Abstract. We consider the non-stationary Heun equation, also known as quantum Painlevé VI, which has appeared in different works on quantum integrable models and conformal field theory. We use a generalized kernel function identity to transform the problem to solve this equation into a differential-difference equation which, as we show, can be solved by efficient recursive algorithms. We thus obtain series representations of solutions which provide elliptic generalizations of the Jacobi polynomials. These series reproduce, in a limiting case, a perturbative solution of the Heun equation due to Takemura, but our method is different in that we expand in non-conventional basis functions that allow us to obtain explicit formulas to all orders; in particular, for special parameter values, our series reduce to a single term.