function distancee=dis(X)
distancee=zeros(3,3);
for i=1:3
        for j=1:3

            for k=1:3
            distancee(i,j)=(X(k,i)-X(k,j))^2+distancee(i,j);
            end
            distancee(i,j)=distancee(i,j)^0.5;
        end
    end