cubic_regula_falsi

cubic_regula_falsi
  Use regula_falsi to estimate a root of a cubic function.

  F(1) = -1 (negative)
  F(2) = 5 (positive)
Step 1: alpha = 0.833333
Step 2: alpha = 0.896266
Step 3: alpha = 0.946009
Step 4: alpha = 0.974746
Step 5: alpha = 0.988809
Step 6: alpha = 0.995163
Step 7: alpha = 0.997932
Step 8: alpha = 0.99912
Step 9: alpha = 0.999626
Step 10: alpha = 0.999841
Step 11: alpha = 0.999933
Step 12: alpha = 0.999971
Step 13: alpha = 0.999988
Step 14: alpha = 0.999995
Step 15: alpha = 0.999998
Step 16: alpha = 0.999999
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(   1.324717957244746) = -8.88178e-16 (negative)
  F(                   2) = 5 (positive)

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

cubic_regula_falsi
  Use regula_falsi to estimate a root of a cubic function.

  F(1) = -1 (negative)
  F(2) = 5 (positive)
Step 1: alpha = 0.833333
Step 2: alpha = 0.896266
Step 3: alpha = 0.946009
Step 4: alpha = 0.974746
Step 5: alpha = 0.988809
Step 6: alpha = 0.995163
Step 7: alpha = 0.997932
Step 8: alpha = 0.99912
Step 9: alpha = 0.999626
Step 10: alpha = 0.999841
Step 11: alpha = 0.999933
Step 12: alpha = 0.999971
Step 13: alpha = 0.999988
Step 14: alpha = 0.999995
Step 15: alpha = 0.999998
Step 16: alpha = 0.999999
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(   1.324717957244746) = -8.88178e-16 (negative)
  F(                   2) = 5 (positive)

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