Advances in Mathematical Physics

Volume 2014 (2014), Article ID 595848, 8 pages

http://dx.doi.org/10.1155/2014/595848

Research Article

A Shifted Jacobi-Gauss Collocation Scheme for Solving Fractional Neutral Functional-Differential Equations

A. H. Bhrawy1,2 and M. A. Alghamdi¹

¹Department of Mathematics, Faculty of Science, King Abdulaziz University, Jeddah 21589, Saudi Arabia

²Department of Mathematics, Faculty of Science, Beni-Suef University, Beni-Suef 62511, Egypt

Received 4 March 2014; Revised 29 March 2014; Accepted 29 March 2014; Published 15 April 2014

Academic Editor: Xiao-Jun Yang

Copyright © 2014 A. H. Bhrawy and M. A. Alghamdi. 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 shifted Jacobi-Gauss collocation (SJGC) scheme is proposed and implemented to solve the fractional neutral functional-differential equations with proportional delays. The technique we have proposed is based upon shifted Jacobi polynomials with the Gauss quadrature integration technique. The main advantage of the shifted Jacobi-Gauss scheme is to reduce solving the generalized fractional neutral functional-differential equations to a system of algebraic equations in the unknown expansion. Reasonable numerical results are achieved by choosing few shifted Jacobi-Gauss collocation nodes. Numerical results demonstrate the accuracy, and versatility of the proposed algorithm.

1. Introduction

Fractional differential equations (FDEs) have drawn the interest of many researchers in recent years [1–6], due to their useful applications in many fields of science. In fact, we may observe several applications in electrochemistry, viscoelasticity, electromagnetic, control, plasma physics, porous media, fluctuating environments, dynamical processes, and so on. In consequence, fractional differential equations are gaining much attention from the researchers. For some recent developments on this subject, see [7–15].

In the last decade or so, comprehensive research has been accomplished on the development of numerical algorithms which are numerically stable for both linear and nonlinear FDEs. Tripathi et al. [16] presented a new operational matrix of hat functions to solve linear FDEs. The spectral tau method was proposed in [17] to achieve an accurate solution of linear and nonlinear FDEs subject to multipoint conditions. In [18], the author proposed Bernstein polynomial to design a numerical algorithm for fractional Riccati equations. The authors of [19] investigated the spline collocation method for approximating the solution of nonlinear FDEs. Furthermore, the author of [20] transformed the time-dependent space FDE with variable coefficients into a system of ordinary differential equations, which is then solved by a standard numerical method. Baleanu et al. [21] developed the generalized Laguerre spectral tau and collocation approximations to solve FDEs on the half line. In [22], Ma and Huang developed spectral collocation method for solving linear fractional integrodifferential equations. Yang and Huang [23] analyzed and developed the Jacobi collocation scheme for pantograph integrodifferential equations with fractional orders in finite interval. In [24] Yin et al. proposed a new fractional-order Legendre function with spectral method to solve partial FDEs; based on the operational matrix of these functions, the same authors developed their approach in combination with variational iteration formula to solve a class of FDEs; see [25]. More recently, the Jacobi Galerkin method was extended in [26] to solve stochastic FDEs.

Polynomial approximations can be quite useful for expressing the solution of a differential equation. One such approach would be the spectral methods. An advantage of a spectral collocation method is that it gives high accurate solutions with relatively fewer spatial grid nodes when compared with other numerical techniques. In [27], the Jacobi rational collocation scheme was proposed and developed to solve generalized pantograph equations. In [28], the authors extended the application of Jacobi-Gauss-Lobatto collocation approximation to solve nonlinear Schrödinger equations. Also, the generalized Laguerre-Legendre collocation method has been successfully applied to initial-boundary value problems [29]. In [30], approximate solutions of nonlinear Klein-Gordon and Sine-Gordon equations were provided using the Chebyshev tau meshless scheme. For some recent developments on spectral methods, see [31–34].

