Constructing Krall–Hahn orthogonal polynomials

Given a sequence of polynomials (pn)n(pn)n, an algebra of operators AA that acts in the linear space of polynomials and an operator Dp∈ADp∈A with Dp(pn)=θnpnDp(pn)=θnpn, where θnθn is any arbitrary eigenvalue, we construct a new sequence of polynomials (qn)n(qn)n by considering a linear combination of m+1m+1 consecutive pnpn: qn=pn+∑j=1mβn,jpn−j. Using the concept of a DD-operator, we determine the structure of the sequences βn,jβn,j, j=1,…,mj=1,…,m, such that the polynomials (qn)n(qn)n are eigenfunctions of an operator in the algebra AA. As an application, from the classical discrete family of Hahn polynomials, we construct orthogonal polynomials (qn)n(qn)n that are also eigenfunctions of higher-order difference operators.