A=xlsread（‘d:/2.xlsx‘）;
x=A（:1）;
y=A（:2）;
N=length（y）;
n=length（y）;
Sk=zeros（size（y））;
UFk=zeros（size（y））;
s=0;
for?i=2:n
????for?j=1:i
????????if?y（i）>y（j）
????????????s=s+1;
????????else
????????????s=s+0;
????????end;
????end
????Sk（i）=s;
????E=i*（i-1）/4;
????Var=i*（i-1）*（2*i+5）/72;
????UFk（i）=（Sk（i）-E）/sqrt（Var）;
end;
y2=zeros（size（y））;
Sk2=zeros（size（y））;
UBk=zeros（size（y））;
s=0;
for?i=1:n
????y2（i）=y（n-i+1）;
end;
for?i=2:n
????for?j=1:i
????????if?y2（i）>y2（j）
????????????s=s+1;
????????else
????????????s=s+0;
????????end;
????end;
????Sk2（i）=s;
????E=i*（i-1）/4;
????Var=i*（i-1）*（2*i+5）/72;
????U