Neutral functional-differential equations play an important role in the mathematical modeling of several phenomena. It is well known that most of delay differential equations cannot be solved exactly. Therefore, numerical methods would be presented and developed to get approximate solutions of these equations. In this direction, Ishiwata and Muroya [35] applied the rational approximation scheme for solving a class of delay differential equations. In [36], Chen and Wang implemented the variational iteration scheme to obtain an analytical solution of the neutral functional-differential equation. Very recently, Heydari et al. [37] proposed a new numerical algorithm based on the operational matrix formulation of Chebyshev cardinal functions for solving delay differential equations arising in electrodynamics. In this paper we propose a numerical solution for a new class of delay differential equations, namely, fractional neutral functional-differential equations (FNFDEs) with proportional delay.

The main aim of this paper is to design a suitable way to approximate a new class of functional-differential equations with fractional orders on the interval using spectral collocation method. The spectral shifted Jacobi-Gauss collocation (SJGC) approximation is proposed to obtain the numerical solution . The SJGC approximation, which is more reliable, is employed to obtain approximate solution of FNFDEs with leading fractional order and initial conditions. We choose the nodes of the shifted Jacobi-Gauss interpolation on as suitable collocation nodes. The Legendre and Chebyshev collocation approximations can be obtained as special cases from our general approach. Finally, the validity and effectiveness of the method are demonstrated by solving two numerical examples. Numerical examples are presented in the form of tables and graphs to make comparisons with the results obtained by other methods and with the exact solutions more easier.

In the next section, we present an overview of shifted Jacobi polynomials and fractional calculus needed hereafter. Section 3 is devoted to present and implement the collocation scheme for solving FNFDEs with proportional delay using Jacobi polynomials. In Section 4, we introduce two numerical examples demonstrating the high accuracy and efficiency of the present numerical algorithm.

2. Preliminaries

Here, we state some preliminaries of fractional calculus [38] and some relevant properties of Jacobi polynomials. The most commonly used definition of fractional integral is the Riemann-Liouville operator.

Definition 1. The Riemann-Liouville fractional integral operator of order is defined as

Definition 2. The Caputo fractional derivatives of order are defined as

where is the classical differential operator of order .

Also

where and , while and are the floor and ceiling functions, respectively,

The Caputo’s fractional differentiation is a linear operation, similar to the integer-order differentiation

where and are constants.

Let and let be the Jacobi polynomial of degree ; then we get

Besides,

Let ; then we define the weighted space as usual, equipped with the following inner product and norm:

The set of Jacobi polynomials forms a complete -orthogonal system, and

Let us define the shifted Jacobi polynomial of degree by , , and thanks to (6) and (7), yield

The shifted Jacobi polynomials of degree on the interval are given by

Since the analytic form of is given by (12), with the use of (4), (5), and (12), we obtain

Now, approximate by terms of shifted Jacobi series; we have

where is given from (18) with , and

Employing (13)–(15) we get

where , and

A function may be expressed in terms of shifted Jacobi polynomials as

In practice, only the first terms shifted Jacobi polynomials are considered. Then we have

Next, let ; then we define the weighted space in the usual way, with the following inner product and norm:

The set of shifted Jacobi polynomials forms a complete -orthogonal system. Moreover, and due to (9), we have

3. Shifted Jacobi Collocation Approximation for FNFDEs

In this section, we propose the shifted Jacobi collocation method with the Jacobi-Gauss quadrature nodes to solve numerically the following FNFDEs with proportional delay:

with the initial conditions

Here, and () are given analytical functions, , and denote given constants with . By using the shifted Jacobi-Gauss collocation method [39], we can approximate the fractional neutral functional-differential equations with proportional delays, without any artificial boundary and variable transformation. Let us first introduce some basic notation that will be used in the sequel.

Now we introduce the Jacobi-Gauss-Lobatto quadratures in two different intervals , and . Denoting by , , and , , the nodes and Christoffel numbers of the standard (shifted) Jacobi-Gauss-Lobatto quadratures on , respectively. Therefore, we can deduce that

Let be the set of all polynomials of degree ; then, for any , we have

