peroxide::numerical::runge_kutta
pub fn rk4<F>( init_param: f64, init_value: Vec<f64>, f: F, step: f64, num: usize) -> Matrix where F: Fn(f64, Vec<f64>) -> Vec<f64> + Copy,