Let ℓj be the Lagrange basispolynomials for the nodes. Necessity and existence give uniquely Aj=∫−11ℓj(t)dt, because every polynomial of degree at most n−1 is ∑jp(tj)ℓj(t).
For degp≤2n−1, divide p=qpn+r with degq≤n−1, degr≤n−1. Orthogonality makes ∫qpn=0, and p(tj)=r(tj). Apply part (i) to r. This is Gaussian quadrature.