 This method does not require any Newton-Raphson iterations or evaluations of Bessel functions as other methods do.As shown in the paper, the method was able to compute the nodes at a problem size of one billion in 11 seconds.Since the nodes near the endpoints of formula_4 become very close to each other at this problem size, this method of computing the nodes is sufficient for essentially any practical application in double-precision floating point.