# Morris-Lecar model Methods Chapter with Kca
dv/dt = ( I - ICa-(gk*w+gkca*z)*(V-VK)-gl*(V-Vl))/c
dw/dt = phi*(winf(V)-w)/tauw(V)
dca/dt=eps*(-mu*Ica-ca)
Ica=gca*minf(V)*(V-Vca)
z=Ca/(Ca+Ca0)
v(0)=-41
w(0)=.16
ca(0)=26.8
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=120,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 mu=.3,ca0=18,eps=0.002,gkca=.75
aux zbar=z
aux icaa=ica
@ xp=zbar,yp=v,xlo=.3,xhi=.7,ylo=-70,yhi=45,total=2000,dt=2
@ meth=gear,toler=1e-5
@ dtmax=5,dtmin=1e-10,bound=1000
done