{
 "metadata": {
  "name": "",
  "signature": "sha256:885ead9482b851b5c86ff5d94634cc23f6d98ffd24283263729cfd9b231402c7"
 },
 "nbformat": 3,
 "nbformat_minor": 0,
 "worksheets": [
  {
   "cells": [
    {
     "cell_type": "heading",
     "level": 1,
     "metadata": {},
     "source": [
      "Chapter 8: Electric Traction Motors"
     ]
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 8.1, Page 237"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "import numpy as np\n",
      "import matplotlib.pyplot as plt\n",
      "%matplotlib inline\n",
      "\n",
      "#Variable declaration\n",
      "v=230;#in volts\n",
      "rm=0.3;#in ohms\n",
      "Ia=np.array([5,10,15,20,25,30,35,40]);#in amperes\n",
      "T=[20,50,100,155,215,290,360,430];#\n",
      "N = []\n",
      "#Calculations&Results\n",
      "for i in range(0,8):\n",
      "    eb= v-(Ia[i])*rm;#\n",
      "    N.append((9.55*eb*Ia[i])/(T[i]));#\n",
      "    print \"speed in rpm is for current \",(Ia[i]),\"  amperes \",(round(N[i])),\" RPM\"\n",
      "    \n",
      "plt.plot(Ia,N)\n",
      "plt.xlabel(\"ARMATURE CURRENT ,Ia IN AMPS\")\n",
      "plt.ylabel(\"SPEED ,N IN RPM\")\n",
      "plt.title(\"Speed-Armature current characteristic\")\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "speed in rpm is for current  5   amperes  546.0  RPM\n",
        "speed in rpm is for current  10   amperes  434.0  RPM\n",
        "speed in rpm is for current  15   amperes  323.0  RPM\n",
        "speed in rpm is for current  20   amperes  276.0  RPM\n",
        "speed in rpm is for current  25   amperes  247.0  RPM\n",
        "speed in rpm is for current  30   amperes  218.0  RPM\n",
        "speed in rpm is for current  35   amperes  204.0  RPM\n",
        "speed in rpm is for current  40   amperes  194.0  RPM\n"
       ]
      },
      {
       "metadata": {},
       "output_type": "pyout",
       "prompt_number": 2,
       "text": [
        "<matplotlib.text.Text at 0xb14178ac>"
       ]
      },
      {
       "metadata": {},
       "output_type": "display_data",
       "png": "iVBORw0KGgoAAAANSUhEUgAAAYcAAAEZCAYAAAB8culNAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XucE9X9//HXchdQ8Y4glxWKiDdQwHol3kCsipdy0bZK\nqdZfvUBVvICtrLeiWG+taPXrBURRFrwUrSCiRsVyEQUVkCroiquAgKBQlC67+f3xmTGzs0k2WZJM\nsnk/H488ksxMZj6Z3cxnzjkz54CIiIiIiIiIiIiIiIiIiIiIiIiIiIhI4MLA74IOQjKuDDgp6CDS\n5EHgT3X4XHtgM1CU3nBE0utY4N/AJmADMAfoGUAcbwDDalmmGKgCHsh8ODVUAfsHsN1ctCP74nPg\nxDTGUldhsncyUkZufOe80SDoAIRdgJeA+4DdgLbATcC2IINK4AJgCTAYaJJguUYZ2v6OnOkF9f/e\n0Pe+iPScsQZ91ruj3yOyg9tP5e8ZIfj9JZKSnsDGBPOHAu8Af8dKFh9T/QxoV+BR4GugHLiF6j+a\nYcAy4FtgJlacdp0CLHfW+3dqP5MrAlYA5wCfAOf65lcBlwKfAiuBPk5M1wDfODGeBZzmfH4DcL3n\n872Budj++NqJqbEz7y1n/Vuw6oBBzr55O0YM7hn1BKza4WXncycCbYBnnXg+A65I8H13Au7Czjo3\nOdtqBoSAL33LlhH9u5QA04BJwHfYPg0Dt2F/y61OjF2BV539sBwY6FnfBGA8duLwPTDP8738+8L7\nOa+Lsb/998BSoLsz/XPgauAD53s9AzR15rVytvkN9j/zInbC4goDt3q+Ryfgt57trAR+74tjALDY\n2RcrgH7OvtgO/OB8h785y9a2T7x/z5Ocabc48/d0Yt/ofP4t7H92ElDpxLsZGAl0xPah+1vZHXgc\n+Mr53s8jErCdgfXYP/mpWOnBayhQAYzAzkAHYT/oVs7857EfzE7AXsB8oj/OAdiB+gDsR3AD9qMG\n+yF9jx3oGwJ/dLaTqFrpOOzH1Qz7cU/3za8CXnFia4odRCuwOuGGwEXOd30KaAF0w36wHZzPH44l\niAbOtGXO9/au31uVMpTak8Mm4Cjn/U7Ae048jbAqspVA3zjfdzzwOrCvE9PPsdJSiJrJwVtVUwL8\nDzjTed8MO6iWAQc669rVWceFzvvuwDpnvhv7euzkoSHwJPB0nO8Zy0AsMR/hvO9E9MSgDEs2rbH/\nt2XAJc683YGznZhbAqVUP1D6v0cjLNkXO/OPB/4L9HDe98b+Bm4bRxvs/xFqVmO2oPZ94v17NsUO\n6Dc778div4WGzuMYz7r9VWkdqZ4c/oXt312d73QcIjmgK/ZP/iV2MP0nsLczbyh2NuM1H/g1sA/w\nI/ZDdp2HHdAAZlD9x9cA++G2x6qH/u1b75ckTg6PAJOd192xA+BenvlV2IHTFcIO/m5xfmdnmV6e\nZRZiSSyWPwLP+dafanKY4Jl3JPCFb/lRwGMxtt3Aif2QGPNC1J4cwr75bzjTXYOxM1uvh4AbndcT\ngIc98/pjpUZXbcnhFeKXij4Hzve8vwM7qMbSHTuTdvm/RyzPA8Od1w9hpa9Y3qB6STWZfTLBN9+b\nHG4CXsASoV+i5LAvVrLYNU6cBUltDrlhOVY0bwccjJ1d3euZ708OXzjLtMeqXVZjRemNwD+IHrA7\nYG0Z7rwNzvS22A+i3Lde7wFvKVZK2Iydge0E/BKY6sxfjJ1Beg8y/nXgbNOtW/7BeV7rmf8DdsYI\n0AWrFliNVUHcBuxB3UWo/h07YPtto+cximgi9toTS7or67ht/76F6vumA5asvLGcjyV8N3b/fmqZ\nwvb3I3Hsa+Ksuzl2QC7D/gZvYgdNb329/2/cHyuJbMC+x2lE/261xeFtd0hmn/i3jSe2O7Fqq1nO\nNq9LsF2vdlgC/C7J5QuCkkPu+Q8wEUsSrra+ZTpgCeNLrOF6D6x6YDfsh+ye7a7Cqph28zxaYPX6\nq7EfhavI9/4g7Ex/Z6wq6mys8fwh57Pu5y/0xbYjjYwPYlUcnZ3vcQOJ/0f/ix3MXK1jLOONZxV2\nBundH7sAp8f43HqsVNY5ie02pHoJyr/deLG86YtlZ+CyGJ+riy+JHXttrsaSdG/sb9CHmg3P3u/R\nFGvDGYcl2d2wNgF3+URx+PfRju6TLVh7QiesSu8q4IQ42/L6EqtOU8nBQ8kheAdg/8RuAmiHVQ3N\n9SyzN1ZMb4zVJXfFfoBrsLOku7EfUQPsh3G887l/AKOxun2wf363ge9lLAGcjdWxDif2wdV1Idbw\nfTBwmPM4xnk+OMHnUtESK6lsxb7jH3zz11K9yuAD7Dschp3ll/iW91+dssBZ/7VYSaghFnusy4ar\nsOqmu7FSVkOsrrsJ1pjeDDtDboy1YTSNsQ4/bzwvYQfhXzvraIxVt3WNE7uff1/4PYIdKA931tWZ\n6hcjxNMSK0l8hx0wx8RYxhtbE+exHttn/anehvMoVio+Efv/bEu0zcH/HeqyT7zTTse+ZxHWnlbp\nxBRrW16rsSrYB7D2ssZEf0MFS8kheJuxovR87MxnLvAhdgbnmg/8DGucuwW7Ssi9wukC7MfpXpE0\nlehB/gWsPvkZ7Mf+EXalCNiPeSBwu/O6M3Z/RSxtsR/3vdhVLO7jfewKqAsSfD//GVuiM7iRWDXC\n91h9+zO+5UuwUtVGrIrrE6y+eTZW4nrbt3zE974KO4B0x65UWudsZ5cE8XwEvItVmYzFDjzfYVdl\nPYJVH22henWHf7ve6a4t2EF0CFYKXO2sv4ln2UT7roTq+8JvGlYtNxnbn89R82KHWPHeiyXO9Vib\n1Ixa4tiMnViUYv9/52FtZq53seRwD9aYHCaapO5zYv/W2W5d94k7rTN2pdNmJ/bxWEkEZz1/wvbX\nVTG+x2+w9r7lWCIZjmRcGXawW4SduYH9Y5c70xZhZxuuUdgVNsuJfxVJIRlKzUZXEZG89zlWPPUa\nQzR7e3XDGjobY1cTrEClm6EoOYhIlmXrwFtbXaFrAHatcQVW4liBNYwVsnhVFCIiGZON5BDB6oQX\nYndsuq7AGhQfJXpDVxuqXwJYTs0rdQrNRNQ4JiJZlo3kcAx2t2R/7JK047BLFouxhsHVxL9JBnTW\nLCKSdZnqHM1rtfO8DrtzsjfV69AfwfpvAbtCwXut/X74bgDr1KlTZOXKut6XJCJSsFaSwr0vmS45\nNMeuvwe7+aovdmmg93r6s51pYH31DMEuXSvGLt9c4FmWlStXEolE8vYxZsyYwGNQ/MHHUWixK/7g\nHyS+L6aGTJcc9iHaaVcjrMO1WcATWJVSBLuaye30axl2vfQyrMfGS4lRrfTpp/Czn2U0bhGRgpbp\n5PA50W6CvRLdNPUX5xHXL38J8+bBTjvtSGgiIhJPXt5DcNBBcPnlQUdRN6FQKOgQdojiD04+xw6K\nP9/k48hIkc2bI/TuDSNHwrDaBrUUERGKiooghWN+XiaHSCTCsmXQpw/Mng2HHRZ0SCIiuS3V5JCX\n1UoA3brBffdZ+8N36oVdRCSt8rbk4Lr0Uli7FqZNg6J8/DYiIllQMCUH1z33wKpVcO+9tS8rIiLJ\nycdz7WolB4CyMjjySHj+eTj66GCCEhHJZQVXcgDo2BEefRSGDIF164KORkQk/9WLkoNr9GhYuBBm\nzICGDbMclYhIDivIkoPr5puhogJuuSXoSERE8lu9KjkArFkDRxwBjz0G/frFXUxEpKAUzE1wibz1\nFgwaBO++C+3aJVxURKQgFHS1kuv44+GqqyxB/O9/QUcjIpJ/6mXJAaCqCs4+G4qLdQ+EiIhKDo4G\nDWDCBJg+HaZODToaEZH8Um9LDq7337eG6Tlz4IADMhiViEgOU8nB5/DD4bbbrIO+rVuDjkZEJD/U\n+5KDfQAuvNA65pswQR30iUjhUckhhqIiePBBeO8962ZDREQSy0ZyKAM+BBYBC5xpuwOvAp8As4BW\nnuVHAZ8Cy4G+6QqiRQt49lkYNQoWLUrXWkVE6qdsJIcIEAJ6AL2daddjyaEL8JrzHqAbMNh5PhV4\nIJ0xHnAA3H+/tT9s2pSutYqI1D/Zqlby13OdCUx0Xk8EznJeDwCeBiqwEscKogklLQYPhtNOg6FD\nrS1CRERqylbJYTawELjYmbYPsNZ5vdZ5D9AGKPd8thxom+6A/vpXWL0a7ror3WsWEakfGmVhG8cA\nq4G9sKqk5b75EecRT9rP75s2hdJS6N3bBgk67rh0b0FEJL9lIzmsdp7XAc9j1URrgdbAGmBf4Btn\nma8Ab1d5+znTqikpKfnpdSgUIhQKpRxUhw52Wet559lVTPvsU+tHRETyRjgcJhwO1/nzmb7ivznQ\nENgMtMCuTLoJOBnYANyBNUa3cp67AZOxBNIWq47qTPXSQ8r3OSTy5z/DO+/Aq69qgCARqb9yrcvu\nYqy0AFZKeQoYi13KWgq0xxqeBwHu9UOjgWHAdmAE8IpvnWlNDpWV0LcvHHUU3Hpr2lYrIpJTci05\nZEJakwPAN99YNxsPPQS/+EVaVy0ikhOUHOpozhw491xYsMDaI0RE6hN1n1FHxx4L114LAwfCtm1B\nRyMiEiyVHKqtGM45B9q2tTupRUTqC5UcdkBRETz+OMycCU8/HXQ0IiLBUckhhsWL4ZRT4K234MAD\nM7opEZGsUMkhDbp3h9tvtwbqLVuCjkZEJPtUcoi7ERg2DCoqYNIkDRAkIvlNJYc0KSqC8ePhww/t\n/gcRkUKSj+fDWSk5uD75BI45BmbMgJ49s7ZZEZG0Uskhzbp0sSFGBw6Eb78NOhoRkexQySFJV14J\nK1bAP/8JDZRSRSTPqOSQIXfcARs2wLhxQUciIpJ5Kjmk4MsvoVcveOYZqMMQEiIigVHJIYPatYMn\nnoDzz7dhRkVE6islhxT17QuXXAJDhsD27UFHIyKSGapWqoPKSjjtNOjRw+6kFhHJdRrPIUvWrYMj\njrDeW888M+hoREQSU3LIorlzYcAAmDcP9t8/6GhEROJTg3QWHXUU3HCD3SD3449BRyMikj7ZSA4N\ngUXAi877EqDcmbYI6O9ZdhTwKbAc6JuF2HbY8OFQXAx//GPQkYiIpE82ksMIYBng1gVFgLuBHs5j\nhjO9GzDYeT4VeCBL8e2QoiJ47DF4/XV48smgoxERSY9MH3z3A04DHiFa11VE7HqvAcDTQAVQBqwA\nemc4vrTYZReYNs262Fi6NOhoRER2XKaTwz3ANUCVZ1oEuAL4AHgUaOVMb4NVN7nKgbYZji9tDj0U\n7rzTBgjavDnoaEREdkyjDK77dOAbrF0h5Jn+IHCz8/oW4C7gd3HWEfOypJKSkp9eh0IhQjnSl8XQ\noTBnDvz+9zB5sgYIEpHghMNhwuFwnT+fycPXX4DfANuBZsAuwLPABZ5lOmIN1YcA1zvT3NvKZgJj\ngPm+9ebMpayx/PCDXcV08cVw2WVBRyMiYnL1Poc+wEjgDGBfwO2Z6EqgF3A+1hA9GWtnaAvMBjpT\ns/SQ08kBrGvvo4+Gl16C3nnRaiIi9V2qySGT1UpeRUQP8uOAw5z3nwOXONOXAaXO83bgUuJUK+W6\nzp1taNFBg+C992CPPYKOSEQkNflYK57zJQfXyJGwbJmVIDRAkIgESXdI55CxY+H77+1ZRCSfqOSQ\nYV99BT172g1yJ50UdDQiUqhUcsgxbdtaYvj1ry1RiIjkAyWHLDjpJLusdcgQqKgIOhoRkdqpWilL\nqqrg9NPhwAPhrruCjkZECo2qlXJUgwYwaRI8+6z1wyQikstUcsiyhQuhf3/rZuOAA4KORkQKhUoO\nOa5nT7jtNjjnHNiyJehoRERiU8khAJEIDBsG27bBU0+pgz4RyTyVHPJAURGMH293T48fH3Q0IiI1\n5eM5a96XHFwrV1oPrv/8pz2LiGRKOntlvRrr+C7WMu5Qn0GoN8kB4MUX4dJLrYO+vfcOOhoRqa/S\nmRyqsNHaZgDbYsy/KaXI0qdeJQeAG26AefPglVegUbb6yRWRgpLO5NAdOA/oB7yPje/8GtWH/AxC\nvUsOlZXQr5+N/fCXvwQdjYjUR5kY7KcIOApLFCcD1wHT6xJcmtS75ACwbh0ccQT8/e8wYEDQ0YhI\nfZOJq5X2AnoAhwLlwLo6RSYJ7bUXlJba8KIrVgQdjYgUukRZ5HfAIKApMA2YCqzNRlC1qJclB9f4\n8fDwwzB3LjRvHnQ0IlJfpLtBegnwRYx5EeDMlCJLn3qdHCIR6967USOYMEE3yIlIeqQzOYSoPoaz\nd9kI8GYqgaVRvU4OAP/9Lxx5JFxxBVxySe3Li4jUJt0N0j2AzlgJ4uM6xtQQWIi1V5wB7A5MAToA\nZVjV1SZn2VHAMKASGA7MirG+ep8cAD75BI45Bl5+GXr1CjoaEcl36WyQvhE7iJ8DvAz8vo4xjQCW\nES2FXA+8CnTBLo293pneDRjsPJ8KPFBLfPValy7w0EMwcCCsXx90NCJSaBIdfIcQvdehJ3VLDvsB\npwGPEM1YZwITndcTgbOc1wOweykqsBLFCqB3HbZZb5xzjiWHX/3K7oUQEcmWRMlhG7DVeb2hlmXj\nuQe4huo3zu1D9Kqntc57gDZY1ZOrHGhbh23WK2PHwo8/ws03Bx2JiBSSRJ017A+8GOd9MlcrnQ58\nAyzCGrdjiVC90TvW/BpKSkp+eh0KhQiF4q0+/zVqBFOm2DgQRx4Jp50WdEQikg/C4TDhcLjOn6/t\naqV4krla6S/Ab4DtQDNgF+A5oJez7jXAvsAbQFeibQ+3O88zgTHAfP+2C6FB2m/OHDj3XOuDqbg4\n6GhEJN9kovuMWHoDC1JYvg8wErtaaRxWTXUHlhBaOc/dgMnOutsCs7ErpfyZoCCTA8C999o41O+8\nA82aBR2NiOSTdF6t1AA4F7gWa1QGa5ieBTxch9jcI/rtwCnAJ8CJREsKy4BS53kGcCmJq5wKzogR\n0Lmz3f8gIpJJibLII0AxVkLoA6zGqn9uAF7IfGhxFWzJAWDzZuu99ZprbKhREZFkpLNaaQnW2V4V\n1mawBuiEVQkFqaCTA8DHH8Pxx8OsWdCjR9DRiEg+SGe1UgXRS1B/BD4n+MQgwIEHwv33WwP1xo1B\nRyMi9VGiLPIDdiOaqxOw0nkdwUoVQSj4koPrj3+07r2nT4cGBXsvuYgkI53VSh1r+WxZshtJMyUH\nR0UFnHAC9O9vQ42KiMSTrUtZg6Tk4PHVV9Yx38SJcMopQUcjIrlKyaEAhcMwZAgsWADt2wcdjYjk\nokwMEyo5LhSCq66yTvq2bQs6GhGpD1RyqCciEevFtU0bG2pURMQrndVKb8SZ7h6ZT0x2I2mm5BDH\nd99Z+8ONN9pQoyIirnQmh56e1+7R+OfAdVhvqz1rfCI7lBwS+OgjOPFEeP11OOSQoKMRkVyRqQbp\nEPAnYCfgVqzvo6AoOdTiySdt/Id334Vddw06GhHJBelODqdifSn9D0sK8aqasknJIQmXXQZffw3P\nPQdF+diyJCJplc7k8C6wF/BXYK4zzXtUfj/V4NJEySEJ27ZZ/0vnngvXXht0NCIStHQmh7DzHO9I\nfEKyG0kzJYckffmlNVA/84xd7ioihUs3wUk1r74KF14ICxfaZa4iUpjSmRzOJfFgO88lu5E0U3JI\n0W23wYwZ8MYb0Lhx0NGISBDSmRwmkDg5/DbZjaSZkkOKqqpgwAAbRe6ee4KORkSCoGoliWnjRujZ\nE8aOhUGDgo5GRLJNyUHiWrQI+vaFt96yAYNEpHDkUsd7zYD5wGJgGTDWmV4ClAOLnEd/z2dGAZ8C\ny4G+GYytIPXoAePGWR9MmzcHHY2I5LJMlxyaA1uBRsAcYCRwErAZuNu3bDdgMtALaAvMBroQHarU\npZLDDrr4Yvj+e7vEVTfIiRSGVEsOjRLM60P1Bml3pe60t5JY/1bnuQnQEHBHPI4V4ADgaWzs6jJs\niNLewLwktiMp+Pvf4Zhj4L77bKhRERG/RMnhGmJfrXQosB92sK9NA+xO6k7Ag8BS4JfAFcAFwELg\namAT0IbqiaAcK0FImjVrBtOmwc9/bo3Uxx4bdEQikmsSJYfTfe+PAf4MrAYuT3L9VUB3YFfgFawD\nvweBm535twB3Ab+L8/mY9UclJSU/vQ6FQoR0+2/KiothwgQbQW7hQmjdOuiIRCSdwuEw4XC4zp9P\npv7pZKxHVoDbgFfruK0/Az9gfTW5OgIvAocA1zvTbneeZwJjsEZtL7U5pNGYMfDmmzB7NjRKdKog\nInktnVcrnY51uHc1dmAPkVpi2BNo5bzeCTgFuzrJe456NvCR83o6MARrnygGfgYsSGF7Ugc33mjV\nTKNHBx2JiOSSRFmkCqv3/yDGvAhwZi3rPgSYiCWgBsAk4E7gCayqKQJ8DlwCrHU+MxoYBmwHRmBV\nUTW2rZJDem3YAEccAXffbZe5ikj9k86b4ELOcyTGchHgzVQCSyMlhwx49134xS9gzhzo0iXoaEQk\n3XSHtNTZQw/B/ffDvHnQokXQ0YhIOmUjOUzE7l8YDyypw+d3lJJDhkQi8NvfwvbtMGmSbpATqU+y\n0X3GeOA17D4FqUeKiuCBB+Cjj+DBB4OORkSClI/nhio5ZNiKFXD00TB9ut0oJyL5L5c63pM81bkz\nPPKIde29bl3Q0YhIEFRykLhGj7armGbOhIbJdJYiIjlLJQdJm5tvtlHkxowJOhIRybbaksOJ2FjR\ny5zHNOCETAcluaFRI3j6aXjiCXjxxaCjEZFsSpQcfgE8ivV9dD7wK+BlZ9ovMh+a5IK994YpU+Ci\ni+Czz4KORkSyJVH905vAcGp2n3EocD9wfKaCqoXaHAJw//3WSD1nDrRsGXQ0IpKqdLY57EPsfpU+\nBPZOLSzJd5ddZuM+dO0Kjz0GlZVBRyQimZQoOWyt4zyph4qKrPTw7LM2DkT37vDyy3ZXtYjUP4mK\nGN8RfyjQ44h2x51tqlYKWCRiDdTXXQdt2sC4cdarq4jkrkz0yhqLemUVtm+3KqaSEjjhBLjtNujY\nMeioRCSWdCaHXbHSQywdgC+SDyutlBxyzJYtcNdd8Le/Wcd9o0fD7rsHHZWIeKWzQTrsef2ab94L\nyYck9V3Llnaj3NKllii6drVksW1b0JGJSF0le4e0zgOlVq1bwz/+YWNSv/22JYnJk+0uaxHJL+o+\nQ9LuwAPhhRdg4kS47z7o3RveeCPoqEQkFY0SzNsLuAqro/K+dueJJHT88Taq3NSpdod1165wxx1w\n8MFBRyYitUlUcngE2Blo6XvdEvi/JNbdDJgPLMb6ZRrrTN8deBX4BJhF9UtiRwGfAsuBvsl+Ccld\nRUXW9ffHH0O/fnDSSZYovvoq6MhEJJFMd9ndHLthrhEwBxgJnAmsB8YB1wG7AdcD3YDJQC+gLTAb\n6AL4a6x1tVIe27TJSg8PPwx/+ANcey3sskvQUYnUf+m8WqnU8/oO37xZSa7fvZO6CdAQ2Iglh4nO\n9InAWc7rAcDTQAVQBqwAeie5HckTrVrB2LGweDGUl0OXLjB+PFRUBB2ZiHglSg4/87z2V/Ek2+bQ\nAKtWWgu8ASzF+mxa68xf67wHaAOUez5bjpUgpB5q18664XjlFbvb+qCD4Lnn1B2HSK5I1CCdDlVA\nd+yGuleoORZExHnEE3NeSUnJT69DoRChUGhHYpQAHXaYjTT36qtwzTV2f8Sdd9oY1iJSd+FwmHA4\nXOfPJ6p/Wo6N41AEPOW8xvO+a4rb+jPwA3AR1jXHGmBfrETRFWt3ALjdeZ4JjMEatb3U5lBPVVXB\nk0/Cn/4EvXpZ9VOXLkFHJVI/pLP7jDDRM/ciap7F1zYi3J7AdmATsBNWcrgJ6AdswNoxrseuVvI2\nSPcm2iDdOcZ2lRzquR9+sK44/vpXGDwYbrzRBh0SkbpLZ3LYUYdgDc4NnMck4E7sUtZSoD3W8DwI\nSyAAo4FhWFIZgSUUPyWHArF+vXXmN2kSXHmlPZo3DzoqkfyUzuTQBTuYd8YG+BkJ5MLV6UoOBWbl\nSrjhBhuF7qabYOhQaNgw6KhE8ks6k8Mc7Mz/beAM4CjgnB0JLk2UHArUggXWaL1hg40h0b+/3WQn\nIrVLZ3JYjF1p5FoE9KhbWGml5FDAIhF46SUbaGjffTXQkEiy0nkTXDPgcOdxBNao7L4+vO4hitRd\nURGccQZ8+KE1Vp9xBvzqV1BWFnRkIvVLslcruct639d2tVKmqOQgP9FAQyLJyaWrlTJFyUFqWLPG\nhit97jmrcrrsMmjWLOioRHJHOquVemE3qbkuBKYDf0OD/0iOiTXQ0FNPaaAhkbpKlEUWAScB3wLH\nA1OAy7FG6a7ALzMeXWwqOUit3nrLrmzavt264zjxxKAjEglWOquVPgAOc16PB9YBJTHmZZuSgyQl\nEoFp02DUKOuGY9w4DTQkhSud1UoNgcbO65OxPpBcme6wT2SHFRXBwIGwbBmceqoNNNSnj3URvnZt\n7Z8XKWSJksPTwJtYO8NW7GY4sK68N8X7kEiuadIEhg+HVatg5EiYOxcOOMCqmh56CNatCzpCkdxT\nWxHjKKA1NrjPf51pXbChQt/PYFyJqFpJdtgPP1hX4aWlMGMG9O5tw5mefTbssUfQ0Ymkny5lFUnR\n1q3w8sswZQrMmmVjSQwaBGedBbvtFnR0Iumh5CCyA7Zsse45SkvhtdfguOPsTuwzz4Rddw06OpG6\nU3IQSZPvv7chTEtLIRyGE06wEsUZZ8DOOwcdnUhqlBxEMmDTJpg+3aqe5syBk0+2RHH66dCiRdDR\nidROyUEkw779Fl54wUoUc+dCv35W9dS/vwYjktyl5CCSRevXw/PPW4li4UI47TQrUZx6qvp2ktyi\n5CASkG++sY7/pkyBxYutymnQIOjbF5o2DTo6KXRKDiI5YM0a67qjtBSWLLGrnQYPtru0mzQJOjop\nROnsPiMd2mHdbiwFlgDDneklQDnWud8ioL/nM6OAT4HlQN8MxyeSEa1bw+WXWweAH30EPXrALbdA\nmzZw0UV2P0VFRdBRisSX6ZJDa+exGLur+j3gLGAQsBm427d8N2Ay1l14W2A2dke2t+NllRwkb61a\nBVOnWonAKXMxAAAPSUlEQVTis8/gnHOs6qlPH2ikHsskg3Kt5LAGSwwAW4CPsYM+xA5yANanUwVQ\nBqwAemc2RJHsad8err4a5s+HBQugUycbnKhtW7j0UhuPorIy6ChFMp8cvDpiY0HMc95fgXX9/SjQ\nypnWBqtucpUTTSYi9UpxMVx7rV3l9M470K4djBhhz8OH2/0UGqxIgpKtgmxLYBowAitBPAjc7My7\nBbgL+F2cz9aoQyopKfnpdSgUIhQKpS9SkQB07mzjTowaBf/5j1U9/eEPsHGjdTs+eDAceaR1Qy6S\njHA4TDgcrvPns/Gv1hh4CZgB3BtjfkfgReAQ4Hpn2u3O80xgDDDfs7zaHKRgLFtm7RNTplgHgQMH\nWhtFr15KFJKaXLuUtQiYCGwArvRM3xdY7by+EmuAPp9og3Rvog3SnaleelBykIITidglsVOnWqLY\nts2SxMCB0LOnEoXULteSw7HAW8CHRA/wo4HzgO7OtM+BS4C1nvnDgO1YNdQrvnUqOUhBi0Ts8tjS\nUntUVFiiGDQIDj9ciUJiy7XkkAlKDiKOSAQ++CBaoohEoiWKHj2UKCRKyUGkQEUi1m2HW6IoKoqW\nKA47TImi0Ck5iAiRCLz/fvSGu0aNoiWKQw9VoihESg4iUk0kAu+9Fy1RNG0aLVEcfLASRaFQchCR\nuCIRePddSxJTp9r4E26iOOigoKOTTFJyEJGkRCLWhYdbothll+h9FN26BR2dpJuSg4ikrKrK+nty\nSxStWkVLFF27Bh2dpIOSg4jskKoqG/506lR77LFHtDH7gAOCjk7qSslBRNKmqgr+/W8rUUybBnvt\nFU0UXboEHZ2kQslBRDKistJ6j5061RJF69bRRNG5c9DRSW2UHEQk4yorrUvx0lJ49lkb4c5NFJ06\nBR2dxKLkICJZVVlpw6GWlsJzz9l4FG6iKC4OOjpxKTmISGC2b6+eKDp0iCaKjh2Djq6wKTmISE7Y\nvh3CYUsUzz9vVU+HHWY32x18sD23bw8NsjkeZQFTchCRnFNRYX09LV1q41K4z99/bzfceRPGwQdb\nIlG3Huml5CAieWPjRhvtzpswli61wYz8CeOgg2DvvZU06krJQUTy3rp1liT8JY2GDS1JeBPGQQfZ\njXqSmJKDiNRLkQisWVM9YbiP5s1rljS6dYNddw066tyh5CAiBSUSgfLy6gljyRL4+GPYbbeaVVPd\nukGLFkFHnX25lhzaAU8Ae2PjRT8M/A3YHZgCdADKgEHAJuczo7AxpCuB4cAs3zqVHESkVlVVUFZW\ns3rqP/+xu7v9JY2uXaFZs6CjzpxcSw6tncdioCXwHnAW8FtgPTAOuA7YDbge6AZMBnoBbYHZQBeg\nyrNOJQcRqbPKSli5smb11IoVdgOfv6TRpQs0aRJ01Dsu15KD3wvA/c6jD7AWSx5hoCtWaqgC7nCW\nnwmUAPM861ByEJG0q6iATz+tWT31xRew//7WI2379nZjX/v20Ue+XEGVanJolLlQaugI9ADmA/tg\niQHneR/ndRuqJ4JyrAQhIpJRjRtbe4R/oKNt26wq6pNPYNUqSxZvv23Pq1bBli1W4nCThT95tGuX\nn9VV2UoOLYFngRHAZt+8iPOIR8UEEQlM06Zw6KH2iGXrVksS3kc4HH1dXm4N47ESh/t+jz1yr/SR\njeTQGEsMk7BqJYhWJ60B9gW+caZ/hTViu/ZzplVTUlLy0+tQKEQoFEpzyCIiyWne3Bqz442YV1Vl\nl+C6yeKLL6x94/XXo+9//DFx8thvv9TbPcLhMOFwuM7fK9O5qgiYCGwArvRMH+dMuwNriG5F9Qbp\n3kQbpDtTvfSgNgcRqVc2b4Yvv4xWVXkTyapV8PXXsOee8ZNH+/ZWOklU+si1BuljgbeAD4ke4EcB\nC4BSoD01L2UdjV3Kuh2rhnrFt04lBxEpKJWVsHp19eThTySVldUThz95FBfnVnLIBCUHERGf776L\nXepwX5eXKzmIiIhPqtVK6kldRERqUHIQEZEalBxERKQGJQcREalByUFERGpQchARkRqUHEREpAYl\nBxERqUHJQUREalByEBGRGpQcRESkBiUHERGpQclBRERqUHIQEZEalBxERKQGJQcREalByUFERGpQ\nchARkRoynRweA9YCH3mmlQDlwCLn0d8zbxTwKbAc6Jvh2EREJI5MJ4fHgVN90yLA3UAP5zHDmd4N\nGOw8nwo8kIX4si4cDgcdwg5R/MHJ59hB8eebTB983wY2xpgea5DrAcDTQAVQBqwAemcssoDk+z+Y\n4g9OPscOij/fBHVmfgXwAfAo0MqZ1garbnKVA22zHJeIiBBMcngQKAa6A6uBuxIsG8lKRCIiknUd\nqd4gHW/e9c7DNRM4MsZnVmBJQw899NBDj+QfK8gxHameHPb1vL4SmOy87gYsBppgJYuVxG6bEBGR\nPPc08DXwP+BLYBjwBPAh1ubwArCPZ/nRWHZbDvTLaqQiIiIiIlJ/lGGljkXAgmBDSUqsmwB3B14F\nPgFmEb1aK9ckcwOj/x6WXNIOeANYCiwBhjvT82X/x4u/hPz4GzQD5mNVxcuAsc70fNn/8eIvIT/2\nP0BDLMYXnff5su/r5HPsC+aL47Ab/bwH2HHAtc7r64Dbsx1UkmLFPga4KphwUtYauyIOoCXwH+BA\n8mf/x4s/n/4GzZ3nRsA84FjyZ/9D7Pjzaf9fBTwFTHfep7Tv8/EO5HxqpI51E+CZwETn9UTgrKxG\nlLxUbmDMRWuwsz6ALcDH2H0z+bL/48UP+fM32Oo8N8HOYjeSP/sfYscP+bH/9wNOAx4hGm9K+z7f\nkkMEmA0sBC4OOJa62gerrsF53ifBsrko1g2Mua4jVgqaT37u/45Y/POc9/nyN2iAJbi1RKvI8mn/\nx4of8mP/3wNcA1R5puXTvk+ZexnsXtgf7bgAY0lWR6pXzfjPxr/NXigp60j12PfGzkKKgFuxH0eu\nawm8R/QsKZ/2P1j8C4nGn49/g12xxHYC+bf/IRp/iPzY/6cD453XIaJtDint+3wrOax2ntcBz5Of\nfS+txeqTwZLdNwHGkqpviN5Q8wi5v/8bA88Ck7DLpiG/9r8b/5NE48+3vwHAd8C/gCPIr/3vcuPv\nSX7s/6OxKqTPsdsJTsR+Aynt+3xKDs2BnZ3XLbAuvePdeZ3LpgMXOq8vJPqjzwfeGxjPJrf3fxF2\nVrcMuNczPV/2f7z48+VvsCfRKpedgFOwK2fyZf/Hi7+1Z5lc3f+jsavdioEhwOvAb8iffZ+yYqwq\naTF2ad+oYMNJiv8mwN9iV1vNJvcvJ0v1BsZccyxW37qY6pcd5sv+jxV/f/Lnb3AI8D4W/4dY/Tfk\nz/6PF3++7H9XH6JXK+XLvhcRERERERERERERERERERERERERERFxnYVdW3+AZ1pH4AfsOvsl2F2h\n7s2TIWf533mW7+5Mu9ozrRF2V7vb7fENRK/dr/S8vgJ4HDjXF9eWJGP5zrOuRdidoX4tgYewQaUW\nYn3m9Cb2cLYlnu8xAfjMWe/7VO+6JYwNULUYmIuNaugqI9rN/CKiN7RNwLp+buK83xO7o/Vgz7Ib\nPNucFeO7eJWRWo/FYeBwz2eneeb9Evs7xHMvFru3A7qh2N/9JM809//pHM823f00B+jiTD+d6D0F\nS4HfJ/81RCQbpmA3zJR4pnUketBsALxG9OAdwg58r3iWvwM7mHm7Ne6PjQn+SYxtbva9f5zowcS/\nTG2xTKd2zwC3ed53xHqv7EDN5ODtntkbV8i37BtED7RDifZlA/G7mZ+AHZT/n/PeTQ5esfZFPKl2\nZ++NuQxLQgc6788lfnJo4Cw7C9sPrguxG8P+zzNtCnbQd7+Dd5sXA//EThy+Bto40xsTTRoSoHzq\nPkMyqyVwJHA5MDjOMlXYIEudnPcR4AugKdEOyfoBM6h+VjkEeBA7qByVRCzJdInsjyWZz3XCSgl/\n8kwrA15O8NmiGK/n+bbrFWterHVHgPuwcdQT/Q7r0j3081ipaAnJ9V4cAe7CSnS1bTOEJYHHgPN8\n897G9m8j7P+pk7NsLG8DnbEucRoS7QSugtgnEZJlSg7iGoCd3a/CqoAOj7FMM+x2/CXOe/cgMg0Y\niB343we2+T5zIpYwSql5QKkrfyxgVT3eaqVi32cOwqouIju47VN924XovvDPK8LOmN2YRnjmrcKq\nVy5IQ0xew7BO4nphI8glU6KYiv3NO9USy3lYieBFrMTV0DMvgo001g/r+C1WSc7dT2dgpc6NznJf\nAJOB88mP8RLqPSUHcZ2HHSBwnr0H8U7YgW0N1jPuy77PTgUGOZ952jfvdKyu+X9YXzRnkfj/LtaB\nyTstUSxvY+MeuA9/NU2ig168ee70IuBObES2Z4A/eJYpwkbc+oyaI4VFsLNtN6b7fPPGYv32pPO3\nOIJo+8d+wM+S+Ewl9v0S9VnWBKsifBH4LzY+hn+YzCnY/8EQav4vuPtpEXYiMdKZfjHWVrHAmfZY\nEvFKhik5CNiZ5QlYL6CfYwerQZ75K7EDWyegK3ZW6rUWO/ifjLUDQPSgeh7Wo+Xn2LgKuxO7odi1\nAdjNF9v6FGJJZBlwGLH/7/3bBdgDK0WBfZ+RWGP9SOBGz3IR7Ix3f6yR/BqStwI7kMeryktVCDvQ\n/hy7OGAxVu1XmwjWrfPxWI+esfTDOmtbgv09jyN6EuGe7b+LNarvAXwaYxvnY3+/c4CvPPOWYA3d\np1DzggQJgJKDgF2d8gTWOFsMtCf64/fagNVL/yXGOm7ExqV1R54qAnbBehd1uw8uxto0ElUthbED\nZWPn/VCsy2G/RLHEsxKri7/JM60jVj2yBSuJnOBM3x07GM7xLOseAO/HvtNRMeb9GSsdtY8xz8+d\nfhvRs+hUvEb1LrzB9vlG4Ecsef48hfVtx0YQu4rYJanzsCvTij2PU7Aurb2ux7qNjsW/L1pQvWG7\nB9YOJAFTchCwKoDnfdOedaa7A5u4XsAan3v75s2leh1zBDtIvoY1MrqmY/XNjT3Lef0Lqx56j2j1\nw3W+9caLxd/mEOtKn4uwbpZXYFccPU506MQLsIP7IifuEqpXTXm3fSs1Sw9gB+X7qF49421zmBDj\nM8uc71tblZpXA6z09K1vuZlYg/AyrMpqbpzPx/Mo1dsRXM2xZPkvz7StWPI8g+r/CzOBN+Os3/99\nirCS1nJs/4zBTghERKQODgL+GnQQIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIpn2/wHZncU8TQuF\nCgAAAABJRU5ErkJggg==\n",
       "text": [
        "<matplotlib.figure.Figure at 0xb561ffac>"
       ]
      }
     ],
     "prompt_number": 2
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 8.2, Page 238"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "import matplotlib.pyplot as plt\n",
      "%matplotlib inline\n",
      "\n",
      "#Variable declaration\n",
      "v=600;#in volts\n",
      "rm=0.8;#in ohms\n",
      "N1=600.;#\n",
      "Ia=[20,40,60,80];#in amperes\n",
      "EMF=[215,381,485,550]\n",
      "T = []\n",
      "N =[]\n",
      "eb =[]\n",
      "#Calculations&Results\n",
      "for i in range(0,4):\n",
      "    eb.append(v-(Ia[i])*rm);#\n",
      "    N.append((N1/EMF[i])*eb[i]);#\n",
      "    T.append((9.55*eb[i]*Ia[i])/(N[i]));#\n",
      "    print \"for current \",(Ia[i]),\" amperes, speed is\",(round(N[i])),\" RPM and Torque in N-m is \",round((T[i]),1)\n",
      "plt.plot(T,N)\n",
      "plt.xlabel(\"TORQUE ,T IN Nm\")\n",
      "plt.ylabel(\"SPEED ,N IN RPM\")\n",
      "plt.title(\"Speed-torque curve\")\n",
      "#answers differ due to rounding-off errors"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "for current  20  amperes, speed is 1630.0  RPM and Torque in N-m is  68.4\n",
        "for current  40  amperes, speed is 894.0  RPM and Torque in N-m is  242.6\n",
        "for current  60  amperes, speed is 683.0  RPM and Torque in N-m is  463.2\n",
        "for current  80  amperes, speed is 585.0  RPM and Torque in N-m is  700.3\n"
       ]
      },
      {
       "metadata": {},
       "output_type": "pyout",
       "prompt_number": 3,
       "text": [
        "<matplotlib.text.Text at 0xb14fb38c>"
       ]
      },
      {
       "metadata": {},
       "output_type": "display_data",
       "png": "iVBORw0KGgoAAAANSUhEUgAAAZAAAAEZCAYAAAC5AHPcAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3X2c1XP+//HHVKKSSui6plKmkqKIJR0kLVa5lhXFknUR\ndr8lF6txsbS2r8vfF/v7VShqxa6wyEXbodUqVxFj0vTdyqQiIqKkmd8fr/dxPnPmzDRz5nzO53PO\nPO+327mdz9V8Pq8ZOq/zvgYRERERERERERERERERERERERERERHJsIZBByASAlHs38J7AcchklUa\nBB2A1GtHAYuBr4EvgX8BAwOIo9y9qlIGdMtQLCJZo1HQAUi9tRfwD2AcMBfYHRgMbA8yqGrkpfhz\njYCf0hlIhjTAEqeISOgMBDZXc34M8AZwP1ZC+Rg41nO+BTAd+AwoBW6lYon6QqAI+AqYD3T2nDse\nKHb3vR+rwrqoijhexz5IvwO+Bc50xy8GVmIlp2eAdp6fKQMuc+dXuWMTPLFeSMVSTeLzxwCLPPsF\nwCvuWcWeGJLZG3gYWIf97k9Xcc9YnLEYHgEeBJ53v+tEYD0V/6anAu+77QbAJKAE2AQ8AbSqJi4R\nkbRpjn3wPAIMp/KHzxhgB3AV1j5xFvaB39Kdfxr7wGsC7AssAS5x50ZgH94HYB90N2DJCGAfYAtw\nmrvv1e45F1YTa2IV1rHAF0B/oDFwH/BawvUvuVh3d7/fBqA30BSYnXDPhQnPH0P8w74Z8Clwgftd\n+rtn96oi1ueBOViCbYSV6hLvmez3egT7+x7h9nfHksNQz/VPYokF7L/LYqA9sBvwkPu9REQyogD7\ntvwp9iH+DLCfOzcG+xbttQQ4D2gDbAP28JwbBfzTbb9IxQ/kBsBWrBRyPvbB5/UptUsg04Epnv1m\nwI/ESzllQMRzfgZwu2e/BzVPIGdjpSCvvwA3JYmzHbATSx6JvPeMSUwgjyScvxX7XcES/ndAJ7df\nRMUSYTvsb6B21XpE/7ElSMXAWOxD6UDs2+w9nvOJCWSNu6Yz9q13PVYNthn7Bryvu64LcK/n3Jfu\neAfsg6404b6ferY/wqqqvgWOrCLudi6WmK3uGR2quGe7hP21Vdw3mS7AIOK/y2bgXCyJJuqEVVt9\nU4v7x5QnxAhWojgNK2WdBrzjuSYfKwXGYirC2nqSxSU5So3oEhYrgEeJV0NBxQ9ksA/TZ7APse1A\na5I39K7Fvj3PSXKuB/Fv0WCN4979PjWI9TPsAzSmmYvFm/C8vbrWU7ENxrsNloCaefbberbXYtVj\nw2oQ16dYG0gLKieRrVj1WbJnVOVjLFH+Ekta3iqqtVjy/3cN7iMiklYHAL8jniQ6Ye0Uf3H7Y7Bq\nrfFYaeNM7EMx1lYyDyutNMdK0t2Bo925kcByrM0B7AM11vAcawM5FfsCdRW7bgNZjzW8xxwHfA70\nw9oK7qViNVNilddwd49e2If4YwnX3IZVYzUB9sfab2LVTc2B1VjV3W7udShW/ZfMP4DHsfaX3Yj/\nTXpi1X79sKq/h6hchXVrkvtNdLF9jyWnmKvd8Vgy3Bc4pYqYRETSqj3Wc6cUq1svxRrF93Tnx2Dj\nQmK9sIqp2KC7F/AA9q37a+BdrKE95jzgAyzprAWmec6dgJV4Yr2wEtsgEo3DSh2bgTM8x0qwqqtn\n3e8Ts5PK40auxZJIKfbN3fvh3RprdN+CJY7JVExIPbHE8DnW8eBV4KAqYm2FJYMNWHXWU55z12MN\n8GuAXyfE+TBwS5L7dXLXPZdwPA+4BvvvsgX7W9xWRUwitTYD2Ih9E4w5DFiKjfh9C/smFXMd9s2r\nmIrF9QHuHiuxb3pSP4yhcqNvLtHgRJFqDAYOpmICiWLf/sDqVRe67d7AMqzInY99m4kN3FqKJR6A\nF7DqAMl9Y1ACEQk1P3thLaLyQLH1xLsYtiTe6DgCa/DcgdX3lmA9T9phdcBL3XUzsfptyX27ml4k\n2+Xy7yb1RKZ7YU3C6rWnYskrNmipPfCm57pSrHF1BxW7XK6jcs8cyU2Puleu0kSmkvUyPQ5kOtar\npjPWADcjw88XEZE0yXQJ5DDiPWmeIt4zZh0V++J3xEoe69y293ji4DIAunfvXr5q1apkp0REpGqr\nsO7jtZbpEkgJMMRtHwt84rafBc7BRrx2xQZ7LcW6Im7B2kPygNFY//9KVq1aRXl5eehfkydPDjyG\nXIkzG2JUnIoz7C9sDFVK/CyBzMGSxT5YX/2bsFHG/4MNvvqB+KjjImxK79h0CJcRb2S8DOvX3gTr\nhTXfx5hFRKSG/Ewgo6o4PqiK47dTccK5mHeAvmmJSERE0kaTKWZYJBIJOoQayYY4syFGUJzppjjD\nI9VV1sKo3NXniYhIDeXl5UGKuUAlEBERSYkSiIiIpEQJREREUqIEIiIiKVECERGRlCiBiIhISpRA\nREQkJUogIiKSEiUQERFJiRKIiIikRAlERERSogQiIiIpUQIREZGUKIGIiEhKlEBERCQlSiAiIpIS\nJRAREUmJEoiIiKRECURERFKiBCIiIilRAknBN9/A8OHw/fdBRyIiEhw/E8gMYCOwPOH4lcDHwIfA\nnzzHrwNWAsXAMM/xAe4eK4F7/Qq2Nlq0gFat4JZbgo5ERCQ4eT7eezDwHTAT6OuOHQNcD5wI7AD2\nBb4AegOzgUOBDsCrQA+gHFgKXOHeXwDuA+YneV55eXm5T79KZRs2QN++sHAhHHhgxh4rIpJWeXl5\nkGIu8LMEsgjYnHDst8AdWPIASx4AI4A57vhqoAQYBLQDmmPJAywZjfQt4lpo29ZKIOPGQVlZ0NGI\niGRepttAegBHA28CUWCgO94eKPVcV4qVRBKPr3PHQ2HcONi5E6ZPDzoSEZHMaxTA81oBh2PVVXOB\nbum6eWFh4c/bkUiESCSSrlsn1aAB/OUvcPzxMGIE7Lefr48TEamzaDRKNBpNy738bAMByAeeI94G\n8iIwBXjN7ZdgyeQ3bn+Ke58PTAbWAAuBXu74KGAIcGmSZ2W0DcRrwgRrE5k1K5DHi4ikLKxtIMnM\nA4512z2BxsAm4FngHLffFavqWgpsALZg7SF5wGh3j1ApLIRFi2DBgqAjERHJHD8TyBxgMZYoPgXG\nYl17u2HdcucA57tri7DqrCKslHIZ1gMLtz0N68ZbQvIeWIFq1gzuvx9++1vYti3oaEREMsPvKqxM\nCqwKK+a006xr7803BxqGiEiN1aUKSwkkjUpLoX9/eOMNOOCAQEMREamRbGoDyWkdO8KNN8Kll0LA\nuUxExHdKIGl2xRU2V5Z6ZIlIrlMVlg/efhtOPhk++ghatw46GhGRqqkNxIQmgQCMH2+z9U6bFnQk\nIiJVUwIxoUogW7ZA794wZw4MHhx0NCIiyakRPYT22gvuucfmy/rxx6CjERFJPyUQH51+OnTtClOn\nBh2JiEj6qQrLZ6tXw8CBsGQJdO8edDQiIhWpCivE8vNtssXLL9fYEBHJLUogGfC738G6dTB3btCR\niIikj6qwMmTxYjjzTBsb0rJl0NGIiBh14zWhTiBgPbIaNoQHHgg6EhERowRiQp9ANm+2sSHz5sGg\nQUFHIyKiRvSs0aqVdekdNw5++inoaERE6kYJJMPOPRf23RfuvTfoSERE6kZVWAFYuRKOOALefRc6\ndw46GhGpz1SFlWV69ICrroIrrww6EhGR1CmBBGTiRFixwhrURUSykaqwAhSNwujRUFQEzZsHHY2I\n1EfqxmuyLoEAjBljvbPuvjvoSESkPlICMVmZQDZtgj594MUX4ZBDgo5GROobNaJnsX32gTvusLEh\nO3cGHY2ISM35mUBmABuB5UnO/R4oA/b2HLsOWAkUA8M8xwe4e6wEcnL0xNix0KQJPPhg0JGIiNSc\nnwnkYWB4kuOdgOOBNZ5jvYGz3ftw4AHiRaoHgYuAHu6V7J5ZLS8PHnoIbr4ZPvss6GhERGrGzwSy\nCNic5PhdwMSEYyOAOcAOYDVQAgwC2gHNgaXuupnASB9iDVzv3laNdfXVQUciIlIzmW4DGQGUAh8k\nHG/vjseUAh2SHF/njuekG26w0ekvvBB0JCIiu9Yog89qClyPVV/FpLUXWGFh4c/bkUiESCSSztv7\nrkkTm+p93DhbN6Rp06AjEpFcE41GiUajabmX391484HngL7u9SrwvTvXEStRDALGumNT3Pt8YDLW\nTrIQ6OWOjwKGAJcmeVZWduNNZtQo6NIFpkzZ9bUiInWRLd14lwNtgK7uVQocgvXUehY4B2jszvXA\n2j02AFuwJJMHjAZyfvKPu++G6dPhww+DjkREpGp+JpA5wGKgJ/Ap8VJGjLe4UATMde8vApd5zl8G\nTMO68ZZgpZOc1rYt3HKLVWWVlQUdjYhIchqJHlJlZfCLX8BFF8HFFwcdjYjkKk1lYnIqgQB88AEM\nHQrLl0ObNkFHIyK5SAnE5FwCAZgwAdavh8ceCzoSEclFSiAmJxPI1q022eK0aVYaERFJJ78SyO+x\nhuxk15RjI8rDJCcTCMBzz8Hvf29VWnvsEXQ0IpJL/EogZcD7WK+o7UnO35zKA32UswkE4LTToG9f\nmy9LRCRd/Eog/bGBeycA72LdchdgiSWMcjqBlJZC//7wxhtwwAFBRyMiucLvNpA84AgsmQwFrsUG\n/oVNTicQgHvvtTXU//lPm8FXRKSu/B6Jvi9wMHAQNnr8i1QeJHV3xRWwZQvMmhV0JCIi1Wedi4Cz\ngN2Bp4AnsWlHwirnSyAAb78NJ59sky22bh10NCKS7fxsRP+Qigs/xZQDp6TyQB/ViwQCMH68de+d\nPj3oSEQk2/mVQCJUnK/Ke2058FoqD/RRvUkgW7bYAlRz5sDgwUFHIyLZzM9G9IOB/bGSyMepPCCD\n6k0CAXjqKbjpJli2DBo3DjoaEclWfjWi3wQ8AZwGvABcksoDxB+nnw5du8LUqUFHIiL1VXVZpwgY\niC0A1Rp4ye2HVb0qgQCsXg0DB8KSJdC9e9DRiEg28qsEsp346oFf7uJaCUB+PkycCJdfDvUsd4pI\nCFSXdb4BXvfsDwYWuW31wgqJHTtgwAC4/no455ygoxGRbONnL6yqqBdWiCxeDGecAUVF0LJl0NGI\nSDYJYjr3w7A1y8Ok3iYQsOVvGzaEBx4IOhIRySZ+JZAGwKlAd6wb7wtYI/rtwH7YZIthUq8TyObN\nNjZk3jwYNCjoaEQkW/iVQKYBXbGSxhBgPVAA3ADMS+VhPqvXCQTg8cfhz3+26U4aNQo6GhHJBn4l\nkA+xCRTLgD2ADVhp5MtUHpQB9T6BlJfDsGEwfLgtQCUisit+JZD3sJHoVe2HTb1PIAAlJXD44fDO\nO9ClS9DRiEjY+ZVAfgBKPPvdgVVuuxwrnYSJEohz663w1lvwzDNaN0REqudXAsnfxc+u3sX5GcBJ\nwOdAX3fsz8DJwI9YMhqLjTcBuA64ENgJjAdedscHAI9g1WgvAFdV8TwlEGf7dujXD+64A049Neho\nRCTMgujGWxODge+AmcQTyPHEl8Wd4o5NAnoDs4FDgQ7Aq0APrKSzFLjCvb8A3AfMT/I8JRCPaBRG\nj7axIc2bBx2NiISV3ysSpmoRsDnh2CvE11RfAnR02yOwNdd3YCWbEmAQ0A5oTnzMyUxgpG8R55BI\nBI47zmbsFRHxQ5DzW12IlSgA2mPL5caUYiWRxOPr3HGpgalTYfZsePfdoCMRkVwU1GiBG7B2kNnp\nvGlhYeHP25FIhEgkks7bZ5199oEpU2yU+ptv2kh1EanfotEo0Wg0Lfeqrt5rYRXHYw0Nx9bg/vnA\nc8TbQADGABcDxwHb3LFJ7j3WLjIfmIwtp7sQ6OWOj8IGNV6aLC61gVRWXm7VWWeeCVdcEXQ0IhI2\ndWkDqa4EMsGzHftkPhy4FutZlYrh7r5DiCcPgGex0shdWBVVD6zdoxzYgrWHLAVGY43oUkN5efDQ\nQ7b07amnQgdVAIpImtQ060SAG4EmwG3AizX4mTlYotgH2IiVKK4DGgNfuWv+DVzmtq/H2kV+wrrq\nvuSOx7rxNsHaTMZX8TyVQKpx442wYgU8+WTQkYhImPjZjXc48faK26i6WisMlECq8cMP0Lcv3Hcf\nnHhi0NGISFj4lUDeAvYFpmIlBYhXZQGErW+PEsguvPyyNah/9BE0bRp0NCISBn4lkKh7r+pT+ZhU\nHugjJZAaGDXK5siaMmXX14pI7gvrSPRMUwKpgQ0brCpr4UI48MCgoxGRoPmVQE6n6tIHwN9TeaCP\nlEBq6KGHYNYsWLQIGgQ5lFREAudXAnmE6hPI2FQe6CMlkBoqK4Mjj4SxY+GSS4KORkSCpCosowRS\nCx98AEOHwvLl0KZN0NGISFCUQIwSSC1NmADr18NjjwUdiYgERQnEKIHU0tat0KcPTJtmpRERqX/C\nOp27hFyzZnD//XDZZbBt266vFxHxqi7rDKFiI3rs2tix132JKHUqgaTo9NOtS+/NNwcdiYhkml9V\nWP8geS+sg7CFoMI2ObgSSIpKS6F/f/jXv6CgIOhoRCSTMtUGciTwB6Al8EdsmvYwUQKpg3vvhaef\ntgGGebnUMiYi1fI7gQzFZuIFSxyvpPKgDFACqYOdO+Gww2D8eLjggqCjEZFM8SuBnIzNxPs1cDu2\nxnmYKYHU0dtvw8kn22SLrVsHHY2IZIJfCaQMW4/8/STnyoFTUnmgj5RA0mD8eOveO3160JGISCb4\nlUAi7r08yXXlwGupPNBHSiBpsGUL9O4Nc+bYKoYikts0kNAogaTJ3/4Gf/gDLFsGjRsHHY2I+CnT\nAwkfBR4ENBl4jjrtNOjWDaZODToSEQmzVLLOYUBn9z4xveHUiUogabR6NQwcCEuWQPfuQUcjIn5R\nFZZRAkmzO++EBQtg/nyNDRHJVZoLS3xxzTU2W+8TTwQdiYiEUS59r1QJxAeLF8MZZ0BREbRsGXQ0\nIpJuqsIySiA+GTcOGjaEBx4IOhIRSTc/q7COxdY+L3Kvp4BjanjvGcBGYLnn2N7YVCifAC9j82rF\nXAesBIqBYZ7jA9w9VgL31vDZkkZTpsC8edagLiISU10COQmYjk2aeC7wa+AFd+ykGtz7YWB4wrFJ\nWALpCSxw+wC9gbPd+3DgAeIZ8UHgIqCHeyXeU3zWqpV16b3kEtixI+hoRCQsqksgE4GRWCJYBryH\nlSpGAtfW4N6LgM0Jx07BxpHg3ke67RHAHGAHsBooAQYB7YDmwFJ33UzPz0gGjRoF++1ns/aKiED1\nCaQNyefB+gDYL8XntcGqtXDvbdx2e2zerZhSoEOS4+vcccmwvDx48EGrzlqzJuhoRCQMGlVz7vsU\nz9VUOckXrEpZYWHhz9uRSIRIJJLO29d7++8PV10FV14JzzyjsSEi2SgajRKNRtNyr+o+Ar6h6mVr\nB1OxAbwq+VgbSl+3X4xN0rgBq55aCBQQbwuZ4t7nA5OBNe6aXu74KGyp3UuTPEu9sDJg+3bo1w/u\nuANOPTXoaESkrurSC6u6EsiIas6lOkvSs8AFwJ/c+zzP8dnAXVgVVQ+s3aMc2IK1hywFRgP3pfhs\nSYPdd4eHHoLRo2HoUGjePOiIRCQo1WWdFlgpJJkuWOmgOnOw0sI+WHvHTcAzwFxsLq3VwFnYglUA\n1wMXAj8BVwEvueMDgEeAJlgvsPFVPE8lkAwaO9YGFt59d9CRiEhd+DWQ8D3gYLe9ADiuinNhoQSS\nQZs2QZ8+8OKLcMghQUcjIqnKxFxYe6dyc8ld++xjPbLGjbP11EWk/tFkipKyMWOgaVNNcSJSX1XX\niL4v8DusaOPdjp2Tei4vzxrUBw+2Rag6aISOSL1SXb1XIfFxGnlUHrNxsx8B1YHaQAJy442wYgU8\n+WTQkYhIbWk2XqMEEpAffoC+feG+++DEE4OORkRqw69G9Lme7T8lnHs5lYdJbmrSxNpBLr8ctm4N\nOhoRyZTqEkgPz/awhHNqA5EKhg2DI46AW24JOhIRyRT1wpK0uesumDEDli/f9bUikv2q64XVBDgE\nqxuLbePZF6mgbVu49VYbG/Kvf0EDfT0RyWnVNZxEqb4XVk1XJswUNaKHQFkZHHmkTXVyySVBRyMi\nu6JeWEYJJCQ++MAmWly+HNq02fX1IhIcv3ph9cQmP/wImxhRw8SkRg46CC64AC68ED75JOhoRMQv\n1SWQGcA/gNOBd4H7MxKR5ISbb4aCAjj6aDj8cOvm++WXQUclIulUXbFlGdDfsx/GGXi9VIUVQj/9\nBK+8AjNn2sy9xxxja4mcdJKtLSIiwfKrDaQYONdz3eNuP9ag/m4qD/SREkjIbdkCTz0Fs2ZZ+8iZ\nZ8L551sJRcvjigTDrwQSpWLPq8SeWOqFJSlbswYef9xKJjt3WqnkvPOgW7egIxOpX9QLyyiBZKHy\ncnj7bUskTzwBBxxgpZIzz7QVD0XEX34lkEOBUmC9278Aa1Bfjc3U+1UqD/SREkiW+/FHmD/fkskr\nr8AJJ1jJZPhw2G23oKMTyU1+Lml7HJYojgaeAK7AGtILgDNSeaCPlEByyObNMHeutZd88gmcc46V\nTAYMUHuJSDr5lUDeB/q57f8BvsBKHonnwkIJJEetWgWPPWYlk913t1LJr38NnTsHHZlI9vNrIGFD\nIFZxMBRY6DlX3RxaImnVvTtMngwlJTBtmjXAH3wwHHssPPIIfPtt0BGK1E/VZZ0bgJOATUAnYABQ\nhk3z/ghwpN/B1ZJKIPXItm3w/PNWKnntNRtXMnq0TaHSSF9vRGrMz15YRwBtsQWkYksF9QT2RONA\nJCQ2bbIeXDNnwtq1cO651l7SL2yVrCIhlI3deK8DzsNKNMuBsUAzrKG+C9bT6yzga8/1FwI7gfEk\nXxFRCURYscIa3mfNghYtLJGcey60bx90ZCLhlG0JJB/4J9AL2I4ljReAPlh12Z3AtUArYBLQG5iN\ndSvuALyKlYLKEu6rBCI/KyuDRYusVPL003DooZZMRo6EZs2Cjk4kPPxqRPfLFmAH0BRrjG8KfAac\nAjzqrnkUGOm2R2CzAe/ASiYlwGGZC1eyUYMGMGQITJ8OpaW2Psns2dCxI4wZAwsW2Ah4EUldEAnk\nK+C/gbVY4vgaeAVoA2x012x0+wDtsQGNMaVoanmphaZNbRzJ88/Dxx9b28iECZCfD5MmQVFR0BGK\nZKcg+qt0B67GqrK+AZ7E2kO8yqm8AmLi+UoKCwt/3o5EIkQikdSjlJzUti1cc429PvzQ2kqOPx7a\ntbNeXKNGwX77BR2liH+i0SjRaDQt9wqiDeRs4HjgN25/NHA4cCw2QeMGoB027qQAawcBmOLe5wOT\ngSUJ91UbiKRk505YuNDaS559Fo46ytpLfvUraNIk6OhE/JVtbSDFWMJoggU9FCgCnsPm28K9z3Pb\nzwLnAI2Brtg4lKUZjFdyXMOGNn5k5kxrLzn7bBuw2KEDXHwxvP66NcqLSEVBdeOdiCWJMmw8yW+A\n5sBcoDOVu/Fej3Xj/Qm4CngpyT1VApG0WrfOGt5nzoTvvrPp5kePhp49g45MJH2yrRuvX5RAxBfl\n5fD++5ZIZs+Grl0tkZx9NrRuHXR0InWjBGKUQMR3yZboPf98OPFELdEr2UkJxCiBSEYlLtF71llW\nMtESvZJNlECMEogEJnGJ3vPPtzaTrl2DjkykekogRglEApe4RG9BgZVKtESvhJUSiFECkVDxLtH7\n6qswbJiVTE44QUv0SngogRglEAkt7xK9K1fa1CqjR2uJXgmeEohRApGskGyJ3vPOg06dgo5M6iMl\nEKMEIlmlvBwWL7ZSyZNPQv/+lkxOPx2aNw86OqkvlECMEohkrWRL9J5/Phx3nJboFX8pgRglEMkJ\nmzbBX/9qJZNPP40v0XvQQUFHJrlICcQogUjOKS629pJZs6wbcGyJ3nbtgo5McoUSiFECkZyVuETv\nYYdZe4mW6JW6UgIxSiBSL3z/va1bMmuWNcKPGGElk0jElvIVqQ0lEKMEIvXOhg0wZ44lk88/h0GD\nbPR7r172XlAAe+4ZdJQSZkogRglE6rUVK2DZMms3KS629d8/+cSmnE9MKr162fK+GsQoSiBGCUQk\nQVkZrF1rySSWVGIJZvv2iokl9t6tm6ZaqU+UQIwSiEgtfPllxdJK7H3dOksi3tJKQQEccADstVfQ\nUUu6KYEYJRCRNNi2zebr8pZWYtVhrVolrw5r107VYdlKCcQogYj4qKzMBjYmJpbiYvjhh+SJpXt3\nVYeFnRKIUQIRCchXX1VOKsXFlnDy8yu3sxQUqDosLJRAjBKISMhs2wYlJZXbWVassJH1ie0svXpB\n+/aqDsskJRCjBCKSJcrKoLQ0eXXY1q1VV4c1bhx05LlHCcQogYjkgM2b40nFm1jWroUuXZJXh7Vo\nEXTU2SsbE0hLYBrQBygHxgIrgSeALsBq4Czga3f9dcCFwE5gPPByknsqgYjksO3bK1eHxV577VW5\nOqygADp2VHXYrmRjAnkUeA2YATQCmgE3AJuAO4FrgVbAJKA3MBs4FOgAvAr0BMoS7qkEIlIPlZXZ\n2JVk1WHffWfjVxLbWfbfX9VhMdmWQFoA7wHdEo4XA0OAjUBbIAoUYKWPMuBP7rr5QCHwZsLPK4GI\nSAVff528OmzNGujcuXI7S0GBNe7XJ9mWQPoDfwGKgH7AO8DVQClW6ojF9ZXbvx9LFo+7c9OAF4G/\nJdxXCUREauTHH6uuDttzz+TVYZ065WZ1WF0SSBCLZTYCDgGuAN4C7sGqqrzK3asqSc8VFhb+vB2J\nRIhEInUIU0RyVePG0Lu3vbzKy606zJtYnnnGtrdsqVwdVlAAPXrA7rsH83ukIhqNEo1G03KvIPJp\nW+DfQFe3fxRWTdUNOAbYALQDFmJVWLHkMsW9zwcmA0sS7qsSiIj45ptvkleHrV5tpZPE6rADD4Tm\nzYOOeteyrQoL4HXgN8AnWHtGU3f8S6ytYxLWU8vbiH4Y8Ub0/alcClECEZGM+/FHWLWqcmK58kpb\nNTLssjGB9MPaMhoDq7BuvA2BuUBnKnfjvR7rxvsTcBXwUpJ7KoGIiNRSNiYQPyiBiIjUUl0SiFZQ\nFhGRlCiFGX+TAAAHpUlEQVSBiIhISpRAREQkJUogIiKSEiUQERFJiRKIiIikRAlERERSogQiIiIp\nUQIREZGUKIGIiEhKlEBERCQlSiAiIpISJRAREUmJEoiIiKRECURERFKiBCIiIilRAhERkZQogYiI\nSEqUQEREJCVKICIikhIlEBERSYkSiIiIpCTIBNIQeA94zu3vDbwCfAK8DLT0XHsdsBIoBoZlMEYR\nEalCkAnkKqAIKHf7k7AE0hNY4PYBegNnu/fhwANkcckpGo0GHUKNZEOc2RAjKM50U5zhEdQHcUfg\nRGAakOeOnQI86rYfBUa67RHAHGAHsBooAQ7LVKDpli3/U2VDnNkQIyjOdFOc4RFUArkbmACUeY61\nATa67Y1uH6A9UOq5rhTo4HeAIiJSvSASyMnA51j7R14V15QTr9qq6ryIiNQztwOfAv8B1gNbgVlY\nA3lbd007tw/WFjLJ8/PzgUFJ7ltCPPHopZdeeulVs1cJWWoI8V5YdwLXuu1JwBS33RtYBjQGugKr\nqLrkIiIi9cQQ4Fm3vTfwKsm78V6PZcli4IRMBigiIiIiIlLBcKxkspJ4FVhQZmA9yJZ7joVxgGQn\nYCHwEfAhMN4dD1usewBLsCrMIuCOkMYJ2TEwdjXwARbnUncsjHG2BJ4CPsb+uw8KYZwHYH/H2Osb\n7N9R2OKMPfcj7HNpNrB7SOPMuIZY1VY+sBv2QdMrwHgGAwdTMYHcCUx029dSuW1nNyz+EjLXK64t\n0N9t7wmswP5uYYy1qXtvBLwJHBXSOH8HPE68SjaMMf4H++DwCmOcjwIXuu1GQIuQxhnTAOsQ1Inw\nxZkP/C+WNACeAC4IYZyBOALrlRWT2GMrCPlUTCDFxMe0tCXeu+w6KpaY5gOH+x1cFeYBQwl3rE2B\nt4A+hC/Ojlj73THESyBhixEsgbROOBa2OFtgH3iJwhan1zBgkdsOW5x7Y18QW2HJ+Dng+HTFme2Z\npQPWJTgmjIMMwz5AMh8rNS0hnLE2wL4RbSRe7Ra2OLNlYGw5lujeBi52x8IWZ1fgC+Bh4F3g/wHN\nQhin1znYbBkQvji/Av4bWAt8BnyNVV2lJc5sTyDlQQdQS7F+19Wdz6Q9gb9h85J9mySWMMRahlW3\ndQSOxr7lJ8YRZJzZNDD2SOzLwi+By7Eq18Q4go6zEXAINufdIdg4scRahTDEGdMY+BXwZBVxBB1n\nd+Bq7Itie+zf/HlJ4kgpzmxPIOuweseYTlTMnmGwkYoDJD9324mxd3THMmU3LHnMwqqwILyxgjVS\nPg8MIFxx/gKbx+0/2LfQY7G/aZhijFnv3r8AnsbmlAtbnKXu9ZbbfwpLJBsIV5wxvwTewf6mEL6/\n50BgMfAl8BPwd6zqP6x/z4xqhA0szMe+CQTdiA6V20DCOEAyD5iJVb14hS3WfYj3DmkCvA4cF8I4\nY8I8MLYp0NxtNwPewOruwxYn2H/nnm670MUYxjgB/oo1SseELc5+WE/LJu55j2Klz7DFGZhfYo1E\nJVgDUJDmYPWMP2JtM2MJ5wDJo7CqoWXEuyEOD2GsfbF68GVY99MJ7njY4owJ88DYrtjfcRn2gRL7\ntxK2OME+9N4C3se+MbcIaZzNgE3EEzOEM86JxLvxPorVPoQxThERERERERERERERERERERERERER\nERHJXa2Jj2FZj41efg8bN9IJeAbr414C3IP1fQeIYKPa38OmCr8t4b4jsfEHRdjYkxGec1FsJHxM\nPvFBo977xl7HVhP/3901K7H5iWI/kzhx3SPA6Z7nv+U5NxCbHyxRPjbO5wrPsf9DxQFvIiICTMam\nU49ZSvzDsgEwDRt5C/ZBHxstvge25kQsKfTDPtC7uP18t9/X7S/EptXAc96bQJ6l9ryj15N5GDjN\nbUexNT2Gu/3qEsgGLIHGEuf9KIFImmX7XFgiMbHpFo4DfsBG3IJ9E78GW19ij4Sf2YaNzO7m9v8L\n+COwxu2vxhax+q9axlAbtfmZcmAqcEMNrv0CWEDFpBGbFC8K3IWVZj4GDsXmxvoEuLUW8Ug9pwQi\nuaYPNrmd17fYdNY9Eo7vjU0oWOT2eyf52XfcPWtiMBWrsLrW8Odq49/YVDkRdj2b651Y8kv8d14O\nbMcSx4NYdd+lwIHAGGztCJFdUgKRXFOTaakHYyWPT7GZiD/axT1j1UDJ7u09tgibLj32+s+ugk3R\nbcCNNbjuP9g6L+cmORerbvvQvTZiiel/gc5piFHqASUQyTVFVGzoBtgL+1AscfuLsDVG+mDtC508\nPzsw4WcHYI3pYFNie5eE3RubTC+TyrF2jybUbEW727FZV/OoWF223b2XebZj+w3rHqbUB0ogkmsW\nYFOXj3b7DbEV2R7G2jy8VgP3An9w+1OxWWq9jejXAve5/SgVF+O5APhnLeObiVUd1dVtWGy7qsZa\ngSXGX9XgWpFaUQKRXOH9cDwVOBNrFF4BfI9NUR27znvtQ1ivpo5Y991rsV5RK9xrEvGus/8Xa095\nH6sCa4olndh9E9tAYr2nvPpScYGeXa0GV5UXiS8ClIz3nn/Efr+qrlNiERFJszuAudjCZemwF/BE\nmu4lIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIlIj/x90qdUelP+JNgAAAABJRU5ErkJggg==\n",
       "text": [
        "<matplotlib.figure.Figure at 0xb131cf4c>"
       ]
      }
     ],
     "prompt_number": 3
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 8.3, Page 239"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "import math\n",
      "\n",
      "#Variable declaration\n",
      "N1=640;# in rpm\n",
      "I1=15;# in A\n",
      "\n",
      "#Calculations\n",
      "I2=math.sqrt((2)*math.sqrt(2)*I1**2);\n",
      "N2=round((2*I1*N1)/I2);\n",
      "\n",
      "#Results\n",
      "print \"Current drawn,I2(A) = %.3f\"%I2\n",
      "print \"Motor speed,N2(rpm) = %.f\"%N2"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "Current drawn,I2(A) = 25.227\n",
        "Motor speed,N2(rpm) = 761\n"
       ]
      }
     ],
     "prompt_number": 91
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 8.4, Page 243"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "#Variable declaration\n",
      "n1=700;#rpm\n",
      "n2=750;#rpm\n",
      "rm=0.3;#in ohms\n",
      "v=500;#in volts\n",
      "ib=50;#amperes\n",
      "\n",
      "#Calculations\n",
      "eb1=v-(ib*rm);#in volts\n",
      "eb2=eb1;#\n",
      "N=((v-(2*(ib*rm)))/((eb1/n1)+(eb2/n2)));#\n",
      "pdv1=((eb1/n1)*N)+ib*rm;#in volts\n",
      "pdv2=((eb1/n2)*N)+ib*rm;#in volts\n",
      "\n",
      "#Results\n",
      "print \"speed in rpm is %.f\"%N\n",
      "print \"PD across machine 1 in volts is %.f\"%pdv1\n",
      "print \"PD across machine 2 in volts is %.f\"%pdv2"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "speed in rpm is 351\n",
        "PD across machine 1 in volts is 258\n",
        "PD across machine 2 in volts is 242\n"
       ]
      }
     ],
     "prompt_number": 92
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 8.5, Page 243"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "import math\n",
      "\n",
      "#Variable declaration\n",
      "V=500;# in volts\n",
      "Vm=40;# in kmph\n",
      "Ft=1800;# in N\n",
      "Rm=0.4;# in ohm\n",
      "Lm=3200;# losses per motor in watt\n",
      "\n",
      "#Calculations&Results\n",
      "Mo=(Ft*Vm*1000)/3600;\n",
      "Cl=3200;# consatant losses in watt\n",
      "# formuls: Mi=Po+Cl+C_losses\n",
      "#C_losses=I^2*Rm\n",
      "#Mi=V*I\n",
      "#I1=(V+sqrt(V^2-(4*Rm*(Mo+Cl))))/(2*Rm);leaving as gives a very high value\n",
      "I1=(V-math.sqrt(V**2-4*Rm*(Mo+Cl)))/(2*Rm);\n",
      "print \"Current drawn by each motor,(A) = %.2f\"%I1\n",
      "It=I1*2;\n",
      "print \"Total current drawn,(A) = %.1f\"%It"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "Current drawn by each motor,(A) = 48.26\n",
        "Total current drawn,(A) = 96.5\n"
       ]
      }
     ],
     "prompt_number": 93
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 8.6, Page 244"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "import math\n",
      "\n",
      "#Variable declaration\n",
      "Ft=35300.;# in N\n",
      "V=48.;# in kmph\n",
      "\n",
      "#Calculations&Results\n",
      "Po=((Ft*V*1000.)/3600)*10**-3;\n",
      "Ft1=55180.;#in N\n",
      "Pd=Po*math.sqrt(Ft1/Ft);\n",
      "print \"part (a)\"\n",
      "print \"power delivered(kW) = %.1f\"%Pd\n",
      "Pd1=Po*(Ft1/Ft);\n",
      "print \"part (b)\"\n",
      "print \"power delivered(kW) = %.1f\"%Pd1"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "part (a)\n",
        "power delivered(kW) = 588.5\n",
        "part (b)\n",
        "power delivered(kW) = 735.7\n"
       ]
      }
     ],
     "prompt_number": 97
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 8.7, Page 244"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "#Variable declaration\n",
      "Ia=[60,120,180,240,300,360];# in amperes\n",
      "sp1=[80,50,45,42,38,35];#in kmph\n",
      "tf1=[1.7,5,10,14,16,20];#innewtons\n",
      "d1=0.85;#in meters\n",
      "d2=0.9;#in meters\n",
      "y1=71./21;#\n",
      "y2=74./19;#\n",
      "\n",
      "#Calculations&Results\n",
      "for i in range(0,6):\n",
      "    V=((d2/d1)*(y1/y2))*sp1[i];#in kmph\n",
      "    tf2=((d1/d2)*(y2/y1))*(tf1[i]);#in newtons\n",
      "    print \"for armature current\",(Ia[i]),\"amperes , speed is \",(V),\" kmph and tractive effor in thousand newtons is \",round((tf2),2)\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "for armature current 60 amperes , speed is  73.5316829434  kmph and tractive effor in thousand newtons is  1.85\n",
        "for armature current 120 amperes , speed is  45.9573018397  kmph and tractive effor in thousand newtons is  5.44\n",
        "for armature current 180 amperes , speed is  41.3615716557  kmph and tractive effor in thousand newtons is  10.88\n",
        "for armature current 240 amperes , speed is  38.6041335453  kmph and tractive effor in thousand newtons is  15.23\n",
        "for armature current 300 amperes , speed is  34.9275493981  kmph and tractive effor in thousand newtons is  17.41\n",
        "for armature current 360 amperes , speed is  32.1701112878  kmph and tractive effor in thousand newtons is  21.76\n"
       ]
      }
     ],
     "prompt_number": 98
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 8.8, Page 245"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "import numpy as np\n",
      "\n",
      "#Variable declaration\n",
      "n1=500;#in rpm\n",
      "d1=90;#in cm\n",
      "d2=86;#in cm\n",
      "v=600;#in. volts\n",
      "vd=0.1;#drop\n",
      "\n",
      "#Calculations\n",
      "eb1=v-(vd*v);#in volts\n",
      "A = np.array([[90,-86],[1,1]])\n",
      "B = np.array([240,600])\n",
      "Eb1 = np.linalg.solve(A, B)\n",
      "N1=n1*(Eb1[0]-(vd*v))/(v-(vd*v));#\n",
      "N2=N1*(d1/d2);#\n",
      "\n",
      "#Results\n",
      "print \"speed in rpm is %.f\"%N1\n",
      "print \"speed in rpm is %.f\"%N2\n",
      "#N2 is calculated wrong in the book'''"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "speed in rpm is 217\n",
        "speed in rpm is 217\n"
       ]
      }
     ],
     "prompt_number": 83
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 8.9, Page 245"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "#Variable declaration\n",
      "ia=350;#A\n",
      "ib=305;#A\n",
      "v=600;#V\n",
      "\n",
      "#Calculations&Results\n",
      "pa=(v*ia)/1000;#kW\n",
      "pb=(v*ib)/1000;#kW\n",
      "print \"(i) When motors are connected in parallel and train speed is 40kmph\"\n",
      "print \"power input to motor A is,(kW)= %.f\"%pa\n",
      "print \"power input to motor B is,(kW)= %.f\"%pb\n",
      "fta=1625;#kg\n",
      "ftb=1480;#kg\n",
      "print \"tractive effor of motor A is,(kg)= %.f\"%fta\n",
      "print \"tractive effor of motor B is,(kg)= %.f\"%ftb\n",
      "print \"(ii) When motors are connected in series and current is 400A\"\n",
      "rm=0.08;#ohm\n",
      "i=400;#A\n",
      "eba=v-(i*rm);#V\n",
      "abb=eba;#V\n",
      "va=38.5;#V\n",
      "vb=36.7;#V\n",
      "vx=((v-2*(i*rm))*((va*vb)/(va+vb)))/eba;#\n",
      "Va=((eba/va)*vx)+(i*rm);#V\n",
      "Vb=((eba/vb)*vx)+(i*rm);#V\n",
      "pa1=(Va*i)/1000;#kW\n",
      "pb1=(Vb*i)/1000;#kW\n",
      "print \"power input to motor A is,(kW)= %.1f\"%pa1\n",
      "print \"power input to motor B is,(kW)= %.1f\"%pb1\n",
      "fta1=1960;#kg\n",
      "ftb1=2060;#kg\n",
      "print \"tractive effor of motor A is,(kg)= %.f\"%fta1\n",
      "print \"tractive effor of motor B is,(kg)= %.f\"%ftb1"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "(i) When motors are connected in parallel and train speed is 40kmph\n",
        "power input to motor A is,(kW)= 210\n",
        "power input to motor B is,(kW)= 183\n",
        "tractive effor of motor A is,(kg)= 1625\n",
        "tractive effor of motor B is,(kg)= 1480\n",
        "(ii) When motors are connected in series and current is 400A\n",
        "power input to motor A is,(kW)= 117.4\n",
        "power input to motor B is,(kW)= 122.6\n",
        "tractive effor of motor A is,(kg)= 1960\n",
        "tractive effor of motor B is,(kg)= 2060\n"
       ]
      }
     ],
     "prompt_number": 99
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 8.10, Page 254"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "#Variable declaration\n",
      "f=50;#hz\n",
      "t=0.5;#in meter\n",
      "s=0.25;#\n",
      "\n",
      "#Calculations\n",
      "vs=2*f*t*(3600./1000);#kmph\n",
      "vc=vs*(1-s);#kmph\n",
      "\n",
      "#Results\n",
      "print \"linear synchronous velocity in kmph is %.f\"%vs\n",
      "print \"vehicle speed in kmph is %.f\"%vc"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "linear synchronous velocity in kmph is 180\n",
        "vehicle speed in kmph is 135\n"
       ]
      }
     ],
     "prompt_number": 100
    }
   ],
   "metadata": {}
  }
 ]
}