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

% the matrix from exercise 1.8.2
A = [ 0 4 1 ; 1 1 3 ; 2 -2 1 ]; 
b = [ 9 6 -1 ]'; 

[L,U,P] = lu(A); 
P, L, U 
P' * L * U 

[K,V] = lu(A)
K, V
K*V