Main /
## ApproximationAndError< Taylor series redux | Home Page | Calculus > Given a series that is known to converge but for which an exact answer is not known, how does one find a good approximation to the true value? One way to get an approximation is to add up some number of terms and then stop. But how many terms are enough? How close will the result be to the true answer? That is the motivation for this module. ## Error definedGiven a convergent series (:latex:) \begin{equation*} s = \sum_{n=0}^\infty a_n. \end{equation*} (:latexend:) Recall that the partial sum ($s_k$) is the sum of the terms up to and including ($a_k$), i.e., (:latex:)
\begin{align*}
s_k &= a_0+a_1+a_2+\ldots+a_k Then the (:latex:)
\begin{align*}
E_k &= s - s_k In other words, the error is the sum of all the terms from the infinite series which were not included in the partial sum. ## Alternating series error boundFor a decreasing, alternating series, it is easy to get a bound on the error ($E_k$): (:latex:) \begin{equation*} |E_k| \leq |a_{k+1}|. \end{equation*} (:latexend:) In other words, the error is bounded by the next term in the series. ## NoteIf the series is strictly decreasing (as is usually the case), then the above inequality is strict. To see why the alternating bound holds, note that each successive term in the series overshoots the true value of the series. In other words, if ($S$) is the true value of the series, (:latex:)
\begin{align*}
a_0 &> S The above figure shows that if one stops at ($a_0 - a_1+a_2-a_3$), then the error ($E_3$) must be less than ($a_4$). ## ExampleWhat is the minimum number of terms of the series (:latex:) \begin{equation*} \sum_{n=1}^\infty (-1)^{n+1} \frac{1}{n^2} \end{equation*} (:latexend:) one needs to be sure to be within ($\frac{1}{100}$) of the true sum? The goal is to find ($k$) so that ($|E_k| \leq \frac{1}{100}$). Since ($|E_k| \leq |a_{k+1}|$), the question becomes for which value of ($k$) is ($|a_{k+1}| \leq \frac{1}{100}$)? If ($k=9$), then ($|a_{k+1}| = \frac{1}{100}$), and so by the alternating series error bound, ($|E_9| \leq \frac{1}{100}$). Thus 9 terms are required to be within ($\frac{1}{100}$) of the true value of the series. ## Integral test for error boundsAnother useful method to estimate the error of approximating a series is a corollary of the integral test. Recall that if a series ($\sum f(n)$) has terms which are positive and decreasing, then (:latex:) \begin{equation*} \int_{m+1}^\infty f(x)dx < \sum_{n=m+1}^\infty f(n) < \int_{m}^\infty f(x)dx. \end{equation*} (:latexend:) But notice that the middle quantity is precisely ($E_m$). So (:latex:) \begin{equation*} \int_{m+1}^\infty f(x)dx < E_m < \int_{m}^\infty f(x)dx. \end{equation*} (:latexend:) This bound is nice because it gives an upper bound ## ExampleHow many terms of the series (:latex:) \begin{equation*} \sum_{n=1}^\infty \frac{1}{n^3} \end{equation*} (:latexend:) must one add up so that the Integral bound guarantees the approximation is within ($\frac{1}{100}$) of the true answer? (:toggle hide show="Answer" box2:) If one adds up the first ($m$) terms, then by the integral bound, the error ($E_m$) satisfies (:latex:)
\begin{align*}
E_m &< \int_m^\infty \frac{dx}{x^3} Setting ($\frac{1}{2m^2} \leq \frac{1}{100}$) gives that ($m^2 \geq 50$), so ($m \geq 8$). Thus, ($m=8$) is the minimum number of terms required so that the Integral bound guarantees we are within ($\frac{1}{100}$) of the true answer. ## NoteIf you actually compute the partial sums using a calculator, you will find that 7 terms actually suffice. But remember, we want the ## Taylor approximationsRecall that the (:latex:)
\begin{align*}
f(x) &= \sum_{n=0}^\infty \frac{f^{\left(n\right)}(0)}{n!} x^n The (:latex:)
\begin{align*}
f(x) &\approx \sum_{n=0}^N \frac{f^{\left(n \right)}(0)}{n!} x^n This is a good approximation for ($f(x)$) when ($x$) is close to 0. How good an approximation is it? That is the purpose of the last error estimate for this module. As in previous modules, let ($E_N(x)$) be the error between the Taylor polynomial and the true value of the function, i.e., (:latex:) \[ f(x) = \sum_{n=0}^N \frac{f^{\left(n \right)}(0)}{n!} x^n + E_N(x). \] (:latexend:) Notice that the error ($E_N(x)$) is a function of ($x$). In general, the further away ($x$) is from ($0$), the bigger the error will be. A first, weak bound for the error is given by (:latex:) \[ E_N(x) \leq C x^{N+1} \] (:latexend:) for some constant ($C$) and ($x$) sufficiently close to 0. In other words, ($E_N(x)$) is ($O(x^{N+1})$). A stronger bound is given in the next section. ## Taylor remainder theoremThe following gives the precise error from truncating a Taylor series: Taylor remainder theorem The error ($E_N(x)$) is given precisely by (:latex:) \begin{equation*} E_N(x) = \frac{f^{\left(N+1\right)}(t)}{(N+1)!}x^{N+1}, \end{equation*} (:latexend:) for some ($t$) between 0 and ($x$), inclusive. So if ($x<0$), then ($x \leq t \leq 0$), and if ($x>0$), then ($0 \leq t \leq x$). ## ExampleConsider the case when ($N=0$). The Taylor remainder theorem says that (:latex:) \begin{equation*} f(x) = f(0) + f'(t)x, \end{equation*} (:latexend:) for some ($t$) between 0 and ($x$). Solving for ($f'(t)$) gives (:latex:) \begin{equation*} f'(t) = \frac{f(x)-f(0)}{x-0}, \end{equation*} (:latexend:) for some ($0<t<x$) if ($x>0$) and ($x<t<0$) if ($x<0$), which is precisely the statement of the Mean value theorem. Therefore, one can think of the Taylor remainder theorem as a generalization of the Mean value theorem. ## Taylor error boundAs it is stated above, the Taylor remainder theorem is not particularly useful for actually finding the error, because there is no way to actually find the ($t$) for which the equation holds. There is a slightly different form which gives a bound on the error: Taylor error bound (:latex:) \begin{equation*} |E_N(x)| \leq \frac{C}{(N+1)!} |x|^{N+1}, \end{equation*} (:latexend:) where ($C$) is the maximum value of ($|f^{\left(N+1\right)}(t)|$) over all ($t$) between 0 and ($x$), inclusive. ## ExampleEstimate ($\sqrt{e}$) using (:latex:) \begin{equation*} e^{1/2} \approx 1 + \frac{1}{2} + \frac{(1/2)^2}{2!} + \frac{(1/2)^3}{3!} \approx 1.64, \end{equation*} (:latexend:) and bound the error. (:toggle hide show="Answer" box1:) The function is ($f(x) = e^x$), and the approximating polynomial used here is (:latex:) \begin{equation*} e^x \approx 1 + x + \frac{x^2}{2!} + \frac{x^3}{3!}. \end{equation*} (:latexend:) Then according to the above bound, (:latex:) \begin{equation*} |E_3(x)| \leq \frac{C}{4!}|x|^4, \end{equation*} (:latexend:) where ($C$) is the maximum of ($f^{\left(4 \right)}(t) = e^t$) for ($0 \leq t \leq x$). Since ($e^t$) is an increasing function, ($C = e^x$). Thus, (:latex:) \begin{equation*} |E_3(x)| \leq \frac{e^x}{4!}x^4. \end{equation*} (:latexend:) Thus, (:latex:) \begin{equation*} |E_3(1/2)| \leq \frac{e^{1/2}}{4!}(1/2)^4 < \frac{1}{100}. \end{equation*} (:latexend:) < Taylor series redux | Home Page | Calculus > |

Page last modified on August 22, 2013, at 01:00 PM