key: cord-0838729-38aewt2p
authors: Zeb, Anwar; Alzahrani, Ebraheem; Erturk, Vedat Suat; Zaman, Gul
title: Mathematical Model for Coronavirus Disease 2019 (COVID-19) Containing Isolation Class
date: 2020-06-25
journal: Biomed Res Int
DOI: 10.1155/2020/3452402
sha: e8cee2eeb61285eb69783efc45784cd56dddfe85
doc_id: 838729
cord_uid: 38aewt2p

The deadly coronavirus continues to spread across the globe, and mathematical models can be used to show suspected, recovered, and deceased coronavirus patients, as well as how many people have been tested. Researchers still do not know definitively whether surviving a COVID-19 infection means you gain long-lasting immunity and, if so, for how long? In order to understand, we think that this study may lead to better guessing the spread of this pandemic in future. We develop a mathematical model to present the dynamical behavior of COVID-19 infection by incorporating isolation class. First, the formulation of model is proposed; then, positivity of the model is discussed. The local stability and global stability of proposed model are presented, which depended on the basic reproductive. For the numerical solution of the proposed model, the nonstandard finite difference (NSFD) scheme and Runge-Kutta fourth order method are used. Finally, some graphical results are presented. Our findings show that human to human contact is the potential cause of outbreaks of COVID-19. Therefore, isolation of the infected human overall can reduce the risk of future COVID-19 spread.

Mathematical models are useful to understand the behavior of an infection when it enters a community and investigate under which conditions it will be wiped out or continued. Currently, COVID-19 is of great concern to researches, governments, and all people because of the high rate of the infection spread and the significant number of deaths that occurred. In December 2019, coronavirus first reported in Wuhan, China, is an infectious disease caused by a newly discovered coronavirus. The virus that causes COVID-19 is mainly transmitted through droplets generated when an infected person coughs, sneezes, or exhales. These droplets are too heavy to hang in the air and quickly fall on floors or surfaces. Coronavirus-confirmed cases reached nearly four million in 187 countries, and approximately 295,000 people have lost their lives due to this virus.

According to figures collated by Johns Hopkins University, the largest cases occurred at the US. Noting that more than 77,000 deaths happened, it also has the world's highest death toll (see the Figure 1 ) [1] .

Researchers have been tracking the spread of the virus, have mobilized to speed innovative diagnostics, and are working on a number of vaccines to protect against COVID-19. Cao et al. [2, 3] studied the clinical features of coronavirus and discussed the short-term outcomes of 18 patients and 102 patients with COVID-19 in intensive care units. For the demographic details of 102 patients, see Table 1 [3] . Coronaviruses are typically transmitted from person to person through respiratory droplets and close contact. The majority of the transmission is happening through respiratory droplets that we may inhale from close contact with one another [4] . A modified SIR epidemic model is presented in [4] to project the actual number of infected cases and the specific burdens on isolation wards and intensive care units. Nesteruk [5] developed an SIR (susceptible, infected, and recovered) epidemic model and discussed statistically the parameters used in the proposed model and showed how to control this infection.

Unfortunately, the number of coronavirus victims is expected to be much higher than that predicted on February 10, 2020, since 12289 new cases (not previously included in official counts) have been added two days later. Further research should focus on updating the predictions with the use of up-to-date data and using more complicated mathematical models. Currently, there are no licensed vaccines or therapeutic agents for coronavirus prevention or treatment though research studies into potential antivirals and vaccine candidates are underway in a number of countries. Vaccine testing, development, and distribution are typically a much longer process than drug development, and it is not likely that a vaccine for COVID-19 will be ready before 2021 at the earliest. The virus can easily spread in dense places. Social distancing or low contact rate refers to measures that are taken to increase the physical space between people to slow down the spread of the virus. Batista [6] studied the logistic growth regression model which is used for the estimation of the final size of the coronavirus epidemic. Several researchers developed different models of COVID-19 and studied dynamical behavior (see for instance [4] [5] [6] [7] ).

From the above discussion, it was concluded that human to human contact is the potential cause of outbreaks of COVID-19. Therefore, isolation of the infected human overall can reduce the risk of future COVID-19 spread. In order to do this, we divided the total population into five compartments: susceptible, exposed, infected, isolated, and recovered from the disease. This study will lead to the mathematical model formulation in which the interaction of the exposed population and infected population occurred to the susceptible populations. The infected individuals, the individuals showing no symptoms apparently but have the disease in weak form inside their bodies, must be sent to isolated class in different rates. The local stability and global stability of model is discussed, by using the approach of the basic reproductive. For the numerical solution of the proposed model, the nonstandard finite difference (NSFD) scheme and Runge-Kutta fourth order method are used. Finally, some graphical results are presented. Our findings show that human to human contact is the potential cause of outbreaks of COVID-19. β Rate at which susceptible population moves to infected and exposed class π Rate at which exposed population moves to infected one γ Presents the rate at which exposed people take onside as isolated σ Shows the rate at which infected people were added to isolated individual θ Rate at which isolated persons recovered μ Natural death rate plus disease-related death rate 2 BioMed Research International

The paper is organised as follows: Section 2 is related to the model formulation keeping in mind the assumptions that exposed and infected individuals are making contacts with susceptible individuals at the same rate. Section 3 is concerned with the local stability and existence of positive equilibrium solution. Some numerical simulations are executed to illustrate the analytical results in Section 4. Finally, conclusions are presented in Section 5. Time (day)

