28 February 2022 09:03:01 PM

DISK01_RULE:
  C++ version
  Test the DISK01_RULE library.

TEST01
  DISK01_RULE can compute a rule Q(f) for the unit disk
  using NT equally spaced angles and NR radial distances.

  NT = 8
  NR = 4

  Estimate integrals I(f) where f = x^e[0] * y^e[1].

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

   0   0         3.14159         3.14159
   0   2        0.785398        0.785398
   0   4        0.392699        0.392699
   0   6        0.245437        0.245437
   2   2          0.1309          0.1309
   2   4       0.0490874       0.0490874
   2   6       0.0245437        0.019635
   4   4       0.0147262        0.019635
   4   6      0.00613592      0.00818123
   6   6       0.0021914      0.00350624

DISK01_RULE:
  Normal end of execution.

28 February 2022 09:03:01 PM