# Morris Lecar Model as in book and my notes
#  
#
p Iapp=100
p gc=0
p C=20
p Vk=-84
p gk=8
p Vca=120
p gca=4.4
p Vlk=-60
p glk=2
p phi=.04
p v1=-1.2
p v2=18
p v3=2
p v4=30

aux diff=v-vv

i v=-50
i w=.5

i vv=20
i ww=.2

minf(v)=.5*(1+tanh((v-v1)/v2))
winf(v)=.5*(1+tanh((v-v3)/v4))
tau(v)=1/cosh((v-v3)/(2*v4))

v'=(-gca*minf(v)*(v-Vca)-gk*w*(v-Vk)-glk*(v-Vlk)+Iapp+gc*(vv-v))/C
w'=phi*(winf(v)-w)/tau(v)

vv'=(-gca*minf(vv)*(vv-Vca)-gk*ww*(vv-Vk)-glk*(vv-Vlk)+Iapp+gc*(v-vv))/C
ww'=phi*(winf(vv)-ww)/tau(vv)

@ total=600,dt=.2,xhi=600,ylo=-70,yhi=40

done