# Modeling the Effects of Vaccination and Treatment on Tuberculosis Transmission Dynamics

Authors: Ashenafi Kelemu Mengistu; Peter J. Witbooi
Journal: Journal of Applied Mathematics
 (2019)

---

## Abstract

In this study, we develop and analyze a deterministic mathematical model for tuberculosis (TB) transmission dynamics. The model includes vaccination for newborns and treatment for both high-risk latent and active TB patients. The stability of disease-free equilibrium point is discussed in detail. In the numerical simulation, the model parameters are estimated using reported TB incidence data in Ethiopia from the years 2003 to 2017, andR
0 is calculated as R
0
≈
2.13. Finally, the sensitivity indices of R
0 with respect to the model parameters are performed, and their corresponding graphical results are presented. Our results quantify the positive influence of vaccination and the treatment for high-risk latent and active TB patients on the control of tuberculosis.


---

## Body

## 1. Introduction

Tuberculosis (TB) is an airborne and highly infectious disease caused byMycobacterium tuberculosis. A susceptible individual is infected with the bacteria when he or she inhales the TB germs, which are released into the air when infected individuals cough, sneeze, spit or talk [1]. The fight against TB involves vaccination with Bacillus Calmette-Guérin (BCG), screening of those at high risk, early detection, and treatment of cases [2]. In most TB endemic countries, BCG vaccination is recommended for tuberculosis prevention and is usually administered shortly after birth to prevent TB in infants [3, 4].Globally, tuberculosis remains a major global health problem, and it is one of the top 10 causes of death and the leading cause of a single infectious agent (above HIV/AIDS). In 2017, an estimated 10 million people developed TB, and 1.3 million died from the disease (including 300,000 deaths among HIV-positive people) [5]. In Ethiopia, TB is still a serious public health challenge and one of the leading causes of morbidity and mortality [6]. According to the WHO report, Ethiopia is one of the 30 high-burden countries, and there were an estimated 172,000 (164 per 100,000 populations) incident cases of TB in 2017. According to the same report, there were an estimated 25,000 deaths (24 per 100,000) due to TB, excluding HIV related deaths.Mathematical models and computer simulations are inexpensive, easy to manage, relatively fast, and quite productive experimental tools. They have been widely used to examine, explain, and predict the dynamics of infectious disease transmission. Starting from the first mathematical model for TB by Waaler et al. [7], different mathematical models for tuberculosis have been formulated, analyzed, and utilized. These models apply to different types of populations such as a city or country, a school, prison, or refugee camp, see for instance, [8, 9] Also, different models would focus on various factors such as progression rate, treatment, vaccination, immigration, etc. see [10–12].For many years mathematical models have been applied to study the transmission dynamics of TB. For example, Zhao et al. [13] investigated the role of age on the transmission of TB in Mainland China and found that the BCG vaccine is useful only for younger people. They also recommended that the DOTS program should be more focused on the senior-aged group, and more attention should be given to people with latent TB in the middle-aged group. Choi et al. [14] introduced three control mechanisms: distancing, case finding, and case holding into the SEIL model in South Korea. They showed that distancing control is the most effective prevention mechanism of all. On the other hand, Moualeu et al. [15] developed a model for the transmission dynamics of TB and applied to the data for Cameroon. They identified that the combined effort that is education and chemoprophylaxis might lead to a reduction of 80% in the number of infected people in 10 years. Kim et al. [16] developed a mathematical model for TB and fitted to the Philippine data. Their result showed that applying a combination of distancing and case finding control strategies has significant potential for curtailing the spread of TB in the Philippines. Therefore the transmission dynamics of TB by mathematical models are important for proposing the best mechanisms to control the spread of TB. The purpose of this study is to develop a suitable TB dynamics model and calibrating it for Ethiopia.


## 2. Materials and Methods

### 2.1. Model Formulation

The model that we present in this paper is the modification of the models in [14, 16] by including TB vaccination for newborns. The homogeneously mixing total population at time t, of size N
t, is divided into five subclasses: susceptible S
t, Vaccinated V
t), high-risk latently infected E
t, Infectious (or active TB) I
t, and low-risk latent L
t.We assume that the recruitment rate into the population isΛ and some portion it, ε
Λ, will receive a vaccination at birth where 0
≤
ε
≤
1. The natural death rate (any death which is not due to TB) is assumed to be the same for each class and denoted by μ. The mortality due to the TB disease will happen only in the I-class with a rate δ. The efficiency of the BCG vaccine is not complete. Hence, it is assumed that some portion of vaccinated individuals will be susceptible to bacteria with a rate of θ.Susceptible individuals can be infected with TB through the transmission coefficientβ. The treatment rate for the E
-class is denoted by α
. It is assumed that the untreated portion of the E-class will develop active TB at the rate k. If treatment is administered for the I-class with a rate r, then some of them will complete their treatment correctly at a rate 1
−
p
r for 0
≤
p
≤
1. The recovered individuals are moved to the L-class because treatment cannot eradicate the TB bacteria from the body of the patients. Hence, recovered and low-risk latently infected individuals are classified into a single class of low-risk latent individuals. It is assumed that there is no permanent immunity to tuberculosis hence some of the recovered (low-risk latent) individuals can lose their immunity and become high-risk-latently infected, with the relapse rate σ. We further assume that all parameters to be used in this model are nonnegative.Given these assumptions, we have the following flow diagram, which describes the interaction between classes (Figure1).Figure 1
Flow diagram of the TB transmission model.Based on our definitions, assumptions, and interrelations between the variables, the system of ODE that describes the dynamics of TB is formulated as follows,(1)
d
S
d
t
=
1
−
ε
Λ
+
θ
V
−
β
S
I
−
μ
S
,
d
V
d
t
=
ε
Λ
−
θ
+
μ
V
,
d
E
d
t
=
β
S
I
+
p
r
I
+
σ
L
−
k
+
α
+
μ
E
,
d
I
d
t
=
k
E
−
μ
+
r
+
δ
I
,
d
L
d
t
=
1
−
p
r
I
+
α
E
−
μ
+
σ
L
,
N
t
=
S
t
+
V
t
+
E
t
+
I
t
+
L
t
.


### 2.2. Basic Properties

#### 2.2.1. Positivity of Solutions

