We describe a strategy for rigorous arbitrary-precision evaluation of Legendre polynomials on the unit interval and its application in the generation of Gauss-Legendre quadrature rules. Our focus is on making the evaluation practical for a wide range of realistic parameters, corresponding to the requirements of numerical integration to an accuracy of about 100 to 100 000 bits. Our algorithm com...

We extend a collocation method for solving a nonlinear ordinary differential equation ODE via Jacobi polynomials. To date, researchers usually use Chebyshev or Legendre collocation method for solving problems in chemistry, physics, and so forth, see the works of Doha and Bhrawy 2006, Guo 2000, and Guo et al. 2002 . Choosing the optimal polynomial for solving every ODEs problem depends on many f...

The main aim of this article is to generalize the Legendre operational matrix to the fractional derivatives and implemented it to solve the nonlinear multi-order fractional differential equations. In this approach, a truncated Legendre series together with the Legendre operational matrix of fractional derivatives are used. The main characteristic behind the approach using this technique is that...

A new algorithm for solving the general nonlinear third-order differential equation is developed by means of a shifted Jacobi-Gauss collocation spectral method. The shifted Jacobi-Gauss points are used as collocation nodes. Numerical examples are included to demonstrate the validity and applicability of the proposed algorithm, and some comparisons are made with the existing results. The method ...

properties of the hybrid of block-pulse functions and lagrange polynomials based on the legendre-gauss-type points are investigated and utilized to define the composite interpolation operator as an extension of the well-known legendre interpolation operator. the uniqueness and interpolating properties are discussed and the corresponding differentiation matrix is also introduced. the appli...

In this paper, we propose an iterative spectral method for solving differential equations with initial values on large intervals. In the proposed method, we first extend the Legendre wavelet suitable for large intervals, and then the Legendre-Guass collocation points of the Legendre wavelet are derived. Using this strategy, the iterative spectral method converts the differential equation to a s...

Solutions of partial differential equations with coordinate singularities often have special behavior near the singularities, which forces them to be smooth. Special treatment for these coordinate singularities is necessary in spectral approximations in order to avoid degradation of accuracy and efficiency. It has been observed numerically in the past that, for a scheme to attain high accuracy,...

Abstract. Existing analysis shows that when the Gauss Runge-Kutta (GRK) (also called Legendre-Gauss collocation) formulation with s Gaussian nodes is applied to ordinary differential equation initial value problems, the discretization has order 2s (super-convergent) [8]. However, for time-dependent partial differential equations (PDEs) with boundary conditions, super-convergence is only observe...

problem (OCP) for systems governed by Volterra integro-differential (VID) equation is considered. The method is developed by means of the Legendre wavelet approximation and collocation method. The properties of Legendre wavelet together with Gaussian integration method are utilized to reduce the problem to the solution of nonlinear programming one. Some numerical examples are given to confirm t...