{ "cells": [ { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXIAAAEKCAYAAAAPVd6lAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjMsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+AADFEAAAgAElEQVR4nOydd3ib1dn/P0fLkmecxDHZE0JC9iQhhJCE0QQIUCibN/zahlJKoQMCFFpTKC2jtFDoYKZQIFBGeBkvJZSYDRlk7+kMO45jJ97a5/fH40eWLMm2PIjl3J/r8mVJz9HRLQe+uvU997mP0lojCIIgJC+WYx2AIAiC0DpEyAVBEJIcEXJBEIQkR4RcEAQhyREhFwRBSHJEyAVBEJKcNhFypdTPlFIblVIblFIvK6WcbTGvIAiC0DStFnKlVG/gp8AErfUIwApc3tp5BUEQhObRVtaKDXAppWxAKlDYRvMKgiAITWBr7QRa6wNKqYeBvUAt8IHW+oPGntO9e3c9YMCAFr1edXU1aWlpLXpueyJxJYbElRgSV2J01LigdbGtWrXqsNY6p+HjrRZypVQ2MA8YCBwF/q2Uulpr/a8G4xYACwByc3N5+OGHW/R6VVVVpKenty7odkDiSgyJKzEkrsToqHFB62I788wzC2Je0Fq36ge4FHgm7P61wF8be8748eN1S1m2bFmLn9ueSFyJIXElhsSVGB01Lq1bFxuwUsfQ1LbwyPcCpyqlUpVSCpgFbG6DeQVBEIRm0Goh11p/DbwGfAOsr5vzydbOKwiCIDSPVnvkAFrr3wC/aYu5BEGIjc/nY//+/bjd7mMdSoisrCw2b+54X8A7alzQvNicTid9+vTBbrc3a842EXJBENqf/fv3k5GRwYABAzBczGNPZWUlGRkZxzqMKDpqXNB0bFprSktL2b9/PwMHDmzWnLJFXxCSBLfbTbdu3TqMiAvtg1KKbt26JfTNS4RcEJIIEfHjg0T/nTuNkG8r3cZbW9461mEIwnHJueeeS5cuXTjvvPOirl1yySXs2rULgPvvvz/0uNfrZfr06fj9/m8tznB+8IMfsGnTprjXvV4vc+bMYdasWdx8883fYmSJ02mEfP6S+Vz4yoXsKNtxrEMRhOOOW2+9lRdeeCHq8Y0bNxIIBBg0aBAQKeQOh4NZs2bxyiuvfGtxhvP0008zfPjwuNcdDgfvvfce//3vf3n00Ue/xcgSp9MI+aYS45NVhFwQ2ocVK1YwatQo3G431dXVnHLKKaGMdtasWTEX8F588UXmzZsHwO23305tbS1jxozhqquuAuDCCy/kxRdfTCiORYsWMW/ePM4991yGDh3KPffcE7r2yCOPMGLECCZPnsyf//xnwNgSP3fuXEaPHs2IESNCHxwzZsxg5cqVALz//vuMGzeO0aNHM2fOHADefvttJk+ezNixY5k9ezbFxcUAlJWVceGFFzJq1ChOPfVU1q1bl1D87UGnqVrxB42vZ6U1pcc4EkHonEycOJELLriAu+66i9raWq6++upGM1qAzz//nCuuuAKAP/zhDzz++OOsWbMmdH3EiBGsWLEi4ViWL1/Ohg0bSE1NZeLEicydOxelFM899xxff/01FRUVzJ49mzPOOINdu3bRq1cv3n33XQDKy8sj5iopKeH666/nk08+oX///pSVlQEwbdo0vvrqK5RSPP300zz44IP88Y9/5De/+Q1jx45lyZIlfPTRR1x77bUR7+lY0PmEvFaEXOj83HILtLV2jBkDdUlsXH79618zceJEnE4njz32GDU1NY2OLyoqIicnqsdTCKvVisPhSLhc8KyzzqJbt24AXHzxxXz22WcopbjoootIS0sjGAxy8cUX8+mnn3Luuefyy1/+koULF3Leeedx+umnR8z11Vdfcfrpp9O/f38AunbtChjlnpdddhlFRUV4vd5QKeBnn33G66+/DsDMmTMpLS2lvLycrKysZsff1nQaa8VuNQrnJSMXhPajrKyMqqoqKisrm1Ue53K5mhzn8XhwOiPPonniiScYM2YMY8aMobAwuit2w6oOpZTZ6ymKk046iVWrVjFy5EjuuOMOfvvb30Zcj/e8m266iZ/85CesX7+ef/zjH6H3EWv8sa4m6hQZudYat9/4I1d6K49xNILQ/jSVObcXCxYs4N5772X37t0sXLiQ3//+942OHzZsGDt27MBsW2232/H5fKEdi6WlpeTk5ETtYLzxxhu58cYb4867dOlSysrKcLlcLFmyhGeffRaLxcL8+fO5/fbbqa6u5s033+SFF16gsLCQrl27cvXVV5Oens6iRYsi5poyZQo33XQTBQUFIWula9eulJeX07t3bwD++c9/hsZPnz6dF198kbvvvpv8/Hy6d+9OZmZmc/+E7UKnEHJvwBuyVmp8jX/VEwShZTz//PPYbDauvPJKAoEAU6dO5eOPPw7ZFVu2bKGqqoo+ffrwzDPPcM455zB37lzy8/OZPXs2YHwQjBo1inHjxvHiiy+ybNmy0OJiIkybNo1rrrmGHTt2cOWVVzJhwgQA5s+fz6RJkwgGgyxYsICxY8fyn//8h1tvvRWLxYLdbudvf/tbxFw5OTn8/e9/58ILL+TQoUOMHTuWd955h7y8PC699FJ69+7Nqaeeyu7duwHIy8vjuuuuY9SoUaSmpkaI/DEjVkvE9v5p6za2h6sPa/LQ5KGvfuPqFs/dGjpq20yJKzE6clybNm061mFEUVFR0ej1mpoaPXnyZO33+2Nev+iii/SWLVsSes3nnntO33jjja2KKx4///nP9dGjR1v03ObS3Nhi/XvTjm1sjznVvurQbcnIBaHj4HK5uOeeezhw4EDUNa/Xy4UXXsjQoUOPQWTRXHHFFbz99tv4fL5jHUrCdAprpcpbFbotQi4IHYtzzjkn5uMOh4Nrr7024fnmz5/P/PnzWxlVNC+//HKbz/lt0Sky8nDxrvZWNzJSEASh89EphNwb8AJgs9gkIxcE4bijTYRcKdVFKfWaUmqLUmqzUmpKW8zbXDx+DwDZzmwRckEQjjvayiN/FHhfa32JUsoBpLbRvM3CEzCEPCMlI3RbEATheKHVGblSKhOYDjwDoLX2aq2PtnbeRDCtlQxHRui2IAjC8UJbWCuDgBLgOaXUaqXU00qptDaYt9mY1kpmSmbotiAI3w5r1qxhypQpnHLKKYwaNSqqLe2x7keen58fs096c5gzZw5Hj0bnpXl5eTz88MMJzZWent6iGJpDW1grNmAccJPW+mul1KPA7cDd4YOUUguABQC5ubnk5+e36MWqqqqinrum2Oge5K3yUuOpafHcrSFWXB0BiSsxOnJcWVlZVFZ2rBYUgUCAYDDIX//6V4YMGUJRURHTp09n6tSpdOnShc2bN+PxeMjJyaGyspL777+fm266KfT8adOmsWjRIi677LI2j8v8W9XU1OD3+1v0tzM/lBo+1+PxYLfbE56zsrIyIrbGcLvdzf9vMdYuoUR+gBOAPWH3Twfebew5bb2z85lvntHkoS/792XaeZ+zxXO3ho68I7AjInElRkfY2bl8+XI9cuRIXVtbq6uqqvTw4cP1V199FTVu1KhRetu2bVprre+44w793HPPaa21XrhwobZYLHr06NH6yiuv1FprvWbNGv2d73wnoTgmTZqkN2zYELp/xhln6JUrV+qqqip93XXX6QkTJuhRo0bpJUuWaK2Nv93cuXO11lqXlpbqefPm6ZEjR+rJkyfrtWvXaq21rqys1PPnz9cjRozQI0eO1K+99prWWuv+/fvrkpISrbXW9913nz7ppJP0rFmz9OWXX64feughrbXWTz75ZOg1L774Yl1dXa211nrXrl361FNP1RMmTNB33XWXTktL01obOzsffPBBPWHCBD1y5Ej961//Oub7TGRnZ6szcq31QaXUPqXUUK31VmAWEP/8pHbAtFMyHBl4/B601se8G5kgtCe3vH8Law62bR/bMSeM4c/nxu/G1Zx+5MuXL8fr9TJ48GCgffqRX3755bz66qvcc889FBUVUVhYyPjx47nzzjuZOXMmzz77LPv27WPWrFmhHi8m8XqJ33vvvWRlZbF+/XoAjhw5EvG8VatWsXjxYlavXo3f72fcuHGMHz8eMNro/vCHPwTgrrvu4plnnuGmm27i5ptv5oYbbuDaa6/liSeeCM313//+l+3bt7N8+XK01lxwwQV88sknTJ8+PaG/QzhtVUd+E/CiUmodMAa4v4nxbUposTMlA40moAPf5ssLwnHDr3/9a5YuXcrKlSu57bbbIq4VFRVxzTXX8Nxzz2GxWEKPNbcfeXP53ve+x7///W8AXn31VS699FIAPvjgA/7whz8wZswY5s6di9vtZu/evRHP/eyzz7jmmmuAyF7iH374YUS3xezs7Ijnffrpp1x00UWkpqaSmZnJBRdcELq2YcMGTj/9dEaOHMmLL77Ixo0bgcgPMfM1AT766CM++OADxo4dy7hx49iyZQvbt29v9vuPRZuUH2qt1wAT2mKulmCWHGamGK0kPX4PNken6D4gCDFpLHNuT8x+5D6fL6LPeEVFBXPnzuW+++7j1FNPDT3emn7kTz31FADvvfcevXr1Cl3r3bs33bp1Y926dbzyyiv84x//AAyb+PXXX2fo0KERB1WYR7SZYxpi9jJv6lt8vOvz589nyZIljB49mkWLFkX42rGeo7Xmjjvu4Prrr2/09RKhU+3szHAY/3BSSy4I7YPZj/yqq65i4cKFgFF9ctFFF3HttdeGsmMTsx+5idmP3KSxfuRr1qxhzZo1ESJucvnll/Pggw9SXl7OyJEjAaOny1/+8peQWK9evTrqeWYvcSCil/jZZ5/N448/HhrX0FqZPn06b775JrW1tVRWVvL222+HrlVWVtKzZ098Pl/E+aOnnXYaixcvBoh4fNasWTz77LNUVRk9og4cOMChQ4eiYk2ETiHkHr8Hi7KQajf2ITVWSz7t2Wnc/dHdca8LghCb8H7kt99+OytWrODjjz/m1Vdf5ZNPPmHRokWhU31MH9zsR25i9iM3D19uaT/ySy65hMWLF/O9730v9Njdd9+Nz+dj1KhRTJ48mbvvjv7/PC8vj5UrVzJq1Chuv/32UC/xu+66iyNHjjBixAhGjx7NsmXLIp43btw4LrvsMsaMGcN3v/vdiOPi7r33XiZPnsxZZ53FySefHHr80Ucf5YknnmDixIkR54TOmjWLK6+8kilTpjBy5EguueSS1lcjxVoBbe+ftq5aue2D27TzPmeoemXPkT0xn7v36N5Q3/K2piNXO3REJK7E6AhVK7E4Fv3I2yKuY4n0I4+DJ+AhxZpCijUFiJ+R76vY922GJQjHPcnUjzyZ6RQrgh6/B4fVgcPqMO7H8cgPVdf7UN6ANzReEIT2o637kQvRdIqM3BvwkmJLIcVmZOTxtumXVJeEbh+pPRJzjCB0ZHScE9+FzkWi/86dQshNa8XMsONZK2W1ZaHbR9wi5EJy4XQ6KS0tFTHv5GitKS0tjSrJbIxOYa2YNonpkcezViq99SvDkpELyUafPn3Yv38/JSUlTQ/+lnC73QkJzrdFR40Lmheb0+mkT58+zZ6zUwi5J+BplrUSfrZnuac85hhB6KjY7XYGDhx4rMOIID8/n7Fjxx7rMKLoqHFB+8TWOawVf/OslXAhr/XVfiuxCYIgtDedQsgTsVYsynjLciScIAidhU4h5A2tlcYy8pxUo4FPrV8yckEQOgedQsjNjDxUR96IR56TZgi5ZOSCIHQWOoWQ+wK+Zlkrbr+b7qndARFyQRA6D51CyP1BPzaLrcmM3O13k5WShUKJkAuC0GnoFELuC/oihNwX9MUc5/a7cdlduOwuEXJBEDoNnULIzYzcbjV6GvsC8YXcaXOSak+V8kNBEDoNbSbkSimrUmq1UuqdtpqzufiDfuwWOzaLLXQ/Fm6/mxRrCqn2VGr8kpELgtA5aMuM/GZgcxvO12zMjNyqrEDj1orT5sRla9xa+WLfF8z858yI3iyCIAgdlTYRcqVUH2Au8HRbzJcovoDhkSulsFvszbJWGhPyx5c/zrI9y3hp/UvtFbIgCEKb0VYZ+Z+B24BgG82XEKa1AmC32mNm5EEdxBvwNssjNzPxfeVyEIUgCB2fVjfNUkqdBxzSWq9SSs1oZNwCYAFAbm5uxDl+iVBVVRX1XI/PQ1FhEfn5+aigYs/ePdFj6mrLC/cW4q50czR4NG4Mu4p3AbBi+wry7c2LM1ZcHQGJKzEkrsSQuBKnXWKLdf5bIj/A74H9wB7gIFAD/Kux57T1mZ3O+5z6tg9u01pr3e2BbvrH7/w4akxZTZkmD/2nL/+kL3j5Aj36b6PjvkafR/po8tBnLjqzVXF1BCSuxJC4EkPiSpzWxEZ7ndmptb5Da91Haz0AuBz4SGt9dWvnTQRzsRPiWytmRt6UR661Dh0JJ4dPCIKQDCR9HbnW2vDI62rI7RZ7zPJDt98NGELutDnjNs2q9deGmm5J1YogCMlAmx4sobXOB/Lbcs6mCOgAQCgjt1lsMTPycCFPsaY02bPcbrHLKUKCICQFSZ+Rm9l3hLUSo/ywoZA3dYpQ36y+VHor424uEgRB6Ch0PiG3xPbIw4XcYXU0mZGfkH4CIF0SBUHo+HQaIQ+vI2/KI0+xpcRtdWsKeW5absR9QRCEjkrSC7lpo0R45I1YKynWFFKsKQR1kEAwEDWuYUZe7a1ul7gFQRDaiqQX8pZaKxD7AIqGGXm1T4RcEISOTacR8lD5YXMWO+vO9oy14BkS8vSmrZUFby/g6je+1ZJ5QRCEKJJeyM3sOzwjj+WRm6KdYksJZeSxFjyba61UeCp46puneHH9i1R6Klv5LgRBEFpO0gt5Q2slXh25KdpNne3Z3MXOrYe3hm5vL9ve0vAFQRBaTacT8njWiinuDqujSWtFochJywHie+QF5QWh29IlURCEY0mnEfJQ+WGcxU4zI7db7E1aK2mONDIcGUB8a8XsxwKwv2J/K96BIAhC60h6IY9VfhjLIzfHNcdaSXekk+5ID92PRXFVceh2aW1pK96BIAhC62jTXivHguZaK6GM3Np4Rl7jqyHVnorL7gLiWyslNSV0c3XDG/AaPVmcrX8vgiAILaHTCHl498NY1or5mN1ib9QjN4+DsygLqfbUuBl5uaecLs4uhpC7RcgFQTh2JL2Qxyo/jJeRm+d6NmateAKe0PV0R3qj5YeZKZkEdED6lguCcExJeo88VvlhLI/cG/CGLJXGrBUzIweMAyj8sZtmVXgqyHJmke3M5qj7aKMxriteR1Afk+NMBUE4Duh0Qh7vhCBfwBeqbGnMWvH4PSEhd9qcoR2hDSl3l5OZkkkXZ5dG+5Yv2bKE0X8fzf2f3p/AuxIEQWg+rRZypVRfpdQypdRmpdRGpdTNbRFYc4lZfhjHWjEzcdM6iZeRm0LvsrniCrlprWS7shu1VpbuXArAJwWfNPctCYIgJERbeOR+4Bda62+UUhnAKqXUUq31pjaYu0kalh/GzciDvtCCaGNNs8KtFafNSa0v9pFwFZ4KslKycFgdjWbk5sahzYc3N/ctCYIgJERbHL5cpLX+pu52JbAZ6N3aeZtLSzzyRq2VsMXOeNaK1ro+I3dmU+2rjvnhAfVCXlRZJD65IAjtQpt65EqpAcBY4Ou2nLcxYpUfBnUwSjR9QV/Ci50uuyvmIc1uvxtf0EdmSiaZKZkA1ASiF0W11hQcLUChCOgAJdUlLX2bgiAIcWmz8kOlVDrwOnCL1roixvUFwAKA3Nxc8vPzW/Q6VVVVEc9dd3AdACuXr6TQVci+AqPvyYfLPsRhcYTGHTh4AG+tl/z8/FAlyqZtm8ivjYyjsqaS0kOl5OfnU3mkktLa0qhYy7xlABwsOBjy5ksro8dV+iqp9FZySuYpbKzYyNv5bzMkfUiL3ndLafj36ihIXIkhcSVGR40L2ie2NhFypZQdQ8Rf1Fq/EWuM1vpJ4EmACRMm6BkzZrTotfLz8wl/7q7Vu2ArTJs6jX5Z/Vjx+QrYA1OnTQ1tswfoUtSFKlsVM2bMMDLxz6HvgL7MOD0yjuBXQQb1HcSMGTPoW9aXogNFNIx1e+l2+BImjpyI1WKFbWBJsUSN23hoI3wBs06excblG+k9tDczTmzZ+24pDf9eHQWJKzEkrsToqHFB+8TWFlUrCngG2Ky1fqT1ISVGLI88/HETb8AbYb9AfI88ZK3YXDEXOyu9Rv/xjJSM0IeFOxDtpZvVLCN6jACgqKoo7vtYuHQhE5+aGPdQaEEQhHi0hUd+GnANMFMptabuZ04bzNsszKqV8MOXwx8PjQvzyJVSOKyOqKoVrXVE+WG8xc4an2HNpNpTSbOnAVAbiBZ8s5plWM4wAA5WHYz5Hqq91Tz4xYOsLFwZKlcUBEFoLq22VrTWnwGqDWJpEbHO7ASiqki8AW/oGhgLng2zX/N+UxuCwoXcnNMdjB5XVmt46T3Te+KyueKWKW44tCF0e1XRKuaeNDfmOEEQhFh0yp2d4Y+b+AL1GTkYm4IaWitmhh5efljrr0VrHTEuXMhNayVmRl5nrWS7suni7BJ3K394jbnUmwuCkCidTsjN3w2tlXCPHIxa8obWSvgBzWB45EEdjPpQiLBWHI1bKwpFVkqWIeSe2EJunjA0pc8UOW1IEISESXohD7WnbbCQGctaCc/IY1krZoYebq0AUfaKKeQum6vRjLystowsZxZWi7XRjLywspBurm4M7jqYfRXxhVxrzT/X/JO95XvjjhEE4fgj6YXczJatygo0vtgZ7pGnWONn5OGLnUDUpiCzkiV8sTNe1Uq2MxugUSEvqiqiV0Yv+mb2pbCykEAwEHPcqxtfZf5b87nolYtiXhcE4fikUwi5VVkxqiDrM/JY5YcRHrktJSojj7JW6k4JipeRp9pTcVgd2Cy2uB55tssQ8ixnVqMZuSnk/qCf4urimOPe2f4OAN8UfdNk61xBEI4fkl7IfQFfhPcd8sgbWCsNFzsdVkezFjshvpA7bU6UUqTZ02JWrRx1H6WLswsAXVIat1Z6ZvSkb1ZfIP5hzuHVLRsPbYw5RhCE44+kEvKCAli/PiviMX/QHxJviG+tNCw/bMxaaeiRN9wUZJ7raX4LSHekx8zIKz2VZDgygHprpWEFTFAHOVh1kJ7pPemR1gMgZk+WoA6yrXQbF51s2CobS+ILebW3mm+Kvol7XRCEzkVSCfmdd8Itt4xh1ar6x6KEPM5iZ/iGIGjeYqfLFt9aMa8BpDnSYnrkld5KMlLqhdwf9IeyeZNydzkBHSAnNYfuqd0BOFxzOGqu/RX7qfHVcPbgs0mzp7GpJH6X4HP+dQ7jnxzPki1L4o4RBKHzkFRC/thj4HAEefbZ+sfiZeSNbdGHuvJDf/MWOxsKea2/llR7auh+czNyMA5tDsfcNNTV1TUk5CU10Rn5ttJtAAztNpQTu53I9rLtUWMA9hzdw+f7Pgdg0ZpFMccIgtC5SCoh79YNRo4s5+OP6x9rWI2itCHqu/c27pGnWJte7IxXtWJaKybpjvT4GXkDIW/ok4dvGspwZOCwOmJm5AcqDgDQN6svfTP7hu435Kv9XwEw9oSxfLHviygrp+GcVd6quNcFQUgOkkrIAUaPPsrGjVBaatxvmJG/utgQ9VsX+qipczG01lGCH6vXSsPFzsaqVsKFPM2eRm0wUuy9AS/egDfCWoEYQl63bb+rqytKKbqndo8p5GbDrZ7pPemd0Tvugui64nXYLDauGnkVJTUlHPZGzwWw5uAaBj02iIlPTYzZPEwQhOQh6YT8pJOMzoPr1xv3Gwr5/75piHV5pY+3364fA0Rk5NqfQnmll0BYyXa8jLwpIY+VkVd66jokNpGRm9aKWW8eT8gPVh0kw5FBmiONPpl9KK0tjdkHZm3xWoZ1H8apfU4FYHtlbAvmHyv/gTfgZcvhLby/4/2YYwRBSA6STsj79zfS7M11LUnCz+Lctw/27DRup2X4eb9On0wLxRzn8cCH7zsoKvFw1131c8fb2RmrasXM1sFY7GzokZuWRZMZeZi1AvGFvKiqiJ4ZPQHok9kHIKa9srlkM6f0OIWRuSMBKKgpiBqjteb9ne8z58Q5pDvS+b8d/xc1xuSl9S8x5Zkpcni0IHRgkk7Ic3I8pKfXC3l4Rr5sGRA0bp8yysfSuo6wppCbGfnixXC01IGyeXn8caioO8+o4WJn1VFDrP/1iht/2Npp1GKnPUZG7o3MyLOcRtlkSzPyosoieqZHCnlDeyUQDLCvYh8DuwwkMyWTnNQcDrijxX5v+V72HN3DnCFzmDVwFh/u+jBqDBgVNQveXsBX+7/ihndvaNRv9wf9cXekCoLQviSdkCsFw4bFFvK1a8FhN7LuQUN8HDgARUVh/VjqPPJ//xsyUlNIzfRQVQVL6qr0Glord99u/P7oYzd/+Ut9DA2tlaojaVR53aGYIMxaqcvIs1IMIS93R1atHKk9gtPmDGX43V3dY1atxMrIGwp5UVUR/qCf/ln9ARjcdTCFtYVRc5mli6NPGM3EXhPZeWRnVFwA72x7h2pfNT8a/yM2lWyKm5Xvr9jPoEcH0e/P/dhZtjPmGEEQ2o+kE3KIL+Tr1sHJJxpi3XeAId6rVkVm5B4PfPghnDTYQQAvubmELBhzsdNhdbBjByx5zRDyIcNq+cMfwFtX5FLjqyHVZgj51q3w4nPpBC1epp3hY19dz6uGGbmn2olF21nyfjmHwxLuI+4jdHV1Dd3vntqdI7VHosoniyqLOCHtBAB6Z/YGooW84Khho/TL6gfA4OzGhXxY92GM7TnW+NsVr4sa99bWt+iV0YsHznoAm8XGe9vfixoD8MsPfsnBqoMUVRZx17K7Yo4B+KTgEy54+QJe2fBK3DGCICROmwi5UupcpdRWpdQOpdTtbTFnYwwbBgcOGJaIL1BfjbJuHZxysiHqPXv7sVhg5cr6XZ4Oq4NVqwyPfGA/o478rLM1H3wAgYCRkTusDizKwmuvAQEHCsXEqW4OHarP3MMz8oULwRIwGmdVe6tZuNAYE56RBwJwwQWKYE0Wy9eWM3Mm1NZZ6mW1ZSFbJRCApf/bHY1m9nlH6oWC/EQAACAASURBVD8UPJVU+6pDGXmKSifVksV7n+/nQJhzUlBuCHn/LnUZefZgSjwlVNVGVqVsPryZnNQcuqV2Y8wJYwCjiqUhX+3/iun9p5OZksm0ftP4z87/RI0pqS7hjc1v8JNJP+HWqbfyyoZXYp6EVFhZyPkvn8/b297mqjeuYlNF7A1Ntb5aFi5dyBWvX8HuI7tjjhEEIZK2OLPTCjwBfAcYDlyhlBre2nkbY+hQ4/e2bfUZeXExHDoEI4Ybom6x+Rg2DFasiFzs/NzYK8Og/g40mrPPCVBaCt98Yyx2mrbKm2/CxIkKp81Jr35ueveG5583nmsudu7eDf/7vzDrdKOV7Q9/XM3ixbBpU2RG/uST8Omn0CMrizPOKWf9erj7bmOu8Iz8rrvgyw9zAFixqYTZs6GkpP6IuJ7pPTlyBKZNg5qDvfnkm0KGDq3/gDEz8v5Z/dm9G175+2CCBOk5fA+33QZH6+z5zYc3MzxnOFrDymU9SfH34HfPruGmm+qrgYoqi9hXsY9JvSbh8UCP6lmsLV7LL+4qY+lS8NWV6b+7/V18QR/XjLqGeQOvRaP5+VNv8M47hq1l8uzqZ6nwVLD8Byvo6szhsXX/Ytky2LMHwq33O/57Bw9+8SBvbn6Tac+cwWcrj7JjR/3rgbFYe98n95H7cC4z/zmT7SV7OHiQULlpOMsPLOfl9S9TUn2YigrjA7QRqx+tjW9ejY0RhI5Gq496AyYBO7TWuwCUUouBeUD8PeStxBTyrVvrhXxdnTMweoQdvjB88YkT4d13wRuWkX/+OZx4InTLMhY0T5/hBWwsXQruUW5SrCns3w/Ll8P998MOixOPv5arroI//hEOFgdx+92k2lN58knDsz93ZhrvfQxXXVfFc4/BfffBlJsNIQ+607n7bjjjDKjonkl6RgULFsCf/gRXXmlk5AO6DOCNN+APf4A5P+nOe8ADfznMrd+Diy+GvOcMIc9QJ3DuubB6NYy4pBfWYYU4dhtjHnkE9g7ZSzdXN5b8O40f/xgCvYbAcBg3aycPPzyU556De+/VbC7fzMzcy5g1C5YtU6R8fwzVXVbz1OPw+ONw3nkw/qrlAKx+ZzJ950FJ6nS4Dv78xqc88rt55OTAZZfBiv7/IZ1cfn7VaD771ALXD+Pl3f/m5Z//GIDBg2HqaUHeHvgsPWpmcuGkCZQM+ikls+9k5hVroXg03brBqFHQfXABr/f6G90P/D9qPrmewiumcvo9t8M7f8dqhX79YOBAKB74KBv73k1a8Szyy1Zy0spz4Knl4MkiOxtyc6FrN82eE2+lcMAfjf8warPhlTdgzwwcDujSxfhJy66mYMjtHO3zCtqdBV/egl5xOjarlcxMyMiAlBRQ6cUcGfAc7rSt2A6PxbntaoLVXUOL4A6H8WNJO0J17lICQY21cAqqoh9+PwSDYLWCzQZ2O5Bagte1D6p7YKnqQzBofCPTGiyW+h9UkIClFl/1FGw2Yx7zQ6au1Q9K1d8Oxxyndf2Pifmc8B9LI2ldvA82r3cKdnvsa8fiw9B8TZ9vaty4jjW3396FGTPads62EPLeQPhpCPuByW0wbxTPrn6WN7e9yWtTZmCxwJYt4Bvow2Vz1Qv5yDohD/iYNAkWLYK9B4yM3GYxMvLzzquvYMnq6mHMmFSWLoWBpxgZuZnhXnQRPP6eC7ffzfXXwIMPwr8WGwuiDpXKX5+B88+H/j2NjNyRXs1PfmKMO+ESo/zwqScyKC01hPYX67Mo95Tzrwfgrbfg+uvhyJVHGOgcy//8ECZNgt/c2p33noNeQw6zaBFcfjncdf8h6A933pzL9m/g9dfh9WAvlu1expZ8uPpq+NnPoPvNBXis/bl6IZx2Gjz69BAmvALf/eEO/nwj3HIL3PDLQ3DrEV7/+zC67TCEe9eQ0Tyx4jH27vfzj7/ZePRReKf2a5hq418Pj+WC78APfjSJS1amcP1DnzDDO4+XXoInnwrivXkpbJ9DySELv/wlHBx2Cc/v+R1vfVTCttU5fPYZ/O+6fMoH7abr8vuYdQaMnfIjflX2W8665y+cr59m5UrYsAHerc4jqBVDC/MYc3ZfNqb+hI8nPMavz1uAf984du2CTYc2s6nX7WQfOp8zDr5F0PYp73SfxcjfXMf39OscOKAoKYFVrocoHPBHBh7+ESOC1/Bl9x9SPn8u8+3vkl0+g6NH4WB1Efk9z6cibTWDai6jNmMvRXN+woA5/2aO53lURT8qKoNsSv0ba3PuwG+tJMXfg4rBi7BNvJPh1TcwzvMLXIETKNFbWed6jJ0ZiwhY6r8a9PCeysnuKxni/y7+gI9dtnfY4XqJQ84vQmMyfUPo55lLP++5pOocKi17Kbat4JD9a0ocK/BZKrH7u5ATGEuOfxw5/rEobFRb9lNu3Um5dScV1h14VSWZgUF0CQwl2z8MR7ALAUstNZaD1FiKqLEWUWMpwoKdjMAA0gN9cQaMRm0eVU6tOkSttZhaSzEKK45gFinBbjiD3UkJdiVoceNV5XgsR/GqcnyWSixuGy57Dxw6E7tOx6pT8KlKvJYK/KqagPIQxIfCgsKKRduM39hQ2oZCEVBegvjQyo/SVhTW0HilrViwotFo5SdIAAgCFix1182xAJogqCDVtVW4XE5A1V1XoI3fGvOTUAONfdrUfTpqVXdP1b2GDj3XmEvXz9kMKtPmAWOaPb45tIWQxzp4OepdKaUWAAsAcnNzyc/PT/iF3t7+Nh8f+pgvv8znhBMm8+mnlZRll5Fpy+SDDw7StWs2Wzcb3sn2XduZbFkJTODN/90KCj7N383hw9C9+1YKdhk2xLJPl3HyyRN47bU+6Iv3EfQFefaFI/Tv7+DgwRXgg4IDBRzOzGfIkPH8/ZnD8F1Y9kEFJSVw+ulr2b7J2HTz6defMnlyLSkpp/LGuwex9rPyp4dtzJxZTEXFZnyVPorcRaxZk88NN3QnL28ElnOOcOjjLFLtXn7xi1UUbDYWMD9f/Tnn9+rKDTf04W8rD0J/2LW+G7/97ToyM8vw7/JTWFnIV19/xI03WujatR//1HuwlZ7IjTdu56KLCqkoCuKyuMhfn8+oIaPIy4PnP9nFIuDC0zL4/l2fkZ7up7DYjifg4d2vX+D00wcyYYKFn678iJrgiTz64mq6dzc+CE9OP5n3t7zDRePP58YbYfrl2/jphlJ+dv5ALhho/HvuqOrHot1B8g8+wHnjzmPcOKjYdD9fH0nnpbu7kWI1xp25cTr/PfQvfnTqPE4+OYM91Xv4/srnubTPd/nxzJ3ATs71z2bN8n/yWs21/GXWX5g5S/PTNT8lvcbB0/Pm09Vh9Grove+H/G3X35g2+CYunXoJ7xa9y7+3PczMnJn8avqlWJSXMu+9/Hztz3nefS63jL2F/rZM3tj+J7z+Ku4ffh9Tuk1Ba81bBW/x5P4n+WfqcM4eeDYbKzayrWobE7IncNOQm+iX2o+dVTt5ed/LLDv0CJszHyUnJYdCdyF2Zees3Fmc3/N8HBYHy8uW899D/+Wz6p/yGT8N/Xc8MG0g38/5Pv1S+3HIc4iVR1byzZG/skE/GhpjVVYGpw1mQuaZ9Ejpwb7KfRR4CthU9Rg+Xe8zpdvS6e3qzUnOAbhsLgprC9lX8x7bvM+HxtiUjW6ObnR1dGWAoxt+7eegew27PEupDlQD4LA4yLZnk+3Ipq89myBBqvzFVPp3UOgrp9JXicPiIM2WRrotnSxrGi6ri2rfUWr1PsoDbmoDtXiDXlxWF2m2NJxWJ6nKjk3Z0GgCOkCQoPFb1/+2W+w4lA2LshDUwfqxOkiAIEEdNORYWYyzB+rEOKADaDRBbYwBsChDrG1dIGip+wBAo7UOjVVKhURZxZSvekyBDv3WOvT80Bxh8zUHZ+aEFulf44Fq3aofYArwn7D7dwB3NPac8ePH65bw0/d+qtPuTdNaaz1njtajRmk97h/j9HkvnafHjNH6nHO0DgaDmjz03R/drT0erVNStL78tk80eeibH12qQett27R+etXTmjx0wdECvXSp8cVzyp8v0sMeG6GtVq3vuMN4zeFPDNfffeW7Wmut//hHrckq0OShe5//tB4+XOtgUOuv9n2lyUO/t+09rbXWd96pNd/5iWZhts7O1rqgwJjr2jev1f3+1C/0fv76d68x1xW/1evXG4/V+mo1eej7Pr4vNO7Hr/5Gqzyl9+73hR77y9d/0eShi6uKtdbG+077XZq++f9ujvibnfjwifo7//pO6P4Ty5/Q5KH3l+8PPbbu4DpNHvqldS9prbUOBAM68/eZ+oZ3boiY6+6P7taWeyz6aO1RrbXW9358ryYPfbDyYGhMMBjUgx4dpM/917laa63Lasp0yr0p+sZ3b4yY66m3n9LkoR/87EGttdYXLr5QZ9yfoUuqSyLGmf9Of13+V31P/j2aPPSi1YsixgSCAT3v5XmaPPRpz5ymyUOf88I52u1zR4w7VHVIn/r0qZo8NHnoQY8O0msPro0Ys2zZMr2zbKc+54VztOs+lx7515H62W+e1cFgUDdke+l2feO7N+pLX71UP/T5QxF/h3A2FG/QD33+kH7868f1pkObYo6p8lTpj3Z9pN/c/KZevn+5rvHWRMWltdZev1evL16vNxRv0IerD8ecS2uty93luqiySJfVlMWM3cTtc0f9nWIRbw4zro5GR41L69bFBqzUMTS1LTLyFcCJSqmBwAHgcuDKNpg3CpvFRhDjk/fkk40NQEMCPizY2LQJzj7b+HS0W+xGkywHjB0LG7f6YCxsWOugVy8YMgS+Xm945N6Al2nTICsLCg54sGU5CQQM3xmMmnKzvvz734f7n6yhFDiwO5U/32N4i+YBzOZuzt/+Ft6+p5Id/gw+/NDwdsGoJQ+v1/7u1Uf48cNwxy1dGTGi/vXSHekRm4KCrmK6pXajb+/6f65eGb0AoxqkR1oPymrLqPZVh2rIQ+NcvdhRtiN0f3PJZjIcGaHnAwztPhS7xc7a4rVcMfIKth7eSoWngsm9Ix2yM/qfwb2f3MsX+77gOyd+h3e3v8vEXhPJTc8NjVFKcfHJF/Po149SWlPK82ufxxPw8INxP4iYa0j6EM4adBa//+z3aDRLtizhdzN/F+oAaXLd2Ot4Yd0L/Pg9w3O/YsQVXDv62ogxFmVh8SWLufO/d/Lhrg/52ak/4/5Z94c2dpnkpOXw2XWfkb8nH7ffzexBs6PGAAzKHsT7VzfdtmBI1yE8PufxJsed0uMUTulxSqNj0hxpnDnwzCbnslvtjOgxoslxmSmZZKZkNjku1vuPhYplwgsdhlZXrWit/cBPgP8Am4FXtdbtcnyNzWIjoI3dg0OHGhUItR4/nho7Xq+xYAbGf+zmJqBJk2DbDsMaWPONnRkzDPE1PXKP34PTCVddBUWH3BTuTWHiRBg/3pjLZXOFhDwrCx54xPA/r7k8lUsuMcaYBzBX+4yvqVYrnDiikkF9Mhg3rj7+rJQsKjwVoa+BZsMsc3u+SU5qDodr64W8uLqY3LTciDGmEJvb9BuWHpr0cfVh99Hdobr0zYc3MyxnWMT/mA6rg+E5w1lbvBaArw98bfztek+KmGtK3ynYLDY+LviYkuoSvt7/NXNPnEtD/mfM/+AL+rj/0/v501d/4vR+p4fKHMN59FzDSlj44UIm9Z7EL6b8ImqMRVn4v6v+j7/O+SsvXPQCz1/0fExRcdqcPHLOI6y7YR2PnPNIqPqoIVaLlVmDZjH3pLnNFjFB6Oi0RUaO1vo9IPZukTYkXMhPPtl4rMbtp7JuW74p5DaLLSRc55wDj/3HEPLSYgez66xKs8OhWZp4++3w1O88+N1OHnigvgrAaXNGtLE9eUQtrDCE3MQ8gDm8JWylp/5QCZMsZxYaTbW3moyUjKjt+SbdU7tHnBJ0qPpQ6PQgk94Zxqagwkpjw0946WHEOFdv/EE/BUcLGNx1MJsPb+asQWfRkNEnjGbpzqVorfl6/9dkpmQytPvQiDGp9lSm9p3Kki1LGNhlIBrN3JOihXxEjxFcMvwSHvnqESzKwkvffSlqDMCwnGGsWrCKL/d/ybyh8+IKq8vu4oaJN8S8JghCGwn5t4XVYg0J+bhxRrlUjdvPkUobKSn1ZYmmtQIwezakZfqoBlwp9lAWHcrI63Zz9u0LI8a46ZaSxZlh33CdNmeosRVEHrxsEsrIvdWhx8J7kZuEtul7yslIyQjNG76zEwwLoLiq/gDm4upiJvSaEDHmhHRjl2dIyOsy8gFdBkSM6+00BH9H2Q66OLtQWFkY86v55N6TeX7t8+w6sovP933OxF4TsajoL2z/b8z/Y/5b8/nRuz9iXM9xjO85PmoMwDMXPMOY3DGM7TmWaf2mxRwDMDB7IAOzB8a9LghC0yTVFn1zK35QB0lPh5Ejodbr42CRjdNOM+p4IdJacTjgqmuMrPuWmxxk1Gmrmf2FHy7hDbrJTIv8Su6yuyK6H5pCHt790GlzolBNZuSmZ2n65OGnA4XTPTWy38qh6kP0SI3MyO1WOz3SenCgss5aOVpAmj0taq7ernohX3/I2O0zssdIGjJ70GwAnvrmKdYfWs93hnwnagzA5SMuZ3LvyThtTh6Y/UBc7zQzJZNfTf8Vc06cE/O6IAhtR1Jl5KaQ+4N+HFYH06fDWq8fT4mdWbPqx9kt9oheJadN9/HkEvjBdZFndgIRhyp4Ap6Q5WLisrkirJVYGblSqq4Uq4mM3FmfkUMjQu6q74Do9rup8FRELCia9M7oHcrI95TvoX+X/lHC2tXRlcyUTDYc2hDy5kfljoqa68SuJzI4ezAPfP4AAPNOnhc1BowPwC++/wXegDeuDy0IwrdLUmbkpkjPnw9Y/KBtXHFF5Ljww5dDW/TDTghq6JGDIZoNxclli52Rhws5gNPqjM7I41krDTJys1e5SU5aDjW+Gmp8NRyqPgQQ5ZGD0Txrb/lewMjIG/rjYHzITOo9iS/3f8mX+7/khPQTQrZMw3F5M/JQKK4dfS1Dug6JGmNiURYRcUHoQCSVkFuVsXvLFPJx4wz/+4rv2RgYZrParfUeOUQ2zTJp6JGDkZ03zMhT7akh8YYwa8XmihgXnpFrrY2MPMZiJ0Rm5F2cXbBarBHjzBK8wzWHQ155w6oVMLLonUd2EtRBCspjCznA1D5TWVu8ljc2v8HsQbPj2iFXj7qasoVlLJq3KOZ1QRA6Jkkl5GZGHnGAgcVPrxMimyrYLfbYGbk1LCOP4ZHHzMjtsa0Vs3bcxGmpz8g9AQ/+oD9uRl7hMU6yKKsti7JVwCg/BEPIG8vIT+p2EjW+GlYXraastoyTup0UNQbg7MFnh+I6d/C5MceYdHF2kZphQUgyktYjN2l4ZicY2XZERh5sJCNv6JHboj1yf9BvtMu12kNCHiX4VleoaqXhoRImoYy8zlpp2IvcJCIjrzYy8nhCDvD65tcBo4QwFlP7TuX68ddTWlvKpadcGnOMIAjJS6cUcrvVHlmN0gyPPKiDMRfwTC+81l8bEvJUe2pUaV64R97wUAmTNHsaVmWNsFZiZuRpRkZeXFXMvnKjH1n4TkyTod2MesuXN7wMxK5GAcP//vt5f495TRCE5CeprBXTSzaFXGujcU6UkDewVszsPNxaMTNyU8jN31FVK3VlhuaCZ7WvOmqhEyI98ngZuVKKzJTMiMXOWELeN7MvYJytWVBeQM/0njE3y/TK6EXvjN7sObqHPpl9Qh8AgiAcXySVkIc88rpNQaagx8rIw60Vb8CLVVkjsmhTGM3FTtNiiVW1AoR88obndZqEZ+SmB94wIwfDgzY3AoWfDhTxmnYXuWm57Dm6h4LygqhNPiZKKS4eZjSFuWTYJTHHCILQ+Ulqa8X8HW6ZmPfN/ihgeOTh/jhEZ+Tm+CiPvC4jN73xeEIe4ZHXWSuxmhblpOVQUlNCUAfjZuRg9EwpKC+g4GgBE3tPjDkG4P5Z9zP2hLHifQvCcUxSZuSmgJv2SazFzoYeebitYo6B+kzczMxjlR9CvbVS46sJ9VYJx2Vx1XvkcawVMCpSSqpLqPRUEtTBuEI+oMsAdh7Zyb6KfXHLCsFoD3Dd2OtCbQIEQTj+SCohb1hHnoi10jBrtygLNostKiNvylppzCOv9dcSCAbiLnZCfUYer2GWyaAug9h1ZBfegJdh3YfFHCMIggBJJuQN68hD1oq18TpyX8AXc7EwxZoS5ZG31FpxWp2h683JyIuqjJOJY+2yBCKaZI3rOS7mGEEQBEhSIU84Iw96ozxyiLRg4mXkMa0VRwxrxWoIfrWvOpSRx7I7clJz8AQ8bD28FTC22cdi9qDZZDuzObn7yc06SEAQhOOXpF7sNMW6OR55LCFPsaU06ZEnUrUCRk/ySk8lLpsrKi6orxFfc3ANUN9XvCFZziw23biJNHua7LQUBKFRWpWRK6UeUkptUUqtU0q9qZTq0vSzWk7DOvLGqlYabtGPm5EHm/DIY1krtmghT7Uaj1V6KqnwVMQ9Zsvcobn64GpSrClxFzvBsF1i2TOCIAjhtNZaWQqM0FqPArZhHLzcbjS7jtwSvdgZMyO3hmXkcTzyhtZKtTf2Yme6zbBRjriPUOaOX1Zobvb5cv+X9M7sLdm2IAitplVCrrX+oO7MToCvgD6tDyk+CXnkwcarVqB5Hnm4taK1juuRZ9iMzPlI7ZFG68PNzT3+oD+urSIIgpAIbbnY+f+A/2vD+aJoaR25LxC9IQjqPPJAEx55mLXiC/oI6EDMjDzTblgpZbVljQp5Rkr9Cfan5DR+srogCEJzaHKxUyn1IRCrRu5XWuu36sb8CvADLzYyzwJgAUBubi75+fkJB7upYhMA36z5Bsc+B1srjcqPLZu2kF9cP1/h/kK8fm/oNUrKSnBYHFGv6a5yc9B9kPz8fNYUGYuPq1eupshZFDHOpmxs3bmVD3wfGPMXFJIfiJxLuQ2LZOXGlRQdKSJXx3+PudZcCikkuzq7RX+HRKiqqmr312gJEldiSFyJ0VHjgvaJrUkh11rPbuy6Uup/gPOAWVpr3cg8TwJPAkyYMEHPmDEjsUiBjMIMWA2njDiFGUNn4NzvhG9g7OixzBhSP18++QT3BjnjjDOMY9i2u+ie2p2Gr5mzx6ggmTFjBptXbIZtMGPajKja7tSvUsnplcPYSWPhcxg1bBQzxkfOtWzZMhxWB9m9sqk+UM3wgcOjXs9kydglLN6wmFtOvSXmN4W2JD8/P24cxxKJKzEkrsToqHFB+8TW2qqVc4GFwAVa65qmxreWZnvkdX64ab00WrVSZ8GY1kqsI8xcNlfo6DUg5hZ9pRTZzmyKq4up8lY1Wo3SL6sft512W7uLuCAIxwet9cgfBzKApUqpNUqpdm16nUgdefj1xurIo5pmWaN3gKbaU6n118Y9r9Okq6srO8t2hm4LgiB8G7RqQ5DWOv4Jve1As+vI67bsewNe0khrNCNvaos+1B335mtayLNd2ew8IkIuCMK3S1Ju0W9OHTkkZq24/W7sFnvUyT9Qb62YB0fEE/Luqd0prCwEYh+WLAiC0B4kpZA3VX5oZuSmtRKrHzk02BAU47xOk1R7aoRHHk/Izc0+AH2z+sYcIwiC0NYktZAnkpE3Z0NQrIVOMA6IqPBUhLoaxtt+3y+rX+h2n8x23RslCIIQIqmEPF4/8niHRoSfxxk3Iw/zyGMtdILRwKrcU85R99HQ/ViEdymM96EgCILQ1iRl98OG/cibslaa5ZEH4mfkWSlZlLvDhDwltpCf3u90xvUcx/fHfj+h9yUIgtAaklLImyo/DLdWtNbNa2Prj++RZ6XUZ+QumyvuuIyUDFYtWNWCdyYIgtBykspaafbhy2EZuTkmXkbuC/oI6mCjHnmWMwt/0E9RVVFcW0UQBOFYkVRCHq+OPN6GIG/AG7JOYgm5WX3i9ruNqpV4HnmdlVJQXkAXZ7u2XBcEQUiYpBLyltSRN0fIa321uP3u+NZKXRa+t3yvCLkgCB2OpBTyROrIzTGxhNzsNW7WiMfqoQL1Gfne8r1xFzoFQRCOFUkl5BZlQaGaLD+MlZE3HAP1GXmNr4Zqb3XMAyMgstxQMnJBEDoaSSXkYIh5W3vkZkOseDs2w7NwycgFQehoJJ2QW5U1VEcet/wwzFppTMjDT/+p9lXHtVbMA5OBqF7lgiAIx5qkFPKGGbm549Mk0cVO0yOPl5HnpOWEbksPFUEQOhpJL+Q2iy3qJPrmZuSmcFd5q3D73XEz8vCOiP2z+rf+TQiCILQhSSfkFixRQt6Q5mbkZtXK4ZrDQPyuhgAXD7sYgKl9p7YiekEQhLYnqbboQ51HXldH7gv6Ygp5+GKn6aM3lpGbQh6vagXgxYtfpLSmtNExgiAIx4I2yciVUr9USmmlVPe2mK8xYlkrDUnUWimpLom4HwunzUnvzN6tC14QBKEdaLWQK6X6AmcBe1sfTtM0FPJYfcYTXew8XFuXkcfxyAVBEDoybZGR/wm4DdBtMFeTtDQjjyX4ZvlhczxyQRCEjkqrhFwpdQFwQGu9to3iaZLmeOSmaDe1IciiLKRYU5rlkQuCIHRUmlzsVEp9CMTaBfMr4E7g7Oa8kFJqAbAAIDc3l/z8/OZHGY6GooNF5Ofnc6DwAAFvIOZcFizs2L0DT7HRb/ybld9Q7CqOGudQDvaX7gdg87rNsKdlYVVVVbX8PbUjEldiSFyJIXElTrvEprVu0Q8wEjiEIX17AD+GT35CU88dP368bimDHhqkL1x8odZa68tfu1yf9JeTYo5LuTdF3/bBbfrpVU9r8tB7j+6NOa73H3vrjPszNHnoDcUbWhzXsmXLWvzc9kTiSgyJKzEkrsRpTWzASh1DU1tcfqi1Xg+E9q4rpfYAE7TWh1s6Z3NojkcOhk/uCzbe/RCMg5QPVB4AIN2R3g4RC4IgtC9JtyGoYa+VeEKeYjWOcWvMASnbRQAADOpJREFUI4fIbobZruw2jlYQBKH9abMNQVrrAW01V2M0p/wQjLpvt9/dbCG3KisZjox2iFgQBKF9ScqMvDnWitPmxB1ovpBnu7KjerYIgiAkA0kn5OG9VuKVH0J9Ru7xG1Ur8cZlO7MjfguCICQbSSfk4XXkjWXkLrsLt99Nrb8Wl80VN9s229LKyT+CICQrSSnkER55jCPcwMjIa3211PpqQzs4YzGgywBAhFwQhOQl6YW8KWulxlcTalcbi+n9p3PWoLP4+ZSft0u8giAI7U3StbENP7OzsfJDp81JaU2pYa00kpH3yujFB9d80C6xCoIgfBt02ozcZav3yKUZliAInZmkE3KbsoUOi2iqjrzWX+eRN2KtCIIgJDvJJ+QWW2jbfXM88qasFUEQhGQn+YRc2UKbfJpTRy4ZuSAInZ2kFPLmWCumR17jq5GMXBCETk3yCXmC1kpT5YeCIAjJTvIJeZi10pSQAxx1H5WqFUEQOjVJKeSmtdJUHTnAEfcRycgFQejUJJ+Q11krWmu8AS8ptpSY48J9cfHIBUHozCSfkCsjA/cH/XgD3rjtac2MHJCMXBCETk2rhVwpdZNSaqtSaqNS6sG2CKoxTCF3+90EdKB5Qi4ZuSAInZhW9VpRSp0JzANGaa09SqkeTT2ntZieeJW3Coh/YESaPS10WzJyQRA6M63NyG8A/qC19gBorQ+1PqTGMTPyal81EF/Is5xZodtpjrSYYwRBEDoDrRXyk4DTlVJfK6U+VkpNbIugGsMUcjMjT7HGXuwM7y/e1dW1vcMSBEE4ZjRprSilPgROiHHpV3XPzwZOBSYCryqlBmmtdYx5FgALAHJzc8nPz29RwAGvcTrQp19/CsDunbvJr4me66D7YOj2ni17yD/YstdrLlVVVS1+T+2JxJUYEldiSFyJ0y6xaa1b/AO8D8wIu78TyGnqeePHj9ct5c7Fd2ry0P/Z8R9NHnrR6kUxxx2pPaLJQ5OHXndwXYtfr7ksW7as3V+jJUhciSFxJYbElTitiQ1YqWNoamutlSXATACl1EmAAzjcyjkbpaG1Es8jz0zJDN0Wa0UQhM5Ma08IehZ4Vim1AfAC/1P3qdFumFUr1d7GFzstqv4zqmdGz/YMSRAE4ZjSKiHXWnuBq9solmbR3Iwc4Onzn6bcUx4h6oIgCJ2NpDuzMxEh//64738rMQmCIBxLki5VbW4duSAIwvFC8gl5g52d8ZpmCYIgHC8kn5Cr5i12CoIgHC8knZCbR7uJtSIIgmCQdEJuVVageYudgiAIxwNJJ+Sy2CkIghBJ0gm5aa1IRi4IgmCQdEIu1oogCEIkSSfkUrUiCIIQSdIJuWmtVHorgfj9yAVBEI4Xkk7IHRYjA6/0GEIuGbkgCMc7SSfkdmVHoaj112JRFqwW67EOSRAE4ZiSdEKulMJpcwL1NosgCMLxTNIJORAS8lR76jGORBAE4diTlELusrsAEXJBEARopZArpcYopb5SSq1RSq1USk1qq8Aaw2UzhDzNkfZtvJwgCEKHprUZ+YPAPVrrMcCv6+63O2KtCIIg1NNaIdeAecpxFlDYyvmahVgrgiAI9bT2qLdbgP8opR7G+FCY2vqQmsYUcBFyQRAEUE0deq+U+hA4IcalXwGzgI+11q8rpb4HLNBaz44zzwJgAUBubu74xYsXtyjgqqoqfr/n93xR+gWndTuN+0bc16J52pqqqirS09OPdRhRSFyJIXElhsSVOK2J7cwzz1yltZ4QdUFr3eIfoJz6DwMFVDTneePHj9ctZdmyZfqaN67R5KGvfP3KFs/T1ixbtuxYhxATiSsxJK7EkLgSpzWxASt1DE1trUdeCJxRd3smsL2V8zWLLs4uAPRI7fFtvJwgCEKHprUe+Q+BR5VSNsBNnXXS3piNsrqldvs2Xk4QBKFD0yoh11p/Boxvo1iajZmRW1RS7mcSBEFoU1qbkR8TfjblZxx1H2X+mPnHOhRBEIRjTlIKeao9lYfOfuhYhyEIgtAhEG9CEAQhyREhFwRBSHJEyAVBEJIcEXJBEIQkR4RcEAQhyREhFwRBSHJEyAVBEJIcEXJBEIQkp8k2tu3yokqVAAUtfHp34HAbhtNWSFyJIXElhsSVGB01LmhdbP211jkNHzwmQt4alFIrdax+vMcYiSsxJK7EkLgSo6PGBe0Tm1grgiAISY4IuSAIQpKTjEL+5LEOIA4SV2JIXIkhcSVGR40L2iG2pPPIBUEQhEiSMSMXBEEQwuiwQq6UOlcptVUptUMpdfv/b+/8Q60oojj++eKPXpakJqWZVEoEFZQvEfslgWElof0ijEjJIMSEooIEQSQI+kFBRRRUkkmUVGqPUNJMCgLN0ucvlJ6KlD/yQT/UKCjt9MfMteW6e99N7+7eyzsfuOzsnLNvDmdmz5ud2Z1JkZ8haUmUr5d0cQE2jZS0VtIOSdslPZqic5Okw5I6429+3nbFcvdK2hrL/DZFLkmvRH9tkdRegE2XJfzQKemIpMeqdArxl6SFkrolbUvkDZG0WlJXPA7OuHZG1OmSNKMAu16QtDPW0zJJgzKurVnnOdi1QNL+RF1Nzri25r2bg11LEjbtldSZcW2e/kqNDYW1sbQdmcv+AX2A3cAooD+wGbi8Smc28EZMTwOWFGDXcKA9pgcC36fYdRPwaQk+2wsMrSGfDKwEBIwH1pdQpz8R3oMt3F/ABKAd2JbIex6YG9NzgedSrhsC7InHwTE9OGe7JgF9Y/q5NLvqqfMc7FoAPFlHPde8dxttV5X8RWB+Cf5KjQ1FtbFm7ZGPA3aZ2R4z+wv4AJhapTMVWBTTHwETJSlPo8zsoJltjOmjwA5gRJ5lNpCpwLsWWAcMkjS8wPInArvN7FQ/BDstzOwr4Jeq7GQbWgTckXLpLcBqM/vFzH4FVgO35mmXma0ys2PxdB1wYaPKOx276qSeezcXu+L9fy/wfqPKq5casaGQNtasgXwE8GPifB8nB8wTOrHRHwbOLcQ6IA7ljAHWp4ivlbRZ0kpJVxRkkgGrJH0n6eEUeT0+zZNpZN9gZfgL4HwzOwjhRgTOS9Ep228zCU9SafRU53kwJw75LMwYJijTXzcCh8ysK0NeiL+qYkMhbaxZA3laz7r69Zp6dHJB0tnAx8BjZnakSryRMHxwFfAqsLwIm4DrzawduA14RNKEKnmZ/uoPTAE+TBGX5a96KdNv84BjwHsZKj3VeaN5HRgNXA0cJAxjVFOav4D7qN0bz91fPcSGzMtS8v6Xz5o1kO8DRibOLwQOZOlI6gucw6k9Cv4vJPUjVNR7Zra0Wm5mR8zs95heAfSTNDRvu8zsQDx2A8sIj7hJ6vFpXtwGbDSzQ9WCsvwVOVQZXorH7hSdUvwWJ7xuB+63OJBaTR113lDM7JCZHTezf4A3M8ory199gbuAJVk6efsrIzYU0saaNZBvAC6VdEnszU0DOqp0OoDK7O49wBdZDb5RxDG4t4EdZvZShs6wyli9pHEEH/+cs11nSRpYSRMmy7ZVqXUA0xUYDxyuPPIVQGZPqQx/JUi2oRnAJyk6nwGTJA2OQwmTYl5uSLoVeAqYYmZ/ZOjUU+eNtis5p3JnRnn13Lt5cDOw08z2pQnz9leN2FBMG8tjBrdBs8CTCTO/u4F5Me9pQuMGaCM8qu8CvgFGFWDTDYRHni1AZ/xNBmYBs6LOHGA7YbZ+HXBdAXaNiuVtjmVX/JW0S8Br0Z9bgbEF1eMAQmA+J5FXuL8I/0gOAn8TekAPEeZU1gBd8Tgk6o4F3kpcOzO2s13AgwXYtYswZlppY5W3sy4AVtSq85ztWhzbzhZCgBpebVc8P+nezdOumP9OpU0ldIv0V1ZsKKSN+ZedjuM4LU6zDq04juM4deKB3HEcp8XxQO44jtPieCB3HMdpcTyQO47jtDgeyJ1egaRBkmbXkJ8p6UtJfWrofJ61ep3jlIkHcqe3MIiwYmYWM4GlZna8hs7iHv6G45SCB3Knt/AsMDquRf1Civx+4ld3koZL+irqbpN0Y9TpIHyl6jhNhX8Q5PQK4op0n5rZlSmy/sAPZjYsnj8BtJnZM3GoZYCFpUmR1AWMN7OilhFwnB7pW7YBjtMEDAV+S5xvABbGRZCWm1lyx5luwqffHsidpsGHVhwH/iSs3QOc2LxgArAfWCxpekK3Leo7TtPggdzpLRwlbMF1EhZ2ZekjqQ1A0kVAt5m9SVjRrj3mCxhG2DLMcZoGD+ROryCOaX8dJy/TJjtXEVawg7CPaKekTcDdwMsx/xpgnf23DZvjNAU+2ek4gKQxwONm9kANnZeBDjNbU5xljtMz3iN3HMDMNgFra30QRNi53YO403R4j9xxHKfF8R654zhOi+OB3HEcp8XxQO44jtPieCB3HMdpcTyQO47jtDgeyB3HcVqcfwEqTRWzfXmUiQAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "from scipy.integrate import odeint\n", "import matplotlib.pyplot as plt\n", "import numpy as np\n", "\n", "#definir uma função para o sistema mecânico\n", "def mecanico(y, t, m, b, k, u):\n", " x1, x2 = y\n", " dydt = [x2, -(b/m)*x2 - (k/m)*x1 + (u/m)]\n", " return dydt\n", "\n", "m = 1\n", "b = 0.9\n", "k = 100\n", "u = 10\n", "\n", "y0 = [1.0, 0.0]\n", "t = np.linspace(0, 20, 1001)\n", "sol = odeint(mecanico, y0, t, args=(m, b, k, u))\n", "\n", "plt.plot(t, sol[:, 0], 'b', label='x1(t) - posição')\n", "plt.plot(t, sol[:, 1], 'g', label='x2(t) - velocidade')\n", "plt.legend(loc='best')\n", "plt.xlabel('t (s)')\n", "plt.grid()\n", "plt.show()" ] } ], "metadata": { "celltoolbar": "Raw Cell Format", "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.7.6" } }, "nbformat": 4, "nbformat_minor": 4 }