#yildirim1

FIX: L_e

#reactions
r_a1:
 {1.0}A = $pool
  cell * A * (gamma_A + mu)
r_a2:
 {1.0}A = $pool
  cell * beta_A * B * (A / (K_A + A))
r_a3_l1:
 {1.0}L = {1.0}A
  cell * alpha_A * B * (L / (K_L + L))
r_b1:
 {1.0}B = $pool
  cell * B * (gamma_B + mu)
r_b2_i2:
 {1.0}I2 = {1.0}B
  cell * I2 / tau_B
r_i1:
 $pool = {1.0}I1
  cell * alpha_M * ((K_1 * exp(mu * tau_M * -2) * pow(A, 2) + 1) / (K + K_1 * exp(-2 * mu * tau_M) * pow(A, 2)))
r_i2:
 $pool = {1.0}I2
  cell * alpha_B * M * exp(-(mu * tau_B))
r_i3:
 $pool = {1.0}I3
  cell * alpha_P * M * exp(-1 * mu * (tau_B + tau_P))
r_l2:
 {1.0}L = $pool
  cell * L * (gamma_L + mu)
r_l3:
 {1.0}L = $pool
  cell * beta_L1 * P * (L / (K_L1 + L))
r_l4:
 $pool = {1.0}L
  cell * alpha_L * P * (L_e / (K_Le + L_e))
r_m1:
 $pool = {1.0}M
  cell * gamma_0
r_m2:
 {1.0}M = $pool
  cell * M * (gamma_M + mu)
r_m3_i1:
 {1.0}I1 = {1.0}M
  cell * I1 / tau_M
r_p1:
 {1.0}P = $pool
  cell * P * (gamma_P + mu)
r_p2_i3:
 {1.0}I3 = {1.0}P
  cell * I3 / (tau_B + tau_P)


#parameters
K = 7200.0
K_1 = 25200.0
K_A = 1.95
K_L = 0.97
K_L1 = 1.81
K_Le = 0.26
alpha_A = 17600.0
alpha_B = 0.0166
alpha_L = 2880.0
alpha_M = 0.000997
alpha_P = 10.0
beta_A = 21500.0
beta_L1 = 2650.0
gamma_0 = 7.25e-07
gamma_A = 0.52
gamma_B = 0.000833
gamma_L = 0.0
gamma_M = 0.411
gamma_P = 0.65
mu = 0.0226
tau_B = 2.0
tau_M = 0.1
tau_P = 0.83
L_e = 0.08
cell = 1.0


#initial values
A = 0.038
B = 0.0
I1 = 0.0
I2 = 0.0
I3 = 0.0
L = 0.372
M = 0.000626
P = 0.0149


#assignment rules