Numerical analysis of a bi-modal covid-19 SITR model

This study presents a structure preserving nonstandard finite difference scheme to analyze a susceptible-infected-treatment-recovered (SITR) dynamical model of coronavirus 2019 (covid-19) with bimodal virus transmission in susceptible population. The underlying model incorporates the possible treatment measures as the emerging scenario of covid-19 vaccines. Keeping in view the fact that the real time data for covid-19 is updated at discrete time steps, we propose a new structure preserving numerical scheme for the proposed model. The proposed numerical scheme produces realistic solutions of the complex bi-modal SITR nonlinear model, converges unconditionally to steady states and reflects dynamical consistency with continuous sense of the model. The analysis of the model reveals that the model remains stable at the steady state points. The basic reproduction number [Formula: see text] falls less than 1 when treatment rate is increased and disease will die out. On the other hand, it predicts that human population may face devastating effects of pandemic if the treatment measures are not strictly implemented.