function [L,U] = min_lu(A)
%
%  [L,U]=minlu(A) utfoerer en LU-faktorisering
%        uten pivotering av matrisen A.
%
[n,m] = size(A); 
if n~=m 
    error('matrisen er ikke kvadratisk');
end
L = zeros(n,n);   
U = zeros(n,n);

for k=1:n
    L(k,k)=1; %Doolittle factorization.
    U(k,k) = (A(k,k)-L(k,1:k-1)*U(1:k-1,k))/L(k,k);
    for i=k:n
        L(i,k) = (A(i,k)-L(i,1:k-1)*U(1:k-1,k))/U(k,k);
        U(k,i) = (A(k,i)-L(k,1:k-1)*U(1:k-1,i))/L(k,k);
    end
end
        


