# Set-Valued Control Approach Applied to a COVID-19 Model with Screening and Saturated Treatment Function

Authors: Mohamed Elhia; Lahoucine Boujallal; Meryem Alkama; Omar Balatif; Mostafa Rachik
Journal: Complexity
 (2020)

---

## Abstract

The purpose of this paper is modelling and controlling the spread of COVID-19 disease in Morocco. A nonlinear mathematical model with two subclasses of infectious individuals is proposed. The population is divided into five classes, namely, susceptible (S), exposed (E), undiagnosed infectious (Inc), diagnosed patients (Ic), and removed individuals. To reflect the real dynamic of the COVID-19 transmission in Morocco, the real reported data are used for estimating model parameters. Two controls representing screening effort and limited treatment are considered. Based on viability theory and set-valued analysis, a Lyapunov function is constructed such that both exposed and infected populations are decreased to zero asymptotically. The corresponding controls are derived via a continuous selection of adequately designed feedback map. Numerical simulations are presented with three scenarios (cases when each control is used alone and the case when two controls are combined). Our results show that when only one control is to be applied, screening is the most effective in decreasing the number of people in the three infected compartments, whereas combining both controls is found to be highly effective and leads to a significant improvement in the epidemiological situation of Morocco. To the best of our knowledge, this work is the first one that applies the set-valued approach to a controlled COVID-19 model which agrees with the observed cases in Morocco.


---

## Body

## 1. Introduction

Coronavirus disease 2019, known as COVID-19, is a disease caused by a novel betacoronavirus named SARS-CoV-2 that affects the lower respiratory tract [1]. It is a fatal disease for patients who developed various complications including organ failure, septic shock, and pulmonary oedema (see [1, 2]). Since its appearance in December 2019 in Wuhan, Hubei Province, China, the virus has extended internationally [3], and the number of reported cases has accelerated around the world. This is why, on 20th January 2020, the World Health Organization (WHO) declared COVID-19 as a public health emergency of international concern [4].Given the extent of the COVID-19 disease and the urgency to deal with it, mathematical models are a powerful tool for the management and control of this pandemic. Nowadays, the contribution of the epidemiological models in the fight against infectious diseases is indisputable. Indeed, they provide important information on the dynamics of disease transmission, give an overview of the characteristics of an epidemic, and allow predicting its future evolution and evaluating different intervention strategies to find the best control program. For more details on the use of mathematical models in epidemiology, we refer the interested reader to [5–10].Several mathematical modelling studies have been proposed, which aim to understand the dynamics of this disease and provide appropriate responses to the challenges that it presents. For instance, in [11], a mathematical model considering susceptible, exposed, infected, asymptotic, quarantine/isolation, and recovered classes as in the case of COVID-19 disease was developed. Elasticity and sensitivity analyses are performed, and the global stability for the proposed model is studied. The study in [12] proposed a fractional model to describe the dynamics of COVID-19 and parametrized the model using the available infection cases. In [13], the authors proposed a compartmental model taking hospitalized and asymptomatic cases as extra compartments for COVID-19 disease. The local stability of the disease-free equilibrium in terms of the basic reproduction number is investigated, and the sensitivity analysis of the model with respect to the variation of each one of its parameters is provided. Using a generalized SEIR model and based on the public data, Peng et al. [14] estimated the key epidemic parameters and made predictions on the inflexion point and possible ending time for 5 different regions in China. Based on the collected epidemic data in China, the authors proposed in [15] an extension of the classical SEIR model. They divided the population into 6 different categories and established the time series models based on different mathematical formulas according to the variation law of the original data. Extensions of the classical SIR model by adding time delays were also proposed in [16, 17].To contain the spread of COVID-19, almost all the countries affected by this pandemic have implemented severe control measures such as complete lockdown, social distancing, isolation, early detection, and the implementation of a robust system to trace contacts. In this context, some authors proposed optimal control frameworks for determining the optimal strategies that can mitigate the negative impact of this disease. The study in [18] considered a mathematical model including a quarantine class and governmental intervention measures such as lockdown, media coverage on social distancing, and improvement of public hygiene. To reduce the infected individuals as well as to minimize the cost of implementing government control measures, an optimal control problem is formulated and solved, where the control considered represents the awareness due to media coverage. The authors in [19] considered two control SEIR-type models describing the spread of COVID-19 in a human population. They used one control representing the quarantine intensity and tried to find the optimal solution by applying Pontryagin’s maximum principle. In [20], a model with five compartments was proposed to highlight the role of infected people without symptoms, infected individuals with mild symptoms, and the complicated cases. Then, an optimal control problem was formulated to lessen the infected individuals by carrying out awareness campaigns, diagnosis, surveillance of airports, and quarantine of infected. In [21], the authors proposed stochastic and deterministic models to investigate the transmission mechanism of 2019-nCoV from 15 January to 5 February 2020 in Hubei province. Stochastic and deterministic analyses were performed, and optimal strategies were provided for the deterministic model. For more examples of works that address the problem of COVID-19 control, we refer to [22–28].According to Rismanbaf [29], there is still no vaccine or definitive treatment against the SARS-CoV-2 virus. However, new potential therapies are suggested based on in vitro studies, virtual screenings, and records of their effectiveness on earlier strains of coronavirus, SARS, and MERS. Furthermore, it is well known that lockdown has a colossal economic and social impact; this is why we have to think of alternative solutions. Inspired by the experience of South Korea, which has shown that mass testing is crucial to control COVID-19 disease, a solution to be explored would, therefore, be the extension of screening capacities. In this context, we propose a mathematical model that describes the COVID-19 transmission dynamics, and we investigate the impact of implementing additional treatment protocol and/or the multiplication of the number of screening tests on the spread of this disease in Morocco. Novelty and the main contributions are as follows.Firstly, based on the collected epidemic data and analysis of the actual situation in Morocco, we propose an extension of the classical susceptible-exposed-infectious-removed (SEIR) model. In our model, infectious people are divided into two subclasses depending on whether they are diagnosed positive and isolated or they are still undiagnosed. Considering these two categories (diagnosed and undiagnosed patients) will allow us to emphasize the importance of detection effort to fight against the COVID-19 pandemic. Furthermore, by estimating the model parameters using the real data, we obtained a more realistic model that reflects the dynamic of this disease in Morocco. Note that this model could be adapted easily to other countries.Then, we consider two controls representing mass screening and limited treatment. For the treatment control, we use a saturated function to take into account the constraints on the health offered in this period of crisis such as insufficient human and material resources and the limitations related to the effectiveness of treatment. To investigate the impact of these intervention strategies on the spread of COVID-19, we formulate a control problem which we solve by a new direct approach. It is important to note that most of the control frameworks indicated above (paragraph 4) used an indirect approach based on the maximum principle of Pontryagin [30], which requires to design a suitable objective functional, prove the existence of an optimal solution, and solve the optimality system consisting of the state and adjoint systems that can be hard to solve. Motivated by this, we propose an alternative approach based on viability theory [31] and set-valued analysis [32]. Our approach enables the user to establish explicit formulas for the feedback controls. More precisely, we construct an appropriate Lyapunov function such that both exposed and infectious individuals decreased to zero asymptotically. The expression of our controls is obtained via a continuous selection of adequately designed feedback map.Finally, to show the effectiveness of our theoretical results, we provide numerical simulations with different scenarios. To our knowledge, the present work is the first one that proposes a controlled dynamical model for COVID-19, which is fitted to the existing epidemic data in Morocco, and where the control problem is resolved using a set-valued approach.The remainder of the paper is structured as follows: In Section2, we introduce the compartmental model of COVID-19. Section 3 will be devoted to treat the control problem of COVID-19 by using a set-valued approach. In Section 4, we present numerical simulations with three scenarios in order to show the efficiency of our control approach. Conclusions are given in Section 5.


