pub fn euler( f: impl Fn(f64, &[f64]) -> Vec<f64>, t0: f64, y0: &[f64], dt: f64, steps: usize, ) -> Vec<Vec<f64>>
Forward Euler method. Returns list of state vectors at each step.