#decroly1 #reactions r1: $pool = {1.0}alpha v_Km1 r2: {1.0}alpha = {50.0}beta sigma1 * alpha * (1 + alpha) * pow(1 + beta, 2) / (L1 + pow(1 + alpha, 2) * pow(1 + beta, 2)) r3: {1.0}beta = {0.02}gamma sigma2 * beta * (1 + d * beta) * pow(1 + gamma, 2) / (L2 + pow(1 + d * beta, 2) * pow(1 + gamma, 2)) r4: {1.0}gamma = $pool ks * gamma #parameters r1_v_Km1 = 0.45 r2_L1 = 500000000.0 r2_sigma1 = 10.0 r3_L2 = 100.0 r3_d = 0.0 r3_sigma2 = 10.0 r4_ks = 1.99 cell = 1.0 #initial values alpha = 29.19988 beta = 188.8 gamma = 0.3367 #assignment rules