# Mathematical Approach to Investigate Stress due to Control Measures to Curb COVID-19

Authors: James Nicodemus Paul; Silas Steven Mirau; Isambi Sailon Mbalawata
Journal: Computational and Mathematical Methods in Medicine
 (2022)

---

## Abstract

COVID-19 is a world pandemic that has affected and continues to affect the social lives of people. Due to its social and economic impact, different countries imposed preventive measures that are aimed at reducing the transmission of the disease. Such control measures include physical distancing, quarantine, hand-washing, travel and boarder restrictions, lockdown, and the use of hand sanitizers. Quarantine, out of the aforementioned control measures, is considered to be more stressful for people to manage. When people are stressed, their body immunity becomes weak, which leads to multiplying of coronavirus within the body. Therefore, a mathematical model consisting of six compartments, Susceptible-Exposed-Quarantine-Infectious-Hospitalized-Recovered (SEQIHR) was developed, aimed at showing the impact of stress on the transmission of COVID-19 disease. From the model formulated, the positivity, bounded region, existence, uniqueness of the solution, the model existence of free and endemic equilibrium points, and local and global stability were theoretically proved. The basic reproduction number (R0) was derived by using the next-generation matrix method, which shows that, when R0<1, the disease-free equilibrium is globally asymptotically stable whereas when R0>1 the endemic equilibrium is globally asymptotically stable. Moreover, the Partial Rank Correlation Coefficient (PRCC) method was used to study the correlation between model parameters and R0. Numerically, the SEQIHR model was solved by using the Rung-Kutta fourth-order method, while the least square method was used for parameter identifiability. Furthermore, graphical presentation revealed that when the mental health of an individual is good, the body immunity becomes strong and hence minimizes the infection. Conclusively, the control parameters have a significant impact in reducing the transmission of COVID-19.


---

## Body

## 1. Introduction

Coronavirus disease-2019 (COVID-19) is an infectious disease caused by a newly discovered coronavirus named severe acute respiratory syndrome coronavirus-2 (SARS-COV-2). This pandemic originated in Wuhan, China, with the first case reported in December 2019, and has spread to other parts of the world in early 2020 as discussed in [1, 2]. When the total confirmed cases globally were 125,260 and 4613 deaths in 24 hours, the World Health Organization (WHO), on 12th March 2020, announced COVID-19 disease as the World pandemic as presented in [3]. The global leaders were greatly bothered by this disease due to its fast spread from one person to another and its social and economic impact on their respective countries. The WHO and country leaders focused on finding ways to reduce the transmission of the disease by introducing some measures such as lockdown, quarantine, closing borders, travel bans, and isolation centers [4–7].Mathematical models are essential tools in evaluating various transmission and control intervention programs for infectious diseases. There are a number of mathematical models on COVID-19 pandemic developed from the start of this human disturbing disease [8]. The author discusses social isolation measures taken by the government in Brazil to fight COVID-19 disease, and also, the protection of health workers was discussed by Masandawa et al. [9]. In their study, Khan et al. [10] discussed on isolation and quarantine as the best ways of fighting COVID-19 disease. At the beginning of the pandemic, many ways to fight the disease were introduced by many scientists. Ullah and Khan [11] explored several ways such as social distancing, self-isolation, quarantine, and hospitalization and concluded that these are the best ways to curb COVID-19. Dos Santos [12], in his review, studied on the effects of social distancing and social isolation measures as good ways to fight against the disease, although they lead to social stress which influences the spread of the disease. Early 2020, some mathematical models for COVID-19 were developed and published. Dhanwant and Ramanathan [13] explored how the social distancing in India during the pandemic helped to reduce the transmission of COVID-19 disease because the spread of the disease is by social contact. Khan and Atangana [14] discussed the interaction of bats and unknown hosts and then the interaction of people in the seafood market where there is enough source of infection. They presented their results graphically to show how they can minimize the infections.Ashcroft et al. [15] discussed the impact of quarantine on the transmission of COVID-19 disease showing that many countries impose quarantine to ensure the exposed people and those from abroad are isolated for a specific period of time to prevent the spread of the disease. However, these measures took a larger economical toll and affected the health of the isolated individuals. Prati [16], from Italy, discussed the national impacts of quarantine psychologically. In their online survey of 1569 people living in Italy, they found that there are psychosocial factors that influence the disease, such as media exposure to COVID-19 outbreak, financial loss, higher worry, and negative attitude towards quarantine leading to psychological impacts. The mental health of public and healthcare professionals is affected by the pandemic, especially during quarantine time where hypervigilance arises because of fear and anxiety [17, 18].Stress affects many quarantined people whereby their immune system is disturbed by the COVID-19, and this is most likely because during quarantine, people are isolated from their families and community members, so they develop fear, and later, the body becomes stressed which affects their immune system. When the immune system is disturbed, it fails to fight against the intruders, which leads to the fast spread of COVID-19 throughout the body. Therefore, this research is aimed at formulating a deterministic model to explore the impact of stress in quarantine to the human population. The model has six compartmental classes (Susceptible, Exposed, Quarantine, Infectious, Hospitalized, and Recovered). The model is extended from the model given in [19] by incorporating the hospitalized class and introducing a stress parameter in a quarantined and infected class.The introduction of this work presented in Section1. Model formulation, discussion of its compartments, and parameters are presented in Section 2. Section 3 contains the discussion of the model analysis theoretically, which includes positivity and bounded regions, the existence, uniqueness of the model, reproduction number, and local and global stability of the COVID-19 disease. Section 4 deals with a discussion of numerical simulation for the model, including sensitivity analysis, numerical solutions, PRCC results, parameter identifiability, and model fitting by the least square method are presented. Section 5 concludes this work and contains the possible extension of this model.


## 2. Model Formulation

In this study, a mathematical model for COVID-19 was formulated based on realistic assumptions. The total population modelNt is divided into six human subclasses, namely, susceptible St (those who are at risk to contact COVID-19 infection), exposed Et (the population which is infected but not infectious), quarantined Qt (those who contacted a COVID-19-infected individual but did not develop any symptom), infectious It (those who have COVID-19 symptoms and are capable of spreading the disease), hospitalized Ht (infectious individuals admitted to a healthcare facility (active cases)), and recovered Rt (those recovered from the COVID-19).


### 2.1. Model Assumptions

By presenting an infectious disease with a mathematical model, the following assumptions for theSEQIHR model are considered based on the characteristics of COVID-19 disease in this work.
(i)
All members of the population can have an equal chance of getting a COVID-19 disease(ii)
Stress in a quarantined class is higher than that in an infectious class(iii)
Population from outside the country were taken directly to quarantine class(iv)
All compartments have an equal natural death rate(v)
The death due to the disease may be only in two variables (infectious and hospitalized)(vi)
Recruitment rates (Newborns) are assumed to be susceptible(vii)
Individuals are equally likely to be infectious to the infected individuals when coming into contact(viii)
Infected individuals are identified early and isolated (hospitalized) immediately for treatment(ix)
There are only two options, a patient will either recover or die, which means no treatment failure(x)
The population differs within a given time step where recruitment and leaving rates differ


### 2.2. Model Compartments and Dynamics

From Section2.1, the variable and parameter descriptions are to be presented by the following SEQIHR model compartmental diagram:LetN be the total population divided; the transmission dynamics of COVID-19 disease in a population are shown in Figure 1. Consider Table 1, showing the variables and their descriptions:Figure 1
Schematic flow diagram showing dynamics of COVID-19.Table 1
Model variable description.
VariablesDescriptionStNumber of susceptible population at timetQtNumber of quarantined population at timetEtNumber of exposed population at timetItNumber of infectious population at timetHtNumber of hospitalized population at timetRtNumber of recovered population at timetThe total populationNt is given by the mathematical equation:
(1)Nt=St+Qt+Et+It+Ht+Rt,where t∈0,t and t>0.Table2 shows the model parameters and their descriptions:Table 2
Model parameter description.
ParametersDescriptionβContact rate (effective transmission rate)μRecruitment rate to the susceptible populationbHuman natural death rateϕQuarantined population from infected countriesρRecovered population rate back to susceptible classαPopulation rate after the quarantined period to susceptibleνRecovery rate from infected populationλRecovered rate from hospitalized populationωProportion of exposed population with contradicting symptomsqProgression rate from exposed to hospitalized and infectious classesθProportion of exposed population with no symptomsη1Stress to the infected population from quarantineη2Stress to the hospitalized population from infectious classδDeath due to the disease from infectious and hospitalized classesThe model parameters found in Equation (2) are described in Table 2. By considering the SEQIHR model with six compartments in Figure 1, the following are the transmission phases:The Susceptible class,St, increases by the addition of a recruitment rate, μ. It also decreases by infection, if contacted with an infected individual at the rate of β and natural death at the rate of b.An individual enters the Exposed class,Et, after direct contact with an infected person, with an infection rate β. Furthermore, taking to hospital decreases the rate of ω and q for individuals with multiple symptoms. Individuals with clear (direct) symptoms are presented by (1-ω), and those with no symptoms are represented by θ, and then, the population in Et class is diminished by a leaving rate of natural death b.The Quarantined class,Qt, was formed by individuals’ progress from exposed, θ, and population rate from outside the country, ϕ, also decreased by individuals with no symptoms after being tested. Then, the individuals with no infections going back to the susceptible class are presented by α and those infected individuals after being tested, going to the infectious class at the rate of η1 and diminishing by leaving a natural death rate of b.The Infectious class,It, are individuals who progress from exposed at the rate of (1-ω) and then quarantine at the rate of η1. Additionally, this class decreases with hospitalized individuals at the rate of η2 and recovered at the rate of ν. Not only that but also it diminished by leaving the rate of natural death, b, and death due to the disease, δ.The Hospitalized class,Ht, are individuals who progress from the exposed class at the rate of ω and q, then from the infectious compartment at the rate of η2. Individuals recovered from hospitalized class at the rate of λ and diminished by leaving the rate of natural death, b, and death due to COVID-19 disease, δ.The Recovered class,Rt, are individuals who progress from a hospitalized compartment at the rate of λ and infectious rate ν and decrease by recovered individuals who are going back to the susceptible population ρ and then diminished by the leaving rate of b.


