Calculus/Taylor series

From Wikibooks, open books for an open world
Jump to: navigation, search

Taylor Series[edit]

As the degree of the Taylor series rises, it approaches the correct function. sin(x) and Taylor approximations, polynomials of degree 1, 3, 5, 7, 9, 11 and 13.

The Taylor series of an infinitely, often differentiable real, (or complex) function f defined on an open interval (a-r, a+r) is the power series 
\sum_{n=0}^{\infin} \frac{f^{(n)}(a)}{n!} (x-a)^{n}

Here, n! is the factorial of n and f (n)(a) denotes the nth derivative of f at the point a. 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 called analytic. 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 a power series converges to the function; the coefficients in that power series are then necessarily the ones given in the above Taylor series formula.

If a = 0, the series is also called a Maclaurin series.

The importance of such a power series representation is threefold. 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 approximate values of the function near the point of expansion.

Around zero, the function looks very flat. The function e-1/x² is not analytic: the Taylor series is 0, although the function is not.

Note that there are examples of infinitely often differentiable functions f(x) whose Taylor series converge, but are not equal to f(x). For instance, for the function defined piecewise by saying that f(x) = exp(−1/x²) if x ≠ 0 and f(0) = 0, all the derivatives are zero at x = 0, so the Taylor series of f(x) is zero, and its radius of convergence is infinite, even though the function most definitely is not zero. This particular pathology does not afflict complex-valued functions of a complex variable. Notice that exp(−1/z²) does not approach 0 as z approaches 0 along the imaginary axis.

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) = exp(−1/x²) can be written as a Laurent series.

The Parker-Sockacki theorem is a recent advance in finding Taylor series which are solutions to differential equations. This theorem is an expansion on the Picard iteration.

Derivation/why this works[edit]

If a function f(x) is written as a infinite power series, it will look like this:


where a is half the radius of convergence and c0,c1,c2,c3,c4... are coefficients. If we substitute a for x:


If we differentiate:


If we substitute a for x:


If we differentiate:


If we substitute a for x:




where f0(x)=f(x) and f1(x)=f´(x) and so on.We can actually go ahead and say that the power approximation of f(x) is:


<this needs to be improved>

List of Taylor series[edit]

Several important Taylor series expansions follow. All these expansions are also valid for complex arguments x.

Exponential function and natural logarithm:

e^{x} = \sum^{\infin}_{n=0} \frac{x^n}{n!}\quad\mbox{ for all } x
\ln(1+x) = \sum^{\infin}_{n=1} \frac{(-1)^{n+1}}n x^n\quad\mbox{ for } \left| x \right| < 1

Geometric series:

\frac{1}{1-x} = \sum^{\infin}_{n=0} x^n\quad\mbox{ for } \left| x \right| < 1

Binomial series:

(1+x)^\alpha = \sum^{\infin}_{n=0} C(\alpha,n) x^n\quad\mbox{ for all } \left| x \right| < 1\quad\mbox{ and all complex } \alpha

Trigonometric functions:

\sin x = \sum^{\infin}_{n=0} \frac{(-1)^n}{(2n+1)!} x^{2n+1}\quad\mbox{ for all } x
\cos x = \sum^{\infin}_{n=0} \frac{(-1)^n}{(2n)!} x^{2n}\quad\mbox{ for all } x
\tan x = \sum^{\infin}_{n=1} \frac{B_{2n} (-4)^n (1-4^n)}{(2n)!} x^{2n-1}\quad\mbox{ for } \left| x \right| < \frac{\pi}{2}
\sec x = \sum^{\infin}_{n=0} \frac{(-1)^n E_{2n}}{(2n)!} x^{2n}\quad\mbox{ for } \left| x \right| < \frac{\pi}{2}
\arcsin x = \sum^{\infin}_{n=0} \frac{(2n)!}{4^n (n!)^2 (2n+1)} x^{2n+1}\quad\mbox{ for } \left| x \right| < 1
\arctan x = \sum^{\infin}_{n=0} \frac{(-1)^n}{2n+1} x^{2n+1}\quad\mbox{ for } \left| x \right| < 1

