1 files changed, 59 insertions, 0 deletions
diff --git a/260/CH11/EX11.8/11_8.sce b/260/CH11/EX11.8/11_8.sce
new file mode 100644
index 000000000..27febdf1d
--- /dev/null
+++ b/260/CH11/EX11.8/11_8.sce
@@ -0,0 +1,59 @@
+//Eg-11.8
+//pg-485
+
+clear
+clc
+
+a = 0;
+b = 1;
+h = b-a;
+
+deff('out = func(in)','out = 1/(1+in^2)')
+
+//From equations [30],[31],[32] & [33]
+
+//Please note that the subscripts(i&j) we use here are different from that used in 
+//text book i.e they are increased by 1, because we cant give the index zero in   //scilab. Therefore,
+
+imax = 6;
+jmax = 6;
+
+I(1,1) = h/2*(func(a) + func(b));
+
+I(2,1) = 1/2*(I(1,1) + h*func(a+h/2));
+
+I(3,1) = 1/2*(I(2,1) + h/2*(func(a+h/4) + func(a+3*h/4)));
+
+//From equation [33]
+
+sum1 = 0;
+
+for(j = 1:2:(2^3-1))    //Since we have to consider the odd terms only.
+    sum1 = sum1 + func(a+j*h/2^3);
+end
+
+
+I(4,1) = 1/2*(I(3,1) + h/2^2*sum1);
+
+//Similarly
+
+sum2 = 0;
+
+for(j = 1:2:(2^4-1))
+    sum2 = sum2 + func(a+j*h/2^4);
+end
+
+I(5,1) = 1/2*(I(4,1) + h/2^3*sum2);
+
+for(j = 2:5)
+    for(i = 1:imax-j)
+        I(i,j) = (4^(j-1)*I(i+1,j-1) - I(i,j-1))/(4^(j-1)-1);
+    end
+end
+
+printf(' The complete Romberg tableau is as follows\n')
+
+disp(I)
+
+printf('\n  Therefore, the value of the integral is %f\n',I(1,5))
+