[sympy] Re: [sage-support] Fwd: [sage-release] Bug in limit?

2015-11-23 Thread William Stein
On Mon, Nov 23, 2015 at 12:17 PM, Dima Pasechnik wrote: > > > On Monday, 23 November 2015 19:38:39 UTC, Ondrej Certik wrote: >> >> On Mon, Nov 23, 2015 at 12:20 PM, Dima Pasechnik wrote: >> > >> > >> > On Monday, 23 November 2015 00:43:02 UTC, William wrote: >> >> >> >> This definitely looks like

[sympy] Re: [sage-support] Fwd: [sage-release] Bug in limit?

2015-11-23 Thread Ondřej Čertík
On Mon, Nov 23, 2015 at 12:20 PM, Dima Pasechnik wrote: > > > On Monday, 23 November 2015 00:43:02 UTC, William wrote: >> >> This definitely looks like a bug. In the meantime, a workaround is to >> use sympy: >> >> sage: var('m a0') >> (m, a0) >> sage: x=2/5*((3/4)^m - 1)*(a0 - 100) + 1/5*(3*

Re: [sympy] Re: [sage-support] Fwd: [sage-release] Bug in limit?

2015-11-23 Thread Aaron Meurer
(note, I am not on the sage list or gms list, so this probably won't make it there unless someone forwards it) SymPy's limit primarily uses the Gruntz algorithm, which is fairly capable. I'm not an expert on it, so others will be able to comment in more detail, but as far as I know, it's mostly re

[sympy] Re: [sage-support] Fwd: [sage-release] Bug in limit?

2015-11-22 Thread William Stein
This definitely looks like a bug. In the meantime, a workaround is to use sympy: sage: var('m a0') (m, a0) sage: x=2/5*((3/4)^m - 1)*(a0 - 100) + 1/5*(3*(3/4)^m + 2)*a0;x 2/5*((3/4)^m - 1)*(a0 - 100) + 1/5*(3*(3/4)^m + 2)*a0 sage: limit(x, m=oo) [BAD] sage: limit(x, m=oo, algorithm='sympy