Hyperbolic functions:

\sinh x = \sum^{\infin}_{n=0} \frac{1}{(2n+1)!} x^{2n+1}\quad\mbox{ for all } x
\cosh x = \sum^{\infin}_{n=0} \frac{1}{(2n)!} x^{2n}\quad\mbox{ for all } x
\tanh x = \sum^{\infin}_{n=1} \frac{B_{2n} 4^n (4^n-1)}{(2n)!} x^{2n-1}\quad\mbox{ for } \left| x \right| < \frac{\pi}{2}
\sinh^{-1} x = \sum^{\infin}_{n=0} \frac{(-1)^n (2n)!}{4^n (n!)^2 (2n+1)} x^{2n+1}\quad\mbox{ for } \left| x \right| < 1
\tanh^{-1} x = \sum^{\infin}_{n=0} \frac{1}{2n+1} x^{2n+1}\quad\mbox{ for } \left| x \right| < 1

Lambert's W function:

W_0(x) = \sum^{\infin}_{n=1} \frac{(-n)^{n-1}}{n!} x^n\quad\mbox{ for } \left| x \right| < \frac{1}{e}

The numbers Bk appearing in the expansions of tan(x) and tanh(x) are the Bernoulli numbers. The C(α,n) in the binomial expansion are the binomial coefficients. The Ek in the expansion of sec(x) are Euler numbers.

Multiple dimensions[edit]

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

\sum_{n_1=0}^{\infin} \cdots \sum_{n_d=0}^{\infin}
\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}


The Taylor series is named for mathematician Brook Taylor, who first published the power series formula in 1715.

Constructing a Taylor Series[edit]

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 (such as those above) 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. The use of computer algebra systems to calculate Taylor series is common, since it eliminates tedious substitution and manipulation.

Example 1[edit]

Consider the function

f(x)=\ln{(1+\cos{x})} \,,

for which we want a Taylor series at 0.

We have for the natural logarithm

\ln(1+x) = \sum^{\infin}_{n=1} \frac{(-1)^{n+1}}{n} x^n = x - {x^2\over 2}+{x^3 \over 3} - {x^4 \over 4} + \cdots \quad\mbox{ for } \left| x \right| < 1

and for the cosine function

\cos x = \sum^{\infin}_{n=0} \frac{(-1)^n}{(2n)!} x^{2n} = 1 -{x^2\over 2!}+{x^4\over 4!}- \cdots \quad\mbox{ for all } x\in\mathbb{C}.

We can simply substitute the second series into the first. Doing so gives

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

Expanding by using multinomial coefficients gives the required Taylor series. Note that cosine and therefore f are even functions, meaning that f(x)=f(-x), hence the coefficients of the odd powers x, x^3, x^5, x^7 and so on have to be zero and don't need to be calculated. The first few terms of the series are

\ln(1+\cos x)=\ln 2-{x^2\over 4}-{x^4\over 96}-{x^6\over 1440} -{17x^8\over 322560}-{31x^{10}\over 7257600}-\cdots

The general coefficient can be represented using Faà di Bruno's formula. However, this representation does not seem to be particularly illuminating and is therefore omitted here.

Example 2[edit]

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^\infty_{n=0} {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{align} 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_4x^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_1x - {c_1 \over 2}x^3 + {c_1 \over 4!}x^5 + c_2x^2 - {c_2 \over 2}x^4 + {c_2 \over 4!}x^6 + c_3x^3 - {c_3 \over 2}x^5 + {c_3 \over 4!}x^7 +\cdots \end{align}

Collecting the terms up to fourth order yields

=c_0 + c_1x + \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


Generalized Mean Value Theorem[edit]