summaryrefslogtreecommitdiff
path: root/854/CH7
diff options
context:
space:
mode:
authorpriyanka2015-06-24 15:03:17 +0530
committerpriyanka2015-06-24 15:03:17 +0530
commitb1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (patch)
treeab291cffc65280e58ac82470ba63fbcca7805165 /854/CH7
downloadScilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.gz
Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.bz2
Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.zip
initial commit / add all books
Diffstat (limited to '854/CH7')
-rwxr-xr-x854/CH7/EX7.1/Example7_1.sce32
-rwxr-xr-x854/CH7/EX7.2/Example7_2.sce40
-rwxr-xr-x854/CH7/EX7.3/Example7_3.sce23
-rwxr-xr-x854/CH7/EX7.4/Example7_4.sce28
-rwxr-xr-x854/CH7/EX7.5/Example7_5.sce13
5 files changed, 136 insertions, 0 deletions
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