rtoy / maxima

A Clone of Maxima's repo
Other
0 stars 0 forks source link

pade interfered with by taylor #1564

Closed rtoy closed 3 months ago

rtoy commented 3 months ago

Imported from SourceForge on 2024-07-04 19:42:02 Created by macrakis on 2003-05-08 20:50:59 Original: https://sourceforge.net/p/maxima/bugs/309


(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.

rtoy commented 3 months ago

Imported from SourceForge on 2024-07-04 19:42:04 Created by robert_dodier on 2006-04-09 20:43:42 Original: https://sourceforge.net/p/maxima/bugs/309/#c8d5


rtoy commented 3 months ago

Imported from SourceForge on 2024-07-04 19:42:07 Created by dgildea on 2008-08-06 22:53:21 Original: https://sourceforge.net/p/maxima/bugs/309/#11d7


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)]

rtoy commented 3 months ago

Imported from SourceForge on 2024-07-04 19:42:11 Created by dgildea on 2008-08-06 22:53:21 Original: https://sourceforge.net/p/maxima/bugs/309/#5649