k1 = d k2 = beta k3 = rho k4 = psi_1 k5 = delta k6 = mu_1 k7 = s_1 k8 = s_2 k9 = mu_2 k10 = s_1_prime k11 = s_2_prime k12 = psi_2 k13 = p_1 k14 = delta_1 k15 = p_2 k16 = p_3 k17 = delta_2 k18 = q_1 k19 = q_2 k20 = r_e k21 = alpha k22 = p k23 = s_3 k24 = s_4 k25 = c k26 = k k27 = d_a k28 = q k29 = compartment