blob: ac713b96e321076cc6b6b2eebd2a9b06b1bb5a88 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
|
clc;
clear;
function y=f(x)
y=x^3 - 13*x - 12
endfunction
x1t=-3;
x2t=-1;
x3t=4;
x0=4.5;
x1=5.5;
x2=5;
disp(0,"iteration:")
disp(x2,"xr:")
disp("---------------------------------------------")
for i=1:4
h0=x1-x0;
h1=x2-x1;
d0=(f(x1)-f(x0))/(x1-x0);
d1=(f(x2)-f(x1))/(x2-x1);
a=(d1-d0)/(h1+h0);
b=a*h1+d1;
c=f(x2);
d=(b^2 - 4*a*c)^0.5;
if abs(b+d)>abs(b-d) then x3=x2+((-2*c)/(b+d));
else x3=x2+((-2*c)/(b-d)); end
ea=abs(x3-x2)*100/x3;
x0=x1;
x1=x2;
x2=x3;
disp(i,"iteration:")
disp(x2,"xr:")
disp(ea,"ea(%):")
disp("---------------------------------------------")
end
|
