Solution PDF

Lecture Notes in Mathematics
A First Course in Quasi-Linear Partial

Differential Equations
for Physical Sciences and Engineering
Solution Manual
Marcel B. Finan
Arkansas Tech University
All
c Rights Reserved
March 12, 2018

2
Preface
This manuscript provides the complete and detailed solutions to A First
Course in Partial Differential Equations for Physical Sciences and Engineer-
ing. Distribution of this book in any form is prohibited.
Marcel B Finan
August 2009
Contents
Preface . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2
Solutions to Section 1 . . . . . . . . . . . . . . . . . . . . . . . . . 4
Solutions to Section 6.1 . . . . . . . . . . . . . . . . . . . . . . 36
Solution to Section 6.2 . . . . . . . . . . . . . . . . . . . . . . 39
3
4 CONTENTS
Solutions to Section 1
Problem 1.1
Classify the following equations as either ODE or PDE.
t2
(a) (y 000 )4 + (y 0 )2 +4
= 0.
∂u y−x
(b) ∂x
+ y ∂u
∂y
= y+x
.
(c) y 00 − 4y = 0.
Solution.
(a) ODE with dependent variable y and independent variabe x.
(b) PDE with dependent variable u and independent variabes x and y.
(c) ODE with dependent variable y and independent variabe x
Problem 1.2
Write the equation
uxx + 2uxy + uyy = 0
in the coordinates s = x, t = x − y.
Solution.
We have
ux =us sx + ut tx = us + ut
uxx =uss sx + ust tx + ust sx + utt tx = uss + 2ust + utt
uxy =uss sy + ust ty + ust sy + utt ty = −ust − utt
uy =us sy + ut ty = −ut
uyy = − ust sy − utt ty = utt .
Substituting these expressions into the given equation we find
uss = 0
Problem 1.3
Write the equation
uxx − 2uxy + 5uyy = 0
in the coordinates s = x + y, t = 2x.
SOLUTIONS TO SECTION 1 5
Solution.
We have
ux =us sx + ut tx = us + 2ut
uxx =uss sx + ust tx + 2ust sx + 2utt tx = uss + 4ust + 4utt
uxy =uss sy + ust ty + 2ust sy + 2utt ty = uss + 2ust
uy =us sy + ut ty = us
uyy =uss sy + ust ty = uss .
Substituting these expressions into the given equation we find
uss + utt = 0
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 nonhomo-
geneous:
(a) uux + x2 uyyy + sin x = 0.
2
(b) ux + ex uy = 0.
(c) utt + (sin y)uyy − et cos y = 0.
Solution.
(a) Order 3, non-linear.
(b) Order 1, linear, homogeneous.
(c) Order 2, linear, non-homogeneous
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 nonlinear 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).
Solution.
(a) Linear, homogeneous, order 3.
6 CONTENTS
(b) Linear, non-homogeneous, order 3. The inhomogeneity is − sin y.

(c) Non-linear, order 2. The non-linear term is ex uux .
(d) Non-linear, order 3. The non-linear terms are ux uxxy and ex uuy .
(e) Linear, non-homogeneous, order 2. The inhomogeneity is f (x, y, t)
Problem 1.6
Which of the following PDEs are linear?
(a) Laplace’s 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 .
Solution.
(a) Linear.
(b) Linear.
(c) Non-linear where all the terms are non-linear.
(d) Non-linear with non-linear term 6uux
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) Burger’s equation for h(x, t) :
ht + hhx = νhxx .
Solution.
(a) PDE, linear, second order, homogeneous.
(b) PDE, linear, second order, homogeneous.
(c) PDE, quasi-linear (non-linear), fourth order.
(d) ODE, linear, second order, non-homogeneous.
(e) PDE, linear, second order, non-homogeneous.
(f) PDE, quasilinear (non-linear), second order
Problem 1.8
Write down the general form of a linear second order differential equation of
a function in three variables.
Solution.
A(x, y, z)uxx + B(x, y, z)uxy + C(x, y, z)uyy + E(x, y, z)uxz + F (x, y, z)uyz +
G(x, y, z)uzz +H(x, y, z)ux +I(x, y, z)uy +J(x, y, z)uz +K(x, y, z)u = L(x, y, z)
Problem 1.9
Give the orders of the following PDEs, and classify them as linear or non-
linear. If the PDE is linear, specify whether it is homogeneous or non-
homogeneous.
(a) x2 uxxy + y 2 uyy − log (1 + y 2 )u = 0
(b) ux + u3 = 1
(c) uxxyy + ex ux = y
(d) uuxx + uyy − u = 0
(e) uxx + ut = 3u.
Solution.
(a) Order 3, linear, homogeneous.
(b) Order 1, non-linear.
(c) Order 4, linear, non-homogeneous
(d) Order 2, non-linear.
(e) Order 2, linear, homogeneous
8 CONTENTS
Problem 1.10
Consider the second-order PDE
uxx + 4uxy + 4uyy = 0.
Use the change of variables v(x, y) = y − 2x and w(x, y) = x to show that

uww = 0.
Solution.
Using the chain rule we find
ux = − 2uv + uw
uxx =4uvv − 4uvw + uww
uy =uv
uyy =uvv
uxy = − 2uvv + uvw .
Substituting these into the given PDE we find uww = 0
Problem 1.11
Write the one dimensional wave equation utt = c2 uxx in the coordinates
v = x + ct and w = x − ct.
Solution.
We have
ut =cuv − cuw
utt =c2 uvv − 2c2 uwv + c2 uww
ux =uv + uw
uxx =uvv + 2uvw + uww .
Substituting we find uvw = 0
Problem 1.12
Write the PDE
uxx + 2uxy − 3uyy = 0
in the coordinates v(x, y) = y − 3x and w(x, y) = x + y.
Solution.
We have
ux = − 3uv + uw
uxx = − 3(−3uv + uw )v + (−3uv + uw )w = 9uvv − 6uvw + uww
uxy = − 3uvv + uvw − 3uvw + uww = −3uvv − 2uvw + uww
uy =uv + uw
uyy =(uv + uw )v + (uv + uw )w = uvv + 2uvw + uww .
Substituting into the PDE we find uvw = 0
Problem 1.13
Write the PDE
aux + buy = 0, a 6= 0
in the coordinates s(x, y) = bx − ay and t(x, y) = x.
Solution.
According to the chain rule for the derivative of a composite function, we
have
ux =us sx + ut tx = bus + ut
uy =us sy + ut ty = −aus .
Substituting these into the given equation to obtain
abus + aut − abus = 0
or
aut = 0
and since a 6= 0 we obtain
ut = 0
Problem 1.14
Write the PDE
ux + uy = 1
in the coordinates s = x − y and t = x.
10 CONTENTS
Solution.
ux = us sx + ut tx = us + ut
uy = us sy + ut ty = −us .
Substituting these into the PDE to obtain ut = 1
Problem 1.15
Write the PDE
aut + bux = u, b 6= 0
in the coordinates v = ax − bt and w = x.
Solution.
We have ut = −buv and ux = auv + uw . Substituting we find uw = 1b u
Problem 2.1
Determine a and b so that u(x, y) = eax+by is a solution to the equation
uxxxx + uyyyy + 2uxxyy = 0.
Solution.
We have uxxxx = a4 eax+by , uyyyy = b4 eax+by , and uxxyy = a2 b2 eax+by . Thus,
substituting these into the equation we find
(a4 + 2a2 b2 + b4 )eax+by = 0.
Since eax+by 6= 0, we must have a4 + 2a2 b2 + b4 = 0 or (a2 + b2 ) = 0. This is

true only when a = b = 0. Thus, u(x, y) = 1
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.
Solution.
Substituting into the differential equation we find
tX 00 T − XT 0 = 0
or
X 00 T0
= .
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
12 CONTENTS
Problem 2.3
Consider the initial value problem
xux + (x + 1)yuy = 0, x, y > 1
u(1, 1) = e.
xex
Show that u(x, y) = y
is the solution to this problem.
Solution.
x x x xex
We have xux +(x+1)yuy = y (e +xe )+(x+1)y − y2 = 0 and u(1, 1) = e
Problem 2.4
Show that u(x, y) = e−2y sin (x − y) is the solution to the initial value prob-
lem
ux + uy + 2u = 0, x, y > 1
u(x, 0) = sin x.
Solution.
We have ux +uy +2u = e−2y cos (x − y)−2e−2y sin (x − y)−e−2y cos (x − y)+
2e−2y sin (x − y) = 0 and u(x, 0) = sin x
Problem 2.5
Solve each of the following differential equations:
(a) du
dx
= 0 where u = u(x).
(b) ∂u
∂x
= 0 where u = u(x, y).
Solution.
(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 function of
y
Problem 2.6
Solve each of the following differential equations:
2
(a) ddxu2 = 0 where u = u(x).
∂2u
(b) ∂x∂y = 0 where u = u(x, y).
Solution.
(a) The general solution to this equation is u(x) = C1 x + C2 where C1 and
C2 are arbitrary constants.
(b) We have uy =R f (y) where f is an arbitrary differentiable function of y.
Hence, u(x, y) = f (y)dy + g(x)
Problem 2.7
Show that u(x, y) = f (y + 2x) + xg(y + 2x), where f and g are two arbitrary
twice differentiable functions, satisfy the equation
uxx − 4uxy + 4uyy = 0.
Solution.
Let v(x, y) = y + 2x. Then
ux =2fv (v) + g(v) + 2xgv (v)
uxx =4fvv (v) + 4gv (v) + 4xgvv (v)
uy =fv (v) + xgv (v)
uyy =fvv (v) + xgvv (v)
uxy =2fvv (v) + gv (v) + 2xgvv (v).
Hence,
uxx − 4uxy + 4uyy =4fvv (v) + 4gv (v) + 4xgvv (v)
−8fvv (v) − 4gv (v) − 8xgvv (v)
+4fvv (v) + 4xgvv (v) = 0
Problem 2.8
Find the differential equation whose general solution is given by u(x, t) =
f (x−ct)+g(x+ct), where f and g are arbitrary twice differentiable functions
in one variable.
Solution.
Let v = x − ct and w = x + ct. We have
ux =fv vx + gw wx = fv + gw
uxx =fvv vx + gww wx = fvv + gww
ut =fv vt + gw wt = −cfv + cgw
utt = − cfvv vt + cgww wt = c2 fvv + c2 gww
Hence, u satisfies the wave equation utt = c2 uxx
14 CONTENTS
Problem 2.9
Let p : R → R be a differentiable function in one variable. Prove that
ut = p(u)ux
has a solution satisfying u(x, t) = f (x + p(u)t), where f is an arbitrary

differentiable function. Then find the general solution to ut = (sin u)ux .
Solution.
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 t−interval 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
this will imply that 1 = 1 − tfv pu = 0, a contradiction. Hence, we must have
1 − tfv pu 6≡ 0. In this case,
fv p
ut = .
1 − tfv pu
Likewise,
ux = fv · (1 + pu ux t)
or
fv
ux = .
1 − tfv pu
It follows that ut = p(u)ux .
If ut = (sin u)ux then p(u) = sin u so that the general solution is given by
u(x, t) = f (x + t sin u)
where f is an arbitrary differentiable function in one variable
Problem 2.10
Find the general solution to the pde
uxx + 2uxy + uyy = 0.
Hint: See Problem 1.2.

Solution.
Using Problem 1.2, we found uss = 0. Hence, u(s, t) = sf (t) + g(t) where
f and g are arbitrary differentiable functions. In terms of x and y we find
u(x, y) = xf (x − y) + g(x − y)
Problem 2.11
Let u(x, t) be a function such that uxx exists and u(0, t) = u(L, t) = 0 for all
t ∈ R. Prove that Z L
uxx (x, t)u(x, t)dx ≤ 0.
0
Solution.
Using integration by parts, we compute
Z L Z L
L
uxx (x, t)u(x, t)dx = ux (x, t)u(x, t)|x=0 − u2x (x, t)dx
0 0
Z L
=ux (L, t)u(L, t) − ux (0, t)u(0, t) − u2x (x, t)dx
0
Z L
=− u2 (x, t)dx ≤ 0
0
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]
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
sin nx
u(x, 0) = 1 + .
n
(c) Find sup{|un (x, 0) − 1| : x ∈ R}.
(d) Find sup{|un (x, t) − 1| : x ∈ R, t > 0}.
(e) Show that the problem is ill-posed.
16 CONTENTS
Solution.
(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 .
2
en t
(d) We have sup{|un (x, t) − 1| : x ∈ R} = n
.
2
en t
(e) We have limt→∞ sup{|un (x, t) − 1| : x ∈ R, t > 0} = limt→∞ n
= ∞.
Hence, the solution is unstable and thus the problem is ill-posed
Problem 2.13
Find the general solution of each of the following PDEs by means of direct
integration.
(a) ux = 3x2 + y 2 , u = u(x, y).
(b) uxy = x2 y, u = u(x, y).
(c) uxtt = e2x+3t , u = u(x, t).
Solution.
(a) u(x, y) = x3 + xy 2 + f (y), where f is an arbitrary differentiable function.
3 2
(b) u(x, y) = x 6y + F (x) + g(y), where F (x) = f (x)dx and g(y) is an
R
arbitrary differentiable function.

1 2x+3t
R R
(c) u(x, t) = 18 e + t f1 (x)dx + f2 (x)dx + g(t)
Problem 2.14
uxx + 4uxy + 4uyy = 0.
(a) Use the change of variables v(x, y) = y − 2x and w(x, y) = x to show

that uww = 0.
(b) Find the general solution to the given PDE.
Solution.
(a) Using the chain rule we find
ux = − 2uv + uw
uy =uv
uyy =uvv
uxy = − 2uvv + uvw .
Substituting these into the given PDE we find uww = 0.

(b) Solving the equation uww = 0 we find uw = f (v) and u(v, w) = wf (v) +
g(v). In terms of x and y the general solution is u(x, y) = xf (y − 2x) + g(y −
2x)
Problem 2.15
Derive the general solution to the PDE
utt = c2 uxx
by using the change of variables v = x + ct and w = x − ct.
Solution.
We have
ut =cuv − cuw
utt =c2 uvv − 2c2 uwv + c2 uww
ux =uv + uw
uxx =uvv + 2uvw + uww
Substituting we find uvw = 0 and solving Rthis equation we find uv = f (v)

and u(v, w) = F (v) + G(w) where F (v) = f (v)dv.
Finally, using the fact that v = x + ct and w = x − ct; we get d’Alembert’s
solution to the one-dimensional wave equation:
u(x, t) = F (x + ct) + G(x − ct)
where F and G are arbitrary differentiable functions

18 CONTENTS
Problem 3.1
Solve the IVP: y 0 + 2ty = t, y(0) = 0
Solution.
2
R
Since p(t) = 2t, we find µ(t) = e 2tdt = et . Multiplying the given equation
2
by et to obtain
2 0 2
et y = tet
Integrating both sides with respect to t and using substitution on the right-
hand integral to obtain
2 1 2
et y = et + C
2
2
Dividing the last equation by et to obtain
2 1
y(t) = Ce−t +
2
Since y(0) = 0, we find C = − 12 . Thus, the unique solution to the IVP is

given by
1 2
y = (1 − e−t )
2
Problem 3.2
Find the general solution: y 0 + 3y = t + e−2t
Solution.
Since p(t) = 3, the integrating factor is µ(t) = e3t . Thus, the general solution
is
y(t) = e−3t Re3t (t + e−2t )dt + Ce−3t
R
= e−3t 3t t
3t (te + e )dt + Ce
−3t
= e−3t e9 (3t − 1) + et + Ce−3t

= 3t−1
9
+ e−2t + Ce−3t
Problem 3.3
Find the general solution: y 0 + 1t y = 3 cos t, t > 0
Solution. R dt
Since p(t) = 1t , the integrating factor is µ(t) = e t = eln t = t. Using the
method of integrating factor we find
1
3t cos tdt + Ct
R
y(t) = t
= 3t (t sin t + cos t) + Ct
= 3 sin t + 3 cos
t
t
+ Ct
Problem 3.4
Find the general solution: y 0 + 2y = cos (3t).
Solution.
We have p(t) = 2 so that µ(t) = e2t . Thus,
Z
−2t
y(t) = e e2t cos (3t)dt + Ce−2t
But
e2t
Z Z
2t 2
e cos (3t)dt = sin (3t) − e2t sin (3t)dt
3 3
e2t 2 e2t
Z
2
= sin (3t) − (− cos (3t) + e2t cos (3t)dt)
3 3 3 3
e2t
Z
13
e2t cos (3t)dt = (3 sin (3t) + 2 cos (3t))
9 9
e2t
Z
e2t cos (3t)dt = (3 sin (3t) + 2 cos (3t))
13
Hence,
1
y(t) = (3 sin (3t) + 2 cos (3t)) + Ce−2t
13
Problem 3.5
Find the general solution: y 0 + (cos t)y = −3 cos t.
Solution.
Since p(t) = cos t we have µ(t) = esin t . Thus,
Z
− sin t
y(t) =e esin t (−3 cos t)dt + Ce− sin t
= − 3e− sin t esin t + Ce− sin t

=Ce− sin t − 3
20 CONTENTS
Problem 3.6
Given that the solution to the IVP ty 0 + 4y = αt2 , y(1) = − 13 exists on the
interval −∞ < t < ∞. What is the value of the constant α?
Solution.
4
Solving this equation with the integrating factor method with p(t) = t
we
find µ(t) = t4 . Thus,
Z
1 C
y = 4 t4 (αt)dt + 4
t t
α 2 C
= t + 4
6 t
Since the solution is assumed to be defined for all t, we must have C = 0.
On the other hand, since y(1) = − 13 we find α = −2
Problem 3.7
Suppose that y(t) = Ce−2t + t + 1 is the general solution to the equation
y 0 + p(t)y = g(t). Determine the functions p(t) and g(t).
Solution. R
The integrating factor is µ(t) = e2t . Thus, p(t)dt = 2t and this implies that
p(t) = 2. On the other hand, the function t + 1 is the particular solution
to the nonhomogeneous equation so that (t + 1)0 + 2(t + 1) = g(t). Hence,
g(t) = 2t + 3
Problem 3.8
Suppose that y(t) = −2e−t + 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).
Solution.
First, we find y0 : y0 = y(0) = −2 + 1 + 0 = −1. Next, we find g(t) : g(t) =
y 0 + y = (−2e−t + et + sin t)0 + (−2e−t + et + sin t) = 2e−t + et + cos t − 2e−t +
et + sin t = 2et + cos t + sin 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?
Solution. R
The integrating factor is µ(t) = e (1+cos t)dt = et+sin t . Thus, the general solu-
tion is
Z
−(t+sin t)
y(t) =e et+sin t (1 + cos t)dt + Ce−(t+sin t)
=1 + Ce−(t+sin t)
Since y(0) = 3, we find C = 2 and therefore y(t) = 1 + 2e−(t+sin t) . Finally,
lim y(t) = lim (1 + 2e− sin t e−t ) = 1

