12 January 2014 12:04:39 AM
CIRCLE_MONTE_CARLO_PRB
C++ version
Test the CIRCLE_MONTE_CARLO library.
TEST01
Use CIRCLE01_SAMPLE to estimate integrals
along the circumference of the unit circle in 2D.
N 1 X^2 Y^2 X^4 X^2Y^2 Y^4 X^6
1 6.28319 3.76091 2.52228 2.25116 1.50975 1.01253 1.34747
2 6.28319 5.40608 0.877102 4.74323 0.662853 0.214249 4.23909
4 6.28319 2.41889 3.8643 1.8059 0.612981 3.25132 1.54901
8 6.28319 3.35786 2.92533 2.66522 0.692635 2.23269 2.31146
16 6.28319 4.48599 1.7972 3.68707 0.798921 0.998278 3.22227
32 6.28319 3.11029 3.1729 2.21406 0.896234 2.27666 1.75973
64 6.28319 3.44005 2.84314 2.64112 0.79893 2.04421 2.21612
128 6.28319 2.83891 3.44428 2.13999 0.69892 2.74536 1.8102
256 6.28319 3.03401 3.24917 2.25162 0.782389 2.46678 1.88661
512 6.28319 3.08356 3.19962 2.32384 0.759717 2.43991 1.95132
1024 6.28319 3.19346 3.08972 2.39553 0.797935 2.29179 1.99423
2048 6.28319 3.20166 3.08152 2.41081 0.790848 2.29067 2.00347
4096 6.28319 3.15973 3.12346 2.37218 0.787543 2.33592 1.97909
8192 6.28319 3.12068 3.16251 2.33731 0.783364 2.37914 1.94737
16384 6.28319 3.1149 3.16829 2.331 0.783893 2.3844 1.94203
32768 6.28319 3.12638 3.1568 2.34057 0.785816 2.37099 1.94945
65536 6.28319 3.14116 3.14202 2.35689 0.78427 2.35775 1.96457
Exact 6.28319 3.14159 3.14159 2.35619 0.785398 2.35619 1.9635
CIRCLE_MONTE_CARLO_PRB
Normal end of execution.
12 January 2014 12:04:39 AM