Research Article | Open Access

Aman Ullah, Thabet Abdeljawad, Shabir Ahmad, Kamal Shah, "Study of a Fractional-Order Epidemic Model of Childhood Diseases", *Journal of Function Spaces*, vol. 2020, Article ID 5895310, 8 pages, 2020. https://doi.org/10.1155/2020/5895310

# Study of a Fractional-Order Epidemic Model of Childhood Diseases

**Guest Editor:**Erdal Karapinar

#### Abstract

In this article, we discuss the existence and uniqueness of the solution of the fractional-order epidemic model of childhood diseases by using fixed point theory. The technique of natural transform coupled with the Adomian decomposition is used to find the solution of the proposed model. At the end of the article, the model is demonstrated with appropriate numerical and graphical description.

#### 1. Introduction

The immunity of the children is very weak; therefore, all the children are exposed to various infectious and noninfectious diseases during childhood [1]. The most common childhood diseases are either viral or bacterial infections or allergic and immunologic diseases [2]. “Respiratory Syncytial Virus” (RSV) causes bronchiolitis, which affects baby children up to 1 year. “Coxsackievirus” causes hand, foot, and mouth disease, which has no proper treatment for the infection and needs proper care. The infection due to “rotavirus” causes vomiting, high fever, and diarrhea in children, often leading to critical problems with dehydration. The “varicella” virus causes “chickenpox.” The “rubeola virus” causes measles, also known as “German measles,” leading to serious complications, even death. Mumps is another viral disease which initially starts with flu and then causes parotitis. Viruses or bacteria cause meningitis, which have symptoms like headache, stiff neck, fever, and malaise. It is an inflammation of the tissue surrounding the brain and spinal cord “meninges.” A strep infection causes “scarlet fever,” which commonly appears after the throat infection. There are many other diseases that are found in children in which it is not possible to list them all. The term modeling is derived from the Latin word “modulus” [3] which means that it is without restraint an available application that permits using mathematics to create model interactivity. We use the concept of mathematical modeling [4], for this divided the population into two major classes. The first one is a mature class while the other is a premature class. The premature class is then converted to a mature class with a constant time, known as maturation delay. To analyze the dynamics of childhood diseases, many of the diseases cannot be spread rapidly in the human body; for a certain disease, it will take a particular time called “latent period.” All over the world, proper vaccination is suggested to control all the mentioned childhood diseases. For this purpose, WHO launched the immunization program all over the world in 1974 [5]. Singh et al. [2] provided a mathematical model to formulate the procedure and present a numerical solution of vaccination in various diseases during childhood. The SIR model is presented by Makinde [6], as given below:

Arafa et al. rearranged the model (1) using the following relations, , and modified it as follows:

The model shows 100 percent efficiency of the vaccination. The natural death rate and the birth rate are denoted by and ; however, the rate of mortality of the childhood disease is very low. A fraction of the vaccinated population at birth is represented by the parameter , where . Assume that the rest of the population is susceptible. A susceptible individual suffers from the disease through a contact with infected individuals at rate .

In the last few decades, the research in applied mathematics particularly in modern calculus and fractional calculus has attracted scientists and researchers [7–10]. It should be kept in mind that in classical calculus, DEs have been extended to fill the research gap in many branches of pure and applied mathematics. Therefore, classical calculus is extended to modern calculus, and ODEs are extended to fractional-order DEs (FODEs). Fractional calculus has many applications in physical and natural sciences [11–18]. Different real-world problems have been modeled by FODEs and the system of FODEs [19–21]. Due to the degree of freedom, fractional calculus is used to solve nonlinear problems [22–28]. Various researchers and scientists discussed the existence and uniqueness of the solutions of different phenomena [12, 23, 29, 30] and applied different techniques to solve these problems [17, 20, 21, 31]. Haq et al. presented model (2) in a fractional order as

It is very essential to find the solution of (3) with an efficient technique. Many researchers suggested different analytical and numerical methods for finding the solution of nonlinear problems. In this paper, we will use natural transform coupled with the Adomian decomposition method (NTADM) for finding an approximate solution of a nonlinear problem. This method can also be used for solving PDEs in which no perturbation or linearization is required. This method does not need any extra memory or any additional parameters.

