#
v'=-(gl*(v-vl)+gk*w*(v-vk)+gna*x*(v-vna))+i
w'=(winf(v)-w)/tauw(v)
# x'=(xinf(v)-x)/taux(v)
x=xinf(v)
par i=0,gk=10,gl=.1,gna=2.5,vk=-85,vl=-70,vna=120
xinf(v)=.5*(1+tanh((v-vxt)/vxs))
taux(v)=1/(phix*cosh((v-vxt)/(2*vxs)))
winf(v)=.5*(1+tanh((v-vwt)/vws))
tauw(v)=1/(phiw*cosh((v-vwt)/(2*vws)))
par vxt=-52,vxs=15,vwt=-65,vws=20,phix=1,phiw=.1
done