One of the quantities usually required when solving the equation of radiative transfer is the intensity of radiation emerging from the surface of the medium under consideration. For multi-dimensional situations however, the methods presented to date have been numerical, and these first calculate the so-called source function Sv (r, Ω) as a function of position r, angle Ω and frequency v. This is generally the most difficult part of the exercise since an integro-difierential equation must be solved. The emergent intensity is then determined by solving a relatively simple first order differential equation by any of the well known numerical integration schemes. However, if the emergent intensity is required at a large number of angles, frequencies, and positions on the surface of the medium, and this is usually the case, the amount of computing needed may be considerable.