5 September 2021   5:06:15.680 PM

circle_rule_test():
  FORTRAN90 version
  Test circle_rule().
 
test01():
  circle_rule() computes 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.96350         1.96350    
    2    2    0.785398        0.785398    
    2    4    0.392699        0.392699    
    2    6    0.245437        0.196350    
    4    4    0.147262        0.196350    
    4    6    0.736311E-01    0.981748E-01
    6    6    0.306796E-01    0.490874E-01
 
test01():
  circle_rule() computes 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.96350         1.96350    
    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.736311E-01    0.736311E-01
    6    6    0.306796E-01    0.306796E-01

circle_rule_test():
  Normal end of execution.

 5 September 2021   5:06:15.680 PM