Mathematical modelling of the COVID-19 pandemic with demographic effects

In this paper, a latent infection susceptible–exposed–infectious–recovered model with demographic effects is used to understand the dynamics of the COVID-19 pandemics. We calculate the basic reproduction number (R0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${R}_{0}$$\end{document}) by solving the differential equations of the model and also using next-generation matrix method. We also prove the global stability of the model using the Lyapunov method. We showed that when the R0<1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${R}_{0}<1$$\end{document} or R0≤1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${R}_{0}\le 1$$\end{document} and R0>1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${R}_{0}>1$$\end{document} or R0≥1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${R}_{0}\ge 1$$\end{document} the disease-free and endemic equilibria asymptotic stability exist theoretically. We provide numerical simulations to demonstrate the detrimental impact of the direct and latent infections for the COVID-19 pandemic.

treatments like antibiotic drugs are used, and in severe cases, supportive treatment like a breathing machine is given to protect vital functions of the organs. The virus infects all ages of humans, but the higher risk is more on adult individuals with severe illness relating to respiratory diseases, organ diseases, and blood diseases [7]. The basic reproduction number (R 0 ) is a critical threshold quantity associated with viral transmissibility, and it has been used to understand the transmission of the COVID-19. Epidemiological R 0 is a value used to describe the contagiousness of the pathogen, and it is estimated using incidence data during the first phase of a disease outbreak. It describes the number of people on average that would be infected from a case introduced into a population. The initial COVID-19 pandemic R 0 , according to the WHO, was estimated to be between ranges of 1.4 and 2.5 [6]. That is, one infected person will infect an average of 2 persons in his/her lifetime. In the first phase of the epidemic, Zhao et al. [12] estimated the average R 0 for COVID-19, from 3.3 to 5.5, and Read et al. [13] estimated to range between 3.6 and 4.0.
Stability analysis, which has a direct relationship with R 0 , is also another way to understand infectious disease. It is believed that when R 0 is above unity, the disease will persist, and the stability is endemic, and when R 0 is less than unity, the disease will die out, and the stability is disease-free. The analysis is done by partitioning the state of individuals in the population into different compartments. For instance, since COVID-19 has an incubation period, the population can be divided into those who are capable of being affected by the virus, call the susceptible (S) compartment. When a visibly infected (I, i.e. a person confirmed to have the virus) individual is identified, from the S compartment; the infectious person contact and contact-contact form an exposed (E) compartment; and those overcoming the illness of the virus and get well form the recovered (R) compartment. The SEIR is interpreted using differential equations, where differential equation techniques and simple algebraic methods are used to study the dynamic of the disease. However, the R 0 can also be calculated using the differential equations model at the state when the disease is free from the population (disease-free state).
In this study, a deterministic four-compartment SEIR model is considered to inspect the stability analysis of the COVID-19 pandemic using differential equation techniques. That is, contrary from traditional SEIR model, where an individual in the E compartment is infected but not infectious, we consider E as another infection transition point but not visibly infectious [14,15]. This is done by formulating four nonlinear differential equations and provides theoretical and numerical analysis of the model. Our results show that, theoretically, the disease-free and endemic equilibria of the model are locally and globally asymptotically stable and the direct and the rate of infection transmission from an individual after exposure to the virus are detrimental for the COVID-19 pandemic.

Model framework
In this section, we describe an epidemic transmission SEIR model with demographic changes. The model is used in epidemiology to compute the amount of susceptible, exposed, visibly infected, recovered people in a population (N). Since the asymptomatic and pre-symptomatic people can transmit the virus but their symptoms are not visible, they are grouped into the E compartment and infection from E is referred to as latent infection. This model is used under the following assumptions: • The population is constant but large. • The only way a person can leave the susceptible state (S) is to become infected either from the exposed (E) or from visibly infected (I) state or die of natural causes. • The only way a person can leave the E state is to show signs and symptoms of the illness or die of natural death. • The only way a person can leave the I state is to recover from the disease or die from natural death or die as a result of the disease. • A person who recovered (R) from the illness received permanent immunity. • Age, sex, social status, and race do not affect the probability of being infected. • The member of the population has the same contacts with one another equally. • All births are into the susceptible state, and it is assumed that the birth and natural death rates are equal.
The transmission is measured at Sβ(I + κE)/N , where β is the direct transmission rate, and κ is the proportional rate constant when an uninfected individual comes into contact with an individual from state E. We assume natural birth and death rate to be measured at an equal rate µ and induced death rate measured at δ . The rate for an individual to move from state E to state I is measured at rate σ , , and the rate of recovery is measured at γ . Figure 1 represents the latent infection SEIR model, which is described using the system of nonlinear ordinary differential equations where S(t) = S, E(t) = E, I(t) = I and R(t) = R denote the number of susceptible, exposed, infectious, and remove individuals at time t , respectively, and (1) is subjected to the initial condition For simplicity system (1) is reduced to a proportional framework given as

