Using quantum physics, you can determine the f eigenvalues and matching eigenvectors for systems in which the energies are degenerate. Take a look at this unperturbed Hamiltonian:

In other words, several states have the same energy. Say the energy states are f-fold degenerate, like this:

How does this affect the perturbation picture? The complete Hamiltonian, H, is made up of the original, unperturbed Hamiltonian, H0, and the perturbation Hamiltonian,

In zeroth-order approximation, you can write the eigenfunction

as a combination of the degenerate states

Note that in what follows, you assume that

if m is not equal to n. Also, you assume that the

are normalized — that is,

Plugging this zeroth-order equation into the complete Hamiltonian equation, you get

Now multiplying that equation by

gives you

Using the fact that

if m is not equal to n gives you

Physicists often write that equation as

where

And people also write that equation as

where E(1)n = En – E(0)n. That's a system of linear equations, and the solution exists only when the determinant to this array is nonvanishing:

The determinant of this array is an fth degree equation in E(1)n, and it has f different roots,

Those f different roots are the first-order corrections to the Hamiltonian. Usually, those roots are different because of the applied perturbation. In other words, the perturbation typically gets rid of the degeneracy.
So here's the way you find the eigenvalues to the first order — you set up an f-by-f matrix of the perturbation Hamiltonian,

Then diagonalize this matrix and determine the f eigenvalues

and the matching eigenvectors:

Then you get the energy eigenvalues to first order this way:

And the eigenvectors are
