wiggle_regula_falsi

wiggle_regula_falsi
  Use regula_falsi to estimate a root of the wiggle equation.

  F(8) = -3.13768
  F(0) = 3.1609
Step 1: alpha = 0.501843
Step 2: alpha = 0.387454
Step 3: alpha = 0.591673
Step 4: alpha = 0.446437
Step 5: alpha = 0.29934
Step 6: alpha = 0.293545
Step 7: alpha = 0.416768
Step 8: alpha = 0.545387
Step 9: alpha = 0.904422
Step 10: alpha = 0.993211
Step 11: alpha = 0.999597
Step 12: alpha = 0.999976
Step 13: alpha = 0.999999
Step 14: alpha = 1
Step 15: alpha = 1
Step 16: alpha = 1
Step 17: alpha = 1
Step 18: alpha = 1
Step 19: alpha = 6.02679e-14

  After 19 iterations:
  F(1.03868) = -1.72085e-14
  F(1.03868) = 1.08913e-13

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

wiggle_regula_falsi
  Use regula_falsi to estimate a root of the wiggle equation.

  F(8) = -3.13768
  F(0) = 3.1609
Step 1: alpha = 0.501843
Step 2: alpha = 0.387454
Step 3: alpha = 0.591673
Step 4: alpha = 0.446437
Step 5: alpha = 0.29934
Step 6: alpha = 0.293545
Step 7: alpha = 0.416768
Step 8: alpha = 0.545387
Step 9: alpha = 0.904422
Step 10: alpha = 0.993211
Step 11: alpha = 0.999597
Step 12: alpha = 0.999976
Step 13: alpha = 0.999999
Step 14: alpha = 1
Step 15: alpha = 1
Step 16: alpha = 1
Step 17: alpha = 1
Step 18: alpha = 1
Step 19: alpha = 6.02679e-14

  After 19 iterations:
  F(1.03868) = -1.72085e-14
  F(1.03868) = 1.08913e-13

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