普氏分析法的计算过程，使用的误差，matlab

syms a b c d e f A B C D E F
X=[a,b;c,d;e,f];
Y=[A,B;C,D;E,F];
% X=[1,2;3,3;4,5]
% Y=[1,2;3,4;4,2]
muX = mean(X,1)
muY = mean(Y,1)
X0 = X - repmat(muX, 3, 1)
Y0 = Y - repmat(muY, 3, 1)

ssqX = sum(X0.^2,1)
ssqY = sum(Y0.^2,1)

ssqX = sum(ssqX)
ssqY = sum(ssqY)

normX = sqrt(ssqX)
normY = sqrt(ssqY)

muX =

[ a/3 + c/3 + e/3, b/3 + d/3 + f/3]

muY =

[ A/3 + C/3 + E/3, B/3 + D/3 + F/3]

X0 =

[ (2a)/3 - c/3 - e/3, (2b)/3 - d/3 - f/3]
[ (2c)/3 - a/3 - e/3, (2d)/3 - b/3 - f/3]
[ (2e)/3 - c/3 - a/3, (2f)/3 - d/3 - b/3]

Y0 =

[ (2A)/3 - C/3 - E/3, (2B)/3 - D/3 - F/3]
[ (2C)/3 - A/3 - E/3, (2D)/3 - B/3 - F/3]
[ (2E)/3 - C/3 - A/3, (2F)/3 - D/3 - B/3]

ssqX =

[ (a/3 + c/3 - (2e)/3)^2 + (a/3 - (2c)/3 + e/3)^2 + (c/3 - (2a)/3 + e/3)^2, (b/3 + d/3 - (2f)/3)^2 + (b/3 - (2d)/3 + f/3)^2 + (d/3 - (2b)/3 + f/3)^2]

ssqY =

[ (A/3 + C/3 - (2E)/3)^2 + (A/3 - (2C)/3 + E/3)^2 + (C/3 - (2A)/3 + E/3)^2, (B/3 + D/3 - (2F)/3)^2 + (B/3 - (2D)/3 + F/3)^2 + (D/3 - (2B)/3 + F/3)^2]

ssqX =

(a/3 + c/3 - (2e)/3)^2 + (a/3 - (2c)/3 + e/3)^2 + (c/3 - (2a)/3 + e/3)^2 + (b/3 + d/3 - (2f)/3)^2 + (b/3 - (2d)/3 + f/3)^2 + (d/3 - (2b)/3 + f/3)^2

ssqY =

(A/3 + C/3 - (2E)/3)^2 + (A/3 - (2C)/3 + E/3)^2 + (C/3 - (2A)/3 + E/3)^2 + (B/3 + D/3 - (2F)/3)^2 + (B/3 - (2D)/3 + F/3)^2 + (D/3 - (2B)/3 + F/3)^2

normX =

((a/3 + c/3 - (2e)/3)^2 + (a/3 - (2c)/3 + e/3)^2 + (c/3 - (2a)/3 + e/3)^2 + (b/3 + d/3 - (2f)/3)^2 + (b/3 - (2d)/3 + f/3)^2 + (d/3 - (2b)/3 + f/3)²⁾(1/2)

normY =

((A/3 + C/3 - (2E)/3)^2 + (A/3 - (2C)/3 + E/3)^2 + (C/3 - (2A)/3 + E/3)^2 + (B/3 + D/3 - (2F)/3)^2 + (B/3 - (2D)/3 + F/3)^2 + (D/3 - (2B)/3 + F/3)²⁾(1/2)