# two state process for the Morris-Lecar K channel
v=vhold + (vclamp-vhold)*heav(t-ton)*heav(toff-t)
par vhold=-100,vclamp=0,ton=20,toff=80
par v3=12,v4=17.4,phi=.06666,ek=-84
ninf=.5*(1+tanh((v-v3)/v4))
tau=1/(phi*cosh((v-v3)/(2*v4)))
a=ninf/tau
b=(1-ninf)/tau
markov z 2
{0} {a}
{b} {0}
aux ik=z*(v-ek)
@ total=100,dt=.1,meth=euler
@ yp=ik,xhi=100,yhi=25,ylo=-5
done