Vediamo un algoritmo simile alla risoluzione della matrice triangolare inferiore, solo che in questo caso le sostituzioni avverranno all'indietro.
function x=solveUT(A,b) %SOLVEUT Risolve il sistema lineare Ax=b con A matrice triangolare % superiore % % x=SOLVEUT(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 SOLVEUTCS, SOLVELT n=length(b); x=b; for i=n:-1:1 if A(i,i)==0 disp('Matrice non risolvibile') break end for j=i+1:n x(i)=x(i)-A(i,j)*x(j); end x(i)=x(i)/A(i,i); end