v'= -gl*(v-vl)-gca*m*(v-vca)-gk*n*(v-vk)+i
m=am(v)/(am(v)+bm(v))
n'=an(v)*(1-n)-bn(v)*n
init v=-70
am(v)=exp((v-v1)/v2)
bm(v)=exp(-(v-v3)/v4)
an(v)=exp((v-v5)/v6)
bn(v)=exp(-(v-v7)/v8)
par vl=-70,vk=-90,vca=150,i=0,gca=1,gk=2.5,gl=.1
par v1=-8,v2=8,v3=-11,v4=8
par v5=16,v6=8,v7=-26,v8=8
@ meth=cvode,toler=1e-5,atol=1e-4
@ xp=V,yp=n,xlo=-80,xhi=30,total=50,ylo=-.125,yhi=1.125
done