#13526: simple incorrect limit
-----------------------------------+----------------------------------------
Reporter: pang | Owner: burcin
Type: defect | Status: new
Priority: minor | Milestone:
Component: symbolics | Resolution:
Keywords: symbolic, limit | Work issues:
Report Upstream: N/A | Reviewers:
Authors: | Merged in:
Dependencies: | Stopgaps:
-----------------------------------+----------------------------------------
Comment (by kcrisman):
Well, fair enough.
{{{
Maxima 5.28.0 http://maxima.sourceforge.net
using Lisp SBCL 1.0.55.0-abb03f9
(%i2) limit( (3^n+(-2.0)^n)/(3^(n+1)+(-2)^(n+1)),n,inf);
rat: replaced -2.0 by -2/1 = -2.0
rat: replaced -2.0 by -2/1 = -2.0
rat: replaced 3.141592653589793 by 80143857/25510582 = 3.141592653589793
rat: replaced 3.141592653589793 by 80143857/25510582 = 3.141592653589793
rat: replaced -.693147180559945 by -13614799/19642003 = -.693147180559946
rat: replaced -.693147180559945 by -13614799/19642003 = -.693147180559946
rat: replaced -.693147180559945 by -13614799/19642003 = -.693147180559946
rat: replaced -.693147180559945 by -13614799/19642003 = -.693147180559946
rat: replaced -.693147180559945 by -13614799/19642003 = -.693147180559946
rat: replaced .6931471805599453 by 13614799/19642003 = 0.693147180559946
(%o2) 0
(%i3) limit( (3^n+(-2)^n)/(3^(n+1)+(-2)^(n+1)),n,inf);
(%o3) 0
}}}
while
{{{
sage: [l.subs(n=k) for k in [1.,10.^1,10.^2,10.^4,10.^6,10.^8,10.^10]]
[0.0769230769230769, 0.343080200343806, 0.333333333333333,
0.333333333333333, 0.333333333333333, 0.333333333333333,
0.333333333333333]
}}}
--
Ticket URL: <http://trac.sagemath.org/sage_trac/ticket/13526#comment:3>
Sage <http://www.sagemath.org>
Sage: Creating a Viable Open Source Alternative to Magma, Maple, Mathematica,
and MATLAB
--
You received this message because you are subscribed to the Google Groups
"sage-trac" group.
To post to this group, send email to [email protected].
To unsubscribe from this group, send email to
[email protected].
For more options, visit this group at
http://groups.google.com/group/sage-trac?hl=en.
