The math optional, made finite. Daily Practice

Wave equation

At a Glance

Why This Chapter Matters

The wave equation appears in 6 of the last 11 years and is essentially always a 15- or 20-mark question — the highest mark range in Paper 2’s PDE cluster. Every question uses the same algorithm: separate variables, quantise the spatial modes using the fixed-end boundary conditions, apply the two initial conditions to determine the Fourier coefficients, and write the series. The variation between questions lies entirely in the initial conditions (zero vs. non-zero displacement, zero vs. non-zero velocity, polynomial vs. trigonometric IC), which determines which coefficients survive and how much integration is required. Mastering the algorithm once covers all seven appearances.

Minimum Theory

Wave equation. The transverse displacement u(x,t)u(x,t) of a tightly stretched string of length LL satisfies utt=c2uxx,0<x<L,  t>0,u_{tt}=c^2 u_{xx},\qquad 0<x<L,\;t>0, with c2=T/μc^2=T/\mu (tension/mass-per-unit-length). Boundary conditions u(0,t)=u(L,t)=0u(0,t)=u(L,t)=0 (fixed ends). Initial conditions: u(x,0)=f(x)u(x,0)=f(x) (initial shape) and ut(x,0)=g(x)u_t(x,0)=g(x) (initial velocity).

Separation of variables. Set u=X(x)T(t)u=X(x)T(t). The BCs force X(0)=X(L)=0X(0)=X(L)=0 and the separation constant must be negative (call it λ2-\lambda^2) to get oscillatory solutions. The normalised modes are Xn(x)=sin(nπx/L)X_n(x)=\sin(n\pi x/L), n=1,2,n=1,2,\ldots, with corresponding time factors Tn(t)=Ancos(nπct/L)+Bnsin(nπct/L)T_n(t)=A_n\cos(n\pi ct/L)+B_n\sin(n\pi ct/L). General solution: u(x,t)=n=1sinnπxL ⁣[AncosnπctL+BnsinnπctL].u(x,t)=\sum_{n=1}^\infty\sin\frac{n\pi x}{L}\!\left[A_n\cos\frac{n\pi ct}{L}+B_n\sin\frac{n\pi ct}{L}\right].

IC rule. Apply u(x,0)=f(x)u(x,0)=f(x): An=2L0Lf(x)sinnπxLdxA_n=\dfrac{2}{L}\displaystyle\int_0^L f(x)\sin\frac{n\pi x}{L}\,dx.

Apply ut(x,0)=g(x)u_t(x,0)=g(x): nπcLBn=2L0Lg(x)sinnπxLdx\dfrac{n\pi c}{L}B_n=\dfrac{2}{L}\displaystyle\int_0^L g(x)\sin\frac{n\pi x}{L}\,dx, i.e. Bn=2nπc0Lg(x)sinnπxLdxB_n=\dfrac{2}{n\pi c}\displaystyle\int_0^L g(x)\sin\frac{n\pi x}{L}\,dx.

Zero IC simplifications: If f(x)=0f(x)=0 then An=0A_n=0 for all nn (pure sin\sin in time). If g(x)=0g(x)=0 then Bn=0B_n=0 (pure cos\cos in time). If ff is already one eigenmode asin(n0πx/L)a\sin(n_0\pi x/L), no series is needed — only An0=aA_{n_0}=a survives.

IBP for polynomial ICs. For f(x)=x(Lx)f(x)=x(L-x) on (0,L)(0,L): 0Lx(Lx)sin(nπx/L)dx=2L3n3π3[1(1)n]\int_0^L x(L-x)\sin(n\pi x/L)\,dx=\dfrac{2L^3}{n^3\pi^3}[1-(-1)^n], so only odd nn contribute (1/n31/n^3 decay). For g(x)=λx(Lx)g(x)=\lambda x(L-x): the BnB_n coefficient is 8λL3n4π4c\dfrac{8\lambda L^3}{n^4\pi^4 c} for odd nn. For f(x)=xf(x)=x on (0,L)(0,L): An=2L(1)n+1nπA_n=\dfrac{2L(-1)^{n+1}}{n\pi} (slow 1/n1/n decay, alternating). For constant g=1g=1: Bn=4Ln2π2cB_n=\dfrac{4L}{n^2\pi^2 c} for odd nn only.

