# sintayl.ode
# first 5 terms of the Taylor series for sin
z1=x
z2=z1-x^3/6
z3=z2+x^5/120
z4=z3-x^7/5040
z5=z4+x^9/362880
x'=2*pi
init x=-3.1415926
aux y=sin(x)
aux y[1..5]=z[j]
@ total=1,dt=.005
@ xlo=-3.15,xhi=3.15,ylo=-1,yhi=1
@ nplot=6,xp=x,yp=y
@ yp[2..6]=y[j-1],xp[j]=x
done