clear;

A = rand(5,5);
n = size(A,1); % Groesse der Matrix

L = eye(n); % Wir starten mit einer Einheitsmatrix fuer L
U = zeros(n,n); % und mit einer Nullmatrix für U.

for i = 1:n
    for k = i:n
        s = 0;
        for j = 1:i-1
            s = s + L(i,j)*U(j,k);
        end
        U(i,k) = A(i,k) - s;
    end
    
    for k = i+1:n
        s = 0;
        for j = 1:i-1
            s = s + L(k,j)*U(j,i);
        end
        L(k,i) = (A(k,i) - s)/U(i,i);
    end
    
end



