1 files changed, 33 insertions, 0 deletions
diff --git a/3793/CH10/EX10.5/exp_10_5.sce b/3793/CH10/EX10.5/exp_10_5.sce
new file mode 100644
index 000000000..62c8e4af8
--- /dev/null
+++ b/3793/CH10/EX10.5/exp_10_5.sce
@@ -0,0 +1,33 @@
+clear;
+clc;
+s=30;
+v=11;
+a=complex(-.5,.866);
+A=[1 1 1;1 a^2 a;1 a a^2];
+xs=complex(0,.5);
+xn=complex(0,.3);
+xz=complex(0,.08);
+Z=5;
+X=(Z*s)/v^2;
+Xn=complex(0,3*X);
+Ea=complex(1,0);
+Ia1=Ea/(Xn+xn+xs+xz);
+Ia=3*Ia1;
+Va0=-Ia1*(xz+Xn);
+Va1=Ea-(Ia1*xs);
+Va2=-Ia1*xn;
+V1=[Va0;Va1;Va2];
+V=A*V1;
+Vab=V(1,1)-V(2,1);
+Vbc=V(2,1)-V(3,1);
+Vca=V(3,1)-V(1,1);
+Ibase=(s*1000)/(sqrt(3)*v);
+If=Ia1*Ibase;
+Vab1=Vab*(v/sqrt(3));
+Vbc1=Vbc*(v/sqrt(3));
+Vca1=Vca*(v/sqrt(3));
+mprintf("Sub transient fault current is %.3f A \n",imag(If));
+mprintf("Actual Line Voltages are Vab=%f+%f Kv   Vbc=%f%f Kv  Vca=%f+%f Kv",real(Vab1),imag(Vab1),real(Vbc1),imag(Vbc1),real(Vca1),imag(Vca1));
+
+
+