ODE Assign 5 Soln
ODE Assign 5 Soln
ODE Assign 5 Soln
(For calculations of Particular Integrals by operator method, see Simmons books, page 161,
section 23 of the chapter Second order linear equations.)
d2 u du
2
+ −12u = 0 =⇒ m2 +m−12 = 0 =⇒ m = −4, 3 =⇒ u(t) = Ae−4t +Be3t = y(et ).
dt dt
The general solution is thus
A
y(x) = + Bx3 .
x4
(ii) Using the substitution x = et , the given equation reduces to,
d2 u du
2
+ 4 + 13u = 0 =⇒ m2 + 4m + 13 = 0 =⇒ m = −2 ± 3i.
dt dt
Thus
u(t) = e−2t (A cos 3t + B sin 3t) = y(et ).
The general solution is
1
y(x) = 2
A cos(3 ln x) + B sin(3 ln x) .
x
d2 u du
2
− 2 + u = 0 =⇒ m2 − 2m + 1 = 0 =⇒ m = 1, 1 =⇒ u(t) = et (A + Bt) = y(et )
dt dt
The general solution is thus
y(x) = ex (A + B ln x).
d
2. (Higher order Cauchy-Euler equations) Let us denote D = dx and D = dtd where x = et .
Show that
xD = D, x2 D2 = D(D − 1), x3 D3 = D(D − 1)(D − 2).
Hence conclude that (x3 D + ax2 D2 + bxD + c)y = 0, x > 0 is transformed into constant
coefficients ODE [D(D − 1)(D − 2) + aD(D − 1) + bD + c]y = 0 by the substitution
x = et .
Solution:
Given x = et , so dx
dt
= et = x. Now, by chain rule, dxd
= dtd dx
dt
= e−t dtd . Thus xD =
D. Differentiating this with respect to x, we have xD2 + D = D2 dxdt
= D2 e−t , =⇒
x2 D2 + xD = D2 , =⇒ x2 D2 = D2 − D = D(D − 1).
Differentiating x2 D2 = D2 − D with respect to x, we have x2 D3 + 2xD2 = [D3 − D2 ]e−t ,
=⇒ x3 D3 = D3 − D2 − 2D(D − 1) = D(D − 1)(D − 2).
3. Find a particular solution of each of the following equations by operator methods and
hence find its general solution:
(i) y 00 + 4y = 2 cos2 x + 10ex (ii)(T) y 00 + y = sin x + (1 + x2 )ex
(T) (iii) y 00 − y = e−x (sin x + cos x) (iv) y 000 − 3y 00 − y 0 + 3y = x2 ex
Solution:
(i) Characteristic equation m2 + 4 = 0 =⇒ m = ±2i. Hence homogeneous solution
yh = A cos 2x + B sin 2x. Now r(x) = 2 cos2 x + 10ex = cos 2x + 1 + 10ex . Let D ≡ d/dx
and yp be the particular solution. Then
1 1
2
1 = 2
e0x = 1/4.
D +4 D +4
1 1
10ex = 10 2 ex = 2ex .
D2+4 1 +4
1 1 2ix 1 2ix
2
e2ix = x e =x e = xe2ix /4i.
D +4 2D 2.2i
Taking real part
1
cos 2x = x sin 2x/4.
D2 +4
Adding, we get the particular solution as
x sin 2x 1
yp = + + 2ex .
4 4
ex 1 2 ex
= (x + 1) = (1 − D − D2 /2 + (D + D2 /2)2 + · · · )(x2 + 1)
2 1 + D + D2 /2 2
ex ex ex
= (1−D−D2 /2+D2 /2+· · · )(x2 +1) = (1−D+D2 /2+· · · )(x2 +1) = (1+x2 −2x+1)
2 2 2
Thus the general solution is
x2 x
x cos x
y = A cos x + B sin x − + 1−x+ e
2 2
Then the particular solution is obtained by adding the real and imaginary parts:
e−x (cos x − 3 sin x)
yp (x) =
5
Thus the general solution is
e−x (cos x − 3 sin x)
y = Aex + Be−x +
5
(iv) Characteristic equation m3 − 3m2 − m + 3 = 0 =⇒ m = −1, 1, 3. Hence homoge-
neous solution yh = Ae−x + Bex + Ce3x . Now r(x) = x2 ex . Let D ≡ d/dx and yp be the
particular solution. Then
1 1 1
x2 e x = x2 ex = ex 3 x2
D3 2
− 3D − D + 3 3
(D − 1) − 4(D − 1) D − 4D
1 x 1 x 1 1 e x x3 x
= ex x 2
= e (1 + D 2
/4 + · · · )x 2
= e (x 2
+ ) = − ( + ).
−4D(1 − D2 /4) −4D −4D 2 4 3 2
So the particular integral is
x x3
x
yp (x) = −e + .
8 12
Thus the general solution is
x x3
−x x 3x x
y = Ae + Be + Ce −e + .
8 12
4. Solve y 00 + y 0 − 2y = ex .
Solution: Characteristic equation of the homogeneous part is: m2 + m − 2 = 0, m =
1, −2. Solution for the homogeneous part: c1 ex + c2 e−2x .
Particular integral:
1 x 1 x xex xex
e = x e = = .
D2 + D − 2 2D + 1 (2.1 + 1) 3
General solution:
x −2x xex
c1 e + c2 e + .
3
5. Solve by using operator method (D2 + 9)y = sin 2x cos x.
Solution:
Characteristic equation of the homogeneous part is: m2 + 9 = 0, m = ±3i. Solution for
the homogeneous part: c1 cos 3x + c2 sin 3x.
Particular integral:
1 1
sin 2x cos x = (sin 3x + sin x).
D2 +9 2
2(D + 9)
Now
1
eix = eix /(i2 + 9) = (cos x + i sin x)/8.
+9D2
Taking the imaginary part,
1
(sin x) = sin x/16.
2(D2 + 9)
Now
1 3ix xe3ix
e = = (x cos 3x + ix sin 3x)/6i.
D2 + 9 2.3i
Taking the imaginary part,
1
(sin 3x) = −(x cos 3x)/12.
2(D2 + 9)
General solution:
1 1
P.I = 1 + x + x2 = 1 + x + x2
D2 − 6D + 9 2
9(1 + (D − 6D)/9)
1
= [1 − (D2 − 6D)/9 + (D2 − 6D)2 /81 − · · · ](1 + x + x2 ).
9
1
= [1 + 2D/3 + D2 /3 + · · · ](1 + x + x2 ) = (1 + x + x2 + 2/3 + 4x/3 + 2/3)
9
1
= (7/3 + 7x/3 + x2 ).
9
7. Find P.I: y 00 + 9y = x cos x.
Solution:
Consider
1 1 1 1
xeix = eix x = eix 2 x = eix x
D2 +9 2
(D + i) + 9 D + 2iD + 8 2
8(1 + D /8 + iD/4)
1 1
= eix (1 − D2 /8 − iD/4)x = eix (x − i/4).
8 8
Taking the real part:
1 x cos x sin x
x cos x = + .
D2 +9 8 32
8. (T) Solve x2 y 00 − 2xy 0 − 4y = x2 + 2 log x, x > 0.
Solution:
Apply the transformation x = et the equation reduces to y 00 − 3y − 4 = e2t + 2t.
Solution of the homogeneous part c1 e4t + c2 e−t .
1 1
Particular integral: D2 −3D−4
(e2t + 2t) = −e2t /6 + D2 −3D−4
2t = − 61 e2t − 12 (t − 3/4).
1 1 1 −1
2t = 2 t = − [1+(D2 /4−3D/4+· · · )]t = (t−3/4).
D2 − 3D − 4 2
−4(1 − D /4 + 3D/4) 2 2
9. (T) (Higher order variation of parameter) Consider the n-th order linear equation
n
X
(n)
y + ai (x)y (i) = y (n) + an−1 (x)y (n−1) + · · · + a0 (x)y = r(x).
1
(n−2)
vi0 yi
P
if = 0.
X (n−1)
X (n)
yp(n) = vi0 yi + vi yi − − − − − −(n + 1).
Then n
X X (n−1)
yp(n) + ai (x)yp(i) = vi0 yi .
1
10. (i) Let y1 (x), y2 (x) are two linearly independent solutions of y 00 +p(x)y 0 +q(x)y = 0. Show
that φ(x) = αy1 (x) + βy2 (x) and ψ(x) = γy1 (x) + δy2 (x) are two linearly independent
solutions if and only if αδ 6= βγ.
(ii) Show that the zeros of the functions a sin x + b cos x and c sin x + d cos x are distinct
and occur alternately whenever ad − bc 6= 0.
Solution:
(i) We have W (φ, ψ) = (αδ − βγ)W (y1 , y2 ). Since y1 , y2 are fundamental solutions,
W (y1 , y2 ) 6= 0. If αδ 6= βγ, then W (φ, ψ) 6= 0. Conversely if W (φ, ψ) 6= 0, then
αδ 6= βγ.
(ii) We know sin x, cos x are independent solutions of y 00 + y = 0. So by part (i) a sin x +
b cos x and c sin x + d cos x are independent solutions whenever ad − bc 6= 0. Hence the
result follows from Sturm Separation theorem ( Simmons, page 190, Theorem A).
11. (T) Show that any nontrivial solution u(x) of u00 + q(x)u = 0, q(x) < 0 for all x, has at
most one zero.
Solution:
Consider the equation z 00 = 0. Then z = 1 is a solution of the equation. By Strum
comparison theorem, between two zeros of u(x) there must be at least one zero of z(x).
But z = 1 has no zero. Hence u(x) can have at most one zero.
12. Let u(x) be any nontrivial solution of u00 + [1 + q(x)]u = 0, where q(x) > 0. Show that
u(x) has infinitely many zeros.
Solution:
Consider
v 00 + v = 0, u00 + 1 + q(x) u = 0
13. Let u(x) be any nontrivial solution of u00 + q(x)u = 0 on a closed interval [a, b]. Show
that u(x) has at most a finite number of zeros in [a, b].
Solution:
Suppose, on the contrary, u(x) has infinite number of zeros in [a, b]. It follows that there
exists x0 ∈ [a, b] and a sequence of zeros xn 6= x0 such that xn → x0 . Since u(x) is
continuous and differentiable at x0 , we have
u(xn ) − u(x0 )
u(x0 ) = lim u(xn ) = 0, u0 (x0 ) = lim =0
xn →x0 xn →x0 xn − x0
By uniqueness theorem, u ≡ 0 which contradicts the fact that u is nontrivial.
x2 y 00 + xy 0 + (x2 − p2 )y = 0, x > 0.
1/4 − p2
u00 + (1 + )u = 0.
x2
2
Given p ≥ 0, we can choose x0 large enough such that 1+ 1/4−p x2
> 1/4 for all x ∈ (x0 , ∞).
00 1
Compare Jp with sin(x/2) which is solution of v + 4 v = 0 in (x0 , ∞). Clearly sin(x/2) has
infinitely many zeros in (x0 , ∞). By Sturm comparison theorem, between two consecutive
zeros of sin(x/2) there is a zero of Jp . Hence Jp has infinitely many zero in (x0 , ∞).
15. (T) Consider u00 + q(x)u = 0 on an interval I = (0, ∞) with q(x) ≥ m2 for all t ∈ I.
Show any non trivial solution u(x) has infinitely many zeros and distance between two
consecutive zeros is at most π/m.
Solution: Compare u(x) with sin mx which is a solution of v 00 + m2 v = 0. By Sturm
comparison theorem, between two consecutive zeros of v(x) = sin(mx) there is a zero of
u(x). Hence u(x) has infinitely many zero in (x0 , ∞).
Let u(a) = 0. We will show that u(x) has a zero in (a, a + π/m]. Consider v(x) =
sin(mx − ma) which is a solution of v 00 + m2 v = 0. Clearly v(a) = v(a + π/m) = 0. Hence
by Sturm comparison theorem, there exists at least one zero of u(x) in (a, a + π/m).
Hence distance between two consecutive zeros of u(x) is at most π/m.
16. Consider u00 + q(x)u = 0 on an interval I = (0, ∞) with q(x) ≤ m2 for all t ∈ I. Show
that distance between two consecutive zeros is at least π/m.
Solution:
Suppose u(a) = 0 and u(b) be two consecutive zeros. Consider v(x) = sin(mx − ma)
which is a solution of v 00 + m2 v = 0. By Sturm comparison theorem, there exists a zero
of v(x) in (a, b). But we know that v(a) = 0 and next zero of v is at a + π/m. So
b > a + π/m.
x2 y 00 + xy 0 + (x2 − p2 )y = 0, x > 0.
Show that (i) If 0 ≤ p ≤ 1/2, then every interval of length π has at least contains at
least one zero of Jp .
(ii) If p = 1/2 then distance between consecutive zeros of Jp is exactly π.
(iii) If p > 1/2 then every interval of length π contains at most one zero of Jp .
Solution: The normal form of Bessel equation is
1/4 − p2
u00 + (1 + )u = 0.
x2
The zeros of Jp and u(x) are same.
(i)Apply exercise 15 with m = 1.
(ii) Clear from normal form.
(iii) Apply exercise 16 with m = 1.
18. Let y(x) be a non-trivial solution of y 00 + q(x)y = 0. Prove that if q(x) > k/x2 for some
k > 1/4 then y has infinitely many positive zeros. If q(x) < 4x12 then y has only finitely
many positive zeros.
Solution:
Consider the Cauchy-Euler equation y 00 + ky x2
= 0. With x = et , it transforms into
y 00 − y 0 + ky = 0. So characteristic equation m2 − m + k = 0. So 1 − 4k = 0 implies two
equal real roots and so the solution has finitely many zeros. If 1 − 4k < 0 then complex
conjugate roots and solution look like xm sin(βx) and it has infinitely many zeros. Rest
follows from Sturm comparison theorem.
Note: This service is not intended for secure transactions such as banking, social media, email, or purchasing. Use at your own risk. We assume no liability whatsoever for broken pages.
Alternative Proxies: