Research Article  Open Access
Application of Rational Second Kind Chebyshev Functions for System of Integrodifferential Equations on SemiInfinite Intervals
Abstract
Rational Chebyshev bases and Galerkin method are used to obtain the approximate solution of a system of highorder integrodifferential equations on the interval [0,∞). This method is based on replacement of the unknown functions by their truncated series of rational Chebyshev expansion. Test examples are considered to show the high accuracy, simplicity, and efficiency of this method.
1. Introduction
In recent years, there has been a growing interest in the system of integrodifferential equations (IDE), which arise frequently in many applied areas which include engineering, mechanics, physics, chemistry, astronomy, biology, economics, potential theory, electrostatics, and so forth [1–8]. The systems of integrodifferential equations are generally difficult to solve analytically, thus finding efficient computational algorithms for obtaining numerical solution is required.
There are various techniques for solving systems of IDE, for example, operational Tau method [9, 10], Adomian decomposition method [11], Galerkin method [12], rationalized Haar functions method [13], He’s homotopy perturbation method (HPM) [14, 15], and Ghebyshev polynomial [16].
A number of problems arising in science and engineering are set in semiinfinite domains. One can apply different spectral methods that are used to solve problems in semiinfinite domains. The first approach is using Laguerre polynomials [17–20]. The second approach is replacing semiinfinite domain with interval by choosing , sufficiently large. This method is named domain truncation [21]. The third approach is reformulating original problem in semiinfinite domain to singular problem in bounded domain by variable transformation and then using the Jacobi polynomials to approximate the resulting singular problem [22]. The fourth approach of spectral method is based on rational orthogonal functions.
Boyd [23] defined a new spectral basis, named rational Chebyshev functions on the semiinfinite interval, by mapping to the Chebyshev polynomials. Guo et al. [24] introduced a new set of rational Legendre functions which are mutually orthogonal in . They applied a spectral scheme using the rational Legendre functions for solving the Kortewegde Vries equation on the half line. Boyd et al. [25] applied pseudospectral methods on a semiinfinite interval and compared rational Chebyshev, Laguerre, and mapped Fourier sine.
The authors of [26–29] applied spectral method to solve nonlinear ordinary differential equations on semiinfinite intervals. Their approach was based on a rational Tau method. They obtained the operational matrices of derivative and product of rational Chebyshev and Legendre functions and then they applied these matrices together with the Tau method to reduce the solution of these problems to the solution of system of algebraic equations.
Zarebnia and Ali Abadi [30] used SincCollocation method for solving system of nonlinear secondorder integrodifferential equations of the Fredholm type. Rational second (third) kind Chebyshev (RSC) functions, for the first time, were proposed by Tavassoli Kajani and Ghasemi Tabatabaei [31] to find the numerical solution of LaneEmden equation.
This paper outlines the application of rational second kind Chebyshev functions and Galerkin method to the following system of linear highorder integrodifferential equations on the interval . Two problems of such equations are solved to make clear the application of the proposed method. One has where .
2. Properties of RSC Functions
In this section, we present some properties of RSC functions.
2.1. RSC Functions
The second kind Chebyshev polynomials , are orthogonal in the interval with respect to the weight function and we find that satisfies the recurrence relation [32] The RSC functions are defined by [31, 33] thus RSC functions satisfy
2.2. Function Approximation
Let denotes a nonnegative, integrable, realvalued function over the interval . We define where is the norm induced by the scalar product Thus denote a system which is mutually orthogonal under (2.6), that is, where is the Kronecker delta function. This system is complete in ; as a result, any function can be expanded as follows: with The ’s are the expansion coefficients associated with the family . If the infinite series in (2.8) is truncated, then it can be written as where and .
We can also approximate the function in as follows: where is an matrix that Moreover, from recurrence relation in (2.3) we have
2.3. Product Integration of the RSC Functions
We also use the matrix as follows: To illustrate the calculation we choose and , then we obtain
2.4. Operational Matrix of Derivative
The derivative of the vector defined in (2.10) can be approximated by where is named the operational matrix of derivative. Differentiating (2.3) we get By using (2.17) the matrix can be calculated. The matrix is a lower Hessenberg matrix and can be expressed as , where is a tridiagonal matrix which is obtained from and the elements of matrix are obtained from
2.5. The Product Operational Matrix
The following property of the product of two rational Chebyshev vectors will also be used: where is called product operational matrix for the vector . Using (2.20) and the orthogonal property, the elements of the matrix can be calculated from where is given by
3. Solving System of Integrodifferential Equations over SemiInfinite Interval
Consider the following system of integrodifferential equations: Using (2.10) and (2.11) to approximate , , , and when and , we have According to the operational matrix of derivative we can approximate as With substituting these approximations in (3.1) we have Then using (2.20) we obtain which can be simplified using (2.14) By solving this linear system of algebraic equations we can find vectors , , and then approximate the solutions as
4. Numerical Examples
Example 4.1. Consider the following system of linear integrodifferential equations:
The exact solution of this example is and .
We solved Example 4.1 using the present method with , and we obtained and , which imply
that are the exact solutions.
Example 4.2. Next, consider the following system of integrodifferential equations with the exact solution and : We solved this example by using the method described in Section 3 for and . Results are shown in Table 1 and Figures 1 and 2. The errors for large values of are shown in Table 2. It is seen that the proposed method provides accurate results even for large values of .

