function LR=LRdecomp(A,how)

n=size(A,1);
if how==0
    % zeilenorientierte LR-Zerlegung
    for i=1:n-1
        for j=i+1:n
            A(j,i)=A(j,i)/A(i,i);
            A(j,i+1:n)=A(j,i+1:n)-A(j,i)*A(i,i+1:n);
        end
    end
else
    % spaltenorientierte LR-Zerlegung
    for i=1:n-1
        A(i+1:n,i)=A(i+1:n,i)/A(i,i);
        for k=i+1:n
            A(i+1:n,k)=A(i+1:n,k)-A(i,k)*A(i+1:n,i);
        end
    end
end

LR=A;