# Morris-Lecar model Methods Chapter with Kca and h for K-DR
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)
ds/dt=eps*(sinf(v)-s)/tausca
Ica=(gca*minf(V)+gsca*s)*(V-Vca)
z=Ca/(Ca+Ca0)
v(0)=-40.28
w(0)=.0557
# scaled calcium by 480 
ca(0)=2.099
s(0)=.156
minf(v)=.5*(1+tanh((v-v1)/v2))
winf(v)=.5*(1+tanh((v-v3)/v4))
sinf(v)=.5*(1+tanh((v-v5)/v6))
tauw(v)=1/cosh((v-v3)/(2*v4))
param vk=-84,vl=-60,vca=120
param i=65,gk=8,gl=2,c=20
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 zbar=z
aux icaa=ica
@ method=gear,total=3750,xplot=zbar,yplot=s
@ dt=2.5,dtmin=1e-10,dtmax=5,toler=.01
@ xlo=.5,xhi=.75,ylo=0,yhi=.25,bound=10000
done