Abstract
This paper investigates a new mathematical SQIR model for COVID-19 by means of four dimensions; susceptible, quarantine, infected and recovered (SQIR) via Non-linear Saturated Incidence Rate. First of all the model is formulated in the form of differential equations. Disease-free, endemic equilibriums and Basic Reproduction Number are found for the said model. Local Stability is analyzed through Jacobean Matrix while Lyapunov Function is constructed for the study of Global Stability of the Model. Using nonstandard finite difference method, numerical results are simulated. By Simulation, we mean how protection, exposure, death and cure rates affect the Susceptible, Quarantined, Infected and recovered population with the passage of time.