23 August 2012 10:18:52.907 AM GM_RULE_PRB FORTRAN77 version Test the GM_RULE library. TEST01 SIMPLEX_UNIT_TO_GENERAL maps points in the unit simplex to a general simplex. Here we consider a simplex in 2D, a triangle. The vertices of the general triangle are: 1.0000 1.0000 3.0000 1.0000 2.0000 5.0000 ( XSI ETA ) ( X Y ) 0.0000 1.0000 0.0000 1.0000 1.0000 3.0000 0.0000 1.0000 0.0000 2.0000 1.0000 5.0000 0.8679 2.7613 0.0255 1.1019 0.1383 1.4872 0.2106 1.8425 0.2027 1.7353 0.3299 2.3197 0.1128 1.9149 0.6893 3.7572 0.6425 2.4831 0.1981 1.7923 0.8450 2.7044 0.0145 1.0580 0.3465 2.3242 0.6312 3.5247 0.0242 1.3410 0.2926 2.1704 0.3726 1.7706 0.0254 1.1014 0.4083 1.8926 0.0761 1.3046 TEST02 SIMPLEX_UNIT_TO_GENERAL maps points in the unit simplex to a general simplex. Here we consider a simplex in 3D, a tetrahedron. The vertices are: 1.0000 1.0000 1.0000 3.0000 1.0000 1.0000 1.0000 4.0000 1.0000 1.0000 1.0000 5.0000 ( XSI ETA MU ) ( X Y Z ) 0.0000 1.0000 0.0000 1.0000 0.0000 1.0000 1.0000 3.0000 0.0000 1.0000 0.0000 1.0000 0.0000 1.0000 1.0000 4.0000 0.0000 1.0000 0.0000 1.0000 0.0000 1.0000 1.0000 5.0000 0.6530 2.3060 0.0192 1.0575 0.0802 1.3209 0.1227 1.2455 0.3794 2.1383 0.1895 1.7579 0.4363 1.8726 0.0636 1.1908 0.3885 2.5542 0.1183 1.2365 0.0365 1.1094 0.0293 1.1174 0.0138 1.0277 0.1341 1.4024 0.3020 2.2079 0.0208 1.0415 0.0237 1.0711 0.2865 2.1460 0.2890 1.5780 0.0197 1.0590 0.4669 2.8677 0.0792 1.1585 0.5366 2.6099 0.1496 1.5985 0.0966 1.1933 0.5111 2.5332 0.0596 1.2384 0.3663 1.7327 0.0599 1.1797 0.2031 1.8125 TEST03 GM_RULE_SIZE returns POINT_NUM, the number of points associated with a Grundmann-Moeller quadrature rule for the unit simplex of dimension DIM_NUM with rule index RULE and degree of exactness DEGREE = 2*RULE+1. DIM_NUM RULE DEGREE POINT_NUM 2 0 1 1 2 1 3 4 2 2 5 10 2 3 7 20 2 4 9 35 2 5 11 56 3 0 1 1 3 1 3 5 3 2 5 15 3 3 7 35 3 4 9 70 3 5 11 126 5 0 1 1 5 1 3 7 5 2 5 28 5 3 7 84 5 4 9 210 5 5 11 462 10 0 1 1 10 1 3 12 10 2 5 78 10 3 7 364 10 4 9 1365 10 5 11 4368 TEST04 GM_RULE_SET determines the weights and abscissas of a Grundmann-Moeller quadrature rule for the DIM_NUM dimensional simplex, using a rule of in index RULE, which will have degree of exactness 2*RULE+1. Here we use DIM_NUM = 3 RULE = 2 DEGREE = 5 POINT W X Y Z 1 0.304762 0.125000 0.125000 0.125000 2 0.304762 0.375000 0.125000 0.125000 3 0.304762 0.625000 0.125000 0.125000 4 0.304762 0.125000 0.375000 0.125000 5 0.304762 0.375000 0.375000 0.125000 6 0.304762 0.125000 0.625000 0.125000 7 0.304762 0.125000 0.125000 0.375000 8 0.304762 0.375000 0.125000 0.375000 9 0.304762 0.125000 0.375000 0.375000 10 0.304762 0.125000 0.125000 0.625000 11 -0.578571 0.166667 0.166667 0.166667 12 -0.578571 0.500000 0.166667 0.166667 13 -0.578571 0.166667 0.500000 0.166667 14 -0.578571 0.166667 0.166667 0.500000 15 0.266667 0.250000 0.250000 0.250000 TEST05 GM_RULE_SET determines the weights and abscissas of a Grundmann-Moeller quadrature rule for the DIM_NUM dimensional simplex, using a rule of in index RULE, which will have degree of exactness 2*RULE+1. In this test, we compute various rules, and simply report the number of points, and the sum of weights. DIM_NUM RULE POINT_NUM WEIGHT SUM 2 0 1 1.000000000000000 2 1 4 0.9999999999999998 2 2 10 0.9999999999999998 2 3 20 1.000000000000001 2 4 35 0.9999999999999998 2 5 56 1.000000000000006 3 0 1 1.000000000000000 3 1 5 1.000000000000000 3 2 15 0.9999999999999993 3 3 35 1.000000000000002 3 4 70 0.9999999999999917 3 5 126 0.9999999999999940 5 0 1 1.000000000000000 5 1 7 0.9999999999999998 5 2 28 0.9999999999999983 5 3 84 0.9999999999999918 5 4 210 0.9999999999999998 5 5 462 1.000000000000103 10 0 1 1.000000000000000 10 1 12 1.000000000000000 10 2 78 0.9999999999999898 10 3 364 0.9999999999999676 10 4 1365 0.9999999999992277 10 5 4368 1.000000000006103 TEST06 GM_RULE_SET determines the weights and abscissas of a Grundmann-Moeller quadrature rule for the DIM_NUM dimensional simplex, using a rule of in index RULE, which will have degree of exactness 2*RULE+1. In this test, we write a rule to a file. Here we use DIM_NUM = 3 RULE = 2 DEGREE = 5 Wrote rule 2 to "gm2_3d_w.txt" and "gm2_3d_x.txt". TEST07 GM_RULE_SET determines the weights and abscissas of a Grundmann-Moeller quadrature rule for the DIM_NUM dimensional simplex, using a rule of in index RULE, which will have degree of exactness 2*RULE+1. In this test, look at all the monomials up to some maximum degree, choose a few low order rules and determine the quadrature error for each. Here we use DIM_NUM = 5 Rule Order Quad_Error F(X) = X1^0 * X2^0 * X3^0 * X4^0 * X5^0 0 1 0.00000 1 7 0.222045E-15 2 28 0.166533E-14 3 84 0.810463E-14 F(X) = X1^1 * X2^0 * X3^0 * X4^0 * X5^0 0 1 0.111022E-15 1 7 0.333067E-15 2 28 0.222045E-15 3 84 0.444089E-15 F(X) = X1^0 * X2^1 * X3^0 * X4^0 * X5^0 0 1 0.111022E-15 1 7 0.333067E-15 2 28 0.222045E-15 3 84 0.888178E-15 F(X) = X1^0 * X2^0 * X3^1 * X4^0 * X5^0 0 1 0.111022E-15 1 7 0.333067E-15 2 28 0.222045E-15 3 84 0.155431E-14 F(X) = X1^0 * X2^0 * X3^0 * X4^1 * X5^0 0 1 0.111022E-15 1 7 0.333067E-15 2 28 0.222045E-15 3 84 0.155431E-14 F(X) = X1^0 * X2^0 * X3^0 * X4^0 * X5^1 0 1 0.111022E-15 1 7 0.333067E-15 2 28 0.222045E-15 3 84 0.444089E-15 F(X) = X1^2 * X2^0 * X3^0 * X4^0 * X5^0 0 1 0.416667 1 7 0.222045E-15 2 28 0.133227E-14 3 84 0.111022E-13 F(X) = X1^1 * X2^1 * X3^0 * X4^0 * X5^0 0 1 0.166667 1 7 0.111022E-15 2 28 0.111022E-15 3 84 0.176525E-13 F(X) = X1^0 * X2^2 * X3^0 * X4^0 * X5^0 0 1 0.416667 1 7 0.222045E-15 2 28 0.133227E-14 3 84 0.988098E-14 F(X) = X1^1 * X2^0 * X3^1 * X4^0 * X5^0 0 1 0.166667 1 7 0.111022E-15 2 28 0.111022E-15 3 84 0.169864E-13 F(X) = X1^0 * X2^1 * X3^1 * X4^0 * X5^0 0 1 0.166667 1 7 0.111022E-15 2 28 0.444089E-15 3 84 0.157652E-13 F(X) = X1^0 * X2^0 * X3^2 * X4^0 * X5^0 0 1 0.416667 1 7 0.222045E-15 2 28 0.133227E-14 3 84 0.577316E-14 F(X) = X1^1 * X2^0 * X3^0 * X4^1 * X5^0 0 1 0.166667 1 7 0.111022E-15 2 28 0.444089E-15 3 84 0.146549E-13 F(X) = X1^0 * X2^1 * X3^0 * X4^1 * X5^0 0 1 0.166667 1 7 0.111022E-15 2 28 0.111022E-15 3 84 0.126565E-13 F(X) = X1^0 * X2^0 * X3^1 * X4^1 * X5^0 0 1 0.166667 1 7 0.111022E-15 2 28 0.111022E-15 3 84 0.122125E-13 F(X) = X1^0 * X2^0 * X3^0 * X4^2 * X5^0 0 1 0.416667 1 7 0.222045E-15 2 28 0.133227E-14 3 84 0.133227E-14 F(X) = X1^1 * X2^0 * X3^0 * X4^0 * X5^1 0 1 0.166667 1 7 0.111022E-15 2 28 0.111022E-15 3 84 0.766054E-14 F(X) = X1^0 * X2^1 * X3^0 * X4^0 * X5^1 0 1 0.166667 1 7 0.111022E-15 2 28 0.111022E-15 3 84 0.632827E-14 F(X) = X1^0 * X2^0 * X3^1 * X4^0 * X5^1 0 1 0.166667 1 7 0.111022E-15 2 28 0.444089E-15 3 84 0.577316E-14 F(X) = X1^0 * X2^0 * X3^0 * X4^1 * X5^1 0 1 0.166667 1 7 0.111022E-15 2 28 0.111022E-14 3 84 0.399680E-14 F(X) = X1^0 * X2^0 * X3^0 * X4^0 * X5^2 0 1 0.416667 1 7 0.111022E-15 2 28 0.222045E-15 3 84 0.732747E-14 F(X) = X1^3 * X2^0 * X3^0 * X4^0 * X5^0 0 1 0.740741 1 7 0.222045E-15 2 28 0.222045E-14 3 84 0.355271E-14 F(X) = X1^2 * X2^1 * X3^0 * X4^0 * X5^0 0 1 0.222222 1 7 0.444089E-15 2 28 0.133227E-14 3 84 0.172085E-13 F(X) = X1^1 * X2^2 * X3^0 * X4^0 * X5^0 0 1 0.222222 1 7 0.444089E-15 2 28 0.133227E-14 3 84 0.166533E-13 F(X) = X1^0 * X2^3 * X3^0 * X4^0 * X5^0 0 1 0.740741 1 7 0.222045E-15 2 28 0.166533E-14 3 84 0.421885E-14 F(X) = X1^2 * X2^0 * X3^1 * X4^0 * X5^0 0 1 0.222222 1 7 0.222045E-15 2 28 0.133227E-14 3 84 0.155431E-13 F(X) = X1^1 * X2^1 * X3^1 * X4^0 * X5^0 0 1 0.555556 1 7 0.444089E-15 2 28 0.444089E-15 3 84 0.444089E-15 F(X) = X1^0 * X2^2 * X3^1 * X4^0 * X5^0 0 1 0.222222 1 7 0.222045E-15 2 28 0.888178E-15 3 84 0.132117E-13 F(X) = X1^1 * X2^0 * X3^2 * X4^0 * X5^0 0 1 0.222222 1 7 0.222045E-15 2 28 0.133227E-14 3 84 0.137668E-13 F(X) = X1^0 * X2^1 * X3^2 * X4^0 * X5^0 0 1 0.222222 1 7 0.222045E-15 2 28 0.133227E-14 3 84 0.137668E-13 F(X) = X1^0 * X2^0 * X3^3 * X4^0 * X5^0 0 1 0.740741 1 7 0.00000 2 28 0.888178E-15 3 84 0.466294E-14 F(X) = X1^2 * X2^0 * X3^0 * X4^1 * X5^0 0 1 0.222222 1 7 0.222045E-15 2 28 0.133227E-14 3 84 0.132117E-13 F(X) = X1^1 * X2^1 * X3^0 * X4^1 * X5^0 0 1 0.555556 1 7 0.444089E-15 2 28 0.888178E-15 3 84 0.444089E-15 F(X) = X1^0 * X2^2 * X3^0 * X4^1 * X5^0 0 1 0.222222 1 7 0.222045E-15 2 28 0.888178E-15 3 84 0.107692E-13 F(X) = X1^1 * X2^0 * X3^1 * X4^1 * X5^0 0 1 0.555556 1 7 0.444089E-15 2 28 0.888178E-15 3 84 0.444089E-15 F(X) = X1^0 * X2^1 * X3^1 * X4^1 * X5^0 0 1 0.555556 1 7 0.444089E-15 2 28 0.888178E-15 3 84 0.666134E-15 F(X) = X1^0 * X2^0 * X3^2 * X4^1 * X5^0 0 1 0.222222 1 7 0.222045E-15 2 28 0.222045E-15 3 84 0.788258E-14 F(X) = X1^1 * X2^0 * X3^0 * X4^2 * X5^0 0 1 0.222222 1 7 0.222045E-15 2 28 0.888178E-15 3 84 0.910383E-14 F(X) = X1^0 * X2^1 * X3^0 * X4^2 * X5^0 0 1 0.222222 1 7 0.222045E-15 2 28 0.133227E-14 3 84 0.843769E-14 F(X) = X1^0 * X2^0 * X3^1 * X4^2 * X5^0 0 1 0.222222 1 7 0.222045E-15 2 28 0.888178E-15 3 84 0.788258E-14 F(X) = X1^0 * X2^0 * X3^0 * X4^3 * X5^0 0 1 0.740741 1 7 0.00000 2 28 0.555112E-15 3 84 0.466294E-14 F(X) = X1^2 * X2^0 * X3^0 * X4^0 * X5^1 0 1 0.222222 1 7 0.222045E-15 2 28 0.133227E-14 3 84 0.888178E-14 F(X) = X1^1 * X2^1 * X3^0 * X4^0 * X5^1 0 1 0.555556 1 7 0.444089E-15 2 28 0.888178E-15 3 84 0.444089E-15 F(X) = X1^0 * X2^2 * X3^0 * X4^0 * X5^1 0 1 0.222222 1 7 0.222045E-15 2 28 0.133227E-14 3 84 0.555112E-14 F(X) = X1^1 * X2^0 * X3^1 * X4^0 * X5^1 0 1 0.555556 1 7 0.444089E-15 2 28 0.888178E-15 3 84 0.222045E-15 F(X) = X1^0 * X2^1 * X3^1 * X4^0 * X5^1 0 1 0.555556 1 7 0.444089E-15 2 28 0.444089E-15 3 84 0.222045E-15 F(X) = X1^0 * X2^0 * X3^2 * X4^0 * X5^1 0 1 0.222222 1 7 0.222045E-15 2 28 0.888178E-15 3 84 0.155431E-14 F(X) = X1^1 * X2^0 * X3^0 * X4^1 * X5^1 0 1 0.555556 1 7 0.444089E-15 2 28 0.888178E-15 3 84 0.777156E-15 F(X) = X1^0 * X2^1 * X3^0 * X4^1 * X5^1 0 1 0.555556 1 7 0.444089E-15 2 28 0.888178E-15 3 84 0.555112E-15 F(X) = X1^0 * X2^0 * X3^1 * X4^1 * X5^1 0 1 0.555556 1 7 0.444089E-15 2 28 0.888178E-15 3 84 0.222045E-15 F(X) = X1^0 * X2^0 * X3^0 * X4^2 * X5^1 0 1 0.222222 1 7 0.222045E-15 2 28 0.888178E-15 3 84 0.310862E-14 F(X) = X1^1 * X2^0 * X3^0 * X4^0 * X5^2 0 1 0.222222 1 7 0.222045E-15 2 28 0.888178E-15 3 84 0.333067E-14 F(X) = X1^0 * X2^1 * X3^0 * X4^0 * X5^2 0 1 0.222222 1 7 0.222045E-15 2 28 0.888178E-15 3 84 0.488498E-14 F(X) = X1^0 * X2^0 * X3^1 * X4^0 * X5^2 0 1 0.222222 1 7 0.222045E-15 2 28 0.888178E-15 3 84 0.555112E-14 F(X) = X1^0 * X2^0 * X3^0 * X4^1 * X5^2 0 1 0.222222 1 7 0.222045E-15 2 28 0.888178E-15 3 84 0.555112E-14 F(X) = X1^0 * X2^0 * X3^0 * X4^0 * X5^3 0 1 0.740741 1 7 0.00000 2 28 0.222045E-15 3 84 0.288658E-14 F(X) = X1^4 * X2^0 * X3^0 * X4^0 * X5^0 0 1 0.902778 1 7 0.117188 2 28 0.222045E-14 3 84 0.111022E-13 F(X) = X1^3 * X2^1 * X3^0 * X4^0 * X5^0 0 1 0.611111 1 7 0.937500E-01 2 28 0.444089E-15 3 84 0.510703E-14 F(X) = X1^2 * X2^2 * X3^0 * X4^0 * X5^0 0 1 0.416667 1 7 0.203125 2 28 0.266454E-14 3 84 0.222045E-15 F(X) = X1^1 * X2^3 * X3^0 * X4^0 * X5^0 0 1 0.611111 1 7 0.937500E-01 2 28 0.444089E-15 3 84 0.510703E-14 F(X) = X1^0 * X2^4 * X3^0 * X4^0 * X5^0 0 1 0.902778 1 7 0.117188 2 28 0.199840E-14 3 84 0.104361E-13 F(X) = X1^3 * X2^0 * X3^1 * X4^0 * X5^0 0 1 0.611111 1 7 0.937500E-01 2 28 0.111022E-14 3 84 0.510703E-14 F(X) = X1^2 * X2^1 * X3^1 * X4^0 * X5^0 0 1 0.166667 1 7 0.312500E-01 2 28 0.421885E-14 3 84 0.222045E-14 F(X) = X1^1 * X2^2 * X3^1 * X4^0 * X5^0 0 1 0.166667 1 7 0.312500E-01 2 28 0.355271E-14 3 84 0.244249E-14 F(X) = X1^0 * X2^3 * X3^1 * X4^0 * X5^0 0 1 0.611111 1 7 0.937500E-01 2 28 0.111022E-14 3 84 0.610623E-14 F(X) = X1^2 * X2^0 * X3^2 * X4^0 * X5^0 0 1 0.416667 1 7 0.203125 2 28 0.199840E-14 3 84 0.222045E-15 F(X) = X1^1 * X2^1 * X3^2 * X4^0 * X5^0 0 1 0.166667 1 7 0.312500E-01 2 28 0.421885E-14 3 84 0.222045E-14 F(X) = X1^0 * X2^2 * X3^2 * X4^0 * X5^0 0 1 0.416667 1 7 0.203125 2 28 0.199840E-14 3 84 0.222045E-14 F(X) = X1^1 * X2^0 * X3^3 * X4^0 * X5^0 0 1 0.611111 1 7 0.937500E-01 2 28 0.444089E-15 3 84 0.466294E-14 F(X) = X1^0 * X2^1 * X3^3 * X4^0 * X5^0 0 1 0.611111 1 7 0.937500E-01 2 28 0.444089E-15 3 84 0.510703E-14 F(X) = X1^0 * X2^0 * X3^4 * X4^0 * X5^0 0 1 0.902778 1 7 0.117188 2 28 0.111022E-14 3 84 0.821565E-14 F(X) = X1^3 * X2^0 * X3^0 * X4^1 * X5^0 0 1 0.611111 1 7 0.937500E-01 2 28 0.111022E-14 3 84 0.344169E-14 F(X) = X1^2 * X2^1 * X3^0 * X4^1 * X5^0 0 1 0.166667 1 7 0.312500E-01 2 28 0.266454E-14 3 84 0.155431E-14 F(X) = X1^1 * X2^2 * X3^0 * X4^1 * X5^0 0 1 0.166667 1 7 0.312500E-01 2 28 0.355271E-14 3 84 0.177636E-14 F(X) = X1^0 * X2^3 * X3^0 * X4^1 * X5^0 0 1 0.611111 1 7 0.937500E-01 2 28 0.444089E-15 3 84 0.388578E-14 F(X) = X1^2 * X2^0 * X3^1 * X4^1 * X5^0 0 1 0.166667 1 7 0.312500E-01 2 28 0.266454E-14 3 84 0.155431E-14 F(X) = X1^1 * X2^1 * X3^1 * X4^1 * X5^0 0 1 1.33333 1 7 0.625000E-01 2 28 0.177636E-14 3 84 0.888178E-15 F(X) = X1^0 * X2^2 * X3^1 * X4^1 * X5^0 0 1 0.166667 1 7 0.312500E-01 2 28 0.266454E-14 3 84 0.222045E-14 F(X) = X1^1 * X2^0 * X3^2 * X4^1 * X5^0 0 1 0.166667 1 7 0.312500E-01 2 28 0.222045E-14 3 84 0.155431E-14 F(X) = X1^0 * X2^1 * X3^2 * X4^1 * X5^0 0 1 0.166667 1 7 0.312500E-01 2 28 0.222045E-14 3 84 0.155431E-14 F(X) = X1^0 * X2^0 * X3^3 * X4^1 * X5^0 0 1 0.611111 1 7 0.937500E-01 2 28 0.444089E-15 3 84 0.310862E-14 F(X) = X1^2 * X2^0 * X3^0 * X4^2 * X5^0 0 1 0.416667 1 7 0.203125 2 28 0.133227E-14 3 84 0.111022E-14 F(X) = X1^1 * X2^1 * X3^0 * X4^2 * X5^0 0 1 0.166667 1 7 0.312500E-01 2 28 0.266454E-14 3 84 0.888178E-15 F(X) = X1^0 * X2^2 * X3^0 * X4^2 * X5^0 0 1 0.416667 1 7 0.203125 2 28 0.133227E-14 3 84 0.111022E-14 F(X) = X1^1 * X2^0 * X3^1 * X4^2 * X5^0 0 1 0.166667 1 7 0.312500E-01 2 28 0.266454E-14 3 84 0.155431E-14 F(X) = X1^0 * X2^1 * X3^1 * X4^2 * X5^0 0 1 0.166667 1 7 0.312500E-01 2 28 0.266454E-14 3 84 0.888178E-15 F(X) = X1^0 * X2^0 * X3^2 * X4^2 * X5^0 0 1 0.416667 1 7 0.203125 2 28 0.133227E-14 3 84 0.222045E-15 F(X) = X1^1 * X2^0 * X3^0 * X4^3 * X5^0 0 1 0.611111 1 7 0.937500E-01 2 28 0.444089E-15 3 84 0.310862E-14 F(X) = X1^0 * X2^1 * X3^0 * X4^3 * X5^0 0 1 0.611111 1 7 0.937500E-01 2 28 0.111022E-14 3 84 0.388578E-14 F(X) = X1^0 * X2^0 * X3^1 * X4^3 * X5^0 0 1 0.611111 1 7 0.937500E-01 2 28 0.444089E-15 3 84 0.388578E-14 F(X) = X1^0 * X2^0 * X3^0 * X4^4 * X5^0 0 1 0.902778 1 7 0.117188 2 28 0.888178E-15 3 84 0.199840E-14 F(X) = X1^3 * X2^0 * X3^0 * X4^0 * X5^1 0 1 0.611111 1 7 0.937500E-01 2 28 0.111022E-14 3 84 0.344169E-14 F(X) = X1^2 * X2^1 * X3^0 * X4^0 * X5^1 0 1 0.166667 1 7 0.312500E-01 2 28 0.266454E-14 3 84 0.210942E-14 F(X) = X1^1 * X2^2 * X3^0 * X4^0 * X5^1 0 1 0.166667 1 7 0.312500E-01 2 28 0.266454E-14 3 84 0.310862E-14 F(X) = X1^0 * X2^3 * X3^0 * X4^0 * X5^1 0 1 0.611111 1 7 0.937500E-01 2 28 0.444089E-15 3 84 0.199840E-14 F(X) = X1^2 * X2^0 * X3^1 * X4^0 * X5^1 0 1 0.166667 1 7 0.312500E-01 2 28 0.266454E-14 3 84 0.244249E-14 F(X) = X1^1 * X2^1 * X3^1 * X4^0 * X5^1 0 1 1.33333 1 7 0.625000E-01 2 28 0.177636E-14 3 84 0.888178E-15 F(X) = X1^0 * X2^2 * X3^1 * X4^0 * X5^1 0 1 0.166667 1 7 0.312500E-01 2 28 0.222045E-14 3 84 0.210942E-14 F(X) = X1^1 * X2^0 * X3^2 * X4^0 * X5^1 0 1 0.166667 1 7 0.312500E-01 2 28 0.222045E-14 3 84 0.399680E-14 F(X) = X1^0 * X2^1 * X3^2 * X4^0 * X5^1 0 1 0.166667 1 7 0.312500E-01 2 28 0.222045E-14 3 84 0.399680E-14 F(X) = X1^0 * X2^0 * X3^3 * X4^0 * X5^1 0 1 0.611111 1 7 0.937500E-01 2 28 0.444089E-15 3 84 0.344169E-14 F(X) = X1^2 * X2^0 * X3^0 * X4^1 * X5^1 0 1 0.166667 1 7 0.312500E-01 2 28 0.155431E-14 3 84 0.222045E-15 F(X) = X1^1 * X2^1 * X3^0 * X4^1 * X5^1 0 1 1.33333 1 7 0.625000E-01 2 28 0.177636E-14 3 84 0.222045E-15 F(X) = X1^0 * X2^2 * X3^0 * X4^1 * X5^1 0 1 0.166667 1 7 0.312500E-01 2 28 0.222045E-14 3 84 0.777156E-15 F(X) = X1^1 * X2^0 * X3^1 * X4^1 * X5^1 0 1 1.33333 1 7 0.625000E-01 2 28 0.444089E-15 3 84 0.133227E-14 F(X) = X1^0 * X2^1 * X3^1 * X4^1 * X5^1 0 1 1.33333 1 7 0.625000E-01 2 28 0.177636E-14 3 84 0.133227E-14 F(X) = X1^0 * X2^0 * X3^2 * X4^1 * X5^1 0 1 0.166667 1 7 0.312500E-01 2 28 0.155431E-14 3 84 0.111022E-14 F(X) = X1^1 * X2^0 * X3^0 * X4^2 * X5^1 0 1 0.166667 1 7 0.312500E-01 2 28 0.155431E-14 3 84 0.188738E-14 F(X) = X1^0 * X2^1 * X3^0 * X4^2 * X5^1 0 1 0.166667 1 7 0.312500E-01 2 28 0.155431E-14 3 84 0.188738E-14 F(X) = X1^0 * X2^0 * X3^1 * X4^2 * X5^1 0 1 0.166667 1 7 0.312500E-01 2 28 0.155431E-14 3 84 0.188738E-14 F(X) = X1^0 * X2^0 * X3^0 * X4^3 * X5^1 0 1 0.611111 1 7 0.937500E-01 2 28 0.00000 3 84 0.388578E-14 F(X) = X1^2 * X2^0 * X3^0 * X4^0 * X5^2 0 1 0.416667 1 7 0.203125 2 28 0.666134E-15 3 84 0.111022E-15 F(X) = X1^1 * X2^1 * X3^0 * X4^0 * X5^2 0 1 0.166667 1 7 0.312500E-01 2 28 0.155431E-14 3 84 0.888178E-15 F(X) = X1^0 * X2^2 * X3^0 * X4^0 * X5^2 0 1 0.416667 1 7 0.203125 2 28 0.666134E-15 3 84 0.199840E-14 F(X) = X1^1 * X2^0 * X3^1 * X4^0 * X5^2 0 1 0.166667 1 7 0.312500E-01 2 28 0.888178E-15 3 84 0.188738E-14 F(X) = X1^0 * X2^1 * X3^1 * X4^0 * X5^2 0 1 0.166667 1 7 0.312500E-01 2 28 0.888178E-15 3 84 0.111022E-14 F(X) = X1^0 * X2^0 * X3^2 * X4^0 * X5^2 0 1 0.416667 1 7 0.203125 2 28 0.666134E-15 3 84 0.177636E-14 F(X) = X1^1 * X2^0 * X3^0 * X4^1 * X5^2 0 1 0.166667 1 7 0.312500E-01 2 28 0.222045E-15 3 84 0.188738E-14 F(X) = X1^0 * X2^1 * X3^0 * X4^1 * X5^2 0 1 0.166667 1 7 0.312500E-01 2 28 0.222045E-15 3 84 0.188738E-14 F(X) = X1^0 * X2^0 * X3^1 * X4^1 * X5^2 0 1 0.166667 1 7 0.312500E-01 2 28 0.222045E-15 3 84 0.244249E-14 F(X) = X1^0 * X2^0 * X3^0 * X4^2 * X5^2 0 1 0.416667 1 7 0.203125 2 28 0.666134E-15 3 84 0.111022E-14 F(X) = X1^1 * X2^0 * X3^0 * X4^0 * X5^3 0 1 0.611111 1 7 0.937500E-01 2 28 0.00000 3 84 0.777156E-15 F(X) = X1^0 * X2^1 * X3^0 * X4^0 * X5^3 0 1 0.611111 1 7 0.937500E-01 2 28 0.00000 3 84 0.122125E-14 F(X) = X1^0 * X2^0 * X3^1 * X4^0 * X5^3 0 1 0.611111 1 7 0.937500E-01 2 28 0.00000 3 84 0.00000 F(X) = X1^0 * X2^0 * X3^0 * X4^1 * X5^3 0 1 0.611111 1 7 0.937500E-01 2 28 0.444089E-15 3 84 0.166533E-14 F(X) = X1^0 * X2^0 * X3^0 * X4^0 * X5^4 0 1 0.902778 1 7 0.117188 2 28 0.444089E-15 3 84 0.577316E-14 GM_RULE_PRB Normal end of execution. 23 August 2012 10:18:52.914 AM