time'=1
markov z 2
{0} {r}
{1000} {0}
r=r0+r1*sin(omega*t)
par r0=.2,r1=.19,omega=.05
@ meth=euler,dt=.05,maxstor=100000,bound=100000,total=2000
done