key: cord-0715819-j4d1knwz
authors: Din, Anwarud; Li, Yongjin; Khan, Tahir; Zaman, Gul
title: Mathematical analysis of spread and control of the novel corona virus (COVID-19) in China
date: 2020-09-23
journal: Chaos Solitons Fractals
DOI: 10.1016/j.chaos.2020.110286
sha: 9d825fd853f99e15d815abdf487a9526c01f3a40
doc_id: 715819
cord_uid: j4d1knwz

Number of well-known contagious diseases exist around the world that mainly include HIV, Hepatitis B, influenzas etc., among these, a recently contested coronavirus (COVID-19) is a serious class of such transmissible syndromes. Abundant scientific evidence the wild animals are believed to be the primary hosts of the virus. Majority of such cases are considered to be human-to-human transmission, while a few are due to wild animals-to-human transmission and substantial burdens on healthcare system following this spread. To understand the dynamical behavior such diseases, we fitted a susceptible-infectious-quarantined model for human cases with constant proportions. We proposed a model that provide better constraints on understanding the climaxes of such unseen disastrous spread, relevant consequences, and suggesting future imperative strategies need to be adopted. The main features of the work include the positivity, boundedness, existence and uniqueness of solution of the model. The conditions were derived under which the COVID-19 may extinct or persist in the population. Sensitivity and estimation of those important parameters have been carried out that plays key role in the transmission mechanism. To optimize the spread of such disease, we present a control problem for further analysis using two control measures. The necessary conditions have been derived using the Pontryagin’s maximum principle. Parameter values have been estimated from the real data and experimental numerical simulations are presented for comparison as well as verification of theoretical results. The obtained numerical results also present the verification, accuracy, validation, and robustness of the proposed scheme.

On the basis of genetic properties the coronaviranae family comprises of four genera, which include genus Betacoronavirus, genus Gammacoronavirus and genus Deltacoronavirus etc., [1] . RNA genome of the coronavirus is the largest RNA among the all reported virus RNA ranging from 26 to 32 kb [2] . The coronavirus can infect aves reptiles mammals including homo sapiens. The coronavirus mainly infect subclinical [1] [2] [3] . The severe respiratory syndrome (SARS-CoV) and Middle East respiratory syndrome coronavirus (MERS-CoV) are zoonotic pathogens which are related to the genus of Betacoronavirus causing severe respiratory infections in man. The COVID-19 can rapidly adapt to the new host and ecological niches and attain the tendency for genetic recombination and strengthen the inherently mutation rate of typical RNA viruses [4, 5] . The typical structural representation of the COVID-19 is shown in Figure 1 highlighting various segments of proteins namely E-Protein, S-Protein, M-Protein, and HE-Protein [3] . [3] .

Not yet the definite source of this virus is still under consideration. Apart, the infections in family and health workers clusters confirm the spread from human to human, while the mood of transmission is still open for confirmation. However, on January 21, 2020, WHO approved human to human transmission [3, 5, 6] .

Till March 16, 2020, a total of 81009 confirm cases were reported in Mainland China due to pneumonia/COVID-19 and the same time new patients were also reported worldwide including Thailand, South Korea, Japan, United Kingdom, USA, India etc. The COVID-19 appeared to cause mild from viral pneumonia with low capability from person to person transmission which is different from SARS-CoV [4] . The attitude behind this maybe because of the recombination within receptor-binding of glycoprotein of the virus originated form the snake, wildlife or from both. Accordingly, the COVID-19 will likely be attenuated upon infection to humans. Still the concern about the adaptation in humans may acquire more competence to efficiently replicate and provide feasible route transmit rapidly via close contacts with the infected bodies [7] .

In order to identify the mechanism from which CoVs infected patients are passing, we developed a flow diagram in Figure 2 . In the first step, an infected patient is identified. Further, we will use any traceable app to track back the record of locations recently visited by the positively diagnosed patients. For this purpose, CDR analysis can be used with minimum 2 − 3 weeks of records. Followed by this, the awareness about this disease will be increased among the public. For this aim, we shall spread the news of the result of the second step through social media platforms such as apps. It will also lead to further trace the chain of the infected people. Furthermore, we will use any easiest and possible method to get the chain of people who interacted with patient and the places visited by them. The identified confirmed coronavirus infected patients will comply with strict health and security monitoring. Finally, the same strategy shall be used for all infected patients.

