%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%   my scripts for bits and pieces of the Tridiagonal demonstration
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function  seeLU(A)
% SEELU - just shows side by side plots of a LU decomposition of A
%   
  [L,U]=lu(A);
  figure
  subplot(1,2,1); spy(L); title('L');
  subplot(1,2,2); spy(U); title('U');
  disp(sprintf('\n Here are the first 5 rows of L and U\n'));
  L=full(L(1:5,1:5))
  U=full(U(1:5,1:5))
  return