Published online by Cambridge University Press: 25 April 2016
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.