Sobolev orthogonality of polynomial solutions of second-order partial differential equations

Given a second-order partial differential operator L with nonzero polynomial coefficients of degree at most 2, and a Sobolev bilinear form ( P, Q) S = N i=0 i j=0 u(i, j),. i- j x. j y P. i- j x. j y Q , N 0, where u(i, j), 0 j i N, are linear functionals defined on the space of bivariate polynomials, we study the orthogonality of the polynomial solutions of the partial differential equation L[ p] =.n,m p with respect to (center dot, center dot) S, where.n,m are eigenvalue parameters depending on the total and partial degree of the solutions. We show that the linear functionals in the bilinear form must satisfy Pearson equations related to the coefficients ofL. Therefore, we also study solutions of the Pearson equations that can be obtained from univariatemoment functionals. In fact, the involved univariate functionals must satisfy Pearson equations in one variable. Moreover, we study polynomial solutions of L[ p] =.n, m p obtained from univariate sequences of polynomials satisfying second-order ordinary differential equations.