Abstract:The aim of this paper is to describe a Matlab package for computing the simultaneous Gaussian quadrature rules associated with a variety of multiple orthogonal polynomials. Multiple orthogonal polynomials can be considered as a generalization of classical orthogonal polynomials, satisfying orthogonality constraints with respect to $r$ different measures, with $r \ge 1$. Moreover, they satisfy $(r+2)$--term recurrence relations. In this manuscript, without loss of generality, $r$ is considered equal to $2$. The so-called simultaneous Gaussian quadrature rules associated with multiple orthogonal polynomials can be computed by solving a banded lower Hessenberg eigenvalue problem. Unfortunately, computing the eigendecomposition of such a matrix turns out to be strongly ill-conditioned and the \texttt{Matlab} function \texttt{balance.m} does not improve the condition of the eigenvalue problem. Therefore, most procedures for computing simultaneous Gaussian quadrature rules are implemented with variable precision arithmetic. Here, we propose a \texttt{Matlab} package that allows to reliably compute the simultaneous Gaussian quadrature rules in floating point arithmetic. It makes use of a variant of a new balancing procedure, recently developed by the authors of the present manuscript, that drastically reduces the condition of the Hessenberg eigenvalue problem.

A Golub-Welsch version for simultaneous Gaussian quadrature

A Golub-Welsch version for simultaneous Gaussian quadrature

Computational aspects of simultaneous Gaussian quadrature

A MATLAB package computing simultaneous Gaussian quadrature rules for Multiple Orthogonal Polynomials

On Lebesgue Integral Quadrature

When Lanczos Iterations Generate Symmetric Quadrature Nodes?

Gaussian quadrature rules for composite highly oscillatory integrals

Complementary Romanovski-Routh polynomials and their zeros

Error bounds for Gauss–Lobatto quadrature of analytic functions on an ellipse

A note on Laguerre truncated polynomials and quadrature formula

Complex generalized Gauss-Radau quadrature rules for Hankel transforms of integer order

Positive quadrature formulas on an arc of the unit circle

Extended homogeneous bivariate orthogonal polynomials: symbolic and numerical Gaussian cubature formula

Gram quadrature: numerical integration with Gram polynomials

A Riemann-Hilbert Approach to Skew-Orthogonal Polynomials of Symplectic Type

Applications of multiple orthogonal polynomials with hypergeometric moment generating functions

Dual forms of the orthogonality relations of some classical q-orthogonal polynomials

A Note on the Gauss-Jacobi Quadrature Formulae for Singular Integral Equations of the Second Kind

An evolution of matrix-valued orthogonal polynomials

A Riemann--Hilbert Approach to Skew-Orthogonal Polynomials of Symplectic Type

Error bounds for Gauss–Jacobi quadrature of analytic functions on an ellipse