# Deterministic Epidemic Models for Ebola Infection with Time-Dependent Controls

Authors: Eric Okyere; Johnson De-Graft Ankamah; Anthony Kodzo Hunkpe; Dorcas Mensah
Journal: Discrete Dynamics in Nature and Society
 (2020)

---

## Abstract

In this paper, we have studied epidemiological models for Ebola infection using nonlinear ordinary differential equations and optimal control theory. We considered optimal control analysis of SIR and SEIR models for the deadly Ebola infection using vaccination, treatment, and educational campaign as time-dependent control functions. We have applied indirect methods to study existing deterministic optimal control epidemic models for Ebola virus disease. These methods in optimal control are based on Hamiltonian function and Pontryagin’s maximum principle to construct adjoint equations and optimality systems. The forward-backward sweep numerical scheme with the fourth-order Runge–Kutta method is used to solve the optimality system for the various control strategies. From our numerical illustrations, we can conclude that effective educational campaigns and vaccination of susceptible individuals as well as effective treatments of infected individuals can help reduce the disease transmission.


---

## Body

## 1. Introduction

The re-emergence of the Ebola virus disease in2014−2016 in West Africa has been classified as the largest outbreak since the disease was first discovered in DRC in 1976 [1]. This highly infectious and deadly disease has claimed many lives and caused huge economic burden in the affected West Africa Countries. The recent 2018−2019 outbreak in eastern DRC is a very complex situation due to insecurities which is seriously affecting public health workers response activities [1].Mathematical modeling of epidemics has contributed significantly in understanding the dynamical behaviour and control of infectious diseases [2]. The complex dynamics of Ebola virus disease has attracted the attention of many researchers who are interested in epidemiological modeling (see, e.g., [3–20]).System of nonlinear equations that incorporates optimal control dynamics are key mathematical tools that are used in compartmental modeling to understand the spread of infectious diseases. Sharomi and Malik [21] conducted a comprehensive survey on optimal control modeling of several infectious diseases. Zakary et al. [22] formulated the epidemic model for controlling Ebola outbreak. A mathematical model for Ebola disease with time-dependent controls is proposed and numerically analyzed in [23]. Ahmad et al. [24] proposed and studied optimal control epidemic model for Ebola virus infection. Area et al. [25] introduced and studied the deterministic epidemic model for Ebola disease that incorporates optimal control dynamics with vaccination control. Grigorieva et al. [26] developed and analyzed a controlled dynamical model for Ebola virus infection.Our present work is motivated by indirect methods in optimal control theory. These methods are based on Pontryagin’s maximum principle and Hamiltonian function to construct adjoint equations and optimality systems for optimal control problems. Indirect methods have widely been applied by several authors, see the detailed survey on optimal modeling of infectious diseases by the authors in [21] and references therein. Recently, indirect methods in optimal control have also been used to study the dynamical behaviour and control of Zika virus disease (see, e.g., [27–30]). We are also inspired by the forward-backward numerical scheme with the fourth-order Runge–Kutta method for optimal control problems described in [31].Rachah and Torres [32, 33] formulated and analyzed SIR and SEIR optimal control models for Ebola infection using direct methods. In both papers, they formulated objective functionals with dynamical state constrained equations. In [32], they proposed one optimal control strategy for the SIR model and the same authors in [33] proposed three different optimal control strategies for the SEIR model. After their models formulations, they applied the ACADO solver which is an automatic control and dynamic optimization tool to perform their numerical simulations. In our present work, we will apply indirect methods to study the mathematical models proposed in [32, 33]. In this study, we will also propose one additional optimal control strategy for the SIR model. We will analytically construct Hamiltonian function and optimality system for the various optimal control strategies. We will then apply the forward-backward sweep method with the fourth-order Runge–Kutta method to perform numerical simulations in Matlab.The rest of the paper is organized in three (3) more sections. In Section2, we will introduce the classic SIR epidemic model to describe the dynamical behaviour of Ebola infection. We will then consider optimal control problems for the SIR model with two control strategies. In Section 3, we will present the basic SEIR epidemic model to describe the transmission dynamics of Ebola virus disease. We will further consider SEIR optimal control problems for this infectious disease with three different control strategies. In all these sections, we will formulate Hamiltonian functions and then apply Pontryagin’s maximum principle to construct adjoint equations and optimality system for the various optimal control strategies. We will also perform numerical simulations for the optimality systems. Finally, we will present the conclusion in Section 4.


## 2. SIR Model

