# 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)=-16
w(0)=0.014915
minf(v)=.5*(1+tanh((v-va)/vb))
winf(v)=.5*(1+tanh((v-vc)/vd))
tauw(v)=1/cosh((v-vc)/(2*vd))
param vk=-84,vl=-60,vca=120
param i=0,gk=8,gl=2,c=20
param va=-1.2,vb=18
# Uncomment the ones you like!!
par1-3 vc=2,vd=30,phi=.04,gca=4.4
set typeII {vc=2,vd=30,phi=.04,gca=4.4}
set typeI {vc=12,vd=17.4,phi=.06666667,gca=4}
set typeIb {vc=12,vd=17.4,phi=.23,gca=4}
#par4-6 vc=12,vd=17.4,phi=.06666667,gca=4
#par7-8 vc=12,vd=17.4,phi=.23,gca=4
param s1=0,s2=0,t1=50,t2=55,t3=500,t4=550
# double pulse stimulus 
s(t)=s1*heav(t-t1)*heav(t2-t)+s2*heav(t-t3)*heav(t4-t)
@ total=150,dt=.25,xlo=-75,xhi=75,ylo=-.25,yhi=.5,xp=v,yp=w
done