Stability and numerical analysis via non-standard finite difference scheme of a nonlinear classical and fractional order model

In this paper, we develop a new mathematical model for an in-depth understanding of COVID-19 (Omicron variant). The mathematical study of an omicron variant of the corona virus is discussed. In this new Omicron model, we used idea of dividing infected compartment further into more classes i.e asymptomatic, symptomatic and Omicron infected compartment. Model is asymptotically locally stable whenever [Formula: see text] and when [Formula: see text] at disease free equilibrium the system is globally asymptotically stable. Local stability is investigated with Jacobian matrix and with Lyapunov function global stability is analyzed. Moreover basic reduction number is calculated through next generation matrix and numerical analysis will be used to verify the model with real data. We consider also the this model under fractional order derivative. We use Grunwald–Letnikov concept to establish a numerical scheme. We use nonstandard finite difference (NSFD) scheme to simulate the results. Graphical presentations are given corresponding to classical and fractional order derivative. According to our graphical results for the model with numerical parameters, the population’s risk of infection can be reduced by adhering to the WHO’s suggestions, which include keeping social distances, wearing facemasks, washing one’s hands, avoiding crowds, etc.