← 2022 Paper 1
UPSC 2022 Maths Optional Paper 1 Q8a-ii — Step-by-Step Solution 10 marks · Section B
Euler-Cauchy equation · ODEs · asked 8× in 13 yrs · Read the full method →
Question
Solve ( 3 x + 2 ) 2 d 2 y d x 2 + 5 ( 3 x + 2 ) d y d x − 3 y = x 2 + x + 1 (3x+2)^2\dfrac{d^2 y}{dx^2}+5(3x+2)\dfrac{dy}{dx}-3y=x^2+x+1 ( 3 x + 2 ) 2 d x 2 d 2 y + 5 ( 3 x + 2 ) d x d y − 3 y = x 2 + x + 1 .
Technique
Substitution u = 3 x + 2 u=3x+2 u = 3 x + 2 converts to Cauchy-Euler form; then u = e t u=e^t u = e t converts to constant-coefficient; CF from char roots 1 / 3 , − 1 1/3,-1 1/3 , − 1 ; PI for each term in RHS via 1 / f ( D ) 1/f(D) 1/ f ( D ) acting on e a t e^{at} e a t .
Solution
Setup. Cauchy-Euler-like equation with ( 3 x + 2 ) (3x+2) ( 3 x + 2 ) playing the role of x x x . Substitute u = 3 x + 2 u=3x+2 u = 3 x + 2 (so d u = 3 d x du=3\,dx d u = 3 d x , x = ( u − 2 ) / 3 x=(u-2)/3 x = ( u − 2 ) /3 ).
Step 1 — Convert derivatives
d y d x = d y d u ⋅ d u d x = 3 d y d u \dfrac{dy}{dx}=\dfrac{dy}{du}\cdot\dfrac{du}{dx}=3\dfrac{dy}{du} d x d y = d u d y ⋅ d x d u = 3 d u d y .
d 2 y d x 2 = d d x ( 3 d y d u ) = 3 ⋅ d 2 y d u 2 ⋅ 3 = 9 d 2 y d u 2 \dfrac{d^2 y}{dx^2}=\dfrac{d}{dx}\!\left(3\dfrac{dy}{du}\right)=3\cdot\dfrac{d^2 y}{du^2}\cdot 3=9\dfrac{d^2 y}{du^2} d x 2 d 2 y = d x d ( 3 d u d y ) = 3 ⋅ d u 2 d 2 y ⋅ 3 = 9 d u 2 d 2 y .
Substitute:
u 2 ⋅ 9 y u ′ ′ + 5 u ⋅ 3 y u ′ − 3 y = x 2 + x + 1 u^2\cdot 9 y''_u+5u\cdot 3y'_u-3y=x^2+x+1 u 2 ⋅ 9 y u ′′ + 5 u ⋅ 3 y u ′ − 3 y = x 2 + x + 1 ,
9 u 2 y u ′ ′ + 15 u y u ′ − 3 y = x 2 + x + 1 9u^2 y''_u+15u y'_u-3y=x^2+x+1 9 u 2 y u ′′ + 15 u y u ′ − 3 y = x 2 + x + 1 .
Divide by 3:
3 u 2 y u ′ ′ + 5 u y u ′ − y = x 2 + x + 1 3 3u^2 y''_u+5u y'_u-y=\dfrac{x^2+x+1}{3} 3 u 2 y u ′′ + 5 u y u ′ − y = 3 x 2 + x + 1 .
Express RHS in terms of u u u : x = ( u − 2 ) / 3 x=(u-2)/3 x = ( u − 2 ) /3 , so x 2 = ( u − 2 ) 2 / 9 x^2=(u-2)^2/9 x 2 = ( u − 2 ) 2 /9 , x 2 + x + 1 = ( u − 2 ) 2 / 9 + ( u − 2 ) / 3 + 1 x^2+x+1=(u-2)^2/9+(u-2)/3+1 x 2 + x + 1 = ( u − 2 ) 2 /9 + ( u − 2 ) /3 + 1 .
( u − 2 ) 2 / 9 = ( u 2 − 4 u + 4 ) / 9 (u-2)^2/9=(u^2-4u+4)/9 ( u − 2 ) 2 /9 = ( u 2 − 4 u + 4 ) /9 .
( u − 2 ) / 3 = ( 3 u − 6 ) / 9 (u-2)/3=(3u-6)/9 ( u − 2 ) /3 = ( 3 u − 6 ) /9 .
1 = 9 / 9 1=9/9 1 = 9/9 .
Sum: ( u 2 − 4 u + 4 + 3 u − 6 + 9 ) / 9 = ( u 2 − u + 7 ) / 9 (u^2-4u+4+3u-6+9)/9=(u^2-u+7)/9 ( u 2 − 4 u + 4 + 3 u − 6 + 9 ) /9 = ( u 2 − u + 7 ) /9 .
RHS becomes ( u 2 − u + 7 ) / ( 9 ⋅ 3 ) = ( u 2 − u + 7 ) / 27 (u^2-u+7)/(9\cdot 3)=(u^2-u+7)/27 ( u 2 − u + 7 ) / ( 9 ⋅ 3 ) = ( u 2 − u + 7 ) /27 .
So
3 u 2 y ′ ′ + 5 u y ′ − y = u 2 − u + 7 27 . ( ⋆ ) 3u^2 y''+5u y'-y=\dfrac{u^2-u+7}{27}.\qquad(\star) 3 u 2 y ′′ + 5 u y ′ − y = 27 u 2 − u + 7 . ( ⋆ )
Step 2 — Convert to constant-coefficient ODE via u = e t u=e^t u = e t
Let u = e t u=e^t u = e t , t = ln u t=\ln u t = ln u . Standard Euler-Cauchy:
u y ′ = D y u y'=Dy u y ′ = D y where D = d / d t D=d/dt D = d / d t .
u 2 y ′ ′ = D ( D − 1 ) y u^2 y''=D(D-1)y u 2 y ′′ = D ( D − 1 ) y .
Substitute:
3 D ( D − 1 ) y + 5 D y − y = e 2 t − e t + 7 27 3D(D-1)y+5Dy-y=\dfrac{e^{2t}-e^t+7}{27} 3 D ( D − 1 ) y + 5 D y − y = 27 e 2 t − e t + 7 ,
[ 3 D 2 − 3 D + 5 D − 1 ] y = e 2 t − e t + 7 27 [3D^2-3D+5D-1]y=\dfrac{e^{2t}-e^t+7}{27} [ 3 D 2 − 3 D + 5 D − 1 ] y = 27 e 2 t − e t + 7 ,
[ 3 D 2 + 2 D − 1 ] y = e 2 t − e t + 7 27 [3D^2+2D-1]y=\dfrac{e^{2t}-e^t+7}{27} [ 3 D 2 + 2 D − 1 ] y = 27 e 2 t − e t + 7 .
Factor: 3 D 2 + 2 D − 1 = ( 3 D − 1 ) ( D + 1 ) 3D^2+2D-1=(3D-1)(D+1) 3 D 2 + 2 D − 1 = ( 3 D − 1 ) ( D + 1 ) .
Step 3 — Complementary function
Roots of 3 D 2 + 2 D − 1 = 0 3D^2+2D-1=0 3 D 2 + 2 D − 1 = 0 : D = ( − 2 ± 4 + 12 ) / 6 = ( − 2 ± 4 ) / 6 = 1 / 3 , − 1 D=(-2\pm\sqrt{4+12})/6=(-2\pm 4)/6=1/3,-1 D = ( − 2 ± 4 + 12 ) /6 = ( − 2 ± 4 ) /6 = 1/3 , − 1 .
So CF = C 1 e t / 3 + C 2 e − t = C 1 u 1 / 3 + C 2 / u = C 1 ( 3 x + 2 ) 1 / 3 + C 2 / ( 3 x + 2 ) =C_1 e^{t/3}+C_2 e^{-t}=C_1 u^{1/3}+C_2/u=C_1(3x+2)^{1/3}+C_2/(3x+2) = C 1 e t /3 + C 2 e − t = C 1 u 1/3 + C 2 / u = C 1 ( 3 x + 2 ) 1/3 + C 2 / ( 3 x + 2 ) .
Step 4 — Particular integral
RHS = 1 27 ( e 2 t − e t + 7 ) \dfrac{1}{27}(e^{2t}-e^t+7) 27 1 ( e 2 t − e t + 7 ) .
PI for e 2 t e^{2t} e 2 t : f ( D ) = 3 D 2 + 2 D − 1 f(D)=3D^2+2D-1 f ( D ) = 3 D 2 + 2 D − 1 , f ( 2 ) = 12 + 4 − 1 = 15 f(2)=12+4-1=15 f ( 2 ) = 12 + 4 − 1 = 15 . PI = e 2 t 15 =\dfrac{e^{2t}}{15} = 15 e 2 t .
PI for − e t -e^t − e t : f ( 1 ) = 3 + 2 − 1 = 4 f(1)=3+2-1=4 f ( 1 ) = 3 + 2 − 1 = 4 . PI = − e t 4 =-\dfrac{e^t}{4} = − 4 e t .
PI for 7 7 7 (constant): f ( 0 ) = − 1 f(0)=-1 f ( 0 ) = − 1 . PI = 7 − 1 = − 7 =\dfrac{7}{-1}=-7 = − 1 7 = − 7 .
Combine:
y p = 1 27 [ e 2 t 15 − e t 4 − 7 ] = e 2 t 405 − e t 108 − 7 27 . y_p=\dfrac{1}{27}\!\left[\dfrac{e^{2t}}{15}-\dfrac{e^t}{4}-7\right]=\dfrac{e^{2t}}{405}-\dfrac{e^t}{108}-\dfrac{7}{27}. y p = 27 1 [ 15 e 2 t − 4 e t − 7 ] = 405 e 2 t − 108 e t − 27 7 .
Convert back: e 2 t = u 2 = ( 3 x + 2 ) 2 e^{2t}=u^2=(3x+2)^2 e 2 t = u 2 = ( 3 x + 2 ) 2 , e t = u = 3 x + 2 e^t=u=3x+2 e t = u = 3 x + 2 .
y p = ( 3 x + 2 ) 2 405 − ( 3 x + 2 ) 108 − 7 27 . y_p=\dfrac{(3x+2)^2}{405}-\dfrac{(3x+2)}{108}-\dfrac{7}{27}. y p = 405 ( 3 x + 2 ) 2 − 108 ( 3 x + 2 ) − 27 7 .
Step 5 — General solution
Answer
y = C 1 ( 3 x + 2 ) 1 / 3 + C 2 3 x + 2 + ( 3 x + 2 ) 2 405 − 3 x + 2 108 − 7 27 . \boxed{\;y=C_1(3x+2)^{1/3}+\dfrac{C_2}{3x+2}+\dfrac{(3x+2)^2}{405}-\dfrac{3x+2}{108}-\dfrac{7}{27}.\;} y = C 1 ( 3 x + 2 ) 1/3 + 3 x + 2 C 2 + 405 ( 3 x + 2 ) 2 − 108 3 x + 2 − 27 7 .