#solve NLP m = GEKKO(remote=False) x,y = m.Array(m.Var,2,lb=1,ub=3) y.LOWER=2 m.Minimize(x*y**2-x**2-y**2) m.Equations([x+y>=4,x*y<=5]) m.solve(disp=False) xopt = x.value[0]; yopt = y.value[0] print('x:', xopt,'y:', yopt,'obj:',-m.options.objfcnval)