Abstract
Fluxes and intensifies of light scattered by a model atmosphere are computed by a spherical harmonics approximation and by an iterative method of solving the radiative transfer equation and are compared. The large differences in the net fluxes and intensifies reported by Dave and Armstrong (1974) for the two methods are reduced here by making a few changes in the iterative routine. Decreasing the polar angle increment from 2° to 1° in the iterative method of computing the source function does not improve the results as suggested by Dave and Armstrong.
