# Corrigendum to “Transmission Model of Hepatitis B Virus with the Migration Effect”

Authors: Muhammad Altaf Khan; Saeed Islam; Muhammad Arif; Zahoor ul Haq
Journal: BioMed Research International
 (2016)

---

## Body

## 1. Introduction

Hepatitis B is one of the major public health problems in the world. It is an infection that causes the liver disease. The aim of this note is to provide corrections to the basic reproduction numberR
0 in [1] as described in the letter to the editor [2, 3]. In view of [2, 3], the model is reconstructed and all their mathematical results are computed.


## 2. Model Correction

In view of [2, 3], we decided that the product of F and V
-
1, that is, F
V
-
1, is incorrect. If the product F
V
-
1 is incorrect, then the obtained basic reproduction number in [1] is incorrect. The incorrectness of the basic reproduction number affects the stability (local and global), which is to be addressed again. Further, in light of [2], the model is valid biologically. By observing all these assumptions in [2, 3], we decide to reformulate the model by making the following changes to the model revised; we add the parameter α, a rate from migrated class to susceptible, and ξ is the rate of flow from exposed to migrated classes to the model published in [1]. The complete flow diagram of our new improved model can be seen in Figure 1. For complete details of each class with their parameters description, we refer the reader to [1], except α and ξ. Adding the above parameters we obtain the following improved model:(1)
d
S
t
d
t
=
δ
π
1
-
η
C
t
-
δ
+
p
S
t
-
β
A
t
+
κ
C
t
S
t
+
δ
o
V
t
+
α
M
t
,
d
E
t
d
t
=
β
A
t
+
κ
C
t
S
t
-
δ
+
γ
1
+
ξ
E
t
+
δ
π
η
C
t
+
μ
1
M
t
,
d
A
t
d
t
=
γ
1
E
t
-
δ
+
γ
2
A
t
+
μ
2
M
t
,
d
C
t
d
t
=
q
γ
2
A
t
-
δ
+
γ
3
C
t
,
d
V
t
d
t
=
γ
3
C
t
+
1
-
q
γ
2
A
t
-
δ
o
+
δ
V
t
+
δ
1
-
π
+
p
S
t
,
d
M
t
d
t
=
ξ
E
t
-
μ
1
+
μ
2
+
δ
+
α
M
t
,being subject to the initial conditions(2)
S
0
≥
0
,
E
0
≥
0
,
A
0
≥
0
,
C
0
≥
0
,
V
0
≥
0
,
M
0
≥
0
.The model published in [1] should be replaced by (1) as well as the results. Assume S
(
t
)
+
E
(
t
)
+
A
(
t
)
+
C
(
t
)
+
V
(
t
)
+
M
(
t
)
=
1; we obtain the following reduced model for (1):(3)
d
S
t
d
t
=
δ
π
+
δ
o
-
δ
π
η
+
δ
o
C
t
-
δ
+
p
+
δ
o
S
t
-
β
A
t
+
κ
C
t
S
t
-
δ
o
E
t
+
A
t
+
M
t
+
α
M
t
,
d
E
t
d
t
=
β
A
t
+
κ
C
t
S
t
-
δ
+
γ
1
+
ξ
E
t
+
δ
π
η
C
t
+
μ
1
M
t
,
d
A
t
d
t
=
γ
1
E
t
-
δ
+
γ
2
A
t
+
μ
2
M
t
,
d
C
t
d
t
=
q
γ
2
A
t
-
δ
+
γ
3
C
t
,
d
M
t
d
t
=
ξ
E
t
-
μ
1
+
μ
2
+
δ
+
α
M
t
.Let (4)
Γ
=
S
,
E
,
A
,
C
,
M
∈
R
+
5
∣
S
t
≤
δ
π
+
δ
o
δ
+
δ
o
+
p
,
S
+
E
+
A
+
C
+
M
≤
δ
π
+
δ
o
δ
+
δ
o
. Here Γ is a positively invariant set. All the solutions inside Γ are our main focus.Figure 1
The complete flow diagram of hepatitis B virus transmission model.