### 2.3. Model Equations

Based on the assumptions made and the relationship that exists between the variables shown in Figure1, the system of six ordinary differential equations is formed as in
(2)dSdt=μ+αQ+ρR−βIS−bS,(3)dEdt=βIS−θE−qE−bE,(4)dQdt=ϕ+θE−αQ−η1Q−bQ,(5)dIdt=1−ωqE+η1Q−η2I−νI−δI−bI,(6)dHdt=ωqE+η2I−λH−δH−bH,(7)dRdt=λH+νI−ρR−bR.


## 2.1. Model Assumptions

By presenting an infectious disease with a mathematical model, the following assumptions for theSEQIHR model are considered based on the characteristics of COVID-19 disease in this work.
(i)
All members of the population can have an equal chance of getting a COVID-19 disease(ii)
Stress in a quarantined class is higher than that in an infectious class(iii)
Population from outside the country were taken directly to quarantine class(iv)
All compartments have an equal natural death rate(v)
The death due to the disease may be only in two variables (infectious and hospitalized)(vi)
Recruitment rates (Newborns) are assumed to be susceptible(vii)
Individuals are equally likely to be infectious to the infected individuals when coming into contact(viii)
Infected individuals are identified early and isolated (hospitalized) immediately for treatment(ix)
There are only two options, a patient will either recover or die, which means no treatment failure(x)
The population differs within a given time step where recruitment and leaving rates differ


## 2.2. Model Compartments and Dynamics

From Section2.1, the variable and parameter descriptions are to be presented by the following SEQIHR model compartmental diagram:LetN be the total population divided; the transmission dynamics of COVID-19 disease in a population are shown in Figure 1. Consider Table 1, showing the variables and their descriptions:Figure 1
Schematic flow diagram showing dynamics of COVID-19.Table 1
Model variable description.
VariablesDescriptionStNumber of susceptible population at timetQtNumber of quarantined population at timetEtNumber of exposed population at timetItNumber of infectious population at timetHtNumber of hospitalized population at timetRtNumber of recovered population at timetThe total populationNt is given by the mathematical equation:
(1)Nt=St+Qt+Et+It+Ht+Rt,where t∈0,t and t>0.Table2 shows the model parameters and their descriptions:Table 2
Model parameter description.
ParametersDescriptionβContact rate (effective transmission rate)μRecruitment rate to the susceptible populationbHuman natural death rateϕQuarantined population from infected countriesρRecovered population rate back to susceptible classαPopulation rate after the quarantined period to susceptibleνRecovery rate from infected populationλRecovered rate from hospitalized populationωProportion of exposed population with contradicting symptomsqProgression rate from exposed to hospitalized and infectious classesθProportion of exposed population with no symptomsη1Stress to the infected population from quarantineη2Stress to the hospitalized population from infectious classδDeath due to the disease from infectious and hospitalized classesThe model parameters found in Equation (2) are described in Table 2. By considering the SEQIHR model with six compartments in Figure 1, the following are the transmission phases:The Susceptible class,St, increases by the addition of a recruitment rate, μ. It also decreases by infection, if contacted with an infected individual at the rate of β and natural death at the rate of b.An individual enters the Exposed class,Et, after direct contact with an infected person, with an infection rate β. Furthermore, taking to hospital decreases the rate of ω and q for individuals with multiple symptoms. Individuals with clear (direct) symptoms are presented by (1-ω), and those with no symptoms are represented by θ, and then, the population in Et class is diminished by a leaving rate of natural death b.The Quarantined class,Qt, was formed by individuals’ progress from exposed, θ, and population rate from outside the country, ϕ, also decreased by individuals with no symptoms after being tested. Then, the individuals with no infections going back to the susceptible class are presented by α and those infected individuals after being tested, going to the infectious class at the rate of η1 and diminishing by leaving a natural death rate of b.The Infectious class,It, are individuals who progress from exposed at the rate of (1-ω) and then quarantine at the rate of η1. Additionally, this class decreases with hospitalized individuals at the rate of η2 and recovered at the rate of ν. Not only that but also it diminished by leaving the rate of natural death, b, and death due to the disease, δ.The Hospitalized class,Ht, are individuals who progress from the exposed class at the rate of ω and q, then from the infectious compartment at the rate of η2. Individuals recovered from hospitalized class at the rate of λ and diminished by leaving the rate of natural death, b, and death due to COVID-19 disease, δ.The Recovered class,Rt, are individuals who progress from a hospitalized compartment at the rate of λ and infectious rate ν and decrease by recovered individuals who are going back to the susceptible population ρ and then diminished by the leaving rate of b.


## 2.3. Model Equations

Based on the assumptions made and the relationship that exists between the variables shown in Figure1, the system of six ordinary differential equations is formed as in
(2)dSdt=μ+αQ+ρR−βIS−bS,(3)dEdt=βIS−θE−qE−bE,(4)dQdt=ϕ+θE−αQ−η1Q−bQ,(5)dIdt=1−ωqE+η1Q−η2I−νI−δI−bI,(6)dHdt=ωqE+η2I−λH−δH−bH,(7)dRdt=λH+νI−ρR−bR.


## 3. Model Analysis

In this section, positivity, boundedness, derived equilibrium states, basic reproduction number, and stability analysis are discussed.


### 3.1. Positivity of the Model

For the model equations to be epidemiologically, we need to prove that the state variables are nonnegative∀t≥0.Theorem 1.
Let the initial data set beS,E,Q,I,H,R0>0. The solution set of the model system 2 is positive ∀t>0.Proof.
From the system of model Equation (2), consider the first equation:
(8)dSdt=μ+αQ+ρR−βI+bS.
By considering the negative term, by ignoring the rest, Equation (8) is reduced to
(9)dSdt≥−βI+bS.
This is the first-order linear differential inequality which can be solved by a separable methody′=fxgy (where S≥0) resulting in
(10)∫S0StdSS≥−∫0tβI+bdt,St≥S0e−βI+bt,in the absence of COVID-19 disease,
(11)St≥S0e−bt,S0≥0,∀t>0.By applying the same procedure to the remaining equations, the results areE0≥0, Q0≥0, I0≥0, H0≥0, and R0≥0. Therefore, the set of solutions St,Et,Qt,It,Ht, and Rt of the model is positive ∀t>0.


### 3.2. Invariant (Boundedness) Region

TheSEQIHR model is represented by differential equations in system 2, which is to be analyzed in a feasible region Ω, and all state variables and parameters of the model are assumed to be positive ∀t≥0. The bounded region is obtained through the following theorem.Theorem 2.
The setΩ is positively invariant and attracts all solutions in ℝ+6.Proof.
SinceNt=St+Qt+Et+It+Ht+Rt, then the derivative of Nt is given as
(12)dNdt=dSdt+dEdt+dQdt+dIdt+dHdt+dRdt.
Substituting model Equation (2) to Equation (12) gives
(13)dNdt=μ+αQ+ρR−βIS−bS+βIS−θE−qE−bE+ϕ+θE−αQ−η1Q−bQ+1−ωqE+η1Q−η2I−νI−δI−bI+ωqE+η2I−λH−δH−bH+λH+νI−ρR−bR.
Further simplification leads to(14)dNdt=μ+ϕ−S+E+Q+I+H+Rb−δI+δR,implying that
(15)dNdt=μ+ϕ−Nb−δI+δR.
For disease-free,δI⟹0 and δR⟹0, then dN/dt≤μ+ϕ−Nb.
This is a first-order linear differential equation, by a separable method.(16)∫N0NtdNμ+ϕ−bN≤∫0tdt,⟹Nt≥μ+ϕb1−e−bt+N0e−bt.Whent=0, N0≥0. When t⟹∞ then, N∞≤μ+ϕ/b.The Invariant region is given byΩ=S,E,Q,I,H,R∈ℝ+6:0≤Nt≤μ+ϕ/b. The SEQIHR model is biologically and epidemiologically meaningful; thus, we can consider the flow generated by the model for analysis.


### 3.3. Existence and Uniqueness of the Solution

From the first-order differential equation given in the form:y′=ft,y,yt0=y0. The following questions will be of interest:
(i)
Under what conditions can we say that a solution to the equationy′ exists?(ii)
Under what conditions does a unique solution exist to the equationy′?Consider the following equations to answer the question.(17)f1=μ+αQ+δR−βIS−bS,(18)f2=βIS−θE−qE−bE,(19)f3=ϕ+θE−α+η1+bQ,(20)f4=1−ωqE+η1Q−η2+ν+δ+bI,(21)f5=ωqE+η2I−λ+δ+bH,(22)f6=λH+νI−ρ+bR.


#### 3.3.1. Uniqueness of Solution

Theorem 3.
Let us useD to denote the domain:
(23)t−t0≤a,y−y0≤b,y=y1,y2,⋯,yn,y0=y10,y20,⋯,yn0.Proof.
Supposeft,y satisfies the Lipschitz condition; therefore,
(24)ft,y2−ft,y1≤kx2−x1.
Whenever the pointst,x1 and t,x2 belong to the domain D and k is used to represent the positive constant, then, there exist a constant α>0 and a unique solution yt of system 11 in the interval t−t0≤α. It is essential to note that condition (24) is satisfied by ∂fi/∂yj,i,j=1,2,3,⋯,n to be continuous and bounded in domain D.Ifft,y has a continuous partial derivative ∂fi/∂yj on a bounded closed convex domain ℝ (i.e., the convex set of real numbers), where ℝ is used to denote real numbers; then it satisfies a Lipschitz condition in ℝ. Our interest is in the domain:
(25)1≤ε≤ℝ.Therefore, we look for the bounded solution of the form0<ℝ<∞.


