Modelling Drug Abuse Epidemics in the Presence of Limited Rehabilitation Capacity

The abuse of drugs is now an epidemic globally whose control has been mainly through rehabilitation. The demand for drug abuse rehabilitation has not been matched with the available capacity resulting in limited placement of addicts into rehabilitation. In this paper, we model limited rehabilitation through the Hill function incorporated into a system of nonlinear ordinary differential equations. Not every member of the community is equally likely to embark on drug use, risk structure is included to help differentiate those more likely (high risk) to abuse drugs and those less likely (low risk) to abuse drugs. It is shown that the model has multiple equilibria, and using the centre manifold theory, the model exhibits the phenomenon of backward bifurcation whose implications to rehabilitation are discussed. Sensitivity analysis and numerical simulations are performed. The results show that saturation in rehabilitation will in the long run lead to the escalation of drug abuse. This means that limited access to rehabilitation has negative implications in the fight against drug abuse where rehabilitation is the main form of control. This suggests that increased access to rehabilitation is likely to lower the drug abuse epidemic.