k1 = beta
k2 = gamma
k3 = beta0
k4 = mu
k5 = a
k6 = b
k7 = n
k8 = alpha
k9 = alpha0
k10 = Q0
k11 = delta_0
k12 = theta_0
k13 = delta_1
k14 = k1
k15 = n1
k16 = compartment