function [x,y]=euler_implicitef1(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)/(1-h);
  endfor
endfunction