Triangular IC. For a string plucked at x=ax=a to height hh: bn=2hL2π2n2a(La)sinnπaLb_n=\dfrac{2hL^2}{\pi^2 n^2 a(L-a)}\sin\dfrac{n\pi a}{L}. Harmonics for which sin(nπa/L)=0\sin(n\pi a/L)=0 are absent (nodes at the pluck point).

Standing wave modes n=1,2,3: \bullet = fixed node at endpoints; \circ = interior node; each mode oscillates independently with frequency f_n=nc/(2L).

Question Archetypes

All wave-equation questions in the corpus share the same archetype.

ArchetypeYou are seeing this when…
vibrating-stringa string of length LL with fixed ends and given initial displacement/velocity; find u(x,t)u(x,t)

vibrating-string (7 question(s); 2013, 2014, 2014, 2017, 2020, 2021, 2023)

Recognition Cues

Solution Template

  1. Write the general sine-series solution u=n1sin(nπx/L)[Ancosωnt+Bnsinωnt]u=\sum_{n\ge1}\sin(n\pi x/L)[A_n\cos\omega_n t+B_n\sin\omega_n t] where ωn=nπc/L\omega_n=n\pi c/L.
  2. Set t=0t=0: Ansin(nπx/L)=f(x)\sum A_n\sin(n\pi x/L)=f(x). Compute AnA_n by Fourier projection or by inspection if ff is already a finite sum of modes.
  3. Differentiate in tt and set t=0t=0: Bnωnsin(nπx/L)=g(x)\sum B_n\omega_n\sin(n\pi x/L)=g(x). Compute Bn=2ωnL0Lgsin(nπx/L)dxB_n=\frac{2}{\omega_n L}\int_0^L g\sin(n\pi x/L)\,dx.
  4. For polynomial ICs: use two IBPs (bracket terms vanish at fixed ends); collect odd/even nn cases.
  5. Write the final boxed series. Note which harmonics vanish (e.g. nn divisible by 3 if plucked at L/3L/3).

Worked Example(s)

2013 Paper 2, 2013-P2-Q6c (20 marks)

String length ll, fixed ends, released from rest; initial velocity λx(lx)\lambda x(l-x).

f(x)=0f(x)=0 → all An=0A_n=0. Apply g(x)=λx(lx)g(x)=\lambda x(l-x): Bnnπcl=2λl0lx(lx)sinnπxldx=2λl2l3[1(1)n]n3π3.B_n\cdot\frac{n\pi c}{l}=\frac{2\lambda}{l}\int_0^l x(l-x)\sin\frac{n\pi x}{l}\,dx=\frac{2\lambda}{l}\cdot\frac{2l^3\cdot[1-(-1)^n]}{n^3\pi^3}.

For even nn: Bn=0B_n=0. For odd nn: Bn=8λl3n4π4cB_n=\dfrac{8\lambda l^3}{n^4\pi^4 c}.

  u(x,t)=8λl3π4cn odd1n4sinnπxlsinnπctl.  \boxed{\;u(x,t)=\frac{8\lambda l^3}{\pi^4 c}\sum_{n\text{ odd}}\frac{1}{n^4}\sin\frac{n\pi x}{l}\sin\frac{n\pi ct}{l}.\;}

The 1/n41/n^4 decay signals a quadratic IC; the odd-harmonics-only pattern follows from x(lx)x(l-x) being even about l/2l/2.


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

String length π\pi, c=1c=1; zero velocity; initial shape k(sinxsin2x)k(\sin x-\sin 2x).

