# mlII.ode

dv/dt = ( I - gca*minf(V)*(V-Vca)-gk*w*(V-VK)-gl*(V-Vl)+s(t))/c
dw/dt = phi*(winf(V)-w)/tauw(V)
v(0)=-60.855
w(0)=0.014915
minf(v)=.5*(1+tanh((v-v1)/v2))
winf(v)=.5*(1+tanh((v-v3)/v4))
tauw(v)=1/cosh((v-v3)/(2*v4))

param i=0,vk=-84,vl=-60,vca=120
param gk=8,gl=2,c=20
param v1=-1.2,v2=18
param v3=2,v4=30,phi=.04,gca=4.4

# plot parameters: 
@ total=150,dt=.25,ylo=-75,yhi=45,xlo=0,xhi=150,xp=t,yp=va=4.4

done