I have been working on an optimal control problem using Pyomo + IPOPT, but I just cannot find an optimal solution. It is a 3 phase trajectory problem with an objective of minimizing mass. I have formulated the objective by specifying the that the final mass of the third phase is maximized: