# Mathematical Model for Coronavirus Disease 2019 (COVID-19) Containing Isolation Class

Authors: Anwar Zeb; Ebraheem Alzahrani; Vedat Suat Erturk; Gul Zaman
Journal: BioMed Research International
 (2020)

---

## Abstract

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.


---

## Body

## 1. Introduction

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 Figure1) [1].Figure 1
BBC News up to 8 May [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.Table 1
Parameters and description.
SymbolsDescriptionSSusceptible populationEExposed populationIInfected populationQIsolated populationRRecovered populationβ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 rateUnfortunately, 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–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.The paper is organised as follows: Section2 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.


## 2. Model Formulation

In this section, we develop the mathematical model by taking into account the above assumptions.(1)dStdt=A−μSt−βNStEt+It,dEtdt=βNStEt+It−πEt−μ+γEt,dItdt=πEt−σIt−μIt,dQtdt=γEt+σIt−θQt−μQt,dRtdt=θQt−μRt,where the parameters and variables used are described in Table 1.As the first four equations are independent ofRt, so omit without generality the last equation for Rt and the modified system (1) becomes
(2)dStdt=A−μSt−βNStEt+It,dEtdt=βNStEt+It−πEt−μ+γEt,dItdt=πEt−σIt−μIt,dQtdt=γEt+σIt−θQt−μQt.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
(3)dsdt=μ−μs−βNse+i,(4)dedt=βNse+i−πe−μ+γe,(5)didt=πe−σi−μi,(6)dqdt=γe+σi−θq−μq,with the initial conditions
(7)s0=s0≥0,e0=e0≥0,i0=i0≥0,q0=qo≥0.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).Lemma 1.
Under the initial conditions (7), all the solutions st,et,it,qt of system (3) remain nonnegative for t≥0Proof.
By the initial conditions (7), it was discovered that
(8)dsdts=0=μ>0,dedte=0=βNsi≥0,didti=0=πe≥0,dqdtq=0=γe+σi≥0.


## 3. Local Stability and Existence of Positive Equilibrium Point

The existence of unique positive equilibrium and stability of system (3) depends on the basic reproductive number R0 on free equilibrium point (FEP) C0, which is determined with the help of the next generation matrix method [8]. The free coronavirus equilibrium point is C0=1,0,0,0.Consider the following matrices for finding the basic reproductive numberR0:
(9)F=βNse+i0,V=πe+μ+γe−πe+σi+μi.Now Jacobian ofF and V at C0 are
(10)F=βNβN00,V=π+μ+γ0−πσ+μ.The dominant eigenvalue ofFV−1 represents R0=ρFV−1, which is
(11)R0=βNσ+μ+ππ+μ+γσ+μ.Theorem 2.
The system (3) is locally stable related to virus-free equilibrium point C0, R0<1 and unstable if R0>1.Proof.
For local stability atC0, the Jacobian of system (3) is
(12)J=E0=−μ−βN−βN00βN−π−μ−γβN00π−σ+μ00γσ−θ+μ,which follows the eigenvalues , λ2<0, λ3<0, and λ4<0, if R0<1. So the system (3) is locally stable if R0<1 and unstable if R0>1.Theorem 3.
There exists a unique positive virus equilibrium pointC∗s∗,e∗,i∗,q∗ for system (3), if R0>1.Proof.
By letting the right hand sides of all equations of system (3) to zero, as
(13)μ−μs−βNse+i=0,βNse+i−πe−μ+γe=0,πe−σi−μi=0,γe+σi−θq−μq=0implies that
(14)s∗=1R0,e∗=σ+μπi∗,i∗=πμR0−1βNπ+σ+μ,q∗=γσ+μ+πσπθ+μi∗.
From the value ofi∗, it is obvious that all the values of s∗,e∗,q∗ are positive if R0>1.Theorem 4.
IfR0<1, then the system (3) is globally stable.Proof.
For the proof of this theorem, first, we construct the Lyapunov functionL as
(15)L=e−e0+βNσ+μi−i0.
Differentiating equation (15) with respect to time and keeping the reality in mind that R0<1 and 0<s<1, we obtained
(16)L′=é+βNσ+μí,L′=βNse+βNsi−π+μ+γe+βNσ+μπe−βNσ+μσ+μi,L′=βNse+βNsi−π+μ+γe+βNσ+μπe−βNi,≤βNe−π+μ+γe+βNσ+μπe=R0−1e.Therefore, if R0<1, then L′<0, which implies that the system (3) is globally stable for R0<1.


## 4. Numerical Method and Results

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:
(17)y′k=ft,y1,y2,⋯,yn,where k=1,2,⋯,n, then, by the NFSD method
(18)y1′=y1,k+1−y1,kh,y2′=y2,k+1−y2,kh,⋮yn′=yn,k+1−yn,kh.Now, using the NSFD method for numerical solution of system (3), it follows that
(19)sk+1=hμ+sk1+hβNkik+hβNkek+hμsk+1,ek+1=hβNksk+1ik+hβNksk+1ek+ek1+hπ+μ+γ,ik+1=hπek+1+ik1+hσ+μ,qk+1=hγek+1+hσik+1+qk1+hθ+μ.We assume the parameters of the system (3) shown in Table 1 [4].Figure2 shows the solutions for St, Et, It, Qt, and Rt obtained by NSFD, RK4, and ode45 for R0>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 be controlled and so the current pandemic; otherwise, no control is possible. For R0<1, when the contact rate becomes smaller, then, the current infectious disease may be controlled (see Figure 3).Figure 2
Plots present the susceptible, exposed, infected, isolated, and recovered population whenR0>1.Figure 3
Plots present the susceptible, exposed, infected, isolated and recovered population whenR0<1.


## 5. Conclusion

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.



---

*Source: 3452402-2020-06-29.xml*