Translator Disclaimer
2017 Analysis of a Heroin Epidemic Model with Saturated Treatment Function
Isaac Mwangi Wangari, Lewi Stone
J. Appl. Math. 2017: 1-21 (2017). DOI: 10.1155/2017/1953036


A mathematical model is developed that examines how heroin addiction spreads in society. The model is formulated to take into account the treatment of heroin users by incorporating a realistic functional form that “saturates” representing the limited availability of treatment. Bifurcation analysis reveals that the model has an intrinsic backward bifurcation whenever the saturation parameter is larger than a fixed threshold. We are particularly interested in studying the model’s global stability. In the absence of backward bifurcations, Lyapunov functions can often be found and used to prove global stability. However, in the presence of backward bifurcations, such Lyapunov functions may not exist or may be difficult to construct. We make use of the geometric approach to global stability to derive a condition that ensures that the system is globally asymptotically stable. Numerical simulations are also presented to give a more complete representation of the model dynamics. Sensitivity analysis performed by Latin hypercube sampling (LHS) suggests that the effective contact rate in the population, the relapse rate of heroin users undergoing treatment, and the extent of saturation of heroin users are mechanisms fuelling heroin epidemic proliferation.


Download Citation

Isaac Mwangi Wangari. Lewi Stone. "Analysis of a Heroin Epidemic Model with Saturated Treatment Function." J. Appl. Math. 2017 1 - 21, 2017.


Received: 12 May 2017; Accepted: 26 July 2017; Published: 2017
First available in Project Euclid: 19 September 2017

zbMATH: 07037462
MathSciNet: MR3696024
Digital Object Identifier: 10.1155/2017/1953036

Rights: Copyright © 2017 Hindawi


This article is only available to subscribers.
It is not available for individual sale.

Vol.2017 • 2017
Back to Top