Module 2 - First Order ODE
Module 2 - First Order ODE
Equations
Introduction
In this module, we will discuss several methods to find solutions of first-order differential
equations. After finishing this module, students are expected to be familiar with the different
forms of a first-order differential equation and how to find solutions for them.
Topic Outcomes
2. Show the solution of a variable separable DE using the correct method, expressed in
its simplest form.
4. Show the solution of a DE with homogeneous coefficients using the correct method,
expressed in its simplest form.
6. Show the solution of an exact DE using the correct method, expressed in its simplest
form.
8. Show the solution of a linear DE using the correct method, expressed in its simplest
form.
9. Recognize that a 1st order DE has an integrating factor by grouping terms to give
exact correct differentials.
10. Show solution of a DE by using integrating factors, expressed in its simplest form.
11. Recognize that a non-exact DE can be made exact by means of the use of a correct
integrating factor.
13. Recognize that terms in a DE can be used to make substitution by making appropriate
substitution.
14. Show solutions of a DE by making substitutions out of its terms to make use of other
methods.
16. Show the solution of a Bernoulli DE using the correct method, expressed in its simplest
form.
17. Recognize that coefficients in a DE are linear in two variables by finding the intersection
of the lines or making appropriate substitution correctly.
18. Show the solution of a DE with coefficients linear in two variables using the correct
method, expressed in its simplest form.
Separation of Variables
then the equation is said be variable separable. The solution can be immediately obtained
by integrating both sides of (2).
Z Z
A (x) dx + B (y) dy = c
Another definition:
A first-order differential equation in the form
dy
= g(x)h(y)
dx
is said to have separable variable.
Example 1.
dy
(1 − x) = y2
dx
multiplying both sides by dx,
(1 − x)dy = y 2 dx
dy dx
=
y2 1−x
Integrating both sides of the equation gives,
Z Z
dy dx
2
=
y 1−x
1
− = − ln |1 − x| + c
y
1
− = − ln |1 − x| + ln c
y
Simplifying gives the final result as
y ln |c(1 − x)| = 1
Example 2.
xdx dy
+ 3 =0
ex2 y
2 dy
xe−x dx + 3 = 0
y
Integrating both sides of the equation,
Z Z
−x2 dy
xe dx + =c
y3
1 2 1
− e−x − y −2 = c
2 2
2
e−x + y −2 = −2c
2
e−x + y −2 = c
Example 4.
1 1
(2y dx) = (3x dy)
xy xy
2 dx 3 dy
=
x y
Integrating both sides,
2 dx 3 dy
=
x y
Z Z
dx dy
2 =3 +c
x y
2 ln |x| = 3 ln |y| + ln |c|
ln |x2 | = ln |y 3 | + ln |c|
ln |x2 | = ln |cy 3 |
x2 = cy 3
Example 5.
1 2x
2x
1
2x
ye dx = (4 + e )dy
y(4 + e ) y(4 + e2x )
e2x dx dy
2x
=
4+e y
Integrating both sides, we have
e2x dx
Z Z
dy
= +c
4 + e2x y
1
ln |4 + e2x | = ln |y| + ln |c|
2
ln |4 + e2x | = 2 ln |y| + 2 ln |c|
ln |4 + e2x | = ln |y 2 | + ln |c2 |
ln |4 + e2x | = ln |c2 y 2 |
4 + e2x = c2 y 2
Notice that c2 can be written as just c, since they are both arbitrary, but we should
note that the left side of the equation is always greater than 0. Therefore, if we just write c,
we will be restricted to values of c ≥ 0.
Example 6.
dr = b cos θ dr + br sin θ dθ
dr − b cos θ dr = br sin θ dθ
(1 − b cos θ) dr = br sin θ dθ
dr b sin θdθ
=
r 1 − b cos θ
Integrating both sides of the equation,
Z Z
dr b sin θdθ
= +c
r 1 − b cos θ
ln |r| = ln |1 − b cos θ| + ln |c|
ln |r| = ln |c(1 − b cos θ)|
r = c(1 − b cos θ)
Example 7.
x dx dy
− =0
x+2 y
xdx (y + 1)dy
− =0
e−x y3
xex dx − y −3 (y + 1)dy = 0
Solve the particular solution given the initial condition equation 2xyy 0 = 1 + y 2 , when x = 2
and y = 3.
Solution:
Divide the given equation by x(1 + y 2 ),
2xyy 0 = 1 + y 2
dy
2xy = 1 + y2
dx
2y dy dx
2
=
1+y x
Integrate both sides of equation
Z Z
2ydy dx
2
=
1+y x
2
ln |1 + y | = ln |x| + ln |c|
ln |1 + y 2 | = ln c|x|
1 + y 2 = cx
1 + y 2 = cx
1 + 32 = c(2)
c=5
1 + y 2 = 5x
Equations with Homogeneous Coefficients
Polynomials in which all the terms are of the same degree, such as
x2 + 2xy + y 2 x3 + y 3 x4 y + 8y 5
are called homogeneous polynomials.
Let us extend the concept of homogeneity and define what a homogeneous function is.
Definition
The function f (x, y) is said to be homogeneous of degree k in x and y if, and only if,
f (λx, λy) = λk f (x, y)
Theorem 0.0.1. If M (x, y) and N (x, y) are both homogeneous and of the same degree, then
M (x, y)
the function is homogeneous of degree zero.
N (x, y)
Example 1.
Determine if the function f (x, y) is homogeneous and if it is, state the degree of the function.
y x4
f (x, y) = 2y 3 e x −
x + 3y
Solution:
Original function
y x4
f (x, y) = 2y 3 e x −
x + 3y
f (λx, λy),
(λx)4 x4
λy y
3 3 3 x
f (λx, λy) = 2(λy) e λx − = λ · 2y e −
λx + 3λy x + 3y
Now,
f (λx, λy) = λ3 f (x, y)
Therefore, f (x, y) is homogeneous of degree 3.
Equations with Homogeneous Coefficients
Suppose the coefficients M and N in the differential equation of the first order
M (x, y) + N (x, y) = 0
are both homogeneous and are of same degree in x and y. We can put the equation in the
form
dy M (x, y)
+ = 0.
dx N (x, y)
By Theorems (1) and (2), we can put it in the form
dy y
+f =0
dx x
From this we may introduce a new variable v such that y = vx. Then the equation
becomes
dv
x + v + f (v) = 0
dx
in which the variable are now separable.
The method would have been equally successful if we used x = vy to obtain an equation
in y and v. However, it is sometimes easier to substitute for the variable whose differential
has the simpler coefficients.
Example 2.
Solve the equation
2(2x2 + y 2 )dx − xydy = 0
Solution:
From M dx + N dy = 0, both M and N are homogeneous functions of degree 2, let us put
y = vx
dy = vdx + xdv
Substitute to the original equation,
2(2x2 + v 2 x2 )dx − x2 v(vdx + xdv) = 0
(4x2 + 2v 2 x2 − v 2 x2 )dx − vx3 dv = 0
(4x2 + v 2 x2 )dx − vx3 dv = 0
x2 (4 + v 2 )dx − vx3 dv = 0
By separation of variables, divide the equation by (4 + v 2 )x3 ,
dx vdv
− =0
Z x Z 4 + v2
dx vdv
− =c
x 4 + v2
1
ln |x| − ln (4 + v 2 ) = ln |c|
2
2 ln |x| − ln (v 2 + 4) = ln |c|
2
x
ln 2 = ln |c|
v + 4
2
ln x 2
e 4+v = eln |c|
x2
=c
4 + v2
x2 = c(4 + v 2 )
y
Since v = ,
x
2
2 y
x =c 4+
x
2
y
x2 =c 4+
x2
2
4x + y 2
x2 =c
x2
x4 = c(4x2 + y 2 )
Example 3.
x = vy
dx = vdy + ydv
x
Since v = , then
y
2
x2
6 ln |y| − = ln |c|
y2
6y 2 ln |y| − x2 = y 2 ln |c|
6y 2 ln |y| − y 2 ln |c| = x2
y 2 (ln (y 6 ) − ln |c|) = x2
6
y
x2 = y 2 ln
c
x = y ln |cy 6 |
2 2
Example 4.
y
Since v = , then
x
x
cos v = ln
c
y x
cos = ln
x c
y
cos = ln |cx|
x
Example 5.
Solution:
Since M and N are homogeneous function of degree 1, then let
y = vx
dy = vdx + xdv
xdx + x2 ln |v|dv = 0
dx
+ ln |v| = 0
Z xZ
dx
+ ln |v| = 0
x
ln |x| + v ln |v| − v = ln |c|
y
Since v = , then
x
ln |x| + v ln |v| − v = ln |c|
y y y
ln |x| + ln − = ln |c|
x x x
y
x ln |x| + y ln − y = cx
x
x ln |x| + y ln |y| − y ln |x| − y = cx
(x − y) ln |x| + y ln |y| = cx + y
Example 6.
Solve the equation
y 2 dy = x(xdy − ydx)ex/y
Solution:
Express the given equation as
y 2 dy = x(xdy − ydx)ex/y
(y 2 − x2 ex/y )dy = −xyex/y dx
Since both coefficients are homogeneous and of degree 2, let
x
x = vy; v=
y
dx = vdy + ydv
x
Since v = , then
y
ln |y| − ln |c| = ev (v − 1)
y x/y x
ln = e
−1
c y
x/y x − y
y
ln = e
c y
y
y ln = (x − y)ex/y
c
y ln |cy| = (x − y)ex/y
Example 7.
Solve the equation
2y
xdx + sin (ydx − xdy) = 0
x
Solution:
Rewrite the given equation as
2 y 2 y
x + y sin dx − x sin dy = 0
x x
Let
y
y = vx; v=
x
dy = vdx + xdv
x + vx sin v dx − x sin2 v(vdx + xdv) = 0
2
x2 (v + 1)2 = c(v − 2)
y
Since v = , then
x
y 2 y
x2 +1 =c −2
x x
2
x(y + x) = c(y − 2x)
Example 9.
Solution:
Let
y
y = vx; v=
x
dy = vdx + xdv
2vt3 dt + t4 (v 2 + 1)dv = 0
dt v 2 + 1
2 + dv = 0
t v
dt dv
2 + vdv + =0
t v
1
2 ln |t| + v 2 + ln |v| = ln |c|
2
s
Since v = , then
t
2
1 s s
2 ln |t| + + ln = ln |c|
2 t t
2
1 s
2 ln |t| + + ln |s| − ln |t| = ln |c|
2 t2
s2
ln |s| + ln |t| − ln |c| + 2 = 0
2t2
st s
ln + 2 = 0
c 2t
st
2t2 ln + s2 = 0
c
2 2
st
s = −2t ln
c
2 2
s = −2t ln |cst|
Exact Equations
∂f ∂f
dx + dy = 0 (2)
∂x ∂y
Definition
A differential expression M (x, y)dx + N (x, y)dy is an exact differential in a region R of
the xy-plane if it corresponds to the differential of some function f (x, y) defined in R. A
first-order differential equation in the form
is said to be an exact equation if the expression on the left-hand side is an exact differential.
∂M ∂N
=
∂y ∂x
Method of Solution
Given an equation M (x, y)dx + N (x, y)dy = 0, which is exact. Then there exist a function
F (x, y) for which
∂F
= M (x, y)
∂x
Solving for F (x, y), Z
F (x, y) = M (x, y)dx + T (y)
∂F
where the arbitrary function T (y) is the constant of integration. And assuming = N (x, y)
∂y
Z
∂
M (x, y) dx + T 0 (y) = N (x, y)
∂y
Giving us Z
0 ∂
T (y) = N (x, y) − M (x, y) dx
∂y
Finally integrating to find T (y) with respect to y.
F (x, y) = c
and Z
0 ∂
T (x) = M (x, y) − N (x, y) dy
∂x
The implicit solution is still
F (x, y) = c
Example 1.
Solve the equation
(x + 2y)dx + (2x + y)dy = 0
Solution:
Test for exactness,
∂M ∂N
= 2; =2
∂y ∂x
∂M ∂N
Since = , therefore the equation is exact. Hence, its solution is F = c, where
∂y ∂x
∂F ∂F
= M; =N
∂x ∂y
Using
∂F
= M = x + 2y
∂x
Solving for F by integrating both sides with respect to x, holding y as constant,
Z Z
∂F = (x + 2y)∂x
x2
F = + 2xy + T (y)
2
To solve for T (x), find the derivative of F with respect to y
∂F
= 2x + T 0 (y)
∂y
∂F
Since = N , then
∂y
2x + T 0 (y) = 2x + y
T 0 (y) = y
Z Z
0
T (y)dy = ydy
y2
T (y) =
2
Therefore,
x2 y2
F = + 2xy +
2 2
And the set of solutions of the given equation is defined by
x2 y2
+ 2xy + =c or x2 + 4xy + y 2 = c
2 2
Example 2.
∂F = x2 + 2y∂y
2y 2
F = x2 y + + T (x)
2
F = x2 y + y 2 + T (x)
∂F
Solve for T (x), find the derivative of F with respect to x that is, = 2xy − T 0 (x)
∂x
∂F
Since = N then,
∂x
2xy + T 0 (x) = 2xy − 3x2
T 0 (x) = −3x2
T (x) = −x3
x2 y + y 2 − x3 = C
Example 3.
Solve the equation
Solution:
∂M ∂N
= −2 sin (2y) − 6x2 y = −2 sin (2y) − 6x2 y
∂y ∂x
Hence, the equation is exact and its solution is F = c, where
∂F
= cos (2y) − 3x2 y 2
∂x
and
∂F
= cos (2y) − 2x sin (2y) − 2x3 y
∂y
Solving for F ,
∂F
= cos (2y) − 3x2 y 2
∂x
∂F = cos (2y) − 3x2 y 2 ∂x
∂F
= x2 − x sec2 y + T 0 (y)
∂y
∂F
Since = N , then
∂y
Therefore,
F = x2 y − x tan y + c1
and the set of solutions of the given equation is defined by
x2 y − x tan y = c
Example 5.
Solve the equation
cos x cos y − cot x dx − sin x sin ydy = 0
Solution:
∂M ∂N
= −x cos x sin y = −x cos x sin y
∂y ∂x
Hence, the equation is exact and its solution is F = c, where
∂F ∂F
= cos x cos y − cot x and = − sin x sin y
∂x ∂y
Solving for F using
∂F
=N
∂y
we have
∂F
= − sin x sin y
∂y
∂F = − sin x sin y∂y
F = sin x cos y + T (x)
∂F
Solving for T (y) find , we have
∂x
∂F
= cos x cos y + T 0 (x)
∂x
∂F
Since = M , then
∂x
cos x cos y + T 0 (x) = cos x cos y − cot x
T 0 (x) = − cot x
T (x) = − ln | sin x|
Therefore,
or
sin x cos y = ln |c sin x|
Example 6.
Solve the equation
r + sin θ − cos θ dr + r sin θ + cos θ dθ = 0
Solution:
∂M ∂N
= cos θ + sin θ = cos θ + sin θ
∂θ ∂r
Hence, the equation is exact and its solution is F = c, where
∂F ∂F
= r sin θ + cos θ and = r + sin θ − cos θ
∂θ ∂r
Solving for F using
∂F
=N
∂θ
∂F = r sin θ + cos θ ∂θ
F = r − cos θ + sin θ + T (r)
∂F
Solving for T (r) obtain ,
∂r
∂F
= − cos θ + sin θ + T 0 (r)
∂r
∂F
Since = M , then
∂r
∂F
Since = M , then
∂x
y cos (xy) + T 0 (x) = 2x + y cos (xy)
T 0 (x) = 2x
T (x) = x2
Therefore,
F = sin (xy) + x2 = c
and the set of solutions of the given equation is defined by
sin (xy) + x2 = c
Example 8.
Solve the equation
2 2
2x 3x + y − ye−x dx + x2 + 3y 2 + e−x dy = 0
Solution:
∂M 2 ∂N 2
= 2x − 2xe−x = 2x − 2xe−x
∂y ∂x
∂M ∂N
Since = , therefore the equation is exact. Hence, its solution is F = c, where
∂y ∂x
∂F 2
= M = 2x 3x + y − ye−x
∂x
∂F 2
= N = x2 + 3y 2 + e−x
∂y
Solving for F using
∂F
=N
∂y
then
∂F 2
= x2 + 3y 2 + e−x
∂y
2
∂F = (x2 + 3y 2 + e−x )∂y
2
F = x2 y + y 3 + yex + T (x)
Solving for T (x) get the derivative of F with respect to x we have,
∂F 2
= 2xy − 2xye−x + T 0 (x)
∂x
∂F
Since = M then,
∂x
2 2
2xy − 2xye−x + T 0 (x) = 2x 3x + y − ye−x
2 2
2xy − 2xye−x + T 0 (x) = 6x2 + 2xy − 2xye−x
T 0 (x) = 6x2
T (x) = 2x3
Therefore,
2
F = x2 y + y 3 + yex + 2x3
and the set of solutions of the given equation is defined by
2
x2 y + y 3 + yex + 2x3 = c
Example 9.
Solve the equation
3y(x2 − 1)dx + (x3 + 8y − 3x)dy = 0
when x = 0 and y = 1.
Solution:
∂M ∂N
= 3(x2 − 1) = 3(x2 − 1)
∂y ∂x
∂M ∂N
Since = , therefore the equation is exact. Hence, its solution is F = c, where
∂y ∂x
∂F ∂F
= M = 3y(x2 − 1) = N = x3 + 8y − 3x
∂x ∂y
Solving for F using
∂F
=M
∂x
then
∂F
= 3y(x2 − 1)
∂x
∂F = 3y(x2 − 1)∂x
3
x
F = 3y − x + T (y)
3
Solving for T (x), find the derivative of F with respect of y
3
∂F x
=3 − x + T 0 (y)
∂y 3
∂F
Since = N , then
∂y
3
x
3 − x + T 0 (y) = x3 + 8y − 3x
3
T 0 (y) = 8y
T (y) = 4y 2
Therefore,
x3
F = 3y − x + 4y 2
3
and the set of solutions of the given equation is defined by
3
x
3y − x + 4y 2 = c
3
x3 y − 3xy + 4y 2 = c
When x = 0 and y = 1
x3 y − 3xy + 4y 2 = c
c=4
x3 y − 3xy + 4y 2 = 4
x(x2 − 3) = 4(1 − y 2 )
Linear Equation of Order One
Standard Form
By dividing both sides of the linear equation by a1 (x), we obtain a more useful form,
the standard form, of a linear equation,
dx
+ P (x)y = Q(x)
dy
Suppose a linear equation in the standard form, there exist an integrating factor v(x) >
0, a function of x alone. Then the equation
dy
v(x) + P (x)y = v(x)Q(x)
dx
must be an exact equation. The equation can be written in the form
M (x, y) dx + N (x, y) dy = 0
with
M (x, y) = v(x)P (x)y − v(x)Q(x)
in which v(x), P (x) and Q(x) are functions of x alone.
P (x)dx dy
R R R
P (x)dx P (x)dx
e + P (x)ye = Q(x)e
dx
The left member of the equation is actually the derivative of the function
R
P (x)dx
ye
P (x)dx dy
R R R
P (x)dx P (x)dx
e + P (x)ye = Q(x)e
dx
is exact.
From the previous discussion, after determining the integrating factor v(x), we can
now find the solution of a first-order linear equation
dy
+ P (x)y = Q(x)
dx
and from
P (x)dx dy
R R R
P (x)dx P (x)dx
e + P (x)ye = Q(x)e
dx
From which, we have the general solution
R
Z R
P (x)dx
ye = Q(x)e P (x)dx dx + c
Example 1.
y − 1 + x cot x = c csc x
y + x cot x + c csc x = 1
Example 4.
y sin x = x + c
Example 5.
It has been found out that integrating factor is an aid in finding a solution of linear
equations. At present we are concerned enough to enable us to find integrating factors by
inspection. The ability to do this depends largely upon the recognition of certain common
exact differentials and upon experience.
d(xy) = x dy + y dx (1)
x y dx − x dy
d = (2)
y y2
y x dy − y dx
d = (3)
x x2
x y dx − x dy
d Tan−1 = (4)
y x2 + y 2
In this lecture, some equations when terms are reduced to exact differentials becomes
a more recognizable differential equation. Hence, solving them becomes easier.
Example 1.
Solve the equation
y(2xy + 1)dx − xdy = 0
Solution:
y(2xy + 1)dx − xdy = 0
2xy 2 dx + ydx − xdy = 0
1 2 1
2
2xy dx + ydx − xdy = 0 · 2
y y
ydx − xdy
+ 2xdx = 0
y2
x
d + 2xdx = 0
y
Z Z
x
d + 2xdx = c
y
x
+ x2 = c
y
x + x2 y =cy
x(1 + xy) = cy
Example 2.
Solve the equation
y(y 3 − x)dx + x(y 3 + x)dy = 0
Solution:
y(y 3 − x)dx + x(y 3 + x)dy = 0
y 3 (ydx + xdy) − x(ydx − xdy) = 0
x ydx − xdy
ydx + xdy − =0
y y2
x ydx − xdy
d(xy) − =0
y y2
Z
x ydx − xdy
Z
d(xy) − =c
y y2
( xy )2
xy − =c
2
x2
2xy − 2 = c
y
2xy − x2 = cy 2
3
Example 3.
Solve the equation
(x3 y 3 + 1)dx + x4 y 2 dy = 0
Solution:
(x3 y 3 + 1)dx + x4 y 2 dy = 0
x3 y 2 (xdy + ydx) + dx = 0
1 3 2 1
x y (xdy + ydx) + dx = 0 ·
x x
dx
x2 y 2 (xdy + ydx) + =0
Z Z x
dx
(xy)2 (xdy + ydx) + =c
x
(xy)3
+ ln |x| = ln |c|
3
x3 y 3 + 3 ln |x| = 3 ln |c|
x3 y 3 = −3 ln |cx|
Example 4.
Solve the equation
x4 (ydx + xdy) + y 2 (xdy − ydx) = 0
Solution:
x4 (ydx + xdy) + y 2 (xdy − ydx) = 0
1 4 2 1
4
x (ydx + xdy) + y (xdy − ydx) = 0 · 4
x x
2
y xdy − ydx
(xdy + ydx) + 2 =0
x x2
y 2 y
d(xy) + d =0
x x
Z Z 2
y y
d(xy) + d =c
x x
1 y 3
xy + · =c
3 x
y3
3xy + 3 = c
x
3x y + y 3 = cx3
4
y(3x4 + y 2 ) = cx3
Example 5.
Example 7.
Example 9.
Recall from Linear Equation, the form y 0 + P (x)y = Q(x) can be transformed into an
exact equation when we multiply the equation by an integrating factor. The same basic idea
sometimes apply for a non-exact differential equation M (x, y)dx + N (x, y)dy = 0. That is,
it is sometimes possible to find an integrating factor µ(x, y) so that after multiplying, the
left-hand side of
µ(x, y)M (x, y) dx + µ(x, y)N (x, y) dy = 0 (1)
is an exact differential. In an attempt to find µ(x, y), we return to the criterion for exactness.
The equation (1) is exact if and only if
∂ ∂
µ(x, y)M (x, y) = µ(x, y)N (x, y)
∂y ∂x
By the product rule of differentiation (subscript notation means partial derivatives)
µMy + M µy = µNx + N µx
or
µx N − µy M = (My − Nx ) µ (2)
Although M , N , My and Nx are known functions of x and y, the difficulty here in determining
the unknown µ(x, y) from (2) is that we must solve a partial differential equation. Since we
are not yet prepared for that, we make simplifying assumptions.
Suppose µ is a function of only one variable, e.g., µ depends only on x. In this case
dµ
µx = and µy = 0
dx
so (2) can be written as
dµ My − Nx
= ·µ (3)
dx N
My − Nx
We still at an impasse of the quotient depends on both x and y. However,
N
My − Nx
if all obvious algebraic simplifications are made, the quotient turns out to depend
N
solely on the variable x, then (3) is a first-order ordinary differential equation. We can finally
determine µ because (3) is separable and linear. It follows that
R My −Nx
dx
µ(x) = e N
In the like manner, it follows from (2) if µ depends only on the variable y, then
dµ My − Nx
=− ·µ (4)
dy M
My − Nx
In this case, if − is a function of y alone, then we can solve (4) for µ to be
M
R My −Nx
− dy
µ(y) = e M
M (x, y) dx + N (x, y) dy = 0
My − Nx R My −Nx
dx
• If is a function of x alone, the integrating factor is e N .
N
My − Nx R My −Nx
• If − is a function of y alone, the integrating factor is e − M dy .
M
Example 1.
Solve the equation
(x2 + y 2 + 1)dx + x(x − 2y)dy = 0 (1)
Solution:
Let,
M (x, y) = x2 + y 2 + 1
My = 2y
Therefore,
−y 2 x−1 − x−1 + x + y = c
−y 2 x−1 − x−1 + x + y = c
or
−x2 + y 2 − xy + 1 = cx
Example 2.
Solve the equation
2y(x2 − y + x)dx + (x2 − 2y)dy = 0 (1)
Solution:
Let
M (x, y) = 2y(x2 − y + x)
My = 2x4 − 4y + 2x
N (x, y) = x2 − 2y
Nx = 2x
My − Nx 2x4 − 4y + 2x − 2x
=
N x2 − 2y
My − Nx 2(x2 − 2y)
=
N x2 − 2y
My − Nx
=2
N
My − Nx
= µ(x)
N
Solving for the integrating factor, R
2dx
e = e2x
Multiply the integrating factor to (1)
Simplifying,
2ye2x (x2 − y + x)dx + e2x (x2 − 2y)dy = 0 (2)
Test for exactness of Equation (2),
Therefore,
x2 ye2x − y 2 e2x + c = c
x2 ye2x − y 2 e2x = c
or
y(x2 − y) = ce−2x
Example 3.
Solve the equation
y(2x − y + 1)dx + x(3x − 4y + 3)dy = 0 (1)
Solution:
Let
M (x, y) = y(2x − y + 1)
∂M
= 2x − 2y + 1
∂y
N (x, y) = x(3x − 4y + 3)
∂N
= 6x − 4y + 3
∂x
Using the formula,
My − Nx 2x − 2y + 1 − (6x − 4y + 3)
=
M y(2x − y + 1)
My − Nx −2(2x − y + 1)
=
M y(2x − y + 1)
My − Nx 2
=−
M y
My − Nx
= µ(y)
M
Solving for the integrating factor,
− y2 dx
R
e− = e2 ln y
− y2 dx
R
e− = y2
Simplifying,
y 3 (2x − y + 1)dx + xy 2 (3x − 4y + 3)dy = 0]y 2 (2)
Test for exactness of Equation (2),
M (x, y) = y 3 (2x − y + 1)
∂M
= 6xy 2 − 4y 3 + 3y 2
∂y
N (x, y) = xy 2 (3x − 4y + 3)
∂N
= 6xy 2 − 4y 3 + 3y 2
∂x
∂M ∂N ∂F
Since = , therefore (2) is exact. Find the solution F = c, using = M,
∂y ∂x ∂x
Z Z
∂F = y 3 (2x − y + 1)∂x
F = x2 y 3 − xy 4 + xy 3 + T (y)
∂F
= 3x2 y 2 − 4xy 3 + 3xy 2 + T 0 (y)
∂y
∂F
Since = N (x, y), then
∂y
Therefore,
F = x2 y 3 − xy 4 + xy 3 + T (y)
F = x2 y 3 − xy 4 + xy 3 + c
x2 y 3 − xy 4 + xy 3 + c = c
x2 y 3 − xy 4 + xy 3 = c
or
xy 3 (x − y + 1) = c
Example 4.
Solve the equation
y(4x + y)dx − 2(x2 − y)dy = 0 (1)
Solution:
Let
M (x, y) = y(4x + y)
∂M
= 4x + 2y
∂y
N (x, y) = 2(x2 − y)
∂N
= −4x
∂x
Using the formula,
My − Nx 4x + 2y − (−4x)
=
M y(4x + y)
My − Nx 2(4x + y)
=
M y(4x + y)
My − Nx 2
=
M y
My − Nx
= µ(y)
M
Solving for the integrating factor,
2
R
e− y
dx
= e−2 ln y
2
R
e− y
dx
= y −2
Simplifying,
y −1 (4x + y)dx − 2y −2 (x2 − y)dy = 0 (2)
Test for exactness of (2),
M (x, y) = y −1 (4x + y)
∂M 4x
=− 2
∂y y
N (x, y) = −2y −2 (x2 − y)
∂N 4x
=− 2
∂x y
∂M ∂N ∂F
Since = , therefore (2) is exact. Find the solution F = c, using = M,
∂y ∂x ∂x
Z Z
∂F = y −1 (4x + y)∂x
2x2
F = + x + T (y)
y
∂F 2x
= − 2 + T 0 (y)
∂y y
∂F
Since = N (x, y), then
∂y
2x2
− + T 0 (y) = −2y −2 (x2 − y)
y2
2x2 2x2 2
− 2 + T 0 (y) = − 2 +
y y y
2
T 0 (y) =
y
T (y) = 2 ln |y|
T (y) = ln (y 2 )
Therefore,
2x2
F = + x + T (y)
y
2x2
F = + x + ln (y 2 )
y
2x2
+ x + ln (y 2 ) = c
y
2x2
+ x + ln (y 2 ) = c
y
or
2x2 + xy + y ln (y 2 ) = cy
Example 5.
Let
M (x, y) = (xy + 1)
∂M
=x
∂y
N (x, y) = x(x + 4y − 2)
∂N
= 2x + 4y − 2
∂x
Using the formula,
My − Nx x − 2x + 4y − 2
=
N x(x + 4y − 2)
My − Nx −(x + 4y − 2)
=
N x(x + 4y − 2)
My − Nx 1
=−
N x
My − Nx
= µ(x)
N
Solving for the integrating factor,
− x1 dx
R
e = e− ln |x|
− x1 dx
R
e = x−1
Simplifying,
x−1 (xy + 1)dx + (x + 4y − 2)dy = 0 (2)
Test for exactness of (2),
x + T 0 (y) = x + 4y − 2
T 0 (y) = 4y − 2
T (y) = 2y 2 − 2y
Therefore,
F = xy + ln |x| + T (y)
F = xy + ln |x| + 2y 2 − 2y
xy + ln |x| + 2y 2 − 2y = C
Example 6.
Solution:
Let
M (x, y) = 2y 2 + 3xy − 2y + 6x
∂M
= 4y + 3x − 2
∂y
N (x, y) = x(x + 2y − 1)
∂N
= 2x + 2y − 1
∂x
Using the formula,
My − Nx 4y + 3x − 2 − (2x + 2y − 1)
=
N x(x + 2y − 1)
My − Nx x + 2y − 1
=
N x(x + 2y − 1)
My − Nx 1
=
N x
My − Nx
= µ(x)
N
Solving for the integrating factor,
1
R
e x
dx
= eln |x|
1
R
dx
e x =x
F = x3 y + x2 y 2 − x2 y + T (x)
∂F
= 3x2 y + 2xy 2 − 2xy + T 0 (x)
∂x
∂F
Since = M (x, y), then
∂x
3x2 y + 2xy 2 − 2xy + T 0 (x) = x(2y 2 + 3xy − 2y + 6x)
3x2 y + 2xy 2 − 2xy + T 0 (x) = 2xy 2 + 3x2 y − 2xy + 6x2
T 0 (x) = 6x2
T (x) = 2x3
Therefore,
F = x3 y + x2 y 2 − x2 y + T (x)
F = x3 y + x2 y 2 − x2 y + 2x3
x3 y + x2 y 2 − x2 y + 2x3 = c
x2 (xy + y 2 − y + 2x) = c
Example 7.
Solution:
Let
M (x, y) = y(y + 2x − 2)
∂M
= 2y + 2x − 2
∂y
N (x, y) = −2(x + y)
∂N
= −2
∂x
Using the formula,
My − Nx 2y + 2x − 2 − (−2)
=
N −2(x + y)
My − Nx 2(x + y)
=
N −2(x + y)
My − Nx
= −1
N
My − Nx
= µ(x)
N
Solving for the integrating factor,
R
−dx
e = e−x
Simplifying,
e−x y(y + 2x − 2)dx − 2e−x (x + y)dy = 0 (2)
Test for exactness of (2),
M (x, y) = ye−x (y + 2x − 2)
∂M
= 2ye−x + 2xe−x − 2e−x
∂y
N (x, y) = −2e−x (x + y)
∂N
= 2xe−x − 2e−x + 2ye−x
∂x
∂M 0 ∂N 0 ∂F
Since = , therefore (2) is exact. Find the solution F = c, using = M,
∂y ∂x ∂x
Z Z
∂F = ye−x (y + 2x − 2)∂x
Therefore,
−y 2 e−x − 2xye−x + c = c
−y 2 e−x − 2xye−x = c
−ye−x (y + 2x) = c
y(y + 2x) = cex
Example 8.
Solve the equation
y 2 dx + (3xy + y 2 − 1)dy = 0 (1)
Solution:
Let
M (x, y) = y 2
∂M
= 2y
∂y
N (x, y) = 3xy + y 2 − 1
∂N
= 3y
∂x
Using the formula,
My − Nx 2y − 3y
=
M y2
My − Nx 1
=−
M y
My − Nx
= µ(y)
M
Solving for the integrating factor,
− y1 dx
R
e− = eln |y|
− y1 dx
R
e− =y
Multiply the integrating factor to (1)
y[y 2 dx + (3xy + y 2 − 1)dy] = 0 · y
Simplifying,
y 3 dx + y(3xy + y 2 − 1)dy = 0 (2)
Test for exactness of (2),
M (x, y) = y 3
∂M
= 3y 2
∂y
N (x, y) = y(3xy + y 2 − 1)
∂N
= 3y 2
∂x
∂M ∂N ∂F
Since = , therefore (2) is exact. Find the solution F = c, using = M,
∂y ∂x ∂x
Z Z
∂F = y 3 ∂x
F = xy 3 + T (y)
∂F
= 3xy 2 + T 0 (y)
∂y
∂F
Since = N (x, y), then
∂y
F = xy 3 + T (y)
y4 y2
F = xy 3 + −
4 2
And the set of solutions of (1) is defined by
y4 y2
xy 3 + − =c
4 2
y4 y2
xy 3 + − =c
4 2
y 2 (4xy + y 2 − 2) = 4c
y 2 (4xy + y 2 − 2) = c
Example 9.
Solution:
Let
M (x, y) = 2y(x + y + 2)
∂M
= 2x + 4y + 4
∂y
N (x, y) = y 2 − x2 − 4x − 1
∂N
= −2x − 4
∂x
Using the formula,
My − Nx 2x + 4y + 4 − (−2x − 4)
=
M 2y(x + y + 2)
My − Nx 4(x + y + 2)
=
M 2y(x + y + 2)
My − Nx 2
=
M y
My − Nx
= µ(y)
M
Solving for the integrating factor,
2
R
e− y
dx
= e−2 ln |y| = y −2
Simplifying,
2y −1 (x + y + 2)dx + y −2 (y 2 − x2 − 4x − 1)dy = (2)
Test for exactness of (2),
M (x, y) = 2y −1 (x + y + 2)
∂M
= −2xy −2 − 4y −2
∂y
N (x, y) = y −2 (y 2 − x2 − 4x − 1)
∂N
= −2xy −2 − 4y −2
∂x
∂M ∂N ∂F
Since = , therefore (2) is exact. Find the solution F = c, using = M,
∂y ∂x ∂x
Z Z
∂F = 2y −1 (x + y + 2)∂x
F = x2 y −1 + 2x + 4xy −1 + T (y)
∂F
= −x2 y 2 − 4xy −2 + T 0 (y)
∂y
∂F
Since = N (x, y), then
∂y
F = x2 y −1 + 2x + 4xy −1 + T (y)
1
F = x2 y −1 + 2x + 4xy −1 + y +
y
An equation in the form M (x, y)dx + N (x, y)dy = 0 may not yield at once (or at all)
to the methods discussed in the previous discussions. Even then the usefulness of these
methods will not be exhausted. It may be possible by some change of variables to transform
the equation into a type which we know how to solve.
The term (x + 2y) appears twice in the equation and thus it attracts attention, hence we
may put
x + 2y = u
and, because no other terms of x and ystands out, we retain any one of the two.
Another equation,
(1 + 3x sin y)dx − x2 cos y dy = 0
the presence of both sin y and its differential cos y dy, and the fact that the variable y appears
in the equation in no other manner, leads us to put
u = sin y du = cos y dy
u = Ax + By + C
wherein B 6= 0.
Example 1.
Solve
(3 tan x − 2 cos y) sec2 xdx + tan x sin ydy = 0 (1)
Solution:
Let
u = tan x
du = sec2 xdx
M (x, y) = 3u − 2 cos y
My = 2 sin y
and
N (x, y) = u sin y
Nx = sin y
Then,
u3 − u2 cos y = c (3)
Since u = tan x, then (3) becomes
u = sin y
du = cos ydy
Simplifying,
x2 v(3 + v)dx + 2x3 dv = 0 (3)
Using separation of variables, divide (3) by x3 v(3 + v)
dx 2dv
+ =0
x v(v + 3)
Z Z
dx 2dv
+ =c
x v(v + 3)
2 2
ln |x| + ln |v| − ln |v + 3| = ln |c|
3 3
3 ln |x| + 2 ln |v| − 2 ln |v + 3| = 3 ln |c|
ln (u2 ) + ln |x3 | − ln (u + 3x)2 = ln |c|
u 2 x3
=c
(u + 3x)2
x3 sin2 y
=c
(3x + sin y)2
x3 sin4 y = c(3x + sin y)2
Example 4.
Solve
dy
= (9x + 4y + 1)2 (1)
dx
Solution:
We can write (1) as
dy = (9x + 4y + 1)2 dx (2)
Let
u = 9x + 4y + 1
du = 9dx + 4dy
1
dx = (du − 4dy)
9
Substituting these to (2), we have
1
dy = u2 (du − 4dy)
9
9dy = u2 (du − 4dy)
(9 + 4u2 )dy = u2 du
Separating the variables, gives
1 2 2 1
(9 + 4u )dy = u du ·
4u2 + 9 4u2 + 9
u2 du
dy = 2
4u + 9
u2 du
Z Z
dy =
4u2 + 9
1 9 2 −1 2u
y = u− Tan +c
4 16 3 3
We have,
1 3 2u
y = u − Tan−1 +c (3)
4 8 3
Since u = 9x + 4y + 1, then (3) becomes
1 3 −1 2(9x + 4y + 1)
y = (9x + 4y + 1) − Tan +c
4 8 3
−1 2(9x + 4y + 1) 18x + 2
Tan = +c
3 3
2
(9x + 4y + 1) = tan (6x + c)
3
2(9x + 4y + 1) = 3 tan (6x + c)
Example 5.
dy
= 1 + 6xex−y (1)
dx
Solution:
Let
u=x−y
du = dx − dy
dy
= 1 + 6xeu
dx
dx − du
= 1 + 6xeu
dx
dx − du = (1 + 6xeu )dx
[1 − (1 + 6xeu )]dx = du
We have,
−6xeu dx = du (2)
Diving (2) by eu to separate the variables becomes
du
−6xdx = u
Z Ze
−6xdx = e−u du
−3x2 = −e−u + c
−3x2 = −e−(x−y) + c
−3x2 = −ey−x + c
3x2 + c = ey−x
Example 6.
Solve
(x + 2y − 1)dx − (x + 2y − 5)dy = 0 (1)
Solution:
Let
u = x + 2y
du = dx + 2dy
dx = du − 2dy
We now have,
u − 4 ln |u + 3| − y = C (2)
Since u = x + 2y, then substituting to (2) we have
Solve,
(x + 2y − 1)dx + [2(x + 2y) − 3]dy = 0 (1)
Solution:
Let
u = x + 2y
du = dx + 2dy
dx = du − 2dy
We now have,
1 2
(u) − u − y = c (2)
2
Since u = x + 2y, then (2) becomes
1
(x + 2y)2 − (x + 2y) − y = c
2
(x + 2y)2 − 2(x + 2y) − 2y = c
(x + 2y)2 − 2x − 6y = c
(x + 2y − 1)2 = 2y + c
Example 8.
Solve,
y(x tan x + ln |y|)dx + tan xdy = 0 (1)
Solution:
We can rewrite (1) as
dy
(x tan x + ln |y|)dx + tan x =0 (2)
y
Let
u = ln |y|
dy
du =
y
M (x, u) = x tan x + u
∂M
=1
∂u
and
N (x, u) = tan x
∂N
= sec2 x
∂x
Using the formula,
Mu − Nx 1 − secx
=
N tan x
Mu − Nx − tan2 x
=
N tan x
Mu − Nx
= − tan x
N
Integrating,
−x cos x + sin x + u sin x = c (4)
Since u = ln |y|, then (4) gives us the result as
or
ln |y| = x cot x + c csc x − 1
Example 9.
Solve
4(3x + y − 2)dx − (3x + y)dy = 0 (1)
Solution:
Let
u = 3x + y
y = u − 3x
dy = du − 3dx
The equation
dy
+ P (x)y = Q(x)y n (1)
dx
may be put in the form
dy
y −n
+ P (x)y 1−n = Q(x) (2)
dx
Now considering y 1−n and its derivative
d dy
y 1−n = (1 − n) y −n
dx dx
Now, let
dz dy
z = y 1−n and = (1 − n) y −n
dx dx
Substituting this in (2), we will have
dz
+ (1 − n) P (x)z = (1 − n) Q(x)
dx
which is a linear equation.
The solution is
R
Z R
(1−n)P (x)dx (1−n)P (x)dx
ze = (1 − n) Q(x)e dx + c
In terms of x and y,
R
Z R
1−n (1−n)P (x)dx (1−n)P (x)dx
y e = (1 − n) Q(x)e dx + c
y −2 x3 = −x + c
x3 = (c − x)y 2
Alternate Solution
2x3 y 0 = y(y 2 + 3x2 )
y(y 2 + 3x2 )
y0 =
2x3
3
0 y 3x2
y = 3 + 3y
2x 2x
0 3 1 3
y − y = 3 ·y
2x 2x
3 dx
y −3 dy − y −2 dx = 3 (1)
2x 2x
Let
z = −y 2
dz = 2y −3 dy
x3 z = x + c (3)
1
Since z = − 2 , then (3)
y
3 1
x − 2 =x+c
y
x3 = −xy 2 + cy 2
x3 = y 2 (c − x)
Example 2.
Solution:
1 ln |y|
n=3 P (y) = − Q(y) = −
2y 2y
We have the reduced polynomials,
1 1 ln |y| ln |y|
xr = x1−3 = x−2 Pr (y) = −2 · − = Qr (y) = −2 · − =
2y y 2y y
The integrating factor is
1
R
dy
e y = eln |y| = y
The general solution is,
Z
−2 ln |y|
x y= ydy + c
y
Z
−2
x y= (ln |y|) dy + c
x−2 y = y ln |y| − y + c
y = x2 y ln |y| − x2 y + cx2
Example 3.
Solution:
y 0 = y − xy 3 e−2x
y 0 − y = −xe−2x y 3
y −2 e2x = x2 + c
e2x = y 2 x2 + c
Example 4.
Solution:
1 2 3 1
6y dx − x(2x + y)dy = 0 ·
6y 2 dy 6y 2 dy
dx 1 1
− x = 2 x4
dy 6y 3y
dx
This is a Bernoulli D.E. of the form + P (y)x = Q(y)xn with
dy
1 1
n=4 P (y) = − Q(y) =
6y 3y 2
The reduced polynomials are
1 1
xr = x−3 Pr (y) = Qr (y) = −
2y y2
The integrating factor is,
1
R 1 1
dy
e 2y = e 2 ln |y| = y 2
The general solution is
Z
−3 1 1 1
x y = 2 − 2 y 2 dy + c
y
Z
1 3
x−3 y 2 = −y − 2 dy + c
1 1
x−3 y 2 = 2y − 2 + c
1
y = x3 2 + cy 2
Example 5.
x2 y 2 = 2y 3 + c
Example 6.
(y 4 − 2xy)dx + 3x2 dy = 0
1 4 2 1
(y − 2xy)dx + 3x dy = 0 ·
3x2 dx 3x2 dx
dy 2 1
− y = − 2 y4
dx 3x 3x
This is a Bernoulli D.E. with
2 1
n=4 P (x) = − Q(x) = −
3x 3x2
The reduced polynomials are
2 1
yr = y −3 Pr (x) = Qr (x) =
x x2
The integrating factor is
2
R
e x
dx
= e2 ln |x| = x2
The general solution is
Z
−3 2 1
y x = x2 dx + c
x2
Z
−3 2
y x = dx + c
y −3 x2 = x + c
x2 = y 3 (x + c)
Example 7.
Solution:
2xyy 0 = y 2 − 2x3
1 0 2 3 1
2xyy = y − 2x
2xy 2xy
y 1
y0 = − x2
2x y
dy y
− = −x2 y −1
dx 2x
This is a Bernoulli D.E. with
1
n = −1 P (x) = − Q(x) = −x2
2x
The reduced polynomials are
1
yr = y 2 Pr (x) = − Qr (x) = −2x2
x
The integrating factor is
R
− x1 dx 1
e = e− ln |x| =
x
The general solution is
y2
Z
2
1
= −2x dx + c
x x
y2
Z
= −2x dx + c
x
y2
= −x2 + c
x
y2 = (c − x2 )x
Alternate Solution
2xyy 0 = y 2 − 2x3
1 0 2 3 1
2xyy = y − 2x
2xy 2xy
y 1
y0 = − x2
2x y
dy y 1
− = − x2
dx 2x y
dy y 1 2
y − = − x y
dx 2x y
We have,
dy 1
y − = −x2 (1)
dx 2x
Let
z = −y 2
dz = −2ydy
zx−1 = x2 + c
y2
− = x2 + c
x
y 2 = x(c − x2 )
Example 8.
Now,
dy 2y x2
y2 + = (1)
dx 3x 3
Let
z = y3
dz = 3y 2 dy
Substitute to (1) becomes
1 dy 2 x2
+ z=
3 dx 3x 3
dy 2
+ z = x2 (2)
dx x
Solving the integrating factor we have,
2dx
R
e x = e2 ln x = x2
Substituting the integrating factor to (2) becomes
dz
x2 + 2xz = x4
dx
x2 dz + 2xzdx = x4 dx
d(x2 z) = x4 dx
Z Z
d(x z) = x4 dx
2
x5
x2 z = +c
5
x5
x2 y 3 = +c
5
5x2 y 3 = x5 + C
Example 9.
y2
Z
x 1
= dx + c
x3 2 x3
y2
Z
1
= dx
x3 2x2
y2 1
=− +c
x3 2x
2y 2 = cx3 − x2
Alternate Solution
(x2 + 6y 2 )dx − 4xydy = 0
1 2 2 1
(x + 6y )dx − 4xydy = 0 ·
−4xydx −4xydx
dy x 3y
− − =0
dx 4y 2x
dy 3y x
y − = y
dx 2x 4y
We have,
dy 3 x
y − y2 = (1)
dx 2x 4
Let
z = −y 2
dz = −2ydy
Substituting to (1) becomes
1 dz 3 x
− + z=
2 dx 2x 4
Now,
dz 3 1
− z=− x (2)
dx x 2
Solving for integrating factor we have
− x3 dx
R
e = e3 ln x = x−3
Substituting the integrating factor to (2), gives us
1
x−3 dz − 3x−4 dx = − x−2 dx
2
1
d(x3 z) = − x−2 dx
Z Z2
1
d(x−3 z) = − x−2 dx
2
1 1
x−3 z = · + c
2 x
y2 1 1
− 3 = · +c
x 2 x
y2 1
− 3 = +c
x 2x
2y 2 = −x2 + cx3
Equations with Coefficients Linear in the Two Variables
a1 x + b1 y + c1 = 0 a2 x + b 2 y + c 2 = 0 (2)
These lines may be parallel or they may intersect. There will not be two lines if a1 = b1 = 0
or a2 = b2 = 0, but equation (1) will then be linear in one of its variables.
will change the equations (2) into equation of line through the origin of the uv -coordinate
system, namely,
a1 u + b 1 v + c 1 = 0 a2 u + b 2 v = 0 (4)
Therefore, since dx = du and dy = dv, the change of variables
x=u+h y =v+k
where (h, k) is the point of intersection of the lines in (2), will transform (1) into
a2 x + b2 y = k(a1 x + b2 y)
Solution:
The lines
y−2=0 and x−y−1=0
intersect at (3, 2). Let
x=u+3
dx = du
and
y =v+2
dy = dv
v = ru
dv = rdu + udr
Alternate Solution
Equation (2) can be written as
vdu − (u − v)dv = 0
du u
− = −1
dv v
This is linear in u, with
1
P (v) = − Q(v) = −1.
v
The integrating factor is
R
− v1 dv −1 1
e = e− ln v = ev = .
v
The general solution is
Z
u 1
= −1 dv + c
v v
u
= − ln |v| − ln |c|
v
u = −v ln |cv|
Since u = x − 3 and v = y − 2, then
x − 3 = (2 − y) ln |c(y − 2)|
Example 2.
Solve the equation
(x − 4y − 9)dx + (4x + y − 2)dy = 0 (1)
Solution:
The lines
x − 4y − 9 = 0 and 4x + y − 2 = 0
intersect at (1, −2). Let
x=u+1 dx = du
and
y =v−2 dy = dv
Then (1) becomes
(u − 4v)du + (4u + v)dv = 0 (2)
Let
v = ru
dv = rdu + udr
substitute v and dv to transform (2) into
(u − 4ru)du + (4u + ru)(rdu + udr) = 0
u(1 + r2 )du + u2 (4 + r)dr = 0
du r+4
+ 2 dr = 0
Z u Z r +1
du r+4
+ dr = c
u r2 + 1
1
ln |u| + ln |r2 + 1| + 4 Tan−1 r = c
2
ln |u2 (r2 + 1)| + 8 Tan−1 r = c
v
Since r = , then
u
2 v 2 + u2
−1 v
ln u
2
+ 8 Tan =c
u u
and since u = x − 1, v = y + 2 then
2
2 v + u2
ln u + 8 Tan−1 v = 0
u2 u
x−1
ln |(y + 2)2 + (x − 1)2 | + 8 Tan−1 =0
y+2
Example 3.
Solution:
The lines
2x − y = 0 and 4x + y − 6 = 0
lines intersect at (1, 2). Let
x=u+1 dx = du
and
y =v+2 dy = dv
Then (1) becomes
(2u − v)du + (4u + v)dv = 0 (2)
Let
v = ru
dv = rdu + udr
(v + u)3
=c
(v + 2u)2
3
(y − 2) + (x − 1)
2 =c
(y − 2) + 2(x − 1)
(x + y − 3)3
=c
(2x + y − 4)2
(x + y − 3)3 = c(2x + y − 4)2
Example 4.
Solution:
The lines
x − 4y − 3 = 0 and x − 6y − 5 = 0
intersect at (−1, −1). Let
x=u−1 dx = du
and
y =v−1 dy = dv
Then (1) becomes
(u − 4v)du − (u − 6v)dv = 0 (2)
Let
v = ru
dv = rdu + udr
Solution:
The lines
2x + 3y − 5 = 0 and 3x − y − 2 = 0
intersect at (1, 1). Let
x=u+1 dx = du
and
y =v+1 dy = dv
Then (1) becomes
(2u + 3v)du + (3u − v)dv = 0 (2)
Let
v = ru
dv = rdu + udr
substitute v and dv to transform (2) into
(2u + 3ru)du + (3u − ru)(rdu + udr) = 0
(2 + 6r − r2 )udu + u2 (3 − r)dr = 0
du 3−r
+ dr = 0
u Z 2 + 6r − r2
r−3
Z
du
+ dr = c
u 2 + 6r − r2
1
ln |u| + ln |2 + 6r − r2 | = ln |c|
2
ln |u2 (2 + 6r − r2 )| = ln |c|
u2 (2 + 6r − r2 ) = c
v
Since r = , then
u
2
2 v v
u 2+6 − =c
u u
2u2 + 6uv − v 2 = c
and since u = x − 1, v = y − 1 then
2(x − 1)2 + 6(x − 1)(y − 1) − (y − 1)2 = c
2x2 − y 2 + 6xy − 10x − 4y = c
Example 6.
Solution:
Let
u=x+y
du = dx + dy
Since u = x + y, then
x − 2u + 3 ln |u + 2| = c
x − 2(x + y) + 3 ln |(x + y) + 2| = c
x + 2y + c = 3 ln |x + y + 2|
Example 7.
Solve the equation
(x − 2)dx + 4(x + y − 1)dy = 0 (1)
Solution:
The lines
x−2=0 and 4(x + y − 1) = 0
intersect at (2, −1). Let
x=u+2 dx = du
and
y =v−1 dy = dv
Then (1) becomes
udu + 4(u + v)dv = 0 (2)
Let
u = rv
du = rdv + vdr
substitute v and dv to transform (2) into
rv(rdv + vdr) + 4(rv + v)dv =0
v(r2 + 4r + 4)dv + rv 2 dr =0
dv rdr
+ =0
v (r + 2)2
Z Z
dv rdr
+ =c
v (r + 2)2
2
ln |v| + ln |r + 2| + = − ln |c|
r+2
2
= − ln cv(r + 2)
r+2
u
Since r = , then
v
2 u
u = − ln cv
+ 2
+2 v
v
2v
= − ln c(u + 2v)
u + 2v
and since u = x − 2, v = y + 1 then
2(y + 1)
= − ln |c[(x − 2) + 2(y + 1)]|
(x − 2) + 2(y + 1)
2(y + 1) = −(x + 2y) ln |c(x + 2y)|
Example 8.
Solve the equation
(x − 1)dx − (3x − 2y − 5)dy = 0 (1)
Solution:
The lines
x−1=0 and 3x − 2y − 5 = 0
intersect at (1, −1). Let
x=u+1 dx = du
and
y =v−1 dy = dv
Then (1) becomes
udu − (3u − 2v)dv = 0 (2)
Let
u = rv du = rdv + vdr
substitute v and dv to transform (2) into
rv(rdv + vdr) − (3rv − 2v)dv = 0
rv 2 dr + (r2 v − 3rv + 2v)dv = 0
rdr dv
2
+ =0
Z r − 3r + 2 Z v
rdr dv
2
+ =c
r − 3r + 2 v
2 ln |r − 2| − ln |r − 1| + ln |v| = ln |c|
v(r − 2)2
ln = ln |c|
r−1
v(r − 2)2
=c
r−1
v(r − 2)2 = c(r − 1)
u
Since r = , then
v
u 2 u
v −2 =c −1
v v
2
(u − 2v) = c(u − v)
and since u = x − 1, v = y + 1 then
2
(x − 1) − 2(y + 1) = c [(x − 1) − (y + 1)]
(x − 2y − 3)2 = c(x − y − 2)
Example 9.
Solve the equation
(2x − 3y + 4)dx + 3(x − 1)dy = 0 (1)
Solution:
The lines 2x − 3y + 4 = 0 and 3(x − 1) = 0 intersect at (1, 2). Let
x=u+1 dx = du
and
y =v+2 dy = dv
Then (1) becomes
(2u − 3v)du + 3udv = 0 (2)
Let
v = ru du = rdv + vdr
substitute v and dv to transform (2) into
v
Since r = , then
u
2 ln |u| + 3r = c
v
2 ln |u| + 3 =c
u
and since u = x − 1, v = y − 2 then
v
2 ln |u| + 3 =c
u
y−2
2 ln |x − 1| + 3 =c
x−1
Bibliography
[Boyce and Diprima, 2001] Boyce, W. and Diprima, R. (2001). Elementary Differential
Equations with Boundary Value Problems. John Wylie & Sons, Inc.
[Bronson and Costa, 2009] Bronson, R. and Costa, G. (2009). Schaum’s Outline of
Differential Equations. McGraw-Hill Education, 3rd edition.
[Lebl, 2019] Lebl, J. (2019). Notes on Diffy Qs: Differential Equations for Engineers.
Independently Published.
[Rainville et al., 2013] Rainville, D., Bedient, P., and Bedient, R. (2013). Elementary
Differential Equations: Pearson New International Edition. Pearson Education, Limited.
[Zill, 2012] Zill, D. (2012). A First Course in Differential Equations with Modeling
Applications. Cengage Learning.