In this section, we introduce the classic SIR epidemic model to describe the dynamical behaviour of Ebola infection. The model assumes constant population size with no vital dynamics [2]. The population is divided into three different classes with St,It, and Rt representing Susceptible, Infected, and Recovered individuals, respectively. As in [32], the nonlinear dynamical system describing Ebola infection is given by(1)dStdt=−νStItN,S0=S0≥0,dItdt=νStItN−δIt,I0=I0≥0,dRtdt=δIt,R0=R0≥0,with St+It+Rt=N, where ν is the infection rate and δ is the recovery rate.In this study, we will work with proportional quantities instead of the actual populations by scaling each state variable (S,I,R) by the total population. For this purpose, we introduce new state variables (s,i,r), which are expressed in terms of the original state variables and the total population given as follows:(2)st=StN,it=ItN,rt=RtN.Differentiating the scaling equation (2) with respect to time t and using the unscaled model problem (1), we obtain the scaled SIR model describing the transmission dynamics of Ebola as follows:(3)dstdt=−νstit,s0=s0≥0,ditdt=νstit−δit,i0=i0≥0,drtdt=δit,r0=r0≥0,with st+it+rt=1, where the new state variables st,it, and rt represent proportions of Susceptible, Infected, and Recovered individuals, respectively.As in [32], we will use this scaled SIR model (3) in Sections 2.1 and 2.2 to formulate optimal control models for Ebola infection with two different optimal control strategies.


### 2.1. SIR Model with Optimal Control Strategy 1

In this section, we consider a controlled dynamical system for the scaled SIR model (3), using vaccination as time-dependent control function, ηt. For this optimal control problem, we want to reduce the number of infected individuals and the cost of vaccination. Therefore, in this control strategy, we minimize the objective functional Jη given by(4)Jη=∫0tfit+B2η2tdt,subject to(5)dstdt=−νstit−ηtst,s0=s0≥0,ditdt=νstst−δit,i0=i0≥0,drtdt=δit+ηtst,r0=r0≥0.where the control set is given as follows:(6)G1=η:ηtis lebesgue measurable,0≤ηt≤1,t∈0,tf,and the positive parameter B is the weight on cost of vaccination.The Hamiltonian function corresponding to the objective functional (4) and the dynamic constrained state system (5) is given by(7)H=it+B2η2t+φ1−νstit−ηtst+φ2νstit−δit+φ3δit+ηtst.We then apply the Maximum Principle proposed by authors in [34] to determine an optimal solution as follows.Suppose thatw,η is an optimal solution for a controlled dynamical system, then there exist adjoint vector function φ=φ1,φ2,…,φn, which satisfy the system below:(8)dwdt=∂Ht,w,η,φ∂φ,0=∂Ht,w,η,φ∂η,dφdt=−∂Ht,w,η,φ∂w.Following the constructed Hamiltonian function (7) and equation (8), we present the adjoint equations and the control characterisation as follows.Theorem 1.
Letη∗ be an optimal control and optimal state solutions s∗,i∗,andr∗ of the corresponding controlled dynamical system (4)-(5) that minimize Jη over G1. Then, there exist adjoint variables φi for i=1,2,3, which satisfy the system below:(9)dφ1dt=νi∗tφ1−φ2+η∗tφ1−φ3,dφ2dt=−1+νs∗tφ1−φ2+δφ2−φ3,dφ3dt=0,with transversality conditions φ1tf=0,φ2tf=0,andφ3tf=0, and the control η∗ satisfies the optimality condition:(10)η∗t=minmax0,s∗tBφ1−φ3,1.Proof.
To derive the system of equations describing the adjoint variables and transversality conditions, we apply Pontryagin’s Maximum Principle and the Hamiltonian function (7) as follows:(11)dφ1dt=−∂H∂s,dφ2dt=−∂H∂i,dφ3dt=−∂H∂r,with(12)φitf=0,i=1,2,3.
By using the property of the control spaceG1 and solving the differential equation given by(13)∂H∂η=0,on the interior of the control set, the characterisation of the optimal control is obtained as given by equation (10).


### 2.2. SIR Model with Optimal Control Strategy 2

In this section, we present a controlled dynamical system by incorporating treatment control,η1t, and educational campaign control, η2t, into the scaled SIR model (3). For this optimal control problem, we want to reduce the number of infected individuals, the cost of educational campaign, and cost of treatment. Therefore, in this control strategy, we minimize the objective functional given as(14)Jη1,η2=∫0tfC1it+C22η12t+C32η22tdt,which subjects to(15)dstdt=−νstit−η2tst,s0=s0≥0,ditdt=νstit−δit−η1tit,i0=i0≥0,drtdt=δit+η1tit+η2tst,r0=r0≥0,where the control set is given by(16)G2=η1,η2:ηitis lebesgue measurable,0≤ηit≤1,t∈0,tf,fori=1,2,and C1 is a positive constant to keep a balance in the size of it. The positive parameters C2 and C3 are the weight on cost for treatments and educational campaigns.The Hamiltonian function corresponding to equations (14) and (15) is given as(17)H=C1it+C22η12t+C32η22t+φ1−νstit−η2tst+φ2νstit−δit−η1tit+φ3δit+η1tit+η2tst.We then apply the Maximum Principle proposed by authors in [34] to determine an optimal solution.Assume thatw,η is an optimal solution for a controlled dynamical system; then, there exist adjoint vector function φ=φ1,φ2,…,φn, which satisfy the following equations:(18)dwdt=∂Ht,w,η,φ∂φ,0=∂Ht,w,η,φ∂η,dφdt=−∂Ht,w,η,φ∂w.Using the formulated Hamiltonian function (17) and equation (18), the adjoint equations and the control characterisation are presented in the following theorem.Theorem 2.
Letη1∗ and η2∗ be optimal control pair and optimal state solutions s∗,i∗,andr∗ of the corresponding controlled dynamical systems (14)-(15) that minimize Jη1,η2 over G2. Then, there exist adjoint variables φi for i=1,2,3, satisfying(19)dφ1dt=νi∗tφ1−φ2+η2∗tφ1−φ3,dφ2dt=−C1+νs∗tφ1−φ2+δφ2−φ3+η1∗tφ2−φ3,dφ3dt=0,with transversality conditions φ1tf=0,φ2tf=0,andφ3tf=0 and the control pair η1∗ and η2∗ satisfies the optimality conditions:(20)η1∗t=minmax0,i∗tC2φ2−φ3,1,η2∗t=minmax0,s∗tC3φ1−φ3,1.Remark 1.
We remark that the proof of Theorem2 and subsequent theorems in the next section of this study are similar to that of Theorem 1, and we therefore omit their proofs.


