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

---
 854/CH7/EX7.1/Example7_1.sce | 32 ++++++++++++++++++++++++++++++++
 854/CH7/EX7.2/Example7_2.sce | 40 ++++++++++++++++++++++++++++++++++++++++
 854/CH7/EX7.3/Example7_3.sce | 23 +++++++++++++++++++++++
 854/CH7/EX7.4/Example7_4.sce | 28 ++++++++++++++++++++++++++++
 854/CH7/EX7.5/Example7_5.sce | 13 +++++++++++++
 5 files changed, 136 insertions(+)
 create mode 100755 854/CH7/EX7.1/Example7_1.sce
 create mode 100755 854/CH7/EX7.2/Example7_2.sce
 create mode 100755 854/CH7/EX7.3/Example7_3.sce
 create mode 100755 854/CH7/EX7.4/Example7_4.sce
 create mode 100755 854/CH7/EX7.5/Example7_5.sce

(limited to '854/CH7')

diff --git a/854/CH7/EX7.1/Example7_1.sce b/854/CH7/EX7.1/Example7_1.sce
new file mode 100755
index 000000000..ce89f5a75
--- /dev/null
+++ b/854/CH7/EX7.1/Example7_1.sce
@@ -0,0 +1,32 @@
+//clear//
+//Caption: Derivation of capacitance of a parallel plate capacitor
+//Example7.1
+//page 177
+clc;
+x = sym('x');
+d = sym('d');
+Vo = sym('Vo');
+e = sym('e');
+ax = sym('ax');
+A = sym('A');
+B = sym('B');
+S = sym('S');
+V = integ(A,x)+B;
+V = limit(V,A,Vo/d);
+V = limit(V,B,0);
+disp(V,'Potential in Volts V =')
+E = -diff(V,x)*ax;
+disp(E,'Electric Field in V/m E =')
+D = e*E;
+DN = D/ax;
+disp(D,'Electric Flux Density in C/square metre D =')
+Q = -DN*S;
+disp(Q,'Charge in Coulombs Q =')
+C = Q/Vo;
+disp(C,'Capacitance of parallel plate capacitor C =')
+//Result
+//Potential in Volts V =   Vo*x/d   
+//Electric Field in V/m E =   -ax*Vo/d   
+//Electric Flux Density in C/square metre D =  -ax*e*Vo/d   
+//Charge in Coulombs Q =   e*Vo*S/d   
+//Capacitance of parallel plate capacitor C =  e*S/d   
diff --git a/854/CH7/EX7.2/Example7_2.sce b/854/CH7/EX7.2/Example7_2.sce
new file mode 100755
index 000000000..facd26ba8
--- /dev/null
+++ b/854/CH7/EX7.2/Example7_2.sce
@@ -0,0 +1,40 @@
+//clear//
+//Caption: Capacitance of a Cylindrical Capacitor
+//Example7.2
+//page 179
+clc;
+A = sym('A');
+B = sym('B');
+r = sym('r');
+ar = sym('ar');
+ruo = sym('ruo');
+a = sym('a');
+b = sym('b');
+L = sym('L');
+Vo = sym('Vo');
+V = integ(A/r,r)+B;
+disp(V,'Potential V = ')
+V = limit(V,A,Vo/log(a/b));
+V = limit(V,B,-Vo*log(b)/log(a/b));
+disp(V,'Potential V by substitute the values of constant A & B = ')
+V = Vo*log(b/r)/log(b/a);
+E = -diff(V,r)*ar;
+disp(E,'E = ');
+E = limit(E,r,a);
+disp(E,'E at r =a is =')
+D = e*E;
+DN = D/ar;
+disp(DN,'DN =')
+S = float(2*%pi*a*L); //area of cylinder
+Q = DN*S
+disp(Q,'Q =')
+C = Q/Vo;
+disp(C,'Capacitance of a cylindrical Capacitor C =')
+//Result
+// Potential V =  B+log(r)*A   
+// Potential V by substitute the values of constant A & B =(log(r)-log(b))*Vo/log(a/b)   
+// E = ar*Vo/(log(b/a)*r)   
+// E at r =a is =  ar*Vo/(a*log(b/a))   
+// DN = e*Vo/(a*log(b/a))   
+// Q = 6.283185306023805*e*Vo*L/log(b/a)   
+// Capacitance of a cylindrical Capacitor C = 6.283185306023805*e*L/log(b/a)  
diff --git a/854/CH7/EX7.3/Example7_3.sce b/854/CH7/EX7.3/Example7_3.sce
new file mode 100755
index 000000000..32b9bb618
--- /dev/null
+++ b/854/CH7/EX7.3/Example7_3.sce
@@ -0,0 +1,23 @@
+//clear//
+//Caption: Program to Determine the electric field of a two infinite radial planes with an interior angle alpha
+//Example 7.3
+//page 180
+clc;
+phi = sym('phi');
+A = sym('A');
+B = sym('B');
+Vo = sym('Vo');
+alpha = sym('alpha');
+aphi = sym('aphi');
+r = sym('r');
+V = integ(A,phi)+B;
+disp(V,'V =');
+V = limit(V,B,0);
+V = limit(V,A,Vo/alpha);
+disp(V,'Potential V after applying boundary conditions =')
+E = -(1/r)*diff(V,phi)*aphi;
+disp(E,'E =')
+//Result
+// V =  B+phi*A   
+// Potential V after applying boundary conditions =  phi*Vo/alpha   
+// E =    -aphi*Vo/(alpha*r)   
diff --git a/854/CH7/EX7.4/Example7_4.sce b/854/CH7/EX7.4/Example7_4.sce
new file mode 100755
index 000000000..a74cc6f57
--- /dev/null
+++ b/854/CH7/EX7.4/Example7_4.sce
@@ -0,0 +1,28 @@
+//clear//
+//Caption: Derivation of capacitance of a spherical capacitor
+//Example7.4
+//page 181
+clc;
+a = sym('a');
+b = sym('b');
+Vo = sym('Vo');
+r = sym('r');
+e = sym('e');
+V = Vo*((1/r)-(1/b))/((1/a)-(1/b));
+disp(V,'V =')
+E = -diff(V,r)*ar;
+disp(E,'E =')
+D = e*E;
+DN = D/ar;
+disp(DN,'DN =')
+S = float(4*%pi*r^2); //area of sphere
+Q = DN*S;
+disp(Q,'Q =')
+C = Q/Vo;
+disp(C,'Capacitance of a spherical capacitor =')
+//Result
+//V =  (1/r-1/b)*Vo/(1/a-1/b)   
+//E =    ar*Vo/((1/a-1/b)*r^2)   
+//DN =   e*Vo/((1/a-1/b)*r^2)   
+//Q =    12.56637060469643*e*Vo/(1/a-1/b)   
+//Capacitance of a spherical capacitor = 12.56637060469643*e/(1/a-1/b)
diff --git a/854/CH7/EX7.5/Example7_5.sce b/854/CH7/EX7.5/Example7_5.sce
new file mode 100755
index 000000000..aa98fb99a
--- /dev/null
+++ b/854/CH7/EX7.5/Example7_5.sce
@@ -0,0 +1,13 @@
+//clear//
+//Caption: Potential in spherical coordinates as a function of teta V(teta)
+//Example7.5
+//page 182
+clc;
+teta = sym('teta');
+A = sym('A');
+B = sym('B');
+V = integ(A/float(sin(teta)),teta)+B;
+disp(V,'V = ')
+//Result
+//V =  B+(log(cos(teta)-1)/2-log(cos(teta)+1)/2)*A  
+//Equivalent to V = B+log(tan(teta/2))*A
-- 
cgit