Taylor series

"Series expansion" redirects here. For other notions of the term, see series.

As the degree of the Taylor series rises, it approaches the correct function. This image shows $\sin x$ and Taylor approximations, polynomials of degree 1, 3, 5, 7, 9, 11 and 13.

The exponential function (in blue), and the sum of the first n+1 terms of its Taylor polynomial at 0 (in red).

In mathematics, the Taylor series is a representation of a function as an infinite sum of terms calculated from the values of its derivatives at a single point. It may be regarded as the limit of the Taylor polynomials. Taylor series are named in honour of English mathematician Brook Taylor. If the series uses the derivatives at zero, the series is also called a Maclaurin series, named after Scottish mathematician Colin Maclaurin.

Definition

The Taylor series of a real or complex function f(x) that is infinitely differentiable in a neighborhood of a real or complex number a, is the power series

$f(a)+{\frac {f'(a)}{1!}}(x-a)+{\frac {f''(a)}{2!}}(x-a)^{2}+{\frac {f^{(3)}(a)}{3!}}(x-a)^{3}+\cdots$

which in a more compact form can be written as

$\sum _{n=0}^{\infty }{\frac {f^{(n)}(a)}{n!}}(x-a)^{n}\,,$

where n! is the factorial of n and f ⁽ⁿ⁾(a) denotes the nth derivative of f at the point a; the zeroth derivative of f is defined to be f itself and (x − a)⁰ and 0! are both defined to be 1.

Often f(x) is equal to its Taylor series evaluated at x for all x sufficiently close to a. This is the main reason why Taylor series are important.

In the particular case where $a=0$, the series is also called a Maclaurin series.

Examples

The Maclaurin series for any polynomial is the polynomial itself.

The Maclaurin series for $(1-x)^{-1}$ is the geometric series

$1+x+x^{2}+x^{3}+\cdots$

so the Taylor series for $x^{-1}$ at $a=1$ is

$1-(x-1)+(x-1)^{2}-(x-1)^{3}+\cdots .$

By integrating the above Maclaurin series we find the Maclaurin series for $-\ln(1-x)$, where $\ln$ denotes the natural logarithm:

$x+{\frac {x^{2}}{2}}+{\frac {x^{3}}{3}}+{\frac {x^{4}}{4}}+\cdots$

and the corresponding Taylor series for $\ln(x)$ at $a=1$ is

$(x-1)-{\frac {(x-1)^{2}}{2}}+{\frac {(x-1)^{3}}{3}}-{\frac {(x-1)^{4}}{4}}+\cdots .$

The Maclaurin series for the exponential function $e^{x}$ at $a=0$ is

$1+{\frac {x^{1}}{1!}}+{\frac {x^{2}}{2!}}+{\frac {x^{3}}{3!}}+{\frac {x^{4}}{4!}}+{\frac {x^{5}}{5!}}+\cdots \qquad =\qquad 1+x+{\frac {x^{2}}{2}}+{\frac {x^{3}}{6}}+{\frac {x^{4}}{24}}+{\frac {x^{5}}{120}}+\cdots \ .$

The above expansion holds because the derivative of $e^{x}$ is also $e^{x}$ and $e^{0}$ equals 1. This leaves the terms $(x-0)^{n}$ in the numerator and n! in the denominator for each term in the infinite sum.

Convergence

The Taylor series need not in general be a convergent series, but often it is. The limit of a convergent Taylor series need not in general be equal to the function value f(x), but often it is. If f(x) is equal to its Taylor series in a neighborhood of a, it is said to be analytic in this neighborhood. If f(x) is equal to its Taylor series everywhere it is called entire. The exponential function $e^{x}$ and the trigonometric functions sine and cosine are examples of entire functions. Examples of functions that are not entire include the logarithm, the trigonometric function tangent, and its inverse arctan. For these functions the Taylor series do not even converge if x is far from a.

The sine function (blue) is closely approximated by its Taylor polynomial of degree 7 (pink) for a full period centered at the origin.

A Taylor series can be used to calculate the value of an entire function in every point, if the value of the function, and of all of its derivatives, is known at a single point. Uses of the Taylor series for entire functions include:

1. The partial sums (the Taylor polynomials) of the series can be used as approximations of the entire function. These approximations are good if sufficiently many terms are included.
2. The series representation simplifies many mathematical proofs.

Pictured on the right is an accurate approximation of sin(x) around the point a = 0. The pink curve is a polynomial of degree seven:

$\sin \left(x\right)\approx x-{\frac {x^{3}}{3!}}+{\frac {x^{5}}{5!}}-{\frac {x^{7}}{7!}}.$

The error in this approximation is no more than ${\tfrac {|x|^{9}}{9!}}$. In particular, for $|x|<1$, the error is less than 0.000003.

History

The Pythagorean philosopher Zeno considered the problem of summing an infinite series to achieve a finite result, but rejected it as an impossibility: the result was Zeno's paradox. Later, Aristotle proposed a philosophical resolution of the paradox, but the mathematical content was apparently unresolved until taken up by Democritus and then Archimedes. It was through Archimedes's method of exhaustion that an infinite number of progressive subdivisions could be performed to achieve a finite trigonometric result. Liu Hui independently employed a similar method several centuries later.

In the 14th century, the earliest examples of the use of Taylor series and closely-related methods were given by Madhava of Sangamagrama. Though no record of his work survives, writings of later Indian mathematicians suggest that he found a number of special cases of the Taylor series, including those for the trigonometric functions of sine, cosine, tangent, and arctangent. The Kerala school of astronomy and mathematics further expanded his works with various series expansions and rational approximations until the 16th century.

In the 17th century, James Gregory also worked in this area and published several Maclaurin series. It was not until 1715 however that a general method for constructing these series for all functions for which they exist was finally provided by Brook Taylor, after whom the series are now named.

The Maclaurin series was named after Colin Maclaurin, a professor in Edinburgh, who published the special case of the Taylor result in the 18th century.

Properties

The function e^−1/x² is not analytic at x = 0: the Taylor series is identically 0, although the function is not.

If this series converges for every x in the interval (a − r, a + r) and the sum is equal to f(x), then the function f(x) is said to be analytic in the interval (a − r, a + r). If this is true for any r then the function is said to be an entire function. To check whether the series converges towards f(x), one normally uses estimates for the remainder term of Taylor's theorem. A function is analytic if and only if it can be represented as a power series; the coefficients in that power series are then necessarily the ones given in the above Taylor series formula.

The importance of such a power series representation is at least fourfold. First, differentiation and integration of power series can be performed term by term and is hence particularly easy. Second, an analytic function can be uniquely extended to a holomorphic function defined on an open disk in the complex plane, which makes the whole machinery of complex analysis available. Third, the (truncated) series can be used to compute function values approximately (often by recasting the polynomial into the Chebyshev form and evaluating it with the Clenshaw algorithm).

Fourth, algebraic operations can often be done much more readily on the power series representation; for instance the simplest proof of Euler's formula uses the Taylor series expansions for sine, cosine, and exponential functions. This result is of fundamental importance in such fields as harmonic analysis.

Note that there are examples of infinitely differentiable functions f(x) whose Taylor series converge, but are not equal to f(x). For instance, for the function defined pointwise by saying that f(x) = e^−1/x² if x ≠ 0 and f(0) = 0, all the derivatives are zero at x = 0, so the Taylor series of f(x) at 0 is zero everywhere, even though the function is nonzero for every x ≠ 0. This particular pathology does not afflict Taylor series in complex analysis. There, the area of convergence of a Taylor series is always a disk in the complex plane (possibly with radius 0), and where the Taylor series converges, it converges to the function value. Notice that e^−1/z² does not approach 0 as z approaches 0 along the imaginary axis, hence this function is not continuous as a function on the complex plane.

Since every sequence of real or complex numbers can appear as coefficients in the Taylor series of an infinitely differentiable function defined on the real line, the radius of convergence of a Taylor series can be zero (Exercise 12 on page 418 in Walter Rudin, Real and Complex Analysis. McGraw-Hill, New Dehli 1980, ISBN 0-07-099557-5). There are even infinitely differentiable functions defined on the real line whose Taylor series have a radius of convergence 0 everywhere (Exercise 13, same book).

Some functions cannot be written as Taylor series because they have a singularity; in these cases, one can often still achieve a series expansion if one allows also negative powers of the variable x; see Laurent series. For example, f(x) = e^−1/x² can be written as a Laurent series.

The Parker-Sochacki method is a recent advance in finding Taylor series which are solutions to differential equations. This algorithm is an extension of the Picard iteration.

List of Taylor series of some common functions

The cosine function.

An 8th degree approximation of the cosine function in the complex plane.

The two above curves put together.

Several important Maclaurin series expansions follow. All these expansions are valid for complex arguments $x$.

Exponential function:

$\mathrm {e} ^{x}=\sum _{n=0}^{\infty }{\frac {x^{n}}{n!}}=1+x+{\frac {x^{2}}{2!}}+{\frac {x^{3}}{3!}}+\cdots \quad {\mbox{ for all }}x$

Natural logarithm:

$\ln(1-x)=-\sum _{n=1}^{\infty }{\frac {x^{n}}{n}}\quad {\mbox{ for }}|x|<1$

Finite geometric series:

${\frac {1-x^{m+1}}{1-x}}=\sum _{n=0}^{m}x^{n}\quad {\mbox{ for }}x\not =1{\mbox{ and }}m\in \mathbb {N} _{0}$

Infinite geometric series:

${\frac {1}{1-x}}=\sum _{n=0}^{\infty }x^{n}\quad {\mbox{ for }}|x|<1$

Variants of the infinite geometric series:

${\frac {x^{m}}{1-x}}=\sum _{n=m}^{\infty }x^{n}\quad {\mbox{ for }}|x|<1{\mbox{ and }}m\in \mathbb {N} _{0}$

${\frac {x}{(1-x)^{2}}}=\sum _{n=1}^{\infty }nx^{n}\quad {\mbox{ for }}|x|<1$

Square root:

${\sqrt {1+x}}=\sum _{n=0}^{\infty }{\frac {(-1)^{n}(2n)!}{(1-2n)n!^{2}4^{n}}}x^{n}\quad {\mbox{ for }}|x|<1$

Binomial series (includes the square root for α = 1/2 and the infinite geometric series for α = −1):

$(1+x)^{\alpha }=\sum _{n=0}^{\infty }{\alpha \choose n}x^{n}\quad {\mbox{ for all }}|x|<1{\mbox{ and all complex }}\alpha$

with generalized binomial coefficients

${\alpha \choose n}=\prod _{k=1}^{n}{\frac {\alpha -k+1}{k}}={\frac {\alpha (\alpha -1)\cdots (\alpha -n+1)}{n!}}$

Trigonometric functions:

$\sin x=\sum _{n=0}^{\infty }{\frac {(-1)^{n}}{(2n+1)!}}x^{2n+1}\quad =x-{\frac {x^{3}}{3!}}+{\frac {x^{5}}{5!}}-\cdots {\mbox{ for all }}x$

$\cos x=\sum _{n=0}^{\infty }{\frac {(-1)^{n}}{(2n)!}}x^{2n}\quad =1-{\frac {x^{2}}{2!}}+{\frac {x^{4}}{4!}}-\cdots {\mbox{ for all }}x$

$\tan x=\sum _{n=1}^{\infty }{\frac {B_{2n}(-4)^{n}(1-4^{n})}{(2n)!}}x^{2n-1}\quad =x+{\frac {x^{3}}{3}}+{\frac {2x^{5}}{15}}+\cdots {\mbox{ for }}|x|<{\frac {\pi }{2}}$

where the Bs are Bernoulli numbers.

$\sec x=\sum _{n=0}^{\infty }{\frac {(-1)^{n}E_{2n}}{(2n)!}}x^{2n}\quad {\mbox{ for }}|x|<{\frac {\pi }{2}}$

$\arcsin x=\sum _{n=0}^{\infty }{\frac {(2n)!}{4^{n}(n!)^{2}(2n+1)}}x^{2n+1}\quad {\mbox{ for }}|x|<1$

$\arctan x=\sum _{n=0}^{\infty }{\frac {(-1)^{n}}{2n+1}}x^{2n+1}\quad {\mbox{ for }}|x|\leq 1$

Hyperbolic functions:

$\sinh x=\sum _{n=0}^{\infty }{\frac {x^{2n+1}}{(2n+1)!}}\quad {\mbox{ for all }}x$

$\cosh x=\sum _{n=0}^{\infty }{\frac {x^{2n}}{(2n)!}}\quad {\mbox{ for all }}x$

$\tanh x=\sum _{n=1}^{\infty }{\frac {B_{2n}4^{n}(4^{n}-1)}{(2n)!}}x^{2n-1}\quad {\mbox{ for }}|x|<{\frac {\pi }{2}}$

$\mathrm {arcsinh} (x)=\sum _{n=0}^{\infty }{\frac {(-1)^{n}(2n)!}{4^{n}(n!)^{2}(2n+1)}}x^{2n+1}\quad {\mbox{ for }}|x|<1$

$\mathrm {arctanh} (x)=\sum _{n=0}^{\infty }{\frac {x^{2n+1}}{2n+1}}\quad {\mbox{ for }}|x|<1$

Lambert's W function:

$W_{0}(x)=\sum _{n=1}^{\infty }{\frac {(-n)^{n-1}}{n!}}x^{n}\quad {\mbox{ for }}|x|<{\frac {1}{\mathrm {e} }}$

The numbers B_k appearing in the summation expansions of tan(x) and tanh(x) are the Bernoulli numbers. The E_k in the expansion of sec(x) are Euler numbers.

Calculation of Taylor series

Several methods exist for the calculation of Taylor series of a large number of functions. One can attempt to use the Taylor series as-is and generalize the form of the coefficients, or one can use manipulations such as substitution, multiplication or division, addition or subtraction of standard Taylor series to construct the Taylor series of a function, by virtue of Taylor series being power series. In some cases, one can also derive the Taylor series by repeatedly applying integration by parts. Particularly convenient is the use of computer algebra systems to calculate Taylor series.

First example

Compute the 7^th degree Maclaurin polynomial for the function

$f(x)=\ln \cos x,\quad x\in (-\pi /2,\pi /2).$

First, rewrite the function as

$f(x)=\ln(1+(\cos x-1))\,.$

We have for the natural logarithm (by using the big O notation)

$\ln(1+x)=x-{\frac {x^{2}}{2}}+{\frac {x^{3}}{3}}+{\mathcal {O}}(x^{4})$

and for the cosine function

$\cos x-1=-{\frac {x^{2}}{2}}+{\frac {x^{4}}{24}}-{\frac {x^{6}}{720}}+{\mathcal {O}}(x^{8}).$

The latter expansion has no constant term on the right-hand side, which enables us to substitute the second series into the first one and to easily omit terms of higher order than the 7^th degree by using the big O notation:

${\begin{aligned}f(x)&=\ln(1+(\cos x-1))\\&={\bigl (}\cos x-1{\bigr )}-{\frac {1}{2}}{\bigl (}\cos x-1{\bigr )}^{2}+{\frac {1}{3}}{\bigl (}\cos x-1{\bigr )}^{3}+{\mathcal {O}}{\bigl (}(\cos x-1)^{4}{\bigr )}\\&={\biggl (}-{\frac {x^{2}}{2}}+{\frac {x^{4}}{24}}-{\frac {x^{6}}{720}}+{\mathcal {O}}(x^{8}){\biggr )}-{\frac {1}{2}}{\biggl (}-{\frac {x^{2}}{2}}+{\frac {x^{4}}{24}}+{\mathcal {O}}(x^{6}){\biggr )}^{2}+{\frac {1}{3}}{\biggl (}-{\frac {x^{2}}{2}}+{\mathcal {O}}(x^{4}){\biggr )}^{3}+{\mathcal {O}}(x^{8})\\&=-{\frac {x^{2}}{2}}+{\frac {x^{4}}{24}}-{\frac {x^{6}}{720}}-{\frac {x^{4}}{8}}+{\frac {x^{6}}{48}}-{\frac {x^{6}}{24}}+{\mathcal {O}}(x^{8})\\&=-{\frac {x^{2}}{2}}-{\frac {x^{4}}{12}}-{\frac {x^{6}}{45}}+{\mathcal {O}}(x^{8}).\end{aligned}}$

Since the cosine is an even function, the coefficients for all the odd powers x, x^{3, x5, x7, . . . have to be zero.}

Second example

Suppose we want the Taylor series at 0 of the function

$g(x)={\frac {e^{x}}{\cos x}}\,.$

We have for the exponential function

$e^{x}=\sum _{n=0}^{\infty }{x^{n} \over n!}=1+x+{x^{2} \over 2!}+{x^{3} \over 3!}+{x^{4} \over 4!}+\cdots$

and, as in the first example,

$\cos x=1-{x^{2} \over 2!}+{x^{4} \over 4!}-\cdots$

Assume the power series is

${e^{x} \over \cos x}=c_{0}+c_{1}x+c_{2}x^{2}+c_{3}x^{3}+\cdots$

Then multiplication with the denominator and substitution of the series of the cosine yields

${\begin{aligned}e^{x}&=(c_{0}+c_{1}x+c_{2}x^{2}+c_{3}x^{3}+\cdots )\cos x\\&=\left(c_{0}+c_{1}x+c_{2}x^{2}+c_{3}x^{3}+c_{4}x^{4}+\cdots \right)\left(1-{x^{2} \over 2!}+{x^{4} \over 4!}-\cdots \right)\\&=c_{0}-{c_{0} \over 2}x^{2}+{c_{0} \over 4!}x^{4}+c_{1}x-{c_{1} \over 2}x^{3}+{c_{1} \over 4!}x^{5}+c_{2}x^{2}-{c_{2} \over 2}x^{4}+{c_{2} \over 4!}x^{6}+c_{3}x^{3}-{c_{3} \over 2}x^{5}+{c_{3} \over 4!}x^{7}+\cdots \end{aligned}}$

Collecting the terms up to fourth order yields

$=c_{0}+c_{1}x+\left(c_{2}-{c_{0} \over 2}\right)x^{2}+\left(c_{3}-{c_{1} \over 2}\right)x^{3}+\left(c_{4}+{c_{0} \over 4!}-{c_{2} \over 2}\right)x^{4}+\cdots$

Comparing coefficients with the above series of the exponential function yields the desired Taylor series

${\frac {e^{x}}{\cos x}}=1+x+x^{2}+{2x^{3} \over 3}+{x^{4} \over 2}+\cdots$

Taylor series as definitions

Classically, the above functions are defined by some property that holds for them. For example, the exponential function is defined as the function that is equal to its own derivative. However, in computable analysis, functions must be defined by algorithms rather than properties, so the above Taylor expansions are used as primary definitions rather than derived results. This is also likely to be the case in software implementations of the functions.

Using Taylor series, one may define analytical functions of matrices and operators, such as matrix exponential or matrix logarithm.

Taylor series for several variables

The Taylor series may also be generalized to functions of more than one variable with

$T(x_{1},\cdots ,x_{d})=\sum _{n_{1}=0}^{\infty }\cdots \sum _{n_{d}=0}^{\infty }{\frac {\partial ^{n_{1}}}{\partial x_{1}^{n_{1}}}}\cdots {\frac {\partial ^{n_{d}}}{\partial x_{d}^{n_{d}}}}{\frac {f(a_{1},\cdots ,a_{d})}{n_{1}!\cdots n_{d}!}}(x_{1}-a_{1})^{n_{1}}\cdots (x_{d}-a_{d})^{n_{d}}.$

For example, for a function that depends on two variables, x and y, the Taylor series to second order about the point (a, b) is:

$f(x,y)\,$

$\approx f(a,b)+f_{x}(a,b)(x-a)+f_{y}(a,b)(y-b)\,$

$+{\frac {1}{2!}}\left[f_{xx}(a,b)(x-a)^{2}+2f_{xy}(a,b)(x-a)(y-b)+f_{yy}(a,b)(y-b)^{2}\right].$

A second-order Taylor series expansion of a scalar-valued function of more than one variable can be compactly written as

$T(\mathbf {x} )=f(\mathbf {a} )+\nabla f(\mathbf {a} )^{T}(\mathbf {x} -\mathbf {a} )+{\frac {1}{2!}}(\mathbf {x} -\mathbf {a} )^{T}\nabla ^{2}f(\mathbf {a} )(\mathbf {x} -\mathbf {a} )+\cdots$

where $\nabla f(\mathbf {a} )$ is the gradient and $\nabla ^{2}f(\mathbf {a} )$ is the Hessian matrix (not to be confused with the Laplacian, which sometimes has the same notation). Applying the multi-index notation the Taylor series for several variables becomes

$T(\mathbf {x} )=\sum _{|\alpha |\geq 0}^{}{{\frac {\mathrm {D} ^{\alpha }f(\mathbf {a} )}{\alpha !}}(\mathbf {x} -\mathbf {a} )^{\alpha }}$

in full analogy to the single variable case.

See also

• Taylor theorem
• Laurent series
• Holomorphic functions are analytic — a proof that a holomorphic function can be expressed as a Taylor power series
• Newton's divided difference interpolation
• Madhava of Sangamagrama (credited with the first use of "Taylor" series)
• Difference engine

References

• Thomas, George B. Jr.; Finney, Ross L. (1996). Calculus and Analytic Geometry (9th ed.). Addison Wesley. ISBN 0-201-53174-7.{{cite book}}: CS1 maint: multiple names: authors list (link)
• Greenberg, Michael (1998). Advanced Engineering Mathematics (2nd ed.). Prentice Hall. ISBN 0-13-321431-1.

External links

• Weisstein, Eric W. "Taylor Series". MathWorld.
• Madhava of Sangamagramma
• Taylor Series Representation Module by John H. Mathews
• "Discussion of the Parker-Sochacki Method"
• Why so much fuss about Taylor Series Expansion?
• Another Taylor visualisation - where you can choose the point of the approximation and the number of derivatives