From b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b Mon Sep 17 00:00:00 2001
From: priyanka
Date: Wed, 24 Jun 2015 15:03:17 +0530
Subject: initial commit / add all books

---
 845/CH2/EX2.4/Ex2_4.sce | 57 +++++++++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 57 insertions(+)
 create mode 100755 845/CH2/EX2.4/Ex2_4.sce

(limited to '845/CH2/EX2.4/Ex2_4.sce')

diff --git a/845/CH2/EX2.4/Ex2_4.sce b/845/CH2/EX2.4/Ex2_4.sce
new file mode 100755
index 000000000..da626dcae
--- /dev/null
+++ b/845/CH2/EX2.4/Ex2_4.sce
@@ -0,0 +1,57 @@
+//Example 2.4
+clc
+clear
+
+function [root] = FalsePosition(fun,x,tol,maxit)
+// FalsePosition: Computes roots of the function in the given range using False Position Method
+//// Input: FalsePosition(fun,x,tol,maxit)
+// fun = function handle
+// x = range in between sign change is evident
+// tol = Maximum error between iterations that can be tolerated
+// maxit = Maximum number of iterations
+//// Output: [root]
+// Root: Root of the given function in defined range
+
+if fun(x(1)) > 0 then
+    xu = x(1);    xl = x(2);
+else
+    xu = x(2);    xl = x(1);
+end
+
+Ea = 1;
+iter = 1;
+
+while(1)
+    xr(iter) = xl(iter) - ((xu(iter)-xl(iter)) / (fun(xu(iter))-fun(xl(iter))) * fun(xl(iter)));
+    if fun(xr(iter)) > 0 then
+        xu(iter+1) = xr(iter);
+        xl(iter+1) = xl(iter);
+    elseif fun(xr(iter)) < 0 then
+        xl(iter+1) = xr(iter);
+        xu(iter+1) = xu(iter);
+    else
+        break
+    end
+    
+    if iter>1 then
+        Ea(iter) = 100 * abs((xr(iter) - xr(iter-1)) / xr(iter));
+    end
+   
+    if Ea(iter) < tol | iter == maxit then
+        break
+    end
+    iter = iter + 1;
+end
+root = xr(iter);
+endfunction
+
+function f = fun4(x)
+    f = x.*log10(x) - 1.2;
+endfunction
+
+clc
+x = [2 3];
+tol = 1e-4;
+maxit = 2;
+root = FalsePosition(fun4,x,tol,maxit);
+disp(round(root*10^4)/10^4,"root = ")
-- 
cgit