Generating orthogonal matrix polynomials satisfying second order differential equations from a trio of triangular matrices (Q1040857)

It is known that the construction of orthogonal matrix polynomials satisfying a second order differential equation is closely connected with the symmetry of a second order differential operator. The operator has differential coefficients which are matrix polynomials \(F_2,F_1\) and \(F_0\) of degrees not bigger than 2,1 and 0 respectively. Sufficient conditions for the symmetry include a matrix first order differential equation that is equivalent to the fact that the weight matrix \(W\) can be factorized in the form \(W(t)=\rho(t)T(t)T^*(t),\) where \(\rho\) is a scalar function and \(T\) is a matrix function which satisfies a first order differential equation of the form \(T'=FT\). When \(\rho\) is one of the classical weights (Hermite, Laguerre or Jacobi) then \(F\) depends on two matrices \(A,B.\) The author provides a new method to find examples of symmetric operators when \(A,B\) and \(F_0\) are simultaneously triangularizable.