## 3. Equilibria and Basic Reproduction Number

System (3) has the disease-free equilibrium at D
o
=
(
S
o
,
0,0
,
0,0
); that is, S
o
=
δ
π
+
δ
o
/
δ
+
δ
o
+
p. At D
1
=
(
S
∗
,
E
∗
,
A
∗
,
C
∗
,
M
∗
), the endemic equilibria of system (3) are given by(5)
S
∗
=
δ
+
γ
1
+
ξ
δ
+
γ
2
δ
+
γ
3
α
+
δ
+
μ
1
+
μ
2
-
δ
π
η
q
γ
2
ξ
μ
2
+
γ
1
α
+
δ
+
μ
1
+
μ
2
+
ξ
μ
1
δ
+
γ
2
δ
+
γ
3
β
δ
+
γ
3
+
q
κ
γ
2
ξ
μ
2
+
γ
1
α
+
δ
+
μ
1
+
μ
2
,
E
∗
=
δ
+
γ
2
α
+
δ
+
μ
1
+
μ
2
A
∗
γ
1
α
+
δ
+
μ
1
+
μ
2
+
μ
2
ξ
,
C
∗
=
q
γ
2
A
∗
δ
+
γ
3
,
M
∗
=
ξ
δ
+
γ
2
A
∗
ξ
μ
2
+
γ
1
α
+
δ
+
μ
1
+
μ
2
.The basic reproduction number R
0 for system (3) can be obtained by using the method in [4]. We obtain the following matrices: (6)
F
=
0
β
S
o
β
κ
S
o
0
0
0
0
0
0
0
0
0
0
0
0
0
,
V
=
δ
+
ξ
+
γ
1
0
-
δ
π
η
-
μ
1
-
γ
1
δ
+
γ
2
0
-
μ
2
0
-
q
γ
2
δ
+
γ
3
0
-
ξ
0
0
α
+
δ
+
μ
1
+
μ
2
.The basic reproduction number R
0 for (3) is given by (7)
R
0
=
β
S
o
δ
+
γ
3
+
κ
q
γ
2
+
δ
π
η
q
γ
2
μ
2
ξ
+
γ
1
α
+
δ
+
μ
1
+
μ
2
+
μ
1
ξ
δ
+
γ
2
δ
+
γ
3
δ
+
γ
1
+
ξ
δ
+
γ
2
δ
+
γ
3
α
+
δ
+
μ
1
+
μ
2
.Define(8)
R
1
=
β
δ
π
+
δ
o
δ
+
γ
3
+
κ
q
γ
2
+
δ
π
η
q
γ
2
δ
+
δ
o
+
p
γ
1
δ
+
γ
1
+
ξ
δ
+
γ
2
δ
+
γ
3
δ
+
δ
o
+
p
,
R
2
=
β
δ
π
+
δ
o
δ
+
γ
3
+
κ
q
γ
2
ξ
μ
2
+
ξ
δ
+
δ
o
+
p
δ
π
η
q
γ
2
μ
2
+
μ
1
δ
+
γ
2
δ
+
γ
3
δ
+
γ
1
+
ξ
δ
+
γ
2
δ
+
γ
3
δ
+
δ
o
+
p
α
+
δ
+
μ
1
+
μ
2
,where S
o
=
δ
π
+
δ
o
/
δ
+
δ
o
+
p and R
0
=
R
1
+
R
2.


## 4. Stability of Disease-Free Equilibrium (DFE)

