Discrete Dynamics in Nature and Society
Volume 2020 (2020), Article ID 8709393, 11 pages
https://doi.org/10.1155/2020/8709393

Research Article

Existence Theory and Novel Iterative Method for Dynamical System of Infectious Diseases

Gauhar Ali,¹ Ghazala Nazir,¹ Kamal Shah,¹ and Yongjin Li²

¹Department of Mathematics, University of Malakand, Chakdara, Dir (Lower), Khyber Pakhtunkhwa, Pakistan
²Department of Mathematics, Sun Yat-sen University, Gangzhou, China

Correspondence should be addressed to Yongjin Li; stslyj@mail.sysu.edu.cn

Received 24 February 2020; Revised 30 April 2020; Accepted 13 May 2020; Published 1 June 2020

Academic Editor: Emilio Jiménez Macías

Copyright © 2020 Gauhar Ali 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.

This manuscript is devoted to investigate qualitative theory of existence and uniqueness of the solution to a dynamical system of an infectious disease known as measles. For the respective theory, we utilize fixed point theory to construct sufficient conditions for existence and uniqueness of the solution. Some results corresponding to Hyers–Ulam stability are also investigated. Furthermore, some semianalytical results are computed for the considered system by using integral transform due to the Laplace and decomposition technique of Adomian. The obtained results are presented by graphs also.

1. Introduction

In recent time, the subject of fractional calculus has got much attention from the researcher. This is due to large numbers of applications in various disciplines of science and engineering where the concept of derivatives and integrals is frequently used. Normally, the integer-order derivative does not explore the dynamics of real-world problems more comprehensively as compared to the fractional derivative. Also, the fractional differential operator is global and poses greater degree of freedom, while the ordinary differential operator is local and often cannot explain the memory and hereditary process of the real-world problem more efficiently. In fact, fractional derivative is definite integral over a domain; therefore, it has been defined in many ways. In this regard, a variety of definitions have been introduced by different researchers in the literature. Some famous definitions are those given by Riemann and Liouville, Caputo, Hadamard, and so on (see [1]). The definitions of Riemann–Liouville and Caputo have been very well used in applied problems. These definitions involve singular kernel which often causes difficulty in dealing some problems. Therefore, Caputo and Fabrizio in 2015 introduced a new concept about fractional-order derivatives based on nonsingular kernel (see [2–4]). Some remarkable merits of Caputo–Fabrizio fractional derivative and integral and their applications were given by many researchers in previous few years (see, for detail, [5–13]). The concerned fractional integral of a function is the fractional average of the function itself and its fractional integral in Riemann–Liouville sense. Moreover, in some articles, it has been displayed that the derivative has some constructive applications in thermal science, material sciences, and so on (see, for detail, [11, 12, 14, 15]). Since differential operators have greater degree of freedom, therefore, to find the exact solution to each and every problem is quite difficult. In this regard, great motivation has been observed in the last two decades to establish best tools to handle such problems. One of the important techniques is to find analytical approximate solutions to many nonlinear problems of FODEs. For this purpose, the usual decompositions, perturbation, and integral transform methods were greatly utilized to investigate ordinary differential and integral equations. The mentioned techniques have been very well explored for fractional differential equations (see [16–25]). One of the powerful methods which has been used very frequently in the past is due to the Laplace Adomian decomposition method. For usual FODEs, the mentioned method has been used very regularly in the literature. However, to the best of our information, the aforesaid method is very rarely used for FODEs involving nonsingular kernel, see [26, 27].

Here, we remark that mathematical models are the powerful tools to study various dynamical problems of physical and biological sciences. The concept was initiated by Bernoulli in 1776. However, a formal mathematical model of three compartments was constructed in 1927 by McKendrick and his co-author called the SIR mode. Later on, the subject of mathematical modeling was extended to infectious diseases. Because mathematical models of the biological problem have become powerful tools to understand various infectious diseases, the proper method has to be developed to control the disease or minimize its transmission in the society. From ancient time, measles disease is one of the most threatful diseases. It was a big reason for children death in the past. This dangerous disease was caused by germs called morbilliform. Measles disease spreeds, when an infected person coughs or sneezes, because its virus can live for up to two hours in an airspace where the infected person coughed or sneezed. Measles-infected individuals can transfer their germs to the other people 4–8 days before and after the skin eruptions start. It transmits a disease in young babies up to 30–40 million every year. Measles appears once and is present for a long time in life for immunity. Symptoms of the disease include runny nose, high temperature, coughing, and spots on the whole body and in highly complicated cases, ear infections, diarrhea, and pneumonia. The vaccination of measles has been used to control disease in kids. Some vaccinated individuals could remain susceptible individuals when vaccination gets failed. Worldwide vaccination reduced death caused by measles between 2000 and 2017. However, measles disease is still familiar disease in highly developing countries of Asia and Africa particularly due to the lack of proper treatment of this infectious disease. For this purpose, a massive research has been carried out to enhance the understanding of the virus of measles dynamics in various areas. For example, the authors in [28–33] discussed the global stability of the model with five compartments as

