function Y=resout_euler_explicite(N,T,y0,fcn)
% 
%
%
%	resoud_euler_explicite : résolution d'une équation différentielle par la méthode d'Euler explicite.
%
% *********************************************************
%
%    Y=resout_euler_explicite(N,T,y0,fcn) calcule les valeurs aux instants 0,h,2h,...,Nh=T par la méthode
%    d'Euler explicite pour l'équation différentielle y't)=f(t,y(t)) sur [0,T] et y(0)=y0.
%	
%	* variables d'entrées :
%      N : nombre de valeurs calculées ; 
%      T : borne supérieure de l'intervalle ;
%      y0 : valeur de y à t=0 ;
%      fcn : chaine de caractères (de type 'f(t,y)') représentant la fonction f.
%
%	* variables de sortie :
%     Y : valeurs calculée aux instants 0,h,2h,...,Nh=T.
%
% 
% 
% ************ Fonctions auxiliaires utilisées ************
%
%       aucune
%
% *********************************************************
%



% Contrôles d'entrée

% nombre d'arguments
if nargin~=4
   error('nombre d''arguments de la fonction incorrect');
end
   
% Corps d'algorithme
h=T/N;
Y=zeros(1,N+1);
Y(1)=y0;
auxiY=y0;
for i=1:N
   t=h*(i-1);
   y=auxiY;
   auxiY=auxiY+h*eval(fcn);
   Y(i+1)=auxiY;
end




