{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# 301 古典制御\n", "## 読み込み" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "from pylab import *\n", "import control as matlab\n", "import matplotlib.pyplot as plt\n", "%matplotlib inline\n", "import pandas as pd\n", "import warnings\n", "warnings.filterwarnings(\"ignore\")" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## パラメータ定義" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "TIME_MAX = 10 # min\n", "\n", "V = 10\n", "n = 10\n", "K = 1000\n", "K_1 = 100\n", "a = 100\n", "R_a = 8\n", "J_a = 0.02\n", "D_a = 0.01\n", "K_b = 0.5\n", "K_t = 0.5\n", "N_1 = 25\n", "N_2 = 250\n", "N_3 = 250\n", "J_K = 1\n", "D_L = 1\n", "\n", "K_pot = 0.318\n", "K_m = 2.083\n", "a_m = 1.71\n", "K_g = 0.1" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## 伝達関数定義" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "\n", " 66.24\n", "--------------------\n", "s^2 + 1.71 s + 66.24\n", "\n" ] } ], "source": [ "TF_Potentiometer = matlab.tf([K_pot],[1])\n", "TF_Preamplifer = matlab.tf([K],[1])\n", "#TF_PowerAmplifier = matlab.tf([K_1], [1, a])\n", "TF_PowerAmplifier = matlab.tf([1], [1])\n", "TF_Moter = matlab.tf([K_m], [1, a_m, 0])\n", "TF_Gear = matlab.tf([K_g], [1])\n", "\n", "sysPI = TF_Potentiometer * matlab.feedback(TF_Preamplifer*TF_PowerAmplifier*TF_Moter*TF_Gear, TF_Potentiometer)#If u can not make a closed loop tf, use it\n", "print(sysPI)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## 応答のグラフ化" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY4AAAEWCAYAAABxMXBSAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvNQv5yAAAIABJREFUeJzt3Xm8lHXd//HXGzgssitHQUADQRRX1HArO5UWlmm7kC1ayV231n23qmXqT6uHd3bncmcpGi65ZWrmvqSeNBUVU0ElBQ4qR3ZZDwc42+f3x/eaGA4zZ66ZMxdzZs7n+XjM45q5lu98r+Ewn/nuMjOcc865uHqUOgPOOefKiwcO55xzefHA4ZxzLi8eOJxzzuXFA4dzzrm8eOBwzjmXFw8crixJul7Sz6PnNZLqk0h7R16bI91TJD1S7HSLrRj/FuVyr92ZBw6XlaQvSZotqUHSUkkPSvpAEdJN5Mu1kpnZzWb2sVLnY0foTvdarjxwuIwkfR+4DPglsBuwB/A74KQd8N69kn6PUrxXKSjw/+euqPwPym1H0mDgQuAMM7vLzDaaWbOZ3WtmP4rO6SPpMklLosdlkvpEx2ok1Uv6gaQVUWnltOjYdOAU4MdRSebeaP9bks6SNAfYKKmXpH0l1UpaK+k1SSfGzP/uku6UtFLSIknfTTt2gaQ7JN0kaT1wapZkhkl6VNIGSX+XtGdaGvtEx1ZLekPSF9tdO1TS/dG1z0naK+3ayyUtlrRe0ouSPpiW502Sdk47d5KkVZKqJJ0q6R9px46S9IKkddH2qLRjtZJ+IelpoBEYm+EzOlvSwiiPr0v6TNqxUyX9Q9KvJa2JPsPj046fJmledG2dpP/I8u/wI0l3ttv3f5IuS3ufuiidRZJOSX//6LkkXRr9Ha2TNEfS/ln+zdyOYmb+8Mc2D2AK0AL06uCcC4FZwK5ANfAMcFF0rCa6/kKgCvgE4QtsaHT8euDn7dJ7C3gZGA30i65bAPwE6A18BNgATGifRvR+9dHzHsCLwHnRdWOBOuDj0fELgGbg09G5/TLc2/XRex0D9AEuB/4RHesPLAZOA3oBhwCrgP3Srl0NTI6O3wzclpb2l4FdomM/AJYBfaNjjwOnp517CXBV9PzUtDzsDKwBvhKlMy16vUt0vBZ4B9gvOl6V4R6/AOwefQYnAxuBEWnv1QycDvQEvg0sARQd/ySwFyDgQ9G/7SEZ/i1GROkOiV73AlYAh0af4/q0f88RaZ9h+r1+PPr3HBK9376pfPqjdA8vcbhMdgFWmVlLB+ecAlxoZivMbCXw/whfZCnN0fFmM3sAaAAm5HjfK8xssZltAo4ABgAXm1mTmT0O3Ef4kuzI+4FqM7swuq4OuAaYmnbOs2Z2t5m1Re+Vyf1m9qSZbQF+ChwpaTRwAvCWmV1nZi1m9k/gTuDzadfeZWbPR5/fzcDBqQNmdpOZvRdd+7+EwJT6XG5J3Z8kRXm+JUPePgnMN7M/RuncCvwL+FTaOdeb2WvR8eb2CZjZn81sSfQZ/AmYTwh2KW+b2TVm1grcQPhi3y269n4zW2jB34FHgA9meI+lwJOEIAXhB8kqM3sxet0G7C+pn5ktNbPXMtxrMzAQ2IcQuOZF6boS8sDhMnmPUFXTUf3/7sDbaa/fjvb9O412gaeREAg6srhd+ovNrK3de4zMkcaewO5R9dZaSWsJpZbdsrxPzryYWQOhFLF7lP7h7dI/BRiedu2ytOfb3HdUfTcvqnZZCwwGhkWH7yAEqN0JpR0DnsqQt/afPWz/2XR4j5K+KunltHvYPy0f29yDmTVGTwdE1x4vaVZUVbeWUKJMvzbdDYRSFtH2j1GaGwklnW8BS6OqvX3aXxz9YPgtcCWwXNIMSYM6ujeXPA8cLpNngc2E6pxslhC+RFP2iPbFkW1K5vT9S4DR2rZhdw/g3RxpLwYWmdmQtMdAM/tEjPdPNzr1RNIAQvXQkij9v7dLf4CZfTtXglF7xlnAFwnVdkOAdYQqGMxsLeHX+xeBLwG3mlmmvLb/7GH7zybrPUbtNdcAZxKqt4YAr6bykeMe+hBKWL8GdouufaCDa+8GDozaJU4glMBCBs0eNrPjCKWZf0V52o6ZXWFmhxKq3vYGfpQrny5ZHjjcdsxsHaGN4EpJn5a0U9RAe7ykX0Wn3QqcK6la0rDo/JtivsVyMjTYtvMcoX78x9F71xCqYm7Lcd3zwHqFhvZ+knpK2l/S+2PmLeUTkj4gqTdwEfCcmS0mVJftLekrUb6qJL1f0r4x0hxIaPtZCfSSdB7Q/tfzLcBXgc+RuZoKwhf13grdpXtJOhmYGOUtjv6EwLISQmM3ocQRR29C9dpKoCVqNM/addbMNhNKUrcAz5vZO9F77ibpREn9gS2EqszW9tdHn+3hkqoIfw+bM53ndiwPHC4jM/sN8H3gXMKXxGLCL9S7o1N+DswG5gBzgX9G++L4AzAxqia5O9MJZtYEnAgcT2h8/h3wVTP7V458txICzMHAoujaawlVQvm4BTifUEV1KKE6CjPbQPiinEr45b8M+B/Cl2kuDwMPAm8SqpY2s32V0j3AeGC5mb2SKREze4/w6/0HhGrFHwMnmNmqODdmZq8D/0soWS4HDgCejnntBuC7wO2EBvkvRXnuyA3Re/wxbV+PKP9LCJ/xh4D/zHDtIEJJZA3hM3uPUNpxJaTMJWHnnCsOSXsQqqKGm9n6UufHdZ6XOJxziYnaqL5P6JLsQaNCVPSoWedc6UTtF8sJVUxTSpwdV0ReVeWccy4vXlXlnHMuLxVVVTVkyBAbN25cqbORiI0bN9K/f/9SZyMxfn/lze+vfL344ourzKw6n2sqKnDstttuzJ49u9TZSERtbS01NTWlzkZi/P7Km99f+ZLUfhaCnLyqyjnnXF4SK3FImkkYpLTCzLYblSrpR0SDqqJ87EuYnG61pLcIs5O2Ai1mdlhS+XTOOZefJEsc19NBFzwzu8TMDjazg4FzCPP/rE475cPRcQ8azjnXhSQWOMzsScJUAnFMI8x95JxzrotLdByHpPcB92Wqqko7ZyegHhiXKnFIWkSYm8aAq81sRgfXTwemA1RXVx96++23Fy3/XUlDQwMDBuSalbx8+f2VN7+/8vXhD3/4xXxrdrpCr6pPAU+3q6Y62syWSNoVeFTSv6ISzHaioDIDYMKECVapPR8quVcH+P2VO7+/7qUr9KqaSrtqKjNbEm1XAH9h25XJnHPOlVBJA4ekwYTplP+atq+/pIGp54QprF8tTQ7LUGsr/OEP0NRU6pw45ypUkt1xbyUsXD9MUj1hbYMqADO7KjrtM8Aj0TKSKbsBfwlLLtMLuMXMHkoqnxXnqafgm9+E6mo48cRS58Y5V4ESCxxmNi3GOdcTuu2m76sDDkomV93AggVhuyTuKq7OOZefrtDG4Ypp0aKwXbastPlwzlUsDxyVpq4ubD1wOOcS4oGj0njgcM4lzANHpfHA4ZxLmAeOSrJ+PaxaFZ574HDOJcQDRyVJNYyPGRMChy8L7JxLgAeOSpKqpjrqKNiyBdatK21+nHMVyQNHJUkPHODVVc65RHjgqCR1dTBkCOy7b3jtgcM5lwAPHJWkrg7GjoXhw8NrDxzOuQR44KgkHjicczuAB45K0doKb70VAseQIdC7twcO51wiPHBUiiVLwlTqY8eCFEodHjiccwnwwFEpUj2qxo4N2+HDYenS0uXHOVexPHBUikyBw0sczrkEeOCoFIsWQY8esMce4bUHDudcQjxwVIq6uhA0qqrC6+HDYeVKaGkpbb6ccxXHA0elSHXFTRk+PMxVtXJl6fLknKtIHjgqRabAAV5d5ZwrOg8claClBZYvh1Gjtu7zwOGcS0higUPSTEkrJL2a5XiNpHWSXo4e56UdmyLpDUkLJJ2dVB4rxsaNYTtw4NZ9HjiccwlJssRxPTAlxzlPmdnB0eNCAEk9gSuB44GJwDRJExPMZ/lrbAzb/v237vPA4ZxLSGKBw8yeBFYXcOlkYIGZ1ZlZE3AbcFJRM1dpUiWOnXbauq9fPxg82AOHc67oepX4/Y+U9AqwBPihmb0GjAQWp51TDxyeLQFJ04HpANXV1dTW1iaX2xJqaGjIem/9Fyzg/cCrixaxKu2cyYMG0TBnDq+XwWfS0f1VAr+/8lbp95evUgaOfwJ7mlmDpE8AdwPjAWU4N+saqGY2A5gBMGHCBKupqUkgq6VXW1tL1nvr3RuA/Q8/HNLPGTuWndra2LUMPpMO768C+P2Vt0q/v3yVrFeVma03s4bo+QNAlaRhhBLG6LRTRxFKJC6bTG0c4KPHnXOJKFngkDRckqLnk6O8vAe8AIyXNEZSb2AqcE+p8lkWMrVxgE906JxLRGJVVZJuBWqAYZLqgfOBKgAzuwr4PPBtSS3AJmCqmRnQIulM4GGgJzAzavtw2aQCR6YSx4YN4Xj7Y845V6DEAoeZTctx/LfAb7McewB4IIl8VaRsVVU77xy2a9d64HDOFY2PHK8E2aqqUgMCN2zYsflxzlU0DxyVIFtVlQcO51wCPHBUgsZG6Nnz391y/23AgLD1wOGcKyIPHJVg48ZQTaV2Q2C8xOGcS4AHjkqQrdeUBw7nXAI8cFSCxsaOA0dDw47Nj3OuonngqASpqqr2vMThnEuAB45KkK2qqn//0O7hgcM5V0QeOCpBtqoqKfSs8sDhnCsiDxyVIFtVFYTqqmIFDss6SbFzrhvxwFEJOpqLqliBo7ERDj0ULr6482k558qaB45KkK2qCooXOM49F156CXwxG+e6PQ8claCjqqpitHE88wxcdlkYnb5wYefScs6VPQ8clSDJqqpNm+C002CPPeCMM+Ctt6C5ufD0nHNlr9RrjrvOam4Oj44CR2cGAF58Mbz5Jvztb/DOO9DSErZ77VV4ms65suYljnKXWosjqV5VTz0FRxwBH/0ojBsX9i1YUHh6zrmy54Gj3GWbUj2ls4Fj2TIYOTI8T5UyvJ3DuW7NA0e5y7b6X8rAgaGdoqWlsPSXLQtL0AKMGAH9+nmJw7luzgNHuYtT4oDC2jm2bIE1a7YGDilUV3ngcK5b88BR7rItG5vSmYkOV6wI21TggFBd5YHDuW4tscAhaaakFZJezXL8FElzosczkg5KO/aWpLmSXpY0O6k8VoQ4VVVQWOBYtixsd9tt675x46CuDtra8k/POVcRkixxXA9M6eD4IuBDZnYgcBEwo93xD5vZwWZ2WEL5qwxxq6o6EzjSSxzjxoUqrHffzT8951xFSCxwmNmTwOoOjj9jZmuil7OAUUnlpaLlqqrqzLrj2QIHeHWVc91YVxkA+A3gwbTXBjwiyYCrzax9aeTfJE0HpgNUV1dTW6FzKTU0NGS8txEvvcQE4Nk5c9iydOl2xwfMn89hwKvPPsuqXvn9c+85axZjgL/Pm4dFXXD7rlzJEcAb99/P0vZrnHdCtvurFH5/5a3S7y9vZpbYA3gf8GqOcz4MzAN2Sdu3e7TdFXgFOCbO++29995WqZ544onMBy67zAzMVq/OfHz+/HD8xhvzf9MzzjAbOnTbfS0tZlVVZmedlX96Hch6fxXC76+8VfL9AbMtz+/2kvaqknQgcC1wkpm9l9pvZkui7QrgL8Dk0uSwDCTZqyp9DEdKz54wdqxXVTnXjZUscEjaA7gL+IqZvZm2v7+kgannwMeAjD2zHKFXVc+e0Lt35uPFDhzgXXKd6+YSa+OQdCtQAwyTVA+cD1QBmNlVwHnALsDvFOrKWyz0oNoN+Eu0rxdwi5k9lFQ+y15qZtxs7Q39+kGPHoUHjve/f/v948bBk0+GFQGL2M7hnCsPiQUOM5uW4/g3gW9m2F8HHLT9FS6jjtbigPDFXuh8VcuXZy5xjBsXRqKvWLHtGA/nXLfgI8fLXUer/6UUEjgaGsIjW1UVeHWVc92UB45y19EiTimFBI7ly8M2W4kDOhc4XnkFjj4a/vd/w3xYzrmy4YGj3OWqqoLCAkem6UZSRkVjNZcsyS/NdLffHpak/eEPYdQo9rjppsLTcs7tUB44yl2cqqpC1h3vqMSx004waBBkGHAY26xZcOih8NJLcNRRjLnuus6vje6c2yE8cJS7uFVV+U6rnmm6kXQjRhQeOFpb4fnnw8qCBx8MP/whamuD554rLD3n3A7lgaPcJVlV1aMHVFdnPj58+Nbgkq/XXw+B7IgjwusjjsCkUHXlnOvyCgockk4odkZcgZLqVbVsGQwbFgYXZtKZEsesWWGbChyDB7NxzBgPHM6ViUJLHBlGhbmSSLJXVbZqKtgaOMKcYvmZNQt22WVrt15g/X77wbPPhmos51yXVlDgMLPzi50RV6C4VVVbtkBzc/x0s003kjJiRCjtFLIk7axZobSRNup83f77w/r1oRrLOdel5Rw5Lumrmfab2Y3Fz47LS3NzeMQpcUAodey8c7y0ly2DffbJfjwVVJYu3Zp+HGvXhuAwbduJBdbtt1948swzcMAB8dNzzu1wcUoc7097fBC4ADgxwTy5uHItG5uS70SHZiFwdDSdyIgRYZtvO8cLL4Rtqn0jsnn33WHXXeHpp/NLL5OXX/auvc4lKGeJw8y+k/5a0mDgj4nlyMWXa0r1lHwDx7p10NSUu6oK8g8cs2aFKqr2kydKYSR5ZxvIn38eDj8cJkyAv/41bJ1zRVVIG0cjML7YGXEFyLXeeEq+gSPXGA7YGjjy7ZI7axZMnAiDB29/7KijYOHCrYMP82UG3/9+6A22ejVMngz33VdYWs65rHIGDkn3SronetwHvAH8NfmsuZziVlXlu+54nMAxdGhYAySfEofZ1obxTI46KmwLLXXceWeo6vrFL2D27DCn1oknhlKIc65o4kyr/uu05y3A22ZWn1B+XD7yLXHE7QHV0TxVKVIILPkEjrq6UBI4/PDMxw85JASjZ56Bz3wmfroQeo2ddRbsvz98/evQqxfU1sLIkfC734XSh3OuKOK0cfx9R2TEFSCpNo6O5qlKl+8gwDejhR733Tfz8b59w/xVqQGC+fjtb0NgevjhEDQg3Pcpp8ANN8Bll8GQIfmn65zbTqEjx2cUOyOuAEn1qlq5Mkw3MnRox+eNGJFfG0ddXdimDfzbziGHhF5RbW3x0zWDSy+FY4+Fj31s22Onnw6bNsHNN8dPL5M5c+AjH4GaGl+HxHV7hY4cv7qouXCFSapxfM2aEDR65PjzyLeqauHCsJRtRyWZSZNCldrChfHTfe01ePddmDp1+2OHHBIeM2YUNsq9oQG+972Qxpw5YR2Rgw6Cq64qLD3nKkDswCFpkKSBAGb2YnJZcrHFrarq0ydU38QNHKtXxxsoOGIEvPde6Lobx8KFMHZsx+uUT5oUti+9FC9NgEceCdv2pY2U008PX/qpMST5+I//gMsvD2m8+SbMnRu6DX/72yGgFIMZtLQUJy3ndoA4vaoOkzQXmAO8KukVSYfGSVzSTEkrJL2a5bgkXSFpgaQ5kg5JO/Y1SfOjx9fi3lC3EreqKt91x/MJHBC/++zChR1XUwHst18IcvkEjocfDu0mo0dnPv6lL4Xges018dME+Mc/4JZb4Kc/hd//Pnwmo0bBQw/BGWeEgHLvvfmlmbJ4MZx6Khx2WGh7GTgwpPn224Wl59wOFKfEMRP4TzN7n5ntCZwBXBcz/euBKR0cP54wJmQ8MB34PYCknYHzgcOBycD5knJUuHdDcauqINnAEae6yiy0cYwd2/F5ffqEnlFxA8emTfDkk9lLGxAWnZo6FW69Nf5n0NoK//VfoVfW2Wdve6xHj7Dk7aRJcNpp+a+E+NhjoerrjjvCmJOvfAVOPjkEtnHj4FvfKmwOMIB58+DnP4fvfAf++7/hxz8O77N+fWHpOZdBnO64G8zsqdQLM/uHpFj/+8zsSUnv6+CUk4AbzcyAWZKGSBoB1ACPmtlqAEmPEgLQrXHet9vYuDFMe15VlfvcfAPH3nvnPi99vqpcli0LX/K5ShwQvpDvuy8Em46qtQCeego2b4aPf7zj86ZPh5kzQ/CYPj13Hq6/Hv75z9Conikw9+kT0jrkEPjqV0N1Wa42IbMQcM46K8wDdtdd245sv+gi+NWvQvfhxx+H224L6eeyeXMIOtdcE6rSpFCKaW0N3ZS3bAl/IzU1ofT1uc/Fm19sxYqQj2efDR0C6upC1eSAAeH63XcPQf6AAxjQ1ARHHhk+l45s2hRKVanH2rUhSG7aFHrV9e8fAv3w4eGHyfDhYU2Y/v2z/y1s2hTytWZNmPVg/fqwr6UlfAZVVSHtvn1D3gcNCvnfaafQ5tanT0i7R4+tn1lTU0gjegx4441wXvpEoVJIu6oqdCPv0ydsq6pCqblnz615bmsLabe0bJ1jrrk57Gtt3bYzSI8e4ZFKI/2RyieEv6fUI5VG6n3a2rZvg0tdm9qOGJHfPHNp4gSO5yVdTfjSNuBkoDZVrWRm/yzonYORwOK01/XRvmz7XbrUWhy5vlwhv8CxZk3xSxypxu64geO668Iv+ZE5/tkfeST8Zz3mmI7PmzwZDjwwNJLnChzr1sFPfhLaMtpNxriNCRPgiivgm9+En/0sDDzsyAUXwIUXwuc/H+4vNTAzZfRo+L//C8dPOSUMlDz33DAavv25EL7cbrwxpLl4cTj/iitCYNh993BOS0v44r/33hCoTjstVImdcEK4v8MPD++7aVP4Ap8zJwyi/Mc/QqcDCH9j48eHEf/DhoW/uw0bwhf/449DUxOHQSgp7b13+DcbNizMDrBxYzh35UpYtCj730q/fiEAZutw0K9f+Bvu0yd8MTc1hXxs3Bi+6BN2WOLvUAK33w5f+EJBl8YJHAdH2/ZTqR9FCCQfKeidg0zfeNbB/u0TkKYTqrmorq6mtra2E9npuhoaGra7t73nz2eXqiqejXHPBzY303PJEl7KdW5rKx9au5a316/nrRznqrWVYyTenjWLtzqaSRfY7aGH2Bd4btUqNmVIN/3+BgOTgLk33sh7Rx7ZYbqH3XUXzfvtxysxGr5H1tQw/oormD1jBg0dlKjGXHste65YwYsXXsiGv+cYxjR2LHt/8pPs/stfMq+1leVTMtfMjrjmGrjlFpYefzxvfPvbYWR7B3pdeSV7X3opu55/Pk2XXso706ax7sADae3blx5btrDbY4+x26OP0nvtWtbvuy91v/41aw+Nmh7ffHPrmJmUT3wCjj+eQa+9xvCHH2aXxx+nz+23Z3zvlv79WT9xImtPP501hxxCw/jxWJYFvdTaSr/6enq+/jrDlixhp7feoveSJVTNm0evxkZa+/altV8/mgcNYvNBB7F5yhQ2jxjB5uHD2bzrrjQPGkRbnz7hF7AZPZqa6NXQQO81a+j93nv0Xr2aqrVr6b12LT02b6ZHUxM9Wlpoq6qitW9f2vr0oWXgQJoHDqRl4EBa+venpX9/2vr0wXr2xHr0QK2t9GxqoseWLfTctImejY3h0dQU0mxu3hqwJNqqqrCqKlp796Yteo+Nra30GTAgfA7RDzWZoZYW1NJCj5YW1NREj+Zm1NoaHumlCCnkJ3q09eyJ9eoV0uvRI6yAKYEZMoO2NtTWtjWd1DZ1PMqDwdbre/TAotLKv9NL/ahMK50o2q5rbWVLod+XZpboA3gf8GqWY1cD09JevwGMAKYBV2c7L9tj7733tkr1xBNPbL/zS18y22uveAl85jNm+++f+7z33gt/VpddFi/dXXc1mz4993nnnWcmmW3enPHwNve3fn0498ILO07z3XdDXv/nf+Lldc0as759zb71reznLFlittNOZlOnxkvTzKypyezYY82qqsza/zs1Npr94Achn6eeatbaGj9dM7NZs8yOOy69UiI8qqrMPvc5swcfNGtryy/NlPp6szvvNLvqKrMbbzT785/N5swxa2nJO6mMf58VpJLvD5hteX6vxylxJOke4ExJtxEawteZ2VJJDwO/TGsQ/xhwTqky2WXFWTY2ZdCgeA2kq1eHba7Bfylxx3IsXBiqRHLVgUOokhg/PncD+aOPhm1HDePphgwJjdA33wyXXJK5+ufnPw/VIBddFC9NCFUnf/5zmGvrxBPhy1+GT386fN4//CG8/TZLPvUpdr/22tztIO0dfniojnv5ZaivD1Uzra3hnocNyy+t9kaOhM9+tnNpuG4p0cAh6VZCQ/cwSfWE6q4qADO7CngA+ASwgDDr7mnRsdWSLgJS9Q8XWtRQ7tLEWTY2ZciQUHefy5o1YRt3wae4047E6YqbbtKk3FOPPPRQWMPjwAPjpzt9epiC5E9/gm98Y/s8zpgRxmyMGxc/TQif74MPhkBxww2h+y6EvNXW8qYZu2dbvz2Ogw8OD+e6gEQDh5l10LIIUTHpjCzHZhK6Arts4iwbmzJkSPgF3NbW8a/eVIkjn8CRakTtSF1d+DUe16RJ4cs9W9fg5uYQOD796fx+xR95ZBgrcskl4dd2esnqZz8LDe0/+1n89NLtuWcoeWzaFEpDjY2hoTs14aJzFSLOAMCdJP1M0jXR6/GSTkg+ay6nxsb8AodZ7uqqQgLHsmUdzy21YUPo2plrDEe6VDfUl1/OfPzpp0NXzk99Kn6aEBoLL7889PD55CdDT6LW1jDI79Zbw9iHVG+xQvXrF4Lk1KlbJ1x0roLE+al2HbAFSHVvqQd+nliOXHybNsUPHKmFk9au7fi8QgJHS0voR59NnMkN20tNPZKt99F994W2heOOi59mykc/GoLEc8+FEsuUKfDLX4Zuteedl396znUzcQLHXmb2K6AZwMw2kbm7rNvR8i1xQPzAEbdxPPXrvKPR0/mM4UgZNiy0D9x9d+bj990XBrQVOICJz342jKV47LEwiPAPfwgD6OI03jvXzcUJHE2S+hF1GZa0F6EE4kpt06ZQLRJH3MCxZk3obRRnNDpsrX5KlSoyKaTEAWHw3bPPwltvbbt//nx44438q6naS434nj07LP7knIslTuA4H3gIGC3pZuAx4MeJ5srFU0iJI1fPqrjzVKWkgkFHa1QsXBhKMPkupHTyyWH7pz9tuz+1jvgJRWhqO+64MG2Gcy62nIHDzB4FPgucSph25DAzq002Wy4ns2RKHPkGjsGDw1xC8+dnPyffrrh4Fr3bAAAVPElEQVQpY8aEaTRubTdF2X33hZ5RY8bkn6ZzrtPi9Ko6GthsZvcDQ4CfSNoz8Zy5jm3ZEoJHqQMHhDEPuUochQQOCNVVr7wSZn2FUGJ68snilDaccwWJU1X1e6BR0kHAj4C3gRsTzZXLbdOmsI1bVTVoUNjGaeOI2zCe0lHg2Lw5TIY3fnx+aaZ84QthnMZtt22dXbalxQOHcyUUJ3C0RAP1TgKuMLPLgQK7sriiSS3iFLfE0bNn6IGUVIlj8eKtwSzdvHlhnMQBB+SXZsqIEaH31E03hXEXF10UekTlmPzQOZecOIFjg6RzgC8D90vqSTRtiCuhfEscEKqrOgocZoUHDgiD6tqbOzdsCw0cEAbS1dWFKbyvvDIsTNSZ6Tucc50SZ1jrycCXgG+Y2TJJewCXJJstl1MqcMQtcUDu+aoaG8MEf/kGjlQ11Pz5Yc2GdHPnhrERhVZVQVh8aP78sFJeZwKQc64ocgYOM1sG/Cbt9Tt4G0fppaqqilniyHfUeEqqxJGpnWPOnLAeeGem3ujfP6yM55zrEuL0qvqspPmS1klaL2mDJF/AuNQKLXF0FDhSM+Pm2zg+dGgINpkCx9y5+c1e65zr8uK0cfwKONHMBpvZIDMbaGaDks6Yy6GQEsfgwcmUOCBzz6r33gtTrnv1knMVJU7gWG5m8xLPictPEiWOYgeOYjSMO+e6nDgVz7Ml/Qm4m7Q5qszsrsRy5XIrtI1j3brQe0oZ5qnsbOC47bYwMDE1UaAHDucqUpzAMYiwOl/6+pwGeOAopUJLHG1tYQ2KTLPK5rv6X7rx40PaixbBPvuEfXPmwC67dH59C+dclxKnV9VpOyIjLk/5DgCEbacdyRQ4Vq8Os+LmU4pJSe9ZlQocc+eG0kam0o1zrmzF6VU1StJfJK2QtFzSnZJG7YjMuQ4UMgAw12JOqcF/hXzRt++S29YGr77q1VTOVaC4KwDeA+wOjATujfa5UmpsDF/wvXvHvybXRIeFjBpP2WWXEJhSgeOtt8Ka6B44nKs4cQJHtZldZ2Yt0eN6oDpO4pKmSHpD0gJJZ2c4fqmkl6PHm5LWph1rTTt2T+w76i5Sy8bmUzpIMnBI2/as8oZx5ypWnMbxVZK+TFiLA2Aa0MEC00E0p9WVwHGEdcpfkHSPmb2eOsfMvpd2/neASWlJbDKzg2Pkr3vKZy2OlFyBY80aGNWJWshx4+Dpp8M06nPmhH2+SJJzFSdOiePrwBeBZdHj89G+XCYDC8yszsyagNsIM+xmM42twcnlks/qfym5VgFcvTr/UePpJk+G+voQQC64ICwrO2BA4ek557qkOL2q3gFOLCDtkcDitNf1wOGZTowWhhoDPJ62u6+k2UALcLGZ3Z3l2unAdIDq6mpqa2sLyGrX19DQsM29TXz7bQaY8Xwe96uWFj4ELHrpJd7OcN0HVq5kaWMjCwv9DCdNov/MmQyeO5dBr77K+v33Z0nMtNrfX6Xx+ytvlX5/+coZOCSNBS4HjiCM33gW+J6Z1eW6NMM+y3LuVOAOM2tN27eHmS2J3v9xSXPNbOF2CZrNAGYATJgwwWpqanJkqzzV1tayzb395jcwbBh53+9OOzFm6FDGtL+uuRkaGxl94IGMLtJnOBzYO+a5291fhfH7K2+Vfn/5ilNVdQtwOzCC0LPqz8SrUqoHRqe9HgUsyXLu1PZpmtmSaFsH1LJt+4crpI0Dsk87ktpXaOO4c67biBM4ZGZ/TOtVdRPZSw7pXgDGSxojqTchOGzXO0rSBGAooSST2jdUUp/o+TDgaOD19td2a42NxQ0cnZluxDnXrcTpVfVE1JX2NkLAOJmwEuDOAGa2OtNFZtYi6UzgYaAnMNPMXpN0ITDbzFJBZBpwW7Q8bcq+wNWS2gjB7eL03liOUOLYZZf8r8sVODrTOO6c6xbirgAI8B/t9n+dEEjGZrvQzB4AHmi377x2ry/IcN0zgA8A6EhnShwrVmy/f/nysK2ONUTHOdeNxelVNWZHZMTlKTUAMF9DhsCbb26/v74+bEeP3v6Yc86liTNX1RckDYyenyvpLkneUF1qhTaOZ1vMafHiMH3JsGGdz5tzrqLFaRz/mZltkPQB4OPADcBVyWbL5VTIAEDY2sZh7fo3LF4cRo33iPMn4ZzrzuJ8S6TGVnwS+L2Z/RXIY2Y9V3RmneuO29KydVr2lMWLvZrKORdLnMDxrqSrCdOOPBB1k/WfpaXU1BSmLS+0xAHbV1d54HDOxRQnAHyR0KV2ipmtBXYGfpRorlzHCln9LyXTfFWtrfDuux44nHOx5AwcZtYIrAA+EO1qAeYnmSmXQyHrjadkKnEsXx6qrzxwOOdiiNOr6nzgLOCcaFcVcFOSmXI5dKbEkWkVwMXRXJQeOJxzMcSpqvoMYXbcjfDvOaQyLFjtdphC1htPyVTiSAWOzqzF4ZzrNuIEjqZoOhADkNQ/2Sy5nApZbzylo8DhJQ7nXAxxAsftUa+qIZJOB/4GXJtstlyHOlPiyFZV1a+fT3DonIslzpQjv5Z0HLAemACcZ2aPJp4zl11nShx9+4ZHeq+qVFfcfNYvd851W3EmOSQKFI9CWEtc0ilmdnOiOXPZdaZxHLafdqS+3qupnHOxZa2qkjRI0jmSfivpYwrOBOoIYztcqXSmOy7ArrvCO+9sfe2D/5xzeeioxPFHYA1hgaVvEgb99QZOMrOXd0DeXDadLXF88INwww1huVgJli71wOGci62jwDHWzA4AkHQtsIqwDviGHZIzl11nSxzHHgu/+x089xzssUeYvsS74jrnYuqoV1Vz6omZtQKLPGh0EZ0tcdTUhFlw//Y374rrnMtbRyWOgyStj54L6Be9FmBmNijx3LnMGhtDFVOfPoVdP3QoHHYYPPoo7LNP2OeBwzkXU9YSh5n1NLNB0WOgmfVKe+5Bo5RSU6p3pvvssceGqqrXXguvPXA452Ly6dHLUaHrjac79tgwK+7NN8PAgVsHBjrnXA6JBg5JUyS9IWmBpLMzHD9V0kpJL0ePb6Yd+5qk+dHja0nms+wUut54uiOPDMFn0SIvbTjn8pJY4JDUE7gSOB6YCEyTNDHDqX8ys4Ojx7XRtTsD5wOHA5OB8yUNTSqvZacYJY6+fUO3XPDA4ZzLS5IljsnAAjOrM7Mm4DbgpJjXfhx41MxWm9kawqj1KQnls/wUo8QBoboKvCuucy4vsaYcKdBIYHHa63pCCaK9z0k6BngT+J6ZLc5y7chMbyJpOjAdoLq6mtra2s7nvAtqaGj4970duGQJPZubeamT9zpgyBAOAxa1tvJ2iT+39PurRH5/5a3S7y9fSQaOTF1+rN3re4FbzWyLpG8BNwAfiXlt2Gk2A5gBMGHCBKupqSk4w11ZbW0t/763Pn1g8GA6fa/HHAPLlzNm2jTG7LVXZ7PYKdvcXwXy+ytvlX5/+UqyqqoeSK88HwUsST/BzN4zsy3Ry2uAQ+Ne262luuN2Vo8ecO65UOKg4ZwrL0kGjheA8ZLGSOoNTAXuST9B0oi0lycC86LnDwMfkzQ0ahT/WLTPQWgcL0Ybh3POFSCxqioza4lm030Y6AnMNLPXJF0IzDaze4DvSjoRaAFWA6dG166WdBEh+ABcaGark8pr2SlWicM55wqQZBsHZvYA8EC7feelPT8HOCfLtTOBmUnmr2wVozuuc84VyEeOl6Nidcd1zrkCeOAoN2Ze4nDOlZQHjnLT3BzWz/ASh3OuRDxwlJvUIk5e4nDOlYgHjnKTWsTJSxzOuRLxwFFuOrv6n3POdZIHjnLT2fXGnXOukzxwlBsvcTjnSswDR7nxxnHnXIl54Cg33jjunCsxDxzlxksczrkS88BRbrzE4ZwrMQ8c5cZLHM65EvPAUW68xOGcKzEPHOXGSxzOuRLzwFFuUiWOvn1Lmw/nXLflgaPcpFb/k0qdE+dcN+WBo9z4WhzOuRLzwFFu3nsPdt651LlwznVjHjjKzcqVUF1d6lw457qxRAOHpCmS3pC0QNLZGY5/X9LrkuZIekzSnmnHWiW9HD3uSTKfZcUDh3OuxBILHJJ6AlcCxwMTgWmSJrY77SXgMDM7ELgD+FXasU1mdnD0ODGpfJYdDxzOuRJLssQxGVhgZnVm1gTcBpyUfoKZPWFm0cAEZgGjEsxP+TODVas8cDjnSqpXgmmPBBanva4HDu/g/G8AD6a97itpNtACXGxmd2e6SNJ0YDpAdXU1tbW1nclzl9XQ0MBT99/PB5ubWbBuHfUVdp8NDQ0V+28Hfn/lrtLvL19JBo5MAw0s44nSl4HDgA+l7d7DzJZIGgs8LmmumS3cLkGzGcAMgAkTJlhNTU2nM94V1dbW8sFRoUA27ogjGFdh91lbW0ul/tuB31+5q/T7y1eSVVX1wOi016OAJe1PknQs8FPgRDPbktpvZkuibR1QC0xKMK/lYeXKsPWqKudcCSUZOF4AxksaI6k3MBXYpneUpEnA1YSgsSJt/1BJfaLnw4CjgdcTzGt58MDhnOsCEquqMrMWSWcCDwM9gZlm9pqkC4HZZnYPcAkwAPizwhQa70Q9qPYFrpbURghuF5uZBw4PHM65LiDJNg7M7AHggXb7zkt7fmyW654BDkgyb2XJA4dzrgvwkePlZOXKsA6Hr8XhnCshDxzlxAf/Oee6AA8c5cQDh3OuC/DAUU5WroRhw0qdC+dcN+eBo5z4dCPOuS7AA0c58aoq51wX4IGjTPTYvDms/ueBwzlXYh44ykTV2rXhiQcO51yJeeAoE73XrQtPPHA450rMA0eZ8BKHc66r8MBRJjxwOOe6Cg8cZaLKq6qcc12EB44y0XvtWqiqgkGDSp0V51w354GjTFStXRtGjSvTworOObfjeOAoE1Xr1nk1lXOuS/DAUSaq1q71wOGc6xI8cJQJL3E457oKDxxloreXOJxzXYQHjnLQ1ESvjRs9cDjnugQPHOVg1aqw9cDhnOsCEg0ckqZIekPSAklnZzjeR9KfouPPSXpf2rFzov1vSPp4kvns8lauDFsPHM65LiCxwCGpJ3AlcDwwEZgmaWK7074BrDGzccClwP9E104EpgL7AVOA30XpdU8eOJxzXUiSJY7JwAIzqzOzJuA24KR255wE3BA9vwP4qCRF+28zsy1mtghYEKXXPaUChy8b65zrAnolmPZIYHHa63rg8GznmFmLpHXALtH+We2uHZnpTSRNB6YDVFdXU1tbW4y8dylDli5ltyOPZOGCBbSsWFHq7CSioaGhIv/tUvz+ylul31++kgwcmebGsJjnxLk27DSbAcwAmDBhgtXU1OSRxTJRU0PtIYdQkfcWqa2t9fsrY35/3UuSVVX1wOi016OAJdnOkdQLGAysjnmtc865EkgycLwAjJc0RlJvQmP3Pe3OuQf4WvT888DjZmbR/qlRr6sxwHjg+QTz6pxzLqbEqqqiNoszgYeBnsBMM3tN0oXAbDO7B/gD8EdJCwgljanRta9Juh14HWgBzjCz1qTy6pxzLr4k2zgwsweAB9rtOy/t+WbgC1mu/QXwiyTz55xzLn8+ctw551xePHA455zLiwcO55xzefHA4ZxzLi8KvV8rg6QNwBulzkdChgGrSp2JBPn9lTe/v/I1wcwG5nNBor2qSuANMzus1JlIgqTZlXpv4PdX7vz+ypek2fle41VVzjnn8uKBwznnXF4qLXDMKHUGElTJ9wZ+f+XO76985X1vFdU47pxzLnmVVuJwzjmXMA8czjnn8lJRgUPSJZL+JWmOpL9IGlLqPBWDpCmS3pC0QNLZpc5PMUkaLekJSfMkvSbpv0qdp2KT1FPSS5LuK3Veik3SEEl3RP/v5kk6stR5KiZJ34v+Ll+VdKukvqXOU2dImilphaRX0/btLOlRSfOj7dBc6VRU4AAeBfY3swOBN4FzSpyfTpPUE7gSOB6YCEyTNLG0uSqqFuAHZrYvcARwRoXdH8B/AfNKnYmEXA48ZGb7AAdRQfcpaSTwXeAwM9ufsDzE1NLmqtOuB6a023c28JiZjQcei153qKICh5k9YmYt0ctZhJUDy91kYIGZ1ZlZE3AbcFKJ81Q0ZrbUzP4ZPd9A+OLJuL58OZI0CvgkcG2p81JskgYBxxDW1cHMmsxsbWlzVXS9gH7RCqU7UeYrkZrZk4S1j9KdBNwQPb8B+HSudCoqcLTzdeDBUmeiCEYCi9Ne11NBX6zpJL0PmAQ8V9qcFNVlwI+BtlJnJAFjgZXAdVFV3LWS+pc6U8ViZu8CvwbeAZYC68zskdLmKhG7mdlSCD/kgF1zXVB2gUPS36L6xvaPk9LO+SmhCuTm0uW0aJRhX8X1oZY0ALgT+G8zW1/q/BSDpBOAFWb2YqnzkpBewCHA781sErCRGNUc5SKq6z8JGAPsDvSX9OXS5qprKLu5qszs2I6OS/oacALwUauMQSr1wOi016Mo8+Jye5KqCEHjZjO7q9T5KaKjgRMlfQLoCwySdJOZVcqXTz1Qb2apEuIdVFDgAI4FFpnZSgBJdwFHATeVNFfFt1zSCDNbKmkEsCLXBWVX4uiIpCnAWcCJZtZY6vwUyQvAeEljJPUmNM7dU+I8FY0kEerI55nZb0qdn2Iys3PMbJSZvY/w7/Z4BQUNzGwZsFjShGjXR4HXS5ilYnsHOELSTtHf6UepoMb/NPcAX4uefw34a64Lyq7EkcNvgT7Ao+HfmVlm9q3SZqlzzKxF0pnAw4ReHTPN7LUSZ6uYjga+AsyV9HK07yfRevWu6/sOcHP0o6YOOK3E+SkaM3tO0h3APwlV3y9R5lOPSLoVqAGGSaoHzgcuBm6X9A1CsPxCznQqozbHOefcjlJRVVXOOeeS54HDOedcXjxwOOecy4sHDuecc3nxwOGccy4vldYd17kdQtIuhAnhAIYDrYTpNwAazeyokmTMuR3Au+M610mSLgAazOzXpc6LczuCV1U5V2SSGqJtjaS/S7pd0puSLpZ0iqTnJc2VtFd0XrWkOyW9ED2OLu0dONcxDxzOJesgwnocBxBGyO9tZpMJ06x/JzrncuBSM3s/8DkqcAp2V1m8jcO5ZL2QmrJa0kIgNS33XODD0fNjgYnRNDkQJkMcGK1P4lyX44HDuWRtSXvelva6ja3//3oAR5rZph2ZMecK5VVVzpXeI8CZqReSDi5hXpzLyQOHc6X3XeAwSXMkvQ6U9YzOrvJ5d1znnHN58RKHc865vHjgcM45lxcPHM455/LigcM551xePHA455zLiwcO55xzefHA4ZxzLi//H7J0bBa/IcT+AAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "yout, trange = matlab.step(sysPI, arange(0, TIME_MAX, 0.1))\n", "plt.plot(trange, yout, 'r')\n", " \n", "plt.xlim([-2,TIME_MAX])\n", "plt.xlabel('Time')\n", "plt.ylabel('Response p.u.')\n", "plt.title('Controller behavior analysis')\n", "\n", "plt.grid()\n", "plt.show()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## ナイキスト線図・ボード線図" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXYAAAEICAYAAABLdt/UAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvNQv5yAAAIABJREFUeJzt3Xd4lFXaBvD7JPQmoYYO0gQhgISAoBCwK8padkUUsbJYEAWsILK6rnVd+7qsgo0FUcAV1KVoIkVa6EhAQIr0DiEESMj5/rjzfpOhpU3mnXnn/l3XXJCZIXOMyZ0z533Oc4y1FiIi4h1Rbg9AREQCS8EuIuIxCnYREY9RsIuIeIyCXUTEYxTsIiIeo2AXKSBjzAfGmGfdHofI2RjVsUu4MsZsAlAWwPnW2vSc++4DcIe1NtHFoSFnLJsA3Getnen2WCSyaMYu4a4EgEFuD0IklCjYJdy9BmCoMaZy7juNMe8ZY/5+yn1TjDGP5vy9nTFmiTEmzRjzhTFmvDHmrzmP3WWMmXPKv7XGmCY5f/8413OrGWOmGmMOGmP2G2NmG2OijDGfAagPYIox5ogx5oli+wqInELBLuEuBUAygKGn3P8JgNuMMVEAAxjAZQDGGWNKAfgawGcAqgD4EsDNhXz9IQC2AqgOoCaAZwBYa21fAFsAXG+trWCtfbWQn1+kwBTs4gUjAAw0xlR37rDWLgRwCAxzAOgNINlauwtAJwAlAbxprc201n4FYFEhXzsTQC0ADXI+12yrC1fiMgW7hD1r7SoAUwE8dcpDnwC4I+fvd4AzdACoDWDbKQG8uZAv/xqA9QCmG2N+M8acOgaRoFOwi1c8B+B+AHVy3fc5gF7GmDYAWoDLLwCwA0AdY4zJ9dz6uf6eDqCc84ExJvZsL2qtTbPWDrHWng/gegCDjTHOuwTN3MUVCnbxBGvtegBfAHgk131bwSWWzwBMtNZm5Dw0D0AWgEeMMSWMMTcBSMj16ZYDuNAY09YYUwbAyLO9rjGmpzGmSc4vicMATubcAGAXgPMD8d8nUhAKdvGS5wGUP+W+TwC0hm8ZBtbaEwBuAnAXgAMAbgUwKdfjv+Z8rpkA1gHwq5A5RdOc5x0Bf2G8b61NznnsJQDDcypmTr24K1JstEFJPM0Y0xVckmlorc0+x/M+BrDVWjs8WGMTKS6asYtnGWNKgpuXPjxXqIt4jYJdPMkY0wLAQbAU8U2XhyMSVFqKERHxGM3YRUQ8poQbL1qtWjXbsGFDN15aRCRsLV68eK+1tnpez3Ml2Bs2bIiUlBQ3XlpEJGwZY/K1Q1pLMSIiHqNgFxHxGAW7iIjHKNhFRDxGwS4i4jEKdhERj1Gwi4h4jIJdRM7q4EHg/feBhQvdHokUhIJdRAAA2dnAihXAW28Bk3K601sLPPww0Lkz71drqfCgYBeJcP/6F3DTTUD16kCbNsCjjwJf5xwiGBMDrFsH9OzJ+++4A0hPd3e8kjdXWgqISPBlZwOrVwPJycC2bcBLL/H+ceOAzZuBXr2AxESgWzegQQPfv2vcmDP4l18Ghg9nsDvBL6FJwS7icd99B4wZw0Dfu5f3NWkCvPACUKIEMHUqUKHCuT9HVBTwzDNAfDxQowbvsxbwOw5cQoaWYkQ8wlogNZUXO//0J2DXLt6/di0vfl53HQN+40Yur5TImdblFeq5XXkl0LYt/z5wIH85ZOtsqpCjGbtImFuzBnjuOc7Id+/mfXXrcnmlZk0G8GOPBfY1T54EDh8G3nsPWLQI+PRToHLlwL6GFJ5m7CJhwlqG+AcfAL17A198wftLlgTmzuVs+qOPgA0bgC1bgIQEPl6iGKZv0dHAJ58A77wDfP89X2vVqsC/jhSOZuwiIS4rC+jblzPynTt5X+3avMgJAOefD/z+e/DXu41hKWS7dsAttwCXX85lnrJlgzsOOZ2CXSREWAusX88AT0oCypcH/v1vzrj37AF69GDVSmIiL346Qe72BcwuXYAlS4BffmGoW8t19+hod8cVyQIW7MaYaAApALZZa3sG6vOKRILnnwdGjWIZIgDExgJ/+IPv8Zkz3RlXftWqxRvAuvgJE7hUVD3PQ9ykOARyxj4IQCqASgH8nCKeYS3w22+ckScnA/PnAytXAmXK8HbJJUD37pyRN2vm/ky8sMqXB+bNA9q3ByZOBDp0cHtEkScgwW6MqQvgOgAvAhgciM8pEu6s5S0qCpgyBXjoIa6FA6wFT0wEDhzgTPeJJ1wdakD17Qu0asXdrJdcwvLLe+91e1SRJVBVMW8CeALAWStajTH9jTEpxpiUPXv2BOhlRUKHtbx4OGYMcOed3L05dSofq10b6NSJIbd6NS+CfvGFb/nCa9q1A1JS+MvrvvuAZcvcHlFkKfKM3RjTE8Bua+1iY0zi2Z5nrR0FYBQAxMfHq5WQeMLx40Dp0sCOHQzuLVt4f7VqDLUqVfhx+/Zcd44kVaty1+uPP/o2NR07xmUnKV6BmLF3AXCDMWYTgPEAehhjPg/A5xUJOVu2sH777ruBRo24vALwYmdiIuu6V67krs8vv+RSRCSLjgauuIJ//+kn9p1JTnZ1SBGhyDN2a+3TAJ4GgJwZ+1Br7R1F/bwioeDAAXY4BFhumJTEv1epwjry7t35sTEMfDm7GjWASpVY7/7qq9wNG64XiEOd6thFctm61Ve1kpQEpKVx9m0McPPN7IDYvTsvDkZp33aBtGgBLFgA3HUXMGQIWxF8+CGraCSwAhrs1tpkAMmB/JwixWnbNvZTKVGCteTPPcf7Y2I4I09MBDIzgVKlfMsuUniVKrEE8pVXgGHDgK5dgQcecHtU3qMZu0SU7du51puUxFn5unWsue7UCbjqKgZPYiIQF6cZeXExBnjqKa69t2vH+w4eVBOxQFKwi6ft2ME/a9UCfv6Z298B4LzzOFscMACoX5/3dezImwRH+/b8c8sW9nl/6CHg2Wf1CzUQ9CUUT9mxAxg/noF9wQWsH3/nHT520UXA668DixcD+/YB33wDDB7M54h7qlcHrrkGGDkSuOEGzt6laIx14XTa+Ph4m5KSEvTXFe/ZtYvr5BddxMZTVasyGCpW5Iw8MZGhceGFbo9UzsVa4J//BAYN4sauyZOB1q3dHlXoMcYsttbG5/U8LcVIWNm9m2vkTtVKairQsiU7C0ZFscqifn2u3RZHH3IpHsYADz7Iw7RvuQX4+9+Bjz92e1ThS9/6EtL27OHa+A038Id/8GBg7Fge53bppSydS0z0Pf/mm90aqQSC0wK4YkV+vHMnd/Hql3TBaClGQsqBA9yC7tSSO6fyrF/PXYvLl3Nbevv2+mH3usxM/n+uWpV9dZxDtCNZfpdidPFUXLV3LzBpErBpEz+eMYNvxUeP5kXNl15iOaJTudKmDStXFOreV7Iku14uWMCAX7DA7RGFDwW7BNWxY7wwNmgQQ7p6dS6fTJ7Mx6+8kksvBw8C06ax3rlTJ/6QS+S54w5+P5QowYvh//632yMKD1qKkWK1fz8wezY7+l11Fbfox8RwJ2eXLr6j3jp04H0iZ7J/P9CnD6ugFiyI3O8VVcWIa/73P862k5O5Jm4tdxledRUviqWksJIlUn84peCqVAG+/ZYBX6oUcPgwcOgQUK+e2yMLTVqKkSI5eJCnA73+uu++f/wD+OAD/jA+/zxn7FOm+B5v21ahLgUXHe07Q3XgQO5dcLptij8txUiBpaRwd2dyMrB0KTcGlSvHGvPy5blhqFo1HkAhUhzWrgVuvJF/vvoqy2AjoQWwqmIkIA4d4lvgxx9nAy2Aa5zvvstllREjGPD79vnar9apo1CX4tW8Ob8Pb7wRGDoU6N0bOHLE7VGFDq2xy2m2bgXefptvc5cs4Yy8VCngsstYgnjXXTycWEeciZsqVuQpVa+9xuW/Awe4cU20FBPx0tKAOXM46+7YkSfLb9kCNGnCMkOnauXii4GyZV0erMhZpKUx6LOzgfnzgc6d3R5R8VBVjJyVtTzk4McfuV5+8iTrxJ98ksFevz6XYBTkEi6cFgSjRvHgjhEjeGhKpLYAVrB73JEj3OCRnAwcP87mSsYAM2cyzJ96ijPyzp15AdShUJdw1K8f196ff56Tls8/951ZG0kU7B710Ue8LVoEZGVx516PHpytG8Nt+tHRbo9SJLDKlmU7ioQE7m6Oj+eu5rg4t0cWXBH6RsU70tPZX2XYMHY7dCoDnJODHn+cm4WcLfpOSZhCXbzKGC7HJCdzmTESq2V08TRMzZ3LNfGFC9kFLzqas5PPP+eFTxEBTpzwbYb79lv2IgrnvkOqY/eIjAxe5Bwxgk2Qvv6a95crx0B/7DHg++9Z6jV/vkJdJDcn1FesAHr2BC6/nP1mvE5r7CHq0CEeLjF/PmcdUVG+w38BnhCkNqYi+RMXx3ez99/PVgQTJ7Kc16s0Y3fZsWM86m3kSFan9O/P+ytV4k7ORx4Bpk5l86OFC4E//MHN0YqEr9tvZ9FA6dJ89/vhh26PqPhoxh5kJ0/6Llzeey+PeTt+nDPydu2A88/nY8YA333n3jhFvKhNG5ZB3n672yMpXgr2Ynb8OGfazlFvqanA778z3Js3Bx56iDP1Sy8FKld2ebAiEcBpAexsXvr+e+DCC32ndHmBgj3AnPXwEiV42sugQbwAagzb1fbuDRw9yp1yTzzh9mhFIpMT6hkZfOecmclzVXv0cHdcgaI19iI6cYKlhy++yCvulSvzY4CzgD//mZUse/eyodYbb/i2P4uIu8qW5TvpGjV4GMzrr3MTX7jTjL2AMjO5KahyZeCXX7jD7ehRPhYXx6vuzmEAnTt7txmRiFc0a8YKs3vu4Ya+hQuB//wnvA9MD+OhB0dmJrB4MVvYJiezE+KAAey50rQpg7xrV6BbN6BqVbdHKyKFUaECl2ISEnjuQDiHOqBgP01WFvuRN2zIt2TNmwMbN/KxVq24HtezJz8uVQp4803XhioiAWQMD+1wLF7MLOjVy70xFVbEB3tWFte+naqV2bOB2Fhg3Tr+j37mGXaH69rVt8QiIt73t78BkyYBw4dzn0k49VeKuGDPygKWLeMuTmNYbjhqFB9r2RK4806WHzpdEO+7z9XhiohLxo4FHnwQ+OtfWfs+dixLJcOB54PdCXJnjXz2bJ62sno10KIFcPfdPPKtWzegZk23RysioaJMGba+7tgRGDgQ6NCBGVKvntsjy1uRg90YUw/ApwBiAWQDGGWtfauon7ewTp5kkMfG8lDlb7/1bcNv3hzo0wfo3p2PAewX4eWeESJSeMawZLlNG+C994BatdweUf4EYsaeBWCItXaJMaYigMXGmBnW2tUB+Nx5ys5mkDtr5LNmsYHWSy/xdKBu3YBx4/hnuPxPEZHQknsCuGsX8O677Lgaqi2Aixzs1todAHbk/D3NGJMKoA6AYgn27Gy24ExPB7p04Zb9iy/mRqGmTYFbb+UaubODrHJl7vYUEQmEyZO57v7TT8CECVwdCDUBXWM3xjQE0A5AsTSUnTSJFzMPHGCY//wzd45NmcJdns7yiohIcRkwgN1X77uPRRhffcU8CiUBaylgjKkAYCKAR621h8/weH9jTIoxJmXPnj2Feo2YGF+Ph5o1OXsHeCqKQl1EgqVPH56VUKYMl3mnTnV7RP4CEuzGmJJgqI+11k4603OstaOstfHW2vjqhSwI796d9eUDB7LsyAl5L/R2EJHwEhfHMsi77w691iFFDnZjjAHwEYBUa+0bRR/SucXEAG+/7WuSv3QpT0RJTi7uVxYR8RcTA/zrX5xoHj/OCprNm90eVWBm7F0A9AXQwxizLOd2bQA+7zkZwz8PHeKae/fuwB//GBpfVBGJPKmpwPjxXHefOdPdsRQ52K21c6y1xlobZ61tm3ML2tk/iYn8gv7lL6xZv+AC4JVXgvXqIiLUti2XZmJjgauuAl591b1lYk/0Yy9bljWla9bwAOhjx3i/tVp/F5HgadqUF1VvuQV48km2AXaDp1oK1K/P1ptOmH/zDfDWW7y1bu3u2EQkMlSowCWZTp1YsecGT8zYT+Wsv2dkAMuX85DogQOB/fvdHZeIRAZjgMce4/4aa3lE5tdfB+/1PRnsjt69WR45YADw/vs8KWXsWLdHJSKRJD0dmDcPuPFGYNiw4CwPezrYAZYhvfsuyyJbt2aTMEBr7yISHOXKccYeHc0e7+vXF/9rej7YHXFxwI8/An378uO33wZuuw34/Xd3xyUi3nbiBBsSNmsGTJzIC6zFLWKCHeC6l7P+npXFNa8LLmBDn4wMd8cmIt4xaxYbEh4/zrYDSUnAypXATTcF5/UjKthzGzKE5ZHXXgs8+yxPT/rhB7dHJSLhLCWFNezduvFQn3XreH+TJsE9Wi9igx0AGjQAvvySSzQVK/p6K2v9XUQK4vBhzsY7dOAh2K+9BmzYALRq5c54PFXHXljdu/OwDqep2NChXKoZOZK9IEREzuTIEdatV6zIcP/LX4BHH2VbXzdF9Iw9t9ydIk+cAN55hxc7Ro3yVdKIiADAtm0so27QANi3j9fuZszgDni3Qx1QsJ/GGIb6kiU87PrPf+bbqxUr3B6ZiLhtzx5en2vcGBg9mpV1ztKtU5gRChTsZ9G2LY++GjeOb7HKlXN7RCLipt27GehvvslA//VX7pGpVs3tkZ1OwX4OxnD36tq1vKoN8H/o3/7mazQmIt6Vns6usQBQowbX0H/5BRgzBmjY0NWhnZOCPR+cMqVjx1iXOmwYe0B8840qaES86PhxLsk2bsyOsc5Gxsce496XUKdgL4AyZXig9owZ/HuvXsDVV2v3qohXZGVx7bxZM+CRRxjis2cD9eq5PbKCUbAXwuWXszzyzTeBjRtZ7iQi4W/nTuDBB4GaNYHp07ljNNTOM80PBXshlSzJxj6pqax1z8pi7+XRo4HsbLdHJyL5YS0wZQrw0EP8uG5d7h5dsAC44orQqnQpCAV7ETnr7/v2AUePAvfeC3TsyDadIhK6fvyRs/EbbgCmTQP27uX9rVqFb6A7FOwBUrMm1+LGjgW2b+c3zJ138rBtEQkdW7ZwOfWyy4CtW7kJMTU1NMsWC0vBHkDGAH36sDzy6aeBVatU/y4SKo4e5Z9VqgC7dvEa2bp1wP33+/pEeYWCvRhUqMBa94UL+Q1z+DDQtSswdarKI0WC7ddfuf+kfXteC6tQgTvJBw1idZsXKdiLUYmcFmtbt3LX2vXXA9ddxxm9iBSvLVuA++5jS+5vvuHRdCdO8LFwX0PPi4I9CFq25Azh738H5s7lEX2PPw5kZro9MhFvcd4Rr1vHk4o++wx4+GHgt9/4LjpSlkYV7EFSqhQweDDfFvbty/X3EmqaLFIkmZnAnDnsqtixIxt0AUCjRmzgt24d19Jr1nR3nMGmaAmymjWBjz7iN6QxwObNQL9+PMWpSxfvrvmJBFr//sAXX/AaVlQUg71ZMz5WogTPNY5UCnaXOFfhN27kmvvll/O+Nm2AhAQ2G/JS+ZVIYaWlcQfo9Ok8znLGDE6KKlXiuaJXXQX06KFDcXIz1oUyjfj4eJuSkhL01w1VR47wm3bhQt6WLWMj/7JledB2UhLDvmNH/lm7ttsjFil+337LI+bmzmU1S7lyPO1s/PjIbeNhjFlsrY3P63masYeAChV4XqJzgrm1vqv2lSpxk9Prr/ObG+DFV+fgjzVrGPShcGqLSGFt387JzfTpwPDhLDhIT+cyy5AhnJV37gyULu32SMODZuxhIiODM/mFC/n3p57i/RdeyF1zLVpwNp+QAFx6qXuH6Irk1969wEsvMcxXreJ9NWuy39K117o7tlCV3xm7gj3MTZ/OhkULFjD09+xh1c2nn3Lm//TTnOEnJPCwEK/X70ro2b8fWLTId0tMZF/zw4cZ5F26cEZ+5ZVAXJy+R89FSzER4soreQMY5Js3+5Zsdu3iYQHOVuqYGAb8Y4/xB0kk0NLT2fq2cWN+P7ZvDyxd6nu8eXO+owS4fHjwoJZXioOC3UOM8T+uKzaW6/OrV/suzC5YwCoDgB//6U/+F2bbtgUqVnRl+BKG5s0Dfv6Z13yWLuWxcW3a8DB4Y4CePX3fY+3bA+ed5//vFerFQ0sxEWzpUuDllxnwmzb57p8/n0E/Zw4rExo1As4/n3/Wr++9hklybunpvEi/ahUDfPNm4Kuv+Fjv3qwlr12byygdOgCdOmmNvLhoKUby1K4dfygBLtssWgSsXOnb5LF0qX81DsCNIFu3ArVqMfSXLGHgO+EfG6s10nCVlsYAT03lLLtMGZbbjhjh26pfpgwvzKenA+XLA6+8wuW+6tXdHbv4C8iM3RhzNYC3AEQD+NBa+/K5nq8ZexCMHMlbEZ08ySDfuNF3GzmSAT9kCPDGG/7Pr1gROHCAB5CMGcMWCrGxvEgWG8uZnfOLQ9xx5Ah3ZpYpA/z0EytTVq/2P7t32TIuqSQn85yBli1ZgdW0qe9wGQm+oFXFGGOiAfwK4AoAWwEsAnCbtXb12f6Ngj0IjAlKj+CMDL41/+03hv6BA6xDBtgqYdw4/2ZnDRvyec7jq1cDNWqwR3ZMDC+uOceUzZvHYwYrVeLabKVK/MWhYCmY334DvvuOG30WLwbWr2cL6Wuv5SlCQ4cyuFu2ZNlsy5a8+KleRqEnmEsxCQDWW2t/y3nh8QB6AThrsIt3lC3Lk9wvuOD0xz75BPj4Y4b9zp1c7nHapgIM+V27eEtNZVlcq1YM9uxsHjOYmur/OVu3ZlXPkSN8t5CWxqCPiuLvspo1GU5ZWXzus88W53996Lv1VmDCBP69bl2ugffty9JXgFvxlyxxb3xSPAIR7HUA5HoTh60AOp76JGNMfwD9AaB+/foBeFk5zciRbDLjcBa7n3suIMsyhWEMZ+OVKgGVKzPEv/+eQe/8Uti7l6G+fz8fr1aNvwzOdCj4ypXAPff43xcV5Qv29HTfUsOpFRiRqHJl/pn7F/AFF/DNXFaWZuVeFYilmD8CuMpae1/Ox30BJFhrB57t32gpJgiCtBQDcKnl99/91+G3bGF4O7d9+848nAoVeOGtSpXTbzExDKZTl2IqVOCFu3LlGFi6WHt22dk8qHnqVJYlrlrluxheujRD3lmGcZZimjRR5VOoCuZSzFYA9XJ9XBfA9gB8XgkR1nImnTu4nTX1jRsZ6idP+p4fHc23/bVqMSS6dOGF09wXUWvW5K18eff+uyJBVBRwzTW8AcCxY6w1X7WK735Wr2bgjxvn+zclSvACt3PB9MIL+fdmzRT44SIQwb4IQFNjTCMA2wD0BtAnAJ9XiuK55wr0dGuBDRv4A587tJ1bRob/82NjWeLYubN/uWOjRgx1vcUPTWXKcKNQ+/b+9x85wvbRq1f7bsuWARMn+t5plSzJGX1cnP9NJa6hJ1DljtcCeBMsdxxtrX3xXM/XUoz7du9m3bqzI3XhQq5xOypV8gV27tBu1IgXPSPliLFIl5HB2vbVqznDX7mSm5S2bvU9p1o1BnybNkB8PHeZNm6ssC8OagIm/y89nZUPuUPc2WkaFcVKFKczZNu2/KGMidEPppzd/v2+kM99O3aMj8fEMOQ7dOD3VYcOOkcgEBTsESory783zMKFXF5x1sAbNPCFuNO/Q+vcEghZWVy/d94JOjuZne+92rV9Id+hA4Nfpx4VjII9AljL6pPcDb4WL/Z1c6xc2RfgHTvyhynSDvUVd+U+R8Bp2/vrr77HmzfnsZBXXsnTkdSA7twU7B7k9LXOPRvfvZuPlS7N3i+5Z+Pqvy6h6OBBICWF38tz5rBtwdGjvODeuTNbSl91Fb+fo6LcHm1oUbCHuWPHfDMd57Zune/x3CcmJSTw4lWpUu6NV6Swjh9nyeW0aTw4xunfXq0acMUVvjMHtEavYA8r1rLywFlOWbgQWL7ct5GkVi1fv/SOHc/c11rEK3btAmbO9AX9rl28v1Ur30lLl17KzWmRRsEeJhYsAB55hGEOcI0xPt4X5AkJQJ067o5RxC3Wstpm+nQG/ezZ7DdUpgzQtatv2aZly8hYdlSwh7gdO3ge6SefcIPH8OG8eNS8uboXipxNejowa5ZvNu80iatTx7dkc/nlXMbxIgV7iDpxAnjrLeD557m2OHgwMGyYqgFECmPLFmDGDAb9zJlsHmcMlyv79AEGDPDWkk1+g13XnIPou+/YSvaJJ4Bu3Vjz+/LLCnWRwqpfn+2dJ0wA9uzhsY4jR/Jd7+DBrAz74AP/dtGRQMEeBOvW8VDf667jx99+y257TZu6Oy4RL4mO5rWpESMY8MnJbIHxwAPsYvnpp/7N6rxMwV6M0tKAJ59kd7xZs4DXXuNOPB30K1L8unXjxdbvvuNmvX79WBacu7GZVynYi0F2NmcHzZoBr74K3H47d9sNHapac5FgMoYti1NSgC+/5M/mLbdwF/b//ufdgFewB9iiRdw9168fUK8e3xKOGcPKFxFxR1QUA33VKh7XuG8fA9+Z1XuNgj1Adu3iRZyEBHZOHDOGod7xtEMCRcQt0dGcdK1dC7z3Hg/27toVuPpq9lnyCgV7EZ04wUOVmzXj8svQoVx2uesu9bkQCVWlSgEPPshgf/VVvtOOjwduvpndUcOdoqcIpk3j4QJDhnD5ZdUqXiCtVMntkYlIfpQrBzz+OE8JGzmSNfGtWgF33smTxMKVgr0QNmwAevXi27esLGDKFF55b97c7ZGJSGFUqsTTJDdu5LvuL7/kz/MDDwDbtrk9uoJTsBfAkSPAM8+wL8UPP3Bz0apVrFGPhD4VIl5XtSqXZjZsAPr3Bz76iJuchg4F9u51e3T5p2DPB2uBsWP5G/yll4Bbb+U6+pNPsg+6iHhL7dq8uLp2LX/e//EPbnYaMQI4dMjt0eVNwZ6HxYuBSy4B7riD7XN//pkXSdUbWsT7GjVieeSqVVx6feEF3vfKK76TykKRgv0s9uzhW7EOHdgS4MMP2Vr34ovdHpmIBFuLFlx3X7KEGfDUUzz0/d132cwv1CjYT5GZye6LTZuyFv3RR7nscu+9Kl8UiXTt2rHX05w5LHEeOJB/jhnjOxgnFCiqcpkxA2jblmHesSMb/L/xBvtMiIg4unRhk7Fp04AaNYB77mGZ5IQJbFvgNgU7WK8lbbOyAAAQC0lEQVR6441s0n/sGPD11+wj0aKF2yMTkVBlDDNj4UJg0iQexn3rrcBFF3FW72YfmogO9vR0nlzUsiVPY3nxRfZI79VL5Ysikj/GcGK4fDnw+ecsi+7Zk7P6pCR3xhSRwW4tMH48ezS/+CK3Ea9dyxr1MmXcHp2IhKPoaHZyTU0F/vUvnu7UowdwxRW+M42DJeKCfdkydnS77TagenV2dhs7Fqhb1+2RiYgXlCzJirr163mNbvlyXrP7wx94HkMwREyw793L8w/bt+dv1FGj2PjnkkvcHpmIeFGZMsBjj/Ea3gsv8NS0du2C04PG88GelQW88w7LFz/8EHj4YZYv3n8/3zqJiBS3lSt5LF+vXjynobh5Oth//JG/IR95hDP15ctZox4T4/bIRCQSrFsHdOoEfPUVe0t99RWXaoqbJ4N90yaelnLZZbxCPXEia9QvvNDtkYlIpJgyhT3ed+5kvfuTTwav2s5TwX70KFtvtmjBNrrPP8+m+TfdpPJFEQmO7Gzm0A03sDPk4sXA5ZcHdwwlgvtyxcNavsUZMgT4/Xegd2+23gzGWpaIiOPAATYM/O47nqL2/vtA2bLBH0fYz9hXrAC6dwf+9CegShXgp5+AceMU6iISXCtWcOllxgwG+ujR7oQ6UMRgN8a8ZoxZY4xZYYyZbIwJWleVffuAhx7ixdGVK4F//pNvebp2DdYIREToP//hRdKMDPaQeeABd5d/izpjnwGglbU2DsCvAJ4u+pDOLSuLvw2bNQM++IBfwHXrWKOu8kURCabMTNaq3347Z+tLlvD8Y7cVaY3dWjs914fzAdxStOGc2+rV3DG6YgV3j779NhAXV5yvKCJyZrt2cQl41iyWVL/+enBKGfMjkBdP7wHwRQA/32l++omhDrA1wLBhXF9PTATatNGMXUSCY/589pg6cAD47DNeMA0lxubRW9IYMxNA7BkeGmat/W/Oc4YBiAdwkz3LJzTG9AfQHwDq16/ffvPmzYUa8I4dDPjkZHZO+/VX3n/eeVxfd4I+Lk5BLyKBZS3bkQwcyP5SkydzUhksxpjF1tr4PJ+XV7Dn44X6ARgA4DJrbb5OAYyPj7cpKSlFel3H9u0M+qQkhv26dbw/JoZBn5joC3qdgCQihXXsGAs2Ro/m+adjx7ISL5iCEuzGmKsBvAGgm7V2T37/XSCD/VTbtjHgndv69bw/Jobr8k7Qt26toBeR/NmyhUsvKSk8w2HkSHdWBIIV7OsBlAawL+eu+dbaAXn9u+IM9lNt3eof9Bs28P4qVfyDvlUrBb2InO6HH7jp8cQJ4NNP2cjLLUFbiimMYAb7qX7/3T/onRaaVav6gr57d56qpKAXiVzWstLlqad4KM+kSUDz5u6OScGeT5s3+y7GJicDGzfy/mrVTg969ZsRiQxpaTyg+quv2FBw9GigYkW3R6VgL7RNm/yrbpzinerVGfRO1U2LFgp6ES9au5aNA9esYavdoUND52c9v8HuiSZggdSwIW/9+vHjTZt8IZ+UxN/gAFCjhn/QX3BB6PzPF5HC+e9/gTvv5Eaj6dPZ+jscacZeANYy6J3SyqQkXpwFgJo1fRdiExO5FqegFwkPJ0+y0uWvf+WhPBMnAg0auD2q02kpJgis5Zp87qDfto2Pxcb6B32zZgp6kVC0fz97vfzvf8Ddd7MXVZkybo/qzBTsLrCW5ZTOhdikJG6gAoBatfyDvmlTBb2I25YvB268ke+833kH6N8/tH8uFewhwFpukMod9Dt28LHatf2DvkmT0P6GEvGazz9nkMfEcOmlUye3R5Q3BXsIspYtD3IH/c6dfKxOHV/Id+8OnH++gl6kOGRmstLl7bfZdmTCBF4jCwcK9jBgLZuYOSGfnMxWoAAbDOUO+kaNFPQiRbVzJ1vtzp4NPPooj9AMlVa7+aFgD0PWsoY2d9Dv3s3H6tXzhXxiIksyFfQi+TdvHjcbHTgAfPgh0KeP2yMqOAW7B1jLTRJOyCcnA3tyWq3Vr+8LeSfoReR01vK0tUGDOEGaPDl8D+hRsHuQtTxFKnevm717+ViDBv5BH4o1uCLBlpEBPPgg8PHHwDXXsNVuTIzboyo8BXsEyM4+Pej35fTZbNTIv+qmfn23Rinijs2b2RpgyRJgxAjguefCv7Gfgj0CZWcDv/ziH/T79/Ox88/3D/p69dwapUjxmzmTrXYzM1nWeP31bo8oMBTsguxsYNUq38XYn37ihSMAaNzYv+qmTh0XByoSINay0uWZZ9iob/Jkbgb0CgW7nCY7G1i50j/oDx7kY02a+EK+WzcFvYSftDS2BJg4kSWNH30EVKjg9qgCS8EueTp5kkHvVN3MmuUL+qZN/YO+dm03RypybmvWcD197VrO2AcP9mY5sIJdCuzkSWDFCv+gP3SIjzVr5qu66daNvW9EQsHXX7PVbunSwBdfAD16uD2i4qNglyI7eRJYtsx3IXbWLODwYT7WvLl/0MfGujhQiUgnT7La5W9/Azp04BKM14sCFOwScCdPAkuX+gd9Whofa9HCdzG2W7fw6b0h4WnfPrbanTYNuO8+dmYM1Va7gaRgl2KXleUL+qQk9t84coSPtWzpH/Q1arg4UPGUpUu5nr59O/Duu8D997s9ouBRsEvQZWVxM4gT9HPm+IL+wgv9g756dRcHKmHrs8/YardqVS69dOzo9oiCS8EursvMPD3o09P5WKtW/kFfrZqLA5WQd+IEMGQIZ+jduvEiaSQu9ynYJeRkZgKLF/uqbubMAY4e5WOtW/vKK7t25YxMBODhNH/8IzB3LssYX3kFKFHC7VG5Q8EuIe/ECSAlxXcxdu5cX9DHxfkHfZUqLg5UXDN3LlvtHj7MDUe9e7s9Incp2CXsnDgBLFrkH/QZGdxoEhfnK6/s2jW8O/RJ3qzlodKPPspOpZMn811dpFOwS9g7fvz0oD92jEHfpo0v6C+9VEHvJRkZwIABwKefAtddxyZelSu7ParQoGAXzzl+HFi40Hcx9uefeZ8xQLt2vouxl16qIAhXmzaxlHHpUmDkSODZZ8O/1W4gKdjF844dY9A7F2PnzWPQR0WdHvTnnefyYCVP06cDt93GjXCffw707On2iEKPgl0izrFjwIIF/kF/4gSD/qKL/IO+UiWXByv/z1rg5ZeBYcO432HyZHYbldMp2CXiZWT4B/38+b6gb9/eP+grVnR5sBHq8GHgrrsY5r1785Dp8uXdHlXoUrCLnCIjg7N452Ls/PmsrY+O9gV99+5Aly4K+mBYswa48UZg3TrgtddYAePFVruBpGAXycPRo/5Bv2CBL+jj4/2D3msHNrht8mS22i1bFpgwgV9ryZuCXaSA0tN9QZ+UxAuzWVnc5Rgf7yuv7NxZQV9YJ08Cw4dzTT0hgf1e6tZ1e1ThQ8EuUkTp6SypdIJ+0SJf0Hfo4B/0WhfO2759rHqZMYONvN5+m4djSP4FNdiNMUMBvAagurV2b17PV7BLODpyhEHvXIxdtIgz0JIlOft0LsZ27gyUK+fyYEPMkiWsT9+xA3jvPfZQl4ILWrAbY+oB+BDABQDaK9glUqSlcTess0afkuIL+o4d/YO+bFl3x+qmTz7hTtJq1bj0kpDg9ojCVzCD/SsALwD4L4B4BbtEqrQ0dqzMHfTZ2UCpUv5Bf/HFkRH0J06wG+N773HZavx4HbhSVEEJdmPMDQAus9YOMsZswjmC3RjTH0B/AKhfv377zZs3F/p1RcLB4cP+Qb94sS/oO3XyVd106uS9Y922b2er3Z9/BoYOBV56KXJb7QZSwILdGDMTwJmOKh4G4BkAV1prD+UV7Llpxi6R6NAh/6BfsoRBX7q0f9B37BjeQT9nDkM9LY2tdm+91e0ReUexz9iNMa0B/AAgp4M26gLYDiDBWrvzXP9WwS4CHDzoC/qkJDa+spZBf/HF/kEfDtUj1vKEo8GDgYYNWaveqpXbo/KWoJc7asYuUjQHD/JAcKfqZtkyhmWZMgx6p7wyISH0gv7oUV4g/ewz4Prr2XJXHTYDL7/BrlUvkRBRuTJD8frr+fGBA/5B/9xzvqDv3Nk/6EuVcm/cGzeylHH5cuD559nMS6123aUNSiJhYv9+YNYs3xr98uW8v2xZtj1wqm46dAhe0E+bxk1H1gJjxwLXXhuc141U2nkq4nH79vkH/YoVvL9cudODvmTJwL52djbbAgwfznX0yZOBxo0D+xpyOgW7SITZu9c/6Feu5P3lygGXXOK7GNu+fdGC/vBhoF8/4OuvOVv/97/VUiFYFOwiEW7PHl/QJyUBv/zC+8uX9w/6iy7Kf9CvXs1Wuxs2AK+/DgwapFa7waRgFxE/u3f7B/3q1by/QgUGvXMx9qKLzryZaOJEHopRrhxb7XbrFsTBCwAFu4jkYdcuBr1TdZOayvsrVvQP+rg4VuS88go3Un35pVrtukXBLiIFsmuXb30+OZknHAFcpsnMZJ36m2+GXg19JFEdu4gUSM2a3P7vtADYuZMBP3cuN0j16ePq8KQAFOwickaxsTxgundvt0ciBaX9YSIiHqNgFxHxGAW7iIjHKNhFRDxGwS4i4jEKdhERj1Gwi4h4jIJdRMRjXGkpYIzZA2BzMX36agDyPJ4vAujroK+BQ18H8sLXoYG1tnpeT3Il2IuTMSYlP70UvE5fB30NHPo6UCR9HbQUIyLiMQp2ERGP8WKwj3J7ACFCXwd9DRz6OlDEfB08t8YuIhLpvDhjFxGJaAp2ERGP8WSwG2NeM8asMcasMMZMNsZUdntMbjDG/NEY84sxJtsYExFlXg5jzNXGmLXGmPXGmKfcHo8bjDGjjTG7jTGr3B6LW4wx9YwxScaY1JyfhUFujykYPBnsAGYAaGWtjQPwK4CnXR6PW1YBuAnALLcHEkzGmGgA7wG4BkBLALcZY1q6OypXfAzgarcH4bIsAEOstS0AdALwUCR8L3gy2K210621WTkfzgcQkWeqW2tTrbVr3R6HCxIArLfW/matPQFgPIBeLo8p6Ky1swDsd3scbrLW7rDWLsn5exqAVAB13B1V8fNksJ/iHgDfuz0ICao6AH7P9fFWRMAPs5ybMaYhgHYAFrg7kuIXtodZG2NmAog9w0PDrLX/zXnOMPCt2Nhgji2Y8vN1iEDmDPeprjeCGWMqAJgI4FFr7WG3x1PcwjbYrbWXn+txY0w/AD0BXGY9XKyf19chQm0FUC/Xx3UBbHdpLOIyY0xJMNTHWmsnuT2eYPDkUowx5moATwK4wVp71O3xSNAtAtDUGNPIGFMKQG8A37g8JnGBMcYA+AhAqrX2DbfHEyyeDHYA7wKoCGCGMWaZMeYDtwfkBmPMjcaYrQAuBvCtMWaa22MKhpwL5w8DmAZeLJtgrf3F3VEFnzFmHIB5AJobY7YaY+51e0wu6AKgL4AeOVmwzBhzrduDKm5qKSAi4jFenbGLiEQsBbuIiMco2EVEPEbBLiLiMQp2ERGPUbCLiHiMgl1ExGP+D7eTCcQ3l+hHAAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZEAAAEOCAYAAABIESrBAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvNQv5yAAAIABJREFUeJztnXeYXFX5+D/vzPaSTS+kbRKSYAoQUiihJIoSEKQ3pSkQUUGx/MSICig8+FWQDgqIAVEiIFICIZQklNBSIL2waaT3ZHeT3Z2dmff3x72bTJad3ZnZ6ft+npxn7n3POe95957ceec9555zRVUxDMMwjFjwpNoAwzAMI3MxJ2IYhmHEjDkRwzAMI2bMiRiGYRgxY07EMAzDiBlzIoZhGEbMmBMxDMMwYsaciGEYhhEz5kQMwzCMmDEnYhiGYcRMTqoNSDSdO3fW8vLyQ2T79u2juLg4JlkyiFe7seiJtE5L5ZrLj+ZaN5anqk/i2Xa0elLVJ+Hkdq9EVyfWfmmtvDXXZt68eTtUtUskZbPWiYjIWcBZPXv25K677jokr7q6mpKSkphkySBe7caiJ9I6LZVrLj+aa91Ynqo+iWfb0epJVZ+Ek9u9El2dWPultfLWXJvx48evi7iwqmZ1GjlypDZm5syZMcuSQbzajUVPpHVaKtdcfjTXurE8VX0Sz7aj1ZOqPgknt3slujqx9ktr5a25NsBcjfA7VjRLd/ENiUSuffrppw/Js19Xra9jkUjy9KSiT9bsDfDeRj+nH1ZPl/Z2r7SmToZGIvNUdVREhSP1NpmaLBJJTB2LRJKnJ9l9smjDHh12y+va96ap+vU/vqb76/wt6rF7JfZymR6J2NNZhmEcYPmWSi77+8e0K8jld2cO4fPdQa57eh4+fzDVphlpig1nRSlLBtkcoofLs+Gs1pdvbZ9UUsQfP6nBK8KkYwvoWuRhekU1z1QIY7p7ue6ofDwidq9EWceGszI82XBWYurYcFby9CSjT6ZMfVtH3/6mjvzDm1qxreqQOn97p0L73jRVb3p+gQaDQbtXoqyT7cNZWfuIr2EYkbFh937+NKcW9eYyZeJxDOhy6K/XiScPoLLGz4MzK2hXmMvxhdk5emHEhjkRw2jDbN5bw7cf+5gav/LcxDEM6lbaZLmff2MQlbX1PPruanYOzGX8+CQbaqQtNicSpSwZZPM4b7g8mxNpfflo+2T13gAPzK+jxq/cMEwZ2qP56x9U5bFFdXy4KcD1R+czqntOWN3JIpvvFZsTSZNkcyKJqWNzIsnTk4g+eW7ueh1482t6wp1v65KNeyO+/j5/QE++4zUdffuburfGF7WN8Sab75VMmROxR3wNow3hDwS57ZUl/OK5BYzs04FXbjiRIYe1i7h+rtfDlUPz2F5dx93TVyTQUiNTMCdiGG2EKp9yxROf8I/Za/nu2HL+efUYOhbnRa2nf5mXK47ry1MfrWPB+j0JsNTIJMyJGEYbYPHGvdz2YQ1z1+3mrguP4pazhpLjjf32//lpg+lams+kFxbhD9hCxLaMORHDyGJWbKni+n/P56wH3ycQhGe/fzwXjOzVar3tCnK59ayhLN1cyeQP1rbeUCNjsaezopQlg2x+4iRcnj2d1fryoeXWVQZ4ZVU9c7cGKPDC1/rkclJXH907xG8reFXl3vl1LN8V4OYRSp/Odq/EUs6ezkrzZE9nJaaOPZ2VPD2Rlp8xY4bOX7dLr548R/veNFWH/e51vXv6ct1VXdeinmieBAqVfbFznx7xm2l6zt3TIrIx3mTzvZIpT2fZYkPDyGACQWXeut1MX7KFl+bVsGP6B5QV5vLTUwdx1dhyygpzE9p+745F3HjqQO6ctpzpS7Zw2tDuCW3PSD/MiRhGhlFZW8+8tbt5Y+kW3ly6lR3VPvK8Ho7o6OEXpw/hm0f2oLQgsc4jlO+d2I+n31/JLS8tYezhnSnJt6+VtoT1tmGkMfWBIGv3BvjnR+v47Is9LNiwh1Xbq1GF4jwv44/oymlDuzNucBfmfTSbcWP6JN3GhrUjd3xcy60vL+HPFxyJiCTdDiM1ZJQTEZFi4GHAB8xS1X+l2CTDiAuVtfWs2b6PNTv2sXp7Nat3OMertldTWx8EFtOpOI+je7fnW0cdxog+7Rld3pGCXG+qTQfg8PZebhh/OPfPqGBIj3Z878R+qTbJSBIpdyIi8gRwJrBNVYeFyCcA9wFe4HFV/SNwHvC8qr4iIv8BzIkYaYuqUlXnZ/c+H7v2+di938e2yjo+qvDx+s6FbKmsZcveWrZU1rJnf/2Beh6BXh2K6N+lmOP6dyK3ahPfOe0EenUoTOtf+DeeOojlW6q4/dWlHN61hJMHdUm1SUYSSLkTASYDDwJPNQhExAs8BHwd2ADMEZGXgV7AIrdYILlmGtmKqhIIKr5AkLr6IHX+IHX+gPNZH6TWH2C/L8C8rX52f7qBGl+Q/T4/Nb4AVXV+qmrrqar1u8k53lNTz+59PvzBLz9CL0Cnkm30KCugV4ciRpd35LD2hfTvUsyALsX07lhEfs7BCGPWrG307liUxCsSGx6PcM/FR3P+Ix9w/b/n89L1J9Kvc3GqzTISTMqdiKq+KyLljcRjgApVXQ0gIlOAs3EcSi/gMxK8UPIfs9ewYk09Fd7VTdgcnzaUphWtWlPPSs+qpuuEaVubKLN6tY+lVESlZ9VqH4uDnx+S33CoetDmNWt8zK9feSBDQ3Qqyrp1Pj6uXX6wjvOPYFBZv76Od6uWoqiTr8qGjXW8vWcxQVWc710lGIRNm+t4ZdsCV65s3lLLcxvnEwgqAfcRw0BQ8Qed/EDQqecPBg/I/QHFHwweOK4PBN3kOI76QDDyPv10wSGnBbkeSgtyKc3PobQgh9KCXLqWFtChOJcORXl0LM478Nm+KJeu7QpYNv8jTv1qdu6lXpyfw2NXjOJbD77PtU/N5YUfnkC7JE7yG8kn6sWG4sTT3wH6q+rvRaQP0F1VP4nZCMeJTG0YzhKRC4AJqnqNe345cCxwE07UUgu8H25OREQmAhMBunXrNnLKlCmH5Eey2G3im/vwWawTMQ2DLAdHWxRBHLk4+SIN5dy8ELmq4pEGmeA5UDaIRzx4xBnmUQ2S4/EgglMe8AruuZO8B44Fr4DXc1DmFcHrgRwP5AjkeNxzcWR5HiHXC7keyA05zvcKgboa2pcWke91zvO8Tv1oSZet4CPJi3Vh7rKdAe6aW8uwzl5+cozzWt1EYIsNw8vTdrEh8AjOUNMy97wDMCdaPY10lgOLQ84vxJkHaTi/HHggFt2xLjasrPHpa2/O0MoaX5OpqrY+Lqm6iTTtzRlNyhvSvrqm0/46/4FU4/Pr9LdmaI3PHzbV1h9MdfUBrasP6Jtvz1CfP6A+f0DrQ5I/EFR/IKiBQFCDwaDOmDHjS9ewpWsczfUPJ7fFhrGXS8Riw+ZkT324VvveNFXvfG1Zs3a1BltsGF6erMWGsUQi81X1GBH5VFVHuLIFqnpUVIoO1VnOoZHI8cCtqnqaez7JdXh3RqHTtj1Jw19X4fJs25PWl0+nSKSByUvqmLXez8Qj8znhsPiPnmfzvZIpkUgs8wr17sS3AohIFyDe23jOAQaKSD8RyQMuAV6OcxuGYSSYy76SxxEdPTy+qI45W/ypNsdIBJGGLA0JZz7kZZxJ7juAFcCF0eoJ0fcMsBmod3Ve7crPAFYCq4CbY9Vve2clpo7tnZU8PZk6nNVAVW29nv/wbO0/6VWdumBTszZGSzbfK1k7nAUgIkcAX8OZ+3xbVZe1ypMlABvOSs8QPVyeDWe1vnw6Dmc1UONX7plXS8WeINcdlc+Y7vEZ2srmeyVThrOiiRg6Npci1ZPsZJFIYupYJJI8PZkeiTRQVVuvFzziRCSvLNjYYvlIyOZ7JesiERFZgzMPIkAfYLd73B74QlXTap8Di0TS89dVuDyLRFpfPp0jkQYOiUiOzGdMj9ZFJNl8r2RdJNKQgL8CZ4Scnw7cHa2eZCWLRBJTxyKR5OnJlkikgeqQiOTlz1oXkWTzvZIpkUgsT2eNVtXXQpzQNOCUGPQYhtEGKc7PYfJ3xzCyTwd+POVTnrTX62Y0sawTmQ68BzyNM7x1GXCyums60gUbzkrPED1cng1ntb58JgxnhVLnV/66sI5PtwU4rTyHiwfnRb2yPZvvlWwezuqIs7vup266D5tYjyvZHKKHy7PhrNaXz5ThrFD8gaDe8tJi7XvTVL3un3O1xuePqn423yuZMpwV9ayWqu4CfhJtPcMwjMZ4PcItZw2hV4dC7nhtGVsf+4jHrhhFp5L8VJtmREgsw1kz4cvbz6rqV+NlVDyw4az0DNHD5dlwVuvLZ9pwVmPmbPHz6MI6OhYIPxtZQLfilqdss/leyebhrJEhaSzwF+BP0epJVrLhrMTUseGs5OnJ5uGsxsxdu0uPvm26Hn3bdJ1dsb3F8tl8r2TKcFbUT2ep6ryQNFtVf4azTbthGEarGNm3Ay/8cCwdi/O4/O+f8Ph7qxt+vBppStROREQ6hqTOInIa0D0BthmG0Qbp17mYF380llO/0pXbX13Gjf/5jBp7uU/aEsucSOjKdT+wBvi9qr4ff/Nix+ZE0nOcN1yezYm0vnymz4k0JqjKq6vreeHzenqVerhhRD5diw793ZvN90o2z4kUNCHLj1ZPspLNiSSmjs2JJE9PW5oTadKG5Vt1+C2v65G3Ttd3VmxLSLvpeK9k7ZwI8EETsg9j0GMYhtEi4wZ35ZUbTqRHWQFX/uMT7nvrcwJBmydJFyJeJyIi3YGeQKGIjODga7XbAUUJsM0wDAOAvp2KeeGHJ3Dz/xZzz1sr+WTtTu65+OhUm2UQhRMBTgOuAnrhPNbbQBXw6zjaZBiG8SWK8nL4y0VHcfyATvzupcWccd97fPcID+NSbVgbJ5aJ9fNV9b8Jsidu2MR6ek4WhsuzifXWl8+2ifXm2FgV5KEFtWyuDnJm/zzOOTwXrye6fbdCScd7Jesm1oHL3M+fAz9rnCLVk+xkE+uJqWMT68nT09Yn1sOxv86vl93/uva9aape+MgHumnP/ph1peO9ko0T68XuZwlQ2kQyDMNIGoV5Xq4ens89Fx/F4k17mXDve0xbtDnVZrU5Ip4TUdW/uZ+3Jc4cwzCM6Dh3RC+O7t2BG6d8yg/+NZ+LRvXilrOGUpwfn/e4G80T9VUWkS7AtUB5aH1V/V78zDIMw4icfp2Lef4HJ3DfW5/z0KwKPlmzi3svGcHRvdun2rSsJ5Z1Ii8BZcBbwKshyTAMI2Xkej384rTBTLn2OOoDyvmPfMADb9uakkQTS7xXpKo3xd2SCBCR/sDNQJmqXpAKGwzDSG+O7d+J135yEr99cTF3v7mSWSu385eLjqJvp+KWKxtRE0skMlVEzoi2kog8ISLbRGRxI/kEEVkhIhUi8qvmdKjqalW9Otq2DcNoW5QV5nL/pSO49+KjWbm1itPve49/fbzOdgROALE4kZ/gOJIaEakUkSoRqYyg3mRgQqhARLzAQ8DpwBDgUhEZIiLDRWRqo9Q1BlsNw2jDnDOiJ9NvPJlj+nTg5v8t5qp/zGFrZW2qzcoqYnmfSKmqelS1UFXbueftIqj3LrCrkXgMUOFGGD5gCnC2qi5S1TMbpW3R2moYhnFY+0Ke+t4Yfn/2UD5es5Nv3PMuLy/YlGqzsoZYVqwf04R4L7BOVf0t1C0HpqrqMPf8AmCCql7jnl8OHKuq14ep3wm4A/g68Liq3hmm3ERgIkC3bt1GTpky5ZD8dF2FG+9203EVbrg8W7He+vJtacV6rO1u2Rfk0YV1rN4bZEx3L5cPyac0T9LyXsmUFeuxTKw/DBwDLHLPhwMLgE4icp2qvhGFrqb2KQjr1VR1J3BdS0pV9VER2QyclZOTMzIKewzDyGK6F3u4+dgCXl1Tz0sV9SzftZ8rh+Yz2ObcYyfSpe0NCWfIaWjI+RDgH0B/4LMW6pYDi0POjwemh5xPAiZFa1NzybY9SUwd2/YkeXps25PwtKbdpZv26un3vqt9b5qqF987TXdV1yWk7Wzf9iSW4azPVPXopmRN5TUqV86hw1k5wErga8BGYA7wbVVdEpVRTbdlGzCmYYgeLs+Gs1pf3oazoscfdN6e+NIqHyW5Hq4cmsfIbpEN0Nhwlkuk3qYhAf8BHgFOcdPDwLNAPjCnmXrPAJuBemADcLUrPwPHkawCbo7WnpaSRSKJqWORSPL0WCQSnni1++RLbx2ISn78zPyIohKLRGKPRAqBHwIn4sxpvO86klqchYjVUSlMEBaJWCQSbywSiU6WDOLZJwVFxbyyqp6pq+spzoXLh+Qzunv4qMQiEZdIvU2mJotEElPHIpHk6bFIJDyJ6JPFG/foGfc5UckPnp6r26tqW9W2RSKNEJGBwJ04E+oFIc6of1SKEoxFIhaJxBuLRKKTJYNE9Yk/qExzn+AqyIHvfCWf43p4EZGwdWK1sc1FIjjDV18DFgJ9gVuB26LVk6xkkUhi6lgkkjw9FomEJ9F9snJLpZ794Pva96apevXkT3Tznpqo27ZIpBEiMk9VR4rIIlUd7sreU9WTolKUYCwSsUgk3lgkEp0sGSSjT4KqvLHWz38/95HjgYsH53FKrxz27dtnkQjEFInMxtku5QXgeuBcYEW0epKVLBJJTB2LRJKnxyKR8CSzT9Zsr9aL//aB9r1pql7ytw91ytS346I70yORWDZgvBEoAn4MjAQuB66MQY9hGEbGUN65mH9fcxx3njecxRv38tvZNTz27mr8gWCqTUspUQ9nZQo2nGXDWfHGhrOikyWDVPXJ7togTyzcz6JdQr92Hr43PJ/epU3/JrfhrIPDWC83lyLVk+xkw1mJqWPDWcnTY8NZ4UnlvTJjxgx9ZcFGPeb3b+iASa/qn19frjU+f9S6M304K5oNGI8H1uOsPP+YpjdPNAzDaBOICGceeRhjB3Tm9leX8eDMCl5btJk7zxvOsf07pdq8pBHNnEh34NfAMOA+nO3Yd6jqO6r6TiKMMwzDSHc6FOdx90VH8c+rx+ALBLn40Y/49f8WUVlbn2rTkkJMcyIikg9cCvwZ+L2qPhBvw1qLzYnYnEi8sTmR6GTJIN3ulTq/8kKFjzfW+inLFy4fksfg4lqbE2lIOJssngc8h7Pj7m+BntHoSHayOZHE1LE5keTpsTmR8KTrvfLZF7v1tHve0b43TdXz/jJNt+ytCVs20+dEIh7OEpEngQ9wXkh1m6qOVtU/qOrG6HycYRhGdnNU7/a8csOJ/HLCYBZuD3Dq3e/w9EfrCAaz72nYaCbWLwf2AYOAH4fsISOAagTvWTcMw2gr5Ho9/HDc4XTc9wUvbyriNy8u5qXPNnLnecM5vGtpqs2LGxFHIqrqUdVSN7ULSaXmQAzDMJqme7GHf11zLH+64EhWbq3mjPve5963VlLnD6TatLgQy4p1wzAMIwpEhItG9eatn53CacO6c+9bn3PGfe/xyZpdqTat1diK9ShlySDdnjiJpZw9nRU/PfZ0Vngy9V5ZuN3Pk0t87KxVxnZXvj20mOJcabFeNPK0fDorE5M9nZWYOvZ0VvL02NNZ4cnke6W6tl7/8MoSLb9pqo66/U2dumCTBoPBFutl7NNZhmEYRvwozs/hN2cO4ZbjC+jWLp8f/Xs+1zw5l417alJtWlSYEzEMw0gh5WVeXvzhWH7zza/wwaqdfP0v7/D399cQzJCphmge8TUMwzASQI7XwzUn9ee0od353UuL+cPUpZS389B10F6G9SxLtXnNklGRiIicIyKPichLIvKNVNtjGIYRT3p3LOKJq0bz4LdHsKtW+daD73P71KXsq/On2rSwJM2JiMgTIrJNRBY3kk8QkRUiUiEiv2pOh6q+qKrXAlcBFyfQXMMwjJTQsDvwnScVcvHoPjz+/hq+cc+7zFy+LdWmNUkyI5HJwIRQgYh4gYeA04EhwKUiMkREhovI1Eapa0jV37j1DMMwspLiXOHO84bz3HXHU5Tn5buT5/Cjf81nT216vUkxaXMiqvquiJQ3Eo8BKlR1NYCITAHOVtU7gTMb6xBnr5U/AtNUdX5iLTYMw0g9o8s78uqPT+LRd1dx/4wKZhCkuv06vj2mDx5P6l/rlNTFhq4Tmaqqw9zzC4AJqnqNe345cKyqXh+m/o9x3uc+B/hMVf8aptxEYCJAt27dRk6ZMuWQfFtA1fo6ttgweXpssWF4svleaUq+ZZ/zWt6Ve4XD23u4amg+vdzX8raJxYZAObA45PxC4PGQ88uBB+LZpi02TEwdW2yYPD222DA82XyvhJPPmDFD/ztvvY5wX8v7x2nLtMbnT9liw1RHIscDt6rqae75JNex3RmHtmzbkzT8dRUuzyKR1pe3SCS5elIRiYTKq3zKf1b4eH+jny6FwkX9g4zunfxIJNWP+M4BBopIPxHJAy4BXk6xTYZhGGlPaZ5wzfB8bhpdgFfgoSXCXxfUsrcuyYsUIw1ZWpuAZ4DNQD2wAbjalZ8BrARWATfHu10bzkpMHRvOSp4eG84KTzbfK9HIa+v9+pPHpuvAX7+mw295Xf/98Tp9e8aMiOxrCtJ1OCuZ2HBWeobo4fJsOKv15W04K7l6Uj2c1ZS8kiKeXFLHit1BBrRTfjmmmPyc6J/gStuJ9VQki0QSU8cikeTpsUgkPNl8r8QqDwaD+uycL/TKB16PyL6mwCIRi0Ri1WORSHgsEolOlgyy+V6x94mkSbJIJDF1LBJJnh6LRMKTzfdKa+WtuTbY+0QMwzCMZJC1W8G7w1k/AipF5PNG2WXA3ghknYEdibGwWZqyJVl6Iq3TUrnm8iO9/k3JU9UnTdmSLD2p6pNwcrtXoqsTa7+0Vt6aPukbcclIQ5ZMTMCjkcrDyCIO6ZJhdzL0RFqnpXLN5Ud6/ZuSp6pPUtkvqeqTaPrK7pX490tr5cnqk2wfznolCnm4sqkgXrbEoifSOi2Vay4/mutv/ZK6Pgkntz6Jrk6s/RIveULJ2qez4oGIzNVIn1AwkoL1SXpi/ZJ+JKtPsj0SaS2PptoA40tkbZ+ISLmIqIhk4lxl1vZLBpOUPrFIxDBaiYisBboBAZxtfT4ArlPV9VHqKQfWALmqmr7vQzWMECwSMYz4cJaqlgA9gK3AAym2xzCSgjkRw4gjqloLPI/zumdEpExEnhKR7SKyTkR+IyIeN88rIneJyA4RWQ18M1SXW/fvIrJZRDaKyO3uK6UNI23IxLFXw0hbRKQIuBj4yBU9gPP8fn+gE/AGzm7WfweuxXkN9AhgH/DfRuqexIlqDgeKganAeuBvCf0jDCMKbE4kRkSkP3AzUKaqF6TanraKiBQDDwM+YJaq/isFNqzFWdjlB0qAbcBpwFJgPzBCVZe6Zb8PXKqq40RkBvCsuq95FpFvANOBXByH8wXQXlVr3PxLgYmqOj6Jf17U2L2RnojIOTjRblfgIVV9Ix562+Rwlog8ISLbRGRxI/kEEVkhIhUi8qvmdKjqalW9OrGWtk2i7J/zgOdV9VrgW0k39iDnqGp7IB+4HngH6AXkAetCyq0DerrHh+FEFqF5DfTFcSabRWSPiOzBiUC6Jsb85ommT+zeSB5R9suL7n1yFU60HBfapBMBJgMTQgXuWPNDwOk449mXisgQERkuIlMbpZTcyG2IyUTYPzhf1A1fxIEk2tgkqhpQ1RdcW47DeVordAuJPsBG93gz0LtRXgPrgTqgs6q2d1M7VR2aOOubZTKR94mRPCYTfb/8xs2PC23Siajqu8CuRuIxQIX7K8oHTAHOVtVFqnpmo7Qt6Ua3IaLpH5y3ZPZyy6T8/7M4nA10ABYDzwJ3iEipiPQFfgY0vJvgWeDHItJLRDoAB6JfVd2MM39yt4i0ExGPiAwQkVOS+gcdtCeaPjGSRDT94v7f/D9gmqrOj5cNKb/p0oieHDq0sIGDww5fQkQ6ichfgREiMinRxhlh++cF4HwReYTUbsfxiohUA5XAHcCVqroEuAFn0nw18D7wb+AJt85jOHMgC4D5OH9LKFfgDIctBXbjPPXVI7F/RlQ02Sd2b6SccPfKDcCpwAUicl28GrOnsw7S1Dskwz51oKo7gbh1hNEiTfaPqu4DvptsYxoZUd5M3m7gsjB5fuCnbmrgoZD8vcAP3JSOhOsTuzdSS7h+uR+4P96NWSRykA0cOj7dC9iUIluML2P9k35Yn6QnSe0XcyIHmQMMFJF+IpIHXAK8nGKbjINY/6Qf1ifpSVL7pU06ERF5BvgQGCwiG0Tkando4XqcMeplOM/vL0mlnW0V65/0w/okPUmHfrHFhoZhGEbMtMlIxDAMw4gPWf90VufOnbW8vPwQ2b59+yguLo5Jlgzi1W4seiKt01K55vKjudaN5anqk3i2Ha2eVPVJOLndK9HVibVfWitvzbWZN2/eDlXtElHhZLyDN54JZ3XmCqAC+FVL5UeOHKmNmTlzZsyyZBCvdmPRE2mdlso1lx/NtW4sT1WfxLPtaPWkqk/Cye1eia5OrP3SWnlrrg1RvJ89o4azbJsFwzCM9CLThrMOLOcHEJGGbRaWxruh7VV17KkLsq2qFnHX7kjIEp6GQ3GFB885UB45WEeaKXuIXgF/UKkPBA/UkZC2RZpaR2QYhpEaMs2JNLWc/9hENHTyn2ZSUx+AmW8nQn3LvDEtbJbjfA51MA3OqOHYIxAIBMidOd2VO7KDdRrOwRMi89XVUjRnZojMyfeIHDj2eoTq6hruXTIbj4DXI4gIXhE8HqfM3t21PLV2jlv+0DLbt9cydfsCt7yQ4xE2b6rjveql5Hgc/c6nhy/W+Vguqw7I13xRz5ZPvsDrEXK9HnK8Qo7HQ6734Hmu1+MmIc/rIS/HOW/4zM/xkOf14PGYQzaM1pJRj/iKyIXAaap6jXt+OTBGVW9oVG4iMBGgW7duI6dMmXKInurqakpKSpqVvbO+nv21deTn5wOH7n/S+JJpyIE2lrnlG19ldQsfUt49r/P5yMvLO9CONtITWrax/qAe0I7PV09Obm6IPERHiK5Qma++npycHDSkTtDND+rBsvV+Px5PDooS1EPLBNVxYHi8bh23jNuGPxAE8RDQhvJKIKgo4tR1U6LxCuTcCtQNAAAb+klEQVR6INcLuR5xjj2Q55UDn3leyPNAvleQYD0lBXnke53z/BwoyBEKvFCYI4ccF+Y4jrMpmvr/1xyRlm+pXHP54fIiuVeisTHexKvdWPQkul9aK2/NtRk/fvw8VR0VSdlMi0QiWs6vqo8CjwKMGjVKx40bd0j+rFmzaEk2Lky5ZBCvdmPRE2mdlso1lx/J9QcIBpUZs2Yx9sSTqQ8GCQSUd95/n2OPOx5/wBnyaxj68wcUfzCIz+981gecY+fTOa8PBKnzB/EFgtTVh34GqKsPUusPUlsfoLa+4TxAVX2A2rog+30Bqmv81AXrv/QjIhwl+TmUFebSrjCXdgXOccfiPKp3+hgxpC+dS/LoWJxHp+J8upTm07kkr8nhylT1STh5pLJkkM33SmvlyeqTTHMiB5bz47yT4RLg26k1yUgUDUNdhXleCnFeLd4+30OPssKU2DNr1ixOOeUU6vxBanwB9vn87PcFqKr1s6/OT3VDqvVTWVvP3pp6Kmv87mc963bu59P1e9hZXc/U1V+exsvL8dCjrIAeZQUc1r6Qw8oK6dWhkKrdAY7e76N9UV4K/mrDaJ6MciKq6heRhuX8XuAJtW0WjCQiIhTkeinI9dKhOLYv9RkzZ3LMsWPZuc/HzmofO6vr2FpZy+a9tWzaW8vmPTV8tGonW6vqCLhjind8/CZdSvM5vEsJA7uVMKhbKaPLOzKoW4k9bGGklIxyIgCq+hrwWqrtMIxY8YjQviiP9kV5DGhmOZc/EGTz3lpeePsDCrv14/Ot1Xy+rZoX5m+kus4PQOeSPI7r34kTBnTGsy+IqppTMZJKxjkRw2gr5Hg99O5YxFFdchh38oADclVlw+4aPly9kw9X7eSDVTuYunAzAPctnMG5I3ry3bH96FKanyrTjTaEORHDyDBEhN4di+jdsYiLRvVGVVmzYx+Tp33IJi3jkXdW8fj7azj/mF5MPLk//TqnZpsYo21gTsQwMhwRoX+XEr7aJ5dx40axZsc+Hn13Nf+dv4Epc75gwtDuXHfKgJYVGUYMmBMxjCyjX+di7jxvOD/9+kAmz17LPz9ax7TFWzipZw4nn6y2yNKIKxm1d5ZhGJHTtbSAX044gg8nfY3rThnAexv9/O7lxWTSAmMj/TEnYhhZTkl+Dr86/Qi+2S+Xpz/6gj9PX5Fqk4wswoazDKONcMGgXNp37cHDs1ZRWpDLD8bZPInResyJGEYbQUT4w9nDqK7z83+vL6ddYQ7fObZvqs0yMhxzIobRhvB4hLsuPIrqWj+/eXExJfk5nH10z1SbZWQwNidiGG2MXK+Hh75zDGPKO/LzZxfw9rKtqTbJyGDMiRhGG6Qg18vjV47iiB6l/PL5hQTtiS0jRsyJGEYbpbQgl++fPICd+3ys3htMtTlGhtLsnIiI9MLZbv0k4DCgBlgMvApMU1X7n2cYGcxJAzvjEVi4PZBqU4wMJWwkIiL/AJ4AfMD/AZcCPwTeAiYA74vIyckw0jCMxNC+KI9j+nRgkTkRI0aai0TuVtXFTcgXAy+ISB7QJzFmGYaRLMYf0ZU/T9/N9qo62/nXiJqwkUgYBxKa71PVivibZBhGMjllkPNSk3dWbk+xJUYm0uLEuogsEpGFjdJ7InKPiHSKt0EicquIbBSRz9x0RkjeJBGpEJEVInJavNs2jLbI0MPa0T5fmLViW6pNMTKQSBYbTgMCwL/d80vcz0pgMnBW/M3iHlW9K1QgIkPctofiTPK/JSKDVNUGcw2jFYgIwzt7eXfldvyBIDlee2jTiJxInMhYVR0bcr5IRGar6lgRuSxRhjXB2cAUVa0D1ohIBTAG+DCJNhhGVnJkFy/vbazj0/V7GF3eMdXmGBlEJD85SkTk2IYTERkDlLin/oRYBde7w2ZPiEgHV9YTWB9SZoMrMwyjlQzt5MXrsSEtI3qkpXcLiMhonEd9GxxHFXANsAT4pqo+G3WjIm8B3ZvIuhn4CNgBKPAHoIeqfk9EHgI+VNWnXR1/B15T1f82oX8iMBGgW7duI6dMmXJIfnV1NSUlJTHJkkG82o1FT6R1WirXXH4017qxPFV9Es+2o9WTrD55YImXGj/8fmxhs3XsXomuTqz90lp5a67N+PHj56nqqIgKq2pECSgD2kdaPh4JKAcWu8eTgEkhedOB41vSMXLkSG3MzJkzY5Ylg3i1G4ueSOu0VK65/GiudWN5qvoknm1HqydZffLwzArte9NU3bK3ptk6dq9EVyfWfmmtvDXXBpirEX5PR/J0Vjf3V/8UVd0jIkNE5OrofVtkiEiPkNNzcdalALwMXCIi+SLSDxgIfJIoOwyjrTH+CPdR3xX2qK8ROZHMiUzG+dV/mHu+ErgxUQYBf2p4rBgYD/wUQFWXAM8CS4HXgR+pPZllGHFjcLdSurcrYKbNixhREMnTWZ1V9VkRmQSgqn4RSdiXt6pe3kzeHcAdiWrbMNoyIsL4I7owdcFm6gNBcu1RXyMCIvlfss9dVKgAInIcsDehVhmGkRLGDe5KVZ2feet2p9oUI0OIxIn8DGc+YoCIzAaeAm5IqFWGYaSEsYd3JtcrNqRlREyLTkRV5wOnACcA3weGqurCRBtmGEbyKcnPYXR5R2Ytt8l1IzLCzomIyHlhsgaJCKr6QoJsMgwjhYwf3JU7XlvGpj01qTbFyACam1hv2BOrK04UMsM9Hw/MAsyJGEYWMm5wF+54bRmzVmw/8EimYYSjua3gv6uq38WZUB+iquer6vk4GyAahpGlHN61hJ7tC21exIiISCbWy1V1c8j5VmBQguwxDCPFNDzqO7tiB/5g89siGUYkTmSWiEwXkatE5Eqc96vPTLBdhmGkkBMGdGa/L8D6qmCqTTHSnBYXG6rq9SJyLtDwPvVHVfV/iTXLMIxUMrxnGQBr95oTMZqnuaezxN2IC9dpfMlxhJYxDCN76NWhkLLCXNZWmhMxmqe54ayZInKDiPQJFYpInoh8VUSeBK5MrHmGYaQCEWF4zzJzIkaLNOdEJuC8FvcZEdkkIktFZA3wOXApzitsJyfBRsMwUsCwnmVsqApS57d9To3whB3OUtVa4GHgYRHJBToDNaq6J1nGGYaROob1bEdAYeWWaob3Kku1OUaaEtE2napar6qbzYEYRtuhYXJ98Sbbb9UIj+31bBhGk/TpWERRDizaaE7ECI85EcMwmkRE6NvOw2JzIkYzRORERKSviJzqHheKSGlrGhWRC0VkiYgERWRUo7xJIlIhIitE5LQQ+QRXViEiv2pN+4ZhREZ5mZflm6vw+e0pLaNpInnH+rXA88DfXFEv4MVWtrsYOA94t1FbQ4BLcPbnmoAzqe8VES/wEHA6MAS41C1rGEYC6dvOgy8QZOXWqlSbYqQpkUQiPwLGApUAqvo5zs6+MaOqy1R1RRNZZwNTVLVOVdcAFcAYN1Wo6mpV9QFT3LKGYSSQ8nbOV8QSm1w3whCJE6lzv7gBEJEc3FflJoCewPqQ8w2uLJzcMIwE0rVIKM3Pscl1IyzS0q4lIvInYA9wBc5rcX8ILFXVm1uo9xbQvYmsm1X1JbfMLOAXqjrXPX8I+FBVn3bP/w68huPsTlPVa1z55cAYVW3yNb0iMhGYCNCtW7eRU6ZMOSS/urqakpKSmGTJIF7txqIn0jotlWsuP5pr3Vieqj6JZ9vR6klVnzTIH1zqxReA3x1fGLas3Suxl4v0/3608tZcm/Hjx89T1VEtlwRUtdmE8wV+LfAcztzItbjOp7UJ5+VWo0LOJwGTQs6nA8e7aXq4cs2lkSNHamNmzpwZsywZxKvdWPREWqelcs3lR3OtG8tT1SfxbDtaPanqkwb5H15ZogNvfk19/kDYsnavxF4u0v/70cpbc22AuRrh93gk71gPqupjqnohzq/7j91GEsHLwCUiki8i/YCBwCfAHGCgiPQTkTycyfeXE2SDYRghDO9Vhs8f5POt1ak2xUhDInk6a5aItBORjsBnwD9E5C+taVREzhWRDTgRxqsiMh1AVZcAzwJLgdeBH6lqQFX9wPU4kcky4Fm3rGEYCWZYw8p1mxcxmqDF94kAZapaKSLXAP9Q1VtEZGFrGtUwW8u7eXcAdzQhfw1nfsQwjCTSr1MxJfk5LN60l4vonWpzjDQjkqezckSkB3ARMDXB9hiGkWZ4PMKQw9rZE1pGk0TiRH6PM4xUoapzRKQ/znbwhmG0EYYdVsayzZX4A7Zy3TiUSF6P+xzOk1kN56uB8xNplGEY6cXwXu2onR2kYrtNrhuH0qITEZEC4GqcrUgKGuSq+r0E2mUYRhrRsC38og176ZJiW4z0IpLhrH/iLBo8DXgHZ+8s20jHMNoQ/TqXUJTnZcmmylSbYqQZkTiRw1X1t8A+VX0S+CYwPLFmGYaRTng9wpAeNrlufJlInEi9+7lHRIYBZUB5wiwyDCMtGdazjKWbKgkmbK2xkYlE4kQeFZEOwG9xVokvBf6UUKsMw0g7hvcso6Y+wOZqcyLGQSJ5Outx9/AdoH9izTEMI10Z3suZXF9bGUixJUY6EcnTWfk4j/SWh5ZX1d8nzizDMNKNAV1KKMj1sLbS1ooYB4lk25OXgL3APKAuseYYhpGuNEyur6u0J7SMg0TiRHqp6oSEW2IYRtozvGcZUz7ZQyCoeD2SanOMNCCSifUPRMQe6TUMg2E9y6gLwJodtnLdcAjrRERkkbtb74nAfBFZISILQ+SGYbQxRvbtAMDz8zam2BIjXWhuOOvMpFlhGEZG0L9LCWMPy+Hv76/molG96N8lNa8qNtKH5oaztgLnAv8PmABsVNV1DSkp1hmGkXZcNDiPghwvt76ylMS95NTIFJpzIk8Co4BFwOnA3fFqVEQuFJElIhIUkVEh8nIRqRGRz9z015C8ke5QWoWI3C8iNqtnGCmgLF/46dcH8e7K7byxdGuqzTFSTHNOZIiqXqaqfwMuAE6KY7uLgfOAd5vIW6WqR7vpuhD5IzjveB/oJntizDBSxBXH92Vwt1J+/8pSautt8WFbpjkn0rBnFu47zuOGqi5T1RWRlnffrNhOVT9UJ35+CjgnnjYZhhE5OV4Pt509lI17anhk1qpUm2OkkOacyFEiUummKuDIhmMRSeRqo34i8qmIvCMiDdFPT2BDSJkNrswwjBRxXP9OfOuow3jknVVs22+r2NsqkqiJMRF5C+c9JI25WVVfcsvMAn6hqnPd83ygRFV3ishI4EWcl2ENBu5U1VPdcicBv1TVs8K0PRFn6Itu3bqNnDJlyiH51dXVlJSUxCRLBvFqNxY9kdZpqVxz+dFc68byVPVJPNuOVk+q+iScPFS2uzbIpPdqGFim/HyM3SuxlIvm2kcjb821GT9+/DxVHdVySUBVU5aAWcColvKBHsDyEPmlwN8iaWPkyJHamJkzZ8YsSwbxajcWPZHWaalcc/nRXOvG8lT1STzbjlZPqvoknLyx7K+zKrTvTVN1xrKtzRuYALL5XmmtvDXXBpirEX6PR7JiPWmISBcR8brH/XEm0Fer6magSkSOc5/KugJnTy/DMFLMd8f2o3uxcNsrS6jz2yR7WyMlTkREzhWRDcDxwKsiMt3NOhlYKCILgOeB61R1l5v3A+BxoAJYBUxLstmGYTRBXo6Hy76Sz9qd+/nxM5+yvcr2aW1LRLIBY9xR1f8B/2tC/l/gv2HqzAWGJdg0wzBiYFhnL786/Qj+8sZKvnb3LCad8RUuHtUbj23SmPWk1XCWYRiZy3WnDGDajSfxlR7tmPTCIi5+9EMqtlWl2iwjwZgTMQwjbgzoUsKUicfxpwuOZOXWak6/7z3+8sYKW5CYxaRkOMswjOxFRLhoVG++ekRXbp+6lPtnVPDUR+v4+le6MWFYd8Ye3pmCXG+qzTTihDkRwzASQueSfO69ZAQXj+7Ds3PX8/qSLTw3bwPFeV7GHdGVCUO7M25wF0oLclNtqtEKzIkYhpFQjh/QieMHdMLnD/Lh6p28vngLby7dwqsLNyMC/TsXM6xnGcMOK2Noz3YM7VFGWZE5lkzBnIhhGEkhL8fDKYO6cMqgLtx+zjDmrdvN7IodLNlUyZw1u3jps00HyvZsX0h55yL6dCyiVwfns0/HInp3LKJDUS62iXf6YE7EMIyk4/UIY/p1ZEy/jgdkO6vrWLKpksWb9rJ8cxVf7NrPm0u3sqPad0jdPK+HLqX5dCnNR+pqeWP3IrqW5tOxOI/2RXl0KMqlfWEe7YtyaV+US0l+jjmdBGJOxDCMtKBTST4nD+rCyYO6HCLfV+dnw+4avti1n/W79rOtqo5tVbVsr6pjze4g6xZvYdc+XxitjsMqyc+htCCHkvwc2hXkOscFORTl5VCU52X7Zh/LZRXFeV4K83IozPVSmOehIMdLQZ6XwlwvBbleCnI95Od4yc/xEAjaC7nAnIhhGGlOcX4Og7uXMrh76ZfyZs2axbhx46gPBNmzv549+33sqaln9z7nc89+H5U1fqpq66mq81NV6xxvraqlYruf/b4A++ucz5dXLY/aNu9br5Gf4yEvx0Oe1/10jxvk1ZU1PLH6E/K8Qq7XQ67XQ45XyPU4n1s21/Fe9dJDZDkeYd1aHxXe1Xg9zrnX48Hrgc831rPn0414XLlHBK9HWLLNDyu24fU45xW7A4yLw/VvCXMihmFkPLkhQ1yxMGPmTI4bexL7fQH21fmprQ9SUx+gtj7gfPqcT58/SJ0/SJ0/wPLPV9GjVx/q6oP4AkF8fifVhRz7/EECCntr6vEHgtQHgtQHFJ8/iD8YxB9Qan1+PtryBfVBxR8IckiAs3JZ0wYv+qxp+fw5Bw7b5QnXnBvT5YgKcyKGYbR5PCLu0FYOnUsic0SzgusZN+6IlsvNmsW4cWNbyB934DwYVAKqzJj1DieMPZFAUKkPKEFV/EFl9gcfMnrMsQSCQfxBJRBUgkGYM3cuR404xikXUBYtDONo4ow5EcMwjDTC4xE8CPleaXINTdciD/06F39JvrPCy8i+HQ6c161PzoJO2/bEMAzDiBlzIoZhGEbMZPVwloicBewQkXWNssqAvRHIOgM7EmReczRlS7L0RFqnpXLN5Ud6/ZuSp6pPmrIlWXpS1Sfh5HavRFcn1n5prbw1fdI34pKRvgIxExPwaKTyMLKIXxGZDLuToSfSOi2Vay4/0uvflDxVfZLKfklVn0TTV3avxL9fWitPVp9k+3DWK1HIw5VNBfGyJRY9kdZpqVxz+dFcf+uX1PVJOLn1SXR1Yu2XeMkTirgey2gCEZmrqqNSbYdxEOuT9MT6Jf1IVp9keyTSWh5NtQHGl7A+SU+sX9KPpPSJRSKGYRhGzFgkYhiGYcSMORHDMAwjZsyJGIZhGDFjTiRGRKS/iPxdRJ5PtS1tGREpFpEnReQxEflOqu0x7N5IV0TkHPc+eUlEvhEvvW3SiYjIEyKyTUQWN5JPEJEVIlIhIr9qToeqrlbVqxNradskyv45D3heVa8FvpV0Y9sI0fSJ3RvJI8p+edG9T64CLo6XDW3SiQCTgQmhAhHxAg8BpwNDgEtFZIiIDBeRqY1S1+Sb3KaYTIT9A/QC1rvFAkm0sa0xmcj7xEgek4m+X37j5seFrN47Kxyq+q6IlDcSjwEqVHU1gIhMAc5W1TuBM5NrYdsmmv4BNuA4ks9ouz+KEk6UfbI0uda1XaLpFxFZBvwRmKaq8+Nlg910B+nJwV+04Hw59QxXWEQ6ichfgREiMinRxhlh++cF4HwReYT02o6jLdBkn9i9kXLC3Ss3AKcCF4jIdfFqrE1GImGQJmRhV2Kq6k4gbh1htEiT/aOq+4DvJtsYAwjfJ3ZvpJZw/XI/cH+8G7NI5CAbgN4h572ATSmyxfgy1j/ph/VJepLUfjEncpA5wEAR6SciecAlwMsptsk4iPVP+mF9kp4ktV/apBMRkWeAD4HBIrJBRK5WVT9wPTAdWAY8q6pLUmlnW8X6J/2wPklP0qFfbANGwzAMI2baZCRiGIZhxAdzIoZhGEbMmBMxDMMwYsaciGEYhhEz5kQMwzCMmDEnYhiGYcSMOREjIxGRgIh8FpLKU21TPBGRESLyeCt1TBaRC0LOLxWRm1tvHYjI9SJi280YtneWkbHUqOrR4TJFJMdddJWp/Bq4vbGwlX/XBOK3d9ITwGzgH3HSZ2QoFokYWYOIXCUiz4nIK8Abruz/icgcEVkoIreFlL3ZfWnPWyLyjIj8wpXPEpFR7nFnEVnrHntF5M8hur7vyse5dZ4XkeUi8i8RETdvtIh8ICILROQTESkVkfdE5OgQO2aLyJGN/o5S4EhVXeCe3yoij4rIG8BTIlLu6pnvphPcciIiD4rIUhF5FegaolOAo4H5InJKSAT3qdtec9fqCle2QET+CaCq+4G1IjImHn1nZC4WiRiZSqGIfOYer1HVc93j43G+gHeJ8wrQgTjvVxDgZRE5GdiHs5/QCJx7YD4wr4X2rgb2qupoEckHZrtf6rh6huJscjcbGCsinwD/AS5W1Tki0g6oAR7HebPcjSIyCMhX1YWN2hoFLG4kGwmcqKo1IlIEfF1Va0VkIPCMW+dcYDAwHOiG816PJ0JsXKCq6jrMH6nqbBEpAWqbuVY7gZuBsaq6Q0Q6htg0FzgJ+KSFa2dkMeZEjEwl3HDWm6q6yz3+hps+dc9LcL4oS4H/ub+mEZFINqf7BnBkyBxDmavLB3yiqhtcXZ8B5cBeYLOqzgFQ1Uo3/zngtyLy/4Dv4byZrjE9gO2NZC+rao17nAs86EY0AWCQKz8ZeEZVA8AmEZkRUn8CMM09ng38RUT+BbygqhtcJ9LUtToK5/XDO9y/Y1eIzm3AEU1fLqOtYE7EyDb2hRwLcKeq/i20gIjcSPh3xfg5OMxb0EjXDao6vZGucUBdiCiAc19JU22o6n4ReRPnDYAX4UQQjalp1DYc+nf9FNiK8wXvAWpDm2jqj8JxEOe7NvzRHe46A/hIRE4l/LX6cTM6C1xbjTaMzYkY2cx04HvukA0i0lNEugLvAueKSKE7H3BWSJ21OENHABc00vUDEcl1dQ0SkeJm2l4OHCYio93ypSLS8KPtcZwJ7jmNftk3sAw4vBndZThRThC4HPC68neBS9z5mx7AeLftMiDHfVkUIjJAVRep6v/hDEkdQfhr9TZwkYh0cuWhw1mD+PKwm9HGsEjEyFpU9Q0R+QrwoTvXXQ1cpqrzReQ/OO9lXwe8F1LtLuBZEbkcCB0OehxnmGq+O0m9HTinmbZ9InIx8ICIFOL8Yj8VqFbVeSJSSZgnm1R1uYiUiUipqlY1UeRh4L8iciEwk4NRyv+ArwKLgJXAO67868BbIfVvFJHxOFHTUpx3bteFuVZLROQO4B0RCeAMd13l6hkL3IbRprGt4I02j4jcivPlfleS2jsMmAUc4UYTTZX5KVClqq1aK+Lqehx4XFU/aq2uEJ0jgJ+p6uXx0mlkJjacZRhJRESuAD4Gbg7nQFwe4dC5lphR1Wvi6UBcOgO/jbNOIwOxSMQwDMOIGYtEDMMwjJgxJ2IYhmHEjDkRwzAMI2bMiRiGYRgxY07EMAzDiBlzIoZhGEbM/H8dbwpU3TDd4AAAAABJRU5ErkJggg==\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAfEAAAEWCAYAAAB2c65HAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvNQv5yAAAH5lJREFUeJzt3XuYXXV97/H3dxKSNE64xMBQEshAAYEkAmYECYgzxKaKiqWVFkGwKs6jtnir1tacoqDUqmDV9hyRoraWaGyxXlBaRGXwWIklw6XhapATMJAICISMCSFhvuePvQfHOJPZSWZm7V/yfj3Pfmbv316XD5ud+cxae+21IjORJEnlaak6gCRJ2jGWuCRJhbLEJUkqlCUuSVKhLHFJkgpliUuSVChLXBpFEfHBiLiy6hySdg+WuDSMiFgVERsjoi8ifh4RX4iI1qpzSdIAS1zatldlZivwAuCFwP+qOI8kPcsSlxqQmQ8C/wHMjYgDIuKbEfFYRNwbEW8ebr6IeFFE/CginoiI2yKic9xCS9rlWeJSAyLiQOBU4Bbgy8Bq4ADgNcDfRMTCIeaZCXwb+DAwHXgP8NWI2He8ckvatVni0rZ9PSKeAH4I3ABcDpwEvC8zn8rMW4ErgHOGmPd1wDWZeU1m9mfmdcByan8MSNJOm1h1AKnJ/X5mfnfgQUQcDzyWmesHTXM/0DHEvLOBMyLiVYPG9gCuH5OkknY7lri0fR4CpkfEtEFFfhDw4BDT/gz4l8wc9jNzSdoZ7k6XtkNm/gz4EfCRiJgSEc8H3gQsGWLyK4FXRcTvRcSE+vSdETFrPDNL2nVZ4tL2ey3QTm2r/GvAB+qfd/+aeuG/Gng/8Ai1LfP34r87SaMkMrPqDJIkaQe4RSBJUqEscUmSCmWJS5JUKEtckqRCFfU98RkzZmR7e3vVMYb1y1/+kuc85zlVxxhRKTmhnKyl5IRyspaSE5o/a29v76OZucOn++3t7d1v4sSJVwBzceNvPPUDt2/ZsuW8+fPnPzzUBEWVeHt7O8uXL686xrB6enro7OysOsaISskJ5WQtJSeUk7WUnND8WSPi/p2Zf+LEiVfsv//+R+67776Pt7S0+JWmcdLf3x+PPPLIUWvXrr0COG2oafyLSpI0krn77rvvkxb4+Gppacl99913HbU9IENPM455JEllarHAq1F/3YftaktckqRCWeKSpKY3YcKE+UccccRRhx122JyXv/zlh6xfv74FYOrUqcfu6DI//elPP3fVqlV7jF7Kob373e8+4IILLmjb2WmGYolLkpre5MmT++++++47V65ceccee+yRl1566Q4fbT/gyiuvnPHAAw+MeYmPJUtc2glLViyh/ZPt9K7ppf2T7SxZMdTFzKTdy2U3XTb9gEsPmNdyYcv8Ay49YN5lN102fTSXf9JJJ/Xde++9kwePrVu3ruWEE044/Kijjjry8MMPP+rKK6/cG+Cee+6ZdMghh8w588wzZx966KFzTjzxxMP6+vriC1/4wj6333771HPPPfeQI4444qi+vr4YvLzjjjvueW9605sO7OjoeN4hhxwy54Ybbpi6aNGi35k9e/bct7/97QcMTPfBD36w7bDDDptz2GGHzbnooov2Gxh/3/vet397e/vcBQsWHL5y5cpns95xxx2TX/ziFx82Z86cI+fPn/+8W265ZcrOvBaVfsUsIt4FnAcksAJ4Q2Y+VWUmqVFLViyh++puNmzeAG1w/7r76b66G4Cz551dcTqpGpfddNn0d33nXbOf2vJUC8CavjWT3vWdd80GeMsL3/LYzi5/8+bNXHvttXsuWrToycHjU6dO7f/2t7997/Tp0/vXrFkz8fjjjz/irLPOegLggQcemHLllVfet2DBgvtPPfXUQ774xS/u87a3ve2xz3zmM/tdcsklPzv55JM3DLWuSZMm9S9fvvyeD33oQ/udccYZh95000137bffflva29vnvf/97//5ypUrJ3/pS196bm9v712Zyfz5849cuHDh+v7+/vja1742fcWKFXdu3ryZY4455qhjjz12A8B55503+/LLL79/3rx5m77//e8/561vfetBy5Yt+8mOvh6VlXhEzATeDhyVmRsj4l+BM4F/qiqTtD0Wf29xrcAH2bB5A4u/t9gS127roh9cNHOgwAc8teWplot+cNHMnSnxTZs2tRxxxBFHARx//PHr3/GOdzw6+Pn+/v545zvfOWvZsmWtLS0tPPzww5NWr149EWDmzJmbFixYsBHg2GOP3bBq1arJv7mG33T66ac/AXD00UdvPPTQQzfOnj17M8CBBx646b777pvU09PTeuqppz6x55579gO84hWvePz666+f1t/fz6mnnvrEtGnT+gEWLVr0BNT2Ftxyyy2tZ5xxxu8MrOPpp5+O31xz46o+2ctE4LciYjMwldr1maUiPLDuge0al3YHa/vWTtqe8UYNfCY+3POf/exnp//iF7+YuGLFirsmT56cM2fOnLdx48YWgEmTJj379bgJEybkwPhIpkyZkgAtLS1Mnjz52WW0tLSwZcuWbV7KO+I3u/mZZ55h2rRpW7b137G9KvtMPDMfBC4BHgDWAOsy8ztV5ZG210F7HbRd49LuYP/W/Z/envHRsm7dugkzZszYPHny5Lz66qunPfTQQyP+0dDa2vrMunXrJuzoOk855ZS+a665Zu/169e3PPnkky3XXHPNPl1dXetPOeWUvm9/+9t79/X1xeOPP95y3XXX7Q0wffr0/lmzZj39+c9/fh+A/v5+brzxxt/a0fVDtbvT9wFeDRwMPAH8W0S8LjOv3Gq6bqAboK2tjZ6envGO2rC+vr6mzjeglJzQ3Fk/8bxPcP+6++nPfmZNnsUlh19CS7Qwe6/ZTZsZmvs1HayUnFBW1rF2wckXPDj4M3GAKROn9F9w8gUPjuV6zzvvvMde/vKXHzp37twj58yZs+Hggw8e8fiqc88999Hzzz9/9nvf+97+5cuX39Xa2rpdJ7Q56aSTNpx11lm/eMELXnAkwDnnnPPIiSeeuBHg9NNPf2zu3LlzZs6cuem4447rG5jny1/+8n1vfvObZ3/0ox/97S1btsTpp5/+2AknnLBxe/97B2xzd8BYiogzgJdl5pvqj88FXpSZbxtuno6OjvTc6TuvlJzQ/FmXrFjC4u8t5vy28/n7n/89Fy+8uOk/D2/213RAKTmh+bNGRG9mduzo/Lfddtuqo48++tGRp6y57KbLpl/0g4tmru1bO2n/1v2fvuDkCx4cjYPadle33XbbjKOPPrp9qOeq/Ez8AeBFETEV2AgsBJq3oaUhnD3vbM6edzY9PT2seu2qquNITeEtL3zLY5b2+KjyM/EfA1cBN1P7elkLcHlVeSRJKk2lR6dn5geAD1SZQZI0ov7+/v7wIijjr7+/P6hdV3xInrFNkjSS2x955JG96oWicVK/nvhewO3DTVP198QlSU1uy5Yt561du/aKtWvXzsWNv/HUD9y+ZcuW84abwBKXJG3T/PnzHwZOqzqHfpN/UUmSVChLXJKkQlnikiQVyhKXJKlQlrgkSYWyxCVJKpQlLklSoSxxSZIKZYlLklQoS1ySpEJZ4pIkFcoSlySpUJa4JEmFqrTEI2LviLgqIu6OiLsi4oQq80iSVJKqL0X6KeA/M/M1ETEJmFpxHkmSilFZiUfEnsDJwJ8AZObTwNNV5ZEkqTSRmdWsOOIY4HLgTuBooBd4R2b+cqvpuoFugLa2tvlLly4d76gN6+vro7W1teoYIyolJ5STtZScUE7WUnJC82ft6urqzcyOqnNo9FVZ4h3AMuDEzPxxRHwKeDIz/3q4eTo6OnL58uXjlnF79fT00NnZWXWMEZWSE8rJWkpOKCdrKTmh+bNGhCW+i6rywLbVwOrM/HH98VXACyrMI0lSUSor8cxcC/wsIp5XH1pIbde6JElqQNVHp58PLKkfmX4f8IaK80iSVIxKSzwzbwX8nEaSpB3gGdskSSqUJS5JUqEscUmSCmWJS5JUKEtckqRCWeKSJBXKEpckqVCWuCRJhbLEJUkqlCUuSVKhLHFJkgpliUuSVChLXJKkQlnikiQVyhKXJKlQlZd4REyIiFsi4ltVZ5EkqSSVlzjwDuCuqkNIklSaSks8ImYBrwCuqDKHJEklisysbuURVwEfAaYB78nMVw4xTTfQDdDW1jZ/6dKl4xtyO/T19dHa2lp1jBGVkhPKyVpKTignayk5ofmzdnV19WZmR9U5NPomVrXiiHgl8HBm9kZE53DTZeblwOUAHR0d2dk57KSV6+npoZnzDSglJ5STtZScUE7WUnJCWVm1a6lyd/qJwGkRsQpYCpwSEVdWmEeSpKJUVuKZ+VeZOSsz24Ezge9n5uuqyiNJUmma4eh0SZK0Ayr7THywzOwBeiqOIUlSUdwSlySpUJa4JEmFssQlSSqUJS5JUqEscUmSCmWJS5JUKEtckqRCWeKSJBXKEpckqVCWuCRJhbLEJUkqlCUuSVKhLHFJkgpliUuSVChLXJKkQlnikiQVqrISj4gDI+L6iLgrIu6IiHdUlUWSpBJNrHDdW4A/z8ybI2Ia0BsR12XmnRVmkiSpGJVtiWfmmsy8uX5/PXAXMLOqPJIklSYys+oMREQ78ANgbmY+udVz3UA3QFtb2/ylS5eOe75G9fX10draWnWMEZWSE8rJWkpOKCdrKTmh+bN2dXX1ZmZH1Tk0+hoq8YiYkJnPjEmAiFbgBuDizPz3bU3b0dGRy5cvH4sYo6Knp4fOzs6qY4yolJxQTtZSckI5WUvJCc2fNSIs8V1Uo7vT742Ij0fEUaO58ojYA/gqsGSkApckSb+u0RJ/PvAT4IqIWBYR3RGx586sOCIC+BxwV2Z+YmeWJUnS7qihEs/M9Zn5j5m5APgL4APAmoj454g4dAfXfSJwDnBKRNxav526g8uSJGm309BXzCJiAvAK4A1AO3ApsAR4MXANcPj2rjgzfwjE9s4nSZJqGv2e+ErgeuDjmfmjQeNXRcTJox9LkiSNZMQSr2+F/1NmXjTU85n59lFPJUmSRjTiZ+L1r5Z1jUMWSZK0HRrdnf6jiPgH4CvALwcGB864JkmSxl+jJb6g/nPwLvUEThndOJIkqVENlXhmujtdkqQm0/BVzCLiFcAcYMrA2HAHu0mSpLHX0MleIuIy4I+B86l9t/sMYPYY5pIkSSNo9LSrCzLzXODxzLwQOAE4cOxiSZKkkTRa4hvrPzdExAHAZuDgsYkkSZIa0ehn4t+KiL2BjwM3Uzsy/YoxSyVJkkbU6NHpH6rf/WpEfAuYkpnrxi6WJEkayfYcnb6A2sVPJtYfk5lfHKNckiRpBI1exexfgN8BbgWeqQ8nYIlLklSRRrfEO4CjMjPHMowkSWpco0en3w7sP9orj4iXRcQ9EXFvRPzlaC9fkqRdWaNb4jOAOyPiv4FNA4OZedqOrrh+idP/DfwusBq4KSK+mZl37ugyJUnanTRa4h8cg3UfB9ybmfcBRMRS4NWAJS5JUgOiqo+5I+I1wMsy87z643OA4zPzz7aarhvoBmhra5u/dOnScc/aqL6+PlpbW6uOMaJSckI5WUvJCeVkLSUnNH/Wrq6u3szsqDqHRt82t8Qj4oeZeVJErKd2NPqzTwGZmXvuxLpjiLHf+IsiMy8HLgfo6OjIzs7OnVjl2Orp6aGZ8w0oJSeUk7WUnFBO1lJyQllZtWvZZoln5kn1n9PGYN2r+fXzr88CHhqD9UiStEtq9Hvi04cYXp+Zm3di3TcBh0XEwcCDwJnAWTuxPEmSdiuNHth2M7Wt5sep7QbfG1gTEQ8Db87M3u1dcWZuiYg/A64FJgCfz8w7tnc5kiTtrhot8f8EvpaZ1wJExCLgZcC/Av8HOH5HVp6Z1wDX7Mi8kiTt7ho92UvHQIEDZOZ3gJMzcxkweUySSZKkbWp0S/yxiHgfMPD9rj8GHq+fsKV/TJJJkqRtanRL/CxqR49/HfgGcFB9bALwR2MTTZIkbUuj1xN/FDh/mKfvHb04kiSpUY1+xWxf4C+AOcCUgfHMPGWMckmSpBE0ujt9CXA3cDBwIbCK2ve8JUlSRRot8edm5ueAzZl5Q2a+EXjRGOaSJEkjaPTo9IEzs62JiFdQOz3qrLGJJEmSGtFoiX84IvYC/hz4e2BP4F1jlkqSJI2o0aPTv1W/uw7oGrs4kiSpUY0enX4wta+YtQ+eJzNPG5tYkiRpJI3uTv868DngajxDmyRJTaHREn8qMz89pkkkSdJ2abTEPxURHwC+A2waGMzMm8cklSRJGlGjJT4POAc4hV/tTs/6Y0mSVIFGS/x04JDMfHo0VhoRHwdeBTwN/BR4Q2Y+MRrLliRpd9HoGdtuA/YexfVeB8zNzOcDPwH+ahSXLUnSbqHRLfE24O6IuIlf/0x8h75ilpnfGfRwGfCaHVmOJEm7s8jMkSeKeMlQ45l5w04HiLga+EpmXjnM891AN0BbW9v8pUuX7uwqx0xfXx+tra1VxxhRKTmhnKyl5IRyspaSE5o/a1dXV29mdlSdQ6OvoRLfoQVHfBfYf4inFmfmN+rTLAY6gD/IBoJ0dHTk8uXLRzfoKOrp6aGzs7PqGCMqJSeUk7WUnFBO1lJyQvNnjQhLfBe1zd3pEbGe2lHov/EUkJm553DzZuZLR1j264FXAgsbKXBJkvTrtlnimTltLFYaES8D3ge8JDM3jMU6JEna1TV6dPpo+wdgGnBdRNwaEZdVlEOSpGI1enT6qMrMQ6tYryRJu5KqtsQlSdJOssQlSSqUJS5JUqEscUmSCmWJS5JUKEtckqRCWeKSJBXKEpckqVCWuCRJhbLEJUkqlCUuSVKhLHFJkgpliUuSVChLXJKkQlnikiQVyhKXJKlQlZZ4RLwnIjIiZlSZQ5KkElVW4hFxIPC7wANVZZAkqWRVbon/HfAXQFaYQZKkYkXm+HdoRJwGLMzMd0TEKqAjMx8dZtpuoBugra1t/tKlS8cv6Hbq6+ujtbW16hgjKiUnlJO1lJxQTtZSckLzZ+3q6urNzI6qc2j0jVmJR8R3gf2HeGox8H5gUWauG6nEB+vo6Mjly5ePbtBR1NPTQ2dnZ9UxRlRKTignayk5oZyspeSE5s8aEZb4LmriWC04M1861HhEzAMOBm6LCIBZwM0RcVxmrh2rPJIk7WrGrMSHk5krgP0GHm/PlrgkSfoVvycuSVKhxn1LfGuZ2V51BkmSSuSWuCRJhbLEJUkqlCUuSVKhLHFJkgpliUuSVChLXJKkQlnikiQVyhKXJKlQlrgkSYWyxCVJKpQlLklSoSxxSZIKZYlLklQoS1ySpEJZ4pIkFaqyEo+I8yPinoi4IyI+VlUOSZJKNbGKlUZEF/Bq4PmZuSki9qsihyRJJatqS/ytwN9m5iaAzHy4ohySJBUrMnP8VxpxK/AN4GXAU8B7MvOmYabtBroB2tra5i9dunTccm6vvr4+Wltbq44xolJyQjlZS8kJ5WQtJSc0f9aurq7ezOyoOodG35jtTo+I7wL7D/HU4vp69wFeBLwQ+NeIOCSH+IsiMy8HLgfo6OjIzs7OsYq803p6emjmfANKyQnlZC0lJ5STtZScUFZW7VrGrMQz86XDPRcRbwX+vV7a/x0R/cAM4JGxyiNJ0q6mqs/Evw6cAhARhwOTgEcryiJJUpEqOTod+Dzw+Yi4HXgaeP1Qu9IlSdLwKinxzHwaeF0V65YkaVfhGdskSSqUJS5JUqEscUmSCmWJS5JUKEtckqRCWeKSJBXKEpckqVCWuCRJhbLEJUkqlCUuSVKhLHFJkgpliUuSVChLXJKkQlnikiQVyhKXJKlQlZR4RBwTEcsi4taIWB4Rx1WRQ5KkklW1Jf4x4MLMPAa4oP5YkiRth6pKPIE96/f3Ah6qKIe0U5asWEL7J9vpXdNL+yfbWbJiSdWRJO1GJla03ncC10bEJdT+kFhQUQ5phy1ZsYTuq7vZsHkDtMH96+6n++puAM6ed3bF6STtDiIzx2bBEd8F9h/iqcXAQuCGzPxqRPwR0J2ZLx1mOd1AN0BbW9v8pUuXjkne0dDX10dra2vVMUZUSk5o7qwrHl7B0888DcCsybNYvWk1AJMmTGLefvOqjLZNzfyaDlZKTmj+rF1dXb2Z2VF1Do2+MSvxba40Yh2wd2ZmRASwLjP3HGm+jo6OXL58+dgH3EE9PT10dnZWHWNEpeSE5s7acmELSe3fzyWHX8J7fvIeAIKg/wP9VUbbpmZ+TQcrJSc0f9aIsMR3UVV9Jv4Q8JL6/VOAlRXlkHbYQXsdtF3jkjTaqirxNwOXRsRtwN9Q310uleTihRczdY+pvzY2dY+pXLzw4ooSSdrdVHJgW2b+EJhfxbql0TJw8Nri7y0GYPZes7l44cUe1CZp3FR1dLq0Szh73tmcPe9senp6WPXaVVXHkbSb8bSrkiQVyhKXJKlQlrgkSYWyxCVJKpQlLklSoSo5Y9uOiohHgPurzrENM4BHqw7RgFJyQjlZS8kJ5WQtJSc0f9bZmblv1SE0+ooq8WYXEctLOLVhKTmhnKyl5IRyspaSE8rKql2Lu9MlSSqUJS5JUqEs8dF1edUBGlRKTignayk5oZyspeSEsrJqF+Jn4pIkFcotcUmSCmWJS5JUKEt8CBExPSKui4iV9Z/7DDPdxyLijoi4KyI+HRFRH58UEZdHxE8i4u6I+MP6+Lsj4s6I+J+I+F5EzB60rGci4tb67ZtNkHVyRHwlIu6NiB9HRPugZf1VffyeiPi9ccrZU1/fwGu0X3387waN/SQinhi0rKpe0+Gy/klEPDJo/LxBy3p9fX0rI+L1FedsxvfpcFmb6n066PlvRsTtgx5/ZVD2VRFxa328PSI2DnruskZySs/KTG9b3YCPAX9Zv/+XwEeHmGYB8F/AhPrtRqCz/tyFwIfr91uAGfX7XcDU+v23Al8ZtLy+Jsv6NuCy+v0zB7ICRwG3AZOBg4GfAhPGIWcP0DHCOs4HPt8Er+mQWYE/Af5hiPHpwH31n/vU7+9TYc5mfJ8Ol7Wp3qf15/8A+BJw+zDruBS4oH6/fbjpvHlr5OaW+NBeDfxz/f4/A78/xDQJTAEmUftFsQfw8/pzbwQ+ApCZ/Zn5aP3+9Zm5oT7NMmBWs2bdarlXAQvrWxuvBpZm5qbM/H/AvcBx45CzEa8Fvrwd0w9nPLIO9nvAdZn5WGY+DlwHvKyqnE36Pm1kuZW/TyOiFXg38OGhFl7P9keMzvtUssSH0ZaZawDqP/fbeoLMvBG4HlhTv12bmXdFxN71ST4UETdHxL9FRNsQ63gT8B+DHk+JiOURsSwihvrFMd5ZZwI/q8+/BVgHPHfweN3q+tiY5Rw0yRfquxz/eojdl7OpbXF9f9DwuL+mDWT9w/pu6qsi4sD6WFO+pnWVv09HyNps79MPUdvS3rD1fHUvBn6emSsHjR0cEbdExA0R8eIGMkrPmlh1gKpExHeB/Yd4anGD8x8KHMmvtlKui4iTgTvrY/+Vme+OiHcDlwDnDJr3dUAH8JJBizwoMx+KiEOA70fEisz8aYVZh/qFntsYH7OcmfkD4OzMfDAipgFfrWf84qDZzwSuysxnBo2N+2s6QtargS9n5qaIeAu1Lb1TaNLXtFnepyNkbZr3KfAkcGhmvmvwZ/Nb2Xpv0Rpqr+kvImI+8PWImJOZTzaSRdptSzwzXzrccxHx84j47cxcExG/DTw8xGSnA8sys68+z38ALwL+L7W/wr9Wn+7fqG3NDCz7pdR+WbwkMzcNyvNQ/ed9EdEDHEvtc7yqsq4GDgRWR8REYC/gsUHjA2YBA9nHKucPMvPB+jrWR8SXqO0a3brE/3Twwip6TYfNmpm/GDT/PwIfrd9fDXQOem4Wtc+AK3tNm+x9uq2szfQ+XQ/Mj4hV1H637hcRPZnZWZ9uIrXPy+cPLKj+2m6q3++NiJ8ChwPLh8soDebu9KF9Exg4Qvj1wDeGmOYB4CURMTEi9qC2tXJXZia1La7O+nQLqW3xEhHHAp8FTsvMZ385RMQ+ETG5fn8GcOLAPFVl3Wq5rwG+X5/+m8CZUTsq+GDgMOC/xzJn/fEMgPr4K4HBR/4+j9oBYTcOGqvkNd1W1nopDDgNGNgFey2wqJ55H2BRfayqnE31Ph3h/3/TvE8z8zOZeUBmtgMnAT8ZKPC6lwJ3Z+bqgYGI2DciJtTvH1LPeV8DOaWaHOcj6Uq4UftM7XvAyvrP6fXxDuCK+v0J1H7R3UXtF9knBs0/G/gB8D/1+Q+qj3+X2gEwt9Zv36yPLwBWUDuadgXwpibIOoXalvm91H75HTJonsXUtr7uAV4+1jmB5wC99Yx3AJ9i0JHGwAeBv91qfZW8ptvKSu0Awjvqma4Hjhi0zjfWX+t7gTdUnLOp3qcjZG2a9+lWy2lnq6POgX8C3rLV2B8Oek/cDLxqLH6nedt1b552VZKkQrk7XZKkQlnikiQVyhKXJKlQlrgkSYWyxCVJKpQlLm1D/OqqXbdHxNXxq1PV7siyVg1831mSRoMlLm3bxsw8JjPnUjsT2J+ONIMkjRdLXGrcjQy6iEZEvDciboraBU0uHDT+9Yjojdr1prsrSSppt2CJSw2onxpzIbXTchIRi6idIvM44Bhq58w+uT75GzNzPrWzfL09Ip5bQWRJuwFLXNq234qIW4FfANOpXesbauc3XwTcQu10mUdQK3WoFfdt1K7FfeCgcUkaVZa4tG0bM/MYaueYn8SvPhMP4CP1z8uPycxDM/NzEdFJ7UIXJ2Tm0dRKfkoVwSXt+ixxqQGZuQ54O/Ce+pWrrgXeGBGtABExMyL2o3YpzMczc0NEHEHtEpWSNCZ22+uJS9srM2+p7yY/MzP/JSKOBG6MCIA+4HXAfwJviYj/oXb1rGWVBZa0y/MqZpIkFcrd6ZIkFcoSlySpUJa4JEmFssQlSSqUJS5JUqEscUmSCmWJS5JUqP8PLAA8NI0f8Z0AAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZAAAAEWCAYAAABIVsEJAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvNQv5yAAAHU1JREFUeJzt3XucHGWd7/HPdyYhoAG5JNwSICHEo4FlQecgZncFwUtwdwH3iAZEUdEcFdSzq+7ievag7nJWRT1HRdQs8AJURLxBZFFA5OZCMBEQSRAZIoaRrAS5JYCQZH77Rz0dmkknqamerurq+b5fr3l191PV3b+ZqZlvP09VPaWIwMzMbLT6qi7AzMzqyQFiZmaFOEDMzKwQB4iZmRXiADEzs0IcIGZmVogDxKwgSedL+peq6zCrigPEbAsk3SfpKUlrJf1nCo3JVddl1g0cIGZb99cRMRk4CDgY+EjF9Zh1BQeIWU4R8Z/AlWRBsglJ75I0KOlhSYsk7ZnaZ0gKSROa1r1O0jtHPPcuSWskLZf0ktQekvZrWm/jsJmkKZIul/Roes8bJflv2krjjc0sJ0nTgaOAwRbLjgD+FXgjsAfwW+DinK97HPAx4K3ADsDRwB9yPPWDwBAwFdgN+EfAcxNZaSZsfRWzce9SSQFMBn4CnN5inTcD50XErQCSPgI8ImlGjtd/J/DpiFiSHm8SUJuxjiys9omIQeDGnM8zGxPugZht3bERsT1wOPAiYEqLdfYk63UAEBFryXoR03K8/l7AvQXqOpMsbK6StELSaQVew6wwB4hZThFxPXA+8JkWix8A9mk8kPR8YBfgd8ATqfl5Tevv3nT/fmDWZt72yc09LyLWRMQHI2Jf4K+Bv5N0ZK5vxmwMOEDMRuf/A6+WNHJH+kXA2yUdJGkS8H+BWyLivohYTRYkJ0rql/QOnhsY5wAfkvRSZfaT1Aij24ET0vPmAYc1niTpr9K6Ah4HNqQvs1I4QMxGIYXBhcA/jWi/JrV9F1hFFhDzm1Z5F/BhsmGt/YGbmp77beAMshBaA1wK7JwWf4Csd/Eo2X6WS5teczbwY2AtcDNwdkRc1/53aZaPfEEpMzMrwj0QMzMrxAFiZmaFOEDMzKwQB4iZmRXS02eiT5kyJWbMmFF1GWZmtfLzn//8oYiYurX1ejpAZsyYwdKlS6suw8ysViT9dutreQjLzMwKcoCYmVkhDhAzMyvEAWJmZoU4QMzMrBAHiJmZFeIAMTOzQhwgZiV7Yt0TnHXbWdz50J1Vl2LWFgeIWcmeWv8UX73jqyz/w/KqSzFriwPErGRCAAzHcMWVmLXHAWJWsj5lf3YOEKs7B4hZyRoBEvhqoFZvDhCzkkkewrLe4AAxK5n3gVivcICYlawxhGVWd96SzUrmHoj1CgeIWcl8FJb1CgeIWcl8FJb1CgeIWck8hGW9wgFiVrLGYbwR7oFYvTlAzEq2cR8I7oFYvTlAzErWGMJyD8TqzgFiVrLGEJZZ3TlAzMysEAeImZkV4gAxM7NCHCBmZlZIZQEiaS9J10q6S9IySR9I7TtLulrSPel2p9QuSV+QNCjpDkkvqap2MzOrtgeyHvhgRLwYOBQ4RdIc4DTgmoiYDVyTHgMcBcxOXwuAL5dfspmZNVQWIBGxKiJuTffXAHcB04BjgAvSahcAx6b7xwAXRmYxsKOkPUou28zMkq7YByJpBnAwcAuwW0SsgixkgF3TatOA+5ueNpTaRr7WAklLJS1dvXp1J8s2MxvXKg8QSZOB7wL/KyIe39KqLdo2OZU3IhZGxEBEDEydOnWsyjQzsxEqDRBJE8nC4xsR8b3U/PvG0FS6fTC1DwF7NT19OvBAWbWamdlzVXkUloBzgbsi4nNNixYBJ6X7JwGXNbW/NR2NdSjwWGOoy6yOfD0Qq7sJFb73nwFvAX4p6fbU9o/AJ4FLJJ0MrASOS8uuAF4HDAJPAm8vt1yzsdOnPl8PxGqvsgCJiJ/Ser8GwJEt1g/glI4WZVYSIc/Ga7VX+U50s/FIkoewrPYcIGYVEPIQltWeA8SsAn3qcw/Eas8BYlaBPvV5H4jVngPErCIewrK6c4CYVcBDWNYLHCBmFfBhvNYLHCBmFRiOYfrkPz+rN2/BZhUIwgFitect2KwCwzGMNjsRg1k9OEDMKhARZPOJmtWXA8SsAsN4H4jVn7dgswpEhIewrPYcIGYV8FFY1gu8BZuVLCIIvA/E6s8BYlayxhQm7oFY3XkLNitZI0AmqMoLgpq1zwFiVrL1sR5wD8Tqz1uwWckaPZB+9VdciVl7HCBmJdsQGwD3QKz+vAWblWx4OPVA+twDsXpzgJiVrLEPxENYVncOELOK+Ex0qzsHiFnJGvs+hvElba3eHCBmJetLf3a+JrrVnQPErGR9fQ4Q6w0OELOSuQdivcIBYlayjftAHCBWcw4Qs5I1AqRxQqFZXTlAzErmM9CtV1S6JUs6T9KDku5sattZ0tWS7km3O6V2SfqCpEFJd0h6SXWVmxXXOP/DQ1hWd1V/FDofmDei7TTgmoiYDVyTHgMcBcxOXwuAL5dUo9nYSucPRkS1dZi1qdIAiYgbgIdHNB8DXJDuXwAc29R+YWQWAztK2qOcSs3GTuMorMABYvVWdQ+kld0iYhVAut01tU8D7m9abyi1PYekBZKWSlq6evXqjhdrNlqNS9m6B2J1140BsjmtJg7a5C8wIhZGxEBEDEydOrWEssxGp7EPxD0Qq7tuDJDfN4am0u2DqX0I2KtpvenAAyXXZta2Rg/ErO66MUAWASel+ycBlzW1vzUdjXUo8FhjqMvMzMo3oco3l/RN4HBgiqQh4HTgk8Alkk4GVgLHpdWvAF4HDAJPAm8vvWAzM9uo0gCJiOM3s+jIFusGcEpnKzIzs7y6cQjLzMxqwAFiZmaFOEDMzKwQB4iZmRXiADEr2YbhbBp3nw9idecAMSvZM8PPALBN3zYVV2LWHgeIWcme2ZACpN8BYvXmADEr2brhdQBM6Kv0NCyztjlAzErWuCKhZ+O1ussVIJL6O12I2XjR6HmsH15fcSVm7cnbAxmUdKakOR2txmwcmKAUIOEAsXrLGyAHAr8GzpG0OF20aYcO1mXWsxo9kMbhvGZ1lStAImJNRPxbRMwF/p5s1txVki6QtF9HKzTrMY19IBvCAWL1lnsfiKSjJX0f+DzwWWBf4Adk06ybWU79aZficAxXXIlZe/IeR3gPcC1wZkTc1NT+HUmvGPuyzHqXz0C3XrHVAElHYJ0fEZ9otTwi3j/mVZmZWdfb6hBWRGwAXllCLWZmViN5h7BuknQW8C3giUZjRNzakarMxoHAJxJaveUNkLnptnkYK4AjxrYcs/Fh+2225+GnHq66DLO25AqQiPAQltkYmrnDTO57/L6qyzBrS+7Z3CT9JbA/sG2jbXM71s1sy9auW8uU7aZUXYZZW/KeB/IV4E3A+wABxwH7dLAus5724JMPssfkPaouw6wteacymRsRbwUeiYiPAy8H9upcWWa97ekNT/t6IFZ7eQPkqXT7pKQ9gXXAzM6UZNb71g2v8xUJrfby7gO5XNKOwJnArWRHYJ3TsarMeljjOiCNObHM6irvUVj/nO5+V9LlwLYR8VjnyjLrXRtiAxM0gbXr1lZdillbRnMU1lxgRuM5koiICztUl1lP+sr193Lg9Bew/Tbbs+aZNQDcdO9D3DH0GO8+bFbF1ZmNTt6jsL4GfAb4c+C/p6+BDtZl1pMOnP4CTr3oNh55+hF22XYXbrr3IU696DYOnP6CqkszG7W8PZABYE74Is5mbZk7awqfe9Mc3ncznHvnuVy4cn/OOuFg5s7yOSFWP3n34t0J7N7JQvKSNE/S3ZIGJZ1WdT1mo3XY7D033j/xZXs7PKy28vZApgDLJf0MeLrRGBFHd6SqzUhTy38JeDUwBCyRtCgilpdZh1k7bl7xByau/Bwnvmxvvn7LSg6dtYtDxGopb4B8rJNFjMIhwGBErACQdDFwDOAAsVpo7PNoDFsdOmuX5zw2q5O8h/Fe3+lCcpoG3N/0eAh4WfMKkhYACwD23nvv8iozy2HhDSt4z+H7bgyLubOm8J7D92XhDSscIFY7W9wHIumn6XaNpMebvtZIerycEp9bUou25+zYj4iFETEQEQNTp04tqSyzfBa8Yl++fN0Kbrr3ISDrkXz5uhUseMW+FVdmNnpb7IFExJ+n2+3LKWerhnjuHFzTgQcqqsVs1ObOmsJZxx/E3K9l53yc2v8dD19ZbeUawpK0c4vmNRGxbozr2ZolwGxJM4HfAfOBE0quwawtc/eZvPG+j8KyOst7GO+twGrg18A96f5vJN0q6aWdKm6kiFgPnApcCdwFXBIRy8p6f7OxcNPK7KrQS/c8ka/fsnLjcJZZ3eQNkB8Br4uIKRGxC3AUcAnwXuDsThXXSkRcEREvjIhZEXFGme9t1q7GUVjrJu3EwO79nHXCwZx60W0OEaulvAEyEBFXNh5ExFXAKyJiMTCpI5WZ9aA7hh7jrPkHMnHDU9A/MdsncsLB3DHkuUmtfvKeB/KwpH8ALk6P3wQ8kk7sG+5IZWY96N2HzYLhYVj/R3h+dpTg3FlTvB/EailvD+QEsiOeLgUuA/ZObf3AGztTmlmP6kt/duHPXlZveU8kfIjseuitDI5dOWbjRN9E2PBM1VWYtSXvYbxTgb8H9ge2bbRHxBEdqsust02YBOuf3vp6Zl0s7xDWN4BfkV0H/ePAfWTnZJhZEc+fCmtWVV2FWVvyBsguEXEusC4iro+IdwCHdrAus942eVdYu7rqKszakjdAGmecr5L0l5IOJtupbmZFPLwCdp5RdRVmbcl7GO+/SHoB8EHgi8AOwN92rCqzXvfEathhWtVVmLUl71FYl6e7jwGv7Fw5ZuNJq8mlzeoj71FYM8kO453R/Jyyr0hoZmbdI+8Q1qXAucAP8JnnZmZG/gD5Y0R8oaOVmI0XEVtfx6wG8gbI5yWdDlwFbDz7KSJu7UhVZr1seEN229dfbR1mbcobIH8CvAU4gmeHsCI9NrPRiBQgynsUvVl3yhsgrwf2jQhP3mPWLvdArEfk/Qj0C2DHThZiNm4Mp/Ny+yZWW4dZm/L2QHYDfiVpCc/dB+LDeM1Ga8P67LbfAWL1ljdATu9oFWbjifeBWI/Ieyb69Z0uxGzc6Et/dsPrq63DrE1bDBBJa8iOttpkERARsUNHqjLrZRMmZbe+HojV3BYDJCK2L6sQs3GjPwXIhnVbXs+sy3kQ1qxsjcN3PYRlNecAMSubPAuv9QYHiJmZFeIAMTOzQhwgZmZWiAPEzMwKcYCYmVkhlQSIpOMkLZM0LGlgxLKPSBqUdLek1za1z0ttg5JOK79qMzNrVlUP5E7gb4AbmhslzQHmA/sD84CzJfVL6ge+BBwFzAGOT+uamVlF8k6mOKYi4i4AbXo8/DHAxRHxNPAbSYPAIWnZYESsSM+7OK27vJyKzcaQL2lrPaLb9oFMA+5vejyU2jbXvglJCyQtlbR09erVHSvUrLBIF/X0bLxWcx3rgUj6MbB7i0UfjYjLNve0Fm1B66Br+TEuIhYCCwEGBgb8Uc+6T6MH4gCxmutYgETEqwo8bQjYq+nxdOCBdH9z7Wb1srEHUm0ZZu3qto9Ai4D5kiZJmgnMBn4GLAFmS5opaRuyHe2LKqzTrA3ugVhvqGQnuqTXA18EpgL/Lun2iHhtRCyTdAnZzvH1wCkR2eXbJJ0KXAn0A+dFxLIqajdrW6MH4i6I1VxVR2F9H/j+ZpadAZzRov0K4IoOl2bWeY0AaUzrblZT7kOblW3Y10S33uAt2KxsG3eiuwdi9eYAMSubzwOxHuEt2Kxs3gdiPcIBYla2jftAfBSW1ZsDxKxsjeDwnFhWcw4Qs7I1dp5vPB/ErJ4cIGZl60t/do2hLLOacoCYla0vnb8bDhCrNweIWdkaQ1jD66utw6xNDhCzsjUO3/UQltWcA8SsbN6Jbj3CAWJWNvdArEc4QMzKJgFyD8RqzwFiVgU5QKz+HCBmVVAfG69MaFZTDhCzKqjPPRCrPQeIWRUcINYDHCBmVXCAWA9wgJhVQX2ejddqzwFiVgX1+TwQqz0HiFkVInxBKas9B4hZJcLXRLfa8xZsVoUYdg/Eas8BYlaFGAYcIFZvDhCzKoSHsKz+vAWbVSGGHSBWe96CzargfSDWAxwgZpXwEJbVn7dgsyp4J7r1gEoCRNKZkn4l6Q5J35e0Y9Oyj0galHS3pNc2tc9LbYOSTquibrMx5SEsq7mqeiBXAwdExIHAr4GPAEiaA8wH9gfmAWdL6pfUD3wJOAqYAxyf1jWrMQeI1VslARIRV0XE+vRwMTA93T8GuDgino6I3wCDwCHpazAiVkTEM8DFaV0zM6tIN+wDeQfww3R/GnB/07Kh1La59k1IWiBpqaSlq1ev7kC5ZmYGMKFTLyzpx8DuLRZ9NCIuS+t8FFgPfKPxtBbrB62DruVc2BGxEFgIMDAw4Pmyzcw6pGMBEhGv2tJySScBfwUcGbHxwghDwF5Nq00HHkj3N9duZmYVqOoorHnAPwBHR8STTYsWAfMlTZI0E5gN/AxYAsyWNFPSNmQ72heVXbeZmT2rYz2QrTgLmARcrexQxsUR8e6IWCbpEmA52dDWKRGxAUDSqcCVQD9wXkQsq6Z0MzODigIkIvbbwrIzgDNatF8BXNHJuszMLL9uOArLzMxqyAFiZmaFOEDMyhY+utx6gwPErGyNAPFcWFZzDhCz0jUCxH9+Vm/egs3KFsPZrXsgVnMOELOybQwQ//lZvXkLNiubA8R6hLdgs7I5QKxHeAs2K1sjQHxBKas5B4hZ2cJHYVlv8BZsVjYPYVmP8BZsVjYHiPUIb8FmZfMQlvUIb8FmZfOJhNYjHCBmpfNcWNYbHCBmZeufCHOOhZ1mVl2JWVuquqSt2fi13U7wxguqrsKsbe6BmJlZIQ4QMzMrxAFiZmaFOEDMzKwQB4iZmRXiADEzs0IcIGZmVogDxMzMClE0JnbrQZJWA7/NufoU4KEOltOObq4Nurs+11aMayumV2rbJyKmbm2lng6Q0ZC0NCIGqq6jlW6uDbq7PtdWjGsrZrzV5iEsMzMrxAFiZmaFOECetbDqAragm2uD7q7PtRXj2ooZV7V5H4iZmRXiHoiZmRXiADEzs0LGbYBI2lnS1ZLuSbc7bWa9T0taJukuSV+QOn8d0lHUtrekq1JtyyXN6HRto6kvrbuDpN9JOqtbapN0kKSb0+/1Dklv6nBN8yTdLWlQ0mktlk+S9K20/Jayfo85a/u7tG3dIekaSft0S21N671BUkgq7fDZPLVJemP62S2TdFG31Jb+b1wr6bb0e31d4TeLiHH5BXwaOC3dPw34VIt15gL/AfSnr5uBw7uhtrTsOuDV6f5k4Hnd8rNrWvfzwEXAWd1SG/BCYHa6vyewCtixQ/X0A/cC+wLbAL8A5oxY573AV9L9+cC3SvpZ5antlY3tCnhPN9WW1tseuAFYDAx0S23AbOA2YKf0eNcuqm0h8J50fw5wX9H3G7c9EOAYoHFd0QuAY1usE8C2ZL+IScBE4PfdUJukOcCEiLgaICLWRsSTJdSWqz4ASS8FdgOuKqkuyFFbRPw6Iu5J9x8AHgS2etZtQYcAgxGxIiKeAS5ONW6u5u8AR5bR081TW0Rc27RdLQaml1BXrtqSfyb70PDHkurKW9u7gC9FxCMAEfFgF9UWwA7p/guAB4q+2XgOkN0iYhVAut115AoRcTNwLdkn1FXAlRFxVzfURvYp+lFJ30td0TMl9ZdQW676JPUBnwU+XFJNuWtrJukQsg8I93aonmnA/U2Ph1Jby3UiYj3wGLBLh+oZbW3NTgZ+2NGKnrXV2iQdDOwVEZeXVFNDnp/bC4EXSvoPSYslzeui2j4GnChpCLgCeF/RN5tQ9Il1IOnHwO4tFn005/P3A17Ms5+6rpb0ioi4oerayH53fwEcDKwEvgW8DTi33drGqL73AldExP1j/WF6DGprvM4ewNeAkyJieCxqa/U2LdpGHjufZ51OyP2+kk4EBoDDOlpR01u2aNtYW/qA8v/Itvmy5fm5TSAbxjqc7P/HjZIOiIhHu6C244HzI+Kzkl4OfC3VNuq/gZ4OkIh41eaWSfq9pD0iYlX6R9Kqi/l6YHFErE3P+SFwKNmYa9W1DQG3RcSK9JxLU21jEiBjUN/Lgb+Q9F6y/TPbSFobEZvdGVpibUjaAfh34H9HxOJ2a9qCIWCvpsfT2XTIoLHOkKQJZMMKD3ewptHUhqRXkYXzYRHxdAl15alte+AA4Lr0AWV3YJGkoyNiacW1NdZZHBHrgN9IupssUJZ0QW0nA/MgG2WRtC3ZRIujHmYbz0NYi4CT0v2TgMtarLMSOEzSBEkTyT59lTGElae2JcBOkhpj90cAy0uoDXLUFxFvjoi9I2IG8CHgwrEIj7GoTdI2wPdTTd/ucD1LgNmSZqb3nZ9qbNZc8xuAn0Taw1l1bWmY6KvA0SWO42+1toh4LCKmRMSMtI0tTjV2Ojy2WltyKdkBCEiaQjaktaJLalsJHJlqezHZft7Vhd6tjCMDuvGLbIz5GuCedLtzah8Azolnj2j4KlloLAc+1y21pcevBu4AfgmcD2zTTfU1rf82yjsKK8/v9URgHXB709dBHazpdcCvyfazfDS1fYLsHx7pD/jbwCDwM2DfMn5WOWv7MdmBI42f06JuqW3EutdR0lFYOX9uAj6X/m/8EpjfRbXNITu69Bfpd/qaou/lqUzMzKyQ8TyEZWZmbXCAmJlZIQ4QMzMrxAFiZmaFOEDMrPYknSfpQUl3jtHr/UjSo5IuH9H+jTRR4Z3pPSeO4jX/Nk2seKekb6bzL0au8zZJqyXdnr7e2bRsQ1P7yENzR03SLmlSxbUqONmpA8RsjDX9od8p6QeSdmzjte5L5xHYlp1POjlujJwJvKVF+zeAFwF/AmwHvLPFOpuQNA14P9mhxgeQnSIwfzOrfysiDkpf5zS1P9XUfnTeb2QL/gj8E9l5WoU4QMzGXuMP/QCyM8pPqbqgXhfZ9ELPOXtf0qzUk/i5pBslvWgUr3cNsKZF+xWRkJ2zM5rJJScA26XZBp5HG5MYNpP0UknXp+/zyjQDw1ZFxBMR8VPamIjSAWLWWTfTNJmdpA9LWpKuw/DxpvZL0z+AZZIWVFJp71kIvC8iXkr2KfvssXrhNHT1FuBHedaPiN8BnyE7C3wV8FhEbG6W6v+Rto/vSGqelmRbSUvT5IzHNtXxReAN6fs8Dzij2Hc1ej09F5ZZlZTNjnwkaX4ySa8hmw/pELIzlRc1Tc75joh4WNJ2wBJJ342IP1RVe91Jmkx2PZ9vN03mOSkt+xuyM7NH+l1EvDbnW5wN3BARN6bXPBQ4p9WKEXGAsgubHQPMBB5NdZ0YEV8fsfoPgG9GxNOS3k02zf8RadneEfGApH2Bn0j6Jdkw2gFkE71CNjS2KtX0PuB/tihpSUS8Pef3uUUOELOxt52k24EZwM+Bq1P7a9LXbenxZLJAuQF4v6TXp/a9UrsDpLg+4NGIOGjkgoj4HvC9oi8s6XSy68ds/Occ2YScB2zhaa8CfhMRq9NrfI8s4J4TICM+NPwb8KmmZQ+k2xWSriObiftuYFlEvHzkG0bEF8l6Jx3jISyzsfdU+se1D9m1Rhr7QAT8a9OO0P0i4lxJh5P9g3l5RPwpWcBscoSO5RcRj5PNgnscgDJ/2u7rpqOiXgscH6Ob/nwlcKik5ynrKhxJi4lZR+y/OLqxjqSdJDV6UFOAPyObZ+tuYKqyadmRNFHS/qP/zopxgJh1SEQ8RnbkzYfSWPWVwDvS8AqSpknalWz69kci4sm0o/fQyoquKUnfJNvf9N8kDUk6GXgzcLKkXwDLaH1Fw8293o1kE1wemV6vMbT1FbKrbN6cjrT7P3leLyJuIbva5K1kkyv2ke2jQdInJDWOqnp/2g/2C7Jt522p/cXA0tR+LfDJiFge2VUH3wB8Ki27naxnk/f7vI9s0se3pe9zTt7nAp5M0WysKbvuyeSmxz8ALomIr0n6AM8e+rmWbGbgIbLpv6eRPlECH4uI69If+EBEPFTm92CWhwPEzMwK8RCWmZkV4gAxM7NCHCBmZlaIA8TMzApxgJiZWSEOEDMzK8QBYmZmhfwXkE8roYwSjHIAAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "a= matlab.nyquist(sysPI)\n", "plt.title('Nyquist')\n", "plt.show()\n", "\n", "a=matlab.bode(sysPI)\n", "plt.title('Bode')\n", "plt.show()\n", "\n", "a = matlab.pole(sysPI)\n", "\n", "x = []\n", "y = []\n", "for tmp in a:\n", " x.append(tmp.real)\n", " y.append(tmp.imag)\n", "plt.scatter(x, y, c=\"green\", marker = \"o\", label='Plant model')\n", "\n", "plt.xlabel(\"Real\")\n", "plt.ylabel(\"Imaginary\")\n", "\n", "art = []\n", "lgd = plt.legend(bbox_to_anchor=(1.05, 1), loc=2, borderaxespad=0.)\n", "art.append(lgd)\n", "plt.grid()\n", "plt.title('Pole')\n", "plt.show()\n", "\n", "a=matlab.rlocus(sysPI)\n", "plt.title('Rloucus')\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python [conda env:py3.6]", "language": "python", "name": "conda-env-py3.6-py" }, "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.4" }, "varInspector": { "cols": { "lenName": 16.0, "lenType": 16.0, "lenVar": 40.0 }, "kernels_config": { "python": { "delete_cmd_postfix": "", "delete_cmd_prefix": "del ", "library": "var_list.py", "varRefreshCmd": "print(var_dic_list())" }, "r": { "delete_cmd_postfix": ") ", "delete_cmd_prefix": "rm(", "library": "var_list.r", "varRefreshCmd": "cat(var_dic_list()) " } }, "types_to_exclude": [ "module", "function", "builtin_function_or_method", "instance", "_Feature" ], "window_display": false } }, "nbformat": 4, "nbformat_minor": 2 }