“…These are a generalization of orthogonal polynomials which arises naturally in Hermite-Padé approximation of a system of (Markov) functions [9,10,22]. In particular, they satisfy orthogonality conditions with respect to several positive measures [2,25,29]. Some of their applications are situated in diophantine number theory, rational approximation, spectral and scattering problems for higher-order difference equations and some associated dynamical systems; see, e.g., [5,11,18,27].…”