/ / Article

Research Article | Open Access

Volume 2019 |Article ID 3913840 | https://doi.org/10.1155/2019/3913840

M. Valizadeh, Y. Mahmoudi, F. Dastmalchi Saei, "Application of Natural Transform Method to Fractional Pantograph Delay Differential Equations", Journal of Mathematics, vol. 2019, Article ID 3913840, 9 pages, 2019. https://doi.org/10.1155/2019/3913840

# Application of Natural Transform Method to Fractional Pantograph Delay Differential Equations

Accepted30 Sep 2019
Published03 Nov 2019

#### Abstract

In this paper, a new method based on combination of the natural transform method (NTM), Adomian decomposition method (ADM), and coefficient perturbation method (CPM) which is called “perturbed decomposition natural transform method” (PDNTM) is implemented for solving fractional pantograph delay differential equations with nonconstant coefficients. The fractional derivative is regarded in Caputo sense. Numerical evaluations are included to demonstrate the validity and applicability of this technique.

#### 1. Introduction

Delay differential equations (DDEs) occur when the rate of change of process is specified by its state at a certain past state known as its history. These types of differential equations occur in traffic models, control systems, population dynamics, and many natural phenomena [13]. The expression of delay in mathematical modeling of real-life phenomena leads to a precise and acceptable description of the dynamics of the model. DDEs are complex in nature because the history of the system over the delayed interval is given as an initial condition. Due to this reason, these equations are difficult to solve analytically; hence, a numerical method is required. Pantograph delay equations are one of the well-known functional differential equations with proportional delay and often appear in many scientific models such as mechanics and electrodynamics [4, 5]. In 1851, it was the first time that a device named “pantograph” was used in construction of the electric locomotive after which this name was originated from that time. Pantograph was modeled mathematically in 1971 (see [6] for more details). Thereafter, pantograph delay equation was studied by many authors and solved by several numerical methods. The most important of them are collocation method [7], spline method [8], Runge–Kutta method [9], homotopy perturbation method [10], and Adomian decomposition method [11]. In this paper, we intend to apply the natural transform method to solve fractional delay differential equations with variable coefficients. To solve differential and integral equations, several integral transforms such as Fourier, Laplace, and Sumudu are used [1216]. The natural transform method is a new integral transform was introduced by Khan and Khan [17], and its properties were described. Belgacem used the natural transform method to solve Bessel’s differential equation, Maxwell’s equation, and nonlinear Klein–Gordon equations (see [18, 19] and its references). A table for some natural transformation properties on certain functions was presented by Khan and Khan in [17]. Moreover, this method is applied to find the solution of diffusion equations [20]. Finally, they proved that the natural transform method converges to the Laplace and Sumudu transforms.

In this paper, a novel technique is applied to find an approximate solution for the fractional pantograph delay differential equation of the following form:where , and are the given analytic functions, and denotes the Caputo fractional differential operator. The proposed method is a combination of natural transform method (NTM), Adomian decomposition method (ADM) [21], and coefficient perturbation method (CPM) [22], which is called “perturbed decomposition natural transform method.” Unlike existing analytical and numerical methods, this method gives an approximate and/or exact solution for problem (1) by a simple calculation and in an elegant way.

The structure of this paper is as follows: In Section 2, a brief description of the fractional calculus, natural transform method, and Adomian decomposition method are provided. The existence of unique solution for problem (1) is investigated in Section 3. In Section 4, we introduce our new method (called briefly, PDNTM) for obtaining the approximate solution for problem (1). Numerical examples illustrate efficiency and applicability of the proposed method in Section 5, and conclusions are drawn in Section 6.

#### 2. Preliminaries

In this section, we recall some preliminary results which will be needed throughout the paper.

##### 2.1. Fractional Calculus

There are some kinds of definitions for fractional derivatives and integrals, such as Riemann–Liouville integral [23], which are described as follows.

For , the Riemann–Liouville fractional integral of order α is defined as

For , set , the identity operator. Let for some and . Then,

Let . The operator defined byis called the Riemann–Liouville fractional differential operator of order α. For , set , the identity operator. The Caputo differential operator of order n is defined bywhenever exists, where denotes the Taylor polynomial of degree for the function f around the point In the case , define . Under the above conditions, it is easy to show that

More details can be known by referring to [23, 24].

##### 2.2. Natural Transform Method (NTM)

In this section, we review some properties of the atural transform method that will be used in later sections. For the function defined in , piecewise function continues in every finite interval and whereThe natural transform of is defined as follows:With the above conditions, the right-hand side integral in (39) is convergent. By taking in (39), the Laplace transform of ,and for , the Sumudu transform of ,will be obtained. The inverse natural transform of is defined in [19] by