where represents the susceptible, represents the vaccinated, represents the exposed, represents the infected, and represents the recovered individuals. The description of the parameters is given in the analytical section. These models have been investigated corresponding to ordinary and usual fractional-order derivatives. Furthermore, the researchers have given the global and local dynamics by computing the basic reproductive numbers. Here, we investigate the given model under the nonsingular fractional derivative of Caputo and Fabrizio (CFFD) from other perspectives including the qualitative analysis by using fixed point approach. Further stability is a required aspect in the dynamical problem. Since we are going to derive approximate solutions, therefore, in this regard, Hyers–Ulam-type stability results are investigated. The mentioned stability has been very well studied for the general problem of FODEs, see, for detail, [5, 34–36]. Also, the analytical results are investigated through the Laplace Adomian decomposition method. We considered model (1) under the CFFD with fractional order as

under the initial conditions

We obtain the solution in the form of series for the considered problem. Also, we display the results against different values of fractional order . Also, we provide results about the existence and uniqueness of the solution for the concerned model by using fixed point theorems due to Schauder and Banach. Here, we remark that we use the Laplace Adomian decomposition method because this method is easy and efficient and less expensive. Furthermore, the mentioned method does not require any predefine step size or controlling parameter which are needed by RK methods or the homotopy method, respectively, for detail, see [37–39]. Furthermore, the convergence of the method has been showed in many papers of the proposed method, for instance, see [27].

2. Preliminaries

Definition 1. (see [3]). Let , , and ; then, the CFFD is recalled as

is the normalization function defined as and satisfies the conditions . If the function fails to exist in , then the derivative can be redefined as

Definition 2. (see [3]). The integral of fractional order of a function is defined by

Lemma 1. (see [40]). The solution of the differential equation involving CFFD for as

is given by

Definition 3. (see [34, 35]). The solution of the differential equation

is Hyers–Ulam stable; if there exists such that, for the inequality

there exists a unique solution with a constant such that

Furthermore, if there exists nondecreasing function such that (11) may be written as

then the concerned solution of problem (9) is generalized Hyers–Ulam stable.

Theorem 1 (see [40]). Let be a convex subset of Banach space , with operators and with

(1) for all (2) is the condensing operator(3) is continues and compact

Then, there exists at least one solution such that

Definition 4. (see [4, 14]). The Laplace transform of CFFD is given as

3. Existence and Stability Results for the Considered Model

In this part of the manuscript, we determine existence results for model (2) using the fixed point theorem due to Banach. In this regard, we first define the following functions:

Then, we write some notions for easiness as

and

Using (16), system (2) can be written as

In view of Lemma 1, problem (18) can be converted to the given integral equation as

Now, let , and let ; we define the Banach space under the norm

The assumptions given in the following hold true:

  Under the continuity of , there exists such that

  There exists positive constant such that .

Theorem 2. Under the assumptions , problem (19) has at least one solution if ; consequently, the considered system (2) has at least one solution.

Proof. Let be a closed convex set. Now, we define the operators from (19) as

To derive the required results, we first prove that operator is a contraction.
For any , one has

This shows that is a contraction. Now, to show that is bounded and equicontinuous, the continuity of implies that is continuous. For any , we have

This shows that is bounded; for equicontinuity, let , and we have

which implies that as . So, is uniformly continuous and bounded. Thus, by Arzelá–Ascoli theorem, is relatively compact and so is completely continuous. Thus, by Theorem 1, problem (18) has at least one solution. Consequently, the considered model (2) has at least one solution.

Next, we define the operator by

Theorem 3. Under assumption , the operator as defined in (26) is a contraction; then, problem (18) has a unique solution with the condition , and consequently, our proposed system (2) has a unique solution.

Proof. Let ; then, from (26), one has

This shows that is a contraction. Therefore, problem (18) has a unique solution. Hence, our considered system (2) has a unique solution.

Remark 1. Next, for stability analysis, we consider a small perturbation , such that , depends only on the solution.