### 2.3. Numerical Simulations and Discussion

This section deals with numerical solutions of optimal control problems formulated in Sections2.1 and 2.2. Analytical solution of the nonlinear system of differential equations with optimal control is a very hard task in mathematical modeling. Therefore, in this study, we apply the efficient forward-backward sweep numerical scheme with the fourth-order Runge–Kutta method to solve the optimality systems formulated for the two optimal control strategies. This numerical scheme has extensively been described by the authors in [31], in their mathematical modeling textbook, which concerns with the applications of differential equations and optimal control theory. For our numerical simulations, we have adapted the same model parameter values (ν=0.2andδ=0.1) as in the work by the authors in [32] and initial conditions: s0=0.95,i0=0.05,andr0=0. For the positive parameters in the objectives functionals, we have assumed that B=5,C1=1,C2=5, and C3=5. In Figure 1, it is clear that there is a significant decrease in the number of susceptible individuals with control strategies 1 and 2 than without control. Figure 2 shows solution paths for infected individuals with control strategies 1 and 2 and without control. This plot shows significant decrease in the infected individuals with control strategies 1 and 2 than without control. In Figure 3, there is a rapid increase in recovered individuals with control strategies 1 and 2 than without control. Figures 4 and 5 represent control functions for optimal control strategies 1 and 2, respectively.Figure 1
Solution paths for Susceptible individuals with two control strategies and without control,ν=0.2 and δ=0.1.Figure 2
Solution paths for Infected individuals with two control strategies and without control,ν=0.2 and δ=0.1.Figure 3
Solution paths for Recovered individuals with two control strategies and without control,ν=0.2 and δ=0.1.Figure 4
Optimal control functionηt for strategy 1.Figure 5
Optimal control functionsη1t and η2t for strategy 2.


## 2.1. SIR Model with Optimal Control Strategy 1

In this section, we consider a controlled dynamical system for the scaled SIR model (3), using vaccination as time-dependent control function, ηt. For this optimal control problem, we want to reduce the number of infected individuals and the cost of vaccination. Therefore, in this control strategy, we minimize the objective functional Jη given by(4)Jη=∫0tfit+B2η2tdt,subject to(5)dstdt=−νstit−ηtst,s0=s0≥0,ditdt=νstst−δit,i0=i0≥0,drtdt=δit+ηtst,r0=r0≥0.where the control set is given as follows:(6)G1=η:ηtis lebesgue measurable,0≤ηt≤1,t∈0,tf,and the positive parameter B is the weight on cost of vaccination.The Hamiltonian function corresponding to the objective functional (4) and the dynamic constrained state system (5) is given by(7)H=it+B2η2t+φ1−νstit−ηtst+φ2νstit−δit+φ3δit+ηtst.We then apply the Maximum Principle proposed by authors in [34] to determine an optimal solution as follows.Suppose thatw,η is an optimal solution for a controlled dynamical system, then there exist adjoint vector function φ=φ1,φ2,…,φn, which satisfy the system below:(8)dwdt=∂Ht,w,η,φ∂φ,0=∂Ht,w,η,φ∂η,dφdt=−∂Ht,w,η,φ∂w.Following the constructed Hamiltonian function (7) and equation (8), we present the adjoint equations and the control characterisation as follows.Theorem 1.
Letη∗ be an optimal control and optimal state solutions s∗,i∗,andr∗ of the corresponding controlled dynamical system (4)-(5) that minimize Jη over G1. Then, there exist adjoint variables φi for i=1,2,3, which satisfy the system below:(9)dφ1dt=νi∗tφ1−φ2+η∗tφ1−φ3,dφ2dt=−1+νs∗tφ1−φ2+δφ2−φ3,dφ3dt=0,with transversality conditions φ1tf=0,φ2tf=0,andφ3tf=0, and the control η∗ satisfies the optimality condition:(10)η∗t=minmax0,s∗tBφ1−φ3,1.Proof.
To derive the system of equations describing the adjoint variables and transversality conditions, we apply Pontryagin’s Maximum Principle and the Hamiltonian function (7) as follows:(11)dφ1dt=−∂H∂s,dφ2dt=−∂H∂i,dφ3dt=−∂H∂r,with(12)φitf=0,i=1,2,3.
By using the property of the control spaceG1 and solving the differential equation given by(13)∂H∂η=0,on the interior of the control set, the characterisation of the optimal control is obtained as given by equation (10).


