The math optional, made finite. Daily Practice

Second-order linear PDEs with constant coefficients (CF, PI)

At a Glance

Why This Chapter Matters

This atom appears in every year except one since 2013, always worth 10 or 15 marks. The method — factor the operator, write the CF, compute a PI term by term — is completely algorithmic once you know the four PI rules. The only genuine difficulty is detecting resonance (where the naive exponential rule gives 0/00/0) and applying the xx-multiplier correction. Mastering that single trap, plus the shift identity for the xnsin/cosx^n\sin/\cos RHS, covers every hard case in the corpus.

Minimum Theory

Operator notation and factorisation. Write the PDE as F(D,D)z=R(x,y)F(D,D')z=R(x,y) where D=/xD=\partial/\partial x, D=/yD'=\partial/\partial y. The characteristic auxiliary equation sets D=mD=m and D=1D'=1 (or D/D=mD/D'=m): solving F(m,1)=0F(m,1)=0 for mm gives the slopes of the characteristic families. Each simple linear factor (DmD)(D-mD') contributes an arbitrary function ϕ(y+mx)\phi(y+mx) to the complementary function (CF). A repeated factor (DmD)k(D-mD')^k contributes ϕ1(y+mx)+xϕ2(y+mx)++xk1ϕk(y+mx)\phi_1(y+mx)+x\phi_2(y+mx)+\cdots+x^{k-1}\phi_k(y+mx).

PI rules (four cases).

