03 March 2024 07:52:58 AM

zero_itp_test():
  C version
  Test zero_itp(), which seeks a root of a function F(X)
  in an interval [A,B].

  f_01(x) = sin ( x ) - x / 2

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

        1.000000        1.895494        2.000000
    3.414710e-01   -2.220446e-14   -9.070257e-02
  Number of calls to F = 8

  f_02(x) = 2 * x - exp ( - x )

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

        0.000000        0.351734        1.000000
   -1.000000e+00    6.363410e-11    1.632121e+00
  Number of calls to F = 8

  f_03(x) = x * exp ( - x )

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

       -1.000000        0.000000        0.500000
   -2.718282e+00    2.513422e-12    3.032653e-01
  Number of calls to F = 8

  f_04(x) = exp ( x ) - 1 / ( 100 * x * x )

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

        0.000100        0.095345       20.000000
   -9.999990e+05    3.377484e-07    4.851652e+08
  Number of calls to F = 31

  f_05(x) = (x+3) * (x-1) * (x-1)

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

       -5.000000       -3.000000        2.000000
   -7.200000e+01    0.000000e+00    5.000000e+00
  Number of calls to F = 2

  User has requested additional verbose output.
  step   [a,    b]    x    f(x)

0  [1,2]  f(1.5)=-0.125
1  [1.5,2]  f(1.56515)=0.268999
2  [1.5,1.56515]  f(1.52152)=0.000827384
3  [1.5,1.52152]  f(1.52128)=-0.00056417

  f_06(x) = x^3 - x - 2

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

        1.000000        1.521402        2.000000
   -2.000000e+00    1.315442e-04    4.000000e+00
  Number of calls to F = 4

zero_itp_test():
  Normal end of execution.

03 March 2024 07:52:58 AM