Function solve_initial_condition_sensitivity_with 

pub fn solve_initial_condition_sensitivity_with<Sol, S, F>(
    rhs: F,
    y0: &[S],
    t0: S,
    tf: S,
    options: &SolverOptions<S>,
) -> Result<StateTransitionResult<S>, SolverError>
where
    Sol: Solver<S>,
    S: Scalar,
    F: Fn(S, &[S], &mut [S]),

Closure-shaped convenience wrapper around solve_initial_condition_sensitivity.

Wraps rhs(t, y, dydt) into an internal OdeProblem (forward-FD Jacobian via the OdeSystem default) and forwards. Suitable for one-shot analyses and REPL-style scripts. For an analytical J_y or an autodiff-derived J_y, implement OdeSystem directly (with a jacobian override or by wrapping in [crate::AutodiffJacobianSystem]) and call solve_initial_condition_sensitivity.