A numerical program, which combines Monte Carlo method with geometrical ray trace used to study the behaviour of multi photon excitation (MPE) fluorescence microscopic imaging through turbid media, is presented. The mathematical model used to microscopically image multi layered tissues is presented. The basic skeleton of the simulation program is demonstrated with flow chart. A skill is introduced to improve the traditional Monte Carlo method, which enhances the capability of simulation of the converging light. The prominent effect on microscopic imaging caused by scattering characteristics of turbid media is shown through preliminary simulation experiments.