We set

and the inner product and norm are defined as

Obviously,

Thus, for any , the norms and coincide.

Associating with this quadrature rule, we denote by the Jacobi-Gauss interpolation operator

The shifted Jacobi-Gauss collocation method for solving (22) and (23) is to seek , such that

We now derive the algorithm for solving (22) and (23). To do this, let

We first approximate and , , using (31). By substituting this approximation in (22), we get

Making use of (16), we deduce that

Also, by substituting (31) in (23) we obtain

Now, we collocate (33) at the shifted Jacobi-Gauss interpolation points, yielding

Next (34), after using (10), can be written as

Finally, relations (35) and (36) generate set of algebraic equations which can be solved for the unknown coefficients , , by using any standard solver technique.

4. Numerical Results

In this section, two fractional neutral functional-differential equations with proportional delays are solved by the SJGC method. We implement the method presented in this paper for these two examples to demonstrate the accuracy and capability of the proposed algorithm.

Example 1. Consider the following FNFDEs with proportional delay:

where

and the exact solution is given by .

Table 1 lists the results obtained by the shifted Jacobi collocation method in terms of absolute errors at with (first kind shifted Chebyshev collocation method), (shifted Legendre collocation method), and (second kind shifted Chebyshev collocation method). In the case of , the approximate solution by the presented method is shown in Figure 1, to make it easier to compare with the analytic solution.

Table 1: Absolute errors using SJGC method at for Example 1.

0.5

1.0

1.5

2.0

2.5

3.0

3.5

4.0

4.5

5.0

Figure 1: Graph of exact solution and approximate solution for at for Example 1.

Example 2. Consider the following FNFDEs with proportional delay:

subject to

with exact solution .

In Table 2, we list the absolute errors obtained by the shifted Jacobi collocation method, with several values of and at . It is clear that, for all Jacobi polynomials parameters, the results are stable. Meanwhile, Figure 2 presents the SJGC solution with at and exact solution, which are found to be in excellent agreement.

Table 2: Absolute errors using SJGC method at for Example 2.

0.1

0.2

0.3

0.4

0.5

0.6

0.7

0.8

0.9

1.0

Figure 2: Graph of exact solution and approximate solution for at for Example 2.

5. Conclusion

In this paper, we have proposed a numerical algorithm to solve a class of fractional delay differential equations. The Jacobi collocation approximation was developed to solve this problem. A number of collocation techniques can be obtained as special cases from the proposed technique. Numerical results were given to demonstrate the accuracy and applicability of the presented method.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgment

This paper was funded by the Deanship of Scientific Research DSR, King Abdulaziz University, Jeddah. The authors, therefore, acknowledge and express their thanks to DSR technical and financial support.

References

F. Flandoli and C. A. Tudor, “Brownian and fractional Brownian stochastic currents via Malliavin calculus,” Journal of Functional Analysis, vol. 258, no. 1, pp. 279–306, 2010.

R. L. Magin, C. Ingo, L. Colon-Perez, W. Triplett, and T. H. Mareci, “Characterization of anomalous diffusion in porous biological tissues using fractional order derivatives and entropy,” Microporous and Mesoporous Materials, vol. 178, pp. 39–43, 2013.

A. H. Bhrawy and M. A. Alghamdi, “A shifted Jacobi-Gauss-Lobatto collocation method for solving nonlinear fractional Langevin equation involving two fractional orders in different intervals,” Boundary Value Problems, vol. 2012, article 62, 2012.

M. D. Choudhury, S. Chandra, S. Nag, S. Das, and S. Tarafdar, “Forced spreading and rheology of starch gel: viscoelastic modeling with fractional calculus,” Colloids and Surfaces A, vol. 407, pp. 64–70, 2012.

A. Atangana and N. Bildik, “The use of fractional order derivative to predict the groundwater flow,” Mathematical Problems in Engineering, vol. 2013, Article ID 543026, 9 pages, 2013.

A. Ghomashi, S. Salahshour, and A. Hakimzadeh, “Approximating solutions of fully fuzzy linear systems: a financial case study,” Journal of Intelligent and Fuzzy Systems, vol. 26, pp. 367–378, 2014.

F. Gao, X.-J. Lee, H. Tong, F. Fei, and H. Zhao, “Identification of unknown parameters and orders via cuckoo search oriented statistically by differential evolution for noncommensurate fractional-order chaotic systems,” Abstract and Applied Analysis, vol. 2013, Article ID 382834, 19 pages, 2013.

F. Gao, X. Lee, F. Fei, H. Tong, Y. Deng, and H. Zhao, “Identification time-delayed fractional order chaos with functional extrema model via differential evolution,” Expert Systems with Applications, vol. 41, pp. 1601–1608, 2014.

J. Deng and L. Ma, “Existence and uniqueness of solutions of initial value problems for nonlinear fractional differential equations,” Applied Mathematics Letters, vol. 23, no. 6, pp. 676–680, 2010.

Y. Zhao, D.-F. Cheng, and X.-J. Yang, “Approximation solutions for local fractional Schrödinger equation in the one-dimensional Cantorian system,” Advances in Mathematical Physics, vol. 2013, Article ID 291386, 5 pages, 2013.

Y. Zhao, D. Baleanu, C. Cattani, D.-F. Cheng, and X.-J. Yang, “Maxwell's equations on Cantor sets: a local fractional approach,” Advances in High Energy Physics, vol. 2013, Article ID 686371, 6 pages, 2013.

A.-M. Yang, Y.-Z. Zhang, C. Cattani et al., “Application of local fractional series expansion method to solve Klein-Gordon equations on Cantor sets,” Abstract and Applied Analysis, vol. 2014, Article ID 372741, 6 pages, 2014.

B. Jacobs and C. Harley, “Two hybrid methods for solving two-dimensional linear time-fractional partial differential equations,” Abstract and Applied Analysis. In press

Y. Z. Zhang, A. M. Yang, and X.-J. Yang, “1-D heat conduction in a fractal medium: a solution by the local fractional Fourier series method,” Thermal Science, vol. 17, no. 3, pp. 953–956, 2013.

A. Atangana and S. B. Belhaouari, “Solving partial differential equation with space- and time-fractional derivatives via homotopy decomposition method,” Mathematical Problems in Engineering, vol. 2013, Article ID 318590, 9 pages, 2013.

M. P. Tripathi, V. K. Baranwal, R. K. Pandey, and O. P. Singh, “A new numerical algorithm to solve fractional differential equations based on operational matrix of generalized hat functions,” Communications in Nonlinear Science and Numerical Simulation, vol. 18, no. 6, pp. 1327–1340, 2013.

A. H. Bhrawy and M. M. Al-Shomrani, “A shifted Legendre spectral method for fractional-order multi-point boundary value problems,” Advances in Difference Equations, vol. 2012, article 8, 2012.

Ş. Yüzbaşı, “Numerical solutions of fractional Riccati type differential equations by means of the Bernstein polynomials,” Applied Mathematics and Computation, vol. 219, no. 11, pp. 6328–6343, 2013.

A. Pedas and E. Tamme, “Numerical solution of nonlinear fractional differential equations by spline collocation methods,” Journal of Computational and Applied Mathematics, vol. 255, pp. 216–230, 2014.

A. H. Bhrawy, “A new numerical algorithm for solving a class of fractional advection-dispersion equation with variable coefficients using Jacobi polynomials,” Abstract and Applied Analysis, vol. 2013, Article ID 954983, 9 pages, 2013.

D. Baleanu, A. H. Bhrawy, and T. M. Taha, “Two efficient generalized Laguerre spectral algorithms for fractional initial value problems,” Abstract and Applied Analysis, vol. 2013, Article ID 546502, 10 pages, 2013.

X. Ma and C. Huang, “Spectral collocation method for linear fractional integro-differential equations,” Applied Mathematical Modelling, vol. 38, no. 4, pp. 1434–1448, 2014.

