Let $S$ be the unilateral shift operator on $\mathcal l^2$ (which shifts one place to the right) and $S^*$ its adjoint, the backward shift (which shifts one place to the left). I've been trying to find the spectrum of $T=S+S^*$.
Since $T$ is self-adjoint, any eigenvalues would be real, and I've shown that no $|\lambda|\ge 2$ can be an eigenvalue. The case $|\lambda|< 2$ corresponds to the roots of $ t^2=\lambda t-1 $ being complex (this is the characteristic polynomial corresponding to the recurrence relation $a_n=\lambda a_{n-1}+a_{n-2}$). These complex roots have absolute value 1, and this implies that $a_n$ does not converge to zero.
So I've shown that $T$ has no eigenvalues. But I'm still looking for an elementary way to find that the spectrum is $[-2,2]$, as answered by Joel below. Although I appreciate Joel's answer below and I'm sure it'll be valuable for many in the community, unfortunately I only know basic Hilbert space theory. Any hints would be very helpful!
You know that $\|S\| \le 1$ which gives $\|S^{\star}\|=\|S\| \le 1$ and, hence, $\|S+S^{\star}\| \le 2$. However $S+S^{\star}$ is selfadjoint and, so, its norm and spectral radius are the same, and its spectrum is real. That means $\sigma(S+S^{\star})\subseteq[-2,2]$. So that part is fairly straightforward.
Suppose $\lambda \in [-2,2]$. You want to show that $\lambda$ is definitely in the spectrum. So, try to solve $$ (S+S^{\star}-\lambda I)x = y $$ and see if something goes wrong for some value of $y=\{ y_0,y_1,\cdots\}$. The algebra is a mess, unless you use a power series. $$ x(z) = \sum_{n=0}^{\infty}x_{n}z^{n},\;\;\; y(z)=\sum_{n=0}^{\infty}y_{n}z^{n}. $$ These power series converge in $|z| < 1$ because the coefficients are square integrable. What's nice, though, is that $$ (Sx)(z) = zx(z),\;\;\;\; (S^{\star}x)(z) = \frac{x(z)-x(0)}{z}. $$ Now, the equation you want to solve is $$ zx(z)+\frac{x(z)-x(0)}{z}-\lambda x(z)=y(z). $$ So you don't have to use full-blown function theory, but the power series is really helpful. Here you've got $$ (z+\frac{1}{z}-\lambda)x(z)=y(z)+\frac{x(0)}{z} \\ (z^{2}+1-\lambda z)x(z) = zy(z)+x(0). $$ The roots of $z^{2}-\lambda z+1$ are $$ \frac{\lambda}{2}\pm i\sqrt{1-\frac{\lambda^{2}}{4}} $$ For $-2 \le \lambda \le 2$, the modulus of the above is $$ \frac{\lambda^{2}}{4}+\left(1-\frac{\lambda^{2}}{4}\right)=1. $$ Here's where it's best to take a specific $y$ and show that you can't get a solution. The obvious choice is $y(z)\equiv 1$. Then the power series equation gives you $$ x(z) = \frac{z+x(0)}{z^{2}-\lambda z+1}. $$ You've got two distinct roots of the denominator on the unit circle for $-2 < \lambda < 2$, and there's no way for the numerator to dampen the effect of both roots. So $|x(re^{i\theta})| \sim \frac{1}{1-r}$ if $e^{i\theta}$ is an undamped root of the denominator. But look at the growth rate of $x(z)$ just by using the Cauchy-Schwarz inequality: \begin{align} |x(z)| & \le \sum_{n=0}^{\infty}|x_n||z|^{n} \\ & \le \left(\sum_{n=0}^{\infty}|x_n|^{2}\right)^{1/2}\left(\sum_{n=0}^{\infty}|z|^{2n}\right)^{1/2} \\ & = \|x\|\frac{1}{\sqrt{1-|z|^{2}}} \\ & \le \|x\|\frac{1}{\sqrt{1-|z|}}. \end{align} That's a contradiction for any $-2 < \lambda < 2$. So $(-2,2)\subseteq\sigma(S+S^{\star})$. Because spectrum is closed, then $[-2,2]\subseteq \sigma(S+S^{\star})$.