← 2022 Paper 2

UPSC 2022 Maths Optional Paper 2 Q7a — Step-by-Step Solution

15 marks · Section B

Second-order linear PDEs with constant coefficients (CF, PI) · PDEs · asked 12× in 13 yrs · Read the full method →

Question

Find general solution of (D2+DD6D2)z=x2sin(x+y)(D^2+DD'-6D'^2)z=x^2\sin(x+y) where D=/xD=\partial/\partial x, D=/yD'=\partial/\partial y.

Technique

Factor operator f(D,D)=(D+3D)(D2D)f(D,D')=(D+3D')(D-2D'); CF ϕ1(y+2x)+ϕ2(y3x)\phi_1(y+2x)+\phi_2(y-3x); PI for x2ei(x+y)x^2 e^{i(x+y)} via shift theorem; extract imaginary part.

Solution

Step 1 — Factor operator

D2+DD6D2=(D+3D)(D2D)D^2+DD'-6D'^2=(D+3D')(D-2D')? Check: (D+3D)(D2D)=D22DD+3DD6D2=D2+DD6D2(D+3D')(D-2D')=D^2-2DD'+3D'D-6D'^2=D^2+DD'-6D'^2 ✓.

Step 2 — Complementary function

CF from each factor:

zc=ϕ1(y+2x)+ϕ2(y3x)z_c=\phi_1(y+2x)+\phi_2(y-3x).

Step 3 — Particular integral

For RHS =x2sin(x+y)=x^2\sin(x+y), write sin(x+y)=ei(x+y)\sin(x+y)=\Im e^{i(x+y)}.

Particular for x2ei(x+y)x^2 e^{i(x+y)}: this has the form x2eax+byx^2\cdot e^{ax+by} with a=b=ia=b=i but multiplied by polynomial.

Standard trick: Use the operator inverse 1f(D,D)\dfrac{1}{f(D,D')} acting on polynomial times exponential.

1f(D,D)[x2ei(x+y)]=ei(x+y)1f(D+i,D+i)[x2]\dfrac{1}{f(D,D')}[x^2 e^{i(x+y)}]=e^{i(x+y)}\cdot\dfrac{1}{f(D+i,D'+i)}[x^2].

(Standard shift: f(D,D)[eax+byP(x,y)]=eax+byf(D+a,D+b)[P]f(D,D')[e^{ax+by}P(x,y)]=e^{ax+by}f(D+a,D'+b)[P].)

Compute f(D+i,D+i)f(D+i,D'+i): (D+i)2+(D+i)(D+i)6(D+i)2(D+i)^2+(D+i)(D'+i)-6(D'+i)^2 =D2+2iD1+DD+iD+iD16(D2+2iD1)=D^2+2iD-1+DD'+iD+iD'-1-6(D'^2+2iD'-1) =D2+DD6D2+2iD+iD+iD12iD+(11+6)=D^2+DD'-6D'^2+2iD+iD+iD'-12iD'+(-1-1+6) =D2+DD6D2+3iD11iD+4=D^2+DD'-6D'^2+3iD-11iD'+4.

Hmm careful: 1+(1)+(6)(1)=11+6=4-1+(-1)+(-6)(-1)=-1-1+6=4. ✓

2iD+iD12iD+iD2iD+iD-12iD'+iD': let me redo coefficient of DD and DD':

Coefficient of DD: 2i+i=3i2i+i=3i. Coefficient of DD': i12i=11ii-12i=-11i.

So f(D+i,D+i)=D2+DD6D2+3iD11iD+4f(D+i,D'+i)=D^2+DD'-6D'^2+3iD-11iD'+4.

Acting on x2x^2: only DD-terms produce nonzero results.

D[x2]=2xD[x^2]=2x, D2[x2]=2D^2[x^2]=2, D[x2]=0D'[x^2]=0, D2[x2]=0D'^2[x^2]=0, DD[x2]=0DD'[x^2]=0.

f(D+i,D+i)[x2]=2+00+3i(2x)11i(0)+4(x2)=2+6ix+4x2f(D+i,D'+i)[x^2]=2+0-0+3i(2x)-11i(0)+4(x^2)=2+6ix+4x^2.

We need 1/f(D+i,D+i)[x2]1/f(D+i,D'+i)[x^2], i.e., solve f(D+i,D+i)[P]=x2f(D+i,D'+i)[P]=x^2 for polynomial PP.

Let P=Ax2+Bx+CP=Ax^2+Bx+C (no yy since RHS is yy-free).

f(D+i,D+i)[P]=2A+3i(2Ax+B)11i0+4(Ax2+Bx+C)f(D+i,D'+i)[P]=2A+3i(2Ax+B)-11i\cdot 0+4(Ax^2+Bx+C) =2A+6iAx+3iB+4Ax2+4Bx+4C=2A+6iAx+3iB+4Ax^2+4Bx+4C =4Ax2+(6iA+4B)x+(2A+3iB+4C)=4Ax^2+(6iA+4B)x+(2A+3iB+4C).

Match to x2x^2:

So P=x243ix81332P=\dfrac{x^2}{4}-\dfrac{3ix}{8}-\dfrac{13}{32}.

Step 4 — Particular integral (complex form)

zp(C)=ei(x+y) ⁣[x243ix81332]z_p^{(\mathbb C)}=e^{i(x+y)}\!\left[\dfrac{x^2}{4}-\dfrac{3ix}{8}-\dfrac{13}{32}\right].

Step 5 — Take imaginary part

ei(x+y)=cos(x+y)+isin(x+y)e^{i(x+y)}=\cos(x+y)+i\sin(x+y).

zp(C)=(cos(x+y)+isin(x+y)) ⁣[x2413323ix8]z_p^{(\mathbb C)}=(\cos(x+y)+i\sin(x+y))\!\left[\dfrac{x^2}{4}-\dfrac{13}{32}-\dfrac{3ix}{8}\right].

Real part = cos(x+y)(x2/413/32)+sin(x+y)(3x/8)\cos(x+y)\cdot(x^2/4-13/32)+\sin(x+y)\cdot(3x/8).

Imaginary part = sin(x+y)(x2/413/32)cos(x+y)(3x/8)\sin(x+y)\cdot(x^2/4-13/32)-\cos(x+y)\cdot(3x/8).

We want zp(C)\Im z_p^{(\mathbb C)} (since RHS is sin\sin):

zp= ⁣(x241332)sin(x+y)3x8cos(x+y).z_p=\!\left(\dfrac{x^2}{4}-\dfrac{13}{32}\right)\sin(x+y)-\dfrac{3x}{8}\cos(x+y).

Step 6 — General solution

Answer

  z=ϕ1(y+2x)+ϕ2(y3x)+ ⁣(x241332)sin(x+y)3x8cos(x+y).  \boxed{\;z=\phi_1(y+2x)+\phi_2(y-3x)+\!\left(\dfrac{x^2}{4}-\dfrac{13}{32}\right)\sin(x+y)-\dfrac{3x}{8}\cos(x+y).\;}
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.