5 files changed, 73 insertions, 0 deletions
diff --git a/1931/CH13/EX13.1/1.sce b/1931/CH13/EX13.1/1.sce
new file mode 100755
index 000000000..ae3f17729
--- /dev/null
+++ b/1931/CH13/EX13.1/1.sce
@@ -0,0 +1,15 @@
+clc
+clear
+//INPUT DATA
+a=3.61*10^-10//lattice constant of copper which is Fcc crystal in m
+x=1*10^-18//average displacement of the electrons relative to the nucleus in m
+z=29//atomic number of copper
+n=4//number of atoms per unit cell in FCC crystal
+e=1.6*10^-19//charge of electron in coulombs
+
+//CALCULATION
+ne=((n*z)/(a*a*a))//number of electrons in electrons/m^3 
+P=(ne*e*x)/(10^-7)//The electron polarisation in C/m^2 *10^-7
+
+//OUTPUT
+printf('The electron polarisation is %3.3f*10^-7 in C/m^2',P)
diff --git a/1931/CH13/EX13.2/2.sce b/1931/CH13/EX13.2/2.sce
new file mode 100755
index 000000000..f8f8d2891
--- /dev/null
+++ b/1931/CH13/EX13.2/2.sce
@@ -0,0 +1,18 @@
+clc
+clear
+//INPUT DATA
+rp=11.7//relative permittivity of silicon
+N=4.82*10^28//number of atoms per unit volume in atoms/m^3
+ro=8.85*10^-12//permittivity of free space
+E=10^4//E in Vm^-1
+e=1.6*10^-19//charge of electron in coulombs
+Z=14//atomic number of silicon 
+
+//CALCULATION
+z=((ro*(rp-1))/N)//electronic polarisability in Fm^2
+y=(z*E)/10^-35//The dipole moment of each atom in a field of 10^4 Vm^-1 in Cm^-3 *10^-35
+x=((y*10^-35)/(Z*e))/10^-18//The effective distance at this field strength between the centre of the electron cloud in each atom and the nucleus in m*10^-18
+
+//OUTPUT
+printf('The dipole moment of each atom in a field is %3.4f*10^-35 Cm^-3 \n The effective distance at this field strength between the centre of the electron cloud in each atom and the nucleus is %3.2f*10^-18 m',y,x)
+
diff --git a/1931/CH13/EX13.3/3.sce b/1931/CH13/EX13.3/3.sce
new file mode 100755
index 000000000..994e2bf54
--- /dev/null
+++ b/1931/CH13/EX13.3/3.sce
@@ -0,0 +1,13 @@
+clc
+clear
+//INPUT DATA
+d=9.8*10^26//density of hydrogen gas in atoms/m^3
+r=0.50*10^-10//radius of the hydrogen atom in m
+ro=8.85*10^-12//permittivity of free space
+
+//CALCULATION
+z=(4*3.14*ro*r*r*r)/10^-41//electronic polarisability in Fm^2
+rp=(((d*z*10^-41)/ro)+1)//The relative permittivity in hydrogen gas
+
+//OUTPUT
+printf('The electronic polarisability is %3.3f*10^-41 Fm^2 \n The relative permittivity in hydrogen gas is %3.4f',z,rp)
diff --git a/1931/CH13/EX13.4/4.sce b/1931/CH13/EX13.4/4.sce
new file mode 100755
index 000000000..0131ff132
--- /dev/null
+++ b/1931/CH13/EX13.4/4.sce
@@ -0,0 +1,15 @@
+clc
+clear
+//INPUT DATA
+z=1.75*10^-40//electronic polarisability in Fm^2
+d=1.8*10^3//density of argon atom in Kg/m^3
+Z=39.95//atomic weight of argon
+NA=6.025*10^26//Avagadro number in mole^-1
+ro=8.85*10^-12//permittivity of free space
+
+//CALCULATION
+N=((NA*d)/Z)//The number of atoms/unit volume in atoms/m^3 
+rp=(((N*z)/ro)+1)//The static dielectric constant of solid argon
+
+//OUTPUT
+printf('The static dielectric constant of solid argon is %3.5f',rp)
diff --git a/1931/CH13/EX13.5/5.sce b/1931/CH13/EX13.5/5.sce
new file mode 100755
index 000000000..bcae4482a
--- /dev/null
+++ b/1931/CH13/EX13.5/5.sce
@@ -0,0 +1,12 @@
+clc
+clear
+//INPUT DATA
+er=4.94//static dielecric constant of a material
+n=2.69//n is the index of friction
+
+//CALCULATION
+x=((er-1)*(n+2))/((er+2)*(n-1))-1//Ratio between ionic and electronic polarisability of this material
+y=1/x//Ratio between electronic and ionic polarisability of this material
+
+//OUTPUT
+printf('Ratio between electronic and ionic polarisability of this material is %3.4f ',y)