summaryrefslogtreecommitdiff
path: root/Chemical_Engineering_Thermodynamics___by_S._Sundaram/ch9_1.ipynb
blob: cea268a9dfe9535d226f91efd8eabf57284e9128 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
{
 "metadata": {
  "name": "",
  "signature": "sha256:afa9dce93689170d63e47c0994dac86a6544b2fa2de4c3b3095b5b07cd969902"
 },
 "nbformat": 3,
 "nbformat_minor": 0,
 "worksheets": [
  {
   "cells": [
    {
     "cell_type": "heading",
     "level": 1,
     "metadata": {},
     "source": [
      "Chapter 9 : Fluid Flow in Pipes and Nozzles"
     ]
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 9.1  Page No : 147"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "import math \n",
      "\n",
      "#Given\n",
      "R = 848.0;#gas consmath.tant in m Kgf/Kgmole K\n",
      "M = 29.0;#molecular weight of air\n",
      "g = 9.81;\n",
      "T1 = 90+273.0;#initial temperature in K\n",
      "y = 1.4;#gamma = Cp/Cv\n",
      "W = 800/3600.0;#Mass rate of air in Kg/sec\n",
      "P1 = 3.5;#initial pressure in atm\n",
      "d = 2.5;#diameter of the pipe in cm\n",
      "\n",
      "#To find out the pressure at the final point\n",
      "v1 = (R*T1)/(M*P1*1.033*10**4);#specific volume in cubic meter/Kg\n",
      "u1 = (W*v1)/(math.pi*(d**2*(10**-4))/4);#inital velocity in m/sec\n",
      "#Assume final temperature as\n",
      "T2 = [300,310];\n",
      "#Assume specific heat capacity in J/KgK corresponding to the above temperature as\n",
      "Cp = [2987.56,2983.56];\n",
      "us = []\n",
      "u2 = []\n",
      "for i in range(0,2):\n",
      "    us.append((g*y*R*T2[i]/M)**(1/2));#sonic velocity attained in m/sec\n",
      "    u2.append(((u1**2)-((2*g*Cp[i]/M)*(T2[i]-T1)))**(1/2));#From equation 9.18 & 9.19 (page no 170)\n",
      "if us[i]-u2[i] <= 1: \n",
      "    u2 = u2[i];\n",
      "    Tnew2 = T2[i];\n",
      "v2 = u2*(math.pi/4)*(d**2/10**4)*(1/W);\n",
      "P2 = (P1*v1*Tnew2)/(T1*v2);\n",
      "print \"The pressure at the final point is %f atm\"%(P2);\n",
      "#end\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "The pressure at the final point is 397.261807 atm\n"
       ]
      }
     ],
     "prompt_number": 1
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 9.2  Page No : 151"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "import math \n",
      "\n",
      "#Given\n",
      "A1 = 0.002 #inlet area in sq meter\n",
      "A2 = 0.0005 #throat area in sq meter\n",
      "P1 = 1.3*10**4 #inlet pressure in Kgf/sq m\n",
      "P2 = 0.7*10**4 #throat pressure in Kgf/sq m\n",
      "g = 9.81\n",
      "v = 12*10**(-4) #specific volume in cubic m /Kg\n",
      "\n",
      "#To find out the mass rate of alcohol\n",
      "u2 = ((v*2*g*(P1-P2))/(1-((A2/A1)**2)))**(0.5) #throat velocity in m/sec\n",
      "W = (u2*A2)/v\n",
      "print \"The mass rate of alcohol is \",\n",
      "print \"%.4f\" %W,\n",
      "print \" Kg/sec\"\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "The mass rate of alcohol is  5.1147  Kg/sec\n"
       ]
      }
     ],
     "prompt_number": 2
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 9.3  Page No : 153"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "%matplotlib inline\n",
      "import math\n",
      "import matplotlib.pyplot as plt\n",
      "import numpy\n",
      "\n",
      "#Given\n",
      "P1 = 50.0;#initial pressure in Kgf/sq m\n",
      "T1 = 45+273.0;#initial temperature in K\n",
      "g = 9.81;\n",
      "y = 1.35;#gamma\n",
      "R = 848.0;#gas consmath.tant in m Kgf/Kgmole K\n",
      "M = 29.0;#molecular weight of air\n",
      "d = 1.0;#pipe diameter in cm\n",
      "\n",
      "#(i)To plot velocity,specific volume,mass velocity against P2/P1\n",
      "#(ii)To calculate the critical pressure,critical mass velocity and mass rate of flow\n",
      "#(i)Plotting of graph\n",
      "V1 = (R*T1)/(M*P1*1.033*10**4);#initial volume of the gas in cubic m/Kg\n",
      "#P3 = P2/P1 (say)\n",
      "#Assume P3 values as\n",
      "P3 = [1.0,0.8,0.6,0.4,0.2,0.1,0.05];\n",
      "G = [0,0,0,0,0,0,0];\n",
      "u2 = []\n",
      "v2 = []\n",
      "G = []\n",
      "for i in range(0,7):\n",
      "    u2.append((((2*g*y*R*T1)/((y-1)*M))*(1-(P3[i]**((y-1)/y))))**(1/2.0));#final velocity in m/sec\n",
      "for i in range(0,7):\n",
      "     v2.append(V1/(P3[i]**(1/y)));#final specific volume in cubic meter/Kg\n",
      "for i in range(0,7):\n",
      "    G.append(u2[i]/v2[i]);#Mass velocity in Kg/sq m sec\n",
      "\n",
      "plt.plot(P3,u2,\"o-\")\n",
      "plt.title(\"Velocity vs P2/P1\")\n",
      "plt.xlabel(\"P2/P1\")\n",
      "plt.ylabel(\"Velocity\")\n",
      "plt.show()\n",
      "\n",
      "plt.plot(P3,G,\"+-\")\n",
      "plt.title(\"Mass velocity vs P2/P1\")\n",
      "plt.xlabel(\"P2/P1\")\n",
      "plt.ylabel(\"Mass Velocity\")\n",
      "plt.show()\n",
      "\n",
      "P_3 = [1.0,0.8,0.6,0.4,0.2,0.1,0];\n",
      "plt.plot(P_3,v2,\"*-\")\n",
      "plt.title(\"Sp. volume vs P2/P1\")\n",
      "plt.xlabel(\"P2/P1\")\n",
      "plt.ylabel(\"Specific Volume\")\n",
      "plt.show()\n",
      "\n",
      "#(ii)Calculation of critical pressure,critical mass velocity and mass rate of flow\n",
      "#From equation 9.37(page no 181)\n",
      "P2 = P1*(2/(y+1))**(y/(y-1));\n",
      "print \"The critical pressure is %f atm\"%(P2);\n",
      "#From equation a (page no 183)\n",
      "u2 = (((2*g*y*R*T1)/((y-1)*M))*(1-((P2/P1)**((y-1)/y))))**(1/2.0);\n",
      "print \" The critical velocity is %f m/sec\"%(u2);\n",
      "#From equation b (page no 183)\n",
      "v2 = ((R*T1)/(M*P1*1.033*10**4))/((P2/P1)**(1.0/y));\n",
      "print \" The critical specific volume is %f cubic meter/Kg\"%(v2);\n",
      "#From relation c (page no 183)\n",
      "Gnew = float(u2)/v2;\n",
      "print \" The critical mass velocity is %f Kg/sq meter sec\"%(Gnew);\n",
      "W = Gnew*(math.pi/4.0)*(d/(100.0))**2;\n",
      "print \" Mass rate of flow through nozzle is %f Kg/sec\"%(W);\n",
      "#end\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "metadata": {},
       "output_type": "display_data",
       "png": "iVBORw0KGgoAAAANSUhEUgAAAYkAAAEZCAYAAABiu9n+AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xl8E3X+x/FXaaEFUREPThUohyKi4IWKEFxOUUF0EQQv\nFHRViq4iArqUdRU8FrVV1BUPWK5lxQOt0AKSX3EVPBZXFBGoouIK6gpyCAhtfn98piQtbUnbJJPj\n/Xw88shkMpl8Osp88r1BRERERERERERERERERERERERERERERCKoGVAE1KjmeXY45xIRkSizCJhY\nxv5+wPdUnACaEZokEegl4P4Qnq8qMoF9WPLaCvwL6OS81xd4x9n/PfAcULfU52sBPwKHAV5gt3Ou\nH4H5QEPnuG7AMmAb8FU4/hCJbaH8hyVSVS8BQ8vYfzUwE0sCicYHzAEOB47FksIrzntHAH8GGgEn\nA02AR0p9vguwCtjlnOtW51ytgXrAY85xO4FpwOgw/R0iItVWG/sle0HAvqOwX7+nAknAPcAG4Cfg\nH877cHBJojGwAPgfsB64MeCcNYBxznm2Ax9iN1icc6QDI4DfgL3YL+8FwF3Ay6VizgIeL+NvGQP8\ns9S+J5wHwHVAgfP9XwJXlXEOsJLE3wNen+LEWL+MYy8DPim1bwpwu7O9DBgW8N6twOpSx3dHJQkR\niWJ/w6pNit0E/NvZHgW8iyWAmsAzwGznvWaUTBL5wJNYdctpwA9YlQrYr+VPgFbO6/b4b7pFQAtn\n+0Xsl3qxhtgv7iOd1ynAFqBDGX/HCdiv9+Lqn2Tgv8DZWNXPLwHf3wBoW8Y5oGSSSMVKChvLOfZx\n/Nej2OcB37MMuMHZPgZ4G5he6nglCRGJaudjdey1nNf/wpID2A3vwoBjG2G/9mtQMkkcD+zHbsbF\nHsRu+gBfAJeU8/2lk0TpNomF+EslFwOfVvC3LMeqygB6YCUXnLi2AgOw0lNFMrHSzFYsIS2h7KTU\nA/gZaBmwLx0rRRXzYolrK7AJSz5HlzqPkoSUSW0SEi3+hVUlXYbd5M7C/+v4ROBV7Ca3FViDJYMG\npc7RGLth7grY942zH6ApVtVTFdPxt5sMpWRVUGmzgcHO9lXALGd7F3AlcDNWungTaFPBeYqr1Rpg\nN/FVpd7v5Jz7cvyJCOAi4K2A1z5gpHOuplgC+18F3ytygJKERJMZwDXYTXgR1hMH7EbfG7vJFT/q\nYD17Av0Xqz4K7OlzAvCds/0tJX9xl8dXxr7XseqpdljvolllHFPsZcCDtXf0p2RVUB7QE6vCWkvJ\nKrbSMSRV8B0dnJiuw6qTApVOEiJVpiQh0WQGVn1yIyXrzJ/Bqo1OcF4fC1xaxue/xdouJmH1+O2x\nBtuZzvvTsGqkltgNOLBNItAW/FVPxXZjXUdnAyuxapvy/IhV8byENU5/4ew/DuvWexjWvXUXUFjO\nOSpKEO2wJHobByeDOlgprHTiKO98SUAa1taThF23WuUcKyLiumVYVUjNgH1JwB3YL+/tWNXKX5z3\nmmE32uIfPE2AN5xzbMB6KxWrAYzHbtzbsZt9cVVUIf7E0BKr2tmKv9spQGes7eLaIP6Ooc6xdwbs\na4glj23Oud8GTirn8xOwpFmWF7Dqth0Bj+LeShdjPbICle7dFMjjxFmEXYMiJy6RiGiD/WMrfvwC\nZGC/3hYD67Did72Az4zFGt3WYsVykWhxPCV7LkWjp7A2D5GYUwOrQz4eeBi429k/BpjsbLcFPsZ+\nRTbDfgmqSkyiQQ2sq+k0twM5hOEc3KAvEhN6Yl0DwUoJxf8jFzfggZUixgR8ZhH+qQhE3HIYNk5i\nNf7BdyIJISWC3zUIm2YALEFscba34E8YjYEVAZ/ZhP5RivuivYpJJGwiVZVTCxvEVHq6ArCufmV1\nOQx8X0REXBCpkkQf4CP8/d63YNVMm7HRsz84+7/D2iyKNcXfxx2A9PR0X0FBVcdDiYgkrAKCGydU\nQqRKEoPxVzWBddEr7kZ4LfBawP5BWMmjOTb3zPuBJyooKMDn8+nh8zFhwgTXY4iWh66FroWuRcUP\nbCaDSotESeIwbEqB4QH7JgPzsEnHNgIDnf1rnP3F0y7cgqqbRERcE4kksQubeTLQz1jiKMuDzkNE\nRFymMQgxzOPxuB1C1NC18NO18NO1qL6K5oeJVj6nfk1ERIKUlJQEVbjnqyQhIiLlUpIQEZFyKUmI\niEi5lCRERKRcShIiIlIuJQkRESmXkoSIiJRLSUJERMqlJCEiIuVSkhARkXIpSYiISLmUJEREpFxK\nEiIiUi4lCRERKZeShIiIlEtJQkREyqUkISIi5YrEGtdRIycnn6ysPPbuTSE1dT8ZGT3p27eL22GJ\niESthEkSOTn5jBqVS0HBAwf2FRSMB1CiEBEpR8JUN2Vl5ZVIEAAFBQ+Qnb3YpYhERKJfwiSJvXvL\nLjTt2ZMc4UhERGJHJJJEPeBl4HNgDXAOUB9YDKwD8pxjio0F1gNrgZ6hCiI1dX+Z+9PSCkP1FSIi\ncScSSeIJ4C3gZKA9dvO/B0sSrYGlzmuAtsCVznNvYGqoYszI6El6+vgS+9LSxjF8eI9QnF5EJC4l\nhfn8RwKrgBal9q8FugJbgIaAFzgJK0UUAQ85xy0CMoEVAZ/1+Xy+KgWTk5NPdvZi9uxJJjW1kN9+\n60GtWl14/XVIS6vSKUVEYkJSUhJU4Z4f7iRxOvAsVs10GvARcDuwCTgqIIafndfZWEKY5bw3DVgI\nzA84Z5WTRGn798NVV8Hu3TB/PtSqFZLTiohEnaomiXB3gU0BOgK3AR8Aj+OvWirmcx7lOei9zMzM\nA9sejwePx1O14FJg1iz4/e9h8GCYOxdq1qzSqUREoorX68Xr9Vb7POEuSTQE3gOaO687Y1VKLYBu\nwGagEbAMq24qTiCTnedFwARgZcA5Q1aSKLZ3LwwYAEccATNnQrI6PIlInKlqSSLcDdebgW+xBmqA\n7sBnwBvAtc6+a4HXnO0FwCCgFpZYWgHvhzlGUlOtuumnn2DYMCgqCvc3iojEhnCXJMDaIqZhN/4C\n4HogGZgHnABsBAYC25zjxwHDgP3AKCC31PlCXpIo9uuv0KcPtGkDzzwDNRJmFImIxLtobbgOh7Al\nCYAdO6BXL+jYEbKzISkWr5CISCnRWt0Ucw4/HBYuhJUr4a67IIz5SEQk6ilJlOHIIyE3F95+G+69\nV4lCRBJXwswCW1n168PixeDx2EC7++5zOyIRkchTkqjAMcfA0qXQtav1gLr7brcjEhGJLCWJQ2jQ\nwJ8oatWCVq20cJGIJA4liSA0aWLtE2eemY/Pl8tPP2nhIhFJDGq4DtIJJ0CbNnklEgRo4SIRiW9K\nEpWQnKyFi0QksShJVIIWLhKRRKMkUQllLVxUo8Y4GjbsQaHyhIjEoVicdCKs03IcSuDCRWlphVx5\nZQ9mzOjCvn0wYwa0KL28kohIFNDcTS4qKoInnoAHH7THjTdqzicRiS5KElHgs8/g6quhcWOYNg0a\nNnQ7IhERown+osApp8CKFdChA5x+uq1RISISy1SSCJMVK+Caa6BTJ8jKgnr13I5IRBKZShJRplMn\nWLXKph4/7TSb2kNEJNaoJBEBublwww1wxRUwaRLUru12RCKSaFSSiGK9esEnn8Dmzbbi3Ycfuh2R\niEhwVJKIsLlzISMDbrsNxo6FmjXdjkhEEoG6wMaQ776DYcNg61b4+9+hTRu3IxKReKfqphjSpAks\nWgTXXQfnnw9PPmkD8kREoo1KEi5bt866yh5+OLz4IjRt6nZEIhKPVJKIUa1bwzvv2Mp3HTvCrFkQ\nRzlQRGJcJEoSG4HtQCGwDzgbqA/8AzjReX8gsM05fiwwzDk+A8grdb64KkkE+ugjm9ajXTt4+mk4\n+mi3IxKReBHNJQkf4AE6YAkC4B5gMdAaWOq8BmgLXOk89wamRijGqHDGGZYomjaF9u1h4UK3IxKR\nRBepG3Dp7HUpMN3Zng70d7b7AXOwEsdGYAP+xJIQateGKVNg5ky4+WZ77NzpdlQikqgiVZJYAnwI\nDHf2NQC2ONtbnNcAjYFNAZ/dBDSJQIxRp1s3G4C3Z49NFvjuu25HJCKJqOxFm0PrfOB74Fisimlt\nqfd9zqM8B72XmZl5YNvj8eDxeKobY1Q68kh46SV49VUYMMDGVmRmQq1abkcmItHO6/Xi9XqrfZ5I\nd4GdAOzEShQeYDPQCFgGnIS/bWKy87zI+czKgHPEbcN1RbZsgeHD4ZtvYPjwfBYsyGPv3hRSU/eT\nkdGTvn27uB2iiESxqjZch7skUQdIBnYAhwE9gYnAAuBa4CHn+TXn+AXAbGAKVs3UCng/zDHGhAYN\n4PXX4fbb88nIyKWo6IED7xUU2LrbShQiEmrhbpNoACwHPsZKA29iXVonAz2AdcCF+EsOa4B5zvNC\n4BYqropKKElJsHZtXokEAVBQ8ADZ2YtdikpE4lm4SxJfAaeXsf9noHs5n3nQeUgZ9u4t+z/Zpk3J\n+HxaW1tEQithxiDEi9TU/WXu//rrQs44w6qkErDJRkTCREkixmRk9CQ9fXyJfenp45g9uwf33We9\nnzp2tB5RmjRQRKorFisnErJ3U6CcnHyysxezZ08yaWmFjBzZ40Cjtc8HCxbAxIlQWAh/+hNcdhnU\n0M8BkYSm9SSkBJ8P3nzTksXevZYsLr9cyUIkUSlJSJl8PnjrLUsWv/4K991na20nJ7sdmYhEkpKE\nVMjns4WOJk6E7dstWQwcqGQhkiiUJCQoPh/k5Vmy2LoV7r0XBg1SshCJd0oSUik+HyxZYsnixx8t\nWQweDCmRmM1LRCJOSUKqxOeDt9+2ZLF5M4wfD0OGKFmIxBslCakWnw+8Xhtn8d13liyGDoWaNd2O\nTERCQUlCQsbrtZLFN99Ysrj6aiULkVinJCEhl59vyeKrr2DcOLjmGq1lIRKrlCQkbN55x5LF+vWW\nLK67TslCJNZUNUlo/K0cUufOsHgxzJ4N8+dDq1bwzDM2kltE4puShATtvPMgNxfmzrXZZlu1gqlT\nlSxE4pmShFTauefCwoXwz39CTg60bAlPPgl79rgdmYiEmpKEVNk551iSmD/fpvxo2RKysmD3brcj\nE5FQUZKQajv7bJtx9rXXbBR3ejo8/riShUg8UJKQkDnzTFvL4s03baxFixYwZYrNPisisUlJQkKu\nY0crVSxcaN1n09Phr3+FXbvcjkxEKktJQsLm9NPhlVesR9R771myeOQRJQuRWKIkIWHXvj28/LKN\ntfjgA6uGeugh2LnT7chE5FCUJCRiTj0V5s2zWWdXrbJkMWkS7NjhdmQiUp5IJIlkYBXwhvO6PrAY\nWAfkAfUCjh0LrAfWAj0jEJu44JRTbECe1wurV1s11AMP2Ip5IhJdIpEkRgFrgOIJl+7BkkRrYKnz\nGqAtcKXz3BuYGqH4xCVt29pUH/n58Pnnlizuvx9++cXtyESkWLhvwk2Bi4Bp+CeWuhSY7mxPB/o7\n2/2AOcA+YCOwATg7zPFJFDjpJJg503pCrV9vg/ImToRt29yOTESCSRIfAbcCR1Xh/I8Bo4GigH0N\ngC3O9hbnNUBjYFPAcZuAJlX4TolRbdrAjBnw7rs2PXnLlrYI0tatbkcmkriCWaRyEHA98AHwIfAi\n1pZwqPm6LwZ+wNojPOUc4zvEecp8LzMz88C2x+PB4ynv9BKLWrWCl16CDRusraJVK7jlFrj9dqhf\n3+3oRGKD1+vF6/VW+zyVmVu8BnbjfxorGbwAPAH8XM7xDwJXA/uBNOAI4BXgLCxpbAYaAcuAk/C3\nTUx2nhcBE4CVpc6r9SQSTEEBPPigDdD7wx/gjjvg6KPdjkoktoR70aHTsNJEHyAXmA10BoYCpwfx\n+a7AXcAlwMPA/4CHsMRQz3lu65z3bKyaaQnQkoNLE0oSCeqrryxZvPIK3HQT/PGPsHJlPllZeezd\nm0Jq6n4yMnrSt28Xt0MViTpVTRLBVDd9BPyCNT6PAYpXD1gBnF+J7yq+s08G5gE3YA3UA539a5z9\na7DSxy0cukpLEkjz5vDcc7bu9qRJ0Lx5PikpuWzb9sCBYwoKxgMoUYiESDBZpQXwZal9zYGvQh9O\nUFSSEAC6dLmX5cv/ctD+Xr3uY9Gi+12ISCR6hXP50peD3CcSUTVqlF0Q/v77ZPQ7QiQ0KqpuOhlr\nJ6gHDMAykA9rgE4Lf2giFUtN3V/m/oKCQtq1g4wMGDoUDjsswoGJxJGKShJtsIbmI53ni53njsDw\n8IcmUrGMjJ6kp48vsS89fRxz5/YgO9umKj/xRBg9GjZudCdGkVgXTP3UucB74Q6kEtQmIQfk5OST\nnb2YPXuSSUsrZOTIHiUarb/8EqZOhRdfhC5drHTh8UBSpWtmRWJbOLrAjsG6qWaX8Z4PyKjsl4WI\nkoRU2s6dNvVHVhYkJ1uyGDIE6tRxOzKRyAhHkrgEm7n1Okp2RS1um5hexmciQUlCqszng6VL4Ykn\nYMUKGDYMbr0VTjjB7chEwivcg+miiZKEhMSGDfDUUzZfVLduVrq44AJVRUl8CmcX2MWUXPOhPjbq\nWiSmtWwJjz1mjdrdusGIEdChA7zwAuze7XZ0ItEhmKzyMQdPvVHWvkhRSULCoqjIlljNyrJlVm+8\n0eaKOv54tyMTqb5wliQKgRMDXjej5NTfInGhRg3o1Qtycmxti1274LTTYOBAe63fJpKIgskqvYG/\nAfnO6y7ACGyWVjeoJCERs307TJ8O2dlQt661WwwaBGkaTioxJtwN18cC5zjbK4CfKvtFIaQkIRFX\nVAS5uVYV9e9/w/DhVhXVRMtiSYwIZ3UTwHlAN+dxbmW/RCTW1agBffrYKO78fFtatV07K1W8956q\noiR+BZNVJmMLBc1yjh+ErVA3NoxxVUQlCYkKv/xiI7mzs23FvIwMa79ITXU7MpGDhbO6aTXWk6nQ\neZ2M9W46tbJfFiJKEhJVCguthJGVBZ98Ygsi3XwzNGrkdmQifuGsbvJRcpxEPbQYkMgByclw8cWQ\nlwdvvw0//ght29q0HytLL74rEmOCySqDsSonr/O6K7bc6NwwxXQoKklI1Nu61aqinnwSjj3WqqJ+\n/3uoVcvtyCRRhbt3U2OsXcIHvA9sruwXhZCShMSMwkIbd5GVBZ99ZtVQN98MDRq4HZkkmnAkiTM4\neGI/Avb9u7JfFiJKEhKTPv3UGrnnzYNLLoGRI+Gss9yOShJFOJKEl4rbHrpV9stCRElCYtrPP8Pz\nz9vkgo0bW1XU5ZdDzZpuRybxTLPAisSY/fvhjTesKmrdOhucN2IEHHec25FJPApn76bDgPuA55zX\nrbClTEWkGlJS4LLLYNky60L79dfQpg1cd52N6haJBsEkiReB37BR1wD/BR4IW0QiCah9e3juOVvj\n4uSToX9/6NzZ2i/27XM7OklkwSSJdGwZ09+c17uCPHcasBIbeLcGmOTsr4+tUbEOyKPkGIyxwHpg\nLdAzyO8RiRtHHw1jxtja3HfcYV1oW7SASZPgJzdnTJOEFUyS2AsErgSc7uw7lD1Y4/bpQHtnuzM2\nxmIx0BpY6rwGaAtc6Tz3BqYGGZ9I3ElJscbs/HxYsADWr4dWreCGG+Djj92OThJJRTfhqdhNPRNY\nCDQFZgNvA2OCPP+vznMtbDqPrcCl+NfHng70d7b7AXOAfcBGYANwdpDfIxK3ilfLW7cO0tNtdHfX\nrjB/vjV+i4RTRS3dt2O/7Btj1ULfYmMjKjNVeA3nM+nA08DdWKI4KuD7f3ZeZzvnnuW8Nw1LTvNL\nnVO9mySh7dsHr75qvaK++QZuvdVW0VuxIp+srDz27k0hNXU/GRk96du3i9vhSpSoau+mlAree9x5\nNMNmfh0EDMFKE3OwNoVDKcKqm47E1sUuPbbCR8VjMcp8LzMz88C2x+PB4/EEEYpIfKhZ02abHTgQ\nPvrIBuidcEI+ycm57Njh71NSUDAeQIkiQXm9Xrxeb7XPU9ms0gHr7XQqVn1UGfcBu4EbAQ82tUcj\nYBlwEv62icnO8yJgAtb4HUglCZFSunW7F6/3Lwft79XrPhYtut+FiCTahHOcRArWjjAbu3GvBQYE\n8blj8Pdcqg30AFYBC4Brnf3XAq852wuw0kotoDk2HuP9IL5HJOH5fGVXCnz9dbK60Eq1VFTd1BO7\naffFbtZzsLWtdwZ57kZYw3QN5/F3rDfTKmAecAPWQD3QOX6Ns38NsB+4BU1JLhKU1NSyW7C3bCmk\nZUvrTnvjjbZOt0hlVFT0eBtLDPOxxuVooeomkVJycvIZNSqXggJ/m0R6+jieeKI3xx3XhUcesbUu\nbrrJ5orSLLSJR3M3iSS4nJx8srMXs2dPMmlphYwc2aNEo/WGDTBlCsyZY43ed94JrVu7GLBElJKE\niATlhx9sBtqpU+GCC+Duu6FTJ7ejknBTkhCRStm1ywbpTZkCxx8Po0dD375QQ/McxCUlCRGpkv37\n4eWX4eGHYc8eSxZXXQWpqW5HJqGkJCEi1eLzwdKl8MgjtoreqFHW0H3kkW5HJqEQznESIpIAkpKg\ne3fIzbV1uf/zH5uBdvRo+O47t6MTtyhJiMhBTj8dZs2yxY/27YNTT4Xrr4fPPnM7Mok0JQkRKdeJ\nJ8Ljj1v32ZYt4Xe/g0susSnMVeubGNQmISJB270bZsyARx+1BZLuvhv69YPkys7kJhGnhmsRiZjC\nQnj9dXjoIfj5Z7jrLrjmGqhd2+3IpDxKEiIScT4fLF9uPaI++ABGjoQ//AHq13c7MilNvZtEJOKS\nkqBLF3jjDes+W9x2cfvt8PXXbkcnoaAkISIhccop8OKLsHo11KoFHTvCkCFakzvWKUmISEg1aWKj\nt7/8Ek47zab66NXLShqqKY49apMQkbDauxdmz7Z2i7Q06xF1xRWQUtFqNhJyargWkahWVARvvWWl\njG+/hT/+EYYNg8MOczuyxKAkISIx4733rGTxzjvWG+q22+DYY92OKr6pd5OIxIxzz4VXXrEksXkz\ntGkDt9wCBQVuRyalKUmIiGtat4Znn4XPP7exFeecY6vmffCB25FJMVU3iUjU2LEDnn/eFkJKT7dG\n7t69bTyGVI/aJEQkbuzbB//4h7VbFBXZdOWDBtn4C6kaJQkRiTs+H+TlWY+odevgjjtg+HA4/HC3\nI4s9argWkbiTlOQfiPfaa/D++9C8OYwdC99/73Z0iSHcSeJ4YBnwGfApkOHsrw8sBtYBeUC9gM+M\nBdYDa4GeYY5PRGLEGWfA3LmWKHbuhLZtrVTxxRduRxbfwl3d1NB5fAzUBT4C+gPXAz8BDwNjgKOA\ne4C2wGzgLKAJsARoDRQFnFPVTSLCTz/BU0/Z47zzrJH7vPPsvZycfLKy8ti7N4XU1P1kZPSkb98u\n7gbsslhpk3gNeNJ5dAW2YEnEC5yElSKKgIec4xcBmcCKgHMoSYjIAb/+Ci+9BH/9KzRqBN265TNn\nTi4FBQ8cOCY9fTxPPNEroRNFLLRJNAM6ACuBBliCwHlu4Gw3BjYFfGYTVqIQESlTnTo2EO+LLyAj\nA7Ky8kokCICCggfIzl7sUoSxLVJTbNUF5gOjgB2l3vM5j/Ic9F5mZuaBbY/Hg8fjqXaAIhLbUlJs\nIN5TT6WQn3/w+3v2JNYaq16vF6/XW+3zRCJJ1MQSxN+x6ibwVzNtBhoBPzj7v8Mau4s1dfaVEJgk\nREQCpaXtL2d/YYQjcVfpH9ATJ06s0nnCXd2UBDwPrAEeD9i/ALjW2b4Wf/JYAAwCagHNgVbA+2GO\nUUTiSEZGT9LTx5fYV7PmOAYP7uFSRLEt3A3XnYF84BP81UZjsRv/POAEYCMwENjmvD8OGAbsx6qn\nckudUw3XIlKhnJx8srMXs2dPMmlphTRs2IPly7uwZImNs0hEsdK7KRSUJESk0qZOhUmTbAT3ySe7\nHU3kVTVJaG0oEUkIt9xi03lceKEtftShg9sRxQYlCRFJGFdfbSvh9e4Nr77qH3wn5dPcTSKSUAYM\ngBkzoF8/WLLE7Wiin5KEiCScXr1sZbyrroLXX3c7muim6iYRSUgXXGBtExdfbFN7DB7sdkTRSUlC\nRBLWmWdalVOvXrYq3ogRbkcUfZQkRCShtWsH//d/0KOHJYo773Q7ouiiJCEiCa9lS8jPh+7dLVFM\nmKB1tYvF4mXQYDoRCYstW6BnT0sWjz4aX4lCI65FREJg61bo0wfat4enn4bkOJk8VklCRCREduyw\ncRQNG8L06VCzptsRVV8sLDokIhITDj8ccnJg+3a44grYs8ftiNyjJCEiUobatW3AXVqajaXYudPt\niNyhJCEiUo5atWD2bDjxRGvQ3rbt0J+JN0oSIiIVSE6G556Ds8+Gbt3gxx/djiiylCRERA6hRg14\n7DGrdurSBb47aFHl+KXBdCIiQUhKgvvvt0btCy6w6TxatHA7qvBTkhARqYS774YjjoCuXSE3F9q2\ndTui8FKSEBGppJtvhrp14Xe/s66yHTu6HVH4KEmIiFTB0KElV7k7/3y3IwoPNVyLiFTRZZfBzJnQ\nvz8sXux2NOGhJCEiUg09e1pJYsiQ+FzlTtVNIiLV1LkzLFzoH5k9ZIjbEYVOuEsSLwBbgNUB++oD\ni4F1QB5QL+C9scB6YC3QM8yxiYiEzBlnWLfYMWPg2WfdjiZ0wp0kXgR6l9p3D5YkWgNLndcAbYEr\nnefewNQIxCciEjKnnGKr3E2ebOtRxINw34SXA1tL7bsUmO5sTwf6O9v9gDnAPmAjsAE4O8zxiYiE\nVHo6LF8O06bZCnexvrKBG7/UG2BVUDjPDZztxsCmgOM2AU0iGJeISEg0bWrLob7+uq2ZHcuJwu2G\na5/zqOj9g2RmZh7Y9ng8eDyekAYlIlJdxx0Hy5bBRRfBiBHwzDORXeXO6/Xi9XqrfZ5IrEzXDHgD\nONV5vRbwAJuBRsAy4CT8bROTnedFwARgZanzaWU6EYkZO3faKnfHHQczZri3yl0srUy3ALjW2b4W\neC1g/yCgFtAcaAW8H/HoRERCqG5dm7pj5064/PLYW+Uu3EliDvAu0Ab4FrgeKyn0wLrAXoi/5LAG\nmOc8LwTGAFkYAAAG+ElEQVRuoeKqKBGRmJCWZqvc1akDffvG1ip3kahuCjVVN4lITCoshJtugs8+\ng7fegqOOitx3x1J1k4hIQipe5a5TJ1vl7ocf3I7o0JQkREQiKCkJpkyxxuwuXWDTpkN/xk1ud4EV\nEUk4SUkwcWLJVe7S092OqmxKEiIiLrnrLksUXbtCXl50rnKnJCEi4qKbbvKvcvfmmzZRYDRRkhAR\ncdmQIbbKXZ8+1lW2c2e3I/JTw7WISBTo3x9mzYIBA6JrlTslCRGRKNGjh5Ukhgyx1e6igaqbRESi\nSOfOsGiRjczetQuGDnU3HiUJEZEo07EjLF1q62fv3Ak33+xeLEoSIiJRqG1bW5Oie3fYsQNGj3Yn\nDiUJEZEo1aKFrXLXvTts3w5//rMNxIskTfAnIhLlfvgBevWyQXePPVa1RFHVCf6UJEREYsC2bbbK\nXdu28OyzlV/lTrPAiojEsXr1bOqOjRuti+xvv0Xme5UkRERiRN26NnXHr7/aoLvdu8P/nUoSIiIx\nJC0N5s+3iQEjscqdkoSISIypWRNmzrTpxXv0gK1bw/ddShIiIjEoORn+9jc47zzweGDLlvB8j5KE\niEiMSkqCRx+Fyy6z7rHffhv679BgOhGRGJaUBJmZ1kbRpYvNINuyZejOryQhIhIH7rzTEoXHYxME\ntmsXmvNGY3VTb2AtsB4Y43IsIiIxY8QIePhhm8bjww9Dc85oSxLJwJNYomgLDAZOdjWiKOb1et0O\nIWroWvjpWvgl4rW46iobkX3RRTbvU3VFW5I4G9gAbAT2AXOBfm4GFM0S8R9AeXQt/HQt/BL1WvTr\nB7Nnw+WXw5//nE+vXvdW+VzRliSaAIHt85ucfSIiUgndu8Po0flMnJhLXt5fqnyeaEsSmrlPRCRE\nlizJo6jogWqdI9pmge0EZGJtEgBjgSLgoYBjNgDpkQ1LRCTWpQMF0XbPr7QUoABoBtQCPkYN1yIi\nEqAP8AVWYhjrciwiIiIiIhJrghlUl+W8/x+gQ4TicsOhrsUQ7Bp8AvwLaB+50CIu2MGWZwH7gQGR\nCMolwVwLD7AK+BTwRiQqdxzqWhwDLMKqsD8FrotYZJH1ArAFWF3BMXFx30zGqpuaATUpu23iIuAt\nZ/scYEWkgouwYK7FucCRznZvEvtaFB/3NvAmcHmkgouwYK5FPeAzoKnz+phIBRdhwVyLTGCSs30M\n8D/ic1qiC7Abf3lJotL3zWjrAlssmEF1lwLTne2V2D+IBhGKL5KCuRbvAb842yvx3xTiTbCDLUcC\nLwM/RiyyyAvmWlwFzMfGGwH8FKngIiyYa/E9cISzfQSWJPZHKL5IWg5UtLpEpe+b0ZokghlUV9Yx\n8XhzrOwAwxvw/1KIN8H+f9EPeNp5Ha9jb4K5Fq2A+sAy4EPg6siEFnHBXIvngFOA/2LVLKMiE1rU\nqfR9M1qLW8H+wy7d5zcebwiV+Zu6AcOA88MUi9uCuRaPA/c4xyYRfWOBQiWYa1ET6Aj8DqiDlThX\nYPXR8SSYazEOq4byYAMGFgOnATvCF1bUqtR9M1qTxHfA8QGvj8dfZC7vmKbOvngTzLUAa6x+DmuT\nCONihq4K5lqcgVU3gNU998GqIBaEPbrICuZafItVMe12HvnYjTHekkQw1+I8oHjocQHwFdAGK2El\nkri5bwYzqC6wAaYT8dtYG8y1OAGrk+0U0cgir7KDLV8kfns3BXMtTgKWYA27dbDGzLaRCzFigrkW\nU4AJznYDLInUj1B8kdaM4BquY/6+WdagupucR7Ennff/gxWr49WhrsU0rCFulfN4P9IBRlAw/18U\ni+ckAcFdi7uwHk6rgYyIRhdZh7oWxwBvYPeK1Vijfjyag7W7/IaVJIeRuPdNERERERERERERERER\nEREREREREREREfErxMacrAbmAbWxkarLsHEHn3LwuINOwN+ArtiEi6uANcCfnPePdj6/A8gOb/gi\nIhJOgXP6zATuABoCpzv76mKDtwJH9k4ELsOSxBvOvjrAOmz65jrYvFo3oSQhMSJaZ4EViSbvAC2B\nzdiUDwA7gc+BxgHHXYhNgxE4gdqvwEfO53/FFoXaG+Z4RUJGSUKkYinYlA+flNrfDCsdrHReH4NN\nJFh6VtGjsWqoTwP2xeNsxRKnonUWWBG31cbaFMBmT30+4L262KJGo7ASBUBPIDfgmAuAfwNF2Ipo\nn4czWJFwUZIQKdtuyl7/tya22ttM4LWA/b2Bvwa8Xg5cErboRCJE1U0iwUvCShRrsMWNAve3x2bV\nDPY8IjFBJQmRspXVbnA+MBRrnyiuihqHraW9KuA4XzmfB1uH+XBs3YN+WDXV2uqHKyIi0Wo8MNDt\nIERERERERERERERERERERERERERERERERESkfP8Pe6pb7AKfc/YAAAAASUVORK5CYII=\n",
       "text": [
        "<matplotlib.figure.Figure at 0x104438150>"
       ]
      },
      {
       "metadata": {},
       "output_type": "display_data",
       "png": "iVBORw0KGgoAAAANSUhEUgAAAZUAAAEZCAYAAABfKbiYAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XeYVOXZx/EvVToIKETaqmBBkaaCUXSxZRUDxkIVwd03\nRonBEg02ZNVoLCHGEokiUqSIQSBoBERxTGwUAQEJoSgiKCAiRZG68/5xP+MMy+zu7O6cOVN+n+ua\na2bOnJm59yjnnqfdB0REREREREREREREREREREREREREREREpJQKgOPK+RnLgXPjEIuISEZbB+wF\nGhTavhg7WTdPdEBlEI+kEikfeCmOn1cWA4GDwC5gB/bfo5t7rTMwB/gW2AK8AjSO8hn/A1oBY7D/\nxrvce94ETnT7nArMBr7BjqOkkYp+ByAZKQh8BvSJ2NYGqO5eE/+8D9QG6gGjsORRz93+DrRwt13A\n6ELvPR47p6zG/js+6j6rKZaIxrj99gEvA3ne/Rkikkk+B+4B5kds+zNwN4e2VLphv5Z3AOuBYRH7\nVwPGA1uB79xnHe1eGwisBXZiyatvlBiOAXYDR0Zsa4/9eq7knucCK4BtwCwObUFFtlTqAuOwE+c6\n97dViNj31+5zdgKfAu3c9nXABUAO9qt+H3ayXgxcBSwsFPNtwPQof0svYEGhbbcC/3SPL3XfuxPY\nAPw+ymeAHbf/RDyv6f7ODlH27eA+L9Jg4K/u8WjggYjXumF/W6SWqKUiInHwOXYyXQmchJ3Ev8RO\n2pFJ5TzgFPe4DbAJ6OGe/waYgSWXClhCqI2dCHdgXTAAjYDWRcTxNvB/Ec8fB551j3tgv7hPxH59\n34P9ig+JTCrjgGnuu1tgXUC57rWrsRN5R/f8+Ii/73PgfPd4mPuckKpYt9FJEdsWA7+K8ndUx07w\nLSO2LQB6usdfA2e7x3WxYxXNQMJJpTJwM3Ysa0fZ9xbgg0LbZgEXucejgQfd41rARODdQvsrqYhI\nXISSyj3Aw9gv9dlYciluTOWvwF/c4+uwk3ybQvvUxFouV2An2+LkYYkFLDGtB85xz2cSTgxgieUH\noJl7HkoqlbBWRuTJ/3rgHfd4NvC7Ir4/Mqnkc/iYygjgj+7xKViLqUoRn/USMNQ9boUlmWru+Rcu\npjpFvDdkILAfO37fYEnj/Cj7nYYlvLMjttXAWo2h+MYAP7rP+hprYR1b6HOUVEQkLkIn0+bYCW8S\n0A/7dRyZVDphJ+ctwHbsJDXWvVYZuA/r1tmI9d9Xdq9djA0Mfwe8TniAuLAjsS6wxliraF3Eayuw\n7prvIm4/YAPWEE4qjdzjyASWA6xyjz/Fup+KOw4QPal0xrrvAB7BkkxRLnExg7V6xka8djp2Ut8G\nBCL+hsIGcmj3VzQtsZZXv0Lbf0m4uw0O7/4q6rOUVNKMBurFT+uxk+YlwNQor0/EToZNCQ8Uh/6f\nPYCdtE4Bfg5cBlzrXnsTSyyNsS62kUV8/3du317YuMukQrFdjyWe0K0m8FGhz9iK/brPitjWHDvx\ngnXrtaRk0SYofISNs5yLTWoobnbYW8BRQFugN3bsQhYCl7vXp2OD72XRApsB9gAwodBrlwJvFNpW\nAck4Siritzzs1/qPUV6rhZ349wFnYif+0Mk3G+v6qoS1KPZj02GPxsZDarptP7jtRZkIDACu5NAT\n8d+xiQOh8Zi62PhIYQexk/RDLt4W2CD5ePf6C8Dt2MB2BSzBROve24QlpsIn4peAZ7BjUHgMI9J+\n4B/YhIcjsZM/WHdUPxd/aLpwccejKE2AuS6W56O8ngP8K+J5SQmlGjZuBHCEu4mIlElkt0+kytgJ\nL3TSvRLrktoJvAY8RXgwuzfWCvkeOyH/FfuR1Bjr4tmOJaS5HDreUVg19/nLorx2DbCU8OyzFyJe\nO0h4oL4edvLf4va7l0NPqr9xse5yn9fWbY88DvWxrqdtHDrrq7n7rsiZb0U5B+tOejpiWxVsfGib\n+zvmYS27aAYA/y7itWHus3dF3EKzv07l8ONXXPdXlvusAuxvKyDczSdSpBeBzRz6P9vjwH+BT7Du\njroRr92FzbZZiXVdhHR0n7EaeDJi+xHAZLf9I+wXoki6Cc3sOt7vQIrxB2zMR8RTXbCpi5FJ5SLC\nXW6PEP4fsTWwBPtVlQWsIfxLbz7W9QHWZ5vjHg8iPP2zF7aYSiTd3IaNlySzqyl6MoRIXGURvVsB\nbL59qN/5LmBIxGuzsBkqP8NaNiG9sb7u0D6d3OPK2BRIkXSyDusia1vCfiJJo3LJu3gml/Bsm2M4\ndFbNBmxgcD/hWTRgU0ebuMdNsJk1YDOBdmD90ts8ilck0bL8DkCktPya/XUPNptlYkk7iohI6vCj\npTIQm9N+QcS2jYRXKoOtS9jgtjeNsj30nubAV9jfUZcorZTjjz8+uHbt2jiFLiKSMdYS2xqrQyS6\npZID3IGtI9gTsX0GNl5SFSvl0AoboN+EzXzphA3c9ye8ancGNgUSrPje20Sxdu1agsGgbsEgw4YN\n8z2GZLnpWOhY6FgUf6OMMw69bKlMwkpfNMTGPoZhA/JVCS/M+hCbxbUCW0C2AhsfGUR4kdsgrI5Q\ndWz21yy3fRS2NmA1Voeot4d/i4iIxMDLpNInyrYXi9n/YXcr7GMOLxoIVsSvZ5TtIiLiE5VpySDZ\n2dl+h5A0dCzCdCzCdCzKLxMKvgVd/6CIiMSoQoUKUIYcoZaKiIjEjZKKiIjEjZKKiIjEjZKKiIjE\njZKKiIjEjZKKSIYLBPyOQNKJkopIhlNSkXjSOhWRNHfwIGzeDBs3wldf2S30eONG+OYbmDoVsrL8\njlSSSVnXqSipiKSoYBC2bz80QUS737IFGjSAY46BJk3sft8++PZbqF0bJk6EatXg1FPhzjvhyiv9\n/sskGZQ1qfh5kS4RKcKPPxafKEK3KlXCiSJ037o1XHhheFvjxrZfUVq1gt/+Fh55BH79a5g3D4YM\nsUQkUlpKKiIJdOCAdUUVThCFk8YPPxyaKEL3HTvafehWq1Z84jrqKBg+HG69Ff74RzjxRPjd7+x5\nnTrx+Q7JDOr+EomDYBC++6741kVo/KJhw0MTRbTk0aABVEjQv85AAArXUVy7FvLz4c034Y47rCVT\nvXpi4pHkoDGVoimpyGGinUiLsnt30S2KyNbGEUcUnyiaNIFGjYrviko2y5fDffdZl9i990JeHlSt\n6ndUkghKKkVTUpHD5OfbSTJyVlRR3VE//lh8ogi9VrOm33+VdxYssOO1erUdu379oFIlv6MSLymp\nFE1JRX6yfz8MGwZPPw179lhXVFGJInRfv37iuqKS3bvvwj33wLZt8OCDcMUVOjbpSkmlaEoqAsCE\nCfCHP0CNGrBmDQwdChUrWjeYrs0Uu2AQZs2y5FKhAjz0EPziF0ou6UZJpWhKKsKECXDLLdaFM3gw\n3H+/deNI2RUU2KLJoUOtxffQQ3DuuX5HJfGidSoiUezaZTOX5s+HOXOgXTu/I0ofFSvCVVfBr34F\n48fDgAE2FfmPf4TTT/c7OvGLan9J2lqwADp0sFlZH398aEJRd1f8VKpkCeV//4Pu3aFHD1uV/+mn\nfkcmflBSkbRTUACPPw7dusHDD8PIkYfPzFJSib+qVWHQIJsh1rkzdO0K114Ln33md2SSSEoqkla+\n/hpycmD6dGupXH213xFlnho1bMHkmjVw3HFwxhlw4402PVvSn5KKpI033rDurrPOsqmvLVr4HVFm\nq1PHJkP8739WuLJNG7j9dti61e/IxEtKKpLy9u61mV033ACTJ9vMrsqagpI0GjaExx6z1fk//miD\n+ffdBzt2+B2ZeEFJRVLaypXQqRN8+SUsWaIprcnsmGPgb3+DhQth/Xqrjvzoo1YGR9KHkoqkpGAQ\nXngBunSx/vopU2zluyS/Y4+FMWOsi3LhQmjZEp55xlqckvq0+FFSzvbtcP311kp5+WW7foikrkWL\nbFHqihVWQqd/f3VfJoOyLn70sqXyIrAZWBaxrT4wB1gFvAnUi3jtLmA1sBK4OGJ7R/cZq4EnI7Yf\nAUx22z8CNCybAd5/39abNGpkCxqVUFJfhw42yWL8eGvBnHoqvPKKTQ2X1ONlUhkN5BTadieWVE4A\n3nbPAVoDvdx9DvAs4Qw5AsgDWrlb6DPzgG/dtieAR734IyQ5HDwIDzxgi+qeftpu1ar5HZXE0znn\n2CUJnnrKBvY7dIB//cu6OiV1eN39lQW8BrRxz1cC52EtmMZAADgJa6UUEE4Ms4B84AtgLnCy294b\nyAZucPsMA+Zh5Wa+Bo6KEoO6v1Lc+vVwzTV2HZJx46x6sKS3YNDWGt17L9Sta4tYtWA1sZKx+yua\nRlhCwd03co+PATZE7LcBaBJl+0a3HXf/pXt8ANiBda9JGnn1VasjdemldhVCJZTMUKGC1RRbutRW\n6eflwUUXWZenJDc/Z38F3U3kMLt3w29+Y6XqX3sN7rxTF4XKRJUqWSt15UorXnnFFXD55bBsWcnv\nFX8keo5FqNtrE/AzYIvbvhFoFrFfU6yFstE9Lrw99J7mwFfY31EX2BbtS/MjapxnZ2eTrXZ0Ulu6\nFHr3hvbtYfFiW5ktma1KFfuRce218Pe/w4UX2u3++21KspRfIBAgEAj4HUaJsjh09tdjwBD3+E7g\nEfe4NbAEqAocC6wl3Jc3D+jknr9BeKB+EDaIDzbW8nIRMQQlNRQUBINPPRUMNmwYDI4da89Fotm5\nMxh88MFgsEGDYPDXvw4G16/3O6L0Qxl7krwcqJ+EDco3xFoo9wH/BF7BWhjrgJ7Adrf/3UAuNj5y\nMzDbbe8IjAGqY0llsNt+BPAS0B6bBdbbfWZh7vhIMtu6FXJz7brwkybZamuRkmzbZhWpn3vOyu/f\ndRccfbTfUaUHXfmxaEoqSe7tt+2E0LevXeCpalW/I5JUs2mTzRCbMMFqwN1+Oxx5pN9RpbZUmf0l\n8pP9++2X5bXXwosv2toEJRQpi8aNbX3LokWWYE44wZLM99/7HVnmUVIRX6xda4vdli61wfiLLy75\nPSIladECRo2C996z/7datYInn4Q9e/yOLHMoqUjCTZhgVwbs2xdef1194BJ/J55odeFmzoS33rKW\nywsvwIEDfkeW/jSmIgmzaxfcdBPMm2eD8e3b+x2RZIoPP4R77oENG2wacq9eUFE/qYulMRVJagsX\nWi2nKlXg44+VUCSxzjoL5s6FESOsO6x9e5gxQ3XFvKCWiniqoACGD7dpn888Az17+h2RZLpg0Ko0\n3Hsv1KgBDz0EF1xgrwUCqjEWoinFRVNS8cnXX9tU4R9+gIkTdc14SS4FBXb56fvug2bNLLnMng0R\nBTgymrq/JKm88YZ1d511ll3hTwlFkk3FitCnj10crG9fG2d5/32/o0p9aqlIXO3dC0OGwNSpdtEl\nXTNekl0gYLctW2zMZehQSzjZ2ZndFabur6IpqSTIypX2y++442DkSF0zXlJP06ZWsPKyy/yOxH/q\n/hLfBIO24KxLFyuRMWWKEoqkpvbtrbqDlF2iS99Lmtm+Ha6/3lopgQCccorfEYmU3Q03QL9+1hWm\nRbllo5aKlNn770O7dtCokV2RTwlFUl23bnYRsJde8juS1KUxFSm1gwdt+uWzz8Lzz0P37n5HJBI/\n//433HgjLF9ulzXOVGUdU1H3l5TKl1/a5V0rVbKV8bpmvKSbLl1g3z4rJ9S5s9/RpB51f0nMpk6F\n00+HnByYM0cJRdJThQp2wTgN2JdNJjTu1P1VTrt3w223WSKZOBE6dfI7IhFvbdwIp55qBShr1vQ7\nGn9oSrF4YulSOOMMqzC8aJESimSGJk3g7LNteryUjpKKRBUMWgHICy6wFfLjx0Pdun5HJZI46gIr\nG3V/yWG2brV/UF99Zdc9adXK74hEEm/fPis0+d57mflvQN1fEhdz59rak5NOgg8+yMx/TCIAVava\nTMfRo/2OJLWopSIA7N9vJcDHjbN/RLpmvAh8+qn9W/jiC6icYQsw1FKRMvvsMzjnHPjkE1i8WAlF\nJOSUU6wLbPZsvyNJHUoqGW7CBJvR1acPvP666h2JFKYB+9JR91eG2rULbrrJVg1PmqRrxosUZedO\naN4cVq3KrB9d6v6SEgUCdr9woV2VsUoVK7WihCJStDp1oEcPm1YvJVNLJYMMGwa1asHjj9salJ49\n/Y5IJDVkYpFJFZSUYm3dar+0Gje2MvVZWX5HJJI6QkUm589XVYmSqPsrzQUC1kL5+c9tltcFF8CY\nMeGuMBEpWYUKcN11GrCPhV8NubuAa4ACYBlwHVATmAy0ANYBPYHtEfvnAgeBwcCbbntHYAxQDXgD\nuDnKd2V899fzz8Nzz8Gll8KDD/odjUhq2rgR2rSxyz9kQpHJVBqozwJ+DXQA2gCVgN7AncAc4ATg\nbfccoDXQy93nAM8S/kNHAHlAK3fLScQfkErWrIG777Yr2VWq5Hc0IqmrSRM46yx49VW/I0lufiSV\nncB+oAY2plMD+AroDox1+4wFLnePewCT3HvWAWuATsDPgNrAfLffuIj3CHDggJWZuO8+aN0asrP9\njkgkteXlwahRfkeR3PxIKtuA4cB6LJlsx1oojYDNbp/N7jnAMcCGiPdvAJpE2b7RbRfnT3+C2rVt\nPQooqYiU12WXwcqVsHq135EkLz9mfx0P3IJ1g+0A/oGNr0QKultc5Ofn//Q4Ozub7Aw4uy5YYNOG\nFy2CipqOIRIXkUUmH37Y72jiKxAIEIjDDB4/Bup7ARcB/+ee9wc6A+cDXYFNWNfWO8BJhMdWHnH3\ns4BhwBdun5Pd9j7AecANhb4v4wbqd++2xY333w+9evkdjUh6Wb4cfvGL9C8ymUoD9SuxJFIdC/hC\nYAXwGjDA7TMAmO4ez8AG8qsCx2ID8vOx5LMTG1+pgCWn0Hsy2pAh0LGjEoqIF049FZo2hTffLHnf\nTORHnv0EG1RfiE0pXgQ8jw26v4LN5lqHTSkGSzivuPsDwCDCXWODsCnF1bEpxbMSEH9Smz0b/vlP\nqzgsIt4IDdhfeqnfkSSfWJo2bbC1JKkqY7q/vv0W2raFsWNtkaOIeGPHDmjRwgbsjzrK72i84WX3\n1whgAdYq0FXKk1QwCDfcAFdfrYQi4rW6da3I5Esv+R1J8oklqZwD9AOaY11VkwBdxinJTJgAK1bY\nNGIR8V5urnWBZUhHSMxK07SpjC0ufAqbClwRuBtI9vWlad/9tX69Dcy/+abK2IskSjAIJ5xgP+jO\nPNPvaOLPy+6vtsATwH+xab+XYdN4u7rt4qOCAhgwAH7/eyUUkUQKFZnUCvtDxZKF3gVGAVOA3YVe\nuxabyZXM0rqlMnw4TJsG776r2l4iiRYqMrlhA9So4Xc08eVlS2UaljgiE0qoGnCyJ5S0tmyZjaGM\nG6eEIuKHUJHJKVP8jiR5xJJUBkTZdl28A5HS2bvXykU89hgcd5zf0YhkrtxcXWclUnFNmz5AX6AL\n8J+I7bWx65qkysTVtOz+GjIEVq2CqVMz5/KmIslo3z5bYf/BB9Cypd/RxE9Zu7+Ke0MLrCzKI8CQ\niH13YaviD5T2y3ySdknl3/+G3r1t1Xy6LrwSSSW33QbVq8NDD/kdSfx4kVTSRVollZ07bdX8009b\nGW4R8d/y5ZCTY0Um02V804uB+vfd/fdY6yTytrO0XyTxMXgwXHyxEopIMjn1VBu0nz3b70j8V1xB\nybPdfa1EBCIle/VVeO89WLLE70hEpLDQgH2mF5mMpWnTGasQHGqd1MEWP87zKqg4S4vur6+/tsWN\n06bZFEYRSS7pVmTSy3Uqf8e6wEJ+cNskQYJBK7V9/fVKKCLJqm5d6N4dxo/3OxJ/xXqRroKIxweB\nNBmKSg3PPQdbtsDQoX5HIiLFCV1nJQ06R8oslqTyOTAYqIJdffFm4DMvg5KwVavg3nvt10+VKn5H\nIyLFOfdc2LMHFizwOxL/xJJUbsAG7TcCG7Axluu9DErM/v3Qvz/k58NJJ/kdjYiUJFRkMpNX2Gud\nShK7/35bpTtzJlSMtaNSRHy1YQOcdlrqF5n0cqC+GVZU8ht3exVoWtovktKZPx+efdZ+8SihiKSO\npk2hc2dbApCJYjldjQZmAMe422tum3jkhx+sWOQzz9iCKhFJLaEB+0wUS9PmE+xCXSVtS1Yp1/31\n299aORZd/1okNaVDkUkvu7++Bfpj04grA9cAW0v7RRKbmTPh9dettpeIpKaqVa23YcwYvyNJvFiy\nUBbwNDbrC+AD4HfAeo9iireUaals3WrFIsePh65d/Y5GRMoj1YtMqkpx0VIiqQSDcNVVkJVllwgW\nkdR35pk2i/OSS/yOpPTKmlSKKyhZXAdMEFsQKXHy0ku20HHCBL8jEZF4CQ3Yp2JSKavistBALHlE\n7hd0j4PAWO/Ciqukb6msWwdnnAFvvWXdXyKSHlK5yGQiur9qYsUkU01SJ5WDB+H886FbN/jDH/yO\nRkTi7dprrcL4rbf6HUnpeDn76+dY6fuV7nk74NnSfpFE98QTNp7y+9/7HYmIeCE3N7OKTMaSVP4K\n5BCeRrwEOK+c31sPmAL8F0tYnYD6wBxgFfCm2yfkLmA1ltgujtjeEVjmXnuynDEl3NKl8OijMG5c\nas4OEZGSnXeeFZlcuNDvSBIj1gIghacPHyjn9z4JvIFd7Os0LFnciSWVE4C33XOA1kAvd5+DtZJC\nTbIRQB7Qyt1yyhlXwuzZY/PYH3/cZnyJSHoKFZnMlBX2sSSV9YQvLVwVuB1rYZRVXaALEKrjeQDY\nAXQnPPg/FrjcPe4BTAL2A+uANVjL5mdAbWC+229cxHuS3tChttJ2wAC/IxERrw0YAK+8Art3+x2J\n94pLKqGrd9wI/BZogpW/b++el9WxWGHK0cAiYCQ2CaARsNnts9k9B6s3tiHi/RtcLIW3b3Tbk14g\nABMnwvPP268YEUlvmVRksrh1KhuxQpKTsNIsBcXsW9rv7ADcBCzAxmzuLLRPkPB05nLLz8//6XF2\ndjbZ2dnx+uhS27EDBg6EkSOhYUPfwhCRBMvNhb/9za6RlIwCgQCBQKDcn1Pc7+SGwFVAb2y8YgqW\nYD4q53c2Bj7EWiwA52AD8ccBXYFNWNfWO8BJhBPOI+5+FjAM+MLtc7Lb3gebQHBDoe9LqinFAwbY\nNRZGjPA7EhFJpFCRyQ8/hOOP9zuaknkxpXgr8HcgGzgDu6zwE8Ba4OFSRxi2CfgSG5AHuBD4FCup\nHxphGABMd49nYImtKpaIWmHjKJuAndj4SgWs6GXoPUlpyhSrWvrnP/sdiYgkWtWq0K8fjE7zC4eU\nJgvVBq4AbsNaEkeX43vbAi9giWItcB1WBfkVoDk2IN8T2O72vxvIxQb1bwZmu+0dgTFAdWw2WbTS\nMUnRUpkyxUraz5gBnTr5HY2I+GHZMivZkgpFJr1aUV8d+CXWUjgb63qaBLxF+acVJ4rvSSUYhFat\n7FfK/ff7GoqI+CxVikx60f01EZtO3BOYgJXAH4AlllRJKElhxAj48Ue4916/IxERv+Xm2mXC01Vx\nWWgAMBXYlaBYvOJbSyUQgNmz4cknLakMG2bbs7PtJiKZJ1WKTOp6KkXztftr+HArz3DiiRAxs1lE\nMlj//tCxI9xyi9+RFM3LgpJSRvv2WcHIO+7wOxIRSSZ5edYFlgRziOJOScVDkybBySdDhw7q7hKR\nsHPPhR9+SM8ik7EklZ5AHfd4KDANWxEvxSgosGKRoWukKKmISEjFilZkMh0H7GNJKkOxRYbnABcA\no7DqwFKMmTNtsdOFF/odiYgko4EDYfLk9CsyGUtSOejuL8OKP76OLVqUYjz2mLVSVDBSRKJp2tQW\nQk+d6nck8RVLUtkIPI9d0+RfQLUY35exPvoI1q+Hq67yOxIRSWahAft0Esvv6JrYxa+WYldY/BnQ\nBrs6YypI+JTiK6+Erl3hppsS+rUikmL27oVmzZKzyKSX61SOx1ore7AqwqdhF9HaXtybkkhCk8qq\nVXDOOfD551CzZsK+VkRS1K23Qq1a8OCDfkdyKC/XqUzFyrK0BJ4DmmIlXCSK4cNh0CAlFBGJTW4u\njBkDBw+WuGtKiCWpFGBJ5QrgaeAOrAtMCtm0Cf7xD6tGLCISizZtoHFjmDPH70jiI5aksg/oC1yL\nzfyC8KWGJcLTT0Pfvsldz0dEkk9eHowa5XcU8RFLf9kp2NUUP8DK3h+HLYh8pLg3JZGEjKns2gXH\nHgvz58Nxx3n+dSKSRrZvh6wsWLMmeS4zroKSRUtIUnniCZg3D15+2fOvEpE01L8/nH463Hyz35EY\nL5PKCdjlg1tjF+0CCGItllTgeVLZv9+mA06bZpVHRURK6513LKF88klyLJr2cvbXaOxa9Qew69WP\nxS7aJc7LL8MJJyihiEjZnXeeFZn8+GO/IymfWJJKdezywRWAL4B8oJuHMaWUYDBckkVEpKxCRSZT\nfcA+lqSyB6gErAFuwqYWaxWGM2sWVKoEF13kdyQikuoGDEj9IpOxJJVbgBrAYOB04BrsUsOCCkeK\nSPw0a2ZFJqdN8zuSssuEU6FnA/Xz50PPnnat6SpauSMicfCPf8CIETB3rr9xeDH76zVslle0fYJA\n99J+mU88SypXXw1dusDgwZ58vIhkoL17rSz+vHn+rnnzIql8A2zAFjzOK7R/EHi3tF/mE0+Sypo1\ncNZZsG6d6nyJSHzdcgvUru1vkUkvkkpl4CKgD1bq/l9Ygvm0DPH5yZOkcuONVo7lgQfi/tEikuGW\nLoVu3exHa6VK/sTgxTqVA8BMrOZXZ2z217vYDLCMtnmzzdDQ9VJExAunnZa6RSZLmv1VDbgSGA/8\nFngSSOF5CfHxzDPQuzccfbTfkYhIusrNTc2rQhbXtHkJKyb5BjAZWJaQiOIvrt1f339vhSM//BBa\ntozbx4qIHMLvIpNedH/1A1oBN2MVindF3HaWPsTDVAIWY7PMAOoDc4BV2KWK60Xsexd2KeOVwMUR\n2ztiyW411ory3Asv2KWClVBExEv16sFll8GEFCuKVVxSqQjULuJWJw7ffTOwAptJBnAnllROAN52\nz8EKWfbAkrw4AAALoElEQVRy9znAs4Sz5wggD0t+rdzrntm/H/7yF7jjDi+/RUTEhK6zksAropdb\nLCvqvdAUuBR4gXCC6I4Vq8TdX+4e98Bmne0H1mETBjphV5+sDcx3+42LeI8nJk+2FsoZZ3j5LSIi\n5rzzrMs9lYpM+pVUnsAuS1wQsa0RsNk93uyeAxyDrZcJ2QA0ibJ9o9vuCRWOFJFECxWZTKUB+8o+\nfOdlwBZsPCW7iH2ChLvFyi0/P/+nx9nZ2WRnF/W1RZs92+5/8Yv4xCQiEouBA6FdOxg+HKpXL3H3\nMgsEAgQCgXJ/jh+1vx4G+mPrYKph4zNTgTOwJLMJ69p6BziJ8NhK6PLFs4BhWBn+d4CT3fY+wHnY\npY8jxWX21/nn2xS/a64p90eJiJRKTo5dGbJfv8R9p5cX6Yq3u4FmwLFAb2AulmRmEK5+PACY7h7P\ncPtVde9phY2jbMJmoXXC/vD+Ee+JqwULbFpfr15efLqISPHy8lKnC8yvMZVIoWbEI1hZmFXA+YRb\nJiuAV9z9TGBQxHsGYYP9q7EB/FleBHjHHXDbbapELCL+6N7dSrd89pnfkZRMpe9LsGaNlUzYsgVq\n1YpjVCIipXDLLVCnTuLqDaZS91dKGTkSTj9dCUVE/HXddTB6NBw86HckxfNj9ldKCATsVrUq/Oc/\nEJpAlp1tNxGRRGrbFho1grfeSu5ZqOr+ikF+fjipiIj4ZcQI+7E7ebL336XuLxGRNNenj62Z+/Zb\nvyMpmpJKDNTdJSLJIBWKTKr7S0QkhcydC7feCkuWQAUPz+Dq/hIRyQDZ2bBrFyxa5Hck0SmpiIik\nkGQvMqnuLxGRFLN+PbRvDxs2eFdkUt1fIiIZonlzu67TtGl+R3I4JRURkRSUm2tXhUw26v4SEUlB\ne/dC06Ywfz4ce2z8P1/dXyIiGeSII6BvX6sHlkzUUhERSVGffAK//CV8/jlUqhTfz1ZLRUQkw7Rt\nC0cfDW+/7XckYUoqIiIpLNkG7NX9JSKSwrZvh6wsWLsWGjSI3+eq+0tEJAPVqwfduiVPkUklFRGR\nFJeXZ11gydApo6QiIpLisrNh505YvNjvSJRURERSXqjIZDIM2GugXkQkDcS7yKQG6kVEMljz5nD6\n6f4XmVRSERFJE3l5/l9nRd1fIiJpIp5FJtX9JSKS4UJFJseM8S8GtVRERNJIvIpMqqUiIiK0bQtH\nHeVfkUk/kkoz4B3gU2A5MNhtrw/MAVYBbwL1It5zF7AaWAlcHLG9I7DMvfakp1GLiKQIPwfs/ej+\nauxuS4BawMfA5cB1wFbgMWAIcCRwJ9AamAicATQB3gJaAUFgPnCTu38DeAqYVej71P0lIhnlu+9s\noL48RSZTqftrE5ZQAL4H/osli+7AWLd9LJZoAHoAk4D9wDpgDdAJ+BlQG0soAOMi3iMikrGOPNKK\nTE6cmPjv9ntMJQtoD8wDGgGb3fbN7jnAMcCGiPdswJJQ4e0b3XYRkYyXm+tPF5ifSaUW8CpwM7Cr\n0GtBdxMRkTLo2tWutbJoUWK/t3Jiv+4nVbCE8hIw3W3bjI21bMK6tra47Ruxwf2QplgLZaN7HLl9\nY7Qvy8/P/+lxdnY22dnZ5QxfRCS5hYpMvvgidOhQ8v6BQIBAIFDu7/VjoL4CNmbyLXBrxPbH3LZH\nsQH6ehw6UH8m4YH6llhLZh42e2w+8C80UC8i8pPyFJlMpYH6s4FrgK7AYnfLAR4BLsKmFJ/vngOs\nAF5x9zOBQYS7xgYBL2BTitdweEIREclYoSKT06eXvG+8aEW9iEgamzwZRo6Et94q3fvK2lJRUhER\nSWN790KTJrBwIWRlxf6+VOr+EhGRBEl0kUm1VERE0tySJdCjB3z2WexFJtVSERGRqNq1g4YNYe5c\n779LSUVEJAPk5sKoUd5/j7q/REQyQKjI5GefQf36Je+v7i8RESnSkUfCpZd6X2RSSUVEJEPk5Xnf\nBaakIiKSIbp2tW6wxYu9+w4lFRGRDBEqMulla0UD9SIiGeSLL6BjRysyWa1a0ftpoF5ERErUooWV\nwp82zZvPV1IREckweXneXRVS3V8iIhlmzx5o2rT4IpPq/hIRkZhUqwZ9+nhTZFItFRGRDBQqMvn5\n5zYrrDC1VEREJGbt2kGDBvD22/H9XCUVEZEM5cWAvbq/REQyVHFFJtX9JSIipeJFkUklFRGRDJab\nG98uMCUVEZEMdv75sG1b/IpMKqmIiGSwUJHJeLVWNFAvIpLhohWZ1EC9iIiUSajI5PTp5f8sJRUR\nEYnbgL26v0RE5Kcikx9/bC0XdX+JiEiZxavIZDoklRxgJbAaGOJzLCIiKSs3F0aPhoKCsn9GqieV\nSsAzWGJpDfQBTvY1oiQWCAT8DiFp6FiE6ViEZfqxaN/eyrXMnVv2z0j1pHImsAZYB+wHXgZ6+BlQ\nMsv0fzCRdCzCdCzCdCystTJqVNnfn+pJpQnwZcTzDW6biIiUQd++MHNm2d+f6klF07pEROKofn24\n5JKyvz/VpxR3BvKxMRWAu4AC4NGIfdYAxyc2LBGRVHc8sDbVc0SpVQbWAllAVWAJGqgXEZFyuAT4\nH9YiucvnWERERERERA4VyyLIp9zrnwDtExSXH0o6Fv2wY7AUeB84LXGhJVysi2PPAA4AVyQiKJ/E\nciyygcXAciCQkKj8UdKxaAjMwrrUlwMDExZZYr0IbAaWFbNPppw3D1EJ6/7KAqoQfWzlUuAN97gT\n8FGigkuwWI7FWUBd9ziHzD4Wof3mAq8DVyYquASL5VjUAz4FmrrnDRMVXILFcizygT+5xw2Bb7Ex\n3HTTBUsURSWVUp83U31KcUgsiyC7A2Pd43nYP6BGCYovkWI5Fh8CO9zjeYRPIukm1sWxvwOmAN8k\nLLLEi+VY9AVexdZ7AWxNVHAJFsux+Bqo4x7XwZLKgQTFl0j/Ab4r5vVSnzfTJanEsggy2j7peDIt\n7YLQPMK/RNJNrP9f9ABGuOfpuvYplmPRCqgPvAMsBPonJrSEi+VYjAROAb7Cun1uTkxoSafU5810\nac7FeiIoPOc6HU8gpfmbugK5wNkexeK3WI7FX4E73b4VSP21W0WJ5VhUAToAFwA1sBbtR1h/ejqJ\n5VjcjXWLZWMLNuYAbYFd3oWVtEp13kyXpLIRaBbxvBnhJnxR+zR129JNLMcCbHB+JDamUlzzN5XF\nciw6Yt0fYH3nl2BdIjM8jy6xYjkWX2JdXj+627+xE2m6JZVYjsXPgYfc47XA58CJWAsuk2TKefMw\nsSyCjBxw6kz6Dk7HciyaY33KnRMaWeKVdnHsaNJ39lcsx+Ik4C1sILsGNnjbOnEhJkwsx+IvwDD3\nuBGWdOonKL5EyyK2gfp0Pm9GFW0R5G/cLeQZ9/onWDM/XZV0LF7ABh4Xu9v8RAeYQLH8fxGSzkkF\nYjsWt2MzwJYBgxMaXWKVdCwaAq9h54pl2CSGdDQJGzfah7VUc8nc86aIiIiIiIiIiIiIiIiIiIiI\niIiIiIiISDI5iK37WQa8AlTHViO/g639WM7haz86A88D52FFPhcDK4D73OsN3Pt3AU97G76IiCST\nyLpQ44FbgcZAO7etFrbgLnL19v3Ar7Ck8prbVgNYhZUkr4HVZvsNSiqSAtKlSrFIsnkPaAlswsqA\nAHwP/Bc4JmK/87HSKJFF+3YDH7v378YupLbX43hF4kJJRST+KmNlQJYW2p6FtT7muecNseKVhSvf\nNsC6xZZHbEvHitqShtKlSrFIMqiOjYmAVfgdFfFaLexCYDdjLRaAi4HZEft0ARYBBdhVB//rZbAi\nXlBSEYmfH4l+De8q2BUVxwPTI7bnAMMjnv8H+KVn0YkkgLq/RLxVAWuxrMAuCBa5/TSs8musnyOS\n9NRSEYmfaOMeZwPXYOMroa6xu4FvIp6H3lvUuMk6oDZ27Y8eWLfZyvKHKyIi6eIeoKffQYiIiIiI\niIiIiIiIiIiIiIiIiIiIiIiIiIiIJKP/B5SqERgs/7VDAAAAAElFTkSuQmCC\n",
       "text": [
        "<matplotlib.figure.Figure at 0x104489650>"
       ]
      },
      {
       "metadata": {},
       "output_type": "display_data",
       "png": "iVBORw0KGgoAAAANSUhEUgAAAYwAAAEZCAYAAACEkhK6AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xl8FPX9x/HXGg6JXCoockgQkUvBE1GsrkoFKYottqCi\nAlKxrYVS65GWatpa+VkrKNoqgiIWFSxFRUtBrW7RigiKKCA0gEhAQeSSKFEg+/vjM2s2m00ym2R2\n9ng/H495ZO75Mq3z2e8NIiIiIiIiIiIiIiIiIiIiIiIiIiIiWS/H7wSI+CyE/Xew3Od0iKS8Q/xO\ngGSVc4A3gd3ADuAN4HRfUwRhZ8kEw4GDwF5gDxYEv+cc6w28jL33z4BngFZx7rEW6AQ8Dnzt3GsH\n8BLQ2TnnRGAhsB0orfN/haQsBQxJlqbAi8D9wOFAG+B32EdJ6s5/gSZAc+BRLDA0d5aHgfbOsheY\nHnNtR+ybUIgF0bude7XFgszjznnfALOA67z7Z4hINjsd2FXF8eHYx+4BLAfyIXCBy3v/C/hZzL4V\nwGXO+tnAUue+bwNnRZ33GjDSWS8A/hZ1LA/7BR35YRUC/uCkcy8wD2gBPIn9on8b+xhHdKHsV/0a\n4IeVpH+Ik75o44DnnfUBwCrgC2AzcFMl9xkOvB61fZiT/lPjnHuqc79oY4D7nPXpwO+jjn0P+zdH\nOx7lMLKKchiSLGux4pLHgf5YLiNWL2AdcCRwBzC3kvNiPQVcEbXdDTgW+CdwhPP3Pmd9orMd775u\niqaGAMOwHFJHYDH2S/4ILMjd4Zx3GBYsZgItgaHAX4Guce45DyvuOT5q35VYIMK5//VYLq078KqL\ndNYDRmEf+cI4x88FVsbsG4C9m4iA87cxcBXwrovniojUiS7YL9ciYD/2C/oo59hwYEvM+Uuwj3N1\nmgDFQDtn+4/ANGf9auCtmPPfBK511hPJYbwG5Ecd/zPlP7ADKas8HwIsinnuFOD2Sv4NfwN+66x3\nwn79H+psf0xZwKjKcOy97sLqF94kfi6tB5br6RO1Lxf4HKjvbD8O7HPu9SnwHNAh5j7KYWQZ5TAk\nmdYAI7AP+4lAa8qKQKBiwPjYOac6e7EPdySXMZSyX+etgU01vG8826LWS7Cy/ejtxs56e+BM7IMb\nWa4Ejq7kvtG5pCuBZ537AQzGfv1vxIrFeleRvrew3FNLrCguNjdyPDAfK376b9T+C53t/c52GLjH\nudcxWPHeR1U8V7KAAob4ZS0wAwscEW1izmlPxSBSmaexD+5Z2C/z15z9Wyhfr1DVfYuxX9oR8VoR\nRauqCGsT8B/sgxtZmlCxriXiFewj3xMLeE9FHVuGfbBbYr/0n6kmXZVpjxWT/Z6ygBoxAAsk0QKI\nRFHAkGTpDPySsqDQDvvAL4465yjsl299rIK4MxU/YpWZj30Qf4e14Inef4LzrHpYUVEXrMVWrPew\nsv12QDPKFz9FBCpZj/VP57nDsH9PfeAM59nx7Af+jhVzHY592HGuu8pJT6TJ7MEqnluZNlhu40Hg\nkTjH+xO//qIyhwINnPWGziIZTgFDkmUvVkSzBPslvxh4n/ItfpZg5ffbsdZIl1PWsuohZ6nMN1gl\n+YWU/3W+E6tbuAkro/+Vs70zzj1eAWY76VoKvEDFXEQ4Zr2y43uBi7DcwhasHmACZR/ZeJ5y0v93\nytcNDMOKg/ZgdRlXVXJ9VX1KRmF1EAVO2vZS1krqROx/k80u75UHfIVVmoexuo4PKzlXxLX+WLl1\nIXBrnONdsA9HCRWbCuZjTQk/wP5D0i+YzDac8k1CJXluAf7P70RIdsvBmkjmYdnq96jYpLAl1j7/\nTsoHjDxgA2VBYjZlrVokMw1HAcMvkeI/kSp5WSQVaVO/ESufnQUMijlnO1ahtz9m/xfOvlys3DkX\n95Wfkp4yaYiOdPN3rBGCSJW8DBhtsPb2EZup2AqmMjuBe7GWJp9gPXRfqdPUSaqZgVU4i0iK8jJg\n1ObXYkfgF1jRVGvKepqKiIhP6nl47y2U9bzFWd9cybmxTsd6qe5wtudinZDKtR3v2LFjeP369bVM\npohI1llP+aFoXPEyh7EMayKZhzUlHIKNmRNPbJvvNVhv1kbOsb7A6tiL1q9fTzgc1hIOc8cdd/ie\nhlRZ9C70LvQuql6wUpyEeZnDOADciI2bn4MNoPYhMNo5PgXrSbsUGyOnFBiLDRy3AngCCzql2KBn\n8TobiYhIkngZMMCGnf5XzL4pUetbKV9sFe1PziIiIilAPb0zRDAY9DsJKUPvoozeRRm9i9pL98HF\nwk55nIiIuBQIBKAG33/lMERExBUFDBERcUUBQ0REXFHAEBERVxQwRETEFQUMERFxRQFDRERcUcAQ\nERFX0j5gqOOeiEhypH3AmDv3Jb+TICKSFdI+YOTnL6J794FMmTLT76SIiGQ0r0er9VxJSSl33XUj\ngwf38zspIiIZLe1zGJ9/vo9AIBAZTEtERDyS9gGjW7eLKSws8jsZIiIZz+uA0R+bbrUQuDXO8S7A\nYqAEuCnmWHNgDjZL32psytYKPvqoH1deOaqu0isiIpXwMmDkAA9iQaMbcAXQNeacHcDPgT/Huf5+\nYL5zTQ8scFRw1VUwbVodpVhERCrlZcDoBawDNgL7gVnAoJhztmPzdu+P2d8M+A7wmLN9ANgT7yGj\nR1vA2B97BxERqVNeBow2QHTlwmZnnxsdsGAyHXgXmArkxjuxe3c4/nh44YVapFRERKrlZbPa2nTB\nrgecCtwILAXuA24Dbo89saCggBYt4Oab4Ygjgpq3V0QkRigUIhQK1fo+XrZF7Q0UYHUYAPlAKXB3\nnHPvAIqBe53tVlhleAdn+xwsYAyMuS4cDocpKYF27WDxYsttiIhI5VJxTu9lQCcgD2gADAHmVXJu\nbMK3YsVZJzjbfYFVlT3o0ENh+HB45JFapFZERKrkdW+3i7HipBzgUWACMNo5NgXLSSwFmmK5j71Y\ni6pioCcwDQs264ERVKz4DkcGHywshD59oKgIGjb08F8kIpLmaprDSPfu0eHo0Wq/+10YMQKuvNLH\nFImIpLhULJJKuhtugIcf9jsVIiKZKaNyGPv3Q/v28PLL1txWREQqUg4DqF8frrsOpkzxOyUiIpkn\no3IYAJs2wSmn2N/DDvMpVSIiKUw5DMexx8LZZ8Ps2X6nREQks2RcwABVfouIeCEjA0b//rBtG7z7\nrt8pERHJHBkZMHJy4PrrVfktIlKXMq7SO+LTT6FbN/j4Y2jaNMmpEhFJYar0jnHMMdC3Lzz5pN8p\nERHJDBkbMMAmV3roIagkEyIiIgnI6IBxwQWwbx+89ZbfKRERSX8ZHTAOOcRyGWpiKyJSexlb6R3x\n+ec2qdKGDXDEEUlKlYhIClOldyVatICBA2HGDL9TIiKS3jI+hwHwxhs2KOGaNRBI93+xiEgtpWoO\noz+wBigEbo1zvAs2d3cJcFOc4znAcuCF2iSiTx8bybYO5kAXEclaXgaMHOBBLGh0A64AusacswP4\nOfDnSu4xFlgN1KphbCCg8aVERGrLy4DRC1gHbAT2A7OAQTHnbAeWOcdjtQUGYPN617og6eqr4aWX\nbIwpERFJnJcBow1QFLW92dnn1iTgZqC0LhLTrBkMHgyPPVYXdxMRyT71PLx3bYqRBgKfYfUXwapO\nLCgo+HY9GAwSDFZ++g03wA9/CLfcYgMUiohkg1AoRKgOKnG9bDPUGyjA6jAA8rHcwt1xzr0DKAbu\ndbbvAq4GDgCHAk2BfwDXxFznqpVUtNNPhz/8AS6+OKHLREQyRiq2kloGdALygAbAEGBeJefGJvzX\nQDugAzAUeJWKwaJGVPktIlIzXgaMA8CNwEKspdNs4ENgtLMAtMLqOcYB44FNQOM496qz4QOHDoXX\nX4eiourPFRGRMunejS3hIimAn/0MWraEqOoPEZGsUdMiqawMGO+/DwMGwMaNUM/Lan8RkRSUinUY\nKatHD2jfHl580e+UiIikj6wMGKDKbxGRRGVlkRTYxErHHgtLlsBxx9VxqkREUpiKpBLUqBFccw08\n8ojfKRERSQ9Zm8MAWLsWzj0XNm2Chg3rMFUiIilMOYwa6NwZTjwRnn3W75SIiKS+rA4YoMpvERG3\nsrpICuCbb6zy+7XXoGvsbB0iIhlIRVI11KCBTd86ZYrfKRERSW1Zn8MA6/F92mk2vlRubu0TJSKS\nypTDqIW8POjdG555xu+UiIikLgUMhyq/RUSqpoDhGDAAPvkEli/3OyUiIqlJAcORkwOjRqnyW0Sk\nMqr0jrJli3Xk27QJmjSps9uKiKSUVK707g+sAQqBW+Mc7wIsBkqAm6L2twNeA1YBK4Ex3iYT2rSB\n88+Hp57y+kkiIunH6xxGDrAW6AtsAZYCV2BTtUa0BNoDlwG7gHud/a2c5T1s2tZ3nHOir63THAbA\nSy/BLbdYXUYg3fNfIiJxpGoOoxewDtgI7AdmAYNiztkOLHOOR9uKBQuAYixQtPYqoRF9+8LevfD2\n214/SUQkvbgJGIcBvwWmOtudgIEu798GKIra3uzsS1QecAqwpAbXJuSQQ2D0aDWxFRGJ5WZG6+lY\ncdDZzvYnwBzAzQSndVFe1Nh53lgsp1FOQUHBt+vBYJBgMFjrB44YAZ06wcSJcPjhtb6diIivQqEQ\noVCo1vdxU4b1DnAasBz7lQ+wAujp4treQAFW8Q2QD5QCd8c59w4sINwbta8+Fpj+BdwX55o6r8OI\nuPJKOPNMGDvWk9uLiPjGyzqMr4FGUdsdnX1uLMOKsPKABsAQYF4l58YmPgA8CqwmfrDwVKTnt0fx\nSEQk7bgJGAXAAqAt8BTwKvGbx8ZzALgRWIh9+GdjldejnQWsJVQRMA4YD2zCiqH6AMOA87HczXLK\nciqe+853rJXUokXJeqKISGpzmyVpgRUvAbwFfO5NchLmWZEUwOTJsHgxPP20Z48QEUm6mhZJub2g\nJ1asVI+yiuy5iT7MA54GjF27oEMH+N//4KijPHuMiEhSeRkwpgMnYT2uS6P2j0j0YR7wNGAAjBxp\nc3/f6rYQTkQkxXkZMFYD3ambJrJ1zfOA8fbbMHQorFtnfTRERNKdl62klgLdEr1xpjjjDGjeHF5+\n2e+UiIj4y02ECWJNYbdS1pw2DPTwKE2J8DyHATbk+YIF8Oyznj9KRMRzXhZJrceavK6kfB3GxkQf\n5oGkBIy9e+HYY+GDD6BtW88fJyLiKS+LpD7DchgbsCARWbJGkyZWj/Hoo36nRETEP24izF+B5sAL\nwDfOvjBZ0Kw22ooVMHAgfPQR1HMzApeISIryMoeRiwWKi7BRagcClyT6oHTXs6cVR82f73dKRET8\nke5TBCUthwEwYwbMnq2gISLpzeuOe9EiX+iRiT7MA0kNGPv2Qbt2sHSp9QAXEUlHXhZJ/RMbYvxF\n4N9AM+DLRB+UCRo1gquvhqlTqz9XRCTT1KRI6hDgv8BZdZyWmkhqDgNgzRoIBmHTJmjQIKmPFhGp\nE8mc0/sEoGUNrssIXbpA167w3HN+p0REJLncBIxiYK+zfIE1r83qofgikyuJiGQTNwGjMdDEWZpi\nM+j9w+X9+wNrgELiB5kuwGKgBLgpwWt98/3vw6pVVjwlIpItqirDOo2qR6h9t5p75wBrgb7AFmwQ\nwyuwGfciWgLtgcuAXZTN5+3mWvChDiMiPx++/homTvTl8SIiNVbTOoyq+izfS9UB4/xq7t0LWEfZ\nMCKzgEGU/+hvd5bv1eBaX/34x9CrF/zxj9Z6SkQk01UVMIK1vHcbbK7uiM3AmUm4NimOOw5OPx2e\neSbMhx/ew4QJN0eitohIRnIzKlID4CfAuc52CHgY2F/NdbUpK0rFyZoquOEGuPXWhXz66aecccZL\nDB7cz+8kiYh4xk3AeMg57y9YmdfVzr5R1Vy3BWgXtd0Oyym44fragoKCb9eDwSDBYNDlI2pnypSZ\n3H//LNav78nBgxPJzx/P7bc/wJgxQxk9elhS0iAi4kYoFCIUCtX6Pm7KUN6n4mRJ8fbFqodVXF8I\nfAK8TfyKa4ACrNnuvQle61uldzgcZs6cBYwevYhduybQrl0+Eyeex+DB/VQ0JSIpzYtK74gDwPFY\nJTRAR2efm+tuBBZirZ4exT74o53jU4BWWAuoptjkTGOx6WCLK7k2ZQQCAQKBAAcOlHDYYb9k27bS\nb/eJiGQiN1+3C7EBCD9ytvOAEcCrHqUpEb7lMAAmTJjKCSccS69eF9Gt20uMGFHE5MnVldSJiPjL\ni9Fq5wNPAc9hFdxdnP1rsY52qcDXgBHtiSfgz3+2kWwbNvQ7NSIilfNiLKlHsMmSPgL+BhyHFQul\nSrBIKVdfbU1tf/c7v1MiIuINNxHmMGyGvaHYCLXzgaeBlzxMl1spk8MA2LYNevSA55+H3r39To2I\nSHxeTqAUrScwAzgJq4z2W0oFDIC//x3Gj4flyyE31+/UiIhU5GXAaAX8CMthHAPMxnIYKxJ9mAdS\nLmAAXHEFtGoFkyb5nRIRkYq8CBjXY0GiCzY67dPYyLKp9IVOyYCxYwecdBI8/TScd57fqRERKc+L\ngDEdayX1KnCwZsnyXEoGDIAXXoCxY2HFCmjSxO/UiIiUSVYdRqpJ2YABMHKkTeOqyZZEJJUoYKSg\nPXusaGrqVOincQlFJEUkc05vcalZM3jsMRg1Cnbv9js1IiK14yZgnIWN9RTRlBSbmyKV9e0Ll15q\n9RkiIunMTZbkPeBUbHBAsP4Xy4BTvEpUAlK6SCqiuBhOPhnuvRcGDfI7NSKS7bwukiqNWj9IanTa\nSxuNG8Pjj8NPfgKff+53akREasZNwPgIGAPUx2bfGwts8DJRmeicc+DKKy1opEGmSESkAjcB4wag\nDzYL3magN9apTxJ0552wahXMnu13SkREEqdmtUm2dCkMHAjvvQfHHON3akQkG3nRD+NW4G7ggTjH\nwlgxVXX6A/dhdR7TnPvFmgxcDHwFDAeWO/vzgWFY/ckH2KRNX8emI90CBsBvf2sBY9480AR9IpJs\nXlR6r3b+voO1ioos7zhLdXKAB7Gg0Q2bk7trzDkDsOlfO2HFXA85+/OAH2OtsyIj4w518cy08Nvf\nQlGRVYSLiKSLqub0/hHwAtAcyyUkqhc2D/hGZ3sWMIjyc3Nfig2XDrDEedbRwBfYLH+5WKusXKwO\nJSM0aGAz9F14oS3HHut3ikREqldVDuM0oDUwEjgizlKdNkBR1PZmZ5+bc3YC9wKbgE+A3cArLp6Z\nNnr0gHHj4LrroLS0+vNFRPxWVcB4GPg30JmyYqjIsszFvd1WLsQrR+sI/AIrmmoNNAaucnm/tHHL\nLfDFFxqcUETSQ1VFUpOd5WGsaW2itgDtorbbYTmIqs5p6+wLAm8CO5z9c4GzgSdjH1JQUPDtejAY\nJBgM1iCp/qhXD2bMsD4a/fpBx45+p0hEMlEoFCIUCtX6PlXVkjfF6hKOJH5uYWc1964HrAUuxIqV\n3sYqvqPrMAYANzp/e2N1Jb2Bk4GZwBlACfC4c/1fYp6Rlq2kYk2aBHPnQigEOepDLyIe86JZ7T+B\n72GV1vG+yh1c3P9iyprVPgpMAEY7x6Y4fyMtqb7Ems6+6+y/BbgWa1b7LjAKqwiPlhEBo7QUzj/f\nBim86Sa/UyMimU7zYaS5DRugVy94/XXoGtv4WESkDnk5+OD3seauEc2ByxJ9kFTtuONs6JBrroED\nB/xOjYhIRW4izAqgZ8y+97B6Br9lTA4DbFDCfv3g3HNh/Hi/UyMimcrLHEa8m6pq1gOBADz6KEye\nbEOHiIikEjcB4x1gItY34nhgEu6GBpEaaNcO7rkHrr0Wvo4dOUtExEduAsbPsdZJs7HhPUqAn3mZ\nqGx3zTWQlwe//73fKRERKZNIGdZhWNPXVJJRdRjRtm61aV2ffx7O1AzqIlKHvKzDOBsbuXaNs90T\n+GuiD5LEtGpldRnXXgv79vmdGhERdxHmbeBy4HngFGffKqC7V4lKQMbmMCKGDoXWrWHiRL9TIiKZ\nwsscBtiosdHUUyBJ/vIXm9J10SK/UyIi2c5NwNiEzekN0AD4FeXHgxIPHXmkjWY7fDgUF/udGhHJ\nZm6yJC2B+4G+zvkvYdOz7qjqoiTJ+CKpiBEj4NBD4aGHqj9XRKQqGksqw+3ZY5MuTZ0KF13kd2pE\nJJ15WYfREZuq9XNgO1b5fVyiD5LaadYMpk2DUaNg926/UyMi2chNwHgKeAY4Bpv97u/A014mSuL7\n7ndh4ED4xS/8TomIZCM3WZL3gR4x++INSOiHrCmSiiguhp49bdKlSy/1OzUiko68rMO4G9hNWa5i\nCHA48Cdnu7qZ97yUdQEDbM6MIUPg/fehRQu/UyMi6cbLgLGR+DPu4eyvqj6jP2Uz7k3Dgk+sydjM\nfF8Bw4Hlzv7mzjXdneeMBN6KfX42Bgywmfk2b7Y+GiIiiUjFVlI52JzefYEtwFKqntP7TKz5bm/n\n2AzgP8Bj2PzghwF7Yp6RtQFj3z449VQoKLDchoiIW160kuqFVXRHXAvMw3IER7i4dy9gHZZD2Y+N\ndDso5pxLscAAsATLVRwNNAO+gwULsJ7lscEiqzVqBDNmwJgx8OmnYW677U9ka/AUkeSoKmBMASIz\nMpwL/B/2cf8CeMTFvdsARVHbm5191Z3TFuiANeGdDrwLTAVyXTwzq/TqBT/+MVxyyUL++tdPmTv3\nJb+TJCIZrF4Vxw6hrEJ7CBZA/uEsK1zc2+3P3dhsUdhJ16lYcdVSrB7kNuD22IsLCgq+XQ8GgwSD\nQZePTX9Tpsxk7txZrF/fk2++mUh+/nhuv/0BxowZyujRw/xOnoikiFAoRCgUqvV9qirDWomNTrsf\nq4u4HqtTAHej1fYGCrCKb4B8oJTyFd8PAyGsuApsCPXznHQtxnIaAOdgAWNgzDOytg4DIBwOM2fO\nAsaMWcTWrRNo2DCfMWPO4667+lGvXrp34hcRr3hRh/E0FiDmYS2YXnf2d8Ka2VZnmXNuHjZo4RDn\nXtHmAdc4672d+24DtmJFVSc4x/piQUqiBAIBAoEAX35ZQteuvyQQ2Mezzwbo3DnA5Mmwd6/fKRSR\nTFJVwPgjcBNWj3AOljsAi0o/d3HvA1iR0kJsAqbZWAup0c4CMB/YgFWOTwF+GnX9z4EnseKvHsBd\nLp6ZdQoLi5g+vT+rVt3LzJkXc911RcycCW+8YdO8/upX8PHHfqdSRDJBupdbZHWRVHU+/hgeeACm\nT4cLLoBx4+CssyCQ7v+ri0itpGI/jGRQwHBh714LGvffbz3Dx42DwYOhfn2/UyYiflDAkGodPAgv\nvmjjUG3YADfeaM1yDz/c75SJSDJ5PUWrZICcHBg0CEIheO45WLkSOna0wFFY6HfqRCTVKWBkqVNP\nhSeesKDRvDn06WOj3776KijTJiLxqEhKAPjqK5g5E+67Dxo0sDk3rrgCGjb0O2UiUtdUhyF1IhyG\nhQstcKxYAT/5CdxwAxx1lN8pE5G6ojoMqROBAPTvDwsWwCuv2BDqnTvb1LArV/qdOhHxkwKGVKp7\nd3jkEfjf/6wT4EUX2TJ/PpSWVnu5iGQYFUmJa19/bRM2TZoEJSUwdixccw3kahxhkbSiOgxJmnAY\n/vMfCxxvvml9OX72M2gTO3i9iKQk1WFI0gQCEAzC88/D4sVQXAwnnQTDhsE77/idOhHxinIYUid2\n74Zp02zsqvbtbfiRSy+1zoIiklpUJCUp4cABmDvXiqs++8ymkB0xApo29TtlIhKhIilJCfXqwY9+\nZEVVTz5pdRwdOsAvfwkbN/qdOhGpDQUM8Uzv3taqavlyK5o67TS4/HL47381/IhIOlKRlCRNcTE8\n/rgNs3744VbPcfnlGmZdJNlStUiqPzZPdyFwayXnTHaOr8DmEI+WAywHXvAqgZI8jRvbyLhr1sD4\n8dYp8Ljj4O67YedOv1MnItXxMmDkAA9iQaMbcAXQNeacAcDx2Nzf1wMPxRwfi03vqmxEBsnJsRZU\nr70G8+bB6tU2zPpPf2q9ykUkNXkZMHphc3VvBPYDs4BBMedcCsxw1pcAzYGjne22WECZRvoXnUkl\nTjkFZsywoHHkkXDOOTBwIPz736rnEEk1XgaMNkBR1PZmZ5/bcyYBNwMatSgLHHMM/OEPNg/5ZZdZ\nc9yTT7apZUtK/E6diADU8/Debn8fxuYeAsBA4DOs/iJY1cUFBQXfrgeDQYLBKk+XFNeokY2Me911\n8PLL1p8jP79smPWjj67+HiJSXigUIhQK1fo+Xhb19AYKsDoMgHwst3B31DkPAyGsuAqsgjwIjAGu\nBg4AhwJNgX8A18Q8Q62kssDq1day6pln4Pvft9ZVJ53kd6pE0lcqtpJahlVm5wENgCHAvJhz5lEW\nBHoDu4GtwK+BdkAHYCjwKhWDhWSJbt1gyhSbd7xjR+jXD/r2hX/+U8OsiySTlwHjAHAjsBBr6TQb\n+BAY7SwA84ENWOX4FOCnldxL2QihRQv4zW+sx/jw4XD77RZMHnoIvvzSzgmHw9x2259QzlOk7qV7\n6yMVSWWxcBhef93qOd54w+o+8vIWcPPNC5k+vT+DB/fzO4kiKUmDD0pWu/POmUyaNIvdu3tSWnon\nRxwxnkaNVjBq1FB+85th6k0uEkUBQ7JaOBxmzpwFjBu3iC1bJtC0aT6dO5/Hrl392Lw5wPHHw4kn\n2tK9u/3t0EHDr0t2qmnA8LJZrUjSBAIBAoEAX3xRQrduv6SoqJRbbw0weHCAr76y4UhWrYKVK2Hq\nVFvfvh26dKkYSNq1s0miRKQ8BQzJGIWFRUyf3p8f/OAi5s59icJC6xOamwunnmpLtC++sCa7kUDy\n8sv2t7i4LHhE/p54ovUBUSCRbJbu//dXkZTUuZ07y4LIypW2/sEHdiw2N9K9uw1pIpJOVIch4qFw\nGLZtqxhIVq6Eww6rGEi6ddMsg5K6FDBEfBAOQ1FRWfCI/P3wQ+s3Elus1bWrDX8i4icFDJEUcvCg\ndTCMzY0UFkLbthUDyQknQIMGfqdasoUChkga2L8f1q0rH0RWrrRReo87rmL9SMeONk+6SF1SwBBJ\nYyUlsHbQWkdEAAAJCklEQVRtxUCydSt07lw+N9K9O7RvD4dUM7BPOBwmP/8eJky4OfKBEAEUMEQy\nUnGx1YfEVrbv3m0V67GBpHXrsqa/c+YsYORIDZMiFSlgiGSR3bstcMQGkm++gSOPnMmOHbOoX78n\n27ffSfv242nYcAXjxg3lhhuG+Z10SQEKGCLC9u2wcmWYv/1tAc88s4gvv5xAvXr5HHroeZSU9OOo\nowK0agWtWllHxMh67NKkiTopZjINDSIitGwJ558fYMeOAHPmlA2TMn16gEsuCfDZZ1YvEr0UFtqo\nv9H7Dh6MH0hig8zRR8Ohh/r9r5ZkUcAQyUDxhklp0MCa9LZtW/31xcXWUTE2uLzzTvntbdts6JXK\ncirRAaZlSw32mO7SPdOpIikRH4XDsGtXxcASHVAi6zt32jAqlQWX6ADTvHndFompxVh5qVyH0R+4\nD8gBplF+Tu+IycDFwFfAcGA5NkXrE8BR2Ix7jzjnRVPAEEkTBw5YHUtlwSU6wJSUVF3HEh1gcnOr\nf7ZajJWXqgEjB1gL9AW2AEuBK7CpWiMGYFO5DgDOBO7H5vdu5SzvAY2Bd4DLYq5VwBDJQPv2xS8S\ni7c0aFB5QFm+fCbz588CerJhw5106jSe+vVXMGbMUEaPzt4WY6la6d0Lm697o7M9CxhE+Y/+pcAM\nZ30J0Bw4GtjqLADFzjWtY64VkQzUqBHk5dlSlXDYhqmPF0gKC+HTT6/i4MEjKSpaBAQoLCylWbMb\nmTSpH9OnQ7NmVvwV/beq9aZNq+8wmcm8DhhtgKKo7c1YLqK6c9oC26L25QGnYAFFRASweo7Ix7xz\n57hnMGdOgJEjS2jXzlqM3XNPgD59AuzZA3v2WJ+W6PVPPom/f88eawzQuLH7ABNvPTfX3ybLtSmV\n8TpguE1Z7OuLvq4xMAcYi+U0yikoKPh2PRgMEgwGE0qgiGS2eC3GunWr2b0OHoS9eysGkuj1zz+3\n8cIqO75/f80CTfR6ogNVhkIhQqEQAKtXr6vZPx7v6zB6AwVYxTdAPlBK+Yrvh4EQVlwFsAY4D8th\n1AdeBP6FVZzHUh2GiKSVb74pCyKVBZXq1uvVSzzYvPLKTGbPnkVpaU/Wr78LUrAOYxnQCStS+gQY\nglV6R5uHVXrPwgLMbixYBIBHgdXEDxYiImmnQQPrk9KyZc2uD4etUUB1QWXr1vL7d+++ip07j2Tn\nzkU1TrvXAeMAFgwWYi2mHsUqrUc7x6cA87EWUuuAL4ERzrE+wDDgfayZLVgOZYHHaRYRSVmBgNWD\n5ObCMcckdOW39Tl799bw2TW7LGWoSEpExKUJE6ZywgnHcvnl/SEF+2F4TQFDRCRBNe2HkcUtikVE\nJBEKGCIi4ooChoiIuKKAISIirihgiIiIKwoYIiLiigKGiIi4ooAhIiKuKGCIiIgrChgiIuKKAoaI\niLiigCEiIq4oYIiIiCsKGCIi4orXAaM/NuVqIXBrJedMdo6vAE5J8FoREUkSLwNGDvAg9uHvhk3N\n2jXmnAHA8dg0rtcDDyVwrUSJTPAuehfR9C7K6F3UnpcBoxc27epGYD82Z/egmHMuBWY460uA5kAr\nl9dKFP3HUEbvoozeRRm9i9rzMmC0AYqitjc7+9yc09rFtSIikkReBgy3c6em+zSxIiJSS72BBVHb\n+VSsvH4YGBq1vQY42uW1YMVWYS1atGjRktCyjhRTD1gP5AENgPeIX+k931nvDbyVwLUiIpJBLgbW\nYtEs39k32lkiHnSOrwBOreZaERERERGR2qtNB8BMU927uAp7B+8D/wV6JC9pSee2c+cZwAHgB8lI\nlE/cvIsgsBxYCYSSkip/VPcuWmB1pO9h72J40lKWXI8B24APqjgn476bOVixVB5Qn+rrQs6krC4k\n07h5F2cBzZz1/mT3u4ic9yrwIjA4WYlLMjfvojmwCmjrbLdIVuKSzM27KAAmOOstgB1YvWmm+Q4W\nBCoLGAl/N9NhLKmadgA8OknpSyY372IxsMdZX0LZByLTuO3c+XNgDrA9aSlLPjfv4krgH1ifJoDP\nk5W4JHPzLj4FmjrrTbGAcSBJ6Uum14FdVRxP+LuZDgGjph0AM/FD6eZdRLuOsl8Qmcbt/y8GUTbk\nTDgJ6fKDm3fRCTgCeA1YBlydnKQlnZt3MRXoDnyCFcWMTU7SUk7C3810yIa5/Y88tgNgJn4cEvk3\nnQ+MBPp4lBa/uXkX9wG3OecGyNxOom7eRX2sFeKFQC6WE30LK7/OJG7exa+xoqog0BF4GegJ7PUu\nWSkroe9mOgSMLUC7qO12lGWrKzunrbMv07h5F2AV3VOxOoyqsqTpzM27OA0rkgArq74YK6aY53nq\nksvNuyjCiqH2Ocsi7COZaQHDzbs4G/ijs74e+AjojOW8sklGfjdr0wEw07h5F8diZbi9k5qy5Eu0\nc+d0MreVlJt30QV4BasUzsUqQrslL4lJ4+ZdTATucNaPxgLKEUlKX7Ll4a7SO6O+m7XpAJhpqnsX\n07BKvOXO8nayE5hEbv5/EZHJAQPcvYtfYS2lPgDGJDV1yVXdu2gBvIB9Kz7AGgRkoqexeppvsBzm\nSLL3uykiIiIiIiIiIiIiIiIiIiIiIiIiIiIiyXUQ69fyAfAM0AjrJfsa1rdhJRX7NvQGHgHOwwaE\nXA6sBm53jh/pXL8XeMDb5IuISLJEjzM0ExgHtAJOdvY1xjqLRfcq/h3wfSxgvODsywX+hw07nYuN\n9TUaBQxJA+kwWq1IqnkDOB7Yig09AVAMfAi0jjrvAmw4jugB3r4C3nGu/wqb5Oprj9MrUicUMEQS\nUw8beuL9mP15WK5hibPdAhvoMHYE1COxoqqVUfsycWRlyUDpMFqtSCpohNVBgI30+mjUscbYJE1j\nsZwGwEXAwqhzvgO8C5Ris7196GViRbyggCHizj7iz3lcH5vJbibwXNT+/sC9UduvA5d4ljqRJFCR\nlEjNBbCcxmpssqbo/T2wEUDd3kck5SmHIeJOvHqGPsAwrD4jUlz1a2z+8OVR54UruR5s7ukm2NwN\ng7CirDW1T66IiKSD3wA/8jsRIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIiIuKP/wchKwZgozaTFAAA\nAABJRU5ErkJggg==\n",
       "text": [
        "<matplotlib.figure.Figure at 0x1045fb150>"
       ]
      },
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "The critical pressure is 26.842561 atm\n",
        " The critical velocity is 323.738941 m/sec\n",
        " The critical specific volume is 0.028541 cubic meter/Kg\n",
        " The critical mass velocity is 11343.108938 Kg/sq meter sec\n",
        " Mass rate of flow through nozzle is 0.890886 Kg/sec\n"
       ]
      }
     ],
     "prompt_number": 3
    },
    {
     "cell_type": "heading",
     "level": 2,
     "metadata": {},
     "source": [
      "Example 9.4  Page No : 156"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "import math \n",
      "\n",
      "#Given\n",
      "A1 = 0.1;#Inlet area in sq meter\n",
      "u1 = 60.0;#inlet velocity in m/sec\n",
      "v1 = 0.185;#inlet specific volume in cubic meter/Kg\n",
      "H1 = 715.0;#inlet enthalpy in Kcal/Kg\n",
      "H2 = 660.0;#exit enthalpy in Kcal/Kg\n",
      "v2 = 0.495;#exit specific volume in cubic meter/Kg\n",
      "g = 9.81\n",
      "\n",
      "#To calculate the area at exit of nozzle and hence decide the type of the nozzle\n",
      "#From the first law\n",
      "u2 = ((u1**2)-(2*g*(H2-H1)*427))**(1/2);\n",
      "W = (u1*A1)/v1;#Mass rate of gas in Kg/sec\n",
      "A2 = (W*v2)/u2;#Area at exit of nozzle\n",
      "if(A2 < A1):\n",
      "    print \"The nozzle is convergent\";\n",
      "else:\n",
      "    print \"The nozzle is divergent\";\n",
      "#end\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "The nozzle is divergent\n"
       ]
      }
     ],
     "prompt_number": 4
    }
   ],
   "metadata": {}
  }
 ]
}