6 files changed, 96 insertions, 0 deletions
diff --git a/3487/CH4/EX4.1/Ex4_1.sce b/3487/CH4/EX4.1/Ex4_1.sce
new file mode 100644
index 000000000..915af16c8
--- /dev/null
+++ b/3487/CH4/EX4.1/Ex4_1.sce
@@ -0,0 +1,17 @@
+//Chapter 4,Example 4.1 Page 144
+clc
+clear
+//Determine the voltge when S=2 cm 
+S = 0.2 // cm
+Vb = 24.22*S+6.08*sqrt(S)
+printf (" Vb when S = 2 cm is %f kV \n ",Vb)
+//Determine the voltge when S=1.5 cm 
+S = 1.5 // cm
+Vb = 24.22*S+6.08*sqrt(S)
+printf (" Vb when S = 1.5 is %f kV \n ",Vb)
+b = 75
+t = 35
+D = (3.92*b)/(273+t)
+printf (" Air density correction factor= %f \n ",D)
+
+//Answer may vary due to round off error
diff --git a/3487/CH4/EX4.2/Ex4_2.sce b/3487/CH4/EX4.2/Ex4_2.sce
new file mode 100644
index 000000000..01741bebc
--- /dev/null
+++ b/3487/CH4/EX4.2/Ex4_2.sce
@@ -0,0 +1,12 @@
+//Chapter 4,Example 4.2 Page 145
+clc
+clear
+// Determine the potential difference
+AP = 8^2/4 // Area of plate
+d = 4 // mm
+FA = 0.2*9.8*10^-3 // Force of attraction
+V = sqrt(FA*2*36*16*10^-6/(10^-9*16*10^-4))
+printf (" Area of plate =  %f π sq.cm \n ",AP)
+printf (" V = %f V \n ",V)
+
+//Answers may vary due to round off error
diff --git a/3487/CH4/EX4.3/Ex4_3.sce b/3487/CH4/EX4.3/Ex4_3.sce
new file mode 100644
index 000000000..d22e042dd
--- /dev/null
+++ b/3487/CH4/EX4.3/Ex4_3.sce
@@ -0,0 +1,16 @@
+//Chapter 4,Example 4.3 Page 145
+clc
+clear
+d = 1 //mm
+V = 10^3 // V
+F = 5*10^-3 // pull between the plates in N
+E = 1/(36) // epselon
+A = 10^2/4 // Area of the plate 
+d1 = sqrt((1/(2*F))*E*10^-9*V^2*A*10^-4) // calculation done in the text is wrong
+d21 = 1/(d1*10^4)
+d22 = 1/(d1*10^4+d)
+C = (V*E*10^-9*A*10^-4)*(d21-d22)
+printf (" d = %f mm \n ",d1*10^4)
+printf (" charge in capacitance = %f pF \n ",C*10^12)
+
+//Answers may vary due to round off error
diff --git a/3487/CH4/EX4.4/Ex4_4.sce b/3487/CH4/EX4.4/Ex4_4.sce
new file mode 100644
index 000000000..80341e37d
--- /dev/null
+++ b/3487/CH4/EX4.4/Ex4_4.sce
@@ -0,0 +1,13 @@
+//Chapter 4,Example 4.4 Page 145
+clc
+clear
+Imin = 2*10^-6 // A
+Imax = 35*10^-6 // A
+V = 15*10^4 // V
+w = 2*%pi*1500/60
+Cm = sqrt(2)*Imin/(V*w) 
+Ipeak = 2*250/15
+printf (" Cm = %f pF \n ",Cm*10^13)
+printf (" At 250 kV, the current indicated will be = %f μA \n ",Ipeak)
+
+//Answers vary due to round off error
diff --git a/3487/CH4/EX4.5/Ex4_5.sce b/3487/CH4/EX4.5/Ex4_5.sce
new file mode 100644
index 000000000..8750d4122
--- /dev/null
+++ b/3487/CH4/EX4.5/Ex4_5.sce
@@ -0,0 +1,15 @@
+//Chapter 4,Example 4.5 Page 146
+clc
+clear
+V1 = 150*10^3 // V
+PD = 1200 // potential divider ratio
+I = 10^-6 // A
+t = 8 // sec
+V = V1/PD
+R = V/I
+C = t*10^6/R
+printf (" V = %f V \n ",V)
+printf (" R = %f M Ω \n ",R*10^-6)
+printf (" C = %f μF \n ",C)
+
+// Answers provided in the textbook are wrong
diff --git a/3487/CH4/EX4.6/Ex4_6.sce b/3487/CH4/EX4.6/Ex4_6.sce
new file mode 100644
index 000000000..fc60543fd
--- /dev/null
+++ b/3487/CH4/EX4.6/Ex4_6.sce
@@ -0,0 +1,23 @@
+//Chapter 4,Example 4.6 Page 146
+clc
+clear
+i = 8*10^3 // i(t)
+V0 = 8 // V0(t)
+I = 8*10^3 // A
+rcI = 10^10 // rate of change of current in A/sec
+R = 8*10^3 // ohm
+RCbyM = i/V0 // R*C/M
+t = I/rcI // 1/4 of cycle
+T = t*4 
+f = 1/T
+CR = 5/f
+M = CR/RCbyM
+C = CR/R
+printf (" Time for 1/4 cycle = %e sec \n ",t)
+printf (" Full time = %e sec \n ",T)
+printf (" f = %e Hz \n ",f)
+printf (" M = %e H \n ",M)
+printf (" C = %e F \n ",C)
+printf (" R = %e Ω \n ",R)
+
+// Answers may vary due to round off error