#yang



#reactions
R1:
 $pool = {1.0}x1
  cell * K15 * x15 * lin * (1 + x4 / KI19 + x2 / KI20 + x13 / KI21 + x11 / KI22) / (lin + k15 * (1 + x1 / ks))
R10:
 {1.0}x1 = {1.0}x10
  cell * K21 * x21 * x1 / (x1 + k21 * (1 + x5 / ki7 + x3 / ki8 + x7 / ki11 + x11 / ki12 + x10 / ks))
R11:
 {1.0}x10 = {1.0}x11
  cell * K24 * x24 * x10 / (x10 + k24 * (1 + x11 / ks))
R12:
 {1.0}x10 = {1.0}x12
  cell * K21 * x21 * x10 / (x10 + k21 * (1 + x5 / ki7 + x3 / ki8 + x7 / ki11 + x11 / ki12 + x12 / ks))
R13:
 {1.0}x12 = {1.0}x13
  cell * K22 * x22 * x12 / (x12 + k22 * (1 + x13 / ks))
R14:
 {1.0}x13 = {1.0}x14
  cell * K23 * x23 * x13 / (x13 + k23 * (1 + x5 / ki14 + x11 / ki15 + x14 / ks))
R16:
 $pool = {1.0}x16
  cell * a24 * x7 * x7 / (x7 * x7 + KI24 * KI24)
R17:
 {1.0}x17 = $pool
  cell * ki17 * x2 * x17
R18:
 {1.0}x9 = $pool
  kd9 * x9 * cell
R19:
 {1.0}x20 = $pool
  cell * ki5 * x6 * x20
R2:
 {1.0}x1 = {1.0}x2
  cell * K16 * x16 * x1 / (x1 + k16 * (1 + x2 / ks))
R20:
 {1.0}x20 = $pool
  ki4 * x2 * x20 * cell
R21:
 $pool = {1.0}x21
  cell * KI23 * x13 * x21
R22:
 {1.0}x22 = $pool
  cell * K22 * x22 * x12 / ((x12 + k22) * 129)
R23:
 {1.0}x12 = $pool
  cell * x12 * kd12
R24:
 {1.0}x13 = $pool
  cell * kd13 * x13
R25:
 {1.0}x8 = $pool
  cell * kd8 * x8
R26:
 {1.0}x3 = $pool
  kd3 * x3 * cell
R27:
 {1.0}x2 = $pool
  kd2 * cell * x2
R28:
 {1.0}x16 = $pool
  cell * kd16 * x16
R29:
 {1.0}x11 = $pool
  kd11 * x11 * cell
R3:
 {1.0}x2 = {1.0}x3
  cell * K24 * x24 * x2 / (x2 + k24 * (1 + x3 / ks))
R30:
 {1.0}x21 = $pool
  cell * ki9 * x12 * x21
R31:
 {1.0}x21 = $pool
  cell * ki10 * x10 * x21
R32:
 {1.0}x21 = $pool
  cell * ki6 * x21 * x2
R34:
 {1.0}x1 = $pool
  0.1 * x1 * cell
R4:
 {1.0}x1 = {1.0}x4
  cell * K17 * x17 * x1 / (x1 + k17 * (1 + x4 / ki18 + x3 / ki16 + x4 / ks))
R5:
 {1.0}x4 = {1.0}x5
  cell * K24 * x24 * x4 / (x4 + k24 * (1 + x5 / ks))
R6:
 {1.0}x1 = {1.0}x6
  cell * K18 * x18 * x1 / (x1 + k18 * (1 + x7 / ki3 + x6 / ks))
R7:
 {1.0}x6 = {1.0}x7
  cell * K19 * x19 * x6 / (x6 + k19 * (1 + x1 / ki1 + x3 / ki2 + x7 / ks))
R8:
 {1.0}x6 = {1.0}x8
  cell * K20 * x20 * x6 / (x6 + k20 * (1 + x8 / ks))
R9:
 $pool = {1.0}x9
  kd8 * x8 * cell


#parameters
K15 = 3600.0
K16 = 1000.0
K17 = 1000.0
K18 = 1000.0
K19 = 3000.0
K20 = 1599.0
K21 = 5000.0
K22 = 125.0
K23 = 150.0
K24 = 500.0
KI19 = 500.0
KI20 = 200.0
KI21 = 500.0
KI22 = 500.0
KI23 = 0.053
KI24 = 2.3e-05
a24 = 0.15
k15 = 2600.0
k16 = 70.0
k17 = 50.0
k18 = 50.0
k19 = 160.0
k20 = 4.0
k21 = 5.0
k22 = 20.0
k23 = 3.9
k24 = 70.0
kd11 = 0.001
kd12 = 0.07
kd13 = 0.01
kd16 = 0.01
kd2 = 0.05
kd3 = 0.01
kd8 = 0.1
kd9 = 0.001
ki1 = 0.3
ki10 = 0.01
ki11 = 15.0
ki12 = 6.3
ki14 = 0.2
ki15 = 0.86
ki16 = 10.0
ki17 = 10.0
ki18 = 10.0
ki2 = 30.0
ki3 = 30.0
ki4 = 0.6
ki5 = 0.1
ki6 = 0.01
ki7 = 30.0
ki8 = 4.0
ki9 = 0.175
ks = 500.0
lin = 12.0
cell = 1.0


#initial values
x1 = 0.001
x10 = 0.001
x11 = 0.001
x12 = 0.001
x13 = 0.001
x14 = 0.001
x15 = 1.5
x16 = 1.5
x17 = 0.5
x18 = 0.2
x19 = 0.5
x2 = 0.001
x20 = 0.5
x21 = 5.0
x22 = 0.76
x23 = 0.07
x24 = 0.8
x25 = 0.0
x3 = 0.001
x4 = 0.001
x5 = 0.001
x6 = 0.001
x7 = 0.001
x8 = 0.001
x9 = 0.001


#assignment rules