#### 3.3.2. Existence of a Solution

Theorem 4.
Let the domain be denoted byD, also defined in Equation (23), such that Equations (24) and (25) hold. Then, the existing solution of Equation (2) is bounded in the domain D.Proof.
From Equations (10)–(15), by showing that ∂fi/∂yj,i,j=1,2,3,4,5,6, then, Equations (1) and (2) are continuous and bounded. That means the partial derivatives are continuous and bound. Consider the exploration of the partial derivatives for all model equations.
From Equation (17), we obtain the following system of equations:
(26)∂f1∂S=−βI−b,∂f1∂S=−βI−b<∞,∂f1∂E=0,∂f1∂E=0<∞,∂f1∂Q=α,∂f1∂Q=α<∞,∂f1∂I=−βS,∂f1∂I=−βS<∞,∂f1∂H=0,∂f1∂H=0<∞,∂f1∂R=ρ,∂f1∂R=ρ<∞.
Similarly, from Equation (18), we obtain the following system of equations:
(27)∂f2∂S=βI,∂f2∂S=βI<∞,∂f2∂E=−θ+q+b,∂f2∂E=−θ+q+b<∞,∂f2∂Q=0,δf2∂Q=0<∞,∂f2∂I=βS,∂f2∂I=βS<∞,∂f2∂H=0,∂f2∂H=0<∞,∂f2∂R=0,∂f2∂R=0<∞.The same procedures are taken for Equations (19), (20), (21), and (22). Therefore, all partial derivatives are continuous and bounded; hence, from Theorem 4, it is concluded that there exists a unique solution of the model in Equation (2) in the domain region D.


### 3.4. Existence of Disease-Free Equilibrium Point (DFE)

The disease-free equilibrium point, obtained when the infected components are zero, can be done by setting the right-hand side of the equation equal to zero, as in(28)dSdt=dEdt=dQdt=dIdt=dHdt=dRdt=0.When there is no disease, then,E=0, Q=0, I=0, H=0, and R=0.By considering each model equation from Equation (2),
(29)μ+αQ+ρR−βIS−bS=0givesμ−bS=0,then,S=μb.In addition to the second model equation,(30)βIS=θ+q+bE=0givesE=βISθ+q+b,butβ=0,then,E=0.Similarly, for the third model equation,(31)ϕ+θE−αQ−η1Q−bQ=0givesQ=ϕ+θEα+η1+b,but,E=0,then,
(32)Q=ϕα+η1+b.By considering the same procedures for the fourth, fifth, and sixth model equations, the following is obtained:(33)I=0,H=0,R=0.Therefore,(34)E0=S0,E0,Q0,I0,H0,R0=μb,0,ϕα+η1+b,0,0,0T.Equation (34) represents the state in which there is no infection and is known as the disease-free equilibrium point.


### 3.5. Basic Reproduction Number (R0)

The basic reproduction numberR0 is the midpoint number of infections caused by an infectious individual during the entire period of infectiousness [20]. In an epidemiology study, the basic reproduction number is a nondimensional quantity that sets the threshold during the study, both for predicting the outbreak and for evaluating the control strategies. Additionally, R0 analyzes the equilibrium stability, R0<1, which means that infectious individuals will cause less than one secondary infection and die out. Every infectious individual infects more than one secondary infection when R0>1, and the disease spreads to the population. In the SEQIHR model, the basic reproduction number is computed by using the next-generation matrix approach [21] and then obtained by taking the dominant eigenvalues (Spectral radius). Let Fix be the rate of new infection in compartment i and Vi be the rate of transfer of individuals into compartment i by all means other than the epidemic. The important thing is to obtain the disease-free equilibrium point E0. Thus, the computed matrices F and V which are n×n matrices, where n represents the infected classes, defined by: F=∂Fi/∂xjE0 and V=∂Vi/∂xjE0, where 1≤i, j≤n, F are nonnegative, and V is a nonsingular n-matrix (the matrix with inverse belongs to the class of positive matrices). Since F is nonnegative and V is a nonsingular matrix, then V−1 and FV−1 are nonnegative. Therefore, the next-generation matrix FV−1 is computed as defined by [22].Note that the basic reproduction number is defined as the spectral radius (dominant eigenvalue) of the matrixFV−1 [23], that is,
(35)R0=ρFV−1,FV−1=∂Fi∂xjE0∂Vi∂xjE0−1.where F is the rate of new infection in compartment I. The new forces of infection are
(36)dEdt=βIS−θ+q+bE,dQdt=ϕ+θE−α+η1+bQ,dIdt=1−ωqE+η1Q−η2+ν+δ+bI,dHdt=ωqE+η2I−λ+δ+bH.From Equation (36), when I and S meet, we obtain the following:
(37)Fi=f1f2f3f4=βIS000,From Equation (37), the Jacobian matrix of disease-free equilibrium (DFE) is given by
(38)FJ=00βS0000000000000.The partial derivative of Equation (38) with respect to E,Q,I, and H is given as
(39)V=b+θ+q000−θα+b+η100ω−1q−η1b+δ+η2+v0−ω0−η2b+δ+λ,Given thatA1=b+α+η1 and A2=b+δ+η2+ν. The inverse of V is
(40)V−1=1b+θ+q000θA1b+θ+q1A100−αω+α+b−ω+b+η1θ−η1ω+η1A1A2b+θ+qη1A1A21A20A1A2ω+η2−αω+α+b−ω+b+η1θ−η1ω+η1A1A2b+δ+λb+θ+qη1η2A1A2b+δ+λη2A2b+δ+λ1b+δ+λ.The product matrixFV−1 is given by:
(41)FV−1=βμη1θ−qω−1α+b+η1bα+b+η1b+θ+qb+δ+η2+νβη1μbα+b+η1b+δ+η2+νβμbb+δ+η2+ν0000000000000,Eigenvalues=0,0,0,βμ−bqω+bq+η1θ−αqω+αq−η1qω+η1qbα+b+η1b+θ+qb+δ+η2+ν,λ1=0,λ2=0,λ3=0,λ4=βμη1θ+q−ω+q−qω−1α+bbα+b+η1b+θ+qb+δ+η2+ν.Thus, the basic reproduction number becomes(42)R0=βμη1θ−ωq+q+q1−ωα+bbα+b+η1b+θ+qb+δ+η2+ν.


### 3.6. Existence of Endemic Equilibrium Point

Endemic equilibrium points are the steady-state solutions whereby the disease persists in the population [24]. The stability analysis of the endemic equilibrium point describes the long-term dynamics of COVID-19 in the population [25]. By solving all systems of differential equations from the model Equation (2), all derivatives are equal to zero (solve for all variables simultaneously).Theorem 5.
The endemic equilibrium point of model Equation (2) is locally asymptotically stable in the region Ω if R0<1 and unstable if R0>1.Proof.
At the endemic equilibrium pointS=S∗,E=E∗,Q=Q∗,I=I∗,H=H∗, and R=R∗. The variables are given as follows:
(43)S∗=μ+αQ+ρRb+βIE∗=βISb+θ+qQ∗=ϕ+θEb+α+η1H∗=Eωq+η2Ib+δ+λR∗=Hλ+νIb+ρ,I∗=μϕη1+A1A2R0−1+R0αQ∗+ρR∗A1A2βμA1A2−ϕη1.


### 3.7. Local Stability of the Disease-Free Equilibrium

The eigenvalues, which are determined by finding the partial derivatives of the vector-valued function, are used to study the local stability of the disease-free equilibrium. If the Jacobian matrix evaluated at that point has negative eigenvalues, the equilibrium point is asymptotically stable. The Routh-Hurwitz criterion in [26] will be utilized to demonstrate the local stability in this work.Theorem 6.
The disease-free equilibrium pointE0 is locally asymptotically stable if R0<1, and it is unstable when R0>1.Proof.
The linearization of the system of model 2 is done by computing its Jacobian matrix to prove this theorem. At the disease-free equilibrium point, the partial derivatives of each equation in the system for state variablesS,E,Q,I,H,R, which are used to generate the Jacobian matrix JE0 as in 21.
(44)JE0=−b0α−βS0ρ0−b−θ−q0βS000θ−α−b−η10000q1−ωη1−b−δ−η2−ν000ωq0η2−b−δ−λ0000νλ−b−ρ.
At a disease-free equilibrium,(45)S=μb,I=0.
The disease-free equilibrium will be asymptotically stable if the eigenvalues ofJE0<0.
(46)−b0α−βμb0ρ0−b−θ−q0βμb000θ−α−b−η10000q1−ωη1−b−δ−η2−ν000ωq0η2−b−δ−λ0000νλ−b−ρ=0.
From matrix (46), it is clear that the first, second, and third eigenvalues are
(47)λ1=−b,λ2=−b−ρ,andλ3=−b−δ−λ.
Then matrix (46) reduces to a 3×3 matrix after the cancellation of the respective rows and columns used to obtain the first, second, and third eigenvalues as shown in
(48)JE0=−b−θ−q0βμbθ−α−b−η10q1−ωη1−b−δ−η2−ν.
The characteristic polynomial of matrix (48) is given in the form
(49)Zλ=λ3+a1λ2+a2λ+a3,where
(50)a1=α+3b+δ+η1+η2+θ+ν+q,a2=αδ+αη2+αθ+αν+3b2+2αb+2bδ+2bη1+2bη2+2bθ+2bν−βμqb+2bq+δη1+δθ+η1θ+η2θ+η1ν+η1η2+θν+αq+δq+η1q+η2q+νq,a3=αδθ+αη2θ+αθν+b3+αb2+b2δ+b2η1+b2η2+b2θ+b2ν+b2q+αbδ+αbη2+αbθ+αbν−βη1θμb+bδη1+bδθ+bη1θ+bη2θ+bη1ν+bη1η2+bθν+αβμqωb−αβμqb+αbq+βη1μqωb−βη1μqb+bδq+bη1q+bη2q+bνq+δη1θ+η1θν+η1η2θ+αδq+αη2q+ανq+βμqω−βμq+δη1q+η1νq+η1η2q.
However,a1>0, a2>0, and a3>0, condition a1a2−a3>0.
(51)a1a2−a3=−βμqωα+b+M1M2−βμq+M3+η2M4+M5−βμq+η1M6+M7+βμθ−qωb,M1M2+M3+η2M4+M5+η1βθμ+M6+M7b>βμqωα+b+βμη2q+βμM1q+μqωβb,where
(52)M1=2b+δ+θ+ν+q,M2=4b3+2b22α+δ+θ+ν+q+bα+δ+να+θ+q,M3=bη122b+δ+η2+θ+ν+q,M4=8b3+bα+θ+qα+2δ+θ+2ν+q,M5=b26α+4δ+6θ+4ν+6q+bη2α+2b+θ+q,M6=8b3+b24α+6δ+θ+ν+6q+bδ+θ+ν+q2α+δ+θ+ν+q,M7=bη22α+3b+δ+θ+ν+q+η2.Hence, the conditiona1a2−a3>0 is satisfied. The Routh-Hurwitz criterion states that all elements of a system’s characteristic polynomial must be negative in order for the system to be stable [27]. The disease is asymptotically stable because the eigenvalues are negative, and the Routh-Hurwitz requirements are satisfied.


