Product of Chebyshev polynomials
Chebyshev polynomials satisfy a lot of identities, much like trig functions do. This point will look briefly at just one such identity.
Chebyshev polynomials Tn are defined for n = 0 and 1 by
T0(x) = 1
T1(x) = x
and for larger n using the recurrence relation
Tn+1(x) = 2xTn(x) - Tn-1(x)
This implies
T2(x) = 2xT1(x) - T0(x) = 2x2 - 1
T3(x) = 2xT2(x) - T1(x) = 4x3 - 3x
T4(x) = 2xT3(x) - T2(x) = 8x4 - 8x2 + 1
and so forth.
Now for the identity for this post. If m a n, then
2 TmTn = Tm+n + Tm-n.
In other words, the product of the mth and nth Chebyshev polynomials is the average of the (m + n)th and (m - n)th Chebyshev polynomials. For example,
2 T3(x) T1(x) = 2 (4x3 - 3x) x = T4(x) + T2(x)
The identity above is not at all apparent from the recursive definition of Chebyshev polynomials, but it follows quickly from the fact that
Tn(cos I) = cos nI.
Proof: Let I = arccos x. Then
2 Tm(x) Tn(x)
= 2 Tm(cos I) Tn(cos I)
= 2 cos mI cos nI
= cos (m+n)I + cos (m-n)I
= Tm+n(cos I) + Tm-n(cos I)
= Tm+n(x) + Tm-n(x)
You might object that this only shows that the first and last line are equal for values of x that are cosines of some angle, i.e. values of x in [-1, 1]. But if two polynomials agree on an interval, they agree everywhere. In fact, you don't need an entire interval. For polynomials of degree m+n, as above, it is enough that they agree on m + n + 1 points. (Along those lines, see Binomial coefficient trick.)
The close association between Chebyshev polynomials and cosines means you can often prove Chebyshev identities via trig identities as we did above.
Along those lines, we could have taken
Tn(cos I) = cos nI
as the definition of Chebyshev polynomials and then proved the recurrence relation above as a theorem, using trig identities in the proof.
Forman Acton suggested in this book Numerical Methods that Work that you should think of Chebyshev polynomials as "cosine curves with a somewhat disturbed horizontal scale."