Hi,

I recently found an integral that returns bad result. I already reported it to
Maxima tracker, but before they fix it, is there some place (ticket?) in which
such buggy integrals, limits, etc are listed so when bugs are fixed upstream
they can be doctested?

For example I noticed that

sage: integrate(s^2*exp(-(a+b)*s^2 ), s)(s=1,b=0).full_simplify()
-1/2*gamma(3/2, a)/a^(3/2)
sage: integrate(s^2*exp(-a*s^2), s)(s=1).full_simplify()
1/4*(sqrt(pi)*e^a*erf(sqrt(a)) - 2*sqrt(a))*e^(-a)/a^(3/2)

are off by sqrt(pi)/(4*a^(3/2)) (used s=1 to simplify the difference between
expressions).

Regards,
Andrzej.

-- 
To post to this group, send an email to sage-devel@googlegroups.com
To unsubscribe from this group, send an email to 
sage-devel+unsubscr...@googlegroups.com
For more options, visit this group at http://groups.google.com/group/sage-devel
URL: http://www.sagemath.org

Reply via email to