%-----------------------------------------------------------------------
%             -----------                             --------
%             script file                             testlu.m
%             -----------                             --------
%-----------------------------------------------------------------------
%              ----------------------------------------
%              Metodos Numericos para Sistemas Lineares
%              ----------------------------------------
%
%
%
%       --------------       ---------------------------
%       Programa Teste   -   Procedimento  lu  do Matlab
%---------------------       ---------------------------
%
%
%       --------------------
%       Variaveis de Entrada
%---------------------------
%       A   :  Matriz do Sistema Linear
%
%       b   :  Vetor do Lado Direito
%
%
%       ------------------
%       Variaveis de Saida
%-------------------------
%
%                  U   :  Matriz Triangular Superior
%
%                  L   :  Matriz Triangular Inferior
%
%                  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  lu  do Matlab
%---------------------------------------------
        [L,U,P] = lu(A);
%
%
%       Calculo do vetor d = P.b
%-------------------------------
        d = P*b;
%
%       Resolucao do Sistema Triangular Inferior
%-----------------------------------------------
        [y] = LTriSol(L,d);
%
%       Resolucao do Sistema Triangular Superior
%----------------------------------------------
        [x] = UTriSol(U,y);
%
%
%
        disp('      ')
        disp('      ')
        disp('      ')
        disp('     Sistemas Lineares - Decomposicao  LU')
        disp('     ------------------------------------')
        disp('          ---------------------------')
        disp('          Procedimento  lu  do Matlab')
        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/
%--------------------------------------------------
%-----------------------------------------------------------------------

