Mathematical analysis of COVID-19 by using SIR model with convex incidence rate.

This paper is about a new COVID-19 SIR model containing three classes; Susceptible S(t), Infected I(t), and Recovered R(t) with the Convex incidence rate. Firstly, we present the subject model in the form of differential equations. Secondly, "the disease-free and endemic equilibrium" is calculated for the model. Also, the basic reproduction number R 0 is derived for the model. Furthermore, the Global Stability is calculated using the Lyapunov Function construction, while the Local Stability is determined using the Jacobian matrix. The numerical simulation is calculated using the Non-Standard Finite Difference (NFDS) scheme. In the numerical simulation, we prove our model using the data from Pakistan. "Simulation" means how S(t), I(t), and R(t) protection, exposure, and death rates affect people with the elapse of time.

Introduction

In December 2019, a new kind of virus named “corona” was reported to badly affect the Chinese city of Wuhan. The said virus and its resultant outbreak hit the city of Wuhan first and later affected almost the whole world. It took hundreds of thousands of lives worldwide. It is hard to take a single point of view on this virus’s origin. It may be due to a seafood market exchange, or the people’s migration from one place to another, or the transmission from animals to humans; It may also be due to human-to-human interactions. So far, the virus has devastated almost everything around the world. Social life, health, economy, education - generally, each segment of human life has been severely affected. Health researchers, governmental policymakers, and health care authorities are puzzled in combating the deadly outbreak. They all have their point of view on the situation. They are trying hard to, at least, minimize the number of deaths caused by the outbreak. The people infected in the coronavirus pandemic experience mild respiratory problems. Fever, dry cough, throat infection, and fatigue are the symptoms of this disease. People may also have the symptoms as follows; nasal infection, aches, and sore throat. Mathematical modeling plays an important role in describing the epidemic of infectious diseases. The purpose of mathematical modeling is to represent different types of a real-world situation in mathematical language. A number of mathematical model are studies in the pervious literature [1], [2], [8], [10]. Also SARS-CoV-2 is study by many researchers in current research literature [9], [11], [12], [14]. We will study SARS-CoV-2 by developing SEIR model later on in this work. Recently many authors have established numerious models for COVID-19 under different concept of fractional calculus. In this regards very useful models have been established, we refer some as [20], [21], [22], [23], [24], [25]. To find out the different dynamics of a disease and therefore to overcome it at an early stage, mathematical modeling plays an important role there [26], [27], [28]. The area dedicated to the investigation of biological pandemic and also epidemic models for recent diseases SARS-CoV-2 of research. Numerous examples of mathematical models for this pandemic are found in the current study [16], [17]. To understand the stability theory, existence theory, and theory of reform SARS-CoV-2 [6], [7], [13], [19], can be model, and its outcome can be predicted. Plan of prevention is also possible. In addition, one can find a possible lock-down strategy. Especially impressed with the excellent features of the SEIR model using non-linear saturated incidence rates [3], [4], [5], [6], [8], [15].

Model formulation

In this section of manuscript, we formulate our new model for in the form of following system (1). We take whole population N(t) into three classes S(t), I(t) and R(t), which represent Susceptible, Infected and Recovered compartment in the form of differential equations given below (1),

For above system (1) is presented in the form of flow chart as.

In Table 1 , we describe parameters used in system (1). In system (1), add all equations, implies

Table 1

Physical Interpretation of parameters of the system.

Parameters	The physical Description
S(t)	Susceptible compartment
I(t)	Infected compartment
R(t)	Recovered compartment
d_0	Death due to corona
μ	Natural death
b	Birth rate
β	Protection rate
k	Constant rate
α	Isolation rate
γ	Recovery rate

Here represent whole population as . We get

With

If and only if

From the equation of the system (1), it show

Which implies that if , where . We can get

Equilibria

For the system (1), we suppose the existence of equilibrium. Disease free equilibrium is exist for some values of the variables used in (1), which is denoted by .

Endemic equilibria

Expression for The basic reproductive number

In epidemiology there is most important parameter, which give us idea about how the disease is flow in the whole population. From , we look how the disease id spread in population and we can control it from this. The method of finding is below let , then from system (1),whereand

Jacobian of isand Jacobian of is

Hence

We have

From this, we get is

To computes the basic reproduction number we obtained from the parameters used in Table 2 [18], which show that the COVID-19 that occurred in Pakistan is well controlled by Pakistan government. We have the following theorem on the basis of (3).

Table 2

Description of parameters and their values [18].

Parameters	Physical description	Numerical value
S(t)	Susceptible compartment	220 in millions
I(t)	Infected compartment	0 in million
R(t)	Recovered compartment	0 in million
d_0	Death due to corona	0.02
μ	Natural death	0.0062
b	Birth rate	10.7
β	Protection rate	0.009, 0.0009
k	Constant rate	0.00761
α	Isolation rate	0.009, 0.0009
γ	Recovery rate	0.0003

