I tried to solve a 2D optimal lunar soft landing problem by GEKKO.
The state and control variables as well as equations are listed as follows
