COMPLEX ANALYSIS: Maple Worksheets, 2009
(c) John H. Mathews Russell W. Howell
mathews@fullerton.edu howell@westmont.edu

Complimentary software to accompany the textbook:

COMPLEX ANALYSIS: for Mathematics & Engineering, 4th Ed, 2001, ISBN: 0-7637-1425-9
Jones and Bartlett Publishers, Inc., 40 Tall Pine Drive, Sudbury, MA 01776
Tele. (800) 832-0034; FAX: (508) 443-8000, E-mail: mkt@jbpub.com, http://www.jbpub.com/

Check out the new Complex Analysis Projects page.


CHAPTER 7 TAYLOR and LAURENT SERIES

Section 7.2 Taylor series Representations

In Section 4.4 we saw that functions defined by power series have derivatives of all orders (Theorem 4.16). In Section 6.5 we saw that analytic functions also have derivatives of all orders (Corollary 6.2). It seems natural, therefore, that there would be some connection between analytic functions and power series. As you might guess, the connection exists via the Taylor and Maclaurin series of analytic functions.

 

 

Definition 7.2: Taylor series

 

If f(z) i s analytic at z = alpha , the series sum(f^`(k)`*`(`*alpha*`)`*(z-alpha)^k/k!,k = 0 .. i... is called the Taylor series for f(z) centered around alpha . When the center is alpha = 0 , the series is called the Maclaurin series for f(z) .

 

Theorem 7.4 (Taylor's Theorem)

 

Suppose f(z) is analytic in a domain G, and that D[R](alpha) = {`z:`*abs(z-alpha) < R} is contained in G. Then the Taylor series converges to f(z) for all z in D[R](alpha) ; that is

 

f(z) = sum(f^`(k)`*`(`*alpha*`)`*(z-alpha)^k/k!,k =... for all z*epsilon*D[R](alpha) .

Furthermore, this representation is valid in the largest disk with center z = alpha t hat is contained in G, and the convergence is uniform on any closed subdisk D[r](alpha) = {`z:`*abs(z-alpha) <= r} for `0 < r` < R .

 

 

Corollary 7.3 Suppose f(z) is analytic in the domain G that contains the point z = alpha . Let z[0] be a singular point of minimum distance to alpha . If abs(z[0]-alpha) = R , then
(i) the Taylor series converges to f(z) on all of D[R](alpha) , and
(ii) if R < S , the Taylor series does not converge to f(z) on all of D[S](alpha) .

Example 7.3, Page 272. Show that 1/((1-z)^2) = sum((n+1)*z^n,n = 0 .. infinity) , is valid for z in D : abs(z) < 1 .

> f:='f': p:='p': P:='P': s:='s': t:='t': z:='z': Z:='Z':
f := z -> 1/(1-z)^2:
t := taylor(f(Z), Z=0, 10):
s := subs(Z=z,t):
p := z -> subs(Z=z,convert(t, polynom)):
`f(z) ` = f(z);
`f(z) ` = s;
P[9](z) = p(z);

`f(z)  ` = 1/((1-z)^2)

`f(z)  ` = series(1+2*z+3*z^2+4*z^3+5*z^4+6*z^5+7*z...

P[9](z) = 1+2*z+3*z^2+4*z^3+5*z^4+6*z^5+7*z^6+8*z^7...

Or we could use Maple's "unapply" procedure.

> f:='f': p:='p': P:='P': s:='s': t:='t': z:='z': Z:='Z':
f := z -> 1/(1-z)^2:
s := taylor(f(z), z=0, 10):
p:=unapply(convert(taylor(f(z),z=0,10),polynom),z):
`f(z) ` = f(z);
`f(z) ` = s;
P[9](z) = p(z);

`f(z)  ` = 1/((1-z)^2)

`f(z)  ` = series(1+2*z+3*z^2+4*z^3+5*z^4+6*z^5+7*z...

P[9](z) = 1+2*z+3*z^2+4*z^3+5*z^4+6*z^5+7*z^6+8*z^7...

Sum up the terms to verify that we have things right.

> f:='f': n:='n': S:='S': z:='z':
f := z -> 1/(1-z)^2:
S9 := z -> sum((n+1)*z^n, n=0..9):
S := z -> sum((n+1)*z^n, n=0..infinity):
`f(z) ` = f(z);
s[9](z),` = `,Sum((n+1)*z^n, n=0..9) = S9(z);
s[infinity](z),` = `,Sum((n+1)*z^n, n=0..infinity) = S(z);

`f(z) ` = 1/((1-z)^2)

s[9](z), ` = `, Sum((n+1)*z^n,n = 0 .. 9) = 1+2*z+3...

s[infinity](z), ` = `, Sum((n+1)*z^n,n = 0 .. infin...

 

Example 7.4, Page 273. (a) Show that 1/(1-z^2) = sum(z^(2*n),n = 0 .. infinity) , is valid for z in D : abs(z) < 1 .

> f:='f': p:='p': P:='P': s:='s': t:='t': z:='z': Z:='Z':
f := z -> 1/(1-z^2):
t := taylor(f(Z), Z=0, 20):
s := subs(Z=z,t):
p := z -> subs(Z=z,convert(t, polynom)):
`f(z) ` = f(z);
`f(z) ` = s;
P[18](z) = p(z);

`f(z)  ` = 1/(1-z^2)

`f(z)  ` = series(1+1*z^2+1*z^4+1*z^6+1*z^8+1*z^10+...

P[18](z) = 1+z^2+z^4+z^6+z^8+z^10+z^12+z^14+z^16+z^...

Or we could use Maple's "unapply" procedure.

> f:='f': p:='p': P:='P': s:='s': t:='t': z:='z': Z:='Z':
f := z -> 1/(1-z^2):
s := taylor(f(z), z=0, 20):
p:=unapply(convert(taylor(f(z),z=0,20),polynom),z):
`f(z) ` = f(z);
`f(z) ` = s;
P[18](z) = p(z);

`f(z)  ` = 1/(1-z^2)

`f(z)  ` = series(1+1*z^2+1*z^4+1*z^6+1*z^8+1*z^10+...

P[18](z) = 1+z^2+z^4+z^6+z^8+z^10+z^12+z^14+z^16+z^...

Sum up the terms to verify that we have things right.

> f:='f': n:='n': s:='s': S:='S': z:='z':
f := z -> 1/(1-z^2):
S9 := z -> sum(z^(2*n), n=0..9):
S := z -> sum(z^(2*n), n=0..infinity):
`f(z) ` = f(z);
s[18](z),` = `,Sum(z^(2*n), n=0..9) = S9(z);
s[infinity](z),` = `,Sum(z^(2*n), n=0..infinity) = S(z);

`f(z) ` = 1/(1-z^2)

s[18](z), ` = `, Sum(z^(2*n),n = 0 .. 9) = 1+z^2+z^...

s[infinity](z), ` = `, Sum(z^(2*n),n = 0 .. infinit...

 

(b) Show that 1/(1+z^2) = sum((-1)^n*z^(2*n),n = 0 .. infinity) , is valid for z in D : abs(z) < 1 .

> f:='f': p:='p': P:='P': s:='s': t:='t': z:='z': Z:='Z':
f := z -> 1/(1+z^2):
t := taylor(f(Z), Z=0, 20):
s := subs(Z=z,t):
p := z -> subs(Z=z,convert(t, polynom)):
`f(z) ` = f(z);
`f(z) ` = s;
P[18](z) = p(z);

`f(z)  ` = 1/(z^2+1)

`f(z)  ` = series(1-1*z^2+1*z^4-1*z^6+1*z^8-1*z^10+...

P[18](z) = 1-z^2+z^4-z^6+z^8-z^10+z^12-z^14+z^16-z^...

Or we could use Maple's "unapply" procedure.

> f:='f': p:='p': P:='P': s:='s': t:='t': z:='z': Z:='Z':
f := z -> 1/(1+z^2):
s := taylor(f(z), z=0, 20):
p:=unapply(convert(taylor(f(z),z=0,20),polynom),z):
`f(z) ` = f(z);
`f(z) ` = s;
P[18](z) = p(z);

`f(z)  ` = 1/(z^2+1)

`f(z)  ` = series(1-1*z^2+1*z^4-1*z^6+1*z^8-1*z^10+...

P[18](z) = 1-z^2+z^4-z^6+z^8-z^10+z^12-z^14+z^16-z^...

 

Sum up the terms to verify that we have things right.

> f:='f': n:='n': s:='s': S:='S': z:='z':
f := z -> 1/(1+z^2):
S9 := z -> sum((-1)^n * z^(2*n), n=0..9):
S := z -> sum((-1)^n * z^(2*n), n=0..infinity):
`f(z) ` = f(z);
s[18](z),` = `,Sum((-1)^n * z^(2*n), n=0..9) = S9(z);
s[infinity](z),` = `,Sum((-1)^n * z^(2*n), n=0..infinity) = S(z);

`f(z) ` = 1/(z^2+1)

s[18](z), ` = `, Sum((-1)^n*z^(2*n),n = 0 .. 9) = 1...

s[infinity](z), ` = `, Sum((-1)^n*z^(2*n),n = 0 .. ...

 

 

Theorem 7.5 (Uniqueness of power series) Suppose that in some disk D[r](alpha) we have

 

f(z) = sum(a[n]*(z-alpha)^n,n = 0 .. infinity) = sum(b[n]*... . Then a[n] = b[n] for n = 0, 1, 2, `...` .


Example 7.5, Page 274. Find the Maclaurin series of f(z) = sin(z)^3 .

> f:='f': p:='p': P:='P': s:='s': t:='t': z:='z': Z:='Z':
f := z -> sin(z)^3:
t := taylor(f(Z), Z=0, 14):
s := subs(Z=z,t):
p := z -> subs(Z=z,convert(t, polynom)):
`f(z) ` = f(z);
`f(z) ` = s;
P[13](z) = p(z);

`f(z)  ` = sin(z)^3

`f(z)  ` = series(1*z^3-1/2*z^5+13/120*z^7-41/3024*...

P[13](z) = z^3-1/2*z^5+13/120*z^7-41/3024*z^9+671/6...

 

Sum up the terms to verify that we have things right.

> An:='An': Bn:='Bn': f:='f': n:='n': p:='p': q:='q': s:='s': z:='z':
f := z -> sin(z)^3:
An := (-1)^n * 3/(4*(2*n+1)!):
Bn := - (-1)^n * 3 * 9^n/(4*(2*n+1)!):
`f(z) ` = f(z);
p13 := sum(An* z^(2*n+1), n=0..6):
q13 := sum(Bn* z^(2*n+1), n=0..6):
Sum(An* z^(2*n+1), n=0..6) = p13;
Sum(Bn* z^(2*n+1), n=0..6) = q13;
`Partial sum `, s[13](z) =
sum(An* z^(2*n+1), n=0..6) +
sum(Bn* z^(2*n+1), n=0..6);
s[infinity](z),` = `,
Sum(An* z^(2*n+1)+Bn* z^(2*n+1), n=0..infinity) =
sum(An* z^(2*n+1), n=0..infinity) +
sum(Bn* z^(2*n+1), n=0..infinity);

`f(z) ` = sin(z)^3

Sum(3/4*(-1)^n/(2*n+1)!*z^(2*n+1),n = 0 .. 6) = 3/4...

Sum(-3/4*(-1)^n*9^n/(2*n+1)!*z^(2*n+1),n = 0 .. 6) ...

`Partial sum  `, s[13](z) = z^3-1/2*z^5+13/120*z^7-...

s[infinity](z), ` = `, Sum(3/4*(-1)^n/(2*n+1)!*z^(2...

 

 

Theorem 7.6 Let f(z) and g(z) have the power series representations
f(z) = sum(a[n]*(z-alpha)^n,n = 0 .. infinity) for z*epsilon*D[r[1]](alpha) and g(z) = sum(b[n]*(z-alpha)^n,n = 0 .. infinity) for z*epsilon*D[r[2]](alpha) .
If
r = min(r[1],r[2]) , and beta is any complex constant, then
(i) beta*f(z) = sum(beta*a[n]*(z-alpha)^n,n = 0 .. infi... for z*epsilon*D[r[1]](alpha) ,
(ii) f(z)+g(z) = sum((a[n]+b[n])*(z-alpha)^n,n = 0 .. in... for z*epsilon*D[r](alpha) , and
(iii) f(z)*g(z) = sum(c[n]*(z-alpha)^n,n = 0 .. infinity)... for z*epsilon*D[r](alpha) , where c[n] = sum(a[k]*b[n-k],k = 0 .. n) .

Identity (iii) is known as the Cauchy product of the series for f(z) and g(z) .

 

End of Section 7.2.