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