m = GEKKO(remote=False) x = m.Array(m.Var,2,lb=0,ub=10,integer=True) m.Minimize(c@x) Ax = A_u@x m.Equations([Ax[i]<=b_u[i] for i in range(len(b_u))]) m.options.SOLVER=1 m.solve(disp=False) #Print the results print("Objective function: ", m.options.objfcnval) print("Optimal solution: ", x)