kill
(all);
M: matrix(
[6,
sqrt
(2),
(2),0],
[
(2),5,1,0],
(2),1,5,0],
[0,0,0,4]
);
eigenvalues(M);
eigenvectors(M);