Abstract: A direct method is described for the numerical solution of the spherical harmonies approximation to the equation of transfer for plane-parallel atmospheric models with arbitrary nonhomogeneities in the vertical direction, ...