To properly know the mechanism of transmission and hence control the COVID-19 mathematically, we will formulate a model with the help of available literature. The current study is actually divided into three main parts; the statistic, dynamic and control. In the first part, we just studied the qualitative aspects of the disease and estimated the key rates form the real data. In the dynamics part, the authors have tried to answer questions like when the disease will dies out of the Hubie province? When does it will persist in the population? Which parameters are more responsible for the disease spreading? Mathematically, we will calculate the possible equilibria of the model and its stability analysis will be carried out using the methods of linearization, the Lyapunov theory and geometrical approach. The last part is devoted to control theory of COVID-19. In this part, we have chosen two control measures and formulated a control problem for further mathematical analysis. To obtain an optimal control problem along with characterization of control variables, using the Pontrygin's principle. The obtained analytical results and real data will be compared with numerical solutions using standard numerical method.

Formal organization of the manuscript are as: In the 2nd section, the statistics analysis of Coronavirus (COVID-19) will be presented. The data to be used for the study is till March 2, 2020. We proposed the Coronavirus (COVID-19) model in section 3 using the previous knowledge on epidemic modeling. Further, we deliberate the properties of existence, positivity, boundedness of solution to the problem and biologically feasibility. In section 4, the reproductive number will be calculated and it will be used for analyzing both local and global behavior of the disease free state. The dynamics with local properties as well as global properties around endemic state will be carried out in section 5. Also, local sensitivity analysis will be a part of this section. In section 6, we will present numerical simulation for supporting theoretical results related to dynamical behavior of the problem under discussion. Section 7 is devoted to the the formulation of the optimal control problem, development of suitable objective functional, existence of control variables, optimality condition and characterization of the control variables. The numerical results will be derived both for control and without control problems and a comparison may be presented between real data and the obtained findings. Finally, section 8 and 9 will be presented which are discussion and conclusion of the study, respectively.

The 2019 novel coronavirus( COVID-19) was identified by using the next-generation sequencing. Until 28th January,2020, about 5900 confirmed cases and about 9000 suspected cased were reported via Chinese authorities due to COVID-19 in various region of Mainland China [7] . In the same time various cases of COVID-19 have been reported in other countries like Italy, Iran, Japan, Spain, USA, Germany, South Korea and Singapore etc. The medical worker and family member were also reported infected causing human to human transmission upon exposure to infected person of COVID-19. The usual symptoms of novel COVID in an infected subject are high fever, some had dyspnea, while the chest radiographs reveals invasive lesions in lungs [5] [6] [7] . The epidemiological curve in Figure 3 for the worldwide case reported and upon data implementation, can suggests the trend of infected, death and recovered cases from COVID-19.

On December 31, 2019 Wuhan health authority release a public notice about the outbreak after confirmation, with no evidence from human to human transmission of the virus and suggest that the disease is preventable and controllable after implementing proper precautionary measurements while going out. The WHO was also formally informed about the outbreak shortly [7] . Later on January 20, 2020 scientists from the Chinas National Health Commission expressed their potential from one human to another upon the confirmation of two cases in Guangdong, upon the visit of infected family visited to Wuhan. Later the Wuhan authorities also announced new measurements for the control and prevention of this epidemic by canceling Chinese New Year celebrations, and start examine of body temperature upon traveling in public transport on 14th January, 2020. Beside all these for the general public safeguard a quarantine period was introduce on January 23rd, 2020 and lock-down Wuhan was initiated [7, 8] . On January 26th, 2020, a high level tasked force was initiated for the prevention and control of newly rise epidemic by COVID-19, which decided to implement new measurements for the control and prevention aiding extension of spring festival holidays [9] . Besides on 1st February, Huanggang, Hubei authorities permitted single person for purchasing necessary household items on alternative days. Then same measurements was implemented in other parts of the Mainland China. Till 20th of February, 2020, the number of reported cases were 75465 in China according to NRS [10, 11] .

