International Journal of Mathematics and Mathematical Sciences

Volume 2013 (2013), Article ID 486013, 9 pages

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

## Fibonacci Collocation Method for Solving High-Order Linear Fredholm Integro-Differential-Difference Equations

Department of Mathematics, Celal Bayar University, Muradiye, 45140 Manisa, Turkey

Received 8 May 2013; Accepted 27 June 2013

Academic Editor: Irena Lasiecka

Copyright © 2013 Ayşe Kurt 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 collocation method based on the Fibonacci polynomials is introduced for the approximate solution of high order-linear Fredholm integro-differential-difference equations with the mixed conditions. The proposed method is analyzed to show the convergence of the method. Some further numerical experiments are carried out to demonstrate the method.

#### 1. Introduction

The integro-differential-difference equations (IDDEs) have been developed very rapidly in recent years. This is an important branch of mathematics which has a lot of interest in many application fields such as engineering, mechanics, physics, astronomy, chemistry, biology, economics, and potential theory, electrostatics [1–14]. Since some IDDEs are hard to solve numerically, they are solved by using the approximated methods. Several numerical methods were used such as the successive approximations, Adomian decomposition, Haar Wavelet, and Tau and Walsh series methods [15–20]. Additionally the Monte Carlo method for linear Fredholm integro-differential-difference equation has been presented by Farnoosh and Ebrahimi [21] and the Direct method based on the Fourier and block-pulse method functions by Asady et al. [22].

Since the beginning of 1994, the Taylor and Chebyshev matrix methods have also been used by Sezer et al. to solve linear differential, Fredholm integral, and Fredholm integro-differential equations [23–35]. Lately, the Fibonacci collocation method has been used to find the approximate solutions of differential, integral, and integro-differential equations [36].

In this study, we consider the approximate solution of the th-order Fredholm integro-differential-difference equations, where , are the integer,, , under the mixed conditions where , ,, and are functions defined on , ; , , , and are suitable constants.

Our aim is to obtain an approximate solution of (1) expressed in the truncated Fibonacci series form: where , , are the unknown Fibonacci coefficients. Here is positive integer such that and , , are the Fibonacci polynomials defined by

#### 2. Fundamental Matrix Relations

Firstly, we can write the Fibonacci polynomials in the matrix form as follows: where If is even, if is odd,

Let us show (1) in the following form: where

##### 2.1. Matrix Relations for the Differential Part

Firstly, we consider the solution and its th derivate in the matrix form: Then, from relations (5) and (11), we can obtain the following matrix form: Similar to (13), from relations (5), (11), and (12), we can find matrix form as To find the matrix in terms of the matrix, we can use the following relation: where Subsequently, by substituting the matrix form (15) into (14), we obtain the matrix relations

##### 2.2. Matrix Relations for the Difference Part

If we put in the relation (11), we have the matrix form It is seen that the relation between the matricesandis where By using the relations (15) and (19), we can get Thus from (14) and (21), we can find By using the expressions (14) and (22), we obtain the matrix form

##### 2.3. Matrix Relations for the Integral Part

Let us find the matrix relation for the Fredholm integral part in (9). The kernel function can be shown by the truncated Fibonacci series, and the truncated Taylor series, where The expressions (24) and (25) can be put in matrix forms as From (11), (27), and (28) we can obtain Thus By substituting the matrix forms (22) and (27) into the integral part in (9), we can have the matrix relation as follows: so that From (5) and (32), we have where If we substitute the matrix relation (5) into (31), we have the matrix form

##### 2.4. Matrix Relations for the Conditions

The corresponding matrix form for the conditions (2) can be shown, by means of (17), as

#### 3. Method of Solution

