Research Article  Open Access
Numerical and Analytical Study for FourthOrder IntegroDifferential Equations Using a Pseudospectral Method
Abstract
A numerical method for solving fourthorder integrodifferential equations is presented. This method is based on replacement of the unknown function by a truncated series of wellknown shifted Chebyshev expansion of functions. An approximate formula of the integer derivative is introduced. The introduced method converts the proposed equation by means of collocation points to system of algebraic equations with shifted Chebyshev coefficients. Thus, by solving this system of equations, the shifted Chebyshev coefficients are obtained. Special attention is given to study the convergence analysis and derive an upper bound of the error of the presented approximate formula. Numerical results are performed in order to illustrate the usefulness and show the efficiency and the accuracy of the present work.
1. Introduction
The integrodifferential equation (IDE) is an equation that involves both integrals and derivatives of an unknown function. Mathematical modeling of reallife problems usually results in functional equations, like ordinary or partial differential equations, and integral and integrodifferential equations, stochastic equations. Many mathematical formulations of physical phenomena contain integrodifferential equations; these equations arise in many fields like physics, astronomy, potential theory, fluid dynamics, biological models, and chemical kinetics. Integrodifferential equations; are usually difficult to solve analytically; so, it is required to obtain an efficient approximate solution [1–5]. Recently, several numerical methods to solve IDEs have been given such as variational iteration method [6, 7], homotopy perturbation method [8, 9], spline functions expansion [10, 11], and collocation method [12–15].
Chebyshev polynomials are wellknown family of orthogonal polynomials on the interval that have many applications [4, 6, 8, 13]. They are widely used because of their good properties in the approximation of functions. However, with our best knowledge, very little work was done to adapt these polynomials to the solution of integrodifferential equations. Orthogonal polynomials have a great variety and wealth of properties. Some of these properties take a very concise form in the case of the Chebyshev polynomials, making Chebyshev polynomials of leading importance among orthogonal polynomials. The Chebyshev polynomials belong to an exclusive band of orthogonal polynomials, known as Jacobi polynomials, which correspond to weight functions of the form and which are solutions of SturmLiouville equations [16].
In this work, we derive an approximate formula of the integral derivative and derive an upper bound of the error of this formula, and then we use this formula to solve a class of twopoint boundary value problems (BVPs) for the fourthorder integrodifferential equations as under the boundary and initial conditions where are known functions and , and are suitable constants. Several numerical methods to solve the fourthorder integrodifferential equations have been given such as Chebyshev cardinal functions [17], variational iteration method [7], and others.
2. Some Basic Properties and Derivation of an Approximate Formula of the Derivative for Chebyshev Polynomials Expansion
The Chebyshev polynomial of the first kind is a polynomial in of degree , defined by the relation The Chebyshev polynomials of degree of the first kind have precisely zeros and local extrema in the interval . The zeros of are denoted by The Chebyshev polynomials can be determined with the aid of the following recurrence formula [18]: The analytic form of the Chebyshev polynomials of degree is given by where denotes the integer part of . The orthogonality condition is In order to use these polynomials on the interval , we define the so called shifted Chebyshev polynomials by introducing the change of variable . The shifted Chebyshev polynomials are denoted by and defined as .
The function , which belongs to the space of square integrable in , may be expressed in terms of shifted Chebyshev polynomials as where the coefficients are given by
In practice, only the first terms of shifted Chebyshev polynomials are considered. Then, we have that
Lemma 1. The analytic form of the shifted Chebyshev polynomials of degree is given by
Proof. Since we have , then by substituting in (6), we can obtain that Now, we put in (12) we obtain the desired result (11).
The main approximate formula of the derivative of , and is given in the following theorem.
Theorem 2. Let be approximated by shifted Chebyshev polynomials as (10), and also suppose that is integer; then, where is given by
Proof. Since the differential operator is linear, we can obtain that Since is a constant, and Then, we have that and for , and by using (16), we get that A combination of (17), (18), and (14) leads to the desired result and completes the proof of the theorem.
3. Error Analysis
In this section, special attention is given to study the convergence analysis and evaluate the upper bound of the error of the proposed formula.
Theorem 3 (Chebyshev truncation theorem; see [18]). The error in approximating by the sum of its first terms is bounded by the sum of the absolute values of all the neglected coefficients. If then for all , all , and all .
Theorem 4. The derivative of order for the shifted Chebyshev polynomials can be expressed in terms of the shifted Chebyshev polynomials themselves in the following form: where
Proof. We use the properties of the shifted Chebyshev polynomials [18] and expand in (18) in the following form:
where can be obtained using (9), and; then,
At , we find that = = ; also, at any and using the formula (10), we can find that
employing (18) and (23) gives
where After some lengthy manipulation, can be put in the following form:
and this completes the proof of the theorem.
Theorem 5. The error in approximating by is bounded by
Proof. A combination of (8), (10), and (21) leads to but ; so, we can obtain that and subtracting the truncated series from the infinite series, bounding each term in the difference, and summing the bounds complete the proof of the theorem.
4. Procedure Solution for the FourthOrder IntegroDifferential Equation
In this section, we will present the proposed method to solve numerically the fourthorder integrodifferential equation of the form in (1). The unknown function may be expanded by finite series of shifted Chebyshev polynomials as in the following approximation: and approximated formula of its derivatives can be defined in Theorem 2. From (1), (32), and Theorem 2, we have that We now collocate (33) at points , as For suitable collocation points, we use roots of shifted Chebyshev polynomial . The integral terms in (34) can be found using composite trapezoidal integration technique as where ,, for an arbitrary integer . So, by using (34) and (35), we obtain Also, by substituting (32) in the boundary conditions (2), we can obtain equations as follows: Equation (36), together with equations of the boundary conditions (37), give of system of algebraic equations which can be solved, for the unknowns , using conjugate gradient method or Newton iteration method.
5. Numerical Results
In this section, to verify the validity and the accuracy and support our theoretical discussion of the proposed method, we give some computations results of numerical examples.
Example 6. Consider the nonlinear fourthorder integrodifferential equation as in (1) and (2) with ; then, the integrodifferential equation will be
subject to the boundary conditions
The exact solution of this problem is [7].
We apply the suggested method with and approximate the solution as follows:
From (38), (40), and Theorem 2, we have that
We now collocate (41) at points, , as
For suitable collocation points we use roots of shifted Chebyshev polynomial . The integral terms in (42) can be found using composite trapezoidal integration technique as
where ,, for an arbitrary integer . So, by using (43) and (42), we obtain
Also, by substituting (40) in the boundary conditions (39), we can obtain four equations as follows:
where and .
Equation (44), together with four equations of the boundary conditions (45), represent, a nonlinear system of six algebraic equations in the coefficients ; by solving it using the Newton iteration method, we obtain
The behavior of the approximate solution using the proposed method with , the approximate solution using variational iteration method (VIM), and the exact solution are presented in Figure 1. Table 1 shows the behavior of the absolute error between exact solution and approximate solution using the presented method at and . From Figure 1 and Table 1, it is clear that the proposed method can be considered as an efficient method to solve the nonlinear integrodifferential equations. Table 1 indicates that as increases the errors decrease more rapidly; hence, for better results, using number is recommended. Also, we can conclude that the obtained approximated solution is in excellent agreement with the exact solution.

