%
% 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 = [-4,-5,0,0,0];
a = [ 1 1 1 0 0; 1 2 0 1 0; 1 4 0 0 1 ];
b = [5;7;10];
[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) )