t→∞ t→∞
Problem 3.10
Solve the initial value problem ty 0 = y + t, y(1) = 7
Solution.
Rewriting the equation in the form
1
y0 − y = 1
t
we find p(t) = − 1t and µ(t) = 1t . Thus, the general solution is given by
y(t) = t ln |t| + Ct
But y(1) = 7 so that C = 7. Hence,
y(t) = t ln |t| + 7t
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 = be−λt
has the property that y → 0 as t → ∞. Hint: Consider the cases a = λ and

a 6= λ separately.
Solution.
Since p(t) = a we find µ(t) = eat . Suppose first that a = λ. Then
y 0 + ay = be−at
22 CONTENTS
and the corresponding general solution is
y(t) = bte−at + Ce−at
Thus,
limt→∞ y(t) = limt→∞ ( ebtat + eCat )
= limt→∞ aebat = 0
Now, suppose that a 6= λ then
b −λt
y(t) = e + Ce−at
a−λ
Thus,
lim y(t) = 0
t→∞
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)
Solution.
Substituting into the equation we find
u0 − u = −et , u(0) = 1
Solving this equation by the method of integrating factor with µ(t) = e−t we
find
u(t) = −tet + Cet
Since u(0) = 1, C = 1 and therefore u(t) = −tet + et . Finally, we have
y(t) = (−tet + et )−1
Problem 3.13
Solve the initial-value problem ty 0 + 2y = t2 − t + 1, y(1) = 1
2
Solution.
Rewriting the equation in the form
2 1
y0 + y = t − 1 +
t t
2
Since p(t) = t
we find µ(t) = t2 . The general solution is then given by
t2 t 1 C
y(t) = − + + 2
4 3 2 t
1 1
Since y(1) = 2
we find C = 12
. Hence,
t2 t 1 1
y(t) = − + +
4 3 2 12t2
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.
Solution.
Since p(t) = − 1t we find µ(t) = 1t . Thus,
0 R 0
1 t sin s
t
y = 0 s
1
t
y(t) = Si(t) + C
y(t) = tSi(t) + Ct
Since y(1) = 3, C = 3 − Si(1). Hence, y(t) = tSi(t) + (3 − Si(1))t

24 CONTENTS
Problem 4.1
Solve the (separable) differential equation
2 −ln y 2
y 0 = tet .
Solution.
At first, this equation may not appear separable, so we must simplify the
right hand side until it is clear what to do.
2 −ln y 2
y 0 =tet

1
t2 ln
=te · e y2
2 1
=tet · 2
y
t 2
= 2 et .
y
Separating the variables and solving the equation we find
2
y 2 y 0 =tet
Z Z
1 3 0 2
(y ) dt = tet dt
3
1 3 1 t2
y = e +C
3 2
3 2
y = et + C
3
2
Problem 4.2
t2 y − 4y
y0 = .
t+2
Solution.
Separating the variables and solving we find
y 0 t2 − 4
= =t−2
y t+2
Z Z
0
(ln |y|) dt = (t − 2)dt
t2
ln |y| = − 2t + C
2
t2
y(t) =Ce 2 −2t
Problem 4.3
ty 0 = 2(y − 4).
Solution.
y0 2
=
y−4 t
Z Z
0 2
(ln |y − 4|) dt = dt
t
ln |y − 4| = ln t2 + C
y−4
ln | 2 | =C
t
y(t) =Ct2 + 4
Problem 4.4
y 0 = 2y(2 − y).
Solution.
Separating the variables and solving (using partial fractions in the process)
26 CONTENTS
we find
y0
=2
y(2 − y)
y0 y0
+ =2
2y 2(2 − y)
Z Z Z
1 0 1 0
(ln |y|) dt − (ln |2 − y|) dt = 2dt
2 2

y
ln =4t + C
2 − y

y 4t
2 − y =Ce

2Ce4t
y(t) =
1 + Ce4t
Problem 4.5
Solve the IVP
4 sin (2t)
y0 = , y(0) = 1.
y
Solution.
yy 0 =4 sin (2t)
(y 2 )0 =8 sin (2t)
Z Z
2 0
(y ) dt = 8 sin (2t)dt
y 2 = − 4 cos (2t) + C
p
y(t) = ± C − 4 cos (2t).
Since y(0) = 1, we find C = 5 and hence

p
y(t) = 5 − 4 cos (2t)
Problem 4.6
Solve the IVP:
π
yy 0 = sin t, y( ) = −2.
2
Solution.
0
y2
Z Z
dt = sin tdt
2
y2
= − cos t + C
2
y 2 = − 2 cos t + C.
p
Since y( π2 ) = −2, we find C = 4. Thus, y(t) = ± (−2 cos t + 4). From
y( π2 ) = −2, we have
p
y(t) = − (−2 cos t + 4)
Problem 4.7
Solve the IVP:
y 0 + y + 1 = 0, y(1) = 0.
Solution.
(ln (y + 1))0 = − 1
ln (y + 1) = − t + C
y + 1 =Ce−t
y(t) = Ce−t − 1.
Since y(1) = 0, we find C = e. Thus, y(t) = e1−t − 1
Problem 4.8
Solve the IVP:
y 0 − ty 3 = 0, y(0) = 2.
28 CONTENTS
Solution.
Z Z
0 −3
y y dt = tdt
Z −2 0
y t2
dt = + C
−2 2
1 t2
− 2 = +C
2y 2
1
y2 = 2 .
−t + C
q
1 4
Since y(0) = 2, we find C = 4
. Thus, y(t) = ± −4t2 +1
. Since y(0) = 2, we
have y(t) = √ 2
−4t2 +1
Problem 4.9
Solve the IVP:
π
y 0 = 1 + y 2 , y( ) = −1.
4
Solution.
y0
=1
1 + y2
arctan y =t + C
y(t) = tan (t + C).
Since y( π4 ) = −1, we find C = π2 . Hence, y(t) = tan (t + π2 ) = − cot t
Problem 4.10
Solve the IVP:
1
y 0 = t − ty 2 , y(0) = .
2
Solution.
y0
=−t
y2 − 1
y0 y0
− = − 2t
y−1 y+1

y − 1
ln = − t2 + C
y + 1
y−1 2
=Ce−t
y+1
2
1 + Ce−t
y(t) = .
1 − Ce−t2
Since y(0) = 12 , we find C = − 31 . Thus,

2
3 − e−t
y(t) =
3 + e−t2
Problem 4.11
Solve the equation 3uy + uxy = 0 by using the substitution v = uy .
Solution.
Letting v = uy we obtain 3v + vx = 0. Solving this ODE by the method of
separation of variables we find
vx
=−3
v
ln |v(x, y)| = − 3x + f (y)
v(x, y) =f (y)e−3x .
f (y)e−3x dy = F (y)e−3x + G(x) where F (y) =

