#hodgkin



#reactions
v_vH:
 $pool = {1.0}h
  alphah * (1 - h) - betah * h
v_vM:
 $pool = {1.0}m
  alpham * (1 - m) - betam * m
v_vN:
 $pool = {1.0}n
  alphan * (1 - n) - betan * n
v_vV:
 $pool = {1.0}V
  barGna * pow(m, 3) * h * (ENa - V) + barGk * pow(n, 4) * (EK - V) + barGleak * (ELeak - V) + extInput


#parameters
EK = -72.0
ELeak = -49.387
ENa = 55.0
Q10 = 3.0
T = 6.3
ah = 0.07
am = 0.1
an = 0.01
barGk = 36.0
barGleak = 0.3
barGna = 120.0
bh = 1.0
bm = 4.0
bn = 0.125
delay1 = 0.0
dur1 = 60.0
dur2 = 0.0
kalphah = 20.0
kalpham = 10.0
kalphan = 10.0
kbetah = 10.0
kbetam = 18.0
kbetan = 80.0
stim1 = 6.0
stim2 = 0.0
t1 = 5.0
valphah = -60.0
valpham = -36.0
valphan = -50.0
vbetah = -30.0
vbetam = -60.0
vbetan = -60.0
default_compartment = 1.0


#initial values
V = -60.0
h = 0.587915337258
m = 0.0580917142295
n = 0.321283379998


#assignment rules
!F t2 = t1 + dur1
!F t3 = t1 + dur1 + delay1
!F t4 = t3 + dur2
!F alphah = pow(Q10, (T - 6.3) / 10) * ah * exp(-(V - valphah) / kalphah)
!F betan = pow(Q10, (T - 6.3) / 10) * bn * exp(-(V - vbetan) / kbetan)
!F betah = pow(Q10, (T - 6.3) / 10) * bh / (1 + exp(-(V - vbetah) / kbetah))
!F alphan = piecewise(pow(Q10, (T - 6.3) / 10) * 0.1, eq(V, valphan), pow(Q10, (T - 6.3) / 10) * an * (V - valphan) / (1 - exp(-(V - valphan) / kalphan)))
!F betam = pow(Q10, (T - 6.3) / 10) * bm * exp(-(V - vbetam) / kbetam)
!F extInput = piecewise(stim1, and(leq(t1, time), lt(time, t2)), stim2, and(leq(t3, time), lt(time, t4)), 0)
!F alpham = piecewise(pow(Q10, (T - 6.3) / 10) * 0.1, eq(V, valpham), pow(Q10, (T - 6.3) / 10) * am * (V - valpham) / (1 - exp(-(V - valpham) / kalpham)))