#file morris_lecar.ode
#Morris-Lecar model of excitable barnacle muscle fiber
#adapted from Morris and Lecar (1981) Biophysical Journal 35 pp. 193-231
#Figures 8.6 and 8.7

par C=20 
par gbarCa=4.4 
par ECa=120
par gbarK=8 
par EK=-84
par gleak=2 
par Eleak=-60
par v1=-1.2
par v2=18 
par v3=2
par v4=30
par phi=0.04 
par tau=0.8

Iapplied=0

#model equations
m_inf=(0.5*(1+tanh((V-v1)/v2)))
w_inf=(0.5*(1+tanh((V-v3)/v4)))

V' = (1/C)*(gbarCa*m_inf*(ECa-V) + gbarK*w*(EK-V) + gleak*(Eleak-V)+Iapplied)
w' = phi*(w_inf-w)/(tau)


# initial condition
init V=-35, w=0.0149

# final time
@ total=100
# set the plotting window size:
@ xlo=0, xhi=100, ylo=-60, yhi=20

# run the simulation in the GUI by selecting (I)nitialconds|(G)o

#generate a phase plane by choosing (V)iewaxis|(2)D and assigning the axis
#the nullclines can be generated with (N)ullclies|(N)ew
#a direction field can be generated with (D)ir.field/flow|(D)irect field

done