#ataullahkhanov1 #reactions U1: $pool = {1.0}I cell * P * J U2: {3.0}I + {1.0}E = $pool cell * W2 * I * T U3: $pool = {1.0}E cell * W3 * pow(T, 0.52) * pow(M, 0.41) U6_plus_3U7: {1.0}E = $pool cell * 2 * U U6_plus_U7_minus_U8: $pool = {1.0}A cell * U * (1 - W * pow(T, n) * pow(M, k)) #parameters J = 100.0 P = 0.121 U = 0.02 W = 0.01 W2 = 0.2 W3 = 13.48 k = -1.0 n = 1.2 cell = 1.0 #initial values A = 1.11 E = 2.1 I = 10.0 #assignment rules !F T = (A + 3.0 * E - pow(6.0 * A * E - 3.0 * pow(E, 2.0) + pow(A, 2.0), 0.5)) / 6.0 !F M = (7.0 * A - 3.0 * E - pow(6.0 * A * E - 3.0 * pow(E, 2.0) + pow(A, 2.0), 0.5)) / 6.0