Y = lu(A)

function Get_LU_Matrix
    A = [1,2,2; 4,4,2; 4,6,4];
    Y = A;
    n = size(A,2);

    for k=1:n-1
        if A(k, k) == 0
           break
        end

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

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

    disp(Y);
end