← 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\displaystyle\int_{0}^{0\cdot8} f(x)\,dx, where

f(x)=0⋅2+25x−200x2+675x3−900x4+400x5f(x) = 0\cdot2 + 25x - 200x^2 + 675x^3 - 900x^4 + 400x^5

by using single application of Simpson’s 38\dfrac{3}{8} rule.

(The raised dots are the Indian decimal point: 0⋅8=0.80\cdot8 = 0.8, 0⋅2=0.20\cdot2 = 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=3n = 3 and h=(b−a)/3h = (b-a)/3 — not a free choice. Compute the four ordinates and apply 3h8[f0+3f1+3f2+f3]\frac{3h}{8}\bigl[f_0+3f_1+3f_2+f_3\bigr]. Then, because the integrand is an explicit polynomial, integrate it exactly and report the absolute and relative error: the 3/83/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 −3h580f(4)(ξ)-\frac{3h^5}{80}f^{(4)}(\xi) is what converts a rote answer into a full-marks one.

Solution

Step 1 — The rule, and what “single application” fixes

Simpson’s 3/83/8 rule replaces ff on [x0,x3][x_0,x_3] by the cubic through four equally spaced points x0,x1,x2,x3x_0,x_1,x_2,x_3 (Newton–Cotes, n=3n=3) and integrates that cubic exactly:

∫x0x3f(x) dx  ≈  3h8[f0+3f1+3f2+f3],h=x3−x03,\int_{x_0}^{x_3} f(x)\,dx \;\approx\; \frac{3h}{8}\Bigl[f_0 + 3f_1 + 3f_2 + f_3\Bigr],\qquad h = \frac{x_3-x_0}{3},

with error term

E=−3h580 f(4)(ξ)=−(b−a)56480f(4)(ξ),ξ∈(x0,x3).(1)E = -\frac{3h^5}{80}\,f^{(4)}(\xi) = -\frac{(b-a)^5}{6480}f^{(4)}(\xi),\qquad \xi \in (x_0,x_3). \tag{1}

A single application therefore uses exactly n=3n=3 sub-intervals and 44 ordinates over the whole range [0, 0.8][0,\,0.8]. Hence

a=0,b=0.8,h=b−a3=0.83=415=0.2666667.a = 0,\quad b = 0.8,\quad h = \frac{b-a}{3} = \frac{0.8}{3} = \frac{4}{15} = 0.2666667 .

A useful arithmetic shortcut, worth taking in the hall: the multiplier is

3h8=38⋅b−a3=b−a8=0.88=0.1exactly.\frac{3h}{8} = \frac{3}{8}\cdot\frac{b-a}{3} = \frac{b-a}{8} = \frac{0.8}{8} = 0.1 \quad\text{exactly}.

Step 2 — The four ordinates

The nodes are xi=ihx_i = ih, i.e. 0, 415, 815, 450,\ \tfrac{4}{15},\ \tfrac{8}{15},\ \tfrac{4}{5}. Evaluating f(x)=0.2+25x−200x2+675x3−900x4+400x5f(x) = 0.2 + 25x - 200x^2 + 675x^3 - 900x^4 + 400x^5 (all values carried to 77 decimals; the exact rationals are given because they make the check reproducible):

iixix_i (exact)xix_ifif_i (exact)fif_i
00000.00000000.000000015\tfrac150.20000000.2000000
11415\tfrac{4}{15}0.26666670.26666674351930375\tfrac{43519}{30375}1.43272431.4327243
22815\tfrac{8}{15}0.53333330.533333310592330375\tfrac{105923}{30375}3.48717693.4871769
3345\tfrac{4}{5}0.80000000.800000029125\tfrac{29}{125}0.23200000.2320000

Sample working for f1f_1 (the one candidates get wrong), using x1=415x_1 = \tfrac4{15}:

25x1=203,200x12=1289,675x13=675⋅643375=645,900x14=1024225,400x15=1638430375,25x_1 = \tfrac{20}{3},\quad 200x_1^2 = \tfrac{128}{9},\quad 675x_1^3 = \tfrac{675\cdot 64}{3375} = \tfrac{64}{5},\quad 900x_1^4 = \tfrac{1024}{225},\quad 400x_1^5 = \tfrac{16384}{30375}, f1=15+203−1289+645−1024225+1638430375=4351930375=1.4327243.f_1 = \tfrac15 + \tfrac{20}{3} - \tfrac{128}{9} + \tfrac{64}{5} - \tfrac{1024}{225} + \tfrac{16384}{30375} = \tfrac{43519}{30375} = 1.4327243 .

And for the endpoint,

f3=0.2+20−128+345.6−368.64+131.072=0.232.f_3 = 0.2 + 20 - 128 + 345.6 - 368.64 + 131.072 = 0.232 .

Step 3 — Apply the rule

I3/8=3h8[f0+3(f1+f2)+f3]=0.1[0.2+3(1.4327243+3.4871769)+0.232]I_{3/8} = \frac{3h}{8}\Bigl[f_0 + 3(f_1+f_2) + f_3\Bigr] = 0.1\Bigl[0.2 + 3(1.4327243 + 3.4871769) + 0.232\Bigr] =0.1[0.2+3(4.9199012)+0.232]=0.1[0.2+14.7597036+0.232]=0.1×15.1917036,= 0.1\bigl[0.2 + 3(4.9199012) + 0.232\bigr] = 0.1\bigl[0.2 + 14.7597036 + 0.232\bigr] = 0.1\times 15.1917036 ,  I3/8=1.5191704 (exactly 2563616875=1.5191703703‾).\boxed{\,I_{3/8} = 1.5191704\,}\qquad\left(\text{exactly } \frac{25636}{16875} = 1.5191703\overline{703}\right).

Step 4 — The exact value (the integrand is a polynomial, so this is free)

∫00.8 ⁣f dx=[0.2x+12.5x2−2003x3+6754x4−180x5+2003x6]00.8\int_0^{0.8}\! f\,dx = \left[0.2x + 12.5x^2 - \frac{200}{3}x^3 + \frac{675}{4}x^4 - 180x^5 + \frac{200}{3}x^6\right]_0^{0.8} =0.16+8−34.1333333+69.12−58.9824+17.4762667=1.6405333,= 0.16 + 8 - 34.1333333 + 69.12 - 58.9824 + 17.4762667 = 1.6405333 ,

exactly 30761875=1.6405333‾\dfrac{3076}{1875} = 1.640533\overline{3}.

Step 5 — Error, absolute and relative

E=Iexact−I3/8=30761875−2563616875=204816875=0.1213630,E = I_{\text{exact}} - I_{3/8} = \frac{3076}{1875} - \frac{25636}{16875} = \frac{2048}{16875} = 0.1213630 , εrel=EIexact=2048/168753076/1875=5126921=0.0739777=7.398%.\varepsilon_{\text{rel}} = \frac{E}{I_{\text{exact}}} = \frac{2048/16875}{3076/1875} = \frac{512}{6921} = 0.0739777 = 7.398\% .

Why the error is there, and why it is this size. The 3/83/8 rule integrates the interpolating cubic exactly, so it is exact for every polynomial of degree ≤3\le 3. Our integrand is quintic; the terms −900x4-900x^4 and +400x5+400x^5 are precisely what the rule cannot see, and they carry the whole error.

Quantitatively, from (1) with

f(4)(x)=−21600+48000x(affine),f^{(4)}(x) = -21600 + 48000x \quad(\text{affine}),

the crude bound uses max⁡[0,0.8]∣f(4)∣=max⁡{21600, 16800}=21600\max_{[0,0.8]}|f^{(4)}| = \max\{21600,\,16800\} = 21600:

∣E∣≤3h580 (21600)=380(415)5(21600)=20481875=1.0922667,|E| \le \frac{3h^5}{80}\,(21600) = \frac{3}{80}\Bigl(\tfrac{4}{15}\Bigr)^{5}(21600) = \frac{2048}{1875} = 1.0922667 ,

which is satisfied but slack by exactly the factor 21600/2400=921600/2400 = 9. The sharp statement is available here because f(4)f^{(4)} is affine and the 3/83/8 rule is symmetric about the midpoint of [a,b][a,b]: its degree-44 Peano kernel KK is symmetric, so ∫t K(t) dt=xˉ∫K(t) dt\int t\,K(t)\,dt = \bar{x}\int K(t)\,dt with xˉ=a+b2\bar x = \tfrac{a+b}{2}, and for an affine f(4)f^{(4)} the mean-value point is exactly the midpoint,

ξ=a+b2=0.4,f(4)(0.4)=−21600+19200=−2400,\xi = \frac{a+b}{2} = 0.4,\qquad f^{(4)}(0.4) = -21600 + 19200 = -2400, E=−3h580f(4)(0.4)=−380(415)5(−2400)=204816875=0.1213630,E = -\frac{3h^5}{80}f^{(4)}(0.4) = -\frac{3}{80}\Bigl(\frac{4}{15}\Bigr)^{5}(-2400) = \frac{2048}{16875} = 0.1213630 ,

in exact agreement with the directly computed error. ■\qquad\blacksquare

Interpretation for the answer script (one line, worth writing): h=0.2667h = 0.2667 is a coarse step for an integrand that swings from 0.20.2 up to 3.493.49 and back down to 0.2320.232 on this range; a single application of any four-point rule cannot track that, and ≈7.4%\approx 7.4\% error is the honest price. Composite application (two 3/83/8 panels, h=0.1333h=0.1333) would cut the error by the expected factor ≈25=32\approx 2^5 = 32.

Answer

  ∫00.8 ⁣f dx  ≈3/8  3h8[f0+3f1+3f2+f3]=1.5191704,h=0.83=0.2666667  \boxed{\;\int_0^{0.8}\!f\,dx \;\overset{3/8}{\approx}\; \frac{3h}{8}\bigl[f_0+3f_1+3f_2+f_3\bigr] = 1.5191704,\quad h=\tfrac{0.8}{3}=0.2666667\;}   exact value =30761875=1.6405333,∣E∣=204816875=0.1213630,εrel=7.398%  \boxed{\;\text{exact value } = \tfrac{3076}{1875} = 1.6405333,\qquad |E| = \tfrac{2048}{16875} = 0.1213630,\qquad \varepsilon_{\text{rel}} = 7.398\%\;}
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 — 14 years, mapped. Get the take-away PDF free.