pub fn rk45(
f: impl Fn(f64, &[f64]) -> Vec<f64>,
t0: f64,
y0: &[f64],
t_end: f64,
tol: f64,
h_min: f64,
h_max: f64,
) -> (Vec<f64>, Vec<Vec<f64>>)Expand description
Dormand-Prince RK45 adaptive step integrator. Returns (time_points, state_vectors).