#### 2. Preliminaries

*Definition 1 (see [32]). *Let be a continuous function on , a fractional integral of order in the Riemann-Liouville sense corresponding to defined by
where integral on right side is pointwise defined on .

*Definition 2 (see [32]). *Let be a continuous function on . The Caputo fractional derivative may be expressed as
where and represents the integer part of .

*Definition 3 (see [33]). *The natural transform (abbreviated as NT) of the function is presented as and with variables and is defined as
Now, the NT with the Heaviside function is defined as

*Definition 4 (see [34]). *The NT of the arbitrary derivative in Caputo’s sense of is defined as

#### 3. Existence and Uniqueness

In this section, to address the existence and uniqueness of the proposed model, fixed-point theory is used. The existence of the model is guaranteed by Schauder’s fixed-point theorem, while the uniqueness is guaranteed by Banach’s contraction theorem.

System (2) can be generalized by involving a fractional derivative in the Caputo sense for as follows: with initial conditions , where are positive parameters. Assume that

System (9) becomes

Applying a fractional integral and using initial conditions, we have

Let

So, system (12) becomes

Consider a Banach space , with a norm

Let be a mapping defined as

Further, we impose some conditions on a nonlinear function as follows:

(C1). There exist constants and such that

(C2). There exists constant such that for each , such that

Theorem 5. *Suppose that condition (C1) holds. Then, the system (11) has at least one solution.*

*Proof. *To show that operator is bounded, let , where
is a closed and convex subset of . Now, take
It means that , which shows that is bounded. Next, to show that is completely continuous, let and take
This shows that as . Hence, operator is completely continuous by the Arzelá-Ascoli theorem. Thus, the given system (11) has at least one solution by Schauder’s fixed-point theorem.

Next, the Banach fixed-point theorem was used to show that the system (11) has unique solution.

Theorem 6. *Suppose that condition (C2) holds. Then, the system (11) has unique solution.*

*Proof. *Let . Take
Hence, is the contraction. Using the Banach contraction theorem, the system (11) has unique solution.

#### 4. Approximate Solution

In this part, the natural transform method is used to find the approximate solution of the system (9). By using NT and initial conditions, we get

Applying inverse NT, we get

Consider the series solution of the system (9) to be while the nonlinear term is represented as where each is the Adomian polynomial and is defined as

The first three terms of are given below:

Therefore, equation (24) can be written as

Thus, we establish the series solution as where , , and .

Similarly, we can calculate the other terms.

#### 5. Numerical Results

Using , , , , and , corresponding to these numerical values, we present the approximate solution (30) through MATLAB for ten terms in Figures 1–3.

In Figures 1–3, we have presented the approximate solutions of different compartments against the given data and corresponding to different fractional orders. Further, the obtained results show that the solution is continuously dependent on the time-fractional derivative and values of the parameters. Since the population of the susceptible class decreases, however, that of the recovered class gradually increases with time due to vaccination. As shown in Figure 1, the susceptibility is decreasing for some initial time which results in infection decrease as in shown Figure 2. This is because of vaccination which reduces infection, and hence, the recovery from disease is increasing as shown in Figure 3. The increase and decrease of different compartments are different corresponding to the fractional-order derivative. The concerned process for some initial time is faster on a smaller fractional order, and then, the process reverses and becomes slow at the same smaller fractional order. Hence, the fractional calculus presents global dynamics of the decay and the growth process in the biological models. Here, it is interesting that the number of peaks and their respective amplitudes are similar between models; however, there are variations in the timing of these peaks. The transient oscillations of the considered model under fractional orders are more stretched out than those of its corresponding integer-order model whose solutions experience longer interepidemic times. Hence, we concluded that depending on the particular values of the parameters, the fractional-order model either converges to the classical model and fits data, similarly, or fits the data better and outperforms the classical model correspondingly.

#### 6. Conclusion

