A Padé extrapolated inverse power method for coupled Schrödinger-like equations applied to the two-body relativistic bound state problem in quantum electrodynamics (Q1339949)

The author extends the inverse power method combined with Padé extrapolated bound state eigenvalues from one Schrödinger-like equation to two and four coupled Schrödinger-like equations. An inversion technique for large banded and blocked matrices is presented analogous to that in textbooks for unblocked and banded matrices. It is shown that this method, in conjunction with logarithmic scaling and Padé extrapolation techniques, can be used to obtain numerical solutions for the positronium and muonium spectrum that agree with perturbation theory through order.