Weights in Golub-Welsch

Compute the quadrature weights by full Golub-Welsch algorithm.

Reason: for large n the explicit formula for computing the weights will fail due to the division by a polynomial of high degree.

Add this to the docstring.