triangular.m

来自「programas matlab metodos numericos」· M 代码 · 共 15 行

function [x]=triangular(A)
n=size(A); %tama駉 de la matrix
c=min(n);
for i=1:c-1 % busca la matriz U y el vector L*v (la ultima columna) 
	[a,b]=max(abs(A(i:n(1),i)));
   if a~=0
   % da la posicion del maximo pivote (b) y (a) es su valor
	A([i,b+i-1],i:n(2))=A([i+b-1,i],i:n(2)); % intercambia la fila del maximo con la i-esima
   filas=i+1:n(1);
  	A(i,filas:n(2))=A(i,filas:n(2))/a; %normaliza la fila del pivote
	A(filas,filas:n(2))=A(filas,filas:n(2))-A(filas,i)*A(i,filas:n(2)); % triangula
   end
end
x=triu(A,1)+eye(n);
% resolucion del sistema