test_gaulegf3D.c running 
calling gaulegf(-1.000000, 2.000000, xx, wx, n) 
n=1, equivalent h=3.000000 
xx[1]=0.500000, wx[1]=3.000000 

calling gaulegf(-2.000000, 1.500000, yy, wy, n) 
n=1, equivalent h=3.500000 
yy[1]=-0.250000, wy[1]=3.500000 

calling gaulegf(-1.500000, 0.500000, zz, wz, n) 
n=1, equivalent h=2.000000 
zz[1]=-0.500000, wz[1]=2.000000 

calling gaulegf(-1.000000, 2.000000, xx, wx, n) 
n=2, equivalent h=1.500000 
xx[1]=-0.366025, wx[1]=1.500000 
xx[2]=1.366025, wx[2]=1.500000 

calling gaulegf(-2.000000, 1.500000, yy, wy, n) 
n=2, equivalent h=1.750000 
yy[1]=-1.260363, wy[1]=1.750000 
yy[2]=0.760363, wy[2]=1.750000 

calling gaulegf(-1.500000, 0.500000, zz, wz, n) 
n=2, equivalent h=1.000000 
zz[1]=-1.077350, wz[1]=1.000000 
zz[2]=0.077350, wz[2]=1.000000 

calling gaulegf(-1.000000, 2.000000, xx, wx, n) 
n=3, equivalent h=1.000000 
xx[1]=-0.661895, wx[1]=0.833333 
xx[2]=0.500000, wx[2]=1.333333 
xx[3]=1.661895, wx[3]=0.833333 

calling gaulegf(-2.000000, 1.500000, yy, wy, n) 
n=3, equivalent h=1.166667 
yy[1]=-1.605544, wy[1]=0.972222 
yy[2]=-0.250000, wy[2]=1.555556 
yy[3]=1.105544, wy[3]=0.972222 

calling gaulegf(-1.500000, 0.500000, zz, wz, n) 
n=3, equivalent h=0.666667 
zz[1]=-1.274597, wz[1]=0.555556 
zz[2]=-0.500000, wz[2]=0.888889 
zz[3]=0.274597, wz[3]=0.555556 

calling gaulegf(-1.000000, 2.000000, xx, wx, n) 
n=4, equivalent h=0.750000 
xx[1]=-0.791704, wx[1]=0.521782 
xx[2]=-0.009972, wx[2]=0.978218 
xx[3]=1.009972, wx[3]=0.978218 
xx[4]=1.791704, wx[4]=0.521782 

calling gaulegf(-2.000000, 1.500000, yy, wy, n) 
n=4, equivalent h=0.875000 
yy[1]=-1.756989, wy[1]=0.608746 
yy[2]=-0.844967, wy[2]=1.141254 
yy[3]=0.344967, wy[3]=1.141254 
yy[4]=1.256989, wy[4]=0.608746 

calling gaulegf(-1.500000, 0.500000, zz, wz, n) 
n=4, equivalent h=0.500000 
zz[1]=-1.361136, wz[1]=0.347855 
zz[2]=-0.839981, wz[2]=0.652145 
zz[3]=-0.160019, wz[3]=0.652145 
zz[4]=0.361136, wz[4]=0.347855 

calling gaulegf(-1.000000, 2.000000, xx, wx, n) 
n=5, equivalent h=0.600000 
xx[1]=-0.859270, wx[1]=0.355390 
xx[2]=-0.307704, wx[2]=0.717943 
xx[3]=0.500000, wx[3]=0.853333 
xx[4]=1.307704, wx[4]=0.717943 
xx[5]=1.859270, wx[5]=0.355390 

calling gaulegf(-2.000000, 1.500000, yy, wy, n) 
n=5, equivalent h=0.700000 
yy[1]=-1.835815, wy[1]=0.414622 
yy[2]=-1.192321, wy[2]=0.837600 
yy[3]=-0.250000, wy[3]=0.995556 
yy[4]=0.692321, wy[4]=0.837600 
yy[5]=1.335815, wy[5]=0.414622 

