Approximate method of optical energy distribution calculation in multiple scattering mediums

Approximate approach of transport equation solution for Henyey-Greenstein anisotropic phase function is submitted. Propagating light was presented as a sum of collimated and diffusive components, that allows to lead in fast and slow variables and decay transport equation in two parts. An equation for collimated flux results in analytical solution and describes Bouguer propagation. And integro-differential equation for diffuse part simplifies and transforms into engaging partial equations system. Each rank of the system brings to initial value problem solving by Cauchy method with corresponding ranks derivatives joining boundary conditions. It was shown the fast convergence of proposed approximate solution by its comparison with results of numerical simulations by Monte Carlo method for the same multiple scattering medium.