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