%-----------------------------------------------------------------------
%             -------------                           ---------
%             function file                           aprimor.m
%             -------------                           ---------
%-----------------------------------------------------------------------
%                             --------------------------
%                             Metodo de Euler Aprimorado
%                             --------------------------
%
%	Problema de Valor Inicial:   Y'(t)  =  F(t,Y(t))  ,  t > 0
%   --------------------------
%                                     Y(0)  =  Yo .
%
%
%       Variaveis de Entrada
%---------------------------
%	fname : Nome da Function que contem a funcao  F(t,Y(t))
%
%	 Yold : Aproximacao da Solucao no Tempo t(j)
%
%           h : Passo no Tempo  t(j)
%
%           t : Tempo t(j)
%
%
%       Variavel de Saida
%------------------------
%	 Ynew : Aproximacao da Solucao no Tempo t(j+1)  =  t(j) + h
%
        function [Ynew]  =  aprimor(fname,Yold,h,t)
%--------------------------------------------------
%
%	Preditor  -  Metodo de Euler
%--------------- ---------------
%
	FYold  =  feval(fname,t,Yold);
%
	Ystar  =  Yold  +  h*FYold;
%
%	Corretor
%-----------
%
	FYstar  =  feval(fname,(t+h),Ystar);
%
	  Ynew  =  Yold  +  h*( FYold  +  FYstar )*0.5;
%
%-----------------------------------------------------------------------
%-----------------------------------------------------------------------
%	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
%-----------------------------------------------------------------------