The population of exposed NSDF RK4 ODE45 Figure 2 : Plots present the susceptible, exposed, infected, isolated, and recovered population when R 0 > 1. 

The population of recovered × 10 −4 Figure 3 : Plots present the susceptible, exposed, infected, isolated and recovered population when R 0 < 1.

BioMed Research International

In this section, we develop the mathematical model by taking into account the above assumptions.

where the parameters and variables used are described in Table 1 .

As the first four equations are independent of RðtÞ, so omit without generality the last equation for RðtÞ and the modified system (1) becomes

For system (2) , let N = A/μ, s = S/N, e = E/N, i = I/N, and q = Q/N, and rescale the system (2) to get the normalized form

with the initial conditions

In the remaining sections, we will discuss the local and global stability of the proposed model with initial conditions. First, a result is observed for the positivity and boundedness of the solution of system (3). (7), all the solutions ðsðtÞ, eðtÞ, iðtÞ, qðtÞÞ of system (3) remain nonnegative for t ≥ 0

Proof. By the initial conditions (7), it was discovered that

The existence of unique positive equilibrium and stability of system (3) depends on the basic reproductive number R 0 on free equilibrium point (FEP) C 0 , which is determined with the help of the next generation matrix method [8] . The free coronavirus equilibrium point is C 0 = ð1, 0, 0, 0Þ. Consider the following matrices for finding the basic reproductive number R 0 :

Now Jacobian of F and V at C 0 are

The dominant eigenvalue of FV −1 represents R 0 = ρ ðFV −1 Þ, which is

Theorem 2. The system (3) is locally stable related to virus-free equilibrium point C 0 , R 0 < 1 and unstable if R 0 > 1.

Proof. For local stability at C 0 , the Jacobian of system (3) is

which follows the eigenvalues , λ 2 < 0, λ 3 < 0, and λ 4 < 0, if R 0 < 1. So the system (3) is locally stable if R 0 < 1 and unstable if R 0 > 1.

There exists a unique positive virus equilibrium point C * ðs * , e * , i * , q * Þ for system (3), if R 0 > 1.

Proof. By letting the right hand sides of all equations of system (3) to zero, as

implies that

From the value of i * , it is obvious that all the values of ðs * , e * , q * Þ are positive if R 0 > 1. Proof. For the proof of this theorem, first, we construct the Lyapunov function L as

Differentiating equation (15) with respect to time and keeping the reality in mind that R 0 < 1 and 0 < s < 1, we obtained

Therefore, if R 0 < 1, then L ′ < 0, which implies that the system (3) is globally stable for R 0 < 1.

The NSFD method is used for the numerical solution of the proposed model (3) . Basically, NSFD is an iterative method in which we get closer to solution through iteration [9, 10] . Let nonstandard ODEs be given as follows:

where k = 1, 2, ⋯, n, then, by the NFSD method

Now, using the NSFD method for numerical solution of system (3), it follows that s k+1 = hμ + s k 1 + hβN k i k + hβN k e k + hμs k + 1 ð Þ , e k+1 = hβN k s k+1 i k + hβN k s k+1 e k + e k 1 + h π + μ + γ ð Þ ,

We assume the parameters of the system (3) shown in Table 1 [4] . Figure 2 shows the solutions for SðtÞ, EðtÞ, IðtÞ, QðtÞ, and RðtÞ obtained by NSFD, RK4, and ode45 for R 0 > 1, which show that it is unstable and will never become stable because of contact rates of infected people to susceptible people. So according to the act of Governments to keeping people within specified bounds (may be their home (stay at home and save your lives), offices, etc.), the contact rate will 6 BioMed Research International be controlled and so the current pandemic; otherwise, no control is possible. For R 0 < 1, when the contact rate becomes smaller, then, the current infectious disease may be controlled (see Figure 3 ).

In this work, we presented that isolation of the infected human overall can reduce the risk of future COVID-19 spread. Our model shows that the coronavirus spreads through contact and describes how fast something changes by counting the number of people who are infected and the likelihood of new infections. Those new infections are what induce the epidemic. For this reason, we think that this research may lead to better guessing of the spread of this pandemic in the future. This paper is devoted to implement the coronavirus mathematical model containing isolation class. The reproductive number-related stability is discussed, which showed the impact of interaction of infected people to susceptible population and proved graphically and analytically that if we control this contact rate, the control of the current disease is possible, otherwise. State and territory governments have different restrictions in place for public gatherings. Therefore, citizens need to follow the directions from time to time to minimize the health risk.

The authors confirm that the data supporting the findings of this study are available within the article.

Internet source national health agencies, report up to 08

Clinical features and short-term outcomes of 18 patients with corona virus disease 2019 in intensive care unit

Clinical features and short-term outcomes of 102 patients with Coronavirus disease

Breaking down of the healthcare system: mathematical modelling for controlling the novel coronavirus (2019-nCoV) outbreak in Wuhan, China

Statistics-based predictions of coronavirus epidemic spreading in Mainland China

Estimation of the final size of the coronavirus epidemic by SIR model

Mathematical predictions forcoronavirus as a global pendamic

Further notes on the basic reproduction number

Nonstandard finite difference methods for differential equations

Calculation of denominator functions for nonstandard finite difference schemes for differential equations satisfying a positivity condition

The authors declare that there is no conflict of interest regarding the publication of this paper.

This project was funded by the Deanship of Scientific Research at King Abdulaziz University, Jeddah. The authors, therefore, gratefully acknowledge DSR for the technical and financial support.

The authors equally contributed in preparing this manuscript.