diff options
author | priyanka | 2015-06-24 15:03:17 +0530 |
---|---|---|
committer | priyanka | 2015-06-24 15:03:17 +0530 |
commit | b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (patch) | |
tree | ab291cffc65280e58ac82470ba63fbcca7805165 /845/CH10/EX10.4 | |
download | Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.gz Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.bz2 Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.zip |
initial commit / add all books
Diffstat (limited to '845/CH10/EX10.4')
-rwxr-xr-x | 845/CH10/EX10.4/Ex10_4.sce | 29 |
1 files changed, 29 insertions, 0 deletions
diff --git a/845/CH10/EX10.4/Ex10_4.sce b/845/CH10/EX10.4/Ex10_4.sce new file mode 100755 index 000000000..4840d8c21 --- /dev/null +++ b/845/CH10/EX10.4/Ex10_4.sce @@ -0,0 +1,29 @@ +//Example 10.4
+
+clc
+clear
+
+h = 1/3;
+x = 0:h:1;
+y = 0:h:1;
+m = length(y);
+n = length(x);
+u = zeros(m,n);
+u(m,2:n-1) = 1;
+
+kf = 5;
+tab = zeros(kf,(m-2)*(n-2));
+for k = 1:kf
+ row = [];
+ for j = 2:n-1
+ for i = 2:m-1
+ constant = 10/9* (5 + 1/9*(i-1)^2 + 1/9*(j-1)^2);
+ u(j,i) = (u(j,i-1) + u(j,i+1) + u(j-1,i) + u(j+1,i) + constant) / 4;
+ end
+ row = [row u(j,2:m-1)];
+ end
+ row = round(row*10^4)/10^4;
+ tab(k,:) = row;
+end
+mprintf("%4s %9s %9s %9s %9s",'r','u11','u21','u12','u22')
+disp([(1:k)' tab])
|