# Morris Lecar Model of
# Membrane Potential Oscillations

dV/dt = (-gca*minf*(V-Vca)-gk*W*(V-Vk)-gL*(V-VL)-Ikca)/C
dW/dt = (Winf-W)/tau
dCai/dt = rho*(-alpha*gca*minf*(V-Vca) - vpump*Cai)

tau = phi/cosh((V-v3)/(2*v4))
minf = 0.5*(1+tanh((V-v1)/v2))
winf = 0.5*(1+tanh((V-v3)/v4))
Ikca = gkca*Cai*(V-Vk)/(Kkca+Cai)

par C=5300, phi=28.6
par Vk=-75, gk=2700
par Vca=25, gca=1000
par VL=-75, gL=150
par v1=-20, v2=24
par v3=-16, v4=11.2
par Kkca=5, gkca=2000
par rho=0.001, alpha=0.0000045
par vpump=.1

done