Power Series and Taylor Expansions

A structured guide to representing functions with power series, determining convergence, estimating approximation error, integrating series, and solving differential equations by coefficient matching.

The Basic Structure of

A centered at aa has the form

∑n=0∞cn(x−a)n=c0+c1(x−a)+c2(x−a)2+⋯ .\sum_{n=0}^{\infty}c_n(x-a)^n=c_0+c_1(x-a)+c_2(x-a)^2+\cdots.

When a=0a=0, the series is a Maclaurin . For a fixed value of xx, the expression becomes an ordinary infinite series of numbers, so convergence depends on xx.

The basic model is the geometric series:

∑n=0∞rn=11−r,∣r∣<1.\sum_{n=0}^{\infty}r^n=\frac{1}{1-r},\qquad |r|<1.

Replacing rr with x−ax-a gives

∑n=0∞(x−a)n=11−(x−a),∣x−a∣<1.\sum_{n=0}^{\infty}(x-a)^n=\frac{1}{1-(x-a)},\qquad |x-a|<1.

This identity is the main starting point for constructing many other series.

Takeaway: Treat a as an infinite polynomial, but always determine where it converges before using it to represent a function.

Convergence and Endpoint Testing

For a series centered at aa, the RR describes the interior region where convergence is guaranteed:

∣x−a∣<R.|x-a|<R.

There are three possibilities:

  • R=0R=0: the series converges only at x=ax=a.

  • R=∞R=\infty: the series converges for every real xx.

  • 0<R<∞0<R<\infty: the series converges inside the interval and diverges outside it.

To find RR, apply the ratio test to consecutive terms. For

∑n=0∞cn(x−a)n,\sum_{n=0}^{\infty}c_n(x-a)^n,

consider

L=lim⁡n→∞∣cn+1(x−a)n+1cn(x−a)n∣.L=\lim_{n\to\infty}\left|\frac{c_{n+1}(x-a)^{n+1}}{c_n(x-a)^n}\right|.

The ratio test gives convergence when L<1L<1 and divergence when L>1L>1. Solve the resulting inequality, usually of the form ∣x−a∣<R|x-a|<R, and then test x=a−Rx=a-R and x=a+Rx=a+R separately.

For example, consider

∑n=0∞(x−2)n(n+1)3n.\sum_{n=0}^{\infty}\frac{(x-2)^n}{(n+1)3^n}.

The ratio limit is ∣x−2∣/3|x-2|/3, so the interior condition is ∣x−2∣<3|x-2|<3, or −1<x<5-1<x<5. At x=−1x=-1, the resulting alternating harmonic series converges; at x=5x=5, the harmonic series diverges. Therefore,

R=3,[−1,5).R=3,\qquad [-1,5).

Takeaway: The ratio test usually determines only the open interval. Endpoint tests determine whether the final interval uses parentheses or brackets.

Differentiating and Integrating Series

Inside the , a can be differentiated and integrated term by term. If

f(x)=∑n=0∞cn(x−a)n,f(x)=\sum_{n=0}^{\infty}c_n(x-a)^n,

then

f′(x)=∑n=1∞ncn(x−a)n−1f'(x)=\sum_{n=1}^{\infty}nc_n(x-a)^{n-1}

and

∫f(x) dx=C+∑n=0∞cnn+1(x−a)n+1.\int f(x)\,dx=C+\sum_{n=0}^{\infty}\frac{c_n}{n+1}(x-a)^{n+1}.

The differentiated and integrated series have the same as the original series, although endpoint behavior can change.

Starting with

11−x=∑n=0∞xn,∣x∣<1,\frac{1}{1-x}=\sum_{n=0}^{\infty}x^n, \qquad |x|<1,

differentiation gives

1(1−x)2=∑n=0∞(n+1)xn,\frac{1}{(1-x)^2}=\sum_{n=0}^{\infty}(n+1)x^n,

while integration gives

−ln⁡(1−x)=∑n=1∞xnn.-\ln(1-x)=\sum_{n=1}^{\infty}\frac{x^n}{n}.

Similarly,

11+x=∑n=0∞(−1)nxn\frac{1}{1+x}=\sum_{n=0}^{\infty}(-1)^nx^n

leads, after integration, to

ln⁡(1+x)=∑n=1∞(−1)n+1xnn,−1<x≤1.\ln(1+x)=\sum_{n=1}^{\infty}(-1)^{n+1}\frac{x^n}{n}, \qquad -1<x\leq 1.

