← 2015 Paper 1

UPSC 2015 Maths Optional Paper 1 Q8d — Step-by-Step Solution

13 marks · Section B

Euler-Cauchy equation · ODEs · asked 8× in 13 yrs · Read the full method →

Question

Solve x4y(4)+6x3y(3)+4x2y2xy4y=x2+2cos(logex).x^4 y^{(4)}+6x^3 y^{(3)}+4x^2 y''-2x y'-4y=x^2+2\cos(\log_e x).

Technique

Euler equation ⇒ substitute x=etx=e^t to get constant-coefficient operator f(D)=D43D24=(D2)(D+2)(D2+1)f(D)=D^4-3D^2-4=(D-2)(D+2)(D^2+1); CF from the roots; particular integrals handle resonance (e2te^{2t} and cost\cos t are both in the kernel) via tt-multiplication; convert back via t=logxt=\log x.

Solution

Strategy. The LHS is an Euler equation. Substitute x=etx=e^t (so t=logxt=\log x) and use xDx=DtxD_x=D_t, x2Dx2=Dt(Dt1)x^2 D_x^2=D_t(D_t-1), etc.

Step 1 — Convert to constant coefficients

Let D=d/dtD=d/dt. Standard formulae:

Substitute:

D(D1)(D2)(D3)y+6D(D1)(D2)y+4D(D1)y2Dy4y=e2t+2cost.D(D-1)(D-2)(D-3)y+6D(D-1)(D-2)y+4D(D-1)y-2Dy-4y=e^{2t}+2\cos t.

(RHS: x2=e2tx^2=e^{2t}, cos(logx)=cost\cos(\log x)=\cos t.)

Step 2 — Simplify the LHS polynomial in DD

Let me compute the operator coefficient by coefficient.

D(D1)(D2)(D3)D(D-1)(D-2)(D-3): expand step by step.

D(D1)=D2DD(D-1)=D^2-D.

(D2D)(D2)=D32D2D2+2D=D33D2+2D(D^2-D)(D-2)=D^3-2D^2-D^2+2D=D^3-3D^2+2D.

(D33D2+2D)(D3)=D43D33D3+9D2+2D26D=D46D3+11D26D(D^3-3D^2+2D)(D-3)=D^4-3D^3-3D^3+9D^2+2D^2-6D=D^4-6D^3+11D^2-6D.

6D(D1)(D2)=6(D33D2+2D)=6D318D2+12D6D(D-1)(D-2)=6(D^3-3D^2+2D)=6D^3-18D^2+12D.

4D(D1)=4D24D4D(D-1)=4D^2-4D.

2D=2D-2D=-2D.

4=4-4=-4.

Sum:

LHS operator: D43D24D^4-3D^2-4.

Step 3 — Factor the operator

D43D24=(D24)(D2+1)=(D2)(D+2)(D2+1)D^4-3D^2-4=(D^2-4)(D^2+1)=(D-2)(D+2)(D^2+1).

Step 4 — Auxiliary equation; complementary function

Roots: D=2,2,±iD=2,-2,\pm i.

Complementary function in tt:

yc=c1e2t+c2e2t+c3cost+c4sint.y_c=c_1 e^{2t}+c_2 e^{-2t}+c_3\cos t+c_4\sin t.

In xx (using t=logxt=\log x, et=xe^t=x):

yc=c1x2+c2x2+c3cos(logx)+c4sin(logx).y_c=c_1 x^2+\dfrac{c_2}{x^2}+c_3\cos(\log x)+c_4\sin(\log x).

Step 5 — Particular integral for e2te^{2t}

RHS contains e2te^{2t}. Standard formula: yp=e2tf(D)at D=2y_p=\dfrac{e^{2t}}{f(D)}\bigg|_{\text{at }D=2}, where f(D)=D43D24f(D)=D^4-3D^2-4. But f(2)=16124=0f(2)=16-12-4=0 — resonance.

Use f(D)=(D2)(D+2)(D2+1)f(D)=(D-2)(D+2)(D^2+1). With D=2D=2: f(2)=0f(2)=0, single root. Apply 1D2\dfrac{1}{D-2} specially:

yp=e2t(D2)(D+2)(D2+1)=1(D+2)(D2+1)D=2e2tD2shifted.y_p=\dfrac{e^{2t}}{(D-2)(D+2)(D^2+1)}=\dfrac{1}{(D+2)(D^2+1)}\bigg|_{D=2}\cdot\dfrac{e^{2t}}{D-2}\bigg|_{\text{shifted}}.

