The impact of standard and nonstandard finite difference schemes on HIV nonlinear dynamical model

• Published in: Chaos, solitons and fractals solitons and fractals, 2023-08, Vol.173, p.113755, Article 113755
• Main Authors: Li, Shuo, Bukhsh, Imam, Khan, Ihsan Ullah, Asjad, Muhammad Imran, Eldin, Sayed M., El-Rahman, Magda Abd, Baleanu, Dumitru
• Format: Article
• Language: English
• Subjects: HIV mathematical model; Local and global stability; Lyapunov function; Routh-Hurwitz criterion; SFD and NSFD schemes
• ISSN: 0960-0779; 1873-2887
• DOI: 10.1016/j.chaos.2023.113755

Mathematical models are enormously valuable in recognition the characteristics of infectious afflictions. The present study describes and analyses a nonlinear Susceptible-Infected (S·I) type mathematical model for HIV/AIDS. To better comprehend the dynamics of disease diffusion, it is assumed that by giving AIDS patients timely Anti Retroviral Therapy (ART), their transition into HIV infected class is attainable. The ART treatment can reduce or manage the spread of disease among individuals that can extend their life for some more years. For the model, the basic reproduction number is formed which provides a base to study the stability of disease free and endemic equilibria. To understand the entire dynamical behavior of the model, standard finite difference (SFD) schemes such as Runge-Kutta of order four (RK-4) and forward Euler schemes and nonstandard finite difference (NSFD) scheme are implemented. The goal of constructing the NSFD scheme for differential equations is to ensure that it is dynamically reliable, while maintaining important dynamical properties like the positivity of the solutions and its convergence to equilibria of continuous model for all finite step sizes. However, the essential characteristics of the continuous model cannot be properly maintained by the Euler and RK-4 schemes, leading to the possibility of numerical solutions that are not entirely similar to those of the original model. For the NSFD scheme, the Routh-Hurwitz criterion is used to assess the local stability of disease-free and endemic equilibria. To explain the global stability of both the equilibria, Lyapunov functions are offered. To verify the theoretical findings and validate the dynamical aspects of the abovementioned schemes, numerical simulations are also provided. The outcomes offered in this study may be engaged as an effective tool for forecasting the progression of HIV/AIDS epidemic diseases. •A nonlinear Susceptible-Infected (SI) type mathematical model for HIV/AIDS is presented and examined.•The standard finite difference (SFD) schemes, which include forward Euler scheme and Runge-Kutta of order four (RK-4) scheme as well as nonstandard finite difference (NSFD) scheme are designed to understand the entire dynamical behavior of the model.•The Routh-Hurwitz criterion is used to assess the local stability of disease-free and endemic equilibria.•The Lyapunov functions are provided to discuss the global stability of both the equilibria.•To validate the theoretic