## 2. COVID-19 Mathematical Model

### 2.1. Model Construction

The model under consideration in the present paper consists of five coupled nonlinear differential equations which model the dynamic of the COVID-19 disease. We consider that the total populationN is divided into five different compartments:(i)
Susceptible individualsS, people who may be infected by the virus(ii)
Exposed individualsE, infected with the virus but without the typical symptoms of infection(iii)
Undiagnosed infectious individualsInc, who are infectious but not yet confirmed by a test(iv)
Diagnosed infectious individualsIc, who are diagnosed as COVID-19-positive patients and are isolated in hospitals(v)
Removed individualsR, people who moved out either by recovery or deathWe assume that the transmission of COVID-19 occurs following adequate contact between a susceptible and undetected infectiousInc. As the infectious people in Ic are isolated, we assume that they have no contact with the susceptible and do not transmit the disease. Due to the nonlinear contact dynamics in the population, we use the incidence function βInc/NS to indicate successful transmission of COVID-19, where β denotes the infection contact rate. We assume that all newly infected individuals enter into the exposed compartment E, where they reside for an average of k−1 days (k is the rate at which individuals leave the latent class by becoming infectious). A proportion p of the infectious individuals are diagnosed and isolated, and they followed a therapeutic protocol and entered the compartment Ic. The remaining infectious patients are considered as free infectious people as they are not yet detected. Among the infectious who are undiagnosed, some of them are diagnosed and isolated at a rate γ. We have observed from the clinical data that there was an increasing change within time for the number of individuals who progress from Ic to class R; to reflect this, we assume that people in Ic enter into the removed class with a time-dependent variable rt, which is derived from the following logistic function:(1)rt≐rf+r0−rf1+expt−t1/2/Δ,where t1/2 represents the time at which rt reaches its half value and Δ determines the width of rt. The parameters r0 and rf model asymptotic values.A flowchart of our model is given in Figure1.Figure 1
Transfer diagram.The dynamics of the model is governed by the following differential equations:(2)S˙=−βSIncN,E˙=βSIncN−kE,I˙nc=1−pkE−γInc,I˙c=pkE+γInc−rtIc,R˙=rtIc,with initial data(3)S0=S0,E0=E0,Inc0=I0nc,Ic0=I0c,R0=R0.


### 2.2. Positivity and Boundedness of Solutions

As model (2) describes the temporal evolution of human populations, we shall show that the solutions remain nonnegative and bounded. From system (2), we have(4)dSdt|S=0=0,dEdt|E=0=βSIncN≥0,dIncdt|Inc=0=1−pkE≥0,dIcdt|Ic=0=pkE+γInc≥0,dRdt|R=0=rIc≥0.We observe that all the aforementioned rates are nonnegative on bounding planes of the nonnegative cone ofℝ5. Thus, if we begin in the interior of this cone, we shall always remain in this cone as the direction of the vector field is inward on all the bounding planes. So, all the solutions of system (2) are positive. Furthermore, as we denote the total population by Nt, Nt=St+Et+Inct+Ict+Rt, and by adding equations in (2), it is straightforward to prove that all solutions S, E, Inc, Ic, and R are bounded by N0. Hence, the model is well posed in a biologically feasible domain given by the following positive invariant set:(5)S=S,E,Inc,Ic,R|0≤S,E,Inc,Ic,R≤N0.


### 2.3. Parameter Estimation

Let us now estimate parameters of model (2), as well as the unknown initial conditions E0 and I0nc. Firstly, recall that the progression rate, k, from E to the infectious stage is supposed to be equal to the inverse of the mean incubation period [10]. So, to calculate k, it suffices to know the average incubation period of COVID-19. The study in [33] showed that the incubation period ranges from 2 to 14 days with 95% confidence and has a mean of around 5 days. Another analysis [34] estimated the median incubation period to be 5.1 days (95% confidence interval, 4.5 to 5.8 days). A result from [35] indicated a mean incubation period of 5.2 days (95% CI, 4.1 to 7.0). Thus, we set the mean incubation period as 5 days and fix k=1/5 per day. For the remaining parameters and initial conditions E0 and I0nc, they will be estimated by fitting the model with the collected data, communicated by the Ministry of Public Health in Morocco (see Table 1), using the MATLAB routine fminsearch. The resulting estimations are listed in Table 2, and the best model fit to the real data is shown in Figure 2. Note that, for the parameters of the function rt, the best fit to real data gives rf=0.2352, r0=0.0096, t1/2=70.3434, and Δ=10.8262.Table 1
Cumulative daily reported data of confirmed and removed cases form March 20, 2020, to May 12, 2020, reported by the Ministry of Public Health in Morocco and available in [36].
March202122232425262728293031Ic8696115143170225275345402479556617R56791113183437394860April123456789101112Ic65470879191910211120118412751374144815451661R6875105125146161183190206229257295April131415161718192021222324Ic176318882024228325642685285530463209344635683758R329343357379416451468493538566607644April252627282930Ic389740654120425243214423R69675485794310961154May123456789101112Ic456947294903505352195408554857115910606362816418R125414291612183220192200236225102647273329993179Table 2
Parameters and initial conditions.
ParameterDescriptionValueUnitReferencesβInfection rate0.3861/dayFittedkProgression rate fromE to the infectious stage1/51/dayFixedpProportion of the infectious that become diagnosed0.6936DimensionlessFittedγProgression rate fromInc to Ic0.05971/dayFittedrtProgression rate fromIc to R0.01−0.241/dayFittedNTotal population35952000Dimensionless[37]S0Initial number of susceptible individuals35951556DimensionlessCalculatedE0Initial number of exposed individuals111DimensionlessFittedInc0Initial number of undiagnosed infectious individuals247DimensionlessFittedIc0Initial number of diagnosed and isolated individuals86DimensionlessTable1R0Initial number of removed individuals5DimensionlessTable1Figure 2
The reported diagnosed and removed individuals and their best fit to model (2). (a) Real data and fitted curve for Ic. (b) Real data and fitted curve for R.
(a)(b)


