function [L,U]=myLU(A)
%% Funkce implementujici LU rozklad matice, bez pivotace
%% "Intuitivni implementace" - takto by to clovek nejspis pocital rucne

n=length(A);

L=eye(n); %inicializace diagonaly jednickami

U=A; %matici A ulozim do U a pak uz delam upravy primo v U

    for k=1:n-1   % iteruji pres vsechny sloupecky  
        for j=k+1:n % iteruji v kazdem sloupecku pres poddiagonalni prvky
            
            m=U(j,k)/U(k,k); % prvek m_{j,k}   
            
            for l=k:n %% pricteni m-nasobku k-teho radku k j-temu radku
                U(j,l)=U(j,l)-m*U(k,l); 
            end
            
            L(j,k)=m; %% ulozim koeficient m_{j,k} na spravne misto do L
        end
    end
end