9 May 2025 9:00:25.481 PM disk_monte_carlo_test(): FORTRAN90 version Test disk_monte_carlo(). DISK_AREA_TEST DISK_AREA computes the area of a disk with center = (CX,CY) and radius R. ( CX CY ) R Area (-4.444955, -4.603176) 0.064275 0.012979 ( 3.928128, 4.315105) 0.406207 0.518376 ( 1.166186, -2.704491) 0.950446 2.837953 ( 3.720073, -2.212858) 0.509033 0.814032 ( 4.446055, 3.786452) 0.586273 1.079815 ( 1.736302, -0.020314) 0.861931 2.333966 ( 2.432094, -3.492716) 0.927166 2.700631 (-1.128012, 1.600010) 0.779704 1.909895 ( 4.108766, -2.350324) 0.214266 0.144230 (-2.420354, 0.047108) 0.556278 0.972152 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.14265 0.352008 0.415598 0.128031 0.394416E-01 0.151159 2 3.14159 0.175030 2.62186 0.179926E-01 0.141015 2.19122 0.192071E-02 4 3.14159 1.07047 0.575111 0.705048 0.125059 0.308347 0.566252 8 3.14159 1.07188 0.673041 0.560380 0.248453 0.219637 0.340508 16 3.14159 1.11401 0.887099 0.627213 0.187360 0.440322 0.425920 32 3.14159 0.860949 0.915702 0.402123 0.193360 0.451251 0.220924 64 3.14159 0.803142 0.687447 0.399910 0.100252 0.317661 0.247373 128 3.14159 0.815670 0.700538 0.435243 0.109496 0.334707 0.286917 256 3.14159 0.791295 0.784017 0.401528 0.142043 0.371758 0.254697 512 3.14159 0.774548 0.711726 0.395893 0.123802 0.338708 0.253802 1024 3.14159 0.835203 0.735297 0.434244 0.132272 0.358460 0.278785 2048 3.14159 0.794636 0.802520 0.390897 0.139512 0.408691 0.238131 4096 3.14159 0.767911 0.796681 0.383667 0.129696 0.397742 0.240188 8192 3.14159 0.793441 0.792537 0.396929 0.133074 0.396119 0.246598 16384 3.14159 0.784861 0.783852 0.391677 0.130111 0.391987 0.244174 32768 3.14159 0.776626 0.786689 0.387999 0.128888 0.395187 0.242757 65536 3.14159 0.783331 0.788960 0.391964 0.131150 0.395748 0.245342 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 0.389256E-01 0.738759E-01 0.482303E-03 0.915351E-03 0.173722E-02 0.597592E-05 2 3.14159 5.67782 0.150744 13.6876 0.273490 0.723356E-02 37.1215 4 3.14159 4.53426 0.819960 8.15248 1.27704 0.328480 17.2389 8 3.14159 4.62500 1.12691 8.64541 1.36751 0.623689 17.5354 16 3.14159 2.65965 0.708021 3.68707 0.764988 0.280376 5.86674 32 3.14159 3.29508 0.615134 6.77311 0.686644 0.247071 17.7986 64 3.14159 3.93366 0.932487 8.09027 0.982388 0.483021 20.5182 128 3.14159 3.91970 0.833181 8.24976 0.915409 0.465269 20.7575 256 3.14159 3.70273 0.781037 7.89593 0.826740 0.406540 20.8396 512 3.14159 3.81832 0.793448 8.01219 0.878812 0.391673 20.5694 1024 3.14159 3.87309 0.783108 8.24306 0.867402 0.387379 21.2986 2048 3.14159 3.77867 0.817531 7.75368 0.919841 0.420369 19.5522 4096 3.14159 3.89487 0.771933 8.18030 0.869820 0.386445 21.0158 8192 3.14159 3.94454 0.791851 8.25077 0.933799 0.399511 20.9338 16384 3.14159 3.95410 0.786885 8.35440 0.923525 0.392547 21.4164 32768 3.14159 3.95019 0.791114 8.33752 0.925382 0.398061 21.3885 65536 3.14159 3.91652 0.794122 8.20829 0.923718 0.399048 20.9292 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 180.991 315.305 1158.56 2018.34 3516.17 7416.22 2 28.2743 179.361 95.3736 2273.41 748.402 339.814 28829.5 4 28.2743 74.5576 146.057 653.122 37.8063 1323.12 6243.80 8 28.2743 111.106 174.542 1061.58 648.111 1933.04 12104.3 16 28.2743 48.6734 136.633 197.785 182.715 1156.86 1237.21 32 28.2743 81.8975 231.010 509.474 629.874 3396.77 3886.43 64 28.2743 74.5991 176.497 507.054 407.932 2374.48 4562.75 128 28.2743 90.1330 189.094 720.662 575.534 2584.60 7594.18 256 28.2743 92.9480 163.470 762.120 452.774 2017.30 8119.14 512 28.2743 80.3135 190.313 581.963 469.850 2568.32 5591.07 1024 28.2743 92.8593 184.266 701.154 564.058 2393.47 6844.82 2048 28.2743 90.3260 181.411 677.679 514.902 2363.36 6695.53 4096 28.2743 94.3920 180.906 724.849 556.881 2342.65 7217.89 8192 28.2743 90.3731 175.929 674.110 512.716 2224.66 6625.09 16384 28.2743 90.7002 175.961 686.984 512.857 2260.14 6809.57 32768 28.2743 91.5806 177.291 691.157 529.832 2266.89 6804.13 65536 28.2743 92.1312 176.665 700.455 528.243 2266.85 6949.51 DISK_MONTE_CARLO_TEST(): Normal end of execution. 9 May 2025 9:00:25.545 PM