假设x是n个点的坐标

　　x=[111;

　　222;

　　333

　　...;

　　nnn

　　];

　　fori=1:n

　　forj=1:n

　　dis(i,j)=sqrt(sum((x(i,:)-y(j,:)).^2));

　　end;

　　end;