Takeaway: Differentiation creates factors involving the index, while integration increases powers and divides coefficients by the new exponent.

Taylor Expansions and Standard

The centered at aa is

∑n=0∞f(n)(a)n!(x−a)n.\sum_{n=0}^{\infty}\frac{f^{(n)}(a)}{n!}(x-a)^n.

Its degree-
nn is

Tn(x)=∑k=0nf(k)(a)k!(x−a)k.T_n(x)=\sum_{k=0}^{n}\frac{f^{(k)}(a)}{k!}(x-a)^k.

When a=0a=0, the expansion is a . Important examples include

ex=∑n=0∞xnn!,e^x=\sum_{n=0}^{\infty}\frac{x^n}{n!},
sin⁡x=∑n=0∞(−1)nx2n+1(2n+1)!,\sin x=\sum_{n=0}^{\infty}(-1)^n\frac{x^{2n+1}}{(2n+1)!},

and

cos⁡x=∑n=0∞(−1)nx2n(2n)!.\cos x=\sum_{n=0}^{\infty}(-1)^n\frac{x^{2n}}{(2n)!}.

The geometric and logarithmic expansions are also fundamental:

11−x=∑n=0∞xn,∣x∣<1,\frac{1}{1-x}=\sum_{n=0}^{\infty}x^n, \qquad |x|<1,

and

ln⁡(1+x)=∑n=1∞(−1)n+1xnn,−1<x≤1.\ln(1+x)=\sum_{n=1}^{\infty}(-1)^{n+1}\frac{x^n}{n}, \qquad -1<x\leq 1.

Having derivatives of every order does not automatically prove that the function equals its . One must establish that the remainder approaches zero on the interval under consideration.

Takeaway: A gives a finite approximation; the infinite represents the function only where its remainder vanishes.

Building New Series from Known Ones

Known series can be adapted through substitution, algebraic rearrangement, multiplication by powers, differentiation, and integration.

For a direct geometric-series construction,

12+x=12⋅11+x/2=12∑n=0∞(−x2)n=∑n=0∞(−1)nxn2n+1,∣x∣<2.\frac{1}{2+x}=\frac12\cdot\frac{1}{1+x/2} =\frac12\sum_{n=0}^{\infty}\left(-\frac{x}{2}\right)^n =\sum_{n=0}^{\infty}\frac{(-1)^nx^n}{2^{n+1}}, \qquad |x|<2.

Substitution into a known series gives

ex2=∑n=0∞x2nn!e^{x^2}=\sum_{n=0}^{\infty}\frac{x^{2n}}{n!}

and

e−x2=∑n=0∞(−1)nx2nn!.e^{-x^2}=\sum_{n=0}^{\infty}(-1)^n\frac{x^{2n}}{n!}.

A nonzero center is useful when the approximation is intended near a point other than zero. For f(x)=1/xf(x)=1/x centered at a=1a=1, write

1x=11+(x−1)=∑n=0∞(−1)n(x−1)n,∣x−1∣<1.\frac1x=\frac{1}{1+(x-1)} =\sum_{n=0}^{\infty}(-1)^n(x-1)^n, \qquad |x-1|<1.

Thus the is (0,2)(0,2).

Takeaway: Rewrite the target function so that it matches a known template, then carry along the corresponding convergence condition.

Approximation and Error Control

If Tn(x)T_n(x) approximates f(x)f(x), the remainder is

Rn(x)=f(x)−Tn(x).R_n(x)=f(x)-T_n(x).

Taylor’s theorem states that, for some cc between aa and xx,

Rn(x)=f(n+1)(c)(n+1)!(x−a)n+1.R_n(x)=\frac{f^{(n+1)}(c)}{(n+1)!}(x-a)^{n+1}.

If ∣f(n+1)(t)∣≤M|f^{(n+1)}(t)|\leq M on the interval, then

∣Rn(x)∣≤M(n+1)!∣x−a∣n+1.|R_n(x)|\leq\frac{M}{(n+1)!}|x-a|^{n+1}.

For an alternating expansion with decreasing term magnitudes, the first omitted term provides a simpler estimate. For example, using

T5(x)=x−x33!+x55!T_5(x)=x-\frac{x^3}{3!}+\frac{x^5}{5!}

to approximate sin⁡(x)\sin(x), the next term has magnitude ∣x∣7/7!|x|^7/7!. At x=π/18x=\pi/18, this gives an error bound of approximately 9.8×10−109.8\times 10^{-10}.

