Learn how easy it is to sync an existing GitHub or Google Code repo to a SourceForge project! See Demo

Close

#1161 coeff, ratcoef on Taylor series around inf

open
nobody
5
2007-06-29
2007-04-02
Anonymous
No

It seems that coeff() and ratcoef() have trouble
extracting coefficients from Taylor expansions
around infinity:

(%i119) taylor(1/x, x, inf, 1);
1
(%o119)/T/ - + . . .
x
(%i120) ratcoef(taylor(1/x, x, inf, 1), x, -1);
(%o120)/R/ 0
(%i121) coeff(taylor(1/x, x, inf, 1), x, -1);
(%o121)/R/ 0

A work-around is to get rid of the Taylor property, e.g., by wrapping the expression in subst(0, 0, ...).

Discussion

  • Logged In: YES
    user_id=588346
    Originator: NO

    Yes, ratcoeff is a bit confused by taylor expansions at infinity:

    qq: taylor(sum((i-3+100)*x^(i-3),i,0,5),x,inf,10)
    => 102*x^2+101*x+100+99/x+98/x^2+97/x^3 + ...

    makelist(ratcoeff(qq,x,i),i,-3,4) =>
    => [0, 102, 101, 100, 99, 98, 97, 0]

    In other words, it treats the 1/x^2 as the ratcoeff(...,x,2) term.

    So the workaround is to negate the power... but of course that will become exactly wrong when the bug is fixed....

     
  • Robert Dodier
    Robert Dodier
    2007-06-29

    • labels: --> Lisp Core - Polynomials