MTK and JuMP

[ChrisRackauckas]

There’s a way to do this with an arbitrary differential equation stepper:

https://mtk.sciml.ai/dev/tutorials/nonlinear_optimal_control/

It’s on the ControlSystem but easy to map over the ODESystem.