Abstract
This manuscript investigates fractal-fractional order smoking models with relapse and harmonic mean type incidence rate under the Caputo derivative. We derive the existence and unique results about the solution for the considered model via fixed point theory. For the stability of the considered system, Ulam-Hyers (UH) approach is used. We compute the numerical solution by using fractional Adams-Bashforth method. For the simulation of the model, we consider different values of fractional order delta and fractal dimension theta by using some real values of the parameters. The proposed scheme is used to simulate the available data for some smoking community including potential, light, and quit smokers. Various graphical presentations are given to understand the dynamics of the model at various fractional orders.