function[xfvalexitflag]=intprog（fAbAeqbeqlbubMTolXInteger）
options=optimset（‘display‘‘off‘）;
bound=inf;
[x0fval0]=linprog（fAbAeqbeqlbub[]options）;
[xfvalexitflagb]=rec_BranchBound（fAbAeqbeqlbubx0fval0MTolXIntegerbound）;
function[xxfvalexitflagbb]=rec_BranchBound（fAbAeqbeqlbubxvMTolXIntegerbound）
options=optimset（‘display‘‘off‘）;
[x0fval0exitflag0]=linprog（fAbAeqbeqlbub[]options）;
if?exitflag0<=0||fval0>bound
????xx=x;
????fval=v;
????exitflag=exitflag0;
????bb=bound;
????return;
end
ind=find（abs（x0（M）-round（x0（M）））>TolXInteger）;
if?isempty（ind）
????exitflag=1;
????if?fval0????????x0（M）=round（x0（M））;
????????xx=x0;
????????fval=fval0;
????????bb=fval0;
????else
????????xx=x;
????????fval=v;
????????bb=bound;