Optimal control approach of a mathematical modeling with multiple delays of the negative impact of delays in applying preventive precautions against the spread of the COVID-19 pandemic with a case study of Brazil and cost-effectiveness.

As of June 02, 2020. The number of people infected with COVID-19 virus in Brazil was about 529,405, the number of death is 30046, the number of recovered is 211080, and the number is subject to increase. This is due to the delay by a number of countries in general, and Brazil in particular, in taking preventive and proactive measures to limit the spread of the COVID-19 pandemic. So, we propose to study an optimal control approach with delay in state and control variables in our mathematical model proposed by kouidere et al. which describes the dynamics of the transmission of the COVID-19. That the time with delay represent the delay to applying preventive precautions measures. Pontryagin's maximum principle is used to characterize the optimal controls and the optimality system is solved by an iterative method. Finally, some numerical simulations are performed to verify the theoretical analysis using MATLAB.

Introduction

On December 31, 2019, the Chinese authorities announced a dangerous virus outbreak in the Chinese city of Wuhan. The virus was temporarily called the novel Coronavirus nCoV-19, then the International Virus Classification Committee named it (SARS-CoV-2) [19], [20]. This virus belongs to Corona viruses COVID-19 that infect the human respiratory system from seasonal colds. The danger of the novel Coronavirus 2019 lies in the fact that the infected person does not show symptoms of the disease in the incubation period, which ranges from two days to 14 days. The person then shows mild symptoms such as sneezing, a slight rise in temperature, then severe fatigue, a dangerous rise in temperature, difficulty breathing, and failure to breathe that leads to death, especially if the person suffers from a chronic disease such as diabetes, cancer, or HIV disease.

According to the World Health Organization [1], [2], [3]. The novel Coronavirus COVID-19 belongs to the strains of Coronaviruses, which are transmitted from animals to humans and between humans. The most prominent of these is the SARS virus [21], which also appeared for the first time in China and infected about 8098 people and caused the death of 774 people. The source of this virus was bats. Another virus of this strain is the Middle East Syndrome virus (MERS) [4], [5], [6], [7], which infected about 2498 people and resulted in the death of 858 people. The source of this virus was camels.

Compared to the rest of the Coronaviruses, COVID-19 virus is considered the most dangerous and the fastest in spread. Up to June 02, 2020, this virus caused about 6,400,055 infected cases, 378,069 deaths and 2,930,682 recovered cases. Following this huge spread around the world [13], [14], [24]. The WHO declared a state of emergency.

Globally, the number of infected people. In USA, about the infected individuals reached 1,859,772 and the deaths reached 106929. In Brazil, 529,405 people were infected and 30,046 died. In Russia, about 423,741 people were infected and 5037 died. According to the last stats declared on June 02, 2020. The epidemic has moved from China to Europe and to USA and now to South America, which has become the main epidemic area according to the World Health Organization.

On February 26, 2020, Brazilian authorities announced the first case of COVID-19, the first case recorded in South America, of a Brazilian national from Italy. The spread of the COVID-19 virus was initially stable. Then in recent weeks, where the spread of the virus has experienced tremendous acceleration, reaching in June 02, 2020, 3 months after the registration of the first case, approximately 529,405 infected, 30,046 deaths. Brazil has become the second country in the world in terms of infected, after the United States of America.

Where he recorded in State of S ao Paulo about 109,698 infected, followed by State of Rio de Janeiro about 53,388 infected, and in State of Ceará about 48,489 infected, and in State of Amazonas about 41,378 infected.

A large number of mathematical models have been developed to simulate, analyse and understand the Corona virus. In a related research work, Tian-Mu Chen et al. [8] proposed a mathematical model for simulating the phase-based transmissibility of a novel Coronavirus, and Zhi-Qiang Xia et al. [9] Early dynamics of transmission and control of COVID-19: a mathematical modeling study. Sunhwa Choi and Moran Ki proposed [8] Estimating the reproductive number and the outbreak size of Novel Coronavirus disease (COVID-19) using mathematical model in the Republic of Korea [22] Also, many researches have focused on this topic and other related topics [16], [17], [25], [26], [27], [28], [33], [40] . They have found the basic reproduction number and discussed the stability analysis of their model using the basic reproduction number. They found that the model is locally asymptotically stable at disease free equilibrium when and when endemic equilibrium exists and becomes stable. In their research, equilibrium points are also globally asymptotically stable under certain conditions. Besides these, many mathematical models were formulated to investigate the dynamics of infectious diseases with control (For example: [30], [31], [32]).