R R
Hence, u(x, y) = f (y)dy
Problem 4.12
Solve the IVP
(2y − sin y)y 0 = sin t − t, y(0) = 0.
30 CONTENTS
Solution.
Z Z
−3x 0
f (y)e (2y − sin y)y dt = (sin t − t)dt (4.1)
t2
f (y)e−3x y 2 + cos y = − cos t − + C. (4.2)
2
Since y(0) = 0, we find C = 2. Thus,
t2
y 2 + cos y + cos t + =2
2
Problem 4.13
State an initial value problem, with initial condition imposed at t0 = 2,
having implicit solution y 3 + t2 + sin y = 4.
Solution.
Differentiating both sides of the given equation we find
3y 2 y 0 + cos y + 2t = 0, y(2) = 0
Problem 4.14
Can the differential equation
dy
= x2 − xy
dx
be solved by the method of separation of variables? Explain.
Solution.
If we try to factor the right side of the ODE, we get
dy
= x(x − y).
dx
The second factor is a function of both x and y. The ODE is not separable
Problem 5.1
Classify each of the following PDE as linear, quasi-linear, semi-linear, or non-
linear.
(a) xux + yuy = sin (xy).
(b) ut + uux = 0
(c) u2x + u3 u4y = 0.
(d) (x + 3)ux + xy 2 uy = u3
Solution.
(a) Linear (b) Quasi-linear, non-linear (c) Non-linear (d) Semi-linear, non-
linear
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.
Solution.
Let w = 2x−y. Then ux +2uy −u = ex f (w)+2ex fw (w)−2ex fw (w)−ex f (w) =
0
Problem 5.3
√
Show that u(x, y) = x xy satisfies the equation
xux − yuy = u
subject to the constraint
u(y, y) = y 2 , y ≥ 0.
Solution. 1 1 3 1 √
We have xux − yuy = x 2 x y − y 21 x 2 y − 2 = x xy = u. Also, u(y, y) =
3 2 2
y2
Problem 5.4
Show that u(x, y) = cos (x2 + y 2 ) satisfies the equation
−yux + xuy = 0
subject to the constraint
u(0, y) = cos y 2 .
32 CONTENTS
Solution.
We have −yux + xuy = −2xy sin (x2 + y 2 ) + 2xy sin (x2 + y 2 ) = 0. Moreover,
u(0, y) = cos y 2
Problem 5.5
Show that u(x, y) = y − 12 (x2 − y 2 ) satisfies the equation
1 1 1
ux + uy =
x y y
subject to u(x, 1) = 21 (3 − x2 ).
Solution.
We have x1 ux + y1 uy = x1 (−x) + y1 (1 + y) = y1 . Moreover, u(x, 1) = 12 (3 − x2 )
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.
Solution.
Let v = ax + by. We have
d(ax + by)
ux =fv (v) = afv (v)
dx
d(ax + by)
uy =fv (v) = bfv (v).
dy
Hence, by substitution we find 3a − 7b = 0
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.
Solution.
By the chain rule we find
ux = us sx + ut tx = bus + ut
uy = us sy + ut ty = −aus .
Thus,
0 = aux + buy + cu = aut + cu.
This is a first order linear ODE that can be solved using the method of
separation of variables, assuming a 6= 0
Problem 5.8
Solve the partial differential equation
ux + uy = 1
by introducing the change of variables s = x − y and t = x.
Solution.
ux = us sx + ut tx = us + ut
uy = us sy + ut ty = −us .
Substituting these into the PDE to obtain ut = 1. Solving this ODE we
find u(s, t) = t + f (s) where f is an arbitrary differentiable function in one
variable. Now substituting for s and t we find u(x, y) = x + f (x − y)
Problem 5.9
Show that u(x, y) = e−4x f (2x − 3y) is a solution to the first-order PDE
3ux + 2uy + 12u = 0.
Solution.
We have
ux = − 4e−4x f (2x − 3y) + 2e−4x f 0 (2x − 3y)
uy = − 3e−4x f 0 (2x − 3y)
Thus,
3ux + 2uy + 12u = − 12e−4x f (2x − 3y) + 6e−4x f 0 (2x − 3y)
−6e−4x f 0 (2x − 3y) + 12e−4x f (2x − 3y) = 0
34 CONTENTS
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.
Solution.
We have ut = −buv and ux = auv + uw . Substituting we find uw = 1b u and
solving this equation by the method of integrating factor, we find u(v, w) =
1
f (v)e b w where f is an arbitrary differentiable function in one variable. Thus,
x
u(x, t) = f (ax − bt)e b
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.
Solution.
According to the chain rule for the derivative of a composite function, we
have
ux =us sx + ut tx = bus + ut
uy =us sy + ut ty = −aus .
Substituting these into the given equation to obtain
aut = 0
and since a 6= 0 we obtain

ut = 0.
Solving this equation, we find
u(s, t) = f (s)
where f is an arbitrary differentiable function of one variable. Now, in terms

of x and y we find
u(x, y) = f (bx − ay)
Problem 5.12
Write the equation
ut + cux + λu = f (x, t), c 6= 0
in the coordinates v = x − ct, w = x.
Solution.
Using the chain rule, we find
ut =uv vt + uw wt = −cuv
ux =uv vx + uw wx = uv + uw
Substituting these into the original equation we obtain the equation

w−v
cuw + λu(v, w) = f w,
c
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).
Solution.
ut = −cwv
and
ux = wv .
Substitution into the original PDE gives
vwv (v) = Aw(v)

36 CONTENTS
Solutions to Section 6.1

Problem 6.1.1
Find ~a · ~b where ~a =< 4, 1, 41 > and ~b =< 6, −3, −8 > .
Solution.
We have: ~a · ~b = 4(6) + 1(−3) + 14 (−8) = 19
Problem 6.1.2
Find ~a · ~b where ||~a|| = 6, ||~b|| = 5 and the angle between the two vectors is
120◦ .
Solution.
We have: ~a · ~b = ||~a|| · ||~b|| cos 2π
= 6(5) − 21 = −15

3
Problem 6.1.3
If ~u is a unit vector, find ~u · ~v and ~u · w
~ using the figure below.
Solution.
The figure shows that ||~u|| = ||~v || = ||w||
~ = 1 so the triangle is an equi-
lateral triangle. Hence, ~u · ~v = ||~u||||~v || cos 60◦ = 21 . Likewise, ~u · w
~ =
◦ 1
||~u||||~v || cos 120 = − 2
Problem 6.1.4
Find the angle between the vectors ~a =< 4, 3 > and ~b =< 2, −1 > .
Solution.
~a·~b
We have: cos θ = ||~a||||~b||
= 5
√
5 5
= √1 .
5
Hence, θ = cos−1 √1
5
≈ 63◦
Problem 6.1.5
Find the angle between the vectors ~a =< 4, −3, 1 > and ~b =< 2, 0, −1 > .
Solution.
~a·~b −1
We have: cos θ = ||~a||||~b||
= √ 7√
26 5
= √7 .
130
Hence, θ = cos √7
130
≈ 52◦
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.
Solution.
~a·~b
(a) ~a·~b = −40 6= 0 so the vectors are not orthogonal. Also, ||~a||||~b||
= √ −40
√
83 104
6=
0 or −1, the two vectors are not parallel.
(b) ~a · ~b = 0, the vectors are orhtogonal.
(c) ~a · ~b = 0, the vectors are orhtogonal.
(d) ~u = − 23 ~v , the vectors are parallel
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.
Solution.
−→ −→ −→ −→
We have: QP =< −1, −3, 2 > and QR =< 4, −2, −1 > . Since QP · QR = 0,
−→ −→
the vectors QP and QR are orthogonal so that the triangle P QR is a right
triangle at Q
Problem 6.1.8
Find a unit vector that is orthogonal to both ~i + ~j and ~i + ~k.
Solution.
Let ~u =< a, b, c > with a2 + b2 + c2 = 1. We are told that < a, b, c > · <
1, 1, 0 >= 0 and < a, b, c > · < 1, 0, 1 >= 0. These imply a + b = 0 and
a + c = 0. Hence, a2 + (−a)2 + (−a)2 = 1 or a = ± √13 . It follows that either
~u =< √13 , − √13 , − √13 > or ~u =< − √13 , √13 , √13 >
Problem 6.1.9
Find the acute angle between the lines 2x − y = 3 and 3x + y = 7.
Solution.
Since the points (1, −1) and (0, −3) are on the line 2x − y = 3, the vector
~v =< −1, −2 > is a vector on the line. Likewise, using the points (0, 7) and
38 CONTENTS
~ =< 1, −3 > is on the line 3x + y = 7. The angle between

(1, 4), the vector w
the two lines is the angle between the two vectors given by
~v · w~ 5 1
cos θ = =√ √ =√ .
||~v ||||w||
~ 5 10 2
Hence, θ = 45◦
Problem 6.1.10
Find the scalar and vector projections of the vector ~b =< 1, 2, 3 > onto
~a =< 3, 6, −2 > .
Solution.
We have
~a · ~b 9
Comp~a~b = =
||~a|| 7
~a · ~b 27 54 18
Proj~a~b = ~
a =< , ,− >
||~a||2 49 49 49
Problem 6.1.11
If ~a =< 3, 0, −1 >, find a vector ~b such that comp~a~b = 2.
Solution.
~a·~b
Suppose ~b =< x, y, z > . We are told that Comp~a~b = 2 or ||~ a||
= 2. Hence,
√ √ √
~a · ~b = 2||~a|| = 2 10. Thus, 3x − z = 2 10. Hence, ~b =< x, y, 3x − 2 10 >
where x and y are arbitrary numbers
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.
Solution.
~ =< 6 − 0, 12 − 10, 20 − 8 >=< 6, 2, 12 > . The
The displacement vector is D
work done is W = F~ · D
~ = 8(6) + (−6)(2) + (9)(12) = 144 Joules
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.
Solution.
The work done by the force is W = 30(80) cos 40◦ ≈ 1839 ft − lb = 1839 slug
Solution to Section 6.2

Problem 6.2.1
Find the gradient of the function
F (x, y, z) = exyz + sin (xy).
Solution.
We have
Fx (x, y, z) = yzexyz +y cos (xy), Fy (x, y, z) = xzexyz +x cos (xy), Fz (x, y, z) = xyexyz .
Thus,
∇F (x, y, z) = (yzexyz + y cos (xy))~i + (xzexyz + x cos (xy))~j + xyexyz~k
Problem 6.2.2
Find the gradient of the function
y
F (x, y, z) = x cos .
z
Solution.
We have
y x y xy y
Fx (x, y, z) = cos , Fy (x, y, z) = − sin , Fz (x, y, z) = 2 sin .
z z z z z
Thus, y
~ x y
~ xy y
~k
∇F (x, y, z) = cos i − sin j + 2 sin
z z z z 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 .
Solution. √
The level surfaces are spheres centered at (2, 3, −5) and with radius C, C ≥
0
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).
40 CONTENTS
Solution.
The unit vector in the direction of ~a is
1 2
~v = √ ~i + √ ~j.
5 5
We have
ux (x, y) =8x
ux (1, 1) =8
uy (x, y) =2y
uy (1, 1) =2.
Hence,
1 2 12
u~v (1, 1) = ( √ , √ ) · (8, 2) = √
5 5 5
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).
Solution.
The unit vector in the direction of ~a is
1 3
~v = − √ ~i + √ ~j.
10 10
We have
ux (x, y, z) = 2xz − yz, uy (x, y, z) = 3y 2 z 2 − xz, uz (x, y, z) = x2 + 2y 3 z − xy.
Hence,

1 3
u~v (x, y, z) = − √ , 0, √ · (2xz − yz, 3y 2 z 2 − xz, x2 + 2y 3 z − xy)
10 10
1
= √ (3x2 + 6y 3 z − 3xy − 2xz + yz)
10
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.
Solution.
The maximum rate of change is given by ||∇u(0, 2)||. But
∇u(x, y) = y 2 exy~i + (1 + xy)exy~j.
Thus,
∇u(0, 2) = 4~i + ~j
so that the maximum rate of change at (0, 2) is given by
√ √
||∇u(0, 2)|| = 42 + 11 = 17.
The maximum occurs in the direction of the unit vector

∇u(0, 2) 4 1
= √ ~i + √ ~j
||∇u(0, 2)|| 17 17
Problem 6.2.7
Find the gradient vector field for the function u(x, y, z) = ez − ln (x2 + y 2 ).
Solution.
We have
2x
ux (x, y, z) = −
x2 + y 2
2y
uy (x, y, z) = − 2
x + y2
uz (x, y, z) =ez .
Hence,
2x ~ 2y ~
∇u(x, y, z) = − 2 2
i − 2 2
j + ez~k
x +y x +y
42 CONTENTS
Problem 7.1
Solve ux + yuy = y 2 with the initial condition u(0, y) = sin y.
Solution.
dy
We have a = 1, b = y, and f = y 2 . Solving dx
= y we find y = k1 ex . Solving
k2
du
dx
= y 2 = k12 e2x we find u = 21 e2x + f (k1 ) = 21 y 2 + f (k1 ) = 21 y 2 + f (ye−x ).
Using the initial condition u(0, y) = sin y we find sin y − 21 y 2 = f (y). Hence,
u(x, y) = 21 y 2 − 21 y 2 e−2x + sin (ye−x )
Problem 7.2
Solve ux + yuy = u2 with the initial condition u(0, y) = sin y.
Solution.
dy
We have a = 1, b = y, and f = u2 . Solving dx = y we find y = k1 ex . Solving
du
dx
= u2 we find x + u1 = k2 . Thus, u(x, y) = f (ye−x 1
)−x
. Using the initial
condition u(0, y) = sin y we find f (y) = csc y. Hence, u(x, y) = csc (ye1−x )−x
Problem 7.3
Find the general solution of yux − xuy = 2xyu.
Solution.
The system of ODEs is
dy x du
=− , = 2xu.
dx y dx
Solving the first equation, we find x2 + y 2 = k1 . Solving the second equation,

2 2
we find u = k2 ex . Hence, u(x, y) = ex f (x2 + y 2 ) where f is an arbitrary
differentiable function in one variable
Problem 7.4
Find the integral surface of the IVP: xux + yuy = u, u(x, 1) = 2 + e−|x| .
Solution.
The system of ODEs is
dy y du u
= , = .
dx x dx x
Solving the first equation, we find y = k1 x. Solving the second equation, we

y
find u = k2 x. Hence, u(x, y) = xf x where f is an arbitrary differentiable
function in one variable. From the initial condition u(x, 1) = 2 + e−|x| we
1
find f (x) = x(2 + e− |x| ). Hence, the integral surface is
u(x, y) = y(2 + e−| y | )

x
Problem 7.5
1
Find the unique solution to 4ux + uy = u2 , u(x, 0) = 1+x2
.
Solution.
The system of ODEs can be written as
dx dy du
= = 2.
4 1 u
Solving the equation dx 4