## 2.2. SIR Model with Optimal Control Strategy 2

In this section, we present a controlled dynamical system by incorporating treatment control,η1t, and educational campaign control, η2t, into the scaled SIR model (3). For this optimal control problem, we want to reduce the number of infected individuals, the cost of educational campaign, and cost of treatment. Therefore, in this control strategy, we minimize the objective functional given as(14)Jη1,η2=∫0tfC1it+C22η12t+C32η22tdt,which subjects to(15)dstdt=−νstit−η2tst,s0=s0≥0,ditdt=νstit−δit−η1tit,i0=i0≥0,drtdt=δit+η1tit+η2tst,r0=r0≥0,where the control set is given by(16)G2=η1,η2:ηitis lebesgue measurable,0≤ηit≤1,t∈0,tf,fori=1,2,and C1 is a positive constant to keep a balance in the size of it. The positive parameters C2 and C3 are the weight on cost for treatments and educational campaigns.The Hamiltonian function corresponding to equations (14) and (15) is given as(17)H=C1it+C22η12t+C32η22t+φ1−νstit−η2tst+φ2νstit−δit−η1tit+φ3δit+η1tit+η2tst.We then apply the Maximum Principle proposed by authors in [34] to determine an optimal solution.Assume thatw,η is an optimal solution for a controlled dynamical system; then, there exist adjoint vector function φ=φ1,φ2,…,φn, which satisfy the following equations:(18)dwdt=∂Ht,w,η,φ∂φ,0=∂Ht,w,η,φ∂η,dφdt=−∂Ht,w,η,φ∂w.Using the formulated Hamiltonian function (17) and equation (18), the adjoint equations and the control characterisation are presented in the following theorem.Theorem 2.
Letη1∗ and η2∗ be optimal control pair and optimal state solutions s∗,i∗,andr∗ of the corresponding controlled dynamical systems (14)-(15) that minimize Jη1,η2 over G2. Then, there exist adjoint variables φi for i=1,2,3, satisfying(19)dφ1dt=νi∗tφ1−φ2+η2∗tφ1−φ3,dφ2dt=−C1+νs∗tφ1−φ2+δφ2−φ3+η1∗tφ2−φ3,dφ3dt=0,with transversality conditions φ1tf=0,φ2tf=0,andφ3tf=0 and the control pair η1∗ and η2∗ satisfies the optimality conditions:(20)η1∗t=minmax0,i∗tC2φ2−φ3,1,η2∗t=minmax0,s∗tC3φ1−φ3,1.Remark 1.
We remark that the proof of Theorem2 and subsequent theorems in the next section of this study are similar to that of Theorem 1, and we therefore omit their proofs.


## 2.3. Numerical Simulations and Discussion

This section deals with numerical solutions of optimal control problems formulated in Sections2.1 and 2.2. Analytical solution of the nonlinear system of differential equations with optimal control is a very hard task in mathematical modeling. Therefore, in this study, we apply the efficient forward-backward sweep numerical scheme with the fourth-order Runge–Kutta method to solve the optimality systems formulated for the two optimal control strategies. This numerical scheme has extensively been described by the authors in [31], in their mathematical modeling textbook, which concerns with the applications of differential equations and optimal control theory. For our numerical simulations, we have adapted the same model parameter values (ν=0.2andδ=0.1) as in the work by the authors in [32] and initial conditions: s0=0.95,i0=0.05,andr0=0. For the positive parameters in the objectives functionals, we have assumed that B=5,C1=1,C2=5, and C3=5. In Figure 1, it is clear that there is a significant decrease in the number of susceptible individuals with control strategies 1 and 2 than without control. Figure 2 shows solution paths for infected individuals with control strategies 1 and 2 and without control. This plot shows significant decrease in the infected individuals with control strategies 1 and 2 than without control. In Figure 3, there is a rapid increase in recovered individuals with control strategies 1 and 2 than without control. Figures 4 and 5 represent control functions for optimal control strategies 1 and 2, respectively.Figure 1
Solution paths for Susceptible individuals with two control strategies and without control,ν=0.2 and δ=0.1.Figure 2
Solution paths for Infected individuals with two control strategies and without control,ν=0.2 and δ=0.1.Figure 3
Solution paths for Recovered individuals with two control strategies and without control,ν=0.2 and δ=0.1.Figure 4
Optimal control functionηt for strategy 1.Figure 5
Optimal control functionsη1t and η2t for strategy 2.


## 3. SEIR Model