## 2.1. Model Construction

The model under consideration in the present paper consists of five coupled nonlinear differential equations which model the dynamic of the COVID-19 disease. We consider that the total populationN is divided into five different compartments:(i)
Susceptible individualsS, people who may be infected by the virus(ii)
Exposed individualsE, infected with the virus but without the typical symptoms of infection(iii)
Undiagnosed infectious individualsInc, who are infectious but not yet confirmed by a test(iv)
Diagnosed infectious individualsIc, who are diagnosed as COVID-19-positive patients and are isolated in hospitals(v)
Removed individualsR, people who moved out either by recovery or deathWe assume that the transmission of COVID-19 occurs following adequate contact between a susceptible and undetected infectiousInc. As the infectious people in Ic are isolated, we assume that they have no contact with the susceptible and do not transmit the disease. Due to the nonlinear contact dynamics in the population, we use the incidence function βInc/NS to indicate successful transmission of COVID-19, where β denotes the infection contact rate. We assume that all newly infected individuals enter into the exposed compartment E, where they reside for an average of k−1 days (k is the rate at which individuals leave the latent class by becoming infectious). A proportion p of the infectious individuals are diagnosed and isolated, and they followed a therapeutic protocol and entered the compartment Ic. The remaining infectious patients are considered as free infectious people as they are not yet detected. Among the infectious who are undiagnosed, some of them are diagnosed and isolated at a rate γ. We have observed from the clinical data that there was an increasing change within time for the number of individuals who progress from Ic to class R; to reflect this, we assume that people in Ic enter into the removed class with a time-dependent variable rt, which is derived from the following logistic function:(1)rt≐rf+r0−rf1+expt−t1/2/Δ,where t1/2 represents the time at which rt reaches its half value and Δ determines the width of rt. The parameters r0 and rf model asymptotic values.A flowchart of our model is given in Figure1.Figure 1
Transfer diagram.The dynamics of the model is governed by the following differential equations:(2)S˙=−βSIncN,E˙=βSIncN−kE,I˙nc=1−pkE−γInc,I˙c=pkE+γInc−rtIc,R˙=rtIc,with initial data(3)S0=S0,E0=E0,Inc0=I0nc,Ic0=I0c,R0=R0.


## 2.2. Positivity and Boundedness of Solutions

As model (2) describes the temporal evolution of human populations, we shall show that the solutions remain nonnegative and bounded. From system (2), we have(4)dSdt|S=0=0,dEdt|E=0=βSIncN≥0,dIncdt|Inc=0=1−pkE≥0,dIcdt|Ic=0=pkE+γInc≥0,dRdt|R=0=rIc≥0.We observe that all the aforementioned rates are nonnegative on bounding planes of the nonnegative cone ofℝ5. Thus, if we begin in the interior of this cone, we shall always remain in this cone as the direction of the vector field is inward on all the bounding planes. So, all the solutions of system (2) are positive. Furthermore, as we denote the total population by Nt, Nt=St+Et+Inct+Ict+Rt, and by adding equations in (2), it is straightforward to prove that all solutions S, E, Inc, Ic, and R are bounded by N0. Hence, the model is well posed in a biologically feasible domain given by the following positive invariant set:(5)S=S,E,Inc,Ic,R|0≤S,E,Inc,Ic,R≤N0.


## 2.3. Parameter Estimation

Let us now estimate parameters of model (2), as well as the unknown initial conditions E0 and I0nc. Firstly, recall that the progression rate, k, from E to the infectious stage is supposed to be equal to the inverse of the mean incubation period [10]. So, to calculate k, it suffices to know the average incubation period of COVID-19. The study in [33] showed that the incubation period ranges from 2 to 14 days with 95% confidence and has a mean of around 5 days. Another analysis [34] estimated the median incubation period to be 5.1 days (95% confidence interval, 4.5 to 5.8 days). A result from [35] indicated a mean incubation period of 5.2 days (95% CI, 4.1 to 7.0). Thus, we set the mean incubation period as 5 days and fix k=1/5 per day. For the remaining parameters and initial conditions E0 and I0nc, they will be estimated by fitting the model with the collected data, communicated by the Ministry of Public Health in Morocco (see Table 1), using the MATLAB routine fminsearch. The resulting estimations are listed in Table 2, and the best model fit to the real data is shown in Figure 2. Note that, for the parameters of the function rt, the best fit to real data gives rf=0.2352, r0=0.0096, t1/2=70.3434, and Δ=10.8262.Table 1
Cumulative daily reported data of confirmed and removed cases form March 20, 2020, to May 12, 2020, reported by the Ministry of Public Health in Morocco and available in [36].
March202122232425262728293031Ic8696115143170225275345402479556617R56791113183437394860April123456789101112Ic65470879191910211120118412751374144815451661R6875105125146161183190206229257295April131415161718192021222324Ic176318882024228325642685285530463209344635683758R329343357379416451468493538566607644April252627282930Ic389740654120425243214423R69675485794310961154May123456789101112Ic456947294903505352195408554857115910606362816418R125414291612183220192200236225102647273329993179Table 2
Parameters and initial conditions.
ParameterDescriptionValueUnitReferencesβInfection rate0.3861/dayFittedkProgression rate fromE to the infectious stage1/51/dayFixedpProportion of the infectious that become diagnosed0.6936DimensionlessFittedγProgression rate fromInc to Ic0.05971/dayFittedrtProgression rate fromIc to R0.01−0.241/dayFittedNTotal population35952000Dimensionless[37]S0Initial number of susceptible individuals35951556DimensionlessCalculatedE0Initial number of exposed individuals111DimensionlessFittedInc0Initial number of undiagnosed infectious individuals247DimensionlessFittedIc0Initial number of diagnosed and isolated individuals86DimensionlessTable1R0Initial number of removed individuals5DimensionlessTable1Figure 2
The reported diagnosed and removed individuals and their best fit to model (2). (a) Real data and fitted curve for Ic. (b) Real data and fitted curve for R.
(a)(b)