Example 7. Consider the linear fourthorder integrodifferential equation as in (1) and (2) with ; then, the integrodifferential equation will be
subject to the boundary conditions
The exact solution of this problem is [17].
We apply the suggested method with and approximate the solution as follows:
By the same procedure in the previous example, we have
where , and the nodes , . We can write the initialboundary conditions in the form
By using (50) and (51), we obtain a linear system of seven algebraic equations in the coefficients ; by solving it using the conjugate gradient method, we obtain
The behavior of the approximate solution using the proposed method with , the approximate solution using variational iteration method (VIM) and the exact solution are presented in Figure 2. Table 2 shows the behaviour of the absolute error between exact solution and approximate solution using the presented method at and . From this figure, it is clear that the proposed method can be considered as an efficient method to solve the linear integrodifferential equations. Also, we can conclude that the obtained approximate solution is in excellent agreement with the exact solution.

6. Conclusion and Discussion
Integrodifferential equations are usually difficult to solve analytically; so, it is required to obtain the approximate solution. In this paper, we proposed the pseudospectral method using shifted Chebyshev method for solving the integrodifferential equations. The Chebyshev method is useful for acquiring both the general solution and particular solution as demonstrated in examples. Special attention is given to study the convergence analysis and derive an upper bound of the error of the derived approximate formula. From our obtained results, we can conclude that the proposed method gives solutions in excellent agreement with the exact solution and better than the other methods. An interesting feature of this method is that when an integral system has linearly independent polynomial solution of degree or less than , the method can be used for finding the analytical solution. All computations are done using MATLAB 8.
References
 R. P. Agarwal, “Boundary value problems for higher order integrodifferential equations,” Nonlinear Analysis: Theory, Methods & Applications, vol. 7, no. 3, pp. 259–270, 1983. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 E. Babolian, F. Fattahzadeh, and E. G. Raboky, “A Chebyshev approximation for solving nonlinear integral equations of Hammerstein type,” Applied Mathematics and Computation, vol. 189, no. 1, pp. 641–646, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 A. H. Borzabadi, A. V. Kamyad, and H. H. Mehne, “A different approach for solving the nonlinear Fredholm integral equations of the second kind,” Applied Mathematics and Computation, vol. 173, no. 2, pp. 724–735, 2006. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 M. M. Khader, “On the numerical solutions for the fractional diffusion equation,” Communications in Nonlinear Science and Numerical Simulation, vol. 16, no. 6, pp. 2535–2542, 2011. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 N. H. Sweilam, M. M. Khader, and A. M. Nagy, “Numerical solution of twosided spacefractional wave equation using finite difference method,” Journal of Computational and Applied Mathematics, vol. 235, no. 8, pp. 2832–2841, 2011. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 M. M. Khader, “Introducing an efficient modification of the variational iteration method by using Chebyshev polynomials,” Applications and Applied Mathematics, vol. 7, no. 1, pp. 283–299, 2012. View at: Google Scholar  Zentralblatt MATH  MathSciNet
 N. H. Sweilam, “Fourth order integrodifferential equations using variational iteration method,” Computers & Mathematics with Applications, vol. 54, no. 78, pp. 1086–1091, 2007. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 M. M. Khader, “Introducing an efficient modification of the homotopy perturbation method by using Chebyshev polynomials,” Arab Journal of Mathematical Sciences, vol. 18, no. 1, pp. 61–71, 2012. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 N. H. Sweilam, M. M. Khader, and R. F. AlBar, “Homotopy perturbation method for linear and nonlinear system of fractional integrodifferential equations,” International Journal of Computational Mathematics and Numerical Simulation, vol. 1, no. 1, pp. 73–87, 2008. View at: Google Scholar
 M. M. Khader and S. T. Mohamed, “Numerical treatment for first order neutral delay differential equations using spline functions,” Engineering Mathematics Letters, vol. 1, no. 1, pp. 32–43, 2012. View at: Google Scholar
 S. T. Mohamed and M. M. Khader, “Numerical solutions to the second order Fredholm integrodifferential equations using the spline functions expansion,” Global Journal of Pure and Applied Mathematics, vol. 34, pp. 21–29, 2011. View at: Google Scholar
 M. M. Khader and A. S. Hendy, “The approximate and exact solutions of the fractionalorder delay differential equations using Legendre pseudospectral method,” International Journal of Pure and Applied Mathematics, vol. 74, no. 3, pp. 287–297, 2012. View at: Google Scholar
 M. M. Khader, N. H. Sweilam, and A. M. S. Mahdy, “An efficient numerical method for solving the fractional difusion equation,” Journal of Applied Mathematics and Bioinformatics, vol. 1, no. 2, pp. 1–12, 2011. View at: Google Scholar
 N. H. Sweilam, M. M. Khader, and W. Y. Kota, “On the numerical solution of Hammerstein integral equations using Legendre approximation,” International Journal of Applied Mathematical Research, vol. 1, no. 1, pp. 65–76, 2012. View at: Google Scholar
 S. Yousefi and M. Razzaghi, “Legendre wavelets method for the nonlinear VolterraFredholm integral equations,” Mathematics and Computers in Simulation, vol. 70, no. 1, pp. 1–8, 2005. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 J. C. Mason and D. C. Handscomb, Chebyshev Polynomials, Chapman & Hall/CRC, Washington, DC, USA, 2003. View at: MathSciNet
 M. Lakestani and M. Dehghan, “Numerical solution of fourthorder integrodifferential equations using Chebyshev cardinal functions,” International Journal of Computer Mathematics, vol. 87, no. 6, pp. 1389–1394, 2010. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 M. A. Snyder, Chebyshev Methods in Numerical Approximation, PrenticeHall Inc., Englewood Cliffs, NJ, USA, 1966. View at: MathSciNet
Copyright
Copyright © 2013 N. H. Sweilam 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.