Resolução retroativa

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

function x=resret (A,b)
 [l,c]=size(A); // dimensão do argumento/tamanha da matriz A.
 for i=l:-1:1
 soma=0;
 for j=i+1:c
 soma=soma+A(i,j)*x(j);
 end
 x(i)=(b(i)-soma)/A(i,i)
 end
endfunction
A=[2 3 -1;0 -2 -1;0 0 5];
b=[5;-7;15];
x=resret(A,b)