Pde Book
Pde Book
Pde Book
Preface
Partial differential equations are often used to construct models of the most
basic theories underlying physics and engineering. The goal of this book is to
develop the most basic ideas from the theory of partial differential equations,
and apply them to the simplest models arising from the above mentioned
fields.
It is not easy to master the theory of partial differential equations. Unlike
the theory of ordinary differential equations, which relies on the fundamental
existence and uniqueness theorem, there is no single theorem which is central
to the subject. Instead, there are separate theories used for each of the major
types of partial differential equations that commonly arise.
It is worth pointing out that the preponderance of differential equations arising in applications, in science, in engineering, and within mathematics itself,
are of either first or second order, with the latter being by far the most prevalent. We will mainly cover these two classes of PDEs.
This book is intended for a first course in partial differential equations at
the advanced undergraduate level for students in engineering and physical
sciences. It is assumed that the student has had the standard three semester
calculus sequence, and a course in ordinary differential equations.
Marcel B Finan
August 2009
ii
PREFACE
Contents
Preface
. 71
. 75
.
.
.
.
.
.
.
.
.
.
35
36
41
41
50
59
. 65
87
88
93
100
109
120
CONTENTS
16 Fourier Sines Series and Fourier Cosines Series . . . . . . . . .
17 Separation of Variables for PDEs . . . . . . . . . . . . . . . .
17.1 Second Order Linear Homogenous ODE with Constant
Coefficients . . . . . . . . . . . . . . . . . . . . . . .
17.2 The Method of Separation of Variables for PDEs . . . .
18 Solutions of the Heat Equation by the Separation of Variables
Method . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
19 Elliptic Type: Laplaces Equations in Rectangular Domains . .
20 Laplaces Equations in Circular Regions . . . . . . . . . . . . .
. 133
. 140
. 140
. 141
. 147
. 154
. 165
233
Index
289
1 Basic Concepts
The goal of this section is to introduce the reader to the basic concepts and
notations that will be used in the remainder of this book.
We start this section by reviewing the concept of partial derivatives and the
chain rule of functions in two variables.
Let u(x, y) be a function of the independent variables x and y. The first
derivative of u with respect to x is defined by
u(x + h, y) u(x, y)
u
(x, y) = lim
h0
x
h
provided that the limit exists.
Likewise, the first derivative of u with respect to y is defined by
ux (x, y) =
uy (x, y) =
u
u(x, y + h) u(x, y)
(x, y) = lim
h0
y
h
and
2u
ux (x, y + h) ux (x, y)
uyx (x, y) =
(x, y) = lim
h0
yx
h
1
provided that the limits exist.
An important formula of differentiation is the so-called chain rule. If u =
u(x, y), where x = x(s, t) and y = y(s, t), then
u
u x u y
=
+
.
s
x s y s
Likewise,
u
u x u y
=
+
.
t
x t
y t
1
If uxy and uux are continuous then uxy (x, y) = uyx (x, y).
1 BASIC CONCEPTS
Example 1.1
Compute the partial derivatives indicated:
(a) y
(y 2 sin xy)
(b)
2
[ex+y ]2
x2
Solution.
(a) We have
(b) We have
4e2(x+y)
2(x+y)
2
x+y 2
[ex+y ]2 = x
e
= 2e2(x+y) . Thus, x
] = x
2e2(x+y) =
2 [e
x
Example 1.2
Suppose u(x, y) = sin (x2 + y 2 ), where x = tes and y = s + t. Find us (s, t)
and ut (s, t).
Solution.
We have
us (s, t) =ux xs + uy ys = 2x cos (x2 + y 2 )tes + 2y cos (x2 + y 2 )
=[2t2 e2s + 2(s + t)] cos [t2 e2s + (s + t)2 ].
Likewise,
ut (s, t) =ux xt + uy yt = 2x cos (x2 + y 2 )es + 2y cos (x2 + y 2 )
=[2te2s + 2(s + t)] cos [t2 e2s + (s + t)2 ]
A differential equation is an equation that involves an unknown scalar
function (the dependent variable) and one or more of its derivatives. For
example,
d2 y
dy
5 + 3y = 3
(1.1)
2
dx
dx
or
u 2 u 2 u
2 2 + u = 0.
(1.2)
t
x
y
If the unknown function is a function in one single variable then the differential equation is called an ordinary differential equation, abbreviated by
ODE. An example of an ordinary differential equation is Equation (1.1). In
contrast, when the unknown function is a function of two or more independent variables then the differential equation is called a partial differential
equation, in short PDE. Equation (1.2) is an example of a partial differential
equation. In this book we will be focusing on partial differential equations.
1 BASIC CONCEPTS
order
A first order partial differential equation is called quasi-linear if it can
be written in the form
a(x, y, u)ux + b(x, y, u)uy = c(x, y, u).
(1.3)
Note that linear and semi-linear partial differential equations are special cases
of quasi-linear equations. However, a quasi-linear PDE needs not be linear:
A partial differential equation that is not linear is called non-linear. For
example, u2x + 2uxy = 0 is non-linear. Note that this equation is quasi-linear
and semi-linear.
As for ODEs, linear PDEs are usually simpler to analyze/solve than nonlinear PDEs.
Example 1.6
Determine whether the given PDE is linear, quasi-linear, semi-linear, or nonlinear:
(a) xux + yuy = x2 + y 2 .
(b) uux + uy = 2.
(c) utt c2 uxx = f (x, t).
(d) ut + uux + uxxx = 0.
(e) u2tt + uxx = 0.
Solution.
(a) Linear, quasi-linear, semi-linear.
(b) Quasi-linear, non-linear.
(c) Linear, quasi-linear, semi-linear.
1 BASIC CONCEPTS
Example 1.8
Determine whether the equation is homogeneous or non-homogeneous:
(a) xux + yuy = x2 + y 2 .
(b) utt = c2 uxx .
(c) y 2 uxx + xuyy = 0.
Solution.
(a) Non-homogeneous because of x2 + y 2 .
(b) Homogeneous.
(c) Homogeneous
Practice Problems
Problem 1.1
Classify the following equations as either ODE or PDE.
(a) (y 000 )4 +
(b)
u
x
t2
(y 0 )2 +4
+ y u
=
y
= 0.
yx
.
y+x
(c) y 00 4y = 0.
Problem 1.2
Write the equation
uxx + 2uxy + uyy = 0
in the coordinates s = x, t = x y.
Problem 1.3
Write the equation
uxx 2uxy + 5uyy = 0
in the coordinates s = x + y, t = 2x.
Problem 1.4
For each of the following PDEs, state its order and whether it is linear
or non-linear. If it is linear, also state whether it is homogeneous or nonhomogeneous:
(a) uux + x2 uyyy + sin x = 0.
2
(b) ux + ex uy = 0.
(c) utt + (sin y)uyy et cos y = 0.
Problem 1.5
For each of the following PDEs, determine its order and whether it is linear
or not. For linear PDEs, state also whether the equation is homogeneous or
not. For non-linear PDEs, circle all term(s) that are not linear.
(a) x2 uxx + ex u = xuxyy .
(b) ey uxxx + ex u = sin y + 10xuy .
(c) y 2 uxx + ex uux = 2xuy + u.
(d) ux uxxy + ex uuy = 5x2 ux .
(e) ut = k 2 (uxx + uyy ) + f (x, y, t).
1 BASIC CONCEPTS
11
Problem 1.6
Which of the following PDEs are linear?
(a) Laplaces equation: uxx + uyy = 0.
(b) Convection (transport) equation: ut + cux = 0.
(c) Minimal surface equation: (1 + Zy2 )Zxx 2Zx Zy Zxy + (1 + Zx2 )Zyy = 0.
(d) Korteweg-Vries equation: ut + 6uux = uxxx .
Problem 1.7
Classify the following differential equations as ODEs or PDEs, linear or
non-linear, and determine their order. For the linear equations, determine
whether or not they are homogeneous.
(a) The diffusion equation for u(x, t) :
ut = kuxx .
(b) The wave equation for w(x, t) :
wtt = c2 wxx .
(c) The thin film equation for h(x, t) :
ht = (hhxxx )x .
(d) The forced harmonic oscillator for y(t) :
ytt + 2 y = F cos (t).
(e) The Poisson Equation for the electric potential (x, y, z) :
xx + yy + zz = 4(x, y, z)
where (x, y, z) is a known charge density.
(f) Burgers equation for h(x, t) :
ht + hhx = hxx .
Problem 1.8
Write down the general form of a linear second order differential equation of
a function in three variables.
13
2 e
2 2 t
(cos x sin x)
(2 cos x + 2 sin x) = 0,
The idea behind the name is due to the fact that integration is being used to finding
the solution.
Figure 2.1
Example 2.3
Find the general solution of uxy = 0.
Solution.
Integrating first we respect to y we find ux (x, y) = f (x), where f is an
arbitrary Rdifferentiable function. Integrating ux with respect to x we find
u(x, y) = f (x)dx + g(y), where g is an arbitrary differentiable function
Note that the general solution in the previous example involves two arbitrary
functions. In general, the general solution of a partial differential equation
is an expression that involves arbitrary functions. This is in contrast to the
general solution of an ordinary differential equation which involves arbitrary
constants.
Usually, a classical solution enjoys properties such as smootheness (i.e. differentiability) and continuity. However, in the theory of non-linear pdes,
there are solutions that do not require the smoothness property. Such solutions are called weak solutions or generalized solutions. For example,
u(x) = x is a classical solution to the differential equation uu0 = x. In contrast, u(x) = |x| is a generalized solution since it is not differentiable at 0.
In this book, the word solution will refer to a classical solution.
Example 2.4
Show that u(x, t) = t + 12 x2 is a classical solution to the PDE
ut = uxx .
(2.1)
15
Solution.
Assume that the domain of definition of u is D R2 . Since u, ut , ux , utx , uxx
exist and are continuous in D(i.e., u is smooth in D) and u satisfies equation
(2.1), we conclude that u is a classical solution to the given PDE
We next consider the structure of solutions to linear partial differential equations. To this end, consider the linear differential operator L as defined in
the previous section. The defining properties of linearity immediately imply
the key facts concerning homogeneous linear differential equations.
Theorem 2.1
The sum of two solutions to a homogeneous linear differential equation is
again a solution, as is the product of a solution by any constant.
Proof.
Let u1 , u2 be solutions, meaning that L[u1 ] = 0 and L[u2 ] = 0. Then, thanks
to linearity,
L[u1 + u2 ] = L[u1 ] + L[u2 ] = 0,
and hence their sum u1 + u2 is a solution. Similarly, if is any constant, and
u any solution, then
L[u] = L[u] = 0 = 0,
and so the scalar multiple u is also a solution
The following result is known as the superposition principle for homogeneous linear equations. It states that from given solutions to the equation
one can create many more solutions.
Theorem 2.2
If u1 , , un are solutions to a common homogeneous linear partial differential equation L[u] = 0, then the linear combination u = c1 u1 + + cn un is
a solution for any choice of constants c1 , , cn .
Proof.
The key fact is that, thanks to the linearity of L, for any sufficiently smooth
functions u1 , , un and any constants c1 , , cn ,
L[u] =L[c1 u1 + + cn un ] = L[c1 u1 + + cn1 un1 ] + L[cn un ]
= = L[c1 u1 ] + + L[cn un ] = c1 L[u1 ] + + cn L[un ].
17
if 0 < x, y < 1
if 0 < x < 1
if 0 < y < 1.
There are three types of boundary conditions which arise frequently in formulating physical problems:
1. Dirichlet Boundary Conditions: In this case, the dependent function u is prescribed on the boundary of the bounded domain. For example, if
the bounded domain is the rectangular plate 0 < x < L1 and 0 < y < L2 , the
boundary conditions u(0, y), u(L1 , y), u(x, 0), and u(x, L2 ) are prescribed.
The boundary conditions are called homogeneous if the dependent variable
is zero at any point on the boundary, otherwise the boundary conditions are
called nonhomogeneous.
2. Neumann Boundary Conditions: In this case, first partial derivatives are prescribed on the boundary of the bounded domain. For example,
the Neuman boundary conditions for a rod of length L, where 0 < x < L,
are of the form ux (0, t) = and ux (L, t) = , where and are constants.
x
.
(b) Show that sup{|u (x, t) u(x, t)| : x R, t > 0} = 2 . Thus, a small
change in the initial data leads to a small change in the solution. Hence, the
initial value problem is well-posed.
19
Solution.
(a) We have
u
t
2
u
t2
u
x
2 u
x2
2
Thus, tu2
(b) We have
2 u
x2
t
= sin
cos
x
t
sin
= sin
x
t
= cos
sin
x
t
= sin
sin
.
x
u (x, 0)
t
= sin
x
.
t
x
2
sin
: x R, t > 0}
sup{|u (x, t) u(x, t)| : x R, t > 0} = sup{sin
=2
A problem that is not well-posed is referred to as an ill-posed problem. We
illustrate this concept in the next example.
Example 2.6
For x R and t > 0 we consider the initial value problem
utt + uxx =0
u(x, 0) = ut (x, 0) =0.
Clearly, u(x, t) = 0 is a solution to this problem.
(a) Let 0 < << 1 be a very small number. Show that the function u (x, t) =
2 sin x sinh t , where
ex ex
sinh x =
2
is a solution to the problem
utt + uxx =0
u(x, 0) =0
ut (x, 0) = sin
x
Solution.
(a) We have
u
t
2 u
t2
u
x
2 u
x2
2
Thus, tu2 +
(b) We have
2 u
x2
t
= sin
cosh
x
t
sinh
= sin
x
t
= cos
sinh
x
t
= sin
sinh
.
x
u (x, 0)
t
= sin
x
.
x
sup{| u (x, 0) ut (x, 0)| : x R} = sup{ sin
: x R}
t
x
= sup{sin
: x R} =
and
x
t
: x R}
sup{|u (x, t) u(x, t)| : x R} = sup{sinh
sin
t
2
= sinh
.
2
(c) We have
t
lim sup{|u (x, t) u(x, t)| : x R} = lim sinh
= .
t
t
2
21
Practice Problems
Problem 2.1
Determine a and b so that u(x, y) = eax+by is a solution to the equation
uxxxx + uyyyy + 2uxxyy = 0.
Problem 2.2
Consider the following differential equation
tuxx ut = 0.
Suppose u(t, x) = X(x)T (t). Show that there is a constant such that
X 00 = X and T 0 = tT.
Problem 2.3
Consider the initial value problem
xux + (x + 1)yuy = 0, x, y > 1
u(1, 1) = e.
Show that u(x, y) =
xex
y
Problem 2.4
Show that u(x, y) = e2y sin (x y) is the solution to the initial value problem
ux + uy + 2u = 0 for x, y > 1
u(x, 0) = sin x.
Problem 2.5
Solve each of the following differential equations:
= 0 where u = u(x).
(a) du
dx
(b) u
= 0 where u = u(x, y).
x
Problem 2.6
Solve each of the following differential equations:
2
(a) ddxu2 = 0 where u = u(x).
2u
(b) xy
= 0 where u = u(x, y).
Problem 2.12
Consider the initial value problem
ut + uxx = 0, x R, t > 0
u(x, 0) = 1.
(a) Show that u(x, t) 1 is a solution to this problem.
n2 t
(b) Show that un (x, t) = 1 + e n sin nx is a solution to the initial value
problem
ut + uxx = 0, x R, t > 0
23
sin nx
.
n
25
26
(3.1)
where p(t) and g(t) are continuous on the open interval aR < t < b.
Since
p(t) is continuous, it has an antiderivative namely p(t)dt. Let (t) =
R
p(t)dt
e
. Multiply Equation (3.1) by (t) and notice that the left hand side of
the resulting equation is the derivative of a product. Indeed,
d
((t)y) = (t)g(t).
dt
Integrate both sides of the last equation with respect to t to obtain
Z
(t)y = (t)g(t)dt + C
Hence,
or
Z
C
1
(t)g(t)dt +
y(t) =
(t)
(t)
Z R
R
R
p(t)dt
y(t) = e
e p(t)dt g(t)dt + Ce p(t)dt
Notice that the second term of the previous expression is just the general
solution for the homogeneous equation
y 0 + p(t)y = 0
whereas the first term is a solution to the nonhomogeneous equation. That
is, the general solution to Equation (3.1) is the sum of a particular solution of
the nonhomogeneous equation and the general solution of the homogeneous
equation.
Example 3.1
Solve the initial value problem
y0
y
= 4t, y(1) = 5.
t
27
Solution.
We have p(t) = 1t so that (t) = 1t . Multiplying the given equation by the
integrating factor and using the product rule we notice that
0
1
y = 4.
t
Integrating with respect to t and then solving for y we find that the general
solution is given by
Z
y(t) = t 4dt + Ct = 4t2 + Ct.
Since y(1) = 5, we find C = 1 and hence the unique solution to the IVP is
y(t) = 4t2 + t, 0 < t <
Example 3.2
Find the general solution to the equation
2
y 0 + y = ln t, t > 0.
t
Solution.
R 2
The integrating factor is (t) = e t dt = t2 . Multiplying the given equation
by t2 to obtain
(t2 y)0 = t2 ln t.
Integrating with respect to t we find
Z
2
t y = t2 ln tdt + C.
The integral on the right-hand side is evaluated using integration by parts
3
with u = ln t, dv = t2 dt, du = dtt , v = t3 obtaining
t2 y =
t3
t3
ln t + C
3
9
Thus,
y=
t
t C
ln t + 2
3
9 t
28
Example 3.3
Solve
aux + buy + cu = 0
by using the change of variables s = ax + by and t = bx ay.
Solution.
By the Chain rule for functions of two variables, we have
ux =us sx + ut tx = aus + but
uy =us sy + ut ty = bus aut .
Substituting into the given equation, we find
us +
a2
c
u = 0.
+ b2
u(s, t) = f (t)e
cs
a2 +b2
29
Practice Problems
Problem 3.1
Solve the IVP: y 0 + 2ty = t, y(0) = 0.
Problem 3.2
Find the general solution: y 0 + 3y = t + e2t .
Problem 3.3
Find the general solution: y 0 + 1t y = 3 cos t, t > 0.
Problem 3.4
Find the general solution: y 0 + 2y = cos (3t).
Problem 3.5
Find the general solution: y 0 + (cos t)y = 3 cos t.
Problem 3.6
Given that the solution to the IVP ty 0 + 4y = t2 , y(1) = 31 exists on the
interval < t < . What is the value of the constant ?
Problem 3.7
Suppose that y(t) = Ce2t + t + 1 is the general solution to the equation
y 0 + p(t)y = g(t). Determine the functions p(t) and g(t).
Problem 3.8
Suppose that y(t) = 2et + et + sin t is the unique solution to the IVP
y 0 + y = g(t), y(0) = y0 . Determine the constant y0 and the function g(t).
Problem 3.9
Find the value (if any) of the unique solution to the IVP y 0 + (1 + cos t)y =
1 + cos t, y(0) = 3 in the long run?
Problem 3.10
Solve the initial value problem ty 0 = y + t, y(1) = 7.
Problem 3.11
Show that if a and are positive constants, and b is any real number, then
every solution of the equation
y 0 + ay = bet
has the property that y 0 as t . Hint: Consider the cases a = and
a 6= separately.
30
Problem 3.12
Solve the initial-value problem y 0 + y = et y 2 , y(0) = 1 using the substitution
1
u(t) = y(t)
Problem 3.13
Solve the initial-value problem ty 0 + 2y = t2 t + 1, y(1) =
1
2
Problem 3.14
Solve y 0 1t y = sin t, y(1) = 3. Express your answer in terms of the sine
Rt
integral, Si(t) = 0 sins s ds.
31
d
d
d
F (y) G(t) = [F (y) G(t)] = 0.
dt
dt
dt
It follows that
F (y) G(t) = C
which is equivalent to
Z
f (y)y dt =
Z
g(t)dt + C.
As you can see, the result is generally an implicit equation involving a function of y and a function of t. It may or may not be possible to solve this to
get y explicitly as a function of t. For an initial value problem, substitute the
values of t and y by t0 and y0 to get the value of C.
Remark 4.1
If F is a differentiable function of y and y is a differentiable function of t and
given by
both F and y are given then the chain rule allows us to find dF
dt
dF
dF dy
=
dt
dy dt
For separable equations, we are given f (y)y 0 = dF
and we are asked to find
dt
F (y). This process is referred to as reversing the chain rule.
32
Example 4.1
Solve the initial value problem y 0 = 6ty 2 , y(1) =
1
.
25
Solution.
Separating the variables and integrating both sides we obtain
Z 0
Z
y
dt = 6tdt
y2
or
d
dt
Z
1
dt = 6tdt.
y
Thus,
1
= 3t2 + C.
y(t)
1
Since y(1) = 25
, we find C = 28. The unique solution to the IVP is then
given explicitly by
1
y(t) =
28 3t2
Example 4.2
Solve the IVP yy 0 = 4 sin (2t), y(0) = 1.
Solution.
This is a separable differential equation. Integrating both sides we find
Z
Z
d y2
dt = 4 sin (2t)dt.
dt 2
Thus,
y 2 = 4 cos (2t) + C.
Since y(0) = 1, we find C = 5. Now, solving explicitly for y(t) we find
y(t) = 4 cos t + 5.
Practice Problems
Problem 4.1
Solve the (separable) differential equation
2 ln y 2
y 0 = tet
Problem 4.2
Solve the (separable) differential equation
y0 =
t2 y 4y
.
t+2
Problem 4.3
Solve the (separable) differential equation
ty 0 = 2(y 4).
Problem 4.4
Solve the (separable) differential equation
y 0 = 2y(2 y).
Problem 4.5
Solve the IVP
y0 =
Problem 4.6
Solve the IVP:
4 sin (2t)
, y(0) = 1.
y
yy 0 = sin t, y( ) = 2.
2
Problem 4.7
Solve the IVP:
y 0 + y + 1 = 0, y(1) = 0.
Problem 4.8
Solve the IVP:
y 0 ty 3 = 0, y(0) = 2.
33
34
Problem 4.9
Solve the IVP:
Problem 4.10
Solve the IVP:
y 0 = 1 + y 2 , y( ) = 1.
4
1
y 0 = t ty 2 , y(0) = .
2
Problem 4.11
Solve the equation 3uy + uxy = 0 by using the substitution v = uy .
Problem 4.12
Solve the IVP
(2y sin y)y 0 = sin t t, y(0) = 0.
Problem 4.13
State an initial value problem, with initial condition imposed at t0 = 2,
having implicit solution y 3 + t2 + sin y = 4.
Problem 4.14
Can the differential equation
dy
= x2 xy
dx
be solved by the method of separation of variables? Explain.
35
36
(5.1)
(5.2)
then we say that the equation is quasi-linear. The following are examples
of quasi-linear equations:
uux + uy + cu2 = 0
x(y 2 + u)ux y(x2 + u)uy = (x2 y 2 )u.
If Equation (5.1) can be written in the form
a(x, y)ux + b(x, y)uy = c(x, y, u)
(5.3)
then we say that the equation is semi-linear. The following are examples
of semi-linear equations:
xux + yuy = u2 + x2
(x + 1)2 ux + (y 1)2 uy = (x + y)u2 .
If Equation (5.1) can be written in the form
a(x, y)ux + b(x, y)uy + c(x, y)u = d(x, y)
(5.4)
then we say that the equation is linear. Examples of linear equations are:
xux + yuy = cu
(y z)ux + (z x)uy + (x y)uz = 0.
A first order pde that is not linear is said to be non-linear. Examples of
non-linear equations are:
ux + cu2y = xy
37
u2x + u2y = c.
First order partial differential equations are classified as either linear or nonlinear. Clearly, linear equations are a special kind of quasi-linear equation
(5.2) if a and b are functions of x and y only and c is a linear function of u.
Likewise, semi-linear equations are quasilinear equations if a and b are functions of x and y only. Also, semi-linear equations (5.3) reduces to a linear
equation if c is linear in u.
A linear first order partial differential equation is called homogeneous if
d(x, y) 0 and non-homogeneous if d(x, y) 6= 0. Examples of linear homogeneous equations are:
xux + yuy = cu
(y z)ux + (z x)uy + (x y)uz = 0.
Examples of non-homogeneous equations are:
ux + (x + y)uy u = ex
yux + xuy = xy.
Recall that for an ordinary linear differential equation, the general solution
depends mainly on arbitrary constants. Unlike ODEs, in linear partial differential equations, the general solution depends on arbitrary functions.
Example 5.1
Solve the equation ut (x, t) = 0.
Solution.
The general solution is given by u(x, t) = f (x) where f is an arbitrary differentiable function of x
Example 5.2
Consider the transport equation
aut (x, t) + bux (x, t) = 0
where a and b are constants. Show that u(x, t) = f (bt ax) is a solution
to the given equation, where f is an arbitrary differentiable function in one
variable.
Solution.
Let v(x, t) = bt ax. Using the chain rule we see that ut (x, t) = bfv (v) and
ux (x, t) = afv (v). Hence, aut (x, t) + bux (x, t) = abfv (v) abfv (v) = 0
38
Practice Problems
Problem 5.1
Classify each of the following PDE as linear, quasi-linear, semi-linear, or nonlinear.
(a) xux + yuy = sin (xy).
(b) ut + uux = 0
(c) u2x + u3 u4y = 0.
(d) (x + 3)ux + xy 2 uy = u3 .
Problem 5.2
Show that u(x, y) = ex f (2x y), where f is a differentiable function of one
variable, is a solution to the equation
ux + 2uy u = 0.
Problem 5.3
39
Problem 5.6
Find a relationship between a and b if u(x, y) = f (ax+by) is a solution to the
equation 3ux 7uy = 0 for any differentiable function f such that f 0 (x) 6= 0
for all x.
Problem 5.7
Reduce the partial differential equation
aux + buy + cu = 0
to a first order ODE by introducing the change of variables s = bx ay and
t = x.
Problem 5.8
Solve the partial differential equation
ux + uy = 1
by introducing the change of variables s = x y and t = x.
Problem 5.9
Show that u(x, y) = e4x f (2x 3y) is a solution to the first-order PDE
3ux + 2uy + 12u = 0.
Problem 5.10
Derive the general solution of the PDE
aut + bux = u, b 6= 0
by using the change of variables v = ax bt and w = x.
Problem 5.11
Derive the general solution of the PDE
aux + buy = 0, a 6= 0
by using the change of variables s = bx ay and t = x.
40
Problem 5.12
Write the equation
ut + cux + u = f (x, t), c 6= 0
in the coordinates v = x ct, w = x.
Problem 5.13
Suppose that u(x, t) = w(x ct) is a solution to the PDE
xux + tut = Au
where A and c are constants. Let v = x ct. Write the differential equation
with unknown function w(v).
41
Figure 6.1.1
Indeed, the work W is given by the expression
42
where is the angle between the two vectors as shown in Figure 6.1.2.
Figure 6.1.2
The above definition is the geometric definition of the dot product. We
next provide an algebraic way for computing the dot product. Indeed, let
~u = u1~i + u2~j + u3~k and ~v = v1~i + v2~j + v3~k. Then ~v ~u = (v1 u1 )~i + (v2
u2 )~j + (v3 u3 )~k. Moreover, we have ||~u||2 = u21 + u22 + u23 , ||~v ||2 = v12 + v22 + v32
and
||~v ~u||2 =(v1 u1 )2 + (v2 u2 )2 + (v3 u3 )2
=v12 2v1 u1 + u21 + v22 2v2 u2 + u22 + v32 2v3 u3 + u23 .
Now, applying the Law of Cosines to Figure 6.1.3 we can write
||~v ~u||2 = ||~u||2 + ||~v ||2 2||~u|| ||~v || cos .
Thus, by substitution we obtain
v12 2v1 u1 +u21 +v22 2v2 u2 +u22 +v32 2v3 u3 +u23 = u21 +u22 +u23 +v12 +v22 +v32 2||~u|| ||~v || cos
or
||~u|| ||~v || cos = u1 v1 + u2 v2 + u3 v3
so that we can define the dot product algebraically by
~u ~v = u1 v1 + u2 v2 + u3 v3 .
Figure 6.1.3
1 ~
i
2
1 ~
j
2
1 ~
k
2
43
Solution.
We have
1 1
1
1
1
1
1 1
~u ~v = + + 1 = + + = 2.
2 2
2 2
2
2 2 2 2
2
We also have
2
2
2
1
1
1
3
||~u|| =
+
+
=
2
2
2
2
2 2
1
3
1
+
+1= .
||~v ||2 =
2
2
2
2
Thus,
~u ~v
2 2
cos =
=
.
||~u|| ||~v ||
3
Hence,
= cos1
!
2 2
0.34 rad 19.5
3
Remark 6.1.1
The algebraic definition of the dot product extends to vectors with any number of components.
Next, we discuss few properties of the dot product.
Theorem 6.1.1
For any vectors ~u, ~v , and w
~ and any scalar we have
(i) Commutative law: ~u ~v = ~v ~u.
(ii) Distributive law: (~u + ~v ) w
~ = ~u w
~ + ~v w.
~
(iii) ~u (~v ) = (~u) ~v = (~u ~v ).
(iv) Magnitude: ||~u||2 = ~u ~u.
(v) Two nonzero vectors ~u and ~v are orthogonal or perpendicular if and
only if ~u ~v = 0.
(vi)) Two nonzero vectors ~u and ~v are parallel if and only if ~u ~v = ||~u|| ||~v ||.
(vii) ~0 ~v = ~0.
44
Proof.
Write ~u = u1~i + u2~j + u3~k, ~v = v1~i + v2~j + v3~k, and w
~ = w1~i + w2~j + w3~k.
Then
(i) ~u ~v = u1 v1 + u2 v2 + u2 v3 = v1 u1 + v2 u2 + v3 u3 = ~v ~u since product of
numbers is commutative.
(ii) (~u + ~v ) w
~ = ((u1 + v1 )~i + (u2 + v2 )~j + (u2 + v3 )~k) (w1~i + w2~j + w3~k) =
(u1 + v1 )w1 + (u2 + v2 )w2 + (u3 + v3 )w3 = u1 w1 + u2 w2 + u3 w3 + v1 w1 + v2 w2 +
v3 w3 = ~u w
~ + ~v w.
~
(iii) ~u (~v ) = (u1~i+u2~j +u3~k)(v1~i+v2~j +v3~k) = u1 v1 +u2 v2 +u3 v3 =
(u1 v1 + u2 v2 + u3 v3 ) = (~u ~v ).
(iv) ||~u||2 = ~u ~u cos 0 = ~u ~u.
(v) If ~u and ~v are perpendicular then the cosine of their angle is zero and so
the dot product is zero. Conversely, if the dot product of the two vectors is
zero then the cosine of their angle is zero and this happens only when the
two vectors are perpendicular.
(vi) If ~u and ~v are parallel then the cosine of their angle is either 1 or 1.
That is, ~u ~v = ||~u|| ||~v ||. Conversely, if ~u ~v = ||~u|| ||~v || then cos = 1
and this implies that either = 0 or = . In either case, the two vectors
are parallel.
(vii) In 3-D, ~0 =< 0, 0, 0 > and ~v =< a, b, c > so that ~0 ~v = (0 a)~i + (0
b)~j + (0 c)~k = ~0
Remark 6.1.2
Note that the unit vectors ~i, ~j, ~k associated with the coordinate axes satisfy
the equalities
~i ~i = ~j ~j = ~k ~k = 1 and ~i ~j = ~j ~k = ~i ~k = 0.
Example 6.1.2
(a) Show that the vectors ~u = 3~i 2~j and ~v = 2~i + 3~j are perpendicular.
(b) Show that the vectors ~u = 2~i + 6~j 4~k and ~v = 3~i 9~j + 6~k are parallel.
Solution.
(a) We have: ~u ~v = 3(2) 2(3) = 0. Hence ~u is perpendicular to ~v .
(b) We have:
cos =
45
Hence, = so that the two vectors are parallel. Another way to see that
the vectors are parallel is to notice that ~u = 23 ~v
Projection of a vector onto a line
The orthogonal projection of a vector along a line is obtained by taking a
vector with same length and direction as the given vector but with its tail on
the line and then dropping a perpendicular onto the line from the tip of the
vector. The resulting vector on the line is the vectors orthogonal projection
or simply its projection. See Figure 6.1.4.
Figure 6.1.4
Now, if ~u is a unit vector along the line of projection and if ~vparallel is the
vector projection of ~v onto ~u then
~vparallel = (||~v || cos )~u = (~v ~u)~u.
See Figure 6.1.5. Also, the component perpendicular to ~u is given by
~vperpendicular = ~v ~vparallel .
Figure 6.1.5
46
We call Comp~u~v = ~v ~u the the scalar projection of ~v onto ~u. We call the
vector Proj~u~v = ~vparallel the vector projection of ~v onto ~u.
It follows that the vector ~v can be written in terms of ~vparallel and ~vperpendicular
~v = ~vparallel + ~vperpendicular .
Example 6.1.3
Write the vector ~v = 3~i + 2~j 6~k as the sum of two vectors, one parallel,
and one perpendicular to w
~ = 2~i 4~j + ~k.
Solution.
~
=
Let ~u = ||w
w||
~
2 ~
i
21
4 ~
j
21
~vparallel = (~v ~u)~u =
1 ~
k.
21
Then,
6
8
6
32
8
16
~u = ~i + ~j ~k.
21
21
21
21
21
21
Also,
~vperpendicular
32 ~
8 ~
16 ~
i+ 2
j + 6 +
k
=~v ~vparallel = 3 +
21
21
21
79
10
118 ~
= ~i + ~j
k.
21
21
21
Hence,
~v = ~vparallel + ~vperpendicular
Example 6.1.4
Find the scalar projection and vector projection of ~u =< 1, 1, 2 > onto
~v =< 2, 3, 1 > .
Solution.
We have
1(2) + (1)(3) + 2(1)
3
~u ~v
= p
=
2
2
2
||~v ||
14
(2) + 3 + 1
~u ~v ~v
Proj~v ~u =
||~v || ||~v ||
3
2
9
3
= ~v = ~i + ~j + ~k
14
7
14
14
comp~v ~u =
Applications
As pointed out earlier in the section, scalar products are used in Physics.
For instance, in finding the work done by a force applied on an object.
47
Example 6.1.5
A wagon is pulled a distance of 100 m along a horizontal path by a constant
force of 70 N. The handle of the wagon is held at an angle of 35 above the
horizontal. Find the work done by the force.
Solution.
The work done is
W = F d cos 35 = 70(100) cos 35 5734 J
48
Practice Problems
Problem 6.1.1
Find ~a ~b where ~a =< 4, 1, 41 > and ~b =< 6, 3, 8 > .
Problem 6.1.2
Find ~a ~b where ||~a|| = 6, ||~b|| = 5 and the angle between the two vectors is
120 .
Problem 6.1.3
If ~u is a unit vector, find ~u ~v and ~u w
~ using the figure below.
Problem 6.1.4
Find the angle between the vectors ~a =< 4, 3 > and ~b =< 2, 1 > .
Problem 6.1.5
Find the angle between the vectors ~a =< 4, 3, 1 > and ~b =< 2, 0, 1 > .
Problem 6.1.6
Determine whether the given vectors are orthogonal, parallel, or neither.
(a) ~a =< 5, 3, 7 > and ~b =< 6, 8, 2 > .
(b) ~a =< 4, 6 > and ~b =< 3, 2 > .
(c) ~a = ~i + 2~j + ~k and ~b = 3~i + 4~j ~k.
(d) ~a = 2~i + 6~j 4~k and ~b = 3~i 9~j + 6~k.
Problem 6.1.7
Use vectors to decide whether the triangle with vertices P (1, 3, 2), Q(2, 0, 4),
and R(6, 2, 5) is right-angled.
Problem 6.1.8
Find a unit vector that is orthogonal to both ~i + ~j and ~i + ~k.
49
Problem 6.1.9
Find the acute angle between the lines 2x y = 3 and 3x + y = 7.
Problem 6.1.10
Find the scalar and vector projections of the vector ~b =< 1, 2, 3 > onto
~a =< 3, 6, 2 > .
Problem 6.1.11
If ~a =< 3, 0, 1 >, find a vector ~b such that comp~a~b = 2.
Problem 6.1.12
Find the work done by a force F~ = 8~i 6~j + 9~k that moves an object from
the point (0, 10, 8) to the point (6, 12, 20) along a straight line. The distance
is measured in meters and the force in newtons.
Problem 6.1.13
A sled is pulled along a level path through snow by a rope. A 30-lb force
acting at an angle of 40 above the horizontal moves the sled 80 ft. Find the
work done by the force.
50
Figure 6.2.1
The slope of the tangent line T to C at the point P is the rate of change
of z in the direction of ~u. Let Q(x, y, z) be an arbitrary point on C and let
P 0 (x0 , y0 , 0) and Q0 (x, y, 0) be the orthogonal projection of P and Q respec
tively onto the xyplane. Then the vectors P 0 Q0 =< xx0 , y y0 > is paral
lel to ~u so that P 0 Q0 = h~u for some scalar h. Hence, x = x0 + ha, y = y0 + hb
and
f (x0 + ha, y0 + hb) f (x0 , y0 )
f (x, y) f (x0 , y0 )
=
.
h
h
If we take the limit of the above average rate as h 0, we obtain the rate
of change of z(with respect to distance) in the direction of ~u, which is called
the directional derivative of f at (x0 , y0 ) in the direction of ~u. We write
f (x0 + ah, y0 + bh) f (x0 , y0 )
.
h0
h
51
g 0 (0) = lim
f dx f dy
+
= fx (x, y)a + fy (x, y)b.
x dh y dh
(6.2.1)
Example 6.2.1
2
Find u~v (4, 0) if u(x, y) = x + y and ~v = 12 , 23 .
Solution.
We have
1
u~v (4, 0) = ux (4, 0)
+ uy (4, 0)
2
!
3
1
=
2
2
52
Example 6.2.2
Let F (x, y, z) = u(x, y) z. Find F (x, y, z).
Solution.
We have
F (x, y, z) = ux~i + uy~j ~k
Example 6.2.3
Find the gradient vector of f (x, y, z) = (2x 3y + 5z)5 .
Solution.
We have
fx (x, y, z) =10(2x 3y + 5z)4
fy (x, y, z) = 15(2x 3y + 5z)4
fz (x, y, z) =25(2x 3y + 5z)4 .
Thus,
f (x, y, z) = 5(2x 3y + 5z)4 [2~i 3~j + 5~k]
With the notation for the gradient vector, we can rewrite the expression
(6.2.1) for the directional derivative as
f~u (x0 , y0 ) = f (x0 , y0 ) ~u.
This expresses the directional derivative in the direction of ~u as the scalar
projection of the gradient vector onto ~u.
Maximizing the Directional Derivative
Suppose we have a function of two or three variables and we consider all
possible directional derivatives of f at a given point. These give the rates of
change of in all possible directions. We can then ask the questions: In which
of these directions does f change fastest and what is the maximum rate of
change? The answers are provided by the following theorem.
Theorem 6.2.2
The maximum value of the directional derivative of a function f (x, y) or
f (x, y, z) at a point (x, y) or (x, y, z) is ||f || and it occurs in the direction
of the gradient of f at that point.
53
Proof.
We have
f~u (x, y) = f ~u = ||f ||||~u|| cos = ||f || cos ,
where is the angle between f and ~u. The maximum value of cos is 1 and
this occurs when = 0. Therefore the maximum value of f~u is ||f || and it
occurs when = 0, that is, when ~u has the same direction as f
Example 6.2.4
Find the maximum rate of change of the function u(x, y) = 50 x2 2y 2 at
the point (1, 1).
Solution.
The maximum rate of change occurs in the direction of the gradient vector:
u(1, 1) = ux (1, 1)~i + uy (1, 1)~j = 2~i + 4~j.
The maximum rate of change at (1, 1) is
||u(1, 1)|| =
(2)2 + 42 = 2 5
54
Figure 6.2.2
Now, for a function in two variables u(x, y), the equation u(x, y) = C is
called a level curve of u( a level surface of u(x, y, z)). The level curves
u(x, y) = C are just the traces of the graph of u(x, y) in the horizontal plane
z = C projected down to the xyplane.
An important property of the gradient of u is that it is normal to a level
surface of u at every point. To see this, let S be the level surface f (x, y, z) = k
and P0 (x0 , y0 , z0 ) be a point on S. Let C be any curve on S that passes
through P0 . We can describe C in parametric form x = x(t), y = y(t), and
z = z(t). Any point on C satisfies f (x(t), y(t), z(t)) = k. Differentiating both
sides of this equation with respect to t we find by means of the Chain Rule
fx (x, y, z)x0 (t) + fy (x, y, z)y 0 (t) + fz (x, y, z)z 0 (t) = 0
which can be written as f r0~(t) = 0. This means that the gradient is normal
to a level surface (respectively a level curve). See Figure 6.2.3.
Figure 6.2.3
55
Figure 6.2.4
Vector Fields and Integral Curves
In vector calculus, a vector field is a function F~ (x, y) (or F~ (x, y, z) in 3-D
space) that assigns a vector to each point of its domain as shown in Figure
6.2.5.
Figure 6.2.5
56
Creating vector fields manually is very tedious. Thus, vector fields are generally generated using computer softwares such as Mathematica, Maple, or
Mathlab.
Example 6.2.5
The gradient vector of a function is an example of a vector field called the
gradient vector field. Sketch the gradient vector field of the function
u(x, y) = x2 + y 2 .
Describe the level curves of u(x, y).
Solution.
The gradient vector field of the given function is
u(x, y) = 2x~i + 2y~j.
A level curve is defined by the equation
x2 + y 2 = C, C 0.
Thus, level curves are circles centered at the origin. Figure 6.2.6 shows the
gradient vector field as well as some of the level curves.
Figure 6.2.6
57
For example, at the point (1, 2), the corresponding vector in the vector field
is the vector with tail (1, 2) and tip (2, 4)
An integral curve of a vector field is a smooth curve5 such that F~ (x, y)
assigns a tangent vector at each point of . For example, the integral curves
of the vector field F~ (x, y) = y~i x~j are circles centered at the origin. See
Figure 6.2.7.
Figure 6.2.7
If ~r(t) is a parametrization of then r~0 (t) is continuous and r~0 (t) 6= ~0.
58
Practice Problems
Problem 6.2.1
Find the gradient of the function
F (x, y, z) = exyz + sin (xy).
Problem 6.2.2
Find the gradient of the function
F (x, y, z) = x cos
y
z
Problem 6.2.3
Describe the level surfaces of the function f (x, y, z) = (x 2)2 + (y 3)2 +
(z + 5)2 .
Problem 6.2.4
Find the directional derivative of u(x, y) = 4x2 + y 2 in the direction of ~a =
~i + 2~j at the point (1, 1).
Problem 6.2.5
Find the directional derivative of u(x, y, z) = x2 z +y 3 z 2 xyz in the direction
of ~a = ~i + 3~k at the point (x, y, z).
Problem 6.2.6
Find the maximum rate of change of the function u(x, y) = yexy at the point
(0, 2) and the direction in which this maximum occurs.
Problem 6.2.7
Find the gradient vector field for the function u(x, y, z) = ez ln (x2 + y 2 ).
59
(7.1)
(7.2)
60
Hence, the vectors r~0 (t) and ~v are parallel so these two vectors are proportional and this leads to the ODE system
dx
dt
dy
dt
du
dt
f (x, y, u)
(7.3)
or in differential form
dx
dy
du
=
=
.
(7.4)
a
b
f (x, y, u)
By solving the system (7.3) or (7.4), we are assured that the vector ~v is
tangent to the curve which in turn lies in the solution surface S. In our
context, integral curves are called characteristic curves or simply characteristics of the PDE (7.1). We call (7.3) the characteristic equations.
The projection of into the xyplane is called the projected characteristic curve.
Once we have found the characteristic curves, the surface S is the union of
these characteristic curves. In summary, by introducing these characteristic
equations, we have reduced our partial differential equation to a system of
ordinary differential equations. We can use ODE theory to solve the characteristic equations, then piece together these characteristic curves to form a
surface. Such a surface will provide us with a solution to our PDE.
Remark 7.1
dy
= ab one obtains the general solution h(x, y) = k1 where k1
Solving dx
is constant. Likewise, solving du
= fa one obtains the general solution
dx
j(x, y, u) = k2 where k2 is a constant. The constant k2 is a function of
k1 . For the sake of discussion, suppose that h(x, y) = k1 can be expressed as
= fa is being replaced by g(x, k1 ) so that the
y = g(x, k1 ). Then, the y in du
dx
constant in j(x, y, u) = k2 will depend on k1 .
Example 7.1
Find the general solution to aux + buy = 0 where a and b are constants with
a 6= 0.
Solution.
dy
From (7.3) we can write dx
= ab which yields bx ay = k1 for some arbitrary
du
constant k1 . From dx = 0 we find u(x, y) = k2 where k2 is a constant. That
is, u(x, y) is constant on . Since (0, ka1 , k2 ) is on , we have
u(x, y) = u(0,
k1
) = k2
a
61
Example 7.4
Find the general solution of the equation
xux + yuy = xeu , x > 0.
Solution.
We have a(x, y) = x, b(x, y) = y, and f (x, y, u) = xeu . So we have to solve
the system
dy
y du
= ,
= eu .
dx
x dx
62
From the first equation, we can use the separation of variables method to find
y = k1 x for some constant k1 . Solving the second equation by the method of
separation of variables, we find
eu x = k2 .
But k2 = g(k1 ) so that
u
e x = g(k1 ) = g
y
x
u = (1 + y) + f (y x)ex
Example 7.6
Find the general solution to x2 ux + y 2 uy = (x + y)u.
Solution.
Using properties of proportions6 we have
dx
dy
du
dx dy
= 2 =
= 2
.
2
x
y
(x + y)u
x y2
6
If
a
b
c
d
then
ab
b
cd
d .
Also,
a
b
c
d
e
f
a+c+e
b+d+f .
dy
dx
y2
x2
du
(x+y)u
d(xy)
x2 y 2
63
1
x
y1 = k1 .
we find
d(x y)
du
=
u
xy
which implies
u = k2 (x y) = f
1 1
x y
(x y)
Example 7.7
Find the solution satisfying yux + xuy = x2 + y 2 subject to the conditions
u(x, 0) = 1 + x2 and u(0, y) = 1 + y 2 .
Solution.
Solving the equation
have
dy
dx
x
y
du =y 1 (x2 + y 2 )dx
=ydx + x2 y 1 dx
y
=ydx + x2 y 1
dy
x
=ydx + xdy = d(xy).
Hence,
u(x, y) = xy + f (x2 y 2 ).
From u(x, 0) = 1 + x2 we find f (x) = 1 + x, x 0. From u(0, y) = 1 + y 2 we
find f (y) = 1 y, y 0. Hence, f (x) = 1 + |x| and
u(x, y) = xy + |x2 y 2 |
Remark 7.2
The method of characteristics discussed in this section applies as well to any
quasi-linear first order PDE. See Chapter 9.
64
Practice Problems
Problem 7.1
Solve ux + yuy = y 2 with the initial condition u(0, y) = sin y.
Problem 7.2
Solve ux + yuy = u2 with the initial condition u(0, y) = sin y.
Problem 7.3
Find the general solution of yux xuy = 2xyu.
Problem 7.4
Find the integral surface of the IVP: xux + yuy = u, u(x, 1) = 2 + e|x| .
Problem 7.5
Find the unique solution to 4ux + uy = u2 , u(x, 0) =
1
.
1+x2
Problem 7.6
2
Find the unique solution to e2y ux + xuy = xu2 , u(x, 0) = ex .
Problem 7.7
Find the unique solution to xux + uy = 3x u, u(x, 0) = tan1 x.
Problem 7.8
Solve: xux yuy = 0, u(x, x) = x4 .
Problem 7.9
Find the general solution of yux 3x2 yuy = 3x2 u.
Problem 7.10
Find u(x, y) that satisfies yux + xuy = 4xy 3 subject to the boundary conditions u(x, 0) = x4 and u(0, y) = 0.
8 LINEAR FIRST ORDER PDE: THE ONE DIMENSIONAL SPATIAL TRANSPORT EQUATIONS65
8 Linear First Order PDE: The One Dimensional Spatial Transport Equations
Modeling is the process of writing a differential equation to describe a physical situation. In this section we derive the one-dimensional spatial transport
equation and use the method of characteristics to solve it.
Linear Transport Equation for Fluid Flows
We shall describe the transport of a dissolved chemical by water that is traveling with uniform velocity c through a long thin tube G with uniform cross
section A. (The very same discussion applies to the description of the transport of gas by air moving through a pipe.) We assume the velocity c > 0
is in the (rightward) positive direction of the xaxis. We will also assume
that the concentration of the chemical is constant across the cross section
A located at x so that the chemical changes only in the xdirection and
thus the term one-dimensional spatial equation. This condition says that
the quantity of the chemical in a portion of the tube is the same but it is
traveling with time.
Let u(x, t) be a continuously differentiable function denoting the concentration of the chemical (i.e. amount of chemical per unit volume) at position
x and time t. Then at time t, the amount of chemical stored in a section of
the tube between positions a and x (see Figure 8.1) is given by the definite
integral
Z x
Au(s, t)ds.
a
Figure 8.1
66
Taking the derivative of both sides with respect to x and using the Fundamental Theorem of Calculus, we find
u(x, t) = u(x + ch, t + h).
Now, taking the derivative of this last equation with respect to h and using
the chain rule, we find
0 = ut (x + ch, t + h) + cux (x + ch, t + h).
Taking the limit of this last equation as h approaches 0 and using the fact
that ut and ux are continuous, we find
ut (x, t) + cux (x, t) = 0
(8.1)
for all (x, t). This equation is called the transport equation in one-dimensional
space. It is a linear, homogeneous first order partial differential equation.
Note that (8.1) can be written in the form
< 1, c > < ut , ux >= 0
so that the left-hand side of (8.1) is the directional derivative of u(t, x) at
(t, x) in the direction of the vector < 1, c > .
Solvability via the method of characteristics
We will use the method of characteristics discussed in Chapter 7 to solve
(8.1). The characteristic equations are
dt =
dx
du
=
.
c
0
8 LINEAR FIRST ORDER PDE: THE ONE DIMENSIONAL SPATIAL TRANSPORT EQUATIONS67
Solving the first equation, we find x ct = k1 . Solving the second equation
we find
u(x, t) = k2 = f (k1 ) = f (x ct).
One can check that this is indeed a solution to (8.1). Indeed, by using the
chain rule one finds
ut = cf 0 (x ct) and ux = f 0 (x ct).
Hence, by substituting these results into the equation one finds
ut + cux = cf 0 (x ct) + cf 0 (x ct) = 0.
The solution u(x, t) = f (x ct) is called the right traveling wave, since
the graph of the function f (x ct) at a given time t is the graph of f (x)
shifted to the right by the value ct. Thus, with growing time, the function
f (x) is moving without changes to the right at the speed c.
An initial value condition determines a unique solution to the transport equation as shown in the next example.
Example 8.1
2
Find the solution to ut 3ux = 0, u(x, 0) = ex .
Solution.
The characteristic equations lead to the ODEs
1 du
dt
= ,
= 0.
dx
3 dx
Solving the first equation, we find 3t + x = k1 . From the second equation, we
find u(x, t) = k2 = f (k1 ) = f (3t + x). From the initial condition, u(x, 0) =
2
f (x) = ex . Hence,
2
u(x, t) = e(3t+x)
Transport Equation with Decay
Recall from ODE that a function u is an exponential decay function if it
satisfies the equation
du
= u, < 0.
dt
68
(8.2)
8 LINEAR FIRST ORDER PDE: THE ONE DIMENSIONAL SPATIAL TRANSPORT EQUATIONS69
Practice Problems
Problem 8.1
Find the solution to ut + 3ux = 0, u(x, 0) = sin x.
Problem 8.2
Solve the equation aux + buy + cu = 0.
Problem 8.3
Solve the equation ux +2uy = cos (y 2x) with the initial condition u(0, y) =
f (y), where f : R R is a given function.
Problem 8.4
Show that the initial value problem ut + ux = x, u(x, x) = 1 has no solution.
Problem 8.5
Solve the transport equation ut + 2ux = 3u with initial condition u(x, 0) =
1
.
1+x2
Problem 8.6
Solve ut + ux 3u = t with initial condition u(x, 0) = x2 .
Problem 8.7
Show that the decay term u in the transport equation with decay
ut + cux + u = 0
can be eliminated by the substitution w = uet .
Problem 8.8 (Well-Posed)
Let u be the unique solution to the IVP
ut + cux = 0
u(x, 0) = f (x)
and v be the unique solution to the IVP
ut + cux = 0
u(x, 0) = g(x)
70
(9.1)
72
The required surface can be found as the union of integral curves, that is,
curves that are tangent to ~v at every point on the curve. If an integral curve
has a parametrization
~r(t) = x(t)~i + y(t)~j + u(t)~k
then the integral curve (i.e. the characteristic) is a solution to the ODE
system
dy
du
dx
= a(x, y, u),
= b(x, y, u),
= c(x, y, u)
(9.2)
dt
dt
dt
or in differential form
dx
dy
du
=
=
.
a(x, y, u)
b(x, y, u)
c(x, y, u)
(9.3)
dy
y(u2 +x2 )
du
.
(x2 +y 2 )u
Using a
y2
or
Hence, we find
dx
x
u2
dy
y
u2
+ x2
du
(x2 + y 2 )u
dx dy
du
=
.
x
y
u
yu
x
x
f (x2 + y 2 + u2 )
y
74
Practice Problem
Problem 9.1
Find the general solution of the PDE ln (y + u)ux + uy = 1.
Problem 9.2
Find the general solution of the PDE x(y u)ux + y(u x)uy = u(x y).
Problem 9.3
Find the general solution of the PDE u(u2 + xy)(xux yuy ) = x4 .
Problem 9.4
Find the general solution of the PDE (y + xu)ux (x + yu)uy = x2 y 2 .
Problem 9.5
Find the general solution of the PDE (y 2 + u2 )ux xyuy + xu = 0.
Problem 9.6
Find the general solution of the PDE ut + uux = x.
Problem 9.7
Find the general solution of the PDE (y u)ux + (u x)uy = x y.
Problem 9.8
Solve
x(y 2 + u)ux y(x2 + u)uy = (x2 y 2 )u.
Problem 9.9
Solve
1 x2 ux + uy = 0.
Problem 9.10
Solve
u(x + y)ux + u(x y)uy = x2 + y 2 .
(10.1)
Recall that the initial value problem of a first order ordinary differential
equation asks for a solution of the equation which has a given value at a
given point in R. The Cauchy problem for the PDE (10.1) asks for a solution
of (10.1) which has given values on a given curve in R2 . A precise statement
of the problem is given next.
Initial Value Problem or Cauchy Problem
Let C be a given curve in R2 defined parametrically by the equations
x = x0 (t), y = y0 (t)
where x0 , y0 are continuously differentiable functions on some interval I. Let
u0 (t) be a given continuously differentiable function on I. The Cauchy problem for (10.1) asks for a continuously differentiable function u = u(x, y)
defined in a domain R2 containing the curve C and such that:
(1) u = u(x, y) is a solution of (10.1) in .
(2) On the curve C, u equals the given function u0 (t), i.e.
u(x0 (t), y0 (t)) = u0 (t), t I.
(10.2)
We call C the initial curve of the problem, u0 (t) the initial data, and
(10.2) the initial condition or Cauchy data of the problem. See Figure
10.1.
76
Figure 10.1
If we view a solution u = u(x, y) of (10.1) as an integral surface of (10.1),
we can give a simple geometrical statement of the problem: Find a solution surface of (10.1) containing the curve described parametrically by the
equations
: x = x0 (t), y = y0 (t), u = u0 (t), t I.
Note that the projection of this curve in the xyplane is just the curve C.
The following theorem asserts that under certain conditions the Cauchy problem (10.1) - (10.2) has a unique solution.
Theorem 10.1
Suppose that x0 (t), y0 (t), and u0 (t) are continuously differentiable functions
of t in an interval I, and that a, b, and c are functions of x, y, and u with
continuous first order partial derivatives with respect to their argument in
some domain D of (x, y, u)space containing the initial curve
: x = x0 (t), y = y0 (t), u = u0 (t)
where t I. If (x0 (t), y0 (t), u0 (t)) is a point on that satisfies the condition
a(x0 (t), y0 (t), u0 (t))
dx0
dy0
(t) b(x0 (t), y0 (t), u0 (t))
(t) 6= 0
dt
dt
(10.3)
Figure 10.2
It follows that the Cauchy problem has a unique solution if C is nowhere
characteristic.
The unique solution is found as follows: We solve the PDE by the method
of characteristics to obtain the general solution that involves an unknown
function. This unknown function is determined using the initial condition.
We illustrate this process in the next example.
Example 10.1
Solve the Cauchy problem
ux + uy =1
u(x, 0) =f (x).
78
Solution.
The initial curve in R3 can be given parametrically as
: x0 (t) = t, y0 (t) = 0, u0 (t) = f (t).
We have
dy0
dx0
(t) b(x0 (t), y0 (t), u0 (t))
(t) = 1 6= 0
dt
dt
so by the above theorem the given Cauchy problem has a unique solution.
Next we apply the results of the previous section to find the unique solution.
If we solve the characteristic equations in non-parametric form
a(x0 (t), y0 (t), u0 (t))
dy
du
dx
=
=
1
1
1
we find x y = c1 and u x = c2 . Thus, the general solution of the PDE
is given by u = x + F (x y). Using the Cauchy data u(x, 0) = f (x) we find
f (x) = x + F (x) which implies that F (x) = f (x) x. Hence, the unique
solution is given by
u(x, y) = x + f (x y) (x y) = y + f (x y)
Next, if condition (10.3) is not satisfied and is a characteristic curve, i.e.,
dy0
dx0
=b(x0 , y0 , u0 )
dt
dt
du0
dx0
a(x0 , y0 , u0 )
=c(x0 , y0 , u0 )
dt
dt
dy0
du0
c(x0 , y0 , u0 )
=b(x0 , y0 , u0 )
dt
dt
for all point on then the problem has infinitely many solutions. To see
this, pick an arbitrary point P0 = (x0 , y0 , u0 ) on . Pick a new initial curve
0 passing through P0 which is nowhere characteristic in a neighborhood of
P0 . In this case, condition (10.3) is satisfied and the new Cauchy problem has
a unique solution. Since there are infinitely many ways of selecting 0 , we
obtain infinitely many solutions. We illustrate this case in the next example.
a(x0 , y0 , u0 )
Example 10.2
Solve the Cauchy problem
ux + uy =1
u(x, x) =x.
du0
dx0
6= a(x0 , y0 , u0 )
dt
dt
for some points on then the Cauchy problem has no solutions. We illustrate
this case next.
Example 10.3
Solve the Cauchy problem
ux + uy =1
u(x, x) =1.
Solution.
The initial curve in R3 can be given parametrically as
: x0 (t) = t, y0 (t) = t, u0 (t) = 1.
80
We have
a(x0 (t), y0 (t), u0 (t))
dy0
dx0
(t) b(x0 (t), y0 (t), u0 (t))
(t) = 0.
dt
dt
and
du0
dx0
6= a(x0 , y0 , u0 )
= 0.
dt
dt
As in Example 10.1, the general solution to the PDE is given by u = y +
f (x y). Using the Cauchy data u(x, x) = 1 we find f (0) = 1 x, which is
not possible since the LHS is a fixed number whereas the RHS is a variable
expression. Hence, the problem has no solutions
1 = c(x0 , y0 , u0 )
Example 10.4
Solve the Cauchy problem
ux uy =1
u(x, 0) =x2 .
Solution.
The initial curve is given parametrically by
: x0 (t) = t, y0 (t) = 0, u0 (t) = t2 .
We have
a(x0 (t), y0 (t), u0 (t))
dx0
dy0
(t) b(x0 (t), y0 (t), u0 (t))
(t) = 1 6= 0
dt
dt
dx0
dy0
(t) b(x0 (t), y0 (t), u0 (t))
(t) = 1 6= 0
dt
dt
dt
d(x + u)
=
1
x+u
t
we find that (x + u)e = c1 . Now, using the last two fractions we find
u2 x2 = k2 = f ((x + u)et ).
Using the Cauchy data u(x, 0) = 1, we find 1 x2 = f (1 + x) or f (1 + x) =
(1 + x)2 2x(1 + x). Thus, f (x) = x2 2x(x 1). The unique solution is
given by
u2 x2 = (x + u)2 e2t 2(x + u)et [(x + u)et 1]
or
u x = (x + u)e2t 2et [(x + u)et 1] = 2et (x + u)e2t .
This can be reduced further as follows: u + ue2t = x + 2et xe2t =
1e2t
2et
2et + x(1 e2t ) = u = 1+e
2t + x 1+e2t = sech(t) + xtanh(t)
Example 10.6
Solve the initial value problem
uux + uy = 1
where u(x, y) = 0 on the curve y 2 = 2x.
82
Solution.
A parametrization of is
: x0 (t) = 2t2 , y0 (t) = 2t, u0 (t) = 0, t > 0.
We have
a(x0 (t), y0 (t), u0 (t))
dy0
dx0
(t) b(x0 (t), y0 (t), u0 (t))
(t) = 4t 6= 0, t > 0
dt
dt
2u = y (4x y 2 ) 2 .
The solution surface satisfying u = 0 on y 2 = 2x is given by
1
2u = y (4x y 2 ) 2 .
This represents a solution surface only when y 2 < 4x. The solution does not
exist for y 2 > 4x
Practice Problems
Problem 10.1
Solve
(y u)ux + (u x)uy = x y
with the condition u x, x1 = 0, x 6= 0, 1.
Problem 10.2
Solve the linear equation
yux + xuy = u
with the Cauchy data u(x, 0) = x3 , x > 0.
Problem 10.3
Solve
x(y 2 + u)ux y(x2 + u)uy = (x2 y 2 )u
with the Cauchy data u(x, x) = 1, x > 0
Problem 10.4
Solve
xux + yuy = xeu
with the Cauchy data u(x, x2 ) = 0, x > 0.
Problem 10.5
Solve the initial value problem
xux + uy = 0, u(x, 0) = f (x).
Problem 10.6
Solve the initial value problem
ut + aux = 0, u(x, 0) = f (x).
Problem 10.7
Solve the initial value problem
aux + uy = u2 , u(x, 0) = cos x
84
Problem 10.8
Solve the initial value problem
ux + xuy = u, u(1, y) = h(y).
Problem 10.9
Solve the initial value problem
uux + uy = 0, u(x, 0) = f (x)
where f is an invertible function.
Problem 10.10
Solve the initial value problem
1 x2 ux + uy = 0, u(0, y) = y.
Problem 10.11
Consider
xux + 2yuy = 0.
(i) Find and sketch the characteristics.
(ii) Find the solution with u(1, y) = ey .
(iii) What happens if you try to find the solution satisfying either u(0, y) =
g(y) or u(x, 0) = h(x) for given functions g and h?
Problem 10.12
Solve the equation ux + uy = u subject to the condition u(x, 0) = cos x.
Problem 10.13
(a) Find the general solution of the equation
ux + yuy = u.
(b) Find the solution satisfying the Cauchy data u(x, 3ex ) = 2.
(c) Find the solution satisfying the Cauchy data u(x, ex ) = ex .
Problem 10.14
Solve the Cauchy problem
ux + 4uy = x(u + 1), u(x, 5x) = 1.
.
4
Problem 10.16
(a) Find the characteristics of the equation
yux + xuy = 0.
(b) Sketch some of the characteristics.
2
(c) Find the solution satisfying the boundary condition u(0, y) = ey .
Problem 10.17
Consider the equation ux + yuy = 0. Is there a solution satisfying the extra
condition
(a) u(x, 0) = 1
(b) u(x, 0) = x?
If yes, give a formula; if no, explain why.
86
87
(11.1)
(11.3)
89
Example 11.1
Determine whether the equation uxx + xuyy = 0 is hyperbolic, parabolic or
elliptic.
Solution.
Here we are given A = 1, B = 0, and C = x. Since B 2 4AC = 4x, the
given equation is hyperbolic if x < 0, parabolic if x = 0 and elliptic if x > 0
Second order partial differential equations arise in many areas of scientific
applications. In what follows we list some of the well-known models that are
of great interest:
1. The heat equation in one-dimensional space is given by
ut = kuxx
where k is a constant.
2. The wave equation in one-dimensional space is given by
utt = c2 uxx
where c is a constant.
3. The Laplace equation is given by
u = uxx + uyy = 0.
Practice Problems
Problem 11.1
Classify each of the following equation as hyperbolic, parabolic, or elliptic:
(a) Wave propagation: utt = c2 uxx , c > 0.
(b) Heat conduction: ut = cuxx , c > 0.
(c) Laplaces equation: u = uxx + uyy = 0.
Problem 11.2
Classify the following linear scalar PDE with constant coefficents as hyperbolic, parabolic or elliptic.
(a) uxx + 4uxy + 5uyy + ux + 2uy = 0.
(b) uxx 4uxy + 4uyy + 3ux + 4u = 0.
(c) uxx + 2uxy 3uyy + 2ux + 6uy = 0.
Problem 11.3
Find the region(s) in the xyplane where the equation
(1 + x)uxx + 2xyuxy y 2 uyy = 0
is elliptic, hyperbolic, or parabolic. Sketch these regions.
Problem 11.4
Show that u(x, t) = cos x sin t is a solution to the problem
utt
u(x, 0)
ut (x, 0)
ux (0, t)
=
=
=
=
uxx
0
cos x
0
91
Problem 11.6
Show that, for all (x, y) 6= (0, 0), u(x, y) = ln (x2 + y 2 ) is a solution of
uxx + uyy = 0,
and that, for all (x, y, z) 6= (0, 0, 0), u(x, y, z) =
1
x2 +y 2 +z 2
is a solution of
93
Figure 12.1
p
2u
(x)2 + (u)2 2 = T sin (x + x, t) T sin (x, t).
t
u u
=T (sin (x, t))
1+
2
x t
x
=T cos (x, t)
.
x
From Figure 12.1, we have
tan [(x, t)] = slope of tangent line to the string =
Applying the chain rule on this last equation, we find
2u
= cos2 (x, t) 2 .
x
x
Hence,
s
Since 0, we have
1+
u
x
2
2u
2u
3
=
T
cos
(x,
t)
.
t2
x2
u
.
x
95
Hence,
2u
2u
=
T
t2
x2
or
utt (x, t) = c2 uxx (x, t)
(12.1)
(12.2)
The function f (x ct) repesents wave that are moving to the right at a constant speed c and the function g(x + ct) represents waves that are moving to
w(s)ds.
0
Therefore,
1
1
f (x) = (v(x) W (x))
2
c
and
1
1
g(x) = (v(x) + W (x)).
2
c
Finally,
1
u(x, t) = [v(x ct) + v(x + ct) +
2
1
= [v(x ct) + v(x + ct) +
2
1
(W (x + ct) W (x ct))]
cZ
1 x+ct
w(s)ds]
c xct
97
Practice Problems
Problem 12.1
Show that if v(x, t) and w(x, t) satisfy equation (12.1) then v + w is also
a solution to (12.1), where and are constants.
Problem 12.2
Show that any linear time independent function u(x, t) = ax + b is a solution
to equation (12.1).
Problem 12.3
Find a solution to (12.1) that satisfies the homogeneous conditions u(x, 0) =
u(0, t) = u(L, t) = 0.
Problem 12.4
Solve the initial value problem
utt =9uxx
u(x, 0) = cos x
ut (x, 0) =0.
Problem 12.5
Solve the initial value problem
utt =uxx
1
u(x, 0) =
1 + x2
ut (x, 0) =0.
Problem 12.6
Solve the initial value problem
utt =4uxx
u(x, 0) =1
ut (x, 0) = cos (2x).
1 if x < 0
0 if x 0.
Problem 12.8
Solve the initial value problem
utt =c2 uxx
2
u(x, 0) =ex
ut (x, 0) = cos2 x.
Problem 12.9
Prove that the wave equation, utt = c2 uxx satisfies the following properties,
which are known as invariance properties. If u(x, t) is a solution, then
(i) Any translate, u(x y, t) where y is a fixed constant, is also a solution.
(ii) Any derivative, say ux (x, t), is also a solution.
(iii) Any dilation, u(ax, at), is a solution, for any fixed constant a.
Problem 12.10
Find v(r) if u(r, t) =
v(r)
r
Problem 12.11
Find the solution of the wave equation on the real line ( < x < +)
with the initial conditions
u(x, 0) = ex , ut (x, 0) = sin x.
99
Problem 12.12
The total energy of the string (the sum of the kinetic and potential energies)
is defined as
Z
1 L 2
(ut + c2 u2x )dx.
E(t) =
2 0
(a) Using the wave equation derive the equation of conservation of energy
dE(t)
= c2 (ut (L, t)ux (L, t) ut (0, t)ux (0, t)).
dt
(b) Assuming fixed ends boundary conditions, that is the ends of the string
are fixed so that u(0, t) = u(L, t) = 0, for all t > 0, show that the energy is
constant.
(c) Assuming free ends boundary conditions for both x = 0 and x = L, that
is both u(0, t) and u(L, t) vary with t, show that the energy is constant.
Problem 12.13
For a wave equation with damping
utt c2 uxx + dut = 0, d > 0, 0 < x < L
with the fixed ends boundary conditions show that the total energy decreases.
Problem 12.14
(a) Verify that for any twice differentiable R(x) the function
u(x, t) = R(x ct)
is a solution of the wave equation utt = c2 uxx . Such solutions are called
traveling waves.
(b) Show that the potential and kinetic energies (see Exercise 12.12) are
equal for the traveling wave solution in (a).
Problem 12.15
Find the solution of the Cauchy wave equation
utt = 4uxx
u(x, 0) = x2 , ut (x, 0) = sin 2x.
Simplify your answer as much as possible.
Figure 13.1
Assume the heat flows only in the xdirection, with the lateral sides well
insulated, and the only way heat can enter or leave the rod is at either end.
Also we assume that the temperature of the rod is constant at any point of
the cross section. In other words, temperature will only vary in x and we
can hence consider the rod to be a one spatial dimensional rod. We will also
assume that heat energy in any piece of the rod is conserved. That is, the
Figure 13.2
Divide the interval [x, x+x] into n subintervals each of length s using the
partition points x = s0 < s1 < < sn = x + x. Consider the portion Ui of
U of height s. The portion Ui is assumed to be thin so that the temperature
is constant throughout the volume. From the theory of heat conduction, the
quantity of heat Qi in Ui at time t is given by
Qi = cmi u(si1 , t) = cu(si1 , t)Vi
where mi is the mass of Ui , Vi is the volume of Ui and c is the specific
heat, that is, the amount of heat that it takes to raise one unit of mass of
the material by one unit of temperature.
But Ui is a cylinder of height s and area of base A so that Vi = As.
Hence,
Qi = cAu(si1 , t)s.
The quantity of heat in the portion U is given by
Q(t) = lim
n
X
i=1
Qi = lim
n
X
i=1
Z
cAu(si1 , t)s =
x+x
cAu(s, t)ds.
x
x+x
(13.1)
K
is called the diffusivity constant.
where k = c
Equation (13.1) is the one dimensional heat equation which is second order,
The total internal energy in the rod generated by the rods temeprature.
RL
0
u(x, 0)dx =
RL
0
f (x)dx + q0 Lt
E(t) =
0
Practice Problems
Problem 13.1
Show that if u(x, t) and v(x, t) satisfy equation (13.1) then u + v is also a
solution to (13.1), where and are constants.
Problem 13.2
Show that any linear time independent function u(x, t) = ax + b is a solution
to equation (13.1).
Problem 13.3
Find a linear time independent solution u to (13.1) that satisfies u(0, t) = T0
and u(L, T ) = TL .
Problem 13.4
Show that to solve (13.1) with the boundary conditions u(0, t) = T0 and
u(L, t) = TL it suffices to solve (13.1) with the homogeneous boundary
conditions u(0, t) = u(L, t) = 0.
Problem 13.5
Find a solution to (13.1) that satisfies the conditions u(x, 0) = u(0, t) =
u(L, t) = 0.
Problem 13.6
Let (I) denote equation (13.1) together with intial condition u(x, 0) = f (x),
where f is not the zero function, and the homogeneous boundary conditions
u(0, t) = u(L, t) = 0. Suppose a nontrivial solution to (I) can be written in
the form u(x, t) = X(x)T (t). Show that X and T satisfy the ODE
X 00 k X = 0 and T 0 T = 0
for some constant .
Problem 13.7
Consider again the solution u(x, t) = X(x)T (t). Clearly, T (t) = T (0)et .
Suppose that > 0.
(a) Show that X(x) = Aex + Bex , where = k and A and B are
arbitrary constants.
(b) Show that A and B satisfy the two equations A + B = 0 and A(eL
eL ) = 0.
(c) Show that A = 0 leads to a contradiction.
(d) Using (b) and (c) show that eL = eL . Show that this equality leads
to a contradiction. We conclude that < 0.
.
k
2 2
i
x
L
satis-
Problem 13.10
Suppose that a wire is stretched between 0 and a. Describe the boundary
conditions for the temperature u(x, t) when
(i) the left end is kept at 0 degrees and the right end is kept at 100 degrees;
and
(ii) when both ends are insulated.
Problem 13.11
Let ut = uxx for 0 < x < and t > 0 with boundary conditions
R 2 u(0,2t) =
0 = u(, t) and initial condition u(x, 0) = sin x. Let E(t) = 0 (ut + ux )dx.
Show that E 0 (t) < 0.
Problem 13.12
Suppose
ut = uxx + 4, ux (0, t) = 5, ux (L, t) = 6, u(x, 0) = f (x).
Calculate the total thermal energy of the one-dimensional rod (as a function
of time).
Problem 13.13
Consider the heat equation
ut = kuxx
for x (0, 1) and t > 0, with boundary conditions u(0, t) = 2 and u(1, t) = 3
for t > 0 and initial condition u(x, 0) = x for x (0, 1). A function v(x) that
satisfies the equation v 00 (x) = 0, with conditions v(0) = 2 and v(1) = 3 is
called a steady-state solution. That is, the steady-state solutions of the
heat equation are those solutions that dont depend on time. Find v(x).
Example 14.2
For each positive integer n let fn : (0, ) (0, ) be given by fn (x) = nx.
Show that {fn }
n=1 does not converge pointwise in D = (0, ).
One of the weaknesses of this type of convergence is that it does not preserve
some of the properties of the base functions {fn }
n=1 . For example, if each fn
is continuous then the pointwise limit function need not be continuous. (See
Problem 14.1) A stronger type of convergence which preserves most of the
properties of the base functions is the uniform convergence which we define
next.
Let D be a subset of R and let {fn }
n=1 be a sequence of functions defined on
Figure 14.1
Example 14.3
For each positive integer n let fn : [0, 1] R be given by fn (x) = nx . Show
that {fn }
n=1 converges uniformly to the zero function.
Solution.
Let > 0 be given. Let N be a positive integer such that N > 1 . Then for
f
n
n 2
Problem 14.1 shows a sequence of continuous functions converging pointwise
to a discontinuous function. That is, pointwise convergence does not preserve the property of continuity. One of the interesting features of uniform
convergence is that it preserves continuity as shown in the next example.
Example 14.5
Suppose that for each n 1 the function fn : D R is continuous in D.
Suppose that {fn }
n=1 converges uniformly to f. Let a D.
(a) Let > 0 be given. Show that there is a positive integer N such that if
n N then |fn (x) f (x)| < 3 for all x D.
(b) Show that there is a > 0 such that for all |x a| < we have |fN (x)
fN (a)| < 3 .
(c) Using (a) and (b) show that for |x a| < we have |f (x) f (a)| < .
Hence, f is continuous in D since a was arbitrary.
xa n
n xa
Indeed,
lim lim fn (x) = lim f (x)
xa n
xa
Solution.
We have
Z
1
x
x2
dx =
= .
n
2n 1
Hence,
Z
fn (x)dx =
lim
whereas
lim fn (x)dx = 0
x
n
converges
lim
fn (x)dx =
a
lim fn (x)dx.
a n
Theorem 14.1
Suppose that fn : [a, b] R is a sequence of continuous functions that
converges uniformly to f : [a, b] R. Then
Z b
Z b
Z b
lim
fn (x)dx =
lim fn (x)dx =
f (x)dx.
n
a n
Proof.
From Example 14.5, we have that f is continuous and hence integrable. Let
> 0 be given. By uniform convergence, we can find a positive integer N
for all x in [a, b] and n N. Thus, for n N ,
such that |fn (x) f (x) < ba
we have
Z b
Z b
Z b
|fn (x) f (x)|dx < .
f (x)dx
fn (x)dx
a
x2 + n1 .
1
|fn (x) f (x)|
n
q
n
x2 + n1 + x2
1
n
=q
x2 + n1 + x2
1
n
1
n
1
= .
n
(c) Let > 0 be given. Since limn 1n = 0 we can find a positive integer
N such that for all n N we have 1n < . Now the answer to the question
follows from this and part (b)
Even when uniform convergence occurs, the process of interchanging limits and differentiation may fail as shown in the next example.
Example 14.8
Consider the functions fn : R R defined by fn (x) = sinnnx .
(a) Show that {fn }
n=1 converges uniformly to the function f (x) = 0.
(b) Note that {fn }
n=1 and f are differentiable functions. Show that
h
i0
lim fn0 (x) 6= f 0 (x) = lim fn (x) .
Proof.
First, note that the function g is continuous in [a, b] since uniform convergence
preserves continuity. Let c be an arbitrary point in [a, b]. Then
Z x
fn0 (t)dt = fn (x) fn (c), x [a, b].
c
Taking the limit of both sides and using the facts that fn0 converges uniformly
to g and fn converges pointwise to f , we can write
Z x
g(t)dt = f (x) f (c).
c
Proof.
Suppose that fn converges uniformly to f. Let > 0 be given. Then there
is a positive integer N such that |fn (x) f (x)| < 2 for all n N and all
x D. Thus, for n N, we have
sup{|fn (x) f (x)| : x D}
< .
2
cos x
n
cos x
1
| : x R} .
n
n
8
Now apply the squeeze rule for sequences we find that
0 sup{|
lim sup{|
cos x
| : x R} = 0
n
which implies that the given sequence converges uniformly to the zero function on R
L.
Practice Problems
Problem 14.1
Define fn : [0, 1] R by fn (x) = xn . Define f : [0, 1] R by
0 if 0 x < 1
f (x) =
1
if x = 1.
(a) Show that the sequence {fn }
n=1 converges pointwise to f.
(b) Show that the sequence {fn }
n=1 does not converge uniformly to f. Hint:
Suppose otherwise. Let = 0.5 and get a contradiction by using a point
1
(0.5) N < x < 1.
Problem 14.2
Consider the sequence of functions
fn (x) =
nx + x2
n2
sin (nx + 3)
n+1
n+cos x
.
2n+sin2 x
Problem 14.9
Show that the sequence defined by fn (x) = (cos x)n does not converge uniformly on [ 2 , 2 ].
Problem 14.10
Let {fn }
n=1 be a sequence of functions such that
sup{|fn (x)| : 2 x 5}
2n
.
1 + 4n
n=1
Example 15.1
P
n
Show that
n=0 x converges pointwise to a function to be determined for
all all 1 < x < 1.
Solution.
The nth term of the sequence of partial sums is given by
Sn (x) = 1 + x + x2 + + xn =
1 xn+1
.
1x
Since
lim xn+1 = 0 1 < x < 1,
X
n=0
xn =
121
1
.
1x
Thus,
1
1x
P
Likewise, we say that a series of functions
n=1 fn (x) converges uniformly
to a function f if and only if the sequence of partial sums {Sn }
n=1 converges
uniformly to f.
The following theorem provide a tool for uniform convergence.
Theorem 15.1 (Weierstrass M-test)
P
Suppose that for each x in an interval I the series
n=1 fn (x) is well-defined.
Suppose further that
|fn (x)| Mn , x I.
P
P
If n=1 Mn (a scalar series) is convergent then the series
n=1 fn (x) is uniformly convergent.
Example 15.2
P
Show that
n=1
sin (nx)
n2
is uniformly convergent.
Solution.
For all x R, we have
sin (nx) | sin (nx)|
1
2.
n2
2
n
n
P
1
The series
n=1 n2 is convergent being a pseries with p = 2 > 1. Thus, by
Weierstrass M-test the given series is uniformly convergent
In this section we introduce a type of series of functions known as Fourier
series. They are given by
n
n i
a0 X h
f (x) =
+
an cos
x + bn sin
x , LxL
2
L
L
n=1
(15.1)
where an and bn are called the Fourier coefficients. The expression on the
right is called a trigonometric series. Note that we begin the series with
a0
as opposed to simply a0 to simplify the coefficient formula for an that we
2
Figure 15.1
For a T periodic function we have
f (x) = f (x + T ) = f (x + 2T ) = .
Note that the definite integral of a T periodic function is the same over any
interval of length T. By Problem 15.1 below, if f and g are two periodic functions with common period T, then the product f g and an arbitrary linear
combination c1 f + c2 g are also periodic with period T. It is an easy exercise
to show that the Fourier series (15.1) is periodic with fundamental period 2L.
Orthogonality Property
Recall from Calculus that for each pair of vectors ~u and ~v we associate a
scalar quantity ~u ~v called the dot product of ~u and ~v . We say that ~u and ~v
are orthogonal if and only if ~u ~v = 0. We want to define a similar concept
for functions.
Let f and g be two functions with domain the closed interval [a, b]. We define
123
We call < f, g > the inner product of f and g. We say that f and g are
orthogonal if and only if < f, g >= 0. A set of functions is said to be mutually orthogonal if each distinct pair of functions in the set is orthogonal.
Before we proceed any further into computations, we would like to remind
the reader of the following two facts from calculus:
RL
If f (x) is an odd function defined on the interval [L, L] then L f (x)dx =
0.
RL
If f (x) is an even function defined on the interval [L, L] then L f (x)dx =
RL
2 0 f (x)dx.
Example 15.3
Show that the set 1, cos n
x , sin
L
ally orthogonal in [L, L].
n
x
L
: n N , where m 6= n, is mutu-
Solution.
Since the cosine function is even, we have
Z L
Z L
n
n
2L h n iL
x dx = 2
cos
x dx =
sin
x
= 0.
1 cos
L
L
n
L
0
0
L
Since the sine function is odd, we have
Z L
n
1 sin
x dx = 0.
L
L
Now, for n 6= m we have
Z L
Z
m
n
1 L
(m + n)
(m n)
cos
x cos
x dx =
cos
x + cos
x dx
L
L
2 L
L
L
L
1
L
(m + n)
=
sin
x
2 (m + n)
L
L
L
(m n)
+
sin
x
=0
(m n)
L
L
sin
L
and
Z
m
n
x sin
x dx = 0
L
L
n
m
x sin
x dx = 0
cos
L
L
L
The reason we care about these functions being orthogonal is because we will
exploit this fact to develop a formula for the coefficients in our Fourier series.
Now, in order to answer the first question mentioned earlier, that is, which
functions can be expressed as a Fourier series expansion, we need to introduce some mathematical concepts.
A function f (x) is said to be piecewise continuous on [a, b] if it is continuous in [a, b] except possibly at finitely many points of discontinuity within
the interval [a, b], and at each point of discontinuity, the right- and lefthanded limits of f exist. An example of a piecewise continuous function is
the function
x
0x<1
f (x) =
2
x x 1 x 2.
We will say that f is piecewise smooth in [a, b] if and only if f (x) as well
as its derivatives are piecewise continuous.
The following theorem, proven in more advanced books, ensures that a
Fourier decomposition can be found for any function which is piecewise
smooth.
Theorem 15.2
Let f be a 2L-periodic function. If f is a piecewise smooth on [L, L] then
for all points of discontinuity x [L, L] we have
n
n i
a0 X h
f (x ) + f (x+ )
=
+
an cos
x + bn sin
x
2
2
L
L
n=1
125
n
n i
a0 X h
f (x) =
+
an cos
x + bn sin
x .
2
L
L
n=1
Remark 15.1
(1) Almost all functions occurring in practice are piecewise smooth functions.
(2) Given a piecewise smooth function f on [L, L]. The above theorem
applies to the periodic extension F of f where F (x + 2nL) = f (x) (n Z)
and F (x) = f (x) on (L, L). Note that if f (L) = f (L) then F (x) = f (x)
on [L, L]. Otherwise, the end points of f (x) may be jump discontinuities
of F (x).
Convergence Results of Fourier Series
We list few of the results regarding the convergence of Fourier series:
(1) The type of convergence in the above theorem is pointwise convergence.
(2) The convergence is uniform for a continuous function f on [L, L] such
that f (L) = f (L).
P
2
2
(3) The convergence is uniform whenever
n=1 (|an | + |bn | ) is convergent.
(4) If f (x) is periodic, continuous, and has a piecewise continuous derivative,
then the Fourier Series corresponding to f converges uniformly to f (x) for
the entire real line.
(5) The convergence is uniform on any closed interval that does not contain
a point of discontinuity.
Euler-Fourier Formulas
Next, we will answer the second question mentioned earlier, that is, the question of finding formulas for the coefficients an and bn . These formulas for an
and bn are called Euler-Fourier formulas which we derive next. We will assume that the RHS in (15.1) converges uniformly to f (x) on the interval
[L, L]. Integrating both sides of (15.1) we obtain
Z LX
Z L
Z L
h
n
n i
a0
dx +
an cos
x + bn sin
x dx.
f (x)dx =
L
L
L n=1
L
L 2
Since the trigonometric series is assumed to be uniformly convergent, from
Theorem 14.2, we can interchange the order of integration and summation
to obtain
Z L
Z L
Z L h
n
n i
X
a0
f (x)dx =
dx +
an cos
x + bn sin
x dx.
L
L
L
L 2
n=1 L
n
n iL
L
x dx =
sin
x
=0
L
n
L
L
cos
L
and likewise
Z
sin
L
n iL
n
L
x dx =
cos
x
= 0.
L
n
L
L
Thus,
Z
1 L
f (x)dx.
a0 =
L L
To find the other Fourier coefficients, we recall the results of Problems 15.2
- 15.3 below.
Z L
n
m
L if m = n
cos
x cos
x dx =
0 if m 6= n
L
L
L
Z L
n
m
L if m = n
sin
x sin
x dx =
0 if m 6= n
L
L
L
Z L
n
m
sin
x cos
x dx = 0, m, n.
L
L
L
Now, to find the formula for the Fourier
coefficients am for m > 0, we multiply
both sides of (15.1) by cos m
x
and
integrate
from L to L to otbain
L
Z L
Z L
m Z L a
m
n
m
X
0
f (x) cos
x =
cos
x dx +
an
cos
x cos
x dx
L
L
L
L
L 2
L
L
n=1
Z L
n
m
+ bn
sin
x cos
x dx .
L
L
L
Hence,
Z
f (x) cos
L
and therefore
m
x dx = am L
L
1
am =
L
Likewise, we can show that
1
bm =
L
f (x) cos
m
x dx.
L
f (x) sin
m
x dx.
L
1
f (x)dx =
xdx =
a0 =
0
2
Z
1
1 x sin nx cos nx
(1)n 1
an =
x cos nxdx =
+
=
0
n
n2 0
n2
Z
1
x cos nx sin nx
(1)n+1
1
x sin nxdx =
+
=
.
bn =
0
n
n2 0
n
Hence,
X (1)n 1
(1)n+1
f (x) = +
cos (nx) +
sin (nx) < x <
4 n=1
n2
n
Example 15.5
Apply Theorem 15.2 to the function in Example 15.4.
Solution.
Let F be a periodic extension of f of period 2. See Figure 152.
Figure 152
127
= .
2
2
Similar result for x = . Clearly, F is a piecewise smooth function so that
by the previous thereom we can write
,
if x =
n
n+1
X
2
(1) 1
(1)
f
(x),
if
<x<
+
sin
(nx)
=
cos
(nx)
+
4 n=1
n2
n
,
if x = .
2
Taking x = we have the identity
X (1)n 1
+
(1)n =
2
4 n=1
n
2
which can be simplified to
X
n=1
1
2
=
.
(2n 1)2
8
129
Theorem 15.4
A Fourier series of a continuous function f (x) can be differentiated term by
term if f 0 (x) is piecewise smooth. The result of the differentiation is the
Fourier series of f 0 (x).
Practice Problems
Problem 15.1
Let f and g be two functions with common domain D and common period
T. Show that
(a) f g is periodic of period T.
(b) c1 f + c2 g is periodic of period T, where c1 and c2 are real numbers.
Problem 15.2
Show
have
R Lthat for m6= n we
n
x
sin
x
dx = 0 and
(a) L sin m
L
L
RL
(b) L cos m
x sin n
x dx = 0.
L
L
Problem 15.3
Compute
the following
integrals:
RL
2 n
(a) L cos L x dx.
RL
(b) L sin2 n
x dx.
L
RL
n
(c) L cos n
x
sin
x
dx.
L
L
Problem 15.4
Find the Fourier coefficients of
, x < 0
,
0<x<
f (x) =
0,
x = 0,
on the interval [, ].
Problem 15.5
Find the Fourier series of f (x) = x2
1
2
Problem 15.6
Find the Fourier series of the function
1, 2 < x <
0,
< x <
f (x) =
1,
< x < 2.
131
Problem 15.7
Find the Fourier series of the function
1 + x, 2 x 0
f (x) =
1 x, 0 < x 2.
Problem 15.8
Show that f (x) =
1
x
Problem 15.9
Assume that f (x) is continuous and has period 2L. Prove that
Z L+a
Z L
f (x)dx
f (x)dx =
L
L+a
a0 X
+
(an cos nt + bn sin nt), 0 < x < 2.
2
n=1
Find
a3
.
b3
Problem 15.13
Let f (x) = x3 on [, ], extended periodically to all of R. Find the Fourier
coefficients an , n = 1, 2, 3, .
Problem 15.14
Let f (x) be the square wave function
x < 0
f (x) =
0x
extended periodically to all of R. To what value does the Fourier series of
f (x) converge when x = 0?
Problem 15.15
(a) Find the Fourier series of
f (x) =
1 x < 0
2 0x
133
(b) Show that for any odd function f (x) defined on the interval [L, L], we
have
Z L
f (x)dx = 0.
L
Solution.
(a) Since f (x) is even, we have f (x) = f (x) for all x in [L, L]. Thus,
Z 0
Z 0
Z 0
Z L
f (x)dx =
f (x)dx =
f (u)du =
f (x)dx.
L
Z
f (x)dx +
Z
f (x)dx = 2
f (x)dx.
0
Hence,
Z
f (x)dx
f (x)dx =
f (x)dx +
0=
Example 16.3
RL
(a) Find the value of the integral L f (x) sin
RL
(b) Find the value of the integral L f (x) cos
n
x dx when f
L
n
x
dx when f
L
is even.
is odd.
Solution.
x is odd and f is even, we have that f (x) sin n
x
(a) Since the function sin n
L
L
is odd so that
Z L
n
f (x) sin
x dx = 0
L
L
n
x
is
even
and
f
is
odd,
we
have
that
f
(x)
cos
x
(b) Since the function cos n
L
L
is odd so that
Z L
n
f (x) cos
x dx = 0
L
L
Even and Odd Extensions
Let f : (0, L) R be a piecewise smooth function. We define the odd
extension of this function on the interval L < x < L by
f (x)
0<x<L
fodd (x) =
f (x) L < x < 0.
This function will be odd on the interval (L, L), and will be equal to f (x)
on the interval (0, L). We can then further extend this function to the entire
real line by defining it to be 2L periodic. Let f odd denote this extension. We
note that f odd is an odd function and piecewise smooth so that by Theorem
15.2 it possesses a Fourier series expansion, and from the fact that it is odd
all of the a0n s are zero. Moreover, in the interval (0, L) we have
n
f (x) =
bn sin
x .
L
n=1
(16.1)
135
f (x) =
n
a0 X
+
an cos
x .
2
L
n=1
(16.2)
We call (16.2) the Fourier cosine series of f. The coefficients an are given
by
Z
n
2 L
f (x) cos
x dx, n = 0, 1, 2, .
an =
L 0
L
Example 16.4
Graph the odd and even extensions of the function f (x) = x, 0 x 1.
Solution.
We have fodd (x) = x for 1 x 1. The odd extension of f is shown in
Figure 16.1(a). Likewise,
x
0x1
feven (x) =
x 1 x < 0.
Figure 16.1
Example 16.5
Find the Fourier sine series of the function
f (x) =
0 x 2
x .
2
x,
x,
Solution.
We have
2
bn =
"Z
( x) sin nxdx .
x sin nxdx +
0
Z
0
Z
h x
i 2
1 2
x sin nxdx = cos nx +
cos nxdx
n
n 0
0
cos (n/2)
1
=
+ 2 [sin nx]02
2n
n
cos (n/2) sin (n/2)
=
+
2n
n2
( x)
( x) sin nxdx =
cos nx
n
cos nxdx
cos (n/2)
1
2 [sin nx]
2
2n
n
cos (n/2) sin (n/2)
=
+
2n
n2
=
Thus,
4 sin (n/2)
,
n2
and the Fourier sine series of f (x) is
bn =
f (x) =
X
4 sin (n/2)
n=1
n2
X
4(1)n+1
sin nx =
sin (2n 1)x
(2n 1)2
n=1
137
Practice Problems
Problem 16.1
Give an example of a function that is both even and odd.
Problem 16.2
Graph the odd and even extensions of the function f (x) = 1, 0 x 1.
Problem 16.3
Graph the odd and even extensions of the function f (x) = L x for 0 x
L.
Problem 16.4
Graph the odd and even extensions of the function f (x) = 1 + x2 for 0
x L.
Problem 16.5
Find the Fourier cosine series of the function
x,
0 x 2
f (x) =
x, 2 x .
Problem 16.6
Find the Fourier cosine series of f (x) = x on the interval [0, ].
Problem 16.7
Find the Fourier sine series of f (x) = 1 on the interval [0, ].
Problem 16.8
Find the Fourier sine series of f (x) = cos x on the interval [0, ].
Problem 16.9
Find the Fourier cosine series of f (x) = e2x on the interval [0, 1].
Problem 16.10
For the following functions on the interval [0, L], find the coefficients bn of
the Fourier sine expansion.
(a) f (x) = sin 2
x
.
L
(b) f (x) = 1
(c) f (x) = cos L x .
139
Problem 16.11
For the following functions on the interval [0, L], find the coefficients an of
the Fourier cosine expansion.
(a) f (x) = 5 + cos L x .
(b) f (x) = x
(c)
1 0 < x L2
f (x) =
0 L2 < x L.
Problem 16.12
Consider a function f (x), defined on 0 x L, which is even (symmetric)
around x = L2 . Show that the even coefficients (n even) of the Fourier sine
series are zero.
Problem 16.13
Consider a function f (x), defined on 0 x L, which is odd around x = L2 .
Show that the even coefficients (n even) of the Fourier cosine series are zero.
Problem 16.14
The Fourier sine series of f (x) = cos
cos
x
L
bn sin
n=1
where
b1 = 0, bn =
x
L
(n2
for 0 x L is given by
nx
L
, nN
2n
[1 + (1)n ].
1)
x
L
.
Problem 16.15
Consider the function
f (x) =
1 0x<1
2 1 x < 2.
b b2 4ac
.
r1,2 =
2a
We consider the following three cases:
If b2 4ac > 0 then the general solution to (17.1) is given by
y(t) = Ae
b2 4ac
2a
t
b+
+ Be
b2 4ac
2a
t
y(t) = Ae 2a t + Bte 2a t .
If b2 4ac < 0 then
b
4ac b2
r1,2 = i
2a
2a
and the general solution to (17.1) is given by
b
b
4ac b2
4ac b2
2a
t
2a
t
y(t) = Ae
cos
t + Be
sin
t.
2a
2a
141
1 1 + 4
=
.
2
If > 14 then
X(x) = Ae
In this case,
u(x, t) = De
If = 41 then
1+ 1+4
x
2
1+ 1+4
x
2
+ Be
1 1+4
x
2
e + Ee
x
1 1+4
x
2
et .
X(x) = Ae 2 + Bxe 2
and in this case
u(x, t) = (D + Ex)e 2 4 .
If < 14 then
x
2
X(x) = Ae cos
!
p
(1 + 4)
x
x + Be 2 sin
2
!
p
(1 + 4)
x .
2
In this case,
u(x, t) = D0 e
x
+t
2
cos
!
p
(1 + 4)
x
x + B 0 e 2 +t sin
2
(1 + 4)
x
2
Example 17.2
Solve Laplaces equation using the separation of variables method
u = uxx + uyy = 0.
143
Solution.
We look for a solution of the form u(x, y) = X(x)Y (y). Substituting in the
Laplaces equation, we obtain
X 00 (x)Y (y) + X(x)Y 00 (y) = 0.
Assuming X(x)Y (y) is nonzero, dividing for X(x)Y (y) and subtracting
from both sides, we find:
Y 00 (y)
Y (y)
X 00 (x)
Y 00 (y)
=
.
X(x)
Y (y)
The left hand side is a function of x while the right hand side is a function
of y. This says that they must equal to a constant. That is,
Y 00 (y)
X 00 (x)
=
=
X(x)
Y (y)
where is a constant. This results in the following two ODEs
X 00 X = 0 and Y 00 + Y = 0.
The solutions of these equations depend on the sign of .
If > 0 then the solutions are given
+ Be x
X(x) =Ae
Y (y) =Ce
+ De
Y0
X0
=
.
xX
yY
The left hand side is a function of x while the right hand side is a function
of y. This says that they must equal to a constant. That is,
X0
Y0
=
=
xX
yY
where is a constant. This results in the following two ODEs
X 0 xX = 0 and Y 0 yY = 0.
Solving these equations using the method of separation of variable for ODEs
y 2
x2
we find X(x) = Ae 2 and Y (y) = Be 2 . Thus, the general solution is given
by
u(x, y) = Ce
(x2 +y 2 )
2
145
Practice Problems
Problem 17.1
Solve using the separation of variables method
u + u = 0.
Problem 17.2
Solve using the separation of variables method
ut = kuxx .
Problem 17.3
Derive the system of ordinary differential equations for R(r) and () that
is satisfied by solutions to
1
1
urr + ur + 2 u = 0.
r
r
Problem 17.4
Derive the system of ordinary differential equations and boundary conditions
for X(x) and T (t) that is satisfied by solutions to
utt = uxx 2u, 0 < x < 1, t > 0
u(0, t) = 0 = ux (1, t) t > 0
of the form u(x, t) = X(x)T (t). (Note: you do not need to solve for X and
T .)
Problem 17.5
Derive the system of ordinary differential equations and boundary conditions
for X(x) and T (t) that is satisfied by solutions to
ut = kuxx , 0 < x < L, t > 0
u(x, 0) = f (x), u(0, t) = 0 = ux (L, t) t > 0
of the form u(x, t) = X(x)T (t). (Note: you do not need to solve for X and
T .)
Problem 17.6
Find all product solutions of the PDE ux + ut = 0.
= and
T0
kT
= .
n
x, n N.
L
n2 2
kt
L2
, n N.
Cn sin
n=1
n n2 2
x e L2 kt .
L
(18.1)
n
f (x) =
Cn sin
x .
L
n=1
n n2 2
TL T0
Cn sin
x e L2 kt + T0 +
x.
u(x, t) =
L
L
n=1
The Heat Equation with Neumann Boundary Conditions
When both ends of the bar are insulated, that is, there is no heat flow out
of them, we use the boundary conditions
ux (0, t) = ux (L, t) = 0.
In this case, the general form of the heat equation initial boundary value
problem is to find u(x, t) satisfying
ut (x, t) =kuxx (x, t), 0 x L, t > 0
u(x, 0) =f (x), 0 x L
ux (0, t) =ux (L, t) = 0, t > 0.
Since 0 = ux (0, t) = X 0 (0)T (t) we obtain X 0 (0) = 0. Likewise, 0 = ux (L, t) =
X 0 (L)T (t) implies X 0 (L) = 0. We again consider the following three cases:
If = 0 then X(x) = A + Bx. Since X 0 (0) = 0, we find B = 0. Thus,
X(x) = A and T (t) = constant so that u(x, t) = constant which is impossible
if f (x) is not the constant function.
If > 0 then a simple calculation shows that u(x, t) is the trivial solution.
Again, because of the condition
this solution is discarded.
u(x, 0) = f (x),
If < 0 then X(x) = A cos x + B sin x and upon differentiation
with respect to x we find
n n2 2
x e L2 kt .
L
By the superposition principle, the required solution to the heat equation
with Neumann boundary conditions is given by
un (x, t) = Cn cos
u(x, t) =
where
2
Cn =
L
n n2 2
C0 X
Cn cos
+
x e L2 kt
2
L
n=1
f (x) cos
0
n
x dx, n = 0, 1, 2, .
L
Practice Problems
Problem 18.1
Find the temperature in a bar of length 2 whose ends are kept at zero
and lateral
surface insulated if the initial temperature is f (x) = sin 2 x +
x .
3 sin 5
2
Problem 18.2
Find the temperature in a homogeneous bar of heat conducting material of
length L with its end points kept at zero and initial temperature distribution
given by f (x) = Lxd2 (L x), 0 x L.
Problem 18.3
Find the temperature in a thin metal rod of length L, with both ends insulated (so that there is no passage of heat
through the ends) and with initial
sin
L
L
Problem 18.13
Using the method of separation of variables find the solution of the heat
equation
ut = kuxx
satisfying the following boundary and initial conditions:
+ 4 cos 5x
.
(a) ux (0, t) = ux (L, t) = 0, u(x, 0) = cos x
L
L
(b) ux (0, t) = ux (L, t) = 0, u(x, 0) = 5.
Problem 18.14
Find the solution of the following heat conduction partial differential equation
ut = 8uxx , 0 < x < 4, t > 0
u(0, t) = u(4, t) = 0, t > 0
u(x, 0) = 6 sin x, 0 < x < 4.
(19.1)
where the symbol is referred to as the Laplacian. Solutions of this equation are called harmonic functions.
Example 19.1
Show that, for all (x, y) 6= (0, 0), u(x, y) = ax2 ay 2 + cx + dy + e is a
harmonic function, where a, b, c, d, and e are constants.
Solution.
We have
ux
uxx
uy
uyy
=2ax + c
=2a
= 2ay + d
= 2a.
Plugging these expressions into the equation we find uxx + uyy = 0. Hence,
u(x, y) is harmonic
The Laplace equation is arguably the most important differential equation in
all of applied mathematics. It arises in an astonishing variety of mathematical and physical systems, ranging through fluid mechanics, electromagnetism,
potential theory, solid mechanics, heat conduction, geometry, probability,
uxx + uyy = 0
uxx + uyy = 0
u(0, y) = f1 (y),
u(a, y) = f2 (y),
(I)
(II)
u(a, y) = u(x, 0) = u(x, b) = 0
u(0, y) = u(x, 0) = u(x, b) = 0
uxx + uyy = 0
u(x, 0) = g1 (x),
(III)
u(0, y) = u(a, y) = u(x, b) = 0
(IV )
uxx + uyy = 0
u(x, b) = g2 (x),
u(0, y) = u(a, y) = u(x, 0) = 0.
X(x) = 0 and
this leads to
B sin a = 0. We must have B 6= 0 otherwise
the trivial solution. Since B 6= 0, we obtain sin a = 0 or a = n
2 2
where n Z. Solving for we find n = na2 . Thus, we obtain infinitely
many solutions given by
Xn (x) = sin
n
x, n N.
a
Yn (y) = an e
n y
+ bn e
n y
= An cosh
p
p
n y + Bn sinh n y, n N.
X
n=1
Bn sin
n
n
x sinh
y .
a
a
(19.2)
X
n=1
n
n
Bn sinh
b sin
x .
a
a
uxx + uyy = 0
u(0, y) = f1 (y),
Y 00 (y)
Y (y)
X 00 (x)
Y 00 (y)
=
.
X(x)
Y (y)
The left hand side is a function of x while the right hand side is a function
of y. This says that they must equal to a constant. That is,
X 00 (x)
Y 00 (y)
=
=
X(x)
Y (y)
where is a constant. This results in the following two ODEs
X 00 X = 0 and Y 00 + Y = 0.
As far as the boundary conditions, we have for all 0 x a and 0 y b
u(0, y) = f1 (y) = X(0)Y (y)
u(a, y) = 0 = X(a)Y (y) = X(a) = 0
Xn (x) = an e
n x
+ bn e
n x
= An cosh
n
n
x + Bn sinh
x , n N.
b
b
However, this is not really suited for dealing with the boundary condition
X(a) = 0. So, lets also notice that the following is also a solution.
n
n
Xn (x) = An cosh
(x a) + Bn sinh
(x a) , n N.
b
b
Now, using the boundary condition X(a) = 0 we obtain An = 0 for all n N.
Hence, the functions
n
n
un (x, y) = Bn sin
y sinh
(x a) , n N
b
b
Bn sin
n=1
n
n
y sinh
(x a) .
b
b
(19.4)
X
n=1
n
n
Bn sinh a sin
y .
b
b
Since the right-hand side is the Fourier sine series of f1 on the interval [0, b],
the coefficients Bn are given by
Z b
n i1
n h
2
f1 (y) sin
y dy sinh a
.
(19.5)
Bn =
b 0
b
b
Thus, the solution to the Laplaces equation is given by (19.4) with the Bn0 s
calculated from (19.5)
Example 19.3
Solve
uxx + uyy = 0, 0 < x < L, 0 < y < H
u(0, y) = u(L, y) = 0, 0 < y < H
u(x, 0) = uy (x, 0), u(x, H) = f (x), 0 < x < L.
Solution.
Using separation of variables we find
Y 00
X 00
=
= .
X
Y
We first solve
X 00 X = 0
0<x<L
X(0) = X(L) = 0.
2 2
n
x, n N.
L
Y 00 + Y = 0
0<y<H
0
Y (0) Y (0) = 0.
n
n i
X
n h n
x
cosh
y + sinh
y .
u(x, y) =
Bn sin
L
L
L
L
n=1
Yn (y) = An cosh
n
n i
X
n h n
f (x) =
Bn sin
x
cosh
H + sinh
H .
L
L
L
L
n=1
Recall the Fourier sine series of f on [0, L] given by
f (x) =
X
n=1
An sin
n
x
L
where
Z
n
2 L
f (x) sin
x dx.
An =
L 0
L
Thus, the general solution is given by
n
n i
X
n h n
u(x, y) =
Bn sin
x
cosh
y + sinh
y .
L
L
L
L
n=1
with the Bn satisfying
Z
n
n i
n
h n
2 L
Bn
cosh
H + sinh
H =
f (x) sin
x dx
L
L
L
L 0
L
Practice Problems
Problem 19.1
Solve
uxx + uyy = 0
u(a, y) = f2 (y),
uxx + uyy = 0
u(x, 0) = g1 (x),
uxx + uyy = 0
u(x, 0) = u(0, y) = 0,
H
H
y
2
2
subject to
u(0, y) = u(L, y) = 0,
u(x,
H
H
<y<
2
2
H
H
) = f1 (x), u(x, ) = f2 (x), 0 x L.
2
2
Problem 19.6
an z n .
n=0
x
2
4
y , u(2, y) = 7.
3
Problem 19.9
Solve
uxx + uyy = 0, 0 x L, 0 y H
subject to
uy (x, 0) = 0, u(x, H) = 0
y
u(0, y) = u(L, y) = 4 cos
.
2H
Problem 19.10
Solve
uxx + uyy = 0, x > 0, 0 y H
subject to
u(0, y) = f (y), |u(x, 0)| <
uy (x, 0) = uy (x, H) = 0.
x
L
5 sin
3x
.
L
165
sin
u
r
uxx =uxr rx + ux x
sin
sin
= cos urr + 2 u
ur cos
r
r
sin
sin
cos
u
u
+ sin ur + cos ur
r
r
r
cos
uy =ur ry + u y = sin ur +
u
r
uyy =uyr ry + uy y
cos
cos
= sin urr 2 u +
ur sin
r
r
cos
cos
sin
u +
u
.
+ cos ur + sin ur
r
r
r
Substituting these equations into u = 0 we obtain
1
1
urr + ur + 2 u = 0.
r
r
Example 20.1
Find the solution to
u = 0, x2 + y 2 < a2
subject to
(i) Boundary condition: u(a, ) = f (), 0 2.
(ii) Boundedness at the origin: |u(0, )| < .
(iii) Periodicity: u(r, + 2) = u(r, ), 0 2.
(20.1)
R0 (r)
00 ()
R00 (r)
= r2
+r
= .
()
R(r)
R(r)
(20.2)
(20.3)
and
The second equation is known as Eulers equation. Both of these equations
are easily solvable. To solve (20.2), We only have to add the appropriate
boundary conditions. We have (0) = (2) and 0 (0) = 0 (2). The
periodicity of implies that = n2 and must be of the form
n () = A0n cos n + Bn0 sin n, n = 0, 1, 2
The equation in R is of Euler type and its solution must be of the form
R(r) = r . Since = n2 , the corresponding characteristic equation is
( 1)r + r n2 r = 0.
Solving this equation we find = n. Hence, we let
Rn (r) = Cn rn + Dn rn , n N.
167
X
rn (An cos n + Bn sin n).
u(r, ) = C0 +
n=1
X
f () = C0 +
(an An cos n + an Bn sin n)
n=1
f () =
a0 X
+
(an cos n + bn sin n).
2
n=1
It is obvious that an and bn are the Fourier coefficients, and therefore can be
determined by the formulas
Z
1 2
an =
f () cos nd, n = 0, 1,
0
and
1
bn =
f () sin nd, n = 1, 2, .
0
n=1
where
Z 2
a0
1
C0 = =
f ()d
2
2 0
Z 2
1
an
f () cos nd, n = 1, 2,
An = n = n
a
a 0
Z 2
bn
1
f () sin nd, n = 1, 2,
Bn = n = n
a
a 0
Example 20.2
Solve
u = 0, 0 < 2, 1 r 2
subject to
u(1, ) = u(2, ) = sin , 0 < 2.
Solution.
Use separation of variables. First, solving for (), we see that in order to
ensure that the solution is 2periodic in , the eigenvalues are = n2 .
When solving the equation for R(r), we do NOT need to throw out solutions
which are not bounded as r 0. This is because we are working in the
annulus where r is bounded away from 0 and . Therefore, we obtain the
general solution
u(r, ) = (C0 + C1 ln r) +
n=1
But
C0 +
X
[(Cn + Dn ) cos n + (An + Bn ) sin n] = sin
n=1
and
C0 +
X
n=1
169
=0
=0
=0 (n 6= 1)
=1
=0
=0 (n =
6 1)
=1.
() = A cos + B sin .
The condition 0 (0) = 0 implies B = 0. The condition 0 3 = 0 implies
n = (3n)2 , n = 0, 1, 2, . Thus, the angular solution is
n () = An cos 3n, n = 0, 1, 2,
n=0
1
1
cos 9 cos 3
3
9
we obtain D1 a3 = 19 and D3 a9 =
u(a, ) =
1
3
1 r 9
1 r 3
cos 9
cos 3
3 a
9 a
171
Practice Problems
Problem 20.1
Solve the Laplaces equation as in Example 20.1 in the unit disk with u(1, ) =
3 sin 5.
Problem 20.2
Solve the Laplaces equation in the upper half of the unit disk with u(1, ) =
.
Problem 20.3
Solve the Laplaces equation in the unit disk with ur (1, ) = 2 cos 2.
Problem 20.4
Consider
u(r, ) = C0 +
n=1
with
Z 2
1
a0
f ()d
C0 = =
2
2 0
Z 2
an
1
An = n = n
f () cos nd, n = 1, 2,
a
a 0
Z 2
bn
1
Bn = n = n
f () sin nd, n = 1, 2,
a
a 0
Using the trigonometric identity
cos a cos b + sin a sin b = cos (a b)
show that
1
u(r, ) =
2
"
f () 1 + 2
X
r n
n=1
#
cos n( ) d.
Problem 20.5
(a) Using Eulers formula from complex analysis eit = cos t + i sin t show that
1
cos t = (eit + eit ),
2
where i = 1.
(b) Show that
1+2
X
r n
n=1
cos n( ) = 1 +
X
r n
n=1
in()
X
r n
n=1
ein() .
X
r n
n=1
rei()
a rei()
ein() =
rei()
a rei()
n=1
and
ein() =
Hint: Each sum is a geoemtric series with a ratio less than 1 in absolute
value so that these series converges.
(b) Show that
1+2
X
r n
n=1
cos n( ) =
a2 r 2
.
a2 2ar cos ( ) + r2
Problem 20.7
Show that
a2 r 2
u(r, ) =
2
Z
0
a2
f ()
d.
2ar cos ( ) + r2
173
Problem 20.9
The vibrations of a symmetric circular membrane where the displacement
u(r, t) depends on r and t only can be describe by the one-dimensional wave
equation in polar coordinates
1
utt = c2 (urr + ur ), 0 < r < a, t > 0
r
with initial condition
u(a, t) = 0, t > 0
and boundary conditions
u(r, 0) = f (r), ut (r, 0) = g(r), 0 < r < a.
(a) Show that the assumption u(r, t) = R(r)T (t) leads to the equation
1 00 1 R0
1 T 00
=
R +
= .
c2 T
R
rR
(b) Show that < 0.
Problem 20.10
Cartesian coordinates and cylindrical coordinates are shown in Figure 20.1
below.
Figure 20.1
(x,y)
(x,y)
(x,y)
175
177
178
Figure 21.1
In this section we introduce the concept of Laplace transform and discuss
some of its properties.
The Laplace transform is defined in the following way. Let f (t) be defined
for t 0. Then the Laplace transform of f, which is denoted by L[f (t)]
or by F (s), is defined by the following equation
Z
L[f (t)] = F (s) = lim
T
st
f (t)e
0
Z
dt =
f (t)est dt
179
Solution.
(a) Using the definition of Laplace transform we see that
Z
Z T
(sa)t
at
e
dt = lim
e(sa)t dt.
L[e ] =
T
But
(sa)t
dt =
T
1e(sa)T
sa
if s = a
if s 6= a.
1
, s > a.
sa
st
te
0
test est
2
dt =
s
s
=
0
1
, s > 0.
s2
R 2
2
t2 st
If
1 and this implies that 0 et st dt
R s 0 then t st 0 so that e
dt. Since the integral on the right is divergent, by the comparison theo0
rem of improper
integrals,Rthe integral on the left is also divergent. Now, if
R t(ts)
s > 0 then 0 e
dt s dt. By the same reasoning the integral on the
180
left is divergent. This shows that the function f (t) = et does not possess a
Laplace transform
The above example raises the question of what class or classes of functions
possess a Laplace transform. To answer this question we introduce few mathematical concepts.
A function f that satisfies
|f (t)| M eat ,
tC
(21.1)
Solution.
P
tn
tn
n
t
n
(a) Since et =
n=0 n! n! , t n!e . Hence, t is piecewise continuous and
of exponential order at infinity.
(b) Since |tn sin at| n!et , tn sin at is piecewise continuous and of exponential
order at infinity
181
exists as long as s > a. Note that the two conditions above are sufficient, but
not necessary, for F (s) to exist.
In what follows, we will denote the class of all piecewise continuous functions
with exponential order at infinity by PE. The next theorem shows that any
linear combination of functions in PE is also in PE. The same is true for the
product of two functions in PE.
Theorem 21.2
Suppose that f (t) and g(t) are two elements of PE with
|f (t)| M1 ea1 t ,
t C1
and
|g(t)| M2 ea1 t ,
t C2 .
(i) For any constants and the function f (t) + g(t) is also a member of
PE. Moreover
L[f (t) + g(t)] = L[f (t)] + L[g(t)].
(ii) The function h(t) = f (t)g(t) is an element of PE.
We next discuss the problem of how to determine the function f (t) if F (s)
is given. That is, how do we invert the transform. The following result on
uniqueness provides a possible answer. This result establishes a one-to-one
correspondence between the set PE and its Laplace transforms. Alternatively, the following theorem asserts that the Laplace transform of a member
in PE is unique.
Theorem 21.3
Let f (t) and g(t) be two elements in PE with Laplace transforms F (s) and
G(s) such that F (s) = G(s) for some s > a. Then f (t) = g(t) for all t 0
where both functions are continuous.
182
With the above theorem, we can now officially define the inverse Laplace
transform as follows: For a piecewise continuous function f of exponential
order at infinity whose Laplace transform is F, we call f the inverse Laplace
transform of F and write f = L1 [F (s)]. Symbolically
f (t) = L1 [F (s)] F (s) = L[f (t)].
Example 21.3
1
, s > 1.
Find L1 s1
Solution.
1
From Example 21.1(a), we have that L[eat ] = sa
, s > a. In particular, for
1
1
t
1
a = 1 we find that L[e ] = s1 , s > 1. Hence, L
= et , t 0 .
s1
The above theorem states that if f (t) is continuous and has a Laplace transform F (s), then there is no other function that has the same Laplace transform. To find L1 [F (s)], we can inspect tables of Laplace transforms of
known functions to find a particular f (t) that yields the given F (s).
When the function f (t) is not continuous, the uniqueness of the inverse
Laplace transform is not assured. The following example addresses the
uniqueness issue.
Example 21.4
Consider the two functions f (t) = H(t)H(3 t) and g(t) = H(t) H(t 3),
where H is the Heaviside function defined by
1, t 0
H(t) =
0, t < 0
(a) Are the two functions identical?
(b) Show that L[f (t)] = L[g(t).
Solution.
(a) We have
1, 0 t 3
0,
t>3
1, 0 t < 3
0,
t3
f (t) =
and
g(t) =
183
Example 21.5
Find f g where f (t) = et and g(t) = sin t.
Solution.
Using integration by parts twice we arrive at
Z t
(f g)(t) =
e(ts) sin sds
0
t
1
= e(ts) (sin s cos s) 0
2
et 1
=
+ (sin t cos t)
2
2
(21.2)
184
Theorem 21.5
Let f (t), g(t), and k(t) be three piecewise continuous scalar functions defined
for t 0 and c1 and c2 are arbitrary constants. Then
(i) f g = g f (Commutative Law)
(ii) (f g) k = f (g k) (Associative Law)
(iii) f (c1 g + c2 k) = c1 f g + c2 f k (Distributive Law)
Example 21.6
Express the solution to the initial value problem y 0 + y = g(t), y(0) = y0
in terms of a convolution integral.
Solution.
Solving this initial value problem by the method of integrating factor we find
Z t
t
e(ts) g(s)ds = et y0 + (et g)(t)
y(t) = e y0 +
0
P (s) =
1
2
s + a2
1
2
s + a2
.
1
1
So, in this case we have, F (s) = G(s) = s2 +a
2 so that f (t) = g(t) = a sin (at).
Thus,
Z
1 t
1
(f g)(t) = 2
sin (at as) sin (as)ds = 3 (sin (at) at cos (at))
a 0
2a
185
Example 21.8
Solve the initial value problem
4y 00 + y = g(t), y(0) = 3, y 0 (0) = 7
Solution.
Take the Laplace transform of all the terms and plug in the initial conditions
to obtain
4(s2 Y (s) 3s + 7) + Y (s) = G(s)
or
(4s2 + 1)Y (s) 12s + 28 = G(s).
Solving for Y (s) we find
G(s)
12s 28
+
4 s2 + 14
4 s2 + 41
1 2
1 2
3s
1
2
2
=
7
+ G(s)
2
1 2
1 2
2
2
2
4
s + (2
s + 2
s + 21
Y (s) =
Hence,
Z
s
t
t
1 t
sin
g(t s)ds.
y(t) = 3 cos
7 sin
+
2
2
2 0
2
So, once we decide on a g(t) all we need to do is to evaluate the integral and
well have the solution
We conclude this section with the following table of Laplace transform pairs
where H is the Heaviside function defined by H(t) = 1 for t 0 and 0
otherwise.
186
F(s)
H(t) =
1, t 0
0, t < 0
1
s,
tn , n = 1, 2,
et
sin (t)
cos (t)
sinh (t)
cosh (t)
et f (t), with |f (t)| M eat
et H(t)
et tn , n = 1, 2,
et sin (t)
et cos (t)
f (t )H(t ), 0
with |f (t)| M eat
s>0
n!
, s>0
sn+1
1
,
s s >
, s>0
s2 + 2
s
, s>0
s2 + 2
, s > ||
s2 2
s
, s > ||
s2 2
F (s ), s > + a
1
s , s >
n!
, s>
(s)n+1
, s>
(s)2 + 2
s
, s>
(s)2 + 2
s
e
F (s), s > a
H(t ), 0
tf (t)
es
s , s
-F 0 (s)
t
2
s
,
(s2 + 2 )2
s>0
1
,
(s2 + 2 )2
s>0
sin t
1
[sin t
2 3
t cos t]
>0
sF (s) f (0)
s > max{a, 0} + 1
Rt
F (s)
s ,
Table L
s > max{a, 0} + 1
187
Practice Problems
Problem 21.1
R
Determine whether the integral 0
verges, give its value.
1
dt
1+t2
Problem 21.2
R
Determine whether the integral 0
verges, give its value.
t
dt
1+t2
Problem 21.3
R
Determine whether the integral 0 et cos (et )dt converges. If the integral
converges, give its value.
Problem 21.4
Using the definition, find L[e3t ], if it exists. If the Laplace transform exists
then find the domain of F (s).
Problem 21.5
Using the definition, find L[t 5], if it exists. If the Laplace transform exists
then find the domain of F (s).
Problem 21.6
2
Using the definition, find L[e(t1) ], if it exists. If the Laplace transform
exists then find the domain of F (s).
Problem 21.7
Using the definition, find L[(t 2)2 ], if it exists. If the Laplace transform
exists then find the domain of F (s).
Problem 21.8
Using the definition, find L[f (t)], if it exists. If the Laplace transform exists
then find the domain of F (s).
f (t) =
0,
0t<1
t 1,
t1
188
Problem 21.9
Using the definition, find L[f (t)], if it exists. If the Laplace transform exists
then find the domain of F (s).
0t<1
0,
t 1, 1 t < 2
f (t) =
0,
t 2.
Problem 21.10
Let n be a positive integer. Using integration by parts establish the reduction
formula
Z
Z
tn est n
n st
+
tn1 est dt, s > 0.
t e dt =
s
s
Problem 21.11
For s > 0 and n a positive integer evaluate the limits
(a) limt0 tn est
Problem 21.12
Use the linearity property of Laplace transform to find L[5e7t + t + 2e2t ].
Find the domain of F (s).
Problem 21.13
3
.
Find L1 s2
Problem 21.14
Find L1 s22 +
1
s+1
Problem 21.15
2
2
Find L1 s+2
+ s2
.
Problem 21.16
Use Table L to find L[2et + 5].
Problem 21.17
Use Table L to find L[e3t3 H(t 1)].
Problem 21.18
Use Table L to find L[sin2 t].
4
].
s3
Problem 21.23
5
Use Table L to find L1 [ (s3)
4 ].
Problem 21.24
2s
Use Table L to find L1 [ es9 ].
Problem 21.25
h
i
12
.
Using the partial fraction decomposition find L1 (s3)(s+1)
Problem 21.26
h
i
1 24e5s
Using the partial fraction decomposition find L
.
s2 9
Problem 21.27
Use Laplace transform technique to solve the initial value problem
y 0 + 4y = g(t), y(0) = 2
where
0, 0 t < 1
12, 1 t < 3
g(t) =
0,
t3
Problem 21.28
Use Laplace transform technique to solve the initial value problem
y 00 4y = e3t , y(0) = 0, y 0 (0) = 0.
189
190
Problem 21.29
Consider the functions f (t) = et and g(t) = e2t , t 0. Compute f g in
two different ways.
(a) By directly evaluating the integral.
(b) By computing L1 [F (s)G(s)] where F (s) = L[f (t)] and G(s) = L[g(t)].
Problem 21.30
Consider the functions f (t) = sin t and g(t) = cos t, t 0. Compute f g in
two different ways.
(a) By directly evaluating the integral.
(b) By computing L1 [F (s)G(s)] where F (s) = L[f (t)] and G(s) = L[g(t)].
Problem 21.31
Compute t t t.
Problem 21.32
Compute H(t) et e2t .
Problem 21.33
Compute t et et .
191
The time derivatives are transformed in the same way as in the case of
functions in one variable, that is, for example
L(ut )(x, t) = sU (x, s) u(x, 0)
and
L(utt )(x, s) = s2 U (x, s) su(x, 0) ut (x, 0).
The spatial derivatives remain unchanged, for example,
Z
Z
s
ux (x, )e d =
u(x, )es d = Ux (x, s).
Lux (x, t) =
x 0
0
Likewise, we have
Luxx (x, t) = Uxx (x, s).
Thus, applying the Laplace transform to a PDE in two variables x and t we
obtain an ODE in the variable x and with the parameter s.
Example 22.1
Let u(x, t) be the concentration of a chemical contaminant dissolved in a
liquid on a half-infinte domain x > 0. Let us assume that at time t = 0 the
concentration is 0 and on the boundary x = 0, constant unit concentration of
the contaminant is kept for t > 0. The behaviour of this problem is described
by the following mathematical model
u(x, 0) = 0,
u(0, t) = 1,
192
Solution.
Applying Laplace transform to both sides of the equation we obtain
sU (x, s) u(x, 0) Uxx (x, s) = 0
or
Uxx (x, s) sU (x, s) = 0.
This is a second order linear ODE in the variable x and positive parameter
s. Its general solution is
U (x, s) = A(s)e
sx
+ B(s)e
sx
U (x, s) = A(s)e sx .
Next, we apply Laplace transform to the boundary condition obtaining
1
U (0, s) = L(1) = .
s
This leads to A(s) =
1
s
Thus,
u(x, t) = L
1 sx
e
s
2
=
ew dw
x
2 t
Example 22.2
Solve the following initial boundary value problem
u(x, 0) = 0,
u(0,
t) = f (t),
sx
, F (s) = Lf (t).
193
1
sx
u(x, t) = L
F (s)e
= f L1 (e sx ).
It can be shown that
L1 (e
sx
x2
x
)=
e 4t .
4t3
Hence,
t
Z
u(x, t) =
0
x2
x
p
e 4(ts) f (s)ds
4(t s)3
Example 22.3
Solve the wave equation
utt c2 uxx = 0
, x > 0, t > 0
u(x, 0) = ut (x, 0) = 0,
u(0, t) = f (t),
U (x, s) = A(s)e c x .
Next, we apply Laplace transform to the boundary condition obtaining
U (0, s) = L(f (t)) = F (s).
194
U (x, s) = F (s)e c x .
Thus,
x
x
x
f t
u(x, t) = L1 F (s)e c s = H t
c
c
Remark 22.1
Laplace transforms are useful in solving parabolic and some hyperbolic PDEs.
They are not in general useful in solving elliptic PDEs.
195
Practice Problems
Problem 22.1
Solve by Laplace transform
, x > 0, t > 0
ut + ux = 0
u(x, 0) = sin x,
u(0, t) = 0
Hint: Method of integrating factor of ODEs.
Problem 22.2
Solve by Laplace transform
ut + ux = u , x > 0, t > 0
u(x, 0) = sin x,
u(0, t) = 0
Problem 22.3
Solve
ut = 4uxx
u(0, t) = u(1, t) = 0
u(x, 0) = 2 sin x + 6 sin 2x.
Hint: A particular solution of a second order ODE must be found using the
method of variation of parameters.
Problem 22.4
Solve by Laplace transform
ut ux = u , x > 0, t > 0
u(x, 0) = e5x ,
ut + ux = t , x > 0, t > 0
u(x, 0) = 0,
u(0, t) = t2
196
Problem 22.6
Solve by Laplace transform
u(0, t) = t
Problem 22.7
Solve by Laplace transform
utt c2 uxx = 0
, x > 0, t > 0
u(x, 0) = ut (x, 0) = 0,
u(0, t) = sin t,
uxy = 1
, x > 0, y > 0
u(x, 0) = 1,
u(0, y) = y + 1.
Problem 22.10
Solve by Laplace transform
utt = c2 uxx
, x > 0, t > 0
u(x, 0) = ut (x, 0) = 0,
ux (0, t) = f (t),
, x > 0, t > 0
ut + ux = u
u(x, 0) = sin x,
u(0, t) = 0
Problem 22.12
Solve by Laplace transform
u(x, 0) = T,
u(0, t) = 0,
197
198
199
200
n
n
a0 X
f (x) =
+
an cos
x + bn sin
x
2
L
L
n=1
at all points of continuity of f. In the context of Fourier analysis, this is
referred to as the real form of the Fourier series. It is often convenient to
recast this series in complex form by means of Euler formula
eix = cos x + i sin x.
It follows from this formula that
eix + eix = 2 cos x and eix eix = 2i sin x
or
cos x =
eix +eix
2
and sin x =
eix eix
.
2i
a0 X
e L + e L
an
f (x) = +
2
2
n=1
!#
inx
inx
e L e L
+bn
2i
f (x) =
cn e
inx
L
n=
where c0 =
a0
2
an ibn
2
an + ibn
.
2
(23.1)
201
and bn = i(cn cn ).
(23.2)
That is, an and bn can be easily found once we have formulas for cn . In order
to find these formulas, we need to evaluate the following integral
Z
inx
L
imx
L
Z
dx =
i(nm)x
L
dx
i
i(nm)x L
L
L
e
=
i(n m)
L
iL
=
[cos [(n m)] + i sin [(n m)]
(n m)
cos [(n m)] i sin [(n m)]]
=0
if n 6= m. If n = m then
Z
inx
L
inx
L
dx = 2L.
L
inx
which yields the formula for coefficients of the complex form of the Fourier
series:
Z L
inx
1
cn =
f (x)e L dx, n = 0, 1, 2, .
2L L
Example 23.1
Find the complex Fourier coefficients of the function
f (x) = x,
extended to be periodic of period 2.
202
Solution.
Using integration by parts and the fact that ei = ei = 1 we find
Z
1
xeinx dx
cn =
2
Z
i
ix inx
1
inx
e
e
dx
=
2
n
n
1
i in
i in
=
e
+
e
2
n
n
1 in
1 1 in
e
2e
+
2 n2
n
h
i
1
1
(1)n i
=
2i (1)n +
(0) =
2
n
2
n
for n N and for n = 0, we have
1
c0 =
2
xdx = 0
Remark 23.1
It is often the case that the complex form of the Fourier series is far simpler
to calculate than the real form. One can then use (23.2) to find the real form
of the Fourier series. For example, the Fourier coefficients of the real form of
the previous function are given by
an = (cn + cn ) = 0 and bn = i(cn cn ) = n2 (1)n+1 , n N
203
Practice Problems
Problem 23.1
Find the complex Fourier coefficients of the function
f (x) = x,
1x1
0 < x <
2
1
<
x
<
f (x) =
2
2
0
<x<
204
Problem 23.8
Let f (x) = 2 x, 2 < x < 2, be of period 4.
(a) Calculate the coefficients an , bn and cn .
(b) Find the complex Fourier series representation of f.
Problem 23.9
Suppose that the coefficients cn of the complex Fourier series are given by
2
if |n| is odd
in
cn =
0 if |n| is even.
Find an , n = 0, 1, 2, and bn , n = 1, 2, .
Problem 23.10
Recall that any complex number z can be written as z = Re(z) + iIm(z)
where Re(z) is called the real part of z and Im(z) is called the imaginary
part. The complex conjugate of z is the complex number z = Re(z)
iIm(z). Using these definitions show that an = 2Re(cn ) and bn = 2Im(cn ).
Problem 23.11
Suppose that
cn =
i
[einT
2n
T
2
1] if n 6= 0
if n = 0.
Find an and bn .
Problem 23.12
Find the complex Fourier series of the function f (x) = ex on [2, 2].
Problem 23.13
Consider the wave form
205
206
fe (x) =
cn e
inx
T
(24.1)
n=
where
1
cn =
2T
fe (x)e
inx
T
dx.
Let R. Multiply both sides of (24.1) by eix and then integrate both sides
from T to T. Assuming integration and summation can be interchanged
we find
Z T
Z T
X
inx
ix
fe (x)e
dx =
cn
eix e T dx.
T
n=
207
f (x)eix dx.
(24.2)
f () =
The function f is called the Fourier transform of f. We will use the notation
F[f (x)] = f().
Next, it can be shown that
n
f
= 2T cn
T
so that
1 X n inx
f
eT .
fe (x) =
2T n=
T
1 X n inx
T
f
e =
f()eix d
lim
T T
T
n=
so that
1
f (x) =
2
f()eix d
(24.3)
Equation (24.3) is called the Fourier inversion formula and we use the
notation F 1 [f()]. Now, if we make use of Eulers formula, we can write the
Fourier inversion formula in terms of sines and cosines,
Z
Z
i
1
f (x) =
f() cos xd +
f() sin xd
2
2
a superposition of sines and cosines of various frequencies.
Example 24.1
Find the Fourier transform of the function f (x) defined by
ax
e
if x 0
f (x) =
0
if x < 0
for some a > 0.
208
Solution.
We have
f() =
ix
eax eix dx
0
Z
x(a+i)
e
axix
=
e
dx =
(a + i) 0
0
1
=
a + i
f (x)e
dx =
The following theorem lists the basic properties of the Fourier transform
Theorem 24.1
Let f, g, be piecewise continuous functions. Then we have the following
properties:
(1) Linearity: F[f (x) + g(x)] = F[f (x)] + F[g(x)], where and are
arbitrary numbers.
(2) Shifting: F[f (x )] = ei F[f (x)].
(3) Scaling: F[f x ] = F[f (x)].
R
(4) Continuity: If |f (x)|dx < then f is continuous in .
n
(5) Differentiation:RF[f (n) (x)]
F[f (x)].
= (i)
x
1
(6) Integration: F 0 f (s)ds = i F[f (x)].
R
R
1
|f()|2 d.
(7) Parsevals Relation: |f (x)|2 dx = 2
u() =
ex eix dx.
If we differentiate this relation with respect to the variable and then integrate by parts we obtain
209
xex eix dx
u () = i
Z
i
d x2 ix
(e
)e
dx
=
2 dx
Z
i
x2 ix
x2 ix
(e
)e
dx
=
+ i
e
2
Z
i2 x2 ix
=
(e
)e
dx = u()
2
2
0
2
ex dx = ,
we can write
Z
u(0) =
x2
r
dx =
and therefore
r
u() =
= C,
2
e 4
Example 24.3
Prove
F[f (x)] = f().
Solution.
Using a change of variables we find
Z
F[f (x)] =
f (x)e
ix
dx =
f (x)eix dx = f()
Example 24.4
Prove
F[F[f (x)]] = 2f (x).
210
Solution.
We have
1
f (x) =
2
Thus,
Z
f()eix d
2f (x) =
x2
1
e 4 .
(50 ) Gaussians: F 1 [e ] = 4
(60 ) Product: F 1 [f()
g ()] = f (x) g(x).
(70 ) Convolution: F 1 [f g()] = 2(f g)(x).
Remark 24.1
It is important to mention that there exists no established convention of how
to define the Fourier transform. In the literature, we can meet an equivalent
1
in front of the integral.
definition of (24.3) with the constant 12 or 2
There also exist definitions with positive sign in the exponent. The reader
should keep this fact in mind while working with various sources or using the
transformation tables.
211
Practice Problems
Problem 24.1
Find the Fourier transform of the function
1 if 1 x 1
f (x) =
0
otherwise.
Problem 24.2
Obtain the transformed problem when applying the Fourier transform with
respect to the spatial variable to the equation and initial condition
ut + cux = 0
u(x, 0) = f (x).
Problem 24.3
Obtain the transformed problem when applying the Fourier transform with
respect to the spatial variable to the equation and both initial conditions
utt = c2 uxx , x R, t > 0
u(x, 0) = f (x)
ut (x, 0) = g(x).
Problem 24.4
Obtain the transformed problem when applying the Fourier transform with
respect to the spatial variable to the equation and both initial conditions
u = uxx + uyy = 0, x R, 0 < y < L
u(x, 0) = 0
1 if a < x < a
u(x, L) =
0
otherwise
Problem 24.5
Find the Fourier transform of f (x) = e|x| , where > 0.
212
Problem 24.6
Prove that
F[ex H(x)] =
where
H(x) =
Problem 24.7
Prove that
1
1 + i
1 if x 0
0 otherwise.
1
= 2e H().
F
1 + ix
Problem 24.8
Prove
F[f (x )] = ei f().
Problem 24.9
Prove
F[eix f (x)] = f(x ).
Problem 24.10
Prove the following
1
F[cos (x)f (x)] = [f( + ) + f( )]
2
1
F[sin (x)f (x)] = [f( + ) f( )]
2
Problem 24.11
Prove
F[f 0 (x)] = (i)f().
Problem 24.12
Find the Fourier transform of f (x) = 1 |x| for 1 x 1 and 0 otherwise.
Problem 24.13
Find, using the definition, the Fourier transform of
1 a < x < 0
1
0<x<a
f (x) =
0
otherwise
213
214
215
u(x, 0) = f (x)
ut (x, 0) = g(x).
Again, by performing the Fourier transform of u in x, we reduce the PDE
problem into an ODE problem in the variable t:
2 u
= c2 2 u
2
t
u(, 0) = f()
ut (, 0) = g().
General solution to the ODE is
u(, t) = ()eict + ()eict
where and are two arbitrary functions of . Performing the inverse
transformation and making use of the translation theorem, we get the general
solution
u(x, t) = (x ct) + (x + ct)
where F() = and F() = . But
1
1
f ()
g()
() =
2
ic
1
1
f () +
g() .
() =
2
ic
By using the integration property, we find the inverse transforms of and
Z
1
1 x
(x) =
f (x) +
g(s)ds
2
c 0
Z
1
1 x
(x) =
g(s)ds .
f (x)
2
c 0
Application of the translation property then yields directly the DAlambert
solution
Z
1
1 x+ct
u(x, t) = [f (x ct) + f (x + ct)] +
g(s)ds.
2
2c xct
216
217
Using the boundary condition u(, 0) = 0 we find B() = 0. Using the second
boundary condition we find
Z
u(x, L)eix dx
u(, L) =
Z
Z a
a
ix
=
e
dx =
cos (x)dx
a
2 sin (a)
.
=
Hence,
A() sinh (L) =
2 sin (a)
2 sin (a)
.
sinh (L)
Thus,
u(, y) =
2 sin (a)
sinh (y).
sinh (L)
218
Practice Problems
Problem 25.1
Solve, by using Fourier transform
ut + cux = 0
x2
u(x, 0) = e 4 .
Problem 25.2
Solve, by using Fourier transform
ut = kuxx u, x R
x2
u(x, 0) = e .
Problem 25.3
Solve the heat equation
ut = kuxx
subject to the initial condition
u(x, 0) =
1 if x 0
0 otherwise.
Problem 25.4
Use Fourier transform to solve the heat equation
ut = uxx + u,
<x<< t>0
u(x, 0) = f (x).
Problem 25.5
Prove that
e||y eix d =
x2
2y
.
+ y2
Problem 25.6
Solve the Laplaces equation in the half plane
uxx + uyy = 0,
x2 +y 2
u(x, y) = 0.
219
220
Appendix
221
222
APPENDIX
(26.1)
(26.2)
6
16
cos 3t
sin 3t.
73
73
Next, the
characteristic equation r2 r + 1 = 0 has roots r1 = 21 i 23 and
r2 = 12 + i 23 . Thus, the general solution to the homogeneous equation is
1
3
3
yh (t) = e 2 t (c1 cos
t + c2 sin
t).
2
2
The general solution to the differential equation is
1
3
3
6
16
y(t) = yh (t) + yp (t) = e 2 t (c1 cos
t + c2 sin
t) +
cos 3t
sin 3t
2
2
73
73
224
APPENDIX
Example 26.3
Find the general solution of
y 00 + 4y 0 2y = 2t2 3t + 6.
Solution.
We see from the previous two examples that the trial function has usually
the appearance of the nonhomogeneous term g(t). Since g(t) is a quadratic
function, we are going to try yp (t) = At2 + Bt + C. Inserting this into the
differential equation leads to
2At2 + (8A 2b)t + (2A + 4B 2C) = 2t2 3t + 6.
Equating coefficients of like powers of t we find A = 1, B = 25 , and
C = 9. Thus, a particular solution is
5
yp (t) = t2 t 9.
2
We next solve the homogeneous equation.
The characteristic
equation r2 +
yh (t) = c1 e(2
6)t
+ c2 e(2+
6)t
6)t
+ c2 e(2+
6)t
5
t2 t 9
2
Remark 26.1
The same principle used in the previous three examples extends to the case
where g(t) is a product of any two or all three of the three types of functions
discussed above, as the next example illustrates.
Example 26.4
Find a particular solution of
y 00 3y 0 4y = 8et cos 2t.
4
23
4 2t
+ c2 e t +
2t +
e
3
9
3
3t
226
APPENDIX
Example 26.6
Find the general solution of the nonhomogeneous equation
y 00 y 0 2y = 4et .
Solution.
Lets try with yp (t) = Aet . Substituting this into the differential equation
leads to 0Aet = 4et . Thus, A does not exist. Why did the procedure of the
previous examples fail here? The reason is that the function et that appears
in yp is a solution to the homogeneous equation and so cannot possibly be
a solution to the nonhomogeneous equation at the same time. Then comes
the question of how to find a correct form for the particular solution.
We will try to solve a simpler equation with the same difficulty and to use
its general solution to suggest how to proceed with our given equation. The
simpler equation we consider is y 0 + y = 4et . By the method of integrating
factor we find the general solution y(t) = 4tet +cet . The second term is the
solution to the homogeneous equation whereas the first one is the solution
to the nonhomogeneous equation. We conclude from this discussion that a
good guess for the original equation would be yp (t) = Atet . If we insert
this into the differential equation we end up with 3Aet = 4et . Solving
for A we find A = 34 . Thus, yp (t) = 43 tet and the general solution to the
differential equation is y(t) = c1 et + c2 e2t 43 tet
Example 26.7
Find the general solution of the nonhomogeneous equation
y 00 + 2y 0 + y = 2et .
Solution.
The characteristic equation is r2 +2r +1 = 0 with double roots r1 = r2 = 1.
Thus, yh (t) = c1 et + c2 tet . Our trial function can not contain either et or
tet since both are solutions to the homogeneous equation. Thus, a proper
guess is yp (t) = At2 et . Finding derivatives up to order 2 we find yp0 (t) =
2Atet At2 et and yp00 (t) = 2Aet 4Atet + At2 et . Substituting this in
the original equation and collecting like terms we find
2Aet = 2et .
Solving for A we find A = 1 so that yp (t) = t2 et . Hence, the general solution
is given by
y(t) = c1 et + c2 tet + t2 et
Form of yp (t)
tr [An tn + An1 tn1 + + A1 t + A0 ]
tr [An tn + An1 tn1 + + A1 t + A0 ]et
tr et [(An tn + An1 tn1 + + A1 t + A0 ) cos t
+(Bn tn + Bn1 tn1 + + B1 t + B0 ) sin t]
228
APPENDIX
Solution.
First, the characteristic equation is r2 + 1 = 0, with roots r = i, so the
homogeneous solution is yh (t) = c1 sin t + c2 cos t. The trial function for the
particular solution is yp (t) = Aet + Bt3 + Ct2 + Dt + E. Plugging into the
differential equation, we see
Aet + 6Bt + 2C + Aet + Bt3 + Ct2 + Dt + E = et + t3 .
Matching coefficients, we see:
2A = 1, B = 1, C = 0, 6B + D = 0, E = 0
The particular solution is
1
yp (t) = et + t3 6t,
2
and so the general solution is
1
y(t) = c1 sin t + c2 cos t + et + t3 6t
2
When t = 0, this is y(0) = c2 + 21 = 2, so c2 = 32 . The first derivative of the
general solution is y 0 (t) = c1 cos t 23 sin t + 21 et + 3t2 6. At t = 0, y 0 (0) =
c1 + 21 6 = 0, so c1 = 11
. We thus have solution y(t) = 11
sin t + 23 cos t +
2
2
1 t
e + t3 6t
2
(27.1)
(27.2)
230
APPENDIX
and
yp00 = y100 u1 + y10 u01 + y200 u2 + y20 u02 .
In particular, yp00 does not involve u001 and u002 .
Inserting yp , yp0 , and yp00 into equation (27.1) to obtain
[y100 u1 + y10 u01 + y200 u2 + y20 u02 ] + p(t)(y10 u1 + y20 u2 ) + q(t)(u1 y1 + u2 y2 ) = g(t).
Rearranging terms,
[y100 + p(t)y10 + q(t)y1 ]u1 + [y200 + p(t)y20 + q(t)y2 ]u2 + [u01 y10 + u02 y20 ] = g(t).
Since y1 and y2 are solutions to the homogeneous equation, the previous
equation yields our second constraint
u01 y10 + u02 y20 = g(t).
(27.3)
Combining equation (27.2) and (27.3) we find the system of two equations
in the unknowns u01 and u02
y1 u01 + y2 u02 =0
u01 y10 + u02 y20 =g(t).
Since {y1 , y2 } is a fundamental set, the expression W (t) = y1 y20 y10 y2 is
nonzero so that one can find unique u01 and u02 . Using the method of elimination, these functions are given by
(t)g(t)
u01 (t) = y2W
and u02 (t) =
(t)
y1 (t)g(t)
.
W (t)
y1 (t)g(t)
dt.
W (t)
Z
u2 (t) =
et 2et
2
dt = e3t .
t
3e
9
4
4t 0
y +
y = (2t 1)et .
2t 1
2t 1
Z
u2 (t) =
t (2t 1)et
1 3t 1 3t
dt
=
te e .
(2t 1)e2t
3
9
Thus,
1
1
2
1
yp (t) = tet tet et = tet et .
3
9
3
9
The general solution is
2
1
y(t) = c1 t + c2 e2t + tet et
3
9
232
APPENDIX
Example 27.3
Find the general solution to the differential equation y 00 + y 0 = ln t, t > 0.
Solution.
The characterisitc equation r2 + r = 0 has roots r1 = 0 and r2 = 1 so that
y1 (t) = 1, y2 (t) = et , and W (t) = et . Hence,
Z t
Z
e ln t
u1 (t) =
dt = ln tdt = t ln t t
et
Z
Z
Z t
ln t
e
t
t
dt
u2 (t) =
dt = e ln tdt = e ln t +
t
e
t
Thus,
Z t
e
t
dt
yp (t) = t ln t t ln t + e
t
and
Z t
e
t
t
y(t) = c1 + c2 e + t ln t t ln t + e
dt
t
Example 27.4
Find the general solution of
1
y 00 + y =
.
2 + sin t
Solution.
Since the characteristic equation r2 + 1 = 0 has roots r = i, the general
solution of the corresponding homogeneous equation y 00 + y = 0 is given by
yh (t) = c1 cos t + c2 sin t
Since W (t) = 1 we find
Z
sin t
2
dt = t +
dt
u1 (t) =
2 + sin t
2 + sin t
Z
cos t
u2 (t) =
dt = ln (2 + sin t)
2 + sin t
Hence, the particular solution is
Z
2
yp (t) = sin t ln (2 + sin t) + cos t(
dt t)
2 + sin t
and the general solution is
Z
234
Section 2
2.1 a = b = 0.
2.2 Substituting into the differential equation we find
tX 00 T XT 0 = 0
or
T0
X 00
=
.
X
tT
The LHS is a function of x only whereas the RHS is a function of t only.
This is true only when both sides are constant. That is, there is such that
X 00
T0
=
=
X
tT
and this leads to the two ODEs X 00 = X and T 0 = tT.
2.3 We have xux + (x + 1)yuy =
x x
(e
y
+ xe ) + (x + 1)y
x
xe
y2
= 0 and
235
u(1, 1) = e.
2.4 We have ux +uy +2u = e2y cos (x y)2e2y sin (x y)e2y cos (x y)+
2e2y sin (x y) = 0 and u(x, 0) = sin x.
2.5 (a) The general solution to this equation is u(x) = C where C is an
arbitrary constant.
(b) The general solution is u(x, y) = f (y) where f is an arbitrary diferentiable function of y.
2.6 (a) The general solution to this equation is u(x) = C1 x + C2 where
C1 and C2 are arbitrary constants.
(b)
R We have uy = f (y) where f is an arbitrary function of y. Hence, u(x, y) =
f (y)dy + g(x).
2.7 Let v(x, y) = y + 2x. Then
ux
uxx
uy
uyy
uxy
Hence,
uxx 4uxy + 4uyy =4fvv (v) + 4gv (v) + 4xgvv (v)
8fvv (v) 4gv (v) 8xgvv (v)
+4fvv (v) + 4xgvv (v) = 0.
2.8 utt = c2 uxx .
2.9 Let v = x + p(u)t. Using the chain rule we find
ut = fv vt = fv (p(u) + pu ut t).
Thus
(1 tfv pu )ut = fv p.
If 1 tfv pu 0 on any tinterval I then fv p 0 on I which implies that
fv 0 or p 0 on I. But either condition will imply that tfv pu 0 and
236
fv p
.
1 tfv pu
Likewise,
ux = fv (1 + pu ux t)
or
ux =
fv
.
1 tfv pu
Note that we have used the boundary conditions u(0, t) = u(L, t) = 0 and
the fact that u2x (x, t) 0 for all x [0, L].
2.12 (a) This can be done by plugging in the equations.
(b) Plug in.
(c) We have sup{|un (x, 0) 1| : x R} = n1 sup{| sin nx| : x R} = n1 .
(d) We have sup{|un (x, t) 1| : x R} =
en t
.
n
en t
n
= .
237
=cuv cuw
=c2 uvv 2c2 uwv + c2 uww
=uv + uw
=uvv + 2uvw + uww
Section 3
2
3.1 y = 21 (1 et ).
3.2 y(t) =
3t1
9
+ e2t + Ce3t .
3 cos t
t
1
(3 sin (3t)
13
+ Ct .
238
b t
e + Ceat
a
Thus,
lim y(t) = 0.
t2
4
3t + 12 +
1
.
12t2
Section 4
4.1 y(t) =
3 t2
e
2
+C
t2
4.2 y(t) = Ce 2 2t .
4.3 y(t) = Ct2 + 4.
31
239
4.4 y(t) =
2Ce4t
.
1+Ce4t
4.5 y(t) =
p
5 4 cos (2t).
p
4.6 y(t) = (2 cos t + 4).
4.7 y(t) = e1t 1.
4.8 y(t) =
2
.
4t2 +1
4.10 y(t) =
3et
.
3+et2
t2
2
f (y)dy.
= 2.
Section 5
5.1 (a) Linear (b) Quasi-linear, nonlinear (c) Nonlinear (d) Semi-linear, nonlinear.
5.2 Let w = 2x y. Then ux + 2uy u = ex f (w) + 2ex fw (w) 2ex fw (w)
ex f (w) = 0.
3 1
1 1
240
Section 6.1
6.1.1 19.
6.1.2 15.
6.1.3 ~u ~v =
6.1.4 63 .
6.1.5 52 .
1
2
and ~u w
~ = 12 .
241
6.1.6 (a) Neither (b) Orthogonal (c) Orthogonal (d) Parallel.
6.1.7 P QR is a right triangle at Q.
6.1.8 ~u =< 13 , 13 , 13 > .
6.1.9 45 .
6.1.10 Comp~a~b =
9
7
27 54
18
, , 49
49 49
>.
y
z
~i x sin
z
y
z
~j +
xy
z2
sin
y
z
~k.
6.2.3
The level surfaces are spheres centered at (2, 3, 5) and with radius
C, C 0.
6.2.4
12
.
5
6.2.5
1 (3x2
10
Section 7
2y ~
j
x2 +y 2
+ ez~k.
242
1
.
csc (yex )x
2
7.5 u(x, y) =
1
.
(x4y)2 +1y
7.6 u(x, y) =
1
.
2y
ex2 +e 1 y
Section 8
8.1 u(x, t) = sin (x 3t).
c
8.5 u(x, t) =
e3t
.
1+(x2t)2
8.6 u(x, t) = e3t (x t)2 + 19 13 t 91 .
243
8.7 Using the chain rule we find wt = ut et + uet and wx = ux et . Substituting these equations into the original equation we find
wt et u + cwx et + u = 0
or
wt + cwx = 0
8.8 (a) w(x, t) is a solution to the equation follows from the principle of
superposition. Moreover, w(x, 0) = u(x, 0) v(x, 0) = f (x) g(x).
(b) w(x, t) = f (x ct) g(x ct).
(c) From (b) we see that
sup{|u(x, t) v(x, t)|} = sup{|f (x) g(x)|}.
x,t
Thus, small changes in the initial data produces small changes in the solution. Hence, the problem is a well-posed problem.
8.9
u(x, t) =
8.10 u(x, t) = sin
2x3t
2
g t
x
c
e c x
if x < ct
if x ct.
.
Section 9
9.1 u = y+f (y ln (y + u)x) where f is an arbitrary differentiable function.
9.2 u =
f (x+y+u)
xy
244
9.5 x2 + y 2 + u2 = f
y
u
u2
2
Section 10
10.1 u(x, y) =
1xy
,
x+y
x + y 6= 0.
y
x
.
1
.
sec (xay)y
10.8 u(x, y) = h y
(x2
2
1
2
ex1 .
245
10.14 u(x, y) = 1 + 2e 2 e
(4xy)2
2
246
Section 11
11.1 (a) Hyperbolic (b) Parabolic (c) Elliptic.
11.2 (a) Ellitpic (b) Parabolic (c) Hyperbolic.
11.3 The PDE is of hyperbolic type if 4y 2 (x2 + x + 1) > 0. This is true for
all y 6= 0. Graphically, this is the xyplane with the xaxis removed,
The PDE is of parabolic type if 4y 2 (x2 + x + 1) = 0. Since x2 + x + 1 > 0
for all x R, we must have y = 0. Graphically, this is xaxis.
The PDE is of elliptic type if 4y 2 (x2 + x + 1) < 0 which can not happen.
247
11.4 We have
ux (x, t)
uxx (x, t)
ut (x, t)
utt (x, t)
=
=
=
=
sin x sin t,
cos x sin t,
cos x cos t,
cos x sin t.
Thus,
uxx (x, t)
u(x, 0)
ut (x, 0)
ux (0, t)
=
=
=
=
Plugging these expressions into the equation we find uxx + uyy = 0. Similar
argument holds for the second part of the problem.
11.7 Multiplying the equation by u and integrating, we obtain
Z L
Z L
2
u (x)dx =
uuxx (x)dx
0
0
Z L
=[u(L)ux (L) u(0)ux (0)]
u2x (x)dx
0
Z L
2
2
2
= kL u(L) + k0 u(0) +
ux (x)dx
0
248
For > 0, because k0 , kL > 0, the right-hand side is nonpositive and the
left-hand side is nonnegative. Therefore, both sides must be zero, and there
can be no solution other than u 0, which is the trivial solution.
11.8 Substitute u(x, y) = f (x)g(y) into the left side of the equation to obtain
f (x)g(y)(f (x)g(y))xy = f (x)g(y)f 0 (x)g 0 (y). Now, substitute the same thing
into the right side to obtain (f (x)g(y))x (f (x)g(y))y = f 0 (x)g(y)f (x)g 0 (y) =
f (x)g(y)f 0 (x)g 0 (y). So the sides are equal, which means f (x)g(y) is a solution.
11.9 We have
(un )xx = n2 sin nx sinh ny and (un )yy = n2 sin nx sinh ny
Hence, un = 0.
11.10 u(x, y) =
x2 y 2
4
f (x)dx.
249
By substitutition we see that c2 uxx = utt . Moreover,
Z
1
1 x
u(x, 0) = (f (x) + f (x)) +
g(s)ds = f (x)
2
2c x
and
ut (x, 0) = g(x).
11.13 (a) 1 + 4x2 y > 0, (b) 1 + 4x2 y = 0, (c) 1 + 4x2 y < 0.
11.14 u(x, y) = f (y 3x) + g(x + y).
11.15 u(x, y) = f (y 3x) + g(x + y) =
10x2 +y 2 7xy+6
.
6
Section 12
12.1 Let z(x, t) = v(x, t) + w(x, t). Then we have
c2 zxx =c2 vxx + c2 wxx
=vtt + vtt
=ztt .
12.2 Indeed we have c2 uxx (x, t) = 0 = utt (x, t).
12.3 u(x, t) = 0.
12.4 u(x, t) = 12 (cos (x 3t) + cos (x + 3t)).
h
i
1
1
.
12.5 u(x, t) = 21 1+(x+t)
2 + 1+(xt)2
12.6 u(x, t) = 1 +
12.7
1
[sin (2x
8
1
4c
250
Z
ut utt dx +
c2 ux uxt dx
0
L
2
0
2
(ut )2 dx.
The right-hand side is nonpositive, so the energy either decreases or is constant. The latter case can occur only if ut (x, t) is identically zero, which
means that the string is at rest.
12.14 (a) By the chain rule we have ut (x, t) = cR0 (x ct) and utt (x, t) =
c2 R00 (x ct). Likewise, ux (x, t) = R0 (x ct) and uxx = R00 (x ct). Thus,
utt = c2 uxx .
251
(b) We have
1
2
L
2
(ut ) dx =
0
c2 0
[R (x ct)]2 dx =
2
c2
(ux )2 dx.
2
Section 13
13.1 Let z(x, t) = u(x, t) + v(x, t). Then we have
kzxx =kuxx + kvxx
=ut + vt
=zt .
13.2 Indeed we have kuxx (x, t) = 0 = ut (x, t).
13.3 u(x, t) = T0 +
TL T0
x.
L
13.4 Let u be the solution to (13.1) that satisfies u(0, t) = u(L, t) = 0. Let
w(x, t) be the time independent solution to (13.1) that satisfies w(0, t) = T0
and w(L, t) = TL . That is, w(x, t) = T0 + TLLT0 x. From Exercise 13.1,
the function u(x, t) = u(x, t) + w(x, t) is a solution to (13.1) that satisfies u(0, t) = T0 and u(L, t) = TL .
13.5 u(x, t) = 0.
13.6 Substituting u(x, t) = X(x)T (t) into (13.1) we obtain
k
X 00
T0
= .
X
T
k XX = and
T0
T
= .
252
the ODE
X 00 X = 0 whose general
13.7 (a) Letting = k > 0 we obtain
L
+ BeL = 0.
A + B = 0. Likewise,
the
condition
u(L,
t)
=
0
implies
Ae
Hence, A(eL eL ) = 0.
(c) If A = 0 then B = 0 and u(x, t) is the trivial solution which contradicts
the assumption that u is non-trivial.Hence, we
must have
A 6= 0.
L
2L
L
= e
or e
= 1. This equa(d) Using (b) and (c) we obtain
e
>
0.
Hence,
we
must
have
< 0 so that
tion is impossible
since
2L
RL
0
f (x)dx + (1 + 4L)t.
13.13 v(x) = x + 2.
253
13.14 (a) v(x) = TL x.
(b) v(x) = T.
(c) v(x) = x + T.
RL
13.15 (a) E(t) = 0 u(x, t)dx.
(b) We integrate the equation in x from 0 to L :
Z
ut (x, t)dx =
0
u(x, t)dx.
E(t) =
0
We have
dE
=
dt
ut (x, t)dx =
ux |L0
xdx = (7 ) +
+
0
L2
.
2
Hence,
Z
E(t) =
0
L2
f (x)dx + (7 ) +
t.
2
(b) The steady solution (equilibrium) is possible if the right-hand side vanishes:
L2
(7 ) +
=0
2
2
x3
+ C1 x + C2
6
254
Section 14
14.1 (a) For all 0 x < 1 we have limn fn (x) = limn xn = 0. Also,
limn fn (1) = 1. Hence, the sequence {fn }
n=1 converges pointwise to f.
(b) Suppose the contrary. Let = 21 . Then there exists a positive integer N
such that for all n N we have
|fn (x) f (x)| <
1
2
1
2
for all x [0, 1]. Choose (0.5) N < x < 1. Then |fN (x)f (x)| = xN > 0.5 =
which is a contradiction. Hence, the given sequence does not converge uniformly.
14.2 For every real number x, we have
nx + x2
x
x2
=
lim
+
lim
=0
n
n n
n n2
n2
Thus, {fn }
n=1 converges pointwise to the zero function on R.
255
14.3 For every real number x, we have
1
1
fn (x)
.
n+1
n+1
Moreover,
1
= 0.
lim
n
n+1
Applying the squeeze rule for sequences, we obtain
lim fn (x) = 0
Therefore, {fn }
n=1 is not pointwise convergent on [0, 1].
14.5 For 2 x < 0 and 0 < x
we have
x
<
< .
n
n
n
N
256
Thus, the given sequence converges uniformly (and pointwise) to the function
f (x) = x.
(b) Since limn fn0 (x) = 1 for all x [0, 1), the sequence {fn0 }
n=1 converges
pointwise to f 0 (x) = 1. However, the convergence is not uniform. To see
this, let = 12 and suppose that the convergence is uniform. Then there is a
positive integer N such that for n N we have
1
|1 xn1 1| = |x|n1 < .
2
In particular, if we let n = N + 1 we must have xN < 21 for all x [0, 1).
1
But x = 12 N [0, 1) and xN = 12 which contradicts xN < 12 . Hence, the
convergence is not uniform.
14.7 (a) The pointwise limit is
0 if 0 x < 1
1
if x = 1
f (x) =
2
1 if 1 < x 2
(b) The convergence cannot be uniform because if it were f would have to
be continuous.
14.8 (a) Let > 0 be given. Note that
2 cos x sin2 x
1
3 .
|fn (x) | =
2
2
2(2n + sin x) 4n
3
= 0 we can find a positive integer N such that if n N
Since limn 4n
3
then 4n < . Thus, for n N and all x R we have
1
3
|fn (x) |
< .
2
4n
This shows that fn
(b) We have
Z
lim
1
2
Z
fn xdx =
Z
lim fn xdx =
2 n
1
5
dx = .
2
2
257
14.9 We have proved earlier that this sequence converges pointwise to the
discontinuous function
0 if 2 x < 0 and 0 < x 2
f (x) =
1
if x = 0
Therefore, uniform convergence cannot occur for this given sequence.
14.10 (a) Using the squeeze rule we find
lim sup{|fn (x)| : 2 x 5} = 0.
Thus, {fn }
n=1 converges uniformly to the zero function.
(b) We have
Z 5
Z 5
0dx = 0.
fn (x)dx =
lim
n
Section 15.
15.1 (a) We have (f g)(x + T ) = f (x + T )g(x + T ) = f (x)g(x) = (f g)(x).
(b) We have (c1 f +c2 g)(x+T ) = c1 f (x+T )+c2 g(x+T ) = c1 f (x)+c2 g(x) =
(c1 f + c2 g)(x).
15.2 (a) For n 6= m we have
Z
m
n
1 L
(m n)
(m + n)
sin
x sin
x dx =
x cos
x dx
cos
L
L
2 L
L
L
L
1
L
(m + n)
=
sin
x
2 (m + n)
L
L
L
(m n)
sin
x
(m n)
L
L
=0
where we used the trigonometric identiy
1
sin a sin b = [ cos (a + b) + cos (a b)].
2
258
Z
1
bn =
f (x) sin nxdx
Z
Z 0
=
sin nxdx +
sin nxdx
2
= [1 (1)n ]
n
15.5 f (x) = 61 +
15.6 f (x) =
15.7 f (x) =
2
n=1 n
4
n
n=1 (n)2 (1)
cos
4
n=1 (n)2 [1
n
2
cos (nx).
(1)n sin
(1)n ] cos
n
x
2
.
nx
2
259
15.8 Since the sided limits at the point of discontinuity x = 0 do not exist,
the function is not piecewise continuous in [1, 1].
15.9 Define the function
Z
L+a
g(a) =
f (x)dx.
L+a
sin
cos
cos
+
1
sin
15.10 (i) f (x) = 10
n=1
3
n
3
3
n
3
(ii) Using the theorem discussed in class, because this function and its derivative are piecewise continuous, the Fourier series will converge to the function
at each point of continuity. At any point of discontinuity, the Fourier series
will converge to the average of the left and right limits.
(iii)
.
15.11 (a) a0 = 2, an = bn = 0 for n N.
(b) a0 = 4, an = 0, b1 = 1, and bn = 0.
1
(c) a0 = 1, an = 0, bn = n
[(1)n 1], n N.
2L
(d) a0 = an = 0, bn = n
(1)n+1 , n N.
2nx
3
260
15.12 1
15.13 an = 0 for all n N.
15.14
f (0 )+f (0+ )
2
+
2
(b)
n=1 2n1 = 4 .
Section 16
16.1 f (x) = 0.
16.2
16.3
= 0.
n=1
sin (2n1)x
.
2n1
261
16.4
16.5 f (x) =
16.6 f (x) =
16.7 f (x) =
2
n=1 n2 [2 cos (n/2)
2
n
n=1 n2 [(1)
2
n=1 n [1
1] cos nx.
262
16.8 f (x) =
n=1 n
1+(1)n
n2 1
2
x
L
n=1
sin nx.
4[(1)n e2 1]
4+n2 2
cos (nx).
then bn = 0 if n 6= 2 and b2 = 1.
Z
n
2 L
2
bn =
sin
x dx =
[1 (1)n ].
L 0
L
n
(c) If f (x) = cos L x then
Z
2 L
cos
x sin
x dx = 0
b1 =
L 0
L
L
and for n 6= 1 we have
Z
n
2 L
bn =
cos
x sin
x dx
L 0
L
L
Z
x
i
1 2 L h x
=
sin
(1 + n) sin
(1 n) dx
2L 0
L
L
L
x
x
1
L
L
=
cos
(1 + n) +
cos
(1 n)
L
(1 + n)
L
(1 n)
L
0
2n
= 2
[1 + (1)n ].
(n 1)
16.11 (a) a0 = 10 and a1 = 1, and an = 0 for n 6= 1.
2L
n
(b) a0 = L and an = (n)
n N.
2 [(1) 1],
n
2
(c) a0 = 1 and an = n sin 2 , n N.
16.12 By definition of Fourier sine coefficients,
Z
n
2 L
bn =
f (x) sin
x dx
L 0
L
The symmetry around x = L2 can be written as
L
L
f
+x =f
x
2
2
263
for all x R. To use this symmetry it is convenient to make the change of
variable x L2 = u in the above integral to obtain
L
2
Z
bn =
f
L
2
L
n L
+ u sin
+ u du.
2
L 2
L
nu
Since f L2 + u is even in u and for n even sin n
+
u
=
sin
is
L
2
L
odd in u, the integrand of the above integral is odd in u for n even. Since
the intergral is from L2 to L2 we must have b2n = 0 for n = 0, 1, 2,
16.13 By definition of Fourier cosine coefficients,
2
an =
L
Z
0
n
f (x) cos
x dx
L
L
2
can be written as
L
y
2
= f
L
+y
2
L
2
an =
f
L
2
L
n L
+ y cos
+ y dy.
2
L 2
L
Since f L2 + y is odd in y and for n even cos n
is
+ y = cos ny
L
2
L
even in y, the integrand of the above integral is odd in y for n even. Since
the intergral is from L2 to L2 we must have a2n = 0 for all n = 0, 1, 2, .
16.14 sin
16.15 (a)
x
L
1+(1)n
n=2 n2 1
cos
nx
L
.
264
R2
(b) a0 = 22 0 f (x)dx = 3.
(c) We have
2
an =
2
Z
=
f (x) cos
0
1
nx
nx
dx
Z 2
nx
2 cos
dx +
dx
2
2
0
1
nx 1
nx 2
2
2
sin
+2
sin
=
n
2 0
n
2 1
n
2
=
sin
.
n
2
cos
nx
2
is odd in2 x 2.
n
nx
3 X
2
f (x) = +
sin
cos
.
2 n=1
n
2
2
Section 17
17.1 We look for a solution of the form u(x, y) = X(x)Y (y). Substituting in
the given equation, we obtain
X 00 Y + XY 00 + XY = 0.
265
Assuming X(x)Y (y) is nonzero, dividing for X(x)Y (y) and subtract both
00 (x)
sides for XX(x)
, we find:
X 00 (x)
Y 00 (y)
=
+ .
X(x)
Y (y)
The left hand side is a function of x while the right hand side is a function
of y. This says that they must equal to a constant. That is,
Y 00 (y)
X 00 (x)
=
+ = .
X(x)
Y (y)
X(x) =Ae x + Be
Y (y) =Cy + D
X(x) =Ae
+ Be
p
p
Y (y) =C cos ( )y + D sin ( )y
266
X(x) =Ae
Y (y) =Ce
+ Be
()y
+ De
()y
17.2 Lets assume that the solution can be written in the form u(x, t) =
X(x)T (t). Substituting into the heat equation we obtain
X 00
T0
=
.
X
kT
Since X only depends on x and T only depends on t, we must have that
there is a constant such that
X 00
X
= and
T0
kT
= .
267
X(x) =Ae
T (t) =Ce
+ Be
+ De
u(x, t) = k1 e
(x+t)
+ k2 e
(xt)
+ k3 e
If = 0 then
X(x) =Ax + B
T (t) =Ct + D
(x+t)
+ k4 e
(xt)
268
Section 18
18.1 u(x, t) = sin
x
2
2 k t
e 4 + 3 sin
18.2 u(x, t) =
8d
3
18.3 u(x, t) =
18.4 u(x, t) =
1
n=1 (2n1)3
n=1
n=1
Cn sin
sin
1
(4n2 1)
n
x
L
5
x
2
252 k t
e 4 .
(2n1)
x
L
cos
2n
x
L
k(2n1)2 2
t
L2
k 4n2 2 t
L2
e
.
kn2 2 t
e L2 where
4
n = 2, 6, 10,
n
0
n = 4, 8, 12,
Cn =
6
n is odd.
n
269
9
x
L
1
2
81k2 t
e L2 .
n=1 Cn cos
Cn =
n = 1, 5, 9,
n = 3, 7, 11,
n is even
2
n
3
x
L
8
x
L
kn2 2 t
e L2 where
2
n
n
x
L
0
9k 2
t
L2
64k2 t
e L2 .
n n2 2
1+ 2 t
L
u(x, t) =
an cos
x e
.
L
n=0
As t , e
2 2
1+ n 2 t
L
E (t) =2
270
(b) Lets assume that the solution can be written in the form u(x, t) =
X(x)T (t). Substituting into the heat equation we obtain
X 00
T0
= .
X
T
Since X only depends on x and T only depends on t, we must have that
there is a constant such that
X 00
X
= and
T0
T
= .
= and
T0
kT
= .
271
As far as the boundary conditions, we have
u(0, t) = 0 = X(0)T (t) = X(0) = 0
and
u(L, t) = 0 = X(L)T (t) = X(L) = 0.
Note that T is not the zero function for otherwise u 0 and this contradicts
our assumption that u is the non-trivial solution.
Next, we consider the three cases of the sign of .
Case 1: = 0
In this case, X 00 = 0. Solving this equation we find X(x) = ax + b. Since
X(0) = 0 we find b = 0. Since X(L) = 0 we find a = 0. Hence, X 0 and
u(x, t) 0. That is, u is the trivial solution.
Case 2: > 0
n
x, n N.
L
n2 2
kt
L2
, n N.
n n2 2
x e L2 kt , n N
L
272
Cn sin
n=1
n n2 2
x e L2 kt .
L
(13.4)
To determine
the unknown constants Cn we use the initial condition u(x, 0) =
6 sin 9x
in
(13.4)
to obtain
L
X
n
9x
Cn sin
=
x .
6 sin
L
L
n=1
By equating coefficients we find C9 = 6 and Cn = 0 if n 6= 9. Hence, the
solution to the problem is given by
9x 8122 kt
u(x, t) = 6 sin
e L .
L
(b) Similar to (a), we find
2 kt
9 2 kt
3
2
u(x, t) = 3 sin
x e L sin
x e L2
L
L
18.13 u(x, t) = cos
(b) u(x, t) = 5.
x
L
2 kt
e L2 + 4 cos
5x
L
252 kt
e L2 .
Section 19
P
n
n
y
sinh
x
where
B
sin
n
n=1
b
b
Z b
n h
n i1
2
Bn =
f2 (y) sin
y dy sinh
a
.
b 0
b
b
19.1 u(x, y) =
n
n
B
sin
x
sinh
(y
b)
where
n
n=1
a
a
Z a
n
n
2
Bn =
g1 (x) sin
x dx [sinh b ]1 .
a 0
a
a
19.2 u(x, y) =
273
3
sinh
sin x sinh y.
h
n
n i
X
n
An cosh
y + Bn sinh
y sin
x.
L
L
L
n=1
where
2
An =
L
Z
0
and
2
Bn =
L
Z
0
n
(f1 (x) + f2 (x)) sin
xdx
L
1
nH
cosh
2L
n
(f2 (x) f1 (x)) sin
xdx
L
1
nH
sinh
2L
nan (x + iy)n1 .
n=0
n=0
nan (x + iy)n1 .
n=0
274
r cos and y = r sin where r = (x2 + y 2 ) 2 and tan = xy . Using the chain
rule we obtain
ux =ur rx + u x = cos ur
sin
u
r
uxx =uxr rx + ux x
sin
sin
= cos urr + 2 u
ur cos
r
r
sin
sin
cos
u
u
+ sin ur + cos ur
r
r
r
cos
uy =ur ry + u y = sin ur +
u
r
uyy =uyr ry + uy y
cos
cos
= sin urr 2 u +
ur sin
r
r
sin
cos
cos
+ cos ur + sin ur
u +
u
r
r
r
Substituting these equations into (21.1) we obtain the dersired equation.
19.8 u(x, y) = u1 (x, y) + u2 (x, y) + u3 (x, y) + u4 (x, y) where
u1 (x, y) = 0
#
n
(1)n
n
2
sin
x
sinh
y
u2 (x, y) =
n sinh 3n
2
2
2
n=1
1
4(x 2)
4
sinh
u3 (x, y) =
sin
y
3
3
sinh 8
3
"
n
n
X
14(1 (1)n )
u4 (x, y) =
sin
y sinh
x .
3
3
n sinh 2n
3
n=1
19.9
u(x, y) =
4
sinh
L
2H
sinh
x
2H
sinh
(x L)
2H
cos
y
.
2H
275
19.10 u(x, t) =
A0
2
n x
A
e
cos n y where
n
n=1
Z
2 H
f (y)dy
A0 =
H 0
Z
2 H
n
An =
ydy.
f (y) cos
H 0
H
19.11
x
20
sin
L
cosh L H + sinh L H
L
3x
5
sin
3
3
3
L
cosh
H
+
sinh
H
L
L
L
u(x, y) =
1 2
x
2
2 cosh 3y sin 3x
3 sinh 6
Section 20
20.1 u(r, ) = 3r5 sin 5.
20.2 u(r, ) =
n
n=1 r
1(1)n
n2
cos n +
sin n
n
X
1
rn
u(r, ) =
f ()d +
f () [cos n cos n + sin n sin n] d
n
2 0
a
0
n=1
"
#
Z 2
X
1
r n
=
f () 1 + 2
cos n( ) d.
2 0
a
n=1
276
20.5 (a) We have eit = cos t + i sin t and eit = cos t i sin t. The result
follows by adding these two equalities and dividing by 2.
(b) This follows from the fact that
1
cos n( ) = (ein() + ein() ).
2
q
(c) We have |q1 | = ar cos ( )2 + sin ( )2 = ar < 1 since 0 < r < a. A
similar argument shows that |q2 | < 1.
20.6 (a) The first sum is a convergent geometric series with ratio q1 and
sum
X
r n in() ar ei()
e
=
a
1 q1
n=1
=
rei()
a rei()
X
r n
n=1
cos n( ) =1 +
rei()
a rei()
rei()
a rei()
r
r
+ i()
=1 + i()
ae
r ae
r
r
=1 +
a cos ( ) r ai sin ( )
r
+
a cos ( ) r + ai sin ( )
r[a cos ( ) r + ai sin ( )]
=1 +
a2 + 2ar cos ( ) + r2
r[a cos ( ) r ai sin ( )]
+
a2 2ar cos ( ) + r2
a2 r 2
= 2
.
a 2ar cos ( ) + r2
+
277
20.7 We have
1
u(r, ) =
2
"
f () 1 + 2
X
r n
n=1
2
#
cos n( ) d
1
a r2
=
f () 2
d
2 0
a 2ar cos ( ) + r2
Z
f ()
a2 r2 2
=
d.
2
2
a 2ar cos ( ) + r2
0
Z
20.8 u(r, ) = 2
n+1 n sin n
r n .
n=1 (1)
20.9 (a) Differentiating u(r, t) = R(r)T (t) with respect to r and t we find
utt = RT 00 and ur = R0 T and urr = R00 T.
Substituting these into the given PDE we find
1 0
00
2
00
RT = c R T + R T
r
Dividing both sides by c2 RT we find
R00 1 R0
1 T 00
=
+
.
c2 T
R
rR
Since the RHS of the above equation depends on r only, and the LHS depends
on t only, they must equal to a constant .
(b) The given boundary conditions imply
u(a, t) = 0 = R(a)T (t) = R(a) = 0
u(r, 0) = f (r) = R(r)T (0)
ut (r, 0) = g(r) = R(r)T 0 (0).
If = 0 then R00 + 1r R0 = 0 and this implies R(r) = C ln r. Using the condition
R(a) = 0 we find C = 0 so that R(r) = 0 and hence u 0. If > 0 then
T 00 c2 T = 0. This equation has the solution
278
The condition u(r, 0) = f (r) implies that A = f (r) which is not possible.
Hence, < 0
20.10 (a) Follows from the figure and the definitions of trigonometric functions in a right triangle.
(b) The result follows from equation (20.1).
20.11 By the maximum principle we have
min u(x, y) u(x, y) max u(x, y), (x, y)
(x,y)
(x,y)
(x,y)
(x,y)
(x,y)
(x,y)
(x,y)
20.14 We have
2 n
1 n
1 2 n
(r
cos
(n))
+
(r
cos
(n))
+
(r cos (n))
r2
r r
r2 2
=n(n 1)rn2 cos (n) + nrn2 cos (n) rn2 n2 cos (n) = 0
279
1
2
20.15 u(r, ) =
r2
2a2
cos 2.
20.16 u(r, ) = ln 2 + 4
a 3
r
cos 3.
Section 21
21.1 Convergent.
21.2 Divergent.
21.3 Convergent.
21.4
1
,
s3
21.5
1
s2
s > 3.
5s , s > 0.
2
4
s
21.8
es
,
s2
4
s2
2
,
s3
s > 0.
s > 0.
2s
21.9 e s +
n est
21.10 t
1
(es
s2
n
s
e2s ), s 6= 0.
tn1 est dt, s > 0.
5
s+7
1
s2
2
,
s2
s > 2.
21.13 3e2t , t 0.
21.14 2t + et , t 0.
21.15 2(e2t + e2t ), t 0.
21.16
2
s1
+ 5s , s > 1.
280
21.17
es
,
s3
21.18
1
2
21.19
3
,
s2 +36
21.20
s2
,
(s2)2 +9
21.21
2
(s4)3
1
s
s > 3.
s
s2 +4 2
, s > 0.
s > 0.
s > 3.
3
(s4)2
5
,
s4
s > 4.
0,
9(t2)
0t<2
,
t 2.
1 2t
e
20
21.29
et e2t
.
3
21.30
t
2
21.31
t5
.
120
21.32
1
2
14 e2t , t 0.
sin t.
et + 21 e2t .
21.33 t +
et
2
Section 22
et
.
2
281
Rt
0
f (t )H( xc )d.
s
22.12 u(x, t) = L1 Ts e c x + Ts .
2
Section 23
23.1
(1)n i
.
n
23.2 f (x) =
1
2
1
n=1 n
sin
n
2
(einx + einx ).
282
23.3 f (x) =
sinh a
23.4 f (x) =
eix eix
.
2i
23.5 f (x) =
1
2
n=
T+
P1
i inT
n= n [e
1]einx +
P
P
2
2
n inx
23.6 (a) f (x) = 3 + 1
+
n= n2 (1) e
n=1
P
2
4
n
(b) f (x) = 3 +
(1)
cos
nx.
n=1 n2
i inT
n=1 n [e
2
(1)n einx .
n2
23.7 (a)
Z
1
2
a0 =2
12
2
sin xdx = [cos cos ] = 0
2
2
1
2
an =2
1
2
bn =2
1]einx .
8n
4n2
c0 =0
4(1)n n
i( 4n2 )
4i(1)n n
cn =
.
4n2
P
(1)n n 2nix
(b) f (x) = 4
.
n= i(14n2 ) e
cn =
23.8 (a)
Z
1 2
(2 x)dx = 4
a0 =
2 2
Z
n
1 2
an =
(2 x) cos
x dx = 0
2 2
2
Z
n
1 2
4(1)n
bn =
(2 x) sin
x dx =
2 2
2
n
P
P
2(1)n+1 i ( in
2(1)n+1 i ( in
(b) f (x) = 2
e 2 x) +
e 2 x) .
n=1
n=1
n
n
283
4
23.9 an = cn + cn = 0. We have for |n| odd bn = i in
=
|n| even bn = 0.
4
n
and for
23.10 Note that for any complex number z we have z + z = 2Re(z) and
z z = 2iIm(z). Thus,
c n + c n = an
which means that an = 2Re(cn ). Likewise, we have
cn cn = ibn
That is ibn = 2iIm(cn ). Hence, bn = 2Im(cn ).
23.11 an = 2Re(cn ) =
23.12 f (x) = i
n=
1
n
1cos (nT )
.
n
i sin (2in) in x
e 2 .
2in
1 0<t<1
0 1<t<2
bn =
sin nxdx =
0
1 cos n
1 (1)n
=
.
n
n
Hence,
bn =
(d) We have c0 =
a0
2
1
2
2
n
if n is odd
if n is even
284
Section 24
24.1
f() =
24.2
2 sin if 6= 0
2
if = 0.
u
+ ic
u=0
t
u(, 0) = f().
24.3
2 u
= c2 2 u
2
t
u(, 0) = f()
ut (, 0) = g().
24.4
uyy = 2 u
u(, 0) = 0, u(, L) =
24.5
1
i
1
+i
2 sin a
.
2
.
2 + 2
24.6 We have
Z
F[e H(x)] =
ex H(x)eix dx
Z
=
e
0
x(1+i)
ex(1+i)
1
=
.
dx =
1 + i 0
1 + i
285
24.8 We have
Z
f (x )eix dx
Z
i
f (u)eiu du
=e
F[f (x )] =
=e
f()
where u = x .
24.9 We have
ix
F[e
ix
f (x)] =
ix
f (x)e
dx =
f (x)ei()x dx = f( ).
f (x)eix
eix
+ f (x)
2
2
1
= [F[f (x)eix ] + F[f (x)eix ]]
2
1
= [f ( ) + f( + )].
2
24.11 Using the definition and integration by parts we find
Z
0
f 0 (x)eix dx
F[f (x)] =
Z
ix
= f (x)e
+ (i)
f (x)eix dx
=(i)f()
where we used the fact that limx f (x) = 0.
24.12
2
(1
2
cos ).
24.13
2
(1
i
cos a).
286
x
24.14 F 1 [f()] = 12 e 2 .
1
= eax , x 0.
24.15 F 1 a+i
Section 25
25.1 u(x, t) = F 1 [u(, t)] = e
(xct)2
4
25.2
r
t 1 2 (kt+ )
4 ]
u(x, t) =
e F [e
4
r
r
x2
t
=
e
e 4(kt+/4)
4
kt + /4
r
x2
e 4kt+ et .
=
4kt +
25.3 u(x, t) =
1
4kt
R
0
(xs)2
4kt
ds.
25.4
2
u(x, t) =et F 1 [e t ]
1 x2
=et
e 4t .
4t
25.5 We have
Z
||y ix
Z
d =
0
y ix
d +
0
1
(y+ix)
=
e
y + ix
e e
ey eix d
1
(y+ix)
e
y ix
0
1
1
2y
=
+
= 2
.
y + ix y ix
x + y2
287
25.6
Z
1
u(x, y) =
f()e||y eix d
2
2y
1
= f (x) 2
2
x + y2
Z
1
2y
=
f (x)
d.
2
(x )2 + y 2
25.7 utt + ( + )
ut + u = c2 2 u.
25.8 u(x, t) = e(x3t) .
25.9 u(x, t) = e(xkt) .
25.10 u(x, t) =
1
4kt
2
2 (xs)
4kt
es
ds.
288
Index
Boundary value problem, 17
Burgers equation, 11
Cauchy data, 75
Cauchy problem, 75
Characteristic curve, 60
Characteristic direction, 59
Characteristic equation, 140
Characteristic equations, 60, 72
Characteristics, 60
Classical solution, 13
Convection (transport) equation, 11
Convolution, 183
Descriminant, 88
Differential equation, 5
Diffusion equation, 11
Diffusivity constant, 102
Directional derivative, 50
Dirichlet boundary conditions, 17
Dirichlet conditions, 104
Dot product, 41
Eigenvalue problem, 155
Elliptic, 88
Euler equation, 166
Euler-Fourier Formulas, 125
Even extension, 135
Even function, 133
Evolution equation, 177
Exponential order at infinity, 180
First derivative, 4
First order PDE, 36
Forced harmonic oscillator, 11
Fourier coefficients, 121
Fourier cosine series, 135
Fourier inversion formula, 207
Fourier law, 102
Fourier series, 109, 121
Fourier sine series, 135
Fourier transform, 207
Function series, 120
Fundamental period, 122
General solution, 13
Generalized solution, 14
Gradient, 51
Gradient vector field, 56
Harmonic function, 154
Heat equation, 89, 100
Heat source, 103
Helmholtz equation, 155
Homogeneous, 37, 88
Homogeneous linear PDE, 8
Hyperbolic, 88
Ill-posed, 19
Initial curve, 75
Initial data, 75
Initial temperature distribution, 104
Initial value conditions, 18
289
290
Initial value problem, 18
Inner product, 123
integral curve, 57
Integral surface, 13, 59
integral transforms, 177
Integrating factor, 26
Inverse Laplace transform, 182
Korteweg-Vries equation, 11
Lagranges method, 71
Laplace equation, 89, 154
Laplace transform, 178
Laplaces equation, 11
Laplacian, 154
Level curve, 54
Level surface, 54
Linear, 7, 36
linear, 88
Linear differential operator, 8
Linear operator, 8
INDEX
Order, 6
Ordinary differential equation, 5
Orthogonal, 123
Orthogonal projection, 45
Parabolic, 88
Partial differential equation, 5
Piecewise continuous, 124, 180
Piecewise smooth, 124
Pointwise convergence, 109, 120
Poisson Equation, 11
Poisson equation, 155
Projected characteristic curve, 60
Quasi-linear, 7, 36, 88
Right traveling wave, 67
Scalar projection, 46
Semi-linear, 7, 36, 88
Separable, 31
Separation of variables, 31
Method of characteristics, 59, 66, 71 Smooth functions, 8
Method of undetermined coefficients, Solution surface, 13
Specific heat, 101
222
Method of Variation of Parameters, Squeeze rule, 117
Stable solution, 18
229
stationary equation, 177
Minimal surface equation, 11
Strong solution, 13
Mixed boundary condition, 18
Superposition principle, 15
Mutually orthogonal, 123
Neumann boundary conditions, 17, 105 Thermal conductivity, 102
Thermal energy, 100
Non-homogeneous, 37, 88
Thin film equation, 11
Non-homogeneous PDE, 8
Total thermal energy, 103
Non-linear, 7, 36
Transport equation, 18
Nowhere characteristic, 77
Transport equation with decay, 68
Odd extension, 134
Transport equationin 1-D space, 66
Odd function, 133
Trigonometric series, 121
INDEX
Uniform convergence, 110, 121
Vector field, 55
Vector function, 53
Wave equation, 89, 93
wave equation, 11
Weak solution, 14
Weierstrass M-test, 121
Well-posed, 18
291