Function three_compartments

pub fn three_compartments(
    x: &Matrix<f64, Dyn, Const<1>, VecStorage<f64, Dyn, Const<1>>>,
    p: &Matrix<f64, Dyn, Const<1>, VecStorage<f64, Dyn, Const<1>>>,
    t: f64,
    rateiv: Matrix<f64, Dyn, Const<1>, VecStorage<f64, Dyn, Const<1>>>,
    _cov: &Covariates,
) -> Matrix<f64, Dyn, Const<1>, VecStorage<f64, Dyn, Const<1>>>

Analytical for three compartments Assumptions:

• p is a vector of length 5 with k10, k12, k13, k21 and k31 in that order. k10 represents the elimination rate constant from compartment 1 (central). The remaining parameters of the form kij represent the transfer rate constants from the ith to the jth compartment
• rateiv is a vector of length 1 with the value of the infusion rate (only one drug)
• x is a vector of length 3
• covariates are not used