Function finite_difference_jacobian 

pub fn finite_difference_jacobian<F, Func>(
    f: &Func,
    t: F,
    y: &Array1<F>,
    f_eval: &Array1<F>,
    _perturbation_scale: F,
) -> Array2<F>
where
    F: IntegrateFloat,
    Func: Fn(F, ArrayView1<'_, F>) -> Array1<F>,

Calculate finite difference approximation of the jacobian matrix