Sage says that a divergent integral converges

When i use the command N(integrate(e^(-x)/x,x,0,oo)) Sage returns 37.2, but this integral obviously diverges. Is there any way to fix this?

1 answer

Sort by » oldest newest most voted

Thanks for reporting!

Actually, it is a known bug, see trac ticket 14274.

Sympy seems not able to solve it either, though it is able to see that the integral of 1/x diverges at 0.

sage: import sympy
sage: x = sympy.Symbol('x')
sage: from sympy import oo
sage: sympy.integrate(1/x,(x,0,oo))
nan

Asked: 2015-11-10 17:07:47 -0500

Seen: 98 times

Last updated: Nov 11 '15