g(x)=0g(x)=0 → all Bn=0B_n=0. The IC f(x)=ksinxksin2xf(x)=k\sin x-k\sin 2x is already a finite sum of eigenmodes (L=πL=\pi, so basis is sinnx\sin nx): by inspection A1=kA_1=k, A2=kA_2=-k, all others zero.

  u(x,t)=ksinxcostksin2xcos2t.  \boxed{\;u(x,t)=k\sin x\cos t-k\sin 2x\cos 2t.\;}

No Fourier integral required; the solution is a finite sum and stays finite for all tt.


2014 Paper 2, 2014-P2-Q8a (15 marks)

utt=uxxu_{tt}=u_{xx} on (0,1)(0,1); u(x,0)=0u(x,0)=0; ut(x,0)=x2u_t(x,0)=x^2.

L=c=1L=c=1. An=0A_n=0. Compute BnB_n from x2x^2 via two IBPs (using x2sin(nπx)dx=x2cos(nπx)/(nπ)+2xsin(nπx)/(nπ)2+2cos(nπx)/(nπ)3\int x^2\sin(n\pi x)\,dx=-x^2\cos(n\pi x)/(n\pi)+2x\sin(n\pi x)/(n\pi)^2+2\cos(n\pi x)/(n\pi)^3):

  u(x,t)=n=1Bnsin(nπt)sin(nπx),  \boxed{\;u(x,t)=\sum_{n=1}^\infty B_n\sin(n\pi t)\sin(n\pi x),\;}

with the BnB_n as above. Unlike the x(lx)x(l-x) case, both odd and even harmonics contribute (because x2x^2 lacks the midpoint symmetry).


2017 Paper 2, 2017-P2-Q8a (20 marks)

(i) Derive the general solution by separation of variables. (ii) Boundary conditions y(0,t)=y(l,t)=0y(0,t)=y(l,t)=0, yt(x,0)=0y_t(x,0)=0, y(x,0)=asin(πx/l)y(x,0)=a\sin(\pi x/l).

Part (i). Set y=X(x)G(t)y=X(x)G(t); show the separation constant must be negative k2-k^2; write y=(Acoskx+Bsinkx)(Ccosckt+Esinckt)y=(A\cos kx+B\sin kx)(C\cos ckt+E\sin ckt). State why +k2+k^2 (exponential in xx) is rejected. The appropriate solution is the oscillatory form with negative constant.

Part (ii). BCs: A=0A=0, k=nπ/lk=n\pi/l. IC yt=0y_t=0: En=0E_n=0. IC y(x,0)=asin(πx/l)y(x,0)=a\sin(\pi x/l) is already the n=1n=1 mode: C1=aC_1=a, all others zero.

  y(x,t)=asinπxlcosπctl.  \boxed{\;y(x,t)=a\sin\frac{\pi x}{l}\cos\frac{\pi ct}{l}.\;}

The initial shape is a single eigenmode so no Fourier series arises — only the fundamental survives.


2020 Paper 2, 2020-P2-Q8a (20 marks)

String length ll, plucked at x=l/3x=l/3 to height hh, released from rest.

