4 files changed, 76 insertions, 0 deletions
diff --git a/3487/CH1/EX1.1/Ex1_1.sce b/3487/CH1/EX1.1/Ex1_1.sce
new file mode 100644
index 000000000..e70267c0b
--- /dev/null
+++ b/3487/CH1/EX1.1/Ex1_1.sce
@@ -0,0 +1,14 @@
+//Chapter 1,Example 1.1 Page 51
+clc
+clear
+I = 600 // micor amps
+x = 0.5 // distance in cm
+V = 10 // kV
+I2 = 60 // micro amps
+x2 = 0.1 // distance in cm 
+//Calculation 600 = I0*exp(0.5*alpha) and 60 = I0*exp(0.1*alpha)
+alpha = log(600/60)/(0.5-0.1)
+printf("Townsends first ionising coefficient = %f ionizing collisions/cm", alpha)
+
+//Answers may vary due to round of error  
+
diff --git a/3487/CH1/EX1.2/Ex1_2.sce b/3487/CH1/EX1.2/Ex1_2.sce
new file mode 100644
index 000000000..0149dc124
--- /dev/null
+++ b/3487/CH1/EX1.2/Ex1_2.sce
@@ -0,0 +1,27 @@
+//Chapter 1,Example 1.2 Page 52
+clc
+clear
+// Refering the table in example 1.2
+// slope between any two points (log(I/I0)/x)
+// taking the gap between 2 and 2.5 mm
+I1= 1.5*10^-12
+I2= 5.6*10^-12
+I0 = 6*10^-14
+gi1 = log(I1/I0) // gradual increase when gap is 2
+gi2 = log(I2/I0) // gradual increase when gap is 2.5 //claculation in text is wrong
+slope = (gi1-gi2)/0.05
+printf(" Slope = %f \n", -slope) 
+//evaluvating ghama
+e1 = exp(-slope*0.5)
+e2 = exp(-slope*0.5) // -1 is ignored due to the large magnitude
+ghama = (7*10^7-6*e1)/(e2*7*10^7)
+printf(" Ghama for set 1= %e /cm \n ", ghama)
+//Gap between the slope for set 2
+alpha = log(12/8)/0.05
+printf(" Alpha = %e collosions/cm \n", alpha)
+e1 = exp(alpha*0.5)
+e2 = exp(alpha*0.5) // -1 is ignored due to the large magnitude
+ghama = (2*10^5-e1)/(e2*2*10^5)
+printf(" Ghama for set 2= %e colissions/cm \n", ghama)
+
+//Answers may vary due to round of error  
diff --git a/3487/CH1/EX1.3/Ex1_3.sce b/3487/CH1/EX1.3/Ex1_3.sce
new file mode 100644
index 000000000..cb7cec2db
--- /dev/null
+++ b/3487/CH1/EX1.3/Ex1_3.sce
@@ -0,0 +1,21 @@
+//Chapter 1,Example 1.3 Page 53
+clc
+clear
+//employing equation Vb = K*d^n
+//88 = K*4^n --- eq(1)  165 = K*8^n ---eq(2) 
+//dividing eq(2)/q(1)
+Vb1 = 88
+Vb2 = 165
+n1 = 0.6286/0.693
+K1 = Vb1/4^n1
+//135 = K*6^n --- eq(1)  212 = K*10^n ---eq(2) 
+//dividing eq(2)/q(1) 
+Vb1 = 135
+Vb2 = 212
+n2 = 0.4513/0.5128
+K2 = Vb1/6^n2
+n = (n1+n2)/2
+K = (K1+K2)/2
+printf (" n = %f (approx.) K = %f (approx.)",n,K)
+
+//Answer may vary due to round of error 
diff --git a/3487/CH1/EX1.4/Ex1_4.sce b/3487/CH1/EX1.4/Ex1_4.sce
new file mode 100644
index 000000000..6d8b71e40
--- /dev/null
+++ b/3487/CH1/EX1.4/Ex1_4.sce
@@ -0,0 +1,14 @@
+//Chapter 1,Example 1.4 Page 53
+// Determine (pd)min Vbmin
+clc
+clear
+A = 12
+B = 365
+e = 2.718
+ghama = 0.02
+K = 51
+pd = (e/A)*log(1+(1/ghama))
+Vbmin = (B/A)*e*log(K)
+printf (" (pd)min = %f Vbmin = %f Volts",pd,Vbmin)
+
+//Answers may vary due to round of error