In this section, we discuss the local and global stability of system (3) at disease-free equilibrium. We have the following results.Theorem 1.
The disease-free equilibriumD
o of system (3) is locally asymptotically stable if R
0
<
1 and the conditions (
γ
2
+
δ
)
(
γ
3
+
δ
)
>
β
S
o
γ
1 and γ
2
+
δ
>
β
S
o are satisfied; otherwise, they become unstable.Proof.
At the disease-free equilibriumD
o the corresponding Jacobian matrix J
o
(
ζ
) of system (3) is computed as follows: (9)
J
o
ζ
=
-
δ
+
δ
o
+
p
-
δ
o
-
β
S
o
+
δ
o
-
T
2
-
δ
o
+
α
0
-
δ
+
γ
1
+
ξ
β
S
o
β
κ
S
o
+
δ
π
η
μ
1
0
γ
1
-
δ
+
γ
2
0
μ
2
0
0
q
γ
2
-
δ
+
γ
3
0
0
ξ
0
0
-
μ
1
+
μ
2
+
δ
+
α
,where T
2
=
(
δ
π
η
+
β
κ
S
o
+
δ
o
). The root -
(
δ
+
δ
o
+
p
) of J
o
(
ζ
) is clearly negative; the other roots can be obtained from the following equation: (10)
λ
4
+
c
1
λ
3
+
c
2
λ
2
+
c
3
λ
+
c
4
=
0
,where (11)
c
1
=
γ
2
+
δ
+
γ
3
+
δ
+
α
+
δ
+
μ
1
+
μ
2
+
γ
1
+
δ
+
ξ
,
c
2
=
α
+
δ
+
μ
2
γ
1
+
δ
+
ξ
+
μ
1
γ
1
+
δ
+
γ
2
+
δ
γ
3
+
δ
-
β
γ
1
S
o
+
γ
2
+
δ
+
γ
3
+
δ
α
+
δ
+
μ
1
+
μ
2
+
γ
1
+
δ
+
ξ
,
c
3
=
γ
2
+
δ
+
γ
3
+
δ
μ
1
γ
1
+
δ
+
α
+
δ
+
μ
1
+
μ
2
γ
2
+
δ
γ
3
+
δ
-
β
γ
1
S
o
+
γ
2
+
δ
α
+
δ
γ
1
+
δ
+
ξ
+
μ
2
γ
1
+
δ
+
μ
2
ξ
γ
2
+
δ
-
β
S
o
+
γ
1
+
δ
+
ξ
γ
3
+
δ
α
+
δ
+
μ
2
+
γ
2
+
δ
1
-
R
1
,
c
4
=
γ
2
+
δ
γ
3
+
δ
α
+
δ
+
μ
1
+
μ
2
γ
1
+
δ
+
ξ
1
-
R
0
.The eigenvalues of the characteristics equation have negative real part if the Routh-Hurwitz condition is satisfied; that is, c
i
,
i
=
1,2
,
3,4, with c
1
c
2
c
3
-
c
3
2
-
c
1
2
c
4
>
0. Therefore, if R
0
<
1 and (
γ
2
+
δ
)
(
γ
3
+
δ
)
>
β
S
o
γ
1 and γ
2
+
δ
>
β
S
o, we get c
i
>
0 for i
=
1,2
,
3,4. Thus, at D
o the disease-free equilibrium of system (3) is locally asymptotically stable if the above conditions are satisfied.Next, we show the global stability of DFE by using the method in [5].Theorem 2.
The disease-free equilibrium of system (3) is globally asymptotically stable if R
0
<
1.Proof.
ConsiderY
1
=
S, Y
2
=
(
E
,
A
,
C
,
M
)
T, and Y
=
(
Y
1
,
Y
2
). The invariant domain Γ is clearly a compact positive set. We can present the subsystem Y
1
=
G
1
(
Y
1
,
0
)
(
Y
1
-
Y
1
∗
) as (12)
S
˙
=
-
δ
+
δ
o
+
p
S
-
δ
π
+
δ
o
δ
+
δ
o
+
p
.System (3) represents a linear system which is globally asymptotically stable at Y
1
∗
=
(
δ
π
+
δ
o
)
/
(
δ
+
δ
o
+
p
). The hypotheses in [5] are satisfied. The matrix G
2
(
Y
) is given by (13)
G
2
Y
=
-
δ
+
γ
1
+
ξ
β
S
β
κ
S
+
δ
π
η
μ
1
γ
1
-
δ
+
γ
2
0
μ
2
0
q
γ
2
-
δ
+
γ
3
0
ξ
0
0
-
μ
1
+
μ
2
+
δ
+
α
.Following the hypothesis H
3 in [5], for any Y
∈
Γ the above matrix G
2
(
Y
) is Metzler and irreducible. Further, we check the fourth condition H
4 in [5]. There is a maximum in Γ if S
=
(
δ
π
+
δ
o
)
/
(
δ
+
δ
o
+
p
). This corresponds to disease-free equilibrium and the maximum G
2
(
Y
¯
) is given by (14)
G
2
Y
¯
=
-
δ
+
γ
1
+
ξ
β
S
o
β
κ
S
o
+
δ
π
η
μ
1
γ
1
-
δ
+
γ
2
0
μ
2
0
q
γ
2
-
δ
+
γ
3
0
ξ
0
0
-
μ
1
+
μ
2
+
δ
+
α
.The last hypothesis H
5 requires that α
(
G
2
(
Y
¯
)
)
≤
0
. We write G
2
(
Y
¯
) in block form as (15)
G
2
Y
¯
=
B
1
B
2
B
3
B
4
,where (16)
B
1
=
-
δ
+
γ
1
+
ξ
,
B
2
=
β
S
o
,
β
κ
S
o
+
δ
π
η
,
μ
1
,
B
3
=
γ
1
0
ξ
,
B
4
=
-
δ
+
γ
2
0
μ
2
q
γ
2
-
δ
+
γ
3
0
0
0
-
μ
1
+
μ
2
+
δ
+
α
.Clearly B
1 is stable Metzler matrix, so we can write α
(
B
4
-
B
3
B
1
-
1
B
2
)
≤
0
. Let K
=
B
4
-
B
3
B
1
-
1
B
2. (17)
K
=
-
δ
+
γ
2
-
γ
1
β
S
o
δ
+
γ
1
+
ξ
γ
1
β
κ
S
o
+
δ
π
η
δ
+
γ
1
+
ξ
μ
2
-
γ
1
μ
1
δ
+
γ
1
+
ξ
q
γ
2
-
δ
+
γ
3
0
ξ
β
S
o
δ
+
γ
1
+
ξ
ξ
β
κ
S
o
+
δ
π
η
δ
+
γ
1
+
ξ
-
μ
1
+
μ
2
+
δ
+
α
-
ξ
μ
1
δ
+
γ
1
+
ξ
.The characteristics equation of K are given by (18)
λ
3
+
l
1
λ
2
+
l
2
λ
+
l
3
=
0
,where(19)
l
1
=
γ
1
+
δ
+
ξ
γ
2
+
δ
+
γ
3
+
δ
+
α
+
δ
+
μ
2
+
μ
1
γ
1
+
δ
-
β
γ
1
S
o
,
l
2
=
γ
2
+
δ
γ
3
+
δ
1
-
R
1
+
γ
2
+
δ
+
γ
3
+
δ
α
+
δ
+
μ
2
+
μ
1
γ
1
+
δ
γ
+
δ
+
ξ
+
β
S
o
γ
1
+
δ
+
ξ
2
2
γ
1
μ
1
ξ
-
γ
1
+
δ
+
ξ
ξ
μ
2
+
γ
1
α
+
δ
+
μ
1
+
μ
2
γ
1
+
δ
+
ξ
,
l
3
=
2
μ
1
ξ
γ
1
γ
1
+
δ
+
ξ
2
q
γ
2
δ
π
η
+
β
κ
S
o
+
γ
3
+
δ
β
S
o
+
γ
2
+
δ
γ
3
+
δ
×
α
+
δ
+
μ
1
+
μ
2
1
-
R
0
.The Routh-Hurwitz criteria ensure that the above characteristics equation has three negative eigenvalues if R
0
<
1 and the condition l
1
l
2
-
l
3
>
0. All the conditions in [5] are satisfied. Thus, we conclude that system (3) at disease-free equilibrium is globally asymptotically stable.


