← 2026 Paper 2

UPSC 2026 Maths Optional Paper 2 Q5b — Step-by-Step Solution

10 marks · Section B

Newton-Raphson method (convergence, geometric meaning) · Numerical Analysis · asked 7× in 14 yrs · Read the full method →

Question

By using the Newton-Raphson iterative formula, establish the formula

xi+1=13(2xi+Nxi2)x_{i+1} = \frac{1}{3}\left(2x_i + \frac{N}{x_i^2}\right)

to find the cube root of NN. Use this formula, if established, to find the cube root of 63 correct to four decimal places with the initial approximation 3⋅93\cdot9.

Technique

“Establish” means derive, not quote: choose f(x)=x3−Nf(x)=x^{3}-N (whose positive root is N1/3N^{1/3}), write Newton–Raphson xi+1=xi−f(xi)/f′(xi)x_{i+1}=x_i-f(x_i)/f'(x_i), and simplify. The numerical part is then three iterations from x0=3.9x_0=3.9; the marks in the second half are for displaying every iterate and for justifying the four-decimal claim — and the justification that actually holds water is a bracket α<N1/3<β\alpha<N^{1/3}<\beta, not “two iterates agreed”.

Solution

Step 1 — Choose the function and state Newton–Raphson.

To compute N1/3N^{1/3} (N>0N>0) we need the real root of

f(x)=x3−N,f′(x)=3x2.f(x)=x^{3}-N,\qquad f'(x)=3x^{2}.

ff is continuous, f(0)=−N<0f(0)=-N<0, f(x)→+∞f(x)\to+\infty, and f′(x)=3x2>0f'(x)=3x^{2}>0 for x>0x>0: so ff has exactly one positive zero, namely x=N1/3x=N^{1/3}, and it is a simple root (f′(N1/3)=3N2/3≠0f'(N^{1/3})=3N^{2/3}\neq0) — which is what entitles us to expect second-order convergence.

The Newton–Raphson iteration (the tangent at (xi,f(xi))(x_i,f(x_i)) met with the xx-axis) is