## 3. Control by a Set-Valued Approach

In this section, we applied a set-valued approach to the control problem of COVID-19. For more details about this approach, we refer the reader to the works in [38, 39], where this approach is used to deal with the problem of asymptotic null-controllability. Note that, in [39], the proposed approach is used for epidemic models with only one infectious compartment, while the main added value here is to control the spread of the disease and decrease the number of infected people in several infectious compartments to zero asymptotically.


### 3.1. Control Problem Statement

In Morocco, as in several countries around the world, government authorities have opted for severe confinement as a strategic measure to limit the spread of COVID-19 in the country. Of course, this strategy is effective to avoid large losses in human lives, but it can lead to serious economic and social consequences, especially if the duration of the confinement is long. Other measures such as mass screening and treatment can also be implemented to limit the damage and allow a safe return of economic activity. To study the impact of mass screening and treatment on the dynamics of COVID-19 in Morocco, we introduce two controlsu1 and u2 into model (2), where(i)
u1 represents all the efforts which allow the multiplication of the number of screening tests such as the widening of diagnostic sites in professional environments in commercial and industrial units, the detection of contact subjects in family clusters, and the increase of the number of laboratories authorized to carry out COVID-19 tests.(ii)
u2 represents the treatment of people confirmed to be COVID-19 positive and who are placed in isolation in hospitals. It is true that there is not yet a specific drug or vaccine for COVID-19, but some studies have shown that protocols combining several drugs have proven their effectiveness (a review on potential treatments for COVID-19 can be found in [29]). It is often assumed that when the resources for treatment are sufficient, the term of treatment is a linear function to infected individuals. However, in the case of COVID-19, as there are constraints and limitations in terms of medical resources, drugs, hospital beds, adequate isolation spaces, and the effectiveness of treatment, it would be more appropriate to consider a saturated treatment function u2Ic/1+dIc, where the control u2 is the treatment rate and d≥0 is the saturation factor that measures the effect of the infected being delayed for treatment. For examples on the use of the saturated treatment function in epidemiological modelling, we refer to [40–42].Our model with controls is given as follows:(6)S˙=−βSIncN,E˙=βSIncN−kE,I˙nc=1−pkE−γ1+u1Inc,I˙c=pkE+γ1+u1Inc−rtIc−u21+dIcIc,R˙=rtIc+u21+dIcIc,with initial data(7)S0=S0,E0=E0,Inc0=I0nc,Ic0=I0c,R0=R0.The control in system (6) is represented by u≐u1,u2 and has values in the subset(8)U≐0,u1max×0,u2max,where u1max and u2max are positive numbers. The aim of this paper is to reduce the number of exposed, undiagnosed, and diagnosed infectious individuals. Then, before we state our control problem, we set and transform the above nonautonomous system to the following autonomous one:(9)x˙=fx,y,u,y˙=hx,y,u,x0=x0,y0=y0,where x≐S,R,t′,y≐E,Inc,Ic′,x0=S0,R0,0,y0=E0,I0nc,I0c,(10)fx,y,u≐−βSIncNrtIc+u21+dIcIc1,hx,y,u≐βSIncN−kE1−pkE−γ1+u1IncpkE+γ1+u1Inc−rtIc−u21+dIcIc.Therefore, the control problem we have to deal with is stated as follows.Problem: for eachx0,y0∈Γ⊆ℝ+6, find a control u¯ such that(11a)u¯:0,∞⟶U,(11b)lims⟶∞y¯s=0,where x¯,y¯ denotes a solution of system (9) for the control u¯.Remark 1.
(i)
It is easy to see that the solution of system (9) remains nonnegative and bounded for all t≥0; therefore, we do not need to impose (in the context of our approach) this constraint on the state of the system.(ii)
It is also noteworthy to mention that the notion of the invariant set coincides with the viability one when there exists a unique solution of system (9). We refer the interested readers to [31].


### 3.2. Set-Valued Approach

