# Multiple Control Strategies for Prevention of Avian Influenza Pandemic

Authors: Roman Ullah; Gul Zaman; Saeed Islam
Journal: The Scientific World Journal
 (2014)

---

## Abstract

We present the prevention of avian influenza pandemic by adjusting multiple control functions in the human-to-human transmittable avian influenza model. First we show the existence of the optimal control problem; then by using both analytical and
numerical techniques, we investigate the cost-effective control effects for the prevention of transmission of disease. To do this, we use three control functions, the effort to reduce the number of contacts with human infected with mutant avian influenza, the antiviral treatment of infected individuals, and the effort to reduce the number of infected birds. We completely characterized the optimal control and compute numerical solution of the optimality system by using an iterative method.


---

## Body

## 1. Introduction

Highly pathogenic avian influenza is a zoonotic disease caused by H5N1 virus and it has a devastating impact on poultry causing above 90 percent mortality within 48 hours of infection [1]. The disease also spreads to humans through the direct contact with the infected poultry. Despite all advances of medical, avian influenza still poses a significant threat to health and economy of the society. The laboratory confirmed that the fatality rate in human case is about 60 percent [2].Several reports have so far been made of possible coinfection of humans with an H5N1 strain and a human strain [3, 4]. One of the coinfection reports is of a Turkish family cluster. In late December 2005 and early January 2006, a cluster of 8 confirmed highly pathogenic avian influenza (H5N1) cases was detected in Dogubayazit district in eastern Turkey [5]. All confirmed case patients were hospitalized after onset of symptoms. Four of them were died and the other four recovered. Another cluster of 8 cases of HPAI (H5N1) was detected in the village of Kubu Sambilang, Samatra, Indonesia, in April-May 2006 [6]. Only one of them recovered and the remaining seven members died. Though the human-to-human transmission of the disease is rare, its potential to change into an extremely high virulent human-to-human transmittable pandemic strain is the real danger to human health. Due to this threat, some international organizations have developed and began to implement different strategies (biosecurity, isolation, and antiviral treatment) for forestalling the onset of a pandemic. It has been estimated that a future avian influenza pandemic could cause hundreds of billions of dolors due to lost productivity, associated medical policies, and prevention policies.Mathematical models have been widely used to evaluate the effects of control measures of avian influenza [7–9]. In 2007, Iwami et al. [10] proposed a mathematical model of the spread of avian influenza from bird population to human population. They discussed that to minimize the spread of the disease in human population, we must take the measures for the infected human with bird flu to quarantine when mutant bird flu has already occurred. Lee et al. [11] show that the optimal control strategies that rely on the use of antiviral treatment or isolation strategies can reduce the number of clinical cases; therein they emphasize isolation strategies in mitigation of pandemic of influenza particularly when the access to antiviral resources is limited. Later on Ullah et al. [12] proposed the control strategies where they focus on minimizing the impact of influenza by minimizing the vaccine wanning.In this paper, we focus on identifying the optimal control strategies for the model developed in [13]; these controls minimize the impact of avian influenza by isolating the infected individuals, the judicious use of drug supply, and eliminating the infected birds. Recommended medicines like oseltamivir or zanamivir must be taken within two days after the appearance of symptoms. The clinically ill and infectious individuals can be isolated by reasonably effective ways to reduce the transmission of disease, like to educate them to cover their sneeze and cough, not to spit openly, and to avoid the closed contacts with others, sanitizing the rooms or equipments occupied by the patients. The diagnosed infected birds should be culled within a few hours and all the infected dead birds and other contaminated objects (faeces, blood, and feathers) must be destroyed properly (burying, incineration) as soon as possible. We can get the most effective results if we use the isolation, treatment of clinically infectious humans, and elimination of infected birds concurrently.


## 2. Model Frame Work