For a simple-pole resonance, yp=e2tf(2)ty_p=\dfrac{e^{2t}}{f'(2)}\cdot t — standard “shift” formula.

f(D)=D43D24f(D)=D^4-3D^2-4, f(D)=4D36D=2D(2D23)f'(D)=4D^3-6D=2D(2D^2-3). f(2)=22(83)=20f'(2)=2\cdot 2\cdot(8-3)=20.

So yp,1=te2t20y_{p,1}=\dfrac{te^{2t}}{20}.

Step 6 — Particular integral for 2cost2\cos t

RHS contains 2cost2\cos t. Operator f(D)f(D) evaluated at D=±iD=\pm i: f(i)=i43i24=1+34=0f(i)=i^4-3i^2-4=1+3-4=0 — resonance again!

Resonance with D2+1D^2+1 factor. Use the shift method.

For cost=eit\cos t=\Re e^{it}, ypy_p for eite^{it} part: f(i)=0f(i)=0, simple root of D2+1=(Di)(D+i)D^2+1=(D-i)(D+i).

f(D)=4D36Df'(D)=4D^3-6D. f(i)=4i36i=4i6i=10if'(i)=4i^3-6i=-4i-6i=-10i.

Particular for eite^{it}: yp,eit=teitf(i)=teit10i=iteit10y_{p,e^{it}}=\dfrac{te^{it}}{f'(i)}=\dfrac{te^{it}}{-10i}=\dfrac{it e^{it}}{10}.

For 2cost=2(eit)2\cos t=2\Re(e^{it}): yp,2=2(iteit10)=t5(ieit)=t5(icostsint)=t5(sint)=tsint5y_{p,2}=2\Re\bigl(\dfrac{it e^{it}}{10}\bigr)=\dfrac{t}{5}\Re(ie^{it})=\dfrac{t}{5}\Re(i\cos t-\sin t)=\dfrac{t}{5}(-\sin t)=-\dfrac{t\sin t}{5}.

Hmm wait: (icostsint)=(icost)(sint)=0sint=sint\Re(i\cos t-\sin t)=\Re(i\cos t)-\Re(\sin t)=0-\sin t=-\sin t. So yp,2=2t(sint)/10=tsint/5y_{p,2}=2\cdot t\cdot(-\sin t)/10=-t\sin t/5.

Actually let me redo more carefully. 2cost=(2eit)2\cos t=\Re(2e^{it}). Particular for 2eit2e^{it}: 2teitf(i)=2teit10i=2teiti10ii=2iteit10=iteit52\cdot \dfrac{te^{it}}{f'(i)}=\dfrac{2t e^{it}}{-10i}=\dfrac{2t e^{it}\cdot i}{-10i\cdot i}=\dfrac{2it e^{it}}{10}=\dfrac{it e^{it}}{5}.

Take real part: (iteit/5)=(t/5)(ieit)=(t/5)(icostsint)=(t/5)(sint)=tsint/5\Re(it e^{it}/5)=(t/5)\Re(ie^{it})=(t/5)\Re(i\cos t-\sin t)=(t/5)(-\sin t)=-t\sin t/5.

Hmm: but my standard sign-tracking suggests (cost+isint)i=icostsint\Re(\cos t+i\sin t)\cdot i=i\cos t-\sin t, so =sint\Re=-\sin t. Confirmed.

yp,2=tsint5y_{p,2}=-\dfrac{t\sin t}{5}.

Wait — let me double-check by direct substitution. Try y=Atsint+Btcosty=At\sin t+Bt\cos t, compute (D43D24)y(D^4-3D^2-4)y and match to 2cost2\cos t.

Actually easier: yp=tg(t)y_p=t\,g(t) where g(t)=Acost+Bsintg(t)=A\cos t+B\sin t. By the operator factoring f(D)=(D2)(D+2)(D2+1)f(D)=(D-2)(D+2)(D^2+1), the resonance is with the (D2+1)(D^2+1) factor. Reduce: let L1=D2+1L_1=D^2+1 and L2=(D2)(D+2)=D24L_2=(D-2)(D+2)=D^2-4. Then f(D)=L1L2=L2L1f(D)=L_1 L_2=L_2 L_1.

For RHS =2cost=2\cos t, L1(2cost)=(D2+1)(2cost)=2(cost+cost)=0L_1(2\cos t)=(D^2+1)(2\cos t)=2(-\cos t+\cos t)=0. So 2cost2\cos t is in the kernel of L1L_1.

Apply f1f^{-1}: yp=L21L11(2cost)y_p=L_2^{-1}L_1^{-1}(2\cos t). We have L11(2cost)=L_1^{-1}(2\cos t)= particular solution of L1z=2costL_1 z=2\cos t, i.e. z+z=2costz''+z=2\cos t. Standard: z=tsintz=t\sin t. (Check: z=sint+tcostz'=\sin t+t\cos t, z=cost+costtsint=2costtsintz''=\cos t+\cos t-t\sin t=2\cos t-t\sin t. z+z=2costtsint+tsint=2costz''+z=2\cos t-t\sin t+t\sin t=2\cos t ✓.)

Now apply L21=1/(D24)L_2^{-1}=1/(D^2-4): yp=tsintD24y_p=\dfrac{t\sin t}{D^2-4}. Plug D24D^2-4 acting on tsintt\sin t — find ww such that (D24)w=tsint(D^2-4)w=t\sin t.

Try w=(At+B)sint+(Ct+E)costw=(At+B)\sin t+(Ct+E)\cos t. Then w=Asint+(At+B)cost+Ccost(Ct+E)sintw'=A\sin t+(At+B)\cos t+C\cos t-(Ct+E)\sin t =(ACtE)sint+(At+B+C)cost=(A-Ct-E)\sin t+(At+B+C)\cos t.

w=Csint+(ACtE)cost+Acost(At+B+C)sintw''=-C\sin t+(A-Ct-E)\cos t+A\cos t-(At+B+C)\sin t =(CAtBC)sint+(2ACtE)cost=(-C-At-B-C)\sin t+(2A-Ct-E)\cos t =(AtB2C)sint+(2ACtE)cost=(-At-B-2C)\sin t+(2A-Ct-E)\cos t.

w4w=(AtB2C4At4B)sint+(2ACtE4Ct4E)costw''-4w=(-At-B-2C-4At-4B)\sin t+(2A-Ct-E-4Ct-4E)\cos t

Wait, 4w=4[(At+B)sint+(Ct+E)cost]=(4At4B)sint+(4Ct4E)cost-4w=-4[(At+B)\sin t+(Ct+E)\cos t]=(-4At-4B)\sin t+(-4Ct-4E)\cos t.

w4ww''-4w: coefficient of sint\sin t: AtB2C4At4B=5At5B2C-At-B-2C-4At-4B=-5At-5B-2C.

Coefficient of cost\cos t: 2ACtE4Ct4E=2A5Ct5E2A-Ct-E-4Ct-4E=2A-5Ct-5E.

Match to tsintt\sin t: coefficient of sint\sin t is 5At5B2C=t-5At-5B-2C=t, giving 5A=1-5A=1 and 5B2C=0-5B-2C=0.

Coefficient of cost\cos t is 2A5Ct5E=02A-5Ct-5E=0, giving 5C=0-5C=0 and 2A5E=02A-5E=0.

So A=1/5A=-1/5, C=0C=0, BB from 5B20=0B=0-5B-2\cdot 0=0\Rightarrow B=0. EE from 2(1/5)5E=02/5=5EE=2/252(-1/5)-5E=0\Rightarrow -2/5=5E\Rightarrow E=-2/25.

w=t5sint225costw=-\dfrac{t}{5}\sin t-\dfrac{2}{25}\cos t.

Check: derivative of 225cost-\dfrac{2}{25}\cos t is in L21L_2^{-1} kernel? No — it’s a particular contribution.

So yp,2=w=tsint52cost25y_{p,2}=w=-\dfrac{t\sin t}{5}-\dfrac{2\cos t}{25}.

The 2cost25-\dfrac{2\cos t}{25} term can be absorbed into ycy_c (since cost\cos t is in the kernel of L1L_1 and hence of ff). So drop it:

yp,2=tsint5.y_{p,2}=-\dfrac{t\sin t}{5}.

Good — matches my earlier answer.

Step 7 — Total general solution in tt

y(t)=c1e2t+c2e2t+c3cost+c4sint+te2t20tsint5.y(t)=c_1 e^{2t}+c_2 e^{-2t}+c_3\cos t+c_4\sin t+\dfrac{te^{2t}}{20}-\dfrac{t\sin t}{5}.

Step 8 — Back to xx

e2t=x2e^{2t}=x^2, e2t=1/x2e^{-2t}=1/x^2, cost=cos(logx)\cos t=\cos(\log x), sint=sin(logx)\sin t=\sin(\log x), te2t=x2logxte^{2t}=x^2\log x, tsint=sin(logx)logxt\sin t=\sin(\log x)\cdot\log x.

Answer

  y=c1x2+c2x2+c3cos(logx)+c4sin(logx)+x2logx20(logx)sin(logx)5.  \boxed{\;y=c_1 x^2+\dfrac{c_2}{x^2}+c_3\cos(\log x)+c_4\sin(\log x)+\dfrac{x^2\log x}{20}-\dfrac{(\log x)\sin(\log x)}{5}.\;}
We post more of this — worked solutions, CSAT trap breakdowns, guide chapters — a few times a week on Telegram. Free, no sign-in. Join

This solution is part of the Maths Coverage Map — 13 years, mapped. Get the take-away PDF free.