initial commit / add all books

11 files changed, 198 insertions, 0 deletions

diff --git a/3369/CH4/EX4.1/Ex4_1.sce b/3369/CH4/EX4.1/Ex4_1.sce
new file mode 100755
index 000000000..4a9e91ba7
--- /dev/null
+++ b/3369/CH4/EX4.1/Ex4_1.sce
@@ -0,0 +1,17 @@
+//Chapter 4, Exmaple 1, page 139
+//Claculate alpha and No. of electrons emmited
+clc
+clear
+//Claculate (a)alpha
+d2 = 0.01
+d1 = 0.005
+I2 = 2.7*10**-7
+I1 = 2.7*10**-8
+alpha = 1/(d2-d1)*log(I2/I1)
+//(b)number of electrons emmited from cathode per second
+I0 = I1*%e**(-alpha*d1)
+n0 = I0/(1.6*10**-19)
+printf("\n Part (a)\n alpha = %f m^-1",alpha)
+printf("\n Part (b)\n I0 = %e ",I0)
+printf("\n No of electrons emitted = %e electrons/s",n0)
+//Answer may vary due to round off error
diff --git a/3369/CH4/EX4.10/Ex4_10.sce b/3369/CH4/EX4.10/Ex4_10.sce
new file mode 100755
index 000000000..ab65e2563
--- /dev/null
+++ b/3369/CH4/EX4.10/Ex4_10.sce
@@ -0,0 +1,28 @@
+//Chapter 4, Exmaple 10, page 144
+//Claculate (a)Raether's criterion (b)Meek and Lobe's criterion
+clc
+clear
+//(a)Raether's criterion
+// as assumed by Raether and based equation 3.3, 3.50, 4.22 and 4.23
+d = 0.001 // m
+alpha = 10792.2 // m^-1
+p = 101.3 //kPa^-1
+ap = 106.54 // alpha/p Unit: m^-1*kPa^-1
+T = 11253.7 // m^-1*kPa^-1
+B = 273840 // V/m*kPa
+Ep = 58764.81 // E/p Unit:V/m*kPa
+
+ad = 17.7 + log(d)
+E = Ep*p
+Vs = E*d*10^-3 //Voltage breakdown
+printf("\n E = %e V/m",E)
+printf("\n Voltage breakdown = %f kV",Vs)
+
+//(b)Meek and Loeb's criterion
+//Using equation 4.11 and based on 4.24 & 4,25 
+//+ we get Er = 468*10^4 V/m
+Er = 468*10^4 // V/m
+Vs2 = Er*0.001*10^-3
+printf("\n Voltage breakdown = %f kV",Vs2)
+
+// Answers may vary due to round of error
diff --git a/3369/CH4/EX4.11/Ex4_11.sce b/3369/CH4/EX4.11/Ex4_11.sce
new file mode 100755
index 000000000..8efd2abd0
--- /dev/null
+++ b/3369/CH4/EX4.11/Ex4_11.sce
@@ -0,0 +1,13 @@
+//Chapter 4, Exmaple 11, page 146
+//Claculate the first Townsend's ionization coefficient
+clc
+clear
+t = 0.2*10**-6 // transit time of electrons in seconds
+d = 0.05 // m
+ve = d/t
+TC = 35*10**-9 // Time constant
+a = 1/(ve*TC)
+printf("\n Electron drift velocity = %e m/s",ve)
+printf("\n alpha = %e m^-1",a)
+
+// Answers may vary due to round of error
diff --git a/3369/CH4/EX4.12/Ex4_12.sce b/3369/CH4/EX4.12/Ex4_12.sce
new file mode 100755
index 000000000..6e591561e
--- /dev/null
+++ b/3369/CH4/EX4.12/Ex4_12.sce
@@ -0,0 +1,19 @@
+//Chapter 4, Exmaple 12, page 146
+//Travel time and maximum frequency
+clc
+clear
+//(a)Determine the travel time
+Ea = 200*sqrt(2)*10**3/0.1
+x = 1.4*10**-4*2828.4*10**3/(2*%pi*50)
+d = 0.1
+printf("\n Ea = %e V/m",Ea)
+printf("\n x = %f*sin(3.14*t)",x)
+//obtaining t from x
+t = asin(d/x)/3.14
+printf("\n t = %f ms",t) // answer mentioned in the text is wrong
+//(b)Determine the maximum frequency
+k = 1.4*10**-4
+fmax = k*Ea/(2*%pi*d)
+printf("\n fmax = %f Hz",fmax)
+
+//Answer may vary due to round off error
diff --git a/3369/CH4/EX4.2/Ex4_2.sce b/3369/CH4/EX4.2/Ex4_2.sce
new file mode 100755
index 000000000..4c6c80ee8
--- /dev/null
+++ b/3369/CH4/EX4.2/Ex4_2.sce
@@ -0,0 +1,18 @@
+//Chapter 4, Exmaple 2, page 140
+//Claculate electrode space
+clc
+clear
+//based on the values of example 1
+d2 = 0.01
+d1 = 0.005
+I2 = 2.7*10**-7
+I1 = 2.7*10**-8
+a = 1/(d2-d1)*log(I2/I1) // alpha
+//10^9 = %e^a(a*d) 
+//multiplying log on bith sides log(10^9) = a*d
+ad = log(10^9)
+printf("\n a*d = %f ",ad)
+d = ad/a
+printf("\n electrode space = %f m",d) 
+
+//Answers may vary due to round off error
diff --git a/3369/CH4/EX4.3/Ex4_3.sce b/3369/CH4/EX4.3/Ex4_3.sce
new file mode 100755
index 000000000..5858a2127
--- /dev/null
+++ b/3369/CH4/EX4.3/Ex4_3.sce
@@ -0,0 +1,13 @@
+//Chapter 4, Exmaple 3, page 140
+//Claculate size of developed avalanche
+clc
+clear
+a = 4*10**4
+b = 15*10**5
+//Rewriting equation 4.2
+x0=0;x1=0.0005;
+X=integrate('a-b*sqrt(x)','x',x0,x1);
+As = exp(X) // Avelanche size
+printf("\n Avalanche size = %f m",As)
+ 
+//Answers may vary due to round of error
diff --git a/3369/CH4/EX4.4/Ex4_4.sce b/3369/CH4/EX4.4/Ex4_4.sce
new file mode 100755
index 000000000..1db17855f
--- /dev/null
+++ b/3369/CH4/EX4.4/Ex4_4.sce
@@ -0,0 +1,13 @@
+//Chapter 4, Exmaple 4, page 141
+//Claculate distance to produce avalanche
+clc
+clear
+//Rewrite equation 4.2
+//using the values of a and b from previous example
+//convert integartion to quaderatic equation form
+x=poly(0,"x");
+p=59.97-4*10**4*x+7.5*10**5*x^2 // making the polinomial equation
+r= roots(p) //obtaining the roots
+printf("\n %f m or %f m away from the cathode",r(1),r(2))
+
+//Answer may vary due to round of error.
diff --git a/3369/CH4/EX4.5/Ex4_5.sce b/3369/CH4/EX4.5/Ex4_5.sce
new file mode 100755
index 000000000..e78b58fc0
--- /dev/null
+++ b/3369/CH4/EX4.5/Ex4_5.sce
@@ -0,0 +1,11 @@
+//Chapter 4, Exmaple 5, page 141
+//Claculate minimum distance to produce avalanche of size 10^19
+clc
+clear
+//Rewriting equation 4.2 and converting it into quadratic equation
+x=poly(0,"x");
+p=43.75-4*10**4*x+7.5*10**5*x^2 // making the polinomial equation
+r= roots(p) //obtaining the roots
+printf("\n Minimum distance = %f m",r(2)) // other root is disregarded
+
+//Answer may vary due to round of error.
diff --git a/3369/CH4/EX4.7/Ex4_7.sce b/3369/CH4/EX4.7/Ex4_7.sce
new file mode 100755
index 000000000..b388e41e3
--- /dev/null
+++ b/3369/CH4/EX4.7/Ex4_7.sce
@@ -0,0 +1,17 @@
+//Chapter 4, Exmaple 7, page 142
+//Claculate secondary coefficient
+clc
+clear
+//Using equation 3.15
+E = 9*10**3/0.002
+T = 11253.7 // m^-7*kPa^-1
+B = 273840 // V/mkPa
+p = 101.3 // kPa or 1 atm
+d = 0.002 // m
+alpha = p*T*exp(-B*p/E)
+Y = 1/(exp(alpha*d)-1)
+printf("\n E = %e V/m",E)
+printf("\n Alpha = %f m^-1",alpha)
+printf("\n Total secondary coefficient of ionization = %f ",Y)
+
+//Answer may vary due to round off error
diff --git a/3369/CH4/EX4.8/Ex4_8.sce b/3369/CH4/EX4.8/Ex4_8.sce
new file mode 100755
index 000000000..e4c3f69cd
--- /dev/null
+++ b/3369/CH4/EX4.8/Ex4_8.sce
@@ -0,0 +1,31 @@
+//Chapter 4, Exmaple 8, page 143
+//Claculate first and secondary ionization coefficient
+clc
+clear
+//(a)first ionization coefficient
+//Using equation 4.7a
+d1 = 0.005
+a1d1 = log(1.22)
+a1 = a1d1/d1
+
+d2 = 0.01504
+a2d2 = log(1.82)
+a2 = a2d2/d2
+
+d3 = 0.019 // wrong value used in the text
+a3d3 = log(2.22)
+a3 = a3d3/d3
+
+printf("\n Alpha 1 = %f m^-1",a1)
+printf("\n Alpha 2 = %f m^-1",a2)
+printf("\n Alpha 3 = %f m^-1",a3)
+printf("\n From the above results we can understand that ionization mechanism must be acting at d3 ")
+
+//secondary ionization coefficient
+I = 2.22
+e = exp(a1*d3)
+Y = (I-e)/(I*(e-1))
+printf("\n secondary ionization coefficient = %f ",Y)
+
+//Answer may vary due to round off error.
+
diff --git a/3369/CH4/EX4.9/Ex4_9.sce b/3369/CH4/EX4.9/Ex4_9.sce
new file mode 100755
index 000000000..121fb59c9
--- /dev/null
+++ b/3369/CH4/EX4.9/Ex4_9.sce
@@ -0,0 +1,18 @@
+//Chapter 4, Exmaple 9, page 144
+//Claculate distance and voltage
+clc
+clear
+a = 39.8 // alpha
+Y = 0.0354 // corfficient
+p = 0.133 // kPa
+Ep = 12000 // E/P , unit : V/m*kPa
+
+d = (1/a)*(log(1/Y + 1)) // distance
+E = Ep*p
+V = E*d
+
+printf("\n Distance = %f m",d)
+printf("\n E = %f V/m",E)
+printf("\n Volatge  = %f V",V)
+
+//Answers may vary due to round off error