Ode
Ode
Ode
I. Differential equations
Ex. 1
y cos x, (1.1)
y + 4y 0, xdx + ydy 0, (1.2)
2u 2 u
2
c 0 (1.3)
t 2 x 2
2. Partial differential equations: If an equation involves more than one independent variable, so
that partial derivatives enter, it is known as a partial differential
equation. Equation (1.3) is a partial differential equation.
O. D. E. 1
Ex. 2
d2y
x2 + 2y 2 is linear. (1.4)
dx 2
dy
x + y2 is nonlinear. (1.5)
dx
d2y
dx 2 + siny 0 is nonlinear. (1.6)
d2y
dx 2 + y2 0 is quasi-linear. (1.7)
1. Order: The order of a differential equation is the order of the highest derivative.
2. Degree: The degree of a differential equation is the degree of the highest derivative entering,
when the equation has been rationalized and cleared of fractions.
Ex. 3
y + ytanx sin2x is first order and first degree. (1.8)
(y xy)2 k2(1 + y2) is first order and second degree. (1.9)
y2 (1 + y2)3 is second order and second degree. (1.10)
IV. Solutions
A relation, connecting the independent and dependent variables, which satisfies the differential
equation identically, is called a solution of a differential equation.
1. General solution: A solution involving the maximum number of arbitrary constants is called the
general solution.
2. Particular solution: A solution which is derivable from the general solution by assigning fixed
values to the arbitrary constants is called a particular solution.
3. Singular solution: A solution which cannot be obtained from the general solution is called a
singular solution. A singular solution can occur only in the solution of
nonlinear differential equations.
O. D. E. 2
Chapter2 Ordinary Differential Equations of the First Order and
First Degree
M 1 ( x) N ( y)
M 2 ( x)
dx 2
N1 ( y )
dy C
Ex. 1
9yy + 4x 0.
Solution: 9ydy + 4xdx 0 9y2 + 4x2 C
Ex. 2
y 1 + y2.
dy
Solution: dx tan1y x + C y tan(x + C)
1 y2
Ex. 3
y + 5x4y2 0, y(0) 1.
dy 1
Solution: 2
5 x 4 dx x 5 C
y y
1
y(0) 1 1 C The solution is y x5 + 1
Ex. 4
y 2xy, y(0) 1.
dy
Solution: 2xdx lny x2 + C y e x C y C e x
2 2
y
y(0) 1 1 Ce0 C 1 The solution is y e x
2
O. D. E. 3
1. Homogeneous equation
Definition: If f (tx, ty) trf (x, y), then f (x, y) is a homogeneous function of degree r.
Substituting t 1/x into f (tx, ty) trf (x, y), we have
y 1 y y
f (1, ) r f ( x, y ) f ( x, y ) x r f (1, ) x r F ( ) .
x x x x
y
In particular r 0, then f (x, y) F( )
x
If M(x, y)dx + N(x, y)dy 0, when M(x, y) and N(x, y) are of the same degree in x and y,
M/N is a homogeneous function of degree zero, and the differential equation can be written
as
dy M y (2.2)
F( )
dx N x .
Let y ux, equation (2.2) becomes
du dx du
ux F (u ) is separable.
dx x F (u ) u
Ex. 5
2xyy y2 + x2 0.
Solution: Let y ux, the equation becomes
2xux(u + ux) u2x2 + x2 0 2u(u + xu) u2 + 1 0 2xuu + u2 + 1 0
2udu dx
ln(1 + u2) lnx + C 1 + u2 C/x
1 u 2
x
2
1+ x2 + y2 Cx
y C
x x
Ex. 6
y y
(x + ycos )dx xcos dy 0.
x x
dx
Solution: Let y ux (x + uxcosu)dx xcosu(udx + xdu) 0 cosudu
x
y
sinu lnx + C sin lnx C
x
O. D. E. 4
Form: (a1x + b1y + c1)dx + (a2x + b2y + c2)dy 0 (2.3)
a1 b1
(1) If , let x X + , y Y + , equation (2.3) becomes
a 2 b2
[(a1X + b1Y) + (a1 + b1 + c1)]dX + [(a2X + b2Y) + (a2 + b2 + c2)]dY 0. (2.4)
a1 b1 c1 0
We choose
a 2 b2 c2 0
Then equation (2.4) reduces to
(a1X + b1Y) dX + (a2X + b2Y) dY 0
is a homogeneous equation.
a1 b1 c
(2) If k 1 , equation (2.3) becomes
a 2 b2 c2
[k (a2 x + b2 y) + c1] dx + (a2 x + b2 y + c2) dy 0. (2.5)
dv a 2 dx
Let v a2x + b2y dy , equation (2.5) becomes
b2
dv a 2 dx
(kv + c1)dx + (v + c2) 0
b2
a v c2
[kv c1 2 ( v c2 )]dx dv 0 is a separable one.
b2 b2
a1 b1 c1
(3) If k, equation (2.3) becomes
a 2 b2 c2
k(a2x + b2y + c2)dx + (a2x + b2y + c2)dy 0. (2.6)
a. If a2x + b2y + c2 0, we get only a trivial solution.
b. If a2x + b2y + c2 0, equation (2.6) reduces to kdx + dy 0 y kx + C.
Ex. 7
O. D. E. 5
(4x + 3y + 1)dx + (x + y + 1)dy 0.
4 x 3 y 1 0
Solution: x 2, y 3
x y 1 0
Let X x 2, Y y + 3, we have
(4X + 3Y)dX + (X + Y)dY 0
Let Y uX (4X + 3uX)dX + (X + uX)(udX + Xdu) 0
dX 1 u
( 4 4u u 2 )dX X (1 u )du 0 du 0
X (u 2) 2
dX 1 1 1
2
du 0 ln X ln(u 2) C
X u 2 ( u 2) u2
1 x2
ln[ X (u 2)] C ln(2 x y 1) C
u2 2x y 1
Ex. 8
(2x 4y + 5)y + x 2y + 3 0.
1 u
Solution : Let u x 2 y u 1 2 y y
2
1 u
The equation becomes ( 2u 5) u30
2
( 2u 5)(1 u ) 2u 6 0 4u 11 ( 2u 5)u 0
2u 5
( 4u 11)dx ( 2u 5)du 0 dx du 0
4u 11
1 1/ 2 1 1/ 2
dx ( )du 0 dx ( )du C1
2 4u 11 2 4u 11
u 1
x ln(4u 11) C1 8 x 4u ln( 4u 11) 8C1
2 8
8 x 4( x 2 y ) ln[4( x 2 y ) 11] 8C1 4 x 8 y ln( 4 x 8 y 11) C
Ex. 9
(x y 1)dx + (2x 2y 2)dy 0.
Solution: Dividing by x y 1 dx + 2dy 0 x + 2y C
O. D. E. 6
du
a
Let u ax + by + c u a + by y dx , equation (2.7) becomes
b
du
a
du du du
dx f (u ) a bf (u ) a bf (u ) dx is separable.
b dx dx a bf (u )
Ex. 10
y (x + y 7)2. (2.8)
Solution: Let u x + y 7 u 1 + y, equation (2.8) has the result
du 1
u 1 u2 2 dx tan u x + C x + y 7 tan(x + C)
1 u
4. Isobaric equations
O. D. E. 7
Definition: If f (tx, tmy, tm1y) trf (x, y, y), we say that f (x, y, y) is an isobaric function of weight
r.
In particular, if t 1/x, there follows
y y 1
f (1, m , m1 ) r f ( x, y , y ) or
x x x
y y y y
f ( x, y, y ) x r f (1, m , m1 ) x r F ( m , m 1 ). ( 2.9)
x x x x
y y
A differential equation can be put into the form (2.9), then F ( m , m1 ) 0.
x x
y
Solving this for m 1 , and multiplying by xm1, we have
x
y
y x m1( ) (2.10)
xm
y dy du
Let u m
y ux m x m muxm1 , (2.10) is reducible to
x dx dx
du du dx
x m muxm 1 x m1(u ) is separable.
dx (u) mu x
[Note] If the weights m and m 1 are assigned to y and y respectively, the term xaybyc has the
weight a + bm + c(m 1). If an expression is to be isobaric, all its terms must be isobaric
and of the same weight.
Ex. 11
2x3y 1 + 1 4 x 2 y
.
1
Solution: The weights of each term are 3 + (m 1), 0, (0,2 + m), if m 2, every term has the
2
same weight. Let y ux2 y ux2 2ux3, the equation becomes
2x3(x2u 2ux3) 1 + 1 4u 2xu 4u 1 + 1 4u
du dx du dx d ( 1 4u 1) dx
1 4u 1 4u 2 x 1 4u ( 1 4u 1) 2 x 2( 1 4u 1) 2 x
1 4x2 y 1
ln( 1 4u 1) ln x C C
x
Ex. 12
y2 + (1 + xy)y 0.
O. D. E. 8
Solution: weight 2m,m , 1 + m + m , so 2m m m 1.
Let y ux1 y ux1 ux2
x2u2 + (1 + u) (x1u x2u) 0 x2u2 + x1u x2u + x1uu x2u2 0
(1 u )du dx 1 dx
(1 u )u x 1u 1du ln u u ln x C
u x u x
u
ln u C ln y xy C
x
Ex. 13
x3y x2y + y2 0.
Solution: weight 3 + m 1, 2 + m, 2m, so m + 2 2m m 2, Let y ux2, y x2u + 2xu
x3(x2u + 2xu) x2ux2 ux2 xu + 2u u u2
du dx 1 1 dx ux
0 du 0 ln C
u u2 x u u 1 x u 1
y/x xy
C C
y / x 1
2
y x2
u
and from N , we can determine k ( y ).
y
Ex. 14
(x3 + 3xy2) dx + (3x2y + y3) dy 0.
( x 3 3 xy 2 ) (3 x 2 y y 3 )
Solution: Since 6xy, it is exact, and
y x
u x (x3 + 3xy2) dx + k(y) x4/4 + 3x2y2/2 + k(y)
u
y
N 3x2y + k(y) 3x2y + y3, k(y) y4/4
Ex. 15
(sinx coshy) dx (cosx sinhy) dy 0, y(0) 0.
(sin x cosh y ) ( cos x sinh y )
Solution: Since sinxsinhy, it is exact, and
y x
u x sinx coshy dx + k(y) cosx coshy + k(y)
u
y
N cosx sinhy + k(y) cosx sinhy k(y) C u cosx coshy C
dx 1 x
[ Exercises ] 1. dy 0
x y
2 2 y y x y2
2
( x y )dx ( x y )dy
2. 0
x2 y2
dx 2 x
3. dy 0
xy y y3
O. D. E. 10
1 x x
[Answers] 1. y2 + 2Cx C 2 2. log x y tan y C 3. log y C
2 2
2. Integrating factors
If M(x, y) dx + N(x, y) dy 0 is not exact, we multiply by (x, y), M(x, y) dx + N(x, y)dy
0 is exact. The function (x, y) is called an integrating factor.
Ex. 16
ydx xdy 0 (2.11)
O. D. E. 11
M N
Solution : M y , N x 1, 1
y x
M N
y x 1 ( 1) 2
N x x
2
x dx 2
e e 2 ln x e ln x x 2
u Mdx f ( y ) x 2 ydx f ( y )
x x
1
x y f ( y)
u
N x 1 f ( y ) x 2 ( x ) f ( y ) 0 f ( y ) 0
y
y
The solution is x 1 y C1 C
x
Ex. 17
2xydx + (4y + 3x2) dy 0, y(0.2) 1.5.
M / y N / x 2x 6x 2
Solution: Since depends only on y, integrating factor is
M 2 xy y
2
y y dy y2 and u x 2xy3dx + k(y) x2y3 + k(y)
e
u
y
N 3x2y2 + k(y) y(4y + 3x2) k(y) y, and the general solution is
O. D. E. 12
xrys(mydx + nxdy) + xpyq(uydx + vxdy) 0. (2.13)
Multiplying (2.13) by xy and rearranging the terms, we have
(mx+ry+s+1 + ux+py+q+1) dx + (nx+r+1y+s+vx+p+1y+q) dy 0.
M N
This is exact, provided y x , i.e.,
Ex. 18
x4y(3ydx + 2xdy) + x2(4ydx + 3xdy) 0.
Solution: (3x4y2 + 4x2y)dx + (2x5y + 3x3)dy 0
Multiplying xy
(3x+4y+2 + 4x+2y+1)dx + (2x+5y+1 + 3x+3y)dy 0
M 3x+4y+2 + 4x+2y+1, N 2x+5y+1 + 3x+3y
M N
To be exact
y x
3( + 2)x+4y+1 + 4( + 1)x+2y 2( + 5)x+4y+1 + 3( + 3)x+2y
3( + 2) 2( + 5), 4( + 1) 3( + 3) 1, 2
1 6 4 4 3
u xMdx + f (y) x(3x5y4 + 4x3y3)dx + f (y) x y + x y + f (y)
2
u
N 2x6y3 + 3x4y2 + f (y) 2x6y3 + 3x4y2 f (y) 0 f (y) 0
y
1 6 4 4 3
The solution is x y +x y C
2
Ex. 19
2ydx + xdy + xy(3ydx + 2xdy) 0.
Solution: Multiplying xy and rearranging the terms, we get
(2xy+1 + 3x+1y+2)dx + (x+1y + 2x+2y+1)dy 0
For exactness 2( + 1)xy + 3( + 2)x+1y+1 ( + 1)xy + 2( + 2)x+1y+1
O. D. E. 13
2( + 1) + 1, 3( + 2) 2( + 2) 1, 0
u x (2xy + 3x2y2)dx + k(y) x2y + x3y2 + k(y)
u
y
N x2 + 2x3y + k(y) x2 + 2x3y k(y) C
Ex. 20
y y e2x.
Solution: F e dx ex
1 1
y [ Frdx C ] x [ e x e 2 x dx C ] e x [ e x dx C ]
F e
e ( e C ) e Ce x
x x 2x
O. D. E. 14
Ex. 21
y + 2y ex(3sin2x + 2cos2x).
Solution: F e 2dx e2x
1 1
y [ Frdx C ] 2 x [ e 2 x e x (3 sin 2 x 2 cos 2 x )dx C ]
F e
e 2 x [ e 3 x (3 sin 2 x 2 cos 2 x )dx C ] e 2 x ( e 3 x sin 2 x C )
e x sin 2 x Ce 2 x
Ex. 22
y + ytanx sin2x, y(0) 1.
Solution: F e tanxdx elnsecx secx
1 1
[ Frdx C ]
sec x
y [ sec x sin 2 xdx C ]
F
1
cos x[ 2 sin x cos xdx C ] cos x[ 2 sin xdx C ]
cos x
cos x ( 2 cos x C ) 2 cos 2 x C cos x
y(0) 1 1 2cos2(0) + Ccos(0) C 3
y 2cos2x + 3cosx
O. D. E. 15
This is linear in u.
Ex. 23
y Ay By2 (A, B are constants).
Solution: y2y Ay1 B
Let u y1 u y2y, there follows
B
u Au B F e Adx e Ax u e Ax ( e Ax Bdx C ) Ce Ax
A
1
y
B / A Ce Ax
Ex. 24
(ey + x)y 1.
Solution: The equation can be written in the form
dx
x ey which is linear in x.
dy
F e dy ey, x ey[ ey eydy + C] ey(y + C)
O. D. E. 16
Chapter 3 Ordinary Differential Equations of the First Order and Higher
Degree
Definition: f (x, y, y) can be replaced by a product of factors, rational in x and y, each of lower
degree in y, then f (x, y, y) is said to be reducible, otherwise, it is said to be irreducible.
Ex. 1
y2 + (x + y)y + xy 0.
Solution: (y + x) (y + y) 0
y + x 0 or y + y 0
x2
y + C or y Cex
2
Ex. 2
y2 + y2 1.
Solution : y 1 y 2
y 1 y 2 or y 1 y 2
sin1y x + C or cos1y x + C
y sin(x + C) or y cos(x + C)
The two solutions are the same.
O. D. E. 17
Differentiating (3.1), we obtain
f f dy
y x y dx
Since y is no longer present, it may be looked upon as an equation of the first order in the variables
x and y. It may happen that we can integrate this equation. Suppose its solution to be
Ex. 3
4xy2 + 2xy y 0.
dy dy
Solution : Differentiating ( 4 y 2 8 xy ) (2 y 2 x ) y 0
dx dx
dy dy dy dy
( 4 y 2 8 xy ) ( y 2 x ) 0 4 y ( y 2 x ) ( y 2x )0
dx dx dx dx
dy
( 4 y 1)( y 2 x )0
dx
dy dx dy C
1. y 2 x 0 2 0 ln x 2 ln y C1 xy 2 C y
dx x y x
Substituting into the original equation, we have
4C 2 Cx y 0 (y 4C)2 4Cx (y k)2 kx
2. 4y + 1 0 y 1/4, Substitute into the original equation, there follows
x x
y 0 x + 4y 0
4 2
This is the singular solution.
Ex. 4
O. D. E. 18
y + xy x4y2.
dy dy
Solution : Differenti ating y y x 4 x 3 y 2 2 x 4 y
dx dx
dy
(1 2 x 3 y )( 2 y x)0
dx
dy 2dx dy C
1. 2 y x 0 0 x 2 y C y 2
dx x y x
Substituting into the original equation, we have
y + x(C/x2) x4(C/x2)2 xy C2x C
2. 1 2x3y 0 y 1/2x3, The singular solution is
y + x(1/2x3) x4(1/2x3)2 y + 1/2x2 1/4x2 4x2y + 1 0
Ex. 5
O. D. E. 19
2yy2 2xy + y 0.
2 yy 2 y y
Solution : x x yy
2 y 2 y
dx dy 1 y
( y y ) (
dy dy 2 y 2 y
1 dy 1 y
( y y ) (
y dy 2 y 2 y 2
dy 1 dy
( y y ) ( y y )
dy 2 y 2 dy
dy dy dy
1. y y 0
dy y y
Substituting into the original equation, we have
2
C C
2 y 2 x y 0 2C 2 2Cx y 2 0
y y
1 1
2. 1 0 y , substituting into the original equation, we have
2 y 2
2
1 1 1
2 y 2 x y 0 y x0
2 2 2
Ex. 6
y2y2 3xy + y 0.
O. D. E. 20
y 2 y 2 y 1 2 yy 2 1 y 2 y dy
Solution : x 2
3y y 3y 3 3 y dy
2 yy 2 2 yy 2 1 dy y y 2 1 y dy
0 y 2 0
3y 3y 2
dy 3 y dy
y dy C
1. 2 0 y 2 y C y 2 , the general solution is
y dy y
2
C C
y 2 2 3x 2 y 0 C 2 3Cx y 3 0
y y
yy 2 1 1
2. 0 y , the singular solution is
3 y
1 1 1
y 3x y 0 2 y 3x 0
2
y y y
O. D. E. 21
y Cx f (C)
Ex. 7
xy2 yy + 3 0.
Solution: The equation can be written as
3
y xy
y
This is Clairaut equation, the general solution is
C 2x Cy + 3 0
Ex. 8
x2y2 2(xy 4)y + y2 0.
Solution: The equation can be written as
x2y2 2xyy + 8y + y2 0
(xy y)2 + 8y 0
This is Clairaut equation, the general solution is
C 2x2 2C(xy 4) + y2 0
O. D. E. 22
Chapter 4 Ordinary Differential Equations of Higher Order
1. Linear dependence
Definition: A linear combination of n functions y1(x), y2(x),…, yn(x), is an expression of the form
k1y1(x) + k2y2(x) +…+ knyn(x),
where ki’s are constants.
Definition: If (4.2) holds on interval I only for k1 k2 … kn 0, then y1(x), y2(x),…, yn(x) are
said to be linearly independent.
Assume that each of a set of n functions y1(x), y2(x),…, yn(x) possesses n finite derivatives at
all points of an interval I. Then
dy dy dy
k1 1 k 2 2 k n n 0
dx dx dx
2
d y1 2
d y2 d 2 yn
k1 k 2 k n 0
dx 2 dx 2 dx 2
(4.4)
d n 1 y1 d n 1 y 2 d n 1 y n
k1 k 2 k n 0
dx n 1 dx n 1 dx n 1
If the functions y1(x), y2(x),…, yn(x) are linearly dependent, i.e., there exists nontrivial solution of
ki’s, the determinant
O. D. E. 23
y1 y2 yn
y1 y 2 y n
W(y1, y2 ,…, yn; x) (4.5)
y1( n 1) y 2( n 1) y n( n 1)
vanishes identically over I. This determinant is called the Wronskian of the functions y1(x), y2(x),
…, yn(x). If the Wronskian of y1(x), y2(x),…, yn(x) is not identically zero over I., then the
functions y1(x), y2(x),…, yn(x) are linearly independent.
[Note] The vanishing of the Wronskian is necessary but not sufficient condition for linear
dependence of a set of functions.
Ex. 1
Show that ex, e2x are linearly independent on all interval.
ex e2 x
Solution: W(e , e ; x) x
x 2x
exe2x 0
e 2e 2 x
Ex. 2
Show that x, x2 are linearly independent on {x xR, x 0}.
x x2
Solution: W(x, x2; x) 1 2x
x2 0 if x 0.
Principle of superposition: If y1(x), y2(x) are the solutions of the homogeneous linear differential
equation
pn(x)y (x) + pn1(x)y(n1)(x) +…+ p1(x)y(x) + p0(x)y(x) 0
(n)
(4.6)
then k1y1(x) + k2y2(x) is also a solution of (4.6), where k1, k2 are arbitrary constants.
An nth order linear differential equation (4.6) has n linearly independent solutions y1(x),
y2(x),…, yn(x), and its general solution can be expressed as
yh(x) k1y1(x) + k2y2(x) +…+ knyn(x)
where ki’s are arbitrary constants, and yh(x) is called homogeneous solution or complementary
solution of (4.1). Now suppose that one particular solution of (4.1), say yp(x), then the complete
solution of (4.1) is
y(x) yh(x) + yp(x)
O. D. E. 24
Form: any(n)(x) + an1y(n1)(x) +…+ a1y(x) + a0y(x) r(x) (4.7)
or Ly r(x)
where a0, a1, …, an are constants, and
dn d n 1 d n 2 d
L an n a n1 n 1 a n 2 n 2 a1 a0 .
dx dx dx dx
Let yh(x) ex, substitute into the homogeneous equation of (4.7), i.e., r(x) 0, we have
(ann + an1n1 +…+ a1 + a0)ex 0
Since ex 0, we get
ann + an1n1 +…+ a1 + a0 0 (4.8)
Equation (4.8) is called the characteristic equation of the homogeneous equation of (4.7).
[Case 1] If (4.8) has n distinct real roots, say 1, 2,…, n, then
yh(x) k1e 1x k 2 e 2 x k n e n x
and
(Lex) 2( 1)( 3)( 4)…( n)ex
e x
but
L e x 1
L
L( xe 1x ) 0
1
Ex. 3
y 2y y + 2y 0.
O. D. E. 25
Solution: The characteristic equation is
3 22 + 2 0 ( + 1) ( 1) ( 2) 0 1, 1, 2
The general solution is y k1ex + k2ex + k3e2x
Ex. 4
y(5) 3y(4) + 3y(3) y 0.
Solution: The characteristic equation is
5 34 + 33 2 0 2 ( 1)3 0 0, 0, 1, 1, 1
The general solution is y k1 + k2x + (k3 + k4x + k5x2)ex
[Case 3] If (4.8) has simple complex root, say 1 i, then solution
k1e(+i)x + k2e(i)x k1ex(cosx + isinx) + k2ex(cosx isinx)
ex[(k1 + k2)cosx + i(k1 k2)sinx] ex(C1cosx + C2sinx)
where C1 k1 + k2, C2 i(k1 k2)
Hence the solution can be expressed as ex(k1cosx + k2sinx)
Ex. 5
y 2y + 2y 0, y(0) 0.5, y(0) 1, y(0) 2.
Solution: The characteristic equation is
3 22 + 2 0 (2 2 + 2) 0 0, 1 i
The general solution is y(x) k1 + ex(k2cosx + k3sinx)
y(0) 0.5 0.5 k1 + k2
y(0) 1 1 k2 + k3 k1 2.5, k2 2, k3 1
y(0) 2 2 2k3
∴y 2.5 + ex(2cosx + sinx)
Ex. 6
y(7) + 18y(5) + 81y(3) 0.
O. D. E. 26
Solution: The characteristic equation is
7 + 185 + 813 0 3 (2 + 9)2 0 0, 0, 0, 3i, 3i
The general solution is y(x) k1 + k2x + k3x2 + (k4 + k5x)cos3x + (k6 + k7x)sin3x
Ex. 7
y(4) y 4.5e2x.
Solution: The characteristic equation is
4 1 0 (2 1) (2 + 1) 0 1, i
yh(x) c1ex + c2ex + c3cosx + c4sinx
Let yp(x) Ae2x yp 2Ae2x yp 4Ae2x yp 8Ae2x yp(4) 16Ae2x
Substitute in the given equation, there follows
(16A A)e2x 4.5e2x A 0.3
The complete solution is y(x) yh(x) + yp(x) c1ex + c2ex + c3cosx + c4sinx + 0.3e2x
Ex. 8
y 3y + 3y y 30ex.
O. D. E. 27
Solution: The characteristic equation is
3 32 + 3 1 0 ( 1)3 0 1, 1, 1
yh(x) (k1 + k2x + k3x2)ex
Let yp(x) Ax3ex yp A(x3 + 3x2)ex, yp A(x3 + 6x2 + 6x)ex,
yp A(x3 + 9x2 + 18x + 6)ex
Substituting into the problem, we have
A[(x3 + 9x2 + 18x + 6) 3(x3 + 6x2 + 6x) + 3(x3 + 3x2) x3]ex 30ex A 5
The complete solution is y(x) yh(x) + yp(x) (k1 + k2x + k3x2 + 5x3)ex
Ex. 9
y 2y + y x2ex.
Solution: The characteristic equation is
2 2 + 1 0 ( 1)2 0 1, 1
yh(x) (k1 + k2x)ex
Let yp(x) (Ax4 + Bx3 + Cx2)ex yp [Ax4 + (4A + B)x3 + (3B + C)x2 + 2Cx]ex,
yp [Ax4 + (8A + B)x3 + (12A + 6B + C)x2 + (6B + 4C)x + 2C]ex,
Substituting into the problem, we have
[(A 2A + A)x4 + (8A + B 8A 2B + B)x3 + (12A + 6B + C 6B 2C + C)x2
+ (6B + 4C 4C)x + 2C]ex x2ex A 1/12, B 0, C 0
The complete solution is y(x) yh(x) + yp(x) (k1 + k2x + 1/12x4)ex
Ex. 10
y + 4y x2 + cosx.
Solution: The characteristic equation is
2 + 4 0 2i yh(x) k1cos2x + k2sin2x
Let yp(x) (Ax2 + Bx + C) + (Dcosx + Esinx)
yp 2Ax + B Dsinx + Ecosx, yp 2A Dcosx Esinx
Substituting into the problem, we have
4Ax2 + 4Bx + (2A + 4C) + (D + 4D)cosx + (E + 4E)sinx x2 + cosx
A 1/4, B 0, C 1/8, D 1/3, E 0
1 1 1
The complete solution is y(x) yh(x) + yp(x) k1cos2x + k2sin2x + x cosx
2
4 8 3
O. D. E. 28
The homogeneous solution of the equation (4.7) is
n
yh(x) k1y1(x) + k2y2(x) +…+ knyn(x) k i y i ( x )
i 1
u [a
i 1
i
(n)
n yi a n 1 yi( n 1) a1 yi a0 yi ] a n ui y i( n 1) r
i 1
n
u y
i 1
i
( n 1)
i r / an (4.12)
The equations (4.9), (4.10),…,(4.12) gives a system of n equations for the unknown functions
u1, u2,…, un, where we can find u1, u2,…, un.
Ex. 11
O. D. E. 29
y + y tanx.
Solution:2 + 1 0 i yh(x) k1cosx + k2sinx
Let yp(x) u1cosx + u2sinx, then
yp (u1sinx + u2cosx) + (u1cosx + u2sinx)
Set u1cosx + u2sinx 0…………
yp (u1cosx u2sinx) + (u1sinx + u2cosx)
Substitute in the given equation, there follows
[(u1cosx u2sinx) + (u1sinx + u2cosx)] + (u1cosx + u2sinx) tanx
u1sinx + u2cosx tanx…………
From and we get
0 sin x
u1 tan x cos x sinx tanx
W(cos x , sin x )
u1 sinx tanx dx tanx dcosx cosx tanx cosx sec2x dx sinx lnsecx + tanx
cos x 0
u2 sin x tan x sinx u2 sinx dx cosx
W(cos x, sin x )
y(x) yh(x) + yp(x) k1cosx + k2sinx + (sinx ln secx + tanx )cosx + ( cosx) sinx
(k1 ln secx + tanx )cosx + k2sinx
Ex. 12
y + y cotx.
Solution: 3 + 0 0, i yh(x) k1 + k2cosx + k3sinx
Let yp(x) u1 + u2cosx + u3sinx, then
u1 + u2cosx + u3sinx 0
u2sinx + u3cosx 0
u2cosx u3sinx cotx
u1 cotx, u2 cotxcosx, u3 cosx
u1 lnsinx, u2 cosx lncscx cotx, u3 sinx
y(x) yh(x) + yp(x) k1 + k2cosx + k3sinx 1 cosx lncscx cotx + lnsinx
O. D. E. 30
Form: anxny(n)(x) + an1xn1y(n1)(x) +…+ a1xy(x) + a0y(x) r(x) (4.13)
d dt d 1 d
Let x et t
dx dx dt e dt
d d
x
dx dt
d2 t 1 d 1 d
2
d d 2 1 d 1 d
x2 2 x2 x t t e t 2 t
dx dx dx e dt e dt e dt e dt
2
d d d d
2 1
dt dt dt dt
…
d3 3t 1 d 1 d 1 d 2t d 1 1 d 1 d 2t d 1 d 1 d
2 2
x3 3
e t t t e
t t 2 t
e
2t 2
2 t
dx e dt e dt e dt dt e e dt e dt dt e dt e dt
1 d3 1 d2 1 d2 1 d d3 d2 d
e 2 t 2 t 3 2 2t 2 2 t 2 2 2 t 3 3 2 2
e dt e dt e dt e dt dt dt dt
d d d
1 2
dt dt dt
………………………………………………………
dn d d d d d
xn n
1 2 3 n 1
dx dt dt dt dt dt
The transformed equation thus becomes linear with constant coefficients, and y then can be
determined in terms of t.
Ex. 13
x2y 3xy + 4y x2lnx.
Solution: Let x et t lnx, we get
d d dy d2y dy
1 y 3 4 y te 2t
2
4 4 y te 2 t
dt dt dt dt dt
The characteristic equation is 2 4 + 4 0 2, 2 yh (k1 + k2t)e2t
Assume particular solution is yp (At3 + Bt2) e2t
yp [2At3 + (3A + 2B)t2 + 2Bt]e2t
yp [4At3 + (12A + 4B)t2 + (6A + 8B)t + 2B]e2t
[(4A 8A + 4A)t3 + (12A + 4B 12A 8B + 4B)t2 + (6A + 8B 8B)t + 2B]e2t te2t
A 1/6, B 0
1 1
y yh + yp (k1 + k2t + t3)e2t [k1 + k2lnx + (lnx)3]x2
6 6
Ex. 14
x2y 4xy + 6y 7x4sinx.
Solution: Let x et t lnx, we get
O. D. E. 31
d2y dy
2
5 6 y 7e 4 t sin e t 2 5 + 6 0 2, 3 yh k1e2t + k2e3t
dt dt
Let yp u1e2t + u2e3t
u1e2t + u2e3t 0
2u1e2t + 3u2e3t 7e4tsinet
0 e 3t
u1 7e 4 t sin e t 3e 3t 7e2tsinet u1 7etcoset + 7sinet
W(e 2 t , e 3t )
e 2t 0
u2 2e 2t
7e sin e t 7etsinet u2 7coset
4t
W(e 2 t , e 3t )
y yh + yp k1e2t + k2e3t + (7etcoset + 7sinet)e2t + 7e3tcoset
k1e2t + k2e3t + 7e2tsinet k1x2 + k2x3 + 7x2sinx
Ex. 15
O. D. E. 32
dx dy
dt 2 x dt 6 y 2e
t
dx dy
2 3x 3 8 y 1
dt dt
d ( D 2) x ( D 6) y 2e t (i)
Solution: Let D
dt ( 2 D 3) x (3D 8) y 1 (ii)
(i) 2 (ii) x + (D + 4)y 4et + 1………(iii)
(i) (D + 2)(iii) [(D + 6) (D + 2)(D + 4)]y 2et (D + 2)(4et + 1)
(D2 D 2)y 2et (4et + 8et + 2) (D2 D 2)y 10et 2………(iv)
The characteristic equation of (iv) is 2 2 0 1, 2 yh k1et + k2e2t
Assume yp Aet + B yp Aet yp Aet
(iv) (A A 2A)et 2B 10et 2 A 5, B 1
y yh + yp k1et + k2e2t + 5et + 1
Substituting into (iii), we get
x (k1 4k1)et + (2k2 4k2)e2t + (5 20)et 4 + 4et + 1 5k1et 2k2e2t 11et 3
x 5 t 2 2 t 11 t 3
k1 e k 2 e e
y 1 1 5 1
Ex. 16
O. D. E. 33
dx dy
dt x dt y sec t
dx dy
4 5 x 5 4 y 5 sec t
dt dt
d ( D 1) x ( D 1) y sec t (i)
Solution: Let D
dt ( 4 D 5) x (5D 4) y 5 sec t (ii)
(ii) (i) 4 x + Dy sect………(iii)
(D + 1)(iii) (i) (D2 + 1)y sect tant………(iv)
The characteristic equation of (iv) is 2 + 1 0 i yh k1cost + k2sint
Assume yp u1cost + u2sint
u1cost + u2sint 0
u1sint + u2cost sect tant
0 sin t
u1 tan 2 t u1 tan t t
sec t tan t cos t
cos t 0
u 2 tan t u 2 ln | sec t |
sin t sec t tan t
O. D. E. 34
pn(x)y(n)(x) + pn1(x)y(n1)(x) +…+ p1(x)y(x) + p0(x)y(x) r(x) (4.14)
(1) Exact differential equations
Consider a simple case of equation (4.14):
p3y + p2y + p1y + p0y (p3y) p3y + p2y + p1y + p0y (p3y) + (p2 p3)y + p1y + p0y
(p3y) + [(p2 p3)y] (p2 p3)y + p1y + p0y
(p3y) + [(p2 p3)y] + (p1 p2 + p3)y + p0y
(p3y) + [(p2 p3)y] + [(p1 p2 + p3)y] (p1 p2 + p3)y + p0y
(p3y) + [(p2 p3)y] + [(p1 p2 + p3)y] + (p0 p1 + p2 p3)y
If the coefficient of y is identically zero, the other three terms can be integrated, and we have a
new equation of less order. In a similar way, equation (4.1) is said to be exact, if
n
( 1)
i 0
i
pi( i ) ( x ) 0
Ex. 17
(x 1)y + (x + 1)y + y 2x.
Solution: Since p0 p1 + p2 0, the equation is exact.
[(x 1)y] y + (x + 1)y + y 2x
[(x 1)y] + xy + y 2x
[(x 1)y] + (xy) 2x (x 1)y + xy x2 + C1
x x 2 C1
y y
x 1 x 1
This is a linear first order equation.
x
h x 1 dx x ln( x 1)
x 2 C1 h
y e h
e dx C 2
ex
e ( x
x 2
C1 )dx C 2
x 1 x 1
x
e
[ x 2 2 x 2 C1 )e x C 2 ]
x 1
One of the important properties of linear differential equations is the fact that, if one
O. D. E. 35
homogeneous solution of an equation of order n is known, a new linear differential equation of
order n 1, determining the remainder of the solution, can be obtained.
Suppose u(x) is the homogeneous solution, substituting y v(x) u(x) into the original
equation, we have a new equation of less order. We can try emx, cosmx, sinmx, xm for the
homogeneous solution.
Ex. 18
xy + (x 1)y y 0.
Solution: Since ex is a homogeneous solution, let y vex
y vex vex, y vex 2vex + vex
x(v 2v + v)ex + (x 1)(v v)ex vex 0
xv (x + 1)v 0
This is a first order equation in v
x 1
v C e x
dx
C1xex
1
O. D. E. 36
Let p y, the resulting equation is of order n 1 in p.
Ex. 19
(1 + x2)y + y2 + 1 0.
dp dp dx
Solution: Let p y (1 x ) p2 1 0 2 0
2
dx p 1 1 x2
px C x
tan 1 p tan 1 xC C1 p 1
1 px C1 x 1
C1 x 1 C1 1 / C1
y dx C2 dx C 2
C1 x 1 C1 C1 x 1
x 1
(C1 1 / C1 ) ln(C1 x 1) C 2
C1 C1
3 2
d2y dy d 2 y dy
[Exercises] 1. x 2. 0
dx 2 dx dx 2 dx
dp dp dy dp
Let p y, y dx dy dx p dy ,…, the equation can be reduce to an (n 1) order, where p is
Ex. 20
yy + y2 y 0.
dp
Solution: Let p y, y p dy , the equation reduces to
O. D. E. 37
dp dp dy
yp p2 p 0 0 ln( p 1) ln y C
dy p 1 y
dy ydy C1
( p 1) y C1 y y C1 dx 1 dy dx
dx C1 y C1 y
y C1ln(y + C1) x + C2
2 2
d 2 y dy d 2 y dy
[Exercises] 1. y 2. 0
dx 2 dx dx 2 dx
Ex. 21
yy y2 6xy2 0.
Solution: Let y ez y ezz, y ez(z + z2)
e2z(z + z2) z2e2z 6xe2z 0 z 6x 0 z x3 + C1x + C2
3
y e x C1x C2
2 2
d 2 y dy d 2 y dy
[Exercises ] 1. y 2 . y 0
dx 2 dx dx 2 dx
[Answers] 1. y C 2 e C1x 2. y 2 C 2 x C1
4. Isobaric equations
O. D. E. 38
the original equation.
Ex. 22
x2yy 3xyy + x2y2 + 2y2 + x2 0.
Solution: weight: 2 + m + (m 2), 1 + m + (m 1), 2 + 2(m 1), 2m, 2 m 1
dy dy dt du dx du
Let y ux, x e t y x u e t u
dx dt dx dt dt dt
dx d u du t
2
d du
y u 2 e
dt dt dt dt dt
2
d 2 u du du du
e ue 2 e t 3e t ue t u e 2 t u 2u 2 e 2 t e 2 t 0
2t t
dt dt dt dt
2
d 2 u du d du du
u 2
1 0 u 1 0 u t C1
dt dt dt dt dt
u 2
(C t ) 2
udu (C1 t )dt 1 C2
2 2
2
1 y (C1 ln x ) 2
C2
2 x 2
[Exercises] 1. Show that the substitution y u/[Q(x)u], where u du/dx, reduces the nonlinear
first-order equation
y + P(x)y + Q(x)y2 R(x)
to the second-order equation
d 2u Q du
2
P RQu 0
dx Q dx
and also that the substitution y y1 + 1/v, where y1(x) is any known solution of the
given equation, leads to the linear first-order equation
v (P + 2Qy1)v Q.
(The nonlinear form is known as Riccati’s equation.)
2. Use the procedures suggested above to obtain the general solution of the equation
x2y + xy + x2y2 1 in the form (x2 k)/(x2 + k), where k is an arbitrary constant.
O. D. E. 39
Chapter 5 Series Solutions of Differential Equations
A
m 0
m ( x x0 ) m (5.1)
is called a power series, where A0, A1, A2,…, are constants, called the coefficients of the series. x0 is
a constant, called the center of the series, x is a variable.
Ex1. y + y 0
Solution: Assume y Am x
m
, then
m 0
m(m 1) Am x m 2 Am x m 0
m0 m 0
m(m 1) A
m2
m x m 2 Am 2 x m 2 0
m2
[m(m 1) A
m2
m Am 2 ) x m 2 0
m(m 1) Am Am 2 0 m 2,3,4,
A A A A A
A2 0 , A4 2 0 , A6 4 0 ,
2! 4 3 4! 65 6!
A1 A A A A
A3 , A5 3 1 , A7 5 1 ,
3! 5 4 5! 76 7!
x 2
x 4
x 6
x 3
x5 x7
y A0 1 A1 x A0 cos x A1 sin x
2! 4! 6! 3! 5! 7!
Ex2. y + xy + y 0
O. D. E. 40
Solution: Assume y Am x
m
, then
m 0
m(m 1) A
m0
m x m2 mA
m0
m xm A
m0
m xm 0
m(m 1) Am x m 2
m0
(m 2) Am 2 x m 2
m2
A
m2
m2 x m2 0
[m(m 1) A
m2
m ( m 2 1) Am 2 ] x m 2 0
Am 2
Am m 2,3,4,
m
A A A A A0
A2 0 , A4 2 0 , A6 4 ,,
2 4 24 6 246
A0 A
A2 k ( 1) k ( 1) k k 0
2 4 6 ( 2k ) 2 k!
A A A A A1
A3 1 , A5 3 1 , A7 5 ,,
3 5 35 7 357
A1 2 4 6 (2k ) A1 2 k k! A1
A2 k 1 (1) k (1) k (1) k
3 5 7 (2k 1) (2k 1)! (2k 1)!
(1) k
(1) 2 k! 2 k 1
k k
y A0 1 k x 2 k A1 x x
k 1 2 k ! k 1 ( 2 k 1)!
(1) k 2 k
(1) k 2 k k! 2 k 1
A0 k x A1 x
k 0 2 k! k 0 ( 2k 1)!
[Exercises] Obtain the general solutions of the following differential equations in terms of power
series around x 0.
1. y xy 0 2. y + xy 1 x + x2
x 3m
x 3m 1
[Answers]1. y A0 1 A1 x
m 1 2 3 5 6 (3m 1)(3m) m 1 3 4 6 7 (3m)(3m 1)
x2
2. y x 1 Ae 2
O. D. E. 41
A homogeneous second-order linear differential equation
d2y dy
2
a1 ( x ) a0 ( x) y 0 (5.2)
dx dx
1. If a1(x), a0(x) are both analytic at point x x0, the point x x0 is said to be an ordinary point of
the differential equation, then the equation possesses two linearly independent solutions which
are both expressible in the form A
m 0
m ( x x0 ) m .
2. Although a1(x) or a0(x) are not analytic at point x x0, (x x0) a1(x) and (x x0)2 a0(x) are both
analytic at point x x0, the point x x0 is said to be a regular singular point, the equation has at
m 0
3. If a1(x) and a0(x) are not the case mentioned above, the equation is said to be irregular singular
at point x x0, and a nontrivial solution of this type may or may not exist.
O. D. E. 42
(R0 + R1x + R2x2 +…) [r(r 1)A0xr + (r + 1)rA1xr + 1 + (r + 2)(r + 1)A2xr + 2 +…]
+ (P0 + P1x + P2x2 +…) [rA0xr + (r + 1)A1xr + 1 + (r + 2)A2xr + 2 +…]
+ (Q0 + Q1x + Q2x2 +…) [A0xr + A1xr + 1 + A2xr + 2 +…] +… 0
after multiplying and collecting the coefficients of successive power of x,
[R0r(r 1) + P0r + Q0]A0xr
+{[R0(r + 1)r + P0(r + 1) + Q0]A1 + [R1r(r 1) + P1r + Q1]A0}xr + 1
+{[R0(r + 2)(r + 1) + P0(r + 2) + Q0]A2 + [R1(r + 1)r + P1(r + 1) + Q1]A1
+ [R2r(r 1) + P2r + Q2]A0}xr + 2 +… 0 (5.5)
In order to abbreviate this relation, we next define the functions
fi(r + m) Ri(r + m)(r + m 1) + Pi(r + m) + Qi
(5.5) becomes
r r+1 r+2
f0(r)A0x + [f0(r + 1)A1 + f1(r)A0]x + [f0(r + 2)A2 + f1(r + 1)A1 + f2(r)A0]x +… 0
m
f0(r)A0x + [ f 0 ( r m ) Am f i ( r m i ) Am i ] x
r m
r 0
m 1 i 1
The coefficients of all powers of x must vanish independently. The vanishing of the lowest power xr
gives the requirement
f0(r) 0 (5.6)
This equation determines two values of r and is called the indicial equation. The vanishing of the
coefficient of xr + m gives the recurrence formula
m
f 0 ( r m ) Am f i ( r m i ) Ami (5.7)
i 1
which determines each Am in terms of the preceding A’s, and hence in terms of A0, if for each
positive m the quantity f0(r + m) is not zero.
The root of the indicial equation has the following three cases:
1. Distinct roots not differing by an integer (r r1, r2), two independent solutions are
y1 ( x ) x r1 Am x m
m 0
y2 ( x) x r2
A x
m 0
m
m
2. Double roots r r1 r2, only one solution can be obtained. From recurrence formula (5.7), (5.4)
can be written as
y ( x, r ) x r Am ( r ) x m (5.8)
m 0
and
Ly(x,r) f(r)A0xr (r r1)2A0xr (5.9)
The right-hand member of (5.9) vanishes when r r1, and the first solution is
y1(x) y(x,r1) x
r1
A
m 0
m ( r1 ) x m
O. D. E. 43
Since r r1 is a repeated zero of the right-hand member of (5.9), it follows
r L y ( x, r ) 0 L y ( x, r ) 0
r r1 r r r1
Hence a second solution is
y2 ( x) y ( x, r )
r r r1
we have
L[(r r2)y(x,r)] f(r)A0xr (r r1)(r r2)2A0xr (5.11)
and
[( r r2 ) L y ( x, r )] 0 L [( r r2 ) y ( x, r )] 0
r r r2 r r r2
Hence a second solution is
y2 ( x) [( r r2 ) y ( x, r )]
r r r2
m 0
4 ( m r )( m r 1) Am x m r 1 2 ( m r ) Am x m r 1 Am x m r 0
m 0 m 0 m 0
4 ( m r )( m r 1) Am x m r 1 2 ( m r ) Am x m r 1 Am 1 x m r 1 0
m 0 m 0 m 1
[4 r ( r 1) 2 r ] A0 x r 1 {4( m r )( m r 1) 2( m r )] Am Am 1 }x m r 1 0
m 0
O. D. E. 44
The indicial equation is
4r(r 1) + 2r 0 r 1/2, 0
The recurrence formula is
Am 1
Am , m 1,2,3,…
2( m r )( 2m 2 r 1)
r1 1/2
Am 1
Am , m 1, 2, 3,
( 2m )( 2m 1)
A A A A A0
A1 0 , A2 1 0 , A3 0 ,, Am ( 1) m
3! 5 4 5! 7! ( 2m 1)!
xm
y1 ( x ) x ( 1) m
m 0 ( 2m 1)!
r2 0
Am 1
Am , m 1, 2, 3,
( 2m 1)(2m )
A A A A A0
A1 0 , A2 1 0 , A3 0 , , Am ( 1) m
2! 43 4! 6! ( 2m )!
xm
y2 ( x) ( 1) m
m 0 ( 2m )!
xm xm
The solution is y C1 y1 C 2 y 2 C1 x ( 1) m C 2 ( 1) m
m 0 ( 2m 1)! m 0 ( 2m )!
m 0
O. D. E. 45
( m r )( m r 1) A
m 0
m x m r ( m r )( m r 1) Am x m r 1
m 0
3 ( m r ) Am x m r ( m r ) Am x m r 1 Am x m r 0
m0 m 0 m 0
( m r 1)( m r 2) A
m 1
m 1 x m r 1 ( m r )( m r 1) Am x m r 1
m 0
3 ( m r 1) Am 1 x m r 1 ( m r ) Am x m r 1 Am 1 x m r 1 0
m 1 m 0 m 1
[ r ( r 1) r ] A0 x r 1
[( m r ) Am 1 ( m r ) Am ] x m r 1 0
2 2
m 1
m 0 1 x
Let y2 uy1 y2 uy1 + uy1 y2 uy1 + 2uy1 + uy1
Substitute into the given equation
x(x 1)(uy1 + 2uy1 + uy1) + (3x 1)(uy1 + uy1) + uy1 0
[x(x 1)y1 + (3x 1)y1 + y1]u + x(x 1)y1u + [2x(x 1)y1 + (3x 1)y1]u 0
du dx
xu u 0 0 u x 1 u ln x
u x
ln x
y2
1 x
C C ln x
The Solution is y C1 y1 C 2 y 2 1 2
1 x 1 x
Solution: Assume y x Am x
r m
m 0
O. D. E. 46
(m r )(m r 1) A
m 0
m x m r ( m r )( m r 1) Am x m r 1
m 0
( m r ) Am x m r Am x m r 0
m 0 m 0
(m r 1)(m r 2) A
m 1
m 1 x m r 1 ( m r )( m r 1) Am x m r 1
m 0
( m r 1) Am 1 x m r 1 Am 1 x m r 1 0
m 1 m 1
r ( r 1) A0 x r 1 [( m r 2) 2 Am 1 ( m r )( m r 1) Am ] x m r 1 0
m 1
O. D. E. 47