Create ilprog.m

This function is an ILP solver. It uses branch and bound method to solve the ILP. MATLAB linprog is used for LP relaxation