%DUERER  Basic computations on Duerer's magic square.
%
% Thorsten Hansen  2008-04-11
M = magic(4); % Duerer's magic square
M = M(:, [ 1 3 2 4]); % switch columns 2 and 3 to get exactly Duerer's
                      % matrix; this is not really necessary
disp(['Spaltensumme: ' int2str(sum(M))]) % int2str works also on vectors
disp(['Zeilensumme: ' int2str(sum(M'))]) 
disp(['Summe der Hauptdiagonalen: ' int2str(sum(diag(M)))]);

m = mean(sum(M));
sd = std(sum(M));
disp(['Mittelwert der Spaltensumme: ' num2str(m)])
disp(['Standardabweichung der Spaltensumme: ' num2str(sd)])