Solution (source code)

= Solution

The <eigenvector> equations $Q\mathbf e_i=\lambda_i\mathbf e_i$ combine column by column into $QE=ED$. Since $|E|\ne0$, multiplication by the <inverse matrix> gives the <diagonalization>
$$
\boxed{Q=EDE^{-1}}.
$$
The nonzero <determinant> of $E$ ensures that its columns form an <eigenbasis>.