Interpolation approximations based on Gauss–Lobatto–Legendre–Birkhoff quadrature

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 user-oriented implementation of the pseudospectral methods based on the GLLB quadrature nodes for Neumann problems. This approach allows an exact imposition of Neumann boundary conditions, and is as efficient as the pseudospectral methods based on Gauss–Lobatto quadrature for PDEs with Dirichlet boundary conditions.