= dy 1
we find x − 4y = k1 . Solving the equation
dy du 1 1
1
= u2 we find u(x, y) = f (x−4y)−y . Using the initial condition u(x, 0) = 1+x 2
2 1
we find f (x) = 1 + x . Hence, u(x, y) = (x−4y)2 +1−y
Problem 7.6
2
Find the unique solution to e2y ux + xuy = xu2 , u(x, 0) = ex .
Solution.
dx dy du
2y
= = 2.
e x xu
du
Thus, xdx = e2y dy which implies x2 − e2y = k1 . Solving the equation u2
= dy
we find y + u1 = k2 = f (x2 − e2y ). Hence,
1
u(x, y) = .
f (x2 − e2y ) − y
2
Using the initial condition u(x, 0) = ex we find f (x) = e−(x+1) . Hence,
1
u(x, y) = −x2 +e2y −1
e −y
44 CONTENTS
Problem 7.7
Find the unique solution to xux + uy = 3x − u, u(x, 0) = tan−1 x.
Solution.
dx dy du
= = .
x 1 3x − u
dy
Solving the equation dx
x
= 1
we find xe−y = k1 . On the other hand, we have
dx du d(3x − u)
= = =⇒ (3x − u)d(3x − u) = udu.
x 3x − u u
Thus, (3x − u)2 − u2 = k2 = f (xe−y ) which leads to
3 1
u(x, y) = x − f (xe−y ).
2 6x
Using the initial condition, u(x, 0) = tan−1 x we find f (x) = 9x2 −6x tan−1 x.
Hence,
3 9x2 e−2y − 6xe−y tan−1 (xe−y ) 3 3
u(x, y) = x− = x− xe−2y +e−y tan−1 (xe−y )
2 6x 2 2
Problem 7.8
Solve: xux − yuy = 0, u(x, x) = x4 .
Solution.
dy
Solving the equation dx = − xy we find xy = k1 . Since the right-hand side
is 0, u(x, y) = k2 = f (k1 ) = f (xy). But u(x, x) = x4 = f (x2 ). Hence,
f (x) = x2 , x ≥ 0, Hence,
u(x, y) = x2 y 2 , xy ≥ 0
Problem 7.9
Find the general solution of yux − 3x2 yuy = 3x2 u.
Solution.
dy
Solving the equation dx = −3x2 we find y + x3 = k1 . Solving the equation
du
u
= − dy
y
we find uy = k2 = f (k1 ) = f (y + x3 ) where f is a differentiable
function in one variable
Problem 7.10
Find u(x, y) that satisfies yux + xuy = 4xy 3 subject to the boundary condi-
tions u(x, 0) = −x4 and u(0, y) = 0.
Solution.
dy
Solving the equation dx = xy we find y 2 − x2 = k1 . On the other hand,
du = 4y 3 dy so that u(x, y) = y 4 + f (y 2 − x2 ). Since u(x, 0) = −x4 , we
have f (−x2 ) = −x4 or f (x) = −x2 for x ≤ 0. Since u(0, y) = 0 we find
f (y 2 ) = −y 4 so that f (y) = −y 2 for y ≥ 0. Hence, f (x) = −x2 for all x.
Finally,
u(x, y) = y 4 − (y 2 − x2 )2 = 2x2 y 2 − x4
46 CONTENTS
Problem 8.1
Find the solution to ut + 3ux = 0, u(x, 0) = sin x.
Solution.
dt
Solving dx = 13 we find x − 3t = k1 . Solving the equation du dx
= 0 we find
u(x, t) = k2 = f (x − 3t) where f is a differentiable function in one variable.
Since u(x, 0) = sin x, we find sin x = f (x). Hence, u(x, t) = sin (x − 3t)
Problem 8.2
Solve the equation aux + buy + cu = 0.
Solution.
dy
Solving the equation dx = ab we find bx − ay = k1 . Solving the equation du dx
=
c − ac x − ac x
− a u we find u(x, y) = k2 e = f (bx − ay)e where f is a differentiable
function in one variable
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.
Solution.
dy
Solving the equation dx = 2 we find 2x − y = k1 . Solving the equation
du
dx
= cos (y − 2x) = cos k 1 we find
u(x, y) = x cos k1 + k2 = x cos (y − 2x) + g(2x − y)

where g is a differentiable function in one variable.
Since u(0, y) = f (y), we obtain f (y) = g(−y) or g(y) = f (−y). Thus,
u(x, y) = x cos (y − 2x) + f (y − 2x)
Problem 8.4
Show that the initial value problem ut + ux = x, u(x, x) = 1 has no solution.
Solution.
dy
Solving the equation dx = 1 we find x − y = k1 . Solving the equation dudx
=x
1 2
we find u(x, y) = 2 x + f (x − y) where f is a differentiable function of one
2
variable. Since u(x, x) = 1 we find 1 = 21 x2 + f (0) or f (0) = 1 − x2 which
is impossible since f (0) is a constant. Hence, the given initial value problem
has no solution
Problem 8.5
Solve the transport equation ut + 2ux = −3u with initial condition u(x, 0) =
1
1+x2
.
Solution.
dt
Solving the equation dx = 12 we find x − 2t = k1 . Solving the equation
3
du
dx
= − 32 u we find u(x, t) = f (x − 2t)e− 2 x . Since u(x, 0) = 1+x
1
2 we find
3
e2x
f (x) = 1+x2
. Hence,
e−3t
u(x, t) =
1 + (x − 2t)2
Problem 8.6
Solve ut + ux − 3u = t with initial condition u(x, 0) = x2 .
Solution.
dt du
Solving the equation dx = 1 we find x − t = k1 . Solving the equation dt
=
3u + t by the method of integrating factor, we find
1 1 1 1
u(x, t) = − t − + k2 e3t = − t − + f (x − t)e3t .
3 9 3 9
But u(x, 0) = x2 which leads to f (x) = x2 + 19 . Hence,

3t 2 1 1 1
u(x, t) = e (x − t) + − t−
9 3 9
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 = ueλt .
Solution.
Using the chain rule we find wt = ut eλt + λueλt and wx = ux eλt . Substituting
these equations into the original equation we find
wt e−λt − λu + cwx e−λt + λu = 0
or
wt + cwx = 0
48 CONTENTS
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)
where f and g are continuously differentiable functions.
(a) Show that w(x, t) = u(x, t) − v(x, t) is the unique solution to the IVP
ut + cux = 0
u(x, 0) = f (x) − g(x)

(b) Write an explicit formula for w in terms of f and g.
(c) Use (b) to conclude that the transport problem is well-posed. That is, a
small change in the initial data leads to a small change in the solution.
Solution.
(a) w(x, t) is a solution to the equation follows from the principle of super-
position. 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 x
Thus, small changes in the initial data produces small changes in the solution.
Hence, the problem is a well-posed problem
Problem 8.9
Solve the initial boundary value problem
ut + cux = −λu, x > 0, t > 0
u(x, 0) = 0, u(0, t) = g(t), t > 0.

Solution.
dt
Solving dx = 1c we find x − ct = k1 . Solving the equation dudx
= − λc u we find
λ
u(x, t) = f (x − ct)e− c x . From the condition u(0, t) = g(t) we find f (−ct) =
g(t) or f (t) = g − ct . Thus,
x −λx
u(x, t) = g t − e c .
c
This is valid only for x < ct since g is defined on (0, ∞). Also, this expression
will not satisfy u(x, 0) = 0. So we define u(x, t) = 0 for x ≥ ct. That is, the
solution to the initial boundary value problem is
λ
g t − xc e− c x if x < ct
u(x, t) =
0 if x ≥ ct
Problem 8.10
Solve the first-order equation 2ut +3ux = 0 with the initial condition u(x, 0) =
sin x.
Solution.
dt
Solving the equation dx = 23 we find 2x − 3t = k1 . Solving the equation
du
dx
= 0 we find u(x, t) = k2 = f (2x−3t) where f is an arbitrary differentiable
function. Using the initial condition we find f (2x) = sin x or f (x) = sin x2 .

The final answer is u(x, t) = sin 2x−3t

2
Problem 8.11
Solve the PDE ux + uy = 1.
Solution.
dy
Solving the equation dx = 1 we find x − y = k1 . Solving the equation du
dx
=1
we find u(x, y) = x + f (x − y) where f is an arbirary differentiable function
in one variable
50 CONTENTS
Problem 9.1
Find the general solution of the PDE ln (y + u)ux + uy = −1.
Solution.
dx
The characteristic equations are ln (y+u) = dy
1
= −1du
. Using the second and
third fractions we find that y + u = k1 . Now, from the first and second
fractions we have lndxk1 = dy
1
so that x + k2 = y ln k1 . Hence, y ln (y + u) − x =
k2 . Hence, the general solution is given by u = −y + f (y ln (y + u) − x) where
f is an arbitrary differentiable function
Problem 9.2
Find the general solution of the PDE x(y − u)ux + y(u − x)uy = u(x − y).
Solution.
dx dy du
The characteristic equations are given by x(y−u)
= y(u−x)
= u(x−y)
. We have
dx + dy + du du
=
x(y − u) + y(u − x) + u(x − y) u(x − y)
or
d(x + y + u) du
= .
0 u(x − y)
Hence. x + y + u = k1 . On the other hand we have
dx dy du
x
+ y u
= .
−(x − y) x−y
This implies that

dx dy du
+ =−
x y u
or
ln xyu = k
f (x+y+u)
that is xyu = k2 . Hence, the general solution is given by u = xy
where
f is an arbitrary differentiable function
Problem 9.3
Find the general solution of the PDE u(u2 + xy)(xux − yuy ) = x4 .
Solution.
dy
The characteristic equations are ux(udx du
2 +xy) = − yu(u2 +xy) = x4 . From the first
and second fractions we get dx x

= − dyy
. Upon integration we find xy = k1 .
From first and third fractions we get x dx = (u3 + uxy)du or x3 dx = (u3 +
3
4 4
k1 u)du. Integration leads to x4 = u4 + k21 u2 + k2 or x4 − u4 − 2k1 u2 = k2 .
Substituting for k1 we find x4 − u4 − 2xyu2 = k2 . Hence, the general solution
is given by x4 − u4 − 2xyu2 = f (xy) where f is an arbitrary differentiable
function
Problem 9.4
Find the general solution of the PDE (y + xu)ux − (x + yu)uy = x2 − y 2 .
Solution.
dx dy du
The characteristic equations are y+xu
= − x+yu = x2 −y 2
. We have
xdx + ydy − udu du

= 2 .
xy + x2 u 2 2
− xy − y u − ux + uy 2 x − y2
Thus, xdx + ydy − udu = 0 or x2 + y 2 − u2 = k1 . On the other hand,
ydx + xdy + du du
= 2 .
y2 2 2
+ xyu − x − xyu + x − y 2 x − y2
That is, ydx + xdy + udu = 0. Hence, xy + u = k2 . The general solution is
xy + u = f (x2 + y 2 − u2 ) where f is an arbitrary differentiable function
Problem 9.5
Find the general solution of the PDE (y 2 + u2 )ux − xyuy + xu = 0.
Solution.
dy
The characteristic equations are y2dx
+u2
= −xy du
= −xu . Using the last two
dy y
fractions we find y = du
u
which leads to u = k1 . On the other hand, we have
xdx + ydy + udu du

= .
xy 2 2 2
+ xu − xy − xu 2 −xu
Thus, xdx + ydy +udu = 0 or x2 + y 2 + u2 = k2 . The general solution is
x2 + y 2 + u2 = f uy where f is an arbitrary differentiable function
Problem 9.6
Find the general solution of the PDE ut + uux = x.
52 CONTENTS
Solution.
The characteristic equations are
dx dt du
= = .
u 1 x
Solving dxu
= du
x
we find u2 − x2 = k1 . Solving dt1 = d(x+u)
x+u
we find x + u =
t t 2 2
k2 e = e f (u − x ) where f is an arbitrary differentiable function
Problem 9.7
Find the general solution of the PDE (y − u)ux + (u − x)uy = x − y.
Solution.
dx dy du
= = .
y−u u−x x−y
We have
dx + dy + du dx + dy + du
=
y−u+u−x+x−y 0
so that dx + dy + du = 0. Hence, x + y + u = k1 . Likewise,
xdx + ydy + udu xdx + ydy + udu
=
x(y − u) + y(u − x) + u(x − y) 0
so that xdx + ydy + udu = 0. Hence, x2 + y 2 + u2 = k2 . The general solution

is given by
x2 + y 2 + u2 = f (x + y + u)
where f is an arbitrary differentiable function
Problem 9.8
Solve
x(y 2 + u)ux − y(x2 + u)uy = (x2 − y 2 )u.
Solution.
dx dy du
= = 2 .
x(y 2+ u) 2
−y(x + u) (x − y 2 )u
We first note that

dx dy du dx dy du
x y u x
+ y
+ u
= = = .
y2 + u 2
−(x + u) x2 − y 2 0
Thus,
dx dy du
+ + =0
x y u
which gives xyu = k1 . Likewise, we have
xdx + ydy − du xdx + ydy − du
= .
x2 (y 2 2 2 2 2
+ u) − y (x + u) − (x − y )u 0
Thus, xdx+ydy −du = 0 and this implies that x2 +y 2 −2u = k2 . The general
solution is given by
x2 + y 2 − 2u = f (xyu)
where f is an arbitrary differentiable function
Problem 9.9
Solve √
1 − x2 ux + uy = 0.
Solution.
dx dy du
√ = = .
1−x 2 1 0
From the last fraction, we have u(x, y) = k1 . From the first two fractions, we
have y = sin−1 x + k2 = sin−1 x + f (u) where f is a differentiable function
Problem 9.10
Solve
u(x + y)ux + u(x − y)uy = x2 + y 2 .
Solution.
dx dy du
= = 2 .
u(x + y) u(x − y) x + y2
54 CONTENTS
Each of these ratio is equivalent to

ydx + xdy − udu xdx − ydy − udu
=
0 0
or
u2 1
d(xy − 2
) 2
d(x2 − y 2 − u2 )
= .
0 0
Hence,
1 2 u2
(x − y 2 − u2 ) = f (xy − )
2 2
where f is a differentiable function
Problem 10.1
Solve
(y − u)ux + (u − x)uy = x − y
with the condition u x, x1 = 0, x 6= 0, x 6= 1.

Solution.
Let
1
Γ : x0 (t) = t, y0 (t) = , u0 (t) = 0.
t
Then
dy0 dx0 1
− b(x0 , y0 , u0 )
a(x0 , y0 , u0 ) = t − 3 6= 0.
dt dt t
Hence, the problem has a unique solution.
dx dy du
= = .
y−u u−x x−y
We have
dx + dy + du dx + dy + du
=
y−u+u−x+x−y 0
so that dx + dy + du = 0. Hence, x + y + u = c1 . Likewise,
xdx + ydy + udu xdx + ydy + udu
=
x(y − u) + y(u − x) + u(x − y) 0
so that xdx + ydy + udu = 0. Hence, x2 + y 2 + u2 = c2 . The general solution
is given by
x2 + y 2 + u2 = f (x + y + u)
where f is an arbitrary differentiable function. Now, using the Cauchy data
1 1 1 1 2

u x, x = 0 we find f x + x = x + x2 = x + x −2. Hence, f (x) = x2 −2.
2
Hence, the integral surface is described by

