The Central Idea
A number \(r\) is a fixed point of a function \(g\) if
\[ \boxed{g(r)=r}. \]Fixed-point iteration begins with an initial estimate \(x_0\) and repeatedly applies
\[ \boxed{x_{n+1}=g(x_n)}. \]If the sequence converges to a limit \(r\), and \(g\) is continuous at \(r\), then
\[ r=g(r). \]The limiting value is therefore a fixed point of \(g\).
What Is a Fixed Point?
A fixed point is an input that a function sends back to itself.
For example, let
\[ g(x)=\frac{x+6}{2}. \]A fixed point satisfies
\[ x=\frac{x+6}{2}. \]Solving,
\[ 2x=x+6, \] \[ \boxed{x=6}. \]Indeed,
\[ g(6)=\frac{6+6}{2}=6. \]Fixed Points and Equations
Many equations can be rewritten so that their solutions become fixed points.
Suppose we want to solve
\[ f(x)=0. \]If the equation can be rearranged into
\[ \boxed{x=g(x)}, \]then every solution of the original equation that preserves the rearrangement is a fixed point of \(g\).
From Equation to Iteration Rule
Consider
\[ x^3+x-1=0. \]Rearrange:
\[ x^3+x=1, \] \[ x(x^2+1)=1, \] \[ \boxed{x=\frac{1}{x^2+1}}. \]This suggests the iteration
\[ \boxed{ x_{n+1}=\frac{1}{x_n^2+1} }. \]Example: Beginning the Iteration
Use
\[ x_{n+1}=\frac{1}{x_n^2+1} \]with
\[ x_0=0.5. \]Then
\[ x_1 = \frac{1}{0.5^2+1} = \frac{1}{1.25} = 0.8. \]Next,
\[ x_2 = \frac{1}{0.8^2+1} = \frac{1}{1.64} \approx0.609756. \]Continuing the Iteration
Continue:
\[ x_3 \approx \frac{1}{0.609756^2+1} \approx0.729293, \] \[ x_4\approx0.652999, \] \[ x_5\approx0.701033. \]The sequence alternates around the root while gradually moving closer.
The real solution is approximately
\[ \boxed{x\approx0.682328}. \]An Iteration Table
Fixed-point iterations are often organized in a table.
| \(n\) | \(x_n\) |
|---|---|
| \(0\) | \(0.500000\) |
| \(1\) | \(0.800000\) |
| \(2\) | \(0.609756\) |
| \(3\) | \(0.729293\) |
| \(4\) | \(0.652999\) |
| \(5\) | \(0.701033\) |
A longer table makes the convergence easier to observe.
Fixed-Point Iteration Is Repeated Substitution
If
\[ x_{n+1}=g(x_n), \]then
\[ x_1=g(x_0), \] \[ x_2=g(g(x_0)), \] \[ x_3=g(g(g(x_0))), \]and so on.
Each new estimate is created by substituting the previous estimate into the same function.
Graphical Meaning of a Fixed Point
A fixed point satisfies
\[ y=g(x) \]and
\[ y=x \]at the same point.
Therefore, fixed points are the \(x\)-coordinates where the graphs of
\[ y=g(x) \]and
\[ y=x \]intersect.
The Line \(y=x\)
The line
\[ y=x \]plays a special role in fixed-point iteration.
At a fixed point \(r\),
\[ g(r)=r, \]so the point
\[ (r,r) \]lies on both \(y=g(x)\) and \(y=x\).
Cobweb Interpretation
Fixed-point iteration can be visualized graphically.
Starting at \(x_0\):
- move vertically to the graph \(y=g(x)\);
- move horizontally to the line \(y=x\);
- repeat.
The resulting path is often called a cobweb diagram.
If the path contracts toward an intersection, the iteration converges.
Convergence
A fixed-point iteration converges when the sequence
\[ x_0,\quad x_1,\quad x_2,\quad\ldots \]approaches a limiting value \(r\).
Symbolically,
\[ \boxed{x_n\to r}. \]If \(g\) is continuous at \(r\), then the limit must satisfy
\[ r=g(r). \]Divergence
An iteration diverges when the successive values do not approach the desired finite fixed point.
Divergence may appear as:
- unbounded growth;
- persistent oscillation;
- cycling among several values;
- movement toward another fixed point;
- failure of the formula because an operation becomes undefined.
The Same Equation Can Produce Different Iterations
A crucial feature of fixed-point methods is that one equation can often be rearranged in several ways.
Consider
\[ x^2-3=0. \]One rearrangement is
\[ x=\frac{3}{x}, \]giving
\[ x_{n+1}=\frac{3}{x_n}. \]Another is
\[ x=\frac12\left(x+\frac3x\right), \]giving a very different iteration.
A Poor Rearrangement Can Fail
Use
\[ x_{n+1}=\frac{3}{x_n} \]with
\[ x_0=1. \]Then
\[ x_1=3, \] \[ x_2=1, \] \[ x_3=3, \] \[ x_4=1. \]The sequence cycles between \(1\) and \(3\).
It does not converge to \(\sqrt3\).
A Better Rearrangement
Now use
\[ \boxed{ x_{n+1} = \frac12\left(x_n+\frac3{x_n}\right) } \]with
\[ x_0=1. \]Then
\[ x_1=2, \] \[ x_2=1.75, \] \[ x_3\approx1.732143, \] \[ x_4\approx1.732051. \]This iteration converges rapidly to
\[ \boxed{\sqrt3\approx1.732051}. \]Why the Rearrangement Matters
Algebraically equivalent formulas can behave very differently when used as iteration rules.
The equation
\[ x=g(x) \]determines the fixed point, but the shape of \(g\) near that point strongly affects whether nearby iterates move toward or away from it.
The Derivative Test: The Main Local Idea
Near a fixed point \(r\), the magnitude of
\[ g'(r) \]gives important information about local behavior.
Roughly:
\[ \boxed{|g'(r)|<1} \]tends to pull nearby points toward the fixed point, while
\[ \boxed{|g'(r)|>1} \]tends to push nearby points away.
Why \(|g'(r)|<1\) Encourages Convergence
If \(x_n\) is near the fixed point \(r\), then a linear approximation gives
\[ g(x_n)-g(r) \approx g'(r)(x_n-r). \]Since
\[ g(r)=r, \]we obtain
\[ x_{n+1}-r \approx g'(r)(x_n-r). \]Therefore, if
\[ |g'(r)|<1, \]the error is locally multiplied by a factor whose magnitude is less than \(1\).
Positive and Negative Derivatives
If
\[ 0If
\[ -1Example of Alternating Convergence
Recall
\[ g(x)=\frac{1}{x^2+1}. \]Then
\[ g'(x) = -\frac{2x}{(x^2+1)^2}. \]Near the positive fixed point \(r\approx0.682328\), the derivative is negative with magnitude less than \(1\).
This is consistent with the observed pattern of alternating convergence.
When \(|g'(r)|>1\)
If
\[ |g'(r)|>1, \]small errors near the fixed point are typically enlarged rather than reduced.
The fixed point is then locally repelling for the iteration.
A nearby starting value will generally not converge to that fixed point.
The Borderline Case \(|g'(r)|=1\)
If
\[ |g'(r)|=1, \]the simple derivative test does not decide the outcome.
Higher-order behavior may produce convergence, divergence, or neutral behavior.
More analysis is needed.
A Stronger Interval Condition
A useful sufficient condition is available when an interval
\[ I=[a,b] \]satisfies both:
- \(g(x)\in I\) for every \(x\in I\);
- \(|g'(x)|\le k<1\) for every \(x\in I\).
Under these conditions, the iteration has a unique fixed point in the interval and starting values in the interval converge to it.
This is a standard contraction principle.
Mapping an Interval into Itself
The condition
\[ g(I)\subseteq I \]means that once the iteration enters the interval, the update rule keeps future iterates inside it.
This helps prevent the sequence from escaping into an unrelated region.
Contraction
A function acts like a contraction on an interval when it consistently reduces distances.
A derivative bound
\[ |g'(x)|\le k<1 \]provides one common way to establish this behavior.
The smaller \(k\) is, the stronger the contraction tends to be.
Estimated Rate of Convergence
Near a fixed point,
\[ |x_{n+1}-r| \approx |g'(r)|\,|x_n-r|. \]If
\[ |g'(r)|=0.2, \]then the error may shrink by roughly a factor of \(0.2\) per iteration near the fixed point.
If
\[ |g'(r)|=0.9, \]convergence may be much slower.
Choosing a Starting Value
Fixed-point iteration requires an initial guess
\[ x_0. \]Useful starting values may come from:
- a graph;
- a numerical table;
- a bracketing interval;
- physical constraints;
- previous calculations;
- simple nearby values.
Initial Guess and Basin of Attraction
A fixed point may attract some starting values but not others.
The collection of starting values that converge to a particular fixed point is called its basin of attraction.
In elementary work, the practical lesson is simple: a reasonable starting estimate matters.
Example: Solving \(x=\cos x\)
Consider
\[ x=\cos x. \]The equation is already in fixed-point form with
\[ g(x)=\cos x. \]Use the iteration
\[ \boxed{x_{n+1}=\cos x_n}. \]Beginning the Cosine Iteration
Let
\[ x_0=1. \]Then
\[ x_1=\cos(1)\approx0.540302, \] \[ x_2=\cos(0.540302)\approx0.857553, \] \[ x_3\approx0.654290, \] \[ x_4\approx0.793480. \]The values oscillate while gradually narrowing around the fixed point.
Continuing the Cosine Iteration
Further iterations give
\[ x_{10}\approx0.744237, \] \[ x_{20}\approx0.739184, \] \[ x_{30}\approx0.739087. \]The sequence approaches
\[ \boxed{x\approx0.739085}. \]Why the Cosine Iteration Converges
Here,
\[ g'(x)=-\sin x. \]At the fixed point \(r\approx0.739085\),
\[ |g'(r)| = |\sin r| \approx0.674. \]The magnitude is less than \(1\), which is consistent with local convergence.
Residual Checking
If \(x_n\) is intended to solve
\[ f(x)=0, \]evaluate
\[ f(x_n). \]A small value of
\[ |f(x_n)| \]supports the claim that \(x_n\) is close to satisfying the original equation.
Fixed-Point Residual
For the fixed-point equation
\[ x=g(x), \]one natural residual is
\[ \boxed{g(x_n)-x_n}. \]At an exact fixed point,
\[ g(r)-r=0. \]Stopping by Successive Difference
A common stopping condition is
\[ \boxed{|x_{n+1}-x_n|<\varepsilon}. \]This checks whether the sequence is changing by less than a specified tolerance.
It is easy to compute because both values are already available.
Stopping by Residual
Another condition is
\[ \boxed{|g(x_n)-x_n|<\varepsilon}. \]If solving an original equation \(f(x)=0\), we may instead require
\[ \boxed{|f(x_n)|<\varepsilon}. \]These checks ask whether the current estimate nearly satisfies the target equation.
Using More Than One Stopping Condition
A robust implementation may require both
\[ |x_{n+1}-x_n|<\varepsilon \]and
\[ |f(x_{n+1})|<\varepsilon. \]This checks both stability of the iteration and satisfaction of the equation.
Maximum Iteration Safeguard
Fixed-point iteration is not guaranteed to converge for every rule and starting value.
Therefore, computational implementations should also include a maximum iteration count such as
\[ n\le100. \]If convergence has not occurred by then, the algorithm can stop and report failure to meet the tolerance.
Fixed-Point Iteration in Pseudocode
set x = initial guess
repeat
new_x = g(x)
if stopping condition is satisfied
return new_x
x = new_x
until maximum iterations are reached
This simple loop captures the basic structure of the method.
Adding a Residual Check
set x = initial guess
repeat
new_x = g(x)
if abs(new_x - x) < tolerance
and abs(f(new_x)) < tolerance
return new_x
x = new_x
until maximum iterations are reached
The exact stopping logic depends on the problem and the error guarantees required.
Fixed-Point Iteration versus Bisection
Bisection maintains a bracket:
\[ [a_n,b_n]. \]Fixed-point iteration usually maintains a single current estimate:
\[ x_n. \]Bisection offers strong interval guarantees under its assumptions.
Fixed-point iteration can be faster or simpler, but its convergence depends more strongly on the chosen update rule.
Fixed-Point Iteration versus Direct Algebra
If an equation has a simple exact algebraic solution, fixed-point iteration may be unnecessary.
For example,
\[ 2x+3=11 \]is solved immediately by
\[ x=4. \]Numerical iteration is most useful when direct algebra is unavailable, inconvenient, or embedded inside a larger computational process.
Fixed-Point Iteration as a Building Block
Many important numerical methods can be interpreted as specialized fixed-point iterations.
Newton's method, for example, updates by
\[ x_{n+1} = x_n-\frac{f(x_n)}{f'(x_n)}. \]This can be viewed as
\[ x_{n+1}=g(x_n) \]for a particular choice of \(g\).
Common Error: Assuming Any Rearrangement Will Work
Algebraic equivalence does not guarantee numerical convergence.
Two formulas may have the same fixed point but very different iterative behavior.
The rule
\[ x_{n+1}=\frac3{x_n} \]illustrates this danger for \(x^2=3\).
Common Error: Ignoring the Starting Value
A convergent rule near one fixed point may fail from a distant starting value.
Always consider whether the initial guess lies in a region where the update rule is defined and expected to converge.
Common Error: Dividing by Zero
Some iteration formulas contain denominators.
For example,
\[ x_{n+1} = \frac12\left(x_n+\frac3{x_n}\right) \]is undefined when
\[ x_n=0. \]Domain restrictions must be checked before and during iteration.
Common Error: Treating a Repeating Cycle as Convergence
The sequence
\[ 1,\quad3,\quad1,\quad3,\quad\ldots \]does not converge merely because it repeats predictably.
Convergence requires the sequence to approach a single limiting value.
Common Error: Stopping Because Several Digits Look Similar
Visual similarity is not a mathematical stopping rule.
Use a criterion such as
\[ |x_{n+1}-x_n|<10^{-6} \]or
\[ |f(x_n)|<10^{-6}. \]Common Error: Reporting More Precision Than the Iteration Supports
If the stopping tolerance controls only about four decimal places, reporting twelve digits suggests unjustified accuracy.
The reported precision should match the numerical evidence.
Common Error: Forgetting the Original Equation
The purpose of fixed-point iteration is to solve the original problem.
After convergence, substitute the approximation into the original equation or compute its residual.
This guards against algebraic rearrangement errors and unintended fixed points.
A Reliable Fixed-Point Workflow
- Write the target equation clearly.
- Rearrange it into \(x=g(x)\).
- Check domain restrictions introduced by the rearrangement.
- Choose a reasonable initial guess \(x_0\).
- Examine \(g'(x)\) near the expected fixed point when possible.
- Compute \(x_{n+1}=g(x_n)\).
- Record the successive estimates.
- Apply a stated stopping condition.
- Impose a maximum iteration count.
- Verify the final approximation in the original equation.
Example: Full Workflow for \(x=\cos x\)
Solve
\[ x=\cos x. \]The equation already has the form
\[ x=g(x) \]with
\[ g(x)=\cos x. \]Choose
\[ x_0=1. \]Iterate:
\[ x_{n+1}=\cos x_n. \]Continue until
\[ |x_{n+1}-x_n|<10^{-6}. \]Verifying the Cosine Solution
The iteration gives approximately
\[ x\approx0.739085. \]Check:
\[ \cos(0.739085)\approx0.739085. \]Equivalently, for
\[ f(x)=\cos x-x, \]the residual is close to zero.
Practice 1: Identify the Fixed Point
Find the fixed point of
\[ g(x)=\frac{x+8}{2}. \]Practice 2: Generate Iterates
Let
\[ x_{n+1}=\frac{x_n+8}{2}, \qquad x_0=0. \]Compute \(x_1,x_2,x_3,x_4\).
Practice 3: Rearrange into Fixed-Point Form
Rewrite
\[ x^3+x-2=0 \]in the form
\[ x=g(x) \]by isolating the linear \(x\)-term.
Practice 4: Detect a Cycle
Use
\[ x_{n+1}=\frac4{x_n}, \qquad x_0=1. \]Compute the first four iterates and describe the behavior.
Practice 5: Derivative Test
A fixed point \(r\) satisfies
\[ |g'(r)|=0.35. \]What does this suggest about local convergence?
Practice 6: Repelling Fixed Point
A fixed point satisfies
\[ |g'(r)|=1.4. \]What does this suggest about nearby iteration?
Practice 7: Stopping Condition
Suppose
\[ x_n=2.718281 \]and
\[ x_{n+1}=2.718282. \]Does the condition
\[ |x_{n+1}-x_n|<10^{-5} \]hold?
Practice 8: Fixed-Point Residual
If
\[ g(x)=\cos x \]and
\[ x=0.74, \]write the fixed-point residual.
Practice 9: Domain Check
Why is
\[ x_0=0 \]invalid for the iteration
\[ x_{n+1} = \frac12\left(x_n+\frac5{x_n}\right)? \]Practice 10: Compare Methods
State one major advantage of bisection and one major advantage of fixed-point iteration.
Practice Answers
Practice 1
Solve
\[ x=\frac{x+8}{2}. \]Then
\[ 2x=x+8, \] \[ \boxed{x=8}. \]Practice 2
\[ x_1=\frac{0+8}{2}=4, \] \[ x_2=\frac{4+8}{2}=6, \] \[ x_3=\frac{6+8}{2}=7, \] \[ x_4=\frac{7+8}{2}=7.5. \]Therefore,
\[ \boxed{ x_1=4,\quad x_2=6,\quad x_3=7,\quad x_4=7.5 }. \]Practice 3
From
\[ x^3+x-2=0, \]isolate \(x\):
\[ x=2-x^3. \]Therefore, one fixed-point form is
\[ \boxed{x=2-x^3}. \]Whether the associated iteration converges is a separate numerical question.
Practice 4
\[ x_1=4, \] \[ x_2=1, \] \[ x_3=4, \] \[ x_4=1. \]The sequence cycles between \(1\) and \(4\) rather than converging to the positive solution \(2\).
Practice 5
Since
\[ |g'(r)|=0.35<1, \]the fixed point is locally attracting, so nearby starting values are expected to converge under appropriate smoothness conditions.
Practice 6
Since
\[ |g'(r)|=1.4>1, \]the fixed point is locally repelling, so nearby errors tend to grow.
Practice 7
\[ |2.718282-2.718281| = 0.000001 = 10^{-6}. \]Since
\[ 10^{-6}<10^{-5}, \]the condition is satisfied.
Practice 8
The fixed-point residual is
\[ \boxed{\cos(0.74)-0.74}. \]Practice 9
The iteration contains the term
\[ \frac5{x_n}. \]If
\[ x_0=0, \]the first update requires division by zero, so the iteration is undefined.
Practice 10
Bisection provides a strong bracket and predictable error reduction when its conditions hold.
Fixed-point iteration can use a very simple one-value update and may converge efficiently when a good iteration function is chosen.
What to Remember
- A fixed point satisfies \(g(r)=r\).
- Fixed-point iteration uses \(x_{n+1}=g(x_n)\).
- Solutions of an equation can often be found by rewriting it as \(x=g(x)\).
- Different rearrangements of the same equation can behave very differently.
- Local convergence is often associated with \(|g'(r)|<1\).
- If \(|g'(r)|>1\), the fixed point is typically locally repelling.
- The case \(|g'(r)|=1\) requires additional analysis.
- Starting values and domain restrictions matter.
- Stopping conditions should use numerical criteria, not visual judgment.
- The final approximation should be checked against the original equation.
Where We Go Next
Fixed-point iteration shows how a carefully chosen update rule can move a single estimate toward a solution. The next method uses local slope information to construct a particularly powerful update rule.
The next lesson, Newton's Method, uses tangent lines and derivatives to generate approximations that often converge much faster than basic fixed-point iteration.