Solution
= Solution
Let $\ell_i$ be the degree-$n-1$ cardinal polynomial with $\ell_i(c_j)=\delta_{ij}$. Exactness applied to the nonnegative polynomial $\ell_i^2$ gives
$$
b_i=\sum_jb_j\ell_i(c_j)^2
=\int_a^bw(x)\ell_i(x)^2\,dx>0,
$$
because $w>0$ and $\ell_i$ is nonzero. Thus every <Gaussian quadrature> weight is positive.