In this section, the compartmental model that divides the human and birds populations into two different classes presents the optimal control problem for the transmission dynamics of avian influenza. Our main aim is to show that it is possible to implement the time-dependent antiavian influenza control techniques while minimizing the cost of such measures. In order to do this, we introduce three time-dependent control functionsu
1, u
2, and u
3. The control u
1 represents the effort to reduce the number of contacts with human infected by mutant avian influenza, u
2 represents the fraction of clinically infectious cases treated with antiviral per unit of time, and u
3 is the effort to reduce the number of infected birds. Note that the controls are fully effective when u
i
=
1 for i
=
1,2
,
3, while there is no control if u
i
=
0. We divide the total human population at time t denoted by N
h into five distinct epidemiological subclasses which are susceptible X
h, exposed E
h, clinically ill and infectious I
h, treated T
h, and recovered R
h and the birds population N
b into three distinct subclasses which are susceptible X
b, exposed E
b, and infected I
b. Taking into account the assumptions above, the dynamics of the control problem is given by
(1)
d
X
h
d
t
=
Λ
-
μ
h
X
h
-
(
1
-
u
1
)
α
1
X
h
I
h
-
α
2
X
h
I
b
+
ϵ
h
R
h
,
d
E
h
d
t
=
(
1
-
u
1
)
α
1
X
h
I
h
+
α
2
X
h
I
b
-
(
μ
h
+
ϕ
h
)
E
h
,
d
I
h
d
t
=
ϕ
h
E
h
-
(
ρ
h
+
β
h
+
μ
h
+
u
2
)
I
h
,
d
T
h
d
t
=
ρ
h
I
h
-
(
γ
h
+
μ
h
)
T
h
,
d
R
h
d
t
=
γ
h
T
h
+
u
2
I
h
-
(
μ
h
+
ϵ
h
)
R
h
,
d
X
b
d
t
=
π
-
μ
b
X
b
-
(
1
-
u
3
)
α
3
X
b
I
b
,
d
E
b
d
t
=
(
1
-
u
3
)
α
3
X
b
I
b
-
(
μ
b
+
ϕ
b
)
E
b
,
d
I
b
d
t
=
ϕ
b
E
b
-
(
β
b
+
μ
b
)
I
b
,
with the initial conditions
(2)
X
h
(
0
)
≥
0
,
E
h
(
0
)
≥
0
,
I
h
(
0
)
≥
0
,
T
h
(
0
)
≥
0
,
R
h
(
0
)
≥
0
,
X
b
(
0
)
≥
0
,
E
b
(
0
)
≥
0
,
I
b
(
0
)
≥
0
.The human population is recruited at a constant birth rateΛ, ϵ
h is the rate of immunity loss, ϕ
h is the progression rate from E
h to I
h, ρ
h represents the treatment rate of human, μ
h and μ
b are the natural death rates of humans and birds, respectively, β
h and β
b are disease induced death rates in humans and birds, respectively, π represents recruitment rate of birds population, γ
h shows recovery due to treatment, ϕ
b is progression rate of birds from E
b to I
b, and α
1, α
2, and α
3 are effective contact rates between X
h and I
h, X
h and I
b, and X
b and I
b, respectively.The objective of our control problem is to minimize the number of clinically infectious humans and infected birds and the cost of implementing the control by using possible minimal control variablesu
i
(
t
) for i
=
1,2
,
3. We use the Lebesgue measurable control and define our objective functional as
(3)
J
(
u
1
,
u
2
,
u
3
)
=
∫
0
t
end
‍
(
A
1
I
h
+
A
2
I
b
+
1
2
(
C
1
u
1
2
+
C
2
u
2
2
+
C
3
u
3
2
)
)
d
t
,
where A
1 and A
2 are positive constants to keep a balance in the size of I
h and I
b. The square of the control variables reflects the severity of the side effects of the controls. C
1, C
2, and C
3 are positive weight parameters such that 0
<
C
1
<
Λ
/
μ
h, 0
<
C
2
<
Λ
/
μ
h, and 0
<
C
3
<
π
/
μ
b. The objective of the optimal control problem is to seek optimal control functions (
u
1
*
(
t
)
,
u
2
*
(
t
)
,
u
3
*
(
t
)
) such that
(4)
J
(
u
1
*
,
u
2
*
,
u
3
*
)
=
min
⁡
(
u
1
,
u
2
,
u
3
)
∈
U
{
J
(
u
1
,
u
2
,
u
3
)
∣
(
u
1
,
u
2
,
u
3
)
∈
U
}
,
where the control set is defined as
(5)
U
=
{
u
=
(
u
1
,
u
2
,
u
3
)
∣
u
i
is
Lebesgue
measurable
on
⁡
[
0,1
]
,
a
a
a
0
≤
u
i
(
t
)
≤
1
,
t
∈
[
0
,
t
end
]
,
for
i
=
1,2
,
3
}
,
subject to the system (1) with appropriate initial conditions. Pontryagin’s maximum principle [14] is used to solve this optimal control problem and the derivation of the necessary conditions. First we prove the existence of an optimal control for the system (1) and then derive the optimality system.


