An efficient computational analysis for stochastic fractional heroin model with artificial decay term

Heroin addiction is a continuously progressing phenomenon that represents a major problem for world's the public health; this indicates that the development of new methodologies to address the issue at the international level is a crucial priority. To study its transmission dynamics, a new stochastic fractional delayed heroin model based on stochastic fractional delay differential equations (SFDDEs) was developed to focus on the positive aspects of randomness and memory effects. The positive, boundedness, existence, and uniqueness of the model were studied rigorously. The equilibria (i.e., heroin-free equilibrium and the present equilibrium, which gives a clue about both eradication and persistence cases), reproduction number, and sensitivity of parameters were analyzed. The local and global stability of the new model was studied around its steady states. Also, well-known theorems are presented to investigate the extinction and persistence of heroin. The Grunwald-Letnikove non-standard finite difference (GL-NSFD) method was used for the efficient computational analysis of the stochastic fractional delayed model. For the dynamical consistency of the model, the positivity and boundedness of an efficient method were studied rigorously. The given study focuses on delay strategies and fractional calculus that could be useful in formulating specific measures for regulating addiction. Moreover, the simulated results support the theoretical analysis of the model and validate it.