Ask Your Question
1

improper integral error?

asked 10 years ago

rafarob gravatar image

Is there an alternative way to calculate this improper integral?

integral(log(x)/(x^2 - 1),x,0,infinity)

the integral is convergent to pi^2/4 but Sage says divergent.

Preview: (hide)

2 Answers

Sort by » oldest newest most voted
2

answered 10 years ago

rws gravatar image

Always try the other available "algorithms":

sage: integral(log(x)/(x^2 - 1),x,0,infinity, algorithm='mathematica_free')
-1/2*log(x + 1)*log(x) - 1/2*polylog(2, -x) - 1/2*polylog(2, -x + 1)

See http://sagemath.org/doc/reference/cal...

Preview: (hide)
link

Comments

Though this algorithm requires the internet. Interestingly, using algorithm='sympy' gives AttributeError: 'NonElementaryIntegral' object has no attribute '_sage_'.

kcrisman gravatar imagekcrisman ( 10 years ago )

Unfortunately, the mathematica_free algorithm will not work from within Sage Math Cloud with current connection settings.

calc314 gravatar imagecalc314 ( 10 years ago )

But the correct answer is a numerical value or exact expression pi^2/4. So, algorithm = 'mathematica_free' doesn't work.

rafarob gravatar imagerafarob ( 10 years ago )
1

answered 10 years ago

calc314 gravatar image

With the algorithm set to mathematica_free, Sage is returning the antiderivative even though you are asking for the improper integral. You can get the improper integral by doing the following:

f=integral(log(x)/(x^2 - 1),x,algorithm='mathematica_free')
limit(f,x=oo)-limit(f,x=0)

This returns pi^2/4.

Preview: (hide)
link

Your Answer

Please start posting anonymously - your entry will be published after you log in or create a new account.

Add Answer

Question Tools

1 follower

Stats

Asked: 10 years ago

Seen: 1,554 times

Last updated: Mar 11 '15