Abstract
Malaria, a life-threatening disease responsible for millions of deaths worldwide, remains a major public health issue, especially in under-resourced regions. It is caused by Plasmodium parasites, transmitted through mosquito bites, and disproportionately affects vulnerable groups like children and pregnant women. To improve understanding and management of malaria transmission, we investigated different mathematical models, traditionally based on integer-order derivatives. In this study, we introduced a novel approach using a fractional-order mathematical model to evaluate how treatment strategies impact malaria's spread. Initially, we modeled limited treatment scenarios with integer-order nonlinear differential equations. However, recognizing the complexity of malaria dynamics, we enhanced the model with fractional-order derivatives and power laws to capture a more detailed picture of disease behavior. The research established conditions for solution existence and uniqueness within the fractional framework and assessed the stability of the endemic equilibrium using the Lyapunov function technique. A sensitivity analysis of the basic reproduction number identified key factors influencing malaria transmission. Using the fractional Adams-Bashforth-Moulton method, we simulated various scenarios to explore the effects of model parameters and fractional-order values. Visual tools like surface and contour plots helped illustrate the findings. The results showed that improving treatment strategies and implementing preventive measures, such as mosquito control and timely medication, significantly reduced malaria cases. On the other hand, factors like increased mosquito contact and ineffective treatments aggravated the disease's impact. This study provided valuable insights into malaria dynamics, highlighting the critical need for sustained efforts in treatment and prevention to mitigate its devastating effects on communities.