clc,clear all
X=[];
[row,col]=size(X);
Xchu=[];
for k=1:col
    tmp=X(:,k)./X(:,1);
    Xchu=[Xchu,tmp];
end
Xchu
delta=[];
for k=2:row
    tmp=abs(Xchu(k,:)-Xchu(1,:));
    delta=[delta;tmp];
end
M=max(max(delta));
m=min(min(delta));
r=[];
kk=0.5;
for i=1:row-1
    r1=[];
    for j=1:col
        tmp=(m+kk*M)/(delta(i,j)+kk*M);
        r1=[r1,tmp];
    end
    r=[r;r1];
end
r
R=mean(r,2)