lambert_regula_falsi

lambert_regula_falsi
  Use regula_falsi to estimate solution to a lambert equation.

  F(0) = -5
  F(2) = 9.77811
Step 1: alpha = 0.661662
Step 2: alpha = 0.272833
Step 3: alpha = 1.29581
Step 4: alpha = 0.462974
Step 5: alpha = 0.152919
Step 6: alpha = 0.143686
Step 7: alpha = 0.0266768
Step 8: alpha = 0.00330483

  After 8 iterations:
  F(1.32672) = -3.67953e-06
  F(1.32672) = -3.32852e-10

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

lambert_regula_falsi
  Use regula_falsi to estimate a root of the lambert function.

  F(0) = -5
  F(2) = 9.77811
Step 1: alpha = 0.661662
Step 2: alpha = 0.727167
Step 3: alpha = 0.825238
Step 4: alpha = 0.909529
Step 5: alpha = 0.959302
Step 6: alpha = 0.982968
Step 7: alpha = 0.993098
Step 8: alpha = 0.997241
Step 9: alpha = 0.998903
Step 10: alpha = 0.999565
Step 11: alpha = 0.999827
Step 12: alpha = 0.999932
Step 13: alpha = 0.999973
Step 14: alpha = 0.999989
Step 15: alpha = 0.999996
Step 16: alpha = 0.999998
Step 17: alpha = 0.999999
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 = 3.29798e-16

  After 41 iterations:
  F(1.32672) = -1.77636e-15
  F(1.32672) = 0

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

lambert_regula_falsi
  Use regula_falsi to estimate a root of the lambert function.

  F(0) = -5
  F(2) = 9.77811
Step 1: alpha = 0.661662
Step 2: alpha = 0.727167
Step 3: alpha = 0.825238
Step 4: alpha = 0.909529
Step 5: alpha = 0.959302
Step 6: alpha = 0.982968
Step 7: alpha = 0.993098
Step 8: alpha = 0.997241
Step 9: alpha = 0.998903
Step 10: alpha = 0.999565
Step 11: alpha = 0.999827
Step 12: alpha = 0.999932
Step 13: alpha = 0.999973
Step 14: alpha = 0.999989
Step 15: alpha = 0.999996
Step 16: alpha = 0.999998
Step 17: alpha = 0.999999
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 = 3.29798e-16

  After 41 iterations:
  F(1.32672) = -1.77636e-15
  F(1.32672) = 0

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