replaced L,.~ proximation has always imation is called is ing Llu,, real a Assume AMS 1. (Received Communi~dtcd. in In Konrad-Zuse-Zentrum WOLFRAM tions, Taylor Printed @ Reprints Variables, Complex Photocopying of Ordinary the If this 1'34 the Taylor TAYLOR a a of ?:a. of Computer Inverse note in the Gordon the nes ucin, and functions, inverse point available directly from directly available February converging haan, the by the we present a simple and efficient algorithm efficient and simple a present we case. Taylor permitted ntd States United Polynomials off KOEPF Taylor R. xo. ler system Algebra and Breach Science Publishers Science Breach and POLYNOMIALS 3mn, function f 24, T? 1994, Locally Then the Then Differential of is Gilbcrt po!ymmia!s; 1993) Functions, and to explicitly by order hsn4, Vol. fiir f license only license of (x) f lnformationstechnik, Heilbronner 25. T,(f,x,O) near of America ?n-n4 f the O((x : DERlVE polynomial pp. (x) I+ in solutions publisher given 23-33 I xo, := of We OF - when R xo)"). order is - however, that of of then mention EXPLICIT FUNCTIONS e-'lxZ 0 n ordinary n demonstrate S.A. times differentiable times for Implicit Functions, Implicit + n (1) off. ax all is that the differential to calculate Taylor to The an n otherwise if of Often E Taylor the in algorithm N, x function general, efficiency of Solutions of Solutions f which it equations. O Str: polynomial represents a represents in 0 7071 10, for however, polynomials shows an interval an the We an iterative an algorithms give T, 1 ht this that global approx- global Berlin (1) implementations of gives implicit func- implicit I has in contain- calcula- practice. an is to be not ap-

Downloaded By: [Koepf, Wolfram] At: 16:20 23 October 2007 Downloaded By: [Koepf, Wolfram] At: 16:20 23 October 2007 2. for In representation We [I], following functions memory overflow or time or overflow memory function n of ail personal patible computations pute; most and cosx An This 24 The applications, functions often functions applications, The of the x2 f, TAYLOR [2], mention example so - is f Computer i.e. A!get;ia 1 DERIVE iayior calculation on. as the as and at a TAYLOR((COS(x) that calculations DERIVE h point the Another formula [3]) f of that f POLYNOMIALS poiynomiai functional "(0) may be may we'll of '(0) function system that Algebra s available is input computer, the in : : of for = o However, do. = example type the Taylor xo of of Taylor x x lim lim :hat -0 -0 f TAYLOR(fjw,xO,n) expression Taylor interest. systems available. In this In available. systems with @i(xo) - is case searched fl'(x) f etitos Further restrictions 1)/xA2,x, f '(x) is series n on and with (x) is polynomials especialljl respect of OF are function polynomials for = = = explicitly hc n many in which x-0 lim x-0 the sinxlx 0: given only implicitly only given lim for IMPLICIT symbolic by h other the 8) W. x to which - is following cos - DERIVE (2). KOEPF 2x easy too h variable the 3628800 - for (x generates the given cos x x2 at x k. f - iarge, f to FUNCTIONS explicit hand the x DERIVE output after output is sin 0) + use, restrictions cases given +---+--- functions note, x at origin x3 (x2 02 70 24 720 40320 DERIVE having = xo x6 running functions - 0, . For x. one can one may we use we does by = 2) sin Taylor /Expand1 0 an an strength x4 also with not may not apply. x example =n DERIVE pniynnmia! be - find is an .I2 support faii every -- implemented valid If loih (see algorithm 3' 1 enough a by 2 1 the closed we [j], IBM for reasons all of order a get x special for Com- order form com- = the xo. all of in n Downloaded By: [Koepf, Wolfram] At: 16:20 23 October 2007 produce of further and To chain and by impiiciily are where arid where which a tion g function Firre A get the so on, so corresponding thus in for typical guarantees theorem rule to rule the two-dimensional , chain two-dimensional F(xo, for each for y terms we each the = given inductively. ihr higher g g(x) present Taylor yo) example produce for (xo,)ioj of function. (xo,yo) = x is which ALR POLYNOMIALS TAYLOR 0. and explicit polynomial considered ar! with is By iterative F!x.g(x,)) y, with the Differzntiating under certain weak conditions weak certain under this and xi functions equation xi + of procedure to by of y; + procedure g be y; order taking = = successively, we 0 1 = locally. a function a of in and get 1 n. the and a neighborhood a the iilc we Of y, to nt circle unit icii~ii~ig produce limits yo < generate the generate IMPLICIT we of 0 > the 0 the for may the ~yudiiuii function a variable y list FUNCTIONS function differentiate 4 of the Taylor of yo a local given the and x. polynomials first The existence x point now + implicit n xo, derivatives (xo, Fl we of by yo) of func- may such the for an ? LJ c Downloaded By: [Koepf, Wolfram] At: 16:20 23 October 2007 g",. Thus DERIVE the ITERATES(DIF(g_, the ITERATES(DIF(g_, IMPLICIT_TAYLOR_YPRIME(f,x, the calculates implement DERIVE lation The and and with 26 As command equation ..,g(n). fourth order fourth iteratively (xo, an of (like Taylor yo) xml ecnie again consider we example the our If F(x,y) = other Computer other yprime y)*yprirne y)*yprime Taylor given (0,l). polynomials Taylor function = algorithmic By is polynomial 0 + + the polynomial thus polynomial the using DIF(g_,x), y, DIF(g_,x), xO, IMPLICIT-TAYLOR(f,x,y,xO,yO,n), is , first of above yo, Algebra limits where implicit n, g-, procedure into procedure g-, yprirne) of procedure yprime, yprirne, W. order systems) functions. := KOEPF is n n given - - aux then n 1)) 1). we get we of represents does DERIVE. by the On the o directly not the eiaie list derivative ucin y function other The the DERIVE function list hand, we hand, given support of = is g(x) derivatives by produced can the given calcu- easily g', by by Downloaded By: [Koepf, Wolfram] At: 16:20 23 October 2007 DERIVE lNVERSETAYLOR(g,y,x, the near case The nomial 3. the The cases limit list results IMPLICIT_TAYLOR(f,x, !MPL!C!T-?.AVLnR_YPRIME(f.x. IMPLICIT-TAYLOR-AUX(f,x, DERIVE We Taylor of iterated TAYLOR last example is equivalent is example last local of (O,g(O)). does not does one mention example in in input an function oh x both the can inverse implicit polynomial ii in limit desired POLYNOMIALS OF An find exist. that calculations. and TAYLOR function application y,xO, function: a for n) closed IMPLICIT-TAYLOR-AUX y, Taylor := y,xO, yo, of implicitly then IF(LIM(DIF(g,y),y,O) n) g-l y. POLYNOMIALS yo, 1 g-l(y) IMPLICIT-TAYLOR(g "Taylor xO, := form representation form polynomjal It is It again, an algorithm an again, n, yo. of at to n, the the the expansion given the -DIF(f,x)/DIF(f, of INVERSE FUNCTIONS general algorithm general local point explicit example explicit g(x) DERIVE algebraic functions, algebraic where OF = solution of in g(0j. 0, - inverse does not exist" not inverse does MLCT FUNCTIONS IMPLICIT upt after output x,x,y,LIM(g,y,O),O, a neighborhood a we y)) of [4] This is used of the essential is the procedure available with available (2). Taylor of (1) I Expand equation the Here to i.e, as calculate series. last section last n) of if the order the ( a is x F(x,y) two = covered which 0 g'. dimensional is is generates of Here a in a special by taking many poly- the is a g we lazy evaluation. lazy does the does the to InverseTaylor[g-,x-, VECTOR(INVERSE-TAYLOR(x*EXP(x), order implemented the Thus converge simplifying the of We x One approximations We 28 INVERSE-TAVLOR(EXP(X~~ INVERSE_TAYLOR((#eAx INVERSE-TAYLOR(X~~, DERIVE to = We work dynamically work convergence have function first five first order consider get -1 get may only n eto that mention input of for example for that the to job to n. s DERIVE'S capabilities use graphical value in redo the redo g-l(y). Taylor are Taylor desired. This method This desired. If the in the we y-, by of some -e-' led calculated already at example. last nJ interval polynomial the ltig both plotting decide polynomials x, As the whole calculation. With calculation. whole in := y,5) + Computer z #e" Thylor InverseSeries[Series[g,(x, the - an calculation -0.367879. 1) to (-e-l,e-') lutain Figure illustration - order x y x x)/2,x,y,5) calculate series 5) of x, The TI,. Algebra y, g as order k), of in and It follows from It of . one function W. of k, .,z inverting memory, 1,s). can the Taylor n the "Taylor "Taylor "Taylor DERIVE KOEPF the may systems, of of Taylor e expect we inverse g-' g-', calculated our method, our 0, store 1 g(x) expansion expansion expansion oyoil fr inverse for polynomials n)], the output to i.e. shows that however, y] explore polynomial e.g. := Taylor one complex the function after can be calculated be can xeX the in that of of of derivatives may Taylor MATHEMATICA inverse inverse inverse polynomial in graph has the the is work analysis principle it principle does static of polynomials a quality Taylor does does not does does of local order with not and g with that not not together of of minimum streams and streams polynomials exceed functions limits n by exist". exist". exist", order is [6], the respect the + DERIVE possible of 1 where radius up given later, with g-l. n e-I. to of to at is

Downloaded By: [Koepf, Wolfram] At: 16:20 23 October 2007 Downloaded By: [Koepf, Wolfram] At: 16:20 23 October 2007 rather than the than rather analytic Then-if before. complex nth differential equation differential and We 4. EQUATIONS consider Taylor TAYLOR initial Obviously in function FIGURE the data polynomial. a neighborhood the solution POLYNOMIAL usual soiution TAYLOR 1. g The ih nltc right analytic with z). g(x) A graph To of POLITOMIALS typical of of find an of SOLU - the g(x) (3)-(4) initial value initial -- its situation initial - = y' --. Taylor y(xo) I xeX IONS = is hand -.-. and F(x, point OF n = coefficients, we is the times priiblcn; Yo. ~r -r y) IMPLICIT an side first xo r. r~Ra nta value initial (we differentiable-we F. 5 -.r Taylor give:: In I hold F!JNCT!ONS vnvcn A"-,-" this polynomials use 5.j on case the case problem an the to UIFFERENTIAL *I explicit first explicit use same of search g-I. the solution (3)-(4) method symbol o its for order for g (3) 29 as is x a Downloaded By: [Koepf, Wolfram] At: 16:20 23 October 2007 order Taylor bian. differentiai For uses iiviPiiCii~i~fiOFi~A'di;(i,x, DSOLVEl The an DSOLVEl_TAYLOR(y/x,x, DERIVE and We appiication the DERIVE with this ytm, with systems, eak further remark polynomial approximations polynomial input algorithmic construction algorithmic TAYLOR(f,x, approach. the equations function derivatives of the y,xO,yO, We the y,xO, ih oyoii coefficients poiynomiai with that aigoriiii~ii y,xO, give only list yo, nj yii, h sm method same the := ihe 5) n, ifrne that difference ~TERATES(DIF(S_, fdlowirig of fur we the DERIVE - 'OX, iiiiplicit xo W. again Taylor KOEPF examples. output refer fiiii~iiiiiis here series solution can yjrf after to t e ple t sle explicit solve to applied be dFk/ay (of i3ii-iy_,nj,d-,i, [I]. / jiieid~ arbii~ai-y hn denotes then of the ww.r.t. 1-1 homogeneous - order) Tqlnr 7)) x rather polynoxia!. the linear Jaco- than first Downloaded By: [Koepf, Wolfram] At: 16:20 23 October 2007 ITERATES(DIF(~-,u)*f DSOLVE2-TAYLOR(f, To DERIVE the and I., by defining and we and Let the Let the In 5. DIFFERENTIAL and get the the limit nth iteratively TAYLOR g(k)(~~) chain case v ., given right given order can function for equatim that ...-. of u apply rule to rule -t Taylor POLYNOMIAL SOLUTIONS is an (k yl, x, :Iftf.r: TAYLOR EQUATIONS expiicit of 2 + y, the same the hand side hand DIF(g_,y)*u get u, y 2) g polynomial xO, -t p-n..mh we -,Iu,,i. YO, yo, second POLYNOMIALS have yl and method , X'!J.X) of n) + ,,iLL.Li,..,aeiL. AiCferenti~hl~ DIF(g-,x),g_,f, : to (5) = of x order DSOLVE2_TAYLORPUX(f, -t evaluate the = F(x,y, as before after before as xo, F(X,~(X)?~'(X)) solution differentiai which yields which OF u) n g(k)(x) OF - IMPLICiT he Again; 2)) of HIGHER ORDER HIGHER a function (5)-(6) equation at calculation we the the FUNCTIONS x, iterative!:, result. Y, point has u, of xO, the three the the YO, x This of = Y form g'k)(xo) diff~r~ntjte 1 xo, n. is done i.e. varia'hies we (k by > take the the 2). x, Downloaded By: [Koepf, Wolfram] At: 16:20 23 October 2007 with differential We gives the with This differeiitial DERIVE We DSOLVEZ-TAYLG~~EXP~U)~~~~ 3 DERIVE 2 remark DERIVE present initial initial method-as the input input Taylor data equations conditions that ecjiiatioii. the function there polynomial following in the of For is yJJ(xj an higher first the examples y(x0) - obvious solution SIN order + lincar order pjx)yl jx),x,y, = Yo case-works generalization differefitia! of W. (m i.i,O,O, jx) order and KOEPF E + 1,5) N) qjx));(xj yl(xo) n. equdon The independent - DERIVE ex4 DERIVE 12 of = following +---+x. thc = Yl 120 x5 r(xj upt after output upt after output given x3 6 of are technique the some v] /Expand] linearity examples to explicit of the Downloaded By: [Koepf, Wolfram] At: 16:20 23 October 2007 [6] [4] [5] [3] 12) [I] References SIGSAM t Wolfram, St. Comp., A. B-92-21 sruhe, Waialae Avenue, W. Heilelberg, W. W. matical W. Rich, Koepf, Koepf, Koepf, Algorithmic development Algorithmic Koepf, Kucpf, Germany, August Rcckood computing, Fachbereich Mathematik Fachbereich J. Bulletin Power A Examples Rich 1993, new Mathematica. and Suite TAYLOR 27 series 195-213. Cityj loih for algorithm ed. (1993). D. [or 304, CA, by in 1992, Stoutemyer, the optr algebra, computer J. Honolulu, A 1991. 20-32. POLYNOMIALS Calmet and Calmet zlgorithmic Proceedings, Notes Lecture system the der development of Freien DERIVE HI, u dig doing fur oe series. power J. calculation 96816-3236. .Cmbl,Itrainl Conference International Campbell, A. Universitat Journal User Manual, Version Manual, User OF of MLCT FUNCTIONS IMPLICIT of nf algebraic functions algebraic In: Symbolic Berlin, formal Artificiai in Computer Science by 19%. usu.Luet n power and Laurent Puiseux. Computation Computer. inieiiigeiicc 2, Soft 2, in Puiseux Addison-Wesley Warehouse, Inc., Warehouse, 13 737, and (1992), Springer, AISMC-I, sy~~bc!ic series, 581-603. Preprint Berlin- math- series, Publ. Karl- 3660