Integrals which are individually singular, but which may be combined to yield convergent expressions, are needed for computations of relativistic effects and various properties of atomic and quasiatomic systems. As computations become more detailed and precise, more such integrals are required. This paper presents general formulas for the radial parts of the singular and nonsingular (regular) integrals that occur when three-body systems are described using wave functions that include exponentials in all three interparticle coordinates. Our results are compared with those found in the literature for some of the integrals, and are also shown to be consistent with previously reported results for Hylleraas functions (a limiting case in which one of the exponential parameters is set to zero). (C) 2004 American Institute of Physics.