diff options
Diffstat (limited to '845/CH2/EX2.12')
-rwxr-xr-x | 845/CH2/EX2.12/Ex2_12.sce | 19 |
1 files changed, 19 insertions, 0 deletions
diff --git a/845/CH2/EX2.12/Ex2_12.sce b/845/CH2/EX2.12/Ex2_12.sce new file mode 100755 index 000000000..7ff416a86 --- /dev/null +++ b/845/CH2/EX2.12/Ex2_12.sce @@ -0,0 +1,19 @@ +//Example 2.12
+clc
+clear
+
+function [f] = fun12(x)
+ f = x.^3 - x - 1;
+endfunction
+
+x = [0 1 2];
+h = [x(2)-x(1) x(3)-x(2)];
+lamdai = h(2)/h(1);
+deli = 1 + lamdai;
+f = fun12(x);
+
+g = f(1)*lamdai^2 - f(2)*deli^2 + f(3)*(lamdai + deli);
+lamda = -2*f(3)*deli / (g + sqrt(g^2 - 4*f(3)*deli*(f(1)*lamdai - f(2)*deli + f(3))));
+xnew = x(3) + lamda*h(2);
+xnew = round(xnew*10^5) / 10^5;
+disp(xnew,"root = ")
|