This chapter is devoted to the discussion of a distinctive Monte Carlo method for evaluation of angular distribution of outgoing polarized radiation. The algorithm in consideration is based on the modification of N. N. Chentsov method for unknown probability density evaluation via the orthonormal polynomial expansion. A polarization was introduced into a mathematical model of radiation transfer with use of four-dimensional vector of Stokes parameters. Corresponding weighted Monte Carlo algorithm was constructed. Using this method and precise computer simulation, the angular distribution of outgoing radiation was investigated. Special attention was given to the value of polarization impact in the mathematical model of radiation. Algorithm in consideration allows us precisely estimate even a small effect of polarization as well as a deviation of the calculated angular distribution from the Lambertian one.