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

DOI10.1016/J.CAM.2019.03.024MaRDI QIDQ2000638zbMATH OpenOpenAlexWikidataFDO

Authors Vladimir L. Makarov, N. M. Romanyuk

Publication date 28 June 2019

Published in Journal of Computational and Applied Mathematics (Search for Journal in Brave)

Full work available at URL https://arxiv.org/abs/1806.09193

zbMATH Keywords

eigenvalue problemspolynomial coefficientsfunctional-discrete methodsymbolic algorithmexponential convergence ratefourth order Sturm-Liouville problems

Mathematics Subject Classification ID

Boundary eigenvalue problems for ordinary differential equations (34B09) Boundary value problems for linear higher-order PDEs (35G15) Stability and convergence of numerical methods for ordinary differential equations (65L20) Error bounds for numerical methods for ordinary differential equations (65L70) Sturm-Liouville theory (34B24) Numerical approximation of eigenvalues and of other parts of the spectrum of ordinary differential operators (34L16) Numerical solution of eigenvalue problems involving ordinary differential equations (65L15)

Abstract: 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 $(2imes1)$ column vectors and $(2imes2)$ 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, and there are no rounding errors. 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.

Recommendations

• Symbolic algorithm of the functional-discrete method for a Sturm-Liouville problem with a polynomial potential
• Functional-discrete Method (FD-method) for Matrix Sturm-Liouville Problems
• New implementation of the FD-method for Sturm-Liouville problems with Dirichlet-Neumann boundary conditions
• FD-method for Sturm-Liouville problems. Exponential rate of convergence
• The FD-method for solving Sturm-Liouville problems with special singular differential operator

Cites work

• scientific article; zbMATH DE number 3812514 (Why is no real title available?)
• scientific article; zbMATH DE number 3573250 (Why is no real title available?)
• scientific article; zbMATH DE number 490142 (Why is no real title available?)
• scientific article; zbMATH DE number 3253261 (Why is no real title available?)
• scientific article; zbMATH DE number 3353287 (Why is no real title available?)
• scientific article; zbMATH DE number 7693782 (Why is no real title available?)
• A new application of the homotopy analysis method: solving the Sturm-Liouville problems
• Accurate solutions of fourth order Sturm-Liouville problems
• Algorithm 775
• An efficient method for computing eigenelements of Sturm-Liouville fourth-order boundary value problems
• An efficient technique for finding the eigenvalues of fourth-order Sturm-Liouville problems
• Application of the FD-method to the solution of the Sturm-Liouville problem with coefficients of special form
• Computation of eigenvalues of Sturm-Liouville problems using homotopy perturbation method
• How many numerical eigenvalues can we trust?
• Introduction to Numerical Continuation Methods
• Matrix methods for computing eigenvalues of Sturm-Liouville problems of order four
• NIST handbook of mathematical functions
• New properties of the FD-method in its applications to the Sturm-Liouville problems
• Solving frontier problems of physics: the decomposition method
• Super-exponentially convergent parallel algorithm for a fractional eigenvalue problem of Jacobi-type
• Superexponentially convergent algorithm for an abstract eigenvalue problem with applications to ordinary differential equations
• Symbolic algorithm of the functional-discrete method for a Sturm-Liouville problem with a polynomial potential

Cited in

(2)

• Generalization of resonant equations for the Laguerre- and Legendre-type polynomials to equations of the fourth order
• Symbolic algorithm of the functional-discrete method for a Sturm-Liouville problem with a polynomial potential

Describes a project that uses

Uses Software

• SLEUTH
• DLMF

This page was built for publication: Exponentially convergent symbolic algorithm of the functional-discrete method for the fourth order Sturm-Liouville problems with polynomial coefficients

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2000638)