Advances in Mathematical Physics

Advances in Mathematical Physics / 2014 / Article

Research Article | Open Access

Volume 2014 |Article ID 694580 | 5 pages |

A Matrix Method Based on the Fibonacci Polynomials to the Generalized Pantograph Equations with Functional Arguments

Academic Editor: Pavel Kurasov
Received12 Apr 2014
Accepted31 Jul 2014
Published13 Aug 2014


A pseudospectral method based on the Fibonacci operational matrix is proposed to solve generalized pantograph equations with linear functional arguments. By using this method, approximate solutions of the problems are easily obtained in form of the truncated Fibonacci series. Some illustrative examples are given to verify the efficiency and effectiveness of the proposed method. Then, the numerical results are compared with other methods.

1. Introduction

Many phenomena in applied branches that fail to be modeled by the ordinary differential equations can be described by the delay differential equations. Many researchers have studied different applications of those equations in variety of applied sciences such as biology, physics, economy, and electrodynamics (see [14]). Pantograph equations with proportional delays play an important role in this context. The existence and uniqueness of the analytic solutions of the multipantograph equation are investigated in [5]. A numerical approach to multipantograph equations with variable coefficients is also studied in [6]. An extension of the multipantograph equation is known to be the generalized pantograph equation with functional arguments defined as under the mixed conditions where proportional delay-, constant delay-, , , and are real and/or complex coefficients, and the coefficients of th order unknown function- and known are the analytical functions defined in the interval .

In recent years, many researchers have developed different numerical approaches to the generalized pantograph equations as variational iteration method [7], differential transform approach [8], Taylor method [9], collocation method based on Bernoulli matrix [10], and Bessel collocation method [11]. In this study, we investigate a collocation method based on the Fibonacci polynomial operational matrix for the numerical solution of the generalized pantograph equation (1). Even the Fibonacci numbers have been known for a long time; the Fibonacci polynomials are very recently defined to be an important agent in the world of polynomials [12, 13]. Compared to the methods of the orthogonal polynomials, the Fibonacci approach has proved to give more precise and reliable results in the solution of differential equations [14].

This study is organized as follows. In the second part, a short review of the Fibonacci polynomials is presented. A Fibonacci operational matrix for the solution of the pantograph equation is developed in Section 3. Some numerical examples are given in Section 4 to illustrate efficiency and effectiveness of the method.

2. Operational Matrices of the Fibonacci Polynomials

The Fibonacci polynomials are determined by following general formula [12, 13]: with and . Now, we will mention some matrix relations in terms of Fibonacci polynomials.

2.1. Fibonacci Series Expansions

To obtain an expansion form of the analytic solution of the pantograph equation, we use the Fibonacci collocation method as follows.

Suppose that (1) has a continuous function solution that can be expressed in the Fibonacci polynomials Then, a truncated expansion of -Fibonacci polynomials can be written in the vector form where the Fibonacci row vector and the unknown Fibonacci coefficients column vector are given, respectively, by

2.2. Matrix Relations of the Derivatives

The th order derivative of (5) can be written as where , , and is the coefficient vector of the polynomial approximation of th order derivative. Then, there exists a relation between the Fibonacci coefficients as where is operational matrix for the derivative defined by [14] Making use of (7) and (9) yields

3. Solution Procedure for the Pantograph Differential Equations

Let us recall the th order linear pantograph differential equation, The first step in the solution procedure is to define the collocation points in the domain, so that Then, collocating problem (12) at the points in (13) yields The system (14) can, alternatively, be rewritten in the matrix form where Therefore, the th order derivative of the unknown function at the collocation points can be written in the matrix form as or, equivalently, To express the functional terms of (1) as in the form (5), let we put instead of in the relation (18) and then obtain or where are Fibonacci operational matrices corresponding to the coefficients . Therefore, replacing (18) and (20) in (15) gives the fundamental matrix equation for the problem (12) as which corresponds to a system of algebraic equations for the unknown Fibonacci coefficients ,  . In other words, when we denote the expression in the sum by , for and , we get Thus, the augmented matrix of (22) becomes

