Function aumc 

pub fn aumc(
    times: &[f64],
    values: &[f64],
    method: &AUCMethod,
) -> Result<f64, ObservationError>

∫ t·C(t) dt from first to last time point (for MRT = AUMC / AUC).