from gekko import GEKKO m = GEKKO(remote=False) c = [0,1] A = [[-1,1],[3,2],[2,3]] b = [1,12,12] z = m.Array(m.Var,2,integer=True,lb=0) m.qobj(c,x=z,otype='max') m.axb(A,b,x=z,etype='<=') m.options.SOLVER = 1 m.solve() print('Objective: ', -m.options.OBJFCNVAL) print('x: ', z[0].value[0]) print('y: ', z[1].value[0])