Higher-Order Ordinary Differential Equations — Practice Questions & Solutions
> Perfect for: Semester exam preparation, competitive exams (GATE, ESE), and concept reinforcement.
Module 1: Linear Differential Equation of nth Order
Q1. Solve: d 2 y d x 2 − 5 d y d x + 6 y = 0 \frac{d^2y}{dx^2} - 5\frac{dy}{dx} + 6y = 0 d x 2 d 2 y − 5 d x d y + 6 y = 0
Solution:
This is a second-order linear homogeneous ODE with constant coefficients .
Step 1: Write the auxiliary equation by substituting D = d d x D = \frac{d}{dx} D = d x d :
( D 2 − 5 D + 6 ) y = 0 (D^2 - 5D + 6)y = 0 ( D 2 − 5 D + 6 ) y = 0
m 2 − 5 m + 6 = 0 m^2 - 5m + 6 = 0 m 2 − 5 m + 6 = 0
Step 2: Solve the quadratic:
( m − 2 ) ( m − 3 ) = 0 (m - 2)(m - 3) = 0 ( m − 2 ) ( m − 3 ) = 0
m 1 = 2 , m 2 = 3 m_1 = 2, \quad m_2 = 3 m 1 = 2 , m 2 = 3
Step 3: Since roots are real and distinct , the complementary function is:
y c = C 1 e 2 x + C 2 e 3 x \boxed{y_c = C_1 e^{2x} + C_2 e^{3x}} y c = C 1 e
Since f ( x ) = 0 f(x) = 0 f ( x ) = 0 , the particular integral y p = 0 y_p = 0 y p = 0 .
Final Answer: y = C 1 e 2 x + C 2 e 3 x y = C_1 e^{2x} + C_2 e^{3x} y = C 1 e 2 x + C 2
Q2. Solve: d 2 y d x 2 + 4 d y d x + 4 y = 0 \frac{d^2y}{dx^2} + 4\frac{dy}{dx} + 4y = 0 d x 2 d 2 y +
Solution:
Step 1: Auxiliary equation:
m 2 + 4 m + 4 = 0 m^2 + 4m + 4 = 0 m 2 + 4 m + 4 = 0
( m + 2 ) 2 = 0 (m + 2)^2 = 0 ( m + 2 ) 2 = 0
m = − 2 , − 2 (repeated root) m = -2, -2 \quad \text{(repeated root)} m = − 2 , − 2 (repeated root)
Step 2: For real and equal roots , the solution is:
y = ( C 1 + C 2 x ) e − 2 x \boxed{y = (C_1 + C_2 x)e^{-2x}} y = ( C 1 + C 2 x
Q3. Solve: d 2 y d x 2 + 4 y = 0 \frac{d^2y}{dx^2} + 4y = 0 d x 2 d 2 y + 4 y
Solution:
Step 1: Auxiliary equation:
m 2 + 4 = 0 m^2 + 4 = 0 m 2 + 4 = 0
m 2 = − 4 m^2 = -4 m 2 = − 4
m = ± 2 i (complex conjugate roots, α = 0 , β = 2 ) m = \pm 2i \quad \text{(complex conjugate roots, } \alpha = 0, \beta = 2\text{)} m = ± 2 i (complex conjugate roots, α = 0 , β = 2 )
Step 2: For complex roots α ± i β \alpha \pm i\beta α ± i β :
y = e 0 x ( C 1 cos 2 x + C 2 sin 2 x ) = C 1 cos 2 x + C 2 sin 2 x \boxed{y = e^{0x}(C_1 \cos 2x + C_2 \sin 2x) = C_1 \cos 2x + C_2 \sin 2x} y = e 0 x ( C
Q4. Solve: d 3 y d x 3 − 6 d 2 y d x 2 + 11 d y d x − 6 y = 0 \frac{d^3y}{dx^3} - 6\frac{d^2y}{dx^2} + 11\frac{dy}{dx} - 6y = 0 d x 3 d 3 y
Solution:
Step 1: Auxiliary equation:
m 3 − 6 m 2 + 11 m − 6 = 0 m^3 - 6m^2 + 11m - 6 = 0 m 3 − 6 m 2 + 11 m − 6 = 0
Step 2: By trial, m = 1 m = 1 m = 1 is a root. Factorize:
( m − 1 ) ( m 2 − 5 m + 6 ) = 0 (m - 1)(m^2 - 5m + 6) = 0 ( m − 1 ) ( m 2 − 5 m + 6 ) = 0
( m − 1 ) ( m − 2 ) ( m − 3 ) = 0 (m - 1)(m - 2)(m - 3) = 0 ( m − 1 ) ( m − 2 ) ( m − 3 ) = 0
m = 1 , 2 , 3 m = 1, 2, 3 m = 1 , 2 , 3
Step 3: All roots real and distinct:
y = C 1 e x + C 2 e 2 x + C 3 e 3 x \boxed{y = C_1 e^x + C_2 e^{2x} + C_3 e^{3x}} y = C 1 e x + C
Q5. Solve: d 2 y d x 2 − 3 d y d x + 2 y = e 3 x \frac{d^2y}{dx^2} - 3\frac{dy}{dx} + 2y = e^{3x} d x 2 d 2 y
Solution:
Step 1: Auxiliary equation:
m 2 − 3 m + 2 = 0 m^2 - 3m + 2 = 0 m 2 − 3 m + 2 = 0
( m − 1 ) ( m − 2 ) = 0 (m - 1)(m - 2) = 0 ( m − 1 ) ( m − 2 ) = 0
m = 1 , 2 m = 1, 2 m = 1 , 2
Complementary Function:
y c = C 1 e x + C 2 e 2 x y_c = C_1 e^x + C_2 e^{2x} y c = C 1 e x + C
Step 2: Find Particular Integral using operator method:
y p = 1 D 2 − 3 D + 2 e 3 x y_p = \frac{1}{D^2 - 3D + 2} e^{3x} y p = D 2 − 3 D + 2
Substitute D = 3 D = 3 D = 3 :
y p = 1 9 − 9 + 2 e 3 x = 1 2 e 3 x y_p = \frac{1}{9 - 9 + 2} e^{3x} = \frac{1}{2}e^{3x} y p = 9 − 9 + 2 1
Final Answer:
y = C 1 e x + C 2 e 2 x + 1 2 e 3 x \boxed{y = C_1 e^x + C_2 e^{2x} + \frac{1}{2}e^{3x}} y = C 1 e x + C
Q6. Solve: d 2 y d x 2 + y = sin x \frac{d^2y}{dx^2} + y = \sin x d x 2 d 2 y + y
Solution:
Step 1: Auxiliary equation:
m 2 + 1 = 0 ⇒ m = ± i m^2 + 1 = 0 \Rightarrow m = \pm i m 2 + 1 = 0 ⇒ m = ± i
y c = C 1 cos x + C 2 sin x y_c = C_1 \cos x + C_2 \sin x y c = C 1 cos x + C
Step 2: Particular Integral (resonance case — f ( x ) = sin x f(x) = \sin x f ( x ) = sin x matches CF):
y p = 1 D 2 + 1 sin x y_p = \frac{1}{D^2 + 1} \sin x y p = D 2 + 1 1
Since D 2 + 1 = 0 D^2 + 1 = 0 D 2 + 1 = 0 when D 2 = − 1 D^2 = -1 D 2 = − 1 , we use the rule for failure case:
y p = x ⋅ 1 2 D sin x = x 2 ∫ sin x d x = − x 2 cos x y_p = x \cdot \frac{1}{2D} \sin x = \frac{x}{2} \int \sin x \, dx = -\frac{x}{2}\cos x y p = x ⋅ 2 D
Final Answer:
y = C 1 cos x + C 2 sin x − x 2 cos x \boxed{y = C_1 \cos x + C_2 \sin x - \frac{x}{2}\cos x} y = C 1 cos x + C
Q7. Solve: d 2 y d x 2 − 2 d y d x + y = x 2 e x \frac{d^2y}{dx^2} - 2\frac{dy}{dx} + y = x^2 e^x d x 2 d 2 y
Solution:
Step 1: Auxiliary equation:
m 2 − 2 m + 1 = 0 ⇒ ( m − 1 ) 2 = 0 ⇒ m = 1 , 1 m^2 - 2m + 1 = 0 \Rightarrow (m-1)^2 = 0 \Rightarrow m = 1, 1 m 2 − 2 m + 1 = 0 ⇒ ( m − 1 )
y c = ( C 1 + C 2 x ) e x y_c = (C_1 + C_2 x)e^x y c = ( C 1 + C 2
Step 2: Particular Integral:
y p = 1 ( D − 1 ) 2 x 2 e x = e x 1 D 2 x 2 = e x ⋅ x 4 12 y_p = \frac{1}{(D-1)^2} x^2 e^x = e^x \frac{1}{D^2} x^2 = e^x \cdot \frac{x^4}{12} y p = ( D − 1 )
Final Answer:
y = ( C 1 + C 2 x ) e x + x 4 12 e x \boxed{y = (C_1 + C_2 x)e^x + \frac{x^4}{12}e^x} y = ( C 1 + C 2
Q8. Solve: d 3 y d x 3 − 7 d y d x + 6 y = 1 + x + x 2 \frac{d^3y}{dx^3} - 7\frac{dy}{dx} + 6y = 1 + x + x^2 d x 3 d 3 y
Solution:
Step 1: Auxiliary equation:
m 3 − 7 m + 6 = 0 m^3 - 7m + 6 = 0 m 3 − 7 m + 6 = 0
By trial m = 1 m = 1 m = 1 : 1 − 7 + 6 = 0 1 - 7 + 6 = 0 1 − 7 + 6 = 0 ✓
( m − 1 ) ( m 2 + m − 6 ) = 0 (m - 1)(m^2 + m - 6) = 0 ( m − 1 ) ( m 2 + m − 6 ) = 0
( m − 1 ) ( m + 3 ) ( m − 2 ) = 0 (m - 1)(m + 3)(m - 2) = 0 ( m − 1 ) ( m + 3 ) ( m − 2 ) = 0
m = 1 , 2 , − 3 m = 1, 2, -3 m = 1 , 2 , − 3
y c = C 1 e x + C 2 e 2 x + C 3 e − 3 x y_c = C_1 e^x + C_2 e^{2x} + C_3 e^{-3x} y c = C 1 e x +
Step 2: Particular Integral for polynomial:
y p = 1 D 3 − 7 D + 6 ( 1 + x + x 2 ) y_p = \frac{1}{D^3 - 7D + 6}(1 + x + x^2) y p = D 3 − 7 D + 6
Using binomial expansion or long division:
y p = 1 6 ( 1 + 7 D − D 3 6 ) − 1 ( 1 + x + x 2 ) y_p = \frac{1}{6}\left(1 + \frac{7D - D^3}{6}\right)^{-1}(1 + x + x^2) y p = 6 1
= 1 6 ( 1 − 7 D 6 + 49 D 2 36 + ⋯ ) ( 1 + x + x 2 ) = \frac{1}{6}\left(1 - \frac{7D}{6} + \frac{49D^2}{36} + \cdots\right)(1 + x + x^2) = 6 1 ( 1 −
Apply to 1 + x + x 2 1 + x + x^2 1 + x + x 2 :
D ( 1 + x + x 2 ) = 1 + 2 x D(1 + x + x^2) = 1 + 2x D ( 1 + x + x 2 ) = 1 + 2 x
D 2 ( 1 + x + x 2 ) = 2 D^2(1 + x + x^2) = 2
y p = 1 6 [ ( 1 + x + x 2 ) − 7 6 ( 1 + 2 x ) + 49 36 ( 2 ) ] y_p = \frac{1}{6}\left[(1 + x + x^2) - \frac{7}{6}(1 + 2x) + \frac{49}{36}(2)\right] y p = 6 1
= 1 6 [ 1 + x + x 2 − 7 6 − 14 x 6 + 98 36 ] = \frac{1}{6}\left[1 + x + x^2 - \frac{7}{6} - \frac{14x}{6} + \frac{98}{36}\right] = 6 1 [ 1 + x + x
Simplify:
y p = x 2 6 − 4 x 36 + 50 216 = x 2 6 − x 9 + 25 108 y_p = \frac{x^2}{6} - \frac{4x}{36} + \frac{50}{216} = \frac{x^2}{6} - \frac{x}{9} + \frac{25}{108} y p = 6 x
Final Answer:
y = C 1 e x + C 2 e 2 x + C 3 e − 3 x + x 2 6 − x 9 + 25 108 \boxed{y = C_1 e^x + C_2 e^{2x} + C_3 e^{-3x} + \frac{x^2}{6} - \frac{x}{9} + \frac{25}{108}} y = C 1 e
Module 2: Simultaneous Linear Differential Equations
Q9. Solve the simultaneous equations:
d x d t + y = sin t \frac{dx}{dt} + y = \sin t d t d x + y = sin t
d y d t + x = cos t \frac{dy}{dt} + x = \cos t d t d y + x = cos t
Solution:
Step 1: Write in operator form (D = d d t D = \frac{d}{dt} D = d t d ):
( 1 ) D x + y = sin t (1) \quad Dx + y = \sin t ( 1 ) D x + y = sin t
( 2 ) x + D y = cos t (2) \quad x + Dy = \cos t ( 2 ) x + D y =
Step 2: Eliminate y y y . From (1): y = sin t − D x y = \sin t - Dx y = sin t − D x
Substitute into (2):
x + D ( sin t − D x ) = cos t x + D(\sin t - Dx) = \cos t x + D ( sin t − D x ) = cos t
x + cos t − D 2 x = cos t x + \cos t - D^2x = \cos t x + cos t − D 2 x = cos t
( 1 − D 2 ) x = 0 (1 - D^2)x = 0 ( 1 − D 2 ) x = 0
D 2 x − x = 0 ⇒ ( D 2 − 1 ) x = 0 D^2x - x = 0 \Rightarrow (D^2 - 1)x = 0 D 2 x − x = 0 ⇒ ( D 2 − 1 ) x =
Step 3: Solve for x x x :
m 2 − 1 = 0 ⇒ m = ± 1 m^2 - 1 = 0 \Rightarrow m = \pm 1 m 2 − 1 = 0 ⇒ m = ± 1
x = C 1 e t + C 2 e − t x = C_1 e^t + C_2 e^{-t} x = C 1 e t + C 2 e
Step 4: Find y y y from equation (1):
y = sin t − D x = sin t − ( C 1 e t − C 2 e − t ) y = \sin t - Dx = \sin t - (C_1 e^t - C_2 e^{-t}) y = sin t − D x = sin t − ( C 1
y = sin t − C 1 e t + C 2 e − t y = \sin t - C_1 e^t + C_2 e^{-t} y = sin t − C 1 e t + C
Final Answer:
x = C 1 e t + C 2 e − t , y = sin t − C 1 e t + C 2 e − t \boxed{x = C_1 e^t + C_2 e^{-t}, \quad y = \sin t - C_1 e^t + C_2 e^{-t}} x = C 1 e
Q10. Solve:
d x d t + 2 x − 3 y = t \frac{dx}{dt} + 2x - 3y = t d t d x + 2 x − 3 y = t
d y d t − 3 x + 2 y = e 2 t \frac{dy}{dt} - 3x + 2y = e^{2t} d t d y − 3 x + 2 y = e
Solution:
Step 1: Operator form:
( D + 2 ) x − 3 y = t ...(1) (D + 2)x - 3y = t \quad \text{...(1)} ( D + 2 ) x − 3 y = t ...(1)
− 3 x + ( D + 2 ) y = e 2 t ...(2) -3x + (D + 2)y = e^{2t} \quad \text{...(2)} − 3 x + ( D + 2 ) y = e 2 t ...(2)
Step 2: Eliminate y y y . Multiply (1) by ( D + 2 ) (D+2) ( D + 2 ) and (2) by 3 3 3 :
( D + 2 ) 2 x − 3 ( D + 2 ) y = ( D + 2 ) t = 1 + 2 t (D+2)^2 x - 3(D+2)y = (D+2)t = 1 + 2t ( D + 2 ) 2 x − 3 ( D + 2 ) y = ( D
− 9 x + 3 ( D + 2 ) y = 3 e 2 t -9x + 3(D+2)y = 3e^{2t} − 9 x + 3 ( D + 2 ) y = 3 e 2 t
Add:
( D + 2 ) 2 x − 9 x = 1 + 2 t + 3 e 2 t (D+2)^2 x - 9x = 1 + 2t + 3e^{2t} ( D + 2 ) 2 x − 9 x = 1 + 2 t + 3 e
( D 2 + 4 D + 4 − 9 ) x = 1 + 2 t + 3 e 2 t (D^2 + 4D + 4 - 9)x = 1 + 2t + 3e^{2t} ( D 2 + 4 D + 4 − 9 ) x = 1 + 2 t
( D 2 + 4 D − 5 ) x = 1 + 2 t + 3 e 2 t (D^2 + 4D - 5)x = 1 + 2t + 3e^{2t} ( D 2 + 4 D − 5 ) x = 1 + 2 t + 3 e
Step 3: Solve homogeneous part:
m 2 + 4 m − 5 = 0 ⇒ ( m + 5 ) ( m − 1 ) = 0 ⇒ m = 1 , − 5 m^2 + 4m - 5 = 0 \Rightarrow (m+5)(m-1) = 0 \Rightarrow m = 1, -5 m 2 + 4 m − 5 = 0 ⇒ ( m + 5 ) (
x c = C 1 e t + C 2 e − 5 t x_c = C_1 e^t + C_2 e^{-5t} x c = C 1 e t + C
Step 4: Particular Integral:
x p = 1 D 2 + 4 D − 5 ( 1 + 2 t + 3 e 2 t ) x_p = \frac{1}{D^2 + 4D - 5}(1 + 2t + 3e^{2t}) x p = D 2 + 4 D − 5
For 1 + 2 t 1 + 2t 1 + 2 t :
1 − 5 ( 1 − 4 D + D 2 5 ) − 1 ( 1 + 2 t ) = − 1 5 ( 1 + 4 D 5 ) ( 1 + 2 t ) \frac{1}{-5}\left(1 - \frac{4D + D^2}{5}\right)^{-1}(1 + 2t) = -\frac{1}{5}\left(1 + \frac{4D}{5}\right)(1 + 2t) − 5 1 ( 1 −
= − 1 5 ( 1 + 2 t + 8 5 ) = − 1 5 ( 2 t + 13 5 ) = − 2 t 5 − 13 25 = -\frac{1}{5}\left(1 + 2t + \frac{8}{5}\right) = -\frac{1}{5}\left(2t + \frac{13}{5}\right) = -\frac{2t}{5} - \frac{13}{25} = − 5 1 ( 1 + 2 t
For 3 e 2 t 3e^{2t} 3 e 2 t :
3 e 2 t 4 + 8 − 5 = 3 e 2 t 7 \frac{3e^{2t}}{4 + 8 - 5} = \frac{3e^{2t}}{7} 4 + 8 − 5 3 e 2 t =
x p = − 2 t 5 − 13 25 + 3 e 2 t 7 x_p = -\frac{2t}{5} - \frac{13}{25} + \frac{3e^{2t}}{7} x p = − 5 2 t −
Step 5: Find y y y from equation (1):
3 y = ( D + 2 ) x − t 3y = (D+2)x - t 3 y = ( D + 2 ) x − t
Compute ( D + 2 ) x (D+2)x ( D + 2 ) x :
x = C 1 e t + C 2 e − 5 t − 2 t 5 − 13 25 + 3 e 2 t 7 x = C_1 e^t + C_2 e^{-5t} - \frac{2t}{5} - \frac{13}{25} + \frac{3e^{2t}}{7} x = C 1 e t + C
D x = C 1 e t − 5 C 2 e − 5 t − 2 5 + 6 e 2 t 7 Dx = C_1 e^t - 5C_2 e^{-5t} - \frac{2}{5} + \frac{6e^{2t}}{7} D x = C 1 e t − 5 C
( D + 2 ) x = 3 C 1 e t − 3 C 2 e − 5 t − 2 5 − 4 t 5 − 26 25 + 9 e 2 t 7 (D+2)x = 3C_1 e^t - 3C_2 e^{-5t} - \frac{2}{5} - \frac{4t}{5} - \frac{26}{25} + \frac{9e^{2t}}{7} ( D + 2 ) x = 3 C 1 e
3 y = 3 C 1 e t − 3 C 2 e − 5 t − 4 t 5 − 36 25 + 9 e 2 t 7 − t 3y = 3C_1 e^t - 3C_2 e^{-5t} - \frac{4t}{5} - \frac{36}{25} + \frac{9e^{2t}}{7} - t 3 y = 3 C 1 e t −
y = C 1 e t − C 2 e − 5 t − 3 t 5 − 12 25 + 3 e 2 t 7 \boxed{y = C_1 e^t - C_2 e^{-5t} - \frac{3t}{5} - \frac{12}{25} + \frac{3e^{2t}}{7}} y = C 1 e
Module 3: Second-Order Linear ODEs with Variable Coefficients
Q11. Solve: x 2 d 2 y d x 2 + x d y d x − 4 y = 0 x^2\frac{d^2y}{dx^2} + x\frac{dy}{dx} - 4y = 0 x 2 d x 2 d
Solution:
This is a Cauchy-Euler equation . Substitute x = e t x = e^t x = e t (or t = ln x t = \ln x t = ln x ):
x d y d x = D y , x 2 d 2 y d x 2 = D ( D − 1 ) y x\frac{dy}{dx} = Dy, \quad x^2\frac{d^2y}{dx^2} = D(D-1)y x d x d y = D y , x
where D = d d t D = \frac{d}{dt} D = d t d .
Step 1: Transform:
[ D ( D − 1 ) + D − 4 ] y = 0 [D(D-1) + D - 4]y = 0 [ D ( D − 1 ) + D − 4 ] y = 0
( D 2 − D + D − 4 ) y = 0 (D^2 - D + D - 4)y = 0 ( D 2 − D + D − 4 ) y = 0
( D 2 − 4 ) y = 0 (D^2 - 4)y = 0 ( D 2 − 4 ) y = 0
Step 2: Solve:
m 2 − 4 = 0 ⇒ m = ± 2 m^2 - 4 = 0 \Rightarrow m = \pm 2 m 2 − 4 = 0 ⇒ m = ± 2
y = C 1 e 2 t + C 2 e − 2 t y = C_1 e^{2t} + C_2 e^{-2t} y = C 1 e 2 t + C 2
Step 3: Substitute back t = ln x t = \ln x t = ln x , so e t = x e^t = x e t = x :
y = C 1 x 2 + C 2 x 2 \boxed{y = C_1 x^2 + \frac{C_2}{x^2}} y = C 1 x 2 +
Q12. Solve: x 2 d 2 y d x 2 + 3 x d y d x + y = 0 x^2\frac{d^2y}{dx^2} + 3x\frac{dy}{dx} + y = 0 x 2 d x 2 d
Solution:
Step 1: Substitute x = e t x = e^t x = e t :
[ D ( D − 1 ) + 3 D + 1 ] y = 0 [D(D-1) + 3D + 1]y = 0 [ D ( D − 1 ) + 3 D + 1 ] y = 0
( D 2 − D + 3 D + 1 ) y = 0 (D^2 - D + 3D + 1)y = 0 ( D 2 − D + 3 D + 1 ) y = 0
( D 2 + 2 D + 1 ) y = 0 (D^2 + 2D + 1)y = 0 ( D 2 + 2 D + 1 ) y = 0
( D + 1 ) 2 y = 0 (D+1)^2 y = 0 ( D + 1 ) 2 y = 0
Step 2: Repeated root m = − 1 , − 1 m = -1, -1 m = − 1 , − 1 :
y = ( C 1 + C 2 t ) e − t y = (C_1 + C_2 t)e^{-t} y = ( C 1 + C 2 t ) e
Step 3: Substitute t = ln x t = \ln x t = ln x :
y = C 1 + C 2 ln x x \boxed{y = \frac{C_1 + C_2 \ln x}{x}} y = x C 1 + C
Module 4: Solution by Changing the Independent Variable
Q13. Solve: d 2 y d x 2 + tan x d y d x + y cos 2 x = 0 \frac{d^2y}{dx^2} + \tan x \frac{dy}{dx} + y \cos^2 x = 0 d x 2 d 2 y
Solution:
Step 1: Compare with standard form d 2 y d x 2 + P d y d x + Q y = R \frac{d^2y}{dx^2} + P\frac{dy}{dx} + Qy = R d x 2 d 2 y + :
P = tan x , Q = cos 2 x , R = 0 P = \tan x, \quad Q = \cos^2 x, \quad R = 0 P = tan x , Q = cos 2 x , R = 0
Step 2: Choose new independent variable z z z such that:
( d z d x ) 2 = Q = cos 2 x \left(\frac{dz}{dx}\right)^2 = Q = \cos^2 x ( d x d z ) 2 = Q =
d z d x = cos x ⇒ z = sin x \frac{dz}{dx} = \cos x \Rightarrow z = \sin x d x d z = cos x ⇒ z = sin x
Step 3: Transform the equation. We have:
d 2 y d z 2 + P 1 d y d z + Q 1 y = R 1 \frac{d^2y}{dz^2} + P_1 \frac{dy}{dz} + Q_1 y = R_1 d z 2 d 2 y
where:
Q 1 = Q ( d z / d x ) 2 = cos 2 x cos 2 x = 1 Q_1 = \frac{Q}{(dz/dx)^2} = \frac{\cos^2 x}{\cos^2 x} = 1 Q 1 = ( d z / d x )
P 1 = d 2 z d x 2 + P d z d x ( d z / d x ) 2 = − sin x + tan x ⋅ cos x cos 2 x = − sin x + sin x cos 2 x = 0 P_1 = \frac{\frac{d^2z}{dx^2} + P\frac{dz}{dx}}{(dz/dx)^2} = \frac{-\sin x + \tan x \cdot \cos x}{\cos^2 x} = \frac{-\sin x + \sin x}{\cos^2 x} = 0 P 1 =
Step 4: The transformed equation is:
d 2 y d z 2 + y = 0 \frac{d^2y}{dz^2} + y = 0 d z 2 d 2 y + y
Step 5: Solve:
m 2 + 1 = 0 ⇒ m = ± i m^2 + 1 = 0 \Rightarrow m = \pm i m 2 + 1 = 0 ⇒ m = ± i
y = C 1 cos z + C 2 sin z y = C_1 \cos z + C_2 \sin z y = C 1 cos z + C 2 sin z
Final Answer:
y = C 1 cos ( sin x ) + C 2 sin ( sin x ) \boxed{y = C_1 \cos(\sin x) + C_2 \sin(\sin x)} y = C 1 cos ( sin x ) + C
Module 5: Method of Variation of Parameters
Q14. Solve by variation of parameters: d 2 y d x 2 + y = tan x \frac{d^2y}{dx^2} + y = \tan x d x 2 d 2 y +
Solution:
Step 1: Complementary function:
m 2 + 1 = 0 ⇒ m = ± i m^2 + 1 = 0 \Rightarrow m = \pm i m 2 + 1 = 0 ⇒ m = ± i
y c = C 1 cos x + C 2 sin x y_c = C_1 \cos x + C_2 \sin x y c = C 1 cos x + C
So y 1 = cos x y_1 = \cos x y 1 = cos x , y 2 = sin x y_2 = \sin x y 2 = sin x .
Step 2: Compute Wronskian:
W = ∣ y 1 y 2 y 1 ′ y 2 ′ ∣ = ∣ cos x sin x − sin x cos x ∣ = cos 2 x + sin 2 x = 1 W = \begin{vmatrix} y_1 & y_2 \\ y_1' & y_2' \end{vmatrix} = \begin{vmatrix} \cos x & \sin x \\ -\sin x & \cos x \end{vmatrix} = \cos^2 x + \sin^2 x = 1 W =
Step 3: Find particular solution:
y p = − y 1 ∫ y 2 f ( x ) W d x + y 2 ∫ y 1 f ( x ) W d x y_p = -y_1 \int \frac{y_2 f(x)}{W} dx + y_2 \int \frac{y_1 f(x)}{W} dx y p = − y 1 ∫
where f ( x ) = tan x f(x) = \tan x f ( x ) = tan x .
First integral:
∫ sin x ⋅ tan x 1 d x = ∫ sin 2 x cos x d x = ∫ 1 − cos 2 x cos x d x \int \frac{\sin x \cdot \tan x}{1} dx = \int \frac{\sin^2 x}{\cos x} dx = \int \frac{1 - \cos^2 x}{\cos x} dx ∫ 1 s i n x ⋅ t a n x
= ∫ ( sec x − cos x ) d x = ln ∣ sec x + tan x ∣ − sin x = \int (\sec x - \cos x) dx = \ln|\sec x + \tan x| - \sin x = ∫ ( sec x − cos x ) d x = ln ∣ sec x + tan x ∣ −
Second integral:
∫ cos x ⋅ tan x 1 d x = ∫ sin x d x = − cos x \int \frac{\cos x \cdot \tan x}{1} dx = \int \sin x \, dx = -\cos x ∫ 1 c o s x ⋅ t a n x d x = ∫
Step 4: Compute y p y_p y p :
y p = − cos x [ ln ∣ sec x + tan x ∣ − sin x ] + sin x [ − cos x ] y_p = -\cos x [\ln|\sec x + \tan x| - \sin x] + \sin x [-\cos x] y p = − cos x [ ln ∣ sec x + tan x
= − cos x ln ∣ sec x + tan x ∣ + cos x sin x − sin x cos x = -\cos x \ln|\sec x + \tan x| + \cos x \sin x - \sin x \cos x = − cos x ln ∣ sec x + tan x ∣ + cos x sin x −
= − cos x ln ∣ sec x + tan x ∣ = -\cos x \ln|\sec x + \tan x| = − cos x ln ∣ sec x + tan x ∣
Final Answer:
y = C 1 cos x + C 2 sin x − cos x ln ∣ sec x + tan x ∣ \boxed{y = C_1 \cos x + C_2 \sin x - \cos x \ln|\sec x + \tan x|} y = C 1 cos x + C
Q15. Solve by variation of parameters: d 2 y d x 2 − 2 d y d x + y = e x x \frac{d^2y}{dx^2} - 2\frac{dy}{dx} + y = \frac{e^x}{x} d x 2 d 2 y
Solution:
Step 1: CF: m 2 − 2 m + 1 = 0 ⇒ ( m − 1 ) 2 = 0 ⇒ m = 1 , 1 m^2 - 2m + 1 = 0 \Rightarrow (m-1)^2 = 0 \Rightarrow m = 1, 1 m 2 − 2 m + 1 = 0 ⇒ ( m − 1 )
For repeated roots in variation of parameters, we use:
y 1 = e x , y 2 = x e x y_1 = e^x, \quad y_2 = xe^x y 1 = e x , y 2 =
Step 2: Wronskian:
W = ∣ e x x e x e x e x + x e x ∣ = e x ( e x + x e x ) − x e x ⋅ e x = e 2 x W = \begin{vmatrix} e^x & xe^x \\ e^x & e^x + xe^x \end{vmatrix} = e^x(e^x + xe^x) - xe^x \cdot e^x = e^{2x} W =
Step 3: Particular integral:
y p = − e x ∫ x e x ⋅ e x x e 2 x d x + x e x ∫ e x ⋅ e x x e 2 x d x y_p = -e^x \int \frac{xe^x \cdot \frac{e^x}{x}}{e^{2x}} dx + xe^x \int \frac{e^x \cdot \frac{e^x}{x}}{e^{2x}} dx y p = − e x ∫
= − e x ∫ 1 d x + x e x ∫ 1 x d x = -e^x \int 1 \, dx + xe^x \int \frac{1}{x} dx = − e x ∫ 1 d x + x e x ∫
= − x e x + x e x ln ∣ x ∣ = -xe^x + xe^x \ln|x| = − x e x + x e x ln ∣ x ∣
Final Answer:
y = ( C 1 + C 2 x ) e x − x e x + x e x ln ∣ x ∣ = ( C 1 + C 2 x ) e x + x e x ( ln ∣ x ∣ − 1 ) \boxed{y = (C_1 + C_2 x)e^x - xe^x + xe^x \ln|x| = (C_1 + C_2 x)e^x + xe^x(\ln|x| - 1)} y = ( C
Module 6: Cauchy-Euler Equation
Q16. Solve: x 2 d 2 y d x 2 + x d y d x + y = sin ( ln x ) x^2\frac{d^2y}{dx^2} + x\frac{dy}{dx} + y = \sin(\ln x) x 2 d x 2 d
Solution:
Step 1: Substitute x = e t x = e^t x = e t , so t = ln x t = \ln x t = ln x :
[ D ( D − 1 ) + D + 1 ] y = sin t [D(D-1) + D + 1]y = \sin t [ D ( D − 1 ) + D + 1 ] y = sin t
( D 2 + 1 ) y = sin t (D^2 + 1)y = \sin t ( D 2 + 1 ) y = sin t
Step 2: CF: m 2 + 1 = 0 ⇒ m = ± i m^2 + 1 = 0 \Rightarrow m = \pm i m 2 + 1 = 0 ⇒ m = ± i
y c = C 1 cos t + C 2 sin t y_c = C_1 \cos t + C_2 \sin t y c = C 1 cos t + C
Step 3: PI (resonance case):
y p = 1 D 2 + 1 sin t = − t 2 cos t y_p = \frac{1}{D^2 + 1} \sin t = -\frac{t}{2}\cos t y p = D 2 + 1
Step 4: Substitute t = ln x t = \ln x t = ln x :
y = C 1 cos ( ln x ) + C 2 sin ( ln x ) − ln x 2 cos ( ln x ) \boxed{y = C_1 \cos(\ln x) + C_2 \sin(\ln x) - \frac{\ln x}{2}\cos(\ln x)} y = C 1 cos ( ln x )
Q17. Solve: x 2 d 2 y d x 2 − 2 x d y d x − 4 y = x 2 x^2\frac{d^2y}{dx^2} - 2x\frac{dy}{dx} - 4y = x^2 x 2 d x 2 d
Solution:
Step 1: Substitute x = e t x = e^t x = e t :
[ D ( D − 1 ) − 2 D − 4 ] y = e 2 t [D(D-1) - 2D - 4]y = e^{2t} [ D ( D − 1 ) − 2 D − 4 ] y = e 2 t
( D 2 − 3 D − 4 ) y = e 2 t (D^2 - 3D - 4)y = e^{2t} ( D 2 − 3 D − 4 ) y = e 2 t
Step 2: CF: m 2 − 3 m − 4 = 0 ⇒ ( m − 4 ) ( m + 1 ) = 0 m^2 - 3m - 4 = 0 \Rightarrow (m-4)(m+1) = 0 m 2 − 3 m − 4 = 0 ⇒ ( m − 4 ) ( m +
y c = C 1 e 4 t + C 2 e − t y_c = C_1 e^{4t} + C_2 e^{-t} y c = C 1 e 4 t +
Step 3: PI:
y p = 1 D 2 − 3 D − 4 e 2 t = e 2 t 4 − 6 − 4 = − e 2 t 6 y_p = \frac{1}{D^2 - 3D - 4} e^{2t} = \frac{e^{2t}}{4 - 6 - 4} = -\frac{e^{2t}}{6} y p = D
Step 4: Substitute back:
y = C 1 x 4 + C 2 x − x 2 6 \boxed{y = C_1 x^4 + \frac{C_2}{x} - \frac{x^2}{6}} y = C 1 x 4 +
Q18. Solve: ( x 2 D 2 + x D + 1 ) y = ln x ⋅ sin ( ln x ) (x^2D^2 + xD + 1)y = \ln x \cdot \sin(\ln x) ( x 2 D 2 + x D + 1 ) y = ln x
Solution:
Step 1: Substitute x = e t x = e^t x = e t :
( D 2 + 1 ) y = t sin t (D^2 + 1)y = t \sin t ( D 2 + 1 ) y = t sin t
Step 2: CF: y c = C 1 cos t + C 2 sin t y_c = C_1 \cos t + C_2 \sin t y c = C 1 cos t + C
Step 3: PI using complex method:
Consider 1 D 2 + 1 t e i t \frac{1}{D^2 + 1} t e^{it} D 2 + 1 1 t e i t and take imaginary part.
Using 1 D 2 + 1 = 1 ( D + i ) ( D − i ) \frac{1}{D^2+1} = \frac{1}{(D+i)(D-i)} D 2 + 1 1 = and shifting:
1 D 2 + 1 t e i t = e i t 1 ( D + i ) 2 + 1 t = e i t 1 D 2 + 2 i D t \frac{1}{D^2+1} t e^{it} = e^{it} \frac{1}{(D+i)^2 + 1} t = e^{it} \frac{1}{D^2 + 2iD} t D 2 + 1 1
= e i t 1 2 i D ( 1 + D 2 i ) − 1 t = e i t 1 2 i D ( t − 1 2 i ) = e^{it} \frac{1}{2iD}\left(1 + \frac{D}{2i}\right)^{-1} t = e^{it} \frac{1}{2iD}\left(t - \frac{1}{2i}\right) = e i t 2 i D
= e i t 1 2 i ( t 2 2 − t 2 i ) = e i t ( − i t 2 4 + t 4 ) = e^{it} \frac{1}{2i}\left(\frac{t^2}{2} - \frac{t}{2i}\right) = e^{it}\left(-\frac{it^2}{4} + \frac{t}{4}\right) = e i t 2 i 1
= ( cos t + i sin t ) ( t 4 − i t 2 4 ) = (\cos t + i\sin t)\left(\frac{t}{4} - \frac{it^2}{4}\right) = ( cos t + i sin t ) ( 4 t
Imaginary part:
= − t 2 4 cos t + t 4 sin t = -\frac{t^2}{4}\cos t + \frac{t}{4}\sin t = − 4 t 2 cos t +
Step 4: Substitute t = ln x t = \ln x t = ln x :
y = C 1 cos ( ln x ) + C 2 sin ( ln x ) − ( ln x ) 2 4 cos ( ln x ) + ln x 4 sin ( ln x ) \boxed{y = C_1 \cos(\ln x) + C_2 \sin(\ln x) - \frac{(\ln x)^2}{4}\cos(\ln x) + \frac{\ln x}{4}\sin(\ln x)} y = C 1
Module 7: Applications of Differential Equations
Q19. Mechanical Vibrations — Free Damped Motion
A mass of 2 2 2 kg is attached to a spring with constant k = 128 k = 128 k = 128 N/m. The damping coefficient is c = 16 c = 16 c = 16 N·s/m. The mass is displaced 0.1 0.1 0.1 m and released. Find the equation of motion.
Solution:
Step 1: The equation of motion is:
m d 2 x d t 2 + c d x d t + k x = 0 m\frac{d^2x}{dt^2} + c\frac{dx}{dt} + kx = 0 m d t 2 d 2 x +
2 d 2 x d t 2 + 16 d x d t + 128 x = 0 2\frac{d^2x}{dt^2} + 16\frac{dx}{dt} + 128x = 0 2 d t 2 d 2 x +
d 2 x d t 2 + 8 d x d t + 64 x = 0 \frac{d^2x}{dt^2} + 8\frac{dx}{dt} + 64x = 0 d t 2 d 2 x +
Step 2: Auxiliary equation:
m 2 + 8 m + 64 = 0 m^2 + 8m + 64 = 0 m 2 + 8 m + 64 = 0
m = − 8 ± 64 − 256 2 = − 8 ± − 192 2 = − 4 ± 4 3 i m = \frac{-8 \pm \sqrt{64 - 256}}{2} = \frac{-8 \pm \sqrt{-192}}{2} = -4 \pm 4\sqrt{3}i m = 2 − 8 ± 64 − 256
Step 3: Underdamped solution:
x ( t ) = e − 4 t ( C 1 cos 4 3 t + C 2 sin 4 3 t ) x(t) = e^{-4t}(C_1 \cos 4\sqrt{3}t + C_2 \sin 4\sqrt{3}t) x ( t ) = e − 4 t ( C 1 cos 4
Step 4: Apply initial conditions x ( 0 ) = 0.1 x(0) = 0.1 x ( 0 ) = 0.1 , x ′ ( 0 ) = 0 x'(0) = 0 x ′ ( 0 ) = 0 :
At t = 0 t = 0 t = 0 : 0.1 = C 1 0.1 = C_1 0.1 = C 1
x ′ ( t ) = − 4 e − 4 t ( C 1 cos 4 3 t + C 2 sin 4 3 t ) + e − 4 t ( − 4 3 C 1 sin 4 3 t + 4 3 C 2 cos 4 3 t ) x'(t) = -4e^{-4t}(C_1 \cos 4\sqrt{3}t + C_2 \sin 4\sqrt{3}t) + e^{-4t}(-4\sqrt{3}C_1 \sin 4\sqrt{3}t + 4\sqrt{3}C_2 \cos 4\sqrt{3}t) x ′ ( t ) = − 4 e
At t = 0 t = 0 t = 0 : 0 = − 4 C 1 + 4 3 C 2 0 = -4C_1 + 4\sqrt{3}C_2 0 = − 4 C 1 + 4 3
C 2 = C 1 3 = 0.1 3 = 3 30 C_2 = \frac{C_1}{\sqrt{3}} = \frac{0.1}{\sqrt{3}} = \frac{\sqrt{3}}{30} C 2 = 3
Final Answer:
x ( t ) = e − 4 t ( 0.1 cos 4 3 t + 3 30 sin 4 3 t ) m \boxed{x(t) = e^{-4t}\left(0.1 \cos 4\sqrt{3}t + \frac{\sqrt{3}}{30} \sin 4\sqrt{3}t\right) \text{ m}} x ( t ) = e − 4 t ( 0.1 cos 4
Q20. LCR Circuit
An LCR circuit has L = 0.5 L = 0.5 L = 0.5 H, R = 6 R = 6 R = 6 Ω, C = 0.02 C = 0.02 C = 0.02 F, and E ( t ) = 24 sin 10 t E(t) = 24\sin 10t E ( V. Find the charge given , .
Solution:
Step 1: Circuit equation:
L d 2 q d t 2 + R d q d t + q C = E ( t ) L\frac{d^2q}{dt^2} + R\frac{dq}{dt} + \frac{q}{C} = E(t) L d t 2 d 2 q
0.5 d 2 q d t 2 + 6 d q d t + 50 q = 24 sin 10 t 0.5\frac{d^2q}{dt^2} + 6\frac{dq}{dt} + 50q = 24\sin 10t 0.5 d t 2 d 2 q
d 2 q d t 2 + 12 d q d t + 100 q = 48 sin 10 t \frac{d^2q}{dt^2} + 12\frac{dq}{dt} + 100q = 48\sin 10t d t 2 d 2 q
Step 2: CF:
m 2 + 12 m + 100 = 0 m^2 + 12m + 100 = 0 m 2 + 12 m + 100 = 0
m = − 12 ± 144 − 400 2 = − 6 ± 8 i m = \frac{-12 \pm \sqrt{144 - 400}}{2} = -6 \pm 8i m = 2 − 12 ± 144 − 400
q c = e − 6 t ( C 1 cos 8 t + C 2 sin 8 t ) q_c = e^{-6t}(C_1 \cos 8t + C_2 \sin 8t) q c = e − 6 t ( C 1
Step 3: PI:
q p = 1 D 2 + 12 D + 100 48 sin 10 t q_p = \frac{1}{D^2 + 12D + 100} 48\sin 10t q p = D 2 + 12 D + 100
Substitute D 2 = − 100 D^2 = -100 D 2 = − 100 :
q p = 48 sin 10 t 12 D = 4 sin 10 t D = 4 ∫ sin 10 t d t = − 4 cos 10 t 10 = − 2 cos 10 t 5 q_p = \frac{48\sin 10t}{12D} = \frac{4\sin 10t}{D} = 4 \int \sin 10t \, dt = -\frac{4\cos 10t}{10} = -\frac{2\cos 10t}{5} q p = 12 D
Wait — check: D 2 + 100 = 0 D^2 + 100 = 0 D 2 + 100 = 0 at D 2 = − 100 D^2 = -100 D 2 = − 100 , so we need to use the derivative rule:
q p = 48 sin 10 t 12 D = 4 D sin 10 t = − 4 cos 10 t 10 = − 0.4 cos 10 t q_p = \frac{48\sin 10t}{12D} = \frac{4}{D}\sin 10t = -\frac{4\cos 10t}{10} = -0.4\cos 10t q p = 12 D 48
Actually, let's verify: 1 D 2 + 12 D + 100 sin 10 t \frac{1}{D^2 + 12D + 100}\sin 10t D 2 + 12 D + 100 1 sin 10 t with D 2 = − 100 D^2 = -100 :
= 1 12 D sin 10 t = 1 12 ⋅ ( − cos 10 t 10 ) = − cos 10 t 120 = \frac{1}{12D}\sin 10t = \frac{1}{12} \cdot \left(-\frac{\cos 10t}{10}\right) = -\frac{\cos 10t}{120} = 12 D 1 sin 10 t =
So q p = 48 × ( − cos 10 t 120 ) = − 48 cos 10 t 120 = − 2 cos 10 t 5 = − 0.4 cos 10 t q_p = 48 \times \left(-\frac{\cos 10t}{120}\right) = -\frac{48\cos 10t}{120} = -\frac{2\cos 10t}{5} = -0.4\cos 10t q p = 48 × ( −
Step 4: Apply ICs q ( 0 ) = 0 q(0) = 0 q ( 0 ) = 0 , q ′ ( 0 ) = 0 q'(0) = 0 q ′ ( 0 ) = 0 :
q ( t ) = e − 6 t ( C 1 cos 8 t + C 2 sin 8 t ) − 0.4 cos 10 t q(t) = e^{-6t}(C_1 \cos 8t + C_2 \sin 8t) - 0.4\cos 10t q ( t ) = e − 6 t ( C 1 cos 8 t
At t = 0 t = 0 t = 0 : 0 = C 1 − 0.4 ⇒ C 1 = 0.4 0 = C_1 - 0.4 \Rightarrow C_1 = 0.4 0 = C 1 − 0.4 ⇒ C
q ′ ( t ) = − 6 e − 6 t ( C 1 cos 8 t + C 2 sin 8 t ) + e − 6 t ( − 8 C 1 sin 8 t + 8 C 2 cos 8 t ) + 4 sin 10 t q'(t) = -6e^{-6t}(C_1 \cos 8t + C_2 \sin 8t) + e^{-6t}(-8C_1 \sin 8t + 8C_2 \cos 8t) + 4\sin 10t q ′ ( t ) = − 6 e
At t = 0 t = 0 t = 0 : 0 = − 6 C 1 + 8 C 2 ⇒ C 2 = 6 C 1 8 = 6 × 0.4 8 = 0.3 0 = -6C_1 + 8C_2 \Rightarrow C_2 = \frac{6C_1}{8} = \frac{6 \times 0.4}{8} = 0.3 0 = − 6 C 1 +
Final Answer:
q ( t ) = e − 6 t ( 0.4 cos 8 t + 0.3 sin 8 t ) − 0.4 cos 10 t \boxed{q(t) = e^{-6t}(0.4\cos 8t + 0.3\sin 8t) - 0.4\cos 10t} q ( t ) = e − 6 t ( 0.4 cos 8 t + 0.3
Q21. Simple Harmonic Motion — Spring-Mass System
A 4 4 4 kg mass stretches a spring by 0.392 0.392 0.392 m. The mass is set in motion from equilibrium with an upward velocity of 0.4 0.4 0.4 m/s. Find the position function.
Solution:
Step 1: Find spring constant k k k :
m g = k ⋅ Δ L ⇒ 4 × 9.8 = k × 0.392 mg = k \cdot \Delta L \Rightarrow 4 \times 9.8 = k \times 0.392 m g = k ⋅ Δ L ⇒ 4 × 9.8 = k × 0.392
k = 39.2 0.392 = 100 N/m k = \frac{39.2}{0.392} = 100 \text{ N/m} k = 0.392 39.2 = 100 N/m
Step 2: Equation of motion (no damping, no external force):
4 d 2 x d t 2 + 100 x = 0 4\frac{d^2x}{dt^2} + 100x = 0 4 d t 2 d 2 x + 100 x
d 2 x d t 2 + 25 x = 0 \frac{d^2x}{dt^2} + 25x = 0 d t 2 d 2 x + 25 x
Step 3: Solve:
m 2 + 25 = 0 ⇒ m = ± 5 i m^2 + 25 = 0 \Rightarrow m = \pm 5i m 2 + 25 = 0 ⇒ m = ± 5 i
x ( t ) = C 1 cos 5 t + C 2 sin 5 t x(t) = C_1 \cos 5t + C_2 \sin 5t x ( t ) = C 1 cos 5 t + C 2 sin
Step 4: Apply ICs: x ( 0 ) = 0 x(0) = 0 x ( 0 ) = 0 (equilibrium), x ′ ( 0 ) = − 0.4 x'(0) = -0.4 x ′ ( 0 ) = − 0.4 (upward = negative):
0 = C 1 0 = C_1 0 = C 1
x ′ ( t ) = − 5 C 1 sin 5 t + 5 C 2 cos 5 t x'(t) = -5C_1 \sin 5t + 5C_2 \cos 5t x ′ ( t ) = − 5 C 1 sin 5 t + 5
− 0.4 = 5 C 2 ⇒ C 2 = − 0.08 -0.4 = 5C_2 \Rightarrow C_2 = -0.08 − 0.4 = 5 C 2 ⇒ C 2 = − 0.08
Final Answer:
x ( t ) = − 0.08 sin 5 t m \boxed{x(t) = -0.08\sin 5t \text{ m}} x ( t ) = − 0.08 sin 5 t m
Q22. Forced Vibrations — Resonance
Solve: d 2 x d t 2 + 9 x = 6 cos 3 t \frac{d^2x}{dt^2} + 9x = 6\cos 3t d t 2 d 2 x + 9 x with , .
Solution:
Step 1: CF: m 2 + 9 = 0 ⇒ m = ± 3 i m^2 + 9 = 0 \Rightarrow m = \pm 3i m 2 + 9 = 0 ⇒ m = ± 3 i
x c = C 1 cos 3 t + C 2 sin 3 t x_c = C_1 \cos 3t + C_2 \sin 3t x c = C 1 cos 3 t + C
Step 2: PI (resonance — driving frequency = natural frequency):
x p = 1 D 2 + 9 6 cos 3 t = 6 ⋅ t 2 ⋅ 3 sin 3 t = t sin 3 t x_p = \frac{1}{D^2 + 9} 6\cos 3t = 6 \cdot \frac{t}{2 \cdot 3}\sin 3t = t\sin 3t x p = D 2 + 9
(Using 1 D 2 + a 2 cos a t = t 2 a sin a t \frac{1}{D^2+a^2}\cos at = \frac{t}{2a}\sin at D 2 + a 2 1 cos a )
Step 3: General solution:
x ( t ) = C 1 cos 3 t + C 2 sin 3 t + t sin 3 t x(t) = C_1 \cos 3t + C_2 \sin 3t + t\sin 3t x ( t ) = C 1 cos 3 t + C 2
Step 4: Apply ICs:
x ( 0 ) = C 1 = 0 x(0) = C_1 = 0 x ( 0 ) = C 1 = 0
x ′ ( t ) = − 3 C 1 sin 3 t + 3 C 2 cos 3 t + sin 3 t + 3 t cos 3 t x'(t) = -3C_1 \sin 3t + 3C_2 \cos 3t + \sin 3t + 3t\cos 3t x ′ ( t ) = − 3 C 1 sin 3 t
x ′ ( 0 ) = 3 C 2 = 0 ⇒ C 2 = 0 x'(0) = 3C_2 = 0 \Rightarrow C_2 = 0 x ′ ( 0 ) = 3 C 2 = 0 ⇒ C
Final Answer:
x ( t ) = t sin 3 t \boxed{x(t) = t\sin 3t} x ( t ) = t sin 3 t
> Note: The amplitude grows without bound — this is pure resonance .
Q23. Beats Phenomenon
Solve: d 2 x d t 2 + 25 x = 4 cos 4.5 t \frac{d^2x}{dt^2} + 25x = 4\cos 4.5t d t 2 d 2 x + 25 with , .
Solution:
Step 1: CF: m 2 + 25 = 0 ⇒ m = ± 5 i m^2 + 25 = 0 \Rightarrow m = \pm 5i m 2 + 25 = 0 ⇒ m = ± 5 i
x c = C 1 cos 5 t + C 2 sin 5 t x_c = C_1 \cos 5t + C_2 \sin 5t x c = C 1 cos 5 t + C
Step 2: PI:
x p = 1 D 2 + 25 4 cos 4.5 t = 4 cos 4.5 t − 20.25 + 25 = 4 cos 4.5 t 4.75 = 16 cos 4.5 t 19 x_p = \frac{1}{D^2 + 25} 4\cos 4.5t = \frac{4\cos 4.5t}{-20.25 + 25} = \frac{4\cos 4.5t}{4.75} = \frac{16\cos 4.5t}{19} x p = D
Step 3: General solution:
x ( t ) = C 1 cos 5 t + C 2 sin 5 t + 16 19 cos 4.5 t x(t) = C_1 \cos 5t + C_2 \sin 5t + \frac{16}{19}\cos 4.5t x ( t ) = C 1 cos 5 t + C 2
Step 4: Apply ICs:
x ( 0 ) = C 1 + 16 19 = 0 ⇒ C 1 = − 16 19 x(0) = C_1 + \frac{16}{19} = 0 \Rightarrow C_1 = -\frac{16}{19} x ( 0 ) = C 1 + 19 16
x ′ ( t ) = − 5 C 1 sin 5 t + 5 C 2 cos 5 t − 72 19 sin 4.5 t x'(t) = -5C_1 \sin 5t + 5C_2 \cos 5t - \frac{72}{19}\sin 4.5t x ′ ( t ) = − 5 C 1 sin 5 t +
x ′ ( 0 ) = 5 C 2 = 0 ⇒ C 2 = 0 x'(0) = 5C_2 = 0 \Rightarrow C_2 = 0 x ′ ( 0 ) = 5 C 2 = 0 ⇒ C
Final Answer:
x ( t ) = − 16 19 cos 5 t + 16 19 cos 4.5 t = 16 19 ( cos 4.5 t − cos 5 t ) \boxed{x(t) = -\frac{16}{19}\cos 5t + \frac{16}{19}\cos 4.5t = \frac{16}{19}(\cos 4.5t - \cos 5t)} x ( t ) = − 19 16
Using cos A − cos B = − 2 sin A + B 2 sin A − B 2 \cos A - \cos B = -2\sin\frac{A+B}{2}\sin\frac{A-B}{2} cos A − cos B = − 2 sin 2 A + B :
x ( t ) = 32 19 sin 4.75 t sin 0.25 t x(t) = \frac{32}{19}\sin 4.75t \sin 0.25t x ( t ) = 19 32 sin 4.75 t sin 0.25 t
> This shows beats — amplitude modulation at frequency 0.25 0.25 0.25 rad/s.
Q24. RLC Circuit — Transient and Steady-State
In an RLC circuit, L = 1 L = 1 L = 1 H, R = 2 R = 2 R = 2 Ω, C = 0.25 C = 0.25 C = 0.25 F, and E = 50 cos 2 t E = 50\cos 2t E = V. Find .
Solution:
Step 1: Equation:
d 2 q d t 2 + 2 d q d t + 4 q = 50 cos 2 t \frac{d^2q}{dt^2} + 2\frac{dq}{dt} + 4q = 50\cos 2t d t 2 d 2 q
Step 2: CF:
m 2 + 2 m + 4 = 0 ⇒ m = − 2 ± 4 − 16 2 = − 1 ± 3 i m^2 + 2m + 4 = 0 \Rightarrow m = \frac{-2 \pm \sqrt{4-16}}{2} = -1 \pm \sqrt{3}i m 2 + 2 m + 4 = 0 ⇒ m =
q c = e − t ( C 1 cos 3 t + C 2 sin 3 t ) q_c = e^{-t}(C_1 \cos \sqrt{3}t + C_2 \sin \sqrt{3}t) q c = e − t ( C 1
Step 3: PI:
q p = 1 D 2 + 2 D + 4 50 cos 2 t q_p = \frac{1}{D^2 + 2D + 4} 50\cos 2t q p = D 2 + 2 D + 4
Substitute D 2 = − 4 D^2 = -4 D 2 = − 4 :
q p = 50 cos 2 t 2 D = 25 ⋅ sin 2 t 2 = 25 2 sin 2 t q_p = \frac{50\cos 2t}{2D} = 25 \cdot \frac{\sin 2t}{2} = \frac{25}{2}\sin 2t q p = 2 D 50 c
Wait — check: 1 D cos 2 t = sin 2 t 2 \frac{1}{D}\cos 2t = \frac{\sin 2t}{2} D 1 cos 2 t = 2
Actually: q p = 50 2 D cos 2 t = 25 ⋅ sin 2 t 2 = 12.5 sin 2 t q_p = \frac{50}{2D}\cos 2t = 25 \cdot \frac{\sin 2t}{2} = 12.5\sin 2t q p = 2 D 50
Final Answer:
q ( t ) = e − t ( C 1 cos 3 t + C 2 sin 3 t ) + 25 2 sin 2 t \boxed{q(t) = e^{-t}(C_1 \cos \sqrt{3}t + C_2 \sin \sqrt{3}t) + \frac{25}{2}\sin 2t} q ( t ) = e − t ( C
The first term is transient (decays to zero), the second is steady-state .
Q25. Critical Damping Application
A door closer has m = 1 m = 1 m = 1 kg, k = 25 k = 25 k = 25 N/m. Find the damping coefficient c c c for critical damping, and solve if released from x = 0.1 x = 0.1 x = 0.1 m.
Solution:
Step 1: Critical damping condition:
c 2 = 4 m k = 4 × 1 × 25 = 100 c^2 = 4mk = 4 \times 1 \times 25 = 100 c 2 = 4 mk = 4 × 1 × 25 = 100
c = 10 N⋅s/m c = 10 \text{ N·s/m} c = 10 N⋅s/m
Step 2: Equation:
d 2 x d t 2 + 10 d x d t + 25 x = 0 \frac{d^2x}{dt^2} + 10\frac{dx}{dt} + 25x = 0 d t 2 d 2 x +
Step 3: Solve:
m 2 + 10 m + 25 = 0 ⇒ ( m + 5 ) 2 = 0 ⇒ m = − 5 , − 5 m^2 + 10m + 25 = 0 \Rightarrow (m+5)^2 = 0 \Rightarrow m = -5, -5 m 2 + 10 m + 25 = 0 ⇒ ( m + 5 )
x ( t ) = ( C 1 + C 2 t ) e − 5 t x(t) = (C_1 + C_2 t)e^{-5t} x ( t ) = ( C 1 + C 2 t ) e
Step 4: Apply ICs x ( 0 ) = 0.1 x(0) = 0.1 x ( 0 ) = 0.1 , x ′ ( 0 ) = 0 x'(0) = 0 x ′ ( 0 ) = 0 :
0.1 = C 1 0.1 = C_1 0.1 = C 1
x ′ ( t ) = C 2 e − 5 t − 5 ( C 1 + C 2 t ) e − 5 t x'(t) = C_2 e^{-5t} - 5(C_1 + C_2 t)e^{-5t} x ′ ( t ) = C 2 e
0 = C 2 − 5 C 1 ⇒ C 2 = 0.5 0 = C_2 - 5C_1 \Rightarrow C_2 = 0.5 0 = C 2 − 5 C 1 ⇒ C
Final Answer:
x ( t ) = ( 0.1 + 0.5 t ) e − 5 t \boxed{x(t) = (0.1 + 0.5t)e^{-5t}} x ( t ) = ( 0.1 + 0.5 t ) e − 5 t
Quick Reference: Method Selection Guide
Equation Type Method to Use a n y ( n ) + ⋯ + a 0 y = f ( x ) a_n y^{(n)} + \dots + a_0 y = f(x) a n y ( n ) + ⋯ + a (constant coeffs)
Exam Tips & Common Mistakes
Always find CF first — the form of PI depends on it.
Check for resonance — if f ( x ) f(x) f ( x ) matches CF terms, multiply by x x x .
For Cauchy-Euler , remember: x D = δ xD = \delta x D = δ , x 2 D 2 = δ ( δ − 1 ) x^2D^2 = \delta(\delta-1) x where and .
Summary
This unit covers the complete toolkit for solving higher-order ODEs:
Linear ODEs with constant coefficients via auxiliary equations
Simultaneous systems via elimination
Variable coefficient equations via transformation
Cauchy-Euler equations via logarithmic substitution
Variation of parameters for general non-homogeneous equations
Engineering applications in vibrations and circuits
Master these 25 solved problems and you'll be well-prepared for any exam on higher-order differential equations!