NRS updated the record on daily basis for the cumulative confirmed, as well as suspected, dead and new cases. Beside Hubie, each province hold a media talk at 1300hr daily and update about the CoVID-19 and compare with the previous day [10, 11] . From Figure 4a , the curves for epidemics are drawn with the help of Infectious Disease Information System (IDIS) data Chinas National, which need every case of covid and therefore have been reported electronically by authorities immediately after diagnosed. It contain cases which are reported as asymptomatic and data is updated in real time, while individually cases are reported and downloaded after every 2400hr. Figure 4b represents the epidemiological curve of the reported cases of Hubei province. The curves includes the graphical presentation of the infected, death and quarantine cases [5] [6] [7] . 

We discuss transmission dynamics of the novel coronavirus a mathematical model according to the component of the virus transmittal [12] [13] [14] . We subdivide the human host populace (T(t)) in three groups: susceptible individuals S(t) as vulnerable for the infection; infected I(t) i.e., confirmed infected; Q(t) those who have contact history with an infectious individual and may or may not be infected with COVID-19 and whose activities have been restricted for a specific period of time. SIQ models are appropriate to model COVID-19 as quarantine destabilizes the epidemic and lead to sustained oscillations in the dynamics of the disease. The COVID-19 model whose the transmission structure is depicted in the Figure 5 . Moreover, we also put some assumptions before presenting the proposed model in the follow-ing form: a 1 . All parameters as well as variable and state variables are positive as well as non-negative respectively. a 2 . The susceptible agents goes to the infection classes and there is a constant inflow into the susceptible population. a 3 . Initially infected or suspected people move to quarantined class and confirmed cases from quarantined come back to the infected compartment.

Keeping in view the 5 and assumptions a 1 − a 3 we can describe the dynamics of COVID-19 in Hubie in the form of the system of three differential equations:

with initial conditions

where the parameters detail and its associated description are explained in the We discuss the well possedness of the proposed model in the form of the following axioms.

where

Clearly, C ≥ 0 and L preserve the axioms of Metzler matrix, so system (1) is invariant in R 3 + . Proposition 2. The solution of (1) i.e., (S, I, Q) with (2) are positive.

Proof. It could be clearly noted that the solution of the first equation of system (1) becomes

The solution of equation (4) is

∀ t > 0, which shows that S(t) > 0. Similarly the second equation of (1) gives the form

which implies that I(t) ≥ 0. Continuing the same process it is very simple to prove that Q(t) is also positive. Thus (S, I, Q) is nonnegative. Proof. Since

The differentiation of the above equation (6) gives

Clearly, dT dt + d 0 T ≤ Λ. Solving we obtain

which is biologically feasible. Moreover every solution in ∆ remains in ∆ for t ≥ 0.

Proof. It is very much clear that all the states variables of the proposed problem are nonnegative, so the problem as stated by (1)- (2) is well-posed and biologically feasible. From

So every solution of (1) along with (2) in ∆ remains in ∆.

Let

The disease-free state (E 1 ) takes the form

This disease free state is used to calculate the threshold parameter (R 0 ), also called the basic reproduction number i.e., the average of secondary number of infections, see for detail [15] [16] [17] [18] . Moreover, the threshold parameter (R 0 ) is used in the calculating of the endemic state. We follow the Watmough et al., method for the purposes of calculating the threshold parameter, whose detail is given in [14] . We know that I * > 0, so

By rearranging the terms, we can write I in the following way

The term R 0 used in (12) and so called the threshold number or threshold quantity which is given by

Proof. The three eigen-values of the Jacobean matrix at DFE are

Clearly, the first two eigenvalues are negative, whereas, the third eigenvalue is negative only if R 0 < 1. Hence the proof.

Proof. We define a function is given by

Differentiating F(t), we get

Then

The simplification and some re-arrangements gives

Thus, when R 0 < 1, then dF dt < 0. Also, dF dt = 0 iff S(t) = S 1 , I(t) = 0 and Q(t) = 0, which proves the conclusion.

We discuss the endemic state (EE) in this section. To prove the stability of model (1) at endemic equilibrium, let E * 2 is the EE and then by the use of Eq. (10) and Eq. (13) we arrive

where

Clearly, the relations (15) exists if R 0 > 1. (15) exists for the problem as stated by (1) .

The local dynamics at EE (15) is discussed via the following theorem.

Theorem 2. If R 0 > 1, then system (1) is locally asymptotically stable at the endemic equilibrium E * 2 defined in (15) . Otherwise, (1) is unstable.

Proof. Let J * is the Jacobian matrix of (1) at E * 2 (15), then