Now, we return to problems (11a) and (11b), which we will treat in the context of viability theory and set-valued analysis stated in Appendix. Therefore, in order to establish our control method, first we need to extend the dynamics of COVID-19 model (9) with the control dynamics governed by the following differential equation:(12)u˙=v−αu,where v is the new control that has values in the closed subset(13)V≐0,α1u1max×0,α2u2max,with α1 and α2 being positive numbers. If the initial state of dynamical system (12) belongs to U, i.e., u0=u0∈U and v∈V, then the solution of differential equation (12) also belongs to U, i.e., u¯t∈U for all t≥0, since the solution of the differential equation satisfies the following inequality:(14)u0iexp−αit≤u¯it≤uimax−uimax−ui0exp−αit,for i=1,2. According to Definition A.2, a convenient ℝ+3-Lyapunov function related to problems (11a) and (11b) can be given by(15)φy,z=v,y+δ,z,for ally,z∈ℝ3×ℝ3,where ν and δ denote vectors of ℝ+3 whose both coordinates are positive. Indeed, if a C1 function w:ℝ+⟶ℝ+3 satisfies φws,w˙s≤0, for all s, then the real nonnegative function q=δ1w1+δ2w2+δ3w3 is such that(16)q˙s+infiνiδiqs≤0,for alls≥0.As a result,qs⟶0 at infinity. Thereby, both functions wi do so for all i=1,2,3.Then, condition (11b) can be written as viability of the subset(17)Dφ≐x,y,u∈ℝ8|φy,hx,y,u|≤0,under the augmented control system(18)x˙=fx,y,u,y˙=hx,y,u,u˙=v−αu.Ultimately, consider the feedback mapGφ, for each x,y,u∈Dφ:(19)Gφx,y,u≐v∈V|fx,y,u,hx,y,u,v−αu∈TDφx,y,u.Let(20)π1Dφ≐Ωφ,where π1 denotes the mapping x,y,u⟶x,y. Then, we are ready to show the following result.Theorem 1.
Any continuous selection of the mapGφ provides a solution of problems (11a) and (11b) for each x0,y0∈Ωφ∩ℝ+6.Proof.
For such a selectionζ, it follows that(21)fx,y,u,hx,y,u,ζx,y,u−αu∈TDφx,y,u,for allx,y,u∈Dφ,ζx,y,u∈V.
SinceDφ is closed, f, h, and ζ are continuous; then, by using Lemma A.1, we deduce that subset Dφ is locally viable under system (18). Since both f and h have linear growth and ζ is bounded, subset Dφ is viable under system (18).
Now, letx0,y0 belong to subset Ωφ and u0 be such that x0,y0,u0∈Dφ; then, system (18) has solution x¯,y¯,u¯:0,∞⟶Dφ, which satisfies x¯,y¯,u¯0=x0,y0,u0.
Thanks to equation (17), this solution satisfies(22)φy¯s,y¯˙s≤0,for alls≥0.
By takingx0,y0∈ℝ+6 and using Gronwall inequality, we get y¯s⟶0 as s⟶∞. Furthermore, since ζx,y,u∈V, thanks to inequality (14), u¯s∈U for all s≥0. Ending the proof, problems (11a) and (11b) have the solution from Ωφ∩ℝ+6.
In what follows, we will provide the explicit expression of the feedback mapGφ⋅ of equation (19) in order to get its continuous selection as required by Theorem 1. Now, let us set(23)ϕx,y,u≐φy,hx,y,u.
Then, according to equation (15), we get(24)ϕx,y,u=ν,y+δ,hx,y,u,for allx,y,u∈ℝ8,and let(25)Dφ≐x,y,u∈ℝ8|ϕx,y,u≤0,∇xϕ=δ1βIncN0−δ3r′tIc,∇yϕ=ν1+−δ1+δ21−p+δ3pkν2+δ1βSN+1+u1γδ3−δ2ν3−δ3rt+u21+dIc2,∇uϕ=δ3−δ2γInc−δ3Ic1+dIc.
According to LemmaA.2, we can easily get the expression of the contingent cone TDφ⋅ as follows.Lemma 1.
For eachx,y,u∈Dφ, we have(26)θ,λ,Λ∈TDφx,y,u⟺|θ,λ,Λ∈ℝ8,∇xϕx,y,u,θ+∇yϕx,y,u,λ+∇uϕx,y,u,Λ≤0ifϕx,y,u=0.
Hence, we are ready to determine a useful expression of the feedback mapGφ⋅ given by (19). To that end, we need to define the following functions and map:(27)ℓφ≐∇xϕf+∇yϕh−∇uϕαu,(28)mφ≐∇uϕ,(29)Cφx,y,u≐v∈V|ℓφx,y,u+mφx,y,v≤0,for all x,y,u∈ℝ8. It is noteworthy here that the control affine structure of the system involves that function ℓφ is affine in control u, and function mφ depends only from state x,y. Consequently, we are ready to state the following result.Lemma 2.
We have, for eachx,y,u∈Dφ,(30)Gφx,y,u=|V,ifϕx,y,u<0,Cφx,y,u,ifϕx,y,u=0.Proof.
We can easily get expression (30) by using Lemma 1 and considering equations (19) and (29).
Recall that, in Theorem1, we have proved the case where problems (11a) and (11b) have a solution from Ωφ∩ℝ+6. In the next theorem, we will treat these problems from ℝ+6. Thereafter, we need to introduce the following map, for each σ>0:(31)Cφσx,y,u≐v∈V|ℓφx,y,u+mφx,y,v≤−σ,for each x,y,u∈ℝ8.Theorem 2.
Assume that, for someσ>0, the map Cφσ given by (31) has a continuous selection; then, problems (11a) and (11b) have a solution for each x0,y0∈ℝ+6.Proof.
Letζ be such a selection of the map Cφσ. For each x,y,u∈Dφ, this map has values included in Gφ. Then, the selection ζ is also a selection of Gφ. By Theorem 1, problems (11a) and (11b) have, therefore, a solution for each x0,y0∈Ωφ. Now, we have to show that problems (11a) and (11b) have a solution from ℝ+6∖Ωφ.
Letx0,y0 belong to ℝ+6∖Ωφ. Thereby, ϕx0,y0,u>0, for all u, such that x0,y0,u∈ℝ+8. Let such u0 be given, and we readily have(32)ddsϕx¯s,y¯s,u¯s=∇xϕx¯s,y¯s,u¯s,x¯˙s+∇yϕx¯s,y¯s,u¯s,y¯˙s+∇uϕx¯s,y¯s,u¯s,u¯˙s,which yields, according to equations (27) and (28),(33)ddsϕx¯s,y¯s,u¯s=ℓφx¯s,y¯s,u¯s+mx¯s,y¯s,u¯˙s.
Hence,(34)ϕx¯s1,y¯s1,u¯s1=ϕx0,y0,u0+∫0s1ℓφx¯s,y¯s,u¯s+ζx¯s,y¯s,u¯s,mx¯s,y¯sds.
Sinceζ is also a continuous selection of Cφσ,(35)ϕx¯s1,y¯s1,u¯s1≤ϕx0,y0,u0−σs1.
Thereby,(36)x¯s1,y¯s1,u¯s1∈Dφ,fors1≥ϕx0,y0,u0σ.
Letx1≐x¯s1,y1≐y¯s1, and u1≐u¯s1. As a result, the use of Theorem 1 implies that the system,(37)z˙=fz,p,w,zs1=x1,p˙=hz,p,w,ps1=y1,w˙=ζz,p,w−αw,ws1=u1,has a solution z¯,p¯,w¯:s1,∞⟶ℝ+6×U, which satisfies(38)p¯s⟶0,whens⟶∞.
Consequently, the control given byu¯ on 0,s1 and w¯ on s1,∞ leads to the solution of problems (11a) and (11b) from x0,y0.


## 3.1. Control Problem Statement

