16-Jul-2023 18:20:10

square_symq_rule_test():
  MATLAB/Octave version 5.2.0
  Test square_symq_rule().

square_symq_rule_test01():
  Symmetric quadrature rule for a square.
  Precision = 5

     J  W       X       Y

     1       0.0794152        0.508879       0.0170357
     2       0.0794152        0.491121        0.982964
     3         0.13795        0.889436        0.214215
     4         0.13795        0.110564        0.785785
     5        0.139778        0.114811        0.208516
     6        0.139778        0.885189        0.791484
     7        0.285714             0.5             0.5
  Weight sum               1

square_symq_rule_test02():
  Test a quadrature rule for the square of precision P.

  Stated precision of rule    = 5
  Number of quadrature points = 7

  Degree  Maximum error

   0                         0
   1     1.110223024625157e-16
   2     1.665334536937735e-16
   3     5.551115123125783e-17
   4     8.326672684688674e-17
   5     8.326672684688674e-17
   6     0.0004624560744359296
   7      0.001248225890155469

square_symq_rule_test03():
  Check rules of precision p = 0 through 20
  for error in approximating integrals of monomials.

              maximum                   maximum
   p          absolute                  relative
              error                     error
   0                         0                         0
   1                         0                         0
   2     5.551115123125783e-17     1.665334536937735e-16
   3     5.551115123125783e-17      3.33066907387547e-16
   4     1.110223024625157e-16     4.440892098500626e-16
   5     1.665334536937735e-16     4.996003610813204e-16
   6     1.110223024625157e-16      3.33066907387547e-16
   7     1.110223024625157e-16     4.996003610813204e-16
   8     2.220446049250313e-16     4.440892098500626e-16
   9     2.220446049250313e-16     4.440892098500626e-16
  10     5.551115123125783e-17     3.885780586188048e-16
  11     1.110223024625157e-16     5.828670879282072e-16
  12     2.220446049250313e-16     5.412337245047638e-16
  13     2.220446049250313e-16     4.683753385137379e-16
  14     1.110223024625157e-16     8.326672684688674e-16
  15     5.551115123125783e-17     7.216449660063518e-16
  16     6.661338147750939e-16     1.165734175856414e-15
  17     1.249000902703301e-16     1.873501354054952e-15
  18     1.110223024625157e-16     5.828670879282072e-16
  19     5.551115123125783e-17     1.144917494144693e-15
  20     2.220446049250313e-16     2.518818487118324e-15

square_symq_rule_test():
  Normal end of execution.

16-Jul-2023 18:20:10