function [x1,x2]=FAUSSEPOSITIONsimple(init1,init2,f,Niter)
x1 = init1; x2 = init2;
for i=1:Niter 
   delta = f(x2) - f(x1);  
   if delta == 0 
      x1 = x1; x2 = x2;  
   else 
      y  = x2 - (f(x2) .* (x2-x1))./delta ;  
      fy = f(y).* f(x1);  
      if fy < 0
         x1 = x1; x2 = y; 
     else 
         x1 = y; x2 = x2; 
     end
   end  
end 
endfunction      
