blob: 3010f55efd70b80132a31347e0146a8cf3f923dd (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
|
//Example 5.2
//Bisection Method
//Page no. 147
clc;clear;close;
deff('y=f(x)','y=x^x-2*x+2')
x1=0;x2=2;e=0.001;i=0;
printf('Iteration\tx1\t\tx2\t\tz\t\tf(z)\n')
printf('--------------------------------------------------------------------------\n')
while abs(x1-x2)>e
z=(x1+x2)/2
printf(' %i\t\t%f\t%f\t%f\t%f\n',i,x1,x2,z,f(z))
if f(z)*f(x1)>0
x1=z
else
x2=z
end
i=i+1
end
printf('\n\nThe solution of this equation is %g after %i Iterations',z,i-1)
printf('\n\n\nNote : There are computational errors in the answer given by the book for this example')
|
