PHYSICS · Multi-part · Fermi University Physics · Mathematical Methods for Physics · Special functions and orthogonal polynomials · Numerical and computational mathematics
Error, Chance and the Machine, Ch Integration: error expansions, Gauss rules and adaptivity
Problem
The Legendre polynomial of degree 3 is P3(x)=21(5x3−3x), and the three-point Clenshaw-Curtis nodes on [−1,1] are xk=cos(kπ/2) for k=0,1,2.
((a)) Find the three roots of P3 (one of them is x=0) and take them as the nodes of a three-node Gauss-Legendre rule. Using the symmetry x1=−x3, x2=0, write its weights as w1=w3=a and w2=b, and find a and b by demanding exactness for f(x)=1 and f(x)=x2.
((b)) Check that the same a,b also make the rule exact for f(x)=x4, and state the degree of exactness this three-node rule reaches.
((c)) Evaluate the three Clenshaw-Curtis nodes xk=cos(kπ/2), and find their weights p (at x=±1) and q (at x=0) from exactness for f(x)=1 and f(x)=x2. Identify the classical rule this reproduces.
((d)) Test the Clenshaw-Curtis rule of part (c) on f(x)=x4 and state its degree of exactness. Compare, at the same cost of three function evaluations, the degrees of exactness reached by the Gauss and Clenshaw-Curtis rules.