Research Article  Open Access
A Shifted JacobiGauss Collocation Scheme for Solving Fractional Neutral FunctionalDifferential Equations
Abstract
The shifted JacobiGauss collocation (SJGC) scheme is proposed and implemented to solve the fractional neutral functionaldifferential 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 JacobiGauss scheme is to reduce solving the generalized fractional neutral functionaldifferential equations to a system of algebraic equations in the unknown expansion. Reasonable numerical results are achieved by choosing few shifted JacobiGauss 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 timedependent 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 fractionalorder 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 JacobiGaussLobatto collocation approximation to solve nonlinear Schrödinger equations. Also, the generalized LaguerreLegendre collocation method has been successfully applied to initialboundary value problems [29]. In [30], approximate solutions of nonlinear KleinGordon and SineGordon equations were provided using the Chebyshev tau meshless scheme. For some recent developments on spectral methods, see [31–34].
Neutral functionaldifferential 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 functionaldifferential 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 functionaldifferential equations (FNFDEs) with proportional delay.
The main aim of this paper is to design a suitable way to approximate a new class of functionaldifferential equations with fractional orders on the interval using spectral collocation method. The spectral shifted JacobiGauss 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 JacobiGauss 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 RiemannLiouville operator.
Definition 1. The RiemannLiouville 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 integerorder 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 JacobiGauss 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 JacobiGauss collocation method [39], we can approximate the fractional neutral functionaldifferential 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 JacobiGaussLobatto quadratures in two different intervals , and . Denoting by , , and , , the nodes and Christoffel numbers of the standard (shifted) JacobiGaussLobatto 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 JacobiGauss interpolation operator
The shifted JacobiGauss 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 JacobiGauss 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 functionaldifferential 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.

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.

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. View at: Publisher Site  Google Scholar  MathSciNet
 R. L. Magin, C. Ingo, L. ColonPerez, 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. View at: Google Scholar
 A. H. Bhrawy and M. A. Alghamdi, “A shifted JacobiGaussLobatto collocation method for solving nonlinear fractional Langevin equation involving two fractional orders in different intervals,” Boundary Value Problems, vol. 2012, article 62, 2012. View at: Google Scholar
 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. View at: Publisher Site  Google Scholar
 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. View at: Publisher Site  Google Scholar  MathSciNet
 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. View at: Google Scholar
 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 fractionalorder chaotic systems,” Abstract and Applied Analysis, vol. 2013, Article ID 382834, 19 pages, 2013. View at: Publisher Site  Google Scholar  MathSciNet
 F. Gao, X. Lee, F. Fei, H. Tong, Y. Deng, and H. Zhao, “Identification timedelayed fractional order chaos with functional extrema model via differential evolution,” Expert Systems with Applications, vol. 41, pp. 1601–1608, 2014. View at: Google Scholar
 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. View at: Publisher Site  Google Scholar
 Y. Zhao, D.F. Cheng, and X.J. Yang, “Approximation solutions for local fractional Schrödinger equation in the onedimensional Cantorian system,” Advances in Mathematical Physics, vol. 2013, Article ID 291386, 5 pages, 2013. View at: Publisher Site  Google Scholar  MathSciNet
 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. View at: Publisher Site  Google Scholar  MathSciNet
 A.M. Yang, Y.Z. Zhang, C. Cattani et al., “Application of local fractional series expansion method to solve KleinGordon equations on Cantor sets,” Abstract and Applied Analysis, vol. 2014, Article ID 372741, 6 pages, 2014. View at: Publisher Site  Google Scholar  MathSciNet
 B. Jacobs and C. Harley, “Two hybrid methods for solving twodimensional linear timefractional partial differential equations,” Abstract and Applied Analysis. In press. View at: Google Scholar
 Y. Z. Zhang, A. M. Yang, and X.J. Yang, “1D 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. View at: Google Scholar
 A. Atangana and S. B. Belhaouari, “Solving partial differential equation with space and timefractional derivatives via homotopy decomposition method,” Mathematical Problems in Engineering, vol. 2013, Article ID 318590, 9 pages, 2013. View at: Publisher Site  Google Scholar  MathSciNet
 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. View at: Publisher Site  Google Scholar  MathSciNet
 A. H. Bhrawy and M. M. AlShomrani, “A shifted Legendre spectral method for fractionalorder multipoint boundary value problems,” Advances in Difference Equations, vol. 2012, article 8, 2012. View at: Google Scholar
 Ş. 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. View at: Publisher Site  Google Scholar  MathSciNet
 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. View at: Publisher Site  Google Scholar  MathSciNet
 A. H. Bhrawy, “A new numerical algorithm for solving a class of fractional advectiondispersion equation with variable coefficients using Jacobi polynomials,” Abstract and Applied Analysis, vol. 2013, Article ID 954983, 9 pages, 2013. View at: Publisher Site  Google Scholar  MathSciNet
 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. View at: Publisher Site  Google Scholar  MathSciNet
 X. Ma and C. Huang, “Spectral collocation method for linear fractional integrodifferential equations,” Applied Mathematical Modelling, vol. 38, no. 4, pp. 1434–1448, 2014. View at: Publisher Site  Google Scholar  MathSciNet
 Y. Yang and Y. Huang, “Spectralcollocation methods for fractional pantograph delayintegrodifferential equations,” Advances in Mathematical Physics, vol. 2013, Article ID 821327, 14 pages, 2013. View at: Publisher Site  Google Scholar  MathSciNet
 F. Yin, J. Song, Y. Wu, and L. Zhang, “Numerical solution of the fractional partial differential equations by the twodimensional fractionalorder Legendre functions,” Abstract and Applied Analysis, vol. 2013, Article ID 562140, 13 pages, 2013. View at: Publisher Site  Google Scholar  MathSciNet
 F. Yin, J. Song, H. Leng, and F. Lu, “Couple of the variational iteration method and fractionalorder legendre functions method for fractional differential equations,” The Scientific World Journal, vol. 2014, Article ID 928765, 9 pages, 2014. View at: Publisher Site  Google Scholar
 M. Kamrani, “Numerical solution of stochastic fractional differential equations,” Numerical Algorithms, 2014. View at: Publisher Site  Google Scholar
 E. H. Doha, A. H. Bhrawy, D. Baleanu, and R. M. Hafez, “A new Jacobi rationalGauss collocation method for numerical solution of generalized pantograph equations,” Applied Numerical Mathematics, vol. 77, pp. 43–54, 2014. View at: Publisher Site  Google Scholar  MathSciNet
 E. H. Doha, A. H. Bhrawy, M. A. Abdelkawy, and R. A. van Gorder, “JacobiGaussLobatto collocation method for the numerical solution of $1+1$ nonlinear Schrödinger equations,” Journal of Computational Physics, vol. 261, pp. 244–255, 2014. View at: Publisher Site  Google Scholar  MathSciNet
 M. Tatari and M. Haghighi, “A generalized LaguerreLegendre spectral collocation method for solving initialboundary value problems,” Applied Mathematical Modelling, vol. 38, no. 4, pp. 1351–1364, 2014. View at: Publisher Site  Google Scholar  MathSciNet
 W. Shao and X. Wu, “The numerical solution of the nonlinear KleinGordon and SineGordon equations using the Chebyshev tau meshless method,” Computer Physics Communications, vol. 185, no. 5, pp. 1399–1409, 2014. View at: Publisher Site  Google Scholar  MathSciNet
 E. H. Doha and A. H. Bhrawy, “An efficient direct solver for multidimensional elliptic Robin boundary value problems using a Legendre spectralGalerkin method,” Computers and Mathematics with Applications, vol. 64, no. 4, pp. 558–571, 2012. View at: Publisher Site  Google Scholar  MathSciNet
 E. H. Doha, A. H. Bhrawy, and R. M. Hafez, “On shifted Jacobi spectral method for highorder multipoint boundary value problems,” Communications in Nonlinear Science and Numerical Simulation, vol. 17, no. 10, pp. 3802–3810, 2012. View at: Publisher Site  Google Scholar  MathSciNet
 A. H. Bhrawy, E. Tohidi, and F. Soleymani, “A new Bernoulli matrix method for solving highorder linear and nonlinear Fredholm integrodifferential equations with piecewise intervals,” Applied Mathematics and Computation, vol. 219, no. 2, pp. 482–497, 2012. View at: Publisher Site  Google Scholar  MathSciNet
 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. View at: Google Scholar
 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. View at: Publisher Site  Google Scholar  MathSciNet
 X. Chen and L. Wang, “The variational iteration method for solving a neutral functionaldifferential equation with proportional delays,” Computers and Mathematics with Applications, vol. 59, no. 8, pp. 2696–2702, 2010. View at: Publisher Site  Google Scholar  MathSciNet
 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. 1415, pp. 7789–7809, 2013. View at: Publisher Site  Google Scholar  MathSciNet
 I. Podlubny, Fractional Differential Equations, vol. 198 of Mathematics in Science and Engineering, Academic Press, San Diego, Calif, USA, 1999. View at: MathSciNet
 A. H. Bhrawy, “A JacobiGaussLobatto collocation method for solving generalized FitzhughNagumo equation with timedependent coefficients,” Applied Mathematics and Computation, vol. 222, pp. 255–264, 2013. View at: Publisher Site  Google Scholar  MathSciNet
Copyright
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.