Theorem 1.
Let the initial dataS
0
,
V
0
,
E
0
,
I
0 and L
0 be nonnegative. Then the solution set t
>
0.Proof.
Take the second equation of the model (1)(2)
d
V
t
d
t
=
ε
Λ
−
θ
+
μ
V
t
.For simplicity let us writeθ
+
μ
=
φ and ε
Λ
=
λ. Then,(3)
d
V
t
d
t
+
φ
V
t
=
λ
.Multiplying both sides of (3) by exp
φ
t gives(4)
d
V
t
d
t
exp
φ
t
+
φ
V
t
exp
φ
t
=
λ
e
x
p
φ
t
.By the product rule of the derivative we have(5)
d
V
t
d
t
exp
φ
t
+
φ
V
t
exp
φ
t
=
d
d
t
V
t
exp
φ
t
.Hence from (4), we have(6)
d
d
t
V
t
exp
φ
t
=
λ
exp
φ
t
.Integrating both sides of (6) gives(7)
V
t
=
V
0
exp
−
φ
t
+
λ
φ
1
−
exp
−
φ
t
≥
0
.Similarly taking the first equation of (1) gives(8)
d
S
t
d
t
=
1
−
ε
Λ
+
θ
V
t
−
β
I
t
−
μ
S
t
.By letting
1
−
ε
Λ
=
ϕ and β
I
t
−
μ
=
H
t
,we have(9)
d
S
t
d
t
+
H
t
S
t
=
ϕ
+
θ
V
t
.Multiply both sides of (9) by exp
∫
0
t
H
τ
d
τ gives(10)
d
S
t
d
t
exp
∫
0
t
H
τ
d
τ
+
H
t
S
t
exp
∫
0
t
H
τ
d
τ
=
ϕ
exp
∫
0
t
H
τ
d
τ
+
θ
V
t
exp
∫
0
t
H
τ
d
τ
.By the product rule of the derivative, we have(11)
d
S
t
d
t
exp
∫
0
t
H
τ
d
τ
+
H
t
S
t
exp
∫
0
t
H
τ
d
τ
=
d
d
t
S
t
exp
∫
0
t
H
τ
d
τ
.Hence(12)
d
d
t
S
t
exp
∫
0
t
H
τ
d
τ
=
ϕ
exp
∫
0
t
H
τ
d
τ
+
θ
V
t
exp
∫
0
t
H
τ
d
τ
.Integrating both sides of (12) gives(13)
S
t
exp
∫
0
t
H
τ
d
τ
−
S
0
=
ϕ
∫
0
t
exp
∫
0
τ
H
u
d
u
+
∫
0
t
θ
V
u
exp
∫
0
τ
H
u
d
u
.Then(14)
S
t
=
S
0
exp
−
∫
0
t
H
τ
d
τ
+
ϕ
∫
0
t
exp
∫
0
τ
H
u
d
u
exp
−
∫
0
t
H
τ
d
τ
+
∫
0
t
θ
V
u
exp
∫
0
τ
H
u
d
u
∫
0
t
θ
V
u
exp
∫
0
τ
H
u
d
u
≥
0
.Similarly, we can show thatE
t
,
I
t, and L
t are nonnegative.Remark. Note that, in particular, from equation (13) it follows that(15)
l
i
m
t
→
∞
V
t
=
ε
Λ
θ
+
μ
.


#### 2.2.2. Invariant Regions

Theorem 2.
With the nonnegative initial conditions, the feasible region of the model is defined by(16)
Ω
=
S
t
,
V
t
,
E
t
,
I
t
,
L
t
∈
R
+
5
S
t
+
V
t
+
E
t
+
I
t
+
L
t
≤
Λ
μ
.Proof.
The change of total population size is(17)
d
N
t
d
t
=
d
S
t
d
t
+
d
V
t
d
t
+
d
E
t
d
t
+
d
I
t
d
t
+
d
L
t
d
t
=
Λ
−
μ
N
t
−
δ
I
t
≤
Λ
−
μ
N
t
.The inequality (17) implies(18)
N
t
≤
Λ
μ
−
N
0
exp
−
μ
t
.Thus for everyt
>
0, as t
→
∞
,
N
t
≤
Λ
/
μ
.Hence an invariant region for this model is(19)
Ω
=
S
t
,
V
t
,
E
t
,
I
t
,
L
t
∈
R
+
5
S
t
+
V
t
+
E
t
+
I
t
+
L
t
≤
Λ
μ
.


## 2.1. Model Formulation

The model that we present in this paper is the modification of the models in [14, 16] by including TB vaccination for newborns. The homogeneously mixing total population at time t, of size N
t, is divided into five subclasses: susceptible S
t, Vaccinated V
t), high-risk latently infected E
t, Infectious (or active TB) I
t, and low-risk latent L
t.We assume that the recruitment rate into the population isΛ and some portion it, ε
Λ, will receive a vaccination at birth where 0
≤
ε
≤
1. The natural death rate (any death which is not due to TB) is assumed to be the same for each class and denoted by μ. The mortality due to the TB disease will happen only in the I-class with a rate δ. The efficiency of the BCG vaccine is not complete. Hence, it is assumed that some portion of vaccinated individuals will be susceptible to bacteria with a rate of θ.Susceptible individuals can be infected with TB through the transmission coefficientβ. The treatment rate for the E
-class is denoted by α
. It is assumed that the untreated portion of the E-class will develop active TB at the rate k. If treatment is administered for the I-class with a rate r, then some of them will complete their treatment correctly at a rate 1
−
p
r for 0
≤
p
≤
1. The recovered individuals are moved to the L-class because treatment cannot eradicate the TB bacteria from the body of the patients. Hence, recovered and low-risk latently infected individuals are classified into a single class of low-risk latent individuals. It is assumed that there is no permanent immunity to tuberculosis hence some of the recovered (low-risk latent) individuals can lose their immunity and become high-risk-latently infected, with the relapse rate σ. We further assume that all parameters to be used in this model are nonnegative.Given these assumptions, we have the following flow diagram, which describes the interaction between classes (Figure1).Figure 1
Flow diagram of the TB transmission model.Based on our definitions, assumptions, and interrelations between the variables, the system of ODE that describes the dynamics of TB is formulated as follows,(1)
d
S
d
t
=
1
−
ε
Λ
+
θ
V
−
β
S
I
−
μ
S
,
d
V
d
t
=
ε
Λ
−
θ
+
μ
V
,
d
E
d
t
=
β
S
I
+
p
r
I
+
σ
L
−
k
+
α
+
μ
E
,
d
I
d
t
=
k
E
−
μ
+
r
+
δ
I
,
d
L
d
t
=
1
−
p
r
I
+
α
E
−
μ
+
σ
L
,
N
t
=
S
t
+
V
t
+
E
t
+
I
t
+
L
t
.


## 2.2. Basic Properties

### 2.2.1. Positivity of Solutions

