Abstract and Applied Analysis

Volume 2014 (2014), Article ID 816473, 9 pages

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

Research Article

A Pseudospectral Algorithm for Solving Multipantograph Delay Systems on a Semi-Infinite Interval Using Legendre Rational Functions

E. H. Doha,¹ D. Baleanu,2,3,4 A. H. Bhrawy,5,6 and R. M. Hafez⁷

¹Department of Mathematics, Faculty of Science, Cairo University, Giza 12613, Egypt

²Department of Chemical and Materials Engineering, Faculty of Engineering, King Abdulaziz University, Jeddah 21589, Saudi Arabia

³Department of Mathematics and Computer Sciences, Cankaya University, Ankara, Turkey

⁴Institute of Space Sciences, 76900 Magurele-Buchaarest, Romania

⁵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

⁷Department of Basic Science, Institute of Information Technology, Modern Academy, Cairo 11931, Egypt

Received 26 December 2013; Accepted 25 April 2014; Published 19 May 2014

Academic Editor: Robert A. Van Gorder

Copyright © 2014 E. H. Doha 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.

A new Legendre rational pseudospectral scheme is proposed and developed for solving numerically systems of linear and nonlinear multipantograph equations on a semi-infinite interval. A Legendre rational collocation method based on Legendre rational-Gauss quadrature points is utilized to reduce the solution of such systems to systems of linear and nonlinear algebraic equations. In addition, accurate approximations are achieved by selecting few Legendre rational-Gauss collocation points. The numerical results obtained by this method have been compared with various exact solutions in order to demonstrate the accuracy and efficiency of the proposed method. Indeed, for relatively limited nodes used, the absolute error in our numerical solutions is sufficiently small.

1. Introduction

Over the last three decades, the scientists have paid much attention to spectral methods due to their high accuracy (see, for instance, [1–6] and the references therein). On the other hand, spectral methods, in the context of numerical schemes for differential equations, generically belong to the family of weighted residual methods (WRMs) (cf. Finlayson [7]). WRMs represent a particular group of approximation techniques, in which the residuals (or errors) are minimized in a certain way and thereby leading to specific methods including Galerkin, Petrov-Galerkin, collocation, and tau formulations. WRMs are traditionally regarded as the foundation and cornerstone of the finite element, spectral, finite volume, boundary element, and some other methods. Many problems in science and engineering arise in unbounded domains (see, e.g., [8–12]). In general, the use of Jacobi rational functions has the advantage of obtaining the solutions in terms of the Jacobi rational parameters (see, e.g., [13–16]). Moreover, the authors of [17, 18] proposed an efficient collocation schemes based on the operational matrices of rational Legendre and Chebyshev functions for solving problems in the half line.

Systems of multipantograph equations model many real-life phenomena in physics and biology. Numerous numerical and analytical schemes have been presented and developed for solving such multipantograph equations, like variational iteration [19], Bernstein collocation [20], spline [21], homotopy perturbation [22], and Taylor collocation [23, 24] methods. This problem has been the focus of many studies; see, for instance, [21, 25–29].

The main aim of this paper is to develop the spectral-Gauss collocation integration process based on Legendre rational functions for solving systems of linear and nonlinear multipantograph equations on the half line. The implementation of this method reduces the problem to systems of linear and nonlinear algebraic equations. The numerical simulation of multipantograph systems on a semi-infinite domain is investigated. Finally, the accuracy and applicability of the proposed method are demonstrated by test problems. Numerical results are given in which the exponential convergence behaviour is exhibited.

This paper is organized as follows. We present few revelent properties of Legendre rational functions and function approximation in the coming section. In Sections 3 and 4, the pseudospectral algorithms are implemented for solving systems of linear and nonlinear multipantograph equations. Two test examples are introduced in Section 5. Finally, some concluding remarks are given in the last section.

2. Legendre Rational Interpolation

In this section, we present Legendre rational functions and Legendre rational approximation that will be used to construct the Legendre rational-Gauss collocation (LR-GC) method.

2.1. Legendre Rational Functions

The well-known Legendre polynomials are defined on the interval with respect to the weight function . In order to use these polynomials on the interval , we recall the Legendre rational functions by introducing the change of variable . Let the Legendre rational functions be denoted by . Then can be obtained with the aid of the following recurrence formula:

According to the properties of the standard Legendre polynomials, we have

2.2. Function Approximation

Let denote a nonnegative, integrable, real-valued function over the interval . We define

where

