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/


CHAPTER 3 ANALYTIC and HARMONIC FUNCTIONS

Section 3.1 Differentiable Functions

Does the notion of a derivative of a complex function make sense? If so, how should it be defined, and what does it represent? These and other questions will be the focus of the next few sections.

Using our imagination, we take our lead from elementary Calculus and define the
derivate of f at z[0] , written f '( z[0] ), by

 

f '( z[0] ) = Limit((f(z)-f(z[0]))/(z-z[0]),z = z[0]) ,

 

provided that the limit exists. When this happens, we say that the function f is differentiable at z[0] . If we write Delta*z = z-z[0] , then this definition can be expressed in the form

 

f '( z[0] ) = Limit((f(z[0]+Delta*z)-f(z[0]))/(Delta*z),Delta*z =... .

 


Example 3.1, Page 94. Use the limit definition to find the derivative of f(z) = z^3 .

 

> df:='df': dz:='dz': f:='f': z:='z': z0:='z0':
f := z ->z^3:
`f(z) ` = f(z); ` `;
`Form the difference quotient.`;
dfdz := (f(z) - f(z0))/(z - z0):
df/dz = dfdz; ` `;

`f(z) ` = z^3

` `

`Form the difference quotient.`

df/dz = (z^3-z0^3)/(z-z0)

` `

Substitution of z[0] for z is indeterminate.

> subs(z=z0, dfdz);

 

Error, numeric exception: division by zero

The derivative is the limit of the difference quotient.

> Limit(df/dz, z=z0) = limit(dfdz, z=z0);

Limit(df/dz,z = z0) = 3*z0^2

The difference quotient can be simplified before taking a limit.

> df/dz = dfdz;
Q := simplify(dfdz):
df/dz = Q;
`f '(z0) = `, Limit(df/dz, z=z0) = subs(z=z0, Q);

df/dz = (z^3-z0^3)/(z-z0)

df/dz = z^2+z0*z+z0^2

`f '(z0) = `, Limit(df/dz,z = z0) = 3*z0^2

Which is known to be the formula for the derivative of f(z) .

 


Example 3.2, Page 94. Use the limit definition to show that f(z) = conjugate(z) is NOWHERE differentiable.

> f:='f': z:='z':
f := z -> conjugate(z):
`f(z) ` = f(z);

`f(z) ` = conjugate(z)

Form the difference quotient using a change in dx in x .

> dx:='dx': dX:='dX': dy:='dy': dY:='dY':
X:='X': X0:='X0': Y:='Y': Y0:='Y0':
var := {dX='dx',dY='dy',X='x',X0='x0',Y='y',Y0='y0'}:
assume(dX,real); assume(dY,real);
assume(X,real); assume(Y,real);
assume(X0,real); assume(Y0,real);
dFdX := (f(X0 + dX + I*Y0) - f(X0 + I*Y0))/
(X0 + dX + I*Y0 - X0 - I*Y0):
dfdx := subs(var, dFdX):
`f(z0+dx) ` = f(z0+dx),` and `,
`f(z0) ` = f(z0); ` `;
`f(z0+dx) - f(z0)`/dx = dfdx;

`f(z0+dx) ` = conjugate(z0+dx), `  and  `, `f(z0) `...

`  `

`f(z0+dx) - f(z0)`/dx = 1

And this limit is easy to compute:

> lim1 := limit(dfdx , dx=0):
Limit(df/dx, x=0) = lim1;

Limit(df/dx,x = 0) = 1

Form the difference quotient using a change in dy in y .

> dFdY := (f(X0 + I*(Y0+dY)) - f(X0 + I*Y0))/
(X0 + I*(Y0+dY) - X0 - I*Y0):
dfdy := subs(var, dFdY):
`f(z0+Idy) - f(z0)`/`I dy` = dfdy;
`f(z0+Idy) - f(z0)`/`I dy` = simplify(dfdy);

