Abstract:
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 and this is represented by a matrix, for a given beam of rays. An error analysis of the curvature matrix and results using the solution for few representative cases has been presented.