But they did not take into account delay in applying precautionary and preventive measures which may cause a delay in response to the potential threats.

Delay in applying precautionary and preventive measures, such as:

The delay in schools closure in the affected countries.

The delay in closing airports and ports in a timely manner.

The delay in suspending air, land and navigation travels from and to the affected countries.

The delay in applying precautionary, preventive and strict measures when receiving travelers.

The delay in diagnosis of the infected cases, given that symptoms do not appear in the disease’s incubation period.

The tolerance of others in the passage of passengers through airports in some countries.

The delay in diagnosis of infected cases, especially as symptoms do not appear in the disease’s incubation period.

The delay in imposing quarantine among patients increases the spread of the disease.

Perhaps the delay in applying the aforementioned precautionary measures was the main reason that Western Europe turned into the main focus of the global pandemic. This led to panic as the number of the infected individuals and deaths keep rising. Thousands of infected cases and hundreds of deaths are recorded daily in Italy, Spain, Germany, the United States of America and France. In Italy, on Saturday, 21 March 2020, a record number of deaths reached nearly 800 deaths, and more than 6000 infected. The pandemic focus shifted from China to Western Europe and that led to the ease of transmission and spread of the disease in the neighboring countries and the countries of North Africa.

In addition, most of the previous research has focused on continuous modeling. In this research, the statistical data are collected at the time (day, week, month, and year). Treatment and vaccination of some patients are also given at the same time. Consequently, it is more direct, more practical, and more precise to describe a phenomenon using continuous temporal modeling. Thus, differential equations appear like a more natural way of describing epidemiological models. In addition, numerical solutions are used for the discretionary differential equations, which encourages us to use the differential equations directly. and therefore can be easily implemented by mathematicians.

To make the modeling of this phenomenon more realistic, we consider an optimal control problem governed by a system of differential equations with time to delay. We study an optimal control problem with time to delay in the state and control variable in a mathematical model of kouidere et al. [25] and a time delay representing delay in applying preventive precaution and measures.

The paper is organized as follows: In Section 2, we represent our mathematical model , analysis of a mathematical modeling of transmission of novel coronavirus “COVID-19”. In Section 3, a mathematical modeling with optimal control startegy of transmission of novel coronavirus “COVID-19”. In Section 4, we present a optimal control problem with multiple delays in our proposed model where we give some results concerning the existence of the optimal control and we caracterize the optimal controls using the Pontryagin’s maximum principle. Numerical simulations through MATLAB are given in Section 5. Finally we conclude the paper in Section 6.

Mathematical model

We consider a mathematical model of kouidere et al., that describes the of transmission of novel coronavirus “COVID-19”. We divide the population denoted by into six compartments: susceptible population Brazil the infected without symptoms the infected with symptoms the infected with complications the number of people who have been quarantined in hospitals and recovered

Hence, we present the COVID-2019 mathematical model is governed by the following system of differential equation:where and are the initial state.

With

- : Denote the incidence of susceptible in Brazil.

- : natural mortality

- : The rate of people who were infected with the virus by contact with the infected patients without symptoms

- : The rate of people who were infected with the virus by communicating with the infected patients with symptoms

- : The rate of people become normaly infected with symptoms

- : The rate of people have developed a rapid and dangerous development of the disease due to immunodeficiency, old age or children.

- : The rate of people have severe complications such as pulmonary failure.

- : The rate of people with symptoms of mild virus who have been quarantined.

- : The rate of people with serious complications who have been quarantined.

- : The rate of people who recovered from the virus.

- : Mortality rate due to complications.

- The rate of people who died under quarantine in hospitals.

Fig. 1 shows that, if Brazil does not apply preventive precautions to reduce the spread of the COVID-19 pandemic, the number of COVID-19 infected cases is likely to reach millions.

Fig. 1

Accumulated cases of COVID-19 in Brazil after 4 months.

Model basic properties

Positivity of solutions

If and the solutions and of system (15) are positive for all .

It follows from the first equation of system (15) that where The both sides in last inequality are multiplied by

We obtain then

Integrating this inequality from 0 to gives: then

.

Similarly, we prove that and  □

Boudedness of the solutions

The set positively invariant under system (15) with initial conditions and .

Also, one assumes that:

If we take limit we have.

It implies that the region is a postivily invariant set for the system (15). □

