# model for the hh K channel
# 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 scale=40
par ton=5,toff=25,ek=-80
v=vhold+heav(t-ton)*heav(toff-t)*(vnew-vhold)
an=.01*(v+55)/(1-exp(-(v+55)/10))
bn=.125*exp(-(v+65)/80)
markov z 5
{0} {4*an} {0} {0} {0}
{bn} {0} {3*an} {0} {0}
{0} {2*bn} {0} {2*an} {0}
{0} {0} {3*bn} {0} {an}
{0} {0} {0} {4*bn} {0}
aux ik=(z==4)*(v-ek)/scale
aux pot=v
@ meth=euler,dt=.1,total=50
@ yp=ik,yhi=3,ylo=-1,xlo=0,xhi=50,bound=100
done