Function euler_ode 

pub fn euler_ode<F>(
    f: F,
    t0: f64,
    y0: Vec<f64>,
    t_end: f64,
    dt: f64,
) -> Vec<(f64, Vec<f64>)>
where
    F: Fn(f64, &[f64]) -> Vec<f64>,

Euler’s method for ODEs: dy/dt = f(t, y).

Returns (time, state) pairs.