Exictence of solutions

The system (15) that satisfies a given initial condition

has a unique solution.

Let and so the system (15) can be rewritten in the following form:

whereandThe second term on the right-hand side of (8) satisfies

Where thenwhere .

Thus, it follows that the function is uniformly Lipschitz continuous, and the restriction on and we see that a solution of the system exists [18]. □

Stability analysis and sensitivity of the model parameters

In this section, we will study the stability behavior of system (1) at an Disease Free Equilibrium point is denoted by and an Endemic Equilibrium point is denoted by . The first three equations in system (1) are independent of the variables and . Hence, the dynamics of equation system (1) is equivalent to the dynamics of the equation system:

Equilibrium point

The disease free equilibrium

To find the disease free equilibrium point, we equated the right hand side of model (1) to zero, evaluating it at =0 and solving for the noninfected and noncarrier state variables. Therefore, the disease free equilibrium point

The endemic equilibrium

The endemic equilibrium point it occurs when the disease persists in the community. To obtain it, we equate all the model Eq. (1) to zero. then we obtain

where is the basic reproduction number given by:

The basic reproductive number

In our work, the basic reproduction number is defined as the average number of secondary infections produced by an infected individual in a completely susceptible population. To obtain the basic reproduction number, we used the next-generation matrix method formulated in Bentout et al. [34], Bani-Yaghoub et al. [35], Driessche and Watmough [36].

Through the model equations system (2), then by the principle of next generation matrix, we obtained: where

The inverse of V is given by then

Finally, we have

Local stability of disease free equilibrium

The disease free equilibrium point is locally asymptotically stable if and unstable if .

The Jacobian matrix of system (3) at the disease free equilibrium point as follows:

From the jacobian matrix we obtained a characteristic polynomial:

with

we see that the characteristic equation of has an eigenvalue is negative. So, in order to determine the stability of the we discuss the roots of the following equation .

By Routh-Hurwitz criterion, system (2) is locally stable if and .

Obviously we see that a and b to be positive, must be positive, which leads to .

So is stable if and unstable if . □

Numerical simulation

In this section, we present some numerical solutions of system (1) for different values of the parameters. The resolution of system (1) was created using the Gauss-Seidel-like implicit finite-difference method developed by LaSalle [37], presented in Gumel et al. [38] and denoted the GSS1 method. We use the different initial values for each variable of state, and we use the following parameters: we have the Disease free Equilibrium point and . In this case, we obtained the following remarks: Over time, we notice that the number of susceptible people is close to . We also note that the number of the people infected without symptoms and the number of the people infected with symptoms, are close to zero (see Fig. 2 ).

Fig. 2

COVID-19 model with a disease-free equilibrium.

We start by a graphic representation of the COVID-19 disease-free equilibrium and we use the same parameters and different initial values in table1. . From these figures, using the different values of initial variables and we obtained the following remarks: The number of potential individuals increases and approaches the number (see Fig. 2): Also, the number of the asymptomatic infected cases or cases with mild symptoms decreases and converges to zero (see Fig. 2). The number of the infected people with symptoms and carriers of the virus increases at first, after that it decreases and approaches zero (see Fig. 2). Therefore, the solution curves to the equilibrium = when . Hence, model (1) is stable.

According to theorem (3), the Disease free equilibrium point of system (1) is stable on

Also, for the different initial values for each variable of state, and the following parameters: we have the Endemic Equilibrium point and . In this case, we obtained the following remarks: Over time, we notice that the number of susceptible people is close to . We also note that the number of the people infected without symptoms is close to and the number of the people infected with symptoms are close to (see Fig. 3 ).

Fig. 3

COVID-19 model with an epidemic equilibrium.

Therefore the solution curves to the equilibruim point when .

Sensitivity analysis of model parameters

we performed a sensitivity analysis, On the basic parameters. to help us to know the parameters that have a high impact on the breeding number of base .

To achieve a sensitivity analysis of the model (1), we followed the technique described by Karrakchou et al. [39].

This technique develops a formula to obtain the sensitivity index of all the basic parameters, defined as for represents all the basic parameters. The sensitivity index of with respect to the remaining parameters are computed as follows:

are obtained and evaluated at (Table 1 ).

Table 1

The sensitivity of indices.

Parameter symbol	Value	Sensitivity indexes
μ	0.02	−0.0672
α_1	0.8	−0.3273
α_2	0.2	−0.2667
α_3	0.2	−0.4947
β_1	0.35	0.7383
β_2	0.3	0.2613
θ_1	0.1	−0.2473

