Research Article  Open Access
Mathematical Model for Coronavirus Disease 2019 (COVID19) Containing Isolation Class
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 COVID19 infection means you gain longlasting 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 COVID19 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 RungeKutta 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 COVID19. Therefore, isolation of the infected human overall can reduce the risk of future COVID19 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, COVID19 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 COVID19 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. Coronavirusconfirmed 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 COVID19. Cao et al. [2, 3] studied the clinical features of coronavirus and discussed the shortterm outcomes of 18 patients and 102 patients with COVID19 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 uptodate 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 COVID19 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 COVID19 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 COVID19. Therefore, isolation of the infected human overall can reduce the risk of future COVID19 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 RungeKutta 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 COVID19.
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 virusfree 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).
5. Conclusion
In this work, we presented that isolation of the infected human overall can reduce the risk of future COVID19 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 numberrelated 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. View at: Google Scholar
 J. L. Cao, X. R. Hu, W. J. Tu, and Q. Liu, “Clinical features and shortterm outcomes of 18 patients with corona virus disease 2019 in intensive care unit,” Intensive Care Medicine, vol. 46, no. 5, pp. 851–853, 2020. View at: Publisher Site  Google Scholar
 J. L. Cao, W. J. Tu, X. R. Hu, and Q. Liu, “Clinical features and shortterm outcomes of 102 patients with Coronavirus disease 2019 in Wuhan, China,” Clinical Infectious Diseases, 2020. View at: Publisher Site  Google Scholar
 W. Ming, J. V. Huang, and C. J. P. Zhang, “Breaking down of the healthcare system: mathematical modelling for controlling the novel coronavirus (2019nCoV) outbreak in Wuhan, China,” medRxiv and bioRxiv, 2020. View at: Publisher Site  Google Scholar
 I. Nesteruk, “Statisticsbased predictions of coronavirus epidemic spreading in Mainland China,” Innovative Biosystems and Bioengineering, vol. 4, no. 1, pp. 13–18, 2020. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar
Copyright
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.