Faster and more accurate Gauss-Legendre nodes and weights
This algorithm was pointed out to me by Sigurd Naess.
For degrees <=100 I'm still using the old one (however now working with long double
instead of double
internally). This is also quick enough and avoids storing huge coefficient tables.