p=500*rand(6,2)  % Position der Staedte
% Vektoren zu Matrizen aufblaehen fuer gegenseitige Differenzen
x=p(:,1)*ones(1,6)
y=repmat(p(:,2),1,6)
% Andere Variante 
% [x,y]=meshgrid(p(:,1),p(:,2));

% Abstaende bestimmen
d=sqrt((x-x').^2+(y-y').^2)

fid=fopen('tabelle.txt','w'); % Datei zum Schreiben oeffnen

% Erste Zeile mit den Ueberschriften
fprintf(fid,'         %d      %d      %d      %d      %d      %d\n',1:6);

% Jeweils Zeile mit den Abstaenden (muss fuer jede Zeile wiedeholt werden)
% fprintf(fid,'  %d   %6.1f %6.1f %6.1f %6.1f %6.1f %6.1f\n',1,d(:,1));

% Einfacher unter Ausnutzung der vektorisierten Ausgabe 
fprintf(fid,'  %d   %6.1f %6.1f %6.1f %6.1f %6.1f %6.1f\n',[[1:6];d]);

fclose(fid); % Datei schliessen

% save -ASCII Variable d   % Speichern der Variable im ASCII-Format