 Solution

Write $W={\int }_{-1}^{1}w(x)dx$, finite and positive. Exactness on the constant polynomial gives ${\sum }_i{a}_i^{(n)}=W$. Positivity therefore implies the uniform operator bound $|I_n(h)|\le W\Vert h{\Vert }_{\infty }$, and the integral has the same bound.

For a continuous $f$ and $ϵ>0$, the Weierstrass approximation theorem gives a fixed polynomial $Q$ of degree $m$ with $\Vert f-Q{\Vert }_{\infty }<ϵ$. For every $n\ge m$, polynomial exactness gives $I_n(Q)=I(Q)$. ThusSince $ϵ$ is arbitrary, $\boxed{I_n(f)\to I(f)}$ for every continuous $f$. This is convergence of positive quadrature on continuous functions. The argument requires nodes in $[-1,1]$, implicit in a rule defined on all continuous functions there; otherwise their function values would not even be specified by the stated domain. Exactness alone would not provide the norm bound if arbitrarily signed weights were allowed.