initial commit / add all books

6 files changed, 114 insertions, 0 deletions

diff --git a/1673/CH9/EX9.1/9_1.sce b/1673/CH9/EX9.1/9_1.sce
new file mode 100755
index 000000000..b0d6ed66e
--- /dev/null
+++ b/1673/CH9/EX9.1/9_1.sce
@@ -0,0 +1,11 @@
+//standard five point formula

+//example 9.1

+//page 350

+clc;clear;close;

+u2=5;u3=1;

+for i=1:3

+    u1=(u2+u3+6)/4;

+    u2=u1/2+5/2;

+    u3=u1/2+1/2;

+    printf(' the values are u1=%d\t u2=%d\t u3=%d\t\n\n',u1,u2,u3);

+end
\ No newline at end of file
diff --git a/1673/CH9/EX9.2/9_2.sce b/1673/CH9/EX9.2/9_2.sce
new file mode 100755
index 000000000..db3e3a824
--- /dev/null
+++ b/1673/CH9/EX9.2/9_2.sce
@@ -0,0 +1,29 @@
+//solution of laplace equation by jacobi method,gauss-seidel method and SOR method

+//example 9.2

+//page 351

+clc;clear;close;

+u1=0.25;u2=0.25;u3=0.5;u4=0.5;//initial values

+printf('jacobis iteration process\n\n')

+printf('u1\t     u2\t     u3\t     u4\t \n\n')

+printf('%f\t     %f\t     %f\t     %f\t  \n',u1,u2,u3,u4)

+for i=1:7

+    u11=(0+u2+0+u4)/4

+    u22=(u1+0+0+u3)/4;

+    u33=(1+u2+0+u4)/4;

+    u44=(1+0+u3+u1)/4;

+    u1=u11;u2=u22;u3=u33;u4=u44;

+printf('%f\t     %f\t     %f\t     %f\t  \n',u11,u22,u33,u44)    

+end

+printf(' gauss seidel process\n\n');

+u1=0.25;u2=0.3125;u3=0.5625;u4=0.46875;//initial values

+printf('u1\t     u2\t     u3\t     u4\t \n\n')

+printf('%f\t     %f\t     %f\t     %f\t  \n',u1,u2,u3,u4)

+for i=1:4

+    u1=(0+u2+0+u4)/4

+    u2=(u1+0+0+u3)/4;

+    u3=(1+u2+0+u4)/4;

+    u4=(1+0+u3+u1)/4;

+    printf('%f\t     %f\t     %f\t     %f\t  \n',u1,u2,u3,u4)    

+end

+printf('u1\t     u2\t     u3\t     u4\t \n\n')

+printf('%f\t     %f\t     %f\t     %f\t  \n',u1,u2,u3,u4)

diff --git a/1673/CH9/EX9.4/9_4.sce b/1673/CH9/EX9.4/9_4.sce
new file mode 100755
index 000000000..989e9cde8
--- /dev/null
+++ b/1673/CH9/EX9.4/9_4.sce
@@ -0,0 +1,13 @@
+//poisson equation

+//exaample 9.4

+//page 354

+clc;clear;close;

+u2=0;u4=0;

+printf('  u1\t    u2\t   u3\t   u4\t\n\n');

+for i=1:6

+    u1=u2/2+30;

+    u2=(u1+u4+150)/4;

+    u4=u2/2+45;

+    printf(' %0.2f\t   %0.2f\t   %0.2f\t   %0.2f\n',u1,u2,u2,u4);

+end

+printf(' from last two iterates we conclude u1=67   u2=75   u3=75   u4=83\n')
\ No newline at end of file
diff --git a/1673/CH9/EX9.6/9_6.sce b/1673/CH9/EX9.6/9_6.sce
new file mode 100755
index 000000000..42d52bd1f
--- /dev/null
+++ b/1673/CH9/EX9.6/9_6.sce
@@ -0,0 +1,26 @@
+//bender-schmidt formula

+//example 9.6

+//page 362

+clc;clear;close;

+deff('y=f(x)','y=4*x-x^2');

