kepler_regula_falsi

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

  F(0) = 5
  F(10) = -6.08804
Step 1: alpha = 0.549064
Step 2: alpha = 0.317827
Step 3: alpha = 0.669367
Step 4: alpha = 0.106471
Step 5: alpha = 0.143208
Step 6: alpha = 0.0355241
Step 7: alpha = 0.00403346
Step 8: alpha = 0.00014688

  After 8 iterations:
  F(3.79083) = -9.71978e-10
  F(3.79083) = 1.11022e-15

  it <= itmax?  true
  |a-b| <= xtol?  true
  |f(b)| <= ftol?  true
diary lambert_regula_falsi.txt
kepler_regula_falsi

kepler_regula_falsi
  Use regula_falsi to estimate solution to Kepler's equation.

  F(10) = -6.08804
  F(0) = 5
Step 1: alpha = 0.450936
Step 2: alpha = 0.226997
Step 3: alpha = 0.363754
Step 4: alpha = 0.830282
Step 5: alpha = 0.987612
Step 6: alpha = 0.999275
Step 7: alpha = 0.999958
Step 8: alpha = 0.999998
Step 9: alpha = 1
Step 10: alpha = 1
Step 11: alpha = 1
Step 12: alpha = 1
Step 13: alpha = 1
Step 14: alpha = 1
Step 15: alpha = 1
Step 16: alpha = 1
Step 17: alpha = 1
Step 18: alpha = 1
Step 19: alpha = 1
Step 20: alpha = 1
Step 21: alpha = 1
Step 22: alpha = 1
Step 23: alpha = 1
Step 24: alpha = 1
Step 25: alpha = 1
Step 26: alpha = 1
Step 27: alpha = 1
Step 28: alpha = 1
Step 29: alpha = 1
Step 30: alpha = 1
Step 31: alpha = 1
Step 32: alpha = 1
Step 33: alpha = 1
Step 34: alpha = 1
Step 35: alpha = 1
Step 36: alpha = 1
Step 37: alpha = 1
Step 38: alpha = 1
Step 39: alpha = 1
Step 40: alpha = 1
Step 41: alpha = 1
Step 42: alpha = 1
Step 43: alpha = 1
Step 44: alpha = 1
Step 45: alpha = 1
Step 46: alpha = 1
Step 47: alpha = 1
Step 48: alpha = 1
Step 49: alpha = 1
Step 50: alpha = 1

  After 50 iterations:
  F(3.79083) = -2.22045e-16
  F(3.48575) = 0.839447

  it <= itmax?  true
  |xp - xn| <= xtol?  false
  |min(f(xn),f(xp))| <= ftol?  true
diary lambert_regula_falsi.txt
kepler_regula_falsi

kepler_regula_falsi
  Use regula_falsi to estimate solution to Kepler's equation.

  F(10) = -6.08804
  F(0) = 5
Step 1: alpha = 0.450936
Step 2: alpha = 0.226997
Step 3: alpha = 0.363754
Step 4: alpha = 0.830282
Step 5: alpha = 0.987612
Step 6: alpha = 0.999275
Step 7: alpha = 0.999958
Step 8: alpha = 0.999998
Step 9: alpha = 1
Step 10: alpha = 1
Step 11: alpha = 1
Step 12: alpha = 1
Step 13: alpha = 1
Step 14: alpha = 1
Step 15: alpha = 1
Step 16: alpha = 1
Step 17: alpha = 1
Step 18: alpha = 1
Step 19: alpha = 1
Step 20: alpha = 1
Step 21: alpha = 1
Step 22: alpha = 1
Step 23: alpha = 1
Step 24: alpha = 1
Step 25: alpha = 1
Step 26: alpha = 1
Step 27: alpha = 1
Step 28: alpha = 1
Step 29: alpha = 1
Step 30: alpha = 1
Step 31: alpha = 1
Step 32: alpha = 1
Step 33: alpha = 1
Step 34: alpha = 1
Step 35: alpha = 1
Step 36: alpha = 1
Step 37: alpha = 1
Step 38: alpha = 1
Step 39: alpha = 1
Step 40: alpha = 1
Step 41: alpha = 1
Step 42: alpha = 1
Step 43: alpha = 1
Step 44: alpha = 1
Step 45: alpha = 1
Step 46: alpha = 1
Step 47: alpha = 1
Step 48: alpha = 1
Step 49: alpha = 1
Step 50: alpha = 1

  After 50 iterations:
  F(3.79083) = -2.22045e-16
  F(3.48575) = 0.839447

  it <= itmax?  true
  |xp - xn| <= xtol?  false
  |min(f(xn),f(xp))| <= ftol?  true
diary lambert_regula_falsi.txt