Disproportionation of toluene is one of the key processes in the direct production of aromatic products and the intermediate production of polymer resins, synthetic materials and plasticizers. The aim of this research is to model the kinetics of the disproportionation reaction of toluene with zeolite catalyst. In this research, the Langmuir-Hinshelwood kinetic model was selected and used to solve the model and estimate the kinetic parameters of the model from experimental data using MATLAB code. The selected kinetic model was solved using MATLAB programming (using the MATLAB ODE45 library function and the 4th-order Runge-Kutta method), and kinetic parameters were obtained using optimization. In this study, it was used to predict the behavior of toluene and xylene gases in a plug Reactor. The model was presented and after optimizing the parameters related to the reaction model, it was able to predict the rate of conversion of the reactants at the reactor outlet well. Finally, the findings of the study showed that: the reaction model, along with process modeling, were able to predict the rate of conversion of the reactants at the reactor outlet well. It was also found that with increasing temperature and keeping other parameters such as pressure, gas composition percentage and reactor characteristics constant; the rate of conversion of toluene and xylene increases. Increasing the pressure, if other parameters are constant, increases the percentage of xylene and toluene conversion, and increasing the percentage of toluene and xylene at the reactor inlet increases the percentage of toluene conversion and xylene conversion, respectively. Also The average error in the mole fraction of toluene exiting the reactor was 2.96 and the average error in the mole fraction of xylene exiting the reactor was 7.53. Also, a plug reactor was simulated using the selected kinetic model in the HISIS software and the mole fractions of toluene and xylene at the reactor outlet were compared w