3 April 2023   9:06:54.286 AM

square_arbq_rule_test():
  FORTRAN90 version
  Test SQUARE_ARBQ_RULE().

TEST01
  Quadrature rule for a symmetric square.
  Polynomial exactness degree DEGREE =    8

  Number of nodes N =     16

     J  W       X       Y

     1    0.204246       -0.227222        0.870315    
     2    0.736368E-01    0.278678        0.985626    
     3    0.193156        0.921572        0.222410    
     4    0.160791       -0.522943       -0.928226    
     5    0.163512        0.830917        0.843511    
     6    0.310286       -0.608025        0.582595    
     7    0.646316E-01   -0.982255       -0.821127    
     8    0.429943        0.495947E-01   -0.691724    
     9    0.456476        0.591001       -0.261441    
    10    0.472514        0.362659        0.519812    
    11    0.170105       -0.936916        0.215377    
    12    0.870481E-01   -0.885013        0.909038    
    13    0.570953       -0.193466        0.352632E-01
    14    0.120350        0.577245       -0.962256    
    15    0.145236        0.921307       -0.708268    
    16    0.377116       -0.717604       -0.413062    
   Sum     4.00000    
  Area     4.00000    

TEST02
  Get a quadrature rule for the symmetric square.
  Then write it to a file.
  Polynomial exactness degree DEGREE =    8

  Quadrature rule written to file "square08.txt".

TEST03
  Get a quadrature rule for the symmetric square.
  Set up GNUPLOT graphics input.
  Polynomial exactness degree DEGREE =    8
 
  Created vertex file "square08_vertices.txt".
  Created node file "square08_nodes.txt".
  Created command file "square08_commands.txt".

TEST04
  Get a quadrature rule for the symmetric square.
  Test its accuracy.
  Polynomial exactness degree DEGREE =    8
 
  RMS error =   0.422762E-16

SQUARE_ARBQ_RULE_TEST
  Normal end of execution.

 3 April 2023   9:06:54.286 AM