## 3. Existence of Control Problem

In this section, we consider the control system (1) with initial conditions at t
=
0 to show the existence of the control problem. Note that, for the bounded Lebesgue measurable controls and nonnegative initial conditions, nonnegative bounded solutions to the state system exist (see [15]). To prove the existence of the solution of system (1), we can write
(6)
Y
t
=
B
Y
+
F
(
Y
)
,
where(7)
Y
=
[
X
h
E
h
I
h
T
h
R
h
X
b
E
b
I
b
]
,
F
(
Y
)
=
[
Λ
-
(
1
-
u
1
)
α
1
X
h
I
h
-
α
2
X
h
I
b
(
1
-
u
1
)
α
1
X
h
I
h
+
α
2
X
h
I
b
0
0
0
π
-
(
1
-
u
3
)
α
3
X
b
I
b
(
1
-
u
3
)
α
3
X
b
I
b
0
]
,
B
=
[
-
μ
h
0
0
0
ϵ
h
0
0
0
0
-
μ
h
-
ϕ
h
0
0
0
0
0
0
0
ϕ
h
-
ρ
h
-
β
h
-
μ
h
-
u
2
0
0
0
0
0
0
0
ρ
h
-
γ
h
-
μ
h
0
0
0
0
0
0
0
γ
h
+
u
2
-
μ
h
-
ϵ
h
0
0
0
μ
b
0
0
0
0
0
0
0
0
0
0
0
0
0
-
μ
b
-
ϕ
b
0
0
0
0
0
0
0
ϕ
b
-
β
b
-
μ
b
]
,where Y
t denotes derivative of Y with respect to time t. Equation (3) is a nonlinear system with a bounded coefficient. We set
(8)
D
(
Y
)
=
B
Y
+
F
(
Y
)
.
Now,
(9)
F
(
Y
1
)
-
F
(
Y
2
)
=

[
(
1
-
u
1
)
α
1
(
X
h
2
I
h
2
-
X
h
1
I
h
1
)
+
α
2
(
X
h
2
I
b
2
-
X
h
1
I
b
1
)
-
(
1
-
u
1
)
α
1
(
X
h
2
I
h
2
-
X
h
1
I
h
1
)
-
α
2
(
X
h
2
I
b
2
-
X
h
1
I
b
1
)
0
0
0
(
1
-
u
3
)
α
3
(
X
b
2
I
b
2
-
X
b
1
I
b
1
)
-
(
1
-
u
3
)
α
3
(
X
b
2
I
b
2
-
X
b
1
I
b
1
)
0
]