In this paper, we show the existence and uniqueness of the fractional-order childhood disease model applying metric fixed-point theory. We also develop a proper procedure of natural transform, coupled with the Adomian decomposition method, and obtained an approximate solution of the proposed model. In the end, some numeric interpretation and their graphical presentation have been provided.

#### Data Availability

No data were used to support this study.

#### Conflicts of Interest

There does not exist any conflict of interest regarding this paper.

#### Authors’ Contributions

All authors have equal contribution in this research paper; particularly, the second author has significant contributions to revising the manuscript. He also polished the paper and its language style.

#### Acknowledgments

The second author would like to thank Prince Sultan University for funding this work through the research group Nonlinear Analysis Methods in Applied Mathematics (NAMAM) (Group Number RG-DES-2017-01-17).

#### References

- A. A. M. Arafa, S. Z. Rida, and M. Khalil, “Solutions of fractional order model of childhood diseases with constant vaccination strategy,”
*Mathematical Sciences Letters*, vol. 1, no. 1, pp. 17–23, 2012. View at: Publisher Site | Google Scholar - 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. 2, 2016. View at: Publisher Site | Google Scholar - E. Leah,
*Mathematical Models in Biology*, Society for Industrial and Applied Mathematics, 2005. View at: Publisher Site - N. F. Britton,
*Essential Mathematical Biology, Springer Undergraduate Mathematics Series*, Springer Verlag, UK, 2003. View at: Publisher Site - R. H. Henderson, “The expanded programme on immunization of the world health organization,”
*Clinical Infectious Diseases*, vol. 6, Supplement 2, pp. S475–S479, 1984. View at: Publisher Site | Google Scholar - 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. View at: Publisher Site | Google Scholar - K. S. Miller and B. Ross,
*An Introduction to the Fractional Calculus and Fractional Differential Equations, a Wiley-Interscience Publication*, JohnWiley & Sons, New York, NY, USA, 1993. - K. B. Oldham and J. Spanier,
*The Fractional Calculus*, Academic Press, New York, NY, USA, 1974. - K. Shah, T. Abdeljawad, I. Mahariq, and F. Jarad, “Qualitative analysis of a mathematical model in the time of COVID-19,”
*BioMed Research International*, vol. 2020, Article ID 5098598, 11 pages, 2020. View at: Publisher Site | Google Scholar - V. Lakshmikantham and A. S. Vatsala, “Basic theory of fractional differential equations,”
*Nonlinear Analysis: Theory, Methods & Applications*, vol. 69, no. 8, pp. 2677–2682, 2008. View at: Publisher Site | Google Scholar - D. Baleanu, J. T. Machado, and A. C. Luo,
*Fractional Dynamics and Control*, Springer Science and Business Media, 2011. - N. Ali and R. A. Khan, “Existence of positive solution to a class of fractional differential equations with three point boundary conditions,”
*Mathematical Sciences Letters*, vol. 5, no. 3, pp. 291–296, 2016. View at: Publisher Site | Google Scholar - R. S. Adigüzel, Ü. Aksoy, E. Karapinar, and İ. M. Erhan, “On the solution of a boundary value problem associated with a fractional differential equation,”
*Mathematical Methods in the Applied Sciences*, pp. 1–12, 2020. View at: Publisher Site | Google Scholar - E. Karapinar, T. Abdeljawad, and F. Jarad, “Applying new fixed point theorems on fractional and ordinary differential equations,”
*Advances in Difference Equations*, vol. 2019, no. 1, 2019. View at: Publisher Site | Google Scholar - B. Alqahtani, H. Aydi, E. Karapınar, and V. Rakočević, “A solution for Volterra fractional integral equations by hybrid contractions,”
*Mathematics*, vol. 7, no. 8, p. 694, 2019. View at: Publisher Site | Google Scholar - T. Abdeljawad, R. P. Agarwal, E. Karapınar, and P. S. Kumari, “Solutions of the nonlinear integral equation and fractional differential equation using the technique of a fixed point with a numerical experiment in extended b-metric space,”
*Symmetry*, vol. 11, no. 5, p. 686, 2019. View at: Publisher Site | Google Scholar - K. Shah, H. Khalil, and R. A. Khan, “Analytical solutions of fractional order diffusion equations by natural transform method,”
*Iranian Journal of Science and Technology Transactions A: Science*, vol. 42, no. 3, pp. 1479–1490, 2018. View at: Publisher Site | Google Scholar - J. Danane, K. Allali, and Z. Hammouch, “Mathematical analysis of a fractional differential model of HBV infection with antibody immune response,”
*Chaos, Solitons and Fractals*, vol. 136, article 109787, 2020. View at: Publisher Site | Google Scholar - 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. View at: Publisher Site | Google Scholar - D. Kumar, J. Singh, and Sushela, “Application of homotopy analysis transform method to fractional biological population model,”
*Romanian Reports Physics*, vol. 65, no. 1, pp. 63–75, 2013. View at: Google Scholar - F. Haq, K. Shah, G. 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. View at: Google Scholar - C. F. Li, X. N. Luo, and Y. Zhou, “Existence of positive solutions of the boundary value problem for nonlinear fractional differential equations,”
*Computers & Mathematics with Applications*, vol. 59, no. 3, pp. 1363–1375, 2010. View at: Publisher Site | Google Scholar - 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. View at: Publisher Site | Google Scholar - M. Belmekki, J. J. Nieto, and R. Rodríguez-López, “Existence of periodic solution for a nonlinear fractional differential equation,”
*Boundary Value Problems*, vol. 2009, Article ID 324561, 18 pages, 2009. View at: Publisher Site | Google Scholar - M. Benchohra, A. Cabada, and D. Seba, “An existence result for nonlinear fractional differential equations on Banach spaces,”
*Boundary Value Problems*, vol. 2009, 11 pages, 2009. View at: Publisher Site | Google Scholar - A. Atangana and Z. Hammouch, “Fractional calculus with power law: the cradle of our ancestors,”
*The European Physical Journal Plus*, vol. 134, no. 9, 2019. View at: Publisher Site | Google Scholar - M. K. Owolabi and Z. Hammouch, “Mathematical modeling and analysis of two-variable system with noninteger-order derivative,”
*Chaos: An Interdisciplinary Journal of Nonlinear Science*, vol. 29, no. 1, article 013145, 2019. View at: Publisher Site | Google Scholar - M. K. Owolabi and Z. Hammouch, “Spatiotemporal patterns in the Belousov-Zhabotinskii reaction systems with Atangana-Baleanu fractional order derivative,”
*Physica A: Statistical Mechanics and its Applications*, vol. 523, pp. 1072–1090, 2019. View at: Publisher Site | Google Scholar - R. A. Khan and K. Shah, “Existence and uniqueness of solutions to fractional order multi-point boundary value problems,”
*Communications in Applied Analysis*, vol. 19, pp. 515–526, 2015. View at: Google Scholar - F. Haq, K. Shah, G. Rahman, and M. Shahzad, “Existence and uniqueness of positive solution to boundary value problem of fractional differential equations,”
*Sindh University Research Journal*, vol. 48, no. 2, pp. 451–456, 2016. View at: Google Scholar - A. Arara, M. Benchohra, N. Hamidi, and J. J. Nieto, “Fractional order differential equations on an unbounded domain,”
*Nonlinear Analysis: Theory, Methods & Applications*, vol. 72, no. 2, pp. 580–586, 2010. View at: Publisher Site | Google Scholar - I. Podlubny,
*Fractional Differential Equations, Mathematics in Science and Engineering*, Academic Press, New York, NY, USA, 1999. - Z. H. Khan and W. A. Khan, “N-transform properties and applications,”
*NUST Journal of Engineering Sciences*, vol. 1, no. 1, pp. 127–133, 2008. View at: Google Scholar - D. Loonker and P. K. Banerji, “Solution of fractional ordinary differential equations by natural transform,”
*International Journal of Mathematics and Mathematical Sciences*, vol. 12, no. 2, 7 pages, 2013. View at: Google Scholar

#### Copyright

Copyright © 2020 Aman Ullah 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.