1 files changed, 35 insertions, 0 deletions
diff --git a/191/CH3/EX3.1/Example3_1.sce b/191/CH3/EX3.1/Example3_1.sce
new file mode 100755
index 000000000..5a09ca883
--- /dev/null
+++ b/191/CH3/EX3.1/Example3_1.sce
@@ -0,0 +1,35 @@
+//Bisection Method
+clc;
+clear;
+close();
+format('v',9);
+b(1)=1;a(1)=0;k=5;
+deff('[fx]=bisec(x)','fx =(x+1).^2.*exp(x.^2-2)-1');
+x = linspace(0,1);
+plot(x,((x+1).^2).*(exp(x.^2-2))-1);
+//in interval [0,1]changes its sign thus has a root
+//k = no of decimal place of accuracy
+//a = lower limit of interval
+//b = upper limit of interval
+//n = no of iterations required
+n = log2((10^k)*(b-a));
+n = ceil(n);
+disp(n,'Number of iterations : ')
+for i = 1:n-1
+    N(i) = i;
+    c(i) = (a(i)+b(i))/2;
+    bs(i) = bisec(c(i));
+    if (bisec(b(i))*bisec(c(i))<0)
+        a(i+1)=c(i);
+        b(i+1)=b(i);
+    else
+        b(i+1)=c(i);
+        a(i+1)=a(i);
+    end
+end
+N(i+1)=i+1;
+c(i+1) = (a(i+1)+b(i+1))/2;
+bs(i+1) = bisec(c(i));
+ann = [N a b c bs];
+disp(ann , 'The Table : ');
+disp(c(i),'The root of the function is : ')
+\ No newline at end of file