In Morocco, as in several countries around the world, government authorities have opted for severe confinement as a strategic measure to limit the spread of COVID-19 in the country. Of course, this strategy is effective to avoid large losses in human lives, but it can lead to serious economic and social consequences, especially if the duration of the confinement is long. Other measures such as mass screening and treatment can also be implemented to limit the damage and allow a safe return of economic activity. To study the impact of mass screening and treatment on the dynamics of COVID-19 in Morocco, we introduce two controlsu1 and u2 into model (2), where(i)
u1 represents all the efforts which allow the multiplication of the number of screening tests such as the widening of diagnostic sites in professional environments in commercial and industrial units, the detection of contact subjects in family clusters, and the increase of the number of laboratories authorized to carry out COVID-19 tests.(ii)
u2 represents the treatment of people confirmed to be COVID-19 positive and who are placed in isolation in hospitals. It is true that there is not yet a specific drug or vaccine for COVID-19, but some studies have shown that protocols combining several drugs have proven their effectiveness (a review on potential treatments for COVID-19 can be found in [29]). It is often assumed that when the resources for treatment are sufficient, the term of treatment is a linear function to infected individuals. However, in the case of COVID-19, as there are constraints and limitations in terms of medical resources, drugs, hospital beds, adequate isolation spaces, and the effectiveness of treatment, it would be more appropriate to consider a saturated treatment function u2Ic/1+dIc, where the control u2 is the treatment rate and d≥0 is the saturation factor that measures the effect of the infected being delayed for treatment. For examples on the use of the saturated treatment function in epidemiological modelling, we refer to [40–42].Our model with controls is given as follows:(6)S˙=−βSIncN,E˙=βSIncN−kE,I˙nc=1−pkE−γ1+u1Inc,I˙c=pkE+γ1+u1Inc−rtIc−u21+dIcIc,R˙=rtIc+u21+dIcIc,with initial data(7)S0=S0,E0=E0,Inc0=I0nc,Ic0=I0c,R0=R0.The control in system (6) is represented by u≐u1,u2 and has values in the subset(8)U≐0,u1max×0,u2max,where u1max and u2max are positive numbers. The aim of this paper is to reduce the number of exposed, undiagnosed, and diagnosed infectious individuals. Then, before we state our control problem, we set and transform the above nonautonomous system to the following autonomous one:(9)x˙=fx,y,u,y˙=hx,y,u,x0=x0,y0=y0,where x≐S,R,t′,y≐E,Inc,Ic′,x0=S0,R0,0,y0=E0,I0nc,I0c,(10)fx,y,u≐−βSIncNrtIc+u21+dIcIc1,hx,y,u≐βSIncN−kE1−pkE−γ1+u1IncpkE+γ1+u1Inc−rtIc−u21+dIcIc.Therefore, the control problem we have to deal with is stated as follows.Problem: for eachx0,y0∈Γ⊆ℝ+6, find a control u¯ such that(11a)u¯:0,∞⟶U,(11b)lims⟶∞y¯s=0,where x¯,y¯ denotes a solution of system (9) for the control u¯.Remark 1.
(i)
It is easy to see that the solution of system (9) remains nonnegative and bounded for all t≥0; therefore, we do not need to impose (in the context of our approach) this constraint on the state of the system.(ii)
It is also noteworthy to mention that the notion of the invariant set coincides with the viability one when there exists a unique solution of system (9). We refer the interested readers to [31].


## 3.2. Set-Valued Approach

