# Wiener
param kp=0.005,km=0.005
tau=1/(kp+km)
finf=kp/(kp+km)
N=1000
gamma(f)=(kp*(1-f)+km*f)/N
wiener w
f(0)=0.5
f'=-(f-finf)/tau+sqrt(gamma(f))*w
@ total=1000,trans=0,DT=0.01,xlo=0,xhi=1000,ylo=-0.1,yhi=1.1
@ maxstore=1000000,bounds=10000
@ xplot=t,yplot=x
@ njmp=1
done