Classical orthogonal polynomials in two variables: a matrix approach

Abstract Classical orthogonal polynomials in two variables are defined as the orthogonal polynomials associated to a two-variable moment functional satisfying a matrix analogue of the Pearson differential equation. Furthermore, we characterize classical orthogonal polynomials in two variables as the polynomial solutions of a matrix second order partial differential equation.