Positivity of the solution
Assume that system (4) has a global solution corresponding to non-negative initial conditions. Then, the following Lemma confirms that the solution is non-negative at all times. Proof. We use the contradiction: we assuming there exists positive real t 1 , t 2 and t 3 for which one of the conditions hold: . Thus, the solutions of s(t), e(t)andi(t) remain positive for all t > 0. Hence, the positively invariant for the system (4) is The equilibrium points and reproduction number calculations of the model There are two equilibrium points for the system (4), i.e. the disease-free equilibrium (DFE), the state when the disease is absent, and the endemic equilibrium (EE), which is the state when the disease continues to persist in the population. Let the DFE points of the model are denoted as E 0 = s 0 , e 0 , i 0 and represent a system (4) at E 0 as In terms of i 0 , from the last equation of (6), we get Adding the first two equations of (6) and substituting for e 0 , we get Because at the disease-free state no one has the infection then, i 0 = e 0 = 0 . We can see that E 0 = s 0 , e 0 , i 0 = (1, 0, 0).
Also, the EE points are denoted asE * = (s * , e * , i * ) , where s * , e * and i * are calculated by letting s 0 = s * , e 0 = e * , i 0 = i * , and then, the second equation of (6) becomes Multiplying both sides of (7) by σ/(µ + σ )(µ + γ + δ) , we get which implies that Hence also (6) µ − βs 0 i 0 + κe 0 − µs 0 = 0, To understand the stability of the model, we need an expression to estimate the basic reproduction number ( R 0 ). However, from (6) lets 0 = s, e 0 = e, i 0 = i , in terms of i the second equation of (6) becomes By the factorising method, we get either i = 0 or Equating the highest i value from (9) to zero, we get Since the threshold for R 0 is unity, we then assume Equation (11) is justified using the next-generation matrix method defined in [16] as is the spectral radius of the matrix F V −1 and the largest eigenvalue of K is the R 0 . F and V are the matrices associated with the DFE points defined as The largest eigenvalues of K is as (11) given as Hence
Proof The Jacobian matrix of system (3) associated with DFE is given as with characteristic polynomial where ⋋ is an eigenvalue. It is easy to see that for Theorem 1 to satisfy and The proof of Theorem 1 is complete.

the DFE is globally asymptotically stable in Ω.
Proof To prove the global asymptotic stability (GAS) of the DFE, we construct the following Lyapunov function V : Ω → R , where V (s, e, i) = i(t) . Then, the time derivative of V is given as since at the equilibrium points di 0 /dt = de 0 /dt = 0 . Therefore Substituting s 0 and e 0 into (13), we get

Stability analysis of the endemic equilibrium
Theorem 3 If R 0 > 1 , the endemic equilibrium is locally asymptotically stable.
Proof To prove the LAS of the endemic equilibrium, we consider the Jacobian matrix associated with E * , that is Substituting for s * , e * and i * , we get if is an eigenvalue, then where I is a three-dimensional identity matrix which by matrix simplification method we then get where From R 0 , we get hence using (18) we get and Since a > 0, b > 0, c > 0, and ab − c > 0 , according to the Routh-Hurwitz criterion, the endemic equilibrium of system (4) is LAS.
Firstly, we investigate the DFE by assuming β = 0.0533 ; we observe that when R 0 = 0.163 in Fig. 2a, Theorems 1 and 3 are satisfied for the DFE to be asymptotically stable. It is observed that when states E and I are decreasing, the susceptible population approaches unity with increasing time. Also, in Fig. 2b, we observe that when β = 0.533 , and R 0 = 1.63 , Theorems 2 and 4 for the EE of the model to be asymptotically stable are satisfied. It is observed that increasing the I proportion the S proportion declines until at a certain point in time when the I proportion started to decrease and the S proportion then increases. The decrease in the trajectories in the case of the E state is the result of the increase in asymptomatic individual to the visibly infectious state and natural death, whereas for the I state is the results of the increase in the recovered individual and those who might have died of natural or virus death.
From a mathematical point of view, it is easy to see that the EE tends to DFE which is dependent on the decreasing rate of κ and β . We investigate the effect of the direct and latent infection rate numerically by keeping the EE parameter values constant and regulating the degree of κ and β . That is, the positive effect for the direct and latent infection  parameters is noticeable when the magnitude of κ and β decreases, and their output is similar to that of Fig. 2a. It is observed in Fig. 3 that as κ decreases the endemic trajectory patterns are similar to the DFE curves in Fig. 2a. Also, we observe that when β is lower in magnitude, lesser susceptible individuals become infected as the curve tends to increase in proportion. In Figs. 3 and 4, it is easy to see that direct and latent infection transmissions can enhance the persistence of the COVID-19 pandemic.

Conclusion
In this paper, we formulate a latent infection SEIR model to investigate the stability analysis of the COVID-19 pandemic with demographic effects. We use differential equation techniques and simple algebraic procedures to describe the dynamics of the model theoretically. We showed that the model has two equilibrium states, which are disease-free and endemic equilibrium. The stability analyses show that the two equilibria states are locally and globally asymptotically stable theoretically, which are confirmed numerically using epidemiological data of COVID-19 pandemic. From our study, we observe that when κ and β decrease, the infected population also decreases. The biological implication of this is that the direct and latent infections are detrimental to the COVID-19 pandemic. Therefore, isolating exposed and visibly infected individual is an important strategy in controlling the COVID-19 pandemic.