% Resolucion de un sistema triangular superior AX=B por el metodo de
% sustitucion regresiva. El metodo funciona solo si todos los elementos
% diagonales son distintos de cero. Primero se calcula X
% y luego se usa la regla
% X
%
%Sintaxis
% X=backsub(A,B)
%
%Entrada
% - A es una matriz triangular superior invertible de orden n x n (cuadrada)
% - B es una matriz de orden n x 1 (columna)
%Resultados
% - X es la solucin del sistema lineal
%Calculo de la dimencion de B e inicializacion de X
n=length(B);
X=zeros(n,1);
X(n)=B(n)/A(n,n);
for k=n-1:-1:1
X(k)=(B(k)-A(k,k+1:n)*X(k+1:n))/A(k,k);
end
No hay comentarios:
Publicar un comentario