Convergence of Monte Carlo algorithms for reconstructing the scattering phase function with polarization

The problem of reconstructing the atmospheric scattering phase function from groundbased solar almucantar sky brightness observations is considered. A newiterative algorithm for solving this problem was developed as a combination of existing additive and multiplicative methods of refining the single-scattering contribution to the observed brightness with polarization of scattered radiation in the atmosphere. Also, some modifications of these methods were proposed. The objective of this paper is to numerically substantiate the convergence of these methods. For this purpose, an algorithm of Jacobi matrix calculation for the iteration operators of the methods was developed, and calculations were carried out for various parameters of the atmosphere.