# a multiplicative noise ratchet with correlations
x'=-sin(x/2)*(j1+sqrt(d1)*y1*l1)-sin(x+x0)*(j2+sqrt(d2)*y2*l2) +  \
    sqrt(d)*(l1*y1+l2*y2+sqrt(1-l1*l1-l2*l2)*w)
par d=3,j1=.7,j2=1,d1=.3,d2=.3,l1=.3,l2=.3,x0=1.57
wiener  y1,y2,w
@ total=500,meth=euler,nout=10,bound=10000
done