%-----------------------------------------------------------------------
%             -----------                           ----------
%             script file                           tstode11.m
%             -----------                           ----------
%-----------------------------------------------------------------------
%                             ---------------
%                             Metodo de Euler
%                             ---------------
%
%	Problema de Valor Inicial:   Y'(t)  =  F(t,Y(t))  ,  t > 0
%   --------------------------
%                                     Y(0)  =  Yo .
%
%
%
%       -----------------------------------------
%	    Programa Teste para o Procedimento  ode11
%------------------------------------------------
%
%
%	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 = 'queda';
        Tinicial = 0.0;
  	      Tfinal = 6.0;
   	          Y0 = 0.0;
               h = 0.01;
%
%
%
%	Chamada do Procedimento ode11
%------------------------------------
	[T,Y]  =  ode11(fname,Tinicial,Tfinal,Y0,h);
%--------------------------------------------------
%
%
	figure(1), plot(T,Y)
	title('Programa Teste para o Metodo de Euler')
        xlabel('F(t,Y(t))  =  1.0  -  0.5*Y*Y   com  Y(0) = 0')
%
%
%   -----------------------------------------
%	Saida do Resultado em um Arquivo de Dados
%--------------------------------------------
%
%    passo = 15;
%    Tm = T(1:passo:600);
%    Ym = Y(1:passo:600);
%       Nome do Arquivo de Dados
%-------------------------------
	name = 'ode11.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  ode11.dat')
%
%	
%-----------------------------------------------------------------------
%-----------------------------------------------------------------------
%-----------------------------------------------
%	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/
%-----------------------------------------------
