# Prevention of Leptospirosis Infected Vector and Human Population by Multiple Control Variables

Authors: Muhammad Altaf Khan; Saeed Islam; Sher Afzal Khan; Ilyas Khan; Sharidan Shafie; Taza Gul
Journal: Abstract and Applied Analysis
 (2014)

---

## Abstract

Leptospirosis is an infectious disease that damages the liver and kidneys, found mainly in dogs and farm animals and caused by bacteria. In this paper, we present the optimal control problem applied to a dynamical leptospirosis infected vector and human population by using multiple control variables. First, we show the existence of the control problem and then use analytical and numerical techniques to investigate the existence cost effective control efforts for prevention of indirect and direct transmission of this disease. In order to do this, we consider three control functions two for human
and one for vector population. We completely characterize the optimal control problem and compute the numerical solution of the optimality system using an iterative method.


---

## Body

## 1. Introduction

Leptospirosis disease is a globally important infectious disease. The disease is caused by a bacteria which is calledLeptospira. Human as well as cattle is mostly infected from this disease [1]. The human is infected by means of drinking the water in which a rat was found dead, and cattle that drink this water become infectious. The human whose urine is used by other animals and cattle is also infected because the leptospirosis disease germs come out in urine. Those who wade through dirty water are mostly infected from this disease. Weil’s first time described leptospirosis as a unique disease process in 1886, while 30 years before Inada and his colleagues identified the causal organism. The symptoms of leptospirosis are high fever, headache, chills, muscle aches, conjunctivitis (red eyes), diarrhea, vomiting, and kidney or liver problems (which may also include jaundice), anemia, and, sometimes, rash. Symptoms may last from a few days and up to several weeks. Deaths from this disease may occur, but they are rare. For some cases the infections can be mild and without obvious symptom [2–6].Many models have been proposed to represent the dynamics of both human and vector population [7–9]. Pongsuumpun et al. [10] developed mathematical models to study the behavior of leptospirosis disease. They represent the rate of change for both rats and human population. The human population is further divided into two main groups: juveniles and adults. Triampo et al. [11] considered a deterministic model for the transmission of leptospirosis disease [11]. In their work they considered a number of leptospirosis infections in Thailand and showed the numerical simulations. Zaman [12] considered the real data presented in [11] to study the dynamical behavior and role of optimal control theory. The dynamical interaction including local and global stability of leptospirosis infected vector and human population can be found in Zaman et al. [13]. In their work they also presented the bifurcation analysis and presented the numerical simulations for different values of infection rate. In [14], the author presented an epidemic model of malaria, by using three control variables, and obtained their optimal solutions; for more references, see, for example, [14–17].In this paper, we consider the basic model studied in [13] to incorporate some important epidemiological features. We use optimal control theory to reduce the proportion of the infected human and infected vector population by using multiple control. At the long-term level of infected human, every infected human on average causes one further secondary case. Therefore, if we can reduce the number of infected humans further, the disease does less well and increase the recovered human. Here we define the control variables; the first control is to cover all cuts and wear dry, full-cover boots, shoes and long sleeve shirts when handling animals. The second control represents washing hands thoroughly on a regular basis and showering after work. Cleaning up both work place and home is our third control. To do this, we first show the existence of the optimal control system. Then, by solving the optimality system analytically, which consists of the original state system, the adjoint system and their boundary conditions. The real data presented for leptospirosis epidemic in Thailand have been used in numerical simulations. We also conclude by discussing the results of the numerical simulations for our epidemic mathematical model.The structure of the paper is organized as follows. Section2 is devoted to the mathematical formulation and solution of the problem which is further divided into the following subsections. Section 2.1 includes the formulation of the mathematical model followed by Section 2.2 which contains optimal control problem whereas Section 2.3 explains the bifurcation analysis of the control problem. The last Section 2.4 discusses the existence of control problem. Section 3 is devoted to the numerical solution of the optimality system. Finally, the conclusion is presented in Section 4.


## 2. Mathematical Formulation and Solution

### 2.1.  Mathematical Model

In this section, a vector-host epidemic model with direct transmission is presented. The host population at timet is divided into susceptible S
h
(
t
), I
h
(
t
) infected and recovered R
h
(
t
) individuals. The vector (rats) populations at time t is divided into susceptible S
v
(
t
) and infected vector population I
v
(
t
). The total population of human is denoted by N
h and the total population of the vector is denoted by N
v. Thus, N
h
(
t
)
=
S
h
(
t
)
+
I
h
(
t
)
+
R
h
(
t
) and N
v
(
t
)
=
S
v
(
t
)
+
I
v
(
t
). The mathematical representation of the model which consists of the system of nonlinear differential equations with five state variables is given by
(1)
d
S
h
d
t
=
b
1
-
μ
h
S
h
-
β
2
S
h
I
v
-
β
1
S
h
I
h
+
λ
h
R
h
,
d
I
h
d
t
=
β
2
S
h
I
v
+
β
1
S
h
I
h
-
μ
h
I
h
-
δ
h
I
h
-
γ
h
I
h
,
d
R
h
d
t
=
γ
h
I
h
-
μ
h
R
h
-
λ
h
R
h
,
d
S
v
d
t
=
b
2
-
γ
v
S
v
-
β
3
S
v
I
h
,
d
I
v
d
t
=
β
3
S
v
I
h
-
γ
v
I
v
-
δ
v
I
v
.
Here b
1 is the recruitment rate of human population, susceptible human can be infected by two ways of transmission, that is, directly or through infected individuals, and β
1, β
2 are the mediate transmission coefficients. μ
h is the natural mortality rate for human; γ
h is the recovery rate for human from the infections. We assumed that disease may be fatal to some infectious host, so disease related death rate from infected class occurs at human populations at δ
h. The immune human once again is susceptible at constant rate λ
h. b
2 is the recruitment rate for vector population. The infectious vector dies due to disease at vector populations at the rate of δ
v. β
3 represents the disease carrying of susceptible vector per host per unit time; γ
v is the death rate of vector.


### 2.2. Optimal Control Problem

