#clancy2



#reactions
reaction_0000001:
 {1.0}IC3 = {1.0}IC2
  cell * (IC3 * a11 - IC2 * b11)
reaction_0000002:
 {1.0}IC2 = {1.0}IF
  cell * (IC2 * a12 - IF * b12)
reaction_0000003:
 {1.0}IF = {1.0}IM1
  cell * (IF * a4 - IM1 * b4)
reaction_0000004:
 {1.0}IM1 = {1.0}IM2
  cell * (IM1 * a5 - IM2 * b5)
reaction_0000005:
 {1.0}C3 = {1.0}IC3
  cell * (C3 * b3 - IC3 * a3)
reaction_0000006:
 {1.0}C2 = {1.0}C3
  cell * (C2 * b11 - C3 * a11)
reaction_0000007:
 {1.0}C2 = {1.0}IC2
  cell * (C2 * b3 - IC2 * a3)
reaction_0000008:
 {1.0}C1 = {1.0}C2
  cell * (C1 * b12 - C2 * a12)
reaction_0000009:
 {1.0}C1 = {1.0}IF
  cell * (C1 * b3 - IF * a3)
reaction_0000010:
 {1.0}IF = {1.0}O
  cell * (IF * b2 - a2 * O)
reaction_0000011:
 {1.0}O = {1.0}C1
  cell * (-C1 * a13 + O * b13)


#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


#initial values
C1 = 0.0
C2 = 0.0
C3 = 1.0
IC2 = 0.0
IC3 = 0.0
IF = 0.0
IM1 = 0.0
IM2 = 0.0
O = 0.0


#assignment rules
!F a11 = 3.802 / (0.1027 * exp(-v / 17.0) + 0.2 * exp(-v / 150.0))
!F a12 = 3.802 / (0.1027 * exp(-v / 15.0) + 0.23 * exp(-v / 150.0))
!F a2 = 9.178 * exp(v / 29.68)
!F b13 = 0.22 * exp(-(v - 10.0) / 20.3)
!F b3 = 0.0084 + 2.000000e-05 * v
!F a4 = (9.178 * exp(v / 29.68)) / 100.0
!F b4 = 3.7933 * 1.000000e-07 * exp(-v / 7.7)
!F b5 = (3.7933 * 1.000000e-07 * exp(-v / 7.7)) / 50.0
!F a5 = (9.178 * exp(v / 29.68)) / (9.5 * 10000.0)
!F Ena = Rk * Temp / Fara * log(nao / nai)
!F b11 = 0.1917 * exp(-v / 20.3)
!F b2 = (3.802 / (0.1027 * exp(-v / 12.0) + 0.25 * exp(-v / 150.0))) * (9.178 * exp(v / 29.68)) * (3.7933 * 1.000000e-07 * exp(-v / 7.7)) / ((0.22 * exp(-(v - 10.0) / 20.3)) * (0.0084 + 2.000000e-05 * v))
!F b12 = 0.2 * exp(-(v - 5.0) / 20.3)
!F Ina = Gna * (O / (IC3 + IC2 + IF + IM1 + IM2 + C3 + C2 + C1 + O)) * (v - (Rk * Temp / Fara * log(nao / nai))) / scale
!F a3 = 3.7933 * 1.000000e-07 * exp(-v / 7.7)
!F a13 = 3.802 / (0.1027 * exp(-v / 12.0) + 0.25 * exp(-v / 150.0))
!F v = piecewise(vtest,5.0<time<=300.0,vhold)