function MKI=klgm

for i=1:n

    for j=1:n

        for k=1:m

            X(k,i)=B(k,i)/sum(B(:,i));

            X(k,j)=B(k,j)/sum(B(:,j));

            e(k)=abs(X(k,i)-X(k,j));

            MKI(i,j)=1-sum(e)*0.5;

        end

    end   

end

disp('改进的克鲁格曼系数矩阵MKI为')