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

2015-11-24 Thread Dima Pasechnik
On Tuesday, 24 November 2015 11:12:13 UTC, Sergey Kirpichev wrote: > > On Monday, November 23, 2015 at 10:38:39 PM UTC+3, Ondrej Certik wrote: >> >> Note that William's result has one more zero in the answer... Which >> one is correct? >> > > Maxima's, of course. btw, sympy's answer is same:

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

2015-11-24 Thread Sergey Kirpichev
On Monday, November 23, 2015 at 10:38:39 PM UTC+3, Ondrej Certik wrote: > > Note that William's result has one more zero in the answer... Which > one is correct? > Maxima's, of course. btw, sympy's answer is same: In [3]: limit(S(2)/5*((S(3)/4)**m - 1)*(a0 - 100) + S(1)/5*(3*(S(3)/4)**m

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

2015-11-24 Thread Sergey Kirpichev
On Monday, November 23, 2015 at 3:43:02 AM UTC+3, William wrote: > > This definitely looks like a bug. In the meantime, a workaround is to > use sympy: > This is not a maxima bug: Maxima 5.34.1 http://maxima.sourceforge.net using Lisp GNU Common Lisp (GCL) GCL 2.6.12 (a.k.a. GCL)

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

2015-11-24 Thread Dima Pasechnik
On Tuesday, 24 November 2015 11:12:13 UTC, Sergey Kirpichev wrote: > > > On Monday, November 23, 2015 at 3:43:02 AM UTC+3, William wrote: >> >> This definitely looks like a bug. In the meantime, a workaround is to >> use sympy: >> > > > This is not a maxima bug: > > Maxima 5.34.1

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

2015-11-24 Thread Sergey B Kirpichev
On Tue, Nov 24, 2015 at 04:12:52AM -0800, Dima Pasechnik wrote: >If you first set domain to complex (and this is the setting Sage uses to >call maxima), you get the same error as from Sage. I have seen this post. My messages arrive too late due to moderation, so they a little dated.

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: [sage-support] Fwd: [sage-release] Bug in limit?

2015-11-23 Thread Dima Pasechnik
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*(3/4)^m + 2)*a0;x > 2/5*((3/4)^m - 1)*(a0 - 100) +

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

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

2015-11-23 Thread kcrisman
> > > I wonder -- to what extent should we be using maxima by default still > for limits, instead of sympy...? At some point, presumably sympy will > be uniformly better than maxima? > > > I've been wondering about this as well (also integrals) for some time. Unfortunately I haven't had

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

2015-11-23 Thread Dima Pasechnik
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 a bug. In the meantime, a workaround

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:

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

2015-11-23 Thread Nils Bruin
On Monday, November 23, 2015 at 11:20:23 AM UTC-8, Dima Pasechnik wrote: > > the bug is not really in maxima, it's in Sage's interface to maxima: > > (%i7) limit(2/5*((3/4)^m - 1)*(a0 - 100) + 1/5*(3*(3/4)^m + > 2)*a0,m,inf); > (%o7) 40 > > The stack overflow

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,

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

2015-11-22 Thread David Joyner
Forwarded to the correct list -- Forwarded message -- From: G. M.-S. Date: Sun, Nov 22, 2015 at 4:34 PM Subject: [sage-release] Bug in limit? To: sage-rele...@googlegroups.com Hello. This is my first post, please be indulgent. Is the following a bug?