%
% Written by:
% --
% John L. Weatherwax                2009-04-21
%
% email: wax@alum.mit.edu
%
% Please send comments and especially bug reports to the
% above email address.
%
%-----

close all; clc; clear all; 

addpath('../../BookCode');

c = [5,7,10,0,0];
a = [ 1 1 1 -1 0; 1 2 4 0 -1 ];
b = [4;5];
[xsol,ind] = barnes(a,b,c,1.e-6);
i=1; 
x_full = zeros( size(a,2), 1 ); 
for j=ind
  x_full(j) = xsol(i);
  i = i+1;
end

fprintf('Full solution is:\n');
for i=1:size(a,2), 
  fprintf('x(%1.0f)=%8.4f\n', i, x_full(i));
end
fprintf('objective function= %10.6f\n', dot(c,x_full) )