double Sxx = 0, Sxy = 0, Sxz = 0, Syx = 0, Syy = 0, Syz = 0, Szx = 0, Szy = 0, Szz = 0;
for (int i = n + 1; --i >= 1;) {
Point3f aij = centerAndPoints[0][i];
Point3f bij = centerAndPoints[1][i];
Point3f ptA = new Point3f(aij);
ptA.sub(centerAndPoints[0][0]);
Point3f ptB = new Point3f(bij);
ptB.sub(centerAndPoints[0][1]);
Sxx += (double) ptA.x * (double) ptB.x;
Sxy += (double) ptA.x * (double) ptB.y;
Sxz += (double) ptA.x * (double) ptB.z;