  The method proposed by Kinoshita and Nio to compute higher order vacuum
polarization contributions to the Coulomb potential in muonic hydrogen is
generalized to obtain relativistic corrections to their results.
