#12589: series yields wrong result
--------------------------------------------+------------------------
Reporter: dkrenn | Owner: burcin
Type: defect | Status: new
Priority: major | Milestone: sage-7.0
Component: symbolics | Resolution:
Keywords: series, order, symbolics | Merged in:
Authors: | Reviewers:
Report Upstream: N/A | Work issues:
Branch: | Commit:
Dependencies: | Stopgaps: todo
--------------------------------------------+------------------------
Changes (by rws):
* milestone: sage-6.4 => sage-7.0
Old description:
> The following was posted on
> [https://spreadsheets.google.com/pub?key=pCwvGVwSMxTzT6E2xNdo5fA the
> public bug reports from the notebook interface] by Clemens Heuberger on
> 1/4/2011 and on [http://groups.google.com/group/sage-
> support/browse_thread/thread/2949e7c3d0e84ef/852fabf0229526ed sage-
> support]. In the latter other code examples were posted.
>
> {{{f.series(q,2)}}} for the {{{f}}} defined below yields {{{Order(q^2)}}}
> which is incorrect, as {{{f.subs(q=0)}}} equals {{{1}}} (which is
> correct).
>
> {{{
> sage: var('q')
> sage: f=(q^13362120470/((q - 1)*(q^5 - 1)*(q^21 - 1)*(q^85 - 1)*(q^341 -
> 1)*(q^1365 - 1)*(q^5461 - 1)*(q^21845 - 1)*(q^87381 - 1)*(q^349525 -
> 1)*(q^1398101 - 1)*(q^5592405 - 1)*(q^22369621 - 1)*(q^89478485 -
> 1)*(q^357913941 - 1)*(q^1431655765 - 1)*(q^5726623061 - 1)) +
> 1)/(q^7635497409/((q - 1)*(q^5 - 1)*(q^21 - 1)*(q^85 - 1)*(q^341 -
> 1)*(q^1365 - 1)*(q^5461 - 1)*(q^21845 - 1)*(q^87381 - 1)*(q^349525 -
> 1)*(q^1398101 - 1)*(q^5592405 - 1)*(q^22369621 - 1)*(q^89478485 -
> 1)*(q^357913941 - 1)*(q^1431655765 - 1)*(q^5726623061 - 1)) + 1)
> sage: f.series(q,2)
> Order(q^2)
> sage: f.subs(q=0)
> 1
> }}}
New description:
The following was posted on
[https://spreadsheets.google.com/pub?key=pCwvGVwSMxTzT6E2xNdo5fA the
public bug reports from the notebook interface] by Clemens Heuberger on
1/4/2011 and on [http://groups.google.com/group/sage-
support/browse_thread/thread/2949e7c3d0e84ef/852fabf0229526ed sage-
support]. In the latter other code examples were posted.
{{{f.series(q,2)}}} for the {{{f}}} defined below yields {{{Order(q^2)}}}
which is incorrect, as {{{f.subs(q=0)}}} equals {{{1}}} (which is
correct).
{{{
sage: var('q')
sage: f=(q^13362120470/((q - 1)*(q^5 - 1)*(q^21 - 1)*(q^85 - 1)*(q^341 -
1)*(q^1365 - 1)*(q^5461 - 1)*(q^21845 - 1)*(q^87381 - 1)*(q^349525 -
1)*(q^1398101 - 1)*(q^5592405 - 1)*(q^22369621 - 1)*(q^89478485 -
1)*(q^357913941 - 1)*(q^1431655765 - 1)*(q^5726623061 - 1)) +
1)/(q^7635497409/((q - 1)*(q^5 - 1)*(q^21 - 1)*(q^85 - 1)*(q^341 -
1)*(q^1365 - 1)*(q^5461 - 1)*(q^21845 - 1)*(q^87381 - 1)*(q^349525 -
1)*(q^1398101 - 1)*(q^5592405 - 1)*(q^22369621 - 1)*(q^89478485 -
1)*(q^357913941 - 1)*(q^1431655765 - 1)*(q^5726623061 - 1)) + 1)
sage: f.series(q,2)
Order(q^2)
sage: f.subs(q=0)
1
}}}
The pynac ticket implementing an overflow check is
https://github.com/pynac/pynac/issues/121. When this is done this ticket
should doctest it.
--
--
Ticket URL: <http://trac.sagemath.org/ticket/12589#comment:15>
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 unsubscribe from this group and stop receiving emails from it, send an email
to [email protected].
To post to this group, send email to [email protected].
Visit this group at https://groups.google.com/group/sage-trac.
For more options, visit https://groups.google.com/d/optout.
