The purpose of this work is to investigate the boronizing kinetics of AISI M2 steel by using the integral diffusion model with consideration of boride incubation periods. This simulation model was established by solving the differential algebraic equations (DAE) resulting from the integral method in the temperature range of 1173 to 1323 K. By using a particular solution of the obtained DAE system, the values of boron diffusivities in the FeB and Fe2B layers were estimated. The estimated values of activation energies for boron diffusion in AISI M2 steel were respectively 228.06 kJ mol-1 and 212.10 for FeB and Fe2B. Finally, a comparison was made between the simulated thicknesses of FeB and Fe2B layers and the experimental values obtained at 1173, 1223, 1273 and 1323 K for 10 h. The findings of this research work may serve as a tool to simulate the boronizing kinetics of any steel with a microstructure consisting of FeB and Fe2B layers versus the boriding parameters (the time and the temperature).