function x = SECANTEsimple(init1,init2,f,Niter)
x1 = init1; x2 = init2;
for i=1:Niter 
   delta = f(x2) - f(x1);  
   if delta == 0 
      x1 = (x1+x2)/2; x2 = (x1+x2)/2; 
   else 
      c  = (x1.*f(x2) -x2.*f(x1))./delta; x1 = x2; x2 = c;    
   end  
end
x = x2;   
endfunction
