matlab 求点间距离矩阵

clear;
n=input('请输入点的个数：');
x=zeros(1,n);
y=x;
for i=1:n
disp(['请输入第',num2str(i),'个点的坐标']);
x(i)=input('x=?');
y(i)=input('y=?');
end
r=zeros(n,n);
for j=1:n
for k=1:n
r(j,k)=sqrt((x(k)-x(j))^2+(y(k)-y(j))^2);
end
end