rtoy / maxima

A Clone of Maxima's repo
Other
0 stars 0 forks source link

limit bugs #3100

Closed rtoy closed 1 week ago

rtoy commented 1 week ago

Imported from SourceForge on 2024-07-07 13:43:08 Created by andrejv on 2008-08-31 10:31:04 Original: https://sourceforge.net/p/maxima/bugs/1485


ex and ex1 are equivalent (replace 2^(2*N) with 4^N).

(%i1) ex:(%pi*N^(2*N+1)*2^(2*N))/((2*N-1)^(2*N)*(2*%e*N+%e)); (%o1) (%pi*N^(2*N+1)*2^(2*N))/((2*N-1)^(2*N)*(2*%e*N+%e)) (%i2) ex1:(%pi*N^(2*N+1)*4^N)/((2*N-1)^(2*N)*(2*%e*N+%e)); (%o2) (%pi*N^(2*N+1)*4^N)/((2*N-1)^(2*N)*(2*%e*N+%e)) (%i3) limit(ex, N, inf); (%o3) 0 (%i4) limit(ex*%e, N, inf); (%o4) inf (%i5) limit(ex1, N, inf); (%o5) %pi/2 (%i6) limit(ex1/%e, N, inf); (%o6) 0 (%i7) limit(ex1*%e, N, inf); (%o7) inf

The only correct answer is %o5. tlimit is correct for all examples.

Andrej

rtoy commented 1 week ago

Imported from SourceForge on 2024-07-07 13:43:10 Created by dgildea on 2008-08-31 18:25:27 Original: https://sourceforge.net/p/maxima/bugs/1485/#f8b9


Logged In: YES user_id=1797506 Originator: NO

Fixed in limit.lisp rev 1.55.

(%i11) ex:(%pi*N^(2*N+1)*2^(2*N))/((2*N-1)^(2*N)*(2*%e*N+%e)); (%o11) %pi*N^(2*N+1)*2^(2*N)/((2*N-1)^(2*N)*(2*%e*N+%e)) (%i12) ex1:(%pi*N^(2*N+1)*4^N)/((2*N-1)^(2*N)*(2*%e*N+%e)); (%o12) %pi*N^(2*N+1)*4^N/((2*N-1)^(2*N)*(2*%e*N+%e)) (%i13) limit(ex, N, inf); (%o13) %pi/2 (%i14) limit(ex*%e, N, inf); (%o14) %e*%pi/2 (%i15) limit(ex1, N, inf); (%o15) %pi/2 (%i16) limit(ex1/%e, N, inf); (%o16) %e^-1*%pi/2 (%i17) limit(ex1*%e, N, inf); (%o17) %e*%pi/2

rtoy commented 1 week ago

Imported from SourceForge on 2024-07-07 13:43:13 Created by dgildea on 2008-08-31 18:25:27 Original: https://sourceforge.net/p/maxima/bugs/1485/#2e28