19 March 2024 03:56:33 PM

zero_chandrupatla_test():
  C++ version
  zero_chandrupatla() seeks a root of a function f(x)
  in an interval [a,b].

  f_01(x) = x^3 - 2 x - 5

      A                 Z             B
    F(A)              F(Z)          F(B)

               2         2.09455               3
              -1    -3.49042e-06              16
  Number of calls to F = 7

  f_02(x) = 1 - 1/x^2

      A                 Z             B
    F(A)              F(Z)          F(B)

             0.5               1            1.51
              -3    -9.97854e-07        0.561423
  Number of calls to F = 8

  f_03(x) = ( x - 3 )^3

      A                 Z             B
    F(A)              F(Z)          F(B)

               0               3               5
             -27     6.93889e-18               8
  Number of calls to F = 21

  f_04(x) = 6 * ( x - 2 )^5

      A                 Z             B
    F(A)              F(Z)          F(B)

               0               2               5
            -192    -1.51461e-28            1458
  Number of calls to F = 21

  f_05(x) = x^9

      A                 Z             B
    F(A)              F(Z)          F(B)

              -1      3.8147e-06               4
              -1     1.71057e-49          262144
  Number of calls to F = 21

  f_06(x) = x^19

      A                 Z             B
    F(A)              F(Z)          F(B)

              -1      3.8147e-06               4
              -1     1.1162e-103     2.74878e+11
  Number of calls to F = 21

  f_07(x) = x e^(-1/x2)

      A                 Z             B
    F(A)              F(Z)          F(B)

              -1     -0.00827833               4
       -0.367879              -0         3.75765
  Number of calls to F = 8

  f_08(x) = -(3062(1-xi)e^(-x)/(xi+(1-xi)e^(-x)) - 1013 + 1628/x

      A                 Z             B
    F(A)              F(Z)          F(B)

          0.0002         1.03754               2
     8.13781e+06     0.000335259         -438.26
  Number of calls to F = 9

  f_09(x) = e^x - 2 - 0.01/x^2 + 0.000002/x^3

      A                 Z             B
    F(A)              F(Z)          F(B)

          0.0002        0.703205               1
         -0.9998    -4.38667e-08        0.708284
  Number of calls to F = 7

zero_chandrupatla_test():
  Normal end of execution.

19 March 2024 03:56:33 PM