High-Order Shifted Gegenbauer Integral Pseudospectral Method for Solving Differential Equations of Lane-Emden Type

Abstract

We present a novel, high-order, efficient, and exponentially convergent shifted Gegenbauer integral pseudospectral method (SGIPSM) to solve numerically Lane-Emden equations provided with some mixed Neumann and Robin boundary conditions. The framework of the proposed method includes: (i) recasting the problem into its integral formulation, (ii) collocating the latter at the shifted flipped-Gegenbauer-Gauss-Radau (SFGGR) points, and (iii) replacing the integrals with accurate and well-conditioned numerical quadratures constructed via SFGGR-based shifted Gegenbauer integration matrices. The integral formulation is eventually discretized into linear/nonlinear system of equations that can be solved easily using standard direct system solvers. The implementation of the proposed method is further illustrated through four efficient computational algorithms. Moreover, we furnish rigorous error and convergence analyses of the SGIPSM. Five numerical test examples are presented to verify the effectiveness, accuracy, exponential convergence, and numerical stability of the proposed method. The numerical simulations are associated with extensive numerical comparisons with other rival methods in the literature to demonstrate further the power of the proposed method. The SGIPSM is broadly applicable and represents a strong addition to common numerical methods for solving linear/nonlinear differential equations when high-order approximations are required using a relatively small number of collocation points.

arXiv comment26 pages, 5 figures

References in corpus (6)

• An approximation algorithm for the solution of the nonlinear Lane-Emden type equations arising in astrophysics using Hermite functions collocation method
• High-Order Numerical Solution of Second-Order One-Dimensional Hyperbolic Telegraph Equation Using a Shifted Gegenbauer Pseudospectral Method
• High-Order, Stable, And Efficient Pseudospectral Method Using Barycentric Gegenbauer Quadratures
• Lebesgue Constants Arising in a Class of Collocation Methods
• Optimal Control of a Parabolic Distributed Parameter System Using a Barycentric Shifted Gegenbauer Pseudospectral Method
• High-Order Adaptive Gegenbauer Integral Spectral Element Method for Solving Nonlinear Optimal Control Problems

Cited by in corpus (1)

• A Direct Integral Pseudospectral Method for Solving a Class of Infinite-Horizon Optimal Control Problems Using Gegenbauer Polynomials and Certain Parametric Maps