The Modeling Dynamics of HIV and CD4$$^{}$$ T-cells During Primary Infection in Fractional Order: Numerical Simulation

This article is devoted to introduce a numerical treatment using Adams–Bashforth–Moulton method of the fractional model of HIV-1 infection of CD4(^{+}) T-cells. We study the effect of the changing the average number of viral particles N with different sets of initial conditions on the dynamics of the presented model. The fractional derivative is described in Caputo sense. Special attention is given to present the local stability of the proposed model using fractional Routh–Hurwitz stability criterion. Qualitative results show that the model has two equilibria: the disease-free equilibrium and the endemic equilibrium points. We compare our numerical solutions with those numerical solutions using fourth-order Runge–Kutta method (RK4). The obtained numerical results of the proposed model show the simplicity and the efficiency of the proposed method.