Function rk4 

pub fn rk4(
    f: impl Fn(f64, &[f64]) -> Vec<f64>,
    t0: f64,
    y0: &[f64],
    dt: f64,
    steps: usize,
) -> Vec<Vec<f64>>

Classical 4th-order Runge-Kutta.