We developed an efficient Monte-Carlo algorithm to solve dust-scattering radiative transfer problems for continuum radiation. The method calculates the scattered intensities for various anisotropic factors ($g_i$) all at once, while actual photon pack...
We developed an efficient Monte-Carlo algorithm to solve dust-scattering radiative transfer problems for continuum radiation. The method calculates the scattered intensities for various anisotropic factors ($g_i$) all at once, while actual photon packets are tracked following a scattering phase function given by a single anisotropic factor ($g_0$). The algorithm was tested by applying the method to a dust cloud embedding a star at the cloud center and found to provide accurate results within the statistical fluctuation that is intrinsic in Monte-Carlo simulations. It was found that adopting $g_0$ = 0.4 - 0.5 in the algorithm is most efficient. The method would be efficient in estimating the anisotropic factor of the interstellar dust by comparing the observed data with radiative transfer models.