The propagation of polarized light in anisotropic medium like some fibrous tissues (e.g., muscle, skin, tooth) is investigated with Monte Carlo method. In the simulation, the myofibrils and collagen fibers are approximated as long cylinders and the tissue phantom is composed of spherical and cylindrical structures. We present the realization and validation of Monte Carlo program, and then apply the program to simulate and analyze polarization imaging process of muscle. The good agreement between the simulation results and the experimental results validates the assumption of the phantom composition and shows the feasibility of Monte Carlo method. This paper also presents how to describe tissue anisotropy according to two parameters derived from the polarization imaging.