In this section, we introduce the basic SEIR mathematical model to describe the transmission dynamics of Ebola infection. This model also assumes constant population size with no vital dynamics (birth and death rates). The total population is divided into four different classes withSt,Et,Et, and Rt representing Susceptible, Exposed, Infectious, and Recovered individuals, respectively. Following the compartmental modeling concepts in [2], the nonlinear dynamical system describing Ebola virus disease is given by(21)dStdt=−νStItN,S0=S0≥0,dEtdt=νStItN−ρEt,E0=E0≥0,dItdt=ρEt−δIt,I0=I0≥0,dRtdt=δIt,R0=R0≥0,with St+Et+It+Rt=N, where ν is the transmission rate, ρ represent the infectious rate, and δ is the recovery rate.As we carried out in Section2, we will work with proportional quantities instead of the actual populations by scaling each state variable (St,Et,It,Rt) by the total population. We therefore introduce new state variables (st,et,it,rt), which are expressed in terms of the original state variables, and the total population is given by(22)st=StN,et=EtN,it=ItN,rt=RtN.By differentiating the scaling equation (22) with respect to time t and using the unscaled model problem (21), we obtain the scaled SEIR model as follows:(23)dstdt=−νstit,s0=s0≥0,detdt=νstit−ρet,e0=e0≥0,ditdt=ρet−δit,i0=i0≥0,drrdt=δit,r0=r0≥0,with st+et+it+rt=1, where the new state variables st,et,it, and rt represent proportions of Susceptible, Infected, and Recovered individuals, respectively.As in [33], we will use this scaled SEIR model problem (23) in Sections 3.1, 3.2, and 3.3 to construct optimal control models for Ebola infection with three different control strategies.


### 3.1. SEIR Model with Optimal Control Strategy 1

In this section, we consider a controlled dynamical system for the scaled SEIR model (23), using vaccination as time-dependent control function, ηt. As we did in Section 2.1, our main objective for this optimal control problem is that we want to reduce the number of infected individuals and the cost of vaccination. Therefore, in this control strategy, we minimize the objective functional Jη given by(24)Jη=∫0tfit+D2η2tdt,which subjects to(25)dstdt=−νstit−ηtst,s0=s0≥0,detdt=νstit−ρet,e0=e0≥0,ditdt=ρet−δit,i0=i0≥0,drtdt=δit+ηtst,r0=r0≥0,where the control set is given by(26)G3=η:ηtis lebesgue measurable,0≤ηt≤1,t∈0,tf,and the positive constant D is the weight on the cost of vaccination.The Hamiltonian functionH is then given as follows:(27)H=it+D2η2t+φ1−νstit−ηtst+φ2νstit−ρet+φ3ρet−δit+φ4δit+ηtst.We then apply the Maximum Principle proposed by authors in [34] to determine an optimal solution.Given thatw,η is an optimal solution for a controlled dynamical system, there exist adjoint vector function φ=φ1,φ2,…,φn, which satisfy the following equations:(28)dwdt=∂Ht,w,η,φ∂φ,0=∂Ht,w,η,φ∂η,dφdt=−∂Ht,w,η,φ∂w.Using the formulated Hamiltonian function (27) and equation (28), the adjoint equations and the control characterisation are presented in the following theorem.Theorem 3.
Letη∗ be an optimal control and optimal state solutions s∗,e∗,i∗,andr∗ of the corresponding controlled dynamical system (24)-(25) that minimize Jη over G3. Then, there exist adjoint variables φi for i=1,2,3,4, satisfying(29)dφ1dt=νi∗tφ1−φ2+η∗tφ1−φ4,dφ2dt=ρφ2−φ3,dφ3dt=−1+νs∗tφ1−φ2+δφ3−φ4,dφ4dt=0,with transversality conditions φ1tf=0,φ2tf=0,φ3tf=0,andφ4tf=0 and the control η∗t satisfies the optimality condition:(30)η∗t=minmax0,s∗tDφ1−φ4,1.


### 3.2. SEIR Model with Optimal Control Strategy 2

In this section, we present a controlled dynamical system for the scaled SEIR model (23), using vaccination as the time-dependent control function, ηt. Our main objective for this strategy is to reduce the number of infected and exposed individuals as well as the cost of vaccination. Therefore, in this control strategy, we minimize the objective functional Jη given by(31)Jη=∫0tfK1et+K2it+K32η2tdt,which subjects to(32)dstdt=−νstit−ηtst,s0=s0≥0,detdt=νstit−ρet,e0=e0≥0,ditdt=ρet−δit,i0=i0≥0,drtdt=δit+ηtst,r0=r0≥0,where the control set is given as follows:(33)G4=η:ηtis lebesgue measurable,0≤ηt≤1,t∈0,tf,and the parameters K1 and K2 are positive constants associated with exposed individuals and infected individuals, respectively, and K3 is a positive weight constant for control function ηt.The HamiltonianH is given by(34)H=K1et+K2it+K32η2t+φ1−νstit−ηtst+φ2νstit−ρet+φ3ρet−δit+φ4δi+ηtst.We then apply the Maximum Principle proposed by authors in [34] to determine an optimal solution.Assume thatw,η is an optimal solution for a controlled dynamical system; then, there exist adjoint vector function φ=φ1,φ2,…,φn, which satisfy the following equations:(35)dwdt=∂Ht,w,η,φ∂φ,0=∂Ht,w,η,φ∂η,dφdt=−∂Ht,w,η,φ∂w.Using the formulated Hamiltonian function (34) and equation (35), the adjoint equations and the control characterisation are given in the following theorem.Theorem 4.
Letη∗ be an optimal control and optimal state solutions s∗,e∗,i∗,andr∗ of the corresponding controlled dynamical system (31)-(32) that minimize Jη over G4. Then, there exist adjoint variables φi for i=1,2,3,4, satisfying(36)dφ1dt=νi∗tφ1−φ2+η∗tφ1−φ4,dφ2dt=−K1+ρφ2−φ3,dφ3dt=−K2+νs∗tφ1−φ2+δφ3−φ4,dφ4dt=0,with transversality conditions φ1tf=0,φ2tf=0,φ3tf=0,andφ4tf=0 and the control function η∗ is given by(37)η∗t=minmax0,s∗tK3φ1−φ4,1.


