 Solution

If $e_1,\dots ,e_n$ are orthonormal, any relation ${\sum }_i{c}_i{e}_i=0$ gives $c_j=0$ on taking the dot product with $e_j$. Thus the vectors are linearly independent. A linearly independent set of $n$ vectors in the $n$-dimensional vector space ${\mathbb{R}}^n$ is a basis, proving the claim.

Expand both $x$ and $f$ in this orthonormal basis, writing $f_i=e_i\cdot f$. Since $A{e}_i={\lambda }_i{e}_i$, the linear equation becomesComparing the unique basis coefficients gives, when $\mu$ is not an eigenvalue,If $\mu ={\lambda }_1$, the equation is solvable exactly when $e_i\cdot f=0$ for every $i$ with ${\lambda }_i=\mu$. For those indices $a_i$ is arbitrary; for all other indices the displayed formula still applies. Thus repeated eigenvalues require orthogonality to the entire corresponding eigenspace, not merely to one chosen eigenvector. This is the finite-dimensional Fredholm solvability condition for a self-adjoint operator.

For the specific real symmetric matrix, chooseTheir dot products with $f$ are $-1/\sqrt{2}$, $3/\sqrt{2}$ and $3$. At $\mu =2$ this givesOne can check $(A-2I)x=f$ directly. At $\mu =1$, the first coefficient equation would require $0\cdot a_1=-1/\sqrt{2}$, which is impossible. Hence there is no solution when $\mu =1$; the nonzero projection onto its eigenspace is the obstruction.