Optimal control theory is a powerful mathematical tool which makes the decision involving complex dynamical systems [18]. Optimal control method has been used to study the dynamics of the disease; we refer the reader to [19–21]; no such method is used, according to the author’s knowledge, to determine optimal control measure for vector-host epidemic direct transmission. The problem is to minimize the infected human and vector population and to maximize the recovered human population. In the system (1) we have five state variables S
h
(
t
), I
h
(
t
), R
h
(
t
), S
v
(
t
), and I
v
(
t
). In this optimal control problem we use three control variables.(i)
The first controlu
1
(
t
) (Figure 7) represents that human should cover all cuts, abrasions with waterproof dressing, grazes, wear dry clothes, wear full-cover shoes, gloves, and use the shirts with long sleeves when handling the animals.
(ii)
The second controlu
2
(
t
) (Figure 8) shows that after work human should bath or shower regularly and adopt the habit of washing hands regularly.
(iii)
Our third controlu
3
(
t
) (Figure 9) represents cleaning the home and working area.In the human population, the associated force of infections are reduced by factors of(
1
-
u
1
(
t
)
) and (
1
-
u
2
(
t
)
), respectively. We assume that the mortality rate of vector population increases at a rate proportional to u
3
(
t
), where ϵ
1
>
0 and ϵ
2
>
0 are rate constants. Taking into account the extensions and assumptions made above, it follows that the dynamics of the system (1) are governed by the following system of five differential equations:
(2)
d
S
h
d
t
=
b
1
-
μ
h
S
h
-
β
2
S
h
I
v
1
-
u
1
t
d
S
h
d
t
=
-
β
1
S
h
I
h
1
-
u
2
t
+
λ
h
R
h
,
d
I
h
d
t
=
β
2
S
h
I
v
1
-
u
1
t
+
β
1
S
h
I
h
1
-
u
2
t
d
I
h
d
t
=
-
μ
h
I
h
-
δ
h
I
h
-
γ
h
I
h
-
u
2
t
I
h
,
d
R
h
d
t
=
γ
h
I
h
-
μ
h
R
h
-
λ
h
R
h
+
u
1
t
S
h
+
u
2
t
I
h
,
d
S
v
d
t
=
b
2
-
γ
v
S
v
-
β
3
S
v
I
h
-
ϵ
1
u
3
t
S
v
,
d
I
v
d
t
=
β
3
S
v
I
h
-
γ
v
I
v
-
δ
v
I
v
-
ϵ
2
u
3
(
t
)
I
v
with initials conditions
(3)
S
h
(
0
)
≥
0
,
I
h
(
0
)
≥
0
,
R
h
(
0
)
≥
0
,
S
v
(
0
)
≥
0
,
I
v
(
0
)
≥
0
,
where A
1, A
2, B
1, and B
2 are positive constants to keep the balance of the size of the individuals S
h
(
t
), I
h
(
t
), S
v
(
t
), and I
v
(
t
), respectively. For the optimal control problem we consider the control variable u
1
∈
U. Here, u
1, u
2, u
3
∈
U is measurable; 0
≤
u
1
(
t
), u
2
(
t
), u
3
(
t
)
≤
1.


### 2.3. Bifurcation Analysis of the Control Problem

In this subsection, we present the endemic equilibria which is further used in the bifurcation analysis. In order to do this we set the left-hand side of the system (2) equal to zero and use the technique developed in [22], getting
(4)
S
h
*
=
P
1
P
2
(
γ
v
+
β
3
I
h
*
+
ϵ
1
u
3
)
β
2
β
3
b
2
(
1
-
u
1
)
+
P
1
β
1
(
1
-
u
2
)
(
γ
v
+
β
3
I
h
*
+
ϵ
1
u
3
)
,
S
v
*
=
b
2
γ
v
+
β
3
I
h
*
+
ϵ
1
u
3
,
I
v
*
=
β
3
b
2
I
h
*
P
1
(
γ
v
+
β
3
I
h
*
+
ϵ
1
u
3
)
,
R
h
*
=
P
3
I
h
*
β
2
β
3
b
2
1
-
u
1
×
γ
v
+
β
3
I
h
*
+
ϵ
1
u
3
s
s
s
s
s
s
s
s
s
s
s
s
+
P
1
β
1
1
-
u
2
γ
v
+
β
3
I
h
*
+
ϵ
1
u
3
+
u
1
s
s
s
s
s
s
s
×
P
1
P
2
γ
v
+
β
3
I
h
*
+
ϵ
1
u
3
s
s
s
s
s
s
s
×
P
4
β
2
β
3
b
2
1
-
u
1
+
P
1
β
1
1
-
u
2
s
s
s
s
s
s
s
s
s
s
s
s
s
s
×
γ
v
+
β
3
I
h
*
+
ϵ
1
u
3
-
1
,
where
(5)
P
1
=
γ
v
+
δ
v
+
ϵ
2
u
3
,
P
2
=
μ
h
+
δ
h
+
γ
h
+
u
2
,
P
3
=
μ
h
+
λ
h
.
The reproduction number for the control system (2) is given by
(6)
K
u
=
b
1
μ
h
β
2
β
3
b
2
1
-
u
1
P
1
P
2
γ
v
-
ϵ
1
u
3
+
β
1
1
-
u
2
P
2
.
The reproduction number for without control system is given by
(7)
K
0
=
b
1
μ
h
β
2
β
3
b
2
γ
v
P
0
1
P
0
2
+
β
1
P
0
2
,
where P
0
1
=
(
δ
v
+
γ
v
) and P
0
2
=
(
μ
h
+
δ
h
+
γ
h
). In the above expression for the endemic equilibria the infected component I
h
* is nonzero. Using the value of S
h
*, R
h
*, and I
v
* in the first equation of the system (2), we obtained
(8)
f
I
h
*
=
a
I
h
*
2
+
b
I
h
*
+
c
=
0
,
where
(9)
a
=
λ
h
β
3
2
P
1
β
1
(
1
-
u
2
)
,
b
=
λ
h
β
3
2
β
2
b
2
1
-
u
1
+
P
4
b
1
β
1
P
1
1
-
u
2
β
3
b
=
+
λ
h
β
3
P
1
γ
v
+
ϵ
1
u
3
β
1
1
-
u
2
b
=
-
P
2
P
4
β
2
1
-
u
1
β
3
β
2
b
2
-
β
1
1
-
u
2
P
1
β
3
b
=
+
u
1
λ
h
P
1
β
3
,
c
=
-
P
1
P
2
P
4
γ
v
+
ϵ
1
u
3
μ
h
β
1
1
-
u
2
-
u
1
λ
h
1
-
K
u
.
Here the coefficient a is positive always and c depends upon the value of K
u; if the value of K
u
<
1, then c is positive; otherwise negative. The positive solution of the above equation depends upon the value of b and c. For the value of K
u
>
1, the above equation leads to two roots, positive and negative. If we substitute K
u
=
1, then the equation has no positive solution. This is possible if and only if b
<
0. For b
<
0 and K
u
=
1, the equilibria depend upon K
u; then there exists an open interval having two positive roots, that is, I
1
=
-
b
-
b
2
-
4
a
c
/
2
a and I
2
=
-
b
+
b
2
-
4
a
c
/
2
a. For c
>
0 either b
2
<
4
a
c or b
≥
0, then the above have no positive solution.For backward bifurcation, we setb
2
-
4
a
c
=
0 and K
u
=
R
c and solve for the critical value of K
c, which is given by
(10)
R
c
=
1
-
b
2
4
a
P
1
P
2
P
4
γ
v
+
ϵ
1
u
3
μ
h
β
1
1
-
u
2
-
u
1
λ
h
.
It is proved by simulating the set of parameter values presented in Table 1. Figure 1 shows the bifurcation of the control system (2). The high bold black line shows the bifurcation at u
1
=
u
2
=
u
3
=
0, that is, without control. The small black bold line shows the bifurcation at u
1
=
u
2
=
u
3
≠
0, that is, with control. The change occurs in the bifurcation rapidly when we apply the control variable. u
1 is the first control, which is represented by the bold dotted line, the dotted dashed line represents the second control, which shows the change in the control, the very narrow dotted line represent the third control, and bifurcation occurs for the control variable.Table 1
Parameter values used in the control bifurcation.
Notation
Value
b
2
10
β
1
0.048
β
3
0.006
u
2
0.001
λ
h
0.068
γ
h
0.041
b
1
10
β
2
0.0008
u
1
0.541
u
3
0.081
μ
h
0.091
δ
h
0.10
δ
v
0.009
ϵ
2
0.009
ϵ
1
0.00001
γ
v
0.087Figure 1
The plot shows the backward bifurcation for control and without control system.


