Trait fast_ode::DifferentialEquation

pub trait DifferentialEquation<const N: usize> {
    fn ode_dot_y(&self, t: f64, y: &Coord<{ N }>) -> (Coord<{ N }>, bool);
}

Required methods

fn ode_dot_y(&self, t: f64, y: &Coord<{ N }>) -> (Coord<{ N }>, bool)

The differential equation is dy/dt = ode_dot_y(t, y).0 ∀ t ∈ ℝ.

If you do not like your argument, you have two choices:

If you return (Coord::<N>([NAN; N]), b) you get a smaller step, independent of b.

If you return (c, false), the solver will terminate if (t, y) is part of the solution curve. In this case, IvpResult::OdeRequestedExit is returned.

Implementors