27 February 2022 08:12:18 AM

CIRCLE_RULE:
  C++ version
  Test the CIRCLE_RULE library.

TEST01
  CIRCLE_RULE can compute a rule Q(f) for the unit circle
  using NT equally spaced angles.
  Estimate integrals I(f) where f = x^e(1) * y^e(2)
  using 8 points.

  E(1)  E(2)    I(f)            Q(f)

   0   0         6.28319         6.28319
   0   2         3.14159         3.14159
   0   4         2.35619         2.35619
   0   6          1.9635          1.9635
   2   2        0.785398        0.785398
   2   4        0.392699        0.392699
   2   6        0.245437         0.19635
   4   4        0.147262         0.19635
   4   6       0.0736311       0.0981748
   6   6       0.0306796       0.0490874

TEST01
  CIRCLE_RULE can compute a rule Q(f) for the unit circle
  using NT equally spaced angles.
  Estimate integrals I(f) where f = x^e(1) * y^e(2)
  using 32 points.

  E(1)  E(2)    I(f)            Q(f)

   0   0         6.28319         6.28319
   0   2         3.14159         3.14159
   0   4         2.35619         2.35619
   0   6          1.9635          1.9635
   2   2        0.785398        0.785398
   2   4        0.392699        0.392699
   2   6        0.245437        0.245437
   4   4        0.147262        0.147262
   4   6       0.0736311       0.0736311
   6   6       0.0306796       0.0306796

CIRCLE_RULE:
  Normal end of execution.

27 February 2022 08:12:18 AM