perfect.m

programas matlab metodos numericos

function [x]=perfect(A)
%resuleve un solo sistema con el metodo de doble for
n=size(A); %tama駉 de la matrix
for i=1:n(1)-1 % busca la matriz U y el vector L*v (la ultima columna) 
    [a,b]=max(abs(A(i:n(1),i)));
 % 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
    for j=i+1:n(1)
        A(i,i+1:n(2))=A(i,i+1:n(2))/A(i,i); %normaliza la fila del pivote, a no sirve no tiene el signo correspondiente
        A(j,i+1:n(2))=A(j,i+1:n(2))-A(i,i+1:n(2))*A(j,i); % triangula
    end
end
% resolucion del sistema
x(n(1),1)=A(n(1),n(2))/A(n(1),n(1));%el ultimo es el unico no normalizado
for i=n(1)-1:-1:1
    x(i,1)= A(i,n(2))-A(i,i+1:n(1))*x(i+1:n(1),1); % no divide porque esta normalizado
end
%imprime solucion