“…, N. Equation (4) was obtained by employing two separate (but symmetrical) quadrature schemes of order N in the negative and positive ranges of the angular variable ξ . The quadrature nodes and weights for integration in ξ ∈ (0, ∞) are denoted, respectively, by {ξ n } and {w n }, and are defined as in Garcia (2013) by mapping the standard Gauss-Legendre nodes {μ n } and weights {v n } from (−1, 1) into (0, ∞), according to a transformation suggested by Gautschi (1968), viz.…”