Stability analysis of five-grade Leishmania epidemic model with harmonic mean-type incidence rate

In this paper, we discuss the Anthroponotic Cutaneous Leishmania transmission. In addition, we develop a mathematical model for the Anthroponotic Cutaneous Leishmania transmission and consider its qualitative behavior. We derive the threshold number R0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$R_{0}$\end{document} of the model using the next generation method. In the disease-free case, we carry out the local and global stability under the condition R0<1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$R_{0}<1$\end{document}. Moreover, we derive the global stability at the disease-free equilibrium point by utilizing the Castillo-Chavez method. On the other hand, at the endemic equilibrium point, we show the local and global stability to be held under specific conditions and R0>1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$R_{0}>1$\end{document}. We also establish the global stability at the endemic equilibrium point with the help of a geometrical approach, which is a generalization of Lyapunov theory, by using a second additive compound matrix. Finally, we take into account the sensitivity analysis of the threshold number with other parameters. We also discuss several graphs of important parameters.