{ "cells": [ { "cell_type": "markdown", "metadata": { "nbpages": { "level": 0, "link": "[](https://jckantor.github.io/CBE32338/04.03-PID_Control.html)", "section": "" } }, "source": [ "\n", "*This notebook contains material from [CBE32338](https://jckantor.github.io/CBE32338);\n", "content is available [on Github](https://github.com/jckantor/CBE32338.git).*\n" ] }, { "cell_type": "markdown", "metadata": { "nbpages": { "level": 0, "link": "[](https://jckantor.github.io/CBE32338/04.03-PID_Control.html)", "section": "" } }, "source": [ "\n", "< [4.1 Relay Control](https://jckantor.github.io/CBE32338/04.01-Relay-Control.html) | [Contents](toc.html) | [4.10 Lab Assignment: PID Control](https://jckantor.github.io/CBE32338/04.10-Lab-Assignment-PID-Control.html) >

\"Open

\"Download\"" ] }, { "cell_type": "markdown", "metadata": { "nbpages": { "level": 1, "link": "[4.3 PID Control](https://jckantor.github.io/CBE32338/04.03-PID_Control.html#4.3-PID-Control)", "section": "4.3 PID Control" } }, "source": [ "# 4.3 PID Control\n", "\n", "## References\n", "\n", "Karl J. Astrom and Richard Murray (2016). _Feedback Systems: An Introduction for Scientists and Engineers, 2nd Edition_. [Web version available]( http://www.cds.caltech.edu/~murray/amwiki/index.php/Second_Edition)." ] }, { "cell_type": "markdown", "metadata": { "nbpages": { "level": 2, "link": "[4.3.1 PID Simulation](https://jckantor.github.io/CBE32338/04.03-PID_Control.html#4.3.1-PID-Simulation)", "section": "4.3.1 PID Simulation" } }, "source": [ "## 4.3.1 PID Simulation" ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "nbpages": { "level": 2, "link": "[4.3.1 PID Simulation](https://jckantor.github.io/CBE32338/04.03-PID_Control.html#4.3.1-PID-Simulation)", "section": "4.3.1 PID Simulation" } }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Collecting slycot\n", "Installing collected packages: slycot\n", "Successfully installed slycot-0.2.0\n", "Collecting control\n", " Using cached control-0.7.0-py2.py3-none-any.whl\n", "Requirement already satisfied: numpy in /Users/jeff/anaconda3/lib/python3.6/site-packages (from control)\n", "Requirement already satisfied: scipy in /Users/jeff/anaconda3/lib/python3.6/site-packages (from control)\n", "Requirement already satisfied: matplotlib in /Users/jeff/anaconda3/lib/python3.6/site-packages (from control)\n", "Requirement already satisfied: six>=1.10 in /Users/jeff/anaconda3/lib/python3.6/site-packages (from matplotlib->control)\n", "Requirement already satisfied: python-dateutil>=2.0 in /Users/jeff/anaconda3/lib/python3.6/site-packages (from matplotlib->control)\n", "Requirement already satisfied: pytz in /Users/jeff/anaconda3/lib/python3.6/site-packages (from matplotlib->control)\n", "Requirement already satisfied: cycler>=0.10 in /Users/jeff/anaconda3/lib/python3.6/site-packages (from matplotlib->control)\n", "Requirement already satisfied: pyparsing!=2.0.4,!=2.1.2,!=2.1.6,>=2.0.1 in /Users/jeff/anaconda3/lib/python3.6/site-packages (from matplotlib->control)\n", "Installing collected packages: control\n", "Successfully installed control-0.7.0\n" ] } ], "source": [ "!pip install slycot\n", "!pip install control" ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "nbpages": { "level": 2, "link": "[4.3.1 PID Simulation](https://jckantor.github.io/CBE32338/04.03-PID_Control.html#4.3.1-PID-Simulation)", "section": "4.3.1 PID Simulation" } }, "outputs": [], "source": [ "%matplotlib inline\n", "\n", "import numpy as np\n", "import matplotlib.pyplot as plt\n", "import control.matlab as control\n", "\n", "# control constants\n", "Kc = 0.85\n", "tauI = 10000\n", "\n", "# control transfer function\n", "Gc = Kc*control.tf([tauI,1],[tauI,0])\n", "\n", "# model transfer functions\n", "Gp = control.tf([0.31],[16,1])*control.tf([1],[135,1])" ] }, { "cell_type": "code", "execution_count": 5, "metadata": { "nbpages": { "level": 2, "link": "[4.3.1 PID Simulation](https://jckantor.github.io/CBE32338/04.03-PID_Control.html#4.3.1-PID-Simulation)", "section": "4.3.1 PID Simulation" } }, "outputs": [ { "data": { "text/plain": [ "[]" ] }, "execution_count": 5, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXQAAAD8CAYAAABn919SAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4wLCBo\ndHRwOi8vbWF0cGxvdGxpYi5vcmcvpW3flQAAIABJREFUeJzt3Xl4VfW97/H3NzsjJCQMASIEwqw4\nABoBhzrgRO259Xoee489rcdWLZ2OVWvbc2zPvdb2PM89rT0dT2vrU21txbZW0Xqo1lorDrcKDTLP\noDJDghCSEDLt/b1/7BVAGsiGZLOy1/68nmc/e6/fWgu+K0s/Wfz2b62fuTsiIpL5csIuQEREeocC\nXUQkIhToIiIRoUAXEYkIBbqISEQo0EVEIkKBLiISEQp0EZGIUKCLiEREbncbmFkh8ApQEGz/hLvf\na2ZXAPeT/KXQBHzM3Tce788aMmSIV1VV9bhoEZFssnjx4j3uXt7ddt0GOtAKzHL3JjPLA14zs+eA\nB4Dr3H2NmX0G+DfgY8f7g6qqqqipqUnhrxQRkU5mtjmV7boNdE8+7KUpWMwLXh68BgTtpcCOEy9T\nRER6SypX6JhZDFgMjAd+6O4Lzew24FkzOwg0ADPTV6aIiHQnpS9F3T3u7lOBkcB0MzsLuAu41t1H\nAj8Dvt3VvmY2x8xqzKymrq6ut+oWEZGjnNAoF3evBxYA7wemuPvCYNVvgAuPsc+D7l7t7tXl5d32\n6YuIyEnqNtDNrNzMyoLPRcCVwBqg1MwmBptdFbSJiEhIUulDrwAeCfrRc4DH3X2+mX0CeNLMEsA+\n4JY01ikiIt1IZZTLcmBaF+1PAU+loygRETlxKY1yERHpCXcn4dCRSBBP+KFXxxGfD73cSQTv8YST\nSEDCk8vuTjxYTiSSf2bC/fCrc50f/js717mDc3gbD9Y5h5e9cznR2Z5sSx5Dcv/k++F9Oo+vs/3o\n7Tp3vv7ckYwZ0j+tP2cFukiExBPOwfY4B9vitLR3vhK0dBz+3NoRp7U9QVs8QWt7PHhPLrd1JGjt\nSNAefG6Pd7Y7HYnkcnuH0x587oh78j3hR31OvCewOxKau/jc0QMV6CJR5u40t8XZf7CdhpZ2Gg52\n0NTaTmNLB02tHTQF740tHTS3dXCgLU5za/De1kFza5zmtngyxNvjtHUkTroWM8iP5ZCfm3P4PTeH\nvFjylR8z8mI55MaM4rzc5Oecw225OTnkxezQ59wcIxaz5HvncvA68nMsx4jZe5dz7Oh3yAm2yzEj\nJ4f3rLOgPWaGGe/ZpnO90bkcbJOTbOvc3gyMv90eI1h3uL1zn86fm2FdbmPBNqeKAl2kl8QTTn1z\nG+8eaGNPUyvvNrWx90Ab+5rbqG9up765jX3Be/3BdhoOttPQ0kG8m6tXMyjOz6V/QS79CmL0z8+l\nX36MoSWFFA2O0T8/RlFejML8GP3ycinKz0kuv+eVQ0Hue9/zc5Ofk+/JwD3VASS9S4Eu0o14wqlt\nbGFHfQu1DS3sbmihtrGV3Q2t1Da2UNvQyp6mVvY1t3GsbC4pyKWsfx4D++VT1i+f0YP7U1qUx4Ci\n3OR7YR4DgvfiwlyKC3IpCd775ccUtJISBbpkvY54gp37W9iyt5nN7zazdV8zO+oPBq9kgB/dB5yb\nYwwtKaB8QCGjBvfjvKqBDOmfz+DiAgYX5zOofz5DigsY1D+f0qI88mJ6UrWknwJdsoK7s7uhlU11\nTclXbRNv7TnAlr3NbN938D2BnRczKkqLqCgtZMaYQZxWVkRFWSGnlRYxbEAhwwYUMLBfPjk5umqW\nvkWBLpHzblMra3Y2snZXA6t3NrCxNhngB9rih7YpLshlzJD+nD2ilA+cXcHowf2oHNSP0YP7M3xA\nITGFtWQgBbpktF37W1i6tZ5l2+pZvaOBNTsbqG1sPbR+aEkBE4eV8KHqSsaV92dceTHjhhYztKRA\n/dISOQp0yRgH2+Is2bovGeBb61m2dT+7GlqAZJ/2hGElXDxhCJMrBnD68AGcUVHC4OKCkKsWOXUU\n6NJnHWjtYPHmfSx8+10WvrWXZdvqaY8n+7qrBvdjxthBTK0sY0plGZMrBlCYFwu5YpFwKdClz4gn\nnGXb6lmwro5XN9SxYtt+OhJOLMc4e0Qpt1w8hhljBjGtciAD++eHXa5In6NAl1C929TKKxvqWLCu\njlfW17GvuZ0cgymVZcy5ZCwzxw7mvNED6V+g/1RFuqP/S+SU27W/hedW7uTZFTup2bwPdxjcP5/L\nJw3lstOHcsmEIZT10xW4yIlSoMsp0VWInz68hDuumMDlk4Zy9ohSjesW6SEFuqRNS3uc51bu5Dd/\n3crCt/ceCvG7rpzItWdXMH5ocdglikSKAl163crt+/n1X7fwu6U7aGzpYPTgftx5xUQ+cI5CXCSd\nug10MysEXgEKgu2fcPd7LXlXxr8DHwLiwAPu/v10Fit9V3NbB08u3savFm1l9c4GCnJzuPbsCv7h\n/EpmjBmkm3hEToFUrtBbgVnu3mRmecBrZvYccAZQCZzu7gkzG5rOQqVvqm1o4ZHX3+HRN7aw/2A7\nZ542gK9fdyYfnDqC0qK8sMsTySqpzCnqQFOwmBe8HPg08I/ungi2q01XkdL3rNvVyE9ffYvfLd1B\neyLBNZOH84lLxnDe6EFhlyaStVLqQzezGLAYGA/80N0Xmtk44B/M7HqgDvicu2/oYt85wByAUaNG\n9VrhEo4lW/bx3T9t4OX1dRTm5XDj9EpuuWgMVWmeWktEupdSoLt7HJhqZmXAU2Z2Fsk+9RZ3rzaz\nvwceBt7Xxb4PAg8CVFdXa2LBDLVuVyP/+cd1/HH1bgb1z+cLV0/kIzNG645NkT7khEa5uHu9mS0A\nZgPbgCeDVU8BP+vd0qQv2PJuM9/503qeXrqd4vxc7r5qIh+/eAzFunNTpM9JZZRLOdAehHkRcCXw\nDeBpYBbJK/NLgfXpLFROrbrGVr734np+vWgruTFjziVj+fSl43QHp0gflsplVgXwSNCPngM87u7z\nzew1YK6Z3UXyS9Pb0linnCLxhPPYoi3c/4e1NLfFuXF6JbfPmsCwAYVhlyYi3UhllMtyYFoX7fXA\nB9JRlIRjxbb9/NvTK1i2bT8XjhvM1647SzcCiWQQdYQK+w+2859/XMcv39jMkOICvnfjVD445TTd\nDCSSYRToWe73y3dy7zOr2HuglZsvqOLzV09kQKFuCBLJRAr0LNXY0s5Xn1nNk29u45yRpfz84+dz\n1ojSsMsSkR5QoGehxZv3cddvlrJtXzOfu2ICt88aT14sJ+yyRKSHFOhZpCOe4L9e2sgP/ryRitJC\nHv/kBVRX6VZ9kahQoGeJrXubuePXS3hzSz3XTxvBfdedqb5ykYhRoGeBVzfU8dm5b+LA926cynVT\nR4RdkoikgQI9wtydR/7yDl///RrGlxfz05urqRzUL+yyRCRNFOgR1daR4N5nVvKrRVu58oxhfPfG\nqXr+ikjE6f/wCNp7oI1PPbqYRW/v5TOXjeMLV0/SBMwiWUCBHjFrdzVw2yM11Da2qr9cJMso0CNk\n4VvvcusjNfTLj/H4Jy9gamVZ2CWJyCmkQI+IVzfU8Ylf1DCirIhHb5tBRWlR2CWJyCmmQI+AP63e\nzWfmvsm4ocX88tbpDCkuCLskEQmB7vfOcL9fvpNPPbqYMypK+NUnZijMRbKYrtAz2JOLt/HFJ5Zx\n3uiBPPyx8ynRnZ8iWa3bK3QzKzSzRWa2zMxWmdl9R63/gZk1pa9E6crchZu5+7fLuGDcYB65ZbrC\nXERSukJvBWa5e5OZ5QGvmdlz7v6GmVUDGkpxij1es5WvPLWSWacP5UcfOZfCvFjYJYlIH9DtFbon\ndV6B5wUvD+YYvR/4Uhrrk6P8ee1u7pm3gvdNGMKPP3qewlxEDknpS1Ezi5nZUqAWeMHdFwL/DDzj\n7jvTWaActnRrPZ+du4TJFQN44KPnkZ+r77RF5LCUvhR19zgw1czKgKfM7BLgQ8Bl3e1rZnOAOQCj\nRo06+Uqz3Ft1Tdzy879SXlLAwx87X89lEZG/cUKXeO5eDywALgfGAxvN7B2gn5ltPMY+D7p7tbtX\nl5eX97Dc7FTb2MLNP1sEwCO3TKe8REMTReRvpTLKpTy4MsfMioArgcXuPtzdq9y9Cmh29/HpLTU7\nNbV2cMvP/8qexjYe/tj5jBnSP+ySRKSPSuXf7RXAI8GXoDnA4+4+P71lCSQfgfvpRxezZmcjP725\nWs9mEZHj6jbQ3X05MK2bbYp7rSI55L7/XsWrG/Zw/w3ncPmkoWGXIyJ9nIZJ9FFPLdnG3IVb+NSl\n4/hQdWXY5YhIBlCg90HrdjVyz7wVzBgziC9cPTHsckQkQyjQ+5jGlnY+/ehiSgrz+ME/TiM3plMk\nIqnRYOY+xN35lyeXs3lvM4/dNoOhJYVhlyQiGUSXf33Iw//vHZ5dsYsvXTOJGWMHh12OiGQYBXof\nUfPOXv7vs2u4evIw5lwyNuxyRCQDKdD7gD1NrXz2sTcZMbCI+z80BTMLuyQRyUDqQw+Zu3P348uo\nb25n3mfOp7RIzzUXkZOjK/SQ/XbxNl5eX8c97z+dM08rDbscEclgCvQQ7drfwtfnr2Z61SD+6YKq\nsMsRkQynQA+Ju/Plp1bQHk/wzRvOISdH/eYi0jMK9JA8tWQ7f15byxeunkSVnqAoIr1AgR6C2oYW\n7vvv1Zw3eiAfv2hM2OWISEQo0E8xd+crT6+kpT3ON284h5i6WkSklyjQT7Fnlu3ghdW7+fxVExlX\nrqcOi0jvUaCfQnWNrXz1mVVMrSzjtvfpblAR6V0K9FPoq8+s4kBrnPvV1SIiaZDKnKKFZrbIzJaZ\n2Sozuy9on2tm68xspZk9bGa6xfE4/rJxD79fsZPbZ41nwrCSsMsRkQhK5Qq9FZjl7lOAqcBsM5sJ\nzAVOB84GioDb0lZlhosnnK/NX82IsiI+oQdviUiadBvontQULOYFL3f3Z4N1DiwCRqaxzoz225qt\nrN3VyD3Xnk5hXizsckQkolLqQzezmJktBWqBF9x94RHr8oCbgD8cY985ZlZjZjV1dXW9UXNGaWxp\n51t/XEf16IF84OyKsMsRkQhLKdDdPe7uU0lehU83s7OOWP0j4BV3f/UY+z7o7tXuXl1eXt7zijPM\njxZsYk9TG//77ybrsbgiklYnNMrF3euBBcBsADO7FygHPt/rlUXA1r3NPPTq2/z9tBFMqSwLuxwR\nibhURrmUm1lZ8LkIuBJYa2a3AdcAH3b3RHrLzEz/8dxaYjnGF2dPCrsUEckCqUxwUQE8YmYxkr8A\nHnf3+WbWAWwGXg+6Eua5+9fSV2pmWfT2Xn6/Yid3XjmBitKisMsRkSzQbaC7+3JgWhftmu3oGBIJ\n5+vzV1NRWsgnLxkXdjkikiV0p2gazFuynRXb9/Ol2ZMoytcwRRE5NRTovexgW5z7n1/LlMoyrpsy\nIuxyRCSLKNB72WOLtrC7oZWvXHuGZiESkVNKgd6LWtrj/OTlTcwcO4jpYwaFXY6IZBkFei96vGYr\ntY2tfG7WhLBLEZEspEDvJW0dCX68YBPVowdywbjBYZcjIllIgd5LnnxzGzv2t3D7FRN0i7+IhEKB\n3gva4wl+tGAjU0aWcsmEIWGXIyJZSoHeC55esp2tew9y+yxdnYtIeBToPRRPOD9asInJFQO44oyh\nYZcjIllMgd5D85fv4O09B7h91nhdnYtIqBToPZBIOP/1541MHFbMNWcOD7scEclyCvQe+MOqXWyo\nbeKfZ03QXaEiEjoF+klKJJzvv7iBseX9NbWciPQJCvSTtGB9LWt3NfLZy8YT09W5iPQBCvST9PO/\nbGbYgAI+OPW0sEsREQFSm4Ku0MwWmdkyM1tlZvcF7WPMbKGZbTCz35hZfvrL7RveqmvilfV1fGTG\naPJi+p0oIn1DKmnUCsxy9ynAVGC2mc0EvgF8x90nAPuAW9NXZt/yyzc2kxczbpxeGXYpIiKHdBvo\nntQULOYFLwdmAU8E7Y8A/zMtFfYxB1o7eGLxNt5/VgVDSwrDLkdE5JCU+gvMLGZmS4Fa4AVgE1Dv\n7h3BJtuArJie5+ml22ls6eDmC0eHXYqIyHukFOjuHnf3qcBIYDpwRlebdbWvmc0xsxozq6mrqzv5\nSvsAd+cXf9nMmacN4NxRA8MuR0TkPU7oGz13rwcWADOBMjPLDVaNBHYcY58H3b3a3avLy8t7Umvo\nFr69l3W7G7n5gird5i8ifU4qo1zKzaws+FwEXAmsAV4Cbgg2uxn4XbqK7Ct++fpmSovy+B9TNFRR\nRPqe3O43oQJ4xMxiJH8BPO7u881sNfBrM/t3YAnwUBrrDN2u/S38YdUubr14DEX5sbDLERH5G90G\nursvB6Z10f4Wyf70rPDYws0k3PnoDH0ZKiJ9k+6KSUFrR5zHFm1h1qShjBrcL+xyRES6pEBPwR9W\n7mJPUxs3XaCrcxHpuxToKfjF65upGtyPSyZk9igdEYk2BXo3Vm7fz+LN+7jpgio981xE+jQFejce\nW7SFwrwcbjhvZNiliIgclwL9OFra48xftoPZZw6ntCgv7HJERI5LgX4cL62tpaGlg+vP1dW5iPR9\nCvTjmLdkO+UlBVw0bnDYpYiIdEuBfgz7DrSxYF0t1005jVxNYiEiGUBJdQzzl++gPe5cf25WPBVY\nRCJAgX4M85Zs5/ThJUyuGBB2KSIiKVGgd+HtPQdYsqWe66eN0GNyRSRjKNC78NSS7ZjBdVPV3SIi\nmUOBfhR35+kl27lo3BCGl2rOUBHJHAr0oyzevI8te5u5fpquzkUksyjQjzJvyXaK8mLMPmt42KWI\niJwQBfoRDt3qf9Zw+hekMpmTiEjfkcqcopVm9pKZrTGzVWZ2R9A+1czeMLOlZlZjZhk/e9GhW/3V\n3SIiGSiVy9AO4G53f9PMSoDFZvYC8E3gPnd/zsyuDZYvS1+p6TdvyXaGlhRw0fghYZciInLCur1C\nd/ed7v5m8LkRWAOMABzovOumFNiRriJPhUO3+k89jZieey4iGeiEOorNrIrkhNELgTuB583sWyR/\nMVx4jH3mAHMARo0a1YNS0+vQrf7T9GRFEclMKX8pambFwJPAne7eAHwauMvdK4G7gIe62s/dH3T3\nanevLi/vu1O4/W7pjuSt/qfpVn8RyUwpBbqZ5ZEM87nuPi9ovhno/PxbIGO/FK1tbGHxln1ce3ZF\n2KWIiJy0VEa5GMmr7zXu/u0jVu0ALg0+zwI29H55p8YLq3fjDtecqbHnIpK5UulDvwi4CVhhZkuD\nti8DnwC+Z2a5QAtBP3kmen7VbqoG92PisOKwSxEROWndBrq7vwYca9jHeb1bzqnX0NLO65v2cMtF\nY/RkRRHJaFl/p+hLa2tpjztXq7tFRDJc1gf686t2MbSkgGmVZWGXIiLSI1kd6C3tcV5aW8fVZw4j\nRzcTiUiGy+pAf3XDHg62xzW6RUQiIasD/flVuxhQmMvMsYPDLkVEpMeyNtA74gleXLObK84YRl4s\na38MIhIhWZtki97Zy77mdq45c1jYpYiI9IqsDfQ/rtpNYV4Ol0zsu8+XERE5EVkZ6O7O86t2ccmE\ncvrla2YiEYmGrAz05dv2s3N/i0a3iEikZGWgP79qF7Ec44ozhoZdiohIr8naQJ85dhBl/fLDLkVE\npNdkXaBvrG1iU90BZqu7RUQiJusC/flVuwC4arICXUSiJSsDfWplGcNLC8MuRUSkV2VVoNc2trB8\n236umqybiUQkelKZgq7SzF4yszVmtsrM7jhi3e1mti5o/2Z6S+25l9fVAXD5JI1uEZHoSeWumg7g\nbnd/08xKgMVm9gIwDLgOOMfdW82sz6fkgvV1DBtQwBkVJWGXIiLS61KZgm4nsDP43Ghma4ARJOcU\n/Q93bw3W1aaz0J7qiCd4dX0ds88arqnmRCSSTqgP3cyqgGnAQmAi8D4zW2hmL5vZ+b1fXu9ZurWe\nhpYOLlN3i4hEVMoPMjGzYuBJ4E53bzCzXGAgMBM4H3jczMa6ux+13xxgDsCoUaN6rfATtWBdHbEc\n46LxQ0KrQUQknVK6QjezPJJhPtfd5wXN24B5nrQISAB/k5bu/qC7V7t7dXl5eE82XLC+lvNGDaS0\nKC+0GkRE0imVUS4GPASscfdvH7HqaWBWsM1EIB/Yk44ie6q2sYWV2xu4dJIelSsi0ZVKl8tFwE3A\nCjNbGrR9GXgYeNjMVgJtwM1Hd7f0Fa+sT/6euUyBLiIRlsool9eAYw0L+WjvlpMeC9bVMrSkgMkV\nA8IuRUQkbSJ/p2hHPMGrG/Zw6cRyDVcUkUiLfKAv21bP/oPtGq4oIpEX+UDvHK548QQNVxSRaMuK\nQD93VJmGK4pI5EU60GsbW1ixfb+6W0QkK0Q60DuHK146UcMVRST6Ih3oC9bVUl5SwJmnabiiiERf\nZANdwxVFJNtENtAPD1dUd4uIZIfIBvqCdXXkGLxvvAJdRLJDpAP93FEDKe2n4Yoikh0iGeh1ja3B\ncEVdnYtI9ohkoL+yPjkZtMafi0g2iWSgL1hfx5BiPV1RRLJL5AI9nnBeWV/HpRPLycnRcEURyR6R\nC/SlW/ex/2A7l5+u/nMRyS6RC3QNVxSRbJXKnKKVZvaSma0xs1VmdsdR679gZm5mfeL5tBquKCLZ\nKpUr9A7gbnc/A5gJfNbMJkMy7IGrgC3pKzF1Gq4oItms20B3953u/mbwuRFYA4wIVn8H+BLQJyaH\nflnDFUUki51QH7qZVQHTgIVm9kFgu7sv62afOWZWY2Y1dXV1J11oKjqfrqjhiiKSjVIOdDMrBp4E\n7iTZDfMV4P90t5+7P+ju1e5eXV6evq6QI5+uqOGKIpKNUgp0M8sjGeZz3X0eMA4YAywzs3eAkcCb\nZjY8XYV2R09XFJFsl9vdBpZ8mPhDwBp3/zaAu68Ahh6xzTtAtbvvSVOd3dJwRRHJdqlcoV8E3ATM\nMrOlwevaNNd1wl5aV8t5ozVcUUSyV7dX6O7+GnDcTml3r+qtgk5GbWMLK7c38MVrJoVZhohIqCJx\np6gmgxYRiUigazJoEZEIBHpHPHHo6YqaDFpEslnGB/rSrfU0tHRwue4OFZEsl/GBvmBdHbEc4+IJ\nfeLZYCIiocn8QF9fy7mjyigt0nBFEcluGR3oncMV9TAuEZEMD/Rnl+8E0O3+IiJkcKC3tMf58ctv\nUT16oJ6uKCJCBgf6b/66lV0NLdx11UQNVxQRIUMDvaU9zo8WbGR61SAuHDc47HJERPqEjAz0Xy3a\nwu6GVu68aoKuzkVEAhkX6Mmr803MGDOIC8dp7LmISKeMC/S5C7dQ19jKXVdNDLsUEZE+JaMC/WBb\nnAcWbOLCcYOZOVZ95yIiR8qoQH/0jc3sadLVuYhIV7oNdDOrNLOXzGyNma0yszuC9vvNbK2ZLTez\np8ysLJ2FNrd18OOXN3Hx+CGcXzUonX+ViEhGSuUKvQO4293PAGYCnzWzycALwFnufg6wHrgnfWXC\nL1/fzLsH2rjrqgnp/GtERDJWt4Hu7jvd/c3gcyOwBhjh7n90945gszeAkekq8kBrBz955S3eN2EI\n543W1bmISFdOqA/dzKqAacDCo1bdAjzXOyX9rV+8vpm9B9rUdy4ichwpB7qZFQNPAne6e8MR7V8h\n2S0z9xj7zTGzGjOrqaurO6kiy0sK+F/VIzl31MCT2l9EJBuYu3e/kVkeMB943t2/fUT7zcCngCvc\nvbm7P6e6utpramp6UK6ISPYxs8XuXt3ddrkp/EEGPASsOSrMZwP/AlyaSpiLiEh6dRvowEXATcAK\nM1satH0Z+D5QALwQPE/lDXf/VFqqFBGRbnUb6O7+GtDVE7Ce7f1yRETkZGXUnaIiInJsCnQRkYhQ\noIuIRIQCXUQkIhToIiIRkdKNRb32l5nVAZtPcvchwJ5eLCcT6Jizg445O/TkmEe7e3l3G53SQO8J\nM6tJ5U6pKNExZwcdc3Y4FcesLhcRkYhQoIuIREQmBfqDYRcQAh1zdtAxZ4e0H3PG9KGLiMjxZdIV\nuoiIHEdGBLqZzTazdWa20cz+Nex6esNxJt8eZGYvmNmG4H1g0G5m9v3gZ7DczM4N9whOnpnFzGyJ\nmc0PlseY2cLgmH9jZvlBe0GwvDFYXxVm3SfLzMrM7IlgUvU1ZnZB1M+zmd0V/He90sx+ZWaFUTvP\nZvawmdWa2coj2k74vJrZzcH2G4I5Jk5anw90M4sBPwTeD0wGPhxMUp3pjjX59r8CL7r7BODFYBmS\nxz8heM0BHjj1JfeaO0jOTdvpG8B3gmPeB9watN8K7HP38cB3gu0y0feAP7j76cAUksce2fNsZiOA\nzwHV7n4WEANuJHrn+efA7KPaTui8mtkg4F5gBjAduLfzl8BJcfc+/QIuIDlTUufyPcA9YdeVhuP8\nHXAVsA6oCNoqgHXB558AHz5i+0PbZdKL5GTiLwKzSM6CZSRvtsg9+nwDzwMXBJ9zg+0s7GM4weMd\nALx9dN1RPs/ACGArMCg4b/OBa6J4noEqYOXJnlfgw8BPjmh/z3Yn+urzV+gc/o+j07agLTKOmnx7\nmLvvBAjehwabReXn8F3gS0AiWB4M1Lt7R7B85HEdOuZg/f5g+0wyFqgDfhZ0M/3UzPoT4fPs7tuB\nbwFbgJ0kz9tion2eO53oee3V850Jgd7V5BqRGZpzrMm3u9q0i7aM+jmY2d8Bte6++MjmLjb1FNZl\nilzgXOABd58GHODwP8O7kvHHHHQZXAeMAU4D+pPscjhalM5zd451jL167JkQ6NuAyiOWRwI7Qqql\nVwWTbz8JzHX3eUHzbjOrCNZXALVBexR+DhcBHzSzd4Bfk+x2+S5QZmads2cdeVyHjjlYXwrsPZUF\n94JtwDZ3XxgsP0Ey4KN8nq8E3nb3OndvB+YBFxLt89zpRM9rr57vTAj0vwITgm/I80l+ufJMyDX1\nmFnXk2+TPLbOb7pvJtm33tn+T8G35TOB/Z3/tMsU7n6Pu4909yqS5/HP7v4R4CXghmCzo4+582dx\nQ7B9Rl25ufsuYKuZTQqargBWE+HzTLKrZaaZ9Qv+O+885sie5yOc6Hl9HrjazAYG/7K5Omg7OWF/\nqZDiFw/XAuuBTcBXwq6nl47BQGKJAAAArUlEQVTpYpL/tFoOLA1e15LsO3wR2BC8Dwq2N5KjfTYB\nK0iOIAj9OHpw/JcB84PPY4FFwEbgt0BB0F4YLG8M1o8Nu+6TPNapQE1wrp8GBkb9PAP3AWuBlcAv\nSU4oH6nzDPyK5HcE7SSvtG89mfMK3BIc+0bg4z2pSXeKiohERCZ0uYiISAoU6CIiEaFAFxGJCAW6\niEhEKNBFRCJCgS4iEhEKdBGRiFCgi4hExP8He4mQAxeTyt8AAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "t = np.linspace(0,1000)\n", "y,t = control.step(Gp,t)\n", "plt.plot(t,50*y + 22)" ] }, { "cell_type": "code", "execution_count": 7, "metadata": { "nbpages": { "level": 2, "link": "[4.3.1 PID Simulation](https://jckantor.github.io/CBE32338/04.03-PID_Control.html#4.3.1-PID-Simulation)", "section": "4.3.1 PID Simulation" } }, "outputs": [ { "data": { "text/plain": [ "[]" ] }, "execution_count": 7, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAD8CAYAAACMwORRAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4wLCBo\ndHRwOi8vbWF0cGxvdGxpYi5vcmcvpW3flQAAIABJREFUeJzt3Xl4XHd97/H3d0b7ZknW4k3e4lV2\nHNtR7JCEAFmdhQTC5hCWAq25z1PKltvecKHQ0nK5lFu2h5TiQuE2N81CkoIJjg0YlxDSeHe8r7It\nydqtzZasbeZ3/5ixoyhyNJZHOpozn9fzTGbOmV9G36Pj56Mzv/M752fOOURExF8CXhcgIiLxp3AX\nEfEhhbuIiA8p3EVEfEjhLiLiQwp3EREfUriLiPiQwl1ExIcU7iIiPpTi1Q8uKipyM2fO9OrHi4gk\npB07djQ754qHa+dZuM+cOZPt27d79eNFRBKSmZ2KpZ26ZUREfEjhLiLiQwp3EREfUriLiPiQwl1E\nxIcU7iIiPqRwFxHxIc/GuUvy6e4LcaK5k8azPZzvDdHdF+J8X4iu6Ou+UJjstBSy01PITg+Sm5FC\ndloKuRmpTJ6QQX5WKmbm9WaIJASFu8RdKOzYU9PGgboOjjd2crzpHJXN56hpPc+VTNmbnRZkakEm\nU/MzmVaQxbSCTK4qzmH+pFym5mcSCCj4RS5QuEtctHf18fujTWw+1MjvjzTR0tkLQGZqkNnF2Swr\nK+A9y6dxVXEOU/IzyExNITMtSGZq5JGRFiAlEKCzt5/OnsjjXE+Izp5+2s/3Udt2ntNt56lpPc/p\n1vPsrGqj/XzfxZ+fnRZk3qRcFkzKZX5pLoumTmDxlAlkpgW9+pWIeErhLiPW1tXLz7bX8JsDDeyo\naiUUdhRkpfL2+SW8Y0EJ184oYHJexmUdUedlpJKXkRpT247uPo42nONw/VkO13dwqP4sL+yr54mt\n1QAEA8a80lyWluWztGwC15TlM7ckl6CO8CUJmLuS78lXoKKiwuneMompuqWLH790gqe3V9PVG6J8\nch63LIgE+tKyfE/D0zlH49ke9p1uZ3d1G7ur23i1uo2O7n4gcoS/fEYBFTMKuW5mAUun55OVpmMc\nSRxmtsM5VzFcO/2rlpjtO93OD1+sZP3eOgy4b+kU1tw8mwWT8rwu7SIzozQvg9K8DG5dWApEAv9E\ncye7q9vYVdXGtpMtfGfTEZyLHN0vnpJHxcxCVswqZMXMQgqy0zzeCpErpyN3Gdb+2na+9quDvHz8\nDDnpKXxw5XQ+duNMJk/I9Lq0EWs/38euqla2n2xl28kWdle30dMfBmB+aS4rZxeyctZErptVQElu\nhsfVirwm1iN3hbtcUndfiO9uOsraFyspyEplzc2zWb1iesx94omkpz/Enpp2tp5o4ZXKM+w41UpX\nbwiA2cXZrJw1kZWzClk5uzCh/6hJ4lO4yxXZUnmGLzy3l8rmTt537TS+eM9C8rOSp7uiLxRmf20H\nWyrPsOVEC9tOtHC2J9JvP70wi5Wzot04swqZXpil8fcyZhTuMiJnu/v43y8c4vEtVZQVZvL1dy/h\nprlFXpfluVDYcbCugy0nWthSeYatJ1to64oMxSzJTb8Y9NfNLGR+aa7G3MuoUbjLZfuv42f43FO7\naTzbzcdvnMXn75inkSSXEA47jjWdY+uJFraeaGHbyRbq2rsByM1IYdn0Aq6dXsC1MyIjcnLS9XuU\n+FC4y2V5ens1//O5vcyYmMU/vn8pS8vyvS4poTjnqGk9z9YTLeyoamXnqVYON5zFOQgYzJ+Ux/Lp\n+VxTls/SsnyuKs7ReHsZEYW7xCQcdnzrN0f4/uZjvHVuEY8+tNyXJ0y90NHdx+6qNnacamVnVSu7\nqto41/PaePvFUyewtCyfJdPyWTQlj+mFWerOkWFpnLsMq7svxF8+s4dfvlrLgyvK+Or9i0kN6kah\n8ZKXkcrN84q5eV5kovpw2FHZ3Mmr1W28WtPGqzXt/OSPJ+kNRYZg5qSnUD45j/IpkceiKXnMKckh\nPUW3UJDLp3BPUmfO9bDmsR3sONXKI3ct4JM3z9aIj1EWCBhzSnKYU5LDe66dBkSGYB6uP8uB2g72\n13awv7adp7ZVc74vMgwzGDBmTMxiXkku80pzmFuay7zSXGYWZSn05U0NG+5m9q/AvUCjc27xEO8b\n8F3gbqAL+BPn3M54FyrxU9l0jo/9dBv17d3800PLufvqyV6XlLTSU4IsmRbpmrkgFHacPNPJgdoO\njjScvfj49YF6wtFe1IDB1IJMZhXlMGtiFrOKsplZlM3MidlMzs9Q8EtMR+4/Bb4P/Nsl3r8LmBt9\nrAR+EH2Wcai6pYvVa18hFHY8seZ6lk8v8LokGSQYMK4qzuGq4pzXre/uC3G86RxHG85R2dzJieZO\nTjZ3svNU68W+fAAzKM3NYFpBZvSRxdSCTCblZVCSl86kvAwKs9P0Tc3nhg1359yLZjbzTZrcD/yb\ni5yZfcXM8s1ssnOuLk41Spw0n+vhwz/eQk9/mJ/9t7cwrzTX65LkMmSkBlk0ZQKLpkx43XrnHM3n\nejnR3ElVSxc1rV3UtJ6nprWL7ada+eWeOkLh1w+cSAsGKM5NpzQvnaKcdCbmpFGYncbE7Mjridnp\n5GdF7tA5ITOV3IwUnexNMPHoc58KVA9Yromue0O4m9kaYA3A9OnT4/CjJVZnu/v4k59spb6jm8f/\n9HoFu4+YGcW56RRHL6YarD8UpuFsD/Xt3TR2dFPf0U1DRw8NHd00dHRT1dLFzqo2Wrt63/BH4LWf\nAbnpKUzISiUnPZWc9CBZaSnkpKeQlRYkO/qckRokIzUQeU4Jkp4aID0lSFqKkRIIkBoMvO51MGAE\nAxAwI2BGMGAEAkbAwLCLP9sg8h8HDnAOHC76HDlZ7RyEnCPsHOGwi7wOR7q5Qs5Fngc8+sPh6LOj\nPzRgORRdFw5fXH+xTShMXzj6HHqtzetfRz7nQruBn9EX/YxP3zqXd14zZbT+SQDxCfeh/pwP+S/E\nObcWWAuRoZBx+NkSg57+EJ98bAcH687yo49UcO0MdcUkk5RggKn5kRms3kw47Gg/38eZzl7OnOuh\n/XzfxUdHdz8d0ddnuyOTqbR19VLT2kVnT4jO3n66ekOX/OPgF2aQGgiQEjRSAkZKMEBKwEgNvrbu\nwh+tlGCA1ICREjSyU1Mi6wIBUoPGhMzRH24cj3CvAcoGLE8DauPwuRIHobDjc0/t5uXjZ/jW+6/h\nHQtKvC5JxqlAwCjITqMgO405JTnD/w9D6A+F6e4P090Xij7C9PSHLh7R9oZeO7qNHOFGj7bDkSPu\nUPT5wt+IwUfoxoUjebt4RG8W+fYSNCMw+FuAQTAQIBiIPkfbBC0avkG7GLopQYuuj4b1gPWR4A4k\n1IVn8Qj3dcCnzOxJIidS29XfPj445/jKun2s31vPl+5ZyAPLp3ldkvhcSjBATjCg2y2MA7EMhXwC\neDtQZGY1wFeAVADn3D8D64kMgzxGZCjkx0arWLk839t0jP/3ShWffNts/vSts70uR0TGUCyjZR4c\n5n0H/HncKpK4+N2hBr792yO8Z/k0Hlm1wOtyRGSM6VpzH6pv7+bhp19l4eQ8vvbuxRrPLJKEFO4+\nEwo7PvPkLnr6w3z/g8vISNWViiLJSGc9fOZ7m46y5UQL//i+a95whaOIJA8dufvIy8eb+d7vjvLA\n8qkXb0wlIslJ4e4Tzed6+OyTu5lVlM3f3f+G+7uJSJJRuPtAOOx4+OlXaTvfx/cfXE62xhiLJD2F\nuw/8yx8q+f2RJv76noWUT8nzuhwRGQcU7gnuYF0H39x4mLsWT+JD18/wuhwRGScU7gksHHZ88T/2\nkpeZyv9699Uazy4iFyncE9jT26vZWdXGF+5aQEF2mtfliMg4onBPUGfO9fD1Fw6xYlYh79WwRxEZ\nROGeoL7+wiE6e/r5+3fp9gIi8kYK9wS09UQLz+yo4U/fOlszKonIkBTuCaa3P8yXfr6XqfmZfPrW\nOV6XIyLjlK52STA/fukERxrO8aOPVJCVpt0nIkPTkXsCqWnt4nubjnJ7eSm3lZd6XY6IjGMK9wTy\nN+sORJ7vW+RxJSIy3incE8TmQ4389mADn7lt7rCz2IuIKNwTQDjs+MaGQ8yYmMXHb5zldTkikgAU\n7gng+b11HKo/y+dvn0dainaZiAxPSTHO9YXCfOvXh1kwKZd3LpnidTkikiAU7uPcsztqOHmmi4fv\nmE8goCtRRSQ2CvdxrLsvxHc3HWVpWT63LSzxuhwRSSAK93Hs8S1V1LV381d3ztf9Y0Tksijcx6lz\nPf380+Zj3DhnIjfMKfK6HBFJMAr3ceonL53gTGcv//2O+V6XIiIJSOE+DrV19bL2xUpuLy9l2fQC\nr8sRkQQUU7ib2SozO2xmx8zskSHen25mm81sl5ntMbO7419q8vjn31dyrrefh++Y53UpIpKghg13\nMwsCjwJ3AeXAg2ZWPqjZl4CnnXPLgNXAP8W70GTR2NHNT18+wf3XTGHBpDyvyxGRBBXLkfsK4Jhz\nrtI51ws8Cdw/qI0DLiTRBKA2fiUmlx/8/jj9Icdnb9NRu4iMXCzhPhWoHrBcE1030N8AHzKzGmA9\n8BdDfZCZrTGz7Wa2vampaQTl+ltbVy9PbavmvqVTmFmU7XU5IpLAYgn3oQZYu0HLDwI/dc5NA+4G\nHjOzN3y2c26tc67COVdRXFx8+dX63ONbqujqDfFnb53tdSkikuBiCfcaoGzA8jTe2O3yCeBpAOfc\nfwEZgAZnX4buvhA/+eNJbp5XzMLJ6msXkSsTS7hvA+aa2SwzSyNywnTdoDZVwK0AZraQSLir3+Uy\n/GL3aZrP9fDJm3XULiJXbthwd871A58CNgIHiYyK2W9mXzWz+6LNHgb+zMxeBZ4A/sQ5N7jrRi4h\nHHasfbGS8sl53HDVRK/LEREfiGmGZefceiInSgeu+/KA1weAG+NbWvLYfLiR402dfHf1Ut1DRkTi\nQleojgM/fLGSKRMyuPvqyV6XIiI+oXD32O7qNraeaOHjN80iNajdISLxoTTx2L+8WEluRgqrV0z3\nuhQR8RGFu4eqznTxwr46Hlo5g5z0mE5/iIjEROHuoR+9VEkwYHzsxplelyIiPqNw90hrZy9Pb6/m\n/qVTKc3L8LocEfEZhbtHHnvlFN19YdbooiURGQUKdw/0hcI89sop3javmHmluV6XIyI+pHD3wKaD\nDTSd7eHD18/wuhQR8SmFuwce31LFlAkZvGNBideliIhPKdzH2MnmTv5wtJnVK6YTDOhWAyIyOhTu\nY+yJbVUEA8YHrisbvrGIyAgp3MdQT3+In22v4baFJRr+KCKjSuE+hjbsq6els5eHVupEqoiMLoX7\nGPr3LVVML8zipjmapEpERpfCfYwcazzLlhMtfHDldAI6kSoio0zhPkYe31JFatB477XTvC5FRJKA\nwn0MdPeFeHZHDasWT6YoJ93rckQkCSjcx8Dze+ro6O7noZW6Z7uIjA2F+xh4fMsprirOZuWsQq9L\nEZEkoXAfZQdqO9hV1cYHV87Q5NciMmYU7qPs37eeIj0lwHuWT/W6FBFJIgr3UXS+N8TPd9Vyz5LJ\n5GeleV2OiCQRhfso+vWBes719POBCt1HRkTGlsJ9FD278zTTCjK5bqZOpIrI2FK4j5KGjm5eOtrE\nA8um6opUERlzCvdR8vNdpwk7ePdyXZEqImMvpnA3s1VmdtjMjpnZI5do834zO2Bm+83s3+NbZmJx\nzvHszhqWT89nVlG21+WISBIaNtzNLAg8CtwFlAMPmln5oDZzgS8ANzrnFgGfHYVaE8b+2g6ONJzj\nAR21i4hHYjlyXwEcc85VOud6gSeB+we1+TPgUedcK4BzrjG+ZSaW53aeJi0Y4N4lk70uRUSSVCzh\nPhWoHrBcE1030Dxgnpn90cxeMbNV8Sow0fSFwqx79TS3LizR2HYR8UxKDG2GGurhhvicucDbgWnA\nH8xssXOu7XUfZLYGWAMwfbo/b6L1h6NNNJ/rVZeMiHgqliP3GmDgVTjTgNoh2vzCOdfnnDsBHCYS\n9q/jnFvrnKtwzlUUFxePtOZx7dmdpynMTuNt8/y5fSKSGGIJ923AXDObZWZpwGpg3aA2PwfeAWBm\nRUS6aSrjWWgiaD/fx28ONHDfNVNIS9EoUxHxzrAJ5JzrBz4FbAQOAk875/ab2VfN7L5os43AGTM7\nAGwG/tI5d2a0ih6v1u+to7c/zAO6SZiIeCyWPnecc+uB9YPWfXnAawd8PvpIWs/trGFOSQ5XT53g\ndSkikuTUdxAnp850su1kKw8sn6r7touI5xTucfIfu05jBu9aqi4ZEfGewj0OnHM8t/M0N1w1kSn5\nmV6XIyKicI+HnVWtVLV08e5lGtsuIuODwj0OfvlqHWkpAe5cVOp1KSIigML9ioXCjl/treOW+SXk\nZqR6XY6ICKBwv2JbT7TQdLaHe6/RTcJEZPxQuF+h5/fUkpka5JYFJV6XIiJykcL9CvSHwrywr55b\nF5aQlRbT9WAiImNC4X4FXj5+hpbOXt55zRSvSxEReR2F+xV4fk8tuekpugOkiIw7CvcR6u0Ps2Ff\nPbeXl5KRGvS6HBGR11G4j9AfjjbR0d2vUTIiMi4p3Efo+T11TMhM5aY56pIRkfFH4T4C3X0hfnOg\ngVWLJmlSDhEZl5RMI/Cfh5s416MuGREZvxTuI/DLPbVMzE7jLbMnel2KiMiQFO6Xqau3n98dbOSu\nqyeREtSvT0TGJ6XTZdp0sJHzfSHuXaILl0Rk/FK4X6bn99RSkpvOdTMLvS5FROSSFO6X4Wx3H5sP\nN3HPkskEA5onVUTGL4X7ZfjtwQZ6+8PqkhGRcU/hfhnW761n8oQMlpXle12KiMibUrjHqLOnnxeP\nNHHnokkE1CUjIuOcwj1Gmw830tMf5q7Fk7wuRURkWAr3GL2wr56inDQqNEpGRBKAwj0G3X0hNh9q\n5I5FkzRKRkQSQkzhbmarzOywmR0zs0fepN17zcyZWUX8SvTei0ea6OoNqUtGRBLGsOFuZkHgUeAu\noBx40MzKh2iXC3wa2BLvIr22YV89EzJTuV73khGRBBHLkfsK4JhzrtI51ws8Cdw/RLu/A/4B6I5j\nfZ7r7Q/zm4MN3F5eSqruJSMiCSKWtJoKVA9Yromuu8jMlgFlzrnn41jbuPDy8WbOdverS0ZEEkos\n4T7UGUR38U2zAPBt4OFhP8hsjZltN7PtTU1NsVfpoQ376slJT+GmuUVelyIiErNYwr0GKBuwPA2o\nHbCcCywG/tPMTgLXA+uGOqnqnFvrnKtwzlUUF4//6en6Q2F+faCBWxaUkJ6iSbBFJHHEEu7bgLlm\nNsvM0oDVwLoLbzrn2p1zRc65mc65mcArwH3Oue2jUvEY2nqyhZbOXnXJiEjCGTbcnXP9wKeAjcBB\n4Gnn3H4z+6qZ3TfaBXppw756MlIDvG3++P+WISIyUEosjZxz64H1g9Z9+RJt337lZXkvHHZs2FfP\nO+aXkJUW069JRGTc0Ni+S9hZ1Urj2R5WqUtGRBKQwv0SXthXT1owwC0LSrwuRUTksinch+BcpEvm\nrXOLyM1I9bocEZHLpnAfwt7T7ZxuO68uGRFJWAr3IWzYV09KwLi9vNTrUkRERkThPsiFLpnrZ08k\nPyvN63JEREZE4T7IscZzVDZ3cqe6ZEQkgSncB9m4vx4zuFNdMiKSwBTug2zYX8+ysnxK8jK8LkVE\nZMQU7gNUt3Sx73SHRsmISMJTuA/w6wMNANy5SOEuIolN4T7Axn31LJiUy4yJ2V6XIiJyRRTuUU1n\ne9h2qkVdMiLiCwr3qN8ebMA5dcmIiD8o3KM27KtnxsQsFkzK9boUEZErpnAHOrr7ePl4M3cumoTZ\nUFPGiogkFoU7sPlQI30hpy4ZEfENhTuRq1JLctNZVpbvdSkiInGR9OHe3Rdi86Em7lhUSiCgLhkR\n8YekD/cXjzRxvi/EqkWTvS5FRCRukj7cN+5vYEJmKitnF3pdiohI3CR1uPeFwvz2YAO3LiwhNZjU\nvwoR8ZmkTrStJ1poP9/HKo2SERGfSepw37CvnszUIDfPK/a6FBGRuEracA+HHRv31/O2ecVkpAa9\nLkdEJK6SNty3n2ql8WwPdy/RKBkR8Z+kDff1e+tITwlwy4ISr0sREYm7mMLdzFaZ2WEzO2Zmjwzx\n/ufN7ICZ7TGzTWY2I/6lxk847HhhXx1vm1dMTnqK1+WIiMTdsOFuZkHgUeAuoBx40MzKBzXbBVQ4\n55YAzwD/EO9C42lXdSsNHT3coy4ZEfGpWI7cVwDHnHOVzrle4Eng/oENnHObnXNd0cVXgGnxLTO+\nfrWnnjR1yYiIj8US7lOB6gHLNdF1l/IJ4IUrKWo0XeiSuXluMbkZqV6XIyIyKmIJ96HupuWGbGj2\nIaAC+OYl3l9jZtvNbHtTU1PsVcbR7po26tq7uWeJLlwSEf+KJdxrgLIBy9OA2sGNzOw24IvAfc65\nnqE+yDm31jlX4ZyrKC725sKh9XvqSAsGuHVhqSc/X0RkLMQS7tuAuWY2y8zSgNXAuoENzGwZ8EMi\nwd4Y/zLjwznHC/vqeevcIvLUJSMiPjZsuDvn+oFPARuBg8DTzrn9ZvZVM7sv2uybQA7wMzPbbWbr\nLvFxnnq1pp3Tbee5+2qNkhERf4tpkLdzbj2wftC6Lw94fVuc6xoV6/fWkRo0bitXl4yI+FvSXKHq\nnGP93jpumlPEhEx1yYiIvyVNuO893U5N63nuUpeMiCSBpAn3X+2tIyVg3KEuGRFJAkkR7s45Xthb\nz41zisjPSvO6HBGRUZcU4b6/toOqli7uvloXLolIckiKcF+/t45gwLijXOEuIsnB9+F+YZTMDVdN\npCBbXTIikhx8H+67qts4eaaLe3V7XxFJIr4P96e2VpOVFuSeJVO8LkVEZMz4OtzP9fTzyz21vHPJ\nFM24JCJJxdfh/qs9tXT1hnj/dWXDNxYR8RFfh/uT26qZW5LD8un5XpciIjKmfBvuRxrOsquqjQ9c\nV4bZUPONiIj4l2/D/alt1aQGjXcve7MZAUVE/MmX4d7TH+I/dp3m9vJSJuake12OiMiY82W4//ZA\nIy2dvXzguulelyIi4glfhvuT26qYMiGDm+YUeV2KiIgnfBfuNa1dvHSsmfdVlBEM6ESqiCQn34X7\nz7bXAPC+imkeVyIi4h1fhXso7HhmRw03zSliWkGW1+WIiHjGV+H+0rFmTredZ7VOpIpIkvNVuD+9\nrZqCrFRuKy/xuhQREU/5Jtwbz3bz6wP1PLB8GukpQa/LERHxlC/CvT8U5nNP7cbM+OBKdcmIiPgi\n3L+x4RB/PHaGr71rMVcV53hdjoiI5xI+3H+x+zT/8ocTfPQtM3hfhW7tKyICCR7u+2vb+R/P7mHF\nzEK+dG+51+WIiIwbMYW7ma0ys8NmdszMHhni/XQzeyr6/hYzmxnvQgdr7ezlk4/tID8zjUcfWk5q\nMKH/TomIxNWwiWhmQeBR4C6gHHjQzAYfJn8CaHXOzQG+DXwj3oUO1B8K8xdP7KKxo4d//vC1FOfq\nzo8iIgPFcri7AjjmnKt0zvUCTwL3D2pzP/B/o6+fAW61UZwh45sbD/PSsWb+/l2LWVqmWZZERAaL\nJdynAtUDlmui64Zs45zrB9qBifEocLBfvlrLD1+s5MPXz9DcqCIilxBLuA91BO5G0AYzW2Nm281s\ne1NTUyz1vcHE7DTuKC/lr3UCVUTkklJiaFMDDDxEngbUXqJNjZmlABOAlsEf5JxbC6wFqKioeEP4\nx+KGOUXcoPu0i4i8qViO3LcBc81slpmlAauBdYParAM+Gn39XuB3zrkRhbeIiFy5YY/cnXP9ZvYp\nYCMQBP7VObffzL4KbHfOrQN+DDxmZseIHLGvHs2iRUTkzcXSLYNzbj2wftC6Lw943Q28L76liYjI\nSOnKHxERH1K4i4j4kMJdRMSHFO4iIj6kcBcR8SHzaji6mTUBp0b4vxcBzXEsJxFom5ODtjk5XMk2\nz3DOFQ/XyLNwvxJmtt05V+F1HWNJ25wctM3JYSy2Wd0yIiI+pHAXEfGhRA33tV4X4AFtc3LQNieH\nUd/mhOxzFxGRN5eoR+4iIvImEi7ch5usOxGZWZmZbTazg2a238w+E11faGa/MbOj0eeC6Hozs+9F\nfwd7zGy5t1swcmYWNLNdZvZ8dHlWdJL1o9FJ19Oi68d8EvbRYGb5ZvaMmR2K7u+3+H0/m9nnov+u\n95nZE2aW4bf9bGb/amaNZrZvwLrL3q9m9tFo+6Nm9tGhflasEircY5ysOxH1Aw875xYC1wN/Ht2u\nR4BNzrm5wKboMkS2f270sQb4wdiXHDefAQ4OWP4G8O3oNrcSmXwdxngS9lH0XWCDc24BcA2Rbfft\nfjazqcCngQrn3GIitw1fjf/280+BVYPWXdZ+NbNC4CvASiJzV3/lwh+EEXHOJcwDeAuwccDyF4Av\neF3XKGznL4DbgcPA5Oi6ycDh6OsfAg8OaH+xXSI9iMzqtQm4BXieyHSNzUDK4P1NZD6Bt0Rfp0Tb\nmdfbcJnbmwecGFy3n/czr82vXBjdb88Dd/pxPwMzgX0j3a/Ag8APB6x/XbvLfSTUkTuxTdad0KJf\nQ5cBW4BS51wdQPS5JNrML7+H7wB/BYSjyxOBNheZZB1ev11jNgn7KJoNNAE/iXZF/cjMsvHxfnbO\nnQb+D1AF1BHZbzvw936+4HL3a1z3d6KFe0wTcScqM8sBngU+65zreLOmQ6xLqN+Dmd0LNDrndgxc\nPURTF8N7iSIFWA78wDm3DOjkta/qQ0n4bY52K9wPzAKmANlEuiUG89N+Hs6ltjGu255o4R7LZN0J\nycxSiQT7486556KrG8xscvT9yUBjdL0ffg83AveZ2UngSSJdM98B8qOTrMPrt+viNr/ZJOzjXA1Q\n45zbEl1+hkjY+3k/3waccM41Oef6gOeAG/D3fr7gcvdrXPd3ooV7LJN1JxwzMyLz0B50zn1rwFsD\nJx7/KJG++AvrPxI963490H7h61+icM59wTk3zTk3k8h+/J1z7iFgM5FJ1uGN25zQk7A75+qBajOb\nH111K3AAH+9nIt0x15tZVvS4/GByAAAA2ElEQVTf+YVt9u1+HuBy9+tG4A4zK4h+47kjum5kvD4J\nMYKTFncDR4DjwBe9ridO23QTka9fe4Dd0cfdRPoaNwFHo8+F0fZGZNTQcWAvkZEInm/HFWz/24Hn\no69nA1uBY8DPgPTo+ozo8rHo+7O9rnuE27oU2B7d1z8HCvy+n4G/BQ4B+4DHgHS/7WfgCSLnFPqI\nHIF/YiT7Ffh4dNuPAR+7kpp0haqIiA8lWreMiIjEQOEuIuJDCncRER9SuIuI+JDCXUTEhxTuIiI+\npHAXEfEhhbuIiA/9fyjHSrgmK/60AAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# control constants\n", "Kc = 2\n", "tauI = 60\n", "\n", "# control transfer function\n", "Gc = Kc*control.tf([tauI,1],[tauI,0])\n", "\n", "t = np.linspace(0,1000)\n", "\n", "H = Gp*Gc/(1+Gp*Gc)\n", "y,t = control.step(H,t)\n", "plt.plot(t,y)" ] }, { "cell_type": "markdown", "metadata": { "nbpages": { "level": 2, "link": "[4.3.2 PID Implementation](https://jckantor.github.io/CBE32338/04.03-PID_Control.html#4.3.2-PID-Implementation)", "section": "4.3.2 PID Implementation" } }, "source": [ "## 4.3.2 PID Implementation\n", "\n", "### Reference Conditions\n", "\n", "One of the implementation issues for PID control of the Temperature Control Lab is choice of reference conditions. One reason is that the linearization on which PID analysis is based is only valid in some 'neighborhood' of a nominal operating condition. But perhaps a more typical situation in most practical " ] }, { "cell_type": "code", "execution_count": 8, "metadata": { "nbpages": { "level": 2, "link": "[4.3.2 PID Implementation](https://jckantor.github.io/CBE32338/04.03-PID_Control.html#4.3.2-PID-Implementation)", "section": "4.3.2 PID Implementation" } }, "outputs": [ { "data": { "text/plain": [ "(35.294117647058826, 24.0, 0, 0)" ] }, "execution_count": 8, "metadata": {}, "output_type": "execute_result" } ], "source": [ "%matplotlib inline\n", "\n", "import sys\n", "sys.path.append('..')\n", "from TCLab import TCLab, clock, pid\n", "\n", "# ambient and reference values\n", "Tamb = 20\n", "Tref = 50\n", "uref = (Tref - Tamb)/0.85\n", "\n", "# control parameters\n", "b = 1 # setpoint weighting\n", "kp = 0.8 # proportional control gain\n", "ki = kp/60\n", "\n", "# sampling period\n", "tf = 1200 # experiment length (sec.)\n", "h = 1 # sample time (sec.)\n", "\n", "# setpoint function\n", "def Tset(t):\n", " if t <= 900:\n", " return 50\n", " else:\n", " return 35\n", "\n", "\n", "bi = ki*h\n", "\n", "r = Tset(0) - Tref\n", "y = Tamb - Tref\n", "\n", "P = kp*(b*r - y)\n", "I = 0\n", "\n", "uref,P,I,r" ] }, { "cell_type": "code", "execution_count": 9, "metadata": { "nbpages": { "level": 2, "link": "[4.3.2 PID Implementation](https://jckantor.github.io/CBE32338/04.03-PID_Control.html#4.3.2-PID-Implementation)", "section": "4.3.2 PID Implementation" }, "scrolled": false }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAA1UAAAGoCAYAAACubM8oAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4wLCBo\ndHRwOi8vbWF0cGxvdGxpYi5vcmcvpW3flQAAIABJREFUeJzs3Xd81eXd//HXdU7GyQ4JSQiEvRGR\nJSgggmLdVtE66vpVW6qtVettW+3dVupdW9u7t3a32mkdVdxWq9UKceBkyQ4zQEIge+9zrt8f35M9\nyD4Z7+fjkUfOd39OviGcz/e6rs9lrLWIiIiIiIhI17gCHYCIiIiIiMhApqRKRERERESkG5RUiYiI\niIiIdIOSKhERERERkW5QUiUiIiIiItINSqpERERERES6QUmViIjIAGCM8RhjrDEmpYfP+7Qx5ns9\neU4RkaFGSZWIyABmjClt9OUzxlQ0Wr420PF1hzHmmDFmSR9fc5gx5tfGmCP+n+E+Y8zPjTFx3Txv\nryREza7xkTHmut46v4iItE1JlYjIAGatjaz7Ag4DFzda92Sg42uLMSaov13DGBMGpAITgRVANLAE\nKAfm9nR8rVy/138mIiLSO5RUiYgMYsYYtzHm+8aYA8aYXGPMk8aYWP+2acaYWmPMzcaYTGNMnjHm\nJmPM6caY7caYQmPMQ43OdYsxZq0x5hFjTLExZqcxZmmj7XHGmL/7W5iOGGPuM8a4mh37W2NMAXCP\n//qpxph8Y0yOMeYxY0yUf/9ngUTgTX+L0e3GmPOMMfuavb/61ixjzIPGmKeMMc8YY0qAq9t7/624\nGYgHLrfWpllrfdbaY9baH1hr/+O/xsnGmPf8P5utxpjzG8XytDHmF8aYfxtjSowx640xY/2b3/V/\nT/O/n0vr3o8/vuPA7/3n+boxZr//frxgjEnqyr1vFFeQMeZ5Y8xxf9zrjDFTm+2W5F9fYox52xgz\nqjvXFBEZapRUiYgMbt8CPofT4pIC1AAPN9ruBmYBE4AvAb8G7gbO9K//kjFmYaP9lwKf4SQfDwIv\nGWOi/dueBIr851oAXApc3+zYLcBw4P/86+4HRgAnA1OB/waw1n4ByAY+5291+1UH3+/lwGNADPB8\nB95/YyuA16y1Fa1tNMZ4gFeBl4AE/7mfNcaMb7TbF4F7gTggC/hho/cOMNX/fl7yL48DgoHRwO3G\nmAuA7wOXAaOAXOCJDr739ryC0wI3AtiN8zNq7Hrgu/73tbeV7SIi0g4lVSIig9tXgXustUettZU4\nH/KvMsaYRvvcb62tsta+4l/+u7U2z1p7GPgAmNNo3yPW2t9Za2ustX8HMoBz/S0yS4G7rLXl1tos\n4FfA1Y2OPWCt/aO11mutrbDW7rbWrrXWVltrjwG/wEnmuuMda+2//K1MFR18/3XicRKhtpzh//6Q\n//3/G3gLuKrRPmustZustTXAU8DsE8RbBfyP/2dQAVwLPGqt3eqP99vA2caYESc4T5ustbXW2ses\ntaWNfgYL/ElinZestR/6t3/Xf82Erl5TRGSoUf9tEZFByp84jAb+ZYyxjTa5cBIIAK+1Nq/Rtgrg\neLPlyEbLGc0ucwgYCYwFPEBOo3zFBTTurnekWXwjgV8Ci4Ao//7tJTUdUX+NDrz/3GbH5gHJ7Zx7\nJHDYWtv4XIdwWpTqHGv0upymP7vWHPMnYI2vsbZuwVpbaIwp9l+j8ATnapV/rNaDwEqcVkIfYHB+\nBpn+3ep/btbafGNMqT+WnK5cU0RkqFFLlYjIIOX/8J8JnGWtjW305bHWNk8oOqp59boxwFGcD+Wl\nwLBG14m21jYu8GCbHfu/QBkw01obDXwZ58N+W/uXAeF1C8aYYJxudo3VH9OF9/8f4IJmLTiNHfW/\n38bG0JCYtKf5e2lr/VGcBBUAY0wMTsGMjlyjLV/C6QK5HKdb5LS60zfaZ3Sja8bhJIPdTXBFRIYM\nJVUiIoPbH4AHjTGjAYwxicaYi7txvtH+ohNB/vLdY4A3rbUHgY+AnxljoowxLmPMZNN+SfQonESs\n2BgzBrir2fbjOOOz6uwC4owxZ/sTqh9y4v/HOvP+/wzk44yTmmIcCf6CG2cD7wEuY8yd/vd/Dk6y\n8uwJYsBaW0XDeLP2/AP4ijFmpj+5+ymw1t89siOCjVO+ve4rCOfnXInTEhcB/KiV4z5vjFlojAn1\nb19nrc3u4DVFRIY8JVUiIoPbz3BaYNb6K+J9QPfKg7+LM8YqH6eoxGXW2iL/tmuAWJxCCPnAM0B7\nlet+gFNAogh4EaewRGMPAA/4K9bd5m9dugOnIEYGTle7E7W4dfj9+8c0LcPp0rcWKAE+xElENvnH\nG10EXIGToDwEXGWt3X+CGBq/32f97+eSNmJ4FfgJTmGJoziFJa5vbd82/AWny2bd1x9wksUcnJ/X\nNuD9Vo57AqeLYC4wHbixE9cUERnyTNOu4SIiIq0zxtwCXGGtXRHoWERERPoTtVSJiIiIiIh0Q0CS\nKmNMujFmmzFmizFmg3/dauNMPrnF/3VBIGITERERERHpjIB0/zPGpAPzG1dfMsasBkqttT/v84BE\nRERERES6SN3/REREREREuiFQk/9a4E3/ZIyPWGsf9a+/zRhzA7AB+C9rbUHzA40xq4BVAB6PZ96Y\nMc2nDJHBzufz4XLpecBQpHs/NOm+D02670OX7v3Q1F/v+549e3KttQkn2i9Q3f9GWmuPGmMSgbeA\nbwBpOKVcLfA/QLK19qb2zjN16lSblpbW6/FK/5KamsqyZcsCHYYEgO790KT7PjTpvg9duvdDU3+9\n78aYjdba+SfaLyDpoLX2qP97Ns7cJAustcettV5rrQ/4I7AgELGJiIiIiIh0Rp8nVcaYCGNMVN1r\nnNnotxtjkhvtdhmwva9jExERERER6axAjKlKAl40xtRd/ylr7RvGmMeNMbNxuv+lA18NQGwiIiIi\nIiKd0udJlbX2AHBKK+uv7+tYRERERESk+2pqasjIyKCysrJLx8fExLBr164ejqrjPB4PKSkpBAcH\nd+n4QFX/ExERERGRQSIjI4OoqCjGjRuHv0dap5SUlBAVFdULkZ2YtZa8vDwyMjIYP358l87R/+oW\nioiIiIjIgFJZWUl8fHyXEqpAM8YQHx/f5VY2UFIlIiIiIiI9YCAmVHW6G7uSKhERERERkW7QmCoR\nERERERnQ8vLyOPvsswE4duwYbrebhIQEAObOncurr75KYmIi27f3zqxNSqpERERERGRAi4+PZ8uW\nLQCsXr2ayMhI7r77bgDeffddbrvtNm644YZeu76SKhGRIcLrszz58SEmJUQyLCKEFzdnUl5diyfI\njfXvU1XrJSzYzednj2LmqJhWz5NVVMHjHx6iqtYHgMtAebUXT7Cbqlpvk/PVXTci1M0ZkxPYkJ7P\nFxeOJS4ipHffbAC8vi2LihovK+emBDoUERFpZOnSpaSnp/fqNZRUiYh0Q3Wtj/05pbiMYWJCBEHu\nvhmqWlheTY3Xkl1SibUn3h9gZ1Yxz27IAGB0XBhH8iva3PdwfjnfOGtyq9ue35TBhvSCTsf87p5c\nADILK7l24ZgOH3eo2Mv2zKJOX68rYsODqaj2UlXrIzzEjdtlMBiqvT4qa7xtHldcWcPvUvcDMCw8\nhISo0A5dz+0yDAsPIbe0qssxu4xhQkIExsD+7DKC3IYJwyM4XlJFkP/8Pms5kFPGmPhwIkP1X7+I\n9K6Lf/1+p4/xer243e529/nnN5Z0NaRep7+sIiLd8PvU/fxn13EAVs4dxZcWd21+i84oqazha09u\noqSytsvnaJ5Qfe/C6fz0jd3UeJ0M7WhhJfe+sK3N410uw3+dM4XQIBc/es2ZrHF4ZAi5pdUA3Hv+\nNNwuw+H8cv7+4aEmx27PLGr33M2VlFTz3KGO7x9oD721p8+vOX/cMEKD3Kzf5ySuxlCfbC+fmkBx\nZS0bDxUwNj6cX109B5dr4FboEhHpj5RUiYg08tmRQp7bmIG3g80/O48W179+bWsWe7NLeyu0esUV\nNU0SqhnJ0R0+tqSqhqjQYFwuWDxpODuPFnPhrGROGhnDg5fP4oVNmZySEsO7e3Pw+do+z6JJ8Syd\n4gwAXn3JDN5Jy+GmJeN56pPDTE6MYtGk4QCcOi6O4spaNh8uIMoTREFZDTFhnZutPsuUkTyi4++x\nqwrKq8kqcuYoGRUbRmZh08RzQkIEnqC2n6LWeH14rW13n+Z2Zjm/PyFBLiYlRHYhatibXcKG9AKM\naUimGv/6rkvLqX99KK+c7zy/leCg/lf8d3x8BDcvGa+ET2QQ6EqLUiAn/+0JSqpEZMA7nFdOTmnX\nJ+xr7I/vHmzxYfpElkweTlFFDdsyitiW0Tfd1OpcMS+FGxeN6/LxF80aWf96SlIU95w/DYDzT07u\n8DnmjY1j3tg4AL62bFKTbS6X4eYl44Gut+ClpuazbNmsLh/fUceLK/n6k5sYGRvGz66YxS1PbCTP\n3/KWHOPhoStn4+7hD/wPvZnGurQcvnfhdOaMGdalczzyzn5e3ZqFtbBiehIZBeXsPlbS5v7tbQuk\nbRlFjI4Lq+86ubfAS9ShfDzBbiJCghgR48ET3PGEtT3HiyuprvWRXeL83YiLCKWwvJq4CKe1NcTt\nYsbI6B6/3yIyeCmpEpEB7VBeGbf/YzO+Do4r6oj4yBDuXDGFjnyccrsMkxIj8VnLvuzSDo9v6q4o\nTxBBLhcpw8L65oJDQFK0h99fNw9PsAtPsJtfXjUHn7V4rSUkyNUrH7BvO2syVy0Yw6jYrt/Hm5aM\nZ7G/ZXBKUhQ1Xh8Hcsrqk5Pc0iqMgcmJUezPKaW6tp0myAD5ND2fl7cc5bfr9tevKymp5pWMnfXL\nS6cM51vnTuv2tXJLq/j6k5vqC6205frTxnLlqaO7fT0RCbxrrrmG1NRUcnNzSUlJ4Yc//CE333xz\nj15DSZWI9BuZhRX8bt0+5o4ZxoZDBcweHcP6fXlkFJSTHBtGdnYlz2ZuanJMcWUNPusUXkiM8nQ7\nBmPgolnJzB4d2+ljZ6V0/hjpXxoXmIgJ71w3xa4ICXJ1K6ECCHa7mlRqDAlycXJKw/KImIZ/F9M7\n0VW0L01IiODlLUcBCHIbTkmJJSOjjLjEKHZlOS1r7+7JJT13U3un6ZCy6tomCVV0WBDFFS3HJz6z\n4Qjv7Mlpsb4tI2I8fPu8qQS7XPz+nf2Eh7hPOMbyqY8Pc7y4ktvPnoy1lv97aw+H88qb7FNaVUuw\n20VVrZdLThnJF+Yr0RM5kdWrVzdZ/sc//tHr11RSJSIBsftYcYtiCS9syiCjoIKt/i50jSu+Hc4r\np6TcUpnf9AMHQLDb8L0LZzCymx9ORSQwojzBXD53FC9uzuS+i09i9uhYUlNzOGPpLG7/x2YO+//d\nH27l3393fPmM8UxPjubbz20lJiyYgvJqrj51DJsPF7D7WEmnrnc4v5y/rk8n2hPMG9uPARATFkyU\np/XkvLLGyz8+OQw4VScB3t+b2+41/v7hISJCgwhuo8poyrCwfps4iwx2AUmqjDHpQAngBWqttfON\nMXHAM8A4IB240lrb+Zq9ItLvZRVV8J3ntnaqy96wiBCunxjKaQvnttgWExbcJ60KItJ7bjh9HCvn\npRDdKAlxuwwPXXUKXp8lr7S6x7rXRnmCcBlDdFgQxhgev3kB4SFBlFbVEu0J4op5KRwr6vg4zS0Z\nhfzx3QO8tjWryfq/rk/v0PEvbMqsf339aWM5bUI8AMdLKrn/nzub7Pv71P20xWXgt9fOJWVYeAcj\nF5GeEsiWquXW2saPZO4B3rbWPmiMuce//J3AhCYi3bUuLZsnPzqMr5VPQRXVXnwWxg+PYEJCRJNt\nxj+SyfiHryyeFM/Oo8WcOSWRA9s+YUy8PiyIDEYul2mSUNUJ9VdTDI/rvY8sda1JdZUpQ4JMp/7W\npAwLo7Sytr7whcFgDK3+/WssyGWobfR0KS4ihMvmjqpviRoTH87tZ09mf04pPmvx+WyT/RtLzy1j\nf04Z33p2K2EhLQt6jB8ewb3nT+uzufREhpr+1P3v88Ay/+vHgFSUVIn0a2VVtWw5UsjpE+LZfrSI\nQ43GAjy7MYOCsuo2jw1yG755zhTGD49oc586dZXlDnQ/ZBGRHudyGb7YiQmtO+OcGUmcQ9IJ9zuc\nV86dz2ymtKqW0qqWY8RySqr48/sHT9hNOjEqlIX+ljIR6Thj+6pUVeOLGnMQKAAs8Ii19lFjTKG1\nNrbRPgXW2hb1ZY0xq4BVAAkJCfPWrFnTV2FLP1FaWkpkZNfmc5GetSatmp15XuYmudmc7W3RNSch\n3HDttJBWj/UEGTxBnaumpns/NOm+D026751XWWuprG35uW5Hnpe3DnV8svAbTwphfEzPlK/vCt37\ngSkmJoZJkyadeMc2eL1e3O7A/d4B7Nu3j6KiplOjLF++fKO1dv6Jjg1UUjXSWnvUGJMIvAV8A3il\nI0lVY1OnTrVpaWm9HK30N6mpqSxbtizQYQwp2SWV/OClHWQWVhDsbkiEarxN/35MSYpicpLzH6Hb\nGM6ansjELk5o2hrd+6FJ931o0n3vObVeH89vyqCgvKbd/TIKyvnsSBEuQ/0UAm7/XHPnzez43HXd\npXs/MO3atYvp06d3+fj+MPlva+/BGNOhpCog3f+stUf937ONMS8CC4Djxphka22WMSYZyA5EbCJD\nlbWW1LQcxsSHsz2ziOKKhv98P8soqp8Qt3kiVSci1M2dKyYzOk5jnkRE+pMgt4urTj1x98Si8hru\neGYzeaXV+Px/62u8lic+OkxOSRUAnmA3F80a2eq4LZFAysvL4+yzzwbg2LFjuN1uEhISKCkpYcyY\nMRw7dgyXy8WqVau44447evz6fZ5UGWMiAJe1tsT/+nPA/cArwI3Ag/7vL/d1bCJDSY3X12Sulo2H\nCnjorT0nPO7608Zy6ZxR9cvBbkON1xLkMrh6YXJUERHpGzHhwfzlxlPri2FYLN96disHc8tYsyGj\nfr+iihquXjAGtzFKrqTfiI+PZ8uWLYAzT1VkZCR33303WVlZZGVlMXfuXEpKSpg3bx7nnHMOM2bM\n6NHrB6KlKgl40TilvYKAp6y1bxhjPgXWGGNuBg4DXwhAbCJDQmZhBd98ZgsV1d429zn3pKQmE6FG\ne4IJcrs4a1pifbeQOiGdHBslIiL9k8tlCGn0N/47509j/d5cLJbyai8vbMrk5S1H6ydrrpMYFcov\nrp7d5rxcIoGSnJxMcrLTfTUqKorp06eTmZk58JMqa+0B4JRW1ucBZ/d1PCKDlbWWVz47SnZxVYtt\ne7NLqKj2Euw2TSaRLPevmzEymluXTWqRPImIyNAyKjaMK08dXb9cUF7DxwfysNDkwVx2SRU//3da\nq3NkjYkP59yTRvRFuNJfPHJmpw8J9/nAdYKS/199p4sBOdLT09m8eTMLFy7s1nla059KqotINxWW\nV9d329iWUcSf3jvY7v4PXTmbcR0oaS4iIgJw1zlTAOfB3a1PbKofbwuw6XAhmw4Xtnrc8MgQxsY7\n/9/EhYeou7j0udLSUi6//HJ+8YtfEB0d3ePnV1IlMkis+fQIj390qMX65VMTmJjYsgJfyrBwJVQi\nItIlxhjuv/Qk9h0vpdrrw1oormxZXXB7ZhEfHchn9Ss769edNDKan6w8GWOUWA1aXWhRKu/F6n81\nNTVcfvnlXHvttaxcubJXrqGkSmQA2pddwitbjta3SoFTaAIgLiKEuv+nkqI8fPXMiUSE6p+6iIj0\nrMQoD4lRnnb3WTo5gdzSnRSUO5PBF1XUsONoMQ+8touQoIauXnPGDOOcGSee5Fiks6y13HzzzUyf\nPp277rqr166jT1oi/VRuaRWF/v+Emvv12n0cyClrsX7aiCj+9wsthiyKiIgExLCIEB6+anb98lMf\nH+Yfnxzm44P5TfZbvy+XhKhQIkPdlFb3/RyqMnitX7+exx9/nJNPPpnZs53fxR//+MdccMEFPXod\nJVUi/dCR/HJuf3oztW3MCQUQExbMl88YX999wgCzUmL6KEIREZHO+8L8FCYkRDSZ0mPtruNsOlzI\n91/aDkB1eSWLF1czLCIkUGHKALd69er610uWLMHa3k/UlVSJ9BPPb8xgx9FiAI4XV1LrtcRHhhAb\n1rI8rcsYLpk9kmVTE/s6TBERkS4Ldrs4bUJ8k3WTEiP5zdq9VFR7ySur5kgJ/PCfO4iLaJjWIznG\nw01LxqsqrfRbSqpEAiC3tIpjRZX1y/ll1fztg/Qm+7gM3H/JTMbEtyxPKyIiMliMig3jJytnAbA1\no5Db//4h+3PK2N+sm/uwiBCmjXAKGYweFk5MuObEkv5DSZVIHyuprOEbT22mtKq2xbYzJg/nzCkJ\nACRGe5RQiYjIkDIrJZavnBzKpJOm16/bmVXMC5syeazRw8fhkSE8cv38JsUuJPCstQO2qmN3uwgq\nqRLpRT6f5Tfr9pGe2/C0ray6ltKqWuIiQhgZ21A1KTwkiJuXjCc+MrS1U4mIiAwJo6JcLGzURXD2\nmFgKyqrJKXUmsz+SX0FuaTV3PrMZT5C7fj+Xy3DVqaM5dVxcn8cs4PF4yMvLIz4+fsAlVtZa8vLy\n8Hjar2bZHiVVIr1g59Fi8sqqyCyo4K2dx1vd584Vk5kzZlgfRyYiIjKwhAa5uetzU+uXU9Oy+b83\n93Akv6LFvo+8c4DKGm+L9UEuF3ERIQyPDNHDy16SkpJCRkYGOTk5XTq+srKyW0lNd3k8HlJSUrp8\nvJIqkR62PbOIe1/Y1mTdFfNSWDSx4albRGgQI2PD+jo0ERGRAe/MKQlMTIhskjxZ4Kev7+Z4cSU/\neyOtzWOTYzz8/rp5KnjRC4KDgxk/fnyXj09NTWXOnDk9GFHfUlIl0kG/XbePDelN59XILa1meGTT\nkq9l1c4f+clJkSRGeYgND+aqU0fjCXYjIiIi3WOMYXRcyzHHd6yYzBvbj+FrNjTGYvlgXx4AWUWV\n3PzYp1w2ZxSfnz2qL8KVIUJJlcgJHC+uZP2+XN7YfqzV7bmlLSfojQwN4gcXzSA2XHNsiIiI9IVZ\nKbHMSoltddtDb+1h3e5sAPJKq3nio0N4gt24mo39SRkWxvTk6F6PVQYfJVUi7fD5LP/94naOFzvl\nz8+bOYIvzHf62/709TT2HC/h1HFx3LJsQpPjoj3BapkSERHpJ+48ezK3LZ9ERbWXB9/YxfbMYn6z\ndl+L/VwGfnvtXFKGqfqudE7AkipjjBvYAGRaay8yxvwNOBMo8u/y/6y1WwIVnwxd1lp+8vputmUU\n4bWWimovwyJCWDwxnmsWjiHa48yLce8F03htaxaXzR1Vv05ERET6H5fLEOIyhAS5+NqySbzy2VFq\nvU37CabnlbEvu5Q7n95CsNvFsIhg7v/8TIarsIV0QCBbqu4AdgGN21i/Za19LkDxyBC3fl8umQUV\nFFXU8OH+vPr1xsCNp4/l7OlJTfYfHhnKjYvG9XGUIiIi0h2j48L5+vJJLdZnFJRz1zOfUVHjparW\nR2lVLb9Zu48Z/u6ALpfhjMnDSYoOXIU66b8CklQZY1KAC4EHgLsCEYMION37fNaSdryEB1/f3WTb\n1QtGc8kpIwlyuQgLUVc+ERGRwSxlWDhPfHkhVbVe0nPL+e6L29h4qICNhwrq9/n4QB4/WXly/bIx\nRpUEBQDT3dmDu3RRY54DfgJEAXc36v53OlAFvA3cY62tauXYVcAqgISEhHlr1qzps7ilfygtLSUy\nMrLb58mt8PHnbVVU1DasmxjrIjnChScIFiYHEaw/lP1KT917GVh034cm3fehq7/c+605tWSXN3xO\n/vRYLVUtp8BiyaggVozVMIDu6i/3vbnly5dvtNbOP9F+fZ5UGWMuAi6w1n7NGLOMhqQqGTgGhACP\nAvuttfe3d66pU6fatLS25yKQwSk1NZVly5Z1+fi9x0v46EAeO44Ws+NoMS4DGENCZAg/vXyWJgXs\nx7p772Vg0n0fmnTfh67+eu9f2pzJ3z9Mp/FQLJ/PEuw2XDpnFM0fw05MiGTRpOF9GeKA1l/vuzGm\nQ0lVILr/LQYuMcZcAHiAaGPME9ba6/zbq4wxfwXuDkBsMkj5fJay6lp8Fn702i7yy5wy6EZVfkRE\nRKQDLp0zikvnNJ3bavUrO9h4qIBnN2S0eszDV80mKbrhYW14SJC6Cw5SfZ5UWWvvBe4FaNRSdZ0x\nJtlam2WMMcClwPa+jk0GJ5/Pcvezn7E3u7R+XVK0h8/NSGLc8AglVCIiItIlt501iXfScvA2m3F4\nS0Yh2zKK+OYzTQtZJ0aF8usvziE8RLMaDTb96Y4+aYxJAAywBbglwPHIAGat5aUtmRzOq6Csupa9\n2aW4XIawYBfBbhc3LRnHoolqkhcREZGuGx4ZyuXzUlqsXzJ5OPe9soOSypr6dVW1PrJLqvjxv3aR\nENmyguDUEZGcNzO5V+OV3hPQpMpamwqk+l+fFchYZPCoqPbySXo+f3k/vcn6mxaP4/OzR7V+kIiI\niEgPGRkbxh9vaDoM58P9efz4X7v47EgRDdOyNvjPruMkRHlIjGroLhgfGaJWrQFCd0kGlapaL19/\nahM5JU7hyGVTEzh5VAxhIW61TImIiEjAnDYhju9dOJ2iipoW2z4+mM8nB/NZ/cqOJutjwoJ55Pp5\nRITqI3t/pzskg0JmYQXPfHqE3NIqckqqCAtxMykxkq8snUC0R2VORUREJLCMMSycEN/qttljYimu\nqKG0qmGel4LyaooqavifV3cyLCIEgEUT4zljckKfxCudo6RKBqxar4/0vHKstTzx0SE2HS6s33br\nmRNZPi0xgNGJiIiIdExilIf//cIpTdbVdRfccbS4ft2nB/OJiwghxO1qcY7RceF4gt29Hqu0TkmV\nDFi/Xbef/+w6Xr8c5DbcftZkosOCmTsmNoCRiYiIiHTPaRPi+J9LZ1Ls7y74z8+OsvtYCfc8v63V\n/WeOiuEnK0/uyxClESVVMuAcKfHx/Ze2szWzCJdxJtcDWDYtUa1TIiIiMigYY5g9uuEh8chYD396\n7yDVtb4W+x7KL2d7ZhHffu6zVufBCgsO4pZlE0iMall1UHqGkioZMHw+y7bMIv51sJpS43T1Wzpl\nON86d1qAIxMRERHpXZMSo3hax5TgAAAgAElEQVTw8lmtbnvio0M88+kRdmWVtHl88HuGs6cnAWAM\nTEmMIiZc4857ipIqGTBe2pLJX9enU1JqGTHczfcunMHkpMhAhyUiIiISUNcsGMO8scOo9doW24or\na3jw9d18sD+PD/bn1a8fNzyCX109G2NatmxJ5ympkn7rqY8P8+6enPrl3FKnTPqEGBdfPmsyM0fF\nBCo0ERERkX7D7TJMT45uc/vNS8bzWUZDQa/dWSWk55ax6vGNuFtJqsJC3MwdO4yP9ucxe3QsX1k6\noVfiHkyUVEm/4vNZNh0uoKC8hmc+PYyv2QOXiQkRfD6phkWTNOeUiIiISEdcOmcUl84ZVb/88pZM\n/vTeQY4VVbZ5zL7sUgAO55eTHOtpdYqahKhQQoNcTEhQzyElVdKv/HvHMX6Xur9+edHEeK47bWz9\nclK0hw/efzcQoYmIiIgMCpecMpLTJsS3WvSitKqWbz+3tcm6R9450O75Hr5qNpMSh3ZipaRK+oU/\nvLOfjw7k1ZcNnTMmlviIUK5ZMJrEaFWqEREREekpxhiS2vl8dcfZk9lypJBFk+LZmF5ARY23xT7v\n7c2tf33vC1uJCG2aVsweHUt2SRULx8fx+dmjmh8+6CipkoDKLKxg/d5cXtuaVb9uTFw49118Uqsl\nQUVERESkd62YkcSKGU6lwEUTWx9yccZkZ3JigMoaH5U11U22v70rG4BtGUUtjg0NcrNsasKgmqxY\nSZUEjM9nuf+fOzha6PTnPfekJK5eMIbYsGAlVCIiIiL92OkT43n+1kX4rKW0qrbJtsc/PMTa3dn1\ny39672CL47OKKvjS4vG9HmdfCVhSZYxxAxuATGvtRcaY8cDTQBywCbjeWlvd3jlkYDqSX853X9xG\ncUUNPgvDIkJYMT2RS+eManUQpIiIiIj0PyFBLoAWLU43LR5PXEQIXn/FsVpfw9itqhofb+48zgub\nMnlpy9H69baykmXLej/m3hLIlqo7gF1AXf3HnwIPW2ufNsb8AbgZ+H2ggpOeZ61lXVo2b+3MprDc\nGTvlchluOG1sfROziIiIiAxsMeHB3LhoXJvbiytr+OhAPr5GZZ5tyym2BpROJ1XGmElAkrV2fbP1\nZwBHrbX7Wz+yyb4pwIXAA8Bdxpl17Czgi/5dHgNWo6Rq0PD5LJ+k5/PwW3vr1/3u2rmMig3Dpa5+\nIiIiIkPGf184g1pv08qD77zzToCi6Rldaan6BfDdVtZX+Ldd3MFzfBuI8i/HA4XW2roOmRnA4C8T\nMkRUVHu585nN9WOnTpsQx/KpiYyOCw9wZCIiIiISCEFuV5PlgT6evitJ1Thr7dbmK621G4wx4050\nsDHmIiDbWrvRGLOsbnUru7baCGiMWQWsAkhISCA1NbVjUUufy63wseGYl7xKH3sLfBgDw8MMc0Nr\nqM7IJjWja+ctLS3VfR+idO+HJt33oUn3fejSvR+aBvp970pS1d6kQWEdOH4xcIkx5gL/uaJxWq5i\njTFB/taqFOBoawdbax8FHgWYOnWqXTaQR7QNYl6f5fZ/bOZwaTkAUVFw7/nTWDSp9bKcnZGamoru\n+9Ckez806b4PTbrvQ5fu/dA00O+768S7tPCpMeYrzVcaY24GNp7oYGvtvdbaFGvtOOBqYK219lpg\nHXCFf7cbgZe7EJv0A9ZavvXcZxzOLycs2M1Xlk7gO+dN4/SJ8YEOTURERESkx3WlpepO4EVjzLU0\nJFHzgRDgsm7E8h3gaWPMj4DNwJ+7cS4JgMLyah774BBFFTXsPV4KwPWnj+XiU0YGODIRERERkd7T\n6aTKWnscWGSMWQ7M9K9+zVq7tgvnSgVS/a8PAAs6ew7pH/LLqnnio0P8Z9fx+nVfXDhGCZWIiIiI\nDHpdnqfKWrsOp8ueDHG7soq55/mt1E01cMuZE0mICmXumNjABiYiIiIi0gcCOfmvDHA1Xh9/fO8A\nmw4V4LMQHxnCkknDuXBWcqBDExERERHpM0qqpEuOF1fy9q5sXt92DIDQIBcPXzmbYREhAY5MRERE\nRKRvKamSTssrreLrT26iqtaZCfvqBaNZPjVRCZWIiIiIDEldTqqMMSW0nKC3CNgA/Je/8IQMMn95\n/yAf7M+lqtZHYlQos0fHcuX80QS7u1KdX0RERERk4OtOS9VDOBP0PgUYnDmnRgBpwF+AZd0NTvqP\n/LJqPk3P58XNmQC4DNx97lSmJ0cHODIRERERkcDqTlJ1nrV2YaPlR40xH1lr7zfGfLe7gUn/4fNZ\n7nl+K1lFlQCsmJ7ENQtGkxjtCXBkIiIiIiKB152kymeMuRJ4zr98RaNtzbsFygCUUVDOQ2/uobiy\nhuPFVUSHBTFn9DBuXDSW2HCNnxIRERERge4lVdcCvwR+h5NEfQRcZ4wJA27rgdgkQKy1fHwwn1e3\nHmVvdmn9+utPG8t5M1UuXURERESkse5M/nsAuLiNze939bwSeB/sz+PB13cDztipH688mYSoUBKj\n1N1PRERERKS57lT/mwL8Hkiy1s40xswCLrHW/qjHopM+VVJZw3+/uJ3D+eUAzB0Ty7JpiZw0MibA\nkYmIiIiI9F/dqYP9R+BeoAbAWrsVpwKgDEBbjhTyyDsHOJhbhtdnSYr2cO8F01k+NTHQoYmIiIiI\n9GvdGVMVbq39xBjTeF1tN+ORAMgsrOC+l7fj85cXuef8aZw+IR6Xy7R/oIiIiIiIdCupyjXGTMRf\n6c8YcwWQ1SNRSZ85kl/O157cBMDUEVGcMXk4iybG0yxZFhERERGRNnQnqfo68CgwzRiTCRwErjvR\nQcYYD/AuEOq//nPW2vuMMX8DzgSK/Lv+P2vtlm7EJ+2w1vLatizeScupX3fH2ZMZHRcewKhERERE\nRAae7lb/W2GMiQBc1tqSDh5aBZxlrS01xgQD7xtjXvdv+5a19rl2jpUeUOv18d7eXB555wAAxsAj\n188jOSYswJGJiIiIiAw8nU6qjDF3tbEeAGvtQ+0db621QN3kR8H+L00W3EdKq2q57alN5JVWA7Bs\nagLLpyUqoRIRERER6SLj5DidOMCY+/wvpwKnAq/4ly8G3rXWfrkD53ADG4FJwG+ttd/xd/87Hacl\n623gHmttVSvHrgJWASQkJMxbs2ZNp+IfyjYdr2Vrrpf0Ih9BLkiOcHHdjBBC3QNr/FRpaSmRkZGB\nDkMCQPd+aNJ9H5p034cu3fuhqb/e9+XLl2+01s4/0X6dTqrqDzTmTeDyum5/xpgo4Flr7XmdOEcs\n8CLwDSAPOAaE4IzV2m+tvb+946dOnWrT0tK6FP9Qs/d4CXet+ax++b6LZzB/XFwAI+q61NRUli1b\nFugwJAB074cm3fehSfd96NK9H5r66303xnQoqerOPFVjgOpGy9XAuM6cwFpbCKQC51lrs6yjCvgr\nsKAbsUkjh/LK6hOqeWOHsfqSgZtQiYiIiIj0N92p/vc48Ikx5kWcMVGXAY+d6CBjTAJQY60tNMaE\nASuAnxpjkq21WcYZnHUpsL0bsQng81me+PgQHx/Ir1/3laUTGBWr8VMiIp3i88Kmx6DwCCRMhVlX\nOVV+Gqsuhw1/gUkrIHGasy7tDTjyMVgvjJoP5bnOOQB8NeAKhonLIWMDVJeCt6bhvO5QmHQ2HP4Q\nxixyvs+9AcJ74aFYUSZsW+NcJ+MTmPclqKmAzY/DKddAdHLPX1NEZBDpTvW/B/xV+87wr/qStXZz\nBw5NBh7zj6tyAWusta8aY9b6Ey4DbAFu6WpsAkXlNby7N4dnN2QA4HYZfnftXEYqoRKR3lZdDt5q\nCIuF2mrIPwCuIGc5NAqMG/L3O4lKB0SUHoLjO3s56BPI2gIb/c8N96+FsGEQM7rpPjtfhj1vwM6X\n4OJfOUlJ6k8ath94p/Vz71/b9nX3vul83/GS870oA+bf1LX30J7UHzvJXt11KgqgOAuyd8KxbbD0\nWz1/zTrGwLBxUFUKnhgoyQJ3MGHlWb1z311uiEhwXpdmQ8RwKMuF8HjnfcdPdNZHjWiZOJfnQ3AY\nuEOcY6KSWp6/NMf5Xbc+53ffWucaFYXgq23YLzy+6fElx524XN3pRCQigdKdliqstZuATZ08Zisw\np5X1Z3UnFmlQVlXL157aSHGF88f7wlnJfG5GkhIqEel9Ph+88g0oPQ5feAw++BUcSG3YnnIqRCbC\n7tc6fMppJSWQGdXzsXbHuh+3vc1bAy/d2vVzB4c5CVlrMj51vnrbvrcbXucf6N776QjjTyRcQU5C\nDswoKYEjAbjv7hAnhsV3wMyVDesLj8DzX3ZaKodPhm3PwXk/gbGLGvY5ugVe/SZMOddJwI583PZ1\nXEFwxZ+dhDJ9Pfz7u3DyF2DRbb321kSk93QrqZL+p6Syhq89uYniilpiw4OZkRzNjaePIyzEHejQ\nRKS3leXCuz+Hshz/B/NyCO6BCb0rCpyWmY7wVkPePuf1P2+H4kznaX9dUaSMT50P0MblfDjtgDLf\nMaISR3Qh8B5WVQIhEVBb2fbPtSwXIpOc7n4AQR4YfwbsedP52USPhMTpsO8/zn6jF8C+tVCeByHh\n4ImFJXfCnn/7W8TinPsZMdzZp6P3odOMc53IJKelKDLRWV1R6LS69Kb8g87PFOoTqnrxk8Ad3HPX\nshZydre/T10M63/ZtBWxLNeJM+sz5wvgjXthxMkN++Ttd1qo0vzTb7pDICjU+d0BSJjm/HuoKICS\nY/D6Pc69PbbN2b7t2RPH11/FjoUz7nK6su57CxbdDp7oQEcl0meUVA0iR/LLeeWzoxSW1wBw54rJ\nzBurghQi/Vp5PhQecj5sWZ/zgaquKGtIhPMBrKq03VPU2/O6M+6mNxRldP2Y6Rc7H5yP+4fKWh9M\nPAtW3Nf2sY3sTk1lRD+sCNUpMy9vujznura3ASz4ivM1FGz8G2z4a9N1xkV5+EiiVv6x57vDvXVf\nQ+upy+10QzUu5/ey8QMAaEh22tPePidd5oyZe/nrMPo0OM/fwllwCJ67yUlgS7I6fr7+7Ng256HB\nJ486y9YHUy/s0qmiivdARidaKY3xJ6wuZ9xiTEqXrivSHUqqBon03DLueHozPv//BV89c4ISKpH+\nzlsLL33N+VA1fqlTqCCzUz2qT2zhLZA0o+vHv32/84R++JSOd0sK8jhflYVOF6fhU5wucVUlUJHv\nfIgdPqXrMcngMucGJ+HwRDstOFHJYC17PvyEpN4YX7T8v+G0W53fyehRTmtq1AinpS5smNPSGpHg\ndGFtLiwOqsucIiPhw6Esu+U+EQnOwxLjcn7P3UFw9T+atvgNGwtXPe5cs/64RCchsL6ef8+97cgn\nsPmJhoQKnC6kjbuRdsKUkhLI6mTXz9ELnVbNQx/ARQ/DyNldurZIVympGuB8Pssv3t7LZ0cK8VkY\nGethwfh4LjxZlZpEet2nf3I+TJTlNAxG78zcf96qhqfUB991vgd5nCSougxy/PPwDZ/sFHjoiBEn\nOx8Ws7Y655l1pfM0vqsu+LlTAe7Ur3ShAtzYhpfuYKd7W2sD+2Voc7kaqiVGj6xf7Q3qpXHAQSFO\nElUndnTTa9d154vqQJfTtv5NNHofzrla+b2PHtlyv4FaZTFhmvPwpSzbSSYriyG065O4Fh89StTI\nkSfesc7xHU3Hr71xj5Mw95TKQgiJhNoqpwdBnWAPnPFfEDeh564lA5aSqgHsSH457+zJYd1u50lZ\nsNuw+pKTSI5RQQoZxI5uaXi664lx/jPL2uJ/apzT/rE9qboMNj3edF1FQdfOFRzujH8CpyvYwlVO\n1bxnb3SqhV3ya2eMVCDEjYezfxCYa4vIwBAUCsvv7bHT7U1NZVRnuvx+9Hv47OmG5ZqKhrGdPaUs\nt/X17//C6eLcVUkzofRY0/MnzoCYHkwKpU8oqRqgyqpqufvZzyivdgZDX3nqaC46OZlhESEBjkyk\n90SUHoR//i3QYbQtagR87kedOybI4zytLkh3nvDG+lt3gkLg8j8DNnAJlYjIQLDgq07FRWudLpzl\neT137oxP4eNHGpZXrHbGbNVUwGt3NS1c0lPC4+Gap53/B2TAUFI1ABWV1/Clv31CjdcyIsbDwvFx\nXDk/hdAgVfiTfmTvf5xKVitWO11aPn7EKc3sDnG6anTB5NwsCA+BpJOc8xxtZWq8SSu6FXanBIfB\nuDMg/T1nedaVEDuma+eKn9hyXUgPVO4TERnsXK6mXfB6coLsuAlOl+qSY05X7InLG7YtuwcOdaM4\nUNaWhhaqpJnOg7ljW5150p68whmT2lhEAlzwM6eXhvQ7SqoGmIpqL99/eTs1XmfcxteXT2L26F4u\ndytSx+eDA+ucgdQn6ur24W+d7299H8afCVue6vbl3d4qCI13WoOCw2DNjc5g8rjxTkvP5x6AcYu7\nfZ1OG7Ow768pIiK9z+WGeTe2vm3Siu49yNvzJqx7wClicu6PnFa2A+/AWz+AyqKW+5fnwXsPOQ8W\nT8QT65zDHQyTP6eHdH1ASdUA8/hH6RzMLQPgexdOV0IlfWvrM/DxHzp3TO5e56uxU66Bk6/o/OXX\nr2fxWec5/fcBrn7SGdsUFus8SezJ+WxERER605TPwZjT/BVT/V39JpwJN/6z5Zxt2bvgze850wE0\nnlC9IwoPOZNZS69SUjVA7Msu5Uev7SS/zPlHtnLuKBaMV8l06ab8A/D6dzo+D1JdMQVwqj01nvSy\nOZfb6bpQW+Usu0OcFq6gUDjl6i6NE6oJiWlIqMBJourKFCuhEhGRgaa1CZJbWzduiVNpsPDwic+Z\n9i/ngWOdHS9C2hsdj8kY5/jmk5xPOdeZnFxapaRqADheXMl3nt9Kda0zd8WiifF8afH4AEcl/UZ1\nuTN7vbfGmeeoM4586vTd7ixPtDNxa/NywCIiItLzjIEZl3Rs3/FL4V/fgtNugYwNkP5+04eiHdX8\nmB0vOiXlO/UQ0zhTglSXOV0lB+q0AR2gpKqfq6718Yv/7KG61ofbZXjk+nkkRoWe+EAZmLy1Tgnt\nzlj/S9jTiSdQzRkXfOFvzhxLHREU5gwKFhERkf4neRbc/G/n9YxLO59QvXWfU/Vw5ko49cvOug9+\n47SAbX6i63Ed+diZ+7ANxlfT9XP3A0qq+rG3dh7nN2v34vPPJfrt86aSFO0JbFDSe7J3wT/vaOgu\n1xUTzmwoyd1RSSfBsE4eIyIiIv2fMU0nLO6Is77nFKWaemHDWK/TbnXmzursZ5RNf294fWwb/OXc\nNnc9udLAWed07vz9iJKqfqi4sobnN2awLi0HnwVPsIurTx3DookdbEmQnuXzwdannSpzvSlrq/PH\nyh3stB51Rm2VU471zHtU4UdERES6LiwWTrqs6TpPNMy5rvPnipvgTKkyYiakrwfrbXNXn8vX+fP3\nI32eVBljPMC7QKj/+s9Za+8zxowHngbigE3A9dba6rbPNHj96d0DrEvLASA5xsMfrpuHy2UCHNUQ\nU1XS0A3vyCdNJ/7rTe4Q+OIzPTvHhoiIiEggTFzedG6vdmxPTWVZ70bTqwLRUlUFnGWtLTXGBAPv\nG2NeB+4CHrbWPm2M+QNwM/D7AMQXUH96ryGhuuH0sSyeNFwJVV/b/CR88mjL9VMvcCb+600J05RQ\niYiIiAwwfZ5UWWstUFeiLNj/ZYGzgC/61z8GrGYIJVWH88p5aUsm69KcSmxnTkngC/NHBziqIWLH\ni854pjrp653vnmjAn9DGpMCib6hrnYiIiIi0YJwcp48vaowb2AhMAn4L/C/wkbV2kn/7aOB1a+3M\nVo5dBawCSEhImLdmzZo+i7u31Posv9lcRWGVcy+mxbm4epoq/LWltLSUyMjILh8fUpWH2z+pXkh1\nPpP2/bnlNSLHkTbt9i5fQ3pHd++9DEy670OT7vvQpXs/NPXX+758+fKN1tr5J9ovIIUqrLVeYLYx\nJhZ4EZje2m5tHPso8CjA1KlT7bJly3orzD5hreWbz2zBG1JGYqSbW8+cyPxxw4jyaCLTtqSmptLl\n+773LVj7q6broqJg/BkwdomzbFxEpcwnWd3w+p1u3XsZsHTfhybd96FL935oGuj3PaDV/6y1hcaY\nVOA0INYYE2StrQVSgKOBjK0vfLg/j39uPcr+HGfW65sWj2f5tMQARzXI7P4X7F/bsJy3z/kePaqh\nTGhoNCy6AyIT+j4+ERERERnwAlH9LwGo8SdUYcAK4KfAOuAKnAqANwIv93VsfeloYQU//lfDOJ7r\nTxvLeTNHBDCiQaQ0BwoPg68G3n8YvM2KSEYmwpWPdXJGcBERERGR1gWipSoZeMw/rsoFrLHWvmqM\n2Qk8bYz5EbAZaDnQZZAor67l609tAmBsfDi3nDmR6cnRAY5qkKgug+dvgsrihnUjZ8PsaxuW4ycp\noRIRERGRHhOI6n9bgTmtrD8ALOjrePra8xszWLs7m1qvM2Ts1mUTOWlkTICjGuDS3oDdrwIWqkqd\nhCpiOMSOdeZ9WrjKmXxORERERKQXBHRM1VDzwb5c/vZBev3y9y6croSqC4yvBg5/BN4asD5Y/0uo\nKW+605K7YNziwAQoIiIiIkOKkqo+cjC3jJ+8vhuA0yfGc93CsYyJ15xHXZGS8U/Y+1nTlUkzYeFX\nndehkWqZEhEREZE+o6Sql/l8lp/9O42tGYUABLsNX14ynsRoT4AjGyCsdYpN7HrFqdhnLQk5eyEy\nHMYuBmOc8VGzr4PhkwIdrYiIiIgMQUqqepG1ll+v3cf6fbkAhAa5+M0X5yqhao/PC/vXQY1TZp7y\nfNjpLwRZlAGAsV4YtwTOfSBAQYqIiIiINFBS1YtS03L4z67jAFy7cAwXnTKSyFD9yNu140X44Net\nb4tMhIseZvuHH3H6OSv7Ni4RERERkTboE34vKK6sYfXLO0jPc1pbpiRFsXJuCiFBrgBH1g9t/JvT\nEmWdaohUlzrfxy8Fj7+IhzsY3KEw7UKISaHKMxxc+lmKiIiISP+gpKqHVdZ4ue/lHezLdpKDpOhQ\nfrxyphKqOhWFsP9t8NY6lfs2PQ6+2qb7xE+EFT9U4iQiIiIiA4KSqh721MeH6xOqe8+fxoLxcQS5\nlRzU++BXsO/tpuvGLoaldzcse2KUUImIiIjIgKGkqgc9+fEhXtycCcDFpyRz+sR4jDEBjiqAdrwI\nH/ymZUuUccHMy/2V+0JgxqUQHheYGEVEREREuklJVQ84VlTJmzuP8ZI/oTp1XByrlk4McFQBsO9t\nyN3TsJz2esuECmDmSlh0W9/FJSIiIiLSi5RUdVN1rY/frNvLZ0eKAJg3dhg/uHhGgKPqQzWVgIWi\nTHj7/pbb4yfByj86rVJ1hnLrnYiIiIgMOkqquiE1LZuH39qDz0JIkIsbTh/LGZMTAh1W33n7/pbj\no0bNg5RTndfG5cwnpfFRIiIiIjKIKanqou2ZRfzfm05XN0+wiyvmpfD52aMCHFUv8nlh8xNQntuw\nvO9t/7ioUGedJxqWfBNiRwcuThERERGRPqakqpOstWSXVPHzN9MAiIsI4c83zh98Ff4qixvmjAI4\n9AFs+EvL/WZ83kmkRERERESGqD5Pqowxo4G/AyMAH/CotfaXxpjVwFeAHP+u37XW/quv4zuRP79/\nkJe3HK1ffuCymYMvocpJg5dudVqjmjvpMhg2znntDoEJy/owMBERERGR/icQLVW1wH9ZazcZY6KA\njcaYt/zbHrbW/jwAMXXIy1sy6xOqUbFhfGnxOFKGhQc4qh5wdLNT/tz6nOWCdCehCouFoLCG/YaN\ng9O+BkEhgYhSRERERKRf6vOkylqbBWT5X5cYY3YB/X4w0vbMIv703kEAFoyP4/sXDdAKf/kHobqs\n6bp3/heKM5uuc4fAyj9B5BAqvCEiIiIi0gXGWhu4ixszDngXmAncBfw/oBjYgNOaVdDKMauAVQAJ\nCQnz1qxZ0+tx5pT7+O2WKgDiPIZVs0LxBA28suCxBduYuP+vrW6rCo0jI+WS+uVKTyKVYSP6KrRO\nKS0tJTIyMtBhSADo3g9Nuu9Dk+770KV7PzT11/u+fPnyjdba+SfaL2BJlTEmEngHeMBa+4IxJgnI\nBSzwP0Cytfam9s4xdepUm5aW1msxWmv52wfpfHQgj6OFlQD84fp5jIoNO8GR/cT+dbD7NZwfKU63\nvrJcpxtfSKNfWpcLTrkGxi4KQJCdl5qayrJlywIdhgSA7v3QpPs+NOm+D12690NTf73vxpgOJVUB\nqf5njAkGngeetNa+AGCtPd5o+x+BVwMRW53c0io+3J/HC5ucbnEul+HXV8/p3wlV/kEnaQJnfNR7\n/wdVJU338UTDpb+HkEEwFkxEREREpB8IRPU/A/wZ2GWtfajR+mT/eCuAy4DtfR1bncoaL998ZguF\n5TUAXDgrmQtPTmZ0XD9ORPL2w/Nfbig2USdhKpz6lYblYWOVUImIiIiI9KBAtFQtBq4HthljtvjX\nfRe4xhgzG6evWjrw1QDERllVLbc+uYnC8hqGRYRwSkoMN5w+lvCQfjSl16bH4cC6pusqCp2EKnYM\nRCY661zBMO9GSJze9zGKiIiIiAwRgaj+9z7QWpWHgM9JtftYMW/uOE5BWTUAt545kdMnxgcmmNIc\np9R5c94qZxLe5i1SAEGhcP5PIXpk78cnIiIiIiJAgMZU9UdH8sv5znNb8flrOtwSyITK54N/3e0U\nlmjL+KUw98am68LjnC8REREREekzSqqAzMIKvvbkJgDGD49g8aR4zp/ZB+XEa6vhrR9A4aGm631e\nKD0OnhgYvbDlcUEhMOcGiErq/RhFRERERKRdQzqpOpRXxo6jxXywP7d+3d2fm8qY+B4u5JC92+ma\nV5YDxUcb1hcchMMftn3c3Bvg5Ct6NhYREREREelRQzapqqzx8t0Xt1FcUVu/7uGrZvd8QlWcBS/d\nAu3NB3b6bS3niHKHQGRCz8YiIiIiIiI9bkgmVX9df5BXt2ZRXesjKdrDnDGxnDQymkmJXZzFec+b\nsOkxp3z5oQ+bFpGw3jeutuIAACAASURBVKYJVfIpTlnzOuHDYeZKcLm7dm0REREREQmoIZNUeX2W\n17ZlUVRezUubM/FZZ0LfVUsnsGB8o+IO1sK+/0BNhbPsq4HaKqitbPvkGx9zvhdltB9E2DA4+wcQ\nMbx7b0ZERERERPqNAZ1UGSxUl59wP6+1/Hv7Mf6+/iAAwcBpE4Zzx4pJBLtcTc+R/j6se6DrQY2a\nB+c2O94VDFgwbnC5un5uERERERHpdwZ0UhVengF/Pb/dfWp8liP55Uz2WX4CRIcFE+w2RB8JJuix\n1qbLasOsqyA47P+3d+dxdtVlgv8/TypLJamQnRCSQBZqIgQ0GFQE1CC2oqJBR1sYx0G0B+2Wn9rN\nNGrjb5huh1fjtCO92E037QbdjkqzjHEXwYBAIxL2EDCBkFBkJXtSWaue+eOeSipJZb1169at+3m/\nXnnVPd/zPec8ladO1X3u+Z7v6Xpd9IOBQ0tXt5rffvB+kiRJkvqcmi6qkn4HLWC27Wpj47ZdbN/V\nzrYcSAQMHtjA8SOHEHGYYqpj6F+H6e+EN/5RN0UtSZIkqS+p6aJqY+ME2j/6UwB2tydbd5Rm8kvg\nz+54ipd3boMGoBH+1wdezZTxx1UvWEmSJEl9Uk0XVWta2/nAPz7IiCED2bazjS07du+zftTQgcyZ\neSLjhw/mVAsqSZIkSRVQ00UVwK62ZM3mHQA0DuhH44DS1OQN/YLLzpnM+dOPr2Z4kiRJkvq4mi+q\nOrv2PTM4fcLwaochSZIkqY70eFEVEZOAW4ATgHbgpsz8m4gYBXwfmAy8CPx+Zq4/1L5GD+7H/3/R\nabRnMnRgfwsqSZIkST2uGleqdgNXZeajETEMmB8RdwEfBe7OzOsj4vPA54HPHWpHgxrY98G9kiRJ\nktTDevxJtJm5IjMfLV5vBhYCE4A5wM1Ft5uBi3s6NkmSJEk6WlW9pyoiJgNnAr8BxmXmCigVXhHR\n5QwTEXEFcEWxuCMinu6BUNW7jAFeqXYQqgpzX5/Me30y7/XL3Nen3pr3k4+kU2RmpQPp+sARTcC9\nwHWZeUdEbMjMEZ3Wr8/MkYfZxyOZeValY1XvYt7rl7mvT+a9Ppn3+mXu61Ot573Hh/8BRMQA4Hbg\nO5l5R9G8KiLGF+vHA6urEZskSZIkHY0eL6oiIoBvAAsz86udVs0FLiteXwb8oKdjkyRJkqSjVY17\nqs4FPgI8FRGPF21/BlwP3BoRHweWAR88gn3dVJkQ1cuZ9/pl7uuTea9P5r1+mfv6VNN5r9o9VZIk\nSZLUF1TlnipJkiRJ6issqiRJkiSpDDVbVEXEhRHxXEQsjojPVzsedZ+ImBQRv4qIhRGxICI+U7SP\nioi7ImJR8XVk0R4R8bfFz8KTEfHa6n4HKkdENETEYxHxo2J5SkT8psj79yNiYNE+qFheXKyfXM24\ndewiYkRE3BYRzxbn/Rs93+tDRPxx8Xv+6Yj4bkQ0es73PRHxzYhY3fnZosdyjkfEZUX/RRFxWVfH\nUu9xkLz/VfG7/smIuDMiOj9O6QtF3p+LiHd0aq+J9/w1WVRFRAPw98A7gdOASyPitOpGpW60G7gq\nM08FzgY+VeT388DdmdkM3F0sQ+nnoLn4dwVwY8+HrG70GWBhp+UvAzcUeV8PfLxo/ziwPjNPAW4o\n+qk2/Q3ws8x8FfAaSvn3fO/jImIC8GngrMw8HWgALsFzvi/6NnDhfm1HdY5HxCjgWuANwOuBazsK\nMfVa3+bAvN8FnJ6ZrwZ+B3wBoHifdwkwo9jmH4oPWWvmPX9NFlWUTqbFmflCZu4EvgfMqXJM6iaZ\nuSIzHy1eb6b0BmsCpRzfXHS7Gbi4eD0HuCVLHgJGdDzzTLUlIiYC7wa+XiwH8FbgtqLL/nnv+Hm4\nDbig6K8aEhHHAW+m9KgNMnNnZm7A871e9AcGR0R/YAiwAs/5Picz7wPW7dd8tOf4O4C7MnNdZq6n\n9OZ8/zfs6kW6yntm/iIzdxeLDwETi9dzgO9l5o7MXAIspvR+v2be89dqUTUBeKnTckvRpj6mGN5x\nJvAbYFxmroBS4QUcX3Tz56Hv+GvgaqC9WB4NbOj0C7hzbvfkvVi/seiv2jIVWAN8qxj2+fWIGIrn\ne5+XmS8DX6H0GJUVlM7h+XjO14ujPcc99/uejwE/LV7XfN5rtajq6pMp54bvYyKiCbgd+GxmbjpU\n1y7a/HmoMRFxEbA6M+d3bu6iax7BOtWO/sBrgRsz80xgK3uHAXXFvPcRxdCtOcAU4ERgKKUhPvvz\nnK8vB8uz+e9DIuIaSrd7fKejqYtuNZX3Wi2qWoBJnZYnAsurFIsqICIGUCqovpOZdxTNqzqG+RRf\nVxft/jz0DecC742IFyld3n8rpStXI4qhQbBvbvfkvVg/nAOHl6j3awFaMvM3xfJtlIosz/e+723A\nksxck5m7gDuAc/CcrxdHe4577vcRxSQjFwEfzr0PzK35vNdqUfVboLmYIWggpRvb5lY5JnWTYoz8\nN4CFmfnVTqvmAh2z/VwG/KBT+38pZgw6G9jYMaRAtSMzv5CZEzNzMqVz+p7M/DDwK+ADRbf9897x\n8/CBon+v/PRKB5eZK4GXImJ60XQB8Aye7/VgGXB2RAwpfu935N5zvj4c7Tn+c+DtETGyuMr59qJN\nNSQiLgQ+B7w3M1s7rZoLXFLM8jmF0kQlD1ND7/mjVn8fRcS7KH2K3QB8MzOvq3JI6iYRcR7wa+Ap\n9t5b82eU7qu6FTiJ0h/jD2bmuuKP8dco3bDaClyemY/0eODqNhExG/hvmXlRREyldOVqFPAY8J8z\nc0dENAL/Qumeu3XAJZn5QrVi1rGLiJmUJicZCLwAXE7pQz/P9z4uIv4c+BClYUCPAX9A6X4Jz/k+\nJCK+C8wGxgCrKM3i9385ynM8Ij5G6f0AwHWZ+a2e/D50dA6S9y8Ag4C1RbeHMvOTRf9rKN1ntZvS\nrR8/Ldpr4j1/zRZVkiRJktQb1OrwP0mSJEnqFSyqJEmSJKkMFlWSJEmSVAaLKkmSJEkqg0WVJEmS\nJJXBokqS1KtFxDURsSAinoyIxyPiDT18/NkR8aOePKYkqbb0P3wXSZKqIyLeCFwEvLZ4VtEYSs+z\nkiSp1/BKlSSpNxsPvJKZOwAy85XMXB4RsyLi3oiYHxE/j4jxABFxSkT8MiKeiIhHI2JalPxVRDwd\nEU9FxIeKvrMjYl5E3BYRz0bEd4oHjxIRFxZt9wPv7wgmIt5SXC17PCIei4hhPf9fIknqbXz4rySp\n14qIJuB+YAjwS+D7wIPAvcCczFxTFEnvyMyPRcRvgOsz886IaKT04eE7gU8CFwJjgN8CbwCmAz8A\nZgDLgQeAPwUeARYBbwUWF8cckpkXRcQPi/0/UMS2PTN398T/hSSp9/JKlSSp18rMLcAs4ApgDaUC\n5xPA6cBdEfE48EVgYnHVaEJm3llsuz0zW4HzgO9mZltmrqJUkL2uOMTDmdmSme3A48Bk4FXAksxc\nlKVPHv+1U0gPAF+NiE8DIyyoJEngPVWSpF4uM9uAecC8iHgK+BSwIDPf2LlfRBx3kF3EIXa/o9Pr\nNvb+XexyGEdmXh8RPwbeBTwUEW/LzGcP/11Ikvoyr1RJknqtiJgeEc2dmmYCC4GxxSQWRMSAiJiR\nmZuAloi4uGgfFBFDgPuAD0VEQ0SMBd4MPHyIwz4LTImIacXypZ3imZaZT2XmlykNE3xVN32rkqQa\nZlElSerNmoCbI+KZiHgSOA3478AHgC9HxBOUhu2dU/T/CPDpou+DwAnAncCTwBPAPcDVmbnyYAfM\nzO2Uhhv+uJioYmmn1Z8tJrx4AtgG/LT7vlVJUq1yogpJkiRJKoNXqiRJkiSpDBZVkiRJklQGiypJ\nkiRJKoNFlSRJkiSVwaJKkiRJkspgUSVJkiRJZbCokiRJkqQyWFRJkiRJUhksqiRJkiSpDBZVkiRJ\nklQGiypJkiRJKoNFlSRJkiSVwaJKkiRJkspgUSVJkiRJZbCokiRJkqQy9K92AOUYMWJEnnLKKdUO\nQz1s69atDB06tNphqArMfX0y7/XJvNcvc1+femve58+f/0pmjj1cv5ouqsaNG8cjjzxS7TDUw+bN\nm8fs2bOrHYaqwNzXJ/Nen8x7/TL39am35j0ilh5JP4f/SZIkSVIZKlZURcSkiPhVRCyMiAUR8Zmi\n/X9ExMsR8Xjx712dtvlCRCyOiOci4h2Vik2SJEmSukslh//tBq7KzEcjYhgwPyLuKtbdkJlf6dw5\nIk4DLgFmACcCv4yI/5CZbRWMUZIkSZLKUrGiKjNXACuK15sjYiEw4RCbzAG+l5k7gCURsRh4PfDv\nlYpRkiRJUnl27dpFS0sL27dvP+Z9DB8+nIULF3ZjVEensbGRiRMnMmDAgGPaPjKzm0Pq4iARk4H7\ngNOBPwE+CmwCHqF0NWt9RHwNeCgz/7XY5hvATzPztv32dQVwBcDYsWNn3XrrrRWPX73Lli1baGpq\nqnYYqgJzX5/Me30y7/XL3NeepqYmxo0bx/Dhw4mIY9pHW1sbDQ0N3RzZkclMNm7cyKpVq9iyZcs+\n684///z5mXnW4fZR8dn/IqIJuB34bGZuiogbgS8BWXz938DHgK4ycEDFl5k3ATcBTJ8+PXvjLCGq\nrN46O4wqz9zXJ/Nen8x7/TL3tWfhwoVMnDjxmAsqgM2bNzNs2LBujOroDBs2jC1btnDWWYetn7pU\n0dn/ImIApYLqO5l5B0BmrsrMtsxsB/6Z0hA/gBZgUqfNJwLLKxmfJEmSpPKVU1D1BuXGX8nZ/wL4\nBrAwM7/aqX18p27vA54uXs8FLomIQRExBWgGHq5UfJIkSZLUHSp5pepc4CPAW/ebPv1/RcRTEfEk\ncD7wxwCZuQC4FXgG+BnwKWf+kyRJknQ4LS0tzJkzh+bmZqZOncqVV17Jjh07WLt2Leeffz5NTU1c\neeWVFTt+xYqqzLw/MyMzX52ZM4t/P8nMj2TmGUX7e4tZAju2uS4zp2Xm9Mz8aaVikyRJktQ3ZCbv\nf//7ufjii1m0aBGLFi1i27ZtXH311TQ2NvKlL32Jr3zlK4ffURkqPlGFJEmSpPrwnr+7/5i2O9zs\nfz/8/8476Lp77rmHxsZGLr/8cgAaGhq44YYbOPnkk7nuuus477zzWLx48THFdaQqOlGFJEmSJFXS\nggULmDVr1j5txx13HJMnT654MdXBK1WSJEmSusWhrigdSjlTqmdml7P39cTzeDt4pUqSJElSzZox\nYwaPPPLIPm2bNm1i1apVTJ8+vUdisKiSJEmSVLMuuOACWltbueWWW4DS/VlXXXUVV155JYMHD+6R\nGCyqJEmSJNWsiODOO+/ktttuo7m5mdGjR9OvXz+uueYaACZPnsyf/Mmf8O1vf5uJEyfyzDPPdHsM\n3lMlSZIkqaZNmjSJuXPnAvDggw9y6aWXMn/+fGbNmsWLL75Y8eNbVEmSJEnqM8455xyWLl3ao8d0\n+J8kSZIklcGiSpIkSZLKYFElSZIkSWWwqJIkSZKkMlhUSZIkSVIZLKokSZIk1bSWlhbmzJlDc3Mz\nU6dO5corr2THjh3cddddzJo1izPOOINZs2Zxzz33VOT4FlWSJEmSalZm8v73v5+LL76YRYsWsWjR\nIrZt28bVV1/NmDFj+OEPf8hTTz3FzTffzEc+8pGKxOBzqiRJkiR1j396yzFtNqS9Hfod4nrPJ+49\n6Kp77rmHxsZGLr/8cgAaGhq44YYbOPnkk7nuuutoamoCYMaMGWzfvp0dO3YwaNCgY4rzYLxSJUmS\nJKlmLViwgFmzZu3TdtxxxzF58mQWL168p+3222/nzDPP7PaCCrxSJUmSJKm7HOKK0qG0bt7MsGHD\njmnbzCQiumzvsGDBAj73uc/xi1/84piOcTheqZIkSZJUs2bMmMEjjzyyT9umTZtYtWoV06dPp6Wl\nhfe9733ccsstTJs2rSIxWFRJkiRJqlkXXHABra2t3HLLLQC0tbVx1VVX7ZkB8N3vfjd/+Zd/ybnn\nnluxGCyqJEmSJNWsiODOO+/ktttuo7m5mdGjR9OvXz+uueYavva1r7F48WK+9KUvMXPmTGbOnMnq\n1au7PYaK3VMVEZOAW4ATgHbgpsz8m4gYBXwfmAy8CPx+Zq6P0kDIvwHeBbQCH83MRysVnyRJkqS+\nYdKkScydOxeABx98kEsvvZT58+fzxS9+kS9+8YsVP34lr1TtBq7KzFOBs4FPRcRpwOeBuzOzGbi7\nWAZ4J9Bc/LsCuLGCsUmSJEnqg8455xyWLl16wIyAlVSxoiozV3RcacrMzcBCYAIwB7i56HYzcHHx\neg5wS5Y8BIyIiPGVik+SJEmSukN0nmqwYgeJmAzcB5wOLMvMEZ3Wrc/MkRHxI+D6zLy/aL8b+Fxm\nPrLfvq6gdCWLsWPHzrr11lsrHr96ly1btux5iJvqi7mvT+a9Ppn3+mXua8/w4cOZNm1al9OaH6m2\ntjYaGhq6Maqjk5k8//zzbNy4cZ/2888/f35mnnW47Sv+nKqIaAJuBz6bmZsO8Z/d1YoDKr7MvAm4\nCWD69Ok5e/bsbopUtWLevHmY9/pk7uuTea9P5r1+mfvas2TJEnbu3Mno0aOPubDaXMZzqsqVmaxd\nu5YRI0Zw5plnHtM+KlpURcQASgXVdzLzjqJ5VUSMz8wVxfC+juk3WoBJnTafCCyvZHySJEmSyjNx\n4kRaWlpYs2bNMe9j+/btNDY2dmNUR6exsZGJEyce8/aVnP0vgG8ACzPzq51WzQUuA64vvv6gU/uV\nEfE94A3AxsxcUan4JEmSJJVvwIABTJkypax9zJs375ivEvUGlbxSdS7wEeCpiHi8aPszSsXUrRHx\ncWAZ8MFi3U8oTae+mNKU6pdXMDZJkiRJ6hYVK6qKCScONqjygi76J/CpSsUjSZIkSZVQyedUSZIk\nSVKfZ1ElSZIkSWWwqJIkSZKkMlhUSZIkSVIZLKokSZIkqQwWVZIkSZJUBosqSZIkSSqDRZUkSZIk\nlcGiSpIkSZLKYFElSZIkSWWwqJIkSZKkMlhUSZIkSVIZLKokSZIkqQwWVZIkSZJUBosqSZIkSSqD\nRZUkSZIklcGiSpIkSZLKYFElSZIkSWWwqJIkSZKkMlhUSZIkSVIZLKokSZIkqQwVK6oi4psRsToi\nnu7U9j8i4uWIeLz4965O674QEYsj4rmIeEel4pIkSZKk7lTJK1XfBi7sov2GzJxZ/PsJQEScBlwC\nzCi2+YeIaKhgbJIkSZLULSpWVGXmfcC6I+w+B/heZu7IzCXAYuD1lYpNkiRJkrpL/yoc88qI+C/A\nI8BVmbkemAA81KlPS9F2gIi4ArgCYOzYscybN6+y0arX2bJli3mvU+a+Ppn3+mTe65e5r0+1nvee\nLqpuBL4EZPH1fwMfA6KLvtnVDjLzJuAmgOnTp+fs2bMrEqh6r3nz5mHe65O5r0/mvT6Z9/pl7utT\nree9R2f/y8xVmdmWme3AP7N3iF8LMKlT14nA8p6MTZIkSZKORY8WVRExvtPi+4COmQHnApdExKCI\nmAI0Aw/3ZGySJEmSdCwqNvwvIr4LzAbGREQLcC0wOyJmUhra9yLwCYDMXBARtwLPALuBT2VmW6Vi\nkyRJkqTuUrGiKjMv7aL5G4fofx1wXaXikSRJkqRK6NHhf5IkSZLU1xyyqIqIs3oqEEmSJEmqRYe7\nUvXPEbEoIv4iIk7rkYgkSZIkqYYcsqjKzDOBi4A24LaIeDwiPhcRJ/dIdJIkSZLUyx32nqrMfC4z\n/zwzTwMuA0YA90TEAxWPTpIkSZJ6uSOeqCIi+gHHA+OAocCaSgUlSZIkSbXisFOqR8SbgEuBiyk9\nrPd7wB9n5sYKxyZJkiRJvd4hi6qIeAlYRqmQ+vPMXNUjUUmSJElSjTjclarzMnNpRLwOuDEiJhfb\nBJCZ+eoKxydJkiRJvdohi6rMXFq8/FfgT4GngPZKByVJkiRJteKw91QV1mTm3IpGIkmSJEk16EiL\nqmsj4uvA3cCOjsbMvKMiUUmSJElSjTjSoupy4FXAAPYO/0vAokqSJElSXTvSouo1mXlGRSORJEmS\npBp0pA//fSgiTqtoJJIkSZJUg470StV5wGURsYTSPVVOqS5JkiRJHHlRdWFFo5AkSZKkGnVERVWn\n51VJkiRJkjo50nuqJEmSJEldsKiSJEmSpDJYVEmSJElSGSpWVEXENyNidUQ83altVETcFRGLiq8j\ni/aIiL+NiMUR8WREvLZScUmSJElSd6rklapvc+CsgZ8H7s7MZuDuYhngnUBz8e8K4MYKxiVJkiRJ\n3aZiRVVm3ges2695DnBz8fpm4OJO7bdkyUPAiIgYX6nYJEmSJKm7RGZWbucRk4EfZebpxfKGzBzR\naf36zBwZET8Crs/M+4v2u4HPZeYjXezzCkpXsxg7duysW2+9tWLxq3fasmULTU1N1Q5DVWDu65N5\nr0/mvX6Z+/rUW/N+/vnnz8/Msw7X70gf/ltp0UVbl9VeZt4E3AQwffr0nD17dgXDUm80b948zHt9\nMvf1ybzXJ/Nev8x9far1vPf07H+rOob1FV9XF+0twKRO/SYCy3s4NkmSJEk6aj19pWoucBlwffH1\nB53ar4yI7wFvADZm5ooejk3qFi9v2MaTL23grMmj6N8v+M2StZw2fjgnjR7CsrWtPPnyBgBmnDic\nTdt28dL6VgDa2pP1W3cyZtigPfsaOrA//foFm7fvOuxxW3e2sbstOW7wvqd1xz5eWLOF0U2D6Fdc\nF26I4HVTRjGmaVAXe5MkSdKRqlhRFRHfBWYDYyKiBbiWUjF1a0R8HFgGfLDo/hPgXcBioBW4vFJx\nqfZs3bF7T+EB8NLmdp5duemIt39m+SZ+vmAlyzdsZ8KIwby8YRsAk8cM7fZYAV58ZWvx6vl92ieP\nGdppXS8x7/mK/T8ciZ2723jvayYw7fgji+FIc3/SqCEMGdhbRjdLkqS+rmLvOjLz0oOsuqCLvgl8\nqlKxqPLa27Prm+A6efyl9dy/aO1R7/uXC1fts7x58w5uffHJo94PsKegAnq8wDnU8d51xnh+8tSK\ngy53tB3MK1t28PCS0mSbJ40awukThgPwswUraW/fNzPvOmM8qzZtZ/7S9YeNqyf8473PH75T4Why\n/7ZTxx1rSEdsyMAGPjBrIscNHlDxY+0vgH79urodVZIk9TQ/ytVhbd/VxoPPv8KLr7Ry/HEHDhVb\ns3kHdzz6csXjaD6+iYhgZW7lhHHDjmrb363azOimgexqa2fowP6MGDKAK948jUq9Jx3Yvx87d7cf\n8BrghOGNBMGKjduICCaNHEz/hn5cfu5k1mzewcSRg4kIPvmWqSzfuJ22tmTSqFLboWxo3cmutmRs\np+GDfzh7Gi9v2EZbW6mwGjd8EIP6N+zpv27rzu7+1o/Yyo3b+f4jL9G/35Hf2nm43Gcmi1ZvAQ4s\nxitl7hPVu/1zzswTOWF44zFtu3rTDo4/bhCnjj+OaWN732xLkiTVEouqOrSxdRcrNm07oP2uBat4\nduXmAz79PporGYcrUtqz9EZ/QMPRzZHSfHzTnmFq8+atZ/bs1xzV9r3R1P3eyDYOaGDSqCF7liOC\nCSMGH/H+RgwZ2GX7wfYxYsjAg27TE6aObeKcU8Yc1TZHkvslr2xlcVFYVdLDS9by0AvrKlaYH0rH\nxccfPN49BV13DAHNTI4f1siaLTsAGNAv+K9vnsphPgs4Ikcy7HPy6KE0Dmgo/2CSJB0Di6o+Zv9h\neGu37OCOx15mx67SlZIdu9v49aJXyjpGV8PQ+vcL3jHjBE4aPaSLLaSeM2XMUKb0wH1iv3da5YcX\nHsxL61r52dMr2d1+bM8Z3H9oaXcNAV26tnWf5atvO7Zhuvs70mGflR7yOX54I+9/7YTDXjUuV7+g\n4seQJHUvi6oatXn7Lh58fi272vYOK/v359fyZMvGI97Hf+hiGFVDP7ocFjegoR/tmUwaOcT7OKQq\nmzRqCP/1zVOPefuPnzeFlRu3M3bYIFZsPPCq9bH4l4eWsmnbbjKTF17ZytSxQ4kuH0F49A417HP7\nrjaWrSsVcz0x5PNfHlpa8WMA/MGbptDQC37XDh7QwLmnjPEqoCQdhkVVL9WyvpWtO9r2LN/87y+y\nefvuPcuH+2S589/i9oS3vur4PZMXAJwxYfgx34shqbYN7N9vz1Xl/YehHqtr3zOjW/bTlcMN+1y0\najMv7neVrLt95zdLWbtlZ8WHe3ZcfPz6r5dU9kBH4a9/uagqs4SuXrWd21c8VrH9d/47OmHEYM45\nZTRvmDK6Ysc7Wk2N/Y9qCLik6rKoqqKOoXqPLlvPg4v3zor37y+8sk9BdSidZ3oDaBrUwMVnTmBY\nY8/PRiZJ1dA8bhjNRzl5zdHqqeGej7+0gX9//uhnSa2EBxa/wsZtpWfkVWOW0M2tSWsPHfflDdv4\nt0da+LdHWnrkeEfq2vecxlmTR1U7DElHwKKqB23bWZpFb/uudl5cu5WfPb3ysNs0jys+RU4YP6KR\n//jaiXvWDWscsM9Mb5Kk2jZz0ghmThpR7TCA0qRCKzdup3Xn7sN3roCHH36Y179+ZkX2vXLjdv7P\nw8tYuraVE0c0snzD9r1/b3uBNZt3sKF1Fys3ba92KJKOkEVVBby0rpXWnW1s2r6Lr//6BQYWU1gf\n7JO+flEa8vFHs6fRv5gVb1D/frx+yijHsUuSqqaaw8SXDe3XbcNT93css4/2pH+893l+/OSKw3eU\n1GtYVJUhM/eMf7/vd2t4smUj9y9ew/Zd7YfcbuTQgZw9dRQDG/rx7lePZ/xwx0xLkiRJtcqi6ii0\ntycPLVnL+q272NXWzjfuP/SNxM3jmtjdlsyePnbPcI6hg/oz7jgniJAkSZL6CouqQ1i1aTsbWnfx\n8wUreW7VZpYdZHapjtmg2hOufOspDB7QwBumjmJQf4fuSZIkSX2dRVWhYya+x5at54HFa/nd6oMX\nUQDvPOMEguBNp87OUwAAEoxJREFUzWP2mX1PkiSpO7Qd4wO+JfW8ui6qXlrXypMtGw87E1/zuCYa\nIvjEW6YysKGBiSMH+wBcSZJUUV//9RJ+9vRKLnr1iQAMGdjAG6eNdhIrqReqq6KqYzjfLxeu4pnl\nm1i27sArUR0z8f3h7GkMHtDAWZNH+swnSZLUY149cfie2f9a1m/jH+99fu/Ku2DKmKHsamvnTc1j\nmXXySCJg8uihDOzfr0oRS+rTRVVmkgkrNm3nr+/6Hc+u3Nxlv7edOo6hgxqciU+SJFXdOdPGcMcf\nncPt81tY17oTgHufW0PrzjYAlhSPaPnuw8v47sPL9mz3tlNLD6nu3xC8+4zxnDRqCICja6Qe0CeL\nqsWrt/DMik38830vdLm+eVwT/SL4xJunMn7EYJoG9cn/BkmSVKMGNPTjkteftGf5j2afwoqN29i2\ns43N23fzrQeW0K9fsLF1F6s37wDglwtX7enf+baGoYMa+MjZk4v9Bm+cNtpROFI36xPVxO62dpa8\nspUfPL6cpetaD3jIbseQvjkzT+RDr5vkLxJJklRzOo+m+etLztzzesHyjby8fhsAyzds4/ZHX97z\n3gdg6462fYYQ/t09i5kyZiht7cnMSSN4y/SxAJw0aoj3a0nHqOaLql8sWMnf3bO4y3XvPOMEzp46\nmteeNLKHo5IkSeoZM04czowT985E/NFzpwCl2yDmPrGclzeUCq4Fyzftmdm4YwjhsnWtzH1i+Z5t\nO4YQAgxr7M/vv24SQ4pCKwIiHEoodaWmi6qVW9v3KaiGDmpg6tgmPn7eFCaMGOynLZIkqW5FBHNm\nTtinbfXm7WzZvptdbck/FVevdra1s7QotjoPIQS487GX91n+2HmT9zyHc0BDP944bbS3UUjUeFHV\ncVm7X8BfzDmd10waUd2AJEmSerHjhzVy/LDS669+aOae9sWrN/PCmr23T9z7uzU82bKRjjkuOt5z\nffP+F/fZ39/evYgpY4buWR7W2J/LzpkMwJimQYwaOrDbvwepN6rpogrgrMkjufY9M6odhiRJUs06\n5fhhnNJRbQFvn3HCPuuffnkj9y1aQ+be5ZbiPq4l+93LftWtT+x5PevkkXsKq/4NwUVnnMjEkfvO\ntOzshOoLqlJURcSLwGagDdidmWdFxCjg+8Bk4EXg9zNzfTXikyRJ0l6nTxjO6ROG79PWMZSww23z\nW1i5cTsAi1ZvAWD+0n3fyv30qZXs76RRQ3j3q8fvWV6zqa3b4pZ6SjWvVJ2fma90Wv48cHdmXh8R\nny+WP1ed0CRJknQonYcSAlx94av2vH5lyw4eW7aB9uLS1kvrWvnB48vpfFGqY0jhsnWt3Dhv7+yE\nmzfv5L51v2Xofvdq/eHsaTQUO+gXASQnjfKhx+odetPwvznA7OL1zcA8LKokSZJqzpimQfzeaeP2\nafuDN03dZ3nn7nZum9/C+uIBx7D3+VqrN++A4vlbHa6+7ckuj9VxnIZ++z70uDOHGKrSIjsGx/bk\nQSOWAOuBBP4pM2+KiA2ZOaJTn/WZecBc6BFxBXAFwOBxU2Z9+Nob+fCpg3oqdPUCW7Zsoampqdph\nqArMfX0y7/XJvNen3e3J0rVbGTpkb2H06OrdvLS5naBUGC3f0n7U+z17fH9GD95bWG3amRw3MDht\ndANDB1hw9Qa99Zw///zz52fmWYfrV60rVedm5vKIOB64KyKePdINM/Mm4CaAYROac+LEicye7UQV\n9WTevHnMnj272mGoCsx9fTLv9cm816/9c3/Jfut3t7Xz8JJ1vLi2ldFNpUkwVnR66HFnHUMMF2wB\nthx4rPvWsGf2wh2721i+YTsTRw7mo+dMZuR+MxeOH97IsMYBZXxnOpRaP+erUlRl5vLi6+qIuBN4\nPbAqIsZn5oqIGA+srkZskiRJ6r36N/TjnFPGcM4p+7Z3PPS4s5fWtfLTp1ewq23vyKyOIYYd9p+9\nsGX9Nv7njxd2eezzX3U8/TtVbis2bmPH7nbe85oTeUvz2AP6+8Dk+tHjRVVEDAX6Zebm4vXbgb8A\n5gKXAdcXX3/Q07FJkiSp75g0aghXvHnaPm1/NHvanungd7btHUr4oydW7Hn4cfPx+w5D65jN8FfP\ndv2Z/1d/8Tu++ovfdbnuP752Ascf17hP24CGfpw9dZRXvvqQalypGgfcWVTt/YH/k5k/i4jfArdG\nxMeBZcAHqxCbJEmS+rCIYFIXk1l85m3NfOZtzV1u07K+lQXLN+3TtrF1F//y0NJO+913m45pC25/\n9OUu9/m3d7PPg5M76xfwydnTaDjEVa4JIwczZGBvmnOuvvV4JjLzBeA1XbSvBS7o6XgkSZKkQ5k4\ncggTRx5YiP3+6yYddJtla0tDD3e37zsp3DPLN7FsXStw4NDDzv7037qe7bCzt5067oBi7tFl65k5\naQQfPGsS4/e7QtbBYYndz/JWkiRJ6mYnjR7CJ94yrct1qzdvZ3OnByd3dtczq1i4YtMhr1J1DEfs\nGK64v7sXrubuhYeenuCj50xmyMCGg66fOraJ6ScMO+h67cuiSpIkSepB+z84ubNpbzn8tOJL127l\n2ZWbD2j/5TOr9rQfrCbrGJb47QdfPOxxJowYfNCHKy/fUJqkY8qYofzh7Gn7TOBxMJNGDaFxwMEL\nuVpmUSVJkiTVkJNHD+Xk0Qfej/WOGSccdtuFKzbxq+dWc6hH1XbMkPjyhm2H3d+SV7Ye9MHMXXn7\nfg+F7rD8pV3MPuK99D4WVZIkSVKdOHX8cZw6/rhD9vmDN03ZM0PiwXzrgSU88dJG4MDZEve3qz15\nsbh/7BfPdD1kMXe0HXIfvZ1FlSRJkqQ9BvVvYNrYQxdK//PiM45qn4tXb2bx6oNPzPG7hU8f1f56\nG4sqSZIkSRV1yvHDOOVgN5IBja8824PRdL+u7zyTJEmSJB0RiypJkiRJKoNFlSRJkiSVwaJKkiRJ\nkspgUSVJkiRJZbCokiRJkqQyWFRJkiRJUhksqiRJkiSpDBZVkiRJklQGiypJkiRJKoNFlSRJkiSV\nwaJKkiRJkspgUSVJkiRJZbCokiRJkqQyWFRJkiRJUhl6XVEVERdGxHMRsTgiPl/teCRJkiTpUPpX\nO4DOIqIB+Hvg94AW4LcRMTczn+mq/0B2MWbbi7CmV30bqrAhW1+CNc9VOwxVgbmvT+a9Ppn3+mXu\n69Pg1uXVDqEsva0aeT2wODNfAIiI7wFzgC6LqvG5mve9+BdwR2MPhqhqO3XzZmgZVu0wVAXmvj6Z\n9/pk3uuXua9PzdsD+E/VDuOY9baiagLwUqflFuANnTtExBXAFQDNJzSxdeDxrGqLnotQVbdr4FBa\n23rbj656grmvT+a9Ppn3+mXu69O2/gN5ct68aodxzHrbT2xX1VHus5B5E3ATwPTp07P5sz/sibjU\ni8ybN4/Zs2dXOwxVgbmvT+a9Ppn3+mXu61Ot5723TVTRAkzqtDwRqO0BlpIkSZL6tMjMw/fqIRHR\nH/gdcAHwMvBb4D9l5oKD9N8MeCdj/RkDvFLtIFQV5r4+mff6ZN7rl7mvT7017ydn5tjDdepVw/8y\nc3dEXAn8HGgAvnmwgqrwXGae1TPRqbeIiEfMe30y9/XJvNcn816/zH19qvW896qiCiAzfwL8pNpx\nSJIkSdKR6G33VEmSJElSTan1ouqmagegqjDv9cvc1yfzXp/Me/0y9/WppvPeqyaqkCRJkqRaU+tX\nqiRJkiSpqiyqJEmSJKkMNVtURcSFEfFcRCyOiM9XOx51n4iYFBG/ioiFEbEgIj5TtI+KiLsiYlHx\ndWTRHhHxt8XPwpMR8drqfgcqR0Q0RMRjEfGjYnlKRPymyPv3I2Jg0T6oWF5crJ9czbh17CJiRETc\nFhHPFuf9Gz3f60NE/HHxe/7piPhuRDR6zvc9EfHNiFgdEU93ajvqczwiLiv6L4qIy6rxvejIHSTv\nf1X8rn8yIu6MiBGd1n2hyPtzEfGOTu018Z6/JouqiGgA/h54J3AacGlEnFbdqNSNdgNXZeapwNnA\np4r8fh64OzObgbuLZSj9HDQX/64Abuz5kNWNPgMs7LT8ZeCGIu/rgY8X7R8H1mfmKcANRT/Vpr8B\nfpaZrwJeQyn/nu99XERMAD4NnJWZp1N6PuUleM73Rd8GLtyv7ajO8YgYBVwLvAF4PXBtRyGmXuvb\nHJj3u4DTM/PVwO+ALwAU7/MuAWYU2/xD8SFrzbznr8miitLJtDgzX8jMncD3gDlVjkndJDNXZOaj\nxevNlN5gTaCU45uLbjcDFxev5wC3ZMlDwIiIGN/DYasbRMRE4N3A14vlAN4K3FZ02T/vHT8PtwEX\nFP1VQyLiOODNwDcAMnNnZm7A871e9AcGR0R/YAiwAs/5Picz7wPW7dd8tOf4O4C7MnNdZq6n9OZ8\n/zfs6kW6yntm/iIzdxeLDwETi9dzgO9l5o7MXAIspvR+v2be89dqUTUBeKnTckvRpj6mGN5xJvAb\nYFxmroBS4QUcX3Tz56Hv+GvgaqC9WB4NbOj0C7hzbvfkvVi/seiv2jIVWAN8qxj2+fWIGIrne5+X\nmS8DXwGWUSqmNgLz8ZyvF0d7jnvu9z0fA35avK75vNdqUdXVJ1PODd/HREQTcDvw2czcdKiuXbT5\n81BjIuIiYHVmzu/c3EXXPIJ1qh39gdcCN2bmmcBW9g4D6op57yOKoVtzgCnAicBQSkN89uc5X18O\nlmfz34dExDWUbvf4TkdTF91qKu+1WlS1AJM6LU8EllcpFlVARAygVFB9JzPvKJpXdQzzKb6uLtr9\neegbzgXeGxEvUrq8/1ZKV65GFEODYN/c7sl7sX44Bw4vUe/XArRk5m+K5dsoFVme733f24Almbkm\nM3cBdwDn4DlfL472HPfc7yOKSUYuAj6cex+YW/N5r9Wi6rdAczFD0EBKN7bNrXJM6ibFGPlvAAsz\n86udVs0FOmb7uQz4Qaf2/1LMGHQ2sLFjSIFqR2Z+ITMnZuZkSuf0PZn5YeBXwAeKbvvnvePn4QNF\n/1756ZUOLjNXAi9FxPSi6QLgGTzf68Ey4OyIGFL83u/Ived8fTjac/znwNsjYmRxlfPtRZtqSERc\nCHwOeG9mtnZaNRe4pJjlcwqliUoepobe80et/j6KiHdR+hS7AfhmZl5X5ZDUTSLiPODXwFPsvbfm\nzyjdV3UrcBKlP8YfzMx1xR/jr1G6YbUVuDwzH+nxwNVtImI28N8y86KImErpytUo4DHgP2fmjoho\nBP6F0j1364BLMvOFasWsYxcRMylNTjIQeAG4nNKHfp7vfVxE/DnwIUrDgB4D/oDS/RKe831IRHwX\nmA2MAVZRmsXv/3KU53hEfIzS+wGA6zLzWz35fejoHCTvXwAGAWuLbg9l5ieL/tdQus9qN6VbP35a\ntNfEe/6aLaokSZIkqTeo1eF/kiRJktQrWFRJkiRJUhksqiRJkiSpDBZVkiRJklQGiypJkiRJKoNF\nlSSpV4uIayJiQUQ8GRGPR8Qbevj4syPiRz15TElSbel/+C6SJFVHRLwRuAh4bfGsojGUnmclSVKv\n4ZUqSVJvNh54JTN3AGTmK5m5PCJmRcS9ETE/In4eEeMBIuKUiPhlRDwREY9GxLQo+auIeDoinoqI\nDxV9Z0fEvIi4LSKejYjvFA8eJSIuLNruB97fEUxEvKW4WvZ4RDwWEcN6/r9EktTb+PBfSVKvFRFN\nwP3AEOCXwPeBB4F7gTmZuaYokt6RmR+LiN8A12fmnRHRSOnDw3cCnwQuBMYAvwXeAEwHfgDMAJYD\nDwB/CjwCLALeCiwujjkkMy+KiB8W+3+giG17Zu7uif8LSVLv5ZUqSVKvlZlbgFnAFcAaSgXOJ4DT\ngbsi4nHgi8DE4qrRhMy8s9h2e2a2AucB383MtsxcRakge11xiIczsyUz24HHgcnAq4AlmbkoS588\n/munkB4AvhoRnwZGWFBJksB7qiRJvVxmtgHzgHkR8RTwKWBBZr6xc7+IOO4gu4hD7H5Hp9dt7P27\n2OUwjsy8PiJ+DLwLeCgi3paZzx7+u5Ak9WVeqZIk9VoRMT0imjs1zQQWAmOLSSyIiAERMSMzNwEt\nEXFx0T4oIoYA9wEfioiGiBgLvBl4+BCHfRaYEhHTiuVLO8UzLTOfyswvUxom+Kpu+lYlSTXMokqS\n1Js1ATdHxDMR8SRwGvDfgQ8AX46IJygN2zun6P8R4NNF3weBE4A7gSeBJ4B7gKszc+XBDpiZ2ykN\nN/xxMVHF0k6rP1tMePEEsA34afd9q5KkWuVEFZIkSZJUBq9USZIkSVIZLKokSZIkqQwWVZIkSZJU\nBosqSZIkSSqDRZUkSZIklcGiSpIkSZLKYFElSZIkSWX4f9XNMRfi23TgAAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" }, { "name": "stdout", "output_type": "stream", "text": [ "TCLab disconnected successfully.\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAA1UAAAGoCAYAAACubM8oAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4wLCBo\ndHRwOi8vbWF0cGxvdGxpYi5vcmcvpW3flQAAIABJREFUeJzs3Xd81eXd//HXdU7GyQ4JSQiEvRGR\nJSgggmLdVtE66vpVW6qtVettW+3dVupdW9u7t3a32mkdVdxWq9UKceBkyQ4zQEIge+9zrt8f35M9\nyD4Z7+fjkUfOd39OviGcz/e6rs9lrLWIiIiIiIhI17gCHYCIiIiIiMhApqRKRERERESkG5RUiYiI\niIiIdIOSKhERERERkW5QUiUiIiIiItINSqpERERERES6QUmViIjIAGCM8RhjrDEmpYfP+7Qx5ns9\neU4RkaFGSZWIyABmjClt9OUzxlQ0Wr420PF1hzHmmDFmSR9fc5gx5tfGmCP+n+E+Y8zPjTFx3Txv\nryREza7xkTHmut46v4iItE1JlYjIAGatjaz7Ag4DFzda92Sg42uLMSaov13DGBMGpAITgRVANLAE\nKAfm9nR8rVy/138mIiLSO5RUiYgMYsYYtzHm+8aYA8aYXGPMk8aYWP+2acaYWmPMzcaYTGNMnjHm\nJmPM6caY7caYQmPMQ43OdYsxZq0x5hFjTLExZqcxZmmj7XHGmL/7W5iOGGPuM8a4mh37W2NMAXCP\n//qpxph8Y0yOMeYxY0yUf/9ngUTgTX+L0e3GmPOMMfuavb/61ixjzIPGmKeMMc8YY0qAq9t7/624\nGYgHLrfWpllrfdbaY9baH1hr/+O/xsnGmPf8P5utxpjzG8XytDHmF8aYfxtjSowx640xY/2b3/V/\nT/O/n0vr3o8/vuPA7/3n+boxZr//frxgjEnqyr1vFFeQMeZ5Y8xxf9zrjDFTm+2W5F9fYox52xgz\nqjvXFBEZapRUiYgMbt8CPofT4pIC1AAPN9ruBmYBE4AvAb8G7gbO9K//kjFmYaP9lwKf4SQfDwIv\nGWOi/dueBIr851oAXApc3+zYLcBw4P/86+4HRgAnA1OB/waw1n4ByAY+5291+1UH3+/lwGNADPB8\nB95/YyuA16y1Fa1tNMZ4gFeBl4AE/7mfNcaMb7TbF4F7gTggC/hho/cOMNX/fl7yL48DgoHRwO3G\nmAuA7wOXAaOAXOCJDr739ryC0wI3AtiN8zNq7Hrgu/73tbeV7SIi0g4lVSIig9tXgXustUettZU4\nH/KvMsaYRvvcb62tsta+4l/+u7U2z1p7GPgAmNNo3yPW2t9Za2ustX8HMoBz/S0yS4G7rLXl1tos\n4FfA1Y2OPWCt/aO11mutrbDW7rbWrrXWVltrjwG/wEnmuuMda+2//K1MFR18/3XicRKhtpzh//6Q\n//3/G3gLuKrRPmustZustTXAU8DsE8RbBfyP/2dQAVwLPGqt3eqP99vA2caYESc4T5ustbXW2ses\ntaWNfgYL/ElinZestR/6t3/Xf82Erl5TRGSoUf9tEZFByp84jAb+ZYyxjTa5cBIIAK+1Nq/Rtgrg\neLPlyEbLGc0ucwgYCYwFPEBOo3zFBTTurnekWXwjgV8Ci4Ao//7tJTUdUX+NDrz/3GbH5gHJ7Zx7\nJHDYWtv4XIdwWpTqHGv0upymP7vWHPMnYI2vsbZuwVpbaIwp9l+j8ATnapV/rNaDwEqcVkIfYHB+\nBpn+3ep/btbafGNMqT+WnK5cU0RkqFFLlYjIIOX/8J8JnGWtjW305bHWNk8oOqp59boxwFGcD+Wl\nwLBG14m21jYu8GCbHfu/QBkw01obDXwZ58N+W/uXAeF1C8aYYJxudo3VH9OF9/8f4IJmLTiNHfW/\n38bG0JCYtKf5e2lr/VGcBBUAY0wMTsGMjlyjLV/C6QK5HKdb5LS60zfaZ3Sja8bhJIPdTXBFRIYM\nJVUiIoPbH4AHjTGjAYwxicaYi7txvtH+ohNB/vLdY4A3rbUHgY+AnxljoowxLmPMZNN+SfQonESs\n2BgzBrir2fbjOOOz6uwC4owxZ/sTqh9y4v/HOvP+/wzk44yTmmIcCf6CG2cD7wEuY8yd/vd/Dk6y\n8uwJYsBaW0XDeLP2/AP4ijFmpj+5+ymw1t89siOCjVO+ve4rCOfnXInTEhcB/KiV4z5vjFlojAn1\nb19nrc3u4DVFRIY8JVUiIoPbz3BaYNb6K+J9QPfKg7+LM8YqH6eoxGXW2iL/tmuAWJxCCPnAM0B7\nlet+gFNAogh4EaewRGMPAA/4K9bd5m9dugOnIEYGTle7E7W4dfj9+8c0LcPp0rcWKAE+xElENvnH\nG10EXIGToDwEXGWt3X+CGBq/32f97+eSNmJ4FfgJTmGJoziFJa5vbd82/AWny2bd1x9wksUcnJ/X\nNuD9Vo57AqeLYC4wHbixE9cUERnyTNOu4SIiIq0zxtwCXGGtXRHoWERERPoTtVSJiIiIiIh0Q0CS\nKmNMujFmmzFmizFmg3/dauNMPrnF/3VBIGITERERERHpjIB0/zPGpAPzG1dfMsasBkqttT/v84BE\nRERERES6SN3/REREREREuiFQk/9a4E3/ZIyPWGsf9a+/zRhzA7AB+C9rbUHzA40xq4BVAB6PZ96Y\nMc2nDJHBzufz4XLpecBQpHs/NOm+D02670OX7v3Q1F/v+549e3KttQkn2i9Q3f9GWmuPGmMSgbeA\nbwBpOKVcLfA/QLK19qb2zjN16lSblpbW6/FK/5KamsqyZcsCHYYEgO790KT7PjTpvg9duvdDU3+9\n78aYjdba+SfaLyDpoLX2qP97Ns7cJAustcettV5rrQ/4I7AgELGJiIiIiIh0Rp8nVcaYCGNMVN1r\nnNnotxtjkhvtdhmwva9jExERERER6axAjKlKAl40xtRd/ylr7RvGmMeNMbNxuv+lA18NQGwiIiIi\nIiKd0udJlbX2AHBKK+uv7+tYRERERESk+2pqasjIyKCysrJLx8fExLBr164ejqrjPB4PKSkpBAcH\nd+n4QFX/ExERERGRQSIjI4OoqCjGjRuHv0dap5SUlBAVFdULkZ2YtZa8vDwyMjIYP358l87R/+oW\nioiIiIjIgFJZWUl8fHyXEqpAM8YQHx/f5VY2UFIlIiIiIiI9YCAmVHW6G7uSKhERERERkW7QmCoR\nERERERnQ8vLyOPvsswE4duwYbrebhIQEAObOncurr75KYmIi27f3zqxNSqpERERERGRAi4+PZ8uW\nLQCsXr2ayMhI7r77bgDeffddbrvtNm644YZeu76SKhGRIcLrszz58SEmJUQyLCKEFzdnUl5diyfI\njfXvU1XrJSzYzednj2LmqJhWz5NVVMHjHx6iqtYHgMtAebUXT7Cbqlpvk/PVXTci1M0ZkxPYkJ7P\nFxeOJS4ipHffbAC8vi2LihovK+emBDoUERFpZOnSpaSnp/fqNZRUiYh0Q3Wtj/05pbiMYWJCBEHu\nvhmqWlheTY3Xkl1SibUn3h9gZ1Yxz27IAGB0XBhH8iva3PdwfjnfOGtyq9ue35TBhvSCTsf87p5c\nADILK7l24ZgOH3eo2Mv2zKJOX68rYsODqaj2UlXrIzzEjdtlMBiqvT4qa7xtHldcWcPvUvcDMCw8\nhISo0A5dz+0yDAsPIbe0qssxu4xhQkIExsD+7DKC3IYJwyM4XlJFkP/8Pms5kFPGmPhwIkP1X7+I\n9K6Lf/1+p4/xer243e529/nnN5Z0NaRep7+sIiLd8PvU/fxn13EAVs4dxZcWd21+i84oqazha09u\noqSytsvnaJ5Qfe/C6fz0jd3UeJ0M7WhhJfe+sK3N410uw3+dM4XQIBc/es2ZrHF4ZAi5pdUA3Hv+\nNNwuw+H8cv7+4aEmx27PLGr33M2VlFTz3KGO7x9oD721p8+vOX/cMEKD3Kzf5ySuxlCfbC+fmkBx\nZS0bDxUwNj6cX109B5dr4FboEhHpj5RUiYg08tmRQp7bmIG3g80/O48W179+bWsWe7NLeyu0esUV\nNU0SqhnJ0R0+tqSqhqjQYFwuWDxpODuPFnPhrGROGhnDg5fP4oVNmZySEsO7e3Pw+do+z6JJ8Syd\n4gwAXn3JDN5Jy+GmJeN56pPDTE6MYtGk4QCcOi6O4spaNh8uIMoTREFZDTFhnZutPsuUkTyi4++x\nqwrKq8kqcuYoGRUbRmZh08RzQkIEnqC2n6LWeH14rW13n+Z2Zjm/PyFBLiYlRHYhatibXcKG9AKM\naUimGv/6rkvLqX99KK+c7zy/leCg/lf8d3x8BDcvGa+ET2QQ6EqLUiAn/+0JSqpEZMA7nFdOTmnX\nJ+xr7I/vHmzxYfpElkweTlFFDdsyitiW0Tfd1OpcMS+FGxeN6/LxF80aWf96SlIU95w/DYDzT07u\n8DnmjY1j3tg4AL62bFKTbS6X4eYl44Gut+ClpuazbNmsLh/fUceLK/n6k5sYGRvGz66YxS1PbCTP\n3/KWHOPhoStn4+7hD/wPvZnGurQcvnfhdOaMGdalczzyzn5e3ZqFtbBiehIZBeXsPlbS5v7tbQuk\nbRlFjI4Lq+86ubfAS9ShfDzBbiJCghgR48ET3PGEtT3HiyuprvWRXeL83YiLCKWwvJq4CKe1NcTt\nYsbI6B6/3yIyeCmpEpEB7VBeGbf/YzO+Do4r6oj4yBDuXDGFjnyccrsMkxIj8VnLvuzSDo9v6q4o\nTxBBLhcpw8L65oJDQFK0h99fNw9PsAtPsJtfXjUHn7V4rSUkyNUrH7BvO2syVy0Yw6jYrt/Hm5aM\nZ7G/ZXBKUhQ1Xh8Hcsrqk5Pc0iqMgcmJUezPKaW6tp0myAD5ND2fl7cc5bfr9tevKymp5pWMnfXL\nS6cM51vnTuv2tXJLq/j6k5vqC6205frTxnLlqaO7fT0RCbxrrrmG1NRUcnNzSUlJ4Yc//CE333xz\nj15DSZWI9BuZhRX8bt0+5o4ZxoZDBcweHcP6fXlkFJSTHBtGdnYlz2ZuanJMcWUNPusUXkiM8nQ7\nBmPgolnJzB4d2+ljZ6V0/hjpXxoXmIgJ71w3xa4ICXJ1K6ECCHa7mlRqDAlycXJKw/KImIZ/F9M7\n0VW0L01IiODlLUcBCHIbTkmJJSOjjLjEKHZlOS1r7+7JJT13U3un6ZCy6tomCVV0WBDFFS3HJz6z\n4Qjv7Mlpsb4tI2I8fPu8qQS7XPz+nf2Eh7hPOMbyqY8Pc7y4ktvPnoy1lv97aw+H88qb7FNaVUuw\n20VVrZdLThnJF+Yr0RM5kdWrVzdZ/sc//tHr11RSJSIBsftYcYtiCS9syiCjoIKt/i50jSu+Hc4r\np6TcUpnf9AMHQLDb8L0LZzCymx9ORSQwojzBXD53FC9uzuS+i09i9uhYUlNzOGPpLG7/x2YO+//d\nH27l3393fPmM8UxPjubbz20lJiyYgvJqrj51DJsPF7D7WEmnrnc4v5y/rk8n2hPMG9uPARATFkyU\np/XkvLLGyz8+OQw4VScB3t+b2+41/v7hISJCgwhuo8poyrCwfps4iwx2AUmqjDHpQAngBWqttfON\nMXHAM8A4IB240lrb+Zq9ItLvZRVV8J3ntnaqy96wiBCunxjKaQvnttgWExbcJ60KItJ7bjh9HCvn\npRDdKAlxuwwPXXUKXp8lr7S6x7rXRnmCcBlDdFgQxhgev3kB4SFBlFbVEu0J4op5KRwr6vg4zS0Z\nhfzx3QO8tjWryfq/rk/v0PEvbMqsf339aWM5bUI8AMdLKrn/nzub7Pv71P20xWXgt9fOJWVYeAcj\nF5GeEsiWquXW2saPZO4B3rbWPmiMuce//J3AhCYi3bUuLZsnPzqMr5VPQRXVXnwWxg+PYEJCRJNt\nxj+SyfiHryyeFM/Oo8WcOSWRA9s+YUy8PiyIDEYul2mSUNUJ9VdTDI/rvY8sda1JdZUpQ4JMp/7W\npAwLo7Sytr7whcFgDK3+/WssyGWobfR0KS4ihMvmjqpviRoTH87tZ09mf04pPmvx+WyT/RtLzy1j\nf04Z33p2K2EhLQt6jB8ewb3nT+uzufREhpr+1P3v88Ay/+vHgFSUVIn0a2VVtWw5UsjpE+LZfrSI\nQ43GAjy7MYOCsuo2jw1yG755zhTGD49oc586dZXlDnQ/ZBGRHudyGb7YiQmtO+OcGUmcQ9IJ9zuc\nV86dz2ymtKqW0qqWY8RySqr48/sHT9hNOjEqlIX+ljIR6Thj+6pUVeOLGnMQKAAs8Ii19lFjTKG1\nNrbRPgXW2hb1ZY0xq4BVAAkJCfPWrFnTV2FLP1FaWkpkZNfmc5GetSatmp15XuYmudmc7W3RNSch\n3HDttJBWj/UEGTxBnaumpns/NOm+D026751XWWuprG35uW5Hnpe3DnV8svAbTwphfEzPlK/vCt37\ngSkmJoZJkyadeMc2eL1e3O7A/d4B7Nu3j6KiplOjLF++fKO1dv6Jjg1UUjXSWnvUGJMIvAV8A3il\nI0lVY1OnTrVpaWm9HK30N6mpqSxbtizQYQwp2SWV/OClHWQWVhDsbkiEarxN/35MSYpicpLzH6Hb\nGM6ansjELk5o2hrd+6FJ931o0n3vObVeH89vyqCgvKbd/TIKyvnsSBEuQ/0UAm7/XHPnzez43HXd\npXs/MO3atYvp06d3+fj+MPlva+/BGNOhpCog3f+stUf937ONMS8CC4Djxphka22WMSYZyA5EbCJD\nlbWW1LQcxsSHsz2ziOKKhv98P8soqp8Qt3kiVSci1M2dKyYzOk5jnkRE+pMgt4urTj1x98Si8hru\neGYzeaXV+Px/62u8lic+OkxOSRUAnmA3F80a2eq4LZFAysvL4+yzzwbg2LFjuN1uEhISKCkpYcyY\nMRw7dgyXy8WqVau44447evz6fZ5UGWMiAJe1tsT/+nPA/cArwI3Ag/7vL/d1bCJDSY3X12Sulo2H\nCnjorT0nPO7608Zy6ZxR9cvBbkON1xLkMrh6YXJUERHpGzHhwfzlxlPri2FYLN96disHc8tYsyGj\nfr+iihquXjAGtzFKrqTfiI+PZ8uWLYAzT1VkZCR33303WVlZZGVlMXfuXEpKSpg3bx7nnHMOM2bM\n6NHrB6KlKgl40TilvYKAp6y1bxhjPgXWGGNuBg4DXwhAbCJDQmZhBd98ZgsV1d429zn3pKQmE6FG\ne4IJcrs4a1pifbeQOiGdHBslIiL9k8tlCGn0N/47509j/d5cLJbyai8vbMrk5S1H6ydrrpMYFcov\nrp7d5rxcIoGSnJxMcrLTfTUqKorp06eTmZk58JMqa+0B4JRW1ucBZ/d1PCKDlbWWVz47SnZxVYtt\ne7NLqKj2Euw2TSaRLPevmzEymluXTWqRPImIyNAyKjaMK08dXb9cUF7DxwfysNDkwVx2SRU//3da\nq3NkjYkP59yTRvRFuNJfPHJmpw8J9/nAdYKS/199p4sBOdLT09m8eTMLFy7s1nla059KqotINxWW\nV9d329iWUcSf3jvY7v4PXTmbcR0oaS4iIgJw1zlTAOfB3a1PbKofbwuw6XAhmw4Xtnrc8MgQxsY7\n/9/EhYeou7j0udLSUi6//HJ+8YtfEB0d3ePnV1IlMkis+fQIj390qMX65VMTmJjYsgJfyrBwJVQi\nItIlxhjuv/Qk9h0vpdrrw1oormxZXXB7ZhEfHchn9Ss769edNDKan6w8GWOUWA1aXWhRKu/F6n81\nNTVcfvnlXHvttaxcubJXrqGkSmQA2pddwitbjta3SoFTaAIgLiKEuv+nkqI8fPXMiUSE6p+6iIj0\nrMQoD4lRnnb3WTo5gdzSnRSUO5PBF1XUsONoMQ+8touQoIauXnPGDOOcGSee5Fiks6y13HzzzUyf\nPp277rqr166jT1oi/VRuaRWF/v+Emvv12n0cyClrsX7aiCj+9wsthiyKiIgExLCIEB6+anb98lMf\nH+Yfnxzm44P5TfZbvy+XhKhQIkPdlFb3/RyqMnitX7+exx9/nJNPPpnZs53fxR//+MdccMEFPXod\nJVUi/dCR/HJuf3oztW3MCQUQExbMl88YX999wgCzUmL6KEIREZHO+8L8FCYkRDSZ0mPtruNsOlzI\n91/aDkB1eSWLF1czLCIkUGHKALd69er610uWLMHa3k/UlVSJ9BPPb8xgx9FiAI4XV1LrtcRHhhAb\n1rI8rcsYLpk9kmVTE/s6TBERkS4Ldrs4bUJ8k3WTEiP5zdq9VFR7ySur5kgJ/PCfO4iLaJjWIznG\nw01LxqsqrfRbSqpEAiC3tIpjRZX1y/ll1fztg/Qm+7gM3H/JTMbEtyxPKyIiMliMig3jJytnAbA1\no5Db//4h+3PK2N+sm/uwiBCmjXAKGYweFk5MuObEkv5DSZVIHyuprOEbT22mtKq2xbYzJg/nzCkJ\nACRGe5RQiYjIkDIrJZavnBzKpJOm16/bmVXMC5syeazRw8fhkSE8cv38JsUuJPCstQO2qmN3uwgq\nqRLpRT6f5Tfr9pGe2/C0ray6ltKqWuIiQhgZ21A1KTwkiJuXjCc+MrS1U4mIiAwJo6JcLGzURXD2\nmFgKyqrJKXUmsz+SX0FuaTV3PrMZT5C7fj+Xy3DVqaM5dVxcn8cs4PF4yMvLIz4+fsAlVtZa8vLy\n8Hjar2bZHiVVIr1g59Fi8sqqyCyo4K2dx1vd584Vk5kzZlgfRyYiIjKwhAa5uetzU+uXU9Oy+b83\n93Akv6LFvo+8c4DKGm+L9UEuF3ERIQyPDNHDy16SkpJCRkYGOTk5XTq+srKyW0lNd3k8HlJSUrp8\nvJIqkR62PbOIe1/Y1mTdFfNSWDSx4albRGgQI2PD+jo0ERGRAe/MKQlMTIhskjxZ4Kev7+Z4cSU/\neyOtzWOTYzz8/rp5KnjRC4KDgxk/fnyXj09NTWXOnDk9GFHfUlIl0kG/XbePDelN59XILa1meGTT\nkq9l1c4f+clJkSRGeYgND+aqU0fjCXYjIiIi3WOMYXRcyzHHd6yYzBvbj+FrNjTGYvlgXx4AWUWV\n3PzYp1w2ZxSfnz2qL8KVIUJJlcgJHC+uZP2+XN7YfqzV7bmlLSfojQwN4gcXzSA2XHNsiIiI9IVZ\nKbHMSoltddtDb+1h3e5sAPJKq3nio0N4gt24mo39SRkWxvTk6F6PVQYfJVUi7fD5LP/94naOFzvl\nz8+bOYIvzHf62/709TT2HC/h1HFx3LJsQpPjoj3BapkSERHpJ+48ezK3LZ9ERbWXB9/YxfbMYn6z\ndl+L/VwGfnvtXFKGqfqudE7AkipjjBvYAGRaay8yxvwNOBMo8u/y/6y1WwIVnwxd1lp+8vputmUU\n4bWWimovwyJCWDwxnmsWjiHa48yLce8F03htaxaXzR1Vv05ERET6H5fLEOIyhAS5+NqySbzy2VFq\nvU37CabnlbEvu5Q7n95CsNvFsIhg7v/8TIarsIV0QCBbqu4AdgGN21i/Za19LkDxyBC3fl8umQUV\nFFXU8OH+vPr1xsCNp4/l7OlJTfYfHhnKjYvG9XGUIiIi0h2j48L5+vJJLdZnFJRz1zOfUVHjparW\nR2lVLb9Zu48Z/u6ALpfhjMnDSYoOXIU66b8CklQZY1KAC4EHgLsCEYMION37fNaSdryEB1/f3WTb\n1QtGc8kpIwlyuQgLUVc+ERGRwSxlWDhPfHkhVbVe0nPL+e6L29h4qICNhwrq9/n4QB4/WXly/bIx\nRpUEBQDT3dmDu3RRY54DfgJEAXc36v53OlAFvA3cY62tauXYVcAqgISEhHlr1qzps7ilfygtLSUy\nMrLb58mt8PHnbVVU1DasmxjrIjnChScIFiYHEaw/lP1KT917GVh034cm3fehq7/c+605tWSXN3xO\n/vRYLVUtp8BiyaggVozVMIDu6i/3vbnly5dvtNbOP9F+fZ5UGWMuAi6w1n7NGLOMhqQqGTgGhACP\nAvuttfe3d66pU6fatLS25yKQwSk1NZVly5Z1+fi9x0v46EAeO44Ws+NoMS4DGENCZAg/vXyWJgXs\nx7p772Vg0n0fmnTfh67+eu9f2pzJ3z9Mp/FQLJ/PEuw2XDpnFM0fw05MiGTRpOF9GeKA1l/vuzGm\nQ0lVILr/LQYuMcZcAHiAaGPME9ba6/zbq4wxfwXuDkBsMkj5fJay6lp8Fn702i7yy5wy6EZVfkRE\nRKQDLp0zikvnNJ3bavUrO9h4qIBnN2S0eszDV80mKbrhYW14SJC6Cw5SfZ5UWWvvBe4FaNRSdZ0x\nJtlam2WMMcClwPa+jk0GJ5/Pcvezn7E3u7R+XVK0h8/NSGLc8AglVCIiItIlt501iXfScvA2m3F4\nS0Yh2zKK+OYzTQtZJ0aF8usvziE8RLMaDTb96Y4+aYxJAAywBbglwPHIAGat5aUtmRzOq6Csupa9\n2aW4XIawYBfBbhc3LRnHoolqkhcREZGuGx4ZyuXzUlqsXzJ5OPe9soOSypr6dVW1PrJLqvjxv3aR\nENmyguDUEZGcNzO5V+OV3hPQpMpamwqk+l+fFchYZPCoqPbySXo+f3k/vcn6mxaP4/OzR7V+kIiI\niEgPGRkbxh9vaDoM58P9efz4X7v47EgRDdOyNvjPruMkRHlIjGroLhgfGaJWrQFCd0kGlapaL19/\nahM5JU7hyGVTEzh5VAxhIW61TImIiEjAnDYhju9dOJ2iipoW2z4+mM8nB/NZ/cqOJutjwoJ55Pp5\nRITqI3t/pzskg0JmYQXPfHqE3NIqckqqCAtxMykxkq8snUC0R2VORUREJLCMMSycEN/qttljYimu\nqKG0qmGel4LyaooqavifV3cyLCIEgEUT4zljckKfxCudo6RKBqxar4/0vHKstTzx0SE2HS6s33br\nmRNZPi0xgNGJiIiIdExilIf//cIpTdbVdRfccbS4ft2nB/OJiwghxO1qcY7RceF4gt29Hqu0TkmV\nDFi/Xbef/+w6Xr8c5DbcftZkosOCmTsmNoCRiYiIiHTPaRPi+J9LZ1Ls7y74z8+OsvtYCfc8v63V\n/WeOiuEnK0/uyxClESVVMuAcKfHx/Ze2szWzCJdxJtcDWDYtUa1TIiIiMigYY5g9uuEh8chYD396\n7yDVtb4W+x7KL2d7ZhHffu6zVufBCgsO4pZlE0iMall1UHqGkioZMHw+y7bMIv51sJpS43T1Wzpl\nON86d1qAIxMRERHpXZMSo3hax5TgAAAgAElEQVTw8lmtbnvio0M88+kRdmWVtHl88HuGs6cnAWAM\nTEmMIiZc4857ipIqGTBe2pLJX9enU1JqGTHczfcunMHkpMhAhyUiIiISUNcsGMO8scOo9doW24or\na3jw9d18sD+PD/bn1a8fNzyCX109G2NatmxJ5ympkn7rqY8P8+6enPrl3FKnTPqEGBdfPmsyM0fF\nBCo0ERERkX7D7TJMT45uc/vNS8bzWUZDQa/dWSWk55ax6vGNuFtJqsJC3MwdO4yP9ucxe3QsX1k6\noVfiHkyUVEm/4vNZNh0uoKC8hmc+PYyv2QOXiQkRfD6phkWTNOeUiIiISEdcOmcUl84ZVb/88pZM\n/vTeQY4VVbZ5zL7sUgAO55eTHOtpdYqahKhQQoNcTEhQzyElVdKv/HvHMX6Xur9+edHEeK47bWz9\nclK0hw/efzcQoYmIiIgMCpecMpLTJsS3WvSitKqWbz+3tcm6R9450O75Hr5qNpMSh3ZipaRK+oU/\nvLOfjw7k1ZcNnTMmlviIUK5ZMJrEaFWqEREREekpxhiS2vl8dcfZk9lypJBFk+LZmF5ARY23xT7v\n7c2tf33vC1uJCG2aVsweHUt2SRULx8fx+dmjmh8+6CipkoDKLKxg/d5cXtuaVb9uTFw49118Uqsl\nQUVERESkd62YkcSKGU6lwEUTWx9yccZkZ3JigMoaH5U11U22v70rG4BtGUUtjg0NcrNsasKgmqxY\nSZUEjM9nuf+fOzha6PTnPfekJK5eMIbYsGAlVCIiIiL92OkT43n+1kX4rKW0qrbJtsc/PMTa3dn1\ny39672CL47OKKvjS4vG9HmdfCVhSZYxxAxuATGvtRcaY8cDTQBywCbjeWlvd3jlkYDqSX853X9xG\ncUUNPgvDIkJYMT2RS+eManUQpIiIiIj0PyFBLoAWLU43LR5PXEQIXn/FsVpfw9itqhofb+48zgub\nMnlpy9H69baykmXLej/m3hLIlqo7gF1AXf3HnwIPW2ufNsb8AbgZ+H2ggpOeZ61lXVo2b+3MprDc\nGTvlchluOG1sfROziIiIiAxsMeHB3LhoXJvbiytr+OhAPr5GZZ5tyym2BpROJ1XGmElAkrV2fbP1\nZwBHrbX7Wz+yyb4pwIXAA8Bdxpl17Czgi/5dHgNWo6Rq0PD5LJ+k5/PwW3vr1/3u2rmMig3Dpa5+\nIiIiIkPGf184g1pv08qD77zzToCi6Rldaan6BfDdVtZX+Ldd3MFzfBuI8i/HA4XW2roOmRnA4C8T\nMkRUVHu585nN9WOnTpsQx/KpiYyOCw9wZCIiIiISCEFuV5PlgT6evitJ1Thr7dbmK621G4wx4050\nsDHmIiDbWrvRGLOsbnUru7baCGiMWQWsAkhISCA1NbVjUUufy63wseGYl7xKH3sLfBgDw8MMc0Nr\nqM7IJjWja+ctLS3VfR+idO+HJt33oUn3fejSvR+aBvp970pS1d6kQWEdOH4xcIkx5gL/uaJxWq5i\njTFB/taqFOBoawdbax8FHgWYOnWqXTaQR7QNYl6f5fZ/bOZwaTkAUVFw7/nTWDSp9bKcnZGamoru\n+9Ckez806b4PTbrvQ5fu/dA00O+768S7tPCpMeYrzVcaY24GNp7oYGvtvdbaFGvtOOBqYK219lpg\nHXCFf7cbgZe7EJv0A9ZavvXcZxzOLycs2M1Xlk7gO+dN4/SJ8YEOTURERESkx3WlpepO4EVjzLU0\nJFHzgRDgsm7E8h3gaWPMj4DNwJ+7cS4JgMLyah774BBFFTXsPV4KwPWnj+XiU0YGODIRERERkd7T\n6aTKWnscWGSMWQ7M9K9+zVq7tgvnSgVS/a8PAAs6ew7pH/LLqnnio0P8Z9fx+nVfXDhGCZWIiIiI\nDHpdnqfKWrsOp8ueDHG7soq55/mt1E01cMuZE0mICmXumNjABiYiIiIi0gcCOfmvDHA1Xh9/fO8A\nmw4V4LMQHxnCkknDuXBWcqBDExERERHpM0qqpEuOF1fy9q5sXt92DIDQIBcPXzmbYREhAY5MRERE\nRKRvKamSTssrreLrT26iqtaZCfvqBaNZPjVRCZWIiIiIDEldTqqMMSW0nKC3CNgA/Je/8IQMMn95\n/yAf7M+lqtZHYlQos0fHcuX80QS7u1KdX0RERERk4OtOS9VDOBP0PgUYnDmnRgBpwF+AZd0NTvqP\n/LJqPk3P58XNmQC4DNx97lSmJ0cHODIRERERkcDqTlJ1nrV2YaPlR40xH1lr7zfGfLe7gUn/4fNZ\n7nl+K1lFlQCsmJ7ENQtGkxjtCXBkIiIiIiKB152kymeMuRJ4zr98RaNtzbsFygCUUVDOQ2/uobiy\nhuPFVUSHBTFn9DBuXDSW2HCNnxIRERERge4lVdcCvwR+h5NEfQRcZ4wJA27rgdgkQKy1fHwwn1e3\nHmVvdmn9+utPG8t5M1UuXURERESkse5M/nsAuLiNze939bwSeB/sz+PB13cDztipH688mYSoUBKj\n1N1PRERERKS57lT/mwL8Hkiy1s40xswCLrHW/qjHopM+VVJZw3+/uJ3D+eUAzB0Ty7JpiZw0MibA\nkYmIiIiI9F/dqYP9R+BeoAbAWrsVpwKgDEBbjhTyyDsHOJhbhtdnSYr2cO8F01k+NTHQoYmIiIiI\n9GvdGVMVbq39xBjTeF1tN+ORAMgsrOC+l7fj85cXuef8aZw+IR6Xy7R/oIiIiIiIdCupyjXGTMRf\n6c8YcwWQ1SNRSZ85kl/O157cBMDUEVGcMXk4iybG0yxZFhERERGRNnQnqfo68CgwzRiTCRwErjvR\nQcYYD/AuEOq//nPW2vuMMX8DzgSK/Lv+P2vtlm7EJ+2w1vLatizeScupX3fH2ZMZHRcewKhERERE\nRAae7lb/W2GMiQBc1tqSDh5aBZxlrS01xgQD7xtjXvdv+5a19rl2jpUeUOv18d7eXB555wAAxsAj\n188jOSYswJGJiIiIiAw8nU6qjDF3tbEeAGvtQ+0db621QN3kR8H+L00W3EdKq2q57alN5JVWA7Bs\nagLLpyUqoRIRERER6SLj5DidOMCY+/wvpwKnAq/4ly8G3rXWfrkD53ADG4FJwG+ttd/xd/87Hacl\n623gHmttVSvHrgJWASQkJMxbs2ZNp+IfyjYdr2Vrrpf0Ih9BLkiOcHHdjBBC3QNr/FRpaSmRkZGB\nDkMCQPd+aNJ9H5p034cu3fuhqb/e9+XLl2+01s4/0X6dTqrqDzTmTeDyum5/xpgo4Flr7XmdOEcs\n8CLwDSAPOAaE4IzV2m+tvb+946dOnWrT0tK6FP9Qs/d4CXet+ax++b6LZzB/XFwAI+q61NRUli1b\nFugwJAB074cm3fehSfd96NK9H5r66303xnQoqerOPFVjgOpGy9XAuM6cwFpbCKQC51lrs6yjCvgr\nsKAbsUkjh/LK6hOqeWOHsfqSgZtQiYiIiIj0N92p/vc48Ikx5kWcMVGXAY+d6CBjTAJQY60tNMaE\nASuAnxpjkq21WcYZnHUpsL0bsQng81me+PgQHx/Ir1/3laUTGBWr8VMiIp3i88Kmx6DwCCRMhVlX\nOVV+Gqsuhw1/gUkrIHGasy7tDTjyMVgvjJoP5bnOOQB8NeAKhonLIWMDVJeCt6bhvO5QmHQ2HP4Q\nxixyvs+9AcJ74aFYUSZsW+NcJ+MTmPclqKmAzY/DKddAdHLPX1NEZBDpTvW/B/xV+87wr/qStXZz\nBw5NBh7zj6tyAWusta8aY9b6Ey4DbAFu6WpsAkXlNby7N4dnN2QA4HYZfnftXEYqoRKR3lZdDt5q\nCIuF2mrIPwCuIGc5NAqMG/L3O4lKB0SUHoLjO3s56BPI2gIb/c8N96+FsGEQM7rpPjtfhj1vwM6X\n4OJfOUlJ6k8ath94p/Vz71/b9nX3vul83/GS870oA+bf1LX30J7UHzvJXt11KgqgOAuyd8KxbbD0\nWz1/zTrGwLBxUFUKnhgoyQJ3MGHlWb1z311uiEhwXpdmQ8RwKMuF8HjnfcdPdNZHjWiZOJfnQ3AY\nuEOcY6KSWp6/NMf5Xbc+53ffWucaFYXgq23YLzy+6fElx524XN3pRCQigdKdliqstZuATZ08Zisw\np5X1Z3UnFmlQVlXL157aSHGF88f7wlnJfG5GkhIqEel9Ph+88g0oPQ5feAw++BUcSG3YnnIqRCbC\n7tc6fMppJSWQGdXzsXbHuh+3vc1bAy/d2vVzB4c5CVlrMj51vnrbvrcbXucf6N776QjjTyRcQU5C\nDswoKYEjAbjv7hAnhsV3wMyVDesLj8DzX3ZaKodPhm3PwXk/gbGLGvY5ugVe/SZMOddJwI583PZ1\nXEFwxZ+dhDJ9Pfz7u3DyF2DRbb321kSk93QrqZL+p6Syhq89uYniilpiw4OZkRzNjaePIyzEHejQ\nRKS3leXCuz+Hshz/B/NyCO6BCb0rCpyWmY7wVkPePuf1P2+H4kznaX9dUaSMT50P0MblfDjtgDLf\nMaISR3Qh8B5WVQIhEVBb2fbPtSwXIpOc7n4AQR4YfwbsedP52USPhMTpsO8/zn6jF8C+tVCeByHh\n4ImFJXfCnn/7W8TinPsZMdzZp6P3odOMc53IJKelKDLRWV1R6LS69Kb8g87PFOoTqnrxk8Ad3HPX\nshZydre/T10M63/ZtBWxLNeJM+sz5wvgjXthxMkN++Ttd1qo0vzTb7pDICjU+d0BSJjm/HuoKICS\nY/D6Pc69PbbN2b7t2RPH11/FjoUz7nK6su57CxbdDp7oQEcl0meUVA0iR/LLeeWzoxSW1wBw54rJ\nzBurghQi/Vp5PhQecj5sWZ/zgaquKGtIhPMBrKq03VPU2/O6M+6mNxRldP2Y6Rc7H5yP+4fKWh9M\nPAtW3Nf2sY3sTk1lRD+sCNUpMy9vujznura3ASz4ivM1FGz8G2z4a9N1xkV5+EiiVv6x57vDvXVf\nQ+upy+10QzUu5/ey8QMAaEh22tPePidd5oyZe/nrMPo0OM/fwllwCJ67yUlgS7I6fr7+7Ng256HB\nJ486y9YHUy/s0qmiivdARidaKY3xJ6wuZ9xiTEqXrivSHUqqBon03DLueHozPv//BV89c4ISKpH+\nzlsLL33N+VA1fqlTqCCzUz2qT2zhLZA0o+vHv32/84R++JSOd0sK8jhflYVOF6fhU5wucVUlUJHv\nfIgdPqXrMcngMucGJ+HwRDstOFHJYC17PvyEpN4YX7T8v+G0W53fyehRTmtq1AinpS5smNPSGpHg\ndGFtLiwOqsucIiPhw6Esu+U+EQnOwxLjcn7P3UFw9T+atvgNGwtXPe5cs/64RCchsL6ef8+97cgn\nsPmJhoQKnC6kjbuRdsKUkhLI6mTXz9ELnVbNQx/ARQ/DyNldurZIVympGuB8Pssv3t7LZ0cK8VkY\nGethwfh4LjxZlZpEet2nf3I+TJTlNAxG78zcf96qhqfUB991vgd5nCSougxy/PPwDZ/sFHjoiBEn\nOx8Ws7Y655l1pfM0vqsu+LlTAe7Ur3ShAtzYhpfuYKd7W2sD+2Voc7kaqiVGj6xf7Q3qpXHAQSFO\nElUndnTTa9d154vqQJfTtv5NNHofzrla+b2PHtlyv4FaZTFhmvPwpSzbSSYriyG065O4Fh89StTI\nkSfesc7xHU3Hr71xj5Mw95TKQgiJhNoqpwdBnWAPnPFfEDeh564lA5aSqgHsSH457+zJYd1u50lZ\nsNuw+pKTSI5RQQoZxI5uaXi664lx/jPL2uJ/apzT/rE9qboMNj3edF1FQdfOFRzujH8CpyvYwlVO\n1bxnb3SqhV3ya2eMVCDEjYezfxCYa4vIwBAUCsvv7bHT7U1NZVRnuvx+9Hv47OmG5ZqKhrGdPaUs\nt/X17//C6eLcVUkzofRY0/MnzoCYHkwKpU8oqRqgyqpqufvZzyivdgZDX3nqaC46OZlhESEBjkyk\n90SUHoR//i3QYbQtagR87kedOybI4zytLkh3nvDG+lt3gkLg8j8DNnAJlYjIQLDgq07FRWudLpzl\neT137oxP4eNHGpZXrHbGbNVUwGt3NS1c0lPC4+Gap53/B2TAUFI1ABWV1/Clv31CjdcyIsbDwvFx\nXDk/hdAgVfiTfmTvf5xKVitWO11aPn7EKc3sDnG6anTB5NwsCA+BpJOc8xxtZWq8SSu6FXanBIfB\nuDMg/T1nedaVEDuma+eKn9hyXUgPVO4TERnsXK6mXfB6coLsuAlOl+qSY05X7InLG7YtuwcOdaM4\nUNaWhhaqpJnOg7ljW5150p68whmT2lhEAlzwM6eXhvQ7SqoGmIpqL99/eTs1XmfcxteXT2L26F4u\ndytSx+eDA+ucgdQn6ur24W+d7299H8afCVue6vbl3d4qCI13WoOCw2DNjc5g8rjxTkvP5x6AcYu7\nfZ1OG7Ow768pIiK9z+WGeTe2vm3Siu49yNvzJqx7wClicu6PnFa2A+/AWz+AyqKW+5fnwXsPOQ8W\nT8QT65zDHQyTP6eHdH1ASdUA8/hH6RzMLQPgexdOV0IlfWvrM/DxHzp3TO5e56uxU66Bk6/o/OXX\nr2fxWec5/fcBrn7SGdsUFus8SezJ+WxERER605TPwZjT/BVT/V39JpwJN/6z5Zxt2bvgze850wE0\nnlC9IwoPOZNZS69SUjVA7Msu5Uev7SS/zPlHtnLuKBaMV8l06ab8A/D6dzo+D1JdMQVwqj01nvSy\nOZfb6bpQW+Usu0OcFq6gUDjl6i6NE6oJiWlIqMBJourKFCuhEhGRgaa1CZJbWzduiVNpsPDwic+Z\n9i/ngWOdHS9C2hsdj8kY5/jmk5xPOdeZnFxapaRqADheXMl3nt9Kda0zd8WiifF8afH4AEcl/UZ1\nuTN7vbfGmeeoM4586vTd7ixPtDNxa/NywCIiItLzjIEZl3Rs3/FL4V/fgtNugYwNkP5+04eiHdX8\nmB0vOiXlO/UQ0zhTglSXOV0lB+q0AR2gpKqfq6718Yv/7KG61ofbZXjk+nkkRoWe+EAZmLy1Tgnt\nzlj/S9jTiSdQzRkXfOFvzhxLHREU5gwKFhERkf4neRbc/G/n9YxLO59QvXWfU/Vw5ko49cvOug9+\n47SAbX6i63Ed+diZ+7ANxlfT9XP3A0qq+rG3dh7nN2v34vPPJfrt86aSFO0JbFDSe7J3wT/vaOgu\n1xUTzmwoyd1RSSfBsE4eIyIiIv2fMU0nLO6Is77nFKWaemHDWK/TbnXmzursZ5RNf294fWwb/OXc\nNnc9udLAWed07vz9iJKqfqi4sobnN2awLi0HnwVPsIurTx3DookdbEmQnuXzwdannSpzvSlrq/PH\nyh3stB51Rm2VU471zHtU4UdERES6LiwWTrqs6TpPNMy5rvPnipvgTKkyYiakrwfrbXNXn8vX+fP3\nI32eVBljPMC7QKj/+s9Za+8zxowHngbigE3A9dba6rbPNHj96d0DrEvLASA5xsMfrpuHy2UCHNUQ\nU1XS0A3vyCdNJ/7rTe4Q+OIzPTvHhoiIiEggTFzedG6vdmxPTWVZ70bTqwLRUlUFnGWtLTXGBAPv\nG2NeB+4CHrbWPm2M+QNwM/D7AMQXUH96ryGhuuH0sSyeNFwJVV/b/CR88mjL9VMvcCb+600J05RQ\niYiIiAwwfZ5UWWstUFeiLNj/ZYGzgC/61z8GrGYIJVWH88p5aUsm69KcSmxnTkngC/NHBziqIWLH\ni854pjrp653vnmjAn9DGpMCib6hrnYiIiIi0YJwcp48vaowb2AhMAn4L/C/wkbV2kn/7aOB1a+3M\nVo5dBawCSEhImLdmzZo+i7u31Posv9lcRWGVcy+mxbm4epoq/LWltLSUyMjILh8fUpWH2z+pXkh1\nPpP2/bnlNSLHkTbt9i5fQ3pHd++9DEy670OT7vvQpXs/NPXX+758+fKN1tr5J9ovIIUqrLVeYLYx\nJhZ4EZje2m5tHPso8CjA1KlT7bJly3orzD5hreWbz2zBG1JGYqSbW8+cyPxxw4jyaCLTtqSmptLl\n+773LVj7q6broqJg/BkwdomzbFxEpcwnWd3w+p1u3XsZsHTfhybd96FL935oGuj3PaDV/6y1hcaY\nVOA0INYYE2StrQVSgKOBjK0vfLg/j39uPcr+HGfW65sWj2f5tMQARzXI7P4X7F/bsJy3z/kePaqh\nTGhoNCy6AyIT+j4+ERERERnwAlH9LwGo8SdUYcAK4KfAOuAKnAqANwIv93VsfeloYQU//lfDOJ7r\nTxvLeTNHBDCiQaQ0BwoPg68G3n8YvM2KSEYmwpWPdXJGcBERERGR1gWipSoZeMw/rsoFrLHWvmqM\n2Qk8bYz5EbAZaDnQZZAor67l609tAmBsfDi3nDmR6cnRAY5qkKgug+dvgsrihnUjZ8PsaxuW4ycp\noRIRERGRHhOI6n9bgTmtrD8ALOjrePra8xszWLs7m1qvM2Ts1mUTOWlkTICjGuDS3oDdrwIWqkqd\nhCpiOMSOdeZ9WrjKmXxORERERKQXBHRM1VDzwb5c/vZBev3y9y6croSqC4yvBg5/BN4asD5Y/0uo\nKW+605K7YNziwAQoIiIiIkOKkqo+cjC3jJ+8vhuA0yfGc93CsYyJ15xHXZGS8U/Y+1nTlUkzYeFX\nndehkWqZEhEREZE+o6Sql/l8lp/9O42tGYUABLsNX14ynsRoT4AjGyCsdYpN7HrFqdhnLQk5eyEy\nHMYuBmOc8VGzr4PhkwIdrYiIiIgMQUqqepG1ll+v3cf6fbkAhAa5+M0X5yqhao/PC/vXQY1TZp7y\nfNjpLwRZlAGAsV4YtwTOfSBAQYqIiIiINFBS1YtS03L4z67jAFy7cAwXnTKSyFD9yNu140X44Net\nb4tMhIseZvuHH3H6OSv7Ni4RERERkTboE34vKK6sYfXLO0jPc1pbpiRFsXJuCiFBrgBH1g9t/JvT\nEmWdaohUlzrfxy8Fj7+IhzsY3KEw7UKISaHKMxxc+lmKiIiISP+gpKqHVdZ4ue/lHezLdpKDpOhQ\nfrxyphKqOhWFsP9t8NY6lfs2PQ6+2qb7xE+EFT9U4iQiIiIiA4KSqh721MeH6xOqe8+fxoLxcQS5\nlRzU++BXsO/tpuvGLoaldzcse2KUUImIiIjIgKGkqgc9+fEhXtycCcDFpyRz+sR4jDEBjiqAdrwI\nH/ymZUuUccHMy/2V+0JgxqUQHheYGEVEREREuklJVQ84VlTJmzuP8ZI/oTp1XByrlk4McFQBsO9t\nyN3TsJz2esuECmDmSlh0W9/FJSIiIiLSi5RUdVN1rY/frNvLZ0eKAJg3dhg/uHhGgKPqQzWVgIWi\nTHj7/pbb4yfByj86rVJ1hnLrnYiIiIgMOkqquiE1LZuH39qDz0JIkIsbTh/LGZMTAh1W33n7/pbj\no0bNg5RTndfG5cwnpfFRIiIiIjKIKanqou2ZRfzfm05XN0+wiyvmpfD52aMCHFUv8nlh8xNQntuw\nvO9t/7ioUGedJxqWfBNiRwcuThERERGRPqakqpOstWSXVPHzN9MAiIsI4c83zh98Ff4qixvmjAI4\n9AFs+EvL/WZ83kmkRERERESGqD5Pqowxo4G/AyMAH/CotfaXxpjVwFeAHP+u37XW/quv4zuRP79/\nkJe3HK1ffuCymYMvocpJg5dudVqjmjvpMhg2znntDoEJy/owMBERERGR/icQLVW1wH9ZazcZY6KA\njcaYt/zbHrbW/jwAMXXIy1sy6xOqUbFhfGnxOFKGhQc4qh5wdLNT/tz6nOWCdCehCouFoLCG/YaN\ng9O+BkEhgYhSRERERKRf6vOkylqbBWT5X5cYY3YB/X4w0vbMIv703kEAFoyP4/sXDdAKf/kHobqs\n6bp3/heKM5uuc4fAyj9B5BAqvCEiIiIi0gXGWhu4ixszDngXmAncBfw/oBjYgNOaVdDKMauAVQAJ\nCQnz1qxZ0+tx5pT7+O2WKgDiPIZVs0LxBA28suCxBduYuP+vrW6rCo0jI+WS+uVKTyKVYSP6KrRO\nKS0tJTIyMtBhSADo3g9Nuu9Dk+770KV7PzT11/u+fPnyjdba+SfaL2BJlTEmEngHeMBa+4IxJgnI\nBSzwP0Cytfam9s4xdepUm5aW1msxWmv52wfpfHQgj6OFlQD84fp5jIoNO8GR/cT+dbD7NZwfKU63\nvrJcpxtfSKNfWpcLTrkGxi4KQJCdl5qayrJlywIdhgSA7v3QpPs+NOm+D12690NTf73vxpgOJVUB\nqf5njAkGngeetNa+AGCtPd5o+x+BVwMRW53c0io+3J/HC5ucbnEul+HXV8/p3wlV/kEnaQJnfNR7\n/wdVJU338UTDpb+HkEEwFkxEREREpB8IRPU/A/wZ2GWtfajR+mT/eCuAy4DtfR1bncoaL998ZguF\n5TUAXDgrmQtPTmZ0XD9ORPL2w/Nfbig2USdhKpz6lYblYWOVUImIiIiI9KBAtFQtBq4HthljtvjX\nfRe4xhgzG6evWjrw1QDERllVLbc+uYnC8hqGRYRwSkoMN5w+lvCQfjSl16bH4cC6pusqCp2EKnYM\nRCY661zBMO9GSJze9zGKiIiIiAwRgaj+9z7QWpWHgM9JtftYMW/uOE5BWTUAt545kdMnxgcmmNIc\np9R5c94qZxLe5i1SAEGhcP5PIXpk78cnIiIiIiJAgMZU9UdH8sv5znNb8flrOtwSyITK54N/3e0U\nlmjL+KUw98am68LjnC8REREREekzSqqAzMIKvvbkJgDGD49g8aR4zp/ZB+XEa6vhrR9A4aGm631e\nKD0OnhgYvbDlcUEhMOcGiErq/RhFRERERKRdQzqpOpRXxo6jxXywP7d+3d2fm8qY+B4u5JC92+ma\nV5YDxUcb1hcchMMftn3c3Bvg5Ct6NhYREREREelRQzapqqzx8t0Xt1FcUVu/7uGrZvd8QlWcBS/d\nAu3NB3b6bS3niHKHQGRCz8YiIiIiIiI9bkgmVX9df5BXt2ZRXesjKdrDnDGxnDQymkmJXZzFec+b\nsOkxp3z5oQ+bFpGw3jeutuIAACAASURBVKYJVfIpTlnzOuHDYeZKcLm7dm0REREREQmoIZNUeX2W\n17ZlUVRezUubM/FZZ0LfVUsnsGB8o+IO1sK+/0BNhbPsq4HaKqitbPvkGx9zvhdltB9E2DA4+wcQ\nMbx7b0ZERERERPqNAZ1UGSxUl59wP6+1/Hv7Mf6+/iAAwcBpE4Zzx4pJBLtcTc+R/j6se6DrQY2a\nB+c2O94VDFgwbnC5un5uERERERHpdwZ0UhVengF/Pb/dfWp8liP55Uz2WX4CRIcFE+w2RB8JJuix\n1qbLasOsqyA47P+3d+dxdtVlgv8/TypLJamQnRCSQBZqIgQ0GFQE1CC2oqJBR1sYx0G0B+2Wn9rN\nNGrjb5huh1fjtCO92E037QbdjkqzjHEXwYBAIxL2EDCBkFBkJXtSWaue+eOeSipJZb1169at+3m/\nXnnVPd/zPec8ladO1X3u+Z7v6Xpd9IOBQ0tXt5rffvB+kiRJkvqcmi6qkn4HLWC27Wpj47ZdbN/V\nzrYcSAQMHtjA8SOHEHGYYqpj6F+H6e+EN/5RN0UtSZIkqS+p6aJqY+ME2j/6UwB2tydbd5Rm8kvg\nz+54ipd3boMGoBH+1wdezZTxx1UvWEmSJEl9Uk0XVWta2/nAPz7IiCED2bazjS07du+zftTQgcyZ\neSLjhw/mVAsqSZIkSRVQ00UVwK62ZM3mHQA0DuhH44DS1OQN/YLLzpnM+dOPr2Z4kiRJkvq4mi+q\nOrv2PTM4fcLwaochSZIkqY70eFEVEZOAW4ATgHbgpsz8m4gYBXwfmAy8CPx+Zq4/1L5GD+7H/3/R\nabRnMnRgfwsqSZIkST2uGleqdgNXZeajETEMmB8RdwEfBe7OzOsj4vPA54HPHWpHgxrY98G9kiRJ\nktTDevxJtJm5IjMfLV5vBhYCE4A5wM1Ft5uBi3s6NkmSJEk6WlW9pyoiJgNnAr8BxmXmCigVXhHR\n5QwTEXEFcEWxuCMinu6BUNW7jAFeqXYQqgpzX5/Me30y7/XL3Nen3pr3k4+kU2RmpQPp+sARTcC9\nwHWZeUdEbMjMEZ3Wr8/MkYfZxyOZeValY1XvYt7rl7mvT+a9Ppn3+mXu61Ot573Hh/8BRMQA4Hbg\nO5l5R9G8KiLGF+vHA6urEZskSZIkHY0eL6oiIoBvAAsz86udVs0FLiteXwb8oKdjkyRJkqSjVY17\nqs4FPgI8FRGPF21/BlwP3BoRHweWAR88gn3dVJkQ1cuZ9/pl7uuTea9P5r1+mfv6VNN5r9o9VZIk\nSZLUF1TlnipJkiRJ6issqiRJkiSpDDVbVEXEhRHxXEQsjojPVzsedZ+ImBQRv4qIhRGxICI+U7SP\nioi7ImJR8XVk0R4R8bfFz8KTEfHa6n4HKkdENETEYxHxo2J5SkT8psj79yNiYNE+qFheXKyfXM24\ndewiYkRE3BYRzxbn/Rs93+tDRPxx8Xv+6Yj4bkQ0es73PRHxzYhY3fnZosdyjkfEZUX/RRFxWVfH\nUu9xkLz/VfG7/smIuDMiOj9O6QtF3p+LiHd0aq+J9/w1WVRFRAPw98A7gdOASyPitOpGpW60G7gq\nM08FzgY+VeT388DdmdkM3F0sQ+nnoLn4dwVwY8+HrG70GWBhp+UvAzcUeV8PfLxo/ziwPjNPAW4o\n+qk2/Q3ws8x8FfAaSvn3fO/jImIC8GngrMw8HWgALsFzvi/6NnDhfm1HdY5HxCjgWuANwOuBazsK\nMfVa3+bAvN8FnJ6ZrwZ+B3wBoHifdwkwo9jmH4oPWWvmPX9NFlWUTqbFmflCZu4EvgfMqXJM6iaZ\nuSIzHy1eb6b0BmsCpRzfXHS7Gbi4eD0HuCVLHgJGdDzzTLUlIiYC7wa+XiwH8FbgtqLL/nnv+Hm4\nDbig6K8aEhHHAW+m9KgNMnNnZm7A871e9AcGR0R/YAiwAs/5Picz7wPW7dd8tOf4O4C7MnNdZq6n\n9OZ8/zfs6kW6yntm/iIzdxeLDwETi9dzgO9l5o7MXAIspvR+v2be89dqUTUBeKnTckvRpj6mGN5x\nJvAbYFxmroBS4QUcX3Tz56Hv+GvgaqC9WB4NbOj0C7hzbvfkvVi/seiv2jIVWAN8qxj2+fWIGIrn\ne5+XmS8DX6H0GJUVlM7h+XjO14ujPcc99/uejwE/LV7XfN5rtajq6pMp54bvYyKiCbgd+GxmbjpU\n1y7a/HmoMRFxEbA6M+d3bu6iax7BOtWO/sBrgRsz80xgK3uHAXXFvPcRxdCtOcAU4ERgKKUhPvvz\nnK8vB8uz+e9DIuIaSrd7fKejqYtuNZX3Wi2qWoBJnZYnAsurFIsqICIGUCqovpOZdxTNqzqG+RRf\nVxft/jz0DecC742IFyld3n8rpStXI4qhQbBvbvfkvVg/nAOHl6j3awFaMvM3xfJtlIosz/e+723A\nksxck5m7gDuAc/CcrxdHe4577vcRxSQjFwEfzr0PzK35vNdqUfVboLmYIWggpRvb5lY5JnWTYoz8\nN4CFmfnVTqvmAh2z/VwG/KBT+38pZgw6G9jYMaRAtSMzv5CZEzNzMqVz+p7M/DDwK+ADRbf9897x\n8/CBon+v/PRKB5eZK4GXImJ60XQB8Aye7/VgGXB2RAwpfu935N5zvj4c7Tn+c+DtETGyuMr59qJN\nNSQiLgQ+B7w3M1s7rZoLXFLM8jmF0kQlD1ND7/mjVn8fRcS7KH2K3QB8MzOvq3JI6iYRcR7wa+Ap\n9t5b82eU7qu6FTiJ0h/jD2bmuuKP8dco3bDaClyemY/0eODqNhExG/hvmXlRREyldOVqFPAY8J8z\nc0dENAL/Qumeu3XAJZn5QrVi1rGLiJmUJicZCLwAXE7pQz/P9z4uIv4c+BClYUCPAX9A6X4Jz/k+\nJCK+C8wGxgCrKM3i9385ynM8Ij5G6f0AwHWZ+a2e/D50dA6S9y8Ag4C1RbeHMvOTRf9rKN1ntZvS\nrR8/Ldpr4j1/zRZVkiRJktQb1OrwP0mSJEnqFSyqJEmSJKkMFlWSJEmSVAaLKkmSJEkqg0WVJEmS\nJJXBokqS1KtFxDURsSAinoyIxyPiDT18/NkR8aOePKYkqbb0P3wXSZKqIyLeCFwEvLZ4VtEYSs+z\nkiSp1/BKlSSpNxsPvJKZOwAy85XMXB4RsyLi3oiYHxE/j4jxABFxSkT8MiKeiIhHI2JalPxVRDwd\nEU9FxIeKvrMjYl5E3BYRz0bEd4oHjxIRFxZt9wPv7wgmIt5SXC17PCIei4hhPf9fIknqbXz4rySp\n14qIJuB+YAjwS+D7wIPAvcCczFxTFEnvyMyPRcRvgOsz886IaKT04eE7gU8CFwJjgN8CbwCmAz8A\nZgDLgQeAPwUeARYBbwUWF8cckpkXRcQPi/0/UMS2PTN398T/hSSp9/JKlSSp18rMLcAs4ApgDaUC\n5xPA6cBdEfE48EVgYnHVaEJm3llsuz0zW4HzgO9mZltmrqJUkL2uOMTDmdmSme3A48Bk4FXAksxc\nlKVPHv+1U0gPAF+NiE8DIyyoJEngPVWSpF4uM9uAecC8iHgK+BSwIDPf2LlfRBx3kF3EIXa/o9Pr\nNvb+XexyGEdmXh8RPwbeBTwUEW/LzGcP/11Ikvoyr1RJknqtiJgeEc2dmmYCC4GxxSQWRMSAiJiR\nmZuAloi4uGgfFBFDgPuAD0VEQ0SMBd4MPHyIwz4LTImIacXypZ3imZaZT2XmlykNE3xVN32rkqQa\nZlElSerNmoCbI+KZiHgSOA3478AHgC9HxBOUhu2dU/T/CPDpou+DwAnAncCTwBPAPcDVmbnyYAfM\nzO2Uhhv+uJioYmmn1Z8tJrx4AtgG/LT7vlVJUq1yogpJkiRJKoNXqiRJkiSpDBZVkiRJklQGiypJ\nkiRJKoNFlSRJkiSVwaJKkiRJkspgUSVJkiRJZbCokiRJkqQyWFRJkiRJUhksqiRJkiSpDBZVkiRJ\nklQGiypJkiRJKoNFlSRJkiSVwaJKkiRJkspgUSVJkiRJZbCokiRJkqQy9K92AOUYMWJEnnLKKdUO\nQz1s69atDB06tNphqArMfX0y7/XJvNcvc1+femve58+f/0pmjj1cv5ouqsaNG8cjjzxS7TDUw+bN\nm8fs2bOrHYaqwNzXJ/Nen8x7/TL39am35j0ilh5JP4f/SZIkSVIZKlZURcSkiPhVRCyMiAUR8Zmi\n/X9ExMsR8Xjx712dtvlCRCyOiOci4h2Vik2SJEmSukslh//tBq7KzEcjYhgwPyLuKtbdkJlf6dw5\nIk4DLgFmACcCv4yI/5CZbRWMUZIkSZLKUrGiKjNXACuK15sjYiEw4RCbzAG+l5k7gCURsRh4PfDv\nlYpRkiRJUnl27dpFS0sL27dvP+Z9DB8+nIULF3ZjVEensbGRiRMnMmDAgGPaPjKzm0Pq4iARk4H7\ngNOBPwE+CmwCHqF0NWt9RHwNeCgz/7XY5hvATzPztv32dQVwBcDYsWNn3XrrrRWPX73Lli1baGpq\nqnYYqgJzX5/Me30y7/XL3NeepqYmxo0bx/Dhw4mIY9pHW1sbDQ0N3RzZkclMNm7cyKpVq9iyZcs+\n684///z5mXnW4fZR8dn/IqIJuB34bGZuiogbgS8BWXz938DHgK4ycEDFl5k3ATcBTJ8+PXvjLCGq\nrN46O4wqz9zXJ/Nen8x7/TL3tWfhwoVMnDjxmAsqgM2bNzNs2LBujOroDBs2jC1btnDWWYetn7pU\n0dn/ImIApYLqO5l5B0BmrsrMtsxsB/6Z0hA/gBZgUqfNJwLLKxmfJEmSpPKVU1D1BuXGX8nZ/wL4\nBrAwM7/aqX18p27vA54uXs8FLomIQRExBWgGHq5UfJIkSZLUHSp5pepc4CPAW/ebPv1/RcRTEfEk\ncD7wxwCZuQC4FXgG+BnwKWf+kyRJknQ4LS0tzJkzh+bmZqZOncqVV17Jjh07WLt2Leeffz5NTU1c\neeWVFTt+xYqqzLw/MyMzX52ZM4t/P8nMj2TmGUX7e4tZAju2uS4zp2Xm9Mz8aaVikyRJktQ3ZCbv\nf//7ufjii1m0aBGLFi1i27ZtXH311TQ2NvKlL32Jr3zlK4ffURkqPlGFJEmSpPrwnr+7/5i2O9zs\nfz/8/8476Lp77rmHxsZGLr/8cgAaGhq44YYbOPnkk7nuuus477zzWLx48THFdaQqOlGFJEmSJFXS\nggULmDVr1j5txx13HJMnT654MdXBK1WSJEmSusWhrigdSjlTqmdml7P39cTzeDt4pUqSJElSzZox\nYwaPPPLIPm2bNm1i1apVTJ8+vUdisKiSJEmSVLMuuOACWltbueWWW4DS/VlXXXUVV155JYMHD+6R\nGCyqJEmSJNWsiODOO+/ktttuo7m5mdGjR9OvXz+uueYaACZPnsyf/Mmf8O1vf5uJEyfyzDPPdHsM\n3lMlSZIkqaZNmjSJuXPnAvDggw9y6aWXMn/+fGbNmsWLL75Y8eNbVEmSJEnqM8455xyWLl3ao8d0\n+J8kSZIklcGiSpIkSZLKYFElSZIkSWWwqJIkSZKkMlhUSZIkSVIZLKokSZIk1bSWlhbmzJlDc3Mz\nU6dO5corr2THjh3cddddzJo1izPOOINZs2Zxzz33VOT4FlWSJEmSalZm8v73v5+LL76YRYsWsWjR\nIrZt28bVV1/NmDFj+OEPf8hTTz3FzTffzEc+8pGKxOBzqiRJkiR1j396yzFtNqS9Hfod4nrPJ+49\n6Kp77rmHxsZGLr/8cgAaGhq44YYbOPnkk7nuuutoamoCYMaMGWzfvp0dO3YwaNCgY4rzYLxSJUmS\nJKlmLViwgFmzZu3TdtxxxzF58mQWL168p+3222/nzDPP7PaCCrxSJUmSJKm7HOKK0qG0bt7MsGHD\njmnbzCQiumzvsGDBAj73uc/xi1/84piOcTheqZIkSZJUs2bMmMEjjzyyT9umTZtYtWoV06dPp6Wl\nhfe9733ccsstTJs2rSIxWFRJkiRJqlkXXHABra2t3HLLLQC0tbVx1VVX7ZkB8N3vfjd/+Zd/ybnn\nnluxGCyqJEmSJNWsiODOO+/ktttuo7m5mdGjR9OvXz+uueYavva1r7F48WK+9KUvMXPmTGbOnMnq\n1au7PYaK3VMVEZOAW4ATgHbgpsz8m4gYBXwfmAy8CPx+Zq6P0kDIvwHeBbQCH83MRysVnyRJkqS+\nYdKkScydOxeABx98kEsvvZT58+fzxS9+kS9+8YsVP34lr1TtBq7KzFOBs4FPRcRpwOeBuzOzGbi7\nWAZ4J9Bc/LsCuLGCsUmSJEnqg8455xyWLl16wIyAlVSxoiozV3RcacrMzcBCYAIwB7i56HYzcHHx\neg5wS5Y8BIyIiPGVik+SJEmSukN0nmqwYgeJmAzcB5wOLMvMEZ3Wrc/MkRHxI+D6zLy/aL8b+Fxm\nPrLfvq6gdCWLsWPHzrr11lsrHr96ly1btux5iJvqi7mvT+a9Ppn3+mXua8/w4cOZNm1al9OaH6m2\ntjYaGhq6Maqjk5k8//zzbNy4cZ/2888/f35mnnW47Sv+nKqIaAJuBz6bmZsO8Z/d1YoDKr7MvAm4\nCWD69Ok5e/bsbopUtWLevHmY9/pk7uuTea9P5r1+mfvas2TJEnbu3Mno0aOPubDaXMZzqsqVmaxd\nu5YRI0Zw5plnHtM+KlpURcQASgXVdzLzjqJ5VUSMz8wVxfC+juk3WoBJnTafCCyvZHySJEmSyjNx\n4kRaWlpYs2bNMe9j+/btNDY2dmNUR6exsZGJEyce8/aVnP0vgG8ACzPzq51WzQUuA64vvv6gU/uV\nEfE94A3AxsxcUan4JEmSJJVvwIABTJkypax9zJs375ivEvUGlbxSdS7wEeCpiHi8aPszSsXUrRHx\ncWAZ8MFi3U8oTae+mNKU6pdXMDZJkiRJ6hYVK6qKCScONqjygi76J/CpSsUjSZIkSZVQyedUSZIk\nSVKfZ1ElSZIkSWWwqJIkSZKkMlhUSZIkSVIZLKokSZIkqQwWVZIkSZJUBosqSZIkSSqDRZUkSZIk\nlcGiSpIkSZLKYFElSZIkSWWwqJIkSZKkMlhUSZIkSVIZLKokSZIkqQwWVZIkSZJUBosqSZIkSSqD\nRZUkSZIklcGiSpIkSZLKYFElSZIkSWWwqJIkSZKkMlhUSZIkSVIZLKokSZIkqQwVK6oi4psRsToi\nnu7U9j8i4uWIeLz4965O674QEYsj4rmIeEel4pIkSZKk7lTJK1XfBi7sov2GzJxZ/PsJQEScBlwC\nzCi2+YeIaKhgbJIkSZLULSpWVGXmfcC6I+w+B/heZu7IzCXAYuD1lYpNkiRJkrpL/yoc88qI+C/A\nI8BVmbkemAA81KlPS9F2gIi4ArgCYOzYscybN6+y0arX2bJli3mvU+a+Ppn3+mTe65e5r0+1nvee\nLqpuBL4EZPH1fwMfA6KLvtnVDjLzJuAmgOnTp+fs2bMrEqh6r3nz5mHe65O5r0/mvT6Z9/pl7utT\nree9R2f/y8xVmdmWme3AP7N3iF8LMKlT14nA8p6MTZIkSZKORY8WVRExvtPi+4COmQHnApdExKCI\nmAI0Aw/3ZGySJEmSdCwqNvwvIr4LzAbGREQLcC0wOyJmUhra9yLwCYDMXBARtwLPALuBT2VmW6Vi\nkyRJkqTuUrGiKjMv7aL5G4fofx1wXaXikSRJkqRK6NHhf5IkSZLU1xyyqIqIs3oqEEmSJEmqRYe7\nUvXPEbEoIv4iIk7rkYgkSZIkqYYcsqjKzDOBi4A24LaIeDwiPhcRJ/dIdJIkSZLUyx32nqrMfC4z\n/zwzTwMuA0YA90TEAxWPTpIkSZJ6uSOeqCIi+gHHA+OAocCaSgUlSZIkSbXisFOqR8SbgEuBiyk9\nrPd7wB9n5sYKxyZJkiRJvd4hi6qIeAlYRqmQ+vPMXNUjUUmSJElSjTjclarzMnNpRLwOuDEiJhfb\nBJCZ+eoKxydJkiRJvdohi6rMXFq8/FfgT4GngPZKByVJkiRJteKw91QV1mTm3IpGIkmSJEk16EiL\nqmsj4uvA3cCOjsbMvKMiUUmSJElSjTjSoupy4FXAAPYO/0vAokqSJElSXTvSouo1mXlGRSORJEmS\npBp0pA//fSgiTqtoJJIkSZJUg470StV5wGURsYTSPVVOqS5JkiRJHHlRdWFFo5AkSZKkGnVERVWn\n51VJkiRJkjo50nuqJEmSJEldsKiSJEmSpDJYVEmSJElSGSpWVEXENyNidUQ83altVETcFRGLiq8j\ni/aIiL+NiMUR8WREvLZScUmSJElSd6rklapvc+CsgZ8H7s7MZuDuYhngnUBz8e8K4MYKxiVJkiRJ\n3aZiRVVm3ges2695DnBz8fpm4OJO7bdkyUPAiIgYX6nYJEmSJKm7RGZWbucRk4EfZebpxfKGzBzR\naf36zBwZET8Crs/M+4v2u4HPZeYjXezzCkpXsxg7duysW2+9tWLxq3fasmULTU1N1Q5DVWDu65N5\nr0/mvX6Z+/rUW/N+/vnnz8/Msw7X70gf/ltp0UVbl9VeZt4E3AQwffr0nD17dgXDUm80b948zHt9\nMvf1ybzXJ/Nev8x9far1vPf07H+rOob1FV9XF+0twKRO/SYCy3s4NkmSJEk6aj19pWoucBlwffH1\nB53ar4yI7wFvADZm5ooejk3qFi9v2MaTL23grMmj6N8v+M2StZw2fjgnjR7CsrWtPPnyBgBmnDic\nTdt28dL6VgDa2pP1W3cyZtigPfsaOrA//foFm7fvOuxxW3e2sbstOW7wvqd1xz5eWLOF0U2D6Fdc\nF26I4HVTRjGmaVAXe5MkSdKRqlhRFRHfBWYDYyKiBbiWUjF1a0R8HFgGfLDo/hPgXcBioBW4vFJx\nqfZs3bF7T+EB8NLmdp5duemIt39m+SZ+vmAlyzdsZ8KIwby8YRsAk8cM7fZYAV58ZWvx6vl92ieP\nGdppXS8x7/mK/T8ciZ2723jvayYw7fgji+FIc3/SqCEMGdhbRjdLkqS+rmLvOjLz0oOsuqCLvgl8\nqlKxqPLa27Prm+A6efyl9dy/aO1R7/uXC1fts7x58w5uffHJo94PsKegAnq8wDnU8d51xnh+8tSK\ngy53tB3MK1t28PCS0mSbJ40awukThgPwswUraW/fNzPvOmM8qzZtZ/7S9YeNqyf8473PH75T4Why\n/7ZTxx1rSEdsyMAGPjBrIscNHlDxY+0vgH79urodVZIk9TQ/ytVhbd/VxoPPv8KLr7Ry/HEHDhVb\ns3kHdzz6csXjaD6+iYhgZW7lhHHDjmrb363azOimgexqa2fowP6MGDKAK948jUq9Jx3Yvx87d7cf\n8BrghOGNBMGKjduICCaNHEz/hn5cfu5k1mzewcSRg4kIPvmWqSzfuJ22tmTSqFLboWxo3cmutmRs\np+GDfzh7Gi9v2EZbW6mwGjd8EIP6N+zpv27rzu7+1o/Yyo3b+f4jL9G/35Hf2nm43Gcmi1ZvAQ4s\nxitl7hPVu/1zzswTOWF44zFtu3rTDo4/bhCnjj+OaWN732xLkiTVEouqOrSxdRcrNm07oP2uBat4\nduXmAz79PporGYcrUtqz9EZ/QMPRzZHSfHzTnmFq8+atZ/bs1xzV9r3R1P3eyDYOaGDSqCF7liOC\nCSMGH/H+RgwZ2GX7wfYxYsjAg27TE6aObeKcU8Yc1TZHkvslr2xlcVFYVdLDS9by0AvrKlaYH0rH\nxccfPN49BV13DAHNTI4f1siaLTsAGNAv+K9vnsphPgs4Ikcy7HPy6KE0Dmgo/2CSJB0Di6o+Zv9h\neGu37OCOx15mx67SlZIdu9v49aJXyjpGV8PQ+vcL3jHjBE4aPaSLLaSeM2XMUKb0wH1iv3da5YcX\nHsxL61r52dMr2d1+bM8Z3H9oaXcNAV26tnWf5atvO7Zhuvs70mGflR7yOX54I+9/7YTDXjUuV7+g\n4seQJHUvi6oatXn7Lh58fi272vYOK/v359fyZMvGI97Hf+hiGFVDP7ocFjegoR/tmUwaOcT7OKQq\nmzRqCP/1zVOPefuPnzeFlRu3M3bYIFZsPPCq9bH4l4eWsmnbbjKTF17ZytSxQ4kuH0F49A417HP7\nrjaWrSsVcz0x5PNfHlpa8WMA/MGbptDQC37XDh7QwLmnjPEqoCQdhkVVL9WyvpWtO9r2LN/87y+y\nefvuPcuH+2S589/i9oS3vur4PZMXAJwxYfgx34shqbYN7N9vz1Xl/YehHqtr3zOjW/bTlcMN+1y0\najMv7neVrLt95zdLWbtlZ8WHe3ZcfPz6r5dU9kBH4a9/uagqs4SuXrWd21c8VrH9d/47OmHEYM45\nZTRvmDK6Ysc7Wk2N/Y9qCLik6rKoqqKOoXqPLlvPg4v3zor37y+8sk9BdSidZ3oDaBrUwMVnTmBY\nY8/PRiZJ1dA8bhjNRzl5zdHqqeGej7+0gX9//uhnSa2EBxa/wsZtpWfkVWOW0M2tSWsPHfflDdv4\nt0da+LdHWnrkeEfq2vecxlmTR1U7DElHwKKqB23bWZpFb/uudl5cu5WfPb3ysNs0jys+RU4YP6KR\n//jaiXvWDWscsM9Mb5Kk2jZz0ghmThpR7TCA0qRCKzdup3Xn7sN3roCHH36Y179+ZkX2vXLjdv7P\nw8tYuraVE0c0snzD9r1/b3uBNZt3sKF1Fys3ba92KJKOkEVVBby0rpXWnW1s2r6Lr//6BQYWU1gf\n7JO+flEa8vFHs6fRv5gVb1D/frx+yijHsUuSqqaaw8SXDe3XbcNT93css4/2pH+893l+/OSKw3eU\n1GtYVJUhM/eMf7/vd2t4smUj9y9ew/Zd7YfcbuTQgZw9dRQDG/rx7lePZ/xwx0xLkiRJtcqi6ii0\ntycPLVnL+q272NXWzjfuP/SNxM3jmtjdlsyePnbPcI6hg/oz7jgniJAkSZL6CouqQ1i1aTsbWnfx\n8wUreW7VZpYdZHapjtmg2hOufOspDB7QwBumjmJQf4fuSZIkSX2dRVWhYya+x5at54HFa/nd6oMX\nUQDvPOMEguBNp87OUwAAEoxJREFUzWP2mX1PkiSpO7Qd4wO+JfW8ui6qXlrXypMtGw87E1/zuCYa\nIvjEW6YysKGBiSMH+wBcSZJUUV//9RJ+9vRKLnr1iQAMGdjAG6eNdhIrqReqq6KqYzjfLxeu4pnl\nm1i27sArUR0z8f3h7GkMHtDAWZNH+swnSZLUY149cfie2f9a1m/jH+99fu/Ku2DKmKHsamvnTc1j\nmXXySCJg8uihDOzfr0oRS+rTRVVmkgkrNm3nr+/6Hc+u3Nxlv7edOo6hgxqciU+SJFXdOdPGcMcf\nncPt81tY17oTgHufW0PrzjYAlhSPaPnuw8v47sPL9mz3tlNLD6nu3xC8+4zxnDRqCICja6Qe0CeL\nqsWrt/DMik38830vdLm+eVwT/SL4xJunMn7EYJoG9cn/BkmSVKMGNPTjkteftGf5j2afwoqN29i2\ns43N23fzrQeW0K9fsLF1F6s37wDglwtX7enf+baGoYMa+MjZk4v9Bm+cNtpROFI36xPVxO62dpa8\nspUfPL6cpetaD3jIbseQvjkzT+RDr5vkLxJJklRzOo+m+etLztzzesHyjby8fhsAyzds4/ZHX97z\n3gdg6462fYYQ/t09i5kyZiht7cnMSSN4y/SxAJw0aoj3a0nHqOaLql8sWMnf3bO4y3XvPOMEzp46\nmteeNLKHo5IkSeoZM04czowT985E/NFzpwCl2yDmPrGclzeUCq4Fyzftmdm4YwjhsnWtzH1i+Z5t\nO4YQAgxr7M/vv24SQ4pCKwIiHEoodaWmi6qVW9v3KaiGDmpg6tgmPn7eFCaMGOynLZIkqW5FBHNm\nTtinbfXm7WzZvptdbck/FVevdra1s7QotjoPIQS487GX91n+2HmT9zyHc0BDP944bbS3UUjUeFHV\ncVm7X8BfzDmd10waUd2AJEmSerHjhzVy/LDS669+aOae9sWrN/PCmr23T9z7uzU82bKRjjkuOt5z\nffP+F/fZ39/evYgpY4buWR7W2J/LzpkMwJimQYwaOrDbvwepN6rpogrgrMkjufY9M6odhiRJUs06\n5fhhnNJRbQFvn3HCPuuffnkj9y1aQ+be5ZbiPq4l+93LftWtT+x5PevkkXsKq/4NwUVnnMjEkfvO\ntOzshOoLqlJURcSLwGagDdidmWdFxCjg+8Bk4EXg9zNzfTXikyRJ0l6nTxjO6ROG79PWMZSww23z\nW1i5cTsAi1ZvAWD+0n3fyv30qZXs76RRQ3j3q8fvWV6zqa3b4pZ6SjWvVJ2fma90Wv48cHdmXh8R\nny+WP1ed0CRJknQonYcSAlx94av2vH5lyw4eW7aB9uLS1kvrWvnB48vpfFGqY0jhsnWt3Dhv7+yE\nmzfv5L51v2Xofvdq/eHsaTQUO+gXASQnjfKhx+odetPwvznA7OL1zcA8LKokSZJqzpimQfzeaeP2\nafuDN03dZ3nn7nZum9/C+uIBx7D3+VqrN++A4vlbHa6+7ckuj9VxnIZ++z70uDOHGKrSIjsGx/bk\nQSOWAOuBBP4pM2+KiA2ZOaJTn/WZecBc6BFxBXAFwOBxU2Z9+Nob+fCpg3oqdPUCW7Zsoampqdph\nqArMfX0y7/XJvNen3e3J0rVbGTpkb2H06OrdvLS5naBUGC3f0n7U+z17fH9GD95bWG3amRw3MDht\ndANDB1hw9Qa99Zw///zz52fmWYfrV60rVedm5vKIOB64KyKePdINM/Mm4CaAYROac+LEicye7UQV\n9WTevHnMnj272mGoCsx9fTLv9cm816/9c3/Jfut3t7Xz8JJ1vLi2ldFNpUkwVnR66HFnHUMMF2wB\nthx4rPvWsGf2wh2721i+YTsTRw7mo+dMZuR+MxeOH97IsMYBZXxnOpRaP+erUlRl5vLi6+qIuBN4\nPbAqIsZn5oqIGA+srkZskiRJ6r36N/TjnFPGcM4p+7Z3PPS4s5fWtfLTp1ewq23vyKyOIYYd9p+9\nsGX9Nv7njxd2eezzX3U8/TtVbis2bmPH7nbe85oTeUvz2AP6+8Dk+tHjRVVEDAX6Zebm4vXbgb8A\n5gKXAdcXX3/Q07FJkiSp75g0aghXvHnaPm1/NHvanungd7btHUr4oydW7Hn4cfPx+w5D65jN8FfP\ndv2Z/1d/8Tu++ovfdbnuP752Ascf17hP24CGfpw9dZRXvvqQalypGgfcWVTt/YH/k5k/i4jfArdG\nxMeBZcAHqxCbJEmS+rCIYFIXk1l85m3NfOZtzV1u07K+lQXLN+3TtrF1F//y0NJO+913m45pC25/\n9OUu9/m3d7PPg5M76xfwydnTaDjEVa4JIwczZGBvmnOuvvV4JjLzBeA1XbSvBS7o6XgkSZKkQ5k4\ncggTRx5YiP3+6yYddJtla0tDD3e37zsp3DPLN7FsXStw4NDDzv7037qe7bCzt5067oBi7tFl65k5\naQQfPGsS4/e7QtbBYYndz/JWkiRJ6mYnjR7CJ94yrct1qzdvZ3OnByd3dtczq1i4YtMhr1J1DEfs\nGK64v7sXrubuhYeenuCj50xmyMCGg66fOraJ6ScMO+h67cuiSpIkSepB+z84ubNpbzn8tOJL127l\n2ZWbD2j/5TOr9rQfrCbrGJb47QdfPOxxJowYfNCHKy/fUJqkY8qYofzh7Gn7TOBxMJNGDaFxwMEL\nuVpmUSVJkiTVkJNHD+Xk0Qfej/WOGSccdtuFKzbxq+dWc6hH1XbMkPjyhm2H3d+SV7Ye9MHMXXn7\nfg+F7rD8pV3MPuK99D4WVZIkSVKdOHX8cZw6/rhD9vmDN03ZM0PiwXzrgSU88dJG4MDZEve3qz15\nsbh/7BfPdD1kMXe0HXIfvZ1FlSRJkqQ9BvVvYNrYQxdK//PiM45qn4tXb2bx6oNPzPG7hU8f1f56\nG4sqSZIkSRV1yvHDOOVgN5IBja8824PRdL+u7zyTJEmSJB0RiypJkiRJKoNFlSRJkiSVwaJKkiRJ\nkspgUSVJkiRJZbCokiRJkqQyWFRJkiRJUhksqiRJkiSpDBZVkiRJklQGiypJkiRJKoNFlSRJkiSV\nwaJKkiRJkspgUSVJkiRJZbCokiRJkqQyWFRJkiRJUhl6XVEVERdGxHMRsTgiPl/teCRJkiTpUPpX\nO4DOIqIB+Hvg94AW4LcRMTczn+mq/0B2MWbbi7CmV30bqrAhW1+CNc9VOwxVgbmvT+a9Ppn3+mXu\n69Pg1uXVDqEsva0aeT2wODNfAIiI7wFzgC6LqvG5mve9+BdwR2MPhqhqO3XzZmgZVu0wVAXmvj6Z\n9/pk3uuXua9PzdsD+E/VDuOY9baiagLwUqflFuANnTtExBXAFQDNJzSxdeDxrGqLnotQVbdr4FBa\n23rbj656grmvT+a9Ppn3+mXu69O2/gN5ct68aodxzHrbT2xX1VHus5B5E3ATwPTp07P5sz/sibjU\ni8ybN4/Zs2dXOwxVgbmvT+a9Ppn3+mXu61Ot5723TVTRAkzqtDwRqO0BlpIkSZL6tMjMw/fqIRHR\nH/gdcAHwMvBb4D9l5oKD9N8MeCdj/RkDvFLtIFQV5r4+mff6ZN7rl7mvT7017ydn5tjDdepVw/8y\nc3dEXAn8HGgAvnmwgqrwXGae1TPRqbeIiEfMe30y9/XJvNcn816/zH19qvW896qiCiAzfwL8pNpx\nSJIkSdKR6G33VEmSJElSTan1ouqmagegqjDv9cvc1yfzXp/Me/0y9/WppvPeqyaqkCRJkqRaU+tX\nqiRJkiSpqiyqJEmSJKkMNVtURcSFEfFcRCyOiM9XOx51n4iYFBG/ioiFEbEgIj5TtI+KiLsiYlHx\ndWTRHhHxt8XPwpMR8drqfgcqR0Q0RMRjEfGjYnlKRPymyPv3I2Jg0T6oWF5crJ9czbh17CJiRETc\nFhHPFuf9Gz3f60NE/HHxe/7piPhuRDR6zvc9EfHNiFgdEU93ajvqczwiLiv6L4qIy6rxvejIHSTv\nf1X8rn8yIu6MiBGd1n2hyPtzEfGOTu018Z6/JouqiGgA/h54J3AacGlEnFbdqNSNdgNXZeapwNnA\np4r8fh64OzObgbuLZSj9HDQX/64Abuz5kNWNPgMs7LT8ZeCGIu/rgY8X7R8H1mfmKcANRT/Vpr8B\nfpaZrwJeQyn/nu99XERMAD4NnJWZp1N6PuUleM73Rd8GLtyv7ajO8YgYBVwLvAF4PXBtRyGmXuvb\nHJj3u4DTM/PVwO+ALwAU7/MuAWYU2/xD8SFrzbznr8miitLJtDgzX8jMncD3gDlVjkndJDNXZOaj\nxevNlN5gTaCU45uLbjcDFxev5wC3ZMlDwIiIGN/DYasbRMRE4N3A14vlAN4K3FZ02T/vHT8PtwEX\nFP1VQyLiOODNwDcAMnNnZm7A871e9AcGR0R/YAiwAs/5Picz7wPW7dd8tOf4O4C7MnNdZq6n9OZ8\n/zfs6kW6yntm/iIzdxeLDwETi9dzgO9l5o7MXAIspvR+v2be89dqUTUBeKnTckvRpj6mGN5xJvAb\nYFxmroBS4QUcX3Tz56Hv+GvgaqC9WB4NbOj0C7hzbvfkvVi/seiv2jIVWAN8qxj2+fWIGIrne5+X\nmS8DXwGWUSqmNgLz8ZyvF0d7jnvu9z0fA35avK75vNdqUdXVJ1PODd/HREQTcDvw2czcdKiuXbT5\n81BjIuIiYHVmzu/c3EXXPIJ1qh39gdcCN2bmmcBW9g4D6op57yOKoVtzgCnAicBQSkN89uc5X18O\nlmfz34dExDWUbvf4TkdTF91qKu+1WlS1AJM6LU8EllcpFlVARAygVFB9JzPvKJpXdQzzKb6uLtr9\neegbzgXeGxEvUrq8/1ZKV65GFEODYN/c7sl7sX44Bw4vUe/XArRk5m+K5dsoFVme733f24Almbkm\nM3cBdwDn4DlfL472HPfc7yOKSUYuAj6cex+YW/N5r9Wi6rdAczFD0EBKN7bNrXJM6ibFGPlvAAsz\n86udVs0FOmb7uQz4Qaf2/1LMGHQ2sLFjSIFqR2Z+ITMnZuZkSuf0PZn5YeBXwAeKbvvnvePn4QNF\n/1756ZUOLjNXAi9FxPSi6QLgGTzf68Ey4OyIGFL83u/Ived8fTjac/znwNsjYmRxlfPtRZtqSERc\nCHwOeG9mtnZaNRe4pJjlcwqliUoepobe80et/j6KiHdR+hS7AfhmZl5X5ZDUTSLiPODXwFPsvbfm\nzyjdV3UrcBKlP8YfzMx1xR/jr1G6YbUVuDwzH+nxwNVtImI28N8y86KImErpytUo4DHgP2fmjoho\nBP6F0j1364BLMvOFasWsYxcRMylNTjIQeAG4nNKHfp7vfVxE/DnwIUrDgB4D/oDS/RKe831IRHwX\nmA2MAVZRmsXv/3KU53hEfIzS+wGA6zLzWz35fejoHCTvXwAGAWuLbg9l5ieL/tdQus9qN6VbP35a\ntNfEe/6aLaokSZIkqTeo1eF/kiRJktQrWFRJkiRJUhksqiRJkiSpDBZVkiRJklQGiypJkiRJKoNF\nlSSpV4uIayJiQUQ8GRGPR8Qbevj4syPiRz15TElSbel/+C6SJFVHRLwRuAh4bfGsojGUnmclSVKv\n4ZUqSVJvNh54JTN3AGTmK5m5PCJmRcS9ETE/In4eEeMBIuKUiPhlRDwREY9GxLQo+auIeDoinoqI\nDxV9Z0fEvIi4LSKejYjvFA8eJSIuLNruB97fEUxEvKW4WvZ4RDwWEcN6/r9EktTb+PBfSVKvFRFN\nwP3AEOCXwPeBB4F7gTmZuaYokt6RmR+LiN8A12fmnRHRSOnDw3cCnwQuBMYAvwXeAEwHfgDMAJYD\nDwB/CjwCLALeCiwujjkkMy+KiB8W+3+giG17Zu7uif8LSVLv5ZUqSVKvlZlbgFnAFcAaSgXOJ4DT\ngbsi4nHgi8DE4qrRhMy8s9h2e2a2AucB383MtsxcRakge11xiIczsyUz24HHgcnAq4AlmbkoS588\n/munkB4AvhoRnwZGWFBJksB7qiRJvVxmtgHzgHkR8RTwKWBBZr6xc7+IOO4gu4hD7H5Hp9dt7P27\n2OUwjsy8PiJ+DLwLeCgi3paZzx7+u5Ak9WVeqZIk9VoRMT0imjs1zQQWAmOLSSyIiAERMSMzNwEt\nEXFx0T4oIoYA9wEfioiGiBgLvBl4+BCHfRaYEhHTiuVLO8UzLTOfyswvUxom+Kpu+lYlSTXMokqS\n1Js1ATdHxDMR8SRwGvDfgQ8AX46IJygN2zun6P8R4NNF3weBE4A7gSeBJ4B7gKszc+XBDpiZ2ykN\nN/xxMVHF0k6rP1tMePEEsA34afd9q5KkWuVEFZIkSZJUBq9USZIkSVIZLKokSZIkqQwWVZIkSZJU\nBosqSZIkSSqDRZUkSZIklcGiSpIkSZLKYFElSZIkSWX4f9XNMRfi23TgAAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# device initialization\n", "with TCLab() as a:\n", " a.initplot(tf)\n", " for t in clock(tf,h):\n", " r = Tset(t) - Tref\n", " y = a.T1 - Tref\n", " \n", " P = kp*(b*r - y)\n", " v = P + I\n", " \n", " u = max(0,min(200,v + uref))\n", " I += bi*(r-y)\n", " \n", " a.Q1 = u\n", " a.updateplot()\n", " " ] }, { "cell_type": "code", "execution_count": null, "metadata": { "collapsed": true, "nbpages": { "level": 2, "link": "[4.3.2 PID Implementation](https://jckantor.github.io/CBE32338/04.03-PID_Control.html#4.3.2-PID-Implementation)", "section": "4.3.2 PID Implementation" } }, "outputs": [], "source": [] }, { "cell_type": "markdown", "metadata": { "nbpages": { "level": 2, "link": "[4.3.2 PID Implementation](https://jckantor.github.io/CBE32338/04.03-PID_Control.html#4.3.2-PID-Implementation)", "section": "4.3.2 PID Implementation" } }, "source": [ "\n", "< [4.1 Relay Control](https://jckantor.github.io/CBE32338/04.01-Relay-Control.html) | [Contents](toc.html) | [4.10 Lab Assignment: PID Control](https://jckantor.github.io/CBE32338/04.10-Lab-Assignment-PID-Control.html) >

\"Open

\"Download\"" ] } ], "metadata": { "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.6.8" } }, "nbformat": 4, "nbformat_minor": 2 }