ABSTRACT
We formulate a deterministic epidemic model for the spread of Corona Virus Disease (COVID-19). We have included asymptomatic, quarantine and isolation compartments in the model, as studies have stressed upon the importance of these population groups on the transmission of the disease. We calculate the basic reproduction number and show that for the disease dies out and for the disease is endemic. Using sensitivity analysis we establish that is most sensitive to the rate of quarantine and isolation and that a high level of quarantine needs to be maintained as well as isolation to control the disease. Based on this we devise optimal quarantine and isolation strategies, noting that high levels need to be maintained during the early stages of the outbreak. Using data from the Wuhan outbreak, which has nearly run its course we estimate that which while in agreement with other estimates in the literature is on the lower side.
1. Introduction
Coronavirus disease (COVID-19) [Citation23] is a pandemic that has spread around the globe in the early months of 2020. With the first reported cases in Wuhan, China in late November to early December 2019 [Citation23] the disease is prevalent in at least 185 countries, with more than 2.6 million cases reported and around 176,000 fatalities at the time of writing on 24 April 2020 [Citation24].
The causative agent is a coronavirus SARS-CoV-2 [Citation23]. This belongs to a family of viruses that are found in humans and different species of animals including cattle, camels and bats. Animal coronaviruses have in the past caused serious disease outbreaks in human populations as happened in MERS (caused by MERS-CoV) and SARS (caused by SARS-CoV) [Citation2,Citation10,Citation18], it has been reported that both these viruses like SARS-CoV-2 have their origins in bats.
The virus is mainly spread from person to person, through respiratory droplets, the spread is more likely when people are within 6 feet of each other [Citation2]. This has prompted health authorities to stress upon social distancing as a means to control the spread of the disease. In fact, at this time in the absence of a vaccine and any confirmed anti-viral drug, social distancing, quarantine and isolation seem to be the only viable control strategies available [Citation23]. The virus may be also be transmitted by touching infected surfaces and then touching ones own mouth, nose or eyes hence disinfection of surfaces and hands is also recommended as a measure to help prevent the spread of the disease [Citation23].
Symptoms of the disease may appear 2–14 days after exposure and may include fever, cough shortness of breath, chills, muscle pain and loss of taste or smell [Citation11,Citation23]. Many of these are common with influenza, however high persistent fever, dry cough and difficulty in breathing [Citation11] seem to characterize COVID-19. The symptoms may range from very mild (in 80 % of the cases) to severe (in 15 % of the cases) to critical (in 5 % of the cases) [Citation23]. Those at higher risk for severe illness include the elderly (people ages 65 and above) and people with underlying medical conditions which may include chronic lung diseases, diabetes, chronic kidney diseases, serious heart conditions and immunocompromised individuals [Citation26].
With the sudden and overwhelming outbreak of COVID-19, many mathematical models have been proposed to estimate the growth rates and understand the transmission dynamics of the disease. These range from phenomenological [Citation27] and stochastic models [Citation14], which are very useful in the early stages of the outbreak, to mechanistic models that incorporate our understanding to the transmission pathways. The aim of such modelling is twofold, one to provide estimates of the severity of the outbreak by calculating quantities like the growth trends of the epidemic, estimates of the final outbreak size and duration of the outbreak and second to provide insights into efficacy of various control measures.
Most mechanistic models build upon the basic SEIR model, incorporating features particular to the specific disease being modelled. As mentioned above several such models have been proposed to help understand the transmission dynamics of COVID-19. Kucharski et al. [Citation14] have studied a basic SEIR model to look into transmission dynamics within Wuhan and use estimates to project the disease growth in other areas again using the SEIR framework. Chen et al. [Citation3] have proposed a Model involving dynamics of Bats, Hosts, a Reservoir and Humans, within each group they propose an SEIR or SIR model and study the dynamics. Lin et al. [Citation15] have used an extension of the SEIR model with two additional compartments, one incorporating effects of public perception of risk of COVID-19 and a compartment representing cumulative cases including unreported ones. Bhatnagar et al. [Citation16] consider an SEIR model with Quarantine in their work. Aslan et al. [Citation1] consider an SEIR model with a quarantine compartment for susceptibles and isolation for the infected class in their study. They have used the model to forecast the outbreak in Turkey. The also study the effect of various control strategies. Gumel et al. [Citation7] propose an extension of the SEIR model with an asymptomatic class and a hospitalized class. They then extend this further to two population groups, those using masks and those that are not. Pang et al. [Citation19] use an SEIR model with isolation to estimate for three different stages of the outbreak in Wuhan. Perkins and Espana [Citation20] have proposed an extended SEIR model with an asymptomatic class, a hospitalized class and a vaccinated class. They look at time dependent optimal control strategies depending on different values of .
In this study we propose a mechanistic model for the transmission dynamics of COVID-19. Asymptomatic Individuals have been identified as an important source of disease transmission, quarantine of the healthy population and isolating infected individuals as mentioned are the primary mechanisms adopted to date to ‘flatten the disease cure’[Citation2,Citation23], we have incorporated these important pathways in our model. We establish a threshold quantity, the basic reproductive number , the disease dies out if and is endemic in the population when . We calculate the sensitivity of on various model parameters, identifying the parameters to which is most sensitive, we use this information to suggest strategies for controlling the epidemic using techniques from optimal control theory. We also look at the variation in as these parameters are varied over realistic values, this leads to some interesting observations which support the measures taken by health authorities to control the outbreak in most countries. We also use a simplified model to estimate for the outbreak in Wuhan, as the epidemic has nearly run its course there.
2. Model formulation
We formulate a deterministic compartmental model for the transmission dynamics of COVID-19. We extend the basic SEIR model to include dynamics considered important in the transmission of COVID-19. These include an asymptomatic infected group, susceptible individuals who are under quarantine and infected individuals who have been isolated. The total population is divided into nine mutually exclusive compartments where a fraction of the susceptible individuals are quarantined at rate ε though self quarantine or lock-down measures, the quarantined individuals re-enter the susceptible class at rate ξ. We also assume that susceptibles are exposed to the virus at rate or they can have asymptomatic infection at rate of . Since lock-down prevails in most of the countries, exposed and asymptomatic individuals can quarantine themselves, with the quarantine rates being assumed to be and , respectively. The exposed population and the quarantined-but-exposed population is assumed to become infected at the rates σ and ν while the infected individuals either recover to enter at rate γ or are isolated at the rate τ. Some of the the quarantined exposed individuals are also isolated at rate . After a time the isolated individuals recover. Similarly, asymptomatic and quarantined asymptomatic individuals recover after time and , respectively. We assume that individuals have complete immunity after recovery from the disease. The schematic of the transmission pathways is given in Figure . (1) (1) The governing equations are given below (Equation2(2) (2) ). (2) (2) where λ is the force of infection (3) (3) Here β is the effective contact rate, while , and are the relative infectiousness parameters associated with the ,A and classes. The description of variables and parameters of the model are presented in Tables and
Table 1. Description of the variables of the model.
Table 2. Description of the parameters of the model.
3. Basic properties
The described variables in model (Equation2(2) (2) ) exhibit non-negative solutions for all time , i.e. the time series Solutions of model (Equation2(2) (2) ) for non-negative initial conditions are positive for all time .
Lemma 3.1
For any given non-negative initial conditions, there exists a unique solution , respectively, for all Moreover
Proof is presented in the Appendix.
Lemma 3.2
The closed set: is positively invariant.
Proof is attached in the Appendix.
4. Steady state analysis
4.1. Disease-free equilibrium
The model (Equation2(2) (2) ) achieves a Disease-Free Equilibrium (DFE) when there is no infection induced by the disease. i.e. . Mathematically this is found by equating the right-hand side of (Equation2(2) (2) ) to zero. Let represent the DFE of the model. The local stability of the steady state can be found by computing the through next generation operator method [Citation22]. The associated F and V matrices of above model (Equation2(2) (2) ) are given as where , , , , ,
4.2. The basic reproduction number
The basic reproduction number is defined to be the expected number of secondary infections caused by a single infected in a fully susceptible population. Its is calculated by finding the spectral radius ρ of the matrix [Citation22]. (4) (4)
Stability of DFE
Lemma 4.1
The steady state of the model (Equation2(2) (2) ) is locally-asymptotically stable if and unstable if .
The basic implication of this theorem is that small number of the infected individuals will not lead to the large outbreaks and disease will die out in long run. To certify that the disease extinction is independent of the initially available sub-populations in model (Equation2(2) (2) ), it is essential to show that DFE is globally asymptotically stable (GAS).
Lemma 4.2
The Steady state of model is globally asymptotically stable if and unstable
Proof is accessible in the Appendix.
4.3. Endemic equilibrium
The endemic equilibrium is the steady state of the system (Equation2(2) (2) ) in the presence of the infection i.e. .
Let denote the arbitrary endemic equilibrium of the model (Equation2(2) (2) ). (5) (5) Moreover, Let the force of infection can be updated as with Simplifying the system (Equation2(2) (2) ) at this particular fix point, the endemic equilibrium becomes (6) (6)
Lemma 4.3
The model (Equation2(2) (2) ) attains the unique positive endemic equilibrium whenever .
5. Numerical simulations
The Numerical simulation are carried out with help of the ODE-45 Matlab tool. The parameters values are assumed to reflect the real pandemic. The average deaths by the corona virus is reported to be 3% [Citation25].The parameters are summarized in the Table
It is evident from Figure that the model (Equation2(2) (2) ) attains the DFE and Endemic equilibrium when is less than 1 and more than 1, respectively, validating our theoretical results of Lemmas.
In Figure we study the variation in as isolation and quarantine rates are varied. We note that for realistic and particular isolation rate we need a high quarantine rate in order to keep across a range of realistic contact rates . We can also see this effect in the contour plot as we study the effect on as both the parameters are varied. This provides support to the lock down and strict quarantine measures now in force in most countries to control the outbreak.
6. Sensitivity and uncertainty analysis
To be able to suggest the most efficient way of controlling the disease we need to determine the parameters we can control and to which is most sensitive. Knowing this will also inform us about the parameters that need to be estimated accurately in estimating from data. PRCC is an efficient sensitivity analysis method based on sampling. PRCC assigns a value between to for each parameter. Monotonic relationship is measured between two variables in a reliable manner by eliminating the effects (linear) of all parameters excluding the one of interest [Citation17]. Positive PRCC value indicates a positive correlation of the parameter with the disease maintenance, whereas a negative value indicates a negative correlation with the infectiousness of the disease. The parameters studied are:
From our analysis we determine that is most sensitive to β, ε and τ. This points to the effectiveness of quarantine and isolation in controlling the epidemic. In the next section we consider the optimal quarantine and isolation strategies.
Uncertainty analysis enables us to quantify a degree of confidence on the model parameters by constructing a confidence interval containing 95% of the future estimates when the same assumptions are made and the only noise source is observational error. We use the Latin-hypercube sampling-based method to quantify the uncertainty of as a function of the model parameters. We sample independently each of the parameter distributions, and these samples are then randomly permuted forming the input parameter vectors. These samples are then used to calculate the distribution for values of (Figure )
The uncertainty analysis, as shown in Figure yields an estimated value of with 95% CI .
7. Optimal control
Optimal Control Theory uses Pontryagin's principle to reduce cost by finding optimal strategies for control parameters of the system.
Using this principle we affix an adjoint system of differential equations (with terminal conditions) to the state system to evaluate the optimal control strategy for the state system. State system along with the adjoint system is known as the optimality system. The action of the adjoint system is similar to Lagrange multiplier in a sense that it attaches constraint parameters to the state system in order to minimize maximize the functional. Further details regarding Optimal Control and adjoint system can be found in [Citation13,Citation21].
This paper uses optimal control theory to find quarantine (for susceptible and exposed) and isolation strategies that would minimize the total infected population and while using optimum resources. The numerical simulations were performed in MATLAB using fourth-order forward Runge–Kutta method to solve the state system and fourth-order backward Runge–Kutta method to solve the adjoint system (because of the transversality condition). The controls were updated and convergence criteria were applied similar to what M. Imran et al. did in [Citation12].
7.1. Optimal quarantine and isolation using state system
COVID-19 is a pandemic at the moment and no vaccine or cure is available. Being highly contagious, it has spread across the globe within months. The death toll is rising and the infectious population is increasing at an alarming rate. The only ray of hope is to stop the disease from spreading by enforcing social distancing rules strictly. We study the optimal social distancing rules in this paper to be applied and for how long should they be kept in place for the disease to die out. For this purpose we define our control set U to be The Lebesgue measurable quantities are bounded above by , which depends on the amount of resources available for the implementation of the control strategies. The primary objective of this section is to minimize the cost functional The for i = 1, 2 are cost balancing coefficients associated with the isolation, susceptible quarantine and exposed quarantine parameters. Our objective is to find the optimal values for the control parameters such that
Theorem 7.1
There exist unique optimal controls which minimize J over U. Also, there exist an adjoint system of such that the optimal quarantine (susceptible) and isolation is characterized as: (7) (7) (8) (8) The adjoint system is given as The above adjoint system satisfies the transversality conditions (9) (9)
7.2. Optimal control comparisons and strategies
In this section we compare the results of optimal and constant control. We first compare the cost of infection i.e. J for both strategies and then compare the infected population. Finally, we discuss the optimal strategies for the control parameters.
The optimal control cost for each parameter is less then the constant control at all times. It can be observed that for ε the cost reduction is significant compared to the other parameters which is due to the high sensitivity of parameter ε over the system.
Now we compare the total infected population using optimal and constant controls. We note a the optimal time dependent quarantine strategy, ε, yields a significant drop in infected population when compared to its constant counterparts, similarly using the optimal time dependent isolation strategy τ we see a significant drop in the infected population as compared to constant rates of isolation over time. The overall decreased population is consistent with all optimal strategies.
The optimal strategies for each parameter is evaluated for different values of the contact rate β. As β increases, so does the time for which the maximum control should be applied. For both quarantine an isolation the maximum possible quarantine and isolation rates should be maintained to the first few weeks of the disease which can then be reduced over time as depicted in Figures and .
8. Estimation of basic reproduction number
The estimation of basic reproduction number is based upon the estimation of the parameters that define . It has been carried out using nonlinear least-square method by fitting the data of COVID-19 outbreak from the Hubei province in China. We use a simplified SEIR model (Equation10(10) (10) ) to get this estimate. The data is sourced from 2019 Novel Corona virus COVID-19 (2019-nCoV) Data Repository by Johns Hopkins CSSE [Citation5]. It consists of confirmed cases reported with the number of deaths and recovered people on a particular day ranging from 22 January 2020 to 13 April 2020. On a particular day, the number of infected individuals is found by subtracting the deaths and recovered people from the confirmed cases on that day. Figure represents the sourced data. (10) (10) The recruitment rate Π depends upon the demographic region from where data is collected and μ is the average natural death rate of humans.
The Effective contact rate β is the critical parameter we need in order to estimate . It describes the rate at which susceptible and infected individuals interact with each other and the associated chance of getting infected after this interaction. It is difficult to determine this parameter directly; however, an indirect methodology proposed by Cintrón-Arias et al. [Citation4] is used to estimate this parameter for estimating . (11) (11) Calculation of is carried by estimating the parameters involved using the ordinary least squares (OLS) scheme. Here we have assumed that the observed data have constant variance error distribution. Since the observed data is assumed to have constant variance error distribution. where and are observed values and set of parameters involving fitting model, respectively. we can estimate using least squares (12) (12) The numerical simulations are performed using the lsqcurvefit function from the optimization tool in Matlab. This function uses the trust-region-reflective algorithm to minimize the squared residuals.Also the bootstrap technique [Citation6] is used to find the associated confidence intervals.
Table 3. Estimates.
9. Conclusions
In this study we presented a compartmental ordinary differential equation model to understand the transmission dynamics of COVID-19. We incorporated the role of asymptomatic individuals as well as effects of quarantine and isolation of different population groups. We analysed the model to suggest strategies for effective control of the epidemic. We also estimated for the Wuhan outbreak.
After establishing the existence of a threshold quantity such that the outbreak dies out if and the disease is endemic if , we look at the model parameters to which is most sensitive. These turn out to be the contact rate (as expected) and the quarantine (ε) and isolation (τ) rates. We study the variation of with these parameters, our analysis shows that we need a high rate of quarantine along with isolation of infected individuals in order to control the disease. This lends support to the lock down measures adopted across the globe to help combat the outbreak.
Using optimal control methods we suggest quarantine and isolation strategies. The study shows that at the early stages of the outbreak we need the maximum possible levels of quarantine and isolation, with these rates dropping over time, the reduction depending on the contact rate. Further, our analysis reveals that an optimal control rather than a high constant control is preferable. This means that quarantining and isolating according the proposed strategy (as opposed to a high constant control) has a more favourable impact in controlling the disease while keeping the costs.
In the last section we use a simplified model to estimate for the outbreak in Wuhan. Using the disease incidence data we use a simple least square fit, the estimated value comes out to be . The value is within range of various estimates in the literature, it is towards the lower side as we have used a simplified model and the effects of the strict quarantine and isolation measures by the Wuhan authorities are reflected in a lower contacted rate on which we base our estimate of . A follow up study using more elaborate models to estimate is being done and we will report further findings in the future.
The pandemic is still ongoing and in many countries the epidemic cure is yet to 'turn'. Our study is an attempt to use mathematical models to provide health authorities some insights into the disease dynamics to help them in decision making for effective disease control.
Disclosure statement
No potential conflict of interest was reported by the author(s).
References
- H.I. Aslan, M. Demir, M.M. Wise, and S. Lenhart, Modeling COVID-19: Forecasting and analyzing the dynamics of the outbreak in Hubei and Turkey, preprint. doi:10.1101/2020.04.11.20061952.
- Centers for Disease Control and Prevention, Coronavirus disease 2019 (COVID-19), 2020. Available at https://www.cdc.gov/coronavirus/2019-ncov/ (accessed 25 April 2020).
- T. Chen, J. Rui, Q. Wang, Z. Zhao, J. Cui, and L. Yin, A mathematical model for simulating the phase-based transmissibility of a novel coronavirus, Infect. Dis. Poverty 9(1) (2020). doi:10.1186/s40249-020-00640-3.
- A. Cintrón-Arias, C. Castillo-Cháivez, L.M.A. Bettencourt, A.L. Lloyd, and H.T. Banks, The estimation of the effective reproductive number from disease outbreak data, Math. Biosci. Eng. 6 (2009), pp. 261–282. doi:10.3934/mbe.2009.6.261.
- Data.humdata.org, Novel Coronavirus (COVID-19) Cases Data – Humanitarian Data Exchange, 2020. Available at https://data.humdata.org/dataset/novel-coronavirus-2019-ncov-cases (accessed 7 April 2020).
- B. Efron and R.J. Tibshirani, An Introduction to the Bootstrap, Chapman & Hall/CRC, Boca Raton, FL, 1993.
- S. Eikenberry, M. Mancuso, E. Iboi, T. Phan, K. Eikenberry, Y. Kuang, E. Kostelich, and A. Gumel, To mask or not to mask: Modeling the potential for face mask use by the general public to curtail the COVID-19 pandemic, Infect. Dis. Model. 5 (2020), pp. 293–308.
- N. Ferguson, D. Laydon, G. Nedjati Gilani, N. Imai, K. Ainslie, M. Baguelin, S. Bhatia, A. Boonyasiri, Z. Cucunuba Perez, G. Cuomo-Dannenburg, A. Dighe, I. Dorigatti, H. Fu, K. Gaythorpe, W. Green, A. Hamlet, W. Hinsley, L. Okell, S. Van Elsland, H. Thompson, R. Verity, E. Volz, H. Wang, Y.Wang, P. Walker, C. Walters, P. Winskill, C. Whittaker, C. Donnelly, S. Riley, and A. Ghani, Report 9: Impact of non-pharmaceutical interventions (Npis) to reduce COVID19 mortality and healthcare demand, 2020. [online] Hdl.handle.net. Available at http://hdl.handle.net/10044/1/77482 (accessed 2 May 2020).
- W. Fleming and R. Rishel, Deterministic and Stochastic Optimal Control, Springer-Verlag, Berlin, 1975.
- B. Hu, X. Fi, L. Fa, Bat Origin of human coronaviruses. Virol J 12, 221 (2015). https://doi.org/10.1186/s12985-015-0422-1.
- C. Huang, Y. Wang, X. Li, and L. Ren, Clinical features of patients infected with 2019 novel coronavirus in Wuhan, China Lancet 395 (2020), pp. 497–506. doi: 10.1016/S0140-6736(20)30183-5
- M. Imran, A. Khan, A. Ansari, and S. Shah, Modeling transmission dynamics of Ebola virus disease, Int. J. Biomath. 10(04) (2017), p. 1750057. doi: 10.1142/S1793524517500577
- D. Kirschner, S. Lenhart, and S. Serbin, Optimal control of the chemotherapy of HIV, J. Math. Biol.35(7) (1997), pp. 775–792. doi: 10.1007/s002850050076
- A.J. Kucharski, T.W. Russell, C. Diamond, Y. Liu, J. Edmunds, S. Funk, R.M. Eggo, F. Sun, M. Jit, J.D. Munday, and N. Davies, Early dynamics of transmission and control of COVID-19: A mathematical modelling study, Lancet Infect. Dis. 20 (2020), pp. 30144–4. doi:10.1016/s1473-3099. doi: 10.1016/S1473-3099(20)30144-4
- Q. Lin, S. Zhao, D. Gao, Y. Lou, S. Yang, S.S. Musa, M.H. Wang, Y. Cai, W. Wang, L. Yang, and D. He, A conceptual model for the coronavirus disease 2019 (COVID-19) outbreak in Wuhan, China with individual reaction and governmental action, Int. J. Infect. Dis. 93 (2020), pp. 211–216. doi:10.1016/j.ijid.2020.02.058.
- S. Mandal, T. Bhatnagar, N. Arinaminpathy, A. Agarwal, A. Chowdhury, M. Murhekar, R.R.Gangakhedkar, and S. Sarkar, Prudent public health intervention strategies to control the coronavirus disease 2019 transmission in India: A mathematical model-based approach, Indian J. Med. Res. 151(2) (2020), pp. 190–199. doi:10.4103/ijmr.ijmr_504_20.
- S. Marino, I.B. Hogue, C.J. Ray, and D.E. Kirschner, A methodology for performing global uncertainty and sensitivity analysis in systems biology, J. Theor. Biol 254 (2008), pp. 178–196. doi:10.1016/j.jtbi.2008.04.011.
- Oie.int., Questions and answers on the COVID-19: OIE – world organisation for animal health, 2020. Available at https://www.oie.int/scientific-expertise/specific-information-and-recommendations/questions-and-answers-on-2019novel-coronavirus (accessed 25 April 2020).
- L. Pang, S. Liu, X. Zhang, Transission dynamics and control strategies of COVID-19 in Wihan China. J. Biol. Syst. XX (2020), pp. 1–18. doi: 10.1142/S0218339020500096
- A. Perkins and G. Espana, Optimal control of the COVID-19 pandemic with non-pharmaceutical interventions, preprint. doi:10.1101/2020.04.22.20076018.
- L.S. Pontryagin, V.G. Boltyanskii, The Mathematical Theory of Optimal Processes, John Wiley & Sons New York/London1962.
- P. van den Driessche and J. Watmough, Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission, Math. Biosci. 180 (2002), pp. 29–48. doi: 10.1016/S0025-5564(02)00108-6
- Who.int., Coronavirus, 2020. Available at https://www.who.int/health-topics/coronavirus (accessed 25 April 2020).
- Worldometers. info., Coronavirus update (Live): 2,874,022 cases and 200,790 deaths from COVID-19 virus pandemic – worldometer, 2020. Available at https://www.worldometers.info/coronavirus/ (accessed 25 April 2020).
- Worldometers.info., Coronavirus mortality rate (COVID-19) – worldometer, 2020. Available at https://www.worldometers.info/coronavirus/coronavirus-death-rate/#who-03-03-20 (accessed 19 April 2020).
- J. Yang, Y. Zheng, X. Gou, K. Pu, Z. Chen, Q. Guo, R. Ji, H. Wang, Y. Wang, and Y. Zhou, Prevalence of comorbidities and its effects in coronavirus disease 2019 patient: A systematic review and meta-analysis, Int. J. Infect. Dis. 94 (2020), pp. 91–95. doi: 10.1016/j.ijid.2020.03.017
- S. Zhao, Q. Lin, J. Ran, S.S. Musa, G. Yang, W. Wangh, Y. Loue, D. Gao, L. Yang, D. He, and M.H.Wang, Preliminary estimation of the basic reproduction number of novel coronavirus (2019-nCoV) in China, from 2019 to 2020: A data-driven analysis in the early phase of the outbreak, Int. J. Infect. Dis.92 (2020), pp. 214–217. doi:10.1016/j.ijid.2020.01.050.
Appendix
A.1. Proof of Lemma 3.1
Proof.
Adding the right-hand side of model (Equation2(2) (2) ) results in total (A1) (A1) It follows that (A2) (A2) Therefore solutions exists, and are eventually bounded on every finite time interval.
A.2. Proof of Lemma 3.2
Proof.
From (EquationA1(A1) (A1) ) and (EquationA2(A2) (A2) ), it follows that as time approaches to infinity Hence, the set is positively invariant.
A.3. Proof of Lemma 4.2
Proof.
Using comparison theorem, we can construct the equations of the model (Equation2(2) (2) ) infected compartments as (A3) (A3) where .
Since Lemma 4.1 guaranties that is locally asymptotically stable when or, Equivalently saying that all eigenvalues of F−V matrix have negative real parts [Citation22]. Hence the inequality system (EquationA3(A3) (A3) ) is stable whenever . So by comparison theorem, all compartments approaches to 0 as Therefore substituting these values to model (Equation2(2) (2) ) results in , and . Thus attaining the for . Hence the steady state is globally asymptotically stable GAS.
A.4. Proof of Theorem 7.1
Proof.
The proof follows similar steps in [Citation12]. The integrand of our functional J is convex with respect to . Our solution for the model is bounded above by for all t. The model has Lipschitz property with regard to the state variables. The convex property of integrand of the functional J with respect to the set U ensures a unique optimal control, based on the boundedness of the state solutions, and the Lipschitz property of the state system with respect to the state variables [Citation9]. Hence the existence of optimal control is established.
The adjoint system is obtained through the following evaluations: The optimality condition gives: Since are bounded in U by , respectively, optimal control becomes: The uniqueness of optimality system implies the uniqueness of optimal control.