 Solution

Define the Lagrange interpolation polynomialsDistinct interpolation nodes ensure every denominator is nonzero and $L_j(x_i)={\delta }_{ij}$. For a polynomial of degree at most $n-1$, the difference $P(x)-{\sum }_{j}P(x_j)L_j(x)$ has degree at most $n-1$ and vanishes at $n$ distinct points, so it is zero. Integrating this Lagrange interpolation identity proves the required exact quadrature formula. This also covers $n=1$, when the empty product is one and its weight is two.