# Morris-Lecar model Methods Chapter
dv/dt = ( I(t) - gca*minf(V)*(V-Vca)-gk*w*(V-VK)-gl*(V-Vl))/c
dw/dt = phi*(winf(V)-w)/tauw(V)
v(0)=-41.84
w(0)=0.002
minf(v)=.5*(1+tanh((v-v1)/v2))
winf(v)=.5*(1+tanh((v-v3)/v4))
tauw(v)=1/cosh((v-v3)/(2*v4))
param i0=30,ip=30,ton=30,toff=50
param vk=-84,vl=-60,vca=120
param gk=8,gl=2,c=10
param v1=-1.2,v2=18
par4-6 v3=12,v4=17.4,phi=.06666667,gca=4
i(t)=i0+ip*heav(t-ton)*heav(toff-t)
@ total=150,dt=.25,xlo=0,xhi=150,ylo=-60,yhi=50,xp=t,yp=v
done