### 2.4. Existence of Control Problem

We use the bounded Lebesgue measurable control and define our objective functional as(11)
J
(
u
1
,
u
2
,
u
3
)
=
min
⁡
∫
0
T
+
1
2
W
1
u
1
2
+
W
2
u
2
2
+
W
3
u
3
2
A
1
S
h
(
t
)
+
A
2
I
h
(
t
)
+
B
1
S
v
(
t
)
+
B
2
I
v
(
t
)
sssssssssssssss
+
1
2
W
1
u
1
2
+
W
2
u
2
2
+
W
3
u
3
2
d
t
.
Here, in (11), A
1, A
2, B
1, B
2, W
1, W
2, and W
3 represent the weight/balance factors just to keep the balance of individuals in the objective functional. The control set is defined as
(12)
U
=
u
=
u
1
,
u
2
,
u
3
∣
u
i
t
is Lebesgue measureable
,
s
s
s
s
0
≤
u
i
t
≤
1
,
t
∈
0
,
T
,
for
i
=
1,2
,
3
.Note. The details of the control variables u
1, u
2, and u
3 are available in Section 2.2.First, we show the existence for the control system (2). Let S
h
(
t
), I
h
(
t
), R
h
(
t
), S
v
(
t
), and I
v
(
t
) be the state variables with control variables u
1
(
t
), u
2
(
t
), and u
3
(
t
). For existence we consider the control system (2). Then we can write the system (2) in the following form:
(13)
X
′
=
A
X
+
F
X
,
where(14)
X
′
=
S
h
t
I
h
t
R
h
t
S
v
t
I
v
t
,
A
=
-
μ
h
0
λ
h
0
0
0
-
μ
h
-
δ
h
-
γ
h
-
u
2
t
0
0
0
u
1
γ
h
+
u
2
t
-
μ
h
-
λ
h
0
0
0
0
0
-
γ
v
-
ϵ
1
u
3
t
0
0
0
0
0
-
γ
v
-
δ
v
-
ϵ
2
u
3
t
,
F
X
=
b
1
-
β
2
S
h
t
I
v
1
-
u
1
t
-
β
1
S
h
t
I
h
t
1
-
u
2
t
β
2
S
h
t
I
v
1
-
u
1
t
+
β
1
S
h
t
I
h
t
1
-
u
2
t
1
-
u
2
t
0
b
2
-
β
3
S
v
t
I
h
t
β
3
S
v
t
I
h
t
,where X
′ denotes the derivative with respect to time t. The system (13) is a nonlinear system with a bounded coefficient. We set
(15)
G
(
X
)
=
A
X
+
F
(
X
)
.
The second term on the right-hand side of (15) satisfies
(16)
F
(
X
1
)
-
F
(
X
2
)
≤
C
1
S
1
h
t
-
S
2
h
t
+
C
2
I
1
h
t
-
I
2
h
t
ssssissss
+
C
3
R
1
h
t
-
R
2
h
t
+
C
4
S
1
v
t
-
S
2
v
t
ssssssss
+
C
5
I
1
v
t
-
I
2
v
t
≤
C
S
1
h
t
-
S
2
h
t
+
I
1
h
t
-
I
2
h
t
sssssss
+
R
1
h
t
-
R
2
h
t
+
S
1
v
t
-
S
2
v
t
sssissss
+
I
1
v
t
-
I
2
v
t
,
where the positive constant C
=
max
⁡
(
C
1
,
C
2
,
C
3
,
C
4
,
C
5
) is independent of the state variables. Also we have
(17)
G
X
1
-
G
X
2
≤
C
X
1
-
X
2
,
where C
=
C
1
+
C
2
+
C
3
+
C
4
+
C
5
+
M
<
∞. So, it follows that the function G is uniformly Lipschitz continuous. From the definition of control variables and nonnegative initial conditions we can see that a solution of the system (13) exists; see [23]. Now, we consider the control system (2) with the initial conditions (3) to show the existence of the control problem. Note that for bounded Lebesgue measurable controls and nonnegative initial conditions, nonnegative bounded solutions to the state system exist [23]. Let us go back to the optimal control problem (2)-(3). In order to find an optimal solution, first we will find the Lagrangian L and Hamiltonian H for the optical control problem (2)-(3). The Lagrangian for our control problem is given by
(18)
L
=
A
1
S
h
+
A
2
I
h
+
B
1
S
v
+
B
2
I
v
+
1
2
W
1
u
1
2
+
W
2
u
2
2
+
W
3
u
3
2
.
For the minimum value of Lagrangian, we define the Hamiltonian for the control problem:
(19)
H
=
L
S
h
,
I
h
,
S
v
,
I
v
,
u
1
,
u
2
,
u
3
+
λ
1
d
S
h
d
t
+
λ
2
d
I
h
d
t
+
λ
3
d
R
h
d
t
+
λ
4
d
S
v
d
t
+
λ
5
d
I
v
d
t
.
For the existence of our control problem, we state and prove the following theorem.Theorem 1.
There exists an optimal controlu
*
=
(
u
1
*
,
u
2
*
,
u
3
*
)
∈
U such that
(20)
J
u
1
*
,
u
2
*
,
u
3
*
=
min
⁡
u
1
,
u
2
,
u
3
∈
U
J
u
1
,
u
2
,
u
3
,
subject to the control system (2) with the initial conditions (3).Proof.
To prove the existence of an optimal control, we use the result in [24]; the control and the state variable are nonnegative values. In this minimizing of the problem, the necessary convexity of the objective functional in u
1, u
2, and u
3 is satisfied. The set of control variables (
u
1
,
u
2
,
u
3
)
∈
U is also convex and closed by the definition. The optimal system is bounded which determines the compactness needed for the existence of optimal control. The integrand in the objective functional (11) is A
1
S
(
t
)
+
A
2
I
(
t
)
+
B
1
S
v
+
B
2
I
v
+
1
/
2
(
W
1
u
1
2
+
W
2
u
2
2
+
W
3
u
3
2
) which is convex in the control set U. Also we can easily see that there exist a constant σ
>
1 and positive numbers ω
1 and ω
2 such that
(21)
J
u
1
,
u
2
,
u
3
≥
ω
1
u
1
2
+
u
2
2
+
u
3
2
σ
/
2
-
ω
2
,
which is the existence of an optimal control problem.To find the optimal solution, we apply Pontryagin’s maximum principle [25] given by the following.If(
x
,
u
) is an optimal solution for an optimal control problem, then there exists a nontrivial vector function λ
=
(
λ
1
,
λ
2
,
…
,
λ
n
) which satisfies the following inequalities:
(22)
d
x
d
t
=
∂
H
t
,
x
,
u
,
λ
∂
λ
,
0
=
∂
H
t
,
x
,
u
,
λ
∂
u
,
d
λ
d
t
=
-
∂
H
t
,
x
,
u
,
λ
∂
x
.
Now we apply the necessary conditions to the Hamiltonian H in (19).Theorem 2.
Suppose thatS
h
*, I
h
*, R
h
*, S
v
* and I
v
* are the optimal state solutions with associated optimal control variables (
u
1
*
,
u
2
*
,
u
3
*
) for the optimal control problem (2)-(3). Then there exist adjoint variables λ
i, for i
=
1,2
,
…
,
5, satisfying
(23)
d
λ
1
d
t
=
λ
1
-
λ
2
β
2
I
v
1
-
u
1
+
β
1
I
h
1
-
u
2
d
λ
1
d
t
=
+
λ
1
μ
h
-
u
1
λ
3
-
A
1
,
d
λ
2
d
t
=
(
λ
1
-
λ
2
)
β
1
S
h
(
1
-
u
2
)
+
μ
h
λ
2
+
δ
h
λ
2
+
γ
h
λ
2
d
λ
2
d
t
=
+
u
2
λ
2
-
λ
3
γ
h
-
λ
3
u
2
+
β
3
S
v
λ
4
-
λ
5
-
A
2
,
d
λ
3
d
t
=
λ
3
-
λ
1
λ
h
+
μ
h
λ
3
,
d
λ
4
d
t
=
λ
4
-
λ
5
β
3
I
h
+
λ
4
γ
v
+
λ
4
ϵ
1
u
3
-
B
1
,
d
λ
5
d
t
=
λ
1
-
λ
2
β
2
S
h
1
-
u
1
d
λ
5
d
t
=
-
λ
5
μ
v
+
λ
5
γ
v
+
λ
5
δ
v
+
λ
5
ϵ
2
u
3
-
B
2
with transversality conditions (or boundary conditions)
(24)
λ
i
T
=
0
,
i
=
1,2
,
…
,
5
.
Furthermore, optimal controls u
1
*, u
2
*, and u
3
* are given by
(25)
u
1
*
=
max
⁡
min
⁡
λ
2
-
λ
1
β
2
S
h
*
I
v
*
-
λ
3
S
h
*
W
1
,
1
,
0
,
u
2
*
=
max
⁡
min
⁡
λ
2
-
λ
1
β
1
S
h
*
I
h
*
+
λ
2
-
λ
3
I
h
*
W
2
,
1
,
0
,
u
3
*
=
max
⁡
min
⁡
λ
4
ϵ
1
S
v
*
+
λ
5
ϵ
2
I
v
*
W
3
,
1
,
0
.Proof.
To find the adjoint equations and the transversality conditions, we use the Hamiltonian (19). By setting S
h
t
=
S
h
*
t, I
h
t
=
I
h
*
t, R
h
t
=
R
h
*
t, S
v
t
=
S
v
*
t, and I
v
t
=
I
v
*
t and differentiating the Hamiltonian (19) with respect to S
h, I
h, R
h, S
v, and I
v, respectively, we get (23). Then solving the equations ∂
H
/
∂
u
1
=
0, ∂
H
/
∂
u
2
=
0, and ∂
H
/
∂
u
3
=
0 on the interior of the control set and then using the optimality conditions and also the property of control space U, we can derive (25).
Here we call formulas (25) for the characterization of the optimal control. The optimal control and the state are determined by solving the optimality system, which consist of the state system (1), the adjoint system (23), initial conditions at (3), boundary conditions (24), and the characterization of the optimal control (
u
1
*
,
u
2
*
,
u
3
*
) which is given by (25). In addition, the second derivative of the Lagrangian with respect to u
1, u
2, and u
3, respectively, is positive, which shows the minimum of the optimal control u
1
*, u
2
*, and u
3
*. Substituting the values of u
1
*, u
2
*, and u
3
* in the control system (2), we obtain the following system:
(26)
d
S
h
*
d
t
=
b
1
-
μ
h
S
h
*
-
β
2
S
h
*
I
v
*
×
1
-
max
⁡
⁡
min
⁡
⁡
λ
2
-
λ
1
β
2
S
h
*
I
v
*
-
λ
3
S
h
*
W
1
,
1
,
0
-
β
1
S
h
*
I
h
*
1
-
max
⁡
min
⁡
λ
2
-
λ
1
β
1
S
h
*
I
h
*
+
λ
2
-
λ
3
I
h
*
W
2
,
1
,
0
+
λ
h
R
h
*
,
d
I
h
*
d
t
=
β
2
S
h
*
I
v
*
1
-
max
⁡
⁡
min
⁡
⁡
λ
2
-
λ
1
β
2
S
h
*
I
v
*
-
λ
3
S
h
*
W
1
,
1
,
0
-
μ
h
I
h
*
-
δ
h
I
h
*
+
β
1
S
h
*
I
h
*
×
1
-
max
⁡
⁡
min
⁡
λ
2
-
λ
1
β
1
S
h
*
I
h
*
+
λ
2
-
λ
3
I
h
*
W
2
,
1
,
0
-
γ
h
I
h
*
-
max
⁡
min
⁡
λ
2
-
λ
1
β
1
S
h
*
I
h
*
+
λ
2
-
λ
3
I
h
*
W
2
,
1
,
0
I
h
*
,
d
R
h
*
d
t
=
γ
h
I
h
*
-
μ
h
R
h
*
-
λ
h
R
h
*
+
max
⁡
⁡
min
⁡
λ
2
-
λ
1
β
2
S
h
*
I
v
*
-
λ
3
S
h
*
W
1
,
1
,
0
S
h
*
+
max
⁡
⁡
min
⁡
⁡
λ
2
-
λ
1
β
1
S
h
*
I
h
*
+
λ
2
-
λ
3
I
h
*
W
2
,
1
,
0
I
h
*
,
d
S
v
*
d
t
=
b
2
-
γ
v
S
v
*
-
β
3
S
v
*
I
h
*
-
ϵ
1
S
v
*
max
⁡
⁡
min
⁡
⁡
λ
4
ϵ
1
S
v
*
+
λ
5
ϵ
2
I
v
*
W
3
,
1
,
0
,
d
I
v
*
d
t
=
β
3
S
v
*
I
h
*
+
μ
v
I
v
*
-
γ
v
I
v
*
-
δ
v
I
v
*
-
ϵ
2
I
v
*
max
⁡
⁡
min
⁡
⁡
λ
4
ϵ
1
S
v
*
+
λ
5
ϵ
2
I
v
*
W
3
,
1
,
0
with Hamiltonian H
* at (
t
,
S
h
*
,
I
h
*
,
R
h
*
,
S
v
*
,
I
v
*
,
u
1
*
,
u
2
*
,
u
3
*
,
λ
1
,
λ
2
,
λ
3
,
λ
4
,
λ
5
):
(27)
H
*
=
A
1
S
h
*
+
A
2
I
h
*
+
B
1
S
v
*
+
B
2
I
v
*
+
1
2
W
1
max
⁡
⁡
min
⁡
⁡
λ
2
-
λ
1
β
2
S
h
*
I
v
*
-
λ
3
S
h
*
W
1
,
1
,
0
2
+
W
2
×
max
⁡
⁡
min
⁡
⁡
λ
2
-
λ
1
β
1
S
h
*
I
h
*
+
λ
2
-
λ
3
I
h
*
W
2
,
1
,
0
2
+
W
3
×
max
⁡
min
⁡
λ
4
ϵ
1
S
v
*
+
λ
5
ϵ
2
I
v
*
W
3
,
1
,
0
2
+
λ
1
d
S
h
*
d
t
+
λ
2
d
I
h
*
d
t
+
λ
3
d
R
h
*
d
t
+
λ
4
d
S
v
*
d
t
+
λ
5
d
I
v
*
d
t
.


