4
v v=-60.899,w=.014873
a ica,ik
param vk=-84,vl=-60,vca=120
param i=0,gk=8,gl=2, gca=4, c=20
# type II -- v3=2,v4=30,phi=.04 -- like HH 
# type I -- v3=12,v4=17.4,phi=.06666667
param v1=-1.2,v2=18,v3=2,v4=30,phi=.04
user minf 1 .5*(1+tanh((arg1-v1)/v2))
user winf 1 .5*(1+tanh((arg1-v3)/v4))
user lamw 1 phi*cosh((arg1-v3)/(2*v4))
odev  (i+gl*(vl-v)+gk*w*(vk-v)+gca*minf(v)*(vca-v))/c
odew lamw(v)*(winf(v)-w)
# currents
oica gca*minf(v)*(v-vca)
oik gk*w*(v-vk)
done