For the record, that integral is not fixed by http://trac.sagemath.org/sage_trac/ticket/10187 (Update ecl and maxima).
Volker On Nov 8, 3:30 pm, Andrzej Giniewicz <ggi...@gmail.com> wrote: > 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
