Spectroscopy of polarized light scattering is a novel technique prospected to diagnose precancerous lesions. However, a perfect theoretical method to study the mechanism of this technique is absent. With the parallel computing and Stokes-Müller formalism to describe the polarization of photon and media, respectively, a Monte-Carlo program, which can simulate the spectra of polarized light scattering, is developed. The program is designed with a master-slave architecture based on message passing interface (MPI), and a dynamic load balance is obtained via taskpool scheme. It has been validated by the comparison of polarized light backscattering spectra in turbid media between simulation and experiments. Furthermore, results show that this program can provide excellent scalability of task size and slaver number. For long-time required simulations, the computational time of using n slavers can be reduced as almost 1/n times of using 1 slaver.