Numerical simulation for a one-dimensional radiative transfer equation

This numerical simulation has been developed for the diffusion approximation of a one-dimensional stationary radiative transfer equation. Specifically, as the mean free path of the photons tends to zero this simulation depicts the formation of boundary layers and the fact that at the interior of the domain the radiation intensity becomes isotropic solving a Laplace equation. This numerical simulation supports the results obtained analyzing the radiative transfer equation via matched asymptotic expansions.