begin name
hodgkin
end name

begin reactions
vvH {1.0}$EXTERNAL = {1.0}h
vvM {1.0}$EXTERNAL = {1.0}m
vvN {1.0}$EXTERNAL = {1.0}n
vvV {1.0}$EXTERNAL = {1.0}V
end reactions

begin rate equations
vvH = alphah*(1 - h[t]) - betah*h[t]
vvM = alpham*(1 - m[t]) - betam*m[t]
vvN = alphan*(1 - n[t]) - betan*n[t]
vvV = extInput + barGk*n[t]^4*(EK - V[t]) + barGleak*(ELeak - V[t]) + barGna*h[t]*m[t]^3*(ENa - V[t])
end rate equations

begin 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
defaultcompartment = 1.0
end parameters

begin initial conditions
V[0] = Vi
h[0] = hi
m[0] = mi
n[0] = ni
end initial conditions

begin initial values
Vi = -60.0
hi = 0.587915337258
mi = 0.0580917142295
ni = 0.321283379998
end initial values

begin assignment rules
t2 := dur1 + t1
t3 := delay1 + dur1 + t1
t4 := dur2 + t3
alphah := (ah*Q10^((-6.3 + T)/10))/E^((-valphah + V[t])/kalphah)
betan := (bn*Q10^((-6.3 + T)/10))/E^((-vbetan + V[t])/kbetan)
betah := (bh*Q10^((-6.3 + T)/10))/(1 + E^(-((-vbetah + V[t])/kbetah)))
alphan := Piecewise[{{0.1*Q10^((-6.3 + T)/10), V[t] == valphan}}, (an*Q10^((-6.3 + T)/10)*(-valphan + V[t]))/(1 - E^(-((-valphan + V[t])/kalphan)))]
betam := (bm*Q10^((-6.3 + T)/10))/E^((-vbetam + V[t])/kbetam)
extInput := Piecewise[{{stim1, t1 <= t && t < t2}, {stim2, t3 <= t && t < t4}}, 0]
alpham := Piecewise[{{0.1*Q10^((-6.3 + T)/10), V[t] == valpham}}, (am*Q10^((-6.3 + T)/10)*(-valpham + V[t]))/(1 - E^(-((-valpham + V[t])/kalpham)))]
end assignment rules

begin function definitions
end function definitions

begin events
end events

begin process annotations
end process annotations

begin species annotations
end species annotations

begin units
end units