Continued fractions and orthogonal polynomials in several variables

We extend the close interplay between continued fractions, orthogonal polynomials, and Gaussian quadrature rules to several variables in a special but natural setting which we characterize in terms of moment sequences. The crucial condition for the characterization is the commutativity of the multiplication operators on finite polynomial subspaces modulo an ideal. Moreover, starting from the orthogonal polynomials or the three-term recurrence, our method constructs a sequence of moment sequences that provides an approximation of the maximal order for recovering the moment sequence that defines the orthogonality.