## 5. Stability EE

In this section, we determine the local and global stability of (3) at endemic equilibrium.Theorem 3.
IfR
0
>
1, then model (3) at endemic equilibrium is locally asymptotically stable.Proof.
The Jacobian matrixJ
∗ computed at D
1 is given by (20)
J
∗
=
-
z
1
-
δ
o
-
β
S
∗
-
δ
o
-
δ
π
η
-
β
κ
S
∗
-
δ
o
-
δ
o
+
α
z
2
-
δ
+
γ
1
+
ξ
β
S
∗
δ
π
η
+
β
κ
S
∗
μ
1
0
γ
1
-
δ
+
γ
2
0
μ
2
0
0
q
γ
2
-
δ
+
γ
3
0
0
ξ
0
0
-
α
+
δ
+
μ
1
+
μ
2
.The characteristics equation at D
1 is (21)
λ
5
+
g
1
λ
4
+
g
2
λ
3
+
g
3
λ
2
+
g
4
λ
+
g
5
=
0
,where (22)
g
1
=
4
δ
+
z
1
+
α
+
γ
1
+
γ
2
+
γ
3
+
μ
1
+
μ
2
+
ξ
,
g
2
=
α
+
δ
+
μ
1
+
μ
2
γ
1
+
δ
+
γ
2
+
δ
+
γ
3
+
δ
+
ξ
α
+
δ
+
μ
2
+
α
+
δ
+
μ
1
+
μ
2
+
γ
1
+
δ
+
ξ
+
γ
2
+
δ
+
γ
3
+
δ
z
1
+
γ
2
+
δ
γ
1
+
δ
+
ξ
+
γ
3
+
δ
+
γ
3
+
δ
γ
1
+
δ
+
ξ
-
z
3
+
δ
o
z
2
,
g
3
=
α
+
δ
+
μ
1
+
μ
2
γ
1
+
δ
+
ξ
γ
3
+
δ
+
γ
2
+
δ
+
γ
2
+
δ
γ
3
+
δ
-
δ
o
z
2
+
z
2
ξ
α
-
δ
o
-
γ
3
+
δ
δ
o
-
γ
2
+
δ
δ
o
+
γ
2
+
δ
z
1
γ
1
+
δ
+
ξ
+
γ
3
+
δ
+
α
+
δ
+
μ
1
+
μ
2
z
1
γ
1
+
δ
+
ξ
+
γ
3
+
δ
+
γ
2
+
δ
+
γ
3
+
δ
γ
1
+
δ
+
ξ
z
1
-
z
1
μ
1
ξ
+
z
3
-
γ
1
z
2
z
4
-
μ
1
ξ
γ
3
+
δ
+
γ
2
+
δ
+
γ
2
+
δ
γ
3
+
δ
γ
1
+
δ
+
ξ
-
z
5
-
z
3
α
+
δ
+
μ
1
+
μ
2
-
z
3
γ
3
+
δ
-
z
6
,
g
4
=
γ
2
+
δ
γ
3
+
δ
α
+
δ
+
μ
1
+
μ
2
γ
1
+
δ
+
ξ
-
γ
2
+
δ
δ
o
α
+
δ
+
μ
1
+
μ
2
z
2
-
γ
3
+
δ
δ
o
α
+
δ
+
μ
1
+
μ
2
z
2
+
ξ
γ
2
+
δ
z
2
α
-
δ
o
+
z
2
ξ
γ
3
+
δ
α
-
δ
o
-
μ
1
ξ
γ
2
+
δ
z
1
-
μ
1
ξ
γ
3
+
δ
z
1
-
γ
2
+
δ
γ
3
+
δ
δ
o
z
2
-
z
3
α
+
δ
+
μ
1
+
μ
2
z
1
+
γ
2
+
δ
+
γ
3
+
δ
α
+
δ
+
μ
1
+
μ
2
γ
1
+
δ
+
ξ
z
1
-
γ
1
γ
3
+
δ
z
2
z
4
-
z
1
z
5
+
γ
2
+
δ
γ
3
+
δ
α
+
δ
+
μ
1
+
μ
2
z
1
+
γ
2
+
δ
γ
3
+
δ
γ
1
+
δ
+
ξ
z
1
-
μ
2
ξ
z
2
z
4
-
z
3
γ
3
+
δ
z
1
-
z
1
z
6
-
z
1
z
5
-
γ
1
α
+
δ
+
μ
1
+
μ
2
z
2
z
4
+
μ
1
ξ
γ
2
+
δ
γ
3
+
δ
-
α
+
δ
+
μ
1
+
μ
2
z
5
-
z
7
-
z
3
γ
3
+
δ
α
+
δ
+
μ
1
+
μ
2
-
z
6
γ
3
+
δ
,
g
5
=
-
γ
2
+
δ
γ
3
+
δ
δ
o
α
+
δ
+
μ
1
+
μ
2
z
2
+
ξ
γ
2
+
δ
γ
3
+
δ
z
2
α
-
δ
o
+
γ
2
+
δ
γ
3
+
δ
α
+
δ
+
μ
1
+
μ
2
γ
1
+
δ
+
ξ
z
1
-
μ
1
ξ
γ
2
+
δ
γ
3
+
δ
z
1
-
α
+
δ
+
μ
1
+
μ
2
z
5
z
1
-
z
1
z
7
-
z
3
γ
3
+
δ
α
+
δ
+
μ
1
+
μ
2
z
1
-
z
6
z
1
-
γ
1
z
2
q
γ
2
α
+
δ
+
μ
1
+
μ
2
z
8
-
μ
2
ξ
q
γ
2
z
2
z
8
-
γ
1
γ
3
+
δ
α
+
δ
+
μ
1
+
μ
2
z
2
z
4
-
μ
2
ξ
γ
3
+
δ
z
2
z
4
,where z
1
=
(
δ
+
δ
o
+
p
)
+
β
(
A
∗
+
κ
C
∗
), z
2
=
β
(
A
∗
+
κ
C
∗
), z
3
=
β
γ
1
S
∗, z
4
=
δ
o
+
β
S
∗, z
5
=
γ
1
q
γ
2
π
δ
η
+
β
κ
S
∗, z
6
=
β
μ
2
ξ
S
∗, z
7
=
μ
2
ξ
q
γ
2
π
δ
η
+
β
κ
S
∗, and z
8
=
π
δ
η
+
δ
o
+
β
κ
S
∗. The characteristics equation will give negative roots if the following conditions are satisfied: (i)
g
1
,
g
2
,
g
3
,
g
4
,
g
5
>
0,
(ii)
g
1
g
2
-
g
3
>
0,
(iii)
g
4
g
5
g
3
g
1
g
2
-
g
3
-
g
1
2
g
4
>
0,
(iv)
g
5
2
(
2
g
1
g
4
-
g
2
(
g
1
g
2
-
g
3
)
-
g
5
)
>
0
.Conditions (i) and (ii) are easy to satisfy. If conditions (iii) and (iv) are satisfied; then the characteristics equation given above will give negative eigenvalues. Thus, it follows from Routh-Hurwitz criteria that system (3) is locally asymptotically stable at D
1.Next result shows the global stability of the endemic equilibriumD
1 of system (3) for special case when μ
1
=
0 and α
=
0. We have the following result.Theorem 4.
The endemic equilibriumD
1 of system (3) is globally asymptotically stable if condition (30) holds.Proof.
We define the Lyapunov function in the following form:(23)
V
=
w
1
S
-
S
∗
-
S
∗
log
⁡
S
∗
S
+
w
2
E
-
E
∗
-
E
∗
log
⁡
E
∗
E
+
w
3
A
+
w
4
C
+
w
5
M
.The time durative of V is(24)
V
′
=
w
1
1
-
S
∗
S
S
′
+
w
2
1
-
E
∗
E
E
′
+
w
3
A
′
+
w
4
C
′
+
w
5
M
′
.The coefficients w
1
,
…
,
w
5 are positive and will be determined later. At endemic steady state the first equation of (3) is given by(25)
δ
π
+
δ
o
=
δ
π
η
C
∗
+
δ
+
p
S
∗
+
β
A
∗
+
κ
C
∗
S
∗
+
δ
o
X
∗
.Using (25) and the S equation (3) we obtain(26)
w
1
1
-
S
∗
S
S
′
=
w
1
1
-
S
∗
S
δ
π
η
C
∗
+
δ
+
p
S
∗
+
β
A
∗
+
κ
C
∗
S
∗
+
δ
o
X
∗
-
w
1
δ
π
η
C
+
δ
+
p
S
+
β
A
+
κ
C
S
+
δ
o
X
+
w
1
S
∗
S
δ
π
η
C
t
+
δ
+
p
S
t
+
β
A
t
+
κ
C
t
S
t
+
δ
o
X
.Similarly,(27)
w
2
1
-
E
∗
E
E
′
=
w
2
β
A
t
+
κ
C
t
S
t
-
δ
+
γ
1
+
ξ
E
t
+
δ
π
η
C
t
-
w
2
E
∗
E
β
A
t
+
κ
C
t
S
t
+
w
2
β
A
∗
+
κ
C
∗
S
∗
+
w
2
δ
π
η
C
∗
+
δ
π
η
C
t
E
∗
E
.Using the values of (26) and (27) and the last three equations of system (3) and substituting in (24) we obtain(28)
V
′
=
w
1
1
-
S
∗
S
δ
π
η
C
∗
+
δ
+
p
S
∗
+
β
A
∗
+
κ
C
∗
S
∗
+
δ
o
X
∗
-
w
1
δ
π
η
C
-
w
1
δ
+
p
S
-
w
1
β
A
+
κ
C
S
-
w
1
δ
o
X
+
w
1
S
∗
S
δ
π
η
C
t
+
w
1
δ
+
p
S
+
w
1
β
A
+
κ
C
S
+
w
1
δ
o
X
+
w
2
β
A
+
κ
C
S
-
δ
+
γ
1
+
ξ
E
+
δ
π
η
C
-
w
2
E
∗
E
β
A
+
κ
C
S
+
w
2
β
A
∗
+
κ
C
∗
S
∗
+
w
2
δ
π
η
C
∗
+
w
2
δ
π
η
C
E
∗
E
.After implication we obtain(29)
V
′
=
-
q
γ
2
μ
2
ξ
+
γ
1
μ
1
+
μ
2
+
δ
+
α
μ
1
+
μ
2
+
δ
+
α
δ
+
γ
1
+
ξ
δ
+
p
S
-
S
∗
2
S
-
q
γ
2
μ
2
ξ
+
γ
1
μ
1
+
μ
2
+
δ
+
α
μ
1
+
μ
2
+
δ
+
α
δ
+
γ
1
+
ξ
β
S
∗
A
∗
×
-
2
+
S
∗
S
+
A
S
E
∗
A
∗
S
∗
E
-
q
γ
2
μ
2
ξ
+
γ
1
μ
1
+
μ
2
+
δ
+
α
μ
1
+
μ
2
+
δ
+
α
δ
+
γ
1
+
ξ
β
S
∗
κ
C
∗
-
2
+
S
∗
S
+
C
S
E
∗
C
∗
S
∗
E
-
q
γ
2
μ
2
ξ
+
γ
1
μ
1
+
μ
2
+
δ
+
α
μ
1
+
μ
2
+
δ
+
α
δ
+
γ
1
+
ξ
-
2
+
S
∗
S
+
C
C
∗
E
∗
E
-
C
C
∗
S
∗
S
δ
π
η
C
∗
-
δ
+
γ
2
δ
+
γ
3
C
-
q
γ
2
μ
2
ξ
+
γ
1
μ
1
+
μ
2
+
δ
+
α
μ
1
+
μ
2
+
δ
+
α
δ
+
γ
1
+
ξ
δ
o
X
∗
-
1
+
S
∗
S
+
X
X
∗
-
X
X
∗
S
∗
S
,where the constants w
1
,
…
,
w
5 are chosen as w
4
=
(
δ
+
γ
2
), w
3
=
q
γ
2, w
5
=
q
γ
2
μ
2
/
(
μ
1
+
μ
2
+
δ
+
α
), w
2
=
w
1
=
q
γ
2
(
μ
2
ξ
+
γ
1
(
μ
1
+
μ
2
+
δ
+
α
)
)
/
(
μ
1
+
μ
2
+
δ
+
α
)
(
δ
+
γ
1
+
ξ
), and X
=
S
+
E
+
A
+
C
+
M, X
∗
=
S
∗
+
E
∗
+
A
∗
+
C
∗
+
M
∗.
Equation (29) V
′
≤
0 if the following inequalities are satisfied:(30)
-
2
+
S
∗
S
+
C
C
∗
E
∗
E
-
C
C
∗
S
∗
S
≥
0
,
-
2
+
S
∗
S
+
C
S
E
∗
C
∗
S
∗
E
≥
0
,
-
2
+
S
∗
S
+
C
C
∗
E
∗
E
-
C
C
∗
S
∗
S
≥
0
,
-
1
+
S
∗
S
+
X
X
∗
-
X
X
∗
S
∗
S
≥
0
.The endemic equilibrium D
1 of system (3) is said to be globally asymptotically stable if condition (30) holds.