+u=[f(0) f(1) f(2) f(3) f(4)];

+u11=(u(1)+u(3))/2;

+u12=(u(2)+u(4))/2;

+u13=(u(3)+u(5))/2;

+printf(' u11=%0.2f\t  u12=%0.2f\t  u13=%0.2f\t \n',u11,u12,u13)

+u21=(u(1)+u12)/2;

+u22=(u11+u13)/2;

+u23=(u12+0)/2;

+printf(' u21=%0.2f\t  u22=%0.2f\t  u23=%0.2f\t \n',u21,u22,u23)

+u31=(u(1)+u22)/2;

+u32=(u21+u23)/2;

+u33=(u22+u(1))/2;

+printf(' u31=%0.2f\t  u32=%0.2f\t  u33=%0.2f\t \n',u31,u32,u33)

+u41=(u(1)+u32)/2;

+u42=(u31+u33)/2;

+u43=(u32+u(1))/2;

+printf(' u41=%0.2f\t  u42=%0.2f\t  u43=%0.2f\t \n',u41,u42,u43)

+u51=(u(1)+u42)/2;

+u52=(u41+u43)/2;

+u53=(u42+u(1))/2;

+printf(' u51=%0.2f\t  u52=%0.2f\t  u53=%0.2f\t \n',u51,u52,u53)
\ No newline at end of file
diff --git a/1673/CH9/EX9.7/9_7.sce b/1673/CH9/EX9.7/9_7.sce
new file mode 100755
index 000000000..f259cae3c
--- /dev/null
+++ b/1673/CH9/EX9.7/9_7.sce
@@ -0,0 +1,19 @@
+//bender-schimdt's formula and crank-nicolson formula

+//example 9.7

+//page 363

+//bender -schimdt's formula

+clc;clear;close;

+deff('y=f(x,t)','y=exp(-%pi^2*t)*sin(%pi*x)');

+u=[f(0,0) f(0.2,0) f(0.4,0) f(0.6,0) f(0.8,0) f(1,0)];

+u11=u(3)/2;u12=(u(2)+u(4))/2;u13=u12;u14=u11;

+printf(' u11=%f\t  u12=%f\t  u13=%f\t  u14=%f\n\n',u11,u12,u13,u14)

+u21=u12/2;u22=(u12+u14)/2;u23=u22;u24=u21;

+printf(' u21=%f\t  u22=%f\t  u23=%f\t  u24=%f\n\n',u21,u22,u23,u24)

+printf(' the error in the solution is: %f\n\n',abs(u22-f(0.6,0.04)))

+//crank-nicolson formula

+//by putting i=1,2,3,4 we obtain four equation

+A=[4 -1 0 0 ;-1 4 -1 0;0 -1 4 -1;0 0 -1 4];

+C=[0.9510;1.5388;1.5388;0.9510];

+X=A^-1*C;

+printf( ' u11=%f\t  u21=%f\t  u31=%f\t  u41=%f\t\n\n',X(1,1),X(2,1),X(3,1),X(4,1))

+printf(' the error in the solution is: %f\n\n',abs(X(2,1)-f(0.6,0.04)))
\ No newline at end of file
diff --git a/1673/CH9/EX9.8/9_8.sce b/1673/CH9/EX9.8/9_8.sce
new file mode 100755
index 000000000..f5ad169f6
--- /dev/null
+++ b/1673/CH9/EX9.8/9_8.sce
@@ -0,0 +1,16 @@
+//heat equation using crank-nicolson method

+//example 9.8

+//page 364

+clc;clear;close;

+U=0.01878;

+//h=1/2;l=1/8,i=1;

+u01=0;u21=1/8;

+u11=(u21+u01)/6;

+printf(' u11=%f\n\n',u11);

+printf('error is %f\n\n',abs(u11-U));

+//h=1/4,l=1/8,i=1,2,3

+A=[-3 -1 0;1 -3 1;0 1 -3];

+C=[0;0;-1/8];

+X=A^-1*C;

+printf(' u12=%f\n\n',X(2,1));

+printf('error is %f\n\n',abs(X(2,1)-U));