TESTS
skew[i] = calloc(TESTS, sizeof(*skew[i]));
for (j = 0; j < TESTS; j++) {
for (j = 0; j < TESTS; j++) {
aveskew[i] /= TESTS;
sumsq /= TESTS;
for (i = 0; i < TESTS; i++) {