MathLabs

Analysis

Power series and Taylor series

Infinite polynomial-like sums that represent functions locally, such as eˣ = 1+x+x²/2!+…

IntuitionA polynomial that never stops

Some functions are hard to compute directly — square roots, exponentials, sines — but if you let a polynomial grow to infinitely many terms, it can match the function exactly near a point. The classic example is ex=∑n=0∞xnn!\displaystyle e^x = \sum_{n=0}^{\infty} \frac{x^n}{n!}: add up enough terms and you get a value indistinguishable from the true exponential.

Cubic polynomial graph approximating the exponential function near zero, adjustable coefficients.
This cubic 1+x+x22+x361 + x + \dfrac{x^2}{2} + \dfrac{x^3}{6} is the degree-3 Taylor polynomial of exe^x at x=0x=0. Compare its graph with exe^x: they nearly overlap near x=0x=0 and drift apart farther away.

SchoolWriting a function as an infinite sum

Definition: Power series

A power series centered at cc is an expression ∑n=0∞an(x−c)n\displaystyle\sum_{n=0}^{\infty} a_n (x-c)^n, where ana_n are fixed real coefficients and xx is the variable. It is really a "polynomial with infinitely many terms."

∑n=0∞an(x−c)n=a0+a1(x−c)+a2(x−c)2+⋯\sum_{n=0}^{\infty} a_n (x-c)^n = a_0 + a_1(x-c) + a_2(x-c)^2 + \cdots

A power series does not converge for every xx: it always converges on some interval centered at cc, called the interval of convergence, whose half-length is the radius of convergence RR. When xx is inside this interval, the infinite sum adds up to a finite number; when xx is outside, the terms blow up and the sum has no meaning.

an=f(n)(c)n!,f(x)=∑n=0∞f(n)(c)n!(x−c)na_n = \frac{f^{(n)}(c)}{n!}, \qquad f(x) = \sum_{n=0}^{\infty} \frac{f^{(n)}(c)}{n!}(x-c)^n
Behavior by position relative to R
PositionBehavior
∣x−c∣<R|x-c| < RAbsolute convergence
∣x−c∣>R|x-c| > RDivergence
∣x−c∣=R|x-c| = RMust check endpoints separately

UndergraduateTheorems: radius of convergence and term-by-term differentiation

For the power series ∑n=0∞an(x−c)n\displaystyle\sum_{n=0}^{\infty} a_n (x-c)^n, suppose L=lim⁡n→∞∣an+1an∣L = \displaystyle\lim_{n\to\infty} \left| \frac{a_{n+1}}{a_n} \right| exists (as a finite number or +∞+\infty). Then the radius of convergence is R=1LR = \dfrac{1}{L}, with the convention that R=∞R = \infty if L=0L=0 and R=0R = 0 if L=∞L=\infty.

Why is it true?

The ratio test compares consecutive terms of the series to a geometric series: if the ratio of consecutive terms is eventually less than 1 in absolute value, the series behaves like a convergent geometric series and adds up to a finite number.

Proof

Fix x≠cx \neq c and apply the ratio test to the terms bn=an(x−c)nb_n = a_n(x-c)^n of the numerical series ∑bn\sum b_n. We compute

lim⁡n→∞∣bn+1bn∣=lim⁡n→∞∣an+1(x−c)n+1an(x−c)n∣=lim⁡n→∞∣an+1an∣⋅∣x−c∣=L ∣x−c∣\displaystyle \lim_{n\to\infty} \left| \frac{b_{n+1}}{b_n} \right| = \lim_{n\to\infty} \left| \frac{a_{n+1}(x-c)^{n+1}}{a_n(x-c)^n} \right| = \lim_{n\to\infty} \left| \frac{a_{n+1}}{a_n} \right| \cdot |x-c| = L\,|x-c|.

By the ratio test for numerical series, ∑bn\sum b_n converges absolutely when L∣x−c∣<1L|x-c| < 1, i.e. when ∣x−c∣<1/L|x-c| < 1/L, and diverges when L∣x−c∣>1L|x-c| > 1, i.e. when ∣x−c∣>1/L|x-c| > 1/L. So the series converges absolutely for every xx with ∣x−c∣<1/L|x-c| < 1/L and diverges for every xx with ∣x−c∣>1/L|x-c| > 1/L.