On the other hand, in view of (11), the conditions (2) can be taken into account by forming the following matrix equation: where Therefore, the augmented matrix of the specified conditions is Consequently, (23) together with (26) can be written in the new augmented matrix form This form can also be achieved by replacing some rows of the matrix (23) by the rows of (26) or adding those rows to the matrix (23) provided that . Finally, the vector (thereby vector of the coefficients ) is determined by applying some numerical methods designed especially to solve the system of linear equations. On the other hand, when the singular case appears, the least square methods are inevitably available to reach the best possible approximation. Therefore, the approximated solution can be obtained. This would be the Fibonacci series expansion of the solution to the problem (12) with the specified conditions.

3.1. Accuracy of the Results

We can, now, proceed with a short accuracy analysis of the problem in a similar way to [18]. As the truncated Fibonacci series expansion is an approximate solution of (1) with (2), it must satisfy the following equality for , : or When ( is any integer) is prescribed, the truncation limit is increased until the difference at each of the collocation points becomes smaller than the desired value .

4. Numerical Results

In this part, three illustrative examples are given in order to clarify the findings of the previous section. The errors of the proposed method are compared with those of the errors that occurred in the solutions by some other methods in Tables 13 for two sample examples. It is noted here that the number of collocation points in the examples is indicated by the capital letter .

Present methodExponential approach [15]Taylor polynomial approach [16]


Present methodTaylor polynomial approach [16]Taylor method [9]


Present methodTaylor matrix method [6] Boubaker matrix method [17]


Example 1 (see [10, 11]). Consider the following linear pantograph type problem equation: with the initial conditions

The exact solution of this problem is known to be . When the solution procedure in Section 3 is applied to the problem, the solution of the linear algebraic system gives the numerical approximation of the solution to the problem. It is noteworthy that the method reaches the exact solution even for .

Example 2 (see [9, 15, 16]). Now, consider the following equation with variable coefficient given by and the condition

The exact solution is also known to be . A comparison of the absolute errors of the proposed approach Taylor method [16] and the exponential approach [15] is given in Table 1 for . Another comparison of the present method with the methods of Taylor polynomials [9, 16] is also given for and in Table 2. These results show that the Fibonacci approach has better accuracy, at least one decimal place, than the other methods.

Example 3 (see [6, 17]). Finally, let us consider the pantograph equation with variable coefficients and the condition which has the exact solution . Computed results are compared with the results of the Taylor [6] and Boubaker [17] matrix methods in Table 3.

Conflict of Interests

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