We recall some useful properties of the natural transform method from [19] in Table 1.

 Functional form Natural transform form (1) af(t) + bg(t) (2) 1 (3) (4) (5)

Consider the general form of a differential equation aswhere is the linear invertible part, corresponds to the nonlinear term, and is the remaining part. Solving from (12), we have

As L is invertible,

The Adomian decomposition method (ADM) [21] is a powerful tool for solving both linear and nonlinear functional equations. This method decomposes the solution into a rapidly convergent series of solution components and then decomposes the analytic nonlinearity into the series of the Adomian polynomials:where are the well-known Adomian polynomials as follows:

Substituting (15) in (14), we have

Consequently, we write

The convergence analysis of the Adomian decomposition method was studied by Cherruault in [25, 26].

#### 3. Existence and Uniqueness of Solution

In this section, we prove existence and uniqueness of the solution for the fractional delay differential equation (1). This study is based on the Banach contraction principle.

Theorem 1. Assume and are the given analytic functions in (1) and . Ifthen problem (1) has a unique solution .

Proof. The equivalent integral equation representation of problem (1) is of the following form (see [24]):Based on the Banach contraction principle, it is enough to prove that the operator A defined byhas a unique fixed point. It is clear that the operator A maps into itself and thatwhich implies, under our assumption, that A is a contraction mapping. Then, the Banach contraction principle implies that A has a unique fixed point .

#### 4. Main Results

Let us approximate the coefficients and , in equation (1) by polynomials of fractional order α as follows:

Then, problem (1) can be written in perturbed form aswhich is called “perturbed problem” with the exact solution that is an approximation for the exact solution of the main problem (1). This is a first step of our new method that converts the main problem (1) to new one with polynomial coefficients. This strategy will simplify the calculation in the implementation of the proposed method. Similar ideas have been used in [27, 28] for solving boundary value problem and eigenvalue problems.

Theorem 2. Let andfor some constant in (23). Then, we have , i.e.,

Proof. It can be easily proved that is a solution of equation (1) if and only if is a solution for the integral equationwhere is the Riemann–Liouville fractional differential operator of order α (see [24]). Consequently, for perturbed problem (24), we haveSubtracting equation (28) from equation (27) and some manipulation yieldBy taking from both sides of equation (29), we obtainFrom (25), we haveThen, from (30) and (31), we getConsequently,After some manipulation on (33), we obtainFrom the upper bound in (34), it is seen that if , then .

##### 4.1. Method of Solution

Taking the natural transform of both sides of equation (24) yields

Using the properties of the natural transform method in Table 1, we get

Consequently,

By taking the inverse natural transform of (37), we obtain

Now, based on the Adomian decomposition method, we construct the solution in (38) by an infinite convergent series of functions as follows:

Substituting (39) in (38) yieldsin which the components will be obtained by the following recurrence relation:such that

The -term approximation of the solution isso that

Therefore, the solution comes in the form of a convergent series with easily computed components. From Theorem 2 and (23), we know that as .

Remark 1. As is clear, in recurrence relation (41), the terms and can be calculated easily using properties (3) in Table 1, since the term is a fractional-order polynomial for and . This is a main advantage in our proposed method for calculating the components in (43) and making our new method easy to be used and writing a computer program for it.
We summarize the steps of our proposed method (called PDNTM method) as an implementation algorithm as follows (Algorithm 1).

 Input Choose values for parameters ρ, in (23), and M in (44) Step 1. Find the approximation of functions and in (23) Step 2. Compute the components from recurrence relation (41) for Output Construct the approximate solution

#### 5. Numerical Results

In this section, we present some examples to support the accuracy and simplicity of the proposed method. These results in this section indicate the applicability of the method and verification of the theoretical results in the previous section.

Example 1. We consider the following initial value problem:The exact solution for is . According to the proposed method, by taking natural transform to both sides of equation (45), we obtainFrom initial condition , we obtainApplying the inverse of natural transform in (47) yieldsNow, based on (41), we obtainwhich givesThen, all of the components can be obtained as follows:Continuing the above recurrence relation yieldsFor the case , the second term in (52) is equal to zero; therefore, , which is the exact solution of problem (45). For , the approximate solution isIn Figure 1, the convergence of approximate solutions to exact solution as is shown.

Example 2. Consider the following fractional pantograph problem (see [29]):In the case , the exact solution is . According to the proposed method for , we obtainThen, from (43), we obtainwhich is the exact solution. The approximate solutions for different values of α and exact solution are shown in Figure 2 for . Note that as , the numerical solutions tends to the exact solution .