This is exactly the definition of the radius of convergence: the largest RR such that the series converges for all ∣x−c∣<R|x-c| < R. Hence R=1/LR = 1/L. The boundary cases L=0L=0 (ratio always shrinks to 0, so the series converges for every xx, giving R=∞R=\infty) and L=∞L=\infty (ratio blows up for any x≠cx \neq c, so the series converges only at x=cx=c, giving R=0R=0) follow the same argument taking the appropriate limits.

If ∑n=0∞an(x−c)n\displaystyle\sum_{n=0}^{\infty} a_n (x-c)^n has radius of convergence R>0R>0, then f(x)=∑n=0∞an(x−c)nf(x) = \sum_{n=0}^{\infty} a_n(x-c)^n is differentiable on (c−R, c+R)(c-R,\ c+R), and its derivative can be computed term by term: f′(x)=∑n=1∞nan(x−c)n−1\displaystyle f'(x) = \sum_{n=1}^{\infty} n a_n (x-c)^{n-1}, with the same radius of convergence RR.

Why is it true?

A power series behaves like an infinite polynomial, and polynomials can be differentiated term by term; the theorem says this familiar rule survives the passage to infinitely many terms, as long as we stay strictly inside the radius of convergence.

Proof

First we show the differentiated series has the same radius of convergence. Let cn=nanc_n = n a_n be the coefficients of the term-by-term derivative (shifted index). Since n1/n→1n^{1/n} \to 1 as n→∞n \to \infty, we have lim sup⁡n∣nan∣1/n=lim sup⁡n∣an∣1/n\limsup_n |n a_n|^{1/n} = \limsup_n |a_n|^{1/n}, and by the Cauchy–Hadamard formula the two series ∑an(x−c)n\sum a_n(x-c)^n and ∑nan(x−c)n−1\sum n a_n (x-c)^{n-1} have the same radius of convergence RR.

Next, fix any closed subinterval [c−ρ,c+ρ][c-\rho, c+\rho] with ρ<R\rho < R. On this subinterval the series of derivatives ∑nan(x−c)n−1\sum n a_n (x-c)^{n-1} converges uniformly, because its terms are dominated (for nn large) by Mρn−1M \rho^{n-1} for constants coming from a slightly larger radius ρ′∈(ρ,R)\rho' \in (\rho, R), and the Weierstrass M-test applies.

Uniform convergence of the derivative series on [c−ρ,c+ρ][c-\rho, c+\rho] together with pointwise convergence of the original series ∑an(x−c)n\sum a_n(x-c)^n lets us invoke the standard theorem on differentiating a series of functions term by term: the sum f(x)f(x) is differentiable on [c−ρ,c+ρ][c-\rho,c+\rho] and f′(x)f'(x) equals the sum of the derivative series there. Since ρ<R\rho < R was arbitrary, this holds on all of (c−R, c+R)(c-R,\ c+R).

UndergraduateReal-World Applications and Worked Examples

Engineers and scientists rarely evaluate sin⁡\sin, cos⁡\cos, exe^x or  \sqrt{\ } by hand; calculators and computer chips use truncated power series (with a controlled error term) to compute them. Physicists use the small-angle approximation sin⁡x≈x−x36\sin x \approx x - \dfrac{x^3}{6} to linearize pendulum and oscillator equations, and financial models use exponential series to approximate continuous compounding over short time steps.

Example: Radius and interval of convergence with an asymmetric endpoint

Find the radius and interval of convergence of ∑n=1∞(−1)nn⋅4n(x−2)n\displaystyle\sum_{n=1}^{\infty} \frac{(-1)^n}{n\cdot 4^n}(x-2)^n.

Solution

