# cannot calculate limit

No answer from Sage-6.4.1 for this limit :

var('n'); u = (1+sqrt(n))^(-n); limit(u(n=n+1)/u,n=infinity)

Sage 4.7 gives 0 (correct). Why ?

edit retag close merge delete

1
( 2015-02-01 09:41:48 +0200 )edit

Sort by ยป oldest newest most voted

I do not know exactly why, probably a regression in maxima. If you really want to know, you should provide a smaller interval of versions than (4.7, 6.4.1) so that we can have a look at the guilty commit. That said, you can still ask sympy to do that:

sage: limit(u(n=n+1)/u,n=infinity, algorithm="sympy")
0

more

Thanks ! Note that converting to exponentials works. My best interval is (4.7, 5.8), sorry

( 2015-01-30 13:52:04 +0200 )edit

This does not evaluate (returns the question) in 5.2. So it's been a while. Thierry, maybe you can isolate the Maxima command and report upstream?

( 2015-01-30 15:00:35 +0200 )edit

I never used maxima commands (nor maxima itself), but if you give me some hints me how to do that, i will be glad to.work on it.

( 2015-02-01 18:10:51 +0200 )edit

The maxima command is: limit((1+sqrt(n+1))^(-n-1)/(1+sqrt(n))^(-n),n,inf)

Returns zero with Maxima 5.24.0 and is not evaluated when running in http://sagecell.sagemath.org/

( 2015-02-01 23:36:56 +0200 )edit

Thanks @robert.marik - I reported this upstream at https://sourceforge.net/p/maxima/bugs...

( 2015-02-03 15:22:36 +0200 )edit