xi+1=xi−f(xi)f′(xi),f′(xi)≠0.x_{i+1}=x_i-\frac{f(x_i)}{f'(x_i)},\qquad f'(x_i)\neq0 .

Step 2 — Substitute and simplify.

xi+1=xi−xi3−N3xi2=3xi3−xi3+N3xi2=2xi3+N3xi2,x_{i+1}=x_i-\frac{x_i^{3}-N}{3x_i^{2}}=\frac{3x_i^{3}-x_i^{3}+N}{3x_i^{2}}=\frac{2x_i^{3}+N}{3x_i^{2}}, ∴xi+1=13(2xi+Nxi2),xi≠0.■\therefore\quad x_{i+1}=\frac{1}{3}\left(2x_i+\frac{N}{x_i^{2}}\right),\qquad x_i\neq0. \qquad\blacksquare

This is the required formula. Two remarks that cost one line each and earn marks:

Fixed point. With g(x)=13(2x+Nx−2)g(x)=\tfrac13\bigl(2x+Nx^{-2}\bigr),  g(N1/3)=13(2N1/3+N⋅N−2/3)=13⋅3N1/3=N1/3\ g(N^{1/3})=\tfrac13\bigl(2N^{1/3}+N\cdot N^{-2/3}\bigr)=\tfrac13\cdot 3N^{1/3}=N^{1/3} — the exact cube root is the fixed point, as it must be.

Order. g′(x)=13(2−2Nx−3)=23(1−N/x3)g'(x)=\tfrac13\bigl(2-2Nx^{-3}\bigr)=\tfrac23\bigl(1-N/x^{3}\bigr), so g′(N1/3)=0g'(N^{1/3})=0: the iteration is at least second order, i.e. quadratically convergent.

Step 3 — An exact error relation (this replaces hand-waving later).

Let r=N1/3r=N^{1/3}, so N=r3N=r^{3}, and let ei=xi−re_i=x_i-r. Then

xi+1−r=2xi3+r33xi2−r=2xi3−3rxi2+r33xi2.x_{i+1}-r=\frac{2x_i^{3}+r^{3}}{3x_i^{2}}-r=\frac{2x_i^{3}-3r x_i^{2}+r^{3}}{3x_i^{2}} .

The numerator has xi=rx_i=r as a double root; dividing twice by (xi−r)(x_i-r),

2x3−3rx2+r3=(x−r)2(2x+r),2x^{3}-3rx^{2}+r^{3}=(x-r)^{2}(2x+r), ∴ei+1=ei2 2xi+r3xi2 → xi→r  ei2r.(∗)\therefore\quad e_{i+1}=e_i^{2}\,\frac{2x_i+r}{3x_i^{2}}\ \xrightarrow[\ x_i\to r\ ]{}\ \frac{e_i^{2}}{r}. \tag{$\ast$}

Two consequences we will use:

  1. Sign. For xi>0x_i>0 we have 2xi+r>02x_i+r>0, so ei+1>0e_{i+1}>0 whenever ei≠0e_i\neq0: every iterate from x1x_1 onwards overshoots above rr. (Geometrically: ff is convex on x>0x>0, so the tangent line lies below the curve and meets the axis to the right of the root.)
  2. Monotonicity. For xi>rx_i>r, xi−xi+1=xi3−N3xi2>0x_i-x_{i+1}=\dfrac{x_i^{3}-N}{3x_i^{2}}>0, so x1>x2>x3>⋯>rx_1>x_2>x_3>\cdots>r: the iterates decrease monotonically to rr. This gives a rigorous upper bound for rr at every stage — and it is the honest half of the “four decimal places” claim.

Step 4 — The numerical work: N=63N=63, x0=3⋅9x_0=3\cdot9.

xi+1=13(2xi+63xi2).x_{i+1}=\frac13\left(2x_i+\frac{63}{x_i^{2}}\right).

Iteration 1. x0=3.9x_0=3.9, x02=15.21x_0^{2}=15.21, 6315.21=63001521=700169=4.1420118343\dfrac{63}{15.21}=\dfrac{6300}{1521}=\dfrac{700}{169}=4.1420118343,

x1=7.8+4.14201183433=11.94201183433=3.9806706114.x_1=\frac{7.8+4.1420118343}{3}=\frac{11.9420118343}{3}=3.9806706114 .

Iteration 2. x12=15.8457385168x_1^{2}=15.8457385168, 63x12=3.9758323623\dfrac{63}{x_1^{2}}=3.9758323623,

x2=7.9613412229+3.97583236233=11.93717358523=3.9790578617.x_2=\frac{7.9613412229+3.9758323623}{3}=\frac{11.9371735852}{3}=3.9790578617 .

Iteration 3. x22=15.8329014670x_2^{2}=15.8329014670, 63x22=3.9790559002\dfrac{63}{x_2^{2}}=3.9790559002,

x3=7.9581157235+3.97905590023=11.93717162373=3.9790572079.x_3=\frac{7.9581157235+3.9790559002}{3}=\frac{11.9371716237}{3}=3.9790572079 .

| ii | xix_i (10 d.p.) | ∣xi−xi−1∣|x_i-x_{i-1}| | |---|---|---| | 0 | 3.90000000003.9000000000 | — | | 1 | 3.98067061143.9806706114 | 8.067×10−28.067\times10^{-2} | | 2 | 3.97905786173.9790578617 | 1.613×10−31.613\times10^{-3} | | 3 | 3.97905720793.9790572079 | 6.538×10−76.538\times10^{-7} |

The successive differences fall as 8×10−2→1.6×10−3→6.5×10−78\times10^{-2}\to1.6\times10^{-3}\to6.5\times10^{-7} — the number of correct digits roughly doubles each step, exactly the quadratic behaviour predicted by (∗)(\ast) with asymptotic constant 1/r≈0.25131/r\approx0.2513. A fourth iteration would change x3x_3 only in about the thirteenth decimal (by (∗)(\ast), e4≈e32/re_4\approx e_3^{2}/r with e3≈1.1×10−13e_3\approx1.1\times10^{-13}), so the iteration has converged to the precision carried.

Step 5 — Why the answer is correct to four decimal places (the part worth doing properly).

Rounding x3x_3 blindly is not an argument, and neither is ”x2x_2 and x3x_3 agree to six places” — agreement of successive iterates is evidence, not proof (it fails, for instance, for a slowly convergent linear iteration). Here is a proof, and it takes two lines.

Upper bound. By Step 3(2), the iterates decrease monotonically to rr, so

r<x3=3.9790572079.r<x_3=3.9790572079 .

Lower bound. x↦x3x\mapsto x^{3} is strictly increasing on x>0x>0, so it suffices to test one cube:

(3.97905)3=62.999658<63⟹3.97905<r.(3.97905)^{3}=62.999658<63\quad\Longrightarrow\quad 3.97905<r .

Hence

3.97905<r<3.97906,3.97905<r<3.97906,

so r=3.9790…r=3.9790\ldots with the fifth decimal digit 55 followed by more, and

r∈(3.97905, 3.97915) ⟹ r rounded to four decimals=3.9791.r\in(3.97905,\,3.97915)\ \Longrightarrow\ r\ \text{rounded to four decimals}=3.9791 .

(For completeness the upper end of the four-decimal window also checks directly: (3.97915)3=63.004408>63(3.97915)^{3}=63.004408>63, so r<3.97915r<3.97915.)

∴633=3.9791 (correct to four decimal places).■\therefore\quad \sqrt[3]{63}=3.9791\ \text{(correct to four decimal places)}.\qquad\blacksquare

Step 6 — A trap worth naming. A candidate who stops at x1=3.98067x_1=3.98067 and rounds gets 3.98073.9807 — wrong in the third decimal. One iterate is never enough here because e1≈1.6×10−3e_1\approx1.6\times10^{-3}; it is e2≈6.5×10−7e_2\approx6.5\times10^{-7} that first clears the 5×10−55\times10^{-5} needed for four decimals. Two iterations suffice for the stated accuracy; the third is written down to demonstrate convergence.

Answer

  xi+1=xi−xi3−N3xi2=13(2xi+Nxi2);633=3.9791  (4 d.p.),x3=3.9790572079  \boxed{\;x_{i+1}=x_i-\frac{x_i^{3}-N}{3x_i^{2}}=\frac13\left(2x_i+\frac{N}{x_i^{2}}\right);\qquad \sqrt[3]{63}=3.9791\ \ (4\ \text{d.p.}),\quad x_3=3.9790572079\;}
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.