Three Monte Carlo programs of polarized light transport into scattering media: part I

Propagation of light into scattering media is a complex problem that can be modeled using statistical methods such as Monte Carlo. Few Monte Carlo programs have so far included the information regarding the status of polarization of light before and after a scattering event. Different approaches have been followed and limited numerical values have been made available to the general public. In this paper, three different ways to build a Monte Carlo program for light propagation with polarization are given. Different groups have used the first two methods; the third method is original. Comparison in between Monte Carlo runs and Adding Doubling program yielded less than 1 % error.