Three-dimensional steady-state radiative integral transfer equations (RITEs) for a cubic absorbing and isotropically scattering homogeneous medium are solved using the method of "subtraction of singularity". Surface integrals and volume integrals are carried out analytically to eliminate singularities, to assure highly accurate solutions, and to reduce the computational time. The resulting system of linear equations for the incident energy is solved iteratively. Six benchmark problems for cold participating media subjected to various combinations of externally uniform/non-uniform diffuse radiation loads are considered. The solutions for the incident energy and the net heat flux components are given in tabular form for scattering albedos of omega = 0, 0.5 and 1. (C) 2007 Elsevier Ltd. All rights reserved.