#### 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 () was derived by using the next-generation matrix method, which shows that, when , the disease-free equilibrium is globally asymptotically stable whereas when 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 . 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.

#### 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 Section 1. 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 model is divided into six human subclasses, namely, susceptible (those who are at risk to contact COVID-19 infection), exposed (the population which is infected but not infectious), quarantined (those who contacted a COVID-19-infected individual but did not develop any symptom), infectious (those who have COVID-19 symptoms and are capable of spreading the disease), hospitalized (infectious individuals admitted to a healthcare facility (active cases)), and recovered (those recovered from the COVID-19).

##### 2.1. Model Assumptions

By presenting an infectious disease with a mathematical model, the following assumptions for the *SEQIHR* 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 Section 2.1, the variable and parameter descriptions are to be presented by the following *SEQIHR* model compartmental diagram:

Let 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:

The total population is given by the mathematical equation: where and .

Table 2 shows the model parameters and their descriptions:

The 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, , 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 .

An individual enters the Exposed class, , after direct contact with an infected person, with an infection rate . Furthermore, taking to hospital decreases the rate of and 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 class is diminished by a leaving rate of natural death .

The Quarantined class, , 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 and diminishing by leaving a natural death rate of .

The Infectious class, , are individuals who progress from exposed at the rate of (1-) and then quarantine at the rate of . Additionally, this class decreases with hospitalized individuals at the rate of and recovered at the rate of . Not only that but also it diminished by leaving the rate of natural death, , and death due to the disease, .

The Hospitalized class, , are individuals who progress from the exposed class at the rate of and , then from the infectious compartment at the rate of . Individuals recovered from hospitalized class at the rate of and diminished by leaving the rate of natural death, , and death due to COVID-19 disease, .

The Recovered class, , 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 .

##### 2.3. Model Equations

Based on the assumptions made and the relationship that exists between the variables shown in Figure 1, the system of six ordinary differential equations is formed as in

#### 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 .

Theorem 1. *Let the initial data set be . The solution set of the model system 2 is positive .*

*Proof. *From the system of model Equation (2), consider the first equation:
By considering the negative term, by ignoring the rest, Equation (8) is reduced to
This is the first-order linear differential inequality which can be solved by a separable method (where ) resulting in
in the absence of COVID-19 disease,

By applying the same procedure to the remaining equations, the results are , , , , and . Therefore, the set of solutions , and of the model is positive .

##### 3.2. Invariant (Boundedness) Region

The *SEQIHR* 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 . The bounded region is obtained through the following theorem.

Theorem 2. *The set is positively invariant and attracts all solutions in .*

*Proof. *Since , then the derivative of is given as
Substituting model Equation (2) to Equation (12) gives
Further simplification leads to
implying that
For disease-free, and , then .

This is a first-order linear differential equation, by a separable method.

When , . When then, .

The Invariant region is given by . 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: . The following questions will be of interest: (i)Under what conditions can we say that a solution to the equation exists?(ii)Under what conditions does a unique solution exist to the equation ?

Consider the following equations to answer the question.

###### 3.3.1. Uniqueness of Solution

Theorem 3. *Let us use D to denote the domain:
*

*Proof. *Suppose satisfies the Lipschitz condition; therefore,
Whenever the points and belong to the domain D and is used to represent the positive constant, then, there exist a constant and a unique solution of system 11 in the interval . It is essential to note that condition (24) is satisfied by to be continuous and bounded in domain D.

If has a continuous partial derivative 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:

Therefore, we look for the bounded solution of the form .

###### 3.3.2. Existence of a Solution

Theorem 4. *Let the domain be denoted by , also defined in Equation (23), such that Equations (24) and (25) hold. Then, the existing solution of Equation (2) is bounded in the domain .*

*Proof. *From Equations (10)–(15), by showing that 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:
Similarly, from Equation (18), we obtain the following system of equations:

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 .

##### 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

When there is no disease, then, =, =, =, =, and =.

By considering each model equation from Equation (2),

In addition to the second model equation,

Similarly, for the third model equation, then,

By considering the same procedures for the fourth, fifth, and sixth model equations, the following is obtained:

Therefore,

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 ()

The basic reproduction number 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, analyzes the equilibrium stability, , 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 , 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 be the rate of new infection in compartment and be the rate of transfer of individuals into compartment by all means other than the epidemic. The important thing is to obtain the disease-free equilibrium point . Thus, the computed matrices F and V which are matrices, where represents the infected classes, defined by: and , where , , are nonnegative, and is a nonsingular -matrix (the matrix with inverse belongs to the class of positive matrices). Since is nonnegative and is a nonsingular matrix, then and are nonnegative. Therefore, the next-generation matrix is computed as defined by [22].

Note that the basic reproduction number is defined as the spectral radius (dominant eigenvalue) of the matrix [23], that is, where is the rate of new infection in compartment . The new forces of infection are

From Equation (36), when and meet, we obtain the following:

From Equation (37), the Jacobian matrix of disease-free equilibrium (DFE) is given by

The partial derivative of Equation (38) with respect to , and is given as

