This paper presents a mathematical model that studies the importance of treatment at rehabilitation centers intending to show the impact of the rehabilitation centers in minimizing illicit drug usage. We consider the global stability of endemic equilibrium using the properties of Volterra-Lyapunov matrices. Numerical simulations were carried out to show the impact of rehabilitation centers on illicit drug users.