k1 = mu
k2 = l_e
k3 = beta_1
k4 = R0_1
k5 = gamma_1
k6 = beta_2
k7 = R0_2
k8 = gamma_2
k9 = tInf_1
k10 = tInf_2
k11 = sigma
k12 = tImm
k13 = Lambda_1
k14 = Lambda_2
k15 = I_1_frac
k16 = I_2_frac
k17 = S_frac
k18 = R1_frac
k19 = R2_frac
k20 = Rp_frac
k21 = psi
k22 = env