params kp=0.5, km=0.5
params c=2, gl=0.5, gch=1.0, vl=-70, vch=20

v(0)=0

markov n 2
{0}{kp}
{km}{0}

v'=(-gl*(v-vl)-gch*n*(v-vch))/c
aux n=n
@ total=100,trans=0 ,DT=.001,xlo=0,xhi=100,ylo=-80,yhi=0 
@ maxstore=1000000,bounds=10000
done