(i) for (ii)

Lemma 2. The solution of the perturbed problem

satisfies the given relation:

where .

Proof. In view of Lemma 1, the solution of perturb problem (28) is given by

From (30), on using Remark 1, we have

Theorem 4. Under assumption , together with Lemma 2, problem (18) is Hyers–Ulam stable if , which yields that our considered system (2) is Hyers–Ulam stable.

Proof. Let be any solution, and is a unique solution; then,

This implies that

Therefore, the solution of (18) is Hyers–Ulam stable. Hence, the solution of the proposed system (2) is Hyers–Ulam stable.

Remark 2. In the same line, we can also develop the results of generalized Hyers–Ulam, Rassias–Hyers–Ulam stability. The aforementioned stability analysis has been studied for simple mathematical models of biology and physics in [41–43].

4. Derivation of the General Semianalytical Solution to the Considered Model (2)

Here, in this section, we are going to compute series solution for the suggested problem. To receive this goal, taking Laplace transform of (16), we have

Now, assuming the solution in the series form,

Further decomposing the nonlinear terms , etc. in terms of Adomian polynomials,

where the Adomian polynomial can be defined as

In the same way, the other polynomials can be defined.

Hence, in view of (35) and (36), system (34) becomes

Now, equating terms on both sides of (38), we have

Evaluating the Laplace transform in (39), we get

and so on. Therefore, we get the required solution as given by

Note. For the convergence of the proposed method, see the paper [27].

5. Results and Discussion

This part of the manuscript is related to provide numerical results and some discussion about the approximate solution of the considered problem. For this purpose, we apply the Laplace Adomian decomposition methods for the solution. We select appropriate values for parameters, as we take [3]  , , , , , , , , , , , and . In view of these values, we draw the graph of approximate solutions (40) for first ten terms by using Matlab in Figures 1–5 against various fractional-order derivatives.

Figure 1: Plot of the susceptible class at different fractional values of .

Figure 2: Plot of the vaccinated class at different fractional values of .

Figure 3: Plot of the exposed class at different fractional values of .

Figure 4: Plot of the infected class at different fractional values of .

Figure 5: Plot of the recovered class at different fractional values of .

From Figures 1–5, we see that the susceptible class population is decreasing with different rates. It decreases with faster speed when the order is smaller as compared with the larger order. Also, the dynamics of the vaccinated class is increasing with different rates due to the fractional order. Also, the exposed population is increasing up, and hence, on using the vaccine, the density of the population of the infected class is decreasing, while the density of the recovered class increases with the same scenario. From all these figures, we concluded that the fractional derivative with exponential kernel can also be used to provide the global dynamics of the considered model.

6. Conclusion

By fixed point approach, we have investigated the existence of the considered model under CFFD. Also, by using the tools of nonlinear functional analysis, we have established sufficient conditions for Hyers–Ulam-type stability of the approximate solutions of the considered model. Also, we have provided the semianalytical solution to the considered model by the Laplace Adomian decomposition method. The concerned method needs no discretization of data nor extra axillary parameter as needed by homotopy methods on which these methods depend. The proposed method has been utilized extensively for usual fractional differential equations, but in the case of new fractional differential operators, this method has not been properly used. Hence, we concluded that CFFD can also be used as powerful tools to investigate biological models of infectious diseases. Also, the concerned models involving CFFD can be handled easily by using the Laplace Adomian decomposition method.

Data Availability

No data were used to support this study.

Conflicts of Interest

There exist no conflicts of interest regarding this research work.

Authors’ Contributions

All authors have equal contribution to this work.

Acknowledgments

The authors are thankful to the HED for supporting this work under grant no. HEREF-46 and HEC of Pakistan under grant no. NRPU-10039. This work was supported financially by the National Natural Science Foundation of China (11971493).

References

I. Podlubny, Fractional Differential Equations: Mathematics in Science and Engineering, Academic Press, New York, NY, USA, 1999.

A. A. Kilbas, H. Srivastava, and J. Trujillo, “Theory and application of fractional differential equations,” North Holland Mathematics Studies, vol. 204, Elsevier, Amsterdam, Netherlands, 2006.

M. Caputo and M. Fabrizio, “A new definition of fractional derivative with out singular kernel,” Progress in Fractional Differentiation and Applications, vol. 1, pp. 73–85, 2015.

M. Caputo and M. Fabrizio, “Applications of new time and spatial fractional derivatives with exponential kernels,” Progress in Fractional Differentiation and Applications, vol. 2, no. 1, pp. 1–11, 2016.

