Abstract
In this paper, we formulate a. (n+3)-dimensional nonlinear virus dynamicsmodel that considers n-stages of the infected cells and (n+1) distributed time delays. The model incorporates humoral immune response and general nonlinear forms for the incidence rate of infection, the generation and removal rates of the cells and viruses. Under a set of conditions on the general functions, the basic infection reproduction number R-0(M) and the humoral immune response activation number R-1(M) are derived. Utilizing Lyapunov functionals and LaSalle's invariance principle, the global asymptotic stability of all steady states of themodel are proven. Numerical simulations are carried out to confirm the theoretical results. Copyright (C)2016 JohnWiley & Sons, Ltd.