# parameters for Safstrom cell
par ie=0
par tau=20,a=.644,b=6.17,rm=50
dx/dt=(a*b^2*(1-cos(x))+(1+cos(x))*(rm*ie-a*b^2))/(tau*b)
done