The mass matrix for Gauss-Lobatto grid points is usually approximated by GaussLobatto quadrature because this leads to a diagonal matrix that is easy to invert. The exact mass matrix and its inverse are full. We show that the exact mass matrix and its inverse differ from the approximate diagonal ones by a simple rank-1 update (outer product). They can thus be applied to an arbitrary vector in O...

Abstract In this manuscript, we implement a spectral collocation method to find the solution of reaction–diffusion equation with some initial and boundary conditions. We approximate by using two-dimensional interpolating polynomial dependent Legendre–Gauss–Lobatto points. fully show that achieved solutions are convergent exact when number points increases. demonstrate capability efficiency prov...

We consider singularly perturbed linear boundary value problems for ODE's with variable coefficients, but without turning points. Convergence results are obtained for collocation schemes based on Gauss and Lobatto points, showing that highly accurate numerical solutions for these problems can be obtained at a very reasonable cost using such schemes, provided that appropriate meshes are used. Th...

We derive in this paper the asymptotic estimates of the nodes and weights of the Gauss-Lobatto-Legendre-Birkhoff (GLLB) quadrature formula, and obtain optimal error estimates for the associated GLLB interpolation in Jacobi weighted Sobolev spaces. We also present a useroriented implementation of the pseudospectral methods based on the GLLB quadrature nodes for Neumann problems. This approach al...

The Falkner-Skan equation is a nonlinear third-order boundary value problem defined on the semi-infinite interval [0,∞). This equation plays an important role to illustrate the main physical features of boundary layer phenomena. This paper presents a new collocation method for solving the Falkner-Skan equation. The proposed approach is equipped by the orthogonal Chebyshev polynomials that have ...

An algorithm is described for the efficient numerical solution of singular two-point boundary value problems. The algorithm is based on collocation at Gauss points, is applied directly to second order equations and uses a transformation of the independent variable to obtain extra smoothness if needed. Numerical comparisons between a code based on this approach and codes based on other options w...

in this paper, an effective technique is proposed to determine thenumerical solution of nonlinear volterra-fredholm integralequations (vfies) which is based on interpolation by the hybrid ofradial basis functions (rbfs) including both inverse multiquadrics(imqs), hyperbolic secant (sechs) and strictly positive definitefunctions. zeros of the shifted legendre polynomial are used asthe collocatio...

This paper develops a family of preconditioners for pseudospectral approximations of pth-order linear differential operators subject to various types of boundary conditions. The approximations are based on ultraspherical polynomials with special attention being paid to Legendre and Chebyshev polynomial methods based on Gauss–Lobatto quadrature points. The eigenvalue spectrum of the precondition...