pub fn estimate_initial_step<F, Func>( f: &Func, t: F, y: &Array1<F>, dy: &Array1<F>, tol: F, tend: F, ) -> Fwhere F: IntegrateFloat, Func: Fn(F, ArrayView1<'_, F>) -> Array1<F>,
Calculate a safe step size based on function derivatives