This study presents a mathematical model that represents the population growth dynamics of tobacco consumers based on a system of ordinary nonlinear differential equations. The model is used to determine the Basic Reproductive Number (R 0). Two points of equilibrium are found and their local stability is classified. Finally, the Matlab software is used to present numerical simulations using the fourth-order Runge-Kutta method, and it is shown that the solutions approach an asymptotically stable point, under the variation of R 0 .