is the norm induced by the inner product of the space :

Thus, denotes a system which is mutually orthogonal under (6); that is,

where is the Kronecker delta function. This system is complete in . For any function the following expansion holds:

with

2.3. Legendre Rational Interpolation Approximation

We denote by , , the nodes of the standard Legendre-Gauss interpolation on the interval , and , are Christoffel numbers. The nodes of the Legendre rational-Gauss interpolation on the interval are the zeros of , which is denoted by , . It is clear that , and , . We now set . Therefore, making use of the property of the Legendre-Gauss quadrature, we have for any ,

The interpolating function of a smooth function on a semi-infinite interval is denoted by . It is an element of and is defined as

is the orthogonal projection of upon with respect to the inner product (6) and the norm (5). Thus by the orthogonality of Legendre rational functions we have [11]

To obtain the order of convergence of Legendre rational approximation, at first we define the space

where the norm is induced by

and is the Sturm-Liouville operator as follows:

We have the following theorem for the convergence.

Theorem 1. For any and ,

A complete proof of the theorem and discussion on convergence are given in [11].

3. Linear Multipantograph System

In this section, we propose the Legendre rational-Gauss collocation method to solve the following system of linear multipantograph equations:

subject to

Let us first introduce some basic notation that will be used in the sequel. We define the discrete inner product and norm as follows:

Obviously,

The Legendre rational-Gauss collocation method for solving (17) and (18) is to seek , , such that

We derive the algorithm for solving (17)-(18). To do this, let

Now, substitution of (22) into (17) enables us to write

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

Moreover, the initial condition (18)-with the aid of (2)-yields

If we collocate (24) at the Legendre rational roots of , then we get

Thus (26) with relation (25) generate of a set of algebraic equations which can be solved for the unknown coefficients , by using any standard solver technique.

4. Nonlinear Multipantograph System

In this section, we consider the nonlinear multipantograph system of the form

with initial conditions

The Legendre rational-Gauss collocation method for solving (27) and (28) is to seek , such that

is satisfied exactly at the collocation points , . In other words, we have to collocate (29) at the Legendre rational roots , which immediately yields

with (28) written in the form

This constitute a system of nonlinear algebraic equations in the unknown expansion coefficients , which can be solved by using any standard iteration technique, like Newton's iteration method.

5. Numerical Results

In this section, we present two numerical examples in order to show the accuracy of Legendre rational collocation method for solving multipantograph delay system. The following tables and figures contain the values of the exact solutions , the approximate solutions , and the absolute error functions , at the selected points of the given interval.

Example 1. Consider the following linear multipantograph delay system:

with the initial conditions

where

The exact solution of the system is and .

Tables 1 and 2 list the results obtained by the Legendre rational-Gauss collocation method in terms of maximum absolute errors with different values of . The logarithmic graphs of absolute coefficients for Legendre rational functions are shown in Figures 1 and 2. This confirms that the proposed method has reasonable convergence rate.

Table 1: Comparison of the absolute errors for of (32).

LR-GC method

0.0
0.000 × 10⁺⁰⁰
0.000 × 10⁺⁰⁰
2.775 × 10⁻¹⁷

10.0
1.105 × 10⁻⁴
5.638 × 10⁻⁸
3.779 × 10⁻¹⁰

20.0
2.175 × 10⁻⁵
4.347 × 10⁻⁷
2.609 × 10⁻⁹

30.0
3.087 × 10⁻³
6.613 × 10⁻⁷
8.111 × 10⁻⁹

40.0
1.080 × 10⁻²
6.894 × 10⁻⁷
1.604 × 10⁻⁸

50.0
2.036 × 10⁻²
1.644 × 10⁻⁶
4.074 × 10⁻⁸

60.0
3.008 × 10⁻²
6.281 × 10⁻⁶
9.717 × 10⁻⁸

70.0
3.922 × 10⁻²
1.594 × 10⁻⁵
1.190 × 10⁻⁷

80.0
4.756 × 10⁻²
3.053 × 10⁻⁵
8.903 × 10⁻⁸

90.0
5.506 × 10⁻²
4.921 × 10⁻⁵
4.396 × 10⁻⁸

100.0
6.178 × 10⁻²
7.100 × 10⁻⁵
3.624 × 10⁻⁸

Table 2: Comparison of the absolute errors for of (32).

LR-GC method

0.0
8.273 × 10⁻⁷
9.056 × 10⁻¹²
4.393 × 10⁻¹¹

