Implementazione in Matlab

L'algoritmo he implementa il metodo di eliminazione di Gauss per la fattorizzazione della matrice $A$ in $LU$ è il seguente:

function A=fattLU(A)
%FATTLU Fattorizza la matrice A come LU tramite il metodo di
%   eliminazione di Gauss
%
%   A=FATTLU(A)
%
%   I parametri della funzione sono:
%       A -> la matrice quadrata da fattorizzare
%
%   I valori di ritorno sono:
%       A -> la matrice modificata contenente nella parte
%            triangolare superiore U e nella parte strettamente
%            triangolare inferiore L
%
%   See Also SOLVELU, FATTPALU, FATTLDLT, FATTQR
  n=length(A);
  for i=1:n-1
     if A(i,i)==0
        disp('Non è possibile fattorizzare la matrice come LU')
        return
     end
(1)  A(i+1:n,i)=A(i+1:n,i)/A(i,i);
(2)  A(i+1:n,i+1:n)=A(i+1:n,i+1:n)-A(i+1:n,i)*A(i,i+1:n);
  end