### 3.3. SEIR Model with Optimal Control Strategy 3

In this section, we present a controlled dynamical problem by incorporating treatment control,η1t, and educational campaign control, η2t, into the SEIR model (23). As we presented in Section 2.2, our main aim for this controlled problem is to minimize the number of infected individuals and the cost of educational campaigns and treatment. Therefore, in this control strategy, we minimize the objective functional given by(38)Jη1,η2=∫0tfD1it+D22η12t+D32η22tdt,which subjects to(39)dstdt=−νstit−η2tst,s0=s0≥0,detdt=νstit−ρet,e0=e0≥0,ditdt=ρet−δit−η1tit,i0=i0≥0,drtdt=δit+η1tit+η2tst,r0=r0≥0,where the control set is given as follows:(40)G5=η1,η2:ηitis lebesgue measurable,0≤ηit≤1,t∈0,tffori=1,2,where D1 is the weight constant on infected individuals and D2 and D3 are also positive weight parameters associated with controls η1t and η2t, respectively.The Hamiltonian functionH is given as(41)H=D1it+D22η12t+D32η22t+φ1−νstit−η2tst+φ2νstit−ρet+φ3ρet−δit−η1tit+φ4δit+η1tit+η2tst.We then apply the Maximum Principle proposed by authors in [34] to determine an optimal solution.Assume thatw,η is an optimal solution for a controlled dynamical system; then, there exist adjoint vector function φ=φ1,φ2,…,φn, which satisfy the following equations:(42)dwdt=∂Ht,w,η,φ∂φ,0=∂Ht,w,η,φ∂η,dφdt=−∂Ht,w,η,φ∂w.Using the formulated Hamiltonian function (41) and equation (42), the adjoint equations and the control characterisation are presented in the following theorem.Theorem 5.
Letη1∗ and η2∗ be optimal control pair and optimal state solutions s∗,e∗,i∗,andr∗ of the corresponding controlled dynamical system (38)-(39) that minimize Jη1,η2 over G5. Then, there exist adjoint variables φi for i=1,2,3,4, satisfying(43)dφ1dt=νi∗tφ1−φ2+η2∗tφ1−φ4,dφ2dt=ρφ2−φ3,dφ3dt=−D1+νs∗tφ1−φ2+δφ3−φ4+η1∗tφ3−φ4,dφ4dt=0,with transversality conditions φ1tf=0,φ2tf=0,φ3tf=0,andφ4tf=0, and the control pair η1∗ and η2∗ satisfies the optimality condition:(44)η1∗t=minmax0,i∗tD2φ3−φ4,1,η2∗t=minmax0,s∗tD3φ1−φ4,1.


### 3.4. Numerical Simulations and Discussion

As we did in Section2.3, we have used the forward-backward sweep numerical scheme with the fourth-order Runge–Kutta method to solve the optimality systems formulated for the three optimal control strategies. For our numerical simulations, we have adapted the same model parameter values (ν=0.2,δ=0.1,andρ=0.1887) as in the work by the authors in [33] with initial conditions: s0=0.88,e0=0.07,i0=0.05,andr0=0. We have assumed D=5,K1=1,K2=5,K3=5,D1=1,D2=5, and D3=5. Figure 6 shows solution paths for susceptible individuals with control strategies 1,2,and3 and without control. It is clear from the plot that there is a significant decrease in the number of susceptible individuals with control strategies than without control. In Figure 7, there is a rapid decrease in the number of exposed individuals with control strategies 1,2,and3 than without control. A similar effect can be observed in Figure 8 for the infected individuals. In Figure 9, there is a rapid increase in recovered individuals with control strategies 1,2,and3 than without control. Figures 10–12 represent optimal control functions for strategies 1, 2, and 3, respectively.Figure 6
Solution path for Susceptible individuals with three control strategies and without control,ν=0.2, ρ=0.1887, and δ=0.1.Figure 7
Solution paths for Exposed individuals with three control strategies and without control,ν=0.2, ρ=0.1887, and δ=0.1.Figure 8
Solution paths for Infected individuals with three control strategies and without control,ν=0.2, ρ=0.1887, and δ=0.1.Figure 9
Solution paths for Recovered individuals with three control strategies and without control,ν=0.2, ρ=0.1887, and δ=0.1.Figure 10
Optimal control functionηt for strategy 1.Figure 11
Optimal control functionsηt for strategy 2.Figure 12
Optimal control functionsη1t and η2t for strategy 3.


