# chirikov.ode
p'=mod(p-k*sin(x),2*pi)
x'=mod(x+p-k*sin(x),2*pi)
par k=.9716
@ total=50000,meth=discrete,lt=0,maxstor=50000
@ xp=p,yp=x,xlo=0,ylo=0,xhi=6.3,yhi=6.3
done