Their sensitivity indices are in Table 1.

The optimal control problem

Problem statement

As of today June 02, 2020, there is no cure or vaccine for the disease. So, we suggest the following strategies: there are four controls and for that represent consecutively the awareness program through sensitization and prevention, quarantine and health monitoring, diagnosis, monitoring and the last is treatment and psychological support with follow-up at time t. In order to have a realistic and logic model, we need to take in consideration that the movement of controlled individuals from the compartment of Infected without symptoms () to quarantined in hospitals (), and transition by the contact between Infected without symptoms (Respectively Infected with symptoms) with susceptible to susceptible, and the transition of Infected with symptoms (Respectively Infected with complications) to quarantined in hospitals, and the transition of quarantined in hospitals to rocovered is subject to a delay. Thus, the time delay is introduced into the system as follows: at the moment,Thus, the delay is introduced into the system as follows: in time, only a percentage of individuals () that have been treated and controlled time unit ago, that is to say that at the time with are removed to other compartments.

So, the mathematical system with time delay in state and control system of variables is given by the nonlinear retarded system of differential equations:

In addition, for biological reasons, we assume, for that and are nonnegative continuous functions and and .

The set is positively invariant under system (4) with initial conditions, and .

By adding the equations of system (4) we obtain

If we take limit we have.

It implies that the region is a postivily invariant set for the system (4). □

Existence of solutions

The controlled system (4) that satisfies the given initial condition

has a unique solution.

Let and

So the system (4) can be rewritten in the following form:wherewhere and andandThe second term on the right-hand side of (5) satisfiesandwhere and are some positive constants, independent of the state variables and .

Thenwhere .

Thus, it follows that the function is uniformly Lipschitz continuous, From the definition of the controls and and the restriction on and we prove that a solution of the system exists [18]. □

The optimal control: existence and characterization

The problem is to minimize the objective functionalwhere and are the cost coefficients. They are selected to weigh the relative importance of and at time is the final time.

In other words, we seek the optimal controls and such thatwhere is the set of admissible controls defined by

Existence of an optimal control

In this section we introduce a result concerning the existence of optimal control.

Consider the control problem with system (4) .

There exists an optimal control such that

The existence of the optimal control can be obtained using a result by Fleming and Rishel [12], checking the following steps:

The set of controls and corresponding state variables is nonempty. To prove this condition we use a simplified version of an existence result of Boyce and DiPrima ([15], Theorem 7.1.1)

is convex in .

The control set is convex and closed by definition. Take any controls and . then

Additionally, we observe that and

then

Hence, for all and .

The right hand sides of equations of system (4) are continuous, bounded above by linear function in the state and controls variable.

The integrand of the objective functional is clearly convex on .

There exists constants and such that The integrand in the objective functional satisfiesIndeed

.

The state variables being bounded, let and

Therfore, from Fleming and Rishel [12], we conclude that there exists an optimal control.

□

Characterization of the optimal control

In order to derive the necessary conditions for the optimal control, we apply Pontryagin’s maximum principle [11] we have the Hamiltonian at time defined bywhere is the right side of the difference equation of the state variable.

Given the optimal controls and the solutions and of the corresponding state system (4) , there exists adjoint variables and satisfying: where

With the transversality conditions at time : and .

Furthermore, for the optimal controls and are given by

The Hamiltonian is defined as follows:

where:

For the adjoint equations and transversality conditions can be obtained by using Pontryagin’s maximum principle [5], [10], [23] such that

For the optimal controls and can be solved from the optimality condition,

That are

By the bounds in of the controls, it is easy to obtain and are given in (11)–(14) the form of system (4). □

Numerical simulation

In this section, we present the results obtained by numerically solving the optimality system. In our control problem, we have initial conditions for the state variables and terminal conditions for the adjoints. That is, the optimality system is a two-point boundary value problem with separated boundary conditions at times step and . We solve the optimality system by an iterative method with forward solving of the state system followed by backward solving of the adjoint system. We start with an initial guess for the controls at the first iteration and then before the next iteration, we update the controls by using the characterization. We continue until convergence of successive iterates is achieved. A code is written and compiled in Matlab using the following data.

Different simulations can be carried out using various values of parameters. In the present numerical approach, we use the following parameters values taken from [13]:

