# Morris-Lecar model Methods Chapter
dv/dt = ( I - gca*minf(V)*(V-Vca)-gk*w*(V-VK)-gl*(V-Vl)+s(t))/c
dw/dt = phi*(winf(V)-w)/tauw(V)
v(0)=-35
w(0)=0.004
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 i=37.5
param vk=-84,vl=-60,vca=120
param gk=8,gl=2,c=20
param v1=-1.2,v2=18
# Uncomment the ones you like!!
#param_fig1-3 v3=2,v4=30,phi=.04,gca=4.4
#param_fig4-6 v3=12,v4=17.4,phi=.06666667,gca=4
param_fig7-8 v3=12,v4=17.4,phi=.23,gca=4

param s1=0,s2=0,s3=0,t1=117,t2=122,t3=300,t4=305,t5=428,t6=433
# triple pulse stimulus 
s(t)=s1*heav(t-t1)*heav(t2-t)+s2*heav(t-t3)*heav(t4-t)+s3*heav(t-t5)*heav(t6-t)
aux stim=s(t)
@ total=200,dt=.25,xlo=-50,xhi=50,ylo=-.125,yhi=.5,xp=v,yp=w
done