Bivariate orthogonal polynomials in the Lyskova class

Classical orthogonal polynomials in two variables can be characterized as the polynomial solutions of a matrix second-order partial differential equation involving matrix polynomial coefficients. In this work, we study classical orthogonal polynomials in two variables whose partial derivatives satisfy again a second-order partial differential equation of the same type.