A practical approximation procedure is:

  1. Choose the center and the target value.

  2. Select enough terms to meet the desired accuracy.

  3. Compute the polynomial approximation.

  4. Bound the remainder using a derivative bound or the first omitted alternating term.

  5. Report the approximation together with its error guarantee.

Takeaway: The degree of the approximation should be justified by an explicit error estimate, not chosen only by appearance.

Series for Nonelementary Integrals

Term-by-term integration can evaluate definite integrals whose antiderivatives are not elementary. Since

e−t2=∑n=0∞(−1)nt2nn!,e^{-t^2}=\sum_{n=0}^{\infty}(-1)^n\frac{t^{2n}}{n!},

integrating from 00 to xx gives

∫0xe−t2 dt=∑n=0∞(−1)nx2n+1(2n+1)n!=x−x33+x510−x742+⋯ .\int_0^x e^{-t^2}\,dt =\sum_{n=0}^{\infty}(-1)^n\frac{x^{2n+1}}{(2n+1)n!} =x-\frac{x^3}{3}+\frac{x^5}{10}-\frac{x^7}{42}+\cdots.

At x=1x=1, the first four terms give

1−13+110−142.1-\frac13+\frac1{10}-\frac1{42}.

Because the resulting series is alternating with decreasing term magnitudes, the error is at most the next term,

19⋅4!=1216.\frac{1}{9\cdot 4!}=\frac{1}{216}.

The general strategy is to expand the integrand, integrate each term over the desired limits, and then use a convergence or error test to control the approximation.

Takeaway: Series integration converts a difficult integral into a sequence of elementary power integrals while preserving a way to estimate the numerical error.

Power-Series Solutions of Differential Equations

To solve a differential equation by a , assume

y(x)=∑n=0∞anxn.y(x)=\sum_{n=0}^{\infty}a_nx^n.

Then differentiate term by term:

y′(x)=∑n=1∞nanxn−1,y′′(x)=∑n=2∞n(n−1)anxn−2.y'(x)=\sum_{n=1}^{\infty}na_nx^{n-1}, \qquad y''(x)=\sum_{n=2}^{\infty}n(n-1)a_nx^{n-2}.

Substitute these expressions into the differential equation, rewrite all terms with matching powers of xx, and set each coefficient equal to zero.

For y′′+y=0y''+y=0, alignment of powers gives

(n+2)(n+1)an+2+an=0,(n+2)(n+1)a_{n+2}+a_n=0,

so the coefficient recurrence is

an+2=−an(n+2)(n+1).a_{n+2}=-\frac{a_n}{(n+2)(n+1)}.

The even and odd coefficients develop independently:

a2=−a02!,a4=a04!,a6=−a06!,a_2=-\frac{a_0}{2!},\qquad a_4=\frac{a_0}{4!},\qquad a_6=-\frac{a_0}{6!},

and

a3=−a13!,a5=a15!,a7=−a17!.a_3=-\frac{a_1}{3!},\qquad a_5=\frac{a_1}{5!},\qquad a_7=-\frac{a_1}{7!}.

Therefore,

y=a0cos⁡x+a1sin⁡x.y=a_0\cos x+a_1\sin x.

Initial conditions determine the starting coefficients: a0=y(0)a_0=y(0) and a1=y′(0)a_1=y'(0). For y(0)=2y(0)=2 and y′(0)=−3y'(0)=-3, the solution is

y=2cos⁡x−3sin⁡x.y=2\cos x-3\sin x.

Takeaway: Power-series methods turn a differential equation into algebraic relations among coefficients; the recurrence and initial conditions determine the solution.

A Reliable Problem-Solving Workflow

Use this workflow when working with a power or Taylor expansion:

  1. Identify the center and write all powers as (x−a)n(x-a)^n.

  2. Find the with the ratio test or root test.

  3. Test both endpoints separately whenever the radius is finite.

  4. Choose a construction method: geometric series, substitution, differentiation, integration, or Taylor’s formula.

  5. For approximation, choose the smallest polynomial degree that meets the required accuracy.

  6. Estimate the remainder with Taylor’s theorem or an alternating-series bound.

  7. For differential equations, substitute the assumed series, align powers, and derive the coefficient recurrence.

  8. State the interval on which the resulting series represents the function.

The central ideas are connected: convergence determines where a series is valid, term-by-term operations create new representations, Taylor polynomials provide local approximations, error bounds measure their reliability, and coefficient recurrences extend the method to differential equations.