Galerkin-like method for solving linear functional differential equations under initial conditions

In this paper, we present a weighted residual Galerkin method to solve linear functional differential equations. We consider the problem with variable coefficients under initial conditions. Assuming the exact solution of the problem has a Taylor series expansion convergent in the relevant domain, we seek a solution of the given problem in the form of a polynomial having degree $N$ of our choice. Substituting this polynomial with unknown coefficients in the given equation yields an expression linear in these coefficients. We then proceed as in the weighted residual method and take inner product of this expression with monomials up to degree $N$, resulting in $N+1$ linear algebraic equations. Appropriately incorporating the initial conditions and solving the resulting linear system, we obtain the approximate solution to the given problem. Additionally, we present a way of estimating the absolute error of the obtained approximation, which is then used to improve the original approximation through a method called residual correction. We also show that the upper bound for the error of the proposed method depends on the Taylor truncation error of the exact solution. The proposed scheme and the residual correction technique are illustrated in several example problems.

___

  • [1] Ockendon JR, Tayler AB. The dynamics of a current collection system for an electric locomotive. Proceedings of the Royal Society of London A 1971; 322: 447-468. doi: 10.1098/rspa.1971.0078
  • [2] Tayler AB. Mathematical Models in Applied Mechanics. New York, NY, USA: Oxford University Press, 1985.
  • [3] Brunner H. Collocation Methods for Volterra Integral and Related Functional Equations. Cambridge, UK: Cambridge University Press, 2004.
  • [4] Derfel GA, Vogl F. On the asymptotics of solutions of a class of linear functional-differential equations. European Journal of Applied Mathematics 1996; 7 (5): 511-518. doi: 10.1017/S0956792500002527
  • [5] Feldstein A. On neutral functional-differential equations with variable time delays. Mathematical Proceedings of the Cambridge Philosophical Society 1998; 124 (2): 371-398. doi: 10.1017/S0305004198002497
  • [6] Dehghan M, Shakourifar M, Hamidi A. The solution of linear and non-linear systems of Volterra functional equations using Adomian-Padé technique. Chaos Solitons & Fractals 2009; 39 (5): 2509-2521. doi: 10.1016/j.chaos.2007.07.028
  • [7] Feng X. An analytic study on the multi-pantograph delay equations with variable coefficients. Bulletin Mathématiques de la Société des Sciences Mathématiques de Roumanie 2013; 56 (2): 205-215.
  • [8] Yüzbaşı Ş, Karaçayır M. A Galerkin-like approach to solve multi-pantograph type delay differential equations. FILOMAT 2018; 32 (2): 409-422. doi: 10.2298/FIL1802409Y
  • [9] Ishiwata E, Muroya Y. On collocation methods for delay differential and Volterra integral equations with proportional delay. Frontiers of Mathematics in China 2009; 4 (1): 89-111. doi: 10.1007/s11464-009-0004-x
  • [10] Işık OR, Güney Z, Sezer M. Bernstein series solutions of pantograph equations using polynomial interpolation. Journal of Difference Equations and Applications 2012; 18 (3): 357-374. doi: 10.1080/10236198.2010.496456
  • [11] Yüzbaşı Ş, Şahin N, Sezer M. A Bessel collocation method for numerical solution of generalized pantograph equations. Numerical Methods for Partial Differential Equations 2011; 28 (4): 110 5-1123. doi: 10.1002/num.20660
  • [12] Sedaghat S, Ordokhani Y, Dehghan M. Numerical solution of the delay differential equations of pantograph type via Chebyshev polynomials. Communications in Nonlinear Science and Numerical Simulation 2012; 17 (12): 4815-4830. doi: 10.1016/j.cnsns.2012.05.009
  • [13] Yüzbaşı Ş, Sezer M. An exponential approximation for solutions of generalized pantograph-delay differential equations. Applied Mathematical Modelling 2013; 37 (22): 9160-9173. doi: 10.1016/j.apm.2013.04.028
  • [14] Tohidi E, Bhrawy AH, Erfani K. A collocation method based on Bernoulli operational matrix for numerical solution of generalized pantograph equation. Applied Mathematical Modelling 2013; 37 (6): 4283-4294. doi: 10.1016/j.apm.2012.09.032
  • [15] Sezer M, Yalçınbaş S, Gülsu M. A Taylor polynomial approach for solving generalized pantograph equations with nonhomogenous term. International Journal of Computer Mathematics 2011; 85 (7): 1055-1063. doi: 10.1080/00207160701466784
  • [16] Baykuş Savaşaneril N, Sezer M. Hybrid Taylor-Lucas collocation method for numerical solution of high-order pantograph type delay differential equations with variable delays. Applied Mathematics & Information Sciences 2017; 11 (6): 1795-1801. doi: 10.18576/amis/110627
  • [17] Saadatmandi A, Dehghan M. Variational iteration method for solving a generalized pantograph equation. Computers & Mathematics with Applications 2009; 58 (11-12): 2190-2196. doi: 10.1016/j.camwa.2009.03.017
  • [18] Yüzbaşı Ş, İsmailov N. A Taylor operation method for solutions of generalized pantograph type delay differential equations. Turkish Journal of Mathematics 2018; 42: 395-406. doi: 10.3906/mat-1506-71
  • [19] Javadi S, Babolian E, Taheri Z. Solving generalized pantograph equations by shifted orthonormal Bernstein polynomials. Journal of Computational and Applied Mathematics 2016; 303: 1-14. doi: 10.1016/j.cam.2016.02.025
  • [20] Gümgüm S, Baykuş Savaşaneril N, Kürkçü ÖK, Sezer M. Lucas polynomial solution of nonlinear differential equations with variable delays. Hacettepe Journal of Mathematics & Statistics 2019; doi:10.15672/hujms.460975
  • [21] Türkyılmazoğlu M. An effective approach for numerical solutions of high-order Fredholm integro-differential equations. Applied Mathematics and Computation 2014; 227: 384-398. doi: 10.1016/j.amc.2013.10.079
  • [22] Yüzbaşı Ş, Karaçayır M. A Galerkin-like approach to solve high-order integro-differential equations with weakly singular kernel. Kuwait Journal of Science 2016; 43 (2): 61-75.
  • [23] Sezer M, Akyüz Daşçıoğlu A. A Taylor method for numerical solution of generalized pantograph equations with linear functional argument. Journal of Computational and Applied Mathematics 2007; 200 (1): 217-225. doi: 10.1016/j.cam.2005.12.015
  • [24] Yalçınbaş S, Aynıgül M, Sezer M. A collocation method using Hermite polynomials for approximate solution of pantograph equations. Journal of the Franklin Institute 2011; 348 (6): 1128-1139. doi: 10.1016/j.jfranklin.2011.05.003
  • [25] Öztürk Y, Gülsu M. Approximate solution of linear generalized pantograph equations with variable coefficients on Chebyshev-Gauss grid. Journal of Advanced Research in Scientific Computing 2012; 4 (1): 36-51.