initial commit / add all books

5 files changed, 141 insertions, 0 deletions

diff --git a/854/CH3/EX3.1/Example3_1.sce b/854/CH3/EX3.1/Example3_1.sce
new file mode 100755
index 000000000..757c1f4e1
--- /dev/null
+++ b/854/CH3/EX3.1/Example3_1.sce
@@ -0,0 +1,14 @@
+//clear//

+//Caption: Program to find Electric Flux density 'D' of a uniform line charge

+//Example3.1

+//page 54

+clc;

+e0 = 8.854e-12; //free space permittivity in F/m

+rL = 8e-09; //line charge density c/m

+r = 3; // distance in metre

+E = Electric_Field_Line_Charge(rL,e0,r); //electric field intensity of line charge

+D = e0*E;

+disp(D,'Electric Flux Density in Coulombs per square metre D =')

+//Result

+// Electric Flux Density in Coulombs per square metre D =   

+//     4.244D-10  

diff --git a/854/CH3/EX3.2/Example3_2.sce b/854/CH3/EX3.2/Example3_2.sce
new file mode 100755
index 000000000..431d53ecf
--- /dev/null
+++ b/854/CH3/EX3.2/Example3_2.sce
@@ -0,0 +1,28 @@
+//clear//

+//Caption: Program to calculate surface charge density,Flux density, Field Intensity of coaxial cable

+//Example3.2

+//page 64

+clc;

+Q_innercyl = 30e-09; //total charge on the inner conductor in coulombs

+a = 1e-03; // inner radius of coaxial cable in metre

+b = 4e-03; // outer radius of coaxial cable in metre

+L = 50e-02; //length of coaxial cable

+rs_innercyl = Q_innercyl/(2*%pi*a*L);

+rs_outercyl = Q_innercyl/(2*%pi*b*L);

+e0 = 8.854e-12; //free space relative permittivity F/m

+r = sym('r');

+Dr = a*rs_innercyl/r;

+Er = Dr/e0;

+disp(rs_innercyl,'Surface charge density of inner cylinder of coaxial cable in C/square.metre, rs_innercyl=')

+disp(rs_outercyl,'Surface charge density of outer cylinder of coaxial cable in C/square.metre, rs_outercyl=')

+disp(Dr,'Electric Flux Density in C/square.metre Dr=')

+disp(Er,'Electric Field Intensity in V/m Er=')

+//Result

+//Surface charge density of inner cylinder of coaxial cable in C/square.metre, rs_innercyl=   

+//    0.0000095  

+//Surface charge density of outer cylinder of coaxial cable in C/square.metre, rs_outercyl=   

+//    0.0000024  

+//Electric Flux Density in C/square.metre Dr=   

+// 9.5488183337312011E-9/r   

+//Electric Field Intensity in V/m Er=   

+// 1078.47507722286/r  

diff --git a/854/CH3/EX3.3/Example3_3.sce b/854/CH3/EX3.3/Example3_3.sce
new file mode 100755
index 000000000..1cda28ca9
--- /dev/null
+++ b/854/CH3/EX3.3/Example3_3.sce
@@ -0,0 +1,29 @@
+//clear//

+//Caption: Program to calculate the total charge enclosed in a volume at the origin

+//Example3.3

+//page 67

+clc;

+V = 1e-09; //volume in cubic metre

+x = sym('x');

+y = sym('y');

+z = sym('z');

+//Components of Electric Flux Density in cartesian coordinate system

+Dx = exp(-x)*sin(y);

+Dy = -exp(-x)*cos(y);

+Dz = 2*z;

+//Divergence of electric flux density 'D'

+dDx = diff(Dx,x);

+dDy = diff(Dy,y);

+dDz = diff(Dz,z);

+//Total charge enclosed in a given volume

+del_Q = (dDx+dDy+dDz)*V;

+disp(del_Q,'Total charge enclosed in an incremental volume in coulombs, del_Q =')

+//Total Charge enclosed in a given volume at origin (0,0,0)

+del_Q = limit(del_Q,x,0);

+del_Q = limit(del_Q,y,0);

+del_Q = limit(del_Q,z,0);