## 6. Numerical Simulation

This section deals with the numerical solution of model (3). The numerical results for model (3) are presented in Figures 2–8. In this paper, the value for α and ξ is taken between 0 and 1. The rest of the parameters values are taken from [1], except those mentioned in the figures. Figure 1 shows the population behavior of incidentals when μ
1
=
0.1 and μ
2
=
0.1. Figure 3 represents the population of incisiveness when μ
1
=
0.2 and μ
2
=
0.2. In Figures 2 and 3, the individuals of carriers decrease sharply. For μ
1
=
0.001 and μ
2
=
0.02, we present Figure 4. Figure 5 is population behavior of individuals when μ
1
=
0.1 and μ
2
=
0.2. Figures 6 and 7, respectively, represent the population of carriers individuals for different values of parameters. The population behavior of susceptible, exposed, acute, and carriers individuals for different parameters is presented in Figure 8. The numerical results from Figures 2 to 8 show when there is a decrease in the value of suggested parameters, and the population of individuals in the host decreases sharply.Figure 2
Population behavior of individuals whenμ
1
=
0.1, μ
2
=
0.1.Figure 3
Population behavior of individuals whenμ
1
=
0.2, μ
2
=
0.2.Figure 4
Population behavior of individuals whenμ
1
=
0.01,  μ
2
=
0.02.Figure 5
Population behavior of individuals whenμ
1
=
0.1, μ
2
=
0.2.Figure 6
Population behavior of carriers individuals: solid line,π
=
1, p
=
0; dashed π
=
0.1, p
=
0.1.Figure 7
Population behavior of carriers individuals: solid linep
=
0, μ
1
=
0.1, and μ
2
=
0.1; dashed p
=
0.1, μ
1
=
0.01, and μ
2
=
0.01.Figure 8
Population behavior of susceptible, exposed, acute, and carriers individuals: solid lineμ
1
=
0.2, μ
2
=
0.2; dashed μ
1
=
0.002, μ
2
=
0.002.
(a)
(b)
(c)
(d)


## 7. Conclusion

In this corrigendum, we make all the necessary changes to the published paper [1], which are highlighted in the comment papers [2, 3]. We added the parameter α, a rate from migrated class to susceptible, and ξ is the rate of flow from exposed to migrated classes. Further, the basic reproduction number has been investigated. The mathematical results for the revised model are presented successfully.



---

*Source: 9306329-2016-12-20.xml*