Function rk4_step 

pub fn rk4_step<F>(state: &mut PhysicalState, step_size: f64, derivative: F)
where
    F: Fn(f64, &[f64]) -> Vec<f64>,

Runge-Kutta 4th Order ODE Solver Solves dy/dt = f(t, y)