Function is_singular 

pub fn is_singular<F>(matrix: ArrayView2<'_, F>, threshold: Option<F>) -> bool
where
    F: IntegrateFloat,

Check if a matrix is singular or ill-conditioned

Uses condition number estimation to check if a matrix is close to singular, which would cause problems for ODE solvers