Here c=2c=2 and an=(−1)nn⋅4na_n = \frac{(-1)^n}{n\cdot 4^n}. By the ratio test theorem above, L=lim⁡n→∞∣an+1an∣=lim⁡n→∞n4(n+1)=14L = \lim_{n\to\infty} \left| \frac{a_{n+1}}{a_n} \right| = \lim_{n\to\infty} \frac{n}{4(n+1)} = \frac{1}{4}, so R=4R = 4.

This gives the open interval (−2,6)(-2,6) where the series converges absolutely. We must check the two endpoints separately, since the ratio test is inconclusive there.

At x=−2x=-2: (x−2)n=(−4)n(x-2)^n = (-4)^n, so the term becomes (−1)nn⋅4n(−4)n=(−1)n(−1)n4nn⋅4n=1n\frac{(-1)^n}{n\cdot 4^n}(-4)^n = \frac{(-1)^n(-1)^n 4^n}{n \cdot 4^n} = \frac{1}{n}, giving the harmonic series ∑n=1∞1n\displaystyle\sum_{n=1}^{\infty} \frac{1}{n}, which diverges.

At x=6x=6: (x−2)n=4n(x-2)^n = 4^n, so the term becomes (−1)nn⋅4n⋅4n=(−1)nn\frac{(-1)^n}{n\cdot 4^n}\cdot 4^n = \frac{(-1)^n}{n}, the alternating harmonic series ∑n=1∞(−1)nn\displaystyle\sum_{n=1}^{\infty} \frac{(-1)^n}{n}, which converges by the alternating series test.

So the interval of convergence is (−2, 6](-2,\ 6]: closed on the right, open on the left.

Example: Estimating cos(0.2) for a pendulum correction

An engineer needs cos⁡(0.2)\cos(0.2) (radians) to a precision of 10−610^{-6} when modeling the restoring torque of a pendulum. Use the Maclaurin series of cos⁡x\cos x truncated at the x4x^4 term to estimate it, and bound the error.

Solution

The Maclaurin series of cos⁡x\cos x is cos⁡x=∑n=0∞(−1)nx2n(2n)!=1−x22+x424−⋯\cos x = \sum_{n=0}^{\infty} \frac{(-1)^n x^{2n}}{(2n)!} = 1 - \frac{x^2}{2} + \frac{x^4}{24} - \cdots, a special case of f(x)=∑n=0∞f(n)(c)n!(x−c)n\displaystyle f(x) = \sum_{n=0}^{\infty} \frac{f^{(n)}(c)}{n!}(x-c)^n with c=0c=0.

Truncating at the x4x^4 term with x=0.2x=0.2: 1−(0.2)22+(0.2)424=1−0.02+0.0000667≈0.9800671 - \dfrac{(0.2)^2}{2} + \dfrac{(0.2)^4}{24} = 1 - 0.02 + 0.0000667 \approx 0.980067.

Because cos⁡x\cos x's Maclaurin series is alternating with decreasing terms for small xx, the error is bounded by the first omitted term: (0.2)6720≈8.9×10−8\dfrac{(0.2)^6}{720} \approx 8.9\times 10^{-8}, which is already far below the required 10−610^{-6} tolerance.

So the engineer can safely use ≈0.980067\approx 0.980067 instead of evaluating a transcendental cosine directly — exactly what a calculator's internal algorithm does.

For a power series with L=lim⁡n→∞∣an+1/an∣=5L = \lim_{n\to\infty} |a_{n+1}/a_n| = 5, what is the radius of convergence?

What does the Taylor series formula f(x)=∑n=0∞f(n)(c)n!(x−c)n\displaystyle f(x) = \sum_{n=0}^{\infty} \frac{f^{(n)}(c)}{n!}(x-c)^n require of ff at cc?

If ∑n=0∞an(x−c)n\displaystyle\sum_{n=0}^{\infty} a_n (x-c)^n has radius of convergence RR, what is the radius of convergence of the differentiated series f′(x)=∑n=1∞nan(x−c)n−1\displaystyle f'(x) = \sum_{n=1}^{\infty} n a_n (x-c)^{n-1}?

Using degree-4 Maclaurin truncation, engineers approximate cos⁡(0.2)\cos(0.2) with error bounded by which quantity?