Example 3. We consider the following nonpolynomial coefficients problem (see [29]):The exact solution for is . LetThen, the perturbed problem is as follows:The approximate solutions for and exact solution for , , and are shown in Figure 3. Table 2 denotes the behavior of approximate solution for different values of with . These results are consistent with Theorem 2 which shows convergence of the solution of perturbed problem (59) to (57).

 t Exact. Sol App. Sol Abs. Err For 0.0 1.0000000000 1.0000000000 0 0.1 1.1051709181 1.1051873941 0.2 1.2214027582 1.2214234795 0.3 1.3498588076 1.3498811112 0.4 1.4918246976 1.4918490746 0.5 1.6487212707 1.6487484081 0.6 1.8221188004 1.8221489350 0.7 2.0137527075 2.0137858850 0.8 2.2255409285 2.2255774493 0.9 2.4596031112 2.4596436807 1.0 2.7182818285 2.7183270070 For 0.0 1.0000000000 1.0000000000 0 0.1 1.1051709181 1.1051709163 0.2 1.2214027582 1.2214027562 0.3 1.3498588076 1.3498588054 0.4 1.4918246976 1.4918246952 0.5 1.6487212707 1.6487212680 0.6 1.8221188004 1.8221187975 0.7 2.0137527075 2.0137527042 0.8 2.2255409285 2.2255409249 0.9 2.4596031112 2.4596031072 1.0 2.7182818285 2.7182818240 For 0.0 1.0000000000 1.0000000000 0 0.1 1.1051709181 2.7182818285 0.2 1.2214027582 1.1051709181 0.3 1.3498588076 1.3498588076 0.4 1.4918246976 1.4918246976 0.5 1.6487212707 1.6487212707 0.6 1.8221188004 1.8221188004 0.7 2.0137527075 2.0137527075 0.8 2.2255409285 2.2255409285 0.9 2.4596031112 2.4596031112 1.0 2.7182818285 2.7182818285

#### 6. Future Work

In our future work, the PDNTM method will be studied for the general class of problems such as the system of fractional pantograph delay differential equations of the formsubject to the initial conditionswhere and . To this end, we can approximate the coefficients and by polynomials of fractional order as follows:and extend all processes of Section 4.1.

#### 7. Conclusion

In this paper, a new method based on the combination of the natural transform method, Adomian decomposition method, and coefficient perturbation method is introduced to solve the fractional pantograph delay differential equations with nonconstant coefficients. From this method, we obtain an approximate or exact solution for the problem with some simple iterative calculation that is easy in use and computer programming. The numerical results confirm the accuracy of the method.

#### Data Availability

No data were used to support this study.

#### Conflicts of Interest

The authors declare that they have no conflicts of interest.

#### Acknowledgments

The authors appreciate the reviewer’s invaluable advice on improving the paper.

#### References

