use eqsolver::ODESolver;

fn main() {
    // Want to solve: y' = f(t, y) = t*y, starting at t = 0 and ending at t = 2, know that x(0) = 0.2
    let f = |t: f64, y: f64| t * y;
    let (x0, y0) = (0., 0.2);
    let x_end = 2.;
    let step_size = 1e-3;

    let solver = ODESolver::new(f, x0, y0, step_size);
    let solution = solver.solve(x_end).unwrap();

    println!("Solution: {solution}");
}