Solution of radiative transfer equation with spherical-symmetry in partially scattering medium

We have solved the equation of radiative transfer in spherical symmetry with scattering and absorbing medium. We have set the albedo for single scattering to be equal to 0.5. We have set the Planck function constant throughout the medium in one case and in another case the Planck function has been set to vary asr ^–2. The geometrical extension of the spherical shell has been taken as large as one stellar radius. Two kinds of variations of the optical depth are employed (1) that remains constant with radius and (2) that varies asr ^–2. In all these cases the internal source vectors and specific intensities change depending upon the type of physics we have employed in each case.