options = cplexoptimset('cplex');
options.Display = 'on';
options.mip.tolerances.mipgap =  0;
options.parameter2058 = 1;
options.mip.pool.relgap= 0;
options.parameter2104 = 2;
[x, fval, exitflag, output] = cplexbilp (f, Aineq, bineq, Aeq, beq, ...
      [ ], options);