N. Sene, “Stability analysis of electrical RLC circuit described by the Caputo-Liouville generalized fractional derivative,” Alexandria Engineering Journal, 2020.

K. M. Saad, “New fractional derivative with non-singular kernel for deriving Legendre spectral collocation method,” Alexandria Engineering Journal, 2019.

K. M. Saad, D. Baleanu, and A. Atangana, “New fractional derivatives applied to the Korteweg-de Vries and Korteweg-de Vries-Burger’s equations,” Computational and Applied Mathematics, vol. 37, no. 4, pp. 5203–5216, 2018.

K. M. Saad, S. Deni̇z, and D. Baleanu, “On a new modified fractional analysis of Nagumo equation,” International Journal of Biomathematics, vol. 12, no. 3, Article ID 1950034, 2019.

H. M. Srivastava and K. M. Saad, “New approximate solution of the time-fractional Nagumo equation involving fractional integrals without singular Kernel,” Applied Mathematics & Information Sciences, vol. 14, no. 1, pp. 1–8, 2020.

H. A. A. El-Saka, “The fractional-order SIS epidemic model with variable population size,” Journal of the Egyptian Mathematical Society, vol. 22, no. 1, pp. 50–54, 2014.

T. Hernandez, “A fractional calculus approach to the dynamic optimization of biological reactive systems. Part I: fractional models for biological reactions,” Chemical Engineering Science, vol. 117, pp. 217–228, 2014.

Z. Wang, D. Yang, T. Ma, and N. Sun, “Stability analysis for nonlinear fractional-order systems based on comparison principle,” Nonlinear Dynamics, vol. 75, no. 1-2, pp. 387–402, 2014.

A. Atangana and J. F. Gómez-Aguilar, “Fractional derivatives with no-index law property: application to chaos and statistics,” Chaos, Solitons & Fractals, vol. 114, pp. 516–535, 2018.

A. Atangana and B. S. T. Alkahtani, “New model of groundwater flowing within a confine aquifer: application of Caputo-Fabrizio derivative,” Arabian Journal of Geosciences, vol. 9, no. 1, pp. 1–6, 2016.

H. Khan, J. F. Gömez-Aguilar, A. Alkhazzan, and A. Khan, “A fractional order HIV-TB coinfection model with nonsingular Mittag-Leffler law,” Mathematical Methods in the Applied Sciences, vol. 43, pp. 1–21, 2020.

F. Haq, K. Shah, G. ur Rahman, and M. Shahzad, “Numerical solution of fractional order smoking model via Laplace Adomian decomposition method,” Alexandria Engineering Journal, vol. 57, no. 2, pp. 1061–1069, 2018.

A. Ali, K. Shah, and R. A. Khan, “Numerical treatment for traveling wave solutions of fractional Whitham-Broer-Kaup equations,” Alexandria Engineering Journal, vol. 57, no. 3, pp. 1991–1998, 2018.

O. Kiymaz, “An algorithm for solving initial value problems using Laplace Adomian decomposition method,” Applied Mathematical Sciences, vol. 3, no. 30, pp. 1453–1459, 2009.

V. Morales-Delgado, J. F. Gömez-Aguilar, and M. A. Taneco-Hernandez, “Analytical solution of the time fractional diffusion equation and fractional convection-diffusion equation,” Revista Mexicana de Física, vol. 65, no. 1, pp. 82–88, 2018.

J. F. Gömez-Aguilar, H. Yépez-Martínez, R. F. Escobar-Jiménez, V. H. Olivares-Peregrino, J. M. Reyes, and I. O. Sosa, “Series solution for the time-fractional coupled mKdV equation using the homotopy analysis method,” Mathematical Problems in Engineering, vol. 2016, Article ID 7047126, 8 pages, 2016.

K. M. Saad, M. M. Khader, J. F. Gömez-Aguilar, and D. Baleanu, “Numerical solutions of the fractional Fisher’s type equations with Atangana-Baleanu fractional derivative by using spectral collocation methods,” Chaos: An Interdisciplinary Journal of Nonlinear Science, vol. 29, no. 2, Article ID 023116, 2019.

J. F. Gömez-Aguilar and A. Atangana, “Fractional Hunter-Saxton equation involving partial operators with bi-order in Riemann-Liouville and Liouville-Caputo sense,” The European Physical Journal Plus, vol. 132, no. 2, p. 100, 2017.

