function d=diagonale(matrix)
%function d=diagonale(matrix)
%gibt die Elemente auf der Diagonalen der quadratischen Matrix "matrix"
%zurueck.

[rows,cols]=size(matrix);  %Groesse der Matrix

letzter_index=rows*cols;   %Index des Elements unten rechts

v=matrix(1:rows+1:letzter_index)  %Rausgreifen der Elemente auf der Diagonalen