Discrete Dynamics in Nature and Society
Volume 2017 (2017), Article ID 4057089, 7 pages
https://doi.org/10.1155/2017/4057089

Research Article

Numerical Analysis of Fractional Order Epidemic Model of Childhood Diseases

Fazal Haq,¹ Muhammad Shahzad,¹ Shakoor Muhammad,² Hafiz Abdul Wahab,¹ and Ghaus ur Rahman³

¹Department of Mathematics, Hazara University, Mansehra, Khyber Pakhtunkhwa, Pakistan
²Department of Mathematics, Abdul Wali Khan University, Mardan, Pakistan
³Department of Mathematics and Statistics, University of Swat, Khyber Pakhtunkhwa, Pakistan

Correspondence should be addressed to Fazal Haq; fazalhaqphd@gmail.com

Received 26 July 2017; Revised 4 October 2017; Accepted 5 November 2017; Published 3 December 2017

Academic Editor: Stefania Tomasiello

Copyright © 2017 Fazal Haq et al. 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.

Abstract.

The fractional order Susceptible-Infected-Recovered (SIR) epidemic model of childhood disease is considered. Laplace–Adomian Decomposition Method is used to compute an approximate solution of the system of nonlinear fractional differential equations. We obtain the solutions of fractional differential equations in the form of infinite series. The series solution of the proposed model converges rapidly to its exact value. The obtained results are compared with the classical case.

1. Introduction

Childhood diseases are most serious infectious diseases. Measles, poliomyelitis, and rubella are famous among them [1, 2]. Measles is a highly infectious disease, caused by respiratory infection by a Morbillivirus. These diseases normally affect the children, because child population is more prone to the disease as compared to the adults [1]. Therefore, the population can be divided into two major classes: premature and mature populations. Premature population takes a constant time to become mature, which is known as maturation delay. In disease dynamics, a disease cannot spread instantaneously but rather it will take some time in the body which is called latent period for the particular disease. For the control of childhood disease vaccination is a significant strategy being used all over the world. A universal effort to extend vaccination range to all children began in 1974, when the World Health Organization (WHO) founded the Expanded Program on Immunization [3]. Mathematical model plays an important role to comprehend the process of transmission of a disease and provides different techniques to control its propagation. Many mathematicians investigated childhood disease; for instance, Singh et al. [2] studied about vaccination of the childhood diseases. In addition, the authors presented numerical solution of the childhood disease model. Makinde [4] presented a Susceptible-Infected-Recovered model

The authors [1] rearranged model (1) using the relations , , and and obtained a new (SIR) model in the following form:

The description of the preceding model is given below.

The model shows that vaccination is 100 percent efficient and the natural death rate is unequal. Therefore the total population size is not constant. The birth rate is represented by while the rate of mortality of the childhood disease is very low. The parameter, , represents a fraction of vaccinated population at birth, where , considering that the rest of population is susceptible. A susceptible individual suffers from the disease through a contact with infected individuals at rate . Infected individuals recover at a rate . There are many applications of fractional calculus [5, 6]. Mathematicians and researchers used fractional calculus to model real world problems. Fractional calculus has greater degree of freedom; therefore it helps to solve nonlinear problems [7, 8]. Also with the help of fractional derivative interdisciplinary applications can be studied. The nonlinear oscillation of earthquake can be modeled with the help of fractional derivative.

In 1980, Laplace–Adomian Decomposition Method (LADM) was introduced by Adomian, which is an effective method for finding the numerical and explicit solution of a system of differential equations representing physical problems. This method works efficiently for solving several kinds of differential equations. It includes nonlinear boundary value/initial value problems. Moreover, we can use it for the partial differential equations; in addition, it can be used to solve a system of stochastic differential equations also. In this method, no perturbation or linearization is required. The advantage of LADM is that it needs no extra memory and does not require any additional parameters, which wastes time and memory.

Studying the literature review of childhood disease model it has been found that Arafa et al. [1] presented the fractional order model for the childhood disease and provided numerical solution by using homotopy analysis method. The reason behind the use of fractional order differential equations (FDEs) is that FDEs are naturally related to systems that involve memory, which can be found in many biological systems. Also, they show the realistic behavior of infection of disease but at a slower rate. Motivated by the applications of fractional calculus and LADM we explore numerical solution of childhood disease model. The Caputo derivative, which is a modification of the Riemann-Liouville definition [9, 10], is considered as a differential operator in our model.

