Menu

#309 pade interfered with by taylor

closed
5
2008-08-06
2003-05-08
No

(C1) taylor(sin(x),x,0,3);
(D1) x-x^3/6
(C2) pade(d1,2,2);
(D2) [6*x/(x^2+6)]

Fine so far.

(C3) taylor(x,x,0,3);
(D3) +x
(C4) pade(d1,2,2);
(D4) [6/7]

What's this?! Pade is being affected by the Taylor
calculation in C3?! Presumably there is some global
flag being set but not reset....

(C5) taylor(exp(x),x,0,3);
(D5) 1+x+x^2/2+x^3/6
(C6) pade(d1,2,2);
(D6) [6*x/(x^2+6)]

Somehow that restored the thing that was causing the
problem.

Discussion

  • Robert Dodier

    Robert Dodier - 2006-04-09
    • labels: --> Lisp Core - Taylor
     
  • Dan Gildea

    Dan Gildea - 2008-08-06

    Logged In: YES
    user_id=1797506
    Originator: NO

    Fixed in pade.lisp rev 1.8: declare varlist and genvar as special

    (%i4) taylor(sin(x),x,0,3);
    (%o4) x-x^3/6

    (%i5) pade(%,2,2);
    (%o5) [6*x/(x^2+6)]

    (%i6) taylor(x,x,0,3);
    (%o6) +x

    (%i7) pade(%o4,2,2);
    (%o7) [6*x/(x^2+6)]

     
  • Dan Gildea

    Dan Gildea - 2008-08-06
    • status: open --> closed
    • assigned_to: nobody --> dgildea
     

Log in to post a comment.