(i) If there is no positive equilibrium of system.

(ii) If there is a unique positive equilibrium of the model (1) , called the endemic equilibrium.

Local Stability

We reduced our model (1) for local stability. Furthermore, to obtained the result which show “disease free and endemic equilibrium”. For system (1). We reduced and get

Subject to initial condition

For local stability, we have the following theorem.

If , then the system (4) is locally asymptotically stable at the disease free equilibrium .’

At the jacobian matrix is given byThe auxiliary equation of is given bywhereWe haveThe Routh-Hurtwriz criteria is satisfied as and if . which show the system (1) is locally asymptotically stable at .

Furthermore, at the system (4) is locally asymptotically stability analogous to . We are going to prove it in the next theorem.

At , if then system (4) is locally asymptotically stable.

For system (4) jacobian matrix isAfter some operations on matrix , we getWe calculate trace and determinant of andThe determinant of . The real part at “endemic equilibrium” of model (4) has negative. Thus, with condition , we have that the endemic equilibrium of system (4) is locally asymptotically stable.

Global stability

Here, we present Global stability for the system (1). For “global stability of disease-free and endemic equilibrium”, we constructed a function known as Lyapunov function in the following theorem.

If then disease free equilibrium of the system (4) is globally asymptotically stable. Otherwise unstable.

To prove this, we construct a Lyapunov function as followingsuch that are constants. With respect to time t taking derivative of (8) with, we haveWe getLet assume , we get finallyHence “globally asymptotically stable” for system (1) with has reached. Further, We are going to prove a theorem for “global stability of the endemic equilibrium” of model (1).

The endemic equilibrium of model (1) is stable globally asymptotically if .

By constructing Lyapunov function, we prove the above resultTaking derivative with respect to time (9), we getPutting the values from (1) After some arrangement we getThus , the “endemic equilibrium” of the model (1) is “ globally asymptotically stable”, show that .

Numerical results and discussion

In this part of our manuscript, we calculated numerical simulation for model (1) with values used on Table 2. We take data from 1 February 2020 to September corresponding to different compartments involve in the system (1) from Pakistan. Here, we use (NSFD) Non-standard Finite Difference scheme [13], [15], [19] to rewrite the system is

Which is decomposed in Nonstandard Finite Difference scheme as

Just like above Eq. (11), we can write the system (1) in Non-Standard Finite Difference Scheme as

For the real data of Pakistan [18], we testified our model (1) taking the values of parameter of Table 2 from first February 2020 to 20th of September 2020. From Fig. 1, Fig. 2, Fig. 3 , we see that as the susceptibility was decreasing the infection was increasing in first four months but in the month of July and August the infection rate became slow and finally in the last it was nearly become stable. Also during this time the recovery rate was rapid from the infection as in Fig. 3. The concerned simulation was performed for taking the protection parameter . Now by deceasing the protection and isolation rate further up to . We plot the results in the given Fig. 4, Fig. 5, Fig. 6 . We see that the infection rate became slow on reducing the protection and isolation rate. Therefore the recovery is also become slow. From these simulation we observed that protection and isolation rate play significant roles in controlling the infection from further spreading in the community. see Fig. 2,.

Fig. 1

Dynamical behavior in of susceptible population of the considered model.

Fig. 2

Dynamical behavior of infected population of the considered model.

Fig. 3

Dynamical behavior of recovered population of the considered model.

Fig. 4

Dynamical behavior in of susceptible population of the considered model.

Fig. 5

Dynamical behavior of infected population of the considered model.

Fig. 6

Dynamical behavior of recovered population of the considered model.

Conclusion

Conclusion of the numerical results shows the projection of model (1). The output derived from the NCOVID-19 display convex incidence rate. The current manuscript declared the high contiguous rate from infected population to susceptible population. To overcome the pandemic the migration should be strictly prohibited for the sake of saving humanity. Also the immigration of exposed population to infected community increased the infection. Isolation of infected one is the best option to secured the healthy community. It is necessary to judge the spread and model with various parameters for proper supervision. The proper treatment of this pandemic is to keep infected away from healthy people. High internal defense system aids to get healthy soon while the low internal defense system need more attention. This is the only solution to overcome recent outbreak within a short period. The current discussion demonstrate the quick transfer of NCOVID-19. The COVID-19 shared the same properties like SARS having mortality rate of 2 percent. There is no vaccine available in the current time but to isolate was the best option. Also social distancing is the best way to control this deadly various.

Funding

There is no source to support this article financially.

Availability of Data

This is not applicable in this research work.

Declaration of Competing Interest

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.