A nonlinear mathematical model of COVID-19 containing asymptomatic as well as symptomatic classes of infected individuals is considered and examined in the current paper. The largest eigenvalue of the next-generation matrix known as the reproductive number is obtained for the model, and serves as an epidemic indicator. To better understand the dynamic behavior of the continuous model, the unconditionally stable nonstandard finite difference (NSFD) scheme is constructed. The aim of developing the NSFD scheme for differential equations is its dynamic reliability, which means discretizing the continuous model that retains important dynamic properties such as positivity of solutions and its convergence to equilibria of the continuous model for all finite step sizes. The Schur–Cohn criterion is used to address the local stability of disease-free and endemic equilibria for the NSFD scheme; however, global stability is determined by using Lyapunov function theory. We perform numerical simulations using various values of some key parameters to see more characteristics of the state variables and to support our theoretical findings. The numerical simulations confirm that the discrete NSFD scheme maintains all the dynamic features of the continuous model.
