Global stability and modeling with a non-singular kernel for fractional order heroin epidemic model: Insights from different population studies

There has been a worldwide epidemic of heroin that has affected people, families, societies, and cultures across the world. Now, the heroin epidemic has transitioned from heroin abuse to the overuse of synthetic narcotics, which are widely accessible and inexpensively produced. In this work, a novel mathematical approach is applied to investigate the dynamics of the heroin epidemic model and its harmful effect on society with different population data. A heroin model has been constructed with the importance of a non-singular kernel in the sense of a generalized Mittag-Leffler kernel. The well-posedness of the proposed model is proven via fixed-point theory. To examine the heroin model, two equilibrium states have been determined. These equilibrium states are proven to be locally and globally asymptotically stable. To analyze the behavior of heroin, a basic reproduction number and sensitivity analysis are used to determine the impact of different parameters mathematically as well as through simulations. To find the approximate solution, we implement the Toufik–Atangana numerical method at different fractional order values. The sensitivity of the heroin model is carried out, and 3-D graphs show the significance of the parameter involved in the model. Finally, the numerical outcomes are presented with different values of fractional parameters.