Implementazione in Matlab

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