Abstract
The objective of this work is to investigate the qualitative behavior of an Human Immunodeficiency Virus (HIV) dynamics model with two types of cocirculating target cells and under the effect of anti-viral drug therapy. The model takes into account both short-lived infected cells and long-lived chronically infected cells. In the two types of target cells, the drug efficacy is assumed to be different. The incidence rate of virus infection is given by general functional response. We have derived the basic reproduction number which determines the global dynamics of the model. We have established a set of conditions which are sufficient to investigate the global stability of the equilibria of the model. The global stability analysis of the model has been established using the Lyapunov method. Numerical simulations have been performed for the model with a specific form of the incidence rate function. We have shown that the numerical and theoretical results are consistent.