k1 = beta_x k2 = psi k3 = alpha_x k4 = beta_y k5 = alpha_y k6 = alpha_0 k7 = k k8 = alpha_k k9 = compartment