Abstract

We examine possible approximate solutions of both integer and noninteger systems of nonlinear differential equations describing tuberculosis disease population dynamics. The approximate solutions are obtained via the relatively new analytical technique, the homotopy decomposition method (HDM). The technique is described and illustrated with numerical example. The numerical simulations show that the approximate solutions are continuous functions of the noninteger-order derivative. The technique used for solving these problems is friendly, very easy, and less time consuming.

1. Introduction

Tuberculosis, MTB, or TB (short for tubercle bacillus) is a common and in many cases lethal, infectious disease caused by various strains of Mycobacterium, usually Mycobacterium Tuberculosis [1]. Tuberculosis typically attacks the lungs, but can also affect other parts of the body. It is spread through the air when people who have an active TB infection cough, sneeze, or otherwise transmit their saliva through the air [2]. Most infections are asymptomatic and latent, but about one in ten latent infections eventually progresses to active disease which, if left untreated, kills more than 50% of those so infected. Interested reader can find more about this model in [37].

Based on the standard SIRS model, the model population was compartmentalised into the susceptible and the infected which is further broken down into latently infected and actively infected while the recovered subpopulation is ploughed back into the susceptible group due to the possibility of reinfection after successful treatment of the earlier infection. The model monitors the temporary dynamics in the population of susceptible people , TB latently infected people , and TB actively infected people as captured in the model system of ordinary differential equations that follows. subject to the initial conditions where is the total number of new people in the location of interest; is the number of susceptible people in the location; is the number of TB latently infected people; is the number of TB actively infected people; is the probability that a susceptible person is not vaccinated; is the efficient rate of vaccines; is the success rate of latent therapy; is the active TB treatment cure rate; is the TB instantaneous incidence rate per susceptible; is humans natural death rate; is the proportion of infection instantaneously degenerating into active TB; is the TB-induced death rate; and is the breakdown rate from latent to active TB. The equilibrium analysis of the model was studied in [8]. Equation (1) together with (2) does not have an exact solution and is usually solved numerically.

The purpose of this paper is to derive approximate analytical solutions for the standard form as well as the fractional version of (1) together with (2) using the relatively new analytical technique, the homotopy decomposition method (HDM).

The paper is structured as follows. In Section 2, we present the basic ideal of the homotopy decomposition method for solving partial differential equations. We present the application of the HDM for system Tuberculosis disease population dynamics model in Section 3. In Section 4, we present the application of the HDM for system of fractional Tuberculosis disease population dynamics model. The conclusions are then given finally in Section 5.

2. Fundamental Information about Homotopy Decomposition Method

To demonstrate the elementary notion of this technique, we consider a universal nonlinear nonhomogeneous partial differential equation with initial conditions of the following form [913]. focused on the primary condition where is the order of the derivative, where is an identified function, is the common nonlinear differential operator, denotes a linear differential operator, and is the order of the derivative. The procedures first stage here is to apply the inverse operator on both sides of (3) to obtain The multi-integral in (3) can be transformed to so that (3) can be reformulated as Using the homotopy scheme, the solution of the aforementioned integral equation is given in series form as and the nonlinear term can be decomposed as where is an implanting parameter. is the polynomials that can be engendered by The homotopy decomposition method is obtained by the combination of decomposition method with Abel integral and is given by with Relating the terms of same powers of , this gives solutions of various orders. The initial guess of the approximation is that is actually the Taylor series of the exact solution of order . Note that this initial guess insures the uniqueness of the series decompositions [9].

3. Application of the HDM to the Model with Integer-Order Derivative

In this section, we employ this method for deriving the set of the mathematical equations describing the tuberculosis disease population dynamics model.

Resulting from the steps involved in the HDM method, we reach at the following integral equations that are very simple to solve: Integrating the previous, we obtain the following components: For simplicity, let us put In general, we obtain the following recursive formulas: where , , and depend on the fixed set of empirical parameters. It therefore follows that the approximate solution of the system (1) is given as If for instance one supposes that the total number of new people in the location of interest is ; the initial number of susceptible people in the location is ; the initial number of TB latently infected people is ; the initial number of TB actively infected people is ; the probability that a susceptible person is not vaccinated is ; the efficient rate of vaccines is ; the success rate of latent TB therapy is ; the active TB treatment cure rate is ; the TB instantaneous incidence rate per susceptible is ; humans natural death rate is ; the proportion of infection instantaneously degenerating into active TB is ; the TB-induced death rate is ; and the breakdown rate from latent to active TB is , then the following approximate solution is obtained as a result of the first 8 terms of the series decomposition: If in addition we assume that no new person migrates or is born in this area, we obtain the following figures. The approximate solutions of the main problem are depicted in Figures 1, 2, and 3, respectively.

Figure 1 shows that, if there is migration or newborn in the location of interest, the number of susceptible people will vanish as time goes, because of the natural death rate and due to TB. Note that any person that is latently infected is removed from the set of susceptible. Figure 2 indicates that the number of people that are latently infected will increase up to a certain time and then vanish as time goes. The number of susceptible people, will become latently infected since some are not vaccinated against the TB and finally will vanish due to. Figure 3 indicates that the number of TB actively infected people will also vanish because of the natural death rate and the death due to TB.

4. Application of the HDM to the Model with Noninteger-Order Derivative

Fractional calculus has been used to model physical and engineering processes, which are found to be best described by fractional differential equations. It is worth noting that the standard mathematical models of integer-order derivatives, including nonlinear models, do not work adequately in many cases. In the recent years, fractional calculus has played a very important role in various fields such as mechanics, electricity, chemistry, biology, economics, notably control theory, and signal and image processing. Major topics include anomalous diffusion; vibration and control; continuous time random walk; Levy statistics, fractional Brownian motion; fractional neutron point kinetic model; power law; Riesz potential; fractional derivative and fractals; computational fractional derivative equations; nonlocal phenomena; history-dependent process; porous media; fractional filters; biomedical engineering; fractional phase-locked loops, and groundwater problem (see [1421]).

4.1. Properties and Definitions

Definition 1. A real function , , is said to be in the space , , if there exists a real number , such that , where , and it is said to be in space if , .

Definition 2. The Riemann-Liouville fractional integral operator of order , of a function , , is defined as
Properties of the operator can be found in [1416]. We mention only the following: for , , , and ,

Lemma 3. If , , and , , then

Definition 4 (partial derivatives of fractional order). Assume now that is a function of variables also of class on . We define partial derivative of order for respect to the function where is the usual partial derivative of integer-order .

4.2. Approximate Solution of Fractional Version

The system of equations under investigation here is given as Following the discussion presented earlier, we arrive at the following equations: Integrating the previous, we obtain the following components: The remaining terms can be obtained in the same manner. But here only few terms of the series solutions are considered, and the asymptotic solution is given as The following figures show the simulated solutions for different values of the fractional order derivatives. The approximate solutions of the main problem are depicted in Figures 4, 5, 6, 7, 8, and 9, respectively.

The numerical simulations show that the approximate solutions are continuous functions of the noninteger-order derivative. It is worth noting that the standard mathematical models of integer-order derivatives, including nonlinear models, do not work adequately in many cases. It is therefore advisable to use the fractional model for describing this problem.

5. Conclusion

The tuberculosis model was examined for the case of integer- and noninteger-order derivatives. Both systems of nonlinear equations were solved with an iterative analytical model called the homotopy decomposition model method. The basic characters of the relatively new technique are presented in detail. The approximate solutions of the noninteger case are increasing continuous functions of the fractional order derivative. The technique used for solving these problems is friendly, very easy, and less time consuming. The numerical solutions in both cases display the biological behaviour of the real world situation.