Status: New
Owner: ----
Labels: Type-Defect Priority-Medium Integration WrongResult
New issue 3761 by [email protected]: incorrect integral evaluation
http://code.google.com/p/sympy/issues/detail?id=3761
sympy gives:
In [1]: integrate(x**(x**x), (x,0,1))
Out[1]: 1/2
a quick numerical estimation gives
import numpy
dx = numpy.linspace(0, 1, num=10000)
sum((dx ** dx ** dx) /10000)
0.57311425488629053
WolframAlpha agrees:
N[Integral[x**(x**x), {x,0,1}], 25]
0.5731215670436195827588303
I have verified this issue in both the 0.7.2 release and the master branch
in git, both on 32- and 64-bit systems with CPython 2.7. I also tried pypy
on the 64-bit system and found the same issue.
--
You received this message because this project is configured to send all
issue notifications to this address.
You may adjust your notification preferences at:
https://code.google.com/hosting/settings
--
You received this message because you are subscribed to the Google Groups
"sympy-issues" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to [email protected].
To post to this group, send email to [email protected].
Visit this group at http://groups.google.com/group/sympy-issues?hl=en.
For more options, visit https://groups.google.com/groups/opt_out.
