The Spectral Theorem

The following code will calculate a spectral factorization H=UDU^* of a Hermitian matrix H, assuming that the eigenvalues of H are nice enough.