Worldwide, smoking is a common social practice, especially in places like schools and on important occasions. The World Health Organization (WHO) states that smoking is the third leading cause of death worldwide and the most significant avoidable cause of disease. Thus in this study, we present the solution behavior of the fractal-fractional (FF) smoking model. Using an effective numerical integration technique, the discretized system of FF differential equations that results is numerically integrated. By using RK4 to evaluate the numerical solution, we can meet the requirements for both efficiency and accuracy in the given approach. In the hope of providing some recommendations through which to reduce the risks of this bad behavior, the effect of some parameters and external factors affecting the solution behavior of this proposed mathematical model was studied, including the recruitment rate (due to immigration or birth) and the smoking cessation rate. The results show that the implemented technique is a straightforward and efficient tool for modeling the solution for these models.