function?centdist?=?alldist（centers）
%?output:?matrix?of?all?pairwise?distances
%?input:?data?points?（centers）

k?=?size（centers1）;
centdist?=?zeros（kk）;
for?j?=?1:k
????centdist（1:j-1j）?=?calcdist（centers（1:j-1:）centers（j:））;
end
centdist?=?centdist+centdist‘;