sretro

Esta é uma pré-visualização de arquivo. Entre para ver o arquivo original

function [x,r]=sretro(Tsup)
% Rotina para calcular a solução de sistemas triagulares superiores
% Entrada: Tsup = matriz triangular superior aumentada
% Saída: x = vetor solução
% r = vetor resíduo
[nL,nc]=size(Tsup);
n=nL;
 
A = Tsup(:,1:n);
b = Tsup(:,nc);
 
x=zeros(n,1);
x(n)=b(n)/A(n,n);
 
for i=n-1:-1:1
 soma=0;
 for j=i+1:n
 soma=soma+A(i,j)*x(j);
 end
 x(i)=(b(i)-soma)/A(i,i);
end
 
r=b-A*x;