“…The change in light intensity at every location is solved considering the incoming light that is the light from the source plus the scattered light coming from other particles, and the outgoing light, namely the scattered light contribution of that particular location and the loss due to absorption. Here, the wavelength-dependent light interaction properties, the spectral volumetric absorption and scattering coefficients, α λ and σ λ , as well as the scattering phase function, p (Ω′ → Ω), have to be imported to resolve the spectral radiation intensity, I λ,Ω ( s , t ) reaching a given point s ( x ) in space and time t , having a given direction of propagation Ω defined by the polar and azimuthal angles, traveling along distances measured by the spatial parameter s . For the rigorous RTE, isotropic scattering and diffuse reflectance phase functions are usually used. − The accuracy of the solutions is dictated by simplifications made on the scattering spatial distribution function . For symmetry arguments, the six-flux and two-flux ,− approximations are utilized the most.…”