function [e]=countE(x,x1,x2,n);

%分別計算樣本x1中第x個樣本劃入樣本x2時的E值,樣本總數為n
%e為所求E值

xm1=(mean(x1))';        
xm2=(mean(x2))'; 
[n1,nc1]=size(x1);
n2=n-n1;
x11=xm1+(xm1-(x1(x,:))')/(n1-1);
x22=xm2-(xm2-(x1(x,:))')/(n2+1);
n11=n1-1;
n22=n2+1;
e=n11*n22*(x11-x22)'*(x11-x22)/n;