In this paper, a mathematical model is proposed and analysed to assess the impacts of health care providers’ population on the transmission dynamics of COVID-19. The stability theory of differential equations is used to examine a mathematical model. The results of both local and global stability of disease-free equilibrium points were determined by using Routh-Hurwitz criteria and Metzler’s matrix method respectively, which verified that it was asymptotically stable. Also, the endemic equilibrium point was determined by the Lyapunov function which showed that E ∗ was globally asymptotically stable under strict conditions. The findings revealed that non-diagnosed and undetected health care providers seem to contribute to the high spread of COVID-19 in a population. Also, it illustrates that an increase in the number of non-diagnostic testing rates of health care providers may result in high infection rates in the community. Therefore, the particular study recommend that there is a necessity of applying early diagnostic testing to curtail the COVID-19 transmission in the health care providers’ population.