g(x)=0g(x)=0Bn=0B_n=0. The initial shape is the triangular function: f(x)={3hx/l,0xl/3,3h(lx)/(2l),l/3xl.f(x)=\begin{cases}3hx/l,&0\le x\le l/3,\\3h(l-x)/(2l),&l/3\le x\le l.\end{cases}

Fourier sine coefficients (integrate piecewise + IBP on each segment): bn=9hπ2n2sinnπ3.b_n=\frac{9h}{\pi^2 n^2}\sin\frac{n\pi}{3}.

Since sin(nπ/3)=0\sin(n\pi/3)=0 for n=3,6,9,n=3,6,9,\ldots, every third harmonic vanishes (nodes at the pluck point).

  y(x,t)=9hπ2n=11n2sinnπ3sinnπxlcosnπctl.  \boxed{\;y(x,t)=\frac{9h}{\pi^2}\sum_{n=1}^\infty\frac{1}{n^2}\sin\frac{n\pi}{3}\sin\frac{n\pi x}{l}\cos\frac{n\pi ct}{l}.\;}


2021 Paper 2, 2021-P2-Q6a (20 marks)

a2uxx=utta^2u_{xx}=u_{tt} on (0,L)(0,L); u(x,0)=x(Lx)/4u(x,0)=x(L-x)/4; ut(x,0)=0u_t(x,0)=0.

Bn=0B_n=0 (zero velocity). Compute AnA_n from x(Lx)/4x(L-x)/4; using the standard integral 0Lx(Lx)sin(nπx/L)dx=2L3[1(1)n]/(n3π3)\int_0^L x(L-x)\sin(n\pi x/L)\,dx=2L^3[1-(-1)^n]/(n^3\pi^3), only odd nn contribute: An=2L144L3n3π3=2L2n3π3(n odd).A_n=\frac{2}{L}\cdot\frac{1}{4}\cdot\frac{4L^3}{n^3\pi^3}=\frac{2L^2}{n^3\pi^3}\quad(n\text{ odd}).

  u(x,t)=n odd2L2n3π3sin ⁣nπxLcos ⁣anπtL.  \boxed{\;u(x,t)=\sum_{n\text{ odd}}\frac{2L^2}{n^3\pi^3}\sin\!\frac{n\pi x}{L}\cos\!\frac{an\pi t}{L}.\;}

Same integral structure as 2013-P2-Q6c (quadratic IC, odd harmonics, 1/n31/n^3 decay from shape — noting here it’s displacement, not velocity, so decay is 1/n31/n^3 not 1/n41/n^4).


2023 Paper 2, 2023-P2-Q7c (20 marks)

a2uxx=utta^2u_{xx}=u_{tt} on (0,L)(0,L); u(x,0)=xu(x,0)=x; ut(x,0)=1u_t(x,0)=1.

Both ICs non-zero; both AnA_n and BnB_n survive.

AnA_n: Fourier sine series of xx on (0,L)(0,L): An=2L(1)n+1nπA_n=\dfrac{2L(-1)^{n+1}}{n\pi} (slow 1/n1/n decay, alternating sign — because f(L)=L0f(L)=L\ne0).

BnB_n: Fourier sine series of constant 11: nπaLBn=2[1(1)n]nπ\dfrac{n\pi a}{L}B_n=\dfrac{2[1-(-1)^n]}{n\pi}, giving Bn=4Ln2π2aB_n=\dfrac{4L}{n^2\pi^2 a} for odd nn, 00 for even.

  u=n=12L(1)n+1nπsinnπxLcosnπatL+n odd4Ln2π2asinnπxLsinnπatL.  \boxed{\;u=\sum_{n=1}^\infty\frac{2L(-1)^{n+1}}{n\pi}\sin\frac{n\pi x}{L}\cos\frac{n\pi at}{L}+\sum_{n\text{ odd}}\frac{4L}{n^2\pi^2 a}\sin\frac{n\pi x}{L}\sin\frac{n\pi at}{L}.\;}

Common Traps


Marks-Aware Writing

15-mark questions (2014-Q7a, 2014-Q8a): Four steps — write the general solution, apply the zero-velocity or zero-displacement IC (one sentence), compute the Fourier coefficient integral (show the IBP setup), write the final series. For Q7a: identify that the IC is already a finite sum; no integral needed (2 marks saved).

20-mark questions (2013, 2017, 2020, 2021, 2023): The mark scheme expects all five components: PDE + general solution derivation (step 1), BC application (step 2), first IC (step 3), second IC (step 4), assembled answer (step 5). For the 2017 two-part question, Part (i) derivation is worth ~8 marks on its own; state the rejection of the positive constant explicitly.

Practice Set

YearPaper/QMarksOne-line hint
2024P2-Q5a10Not a BVP — verify directly: compute uxx,uyy,uttu_{xx},u_{yy},u_{tt} via chain rule; show uxx+uyy=utt/c2u_{xx}+u_{yy}=u_{tt}/c^2 iff α2=1k2/c2\alpha^2=1-k^2/c^2 (the condition given)

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