We can construct the fundamental matrix equation corresponding for (1). For this aim, we substitute the matrix relations (23) and (35) into (9). So we obtain the matrix equation By using in (37) the collocation points defined by the system of the matrix equations is obtained or shortly the fundamental matrix equation becomes where Therefore, the fundamental matrix equation (40) corresponding for (1) can be written as where Equation (42) corresponds to a system of linear algebraic equations with unknown Fibonacci coefficients. Further, we can express the matrix form (36) conditions where To obtain the solution of (1) under the conditions (2), by replacing the row matrices (44) by the last rows of the matrices (42), we have the new augmented matrix If the last m rows of the (30) are replaced, the augmented matrix of the above system is obtained as If rank , then we can write And so, the matrix (thereby the coefficients) is uniquely determined.

#### 4. Accuracy of Solution

We can check the accuracy of the method. The truncated Fibonacci series in (3) have to be approximately satisfying (1). For each , , or If ( is any positive integer) is prescribed, then the truncation limit is increased until the difference at each of the points becomes smaller than the prescribed .

#### 5. Numerical Examples

In this section, several examples are given to illustrate the applicability of the method and all of them are performed on the computer MATLAB. Also, the absolute errors in tables are the values of at selected points.

*Example 1. *Let us consider the linear Fredholm integro-differential-difference equation given by
with the boundary conditions
and the approximate solution by the truncated Fibonacci series
where , , , , , , , and ,
From (38), the collocation points for are computed
and from (40), the fundamental matrix equation of the problem is
where
The augmented matrix for this fundamental matrix equation is calculated as
From (37), the matrix forms for the boundary conditions are
The new augmented matrix based on conditions can be written as
Solving this system, the unknown Fibonacci coefficients are obtained as
Hence, by substituting the Fibonacci coefficients matrix into (11), we have the approximate solution , which is the exact solution

*Example 2 (see [26]). *Consider the following linear Fredholm integro-differential-difference equation with variable coefficients
with the conditions , , and and the exact solution .

So, , , , , , , and
The fundamental matrix equation of the problem from (40) becomes

For the solutions of the problem are found as, respectively,
We compare the solutions found by the present method for , and the absolute errors in Figure 1 and Table 1. It is seen that when we increase integer , the errors decrease.

*Example 3 (see [30]). *Now consider the linear Fredholm integro-differential-difference equation given by
with the initial conditions , , and and the exact solution .

Here
From (40), the fundamental matrix equation of the problem becomes
The solutions obtained for are compared with the exact solution is which are given in Figure 2 and we compare the absolute errors found by present method for in Figure 3. Also, the numerical solution and absolute errors are compared for in Table 2.

#### 6. Conclusion

The Fibonacci collocation method is used to solve the linear Fredholm integro-differential-difference equations numerically. The obtained numerical results show that the accuracy improves when is increased. Tables and figures indicate that as increases, the errors decrease more rapidly. A considerable advantage of the method is that the Fibonacci polynomial coefficients of the solution are found very easily by using computer programs. This method can also be extended to the system of linear integro-differential-difference equation with variable coefficients, but some modifications are required.

#### References