## 3.1. SEIR Model with Optimal Control Strategy 1

In this section, we consider a controlled dynamical system for the scaled SEIR model (23), using vaccination as time-dependent control function, ηt. As we did in Section 2.1, our main objective for this optimal control problem is that we want to reduce the number of infected individuals and the cost of vaccination. Therefore, in this control strategy, we minimize the objective functional Jη given by(24)Jη=∫0tfit+D2η2tdt,which subjects to(25)dstdt=−νstit−ηtst,s0=s0≥0,detdt=νstit−ρet,e0=e0≥0,ditdt=ρet−δit,i0=i0≥0,drtdt=δit+ηtst,r0=r0≥0,where the control set is given by(26)G3=η:ηtis lebesgue measurable,0≤ηt≤1,t∈0,tf,and the positive constant D is the weight on the cost of vaccination.The Hamiltonian functionH is then given as follows:(27)H=it+D2η2t+φ1−νstit−ηtst+φ2νstit−ρet+φ3ρet−δit+φ4δit+ηtst.We then apply the Maximum Principle proposed by authors in [34] to determine an optimal solution.Given thatw,η is an optimal solution for a controlled dynamical system, there exist adjoint vector function φ=φ1,φ2,…,φn, which satisfy the following equations:(28)dwdt=∂Ht,w,η,φ∂φ,0=∂Ht,w,η,φ∂η,dφdt=−∂Ht,w,η,φ∂w.Using the formulated Hamiltonian function (27) and equation (28), the adjoint equations and the control characterisation are presented in the following theorem.Theorem 3.
Letη∗ be an optimal control and optimal state solutions s∗,e∗,i∗,andr∗ of the corresponding controlled dynamical system (24)-(25) that minimize Jη over G3. Then, there exist adjoint variables φi for i=1,2,3,4, satisfying(29)dφ1dt=νi∗tφ1−φ2+η∗tφ1−φ4,dφ2dt=ρφ2−φ3,dφ3dt=−1+νs∗tφ1−φ2+δφ3−φ4,dφ4dt=0,with transversality conditions φ1tf=0,φ2tf=0,φ3tf=0,andφ4tf=0 and the control η∗t satisfies the optimality condition:(30)η∗t=minmax0,s∗tDφ1−φ4,1.


## 3.2. SEIR Model with Optimal Control Strategy 2

In this section, we present a controlled dynamical system for the scaled SEIR model (23), using vaccination as the time-dependent control function, ηt. Our main objective for this strategy is to reduce the number of infected and exposed individuals as well as the cost of vaccination. Therefore, in this control strategy, we minimize the objective functional Jη given by(31)Jη=∫0tfK1et+K2it+K32η2tdt,which subjects to(32)dstdt=−νstit−ηtst,s0=s0≥0,detdt=νstit−ρet,e0=e0≥0,ditdt=ρet−δit,i0=i0≥0,drtdt=δit+ηtst,r0=r0≥0,where the control set is given as follows:(33)G4=η:ηtis lebesgue measurable,0≤ηt≤1,t∈0,tf,and the parameters K1 and K2 are positive constants associated with exposed individuals and infected individuals, respectively, and K3 is a positive weight constant for control function ηt.The HamiltonianH is given by(34)H=K1et+K2it+K32η2t+φ1−νstit−ηtst+φ2νstit−ρet+φ3ρet−δit+φ4δi+ηtst.We then apply the Maximum Principle proposed by authors in [34] to determine an optimal solution.Assume thatw,η is an optimal solution for a controlled dynamical system; then, there exist adjoint vector function φ=φ1,φ2,…,φn, which satisfy the following equations:(35)dwdt=∂Ht,w,η,φ∂φ,0=∂Ht,w,η,φ∂η,dφdt=−∂Ht,w,η,φ∂w.Using the formulated Hamiltonian function (34) and equation (35), the adjoint equations and the control characterisation are given in the following theorem.Theorem 4.
Letη∗ be an optimal control and optimal state solutions s∗,e∗,i∗,andr∗ of the corresponding controlled dynamical system (31)-(32) that minimize Jη over G4. Then, there exist adjoint variables φi for i=1,2,3,4, satisfying(36)dφ1dt=νi∗tφ1−φ2+η∗tφ1−φ4,dφ2dt=−K1+ρφ2−φ3,dφ3dt=−K2+νs∗tφ1−φ2+δφ3−φ4,dφ4dt=0,with transversality conditions φ1tf=0,φ2tf=0,φ3tf=0,andφ4tf=0 and the control function η∗ is given by(37)η∗t=minmax0,s∗tK3φ1−φ4,1.


## 3.3. SEIR Model with Optimal Control Strategy 3

