kepler_bisection

kepler_bisection
  Use bisection2 to estimate solution to Kepler's equation.

  F(10) = -6.08804 (negative)
  F(0) = 5 (positive)

  After 24 iterations:
  F(3.79084) = -1.52126e-06 (negative)
  F(3.79083) = 2.43377e-08 (positive)

  it <= itmax?  true
  |xn-xp| <= xtol?  true
  min(|f(xn)|,|f(xp)|) <= ftol?  true
diary off