% Musterlösung
%
% Numerische Mathematik, SS 2013
%
% LR_Zerlegung(zeilenorientiert)
%
% Ingolf Busch, Juni 2013

% Input: A     quadratische Matrix


% Output: S  struct mit allen Informationen 
%         S = lr(A)
%
%         L linke untere normierte Dreiecksmatrix
%         R rechte obere Dreiecksmatrix
%         [L,R] = lr(A)


function [L,R] = my_lr(A)

n=size(A,1);

for k=1:(n-1)
    for i = (k+1):n
        A(i,k)       = A(i,k)/A(k,k);
        A(i,(k+1):n) = A(i,(k+1):n)-A(i,k)*A(k,(k+1):n);
    end
end

if(nargout==1)
    L.mode = 'zeile';
    L.L    = eye(n)+tril(A,-1);
    L.R    = triu(A);
else
    L    = eye(n)+tril(A,-1);
    R    = triu(A);
end