Theorem 1.
Let the initial dataS
0
,
V
0
,
E
0
,
I
0 and L
0 be nonnegative. Then the solution set t
>
0.Proof.
Take the second equation of the model (1)(2)
d
V
t
d
t
=
ε
Λ
−
θ
+
μ
V
t
.For simplicity let us writeθ
+
μ
=
φ and ε
Λ
=
λ. Then,(3)
d
V
t
d
t
+
φ
V
t
=
λ
.Multiplying both sides of (3) by exp
φ
t gives(4)
d
V
t
d
t
exp
φ
t
+
φ
V
t
exp
φ
t
=
λ
e
x
p
φ
t
.By the product rule of the derivative we have(5)
d
V
t
d
t
exp
φ
t
+
φ
V
t
exp
φ
t
=
d
d
t
V
t
exp
φ
t
.Hence from (4), we have(6)
d
d
t
V
t
exp
φ
t
=
λ
exp
φ
t
.Integrating both sides of (6) gives(7)
V
t
=
V
0
exp
−
φ
t
+
λ
φ
1
−
exp
−
φ
t
≥
0
.Similarly taking the first equation of (1) gives(8)
d
S
t
d
t
=
1
−
ε
Λ
+
θ
V
t
−
β
I
t
−
μ
S
t
.By letting
1
−
ε
Λ
=
ϕ and β
I
t
−
μ
=
H
t
,we have(9)
d
S
t
d
t
+
H
t
S
t
=
ϕ
+
θ
V
t
.Multiply both sides of (9) by exp
∫
0
t
H
τ
d
τ gives(10)
d
S
t
d
t
exp
∫
0
t
H
τ
d
τ
+
H
t
S
t
exp
∫
0
t
H
τ
d
τ
=
ϕ
exp
∫
0
t
H
τ
d
τ
+
θ
V
t
exp
∫
0
t
H
τ
d
τ
.By the product rule of the derivative, we have(11)
d
S
t
d
t
exp
∫
0
t
H
τ
d
τ
+
H
t
S
t
exp
∫
0
t
H
τ
d
τ
=
d
d
t
S
t
exp
∫
0
t
H
τ
d
τ
.Hence(12)
d
d
t
S
t
exp
∫
0
t
H
τ
d
τ
=
ϕ
exp
∫
0
t
H
τ
d
τ
+
θ
V
t
exp
∫
0
t
H
τ
d
τ
.Integrating both sides of (12) gives(13)
S
t
exp
∫
0
t
H
τ
d
τ
−
S
0
=
ϕ
∫
0
t
exp
∫
0
τ
H
u
d
u
+
∫
0
t
θ
V
u
exp
∫
0
τ
H
u
d
u
.Then(14)
S
t
=
S
0
exp
−
∫
0
t
H
τ
d
τ
+
ϕ
∫
0
t
exp
∫
0
τ
H
u
d
u
exp
−
∫
0
t
H
τ
d
τ
+
∫
0
t
θ
V
u
exp
∫
0
τ
H
u
d
u
∫
0
t
θ
V
u
exp
∫
0
τ
H
u
d
u
≥
0
.Similarly, we can show thatE
t
,
I
t, and L
t are nonnegative.Remark. Note that, in particular, from equation (13) it follows that(15)
l
i
m
t
→
∞
V
t
=
ε
Λ
θ
+
μ
.


### 2.2.2. Invariant Regions

Theorem 2.
With the nonnegative initial conditions, the feasible region of the model is defined by(16)
Ω
=
S
t
,
V
t
,
E
t
,
I
t
,
L
t
∈
R
+
5
S
t
+
V
t
+
E
t
+
I
t
+
L
t
≤
Λ
μ
.Proof.
The change of total population size is(17)
d
N
t
d
t
=
d
S
t
d
t
+
d
V
t
d
t
+
d
E
t
d
t
+
d
I
t
d
t
+
d
L
t
d
t
=
Λ
−
μ
N
t
−
δ
I
t
≤
Λ
−
μ
N
t
.The inequality (17) implies(18)
N
t
≤
Λ
μ
−
N
0
exp
−
μ
t
.Thus for everyt
>
0, as t
→
∞
,
N
t
≤
Λ
/
μ
.Hence an invariant region for this model is(19)
Ω
=
S
t
,
V
t
,
E
t
,
I
t
,
L
t
∈
R
+
5
S
t
+
V
t
+
E
t
+
I
t
+
L
t
≤
Λ
μ
.


## 2.2.1. Positivity of Solutions

Theorem 1.
Let the initial dataS
0
,
V
0
,
E
0
,
I
0 and L
0 be nonnegative. Then the solution set t
>
0.Proof.
Take the second equation of the model (1)(2)
d
V
t
d
t
=
ε
Λ
−
θ
+
μ
V
t
.For simplicity let us writeθ
+
μ
=
φ and ε
Λ
=
λ. Then,(3)
d
V
t
d
t
+
φ
V
t
=
λ
.Multiplying both sides of (3) by exp
φ
t gives(4)
d
V
t
d
t
exp
φ
t
+
φ
V
t
exp
φ
t
=
λ
e
x
p
φ
t
.By the product rule of the derivative we have(5)
d
V
t
d
t
exp
φ
t
+
φ
V
t
exp
φ
t
=
d
d
t
V
t
exp
φ
t
.Hence from (4), we have(6)
d
d
t
V
t
exp
φ
t
=
λ
exp
φ
t
.Integrating both sides of (6) gives(7)
V
t
=
V
0
exp
−
φ
t
+
λ
φ
1
−
exp
−
φ
t
≥
0
.Similarly taking the first equation of (1) gives(8)
d
S
t
d
t
=
1
−
ε
Λ
+
θ
V
t
−
β
I
t
−
μ
S
t
.By letting
1
−
ε
Λ
=
ϕ and β
I
t
−
μ
=
H
t
,we have(9)
d
S
t
d
t
+
H
t
S
t
=
ϕ
+
θ
V
t
.Multiply both sides of (9) by exp
∫
0
t
H
τ
d
τ gives(10)
d
S
t
d
t
exp
∫
0
t
H
τ
d
τ
+
H
t
S
t
exp
∫
0
t
H
τ
d
τ
=
ϕ
exp
∫
0
t
H
τ
d
τ
+
θ
V
t
exp
∫
0
t
H
τ
d
τ
.By the product rule of the derivative, we have(11)
d
S
t
d
t
exp
∫
0
t
H
τ
d
τ
+
H
t
S
t
exp
∫
0
t
H
τ
d
τ
=
d
d
t
S
t
exp
∫
0
t
H
τ
d
τ
.Hence(12)
d
d
t
S
t
exp
∫
0
t
H
τ
d
τ
=
ϕ
exp
∫
0
t
H
τ
d
τ
+
θ
V
t
exp
∫
0
t
H
τ
d
τ
.Integrating both sides of (12) gives(13)
S
t
exp
∫
0
t
H
τ
d
τ
−
S
0
=
ϕ
∫
0
t
exp
∫
0
τ
H
u
d
u
+
∫
0
t
θ
V
u
exp
∫
0
τ
H
u
d
u
.Then(14)
S
t
=
S
0
exp
−
∫
0
t
H
τ
d
τ
+
ϕ
∫
0
t
exp
∫
0
τ
H
u
d
u
exp
−
∫
0
t
H
τ
d
τ
+
∫
0
t
θ
V
u
exp
∫
0
τ
H
u
d
u
∫
0
t
θ
V
u
exp
∫
0
τ
H
u
d
u
≥
0
.Similarly, we can show thatE
t
,
I
t, and L
t are nonnegative.Remark. Note that, in particular, from equation (13) it follows that(15)
l
i
m
t
→
∞
V
t
=
ε
Λ
θ
+
μ
.


## 2.2.2. Invariant Regions

Theorem 2.
With the nonnegative initial conditions, the feasible region of the model is defined by(16)
Ω
=
S
t
,
V
t
,
E
t
,
I
t
,
L
t
∈
R
+
5
S
t
+
V
t
+
E
t
+
I
t
+
L
t
≤
Λ
μ
.Proof.
The change of total population size is(17)
d
N
t
d
t
=
d
S
t
d
t
+
d
V
t
d
t
+
d
E
t
d
t
+
d
I
t
d
t
+
d
L
t
d
t
=
Λ
−
μ
N
t
−
δ
I
t
≤
Λ
−
μ
N
t
.The inequality (17) implies(18)
N
t
≤
Λ
μ
−
N
0
exp
−
μ
t
.Thus for everyt
>
0, as t
→
∞
,
N
t
≤
Λ
/
μ
.Hence an invariant region for this model is(19)
Ω
=
S
t
,
V
t
,
E
t
,
I
t
,
L
t
∈
R
+
5
S
t
+
V
t
+
E
t
+
I
t
+
L
t
≤
Λ
μ
.


## 3. Model Analysis

### 3.1. Disease-Free Equilibrium Point and the Basic Reproduction Number

The disease-free equilibrium point of the model (1) is given by:(20)
P
0
∗
=
S
0
∗
,
V
0
∗
,
0
,
0
,
0
,whereS
0
∗
=
Λ
/
μ
θ
+
μ
1
−
ε
/
θ
+
μ, and V
0
∗
=
ε
Λ
/
θ
+
μ.The basic reproduction number (R
0) is the expected average number of new TB infections caused by a single infected individual when in contact with a completely susceptible population. We obtained R
0 by using the next-generation matrix method given in [17], this amounts to calculating the two matrices M and F, where M is the rate of transfer of individuals into and out of the infected classes and F is the rate of new infections in the compartment. Hence, by the equations we obtain,(21)
F
=
0
β
Λ
1
μ
−
ε
θ
+
μ
0
0
0
0
0
0
0
,(22)
M
=
k
+
α
+
μ
−
p
r
−
σ
−
k
r
+
δ
+
μ
0
−
α
−
1
+
p
r
μ
+
σ
.ThenR
0 is the dominant eigenvalue of the matrix F
M
−
1.Thus we haveR
0
=
k
β
Λ
θ
+
μ
1
−
ε
μ
+
σ
/
μ
θ
+
μ
μ
r
+
δ
+
μ
α
+
μ
+
σ
+
k
r
μ
1
−
p
+
δ
+
μ
μ
+
σ
.The following number,R
g, serves as an indicator for global stability of the disease-free equilibrium point:(23)
R
g
=
k
β
Λ
θ
+
μ
1
−
ε
+
η
μ
θ
+
μ
μ
1
μ
2
−
k
p
r
,withμ
1
=
α
+
μ
+
k, μ
2
=
r
+
δ
+
μ, and η
=
μ
σ
/
k
β
Λ
α
μ
1
+
1
−
p
r
k.Note that sincep
≤
1 and μ
>
0, we are guaranteed that μ
1
μ
2
−
k
p
r
>
0.Theorem 3.
For model (1), the disease-free equilibrium point P
0
∗ is globally asymptotically stable if R
g
<
1.Proof.
We follow a methodology similarly as in the stability analysis of [18, 19].ForR
g
<
1
we have(24)
k
β
Λ
θ
+
μ
1
−
ε
+
η
μ
θ
+
μ
−
μ
1
μ
2
−
k
p
r
<
0
.This can be written as(25)
k
β
Λ
θ
+
μ
1
−
ε
μ
θ
+
μ
+
k
β
Λ
η
μ
θ
+
μ
−
μ
1
μ
2
−
k
p
r
<
0
.By the Archimedean property ofR, there exists γ
0
>
0, for which(26)
k
β
Λ
θ
+
μ
1
−
ε
μ
θ
+
μ
+
γ
0
k
β
+
k
β
Λ
η
μ
θ
+
μ
−
μ
1
μ
2
−
k
p
r
<
0
.Also, we can find a numberγ
1, with 0
<
γ
1
<
μ, and such that(27)
k
β
Λ
θ
+
μ
1
−
ε
μ
θ
+
μ
+
γ
0
k
β
+
k
β
Λ
η
μ
θ
+
μ
σ
+
γ
1
γ
1
−
μ
1
μ
2
−
k
p
r
<
0
.We require an upper bound forS
t. From Remark (15) it follows that there exists t
0 such that  V
t
−
ε
Λ
/
θ
+
μ
<
γ
0 whenever t
>
t
0. Without loss of generality, we can assume that V
t
−
ε
Λ
/
θ
+
μ
<
γ
0 whenever t
>
0
.The latter inequality implies thatε
Λ
/
θ
+
μ
−
γ
0
<
V
t and consequently, that −
V
t
<
−
ε
Λ
/
θ
+
μ
+
γ
0
. Also, we have N
t
≤
Λ
/
μ.Thus for everyt
>
0,(28)
S
t
≤
N
t
−
V
t
<
Λ
μ
−
ε
Λ
θ
+
μ
+
γ
0
=
Λ
μ
θ
+
μ
1
−
ε
θ
+
μ
+
γ
0
.Now takingγ
2
=
γ
1
/
2, we introduce two constants C
0 and C
1 as follows:(29)
C
0
=
μ
2
k
−
α
σ
+
γ
1
k
μ
+
σ
,
(30)
C
1
=
σ
+
γ
2
μ
+
σ
.In particular then,C
0
>
0.Now we define a function:(31)
Q
t
=
E
t
+
C
0
I
t
+
C
1
L
t
.We prove now thatQ
˙
t is negative-definite. Note that we can write(32)
Q
˙
t
=
C
2
E
+
C
3
I
+
C
4
L
,whereC
2
=
C
0
k
−
μ
2
+
C
1
α, C
3
=
β
S
+
p
r
−
C
0
μ
1
+
1
−
p
r
C
1 , and C
4
=
σ
−
C
1
μ
+
σ.ThenC
4
=
−
γ
2
<
0.(33)
C
2
=
−
α
σ
+
γ
1
μ
+
σ
+
α
σ
+
γ
2
μ
+
σ
=
α
γ
2
−
γ
1
μ
+
σ
.Sinceγ
2
<
γ
1, it follows that C
2
<
0.We can writek
C
3 as:(34)
k
C
3
=
k
β
S
+
p
r
−
μ
1
μ
2
+
C
5
,where(35)
C
5
=
μ
1
α
σ
+
γ
1
+
r
k
1
−
p
σ
+
γ
2
μ
+
σ
≤
σ
+
γ
1
μ
+
σ
μ
1
α
+
r
k
1
−
p
=
σ
+
γ
1
k
β
Λ
σ
μ
μ
+
σ
η
.Noting also the upper bound forS
t, we obtain the following inequality(36)
k
C
3
≤
k
Λ
β
μ
θ
+
μ
θ
+
μ
1
−
ε
+
γ
0
β
+
p
r
−
μ
1
μ
2
+
k
Λ
β
μ
θ
+
μ
σ
+
γ
1
σ
η
.Therefore by the inequality (16), it follows thatk
C
3
<
0. This proves that Q
˙
t is negative-definite. Hence, Q
t is a Lyapunov function on Ω. Therefore, by LaSalle’s invariance principle [20], every solution of model (1), with any initial conditions in Ω, approaches P
0
∗ as t
→
∞, whenever R
0
<
1.The results in Theorem3 implies that for any initial size of the subpopulation of the model, TB can be eliminated from the population when R
0
<
1.


