function [L, U]=my_lu(A)
    % trova i fattori LU della fattorizzazione LU di A senza pivoting
    %
    [n, m]=size(A);
    L=eye(n);
    U=A;
    for k=1:n-1
        pivot=U(k, k);
        L(k+1:n, k)=U(k+1:n, k)/pivot;
        colonnaL= L(k+1:n, k);
        rigaU=U(k, k+1:n);
        U(k+1:n, k+1:n)= U(k+1:n, k+1:n)-colonnaL*rigaU;
        U(k+1:n, k)=0;
    end
end