In this article, a shifted Legendre tau method is introduced to get a direct solution technique for solving multi-order fractional differential equations (FDEs) with constant coefficients subject to multi-point boundary conditions. The fractional derivative is described in the Caputo sense. Also, this article reports a systematic quadrature tau method for numerically solving multi-point boundar...

Abstract. A local convergence rate is established for an orthogonal collocation method based on Radau quadrature applied to an unconstrained optimal control problem. If the continuous problem has a sufficiently smooth solution and the Hamiltonian satisfies a strong convexity condition, then the discrete problem possesses a local minimizer in a neighborhood of the continuous solution, and as the...

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

We describe the design and implementation of GAUSS Ð an online algorithm selection system for numerical quadrature. Given a quadrature problem and performance constraints on its solution, GAUSS selects the best (or nearly best) algorithm. GAUSS uses inductive logic programming to generalize a database of performance data; this produces high-level rules that correlate problem features with algor...

We consider the convergence of Gauss-type quadrature formulas for the integral ∫∞ 0 f(x)ω(x)dx, where ω is a weight function on the half line [0,∞). The n-point Gauss-type quadrature formulas are constructed such that they are exact in the set of Laurent polynomials Λ−p,q−1 = { ∑q−1 k=−p akx k}, where p = p(n) is a sequence of integers satisfying 0 ≤ p(n) ≤ 2n and q = q(n) = 2n − p(n). It is pr...

A generalized framework is presented that extends the classical theory of finite-difference summation-by-parts (SBP) operators to include a wide range of operators, where the main extensions are i) non-repeating interior point operators, ii) nonuniform nodal distribution in the computational domain, iii) operators that do not include one or both boundary nodes. Necessary and sufficient conditio...

Tensor products of Gauss-Lobatto quadrature points are frequently used as collocation points in spectral element methods. Unfortunately, it is not known if Gauss-Lobatto points exist in non-tensor-product domains like the simplex. In this work, we show that the n-dimensional tensor-product of Gauss-Lobatto quadrature points are also Fekete points. This suggests a way to generalize spectral meth...

In present paper, a numerical approach for solving Cauchy type singular integral equations is discussed. Lagrange interpolation with Gauss Legendre quadrature nodes and Taylor series expansion are utilized to reduce the computation of integral equations into some algebraic equations. Finally, five examples with exact solution are given to show efficiency and applicability of the method. Also, w...

This paper evaluates a mutual-information integral that appears in EXIT-chart analysis and the capacity equation for the binaryinput AWGN channel. Rapidly-converging series representations are derived and shown to be computationally efficient. Over a wide range of channel signal-to-noise ratios, the series are more accurate than the Gauss-Hermite quadrature and comparable to Monte Carlo integra...