1 files changed, 53 insertions, 0 deletions
diff --git a/1319/CH12/EX12.19/i_19.sce b/1319/CH12/EX12.19/i_19.sce
new file mode 100644
index 000000000..b113ae9f9
--- /dev/null
+++ b/1319/CH12/EX12.19/i_19.sce
@@ -0,0 +1,53 @@
+// To calculate current in each branch using loop analysis.
+
+clc;
+clear;
+
+// MESH Equations for the given network.
+//3*i1-i2+0*i3=11
+//-i1+10*i2-2*i3=0
+//0*i1+-2*i2+5*i3=13
+
+//Voltage supplies are 11V and 13V
+
+R=[3 -1 0;-1 10 -2; 0 -2 5];
+E=[11;0;13];
+
+// Loop Currents
+
+I=inv(R)*E;
+
+i1=I(1);
+i2=I(2);
+i3=I(3);
+
+ia=i1-i2; // Assumed direction from Mesh 1
+ib=i2-i3; // Assumed direction from Mesh 2
+
+disp('A',ib,'ib (through 2 resistor between 7 ohm and 3 ohm resistor) =','A',ia,'ia(through 1 ohm resistor) =','A',i3,'i3 =','A',i2,'i2 =','A',i1,'i1 =','The Calculated Loop Currents are')
+
+disp('The Negative sign indicates that the assumed direction of flow of current should be reveresed')
+
+// To obtain the magnitude of direction.
+
+if(i1<0)
+    i1=abs(i1);
+end
+if(i2<0)
+    i2=abs(i2);
+end
+if(i3<0)
+    i3=abs(i3);
+end
+if(ia<0)
+    ia=abs(ia);
+end
+if(ib<0)
+    ib=abs(ib);
+end
+
+disp('A',i1,'The Current through 2 ohm resistor on the 11V side =')
+disp('A',i2,'The Current through 7 ohm resistor =')
+disp('A',i3,'The Current through 3 ohm resistor on the 13V side =')
+disp('A',ia,'The Current through 1 ohm resistor =')
+disp('A',ib,'The Current through 2 ohm resistor between the 7 and 3 ohm resistors =')