%\begin{verbatim}
%-----------------------------------------------------------------------
%              -------------                   ---------
%              function file                   TridBloco
%              -------------                   ---------
%-----------------------------------------------------------------------
%              Gera um Sistema Tridiagonal por Blocos
%----------------------------------------------------
%
        function [Diag,Sup,Sub,b] = TridBloco(ndim,mblock)
%---------------------------------------------------------
%
        [D] = clock;      rand('seed',D(6))   
%        
%       Diagonal Principal
%-------------------------
        for k = 1:mblock
        Diag(:,:,k) = rand(ndim,ndim);
        end        
%
%       Diagonal Superior
%------------------------
        for k = 1:(mblock-1)
        Sup(:,:,k) = rand(ndim,ndim);
        end                             
%
%       Diagonal Inferior
%------------------------
        for k = 1:(mblock-1)
        Sub(:,:,k) = rand(ndim,ndim);
        end                      
%
%       Vetor do Lado Direito
%----------------------------
        for k = 1:mblock
        b(:,k) = rand(ndim,1);
        end                      
%
%-----------------------------------------------------------------------
%-----------------------------------------------------------------------
%--------------------------------------------------
%       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/
%--------------------------------------------------
%\end{verbatim}
