# Mathematical Modelling of Human African Trypanosomiasis Using Control Measures.

1. IntroductionHuman African trypanosomiasis (HAT), commonly known as sleeping sickness, is a vector-borne tropical disease which is caused by Trypanosoma brucei protozoa species. It is one of the neglected tropical diseases which affect people in sub-Saharan Africa, specifically those living in rural areas. HAT is caused by two species of protozoa which are Trypanosoma brucei gambiense (TBG), which causes the chronic form of HAT in central and western Africa, and Trypanosoma brucei rhodesiense (TBR), which causes the acute form of the disease in eastern and southern Africa [1]. The HAT disease has killed millions of people since the beginning of 20th century and it is transmitted from one individual to another by tsetse flies (genus Glossina); TBG is transmitted by riverine tsetse species, while TBR is transmitted by savanna tsetse species [1]. Rhodesiense HAT is an acute disease that can lead to death if not treated within 6 months, while gambiense HAT is a slow chronic progressive disease which causes death with an average duration of 3 years [2]. The signs and symptoms for both forms of HAT are not specific and their appearances vary from one person to another; at the first stage of HAT, the disease is not severe and the signs and symptoms such as intermittent fever, headache, pruritus, lymphadenopathies, asthenia, anemia, cardiac disorders, endocrine disturbances, musculoskeletal pains, and hepatosplenomegaly may be observed, while in the second stage of HAT, sleep disorders and neuropsychiatric disorders are likely to dominate. The HAT disease can be treated by using drugs such as suramin, eflornithine, melarsoprol, and pentamidine.

The disease is reported to affect about 37 sub-Saharan African countries; it affects much rural areas where there are suitable environments for the tsetse flies to live and reproduce, and the periurban areas can also be affected. The transmission of HAT can occur during human activities such as hunting, farming, as well as fishing [3]. The transmission of HAT needs the reservoir; reservoir is a species that can permanently maintain the pathogen and from which the pathogen can be transmitted to the target population [4]. Rhodesiense HAT is zoonotic which requires a nonhuman reservoir (animals) for maintaining its population, while in gambiense HAT, humans act as key reservoir [4].

Mathematical models have been used to study the transmission and effective control of diseases simply and cheaply with no need of expensive and complicated experiments [5]. So far, different models have been developed and formulated by different researchers. One of the important modelling work on HAT was done by Rogers [6]; the model explained the mathematical framework on transmission of HAT in multiple host populations [6]. Rogers' model was generalized by Hargrove et al [7], and a new parameter which allows the tsetse flies to feed off multiple hosts was introduced. The model compared the effectiveness of two methods used to control HAT: insecticide-treated cattle and the use of trypanocide drugs to treat cattle. They found out that treating cattle with insecticides is more effective and a cheaper approach to control HAT than using trypanocide drugs. Kajunguri [8] developed a model which was based on a constant population with a fixed number of domestic animals, human, and tsetse flies in one of the villages in West Africa. The major findings of their model estimated that the cattle population contributes to about 92% of the total TBR transmission, while the rest 8% is the contribution of human population in transmission of the disease. The study by Kajunguri [8], which also formulated a multihost model, was used to study the control of tsetse flies and TBR in southern Uganda. They found out that the effective application of insecticides brings about a cost-effective method of control and eliminating the disease. They realized that using insecticides for controlling HAT is more effective and efficient in the area where there are few wild hosts.

Due to low mortality rate of the disease and poverty of its sufferers, the efforts toward the control of HAT has reduced. Most attention is given to popular diseases such as HIV/AIDS, tuberculosis, malaria, and ebola, although the disease is still a threat to the lives of sub-Saharan African people. Moreover, very few studies have been carried out on applying optimal control theory to HAT transmission models. In this paper, we use optimal control theory to study the transmission dynamics of HAT diseases by using education, treatment, and insecticides as the control measures.

The rest of this paper is outlined as follows: Section 2 represents the vector-host model and the underlying assumptions. In section 3, the model equilibria and stabilities are determined, whereas in Section 4, the optimal control model is analyzed by modifying the previous one to control the HAT by using control measures (education, insecticides, and treatment). In addition, the numerical simulations for the optimal control model are done in this section and we use the results obtained to compare the efforts of each control measure to control the HAT in Africa. Finally, we provide the conclusion in Section 5.

2. Model Formulation

In this section, the vector-host model as well as the necessary differential equations to describe the transmission of HAT from tsetse fly to human and vice versa are developed. The transmission of HAT in the human population is modelled using four subclasses: Susceptible [S.sub.H], Exposed [E.sub.H], Infectious [I.sub.H], and Recovered [R.sub.H]. The total human population, [N.sub.H], is thus defined by

[N.sub.H] = [S.sub.H] + [E.sub.H] + [I.sub.H] + [R.sub.H]. (1)

The transmission of HAT in the vector (tsetse flies) population is also divided into Susceptible ([S.sub.V]), Exposed ([E.sub.V]), and Infectious ([I.sub.V]). The total population of the tsetse flies, [N.sub.V], is also defined by

[N.sub.V] = [S.sub.V] + [E.sub.V] + [I.sub.V]. (2)

We assume a constant population for both host and vector. It is also assumed that the tsetse fly cannot recover from the disease and the infected tsetse fly remains infectious throughout the rest of its life; there is no disease-induced death rate for tsetse flies and the recruitment rates are assumed to be constant due to birth and immigration.

