1 | initial version |

I've reported this at http://trac.sagemath.org/ticket/18059. I don't have an answer (yet) for why this is happening, because ordinarily it indicates there is something like a complex number as an intermediate result. But I will note that `numerical_integral(lambda r: y.subs(theta=1,phi=1),0,5)`

yields this problem while `z = y.subs(theta=pi/2,phi=pi/2); numerical_integral(z,0,5)`

does not, and yet `numerical_integral(lambda r: z,0,5)`

does but `numerical_integral(lambda r: 0.240000000000000*e^(-1/20*r)/(4*pi - pi^2 + 2*r^2 + 8), 0,5)`

is fine!

Copyright Sage, 2010. Some rights reserved under creative commons license. Content on this site is licensed under a Creative Commons Attribution Share Alike 3.0 license.