08-Jan-2022 10:03:20

square_arbq_rule_test():
  MATLAB/Octave version 9.8.0.1380330 (R2020a) Update 2
  Test square_arbq_rule().

SQUARE_ARBQ_RULE_TEST01
  Symmetric quadrature rule for a square.
  Polynomial exactness degree DEGREE = 8

  Number of nodes N = 16

     J  W       X       Y

     1        0.204246       -0.227222        0.870315
     2       0.0736368        0.278678        0.985626
     3        0.193156        0.921572         0.22241
     4        0.160791       -0.522943       -0.928226
     5        0.163512        0.830917        0.843511
     6        0.310286       -0.608025        0.582595
     7       0.0646316       -0.982255       -0.821127
     8        0.429943       0.0495947       -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.0870481       -0.885013        0.909038
    13        0.570953       -0.193466       0.0352632
    14         0.12035        0.577245       -0.962256
    15        0.145236        0.921307       -0.708268
    16        0.377116       -0.717604       -0.413062
   Sum               4
  Area               4

SQUARE_ARBQ_RULE_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"

SQUARE_ARBQ_RULE_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"

SQUARE_ARBQ_RULE_TEST04
  Get a quadrature rule for the symmetric square.
  Test its accuracy.
  Polynomial exactness degree DEGREE = 8

  RMS error = 4.22762e-17

square_arbq_rule_test():
  Normal end of execution.

08-Jan-2022 10:03:20