06.09.2021 Views

First Semester in Numerical Analysis with Julia, 2020a

First Semester in Numerical Analysis with Julia, 2020a

First Semester in Numerical Analysis with Julia, 2020a

SHOW MORE
SHOW LESS

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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!