Given that and +. The inverse of is

The product matrix is given by:

Thus, the basic reproduction number becomes

##### 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 and unstable if .*

*Proof. *At the endemic equilibrium point , and . The variables are given as follows:

##### 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 point is locally asymptotically stable if , and it is unstable when .*

*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 variables , which are used to generate the Jacobian matrix as in 21.
At a disease-free equilibrium,
The disease-free equilibrium will be asymptotically stable if the eigenvalues of .
From matrix (46), it is clear that the first, second, and third eigenvalues are
Then matrix (46) reduces to a matrix after the cancellation of the respective rows and columns used to obtain the first, second, and third eigenvalues as shown in
The characteristic polynomial of matrix (48) is given in the form
where
However, , , and , condition .
where

Hence, the condition 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 the *SEQIHR* 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 . Theorem 7 is considered.

Theorem 7. *The DFE is globally asymptotically stable if , and unstable if .*

*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
By considering Equation (54), is the vector of the nontransmitting compartment, is the vector of transmitting compartment, and is the vector of disease-free equilibrium point.
We must show that the matrix has real negative eigenvalues and is a Metzler matrix in order for DFE to be globally stable (i.e., the off-diagonal elements of are nonnegative, symbolically denoted by . We can derive equations with and without transmission from model 1, as stated herewith.
The eigenvalues of the matrix are located at the diagonal (-b and -b), and these eigenvalues are real, distinct, and negative. Moreover, matrices and are given by

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 point is asymptotically stable when and unstable when .*

*Proof. *The logarithmic Lyapunov function is used to analyze the stability of the endemic equilibrium and is given in the form
where represents a positive constant, represents some free virus in compartment , and denotes the number of free viruses in compartment at the equilibrium point. Then, model system (2) is now written as follows:
The constants , and are nonnegative constants and the function which is continuous and differentiable. Consider the derivative with respect to each compartment
At the endemic equilibrium point,
where Equation (63) gives

When , , , , , and . Therefore, or zero and the function is negative when .

By following the approach of [29], the largest invariant set in is a singleton set which is the endemic equilibrium point, and using LaSalle [30] invariant principle, is globally asymptotically stable when 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 . The sensitivity index of each partial basic reproduction number for its parameters.

##### 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 the *SEQIHR* 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 is given by

.

From Table 3, it is observed that the parameter is more sensitive since it increases the basic reproductive number by . 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 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 (); then . Other sensitive parameters are with , and which both have . Parameters, and , have the least values which are -0.0717 and -0.9566, respectively, meaning that decreasing and by a certain percentage always decreases by the same percentage and the same thing happens if we increase . 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 by their respective percentages.

From Figure 2, is positively and highly correlated with as the absolute value of its PRCC value is higher than the corresponding value of other parameters. Furthermore, natural death is highly negatively correlated with .

###### 4.1.1. Dynamic Population Simulation with a *SEQIHR* Model

The numerical simulation of the *SEQIHR* 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.

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

From Section 4.1, we observed that the most positive parameter is . The relationship between and is shown in Figure 4(a), whereby we noticed that the increase of stress results in a quick rise of . 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 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.

**(a)**

**(b)**

###### 4.1.3. Simulation of Stability Analysis of *SEQIHR* 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 variables , and varied by considering Figures 5(a) and 7(b), respectively, are illustrated as follows:

**(a)**

**(b)**

**(a)**

**(b)**

**(a)**

**(b)**

###### 4.1.4. Control Parameters

By considering the variation of some control parameters; from Figure 8(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 () increases, the rate of infection increases, and when decreases, then the rate of infected population decreases. Moreover, when increases, the rate of infected and hospitalized increases and then decreases when decreases as in Figures 10(a) and 10(b).

**(a)**

**(b)**

**(a)**

**(b)**

**(a)**

**(b)**

###### 4.1.5. Parameter Identifiability and Model Fitting

The identifiability of parameters is essential to the proposed *SEQIHR* model. Such parameters are , , , , , and . 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
where are the observed data of all compartments, is the number of compartments (i.e., ), and is the solution for all compartments of the *SEQIHR* model. With the initial values of the parameters given as , , , , , and , 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.

#### 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 (), Exposed (), Quarantined (), Infectious (), Hospitalized (), and Recovered () 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 was proved. This exploration suggests that the COVID-19 disease can enter and spread to the human population if provided that the initial human population is close to the infested region. But also, die out when few initial human populations are infected and . 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 the *SEQIHR* 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.

#### Data Availability

Some data used in our numerical simulation (as shown in Table 3) are assumed and others from published articles as cited in this work.

#### Conflicts of Interest

The authors declare that they have no competing interests.

#### Authors’ Contributions

James Nicodemus Paul is responsible for the conceptualization, model formulation, model analysis, and drafting of the manuscript. Silas Steven Mirau and Isambi Sailon Mbalawata are assigned to the model formulation and supervision. All authors read and approved the final manuscript.

#### Acknowledgments

The authors acknowledge African Institute for Mathematical Sciences (AIMS) for the support.