We now gather some well-known definitions and results from the literature, which we will use throughout this paper. For more details, we refer the reader to [1, 7–12].

Definition 1. The Caputo fractional order derivative of a function on the interval is defined by

where and represents the integer part of .

Definition 2. We recall the definition of Laplace transform of Caputo derivative as

2. The Laplace–Adomian Decomposition Method

Using the Caputo fractional derivative system (2) gets the following form:

where while , , , and   are positive parameters and the given initial conditions are , , and . To solve system (5), we use Laplace–Adomian Decomposition Method (LADM) [5, 13]. Moreover, the obtained solution will be compared with the integer order derivative case. Furthermore, we use Laplace transformation to convert the system of differential equations into a system of algebraic equations. Then, the algebraic equations are used to obtain the required solution in form of series. We will discuss the procedure for solving model (5) with given initial conditions. Applying Laplace transform on both sides of model (5), we obtain the following system:

or

Using the initial conditions (7), we obtain the form

Assume that the solutions , , and in the form of infinite series are given by

while the nonlinear term is decomposed as follows:

where each is the Adomian polynomials defined as

The first three polynomials are given by

Substituting (9) and (10) into (8) results in

Matching the two sides of (13) yields the following iterative algorithm:

Taking Laplace inverse of (14) and considering first three terms, we get

3. Numerical Results and Discussion

Using , , , , , , , and , the LADM provides us with an approximate solution in the form of infinite series. Thus we calculate the first four terms of (5) to obtain

For , (16) attains the form

Similarly, we get the following system for :

Now, for , one can obtain the following system:

Similarly, the solution after three terms for is calculated as follows:

From the graphical results it is clear that the result obtained by using LADM is very efficient. It also shows that the presented method can predict the behavior of the variables accurately for the region under consideration. It is also clear that the efficiency of this method can be dramatically increased by increasing the terms. Fractional order derivative provides a greater degree of freedom as compared to integer order derivative. The dynamics of various compartments have been shown in Figure 1, Figure 2, and Figure 3, respectively. In addition, we give a comparison of and LADM in Tables 1 and 2 for , which shows that both the methods agree for short interval of time. The proposed method is better than method as it needs no extra predefined parameter which controls the method.

Table 1: Numerical solution of the proposed model using LADM at classical order .

Time (week)

Table 2: Numerical solution of the proposed model using at classical order .

Time (week)

Figure 1: Plot of approximate solutions of susceptible class corresponding to different fractional values of for .

Figure 2: Plot of approximate solutions of infected class corresponding to different fractional values of for .

Figure 3: Plot of approximate solutions of removed class corresponding to different fractional values of for .

4. Convergence Analysis

Solution (16) is in the form of series, which converges uniformly to the exact solution. To check the convergence of series (16), we use techniques (see [14]). For sufficient conditions of convergence of this method, we give the following theorem by using [14].

Theorem 3 (see [13]). Let be the Banach space and be a contractive nonlinear operator such that, for all .
By using Banach contraction principles has a unique point such that , where . The series given in (16) can be written by applying Adomian Decomposition Method as follows:

and assume that , where ; then, we have

(1);(2)

Theorem 4 (see [15]). Let be an operator from a Hilbert space into and be the exact solution of (5), where . , which is obtained by (9), converges to when , and , .

5. Conclusion

In this paper, we have considered a fractional order childhood disease model. The concern model is investigated for numerical solution by using LADM. The LADM is effective tool to solve nonlinear models and is widely used in engineering and applied mathematics. Also we have provided convergence results for the aforesaid method. It is clear that this method provides good results. One can observe that the method provides excellent numerical solutions for nonlinear fractional order models as compared to other methods like homotopy analysis, homotopy perturbation method, and . Since these methods involve an extra parameter on which the solutions depend, therefore our proposed method needs no parameter and is easy to understand as well as to implement. For initial value problems, our proposed method is better than the PDQM and cubic-spline DQM methods. As PDQM, cubic-spline DQM methods are based on discretization technique which needs extra memory and time as compared to LADM. The method includes a pseudospectral method referred to as the quadrature discretization method (PQDM) and cubic-spline DQM, which are based on nonstandard polynomial basis sets. The aforesaid method is slowly converging as compared to LADM [16].

