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

[L,U] = slu(pascal(4))
[L,U] = slu(pascal(5))

% The pattern is 
% 1) U=L^T
% 2) U starts with a row of all ones and a main diagonal of all ones
% 3) the (i,j)th element of U is computed by U(i,j) = U(i,j-1)+U(i-1,j-1)
% 4) much like pascal's triangle