function iter(a, b, eps)
  xx=x=(a+b)/2
  iter_count = 0
  while true
    xx=x
    x=fi(x)
    iter_count+=1
    if (abs((xx-x)/xx)<eps) break end
    print(x," ",xx)
  end
  return (x,iter_count)
end