% Tests LU factorization with partial pivoting
% B&F 6.5.5a
A = [2 -1 1
     3  3  9
     3  3  5]
[L,U,P]=mylu(A)
fprintf('error=%g\n',norm(P*A - L*U,'inf'));

% B&F 6.5.3c
A = [  1  1 -1  0
       1  1  4  3
       2 -1  2  4
       2 -1  2  3 ]
[L,U,P]=mylu(A)
fprintf('error=%g\n',norm(P*A - L*U,'inf'));


% Sample output:
%  >> test_mylu
%  
%  A =
%  
%       2    -1     1
%       3     3     9
%       3     3     5
%  
%  
%  L =
%  
%      1.0000         0         0
%      0.6667    1.0000         0
%      1.0000         0    1.0000
%  
%  
%  U =
%  
%       3     3     9
%       0    -3    -5
%       0     0    -4
%  
%  
%  P =
%  
%       0     1     0
%       1     0     0
%       0     0     1
%  
%  error=0
%  
%  A =
%  
%       1     1    -1     0
%       1     1     4     3
%       2    -1     2     4
%       2    -1     2     3
%  
%  
%  L =
%  
%      1.0000         0         0         0
%      0.5000    1.0000         0         0
%      0.5000    1.0000    1.0000         0
%      1.0000         0         0    1.0000
%  
%  
%  U =
%  
%      2.0000   -1.0000    2.0000    4.0000
%           0    1.5000    3.0000    1.0000
%           0         0   -5.0000   -3.0000
%           0         0         0   -1.0000
%  
%  
%  P =
%  
%       0     0     1     0
%       0     1     0     0
%       1     0     0     0
%       0     0     0     1
%  
%  error=0

