The H function can be approximated up to an order  as
where  are the zeros of Legendre polynomials  and  are the positive, non vanishing roots of the associated characteristic equation
where  are the quadrature weights given by