It's supposed to return sqrt(2)/pi. I get the result you post in master,
and the correct result in my new branch "gruntz".
I take it you still get nan in the new branch? :/
On 24.01.2012 19:11, Joachim Durchholz wrote:
Am 24.01.2012 17:53, schrieb Tom Bachmann:
Executing the following two commands triggers the problem for me (i.e.
the second result computed is wrong):
gruntz(gamma(x+1)/sqrt(2*pi)- exp(-x)*(x**(x+S(1)/2) +
x**(x-S(1)/2)/12), x, oo)
gruntz(exp(x)*exp(x**2)*(erf(x+1/exp(x))-erf(x)), x, oo)
My (ness01's) branch gruntz contains a fix. Can you verify if this works
for you?
I'm getting this (Linux 64bit):
>>> gruntz(gamma(x+1)/sqrt(2*pi)- exp(-x)*(x**(x+S(1)/2) + x**(x
S(1)/2)/12), x, oo)
∞
>>> gruntz(exp(x)*exp(x**2)*(erf(x+1/exp(x))-erf(x)), x, oo)
nan
No noticeable delays.
Not sure whether that qualifies as "works" or "does not work".
(I never grokked integral solving, so I have no idea what gruntz etc.
are supposed to return.)
--
You received this message because you are subscribed to the Google Groups
"sympy" group.
To post to this group, send email to sympy@googlegroups.com.
To unsubscribe from this group, send email to
sympy+unsubscr...@googlegroups.com.
For more options, visit this group at
http://groups.google.com/group/sympy?hl=en.
