In this manuscript, an epidemic model employed the dynamics of drugs usage among adults. Among smokers, often the desire to quit smoking arises. A large number of smokers attempt to quit, but only a few of them are successful. A non-linear mathematical model is employed to study and assess the dynamics of smoking and its impact on public health in a community. We prove the essential properties, bounded, positivity and well-posed, also local and global stability analysis has been made for the epidemic model. The sensitivity analysis of the model is provided by threshold or reproductive number as well as analyzed qualitatively. We develop an unconditionally convergent nonstandard finite difference scheme by applying Mickens approach φ(h) = h + O(h 2) instead of h to control the spread of bad impact in society. Finally numerical simulations are also established to investigate the influence of the system parameters on the spread of the smoking impact in society.