16 July 2023   8:28:51.651 AM

square_symq_rule_test():
  FORTRAN90 version
  Test square_symq_rule().

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

  Number of nodes N =      7

     J  W       X       Y

     1    0.794152E-01    0.508879        0.170357E-01
     2    0.794152E-01    0.491121        0.982964    
     3    0.137950        0.889436        0.214215    
     4    0.137950        0.110564        0.785785    
     5    0.139778        0.114811        0.208516    
     6    0.139778        0.885189        0.791484    
     7    0.285714        0.500000        0.500000    
  Sum of weight:     1.00000    

test02():
  Get a quadrature rule for the symmetric square.
  Test its accuracy.
  Polynomial precision =    5

  Degree  Maximum error

   0     0.000000000000000    
   1     0.000000000000000    
   2    0.1110223024625157E-15
   3    0.5551115123125783E-16
   4    0.5551115123125783E-16
   5    0.8326672684688674E-16
   6    0.4624560744359296E-03
   7    0.1248225890155455E-02

test02():
  Test the precision of quadrature rules for the unit quadrilateral,
  Check rules of precision p = 0 through 20
  for error in approximating integrals of monomials.

              maximum                   maximum
   p          absolute                  relative
              error                     error

   0     0.000000000000000         0.000000000000000    
   1     0.000000000000000         0.000000000000000    
   2    0.5551115123125783E-16    0.1665334536937735E-15
   3    0.5551115123125783E-16    0.2220446049250313E-15
   4    0.5551115123125783E-16    0.2775557561562891E-15
   5    0.1110223024625157E-15    0.4996003610813204E-15
   6    0.1110223024625157E-15    0.3330669073875470E-15
   7    0.1110223024625157E-15    0.4163336342344337E-15
   8    0.4440892098500626E-15    0.4440892098500626E-15
   9    0.4440892098500626E-15    0.4440892098500626E-15
  10    0.1110223024625157E-15    0.6245004513516506E-15
  11    0.2220446049250313E-15    0.7494005416219807E-15
  12    0.4440892098500626E-15    0.6106226635438361E-15
  13    0.4440892098500626E-15    0.4996003610813204E-15
  14    0.3330669073875470E-15    0.1360023205165817E-14
  15    0.1110223024625157E-15    0.7632783294297951E-15
  16    0.6661338147750939E-15    0.1217775880135718E-14
  17    0.2220446049250313E-15    0.1998401444325282E-14
  18    0.1665334536937735E-15    0.9367506770274758E-15
  19    0.1110223024625157E-15    0.1498801083243961E-14
  20    0.5551115123125783E-15    0.2914335439641036E-14

square_symq_rule_test():
  Normal end of execution.

16 July 2023   8:28:51.652 AM