## 2.1.  Mathematical Model

In this section, a vector-host epidemic model with direct transmission is presented. The host population at timet is divided into susceptible S
h
(
t
), I
h
(
t
) infected and recovered R
h
(
t
) individuals. The vector (rats) populations at time t is divided into susceptible S
v
(
t
) and infected vector population I
v
(
t
). The total population of human is denoted by N
h and the total population of the vector is denoted by N
v. Thus, N
h
(
t
)
=
S
h
(
t
)
+
I
h
(
t
)
+
R
h
(
t
) and N
v
(
t
)
=
S
v
(
t
)
+
I
v
(
t
). The mathematical representation of the model which consists of the system of nonlinear differential equations with five state variables is given by
(1)
d
S
h
d
t
=
b
1
-
μ
h
S
h
-
β
2
S
h
I
v
-
β
1
S
h
I
h
+
λ
h
R
h
,
d
I
h
d
t
=
β
2
S
h
I
v
+
β
1
S
h
I
h
-
μ
h
I
h
-
δ
h
I
h
-
γ
h
I
h
,
d
R
h
d
t
=
γ
h
I
h
-
μ
h
R
h
-
λ
h
R
h
,
d
S
v
d
t
=
b
2
-
γ
v
S
v
-
β
3
S
v
I
h
,
d
I
v
d
t
=
β
3
S
v
I
h
-
γ
v
I
v
-
δ
v
I
v
.
Here b
1 is the recruitment rate of human population, susceptible human can be infected by two ways of transmission, that is, directly or through infected individuals, and β
1, β
2 are the mediate transmission coefficients. μ
h is the natural mortality rate for human; γ
h is the recovery rate for human from the infections. We assumed that disease may be fatal to some infectious host, so disease related death rate from infected class occurs at human populations at δ
h. The immune human once again is susceptible at constant rate λ
h. b
2 is the recruitment rate for vector population. The infectious vector dies due to disease at vector populations at the rate of δ
v. β
3 represents the disease carrying of susceptible vector per host per unit time; γ
v is the death rate of vector.


