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/CH10/EX10.5/Ex10_5.sce | 60 ++++++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 60 insertions(+)
 create mode 100755 845/CH10/EX10.5/Ex10_5.sce

(limited to '845/CH10/EX10.5/Ex10_5.sce')

diff --git a/845/CH10/EX10.5/Ex10_5.sce b/845/CH10/EX10.5/Ex10_5.sce
new file mode 100755
index 000000000..889edee01
--- /dev/null
+++ b/845/CH10/EX10.5/Ex10_5.sce
@@ -0,0 +1,60 @@
+//Example 10.5
+
+clc
+clear
+
+x = 0:4;
+y = 0:4;
+m = length(y);
+n = length(x);
+u = zeros(m,n);
+u(m,:) = x.^3;
+u(:,n) = 16*y';
+u0 = u;
+
+I = ceil(m/2);
+J = ceil(n/2);
+
+u(J,I) = (u0(J-2,I-2) + u0(J-2,I+2) + u0(J+2,I-2) + u0(J+2,I+2)) / 4;
+
+for j = [J-1 J+1]
+    for i = [I-1 I+1]
+        u(j,i) = (u(j-1,i-1) + u(j-1,i+1) + u(j+1,i-1) + u(j+1,i+1)) / 4;
+    end
+end
+
+j1 = [J-1 J J J+1];
+i1 = [I I-1 I+1 I];
+for k = 1:4
+    i = i1(k);
+    j = j1(k);
+    u(j,i) = (u(j,i-1) + u(j,i+1) + u(j-1,i) + u(j+1,i)) / 4;
+end
+disp(u,"u:")
+
+p = m-1;
+q = n-1;
+c = cos(%pi/p) + cos(%pi/q);
+w = 4/(2+sqrt(4-c^2));
+w = round(w*10^3)/10^3;
+
+kf = 10;
+tab = zeros(kf+1,(m-2)*(n-2));
+row = [];
+for j = 2:n-1
+    row = [row u(j,2:m-1)];
+end
+tab(1,:) = row;
+for k = 1:kf
+    row = [];
+    for j = 2:n-1
+        for i = 2:m-1
+            u(j,i) = (u(j,i-1) + u(j,i+1) + u(j-1,i) + u(j+1,i)) *w/4 + (1-w)*u(j,i);
+        end
+        row = [row u(j,2:m-1)];
+    end
+    row = round(row*10^4)/10^4;
+    tab(k+1,:) = row;
+end
+mprintf("\n\n%8s %9s %10s %10s %9s %10s %10s %9s %9s",'u11','u21','u31','u12','u22','u32','u13','u23','u33')
+disp(tab)
-- 
cgit