# American Institute of Mathematical Sciences

November  2020, 25(11): 4449-4477. doi: 10.3934/dcdsb.2020107

## Mathematical analysis of an age structured heroin-cocaine epidemic model

 1 Laboratoire d'Analyse Nonlinéaire et Mathématiques Appliquées, University of Tlemcen, Tlemcen 13000, Algeria 2 Graduate School of System Informatics, Kobe University, 1-1 Rokkodai-cho, Nada-ku, Kobe 657-8501, Japan

* Corresponding author

Received  December 2018 Revised  December 2019 Published  March 2020

This paper is devoted to studying the dynamics of a certain age structured heroin-cocaine epidemic model. More precisely, this model takes into account the following unknown variables: susceptible individuals, heroin users, cocaine users and recovered individuals. Each one of these classes can change or interact with others. In this paper, firstly, we give some results on the existence, uniqueness and positivity of solutions. Next, we obtain a threshold value $r(\Psi'[0])$ such that an endemic equilibrium exists if $r(\Psi'[0]) > 1$. We then show that if $r(\Psi'[0]) < 1$, then the disease-free equilibrium is globally asymptotically stable, whereas if $r(\Psi'[0]) > 1$, then the system is uniformly persistent. Moreover, for $r(\Psi'[0]) > 1$, we show that the endemic equilibrium is globally asymptotically stable under an additional assumption that epidemic parameters for heroin users and cocaine users are same. Finally, some numerical simulations are presented to illustrate our theoretical results.

Citation: Abdennasser Chekroun, Mohammed Nor Frioui, Toshikazu Kuniya, Tarik Mohammed Touaoula. Mathematical analysis of an age structured heroin-cocaine epidemic model. Discrete & Continuous Dynamical Systems - B, 2020, 25 (11) : 4449-4477. doi: 10.3934/dcdsb.2020107
##### References:

show all references

##### References:
Transfer diagram for model (1)
The functions $\beta$, $\theta_1$ and $\theta_2$ with respect to age $a$. (left) $\theta_1 \equiv \theta_2$ as in (40), (right) $\theta_1$ greater than $\theta_2$
The evolution of solution $S$ and $R$ with respect to time $t$ are drawn. The case of disease-free equilibrium with $r(\Psi'[0]) < 1$
The evolution of solutions $i_1$ and $i_2$ with respect to time $t$ and age $a$. The case of disease-free equilibrium with $r(\Psi'[0]) < 1$
The evolution of solution $S$ and $R$ with respect to time $t$ are drawn. The case of endemic equilibrium with $r(\Psi'[0]) > 1$
The evolution of solutions $i_1$ and $i_2$ with respect to time $t$ and age $a$. The case of endemic equilibrium with $r(\Psi'[0]) > 1$
Description of each symbol in model (1)
 Symbol Description $S(t)$ Density of susceptible individuals at time $t$ $i_1(t, \xi_1)$ Density of heroin users at time $t$ and age $\xi_1$ $i_2(t, \xi_2)$ Density of cocaine users at time $t$ and age $\xi_2$ $R(t)$ Density of recovered individuals at time $t$ $A$ Number of all newborns per unit time $\mu$ Natural death rate per capita and unit time $\beta(\xi_i)$ $(i=1, 2)$ Transmission rate for drug users with age $\xi_i$ $(i=1, 2)$ $\theta_1(\xi_1)$ Recovery rate from the consumption of heroin at age $\xi_1$ $\theta_2(\xi_2)$ Recovery rate from the consumption of cocaine at age $\xi_2$ $k_1$ Rate at which an individual recovered from the consumption of heroin becomes a cocaine user $k_2$ Rate at which an individual recovered from the consumption of cocaine becomes a heroin user $\delta_1$ Rate at which a recovered individual becomes a heroin user $\delta_2$ Rate at which a recovered individual becomes a cocaine user
 Symbol Description $S(t)$ Density of susceptible individuals at time $t$ $i_1(t, \xi_1)$ Density of heroin users at time $t$ and age $\xi_1$ $i_2(t, \xi_2)$ Density of cocaine users at time $t$ and age $\xi_2$ $R(t)$ Density of recovered individuals at time $t$ $A$ Number of all newborns per unit time $\mu$ Natural death rate per capita and unit time $\beta(\xi_i)$ $(i=1, 2)$ Transmission rate for drug users with age $\xi_i$ $(i=1, 2)$ $\theta_1(\xi_1)$ Recovery rate from the consumption of heroin at age $\xi_1$ $\theta_2(\xi_2)$ Recovery rate from the consumption of cocaine at age $\xi_2$ $k_1$ Rate at which an individual recovered from the consumption of heroin becomes a cocaine user $k_2$ Rate at which an individual recovered from the consumption of cocaine becomes a heroin user $\delta_1$ Rate at which a recovered individual becomes a heroin user $\delta_2$ Rate at which a recovered individual becomes a cocaine user
 [1] A. M. Elaiw, N. H. AlShamrani, A. Abdel-Aty, H. Dutta. Stability analysis of a general HIV dynamics model with multi-stages of infected cells and two routes of infection. Discrete & Continuous Dynamical Systems - S, 2020  doi: 10.3934/dcdss.2020441 [2] Laurent Di Menza, Virginie Joanne-Fabre. An age group model for the study of a population of trees. Discrete & Continuous Dynamical Systems - S, 2020  doi: 10.3934/dcdss.2020464 [3] Zhihua Liu, Yayun Wu, Xiangming Zhang. Existence of periodic wave trains for an age-structured model with diffusion. Discrete & Continuous Dynamical Systems - B, 2020  doi: 10.3934/dcdsb.2021009 [4] Yolanda Guerrero–Sánchez, Muhammad Umar, Zulqurnain Sabir, Juan L. G. Guirao, Muhammad Asif Zahoor Raja. Solving a class of biological HIV infection model of latently infected cells using heuristic approach. Discrete & Continuous Dynamical Systems - S, 2020  doi: 10.3934/dcdss.2020431 [5] Tin Phan, Bruce Pell, Amy E. Kendig, Elizabeth T. Borer, Yang Kuang. Rich dynamics of a simple delay host-pathogen model of cell-to-cell infection for plant virus. Discrete & Continuous Dynamical Systems - B, 2021, 26 (1) : 515-539. doi: 10.3934/dcdsb.2020261 [6] Eduard Feireisl, Elisabetta Rocca, Giulio Schimperna, Arghir Zarnescu. Weak sequential stability for a nonlinear model of nematic electrolytes. Discrete & Continuous Dynamical Systems - S, 2021, 14 (1) : 219-241. doi: 10.3934/dcdss.2020366 [7] Ting Liu, Guo-Bao Zhang. Global stability of traveling waves for a spatially discrete diffusion system with time delay. Electronic Research Archive, , () : -. doi: 10.3934/era.2021003 [8] Mohamed Dellal, Bachir Bar. Global analysis of a model of competition in the chemostat with internal inhibitor. Discrete & Continuous Dynamical Systems - B, 2021, 26 (2) : 1129-1148. doi: 10.3934/dcdsb.2020156 [9] Mohammad Ghani, Jingyu Li, Kaijun Zhang. Asymptotic stability of traveling fronts to a chemotaxis model with nonlinear diffusion. Discrete & Continuous Dynamical Systems - B, 2021  doi: 10.3934/dcdsb.2021017 [10] Hui Zhao, Zhengrong Liu, Yiren Chen. Global dynamics of a chemotaxis model with signal-dependent diffusion and sensitivity. Discrete & Continuous Dynamical Systems - B, 2020  doi: 10.3934/dcdsb.2021011 [11] Shujing Shi, Jicai Huang, Yang Kuang. Global dynamics in a tumor-immune model with an immune checkpoint inhibitor. Discrete & Continuous Dynamical Systems - B, 2021, 26 (2) : 1149-1170. doi: 10.3934/dcdsb.2020157 [12] Yanhong Zhang. Global attractors of two layer baroclinic quasi-geostrophic model. Discrete & Continuous Dynamical Systems - B, 2021  doi: 10.3934/dcdsb.2021023 [13] Lars Grüne. Computing Lyapunov functions using deep neural networks. Journal of Computational Dynamics, 2020  doi: 10.3934/jcd.2021006 [14] Peter Giesl, Sigurdur Hafstein. System specific triangulations for the construction of CPA Lyapunov functions. Discrete & Continuous Dynamical Systems - B, 2020  doi: 10.3934/dcdsb.2020378 [15] Chao Xing, Jiaojiao Pan, Hong Luo. Stability and dynamic transition of a toxin-producing phytoplankton-zooplankton model with additional food. Communications on Pure & Applied Analysis, 2021, 20 (1) : 427-448. doi: 10.3934/cpaa.2020275 [16] Lu Xu, Chunlai Mu, Qiao Xin. Global boundedness of solutions to the two-dimensional forager-exploiter model with logistic source. Discrete & Continuous Dynamical Systems - A, 2020  doi: 10.3934/dcds.2020396 [17] Pedro Branco. A post-quantum UC-commitment scheme in the global random oracle model from code-based assumptions. Advances in Mathematics of Communications, 2021, 15 (1) : 113-130. doi: 10.3934/amc.2020046 [18] Imam Wijaya, Hirofumi Notsu. Stability estimates and a Lagrange-Galerkin scheme for a Navier-Stokes type model of flow in non-homogeneous porous media. Discrete & Continuous Dynamical Systems - S, 2021, 14 (3) : 1197-1212. doi: 10.3934/dcdss.2020234 [19] Peter Giesl, Zachary Langhorne, Carlos Argáez, Sigurdur Hafstein. Computing complete Lyapunov functions for discrete-time dynamical systems. Discrete & Continuous Dynamical Systems - B, 2021, 26 (1) : 299-336. doi: 10.3934/dcdsb.2020331 [20] Yifan Chen, Thomas Y. Hou. Function approximation via the subsampled Poincaré inequality. Discrete & Continuous Dynamical Systems - A, 2021, 41 (1) : 169-199. doi: 10.3934/dcds.2020296

2019 Impact Factor: 1.27