%\begin{verbatim}
%-----------------------------------------------------------------------
%              -------------                    ---------
%              function file                    tridLDLt
%              -------------                    ---------
%-----------------------------------------------------------------------
%
        function [Diag,Sub] = tridLDLt(Diag,Sub,ndim)
%---------------------------------------------------
%
%       -----------------------------
%       Metodo de Decomposicao de LDL'
%------------------------------------
%
%       Matriz Tridiagonal Simetrica
%-----------------------------------
%
%       Diag    : Vetor que Contem a Diagonal Principal
%       Sub     : Vetor que Contem a Diagonal Inferior
%       ndim    : Dimensao do Sistema
%
%
%       Precisao para verificar a singularidade da Matriz
%--------------------------------------------------------
        epsilon = 1.0e-10;
%
        for i = 2:ndim
%
%       Verifica se a Matriz e' Singular
%       --------------------------------
        if( abs(Diag(i-1)) < epsilon  ),  errormsg(Diag(i-1),i),  end
%
         SubOld = Sub(i);
         Sub(i) = SubOld/Diag(i-1);
%
        Diag(i) = Diag(i) - SubOld*Sub(i);
%
       end
%--------------------------------------------------------------------
%--------------------------------------------------------------------
%-----------------------------------------------------
%       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}
