In this study, we explored novel fractional Visceral Leishmaniasis (VL) disease models with Caputo fractional derivatives. The computational solution of the fractional VL model is derived using the Generalized Euler's method (GEM). The compartment structure of the five-dimensional model includes susceptible, infected, public treatment, private treatment, and recovered classes. This study considers two treatment strategies: public and private. The fractional-order novel model was graphically simulated using MATLAB. Furthermore, the stability of the equilibrium points was investigated. The findings obtained using the proposed method are computationally efficient and easy to implement.
Khan et al. (Sun,) studied this question.