We have developed an efficient Monte Carlo method by which we can evaluate the evolution of comets. There are many poorly known evolutional parameters, and we have investigated the influence of these parameters on the final populations and the inclination distributions of short-period comets. We compare the observed and calculated inclination distributions of different comet populations and obtain a good fit for the inclinations of the Jupiter family comets by assuming a mantle blow-off and a sudden brightening of the comet when its perihelion distance is lowered in a major jump.