+disp(del_Q*1e09,'Total charge enclosed in an incremental volume in  nano coulombs at origin, del_Q =')

+//Result

+//Total charge enclosed in an incremental volume in coulombs, del_Q =    2.0000000000000001E-9   

+//Total charge enclosed in an incremental volume in  nano coulombs at origin, del_Q =   

+//  2.0 

diff --git a/854/CH3/EX3.4/Example3_4.sce b/854/CH3/EX3.4/Example3_4.sce
new file mode 100755
index 000000000..dd647c7a9
--- /dev/null
+++ b/854/CH3/EX3.4/Example3_4.sce
@@ -0,0 +1,27 @@
+//clear//

+//Caption: Program to Find the Divergence of 'D' at the origin

+//Example3.4

+//page 70

+clc;

+x = sym('x');

+y = sym('y');

+z = sym('z');

+//Components of Electric Flux Density in cartesian coordinate system

+Dx = exp(-x)*sin(y);

+Dy = -exp(-x)*cos(y);

+Dz = 2*z;

+//Divergence of electric flux density 'D'

+dDx = diff(Dx,x);

+dDy = diff(Dy,y);

+dDz = diff(Dz,z);

+divD = dDx+dDy+dDz

+disp(divD,'Divergence of Electric Flux Density D in C/cubic.metre, divD =')

+divD = limit(divD,x,0);

+divD = limit(divD,y,0);

+divD = limit(divD,z,0);

+disp(divD,'Divergence of Electric Flux Density D in C/cubic.metre at origin, divD =')

+//Result

+//Divergence of Electric Flux Density D in C/cubic.metre, divD =   

+// 2   

+//Divergence of Electric Flux Density D in C/cubic.metre at origin, divD =   

+// 2   

diff --git a/854/CH3/EX3.5/Example3_5.sce b/854/CH3/EX3.5/Example3_5.sce
new file mode 100755
index 000000000..f8d2cf549
--- /dev/null
+++ b/854/CH3/EX3.5/Example3_5.sce
@@ -0,0 +1,43 @@
+//clear//

+//Caption: Program to verify the Divergence theorem for the field 'D' 

+//Example3.5

+//page 74

+clc;

+x = sym('x');

+y = sym('y');

+z = sym('z');

+//Components of Electric Flux Density in cartesian coordinate system

+Dx = 2*x*y;

+Dy = x^2;

+Dz = 0;

+//Divergence of electric flux density 'D'

+dDx = diff(Dx,x);

+dDy = diff(Dy,y);

+dDz =0;

+divD = dDx+dDy+dDz

+disp(divD,'Divergence of Electric Flux Density D in C/cubic.metre, divD =')

+//Evaluate volume integral on divergence of 'D'

+Vol_int_divD = integ(divD,x);

+Vol_int_divD = limit(Vol_int_divD,x,1)-limit(Vol_int_divD,x,0);

+Vol_int_divD = integ(Vol_int_divD,y);

+Vol_int_divD = limit(Vol_int_divD,y,2)-limit(Vol_int_divD,y,0);

+Vol_int_divD = integ(Vol_int_divD,z);

+Vol_int_divD = limit(Vol_int_divD,z,3)-limit(Vol_int_divD,z,0);

+disp(Vol_int_divD,'Volume Integral of divergence of D, in coulombs vol_int(divD)=')

+//Evaluate surface integral on field D

+Dx = limit(Dx,x,1);

+sur_D = integ(Dx,y);

+sur_D = limit(sur_D,y,2) - limit(sur_D,y,0);

+sur_D = integ(sur_D,z);

+sur_D = limit(sur_D,z,3) - limit(sur_D,z,0);

+disp(sur_D,'Surface Integral of field D, in coulombs sur_int(D.ds)=')

+if(sur_D==Vol_int_divD)

+  disp('Divergence Theorem verified')

+end

+//Result

+// Divergence of Electric Flux Density D in C/cubic.metre, divD =   

+//  2*y   

+//Volume Integral of divergence of D, in coulombs vol_int(divD)=   

+// 12   

+// Surface Integral of field D, in coulombs sur_int(D.ds)=   

+// 12