Now, we return to problems (11a) and (11b), which we will treat in the context of viability theory and set-valued analysis stated in Appendix. Therefore, in order to establish our control method, first we need to extend the dynamics of COVID-19 model (9) with the control dynamics governed by the following differential equation:(12)u˙=v−αu,where v is the new control that has values in the closed subset(13)V≐0,α1u1max×0,α2u2max,with α1 and α2 being positive numbers. If the initial state of dynamical system (12) belongs to U, i.e., u0=u0∈U and v∈V, then the solution of differential equation (12) also belongs to U, i.e., u¯t∈U for all t≥0, since the solution of the differential equation satisfies the following inequality:(14)u0iexp−αit≤u¯it≤uimax−uimax−ui0exp−αit,for i=1,2. According to Definition A.2, a convenient ℝ+3-Lyapunov function related to problems (11a) and (11b) can be given by(15)φy,z=v,y+δ,z,for ally,z∈ℝ3×ℝ3,where ν and δ denote vectors of ℝ+3 whose both coordinates are positive. Indeed, if a C1 function w:ℝ+⟶ℝ+3 satisfies φws,w˙s≤0, for all s, then the real nonnegative function q=δ1w1+δ2w2+δ3w3 is such that(16)q˙s+infiνiδiqs≤0,for alls≥0.As a result,qs⟶0 at infinity. Thereby, both functions wi do so for all i=1,2,3.Then, condition (11b) can be written as viability of the subset(17)Dφ≐x,y,u∈ℝ8|φy,hx,y,u|≤0,under the augmented control system(18)x˙=fx,y,u,y˙=hx,y,u,u˙=v−αu.Ultimately, consider the feedback mapGφ, for each x,y,u∈Dφ:(19)Gφx,y,u≐v∈V|fx,y,u,hx,y,u,v−αu∈TDφx,y,u.Let(20)π1Dφ≐Ωφ,where π1 denotes the mapping x,y,u⟶x,y. Then, we are ready to show the following result.Theorem 1.
Any continuous selection of the mapGφ provides a solution of problems (11a) and (11b) for each x0,y0∈Ωφ∩ℝ+6.Proof.
For such a selectionζ, it follows that(21)fx,y,u,hx,y,u,ζx,y,u−αu∈TDφx,y,u,for allx,y,u∈Dφ,ζx,y,u∈V.
SinceDφ is closed, f, h, and ζ are continuous; then, by using Lemma A.1, we deduce that subset Dφ is locally viable under system (18). Since both f and h have linear growth and ζ is bounded, subset Dφ is viable under system (18).
Now, letx0,y0 belong to subset Ωφ and u0 be such that x0,y0,u0∈Dφ; then, system (18) has solution x¯,y¯,u¯:0,∞⟶Dφ, which satisfies x¯,y¯,u¯0=x0,y0,u0.
Thanks to equation (17), this solution satisfies(22)φy¯s,y¯˙s≤0,for alls≥0.
By takingx0,y0∈ℝ+6 and using Gronwall inequality, we get y¯s⟶0 as s⟶∞. Furthermore, since ζx,y,u∈V, thanks to inequality (14), u¯s∈U for all s≥0. Ending the proof, problems (11a) and (11b) have the solution from Ωφ∩ℝ+6.
In what follows, we will provide the explicit expression of the feedback mapGφ⋅ of equation (19) in order to get its continuous selection as required by Theorem 1. Now, let us set(23)ϕx,y,u≐φy,hx,y,u.
Then, according to equation (15), we get(24)ϕx,y,u=ν,y+δ,hx,y,u,for allx,y,u∈ℝ8,and let(25)Dφ≐x,y,u∈ℝ8|ϕx,y,u≤0,∇xϕ=δ1βIncN0−δ3r′tIc,∇yϕ=ν1+−δ1+δ21−p+δ3pkν2+δ1βSN+1+u1γδ3−δ2ν3−δ3rt+u21+dIc2,∇uϕ=δ3−δ2γInc−δ3Ic1+dIc.
According to LemmaA.2, we can easily get the expression of the contingent cone TDφ⋅ as follows.Lemma 1.
For eachx,y,u∈Dφ, we have(26)θ,λ,Λ∈TDφx,y,u⟺|θ,λ,Λ∈ℝ8,∇xϕx,y,u,θ+∇yϕx,y,u,λ+∇uϕx,y,u,Λ≤0ifϕx,y,u=0.
Hence, we are ready to determine a useful expression of the feedback mapGφ⋅ given by (19). To that end, we need to define the following functions and map:(27)ℓφ≐∇xϕf+∇yϕh−∇uϕαu,(28)mφ≐∇uϕ,(29)Cφx,y,u≐v∈V|ℓφx,y,u+mφx,y,v≤0,for all x,y,u∈ℝ8. It is noteworthy here that the control affine structure of the system involves that function ℓφ is affine in control u, and function mφ depends only from state x,y. Consequently, we are ready to state the following result.Lemma 2.
We have, for eachx,y,u∈Dφ,(30)Gφx,y,u=|V,ifϕx,y,u<0,Cφx,y,u,ifϕx,y,u=0.Proof.
We can easily get expression (30) by using Lemma 1 and considering equations (19) and (29).
Recall that, in Theorem1, we have proved the case where problems (11a) and (11b) have a solution from Ωφ∩ℝ+6. In the next theorem, we will treat these problems from ℝ+6. Thereafter, we need to introduce the following map, for each σ>0:(31)Cφσx,y,u≐v∈V|ℓφx,y,u+mφx,y,v≤−σ,for each x,y,u∈ℝ8.Theorem 2.
Assume that, for someσ>0, the map Cφσ given by (31) has a continuous selection; then, problems (11a) and (11b) have a solution for each x0,y0∈ℝ+6.Proof.
Letζ be such a selection of the map Cφσ. For each x,y,u∈Dφ, this map has values included in Gφ. Then, the selection ζ is also a selection of Gφ. By Theorem 1, problems (11a) and (11b) have, therefore, a solution for each x0,y0∈Ωφ. Now, we have to show that problems (11a) and (11b) have a solution from ℝ+6∖Ωφ.
Letx0,y0 belong to ℝ+6∖Ωφ. Thereby, ϕx0,y0,u>0, for all u, such that x0,y0,u∈ℝ+8. Let such u0 be given, and we readily have(32)ddsϕx¯s,y¯s,u¯s=∇xϕx¯s,y¯s,u¯s,x¯˙s+∇yϕx¯s,y¯s,u¯s,y¯˙s+∇uϕx¯s,y¯s,u¯s,u¯˙s,which yields, according to equations (27) and (28),(33)ddsϕx¯s,y¯s,u¯s=ℓφx¯s,y¯s,u¯s+mx¯s,y¯s,u¯˙s.
Hence,(34)ϕx¯s1,y¯s1,u¯s1=ϕx0,y0,u0+∫0s1ℓφx¯s,y¯s,u¯s+ζx¯s,y¯s,u¯s,mx¯s,y¯sds.
Sinceζ is also a continuous selection of Cφσ,(35)ϕx¯s1,y¯s1,u¯s1≤ϕx0,y0,u0−σs1.
Thereby,(36)x¯s1,y¯s1,u¯s1∈Dφ,fors1≥ϕx0,y0,u0σ.
Letx1≐x¯s1,y1≐y¯s1, and u1≐u¯s1. As a result, the use of Theorem 1 implies that the system,(37)z˙=fz,p,w,zs1=x1,p˙=hz,p,w,ps1=y1,w˙=ζz,p,w−αw,ws1=u1,has a solution z¯,p¯,w¯:s1,∞⟶ℝ+6×U, which satisfies(38)p¯s⟶0,whens⟶∞.
Consequently, the control given byu¯ on 0,s1 and w¯ on s1,∞ leads to the solution of problems (11a) and (11b) from x0,y0.


## 4. Numerical Simulation