### 3.2. Existence of the Endemic Equilibrium Point

In this section, we show the existence of an endemic equilibrium point of the model (1). The endemic equilibrium point is the steady-state where the disease keeps alive in the population that is when at least one of the infected classes of the model is nonzero.Theorem 4.
IfR
0
>
1, then the model (1) has a unique positive endemic equilibrium (37)
P
∗
=
S
∗
,
V
∗
,
E
,
I
∗
,
L
∗with:(38)
S
∗
=
Λ
θ
+
μ
1
−  ε
R
0
θ
+  μ
,
(39)
V
∗
=
ε
Λ
θ
+
μ
,
(40)
E
∗
=
r
+
δ
+  μ
μ
k
β
R
0
−  1
,
(41)
I
∗
=
μ
β
R
0
−  1
,(42)
L
∗
=
μ
k
r
1
−  p
+
α
r
+
δ
+  μ
R
0
−  1
k
β
μ
+  σ
.


## 3.1. Disease-Free Equilibrium Point and the Basic Reproduction Number

The disease-free equilibrium point of the model (1) is given by:(20)
P
0
∗
=
S
0
∗
,
V
0
∗
,
0
,
0
,
0
,whereS
0
∗
=
Λ
/
μ
θ
+
μ
1
−
ε
/
θ
+
μ, and V
0
∗
=
ε
Λ
/
θ
+
μ.The basic reproduction number (R
0) is the expected average number of new TB infections caused by a single infected individual when in contact with a completely susceptible population. We obtained R
0 by using the next-generation matrix method given in [17], this amounts to calculating the two matrices M and F, where M is the rate of transfer of individuals into and out of the infected classes and F is the rate of new infections in the compartment. Hence, by the equations we obtain,(21)
F
=
0
β
Λ
1
μ
−
ε
θ
+
μ
0
0
0
0
0
0
0
,(22)
M
=
k
+
α
+
μ
−
p
r
−
σ
−
k
r
+
δ
+
μ
0
−
α
−
1
+
p
r
μ
+
σ
.ThenR
0 is the dominant eigenvalue of the matrix F
M
−
1.Thus we haveR
0
=
k
β
Λ
θ
+
μ
1
−
ε
μ
+
σ
/
μ
θ
+
μ
μ
r
+
δ
+
μ
α
+
μ
+
σ
+
k
r
μ
1
−
p
+
δ
+
μ
μ
+
σ
.The following number,R
g, serves as an indicator for global stability of the disease-free equilibrium point:(23)
R
g
=
k
β
Λ
θ
+
μ
1
−
ε
+
η
μ
θ
+
μ
μ
1
μ
2
−
k
p
r
,withμ
1
=
α
+
μ
+
k, μ
2
=
r
+
δ
+
μ, and η
=
μ
σ
/
k
β
Λ
α
μ
1
+
1
−
p
r
k.Note that sincep
≤
1 and μ
>
0, we are guaranteed that μ
1
μ
2
−
k
p
r
>
0.Theorem 3.
For model (1), the disease-free equilibrium point P
0
∗ is globally asymptotically stable if R
g
<
1.Proof.
We follow a methodology similarly as in the stability analysis of [18, 19].ForR
g
<
1
we have(24)
k
β
Λ
θ
+
μ
1
−
ε
+
η
μ
θ
+
μ
−
μ
1
μ
2
−
k
p
r
<
0
.This can be written as(25)
k
β
Λ
θ
+
μ
1
−
ε
μ
θ
+
μ
+
k
β
Λ
η
μ
θ
+
μ
−
μ
1
μ
2
−
k
p
r
<
0
.By the Archimedean property ofR, there exists γ
0
>
0, for which(26)
k
β
Λ
θ
+
μ
1
−
ε
μ
θ
+
μ
+
γ
0
k
β
+
k
β
Λ
η
μ
θ
+
μ
−
μ
1
μ
2
−
k
p
r
<
0
.Also, we can find a numberγ
1, with 0
<
γ
1
<
μ, and such that(27)
k
β
Λ
θ
+
μ
1
−
ε
μ
θ
+
μ
+
γ
0
k
β
+
k
β
Λ
η
μ
θ
+
μ
σ
+
γ
1
γ
1
−
μ
1
μ
2
−
k
p
r
<
0
.We require an upper bound forS
t. From Remark (15) it follows that there exists t
0 such that  V
t
−
ε
Λ
/
θ
+
μ
<
γ
0 whenever t
>
t
0. Without loss of generality, we can assume that V
t
−
ε
Λ
/
θ
+
μ
<
γ
0 whenever t
>
0
.The latter inequality implies thatε
Λ
/
θ
+
μ
−
γ
0
<
V
t and consequently, that −
V
t
<
−
ε
Λ
/
θ
+
μ
+
γ
0
. Also, we have N
t
≤
Λ
/
μ.Thus for everyt
>
0,(28)
S
t
≤
N
t
−
V
t
<
Λ
μ
−
ε
Λ
θ
+
μ
+
γ
0
=
Λ
μ
θ
+
μ
1
−
ε
θ
+
μ
+
γ
0
.Now takingγ
2
=
γ
1
/
2, we introduce two constants C
0 and C
1 as follows:(29)
C
0
=
μ
2
k
−
α
σ
+
γ
1
k
μ
+
σ
,
(30)
C
1
=
σ
+
γ
2
μ
+
σ
.In particular then,C
0
>
0.Now we define a function:(31)
Q
t
=
E
t
+
C
0
I
t
+
C
1
L
t
.We prove now thatQ
˙
t is negative-definite. Note that we can write(32)
Q
˙
t
=
C
2
E
+
C
3
I
+
C
4
L
,whereC
2
=
C
0
k
−
μ
2
+
C
1
α, C
3
=
β
S
+
p
r
−
C
0
μ
1
+
1
−
p
r
C
1 , and C
4
=
σ
−
C
1
μ
+
σ.ThenC
4
=
−
γ
2
<
0.(33)
C
2
=
−
α
σ
+
γ
1
μ
+
σ
+
α
σ
+
γ
2
μ
+
σ
=
α
γ
2
−
γ
1
μ
+
σ
.Sinceγ
2
<
γ
1, it follows that C
2
<
0.We can writek
C
3 as:(34)
k
C
3
=
k
β
S
+
p
r
−
μ
1
μ
2
+
C
5
,where(35)
C
5
=
μ
1
α
σ
+
γ
1
+
r
k
1
−
p
σ
+
γ
2
μ
+
σ
≤
σ
+
γ
1
μ
+
σ
μ
1
α
+
r
k
1
−
p
=
σ
+
γ
1
k
β
Λ
σ
μ
μ
+
σ
η
.Noting also the upper bound forS
t, we obtain the following inequality(36)
k
C
3
≤
k
Λ
β
μ
θ
+
μ
θ
+
μ
1
−
ε
+
γ
0
β
+
p
r
−
μ
1
μ
2
+
k
Λ
β
μ
θ
+
μ
σ
+
γ
1
σ
η
.Therefore by the inequality (16), it follows thatk
C
3
<
0. This proves that Q
˙
t is negative-definite. Hence, Q
t is a Lyapunov function on Ω. Therefore, by LaSalle’s invariance principle [20], every solution of model (1), with any initial conditions in Ω, approaches P
0
∗ as t
→
∞, whenever R
0
<
1.The results in Theorem3 implies that for any initial size of the subpopulation of the model, TB can be eliminated from the population when R
0
<
1.


