add books

1 files changed, 2270 insertions, 0 deletions

diff --git a/sample_notebooks/JigneshBhadani/ch5.ipynb b/sample_notebooks/JigneshBhadani/ch5.ipynb
new file mode 100644
index 00000000..189e1172
--- /dev/null
+++ b/sample_notebooks/JigneshBhadani/ch5.ipynb
@@ -0,0 +1,2270 @@
+{
+ "metadata": {
+  "name": "",
+  "signature": "sha256:f2a55784eb064c2602de5db3699903ee6b51018dcf324e42007dffa21eb1e0ba"
+ },
+ "nbformat": 3,
+ "nbformat_minor": 0,
+ "worksheets": [
+  {
+   "cells": [
+    {
+     "cell_type": "heading",
+     "level": 1,
+     "metadata": {},
+     "source": [
+      "Chapter 5 : SINGLE PHASE AC CIRCUITS"
+     ]
+    },
+    {
+     "cell_type": "heading",
+     "level": 2,
+     "metadata": {},
+     "source": [
+      "Example 5.1  Page No : 157"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "import math \n",
+      "from numpy import ones\n",
+      "\n",
+      "def  r2p(x,y):\t\t\t#function to convert recmath.tangular to polar\n",
+      "    polar  =  ones(2)\n",
+      "    polar[0]  =  math.sqrt ((x **2) +(y**2))\n",
+      "    polar[1]  =  math.atan (y/x)\n",
+      "    polar[1]  = (polar [1]*180)/math.pi\n",
+      "    return polar\n",
+      "    \n",
+      "def  p2r(r,theta):\t\t\t#function to convert polar to recmath.tangular\n",
+      "    rect  =  ones(2)\n",
+      "    theta  = ( theta *math.pi) /180\n",
+      "    rect [0] = r* math.cos(theta)\n",
+      "    rect [1] = r* math.sin(theta)\n",
+      "    return rect\n",
+      "\n",
+      "#CALCULATIONS\n",
+      "I1 = r2p(7,-5);\n",
+      "print (I1);\n",
+      "I2 = r2p(-9,6);\n",
+      "I2[1] = I2[1]+(180);\t\t\t#this belongs to quadrant 2 and hence 180 degrees should be added\n",
+      "print (I2);\n",
+      "I3 = r2p(-8,-8);\n",
+      "I3[1] = I3[1]+(180);\t\t\t#this belongs to quadrant 3 and hence 180 degrees should be added\n",
+      "print (I3);\n",
+      "I4 = r2p(6,6);\n",
+      "print (I4);\n",
+      "#note:here direct functions for converson are not available and hence we defined user defined functions for polar to rect and rect to polar conversions\n",
+      "\n",
+      "\n"
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        "[  8.60232527 -45.        ]\n",
+        "[  10.81665383  135.        ]\n",
+        "[  11.3137085  225.       ]\n",
+        "[  8.48528137  45.        ]\n"
+       ]
+      }
+     ],
+     "prompt_number": 1
+    },
+    {
+     "cell_type": "heading",
+     "level": 2,
+     "metadata": {},
+     "source": [
+      "Example 5.2  Page No : 157"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "import math \n",
+      "from numpy import ones\n",
+      "\n",
+      "def  r2p(x,y):\t\t\t#function to convert recmath.tangular to polar\n",
+      "    polar  =  ones(2)\n",
+      "    polar[0]  =  math.sqrt ((x **2) +(y**2))\n",
+      "    polar[1]  =  math.atan (y/x)\n",
+      "    polar[1]  = (polar [1]*180)/math.pi\n",
+      "    return polar\n",
+      "    \n",
+      "def  p2r(r,theta):\t\t\t#function to convert polar to recmath.tangular\n",
+      "    rect  =  ones(2)\n",
+      "    theta  = ( theta *math.pi) /180\n",
+      "    rect [0] = r* math.cos(theta)\n",
+      "    rect [1] = r* math.sin(theta)\n",
+      "    return rect\n",
+      "    \n",
+      "#CALCULATIONS\n",
+      "#for subdivision 1\n",
+      "I1 = p2r(10,60);\n",
+      "I2 = p2r(8,-45);\n",
+      "I3 = I1+I2;\n",
+      "print (I3);\n",
+      "I4 = r2p(I3[0],I3[1]);\n",
+      "print (I4)\n",
+      "#for subdivision 2\n",
+      "I5 = r2p(5,4);\n",
+      "I6 = r2p(-4,-6);\n",
+      "I7 = ones(2)\n",
+      "I7[0] = (I5[0])*(I6[0]);\n",
+      "I7[1] = (I5[1]+I6[1]);\n",
+      "I7[1] = I7[1]-180;\n",
+      "print (I7);\n",
+      "#for subdivision 3\n",
+      "I8 = r2p(-2,-5);\n",
+      "I9 = r2p(5,7);\n",
+      "I10 = ones(2)\n",
+      "I10[0] = I8[0]/I9[0];\n",
+      "I10[1] = I8[1]-I9[1];\n",
+      "I10[1] = I10[1]-180\n",
+      "print (I10);\n",
+      "#note:here direct functions for converson are not available and hence we defined user defined functions for polar to rect and rect to polar conversions\n"
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        "[ 10.65685425   3.00339979]\n",
+        "[ 11.07198956  15.73932193]\n",
+        "[  46.17358552 -135.        ]\n",
+        "[   0.62601269 -161.56505118]\n"
+       ]
+      }
+     ],
+     "prompt_number": 3
+    },
+    {
+     "cell_type": "heading",
+     "level": 2,
+     "metadata": {},
+     "source": [
+      "Example 5.3  Page No : 160"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "import math \n",
+      "\n",
+      "#given i(t) = 5*math.sin(314*t+(2*math.pi/3))&& v(t) = 20*math.sin(314*t+(5*math.pi/6))\n",
+      "#CALCULATIONS\n",
+      "P1 = 2*(math.pi/3);\t\t\t#phase angle of current in radians\n",
+      "P1 = P1*(180/math.pi);\t\t\t#phase angle of current in degrees\n",
+      "P2 = 5*(math.pi/6);\t\t\t#phase angle of voltage in radians\n",
+      "P2 = P2*(180/math.pi);\t\t\t#phase angle of voltage in degrees\n",
+      "P3 = P2-P1;\t\t\t#current lags voltage by P3 degrees\n",
+      "P4 = P3*math.pi/180;\n",
+      "pf = math.cos(P4);\t\t\t#lagging pf\n",
+      "Vm = 20;\t\t\t#peak voltage\n",
+      "Im = 5;\t\t\t#peak current\n",
+      "Z = Vm/Im;\t\t\t#impedance in ohms\n",
+      "R = (Z)*math.cos(P4);\t\t\t#resistance in ohms\n",
+      "Xl = math.sqrt((Z)**2-(R)**2);\t\t\t#reacmath.tance \n",
+      "W = 314;\n",
+      "L = Xl/W;\t\t\t#inductance in henry\n",
+      "V = Vm/math.sqrt(2);\t\t\t#average value of voltage\n",
+      "I = Im/math.sqrt(2);\t\t\t#average value of current\n",
+      "av = (V*I)*math.cos(P4);\t\t\t#average power in watts\n",
+      "print \"thus impedance, resistance, inductance, powerfactor and average power are %d ohms, %1.2f ohms, %g H,%1.3f and %2.1f W respectively\"%(Z,R,L,pf,av);"
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        "thus impedance, resistance, inductance, powerfactor and average power are 4 ohms, 3.46 ohms, 0.00636943 H,0.866 and 43.3 W respectively\n"
+       ]
+      }
+     ],
+     "prompt_number": 4
+    },
+    {
+     "cell_type": "heading",
+     "level": 2,
+     "metadata": {},
+     "source": [
+      "Example 5.4  Page No : 161"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "import math \n",
+      "#Chapter-5, Example 5.4, Page 161\n",
+      "\n",
+      "#INPUT DATA\n",
+      "I = 10.;\t\t\t#given current in A\n",
+      "P = 1000.;\t\t\t#power in Watts\n",
+      "V = 250.;\t\t\t#voltage in volts\n",
+      "f = 25.;\t\t\t#frequency in Hz\n",
+      "\n",
+      "#CALCULATIONS\n",
+      "R = P/((I)**2);\t\t\t#resistance in ohms\n",
+      "Z = V/I;\t\t\t#impedance in ohms\n",
+      "Xl = math.sqrt((Z)**2-(R)**2);\t\t\t#reacmath.tance in ohms\n",
+      "L = Xl/(2*math.pi*f);\t\t\t#inductance in Henry\n",
+      "Pf = R/Z;\t\t\t#power factor,lagging,pf = math.cos(phi)\n",
+      "\n",
+      "# Results\n",
+      "print \"thus impedance, resistance, inductance, reactance and powerfactor are %d ohms, %d ohms, %1.3f H, \\\n",
+      "%2.2f ohms and %1.1f respectively\"%(Z,R,L,Xl,Pf);\n"
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        "thus impedance, resistance, inductance, reactance and powerfactor are 25 ohms, 10 ohms, 0.146 H, 22.91 ohms and 0.4 respectively\n"
+       ]
+      }
+     ],
+     "prompt_number": 5
+    },
+    {
+     "cell_type": "heading",
+     "level": 2,
+     "metadata": {},
+     "source": [
+      "Example 5.5  Page No : 162"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "import math \n",
+      "#Chapter-5, Example 5.5, Page 162\n",
+      "\n",
+      "#INPUT DATA\n",
+      "V = 250.;\t\t\t#supply voltage in volts\n",
+      "f = 50.;\t\t\t#frequency in hz\n",
+      "Vr = 125.;\t\t\t#voltage across resistance in volts\n",
+      "Vc = 200.;\t\t\t#voltage across coil in volts\n",
+      "I = 5.;\t\t\t#current in A\n",
+      "#CALCULATIONS\n",
+      "R = Vr/I;\t\t\t#resistance in ohms\n",
+      "Z1 = Vc/I;\t\t\t#impedance of coil in ohms\n",
+      "#Z1 = math.sqrt((R1)**2+(Xl)**2)------eqn(1)\n",
+      "Z = V/I;\t\t\t#total impedance in ohms\n",
+      "#Z = math.sqrt((R+R1)**2+(Xl)**2)-----eqn(2)\n",
+      "#solving eqn(1)and eqn(2) we get R1 as follows\n",
+      "R1 = (((Z)**2-(Z1)**2)-(R)**2)/(2*R);\t\t\t#in ohms\n",
+      "Xl = math.sqrt((Z1)**2-(R1)**2);\t\t\t#reacmath.tance of coil in ohms\n",
+      "P = ((I)**2*R1);\t\t\t#power absorbed by the coil in Watts\n",
+      "Pt = ((I)**2)*(R+R1);\t\t\t#total power in Watts\n",
+      "\n",
+      "# Results\n",
+      "print \"thus impedance, resistance, reactance are %d ohms, %d ohms, %2.2f ohms respectively\"%(Z1,R,Xl);\n",
+      "print \"power absorbed and total power are %3.1f W and %3.1f W respectively\"%(P,Pt)"
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        "thus impedance, resistance, reactance are 40 ohms, 25 ohms, 39.62 ohms respectively\n",
+        "power absorbed and total power are 137.5 W and 762.5 W respectively\n"
+       ]
+      }
+     ],
+     "prompt_number": 6
+    },
+    {
+     "cell_type": "heading",
+     "level": 2,
+     "metadata": {},
+     "source": [
+      "Example 5.6  Page No : 163"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "import math \n",
+      "#Chapter-5, Example 5.6, Page 163\n",
+      "\n",
+      "#INPUT DATA\n",
+      "V = 240;\t\t\t#supply voltage in volts\n",
+      "Vl = 171;\t\t\t#voltage across inductor in volts\n",
+      "I = 3;\t\t\t#current in A\n",
+      "phi = 37;\t\t\t#power factor laggging in degrees\n",
+      "#CALCULATIONS\n",
+      "Zl = Vl/I;\t\t\t#impedance of coil in ohms\n",
+      "#Zl = math.sqrt((R1)**2+(Xl)**2)------eqn(1)\n",
+      "Z = V/I;\t\t\t#total impedance in ohms\n",
+      "#Z = math.sqrt((R+R1)**2+(Xl)**2)-----eqn(2)\n",
+      "pf = math.cos(phi*math.pi/180);\t\t\t#powerfactor\n",
+      "Rt = pf*Z;\t\t\t#total resistance in ohms\t\t\t#Rt = (R+R1)\n",
+      "#substituting Rt value in eqn(2) we find Xl as follows\n",
+      "Xl = math.sqrt((Z)**2-(Rt)**2);\t\t\t#reacmath.tance of inductor in ohms\n",
+      "#ubstituting Xl value in eqn(1) we find R1 as follows\n",
+      "R1 = math.sqrt((Zl)**2-(Xl)**2);\t\t\t#resistance of inductor in ohms\n",
+      "R = Rt-R1;\t\t\t#resistance of resistor in ohms\n",
+      "print \"Thus resistance of resistor  is %2.2f ohms\"%(R);\n",
+      "print \"Thus resisimath.tance and reacmath.tance of inductor are %2.2f ohms and %2.2f ohms respectively\"%(R1,Xl)"
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        "Thus resistance of resistor  is 33.38 ohms\n",
+        "Thus resisimath.tance and reacmath.tance of inductor are 30.51 ohms and 48.15 ohms respectively\n"
+       ]
+      }
+     ],
+     "prompt_number": 7
+    },
+    {
+     "cell_type": "heading",
+     "level": 2,
+     "metadata": {},
+     "source": [
+      "Example 5.7  Page No : 164"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "import math \n",
+      "#Chapter-5, Example 5.7, Page 164\n",
+      "\n",
+      "#INPUT DATA\n",
+      "V = 100;\t\t\t#supply voltage in volts\n",
+      "#for COIL A\n",
+      "f = 50;\t\t\t#frequency in Hz\n",
+      "I1 = 8;\t\t\t#current in A\n",
+      "P1 = 120;\t\t\t#power in Watts\n",
+      "#for COIL B\n",
+      "I2 = 10;\t\t\t#current in A\n",
+      "P2 = 500;\t\t\t#power in Watts\n",
+      "#CALCULATIONS\n",
+      "#FOR COIL A\n",
+      "Z1 = V/I1;\t\t\t#impedance of coil A in ohms\n",
+      "R1 = P1/(I1)**2;\t\t\t#resistance of coil A in ohms\n",
+      "X1 = math.sqrt(((Z1)**2-(R1)**2));\t\t\t#reacmath.tance of coil A in ohms\n",
+      "#FOR COIL B\n",
+      "Z2 = V/I2;\t\t\t#impedance of coil B in ohms\n",
+      "R2 = P2/(I2)**2;\t\t\t#resistance of coil B in ohms\n",
+      "X2 = math.sqrt(((Z2)**2-(R2)**2));\t\t\t#reacmath.tance of coil B in ohms\n",
+      "#When both COILS A and B are in series\n",
+      "Rt = R1+R2;\t\t\t#total resistance in ohms\n",
+      "Xt = X1+X2;\t\t\t#total reacmath.tance in ohms\n",
+      "Zt = math.sqrt((Rt)**2+(Xt)**2);\t\t\t#total impedance in ohms\n",
+      "It = V/Zt;\t\t\t#current drawn in A\n",
+      "P = ((It)**2)*(Rt);\t\t\t#power taken in watts\n",
+      "print \"Thus current drawn and power taken in watts are %2.2f A and %3.2f W respectively\"%(It,P);"
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        "Thus current drawn and power taken in watts are 4.66 A and 130.12 W respectively\n"
+       ]
+      }
+     ],
+     "prompt_number": 8
+    },
+    {
+     "cell_type": "heading",
+     "level": 2,
+     "metadata": {},
+     "source": [
+      "Example 5.8  Page No : 167"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "import math \n",
+      "#Chapter-5, Example 5.8, Page 167\n",
+      "\n",
+      "#INPUT DATA\n",
+      "R = 100;\t\t\t#resistance in ohms\n",
+      "C = 50*10**-6;\t\t\t#capacitance in F\n",
+      "V = 200;\t\t\t#voltage in Volts\n",
+      "f = 50;\t\t\t#frequency in Hz\n",
+      "#Z = R-(1j)*(Xc)------>impedance\n",
+      "Xc = 1/(2*math.pi*f*C);\t\t\t#capacitive reacmath.tance in ohms\n",
+      "Z = math.sqrt((R)**2+(Xc)**2);\t\t\t#impedance in ohms\n",
+      "I = V/Z;\t\t\t#current in A\n",
+      "pf = R/Z;\t\t\t#power factor ------>math.cos(phi)---->leading\n",
+      "phi = math.acos(0.844);\t\t\t#phase angle in radians\n",
+      "phi = phi*180/math.pi;\t\t\t#phase angle in degrees\n",
+      "Vr = (I)*(R);\t\t\t#voltage across resistor\n",
+      "Vc = (I)*(Xc);\t\t\t#votage across capacitor\n",
+      "print \"Thus impedance, current, powerfactor and phaseangle are %3.2f ohms, %1.2f A, %1.3f and %2.2f degrees respectively\"%(Z,I,pf,phi);\n",
+      "print \"voltage across resistor and capacitor are %d V and %3.2f V respectively\"%(Vr,Vc)\n"
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        "Thus impedance, current, powerfactor and phaseangle are 118.54 ohms, 1.69 A, 0.844 and 32.44 degrees respectively\n",
+        "voltage across resistor and capacitor are 168 V and 107.41 V respectively\n"
+       ]
+      }
+     ],
+     "prompt_number": 9
+    },
+    {
+     "cell_type": "heading",
+     "level": 2,
+     "metadata": {},
+     "source": [
+      "Example 5.9  Page No : 169"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "import math \n",
+      "#Chapter-5, Example 5.9, Page 169\n",
+      "\n",
+      "#INPUT DATA\n",
+      "phi = 40;\t\t\t#phase in degrees\n",
+      "V = 150;\t\t\t#voltage in Volts\n",
+      "I = 8;\t\t\t#current in A\n",
+      "#the applied voltage lags behind the current .That means the current leads the voltage\n",
+      "#hence pf is leading\n",
+      "#CALCULATIONS\n",
+      "pf = math.cos(phi*math.pi/180);\t\t\t#in degrees--->leading\n",
+      "#hence it is a capacitive circuit\n",
+      "pa = V*I*pf;\t\t\t#active power in W\n",
+      "pr = V*I*math.sin(phi*math.pi/180);\t\t\t#reactive power in VAR\n",
+      "print \"Thus active and reactive power are %3.1f W and %3.1f VAR respectively\"%(pa,pr);"
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        "Thus active and reactive power are 919.3 W and 771.3 VAR respectively\n"
+       ]
+      }
+     ],
+     "prompt_number": 10
+    },
+    {
+     "cell_type": "heading",
+     "level": 2,
+     "metadata": {},
+     "source": [
+      "Example 5.10  Page No : 169"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "import math \n",
+      "#Chapter-5, Example 5.10, Page 169\n",
+      "\n",
+      "#INPUT DATA\n",
+      "#given v = 141.4*math.sin(314*t)\n",
+      "P = 700.;\t\t\t#power in Watts\n",
+      "pf = 0.707;\t\t\t#powerfactor------>leading------>math.cos(phi)\n",
+      "Vm = 141.4;\t\t\t#maximum value of supply voltage\n",
+      "#CALCULATIONS\n",
+      "Vr = Vm/(math.sqrt(2));\t\t\t#rms value of supply voltage\n",
+      "I = P/(Vr*pf);\t\t\t#current in A\n",
+      "Z = Vr/I;\t\t\t#impedance in ohms\n",
+      "R = (Z)*(pf);\t\t\t#resistance in ohms\n",
+      "phi = math.acos(pf)*180/math.pi;\t\t\t#angle in degrees\n",
+      "Xc = (Z)*(math.sin(phi));\t\t\t#reacmath.tance in ohms\n",
+      "C = 1/(3.14*7.13);\t\t\t#capacitance in F\n",
+      "print \"Thus resistance and capacitance are %1.2f ohms and %g F respectively\"%(R,C);"
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        "Thus resistance and capacitance are 7.14 ohms and 0.0446664 F respectively\n"
+       ]
+      }
+     ],
+     "prompt_number": 13
+    },
+    {
+     "cell_type": "heading",
+     "level": 2,
+     "metadata": {},
+     "source": [
+      "Example 5.11  Page No : 169"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "import math \n",
+      "#Chapter-5, Example 5.11, Page 169\n",
+      "\n",
+      "#INPUT DATA\n",
+      "V = 200.;\t\t\t#supply voltage in volts\n",
+      "f = 50.;\t\t\t#freq in hz\n",
+      "P = 7000.;\t\t\t#power in Watts\n",
+      "Vr = 130.;\t\t\t#volatge across resistor in volts\n",
+      "P = 7000.;\t\t\t#power in Watts\n",
+      "\n",
+      "#CALCULATIONS\n",
+      "R = ((Vr)**2)/P;\t\t\t#resistance in ohms\n",
+      "I = Vr/R;\t\t\t#current in A\n",
+      "Z = V/I;\t\t\t#total impedance in ohms\n",
+      "Xc = math.sqrt((Z)**2-(R)**2);\n",
+      "C = 1/(2*math.pi*f*Xc);\t\t\t#capacitance in F\n",
+      "pf = R/Z;\t\t\t#power factor------>leading\n",
+      "phi = math.acos(pf);\t\t\t#angle in radians\n",
+      "phi = phi*180/math.pi;\t\t\t#angle in degrees\n",
+      "Vm = V*math.sqrt(2);\t\t\t#maximum value of voltage\n",
+      "#voltage equation v = Vm*math.sin(2*math.pi*f*t)------>282.84*math.sin(314.16*t)\n",
+      "#current leads voltage by phi\n",
+      "#current equation ------>i = 76.155*math.sin(314.16*t+phi)\n",
+      "print \"Thus current, resistance, p.f, capacitance, impedance are %2.2f A , %1.2f ohms, %2.1f , \\\n",
+      "%g F and %1.2f ohms respectively\"%(I,R,pf,C,Z);"
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        "Thus current, resistance, p.f, capacitance, impedance are 53.85 A , 2.41 ohms, 0.6 , 0.00112771 F and 3.71 ohms respectively\n"
+       ]
+      }
+     ],
+     "prompt_number": 14
+    },
+    {
+     "cell_type": "heading",
+     "level": 2,
+     "metadata": {},
+     "source": [
+      "Example 5.12  Page No : 170"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "import math \n",
+      "\n",
+      "#INPUT DATA\n",
+      "C = 50.;\t\t\t#capacitance in uf\n",
+      "R = 100.;\t\t\t#resistance in ohms\n",
+      "V = 200.;\t\t\t#supply voltage in volts\n",
+      "f = 50.;\t\t\t#freq in hz\n",
+      "#CALCULATIONS\n",
+      "Xc = 1/(2*math.pi*f*C*10**-6);\t\t\t#capacitive reacmath.tance in ohms\n",
+      "Z = R-((1j)*Xc);\t\t\t#impedance in ohms\n",
+      "print (Z);\n",
+      "z1 = math.sqrt((R)**2+(Xc)**2);\n",
+      "theta = math.atan(Xc/R);\n",
+      "pf = math.cos(theta);\t\t\t#powerfactor\n",
+      "I = V/z1;\t\t\t#current in A\n",
+      "P = V*I*pf;\t\t\t#power in Watts\n",
+      "print \"Thus current, power factor, power are % 1.2f A ,%1.3f ,%d W respectively\"%(I,pf,P);"
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        "(100-63.6619772368j)\n",
+        "Thus current, power factor, power are  1.69 A ,0.844 ,284 W respectively\n"
+       ]
+      }
+     ],
+     "prompt_number": 15
+    },
+    {
+     "cell_type": "heading",
+     "level": 2,
+     "metadata": {},
+     "source": [
+      "Example 5.13  Page No : 170"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "import math \n",
+      "\n",
+      "#INPUT DATA\n",
+      "C = 0.05;\t\t\t#capacitance in uf\n",
+      "F = 500;\t\t\t#freq in hz\n",
+      "#CALCULATIONS\n",
+      "Xl = 1/(2*math.pi*F*C*10**-6);\t\t\t#capacitive reacmath.tance in ohms\n",
+      "#at resonance Xl = Xc \n",
+      "L = (Xl/(2*math.pi*F));\t\t\t#inductance in H\n",
+      "print \"Thus value of L is %1.2f H\"%(L);"
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        "Thus value of L is 2.03 H\n"
+       ]
+      }
+     ],
+     "prompt_number": 16
+    },
+    {
+     "cell_type": "heading",
+     "level": 2,
+     "metadata": {},
+     "source": [
+      "Example 5.14  Page No : 171"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "import math \n",
+      "\n",
+      "#INPUT DATA\n",
+      "V = 200;\t\t\t#voltage in V\n",
+      "R = 50;\t\t\t#resistance in ohms\n",
+      "L = 0.5;\t\t\t#inductance in Henry\n",
+      "F = 50;\t\t\t#freq in hz\n",
+      "#CALCULATIONS\n",
+      "Xl = 2*math.pi*F*L;\t\t\t#inductive reacmath.tance\n",
+      "Z = (R)+((1j)*Xl)\t\t\t#impedance\n",
+      "print (Z);\n",
+      "z1 = math.sqrt((R)**2+(Xl)**2);\t\t\t#magnitude\n",
+      "theta = math.atan(Xl/R);\t\t\t#angle in radians\n",
+      "I = V/z1;\t\t\t#current in A\n",
+      "P = V*I*math.cos(theta);\t\t\t#power supplied in W\n",
+      "#here capacitive reacmath.tance equals inductive reacmath.tance\n",
+      "#hence Xc = Xl\n",
+      "C = 1/(2*math.pi*F*Xl);\t\t\t#capacitance in uf\n",
+      "r = (V/I)-(R);\t\t\t#additional resistance to be added in series\n",
+      "print \"Thus current and power required are % 1.2f A and %2.2f W respectively\"%(I,P);\n",
+      "print \"Thus additional resistance that neede to be connected in series with R and C to have\\\n",
+      " same current at unity power factor is %1.1f ohms\"%(r);\n"
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        "(50+157.079632679j)\n",
+        "Thus current and power required are  1.21 A and 73.60 W respectively\n",
+        "Thus additional resistance that neede to be connected in series with R and C to have same current at unity power factor is 114.8 ohms\n"
+       ]
+      }
+     ],
+     "prompt_number": 18
+    },
+    {
+     "cell_type": "heading",
+     "level": 2,
+     "metadata": {},
+     "source": [
+      "Example 5.15  Page No : 171"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "import math \n",
+      "\n",
+      "#INPUT DATA\n",
+      "R = 50.;\t\t\t#resistance in ohms\n",
+      "L = 9.;\t\t\t#inductance in Henry\n",
+      "I0 = 1.;\t\t\t#current in A\n",
+      "f = 75.;\t\t\t#ferquency in Hz\n",
+      "#at resonance Xl = Xc \n",
+      "#CALCULATIONS\n",
+      "Xl = 2*math.pi*f*L;\t\t\t#inductive reacmath.tance\n",
+      "Xc = Xl;\t\t\t#capacitive reacmath.tance\n",
+      "C = 1/(2*math.pi*f*Xc);\t\t\t#capacitance in uf\n",
+      "print \"Thus capacitance is %g F\"%(C);"
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        "Thus capacitance is 5.00352e-07 F\n"
+       ]
+      }
+     ],
+     "prompt_number": 19
+    },
+    {
+     "cell_type": "heading",
+     "level": 2,
+     "metadata": {},
+     "source": [
+      "Example 5.16  Page No : 175"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "import math \n",
+      "\n",
+      "#INPUT DATA\n",
+      "R = 10.;\t\t\t#resistance in ohms\n",
+      "L = 0.1;\t\t\t#inductance in Henry\n",
+      "C = 150.;\t\t\t#capacitor in uf\n",
+      "V = 200.;\t\t\t#voltage in V\n",
+      "f = 50.;\t\t\t#frequency in hz\n",
+      "#CALCULATIONS\n",
+      "Xc = 1/(2*math.pi*f*C*10**-6);\t\t\t#Capacitive reacmath.tance in ohms\n",
+      "Xl = (2*math.pi*f*L);\t\t\t#inductive reacmath.tance in ohms\n",
+      "Z = R+((1j)*(Xl-Xc));\t\t\t#impedance in ohms\n",
+      "z1 = math.sqrt((R)**2+(Xl-Xc)**2);\t\t\t#magnitude of Z\n",
+      "I = V/z1;\t\t\t#current in A\n",
+      "pf = R/z1;\t\t\t#power factor----->math.cos(phi)\n",
+      "#As Xl-Xc is inductive,pf is lagging\n",
+      "z2 = math.sqrt((R**2)+(Xl)**2);\t\t\t#impedance of coil in ohms\n",
+      "Vl = I*(z2);\t\t\t#voltage across coil in volts\n",
+      "Vc = I*(Xc);\t\t\t#voltage across capacitor in volts\n",
+      "print \"Thus inductive reacmath.tance, capacitive reacmath.tance, impedance, current, powerfactor are %2.2f ohms, \\\n",
+      "%2.2f ohms, %2.2f ohms, %d A, %1.1f respectively\"%(Xl,Xc,z1,I,pf);\n"
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        "Thus inductive reacmath.tance, capacitive reacmath.tance, impedance, current, powerfactor are 31.42 ohms, 21.22 ohms, 14.28 ohms, 14 A, 0.7 respectively\n"
+       ]
+      }
+     ],
+     "prompt_number": 20
+    },
+    {
+     "cell_type": "heading",
+     "level": 2,
+     "metadata": {},
+     "source": [
+      "Example 5.17  Page No : 176"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "import math \n",
+      "\n",
+      "#INPUT DATA\n",
+      "L = 10;\t\t\t#inductance in milliHenry\n",
+      "C = 5;\t\t\t#capacitor in uf\n",
+      "phi = 50;\t\t\t#phase in degrees-------->lagging\n",
+      "f = 500;\t\t\t#frequency in hz\n",
+      "V = 200;\t\t\t#supply voltage in volts\n",
+      "\n",
+      "#CALCULATIONS\n",
+      "Xc = 1/(2*math.pi*f*C*10**-6);\t\t\t#Capacitive reacmath.tance in ohms\n",
+      "Xl = (2*math.pi*f*L*10**-3);\t\t\t#inductive reacmath.tance in ohms\n",
+      "R = (Xc-Xl)/(math.tan(phi*math.pi/180));\t\t\t#resistance in ohms\n",
+      "Z = math.sqrt((R)**2+(Xc-Xl)**2);\t\t\t#impedance in ohms\n",
+      "I = V/Z;\t\t\t#current in A\n",
+      "Vr = (I)*(R);\t\t\t#voltage across resistance\n",
+      "Vl = (I)*(Xl);\t\t\t#voltage across inductance\n",
+      "Vc = (I)*(Xc);\t\t\t#voltage across capacitance\n",
+      "print \"Thus voltages across resistance, inductance, capacitance are %3.2f volts, %3.2f volts, %3.2f volts respectively\"%(Vr,Vl,Vc);"
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        "Thus voltages across resistance, inductance, capacitance are 128.56 volts, 149.26 volts, 302.47 volts respectively\n"
+       ]
+      }
+     ],
+     "prompt_number": 21
+    },
+    {
+     "cell_type": "heading",
+     "level": 2,
+     "metadata": {},
+     "source": [
+      "Example 5.18  Page No : 176"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "import math \n",
+      "from sympy import Symbol,solve\n",
+      "#Chapter-5, Example 5.18, Page 176\n",
+      "\n",
+      "#INPUT DATA\n",
+      "L = 5;\t\t\t#inductance in Henry\n",
+      "f = 50;\t\t\t#frequency in hz\n",
+      "V = 230;\t\t\t#supply voltage in volts\n",
+      "R = 2;\t\t\t#resistance in ohms\n",
+      "V1 = 250;\t\t\t#voltage across coil in V\n",
+      "\n",
+      "#CALCULATIONS\n",
+      "Xl = (2*math.pi*f*L);\t\t\t#inductive reacmath.tance in ohms\n",
+      "Z1 = math.sqrt((R)**2+(Xl)**2);\t\t\t#impedance of coil in ohms\n",
+      "I = V1/Z1;\t\t\t#current in A\n",
+      "Z = V/I;\t\t\t#total impedance in ohms\n",
+      "#Z = math.sqrt((R)**2+(Xl-Xc)**2) and solving for Xc\n",
+      "Xc = Symbol(\"Xc\");\n",
+      "p = (Xc**2)-3141.58*(Xc)+378004\n",
+      "roots2  =  solve(p);\n",
+      "r2  =  roots2[1];\n",
+      "#Xc cannot be greater than Z\n",
+      "C = 1/(2*math.pi*f*r2);\t\t\t#capacitance in F\n",
+      "print \"Thus value of C that must be present suct that voltage across coil is 250 volts is %g F respectively\"%(C);"
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        "Thus value of C that must be present suct that voltage across coil is 250 volts is 1.05531e-06 F respectively\n"
+       ]
+      }
+     ],
+     "prompt_number": 22
+    },
+    {
+     "cell_type": "heading",
+     "level": 2,
+     "metadata": {},
+     "source": [
+      "Example 5.19  Page No : 178"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "import math \n",
+      "#Chapter-5, Example 5.19, Page 178\n",
+      "\n",
+      "#v = 350*math.cos(3000*t-20)\n",
+      "#i = 15*math.cos(3000*t-60)\n",
+      "#INPUT DATA\n",
+      "L = 0.5;\t\t\t#inductance in Henry\n",
+      "phi = -40;\t\t\t#phase difference between applied voltage and current\n",
+      "#Xl>Xc(P.f is lagging)\n",
+      "w = 3000;\t\t\t#freq in hz\n",
+      "Vm = 350;\t\t\t#peak voltage in volts\n",
+      "Im = 15;\t\t\t#peak current in amps\n",
+      "#CALCULATIONS\n",
+      "Z = Vm/Im;\t\t\t#total impedance in ohms\n",
+      "#Xl-Xc = 0.839*R = X\n",
+      "#Z = math.sqrt((R)**2+(X)**2)\n",
+      "#Z = 1.305*R\n",
+      "R = Z/1.305;\t\t\t#resistance in ohms\n",
+      "X = 0.839*R;\t\t\t#\n",
+      "#X = Xl-Xc\n",
+      "Xl = w*L;\t\t\t#reactive inductance in ohms\n",
+      "Xc = Xl-X;\t\t\t#capacitive reacmath.tance in ohms\n",
+      "C = 1/(w*Xc);\t\t\t#capacitance in uf\n",
+      "print \"Thus resistance and capacitance are %2.2f ohms and %g F respectively\"%(R,C);"
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        "Thus resistance and capacitance are 17.62 ohms and 2.24435e-07 F respectively\n"
+       ]
+      }
+     ],
+     "prompt_number": 23
+    },
+    {
+     "cell_type": "heading",
+     "level": 2,
+     "metadata": {},
+     "source": [
+      "Example 5.20  Page No : 182"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "import math \n",
+      "from numpy.linalg import inv\n",
+      "from scipy.optimize import fsolve\n",
+      "from sympy.solvers import solve\n",
+      "\n",
+      "\n",
+      "#INPUT DATA\n",
+      "R = 10;\t\t\t#resistance in ohms\n",
+      "L = 0.1;\t\t\t#inductance in henry\n",
+      "f = 50;\t\t\t#frequency in hz\n",
+      "#CALCULATIONS\n",
+      "Xl = (2*math.pi*f*L);\t\t\t#inductive reacmath.tance in ohms\n",
+      "Z = R+((1j)*(Xl));\t\t\t#impedance in ohms\n",
+      "Y = inv([[Z]])#[0];\t\t\t#admittance in mho\n",
+      "y = abs(Y);\t\t\t#admittance in mho\n",
+      "print \"admittance is %1.5f mho\"%(y);\n",
+      "\n"
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        "admittance is 0.03033 mho\n"
+       ]
+      }
+     ],
+     "prompt_number": 12
+    },
+    {
+     "cell_type": "heading",
+     "level": 2,
+     "metadata": {},
+     "source": [
+      "Example 5.21  Page No : 182"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "import math \n",
+      "from numpy.linalg import inv\n",
+      "\n",
+      "#INPUT DATA\n",
+      "#CALCULATIONS\n",
+      "Z = 10+((1j)*(5));\t\t\t#impedance in ohms\n",
+      "Y = inv([[Z]]);\t\t\t#admittance in mho\n",
+      "print (Y);\n",
+      "\n"
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        "[[ 0.08-0.04j]]\n"
+       ]
+      }
+     ],
+     "prompt_number": 15
+    },
+    {
+     "cell_type": "heading",
+     "level": 2,
+     "metadata": {},
+     "source": [
+      "Example 5.22  Page No : 182"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "import math \n",
+      "\n",
+      "#INPUT DATA\n",
+      "Z1 = 7.+((1j)*5);\t\t\t#impedance of branch1 in ohms\n",
+      "Z2 = 10.-((1j)*8);\t\t\t#impedance  of branch2 in ohms\n",
+      "V = 230.;\t\t\t#supply voltage in volts\n",
+      "f = 50.;\t\t\t#frequency in hz\n",
+      "#CALCULATIONS\n",
+      "Y1 = 1/(Z1);\t\t\t#admittance of branch1 in mho\n",
+      "Y2 = 1/(Z2);\t\t\t#admittance of branch2 in mho\n",
+      "Y = Y1+Y2;\t\t\t#admittance of combined circuit\n",
+      "print (Y);\n",
+      "g = abs(Y);\t\t\t#conductance in mho;\n",
+      "B = math.atan(Y.imag/Y.real);\t\t\t#susceptance in mho\n",
+      "I = V*(Y);\t\t\t#current\n",
+      "print (I);\t\t\t#total current taken from mains in A\n",
+      "z = math.atan(I.imag/I.real);\n",
+      "pf = math.cos(z);\t\t\t#power factor\n",
+      "print \"thus conductance and susceptance of the circuit is %1.3f mho and %1.3f mho respectively\"%(g,B);\n",
+      "print \"power factor is %1.3f lagging\"%(pf)\n",
+      "\n"
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        "(0.155570204351-0.0187870797627j)\n",
+        "(35.7811470007-4.32102834542j)\n",
+        "thus conductance and susceptance of the circuit is 0.157 mho and -0.120 mho respectively\n",
+        "power factor is 0.993 lagging\n"
+       ]
+      }
+     ],
+     "prompt_number": 14
+    },
+    {
+     "cell_type": "heading",
+     "level": 2,
+     "metadata": {},
+     "source": [
+      "Example 5.23  Page No : 183"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "import math \n",
+      "#Chapter-5, Example 5.23, Page 183\n",
+      "\n",
+      "#INPUT DATA\n",
+      "V = 240.;\t\t\t#voltage in volts\n",
+      "f = 50.;\t\t\t#frequency in Hz\n",
+      "R = 15.;\t\t\t#resisimath.tance in ohms\n",
+      "I = 22.1;\t\t\t#current in A\n",
+      "#CALCULATIONS\n",
+      "G = 1/R;\t\t\t#conductance in mho\n",
+      "#susceptance of the circuit,B = 1/(Xl) = 0.00318/L\n",
+      "#admittance of the circuit,(G-jB) = (0.067-j(0.00318/L))\n",
+      "Y = I/V;\t\t\t#admittance in mho;\n",
+      "#Y = math.sqrt((0.067)**2+(0.00318/L)**2) = 0.092-----eqn(1)\n",
+      "#solving eqn(1) for L we have it as\n",
+      "L = math.sqrt((0.00318)**2/((Y)**2-(G)**2));\t\t\t#inductance in henry\n",
+      "#when current is 34A\n",
+      "I1 = 34;\t\t\t#current in A\n",
+      "Y1 = I1/V;\t\t\t#admittance in mho\n",
+      "#for Y1 we need to find f \n",
+      "f1 = math.sqrt((3.183)**2/((Y1)**2-(G)**2));\t\t\t#frequency in hz\n",
+      "print \"Thus value of frequency when current is 34A is %2.1f Hz\"%(f1);\n",
+      "\n"
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        "Thus value of frequency when current is 34A is 25.5 Hz\n"
+       ]
+      }
+     ],
+     "prompt_number": 29
+    },
+    {
+     "cell_type": "heading",
+     "level": 2,
+     "metadata": {},
+     "source": [
+      "Example 5.24  Page No : 184"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "import math \n",
+      "from numpy.linalg import inv\n",
+      "\n",
+      "#Chapter-5, Example 5.24, Page 184\n",
+      "\n",
+      "#INPUT DATA\n",
+      "L = 0.05;\t\t\t#inductance in henry\n",
+      "R2 = 20.;\t\t\t#resistance in ohms\n",
+      "R1 = 15.;\t\t\t#resistance in ohms\n",
+      "V = 200.;\t\t\t#supply voltage in volts\n",
+      "f = 50.;\t\t\t#frequency in hz\n",
+      "#CALCULATIONS\n",
+      "#for branch 1\n",
+      "Z1 = (R1)+((1j)*(2*math.pi*f*L));\t\t\t#impedance in ohms\n",
+      "Y1 = inv([[Z1]]);\t\t\t#admittance in branch\n",
+      "I1 = V*(Y1);\t\t\t#current in branch\n",
+      "print (I1);\n",
+      "i1 = abs(I1);\t\t\t#magnitude of current\n",
+      "#for branch 2\n",
+      "Y2 = 1/R2;\t\t\t#admittance in branch\n",
+      "I2 = V*Y2;\t\t\t#current in branch\n",
+      "i2 = abs(I2);\t\t\t#magnitude of current\n",
+      "I = I1+I2;\t\t\t#total current in A\n",
+      "i = abs(I);\t\t\t#magnitude of total current\n",
+      "theta = math.atan(I.imag/I.real);\t\t\t#angle in radians\n",
+      "theta = theta*(180)/(math.pi);\t\t\t#angle in degrees\n",
+      "print \"Thus current in branch1,branch2 abd total currents are %1.2f A, %d A, %2.2f A respectively\"%(i1,i2,i);\n",
+      "print \"phase angle of the combination is %2.1f degrees\"%(theta);\n",
+      "\n"
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        "[[ 6.3594338-6.6595835j]]\n",
+        "Thus current in branch1,branch2 abd total currents are 9.21 A, 10 A, 17.66 A respectively\n",
+        "phase angle of the combination is -22.2 degrees\n"
+       ]
+      }
+     ],
+     "prompt_number": 18
+    },
+    {
+     "cell_type": "heading",
+     "level": 2,
+     "metadata": {},
+     "source": [
+      "Example 5.25  Page No : 185"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "import math \n",
+      "from numpy.linalg import inv\n",
+      "\n",
+      "#INPUT DATA\n",
+      "L = 6.;\t\t\t#inductance in millihenry\n",
+      "R2 = 50.;\t\t\t#resistance in ohms\n",
+      "R1 = 40.;\t\t\t#resistance in ohms\n",
+      "C = 4.;\t\t\t#capacitance in uf\n",
+      "V = 100.;\t\t\t#voltage in volts\n",
+      "f = 800.;\t\t\t#frequency in hz\n",
+      "#CALCULATIONS\n",
+      "Xl = (2*math.pi*f*L*10**-3);\t\t\t#inductive reacmath.tance in ohms\n",
+      "Xc = 1/(2*math.pi*f*C*10**-6);\t\t\t#capacitive reacmath.tance in ohms\n",
+      "Y1 = inv([[(R1)+(1j*Xl)]]);\t\t\t#admittance of branch1 in mho\n",
+      "Y2 = inv([[(R2)-(1j*Xc)]]);\t\t\t#admittance of branch2 in mho\n",
+      "I1 = V*(Y1);\t\t\t#current in branch 1\n",
+      "I2 = V*(Y2);\t\t\t#current in branch 2\n",
+      "I = I1+I2;\t\t\t#total curremt in A\n",
+      "theta = (math.atan(I1.imag/I1.real))-math.atan(I2.imag/I2.real);\n",
+      "theta = theta*180/math.pi;\t\t\t#angle in degrees\n",
+      "print \"Thus total current taken from supply is %2.2f\"%(abs(I));\n",
+      "print \"phase angle between currents of coil and capacitor is %2.2f degrees\"%(theta);\n",
+      "\n"
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        "Thus total current taken from supply is 2.61\n",
+        "phase angle between currents of coil and capacitor is -81.86 degrees\n"
+       ]
+      }
+     ],
+     "prompt_number": 21
+    },
+    {
+     "cell_type": "heading",
+     "level": 2,
+     "metadata": {},
+     "source": [
+      "Example 5.26  Page No : 186"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "import math \n",
+      "\n",
+      "#INPUT DATA\n",
+      "Z1 = 10+(1j*15);\t\t\t#impedance in ohms\n",
+      "Z2 = 6-(1j*8);\t\t\t#impedance in ohms\n",
+      "I = 15.;\t\t\t#current in A\n",
+      "#CALCULATIONS\n",
+      "I1 = ((Z2)/(Z1+Z2))*(I);\t\t\t#umath.sing current division rule\n",
+      "I2 = ((Z1)/(Z1+Z2))*(I);\t\t\t#umath.sing current division rule\n",
+      "i1 = abs(I1);\t\t\t#magnitude of current 1\n",
+      "i2 = abs(I2);\t\t\t#magnitdude of current 2\n",
+      "P1 = ((i1)**2)*(Z1*(1));\t\t\t#power consumed by branch 1\n",
+      "P2 = ((i2)**2)*(Z2*(1));\t\t\t#power consumed by branch 2\n",
+      "print \"Thus power consumed by branches 1 and 2 are %3.2f W and %4.1f W respectively\"%(P1.real,P2.real);\n",
+      "\n"
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        "Thus power consumed by branches 1 and 2 are 737.70 W and 1438.5 W respectively\n"
+       ]
+      }
+     ],
+     "prompt_number": 24
+    },
+    {
+     "cell_type": "heading",
+     "level": 2,
+     "metadata": {},
+     "source": [
+      "Example 5.27  Page No : 187"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "import math \n",
+      "from numpy.linalg import inv\n",
+      "#Chapter-5, Example 5.27, Page 187\n",
+      "\n",
+      "#INPUT DATA\n",
+      "V = 200.;\t\t\t#voltage in volts\n",
+      "f = 50.;\t\t\t#frequency in hz\n",
+      "R1 = 10.;\t\t\t#resistance in ohms\n",
+      "L1 = 0.0023;\t\t\t#inductance in henry\n",
+      "R2 = 5.;\t\t\t#resistance in ohms\n",
+      "L2 = 0.035;\t\t\t#inductance in henry\n",
+      "#CALCULATIONS\n",
+      "Xl1 = (2*math.pi*f*L1);\t\t\t#inductive reacmath.tance in branch 1 in ohm\n",
+      "Xl2 = (2*math.pi*f*L2);\t\t\t#inductive reacmath.tance in branch 2 in ohm\n",
+      "Y1 = inv([[10+(1j*7.23)]]);\t\t\t#admittance of branch 1 in mho\n",
+      "Y2 = inv([[5+(1j*10.99)]]);\t\t\t#admittance of branch 2 in mho\n",
+      "Y = Y1+Y2;\t\t\t#total admittance in mho\n",
+      "I1 = V*(Y1);\t\t\t#current through branch1\n",
+      "I2 = V*(Y2);\t\t\t#current through branch2\n",
+      "I = I1+I2;\t\t\t#total current in A\n",
+      "theta = math.atan(I.imag/I.real);\t\t\t#angle in radians\n",
+      "pf_of_combination = math.cos(theta);\t\t\t#powerfactor---->lagging\n",
+      "print \"Thus currents in branch1, branch2 and total current are %2.1f A, %2.1f A and %2.2f A respectively\"%(abs(I1),abs(I2),abs(I));\n",
+      "print \"pf of combination is %1.3f\"%(pf_of_combination);"
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        "Thus currents in branch1, branch2 and total current are 16.2 A, 16.6 A and 31.68 A respectively\n",
+        "pf of combination is 0.631\n"
+       ]
+      }
+     ],
+     "prompt_number": 26
+    },
+    {
+     "cell_type": "heading",
+     "level": 2,
+     "metadata": {},
+     "source": [
+      "Example 5.28  Page No : 189"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "import math \n",
+      "\n",
+      "#INPUT DATA\n",
+      "f = 50.;\t\t\t#freq in hz\n",
+      "V = 100.;\t\t\t#volatge in V\n",
+      "L1 = 0.015;\t\t\t#inductance in branch 1 in henry\n",
+      "L2 = 0.08;\t\t\t#inductance in branch 2 in henry\n",
+      "R1 = 2.;\t\t\t#resistance of branch 1 in ohms\n",
+      "x1 = 4.71;\t\t\t#reacmath.tance of branch 1 in ohms\n",
+      "R2 = 1.;\t\t\t#resistance of branch 2 in ohms\n",
+      "x2 = 25.13;\t\t\t#reacmath.tance of branch 2 in ohms\n",
+      "Z1 = (R1)+(1j*x1);\t\t\t#impedance of branch1 in ohms\n",
+      "Z2 = (R2)+(1j*x2);\t\t\t#impedance of branch1 in ohms\n",
+      "I1 = V/Z1;\t\t\t#current in branch 1 in A\n",
+      "print \"current in branch 1 in A\"\n",
+      "print (I1);\n",
+      "I2 = V/Z2;\t\t\t#current in branch 2 in A\n",
+      "print \"current in branch 2 in A\"\n",
+      "print (I2);\n",
+      "I3 = I1+I2;\t\t\t#total current in A\n",
+      "print \"total current in A\"\n",
+      "print (I3);\n",
+      "#note:Answer for real part of total current given in textbook is wrong.Please check the calculations\n",
+      "\n"
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        "current in branch 1 in A\n",
+        "(7.63822319652-17.9880156278j)\n",
+        "current in branch 2 in A\n",
+        "(0.158098542505-3.97301637316j)\n",
+        "total current in A\n",
+        "(7.79632173903-21.961032001j)\n"
+       ]
+      }
+     ],
+     "prompt_number": 27
+    },
+    {
+     "cell_type": "heading",
+     "level": 2,
+     "metadata": {},
+     "source": [
+      "Example 5.29  Page No : 189"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "import math \n",
+      "from numpy.linalg import inv\n",
+      "\n",
+      "#CALCULATIONS\n",
+      "R = 8;\t\t\t#resistance in ohms\n",
+      "Xc = -(1j)*12;\t\t\t#capacitive reacmath.tance in ohms\n",
+      "Y = (inv([[R]])+inv([[Xc]]));\t\t\t#admittance in mho\n",
+      "print (Y);\n",
+      "\n"
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        "[[ 0.125+0.08333333j]]\n"
+       ]
+      }
+     ],
+     "prompt_number": 28
+    },
+    {
+     "cell_type": "heading",
+     "level": 2,
+     "metadata": {},
+     "source": [
+      "Example 5.30  Page No : 189"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "import math \n",
+      "from numpy.linalg import inv\n",
+      "\n",
+      "#CALCULATIONS\n",
+      "R = 3;\t\t\t#resistance in ohms\n",
+      "Xl = (1j)*4;\t\t\t#inductive reacmath.tance in ohms\n",
+      "Y = (inv([[R]])+inv([[Xl]]));\t\t\t#admittance in mho\n",
+      "print (Y);\n",
+      "\n"
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        "[[ 0.33333333-0.25j]]\n"
+       ]
+      }
+     ],
+     "prompt_number": 29
+    },
+    {
+     "cell_type": "heading",
+     "level": 2,
+     "metadata": {},
+     "source": [
+      "Example 5.31  Page No : 196"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "import math \n",
+      "\n",
+      "#INPUT DATA\n",
+      "R = 10.;\t\t\t#resistance in ohms\n",
+      "L = 10.;\t\t\t#inductance in milli henry\n",
+      "C = 1.;\t\t\t#capacitance in uF\n",
+      "V = 200.;\t\t\t#applied voltage in volts\n",
+      "\n",
+      "#CALCULATIONS\n",
+      "fr = 1/(2*math.pi*(math.sqrt(L*C*10**-3*10**-6)));\t\t\t#resonant frequency in hz\n",
+      "I0 = V/(R);\t\t\t#current at resonance in A\n",
+      "Vr = I0*R;\t\t\t#voltage across resistance in volts\n",
+      "Xl = 2*math.pi*fr*L*10**-3;\t\t\t#inductance in ohms\n",
+      "Vl = I0*Xl;\t\t\t#voltage across inductor in volts\n",
+      "Xc = inv([[2*math.pi*fr*C*10**-6]]);\t\t\t#capacitance in ohms\n",
+      "Vc = I0*Xc;\t\t\t#voltage across capacitor in volts\n",
+      "wr = 2*math.pi*fr\t\t\t#angular resonant frewuency in rad/sec\n",
+      "Q = (wr*L*10**-3)/(R);\t\t\t#quality factor\n",
+      "Bw = (fr/Q);\t\t\t#bandwidth in hz\n",
+      "print \"Thus resonant frequency and current are %4.2f hz and %d A respectively\"%(fr,I0);\n",
+      "print \"voltages across resistance, inductance and capacitance are %d V, %d V and %d V respectively\"%(Vr,Vl,Vc);\n",
+      "print \"bandwidth and quality factor are %3.2f hz and %d respectively\"%(Bw,Q);"
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        "Thus resonant frequency and current are 1591.55 hz and 20 A respectively\n",
+        "voltages across resistance, inductance and capacitance are 200 V, 2000 V and 2000 V respectively\n",
+        "bandwidth and quality factor are 159.15 hz and 10 respectively\n"
+       ]
+      }
+     ],
+     "prompt_number": 31
+    },
+    {
+     "cell_type": "heading",
+     "level": 2,
+     "metadata": {},
+     "source": [
+      "Example 5.32  Page No : 196"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "import math \n",
+      "from numpy.linalg import inv\n",
+      "#Chapter-5, Example 5.32, Page 196\n",
+      "\n",
+      "#INPUT DATA\n",
+      "V = 220.;\t\t\t#applied voltage in volts\n",
+      "f = 50.;\t\t\t#frequency in hz\n",
+      "Imax = 0.4;\t\t\t#maximum current in A\n",
+      "Vc = 330.;\t\t\t#voltage across capacitance in volts\n",
+      "#at resonance condition I0 = 0.4 A\n",
+      "I0 = 0.4\t\t\t#current in A\n",
+      "#CALCULATIONS\n",
+      "Xc = (Vc)/(I0);\t\t\t#capacitive reacmath.tance in ohms\n",
+      "C = inv([[2*math.pi*f*Xc]]);\t\t\t#capacitance in F\n",
+      "#at resonance condition Xc = Xl, hence\n",
+      "L = Xc/(2*math.pi*f);\t\t\t#inductance in henry\n",
+      "R = V/(Imax);\t\t\t#resistance in ohms\n",
+      "print \"Thus resistance, inductance and capacitance are %d ohms, %1.2f H and %g F respectively\"%(R,L,C);"
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        "Thus resistance, inductance and capacitance are 550 ohms, 2.63 H and 3.8583e-06 F respectively\n"
+       ]
+      }
+     ],
+     "prompt_number": 1
+    },
+    {
+     "cell_type": "heading",
+     "level": 2,
+     "metadata": {},
+     "source": [
+      "Example 5.33  Page No : 197"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "import math \n",
+      "\n",
+      "#INPUT DATA\n",
+      "R1 = 5;\t\t\t#resistance of branch1 in ohms\n",
+      "R2 = 2;\t\t\t#resistance of branch2 in ohms\n",
+      "L = 10;\t\t\t#inductance in mH\n",
+      "C = 40;\t\t\t#capacitance in uF \n",
+      "#CALCULATIONS\n",
+      "fr = (1./(2*math.pi*(math.sqrt(L*C*10**-9))))*(math.sqrt(((C*10**-6*(R1)**2)-L*10**-3)/((C*10**-6*(R2)**2)-L*10**-3)));\t\t\t#resonant frequency in hz\n",
+      "print \"Thus resonant frequency is %f hz\"%(fr);\n"
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        "Thus resonant frequency is 240.665502 hz\n"
+       ]
+      }
+     ],
+     "prompt_number": 2
+    },
+    {
+     "cell_type": "heading",
+     "level": 2,
+     "metadata": {},
+     "source": [
+      "Example 5.34  Page No : 197"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "import math \n",
+      "\n",
+      "#INPUT DATA\n",
+      "R = 20;\t\t\t#resistance in ohms\n",
+      "L = 0.2;\t\t\t#inductance in H\n",
+      "C = 100;\t\t\t#capacitance in uF \n",
+      "#resistance will be non-inductive only at reosnant frequency\n",
+      "#CALCULATIONS\n",
+      "fr = (1./(2*math.pi*(math.sqrt(L*C*10**-6))))*(math.sqrt((L-(C*10**-6*(R)**2))/(L)));\t\t\t#resonant frequency in hz\n",
+      "print \"Thus resonant frequency is %2.2f hz\"%(fr);\n",
+      "Rf = (L)/(C*R*10**-6);\t\t\t#non-inductive resistance\n",
+      "print \"Thus value of non-inductive resistance is %d ohms\"%(Rf);"
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        "Thus resonant frequency is 31.83 hz\n",
+        "Thus value of non-inductive resistance is 100 ohms\n"
+       ]
+      }
+     ],
+     "prompt_number": 3
+    },
+    {
+     "cell_type": "heading",
+     "level": 2,
+     "metadata": {},
+     "source": [
+      "Example 5.35  Page No : 198"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "#INPUT DATA\n",
+      "Q = 250;\t\t\t#quality factor\n",
+      "fr = 1.5*10**6;\t\t\t#resonant freq in hertz\n",
+      "\n",
+      "#CALCULATIONS\n",
+      "Bw = (fr)/(Q);\t\t\t#bandwidth in Hz\n",
+      "hf1 = fr+Bw;\t\t\t#half power freq 1\n",
+      "hf2 = fr-Bw;\t\t\t#half power freq 2\n",
+      "print \"Thus bandwidth is %d hz\"%(Bw);\n",
+      "print \"Thus value of half-power frequencies are %g hz and %g hz\"%(hf1,hf2);"
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        "Thus bandwidth is 6000 hz\n",
+        "Thus value of half-power frequencies are 1.506e+06 hz and 1.494e+06 hz\n"
+       ]
+      }
+     ],
+     "prompt_number": 4
+    },
+    {
+     "cell_type": "heading",
+     "level": 2,
+     "metadata": {},
+     "source": [
+      "Example 5.36  Page No : 198"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "import math \n",
+      "\n",
+      "#INPUT DATA\n",
+      "L = 40*10**-3;\t\t\t#inductance in henry\n",
+      "C = 0.01*10**-6;\t\t\t#capacitance in uf\n",
+      "#CALCULATIONS\n",
+      "fr = 1./(2*math.pi*math.sqrt(L*C));\t\t\t#resonant frequency\n",
+      "print \"Thus resonant frequency is %d hz\"%(fr);\n"
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        "Thus resonant frequency is 7957 hz\n"
+       ]
+      }
+     ],
+     "prompt_number": 5
+    },
+    {
+     "cell_type": "heading",
+     "level": 2,
+     "metadata": {},
+     "source": [
+      "Example 5.37  Page No : 198"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "import math \n",
+      "from numpy.linalg import inv\n",
+      "#Chapter-5, Example 5.37, Page 198\n",
+      "\n",
+      "#INPUT DATA\n",
+      "V = 120.;\t\t\t#source voltage in volts\n",
+      "R = 50.;\t\t\t#resistance in ohms\n",
+      "L = 0.5;\t\t\t#inductance in Henry\n",
+      "C = 50.;\t\t\t#capacitance in uF\n",
+      "\n",
+      "#CALCULATIONS\n",
+      "#at Resonance\n",
+      "fr = (1./(2*math.pi*(math.sqrt(L*C*10**-6))));\t\t\t#resonant frequency in hz\n",
+      "I0 = V/R;\t\t\t#current at resonance in A\n",
+      "Vl = (1j)*(I0*L);\t\t\t#voltage developed across inductor in volts\n",
+      "Vc = (-1j)*(I0*L);\t\t\t#voltage developed across capacitor in volts\n",
+      "Q = (inv([[R]]))*(math.sqrt(L/(C*10**-6)));\t\t\t#quality factor\n",
+      "Bw = (fr)/(Q);\t\t\t#Bandwidth in Hz\n",
+      "#given resonance is to occur at 300 rad/sec,then\n",
+      "wr = 300;\t\t\t#wr = (2*math.pi*f*r)------->measured in Hz\n",
+      "#wr = inv(math.sqrt(L*Cn))\n",
+      "Cr = inv([[L*(wr)**2]]);\t\t\t#capacitance required in uF\n",
+      "print \"Thus resonant frequency, current, quality factor and bandwidth are %2.1f Hz, \\\n",
+      "%1.1f A, %d and %2.1f hz respectively\"%(fr,I0,Q,Bw);\n",
+      "print \"New value of capacitance at 300 rad/sec is %g F\"%(Cr)\n",
+      "\n"
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        "Thus resonant frequency, current, quality factor and bandwidth are 31.8 Hz, 2.4 A, 2 and 15.9 hz respectively\n",
+        "New value of capacitance at 300 rad/sec is 2.22222e-05 F\n"
+       ]
+      }
+     ],
+     "prompt_number": 8
+    },
+    {
+     "cell_type": "heading",
+     "level": 2,
+     "metadata": {},
+     "source": [
+      "Example 5.38  Page No : 199"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "import math\n",
+      "#INPUT DATA\n",
+      "Q = 45.;\t\t\t#quality factor\n",
+      "f1 = 600.*10**3;\t\t\t#freq in Hz\n",
+      "f2 = 1000.*10**3;\t\t\t#freq in Hz\n",
+      "#given new resistance is 50% greater than former.let us consider two reismath.tances as R1 = 1 ohm and R2 = 1.5 ohm for ease of calculation.Then\n",
+      "R1 = 1;\t\t\t#resistance in ohm\n",
+      "R2 = 1.5;\t\t\t#resistance in ohm\n",
+      "\n",
+      "#CALCULATIONS\n",
+      "W1 = 2*math.pi*f1;\t\t\t#angular freq 1 in rad/sec\n",
+      "W2 = 2*math.pi*f2;\t\t\t#angular freq 2 in rad/sec\n",
+      "Q = 45;\t\t\t#quality factor\n",
+      "L = (Q*R1)/(W1);\t\t\t#inductance in henry\n",
+      "Q1 = (W2*L)/(R2);\t\t\t#new quality factor\n",
+      "print \"Thus new quality factor is %d\"%(Q1);\n",
+      "\n"
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        "Thus new quality factor is 50\n"
+       ]
+      }
+     ],
+     "prompt_number": 9
+    },
+    {
+     "cell_type": "heading",
+     "level": 2,
+     "metadata": {},
+     "source": [
+      "Example 5.39  Page No : 199"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "import math \n",
+      "from numpy.linalg import inv\n",
+      "\n",
+      "#INPUT DATA\n",
+      "R = 4.;\t\t\t#resistance in ohm\n",
+      "L = 100.*10**-6;\t\t\t#inductance in henry\n",
+      "C = 250.*10**-12;\t\t\t#capacitance in Farads\n",
+      "\n",
+      "#CALCULATIONS\n",
+      "fr = inv([[2*math.pi*math.sqrt(L*C)]]);\t\t\t#resonant frequency in Hz\n",
+      "Q = (inv([[R]]))*(math.sqrt(L/C));\t\t\t#Q-factor\n",
+      "Bw = fr/Q;\t\t\t#bandwidth in Hz\n",
+      "hf1 = fr+Bw;\t\t\t#halfpower freq1 in Hz\n",
+      "hf2 = fr-Bw;\t\t\t#halfpower freq2 in Hz\n",
+      "\n",
+      "print \"Thus resonant freq, Q-factor and new halfpower frequencies are %d hz , %d, %g hz, %g hz respectively\"%(fr,Q,hf1,hf2);\n",
+      "#note:given answers are wrong in textbook.Please check the answers\n",
+      "\n"
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        "Thus resonant freq, Q-factor and new halfpower frequencies are 1006584 hz , 158, 1.01295e+06 hz, 1.00022e+06 hz respectively\n"
+       ]
+      }
+     ],
+     "prompt_number": 10
+    },
+    {
+     "cell_type": "heading",
+     "level": 2,
+     "metadata": {},
+     "source": [
+      "Example 5.40  Page No : 200"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "import math \n",
+      "from numpy.linalg import inv\n",
+      "\n",
+      "#INPUT DATA\n",
+      "R = 10;\t\t\t#resistance in ohm\n",
+      "L = 10**-3;\t\t\t#inductance in henry\n",
+      "C = 1000*10**-12;\t\t\t#capacitance in Farads\n",
+      "V = 20;\t\t\t#voltage in volts\n",
+      "#CALCULATIONS\n",
+      "fr = inv([[2*math.pi*math.sqrt(L*C)]]);\t\t\t#resonant frequency in Hz\n",
+      "Q = (inv([[R]]))*(math.sqrt(L/C));\t\t\t#Q-factor\n",
+      "Bw = fr/Q;\t\t\t#bandwidth in Hz\n",
+      "hf1 = fr+Bw;\t\t\t#halfpower freq1 in Hz\n",
+      "hf2 = fr-Bw;\t\t\t#halfpower freq2 in Hz\n",
+      "print \"Thus resonant freq, Q-factor and new halfpower frequencies are %d hz , %d , %g hz, %g hz respectively\"%(fr,Q,hf1,hf2);\n",
+      "\n"
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        "Thus resonant freq, Q-factor and new halfpower frequencies are 159154 hz , 100 , 160746 hz, 157563 hz respectively\n"
+       ]
+      }
+     ],
+     "prompt_number": 11
+    },
+    {
+     "cell_type": "heading",
+     "level": 2,
+     "metadata": {},
+     "source": [
+      "Example 5.41  Page No : 208"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "import math \n",
+      "\n",
+      "#INPUT DATA\n",
+      "P1 = 1000.;\t\t\t#power1 in watts\n",
+      "P2 = 1000.;\t\t\t#power2 in watts\n",
+      "#CALCULATIONS\n",
+      "#for case(1)\n",
+      "Pt = P1+P2;\t\t\t#total power in watts\n",
+      "phi = math.atan(math.sqrt(3)*((P2-P1)/(P2+P1))*(180/math.pi));\t\t\t#math.since math.tan(phi) = math.sqrt(3)*((P2-P1)/(P2+P1)))\n",
+      "pf = math.cos(phi);\n",
+      "print \"Thus power and powerfactor are %d W ,%d respectively\"%(Pt,pf);\n",
+      "#for case(2)\n",
+      "P3 = 1000;\t\t\t#power3 in watts\n",
+      "P4 = -1000;\t\t\t#power4 in watts\n",
+      "Pt1 = P3+P4;\t\t\t#total power in watts\n",
+      "pf1 = 0;\t\t\t#math.since we cannot perform division by zero in scilab,it doesn't consider it as infinite quantity to yield 90 degree angle and hence powerfactor 0\n",
+      "print \"Thus power and powerfactor are %d W ,%d respectively\"%(Pt1,pf1);\n",
+      "\n"
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        "Thus power and powerfactor are 2000 W ,1 respectively\n",
+        "Thus power and powerfactor are 0 W ,0 respectively\n"
+       ]
+      }
+     ],
+     "prompt_number": 14
+    },
+    {
+     "cell_type": "heading",
+     "level": 2,
+     "metadata": {},
+     "source": [
+      "Example 5.42  Page No : 209"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "import math \n",
+      "\n",
+      "#INPUT DATA\n",
+      "V1 = 400.;\t\t\t#voltage in volts\n",
+      "Z1 = (3.+((1j)*4));\t\t\t#impedance in ohms\n",
+      "#CALCULATIONS\n",
+      "#in star connected system,phase voltage = (line voltage)\n",
+      "Ep = V1/(math.sqrt(3));\t\t\t#voltage in volts\n",
+      "Ip = Ep/Z1;\t\t\t#current in A\n",
+      "ip1 = abs(Ip);\t\t\t#line current in A\n",
+      "theta = math.atan(Ip.imag/Ip.real);\n",
+      "Pt = math.sqrt(3)*V1*ip1*math.cos(theta);\t\t\t#total power consumed in load in W\n",
+      "print \"Thus total power consumed in load is %f W\"%(Pt);\n",
+      "#note:for line current the answer given is 46.02A instead of 46.2 A and hence total power consumed changes\n",
+      "\n"
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        "Thus total power consumed in load is 19200.000000 W\n"
+       ]
+      }
+     ],
+     "prompt_number": 15
+    },
+    {
+     "cell_type": "heading",
+     "level": 2,
+     "metadata": {},
+     "source": [
+      "Example 5.43  Page No : 209"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "#INPUT DATA\n",
+      "V1 = 400;\t\t\t#voltage in volts\n",
+      "Il = 10;\t\t\t#current in A\n",
+      "#CALCULATIONS\n",
+      "#in star connected system,phase current = (line current) = I1\n",
+      "phase_voltage = (V1)/(math.sqrt(3));\t\t\t#voltage in Volts\n",
+      "print \"Thus phase voltage is %1.0f V\"%(phase_voltage);\n",
+      "\n"
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        "Thus phase voltage is 231 V\n"
+       ]
+      }
+     ],
+     "prompt_number": 16
+    },
+    {
+     "cell_type": "heading",
+     "level": 2,
+     "metadata": {},
+     "source": [
+      "Example 5.44  Page No : 209"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "import math \n",
+      "#Chapter-5, Example 5.44, Page 209\n",
+      "\n",
+      "#INPUT DATA\n",
+      "Z1 = (6-((1j)*8));\t\t\t#impedance1 in ohms\n",
+      "Z2 = (16+((1j)*12));\t\t\t#impedance2 in ohms\n",
+      "I1 = (12+((1j)*16));\t\t\t#current in A\n",
+      "#CALCULATIONS\n",
+      "V = I1*Z1;\t\t\t#applied voltage in volts\n",
+      "I2 = V/(Z2);\t\t\t#current in other branch in A\n",
+      "print \"current in other branch in Amps\"\n",
+      "print (I2);\n",
+      "I = I1+I2;\t\t\t#total current in A\n",
+      "print \"total current in Amps\";\n",
+      "print (I);\n",
+      "i1 = abs(I);\t\t\t#magnitude in A\n",
+      "i2 = math.atan(I.imag/I.real);\n",
+      "P = V*i1*math.cos(i2);\t\t\t#power consumed in circuit\n",
+      "print \"Thus voltage applied and power consumed are %d V and %d W respectively\"%(V.real,P.real);\n",
+      "\n"
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        "current in other branch in Amps\n",
+        "(8-6j)\n",
+        "total current in Amps\n",
+        "(20+10j)\n",
+        "Thus voltage applied and power consumed are 200 V and 4000 W respectively\n"
+       ]
+      }
+     ],
+     "prompt_number": 19
+    },
+    {
+     "cell_type": "heading",
+     "level": 2,
+     "metadata": {},
+     "source": [
+      "Example 5.45  Page No : 210"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "import math \n",
+      "\n",
+      "#INPUT DATA\n",
+      "Vl = 415.;\t\t\t#voltage in volts\n",
+      "Z = (4+((1j)*6));\t\t\t#impedance in each phase in ohm\n",
+      "#CALCULATIONS\n",
+      "Ip = Vl/Z;\t\t\t#current in each phase in A\n",
+      "ip1 = abs(Ip);\t\t\t#magnitude of Ip\n",
+      "Il = (math.sqrt(3))*(ip1);\t\t\t#line current in A\n",
+      "phi = math.atan(Ip.imag/Ip.real)\n",
+      "P = (math.sqrt(3))*Vl*Il*math.cos(phi);\t\t\t#power supplied in W\n",
+      "print \"Thus power supplied is %d W\"%(P);\n",
+      "#note:the math.cosfunction of scilab and calculator will differ slightly\n",
+      "\n"
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        "Thus power supplied is 39744 W\n"
+       ]
+      }
+     ],
+     "prompt_number": 20
+    },
+    {
+     "cell_type": "heading",
+     "level": 2,
+     "metadata": {},
+     "source": [
+      "Example 5.46  Page No : 210"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "#INPUT DATA\n",
+      "Vl = 400;\t\t\t#voltage in volts\n",
+      "Il = 20;\t\t\t#current in A\n",
+      "f = 50;\t\t\t#freq in hz\n",
+      "pf = 0.3\t\t\t#power factor\n",
+      "#CALCULATIONS\n",
+      "Ip = Il/math.sqrt(3);\t\t\t#phase current in A\n",
+      "Z = Vl/Ip;\t\t\t#impedance in each phase in ohms\n",
+      "phi = math.acos(0.3);\t\t\t#angle in radians\n",
+      "Zb = Z*(math.cos(phi)+(1j)*math.sin(phi));\t\t\t#impedance connected in each phase\n",
+      "print \"Thus impedance connected in each phase in ohms\";\n",
+      "print (Zb);\n",
+      "\n"
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        "Thus impedance connected in each phase in ohms\n",
+        "(10.3923048454+33.0454232837j)\n"
+       ]
+      }
+     ],
+     "prompt_number": 21
+    },
+    {
+     "cell_type": "heading",
+     "level": 2,
+     "metadata": {},
+     "source": [
+      "Example 5.47  Page No : 210"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "import math \n",
+      "\n",
+      "#INPUT DATA\n",
+      "P1 = 6*10**3;\t\t\t#power in Kw\n",
+      "P2 = -1*10**3;\t\t\t#power in Kw\n",
+      "#CALCULATIONS\n",
+      "P = P1+P2;\t\t\t#total power in Kw\n",
+      "a = math.atan(math.sqrt(3)*((P2-P1)/(P2+P1)));\n",
+      "pf = math.cos(a);\t\t\t#power factor\n",
+      "print \"Thus power and power factor are %d W and %1.2f respectively\"%(P,pf);\n",
+      "\n"
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        "Thus power and power factor are 5000 W and 0.28 respectively\n"
+       ]
+      }
+     ],
+     "prompt_number": 22
+    },
+    {
+     "cell_type": "heading",
+     "level": 2,
+     "metadata": {},
+     "source": [
+      "Example 5.48  Page No : 211"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "import math \n",
+      "\n",
+      "#INPUT DATA\n",
+      "Z = 3-((1j)*4);\t\t\t#impedance in ohms\n",
+      "Vl = 400;\t\t\t#line voltage in volts\n",
+      "#CALCULATIONS\n",
+      "Vp = Vl/(math.sqrt(3));\t\t\t#phase voltage in volts\n",
+      "Ip = Vp/abs(Z);\t\t\t#phase current in Amps\n",
+      "#line current(Il) = phase current(Ip)\n",
+      "Il = Ip;\t\t\t#line current in A\n",
+      "power_factor = math.cos(math.atan(Z.imag/Z.real));\n",
+      "power_consumed = math.sqrt(3)*Vl*Il*power_factor;\n",
+      "print \"Thus power consumed and power factor are %f W and %1.1f respectively\"%(power_consumed,power_factor);\n",
+      "#note:answer computed for power consumed in textbook is wrong.Please check the calculations\n",
+      "\n"
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        "Thus power consumed and power factor are 19200.000000 W and 0.6 respectively\n"
+       ]
+      }
+     ],
+     "prompt_number": 23
+    },
+    {
+     "cell_type": "heading",
+     "level": 2,
+     "metadata": {},
+     "source": [
+      "Example 5.49  Page No : 211"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "import math \n",
+      "\n",
+      "#INPUT DATA\n",
+      "Il = 10.;\t\t\t#current in Amps\n",
+      "Vl = 400.;\t\t\t#line voltage in volts\n",
+      "#CALCULATIONS\n",
+      "Vp = Vl/(math.sqrt(3));\t\t\t#line to neutral voltage\n",
+      "Ip = Il;\t\t\t#phase current in Amps\n",
+      "print \"Thus line to neutral voltage  and phase current are %1.0f V and %d A respectively\"%(Vp,Ip);\n",
+      "\n"
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        "Thus line to neutral voltage  and phase current are 231 V and 10 A respectively\n"
+       ]
+      }
+     ],
+     "prompt_number": 24
+    },
+    {
+     "cell_type": "heading",
+     "level": 2,
+     "metadata": {},
+     "source": [
+      "Example 5.50  Page No : 211"
+     ]
+    },
+    {
+     "cell_type": "code",
+     "collapsed": false,
+     "input": [
+      "import math \n",
+      "\n",
+      "#INPUT DATA\n",
+      "P1 = 2000;\t\t\t#power in watts\n",
+      "P2 = 1000;\t\t\t#power in watts\n",
+      "Vl = 400;\t\t\t#line voltage in volts\n",
+      "#CALCULATIONS\n",
+      "P = P1+P2;\t\t\t#power in Watts\n",
+      "a = math.sqrt(3*(P1-P2)/(P1+P2));\n",
+      "b = math.atan(math.sqrt(a));\n",
+      "power_factor = math.cos(b);\n",
+      "kVA = P/power_factor;\n",
+      "print \"Thus power, power factor and kVA are %d W , %1.3f and %1.2f respectively\"%(P,power_factor,kVA);\n",
+      "#note:computed value for powerfactor and kVA in textbook are wrong.Please check the calculations"
+     ],
+     "language": "python",
+     "metadata": {},
+     "outputs": [
+      {
+       "output_type": "stream",
+       "stream": "stdout",
+       "text": [
+        "Thus power, power factor and kVA are 3000 W , 0.707 and 4242.64 respectively\n"
+       ]
+      }
+     ],
+     "prompt_number": 26
+    }
+   ],
+   "metadata": {}
+  }
+ ]
+}
\ No newline at end of file