This study was supported by Research Projects Center (BAP) of Selcuk University. The authors would like to thank Selcuk University and TUBITAK for their support. Also, the authors would like to thank the editor and referees for their valuable comments and remarks which led to a great improvement of the article. They have denoted here that a minor part of this study was presented orally at the “2nd International Eurasian Conference on Mathematical Sciences and Applications (IECMSA-2013),” Sarajevo, August, 2013.


  1. W. G. Ajello, H. I. Freedman, and J. Wu, “A model of stage structured population growth with density dependent time delay,” SIAM Journal on Applied Mathematics, vol. 52, pp. 855–869, 1992. View at: Google Scholar
  2. Y. Kuang, Delay Differential Equations with Applications in Population Dynamics, Boston, Mass, USA, Academic Press, 1993.
  3. M. Dehghan and F. Shakeri, “The use of the decomposition procedure of Adomian for solving a delay differential equation arising in electrodynamics,” Physica Scripta, vol. 78, no. 6, Article ID 065004, 2008. View at: Publisher Site | Google Scholar | Zentralblatt MATH
  4. J. R. Ockendon and A. B. Tayler, “The dynamics of a current collection system for an electric locomotive,” Proceedings of the Royal Society of London A, vol. 322, pp. 447–468, 1971. View at: Publisher Site | Google Scholar
  5. 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 Site | Google Scholar | Zentralblatt MATH | MathSciNet
  6. M. Sezer, S. Yalcinbas, and N. Sahin, “Approximate solution of multi pantograph equation with variable coefficients,” Journal of Computational and Applied Mathematics, vol. 214, no. 2, pp. 406–416, 2008. View at: Publisher Site | Google Scholar
  7. A. Saadatmandi and M. Dehghan, “Variational iteration method for solving a generalized pantograph equation,” Computers & Mathematics with Applications, vol. 58, no. 11-12, pp. 2190–2196, 2009. View at: Publisher Site | Google Scholar | MathSciNet
  8. Y. Keskin, A. Kurnaz, M. E. Kiris, and G. Oturanc, “Approximate solutions of generalized pantograph equations by the differential transform method,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 8, no. 2, pp. 159–164, 2007. View at: Google Scholar
  9. M. Sezer and A. Akyüz-Daşcıoğlu, “A Taylor method for numerical solution of generalized pantograph equations with linear functional argument,” Journal of Computational and Applied Mathematics, vol. 200, no. 1, pp. 217–225, 2007. View at: Publisher Site | Google Scholar | MathSciNet
  10. 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. Simulation and Computation for Engineering and Environmental Systems, vol. 37, no. 6, pp. 4283–4294, 2013. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
  11. S. Yuzbası, N. Sahin, 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, 2011. View at: Publisher Site | Google Scholar | MathSciNet
  12. S. Falcón and Á. Plaza, “The k-Fibonacci sequence and the Pascal 2-triangle,” Chaos, Solitons & Fractals, vol. 33, no. 1, pp. 38–49, 2007. View at: Publisher Site | Google Scholar | MathSciNet
  13. S. Falcon and A. Plaza, “On K-Fibonacci sequences and polynomials and their derivatives,” Chaos, Solitons & Fractals, vol. 39, no. 3, pp. 1005–1019, 2009. View at: Publisher Site | Google Scholar | MathSciNet
  14. A. B. Koç, M. Çakmak, M. Kurnaz, and K. Uslu, “A new Fibonacci type collocation procedure for boundary value problems,” Advances in Difference Equations, vol. 2013, article 262, 2013. View at: Publisher Site | Google Scholar | MathSciNet
  15. S. Yuzbası and M. Sezer, “An exponential approximation for solutions of generalized pantograph-delay differential equations,” Applied Mathematical Modelling, vol. 37, no. 22, pp. 9160–9173, 2013. View at: Publisher Site | Google Scholar | MathSciNet
  16. M. Sezer, S. Yalçınbas, and M. Gülsu, “A Taylor polynomial approach for solving generalized pantograph equations with nonhomogenous term,” International Journal of Computer Mathematics, vol. 85, no. 7, pp. 1055–1063, 2008. View at: Publisher Site | Google Scholar | MathSciNet
  17. T. Akkaya, S. Yalçinbaş, and M. Sezer, “Numeric solutions for the pantograph type delay differential equation using First Boubaker polynomials,” Applied Mathematics and Computation, vol. 219, no. 17, pp. 9484–9492, 2013. View at: Publisher Site | Google Scholar | MathSciNet
  18. N. Baykus and M. Sezer, “Solution of high-order linear Fredholm integro-differential equations with piecewise intervals,” Numerical Methods for Partial Differential Equations, vol. 27, no. 5, pp. 1327–1339, 2011. View at: Publisher Site | Google Scholar | MathSciNet

Copyright © 2014 Ayşe Betül Koç 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.

1134 Views | 653 Downloads | 2 Citations
 PDF  Download Citation  Citation
 Download other formatsMore
 Order printed copiesOrder

We are committed to sharing findings related to COVID-19 as quickly and safely as possible. Any author submitting a COVID-19 paper should notify us at to ensure their research is fast-tracked and made available on a preprint server as soon as possible. We will be providing unlimited waivers of publication charges for accepted articles related to COVID-19.