I posted this question in MathOverflow but was not lucky with the answers, so wil try here.
Suppose I have a matrix given by a sum
$$A=D+\epsilon B$$
where $D$ is diagonal and $\epsilon$ is small, and I want the eigenvalues of $A$ as power series in $\epsilon$. The first two orders in perturbation theory are well known. Third and higher orders are briefly discussed here. However, the equations become horrible.
I hear that Feynman diagrams are an efficient way to formulate perturbation theory, but I can't find an accessible exposition of this approach. Note that I have in mind the simple matrix setting. I don't want vacuum states, quantum field theory, path integrals, many body, etc. Can someone help?