BioMed Research International
Volume 2020 (2020), Article ID 3452402, 7 pages
https://doi.org/10.1155/2020/3452402

Research Article

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

Anwar Zeb,¹ Ebraheem Alzahrani,² Vedat Suat Erturk,³ and Gul Zaman⁴

¹Department of Mathematics, COMSATS University Islamabad, Abbottabad Campus, Abbottabad 22060, Khyber Pakhtunkhwa, Pakistan
²Department of Mathematics, Faculty of Science, King Abdulaziz University, P.O. Box 80203, Jeddah 21589, Saudi Arabia
³Department of Mathematics, Faculty of Arts and Sciences, Ondokuz Mayis University, 55139 Samsun, Turkey
⁴Department of Mathematics, University of Malakand, Chakdara 18000, Dir (Lower), Khyber Pakhtunkhwa, Pakistan

Correspondence should be addressed to Anwar Zeb; anwar55.ciit@yahoo.com

Received 11 April 2020; Revised 10 May 2020; Accepted 27 May 2020; Published 29 June 2020

Academic Editor: Wen-Jun Tu

Copyright © 2020 Anwar Zeb et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

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.

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 Figure 1) [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.

Symbols
Description

Susceptible population

Exposed population

Infected population

Isolated population

Recovered 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 rate

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–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: 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.

2. Model Formulation

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 , so omit without generality the last equation for and the modified system (1) becomes

For system (2), let , , , , and , 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).

Lemma 1. Under the initial conditions (7), all the solutions of system (3) remain nonnegative for

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

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 on free equilibrium point (FEP) , which is determined with the help of the next generation matrix method [8]. The free coronavirus equilibrium point is .

Consider the following matrices for finding the basic reproductive number :

Now Jacobian of and at are

The dominant eigenvalue of represents , which is

Theorem 2. The system (3) is locally stable related to virus-free equilibrium point , and unstable if .

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

which follows the eigenvalues , , , and , if . So the system (3) is locally stable if and unstable if .

Theorem 3. There exists a unique positive virus equilibrium point for system (3), if .

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

implies that

From the value of , it is obvious that all the values of are positive if .

Theorem 4. If , then the system (3) is globally stable.

Proof. For the proof of this theorem, first, we construct the Lyapunov function as

Differentiating equation (15) with respect to time and keeping the reality in mind that and , we obtained

Therefore, if , then , which implies that the system (3) is globally stable for .

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:

where , then, by the NFSD method

Now, using the NSFD method for numerical solution of system (3), it follows that

We assume the parameters of the system (3) shown in Table 1 [4].

Figure 2 shows the solutions for , , , , and obtained by NSFD, RK4, and ode45 for , 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 , 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 when .

Figure 3: Plots present the susceptible, exposed, infected, isolated and recovered population when .

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.

Data Availability

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

Conflicts of Interest

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

Authors’ Contributions

The authors equally contributed in preparing this manuscript.

Acknowledgments

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.

References

Johns Hopkins University, “Internet source national health agencies, report up to 08 May,” Tech. Rep.

J. L. Cao, X. R. Hu, W. J. Tu, and Q. Liu, “Clinical features and short-term outcomes of 18 patients with corona virus disease 2019 in intensive care unit,” Intensive Care Medicine, vol. 46, no. 5, pp. 851–853, 2020.

J. L. Cao, W. J. Tu, X. R. Hu, and Q. Liu, “Clinical features and short-term outcomes of 102 patients with Coronavirus disease 2019 in Wuhan, China,” Clinical Infectious Diseases, 2020.

W. Ming, J. V. Huang, and C. J. P. Zhang, “Breaking down of the healthcare system: mathematical modelling for controlling the novel coronavirus (2019-nCoV) outbreak in Wuhan, China,” medRxiv and bioRxiv, 2020.

I. Nesteruk, “Statistics-based predictions of coronavirus epidemic spreading in Mainland China,” Innovative Biosystems and Bioengineering, vol. 4, no. 1, pp. 13–18, 2020.

M. Batista, Estimation of the final size of the coronavirus epidemic by SIR model, Research Gate, 2020.

V. A. Okhuese, Mathematical predictions forcoronavirus as a global pendamic, Research Gate, 2020.

P. Van den Driessche and J. Watmough, “Further notes on the basic reproduction number,” in Mathematical Epidemiology, F. Brauer, P. Driessche, and J. Wu, Eds., Lecture Notes in Mathematics, pp. 159–178, 2008.

R. E. Mickens, Nonstandard finite difference methods for differential equations, World Scientific, Singapore, 1994.

R. E. Mickens, “Calculation of denominator functions for nonstandard finite difference schemes for differential equations satisfying a positivity condition,” Numerical Methods for Partial Differential Equations, vol. 23, no. 3, pp. 672–691, 2007.