%-----------------------------------------------------------------------
%             -----------                             -----------
%             script file                             tstLUPivo.m
%             -----------                             -----------
%-----------------------------------------------------------------------
%               ----------------------------------------
%               Metodos Numericos para Sistemas Lineares
%              ------------------------------------------
%              Decomposicao  LU  com Pivoteamento Parcial
%                     
%
%
%
%       --------------       --------------------
%       Programa Teste   -   Procedimento  LUPivo
%---------------------       --------------------
%
%
%       --------------------
%       Variaveis de Entrada
%---------------------------
%       A   :  Matriz do Sistema Linear
%
%       b   :  Vetor do Lado Direito
%
%
%       ------------------
%       Variaveis de Saida
%-------------------------
%
%       LU  :  Matriz que contem a Decomposicao  LU
%
%        P  :  Matriz de Permutacao
% 
%        x  :  Solucao do Sistema Linear
%
%
%       ----------------
%       Dados de Entrada
%-----------------------
%
%       Gera uma Matriz Aleatoria
%--------------------------------
        A = rand(10,10);
%
%       Cria o Vetor do Lado Direito do Sistema
%----------------------------------------------
        b = sum(A');    b = b';     %  Define b como vetor coluna
%
%
%       -------------------------------
%       Chamada do Procedimento  LUPivo
%--------------------------------------
        [LU,P] = LUPivo(A);
%
%
%       Resolucao dos Sistemas Triangulares
%------------------------------------------
        [x] = LUSolve(LU,b,P);
%
%
%
        disp('      ')
        disp('      ')
        disp('      ')
        disp('     Sistemas Lineares - Decomposicao  LU')
        disp('     ------------------------------------')
        disp('              --------------------')
        disp('              Procedimento  LUPivo')
        disp('              --------------------')
        disp('      ')
%
        format long
        disp('    Solucao do Sistema Linear')
        disp('      ')
        disp(x)      
        format short
%
%
%-----------------------------------------------------------------------
%-----------------------------------------------------------------------
%--------------------------------------------------
%       Petronio Pulino
%       Departamento de Matemática Aplicada - D.M.A
%       IMECC - UNICAMP
%       e-mail: pulino@ime.unicamp.br
%       http://www.ime.unicamp.br/~pulino/MS411/
%--------------------------------------------------
%-----------------------------------------------------------------------

