function [x,y]=euler_explicite2(f,ya,a,b,n)
  h=(b-a)/n;
  %initialisation du vecteur des abscisses
  x=[];
  t=a;
  x=[x t];
  %initialisation du vecteur des ordonnees
  y=[ya];
  for i = 1:n
    y=[ y  y(i)+h*f(t,y(i)) ];
    t=t+h;
    x=[x t];
  endfor
endfunction