## 2.2. Optimal Control Problem

Optimal control theory is a powerful mathematical tool which makes the decision involving complex dynamical systems [18]. Optimal control method has been used to study the dynamics of the disease; we refer the reader to [19–21]; no such method is used, according to the author’s knowledge, to determine optimal control measure for vector-host epidemic direct transmission. The problem is to minimize the infected human and vector population and to maximize the recovered human population. In the system (1) we have five state variables S
h
(
t
), I
h
(
t
), R
h
(
t
), S
v
(
t
), and I
v
(
t
). In this optimal control problem we use three control variables.(i)
The first controlu
1
(
t
) (Figure 7) represents that human should cover all cuts, abrasions with waterproof dressing, grazes, wear dry clothes, wear full-cover shoes, gloves, and use the shirts with long sleeves when handling the animals.
(ii)
The second controlu
2
(
t
) (Figure 8) shows that after work human should bath or shower regularly and adopt the habit of washing hands regularly.
(iii)
Our third controlu
3
(
t
) (Figure 9) represents cleaning the home and working area.In the human population, the associated force of infections are reduced by factors of(
1
-
u
1
(
t
)
) and (
1
-
u
2
(
t
)
), respectively. We assume that the mortality rate of vector population increases at a rate proportional to u
3
(
t
), where ϵ
1
>
0 and ϵ
2
>
0 are rate constants. Taking into account the extensions and assumptions made above, it follows that the dynamics of the system (1) are governed by the following system of five differential equations:
(2)
d
S
h
d
t
=
b
1
-
μ
h
S
h
-
β
2
S
h
I
v
1
-
u
1
t
d
S
h
d
t
=
-
β
1
S
h
I
h
1
-
u
2
t
+
λ
h
R
h
,
d
I
h
d
t
=
β
2
S
h
I
v
1
-
u
1
t
+
β
1
S
h
I
h
1
-
u
2
t
d
I
h
d
t
=
-
μ
h
I
h
-
δ
h
I
h
-
γ
h
I
h
-
u
2
t
I
h
,
d
R
h
d
t
=
γ
h
I
h
-
μ
h
R
h
-
λ
h
R
h
+
u
1
t
S
h
+
u
2
t
I
h
,
d
S
v
d
t
=
b
2
-
γ
v
S
v
-
β
3
S
v
I
h
-
ϵ
1
u
3
t
S
v
,
d
I
v
d
t
=
β
3
S
v
I
h
-
γ
v
I
v
-
δ
v
I
v
-
ϵ
2
u
3
(
t
)
I
v
with initials conditions
(3)
S
h
(
0
)
≥
0
,
I
h
(
0
)
≥
0
,
R
h
(
0
)
≥
0
,
S
v
(
0
)
≥
0
,
I
v
(
0
)
≥
0
,
where A
1, A
2, B
1, and B
2 are positive constants to keep the balance of the size of the individuals S
h
(
t
), I
h
(
t
), S
v
(
t
), and I
v
(
t
), respectively. For the optimal control problem we consider the control variable u
1
∈
U. Here, u
1, u
2, u
3
∈
U is measurable; 0
≤
u
1
(
t
), u
2
(
t
), u
3
(
t
)
≤
1.


## 2.3. Bifurcation Analysis of the Control Problem

