matlab线性规划取整,MATLAB求整数线性规划问题

为什么下面的这个程序会解不出来线性规划解呢？求正确的程序代码，感谢大家。

function [x,y]=ILp(f,G,h,Geq,heq,lb,ub,x,id,options)

%整数线性规划分支定界法，可求解纯整数规划和混合整数规划。

%y=minf’*x  s.t.  G*x<=h  Geq*x=heq  x为全整数或混合整数列向量

%用法

%[x,y]=ILp(f,G,h,Geq,heq,lb,ub,x,id,options)

%参数说明

%lb:解的下界列向量(Default:-int)

%ub:解的上界列向量(Default:int)

%x:迭代初值列向量

%id：整数变量指标列向量，1-整数，0-实数(Default:1)

global upper opt c x0 A b Aeq beq ID options;

if nargin<10,options=optimset({});options.Display='off';

options.LargeScale='off';end

if nargin<9,id=ones(size(f));end

if nargin<8,x=[];end

if nargin<7 |isempty(ub),ub=inf*ones(size(f));end

if nargin<6 |isempty(lb),lb=zeros(size(f));end

if nargin<5,heq=[];end

if nargin<4,Geq=[];end