%例8.3.1多總體的均值和方差估計值距離計算

G1=[10   9   3
    8   7   1
    8   4   4
    2   8   0];

G2=[6   2   7
    1   4   5
    2   5   4
    3   1   4];

G3=[4   1   10
    1   2   6
    1   0   11];

[mu1,sigma1,S1]=musig(G1);
[mu2,sigma2,S2]=musig(G2);
[mu3,sigma3,S3]=musig(G3);
[S,S_1]=sigxd(G1,G2,G3);
S1,S2,S3,S,S_1