For obtaining the characteristic polynomail, we set |J * − λA|=0, where A is identity matrix . Thus

Expansion of the determinant in (17) gives us

We noted that one of the eigenvalue of the matrix is λ 1 = −b. For calculating the values of other eigenvalues, we will consider the quadratic equation

By using the values of the endemic equilibrium, we ca write the above equation in the form

Now by Descartes' rules of sign, if R 0 > 1 then equation (18) has two negative real roots or complex with negative real parts. Hence all of the eigenvalues are negative and hence the theorem.

Theorem 3. Assumption of R 0 > 1 guarantees that endemic equilibrium E * 2 = (S * , I * , Q * ) of model (1) is globally asymptotically stable and unstable otherwise.

Proof. Let J 2 and J |2| 2 be the Jacobian matrix and second additive compound matrix of the system containing only the model (1) and using Eq.(10). Then

Next, we will take into account the function P(χ) = P(S,

Thus,

Taking

For simplicity, we write B, such that

Let us assume that (X 1 , X 2 , X 3 ) be a vector in R 3 with the norm · defined by

Considering the Lozinski measure (B) [19, 20] with endowed norm defined above

where g i = X ij + (X ii ) for i, j = 1, 2 and i = j, which implies that

Since, X 12 = σ and X 21 = η, therefore, g 1 and g 2 becomes

Thus, we can also write Eq. (19) as

As a result, we get

Which proves the conclusion, that the proposed model is globally stable around the endemic equilibrium.

The local sensitivity analysis discuss the relation between the threshold quantity (basic reproductive number) and the model parameters. To find these indices, we use the formula Z Π = Π R 0 ∂R 0 ∂Π , which was used by many authors, e.g., see [21, 22] . This analysis provides the changes in the values of the epidemic parameters to the value of the threshold parameter and then to the disease control and spreading. For simplicity, we suppose that

We perform the local sensitivity analysis using Z Π = Π R 0 ∂R 0 ∂Π along with the following parameter's values Λ = 0.03805333333, η = 0.144211141, , γ = 0.00594474, d 0 = 0.007121000000, µ = 0.007121000000, k = 0.027864676, and σ = 0.0052281. We get the following sensitivity indices

It is worthy noted that Eq. (20) demonstrates the indices of η, σ, γ, µ. On the basis of these it is easy to find the essential parameters related to the disease dominance and spreading. The above equations states that γ is proportional to threshold parameter in a direct way. Therefore increasing in this parameter leads to increasing the value of threshold quantity (R 0 ). On the other side η, µ and σ are proportional to (R 0 ) in inverse direction. This indicate that increase will make decrease in the value of threshold quantity. Thus it is easy to formulate a control analysis on the basis of above sensitivity analysis

We carry out the simulation of the model to verify the previous analytical results with the help of graphical representations. We use 4th order Runge-Kutta technique. The application of Rungekutta method of order 4th on the proposed model leads to the following system

Step 1. For S(0) = 0,

Step 2. For i = 1, 2, 3, 4, ..., n − 1,

,

Step 3. For i = 1, 2, 3, ..., n − 1, writing S * (t i ) = S * , I * (t i ) = I * , Q * (t i ) = Q * . We used the available data from section 2 and estimated the values of some important parameters like γ, η, σ, k. Whereas, the values of other parameters such as Λ, d 0 and µ were taken from the available sources. Precisely, the following set of parameters is considered

where its corresponding numerical values and sources are presented in Table 2 For this set of parameter values, the value of R 0 and S 0 were calculated 0.2261469439 and 5.343818752, respectively. The simulation was carried out for this set of data as shown in Figure 6 .

By using the values shown in Table 2 , sample simulation were carried out for susceptible population. We have consider four different initial population of susceptible individual, that is, S 0 = 58.498998, 65.498998, 50.498998, 70.498998 where the population was considered in million and 58.498998 million is the actual population of the Hubie province. Figure 6a shows that all of the susceptible population will approach to 5.343818752 irrespective of its initial values whenever R 0 < 1. In the case of R 0 < 1, each solution curve S(t) almost taking 550 days in order to reach to its equilibrium value S 0 = 5.343818752. It means that if the we wish to eliminate the disease from the community, still it will take almost one and half year. Further, the figures shows that the disease will effect a major portion of the population during the indicated course of outbreak.

The dynamics of infected population in case of R 0 < 1 was plotted in Figure 6b for initial infected population I(t) = 0.001000 actual, 0.01000, 0.051000, 1.01000 millions. It is clear from the plot that the if we increase the initially infected people, it will markedly contribute to infected population at initial stages. However, the infection may completely eliminated from the community after 250 days irrespective of the initial infected population if we kept R 0 < 1. It means that if we able to reduce the basic reproduction rate was decreased to less than one, the disease will dies out from the province after more than eight months. In the same way, we plotted Q(t) for different values of Q(0), i.e., 0.000002 actual, 0.00002, 0.042, 2.000002 millions for R 0 < 1 in Figure 6c . During initial period of the infection, quarantine people may be increased if we increased the initial value of Q(0). However, in long run (almost 400 days) all of the solution curves Q(t) approaches zero irrespective of its initial values. Thus, for elimination of COVID-19 form Hubie it is necessary to keep the value of R 0 < 1. During initial phase of outbreak, more people have close contact history with infected so one must quarantine major portion of the population. However, as the time evolves, the quarantine people will tend to decline and eventually goes to zero. Further, Figures 6a-6c verify our theoretical findings that the disease-free equilibrium is locally and globally asymptotically stable if and only if R 0 < 1. Next, we assumed the real value of Λ, i.e., Λ = 0.3805333333 and use the same values of other parameters as in Table 2 and study the dynamical analysis of the model under discussion numerically. The corresponding values of R 0 and endemic equilibrium was calculated as follows:

The global analysis of all of the classes were presented in Figure 7 for R 0 > 1. In Figure 7a , we have consider again the same initial population of susceptible individual, i.e., 58.498998, 65.498998, 50.498998, 70.498998 millions, however, the value of R 0 is now greater than unity. It is clear from the figure, that the susceptible population will tend to decrease/increase during initial 250 days of the infection and afterward it will show no effect. That is, in long run the susceptible population will stabilize itself and will approaches 23.62985172. Figure 7b shows that due to an increase in the initially infected people a sudden increase may be observed later on, i.e., initial 50-80 days. The infection will increase in first 50 days and afterward decreases gradually. However, using the current epidemic status in the province, it is to be noted that the disease will persist in the population. The disease will attain its endemic status during initial 150 days and after that 1.511070942 million infected people will persist in the community if control measures were not implemented properly. In Figure 7c , we plotted sample simulations of Q(t) for different initial value of the state. All of the curves tend to stabilize itself during initial 300 days. A sudden increase could be observed at early stages and then a gradual decrease in quarantine population irrespective of the initial population. It is clear from the figure that as lim t→∞ Q(t) = 11.19220059 = Q . These figures suggest that the infection will persist in the Hubie province if the current situation remains the same. However, now Chinese data suggest that COVID-19 is well-controlled by Chinese government.

One of the very applied and effective analysis to design the control programme for infectious diseases is optimal control theory, see for detail [12, 20, 22, 23] . The optimal control theory will be utilized to design a control mechanism for COVID-19 on the basis of local sensitivity analysis. Here our main purposes are:

• Reducing vulnerability to COVID-19 by isolating susceptible population and limiting its movement to only non-dangerous zones. This control will be implemented to susceptible population and will be denoted by u 1 (t).

• Reducing COVID-19 infected people via treatment. Our strategy is to imposed this control on infected people and amount of this control at any time t will be represented by u 2 (t).

• Most importantly we will maximize the size of the class Q(t) in order to reduce the spread of the disease with minimum cost.

Here it could be noted that it would make sense if the isolation of the vulnerable population but the treatment of infected agents are taken. For this we modify system (1) by incorporating isolation u 1 (t) for S(t) and treatment u 2 (t) for I(t). So from this it is straight forward that model (1) becomes a special case control model if we let u 1 (t) = 0 (the case of no isolation) and u 2 (t) = 0 (the case of no treatment). There are various types of representations of the control variables in the corresponding system. Majority of the researchers used (1 − u(t))I(t) type of representation (like in [31] ), whereas, other use u(t)I(t) (like in [32] ) or r 0 u(t)I(t) (like in [1] ) etc. In the present study, the authors used the techniques of [32] . For minimization purpose, we formulate the objective functional of the form

subject to

In the objective functional (24), ξ 1 , ξ 2 and ξ 3 are the weights constants. Moreover ξ 4 and ξ 5 are also the positive constants, measures the relative cost on isolation and treatment respectively. Equation (24) has a clear goal: to decrease the number of people who are at risk as well as infected persons and maximizing the quarantined population. Many researchers assumed that the susceptible population has no concern with the objective functional. However, we believed that more susceptible individual means more people at risk. The longer the time that a vulnerable individual stays in the susceptible state, the longer he or she will be exposed to the disease and it is more likely that individual will be infected. By minimizing vulnerability to the infection ultimately lead towards elimination of the disease from the community. That is why we included the susceptible population in the objective functional. A population can be protected from disease in many ways. For example, the number of susceptible individuals can be reduced through immunizations, the contact rate can be reduced through quarantines or public health campaigns, or the removal rate can be increased through better medical treatment of the sick. In this study, we included quarantine as immunization strategy in the model. It is worthy to notice that quarantine individuals do not participate in the total active population, thus, it should be maximized in any control program. We wish to minimize the objective functional, so a term with negative sign will be ultimately maximized and the same was reflected in our objective functional in case of Q(t). Thus, our objective is to find (u * 1 , u * 2 ) such that

subject to the control system (25) and (26), where U, the control set is as

First we show the existence analysis of such controls.

We prove that there is a solution to problem (24)- (26) , which exists. So clearly the control functions are Lebesgue measurable and not negative, therefore with nonnegativity of initial data it is very handy to investigate that a solution with the axiom of bounded-ness and positivity exists, for example see for more detail [24] [25] [26] .

where

Clearly Eq. (29) shows that it is a nonlinear equation and the coefficient are bounded. Now we set

which satisfies

where N = max{(n 1 , n 2 , n 3 )} is free of state variables of (1). We also write

where M = max{N, A } < ∞, which implies that G is uniformly Lipschitz continuous. This from the control functions and non-negativity of state variables, it is easy to conclude existence of solution.

We show the optimal values of the control variable minimize the objective function in the following result.

There exist optimal values of the controls namely u * = (u * 1 , u * 2 ) ∈ U which is the solution of the control problem (24)- (26) .

Proof. It is already shown that the control and state functions are not negative. Further, the set U as stated by equation (28) that the controls functions are closed as well as convex. Moreover, the bounded-ness of system (25) shows that the necessary compactness holds. It is also to be noted that the state and control functions within objective functional indicated the convexity of objective functional. Hence there exist optimal controls (u * 1 , u * 2 ).

To show the characterization of the optimal solution to (24)-(26), we need to define the Lagrangian as well as Hamiltonian for the proposed control problem. Let us symbolize x = (S, I, Q) and u = (u 1 , u 2 ). Then the Lagrangian L is defined by

and the associated Hamiltonian function H becomes

where λ = (λ 1 , . . . , λ 3 ) and g(x, u) = (g 1 (x, u), . . . , g 3 (x, u)) with

Thus, the Hamiltonian function becomes

Next, to find the optimal solution, we use the standard Pontryagin's Maximum Principle [19, 20] , if (x * , u * ) is an optimal solution for the proposed control problem (24)- (26) , then there exists a non-trivial vector function λ such that the Hamiltonian system

the maximality condition

and the transversality condition

hold.

Theorem 5. Let S * , I * , and Q * be the optimal values of states variables associated to the optimal controls (u * 1 , u * 2 ) for problem (24)- (26) . Then, there exist adjoint variables λ i (t), i = 1, 2, 3, satisfying

with terminal (transversality) conditions

Moreover, the optimal control functions u * 1 (t) and u * 2 (t) are represented by

and

Proof. We founded that the optimal values of states and controls by solving analytically the optimal system, which includes the states (25) and the adjoint equations (37), boundary conditions (26) and (38), along with the characterization of optimal control variables (u * 1 , u * 2 ) in the form of (39) and (40).

In order to obtain the numerical solution for the control problem, we will use Runge-Kutta procedure of order four. More precisely, we descritize the original model (1), the control model (25), the adjoint system (37) and will coded the same in MATLAB along with subsidiary conditions (26), (38) and the characterization of the control (39)-(40). The parameter values used for simulating the control and uncontrolled problem is depicted in Table 3 . The simulation for the control and without control problems were carried out and a comparison was presented among real data, without control and with control curves in Figure 8 for all classes. Figure 8a shows the dynamics of susceptible population in three different perspective. If we compare the uncontrolled curve with real data, we reach to the conclusion that our model best fit the realistic scenario. Further, it is clear that without implication of the control strategies, more people will be vulnerable to the disease and eventually a particular portion of the population will be prone to the disease at all the time. However, if we imposed the control strategies upon the susceptible people, the population will tend to decline and after 100 days there will be no susceptible individual in the community.

Again, Figure 8b shows that our proposed model best explain the dynamics of COVID-19 in the Hubie province. Figure suggest that if we ignore the control policy, the infection will tend to grow and the decrease and eventually persist constantly in the Hubie province for long run. However, after 100 days there will be no infection if one impose the control strategies in true spirit. Similarly, the usefulness of the proposed model and the best output of the control upon quarantine population could be observed in Figure 8c . Thus, we can say that the model under consideration is correct up to the mark for explaining the dynamics of COVID-19 in Hubie province and to eliminate the disease, it is most beneficial and effective to adopt the same control strategies as outlined in the current study. Figures 9a and 9b shows the dynamics of the control variables u 1 (t) and u 2 (t). During initial phase of the control policy (up to 60 days), the control u 1 (t) have very low effect upon the susceptible people, afterward, the control tend to increase very fast and finally reach to its maximum value. Similarly, the control u 2 (t) is silent till initial 70 days and then start working with maximum returns. 

Keeping in view the ongoing disastrous spread of the novel coronavirus (COVID-19) across the globe, we have proposed a robust mathematical model that is able of assisting the health regulatory authorities in adopting important safety measures and management strategies for controlling such spreads. In order to test the validity and accuracy of the proposed model, in this paper, we have considered Hubei Province of the People's Republic of China, and is part of the Central China region. We have conducted a detailed analysis of our model, and applied it to study the Wuhan epidemic using publicly reported data. Parameters in the basic reproduction number R 0 of this model comes from three compartments; i.e., from the quarantine individuals, the infected individuals, and susceptible individuals. We have established the global asymptotic stability of the disease-free equilibrium when R 0 < 1, and the global asymptotic stability of the endemic equilibrium when R 0 > 1. Our numerical simulation results demonstrate the application of our model to the COVID-19 outbreak in Wuhan. Our model best fit the reported data very well. Through data fitting, we obtain an estimate of basic reproduction number, R 0 = 2.261469439. Our model predicts the appearance of an epidemic peak, after which the infection level would decrease and approach an endemic state in the long run. We also find that if we use constant transmission rates instead, the model would predict a much higher and unrealistic epidemic peak. This is caused by the fixed transmission rates that do not reflect the impact of on-going disease control measures. It is an indication that using epidemiologically and environmentally dependent transmission rates can potentially generate more practical simulation results. At present, many aspects regarding the pathology, ecology and epidemiology of the novel coronavirus remain unknown [30] , which add challenges to mathematical modeling. Particularly, in our current model, we have employed a bilinear incidence rate based on the law of mass action to represent the human-to-human transmission route. Though, a saturation based incidence rate might better characterize the environmental pathogen, and we hope to investigate it in our future modeling efforts. Given the current development of COVID-19, it is widely speculated that this disease would persist in the human world and become endemic. Our mathematical analysis and numerical simulation results support this speculation. The findings in this study imply that we should be prepare our self to fight the coronavirus infection for a much longer term than that of the current epidemic wave, in order to reduce the endemic burden and potentially eradicate the disease eventually. Among other intervention strategies, new vaccines for the novel coronavirus, which are currently in research and development, could play an important role in achieving that goal. Further, our study strongly support the fact if we imposed the two indicated control measures, we can easily eliminate the COVID-19 in almost 100 days from Wuhan. We emphasize that our data fitting is based on the reported confirmed cases in Wuhan from January 21 to March 2 in 2020. These confirmed cases were determined by the method of nucleic acid testing kits. On February 12, 2020, the national health commission in China started including cases confirmed by another method; i.e, clinical diagnosis, which refers to using CT imaging scans to diagnose patients. This change of criteria led to a surge of confirmed cases on February 12 (with an increase of about 14,000 new cases for Wuhan in a single day). We plan to address the new development of the outbreak data in another piece of work in the near future. We also plan to expand our modeling efforts to the province and country levels beyond the epicenter, the city of Wuhan, and study the spread of the novel coronavirus in larger spatial scales.

In a nutshell, we can say that understanding the dynamics of a fast-moving disease (like COVID-19) in the early stages of a pandemic is indeed a challenging task for the health official and policy makers. A single model is not appropriate, and long-term predictions of such models are likely to be futile. Perhaps, one can use model to short-term data which will enables us to comprehend deeply the existing data as well as to make predictions when data is unavailable. For this purpose, fast-moving models should be developed and analyzed for better understanding of such diseases. The strategy of short-term data fits will definitely help us to control the spread of COVID-19 at earliest. As infectious disease will be always with us, thus, the results of a model could be applied to any next disease if not working for the current situation.

In this article, we have established a model for the transmission dynamic of novel CoVID-2019 by taking into account the classification of different phases of its spread in population. First, we studied the statistics of novel CoVID-2019 in Hubie, China and then we combined different characteristics of the disease and formulated a mathematical model which describe the dynamics of the infection. Further, we presented different mathematical analysis on the proposed model including boundedness, biological feasibility, positivity and existence of equilibria of the proposed model. For analysis purpose, initially we obtained the basic reproduction number and then derived all the disease-free and endemic equilibrium of the model. The local and global stability of the disease-free equilibrium was carried out by using the linearization and Lyapnove functional theory, respectively. The DFE is globally asymptotically stable if R 0 < 1. Similarly, it is shown that the endemic equilibrium is locally, as well as globally asymptotically stable whenever R 0 > 1. The value of R 0 was calculated for the COVID-19 at Hubie and it was founded that the disease will persist for long term if not properly addressed. Finally, the numerical simulation and sensitivity analysis were presented to show the feasibility of the proposed model and effectiveness of the control strategies.

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Fenner's veterinary virology

Crossspecies transmission of the newly identified coronavirus 2019nCoV

WHO R&D Blueprint: informal consultation on prioritization of candidate therapeutic agents for use in novel coronavirus 2019 infection

Replication of MERS and SARS coronaviruses in bat cells offers insights to their ancestral origins

A familial cluster of pneumonia associated with the 2019 novel coronavirus indicating person-to-person transmission: a study of a family cluster

Clinical features of patients infected with 2019 novel coronavirus in Wuhan, China

Novel Coronavirus Disease 2019 (COVID-19): An Emerging Infectious Disease in the 21st Century

Clinical features of patients infected with 2019 novel coronavirus in Wuhan, China

Epidemiological and clinical characteristics of 99 cases of 2019 novel coronavirus pneumonia in Wuhan, China: a descriptive study

Characteristics of and public health responses to the coronavirus disease 2019 outbreak in China

Report of the WHO-China Joint Mission on Coronavirus Disease 2019 (COVID-19)

The transmission dynamic and optimal control of acute and chronic hepatitis B

The stability of dynamical systems

Further notes on the basic reproduction number

Modelling the spread of COVID-19 with new fractal-fractional operators: Can the lockdown save mankind before vaccination

Modeling the dynamics of novel coronavirus (COVID-19) with fractional derivative

Modelling the epidemiology of hepatitisB in NewZealand

Optimal HIV treatment by maximising immune response

Dynamic optimization: the calculus of variations and optimal control in economics and management

Optimal treatment of an SIR epidemic model with time delay

Differential equations, Mathematics in Science and Engineering

Viral dynamics and control of hepatitis B virus (HBV) using an epidemic model

Modeling and control of the hepatitis b virus spreading using an epidemic model

Optimal control strategy of SEIR endemic model with continuous age-structure in the exposed and infectious classes

Optimal Control Theory for Infinite Dimensional Systems

Differential Equations: Classical to Controlled, Mathematics in Science and Engineering

Infectious Diseases of Humans: Dynamics and Control

Optimal vaccination strategies for an SEIR model of infectious diseases with logistic growth

Applied methods of cost-benefit analysis in health care

Multiple transmission pathways and disease dynamics in a waterborne pathogen model

Asymptomatic carriers in transmission dynamics of dengue with control interventions

Optimal control strategies for dengue transmission in Pakistan

This work was supported by the National Natural Science Foundation of China (11971493) and the Fundamental Research Funds for the central Universities (No. 20lgpy137).