Exponentially convergent symbolic algorithm of the functional-discrete method for the fourth order Sturm-Liouville problems with polynomial coefficients

JOURNAL OF COMPUTATIONAL AND APPLIED MATHEMATICS

卷:358

Article

Makarov, Volodymyr L.¹  Romaniuk, Nataliia M.¹ 

[1] Natl Acad Sci Ukraine, Inst Math, Dept Numer Math, 3 Tereshchenkivska Str, UA-01004 Kiev 4, Ukraine

关键词: Fourth order Sturm-Liouville problems;    Eigenvalue problems;    Polynomial coefficients;    Functional-discrete method;    Symbolic algorithm;    Exponential convergence rate;   

DOI  :  10.1016/j.cam.2019.03.024

来源: Elsevier

【 摘 要 】

A new symbolic algorithmic implementation of the functional-discrete (FD-) method is developed and justified for the solution of fourth order Sturm-Liouville problem on a finite interval in the Hilbert space. The eigenvalue problem for the fourth order ordinary differential equation with polynomial coefficients is investigated. The sufficient conditions of an exponential convergence rate of the proposed approach are received. The obtained estimates of the absolute errors of FD-method significantly improve the accuracy of the estimates obtained earlier by I.P. Gavrilyuk, V.L. Makarov and A.M. Popov in 2010. Our algorithm is symbolic and operates with the decomposition coefficients of the eigenfunction corrections in some basis. The number of summands in these decompositions depends on the degree of the potential coefficients and the correction number. Our method uses only the algebraic operations and basic operations on (2 x 1) column vectors and (2 x 2) matrices. The proposed approach does not require solving any boundary value problems and computations of any integrals, unlike the previous variants of FD-method by I.P. Gavrilyuk, V.L. Makarov, A.M. Popov and N.M. Romaniuk in 2010 and 2017. The corrections to eigenpairs are computed exactly as analytical expressions. The numerical examples illustrate the theoretical results. The numerical results obtained with the FD-method are compared with the numerical test results obtained with other existing numerical techniques. (C) 2019 Elsevier B.V. All rights reserved.

【 授权许可】

Free