### 3.8. Global Stability of Disease-Free Equilibrium Point (DFE)

The global stability of theSEQIHR model around the DFE will be proved. The stability result of DFE in epidemiological implication is that minimizing the COVID-19 infection cases will not generate an infection if R0<1. Theorem 7 is considered.Theorem 7.
The DFE is globally asymptotically stable ifR0<1, and unstable if R0>1.Proof.
Using the technique described in [23, 28], the studies examine the global stability of the model 1 disease-free equilibrium point. The format can be used to write the SEQIHR model as in
(53)dXndt=AXn−Xdfe+A1Xi,dXidt=A2Xi.
By considering Equation (54), Xn is the vector of the nontransmitting compartment, Xi is the vector of transmitting compartment, and Xdfe is the vector of disease-free equilibrium point.
(54)Xn=S,RT,Xi=E,Q,I,HT,Xdfe=1b,0,(55)Xn−Xdfe=S−1bR.
We must show that the matrixA has real negative eigenvalues and A2 is a Metzler matrix in order for DFE to be globally stable (i.e., the off-diagonal elements of A2 are nonnegative, symbolically denoted by A2Xij≥∀≠j). We can derive equations with and without transmission from model 1, as stated herewith.
(56)μ+αQ+ρR−βIS−bSλH+νI−ρR−bR=AS−1bR+A1EQIH,βIS−θE−qE−bEϕ+θE−αQ−η1Q−bQ1−ωqE+η1Q−η2I−νI−δI−bIωqE+η2I−λH−δH−bH=A2EQIH,(57)then,A=−b00−b.
The eigenvalues of the matrixA are located at the diagonal (-b and -b), and these eigenvalues are real, distinct, and negative. Moreover, matrices A1 and A2 are given by
(58)A1=0α−βS000νλ,A2=−θ+q+b0βS0θ−α+η1+b001−ωqη1−η2+ν+δ+b0ωq0η2−λ+δ+b.A2 represents a Metzler matrix where its diagonal elements are negatives while the off-diagonal elements are nonnegative.


### 3.9. Global Stability of Endemic Equilibrium Point

The stability analysis explains the behavior of epidemic near the equilibrium points. The logarithmic Lyapunov function was proposed by Korobeinikov and Wake [29] to prove the global stability of endemic equilibrium for SIS, SIR, and SIRS models.Theorem 8.
The endemic equilibrium pointW∗ is asymptotically stable when R0>1 and unstable when R0<1.Proof.
The logarithmic Lyapunov function is used to analyze the stability of the endemic equilibrium and is given in the form(59)W=∑i=16aiXi−xi∗lnXi,
whereai represents a positive constant, Xi represents some free virus in compartment i, and Xi∗ denotes the number of free viruses in compartment i at the equilibrium point. Then, model system (2) is now written as follows:
(60)WS,E,Q,I,H,R=A1S−S∗lnS+A2E−E∗lnE+A3Q−Q∗lnQ+A4I−I∗lnI+A5H−H∗lnH+A6R−R∗lnR
The constantsA1,A2,A3,A4,A5, and A6 are nonnegative constants and the function W which is continuous and differentiable. Consider the derivative with respect to each compartment
(61)dWdt=A11−S∗SdSdt+A21−E∗EdEdt+A31−Q∗QdQdt+A41−I∗IdIdt+A51−H∗HdHdt+A61−R∗RdRdt.dWdt=A11−S∗Sμ+αQ+ρR−βIS−bS+A21−E∗EβIS−θ+q+bE+A31−Q∗Qϕ+θE−α+η1+bQ+A41−I∗I1−ωqE+η2I−η2+ν+δ+bI+A51−H∗HωqE+η2I−λ+δ+bH+A61−R∗RλH+νI−ρ+bR.
At the endemic equilibrium point,(62)μ=βIS+bS−αQ−ρR,βIS=θ+q+bE∗,ϕ+θE=α+η1+bQ∗,1−ωqE+η1Q=η1+ν+δ+bI∗,ωqE+η2I=λ+δ+bH∗,λH+νI=ρ+bR∗.(63)dWdt=A1S−S∗SbS∗−bS+A2E−E∗Eθ+q+bE∗−E+A3Q−Q∗Qα+η1+bQ∗−Q+A4I−I∗Iη1+ν+δ+bI∗−I+A5H−H∗Hλ+δ+bH∗−H+A6R−R∗Rρ+bR∗−R,where Equation (63) gives
(64)dWdt=−bA1S−S∗2S−θ+q+bA2E−E∗2E−α+η1+bA3Q−Q∗2Q−η1+ν+δ+bA4I−I∗2I−λ+δ+bA5H−H∗2H−ρ+bA6R−R∗2R.WhenS⟶S∗, E⟶E∗, Q⟶Q∗, I⟶I∗, H⟶H∗, and R⟶R∗. Therefore, dW/dt≤0 or zero and the function W is negative when WS,E,Q,I,H,R≥0.By following the approach of [29], the largest invariant set in Δ is a singleton set W which is the endemic equilibrium point, and using LaSalle [30] invariant principle, W∗ is globally asymptotically stable when R0>1 in Δ.


## 3.1. Positivity of the Model

For the model equations to be epidemiologically, we need to prove that the state variables are nonnegative∀t≥0.Theorem 1.
Let the initial data set beS,E,Q,I,H,R0>0. The solution set of the model system 2 is positive ∀t>0.Proof.
From the system of model Equation (2), consider the first equation:
(8)dSdt=μ+αQ+ρR−βI+bS.
By considering the negative term, by ignoring the rest, Equation (8) is reduced to
(9)dSdt≥−βI+bS.
This is the first-order linear differential inequality which can be solved by a separable methody′=fxgy (where S≥0) resulting in
(10)∫S0StdSS≥−∫0tβI+bdt,St≥S0e−βI+bt,in the absence of COVID-19 disease,
(11)St≥S0e−bt,S0≥0,∀t>0.By applying the same procedure to the remaining equations, the results areE0≥0, Q0≥0, I0≥0, H0≥0, and R0≥0. Therefore, the set of solutions St,Et,Qt,It,Ht, and Rt of the model is positive ∀t>0.


## 3.2. Invariant (Boundedness) Region

TheSEQIHR model is represented by differential equations in system 2, which is to be analyzed in a feasible region Ω, and all state variables and parameters of the model are assumed to be positive ∀t≥0. The bounded region is obtained through the following theorem.Theorem 2.
The setΩ is positively invariant and attracts all solutions in ℝ+6.Proof.
SinceNt=St+Qt+Et+It+Ht+Rt, then the derivative of Nt is given as
(12)dNdt=dSdt+dEdt+dQdt+dIdt+dHdt+dRdt.
Substituting model Equation (2) to Equation (12) gives
(13)dNdt=μ+αQ+ρR−βIS−bS+βIS−θE−qE−bE+ϕ+θE−αQ−η1Q−bQ+1−ωqE+η1Q−η2I−νI−δI−bI+ωqE+η2I−λH−δH−bH+λH+νI−ρR−bR.
Further simplification leads to(14)dNdt=μ+ϕ−S+E+Q+I+H+Rb−δI+δR,implying that
(15)dNdt=μ+ϕ−Nb−δI+δR.
For disease-free,δI⟹0 and δR⟹0, then dN/dt≤μ+ϕ−Nb.
This is a first-order linear differential equation, by a separable method.(16)∫N0NtdNμ+ϕ−bN≤∫0tdt,⟹Nt≥μ+ϕb1−e−bt+N0e−bt.Whent=0, N0≥0. When t⟹∞ then, N∞≤μ+ϕ/b.The Invariant region is given byΩ=S,E,Q,I,H,R∈ℝ+6:0≤Nt≤μ+ϕ/b. The SEQIHR model is biologically and epidemiologically meaningful; thus, we can consider the flow generated by the model for analysis.


## 3.3. Existence and Uniqueness of the Solution

From the first-order differential equation given in the form:y′=ft,y,yt0=y0. The following questions will be of interest:
(i)
Under what conditions can we say that a solution to the equationy′ exists?(ii)
Under what conditions does a unique solution exist to the equationy′?Consider the following equations to answer the question.(17)f1=μ+αQ+δR−βIS−bS,(18)f2=βIS−θE−qE−bE,(19)f3=ϕ+θE−α+η1+bQ,(20)f4=1−ωqE+η1Q−η2+ν+δ+bI,(21)f5=ωqE+η2I−λ+δ+bH,(22)f6=λH+νI−ρ+bR.


