Quadrature and polynomial chaos

Builds Gauss-Legendre and Gauss-Hermite nodes and weights, and computes the mean, variance and moments of a polynomial from its expansion coefficients.