Numerical solutions of fractional order rabies mathematical model via Newton polynomial

The recent outbreak of rabies virus has affected numerous individuals in the community, highlighting the importance of studying the disease mathematically in epidemiology. In this paper, we develop a mathematical model for the spread of rabies disease using the harmonic mean incidence rate and determine the reproduction number R0 using the next generation matrix approach. The fractional dynamics of the model incorporate the interaction between infected individuals and environmental factors. We analyze the model using both Atangana-Baleanu-Caputo (ABC) and Caputo-Fabrizio (CF) techniques, drawing on both modern and classical approaches. For qualitative investigation using fractional operators, we utilize the Banach fixed-point theory and derive the Hyers-Ulam stability concept. We assign values to the model parameters and utilize the Newton interpolation technique to obtain a numerical scheme. Additionally, we perform sensitivity analysis on the model. In conclusion, our findings provide important insights for the study of epidemics. 2023 THE AUTHORS. Published by Elsevier BV on behalf of Faculty of Engineering, Alexandria University This is an open access article under the CC BY-NC-ND license