Qualitative behavior of a highly non-linear Cutaneous Leishmania epidemic model under convex incidence rate with real data.

In the context of this investigation, we introduce an innovative mathematical model designed to elucidate the intricate dynamics underlying the transmission of Anthroponotic Cutaneous Leishmania. This model offers a comprehensive exploration of the qualitative characteristics associated with the transmission process. Employing the next-generation method, we deduce the threshold value $ R_0 $ for this model. We rigorously explore both local and global stability conditions in the disease-free scenario, contingent upon $ R_0 $ being less than unity. Furthermore, we elucidate the global stability at the disease-free equilibrium point by leveraging the Castillo-Chavez method. In contrast, at the endemic equilibrium point, we establish conditions for local and global stability, when $ R_0 $ exceeds unity. To achieve global stability at the endemic equilibrium, we employ a geometric approach, a Lyapunov theory extension, incorporating a secondary additive compound matrix. Additionally, we conduct sensitivity analysis to assess the impact of various parameters on the threshold number. Employing center manifold theory, we delve into bifurcation analysis. Estimation of parameter values is carried out using least squares curve fitting techniques. Finally, we present a comprehensive discussion with graphical representation of key parameters in the concluding section of the paper.