% This set of function implements MILP1 and MILP3 from the article
% 'Computing elementary flux modes involving a set of target reactions'
% by Laszlo David and Alexander Bockmayr
% (http://nbn-resolving.de/urn:nbn:de:0296-matheon-12621)
% See twoStep_EM.m for details on usage.