D = diag([1,2,3;4,5,6;7,8,9]);
disp(D)
diag([1,1,1,2,2,2])