# Morris-Lecar model Methods Chapter with Kca 
dv/dt = ( stim(t)+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)=-59.8
w(0)=.0
# scaled calcium by 480 
ca(0)=0.02
minf(v)=.5*(1+tanh((v-v1)/v2))
winf(v)=.5*(1+tanh((v-v3)/v4))
tauw(v)=1/cosh((v-v3)/(2*v4))
stim(t)=I0*.25*(1+tanh((t-ton)*shp))*(1+tanh((toff-t)*shp))
param vk=-84,vl=-60,vca=120
param i=0,gk=8,gl=2,c=20
param ton=1000,toff=5000,shp=20,i0=70
param v1=-1.2,v2=18
param v6=24,v5=12,tausca=.05
param v3=12,v4=17.4,phi=.0666666666,gca=4
param mu=.025,ca0=1,eps=0.0005,gkca=1,gsca=1
aux iap=stim(t)
aux icaa=ica
@ method=gear,total=8000
@ dt=2.5,dtmin=1e-10,dtmax=5,toler=.01
@ xlo=0,xhi=8000,ylo=-70,yhi=40,bound=10000

done