function x=dichotomus(a,b,delta,e)
  c=a+(b-a)/2-delta
  d=a+(b-a)/2+delta
  while b-a>2*e
  f(c)
  f(d)
  if f(c)>f(d)
    a=c
    else b=d
  endif
  c=a+(b-a)/2-delta
  d=a+(b-a)/2+delta
endwhile
x=(a+b)/2
endfunction