#   Parameters

par gkc=10
par gi=1800, gkv=1700, gl=7
par EK=-75, EI=100, EL=-40, EC=100
par r=0.27
par I=0
par eps=1
 
kc=3.3/18

am(v)=0.1*(25+v)/(1-exp(-(v+25)/10))
bm(v)=4*exp(-(v+50)/18)

ah(v)=0.07*exp(-(v+50)/20)
bh(v)=1/(1+exp(-(v+20)/10))

an(v)=0.01*(v+20)/(1-exp(-(v+20)/10))
bn(v)=0.125*exp(-(v+30)/80)
tn(v)=1/(230*(an(v)+bn(v)))


minf(v)=am(v)/(am(v)+bm(v))
hinf(v)=ah(v)/(ah(v)+bh(v))
ninf(v)=an(v)/(an(v)+bn(v))


dv/dt=-gi*minf(v)^3*hinf(v)*(v-EI)-( gkv*n^4+gkc*(u/(1+u)) )*(v-EK)-gl*(v-EL)+I
dn/dt=eps*(ninf(v)-n)/tn(v)
du/dt=r*(-minf(v)^3*hinf(v)*(v-EC)-kc*u)



init v=-60, n=0, u=0.2

#adjust defaults for xpp internal variables

@ MAXSTOR=40000,YP=n,XP=v, XLO=-70,XHI=0,YLO=0,YHI=1
@ METH=adams, DT=0.005, T=50, BOUND=1000

done