Dynamics of a stochastic heroin epidemic model with bilinear incidence and varying population size

We investigate a stochastic heroin epidemic model with bilinear incidence and varying population size. Sufficient criteria for the extinction of the drug abusers and the existence of ergodic stationary distribution for the model are established by constructing suitable stochastic Lyapunov functions. By analyzing the sensitivity of the threshold of spread, we obtain that prevention is better than cure. Numerical simulations are carried out to confirm the analytical results.