7 files changed, 77 insertions, 0 deletions
diff --git a/1535/CH13/EX13.1/Ch13Ex1.sci b/1535/CH13/EX13.1/Ch13Ex1.sci
new file mode 100755
index 000000000..4a60c9c23
--- /dev/null
+++ b/1535/CH13/EX13.1/Ch13Ex1.sci
@@ -0,0 +1,9 @@
+// Scilab Code Ex13.1: Electronic Polarizability of atom : Page-287 (2010)
+epsilon_0 = 8.854e-012;    // Absolute electrical permittivity of free space, farad per metre
+R = 0.52e-010;    // Radius of hydrogen atom, angstrom
+n = 9.7e+026;    // Number density of hydrogen, per metre cube
+alpha_e = 4*%pi*epsilon_0*R^3;    // Electronic polarizability of hydrogen atom, farad-metre square
+printf("\nThe electronic polarizability of hydrogen atom = %4.2e farad-metre square", alpha_e);
+
+// Result
+// The electronic polarizability of hydrogen atom = 1.56e-041 farad-metre square
+\ No newline at end of file
diff --git a/1535/CH13/EX13.2/Ch13Ex2.sci b/1535/CH13/EX13.2/Ch13Ex2.sci
new file mode 100755
index 000000000..16929120e
--- /dev/null
+++ b/1535/CH13/EX13.2/Ch13Ex2.sci
@@ -0,0 +1,13 @@
+// Scilab Code Ex13.2: Parallel plate capacitor: Page-287 (2010)
+epsilon_0 = 8.854e-012;    // Absolute electrical permittivity of free space, farad per metre
+A = 100e-004;    // Area of a plate of parallel plate capacitor, metre square
+d = 1e-002;    // Distance between the plates of the capacitor, m
+V = 100;    // Potential applied to the plates of the capacitor, volt
+C = epsilon_0*A/d;    // Capacitance of parallel plate capacitor, farad
+Q = C/V;    // Charge on the plates of the capacitor, coulomb
+printf("\nThe capacitance of parallel plate capacitor = %5.3e F", C);
+printf("\nThe charge on the plates of the capacitor = %5.3e C", Q);
+
+// Result
+// The capacitance of parallel plate capacitor = 8.854e-012 F
+// The charge on the plates of the capacitor = 8.854e-014 C
+\ No newline at end of file
diff --git a/1535/CH13/EX13.3/Ch13Ex3.sci b/1535/CH13/EX13.3/Ch13Ex3.sci
new file mode 100755
index 000000000..921e43b56
--- /dev/null
+++ b/1535/CH13/EX13.3/Ch13Ex3.sci
@@ -0,0 +1,11 @@
+// Scilab Code Ex13.3: Dielectric displacement of medium: Page-288 (2010)
+epsilon_0 = 8.854e-012;    // Absolute electrical permittivity of free space, farad per metre
+epsilon_r = 5.0;    // Dielectric constant of the material between the plates of capacitor
+V = 15;    // Potential difference applied between the plates of the capacitor, volt
+d = 1.5e-003;    // Separation between the plates of the capacitor, m
+// Electric displacement, D = epsilon_0*epsilon_r*E, as E = V/d, so 
+D = epsilon_0*epsilon_r*V/d;    // Dielectric displacement, coulomb per metre square
+printf("\nThe dielectric displacement = %5.3e coulomb per metre square", D);
+
+// Result
+// The dielectric displacement = 4.427e-007 coulomb per metre square 
+\ No newline at end of file
diff --git a/1535/CH13/EX13.4/Ch13Ex4.sci b/1535/CH13/EX13.4/Ch13Ex4.sci
new file mode 100755
index 000000000..ce398350c
--- /dev/null
+++ b/1535/CH13/EX13.4/Ch13Ex4.sci
@@ -0,0 +1,9 @@
+// Scilab Code Ex13.4: Relative dielectric constant : Page-288 (2010)
+epsilon_0 = 8.854e-012;    // Absolute electrical permittivity of free space, farad per metre
+N = 3.0e+028;    // Number density of solid elemental dielectric, atoms per metre cube
+alpha_e = 1e-040;    // Electronic polarizability, farad metre square
+epsilon_r = 1 + N*alpha_e/epsilon_0;    // Relative dielectric constant of the material
+printf("\nThe Relative dielectric constant of the material = %5.3f", epsilon_r);
+
+// Result
+// The Relative dielectric constant of the material = 1.339 
+\ No newline at end of file
diff --git a/1535/CH13/EX13.5/Ch13Ex5.sci b/1535/CH13/EX13.5/Ch13Ex5.sci
new file mode 100755
index 000000000..fea5820fd
--- /dev/null
+++ b/1535/CH13/EX13.5/Ch13Ex5.sci
@@ -0,0 +1,16 @@
+// Scilab Code Ex13.5: Atomic polarizability of sulphur : Page-288 (2010)
+N_A = 6.023e+023;    // Avogadro's number, per mole
+epsilon_0 = 8.854e-012;    // Absolute electrical permittivity of free space, farad per metre
+epsilon_r = 3.75;    // Relative dielectric constant
+d = 2050;    // Density of sulphur, kg per metre cube
+y = 1/3;    // Internal field constant
+M = 32;    // Atomic weight of sulphur, g/mol
+N = N_A*1e+03*d/M;    // Number density of atoms of sulphur, per metre cube
+// Lorentz relation for local fields give
+// E_local = E + P/(3*epsilon_0) which gives
+// (epsilon_r - 1)/(epsilon_r + 2) = N*alpha_e/(3*epsilon_0), solving for alpha_e
+alpha_e = (epsilon_r - 1)/(epsilon_r + 2)*3*epsilon_0/N;    // Electronic polarizability of sulphur, farad metre square
+printf("\nThe electronic polarizability of sulphur = %5.3e farad metre square", alpha_e);
+
+// Result
+// The electronic polarizability of sulphur = 3.292e-040 farad metre square
+\ No newline at end of file
diff --git a/1535/CH13/EX13.6/Ch13Ex6.sci b/1535/CH13/EX13.6/Ch13Ex6.sci
new file mode 100755
index 000000000..1f99790eb
--- /dev/null
+++ b/1535/CH13/EX13.6/Ch13Ex6.sci
@@ -0,0 +1,10 @@
+// Scilab Code Ex13.6: Electronic polarizability from refractive index : Page-289 (2010)
+N = 3e+028;    // Number density of atoms of dielectric material, per metre cube
+epsilon_0 = 8.854e-012;    // Absolute electrical permittivity of free space, farad per metre
+n = 1.6;    // Refractive index of dielectric material
+// As (n^2 - 1)/(n^2 + 2) = N*alpha_e/(3*epsilon_0), solving for alpha_e
+alpha_e = (n^2 - 1)/(n^2 + 2)*3*epsilon_0/N;    // Electronic polarizability of dielectric material, farad metre square
+printf("\nThe electronic polarizability of dielectric material = %4.2e farad metre square", alpha_e);
+
+// Result
+// The electronic polarizability of dielectric material = 3.03e-040 farad metre square 
+\ No newline at end of file
diff --git a/1535/CH13/EX13.7/Ch13Ex7.sci b/1535/CH13/EX13.7/Ch13Ex7.sci
new file mode 100755
index 000000000..2a096a959
--- /dev/null
+++ b/1535/CH13/EX13.7/Ch13Ex7.sci
@@ -0,0 +1,9 @@
+// Scilab Code Ex13.7: Ratio of electronic polarizability to ionic polarizability: Page-289 (2010)
+epsilon_r = 4.9;    // Absolute relative dielectric constant of material, farad per metre
+n = 1.6;    // Refractive index of dielectric material
+// As (n^2 - 1)/(n^2 + 2)*(alpha_e + alpha_i)/alpha_e = N*(alpha_e + alpha_i)/(3*epsilon_0) = (epsilon_r - 1)/(epsilon_r + 2), solving for alpha_i/alpha_e
+alpha_ratio = ((epsilon_r - 1)/(epsilon_r + 2)*(n^2 + 2)/(n^2 - 1) - 1)^(-1);    // Ratio of electronic polarizability to ionic polarizability
+printf("\nThe ratio of electronic polarizability to ionic polarizability = %4.2f", alpha_ratio);
+
+// Result
+// The ratio of electronic polarizability to ionic polarizability = 1.53 
+\ No newline at end of file