pub trait ODEFunction { // Required method fn derivatives(&self, t: f64, y: &[f64; 4]) -> [f64; 4]; }
ODE function trait for RK4
Calculate derivatives dy/dt = f(t, y)