Table 2 and 3 below we will discuss only the negative impact of delay on the evolution of potential people and recovered people from the disease and their response to precautionary and preventive measures.

Table 2

Parameter values used in numerical simulation.

Paramter	Description	Value in d^−1
μ	Natural mortality	0.02
β_1	The rate of people who were infected by contact with the infected without sympt	0.2
β_2	The rate of people who were infected by contact with the infected with sympt	0.1
α_1	The rate of people become normaly infected with symptoms	0.8
α_2	The rate of people have developed a rapid development of the disease	0.4
α_3	The rate of People have severe complications such as pulmonary failure.	0.2
θ_1	The rate of people with symptoms of mild virus who have been quarantined.	0.1
θ_2	The rate of people with serious complications who have been quarantined	0.2
δ_1	Mortality rate due to complications.	0.3
δ_2	The rate of people who died under quarantine in hospitals.	0.08
σ	The rate of people who recovered from the virus	0.08
Λ	Denote the incidence of susceptible.	2,000,000

Table 3

Population values used in numerical simulation Since control and state functions are on different scales, the weight constant value is chosen as follows: and .

Popoluation	Description	Value
S(0)	Population in Brazil	3×10^8
I_W(0)	Infected without symptoms	3×10^5
I(0)	Infected with symptoms	100,000
C(0)	Infected with complications	10,000
H(0)	Infected who have been quarantined in hospitals	60,000
R(0)	Recovered people	20,000

The proposed control strategy in this work helps to achieve several objectives.

Strategy A: sensitization and prevention

For This strategy we use only the optimal control .

From Fig. 4 and Table 2 and 3, Delays in implementing the strategy to guide, educate and protect against COVID-19 in Brazil had very serious consequences after the outbreak of the epidemic after the initial settlement. This has resulted in a large and recently accelerated spread of more than 370,000 cases of COVID 19 virus, and its results will be a delay in early detection, which reduces the number of people healed and has significantly increased the number of deaths in Brazil.

Fig. 4

The evolution of the number of Recovered and infected in quarantine with and without control and with delay.

This delay is not taking into account advice and measures of prevention in a real and serious way:

Not taking the disease seriously in Brazil.

Not to wearing medical masks while keeping the hand clean, while reducing handshake.

Some consider this a regular virus, just like any other seasonal flu.

Not leaving a safe distance between him and the affected person.

In Fig. 2 shows that the number of recovered from the disease decreased negatively, we notice that after 100 days the number of recovered decreased From to after the delay of 10 days, while the decrease in the number of suscepetible increased from to after the delay of 50 days.

Strategy B: diagonsis and treatment

In his strategy, we combined two optimal controls and .

In Fig. 5 and Tables 2 and 3 showed that the number of recovered from the disease decreased negatively, we noticed that after 100 days the number of recovered decreased From to after the delay of 10 days, while the decreased in the number of suscepetible increased from to after the delay of 20 days. We explained these results by the following: whenever the delay in the application of the strategy of surveillance and diagnosis of expatriates in Brazil from endemic countries such as Italy is a model of the fact that the first infected with Brazilian citizen came from Italy, with a diagnosis of his health and placing it under sanitary control for a period of at least 14 days. This has led to the massive spread of the virus in all regions of Brazil.

Fig. 5

The evolution of the number of Recovered and infected in quarantine with and without control and and with delay.

Strategy C: quarantine and treatment

For This strategy we combined two controls and .

In Fig. 6 and Tables 2 and 3 and showed that the number of recovered from the disease decreased negatively, we notice that after 100 days the number of recovered decreased From to after the delay of 10 days, while the decreased in the number of suscepetible increased from to after the delay of 50 days. We clarify these results by the following: The delay in implementing the quarantine strategy for infected with COVID-19, whether inside hospitals or at home in Brazil, and has led to a spread widespread in all regions of Brazil, where the number of people infected now exceeds 370,000 cases of infection of COVID-19.

Fig. 6

The evolution of the number of Recovered and infected in quarantine with and without control and and with delay.

Strategy D: quarantine and diagnosis and monitoring with treatment

To use this optimal strategy, we combine between the optimal controls and

From Fig. 7 , and depending on the values in Table 2, we note that the longer the delay in serious use and rapid interaction to reduce the spread of the virus. Fig. 4 shows that the number of people recovering from the disease decreased negatively and we notice that after 150 days the number of people recovered decreased From to after the delay of 10 days, while the decrease in the number of people recovered increased from to after the delay of 50 days.