- M. K. Kadalbajoo and K. K. Sharma, “Numerical analysis of boundary-value problems for singularly-perturbed differential-difference equations with small shifts of mixed type,”
*Journal of Optimization Theory and Applications*, vol. 115, no. 1, pp. 145–163, 2002. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - C. E. Elmer and E. S. Van Vleck, “Traveling wave solutions for bistable differential-difference equations with periodic diffusion,”
*SIAM Journal on Applied Mathematics*, vol. 61, no. 5, pp. 1648–1679, 2001. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - D. D. Bainov, M. B. Dimitrova, and A. B. Dishliev, “Oscillation of the bounded solutions of impulsive differential-difference equations of second order,”
*Applied Mathematics and Computation*, vol. 114, no. 1, pp. 61–68, 2000. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - M. K. Kadalbajoo and K. K. Sharma, “Numerical analysis of singularly perturbed delay differential equations with layer behavior,”
*Applied Mathematics and Computation*, vol. 157, no. 1, pp. 11–28, 2004. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - M. T. Rashed, “Numerical solution of functional differential, integral and integro-differential equations,”
*Applied Mathematics and Computation*, vol. 156, no. 2, pp. 485–492, 2004. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - K. Maleknejad and Y. Mahmoudi, “Numerical solution of linear Fredholm integral equation by using hybrid Taylor and block-pulse functions,”
*Applied Mathematics and Computation*, vol. 149, no. 3, pp. 799–806, 2004. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - W. Wang and C. Lin, “A new algorithm for integral of trigonometric functions with mechanization,”
*Applied Mathematics and Computation*, vol. 164, no. 1, pp. 71–82, 2005. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - W. Wang, “An algorithm for solving the high-order nonlinear Volterra-Fredholm integro-differential equation with mechanization,”
*Applied Mathematics and Computation*, vol. 172, no. 1, pp. 1–23, 2006. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - C. E. Elmer and E. S. Van Vleck, “A variant of Newton's method for the computation of traveling waves of bistable differential-difference equations,”
*Journal of Dynamics and Differential Equations*, vol. 14, no. 3, pp. 493–517, 2002. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - P. Rai and K. K. Sharma, “Numerical study of singularly perturbed differential-difference equation arising in the modeling of neuronal variability,”
*Computers & Mathematics with Applications*, vol. 63, no. 1, pp. 118–132, 2012. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - P. Rai and K. K. Sharma, “Parameter uniform numerical method for singularly perturbed differential-difference equations with interior layers,”
*International Journal of Computer Mathematics*, vol. 88, no. 16, pp. 3416–3435, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - V. Kumar and K. K. Sharma, “An optimized B-spline method for solving singularly perturbed differential difference equations with delay as well as advance,”
*Neural, Parallel & Scientific Computations*, vol. 16, no. 3, pp. 371–385, 2008. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - M. K. Kadalbajoo and K. K. Sharma, “An exponentially fitted finite difference scheme for solving boundary-value problems for singularly-perturbed differential-difference equations: small shifts of mixed type with layer behavior,”
*Journal of Computational Analysis and Applications*, vol. 8, no. 2, pp. 151–171, 2006. View at Google Scholar · View at MathSciNet - M. K. Kadalbajoo and K. K. Sharma, “Numerical analysis of boundary-value problems for singularly perturbed differential-difference equations: small shifts of mixed type with rapid oscillations,”
*Communications in Numerical Methods in Engineering*, vol. 20, no. 3, pp. 167–182, 2004. View at Publisher · View at Google Scholar · View at MathSciNet - M. H. Reihani and Z. Abadi, “Rationalized Haar functions method for solving Fredholm and Volterra integral equations,”
*Journal of Computational and Applied Mathematics*, vol. 200, no. 1, pp. 12–20, 2007. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - M. H. AliAbadi and E. L. Ortiz, “Numerical treatment of moving and free boundary value problems with the tau method,”
*Computers & Mathematics with Applications*, vol. 35, no. 8, pp. 53–61, 1998. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - J. Pour-Mahmoud, M. Y. Rahimi-Ardabili, and S. Shahmorad, “Numerical solution of the system of Fredholm integro-differential equations by the Tau method,”
*Applied Mathematics and Computation*, vol. 168, no. 1, pp. 465–478, 2005. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - S. M. Hosseini and S. Shahmorad, “Numerical solution of a class of integro-differential equations by the tau method with an error estimation,”
*Applied Mathematics and Computation*, vol. 136, no. 2-3, pp. 559–570, 2003. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - K. Maleknejad and F. Mirzaee, “Numerical solution of integro-differential equations by using rationalized Haar functions method,”
*Kybernetes*, vol. 35, no. 10, pp. 1735–1744, 2006. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - S. M. Hosseini and S. Shahmorad, “Tau numerical solution of Fredholm integro-differential equations with arbitrary polynomial bases,”
*Applied Mathematical Modelling*, vol. 27, no. 2, pp. 145–154, 2003. View at Publisher · View at Google Scholar · View at Scopus - R. Farnoosh and M. Ebrahimi, “Monte Carlo method for solving Fredholm integral equations of the second kind,”
*Applied Mathematics and Computation*, vol. 195, no. 1, pp. 309–315, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - B. Asady, M. Tavassoli Kajani, A. Hadi Vencheh, and A. Heydari, “Solving second kind integral equations with hybrid Fourier and block-pulse functions,”
*Applied Mathematics and Computation*, vol. 160, no. 2, pp. 517–522, 2005. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - M. Sezer, “Taylor polynomial solutions of Volterra integral equations,”
*International Journal of Mathematical Education in Science and Technology*, vol. 25, no. 5, pp. 625–633, 1994. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - S. Yalçınbaş and M. Sezer, “The approximate solution of high-order linear Volterra-Fredholm integro-differential equations in terms of Taylor polynomials,”
*Applied Mathematics and Computation*, vol. 112, no. 2-3, pp. 291–308, 2000. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - M. Sezer and M. Gülsu, “Polynomial solution of the most general linear Fredholm integrodifferential-difference equations by means of Taylor matrix method,”
*Complex Variables*, vol. 50, no. 5, pp. 367–382, 2005. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - M. Gülsu and M. Sezer, “Approximations to the solution of linear Fredholm integrodifferential-difference equation of high order,”
*Journal of the Franklin Institute*, vol. 343, no. 7, pp. 720–737, 2006. View at Publisher · View at Google Scholar · View at MathSciNet - S. Yalçinbaş, “Taylor polynomial solutions of nonlinear Volterra-Fredholm integral equations,”
*Applied Mathematics and Computation*, vol. 127, no. 2-3, pp. 195–206, 2002. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - A. Akyüz-Daşcıoğlu and M. Sezer, “Chebyshev polynomial solutions of systems of higher-order linear Fredholm-Volterra integro-differential equations,”
*Journal of the Franklin Institute*, vol. 342, no. 6, pp. 688–701, 2005. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - S. Yalçinbaş, M. Sezer, and H. H. Sorkun, “Legendre polynomial solutions of high-order linear Fredholm integro-differential equations,”
*Applied Mathematics and Computation*, vol. 210, no. 2, pp. 334–349, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - M. Gülsu and Y. Öztürk, “A new collocation method for solution of mixed linear integro-differential-difference equations,”
*Applied Mathematics and Computation*, vol. 216, no. 7, pp. 2183–2198, 2010. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - S. Yüzbaşı, N. Şahin, and M. Sezer, “Bessel collocation method for solving high-order linear Fredholm integro-differential-difference equations,”
*Journal of Advanced Research in Differential Equations*, vol. 3, no. 2, pp. 1–23, 2011. View at Google Scholar · View at MathSciNet - K. Erdem, S. Yalçinbas, and M. Sezer, “A new collocation method for solution of linear Fredholm integro differential-difference equations using Bernoulli polynomials,” in
*Proceedings of the International Conference on Applied Analysis and Algebra (ICAAA)*, Yildiz Technical University, Istanbul, Turkey, 1992. - A. Akyüz-Daşcıoğlu and M. Sezer, “A Taylor polynomial approach for solving the most general linear Fredholm integro-differential-difference equations,”
*Mathematical Methods in the Applied Sciences*, vol. 35, no. 7, pp. 839–844, 2012. View at Publisher · View at Google Scholar · View at MathSciNet - S. Yalçınbaş, N. Özsoy, and M. Sezer, “Approximate solution of higher order linear differential equations by means of a new rational Chebyshev collocation method,”
*Mathematical & Computational Applications*, vol. 15, no. 1, pp. 45–56, 2010. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - S. Yalçınbaş, M. Aynigül, and T. Akkaya, “Legendre series solutions of Fredholm integral equations,”
*Mathematical & Computational Applications*, vol. 15, no. 3, pp. 371–381, 2010. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - A. Kurt,
*Fibonacci polynomial solutions of linear differential, integral and integro-differential equations [M.S. thesis]*, Graduate School of Natural and Applied Sciences, Mugla University, Mugla, Turkey, 2012.