Y. Yang and Y. Huang, “Spectral-collocation methods for fractional pantograph delay-integrodifferential equations,” Advances in Mathematical Physics, vol. 2013, Article ID 821327, 14 pages, 2013.

F. Yin, J. Song, Y. Wu, and L. Zhang, “Numerical solution of the fractional partial differential equations by the two-dimensional fractional-order Legendre functions,” Abstract and Applied Analysis, vol. 2013, Article ID 562140, 13 pages, 2013.

F. Yin, J. Song, H. Leng, and F. Lu, “Couple of the variational iteration method and fractional-order legendre functions method for fractional differential equations,” The Scientific World Journal, vol. 2014, Article ID 928765, 9 pages, 2014.

M. Kamrani, “Numerical solution of stochastic fractional differential equations,” Numerical Algorithms, 2014.

E. H. Doha, A. H. Bhrawy, D. Baleanu, and R. M. Hafez, “A new Jacobi rational-Gauss collocation method for numerical solution of generalized pantograph equations,” Applied Numerical Mathematics, vol. 77, pp. 43–54, 2014.

E. H. Doha, A. H. Bhrawy, M. A. Abdelkawy, and R. A. van Gorder, “Jacobi-Gauss-Lobatto collocation method for the numerical solution of nonlinear Schrödinger equations,” Journal of Computational Physics, vol. 261, pp. 244–255, 2014.

M. Tatari and M. Haghighi, “A generalized Laguerre-Legendre spectral collocation method for solving initial-boundary value problems,” Applied Mathematical Modelling, vol. 38, no. 4, pp. 1351–1364, 2014.

W. Shao and X. Wu, “The numerical solution of the nonlinear Klein-Gordon and Sine-Gordon equations using the Chebyshev tau meshless method,” Computer Physics Communications, vol. 185, no. 5, pp. 1399–1409, 2014.

E. H. Doha and A. H. Bhrawy, “An efficient direct solver for multidimensional elliptic Robin boundary value problems using a Legendre spectral-Galerkin method,” Computers and Mathematics with Applications, vol. 64, no. 4, pp. 558–571, 2012.

E. H. Doha, A. H. Bhrawy, and R. M. Hafez, “On shifted Jacobi spectral method for high-order multi-point boundary value problems,” Communications in Nonlinear Science and Numerical Simulation, vol. 17, no. 10, pp. 3802–3810, 2012.

A. H. Bhrawy, E. Tohidi, and F. Soleymani, “A new Bernoulli matrix method for solving high-order linear and nonlinear Fredholm integro-differential equations with piecewise intervals,” Applied Mathematics and Computation, vol. 219, no. 2, pp. 482–497, 2012.

A. H. Bhrawy, M. M. Alghamdi, and T. M. Taha, “A new modified generalized Laguerre operational matrix of fractional integration for solving fractional differential equations on the half line,” Advances in Difference Equations, vol. 2012, article 179, 2012.

E. Ishiwata and Y. Muroya, “Rational approximation method for delay differential equations with proportional delay,” Applied Mathematics and Computation, vol. 187, no. 2, pp. 741–747, 2007.

X. Chen and L. Wang, “The variational iteration method for solving a neutral functional-differential equation with proportional delays,” Computers and Mathematics with Applications, vol. 59, no. 8, pp. 2696–2702, 2010.

M. Heydari, G. B. Loghmani, and S. M. Hosseini, “Operational matrices of Chebyshev cardinal functions and their application for solving delay differential equations arising in electrodynamics with error estimation,” Applied Mathematical Modelling, vol. 37, no. 14-15, pp. 7789–7809, 2013.

I. Podlubny, Fractional Differential Equations, vol. 198 of Mathematics in Science and Engineering, Academic Press, San Diego, Calif, USA, 1999.

A. H. Bhrawy, “A Jacobi-Gauss-Lobatto collocation method for solving generalized Fitzhugh-Nagumo equation with time-dependent coefficients,” Applied Mathematics and Computation, vol. 222, pp. 255–264, 2013.