Chapter 8
Integration Using Chebyshev Polynomials
In this chapter we show how Chebyshev polynomials and some of their funda- mental properties can be made to play an important part in two key techniques of numerical integration.
• Gaussian quadrature estimates an integral by combining values of the integrand at zeros of orthogonal polynomials. We consider the special case of Gauss–Chebyshev quadrature, where particularly simple proce- dures follow for suitably weighted integrands. • One can approximately integrate a function by expanding it in a series and then integrating a partial sum of the series. We show that, for Chebyshev expansions, this process — essentially the Clenshaw–Curtis method — is readily analysed and again provides a natural procedure for appropriately weighted integrands.
Although this could be viewed as an ‘applications’ chapter, which in an introductory sense it certainly is, our aim here is primarily to derive further basic properties of Chebyshev polynomials.
8.1 Indefinite integration with Chebyshev series
If we wish to approximate the indefinite integral X h(X)= w(x)f(x)dx, −1 where −1
Suppose, in particular, that the weight w(x) is one of the four functions 1 1 1 w(x)=√ , 1, √ , √ , (8.2) 1 − x2 1 − x 1+x
and that we take fn(x) as the partial sum of the expansion of f(x) in Cheb- yshev polynomials of the corresponding one of the four kinds
Pk(x)=Tk(x),Uk(x),Vk(x),Wk(x). (8.3)
© 2003 by CRC Press LLC Then we can use the fact that (excluding the case where Pk(x)=Tk(x)with k =0) X w(x)Pk (x)dx = Ck(X)Qk(X) − Ck(−1)Qk(−1) −1 where Qk(X)=Uk−1(X),Tk+1(X),Wk(X),Vk(X) (8.4a) and √ √ √ 1 − X2 1 1 − X 1+X Ck(X)=− , , 2 1 , −2 1 , (8.4b) k k +1 k + 2 k + 2
respectively. (Note that Ck(−1) = 0 in the first and fourth cases.) This follows immediately from the fact that if x =cosθ then we have d k cos kθ sin kθ = − , dx sin θ d (k +1)sin(k +1)θ cos(k +1)θ = , dx sin θ d (k + 1 )cos(k + 1 )θ sin(k + 1 )θ = − 2 2 , dx 2 sin θ d (k + 1 )sin(k + 1 )θ cos(k + 1 )θ = 2 2 . dx 2 sin θ In the excluded case, we use d 1 θ = − dx sin θ to give X 1 √ 0 − − − 2 T (x)dx = arccos( 1) arccos X = π arccos X. −1 1 − x
Thus, for each of the weight functions (8.2) we are able to integrate the weighted polynomial and obtain the approximation hn(X) explicitly. Suppose that n n fn(x)= akTk(x)[Pk = Tk]or akPk(x)[Pk = Uk,Vk,Wk]. (8.5) k=0 k=0 Then in the first case n X hn(X)= ak w(x)Tk(x)dx = k=0 −1 √ n 2 1 1 − X = 2 a0(π − arccos X) − ak Uk−1(X), (8.6) k=1 k
© 2003 by CRC Press LLC while in the second, third and fourth cases