function [x,y]=euler_implicite(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
    g=@(z)[y(i-1)+h*f(x(i),z)-z];
    y(i)=fsolve(g,y(i-1));
  endfor
endfunction