Abstract: A method of fitting polynomials to the intensity distribution as a function of optical depth is described for use with the Gauss-Seidel iterative technique for solving the equation of radiative transfer. This technique ...