Home License -- for personal use only.  Not for government,
	academic, research, commercial, or other organizational use.

13-May-2025 13:09:21

disk01_rule_test():
  MATLAB/Octave version 9.11.0.2358333 (R2021b) Update 7
  Test DISK01_RULE().

disk01_rule_test01():
  disk01_rule() computes 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.

13-May-2025 13:09:21