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.

Classical multiple orthogonal polynomials for arbitrary number of weights and their explicit representation

Integral and hypergeometric representations for multiple orthogonal polynomials

Classical discrete multiple orthogonal polynomials: hypergeometric and integral representations

Toda and Laguerre-Freud equations for multiple discrete orthogonal polynomials with an arbitrary number of weights

Multiple Orthogonal Polynomials Associated with an Exponential Cubic Weight

Toda and Laguerre-Freud equations and tau functions for hypergeometric discrete multiple orthogonal polynomials

On the -Classical Orthogonal Polynomials

Toda and Laguerre–Freud equations and tau functions for hypergeometric discrete multiple orthogonal polynomials

Applications of multiple orthogonal polynomials with hypergeometric moment generating functions

Characterization of classical orthogonal polynomials in two variables

On two families of discrete multiple orthogonal polynomials on a star-like set

Discrete orthogonality relations for multi-indexed Laguerre and Jacobi polynomials

Real roots of random orthogonal polynomials with exponential weights

On a New Family of Weighted Least-square Orthogonal Polynomials in Multivariables

Second order differential operators having several families of orthogonal matrix polynomials as eigenfunctions

Mixed Type Multiple Orthogonal Polynomials Associated with the Modified Bessel Functions and Products of Two Coupled Random Matrices

Bidiagonal matrix factorisations related to multiple orthogonal polynomials

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

R$_{II}$ type three term relations for bivariate polynomials orthogonal with respect to varying weights

Laurent Multiple Orthogonal Polynomials on the Unit Circle

On Certain Wronskians of Multiple Orthogonal Polynomials.