Abstract
In this article, a model is developed to depict the dynamics of gonorrhea and human immunodeficiency virus coinfection with disability and mortality risks. Estimating the parameter values and validating the model are done using real-world data on cases of gonorrhea. The properties of solutions, such as boundedness and positivity, uniqueness and existence are investigated. The basic reproduction number is calculated using the next-generation matrix. The stability of Gonorrhea-Free and Present Equilibrium states is also analyzed. The Lyapunov function is used to verify global stability at the Gonorrhea-Free Equilibrium, while the graph-theoretic approach is employed to analyze global stability at the Gonorrhea-Present Equilibrium. Sensitivity indices are calculated to identify the significant parameters transmitting the disease. The Adams Bashforth predictor-corrector scheme is used to simulate the behavior of all classes individually and in combination with the different classes, considering the effect of fractional order ϕ. The relation between various parameters and the basic reproduction number is analyzed and portrayed. The limitations on the values of the parameters to ensure the basic reproduction number is below one are also discussed. This model incorporates aspects that cause disabilities and mortality from both gonorrhea and HIV and estimates the contribution of both to the long-term neurological and reproductive disabilities. The simulation results demonstrate how the course of disabilities is altered when these aspects are modified. This work highlights the importance of early diagnosis and disability prevention efforts.