%\begin{verbatim}
%-----------------------------------------------------------------------
%             -------------                       ---------
%             function file                       solveLDLt
%             -------------                       ---------
%-----------------------------------------------------------------------
%
        function  [x] = solveLDLt(Diag,Sub,b,ndim)
%------------------------------------------------
%
%       ----------------------------------------------
%       Resolucao de um Sistema Tridiagonal  Simetrico
%-----------------------------------------------------
%
%                --------------------------
%                Metodo de Decomposicao LDL'
%------------------------------------------
%
%
%       Resolucao do Sistema Triangular Inferior
%-----------------------------------------------
        x(1,1) = b(1);
%       
        for i = 2:ndim
        x(i,1) = b(i) - x(i-1)*Sub(i);
        end
%
%       Resolucao do Sistema Diagonal
%------------------------------------
        for i = 1:ndim
        x(i) = x(i)/Diag(i);
        end
%
%       Resolucao do Sistema Triangular Superior
%-----------------------------------------------
        for i = (ndim-1):-1:1
        x(i)  = x(i)  -  x(i+1)*Sub(i+1);
        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}
