Calculates the eigenvalues of a special form of Hill's equation that solves it following the Floquet-Fourier-Hill method. The equation is
$\begin{equation} \frac{d^2\phi_{n}}{dy^2} + \left[\lambda_{n} - 2qU(y) \right]\phi_{n} = 0 \end{equation} $
where