% A single M-L oscillator

p phi=.04,gca=4.4,v3=2,v4=30,eca=120,ek=-84,el=-60
p gk=8,gl=2,v1=-1.2,v2=18,c=20,iapp=20

% current functions
minf(v)=.5*(1.+tanh((v-v1)/v2))
ninf(v)=.5*(1.+tanh((v-v3)/v4))
taun(v)=1/cosh((v-v3)/(2*v4))


v'=(-gl*(v-el)-gk*n*(v-ek)-gca*minf(v)*(v-eca)+iapp)/c
n'=phi*(ninf(v)-n)/taun(v)

i v=-50,n=0.

@ dt=.05,total=250,meth=qualrk,xp=t,yp=v,xlo=0,xhi=250,ylo=-80,yhi=60.,bound=500001,maxstor=5000001



done