This is how impulse function or the singularity function with negative exponents behaves currently
>>> plot(SingularityFunction(x, 1, -1)) # can be written as: Piecewise((oo, Eq(x - 1, 0)), (0, True)) [image: Screenshot - 27-07-2019 , 13_55_01.png] <about:invalid#zClosurez> This seems to be correct for the case (*x-1) != 0, *where it comes out to be zero.. But this plot doesn't stand good for the case *(x-1) = 0. *I believe this is because it is discontinuous at this point. I guess there has been an issue with the sympy.plot() in plotting discontinuities. There is an issue #16572 <https://github.com/sympy/sympy/issues/16572> which mentions other similar examples. Currently, I am looking to plot singularity functions with negative order in order to plot the diagram of the *beam*, where the load is in the form of singularity function. PR #17240 <https://github.com/sympy/sympy/pull/17240> has some prior discussions on it. Now, since we are able to plot for uniformly distributed load, varying load and parabolic load, we just want to plot the same for point loads(order = -1) and moment loads(order=-2). This post is just for discussing how can we use sympy.plot() to plot singularity function with negative order. Is the current plot correct? Or what should it be in general? -- You received this message because you are subscribed to the Google Groups "sympy" group. To unsubscribe from this group and stop receiving emails from it, send an email to sympy+unsubscr...@googlegroups.com. To view this discussion on the web visit https://groups.google.com/d/msgid/sympy/5a952e79-ca31-4055-bf97-fdf446d33477%40googlegroups.com.