Computational and Theoretical Analysis of Human Diseases Associated with Infectious Pathogens
View this Special IssueResearch Article  Open Access
On the Mathematical Analysis of Ebola Hemorrhagic Fever: Deathly Infection Disease in West African Countries
Abstract
For a given West African country, we constructed a model describing the spread of the deathly disease called Ebola hemorrhagic fever. The model was first constructed using the classical derivative and then converted to the generalized version using the betaderivative. We studied in detail the endemic equilibrium points and provided the Eigen values associated using the Jacobian method. We furthered our investigation by solving the model numerically using an iteration method. The simulations were done in terms of time and beta. The study showed that, for small portion of infected individuals, the whole country could die out in a very short period of time in case there is not good prevention.
1. Introduction
We shall recall that the filoviruses belong to a virus family named Filoviridae. This virus can cause unembellished hemorrhagic fever including humans and also monkeys [1–3]. In the literature, only two members of this virus family have been distinguished, namely, Marburgvirus and Ebolavirus. However, so far only in the literature five species of Ebolavirus have been identified including: Ivory Coast, Sudan, Zaire, Reston, and Bundibugyo. Amount these five species, Ebola virus is the only member of the Zaire ebolavirus species and the most dangerous, being responsible for largest number of outbreaks [1–4]. Ebola is an unusual nevertheless fatal virus that causes bleeding inside and outside the body [4]. As the virus spreads through the body, it damages the immune system and organs. Ultimately, it causes levels of bloodclotting cells to drop [5]. This leads to severe, uncontrollable bleeding [5].
One of the big challenges faced by West African countries is lack of employments. It is estimated that about 133 million young people (more than 50 percent of the youth population) in Africa are illiterate. Many young people have little or no skills and are therefore largely excluded from productive economic and social life. Those that have some education often exhibit skills irrelevant to current demand in the labour market, in a situation where educational and skill requirements are increasing, resulting in millions of unemployed and underemployed youth. The incidence of youth unemployment in SubSaharan African is estimated to be over 20 percent. Therefore, to provide for their family daily needs, many families in West and Central Africa reply on farming, fishing, and hunting. In the case of hunt, they will kill many wild animals, sometime dry and supply to the rural markets. Some images underpinning this situation are depicted in Figure 1. These pictures show a trade of these bust meats in some rural area in a SubSahara country. Figure 2 shows the life cycle of the ebolavirus and Table 1 shows parameters used for simulations according to some reported data.

