# determ.ode
#
# The deterministic Morris-Lecar model for electrical activity.
# See p. 35 of the Fall et al. text.

v(0)=-60
w(0)=0.01

# applied current
par iapp=150

# conductances
par gca=4, gk=8, gleak=2

# reversal potentials
par vca=120, vk=-84, vleak=-60

# others
par phi=0.04
num c=20,v1=-1.2,v2=18,v3=2,v4=30

# infinity and time constant functions
minf = 0.5*(1+tanh((v-v1)/v2))
winf = 0.5*(1+tanh((v-v3)/v4))
tauw = 1/(phi*cosh((v-v3)/(2*v4)))

# currents
ica = gca*minf*(v-vca)
ik = gk*w*(v-vk)
ileak = gleak*(v-vleak)

# Differential Equations

v' = -(ica+ik+ileak-iapp)/c
w' = (winf-w)/tauw

# auxilary variables
aux minf=minf

# These are the xpp parameters
@ meth=cvode, toler=1.0e-10, atoler=1.0e-10, dt=1, total=1000, maxstor=20000
@ bounds=10000000, xp=t, yp=v
@ xlo=-1, xhi=200, ylo=-75, yhi=75

done