RESUMO
The fundamental goal of this research is to suggest a novel mathematical operator for modeling visceral leishmaniasis, specifically the Caputo fractional-order derivative. By utilizing the Fractional Euler Method, we were able to simulate the dynamics of the fractional visceral leishmaniasis model, evaluate the stability of the equilibrium point, and devise a treatment strategy for the disease. The endemic and disease-free equilibrium points are studied as symmetrical components of the proposed dynamical model, together with their stabilities. It was shown that the fractional calculus model was more accurate in representing the situation under investigation than the classical framework at α = 0.99 and α = 0.98. We provide justification for the usage of fractional models in mathematical modeling by comparing results to real-world data and finding that the new fractional formalism more accurately mimics reality than did the classical framework. Additional research in the future into the fractional model and the impact of vaccinations and medications is necessary to discover the most effective methods of disease control.