# Morris-Lecar model Methods Chapter
dv/dt = ( I - gca*minf(V)*(V-Vca)-gk*w*(V-VK)-gl*(V-Vl))/c
dw/dt = phi*(winf(V)-w)/tauw(V)
ds/dt=sp
dsp/dt=-s/(alpha*alpha)-(1/alpha+1/alpha)*sp
v(0)=31.64
w(0)=0.2316
s(0)=.07355
sp(0)=.0006
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 vk=-84,vl=-60,vca=120
param i=50,gk=8,gl=2,c=20
param v1=-1.2,v2=18
param alpha=5,vt=0,vsyn=0
# Uncomment the ones you like!!
# param_fig1-3 v3=2,v4=30,phi=.04,gca=4.4
param4-6 v3=12,v4=17.4,phi=.06666667,gca=4
#param_fig7-8 v3=12,v4=17.4,phi=.23,gca=4
global 1 {v-vt} {sp=sp+1/(alpha*alpha)}
@ total=75.8, dt=.2, xlo=0,xhi=80,ylo=-45,yhi=35
done