Ref: https://discourse.julialang.org/t/how-to-create-secondorderodeproblem-from-odesystem/95847 It would be nice to be able to build the following 3 systems (from [this section](https://docs.sciml.ai/DiffEqDocs/latest/solvers/dynamical_solve/) of docs): ``` DynamicalODEProblem{isinplace}(f1, f2, v0, u0, tspan, p = NullParameters(); kwargs...) SecondOrderODEProblem{isinplace}(f, du0, u0, tspan, p = NullParameters(); kwargs...) HamiltonianProblem{T}(H, p0, q0, tspan, p = NullParameters(); kwargs...) ``` Directly from an ODESystem.