Incidentally, EVD (Ebola Virus Disease) is believed to take place after an Ebola virus is transmitted to an initial human by contact with an infected animal’s body fluid. On the other hand, humantohuman transmission can take place with direct contact with blood or bodily fluid from an infected person. Fruit bats are considered to be the most likely natural source of the Ebola virus. As initial transmission, the bat drops incompletely eaten fruits and pulp and then lands animals for instance, gorilla feed fallen fruits, and then the hunter kills and sells the body of the affected animal. Now the affected human is in contact with the rest of his family. The chain of transmission can be viewed in the following picture.
Since all physical problems can be modeled via mathematical equation, we aim in this paper to analyze the spread of this deadly disease using mathematical equations. We shall propose a model underpinning the spread of this disease in a given SubSaharan African country in the next section.
2. Mathematical Formulation
Let us consider a West African country with a total number of populations at a given time. Let us assume that the rate of death caused by natural death and other diseases is factored out to be . , , , and are to be rate of infection by Ebola, rate of recovery, rate of susceptibility, and rate of death by Ebola, respectively. In this work, , , , and will describe the susceptible, the infected, the recovery, and the total death populations, respectively. Therefore, the mathematical equation describing the rate of change of susceptible population is given as follows:
The above equation is obtained because is the rate of infectious persons from recovery population turned out to be vulnerable again at the rate and finally the quantity of population that die with natural death and other diseases at the rate . The differential describing the rate of change of infected group is provided as
The above equation can be justified as follows: the total number of persons removed from susceptible group can be mathematically expressed as . However, due to the introduction of medication, a number of individuals will be recovered at a rate of and also a quantity of infected person will die at a rate . The rate of recovery population is described by the following ordinary differential equation:
The above is quite easy to be obtained. Finally the rate of change of death in that given country will be described via the ordinary equation
Thus, the set of mathematical equations involved can then be given as
Although the above equation is novel and valid, we have employed the common Newtonian idea of rate of change which has been intensively disparaged to not describe accurately the real rate of change in time [6–10]. Therefore, we shall consider modifying the above equation using a derivative that can take into account the time scale and also obeys the classical properties of the Newtonian derivative, for instance, chain rule and others. This derivative was recently proposed in the work [11] and defined as for all , . Then, if the limit of the above exists, is said to be differentiable.
Theorem 1 (see [11]). Assuming that is differential and differentiable the on the opened interval , then
Definition 2 (see [11]). Let be given function, and then we propose that the integral of order integral of is
The above operator is the inverse operator of the proposed beta derivative and is called the Atanganabeta integral. Now using the above derivative, we change (5) to
The above equation shall be referred to as Atangana’s Beta Ebola System of Equations (ABESE). We shall in the next section give a detailed analysis of the ABESE.
3. Analysis and Validation
3.1. Validation
The first aspect of this analysis is to verify that the such that . Indeed, by adding equations in system (9), we have
Now, using the linearity of the Betaderivative, we have
Now applying the inverse operator of given in Definition 2, we obtain
3.2. Endemic Equilibrium Points
We shall use only the first three equations of system (9) to find the endemic equilibrium point. Thanks to the Betaderivative that allows us to have that, a beta derivative of a constant is zero because the equilibrium points are obtained here by assuming that the solution of (9) does not depend on the time; then,
After some manipulations, we obtain the following solutions:
According to the survey done in the last months, it was revealed that the disease has a high rate of mortality with a maximum of 90 percent [12]. This implies the recovery rate is very small, and then the existence conditions are true and are conformable to the real world situation.
3.3. EigenValues Solutions
We shall make use of the Jacobian method to find the Eigenvalue associate of the endemic equilibrium points. The Jacobian matrices associate is given as
However, the free disease equilibrium point where , we have the following:
To find the Eigenvalue associate, we solve the following equation:
And the solutions of the above equation are given as
4. Derivation of the Solution
We shall in this section present the derivative of the solution for system (2). Since the system is nonlinear, some analytical techniques such as Laplace transform, Fourier transform, and Green function will not be suitable for this case. Suitable methods for nonlinear equations have been documented in the literature, for instance, homotopy perturbation method and its derivatives [13–15] and variational iteration method and its derivatives [16–18]. However, in this work, we shall use the homotopy decomposition method [19, 20].
In this method, we first apply the inverse operator of defined in (8) on system (5) to obtain
The next move in this method is to assume that, since the system is nonlinear, then the solutions can be obtained in series as
However, introducing the above proposed solution into (19), also making use of the idea of homotopy and after we compare terms of the same power of , we obtain the following:
For any , we have the following:
One of the important parts of any iteration method is to prove the uniqueness and the convergence of the method; we are going to show the analysis underpinning the convergence and the uniqueness of the proposed method for the general solution for .
Theorem 3. Assuming that and are Banach spaces and is contraction nonlinear mapping. If the progression engender by the threedimensional homotopy decomposition method is regarded as then, the following statements hold:(a), with ;(b)For any other greater than 0, is always in the neighborhood of the exact solution ;(c).
Proof (a). The proof of (a) shall be achieved via induction on the natural number . When , we have the following:
However, by hypothesis, we have that has a fixed point, which is the exact solution. Because
if is the exact solution, then,
since is the same as , therefore, we have that
Then,
Since is a contractive nonlinear mapping, we shall have the following inequality:
then, the property is verified for . Assume that the hypothesis is verified for ; we shall prove that it is also verified for . However, at the level , we have the following:
Using the fact that is a nonlinear contractive mapping, we have the following:
Furthermore, using the induction hypothesis, we arrive at
And the proof is completed. The proof of (b) and (c) can be found [19].
In order to show the efficiency and applicability of this method for handling the system nonlinear equations we shall present some numerical solution in Figures 3, 4, 5, 6, and 7.
The model depends on the order of the derivative; according to the prediction, when beta is 1, meaning the model with the classical derivative, the total number of individuals in the given West African country will all die, which is not realistic because we assume that, amount the 1000 persons living in a given West African country, only 900 are susceptible. We are expecting that, in worse case, all the 900 susceptible will die and only those who are not susceptible will live. Now, we observe from the figure that, for beta less than 0.5, the model predicts that, at the time , about 540 will be dead, about 300 will recover from the disease, about 40 will be infected, and only about 10 will be susceptible.
5. Conclusion
The derivative used to model a real world situation is very important. The classical derivative describes the change of rate, but it is an approximation of the real velocity of the object under study; this study has been examined in detail in [6]. The beta derivative is the modification of the classical derivative that takes into account the time scale and also has new parameter that can be considered as the fractional order. We have used the beta derivative to model the spread of the fatal disease called Ebola that has killed many people in the West African countries including Nigeria, Sierra Leone, Guinea, and Liberia since December 2013. We did the investigation of the stable endemic points and presented the Eigenvalues using the Jacobian method. The homotopy decomposition method was used to solve the resulted system of equations. The convergence of the method was presented and some numerical simulations were done for different values of beta. The simulations showed that our model is more realistic for all betas less than 0.5. The model revealed that, for a given population in a West African country, if there were no precaution of recovering, even if the total number of infected populations is very small, the entire population of that country would all die in a very short period of time. Based on the prediction of this paper, we are calling upon more research around this disease; in particular, we are calling for researchers to pay their attention to finding a very good cure or a better prevention, to reduce the risk of contamination.
Conflict of Interests
Both authors declare there is no conflict of interests for this paper.
Acknowledgments
The authors would like to thank the reviewers and Editor for their valuable suggestions. Abdon Atangana would like to thank Claude Leon Foundation for their financial support.
References
 E. Leroy, J. P. Gonzalez, and X. Pourrut, “Ebolavirus and other filoviruses,” Current Topics in Microbiology and Immunology, vol. 315, pp. 363–387, 2007. View at: Publisher Site  Google Scholar
 E. K. Leffel and D. S. Reed, “Marburg and Ebola viruses as aerosol threats,” Biosecurity and Bioterrorism : Biodefense Strategy, Practice, and Science, vol. 2, no. 3, pp. 186–191, 2004. View at: Publisher Site  Google Scholar
 J. H. Kuhn, S. Becker, H. Ebihara et al., “Proposal for a revised taxonomy of the family Filoviridae: classification, names of taxa and viruses, and virus abbreviations,” Archives of Virology, vol. 155, no. 12, pp. 2083–2103, 2010. View at: Publisher Site  Google Scholar
 T. Hoenen, A. Groseth, D. Falzarano, and H. Feldmann, “Ebola virus: unravelling pathogenesis to combat a deadly disease,” Trends in Molecular Medicine, vol. 12, no. 5, pp. 206–215, 2006. View at: Publisher Site  Google Scholar
 D. G. Bausch, H. Feldmann, T. W. Geisbert et al., “Outbreaks of filovirus hemorrhagic fever: time to refocus on the patient,” The Journal of Infectious Diseases, vol. 196, supplement 2, pp. S136–S141, 2007. View at: Publisher Site  Google Scholar
 I. Podlubny, “Geometric and physical interpretation of fractional integration and fractional differentiation,” Fractional Calculus & Applied Analysis, vol. 5, no. 4, pp. 367–386, 2002. View at: Google Scholar
 M. Caputo, “Linear model of dissipation whose Q is almost frequency independentII,” Geophysical Journal of the Royal Astronomical Society, vol. 13, no. 5, pp. 529–539, 1967. View at: Google Scholar
 A. A. Kilbas, H. M. Srivastava, and J. J. Trujillo, Theory and Applications of Fractional Differential Equations, vol. 204 of NorthHolland Mathematics Studies, Elsevier Science B.V., Amsterdam, The Netherlands, 2006. View at: MathSciNet
 A. Atangana and N. Bildik, “The use of fractional order derivative to predict the groundwater flow,” Mathematical Problems in Engineering, Article ID 543026, 9 pages, 2013. View at: Publisher Site  Google Scholar  MathSciNet
 Ş. Yüzbaşı, “A numerical approach to solve the model for HIV infection of CD4+T cells,” Applied Mathematical Modelling, vol. 36, no. 12, pp. 5876–5890, 2012. View at: Publisher Site  Google Scholar  MathSciNet
 A. Atangana and E. F. D. Goufo, “Extension of match asymptotic method to fractional boundary layers problems,” Mathematical Problems in Engineering, vol. 2014, Article ID 107535, 7 pages, 2014. View at: Publisher Site  Google Scholar
 C. M. Fauquet, Virus Taxonomy Classification and Nomenclature of Viruses; 8th Report of the International Committee on Taxonomy of Viruse, Elsevier/Academic Press, Oxford, UK, 2005.
 J. He, “Homotopy perturbation technique,” Computer Methods in Applied Mechanics and Engineering, vol. 178, no. 34, pp. 257–262, 1999. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 Y. Tan and S. Abbasbandy, “Homotopy analysis method for quadratic Riccati differential equation,” Communications in Nonlinear Science and Numerical Simulation, vol. 13, no. 3, pp. 539–546, 2008. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 M. A. Gondal and M. Khan, “Homotopy perturbation method for nonlinear exponential boundary layer equation using Laplace transformation, He’s polynomials and padé technology He’s polynomials and padé technology,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 11, no. 12, pp. 1145–1153, 2010. View at: Google Scholar
 Z. Odibat and S. Momani, “Numerical methods for nonlinear partial differential equations of fractional order,” Applied Mathematical Modelling, vol. 32, no. 1, pp. 28–39, 2008. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 G. C. Wu and D. Baleanu, “Variational iteration method for the Burgers' flow with fractional derivatives—new Lagrange multipliers,” Applied Mathematical Modelling, vol. 37, no. 9, pp. 6183–6190, 2013. View at: Publisher Site  Google Scholar  MathSciNet
 I. Andrianov and J. Awrejcewicz, “Construction of periodic solutions to partial differential equations with nonlinear boundary conditions,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 1, no. 4, pp. 327–332, 2000. View at: Publisher Site  Google Scholar  MathSciNet
 A. Atangana and E. Alabaraoye, “Solving system of fractional partial differential equations arisen in the model of HIV infection of CD4^{+} cells and attractor onedimensional KellerSegel equation,” Advances in Difference Equations, vol. 2013, article 94, 14 pages, 2013. View at: Google Scholar
 M. Y. Ongun, “The Laplace adomian decomposition method for solving a model for HIV infection of CD4^{+}T cells,” Mathematical and Computer Modelling, vol. 53, no. 56, pp. 597–603, 2011. View at: Publisher Site  Google Scholar  MathSciNet
Copyright
Copyright © 2014 Abdon Atangana and Emile Franc Doungmo Goufo. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.