The Lagrange-mesh method is an approximate variational method taking the form of equations on a grid because of the use of a Gauss quadrature approximation. With a basis of Lagrange functions involving associated Laguerre polynomials related to the Gauss quadrature, the method is applied to the Dirac equation. The potential may possess a 1/r singularity. For hydrogenic atoms, numerically exact ...

This paper propose and analyze the two-level stabilized finite element method for the stationary Navier-Stokes equations based on Newon iteration. This algorithm involves solving one small, nonlinear coarse mesh with mesh size H and two linear problems on the fine mesh with mesh size h. Based on local Gauss integration and the quadratic equal-order triangular element ,the two-level stabilized m...

We analyze the consistency, stability, and convergence of an hp discontinuous Galerkin spectral element method. The analysis is carried out simultaneously for acoustic, elastic, coupled elastic–acoustic, and electromagnetic wave propagation. Our analytical results are developed for both conforming and non-conforming approximations on hexahedral meshes using either exact integration with Legendr...

Abstract. Using only fairly simple and elementary considerations essentially from first year undergraduate mathematics we show how the classical Stokes’ theorem for any given surface and vector field in R follows from an application of Gauss’ divergence theorem to a suitable modification of the vector field in a tubular shell around the given surface. The two stated classical theorems are (like...

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 construct and analyze Gauss-type quadrature rules with complex-valued nodes and weights to approximate oscillatory integrals with stationary points of high order. The method is based on substituting the original interval of integration by a set of contours in the complex plane, corresponding to the paths of steepest descent. Each of these line integrals shows an exponentially decaying behavi...

A software package is described for the numerical evaluation of integrals that arise in Bayesian statistical analysis. Diierent types of transformations that can be user selected to precondition the problem are discussed. These all begin with a standardizing transformation, that can be adaptively constructed, and is followed by a multivariate Normal, multivariate Student's t or split-t transfor...

A mesh refinement method is described for solving a continuous-time optimal control problem using collocation at Legendre-Gauss-Radau points. The method allows for changes in both the number of mesh intervals and the degree of the approximating polynomial within a mesh interval. First, a relative error estimate is derived based on the difference between the Lagrange polynomial approximation of ...

A mesh refinement method is described for solving a continuous-time optimal control problem using collocation at Legendre–Gauss–Radau points. The method allows for changes in both the number of mesh intervals and the degree of the approximating polynomial within a mesh interval. First, a relative error estimate is derived based on the difference between the Lagrange polynomial approximation of ...

This paper pursues obtaining Jacobi spectral collocation methods to solve Caputo fractional differential equations numerically. We used the shifted Jacobi–Gauss–Lobatto or Jacobi–Gauss–Radau quadrature nodes as points and derived differentiation matrices for derivatives. With matrices, were transformed into linear systems, which are easier solve. Two types of numerical simulations, results demo...