Pyrolysis kinetics of a hardwood representative, beech (Fagus sylvatica), was investigated by two different kinetic approaches: model-free isoconversional method and model-fitting method. The model-free isoconversional method was used for the determination of apparent kinetic parameters, i.e. the activation energy and pre-exponential factor. The model fitting method was used for the optimization of kinetic parameters of the reaction pathways of three selected reaction mechanisms: one-step, two-step, and three-step one. In both approaches, thermo-gravimetric data were used at five heating rates: 2A degrees C min(-1), 5A degrees C min(-1), 10A degrees C min(-1), 15A degrees C min(-1) and 20A degrees C min(-1). As the most suitable mechanism, the three-step mechanism containing the intermediate degradation step was chosen. This selection was supported by experimental results from the C-13 NMR analysis of solid residues prepared at the key temperatures within the range of 230-500A degrees C. The progress of mass fraction values of each component in this mechanism was simulated. Conclusions from the simulation were confronted with experimental results from the C-13 NMR.