,
|
F
(
Y
1
)
-
F
(
Y
2
)
|
=
|
(
1
-
u
1
)
α
1
(
X
h
2
I
h
2
-
X
h
1
I
h
1
)
+
α
2
(
X
h
2
I
b
2
-
X
h
1
I
b
1
)
|
+
|
(
1
-
u
1
)
α
1
(
X
h
2
I
h
2
-
X
h
1
I
h
1
)
+
α
2
(
X
h
2
I
b
2
-
X
h
1
I
b
1
)
|
+
|
(
1
-
u
3
)
α
3
(
X
b
2
I
b
2
-
X
b
1
I
b
1
)
|
+
|
(
1
-
u
3
)
α
3
(
X
b
2
I
b
2
-
X
b
1
I
b
1
)
|
≤
M
[
|
X
h
2
-
X
h
1
|
+
|
I
h
2
-
I
h
1
|
+
|
X
b
2
-
X
b
1
|
+
|
I
b
2
-
I
b
1
|
]
,
where
(10)
M
=
max
⁡
{
2
[
(
1
-
u
1
)
α
1
Λ
μ
h
+
α
2
π
μ
b
]
,
a
a
a
a
a
a
a
a
2
[
(
1
-
u
3
)
α
3
π
μ
b
+
α
2
Λ
μ
h
]
}
.
Here the constant M is independent of the state variables X
h, I
h, R
h, X
b, and I
b. Also we get |
D
(
Y
1
)
-
D
(
Y
2
)
|
≤
V
|
Y
1
-
Y
2
|, where V
=
max
⁡
{
M
,
∥
B
∥
}, which shows that D is uniformly Lipschitz continuous; hence by definition, the solution of the system exists (see [15]).Let us go back to the optimal control problem (1)–(3). In order to find an optimal solution, first we should find the Lagrangian and Hamiltonian for the optimal control problem. The Lagrangian is
(11)
L
=
A
1
I
h
+
A
2
I
b
+
1
2
(
C
1
u
1
2
+
C
2
u
2
2
+
C
3
u
3
2
)
.
We seek for the minimal value of the Lagrangian. To do this, we define the Hamiltonian H for the control problem, where λ
i, i
=
1,2
,
…
,
8, are the adjoint variables:
(12)
H
=
A
1
I
h
+
A
2
I
b
+
1
2
(
C
1
u
1
2
+
C
2
u
2
2
+
C
3
u
3
2
)
+
λ
1
[
Λ
-
μ
h
X
h
-
(
1
-
u
1
)
α
1
X
h
I
h
-
α
2
X
h
I
b
+
ϵ
h
R
h
]
+
λ
2
[
(
1
-
u
1
)
α
1
X
h
I
h
+
α
2
X
h
I
b
-
(
μ
h
+
ϕ
h
)
E
h
]
+
λ
3
[
ϕ
h
E
h
-
(
ρ
h
+
β
h
+
μ
h
+
u
2
)
I
h
]
+
λ
4
[
ρ
h
I
h
-
(
γ
h
+
μ
h
)
T
h
]
+
λ
5
[
γ
h
T
h
+
u
2
I
h
-
(
μ
h
+
ϵ
h
)
R
h
]
+
λ
6
[
π
-
μ
b
X
b
-
(
1
-
u
3
)
α
3
X
b
I
b
]
+
λ
7
[
(
1
-
u
3
)
α
3
X
b
I
b
-
(
μ
b
+
ϕ
b
)
E
b
]
+
λ
8
[
ϕ
b
E
b
-
(
β
b
+
μ
b
)
I
b
]
.For the existence of the control system (1), we state and prove the following theorem.Theorem 1.
There exists an optimal controlu
*
=
(
u
1
*
,
u
2
*
,
u
3
*
)
∈
U such that
(13)
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 (1) with the initial conditions (2).Proof.
To prove the existence of an optimal control, we use the result in [16–18]. Note that the control and the state variables are non-negative values. In this minimizing problem, the necessary convexity of the objective functional in u
1, u
2, and u
3 is satisfied. The set of all the control variables (
u
1
,
u
2
,
u
3
)
∈
U is also convex and closed by definition. The optimal system is bounded which determines the compactness needed for the existence of an optimal control. In addition the integrand in the functional (3), A
1
I
h
+
A
2
I
b
+
(
1
/
2
)
(
C
1
u
1
2
+
C
2
u
2
2
+
C
3
u
3
2
), is convex on the control set U. Also we can see that there exist a constant ρ
>
1 and positive numbers ω
1 and ω
2 such that
(14)
J
(
u
1
,
u
2
,
u
3
)
≥
ω
1
(
|
u
1
|
2
+
|
u
2
|
2
+
|
u
3
|
2
)
ρ
/
2
-
ω
2
,
because the state variables are bounded, which complete the existence of an optimal control.In order to derive the necessary conditions, we use Pontryagin’s maximum principle [14] as follows. If (
x
,
u
) is an optimal solution of an optimal control problem, then there exists a nontrivial vector function λ
=
(
λ
1
,
λ
2
,
…
,
λ
n
) satisfying the following equations:
(15)
d
x
d
t
=
∂
H
(
t
,
x
,
u
,
λ
)
∂
λ
,
0
=
∂
H
(
t
,
x
,
u
,
λ
)
∂
u
,
d
λ
d
t
=
∂
H
(
t
,
x
,
u
,
λ
)
∂
x
.
We now derive the necessary conditions that optimal control functions and corresponding states must satisfy. In the following theorem, we present the adjoint system and control characterization.Theorem 2.
LetX
h
*, E
h
*, I
h
*, T
h
*, R
h
*, X
b
*, E
b
*, and I
b
* be optimal state solutions with associated optimal control variables (u
1
*, u
2
*, and u
3
*) for the optimal control problem (1)–(3). Then there exist adjoint variables λ
i, for i
=
1,2
,
…
,
8, satisfying
(16)
d
λ
1
d
t
=
(
λ
1
-
λ
2
)
(
1
-
u
1
)
(
α
1
I
h
+
α
2
I
b
)
+
λ
1
μ
h
,
d
λ
2
d
t
=
(
λ
2
-
λ
3
)
ϕ
h
+
λ
2
μ
h
,
d
λ
3
d
t
=
(
λ
1
-
λ
2
)
(
1
-
u
1
)
α
1
X
h
+
(
λ
3
-
λ
4
)
ρ
h
+
λ
3
(
β
h
+
μ
h
+
u
2
)
-
λ
5
u
2
-
A
1
,
d
λ
4
d
t
=
(
λ
4
-
λ
5
)
γ
h
+
λ
4
μ
h
,
d
λ
5
d
t
=
(
λ
5
-
λ
1
)
ϵ
h
+
λ
5
μ
h
,
d
λ
6
d
t
=
(
λ
6
-
λ
7
)
(
1
-
u
3
)
α
3
I
b
+
λ
6
μ
b
,
d
λ
7
d
t
=
(
λ
7
-
λ
8
)
ϕ
b
+
λ
7
μ
b
,
d
λ
8
d
t
=
(
λ
6
-
λ
7
)
(
1
-
u
3
)
α
3
X
b
+
λ
8
(
β
b
+
μ
b
)
-
A
2
with transversality conditions
(17)
λ
i
(
t
end
)
=
0
,
i
=
1,2
,
…
,
8
.
Furthermore the control functions u
1
*, u
2
*, and u
3
* are given by
(18)
u
1
*
=
max
⁡
{
min
⁡
{
(
λ
2
-
λ
1
)
α
1
X
h
*
I
h
*
c
1
,
1
}
,
0
}
,
u
2
*
=
max
⁡
{
min
⁡
{
(
λ
3
-
λ
5
)
I
h
*
c
2
,
1
}
,
0
}
,
u
3
*
=
max
⁡
{
min
⁡
{
(
λ
7
-
λ
6
)
α
3
X
b
*
I
b
*
c
3
,
1
}
,
0
}
.Proof.
To determine the adjoint equations and the transversality conditions, we use the HamiltonianH in (12). The adjoint system results from Pontryagin’s maximum principle [14]:
(19)
d
λ
1
d
t
=
-
∂
H
∂
X
h
,
d
λ
2
d
t
=
-
∂
H
∂
E
h
,
d
λ
3
d
t
=
-
∂
H
∂
I
h
,
d
λ
4
d
t
=
-
∂
H
∂
T
h
,
d
λ
5
d
t
=
-
∂
H
∂
R
h
,
d
λ
6
d
t
=
-
∂
H
∂
X
b
,
d
λ
7
d
t
=
-
∂
H
∂
E
b
,
d
λ
8
d
t
=
-
∂
H
∂
I
b
with λ
i
(
t
end
)
=
0, i
=
1,2
,
…
,
8.
To get the characterization of the optimal control given by (18), solving the equations
(20)
∂
H
∂
u
1
=
0
,
∂
H
∂
u
2
=
0
,
∂
H
∂
u
3
=
0
,
on the interior of the control set, and setting the property of the control space U, we can derive the desired characterization (18).


