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

---
 1673/CH9/EX9.1/9_1.sce | 11 +++++++++++
 1673/CH9/EX9.2/9_2.sce | 29 +++++++++++++++++++++++++++++
 1673/CH9/EX9.4/9_4.sce | 13 +++++++++++++
 1673/CH9/EX9.6/9_6.sce | 26 ++++++++++++++++++++++++++
 1673/CH9/EX9.7/9_7.sce | 19 +++++++++++++++++++
 1673/CH9/EX9.8/9_8.sce | 16 ++++++++++++++++
 6 files changed, 114 insertions(+)
 create mode 100755 1673/CH9/EX9.1/9_1.sce
 create mode 100755 1673/CH9/EX9.2/9_2.sce
 create mode 100755 1673/CH9/EX9.4/9_4.sce
 create mode 100755 1673/CH9/EX9.6/9_6.sce
 create mode 100755 1673/CH9/EX9.7/9_7.sce
 create mode 100755 1673/CH9/EX9.8/9_8.sce

(limited to '1673/CH9')

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));
-- 
cgit