07-Jan-2022 17:56:15

CAUCHY_PRINCIPAL_VALUE_TEST():
  MATLAB/Octave version 9.8.0.1380330 (R2020a) Update 2.
  Test CAUCHY_PRINCIPAL_VALUE.

CAUCHY_PRINCIPAL_VALUE_TEST01:
  CPV of Integral ( -1 <= t <= 1 ) exp(t) / t dt

   N           Estimate             Error

   2         2.112977728449278      0.00152402
   4         2.114501718105385     3.26446e-08
   6         2.114501750751337     1.33671e-12
   8         2.114501750751458      1.4575e-12

CAUCHY_PRINCIPAL_VALUE_TEST02:
  Compute CPV ( Integral ( 1-delta <= t <= 1+delta ) 1/(1-t)^3 dt )
  Try this for delta = 1, 1/2, 1/4.

   N          Estimate                  Exact                  Error

   2       -0.7397260273972603       -0.7363873204868445      0.00333871
   4       -0.7364116807468272       -0.7363873204868445     2.43603e-05
   6       -0.7363867923558026       -0.7363873204868445     5.28131e-07
   8       -0.7363873252791506       -0.7363873204868445     4.79231e-09

   2       -0.3425852498017447       -0.3425632583544805     2.19914e-05
   4       -0.3425633493215302       -0.3425632583544805      9.0967e-08
   6        -0.342563258302464       -0.3425632583544805     5.20165e-11
   8        -0.342563258354488       -0.3425632583544805     7.54952e-15

   2       -0.1678240178701501       -0.1678238552950642     1.62575e-07
   4       -0.1678238554770549       -0.1678238552950642     1.81991e-10
   6       -0.1678238552950588       -0.1678238552950642     5.38458e-15
   8       -0.1678238552950641       -0.1678238552950642     1.38778e-16

CAUCHY_PRINCIPAL_VALUE_TEST
  Normal end of execution.

07-Jan-2022 17:56:15