We present a Monte Carlo simulation program for the propagation of polarized photons in an anisotropic scattering model consisting of poly-dispersed spherical and infinite long cylindrical scatterers. The cylinders are aligned following a Gaussian distribution. Densities and sizes of the spherical and cylindrical scatterers, as well as the orientation of the cylinders are variables for the simulation of different anisotropic media. The good agreement between the simulation and experimental results of polarization imaging confirms the validity of the polarization-dependent Monte Carlo simulation program.