This paper presents a Gaussian quadrature method for the evaluation of the triple integral ( , , ) T I f x y z d xd yd z = ∫∫∫ , where ) , , ( z y x f is an analytic function in , , x y z and T refers to the standard tetrahedral region:{( , , ) 0 , , 1, 1} x y z x y z x y z ≤ ≤ + + ≤ in three space( , , ). x y z Mathematical transformation from ( , , ) x y z space to ( , , ) u v w space maps th...

The use of Gaussian quadrature formulae is explored for the computation of the Macdonald function Kν(x) = ∫ ∞ 0 e cosh t cosh νt dt when x > 0 and ν is complex, ν = α+ iβ. It is shown that Gaussian quadrature with weight function w(t) = exp(−et) on [0,∞] is a viable approach, unless x is small and/or β large, but in combination with Gauss–Legendre quadrature, even in these latter cases.

Abstract We consider solitary-wave solutions of the generalized Burger’s-Fisher equation ∂Ψ ∂t + αΨ δ ∂Ψ ∂x − ∂ 2Ψ ∂x2 = βΨ(1 − Ψδ). In this paper, we present a new method for solving of the generalized Burger’s-Fisher equation by using the collocation formula for calculating spectral differentiation matrix for Chebyshev-Gauss-Lobatto point. To reduce round-off error in spectral collocation met...

Consider a hermitian positive-definite linear functional F, and assume we have m distinct nodes fixed in advance anywhere on the real line. In this paper we then study the existence and construction of nth rational Gauss-Radau (m = 1) and Gauss-Lobatto (m = 2) quadrature formulas that approximate F{f}. These are quadrature formulas with n positive weights and with the n−m remaining nodes real a...

High-order triangle-based discontinuous Galerkin (DG) methods for hyperbolic equations on a rotating sphere are presented. The DG method can be characterized as the fusion of finite elements with finite volumes. This DG formulation uses high-order Lagrange polynomials on the triangle using nodal sets up to 15th order. The finite element-type area integrals are evaluated using order 2N Gauss cub...

Parallel implementation issues of the explicit Verlet{type (VT) and the implicit Runge{Kutta (RK) methods for molecular dynamics integration (MD) on parallel distributed memory processors with the ring topology are described and compared. Both methods are applied to a system of N particles interacting through the Lennard{Jones potential. The RK method is carried out for the two{stage fourth{ord...

1. Introduction. The estimation of quadrature errors for analytic functions has been considered by Davis and Rabinowitz [1]. An estimate for the error of the Gaussian quadrature formula for analytic functions was obtained by Davis [2]. McNamee [3] has also discussed the estimation of error of the Gauss-Legendre quadrature for analytic functions. Convergence of the Gaussian quadratures was discu...

An asymptotic error expansion for Gauss-Legendre quadrature is derived for an integrand with an endpoint singularity. It permits convergence acceleration by extrapolation.