Solution (source code)

= Solution
{c}

At either quadrature node, $2x_j^2=1$, so $u(x_j)=\arctan1=\pi/4$. Thus the <Gaussian quadrature> approximation is
$$
\boxed{\frac\pi2\frac\pi4+\frac\pi2\frac\pi4=\frac{\pi^2}{4}.}
$$
This is the requested quadrature value; the integrand is not a cubic polynomial, so the degree-three exactness established in part (b) does not make this the exact integral.