10.0
1.395 × 10⁻⁵
6.763 × 10⁻⁹
4.703 × 10⁻¹²

20.0
9.412 × 10⁻⁶
8.911 × 10⁻⁹
5.661 × 10⁻¹¹

30.0
6.376 × 10⁻⁵
1.177 × 10⁻⁸
1.242 × 10⁻¹⁰

40.0
1.499 × 10⁻⁴
7.887 × 10⁻⁹
2.575 × 10⁻¹⁰

50.0
2.121 × 10⁻⁴
1.621 × 10⁻⁸
3.371 × 10⁻¹⁰

60.0
2.529 × 10⁻⁴
4.544 × 10⁻⁸
6.993 × 10⁻¹⁰

70.0
2.791 × 10⁻⁴
1.027 × 10⁻⁷
9.146 × 10⁻¹⁰

80.0
2.957 × 10⁻⁴
1.820 × 10⁻⁷
7.991 × 10⁻¹⁰

90.0
3.061 × 10⁻⁴
2.729 × 10⁻⁷
5.234 × 10⁻¹⁰

100.0
3.125 × 10⁻⁴
3.658 × 10⁻⁷
3.478 × 10⁻¹⁰

Figure 1: Logarithmic graph of absolute coefficients of Legendre rational functions at of Example 1.

Figure 2: Logarithmic graph of absolute coefficients of Legendre rational functions at of Example 1.

Example 2. Consider the following nonlinear multipantograph delay system:

with the initial conditions

where

The exact solution of the system is and .

In Tables 3 and 4, we list the absolute errors obtained by the Legedre rational-Gauss collocation method, with different values of . Figures 3 and 4 are plotted to compare the analytic solution with the approximate solution at .

Table 3: Comparison of the absolute errors for of (35).

LR-GC method

0.0
5.551 × 10⁻¹⁷
0.000 × 10⁺⁰⁰
0.000 × 10⁺⁰⁰

1.0
1.554 × 10⁻³
4.799 × 10⁻⁶
8.262 × 10⁻⁸

2.0
2.085 × 10⁻³
2.349 × 10⁻⁶
4.957 × 10⁻⁸

3.0
6.583 × 10⁻³
1.392 × 10⁻⁶
1.069 × 10⁻⁷

4.0
1.408 × 10⁻³
9.788 × 10⁻⁹
1.247 × 10⁻⁷

5.0
5.932 × 10⁻³
2.183 × 10⁻⁵
8.508 × 10⁻⁷

6.0
6.690 × 10⁻³
6.643 × 10⁻⁵
1.169 × 10⁻⁶

7.0
3.789 × 10⁻³
6.919 × 10⁻⁶
1.485 × 10⁻⁶

8.0
1.092 × 10⁻³
6.729 × 10⁻⁵
6.698 × 10⁻⁷

9.0
1.546 × 10⁻³
2.765 × 10⁻⁵
1.382 × 10⁻⁶

10.0
2.137 × 10⁻⁴
2.345 × 10⁻⁵
1.450 × 10⁻⁶

Table 4: Comparison of the absolute errors for of (35).

LR-GC method

0.0
0.000 × 10⁺⁰⁰
1.110 × 10⁻¹⁶
1.665 × 10⁻¹⁶

1.0
6.260 × 10⁻⁴
2.819 × 10⁻⁶
5.391 × 10⁻⁸

2.0
6.997 × 10⁻⁵
2.048 × 10⁻⁵
4.171 × 10⁻⁸

3.0
5.800 × 10⁻⁵
1.042 × 10⁻⁶
9.088 × 10⁻⁸

4.0
6.901 × 10⁻⁴
7.814 × 10⁻⁶
1.541 × 10⁻⁷

5.0
2.427 × 10⁻³
9.524 × 10⁻⁶
2.014 × 10⁻⁷

6.0
3.012 × 10⁻³
1.126 × 10⁻⁶
1.116 × 10⁻⁷

7.0
1.483 × 10⁻³
1.283 × 10⁻⁵
2.553 × 10⁻⁷

8.0
8.598 × 10⁻³
4.171 × 10⁻⁶
5.499 × 10⁻⁷

9.0
1.529 × 10⁻²
6.713 × 10⁻⁶
2.949 × 10⁻⁷

10.0
1.995 × 10⁻²
3.503 × 10⁻⁵
3.311 × 10⁻⁷

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

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

6. Conclusion

