# rotnum.ode
# the rotation number for the std map
f(x)=x+b+a*sin(2*pi*x)
x'=f(x)
b'=b
par a=.15
init b=0
init x=0
aux rho=x/max(t,1)
@ bound=100000
@ total=1000,trans=1000,meth=disc
@ xp=b,yp=rho,xlo=0,ylo=-0.2,xhi=1.001,yhi=1.2
@ rangeover=b,rangestep=200,rangelow=0,rangehigh=1,rangereset=no
done