1 files changed, 48 insertions, 0 deletions
diff --git a/3872/CH4/EX4.2/Ex4_2.sce b/3872/CH4/EX4.2/Ex4_2.sce
new file mode 100644
index 000000000..090b22775
--- /dev/null
+++ b/3872/CH4/EX4.2/Ex4_2.sce
@@ -0,0 +1,48 @@
+// Book - Power System: Analysis & Design 5th Edition
+// Authors - J. Duncan Glover, Mulukutla S. Sharma, Thomas J. Overbye
+// Chapter - 4 : Example 4.2
+// Scilab Version 6.0.0 : OS - Windows
+
+clc;
+clear;
+
+rx=0.03;                     //Radius of conductor x in meter
+ry=0.04;                     //Radius of conductor y in meter
+N=3;
+M=2;
+
+Ddash=[4 4.3;3.5 3.8;2 2.3];               //Equivalent distances in meter to find Dxy 
+Dx=[exp(-1/4)*rx 0.5 2;0.5 exp(-1/4)*rx 1.5;2 1.5 exp(-1/4)*rx];       //Equivalent distances in meter to find Dxx
+Dy=[exp(-1/4)*ry 0.3;0.3 exp(-1/4)*ry];                                //Equivalent distances in meter to find Dyy
+
+Dxyr=1;
+for i=1:N
+    for j=1:M
+        Dxyr=Dxyr*Ddash(i,j)
+    end
+end
+Dxy=nthroot(Dxyr,M*N);
+
+Dxxr=1
+for i=1:N
+    for j=1:N
+        Dxxr=Dxxr*Dx(i,j)
+    end
+end
+Dxx=nthroot(Dxxr,N*N);
+
+Dyyr=1
+for i=1:M
+    for j=1:M
+        Dyyr=Dyyr*Dy(i,j)
+    end
+end
+Dyy=nthroot(Dyyr,M*M);
+
+Lx=2e-7*log(Dxy/Dxx);
+Ly=2e-7*log(Dxy/Dyy);
+L=Lx+Ly;
+
+printf('The value of inductance in conductor x is, Lx=%3.2e H/m per conductor\n',Lx)
+printf('The value of inductance in conductor y is, Ly=%3.2e H/m per conductor\n',Ly)
+printf('The value of total inductance is, L=%3.2e H/m per circuit',L)