#kim4



#reactions
React1:
 $pool = {1.0}R1
  i1
React10:
 {1.0}R2 = $pool
  d_R2 * R2
React11:
 $pool = {1.0}P2
  sP2R2 * R2
React12:
 {1.0}P2 = $pool
  d_P2 * P2
React13:
 $pool = {1.0}P2_prime
  sP2_prime_P2 * P2
React14:
 $pool = {1.0}P2_prime
  s2 / (1 + pow(P1_prime, n))
React15:
 {1.0}P2_prime = $pool
  d_P2_prime * P2_prime
React16:
 $pool = {1.0}P3_prime
  sP3_prime_P2_prime * (pow(P2_prime, n) / (1 + pow(P2_prime, n)))
React17:
 {1.0}P3_prime = $pool
  d_P3_prime * P3_prime
React18:
 $pool = {1.0}P2_prime
  s3 / (1 + pow(P4_prime, n))
React19:
 $pool = {1.0}P4_prime
  sP4_prime_P1_prime * (pow(P1_prime, n) / (1 + pow(P1_prime, n)))
React2:
 {1.0}R1 = $pool
  d_R1 * R1
React20:
 {1.0}P4_prime = $pool
  d_P4_prime * P4_prime
React3:
 $pool = {1.0}P1
  sP1R1 * R1
React4:
 {1.0}P1 = $pool
  d_P1 * P1
React5:
 $pool = {1.0}P1_prime
  sP1_prime_P1 * P1
React6:
 $pool = {1.0}P1_prime
  s1 / (1 + pow(P2_prime, n))
React7:
 $pool = {1.0}P1_prime
  s3 / (1 + pow(P3_prime, n))
React8:
 {1.0}P1_prime = $pool
  d_P1_prime * P1_prime
React9:
 $pool = {1.0}R2
  i2


#parameters
d_P1 = 0.22367
d_P1_prime = 0.37048
d_P2 = 0.22367
d_P2_prime = 0.37048
d_P3_prime = 0.37048
d_P4_prime = 0.37048
d_R1 = 0.23521
d_R2 = 0.23521
n = 9.0
s1 = 0.4
s2 = 0.3
s3 = 0.2
sP1R1 = 0.47305
sP1_prime_P1 = 0.28687
sP2R2 = 0.47305
sP2_prime_P2 = 0.28687
sP3_prime_P2_prime = 0.5
sP4_prime_P1_prime = 0.5
compartment = 1.0


#initial values
P1 = 0.1
P1_prime = 0.1
P2 = 1.0
P2_prime = 1.0
P3_prime = 0.1
P4_prime = 0.1
R1 = 0.01
R2 = 0.1


#assignment rules
!F i1 = piecewise(1.0, and(geq(time, 50.0), leq(time, 100.0)), 0.0)
!F i2 = piecewise(1.0, and(geq(time, 175.0), leq(time, 225.0)), 0.0)