#include #include typedef double R; R rnd(void); int main(){int c = 10000; R xx=0, xy=0, xz=0; R yy=0, yz=0, zz=0; R xs=0, ys=0, zs=0; R a = sqrt(3)/2; int j=c; while (j--){ R x = rnd(), b = a*rnd(), y = b+0.5*x, z = b-0.5*x; xx += x*x; xy += x*y; xz += x*z; yy += y*y; yz += y*z; zz += z*z; xs += x; ys += y; zs += z; if(j<50) printf("%6.3f %6.3f %6.3f\n", x, y, z);} printf("%8.5f %8.5f %8.5f\n", xx/c, xy/c, xz/c); printf("%8.5f %8.5f %8.5f\n", xy/c, yy/c, yz/c); printf("%8.5f %8.5f %8.5f\n\n", xz/c, yz/c, zz/c); printf("%8.5f %8.5f %8.5f\n", xs/c, ys/c, zs/c); return 0;}