COMPLEX ANALYSIS: Maple Worksheets, 2001
(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 6 COMPLEX INTEGRATION

Section 6.5 Integral Representations for Analytic Functions

We now present some major results in the theory of functions of a complex variable. The first result is known as Cauchy's integral formula and shows that the value of an analytic function
f can be represented by a certain contour integral. The `n th` derivative , f^`(n)` `(z)` , will have a similar representation. In Chapter 7 we will show how the Cauchy integral formulae are used to prove Taylor's theorem, and we will establish the power series representation for analytic functions. The Cauchy integral formulae will also be a convenient tool for evaluating certain contour integrals.

 

Theorem 6.10 (Cauchy's Integral Formula)

Let f(z) be analytic in the simply connected domain D , and let C be a simple closed positively oriented contour that lies in D.
If is a point that lies interior to
C , then int(f(z)/(z-z[0]),z = C .. `   `) = 2*pi*i*f(z[0]) .

Example 6.21, Page 245. Show that int(exp(z)/(z-1),z = C .. `   `) = 2*pi*i*exp(1) ,
where
C is the circle C : abs(z) = 2 with positive orientation.

> f :='f': F:='F': z:='z': Z:='Z': Zn:='Zn':
w := exp(z)/(z-1):
print(`Find `,Int(w, z=C..``));
Zn := sort([solve(denom(w)=0, z)]):
`Singularity at `, z[0] = Zn[1];
f := Z -> subs(z=Z, w*(z - Zn[1])):
`Use f(z) ` = f(z);

`Find `, Int(exp(z)/(z-1),z = C .. ``)

`Singularity at  `, z[0] = 1

`Use  f(z) ` = exp(z)

 

The integral of f(z) taken over C is:

> `f(z)` = f(z);
`f(z0)` = f(Zn[1]);
val := 2*Pi*I*f(Zn[1]):
Int(w ,z=C..``) = `2 Pi I f(z0)`;
Int(w ,z=C..``) = val;

`f(z)` = exp(z)

`f(z0)` = exp(1)

Int(exp(z)/(z-1),z = C .. ``) = `2 Pi I f(z0)`

Int(exp(z)/(z-1),z = C .. ``) = 2*I*Pi*exp(1)

 

Example 6.22, Page 245. Show that int(sin(z)/(4*z+pi),z = C .. `   `) = -i*sqrt(2)*pi... ,
where
C is the circle C : abs(z) = 1 with positive orientation.

> f :='f': F:='F': z:='z': Z:='Z': Zn:='Zn':
w := sin(z)/(4*z + Pi):
print(`Find `,Int(w ,z=C..``));
Zn := sort([solve(denom(w)=0, z)]):
`Singularity at `, z[0] = Zn[1];
`Singularity at `, z[0] = evalf(Zn[1]);
f := Z -> subs(z=Z, simplify(w*(z - Zn[1]))):
`Use f(z) ` = f(z);

`Find `, Int(sin(z)/(4*z+Pi),z = C .. ``)

`Singularity at  `, z[0] = -1/4*Pi

`Singularity at  `, z[0] = -.7853981635

`Use  f(z) ` = 1/4*sin(z)

 

The integral of f(z) taken over C is:

> `f(z)` = f(z);
`f(z0) ` = f(Zn[1]);
`f(z0) ` = evalc(f(Zn[1]));
val := 2*Pi*I*f(Zn[1]):
Int(w ,z=C..``) = `2 Pi I f(z0)`;
Int(w ,z=C..``) = val;

`f(z)` = 1/4*sin(z)

`f(z0) ` = 1/4*sin(-1/4*Pi)

`f(z0) ` = -1/8*sqrt(2)

Int(sin(z)/(4*z+Pi),z = C .. ``) = `2 Pi I f(z0)`

Int(sin(z)/(4*z+Pi),z = C .. ``) = -1/4*I*Pi*sqrt(2...

 

Example 6.23, Page 245. Show that int(exp(i*pi*z)/(2*z^2-5*z+2),z = C .. `   `) = 2*p... ,
where
C is the circle C : abs(z) = 1 with positive orientation.
Solution:
We will find that
z[0] = 1/2 the only singularity that lies inside C .

> f :='f': F:='F': z:='z': Z:='Z': Zn:='Zn':
w := exp(I*Pi*z)/(2*z^2 - 5*z + 2):
print(`Find `,Int(w ,z=C..``));
Zn := sort([solve(denom(w)=0, z)]):
`Singularities occur at: ` = Zn;
`The singularity inside C is `,z[0] = Zn[1];
f := Z -> subs(z=Z, simplify(w*(z - Zn[1]))):
`Use f(z) ` = f(z);

`Find `, Int(exp(I*Pi*z)/(2*z^2-5*z+2),z = C .. ``)...

`Singularities occur at: ` = [1/2, 2]

`The singularity inside C is  `, z[0] = 1/2

`Use  f(z) ` = 1/2*exp(I*Pi*z)/(z-2)

 

 

We now state a general result that shows how differentiation under the integral sign can be accomplished. The proof can be found in some advanced texts. See, for instance, Rolf Nevanlinna and V. Paatero, Introduction to Complex Analysis , (Reading, Massachusetts: Addison-Wesley Publishing Company, 1969), Section 9.7.

 

 

Theorem 6.11 (Leibniz's Rule)

Let D be a simply connected domain, and let `I: ` a `` <= `` t `` <= `` b be an interval of real numbers. Let f(z,t) and its partial derivative diff(f(z,t),z) , with respect to z, be continuous functions for all z in D and all t in I.
Then
F(z) = int(f(z,t),t = a .. b) is analytic for z in D, and `F '(z)` = int(diff(f(z,t),z),t = a .. b) .

 

 

We now show how we can generalize Theorem 6.10 to give an integral representation for the `n th` derivative , f^`(n)` `(z)` . Leibniz's rule is used in the proof, and this method of proof will be a mnemonic device for remembering the upcoming theorem.

 

 

Theorem 6.12 (Cauchy's Integral Formulae for Derivatives)

Let f(z) be analytic in the simply connected domain D , and let C be a simple closed positively oriented contour that lies in D .
If z is a point that lies interior to C, then
f^`(n)`*`(z)` = n!/(2*pi*i) int(f(zeta)/((z-zeta)^(n+1)),zeta = C .. ` `) .

 

Example 6.25, Page 247. Show that int(exp(z^2)/((z-i)^4),z) = -4*pi/(3*exp(1)) .
where
C is the circle C : abs(z) = 2 with positive orientation.
Solution:
Since
z = i is a singularity of order n = 4 , multiply the integrand by (z-i)^4 to get the function f(z) .

> f :='f': F:='F': z:='z': Z:='Z': Zn:='Zn':
w := exp(z^2)/(z - I)^4:
print(`Find `,Int(w ,z=C..``));
Zn := sort([solve(denom(w)=0, z)]):
`Singularity of order 4 at `, z[0] = Zn[1];
f := Z -> subs(z=Z, simplify(w*(z - Zn[1])^4)):
`Use f(z) ` = f(z);

`Find `, Int(exp(z^2)/(z-I)^4,z = C .. ``)

`Singularity of order  4  at  `, z[0] = I

`Use  f(z) ` = exp(z^2)

The integral of f(z) taken over C is obtained by applying the
Cauchy integral formula for derivatives.

> d3 := diff(f(z), z$3):
f3 := Z -> subs(z=Z, d3):
`f '''(z) ` = f3(z);
`f '''(z0) ` = f3(Zn[1]);
val := 2*Pi*I/3! * f3(Zn[1]):
Int(w ,z=C..``) = `2 Pi I/3! f '''(z0)`;
Int(w ,z=C..``) = val;

`f '''(z) ` = 12*z*exp(z^2)+8*z^3*exp(z^2)

`f '''(z0) ` = 4*I*exp(-1)

Int(exp(z^2)/(z-I)^4,z = C .. ``) = `2 Pi I/3! f ''...

Int(exp(z^2)/(z-I)^4,z = C .. ``) = -4/3*Pi*exp(-1)...

 

Remark. Sometimes Maple will form the list of values in a different order.

It is always necessary to visually inspect the above results before proceeding.

 

The integral of f(z) taken over C is obtained by applying the Cauchy integral formula.

> `f(z0) ` = f(Zn[1]);
`f(z0) ` = evalc(f(Zn[1]));
val := 2*Pi*I*f(Zn[1]):
Int(w ,z=C..``) = `2 Pi I f(z0)`;
Int(w ,z=C..``) = val;

`f(z0) ` = exp(-1)

`f(z0) ` = exp(-1)

Int(exp(z^2)/(z-I)^4,z = C .. ``) = `2 Pi I f(z0)`

Int(exp(z^2)/(z-I)^4,z = C .. ``) = 2*I*Pi*exp(-1)

 

 

Corollary 6.2 If f(z) is analytic in the domain D , then all derivatives `f '(z)`, `f ''(z)`, `...`, f^`(n)`*`(z)` exists and are analytic in D .

 

Remark. Corollary 6.2 is interesting, as it illustrates a big difference between real and complex functions. It is possible for a real function f(z) to have the property that `f '(z)` exists everywhere in a domain D , but `f ''(z)` exists nowhere. This corollary states that if a complex function f(z) has the property that `f '(z)` exists everywhere in a domain D , then, remarkably, all derivatives of f(z) exist in D .

 

Corollary 6.3 If u(x,y) is a harmonic function at each point `(x,y)` in the domain D , then all partial derivatives diff(u(x,y),x), diff(u(x,y),y), diff(u(x,y),x,x), d... exists and are harmonic functions.

 

End of Section 6.5.