2. Iterate: $x_{n+1} = \mathrm{g}(x_n)$, starting from an estimate $x_1$
The same equation can be rearranged many ways, and the choice matters enormously โ some rearrangements converge and others do not.
Find three ways of writing $x^3 - 5x + 1 = 0$ in the form $x = \mathrm{g}(x)$.
Use $x_{n+1} = \dfrac{x_n^3+1}{5}$ with $x_1 = 0$ to find a root of $x^3 - 5x + 1 = 0$ to $4$ decimal places.
The equation $x^3 - 5x + 1 = 0$ also has a root near $2$. Use $x_{n+1} = \sqrt[3]{5x_n - 1}$ with $x_1 = 2$.
| Behaviour | What you see |
|---|---|
| Converges | The iterates settle towards one value; successive differences shrink |
| Diverges | The values grow without limit, or produce a calculator error |
| Oscillates and converges | Alternates above and below the root, closing in each time |
| Oscillates and diverges | Alternates, but the swings get bigger |
then round. Show at least the first four iterates as evidence.
Try $x_{n+1} = x_n^3 - 4x_n$ โฆ or rather, take the equation $x^3 - 5x + 1 = 0$ rearranged as $x = x^3 - 4x + 1$, with $x_1 = 0.2$.
2. The starting value is too far away
3. The formula becomes undefined โ a negative square root, or division by zero
Explain what happens using $x_{n+1} = \sqrt{x_n - 3}$ with $x_1 = 2$.
Use a change of sign first, to get a starting value close to the root.
The form
Rearrange to $x = \mathrm{g}(x)$.
The iteration
$x_{n+1} = \mathrm{g}(x_n)$ with a starting value.
ANS key
One button press per step.
Full accuracy
Never round between steps.
Why it works
The limit satisfies $L = \mathrm{g}(L)$.
Stop when
Successive values agree to the accuracy wanted.
Verify
Substitute back into $\mathrm{f}(x) = 0$.
One root only
Which one depends on $\mathrm{g}$ and on $x_1$.
Can diverge
A different rearrangement may fail.
Can be undefined
Watch for negative square roots.
Rearrange $x^3 - 7x + 2 = 0$ into the form $x = \mathrm{g}(x)$ in two different ways.
โถ Show solution
$7x = x^3 + 2$, so $x = \dfrac{x^3+2}{7}$
$x^3 = 7x - 2$, so $x = \sqrt[3]{7x-2}$
Use $x_{n+1} = \dfrac{x_n^3+2}{7}$ with $x_1 = 0$ to find $x_2$, $x_3$ and $x_4$.
โถ Show solution
$x_2 = \dfrac{0+2}{7} = 0.285714$
$x_3 = \dfrac{0.023324+2}{7} = 0.289046$
$x_4 = \dfrac{0.024148+2}{7} = 0.289164$
Give the root found in Question 2 to $3$ decimal places.
โถ Show solution
The iterates are settling at $0.28917$, so the root is $\mathbf{0.289}$ to $3$ d.p.
Use $x_{n+1} = \sqrt{\dfrac{10}{x_n + 1}}$ with $x_1 = 2$ to find $x_2$ and $x_3$.
โถ Show solution
$x_2 = \sqrt{10 \div 3} = \sqrt{3.333333} = 1.825742$
$x_3 = \sqrt{10 \div 2.825742} = \sqrt{3.538882} = 1.881195$
An iteration gives $x_4 = 1.4142$ and $x_5 = 1.4142$. What can you conclude?
โถ Show solution
The iterates agree to $4$ decimal places, so the sequence has converged and the root is $1.4142$ to $4$ d.p.
(Strictly, agreement to $4$ d.p. is strong evidence rather than proof; a change of sign check on either side of $1.4142$ would confirm it.)
Explain why $x_{n+1} = \sqrt{x_n - 10}$ with $x_1 = 3$ fails.
โถ Show solution
$x_2 = \sqrt{3 - 10} = \sqrt{-7}$, which is not a real number.
The iteration is undefined at the very first step, so it fails immediately. A starting value of at least $10$ would be needed to keep the expression under the root non-negative.
Use $x_{n+1} = \sqrt[3]{2x_n + 3}$ with $x_1 = 2$ to find the root to $3$ d.p.
โถ Show solution
$x_2 = \sqrt[3]{7} = 1.912931$
$x_3 = \sqrt[3]{6.825862} = 1.896935$
$x_4 = \sqrt[3]{6.793871} = 1.893967$
$x_5 = \sqrt[3]{6.787934} = 1.893415$
$x_6 = \sqrt[3]{6.786831} = 1.893313$
Converging to $\mathbf{1.893}$ to $3$ d.p.
Check: this solves $x^3 - 2x - 3 = 0$; and $1.893^3 - 2(1.893) - 3 = 6.783 - 6.786 = -0.003 \approx 0$ โ
Show that the iteration $x_{n+1} = x_n^2 - 2$ with $x_1 = 3$ diverges.
โถ Show solution
$x_2 = 9 - 2 = 7$
$x_3 = 49 - 2 = 47$
$x_4 = 2209 - 2 = 2207$
$x_5 = 4\,870\,847$
The values grow explosively, so the iteration diverges.
(The equation $x = x^2 - 2$ does have roots, at $x = 2$ and $x = -1$, but this iteration runs away from both unless started between them.)
Explain why two students using the same equation but different rearrangements might get different roots.
โถ Show solution
An iteration converges only to roots that are "attracting" for that particular function $\mathrm{g}$, and only from starting values close enough to them.
A different rearrangement gives a different $\mathrm{g}$, with a different pattern of which roots attract and which repel.
So one student's formula may converge to a root near $0$, while the other's converges to a root near $2$ โ both correct roots of the same equation, reached by different routes. Neither student can claim to have found all the roots without further work.
The equation $x^3 - 6x - 4 = 0$ is to be solved.
(a) Show there is a root between $2$ and $3$. (b) Show that the equation can be written as $x = \sqrt[3]{6x+4}$. (c) Use this iteration with $x_1 = 2.5$ to find the root to $3$ d.p. (d) Show that the rearrangement $x = \dfrac{x^3-4}{6}$ fails from the same starting value, and explain why.
โถ Show solution
(a) $\mathrm{f}(2) = 8 - 12 - 4 = -8 < 0$
$\mathrm{f}(3) = 27 - 18 - 4 = 5 > 0$
Sign change with $\mathrm{f}$ continuous, so a root lies in $(2, 3)$.
(b) $x^3 - 6x - 4 = 0 \;\Rightarrow\; x^3 = 6x + 4 \;\Rightarrow\; x = \sqrt[3]{6x+4}$ โ
(c) $x_1 = 2.5$
$x_2 = \sqrt[3]{19} = 2.668402$
$x_3 = \sqrt[3]{20.010410} = 2.714888$
$x_4 = \sqrt[3]{20.289331} = 2.727444$
$x_5 = \sqrt[3]{20.364666} = 2.730816$
$x_6 = \sqrt[3]{20.384896} = 2.731720$
$x_7 = \sqrt[3]{20.390319} = 2.731962$
The iterates are converging, so the root is $\mathbf{2.731}$ to $3$ d.p.
Verify: $2.731^3 - 6(2.731) - 4 = 20.369 - 16.386 - 4 = -0.017 \approx 0$ โ
(d) Using $x_{n+1} = \dfrac{x_n^3-4}{6}$ from $x_1 = 2.5$:
$x_2 = \dfrac{15.625-4}{6} = 1.937500$
$x_3 = \dfrac{7.273193-4}{6} = 0.545532$
$x_4 = \dfrac{0.162355-4}{6} = -0.639608$
$x_5 = \dfrac{-0.261625-4}{6} = -0.710277$
$x_6 = -0.726388$, $x_7 = -0.730545$ โฆ
The sequence has moved away from the root near $2.731$ and is settling on a different root, near $-0.732$.
Why: the cubic $x^3-6x-4$ has three roots, at approximately $-2.000$, $-0.732$ and $2.731$. This rearrangement makes the root at $2.731$ repelling and the one at $-0.732$ attracting, so it converges to the wrong one. The first rearrangement has the opposite behaviour.
This is precisely the failure NM3 asks you to recognise: a perfectly valid rearrangement can converge to a root other than the one you wanted, or to none at all.