# The Morris-Lecar model
#  A simple membrane oscillator.  
#
params iapp=0.0,phi=.333
param v1=-.01,v2=0.15,v3=0.1,v4=0.145,gca=1.
params vk=-.7,vl=-.5,gk=2.0,gl=.5,cm=1,vca=1
minf(v)=.5*(1+tanh((v-v1)/v2))
ninf(v)=.5*(1+tanh((v-v3)/v4))
lamn(v)= phi*cosh((v-v3)/(2*v4))
ica=gca*minf(v)*(v-vca)
v'=  (iapp+gl*(vl-v)+gk*w*(vk-v)-ica)/cm
w'= (lamn(v)*(ninf(v)-w))
# aux I_ca=ica
b v-v'
b w-w'
@ TOTAL=30,DT=.05,xlo=-.6,xhi=1.2,ylo=-.25,yhi=1.2
@ xplot=v,yplot=w
@ dsmin=1e-5,dsmax=0.05,parmin=-.5,parmax=.5,autoxmin=-.1,autoxmax=.5
@ colormap=5
@ autoymax=.5,autoymin=-.5
set vvst {xplot=t,yplot=v,xlo=0,xhi=100,ylo=-.6,yhi=.5,total=100 \
	dt=.5,meth=qualrk,iapp=.1}
Done