Abstract and Applied Analysis

Volume 2013, Article ID 513808, 9 pages

http://dx.doi.org/10.1155/2013/513808

## Numerical Solution of a Class of Functional-Differential Equations Using Jacobi Pseudospectral Method

^{1}Department of Mathematics, Faculty of Science, King Abdulaziz University, Jeddah 21589, Saudi Arabia^{2}Department of Mathematics, Faculty of Science, Beni-Suef University, Beni-Suef 62511, Egypt^{3}Department of Chemical and Materials Engineering, Faculty of Engineering, King Abdulaziz University, Jeddah 21589, Saudi Arabia^{4}Department of Mathematics and Computer Sciences, Cankaya University, Eskisehir Yolu 29.km, 06810 Ankara, Turkey^{5}Institute of Space Sciences, P.O. Box MG-23, 76900 Magurele-Bucharest, Romania

Received 24 August 2013; Accepted 18 September 2013

Academic Editor: Soheil Salahshour

Copyright © 2013 A. H. Bhrawy 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

The shifted Jacobi-Gauss-Lobatto pseudospectral (SJGLP) method is applied to neutral functional-differential equations (NFDEs) with proportional delays. The proposed approximation is based on shifted Jacobi collocation approximation with the nodes of Gauss-Lobatto quadrature. The shifted Legendre-Gauss-Lobatto Pseudo-spectral and Chebyshev-Gauss-Lobatto Pseudo-spectral methods can be obtained as special cases of the underlying method. Moreover, the SJGLP method is extended to numerically approximate the nonlinear high-order NFDE with proportional delay. Some examples are displayed for implicit and explicit forms of NFDEs to demonstrate the computation accuracy of the proposed method. We also compare the performance of the method with variational iteration method, one-leg -method, continuous Runge-Kutta method, and reproducing kernel Hilbert space method.

#### 1. Introduction

In the last four decades, spectral method has become increasingly popular and been successfully applied in solving all types of differential equations owing to its high order of accuracy (see, for instance, [1–5]). Recently there has been a growing interest in applying spectral methods for the numerical solution of fractional differential equations and delay differential equations (DDEs). Moreover, the principal difficulty in studying NFDEs with proportional delays lies in their special transcendental nature. Thus we propose an efficient technique to solve such differential equations numerically based on shifted Jacobi polynomials.

In the last two decades, some numerical approaches for treating several types of DDEs were presented in [6–12]. In this direction, Zhao et al. [13] studied the stability of Runge-Kutta approach and applied it for solving the NFDEs with proportional delays. Tohidi et al. [12] derived the operational matrix of Bernoulli polynomial and used it with collocation method to solve the pantograph type equation. The one-leg approach was implemented to solve nonlinear NFDEs in [14]. Meanwhile, Trif [15] proposed a direct solution technique by using the Chebyshev Tau operational matrix method for solving the pantograph type equation. More recently, Bhrawy et al. [16] proposed the Legendre pseudospectral algorithm with studying the error analysis and stability of the proposed algorithm for a class of DDEs. Sun and Zhang [17] proposed a compact difference method for solving nonlinear partial DDE. Cordero and Escalante [18] extended the application of segmented tau approach for solving a class of NFDEs, meanwhile, the history-valued NFDEs were approximated in [19]. Isik et al. [20] introduced a new Bernstein collocation method for the numerical solution of DDEs of pantograph-type with retarded case. The existence of solutions of neutral functional-differential equations with proportional delays had been discussed in [21–24].

The aim of this paper is to develop a direct solution technique to approximate the linear high-order NFDEs with proportional delays using the shifted Jacobi polynomials on the interval ; we dedicate the shifted Jacobi-Gauss-Lobatto pseudospectral (SJGLP) method to find the approximate solution . Approximate semianalytical solution with high accuracy can be obtained by selecting a limited number of Gauss-Lobatto collocation points for the linear and nonlinear high-order NFDEs with proportional delays. For suitable collocation points we use the nodes of the shifted Jacobi-Gauss-Lobatto interpolation on in which the nodes are distinct and lie between and , and the two endpoints of the domain ( and ) are used as the first and last collocation points, respectively. The algorithms given in [16] can be obtained as special cases from the proposed algorithms. Finally, the accuracy of the proposed method is showed by test problems. From the results, these algorithms are extremely efficient and accurate for solving NFDEs.

The paper is arranged in the following way. In the next section, some basic properties of Jacobi polynomials which are required in the present paper are given, and in Section 3, the way of constructing the pseudospectral technique for NFDEs with proportional delays is described using the shifted Jacobi polynomials. In Section 4, we investigate the shifted Jacobi-Gauss-Lobatto pseudo-spectral (SJGLP) method for solving nonlinear high-order NFDEs with proportional delays. Some numerical results exhibiting the accuracy of the proposed algorithm are given in Section 5. Finally, a conclusion is given in Section 6.

