from gekko import GEKKO
m = GEKKO(remote=False)
x1,x2,x3,Z = m.Array(m.Var,4,integer=True)
m.Minimize(Z)
m.Equation(x1+x2+x3==16)
m.Equations([Z>=x1,Z>=x2,Z>=x3])
m.options.solver=1
m.solve()
print('x1: ',x1.value[0])
print('x2: ',x2.value[0])
print('x3: ',x3.value[0])
print('Z:  ',Z.value[0])