Implementazione in Matlab

function x=solveLTcs(A,b)
%SOLVEUTCS Risolve il sistema lineare Ax=b con A matrice triangolare
%   superiore utilizzando il metodo column sweep
%
%   x=SOLVEUTCS(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 SOLVEUT
  n=length(b);
  x=b;
  for j=n:-1:1
     if A(j,j)==0
        disp('Matrice non risolvibile')
        break
     end
     x(j)=x(j)/A(j,j);
     for i=1:j-1
        x(i)=x(i)-A(i,j)*x(j);
     end
  end