Solution (source code)

= Solution

For a normalized eigenstate $|n\rangle$ with a <nondegenerate energy eigenvalue>, <first-order nondegenerate perturbation theory> and <second-order nondegenerate perturbation theory> give
$$
\boxed{
E_n(H)=E_n+\langle n|V|n\rangle
+\sum_{m\ne n}
\frac{|\langle m|V|n\rangle|^2}{E_n-E_m}
+O(V^3).}
$$