In this subsection, we present the endemic equilibria which is further used in the bifurcation analysis. In order to do this we set the left-hand side of the system (2) equal to zero and use the technique developed in [22], getting
(4)
S
h
*
=
P
1
P
2
(
γ
v
+
β
3
I
h
*
+
ϵ
1
u
3
)
β
2
β
3
b
2
(
1
-
u
1
)
+
P
1
β
1
(
1
-
u
2
)
(
γ
v
+
β
3
I
h
*
+
ϵ
1
u
3
)
,
S
v
*
=
b
2
γ
v
+
β
3
I
h
*
+
ϵ
1
u
3
,
I
v
*
=
β
3
b
2
I
h
*
P
1
(
γ
v
+
β
3
I
h
*
+
ϵ
1
u
3
)
,
R
h
*
=
P
3
I
h
*
β
2
β
3
b
2
1
-
u
1
×
γ
v
+
β
3
I
h
*
+
ϵ
1
u
3
s
s
s
s
s
s
s
s
s
s
s
s
+
P
1
β
1
1
-
u
2
γ
v
+
β
3
I
h
*
+
ϵ
1
u
3
+
u
1
s
s
s
s
s
s
s
×
P
1
P
2
γ
v
+
β
3
I
h
*
+
ϵ
1
u
3
s
s
s
s
s
s
s
×
P
4
β
2
β
3
b
2
1
-
u
1
+
P
1
β
1
1
-
u
2
s
s
s
s
s
s
s
s
s
s
s
s
s
s
×
γ
v
+
β
3
I
h
*
+
ϵ
1
u
3
-
1
,
where
(5)
P
1
=
γ
v
+
δ
v
+
ϵ
2
u
3
,
P
2
=
μ
h
+
δ
h
+
γ
h
+
u
2
,
P
3
=
μ
h
+
λ
h
.
The reproduction number for the control system (2) is given by
(6)
K
u
=
b
1
μ
h
β
2
β
3
b
2
1
-
u
1
P
1
P
2
γ
v
-
ϵ
1
u
3
+
β
1
1
-
u
2
P
2
.
The reproduction number for without control system is given by
(7)
K
0
=
b
1
μ
h
β
2
β
3
b
2
γ
v
P
0
1
P
0
2
+
β
1
P
0
2
,
where P
0
1
=
(
δ
v
+
γ
v
) and P
0
2
=
(
μ
h
+
δ
h
+
γ
h
). In the above expression for the endemic equilibria the infected component I
h
* is nonzero. Using the value of S
h
*, R
h
*, and I
v
* in the first equation of the system (2), we obtained
(8)
f
I
h
*
=
a
I
h
*
2
+
b
I
h
*
+
c
=
0
,
where
(9)
a
=
λ
h
β
3
2
P
1
β
1
(
1
-
u
2
)
,
b
=
λ
h
β
3
2
β
2
b
2
1
-
u
1
+
P
4
b
1
β
1
P
1
1
-
u
2
β
3
b
=
+
λ
h
β
3
P
1
γ
v
+
ϵ
1
u
3
β
1
1
-
u
2
b
=
-
P
2
P
4
β
2
1
-
u
1
β
3
β
2
b
2
-
β
1
1
-
u
2
P
1
β
3
b
=
+
u
1
λ
h
P
1
β
3
,
c
=
-
P
1
P
2
P
4
γ
v
+
ϵ
1
u
3
μ
h
β
1
1
-
u
2
-
u
1
λ
h
1
-
K
u
.
Here the coefficient a is positive always and c depends upon the value of K
u; if the value of K
u
<
1, then c is positive; otherwise negative. The positive solution of the above equation depends upon the value of b and c. For the value of K
u
>
1, the above equation leads to two roots, positive and negative. If we substitute K
u
=
1, then the equation has no positive solution. This is possible if and only if b
<
0. For b
<
0 and K
u
=
1, the equilibria depend upon K
u; then there exists an open interval having two positive roots, that is, I
1
=
-
b
-
b
2
-
4
a
c
/
2
a and I
2
=
-
b
+
b
2
-
4
a
c
/
2
a. For c
>
0 either b
2
<
4
a
c or b
≥
0, then the above have no positive solution.For backward bifurcation, we setb
2
-
4
a
c
=
0 and K
u
=
R
c and solve for the critical value of K
c, which is given by
(10)
R
c
=
1
-
b
2
4
a
P
1
P
2
P
4
γ
v
+
ϵ
1
u
3
μ
h
β
1
1
-
u
2
-
u
1
λ
h
.
It is proved by simulating the set of parameter values presented in Table 1. Figure 1 shows the bifurcation of the control system (2). The high bold black line shows the bifurcation at u
1
=
u
2
=
u
3
=
0, that is, without control. The small black bold line shows the bifurcation at u
1
=
u
2
=
u
3
≠
0, that is, with control. The change occurs in the bifurcation rapidly when we apply the control variable. u
1 is the first control, which is represented by the bold dotted line, the dotted dashed line represents the second control, which shows the change in the control, the very narrow dotted line represent the third control, and bifurcation occurs for the control variable.Table 1
Parameter values used in the control bifurcation.
Notation
Value
b
2
10
β
1
0.048
β
3
0.006
u
2
0.001
λ
h
0.068
γ
h
0.041
b
1
10
β
2
0.0008
u
1
0.541
u
3
0.081
μ
h
0.091
δ
h
0.10
δ
v
0.009
ϵ
2
0.009
ϵ
1
0.00001
γ
v
0.087Figure 1
The plot shows the backward bifurcation for control and without control system.


## 2.4. Existence of Control Problem