In our model, the recruitment rate of hosts and vectors are represented by [[pi].sub.H] and [[pi].sub.V], respectively. The susceptible host gets the disease when bitten by infectious tsetse fly, and susceptible tsetse fly gets the disease when it bites an infectious human at the rate a. The natural mortality rate for humans and vectors are represented by [[mu].sub.H] and [[mu].sub.V], respectively. The parameter [omega] represents the disease-induced death rate for humans, while [[xi].sub.H] and [[xi].sub.V] are the force of infection for humans and vectors, respectively. The parameter [sigma] represents per capita rate of a vector becoming infectious, and the rest of the parameters are explained in Table 1. Assuming that the transmission per bite from infectious tsetse fly to human is a, then the rate of infection per susceptible human is given by

[[xi].sub.H] = [ap.sub.H][I.sub.V]/[N.sub.V], (3)

and also if we further assume that a is the tsetse-fly biting rate, that is, the average number of bites per tsetse fly per unit, then the rate of infection per susceptible tsetse fly can be represented by

[[xi].sub.V] = [ap.sub.V][I.sub.H]/[N.sub.H]. (4)

From the model diagram in Figure 1, the following differential equations are derived:

[mathematical expression not reproducible] (5)

From system (5), the dimensionless technique is used to [n'.sub.h] [less than or equal to] derive another equivalent differential equation; we denote [S.sub.h] = ([S.sub.H]/[N.sub.H]), [e.sub.h] = ([E.sub.H]/[N.sub.H]), [i.sub.h] = ([I.sub.H]/[N.sub.H]), [r.sub.h] = ([R.sub.H]/[N.sub.H]), [s.sub.v] = ([S.sub.V]/[N.sub.V]), [e.sub.v] = ([E.sub.V]/[N.sub.V]), and [i.sub.v] = ([I.sub.V]/[N.sub.V]) and substitute in system (5), to obtain the following new equivalent equations:

[mathematical expression not reproducible] (6)

Table 1 shows the description of the model parameters and variables.

2.1. Positivity and Boundedness of the Solutions. In this subsection, we show that system (6) is epidemiologically and mathematically well defined in the positive invariant region:

D = {([s.sub.h],[e.sub.h],[i.sub.h],[r.sub.h],[s.sub.v],[e.sub.v],[i.sub.v]) [member of] [R.sup.7.sub.+]: [n.sub.h] [less than or equal to] [[pi].sub.h]/[[mu].sub.h]; [n.sub.v] [less than or equal to] [[pi].sub.v]/[[mu].sub.v]}. (7)

Theorem 1. There exists a domain D in which the solution ([s.sub.h],[e.sub.h],[i.sub.h],[r.sub.h],[s.sub.v],[e.sub.v],[i.sub.v]) is contained and bounded.

Proof. We provide the proof following the idea by Olaniyi and Obabiyi [9]. Given the solution set ([s.sub.h], [e.sub.h], [i.sub.h], [r.sub.h], [s.sub.v], [e.sub.v], [i.sub.v]) with the positive initial conditions ([mathematical expression not reproducible]), we define

[n.sub.h]([s.sub.h],[e.sub.h],[i.sub.h],[r.sub.h]) = [s.sub.h](t)+ [e.sub.h](t)+ [i.sub.h](t)+ [r.sub.h](t) and [n.sub.v]([s.sub.v], [e.sub.v], [i.sub.v]) = [s.sub.v](t) + [e.sub.v](t) + [i.sub.v](t). (8)

The derivatives of [n.sub.h] and [n.sub.v] with respect to time along the solution of system (6) for human and tsetse flies, respectively, are obtained by

[mathematical expression not reproducible] (9)

From these differential equations, it follows that [n'.sub.h] [less than or equal to] [[pi].sub.h] - [[mu].sub.h][n.sub.h] and [n'.sub.v] [less than or equal to] [[pi].sub.v] - [[mu].sub.v][n.sub.v]. We obtain the solutions as follows:

[mathematical expression not reproducible] (10)

By taking the limits of both [n.sub.h] and [n.sub.v] above as t [right arrow] [infinity], we obtain [n.sub.h] [less than or equal to] ([[pi].sub.h]/[[mu].sub.h]) and [n.sub.v] [less than or equal to] ([[pi].sub.v]/[[mu].sub.v]); hence, the solutions are contained in the region D. This implies that all solutions of the human and tsetse fly population are contained in the region D and are nonnegative; this guarantees that the positive invariant region for system (6) exists and is given by

D = {([s.sub.h],[e.sub.h],[i.sub.h],[r.sub.h],[s.sub.v],[e.sub.v],[i.sub.v]) [member of] [R.sup.7.sub.+] : [n.sub.h] [less than or equal to] [[pi].sub.h]/[[mu].sub.h]; [n.sub.v] [less than or equal to] [[pi].sub.v]/[[mu].sub.v]. (11)

3. Model Equilibria and Stability Analysis

In this section, we give the model equilibria, the basic reproduction number, [R.sub.0], and the stabilities at both diseasefree and endemic equilibrium.

3.1. Disease-Free Equilibrium (DFE). The DFE in system (6) is when there are no HAT infections within the human and tsetse fly population. Thus, the existence of the DFE is given by [E.sub.0] = (([[pi].sub.h]/[[mu].sub.h]), 0,0,0, ([[pi].sub.v]/[[mu].sub.v]), 0,0).

3.2. Endemic Equilibrium (EE). The EE is the nontrivial equilibrium point at which the HAT disease persists in both human and tsetse fly population. Thus, the EE is obtained as follows: E* = ([S*.sub.h],[e*.sub.h],[i*.sub.h],[r*.sub.h],[s*.sub.v],[e*.sub.v],[i*.sub.v]), where

[mathematical expression not reproducible] (12)

and the term

[mathematical expression not reproducible]. (13)

3.3. Basic Reproduction Number, [R.sub.0]. The basic reproduction number, [R.sub.0], is defined as the number of secondary infections caused by one infected host or vector in a completely susceptible population [10]. The next-generation matrix approach as done by Van den Driessche and Watmough in [5, 11] is applied to derive

[mathematical expression not reproducible]. (14)

By denoting matrix F = ([partial derivative]F/[partial derivative][x.sub.i]) and V = ([partial derivative]V/[partial derivative][x.sub.i]), where [x.sub.i] = [e.sub.h],[i.sub.h],[e.sub.v],[i.sub.v], the spectral radius of the next-generation matrix [FV.sup.-1] gives the value of [R.sub.0]:

[mathematical expression not reproducible] (15)

The spectral radius [sigma]([FV.sup.-1]) gives

[R.sub.0] = [sigma]([FV.sup.-1]) = [square root of [[a.sub.2][epsilon][p.sub.h][p.sub.v][[pi].sub.h][[pi].sub.v][sigma]]/[[mu].sup.2.sub.v][[mu].sub.h]([epsilon]+[[mu].su b.h])([sigma]+[[mu].sub.v])([sigma]+[[mu].sub.v])([[mu].s ub.h]+[tau]+[omega]). (16)

One infected human in a population of susceptible vectors will cause [R.sub.v] infected vectors; likewise, one infected vector in a population will cause [R.sub.h] infected humans [5]. Therefore, the basic reproduction number can be rewritten as [R.sub.0] = [square root of [R.sub.h][R.sub.v], where [R.sub.h] = (a[epsilon][p.sub.h][[pi].sub.h]/[[mu].sub.h]([[mu].sub.h] + [epsilon])([[mu].sub.h] + [tau] + [omega])) and [R.sub.v] = ([sigma]a[p.sub.v][[pi].sub.v]/[[mu].sup.2.sub.v] + [[mu].sub.v])). Thus, [R.sub.0] can also be defined as the square root of the product of the number of infected humans in the susceptible population caused by one infected tsetse fly in its infectious lifetime and the number of infected tsetse flies caused by one infected human during the infectious period [12].

3.4. Local Stability of Disease-Free Equilibrium (DFE)

Theorem 2. If [R.sub.0] [less than or equal to] 1, the DFE given by [E.sub.0] is locally asymptotically stable in the region defined by (7), and it is unstable when [R.sub.0] > 1.

Proof. The DFE is locally stable if all eigenvalues of Jacobian matrix [mathematical expression not reproducible] are negative. The matrix has all eigenvalues negative only if the trace of [mathematical expression not reproducible] and determinant of [mathematical expression not reproducible]. By linearizing system (6) around [E.sub.0], we obtain the following Jacobian matrix:

[mathematical expression not reproducible] (17)

The trace of matrix [mathematical expression not reproducible] is such that

[mathematical expression not reproducible]. (18)

Using the basic properties of matrix algebra as in [13], it is clear that the eigenvalues [[lambda].sub.1] = -[[mu].sub.h] and [[lambda].sub.2] = -[[mu].sub.v] of the matrix [mathematical expression not reproducible] have negative real parts. The reduced matrix is

[mathematical expression not reproducible] (19)

From matrix [mathematical expression not reproducible], the eigenvalue [[lambda].sub.3] = -([rho] + [[mu].sub.h]) has negative real part. The remaining matrix is further reduced by using the reduction techniques, and we obtain

[mathematical expression not reproducible]. (20)

Using the properties of matrix algebra, the matrix [mathematical expression not reproducible] has eigenvalue -([epsilon] + [[mu].sub.h]) which has negative real part. We further reduce to a 2 x 2 matrix by using the same reduction techniques. The matrix is

[mathematical expression not reproducible] (21)

From the reduced 2 x 2 matrix, the trace is negative and the determinant is

[mathematical expression not reproducible] (22)

Since

[R.sub.0] = [square root of [[a.sup.2][epsilon][p.sub.h][p.sub.v][[pi].sub.h][[pi].sub.v][sigma]]/[[mu].sup.2.sub.v][[mu].sub.h]([epsilon]+[[mu].su b.h])([sigma]+[[mu].sub.v])([[mu].sub.h]+[tau]+[omega])], (23)

then, by letting [R.sub.T] = ([a.sup.2][epsilon][p.sub.h][p.sub.v][[pi].sub.h][[pi].sub.v][sigma]/[[mu].sup.2.sub.v][[mu].sub.v]) ([[mu].sub.h] + [tau] + [omega])), we find our determinant as

[mathematical expression not reproducible]. (24)

The value of [R.sub.T] can be seen to be positive because all the parameters are positive. As a result, the determinant in (24) is positive if and only if [R.sub.T] < 1. Therefore, the DFE is locally stable if [R.sub.T] [less than or equal to] 1.

3.5. Global Stability of Disease-Free Equilibrium (DFE). To show that the DFE is globally stable, we apply Lyapunov's theorem in [5].

Theorem 3. The DFE defined by [E.sub.0] is globally asymptotically stable in the region defined by (7) if [R.sub.0] [less than or equal to] 1. Otherwise unstable if [R.sub.0] > 1.

Proof. We define Lyapunov's function as

[mathematical expression not reproducible], (25)

satisfying system (6), where [k.sub.1],[k.sub.2],[k.sub.3],[k.sub.4],[k.sub.5],[k.sub.6] > 0 are to be determined and [mathematical expression not reproducible]. We first show that V > 0 for all E [not equal to] (([[pi].sub.h]/[[mu].sub.h]), 0,0,0, ([[pi].sub.v]/[[mu].sub.v]), 0,0). It is enough to check that

[mathematical expression not reproducible]. (26)

The function g(m) = m - 1 - ln m such that [mathematical expression not reproducible] has minimum value equal to zero when m = 1; hence, g (m) > 0 for all m > 0. Thus, Lyapunov's function V > 0. The function V is radially unbounded because as [absolute value of m] [right arrow] [infinity], the function g(m) [right arrow] [infinity]. We now take the derivative of V with respect to time and use system (6) to replace the derivatives in the right hand side such that

[mathematical expression not reproducible] (27)

The terms with [r.sub.h] are ignored because if [S.sub.h],[e.sub.h],[i.sub.h] are globally stable, then [r.sub.h] [right arrow] 0 at any time t and the DFE for system (6) is globally stable. Taking [k.sub.1] = [k.sub.2] = (1/[[mu].sub.h] + [epsilon]), [k.sub.4] = [k.sub.5] = (1/[[mu].sub.v] + [sigma]), [k.sub.3] = (1/[epsilon]), and [k.sub.6] = (1/[sigma]), the derivative of V with respect to time becomes

[mathematical expression not reproducible] (28)

The terms (([[pi].sub.h]/[[mu].sub.h][s.sub.h]) + ([[mu].sub.h][s.sub.h][n.sub.h]) - 2) and (([[pi].sub.v]/[[mu].sub.v][s.sub.v]) + ([[mu].sub.v][s.sub.v]/[[pi].sub.v]) - 2) are positive because if we suppose m = ([[??].sub.h]/[[mu].sub.h][s.sub.h]) = ([[pi].sub.v]/[[mu].sub.v][s.sub.v]), we have m + (1/m) - 2 = ([m.sup.2] - 2m + 1/m) = ([(m - 1).sup.2]/m) > 0 for all m > 1 and since [R.sub.v] < 1 and [R.sub.h] [less than or equal to] 1, then V' is negative. Thus, we have V' < 0 for all [E.sub.0] [not equal to] (([[pi].sub.h]/[[mu].sub.h]), 0,0,0, ([[pi].sub.v]/[[mu].sub.v]), 0,0).

Thus, the largest compact invariant set in D is the singleton set [E.sub.0]. Hence, system (6) is globally asymptotically stable.

3.6. Local Stability of Endemic Equilibrium (EE)

Theorem 4. The unique endemic equilibrium defined by E* is locally asymptotically stable in the region defined by (7) if [R.sub.0] > 1, but is unstable if [R.sub.0] [less than or equal to] 1.

Proof. We give the proof of this theorem based on the approach used by Olaniyi and Obabiyi [9, 14]. From the EE points defined in (12), since all values are positive, we express the value of [i*.sub.h] in terms of [R.sub.0] to obtain

[mathematical expression not reproducible] (29)

where

[mathematical expression not reproducible] (30)

Since the basic reproduction number [R.sub.0] = [square root of ([a.sup.2][epsilon][P.sub.h][P.sub.v][[pi].sub.h][[pi].sub.v]h[sigma]/[[mu].sup.2.sub.v][[mu].sub.h]([sigma]+[[mu].sub. v])([[mu].sub.h]+[tau]+[omega])),

If we let

[R.sub.T] = [a.sup.2][epsilon][p.sub.h][p.sub.v][[pi].sub.h][[pi].sub.v][sigma]/[[mu].sup.2.sub.v][[mu].sub.h]([epsilon]+[[mu].sub. h])([sigma]+[[mu].sub.v])([[mu].sub.h] + [tau] + [omega]), (31)

then we can find [i*.sub.h] as

[i*.sub.h] = [[mu].sup.2.sub.v][[mu].sub.h]([epsilon]+[[mu].sub.h])([rho]+[[mu].sub.h])([sigma]+[[mu].sub.v])([[mu].sub.h]+[tau]+[om ega])/B [[R.sub.T]-1]. (32)

The value of B is clearly positive because all parameters are positive. Hence, [i*.sub.h] > 0 if and only if [R.sub.T] > 1, implying that the EE is locally asymptotically stable if [R.sub.T] > 1.

3.7. Global Stability of Endemic Equilibrium (EE). To show the global stability of the EE, we use Lyapunov's theorem together with the following lemma.

Lemma 1. Suppose that [y.sub.1],[y.sub.2],...,[y.sub.n] are n positive numbers. Then, their arithmetic mean is greater than or equal to the geometric mean, that is ([y.sub.1] + [y.sub.2] + ... + [y.sub.n]/n) [greater than or equal to] [([y.sub.1][y.sub.2] ... [y.sub.n]).sup.1/n].

Theorem 5. The EE defined by E* is globally asymptotically stable if [R.sub.0] > 1, otherwise unstable.

Proof. The proof is based on the idea as explained by Martcheva [5]. We define Lyapunov's function as

[mathematical expression not reproducible] (33)

satisfying system (6) with [k.sub.1],[k.sub.2],[k.sub.3],[k.sub.4],[k.sub.5],[k.sub.6] > 0 to be determined. The function V is nonnegative for all ([s.sub.h],[e.sub.h],[i.sub.h],[r.sub.h],[s.sub.v],[i.sub.v]) [not equal to] ([s*.sub.h],[e*.sub.h],[i*.sub.h],[r*.sub.h],[s*.sub.v],[e*.sub.v],[i*.sub.v]) and radially unbounded.

We need to prove that V' < 0 for all ([s.sub.h],[e.sub.h],[i.sub.h],[r.sub.h],[s.sub.v],[e.sub.v],[i.sub.v]) [not equal to] ([s*.sub.h],[e*.sub.h],[i*.sub.h],[r*.sub.h],[s*.sub.v],[e*.sub.v],[i*.sub.v]). We find the derivative of V with respect to time and replace the derivatives [s'.sub.h], [e'.sub.h], [i'.sub.h], [r'.sub.h], [s'.sub.v], [e'.sub.v], [i'.sub.v]). with system (6). We also ignore the Th terms because if [s.sub.h],[e.sub.h],[i.sub.h] are globally stable, then the [r.sub.h] [right arrow] 0 at any time t and EE is globally stable:

[mathematical expression not reproducible] (34)

We now substitute [[pi].sub.h] = a[P.sub.h][i*.sub.v][s*.sub.h] + [[mu].sub.h][s*.sub.h] and [[pi].sub.v] = a[p.sub.v][i*.sub.h][s*.sub.h] + [[mu].sub.v][s*.sub.v] at the endemic equilibrium and then simplify and put similar terms together to obtain

[mathematical expression not reproducible] (35)

We suppose [k.sub.1] = [k.sub.2] and [k.sub.4] = [k.sub.5] and multiply and divide the same equilibrium value to some of the fractions to obtain

[mathematical expression not reproducible] (36)

We choose [k.sub.3] = [k.sub.2]([[mu].sub.h] + [epsilon]/[epsilon]) such that [k.sub.3]([omega] + [tau] + [[mu].sub.h])[i*.sub.h] = [k.sub.2]([epsilon] + [[mu].sub.h])[e*.sub.h] and choose [k.sub.6] = [k.sub.5]([[mu].sub.h] + [sigma]/[sigma]) such that [k.sub.6][p.sub.v][s*.sub.v][i*.sub.v] = [k.sub.5]([sigma] + [[mu].sub.v])[e*.sub.v]. Now, a[p.sub.h][s*.sub.h][i*.sub.v] = ([[mu].sub.h] + [epsilon])[e*.sub.h] and a[p.sub.v][s*.sub.v][i*.sub.v] = ([[mu].sub.v] + [sigma])[e*.sub.v] because [k.sub.1] = [k.sub.2] and [k.sub.4] = [k.sub.5], respectively. We now obtain

[mathematical expression not reproducible] (37)

Suppose [k.sub.6] = ([k.sub.1]a[p.sub.h][s*.sub.h]/[[mu].sub.v]) and [k.sub.3] = ([k.sub.4]a[p.sub.v][s*.sub.v]/[omega] + [tau] + [[mu].sub.h]). We then substitute and simplify to get

[mathematical expression not reproducible]. (38)

From Lemma 1, the terms

[mathematical expression not reproducible] (39)

Therefore, V'< 0 for all ([s.sub.h],[e.sub.h],[i.sub.h],[r.sub.h],[s.sub.v],[e.sub.v],[i.sub.v]) [not equal to] ([s*.sub.h], [e*.sub.h],[i*.sub.h],[r*.sub.h],[s*.sub.v],[e*.sub.v],[i.sub.v]), implying that the endemic equilibrium is globally asymptotically stable if [R.sub.0] > 1.

4. Analysis of Optimal Control Model

In this section, we formulate the optimal control model by modifying system (6) to an optimal control problem. We thus define some linear functions [c.sub.i](t) = 1, for i = 1,2,3. It is important to note that controls are fully effective when [c.sub.i](t) = 1 and not effective when [c.sub.i](t) = 0. The forces of infection [[xi].sub.h] and [xi.sub.v], which correspond to the human and vector population, respectively, are reduced by the factor (1 - [c.sub.1]), where [c.sub.1] measures the level of success obtained due to the effort of educating people on the dangers of exposing their skin, and encouraging them to wear long sleeves and long pants during the day to minimize tsetse fly-human contacts. The factor [c.sub.2] represents the effort of treatment to control the disease, and the factor [c.sub.3] also represents the effort of using insecticides to ensure that the breeding sites of the tsetse fly are minimized. Hence, taking into account the assumptions and extensions made, we try to find the most effective strategy that reduces the HAT infection in the population at a very minimum cost. With the use of bounded Lebesgue measurable control, we define the objective function to be minimized as

[mathematical expression not reproducible] (40)

Thus, the dynamics of the controls that minimizes the objective function is given by

[mathematical expression not reproducible] (41)

subject to the initial conditions [s.sub.h] [greater than or equal to] 0, [e.sub.h] [greater than or equal to] 0, [i.sub.h] [greater than or equal to] 0, [r.sub.h] [greater than or equal to] 0, [s.sub.v] [greater than or equal to] 0, [e.sub.v] [greater than or equal to] 0, and [i.sub.v] [greater than or equal to] 0. The associated effective reproduction number for (41) denoted by [R.sub.E] is obtained as

[R.sub.E] = [square root of [[a.sup.2][epsilon][(1- [c.sub.1]).sup.2][p.sub.h][[pi].sub.h][[pi].sub.v][sigma]]/[c.sub.3][[mu].sub.h]([epsilon]+[[mu].sub.h])[[mu].sup.2.sub .v]([[mu].sub.h]+[c.sub.2][tau]+[omega])] = [square root of [R.sub.c]]. (42)

The goal is to minimize the exposed and infectious human populations ([e.sub.h], [i.sub.h]), the exposed and infectious vector populations ([e.sub.v], [i.sub.v]), and the cost of implementing the control by the use of possible [c.sub.i], i = 1,2,3. The functional objective includes the social cost which relates to the resources that are needed for educating people on personal protection (1/2)[k.sub.1][c.sup.2.sub.1], the application of treatment (1/2)[k.sub.2][c.sup.2.sub.2], and spraying of tsetse fly operations (1/2)[k.sub.3][c.sup.2.sub.3]. The quantities [M.sub.1] and [M.sub.2], respectively, represent the associated cost with minimizing the exposed and infected human population, while [M.sub.3] and [M.sub.4] also represent the cost associated with minimizing the exposed and infected vectors, respectively. The quantity [t.sub.F] is the time period of intervention. As explained in [15], the costs corresponding to [M.sub.1][e.sub.h], [M.sub.2][i.sub.h], [M.sub.3][e.sub.v], and [M.sub.4][i.sub.v] are linear, while the cost control functions (1/2)[k.sub.1][c.sup.2.sub.1], (1/2)[k.sub.2][c.sup.2.sub.2], and (1/2)[k.sub.3][c.sup.2.sub.3] should be nonlinear and take a quadratic form. Therefore, we seek to minimize the objective function over the given time interval [0, [t.sub.F]]. Pontryagin's maximum principle is used to solve this optimal control problem and the derivation of the necessary conditions. The Lagrangian of the optimal control problem is given by

L = ([M.sub.1][e.sub.h] + [M.sub.2][i.sub.h] + [M.sub.3][e.sub.v] + [M.sub.4][i.sub.v] + 1/2[k.sub.1][c.sup.2.sub.1] + 1/2[k.sub.2][c.sup.2.sub.2] + 1/2[k.sub.3][c.sup.2.sub.3]). (43)

To determine the Lagrangian minimum value, we define the Hamiltonian, H, for the control problem as

[mathematical expression not reproducible] (44)

where [mathematical expression not reproducible] are adjoint variables or costate variables. The differential equations of adjoint variables are obtained by taking the partial derivatives of the Hamiltonian equation with respect to the state variables, which gives

[mathematical expression not reproducible] (45)

Theorem 6. Given the optimal controls [c*.sub.1],[c*.sub.2],[c*.sub.3] and the solutions [s.sub.h],[e.sub.h],[i.sub.h],[r.sub.h],[s.sub.v],[e.sub.v],[i.sub.v] of the corresponding state equations (41) and (40) which minimize J ([c.sub.1],[c.sub.2],[c.sub.3]) over the region [OMEGA], then there exist adjoint variables [mathematical expression not reproducible] satisfying

-d[[lambda].sub.i]/dt = [partial derivative]H/[[partial derivative].sub.i], i [member of] {[s.sub.h],[e.sub.h],[i.sub.h],[r.sub.h],[s.sub.v],[e.sub.v],[i.sub.v]}, (46)

and the optimal solution [c*.sub.1],[c*.sub.2],[c*.sub.3] is given by

[c*.sub.1] = min{1, max(0,[[??].sub.1])},

[c*.sub.2] = min{1, max(0,[[??].sub.2])},

[c*.sub.3] = min{1, max(0,[[??].sub.3])}. (47)

Proof. The Pontryagin's maximum principle described in [16, 17] is applied. Corollary 4.1 in [17] shows the existence of an optimal control due to the convexity of the integrand J with respect to [c.sub.1],[c.sub.1],[c.sub.3], and Lipschitz property of the state system with respect to the state variables. By using the optimal conditions

[partial derivative]H/[partial derivative][c.sub.1] = 0,

[partial derivative]H/[partial derivative][c.sub.2] = 0,

[partial derivative]H/[partial derivative][c.sub.3] = 0, (48)

we obtain,

[mathematical expression not reproducible] (49)

Solving (49), we have

[mathematical expression not reproducible] (50)

As stated earlier, the lower and upper boundaries for the control parameters are 0 and 1, respectively. If [[??].sub.1], [[??].sub.2], [[??].sub.3] < 1, then [c.sub.1] = [c.sub.2] = [c.sub.3] = 0 and if [[??].sub.1], [[??].sub.2], [[??].sub.3] > 1, then [c.sub.1] = [c.sub.2] = [c.sub.3] = 1, otherwise [c.sub.1] = [[??].sub.1], [c.sub.2] = [[??].sub.2], [c.sub.3] = [[??].sub.3]. Therefore, for the control parameters [c*.sub.1], [c*.sub.2], [c*.sub.3], we obtain the optimum value of the function J ([c.sub.1],[c.sub.2],[c.sub.3]).

4.1. Optimal Control Simulations. The Octave programming language is used to simulate the optimal control model using the set of parameters obtained from previously reported studies and datasets, which have been cited. Some of these parameters are assumed for the sake of illustrations. Table 2 represents the values of the model parameters used for simulations. The following initial conditions were considered:

[S.sub.h](0) = 30,

[e.sub.h](0) = 7,

[i.sub.h](0) = 2,

[r.sub.h](0) = 0,

[s.sub.v](0) = 40,

[e.sub.v](0) = 10,

[i.sub.v](0) = 3, (51)

and the weight constants were assumed to be

[M.sub.1] = 1,

[M.sub.2] = 2,

[M.sub.3] = 2,

[M.sub.4] = 2,

[K.sub.1] = 2,

[k.sub.2] = 10,

[k.sub.3] = 5. (52)

Figures 2 and 3 represent the control profiles at different values of [c.sub.1],[c.sub.2], and [c.sub.3], while the rest of the plots are the graphs of infectious human and vector population plotted against time in days and they represent the effect of optimal controls [c.sub.1], [c.sub.2], and [c.sub.3] in reducing the number of individuals infected. From Figure 4, we observe that the use of treatment and insecticides only has a significant impact in reducing the number of infectious individuals and they show that this strategy is effective to control tsetse flies and infected human populations. In Figure 5, we observe that the use of education and insecticides reduce the number of infectious individuals but the results depicted in Figure 5(a) shows that this strategy is not effective and efficient to control the infectious human population. From Figure 6, we observe that the use of education and treatment reduces the number of infectious individuals but the results from Figure 6(b) shows that this strategy is not effective and efficient to control the population of infectious tsetse flies. Lastly, the results depicted from Figure 7 shows that the strategy of using education, treatment, and insecticides is very efficient and effective to reduce the number of infected individuals. Therefore, the use of education, treatment, and insecticides simultaneously is very efficient and effective to eliminate HAT in Africa.

In epidemiology, a reproduction number less than unity implies that the disease can be eradicated in the long run.

Hence, choosing suitable parameters for the controls [c.sub.1],[c.sub.2], and [c.sub.3], it was observed that the effective reproduction numbers obtained for Figures 2-7 were 0.0075,0.00013,0.044, 0.000075, 0.0075, 0.00013, 0.044, and 0.000075, respectively. This shows that incorporating all the control measures, that is, educating individuals, giving treatment, and applying insecticides, is an effective method to help reduce the number secondary infections in the population which corresponds with eradicating the disease in the long run.

5. Conclusion

In this paper, we studied and analyzed the model for transmission of HAT, and determined the basic reproduction number. The local and global stabilities of disease-free equilibrium and endemic equilibrium were also proved. For the optimal control model, education, treatment, and insecticides as control measures were used to optimize the objective function defined by Equation (40). The numerical simulations of the optimal control model show that the best strategy to reduce the number of infected individuals is through the use of education, treatment, and insecticides. This is the effective and efficient method to eliminate the disease. Furthermore, the national authorities, nongovernmental organizations (NGOs), and stakeholders must not lose their interest in controlling the disease because neglecting this disease may cause rapid reoccurrence and much effect to the people who are at risk.

https://doi.org/10.1155/2018/5293568

Data Availability

The secondary data supporting this research are from previously reported studies and datasets, which have been cited. The processed data are available in cited references.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The authors thank the support from Next Einstein Initiative of the African Institute for Mathematical Sciences, Ghana (Annual Grant for Graduate Studies) and MasterCard Foundation.

References

[1] K. S. Rock, S. J. Torr, C. Lumbala, and M. J. Keeling, "Quantitative evaluation of the strategy to eliminate human african trypanosomiasis in the democratic republic of Congo," Parasites and Vectors, vol. 8, no. 1, p. 532, 2015.

[2] J. R. Franco, P. P. Simarro, A. Diarra, and J. G. Jannin, "Epidemiology of human african trypanosomiasis," Clinical Epidemiology, vol. 2014, no. 6, pp. 257-275, 2014.

[3] R. Brun, J. Blum, F. Chappuis, and C. Burri, "Human african trypanosomiasis," Lancet, vol. 375, no. 9709, pp. 148-159, 2010.

[4] World Health Organization et al., Control and Surveillance of Human African Trypanosomiasis: Report of a WHO Expert Committee, World Health Organization, Geneva, Switzerland, 2013.

[5] M. Martcheva, An Introduction to Mathematical Epidemiology, Springer, Vol. 61, Springer, Berlin, Germany, 2015.

[6] D. J. Rogers, "A general model for the african trypanosomiases," Parasitology, vol. 97, no. 1, pp. 193-212, 1988.

[7] J. W. Hargrove, R. Ouifki, D. Kajunguri, G. A. Vale, and S. J. Torr, "Modeling the control of trypanosomiasis using trypanocides or insecticide-treated livestock," PLoS Neglected Tropical Diseases, vol. 6, no. 5, p. 1615, 2012.

[8] D. Kajunguri, "Modelling the control of tsetse and african trypanosomiasis through application of insecticides on cattle in Southeastern Uganda," Ph.D. thesis, Stellenbosch University, Stellenbosch, South Africa, 2013.

[9] S. Olaniyi and O. S. Obabiyi, "Mathematical model for malaria transmission dynamics in human and mosquito populations with nonlinear forces of infection," International Journal of Pure and Applied Mathematics, vol. 88, no. 1, pp. 125-156, 2013.

[10] R. H. Chisholm, P. T. Campbell, Y. Wu, S. Y. C. Tong, J. McVernon, and N. Geard, "Implications of asymptomatic carriers for infectious disease transmission and control," Royal Society Open Science, vol. 5, no. 2, p. 172341, 2018.

[11] P. Yan and S. Liu, "Seir epidemic model with delay," ANZIAM Journal, vol. 48, no. 1, pp. 119-134, 2006.

[12] G. T. Azu-Tungmah, "A mathematical model to control the spread of malaria in Ghana," Ph.D. thesis, Kwame Nkrumah University of Science and Technology, Kumasi, Ghana, 2012.

[13] S. A. Pedro, S. Abelman, F. T. Ndjomatchoua, R. Sang, and H. E. Z. Tonnang, "Stability, bifurcation and chaos analysis of vector-borne disease model with application to rift valley fever," PloS One, vol. 9, no. 10, Article ID 108172, 2014.

[14] S. Olaniyi and O. S. Obabiyi, "Qualitative analysis of malaria dynamics with nonlinear incidence function," Applied Mathematical Sciences, vol. 8, no. 78, pp. 3889-3904, 2014.

[15] L. S. Pontryagin, V. G. Boltyanskii, and R. V. Gamkrelidze, Ef Mishchenko the Mathematical Theory of Optimal Processes, Interscience, New York, NY, USA, 1962.

[16] B. Heimann, "Fleming, WH/rishel, RW, deterministic and stochastic optimal control. New York-heidelbergberlin. Springer-verlag. 1975. xiii, 222 s, dm 60, 60," ZAMMJournal of Applied Mathematics and Mechanics, Zeitschrift fur Angewandte Mathematik und Mechanik, vol. 59, no. 9, p. 494, 1979.

[17] S. Davis, S. Aksoy, and A. Galvani, "A global sensitivity analysis for african sleeping sickness," Parasitology, vol. 138, no. 4, pp. 516-526, 2011.

[18] M. Artzrouni and J. P. Gouteux, "A compartmental model of sleeping sickness in central africa," Journal of Biological Systems, vol. 4, no. 4, pp. 459-477, 1996.

Hamenyimana Emanuel Gervas (iD),(1,2,3) Nicholas Kwasi-Do Ohene Opoku (iD),(1,4) and Shamsuddeen Ibrahim (iD)(1)

(1) African Institute for Mathematical Sciences, Biriwa, Cape Coast, Ghana

(2) University of Dar es Salaam, Dar es Salaam, Tanzania

(3) University of Dodoma, Dodoma, Tanzania

(4) University of Cape Coast, Cape Coast, Ghana

Correspondence should be addressed to Hamenyimana Emanuel Gervas; hamenyimana@aims.edu.gh

Received 10 July 2018; Revised 12 September 2018; Accepted 22 October 2018; Published 22 November 2018

Academic Editor: Konstantin Blyuss

Caption: Figure 1: Compartmental model for the transmission of human African trypanosomiasis.

Caption: Figure 2: (a) Control profile when [c.sub.1] = 0, [c.sub.2] [not equal to] 0, and [c.sub.3] [+ or -] 0 at [R.sub.c] = 0.0075. (b) Control profile when [c.sub.1] [not equal to] 0, [c.sub.2] = 0, and [c.sub.3] [not equal to] 0 at [R.sub.c] = 0.00013.

Caption: Figure 3: (a) Control profile when [c.sub.1] [not equal to] 0, [c.sub.2] [not equal to] 0, and [c.sub.3] = 0 at [R.sub.c] = 0.044. (b) Control profile when [c.sub.1] [not equal to] 0, [c.sub.2] [not equal to] 0, and [c.sub.3] [not equal to] 0 at [R.sub.c] = 0.000075.

Caption: Figure 4: Simulations of the model showing the efforts of treatment and insecticides only on infectious individuals at [R.sub.c] = 0.0075.

Caption: Figure 5: Simulations of the model showing the efforts of education and insecticides only on infectious individuals at [R.sub.c] = 0.00013.

Caption: Figure 6: Simulations of the model showing the efforts of education and treatment only on infectious individuals at [R.sub.c] = 0.044.

Caption: Figure 7: Simulations of the model showing the efforts of education, treatment, and insecticides on infectious individuals at [R.sub.c] = 0.000075.

Table 1: The description of model variables and parameters. Variable Description [s.sub.h] Susceptible human population [s.sub.v] Susceptible tsetse fly population [e.sub.h] and [e.sub.v] Exposed human and tsetse fly population, respectively [i.sub.h] and [i.sub.v] Infectious human and tsetse population, respectively [r.sub.h] Recovered human population Parameter Description [[pi].sub.h] Recruitment rate for human population [[pi].sub.v] Recruitment rate for tsetse fly population [P.sub.h] Proportion of bites by the infectious vector on susceptible human population [p.sub.v] Proportion of bites by susceptible vector on an infectious human population a The biting rate of the tsetse flies [sigma] Per capita rate of a vector becoming infectious [epsilon] Per capita rate of human becoming infectious [omega] Disease induced death rate The rate at which the recovered human can become [rho] susceptible again [tau] Recovery rate [mu] Natural death rate [[xi].sub.h] Force of infection for human population [[xi].sub.v] Force of infection for tsetse flies Table 2: Parameters values used for simulations. Parameter Value Reference [[pi].sub.h] 0.000215/day [14] [[pi].sub.v] 0.07/day [14] [p.sub.h] 0.62 [6] [P.sub.v] 0.065 [6, 15] a Varying Assumed [sigma] 0.001 Assumed [epsilon] 0.083 [18] [omega] 0.004 [3] [rho] 0.02 [6] [tau] 0.125 [3] [[mu].sub.h] 0.00044 Assumed [[mu].sub.v] 0.034 [15]

Printer friendly Cite/link Email Feedback | |

Title Annotation: | Research Article |
---|---|

Author: | Gervas, Hamenyimana Emanuel; Opoku, Nicholas Kwasi-Do Ohene; Ibrahim, Shamsuddeen |

Publication: | Computational and Mathematical Methods in Medicine |

Date: | Jan 1, 2018 |

Words: | 7102 |

Previous Article: | The Time-Domain Integration Method of Digital Subtraction Angiography Images. |

Next Article: | A New Approach towards Minimizing the Risk of Misdosing Warfarin Initiation Doses. |

Topics: |