### 3.3.1. Uniqueness of Solution

Theorem 3.
Let us useD to denote the domain:
(23)t−t0≤a,y−y0≤b,y=y1,y2,⋯,yn,y0=y10,y20,⋯,yn0.Proof.
Supposeft,y satisfies the Lipschitz condition; therefore,
(24)ft,y2−ft,y1≤kx2−x1.
Whenever the pointst,x1 and t,x2 belong to the domain D and k is used to represent the positive constant, then, there exist a constant α>0 and a unique solution yt of system 11 in the interval t−t0≤α. It is essential to note that condition (24) is satisfied by ∂fi/∂yj,i,j=1,2,3,⋯,n to be continuous and bounded in domain D.Ifft,y has a continuous partial derivative ∂fi/∂yj on a bounded closed convex domain ℝ (i.e., the convex set of real numbers), where ℝ is used to denote real numbers; then it satisfies a Lipschitz condition in ℝ. Our interest is in the domain:
(25)1≤ε≤ℝ.Therefore, we look for the bounded solution of the form0<ℝ<∞.


### 3.3.2. Existence of a Solution

Theorem 4.
Let the domain be denoted byD, also defined in Equation (23), such that Equations (24) and (25) hold. Then, the existing solution of Equation (2) is bounded in the domain D.Proof.
From Equations (10)–(15), by showing that ∂fi/∂yj,i,j=1,2,3,4,5,6, then, Equations (1) and (2) are continuous and bounded. That means the partial derivatives are continuous and bound. Consider the exploration of the partial derivatives for all model equations.
From Equation (17), we obtain the following system of equations:
(26)∂f1∂S=−βI−b,∂f1∂S=−βI−b<∞,∂f1∂E=0,∂f1∂E=0<∞,∂f1∂Q=α,∂f1∂Q=α<∞,∂f1∂I=−βS,∂f1∂I=−βS<∞,∂f1∂H=0,∂f1∂H=0<∞,∂f1∂R=ρ,∂f1∂R=ρ<∞.
Similarly, from Equation (18), we obtain the following system of equations:
(27)∂f2∂S=βI,∂f2∂S=βI<∞,∂f2∂E=−θ+q+b,∂f2∂E=−θ+q+b<∞,∂f2∂Q=0,δf2∂Q=0<∞,∂f2∂I=βS,∂f2∂I=βS<∞,∂f2∂H=0,∂f2∂H=0<∞,∂f2∂R=0,∂f2∂R=0<∞.The same procedures are taken for Equations (19), (20), (21), and (22). Therefore, all partial derivatives are continuous and bounded; hence, from Theorem 4, it is concluded that there exists a unique solution of the model in Equation (2) in the domain region D.


## 3.3.1. Uniqueness of Solution

Theorem 3.
Let us useD to denote the domain:
(23)t−t0≤a,y−y0≤b,y=y1,y2,⋯,yn,y0=y10,y20,⋯,yn0.Proof.
Supposeft,y satisfies the Lipschitz condition; therefore,
(24)ft,y2−ft,y1≤kx2−x1.
Whenever the pointst,x1 and t,x2 belong to the domain D and k is used to represent the positive constant, then, there exist a constant α>0 and a unique solution yt of system 11 in the interval t−t0≤α. It is essential to note that condition (24) is satisfied by ∂fi/∂yj,i,j=1,2,3,⋯,n to be continuous and bounded in domain D.Ifft,y has a continuous partial derivative ∂fi/∂yj on a bounded closed convex domain ℝ (i.e., the convex set of real numbers), where ℝ is used to denote real numbers; then it satisfies a Lipschitz condition in ℝ. Our interest is in the domain:
(25)1≤ε≤ℝ.Therefore, we look for the bounded solution of the form0<ℝ<∞.


## 3.3.2. Existence of a Solution

Theorem 4.
Let the domain be denoted byD, also defined in Equation (23), such that Equations (24) and (25) hold. Then, the existing solution of Equation (2) is bounded in the domain D.Proof.
From Equations (10)–(15), by showing that ∂fi/∂yj,i,j=1,2,3,4,5,6, then, Equations (1) and (2) are continuous and bounded. That means the partial derivatives are continuous and bound. Consider the exploration of the partial derivatives for all model equations.
From Equation (17), we obtain the following system of equations:
(26)∂f1∂S=−βI−b,∂f1∂S=−βI−b<∞,∂f1∂E=0,∂f1∂E=0<∞,∂f1∂Q=α,∂f1∂Q=α<∞,∂f1∂I=−βS,∂f1∂I=−βS<∞,∂f1∂H=0,∂f1∂H=0<∞,∂f1∂R=ρ,∂f1∂R=ρ<∞.
Similarly, from Equation (18), we obtain the following system of equations:
(27)∂f2∂S=βI,∂f2∂S=βI<∞,∂f2∂E=−θ+q+b,∂f2∂E=−θ+q+b<∞,∂f2∂Q=0,δf2∂Q=0<∞,∂f2∂I=βS,∂f2∂I=βS<∞,∂f2∂H=0,∂f2∂H=0<∞,∂f2∂R=0,∂f2∂R=0<∞.The same procedures are taken for Equations (19), (20), (21), and (22). Therefore, all partial derivatives are continuous and bounded; hence, from Theorem 4, it is concluded that there exists a unique solution of the model in Equation (2) in the domain region D.


## 3.4. Existence of Disease-Free Equilibrium Point (DFE)

The disease-free equilibrium point, obtained when the infected components are zero, can be done by setting the right-hand side of the equation equal to zero, as in(28)dSdt=dEdt=dQdt=dIdt=dHdt=dRdt=0.When there is no disease, then,E=0, Q=0, I=0, H=0, and R=0.By considering each model equation from Equation (2),
(29)μ+αQ+ρR−βIS−bS=0givesμ−bS=0,then,S=μb.In addition to the second model equation,(30)βIS=θ+q+bE=0givesE=βISθ+q+b,butβ=0,then,E=0.Similarly, for the third model equation,(31)ϕ+θE−αQ−η1Q−bQ=0givesQ=ϕ+θEα+η1+b,but,E=0,then,
(32)Q=ϕα+η1+b.By considering the same procedures for the fourth, fifth, and sixth model equations, the following is obtained:(33)I=0,H=0,R=0.Therefore,(34)E0=S0,E0,Q0,I0,H0,R0=μb,0,ϕα+η1+b,0,0,0T.Equation (34) represents the state in which there is no infection and is known as the disease-free equilibrium point.


## 3.5. Basic Reproduction Number (R0)

The basic reproduction numberR0 is the midpoint number of infections caused by an infectious individual during the entire period of infectiousness [20]. In an epidemiology study, the basic reproduction number is a nondimensional quantity that sets the threshold during the study, both for predicting the outbreak and for evaluating the control strategies. Additionally, R0 analyzes the equilibrium stability, R0<1, which means that infectious individuals will cause less than one secondary infection and die out. Every infectious individual infects more than one secondary infection when R0>1, and the disease spreads to the population. In the SEQIHR model, the basic reproduction number is computed by using the next-generation matrix approach [21] and then obtained by taking the dominant eigenvalues (Spectral radius). Let Fix be the rate of new infection in compartment i and Vi be the rate of transfer of individuals into compartment i by all means other than the epidemic. The important thing is to obtain the disease-free equilibrium point E0. Thus, the computed matrices F and V which are n×n matrices, where n represents the infected classes, defined by: F=∂Fi/∂xjE0 and V=∂Vi/∂xjE0, where 1≤i, j≤n, F are nonnegative, and V is a nonsingular n-matrix (the matrix with inverse belongs to the class of positive matrices). Since F is nonnegative and V is a nonsingular matrix, then V−1 and FV−1 are nonnegative. Therefore, the next-generation matrix FV−1 is computed as defined by [22].Note that the basic reproduction number is defined as the spectral radius (dominant eigenvalue) of the matrixFV−1 [23], that is,
(35)R0=ρFV−1,FV−1=∂Fi∂xjE0∂Vi∂xjE0−1.where F is the rate of new infection in compartment I. The new forces of infection are
(36)dEdt=βIS−θ+q+bE,dQdt=ϕ+θE−α+η1+bQ,dIdt=1−ωqE+η1Q−η2+ν+δ+bI,dHdt=ωqE+η2I−λ+δ+bH.From Equation (36), when I and S meet, we obtain the following:
(37)Fi=f1f2f3f4=βIS000,From Equation (37), the Jacobian matrix of disease-free equilibrium (DFE) is given by
(38)FJ=00βS0000000000000.The partial derivative of Equation (38) with respect to E,Q,I, and H is given as
(39)V=b+θ+q000−θα+b+η100ω−1q−η1b+δ+η2+v0−ω0−η2b+δ+λ,Given thatA1=b+α+η1 and A2=b+δ+η2+ν. The inverse of V is
(40)V−1=1b+θ+q000θA1b+θ+q1A100−αω+α+b−ω+b+η1θ−η1ω+η1A1A2b+θ+qη1A1A21A20A1A2ω+η2−αω+α+b−ω+b+η1θ−η1ω+η1A1A2b+δ+λb+θ+qη1η2A1A2b+δ+λη2A2b+δ+λ1b+δ+λ.The product matrixFV−1 is given by:
(41)FV−1=βμη1θ−qω−1α+b+η1bα+b+η1b+θ+qb+δ+η2+νβη1μbα+b+η1b+δ+η2+νβμbb+δ+η2+ν0000000000000,Eigenvalues=0,0,0,βμ−bqω+bq+η1θ−αqω+αq−η1qω+η1qbα+b+η1b+θ+qb+δ+η2+ν,λ1=0,λ2=0,λ3=0,λ4=βμη1θ+q−ω+q−qω−1α+bbα+b+η1b+θ+qb+δ+η2+ν.Thus, the basic reproduction number becomes(42)R0=βμη1θ−ωq+q+q1−ωα+bbα+b+η1b+θ+qb+δ+η2+ν.


## 3.6. Existence of Endemic Equilibrium Point

Endemic equilibrium points are the steady-state solutions whereby the disease persists in the population [24]. The stability analysis of the endemic equilibrium point describes the long-term dynamics of COVID-19 in the population [25]. By solving all systems of differential equations from the model Equation (2), all derivatives are equal to zero (solve for all variables simultaneously).Theorem 5.
The endemic equilibrium point of model Equation (2) is locally asymptotically stable in the region Ω if R0<1 and unstable if R0>1.Proof.
At the endemic equilibrium pointS=S∗,E=E∗,Q=Q∗,I=I∗,H=H∗, and R=R∗. The variables are given as follows:
(43)S∗=μ+αQ+ρRb+βIE∗=βISb+θ+qQ∗=ϕ+θEb+α+η1H∗=Eωq+η2Ib+δ+λR∗=Hλ+νIb+ρ,I∗=μϕη1+A1A2R0−1+R0αQ∗+ρR∗A1A2βμA1A2−ϕη1.


## 3.7. Local Stability of the Disease-Free Equilibrium

The eigenvalues, which are determined by finding the partial derivatives of the vector-valued function, are used to study the local stability of the disease-free equilibrium. If the Jacobian matrix evaluated at that point has negative eigenvalues, the equilibrium point is asymptotically stable. The Routh-Hurwitz criterion in [26] will be utilized to demonstrate the local stability in this work.Theorem 6.
The disease-free equilibrium pointE0 is locally asymptotically stable if R0<1, and it is unstable when R0>1.Proof.
The linearization of the system of model 2 is done by computing its Jacobian matrix to prove this theorem. At the disease-free equilibrium point, the partial derivatives of each equation in the system for state variablesS,E,Q,I,H,R, which are used to generate the Jacobian matrix JE0 as in 21.
(44)JE0=−b0α−βS0ρ0−b−θ−q0βS000θ−α−b−η10000q1−ωη1−b−δ−η2−ν000ωq0η2−b−δ−λ0000νλ−b−ρ.
At a disease-free equilibrium,(45)S=μb,I=0.
The disease-free equilibrium will be asymptotically stable if the eigenvalues ofJE0<0.
(46)−b0α−βμb0ρ0−b−θ−q0βμb000θ−α−b−η10000q1−ωη1−b−δ−η2−ν000ωq0η2−b−δ−λ0000νλ−b−ρ=0.
From matrix (46), it is clear that the first, second, and third eigenvalues are
(47)λ1=−b,λ2=−b−ρ,andλ3=−b−δ−λ.
Then matrix (46) reduces to a 3×3 matrix after the cancellation of the respective rows and columns used to obtain the first, second, and third eigenvalues as shown in
(48)JE0=−b−θ−q0βμbθ−α−b−η10q1−ωη1−b−δ−η2−ν.
The characteristic polynomial of matrix (48) is given in the form
(49)Zλ=λ3+a1λ2+a2λ+a3,where
(50)a1=α+3b+δ+η1+η2+θ+ν+q,a2=αδ+αη2+αθ+αν+3b2+2αb+2bδ+2bη1+2bη2+2bθ+2bν−βμqb+2bq+δη1+δθ+η1θ+η2θ+η1ν+η1η2+θν+αq+δq+η1q+η2q+νq,a3=αδθ+αη2θ+αθν+b3+αb2+b2δ+b2η1+b2η2+b2θ+b2ν+b2q+αbδ+αbη2+αbθ+αbν−βη1θμb+bδη1+bδθ+bη1θ+bη2θ+bη1ν+bη1η2+bθν+αβμqωb−αβμqb+αbq+βη1μqωb−βη1μqb+bδq+bη1q+bη2q+bνq+δη1θ+η1θν+η1η2θ+αδq+αη2q+ανq+βμqω−βμq+δη1q+η1νq+η1η2q.
However,a1>0, a2>0, and a3>0, condition a1a2−a3>0.
(51)a1a2−a3=−βμqωα+b+M1M2−βμq+M3+η2M4+M5−βμq+η1M6+M7+βμθ−qωb,M1M2+M3+η2M4+M5+η1βθμ+M6+M7b>βμqωα+b+βμη2q+βμM1q+μqωβb,where
(52)M1=2b+δ+θ+ν+q,M2=4b3+2b22α+δ+θ+ν+q+bα+δ+να+θ+q,M3=bη122b+δ+η2+θ+ν+q,M4=8b3+bα+θ+qα+2δ+θ+2ν+q,M5=b26α+4δ+6θ+4ν+6q+bη2α+2b+θ+q,M6=8b3+b24α+6δ+θ+ν+6q+bδ+θ+ν+q2α+δ+θ+ν+q,M7=bη22α+3b+δ+θ+ν+q+η2.Hence, the conditiona1a2−a3>0 is satisfied. The Routh-Hurwitz criterion states that all elements of a system’s characteristic polynomial must be negative in order for the system to be stable [27]. The disease is asymptotically stable because the eigenvalues are negative, and the Routh-Hurwitz requirements are satisfied.


## 3.8. Global Stability of Disease-Free Equilibrium Point (DFE)

The global stability of theSEQIHR model around the DFE will be proved. The stability result of DFE in epidemiological implication is that minimizing the COVID-19 infection cases will not generate an infection if R0<1. Theorem 7 is considered.Theorem 7.
The DFE is globally asymptotically stable ifR0<1, and unstable if R0>1.Proof.
Using the technique described in [23, 28], the studies examine the global stability of the model 1 disease-free equilibrium point. The format can be used to write the SEQIHR model as in
(53)dXndt=AXn−Xdfe+A1Xi,dXidt=A2Xi.
By considering Equation (54), Xn is the vector of the nontransmitting compartment, Xi is the vector of transmitting compartment, and Xdfe is the vector of disease-free equilibrium point.
(54)Xn=S,RT,Xi=E,Q,I,HT,Xdfe=1b,0,(55)Xn−Xdfe=S−1bR.
We must show that the matrixA has real negative eigenvalues and A2 is a Metzler matrix in order for DFE to be globally stable (i.e., the off-diagonal elements of A2 are nonnegative, symbolically denoted by A2Xij≥∀≠j). We can derive equations with and without transmission from model 1, as stated herewith.
(56)μ+αQ+ρR−βIS−bSλH+νI−ρR−bR=AS−1bR+A1EQIH,βIS−θE−qE−bEϕ+θE−αQ−η1Q−bQ1−ωqE+η1Q−η2I−νI−δI−bIωqE+η2I−λH−δH−bH=A2EQIH,(57)then,A=−b00−b.
The eigenvalues of the matrixA are located at the diagonal (-b and -b), and these eigenvalues are real, distinct, and negative. Moreover, matrices A1 and A2 are given by
(58)A1=0α−βS000νλ,A2=−θ+q+b0βS0θ−α+η1+b001−ωqη1−η2+ν+δ+b0ωq0η2−λ+δ+b.A2 represents a Metzler matrix where its diagonal elements are negatives while the off-diagonal elements are nonnegative.


## 3.9. Global Stability of Endemic Equilibrium Point

The stability analysis explains the behavior of epidemic near the equilibrium points. The logarithmic Lyapunov function was proposed by Korobeinikov and Wake [29] to prove the global stability of endemic equilibrium for SIS, SIR, and SIRS models.Theorem 8.
The endemic equilibrium pointW∗ is asymptotically stable when R0>1 and unstable when R0<1.Proof.
The logarithmic Lyapunov function is used to analyze the stability of the endemic equilibrium and is given in the form(59)W=∑i=16aiXi−xi∗lnXi,
whereai represents a positive constant, Xi represents some free virus in compartment i, and Xi∗ denotes the number of free viruses in compartment i at the equilibrium point. Then, model system (2) is now written as follows:
(60)WS,E,Q,I,H,R=A1S−S∗lnS+A2E−E∗lnE+A3Q−Q∗lnQ+A4I−I∗lnI+A5H−H∗lnH+A6R−R∗lnR
The constantsA1,A2,A3,A4,A5, and A6 are nonnegative constants and the function W which is continuous and differentiable. Consider the derivative with respect to each compartment
(61)dWdt=A11−S∗SdSdt+A21−E∗EdEdt+A31−Q∗QdQdt+A41−I∗IdIdt+A51−H∗HdHdt+A61−R∗RdRdt.dWdt=A11−S∗Sμ+αQ+ρR−βIS−bS+A21−E∗EβIS−θ+q+bE+A31−Q∗Qϕ+θE−α+η1+bQ+A41−I∗I1−ωqE+η2I−η2+ν+δ+bI+A51−H∗HωqE+η2I−λ+δ+bH+A61−R∗RλH+νI−ρ+bR.
At the endemic equilibrium point,(62)μ=βIS+bS−αQ−ρR,βIS=θ+q+bE∗,ϕ+θE=α+η1+bQ∗,1−ωqE+η1Q=η1+ν+δ+bI∗,ωqE+η2I=λ+δ+bH∗,λH+νI=ρ+bR∗.(63)dWdt=A1S−S∗SbS∗−bS+A2E−E∗Eθ+q+bE∗−E+A3Q−Q∗Qα+η1+bQ∗−Q+A4I−I∗Iη1+ν+δ+bI∗−I+A5H−H∗Hλ+δ+bH∗−H+A6R−R∗Rρ+bR∗−R,where Equation (63) gives
(64)dWdt=−bA1S−S∗2S−θ+q+bA2E−E∗2E−α+η1+bA3Q−Q∗2Q−η1+ν+δ+bA4I−I∗2I−λ+δ+bA5H−H∗2H−ρ+bA6R−R∗2R.WhenS⟶S∗, E⟶E∗, Q⟶Q∗, I⟶I∗, H⟶H∗, and R⟶R∗. Therefore, dW/dt≤0 or zero and the function W is negative when WS,E,Q,I,H,R≥0.By following the approach of [29], the largest invariant set in Δ is a singleton set W which is the endemic equilibrium point, and using LaSalle [30] invariant principle, W∗ is globally asymptotically stable when R0>1 in Δ.


