← 2026 Paper 2
UPSC 2026 Maths Optional Paper 2 Q7b — Step-by-Step Solution
15 marks · Section B
Simpson's 1/3 and 3/8 rules · Numerical Analysis · asked 6× in 14 yrs · Read the full method →
Question
Evaluate the integral ∫00⋅8f(x)dx, where
f(x)=0⋅2+25x−200x2+675x3−900x4+400x5
by using single application of Simpson’s 83 rule.
(The raised dots are the Indian decimal point: 0⋅8=0.8, 0⋅2=0.2.)
Technique
“Single application” is the operative phrase: it means the basic (un-composited) Newton–Cotes rule on four equally spaced nodes, so the number of strips is fixed at n=3 and h=(b−a)/3 — not a free choice. Compute the four ordinates and apply 83h[f0+3f1+3f2+f3]. Then, because the integrand is an explicit polynomial, integrate it exactly and report the absolute and relative error: the 3/8 rule is exact only up to cubics and this integrand is quintic, so there is a genuine error, and quantifying it against the error term −803h5f(4)(ξ) is what converts a rote answer into a full-marks one.
Solution
Step 1 — The rule, and what “single application” fixes
Simpson’s 3/8 rule replaces f on [x0,x3] by the cubic through four equally spaced points x0,x1,x2,x3 (Newton–Cotes, n=3) and integrates that cubic exactly:
∫x0x3f(x)dx≈83h[f0+3f1+3f2+f3],h=3x3−x0,
with error term
E=−803h5f(4)(ξ)=−6480(b−a)5f(4)(ξ),ξ∈(x0,x3).(1)
A single application therefore uses exactly n=3 sub-intervals and 4 ordinates over the whole range [0,0.8]. Hence
a=0,b=0.8,h=3b−a=30.8=154=0.2666667.
A useful arithmetic shortcut, worth taking in the hall: the multiplier is
83h=83⋅3b−a=8b−a=80.8=0.1exactly.
Step 2 — The four ordinates
The nodes are xi=ih, i.e. 0, 154, 158, 54. Evaluating
f(x)=0.2+25x−200x2+675x3−900x4+400x5 (all values carried to 7 decimals; the exact rationals are given because they make the check reproducible):
| i | xi (exact) | xi | fi (exact) | fi |
|---|
| 0 | 0 | 0.0000000 | 51 | 0.2000000 |
| 1 | 154 | 0.2666667 | 3037543519 | 1.4327243 |
| 2 | 158 | 0.5333333 | 30375105923 | 3.4871769 |
| 3 | 54 | 0.8000000 | 12529 | 0.2320000 |
Sample working for f1 (the one candidates get wrong), using x1=154:
25x1=320,200x12=9128,675x13=3375675⋅64=564,900x14=2251024,400x15=3037516384,
f1=51+320−9128+564−2251024+3037516384=3037543519=1.4327243.
And for the endpoint,
f3=0.2+20−128+345.6−368.64+131.072=0.232.
Step 3 — Apply the rule
I3/8=83h[f0+3(f1+f2)+f3]=0.1[0.2+3(1.4327243+3.4871769)+0.232]
=0.1[0.2+3(4.9199012)+0.232]=0.1[0.2+14.7597036+0.232]=0.1×15.1917036,
I3/8=1.5191704(exactly 1687525636=1.5191703703).
Step 4 — The exact value (the integrand is a polynomial, so this is free)
∫00.8fdx=[0.2x+12.5x2−3200x3+4675x4−180x5+3200x6]00.8
=0.16+8−34.1333333+69.12−58.9824+17.4762667=1.6405333,
exactly 18753076=1.6405333.
Step 5 — Error, absolute and relative
E=Iexact−I3/8=18753076−1687525636=168752048=0.1213630,
εrel=IexactE=3076/18752048/16875=6921512=0.0739777=7.398%.
Why the error is there, and why it is this size. The 3/8 rule integrates the interpolating cubic exactly, so it is exact for every polynomial of degree ≤3. Our integrand is quintic; the terms −900x4 and +400x5 are precisely what the rule cannot see, and they carry the whole error.
Quantitatively, from (1) with
f(4)(x)=−21600+48000x(affine),
the crude bound uses max[0,0.8]∣f(4)∣=max{21600,16800}=21600:
∣E∣≤803h5(21600)=803(154)5(21600)=18752048=1.0922667,
which is satisfied but slack by exactly the factor 21600/2400=9. The sharp statement is available here because f(4) is affine and the 3/8 rule is symmetric about the midpoint of [a,b]: its degree-4 Peano kernel K is symmetric, so ∫tK(t)dt=xˉ∫K(t)dt with xˉ=2a+b, and for an affine f(4) the mean-value point is exactly the midpoint,
ξ=2a+b=0.4,f(4)(0.4)=−21600+19200=−2400,
E=−803h5f(4)(0.4)=−803(154)5(−2400)=168752048=0.1213630,
in exact agreement with the directly computed error. ■
Interpretation for the answer script (one line, worth writing): h=0.2667 is a coarse step for an integrand that swings from 0.2 up to 3.49 and back down to 0.232 on this range; a single application of any four-point rule cannot track that, and ≈7.4% error is the honest price. Composite application (two 3/8 panels, h=0.1333) would cut the error by the expected factor ≈25=32.
Answer
∫00.8fdx≈3/883h[f0+3f1+3f2+f3]=1.5191704,h=30.8=0.2666667
exact value =18753076=1.6405333,∣E∣=168752048=0.1213630,εrel=7.398%