Solution (source code)

= Solution

The normalized vectors $u_1=(1,1)^T/\sqrt2$ and $u_2=(-1,1)^T/\sqrt2$ give <eigenvalues> $I_0$ and $3I_0$. Therefore
$$
\boxed{U=\frac1{\sqrt2}\begin{pmatrix}1&-1\\1&1\end{pmatrix},\qquad
D=\begin{pmatrix}I_0&0\\0&3I_0\end{pmatrix}.}
$$
The columns are orthonormal, so $U^{-1}=U^T=U^*$ and direct multiplication verifies $U^{-1}MU=D$. Signs of individual <eigenvectors> are free; the displayed signs also make $\det U=+1$.