next up previous
Next: About this document ... Up: lab_template Previous: lab_template

Subsections


Operations on power series

Purpose

The purpose of this lab is to acquaint you with several useful techniques for working with power series. These techniques include substitution and multiplication and division of series by simple polynomials.

Getting Started

To assist you, there is a worksheet associated with this lab that contains examples and even solutions to some of the exercises. You can copy that worksheet to your home directory with the following command, which must be run in a terminal window, not in Maple.

cp ~bfarr/Powseries_start.mws ~

You can copy the worksheet now, but you should read through the lab before you load it into Maple. Once you have read to the exercises, start up Maple, load the worksheet Powseries_start.mws, and go through it carefully. Then you can start working on the exercises.

Background

The general form of a power series in $x-a$ is given below.

\begin{displaymath}f(x) = \sum_{k=0}^{\infty} a_k (x-a)^k \end{displaymath}

The number $a$ is called the base point of the power series. In this lab, we will consider only the special case $a=0$. Historically, power series have been used most often to approximate functions that do not have simple formulas. One example is the exponential function, whose power series is given below along with some other important power series you have seen.

The most familiar example of a power series is the geometric series.

\begin{displaymath}\frac{1}{1-x} = \sum_{k=0}^{\infty} x^k \mbox{, for $\mid x
\mid < 1$} \end{displaymath}

Once you have a convergent series representation for a function, it can be manipulated in several ways to generate convergent power series representations of related functions. The rest of the Background describes several different techniques for generating power series representations of functions that are related to power series that are already known. The four techniques are substitution, multiplication and division, integration, and differentiation.

Substitution

Suppose you have a power series representation for $f(x)$, and you want a power series representation for $f(ax^p)$ where $a$ is a constant and $p$ is a positive integer The theory for power series says that you can obtain the pwer series you want simply by subsitution. For example, to get a power series for $\exp(2x)$ you just take the series for $\exp(x)$ and substitute $2x$ for $x$, obtaining

\begin{displaymath}\exp(2x) = \sum_{k=0}^{\infty} \frac{(2x)^k}{k!}. \end{displaymath}

Writing out the first few terms gives

\begin{displaymath}\exp(2x) = 1 + 2x +2x^2 + \frac{4 x^3}{3} + \ldots \end{displaymath}

If you use such a substitution, you have to be careful if the series is only valid for a finite interval. For example, suppose you wanted to find the power series for

\begin{displaymath}f(t) = \frac{1}{1+2t}\end{displaymath}

You can obtain the desired series by substitution as

\begin{displaymath}\frac{1}{1+2t} = \sum_{k=0}^{\infty} (-1)^k (2t)^k \end{displaymath}

but you have to be careful because this formula is not valid for all values of $t$. In fact this formula is only valid if $\mid t \mid < 1/2$. The reason for this is that the series for $1/(1-x)$ is only valid if $\mid x \mid < 1$ and when we substitute $-2t$ for $x$, the formula only makes sense if $\mid
2t \mid < 1$.

Multiplication and Division

If you have a power series representation for $f(x)$, and you want the power series for something like $x^2 f(x)$, you can just multiply each term of the series for $f(x)$ by $x^2$. If the leading term for the power series representation of $g(x)$ is $x^k$ for some integer $k > 0$, you can use division to obtain the power series representation for $g(x)/x^n$ for any integer $n \leq k$. Some examples follow.

\begin{displaymath}\frac{x}{1-x} = x +x^2+x^3 + \ldots \end{displaymath}


\begin{displaymath}x^2 \exp(x) = x^2 + x^3 + \frac{x^4}{2} + \ldots \end{displaymath}


\begin{displaymath}\frac{\exp(x)-1}{x} = 1+\frac{x}{2!}+\frac{x^{2}}{3!} +
\ldots \end{displaymath}

Term-by-term integration and differentiation

One of the nice properties of power series is that they can be integrated and differentiated term-by-term. For example,

\begin{displaymath}\frac{d}{dx} \frac{1}{1-x} = \frac{1}{(1-x)^2} \end{displaymath}

so we have the following power series representation.

\begin{displaymath}\frac{1}{(1-x)^2} = 1+2x+3x^2+4x^3+ \ldots \mbox{, for $\mid x
\mid < 1$} \end{displaymath}

Note that the interval of convergence is exactly the same as for the original series, neglecting the behavior at the endpoints.

Integrating a power series term-by-term is very similar, but you may have to include a constant of integration. For example, integrating the power series representation for $\exp(x)$ term by term gives

\begin{displaymath}\int 1+x + \frac{x^2}{2!} + \frac{x^3}{3!} + \ldots \, dx = C +x +
\frac{x^2}{2!} + \frac{x^3}{3!} + \frac{x^{4}}{4!} + \ldots \end{displaymath}

You would have to set $C=1$ to make the right hand side the power series representation for $\exp(x)$.

Exercises

  1. Verify that substitution and multiplication work as described above to generate partial sums of power series (with base point a=0) for the following functions. That is, compare the partial sums of power series up to various orders obtained directly with those obtained by substitution, multiplication, or division.
    1. $f(x) = \exp(2x)$.
    2. $f(x) = x \exp(-x^2)$.
    3. $f(x) = (1-\cos(x))/x^2$
    4. $f(x) = \cos(x^2)$
    5. $f(x) = x/(1-x^2)$
  2. Use substitution followed by integration to find the first four non-zero terms in the power series representation for $\ln(1+x)$. Start with the series for $1/(1-x)$.
  3. In the book it says that power series can also be added term by term. Try the following two Maple commands and verify that the output is the same.
    > 1/2*Taylor(exp(x),x=0,5)+1/2*Taylor(exp(-x),x=0,5);
    > Taylor((exp(x)+exp(-x))/2,x=0,5);
    
    Can you explain why the power series only contains even powers of $x$?

next up previous
Next: About this document ... Up: lab_template Previous: lab_template
William W. Farr
2003-01-28