1. E. Fridman, L. Fridman, and E. Shustin, “Steady modes in relay control systems with time delay and periodic disturbances,” Journal of Dynamic Systems, Measurement, and Control, vol. 122, no. 4, pp. 732–737, 2000. View at: Publisher Site | Google Scholar
2. I. R. Epstein and Y. Luo, “Differential delay equations in chemical kinetics. Nonlinear models: the cross-shaped phase diagram and the Oregonator,” The Journal of Chemical Physics, vol. 95, no. 1, pp. 244–254, 1991. View at: Publisher Site | Google Scholar
3. Y. Kuang, Delay Dierential Equations with Applications in Population Biology, Academic Press, Boston, MA, USA, 1993.
4. C. T. H. Baker and E. Buckwar, “Numerical analysis of explicit one-step methods for stochastic delay differential equations,” LMS Journal of Computation and Mathematics, vol. 3, pp. 315–335, 2000. View at: Publisher Site | Google Scholar
5. J. J. Zhao, Y. Xu, H. X. Wang, and M. Z. Liu, “Stability of a class of Runge–Kutta methods for a family of pantograph equations of neutral type,” Applied Mathematics and Computation, vol. 181, no. 2, pp. 1170–1181, 2006. View at: Publisher Site | Google Scholar
6. J. R. Ockendon and A. B. Tayler, “The dynamics of a current collection system for an electric locomotive,” Proceedings of the Royal Society of London A: Mathematical, Physical and Engineering Sciences, vol. 322, pp. 447–468, 1971. View at: Publisher Site | Google Scholar
7. Y. Muroya, E. Ishiwata, and H. Brunner, “On the attainable order of collocation methods for pantograph integro-differential equations,” Journal of Computational and Applied Mathematics, vol. 152, no. 1-2, pp. 347–366, 2003. View at: Publisher Site | Google Scholar
8. M. Shadia, “Numerical solution of delay differential and neutral differential equations using spline methods,” Assuit University, Asyut, Egypt, 1992, Ph. D. thesis. View at: Google Scholar
9. D. Li and M. Z. Liu, “Runge–Kutta methods for the multi-pantograph delay equation,” Applied Mathematics and Computation, vol. 163, no. 1, pp. 383–395, 2005. View at: Publisher Site | Google Scholar
10. E. Yusufoğlu, “An efficient algorithm for solving generalized pantograph equations with linear functional argument,” Applied Mathematics and Computation, vol. 217, no. 7, pp. 3591–3595, 2010. View at: Publisher Site | Google Scholar
11. D. J. Evans and K. R. Raslan, “The Adomian decomposition method for solving delay differential equation,” International Journal of Computer Mathematics, vol. 82, no. 1, pp. 49–54, 2005. View at: Publisher Site | Google Scholar
12. F. B. M. Belgacem, “Applications of the Sumudu transform to indefinite periodic parabolic equations,” in Proceedings of the 6th International Conference on Mathematical Problems and Aerospace Science, vol. 51–60, Cambridge, UK, 2007. View at: Google Scholar
13. F. B. M. Belgacem, “Sumudu applications to maxwell’s equations,” PIERS Online, vol. 5, no. 4, pp. 355–360, 2009. View at: Publisher Site | Google Scholar
14. M. A. Asiru, “Sumudu transform and the solution of integral equations of convolution type,” International Journal of Mathematical Education in Science and Technology, vol. 32, no. 6, pp. 906–910, 2001. View at: Publisher Site | Google Scholar
15. Q. D. Katabeh and F. B. M. Belgacem, “Applications of the Sumudu transform to fractional differential equations,” Nonlinear Studies, vol. 18, no. 1, pp. 1–15, 2011. View at: Google Scholar
16. T. M. Elzaki and S. M. Elzaki, “On the solution of integro-differential equation system using Elzaki transform,” Global Journal of Mathematics and Mathematical Sciences, vol. 3, no. 1, pp. 13–23, 2011. View at: Google Scholar
17. 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
18. F. B. M. Belgacem and R. Silambarasan, “Maxwell’s equations solutions through the natural transform,” Mathematics in Engineering, Science and Aerospace, vol. 3, no. 3, pp. 313–323, 2012. View at: Google Scholar
19. F. B. M. Belgacem and R. Silambarasan, “Theory of natural transform,” International Journal of Mathematics in Engineering, Science and Aerospace, vol. 3, no. 1, pp. 99–124, 2012. View at: Google Scholar
20. 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
21. G. Adomian, Solving Frontier Problems of Physics: The Decomposition Method, Kluwer Academic Publisher, Boston, MA, USA, 1994.
22. S. A. Pruess, “Solving linear boundary value problems by approximation of coefficients,” Mathematics of Computation, vol. 123, no. 27, 1973. View at: Google Scholar
23. I. Podlubny, Fractional Differential Equations, Academic Press, San Diego, CA, USA, 1999.
24. K. Diethelm, The Analysis of Fractional Differential Equations, Springer, Berlin, Germany, 2010.
25. K. Abbaoui and Y. Cherruault, “Convergence of Adomian’s method applied to nonlinear equations,” Mathematical and Computer Modelling, vol. 20, no. 9, pp. 69–73, 1994. View at: Publisher Site | Google Scholar
26. N. Himoun, K. Abbaoui, and Y. Cherruault, “New results of convergence of Adomian’s method,” Kybernetes, vol. 28, no. 4, pp. 423–429, 1999. View at: Publisher Site | Google Scholar
27. J. Canosa and R. G. De Oliveira, “A new method for the solution of the Schrödinger equation,” Journal of Computational Physics, vol. 5, no. 2, pp. 188–207, 1970. View at: Publisher Site | Google Scholar
28. M. H. Alexander and R. G. Gordon, “New method for constructing solutions to time-dependent perturbation equations,” The Journal of Chemical Physics, vol. 55, no. 10, pp. 4889–4897, 1971. View at: Publisher Site | Google Scholar
29. S. Widatalla, “A comparative study on the stability of Laplace-Adomian algorithm and numerical methods in generalized pantograph equations,” ISRN Computational Mathematics, vol. 2012, Article ID 704184, 6 pages, 2012. View at: Publisher Site | Google Scholar