Implementazione in Matlab risolvere $LUx=b$

Il codice che risolve il problema suddetto è il seguente:

function x=solveLU(A,b)
%SOLVELU Risolve il sistema lineare Ax=b fattorizzando la matrice
%   A come LU ed infine risolvendo i sistemi
%                Ly=b
%                Ux=y
%
%   x=SOLVELU(A,b)
%
%   I parametri della funzione sono:
%       A -> la matrice dei coefficienti del sistema lineare
%       b -> il vettore dei termini noti
%
%   I valori di ritorno sono:
%       x -> il vettore soluzione del sistema lineare
%
%   See Also FATTLU
  A=fattLU(A);
  n=length(b);
  x1=b;
  for i=2:n
     for j=1:i-1
         x1(i)=x1(i)-A(i,j)*x1(j);
     end
  end
  x=solveUT(A,x1);