blob: 549db769f6ef77a50f8bc072ab21426cceb8676e (
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
|
//The regula falsi method for f(x)=3*x+sin(x)-exp(x), starting from 0 and 1
clearglobal()
clc;
fx='3*x+sin(x)-exp(x)' //Define function here
x0=0; // intial value
x1=1;// final vale where root need to bracket
x2=4;//just to start loop
x = x0; fa=eval(fx);
x = x1; fb=eval(fx);
tol=0.0000001;
x = x2; fc=eval(fx);
i=1;
if fa*fb<0 then
while abs(fc)>tol
x = x0; fa=eval(fx);
x = x1; fb=eval(fx);
x2=x1-fb*(x0-x1)/(fa-fb);
x = x2; fc=eval(fx);
X = [i,x0,x1,x2,fc];
disp(X)
if fc*fa < 0 then
x1 = x2;
else x0 = x2;
end;
i=i+1;
end;
end
|