## 3.2. Existence of the Endemic Equilibrium Point

In this section, we show the existence of an endemic equilibrium point of the model (1). The endemic equilibrium point is the steady-state where the disease keeps alive in the population that is when at least one of the infected classes of the model is nonzero.Theorem 4.
IfR
0
>
1, then the model (1) has a unique positive endemic equilibrium (37)
P
∗
=
S
∗
,
V
∗
,
E
,
I
∗
,
L
∗with:(38)
S
∗
=
Λ
θ
+
μ
1
−  ε
R
0
θ
+  μ
,
(39)
V
∗
=
ε
Λ
θ
+
μ
,
(40)
E
∗
=
r
+
δ
+  μ
μ
k
β
R
0
−  1
,
(41)
I
∗
=
μ
β
R
0
−  1
,(42)
L
∗
=
μ
k
r
1
−  p
+
α
r
+
δ
+  μ
R
0
−  1
k
β
μ
+  σ
.


## 4. Numerical Simulations

### 4.1. Estimation of the Model Parameters

In this subsection, we estimate the values of the parameters of the model (1) based on the existing literature and the epidemiological data of Ethiopia in the years between 2003 and 2017. The values of the parameters are summarised in Table 1, and the detailed estimation process of the parameter values is as follows.(a)
According to the World Bank report [21], the average life expectancy of Ethiopia in the years between 2003 and 2017 is 60.93 years. The natural death rate can be calculated as the inverse of life expectancy [8, 22]. Hence we estimate μ
=
0.016.
(b)
The upper limit of the total population in the absence of the disease isΛ
/
μ. Hence it is possible to take Λ as a product of μ and the average population size over the years 2003–2017. Using this formula and the WHO report [23], the average yearly recruitment rate to the population of Ethiopia is 1.4
×
10
6.
(c)
Also, from the WHO report [23], the average TB-induced death rate is approximately δ
=
0.17.
(d)
From the WHO report, the average BCG vaccination coverage in Ethiopia over the years 2003–2017 is 71.5% [24]. Therefore, we take ε
=
0.715.
(e)
According to the data obtained from WHO [25], the treatment success rate for TB in Ethiopia is 83.2%. Hence we estimate 1
−
p
=
0.832
.
(f)
On average, the BCG vaccine significantly reduces the risk of TB by 50% [26]. Therefore, θ
=
0.5.
(g)
The transmission coefficient, progression rate fromE-class to I-class, treatment coverage rate of I-class, the relapse rate from L-class to E-class, and treatment rate of E-class are obtained by fitting the yearly TB incidence data obtained from WHO [23] to the model by using fmincon MATLAB routine. Figure 2 shows the graph of TB incidence data obtained () and the estimated solid curve. The estimated values are β
=
1.646
×
10
−
7, κ
=
0.023, r
=
0.546
,
σ
=
0.0013 and α
=
0.153.
(h)
We calculate the initial number of vaccinated children as the product of the average number of newborns and the vaccination coverage, which isV
0
=
1
×
10
6.
(i)
The initial fraction for the infectious classI
0
=
3.73
×
10
5 is taken from the TB-prevalence in 2003 reported by the WHO [23].
(j)
We estimate the initial value of theE-classand the L-class from the data-fitting process. Thus, the amount of E
0 and L
0 is 16.37% and 30% out of the total population, respectively. This gives E
0
=
1.19
×
10
7 and L
0
=
2.18
×
10
7.
(k)
Finally, the initial number of the susceptible class is found fromS
0
=
N
0
−
V
0
+
E
0
+
I
0
+
L
0
=
3.75
×
10
7
.
(l)
Consequently, using these estimated parameters values we calculated the average value ofR
0 for the year 2003–2017 TB cases in Ethiopia is R
0
≈
2.13.Figure 2
The yearly reported TB incidence data () and the corresponding best fit (solid blue curve) of k
E
t.Table 1
The parameter values of the model (1).
Parameters
Description
Value
Source
N
0
Initial total population
7.25
×
10
7
[23]
S
0
The initial number of susceptible individuals
3.75
×
10
7
Estimated
V
0
The initial number of vaccinated individuals
1
×
10
6
Estimated
E
0
The initial number of high-risk latent individuals
1.19
×
10
7
Fitted
I
0
The initial number of infectious individuals
3.73
×
10
5
[23]
L
0
The initial number of low-risk latent individuals
2.18
×
10
7
Fitted
Λ
Recruitment rate
1.4
×
10
6
Estimated
β
The transmission coefficient
1.646
×
10
−
7
Fitted
ε
Vaccination coverage rate
0.715
[24]
θ
Loss of protection for vaccination
0.5
[26]
μ
The natural death rate
0.016
Estimated
k
Progression rate fromE
to I
0.023
Fitted
r
The treatment rate ofI
0.546
Fitted
1
−
p
Successful treatment rate ofI
0.832
[25]
α
Treatment rate ofE
0.153
Fitted
δ
TB induced death rate
0.17
[23]
σ
The relapse rate
0.0013
Fitted


