Application of the NSFD method in a Malaria model with nonlinear incidence and recovery rates
Tapan Sarkar,Prashant K. Srivastava,Pankaj Biswas
DOI: https://doi.org/10.1140/epjp/s13360-024-05045-x
2024-03-15
The European Physical Journal Plus
Abstract:In this manuscript, we formulate a malaria model with nonlinear incidence and recovery rates to study the impact of the health system's available resources on disease spread and control. A detailed analysis of the model is provided. It is determined that when the basic reproduction number, , is less than unity, the disease may or may not die out due to nonlinear recovery rate. In addition, it is asserted that if the recovery rate is constant, disease eradication is possible for . The global stability of the unique endemic equilibrium point is established using geometric approach. The model system is also examined for Hopf bifurcation, and it is established that certain conditions on the transmission rate from vector to human lead to the emergence of periodic oscillations in the model system. Keeping in mind the nonlinear nature of the hypothesized model, we develop a non-standard finite difference (NSFD) scheme by discretizing the system. It is proven that conservation law and the positivity of the solutions are maintained by the proposed NSFD method for all finite step sizes. Furthermore, the convergence and error bounds of the developed schemes are also explored. To validate the analytical results, numerical simulations using our computational scheme are presented and these results are compared to two well-known standard numerical techniques, viz. the fourth-order Runge–Kutta (RK4) method and the forward Euler method. It is found that conventional numerical schemes fail to accurately capture the dynamics of the continuous model for certain step sizes, resulting in unstable and negative numerical solutions. In contrast, the developed NSFD scheme successfully retains the fundamental mathematical characteristics of the continuous model. It is worth mentioning that the NSFD schemes also capture the backward bifurcation phenomenon. In addition, numerical simulations show that Hopf bifurcation occurs due to transmission rate.
physics, multidisciplinary