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

% Create an A (example on page 28 of the book):
A = [ 5 0 0 ; 2 -4 0 ; 1 2 3 ]; 
% Create a b: 
b = [ 15 -2 10 ].'; 
% the known truth: 
y_truth = [ 3 2 1 ].';

y_approx = rco_forward_sub(A,b); 

norm( y_truth-y_approx )