A Discrete-time Mathematical Model of Smoking Dynamics with Two Sub-populations of Smokers

We analyze smoking dynamics with two sub-populations of smokers using a discrete-time mathematical model with a standard difference scheme. We divide the smokers' populations into beginners and heavy smokers. The boundedness of the solution is obtained. The equilibrium stability is assessed through the Jury stability conditions. We additionally present numerical simulations to validate the analytical results by giving several examples to depict the stability of all equilibriums. The sensitivity analysis of the model's parameters is performed to that can offer recommendations for regulators to reduce the number of smokers.