First Semester in Numerical Analysis with Julia, 2020a
First Semester in Numerical Analysis with Julia, 2020a
First Semester in Numerical Analysis with Julia, 2020a
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
CHAPTER 2. SOLUTIONS OF EQUATIONS: ROOT-FINDING 77<br />
2.6 Fixed-po<strong>in</strong>t iteration<br />
Many root-f<strong>in</strong>d<strong>in</strong>g methods are based on the so-called fixed-po<strong>in</strong>t iteration; a method we<br />
discuss <strong>in</strong> this section.<br />
Def<strong>in</strong>ition 37. Anumberp is a fixed-po<strong>in</strong>t for a function g(x) if g(p) =p.<br />
We have two problems that are related to each other:<br />
• Fixed-po<strong>in</strong>t problem: F<strong>in</strong>d p such that g(p) =p.<br />
• Root-f<strong>in</strong>d<strong>in</strong>g problem: F<strong>in</strong>d p such that f(p) =0.<br />
We can formulate a root-f<strong>in</strong>d<strong>in</strong>g problem as a fixed-po<strong>in</strong>t problem, and vice versa. For<br />
example, assume we want to solve the root f<strong>in</strong>d<strong>in</strong>g problem, f(p) =0. Def<strong>in</strong>e g(x) =x−f(x),<br />
and observe that if p is a fixed-po<strong>in</strong>t of g(x), that is, g(p) =p − f(p) =p, then p is a root<br />
of f(x). Here the function g is not unique: there are many ways one can represent the<br />
root-f<strong>in</strong>d<strong>in</strong>g problem f(p) =0as a fixed-po<strong>in</strong>t problem, and as we will learn later, not all<br />
will be useful to us <strong>in</strong> develop<strong>in</strong>g fixed-po<strong>in</strong>t iteration algorithms.<br />
The next theorem answers the follow<strong>in</strong>g questions: When does a function g have a fixedpo<strong>in</strong>t?<br />
If it has a fixed-po<strong>in</strong>t, is it unique?<br />
Theorem 38. 1. If g is a cont<strong>in</strong>uous function on [a, b] and g(x) ∈ [a, b] for all x ∈ [a, b],<br />
then g has at least one fixed-po<strong>in</strong>t <strong>in</strong> [a, b].<br />
2. If, <strong>in</strong> addition, |g(x) − g(y)| ≤λ|x − y| for all x, y ∈ [a, b] where 0