#### 2. Preliminaries

In this section, we briefly recall some properties of the Jacobi polynomials (, , , ), which satisfy the following relations: The th derivative of Jacobi polynomials of degree can be given by These polynomials are the only polynomials arising as eigenfunctions of the following singular Sturm-Liouville equation: Let ; then we define the weighted space as usual. The inner product and the norm of with respect to the weight function are defined as follows:

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

Let us denote , . By the shifted Jacobi polynomial of degree and by using (1) and (2); then we deduce that

The following inner product and norm are defined on the weighted space subject to the weight function . Moreover, forms a complete -orthogonal system.

According to (5), we get

#### 3. Linear High-Order NFDE with Proportional Delay

In this section, we shall investigate solutions to NFDEs with proportional delays of the form with the initial conditions where and () are given functions, meanwhile, , , , and are constants with (). This class of equations plays an important role in modeling phenomena of the real world.

In the pseudo-spectral methods [25–28], one needs to exactly satisfy the differential equation at specified collocation points in the domain of solution. Generally, the distribution of the collocation nodes can be freely chosen, but an accurate approximations are obtained by selecting the collocation nodes as the zeros of the orthogonal polynomials. For shifted Jacobi polynomials, two commonly used quadrature and collocation nodes, namely, (i) shifted Jacobi-Gauss nodes (in the interior of the domain) and (ii) shifted Jacobi-Gauss-Lobatto nodes (in the interior and at the two endpoints of the domain).

Now, we will present the shifted Jacobi-Gauss-Lobatto type quadratures. Let , , be the nodes of the Jacobi-Gauss-Lobatto interpolation on , and let , , be the corresponding weights. Throughout this paper, we assume that , , stands for the nodes of the shifted Jacobi-Gauss-Lobatto interpolation on the interval . Thus , , and their corresponding weights are , . Let be the set of all polynomials of degree =. One gets for any ,

Let us first introduce some basic notations that will be used in the sequel. We set The discrete inner product and norm are defined by where and are the nodes and the corresponding weights of the shifted Jacobi-Gauss-Lobatto quadrature formula on the interval , respectively. Obviously,

The shifted Jacobi-Gauss-Lobatto pseudo-spectral method for solving (28) and (29) is to seek , such that where the , , are distinct and lie between and , , and . For simplicity in presentation and without loss of generality, assume that . We now derive the collocation algorithm for solving (28) and (29). To do this, consider the solution is approximated by a truncated Jacobi expansion

We first approximate and as (32). By substituting these approximation in (28), we get

Then, by virtue of (8), we deduce that

Also, by substituting (32) in (29), we obtain

To find the solution , we first collocate (20) at the shifted Jacobi roots that yields

Next, (21), after using (7), can be written as

Thus (22) with relation (23) can be written as a linear algebraic system. To do this, let us consider

The matrix system associated with (22) and (23) becomes where the matrices , , , , , and are given explicitly in the following. If we denote that , , , , , and , then the elements , , , and are given by Moreover, the elements of the matrix corresponding to the mixed initial conditions are given by

In the case of , , , and , the linear system (25) can be solved by forming explicitly the LU factorization; that is, . The expense of calculating LU factorization is operations, and the expense of solving the linear system (25), provided that the factorization is known, is .

#### 4. Nonlinear High-Order NFDE with Proportional Delay

In this section, we investigate the shifted Jacobi-Gauss-Lobatto pseudospectral method to numerically approximate the nonlinear high-order NFDE with proportional delay; namely, subject to where , , and are constants with (), meanwhile, is nonlinear in general. This equation is a generalized form of the proportional delay differential equations given in [14, 29, 30] and plays an important role in modeling phenomena of the real world.

The shifted Jacobi-Gauss-Lobatto pseudospectral approximation for (28) is to find such that

Now, we approximate the numerical solution as a truncated series expansion of shifted Jacobi polynomial in the form Accordingly, (30) can be written as

Next, making use of relation (8), thus the high-order derivatives of the proposed solution with proportional delays, , can be expressed explicitly in terms of the shifted Jacobi polynomials and the expansion coefficients ; consequently

Finally, to find the unknown expansion coefficients , , we implement any iteration technique to solve a system of nonlinear algebraic equations resulting from the combination of nonlinear algebraic equations resulting from (33) and linear algebraic equations resulting from (29).

#### 5. Numerical Results

In this section, we will carry out three test examples to study the validity and effectiveness of the proposed method and also show that high accurate solutions are achieved using a few number of the Jacobi Gauss-Lobatto points. Moreover, comparisons with other methods reveal that the present method is accurate and convenient. All the numerical computations have been performed by the symbolic computation software Mathematica 8.0.