calling gaulegf(-1.500000, 0.500000, zz, wz, n) 
n=5, equivalent h=0.400000 
zz[1]=-1.406180, wz[1]=0.236927 
zz[2]=-1.038469, wz[2]=0.478629 
zz[3]=-0.500000, wz[3]=0.568889 
zz[4]=0.038469, wz[4]=0.478629 
zz[5]=0.406180, wz[5]=0.236927 

f13(x,y,z)=x+2.0*y+3.0 
n=1, n^3=1, area=-23.100000 
n=2, n^3=8, area=-23.100000 
n=3, n^3=27, area=-23.100000 
n=4, n^3=64, area=-23.100000 
n=5, n^3=125, area=-23.100000 

f23(x,y,z)=x*x+2.0*y*y+3.0*x*y+4.0*x+5.0*y+6.0; 
n=1, n^3=1, area=149.625000 
n=2, n^3=8, area=229.250000 
n=3, n^3=27, area=229.250000 
n=4, n^3=64, area=229.250000 
n=5, n^3=125, area=229.250000 

f33(x,y,z)=x*x*x+2.0*y*y*y+3.0*x*x*y+4.0*x*y*y+ 
           5.0*x*x+6.0*y*y+7.0*x*y+8.0*x+9.0*y+10.0 
n=1, n^3=1, area=357.656250 
n=2, n^3=8, area=735.875000 
n=3, n^3=27, area=735.875000 
n=4, n^3=64, area=735.875000 
n=5, n^3=125, area=735.875000 
f43(x,y,z)=x*x*x*x+2.0*y*y*y*y+3.0*x*x*x*y+4.0*x*x*y*y+ 
           5.0*x*x*x*y+6.0*x*x*x+7.0*y*y*y+8.0*x*x*y+ 
           9.0*x*y*y+10.0*x*x+11.0*y*y+12.0*x*y+13.0*x+ 
           14.0*y+15.0 
n=1, n^3=1, area=384.070312 
n=2, n^3=8, area=1353.916667 
n=3, n^3=27, area=1403.981250 
n=4, n^3=64, area=1403.981250 
n=5, n^3=125, area=1403.981250 

fe3(x,y,z)=exp(x+y+z) 
n=1, n^3=1, area=16.354816 
n=2, n^3=8, area=41.685504 
n=3, n^3=27, area=43.463700 
n=4, n^3=64, area=43.503579 
n=5, n^3=125, area=43.504073 

phi(x,2,3,xg)*phi(y)*phi(z) 
n=1, n^3=1,  area=3.737549 
n=2, n^3=8,  area=11.074219 
n=3, n^3=27,  area=3.946852 
n=4, n^3=64,  area=7.420179 
n=5, n^3=125,  area=6.834375 
n=6, n^3=216,  area=6.834375 
n=7, n^3=343,  area=6.834375 

phi(x,2,3,xg)*phi(y)*phi(z)*f13(x,y,z) 
n=1, n^3=1, area=-4.111304 
n=2, n^3=8, area=19.490625 
n=3, n^3=27, area=5.754164 
n=4, n^3=64, area=7.821562 
n=5, n^3=125, area=7.821562 
n=6, n^3=216, area=7.821562 
n=7, n^3=343, area=7.821562 
n=8, n^3=512, area=7.821562 
n=9, n^3=729, area=7.821562 

phi(x,2,3,xg)*phi(y)*phi(z)*f23(x,y,z) 
n=1, n^3=1, area=26.630035 
n=2, n^3=8, area=172.942383 
n=3, n^3=27, area=68.528373 
n=4, n^3=64, area=80.048973 
n=5, n^3=125, area=82.506094 
n=6, n^3=216, area=82.506094 
n=7, n^3=343, area=82.506094 
n=8, n^3=512, area=82.506094 
n=9, n^3=729, area=82.506094 

phi(x,2,3,xg)*phi(y)*phi(z)*f33(x,y,z) 
n=1, n^3=1, area=63.655128 
n=2, n^3=8, area=579.227783 
n=3, n^3=27, area=246.519584 
n=4, n^3=64, area=249.232028 
n=5, n^3=125, area=268.466183 
n=6, n^3=216, area=268.466183 
n=7, n^3=343, area=268.466183 
n=8, n^3=512, area=268.466183 
n=9, n^3=729, area=268.466183 

