#alexander2

FIX: A_im

#reactions
r10:
 {1.0}A = $pool
  b1 * A
r11:
 {1.0}A = $pool
  sigma1 * A * R
r1a:
 {1.0}G = $pool
  v_max / (k + G) * G
r1b:
 {1.0}A_im = {1.0}A
  f * (v_max / (k + G)) * G
r2:
 $pool = {1.0}G
  gamma * E
r3:
 $pool = {1.0}R
  beta * A
r4:
 $pool = {1.0}R
  pi1 * E * A
r5:
 $pool = {1.0}E
  lambdaE * A
r6:
 {1.0}A = $pool
  muA * A
r7:
 {1.0}R = $pool
  muR * R
r8:
 {1.0}E = $pool
  muE * E
r9:
 {1.0}G = $pool
  muG * G


#parameters
b1 = 0.25
beta = 200.0
f = 0.0001
gamma = 2000.0
k = 50000000.0
lambdaE = 1000.0
muA = 0.25
muE = 0.25
muG = 5.0
muR = 0.25
pi1 = 0.016
sigma1 = 3e-06
v_max = 125000.0
A_im = 0.0
body = 1.0


#initial values
A = 1.0
E = 0.0
G = 100000000.0
R = 0.0


#assignment rules
!F R0 = f * (v_max / k) * lambdaE * gamma / (mG * mA * muE)
!F mG = muG + v_max / k
!F mA = b1 + muA