begin name
clancy2
end name

begin reactions
reaction0000001 {1.0}IC3 = {1.0}IC2
reaction0000002 {1.0}IC2 = {1.0}IF
reaction0000003 {1.0}IF = {1.0}IM1
reaction0000004 {1.0}IM1 = {1.0}IM2
reaction0000005 {1.0}C3 = {1.0}IC3
reaction0000006 {1.0}C2 = {1.0}C3
reaction0000007 {1.0}C2 = {1.0}IC2
reaction0000008 {1.0}C1 = {1.0}C2
reaction0000009 {1.0}C1 = {1.0}IF
reaction0000010 {1.0}IF = {1.0}O
reaction0000011 {1.0}O = {1.0}C1
end reactions

begin rate equations
reaction0000001 = cell*(-(b11*IC2[t]) + a11*IC3[t])
reaction0000002 = cell*(a12*IC2[t] - b12*IF[t])
reaction0000003 = cell*(a4*IF[t] - b4*IM1[t])
reaction0000004 = cell*(a5*IM1[t] - b5*IM2[t])
reaction0000005 = cell*(b3*C3[t] - a3*IC3[t])
reaction0000006 = cell*(b11*C2[t] - a11*C3[t])
reaction0000007 = cell*(b3*C2[t] - a3*IC2[t])
reaction0000008 = cell*(b12*C1[t] - a12*C2[t])
reaction0000009 = cell*(b3*C1[t] - a3*IF[t])
reaction0000010 = cell*(b2*IF[t] - a2*O[t])
reaction0000011 = cell*(-(a13*C1[t]) + b13*O[t])
end rate equations

begin parameters
Fara = 96485.0
Gna = 23.5
Rk = 8314.0
Temp = 310.0
nai = 15.0
nao = 140.0
scale = 1000.0
vhold = -80.0
vtest = -30.0
cell = 1.0
end parameters

begin initial conditions
C1[0] = C1i
C2[0] = C2i
C3[0] = C3i
IC2[0] = IC2i
IC3[0] = IC3i
IF[0] = IFi
IM1[0] = IM1i
IM2[0] = IM2i
O[0] = Oi
end initial conditions

begin initial values
C1i = 0.0
C2i = 0.0
C3i = 1.0
IC2i = 0.0
IC3i = 0.0
IFi = 0.0
IM1i = 0.0
IM2i = 0.0
Oi = 0.0
end initial values

begin assignment rules
a11 := 3.802/(0.1027/E^(v/17) + 0.2/E^(v/150))
a12 := 3.802/(0.1027/E^(v/15) + 0.23/E^(v/150))
a2 := 9.178*E^(0.03369272237196765*v)
b13 := 0.22/E^(0.04926108374384236*(-10 + v))
b3 := 0.0084 + v/50000
a4 := 0.09178000000000001*E^(0.03369272237196765*v)
b4 := 3.7933*^-7/E^(0.12987012987012986*v)
b5 := 7.5866*^-9/E^(0.12987012987012986*v)
a5 := 0.00009661052631578948*E^(0.03369272237196765*v)
Ena := (Rk*Temp*Log[nao/nai])/Fara
b11 := 0.1917/E^(0.04926108374384236*v)
b2 := (0.000060166489970363645*E^(0.04926108374384236*(-10 + v) - 0.0961774074981622*v))/((0.1027/E^(v/12) + 0.25/E^(v/150))*(0.0084 + v/50000))
b12 := 0.2/E^(0.04926108374384236*(-5 + v))
Ina := (Gna*(v - (Rk*Temp*Log[nao/nai])/Fara)*O[t])/(scale*(C1[t] + C2[t] + C3[t] + IC2[t] + IC3[t] + IF[t] + IM1[t] + IM2[t] + O[t]))
a3 := 3.7933*^-7/E^(0.12987012987012986*v)
a13 := 3.802/(0.1027/E^(v/12) + 0.25/E^(v/150))
v := Piecewise[{{vtest, 5 < t && t <= 300}}, vhold]
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