num::rk4_2nd

Defined in header <ode/ode.hpp>

template <typename Accel = accel_fn, typename State = vec<real>>
requires vector_space<State> && std::invocable<Accel &, const State &, State &>
basic_rk4_2nd_steps<Accel, State> rk4_2nd(Accel accel, State q0, State v0,
                                          ode_params p = {});

Lazy fourth-order Nystrom trajectory for \(\ddot{q} = a(q)\).