We use the bounded Lebesgue measurable control and define our objective functional as(11)
J
(
u
1
,
u
2
,
u
3
)
=
min
⁡
∫
0
T
+
1
2
W
1
u
1
2
+
W
2
u
2
2
+
W
3
u
3
2
A
1
S
h
(
t
)
+
A
2
I
h
(
t
)
+
B
1
S
v
(
t
)
+
B
2
I
v
(
t
)
sssssssssssssss
+
1
2
W
1
u
1
2
+
W
2
u
2
2
+
W
3
u
3
2
d
t
.
Here, in (11), A
1, A
2, B
1, B
2, W
1, W
2, and W
3 represent the weight/balance factors just to keep the balance of individuals in the objective functional. The control set is defined as
(12)
U
=
u
=
u
1
,
u
2
,
u
3
∣
u
i
t
is Lebesgue measureable
,
s
s
s
s
0
≤
u
i
t
≤
1
,
t
∈
0
,
T
,
for
i
=
1,2
,
3
.Note. The details of the control variables u
1, u
2, and u
3 are available in Section 2.2.First, we show the existence for the control system (2). Let S
h
(
t
), I
h
(
t
), R
h
(
t
), S
v
(
t
), and I
v
(
t
) be the state variables with control variables u
1
(
t
), u
2
(
t
), and u
3
(
t
). For existence we consider the control system (2). Then we can write the system (2) in the following form:
(13)
X
′
=
A
X
+
F
X
,
where(14)
X
′
=
S
h
t
I
h
t
R
h
t
S
v
t
I
v
t
,
A
=
-
μ
h
0
λ
h
0
0
0
-
μ
h
-
δ
h
-
γ
h
-
u
2
t
0
0
0
u
1
γ
h
+
u
2
t
-
μ
h
-
λ
h
0
0
0
0
0
-
γ
v
-
ϵ
1
u
3
t
0
0
0
0
0
-
γ
v
-
δ
v
-
ϵ
2
u
3
t
,
F
X
=
b
1
-
β
2
S
h
t
I
v
1
-
u
1
t
-
β
1
S
h
t
I
h
t
1
-
u
2
t
β
2
S
h
t
I
v
1
-
u
1
t
+
β
1
S
h
t
I
h
t
1
-
u
2
t
1
-
u
2
t
0
b
2
-
β
3
S
v
t
I
h
t
β
3
S
v
t
I
h
t
,where X
′ denotes the derivative with respect to time t. The system (13) is a nonlinear system with a bounded coefficient. We set
(15)
G
(
X
)
=
A
X
+
F
(
X
)
.
The second term on the right-hand side of (15) satisfies
(16)
F
(
X
1
)
-
F
(
X
2
)
≤
C
1
S
1
h
t
-
S
2
h
t
+
C
2
I
1
h
t
-
I
2
h
t
ssssissss
+
C
3
R
1
h
t
-
R
2
h
t
+
C
4
S
1
v
t
-
S
2
v
t
ssssssss
+
C
5
I
1
v
t
-
I
2
v
t
≤
C
S
1
h
t
-
S
2
h
t
+
I
1
h
t
-
I
2
h
t
sssssss
+
R
1
h
t
-
R
2
h
t
+
S
1
v
t
-
S
2
v
t
sssissss
+
I
1
v
t
-
I
2
v
t
,
where the positive constant C
=
max
⁡
(
C
1
,
C
2
,
C
3
,
C
4
,
C
5
) is independent of the state variables. Also we have
(17)
G
X
1
-
G
X
2
≤
C
X
1
-
X
2
,
where C
=
C
1
+
C
2
+
C
3
+
C
4
+
C
5
+
M
<
∞. So, it follows that the function G is uniformly Lipschitz continuous. From the definition of control variables and nonnegative initial conditions we can see that a solution of the system (13) exists; see [23]. Now, we consider the control system (2) with the initial conditions (3) to show the existence of the control problem. Note that for bounded Lebesgue measurable controls and nonnegative initial conditions, nonnegative bounded solutions to the state system exist [23]. Let us go back to the optimal control problem (2)-(3). In order to find an optimal solution, first we will find the Lagrangian L and Hamiltonian H for the optical control problem (2)-(3). The Lagrangian for our control problem is given by
(18)
L
=
A
1
S
h
+
A
2
I
h
+
B
1
S
v
+
B
2
I
v
+
1
2
W
1
u
1
2
+
W
2
u
2
2
+
W
3
u
3
2
.
For the minimum value of Lagrangian, we define the Hamiltonian for the control problem:
(19)
H
=
L
S
h
,
I
h
,
S
v
,
I
v
,
u
1
,
u
2
,
u
3
+
λ
1
d
S
h
d
t
+
λ
2
d
I
h
d
t
+
λ
3
d
R
h
d
t
+
λ
4
d
S
v
d
t
+
λ
5
d
I
v
d
t
.
For the existence of our control problem, we state and prove the following theorem.Theorem 1.
There exists an optimal controlu
*
=
(
u
1
*
,
u
2
*
,
u
3
*
)
∈
U such that
(20)
J
u
1
*
,
u
2
*
,
u
3
*
=
min
⁡
u
1
,
u
2
,
u
3
∈
U
J
u
1
,
u
2
,
u
3
,
subject to the control system (2) with the initial conditions (3).Proof.
To prove the existence of an optimal control, we use the result in [24]; the control and the state variable are nonnegative values. In this minimizing of the problem, the necessary convexity of the objective functional in u
1, u
2, and u
3 is satisfied. The set of control variables (
u
1
,
u
2
,
u
3
)
∈
U is also convex and closed by the definition. The optimal system is bounded which determines the compactness needed for the existence of optimal control. The integrand in the objective functional (11) is A
1
S
(
t
)
+
A
2
I
(
t
)
+
B
1
S
v
+
B
2
I
v
+
1
/
2
(
W
1
u
1
2
+
W
2
u
2
2
+
W
3
u
3
2
) which is convex in the control set U. Also we can easily see that there exist a constant σ
>
1 and positive numbers ω
1 and ω
2 such that
(21)
J
u
1
,
u
2
,
u
3
≥
ω
1
u
1
2
+
u
2
2
+
u
3
2
σ
/
2
-
ω
2
,
which is the existence of an optimal control problem.To find the optimal solution, we apply Pontryagin’s maximum principle [25] given by the following.If(
x
,
u
) is an optimal solution for an optimal control problem, then there exists a nontrivial vector function λ
=
(
λ
1
,
λ
2
,
…
,
λ
n
) which satisfies the following inequalities:
(22)
d
x
d
t
=
∂
H
t
,
x
,
u
,
λ
∂
λ
,
0
=
∂
H
t
,
x
,
u
,
λ
∂
u
,
d
λ
d
t
=
-
∂
H
t
,
x
,
u
,
λ
∂
x
.
Now we apply the necessary conditions to the Hamiltonian H in (19).Theorem 2.
Suppose thatS
h
*, I
h
*, R
h
*, S
v
* and I
v
* are the optimal state solutions with associated optimal control variables (
u
1
*
,
u
2
*
,
u
3
*
) for the optimal control problem (2)-(3). Then there exist adjoint variables λ
i, for i
=
1,2
,
…
,
5, satisfying
(23)
d
λ
1
d
t
=
λ
1
-
λ
2
β
2
I
v
1
-
u
1
+
β
1
I
h
1
-
u
2
d
λ
1
d
t
=
+
λ
1
μ
h
-
u
1
λ
3
-
A
1
,
d
λ
2
d
t
=
(
λ
1
-
λ
2
)
β
1
S
h
(
1
-
u
2
)
+
μ
h
λ
2
+
δ
h
λ
2
+
γ
h
λ
2
d
λ
2
d
t
=
+
u
2
λ
2
-
λ
3
γ
h
-
λ
3
u
2
+
β
3
S
v
λ
4
-
λ
5
-
A
2
,
d
λ
3
d
t
=
λ
3
-
λ
1
λ
h
+
μ
h
λ
3
,
d
λ
4
d
t
=
λ
4
-
λ
5
β
3
I
h
+
λ
4
γ
v
+
λ
4
ϵ
1
u
3
-
B
1
,
d
λ
5
d
t
=
λ
1
-
λ
2
β
2
S
h
1
-
u
1
d
λ
5
d
t
=
-
λ
5
μ
v
+
λ
5
γ
v
+
λ
5
δ
v
+
λ
5
ϵ
2
u
3
-
B
2
with transversality conditions (or boundary conditions)
(24)
λ
i
T
=
0
,
i
=
1,2
,
…
,
5
.
Furthermore, optimal controls u
1
*, u
2
*, and u
3
* are given by
(25)
u
1
*
=
max
⁡
min
⁡
λ
2
-
λ
1
β
2
S
h
*
I
v
*
-
λ
3
S
h
*
W
1
,
1
,
0
,
u
2
*
=
max
⁡
min
⁡
λ
2
-
λ
1
β
1
S
h
*
I
h
*
+
λ
2
-
λ
3
I
h
*
W
2
,
1
,
0
,
u
3
*
=
max
⁡
min
⁡
λ
4
ϵ
1
S
v
*
+
λ
5
ϵ
2
I
v
*
W
3
,
1
,
0
.Proof.
To find the adjoint equations and the transversality conditions, we use the Hamiltonian (19). By setting S
h
t
=
S
h
*
t, I
h
t
=
I
h
*
t, R
h
t
=
R
h
*
t, S
v
t
=
S
v
*
t, and I
v
t
=
I
v
*
t and differentiating the Hamiltonian (19) with respect to S
h, I
h, R
h, S
v, and I
v, respectively, we get (23). Then solving the equations ∂
H
/
∂
u
1
=
0, ∂
H
/
∂
u
2
=
0, and ∂
H
/
∂
u
3
=
0 on the interior of the control set and then using the optimality conditions and also the property of control space U, we can derive (25).
Here we call formulas (25) for the characterization of the optimal control. The optimal control and the state are determined by solving the optimality system, which consist of the state system (1), the adjoint system (23), initial conditions at (3), boundary conditions (24), and the characterization of the optimal control (
u
1
*
,
u
2
*
,
u
3
*
) which is given by (25). In addition, the second derivative of the Lagrangian with respect to u
1, u
2, and u
3, respectively, is positive, which shows the minimum of the optimal control u
1
*, u
2
*, and u
3
*. Substituting the values of u
1
*, u
2
*, and u
3
* in the control system (2), we obtain the following system:
(26)
d
S
h
*
d
t
=
b
1
-
μ
h
S
h
*
-
β
2
S
h
*
I
v
*
×
1
-
max
⁡
⁡
min
⁡
⁡
λ
2
-
λ
1
β
2
S
h
*
I
v
*
-
λ
3
S
h
*
W
1
,
1
,
0
-
β
1
S
h
*
I
h
*
1
-
max
⁡
min
⁡
λ
2
-
λ
1
β
1
S
h
*
I
h
*
+
λ
2
-
λ
3
I
h
*
W
2
,
1
,
0
+
λ
h
R
h
*
,
d
I
h
*
d
t
=
β
2
S
h
*
I
v
*
1
-
max
⁡
⁡
min
⁡
⁡
λ
2
-
λ
1
β
2
S
h
*
I
v
*
-
λ
3
S
h
*
W
1
,
1
,
0
-
μ
h
I
h
*
-
δ
h
I
h
*
+
β
1
S
h
*
I
h
*
×
1
-
max
⁡
⁡
min
⁡
λ
2
-
λ
1
β
1
S
h
*
I
h
*
+
λ
2
-
λ
3
I
h
*
W
2
,
1
,
0
-
γ
h
I
h
*
-
max
⁡
min
⁡
λ
2
-
λ
1
β
1
S
h
*
I
h
*
+
λ
2
-
λ
3
I
h
*
W
2
,
1
,
0
I
h
*
,
d
R
h
*
d
t
=
γ
h
I
h
*
-
μ
h
R
h
*
-
λ
h
R
h
*
+
max
⁡
⁡
min
⁡
λ
2
-
λ
1
β
2
S
h
*
I
v
*
-
λ
3
S
h
*
W
1
,
1
,
0
S
h
*
+
max
⁡
⁡
min
⁡
⁡
λ
2
-
λ
1
β
1
S
h
*
I
h
*
+
λ
2
-
λ
3
I
h
*
W
2
,
1
,
0
I
h
*
,
d
S
v
*
d
t
=
b
2
-
γ
v
S
v
*
-
β
3
S
v
*
I
h
*
-
ϵ
1
S
v
*
max
⁡
⁡
min
⁡
⁡
λ
4
ϵ
1
S
v
*
+
λ
5
ϵ
2
I
v
*
W
3
,
1
,
0
,
d
I
v
*
d
t
=
β
3
S
v
*
I
h
*
+
μ
v
I
v
*
-
γ
v
I
v
*
-
δ
v
I
v
*
-
ϵ
2
I
v
*
max
⁡
⁡
min
⁡
⁡
λ
4
ϵ
1
S
v
*
+
λ
5
ϵ
2
I
v
*
W
3
,
1
,
0
with Hamiltonian H
* at (
t
,
S
h
*
,
I
h
*
,
R
h
*
,
S
v
*
,
I
v
*
,
u
1
*
,
u
2
*
,
u
3
*
,
λ
1
,
λ
2
,
λ
3
,
λ
4
,
λ
5
):
(27)
H
*
=
A
1
S
h
*
+
A
2
I
h
*
+
B
1
S
v
*
+
B
2
I
v
*
+
1
2
W
1
max
⁡
⁡
min
⁡
⁡
λ
2
-
λ
1
β
2
S
h
*
I
v
*
-
λ
3
S
h
*
W
1
,
1
,
0
2
+
W
2
×
max
⁡
⁡
min
⁡
⁡
λ
2
-
λ
1
β
1
S
h
*
I
h
*
+
λ
2
-
λ
3
I
h
*
W
2
,
1
,
0
2
+
W
3
×
max
⁡
min
⁡
λ
4
ϵ
1
S
v
*
+
λ
5
ϵ
2
I
v
*
W
3
,
1
,
0
2
+
λ
1
d
S
h
*
d
t
+
λ
2
d
I
h
*
d
t
+
λ
3
d
R
h
*
d
t
+
λ
4
d
S
v
*
d
t
+
λ
5
d
I
v
*
d
t
.