*Example 1. *Consider the first-order NFDE with proportional delay considered in [31]
which has the exact solution .

Table 1 lists the absolute error using Jacobi pseudospectral method for three choices of at in the interval . We compare the errors obtained by the proposed method with variational iteration (VI) method [32], the one-leg method [14, 33] with , and RKHSM method [31]. The graph of analytical solution and approximate solution in long interval for and , is displayed in Figure 1 to make it easier to compare with analytical solution. Moreover, in this case the graph of the error is given in Figure 2. Consequently, we conclude that the approximate solution by Jacobi pseudospectral method agree very well with the exact solution.

*Example 2. *Let us consider the second-order NFDE with proportional delay
which enjoys exact solution .

In Table 2, we introduce the absolute error using the proposed method at with various choices of and . The resulting graph of (35) for the presented method in the case of at and the analytic solution are shown in Figure 3.

*Example 3. *Consider the third-order NFDE with proportional delays
where
The exact solution of the problem is .

Table 3 lists the maximum absolute errors using SJGLP method in the interval at and different choices of and .

In case of Chebyshev polynomials of the second kind (), the graph of exact solution and approximate solution for is plotted in Figure 4 in the interval .

#### 6. Conclusion

In this paper, we have demonstrated the feasibility of SJGLP for solving linear NFDEs with proportional delays. We also have discussed the resulting linear system. Moreover, we have implemented the SJGLP method to numerically approximate the nonlinear high-order NFDE with proportional delay.

All the given examples reveal that the results of SJGLP method are in excellent agreement with the analytical solutions. It is concluded from the aforementioned tables and figures that SJGLP method is an accurate and efficient method to solve NFDEs when compared with those generated by some other methods.

In the future work, we address the Jacobi pseudospectral approximation for the solution of linear and nonlinear delay partial differential equations in two and three dimensions (see, e.g., [17, 34]).

#### References

