Function psi 

pub fn psi(
    equation: &impl Equation,
    subjects: &Data,
    support_points: &ArrayBase<OwnedRepr<f64>, Dim<[usize; 2]>>,
    error_models: &ErrorModels,
    progress: bool,
    cache: bool,
) -> Result<ArrayBase<OwnedRepr<f64>, Dim<[usize; 2]>>, PharmsolError>

Calculate the psi matrix for maximum likelihood estimation.

Parameters

• equation: The equation to use for simulation
• subjects: The subject data
• support_points: The support points to evaluate
• error_model: The error model to use
• progress: Whether to show a progress bar
• cache: Whether to use caching

Returns

A 2D array of likelihoods