`f(z0+Idy) - f(z0)`/`I dy` = (-I*(y0+dy)+I*y0)/(I*(...

`f(z0+Idy) - f(z0)`/`I dy` = -1

And this limit is easy to compute:

> lim2 := limit(dfdy , dy=0):
Limit(df/dy, dy=0) = lim2;

Limit(df/dy,dy = 0) = -1

The two limits are different for any point z[0] .

> Limit(df/dx, dx=0) <> Limit(df/dy, dy=0);
lim1 <> lim2;

Limit(df/dx,dx = 0) <> Limit(df/dy,dy = 0)

1 <> -1

>

Since the two limits are different for any point z = z[0] , diff(f(z),z) does NOT exist for any z = z[0] .

 

 

Definition 3.1: Analytic

 

We say that the complex function f is analytic at the point z[0] provided there is some epsilon >0 such that `f '`(z) exists for all z*epsilon*D[epsilon](z[0]) . In other words, f must be differentiable not only at z[0] , but also at all points in some epsilon -neighborhood of z[0] .

 

Definition 3.1: Entire

 

If f is analytic on the whole complex plane then f said to be entire .

 

The rules for differentiation on page 96.

 

> c:='c': f:='f': F:='F': z:='z':
F := z -> c*f(z):
print(`Derivative of a scalar multiple.`);
`F(z)` = F(z);
`F '(z)` = diff(F(z), z);

`Derivative of a scalar multiple.`

`F(z)` = c*f(z)

`F '(z)` = c*diff(f(z),z)

> a:='a': b:='b': f:='f': g:='g': H:='H': z:='z':
H := z -> a*f(z) + b*g(z):
print(`Derivative of a linear combination multiple.`);
`H(z)` = H(z);
`H '(z)` = diff(H(z), z);

`Derivative of a linear combination multiple.`

`H(z)` = a*f(z)+b*g(z)

`H '(z)` = a*diff(f(z),z)+b*diff(g(z),z)

> f:='f': g:='g': P:='P': z:='z':
P := z -> f(z)*g(z):
print(`Derivative of the product of two functions.`);
`P(z)` = P(z);
`P '(z)` = diff(P(z), z);

`Derivative of the product of two functions.`

`P(z)` = f(z)*g(z)

`P '(z)` = diff(f(z),z)*g(z)+f(z)*diff(g(z),z)

> f:='f': F:='F': g:='g': Q:='Q': z:='z':
Q := z -> f(z)/g(z):
print(`Derivative of the quotient of two functions.`);
`Q(z)` = Q(z);
`Q '(z)` = normal(diff(Q(z), z));

`Derivative of the quotient of two functions.`

`Q(z)` = f(z)/g(z)

`Q '(z)` = (diff(f(z),z)*g(z)-f(z)*diff(g(z),z))/g(...

> f:='f': g:='g': H:='H': z:='z':
H := z -> f(g(z)):
print(`Derivative of the composition of two functions.`);
`H(z)` = H(z);
`H '(z)` = diff(H(z), z);

`Derivative of the composition of two functions.`

`H(z)` = f(g(z))

`H '(z)` = D(f)(g(z))*diff(g(z),z)

 

 

Example 3.3, Page 96. Use the rules to find the derivative of f(z) = z^2+i*2*z+3 ,
and then use this result to find the derivative of
f(z)^4 = (z^2+i*2*z+3)^4 .

> f:='f': g:= 'g': h:='h': w:='w': z:='z':
f := z -> z^2 + I*2*z + 3:
g := z -> z^4:
h := z -> g(f(z)):
`f(z)` = f(z);
`f '(z)` = diff(f(z), z); ` `;
`g(z)` = g(z);
`g '(z)` = diff(g(z), z); ` `;
`h(z) = g(f(z))` = g(f(z));
`h '(z)` = diff(h(z), z);

`f(z)` = z^2+2*I*z+3

`f '(z)` = 2*z+2*I

` `

`g(z)` = z^4

`g '(z)` = 4*z^3

` `

`h(z) = g(f(z))` = (z^2+2*I*z+3)^4

`h '(z)` = 4*(z^2+2*I*z+3)^3*(2*z+2*I)

 

 

Theorem 3.1 (Differentiable implies continuous)

 

If f is differentiable at z[0] , then f is continuous at z[0] .

 

Theorem 3.2 (L'Hopital's Rule)

 

Assume f and g are analytic at z[0] . If we have f(z[0]) = 0 , g(z[0]) = 0 , and `g'(`*z[0]*`)` <> 0 , then

 

Limit(f(z)/g(z),z = z[0]) = Limit(`f '(z)`/`g '(z)`,z = z[0]) .

 

Example about L'Hopital's Rule, Page 98.

> f:='f': g:='g': z:='z':
f := z -> z^2 + z - 1 - 3*I:
g := z -> z^2 - 2*z + 2:
`f(z) ` = f(z);
`g(z) ` = g(z); ` `;
`f(z)/g(z) ` = f(z)/g(z);

`f(z) ` = z^2+z-1-3*I

`g(z) ` = z^2-2*z+2

` `

`f(z)/g(z) ` = (z^2+z-1-3*I)/(z^2-2*z+2)

Direct substitution of z = 1+i is indeterminate.

> f(1+I)/g(1+I);

 

Error, numeric exception: division by zero

Find limit(f(z)/g(z),z = 1+i) .

> Limit(f(z)/g(z), z=1+I) = limit(f(z)/g(z), z=1+I);

Limit((z^2+z-1-3*I)/(z^2-2*z+2),z = 1+I) = 1-3/2*I

Which agrees with the L'Hopital rule computation.

> f1:='f1': g1:='g1': Z:='Z':
f1 := z -> subs(Z=z, diff(f(Z), Z)):
g1 := z -> subs(Z=z, diff(g(Z), Z)):
`f '(z)/g '(z)` = f1(z)/g1(z); ` `;
`f '(1+I)` = f1(1+I);
`g '(1+I)` = g1(1+I); ` `;
`f '(1+I)/g '(1+I)` = f1(1+I)/g1(1+I);

`f '(z)/g '(z)` = (2*z+1)/(2*z-2)

` `

`f '(1+I)` = 3+2*I

`g '(1+I)` = 2*I

` `

`f '(1+I)/g '(1+I)` = 1-3/2*I

 

 

Theorem 3.3 (The first identity theorem)

 

Let f = u+i*v be an analytic function in the domain D . Suppose for all z*epsilon*D that abs(f(z)) = K , where K is a constant.

Then f constant in D .

 

End of Section 3.1.