%-----------------------------------------------------------------------
%             -----------                             ---------
%             script file                             teste01.m
%             -----------                             ---------
%-----------------------------------------------------------------------
%
%
%       --------------
%	Programa teste
%---------------------
%
%
%   --------------------------     ---------------------
%	Sistema Linear Tridiagonal  -  Decomposicao de Crout
%-----------------------------     ---------------------
%
%
%
%	Variaveis
%----------------
%	x	:  Solução do Sistema Linear
%
%
%	diag	: Diagonal Principal da Matriz do Sistema Linear
%	sup	: Diagonal Superior  "    "    "    "       "
%       sub     : Diagonal Inferior  "    "    "    "       "
%	b	: Lado Direito do Sistema
%	ndim	: Dimensao do Sistema
%
%
%	Diag	: Diagonal Principal da  Decomposicao da Crout
%	Sup	: Diagonal Superior  "       "        "    "
%       Sub     : Diagonal Inferior  "       "        "    "
%
%
%       ----------------------------------------
%	Construcao do Sistema Linear Tridiagonal
%-----------------------------------------------
%
	ndim = 10;         %  dimensao do sistema linear
%                             -------------------------
%
%
%	Diagonal Principal
%-------------------------
	for i = 1:ndim
	diag(i) = 4.0;
	end
%
%
%	Diagonal Superior
%-----------------------
	for i = 1:(ndim-1)
	sup(i) = -1.0;
	end
%
%
%	Diagonal Inferior
%-----------------------
	for i = 1:(ndim-1)
	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 de Crout da Matriz Tridiagonal A
%----------------------------------------------------
	[Diag,Sup,Sub] = tridiag(diag,sup,sub,ndim);
%
%
%	OBSERVACAO: A Decomposicao de Crout é armazenada nos
%                   vetores,  Diag, Sup e Sub.
%
%
%
%
%       ---------------------------------------------------------
%	Solucao do Sistema  - Resolucao dos Sistemas Triangulares
%----------------------------------------------------------------
%
	[x] = solvetri(Diag,Sup,Sub,b,ndim);
%
%
%	Mostra a solucao do Sistema Linear na Tela
%-------------------------------------------------
	disp('   ')
	disp('   Solucao do Sistema Linear Tridiagonal')
	disp('   -------------------------------------')
	disp('           ---------------------')
	disp('           Decomposicao de Crout')
	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 Superior
%------------------------
%
%	sup(i) = a(i,i+1)  para  i = 1,2, ... (ndim-1)
%
%
%
%	Diagonal Inferior
%------------------------
%
%	sub(i-1) = a(i,i-1)  para  i = 2,3, ... ndim
%
%
%-----------------------------------------------------------------------
%-----------------------------------------------------------------------