### 4.2. Sensitivity Analysis of the Basic Reproduction Number

To determine the best strategy for reducing human mortality and morbidity due to TB, it is crucial to know the relative importance of the different factors responsible for its transmission and prevalence. Sensitivity analysis using the normalized forward sensitivity index is used to determine the relative importance of various parameters accountable for the disease transmission [27]. The basic reproduction number is the most significant threshold in studying the dynamics of infectious disease [28]. In this section, we calculate the sensitivity indices of R
0 for the model parameters by using the normalized forward sensitivity index.Definition. The normalized forward sensitivity index of R
0 which is differentiable with respect to a given parameter π, is defined by(43)
Γ
τ
R
0
=
∂
R
0
∂
π
×
π
R
0
.Using this formula, we calculate the sensitivity indices ofR
0 with respect to β
,
α
,
r
,
ε
, and p. The values of Γ
β
R
0 and Γ
p
R
0 are positive. This tells us, the total number of infected people can be decreased by reducing the contact rate of active TB infected individuals and the treatment failure rate. On the other hand, Γ
α
R
0, Γ
r
R
0, and Γ
ε
R
0 are negative. This means that TB infection can be controlled by increasing the treatment rate of active and latent TB infected individuals, and vaccination coverage rate. The values of the sensitivity indices for R
0 are summarized in Table 2.Table 2
Sensitivity indices ofR
0.
Parameters
Sensitivity index ofR
0
Corresponding % changes
β
The transmission coefficient
+
1
1%
α
Treatment rate ofE
−
0.801543
+1.25%
r
The treatment rate ofI
−
0.742021
+1.35%
ε
Vaccination coverage rate
−
0.0226732
+44.1%
p
Treatment failure rate
+0.0219702
+45.5%The most sensitive parameter forR
0 is the transmission coefficient β. The next essential parameters are the treatment rates of latent and active TB patients. Table 2 shows that to have 1% decrease in the value of R
0
, it is necessary to decrease the amount of β and p to 1
%, and 45.5
% respectively. While to have 1
% decrease in the value of R
0 it is necessary to increase the amount of ε
,
r, and α to 44.1, 1.35
, and 1.25
% respectively.


### 4.3. Numerical Simulations

In this subsection, the results of the model are simulated using ODE 45 solvers code in MATLAB programming language. The effect of epidemiological parameters on the number of active TB patients is simulated, and their impact is determined. The results of numerical simulations are displayed graphically.As shown in Figure3, for R
g
<
1, for different initial conditions the solution curve of the model converges to the disease-free equilibrium point P
0
∗, this indicates that P
0
∗ is globally asymptotically stable for R
g
<
1 and this agrees with Theorem 3 (Figure 4).Figure 3
The plot shows for different initial values ofS
,
V
,
E
,
I
,  L, the solution of the model converges to the disease-free equilibrium point P
0
∗
=
6.84
×
1
0
7
,
1.61
×
1
0
6
,
0
,
0
,  0 for R
g
=
0.65.Figure 4
The plot shows the effect of transmission coefficient on the total number of infected class.Figures5–8 show the effect of different epidemiological parameters on the number of the infected population. It is found that decreasing transmission coefficient (the contact between susceptible class and active TB patients) and minimizing the failure rate of treatment is essential to reduce the number of infected TB patient population. Also, increasing the treatment and vaccination coverage is vital to decrease the amount of infected TB patient population.Figure 5
The effect of vaccination coverage on the number of infected class.Figure 6
The effect of treatment failure rate on the number of infected class.Figure 7
The effect of treatment coverage on the number of infected class.Figure 8
The effect of treatment coverage of high-risk latent class on the number of infected class.


## 4.1. Estimation of the Model Parameters

In this subsection, we estimate the values of the parameters of the model (1) based on the existing literature and the epidemiological data of Ethiopia in the years between 2003 and 2017. The values of the parameters are summarised in Table 1, and the detailed estimation process of the parameter values is as follows.(a)
According to the World Bank report [21], the average life expectancy of Ethiopia in the years between 2003 and 2017 is 60.93 years. The natural death rate can be calculated as the inverse of life expectancy [8, 22]. Hence we estimate μ
=
0.016.
(b)
The upper limit of the total population in the absence of the disease isΛ
/
μ. Hence it is possible to take Λ as a product of μ and the average population size over the years 2003–2017. Using this formula and the WHO report [23], the average yearly recruitment rate to the population of Ethiopia is 1.4
×
10
6.
(c)
Also, from the WHO report [23], the average TB-induced death rate is approximately δ
=
0.17.
(d)
From the WHO report, the average BCG vaccination coverage in Ethiopia over the years 2003–2017 is 71.5% [24]. Therefore, we take ε
=
0.715.
(e)
According to the data obtained from WHO [25], the treatment success rate for TB in Ethiopia is 83.2%. Hence we estimate 1
−
p
=
0.832
.
(f)
On average, the BCG vaccine significantly reduces the risk of TB by 50% [26]. Therefore, θ
=
0.5.
(g)
The transmission coefficient, progression rate fromE-class to I-class, treatment coverage rate of I-class, the relapse rate from L-class to E-class, and treatment rate of E-class are obtained by fitting the yearly TB incidence data obtained from WHO [23] to the model by using fmincon MATLAB routine. Figure 2 shows the graph of TB incidence data obtained () and the estimated solid curve. The estimated values are β
=
1.646
×
10
−
7, κ
=
0.023, r
=
0.546
,
σ
=
0.0013 and α
=
0.153.
(h)
We calculate the initial number of vaccinated children as the product of the average number of newborns and the vaccination coverage, which isV
0
=
1
×
10
6.
(i)
The initial fraction for the infectious classI
0
=
3.73
×
10
5 is taken from the TB-prevalence in 2003 reported by the WHO [23].
(j)
We estimate the initial value of theE-classand the L-class from the data-fitting process. Thus, the amount of E
0 and L
0 is 16.37% and 30% out of the total population, respectively. This gives E
0
=
1.19
×
10
7 and L
0
=
2.18
×
10
7.
(k)
Finally, the initial number of the susceptible class is found fromS
0
=
N
0
−
V
0
+
E
0
+
I
0
+
L
0
=
3.75
×
10
7
.
(l)
Consequently, using these estimated parameters values we calculated the average value ofR
0 for the year 2003–2017 TB cases in Ethiopia is R
0
≈
2.13.Figure 2
The yearly reported TB incidence data () and the corresponding best fit (solid blue curve) of k
E
t.Table 1
The parameter values of the model (1).
Parameters
Description
Value
Source
N
0
Initial total population
7.25
×
10
7
[23]
S
0
The initial number of susceptible individuals
3.75
×
10
7
Estimated
V
0
The initial number of vaccinated individuals
1
×
10
6
Estimated
E
0
The initial number of high-risk latent individuals
1.19
×
10
7
Fitted
I
0
The initial number of infectious individuals
3.73
×
10
5
[23]
L
0
The initial number of low-risk latent individuals
2.18
×
10
7
Fitted
Λ
Recruitment rate
1.4
×
10
6
Estimated
β
The transmission coefficient
1.646
×
10
−
7
Fitted
ε
Vaccination coverage rate
0.715
[24]
θ
Loss of protection for vaccination
0.5
[26]
μ
The natural death rate
0.016
Estimated
k
Progression rate fromE
to I
0.023
Fitted
r
The treatment rate ofI
0.546
Fitted
1
−
p
Successful treatment rate ofI
0.832
[25]
α
Treatment rate ofE
0.153
Fitted
δ
TB induced death rate
0.17
[23]
σ
The relapse rate
0.0013
Fitted


