pub fn rk4_solve<F>( f: &F, t0: f64, y0: &[f64], t1: f64, steps: usize, ) -> Vec<(f64, Vec<f64>)>where F: Fn(f64, &[f64]) -> Vec<f64>,
Full trajectory [(t, y)], including the initial point, over steps RK4 steps.
[(t, y)]
steps