v'=(I-gl*(V-el)-gk*n*(v-ek)-gna*minf(v)*(v-ena))/C
n'=(ninf(v)-n)/tau(v)
minf(v)=1/(1+exp(-(v-va)/vb))
ninf(v)=1/(1+exp(-(v-vc)/vd))
tau(v)=tau0
par I=40,vc=-45,tau0=1,gl=8,gk=10,gna=20,el=-78,ek=-90,ena=60
par va=-20,vb=15,vd=5
par c=1
@ xp=t,yp=v,xlo=0,xhi=4,ylo=-80,yhi=0
@ bound=1000,dt=.01,meth=qualrk,total=3.66
init v=-8.1692323,n=.510050
done