(x + y + u)2 = x2 + y 2 + u2 + 2
and the unique solution is given by
1 − xy
u(x, y) = , x + y 6= 0
x+y
56 CONTENTS
Problem 10.2
Solve the linear equation
yux + xuy = u,
with the Cauchy data u(x, 0) = x3 , x > 0
Solution.
Let
Γ : x0 (t) = t, y0 (t) = 0, u0 (t) = t3 .
Then
dy0 dx0
− b(x0 , y0 , u0 )
a(x0 , y0 , u0 ) = −t 6= 0.
dt dt
dx dy du
= = .
y x u
Using the first two fractions we find x2 − y 2 = c1 Now, since

du dx + dy
=
u x+y
we can write u = c2 (x + y). Hence, the general solution is given by
u = (x + y)f (x2 − y 2 )
where f is an arbitrary differentiable function. Using the Cauchy data

u(x, 0) = x3 , we find f (x2 ) = x2 , that is f (x) = x. Consequently, the unique
u(x, y) = (x + y)(x2 − y 2 )
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.
Solution.
Let
Γ : x0 (t) = t, y0 (t) = −t, u0 (t) = 1.
Then
dy0 dx0
− b(x0 , y0 , u0 )
a(x0 , y0 , u0 ) = −2t(t2 + 1) 6= 0.
dt dt
dx dy du
= = 2 .
x(y 2+ u) 2
−y(x + u) (x − y 2 )u
We first note that

dx dy du dx dy du
x y u x
+ y
+ u
= = = .
y2 + u 2
−(x + u) x2 − y 2 0
Thus,
dx dy du
+ + =0
x y u
which gives xyu = c1 . Likewise, we have
xdx + ydy − du xdx + ydy − du

= .
x2 (y 2 2 2 2 2
+ u) − y (x + u) − (x − y )u 0
Thus, xdx + ydy − du = 0 and this implies that x2 + y 2 − 2u = c2 . The general

f (xyu) = x2 + y 2 − 2u
u(x, −x) = 1 we find f (x) = −2x − 2. Hence, the unique solution is given by
2xyu + x2 + y 2 − 2u + 2 = 0
Problem 10.4
Solve
xux + yuy = xe−u
with the Cauchy data u(x, x2 ) = 0, x > 0.
Solution.
Let
Γ : x0 (t) = t, y0 (t) = t2 , u0 (t) = 0.
58 CONTENTS
Then
dy0 dx0
− b(x0 , y0 , u0 )
a(x0 , y0 , u0 ) = t2 6= 0.
dt dt
dx dy du
= = −u .
x y xe
Using the first two fractions we find xy = c1 . Using the first and the last
fractions we find dx = eu du or x − eu = c2 . Hence, the general solution is
given by y
x − eu ) = f
x
u(x, x2 ) = 0, we find f (x) = x − 1. Hence, the unique integral surface is
described by
y
− + x − eu + 1 = 0
x
or y
u(x, y) = ln x + 1 −
x
Problem 10.5
Solve the initial value problem
xux + uy = 0, u(x, 0) = f (x).
Solution.
The initial curve in parametric form is
Γ : x0 (t) = t, y0 (t) = 0, u0 (t) = f (t).
Since
dy0 dx0
(t) − b(x0 (t), y0 (t), u0 (t))
a(x0 (t), y0 (t), u0 (t)) (t) = −1 6= 0
dt dt
the initial value problem has a unique solution. The characteristic equations
are
dx dy du
= = .
x 1 0
Hence xe−y = c1 and u(x, y) = c2 = g(xe−y ) where g is a differentiable
function in one variable. Using the Cauchy data u(x, 0) = f (x) we find
g(x) = f (x). Hence, the unique solution is u(x, y) = f (xe−y )
Problem 10.6
ut + aux = 0, u(x, 0) = f (x).
Solution.
The initial curve parametrization is given by
Γ : x0 (w) = w, t0 (w) = 0, u0 (w) = f (w).
Since
dt0 dx0
a(x0 (w), t0 (w), u0 (w)) (w) − b(x0 (w), t0 (w), u0 (w)) (w) = −1 6= 0
dw dw
are
dt dx du
= = .
1 a 0
Thus, x − at = c1 and u(t, x) = c2 = g(x − at) where g is a differentiable
function in one variable. Hence, the unique solution is given by u(t, x) =
f (x − at)
Problem 10.7
aux + uy = u2 , u(x, 0) = cos x
Solution.
Γ : x0 (t) = t, y0 (t) = 0, u0 (t) = cos t.
Since
dy0 dx0
a(x0 (t), y0 (t), u0 (t)) (t) − b(x0 (t), y0 (t), u0 (t)) (t) = −1 6= 0
dt dt
are
dx dy du
= = 2.
a 1 u
60 CONTENTS
Hence, x − ay = c1 and y + u1 = c2 = f (x − ay) where f is a differentiable

function in one variable. Using the Cauchy data u(x, 0) = cos x we find
f (x) = sec x. The unique solution is
1
u(x, y) =
sec (x − ay) − y
Problem 10.8
ux + xuy = u, u(1, y) = h(y).
Solution.
Γ : x0 (t) = 1, y0 (t) = t, u0 (t) = h(t).
Since
dy0 dx0
a(x0 (t), y0 (t), u0 (t)) (t) − b(x0 (t), y0 (t), u0 (t)) (t) = 1 6= 0
dt dt
are
dx dy du
= = .
1 x u
Thus, 2y − x2 = c1 and u = f (2y 2 x
− x )e . Using the Cauchy data u(1, y) =
−1 y+1
h(y), we find f (y) = e h 2 . Hence,
x2 1 x−1

u(x, y) = h y − + e
2 2
Problem 10.9
uux + uy = 0, u(x, 0) = f (x)
where f is an invertible function.

Solution.
Γ : x0 (t) = t, y0 (t) = 0, u0 (t) = f (t).
Since
dy0 dx0
dt dt
are
dx dy du
= = .
u 1 0
Thus, u(x, y) = c1 and x − uy = c2 = g(u). Using the Cauchy data u(x, 0) =
f (x) we find g(f (x)) = x and so g(x) = f −1 (x). Hence, the unique solution
is given by
u(x, y) = f (x − uy)
Problem 10.10
√
1 − x2 ux + uy = 0, u(0, y) = y.
Solution.
Γ : x0 (t) = 0, y0 (t) = t, u0 (t) = t.
Since
dy0 dx0
a(x0 (t), y0 (t), u0 (t)) (t) − b(x0 (t), y0 (t), u0 (t)) (t) = 1 6= 0
dt dt
are
dx dy du
√ = = .
1−x 2 1 0
Thus, y − arcsin x = c1 and u(x, y) = c2 = f (y − arcsin x). Using the Cauchy
data u(0, y) = y we find f (y) = y. Hence, the unique solution is given by
u(x, y) = y − arcsin x
62 CONTENTS
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?
Solution.
dy
(i) The characterisitcs are solutions to the ODE dx = 2yx
. Solving we find
2
y = Cx . Thus, the characteristics are parobolas in the plane centered at the
origin. See figure below.
(ii) The general solution is u(x, y) = f (yx−2 ), and so the solution satisfy-
ing the condition at u(1, y) = ey is
−2
u(x, y) = eyx .
(iii) In the first case, we cannot substitute x = 0 into yx−2 (the argument
of the function f, above) because x−2 is not defined at 0. Similarly, in the
second case, we’d need to find a function f so that f (0) = h(x). If h is not
constant, it is not possible to satisfy this condition for all x ∈ R
Problem 10.12
Solve the equation ux + uy = u subject to the condition u(x, 0) = cos x.
Solution.
The initial curve in R3 can be given parametrically as
Γ : x0 (t) = t, y0 (t) = 0, u0 (t) = cos t.
We have
dy0 dx0
dt dt
so by Theorem 10.1 the given Cauchy problem has a unique solution. The
characteristic equations are
dx dy du
= = .
1 1 u
Thus, x − y = c1 and u = f (x − y)ey . Using the Cauchy data u(x, 0) = cos x,
we find f (x) = cos x. Hence, the unique solution is
u(x, y) = ey cos (x − y)
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 .
Solution.
(a) The characteristic equations in non-parametric form are
dx dy du
= = .
1 y u
Using the first two fractions we find y = C1 ex . Using the first and the last
fractions we find u = C2 ex . Thus, the general solution is given by
f (ye−x , ue−x ) = 0
or
u = ex f (ye−x )
64 CONTENTS
where f is an arbitrary differentiable function.

(b) Letting
Γ : x0 (t) = t, y0 (t) = 3et , u0 (t) = 2
we find
dy0 dx0
a(x0 , y0 , u0 ) − b(x0 , y0 , u0 ) =0
dt dt
and
dx0 du0
2 = c(x0 , y0 , u0 ) 6= a(x0 , y0 , u0 ) =0
dt dt
so the problem has no solutions. Indeed, we want 2 = u(x, 3ex ) = ex f (3ex e−x ) =
ex f (3). This equation is impossible so this Cauchy problem has no solutions.
(c) Letting
Γ : x0 (t) = t, y0 (t) = et , u0 (t) = et
we find
dy0 dx0
a(x0 , y0 , u0 ) − b(x0 , y0 , u0 ) =0
dt dt
dx0 du0
c(x0 , y0 , u0 ) = a(x0 , y0 , u0 ) = et
dt dt
and
dy0 du0
c(x0 , y0 , u0 ) = b(x0 , y0 , u0 ) = e2t
dt dt
so the problem has infinitely many solutions. Indeed, we want ex = ex f (ex e−x ) =⇒
f (1) = 1. In this case, there are infinitely many solutions to this Cauchy prob-
lem, namely, u(x, y) = ex f (ye−x ) where f is an arbitrary function satisfying
f (1) = 1
Problem 10.14
Solve the Cauchy problem
ux + 4uy = x(u + 1), u(x, 5x) = 1.
Solution.
Letting
Γ : x0 (t) = t, y0 (t) = 5t, u0 (t) = 1
we find
dy0 dx0
a(x0 , y0 , u0 ) − b(x0 , y0 , u0 ) = 5 − 4 = 1 6= 0
dt dt
so the problem has a unique solution.

The characteristic equations are dx
1
= dy
4
= du
x(u+1)
. Solving we find 4x−y = c1
x2
and u + 1 = f (4x − y)e . Using the condition u(x, 5x) = 1 we obtain
2
x2
f (x) = 2e− 2 . Thus,
x2 (4x−y)2
u(x, y) = −1 + 2e 2 e− 2
Problem 10.15
π
ux − uy = u, u(x, −x) = sin x, x 6= .
4
Solution.
Letting
Γ : x0 (t) = t, y0 (t) = −t, u0 (t) = sin t
we find
dy0 dx0
a(x0 , y0 , u0 ) − b(x0 , y0 , u0 ) = −1 + 1 = 0
dt dt
and
dx0 du0
c(x0 , y0 , u0 )− a(x0 , y0 , u0 ) = sin t − cos t 6= 0
dt dt
so the problem has no solutions.
dy
The characteristic equations are dx
1
= −1 = duu
. Solving we find x + y = c1
x
and u = f (x + y)e . Using the condition u(x, −x) = sin x we find f (0) =
e−x sin x which is an impossible equation. Hence, the Cauchy problem has
no solutions
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) = e−y .
Solution.
dy
(a) The characteristics satisfy the ODE dx = xy . Solving this equation we find
x2 − y 2 = C. Thus, the characteristics are hyperbolas.
66 CONTENTS
(b)
(c) The general solution to the PDE is u(x, y) = f (x2 − y 2 ) where f is

2
an arbitrary differentiable function. Since u(0, y) = e−y we find f (y) = ey .
2 2
Hence, u(x, y) = ex −y
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.
Solution.
(a) Let
Γ : x0 (t) = t, y0 (t) = 0, u0 (t) = 1.
Then
dy0 dx0
a(x0 , y0 , u0 ) − b(x0 , y0 , u0 ) =0
dt dt
dx0 du0
c(x0 , y0 , u0 ) = a(x0 , y0 , u0 ) =0
dt dt
and
dy0 du0
c(x0 , y0 , u0 ) = b(x0 , y0 , u0 ) =0
dt dt
so the problem has infinitely many solutions.

dy
Solving the ODE dx = y we find the characteristics ye−x = C. Thus, u(x, y) =
f (ye−x ). If u(x, 0) = 1 then we choose f to be any arbitrary differentiable
function satisfying f (0) = 1.
(b) Let
Γ : x0 (t) = t, y0 (t) = 0, u0 (t) = t.
Then
dy0 dx0
a(x0 , y0 , u0 ) − b(x0 , y0 , u0 ) =0
dt dt
and
dx0 du0
6= a(x0 , y0 , u0 )
0 = c(x0 , y0 , u0 ) =1
dt dt
so the problem has no solutions. Indeed, the condition u(x, 0) = x implies
x = f (0) which is impossible
68 CONTENTS
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) Laplace’s equation: ∆u = uxx + uyy = 0.
Solution.
(a) We have A = 1, B = 0 and C = −c2 so that B 2 − 4AC = 4c2 > 0. Thus,
the given equation is of hyperbolic type.
(b) We have A = 0, B = 0 and C = c so that B 2 − 4AC = 0. Thus, the given
equation is of parabolic type.
(c) We have A = 1, B = 0 and C = 1 so that B 2 − 4AC = −4 < 0. Thus,
the given equation is of elliptic type
Problem 11.2
Classify the following linear scalar PDE with constant coefficents as hyper-
bolic, parabolic or elliptic.
(a) uxx + 4uxy + 5uyy + ux + 2uy = 0.
(b) uxx − 4uxy + 4uyy + 3ux + 4u = 0.
(c) uxx + 2uxy − 3uyy + 2ux + 6uy = 0.
Solution.
(a) We have A = 1, B = 4 and C = 5 so that B 2 − 4AC = 16 − 20 = −4 < 0.
Thus, the given equation is of elliptic type.
(b) We have A = 1, B = −4 and C = 4 so that B 2 − 4AC = 16 − 16 = 0.
Thus, the given equation is of parabolic type.
(c) We have A = 1, B = 2 and C = −3 so that B 2 − 4AC = 4 + 12 = 16 > 0.
Thus, the given equation is of hyperbolic type
Problem 11.3
Find the region(s) in the xy−plane where the equation
(1 + x)uxx + 2xyuxy − y 2 uyy = 0
is elliptic, hyperbolic, or parabolic. Sketch these regions.