In this section, we provide some numerical simulations using MATLAB to validate the theoretical results obtained in the previous section. To investigate the impact of controlsu1 and u2 on the spread of COVID-19 in Morocco, we will consider three different control scenarios. The first one is when the control u1 (mass screening) is used alone. The second scenario is when only treatment (control u2) is applied. The last one is when controls u1 and u2 are coupled. For our numerical simulations, we consider parameter’s values and initial conditions given in Table 2. The control period is from March 20, 2020, to May 12, 2020 (tf=53 days). It should be noted that in order to give the continuous selection of feedback map Cφσ⋅, for each scenario, as expressed in equation (31), we first derive expressions of the functions ℓφ and mφ given, respectively, by equations (27) and (28). In what follows, r′t represents the time derivative of rt, and its expression is given by(39)r′t=rf−r0expt−t1/2/ΔΔ1+expt−t1/2/Δ2.Scenario 1.
Control with mass screening.
In this case, we assume that only mass screening,u1, is implemented, whereas the control u2 is set to zero. The functions ℓφ and mφ are given by(40)ℓφx,y,u=−δ1βIncN2S−δ3r′tIc+ν1+δ3p+δ21−p−δ1kβSIncN−kE+ν2+δ1βSN+1+u1γδ3−δ21−pkE−1+u1γInc+ν3−δ3rtpkE+1+u1γInc−rtIc−δ3−δ2γIncα1u1,mφx,y=mφ,10,where mφ,1=δ3−δ2γInc. Then, a continuous selection that provides a solution of problems (11a) and (11b) is given by(41)ζ1x,y,u=minα1u1max,max0,−ℓφ+σmφ,1.
In Figures3(a)–3(c), we depict the evolution of the exposed individuals (E) and infectious individuals in compartments (Inc) and (Ic) over time in both cases with and without control. We observe in these figures that mass screening significantly reduces the number of the infected people. At the end of the 53-day period, the number of the infected individuals in E, Inc, and Ic is reduced by 82.7%, 84.3%, and 63.4%, respectively. Figure 3(d) shows the control u1 as a time-dependent function. We can see, in this figure, that the control u1 should be implemented with full effort over 40 days.Figure 3
Number ofE, Inc, and Ic with control u1 and without control. (a) Exposed individuals E. (b) Undiagnosed infectious individuals (Inc). (c) Confirmed and isolated individuals (Ic). (d) Evolution of the control u1 with δ=ν=0.02,0.2,0.1, α1=0.6, σ=1, and u1max=0.9.
(a)(b)(c)(d)Scenario 2.
Control with treatment.
With this scenario, only treatment is employed as a control strategy, i.e.,u1=0. In this case, expressions of the functions ℓφ and mφ are given by(42)ℓφx,y,u=−δ1βIncN2S−δ3r′tIc+ν1+δ3p+δ21−p−δ1kβSIncN−kE+ν2+δ1βSN+γδ3−δ21−pkE−γInc+ν3−δ3rt+u21+dIc2pkE+γInc−rtIc−Ic1+dIcu2+δ3Ic1+dIcα2u2,mφx,y=0mφ,2,where mφ,2=−δ3Ic/1+dIc. Then, a continuous selection that provides a solution of problems (11a) and (11b) is given by(43)ζ2x,y,u=minα2u2max,max0,−ℓφ+σmφ,2.
Figures4(a)–4(c) display the dynamics of the state variables E, Inc, and Ic with and without control u2. It is clear from these figures that the control u2 effectively reduces the number of people who are diagnosed and isolated, while the exposed and the undiagnosed individuals are unaffected by this control. The result in Figure 4(c) indicates that the number of people detected who are isolated decreases with treatment by 19.5% at the end of the control period. Thus, recovery of infective via treatment can reduce the burden on the health system. In addition, from Figure 4(d), one can note that, in order to reduce the number of the isolated individuals, the control u2 should be intensively used during almost the whole control period.Figure 4
Number ofE, Inc, and Ic with control u2 and without control. (a) Exposed individuals E. (b) Undiagnosed infectious individuals (Inc). (c) Confirmed and isolated individuals Ic. (d) Evolution of the control u1 with δ=ν=0.02,0.2,0.1, d=0.02, α2=0.6, σ=1, and u2max=0.9.
(a)(b)(c)(d)Scenario 3.
Control with screening and treatment.
With this approach, both controlsu1 and u2 are used. Expressions of the functions ℓφ and mφ are given by(44)ℓφx,y,u=−δ1βIncN2S−δ3r′tIc+ν1+δ3p+δ21−p−δ1kβSIncN−kE+ν2+δ1βSN+1+u1γδ3−δ21−pkE−1+u1γInc+ν3−δ3rt+u21+dIc2pkE+1+u1γInc−rtIc−Ic1+dIcu2−δ3−δ2γIncα1u1+δ3Ic1+dIcα2u2,mφx,y=δ3−δ2γInc−δ3Ic1+dIc=mφ,1mφ,2.
Then, a continuous selection that provides a solution of problems (11a) and (11b) is given by(45)ζx,y,u=minα1u1max,max0,−ℓφ+σ2mφ,1minα2u2max,max0,−ℓφ+σ2mφ,2.
Figures5(a)–5(c) show a comparison among the number of the infected individuals in the compartments E, Inc, and Ic in the uncontrolled case and when the two controls are combined. We observe, in Figures 5(a) and 5(b), that the curves are almost similar to the case when only the control u1 is applied (Scenario 1). However, we can note from Figure 5(c) that the result obtained for the confirmed and isolated cases (Ic) is more interesting than that in Figures 3(c) and 4(c). Thus, at the end of the 53 days, the number of people in the compartment Ic, when the controls u1 and u2 are used together, is reduced by 81.6% instead of 63.5% in the first scenario and 19.5% in the second scenario. Furthermore, it is important to note that when two controls are activated, the control u1 should be used with a full effort, while the control u2 has to be applied with its maximum level during 36 days with a gradual reduction at the end of the control period (see Figure 5(d)).
Overall, the numerical results under various interventions show that if we have to choose only one control, the screening control is the most useful because it allows a significant decrease in all the infected compartments. In addition, combining the two controls contributes positively to the improvement of the final results, particularly by reducing the number of people in theIc compartment.Figure 5
Number ofE, Inc, and Ic when both controls u1 and u2 are applied and without control. (a) Exposed individuals E. (b) Undiagnosed infectious individuals (Inc). (c) Confirmed and isolated individuals (Ic). (d) Evolution of the controls u1 and u2 with δ=ν=0.02,0.2,0.1, d=0.02, α=0.6,0.6, σ=1, and u1max=u2max=0.9.
(a)(b)(c)(d)


## 5. Concluding Remarks and Discussion

In this work, we were interested in modelling and controlling the COVID-19 pandemic in Morocco. We have proposed an extension of the SEIR model that is adapted to the epidemiological data in Morocco. Thus, in addition to the susceptible (S), exposed (E), and removed (R) compartments, we have divided infectious patients into two classes. The first class (Ic) involves infectious people who are confirmed by a test and placed in isolation in hospitals. The second class (Inc) contains infectious individuals who are not yet tested. This distinction between detected and undetected people seems more realistic than considering a single compartment that includes all the infectious people. It highlights the role of screening in the dynamics of COVID-19 transmission. Also, it helps to differentiate the isolated infectious who no longer contribute to the transmission of the disease from those who are still free and participate positively to the spread of COVID-19. Furthermore, to have a model that reflects as much as possible the epidemiological situation of Morocco, we have estimated the parameters of our model from official data published by the Ministry of Public Health in Morocco from March 20, 2020, to May 12, 2020.To reduce the number of infected people in the compartmentsE, Inc, and Ic, we formulate a control problem based on the proposed model to which we added two controls representing mass screening and limited treatment. This control problem is solved by a new approach based on results from the viability theory which presents the advantage of directly finding the expression of the considered controls. More precisely, we have proposed a set-valued approach based on viability theory and set-valued analysis that allows deriving controls via continuous selections of the designed feedback map.Numerical simulations are carried out for three different scenarios, namely, the application of mass screening only, the adoption of a control strategy based on treatment only, and finally, the combination of these two intervention measures (screening and treatment). The results obtained have shown that when treatment is applied alone, it reduces only the number of people hospitalized, whereas mass screening has a significant effect on the three infected compartmentsE, Inc, and Ic. Our results show that when screening is implemented alone, the number of E+Inc+Ic at the final time decreased by 71%, while there is only a reduction by 19.5% when treatment is used alone. Hence, we conclude that between applying screening only and treatment only, screening is the most effective control intervention. Furthermore, combining screening and treatment is found to be better since the reduction of the total number of the infected at the final time reaches 83%.



---

*Source: 9501028-2020-12-30.xml*