Function adaptive_jacobian 

pub fn adaptive_jacobian<F, Func>(
    f: &Func,
    t: F,
    y: &Array1<F>,
    f_current: &Array1<F>,
    perturbation_scale: F,
) -> IntegrateResult<Array2<F>>
where
    F: IntegrateFloat,
    Func: Fn(F, ArrayView1<'_, F>) -> Array1<F> + Clone,

Jacobian strategy that uses autodiff when available and falls back to finite differences when not