A mathematical model was developed to represent pyrolysis. The components of primary and secondary pyrolysis reactions were simply lump into different groups and were represented through a set of pseudo-first-order reactions. This study presents an algorithm to estimate the kinetic parameters using Monte-Carlo (MC) simulation. The combination of an analytical reaction model and the MC simulation technique rapidly generates a large number of numerical values. Results show that MC-simulated data and experimental data are in fair agreement. Though the technique developed in this study proved to have potential, more experimental data are needed to check the robustness of the model.