clear?all
x=load（‘are.dat‘）
n=length（x）
sd=std（x）
r=0.2*sd
for?ii=1:2
m=ii+1;
num=zeros（n-m+11）;
for?i=1:n-m+1
for?j=1:n-m+1
if?j~=i
for?k=1:m
d（k）=abs（x（i+k-1）-x（j+