text: 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