8 September 2021 7:13:19.019 PM DISK_MONTE_CARLO_TEST FORTRAN90 version Test the DISK_MONTE_CARLO library. DISK_AREA_TEST DISK_AREA computes the area of a disk with center = (CX,CY) and radius R. ( CX CY ) R Area (-1.377389, 1.427678) 0.894170 2.511827 (-2.692951, 3.238712) 0.390937 0.480136 (-1.980868, 4.219477) 0.879175 2.428292 ( 4.966195, 3.392342) 0.514408 0.831314 ( 1.860457, 1.359361) 0.516244 0.837259 (-0.107215, 1.609617) 0.460673 0.666709 (-0.879364, -0.492748) 0.749650 1.765495 ( 2.869551, -4.339380) 0.256584 0.206829 (-4.990871, 4.912071) 0.042256 0.005610 (-1.896553, -2.349584) 0.331308 0.344837 TEST01 Use DISK01_SAMPLE to estimate integrals in the unit disk. N 1 X^2 Y^2 X^4 X^2Y^2 Y^4 X^6 1 3.14159 0.366748 0.770454 0.428141E-01 0.899425E-01 0.188948 0.499810E-02 2 3.14159 0.389085 0.943470E-01 0.615752E-01 0.125851E-01 0.289394E-02 0.109421E-01 4 3.14159 1.13510 0.207599 0.513835 0.264513E-01 0.474681E-01 0.261947 8 3.14159 0.760603 0.613699 0.562568 0.323150E-01 0.339582 0.481865 16 3.14159 0.861721 0.613562 0.399109 0.910327E-01 0.334486 0.216274 32 3.14159 0.971024 0.653362 0.575740 0.113986 0.321525 0.418641 64 3.14159 0.824500 0.638592 0.399413 0.132138 0.281971 0.243908 128 3.14159 0.795615 0.934208 0.411883 0.157738 0.520588 0.268678 256 3.14159 0.763763 0.784414 0.374409 0.125176 0.384909 0.234469 512 3.14159 0.769155 0.815668 0.368415 0.141798 0.403206 0.220118 1024 3.14159 0.760154 0.792253 0.368916 0.128648 0.408939 0.225226 2048 3.14159 0.795928 0.790131 0.401929 0.132913 0.392164 0.254530 4096 3.14159 0.805292 0.787314 0.411039 0.130482 0.395272 0.259641 8192 3.14159 0.764593 0.803994 0.376983 0.131866 0.406815 0.233586 16384 3.14159 0.794883 0.783260 0.399635 0.129735 0.394531 0.250665 32768 3.14159 0.787061 0.787079 0.393535 0.131700 0.396006 0.245808 65536 3.14159 0.783378 0.786510 0.391457 0.131146 0.393037 0.244640 Exact 3.14159 0.785398 0.785398 0.392699 0.130900 0.392699 0.245437 TEST01 Use DISK01_SAMPLE to estimate integrals in the unit disk. N 1 X^2 Y^2 X^4 X^2Y^2 Y^4 X^6 1 3.14159 1.17174 1.09233 0.437033 0.407413 0.379800 0.163003 2 3.14159 0.871073 1.31193 0.242026 0.362493 0.551058 0.673857E-01 4 3.14159 2.70979 1.31659 4.23969 1.98268 0.952633 7.30609 8 3.14159 3.44958 0.559704 9.60503 0.592573 0.310272 31.0658 16 3.14159 3.95475 0.814599 8.13384 1.01169 0.398758 19.9544 32 3.14159 5.02657 0.417648 12.2058 0.599819 0.142484 33.7929 64 3.14159 3.61536 0.885282 6.92055 1.07818 0.479374 16.5801 128 3.14159 4.11572 0.753432 9.16833 0.852023 0.384410 25.0573 256 3.14159 3.86404 0.747865 7.97112 0.884728 0.342426 19.7820 512 3.14159 4.17433 0.795109 8.93209 0.957347 0.390947 22.8230 1024 3.14159 3.75023 0.772190 7.75929 0.857755 0.370021 19.8001 2048 3.14159 3.84176 0.823591 7.97936 0.942209 0.422462 20.2676 4096 3.14159 3.88096 0.814577 8.07981 0.932142 0.414430 20.5243 8192 3.14159 3.89072 0.797961 8.12942 0.924753 0.403171 20.6715 16384 3.14159 3.92817 0.779070 8.23970 0.901788 0.388694 21.0172 32768 3.14159 3.92116 0.787235 8.24206 0.918455 0.393154 21.0496 65536 3.14159 3.91211 0.783140 8.20211 0.909828 0.390520 20.9178 TEST01 Use DISK01_SAMPLE to estimate integrals in the unit disk. N 1 X^2 Y^2 X^4 X^2Y^2 Y^4 X^6 1 28.2743 221.234 162.697 1731.05 1273.03 936.194 13544.7 2 28.2743 129.466 271.178 616.959 997.010 5080.19 3046.14 4 28.2743 118.306 147.129 1540.41 226.823 2231.46 22286.9 8 28.2743 95.2499 182.494 675.380 498.023 1834.81 5358.36 16 28.2743 64.7202 162.045 336.372 378.035 2310.03 2422.12 32 28.2743 88.3711 144.332 728.115 269.261 1573.07 8281.84 64 28.2743 108.618 185.435 895.948 807.857 2313.85 9143.56 128 28.2743 80.7325 171.504 577.148 499.360 2108.05 5557.30 256 28.2743 89.1339 177.973 661.456 487.932 2385.09 6483.48 512 28.2743 93.0096 167.291 709.297 461.167 2047.15 6995.99 1024 28.2743 91.1356 177.754 693.564 520.713 2299.00 6879.36 2048 28.2743 89.6701 177.669 675.905 516.875 2292.99 6715.88 4096 28.2743 90.4137 175.641 683.111 517.159 2259.59 6762.54 8192 28.2743 92.1832 176.351 702.989 530.374 2268.89 6967.15 16384 28.2743 91.3968 177.616 694.283 521.874 2281.64 6887.30 32768 28.2743 92.0556 175.834 700.002 524.457 2247.70 6940.51 65536 28.2743 91.4567 177.979 690.670 527.972 2291.80 6820.60 DISK_MONTE_CARLO_TEST Normal end of execution. 8 September 2021 7:13:19.100 PM