share/contrib/rand/Reverse

   1 /*
   2 \r(d8) THIS PROGRAM REVERSES THE FUNCTION Y EXPRESSED AS
   3
   4
   5 A POWER SERIES IN X.  CALL IT WITH
   6
   7
   8                       REV(Y,X,N)
   9
  10
  11 WHERE N IS THE TRUNCATION ORDER
  12 */
  13
  14 \rrev(y,x,n):=(y:taylor(y,x,0,n),eq:x = sum(b[i]*y^i,i,1,n),
  15     sol:solve(makelist(coeff(eq,x,i),i,1,n),makelist(b[i],i,1,n)),y:'y,
  16     eq:x = sum(b[i]*y^i,i,1,n),eq:ev(eq,sol))$