% a=[1,2,3;
% 2,3,4];
% B=[2,3,4];
% row=size(a,1);
% D=zeros(row,row);
% for i=1:row
% for j=i+1:row
% D(i,j)=((a(i,1)-a(j,1))^2+(a(i,2)-a(j,2))^2+(a(i,3)-a(j,3))^2)^0.5;
% D(j,i)=D(i,j);
% end
% end
% % C=size(A,1);
% D;
% % for i=1:2
% plot3(a(:,1),a(:,2),a(:,3))
% end
A=ceil(rand(11,3)*40);