
function x=metodo_LU(A,b)

[L,U,P]=lu(A); % fattorizzazione PA=LU
c=P*b; 
y=L\c; % sistema triangolare inferiore
x=U\y; % sistema triangolare superiore

