Implementazione in Matlab

Di seguito presentiamo la function Matlab che risolve un sistema lineare con la matrice dei coefficienti triangolare inferiore

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