summaryrefslogtreecommitdiff
path: root/2078/CH4
diff options
context:
space:
mode:
Diffstat (limited to '2078/CH4')
-rwxr-xr-x2078/CH4/EX4.1/Example4_1.sce13
-rwxr-xr-x2078/CH4/EX4.10/Example4_10.sce12
-rwxr-xr-x2078/CH4/EX4.11/Example4_11.sce10
-rwxr-xr-x2078/CH4/EX4.12/Example4_12.sce27
-rwxr-xr-x2078/CH4/EX4.13/Example4_13.sce37
-rwxr-xr-x2078/CH4/EX4.14/Example4_14.sce24
-rwxr-xr-x2078/CH4/EX4.15/Example4_15.sce34
-rwxr-xr-x2078/CH4/EX4.16/Example4_16.sce11
-rwxr-xr-x2078/CH4/EX4.17/Example4_17.sce11
-rwxr-xr-x2078/CH4/EX4.18/Example4_18.sce12
-rwxr-xr-x2078/CH4/EX4.19/Example4_19.sce17
-rwxr-xr-x2078/CH4/EX4.2/Example4_2.sce13
-rwxr-xr-x2078/CH4/EX4.20/Example4_20.sce17
-rwxr-xr-x2078/CH4/EX4.21/Example4_21.sce31
-rwxr-xr-x2078/CH4/EX4.22/Example4_22.sce13
-rwxr-xr-x2078/CH4/EX4.23/Example4_23.sce32
-rwxr-xr-x2078/CH4/EX4.24/Example4_24.sce36
-rwxr-xr-x2078/CH4/EX4.25/Example4_25.sce11
-rwxr-xr-x2078/CH4/EX4.26/Example4_26.sce21
-rwxr-xr-x2078/CH4/EX4.3/Example4_3.sce13
-rwxr-xr-x2078/CH4/EX4.4/Example4_4.sce24
-rwxr-xr-x2078/CH4/EX4.5/Example4_5.sce20
-rwxr-xr-x2078/CH4/EX4.6/Example4_6.sce20
-rwxr-xr-x2078/CH4/EX4.7/Example4_7.sce10
-rwxr-xr-x2078/CH4/EX4.8/Example4_8.sce12
-rwxr-xr-x2078/CH4/EX4.9/Example4_9.sce15
26 files changed, 496 insertions, 0 deletions
diff --git a/2078/CH4/EX4.1/Example4_1.sce b/2078/CH4/EX4.1/Example4_1.sce
new file mode 100755
index 000000000..8d1bdbec4
--- /dev/null
+++ b/2078/CH4/EX4.1/Example4_1.sce
@@ -0,0 +1,13 @@
+//Exa 4.1
+clc;
+clear;
+close;
+//Given data :
+f=50;//Hz
+d=1*100;//cm
+r=1.25/2;//cm
+r_dash=r*0.7788;//cm
+L=0.4*log(d/r_dash);//mH
+disp(L,"Loop inductance per km(mH)");
+XL=2*%pi*f*L*10^-3;//ohm/Km
+disp(XL,"Reactance of transmission line(ohm/km)");
diff --git a/2078/CH4/EX4.10/Example4_10.sce b/2078/CH4/EX4.10/Example4_10.sce
new file mode 100755
index 000000000..73586bed8
--- /dev/null
+++ b/2078/CH4/EX4.10/Example4_10.sce
@@ -0,0 +1,12 @@
+//Exa 4.10
+clc;
+clear;
+close;
+//Given data :
+r=1.2/2*10;//mm
+rdash=0.7788*r;//mm
+d=3.5*1000;//mm
+L=2*10^-7*log(d/rdash);//H/m
+Lav=1/3*(L+L+L);//H/m
+d=rdash*exp(Lav/(2*10^-7)-1/3*log(2));//mm
+disp(d/1000,"Spacing between adjacent conductors(m)");
diff --git a/2078/CH4/EX4.11/Example4_11.sce b/2078/CH4/EX4.11/Example4_11.sce
new file mode 100755
index 000000000..a4f6f4681
--- /dev/null
+++ b/2078/CH4/EX4.11/Example4_11.sce
@@ -0,0 +1,10 @@
+//Exa 4.11
+clc;
+clear;
+close;
+//Given data :
+r=20;//mm
+rdash=0.7788*r;//mm
+d=7*1000;//mm
+L=10^-7*log(sqrt(3)/2*d/rdash);//H/m
+disp(L*10^3/10^-3,"Spacing between adjacent conductors(mH)");
diff --git a/2078/CH4/EX4.12/Example4_12.sce b/2078/CH4/EX4.12/Example4_12.sce
new file mode 100755
index 000000000..6e1c4e325
--- /dev/null
+++ b/2078/CH4/EX4.12/Example4_12.sce
@@ -0,0 +1,27 @@
+//Exa 4.12
+clc;
+clear;
+close;
+//Given data :
+r=0.9;//cm
+rdash=0.7788*r*10^-2;//m
+daa_dash=sqrt(6^2+6^2);//m
+dbb_dash=7;//m
+dcc_dash=daa_dash;//m
+daa=rdash;//m
+d_adash_adash=rdash;//m
+d_adash_a=daa_dash;//m
+Dsa=(daa*daa_dash*d_adash_adash*d_adash_a)^(1/4);//m
+Dsb=(daa*7)^(1/2);//m
+Dsc=(daa*daa_dash)^(1/2);//m
+Ds=(Dsa*Dsb*Dsc)^(1/3);//m
+dab=sqrt(3^2+0.5^2);//m
+dab_dash=sqrt(3^2+6.5^2);//m
+d_adash_b=sqrt(3^2+6.5^2);//m
+d_adash_bdash=sqrt(3^2+0.5^2);//m
+Dab=(dab*dab_dash*d_adash_b*d_adash_bdash)^(1/4);//m
+Dbc=((dab*dab_dash)^2)^(1/4);//m
+Dca=((6*6)^2)^(1/4);//m
+Dm=(Dab*Dbc*Dca)^(1/3);//m
+L=0.2*log(Dm/Ds);//mH/km
+disp(L,"Inductance per phase(mH/km)");
diff --git a/2078/CH4/EX4.13/Example4_13.sce b/2078/CH4/EX4.13/Example4_13.sce
new file mode 100755
index 000000000..def868d03
--- /dev/null
+++ b/2078/CH4/EX4.13/Example4_13.sce
@@ -0,0 +1,37 @@
+//Exa 4.13
+clc;
+clear;
+close;
+format('v',5)
+//Given data :
+r=5/2;//mm
+rdash=2.176*r*10^-3;//m
+daa_dash=sqrt(6^2+16^2);//m
+dbb_dash=6;//m
+dcc_dash=daa_dash;//m
+dab=8;//m
+dab_dash=sqrt(6^2+8^2);//m
+dbc=8;//m
+dbc_dash=sqrt(6^2+8^2);//m
+dca=16;//m
+dca_dash=6;//m
+Dsa=sqrt(rdash*daa_dash);//m
+Dsb=sqrt(rdash*dbb_dash);//m
+Dsc=sqrt(rdash*dcc_dash);//m
+Ds=(Dsa*Dsb*Dsc)^(1/3);//m
+disp(Ds,"GMD(m) : ");
+Dab=(dab*dab_dash)^(1/2);//m
+Dbc=(dbc*dbc_dash)^(1/2);//m
+Dca=(dca*dca_dash)^(1/2);//m
+Dm=(Dab*Dbc*Dca)^(1/3);//m
+disp(Dm,"Deq or Dm(m) : ");
+L=0.2*log(Dm/Ds);//mH/km
+L=L*10^-3*100;//H(for 100 km line)
+disp(L,"Inductance of 100 km line(H)");
+///Alternate method is given below
+d1=dab;//m
+d2=dca_dash;//m
+L=0.2*log(2^(1/6))*sqrt(d1/rdash)*((d1^2+d2^2)/(4*d1^2+d2^2))^(1/6);//mH
+L=L*10^-3*100;//H(for 100 km line)
+disp(L,"Using Alternate method, Inductance of 100 km line(H)");
+
diff --git a/2078/CH4/EX4.14/Example4_14.sce b/2078/CH4/EX4.14/Example4_14.sce
new file mode 100755
index 000000000..edfcd553f
--- /dev/null
+++ b/2078/CH4/EX4.14/Example4_14.sce
@@ -0,0 +1,24 @@
+//Exa 4.14
+clc;
+clear;
+close;
+//Given data :
+r=5/2;//cm
+rdash=0.7788*r*10^-2;//m
+d=6.5;//m
+s=0.4;//m
+Ds=sqrt(rdash*s);//m
+dab=6.5;//m
+dab_dash=6.9;//m
+d_adash_b=6.1;//m
+d_adash_bdash=6.5;//m
+Dab=(dab*dab_dash*d_adash_b*d_adash_bdash)^(1/4);//m
+Dbc=Dab;//m
+dca=13;//m
+dca_dash=12.6;//m
+d_cdash_a=13.4;//m
+d_cdash_adash=13;//m
+Dca=(dca*dca_dash*d_cdash_a*d_cdash_adash)^(1/4);//m
+Dm=(Dab*Dbc*Dca)^(1/3);//m
+L=0.2*log(Dm/Ds);//mH/km
+disp(L,"Inductance per phase(mH/km)");
diff --git a/2078/CH4/EX4.15/Example4_15.sce b/2078/CH4/EX4.15/Example4_15.sce
new file mode 100755
index 000000000..ac8ecb33f
--- /dev/null
+++ b/2078/CH4/EX4.15/Example4_15.sce
@@ -0,0 +1,34 @@
+//Exa 4.15
+clc;
+clear;
+close;
+//Given data :
+f=50;//Hz
+r=3.5/2;//cm
+rdash=0.7788*r*10^-2;//m
+d=7;//m
+s=40/100;//m
+Ds=sqrt(rdash*s);//m
+dab=7;//m
+dab_dash=7.4;//m
+d_adash_b=6.6;//m
+d_adash_bdash=7;//m
+Dab=(dab*dab_dash*d_adash_b*d_adash_bdash)^(1/4);//m
+Dbc=Dab;//m
+dca=14;//m
+dca_dash=13.6;//m
+d_cdash_a=14.4;//m
+d_cdash_adash=14;//m
+Dca=(dca*dca_dash*d_cdash_a*d_cdash_adash)^(1/4);//m
+Dm=(Dab*Dbc*Dca)^(1/3);//m
+L=0.2*log(Dm/Ds);//mH/km
+XL=2*%pi*f*L*10^-3;//ohm/km
+disp(XL,"Inductive reactance of bundled conductor line(ohm/km)");
+//Equivalent single conductor
+n=2;
+r1=sqrt(n*%pi*r^2/%pi);//m
+r1dash=0.7788*r1*10^-2;//m
+Dm1=(Dab*Dbc*Dca)^(1/3);//m
+L1=0.2*log(Dm1/r1dash);//mH/km
+XL1=2*%pi*f*L1*10^-3;//ohm/km
+disp(XL1,"Inductive reactance with single conductor(ohm/km)");
diff --git a/2078/CH4/EX4.16/Example4_16.sce b/2078/CH4/EX4.16/Example4_16.sce
new file mode 100755
index 000000000..185da699b
--- /dev/null
+++ b/2078/CH4/EX4.16/Example4_16.sce
@@ -0,0 +1,11 @@
+//Exa 4.16
+clc;
+clear;
+close;
+//Given data :
+r=15/2;//mm
+d=1.5*1000;//mm
+l=30;//km
+epsilon_o=8.854*10^-12;//permitivity
+C=%pi*epsilon_o/log(d/r)*l*1000;//F
+disp(C*10^6,"Capacitance of line(micro F)");
diff --git a/2078/CH4/EX4.17/Example4_17.sce b/2078/CH4/EX4.17/Example4_17.sce
new file mode 100755
index 000000000..7fa855846
--- /dev/null
+++ b/2078/CH4/EX4.17/Example4_17.sce
@@ -0,0 +1,11 @@
+//Exa 4.17
+clc;
+clear;
+close;
+//Given data :
+r=2/2;//cm
+d=2.5*100;//cm
+l=100;//km
+epsilon_o=8.854*10^-12;//permitivity
+C=2*%pi*epsilon_o/log(d/r)*l*1000;//F
+disp(C*10^6,"Capacitance of line(micro F)");
diff --git a/2078/CH4/EX4.18/Example4_18.sce b/2078/CH4/EX4.18/Example4_18.sce
new file mode 100755
index 000000000..ac3ad613d
--- /dev/null
+++ b/2078/CH4/EX4.18/Example4_18.sce
@@ -0,0 +1,12 @@
+//Exa 4.18
+clc;
+clear;
+close;
+//Given data :
+r=2/2/100;//m
+d1=3.5;//m
+d2=5;//m
+d3=8;//m
+epsilon_o=8.854*10^-12;//permitivity
+CN=2*%pi*epsilon_o*1000/log((d1*d2*d3)^(1/3)/r);//F
+disp(CN*10^6,"Capacitance of line(micro F)");
diff --git a/2078/CH4/EX4.19/Example4_19.sce b/2078/CH4/EX4.19/Example4_19.sce
new file mode 100755
index 000000000..39b46b9fe
--- /dev/null
+++ b/2078/CH4/EX4.19/Example4_19.sce
@@ -0,0 +1,17 @@
+//Exa 4.19
+clc;
+clear;
+close;
+//Given data :
+f=50;//Hz
+VL=220;//KV
+r=20/2/1000;//m
+d1=3;//m
+d2=3;//m
+d3=6;//m
+epsilon_o=8.854*10^-12;//permitivity
+CN=2*%pi*epsilon_o/log((d1*d2*d3)^(1/3)/r);//F
+disp(CN,"Capacitance per phase per meter line(F)");
+Vph=VL*1000/sqrt(3);//V
+Ic=2*%pi*f*CN*Vph;//A
+disp(Ic*1000,"Charging current per phase(mA) : ");
diff --git a/2078/CH4/EX4.2/Example4_2.sce b/2078/CH4/EX4.2/Example4_2.sce
new file mode 100755
index 000000000..8c7207da9
--- /dev/null
+++ b/2078/CH4/EX4.2/Example4_2.sce
@@ -0,0 +1,13 @@
+//Exa 4.2
+clc;
+clear;
+close;
+//Given data :
+f=50;//Hz
+a=10;//cm^2
+l=500/1000;//km
+r=sqrt(a/%pi);//cm
+d=5*100;//cm
+r_dash=r*0.7788;//cm
+L=0.4*log(d/r_dash)*l;//mH
+disp(L,"Loop inductance per km(mH)");
diff --git a/2078/CH4/EX4.20/Example4_20.sce b/2078/CH4/EX4.20/Example4_20.sce
new file mode 100755
index 000000000..83937d30b
--- /dev/null
+++ b/2078/CH4/EX4.20/Example4_20.sce
@@ -0,0 +1,17 @@
+//Exa 4.20
+clc;
+clear;
+close;
+//Given data :
+f=50;//Hz
+VL=110;//kV
+r=1.05/2;//cm
+d1=3.5;//m
+d2=3.5;//m
+d3=7;//m
+epsilon_o=8.854*10^-12;//permitivity
+CN=2*%pi*epsilon_o/log((d1*d2*d3)^(1/3)*100/r);//F
+disp(CN,"Capacitance per phase per meter line(F)");
+Vph=VL*1000/sqrt(3);//V
+Ic=2*%pi*f*CN*Vph;//A/m
+disp(Ic/10^-3,"Charging current per phase(A/km) : ");
diff --git a/2078/CH4/EX4.21/Example4_21.sce b/2078/CH4/EX4.21/Example4_21.sce
new file mode 100755
index 000000000..31957a768
--- /dev/null
+++ b/2078/CH4/EX4.21/Example4_21.sce
@@ -0,0 +1,31 @@
+//Exa 4.21
+clc;
+clear;
+close;
+//Given data :
+r=2.5/2*10^-2;//m
+VL=132;//KV
+epsilon_o=8.85*10^-12;//permitivity
+f=50;//Hz
+dRRdash=sqrt(7^2+(4+4)^2);//m
+dBBdash=dRRdash;//m
+dYYdash=9;//m
+DSR=sqrt(r*dRRdash);//m
+DSY=sqrt(r*dYYdash);//m
+DSB=sqrt(r*dBBdash);//m
+Ds=(DSR*DSB*DSY)^(1/3);//m
+dRY=sqrt(4^2+(4.5-3.5)^2);//m
+dRYdash=sqrt((9-1)^2+4^2);//m
+dRdashY=sqrt((9-1)^2+4^2);//m
+dRdashYdash=sqrt(4^2+(4.5-3.5)^2);//m
+DRY=(dRY*dRYdash*dRdashY*dRdashYdash)^(1/4);//m
+DYB=((dRY*dRYdash)^2)^(1/4);//m
+DBR=((8*7)^2)^(1/4);//m
+Dm=(DRY*DYB*DBR)^(1/3);//m
+C=2*%pi*epsilon_o/log(Dm/Ds);//F/m
+C=C/10^-3;//F/km
+X=1/(2*%pi*f*C);//ohm
+disp(X/1000,"Capacitive reactance too neutral(kohm) : ");
+Vph=VL*1000/sqrt(3);//Volt
+Ic=2*%pi*f*C*Vph;//A
+disp(Ic,"Charging current(A/km)");
diff --git a/2078/CH4/EX4.22/Example4_22.sce b/2078/CH4/EX4.22/Example4_22.sce
new file mode 100755
index 000000000..1dd15494b
--- /dev/null
+++ b/2078/CH4/EX4.22/Example4_22.sce
@@ -0,0 +1,13 @@
+//Exa 4.22
+clc;
+clear;
+close;
+//Given data :
+d1=8;//m
+d2=6;//m
+epsilon_o=8.854*10^-12;//permitivity
+r=3*5/2*10^-3;//m
+C=4*%pi*epsilon_o/log(2^(1/3)*d1/r*((d1^2+d2^2)/(4*d1^2+d2^2)^(1/3)));//F/m
+C100=C*100*1000*10^6;//microF
+disp(C100,"Capacitance of 100 km line(micro Farad) : ");
+//answer in the textbook is wrong.
diff --git a/2078/CH4/EX4.23/Example4_23.sce b/2078/CH4/EX4.23/Example4_23.sce
new file mode 100755
index 000000000..fdb3e36fd
--- /dev/null
+++ b/2078/CH4/EX4.23/Example4_23.sce
@@ -0,0 +1,32 @@
+//Exa 4.23
+clc;
+clear;
+close;
+//Given data :
+VL=132;//kV
+f=50;//Hz
+r=5/2;//cm
+rdash=0.7788*r*10^-2;//m
+d=6.5;//m
+s=0.4;//m
+epsilon_o=8.854*10^-12;//permitivity
+Ds=sqrt(rdash*s);//m
+dab=6.5;//m
+dab_dash=6.9;//m
+d_adash_b=6.1;//m
+d_adash_bdash=6.5;//m
+Dab=(dab*dab_dash*d_adash_b*d_adash_bdash)^(1/4);//m
+Dbc=Dab;//m
+dca=13;//m
+dca_dash=12.6;//m
+d_cdash_a=13.4;//m
+d_cdash_adash=13;//m
+Dca=(dca*dca_dash*d_cdash_a*d_cdash_adash)^(1/4);//m
+Dm=(Dab*Dbc*Dca)^(1/3);//m
+L=0.2*log(Dm/Ds);//mH/km
+C=2*%pi*epsilon_o/log(Dm/Ds);//F/m
+C=C/10^-3;//F/km
+disp(C,"Capacitance per km(F/km) : ");
+Vph=VL*1000/sqrt(3);//Volt
+Ic=2*%pi*f*C*Vph;//A/km
+disp(Ic,"Charging current per km(A/km) : ");
diff --git a/2078/CH4/EX4.24/Example4_24.sce b/2078/CH4/EX4.24/Example4_24.sce
new file mode 100755
index 000000000..29177b55a
--- /dev/null
+++ b/2078/CH4/EX4.24/Example4_24.sce
@@ -0,0 +1,36 @@
+//Exa 4.24
+clc;
+clear;
+close;
+//Given data :
+VL=132;//kV
+f=50;//Hz
+r=31.8/2;//mm
+rdash=0.7788*r;//mm
+d=10*1000;//mm
+epsilon_o=8.854*10^-12;//permitivity
+disp("One conductor ACSR moose conductor line : ");
+LA=0.2*[log(d/rdash)+1/2*log(2)-%i*0.866*log(2)];//mH/km
+LB=0.2*log(d/rdash);//mH/km
+LC=0.2*[log(d/rdash)+1/2*log(2)+%i*0.866*log(2)];//mH/km
+Lav=(LA+LB+LC)/3;//mH/km
+XL=2*%pi*f*Lav*10^-3;//ohm
+disp(XL,"Inductive reactance per Km per phase(ohm) : ");
+d1=10;//m
+d2=10;//m
+d3=20;//m
+CN=2*%pi*epsilon_o/log((d1*d2*d3)^(1/3)/(rdash*10^-3))/10^3;//F/km
+XC=1/(2*%pi*f*CN*10^6);//ohm
+disp(XC/10^6,"Capacitivetive reactance per Km per phase(Mohm) : ");
+disp("Three conductor bundled line : ");
+S=40/100;//m
+Ds=(rdash*10^-3*S^2)^(1/3);//m
+Deq=(d1*d2*d3)^(1/3);//m
+Ldash=0.2*log(Deq/Ds);//mH/km
+XLdash=2*%pi*f*Ldash*10^-3;//ohm
+disp(XLdash,"Inductive reactance per km per phase(ohm) : ");
+Ds=(r*10^-3*S^2)^(1/3);//m
+Cdash=2*%pi*epsilon_o*10^3/log(Deq/Ds);//microF/km
+XC=1/(2*%pi*f*Cdash)/10^6;//Mohm
+disp(XC,"Capacitivetive reactance per km per phase(Mohm) : ");
+
diff --git a/2078/CH4/EX4.25/Example4_25.sce b/2078/CH4/EX4.25/Example4_25.sce
new file mode 100755
index 000000000..8c9a271e8
--- /dev/null
+++ b/2078/CH4/EX4.25/Example4_25.sce
@@ -0,0 +1,11 @@
+//Exa 4.25
+clc;
+clear;
+close;
+//Given data :
+r=1.5/2;//cm
+d=3*100;//cm
+h=6*100;//cm
+epsilon_o=8.854*10^-12;//permitivity
+C=%pi*epsilon_o/log(d/(1+d^2/4/h^2)^r)*10^3;//F
+disp(C,"Capacitance per km of line(F) : ");
diff --git a/2078/CH4/EX4.26/Example4_26.sce b/2078/CH4/EX4.26/Example4_26.sce
new file mode 100755
index 000000000..cd04b63b3
--- /dev/null
+++ b/2078/CH4/EX4.26/Example4_26.sce
@@ -0,0 +1,21 @@
+//Exa 4.26
+clc;
+clear;
+close;
+//Given data :
+r=2/100;//m
+d1=4;//m
+d2=4;//m
+d3=8;//m
+epsilon_o=8.854*10^-12;//permitivity
+CN=2*%pi*epsilon_o/log((d1*d2*d3)^(1/3)/r);//F
+disp(CN,"Part(i) Capacitance per phase per meter length(F) : ");
+h1=20;//m
+h2=20;//m
+h3=20;//m
+h12=sqrt(20^2+4^2);//m
+h23=sqrt(20^2+4^2);//m
+h31=sqrt(20^2+8^2);//m
+Deq=(d1*d2*d3)^(1/3);//m
+CN=2*%pi*epsilon_o/(log(Deq/r)-log((h12*h23*h31/h1/h2/h3)^(1/3)) );//F
+disp(CN,"Part(ii) Capacitance per phase per meter length(F) : ");
diff --git a/2078/CH4/EX4.3/Example4_3.sce b/2078/CH4/EX4.3/Example4_3.sce
new file mode 100755
index 000000000..b5f51bd49
--- /dev/null
+++ b/2078/CH4/EX4.3/Example4_3.sce
@@ -0,0 +1,13 @@
+//Exa 4.3
+clc;
+clear;
+close;
+//Given data :
+r=1/2;//cm
+d=1*100;//cm
+mu=50;//relative permeability
+r_dash=r*0.7788;//cm
+L_cu=.1+0.4*log(d/r);//mH
+disp(L_cu,"Loop inductance per km of copper conductor line(mH)");
+L_steel=(mu+4*log(d/r))*10^-7*10^3;//mH
+disp(L_steel*10^3,"Loop inductance per km of copper conductor line(mH)");
diff --git a/2078/CH4/EX4.4/Example4_4.sce b/2078/CH4/EX4.4/Example4_4.sce
new file mode 100755
index 000000000..c92095c46
--- /dev/null
+++ b/2078/CH4/EX4.4/Example4_4.sce
@@ -0,0 +1,24 @@
+//Exa 4.4
+clc;
+clear;
+close;
+//Given data :
+r=3;//mm
+d11=r;//mm
+d12=2*r;//mm
+d34=2*r;//mm
+d16=2*r;//mm
+d17=2*r;//mm
+d14=4*r;//mm
+d13=sqrt(d14^2-d34^2);//mm
+d15=d13;//mm
+Ds1=(0.7788*d11*d12*d13*d14*d15*d16*d17)^(1/7);//mm
+Ds2=Ds1;//mm
+Ds3=Ds1;//mm
+Ds4=Ds1;//mm
+Ds5=Ds1;//mm
+Ds6=Ds1;//mm
+Ds7=(2*r*0.7788*d11*d12*d13*2*r*2*r)^(1/7);//mm
+Ds=(Ds1*Ds2*Ds3*Ds4*Ds5*Ds6*Ds7)^(1/7);//mm
+disp(Ds,"Geometric mean radius(mm)");
+//Answer in the book is wrong
diff --git a/2078/CH4/EX4.5/Example4_5.sce b/2078/CH4/EX4.5/Example4_5.sce
new file mode 100755
index 000000000..cfe710f3b
--- /dev/null
+++ b/2078/CH4/EX4.5/Example4_5.sce
@@ -0,0 +1,20 @@
+//Exa 4.5
+clc;
+clear;
+close;
+//Given data :
+r=1.2;//cm
+rdash=0.7788*r;//cm
+d12=0.12*100;//cm
+d11dash=(0.2+1.2)*100;//cm
+d22dash=(0.2+1.2)*100;//cm
+d12dash=(0.2+1.2+0.2)*100;//cm
+d21dash=(1.2)*100;//cm
+Dm=(d11dash*d12dash*d21dash*d22dash)^(1/4);//cm
+d11=0.93456;//cm
+d22=0.93456;//cm
+d12=20;//cm
+d21=20;//cm
+Ds=(d11*d12*d21*d22)^(1/4);//cm
+L=0.4*log(Dm/Ds);//mH/km
+disp(L,"Loop inductance of line(mH/km)");
diff --git a/2078/CH4/EX4.6/Example4_6.sce b/2078/CH4/EX4.6/Example4_6.sce
new file mode 100755
index 000000000..e006068f6
--- /dev/null
+++ b/2078/CH4/EX4.6/Example4_6.sce
@@ -0,0 +1,20 @@
+//Exa 4.6
+clc;
+clear;
+close;
+//Given data :
+r=2/2;//cm
+rdash=0.7788*r;//cm
+d12=0.12*100;//cm
+d11dash=300;//cm
+d12dash=sqrt(300^2+100^2);//cm
+d21dash=d12dash;//cm
+d22dash=d11dash;//cm
+d11=rdash;//cm
+d22=rdash;//cm
+d12=100;//cm
+d21=100;//cm
+Dm=(d11dash*d12dash*d21dash*d22dash)^(1/4);//cm
+Ds=(d11*d12*d21*d22)^(1/4);//cm
+L=0.4*log(Dm/Ds);//mH/km
+disp(L,"Loop inductance of line(mH/km)");
diff --git a/2078/CH4/EX4.7/Example4_7.sce b/2078/CH4/EX4.7/Example4_7.sce
new file mode 100755
index 000000000..84627f55b
--- /dev/null
+++ b/2078/CH4/EX4.7/Example4_7.sce
@@ -0,0 +1,10 @@
+//Exa 4.7
+clc;
+clear;
+close;
+//Given data :
+r=1.24/2;//cm
+rdash=0.7788*r;//cm
+d=2*100;//cm
+L=0.2*log(d/rdash);//mH
+disp(L,"Inductance per phase per km(mH)");
diff --git a/2078/CH4/EX4.8/Example4_8.sce b/2078/CH4/EX4.8/Example4_8.sce
new file mode 100755
index 000000000..0678e2204
--- /dev/null
+++ b/2078/CH4/EX4.8/Example4_8.sce
@@ -0,0 +1,12 @@
+//Exa 4.8
+clc;
+clear;
+close;
+//Given data :
+r=(20/2)/10;//cm
+d1=4*100;//cm
+d2=5*100;//cm
+d3=6*100;//cm
+rdash=0.7788*r;//cm
+L=0.2*log((d1*d2*d3)^(1/3)/rdash);//mH
+disp(L,"Inductance per phase(mH)");
diff --git a/2078/CH4/EX4.9/Example4_9.sce b/2078/CH4/EX4.9/Example4_9.sce
new file mode 100755
index 000000000..8fe1ad62e
--- /dev/null
+++ b/2078/CH4/EX4.9/Example4_9.sce
@@ -0,0 +1,15 @@
+//Exa 4.9
+clc;
+clear;
+close;
+//Given data :
+r=4/2;//cm
+rdash=0.7788*r;//cm
+d=300;//cm
+d3=6*100;//cm
+LA=0.2*[log(d/rdash)+1/2*log(2)-%i*0.866*log(2)];//mH
+disp(LA,"Inductance per km of phase1(mH)");
+LB=0.2*log(d/rdash);//mH
+disp(LB,"Inductance per km of phase2(mH)");
+LC=0.2*[log(d/rdash)+1/2*log(2)+%i*0.866*log(2)];//mH
+disp(LC,"Inductance per km of phase3(mH)");