1.09c - Fixed point iteration and cobweb diagrams

1.09c - Fixed point iteration and cobweb diagrams

Can you turn an awkward cubic into a calculation that improves itself? Under 1.09c, you will rearrange an equation, compute successive approximations and draw the staircase or cobweb that explains their movement. You will also check what the resulting decimal actually establishes.

An equation becomes a recurrence

Try repeatedly taking the cube root of twice your last answer plus 33, starting with 22. The outputs begin 1.9129311.912931, 1.8969351.896935, 1.8939671.893967: they seem to settle. Which equation are they solving?

You already know

A recurrence defines each new term from earlier terms. Rearrangements must preserve the equation and its domain. A real cube root is defined for negative as well as positive inputs.

The number being approached should reproduce itself under the rule.

Fixed point and iterate

A fixed point aa of gg satisfies g(a)=ag(a)=a. An iteration applies the rule once; its output is an iterate. In xn+1=g(xn)x_{n+1}=g(x_n), the subscript counts steps, not powers.

If a sequence converges to aa and gg is continuous there, taking limits in the recurrence gives a=g(a)a=g(a). That implication does not promise convergence.

Fixed point iteration

Rearrange the equation as x=g(x)x=g(x), choose a start x0x_0, then calculate

xn+1=g(xn),n=0,1,2,….x_{n+1}=g(x_n),\qquad n=0,1,2,\ldots.

Check every fixed point against the initial equation and domain.

Not in the formula sheet: learn it.

Begin with a rearrangement that makes the leading power easy to undo.

Worked example 1: constructing the rule

Question. Rearrange x3−2x−3=0x^3-2x-3=0 into a cube-root iteration and calculate x1x_1 from x0=2x_0=2.

Plan. Isolate the cube before taking a real cube root.

x3=2x+3move termsx=2x+33real cube rootxn+1=2xn+33turn equality into a rulex1=73≈1.912931\begin{aligned} x^3&=2x+3 && \text{move terms}\\ x&=\sqrt[3]{2x+3} && \text{real cube root}\\ x_{n+1}&=\sqrt[3]{2x_n+3} && \text{turn equality into a rule}\\ x_1&=\sqrt[3]{7}\approx1.912931 \end{aligned}

Why a cube root? Cubing reverses this step for every real input, so no root is lost.

Answer. Use xn+1=2xn+33x_{n+1}=\sqrt[3]{2x_n+3}; the first iterate is 1.9129311.912931 to six decimal places.

Check. Cubing x1x_1 gives 77, the value of 2x0+32x_0+3, not 2x1+32x_1+3.

Change only the constant from 33 to 55, keeping x0=2x_0=2.

Worked example 2: changing the constant

Question. Construct the same type of iteration for x3−2x−5=0x^3-2x-5=0, starting at x0=2x_0=2, and calculate x1x_1.

Plan. Isolate x3x^3 and evaluate the first output.

x3=2x+5xn+1=2xn+53x1=93≈2.080084\begin{aligned} x^3&=2x+5\\ x_{n+1}&=\sqrt[3]{2x_n+5}\\ x_1&=\sqrt[3]{9}\approx2.080084 \end{aligned}

Why keep the start unchanged? The changed output now comes entirely from the changed constant.

Answer. Use xn+1=2xn+53x_{n+1}=\sqrt[3]{2x_n+5}, with x1=2.080084x_1=2.080084 to six decimal places.

Check. Cubing the first output gives 9=2x0+59=2x_0+5.

The related equation x3−2x−4=0x^3-2x-4=0 has fixed point 22; from 11, its first cube-root iterate is 63≈1.817121\sqrt[3]{6}\approx1.817121.

Equivalent algebra can behave very differently. Rearranging the first cubic as x=(x3−3)/2x=(x^3-3)/2 gives 2,2.5,6.3125,…2,2.5,6.3125,\ldots. Both rules have the same fixed points, but one start does not produce the same sequence.

Watch out

Wrong: “A valid rearrangement must converge.” Right: algebraic equivalence preserves solutions; it does not guarantee a useful sequence.

Make the constant one unit larger again, keeping the same start.

Compute a reproducible sequence

The rule uses the previous unrounded value, even when your written table displays fewer digits. This small distinction prevents avoidable drift.

Worked example 3: using the previous output

Question. Use xn+1=2xn+33x_{n+1}=\sqrt[3]{2x_n+3}, x0=2x_0=2, to find x1,x2,x3x_1,x_2,x_3.

Plan. Retain full precision between successive evaluations.

x1=73=1.912931183…x2=2(1.912931183…)+33=1.896935259…x3=2(1.896935259…)+33=1.893967062…\begin{aligned} x_1&=\sqrt[3]{7}=1.912931183\ldots\\ x_2&=\sqrt[3]{2(1.912931183\ldots)+3}=1.896935259\ldots\\ x_3&=\sqrt[3]{2(1.896935259\ldots)+3}=1.893967062\ldots \end{aligned}

Why retain the hidden digits? The rounded printed output is not the number stored for the next step.

Answer. To six decimal places: 1.912931, 1.896935, 1.8939671.912931,\ 1.896935,\ 1.893967.

Check. The next four displayed values are 1.893415, 1.893313, 1.893294, 1.8932901.893415,\ 1.893313,\ 1.893294,\ 1.893290; the sequence is settling near 1.893291.89329.

Change only the starting value to 1.51.5. The recurrence stays the same.

Worked example 4: a different start

Question. Calculate the first three iterates from x0=1.5x_0=1.5.

Plan. Enter 1.51.5, then repeatedly evaluate the cube root of 22 times ANS plus 33.

x1=63=1.817120593…x2=2x1+33=1.879015742…x3=2x2+33=1.890630836…\begin{aligned} x_1&=\sqrt[3]{6}=1.817120593\ldots\\ x_2&=\sqrt[3]{2x_1+3}=1.879015742\ldots\\ x_3&=\sqrt[3]{2x_2+3}=1.890630836\ldots \end{aligned}

Why use ANS? It carries the full preceding result into every occurrence of the input.

Answer. To six decimal places: 1.817121, 1.879016, 1.8906311.817121,\ 1.879016,\ 1.890631.

Check. These approach the same neighbourhood from below, whereas the start 22 approached from above.

A spreadsheet does the same job: put the start in cell A2, calculate the cube root of 22 times A2 plus 33 in A3, then fill down with a relative reference to the preceding row. Keep a step-number column. Use a real cube-root operation if the argument can be negative; software fractional powers do not always handle negative inputs as intended.

Exam tip

State the recurrence, start and enough successive iterates to show your method. Retain calculator precision, then give the requested final accuracy. Set radians if the iteration contains trigonometric functions.

For comparison, the exact recurrence xn+1=(xn+8)/3x_{n+1}=(x_n+8)/3 from 11 gives 3,11/3,35/93,11/3,35/9. Each output supplies the next input, just as in the cubic calculation.

Continue this lesson.

You’re previewing NovaLearn without an account. Create a free account to read the rest.

Create free account