Updates:
Labels: Integration
Comment #2 on issue 1696 by asmeurer: integral of heaviside
http://code.google.com/p/sympy/issues/detail?id=1696
The problem is the same as in issue 1113. diff(Heaviside(x), x) gives
DiracDelta(x), and integrate() doesn't know how to handle that (rewrite it
in terms of Heaviside). If you change diff(Heaviside(x), x) to return 0,
then it works all of a sudden:
In [1]: diff(Heaviside(x), x)
Out[1]: 0
In [2]: integrate(Heaviside(x), x)
Out[2]: xâ‹…Heaviside(x)
I think the better solution would be to implement what is described in
section 5 of the paper linked to in issue 2128. Or, since Heaviside is
essentially a constant, we could have a simple routine in integrate() that
replaces it with a dummy symbol not depending on x, and then to replace it
back when it is done (I think this would work anyway).
--
You received this message because you are subscribed to the Google Groups
"sympy-issues" group.
To post to this group, send email to sympy-iss...@googlegroups.com.
To unsubscribe from this group, send email to
sympy-issues+unsubscr...@googlegroups.com.
For more options, visit this group at
http://groups.google.com/group/sympy-issues?hl=en.