In this section, we present a controlled dynamical problem by incorporating treatment control,η1t, and educational campaign control, η2t, into the SEIR model (23). As we presented in Section 2.2, our main aim for this controlled problem is to minimize the number of infected individuals and the cost of educational campaigns and treatment. Therefore, in this control strategy, we minimize the objective functional given by(38)Jη1,η2=∫0tfD1it+D22η12t+D32η22tdt,which subjects to(39)dstdt=−νstit−η2tst,s0=s0≥0,detdt=νstit−ρet,e0=e0≥0,ditdt=ρet−δit−η1tit,i0=i0≥0,drtdt=δit+η1tit+η2tst,r0=r0≥0,where the control set is given as follows:(40)G5=η1,η2:ηitis lebesgue measurable,0≤ηit≤1,t∈0,tffori=1,2,where D1 is the weight constant on infected individuals and D2 and D3 are also positive weight parameters associated with controls η1t and η2t, respectively.The Hamiltonian functionH is given as(41)H=D1it+D22η12t+D32η22t+φ1−νstit−η2tst+φ2νstit−ρet+φ3ρet−δit−η1tit+φ4δit+η1tit+η2tst.We then apply the Maximum Principle proposed by authors in [34] to determine an optimal solution.Assume thatw,η is an optimal solution for a controlled dynamical system; then, there exist adjoint vector function φ=φ1,φ2,…,φn, which satisfy the following equations:(42)dwdt=∂Ht,w,η,φ∂φ,0=∂Ht,w,η,φ∂η,dφdt=−∂Ht,w,η,φ∂w.Using the formulated Hamiltonian function (41) and equation (42), the adjoint equations and the control characterisation are presented in the following theorem.Theorem 5.
Letη1∗ and η2∗ be optimal control pair and optimal state solutions s∗,e∗,i∗,andr∗ of the corresponding controlled dynamical system (38)-(39) that minimize Jη1,η2 over G5. Then, there exist adjoint variables φi for i=1,2,3,4, satisfying(43)dφ1dt=νi∗tφ1−φ2+η2∗tφ1−φ4,dφ2dt=ρφ2−φ3,dφ3dt=−D1+νs∗tφ1−φ2+δφ3−φ4+η1∗tφ3−φ4,dφ4dt=0,with transversality conditions φ1tf=0,φ2tf=0,φ3tf=0,andφ4tf=0, and the control pair η1∗ and η2∗ satisfies the optimality condition:(44)η1∗t=minmax0,i∗tD2φ3−φ4,1,η2∗t=minmax0,s∗tD3φ1−φ4,1.


## 3.4. Numerical Simulations and Discussion

As we did in Section2.3, we have used the forward-backward sweep numerical scheme with the fourth-order Runge–Kutta method to solve the optimality systems formulated for the three optimal control strategies. For our numerical simulations, we have adapted the same model parameter values (ν=0.2,δ=0.1,andρ=0.1887) as in the work by the authors in [33] with initial conditions: s0=0.88,e0=0.07,i0=0.05,andr0=0. We have assumed D=5,K1=1,K2=5,K3=5,D1=1,D2=5, and D3=5. Figure 6 shows solution paths for susceptible individuals with control strategies 1,2,and3 and without control. It is clear from the plot that there is a significant decrease in the number of susceptible individuals with control strategies than without control. In Figure 7, there is a rapid decrease in the number of exposed individuals with control strategies 1,2,and3 than without control. A similar effect can be observed in Figure 8 for the infected individuals. In Figure 9, there is a rapid increase in recovered individuals with control strategies 1,2,and3 than without control. Figures 10–12 represent optimal control functions for strategies 1, 2, and 3, respectively.Figure 6
Solution path for Susceptible individuals with three control strategies and without control,ν=0.2, ρ=0.1887, and δ=0.1.Figure 7
Solution paths for Exposed individuals with three control strategies and without control,ν=0.2, ρ=0.1887, and δ=0.1.Figure 8
Solution paths for Infected individuals with three control strategies and without control,ν=0.2, ρ=0.1887, and δ=0.1.Figure 9
Solution paths for Recovered individuals with three control strategies and without control,ν=0.2, ρ=0.1887, and δ=0.1.Figure 10
Optimal control functionηt for strategy 1.Figure 11
Optimal control functionsηt for strategy 2.Figure 12
Optimal control functionsη1t and η2t for strategy 3.


## 4. Conclusion

In this paper, we have studied epidemiological models for Ebola virus disease using nonlinear system of ordinary differential equation and optimal control theory. We have used indirect methods in optimal control to study existing mathematical models proposed by the authors in [32, 33]. Using Pontryagin’s maximum principle and Hamiltonian function, we have derived adjoint equations and optimality system for the various optimal control strategies. From the simulation results, we observed that the SIR model with optimal control strategies shows significant decrease in the proportions of infected and susceptible individuals and a rapid increase in the recovered individuals compared to the SIR model without control. A similar effect was observed in the SEIR model with control strategies where there was a significant decrease in Susceptible, Exposed, and Infected individuals and a rapid increase in the Recovered individuals. Our numerical results are similar to the once generated by the authors in [32, 33] who applied direct methods in optimal control for their mathematical models. Therefore, following the numerical results, we can conclude that effective educational campaigns and vaccination of susceptible individuals as well as effective treatments of infected individuals can help reduce the disease transmission. In this work, we have also analytically and numerically studied one additional optimal control strategy for the SIR epidemic model.



---

*Source: 2823816-2020-07-04.xml*