- C. Canuto, M. Y. Hussaini, A. Quarteroni, and T. A. Zang,
*Spectral Methods in Fluid Dynamics*, Springer Series in Computational Physics, Springer, New York, NY, USA, 1988. View at MathSciNet - L. N. Trefethen,
*Spectral Methods in MATLAB*, vol. 10 of*Software, Environments, and Tools*, SIAM, Philadelphia, Pa, USA, 2000. View at Publisher · View at Google Scholar · View at MathSciNet - E. H. Doha, A. H. Bhrawy, and S. S. Ezz-Eldien, “A new Jacobi operational matrix: an application for solving fractional differential equations,”
*Applied Mathematical Modelling*, vol. 36, no. 10, pp. 4931–4943, 2012. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - A. H. Bhrawy and M. Al-shomrani, “A shifted Legendre spectral method for fractional order multi-point boundary value problems,”
*Advances in Difference Equations*, vol. 2012, article 8, 8 pages, 2012. View at Publisher · View at Google Scholar - A. H. Bhrawy, M. A. 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, 12 pages, 2012. View at Publisher · View at Google Scholar - A. Bellen and M. Zennaro,
*Numerical Methods for Delay Differential Equations*, Numerical Mathematics and Scientific Computation, Oxford University Press, New York, NY, USA, 2003. View at Publisher · View at Google Scholar · View at MathSciNet - M. Z. Liu and D. Li, “Properties of analytic solution and numerical solution of multi-pantograph equation,”
*Applied Mathematics and Computation*, vol. 155, no. 3, pp. 853–871, 2004. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - G. A. Bocharov and F. A. Rihan, “Numerical modelling in biosciences using delay differential equations,”
*Journal of Computational and Applied Mathematics*, vol. 125, no. 1-2, pp. 183–199, 2000. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - M. A. Jafari and A. Aminataei, “Method of successive approximations for solving the multi-pantograph delay equations,”
*General Mathematics Notes*, vol. 8, pp. 23–28, 2012. View at Google Scholar - S. Yüzbaşı, “An efficient algorithm for solving multi-pantograph equation systems,”
*Computers & Mathematics with Applications*, vol. 64, no. 4, pp. 589–603, 2012. View at Publisher · View at Google Scholar · View at MathSciNet - S. Yüzbaşi, N. Şahin, and M. Sezer, “A Bessel collocation method for numerical solution of generalized pantograph equations,”
*Numerical Methods for Partial Differential Equations*, vol. 28, no. 4, pp. 1105–1123, 2012. View at Publisher · View at Google Scholar · View at MathSciNet - E. Tohidi, A. H. Bhrawy, and K. Erfani, “A collocation method based on Bernoulli operational matrix for numerical solution of generalized pantograph equation,”
*Applied Mathematical Modelling*, vol. 37, no. 6, pp. 4283–4294, 2013. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - 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 · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - W.-S. Wang and S.-F. Li, “On the one-leg $\theta $-methods for solving nonlinear neutral functional differential equations,”
*Applied Mathematics and Computation*, vol. 193, no. 1, pp. 285–301, 2007. View at Publisher · View at Google Scholar · View at MathSciNet - D. Trif, “Direct operatorial tau method for pantograph-type equations,”
*Applied Mathematics and Computation*, vol. 219, no. 4, pp. 2194–2203, 2012. View at Publisher · View at Google Scholar · View at MathSciNet - A. H. Bhrawy, L. M. Assas, E. Tohidi, and M. A. Alghamdi, “Legendre-Gauss collocation method for neutral functional-differential equations with proportional delays,”
*Advances in Differential Equations*, vol. 2013, article 63, 16 pages, 2013. View at Publisher · View at Google Scholar - Z.-Z. Sun and Z.-B. Zhang, “A linearized compact difference scheme for a class of nonlinear delay partial differential equations,”
*Applied Mathematical Modelling*, vol. 37, no. 3, pp. 742–752, 2013. View at Publisher · View at Google Scholar · View at MathSciNet - L. F. Cordero and R. Escalante, “Segmented tau approximation for test neutral functional differential equations,”
*Applied Mathematics and Computation*, vol. 187, no. 2, pp. 725–740, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - M. Muslim, “Approximation of solutions to history-valued neutral functional differential equations,”
*Computers & Mathematics with Applications*, vol. 51, no. 3-4, pp. 537–550, 2006. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - O. R. Işik, Z. Güney, and M. Sezer, “Bernstein series solutions of pantograph equations using polynomial interpolation,”
*Journal of Difference Equations and Applications*, vol. 18, no. 3, pp. 357–374, 2012. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - E. K. Ifantis, “An existence theory for functional-differential equations and functional-differential systems,”
*Journal of Differential Equations*, vol. 29, no. 1, pp. 86–104, 1978. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - A. Iserles and Y. Liu, “On neutral functional-differential equations with proportional delays,”
*Journal of Mathematical Analysis and Applications*, vol. 207, no. 1, pp. 73–95, 1997. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - C. He, X. Lv, and J. Niu, “A new method based on the RKHSM for solving systems of nonlinear IDDEs with proportional delays,”
*Abstract and Applied Analysis*, vol. 2013, Article ID 541935, 13 pages, 2013. View at Publisher · View at Google Scholar · View at MathSciNet - E. Ishiwata, “On the attainable order of collocation methods for the neutral functional-differential equations with proportional delays,”
*Computing*, vol. 64, no. 3, pp. 207–222, 2000. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - A. H. Bhrawy and A. S. Alofi, “A Jacobi-Gauss collocation method for solving nonlinear Lane-Emden type equations,”
*Communications in Nonlinear Science and Numerical Simulation*, vol. 17, no. 1, pp. 62–70, 2012. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - 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, 13 pages, 2012. View at Publisher · View at Google Scholar - M. Maleki, I. Hashim, M. Tavassoli Kajani, and S. Abbasbandy, “An adaptive pseudospectral method for fractional order boundary value problems,”
*Abstract and Applied Analysis*, vol. 2012, Article ID 381708, 19 pages, 2012. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - 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. View at Google Scholar - W. Wang, Y. Zhang, and S. Li, “Stability of continuous Runge-Kutta-type methods for nonlinear neutral delay-differential equations,”
*Applied Mathematical Modelling*, vol. 33, no. 8, pp. 3319–3329, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - J. X. Kuang and Y. H. Cong,
*Stability of Numerical Methods for Delay Differential Equations*, Science Press, Beijing, China, 2005. - X. Lv and Y. Gao, “The RKHSM for solving neutral functional-differential equations with proportional delays,”
*Mathematical Methods in the Applied Sciences*, vol. 36, no. 6, pp. 642–649, 2013. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - X. Chen and L. Wang, “The variational iteration method for solving a neutral functional-differential equation with proportional delays,”
*Computers & Mathematics with Applications*, vol. 59, no. 8, pp. 2696–2702, 2010. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - W. Wang, T. Qin, and S. Li, “Stability of one-leg $\theta $-methods for nonlinear neutral differential equations with proportional delay,”
*Applied Mathematics and Computation*, vol. 213, no. 1, pp. 177–183, 2009. View at Publisher · View at Google Scholar · View at MathSciNet - Z. Jackiewicz and B. Zubik-Kowal, “Spectral collocation and waveform relaxation methods for nonlinear delay partial differential equations,”
*Applied Numerical Mathematics*, vol. 56, no. 3-4, pp. 433–443, 2006. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet