#S Coombes Jan 2007
#Morris-Lecar model 
#See the online tutorial of Bard Ermentrout and John Rinzel at
#http://www.math.pitt.edu/~bard/bardware/meth3/meth3.html

#Type I
param gca=1.33,v3=0.1,v4=0.145,phi=0.3333,J=0.0695

#Type II
#param gca=1.1,v3=0,v4=0.3,phi=0.2,J=0.25

param vk=-0.7,vl=-0.5,vca=1
param gk=2,gl=0.5
param v1=-0.01,v2=0.15

v(0)=0
w(0)=0

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


dv/dt=J-gl*(v-vl)-gk*w*(v-vk)-gca*minf(v)*(v-vca)
dw/dt=lamw(v)*(winf(v)-w)

aux Ica=gca*minf(V)*(V-Vca)
aux Ik=gk*w*(V-Vk)

done