function x = DICHOTOMIEsimple(init1,init2,f,Niter)
x = init1; y = init2;
  for i = 1:Niter 
    c  = (x + y)/2; u  = f(x)*f(c); 
    if u <= 0
      x = x; y = c;
    else 
      x = c; y = y; 
    end
  end 
x = (x+y)/2; 
end

