#tyson9 FIX: S #reactions r0: {1.0}X = {1.0}R env * (k0_prime + k0 * Ep) * X r1: $pool = {1.0}X env * k1 * S r2: {1.0}R = $pool env * k2 * R r3: $pool = $pool env * k3 * R * E / (Km3 + E) r4: $pool = $pool env * k4 * Ep / (Km4 + Ep) #parameters Et = 1.0 J3 = 0.05 J4 = 0.05 k0 = 0.4 k0_prime = 0.01 k1 = 1.0 k2 = 1.0 k3 = 1.0 k4 = 0.3 S = 0.0 env = 1.0 #initial values R = 0.0 X = 0.0 #assignment rules !F E = Et - Ep !F Ep = goldbeter_koshland(k3 * R, k4, J3, J4) * Et !F Km4 = J4 * Et !F Km3 = J3 * Et