%-----------------------------------------------------------------------
%             -----------                           ----------
%             script file                           tstode13.m
%             -----------                           ----------
%-----------------------------------------------------------------------
%                             --------------------------
%                             Metodo de Euler Modificado
%                             --------------------------
%
%	Problema de Valor Inicial:   Y'(t)  =  F(t,Y(t))  ,  t > 0
%       --------------------------
%                                     Y(0)  =  Yo .
%
%
%
%       -----------------------------------------
%	Programa Teste para o Procedimento  ode13
%------------------------------------------------
%
%
%	Variaveis de Entrada
%---------------------------
%	   fname : Nome da Function que contem a funcao  F(t,Y(t))
%
%	Tinicial : Tempo Inicial
%
%         Tfinal : Tempo Final
%
%             Y0 : Condicao Inicial
%
%              h : Passo no Tempo
%
%
%	Variavel de Saida
%------------------------
%	 Y : Vetor contendo a Solucao Aproximada
%
%	 T : Vetor contendo os Pontos de Integracao do Tempo
%
%
%	Dados de Entrada
%-----------------------
%
	   fname = 'fluido';
        Tinicial = 0.0;
	  Tfinal = 8.0;
	      Y0 = 2.0;
               h = 0.01;
%
%
%
%	Chamada do Procedimento ode13
%------------------------------------
	[T,Y]  =  ode13(fname,Tinicial,Tfinal,Y0,h);
%--------------------------------------------------
%
%
	figure(1), plot(T,Y)
	title('Programa Teste para o Metodo de Euler Modificado')
        xlabel('F(t,Y(t))  =  1.0  -  0.5*Y*Y   com  Y(0) = 2.0')
%
%
%       -----------------------------------------
%	Saida do Resultado em um Arquivo de Dados
%------------------------------------------------
%
%       Nome do Arquivo de Dados
%-------------------------------
	name = 'ode13.dat';
%
%       Abre um Arquivo Sequencial para Escrita
%----------------------------------------------
	fid = fopen(name, 'w');
%
%
	fprintf(fid,'          %15.10f        %15.10f\n', [ T' ; Y' ]);
%
%
%	Fecha o Arquivo de Dados
%-------------------------------
	fclose(fid);
%
%
	disp('    ')
	disp('    ')
	disp('         Saida da Solucao no Arquivo  ode13.dat')
%
%
%	
	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/MS211
%
%-----------------------------------------------------------------------