phi(x,2,3,xg)*phi(y)*phi(z)*f43(x,y) 
n=1, n^3=1, area=68.356264 
n=2, n^3=8, area=870.272095 
n=3, n^3=27, area=513.642158 
n=4, n^3=64, area=423.687276 
n=5, n^3=125, area=469.739544 
n=6, n^3=216, area=477.269345 
n=7, n^3=343, area=477.269345 
n=8, n^3=512, area=477.269345 
n=9, n^3=729, area=477.269345 

phi(x,2,3,xg)*phi(y)*phi(z)*fe3(x,y,z) 
n=1, n^3=1, area=2.910806 
n=2, n^3=8, area=67.231203 
n=3, n^3=27, area=29.191185 
n=4, n^3=64, area=26.958463 
n=5, n^3=125, area=33.398643 
n=6, n^3=216, area=33.151901 
n=7, n^3=343, area=32.785133 
n=8, n^3=512, area=32.738647 
n=9, n^3=729, area=32.735934 
n=10, n^3=1000, area=32.735841 
n=11, n^3=1331, area=32.735839 
n=12, n^3=1728, area=32.735839 
n=13, n^3=2197, area=32.735839 

phi(x,3,5,xg)*phi(y)*phi(z) 
n=1, n^3=1,  area=4.224479 
n=2, n^3=8,  area=0.013227 
n=3, n^3=27,  area=-0.169675 
n=4, n^3=64,  area=4.056794 
n=5, n^3=125,  area=0.569744 
n=6, n^3=216,  area=4.820546 
n=7, n^3=343,  area=2.094790 

phi(x,3,5,xg)*phi(y)*phi(z)*f13(x,y,z) 
n=1, n^3=1, area=-4.646927 
n=2, n^3=8, area=0.026997 
n=3, n^3=27, area=-6.062543 
n=4, n^3=64, area=-1.341706 
n=5, n^3=125, area=-1.069044 
n=6, n^3=216, area=-1.205567 
n=7, n^3=343, area=-1.205567 
n=8, n^3=512, area=-1.205567 
n=9, n^3=729, area=-1.205567 

phi(x,3,5,xg)*phi(y)*phi(z)*f23(x,y,z) 
n=1, n^3=1, area=30.099414 
n=2, n^3=8, area=0.203978 
n=3, n^3=27, area=0.306898 
n=4, n^3=64, area=23.882732 
n=5, n^3=125, area=10.283017 
n=6, n^3=216, area=11.630855 
n=7, n^3=343, area=13.312732 
n=8, n^3=512, area=13.312732 
n=9, n^3=729, area=13.312732 

phi(x,3,5,xg)*phi(y)*phi(z)*f33(x,y,z) 
n=1, n^3=1, area=71.948161 
n=2, n^3=8, area=0.686301 
n=3, n^3=27, area=-7.071318 
n=4, n^3=64, area=77.464340 
n=5, n^3=125, area=34.107383 
n=6, n^3=216, area=37.089313 
n=7, n^3=343, area=43.346226 
n=8, n^3=512, area=43.346226 
n=9, n^3=729, area=43.346226 

phi(x,3,5,xg)*phi(y)*phi(z)*f43(x,y) 
n=1, n^3=1, area=77.261764 
n=2, n^3=8, area=0.813808 
n=3, n^3=27, area=251.181116 
n=4, n^3=64, area=342.028438 
n=5, n^3=125, area=168.004909 
n=6, n^3=216, area=147.212774 
n=7, n^3=343, area=178.204264 
n=8, n^3=512, area=183.332968 
n=9, n^3=729, area=183.332968 

phi(x,3,5,xg)*phi(y)*phi(z)*fe3(x,y,z) 
n=1, n^3=1, area=3.290028 
n=2, n^3=8, area=0.088732 
n=3, n^3=27, area=-27.457147 
n=4, n^3=64, area=-17.685769 
n=5, n^3=125, area=-15.897891 
n=6, n^3=216, area=1.689445 
n=7, n^3=343, area=-3.896706 
n=8, n^3=512, area=-10.090514 
n=9, n^3=729, area=-8.902912 
n=10, n^3=1000, area=-8.266815 
n=11, n^3=1331, area=-8.190459 
n=12, n^3=1728, area=-8.186035 
n=13, n^3=2197, area=-8.185883