09 May 2025   8:57:26.042 PM
 
circle_integrals_test():
  Fortran90 version
  Test circle_integrals().
 
TEST01
  Use CIRCLE01_SAMPLE to compare exact and
  estimated integrals along the circumference 
  of the unit circle in 2D.

  Number of sample points used is   1024

  If any exponent is odd, the integral is zero.
  We restrict this test to randomly chosen even exponents.

  Ex  Ey     MC-Estimate           Exact      Error

   4   4    0.141798        0.147262        0.55E-02
   4   0     2.33754         2.35619        0.19E-01
   6  10    0.832772E-02    0.862864E-02    0.30E-03
   6   8    0.147493E-01    0.153398E-01    0.59E-03
   6   8    0.147493E-01    0.153398E-01    0.59E-03
   8   8    0.642156E-02    0.671117E-02    0.29E-03
   0   8     1.79076         1.71806        0.73E-01
   6   2    0.239049        0.245437        0.64E-02
   4   8    0.416459E-01    0.429515E-01    0.13E-02
   4   4    0.141798        0.147262        0.55E-02
  10  10    0.144526E-02    0.151001E-02    0.65E-04
   6   8    0.147493E-01    0.153398E-01    0.59E-03
   4   0     2.33754         2.35619        0.19E-01
  10  10    0.144526E-02    0.151001E-02    0.65E-04
   4  10    0.268966E-01    0.276117E-01    0.72E-03
  10   0     1.54934         1.54625        0.31E-02
   6  10    0.832772E-02    0.862864E-02    0.30E-03
   2   4    0.383749        0.392699        0.90E-02
  10  10    0.144526E-02    0.151001E-02    0.65E-04
   4   2    0.380848        0.392699        0.12E-01
 
circle_integrals_test():
  Normal end of execution.
 
09 May 2025   8:57:26.042 PM