From b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b Mon Sep 17 00:00:00 2001
From: priyanka
Date: Wed, 24 Jun 2015 15:03:17 +0530
Subject: initial commit / add all books

---
 1076/CH10/EX10.16/10_16.sce | 58 +++++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 58 insertions(+)
 create mode 100755 1076/CH10/EX10.16/10_16.sce

(limited to '1076/CH10/EX10.16/10_16.sce')

diff --git a/1076/CH10/EX10.16/10_16.sce b/1076/CH10/EX10.16/10_16.sce
new file mode 100755
index 000000000..1ce736714
--- /dev/null
+++ b/1076/CH10/EX10.16/10_16.sce
@@ -0,0 +1,58 @@
+clear
+clc
+
+X=[
+ .25 .25 .05   
+ .2 .2 .05
+ .06 .06 .06
+ .07 .07 .07
+ .1 .1 .3
+ .1 .1 .3
+]
+
+B=[
+    100 11
+    100 11
+    100 11
+    100 11
+    100 220 
+    100 220
+    ]
+V1=11
+V2=220
+S=100
+Xe=3*.03
+//end 9
+X1=(((X(1,1)*B(1,1) *V1/(S*B(1,2))) +(X(3,1)*B(3,1) *V1/(S*B(3,2))) )^-1+((X(2,1)*B(2,1) *V1/(S*B(2,2)))+(X(4,1)*B(4,1) *V1/(S*B(4,2))) +((X(5,1)*B(5,1) *V2/(S*B(5,2)))^-1 +(X(6,1)*B(6,1) *V2/(S*B(6,2)))^-1)^-1)^-1)^-1
+
+X2=(((X(1,2)*B(1,1) *V1/(S*B(1,2))) +(X(3,2)*B(3,1) *V1/(S*B(3,2))) )^-1+((X(2,2)*B(2,1) *V1/(S*B(2,2)))+(X(4,2)*B(4,1) *V1/(S*B(4,2))) +((X(5,2)*B(5,1) *V2/(S*B(5,2)))^-1 +(X(6,2)*B(6,1) *V2/(S*B(6,2)))^-1)^-1)^-1)^-1
+
+X0=((X(3,3)*B(3,1) *V1/(S*B(3,2)))^-1 + ((Xe    *B(4,1) *V1/(S*B(4,2))) + (X(2,3)*B(2,1) *V1/(S*B(2,2))) + (X(4,3)*B(4,1) *V1/(S*B(4,2))) +(((X(5,3)*B(5,1) *V2/(S*B(5,2))))^-1 +(X(6,3)*B(6,1) *V2/(S*B(6,2)))^-1)^-1 )^-1 )^-1
+
+Z1=%i * round(X1*1e3)/1e3
+Z2=%i * round(X2*1e3)/1e3
+Z0=%i * round(X0*1e3)/1e3
+
+//end q11
+
+a=exp(%i * 2 * %pi /3)
+A= [1 1 1 ; 1 a^2 a; 1 a a^2]
+
+
+Ia1=1/(Z1 + (Z2*Z0/(Z2+Z0)))
+Ia1=round(Ia1 *1e3)/1e3
+Ia2=(Z0/(Z2+Z0)) * Ia1 *-1
+Ia2=round(Ia2 *1e3)/1e3
+Ia0=(Z2/(Z2+Z0)) * Ia1 *-1
+Ia0=round(Ia0 *1e3)/1e3
+
+mIa1=[ Ia0 Ia1 Ia2]'
+mIa=A*mIa1
+Ib=round(S*1e7/(sqrt(3) * V2 * 1e3))/10
+Ia=round(abs(mIa)) * Ib
+
+
+mprintf("\nline currents are: in Amperes ")
+mprintf("\nIa= %s", string(round(abs(Ia(1))*10)/10)  )
+mprintf("\nIb= %s", string(round(abs(Ia(2))*10)/10) +'/_'+ string(round(atand(imag(mIa(2)*Ib)/real(mIa(2)*Ib))*100)/100 +180 ))
+mprintf("\nIc= %s", string(round(abs(Ia(3))*10)/10) +'/_'+ string(round(atand(imag(mIa(3)*Ib)/real(mIa(3)*Ib))*100)/100  ))
-- 
cgit