Conflicts of Interest

The authors declared that no conflicts of interest exist regarding this paper.

Authors’ Contributions

All authors have equal contributions and they read and approved the final version of the paper. In the review process Dr. Shakoor Muhammad has provided valuable comments and suggestions and finally he checked and approved the final version. The review process without him was not possible.

References

A. A. M. Arafa, S. Z. Rida, and M. Khalil, “Solutions of fractional order model of childhood diseases with constant vaccination strategy,” Mathematical Science Letter, vol. 1, no. 1, pp. 17–23, 2012.

H. Singh, J. Dhar, H. S. Bhatti, and S. Chandok, “An epidemic model of childhood disease dynamics with maturation delay and latent period of infection,” Modeling Earth Systems and Environment, vol. 2, no. 79, 2016.

R. H. Henderson, “The expanded programme on immunization of the world health organization,” Reviews of Infectious Diseases, vol. 6, no. 2, pp. 475–479, 1984.

O. D. Makinde, “Adomian decomposition approach to a {SIR} epidemic model with constant vaccination strategy,” Applied Mathematics and Computation, vol. 184, no. 2, pp. 842–848, 2007.

C. F. Li, X. N. Luo, Y. Zhou, and J. J, “Existence of, positive solution of the boundary value problem for nonlinear fractinal differentialn equations,” Computer and Mathematics with Applications, vol. 59, pp. 1363–1375, 2010.

F. Haq, K. Shah, G.-U. Rahman, and M. Shahzad, “Numerical analysis of fractional order model of HIV-1 infection of CD4+ T-cells,” Computational Methods for Differential Equations, vol. 5, no. 1, pp. 1–11, 2017.

A. Babakhani and V. Daftardar-Gejji, “Existence of positive solutions of nonlinear fractional differential equations,” Journal of Mathematical Analysis and Applications, vol. 278, no. 2, pp. 434–442, 2003.

M. Benchohra, S. Hamani, and S. K. Ntouyas, “Boundary value problems for differential equations with fractional order and nonlocal conditions,” Nonlinear Analysis: Theory, Methods & Applications, vol. 71, no. 7-8, pp. 2391–2396, 2009.

D. Kumar and J. Singh, “Sushela, Application of homotopy analysis transform method to fractional biological population model,” Romanian Reports Physics, vol. 65, no. 1, pp. 63–75, 2013.

F. Haq, K. Shah, G. Rahman, and M. Shahzad, “Existence and uniqueness of positive solution to boundary value problem of fractional differential equations,” Sindh Univ. Res. Jour, vol. 48, no. 2, pp. 451–456, 2016.

A. Yildirim and Y. Cherruault, “Analytical approximate solution of a SIR epidemic model with constant vaccination strategy by homotopy perturbation method,” Kybernetes, vol. 38, no. 9, pp. 1566–1575, 2009.

Y. Lia, F. Haq, K. Shah, M. Shahzad, and G. Rahman, “Numerical analysis of fractional order Pine wilt disease model with bilinear incident rate,” Journal of Mathematics and Computer Science, vol. 17, pp. 420–428, 2017.

K. Shah, H. Khalil, and R. A. Khan, “Analytical solutions of fractional order diffusion equations by natural transform method,” Iranian Journal of Science & Technology, vol. 2016, pp. 1–13, 2016.

K. Shah, H. Khalil, and R. A. Khan, “Analytical solutions of fractional order diffusion equations by natural transform method,” Iranian Journal of Science & Technology, vol. 2016, 14 pages, 2016.

M. M. Hosseini and H. Nasabzadeh, “On the convergence of Adomian decomposition method,” Applied Mathematics and Computation, vol. 182, no. 1, pp. 536–543, 2006.

K. Leung, B. D. Shizgal, and H. Chen, “The quadrature discretization method ({QDM}) in comparison with other numerical methods of solution of the Fokker-Planck equation for electron thermalization,” Journal of Mathematical Chemistry, vol. 24, no. 4, pp. 291–319, 1998.