Abstract
The integral integral f(r) mod r-R mod -1 mod r-R' mod -1 dr is reduced to a form suitable for numerical computation. The reduction does not depend on the form of the f(r) and consequently the method works equally well for Gaussian and Slater orbitals.