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

---
 1670/CH8/EX8.9/8_9.sce | 33 +++++++++++++++++++++++++++++++++
 1 file changed, 33 insertions(+)
 create mode 100755 1670/CH8/EX8.9/8_9.sce

(limited to '1670/CH8/EX8.9/8_9.sce')

diff --git a/1670/CH8/EX8.9/8_9.sce b/1670/CH8/EX8.9/8_9.sce
new file mode 100755
index 000000000..f1387f600
--- /dev/null
+++ b/1670/CH8/EX8.9/8_9.sce
@@ -0,0 +1,33 @@
+//Example 8.9
+//Rombers Method
+//Page no 275
+clc;clear;close;
+
+a=4;b=5.2;
+h=[0.4,0.2];
+for j=1:2
+    n=(b-a)/h(j)+1
+for i=1:n
+    x(1,i)=a+(i-1)*h(j)
+    y(1,i)=log(x(1,i))
+end
+Q(j)=0;
+for i=1:n
+    if i==1 | i==n then
+        Q(j)=Q(j)+h(j)*(y(1,i))/2
+    else
+        Q(j)=Q(j)+h(j)*(y(1,i))
+    end
+end
+printf('\nx :')
+for k=1:n
+    printf('\t %g',x(1,k))
+end
+printf('\nf(x) :')
+for k=1:n
+    printf('\t%.4f',y(1,k))
+end
+printf('\n\nQ(%i) = %g\n\n\n',j,Q(j))
+end
+R1=4*Q(2)/3-Q(1)/3
+printf('R1 = %g',R1)
\ No newline at end of file
-- 
cgit