## 4. Numerical Results and Discussion

In this section, we present a semi-implicit finite difference method by discretizing the interval[
t
0
,
t
f
] at the points t
i
=
t
0
+
i
l, (
i
=
0,1
,
…
,
n
), where l represents the time step such that t
n
=
t
f. We define the state and adjoint variables X
h, E
h, I
h, T
h, R
h, X
b, E
b, I
b, λ
1, λ
2, λ
3, λ
4, λ
5, λ
6, λ
7, and λ
8 and the controls u
1, u
2, and u
3 in terms of nodal points X
h
i, E
h
i, I
h
i, T
h
i, R
h
i, X
b
i, E
b
i, I
b
i, λ
1
i, λ
2
i, λ
3
i, λ
4
i, λ
5
i, λ
6
i, λ
7
i, λ
8
i, u
1
i, u
2
i, and u
3
i. By combination of forward and backward difference approximation, the method developed by [19] to adopt it in our case is given by
(21)
X
h
i
+
1
-
X
h
i
l
=
Λ
-
μ
h
X
h
i
+
1
-
(
1
-
u
1
i
)
α
1
X
h
i
+
1
I
h
i
-
α
2
X
h
i
+
1
I
b
i
+
ϵ
h
R
h
i
,
E
h
i
+
1
-
E
h
i
l
=
(
1
-
u
1
i
)
α
1
X
h
i
+
1
I
h
i
+
α
2
X
h
i
+
1
I
b
i
-
(
μ
h
+
ϕ
h
)
E
h
i
+
1
,
I
h
i
+
1
-
I
h
i
l
=
ϕ
h
E
h
i
+
1
-
(
ρ
h
+
β
h
+
μ
h
+
u
2
i
)
I
h
i
+
1
,
T
h
i
+
1
-
T
h
i
l
=
ρ
h
I
h
i
+
1
-
(
γ
h
+
μ
h
)
T
h
i
+
1
,
R
h
i
+
1
-
R
h
i
l
=
γ
h
T
i
+
1
+
u
2
i
I
h
i
+
1
-
(
μ
h
+
ϵ
h
)
R
h
i
+
1
,
X
b
i
+
1
-
X
b
i
l
=
π
-
μ
b
X
b
i
+
1
-
(
1
-
u
3
i
)
α
3
X
b
i
+
1
I
b
i
,
E
b
i
+
1
-
E
b
i
l
=
(
1
-
u
3
i
)
α
3
X
b
i
+
1
I
b
i
-
(
μ
b
+
ϕ
b
)
E
b
i
+
1
,
I
b
i
+
1
-
I
b
i
l
=
ϕ
b
E
b
i
+
1
-
(
β
b
+
μ
b
)
I
b
i
+
1
.
To approximate the time derivative of the adjoint variables by the first-ordered backward difference, we use the appropriate scheme as follows:
(22)
λ
1
n
-
i
-
λ
1
n
-
i
-
1
l
=
(
λ
1
n
-
i
-
1
-
λ
2
n
-
i
)
(
1
-
u
1
i
)
(
α
1
I
h
i
+
1
+
α
2
I
b
i
+
1
)
+
λ
1
n
-
i
-
1
μ
h
,
λ
2
n
-
i
-
λ
2
n
-
i
-
1
l
=
(
λ
2
n
-
i
-
1
-
λ
3
n
-
i
)
ϕ
h
+
λ
2
n
-
i
-
1
μ
h
,
λ
3
n
-
i
-
λ
3
n
-
i
-
1
l
=
(
λ
1
n
-
i
-
1
-
λ
2
n
-
i
-
1
)
(
1
-
u
1
i
)
α
1
X
h
i
+
1
+
λ
3
n
-
i
-
1
(
ρ
h
+
β
h
+
μ
h
+
u
2
i
)
-
λ
5
n
-
i
u
2
i
-
λ
4
n
-
i
ρ
h
-
A
1
,
λ
4
n
-
i
-
λ
4
n
-
i
-
1
l
=
(
λ
4
n
-
i
-
1
-
λ
5
n
-
i
)
γ
h
+
λ
4
n
-
i
-
1
μ
h
,
λ
5
n
-
i
-
λ
5
n
-
i
-
1
l
=
(
λ
5
n
-
i
-
1
-
λ
1
n
-
i
-
1
)
ϵ
h
+
λ
5
n
-
i
-
1
μ
h
,
λ
6
n
-
i
-
λ
6
n
-
i
-
1
l
=
(
λ
6
n
-
i
-
1
-
λ
7
n
-
i
)
(
1
-
u
3
i
)
α
3
I
b
i
+
1
+
λ
6
n
-
i
-
1
μ
b
,
λ
7
n
-
i
-
λ
7
n
-
i
-
1
l
=
(
λ
7
n
-
i
-
1
-
λ
8
n
-
i
)
ϕ
b
+
λ
7
n
-
i
-
1
μ
b
,
λ
8
n
-
i
-
λ
8
n
-
i
-
1
l
=
(
λ
6
n
-
i
-
1
-
λ
7
n
-
i
-
1
)
(
1
-
u
3
i
)
α
3
X
b
i
+
1
+
λ
8
n
-
i
-
1
(
β
b
+
μ
b
)
-
A
2
.
The algorithm that describes the approximation method for obtaining the optimal control is as follows.Algorithm 3.
Consider the following.
Step 1. Consider X
h
(
0
)
=
X
h
0, E
h
(
0
)
=
E
h
0, I
h
(
0
)
=
I
h
0, T
h
(
0
)
=
T
h
0, R
h
(
0
)
=
R
h
0, X
b
(
0
)
=
X
b
0, E
b
(
0
)
=
E
b
0, I
b
(
0
)
=
I
b
0, u
1
(
0
)
=
u
2
(
0
)
=
u
3
(
0
)
=
0, and λ
i
(
t
f
)
=
0, i
=
1
,
…
,
8.
Step 2. For i
=
1
,
…
,
n
-
1, do
(23)
X
h
i
+
1
=
X
h
i
+
l
[
Λ
+
ϵ
h
R
h
i
]
1
+
l
[
μ
h
+
(
1
-
u
1
i
)
α
1
I
h
i
+
α
2
I
b
i
]
,
E
h
i
+
1
=
E
h
i
+
l
[
(
1
-
u
1
i
)
α
1
X
h
i
+
1
I
h
i
+
α
2
X
h
i
+
1
I
b
i
]
1
+
l
[
μ
h
+
ϕ
h
]
,
I
h
i
+
1
=
I
h
i
+
l
ϕ
h
E
h
i
+
1
1
+
l
[
ρ
h
+
β
h
+
μ
h
+
u
2
i
]
,
T
h
i
+
1
=
T
h
i
+
l
ρ
h
I
h
i
+
1
1
+
l
[
γ
h
+
μ
h
]
,
R
h
i
+
1
=
R
h
i
+
l
[
γ
T
h
i
+
1
+
u
2
i
I
h
i
+
1
]
1
+
l
[
μ
h
+
ϵ
h
]
,
X
b
i
+
1
=
X
b
i
+
l
π
1
+
l
[
μ
b
+
(
1
-
u
3
i
)
α
3
I
b
i
]
,
E
b
i
+
1
=
E
b
i
+
l
(
1
-
u
3
i
)
α
3
X
b
i
+
1
I
b
i
1
+
l
[
μ
b
+
ϕ
b
]
,
I
b
i
+
1
=
I
b
i
+
l
ϕ
b
E
b
i
+
1
1
+
l
[
β
b
+
μ
b
]
,
λ
1
n
-
i
-
1
=
λ
1
n
-
i
+
l
λ
2
n
-
i
(
1
-
u
1
i
)
[
α
1
I
h
i
+
1
+
α
2
I
h
i
+
1
]
1
+
l
[
(
1
-
u
1
i
)
α
1
I
h
i
+
1
+
α
2
I
h
i
+
1
+
μ
h
]
,
λ
2
n
-
i
-
1
=
λ
3
n
-
i
+
l
λ
3
n
-
i
ϕ
h
1
+
l
λ
2
n
-
i
-
1
[
ϕ
h
+
μ
h
]
,
λ
3
n
-
i
-
1
=
λ
3
n
-
i
+
l
[
(
λ
2
n
-
i
-
1
-
λ
1
n
-
i
-
1
)
(
1
-
u
1
i
)
α
1
X
h
i
+
1
a
a
a
a
a
a
a
a
a
a
+
ρ
h
λ
4
n
-
i
+
u
2
i
λ
5
n
-
i
+
A
1
]
×
(
1
+
l
[
ρ
h
+
β
h
+
μ
h
+
u
2
i
]
)
-
1
,
λ
4
n
-
i
-
1
=
λ
4
n
-
i
-
l
γ
h
λ
n
-
i
1
-
l
[
γ
h
+
μ
h
]
,
λ
5
n
-
i
-
1
=
λ
5
n
-
i
+
l
ϵ
h
λ
1
n
-
i
-
1
1
+
l
[
ϵ
h
+
μ
h
]
,
λ
6
n
-
i
-
1
=
λ
6
n
-
i
+
l
(
1
-
u
3
i
)
I
b
i
+
1
λ
7
n
-
i
1
+
l
[
μ
b
+
(
1
-
u
3
i
)
α
3
I
b
i
+
1
]
,
λ
7
n
-
i
-
1
=
λ
7
n
-
i
+
l
ϕ
b
λ
8
n
-
i
1
+
l
[
ϕ
b
+
μ
b
]
,
λ
8
n
-
i
-
1
=
λ
8
n
-
i
+
l
[
(
λ
7
n
-
i
-
1
-
λ
6
n
-
i
-
1
)
(
1
-
u
3
i
)
α
3
X
b
i
+
1
+
A
2
]
1
+
l
[
β
b
+
μ
b
]
,
u
1
i
+
1
=
max
⁡
{
min
⁡
{
(
λ
2
n
-
i
-
1
-
λ
1
n
-
i
-
1
)
α
1
X
h
i
+
1
I
h
i
+
1
c
1
,
1
}
,
0
}
,
u
2
i
+
1
=
max
⁡
{
min
⁡
{
(
λ
3
n
-
i
-
1
-
λ
5
n
-
i
-
1
)
I
h
i
+
1
c
2
,
1
}
,
0
}
,
u
3
i
+
1
=
max
⁡
{
min
⁡
{
(
λ
7
n
-
i
-
1
-
λ
6
n
-
i
-
1
)
α
3
X
b
i
+
1
I
b
i
+
1
c
3
,
1
}
,
0
}
,
end for.
Step 3. For i
=
1
,
…
,
n
-
1, write X
h
*
(
t
i
)
=
X
h
i, E
h
*
(
t
i
)
=
E
h
i, I
h
*
(
t
i
)
=
I
h
i, T
h
*
(
t
i
)
=
T
h
i, R
h
*
(
t
i
)
=
R
h
i, X
b
*
(
t
i
)
=
X
b
i, E
b
*
(
t
i
)
=
E
b
i, I
b
*
(
t
i
)
=
I
b
i, u
1
*
(
t
i
)
=
u
1
i, u
2
*
(
t
i
)
=
u
2
i, and u
3
*
(
t
i
)
=
u
3
i end for.We have plotted susceptible, exposed, infected, treated, and recovered population of humans and susceptible, exposed, and infected population of birds with and without control by considering parameter values given in Table1, with initial values X
h
(
0
)
=
40, E
h
(
0
)
=
20, I
h
(
0
)
=
10, T
h
(
0
)
=
8, R
h
(
0
)
=
5, X
b
(
0
)
=
70, E
b
(
0
)
=
30, and I
b
(
0
)
=
20. In each of the given figures the undashed and the dashed lines represent the individuals without and with control, respectively. The weight constant of the objective functional A
1
=
0.035, A
2
=
0.0015, C
1
=
0.05, C
2
=
0.08, and C
3
=
0.002 is considered for numerical simulation. Figure 1 shows that the population of the susceptible humans increased after control, in Figures 2 and 3 we see that the number of the exposed and infected humans with control decreased more sharply than that without control, Figure 4 shows that per day clinically reported humans decreased after control, Figure 5 shows that the number of recovered human individuals increased after control, Figure 6 shows that the birds population increased after control, and in Figures 7 and 8 we see that the number of the exposed and infected birds decreased after control. Figure 9 shows the graphs of control variables.Table 1
Parameter values used for numerical simulations.
Notation
Parameters definition
Value
Λ
Recruitment rate of humans
1.7/day
α
1
Effective contact rate betweenX
h and I
h
0.3
α
2
Effective contact rate betweenX
h and I
b
0.15
α
3
Effective contact rate betweenX
b and I
b
0.12
π
Recruitment rate of birds
5/day
ϵ
h
Rate of immunity loss
0.005
μ
h
Natural death rate of humans
0.1/day
μ
b
Natural death rate of birds
0.3/day
γ
h
Recovery due to treatment
0.3
ρ
h
Treatment rate
0.75
β
h
Disease induced death rate in humans
0.4/day
β
b
Disease induced death rate in birds
3.7/day
ϕ
h
Progression rate fromE
h to I
h
0.45
ϕ
b
Progression rate fromE
b to I
b
0.85Figure 1
The plot represents the population of susceptible humans with and without control.Figure 2
The plot represents the population of exposed humans with and without control.Figure 3
The plot represents the population of infected humans with and without control.Figure 4
The plot represents the population of treated humans with and without control.Figure 5
The plot represents the population of recovered humans with and without control.Figure 6
The plot represents the population of susceptible birds with and without control.Figure 7
The plot represents the population of exposed birds with and without control.Figure 8
The plot represents the population of infected birds with and without control.Figure 9
The plot represents the control variablesu
1, u
2, and u
3.


## 5. Conclusion

An optimal control problem of the transmission dynamics of the human-avian influenza disease has been presented. We sought to determine three types of control functions associated with isolation and antiviral treatment of the clinically infectious individuals and elimination of the infected birds. Our control plots indicated that the number of exposed, infected, and hospitalized humans and the number of exposed, infected birds decreased in the optimal system. We developed the necessary conditions for the existence of the optimal control by using Pontryagin’s maximum principle. Using the state and adjoint system together with the characterization of the optimal control, we solved the problem numerically by an efficient numerical method based on optimal control with the estimated parameter values based on avian influenza. The results show that the control practices are very effective in reducing the incidence of infectious population.



---

*Source: 949718-2014-02-24.xml*