function a = mylu_nop(a) 
%
%  MYLU_NOP  performs LU decompostion of A 
%  *without* pivoting, the resulting L and 
%  U returned 'logically' in A (L a unit lower 
%  tridiag etc.) 
%
n = size(a,1);                       % determine matrix size
for k = 1:n-1
    for i = k+1:n                      % row loop for elimination
        m = a(i,k)/a(k,k);             % set multiplier
        a(i,k) = m;                    % store multiplier for L
        for j = k+1:n                  % loop over columns k+1:n       
            a(i,j) = a(i,j)-m*a(k,j);
        end
    end                                  
end
        
