413/CH1/EX1.2/Table2.sce


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 secant 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;// second vale
x2=4; //just to start loop
tol=0.0000001; // tolerence value
x = x0;   fa=eval(fx);
x = x1;   fb=eval(fx);
x = x2;   fc=eval(fx);
i=1;
while abs(fc)>tol
    if abs(fa)<abs(fb) then
    xc=x0;
    x0=x1;
    x1=xc;  
    end
    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)
    x0=x1;
    x1=x2;
    i=i+1;
end;