== FGM(1/2)
N= 1  D_1(A,B_N)=0.300000000 (two resolutions differ by 5.6e-17)  N*D_1=0.30000  ||H-P_N H||=0.041667 (exact |theta|/(12N))  bound=1.041667  OK
N= 2  D_1(A,B_N)=0.153377590 (two resolutions differ by 5.6e-17)  N*D_1=0.30676  ||H-P_N H||=0.020833 (exact |theta|/(12N))  bound=0.520833  OK
N= 4  D_1(A,B_N)=0.077001089 (two resolutions differ by 1.4e-17)  N*D_1=0.30800  ||H-P_N H||=0.010417 (exact |theta|/(12N))  bound=0.260417  OK
N= 8  D_1(A,B_N)=0.038534989 (two resolutions differ by 6.9e-18)  N*D_1=0.30828  ||H-P_N H||=0.005208 (exact |theta|/(12N))  bound=0.130208  OK
N=16  D_1(A,B_N)=0.019271644 (two resolutions differ by 1.0e-17)  N*D_1=0.30835  ||H-P_N H||=0.002604 (exact |theta|/(12N))  bound=0.065104  OK
N=32  D_1(A,B_N)=0.009636336 (two resolutions differ by 1.9e-17)  N*D_1=0.30836  ||H-P_N H||=0.001302 (exact |theta|/(12N))  bound=0.032552  OK
N=64  D_1(A,B_N)=0.004818232 (two resolutions differ by 1.6e-17)  N*D_1=0.30837  ||H-P_N H||=0.000651 (exact |theta|/(12N))  bound=0.016276  OK
== Clayton(2)
N= 1  D_1(A,B_N)=0.171572465 (two resolutions differ by 1.7e-07)  N*D_1=0.17157  ||H-P_N H||=0.180066 (numeric)  bound=1.180066  OK
N= 2  D_1(A,B_N)=0.125230729 (two resolutions differ by 2.1e-08)  N*D_1=0.25046  ||H-P_N H||=0.094417 (numeric)  bound=0.594417  OK
N= 4  D_1(A,B_N)=0.068430084 (two resolutions differ by 1.3e-08)  N*D_1=0.27372  ||H-P_N H||=0.047542 (numeric)  bound=0.297542  OK
N= 8  D_1(A,B_N)=0.035223255 (two resolutions differ by 1.8e-08)  N*D_1=0.28179  ||H-P_N H||=0.023684 (numeric)  bound=0.148684  OK
N=16  D_1(A,B_N)=0.017822913 (two resolutions differ by 2.9e-09)  N*D_1=0.28517  ||H-P_N H||=0.011793 (numeric)  bound=0.074293  OK
N=32  D_1(A,B_N)=0.008959027 (two resolutions differ by 1.3e-09)  N*D_1=0.28669  ||H-P_N H||=0.005880 (numeric)  bound=0.037130  OK
N=64  D_1(A,B_N)=0.004490660 (two resolutions differ by 1.6e-10)  N*D_1=0.28740  ||H-P_N H||=0.002935 (numeric)  bound=0.018560  OK
