hump_regula_falsi

hump_regula_falsi
  Use regula_falsi to estimate a root of the hump function.

  F(10) = -5.9773
  F(0) = 5.17647
Step 1: alpha = 0.4641
Step 2: alpha = 0.468366
Step 3: alpha = 0.496983
Step 4: alpha = 0.623868
Step 5: alpha = 0.664884
Step 6: alpha = 0.729801
Step 7: alpha = 0.809779
Step 8: alpha = 0.883887
Step 9: alpha = 0.936835
Step 10: alpha = 0.968109
Step 11: alpha = 0.984554
Step 12: alpha = 0.992676
Step 13: alpha = 0.996562
Step 14: alpha = 0.998394
Step 15: alpha = 0.999252
Step 16: alpha = 0.999652
Step 17: alpha = 0.999838
Step 18: alpha = 0.999925
Step 19: alpha = 0.999965
Step 20: alpha = 0.999984
Step 21: alpha = 0.999992
Step 22: alpha = 0.999996
Step 23: alpha = 0.999998
Step 24: alpha = 0.999999
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.29955) = -6.21725e-15
  F(1.0934) = 8.48293

  it <= itmax?  true
  |xp - xn| <= xtol?  false
  |min(f(xn),f(xp))| <= ftol?  true
diary kepler_regula_falsi.txt
hump_regula_falsi.m
{Undefined variable "hump_regula_falsi" or class
"hump_regula_falsi.m".
} 
hump_regula_falsi

hump_regula_falsi
  Use regula_falsi to estimate a root of the hump function.

  F(10) = -5.9773
  F(0) = 5.17647
Step 1: alpha = 0.4641
Step 2: alpha = 0.468366
Step 3: alpha = 0.496983
Step 4: alpha = 0.623868
Step 5: alpha = 0.664884
Step 6: alpha = 0.729801
Step 7: alpha = 0.809779
Step 8: alpha = 0.883887
Step 9: alpha = 0.936835
Step 10: alpha = 0.968109
Step 11: alpha = 0.984554
Step 12: alpha = 0.992676
Step 13: alpha = 0.996562
Step 14: alpha = 0.998394
Step 15: alpha = 0.999252
Step 16: alpha = 0.999652
Step 17: alpha = 0.999838
Step 18: alpha = 0.999925
Step 19: alpha = 0.999965
Step 20: alpha = 0.999984
Step 21: alpha = 0.999992
Step 22: alpha = 0.999996
Step 23: alpha = 0.999998
Step 24: alpha = 0.999999
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.29955) = -6.21725e-15
  F(1.0934) = 8.48293

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