5 July 2018 8:23:26.823 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 (-2.815817, 4.563176) 0.829509 2.161685 ( 0.616954, -0.846929) 0.066119 0.013734 (-2.424222, -3.900432) 0.043829 0.006035 ( 1.339657, -4.382728) 0.449539 0.634870 (-0.986937, 2.546735) 0.797287 1.997005 (-4.981616, 3.975041) 0.350752 0.386501 (-4.054552, -4.863831) 0.859097 2.318644 ( 3.408475, -3.768961) 0.007512 0.000177 (-2.396970, 4.124837) 0.113664 0.040588 (-1.483713, 3.228873) 0.267132 0.224183 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 2.41455 0.191431 1.85576 0.147129 0.116647E-01 1.42629 2 3.14159 0.118014 0.546915E-01 0.485103E-02 0.204044E-02 0.952587E-03 0.213620E-03 4 3.14159 1.14909 0.780379 0.727307 0.463617E-01 0.514225 0.540924 8 3.14159 0.664859 0.444275 0.328095 0.823115E-01 0.104648 0.200273 16 3.14159 0.667844 0.675045 0.299232 0.117983 0.322511 0.167483 32 3.14159 0.874306 0.836608 0.388713 0.146414 0.451474 0.205617 64 3.14159 0.732941 0.804583 0.353860 0.111915 0.434813 0.218527 128 3.14159 0.756945 0.816620 0.396898 0.108539 0.430483 0.261179 256 3.14159 0.767384 0.752709 0.377832 0.122918 0.381345 0.233849 512 3.14159 0.746584 0.780966 0.357190 0.125881 0.391175 0.213426 1024 3.14159 0.795306 0.829076 0.401273 0.143614 0.429817 0.254394 2048 3.14159 0.759074 0.776137 0.373097 0.131972 0.377560 0.228936 4096 3.14159 0.783685 0.787455 0.392849 0.128952 0.395614 0.245354 8192 3.14159 0.791588 0.780931 0.397664 0.129168 0.389458 0.249651 16384 3.14159 0.785739 0.782836 0.391171 0.130587 0.392443 0.242876 32768 3.14159 0.787655 0.782700 0.395488 0.129711 0.391104 0.248700 65536 3.14159 0.790243 0.774556 0.397422 0.129430 0.384659 0.249120 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 11.0645 0.191431 38.9686 0.674208 0.116647E-01 137.245 2 3.14159 3.07038 0.546915E-01 3.43411 0.539039E-01 0.952587E-03 4.20327 4 3.14159 7.50138 0.780379 20.7441 1.02540 0.514225 62.9682 8 3.14159 3.47952 0.444275 6.02262 0.528037 0.104648 12.4222 16 3.14159 3.98698 0.675045 7.60476 0.735872 0.322511 16.9888 32 3.14159 2.97603 0.836608 5.80072 0.928332 0.451474 13.9917 64 3.14159 3.57716 0.804583 7.39384 0.827737 0.434813 19.4895 128 3.14159 3.86149 0.816620 7.85101 0.968444 0.430483 19.5214 256 3.14159 4.01683 0.752709 8.38095 0.910712 0.381345 21.1228 512 3.14159 4.01407 0.780966 8.28044 0.929975 0.391175 20.5420 1024 3.14159 3.81316 0.829076 7.86037 0.963054 0.429817 19.6905 2048 3.14159 3.86571 0.776137 7.99053 0.893014 0.377560 20.2063 4096 3.14159 3.82103 0.787455 7.94217 0.903813 0.395614 20.2283 8192 3.14159 3.93695 0.780931 8.29047 0.912874 0.389458 21.1994 16384 3.14159 3.92407 0.782836 8.25722 0.912614 0.392443 21.1283 32768 3.14159 3.91229 0.782700 8.21297 0.910051 0.391104 20.9928 65536 3.14159 3.94683 0.774556 8.31645 0.904323 0.384659 21.2827 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 372.579 44.8494 4909.57 590.994 71.1413 64694.8 2 28.2743 32.7243 162.292 68.7500 188.458 931.549 151.040 4 28.2743 208.039 98.6667 2155.21 984.286 594.652 25628.1 8 28.2743 73.3008 69.1533 382.809 83.2799 417.295 2745.65 16 28.2743 87.1634 191.173 532.407 635.184 2509.51 4179.46 32 28.2743 71.0167 117.424 462.755 346.802 1220.98 3980.62 64 28.2743 79.6134 183.961 649.594 461.808 2448.54 6843.33 128 28.2743 88.5867 159.472 645.808 439.722 2093.85 6148.51 256 28.2743 93.3445 192.772 692.574 596.284 2564.05 6768.89 512 28.2743 92.1468 182.394 670.543 512.893 2375.75 6314.85 1024 28.2743 89.3531 174.978 650.375 486.690 2323.49 6246.40 2048 28.2743 88.8154 182.211 665.173 526.013 2354.31 6522.38 4096 28.2743 88.9383 175.834 669.030 504.416 2233.00 6632.29 8192 28.2743 92.4947 179.022 701.665 537.111 2304.05 6941.95 16384 28.2743 91.8312 178.178 699.149 525.112 2292.65 6911.65 32768 28.2743 91.6166 176.242 694.550 522.806 2254.70 6892.08 65536 28.2743 92.6888 175.473 704.319 524.418 2230.35 6981.70 DISK_MONTE_CARLO_TEST Normal end of execution. 5 July 2018 8:23:26.926 PM