5. Conclusion
The fundamental goal of this paper has been to construct an approximation to the solution of the integrodifferential equations system in a semiinfinite interval. In the above discussion, the Galerkin method with RSC functions, which have the property of orthogonality, is employed to achieve this goal. Advantages of this method is that we do not reform the problem to a finite domain, and with a small value of accurate results are obtained. There is a good agreement between obtained results and exact values that demonstrates the validity of the present method for this type of problems and gives the method a wider applicability.
References
 J. Abdul Jerri, Introduction to Integral Equations with Applications, John Wiley & Sons, New York, NY, USA, 1999. View at: MathSciNet
 J. P. Boyd, “Evaluating of Dawson's integral by solving its differential equation using orthogonal rational Chebyshev functions,” Applied Mathematics and Computation, vol. 204, no. 2, pp. 914–919, 2008. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 P. Linz, Analytical and Numerical Methods for Volterra Equations, vol. 7, Society for Industrial and Applied Mathematics, Philadelphia, Pa, USA, 1985. View at: Publisher Site  MathSciNet
 K. Maleknejad and Y. Mahmoudi, “Numerical solution of linear Fredholm integral equation by using hybrid Taylor and blockpulse functions,” Applied Mathematics and Computation, vol. 149, no. 3, pp. 799–806, 2004. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 Y. Ren, B. Zhang, and H. Qiao, “A simple Taylorseries expansion method for a class of second kind integral equations,” Journal of Computational and Applied Mathematics, vol. 110, no. 1, pp. 15–24, 1999. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 M. T. Rashed, “Numerical solution of functional differential, integral and integrodifferential equations,” Applied Mathematics and Computation, vol. 156, no. 2, pp. 485–492, 2004. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 W. Wang, “An algorithm for solving the highorder nonlinear VolterraFredholm integrodifferential equation with mechanization,” Applied Mathematics and Computation, vol. 172, no. 1, pp. 1–23, 2006. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 S. Yalçınbaş and M. Sezer, “The approximate solution of highorder linear VolterraFredholm integrodifferential equations in terms of Taylor polynomials,” Applied Mathematics and Computation, vol. 112, no. 23, pp. 291–308, 2000. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 S. Abbasbandy and A. Taati, “Numerical solution of the system of nonlinear Volterra integrodifferential equations with nonlinear differential part by the operational Tau method and error estimation,” Journal of Computational and Applied Mathematics, vol. 231, no. 1, pp. 106–113, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 G. Ebadi, M. Y. Rahimi, and S. Shahmorad, “Numerical solution of the system of nonlinear Fredholm integrodifferential equations by the operational tau method with an error estimation,” Scientia Iranica, vol. 14, no. 6, pp. 546–554, 2007. View at: Google Scholar  Zentralblatt MATH  MathSciNet
 A. Khani, M. Mohseni Moghadam, and S. Shahmorad, “Numerical solution of special class of systems of nonlinear Volterra integrodifferential equations by a simple high accuracy method,” Bulletin of the Iranian Mathematical Society, vol. 34, no. 2, pp. 141–152, 2008. View at: Google Scholar  MathSciNet
 K. Maleknejad and M. Tavassoli Kajani, “Solving linear integrodifferential equation system by Galerkin methods with hydrid functions,” Applied Mathematics and Computation, vol. 159, no. 3, pp. 603–612, 2004. View at: Publisher Site  Google Scholar  MathSciNet
 K. Maleknejad, F. Mirzaee, and S. Abbasbandy, “Solving linear integrodifferential equations system by using rationalized Haar functions method,” Applied Mathematics and Computation, vol. 155, no. 2, pp. 317–328, 2004. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 J. Biazar, H. Ghazvini, and M. Eslami, “He's homotopy perturbation method for systems of integrodifferential equations,” Chaos, Solitons & Fractals, vol. 39, no. 3, pp. 1253–1258, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 E. Yusufoğlu, “An efficient algorithm for solving integrodifferential equations system,” Applied Mathematics and Computation, vol. 192, no. 1, pp. 51–55, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 A. AkyüzDaşcıoğlu and M. Sezer, “Chebyshev polynomial solutions of systems of higherorder linear FredholmVolterra integrodifferential equations,” Journal of the Franklin Institute, vol. 342, no. 6, pp. 688–701, 2005. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 B.Y. Guo and J. Shen, “LaguerreGalerkin method for nonlinear partial differential equations on a semiinfinite interval,” Numerische Mathematik, vol. 86, no. 4, pp. 635–654, 2000. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 Y. Maday, B. PernaudThomas, and H. Vandeven, “ShockFitting techniques for solving hyperbolic problems with spectral methods,” Recherche Aerospatiale, vol. 6, pp. 1–9, 1985. View at: Google Scholar  Zentralblatt MATH
 H. I. Siyyam, “Laguerre tau methods for solving higherorder ordinary differential equations,” Journal of Computational Analysis and Applications, vol. 3, no. 2, pp. 173–182, 2001. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 J. Shen, “Stable and efficient spectral methods in unbounded domains using Laguerre functions,” SIAM Journal on Numerical Analysis, vol. 38, no. 4, pp. 1113–1133, 2000. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 J. P. Boyd, Chebyshev and Fourier Spectral Methods, Springer, Berlin, Germany, 2nd edition, 2000. View at: MathSciNet
 B.Y. Guo, “Jacobi spectral approximations to differential equations on the half line,” Journal of Computational Mathematics, vol. 18, no. 1, pp. 95–112, 2000. View at: Google Scholar  Zentralblatt MATH  MathSciNet
 J. P. Boyd, “Orthogonal rational functions on a semiinfinite interval,” Journal of Computational Physics, vol. 70, no. 1, pp. 63–88, 1987. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 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. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 J. P. Boyd, C. Rangan, and P. H. Bucksbaum, “Pseudospectral methods on a semiinfinite interval with application to the hydrogen atom: a comparison of the mapped Fouriersine method with Laguerre series and rational Chebyshev expansions,” Journal of Computational Physics, vol. 188, no. 1, pp. 56–74, 2003. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 K. Parand and M. Razzaghi, “Rational Chebyshev tau method for solving higherorder ordinary differential equations,” International Journal of Computer Mathematics, vol. 81, no. 1, pp. 73–80, 2004. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 H. I. Siyyam, “Laguerre Tau methods for solving higherorder ordinary differential equations,” Journal of Computational Analysis and Applications, vol. 3, no. 2, pp. 173–182, 2001. View at: Publisher Site  Google Scholar
 K. Parand and M. Razzaghi, “Rational Chebyshev tau method for solving Volterra's population model,” Applied Mathematics and Computation, vol. 149, no. 3, pp. 893–900, 2004. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 K. Parand and M. Razzaghi, “Rational legendre approximation for solving some physical problems on semiinfinite intervals,” Physica Scripta, vol. 69, no. 5, pp. 353–357, 2004. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 M. Zarebnia and M. G. Ali Abadi, “Numerical solution of system of nonlinear secondorder integrodifferential equations,” Computers & Mathematics with Applications, vol. 60, no. 3, pp. 591–601, 2010. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 M. Tavassoli Kajani and F. Ghasemi Tabatabaei, “Rational Chebyshev approximations for solving LaneEmde equation of index m,” in Proceeding of the International Conference on Computational and Applied Mathematics, pp. 840–844, Bangkok, Thailand, March 2011. View at: Google Scholar
 M. Abramowitz and I. A. Stegun, Handbook of Mathematical Functions, Dover, New York, NY, USA, 10th edition, 1972.
 M. Dadkhah Tirani, F. Ghasemi Tabatabaei, and M. Tavassoli Kajani, “Rational second (third) kind Chebyshev approximations for solving Volterras population model,” in Proceeding of the International Conference on Computational and Applied Mathematics, pp. 835–839, Bangkok, Thailand, March 2011. View at: Google Scholar
Copyright
Copyright © 2012 M. Tavassoli Kajani 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.