function x=secante(init1,init2,N);

% x=secante(init1,init2,N); 
% 
% Approche via la méthode de la sécante (initiée à init1, init2) le zéro
% réel de la fonction polynomiale x -->  x^5+ x^3/2 +1 entre -1 et -1/2.  

x1=init1;
x2=init2;
for i=1:N 
   y=x1;
   x1=x2;
   delta = x2^5 + x2^3/2 - y^5 -y^3/2; 
   if delta == 0
      x=(x1+x2)./2;
   else 
      x2=x2 - (x2^5 + x2^3/2 + 1)*(x2-y)./delta;
   end
end
x=x2;