RHS RRRuleCondition
eax+bye^{ax+by}eax+byF(a,b)\dfrac{e^{ax+by}}{F(a,b)}F(a,b)0F(a,b)\neq0
eax+bye^{ax+by}, resonantxeax+byFD(a,b)\dfrac{x\,e^{ax+by}}{F_D'(a,b)} (or yy\cdot)F(a,b)=0F(a,b)=0 for one simple factor
sin(ax+by)\sin(ax+by) or cos\cosReplace D2a2D^2\to-a^2, D2b2D'^2\to-b^2, DDabDD'\to-ab; divideNon-zero denominator
eax+byQ(x,y)e^{ax+by}Q(x,y)Shift: F(D,D)[eax+byΦ]=eax+byF(D+a,D+b)ΦF(D,D')[e^{ax+by}\Phi]=e^{ax+by}F(D+a,D'+b)\Phi; find Φ\PhiGeneral — used for xnsinx^n\sin

For a polynomial RR: expand 1/F(D,D)1/F(D,D') as a formal power series in D/DD/D' (or D/DD'/D) and apply to RR; the series terminates.

Boundary conditions. When boundary curves are given, pin down the arbitrary functions ϕi\phi_i by substituting the stated conditions into the general solution.

The CF–PI method: factor the operator, read off the CF, then apply the PI rules for each RHS term

Question Archetypes

Three patterns cover the corpus.

ArchetypeYou are seeing this when…
operator-pde”solve F(D,D)z=R(x,y)F(D,D')z=R(x,y)” — standard CF + PI
operator-pde-boundarysolve and fit boundary conditions (surface through given lines)
verify-pde-solution”verify that u=f(ξ)+g(η)u=f(\xi)+g(\eta) satisfies the PDE” — chain-rule differentiation

operator-pde (10 question(s); 2013, 2014, 2015, 2016, 2017, 2018, 2020, 2021, 2022, 2025)

Recognition Cues

Solution Template

  1. Factor. Write F(D,D)=(DmiD)F(D,D')=\prod(D-m_iD'). Find mim_i by solving F(m,1)=0F(m,1)=0 (for second-order, a quadratic; for third-order, a cubic — group and factor).
  2. CF. zc=ϕi(y+mix)z_c=\sum\phi_i(y+m_ix). For a repeated root mm: add xϕj(y+mx)x\phi_j(y+mx) for each repetition.
  3. PI. For each RHS term, apply the appropriate rule. Check for resonance before applying the exponential rule.
  4. Sum. z=zc+zpz=z_c+z_p.

Worked Example(s)

2014 Paper 2, 2014-P2-Q5a (10 marks)

Solve (2D25DD+2D2)z=24(yx)(2D^2-5DD'+2D'^2)z=24(y-x).

Factor. (2DD)(D2D)(2D-D')(D-2D'). Roots: m=1/2m=1/2 (from 2DD=02D-D'=0) and m=2m=2 (from D2D=0D-2D'=0).

CF. zc=ϕ1(y+12x)+ϕ2(y+2x)z_c=\phi_1(y+\tfrac12x)+\phi_2(y+2x), equivalently f(2y+x)+g(y+2x)f(2y+x)+g(y+2x).

PI. RHS =24(yx)=24(y-x), degree 1. Try cubic ansatz zp=Ax3+Bx2yz_p=Ax^3+Bx^2y; matching coefficients: (2D25DD+2D2)(Ax3+Bx2y)=12Ax+2By10Bx(2D^2-5DD'+2D'^2)(Ax^3+Bx^2y)=12Ax+2By-10Bx. Set =(24y24x)=(24y-24x): B=12B=12, A=3A=3.

  z=f(2y+x)+g(y+2x)+3x3+6x2y.  \boxed{\;z=f(2y+x)+g(y+2x)+3x^3+6x^2y.\;}


2015 Paper 2, 2015-P2-Q5b (10 marks)

Solve (D2+DD2D2)u=ex+y(D^2+DD'-2D'^2)u=e^{x+y}.

Factor. (DD)(D+2D)(D-D')(D+2D'). Roots m=1m=1 and m=2m=-2.

CF. uc=ϕ1(y+x)+ϕ2(y2x)u_c=\phi_1(y+x)+\phi_2(y-2x).

PI. F(1,1)=1+12=0F(1,1)=1+1-2=0: resonance with the factor (DD)(D-D').

Apply (D+2D)1(D+2D')^{-1} first (non-resonant, value 1+2=31+2=3): gives ex+y/3e^{x+y}/3.
Then (DD)1(ex+y/3)(D-D')^{-1}(e^{x+y}/3): solve (DD)u=ex+y(D-D')u=e^{x+y} via characteristics — u=xex+yu=xe^{x+y} (integration along x+y=constx+y=\text{const}). Divide by 3:

up=x3ex+y.u_p=\frac{x}{3}e^{x+y}.

  u=ϕ1(y+x)+ϕ2(y2x)+x3ex+y.  \boxed{\;u=\phi_1(y+x)+\phi_2(y-2x)+\frac{x}{3}e^{x+y}.\;}


2017 Paper 2, 2017-P2-Q5a (10 marks)

Solve (D22DD+D2)z=ex+2y+x3+sin2x(D^2-2DD'+D'^2)z=e^{x+2y}+x^3+\sin2x.

Factor. (DD)2(D-D')^2: repeated root m=1m=1.

CF. zc=ϕ1(y+x)+xϕ2(y+x)z_c=\phi_1(y+x)+x\phi_2(y+x).

PI for ex+2ye^{x+2y}. F(1,2)=(12)2=10F(1,2)=(1-2)^2=1\ne0: PI =ex+2y/1=ex+2y=e^{x+2y}/1=e^{x+2y}.

PI for x3x^3. Only D,D2D,D^2 act (since Dx3=0D'x^3=0): 1/D2x3=x5/201/D^2\cdot x^3=x^5/20.

PI for sin2x\sin2x. D24D^2\to-4, DD0DD'\to0, D20D'^2\to0: F4F\to-4, PI =sin2x/4=-\sin2x/4.

  z=ϕ1(y+x)+xϕ2(y+x)+ex+2y+x520sin2x4.  \boxed{\;z=\phi_1(y+x)+x\phi_2(y+x)+e^{x+2y}+\frac{x^5}{20}-\frac{\sin2x}{4}.\;}


2016 Paper 2, 2016-P2-Q7a (15 marks)

Solve (D32D2DDD2+2D3)z=ex+y(D^3-2D^2D'-DD'^2+2D'^3)z=e^{x+y}.

Factor. Auxiliary m32m2m+2=(m2)(m1)(m+1)m^3-2m^2-m+2=(m-2)(m-1)(m+1): roots m=2,1,1m=2,1,-1.

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

PI. F(1,1)=(12)(11)(1+1)=0F(1,1)=(1-2)(1-1)(1+1)=0: resonance on (DD)(D-D'). Non-vanishing factors at (1,1)(1,1): (12)(1+1)=2(1-2)(1+1)=-2.

zp=xex+y2=x2ex+y.z_p=\frac{xe^{x+y}}{-2}=-\frac{x}{2}e^{x+y}.

  z=ϕ1(y+2x)+ϕ2(y+x)+ϕ3(yx)x2ex+y.  \boxed{\;z=\phi_1(y+2x)+\phi_2(y+x)+\phi_3(y-x)-\frac{x}{2}e^{x+y}.\;}


2018 Paper 2, 2018-P2-Q7a (15 marks)

Solve (2D25DD+2D2)z=5sin(2x+y)+24(yx)+e3x+4y(2D^2-5DD'+2D'^2)z=5\sin(2x+y)+24(y-x)+e^{3x+4y}.

Factor. (D2D)(2DD)(D-2D')(2D-D'); CF =ϕ1(y+2x)+ϕ2(2y+x)=\phi_1(y+2x)+\phi_2(2y+x).

PI for 5sin(2x+y)5\sin(2x+y). D24D^2\to-4, DD2DD'\to-2, D21D'^2\to-1: F=2(4)5(2)+2(1)=0F=2(-4)-5(-2)+2(-1)=0. Resonance! The factor (D2D)(D-2D') annihilates sin(2x+y)\sin(2x+y) (since 2x+y2x+y matches the slope m=2m=2). Trial z=x(Pcos(2x+y)+Qsin(2x+y))z=x(P\cos(2x+y)+Q\sin(2x+y)); matching gives PI =53xcos(2x+y)=-\tfrac53x\cos(2x+y).

PI for 24(yx)24(y-x). Polynomial ansatz (cubic): PI =127x2y127xy2=127xy(xy)=\tfrac{12}{7}x^2y-\tfrac{12}{7}xy^2=\tfrac{12}{7}xy(x-y).

PI for e3x+4ye^{3x+4y}. F(3,4)=100F(3,4)=-10\ne0: PI =110e3x+4y=-\tfrac{1}{10}e^{3x+4y}.

  z=ϕ1(y+2x)+ϕ2(2y+x)53xcos(2x+y)+127xy(xy)110e3x+4y.  \boxed{\;z=\phi_1(y+2x)+\phi_2(2y+x)-\frac{5}{3}x\cos(2x+y)+\frac{12}{7}xy(x-y)-\frac{1}{10}e^{3x+4y}.\;}


2013/2022 Paper 2, 2013-P2-Q6a and 2022-P2-Q7a (15 marks each)

Solve (D2+DD6D2)z=x2sin(x+y)(D^2+DD'-6D'^2)z=x^2\sin(x+y).

Sources:,

Factor. (D+3D)(D2D)(D+3D')(D-2D'); CF =f(y3x)+g(y+2x)=f(y-3x)+g(y+2x).

PI. RHS is x2sin(x+y)x^2\cdot\sin(x+y) — use the shift identity. Replace sin(x+y)=Im(ei(x+y))\sin(x+y)=\operatorname{Im}(e^{i(x+y)}) and find PI of x2ei(x+y)x^2e^{i(x+y)}. The shifted operator F(D+i,D+i)F(D+i,D'+i) has non-zero constant term F(i,i)=i2+i26i2=40F(i,i)=i^2+i^2-6i^2=4\ne0, so no resonance. Expand (4+L)1(4+L)^{-1} with L=3iD11iD+L=3iD-11iD'+\ldots:

Φ=[4+L]1x2=14[x26ix+24+1816]=x243ix81332\Phi=[4+L]^{-1}x^2=\tfrac{1}{4}\bigl[x^2-\tfrac{6ix+2}{4}+\tfrac{-18}{16}\bigr]=\tfrac{x^2}{4}-\tfrac{3ix}{8}-\tfrac{13}{32}.

Take imaginary part after multiplying by ei(x+y)e^{i(x+y)}:

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

  z=f(y3x)+g(y+2x)+ ⁣(x241332) ⁣sin(x+y)3x8cos(x+y).  \boxed{\;z=f(y-3x)+g(y+2x)+\!\left(\frac{x^2}{4}-\frac{13}{32}\right)\!\sin(x+y)-\frac{3x}{8}\cos(x+y).\;}


2020 Paper 2, 2020-P2-Q5d (10 marks)

Solve (D32D2DDD2+2D3)z=e2x+y+sin(x2y)(D^3-2D^2D'-DD'^2+2D'^3)z=e^{2x+y}+\sin(x-2y).

Factor. (DD)(D+D)(D2D)(D-D')(D+D')(D-2D'); CF =ϕ1(y+x)+ϕ2(yx)+ϕ3(y+2x)=\phi_1(y+x)+\phi_2(y-x)+\phi_3(y+2x).

PI for e2x+ye^{2x+y}. F(2,1)=0F(2,1)=0 (since D2DD-2D' vanishes at m=2m=2, b/a=1/2b/a=1/2). Non-vanishing factors: (21)(2+1)=3(2-1)(2+1)=3. PI =xe2x+y/3=xe^{2x+y}/3.

PI for sin(x2y)\sin(x-2y). Replace D21D^2\to-1, D24D'^2\to-4, DD2DD'\to2 (since a=1,b=2a=1,b=-2, ab=2ab=-2). Cubic operator reduces to 3D6D3D-6D': 13D6Dsin(x2y)=115cos(x2y)\dfrac{1}{3D-6D'}\sin(x-2y)=-\dfrac{1}{15}\cos(x-2y).

  z=ϕ1(y+x)+ϕ2(yx)+ϕ3(y+2x)+x3e2x+ycos(x2y)15.  \boxed{\;z=\phi_1(y+x)+\phi_2(y-x)+\phi_3(y+2x)+\frac{x}{3}e^{2x+y}-\frac{\cos(x-2y)}{15}.\;}


2025 Paper 2, 2025-P2-Q5a (10 marks)

Solve (D2+DD2D2)z=ysinx(D^2+DD'-2D'^2)z=y\sin x.

Factor. (DD)(D+2D)(D-D')(D+2D'); CF =f(y+x)+g(y2x)=f(y+x)+g(y-2x).

PI. RHS =ysinx=y\sin x. Trial: zp=y(asinx+bcosx)+(csinx+dcosx)z_p=y(a\sin x+b\cos x)+(c\sin x+d\cos x). Apply the operator: only D2D^2 and DDDD' act non-trivially. Matching coefficients of ysinxy\sin x, ycosxy\cos x, sinx\sin x, cosx\cos x:

a=1a=-1, b=0b=0, c=0c=0, d=1d=-1.

  z=f(y+x)+g(y2x)ysinxcosx.  \boxed{\;z=f(y+x)+g(y-2x)-y\sin x-\cos x.\;}


2021 Paper 2, 2021-P2-Q7a (15 marks)

Solve (D2D23D+3D)z=xy+ex+2y(D^2-D'^2-3D+3D')z=xy+e^{x+2y}.

Factor. D2D23D+3D=(DD)(D+D)3(DD)=(DD)(D+D3)D^2-D'^2-3D+3D'=(D-D')(D+D')-3(D-D')=(D-D')(D+D'-3).

CF. From (DD)(D-D'): ϕ1(y+x)\phi_1(y+x). From (D+D3)(D+D'-3): e3xϕ2(yx)e^{3x}\phi_2(y-x) (verify: D+D3D+D'-3 kills e3xf(yx)e^{3x}f(y-x) since D+DD+D' kills functions of yxy-x and 3e3x3e^{3x} cancels).

PI for ex+2ye^{x+2y}. F(1,2)=(12)(1+23)=(1)(0)=0F(1,2)=(1-2)(1+2-3)=(-1)(0)=0: resonance on (D+D3)(D+D'-3). Apply (DD)1(D-D')^{-1} first (value 1-1): gives ex+2y-e^{x+2y}. Then (D+D3)1ex+2y=xex+2y(D+D'-3)^{-1}e^{x+2y}=xe^{x+2y}. So PI =xex+2y=-xe^{x+2y}.

PI for xyxy. Use (u,v)=(x+y,xy)(u,v)=(x+y,x-y) substitution to make the operator separable; result is a degree-3 polynomial in (x,y)(x,y) (full form in the source).

  z=ϕ1(y+x)+e3xϕ2(yx)+(poly PI)xex+2y.  \boxed{\;z=\phi_1(y+x)+e^{3x}\phi_2(y-x)+(\text{poly PI})-xe^{x+2y}.\;}

Common Traps


operator-pde-boundary (1 question(s); 2023)

Recognition Cues

Solution Template

  1. Solve the homogeneous PDE to get the general solution z=f(ξ1)+xg(ξ1)z=f(\xi_1)+xg(\xi_1) (or sum of fi(ξi)f_i(\xi_i)).
  2. Apply first BC (simpler curve, often x=0x=0 or z=0z=0): determines one arbitrary function.
  3. Apply second BC (substitute into the remaining general solution): get an equation for the second arbitrary function; solve by substitution.
  4. Write the explicit surface.

Worked Example(s)

2023 Paper 2, 2023-P2-Q6a (15 marks)

Find the surface through z=x=0z=x=0 and z1=xy=0z-1=x-y=0 satisfying zxx4zxy+4zyy=0z_{xx}-4z_{xy}+4z_{yy}=0.

Factor. (D2D)2(D-2D')^2: repeated root m=2m=2. General solution: z=f(y+2x)+xg(y+2x)z=f(y+2x)+xg(y+2x).

BC 1 (x=0x=0, z=0z=0): f(y)=0f(y)=0 for all yy, so f0f\equiv0. Now z=xg(y+2x)z=xg(y+2x).

BC 2 (x=yx=y, z=1z=1): xg(3x)=1x\,g(3x)=1, so g(3x)=1/xg(3x)=1/x. Set u=3xu=3x: g(u)=3/ug(u)=3/u.

z=x3y+2x=  3x2x+y.  z=x\cdot\frac{3}{y+2x}=\boxed{\;\frac{3x}{2x+y}.\;}

Verify: zxx4zxy+4zyy=0z_{xx}-4z_{xy}+4z_{yy}=0 ✓; zx=0=0z|_{x=0}=0 ✓; zy=x=3x/3x=1z|_{y=x}=3x/3x=1 ✓.

Common Traps


verify-pde-solution (1 question(s); 2024)

Recognition Cues

Solution Template

  1. Compute uxxu_{xx}, uyyu_{yy}, uttu_{tt} (or whatever derivatives the PDE involves) using ξx,ξy,ξt,ηx,ηy,ηt\xi_x,\xi_y,\xi_t,\eta_x,\eta_y,\eta_t.
  2. Substitute into the PDE; collect terms in f(ξ)+g(η)f''(\xi)+g''(\eta).
  3. The PDE holds iff the coefficient of f+gf''+g'' satisfies the stated algebraic condition.

Worked Example(s)

2024 Paper 2, 2024-P2-Q5a (10 marks)

Show u=f(xkt+iαy)+g(xktiαy)u=f(x-kt+i\alpha y)+g(x-kt-i\alpha y) satisfies uxx+uyy=utt/c2u_{xx}+u_{yy}=u_{tt}/c^2 when α2=1k2/c2\alpha^2=1-k^2/c^2.

With ξ=xkt+iαy\xi=x-kt+i\alpha y and η=xktiαy\eta=x-kt-i\alpha y:

uxx=f+gu_{xx}=f''+g''; uyy=(iα)2f+(iα)2g=α2(f+g)\quad u_{yy}=(i\alpha)^2f''+(-i\alpha)^2g''=-\alpha^2(f''+g''); utt=k2(f+g)\quad u_{tt}=k^2(f''+g'').

uxx+uyy=(1α2)(f+g),uttc2=k2c2(f+g).u_{xx}+u_{yy}=(1-\alpha^2)(f''+g''),\qquad \frac{u_{tt}}{c^2}=\frac{k^2}{c^2}(f''+g'').

Equating: 1α2=k2/c21-\alpha^2=k^2/c^2, i.e. α2=1k2/c2\alpha^2=1-k^2/c^2. \blacksquare

Common Traps


Marks-Aware Writing

10-mark questions (compulsory Q5 slot): Show factorisation (one line), write CF (one line), compute PI with the specific rule named, write the general solution. For resonance: one sentence (”F(a,b)=0F(a,b)=0, so multiply by xx”), then the corrected PI.

15-mark questions (Q6–Q8 slot): Show factorisation with the auxiliary equation, verify the factorisation. Write CF with correct argument form. For multi-term RHS: handle each PI term in a labelled sub-step. Verify the PI by substituting back (one line, or quote sympy if time is short). Box the general solution.

For the shift-identity questions (2013, 2022): explicitly write F(D+i,D+i)F(D+i,D'+i), simplify, identify the constant term F(i,i)0F(i,i)\ne0, expand the geometric series, and extract the real/imaginary part.

Practice Set

YearPaper/QMarksArchetypeOne-line hint
2025P2-Q5a10operator-pdeRHS ysinxy\sin x; undetermined coefficients y(asinx+bcosx)+(csinx+dcosx)y(a\sin x+b\cos x)+(c\sin x+d\cos x)
2024P2-Q5a10verify-pde-solutionChain rule: (iα)2=α2(i\alpha)^2=-\alpha^2; collect f+gf''+g''; match coefficients
2023P2-Q6a15operator-pde-boundaryRepeated root (D2D)2(D-2D')^2; f0f\equiv0 from x=0x=0 BC; g(u)=3/ug(u)=3/u from x=yx=y BC
2022P2-Q7a15operator-pdeSame operator/RHS as 2013; shift identity; F(i,i)=4F(i,i)=4; Im part gives the PI
2021P2-Q7a15operator-pdeFactor (DD)(D+D3)(D-D')(D+D'-3); CF includes e3xϕ2(yx)e^{3x}\phi_2(y-x); double resonance for ex+2ye^{x+2y}
2020P2-Q5d10operator-pdeThree roots; F(2,1)=0F(2,1)=0 (resonance); trig PI via D21D^2\to-1, D24D'^2\to-4, DD2DD'\to2
2018P2-Q7a15operator-pdeThree PIs; trig is resonant (F=0F=0); xx-multiply; polynomial is cubic; exponential F(3,4)=10F(3,4)=-10
2017P2-Q5a10operator-pdeRepeated root (DD)2(D-D')^2; three separate PIs; ex+2ye^{x+2y} non-resonant (F(1,2)=1F(1,2)=1)
2016P2-Q7a15operator-pdeThird-order cubic; m=2,1,1m=2,1,-1; resonance at m=1m=1 for ex+ye^{x+y}; PI =xex+y/2=-xe^{x+y}/2
2015P2-Q5b10operator-pdeFactor (DD)(D+2D)(D-D')(D+2D'); F(1,1)=0F(1,1)=0 resonance; apply non-resonant factor first
2014P2-Q5a10operator-pdeFactor (2DD)(D2D)(2D-D')(D-2D'); polynomial PI cubic ansatz; matching gives zp=3x3+6x2yz_p=3x^3+6x^2y
2013P2-Q6a15operator-pdex2sin(x+y)x^2\sin(x+y); shift identity; F(i,i)=4F(i,i)=4; Im of ei(x+y)Φe^{i(x+y)}\Phi with series

We've mapped all 13 years of this exam. Get new chapters, tools, and solutions as we release them — free.