function [L, U]=simaGauss(A)

U=A;
L=eye(size(A,1));

for i=1:size(A,2)-1
    if A(i,i)==0
        error('A főátlóba 0 került!');
    end
    
    gausstr=eye(size(A,1));
    gausstr(i+1:end,i)= -U(i+1:end,i)./U(i,i);
    U=gausstr*U;
    L(i+1:end,i)=-gausstr(i+1:end,i);
end