Solution.
We have A = 1 + x, B = 2xy, and C = −y 2 so that B 2 − 4AC = 4x2 y 2 +
4y 2 (1 + x) = 4y 2 (x2 + x + 1).
• 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 xy−plane with the x−axis 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 x−axis.
• The PDE is of elliptic type if 4y 2 (x2 + x + 1) < 0 which can not happen
Problem 11.4
Show that u(x, t) = cos x sin t is a solution to the problem
utt = uxx
u(x, 0) = 0
ut (x, 0) = cos x
ux (0, t) = 0
for all x, t > 0.
Solution.
We have
ux (x, t) = − sin x sin t,

uxx (x, t) = − cos x sin t,
ut (x, t) = cos x cos t,
utt (x, t) = − cos x sin t.
Thus,
uxx (x, t) = − cos x sin t = utt (x, t),

u(x, 0) = cos x sin 0 = 0,
ut (x, 0) = cos x cos 0 = cos x,
ux (0, t) = − sin 0 sin t = 0
Problem 11.5
Classify each of the following PDE as linear, quasilinear, semi-linear, or non-
linear.
70 CONTENTS
(a) ut + uux = uuxx

(b) xutt + tuyy + u3 u2x = t + 1
(c) utt = c2 uxx
(d) u2tt + ux = 0.
Solution.
(a) Quasi-linear (b) Semi-linear (c) Linear (d) Nonlinear
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
is a solution of
x2 +y 2 +z 2
uxx + uyy + uzz = 0.

Solution.
We have
2x
ux =
x2
+ y2
2y 2 − 2x2
uxx = 2
(x + y 2 )2
2y
uy = 2
x + y2
2x2 − 2y 2
uyy = 2
(x + y 2 )2
Plugging these expressions into the equation we find uxx + uyy = 0. Similar
argument holds for the second part of the problem
Problem 11.7
Consider the eigenvalue problem
uxx = λu, 0 < x < L
ux (0) = k0 u(0)
ux (L) = −kL u(L)
with Robin boundary conditions, where k0 and kL are given positive numbers
and u = u(x). Can this system have a nontrivial solution u 6≡ 0 for λ > 0?
Hint: Multiply the first equation by u and integrate over x ∈ [0, L].
Solution.
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
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
Problem 11.8
Show that u(x, y) = f (x)g(x), where f and g are arbitrary differentiable
functions, is a solution to the PDE
uuxy = ux uy .
Solution.
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 solu-
tion
Problem 11.9
Show that for any n ∈ N, the function un (x, y) = sin nx sinh ny is a solution
to the Laplace equation
∆u = uxx + uyy = 0.
Solution.
We have
(un )xx = −n2 sin nx sinh ny and (un )yy = n2 sin nx sinh ny
Hence, ∆un = 0
72 CONTENTS
Problem 11.10
Solve
uxy = xy.
Solution.
2
Integrate both sides with respect to y to obtain ux (x, y) = xy2 + f (x). Next,
2 2
integrateR both sides w.r.t. x to obtain u(x, y) = x 4y + F (x) + G(y), where
F (x) = f (x)dx can be any function of x, since f (x) itself is an arbitrary
function of x
Problem 11.11
Classify each of the following second-oder PDEs according to whether they
are hyperbolic, parabolic, or elliptic:
(a) 2uxx − 4uxy + 7uyy − u = 0.
(b) uxx − 2 cos xuxy − sin2 xuyy = 0.
(c) yuxx + 2(x − 1)uxy − (y + 2)uyy = 0.
Solution.
(a) We have A = 2, B = −4, C = 7 so B 2 − 4AC = 16 − 56 = −40 < 0. So
this equation is elliptic everywhere in R2 .
(b) We have A = 1, B = −2 cos x, C = − sin2 x so B 2 − 4AC = 4 cos2 x +
4 sin2 x = 4 > 0. So this equation is hyperbolic everywhere in R2 .
(c) We have A = y, B = 2(x − 1), C = −(y + 2) so B 2 − 4AC = 4(x −
1)2 + 4y(y + 2) = 4[(x − 1)2 + (y + 1)2 − 1]. The equation is parabolic if
(x − 1)2 + (y + 1)2 = 1. It is hyperbolic if (x − 1)2 + (y + 1)2 > 1 and elliptic
if (x − 1)2 + (y + 1)2 < 1
Problem 11.12
Let c > 0. By computing ux , uxx , ut , and utt show that
Z x+ct
1 1
u(x, t) = (f (x + ct) + f (x − ct)) + g(s)ds
2 2c x−ct
is a solution to the PDE

utt = c2 uxx
where f is twice differentiable function and g is a differentiable function.
Then compute and simplify u(x, 0) and ut (x, 0).
Solution.
1 1
ut (x, t) = (cf 0 (x + ct) − cf 0 (x − ct)) + [g(x + ct)(c) − g(x − ct)(−c))
2 2c
c 0 0 1
= (f (x + ct) − f (x − ct)) + (g(x + ct) + g(x − ct))
2 2
c2 00 c
utt = (f (x + ct) + f 00 (x − ct)) + (g 0 (x + ct) − g 0 (c − xt))
2 2
1 0 1
ux (x, t) = (f (x + ct) + f 0 (x − ct)) + [g(x + ct) − g(x − ct)]
2 2c
1 00 1
uxx (x, t) = (f (x + ct) + f 00 (x − ct)) + [g 0 (x + ct) − g 0 (x − ct)].
2 2c
By substitutition we see that c2 uxx = utt . Moreover,
1 x
Z
1
u(x, 0) = (f (x) + f (x)) + g(s)ds = f (x)
2 2c x
and
ut (x, 0) = g(x)
Problem 11.13
yuxx + uxy − x2 uyy − ux − u = 0.
Determine the region D in R2 , if such a region exists, that makes this PDE:
(a) hyperbolic, (b) parabolic, (c) elliptic.
Solution.
We have A = y, B = 1, and C = −x2 . Thus, B 2 − 4AC = 1 + 4yx2 . We have
then (a) 1 + 4x2 y > 0, (b) 1 + 4x2 y = 0, (c) 1 + 4x2 y < 0
Problem 11.14
Consider the second-order hyperbolic PDE
uxx + 2uxy − 3uyy = 0.
Use the change of variables v(x, y) = y − 3x and w(x, y) = x + y to solve the

given equation.
74 CONTENTS
Solution.
We have
ux = − 3uv + uw
uxy = − 3uvv − 2uvw + uww
uy =uv + uw
uyy =(uv + uw )v + (uv + uw )w = uvv + 2uvw + uww
Substituting into the PDE we find uvw = 0. Solving this equation we find
u(v, w) = f (v) + g(w). In terms of x and y we have
u(x, y) = f (y − 3x) + g(x + y)
Problem 11.15
uxx + 2uxy − 3uyy = 0.
u(x, 2x) = 1, ux (x, 2x) = x.
Solution.
From the previous exercise we have
u(x, y) = f (y − 3x) + g(x + y).
From the Cauchy data u(x, 2x) = 1 we find
1 = f (−x) + g(3x). (11.3)
Now from the Cauchy data ux (x, 2x) = x we find
x = −3f 0 (−x) + g 0 (3x). (11.4)
Differentiate (11.3) with respect to x we find
−f 0 (−x) + 3g 0 (3x) = 0 (11.5)
Multiply (11.5) by −3 to obtain
3f 0 (−x) − 9g 0 (3x) = 0 (11.6)

Add (11.4) and (11.6) to obtain x = −8g 0 (3x) or g 0 (x) = − 24

x
. Integrating,
x2
we find g(x) = − 48 + C. Now, from (11.3) we have 1 = f (−x) + g(3x) or
9 2
f (−x) = 1 + 48 x − C = f (x). Thus,
10x2 + y 2 − 7xy + 6
u(x, y) = f (y − 3x) + g(x + y) =
6
76 CONTENTS
Problem 12.1
Show that if v(x, t) and w(x, t) satisfy equation 12.2 then αv + βw is also a
solution to 12.2, where α and β are constants.
Solution.
Let z(x, t) = αv(x, t) + βw(x, t). Then we have
c2 zxx =c2 αvxx + c2 βwxx

=αvtt + βvtt
=ztt
Problem 12.2
Show that any linear time independent function u(x, t) = ax + b is a solution
to equation 12.2.
Solution.
Indeed we have c2 uxx (x, t) = 0 = utt (x, t)
Problem 12.3
Find a solution to 12.2 that satisfies the homogeneous conditions u(x, 0) =
u(0, t) = u(L, t) = 0.
Solution.
Clearly the trivial solution u(x, t) = 0 for all x and t is an answer to the
question
Problem 12.4
utt =9uxx
u(x, 0) = cos x
ut (x, 0) =0.
Solution.
According to Example 12.1, the unique solution is given by
1
u(x, t) = (cos (x − 3t) + cos (x + 3t))
2
Problem 12.5
utt =uxx
1
u(x, 0) =
1 + x2
ut (x, 0) =0.
Solution.
According to Example 12.1 with w(x) = 0, the unique solution is given by

1 1 1
u(x, t) = +
2 1 + (x + t)2 1 + (x − t)2
Problem 12.6
utt =4uxx
u(x, 0) =1
ut (x, 0) = cos (2πx).
Solution.
We have v(x) = 1 and w(x) = cos (2πx). The unique solution is given by
1 x+2t
Z
1
u(x, t) = [2 + cos (2πs)ds]
2 2 x−2t
x+2t
1 1
=1 + sin (2πs)
4 2π x−2t
1
=1 + [sin (2πx + 4πt) − sin (2πx − 4πt)]
8π
Problem 12.7
utt =25uxx
u(x, 0) =v(x)
ut (x, 0) =0
where
1 if x < 0
v(x) =
0 if x ≥ 0.
78 CONTENTS
Solution.
The general solution is given by
1
u(x, t) = (v(x + 5t) + v(x − 5t)).
2
Thus, 

 1if x − 5t < 0 and x + 5t < 0
 1
2
if x − 5t < 0 and x + 5t > 0
u(x, t) = 1
if x − 5t > 0 and x + 5t < 0
 2


0 if x − 5t > 0 and x + 5t > 0
Problem 12.8
utt =c2 uxx
2
u(x, 0) =e−x
ut (x, 0) = cos2 x.
Solution.
We have
Z x+ct
1 2 2 1
u(x, t) = [e−(x+ct) + e−(x−ct) ] + cos2 sds
2 2c x−ct
x+ct
1 2 −(x−ct)2 1 s sin 2s
= [e−(x+ct) +e ]+ +
2 2c 2 4 x−ct
1 2 2 t 1
= [e−(x+ct) + e−(x−ct) ] + + [sin (2x + 2ct) − sin (2x − 2ct)]
2 2 8c
1 2 2 t 1
= [e−(x+ct) + e−(x−ct) ] + + cos (2x) sin (2ct)
2 2 4c
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.
Solution.
Just plug the translated/differentiated/dialated solution into the wave equa-
tion and check that it is a solution
Problem 12.10
v(r)
Find v(r) if u(r, t) = r
cos nt is a solution to the PDE
2
urr + ur = utt .
r
Solution.
We have
v 0 (r) v(r)
ur = cos nt − 2 cos nt
r r
v 00 (r) v 0 (r) v(r)
urr = cos nt − 2 2 cos nt + 2 3 cos nt
r r r
v(r)
ut =−n sin nt
r
v(r)
utt = − n2 cos nt
r
Avoiding the trivial solution u = 0, we cancel cos nt and find from urr + 2r ur =
utt the ODE
v 00 (r) v(r)
= −n2
r r
or
v 00 (r) + n2 v(r) = 0.
Solving this equation we find
v(r) = A cos (nr) + B sin (nr)
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.
Solution.
The general solution is given by
Z x+ct
1 1
u(x, t) = [ex−ct + ex+ct + sin sds]
2 c x−ct
Thus,
1 1
u(x, t) = [ex−ct + ex+ct + (cos (x − ct) − cos (x + ct))]
2 c
80 CONTENTS
Problem 12.12
The total energy of the string (the sum of the kinetic and potential energies)
is defined as
1 L 2
Z
E(t) = (ut + c2 u2x )dx.
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.
Solution.
(a) We have
Z L Z L
dE
(t) = ut utt dx + c2 ux uxt dx
dt 0 0
Z L Z L
2 2 2
= ut utt dx + c ut (L, t)ux (L, t) − c ut (0, t)ux (0, t) − c ut uxx dx
0 0
Z L
2 2
=c ut (L, t)ux (L, t) − c ut (0, t)ux (0, t) + ut (utt − c2 uxx )dx
0
2
=c (ut (L, t)ux (L, t) − ut (0, t)ux (0, t))
since utt − c2 uxx = 0.
(b) Since the ends are fixed, we have ut (0, t) = ut (L, t) = 0. From (a) we
have
dE
(t) = c2 (ut (L, t)ux (L, t) − ut (0, t)ux (0, t)) = 0.
dt
(c) Assuming free ends boundary conditions, that is ux (0, t) = ux (L, t) = 0,
we find dEdt
(t) = 0
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.
Solution.
Using the previous exercise, we find
Z L
dE
(t) = −d (ut )2 dx.
dt 0
The right-hand side is nonpositive, so the energy either decreases or is con-

stant. The latter case can occur only if ut (x, t) is identically zero, which
means that the string is at rest
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).
Solution.
(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 .
(b) We have
L L L
c2 0 c2
Z Z Z
1 2
(ut ) dx = [R (x − ct)]2 dx = (ux )2 dx
2 0 0 2 0 2
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.

82 CONTENTS
Solution.
We have
Z x+2t
1 1
u(x, t) = [(x − 2t)2 + (x + 2t)2 ] + sin 2sds
2 4 x−2t
x+2t
1
=x2 + 4t2 − cos 2s
8 x−2t
2 1
2 1
=x + 4t − cos (2x + 4t) + cos (2x − 4t)
8 8
2 2 1
=x + 4t + sin 2x sin 4t
4
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.
Solution.
Let z(x, t) = αu(x, t) + βv(x, t). Then we have
kzxx =kαuxx + kβvxx

=αut + βvt
=zt
Problem 13.2
Show that any linear time independent function u(x, t) = ax + b is a solution
to equation (13.1).
Solution.
Indeed we have kuxx (x, t) = 0 = ut (x, t)
Problem 13.3
Find a linear time independent solution u to (13.1) that satisfies u(0, t) = T0
and u(L, T ) = TL .
Solution.
Let u(x, t) = ax + b. From the assumptions of the problem we must have
b = T0 and a = TLL−T0 . Thus, u(x, t) = T0 + TLL−T0 x
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.
Solution.
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 + TLL−T0 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
84 CONTENTS
Problem 13.5
Find a solution to (13.1) that satisfies the conditions u(x, 0) = u(0, t) =
u(L, t) = 0.
Solution.
Clearly the trivial solution u(x, t) = 0 for all x and t is an answer to the
question
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 λ.
Solution.
Substituting u(x, t) = X(x)T (t) into (13.1) and then separating the variables
we obtain
X 00 T0
k = .
X T
Since X only depends on x and T only depends on t, there must be a constant
λ such that
00 T0
k XX = λ and T
= λ.
This gives the two ordinary differential equations
X 00 − λk X = 0 and T 0 − λT = 0
Problem 13.7
Consider again the solution u(x, t) = X(x)T (t). Clearly, T (t) = T (0)eλt .
Suppose that λ > 0. √ √
(a) Show that X(x) = Aex α + Be−x α , 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 α −
e−L α ) = 0.
(c) Show that A = 0 leads to a contradiction.
√ √
(d) Using (b) and (c) show that eL α = e−L α . Show that this equality leads
to a contradiction. We conclude that λ < 0.
Solution.
(a) Letting α = λk > 0 we obtain √
the ODE√
X 00 − αX = 0 whose general
solution is given by X(x) = Aex α + Be−x α for some constants A and B.
(b) The condition u(0, t) = 0 implies that X(0) = 0 which √
in turn √implies
A + B = 0. Likewise,
√ √
the condition u(L, t) = 0 implies AeL α
+ Be−L α = 0.
Hence, A(eL α − e−L α ) = 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.
(d) Using (b) and (c) √ we obtain e L α
= e−L α or e2L α = 1. This equation
is impossible
√ since 2L α√> 0. Hence, we must have λ < 0 so that X(x) =
A cos (x −α) + B sin (x −α)
Problem 13.8
Consider the results of the previous exercise. q
−λ
(a) Show that X(x) = c1 cos βx + c2 sin βx where β = k
.
2 2
(b) Show that λ = λn = − knL2π , where n is an integer.
Solution. q
(a)Now, write β = − λk . Then we obtain the equation X 00 + β 2 X = 0 whose
general solution is given by
X(x) = c1 cos βx + c2 sin βx.
(b) Using X(0) = 0 we obtain c1 = 0. Since c2 6= 0 we must have sin βL = 0

2 2
which implies βL = nπ where n is an integer. Thus, λ = − knL2π , where n is
an integer
Problem 13.9
ki2 π 2
Show that u(x, t) = ni=1 ui (x, t), where ui (x, t) = ci e− L2 t sin iπ
P
L
x satis-
fies (13.1) and the homogeneous boundary conditions.
Solution.
ki2 π 2
For each integer i ≥ 0 we have ui (x, t) = ci e− L2 t sin iπ