## 4.2. Sensitivity Analysis of the Basic Reproduction Number

To determine the best strategy for reducing human mortality and morbidity due to TB, it is crucial to know the relative importance of the different factors responsible for its transmission and prevalence. Sensitivity analysis using the normalized forward sensitivity index is used to determine the relative importance of various parameters accountable for the disease transmission [27]. The basic reproduction number is the most significant threshold in studying the dynamics of infectious disease [28]. In this section, we calculate the sensitivity indices of R
0 for the model parameters by using the normalized forward sensitivity index.Definition. The normalized forward sensitivity index of R
0 which is differentiable with respect to a given parameter π, is defined by(43)
Γ
τ
R
0
=
∂
R
0
∂
π
×
π
R
0
.Using this formula, we calculate the sensitivity indices ofR
0 with respect to β
,
α
,
r
,
ε
, and p. The values of Γ
β
R
0 and Γ
p
R
0 are positive. This tells us, the total number of infected people can be decreased by reducing the contact rate of active TB infected individuals and the treatment failure rate. On the other hand, Γ
α
R
0, Γ
r
R
0, and Γ
ε
R
0 are negative. This means that TB infection can be controlled by increasing the treatment rate of active and latent TB infected individuals, and vaccination coverage rate. The values of the sensitivity indices for R
0 are summarized in Table 2.Table 2
Sensitivity indices ofR
0.
Parameters
Sensitivity index ofR
0
Corresponding % changes
β
The transmission coefficient
+
1
1%
α
Treatment rate ofE
−
0.801543
+1.25%
r
The treatment rate ofI
−
0.742021
+1.35%
ε
Vaccination coverage rate
−
0.0226732
+44.1%
p
Treatment failure rate
+0.0219702
+45.5%The most sensitive parameter forR
0 is the transmission coefficient β. The next essential parameters are the treatment rates of latent and active TB patients. Table 2 shows that to have 1% decrease in the value of R
0
, it is necessary to decrease the amount of β and p to 1
%, and 45.5
% respectively. While to have 1
% decrease in the value of R
0 it is necessary to increase the amount of ε
,
r, and α to 44.1, 1.35
, and 1.25
% respectively.


## 4.3. Numerical Simulations

In this subsection, the results of the model are simulated using ODE 45 solvers code in MATLAB programming language. The effect of epidemiological parameters on the number of active TB patients is simulated, and their impact is determined. The results of numerical simulations are displayed graphically.As shown in Figure3, for R
g
<
1, for different initial conditions the solution curve of the model converges to the disease-free equilibrium point P
0
∗, this indicates that P
0
∗ is globally asymptotically stable for R
g
<
1 and this agrees with Theorem 3 (Figure 4).Figure 3
The plot shows for different initial values ofS
,
V
,
E
,
I
,  L, the solution of the model converges to the disease-free equilibrium point P
0
∗
=
6.84
×
1
0
7
,
1.61
×
1
0
6
,
0
,
0
,  0 for R
g
=
0.65.Figure 4
The plot shows the effect of transmission coefficient on the total number of infected class.Figures5–8 show the effect of different epidemiological parameters on the number of the infected population. It is found that decreasing transmission coefficient (the contact between susceptible class and active TB patients) and minimizing the failure rate of treatment is essential to reduce the number of infected TB patient population. Also, increasing the treatment and vaccination coverage is vital to decrease the amount of infected TB patient population.Figure 5
The effect of vaccination coverage on the number of infected class.Figure 6
The effect of treatment failure rate on the number of infected class.Figure 7
The effect of treatment coverage on the number of infected class.Figure 8
The effect of treatment coverage of high-risk latent class on the number of infected class.


## 5. Conclusion

In this research, we studied a tuberculosis transmission dynamics model with vaccination and treatment for both high-risk latent and active TB infected classes. The reproduction number is calculated and the equilibrium points are described. We showed that the disease-free equilibrium pointP
0
∗ is globally asymptotically stable when R
g
<
1, so that the disease dies out. Finally, we showed that an increase in the treatment and vaccination coverage gives rise to a decrease in the number of infected TB patient population.The parameter values of the model are obtained from the existing literature and by fitting the yearly reported TB incidence cases in Ethiopia. We estimate that the basic reproduction number for TB transmission in Ethiopia isR
0
≈
2.13. The implication of this is that TB is still endemic in the country, and more emphasis should be given for the prevention of TB transmission. Using sensitivity analysis, we found that the most influential parameter is the transmission coefficient, followed by the treatment rate of latent TB and active TB infective individuals. Therefore, we propose that the effective strategy to eliminate the TB infection from Ethiopia is reducing the transmission coefficient. This will be done by increasing the isolation of infectious people. The second important strategy is to increase the treatment coverage for latent and active TB infectious individuals with proper follow up so that the failure of treatment cases decreases.



---

*Source: 7463167-2019-12-23.xml*