Fig. 7

The evolution of the number of recovered and hospitalization with and without controls with delay.

This decrease is due to several factors, including (Figs. 6 and 7 ):

The delay in schools closure in Brazil.

The delay in closing airports and ports in a timely manner in Brazil.

The delay in suspending air, land and sea travels from Brazil to affected countries.

The delay closure of schools, universities in Brazil.

The delay in applying precautionary, preventive and strict measures when receiving travelers to Brazil.

The tolerance of others in the passage of passengers through airports in some countries

The delay in diagnosis of infected cases, especially as symptoms do not appear in the disease’s incubation period.

The delay in applying total quarantine among patients increases the spread of the disease in all over Brazil.

Cost-effectiveness analysis

In this section, we analyze the profitability of the previous three scenarios and strategies by comparing these two control strategies to determine the most profitable strategy. Following the method applied in several studies [29], [30], [31], we assess costs using the differential cost-effectiveness (ICER). This ratio used compare the differences between costs and health outcomes of two competing intervention strategies.

The ICER is defined as the quotient of the difference in costs in strategies and by the difference in infected averted in strategies and

Given two competing strategies (1) and (2), where strategy (2) has higher effectiveness than strategy (1) ), the values are calculated as follow: where the total costs and the total cases averted are defined, during a given period for strategy i for by: while is the optimal solution associated to the optimal control Using the simulation results and we ranked, in the Table 3 our control strategies in order of increased numbers of averted infections.

Strategy I is compared with strategy J with respect to increased effectiveness, in reference to Table 4. So (Table 4 ):

Table 4

Total costs and total averted infections for strategies 1–4.

Strategy	Total averted infections (TA)	Total cost (TC)
4	5.24×10^7	1.94×10^6
3	4.85×10^7	3.05×10^6
2	4.82×10^7	1.57×10^6
1	4.52×10^7	5.92×10^5

Since then strategy 3 is less effective than strategy 1. Therefore, strategy 2 is excluded from the set of alternatives.

Next, strategy 2 is compared to strategy 4. The ICER values for strategy 4 and strategy 2 are calculated below:

Since then strategy 4 is less effective than strategy 2. Therefore, strategy 2 is excluded from the set of alternatives.

Next, strategy 2 is compared to strategy 1. The ICER values for strategy 2 and strategy 1 are calculated below:

Since then strategy 1 is less effective than strategy 2. Therefore, strategy 2 is excluded from the set of alternatives.

Therefore, the conclusion is that strategy 1 (awarness campaigns to protect potential individuals infected with the virus, prevent contact with people infected with COVID-19 and with hospital quarantine for the infected) is the most effective strategy as previously mentioned by the proportions.

Conclusion

In this paper, we introduced a mathematical model of kouidere et al, that described the spread of COVID-19 virus in Brazil. We divided the population denoted by into six compartments, susceptible population Brazil the infected without symptoms the infected with symptoms the infected with complications the number of people who have been quarantined in hospitals and recovered . We analyzed a mathematical model and the local and the global stability of COVID-19 free equilibrium and COVID-19 endemic equilibrium are obtained. We also studied the sensitivity analysis of the model parameters to know the parameters that have a high impact on the reproduction number . In addition, in order to minimize the number of infected and infected with several complications and number of people in guarantine. We also introduced four controls which, respectively, represent sensitization, prevention, quarantine, diagnosis, monitoring, treatment and psychological support with follow-up. We also studied the optimal control with delay which represents measuring the extent of interaction with the means of treatment or awareness campaigns. that if preventive and proactive measures are implemented, such as awareness-raising and quarantine campaigns in all overe Brazil, the spread of the COVID-19 epidemic will be reduced, thus the number of people infected with the virus and the number of deaths will be reduced. We applied the results of the control theory and we managed to obtain the characterizations of the optimal controls. The numerical simulation of the obtained results showed the effectiveness of the proposed control strategies.

Funding source

None.

Ethical approval

Approval was not required.

CRediT authorship contribution statement

Abdelfatah Kouidere: Writing - original draft, Writing - review & editing. Driss Kada: Writing - original draft, Writing - review & editing. Omar Balatif: Writing - original draft, Writing - review & editing. Mostafa Rachik: Writing - original draft, Writing - review & editing. Mouhcine Naim: .

Declaration of Competing Interest

The authors declare that they have no conflict of interests.