07-Jan-2022 19:05:12

disk01_rule_test():
  MATLAB/Octave version 9.8.0.1380330 (R2020a) Update 2
  Test DISK01_RULE().

DISK01_RULE_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(1) * y^e(2).

  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.

07-Jan-2022 19:05:12