## 4. Numerical Simulation

A series of numerical results for system (2) of the model equations are presented. The explicit Runge-Kutta fourth-order method is considered for solving the first-order ordinary differential equations of SEQIHR model numerically with a given initial condition. Partial Rank Correlation Coefficient (PRCC) was used to show the sensitivity analysis of the parameters and basic reproduction number. Parameter values from the literature reviews were used, and some were assumed as shown in Table 3. The data is simulated by substituting them in R0. The sensitivity index of each partial basic reproduction number R0 for its parameters.Table 3
Parameter values and sensitivity indices.
ParameterValuesSourcesSensitivity index valueβ0.015Assumed1.0000μ8.94Assumed1.0000ω0.083Assumed-0.0388b0.0104[31]-0.9483α0.85Assumed0.7100ν0.07[32]0.090θ0.2435[33]1.2614η10.85Assumed1.2813η20.65Assumed0.8448δ0.039[34]0.0507q0.099Assumed0.7092


### 4.1. Sensitivity Analysis and Uncertainty

The sensitivity analysis for the endemic threshold tells us the importance of each parameter for the transmission of COVID-19 disease. The information is crucial for the analysis of complex systems. We used the sensitivity analysis of the parameters to determine the strongness of theSEQIHR model predictions for the parameter values. There are usually errors in the data collected and in the initial values assumed for the parameters [35]. The standard equation of a sensitivity index for R0 is given by
(65)ΓLR0=∂R0∂L×LR0,.From Table3, it is observed that the η1 parameter is more sensitive since it increases the basic reproductive number by ΓLR0=1.2813. The increase of this parameter means that in quarantine, people are more stressed such that the immune system decreases in its efficiency, so the virus spreads within the body. The virus causes blood clotting because the virus fights the respiratory system and enters the bloodstream through lung capillaries that are adjacent to the alveolus [36]. The PRCC supported graphically in Figure 2 shows that η1 has impact on the transmission of COVID-19 disease. SEQIHR model shows that hospitalized patients from the infected class are less stressed than those in quarantine, although it also increases the basic reproduction number by (ΓLR0=0.8448); then η1>η2. Other sensitive parameters are θ with ΓLR0=1.2614, β and μ which both have ΓLR0=1. Parameters, ω and b, have the least values which are -0.0717 and -0.9566, respectively, meaning that decreasing ω and b by a certain percentage always decreases R0 by the same percentage and the same thing happens if we increase R0. It shows that if most of the patients are hospitalized, the disease will decrease, and subsequently, there will be no more transmission within the community. Despite the fact that the other parameters have small values, they still increase R0 by their respective percentages.Figure 2
Global sensitivity analysis and PRCC results forR0.From Figure2, η1 is positively and highly correlated with R0 as the absolute value of its PRCC value is higher than the corresponding value of other parameters. Furthermore, natural death b is highly negatively correlated with R0.


#### 4.1.1. Dynamic Population Simulation with aSEQIHR Model

The numerical simulation of theSEQIHR model variables is shown in Figure 3. We observe that the susceptible class declines to acquire the endemic equilibrium level exponentially as people die naturally or due to the disease. The exposed, quarantined, and infected populations both assume a parabolic curve which increases exponentially to a certain maximum point before they decelerate to an endemic level. Hospitalized and recovered populations both assumed a parabolic shape as it increases exponentially to a certain maximum point before decelerating to an endemic point.Figure 3
Dynamic simulation of aSEQIHR model.


#### 4.1.2. Relationship between Most Positive and Negative Parameters with Basic Reproduction NumberR0

From Section4.1, we observed that the most positive parameter is η1. The relationship between η1 and R0 is shown in Figure 4(a), whereby we noticed that the increase of stress results in a quick rise of R0. This shows that individuals isolated from their families are more stressed, which increases the disease infection because stress lowers body immunity. Similarly, the natural death b is the most negative sensitive parameter shown in Figure 4(b), which means that any reduction in it will make the basic reproduction number experience significant exponential retardation.Figure 4
Effect of stressη1 and natural death b on R0.
(a)(b)


#### 4.1.3. Simulation of Stability Analysis ofSEQIHR Model

The numerical simulations for the stability analysis are performed to the analytical results of the model. For the equilibrium point to be globally asymptomatic stable, the model trajectories for the state variables should be originated from different initial values and sometimes converge to a common point and maintain an endemic equilibrium point. But if the model trajectories for state variables remain near the equilibrium point and move together in the long run, this implies that the equilibrium point is globally stable. The six trajectories in each class are represented by different colors as shown by a legend that converge towards equilibrium point as time approaches infinity. The model variablesS,E,Q,I,H, and R varied by considering Figures 5(a) and 7(b), respectively, are illustrated as follows:Figure 5
Global stability of endemic equilibrium for susceptible and exposed human population.
(a)(b)Figure 6
Global stability of endemic equilibrium for quarantined and infected human population.
(a)(b)Figure 7
Global stability of endemic equilibrium for hospitalized and recovered human population.
(a)(b)


#### 4.1.4. Control Parameters

By considering the variation of some control parameters; from Figure8(a), when the contact rate (β) increases, the decaying population rate also increases, but when the contact rate decreases, the decaying rate also decreases. From Figure 8(b), when the contact rate increases, the infections also increase and apply the same when the contact rate decreases and the infections decreases. When stress increases in the susceptible class, the rate of decaying increases, and if the stress decreases, the decaying rate decreases, as shown in Figure 9(a). From Figure 9(b), when stress (η1) increases, the rate of infection increases, and when η1 decreases, then the rate of infected population decreases. Moreover, when η2 increases, the rate of infected and hospitalized increases and then decreases when η2 decreases as in Figures 10(a) and 10(b).Figure 8
Variation of contact rateβ in susceptible and infected population.
(a)(b)Figure 9
Variation ofη1 in susceptible and infected human population.
(a)(b)Figure 10
Variation ofη2 in infected and hospitalized human population.
(a)(b)


#### 4.1.5. Parameter Identifiability and Model Fitting

The identifiability of parameters is essential to the proposedSEQIHR model. Such parameters are β, μ, ω, α, η1, and η2. Parameter identifiability is implemented by using the least square method to minimize the sum of squared differences between the observations and the SEQIHR model [37] and is defined as
(66)SSθ=∑i=1nyi−fxi,θ2,where yi are the observed data of all compartments, i is the number of compartments (i.e., i=1,2,⋯,n), and fxi,θ is the solution for all compartments of the SEQIHR model. With the initial values of the parameters given as β, μ, ω, α, η1, and η2, the least squares identifiabilities are obtained as shown in Table 4 as the initial value and estimated values and used to fit the simulated data as shown in Figure 11. The relationship of initial parameter values and the identifiable values from the least square method is very close.Table 4
Parameter identifiability.
ParameterInitial valuesEstimated valuesβ0.0150.0147μ8.949.8ω0.0830.096b0.01040.013α0.851.1ν0.070.063θ0.24350.224η10.850.72η20.650.51δ0.0390.038q0.0990.098ρ0.0030.0032λ0.0020.0018ϕ0.0010.0009Figure 11
Fitted graphs.


## 4.1. Sensitivity Analysis and Uncertainty

The sensitivity analysis for the endemic threshold tells us the importance of each parameter for the transmission of COVID-19 disease. The information is crucial for the analysis of complex systems. We used the sensitivity analysis of the parameters to determine the strongness of theSEQIHR model predictions for the parameter values. There are usually errors in the data collected and in the initial values assumed for the parameters [35]. The standard equation of a sensitivity index for R0 is given by
(65)ΓLR0=∂R0∂L×LR0,.From Table3, it is observed that the η1 parameter is more sensitive since it increases the basic reproductive number by ΓLR0=1.2813. The increase of this parameter means that in quarantine, people are more stressed such that the immune system decreases in its efficiency, so the virus spreads within the body. The virus causes blood clotting because the virus fights the respiratory system and enters the bloodstream through lung capillaries that are adjacent to the alveolus [36]. The PRCC supported graphically in Figure 2 shows that η1 has impact on the transmission of COVID-19 disease. SEQIHR model shows that hospitalized patients from the infected class are less stressed than those in quarantine, although it also increases the basic reproduction number by (ΓLR0=0.8448); then η1>η2. Other sensitive parameters are θ with ΓLR0=1.2614, β and μ which both have ΓLR0=1. Parameters, ω and b, have the least values which are -0.0717 and -0.9566, respectively, meaning that decreasing ω and b by a certain percentage always decreases R0 by the same percentage and the same thing happens if we increase R0. It shows that if most of the patients are hospitalized, the disease will decrease, and subsequently, there will be no more transmission within the community. Despite the fact that the other parameters have small values, they still increase R0 by their respective percentages.Figure 2
Global sensitivity analysis and PRCC results forR0.From Figure2, η1 is positively and highly correlated with R0 as the absolute value of its PRCC value is higher than the corresponding value of other parameters. Furthermore, natural death b is highly negatively correlated with R0.


### 4.1.1. Dynamic Population Simulation with aSEQIHR Model

The numerical simulation of theSEQIHR model variables is shown in Figure 3. We observe that the susceptible class declines to acquire the endemic equilibrium level exponentially as people die naturally or due to the disease. The exposed, quarantined, and infected populations both assume a parabolic curve which increases exponentially to a certain maximum point before they decelerate to an endemic level. Hospitalized and recovered populations both assumed a parabolic shape as it increases exponentially to a certain maximum point before decelerating to an endemic point.Figure 3
Dynamic simulation of aSEQIHR model.


