Skip to content

num::ode_euler

Defined in header <ode/ode.hpp>

template <typename RHS = ode_rhs_fn, typename State = vec<real>>
requires vector_space<State> && std::invocable<RHS &, real, const State &, State &>
ode_result ode_euler(RHS f, State y0, ode_params p = {},
                     const observer_fn &observer = {});

Integrate \(\dot{y} = f(t, y)\) with fixed-step forward Euler.

observer, if given, is called as observer(t, u) after each step.