25-Mar-2024 22:45:38

cpr_test():
  MATLAB/Octave version 5.2.0
  Test cpr(), which uses the Chebyshev Proxy Rootfinder
  to determine all real zeros of a smooth function
  over an interval [a,b].

bessel_test():
  Seek real roots of the J0 Bessel function.

  Interstitial interpolation error norm = 8.32667e-15
  Number of roots = 12
  Roots computed by CPR:
    2.4048
    5.5201
    8.6537
   11.7915
   14.9309
   18.0711
   21.2116
   24.3525
   27.4935
   30.6346
   33.7758
   36.9171
  Maximum residual at roots = 1.03988e-14
  Exact roots:
    2.4048
    5.5201
    8.6537
   11.7915
   14.9309
   18.0711
   21.2116
   24.3525
   27.4935
   30.6346
   33.7758
   36.9171
  Exact f(roots):
  -1.0239e-16
  -1.4117e-17
   1.3004e-17
  -7.9654e-16
  -1.0152e-15
  -7.3745e-16
   7.5595e-17
  -7.2262e-16
   7.8726e-16
   5.8966e-16
   1.9964e-16
  -8.1131e-16
  f ( froots ):
  -7.6794e-15
  -2.5411e-16
  -2.0287e-15
   5.1830e-15
  -1.7581e-15
  -1.0399e-14
  -4.8156e-15
   6.1710e-15
   8.8970e-15
  -9.1415e-15
   7.0285e-15
   1.2179e-16

newton_test():
  Seek real roots of the Newton example x^3-2x-5=0.

  Interstitial interpolation error norm = 5.68434e-14
  Number of roots = 1
  Maximum residual at roots = 8.88178e-16
  Exact real root   :    2.094551481542328
  Computed real root:    2.094551481542327
  Exact f(roots)    :  1.4210854715202e-14
  Computed f(roots) : -8.881784197001252e-16
  Error             : 1.332267629550188e-15

jenkins_test():
  Seek real roots of the Jenkins polynomial
  p(x) = x^4 + 5.6562x^3 + 5.8854x^2
             + 7.3646x   + 6.1354

  Interstitial interpolation error norm = 6.82121e-13
  Maximum residual at roots = 6.21725e-14
  Number of roots = 2
  Computed roots:
  -4.674054017161708e+00
  -9.999999999999927e-01
  Exact roots:
  -1.000000000000000
  -4.674054017161709
  f (Computed roots ):
  -3.996802888650564e-14
   6.217248937900877e-14
  f(Exact roots):
  -8.881784197001252e-16
  -1.865174681370263e-14

cpr_test():
  Normal end of execution.
25-Mar-2024 22:45:38