  A mechanism is described to symmetrize the ultraspherical spectral method for
self-adjoint problems. The resulting discretizations are symmetric and banded.
An algorithm is presented for an adaptive spectral decomposition of
self-adjoint operators. Several applications are explored to demonstrate the
properties of the symmetrizer and the adaptive spectral decomposition.
