 Solution

For node $x_i$, let ${\ell }_i$ be the Lagrange interpolation polynomial with ${\ell }_i(x_j)={\delta }_{ij}$. It has degree $n$, so ${\ell }_i^2$ has degree $2n$. The Gaussian quadrature with $n+1$ nodes is exact through degree $2n+1$, henceThe integrand is nonnegative. Since ${\ell }_i$ is a nonzero polynomial, it is nonzero on some open subinterval; the positive weight has positive integral there. Consequently $\boxed{a_i>0}$ for every node. This proves the positive weights of Gaussian quadrature property without assuming any coefficient sign beforehand.