Course Handout MAT 201: Partial Differential Equations and Complex Analysis Module 1: Partial Differential Equations
Course Handout MAT 201: Partial Differential Equations and Complex Analysis Module 1: Partial Differential Equations
Course Handout MAT 201: Partial Differential Equations and Complex Analysis Module 1: Partial Differential Equations
Pre requisites
(i) Calculating partial derivatives
Definition 1.2. The order of the highest partial derivative in the equation is called the
order of the PDE.
∂z ∂z ∂2z ∂z 3
For example, z = x ∂x + y ∂y has order 1, and ∂x∂y
= ∂x
has order 2.
Definition 1.3. The degree of the highest order partial derivative is called the degree
of the PDE.
∂z 3
Both equations in the above example have degree 1. But z = ∂x
has degree 3.
(i) If the number of arbitrary constants equals the number of independent variables,
then the PDE is of first order.
(ii) If the number of arbitrary constants is greater than the number of independent
variables, then the PDE is of second or higher order.
(iii) If the PDE is formed by eliminating arbitrary functions, then the order of the
PDE is in general equal to the number of arbitrary functions eliminated.
2.1 Elimination of arbitrary constants
For such problems, we find the partial derivatives of z to eliminate the arbitrary con-
stants.
x2 y2
Example 1. Derive a PDE from the relation 2z = a2
+ b2
.
Solution: The given relation has 2 arbitrary constants a and b. There are 2 independent
variables x, y. So by our rule of thumb, the resulting PDE should be of first order.
Find the first order partial derivatives of z.
∂z 2x ∂z 2y
2 = 2,2 = 2
∂x a ∂y b
Simplifying the above expressions, we find
∂z x ∂z y
p= = 2,q = = 2
∂x a ∂y b
.
Now we can write both arbitrary constants in terms of p, q, x, andy as
1 p 1 q
2
= , 2 =
a x b y
Finally, substitute for the constants in the original relation to find
2z = px + qy
Example 2. Find the differential equation of all planes which are at a constant distance
a from the origin.
Solution: A plane at distance a from the origin will have an equation of the form
lx + my + nz = a, where l, m, n are the direction cosines of the normal from the origin
to the plane.
This relation has 3 arbitrary constants l, m, n. But these constants are the direction
√
cosines, and so satisfy the property : l2 + m2 + n2 = 1, or n = 1 − l2 − m2 . So the
equation of the plane can be written as
√
lx + my + 1 − l2 − m2 z = a (1)
Observe that the relation now has just 2 arbitrary constants, and 2 independent
variables. So the PDE satisfied by this relation should be of first order.
Take the partial derivatives of (1) with respect to x and y. We get the following
equations:
√
l+ 1 − l2 − m2 p = 0 (2)
√
m + 1 − l2 − m2 q = 0 (3)
Next we eliminate the arbitrary constants l, m from equations (1), (2), (3).
Squaring and adding (2) and (3),
l2 + m2 = (1 − l2 − m2 )(p2 + q 2 )
p2 + q 2
l2 + m2 =
1 + p2 + q 2
Hence,
1
1 − l2 − m2 =
1 + p2 + q 2
We can substitute this expression in (2), (3) to find l = − √ p 2 2 , m = − √ q 2 2 .
1+p +q 1+p +q
√
2 2
Finally, substitute for l, m, 1 − l − m in (1)
px qy z
−p −p +p =a
1 + p2 + q 2 1 + p2 + q 2 1 + p2 + q 2
Rearranging this equation gives us the required PDE
p
z = px + qy + a 1 + p2 + q 2
Exercise
Form PDEs from the following equations by eliminating arbitrary constants.
(i) z = ax + by + ab
(ii) z = ax + a2 y 2 + b
Solution: In this relation, the independent variables are x, t and there are 2 arbitrary
functions f, g.
Again, by the rule of thumb, we have to differentiate as many times as the number
of arbitrary functions - in this case 2.
∂2z ∂2z
Evaluate ,
∂x2 ∂y 2
.
∂z ∂z
= f 0 (x + at) + g 0 (x − at), = af 0 (x + at) − ag 0 (x − at)
∂x ∂t
(By applying the chain rule)
∂ 2z 00 00 ∂ 2z
2
= f (x + at) + g (x − at), 2
= a2 f 00 (x + at) + a2 g 00 (x − at)
∂x ∂t
Now by comparing the second order derivatives we can write
∂ 2z 2
2∂ z
= a
∂t2 ∂x2
Note that this is just the wave equation from the examples.
Solution: The given relation has only 1 arbitrary function, so the required PDE should
be of first order.
The partial derivatives of z are:
∂z
p= = (x + y)φ0 (x2 − y 2 ).2x + φ(x2 − y 2 )
∂x
∂z
q= = (x + y)φ0 (x2 − y 2 ).(−2y) + φ(x2 − y 2 )
∂y
(Using product rule and the chain rule)
z
From the original relation, we have φ(x2 − y 2 ) = x+y
. Substituting for φ(x2 − y 2 )
in the partial derivatives and rearranging gives us the following equations:
z
p− = 2x(x + y)φ0 (x2 − y 2 )
x+y
z
q− = −2y(x + y)φ0 (x2 − y 2 )
x+y
Remember that to get the required PDE, we need to eliminate the arbitrary function
φ. It is easy to see that dividing the above 2 equations will eliminate the φ term.
The division gives
p − z/(x + y) x
=−
q − z/(x + y) y
Taking the LCM and simplifying gives
py + qx = z
Example 5. Form the PDE by eliminating arbitrary functions from the given relation:
f (x2 + y 2 , z − xy) = 0
Solution: For questions like these, we can’t evaluate the partial derivatives of z directly.
To solve, let u = x2 + y 2 , v = z − xy, so that the relation is f (u, v) = 0.
Now apply the chain rule to find the partial derivative of the above function with
respect to x.
∂f ∂u ∂x ∂u ∂z ∂f ∂v ∂x ∂v ∂z
+ + + =0
∂u ∂x ∂x ∂z ∂x ∂v ∂x ∂x ∂z ∂x
∂z
If we now substitute for the appropriate partial derivatives of u, v and let p = ∂x
,
then
∂f ∂f
(2x) + (−y + p) = 0 (4)
∂u ∂v
Similarly, we can evaluate the partial derivative with respect to y and get the
equation
∂f ∂f
(2y) + (−x + q) = 0 (5)
∂u ∂v
Now that we have 2 equations involving ∂f , ∂f , we can eliminate these partial
∂u ∂v
derivatives from the equations:
Multiply the first equation by y and the second equation by x. If we then subtract
these 2 new equations, we find
∂f
(−y 2 + py + x2 − qx) = 0
∂v
∂f
Then either ∂v
= 0, or (−y 2 + py + x2 − qx) = 0. But the partial derivative cannot
identically be equal to zero (as f is arbitrary, we cannot say anything specific about
its partial derivatives).
Therefore,
(−y 2 + py + x2 − qx) = 0
Exercise
Form PDEs from the following equations by eliminating arbitrary functions.
(iv) z = xf1 (x + t) + f2 (x + t)
y
(v) z = f x
3.1 Method 1
∂ 2z 1
+ 9x2 y 2 − cos(2x − y) = f (y) (6)
∂x∂y 2
[For the constant of integration, we use arbitrary functions of variable
held fixed]
Integrating the equation (1) with respect to x (keeping y fixed)
∂z 1
+ 3x3 y 2 − sin(2x − y) = xf (y) + g(y) (7)
∂y 4
Integrate equation (2) with respect to y (keeping x fixed)
Z Z
1
z + x3 y 3 − cos(2x − y) = x f (y)dy + g(y)dy + w(x) (8)
4
R R
Then we substitute f (y)dy = u(y) and g(y)dy = v(y)
1
z + x3 y 3 − cos(2x − y) = xu(y) + v(y) + w(x) (9)
4
1
⇒ z= cos(2x − y) − x3 y 3 + xu(y) + v(y) + w(x)
4
where u, v and w are arbitrary function.
∂ 2z
Example 7. Solve = xy
∂x2
Solution: Integrate the given PDE with respect to x (keeping y fixed)
∂z x2
= y + f (y) (10)
∂x 2
Integrate equation (5) with respect to x (keeping y fixed)
x3
Z
z = y + f (y) dx + g(y) (11)
6
1
⇒ z = x3 y + xf (y) + g(y)
6
where f and g are arbitrary function.
∂ 2z
Example 8. Solve 2
+ z = 0, given that when x = 0, z = ey and
∂x
∂z
=1
∂x
Soution: Here z is a function of x only. So we can treat the given PDE
as ordinary differential equation (ODE) but constants can depend on y.
⇒ (D2 + 1)z = 0
Auxiliary equation is m2 + 1 = 0 ⇒ m = i, −i
∂ 2z 1
= sin (2x + 3y) + f (x) (17)
∂x2 3
Integrating the equation (12) with respect to x
Z
∂z 1
= − cos (2x + 3y) + f (x)dx + g(y) (18)
∂x 6
Integrating the equation (13) with respect to x
Z Z Z
1
z = − sin (2x + 3y) + f (x)dxdx + g(y) dx + h(y) (19)
12
1
z = − sin (2x + 3y) + u(x) + xg(y) + h(y) (20)
12
RR
where u(x) = f (x)dxdx is an arbitrary function of x
Exercise
Solve the following equations by direct integration:
∂ 2z x
(i) = +a
∂x∂y y
∂ 2z 2 ∂z ∂z
(ii) = a z given that when x = 0, = a siny and =0
∂x2 ∂x ∂y
∂ 2z
(iii) = e−t cosx
∂x ∂t
∂ 2z ∂z
(iv) 2
= z, given that when y = 0, z = e x
and = e−x
∂y ∂y
∂ 2z
(v) = x2 y
∂x∂y
3.2 Method 2
Definition 3.3. If the unknown function and its partial derivatives ap-
pear to the power of 1 in a PDE then it is called First Order Linear
Partial Differential Equation. Otherwise it is called non-linear PDE.
∂z ∂z ∂ 2z ∂ 2z ∂ 2z
= p, = q, = r, = s, =t
∂x ∂y ∂x2 ∂x∂y ∂y 2
Example 10. (i) px + yq = 3z- First order linear PDE
Case 1. Grouping
We may be able to group two fractions from auxiliary equation and cancel
the common terms to get an equation in only two variables. On integrat-
ing the differential equation in only two variables by any methods, we
shall obtain one of the relations u in the general solution. This method
may be repeated to give another relation v with the help od another two
fractions in auxiliary equation.
y2z
Example 11. Solve p + xzq = y 2
x
Solution:-
y2z
p + xzq = y 2 (21)
x
Equation (16) × x ⇒ y zp + x2 zq = xy 2
2
Step 1 P = y 2 z, Q = x2 z, R = xy 2
dx dy dz dx dy dz
Step 2 Auxiliary equation- = = ⇒ 2 = 2 = 2
P Q R y z xz xy
Step 3 Solve the auxiliary equation by grouping
dx dy
Taking first two fractions 2 = 2
y z xz
⇒ x2 dx = y 2 dy
x3 y3 a
On integrating = + ⇒ x3 − y 3 = a
3 3 3
3 3
Set u = x − y = a
dx dz
Taking first and third fractions 2 = 2
y z xy
⇒ xdx = zdz
x2 z2 b
On integrating = + ⇒ x2 − z 2 = b
2 2 2
Set v = x2 − z 2 = b
Solution:-
Result
R R cos x
cotx dx = dx = log (sinx) + c
sin x
Step 1 P = z, Q = z, R = z 2 + (x + y)2
dx dy dz dx dy dz
Step 2 Auxiliary equation- = = ⇒ = = 2
P Q R z −z z + (x + y)2
Step 3 Solve the auxiliary equation by grouping
dx dy
Taking first two fractions =
z −z
⇒ dx = −dy
On integrating ⇒ x = −y + a = x + y = a
set u = x + y = a
dx dz
Taking first and third fractions = 2
z z + (x + y)2
put x + y = a in the above equation
dx dz
⇒ = 2
z z + a2
z dz 1 2z dz
⇒ dx = 2 ⇒ dx =
z + a2 2 z 2 + a2
Result Z 0
f (x)
dx = log f (x) + c
f (x)
1 0
On integrating ⇒ x = log (z 2 + a2 ) + b
2
⇒ 2x = log (z 2 + a2 ) + b
⇒ 2x − log (z 2 + (x + y)2 ) = b [since x + y = a]
set u = 2x − log (z 2 + (x + y)2 ) = b
Solution:-
Solution:-
dx − dy dy − dz
⇒ =
x2 − y 2 − yz + zx y 2 − z 2 − zx + xy
dx − dy dy − dz
⇒ =
(x + y)(x − y) + z(x − y) (y + z)(y − z) + x(y − z)
dx − dy dy − dz
⇒ =
(x + y)(x − y) + z(x − y) (y + z)(y − z) + x(y − z)
dx − dy dy − dz
⇒ =
(x − y)(x + y + z) (y − z)(x + y + z)
dx − dy dy − dz
⇒ =
(x − y) (y − z)
On integrating, log (x − y) = log (y − z) + log a
⇒ log (x − y) − log (y − z) = log a
x−y x−y
⇒ log = log a ⇒ =a
y−z y−z
x−y
Set u = =a
y−z
Taking the difference between last two and last and first
fractions and equating them
dy − dz dz − dx
=
y 2 − zx − (z 2 − xy) z 2 − xy − (x2 − yz)
dy − dz dz − dx
⇒ 2 =
y − z 2 − zx + xy z 2 − x2 − xy + yz
dy − dz dz − dx
⇒ =
(y + z)(y − z) + x(y − z) (z + x)(z − x) + y(z − x)
dy − dz dz − dx
⇒ =
(y − z)(x + y + z) (z − x)(x + y + z)
dy − dz dz − dx
⇒ =
(y − z) (z − x)
On integrating, log (y − z) = log(z − x) + log b
⇒ log (y − z) − log(z − x) = log b
y−z y−z
⇒ log = log b ⇒ =b
z−x z−x
y−z
Set v = =b
z−x
Step 4 General
equationis φ(u, v) = 0 or u
= f (v)
x−y y−z x−y y−z
i.e. φ , = 0 or =f
y−z z−x y−z z−x
Example 16. Solve y 2 p − xy q = x(z − 2y)
⇒ d(zy) = 2ydy
y2
On integrating zy = 2 + a ⇒ zy − y 2 = a
2
Step 4 General equation is φ(u, v) = 0 or u = f (v)
i.e. φ(x2 + y 2 , zy − y 2 ) = 0 or x2 + y 2 = f (zy − y 2 )
Exercise
(i) xp + yq = 3z
√ √
(ii) p x + q y = z
(iii) y 2 zp + x2 zq = y 2 x
(v) xp − yq = y 2 − x2
(vi) p − q = log (x + y)
(iii) The aboe procedure may be repeated for the second solution
of auxiliary equation v = b
Solution:
dx dy dz
⇒ = =
x(y − z) y(z − x) z(x − y)
Step 3 Solve the auxiliary equation by using Lagrange’s multipliers
1 1 1
Choose P1 = , Q1 = and R1 = such that
x y z
1 1 1
P P1 + QQ1 + RR1 = x(y − z) + y(z − x) + z(x − y)
x y z
= y−z+z−x+x−y
= 0
P2 dx + Q2 dy + R2 dz = 0 ⇒ dx + dy + dz = 0
On integrating, x + y + z = b
⇒ set v = x + y + z = b
dx dy dz
⇒ = =
x(y 2 − z 2 ) y(z 2 − x2 ) z(x2 − y 2 )
Step 3 Solve the auxiliary equation by using Lagrange’s multipliers
1 1 1
Choose P1 = , Q1 = and R1 = such that
x y z
1 1 1
P P1 + QQ1 + RR1 = x(y 2 − z 2 ) + y(z 2 − x2 ) + z(x2 − y 2 )
x y z
= y 2 − z 2 + z 2 − x2 + x2 − y 2
= 0
So the solution u = a is the solution of the differential equation
1 1 1
P1 dx + Q1 dy + R1 dz = 0 ⇒ dx + dy + dz = 0
x y z
On integrating, log x + log y + log z = log a
⇒ log (xyz) = log a ⇒ xyz = a
Set u = xyz = a
Choose P2 = x, Q2 = y and R2 = z such that
P2 dx + Q2 dy + R2 dz = 0 ⇒ x dx + y dy + z dz = 0
x2 y 2 z 2 a
On integrating, + + =
2 2 2 2
⇒ x + y + z = b ⇒ set v = x2 + y 2 + z 2 = b
2 2 2
Solution:
dx dy dz
⇒ = =
x2 (y − z) y 2 (z − x) z 2 (x − y)
Step 3 Solve the auxiliary equation by using Lagrange’s multipliers
1 1 1
Choose P1 = , Q1 = and R1 = such that
x y z
1 1 1
P P1 + QQ1 + RR1 = x2 (y − z) + y 2 (z − x) + z 2 (x − y)
x y z
= x(y − z) + y(z − x) + z(x − y)
= xy − xz + yz − yx + zx − zy
= 0
Result
−1
Z
d 1 1 1
= 2 then 2
dx = − + c
dx x x x x
1 1 1
On integrating, −
− − = −b
x y z
1 1 1 1 1 1
⇒ + + = b ⇒ set v = + + = b
x y z x y z
Step 4 General
equation isφ(u, v) = 0 or u =
f (v)
1 1 1 1 1 1
i.e. φ xyz, + + = 0 or xyz = f + +
x y z x y z
∂z ∂z
Example 20. Solve (mz − ny) + (nx − lz) = ly − mx
∂x ∂y
Solution:-
P1 dx + Q1 dy + R1 dz = 0 ⇒ x dx + y dy + z dz = 0
x2 y 2 z 2 a
On integrating, + + =
2 2 2 2
⇒ x + y + z = a ⇒ set u = x2 + y 2 + z 2 = a
2 2 2
Choose P2 = l, Q2 = m and R2 = n such that
On integrating, lx + my + nz = b ⇒ set v = lx + my + nz = b
Exercise
(ii) z(x + y)p + z(x − y)q = x2 + y 2 (Hint: (x, −y, z) and (y, x, −z) are
the multipliers)
We can solve auxiliary equation of a linear PDE by using case 1 and case
2 for u = a and v = b.
Example 21. Solve (z 2 − 2yz − y 2 )p + (xy + zx)q = xy − zx
Solution:-
dx dy dz
⇒ = =
z 2 − 2yz − y 2 xy + zx xy − zx
Step 3 Solve the auxiliary equation
Calculation of u = a by using Lagrange’s multipliers
Choose P1 = x, Q1 = y and R1 = z such that
P1 dx + Q1 dy + R1 dz = 0 ⇒ x dx + y dy + z dz = 0
x2 y 2 z 2 a
On integrating, + + =
2 2 2 2
⇒ x2 + y 2 + z 2 = a ⇒ set u = x2 + y 2 + z 2 = a
Calculation of v = b by grouping[Taking last two fractions]
dy dz
=
xy + zx xy − zx
dy dz dy dz
⇒ = ⇒ =
x(y + z) x(y − z) (y + z) (y − z)
Exercise
(iii) x2 p + y 2 q = (x + y)z
3.3 Method 3
CHARPIT’S METHOD
This is a general method for solving non linear first order partial differ-
ential equation f (x, y, z, p, q) = 0 with two independent variables. The
complete solution of such an equation contains only two arbitrary con-
stant. (equal to the number of independent variables x and y)
Working Rule for the Solution of f (x, y, z, p, q) = 0 in Charpit’s Method
Step 5 The simplest relation along with the given PDE to determine
p and q
Solution:-
Solution:-
Result Z 0
f (x)
dx = logf (x) + x
f (x)
3.4 Method 4
∂z 0 ∂z 0 ∂ 2z 00
If z = XY , then = X Y, = XY , 2
=X Y
∂x ∂y ∂x
0
Case 2 Consider the ordinary differential equation Y + aY = 0 and
the auxiliary equation is m + a = 0 and m = −a.
⇒ y = c3 e−ay
The required solution is z = XY
√ √
⇒ z = c1 e (1+ 1+a)x
+ c2 e (1− 1+a)x
c3 e−ay
√ √
⇒ z = c1 c3 e (1+ 1+a)x
+ c2 c3 e (1− 1+a)x
e−ay
√ √
⇒ z = k1 e (1+ 1+a)x
+ k2 e (1− 1+a)x
e−ay
where c1 c2 = k1 and c2 c3 = k2
∂u
Example 25. Using the method of separation of variables, solve =
∂x
∂u
2 + u where u(x, 0) = 6e−3x .
∂t
∂u ∂u
Solution:- =2 +u
∂x ∂t
Assume the solution of given Partial differential equation is u = XT
where X is a function of x alone and T that of t alone.
∂u 0 ∂u 0
⇒ = X T, = XT
∂x ∂t
Substitute these values in given PDE
0 0
⇒ X T = 2XT + XT
0 0
⇒ X T − XT = 2XT
0 0
⇒ (X − X)T = 2XT
0 0
X −X T
⇒ =
2X T
Since x and y are independent variables, it can only be true if each side
is equal to the same constant, a (say)
0 0
X −X T
⇒ = a and =a
0
2X T 0
⇒ X − X = 2aX and T − aT = 0
0 0
⇒ X − X − 2aX = 0 and T − aT = 0
0 0
⇒ X − (2a + 1)X = 0 and T − aT = 0
0
Case 1 Consider the ordinary differential equation X − (2a + 1)X = 0
and the auxiliary equation is m − (2a + 1) = 0 ⇒ m = 2a + 1
⇒ X = c1 e(2a+1)x
0
Case 2 Consider the ordinary differential equation T − aT = 0 and the
auxiliary equation is m − a = 0 ⇒ m = a
⇒ T = c2 eat
The required solution is u = XT ⇒ u = c1 e(2a+1)x c2 eat
⇒ u = c1 c2 e(2a+1)x eat ⇒ u = ke(2a+1)x eat where c1 c2 = k
Given that u(x, 0) = 6e−3x [i.e. at t = 0, u = 6e−3x ]
⇒ 6e−3x = ke(2a+1)x ⇒ k = 6 and 2a + 1 = −3 ⇒ a = −2
⇒ T he required solution is u = 6e−3x e−2t ⇒ u = 6e−(3x+2t)
References