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 [email protected]
To unsubscribe from this group, send an email to 
[email protected]
For more options, visit this group at http://groups.google.com/group/sage-devel
URL: http://www.sagemath.org

Reply via email to