x=0 y=0 s=0 c=0 e=1 A=1 d=A*A/2+2 n=1/d len=100 xx=vector(len) yy=vector(len) time=vector(len) *loop t=0 while t < len time[t]=10u*t xx[t]=x yy[t]=y t=t+1 x1=n*(2*e-2*A*c-A*y-A*A/2*x-A*A*s) y1=n*(A*e-A*A*c-A*A/2*y+A*x+2*A*s) s=y1/A c=(e-x1)/A x=x1 y=y1 t end display plot xx time plot yy