Consider the positive definite symmetric matrix:
If is the unique positive definite square root of , then the components of are:
By statement 2 of the notes, has a unique positive definite square root. This positive definite square root is a symmetric matrix that has the same eigenvectors of but whose eigenvalues are the positive square roots of the eigenvalues of M. In that case:
has 8 square roots. The following Mathemtica code can be used to find all possible matrices that are square roots of along with their eigenvalues and eigenvectors:
M = {{16, 0, 0}, {0, 4, 0}, {0, 0, 1}};
A = {{A11, A12, A13}, {A21, A22, A23}, {A31, A32, A33}};
A.A // MatrixForm
unknowns = Flatten[A]
a = Solve[A.A == M, unknowns]
As = Table[A /. a[[i]], {i, 1, 8}]
Do[Print[As[[i]] // MatrixForm], {i, 1, Length[As]}]
Do[Print[Eigensystem[As[[i]]] // MatrixForm], {i, 1, Length[As]}]