# logjul.ode
par r=4.5
s=sign(ran(1)-.5)
y'=.5*(1+s*sqrt(r*r-4*r*y)/r)
init y=.123
@ meth=discrete,total=4000
@ xp=y,yp=y,trans=100,lt=0
done