1 files changed, 33 insertions, 0 deletions
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