Elements of the theory of angular momentum in the Rayleigh-Gans-Debye (RGD) approximation are used to derive an expansion of the phase function in Legendre polynomials for a single spherical particle and polydisperse spherical and randomly oriented spheroidal particles. The numerical calculations are then compared with rigorous results of the T-matrix method.