Lecture 4 Fixed Point Method
Lecture 4 Fixed Point Method
Advantages:
Useful for large, complex problems where direct
methods.
Disadvantages:
Convergence is not always guaranteed.
x 0 1 2
f(x) -1 -1 5
Example 1
x0=(1+2)/2=1.5
x1=ϕ(x0)=ϕ(1.5)=1.35721
x2=ϕ(x1)=ϕ(1.35721)=1.33086
x3=ϕ(x2)=ϕ(1.33086)=1.32588
x4=ϕ(x3)=ϕ(1.32588)=1.32494
x5=ϕ(x4)=ϕ(1.32494)=1.32476
Table
Difference
n x0 x1=ϕ(x0) Update
|x1-x0|
2 1.5 1.35721 x0=x1 0.14279
3 1.35721 1.33086 x0=x1 0.02635
4 1.33086 1.32588 x0=x1 0.00498
5 1.32588 1.32494 x0=x1 0.00094
6 1.32494 1.32476 x0=x1 0.00018
Example 2
Find a root of an equation using Fixed
Point Iteration method.
Solution:
x1=ϕ(x0)=ϕ(0.5)=0.44444
x2=ϕ(x1)=ϕ(0.44444)=0.47929
x3=ϕ(x2)=ϕ(0.47929)=0.45698
x4=ϕ(x3)=ϕ(0.45698)=0.47108
x5=ϕ(x4)=ϕ(0.47108)=0.46209
x6=ϕ(x5)=ϕ(0.46209)=0.46779
x7=ϕ(x6)=ϕ(0.46779)=0.46416
x8=ϕ(x7)=ϕ(0.46416)=0.46647
x9=ϕ(x8)=ϕ(0.46647)=0.465
x10=ϕ(x9)=ϕ(0.465)=0.46593
x11=ϕ(x10)=ϕ(0.46593)=0.46534
x12=ϕ(x11)=ϕ(0.46534)=0.46572
Difference
n x0 x1=ϕ(x0) Update
|x1-x0|
2 0.5 0.44444 x0=x1 0.05556
3 0.44444 0.47929 x0=x1 0.03485
4 0.47929 0.45698 x0=x1 0.02231
5 0.45698 0.47108 x0=x1 0.0141
6 0.47108 0.46209 x0=x1 0.00899
7 0.46209 0.46779 x0=x1 0.0057
8 0.46779 0.46416 x0=x1 0.00363
9 0.46416 0.46647 x0=x1 0.0023
10 0.46647 0.465 x0=x1 0.00146
11 0.465 0.46593 x0=x1 0.00093
12 0.46593 0.46534 x0=x1 0.00059
13 0.46534 0.46572 x0=x1 0.00038
Questions
Find a root of an equation using Fixed Point
Iteration method
MATLAB CODE
% iterative method
f=@(x)x^2 - sin(x) - 0.5; % main
function
xf = @(x) (sin(x)+ 0.5)^(1/2); %
Derivative of function
% Initial guess
a = 0.5;
for i = 1:15
iteration = i+1 % iteration number
b = xf(a)
a = b; % replace b by a
end