Untitled

function [L U]=LUFactorization(A)
[nr nc]=size(A);
L=eye(nr);
if nr==nc %allora possiamo effettuare i calcoli
for i=1:nr
    for j=i+1:nr
        if A(i,i)==0
            disp(['errore: pivot nullo!'])
        else
            L(j,i)=A(j,i)/A(i,i);
            for k=i:nr
                A(j,k)=A(j,k)-L(j,i)*A(i,k)
            end
        end
    end
end
else
    disp(['errore: la matrice dei coefficienti non è quadrata'])
end
U=A;
end