Abstract
A novel method is described for avoiding divergent contributions and delta function contributions in calculations involving singular potentials. The method was announced in a Preliminary Communication (1983, Molec. Phys., 48, 1131); this paper presents the details of the method and extends it to angular problems. The method is illustrated by calculations of relativistic corrections in the hydrogen atom using the Foldy-Wouthuysen hamiltonian; the corrections are correctly calculated to order mc 2 α8 for the 1s 1/2, 2P 1/2 and 2p 3/2 states and also to order mc 2 α6 for the 2s 1/2 state.