We describe a technique to analytically compute the multipole moments of a charge distribution confined to a planar triangle, which may be useful in solving the Laplace equation using the fast multipole boundary element method (FMBEM) and for charged particle tracking. This algorithm proceeds by performing the necessary integration recursively within a specific coordinate system, and then trans...

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...

A new plastic hinge integration method overcomes the problems with nonobjective response caused by strain-softening behavior in force-based beam–column finite elements. The integration method uses the common concept of a plastic hinge length in a numerically consistent manner. The method, derived from the Gauss–Radau quadrature rule, integrates deformations over specified plastic hinge lengths ...

Gauss–Legendre quadrature formulas have excellent convergence properties when applied to integrals ∫ 1 0 f(x) dx with f ∈ C∞[0, 1]. However, their performance deteriorates when the integrands f(x) are in C∞(0, 1) but are singular at x = 0 and/or x = 1. One way of improving the performance of Gauss–Legendre quadrature in such cases is by combining it with a suitable variable transformation such ...

This paper reports numerical analyses of free vibration of laminated composite plate/ shell structures of various shapes, span-to-thickness ratios, boundary conditions and lay-up sequences. The method is based on a novel four-node quadrilateral element, namely MISQ20, within the framework of the first-order shear deformation theory (FSDT). The element is built by incorporating a strain smoothin...

Several algorithms are given and compared for computing Gauss quadrature rules. It is shown that given the three term recurrence relation for the orthogonal polynomials generated by the weight function, the quadrature rule may be generated by computing the eigenvalues and first component of the orthornormalized eigenvectors of a symmetric tridiagonal matrix. An algorithm is also presented for c...

Integration methods for the calculation of the magnetostatic field due to coils Marzia Fontana Abstract Analytical and numerical integration methods of the Biot-Savart law are presented for the three-dimensional computation of the magnetic eld generated by stationary currents. For complex-shaped conductors a nite volume approximation is proposed, based on a composite Gauss-Legendre quadrature o...

In recent years a number of authors have considered an error analysis for quadrature rules of Newton-Cotes type. In particular, the mid-point, trapezoid and Simpson rules have been investigated more recently ([2], [4], [5], [6], [11]) with the view of obtaining bounds on the quadrature rule in terms of a variety of norms involving, at most, the first derivative. In the mentioned papers explicit...

In this correspondence, the authors develop a novel method based on spherical radial cubature and Gauss–Laguerre quadrature rule for non-linear state estimation problems. The proposed filter, referred as cubature quadrature Kalman filter (CQKF) would be able to overcome inherent disadvantages associated with the earlier reported cubature Kalman filter (CKF). The theory and formulation of CQKF h...

The variable-step Störmer-Cowell integrator is a non-summed, double-integration multi-step integrator derived in variable-step form. The method has been implemented with a Shampine-Gordon style error control algorithm that uses an approximation of the local error at each step to choose the step size for the subsequent step. In this paper, the variable-step Störmer-Cowell method is compared to s...