E. Heine in the 19th century studied a system of orthogonal polynomials associated with the weight [x(x-α)(x-β)]-12, x∈[0, α], 0<α<β. A related system was studied by C. J. Rees in 1945, associated with the weight [(1-x2)(1-k2x2)]-12, x∈[-1, 1], k∈(0, 1). These are also known as elliptic orthogonal polynomials, since the moments of the weights may be expressed in terms of elliptic integrals. Such orthogonal polynomials are of great interest because the corresponding Hankel determinant, depending on a parameter k, where 0-1,β∈R, satisfy second order non-linear difference equations. The large n expansion based on the difference equations when combined with known asymptotics of the leading terms of the associated Hankel determinant yields a complete asymptotic expansion of the Hankel determinant. The Painlevé equation is also discussed as well as the generalization of the linear second order differential equation found by Rees.