# Morris Lecar Model of
# Membrane Potential Oscillations

dV/dt = (-gca*minf*(V-Vca)-gk*W*(V-Vk)-gL*(V-VL)+Iapp)/C
dW/dt = (Winf-W)/tau

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

par C=20, phi=25
par Vk=-84, gk=8
par Vca=120, gca=4.4
par VL=-60, gL=2
par v1=-1.2, v2=18
par v3=2, v4=30
par Iapp=0

done