Global Dynamics of Drug Use Model and Its Optimal Control Analysis


Ratchada Viriyapong Sirapong Nongnamtip


        In this paper, we propose and analyze a mathematical model of drug use involving rehabilitation. With our model, we assume that the use of drugs can be initiated by three groups of people; light drug users, heavy drug users and drug users under rehabilitation. In this model, the nonnegativity and boundedness of solutions are verified, and two equilibrium points (drug-free and drug-endemic) are obtained. The basic reproduction number is calculated.  We show that each equilibrium point is stable locally and globally under some conditions. Further, an optimal control problem is applied to the model by adding three control variables; awareness and educational program control, family and friends care control, and rehabilitation campaign control. The numerical simulation of this optimal control model is performed, and the results show that each control alone could already reduce the number of drug users for some certain amount, however, a combination of all three controls gives the best result in reducing overall drug use.  

Keywords: drug use, awareness and educational program, family and friends care, basic reproduction number, optimal control


Binuyo, A. O., & Osuntokun, O. T. (2021). Mathematical modelling of the addiction of drug substance among students in tertiary institutions in Nigeria. Journal of Biomedical Research & Environmental Sciences, 2(9), 851-856.
Butler, G., Freedman, H. I., & Waltman, P. (1986). Uniformly persistent systems. Proceedings of the American Mathematical Society, 96(3), 425-430.
Caulkins, J. P., Gragnani, A., Feichtinger, G., & Trangler, G. (2006). High and low-frequency oscillations in drug epidemics. International Journal of Bifurcation and Chaos, 16(11), 3275-3289.
Castillo-Chaves, C., Feng, Z., & Huang, W. (2001). On the computation of R0 and its role in global stability. Retrieved from;
Everingham, S. S., & Rydell, C. P. (1994). Modeling the Demand of Cocaine. Drug Policy Research Centre.
Freedman, H. I., Ruan, S., & Tang, M. (1994). Uniform persistence and flows near a closed positively invariantset. Journal of Dynamics and Differential Equations, 6(4), 583-600.
Hafiruddin, M., Fatmawati, F., & Miswanto, M. (2019). Mathematical model analysis of a drug. AIP Conference Proceedings, 2192, 0600101-10. http:// 10.1063/1.5139156
Islam, A. I., & Biswas, H. A. (2020). Optimal control strategy applied to dynamic model of drug abuse incident for reducing its adverse effects. medRxiv.
Johnston, L. D., O'Malley, P. M., & Bachman, J. G. (1988). Illicit drug use and smoking and drinking by America's high school students and college students and young adults 1975-1987. U.S.: Department of Health and Human Services, Public Health Service.
Kanato, M., & Leyatikul, P. M. P. (2019). Size estimation of substances users population in Thailand 2019. ONCB Journal, 36(2), 37-48. Retrieved from
Kanyaa, J. K., Osman, S., & Wainaina, M. (2018). Mathematical modelling of substance abuse by commercial drivers. Global Journal of Pure and Applied Mathematics, 14(9), 1149-1165. Retrieved from
Li, J., & Ma, M. (2018). The analysis of a drug transmission model with family education and public health education. Infectious Disease Modelling, 3(2018), 74-84.
Li, M. Y., & Muldowney J. S. (1993). On Bendixson’s criterion. Journal of Differential Equations, 106(1), 27-39.
Li, M. Y., & Muldowney, J. S. (1996). A geometric approach to global-stablilty problems. SIAM Journal on Mathematical Analysis, 27(4), 1070-1083.
Mushayabasa, S., & Tapedzesa, G. (2015). Modeling illicit drug use dynamics and its optimal control analysis. Hindawi Corporation Computational and Mathematical Methods in Medicine, 2015, 383154.
Njagarah, J. B. H, & Nyabadza, F. (2013). Modelling the role of drug barons on the prevalence of drug epidemics. Mathematical Biosciences and Engineering, 10(3), 843-860.
Nyabadza, F., & Hove-Musekwa, S. D. (2010). From heroin epidemics to methamphetamine epidemics: Modelling substance abuse in a South African province. Mathematical Biosciences, 225(2010), 132-140.
Olajide, A. J. (2020). Asymptotic stability of illicit drug dynamics with banditry compartment. Applied Mathematics & Information Sciences, 14(5), 791-800.
Pontryagin, L. S., Boltyanskii, V. G., Gamkrelidze, R. V., & Mishchenko, E. F. (1986). In L. S. Pontryagin (Ed), The Mathematical Theory of Optimal Processes. New York: Gordon and BS Publishers.
Rossi, C. (2004). Operational models for epidemics of problematic drug use: The mover-stayer approach to heterogeneity. Socio-Economic Planning Sciences, 38(1), 73-90.
Samanta, G. P. (2011). Dynamic behaviour for a nonautonomous heroin epidemic model with time delay. Journal of Applied Mathematics and Computing, 35,161–178.
UNODC. (2021). UNODC World Drug Report 2021: pandemic effects ramp up drug risks, as youth underestimate cannabis dangers. Retrieved from nabis-dangers.html
van den Driessche, P., & Watmough, J. (2002). Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission. Mathematical Biosciences, 180(1-2), 29-48. https:
White, E., & Comiskey, C. (2007). Heroin epidemics, treatment and ODE modelling. Mathematical Biosciences, 208(1), 312-324.
World Drug Report. (2021). United Nations. Retrieved from

Research Articles


How to Cite
VIRIYAPONG, Ratchada; NONGNAMTIP, Sirapong. Global Dynamics of Drug Use Model and Its Optimal Control Analysis. Naresuan University Journal: Science and Technology (NUJST), [S.l.], v. 31, n. 3, p. 48-66, aug. 2023. ISSN 2539-553X. Available at: <>. Date accessed: 17 june 2024. doi: