Differential Equations-1
Differential Equations-1
Differential Equations-1
(Semester - I)
US01CMTH02
(Differential Equations -1)
Prepared by Nilesh Y. Patel
Head,Mathematics Department,V.P.and R.P.T.P.Science College
US01CMTH02
UNIT-4
1. Differential equations
1.1. Definition. An equation involving derivatives of a dependent variable with respect
to one or more independent variables and is called a differential equation.
we have to take the square of both the sides, which in turn, will give rise the highest power
dy
10 of the differential coefficient dx .
Exact Differential Equations.
Let M and N be functions of x and y. A differential equation of the form
Mdx + Ndy = 0, (1.1.1)
is said to be exact if the expression on the left hand side of (1.1.1) can be obtained directly
by differentiating some function of x and y.
1.2. Theorem. The necessary and sufficient condition for the differential equation
Mdx + Ndy = 0 (1.2.1)
to be exact is that
∂M ∂N
= . (1.2.2)
∂y ∂x
Proof. Necessity: Suppose that (1.2.1) is exact. So, Mdx + Ndy can be obtained
directly by differentiating some function f = f (x, y). Thus,
∂f ∂f
d[f (x, y)] = Mdx + Ndy ⇒ dx + dy = Mdx + Ndy.
∂x ∂y
So,
∂f ∂f
= M and = N;
∂x ∂y
∂ 2f ∂M ∂ 2f ∂N
= and = . (1.2.3)
∂y∂x ∂y ∂x∂y ∂x
Since, we assume the function to be many times continuously differentiable, (see Proposi-
∂2f ∂2f
tion ??), ∂y∂x = ∂x∂y . As a result, (1.2.3) gives ∂M
∂y
= ∂N
∂x
.
∫ ∂2P
Sufficiency. Let P = Mdx. Then ∂P ∂x (
= M. Hence ∂y∂x = ∂M . This together with
2 2
) ∂y
(1.2.2) gives, ∂N
∂x
= ∂M ∂y
∂ P
= ∂y∂x ∂ P
= ∂x∂y = ∂x∂ ∂P
∂y
. This, on integrating with respect to x,
gives,
∂P
N= + φ(y),
∂y
where φ is a function of y only. Thus we have,
[ ]
∂P ∂P
Mdx + Ndy = dx + + φ(y) dy
∂x ∂y
∂P ∂P
= dx + dy + φ(y)dy
∂x ∂y
= dP + d(F(y)) (where d(F(y)) = φ(y)dy)
= d[P + F(y)],
which shows that Mdx + Ndy = 0, is an exact differential equation.
A working rule for solving an exact differential equation
1. D.E. 3
Let us assume that the differential equation (1.1.1) is exact. The following algorithm
describes a working rule to solve the equation.
(1) Integrate M with respect to x regarding y as a constant.
(2) Integrate with respect to y, the terms in N not involving x.
(3) Add the two expressions obtained in the above two steps and equate the result to an
arbitrary constant. This gives the required solution.
1.3. Example. Solve (x2 − 2xy − y 2 )dx − (x + y)2 dy = 0.
Solution. Here M = x2 − 2xy − y 2 and N = −(x + y)2 = −x2 − 2xy − y 2 . Hence,
∂M ∂N
= −2x − 2y and = −2x − 2y.
∂y ∂x
Thus ∂M∂y
= ∂N
∂x
. Hence the given equation is an exact equation. Integrating M with respect
to x treating y as constant, we obtain
∫ ∫
x3
Mdx = (x2 − 2xy − y 2 )dx = − x2 y − y 2 x. (1.3.1)
3
Also, there is only one term −y 2 in N that does not involve x which, on integration, gives
∫
y3
−y 2 dy = − . (1.3.2)
3
Adding the right hand side terms of (1.3.1) and (1.3.2) and equating the result with a
3 3
constant gives the solution x3 − x2 y − y 2 x − y3 = C, where C is constant of integration.
xdy−ydx
1.4. Example. Solve xdx + ydy + x2 +y 2
= 0.
Solution. First of all we have to simplify the given equation to express it in the form of
(1.1.1).
xdy − ydx
xdx + ydy + =0
x2 + y 2
[ ] [ ]
y x
⇒ x− 2 dx + y + 2 dy = 0.
x + y2 x + y2
Thus,
y x
M=x− and N=y+
x2 + y 2 x2 + y 2
∂M y 2 − x2 ∂N y 2 − x2
⇒ = 2 and = 2 .
∂y (x + y 2 )2 ∂x (x + y 2 )2
Hence the given equation is exact. Now integrating M with respect to x regarding y as
constant, we get,
∫ ∫ [ ]
y
Mdx = x− 2 dx
x + y2
∫ ∫
1
= xdx − y dx
x + y2
2
4
( )
x2 1 −1 x
= − y tan
2 y y
2
( )
x x
= − tan−1 . (1.4.1)
2 y
x
Also, N = y + x2 +y 2
. The only term in N not involving x is y, on integrating it with
y2
respect to y gives 2 . Adding this term to the last term of (1.4.1) and equating the sum
2
to the constant gives the solution x2 − tan−1 xy + y2 = C. That is, x2 − 2 tan−1 xy + y 2 = C,
2
2. Introduction
dy dy
Usually dx is denoted by p in differential equations which involve dx in degree greater
th
than one. Thus the general form of a first order n degree differential equation is
If xp + x + y = 0, then x dx dy
+ y = −x, that is, dx
dy
+ x1 y = −1, which is linear and its
∫ dx
integrating factor is e x = elog x = x. And its solution is
∫
y × I.F. = −1 × I.F.dx + C2
∫
⇒ yx = − xdx + C2
x2
⇒ yx = − + C2
2
⇒ 2xy + x2 − C = 0, (C = 2C2 ). (3.1.2)
dy
If p + 2x = 0, then dx
+ 2x = 0, that is, dy + 2xdx = 0. Hence,
y + x2 − C = 0. (3.1.3)
6
∫
⇒ ye = −(xe − ex dx) + C
x x
⇒ yex = −xex + ex + C
⇒ y = −x + 1 + Ce−x
⇒ y + x − 1 − Ce−x = 0. (3.2.2)
Hence from (3.2.1) and (3.2.2), the required solution is
(y − Cex )(y + x − 1 − Ce−x ) = 0.
3.3. Example. Solve p2 + 2py cot x = y 2 .
Solution. Given equation can be written as
p2 + (2y cot x)p − y 2 = 0.
Solving this for p, we get,
√
−2y cot x ± 4y 2 cot2 x + 4y 2
p=
2
⇒ p = −y cot x ± y cosec x
4. Solvable for y 7
( )
− cos x ± 1
⇒p=y
sin x
( ) ( )
− cos x + 1 − cos x − 1
⇒p=y or p = y
sin x sin x
2 x 2 x
2y sin 2 2y cos 2
⇒p= x x or p = −
2 sin 2 cos 2 2 sin x2 cos x2
x
⇒ p = y tan or p = −y cot x2 .
2
If p = y tan x2 , then
dy
= tan x2 dx
y
( )
⇒ log y = 2 log sec x2 + log C1
⇒ y = C1 sec2 x
2
⇒ y cos2 x
2
= C1
⇒ y(1 + cos x) = 2C1 = C, (say)
⇒ y(1 + cos x) − C = 0.
If p = −y cot x2 , then
dy
= − cot x2 dx
y
( )
⇒ log y = −2 log sin x2 + log C2
C2
⇒y=
sin2 x2
⇒ y sin2 x
2
= C2
⇒ y(1 − cos x) = 2C2 = C, (say)
⇒ y(1 − cos x) − C = 0.
So, the required solution is (y(1 + cos x) − C)(y(1 − cos x) − C) = 0.
npn+1
⇒ xp2 = − +C
n+1
npn−1 C
⇒x=− + 2. (4.1.1)
n+1 p
Substituting this value of x in the given differential equation we get,
−2pnpn−1 2pC
y= + 2 + pn
n+1 p
2C n − 1 n
⇒y= − p . (4.1.2)
p n+1
The required solution is obtained by eliminating p from (4.1.1) and (4.1.2).
p cos p = sin p − y
4. Solvable for y 9
√
1 − cos2 p − y
⇒p=
cos p
√
−1 1 − (C − x)2 − y
⇒ cos (C − x) = (as cos p = C − x)
C −x
[√ ]
1 − (C − x)2 − y
⇒ C − x = cos ,
C −x
which is the required solution.
4.3. Example. Solve y = yp2 + 2px.
Solution. The given equation can be rewritten as
2px
y(1 − p2 ) = 2px or y = . (4.3.1)
1 − p2
Differentiating this with respect to x we get,
(1 − p2 )(2p + 2x(dp/dx)) − 2px(−2p)(dp/dx)
p=
(1 − p2 )2
dp
⇒ p(1 − p2 )2 = 2p(1 − p2 ) + [2x(1 − p2 ) + 4p2 x]
dx
dp
⇒ p(1 − p2 )(1 − p2 − 2) = 2x(1 − p2 + 2p2 )
dx
dp
⇒ − p(1 − p2 )(1 + p2 ) = 2x(1 + p2 )
dx
dp
⇒ p(p2 − 1) = 2x
dx
2dp dx
⇒ =
p(p2 − 1) x
[ ]
1 1 2 dx
⇒ + − dp =
p−1 p+1 p x
⇒ log(p − 1) + log(p + 1) − 2 log p = log x + log C
p2 − 1
⇒ = Cx
p2
⇒ p2 − 1 = p2 Cx
⇒ p2 (1 − Cx) = 1
1
⇒ p2 =
1 − Cx
1
⇒p= √ .
1 − Cx
Substituting this value of p in (4.3.1), we get,
√ √
2x/ 1 − Cx 2x 1 − Cx
y= =
1 − 1−Cx
1
−Cx
√
⇒ 2x 1 − Cx + Cxy = 0,
10
5. Clairaut’s equation
The equation
y = px + f (p) (5.0.1)
is known as Clairaut’s equation. To solve it, we differentiate it with respect to x and get,
dp dp
p=p+x + f ′ (p)
dx dx
dp
⇒ [x + f ′ (p)] =0
dx
dp
⇒ = 0 (hence p = C, a constant) (5.0.2)
dx
or
x + f ′ (p) = 0. (5.0.3)
Now eliminating p from (5.0.1) and (5.0.2) gives
y = Cx + f (C), (5.0.4)
as a solution of (5.0.1). Hence the solution of the Clairaut’s equation is obtained on
replacing p by a constant C. If we eliminate p between (5.0.1) and (5.0.3), we get a solution
which does not contain any arbitrary constant and is not a particular case of (5.0.4). This
type of solution is known as singular solution. Some equations can be reduced to Clairaut’s
form by a suitable substitution.
√ √
Solution. Let x2 = u and y 2 = v, that is,√x = u and y = v. Hence dx = 12 u−1/2 du
and dy = 21 v −1/2 dv. This gives p = dx
dy
= √uv du
dv
. Substituting these values in the given
equation we get,
( ) √
√ u dv 2 u dv √
uv − (u + v − 1) + uv = 0
v du v du
( )2
dv dv
⇒u − (u + v − 1) + v = 0
du du
dv
⇒ uP 2 − (u + v − 1)P + v = 0, where P =
du
⇒ u(P − P ) − v(P − 1) + P = 0
2
⇒ v(P − 1) = uP (P − 1) + P
P
⇒ v = Pu + , (5.3.1)
P −1
which is in Clairaut’s form. Hence by substituting P = C, in (5.3.1), we get the solution as
C
v = Cu + C−1 . Now substituting back the values of u and v, we get the required solution
2 2
y = Cx + C−1 .C
5.4. Example. Solve y 2 (y − px) = x4 p2 .
Solution. Let x = u1 , y = v1 . Then dx = −1 du, dy = −1 dy 2
u2 v2
dv, p = dx = uv2 du
dv
. Substituting
these values in the given equation we get,
[ ] ( )2
1 1 u2 dv 1 1 u4 dv
− = 4 4
v 2 v v 2 du u u v du
( )2
1 u dv 1 dv
⇒ 3− 4 = 4
v v du v du
( )2
dv dv
⇒v−u =
du du
dv
⇒ v = P u + P 2 where P = ,
du
which is in Clairaut’s form and so its solution is v = Cu + C 2 . After substituting the
values of u and v back, the solution becomes C 2 xy + Cy − x = 0.
5.5. Example. Solve e2x (p − 1) + p3 e2y e−x = 0.
ey dy dv
Solution. Let ex = u and ey = v. Then ex dx = du and ey dy = dv. Hence ex dx
= du
,
i.e., p = uv du
dv
. Substituting these values in the given equation, we get,
[ ] [ ]3
2 u dv u dv v 2
u −1 + =0
v du v du u
( )3
u3 dv u2 dv
⇒ −u +
2
=0
v du v du
( )3
dv dv
⇒u −v+ =0
du du
12
dv
⇒ v = P u + P 3, where P = ,
du
which is in Clairaut’s form. Its solution is v = Cu + C 3 or ey = Cex + C 3 .
Hence the differential equation of the family of the orthogonal trajectories is φ(X, Y, − dX dY
)=
0. In usual notations, we see that the differential equation of the family of orthogonal tra-
dy
jectories of the given family of φ(x, y, dx ) = 0 is φ(x, y, − dx
dy
) = 0, so that it is obtained on
replacing dx by − dx
dy
dy
.
6.3. Example. Find the orthogonal trajectories of the semi-cubical parabolas ay 2 = x3 ,
where a is the variable parameter.
Solution. Differentiating the given equation with respect to x, we get
dy dy dy dy
2ay = 3x2 ⇒ 2ay 2 = 3x2 y ⇒ 2x3 = 3x2 y ⇒ 2x = 3y, (6.3.1)
dx dx dx dx
which is the differential equation of the given family. Putting − dx
dy
dy
in place of dx in (6.3.1),
we get,
dx 3
2x(− ) = 3y, ⇒ 2xdx + 3ydy = 0 ⇒ x2 + y 2 = C,
dy 2
which the equation of the family of the orthogonal trajectories of the given semi-cubical
parabolas.
Note: The following example deals with the asymptote. The topic is presently out of the
scope of this text book. If the students have not developed the theory of asymptotes and
the curve tracing, this example may safely, skipped.
6.4. Example. Find the orthogonal trajectories of family of parabolas
y 2 = 4a(x + a), (6.4.1)
where a is the parameter.
Solution. Differentiating (6.4.1) with respect to x, we have
dy
2y = 4a.
dx
Putting this value in (6.4.1), we get,
[ ] ( )2
dy y dy dy dy
2
y = 2y x+ ⇒ y = 2x + y , (6.4.2)
dx 2 dx dx dx
which is the differential equation of the given family of parabolas. Putting dx dy
= − dx
dy
in
(6.4.2), we get the following differential equation of the required orthogonal trajectories of
the given family.
( )2 ( )2
dx dx dy dy
y = −2x + y − ⇒y + 2x − y = 0,
dy dy dx dx
which is same as the differential equation (6.4.2) of the given family of parabolas. Hence
the given family (6.4.1) is self orthogonal. That is, the orthogonal trajectories of the system
belong to the system itself. Hence the required equation of the orthogonal trajectories of
the given family is y 2 = 4C(x + C), where C is the parameter.
♣♣♣♣♣♣♣♣