A highly convergent perturbative method for the solution of systems of coupled equations arising from the Schrödinger equation