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