Wed Oct 8 07:30:27 2025 cube_integrals_test(): python version: 3.10.12 numpy version: 1.26.4 Test cube_integrals(). cube01_monomial_integral_test(): cube01_monomial_integral() computes the integral of a monomial within the interior of the unit cube in 3D. Compare with a Monte Carlo estimate. Number of sample points used is 4192 Ex Ey Ez MC-Estimate Exact Error 0 4 5 0.0344122 0.0333333 0.0011 1 4 4 0.0205156 0.02 0.00052 1 5 0 0.087186 0.0833333 0.0039 3 6 1 0.0191193 0.0178571 0.0013 2 6 6 0.00715501 0.00680272 0.00035 6 6 6 0.00316827 0.00291545 0.00025 5 2 7 0.00694216 0.00694444 2.3e-06 4 0 1 0.0977997 0.1 0.0022 6 0 0 0.141858 0.142857 0.001 1 0 2 0.162358 0.166667 0.0043 2 2 1 0.0562364 0.0555556 0.00068 6 4 6 0.00432144 0.00408163 0.00024 3 4 1 0.0261835 0.025 0.0012 0 0 6 0.137149 0.142857 0.0057 4 1 5 0.0161658 0.0166667 0.0005 0 1 2 0.166055 0.166667 0.00061 5 7 6 0.00324801 0.00297619 0.00027 0 5 4 0.0348004 0.0333333 0.0015 7 2 7 0.00527166 0.00520833 6.3e-05 2 3 6 0.0119708 0.0119048 6.6e-05 cube01_sample_test cube01_sample() samples the unit cube. Sample points in the unit cube. Row: 0 1 2 Col 0 : 0.169125 0.740461 0.187758 1 : 0.358098 0.432466 0.768925 2 : 0.150303 0.940819 0.814896 3 : 0.0949028 0.869298 0.349022 4 : 0.385498 0.687337 0.383212 5 : 0.0234536 0.506725 0.728244 6 : 0.216068 0.159416 0.402048 7 : 0.854395 0.69031 0.67979 8 : 0.403015 0.579655 0.679103 9 : 0.294756 0.961644 0.615026 cube01_volume_test() cube01_volume() returns the volume of the unit cube. cube01_volume() = 1.0 cube_integrals_test(): Normal end of execution. Wed Oct 8 07:30:28 2025