In this paper, a collocation Legendre rational method has been proposed to obtain the approximate solutions of systems of multipantograph delay equations. The derivation of this method is essentially based on Legendre rational functions and Gauss quadrature formula. The main advantage of the developed method is that high accurate solutions were achieved using few numbers of the Legendre rational functions. Additionally, if is increased, it can be seen that approximate solutions obtained by the method are close to the exact solutions.

Conflict of Interests

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

Acknowledgments

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

References

C. Canuto, M. Y. Hussaini, A. Quarteroni, and T. A. Zang, Spectral Methods: Fundamentals in Single Domains, Springer, New York, NY, USA, 2006.

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.

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.

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.

A. B. Koc, M. Cakmak, A. Kurnaz, and K. Uslu, “A new Fibonacci type collocation procedure for boundary value problems,” Advances in Difference Equations, vol. 2013, article 262, 2013.

Y. Simsek, “Unification of the Bernstein-type polynomials and their applications,” Boundary Value Problems, vol. 2013, article 56, 2013.

B. A. Finlayson, The Method of Weighted Residuals and Variational Principles, Academic Press, New York, NY, USA, 1972.

D. Funaro, “Computational aspects of pseudospectral Laguerre approximations,” Applied Numerical Mathematics, vol. 6, no. 6, pp. 447–457, 1990.

B.-y. Guo, “Jacobi approximations in certain Hilbert spaces and their applications to singular differential equations,” Journal of Mathematical Analysis and Applications, vol. 243, no. 2, pp. 373–408, 2000.

J. P. Boyd, Chebyshev and Fourier Spectral Methods, Dover, New York, NY, USA, 2nd edition, 2000.

B.-Y. Guo, J. Shen, and Z.-Q. Wang, “A rational approximation and its applications to differential equations on the half line,” Journal of Scientific Computing, vol. 15, no. 2, pp. 117–147, 2000.

B.-Y. Guo, J. Shen, and Z.-Q. Wang, “Chebyshev rational spectral and pseudospectral methods on a semi-infinite interval,” International Journal for Numerical Methods in Engineering, vol. 53, no. 1, pp. 65–84, 2002.

E. H. Doha, A. H. Bhrawy, R. M. Hafez, and R. A. Van Gorder, “A Jacobi rational pseudospectral method for Lane-Emden initial value problems an astrophysics on a semi-infinite interval,” Computational and Applied Mathematics, 2013.

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, R. M. Hafez, and R. A. Van Gorder, “Jacobi rational-Gauss collocation method for Lane-Emden equations of astrophysical significance,” Nonlinear Analysis: Modelling and Control. In press

Y.-G. Yi and B.-Y. Guo, “Generalized Jacobi rational spectral method on the half line,” Advances in Computational Mathematics, vol. 37, no. 1, pp. 1–37, 2012.

K. Parand, M. Shahini, and M. Dehghan, “Rational Legendre pseudospectral approach for solving nonlinear differential equations of Lane-Emden type,” Journal of Computational Physics, vol. 228, no. 23, pp. 8830–8840, 2009.

K. Parand and M. Razzaghi, “Rational Chebyshev tau method for solving higher-order ordinary differential equations,” International Journal of Computer Mathematics, vol. 81, no. 1, pp. 73–80, 2004.

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.

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.

A. El-Safty and S. M. Abo-Hasha, “On the application of spline functions to initial value problems with retarded argument,” International Journal of Computer Mathematics, vol. 32, pp. 173–179, 1990.

J. Biazar and B. Ghanbari, “The homotopy perturbation method for solving neutral functional-differential equations with proportional delays,” Journal of King Saud University: Science, vol. 24, no. 1, pp. 33–37, 2012.

M. Sezer, S. Yalçinbaş, and N. Şahin, “Approximate solution of multi-pantograph equation with variable coefficients,” Journal of Computational and Applied Mathematics, vol. 214, no. 2, pp. 406–416, 2008.

M. Gülsu and M. Sezer, “A Taylor collocation method for solving high-order linear pantograph equations with linear functional argument,” Numerical Methods for Partial Differential Equations, vol. 27, no. 6, pp. 1628–1638, 2011.

Y. Kuang, Delay Differential Equations with Applications in Population Dynamics, vol. 191, Academic Press, New York, NY, USA, 1993.

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.

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.

Z.-H. Yu, “Variational iteration method for solving the multi-pantograph delay equation,” Physics Letters. A, vol. 372, no. 43, pp. 6475–6479, 2008.

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.