N. Sene and K. Abdelmalek, “Analysis of the fractional diffusion equations described by Atangana-Baleanu-Caputo fractional derivative,” Chaos, Solitons & Fractals, vol. 127, pp. 158–164, 2019.

N. Sene, “Second-grade fluid model with Caputo-Liouville generalized fractional derivative,” Chaos, Solitons & Fractals, vol. 133, Article ID 109631, 2020.

A. A. Alderremy, K. M. Saad, P. Agarwal, S. Aly, and S. Jain, “Certain new models of the multi space-fractional Gardner equation,” Physica A: Statistical Mechanics and Its Applications, vol. 545, Article ID 123806, 2020.

K. Shah, M. A. Alqudah, F. Jarad, and T. Abdeljawad, “Semi-analytical study of Pine Wilt disease model with convex rate under Caputo-Febrizio fractional order derivative,” Chaos, Solitons & Fractals, vol. 135, Article ID 109754, 2020.

K. Shah, T. Abdeljawad, and F. Jarad, “On a nonlinear fractional order model of dengue fever disease under Caputo-Fabrizio derivative,” Alexandria Engineering Journal, 2020.

R. T. Prey and N. A-Hasley, “The clinical significance of measles; a review,” The Journal of Infected Disease, vol. 189, no. 1, pp. 4–16, 2004.

B. P. Prawato, “Stability and simulation of measles transmission model with and without vaccination,” For East Journal of Mathematical Sciences, vol. 102, no. 2, pp. 271–281, 2017.

E. O. Oghre and I. I. Ako, “A mathematical model for measles disease,” For East Journal of Mathematical Sciences, vol. 54, no. 1, pp. 47–63, 2011.

K. Ejima, R. Omori, K. Aihara, and H. Nishiura, “Real-time investigation of measles epidemics with estimate of vaccine efficacy,” International Journal of Biological Sciences, vol. 8, no. 5, pp. 620–629, 2012.

S. Okyere-Siabouh and L. A. Adetunde, “Mathematical model for the study of measles in cape coast metropoli,” International Journal of Modern Biology and Medicines, vol. 4, no. 2, pp. 110–113, 2013.

G. Bolarin, “On the dynamical analysis of a new model for measles infection,” International Journal of Mathematical Trends and Technology, vol. 7, no. 2, pp. 2231–5373, 2014.

M. Sher, K. Shah, M. Fečkan, and R. A. Khan, “Qualitative analysis of multi-terms fractional order delay differential equations via the topological degree theory,” Mathematics, vol. 8, no. 2, p. 218, 2020.

M. Benchohra and S. Bouriah, “Existence and stability results for nonlinear boundary value problem for implicit differential equations of fractional order,” Moroccan Journal of Pure and Applied Analysis, vol. 1, no. 1, pp. 22–37, 2015.

H. Khan, A. Khan, W. Chen, and K. Shah, “Stability analysis and a numerical scheme for fractional Klein-Gordon equations,” Mathematical Methods in the Applied Sciences, vol. 42, no. 2, pp. 723–732, 2019.

S. Paul, S. P. Mondal, and P. Bhattacharya, “Numerical solution of Lotka Volterra prey predator model by using Runge-Kutta-Fehlberg method and Laplace Adomian decomposition method,” Alexandria Engineering Journal, vol. 55, no. 1, pp. 613–617, 2016.

M. Z. Mohamed and T. M. Elzaki, “Comparison between the Laplace decomposition method and adomian decomposition in time-space fractional nonlinear fractional differential equations,” Applied Mathematics, vol. 9, no. 4, pp. 448–458, 2018.

Y. Khan, “An efficient iterated method for mathematical biology model,” Neural Computing and Applications, vol. 23, no. 3-4, pp. 677–682, 2013.

Eiman, K. Shah, M. Sarwar, and D. Baleanu, “Study on Krasnoselskii’s fixed point theorem for Caputo-Fabrizio fractional differential equations,” Advances in Difference Equations, vol. 2020, no. 1, 2020.

E. Ahmed, “On applications of Ulam-Hyers stability in biology and economics,” 2010, https://arxiv.org/abs/1004.1354.

A. Atangana and J. F. Gomez-Aguilar, “Decolonisation of fractional calculus rules: breaking commutativity and associativity to capture more natural phenomena,” The European Physical Journal Plus, vol. 133, p. 166, 2018.

T. Abdeljawad and D. Baleanu, “On fractional derivatives with exponential kernel and their discrete versions,” Reports on Mathematical Physics, vol. 80, no. 1, pp. 11–27, 2017.