# model for the hh Na channel
# due to patlack
# click on Numerics Stochastic Compute
# z=0 to 0   N times  N=10 100 1000 eg
# after done click on  stochastics Mean 
# escape to main menu and redraw to get average
# 
par vhold=-100,vnew=10
par ton=5,toff=25,ena=50
par k1=.24,k2=.4,k3=1.5
v=vhold+heav(t-ton)*heav(toff-t)*(vnew-vhold)
am=.1*(v+40)/(1-exp(-(v+40)/10))
bm=4*exp(-(v+65)/18)
ah=.07*exp(-(v+65)/20)
markov z 5
{0} {3*am} {0} {0} {0}
{bm} {0} {2*am} {0} {k1}
{0} {2*bm} {0} {am} {k2}
{0} {0} {3*bm} {0} {k3}
{0} {0} {0} {0} {ah}
aux cond=(z==3)
aux ina=(z==3)*(v-ena)/40
aux pot=v
@ meth=euler,dt=.1,total=25
@ yp=ina,ylo=-1.5,xlo=0,xhi=15,bound=100
done
bh=1/(1+exp(-(v+35)/10))
init h=.996,m=.0005
m'=am*(1-m)-bm*m
h'=ah*(1-h)-bh*h
aux det=h*m^3