#tyson13

FIX: S

#reactions
r1:
 $pool = {1.0}R
  env * k0 * Ep
r2:
 $pool = {1.0}R
  env * k1 * S
r3:
 {1.0}R = $pool
  env * k2 * R
r4:
 $pool = $pool
  env * k3 * R * E / (J3 + E)
r5:
 $pool = $pool
  env * k4 * Ep / (J4 + Ep)


#parameters
Et = 1.0
J3 = 0.05
J4 = 0.05
k0 = 0.4
k1 = 0.01
k2 = 1.0
k3 = 1.0
k4 = 0.2
S = 0.0
env = 1.0


#initial values
R = 0.0


#assignment rules
!F E = Et - Ep
!F Ep = goldbeter_koshland(k3 * R, k4, J3, J4)