Warning: X does not support locale en_US.UTF-8 14-Jul-2023 14:38:58 prism_witherden_rule_test(): MATLAB/Octave version 9.14.0.2206163 (R2023a) Test prism_witherden_rule(). prism_witherden_rule_test01(): Quadrature rule for the unit prism, Precision p = 5 I W X Y Z 1 0.177864 0.333333 0.333333 0.5 2 0.0561775 0.4873 0.0254001 0.5 3 0.0561775 0.0254001 0.4873 0.5 4 0.0561775 0.4873 0.4873 0.5 5 0.0464154 0.445598 0.108804 0.0644985 6 0.0464154 0.108804 0.445598 0.0644985 7 0.0464154 0.445598 0.445598 0.0644985 8 0.0464154 0.445598 0.108804 0.935501 9 0.0464154 0.108804 0.445598 0.935501 10 0.0464154 0.445598 0.445598 0.935501 11 0.0625186 0.100859 0.798282 0.214787 12 0.0625186 0.798282 0.100859 0.214787 13 0.0625186 0.100859 0.100859 0.214787 14 0.0625186 0.100859 0.798282 0.785213 15 0.0625186 0.798282 0.100859 0.785213 16 0.0625186 0.100859 0.100859 0.785213 Weight Sum 1 pyramid_witherden_rule_test02(): Test the precision of a quadrature rule for the unit prism, Stated precision of rule = 5 Number of quadrature points = 16 Degree Maximum error 0 5.551115123125783e-17 1 2.775557561562891e-17 2 5.551115123125783e-17 3 2.775557561562891e-17 4 5.551115123125783e-17 5 4.163336342344337e-17 6 8.681433560950996e-05 7 0.0002425567824553866 pyramid_witherden_rule_test03(): Test the precision of quadrature rules for the unit prism, Check rules of precision p = 0 through 10 for error in approximating integrals of monomials. maximum maximum p absolute relative error error 0 0 0 1 5.551115123125783e-17 3.33066907387547e-16 2 2.775557561562891e-17 1.665334536937735e-16 3 1.110223024625157e-16 3.33066907387547e-16 4 5.551115123125783e-17 4.163336342344337e-16 5 5.551115123125783e-17 6.245004513516506e-16 6 6.661338147750939e-16 1.748601263784622e-15 7 1.665334536937735e-16 1.332267629550188e-15 8 2.220446049250313e-16 5.854691731421724e-16 9 5.551115123125783e-17 1.554312234475219e-15 10 2.220446049250313e-16 2.441894342980477e-15 prism_witherden_rule_test(): Normal end of execution. 14-Jul-2023 14:38:59