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.024zbMATH Open1415.65171arXiv1806.09193OpenAlexW2810783808WikidataQ128215712 ScholiaQ128215712MaRDI QIDQ2000638FDOQ2000638

Vladimir L. Makarov, N. M. Romanyuk

Publication date: 28 June 2019

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

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.

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

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

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)

Cites Work

• Algorithm 775
• Title not available (Why is that?)
• Introduction to Numerical Continuation Methods
• Solving frontier problems of physics: the decomposition method
• Title not available (Why is that?)
• Title not available (Why is that?)
• Title not available (Why is that?)
• Accurate solutions of fourth order Sturm-Liouville problems
• Super-exponentially convergent parallel algorithm for a fractional eigenvalue problem of Jacobi-type
• New properties of the FD-method in its applications to the Sturm-Liouville problems
• Application of the FD-method to the solution of the Sturm-Liouville problem with coefficients of special form
• Title not available (Why is that?)
• Superexponentially convergent algorithm for an abstract eigenvalue problem with applications to ordinary differential equations
• How many numerical eigenvalues can we trust?
• A new application of the homotopy analysis method: solving the Sturm-Liouville problems
• Computation of eigenvalues of Sturm-Liouville problems using homotopy perturbation method
• An efficient method for computing eigenelements of Sturm-Liouville fourth-order boundary value problems
• Title not available (Why is that?)
• An efficient technique for finding the eigenvalues of fourth-order Sturm-Liouville problems
• Matrix methods for computing eigenvalues of Sturm-Liouville problems of order four
• Title not available (Why is that?)
• 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

Uses Software

• Digital Library of Mathematical Functions
• SLEUTH

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)