1 files changed, 46 insertions, 0 deletions
diff --git a/3793/CH9/EX9.3/exp_9_3.sce b/3793/CH9/EX9.3/exp_9_3.sce
new file mode 100644
index 000000000..e1157047f
--- /dev/null
+++ b/3793/CH9/EX9.3/exp_9_3.sce
@@ -0,0 +1,46 @@
+clear;
+clc;
+function fault3faze (Zbus,nfbuses,loc,elemz,col,locs,Zf);
+    for n=1:nfbuses
+        p=input('number of bus to be faulted');
+        Vf=input('fault bus voltage');
+        If=Vf/(Zbus(p,p)+Zf);
+        mprintf("Bus no.         Fault current\n");
+        mprintf("%2i\",p);
+        mprintf("         %15.4f\",real(If));
+        mprintf("         %15.4f\n",imag(If));
+        for k=1:3
+            V(k)=Vf-Zbus(k,p)*If;
+             mprintf("Bus no.         Bus Voltage\n");
+        mprintf("%2i\",k);
+        mprintf("         %15.4f\",real(V(k)));
+        mprintf("         %15.4f\n",imag(V(k)));
+    end
+    kk=1;
+    for k=1:locs-1
+        add=loc(k+1)-loc(k);
+        for m=1:add
+            j=col(kk);
+            I(k,j)=(V(k)-V(j))/elemz(kk);
+            kk=kk+1;
+            mprintf("Bus no.       Bus No.     Current\n");
+        mprintf("%2i\",k);
+         mprintf("......%10i\",j);
+        mprintf(" %15.4f\",real(I(k,j)));
+        mprintf(" %15.4f\n",imag(I(k,j)));
+    end
+    end
+            
+            
+        
+    end
+endfunction
+nfbuses=2;
+loc=[1 3 4];
+elemz=[.2 .3 .25]*%i;
+col=[2 3 3];
+locs=3;
+Zf=%i*.2;
+Zbus=[.0776 .0448 .0597;.0448 .1104 .0806;.0597 .0806 .2075]*%i;
+
+fault3faze (Zbus,nfbuses,loc,elemz,col,locs,Zf);