Weighted Grid Points

Construct a symbolic nonnegative weighting function w(x) on the interval [-1,1] and this script constructs the family of polynomials which are orthogonal with respect to the weight. Then the extrema of the highest polynomial are computed. Examples: w=1*x^0 gives the Legendre Lobatto points and w=1/sqrt(1-x^2) gives the Chebyshev Lobatto points. Useful for pseudospectral/collocation methods.