[BACK]Return to taylor.gp CVS log [TXT][DIR] Up to [local] / OpenXM_contrib / pari / examples
File: [local] / OpenXM_contrib / pari / examples / Attic / taylor.gp (download)

Revision 1.1.1.1 (vendor branch), Sun Jan 9 17:35:30 2000 UTC (24 years, 5 months ago) by maekawa
Branch: PARI_GP
CVS Tags: maekawa-ipv6, VERSION_2_0_17_BETA, RELEASE_20000124, RELEASE_1_2_3, RELEASE_1_2_2_KNOPPIX_b, RELEASE_1_2_2_KNOPPIX, RELEASE_1_2_2, RELEASE_1_2_1, RELEASE_1_1_3, RELEASE_1_1_2
Changes since 1.1: +0 -0 lines 

Import PARI/GP 2.0.17 beta.

\\ originally contributed by Ilya Zakharevich

\\ sample function
f(x) = sin(x)

\\ plot Taylor polynomials of f of index first + i*step <= ordlim
\\ for x in [xmin,xmax].
plot_taylor(xmin=-5, xmax=5, ordlim=16, first=1, step=1) = 
{
  local(T,Taylor_array,s,t,w,h,dw,dh,cw,ch,gh,h1, extrasize = 0.6);

  default(seriesprecision,ordlim+1);
  ordlim -= first; ordlim \= step; ordlim += first;
  T = f(tt); Taylor_array = vector(ordlim+1);
  forstep(i=ordlim+1, 1, -1,
    T += O(tt^(1 + first + (i-1)*step));
    Taylor_array[i] = truncate(T)
  );

  t = plothsizes();
  w=floor(t[1]*0.9)-2; dw=floor(t[1]*0.05)+1; cw=t[5];
  h=floor(t[2]*0.9)-2; dh=floor(t[2]*0.05)+1; ch=t[6];
  h1=floor(h/1.2);

  plotinit(2, w+2*dw, h+2*dh);
  plotinit(3, w, h1);
  s = plotrecth(3, x=xmin,xmax, f(x), 2+8+16+32);
  gh=s[4]-s[3];

  plotinit(3, w, h);
  plotscale(3,s[1],s[2],s[3]-gh*extrasize/2,s[4]+gh*extrasize/2);
  plotrecth(3, x=xmin,xmax,
    concat(f(x), subst(Taylor_array, tt, x)), 4);
  plotclip(3);
  plotcopy(3, 2, dw, dh);

  plotmove(2,floor(dw+w/2-15*cw),floor(dh/2));
  plotstring(2,"Multiple Taylor Approximations");
  plotdraw([2, 0, 0])
}