09 May 2025   9:16:33.856 PM
 
hypercube_integrals_test():
  FORTRAN90 version
  Test hypercube_integrals().

TEST01
  Estimate monomial integrals using Monte Carlo
  over the interior of the unit hypercube in 3 dimensions.

  Number of sample points used is   4192

  We randomly choose the exponents.

  Ex  Ey  Ez     MC-Estimate      Exact           Error

   2   2   4    0.236710E-01    0.222222E-01    0.14E-02
   3   0   4    0.509324E-01    0.500000E-01    0.93E-03
   4   4   4    0.879102E-02    0.800000E-02    0.79E-03
   4   1   4    0.215206E-01    0.200000E-01    0.15E-02
   0   1   1    0.256391        0.250000        0.64E-02
   4   2   1    0.357342E-01    0.333333E-01    0.24E-02
   0   1   1    0.256391        0.250000        0.64E-02
   0   1   2    0.170540        0.166667        0.39E-02
   3   1   4    0.266270E-01    0.250000E-01    0.16E-02
   0   2   0    0.349367        0.333333        0.16E-01
   0   2   1    0.174016        0.166667        0.73E-02
   2   4   2    0.238475E-01    0.222222E-01    0.16E-02
   2   1   3    0.434707E-01    0.416667E-01    0.18E-02
   3   3   0    0.683415E-01    0.625000E-01    0.58E-02
   4   1   1    0.525364E-01    0.500000E-01    0.25E-02
   0   0   0     1.00000         1.00000        0.00E+00
   1   4   0    0.108274        0.100000        0.83E-02
   2   4   1    0.360418E-01    0.333333E-01    0.27E-02
   4   1   0    0.105096        0.100000        0.51E-02
   2   2   2    0.390717E-01    0.370370E-01    0.20E-02

TEST02
  Estimate monomial integrals using Monte Carlo
  over the interior of the unit hypercube in 6 dimensions.

  Number of sample points used is   4192

  We randomly choose the exponents.

  E1  E2  E3  E4  E5  E6     MC-Estimate      Exact           Error

   3   2   1   1   1   2    0.332563E-02    0.347222E-02    0.15E-03
   3   1   2   2   3   1    0.166509E-02    0.173611E-02    0.71E-04
   3   2   2   0   3   3    0.152871E-02    0.173611E-02    0.21E-03
   3   2   1   0   2   3    0.322913E-02    0.347222E-02    0.24E-03
   1   2   0   0   2   2    0.185337E-01    0.185185E-01    0.15E-04
   1   3   3   0   1   1    0.719987E-02    0.781250E-02    0.61E-03
   2   2   0   3   2   1    0.450177E-02    0.462963E-02    0.13E-03
   3   1   0   0   1   0    0.610756E-01    0.625000E-01    0.14E-02
   1   1   2   3   1   2    0.335555E-02    0.347222E-02    0.12E-03
   1   1   0   2   2   3    0.697709E-02    0.694444E-02    0.33E-04
   1   3   2   0   1   1    0.973418E-02    0.104167E-01    0.68E-03
   0   0   0   0   3   3    0.650871E-01    0.625000E-01    0.26E-02
   0   3   0   3   1   0    0.293357E-01    0.312500E-01    0.19E-02
   2   2   3   1   2   3    0.106218E-02    0.115741E-02    0.95E-04
   3   1   3   3   3   0    0.181580E-02    0.195312E-02    0.14E-03
   0   1   1   3   3   2    0.519527E-02    0.520833E-02    0.13E-04
   0   1   0   3   1   3    0.153086E-01    0.156250E-01    0.32E-03
   0   3   0   3   3   3    0.378686E-02    0.390625E-02    0.12E-03
   3   1   3   0   0   0    0.304355E-01    0.312500E-01    0.81E-03
   0   1   3   3   2   2    0.346021E-02    0.347222E-02    0.12E-04
 
HYPERCUBE_INTEGRALS_TEST():
  Normal end of execution.
 
09 May 2025   9:16:33.860 PM