L
x is a solution
to (13.1). By superposition, u(x, t) is also a solution to (13.1). Moreover,
u(0, t) = u(L, t) = 0 since ui (0, t) = ui (L, t) = 0 for i = 1, · · · , n
86 CONTENTS
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.
Solution.
(i) u(0, t) = 0 and u(a, t) = 100 for t > 0.
(ii) ux (0, t) = ux (a, t) = 0 for t > 0
Problem 13.11
Let ut = uxx for 0 < x < π and t > 0 with boundary conditionsRπ u(0, t) =
0 = u(π, t) and initial condition u(x, 0) = sin x. Let E(t) = 0 (u2t + u2x )dx.
Show that E 0 (t) < 0.
Solution.
Solving this problem we find u(x, t) = e−t sin x. We have
Z π Z π
−2t −2t
E(t) = 2 2
[e sin x + e cos x]dx = e−2t dx = πe−2t .
0 0
Thus, E 0 (t) = −2πe−2t < 0 for all 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).
Solution.
We have Z L
d
u(x, t)dx = ux |L0 + 4L = 1 + 4L.
dt 0
Thus, Z L
E(t) = f (x)dx + (1 + 4L)t
0
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 don’t depend on time. Find v(x).
Solution.
Solving the equation v 00 (x) = 0 we find v(x) = ax + b. Using the conditions
v(0) = 2 and v(1) = 3 we find v(x) = x + 2
Problem 13.14
Consider the equation for the one-dimensional rod of length L with given
heat energy source:
ut = uxx + q(x).
Assume that the initial temperature distribution is given by u(x, 0) = f (x).
Find the equilibrium (steady state) temperature distribution in the following
cases.
(a) q(x) = 0, u(0, t) = 0, u(L, t) = T.
(b) q(x) = 0, ux (0, t) = 0, u(L, t) = T.
(c) q(x) = 0, u(0, t) = T, ux (L, t) = α.
Solution.
Recall that a steady-state solution is a solution that does not depend on time
(i.e. ut = 0.).
(a) We have v 00 (x) = 0 =⇒ v(x) = c1 x + c2 . But v(0) = 0 and v(L) = T so
that c1 = TL and c2 = 0. Thus, the steady-state solution is v(x) = TL x.
(b) We have v(x) = c1 x + c2 with v 0 (0) = 0 and v(L) = T. Thus, c1 = 0 and
c2 = T so that v(x) = T.
(c) We have v(x) = c1 x + c2 with v(0) = T and v 0 (L) = α. Thus, c1 = α and
c2 = T so that v(x) = αx + T
Problem 13.15
Consider the equation for the one-dimensional rod of length L with insulated
ends:
ut = kuxx , ux (0, t) = ux (L, t) = 0.
88 CONTENTS
(a) Give the expression for the total thermal energy of the rod.
(b) Show using the equation and the boundary conditions that the total
thermal energy is constant.
Solution. R
L
(a) E(t) = 0 u(x, t)dx.
(b) We integrate the equation in x from 0 to L :
Z L Z L
ut (x, t)dx = kuxx dx = kux (x, t)|L0 = 0,
0 0
since ux (0, t) = ux (L, t) = 0. The left-hand side can also be written as
d L
Z
u(x, t)dx = E 0 (t).
dt 0
Thus, we have shown that E 0 (t) = 0 so that E(t) is constant
Problem 13.16
Suppose
ut = uxx + x, u(x, 0) = f (x), ux (0, t) = β, ux (L, t) = 7.
(a) Calculate the total thermal energy of the one-dimensional rod (as a func-
tion of time).
(b) From part (a) find the value of β for which a steady-state solution exist.
(c) For the above value of β find the steady state solution.
Solution.
(a) The total thermal energy is
Z L
E(t) = u(x, t)dx.
0
We have
L L
L2
Z Z
dE
= ut (x, t)dx = ux |L0 + xdx = (7 − β) + .
dt 0 0 2
Hence,
L
L2
Z
E(t) = f (x)dx + (7 − β) + t.
0 2
(b) The steady solution (equilibrium) is possible if E 0 (t) = 0 :
L2
(7 − β) + =0
2
2
Solving this equation for β we find β = 7 + L2 .
(c) By integrating the equation uxx + x = 0 we find the steady solution
x3
u(x) = − + C1 x + C2
6
From the condition ux (0) = β we find C1 = β. The steady solution should
also have the same value of the total energy as the initial condition. This
means Z L 3 Z L
x
− + βx + C2 dx = f (x)dx = E(0).
0 6 0
Performing the integration and then solving for C2 we find
1 L L3
Z
L
C2 = f (x)dx + −β .
L 0 24 2
Therefore, the steady-state solution is
1 L L3 x3
Z
L
u(x) = f (x)dx + − β + βx −
L 0 24 2 6
90 CONTENTS
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.
Solution.
(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 = 12 . Then there exists a positive integer N
such that for all n ≥ N we have
1
|fn (x) − f (x)| <
2
for all x ∈ [0, 1]. In particular, we have
1
|fN (x) − f (x)| <
2
1
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
Problem 14.2
Consider the sequence of functions
nx + x2
fn (x) =
n2
defined for all x in R. Show that this sequence converges pointwise to a
function f to be determined.
Solution.
For every real number x, we have
nx + x2 x x2
lim fn (x) = lim = lim + lim =0
n→∞ n→∞ n2 n→∞ n n→∞ n2
Thus, {fn }∞
n=1 converges pointwise to the zero function on R
Problem 14.3
Consider the sequence of functions
sin (nx + 3)
fn (x) = √
n+1
defined for all x in R. Show that this sequence converges pointwise to a
function f to be determined.
Solution.
For every real number x, we have
1 1
−√ ≤ fn (x) ≤ √ .
n+1 n+1
Moreover,
1
lim √ = 0.
n→∞ n+1
Applying the squeeze rule for sequences, we obtain
lim fn (x) = 0
n→∞
for all x in R. Thus, {fn }∞

n=1 converges pointwise to the zero function on R
Problem 14.4
Consider the sequence of functions defined by fn (x) = n2 xn for all 0 ≤ x ≤ 1.
Show that this sequence does not converge pointwise to any function.
Solution.
First of all, observe that fn (0) = 0 for every n in N. So the sequence
{fn (0)}∞
n=1 is constant and converges to zero. Now suppose 0 < x < 1
then n x = n2 en ln x . But ln x < 0 when 0 < x < 1, it follows that
2 n
92 CONTENTS
limn→∞ fn (x) = 0 for 0 < x < 1
Finally, fn (1) = n2 for all n. So,
lim fn (1) = ∞.
n→∞
Therefore, {fn }∞
n=1 is not pointwise convergent on [0, 1]
Problem 14.5
Consider the sequence of functions defined by fn (x) = (cos x)n for all − π2 ≤
x ≤ π2 . Show that this sequence converges pointwise to a noncontinuous
function to be determined.
Solution.
For − π2 ≤ x < 0 and 0 < x ≤ π
2
we have
lim (cos x)n = 0.

n→∞
For x = 0 we have fn (0) = 1 for all n in N. Therefore, {fn }∞

n=1 converges
pointwise to
0 if − π2 ≤ x < 0 and 0 < x ≤ π2

f (x) =
1 if x = 0
Problem 14.6
n
Consider the sequence of functions fn (x) = x − xn defined on [0, 1).
(a) Does {fn }∞n=1 converge to some limit function? If so, find the limit func-
tion and show whether the convergence is pointwise or uniform.
(b) Does {fn0 }∞
n=1 converge to some limit function? If so, find the limit func-
tion and show whether the convergence is pointwise or uniform.
Solution.
(a) Let > 0 be given. Let N be a positive integer such that N > 1 . Then
for n ≥ N
n n

x − x − x = |x| < 1 ≤ 1 < .

n n n N
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 = 21 and suppose that the convergence is uniform. Then there is a
positive integer N such that for n ≥ N we have
1
|1 − xn−1 − 1| = |x|n−1 < .
2
In particular, if we let n = N + 1 we must have xN < 21 for all x ∈ [0, 1).
1
But x = 21 N ∈ [0, 1) and xN = 12 which contradicts xN < 12 . Hence, the
convergence is not uniform
Problem 14.7
xn
Let fn (x) = 1+x n for x ∈ [0, 2].
(a) Find the pointwise limit f (x) = limn→∞ fn (x) on [0, 2].
(b) Does fn → f uniformly on [0, 2]?
Solution.
(a) The pointwise limit is

 0 if 0 ≤ x < 1
1
f (x) = if x = 1
 2
1 if 1 < x ≤ 2
(b) The convergence cannot be uniform because if it were f would have to
be continuous
Problem 14.8
n+cos x
For each n ∈ N define fn : R → R by fn (x) = 2n+sin2 x
.
(a) Show that fn → 21 uniformly.
R7
(b) Find limn→∞ 2 fn (x)dx.
Solution.
(a) Let > 0 be given. Note that
2 cos x − sin2 x

1 ≤ 3 .
|fn (x) − | = 2
2 2(2n + sin x) 4n
3
Since limn→∞ 4n = 0 we can find a positive integer N such that if n ≥ N
3
then 4n < . Thus, for n ≥ N and all x ∈ R we have
1 3
|fn (x) − | ≤ < .
2 4n
94 CONTENTS
This shows that fn → 21 uniformly on R and also on [2, 7].

(b) We have
Z 7 Z 7 Z 7
1 5
lim fn (x)dx = lim fn (x)dx = dx =
n→∞ 2 2 n→∞ 2 2 2
Problem 14.9
Show that the sequence defined by fn (x) = (cos x)n does not converge uni-
formly on [− π2 , π2 ].
Solution.
We have proved earlier that this sequence converges pointwise to the discon-
tinuous 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
Problem 14.10
Let {fn }∞
n=1 be a sequence of functions such that
2n
sup{|fn (x)| : 2 ≤ x ≤ 5} ≤ .
1 + 4n
(a) Show that this sequence converges uniformly
R5 to a function f to be found.
(b) What is the value of the limit limn→∞ 2 fn (x)dx?
Solution.
(a) Using the squeeze rule we find
lim sup{|fn (x)| : 2 ≤ x ≤ 5} = 0.

n→∞
Thus, {fn }∞
n=1 converges uniformly to the zero function.
(b) We have Z 5 Z 5
lim fn (x)dx = 0dx = 0
n→∞ 2 2
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.
Solution.
(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)
Problem 15.2
Show
R Lthat for m6= n we have
(a) −L sin L x sin nπ
mπ

L
x dx = 0 and
RL mπ
nπ

(b) −L cos L x sin L x dx = 0.
Solution.
(a) For n 6= m we have
L
1 L
Z
(m − n)π
Z mπ nπ (m + n)π
sin x sin x dx = − cos x − cos x dx
−L L L 2 −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
96 CONTENTS
(b) For n 6= m we have

L
1 L
Z
(m − n)π
Z mπ nπ (m + n)π
cos x sin x dx = sin x − sin x dx
−L L L 2 −L L L

1 L (m + n)π
= − cos x
2 (m + n)π L
L
L (m − n)π
+ cos x
(m − n)π L −L
=0
where we used the trigonometric identiy
1
cos a sin b = [sin (a + b) − sin (a − b)]
2
Problem 15.3
Compute
RL the following integrals:
2 nπ

(a) −L cos L x dx.
RL
(b) −L sin2 nπ

L
x dx.
RL
(c) −L cos nπ nπ

L
x sin L
x dx.
Solution.
1+cos 2a
(a) Using the trigonometric identity cos2 a = 2
we can write
L
1 L
Z nπ Z
2 2nπ
cos x dx = 1 + cos x dx
−L L 2 −L L
L
1 L 2nπ
= x+ sin x = L.
2 2nπ L −L
(b) Using the trigonometric identity sin2 a = 1−cos 2a

2
we can write
L
1 L
Z nπ Z
2 2nπ
sin x dx = 1 − cos x dx
−L L 2 −L L
L
1 L 2nπ
= x− sin x =L
2 2nπ L −L
sin 2a
(c) Using the trigonometric identity cos a sin a = 2
we can write
L
1 L
Z nπ nπ Z
2nπ
cos x sin x dx = sin x dx
−L L L 2 −L L
L
L 2nπ
=− cos x =0
4nπ L −L
Problem 15.4
Find the Fourier coefficients of

 −π, −π ≤ x < 0
f (x) = π, 0<x<π
0, x = 0, π

on the interval [−π, π].
Solution.
We have
1 π
Z
a0 = f (x)dx = 0
π −π
1 π
Z
an = f (x) cos nxdx
π −π
Z 0 Z π
=− cos nxdx + cos nxdx = 0
−π 0
1 π
Z
bn = f (x) sin nxdx
π −π
Z 0 Z π
=− sin nxdx + sin nxdx
−π 0
2
= [1 − (−1)n ]
n
Problem 15.5
1
Find the Fourier series of f (x) = x2 − 2
on the interval [−1, 1].
98 CONTENTS
Solution.
We have
Z 1
1 1
a0 = (x2 − )dx = −
−1 2 3
Z 1
1
an = (x2 − ) cos nπxdx
−1 2
Z 1
1 1
Z
2
= x cos nπxdx − cos nπxdx
−1 2 −1
1 Z 1 1
2 sin (nπx) sin (nπx) 1 sin (nπx)
=x − 2x dx −
nπ −1 −1 nπ 2 nπ −1
1 Z 1
cos (nπx) cos (nπx)
=2x − 2 dx
(nπ)2 −1 −1 (nπ)2
4
= 2
(−1)n .
(nπ)
Z 1
1
bn = (x2 − ) sin nxdx = 0.
−1 2
Note that bn = 0 because the integrand is odd. Hence,
∞
1 X 4
f (x) = − + (−1)n cos (nπx), − 1 ≤ x ≤ 1
6 n=1 (nπ)2
Problem 15.6
Find the Fourier series of the function

 −1, −2π < x < −π
f (x) = 0, −π < x < π
1, π < x < 2π.

Solution.
From the graph of f (x) we see that f is an odd function on (−2π, 2π). Thus,
nx

f (x) cos 2
is odd so that an = 0 for all n ∈ N. Now,
Z 2π
1
a0 = f (x)dx = 0
2π −2π
Z 2π
1 nx
bn = f (x) sin dx
2π −2π 2
1 2π
Z nx
= f (x) sin dx
π 0 2
1 2π
Z nx
= sin dx
π π 2
2 nx i2π
=− cos
nπ 2 π
2 h nπ i
= cos − (−1)n
nπ 2
Hence,
∞
X 2 h nπ i nx
f (x) = cos − (−1)n sin , − 2π < x < 2π
n=1
nπ 2 2
Problem 15.7
Find the Fourier series of the function

1 + x, −2 ≤ x ≤ 0
f (x) =
1 − x, 0 < x ≤ 2.
Solution.
From the graph of f (x) we see that f is an even function on [−2, 2]. Thus,
100 CONTENTS
nπ

f (x) sin 2
x is odd so that bn = 0 for all n ∈ N. Now,
1 2
Z Z 0 Z 2
1
a0 = f (x)dx = (1 + x)dx + (1 − x)dx = 0
2 −2 2 −2 0
1 2
Z nπ
an = f (x) cos x dx
2 −2 2
Z 2 nπ
= (1 − x) cos x dx
0 2
Z 2 nπ Z 2 nπ
= cos x dx − x cos x dx
0 2 0 2
4 nπ i2
=− cos x
(nπ)2 2 0
4
= [1 − (−1)n ]
(nπ)2
Hence,
∞
X 4 n
nπ
f (x) = [1 − (−1) ] cos x , −2≤x≤2
n=1
(nπ)2 2
Problem 15.8
1
Show that f (x) = x
is not piecewise continuous on [−1, 1].
Solution.
Since the sided limits at the point of discontinuity x = 0 do not exist, the
function is not piecewise continuous in [−1, 1]
Problem 15.9
Assume that f (x) is continuous and has period 2L. Prove that
Z L Z L+a
f (x)dx = f (x)dx
−L −L+a
is independent of a ∈ R. In particular, it does not matter over which interval

the Fourier coefficients are computed as long as the interval length is 2L.
[Remark: This result is also true for piecewise continuous functions].
Solution.
Define the function
Z L+a Z L+a Z −L+a
g(a) = f (x)dx = f (x)dx − f (x)dx.
−L+a 0 0
Using the fundamental theorem of calculus, we have
Z L+a Z −L+a
dg d
= f (x)dx − f (x)dx
da da 0 0
=f (L + a) − f (−L + a) = f (−L + a + 2L) − f (−L + a)
=f (−L + a) − f (−L + a) = 0
Hence, g is a constant function, and in particular we can write g(a) = g(0)

for all a ∈ R which gives the desired result
Problem 15.10
Consider the function f (x) defined by

1 0≤x<1
f (x) =
2 1≤x<
and extended periodically with period 3 to R so that f (x + 3) = f (x) for all

x.
(i) Find the Fourier series of f (x).
(ii) Discuss its limit: In particular, does the Fourier series converge pointwise
or uniformly to its limit, and what is this limit?
(iii) Plot the graphs of f (x) and the limit of the Fourier series.
Solution.
(i) The Fourier series is computed for functions of period 2L. Since this
function has period 3, L = 3/2. By the previous problem, we can compute
the coefficients over any interval of length 3, so we might as well use [0, 3].
102 CONTENTS
Using the formulas for the coefficients, we obtain:

Z L Z 1 Z 3
1 2 10
a0 = f (x)dx = dx + 2dx =
L −L 3 0 1 3
Z L Z 1 Z 3
1 nπx 2 2nπx 2nπx
an = f (x) cos dx = cos dx + 2 cos dx
L −L L 3 0 3 1 3

2 3 2nπ 2nπ
= sin − 0 + 2 sin 2nπ − sin
3 2nπ 3 3

1 2nπ
=− sin
nπ 3
Z L Z 1 Z 3
1 nπx 2 2nπx 2nπx
bn = f (x) sin dx = sin dx + 2 sin dx
L −L L 3 0 3 1 3

2 3 2nπ 2nπ
=− cos − 1 + 2 cos 2nπ − cos
3 2nπ 3 3

1 2nπ
=− − cos +1
nπ 3
Thus, the Fourier series is

∞
10 X 1 2nπ 2nπx 1 2nπ 2nπx
f (x) = + − sin cos − − cos + 1 sin .
3 nπ 3 3 nπ 3 3
n=1
(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)
Problem 15.11
For the following functions f (x) on the interval −L < x < L, determine the
coefficients an , n = 0, 1, 2, · · · and bn , n ∈ N of the Fourier series expansion.
(a) f (x) = 1.
(b) f (x) = 2 + sin πx

L .
1 x≤0
(c) f (x) =
0 x > 0.
(d) f (x) = x.
Solution.
The detailed computations are left for the reader to fill.
(a) a0 = 2, an = bn = 0 for n ∈ N.
(b) a0 = 4, an = 0, b1 = 1, and bn = 0 for n > 1.
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
Problem 15.12
Let f (t) be the function with period 2π defined as
2 if 0 ≤ x ≤ π2

f (t) =
0 if π2 < x ≤ 2π
f (t) has a Fourier series and that series is equal to
∞
a0 X
+ (an cos nt + bn sin nt), 0 < x < 2π.
2
n=1
a3
Find b3 .
Solution.
We have Z π
1 2 2
a3 = 2 cos 3xdx = −
π 0 3π
104 CONTENTS
and
Z π
1 2 2
b3 = 2 sin 3xdx = .
π 0 3π
a3
Thus, b3 = −1
Problem 15.13
Let f (x) = x3 on [−π, π], extended periodically to all of R. Find the Fourier
coefficients an , n = 1, 2, 3, · · · .
Solution.
Since the extension is an odd function, we must have an = 0 for all n ∈ N
Problem 15.14
Let f (x) be the square wave function

−π −π ≤ x < 0
f (x) =
π 0≤x≤π
extended periodically to all of R. To what value does the Fourier series of f (x)
converge when x = 0?
Solution.
f (x) is piecewise smooth function with discontinuity at x = 0. Thus, the Fourier
series of f (x) at x = 0 converges to
f (0− ) + f (0+ ) −π + π
= =0
2 2
Problem 15.15
(a) Find the Fourier series of

1 −π ≤ x < 0
f (x) =
2 0≤x≤π
extended periodically to all of R. Simplify your coefficients as much as possible.

(−1)n+1
(b) Use (a) to evaluate the series ∞
P
n=1 (2n−1) . Hint: Evaluate the Fourier series
at x = π2 .
Solution.
(a) We have
1 π
Z Z 0 Z π
1
a0 = f (x)dx = dx + dx = 3
π −π π −π 0
Z π
1
an = f (x) cos nxdx = 0
π −π
1 π 1 1 (−1)n
Z
bn = f (x) sin nxdx = −
π −π π n n
Thus,
∞
3 2 X sin (2n − 1)x
f (x) = + .
2 π 2n − 1
n=1
(b) By the convergence theorem we have

∞
π 3 2 X sin (2n − 1) π2
f = + .
2 2 π 2n − 1
n=1
This implies
∞
3 2 X (−1)n+1
2= +
2 π 2n − 1
n=1
and this reduces to

∞
X (−1)n+1 π
=
2n − 1 4
n=1

Solution PDF

Uploaded by

Copyright:

Available Formats

Solution PDF

Uploaded by

Document Information

Original Title

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

Solution PDF

Uploaded by

Copyright:

Available Formats

Lecture Notes in Mathematics

A First Course in Quasi-Linear Partial

March 12, 2018

Substituting these expressions into the given equation we find

Substituting these expressions into the given equation we find

(b) Linear, non-homogeneous, order 3. The inhomogeneity is − sin y.

(b) The wave equation for w(x, t) :

(c) The thin film equation for h(x, t) :

(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.

uxx + 4uxy + 4uyy = 0.

Use the change of variables v(x, y) = y − 2x and w(x, y) = x to show that

Substituting these into the given PDE we find uww = 0

Substituting we find uvw = 0

Substituting into the PDE we find uvw = 0

Substituting these into the given equation to obtain

abus + aut − abus = 0

uxxxx + uyyyy + 2uxxyy = 0.

(a4 + 2a2 b2 + b4 )eax+by = 0.

Since eax+by 6= 0, we must have a4 + 2a2 b2 + b4 = 0 or (a2 + b2 ) = 0. This is

xux + (x + 1)yuy = 0, x, y > 1

has a solution satisfying u(x, t) = f (x + p(u)t), where f is an arbitrary

where f is an arbitrary differentiable function in one variable

uxx + 2uxy + uyy = 0.

Hint: See Problem 1.2.

arbitrary differentiable function.

uxx + 4uxy + 4uyy = 0.

(a) Use the change of variables v(x, y) = y − 2x and w(x, y) = x to show

Substituting these into the given PDE we find uww = 0.

by using the change of variables v = x + ct and w = x − ct.

Substituting we find uvw = 0 and solving Rthis equation we find uv = f (v)

u(x, t) = F (x + ct) + G(x − ct)

where F and G are arbitrary differentiable functions

Since y(0) = 0, we find C = − 12 . Thus, the unique solution to the IVP is

= e−3t e9 (3t − 1) + et + Ce−3t

= − 3e− sin t esin t + Ce− sin t

Since y(0) = 3, we find C = 2 and therefore y(t) = 1 + 2e−(t+sin t) . Finally,

lim y(t) = lim (1 + 2e− sin t e−t ) = 1

But y(1) = 7 so that C = 7. Hence,

has the property that y → 0 as t → ∞. Hint: Consider the cases a = λ and

and the corresponding general solution is

y(t) = bte−at + Ce−at

y(t) = (−tet + et )−1

Since y(1) = 3, C = 3 − Si(1). Hence, y(t) = tSi(t) + (3 − Si(1))t

Separating the variables and solving the equation we find

Since y(0) = 1, we find C = 5 and hence

Since y(1) = 0, we find C = e. Thus, y(t) = e1−t − 1

Since y( π4 ) = −1, we find C = π2 . Hence, y(t) = tan (t + π2 ) = − cot t

Since y(0) = 12 , we find C = − 31 . Thus,

f (y)e−3x dy = F (y)e−3x + G(x) where F (y) =

Hence, by substitution we find 3a − 7b = 0

to a first order ODE by introducing the change of variables s = bx − ay and

by using the change of variables v = ax − bt and w = x.

by using the change of variables s = bx − ay and t = x.

Substituting these into the given equation to obtain

and since a 6= 0 we obtain