The Heun operator of Hahn type

The Heun-Hahn operator on the uniform grid is defined. This operator is shown to map polynomials of degree $n$ to polynomials of degree $n+1$, to be tridiagonal in bases made out of either Pochhammer or Hahn polynomials and to be bilinear in the operators of the Hahn algebra. The extension of this algebra that includes the Heun-Hahn operator as generator is described. Biorthogonal rational functions on uniform grids are shown to be related to this framework.