diff options
Diffstat (limited to '413/CH1/EX1.3/Table_3.sce')
-rw-r--r-- | 413/CH1/EX1.3/Table_3.sce | 28 |
1 files changed, 28 insertions, 0 deletions
diff --git a/413/CH1/EX1.3/Table_3.sce b/413/CH1/EX1.3/Table_3.sce new file mode 100644 index 000000000..549db769f --- /dev/null +++ b/413/CH1/EX1.3/Table_3.sce @@ -0,0 +1,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
\ No newline at end of file |