from gekko import GEKKO
m = GEKKO(remote=False)
x1,x2,x3 = m.Array(m.Var,3,lb=0)
m.Minimize(1000-x1**2-2*x2**2-x3**2-x1*x2-x1*x3)
m.Equations([8*x1+14*x2+7*x3==56,
             x1**2+x2**2+x3**2>=25])
m.solve(disp=False)
print(x1.value,x2.value,x3.value)