a ¬= 0`` and ``SingularityFunction(x, a, n + 1)`` if ``n < 0``. - If the node is a multiplication or power node having a SingularityFunction term we rewrite the whole expression in terms of Heaviside and DiracDelta and then integrate the output. Lastly, we rewrite the output of integration back in terms of SingularityFunction. - If none of the above case arises, we return None. Examples ======== >>> from sympy.integrals.singularityfunctions import singularityintegrate >>> from sympy import SingularityFunction, symbols, Function >>> x, a, n, y = symbols('x a n y') >>> f = Function('f') >>> singularityintegrate(SingularityFunction(x, a, 3), x) SingularityFunction(x, a, 4)/4 >>> singularityintegrate(5*SingularityFunction(x, 5, -2), x) 5*SingularityFunction(x, 5, -1) >>> singularityintegrate(6*SingularityFunction(x, 5, -1), x) 6*SingularityFunction(x, 5, 0) >>> singularityintegrate(x*SingularityFunction(x, 0, -1), x) 0 >>> singularityintegrate(SingularityFunction(x, 1, -1) * f(x), x) f(1)*SingularityFunction(x, 1, 0) Nréé)éÿÿÿÿéþÿÿÿ) ZhasrÚfuncrÚargsZ is_positiveÚis_zeroZis_MulZis_PowZrewriterr)ÚfÚxÚaÚnÚexpr©rút/Users/vegardjervell/Documents/master/model/venv/lib/python3.9/site-packages/sympy/integrals/singularityfunctions.pyÚsingularityintegrates+       rN)Zsympy.functionsrrZ sympy.corerZsympy.integralsrrrrrrÚs