A numerical method for obtaining solution of radiative transfer equation in spherically symmetric media with spherical harmonic approximation is presented. The angle derivative is approximated by an orthonormal polynomial ...