#11493: Wrong integral computation
------------------------+---------------------------------------------------
Reporter: kcrisman | Owner: burcin
Type: defect | Status: new
Priority: major | Milestone: sage-4.7.1
Component: calculus | Keywords:
Work_issues: | Upstream: N/A
Reviewer: | Author:
Merged: | Dependencies:
------------------------+---------------------------------------------------
Old description:
> Reported by Nicolas Privault via Rado Kirov:
> {{{
> sage: integrate(x^(0.5)*e^(-x-x^2),x,0,infinity).n()
> 0.886226925452758
>
> sage: integrate(x^(1/2)*e^(-x-x^2),x,0,infinity).n()
> 0.32015709036013135
>
> which give two different results, and the second one is correct.
>
> Strangely again, one also has
>
> sage: integrate(x^(0.5)*e^(-x),x,0,infinity).n()
> 0.886226925452758
> }}}
New description:
Reported by Nicolas Privault via Rado Kirov:
{{{
sage: integrate(x^(0.5)*e^(-x-x^2),x,0,infinity).n()
0.886226925452758
sage: integrate(x^(1/2)*e^(-x-x^2),x,0,infinity).n()
0.32015709036013135
which give two different results, and the second one is correct.
}}}
--
Comment(by kcrisman):
The other thing was correct, I believe, so I'm not sure what the poster
thought was wrong with it.
{{{
sage: numerical_integral(x^(0.5)*e^(-x),0,100000)
(0.88622692749267307, 4.3485198683228606e-07)
}}}
--
Ticket URL: <http://trac.sagemath.org/sage_trac/ticket/11493#comment:2>
Sage <http://www.sagemath.org>
Sage: Creating a Viable Open Source Alternative to Magma, Maple, Mathematica,
and MATLAB
--
You received this message because you are subscribed to the Google Groups
"sage-trac" group.
To post to this group, send email to [email protected].
To unsubscribe from this group, send email to
[email protected].
For more options, visit this group at
http://groups.google.com/group/sage-trac?hl=en.
