Solution (source code)

= Solution

Given $\epsilon>0$, the <Weierstrass approximation theorem> supplies a polynomial $Q$ with $\|f-Q\|_\infty<\epsilon$. For all sufficiently large $n$, the quadrature is exact on $Q$. Positivity of the weights and $\sum_jA_j=\int r$ then give
$$
\left|\int fr-\sum_jA_jf(\omega_j)\right|
\leq2\epsilon\int_{-1}^1r.
$$
Letting $\epsilon\downarrow0$ proves convergence.