## 3. Numerical Results and Discussion

In this section, we present the numerical simulation for the proposed model with and without control system. The proposed model (1) and (2) is solved by the Runge-Kutta order 4 scheme and then we compare the control system with and without control. The optimal strategy is obtained by solving the state system and the adjoint system with the transversality conditions. In our numerical simulations, we first start to solve the control system (2) by using the Runge-Kutta order 4 scheme forward method in time and then solve the state equations and the adjoint system with the backward method previous iterations and the value from the characterization of system (25). The weight constants used in the objective functional are A
1
=
0.001, A
2
=
0.30, B
1
=
100, B
2
=
0.03, W
1
=
100, W
2
=
140, and W
3
=
120. The values of the parameters are given in Table 2.Table 2
Parameter values used in the numerical simulations of the optimal control.
Notation
Parameter description
Value
Reference
b
1
Recruitment rate for human
1.6
[12]
λ
h
Proportionality constant
0.066
[11]
μ
h
Natural death rate of human
4.6 × 10−5
[12]
γ
v
Natural death rate of vector
1.8 × 10−3
[12]
δ
h
Death rate due to disease at human class
1.0 × 10−5
[26]
γ
h
Recovery rate of the infection
2.7 × 10−3
[26]
δ
v
Death rate at infected class of vector
1.82 × 10−3
[11]
b
2
Birth rate of vector
1.9
Assumed
β
2
Transmission betweenS
h and I
v
0.0089
Assumed
β
3
Transmission betweenS
v and I
h
0.0079
Assumed
β
1
Transmission coefficient betweenS
h and I
h
0.00013
Assumed
γ
v
Natural death rate of vector
0.0027
[12]Figures2 and 3 represent the susceptible and infected human population. In these plots the dotted line represents the control system and the bold line shows the system of without control. In Figure 2 the dotted line shows the control in the population of susceptible human and the bold line shows the population of susceptible human with control. The dotted line decreases sharply than that of the bold line. Figure 3 represents the plot for the population of infected human for both the systems with and without control. The bold line shows the population of infected human without control and the dotted line shows the population of infected human with control. The dotted line decreases sharply as compared to the bold line; it means that the population of infected individuals decrease in the control system.Figure 2
The plot shows the population of susceptible human with and without control.Figure 3
The plot shows the population of infected human with and without control.In Figure4 the plot represents the population of recovered individuals of the two systems. The bold line shows the population of the system of without control and the dotted shows the population of the system with control. The population of recovered individuals decrease sharply as compared to the system of without control. Figures 5 and 6 represent the population of susceptible vector and infected vector, respectively. Figure 5 represents the population of the vector and in two systems, with and without control. The dotted line in Figures 5 and 6 shows the population of the vector in the control system. The dotted line in the Figures 5 and 6 shows the population of the vector without control system. The population of susceptible vector in the Figure 5 decreases a little bit than that of the dotted line. The population of infected vector in the Figure 6 decreases sharply than that of the dotted line.Figure 4
The plot shows the population of recovered human with and without control.Figure 5
The plot shows the population of susceptible vector with and without control.Figure 6
The plot shows the population of infected vector with and without control.Figure 7
The plot shows the control variableu
1.Figure 8
The plot shows the control variableu
2.Figure 9
The plot shows the control variableu
3.


## 4. Conclusion

In this paper, we studied the interaction of two nonlinear systems of which one is human and the other one is vector. The theoretical studies for the optimal control problem and their numerical simulation are presented in the paper. We used the optimal control strategies to minimize the infected human, infected vector and to maximize the population of susceptible human. The model which is developed from the numerical simulations of the optimality system showed that the population of infected human, infected vector decreases and the population of susceptible human increases. We also showed that with certain values of control rates there exist their corresponding optimal solutions.



---

*Source: 619035-2014-11-25.xml*