%\begin{verbatim}
%-----------------------------------------------------------------------
%             -----------                             --------
%             script file                             testLDLt
%             -----------                             --------
%-----------------------------------------------------------------------
%       --------------
%       Programa teste
%---------------------
%
%       ------------------------------------     -----------------
%       Sistema Linear Tridiagonal Simetrico  -  Decomposicao  LDL'
%-------------------------------------------     -----------------
%
%       Variaveis
%----------------
%       x       :  Solução do Sistema Linear
%
%       diag    : Diagonal Principal da Matriz do Sistema Linear
%       sub     : Diagonal Inferior   "    "    "    "       "
%       b       : Lado Direito do Sistema
%       ndim    : Dimensao do Sistema
%
%
%       Diag    : Diagonal Principal da  Decomposicao LDL'
%       Sub     : Diagonal Inferior   "       "        "  
%
%
%       ----------------------------------------
%       Construcao do Sistema Linear Tridiagonal
%-----------------------------------------------
%
%       dimensao do sistema linear
%       -------------------------
        ndim = 20;         
%
%
%       Diagonal Principal
%-------------------------
        for i = 1:ndim
        diag(i) = 4.0;
        end
%
%
%       Diagonal Inferior
%------------------------
        for i = 2:ndim
        sub(i) = -1.0;
        end
%
%       Lado Direito do Sistema Linear
%------------------------------------
         b(1) = 3.0;
        for i = 2:(ndim-1)
        b(i) = 2.0;
        end 
        b(ndim) = 3.0;  
%
%       -------------------------------------------------
%       Decomposicao LDL' da Matriz Tridiagonal Simetrica
%--------------------------------------------------------
        [Diag,Sub] = tridLDLt(diag,sub,ndim);
%
%       OBSERVACAO: A Decomposicao LDL'  é armazenada nos
%                   vetores,  Diag, Sup e Sub.
%
%       ---------------------------------------------------------
%       Solucao do Sistema  - Resolucao dos Sistemas Triangulares
%----------------------------------------------------------------
%
        [x] = solveLDLt(Diag,Sub,b,ndim);
%
%       Mostra a solucao do Sistema Linear na Tela
%-------------------------------------------------
        disp('   ')
        disp('   Solucao do Sistema Linear Tridiagonal Simetrico')
        disp('   -----------------------------------------------')
        disp('   ')
        format long, 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/MS211
%--------------------------------------------------
%-----------------------------------------------------------------------
%
%
%
%       -----------------------------------
%       Armazenamento da Matriz Tridiagonal
%------------------------------------------
%
%
%
%       Diagonal Principal
%-------------------------
%
%       diag(i) = a(i,i)  para  i = 1,2, ... , ndim
%
%
%       Diagonal Inferior
%------------------------
%
%       sub(i) = a(i,i-1)  para  i = 2,3, ... , ndim
%
%-----------------------------------------------------------------------
%-----------------------------------------------------------------------
%-----------------------------------------------------
%       Petronio Pulino
%       Departamento de Matematica Aplicada - D.M.A
%       IMECC - UNICAMP
%       e-mail: pulino@ime.unicamp.br
%         http://www.ime.unicamp.br/~pulino/
%-----------------------------------------------------
%\end{verbatim}