### 4.1.2. Relationship between Most Positive and Negative Parameters with Basic Reproduction NumberR0

From Section4.1, we observed that the most positive parameter is η1. The relationship between η1 and R0 is shown in Figure 4(a), whereby we noticed that the increase of stress results in a quick rise of R0. This shows that individuals isolated from their families are more stressed, which increases the disease infection because stress lowers body immunity. Similarly, the natural death b is the most negative sensitive parameter shown in Figure 4(b), which means that any reduction in it will make the basic reproduction number experience significant exponential retardation.Figure 4
Effect of stressη1 and natural death b on R0.
(a)(b)


### 4.1.3. Simulation of Stability Analysis ofSEQIHR Model

The numerical simulations for the stability analysis are performed to the analytical results of the model. For the equilibrium point to be globally asymptomatic stable, the model trajectories for the state variables should be originated from different initial values and sometimes converge to a common point and maintain an endemic equilibrium point. But if the model trajectories for state variables remain near the equilibrium point and move together in the long run, this implies that the equilibrium point is globally stable. The six trajectories in each class are represented by different colors as shown by a legend that converge towards equilibrium point as time approaches infinity. The model variablesS,E,Q,I,H, and R varied by considering Figures 5(a) and 7(b), respectively, are illustrated as follows:Figure 5
Global stability of endemic equilibrium for susceptible and exposed human population.
(a)(b)Figure 6
Global stability of endemic equilibrium for quarantined and infected human population.
(a)(b)Figure 7
Global stability of endemic equilibrium for hospitalized and recovered human population.
(a)(b)


### 4.1.4. Control Parameters

By considering the variation of some control parameters; from Figure8(a), when the contact rate (β) increases, the decaying population rate also increases, but when the contact rate decreases, the decaying rate also decreases. From Figure 8(b), when the contact rate increases, the infections also increase and apply the same when the contact rate decreases and the infections decreases. When stress increases in the susceptible class, the rate of decaying increases, and if the stress decreases, the decaying rate decreases, as shown in Figure 9(a). From Figure 9(b), when stress (η1) increases, the rate of infection increases, and when η1 decreases, then the rate of infected population decreases. Moreover, when η2 increases, the rate of infected and hospitalized increases and then decreases when η2 decreases as in Figures 10(a) and 10(b).Figure 8
Variation of contact rateβ in susceptible and infected population.
(a)(b)Figure 9
Variation ofη1 in susceptible and infected human population.
(a)(b)Figure 10
Variation ofη2 in infected and hospitalized human population.
(a)(b)


### 4.1.5. Parameter Identifiability and Model Fitting

The identifiability of parameters is essential to the proposedSEQIHR model. Such parameters are β, μ, ω, α, η1, and η2. Parameter identifiability is implemented by using the least square method to minimize the sum of squared differences between the observations and the SEQIHR model [37] and is defined as
(66)SSθ=∑i=1nyi−fxi,θ2,where yi are the observed data of all compartments, i is the number of compartments (i.e., i=1,2,⋯,n), and fxi,θ is the solution for all compartments of the SEQIHR model. With the initial values of the parameters given as β, μ, ω, α, η1, and η2, the least squares identifiabilities are obtained as shown in Table 4 as the initial value and estimated values and used to fit the simulated data as shown in Figure 11. The relationship of initial parameter values and the identifiable values from the least square method is very close.Table 4
Parameter identifiability.
ParameterInitial valuesEstimated valuesβ0.0150.0147μ8.949.8ω0.0830.096b0.01040.013α0.851.1ν0.070.063θ0.24350.224η10.850.72η20.650.51δ0.0390.038q0.0990.098ρ0.0030.0032λ0.0020.0018ϕ0.0010.0009Figure 11
Fitted graphs.


## 4.1.1. Dynamic Population Simulation with aSEQIHR Model

The numerical simulation of theSEQIHR model variables is shown in Figure 3. We observe that the susceptible class declines to acquire the endemic equilibrium level exponentially as people die naturally or due to the disease. The exposed, quarantined, and infected populations both assume a parabolic curve which increases exponentially to a certain maximum point before they decelerate to an endemic level. Hospitalized and recovered populations both assumed a parabolic shape as it increases exponentially to a certain maximum point before decelerating to an endemic point.Figure 3
Dynamic simulation of aSEQIHR model.


## 4.1.2. Relationship between Most Positive and Negative Parameters with Basic Reproduction NumberR0

From Section4.1, we observed that the most positive parameter is η1. The relationship between η1 and R0 is shown in Figure 4(a), whereby we noticed that the increase of stress results in a quick rise of R0. This shows that individuals isolated from their families are more stressed, which increases the disease infection because stress lowers body immunity. Similarly, the natural death b is the most negative sensitive parameter shown in Figure 4(b), which means that any reduction in it will make the basic reproduction number experience significant exponential retardation.Figure 4
Effect of stressη1 and natural death b on R0.
(a)(b)


## 4.1.3. Simulation of Stability Analysis ofSEQIHR Model

The numerical simulations for the stability analysis are performed to the analytical results of the model. For the equilibrium point to be globally asymptomatic stable, the model trajectories for the state variables should be originated from different initial values and sometimes converge to a common point and maintain an endemic equilibrium point. But if the model trajectories for state variables remain near the equilibrium point and move together in the long run, this implies that the equilibrium point is globally stable. The six trajectories in each class are represented by different colors as shown by a legend that converge towards equilibrium point as time approaches infinity. The model variablesS,E,Q,I,H, and R varied by considering Figures 5(a) and 7(b), respectively, are illustrated as follows:Figure 5
Global stability of endemic equilibrium for susceptible and exposed human population.
(a)(b)Figure 6
Global stability of endemic equilibrium for quarantined and infected human population.
(a)(b)Figure 7
Global stability of endemic equilibrium for hospitalized and recovered human population.
(a)(b)


## 4.1.4. Control Parameters

By considering the variation of some control parameters; from Figure8(a), when the contact rate (β) increases, the decaying population rate also increases, but when the contact rate decreases, the decaying rate also decreases. From Figure 8(b), when the contact rate increases, the infections also increase and apply the same when the contact rate decreases and the infections decreases. When stress increases in the susceptible class, the rate of decaying increases, and if the stress decreases, the decaying rate decreases, as shown in Figure 9(a). From Figure 9(b), when stress (η1) increases, the rate of infection increases, and when η1 decreases, then the rate of infected population decreases. Moreover, when η2 increases, the rate of infected and hospitalized increases and then decreases when η2 decreases as in Figures 10(a) and 10(b).Figure 8
Variation of contact rateβ in susceptible and infected population.
(a)(b)Figure 9
Variation ofη1 in susceptible and infected human population.
(a)(b)Figure 10
Variation ofη2 in infected and hospitalized human population.
(a)(b)


## 4.1.5. Parameter Identifiability and Model Fitting

The identifiability of parameters is essential to the proposedSEQIHR model. Such parameters are β, μ, ω, α, η1, and η2. Parameter identifiability is implemented by using the least square method to minimize the sum of squared differences between the observations and the SEQIHR model [37] and is defined as
(66)SSθ=∑i=1nyi−fxi,θ2,where yi are the observed data of all compartments, i is the number of compartments (i.e., i=1,2,⋯,n), and fxi,θ is the solution for all compartments of the SEQIHR model. With the initial values of the parameters given as β, μ, ω, α, η1, and η2, the least squares identifiabilities are obtained as shown in Table 4 as the initial value and estimated values and used to fit the simulated data as shown in Figure 11. The relationship of initial parameter values and the identifiable values from the least square method is very close.Table 4
Parameter identifiability.
ParameterInitial valuesEstimated valuesβ0.0150.0147μ8.949.8ω0.0830.096b0.01040.013α0.851.1ν0.070.063θ0.24350.224η10.850.72η20.650.51δ0.0390.038q0.0990.098ρ0.0030.0032λ0.0020.0018ϕ0.0010.0009Figure 11
Fitted graphs.


## 5. Conclusion

COVID-19 pandemic spread rapidly all over the world, which led to the severe human and socioeconomic burden worldwide. In this study, a mathematical model was developed for the transmission of COVID-19 when a human is stressed. The model consists of six compartments: Susceptible (S), Exposed (E), Quarantined (Q), Infectious (I), Hospitalized (H), and Recovered (R) human population. Initially, the model was formulated and some mathematical analyses were presented, including positivity, invariant region, existence, uniqueness of the solution, and stability results for the disease-free equilibrium. The disease-free equilibrium for both local and global is stable when R0<1 was proved. This exploration suggests that the COVID-19 disease can enter and spread to the human population if R0>1 provided that the initial human population is close to the infested region. But also, die out when few initial human populations are infected and R0<1. The basic reproduction number obtained from this study was 2.1692, which shows that the disease is endemic and unique. The most sensitivity indices are summarized in Table 3, and the least positively and negatively sensitive parameters are crucial for the transmission of COVID-19.The numerical results in this study showed that stress affects many quarantined people whereby their immune system is disturbed by the COVID-19, and this is most likely because during quarantine, people are isolated from their families and community members, so they develop fear, and later, the body becomes stressed which affects their immune system. When the immune system is disturbed, it fails to fight against the intruders, which leads to the fast spread of COVID-19 throughout the body. Our graphical presentation illustrated that the control parameters showed a great success on minimizing the spread of COVID-19 in the community.The plan for a future work is to use more detailed and authentic data when having access to COVID-19 data which will be employed in theSEQIHR model. Furthermore, we intend to add a vaccination in our model compartment to implement optimal control strategies and extend to a stochastic model. The limitations of this work are the assumption on an equal death in all compartment while in real situation the infected population has a higher death rate than the susceptible population.



---

*Source: 7772263-2022-01-13.xml*