Spectra of Linear Operators
Floquet-Fourier-Hill method
Computational implementation of the Floquet-Fourier-Hill method described in this paper to compute the eigenvalue spectra of linear operators.
You can find a more detailed description, examples, and code here.