Skip to content

Instantly share code, notes, and snippets.

@prl900
Created July 29, 2020 06:15
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save prl900/8fa04259a7b18224d4ddf0bda4cc4e30 to your computer and use it in GitHub Desktop.
Save prl900/8fa04259a7b18224d4ddf0bda4cc4e30 to your computer and use it in GitHub Desktop.
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 13,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.image.AxesImage at 0x7f36ac165190>"
]
},
"execution_count": 13,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAlIAAADFCAYAAACfIxyhAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvnQurowAAEHhJREFUeJzt3V2MHXd9xvHv08RxSmiETbHlvKhpJCuAKuHQFW+pKorrJlCEo0qpkohqVUXyDa2CikSc9qoXlcINgosKyeLNEgSaBmisiJJahqhCoDSmuJTgBENKExNj07SIN9Uk8OvFTsROamfP/vfs7pyd70c6mjNzXubRiZ/opzlzZlNVSJIkafl+Zb0DSJIkzSoHKUmSpEYOUpIkSY0cpCRJkho5SEmSJDVykJIkSWrkICVJktRoRYNUkhuSPJbkW0n2TyuUNKvshNRnJ7TRpfWCnEkuAL4J7AFOAg8Dt1TVN6YXT5oddkLqsxMagwtX8NrXAN+qqscBknwS2AuctyAXZXNdzCUr2KU0Hf/LT/hZnc2U39ZOaGbZCalv0k6sZJC6HHhy0fpJ4LXPf1KSfcA+gIt5Ea/N7hXsUpqOh+rIarytndDMshNS36SdWMk5Uuea0v7f94RVdaCq5qpqbhObV7A7afDshNRnJ7ThrWSQOglcuWj9CuCplcWRZpqdkPrshDa8lQxSDwM7k/xmkouAm4FD04klzSQ7IfXZCW14zedIVdWzSf4MeAC4APhwVT0ytWTSjLETUp+d0Bis5GRzquqzwGenlEWaeXZC6rMT2ui8srkkSVIjBylJkqRGK/pqT5Ikndt373jD1N/z8vd8aervqZXxiJQkSVIjBylJkqRGDlKSJEmNHKQkSZIaOUhJkiQ1cpCSJElq5CAlSZLUyEFKkiSp0ZKDVJIPJzmT5OuLtm1NcjjJiW65ZXVjSsNhJ6Q+O6Exm+SI1EeBG563bT9wpKp2Ake6dWksPoqdkBb7KHZCI7XkIFVV/wz89/M27wUOdvcPAjdOOZc0WHZC6rMTGrPWc6S2V9UpgG657XxPTLIvydEkR5/hbOPupMGzE1KfndAorPrJ5lV1oKrmqmpuE5tXe3fS4NkJqc9OaJa1DlKnk+wA6JZnphdJmkl2QuqzExqF1kHqEDDf3Z8H7ptOHGlm2Qmpz05oFCa5/MEngC8D1yQ5meQ24C5gT5ITwJ5uXRoFOyH12QmN2YVLPaGqbjnPQ7unnEWaCXZC6rMTGjOvbC5JktTIQUqSJKmRg5QkSVKjJc+RkqT18t073jDV97v8PV+a6vtJkkekJEmSGjlISZIkNXKQkiRJauQgJUmS1MhBSpIkqZGDlCRJUiMHKUmSpEaT/NHiK5N8IcnxJI8kub3bvjXJ4SQnuuWW1Y8rrT87IfXZCY3ZJEekngXeVVWvAF4HvCPJK4H9wJGq2gkc6dalMbATUp+d0GgteWXzqjoFnOru/yjJceByYC/wxu5pB4EHgTtWJaU0IHZC6rMT5+aV9MdhWedIJbkKuBZ4CNjelee5Em2bdjhp6OyE1GcnNDYTD1JJXgx8CnhnVf1wGa/bl+RokqPPcLYlozRIdkLqsxMao4kGqSSbWCjHx6vq093m00l2dI/vAM6c67VVdaCq5qpqbhObp5FZWnd2QuqzExqrSX61F+BDwPGqeu+ihw4B8939eeC+6ceThsdOSH12QmO25MnmwHXAnwD/nuRYt+0vgbuAe5LcBjwB3LQ6EaXBsRNSn53QaE3yq70vAjnPw7unG0caPjsh9dkJjZlXNpckSWrkICVJktTIQUqSJKnRJCeba0Z99443TP09vVKvpFny+N27pv6eV996bOknaTQ8IiVJktTIQUqSJKmRg5QkSVIjBylJkqRGDlKSJEmNHKQkSZIaOUhJkiQ1cpCSJElqtOQgleTiJP+S5N+SPJLkr7vtW5McTnKiW25Z/bjS+rMTUp+d0JhNcmXzs8CbqurHSTYBX0zyj8AfAUeq6q4k+4H9wB2rmFUaCjuxRryS/sywExqtJY9I1YIfd6ubulsBe4GD3faDwI2rklAaGDsh9dkJjdlE50gluSDJMeAMcLiqHgK2V9UpgG657Tyv3ZfkaJKjz3B2WrmldWUnpD47obGaaJCqqp9X1S7gCuA1SX5r0h1U1YGqmququU1sbs0pDYqdkPrshMZqWb/aq6ofAA8CNwCnk+wA6JZnpp5OGjg7IfXZCY3NJL/ae1mSl3T3fxX4feBR4BAw3z1tHrhvtUJKQ2InpD47oTGb5Fd7O4CDSS5gYfC6p6ruT/Jl4J4ktwFPADetYk5pSOyE1GcnNFpLDlJV9TXg2nNsfxrYvRqhpCGzE1KfndCYeWVzSZKkRg5SkiRJjRykJEmSGk1ysrmkEXv87l1Tfb+rbz021feTXsgs/Ht74KnpZ7z+sun2VufnESlJkqRGDlKSJEmNHKQkSZIaOUhJkiQ1cpCSJElq5CAlSZLUyEFKkiSp0cSDVJILknw1yf3d+tYkh5Oc6JZbVi+mNDx2QuqzExqj5RyRuh04vmh9P3CkqnYCR7p1aUzshNRnJzQ6E13ZPMkVwB8CfwP8Rbd5L/DG7v5B4EHgjunG00pc/p4vrXeEDctOSH12QmM16RGp9wHvBn6xaNv2qjoF0C23neuFSfYlOZrk6DOcXVFYaUDshNRnJzRKSw5SSd4KnKmqr7TsoKoOVNVcVc1tYnPLW0iDYiekPjuhMZvkq73rgLcleQtwMXBpko8Bp5PsqKpTSXYAZ1YzqDQgdkLqsxMarSWPSFXVnVV1RVVdBdwMfL6q3g4cAua7p80D961aSmlA7ITUZyc0Ziu5jtRdwJ4kJ4A93bo0ZnZC6rMT2vAm+tXec6rqQRZ+dUFVPQ3snn4kaXbYCanPTmhsvLK5JElSIwcpSZKkRg5SkiRJjZZ1jpSW7/G7d039Pa++9djU31M6n1n49/bAU9PNeP1l0++tpI3JI1KSJEmNHKQkSZIaOUhJkiQ1cpCSJElq5CAlSZLUyEFKkiSpkYOUJElSIwcpSZKkRhNdkDPJd4AfAT8Hnq2quSRbgb8DrgK+A/xxVf3P6sSUhsVOSH12QmO1nCNSv1dVu6pqrlvfDxypqp3AkW5dGhM7IfXZCY3OSr7a2wsc7O4fBG5ceRxpptkJqc9OaMObdJAq4J+SfCXJvm7b9qo6BdAtt53rhUn2JTma5OgznF15YmkY7ITUZyc0SpP+0eLrquqpJNuAw0kenXQHVXUAOABwabZWQ0ZpiOyE1GcnNEoTHZGqqqe65RngM8BrgNNJdgB0yzOrFVIaGjsh9dkJjdWSg1SSS5L82nP3gT8Avg4cAua7p80D961WSGlI7ITUZyc0ZpN8tbcd+EyS555/d1V9LsnDwD1JbgOeAG5avZjSoNgJqc9OaLSWHKSq6nHgVefY/jSwezVCSUNmJ6Q+O6Ex88rmkiRJjRykJEmSGjlISZIkNZr0OlJqdPWtx9Y7wpIeeGr6Ga+/bNfU31OSNiL/fznbPCIlSZLUyEFKkiSpkYOUJElSIwcpSZKkRg5SkiRJjRykJEmSGjlISZIkNZpokErykiT3Jnk0yfEkr0+yNcnhJCe65ZbVDisNhZ2Q+uyExmrSI1LvBz5XVS9n4Q9THgf2A0eqaidwpFuXxsJOSH12QqO05CCV5FLgd4EPAVTVz6rqB8Be4GD3tIPAjasVUhoSOyH12QmN2SRHpK4Gvg98JMlXk3wwySXA9qo6BdAtt53rxUn2JTma5OgznJ1acGkd2Qmpz05otCYZpC4EXg18oKquBX7CMg7PVtWBqpqrqrlNbG6MKQ2KnZD67IRGa5JB6iRwsqoe6tbvZaEwp5PsAOiWZ1YnojQ4dkLqsxMarSUHqar6HvBkkmu6TbuBbwCHgPlu2zxw36oklAbGTkh9dkJjduGEz/tz4ONJLgIeB/6UhSHsniS3AU8AN61ORGmQ7ITUZyc0ShMNUlV1DJg7x0O7pxtHmg12QuqzExorr2wuSZLUyEFKkiSpkYOUJElSo0lPNpekwbr+sl3rHUHSSHlESpIkqZGDlCRJUiMHKUmSpEYOUpIkSY0cpCRJkho5SEmSJDVykJIkSWq05CCV5JokxxbdfpjknUm2Jjmc5ES33LIWgaX1ZiekPjuhMVtykKqqx6pqV1XtAn4b+CnwGWA/cKSqdgJHunVpw7MTUp+d0Jgt96u93cC3q+o/gb3AwW77QeDGaQaTZoSdkPrshEZluYPUzcAnuvvbq+oUQLfcdq4XJNmX5GiSo89wtj2pNEx2QuqzExqViQepJBcBbwP+fjk7qKoDVTVXVXOb2LzcfNJg2Qmpz05ojJZzROrNwL9W1elu/XSSHQDd8sy0w0kDZyekPjuh0VnOIHULvzxcC3AImO/uzwP3TSuUNCPshNRnJzQ6Ew1SSV4E7AE+vWjzXcCeJCe6x+6afjxpmOyE1GcnNFYXTvKkqvop8NLnbXuahV9nSKNjJ6Q+O6Gx8srmkiRJjRykJEmSGjlISZIkNXKQkiRJajTRyeba2K6/bNd6R5AkaSZ5REqSJKmRg5QkSVIjBylJkqRGDlKSJEmNHKQkSZIaOUhJkiQ1cpCSJElq5CAlSZLUyEFKkiSpUapq7XaWfB/4CfBfa7bTNr+OGVdq6Pl+o6pett4h7MTUDD0fDD+jnVieof/3HHo+GH7GiTqxpoMUQJKjVTW3pjtdJjOu3NDzDcksfFZDzzj0fDAbGYdiFj6roWccej6YjYyT8Ks9SZKkRg5SkiRJjdZjkDqwDvtcLjOu3NDzDcksfFZDzzj0fDAbGYdiFj6roWccej6YjYxLWvNzpCRJkjYKv9qTJElq5CAlSZLUaM0GqSQ3JHksybeS7F+r/b6QJFcm+UKS40keSXJ7t31rksNJTnTLLQPIekGSrya5f4gZk7wkyb1JHu0+z9cPLePQ2IkVZ7UTG4ydWFHOQfehy7QhO7Emg1SSC4C/Bd4MvBK4Jckr12LfS3gWeFdVvQJ4HfCOLtd+4EhV7QSOdOvr7Xbg+KL1oWV8P/C5qno58CoWsg4t42DYiamwExuInVixofcBNmonqmrVb8DrgQcWrd8J3LkW+15mzvuAPcBjwI5u2w7gsXXOdQUL/8DeBNzfbRtMRuBS4D/ofrywaPtgMg7tZidWnMtObLCbnVhRpkH3ocuwYTuxVl/tXQ48uWj9ZLdtMJJcBVwLPARsr6pTAN1y2/olA+B9wLuBXyzaNqSMVwPfBz7SHVr+YJJLBpZxaOzEytiJjcdOtBt6H2ADd2KtBqmcY9tgrruQ5MXAp4B3VtUP1zvPYkneCpypqq+sd5YXcCHwauADVXUtC38na/YOz64tO9HITmxYdqLBjPQBNnAn1mqQOglcuWj9CuCpNdr3C0qyiYVyfLyqPt1tPp1kR/f4DuDMeuUDrgPeluQ7wCeBNyX5GMPKeBI4WVUPdev3slCYIWUcGjvRzk5sTHaizSz0ATZwJ9ZqkHoY2JnkN5NcBNwMHFqjfZ9XkgAfAo5X1XsXPXQImO/uz7Pwnfi6qKo7q+qKqrqKhc/t81X1doaV8XvAk0mu6TbtBr7BgDIOkJ1oZCc2LDvRYBb6ABu8E2t4otlbgG8C3wb+ar1PDusy/Q4Lh46/Bhzrbm8BXsrCiXsnuuXW9c7a5X0jvzyRcFAZgV3A0e6z/Adgy9AyDu1mJ6aS105soJudWHHWwfahy7QhO+GfiJEkSWrklc0lSZIaOUhJkiQ1cpCSJElq5CAlSZLUyEFKkiSpkYOUJElSIwcpSZKkRv8HahSghaydn9IAAAAASUVORK5CYII=\n",
"text/plain": [
"<Figure size 720x216 with 3 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"%matplotlib inline\n",
"\n",
"import numpy as np\n",
"import matplotlib.pyplot as plt\n",
"\n",
"\n",
"f = np.zeros((100,100))\n",
"\n",
"f[45:55,30:40] = 15\n",
"f[55:75,50:60] = 25\n",
"f[25:40,45:55] = 10\n",
"\n",
"\n",
"f1 = f[:-20,:-20]\n",
"f2 = f[10:-10,10:-10]\n",
"f3 = f[20:,20:]\n",
"\n",
"fig, axs = plt.subplots(1, 3, figsize=(10, 3))\n",
"\n",
"axs[0].imshow(f1)\n",
"axs[1].imshow(f2)\n",
"axs[2].imshow(f3)"
]
},
{
"cell_type": "code",
"execution_count": 9,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.image.AxesImage at 0x7f36ac53e850>"
]
},
"execution_count": 9,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAP8AAAD8CAYAAAC4nHJkAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvnQurowAAIABJREFUeJztfXusbkd13299j/O474d9zcX3go1wDJSATS0CJSEU6spQhFFTEohSuWlUVDWtQCLlkT9a9Y9KSEUpqEJpXYeUKORBCAiCIhKXJC1RKa+almBzMdjGvvhxr+37vvec8z2mf+yZPWvOt+bMzP4ex+fu9ZOOvn1m771m9mP2/GatNWuRMQYKhaJ96Gx3AxQKxfZAO79C0VJo51coWgrt/ApFS6GdX6FoKbTzKxQthXZ+haKlmKrzE9EdRHSCiH5ARB+cVaMUCsX8QU2dfIioC+D7AG4HcBLANwC8yxhz/+yap1Ao5oXeFOe+GsAPjDEPAQAR/QGAOwFEO393z27TO3QIAPCTB08DAB48cbDev7HfN+eF152qt0+ePAIAGOzmwthHa0x+227SiB3LDu3sqnaYi12/u8sO7UXkOvEjfqwr9OfQBk3uj8FEth0f6zC5AybXtZc3LybLbfPrGvkTWdNhUjxQGic6XACTa+9TcA9i48zYncTKupPXHrSPb0tyY3XZ+xA8J5rcD8DfJ/56Bc+fHWvb273oGzbaxe/NJplAcA2BXPt8d+9e87Ie8gdvHPcv7LHVMwCAp/5mFQCwhkvYMOuTL66AaTr/9QAeY/+fBPBTW1Z26BCe/773AgC+/gv/GQDwljf8XL3/x2++rt7+L+/5T/X2r/3rfwEAePK1/ppG+/3d6vCO3K/ucu+8v1ld9qBXbn0WALDxvw/VZYP9/iENDg293EuVXP7SLZ3zsjYOV2+tWRrXZSuP9evt9WvYE3Xv0ZA9F/5RYp17vFztGO/y56887uVuHLD19uWPDn+RXH3Dg/66uux+ddh5wz32OthHj/yl1XKJ7R/uZc9h3d+o3sXqmMFBJmAoy+3Yd5Vfz+iAb++SvXbemdw9AsIPYy1zIHfu0TUbAIDlx5Z8md/E8DC7T+fs82eDg7suANg4zK593wAAsO+rq3XZ2Vdt+BMH1b3pnfPCxkv+GvoXmFz7Pr7mthN12blf3FNv/+gj++rtj7ziMwCA//jilwIAvma+jFxMM+eXvi4T31siejcRfZOIvjm6eGmK6hQKxSwxzch/EsBx9v8xAI9vPsgYczeAuwHgFa/omy/+3EcAAC//2PsBANfe6L+OP/kLfsbwrnv/eb194HnV13K8xIa0Hhs64L+m3UvV96zPvtArr3263r7y9WsAAMO9bLQ/ODnac3QGfnvMRgmzWrWnf8qPyqNVL5ePZJ21ye9sh42EnF2YA1WFKw8t12Ubh/z1ji3T6LK2BhSUQSrvrjEGtTzZ3u4F3xh+vhvxRyucDvjNpTP+vPUjjvfL9fJRz9XLR1isT96vMXtbebt6nGovT17viLOTM9WzIka/h9f4d5Au+UocexixKchwN5PPh79T1bNaO8wbzCl+dR4f4dcPM7mM1Vz3E9WU+OGP31yXPfvP/Hlvv/Gr9fYHPv4rAIAjP325KrjvfyEX04z83wBwExHdSERLAN4J4AtTyFMoFAtE45HfGDMkon8J4M9QDb2fMMZ8d6tzHl4/jLu+/4sAgIPfr77GT7zOj5p37P1xvf3Ivf6r99Rt1Vdx9/ELddnw2wfqbTcHBlB/ja8c9V/7K496pWJvtzBCslGksz459+ajxGCfl0t27kyMGQxYW/jc2M2tgxGvK89hO0/3bV3suvgl2klsoKzj18MmuW7U6531j3rsbznGjKl0LndcBZAwdveJDRm9s3643tg/Ob8vUSiOGaNYOu3lOpbGdSu83tHSZHvHXB/Cnl/vSrW9zpgUNph+KMJOvDC2zd6bXU9WMi7ezOb5rN5dj1Q3/fIxxjIZu+HP/8mHK/qw74jfPzjulX9fPXVjvX30K+cBAKf/bVXv8L2cEW+NaWg/jDF/CuBPp5GhUCi2B+rhp1C0FFON/KUYXOjj8a8cAwAsvaAq+4mffaje/3v33F5vX/gZRl+61fblR7yJg/bJ9MaZkAwzk3GFnKNtA2bS6Z32+8eCwihGIbtWiTeUphIAOlcSFDICVwdXanFK3DtvFaArMu/n53WudCbKAgSmvMl6w4bZH66o7E7uB1DbvTtcGcemGNzU52RwGzlXyNYmRm7S48OWYHfiUwz+HGqFIJuPdDa4ttVvuufOnzmfqvWf9t3n0gurRnbPsOkVe+bOfBfUxdvL3rHVxyoZF25iikr2TM7896P19hP/qCpfHZ2t5BjhZkSgI79C0VJo51coWoqF0n7TA9atV9T6oYqenPhrr7ncuGW93l79obdxX7mh0mR2L25NzwFguKvik5zKD5g3YM9q6LlXIKfPgZul+zRytrkxSeUDV2JOESUXVO7gxyhzh9veV6tr4F5q3E/AtVds62bYegM2GLmemu4b+di6iNfb5zyZybVWE071O36mJd4bXsbrcPSba8e5PT+4/4J7dzA1cayfezGydgXtFdyK+RRgxN4b518y3MOmm+e4xcI+03X5/eBTnivXVw1aOsMsKV3mOXqIeUJauVdOVNav8dqkn0oMOvIrFC2Fdn6FoqVYKO3HGOheCd1vHaUHgM5ZRtX3MffNZ6ry4QFG3yMLJByNdXQICN1gxwK1DWid0Gyumeaaf1eXSO8BcdUWp98BteXXsC5o6HnDnEUjptjl19YN27pZbkCJU/4h9jTuZhuVa689SvUFubwt0kKl4B4xxyojMd3INMdZKgxbjcidnrrrk9O66ApBPoWz21ybP2TvYO+CswwxJzA2hQjaYN9XvuCIW1j44rN999u+8fpzVf0rfA60NXTkVyhaioWO/AT/FXcumf1d/hNOj0/a4wFg4/qKHfSe8p/CQNHEUI8ufNDkI5rdPxaUNUA4sk/IRKgoqssE2zAQKukk23kw0vF13oJ7LR/lO27kj+l2JMVZJ3K/GsgNRvNuhJ3kuzVky6XYoCbVy2/nYJKd8JE/UOLGfByEuiTlL283f69qBWWMrrFn7lyMN64difv3/MD3k/OvtMz5qSrYxXiQP57ryK9QtBTa+RWKlmKxtH8ELJ+1ir5b7PrjH/vIJz1GqdeP+umAW4PNP1WS7ZYfE9hTBQVXEMapm0/1JJtyQPXXC+QaeVvyL8iy6QuyHFsM4gVw19eU4kyQm+NBSimlpCA3gHCNgUtwKjxY4D48eYA0fat2NJPrphHBFIPfcxc+jEdyikwhnGszMb8GPhPcYJGnMHCrJ9356t6rUCgS0M6vULQUC6X94yWDizdWKt3+oxXdHzJ3Sh6UgwefcGGjgijjXIOfoGoS9QwoZMxO785LHBuscuN3ND+uwqbGCW1J7Y/cg5T1Q5w+5TNH0f8AYPc89myE8uDeCs83tT8LKfdfDul+CLZ9gLkCB/eTuSC7dyTS1pQrOJ9O8oAxLhhrbS0qeOd05FcoWgrt/ApFS5Gk/UT0CQBvBXDKGPNyW3YIwB8CuAHAIwB+3hhzJilrQFj9cVXlleMDW+a/P/2zLGYZj4XnKFXED0ai9dwFlbts1s4W3MEmcC6RGs62JU16iYtsDIJmn7t8BvHznSY99ulOaawjVL0JAtdnds9rWp1YIRiU5xyb3TC2yacjLjUBd1HmFg8pMUjEKhM4d1l5wXRC9BWXZUma/9Ee3/DRit/P+wnPHVCdnO9dlTPy/zcAd2wq+yCALxtjbgLwZfu/QqHYQUiO/MaY/0lEN2wqvhPAG+z2JwH8FYAPJGX1gPXDoQGYr2MP7KLdSWUJH8Hdun1gs2ts9RPE2pdcNmMjXmLULIlGKx4bi5Al2LADt+QE00m1IeYXUZSuSzhHHO15vQUTy0Bu4Po8uT/KthLK0nqhE2N7wTOXlIo599s2LuY/4OJPBAuhArdxdqwQ8q13mbNXdp7NHbH8cFUoZS+Koemc/zpjzBMAYH+PNJSjUCi2CXNX+Gm6LoXiuYmmdv6niOioMeYJIjoK4FTsQJ6ua/n4ceNtmHZVH0sNtXad52IUhLCqtrlPQED1gwrdr+xm6+hV1FW0xH4sUMyk3Jg9vmS64Shvhq+DSdDgsJI8uTG6yiHa+UvkCiv4SqY7wXnCFCIW1ix1n2JTHinklxhHICcGg6X9PLQXdyEfsSlv93S10nVgYwckXbQZmo78XwBwl92+C8DnG8pRKBTbhGTnJ6LfB/BVADcT0Uki+hUAHwZwOxE9COB2+79CodhByNH2vyuy603FtXVMnXyi/2xVdZ3NFT6yLhBqPOsouRHtt0gh+Uo7rvnvh8dVhRFZQl0pm35SCx1DIiCFkehzlhZ6k8xNEF1mU3JLbPdTuuEG9UWmBaJFITbdcOkIY+cn2hubmrjymNt4Pd2IDLdBJOcrk3b+HpsejwO/Yts3YtPgLaAefgpFS7HgAJ5UB/AcXGNjkz/pm8CDdnJPpY6waCG6ntutnxjKx9Zf65gHWMIrS5QbU/ilPq0xRZM7L8cXIbG/9gaMMBJxkU9s9CtRHkqLYmKMQZIrLXSJLUhK2fkLlIui3JJr4BACvwZ2/ojvhWOnQS4HHnaOsYChTTu39FT1MsYUsBJ05FcoWgrt/ApFS7FY2t8xGNm45V27Xp8vXgFzw+V2fmcjdWmsgAw7Py/iNLieFvhCrlyUwnRFFX5mcn+AlJIufw2GLDfHzp+YQqT8C8KDhbp4Exq4CgdyY4q3ugLhnBwIzyHMxcAOzQ97Ly7G4vcz8C+o8wVERAXvmBDGiymth0f8P8snbdx+mym65L7ryK9QtBTa+RWKlmKxtB+oaZezvY92cb7jN3nYotq1MeBUTKawzj/Q8Etr4fkUY9o17TG6OyU1LdpfID8axislN+UDkdiflJs6tuQeJKYIASUvuZ8FsRL4e5Fa5RjY+dfdFIH5ujBtf+ecf7mHq5vkFpj7deRXKFoK7fwKRUuxWNpvqM5AO3Ir9Dg1Gsga+Hp/hqtoKmhGkeNNoq5smaVyS1yBG+xPRisuwbRtLTlvFmG+5tRe8R0osX4kjo1Zt6QckLnQkV+haCkWO/KTX3BTuy4yLdGYrVN2ixvceUDcHTJAQrnjFI3JWO2s3mhdJW64Je6urijlKpxwYY3Jbbw/dQ0c85CbsRa+JKdBo/0xSOeV+C1wUfa8LnfvXeUvgy93oe3qUHUFREBHfoWipdDOr1C0FAtW+GEiEu/GQe9PyWP4c/fb0eokl+mkYu2n1n6nVpDFIK3RJnl/45j4gltDYzh35pwsvyV29tz9JXZ+XpRSgKWmRAkX5Cw325LIwy4fQCrXQiq7MNh7xZV5kaQVzidguM850OTzfh35FYqWQju/QtFS5KTrOg7gdwA8DxWRvdsY87FGKbuMt+W74AQkaC4B2c6fpPq8qslIR0F5YAHgclMUsqk7awP/gSh1LUnXJR2b0van2jXtasQMuamstUWJSwqsJo3lCum6gioS/idBuq66j8jzxh6L6uvSdXUu27JI0hAJOSP/EMD7jDEvBfAaAL9KRC+DpuxSKHY0cgJ4PgHAZee5QEQPALgeDVN2OdRhjZg9v06iiU3plCQFWCKcUpQluPXcMQVYiQ1cOq7Elh2Tm/AvEEeRktE61d5Z2Pkb+DVM3a7IsSkWkWxDKlYC4JW7JYrmyEIl1w86bD0/Z8JBSC8bFq9e/DavhT02Z9+tAL4GTdmlUOxoZHd+ItoD4I8BvNcYc77gPJ+u65Km61IonivIsvMTUR9Vx/+UMeaztjgrZddEui5ro+xaBcVwP4vbf4HF7Q9snPY3RqkEWhakUkop9HLDV0Xa0DT1UzQmQUqhN8MpRNHUpMTvIEXVm2AWU4iE3CL/AUluZFogRlGOXEPtgs5yT8QWtTkFunHuv5S6WFZP6gAiIgC/BeABY8xvsF2askuh2MHIGflfB+AfA/gOEX3blv06qhRdn7bpux4F8I75NFGhUMwDOdr+v0acsJWl7DJ+XbJz2eVr+DmCbKhLzjLgiUp0HbNzZy0I81RE9VIJPvhpKVfREv+BEjSdeuQih8rPQ25D/4KpXYVL5PL3KpHgJQUexouY/b7LLGQD6x7fvTQfO79CobgKoZ1foWgpFh+91znsLFdcq3feqz+DaKXc4ccGNXDJDErqASLa+AL6HgVt+kXm6rkmaB6tqULONc7DrbdpIIzU/qZOPpL4kucfkyutnkwleIlN+7pOlj8gSDLCk8xEps050JFfoWgpFpyuy8cf71qbPnfp5V9N4tlMXSubro9vGsapQUiomQTHzLXNz8JllyNl50+hRGEX843Y6liT2F+CgmcabavQnuTzn1aZi00+LC7Gv2MZ83LvVSgUVw+08ysULcXCw3h1LJ0f2Ui9nN53eOgutpaZh/TKRsyOL8QBmNqlEwX7c5BaETftyrSmKKHfKfrMUTIVy0VqGlTiKswRa0sqknMB6pBgXLHHYl10WJ9xEa/r1bEF90pHfoWipdDOr1C0FAtP2sG1+8AmCrMscxZpoVIqG26U6gv21qLQTQX248ZUuyQZSBPMwoow7f55r/Br6iosuXIXrHKk2Hsl+QFE3h+nzQ+mu0yY5Nru+lHBoj4d+RWKtkI7v0LRUizcvbem447GcGbD3Hf5Kia3UinQfvIoqQIFjMb7c6cI7pRVJYLcWay4S2nwm2i3mzr5TIuctja5T02tBSVIyS3R1gvnpRK4FE0xA1dh33Ae5KN7wWa93lWVlbiU68ivULQUi8/Su+Tce+0XazfLzMtdejkjcJF+U+G4NtXlD54sN7GvferLPC8FV5OwVE1chreqK9d/oISx5NyDecidFxNKvTcxHwjp/Ei5e8956q+gbzBG4JTk9QIfVfgpFIoUtPMrFC3F4rP0Wvoy3Fdxl6UzXts22M9cetk6ZRcOKZWiqzrY/sTs/JJCp4Dqi3K5zBJK3dRtOEWTS6h8k1WMTd2Km8rNPZ+Xl9TVdLpRkpG5YGrj3isXxRcIlXx8+luHuNuYA+0nohUi+joR/V8i+i4R/TtbfoiI7iWiB+3vwfxqFQrFdiOH9q8DeKMx5pUAbgFwBxG9BpqrT6HY0ciJ3msAXLT/9u2fQZNcfeTpS9+G7xruZjyFhS1KBcoosvNLTYm56abs/CVUL4USa0CJZr9Eu71IV+Fp5U670i8mvwl9z2kDny5KCV5icmv53NfFV+ZWxALA0rNVPxrstWWztvMTUdfG7D8F4F5jTHauviBd18WL0iEKhWIbkKXwM8aMANxCRAcAfI6IXp5bQZCu6wXHjRvd3Uct9iUM1vY7W2bJSJv4ms9kBG9in5bO34xpYwZM64XYVG7TeqeVW2Lnz92fU690XmTBWf2eZ9yjmiVwrz62mKeO0Q9g6PJfFMTrdygy9RljzqKi93fA5uoDgK1y9SkUiucmcrT919oRH0S0CuDvAfgeNFefQrGjkUP7jwL4JBF1UX0sPm2M+SIRfRUNcvVtDlEUo/LjWIbSujBWgf2dpTKt6f7tosFNz5mH3OfCPVjksTNQrIopxZgyXPR3aaAAzdH2/z8Atwrlz6A0V59CoXjOQN17FYqWYvHpuixqOz03ZXIXRr6KyX2iUvZ4IOkq6qYQ0VWB04aESmnwm7qVltQ7rZ2/qWusJHeWmvRZauhLLC05lH0OcnkYr/FKxPXdTgfGy+UmKx35FYqWQju/QtFSLH5Vn6X7jtIMD3g/XZe/rzqAOTjUWUsL6hJcK4GyMEq5crOmEAlZRfsbyM2KOlyyam8e++cldxaBO3Ll8nMYE3f3X4wkvVmEm5oyxx6u7ZdWmdar/grC9+rIr1C0FItX+Nmv1si6JXausO9Ph432TUevXEVTjBmk0nVJqb1KQlXF0GR0KlBEzSR7cMn5Ta6nqZt0CaT3KpaurSRcWuK9q+XmhIwTnhl/97ssddfAMueeXSiHAjdfHfkVipZCO79C0VIsXuHnbJd9m6WX0RQpvREQUawV2N7F81L7Y3KlOkps903t/E38D3KOnVcIrNy6SuTO6RqiOR6m9YGQkDFNMrZXBmG8WE/lKe+ckrxe469x+xUKRQra+RWKlmLxSTus7bKz7pJ2eFV776I3gnJX32TWWg53WiTkl6jNn5bqIWN/ihKn0IRizlpuypKSOnZaubOebqSOnbYNHAU+FC5adZCejrn6BlYv5wvg/AAKLCM68isULYV2foWipVi8k4/LuLtacfguc/IxHZmzOG1sJ5aIg0OYIszSpXfq+Ho555UEiciVO0tnHI6SVYzTyp3WOSlWVYmTT8H9EjNBp1aj8jYYmepzV9/RajWn7V7qFrdPR36FoqVY+MjPI5ICm9YsMwUHX9TQWetMlMUryG1I5nFbYSctOMnZ38TOPwvkspMc9pKraOQj7QzZSTRNXMmzpklZ/PxAEThoPn5nn2lj999HRF+0/2u6LoViB6Pks/EeAA+w/zVdl0Kxg5GbsecYgH8A4B5WfCeqNF2wv29PC/J/nXVCZ50wWjH1H43h/4ZU/5meqSi/gf8rAatXPF/az4/hZSlZiJzXBLF6U/JJ+EvJQsb+Ju01kT+pvbFjt3o2m6+3CVL3KedYC9PxfxizvyZg8gO5DJ1B9Wf6plrTX3APckf+jwJ4P8LLaJCu61J+yxQKxVyRk7TjrQBOGWO+1aQCY8zdxpjbjDG3dffsbiJCoVDMATna/tcBeBsRvQXACoB9RPS7sOm6jDFPNEnX5dx3eZTeYKVVLCNvA4ga2IgmdWpX0aKGReTm2pdLbPcNbdWijNg9aGoZmFd7pfOnXanXVFYJrIwgP9+a7xBjVj5eqbbr/H0FU4zkyG+M+ZAx5pgx5gYA7wTwF8aYX4Km61IodjSmsfN/GA3SdW228/NsvMPdLDY582KqgxTGvPoEiN5VgPxlnNbOG/valowi0nlN7fFNMa2dv8S/YNqFTrOUO4t7l8lOgvcyUq97Xzkr5qM9f5+5n0wpijq/MeavUGXp1XRdCsUOh7r3KhQtxcLde13YLrdmecRdFXls8kABVpXzhT+BElCgptHQXCmUuIo2td9KaKrcy5WbmkI0RdOpSxO5s1ACTjs1KQHXZadyPPBjx5NlAVg/cIq+oYbxUigUudDOr1C0FNsexmu8xDT8kVBFbhVTZ+jLovSpPiCsV2pL9rEx5GrHm7YhRlFTWvemVoZcuTnXW4Jp5c7y3qbuXcnzbZrgJTF15f1ktBL2p5JnoCO/QtFSaOdXKFqKbUva4ZIM9C7478+AZ+xl4b1c8oJx34uKavtztcwlVD22X6Jys3QVLaGuTacxTeTm1DXtdKOE6s/rGiRIzx9Mm0/y/qRlSJBLfJrLHdUki1NBdt5N1SgUirZh8WG87Ferd9HaJ/f5T2JnTf4WOcVHzsdNdAXmyhJJmZJCiZ1/Xq6i0yrWShRcTSGNijm+EE3uWRMlYOz8khwPkedvJJbIT0u9d4nnEJzHOkLvsu1H+9XOr1AoMqGdX6FoKRav8LO0auxqjsQxD6L6rownyqJVJJRw4kqqaSlktDFTyp3WpXfW+5vYyEvqbeqC3ERuir43xQzkplax1jZ9AKNd9uJc31A7v0KhSEE7v0LRUizcvdfRIufmG6zkY+BZeutjtosGXw1yc87JldvEBfq5fOwsn9MMVgVKFoEgFB1PXmM3G5j5deRXKNoK7fwKRUuRRfuJ6BEAFwCMAAyNMbcR0SEAfwjgBgCPAPh5Y8yZLQUxbX9N/9nnh2/z+GU1zYk52BS4byadLZq60Ur7JbklrrExufPSpE8rt8SSst1yc5ybSmh7YkVkrcHPGW7d687e8fGyb2RnfbISF8W3pM0lI//fNcbcYoy5zf6v6boUih2MaRR+dwJ4g93+JKrAnh/Y8gzjlXfuCzje7f0p+2e8YdOtUwZQR/ztDBgbSNmiuc9ALP86a9eWaOqyWbJwp6De7P3zkjvtuv2mcufFIhq6OIv5IExkf8Ew684LFrLFfFxqBu0qzr+A3CYZAH9ORN8ionfbsqx0XQqF4rmJ3JH/dcaYx4noCIB7ieh7uRXYj8W7AaB3QLN4KxTPFWR1fmPM4/b3FBF9DsCrkZmuyxhzN4C7AWD5+HHj7PfGhubqnWdUf5Xb9pmMjkB5UlSvwLUySKRQEB6sKDRTCk0VXLlyZ7l6j2MWCs4SudL+Ke+/SN8BtlZersuU3FN3bEpRjfTUlPvGDPdXHaV/tjuxL4WcRJ27iWiv2wbw9wH8DTRdl0Kxo5Ez8l8H4HNE5I7/PWPMl4joG2iQrkuhUDw3kOz8xpiHALxSKC9P12W8/X5kbfc8BxkPWzRmyTzcOYa3tqGdXyqLJvjI1QyX2O6bHiu1oSSM1yyPLaHcOZR8HlOEgmsInn/svWpQryg3h5Xb82LRqrnN302bh3urCkryWaqHn0LRUix8YU8dg9/a7MfMnt+56L9FXOHn7J0ukCewxRdO0g1KSpoSP4ASlIxeTRRG/NgUM5iX3JxRWTo21YaS0b7pvU21q2C0T9n5m7IXx4AD9st8XDD22y4QbmfNlhWkkNORX6FoKbTzKxQtxTZk6a1+R3ts3P6znr9zmuNceoEtXBvFCuxPLAqrtL9E/LSuwlmVJOQ28SVIKduayo1h2oVMAqL+GCW0PlkJ207Z+ROZoMX2xsLWBev1XZlM9XnfqN9nDeOlUChyoZ1foWgpFq/ttzZKWp/87oThiTzN6VhKw6cFUXrj1k3n0FzWri0Rc+lMtWVaG3lTTGt7b+p/MO2UJyG3sSWmqatwiVz73gVUX3IVjqDEVdjwEHcbzV8cHfkVipZCO79C0VIsPGlH7apr2QpfyReG7vKn1S7AJa6XW7QB2IJCSpS3qVvqtLS+RDtOke3UsSXXW4JpHYKk65n16snU/pTzkPA+BuKl97VpsJDADd43wlmtRquaq0+hUGRisSN/x7vzknVX7LJghIFCrymkLyzfdhl/Y3b+Jgq9WdiUUyNwU9fXVF1NMOvFOrmKxFmM9tLz58dOqaQtig2RwJj1zu6arywIcWcZQe9CVRmpe69CoUhBO79C0VIs3L13M+VTzd/sAAAQW0lEQVQM1vDv9Zyle5kvYJ6083PaJtEyTr84fXJrpAOqP+167hiNTq3nLrG3SyhRppWgROkVw7RKw1m4STsIzyE51csBd88dJeQW0HHnvtu94hvO1/Cj47c7a9Uxrm+UuKrryK9QtBTa+RWKliI3XdcBAPcAeDkqQvZPAZxAg3RdsNFFnfZyuM/zoQ6jOZy+UCefA9aJNJhrpajZb6r9jk0R6soayuWYh7a+xAV5FnKl/dNOBRI29uixBQFNoqtBE3KTEaLddCMWhIanqnOBbvq8gshpdho73Gcm5KSQe+jHAHzJGPMSVPH8HoCm61IodjSSIz8R7QPwegD/BACMMRsANoioPF0XUA/N4271dQvW6ke+wHXoL+4BGPmCiuum2SdOVMxM6dUVHVnmFWoqhRJbdep8qbyp5+G0HnwZMe+zFypF/EDENfo5C5kkmz9nEa6nFcTtj8WO4Ik6h7ts33AK8hnb+V8E4DSA3yai+4joHhu/X9N1KRQ7GDmdvwfgVQB+0xhzK4BLKKD4RPRuIvomEX1zdPFiw2YqFIpZI0fhdxLASWPM1+z/n0HV+cvTdb3geE1unFsid1scHPA8qXvB83pHc8YrTDnIpwApqpaiVLF1100WlAiupIHcpu6uJdOCplOIXFo/C5ddDnfsLJWp84pJEHm+zl+Fv0uB0lnYn5QbmZpwm3/3krXzO5ffWdr5jTFPAniMiG62RW8CcD80XZdCsaOR6+H3rwB8ioiWADwE4JdRfTg0XZdCsUORm6X32wBuE3aVp+tyq/ks3d+41vv3ds+x5jD/3DqaKbMMxOyq7rRxxM7v6L5LBLJ5fxIpl80Udc2Razb9lsqd11r3BlOepN2cy41dV66bdcY5tR9IxA03GZ055hMg0fmC6SI/nyenqffzBB5X2Hp++1tPiSlnfmflZB+pUCiuKmjnVyhaioWH8XLU3TknBFF8xaR6gOltykeGOIWUtPlSRNWSoAeNHVlix0hlTbTbJZr0kvOaBAvZfF4qinIJUtOCFBIBNrKmJpJY4fnE6Hsqe660SjU4h4fuYttDm6uv7kcFN1xHfoWipVh8GK9NoboCJV7wdWRr962SsCT3eCAp9YmbwXruqeXupGPntZ4/ZXuf03OaZT6AqP9Igax68GZdhSJ2fho1p1Y68isULYV2foWipVgo7aexj9Y7cmGHovHI/babKnDqE1XSCK6i4nr+HLt5E8XXDJR09drvnNwCuXJzbM4N5Ga5SWfKirYld6VejiyHeT3/iFwpzkRMcVqvPOWuwmxtv3PpBXz/Ge1WO79CociEdn6FoqVYePTemoI7s+Sq50bLJ73P7foRHs3DnsIi+gbRTAuCNUifu6aJFkRX4Zh/QUpuSUgoRyFz2pqi8k3oeSr4BZdREJAk5VqbdL2NyRX8PGYRHZi3p56a8qkr39+fLAsQvK/WutVnq1jXSTy2ToLjVrmqnV+hUKSgnV+haCkWSvtNBxjutjH8Vitav/yEbwKn+oErb9f9pqOZOjoYZvzlZgL7E8vVl0CYOXVSVtCWErkpF1SJ8pZovEtckFOx6TJovxitNhG/rrErcGK6IcqNOdDwhaXuucYCw/BjHd2PrBbcsq2b2yBNIUby++ymBsvP9MI2Z0BHfoWipVi8nX/NuepW353BPv8V617w36LhPv8J652vho/RKndr5IIxUc7totK662hapdQoIkUCjuUIKAkvJsjg54vRhnPs9U1CkSXkmpiNnOukpCjKs6y36fUI4dT4ve2w9HElSlqJUXJFcJ0mLrEgLZDB28hYL3eR7z1bdeH1a6uTApabgI78CkVLoZ1foWgpcpJ23IwqLZfDiwD8GwC/g8J0XaYDjJx93mbeNSwVFzGeQwP/XRrtqbha9yIr41MAiVpyBZlEmYXIq0CoxElS9dS0oGlIqII0USKEsFRFPgExGu2o7axdelNypbBmBVQ/sLcLSlrRDwCQIy4n5PI28vX8khIv9kxr/xH+jjMfl94lf+LgSFXJ0uP9SfkJ5ETvPWGMucUYcwuAvw3gMoDPQdN1KRQ7GqW0/00AfmiM+RGAO1Gl6YL9ffssG6ZQKOaLUm3/OwH8vt0O0nURUTpdl4HX3O6tuPbqw0v17rVrPRfj2vreuYqLcS1ndPGSZT0uSjDAphoMgR9ARjZUqV4pWENyCpEDKZhD07xvtPV+yW9hFvZ40QU5JTemzd8kM0dubArhyiMR45IQqT6XG5260Jb1cjv+aHf1sLssazW3QgwO+peBLrh5DMLfDGSP/DZm/9sA/FG+eITpui5dKjlVoVDMESUj/5sB/B9jzFP2//J0XTceM+bwBgCg//gyAB/IEwDMsv9s8nRd9X4+8m/IihM3QvLc5uJiioJ1zzHFnalZBitLpWOKjdDSqJezgCaFXBZRIjfGIqT2NmkrEK6Fb5BrIZqFWULBvY3mixDeMW5zdwtzonb+Pj+2M1FXcB670XW6O0kpmkDJnP9d8JQf0HRdCsWORlbnJ6JdAG4H8FlW/GEAtxPRg3bfh2ffPIVCMS/kpuu6DODwprJnUJiuiwaEnqP71nbPFR29s745dVgiAJ01e36E6gftksol919+5YkwTtHFK4KsqL2+eZDVyfMz7fGxY7MW2+SiYAFNsr0JV+GYS3eTmAHRaVSqvTHfjO5kWRCZWlKmckiuwuz80S55Suzcfjf2jm096t6rUCgS0M6vULQUiw3jtTwGbrgMAOj+aBeA0H4ZZOy9OKk2D9Y8C7bb6h/7wyjTeJmr6629lfkBJO38sdgBtg3BSrCG9F6yH3O54grADDu/ZP9tmp5qy7Zgk9XDyS2w8ycz685gqKrvbczPQ5oixqYIqQjAqVgIMVdhZ9Vix/bP+so2jvobTZeqC9r9aPXLrysFHfkVipZCO79C0VIsNozXkDA4V2n7nb8Gd8Pt71+vt+n0bn+e/UQNrmHTgnO+6UHij9FkGf/EdepgIsy5iG0HdNBpmRMhttwKxaphEbku32AkKbHkKsohhvbKcVGVNOkxSK6xgtwk1Qf8PY85AUly+XMS5IoBL3LksufQXauEcFfxaDCPBIMW5a4wuewddIk2Nueq9LLYP/Y6l874C167ni8R9A3b/ePqmIsvrvaPBVf2GHTkVyhaigXH7adaEbf6ZPU7ev25eu/4/n1+m7lyjmxsf1qLDJsc9nM9XpkMAwb4rzF3/w1G+9TXXnDZREf+2nNWI45IMQWXFMZLYCclWWBjo2ORXIEZRP0t7D0J7u20cokzPJaySlLYBiyCKX/d84nZ8zmL6wvHcrksln494vJj+fUISstwUZU/0Sn3eEBb2vAncmX1cJctW3Iri5ANHfkVipZCO79C0VIslvaPgY5do3zhxXbN8vc91R8e9jRn6VnPEYdWSdI/y8oiYbwcBeRUn0cC7guxAdKKHb8dUD07NeFKoiC+eiT8k1QmhXxyq7uqf4T2xKirEKIstnovmJqk1vML58dchd19iGalFe5NyjWa3w/+/MSIu4nQaoFykZ/PpxBSGC8OKSYBj0PBsukOrfstTzkXKKrZ9GhwoDp26Wl/8zZuWKu3l364Um+vPb9qfPdp+zLOMoyXQqG4OqGdX6FoKRZL+3sGo8MVv91/XxW+6/xLGCU/6KlN79E99fa4W9Gf4fO9H0D/seV6mwcEcdSRrwrkrsJuutDllGz/SDzW0TKuXR2z1VV1SrExTZyzGY7Wx2gwt892rFUjlnhEcv8NXEVZuYsAyyMfc0tHYPVwvggx5uioLbdfr8lU3BdGZDHUAViYjbzHEriMbPmYZXTmz0/M2Bx1QXb+FuwesHbz0Fn83vhCtsnau/SM1dBf4w8YsPdq5anqoa0d8Q8nNq0jO7UcMfmd0z7cXfcWbyF73h/tBQBceMeF6rjl/HxdOvIrFC2Fdn6FoqVYKO3vXSBc+z8qraSj6h9845/U+//rf7iz3r70fH/e0k3nAQCDB2UnII6a1h3xU4TuGa8d7Vu6OLrZBxNdfnBXvT3YM0n1Yqv+XF1jnjWVrSDsMeuESIkjGnpH58esbLSXWUJOV4/NZTze3K5AEy64Mwf18tWRV+yxCRdUbtEIkqMIdfTOs+nVHn8wl+Hay3M18ojL9TRklK/J5tfApzzuPo93s6kecxWXLEfBClDWruXT/uZxOu/gok4D/n2nyDWM9/nz93/bTolv9m189S0P1tvPfOAF9fYP7qqO+ZmjjwEAnupviPIl6MivULQUix35z67hms/fDwD42a88DgD46HffWO9/3qP+q/X0T/nv0t86/CwA4Mkv+ZH/2ZfLdv7h0WrE7z/iR3vu2rr2wmp/58lVX8hZhPBh5vEAuKLJ2fxdyiQA6LBRpMN0LyM3ikTiDCyd8aOEG/X4aN9/loU4W5FHZl+x36zrY6eMDvn29p/0F+/YVGB/DtJL2RGaKcU4++AMydmoB3szVh/ZQwLTP1esHajau/qQV/LyHA/BGnzbXm43589/4O4pu67+Rb+9dp0fgXnKuLqM+3kIi8eIuZXv/Y5/puf+jlVmn/OKO46VR3z5xeOV3P3HvGLvvr+8ud4+csTfnH/4qm8AAL5w4hXVuev+HqWgI79C0VJo51coWgoyJn/979SVEZ0GcAnA0wurdLG4Blfntel17Ry80Bhzbc6BC+38AEBE3zTG3LbQSheEq/Xa9LquTijtVyhaCu38CkVLsR2d/+5tqHNRuFqvTa/rKsTC5/wKheK5AaX9CkVLsdDOT0R3ENEJIvoBEX1wkXXPEkR0nIj+kogeIKLvEtF7bPkhIrqXiB60vwe3u61NQERdIrqPiL5o/79arusAEX2GiL5nn91rr5Zra4KFdX4i6gL4OIA3A3gZgHcR0csWVf+MMQTwPmPMSwG8BsCv2mv5IIAvG2NuAvBl+/9OxHsAPMD+v1qu62MAvmSMeQmAV6K6xqvl2sphjFnIH4DXAvgz9v+HAHxoUfXP+do+D+B2ACcAHLVlRwGc2O62NbiWY6g6wRsBfNGWXQ3XtQ/Aw7B6Lla+46+t6d8iaf/1AB5j/5+0ZTsaRHQDgFsBfA3AdcaYJwDA/h7ZvpY1xkcBvB9hDJ6r4bpeBOA0gN+2U5p7iGg3ro5ra4RFdv5EHNidByLaA+CPAbzXGHN+u9szLYjorQBOGWO+td1tmQN6AF4F4DeNMbeicjNvD8UXsMjOfxLAcfb/MQCPL7D+mYKI+qg6/qeMMZ+1xU8R0VG7/yiAU9vVvoZ4HYC3EdEjAP4AwBuJ6Hex868LqN6/k8aYr9n/P4PqY3A1XFsjLLLzfwPATUR0IxEtAXgngC8ssP6ZgYgIwG8BeMAY8xts1xcA3GW370KlC9gxMMZ8yBhzzBhzA6rn8xfGmF/CDr8uADDGPAngMSJyC+PfBOB+XAXX1hSLXtX3FlRzyi6ATxhj/v3CKp8hiOinAXwFwHfg58a/jmre/2kALwDwKIB3GGOe3ZZGTgkiegOAXzPGvJWIDuMquC4iugXAPQCWADwE4JdRDYA7/tqaQD38FIqWQj38FIqWQju/QtFSaOdXKFoK7fwKRUuhnV+haCm08ysULYV2foWipdDOr1C0FP8fX+piJJw/QWQAAAAASUVORK5CYII=\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAP8AAAD8CAYAAAC4nHJkAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvnQurowAAIABJREFUeJztfWmwbUd13rfOcO9986CJJ72nyRESKhEEqByIXA5GiAiiAv+IMVTsYCcpUi6XCyhcNvhHKvlhih+JDeUkrlIxBAdsDBiCAg62StixgymQhJBkSWiIItBDw3t683Cnc07nx+7evfqe1be79zn3vHffXl/VrbtP771X9x5695oXGWOgUCjah865HoBCoTg30MmvULQUOvkVipZCJ79C0VLo5FcoWgqd/ApFS6GTX6FoKSaa/ER0BxE9QURPE9GHpzUohUKx8aCmTj5E1AXwJIDbARwEcB+A9xhjHpve8BQKxUahN8G5Pw3gaWPMMwBARF8A8E4A0cnf3b7N9PburX5Q9dHpn6R6//VXHqq3n3p4W7297cbq2ONHt9dtozlG2Ixv08g30ZZhvd050a3O3+XbzNmu3+4zYobGaBnOK9lroFV/DWbOn08D1u66EMY6hq5Al43LtZve+FjHxmtpYeT3g22KYyChjR+b2s+OEccCiOPl1yNde9DmH5kM1i+6488keE6cLr+2jomOlY8LADpnq2O6u1b9EE706+3hFuE5xMY7V/2Ye5nNjWterreffHhrvd2/oboRZ04tAABWjx3F8MyZ2BMKMMnkvwLAc+z3QQD/aN3O9u7F5R/6AABgtFBd4P57/Di/9Qf/td6+84rX19uv/+Pq2Lu/8DN129kr/eSlFfZwbHN3ybct3HS83p67ezcAYPHOk3Xb6MFd9fbiFYN6u7NYzfTusqc12OIfuJmvxrXwgr+NS5f7h98/5tsH26tj+YvWCV5mT3e4fWjp+pdn6QpPd+6liu7KReyjtuS/SvzaV+1HrnvW7zfsqZMnUU/e0bz8AXPjHfL9bD5zWiPbR29xfCyx8a7u9fd+/iU/yGV77XPsfgy2s+cgTO4OeycGO3y/c0erybLCnhOny69tuLV6Zp1lNlZ2PYPLV+rtHQ/OAwB2vf2Fuu3U/9xXbx9/ddUfLfmvFv+QdNg7hisXq3+f8sfe+98/VW//0/1+blzymR0AgO/971cBAA7+we8jF5PI/NLXZWwdIaL3EdH9RHT/8PSZCbpTKBTTxCQr/0EAB9jv/QCeX3uQMeYuAHcBwJU37TQfeOv/AgB8/qNvAwBc9sGn62Pv+OX31dv9b71Yb3/hb6pudvgPLbDDf7lxxMsAvTN2tb7pRN02/z9219vdX6xEC/O3l9VtS69gK9Li+Pdw1B9f7QFgy3PV7Vu6btn3/6IfyzAiAtS0WNNwpx/D1meqlWjpxsW6be7ZhXp7dW81hu5pJq4wlnm4wFfmqpPeWd/Z8j6/wvZfZtzJTkv3jL8H/NrrbbbS9o77Y1cuYXSPVnRXd/vr6p3w4+XchWOJaZn3C4/Vqp2LXKOtni6/BkfLrdpAyPUMtlX7u4wrGzIRcsi4hP6RXkATAAY7PN3OET/I09dU7acf9e9V/1J/3vyeparfB3fUbUuXelrDy/07tOvbFVt/7Ucfqtve+s/fW2+f/LpfRE+eqPr4B5+p3uvDL/tnkMIkK/99AK4jomuIaA7AuwHcPQE9hUIxQzTW9gMAEb0dwMcBdAF82hjzu+sdv33PAfOaN78fAPDCrVZ+ZHJe9yj7GrMv8O6Hqi/s/J1eIXjy7y6tt5cv8l9QJ/912Qo+YLKmW9V6TCeycqVnKXrPs5XbfvH5qto75emu2pXOrRAAMGKKqOE23m93jBZfyQJZ0q6W3VOyfOhkWEdzjC7b7p+oOlnZx2TcF/2Ktcpk51pmj6mLbHugU+Ar+0k2HvscAsUcV2qxPtx45475G7J8wD+TuYPVM1m5mK3KjIsYbPOEaw6L0ZfGy7mQ4RaZSxhauoHilt3b+SOexvI11cq+8LTn0JZf6Tm3rY9sAQCcvsFfV+cEY7zZNPQ6EKareNn3tf8NP6m3F36lOuZz3/kSAOC2tx/GDx5aiT3BAJOw/TDG/DmAP5+EhkKhODdQDz+FoqWYiO0vxfxVB8y+D1ds/0VXHQMArNxzcb3/9FWe/RoxVoysbbR/lLHcOxm7ythJZ0LkbFLALh6r2lcv82zw/EHP6nPzGTfL+Ua/6cxJgUKQbQeKM6fgEmzh1QHj7dwUOGAKwf7x6hq42THwa+D3w46ny0xrIy4iMKWkU3YG/gMc7hLYfm4G43a/kR1bTHnI4cSNEVNUBqZJK0KM+FiZaUyy+XNW3ZlOAS8ucJMeF7+k+xjs53R3jCszA6XkAnuHh/a8iOjTO83M1defrsZ6n1cOnr7RKwR3/mC+3n7dex4BADz5n24EADxyzydw+uhzWWy/rvwKRUuhk1+haCkmUvgVYwR0LDt39u8qdv/s9V7bT9v99tbHvdZ0cV/FXjn7NgB0FmW2z9m+Bzu4Fnv82O5xf+krTKvKvQVFF1Vum7dsXSfwMGRs4Zb1tdCcBeSeco4tN2w/Z5+dHZ+Pi4OP0fXLWWYO7tfgjonRrcURdr2jeeEaAXSsxpqLGDG3Yscqc1Y+EKUEC05wPcKlcfGrx6wmTlTi97vDTeP8OVijSPD82fV0T4+LJsPt/mAnngHe54O/q5zu6h523lOVG/vZy7mJx/fFxeMHvvTqauNfVGL08AHusrk+dOVXKFoKnfwKRUsxW7a/Y2oN6VmrgV/4kddcrjDN8fLFzHHHunIuPOe18suXyqx6zZZxtm5xXOvOg0ECllrQoIeRfIzu8rh2PGALJTY1okgfsSfhNPOcteWsfMcaKpKRbeyYjnSP1vTh2F8TWRLc/eDnBFYEPl63ya83oul2LDhn9QO6NL4/Gh3ppCvG/QZOPIvr31se9CTdj9CyNG6d4AFL3OHLaf57zCWbX0+HOSI5MZSLK2AOXXMnfR/9N1XRfuYbF1XjO5nxUrg+s49UKBQXFGa88gOYr75qux+oVvyVn/MBOHMP+9DaJRZ8MveSDXS5jLkCs1WVf5nrsFPm8hkoqOyq12PKmgFT0gT2Y3tIGJPOr0fgIrh76NbxMaRWVT5GviLxsOJa0bg6PtaKGNscunFH3IoL6NaKObY/pkiUFI0U0UVJ3IlElysUJY4F8Pcs4KTYu+KukdOiCHcSxPa73ZxuwDFaukPOoo37Q4wiPhT8XVmw4cxL+5mL80tyeDd9t1KcL95c3YThN/L9dnTlVyhaCp38CkVLMVO2v7NI2P5Yxe4v/uypqu37ntVfvpzF1fMY7K0VKxO4kkbg3Da5HwBXrNTsqNAWA2f1JNaUt/FY+kDJJg09lk5LUFpJ45XY0rV0ayVd5BqGzAXVsfMx0aQeamTcUoqrGKsfEsS6/To3a+5TENxbyb03du8sWx6IWVyEiIlSEl2e0chqJfn5XNHoxEmuUAwi+dj7vmwVfp3T/uDYs659ERbswAqWc135FYqWQie/QtFSzJTtH80bnP4pmwDj6cqFcXk/S/3EEzQwDbxjiYLEkREX1Lo9yMLK9tfabz4wv5li9YIMswKt4FhOy52WFW+ViYwsuo5djGXRDVyfU2OTrkFw0wXg72lMtEmB03XZjIcJMSoDPlOwb+vEsgJL1yBYUvh5hmn4oxmTxYGNNwVJSHcx91+WtmzZWsD6NgmNGIkaga78CkVLoZNfoWgpkmw/EX0awJ0ADhljbrJtewH8KYCrATwL4F3GmGNJWkNC3+ZNX91T8UxzLA/aCnPi6R9mWWW3jUdixTux/xnLxTPuSs4WHYFNjiJgR+353MlEYvX5uCIIovocXZagOMi171jXDLZ/qnD9cfqCSAX4eyO1jdFIiBNivw1RP5+IqCeON2bR4OKefQd41uEsS8c6iIlq3HLjkoi4fH9BYZQEclb+/wbgjjVtHwZwrzHmOgD32t8KhWITIbnyG2P+hoiuXtP8TgBvstufBfDXAH47SatnsHpx9Yns2sylqzv9563Hc6kze7np26CIWMZWpgiqU0LNjysMOV2SVh4gqaDiX/s6wCZyF2N5ACT6kh0+KB2WohWDc1Hm3AkfrxRsk7ofQVqr9emKSsC1dIW2UMlq22KrbkzJug7d2PmjRGBPyb3v8JRfrvoPd0uPZHJ2fhyBzwh7hzlHOLyyyhq89dEqO3AnL3FvdWz2kSEuM8a8AAD2/6WJ4xUKxXmGDVf4BeW6Tmm5LoXifEFTO/9LRLTPGPMCEe0DcCh2IC/XNX/lAUM23tmVj1pm6bbMCk/d5Gn0bYzygIkIPP5Zts1H3HvtJreHxlJRSaxpUiEUuJUK7Tk2YzdeTldQNEbZ6ASbG6XbGd8fQHA7jrmrOromMq4mdGNuxSKrnxJNchR+DegGz5+7DS8JolzkGlzRjsA9mL2XvJjHwpMVu3/2mkoW4C7HKTRd+e8G8F67/V4AX2tIR6FQnCMkJz8R/QmA7wC4nogOEtG/BvAxALcT0VMAbre/FQrFJkKOtv89kV23lXZGI59EY9kmKugd9rxxGJ3l2ZfBVqF6bCR7a80+cdZXSj7B3XQTUX1RC4BkDRC0xcG4MvqoWXUhkQaAkN1PIMXKR2vpiQfnjyVXhCiim+PXUPt5sN2CyCOKOzG6kYzLkmgS0A2Kp1T/Y++EKC7wJDU8mpRtL19i6za69F3Djdf2KxSKTY6ZBvYY8gqJ3svVp5AnQQzyzbOvpvPKG/La6JG8/fX5EcWLayfhq7z2PAlhfn3hnNiKIwWJcLrSeZG8/qLCLwIxnp+vPtLKnFLMJQKhqh0JWgm6or098n40yZWQxUlJY8/xcajbxhOSRoPThPFypbbzcl07hnlblq5O7TVlDz+FQnEBQie/QtFSzDZ7L1CzNINdFa/MSxoNtnJ7O3PJtbH93MUxWkk2EQTi7P+8girPtS7aYSPKIVFEKLEZM0g2bEnEANL59aXgo2AsMdt7imOUFHOxc4QAmpTCLnCdFjLyBqekdLQJV+DovZPuR8Qfg4uLYo5/nm3YpTWL5aHgxzoF9ymeCdp3zIPhlg5UivO5F/tjdFLQlV+haCl08isULcXsy3VZFt6xLqu7Iy67PNWUyyrLNZmxDLKOrYtpt/vjGXcDEaJJzHhDVj8VPRdlg0s+2dL1pMYbG1cqrj7hq1BCN2DrXWRiLDeAIFrEtPlN039J50uiVFhAhu0XLEP82MAV2Iq3A1Y6jmeu5pr/zsmKsBNBiqIO8w9VKBQXEnTyKxQtxWzZfkN1MgNXiINShS04MjTTEusqsZ5F7F+qOEaOq+i5plsSEdcg0+xaNKIbOzZlZZh0vKnnUHBetK5j4r0LzuutbxmQ+q0jV5XtVygUKczezm9Ru50yd8dU5dys2PBEqqlUgEVR0k2XqiqhfAro5nyZmyS0LEg8OjHdgrRnyb5mTVdSHjZ1UU7RTSktY5WVXXky5vfSO8nmxo515kaBwlpXfoWipdDJr1C0FDNl+2nk03et7Bl37+XpiXonfbuz73N7fLSaqpASKhm1V+DiKtGdpggBlEXtpSCJJjMpKdZANInF2BddQ90Y6SSRXqzITyNBN/DTkDIy8zEKFY555CrP5MvHuzY/hkb1KRSKJHTyKxQtRU65rgMA/gjAK1AxJ3cZYz7RpGSX6QCrVlPZWbHJDbYybT9j9XnGXcf2cxfHaFSf2x9J/OBYS04/cPUVEoMEdvFEkokmUXJA3B3Vdybsz9BCi+WpmrL1giWlpBBH0p055pKb8ntIiRaCS260r1QClhIxRsrIHLlfQQSnrR9Hhr/vrErvMd++7KL6fmKj+qZcpXcA4EPGmFcBeAOAXyeiG6EluxSKTY2cBJ4vAHDVeU4R0eMArkDDkl3ua+eCeHhAQ5CnfNl/weoYfCEdF6dZDTjeJ+BXER73LK72MSSCRIoURiX29qbKvxJbda5HW0GMvnh+rL0kN0CKLm8rCexpyMnUdQpiPh8Jbi0s+uleUlk5ONjOuOVD1Yq/utOM959Akcxva/a9FsB3oSW7FIpNjezJT0TbAfwZgA8YY04WnOfLdZ0+3WSMCoViA5Bl5yeiPqqJ/3ljzFdsc1bJrrFyXY6Fn3dsCstwyiqYBgo/q9xLptviiGVZdWxRU9tuQnlUpIgS3ENjdMXzclj2lL29iZttiQKsQOFXUh+hsYtySgxKKVYjdI2wX0rDFvUJEUSaIE0YU+RxkXWwpyLYP1J1EA0GEpBTsYcAfArA48aY32O7tGSXQrGJkbPy3wrglwE8QkQ/sG2/g6pE1xdt+a4fA/iFjRmiQqHYCORo+/8P4kxScckup22vK5GylESBayQveGDLdA23MWtAKg9AzCZca/t5X2y/pO2NaY6F/clUUzGWO2VFkPZnpNPaCG1/tGBGkxRovN+SseTSRIbrdOQ5pMSYFF0pDVtUlONin0tbx2z7PJu1y3wNAHOHq5d3dc9orM8U1MNPoWgpdPIrFC3FbKP6hkD/lI3quzTUUgJhkoIuK9Ax3OpcghOsPpBmIVMJOErQIKtsEXLcWRsgmmoq0204K9qthJXPFU0aJvCY5r2TCqIAkO9diUu1ENVHrOJuENHKrF5BtF8hdOVXKFqK2Vbp7QKrO6ovlasnPgxyk/svXeDKa+2ajRM4lqCElhQbXvIhntSe3vCj33i8IjG23dQnYCMUfjGUjEs6NqUcLElFlhhjtOSYlJbOxf5P086vUCguTOjkVyhainOWvdex8Dz+OFo51yk7GqaEGnF766qjyY5NKWYmZY1j2CgX1Sa0mqLpNXBssJhThBKXbMHmPw1FspsbXMEdzI2V6WgwdeVXKFoKnfwKRUsxe7a/E6bkGmz3RvI+S+PFK5Q6Fj5qT+dw2ZJ4ui2hWmrjtFap6Lum2u8mdBuKQdHrbaJ1n2YEYAlSqcSaihCpqL4IXcnluiiqj9F17+uQFe3gmXxNh7nE960PzGJ3rM8UdOVXKFoKnfwKRUsxc7bfO+xYZ59F2VWRsy9O6zlc4FFOXNWKse1OLGrPiRApN01OtyT6Loe9Xks/g27y/Fh7iUtuLt0cy0GJCNHEujGp67TUP9C8aEd9EtuUrEwZUXeSJYyPkee6nDtWEVzdpVF9CoUiEzMv19W1LrwuiId4lt6IMsSt+Mn8+gyp0k7JPPkxJGzkRUEkJbb3pvtTCq4CW3ajcl5N3WglTJHjiHFoyfJlJUpLzs0Jiubku8SCeYjVrOA+MKs2k2/NCavCT6FQpKCTX6FoKWYb1UfAyLEylk3hrPxwh2zXJFfhNJVuizclYuxjFWFTbFPjarqJlFBN3UpTfRW5nabY7xRLXSKalNAV2qKpxDLHksyMnDGGJIstZJDOUsgJ7r1cBODbrkrvYKd7sTPojw9PBhEtENH3iOghInqUiP6Dbd9LRPcQ0VP2/578bhUKxblGDtu/DODNxpjXALgZwB1E9AZorT6FYlMjJ3uvAeBK7fTtn0GTWn3kk3S4IgQjZtt39fvWwhX44Mk+ktr+CHtWV+kVqu0CEXYwJiI0Tf8l0G2kzc+I6pOq9DZGE+13U21/QoOftMGnXKdjIkTTKr3CNQSiZ0LbH7gCk2tj7/ucXKV35TKbDs9m8Q1qWCaQpfAjoq7N2X8IwD3GmOxafVquS6E4P5Gl8DPGDAHcTES7AXyViG7K7SAo13XVAYM13ksj9kWjVVaui6Xxciu+5Kk33qEjJh/rvuzJ1b4EwsoyRjeVEmpS232MroRJ6U473dak8frSeAq8J4vSmhVwPSV0xRz/rCQdV4C7VHiAT4dXl7+bpsIvGKAxx1Gx93fA1uoDgPVq9SkUivMTOdr+S+yKDyLaAuAtAH4IrdWnUGxq5LD9+wB8loi6qD4WXzTGfJ2IvoPSWn2Gsfvza9wSscadUVLC5LCFmUqnjQp0SdYTKO1jvf1N3WwnpZvT70bQLbmG85VuiZI2EiS0ro9DwTuRo+1/GMBrhfYjaFCrT6FQnB9Q916FoqWYbTy/8Wy+Y/FNl2kugzJEzNXXVS1tWAqrqEprkxRW04hcaxKJlzPW1Bib0C2Jrss5r4TupP1KJdZS7t05zyxFVzonQte9r9ziJZWvA3yE32i+3JFDV36FoqXQya9QtBSzZ/stWzTcXrEpc0c8L79yseeZeid8u7MMTMzqr6XRABLdLBZyUmeYEhFio+g2ccLZ6HsQoxvb76IcI8teNNozQVfMuMwPTfQrOqIxMbieA2uO7dh+68I0G+Xko1AoLhzMduXv+Kq83dM28eBuv5T2WN7+EXNtdMk+O9wnIJG3P7nap/K+CzRjdIvy/kfoFq3AKcVcifKwhG5uX7x9GsrQ3PNzaAgoCuyJjCeVtz+paBbocmU4z9XfO8PSeO21gT3HNiiwR6FQXHjQya9QtBSzZfsJgHPhtSmKeIz+kCk1pLz9WVF9LhY6kcarseKvQESYOLouRjdFv6mNPJduE1+IUrqp/U1FEwHRHA1N3HPZ/kAcTIkQQjsXR+ogfwCDrcz+f8pG9Tm/GVX4KRSKFHTyKxQtxezt/CshX8LTeHVPdcT2Ub/6n5WAQ9D2F7lcptjJVDXdpm62m0HbX3INuePKoSuh6fWm9qeq9BbQbartd+9VkNCGle4iJgKMbOmuWAq89aArv0LRUujkVyhaitmz/c4xxmonObvC2Rzu1OAq+fLqpKkqvSlWP5oQIaX9ThWLmKa2v4k7bYzuNF2Bm4obk4ox0xA3JOTkYMykG2ThFTJE5ySRqav0DllnjNUfbveD7NsqvXXRjoLlXFd+haKlmP3KbxUXw53V16t3irn0bvPLape5MLr45cZ5+9eMAchY7QuQFUS0ET4BJbkBcvbn0p0GF5GiOykn1NAXIWknF2L4Aab8FbLwArL7b8n18ky+QWz/ljVBbwX3KHvlt7n7HySir9vfWq5LodjEKGH73w/gcfZby3UpFJsYuRV79gP4ZwA+yZrfiapMF+z/n8/pbTRvMJo36Cx30FnuYNQ39R+tUv0XtA8INCCYDuq/KIz9I/ZnhHaprSE7a8j/BXSbosk1cEjnbRRi9zF1DdJ4S2hNOt7UvY3dO75/5P/q91JoM50pXAMbY/C+uWHZObIRbP/HAfwWQim2QbmuM/kjUygUG4qcoh13AjhkjHmgSQfGmLuMMbcYY27pbt/WhIRCodgA5Gj7bwXwDiJ6O4AFADuJ6HOw5bqMMS9kl+siwLgqvVZjyRMWcFYo0KCXsEiSZjihlW2ccbWJW2oOmti9IewvobVRyNG6b4S2PwJnh5+KtaeJL0LG/XBj6/KI160RS5i1kPVOV200TbbfGPMRY8x+Y8zVAN4N4FvGmF+ClutSKDY1JrHzfwwTlOtyK/twB6877pfjwTbBm68kSISv9tzeKtRJL0q6KXAnUc6kqafdFFeRIqToNuEicp5Zk+tJ0Y3sr5/1hIlcx/pN0U1xL2y8Ln9FkJ9/Ua5g7ebG0LaVcMlFk98Y89eoqvRquS6FYpND3XsVipZi5u69ji0a2nh9np8/YPVZIIMLcAjYnViWUhrfn5X+a835sbakcrAh3UZKuqbiRlO6JZimaNIksCdFi6MksIeD95vIEJ3K6y+58gYx+izQzfAyXjYHRi0iFNwLXfkVipZCJ79C0VLMnO2vNZk2GmnI0nV1eIovtunY/ayCBJIdV2LJcmzdUnTVpBVdp8G6lmjdU3SnKW5ImKa4UXK9Bc80maMhZwzC2MUM0iVpvPq+jUcIdhb9iTzdXSl05VcoWgqd/ApFSzHzoh0ufZfp2yq9R1mV3os8n+Rq+QEAjQTWJqHpjibYMGv+r0dL2C/SjTl2FDgMlTiqrEszdt40j81wVBmjOQld6fxJHYJiLHeJCCCcx88xJdp+Bpf+i4vBbt6s3XaJPeqkHgXQlV+haClmXqV3ZL9QriLvyl55tQ/smjb4h+cuT6bpmmGV3saY1I7fVAE2qXIwR1k6qX/BRl2DcFyytFtMsRdTBLtDU1xiivvgc4Cl8eqf8J2tXGyr9B7tjfWZgq78CkVLoZNfoWgpZu/euyaenuff57HIvGKvU2pw994o+0TCfo46JzobViqqL8Xm5riHSqztNGz+EnLY8vUwzei7HGwE3dS9jSn2JhTlsjI5iyeO0wgyAbOovtWdnlj3ZDWF3XzRKr0KhSIJnfwKRUsxW7YfPkLP2Sp5SaIR02jSqtDe1LbLIbn/5py33rEJrW8W3Sb7c8a6EXSb9rsR96DpeU1pJY4tSgwTgVTgg9v2uajM20v70ZVfoWgpdPIrFC1FFttPRM8COIVKVz8wxtxCRHsB/CmAqwE8C+Bdxphj6xIaeZfF0fy4Gprn2pMimqaacRcZ+3PdVXPcXVOOO01y9OWwlSWOOZP22/Te5NItuXex/U2sPTFaHLl0M67BWb2CKr+rfptbvbpL1r1328Ym8/g5Y8zNxphb7G8t16VQbGJMovB7J4A32e3Pokrs+dvrnkH+azhaqL5U84f8cr58mV/6e8f80OoqvUGgA6Mrxf2kYvAb2thFusLKAhRyJyXYiHRbOXQnzUkwTbpN9wsBOI0hZIgO3kvJFyXjfrkVP8jbzzhlnuLOub7XCvKCdyr3FhgAf0lEDxDR+2xbVrkuhUJxfiJ35b/VGPM8EV0K4B4i+mFuB/Zj8T4A6O7RKt4KxfmCrMlvjHne/j9ERF8F8NPILNdljLkLwF0AMH/lAeOUFb3jFbu/fInnjXsn/HACV17ns1gQTRYoS3jRDonXKVBaJV1BS8QNDonuLCP1cugWQHJRLYmenKqSNoGkS27MJZuf5943QXEHRETEYBB+072vQ0EpDoRZfV3Rmzq11zTZfiLaRkQ73DaAtwL4e2i5LoViUyNn5b8MwFeJyB3/x8aYbxLRfSgt16VQKM4bJCe/MeYZAK8R2svLdZGp03eZQcV08Iqj3L2XV+/tLNuKvmtdGcU+7DkpVr/E/lwQ1Zes1RZry0jyMDaGHBfVEpa5SRKRCMSaeCUutZNq+5u40wJlUX2x5+foCn4pqSQ0QDqNF0/ZRfaY2uV3A7T9CoXiAsNsA3tGhI71SHJfqsEur/DrH/esjUcrAAATYklEQVSfyuGCP01c8RPKoY1S+ImrfFMPv6YebU0UfjkrQgndSffPgq60P7eabqxffmgiTVcyb3+k39oLdmE01hbrt+YGCpZzXfkVipZCJ79C0VLMPJ6/49j9HRVPxFn9QST3eNeV+OIujrHSXfZzxgMhoi6XuWAsWeMqvSm6UlqpScWC2LHTFDdiaKLwS9FtWk2XI6GwLXHJlvw4uG1fCswJ/AgidJ2PS1C+ju2vK/IC6Nss2LzCdS505VcoWgqd/ApFSzFztt90XBSSzcjLNfmsSAG3/w+3WxdG5tYYZZ8kt1LB5bKp22nj3AFN6DZFU1fgXLo5NvRJ6UrjjWXcbXK9MTfdkucgiIOxGPxa9Cwo2sF9AnjRjh4rbjOwIgBP7ZULXfkVipZCJ79C0VLMnO1fFzxJAc/ka1maxqyxxF7H2C+JXZxiBFkWctng2Fgm1faXoIS9To13GrSmidzkJvAserTKb+reJq6X16nkEa+ufdRgJuvKr1C0FLNd+bsGQ2vfd2m6eJBC4IYrKORI+NJmIfFVTdZkb7oibRQKVqSpKvwEJNOlFRGLnD/F/AW1wpfb9vnzT71XG81lAD4ZKF/tmauvK18HAIOd1Xyae7maMFH/FwG68isULYVOfoWipZgt2z8kdE+FBnqz4PmUzhE/nAGrROoqlI7mWRtjiVKRepLLZlY11SYpoTgbXJISqiQSb1KFXwwpZajQ7zQr3CajJ0sQE+tcdegpuGlzunUMfiyCNPUc2LHOrXfIXHYDv5ctvr1n59PqXpsno2BG68qvULQUOvkVipYit1zXbgCfBHATKgbmXwF4AqXluoCa7RnstlrKw34Iq7s9OxO4K1pemlf0jbJqjq3n7JcQ1RdlVxMpoVLRd0V0S9CUlZcQExck0aSEFkdKNFnvHD4WjmlE3wmus0VRmxGxri4px11yO+vv5+D9Dq1mP0hxx93geR/Ov8AV+CgQkXJX/k8A+KYx5gZU+fweh5brUig2NZIrPxHtBPCzAH4FAIwxKwBWiKi8XBf86t09XX13VvcwJR4L3OHefs6jia+6Udt8bmBP7LOXULwllUNT9FKL+jVM0e49sU9AbH+TIKCY0tKhJLBHog/ISs2mgT2J85Lp4yDv7wj5Kzi6i+NJb43z+ptyGq9rARwG8BkiepCIPmnz92u5LoViEyNn8vcAvA7AHxpjXgvgDApYfCJ6HxHdT0T3j06faThMhUIxbeQo/A4COGiM+a79/WVUk7+8XNeBA8axr0Prlrjwk3597NLlPgCal+6qFX4D/63igT8S2xcTC8TAHo4mCj+OEqVUgqVOiiYlabxiEK4tyfo2zFlQokgU73NT0UZQ6AXXWKLw40PgCj/LqvP3UlTuZYhUdRqvZfa+szoW3M7v+iUnMk9T4WeMeRHAc0R0vW26DcBj0HJdCsWmRq4/0G8A+DwRzQF4BsCvovpwaLkuhWKTIrdK7w8A3CLsKivX1fFRfQvPV+z+0oEVP5iXvQggVuntRFh9Bscu8tJfQdkjKf1XgZutWIIpYtsvsgwI4oKJiRBSAYgpRMQlswanILHXOa7TmVaVabhk16w4HyvTyo/6rH04fmwQY8+1+e5945alRBXfAFzb78rTRcrX9U75g1f3hlmwKSXmyV0qFIo2QSe/QtFSzL5Wn01EsLLLZuQ9zbX67FAeuXSiYml4sYJYVJ9jB4P9nEWUrAEN3U5FFjJGN+XIknIbbqjBn9gVWKKbcg9GhlWlBJJjVopuQgwKxCj+ChYkw5DeG34+FyFq997YchuICOMPmCfwWN3tO6nnxoJZn74AXfkVipZitis/YXwFEuybAOpqvgAwmi8I7Ml0343ub5LGKcf9d9L0UAn/g5nSzfEvyKV1jugWuQfn0HXvFVcIlyh/E+9VMDeY/b9ub5A2TVd+haKl0MmvULQUM2X7aQR0Fyv+ZrDDlhliirmYkqa23XPWJ2a7FxQ6TaulZrurpmLa+TGxeO6SVFOp6DoOd2xJnoEmUXKIjDeHpZauJ+H3UDRGgW7WvZ2UruQTEos2LanSzH0NrFuvy49RYujXlV+haCl08isULcVM2X7TAQbbrD3Sui7OH/F8EI/q6x9l6b1sBGDvpD+W2/y5FUBK05SKEIum/GpSeCLlkislvABEG3Y08iwVfZeyazeNFhT2p+gWFVrh9056JiXuzDG6kghRAH493IW8a61TQ1Zcg4umzuYfS/ARpP+yrrw8aUcQyce0/XXRjiOuaEe+uUJXfoWipdDJr1C0FDOv0uvYprlDFZsSjepjLo606op2sLaBbCWoLQPsWO4w5MSNlBsm0DDKLZFbsCgBR0o0yYlmk5yaUpFlTaLk1hzrxitlTh6jm3C5deMtqgtYIprE+hWemZRrD/D59mJsd507MnI/uRXAufIOdzDRdplfEBM3rGiwat3lo1YsAbryKxQtxYzde01dbXSFqu9O74hfdvlXdbh9PLCHK1hiK5L7cnfYl1IKlOBf4GSW1Zj/geDSGaxuKY4hEhQjnZ8s/RWjmzuWhnSD1Vgab0FgTypwJ3qfRWLyeZLCL2Z7T0GiG5Tw4n4p1p8ltvLzd9CV6eqc5ZpKTyvIdWHHMH/YKvwYnRR05VcoWgqd/ApFS5FTtON6VGW5HK4F8O8A/BFKy3WNfJXe7lLVtLLfK/z6L8z5Q3nmUrtpgqg/2b3XsV8pJV5gr+2tf2yMZTZjG2m6WYo3wW1YrP6bY/euB8C2efGTErquKSedVqJfCSm6Td2dpTRd0bj9WL4F6VCmdHYVpLm7Ot9fs/uRcQcKP1umy9nwgTCen8sO/ePV9pKdR0FW6wRysvc+YYy52RhzM4DXAzgL4KvQcl0KxaZGKdt/G4D/a4z5EYB3oirTBfv/56c5MIVCsbEo1fa/G8Cf2O2gXBcRJct10ciz+8Nrq42tf7+l3n/2Ks+fdRa5QdXaUCOpuzicRSBwjeTpv6wVIGnb50i4ncYiuYroCixtoDlOiRAckgjQVDQREEspJtHNqleX8EVwdIsyMnOWOiLWrR0r4KM+Y3RT2vygXyE5jenIbDl/twe7bEbeE8ydndftYzSGc9V5dNYeO0rIKgzZK7/N2f8OAF/Kpl6dV5frGp7Rcl0KxfmCkpX/bQC+b4x5yf4uL9d11X6zcokNRHh2AQCweIVf8oh9zV3cP+ATFnbP+C9hrFyX+wJHvQFTEFa1lAdYNG9/U29AwRusJOFoMjdAijuJ+QRIgT8Rr7tGATQRr7yai+CrfY7noDRGQfGWWu0DUrKjXd3HaIF5350ZD8yJJp7lQUJWuce9XHlC2/lDfpDLl1UP0AXCUYFfRcmjeQ88yw9ouS6FYlMja/IT0VYAtwP4Cmv+GIDbiegpu+9j0x+eQqHYKOSW6zoL4KI1bUdQWK6LhoT+0YplWbnIKjWOehZmdY/n9UIWsiSlqtSv3xbjqhOln1Jls2IsdeOc9Qn2OpU/P5kbICFuTCXHv4SCFFmSyJPlIyHRl/qN+RSUBHFJeRO4bZ8Hp43Wd++V6HLRtceUf8sX+xfOsftu7vCyXimoh59C0VLo5FcoWorZRvXNj4BrK3Pflke3AwAWr/LuvfPPe+P78qXMtfHM+DcqqnUV9ndGwQ93lj8/leIqFiHmLAvMnltS7inoQogTdy6jwBotcQlb7tjRkgq3BdGIkosr4K8hKIUVs1hshJgRyzBsxb3Az4Np+6UM0cH9KhEBBY/cTuw58KE7Pw/2zPm4eDq7wRXLAIAdD1XWs0PL+TdTV36FoqXQya9QtBSzZfsXO+g+VrH7gy3WZbefVon3TlffqGUeAXjY823c9bFmmbiGn7PljpViogB3oBCjBSNJN+qUYDxKjlkmRgvjbsUxFjJgPefHHUI43DUGFosIt1dHOUquqBG6ycIUq+tfI8CuMyd7rxNNEnSLtPKxdFvW1Zu7j3MHm0C8crMjcj+4yDPcYjNMn/KDHG3z19C17/CQtcWcz7p2vLwab2AVe4V/WXZ+v2L38SYbUPvVfLlTV36FoqWY7crf8XXEd/3DIwCA4Te8+8CxW/wXrXvcD21lr/UJOORX+8FWbkP1XTjbao8pCQc8sMeuzHxl4QpFbltN5qkXlGkBXba6iNVUY66xwoogVTCOuThz1IqoyCpTRNfZ2/lKyVf7IDVWdUwyAIfRDSvR8ptuz+9F+pUCcNhz4jnvJddZHmwj1YCI5QYIUs2d7o6NkXOJ7tpi5ek4x7l6aTUP5p/z+S1cgk4A6LC5cfLVFTe88P09Vf9nE/7JDLryKxQthU5+haKlmH25Lssqrf7FxQCApZ87Ve/v/GSbH9gZxgbZTe7+y90dJXYxUAJydvG0Y/v8uIYLnO1jA65dNn0Ttw9L/gWcHQ1KLDmbcSJmHWAsIqs3QDzyrC+IEHzYkvsuZ8l5FOTSuGgi3YOAPr8fEUWiY2+L8sjHyqYJY02JPIHfBDsviIt3+3nKOM62C34PUrotwN8Hfj7f79p5TgIuYgy3MDu+FW8D0faS5Xp77jGfA2OwUh27ct1i1c9CWoHuoCu/QtFS6ORXKFqK2Wr7R17b+spffAIA8OJ//Kl699F/6UWA4Q921dvdyyuX4G33ba/bTl3v+eAOT/LhtNuM+3GaWMCLHfyz1zvZGd8Pz5YFLqqMJZaizbjPAE/DNLD23ViEIhdDnJWAs9Q8Wqsu58TFlYjrrGM3w8QSjC53ybXsr1TkJBgrG5do0WA0clh1N3aebq1/nN07V7YqUhBFEnMCC8yWccsOZ79Hkcq69X3kWnm2f3CRfwcXDlaa+aVXeNmFj2Gww1qsjvgHzZ8f9/MY7K1o9FilavrxQr29fLEnvOWqas7M/eVOAMBLp/PXc135FYqWQie/QtFSzJTt375jEW/8J48CAA7fVrE2v//of673/9q/f3+9ffgfez6o/0xlBTh1Hcvue0ZWI88dtyzxDT5ZaP8Rb0VYfGVFY+Hped92DasUzGoHOtYwcA/l7pmW/TZzMrvK3TMdexxooZnjzWCB5TK0jhqjrb5t4aAf19KB6t70jvnHF6tjWEcILnAVP9Msn/IDqnMlctdXsaotd7CR88wtvFiNbWmff469E/J4nbjQ59FqOznfXv3rMkeYAXN6ocVxUYrT7zFW2Dl8cdFm/mWWKINFk7roOU7LZdYFgC0/8k44i1dW1xm8l/yRrIxbmYLxXuLfwb3frt5Nc+fRuu3Mw3s9WV7T8qFKPD7+quq6hl46SEJXfoWipZjpyr/63Dxe+uDVAIAPPlxVAPu3H/WrvfGLMX7j1nvr7a995C0AgJP/5mTdtvw9/yXkCpDVV1b2zq33+9X+1A1MOXjM2kX3MCXPqfVvAy+bNMdWiVW7OgUBHlwJt8rs9E55yFbdgO5hT9dlOO6/7Me1tM9zPc71OZayieeL71plJudC5o7wvjzdnqU7ipR8cgo7rjzkKyHPKrt0heVOmIIrqC4r0Q+uwT+frc9Wz2zpxkV/Dc/4JW5l97iSlvscBLZ9O8QtL7B7e92Sp/ucfwkH28b9KXjl3JWdTGG3UN3HPd/zdHe86/l6++g3rgAAnLzJr/D8vVt40l/P4lur93zlyT1128U/9H295Te/XW//1UdvBQDs+rUfAwCObfH0U9CVX6FoKXTyKxQtBRmzPis21c6IDgM4A+DlmXU6W1yMC/Pa9Lo2D64yxlySc+BMJz8AENH9xphbZtrpjHChXpte14UJZfsVipZCJ79C0VKci8l/1znoc1a4UK9Nr+sCxMxlfoVCcX5A2X6FoqWY6eQnojuI6AkiepqIPjzLvqcJIjpARH9FRI8T0aNE9H7bvpeI7iGip+z/PSla5yOIqEtEDxLR1+3vC+W6dhPRl4noh/bZvfFCubYmmNnkJ6IugP8C4G0AbgTwHiK6cVb9TxkDAB8yxrwKwBsA/Lq9lg8DuNcYcx2Ae+3vzYj3A3ic/b5QrusTAL5pjLkBwGtQXeOFcm3lMMbM5A/AGwH8Bfv9EQAfmVX/G3xtXwNwO4AnAOyzbfsAPHGux9bgWvajmgRvBvB123YhXNdOAP8PVs/F2jf9tTX9myXbfwWA59jvg7ZtU4OIrgbwWgDfBXCZMeYFALD/Lz13I2uMjwP4LYR5cy6E67oWwGEAn7EizSeJaBsujGtrhFlOfil/1aY2NRDRdgB/BuADxpiTqePPdxDRnQAOGWMeONdj2QD0ALwOwB8aY16Lys28PSy+gFlO/oMADrDf+wE8Hzn2vAcR9VFN/M8bY75im18ion12/z4Ah87V+BriVgDvIKJnAXwBwJuJ6HPY/NcFVO/fQWPMd+3vL6P6GFwI19YIs5z89wG4joiuIaI5AO8GcPcM+58aiIgAfArA48aY32O77gbwXrv9XlS6gE0DY8xHjDH7jTFXo3o+3zLG/BI2+XUBgDHmRQDPEdH1tuk2AI/hAri2pph1VN/bUcmUXQCfNsb87sw6nyKI6GcA/C2AR+Bl499BJfd/EcCVAH4M4BeMMUdFIuc5iOhNAH7TGHMnEV2EC+C6iOhmAJ8EMAfgGQC/imoB3PTX1gTq4adQtBTq4adQtBQ6+RWKlkInv0LRUujkVyhaCp38CkVLoZNfoWgpdPIrFC2FTn6FoqX4/6QLT8V2RXLWAAAAAElFTkSuQmCC\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAP8AAAD8CAYAAAC4nHJkAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvnQurowAAIABJREFUeJztfWuwZUd13rfO4965837oNUgzSIAKIyAIWwVScChewljBiEoVLijjIjgu/hAXxKR4JZVKUmXCD4pAnBSJIkOgTGxkDAEEliFgl0MVYIlHACEUCUmgkUbMSBrN6z7POZ0fu3v36ntW3+7e59wzc2evr+rW3Wfv3qt7P3r316vXg4wxUCgU7UPnXDdAoVCcG2jnVyhaCu38CkVLoZ1foWgptPMrFC2Fdn6FoqXQzq9QtBQTdX4ieg0R3UtE9xPRe6fVKIVCsfmgpkY+RNQF8P8A3AjgCIA7AbzJGPOT6TVPoVBsFnoTnPsiAPcbYx4AACL6cwA3A4h2/u7OHaa3fz8AYP7hswCAy59/tj7+6N276u3Dzz1Zb99/7DIAwGiOCRuxbcZfyO43Xf9R66yQ3941AAAMz/pLN332ARz5srCbNIRYllarAmYucj7/rrLd4nF+DWsYlzsUBEgygfDedG1Rfg38vI7Q9phcV5TfIlaX4TySbGHe7q4gKyKXw7U9uB9rrLDEX2P33u4P2j3Pnil7V+pnPZIbFsiwZbtnfdnRLlZgpWokfy+DB8GrsKfRnD9//her9fZV1/i+8eA9+wAAVz+n2vfQw2t4/EnpZRnHJJ3/cgAPs99HALx4w8r278fT3vVOAMCz/sW3AQAf+NLf18f/9XNfVm//5y9/td5+/R//IQDg7GH/BneX/BMfsZeiu1Rd99p+X3bn/f4yF15+HABw5tsX1/uWnjbw5y96ue5l7rEHunrZmpf1YPU1WrrKPxha8m94Z5nJci8Se/Yd9gIPd7AHfaySsXyFr6tz2st1/WrEP0TsefOP3dC+gL1T8v3i9fZOVXUMt/l9XC4Nqu0RO95ncgc7vNzRvK33pG/3kHUG9+HkcoOPMNvsn66Or7D70TvWF6/ByeIdk9+nrn0m7j0BgNWrluvtuQe31dsr9ll3FtlXi3/LmIzBweod2PctP0KdeYUf2MzPdgAI38vOkvB+sP3dw/78Z779kXr701/5Ur395hveAAD48l9/GQDwot/gXXJjTDLnT41lVSGitxHRXUR01/DMWeEUhUJxLjDJnP8GAP/WGPMb9vf7AMAY8x9i52x72iFz5e9Xo/jgeWeq/8cX6uM3vviH9fbP3vOcevvku6qy5gsH6n1PvoSNtif8KOCooRs9AWDPi4/V28t3XFLJvMaP9v2n2Og0P34/OFXjjGN4adWG+YfmvfzLvNy5J7zcwXZBLvva9094uWtPXwEA9B7xcvlo7RiJG+UAwPTGjwNAd7kqM9jn29U74ZnQcBsbcdzAaiJy7TYf7VcP+JGMj+Y9y6AGbFTmrIpfjxtGOstsJN3N5K51wvYBGLHn1DvDnom9Ht7u4Lg9zxzw7888G+2XD/pK+k/2Apnr0eUMa25c7txDXu7gWUtVXXf79335Yn9v+Dvmnutl3/Z1feZDH6q33/o7f1Bvv/GWvwIA/MXvvgoA8O27/xtOnX00i/ZPMvLfCeBqIrqKiOYAvBHAFyeQp1AoZojGIz8AENFNAD6CSpXzcWPMH21UfseBQ+a5N1Vz/uv/8C4AwB23v6g+fvgOPy049B9/Vm//4BPPBwCMbjpR7ztjFR1AOAo48C//NsYClg5WIyAf7dfYPL73uGcR9VyRfUfdXBYAFo5UI8PSlWy0OO5H1cEeNhd183s+4PF5KR/JTrlRc1yXAfi5Mx9pObhcNx/mTGh1ny8Q6CXs6BMoNRm7cHJHbCTsnpFH845tG78uw5SLki5huJuxMcZO3HMwbOrNWYKkSwjuLWuvU+jx45zp8P2DXRX76Kz4a+T6lMFF/rk7Xc/8476RXGez/YHqvVq8iulxzjI9DrvPw51Vvd3dvuzcTzxjWDrs9z/nP50CAHzg9k8BAN78W4/hJz+MvBjrMInCD8aYrwD4yiQyFArFuYFa+CkULcVEtL8U80+/wlz2r94BANhzd0U6Vvf64/1f87T+1FG/5k92nTVGXTmVdkqYIaOgo+3jy1mD/YxSPeapfqCYc+SJL8+tcnpm6TdXvEWUg6Pe+BSCU0xOg4e2vVxB5igoAHTPOo2fP99EPuMS/ebXw6dHrmxgByAQyGApcYc8hRgtjN8b8qw+oPC1wm9Vvo+uPSaytBnYNbiy7H70+JLcghmvi8nl07rumaqRgdIzsB9h0wH7ePiyIt92MoKlZHYPRtv985237yNXPvM2bD/qZfzz3/9fAICP/fHrAQD33fZhLB57eNMVfgqFYgtDO79C0VJMpPArhqFaG3vmhkUAwPC0p9xLj3qqv/AIW4u2VG3lEmajyukkt4K1mt0RM42saTJDl9kGOPlASMVqWs5O55pjZ7rK6WxnlVHf/sZTKk6vudVcbfHGpxCMLkrUNgangQ8sSfk0RtD2R81sBW0/N30N1qoFrXu0vU4up8mCqTCfbnBIcjvsmfApoDNnNn1WmE8b+bvSGX9+3CqTX+/QyuPTBm7B56aAnL5zdE+x1Q1L97ltBsdpNtX60BduBgDsEEtuDB35FYqWQju/QtFSzJb2jzz9GVqTzb0/8vzr1DM9nVm+lFF8+4niRiK9Zdm4xNHbHnOE4bTQaVADOso13tyri2ukLTjFrxcD2F00I1lDW2vS+eeWM3k+jRmOO7pwrzxHj/k5UU+8dfWvbxentm5/TK5IryNyqdZ+M1FcK88hTGNCuW6DecyxZy7J5dOcwGlmedy7jq8cGGmqFnHE416G9bSMXQSfItTTL36PBvJcbM2aTHMzbI6OtyCur3346mqlzNwxFM6QoSO/QtFSzHbk7xoM9lRfpt3fq5weFv/RGX/4Qa+24K6P7gvIFVUDtnbPv9zOkSVw8+Tuo4LveMw9uPb95uvT/I5Juhs+eiVcegMfe3ae9yP3+7gyzK0vx5RHUrvCmAURuY6dCIyHyw1MZ7lJ79r4/nDtfmO5/PzQxNjKZErRFEPjzzewP7AKOb5GHyhb1wQKFcQ/8JtcCevcoPm7xG0zXBuC+xFRJM8frV6ylYv9+XPMbHj1EvZCWsV2587K3N0sxm7yOHTkVyhaCu38CkVLMVvaP6Q6ssvS9ZUHX/8HO+vDi1czX+hHmMntjnFvsyh9mh9fEw7W7i2t45R5FKHEtSKK3SVxCsDNdCOUuKb4GSa5Em0Xffd5uZhJrisTCXsmKtaEe8C3A5sBrjjlXnf1fWbnp+Ty5xQoOK1MPo3qjU/PArkRE2YpalDwrsTaK+zj5rduGuFMvgGgy7z2ao/JmBk281Jc3WNNgZnNwSqfAhzzjVyzZu6DXW7uJMsX68wvqlAoLiRo51coWorZa/tteKbuke0AgJX9jJKx4JdBpF7nxJYTlNQxYk4hOVUXtO4lUXZTJrX8eHRdO1WvsCIhnpflu7WBfKy7npSDp1RvbMHByY1R/Y3kQ76PQeRbjgb3IVhpSVH9SF2BDNc2vnbPzYMTbQyu19lIsGkjD4S6tnfcFLwJdORXKFoK7fwKRUuRpP1E9HEArwVwzBjzPLtvP4DPALgSwEMAftsYcyImIxRY/evaUOmBFvM4ozb7GbVZGTc+iZIdgZry6YI38ohoeCVmyfcJ9CzHuy5JTQU6WURHORLa72S7StrKtf2SfUzENDZAYhpT0/7YFKWgvfV9jLSrE5h9C3VF4N7NwFBJ8BANDJJSz4SV7fLIxiyQzfyj1VKIixodnRoJyHlt/weA16zb914AXzfGXA3g6/a3QqHYQkiO/MaYvyOiK9ftvhnAy+z2JwH8LYD3JGsbETrWJNL5U/Por6uXsyi4LLTW2t5qKOyfYvH1uU+5MApwJV8YasqGtWIOJ8HXOvAjX/cfsoNN4MDB7ij3KTeCrEAuK+vaxusSlWwZ6/ziiNNUESkpywoUoEm5Ed//ejvBwIKyMcWcZLsh2BQA/vnFrpGzmtoEOfJeOfNdvp4fPBPhPvVO+p0rl7LIxk/4xq/YbENzD8/bduQrAJvO+S81xhwFAPv/koZyFArFOcKmK/yCdF1nNV2XQnG+oOk6/y+J6KAx5igRHQRwLFbQGHMLgFsAYP7QIeOUb2v7KxpDzLuq/0sWWovR+jri7q5xJWD1Y7zegH4JVCsWvTWl3BFNWLm5LJ9CFCgSjTBdiNkMJBWMXK5rb2S9XVTYJRRz0SlEav2/RG7iOSSnHtI94PsjU5CYwq4uGs1KLBzn70oijFdgi2LlDnaz6SrzIHTTYABYuL/yjl22ivOo56SApiP/FwG8xW6/BcAXGspRKBTnCMnOT0R/BuBbAJ5NREeI6J8B+CCAG4noPgA32t8KhWILIUfb/6bIoVcW19YxdR6y/uO2ah5ii1H9IPSWS47BqH50/dg5N/HAp5InXYzqJ9bIRZocoe+i3Jj3XUIDn2VLkCs3Za6cWkXg9yM2Y2pielxgy5CcIkSmIGLIL76ykKDNqQjEFJmuJMODdfm2DafGVpG6wbvPVgEOVBXX04lYAwSohZ9C0VLMPm6/XYfsWoXb8mHvw0+nfXOC3Oc2TX2QKovF2k9+zSNrwvW+ho474ro3P5xS0sVGie7Gx0WLt5RcXjampCtZ/3e7CpSPAYSRmSIjYd2u2DXE5LpdqeeQugcxy0MeY8FlEo49U1s0pVAEWIgx5hjE33eek8Ip0Oee6I61KQUd+RWKlkI7v0LRUsyU9tMQ6FuTRUf3e8f8gnwsG6qj1UMesTcWW15wsBDbEgu3JcmN0T7BsSf08WaFXX2pcFuRsqL/eSIsVoDYOr+0Xl7gBBRb13br1iXrzin7gRKz4pgyrfaV53YVKfNuvisVpi02xZRsGbhcbmK8MJ6lebjAsjSz/T0b6mvl8qo/iYrFCHTkVyhaCu38CkVLMVPab7rA2p6K0jivvcEeFu30NKM5gilvkCE1FgVVopvCunUYkXdjU+HYurfkn55l+irJ5ZDsB1JTiJjcVFiyErnS+RENvLhiEZMrrfOXeCamYhII2v7A47JkahKj9S7UXCrBS+xd4tMQl3iGRQcOIgGz6fHqvqoRHZflNyfUXdhkhULRNmjnVyhaitka+cBTbGe0wPPojRjN4fulZBExJI046oKJ46nzN1Nuk3pTx5tOEUrqbxISLKOsuAoxodwic+kcubmeiRnPwcmIZZIOxDman/LIFKAjv0LRUsx25B95E92hXcvkJrujBa7gYMo/l4IrwzTSIViPlXzlc9bbUyN7KjZ9wy9/NqPIGQlzcwCUyC0w2W18byW5sfObyM0xd05dAz9NGPlFm4+cd8UpJbmJO0sD1llifcbaGtRK6/xlfh35FYq2Qju/QtFSzJb2E1NcWN7NI+vOs+yjKxf7BdPemYoz8XXPaCosgX5J04WmYbEC882UTQFHiXIwRdU3krlerjQ1icnInSLE6GoiVFnTqYnkxVjkMVli7lzyXrDzRoJZuRTJN8umwLaBm7MHVJ/1g7kT1rz3IjvH6OTzfh35FYqWQju/QtFS5KTrOgTgUwAuQ0VIbjHGfLRpyq46qIGlzJ1Fz4NWD3j1aGeVafu3ja8MJL2vEokYRpGkDaGQcflRr726ArbNqWUijJdINyPTjZKkHUVpvlJadbekHKPGTVYsYvXy01Lr9JLcDA29dLzIE5NBjLhc4okpXE9n2e/kHoT9UyyM18Gq4o6dGmcZwwhVxjAA8C5jzHMAXA/g7UR0DTRll0KxpZETwPMoAJed5zQR3QPgcjRM2dWpLfyqT+AokoySO0i4gIUxKyex3UKSRMB/GJOjfUxu4nMppYaKyi0IHzZxos6csplyo4NLydp9qkxTWamyqfvRIJQZwO5JIsBnTr21gjPyvq/tZIpvy5xd2LuoIlxA0Zzf5ux7IYDvQFN2KRRbGtmdn4h2AvhLAO80xpwqOK9O1zU6o+m6FIrzBVnr/ETUR9XxP22M+ZzdnZWya326rpGjMvaz013kWq9xs8Xqh42MykJ7RX37Bd9wMdRUTJmW8hPnh4X1ZxNT4hQoj6aaj77k/CaOPSkq3tQBJ9XeErk5UX9T7cmVG5Evmf/G7h1XUNdVRUKNObrvUnjxPpJCTsYeAvAnAO4xxnyYHdKUXQrFFkbOyP8SAL8L4EdE9AO77/2oUnTdZtN3/QLAGzaniQqFYjOQo+3/JuIEqCxlVwcwls73n6pIxxrLRBpLOFBH5GUZfWOZdUUqLq0oxNZjS5JrJNbjN0tuURivujGR4yVym3jX5Wjdc+WWJO2I2WaUxAZI2WYIcmM2EPV2xjq/aAocsQlZs2G8ek/ZwhrGS6FQpKCdX6FoKWaetKN3pvreOFPe3knPbZzhDxBOATp2t5gkAQgoUTa9LjFRTXmxpShmpGyR1l6iiznmvevrn0RuSYCOgkAYybIlnonC+cnn0NAoSgrcEdPmSwleAiRWpEZsmstXulzQm5HLcF0wnOvIr1C0FLON298BBjZUV2fJmuzydUm2SC5lzs0yXZw0fFTTdWknquT8HMXZRsdz6iopK9VVUrZE4SfJSDGDWL0pNA2BJiFyvUaSK7/aaQgKzsCuhTNdpwx3Tm8axkuhUKSgnV+haClmHMbLwMzZDKSnbGguFqqou8T8l1PZRguUZaIPflMPsRJ6nrN/M9B0mtNkilByHyel102VdMJxilHyph6RTSI5Bw1im07Bzc3ZmY3LiCnG6/1U8iDsucVnKBSKCwLa+RWKlmLGSTsIXRtuyK3Z8+Qcgz2en7tyADCaLwjj5YJ1sGAg3EuqngKk1u55mYhppZiIIYbUOn8TGpuzzt+k3obhtqZK1VPtSnliNvWYbGLuDHmdPwgl1xvfF7s3rmx3kYWyY1S/y8J7uejXdXq7gqmVjvwKRUuhnV+haClmnqW3Tipg6clwF6MzjOpzrz1Hc4IAHwmDECPk5wPS3lcbyQSQju5aEjtuUtPYHK++Jpr7Ei+3GEqMaXLlpsySJZnrdye8+qbpiRkkjhGi+8bgysbM3Xmim+7pdfEt1chHoVCkMNuRv+MdEDrL1ZeMx+cPTH3Zl92l6QqWMlNKqYQSrig/e2zESbVFKlPi2DNlE9Si83KPp5hOCUvIkbvROSXHI3WVpOiKyaiLNkz95dkpM3dnfvq05MvWbNjJV4WfQqFIQTu/QtFSzJb2GwADR/er/4PtjC+xJAWdxfHpAHHPplgCD2F3KmJqUcZXjpSfeYq6ligHU+c18QosPS91voRpKEBTKDFnLpGfUERKyWAC+i7VlxMdeiTI5+87r8Mqxrs2TkZJdOKc6L3biOjviej/EtHdRPTv7P79RPQ1IrrP/t+XX61CoTjXyKH9KwBeYYx5AYBrAbyGiK6H5upTKLY0cqL3GgBn7M++/TNokqtvBHStln9go/ZyKk+c6vMpgI1mmhPGS1xfFqiWmEE1UjZFz6ZKqXnZmFlxk5WKnHblUvgZXG/jSMDS+amovw3llmjwU2G8+PMdJXolXxXrn7Bh8fbbCqaZtAMAiKhrY/YfA/A1Y0x2rr4gXddZTdelUJwvyFL4GWOGAK4lor0APk9Ez8utIEjXdfiQcf77TkEx3MZGc5KVGs4xR7LUG69Q2CeMKMnRPiaz6SgzoSNLMgglR1MHm1y502ARTeSWWDSmrAGnYZshyEjG7Y/UKzoERerqrPgDLu9FHQtjlH/ji5b6jDFPoaL3r4HN1QcAG+XqUygU5ydytP0X2xEfRLQA4FUAfgrN1adQbGnk0P6DAD5JRF1UH4vbjDG3E9G30CBXn3NQcMo70T8eCGhOneE0h6fkmn028Z/P2T8Nk9wm9c5S7jTaNU25TY43NTtOHJ+KeXCiveEUgcJ6C55Bjrb/hwBeKOx/AqW5+hQKxXkDNe9VKFqKGYfx8ma9Tsvv1vABYDTP1valkF0F67HJ6URTjXXT4yUa69x2lWispyk35zlMahMgyS15Zk1jJTQsmwwPVnINzquPm7NHslL71asS982gGoVC0TZo51coWoqZB/Oo6b6j/7u8nW3/hOfn3NuvDoHEaVKMXjkaJIXbAppphktMRacxhZCmJtJ5TTXHk05NSuqdxqpKCVWfsF1Nw7s5k9zAEE3K4puxouHKuqjVQJirj5u511l6hTaloCO/QtFSzNyf330Z3der/5T/PA4XmKmv4OgQOPPERnPJcacgj3qRqWgKqVF1moq3HLkppOpNOdhwlIzWuUrFEoVfAYvIsTUR6+W7hZG9aZq4OvMuT9HFRvsez3WxY90LWbDOryO/QtFSaOdXKFqKmdN+Z95LNn4/9+ozHUb7WbRSZwsQTdG1ro6xsvywO6+EvjdVLpXIlVCydp863lRWU7kl8pvchyb3ICIrujafmrok5JsS5TDf7aYLvD+wenmfce2p7WJU4adQKFLQzq9QtBSzT9phtZbOHJGnIaoDEiA09XX0KJYBNaU5FlMwbdb6c855qbJNpwMblc253iZefU1WE0rlTno/mkwbpiw36ZnKmbwtw8N18XV+3g96NvTdwK2UqbZfoVCkoJ1foWgpZp6l19H97uK67KJYF8+PUfVRPzQJBtYZYwQV2PM5jUoZ+TTM+CpS6pIYfpN6psVkTdO7LmXkMw1z5iaeeCm5TU2BmxpmJdpVv2+xd0J4b2JUP4jht6d6oesM16rtVygUKcx05KeRN00c7Ko+by5bLxA3yXVmjoFPc8K810S+0LXCJXb+pKaisfzum+FHnnL8iZUtGe0luU2vIYYmjjvTdF5KmXc3lSu8F1GFn2AHEHMC4srw2rHHOQFthsLPxu7/PhHdbn9rui6FYgujhPa/A8A97Lem61IotjByM/ZcAeAfA7iV7b4ZVZou2P+vT8kxHWA4bzCcN1W23gFVdMn9dfwfjfwfOqb642WjlQjHjfBH7C+2fyP5Ri5nyP8F2EjmeqTalboHJWWbtLHkPm5WvSWQ2hWT1fTeJeoV34lYvQ6sP2Dk/2jo/0b96q+zSpVCfBMUfh8B8G6EM5PidF3DM5quS6E4X5CTtOO1AI4ZY77bpAJjzC3GmOuMMdd1d+5oIkKhUGwCcrT9LwHwOiK6CcA2ALuJ6E9h03UZY45mp+siU2slu2crdT4PUsBDIHGIGtKC9dogWEOJaWx9UqLepuveOfVthGmYDzfR0De9ltR5TTX0EqZpQ9FUroSM++FsVEYd3jfkNX9XtjaTnybtN8a8zxhzhTHmSgBvBPANY8yboem6FIotjUnW+T+I0nRdI0LPWiINt1XDcRCSaBcLWLjCh277n33VooEWbdnoeruEknX+ptZkTazFmo5eEkrYScoSr0RuCROaxmidKCuut0/hPtZx+yPOZ8kwXtyExQUDZWG8uCVsINeWHe60DYhZvgoo6vzGmL9FlaVX03UpFFscat6rULQUM3fsGVn64px0nGMCEPrzB/TJOjhw/+aUX3TScSdGv5qYksYwqYkqIsdTirlZyN3ofL4/py25JsZNlamsbDKeQxNzZy43Fgk4FcaLbddh6/oy1R9xBzjbj5yjXEl4Oh35FYqWQju/QtFSzDyMl3FhvGx0XrfeXx2PaDT7AhdLePUFKwOcnknnpFCyDhw7r+R4ino2kZvT7s1Y3885p0l7N8sWISWjxL6Ev8PSylCkDW66YGLTAjY9dn1mpFl6FQpFLrTzKxQtxTlL2sFDcjkMFjyXr7WXQJ3MIwhrlKBPsTBeIsVL0b7YFCIVmqmh5rjI0CRX7matUuS0YVpyMwxksuXmehFuVFZ4L6Lh5VJBZIQ6nNYfWLe6xabHQ9tneH/JhY78CkVLcQ5G/mrTBeU0PPvoGf/ZDBx+UmG8OJx/A/usSWaWFMmjnoK4djuF9edGjiwlLGKa69o5CtDNYCcl97mh407UbDxxvrjOL7yDMfsUqV7Deid/X4csM+/ck9ULubZrE8N4KRSKCwva+RWKlmL26bqsaWKwVmkxinguuRgAgcIv9dmK0TfnfdWTyybXdlOmwrPwYkthUu+6lFdf0+udVG6OR2Sup2WKvpfWK4AkU9/Iuya+z/x89r72Tvu559rO0Fx+M8J4KRSKCwza+RWKlmL22n7noSeZI3IaFJhJ2nNKPlUxmjxNSp1b17TllmBCE9Wi41tJblOz4wK5YuKYKZiVj7i33/ogIqrtVygUKWjnVyhaiizaT0QPATgNYAhgYIy5joj2A/gMgCsBPATgt40xJzYUZLzG3hnxENPgc8MerimVKFNTY4ykmeWkpqIlGutJ5TY0ZJlYbtPrnbRsySpDpGxNkxPBYMbqSBwX5TY1+LIyuGEPP3/Eeq0zAa4DfGwS7X+5MeZaY8x19rem61IotjAmUfjdDOBldvuTqAJ7vmfDMzrAcLuN2num+u6s7fOft95J35zAP9l+zaSopdUPYVtY2x/bls4vWZ9OHC9iJ01MVJsqrSaVO4UROClXen5TYDfJMF5SeyLr/Pz5utGYv6ONbB0AZgLP9gVyed+oTuza/gReLoHckd8A+CoRfZeI3mb3ZaXrUigU5ydyR/6XGGMeJaJLAHyNiH6aW4H9WLwNALr79zZookKh2AxkdX5jzKP2/zEi+jyAFyEzXZcx5hYAtwDA/KFDprNc0ZSBTTLQO+GbIIbrQtr0UaJP0aQdKZNNSW4J3eRFp2gqWkSpcxWV05Y7zalJk+lZjvwmyuNIG2IJOjaUm6HwdYpvHvOCK8ODWBd2ejxwUbC7+S9oTqLOHUS0y20DeDWAH0PTdSkUWxo5I/+lAD5PlWKhB+B/GmPuIKI7UZquS6FQnDdIdn5jzAMAXiDsL0/XZeDDeK1Uu0bznqaEOfdYMI+1ca6UCuMVi3wqIpXoIMakSuj1pHJL6ijRhE8qt8Qkdxp2CalzGjyT5Jp/7BoT740YPi5jSiW9787EHQiTdrhENnVev4I01Grhp1C0FOfMn7/+QLH1/O5p/y0azbPzRtaxhysEU6NFLIzXpA40qfXnpnKbKBVLlGmSLURTuU2Zw6Ryp8EyEseLlH8ckuUor6LEstAi/EYDAAATDklEQVSJYqw4UP6xftA9uy7EnfrzKxSKFLTzKxQtxcyz9LoY/LXij4Xz4lFJO8tsCmBjkzsbAQDxiLtunZ9H5+VRUF0U1VQGVSYroIWS/UDKLFWSuV6uRDcnDYvFyzR1XpLklqy359QryZVQYlacmkbF7n0qOnNK+SeY6QLeVDdq/ssdd2wbgj4w7xvZO+v3D3ZX+2mlZI5p5RefoVAoLgho51coWoqZ0n4aecoy2F5xMeflB8hUH/DmjNzEMeWPH6P1jlInzX8h7xPLpqj++jbmyi2h6hG5RVQ9N+ZATrtS69pN5TY5noAU3bm4XuE+c6+8zkCoKyLXT039TjdNBrxpPOC1/XU/KhjOdeRXKFoK7fwKRUsxU9pviBko0LjnEo/oy2lOTX+m6NWVpPqJ8wM0NdONyWiiVZ9GvU3OKdCqN76eJkjVm/PMSjwTBXnBKkKTVZGY+S/T7A+tIRC5fpQyVWfQkV+haClmbN5rasVE167vm0gLeGDPmhE0XV9mX0MSRv7GJp2pEbppe90urhCcdEQqUdLFkLseP42yJe1qcp+nGRsgJiPFEiP1uvcxCOMVBLRlJ3bDsHhUwOR05FcoWgrt/ApFSzFb2j+kOsPo0K5V9p/y3x+XcbQ6ML6mz9fuO9x8N4gDYP9FPPmkbKlRU9/c9fQYFUx5ACZoXzIGfA51bTI1SZXN8YKbcGrSyNxZkgmIHp5Z8fVTVD4hV5y2ReoNIgELYbwQvMPMq8/awKweGI4dS0FHfoWipdDOr1C0FLnpuvYCuBXA81CRld8DcC9K03V1PKXpnbJ0Zb/n59y8l9OnjqVEUe8r3lZhiiAGskh56nE01NoXyeW7UzQ31ww30q7o/lyqnjJL5jJKbCCmaUY9C0/MhFyeVqt+dyPTymBKu+JWwpjdC3/32ZR46PrTU1bAMGeeZevJLPdRAHcYY34FVTy/e6DpuhSKLY3kyE9EuwG8FMA/BQBjzCqAVSIqT9cF/zVb22PX+8/K3x/pqxrESY98+aUvsFSWjyzJtf2SETZSR5HcVF0NbAaK1tubHk/Vl8M+Njo+DacpiZEUWMVF5Qpt4e+wZF8SIGC1JvhfnecFc4ZcC07lBRCQM/I/A8BxAJ8gou8T0a02fr+m61IotjByOn8PwK8C+Jgx5oUAzqKA4hPR24joLiK6a3jmTMNmKhSKaSNH4XcEwBFjzHfs78+i6vzl6boOHzLOYceZ7/YWPU9ZucTzpDrrKPxUgTv7UIJ+JU12I+a/0TK5cmOYppJOOqepXElGiU97TG6qLQm5qecbVdIJCj+Oej0+1tYS5Z8kVwgZF+zPuHfueutY/ADQkZV/LoxXXXaatN8Y8xiAh4no2XbXKwH8BJquS6HY0si18PsDAJ8mojkADwB4K6oPh6brUii2KHKz9P4AwHXCobJ0XUBNS7rWJ3nloLfT5Rl7hyxhQddqN3myglS6Lr7G2lnz285Tiu8TtcGCTCBiCpzjgVZgKirKLaCgE3s/puQWeNxlTZMEGUlT4YbhtqTVoCyzcUFuMDVx7xs3020Yxsu9m6NtcsTetYv8y9s/3rd1CZ6vCaiFn0LRUmjnVyhaitl69ZHX3A9tLr7uUxGqv8RCFe10Gk0W4CP22bJFOpGkHTX94ldeEsZLoJtF1DbmNSbR2EkjycbqTZ3X9Hgufc+RK6HEM7Fg6hKl+olnFjRXMEEOppbCtCDWxnpqyt73AUto03vczydcFOzAAzATOvIrFC3FbEd+gyA8FxB3XuAZSt05JevqsRFHdOwogfCBzWpXE9PZaZjbbobcSc10Zy03UTbFIovqjSmHU3YeQiyDwJyd2bgEfcOVdf1oyua9CoXiAoR2foWipZg57e9YCu/WJcX4/AinA06ZIa2bAmXrsW79P6nkAfIVUTnpngoUb2JMAknuNENzlZQV2gqkn0Nj+4FUuzaSua5sfW9TCtJYHZPKjb1rks0HX/vnsfp3sXRdp236u73WlrggfK+O/ApFS6GdX6FoKWZL+7sGgz0V355/rOLwK5ev1ofprG9Ob9F/l9waZ5ftG26TVwkkChmEBBPMLEXvKy6LQ/K+ktZz18tNfWYlLXFkCiFG903IjZqwxqYWEgT6XRT5OLHOH6z2SPexxJw5NjVJ3NtUyLdALtfA23Ba4rvI6g3eFX7vuFy3usWmviP2vmMwPgXYdrQXnJsDHfkVipZCO79C0VLMPmnHyarKlYMV/+k96VX4Ae1jHnxkI5IOF3wBTrmkKKgxrz0jUcw5VpckN8KkRE+tVDKQmDaY75YSQEhyCzwIY0lMipKUCOViU5tJpyZc7nqZY3ILjG3qexupK6qtl+T2I1PPDfZFo06zNoxsPssOM3Hnhj1c899/shK4fLiaPps51fYrFIoEZjryk/FfuP6TVdV8hB/N+bKjef/Z7NiMvp01FtorkpaoVmbxr7bwiTPcDmCQGPIiH1MxBRNHSoFWELZKtEUocLAJkONcJB0X5MfiKjRa80/cr2QuBi43dV2x+526jxG5JDwTMQN17F3iCk5r10IsBn9nRW7v8qX2hXd9Q/35FQpFCtr5FYqWIidpx7NRpeVyeAaAfwPgUyhM12U6fs3eRRsdHfCaud7RubCwRW0SzJQZUWWYZUqBt6BA68UMqkDc/NLtkhRFEbqajEAbkSvSekFuY2/CpnIbxCRoHO04EUW3xKxYMu8OFIoxmw9pChGTK5mr85gBa0JbI3KdQo+b8RLz1w/Ce9ms125akJzCMuRE773XGHOtMeZaAL8GYBHA56HpuhSKLY1S2v9KAD8zxvwcwM2o0nTB/n/9NBumUCg2F6Xa/jcC+DO7HaTrIqKsdF31mv3uSku5cP98fWzpcm932j/hefCazeTbYea9Ma88R5+7S7IpcE2/YlcurXvHQjPZumIUVFz/j0wxSkyBkyaqJZr/ArNi0btuUs/E9dsblCX2zLmHZ0DhBapuhOcXaNeZSSxfu0+tfoy4ybTgedpdZma49h2U7BcArEvEURXqn2J9YK8/0ZnyAsDK/urEgQ11x6e7KWSP/DZm/+sA/EW29Oo8lq7rbMmpCoViE1Ey8v8mgO8ZY35pfxen69p2xSHjAnN2l6qql67yjj1zj/rP5upBpgi0AQtjCj8OxwiCUEfcekqK2x+xutrofC4jGIX42m9qNI8prSTf75GwnePo4mRkjLqSZWFUbgImVa8kN9YuK6NpPAdxTT+S7z75zLhiTmAiXe53z4JuOiVe7N7yd9A5sPGAtr2nfIHly1nc/serfrR20eaG8XoTPOUHNF2XQrGlkdX5iWg7gBsBfI7t/iCAG4noPnvsg9NvnkKh2CzkputaBHBg3b4nMEG6rrU9Fb/icftH86zYsqc5NSXKoDT12iqnm5LCJ0H1A5kRm4JaERWjmAUULJktF8LxkhBZGfUm21sSEix1vCAMl+SMld3WjHqjVD8RcVeKCRGYnQ+FE3PugXC9LncFAPSe9H3GKfrIprTDSP35FQpFAtr5FYqWYqZefaYLDHaFXJabIw72efXp3HHftLWdFVeSsu0Cst90uM7L1vkH1qx4jmlimelkUvMvhJrKigQsCpPLOuoYtIs/qdxwW8hYb5eQosl8Kbwkw20MudOY2P1MtVfInCu+M1hnE+DiNRTYUHAv1c7a+HsVjR3Aqxg6mwG2usXe196Sv+mjK6rVsu0/XLB1yjIl6MivULQU2vkVipZitsE8hkD/ZPW9WTtUuSH1HvNN6Kz57dVLmKnvE86Qwe/jGk9u/CMlPAgiAbuMv8z8l3tJBVQ7EVTDGXQM+fkrTC4LSFJPb6JUn8ldHjfyCAxKnGdajL4LXmz8uvj9EiMPJ+g1j6wbTk3Gw1olw20xuUGiFnYf3f4ST8xYmLaOEN028L5L3OeoXBeAY01+f2pPPWb4w8vy7cHOoa3L75v7pb/pdQAPADt/UNF98w9PVjs+G7EfFqAjv0LRUsxc4be2p/ry7b5zGwDg5D/gWjy/SSvjo/U8c2hY3SMHT6wdLNgoMmBrpM4vmo/2PB/ASAiAGBth3YjPwy3x0d6N4ECovJHAR+BaLo8L3x8fZXJG1doEmQcp5ebKwqgYVXBRWD8Qd4SR/NsDSHL5aC+Y3EpBXflxwN8TiUkF7YoEWg0so4UArZz1cOcx976Fcr204Q6rxGXOPrzeEQtO2z/ZHat3uI21i11P59erEBprd+2rzlnMN2DRkV+haCm08ysULcVs4/aPWCTeG5+o9h3ZUx92dAeQ/Z5XLvI7Qzo5rswKlE9MmeKUdCPOVjnVFxRyAS3ka+i2DVwRxSMMB1Q/sZbN7RZcJQG1lTwAE+ahXK5klgyEFF6KXyB5x/F7y20ouIJKnJrE4ORGpkY+FwO/MFnR6I+zTd5eu214W/naeMI6NrA1mRu/D8H7wZ5fd3F87T6Qy6cDdnOw1yv2aMG/+71HvR388o/3AgD6tmhBkl4d+RWKtkI7v0LRUsw8acd688Nnf3yx3n7o3Z76rB7d4Qs52sZMHPsszddgB6NfVkTvDKPfjLYPF6qy3LyXTyGk9WVObXtnx1cGAk08n07Mj68oxMyH+fqyZN4ZpiQTzH8j4cOcXH6NPO0Zn6a4NWhaEaYCfJvtCzTpvfH28vRRPAJtAKehZ4d7i3zdO4z4DADD7ez5ce85wROPTwuchp5r4vn5Q2FliCPU0Hsq7rJOr+73+7qL/J7bd4Uno2HXw9s7uMSmsjvGUtkxu5Ydj/jCJ6+r7GXm99p+dJuu8ysUigS08ysULcVMaf/8jlU864afAwDozdW+3/ubb9bH33/b7/jCh31ysoV7KguHtd3+W7W2W6Z9/VPWjPIZy/W+3s+8hcTKgYoHu7iAANBl2VBX+YqCpX2BeTD3SnSbjK9yWthj0VdHLoKwkJwBCBM0dFzTuXKbTVPmbGbW2D0INPh2HhKYlbLjnD53T3fCtq5v79r4dISbNvMpgPPQnDvGAk/skOW6+9hl7Rqwa3PZagd7mMabrQzxeI0OgeEXm6oNbBTc+WP+/NUr/bvWe8Rr0uuVIz6lYlOIbSzm5PLB6nr7LNbe4Gk+PmX3mE1IY2Sqz+/53u9WZfs3Hfd13bKv3v4nH/hqvf3Rb70KAHDgv+4EAPz8eP54riO/QtFSzHTkX17t46e/uAwA8N7/fQcA4L+/6bfq42/5xDfq7W++9LJ6e/uXqi/s9+96Vr1PUswBwODKatjc9uOFet/iM/0XeOGB6qu6csAPDc6RAghNJ9eHHAOAuSfYiGEjDC886NOMBXL5SGXTKvFRUzoOeIcerhzts7Bma/usAoyvDXMT1AVB0dhhI+nQF+6dZKznQNWegLFIo7yRWQZnEXPWYcu1day93HbCPr/AFIGt6TvWsvCwf12Xns6iO58Yd/LiysMhi3nv8tmvsAi4C/d5Zrh8KXsXnPkut2Vg7eIONu69mLvmZL1v1195G5Ynrx9X4nHGsu0xf0O2vbYKkL3yBZ8K4+X//u/q7c+/59X1du+l1XmPP9+yo+9pGC+FQpGAdn6FoqUgYwrsASetjOg4gLMAHp9ZpbPFRbgwr02va+vg6caYi3MKzrTzAwAR3WWMuW6mlc4IF+q16XVdmFDar1C0FNr5FYqW4lx0/lvOQZ2zwoV6bXpdFyBmPudXKBTnB5T2KxQtxUw7PxG9hojuJaL7iei9s6x7miCiQ0T0N0R0DxHdTUTvsPv3E9HXiOg++39fStb5CCLqEtH3ieh2+/tCua69RPRZIvqpfXY3XCjX1gQz6/xE1AXwXwD8JoBrALyJiK6ZVf1TxgDAu4wxzwFwPYC322t5L4CvG2OuBvB1+3sr4h0A7mG/L5Tr+iiAO4wxvwLgBaiu8UK5tnIYY2byB+AGAH/Nfr8PwPtmVf8mX9sXANwI4F4AB+2+gwDuPddta3AtV6DqBK8AcLvddyFc124AD8Lqudj+LX9tTf9mSfsvB/Aw+33E7tvSIKIrAbwQwHcAXGqMOQoA9v8l8TPPW3wEwLsR+tlcCNf1DADHAXzCTmluJaIduDCurRFm2flL8q5uCRDRTgB/CeCdxphT57o9k4KIXgvgmDHmu+e6LZuAHoBfBfAxY8wLUZmZt4fiC5hl5z8C4BD7fQWAR2dY/1RBRH1UHf/TxpjP2d2/JKKD9vhBAMfOVfsa4iUAXkdEDwH4cwCvIKI/xda/LqB6/44YY75jf38W1cfgQri2Rphl578TwNVEdBURzQF4I4AvzrD+qYGICMCfALjHGPNhduiLAN5it9+CShewZWCMeZ8x5gpjzJWons83jDFvxha/LgAwxjwG4GEierbd9UoAP8EFcG1NMWuvvptQzSm7AD5ujPmjmVU+RRDRrwP4PwB+BD83fj+qef9tAA4D+AWANxhjnjwnjZwQRPQyAP/SGPNaIjqAC+C6iOhaALcCmAPwAIC3ohoAt/y1NYFa+CkULYVa+CkULYV2foWipdDOr1C0FNr5FYqWQju/QtFSaOdXKFoK7fwKRUuhnV+haCn+Py3vLU8yLV/5AAAAAElFTkSuQmCC\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"H1 = np.fft.fft2(f1)\n",
"H2 = np.fft.fft2(f2)\n",
"H3 = np.fft.fft2(f3)\n",
"\n",
"plt.figure()\n",
"#plt.imshow(H1.real)\n",
"plt.imshow(H1.imag)\n",
"plt.figure()\n",
"#plt.imshow(H2.real)\n",
"plt.imshow(H2.imag)\n",
"plt.figure()\n",
"#plt.imshow(H3.real)\n",
"plt.imshow(H3.imag)"
]
},
{
"cell_type": "code",
"execution_count": 14,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<BarContainer object of 3 artists>"
]
},
"execution_count": 14,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAY0AAAD8CAYAAACLrvgBAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvnQurowAAFGtJREFUeJzt3H+MXlV+3/H3p3hDaVOIgdkVsnFNFicKoMaRLS/SaiMap+BsqsBW0AyqFldF8i5ipayUPwrtH2xBSEurDRJqIWJly4C2/ChkA2ogxFq2QZVYYNhQfi5lWEiYxQIHOyzRBlp7v/3jObN9Zng8czwzeMbwfklXz32+95zrc3X/+Pjcc59JVSFJUo+/t9wDkCQdOwwNSVI3Q0OS1M3QkCR1MzQkSd0MDUlSN0NDktTN0JAkdTM0JEndVi33AJbaqaeeWuvXr1/uYUjSMeWpp57666oam6/dRy401q9fz8TExHIPQ5KOKUn+sqedj6ckSd0MDUlSN0NDktTN0JAkdTM0JEndDA1JUjdDQ5LUzdCQJHUzNCRJ3T5yvwjXx8v6q/5kuYfwkfXa1397uYegFciZhiSpm6EhSepmaEiSuhkakqRuhoYkqZuhIUnqZmhIkroZGpKkboaGJKnbvKGRZFeSt5I8N1S7O8nTbXstydOtvj7J3w0d+8OhPpuSPJtkMslNSdLqx7fzTSZ5PMn6oT7bk7zctu1LeeGSpCPX82dEdgP/Gbh9ulBVvzu9n+QbwDtD7V+pqo0jznMLsAP4HvAgsA14CLgcOFBVZyYZB24AfjfJycA1wGaggKeSPFBVB/ovT5K0lOadaVTVo8D+UcfabOFfAnfOdY4kpwEnVtVjVVUMAuiidvhC4La2fy+wtZ33AmBPVe1vQbGHQdBIkpbJYtc0Pge8WVUvD9XOSPIXSf48yedabQ0wNdRmqtWmj70OUFUHGcxaThmuj+gjSVoGi/0rt5cyc5axF1hXVW8n2QT8cZKzgYzoW+3zcMfm6jNDkh0MHn2xbt26zqFLko7UgmcaSVYB/wK4e7pWVe9X1dtt/yngFeCXGMwS1g51Xwu80fangNOHznkSg8dhP6uP6DNDVd1aVZuravPY2NhCL0mSNI/FPJ76TeAHVfWzx05JxpIc1/Z/EdgA/LCq9gLvJjm3rVdcBtzfuj0ATL8ZdTHwSFv3eBg4P8nqJKuB81tNkrRM5n08leRO4Dzg1CRTwDVVtRMY54ML4L8OXJvkIHAI+HJVTS+iX8HgTawTGLw19VCr7wTuSDLJYIYxDlBV+5NcBzzZ2l07dC5J0jKYNzSq6tLD1P/1iNp9wH2HaT8BnDOi/h5wyWH67AJ2zTdGSdLR4S/CJUndDA1JUjdDQ5LUzdCQJHUzNCRJ3QwNSVI3Q0OS1M3QkCR1MzQkSd0MDUlSN0NDktTN0JAkdTM0JEndDA1JUjdDQ5LUzdCQJHUzNCRJ3QwNSVI3Q0OS1G3e0EiyK8lbSZ4bqn0tyY+SPN22zw8duzrJZJKXklwwVN+U5Nl27KYkafXjk9zd6o8nWT/UZ3uSl9u2fakuWpK0MD0zjd3AthH1G6tqY9seBEhyFjAOnN363JzkuNb+FmAHsKFt0+e8HDhQVWcCNwI3tHOdDFwDfAbYAlyTZPURX6EkacnMGxpV9Siwv/N8FwJ3VdX7VfUqMAlsSXIacGJVPVZVBdwOXDTU57a2fy+wtc1CLgD2VNX+qjoA7GF0eEmSjpLFrGl8Jckz7fHV9AxgDfD6UJupVlvT9mfXZ/SpqoPAO8Apc5zrA5LsSDKRZGLfvn2LuCRJ0lwWGhq3AJ8GNgJ7gW+0eka0rTnqC+0zs1h1a1VtrqrNY2Njc41bkrQICwqNqnqzqg5V1U+BbzJYc4DBbOD0oaZrgTdafe2I+ow+SVYBJzF4HHa4c0mSlsmCQqOtUUz7AjD9ZtUDwHh7I+oMBgveT1TVXuDdJOe29YrLgPuH+ky/GXUx8Ehb93gYOD/J6vb46/xWkyQtk1XzNUhyJ3AecGqSKQZvNJ2XZCODx0WvAV8CqKrnk9wDvAAcBK6sqkPtVFcweBPrBOChtgHsBO5IMslghjHezrU/yXXAk63dtVXVuyAvSfoQzBsaVXXpiPLOOdpfD1w/oj4BnDOi/h5wyWHOtQvYNd8YJUlHh78IlyR1MzQkSd0MDUlSN0NDktTN0JAkdTM0JEndDA1JUjdDQ5LUzdCQJHUzNCRJ3QwNSVI3Q0OS1M3QkCR1MzQkSd0MDUlSN0NDktTN0JAkdTM0JEnd5g2NJLuSvJXkuaHaf0rygyTPJPl2kl9o9fVJ/i7J0237w6E+m5I8m2QyyU1J0urHJ7m71R9Psn6oz/YkL7dt+1JeuCTpyPXMNHYD22bV9gDnVNU/Af43cPXQsVeqamPbvjxUvwXYAWxo2/Q5LwcOVNWZwI3ADQBJTgauAT4DbAGuSbL6CK5NkrTE5g2NqnoU2D+r9mdVdbB9/R6wdq5zJDkNOLGqHquqAm4HLmqHLwRua/v3AlvbLOQCYE9V7a+qAwyCanZ4SZKOoqVY0/g3wEND389I8hdJ/jzJ51ptDTA11Gaq1aaPvQ7Qgugd4JTh+og+kqRlsGoxnZP8e+Ag8K1W2gusq6q3k2wC/jjJ2UBGdK/p0xzm2Fx9Zo9jB4NHX6xbt67/AiRJR2TBM422MP3PgX/VHjlRVe9X1dtt/yngFeCXGMwShh9hrQXeaPtTwOntnKuAkxg8DvtZfUSfGarq1qraXFWbx8bGFnpJkqR5LCg0kmwD/i3wO1X1k6H6WJLj2v4vMljw/mFV7QXeTXJuW6+4DLi/dXsAmH4z6mLgkRZCDwPnJ1ndFsDPbzVJ0jKZ9/FUkjuB84BTk0wxeKPpauB4YE97c/Z77U2pXweuTXIQOAR8uaqmF9GvYPAm1gkM1kCm10F2AnckmWQwwxgHqKr9Sa4Dnmztrh06lyRpGcwbGlV16YjyzsO0vQ+47zDHJoBzRtTfAy45TJ9dwK75xihJOjr8RbgkqZuhIUnqZmhIkroZGpKkboaGJKmboSFJ6mZoSJK6GRqSpG6GhiSpm6EhSepmaEiSuhkakqRuhoYkqZuhIUnqZmhIkroZGpKkboaGJKmboSFJ6mZoSJK6zRsaSXYleSvJc0O1k5PsSfJy+1w9dOzqJJNJXkpywVB9U5Jn27GbkqTVj09yd6s/nmT9UJ/t7d94Ocn2pbpoSdLC9Mw0dgPbZtWuAr5TVRuA77TvJDkLGAfObn1uTnJc63MLsAPY0Lbpc14OHKiqM4EbgRvauU4GrgE+A2wBrhkOJ0nS0TdvaFTVo8D+WeULgdva/m3ARUP1u6rq/ap6FZgEtiQ5DTixqh6rqgJun9Vn+lz3AlvbLOQCYE9V7a+qA8AePhhekqSjaKFrGp+qqr0A7fOTrb4GeH2o3VSrrWn7s+sz+lTVQeAd4JQ5zvUBSXYkmUgysW/fvgVekiRpPku9EJ4RtZqjvtA+M4tVt1bV5qraPDY21jVQSdKRW2hovNkeOdE+32r1KeD0oXZrgTdafe2I+ow+SVYBJzF4HHa4c0mSlslCQ+MBYPptpu3A/UP18fZG1BkMFryfaI+w3k1ybluvuGxWn+lzXQw80tY9HgbOT7K6LYCf32qSpGWyar4GSe4EzgNOTTLF4I2mrwP3JLkc+CvgEoCqej7JPcALwEHgyqo61E51BYM3sU4AHmobwE7gjiSTDGYY4+1c+5NcBzzZ2l1bVbMX5CVJR9G8oVFVlx7m0NbDtL8euH5EfQI4Z0T9PVrojDi2C9g13xglSUeHvwiXJHUzNCRJ3QwNSVK3edc0JGkprb/qT5Z7CB9Zr339tz/0f8OZhiSpm6EhSepmaEiSuhkakqRuhoYkqZuhIUnqZmhIkroZGpKkboaGJKmboSFJ6mZoSJK6GRqSpG6GhiSpm6EhSeq24NBI8stJnh7afpzkq0m+luRHQ/XPD/W5OslkkpeSXDBU35Tk2XbspiRp9eOT3N3qjydZv5iLlSQtzoJDo6peqqqNVbUR2AT8BPh2O3zj9LGqehAgyVnAOHA2sA24Oclxrf0twA5gQ9u2tfrlwIGqOhO4EbhhoeOVJC3eUj2e2gq8UlV/OUebC4G7qur9qnoVmAS2JDkNOLGqHquqAm4HLhrqc1vbvxfYOj0LkSQdfUsVGuPAnUPfv5LkmSS7kqxutTXA60NtplptTdufXZ/Rp6oOAu8ApyzRmCVJR2jRoZHk54DfAf5bK90CfBrYCOwFvjHddET3mqM+V5/ZY9iRZCLJxL59+45g9JKkI7EUM43fAr5fVW8CVNWbVXWoqn4KfBPY0tpNAacP9VsLvNHqa0fUZ/RJsgo4Cdg/ewBVdWtVba6qzWNjY0twSZKkUZYiNC5l6NFUW6OY9gXgubb/ADDe3og6g8GC9xNVtRd4N8m5bb3iMuD+oT7b2/7FwCNt3UOStAxWLaZzkn8A/DPgS0Pl/5hkI4PHSK9NH6uq55PcA7wAHASurKpDrc8VwG7gBOChtgHsBO5IMslghjG+mPFKkhZnUaFRVT9h1sJ0VX1xjvbXA9ePqE8A54yovwdcspgxSpKWjr8IlyR1MzQkSd0MDUlSN0NDktTN0JAkdTM0JEndDA1JUjdDQ5LUzdCQJHUzNCRJ3QwNSVI3Q0OS1M3QkCR1MzQkSd0MDUlSN0NDktTN0JAkdTM0JEndDA1JUrdFhUaS15I8m+TpJBOtdnKSPUlebp+rh9pfnWQyyUtJLhiqb2rnmUxyU5K0+vFJ7m71x5OsX8x4JUmLsxQzjX9aVRuranP7fhXwnaraAHynfSfJWcA4cDawDbg5yXGtzy3ADmBD27a1+uXAgao6E7gRuGEJxitJWqAP4/HUhcBtbf824KKh+l1V9X5VvQpMAluSnAacWFWPVVUBt8/qM32ue4Gt07MQSdLRt9jQKODPkjyVZEerfaqq9gK0z0+2+hrg9aG+U622pu3Prs/oU1UHgXeAUxY5ZknSAq1aZP/PVtUbST4J7Enygznajpoh1Bz1ufrMPPEgsHYArFu3bu4RS5IWbFEzjap6o32+BXwb2AK82R450T7fas2ngNOHuq8F3mj1tSPqM/okWQWcBOwfMY5bq2pzVW0eGxtbzCVJkuaw4NBI8g+T/KPpfeB84DngAWB7a7YduL/tPwCMtzeizmCw4P1Ee4T1bpJz23rFZbP6TJ/rYuCRtu4hSVoGi3k89Sng221dehXwX6vqT5M8CdyT5HLgr4BLAKrq+ST3AC8AB4Erq+pQO9cVwG7gBOChtgHsBO5IMslghjG+iPFKkhZpwaFRVT8EfnVE/W1g62H6XA9cP6I+AZwzov4eLXQkScvPX4RLkroZGpKkboaGJKmboSFJ6mZoSJK6GRqSpG6GhiSpm6EhSepmaEiSuhkakqRuhoYkqZuhIUnqZmhIkroZGpKkboaGJKmboSFJ6mZoSJK6GRqSpG6GhiSp24JDI8npSb6b5MUkzyf5vVb/WpIfJXm6bZ8f6nN1kskkLyW5YKi+Kcmz7dhNSdLqxye5u9UfT7J+4ZcqSVqsxcw0DgK/X1W/ApwLXJnkrHbsxqra2LYHAdqxceBsYBtwc5LjWvtbgB3AhrZta/XLgQNVdSZwI3DDIsYrSVqkBYdGVe2tqu+3/XeBF4E1c3S5ELirqt6vqleBSWBLktOAE6vqsaoq4HbgoqE+t7X9e4Gt07MQSdLRtyRrGu2x0a8Bj7fSV5I8k2RXktWttgZ4fajbVKutafuz6zP6VNVB4B3glBH//o4kE0km9u3btxSXJEkaYdGhkeTngfuAr1bVjxk8avo0sBHYC3xjuumI7jVHfa4+MwtVt1bV5qraPDY2doRXIEnqtajQSPIJBoHxrar6I4CqerOqDlXVT4FvAlta8yng9KHua4E3Wn3tiPqMPklWAScB+xczZknSwi3m7akAO4EXq+oPhuqnDTX7AvBc238AGG9vRJ3BYMH7iaraC7yb5Nx2zsuA+4f6bG/7FwOPtHUPSdIyWLWIvp8Fvgg8m+TpVvt3wKVJNjJ4jPQa8CWAqno+yT3ACwzevLqyqg61flcAu4ETgIfaBoNQuiPJJIMZxvgixitJWqQFh0ZV/U9Grzk8OEef64HrR9QngHNG1N8DLlnoGCVJS8tfhEuSuhkakqRuhoYkqZuhIUnqZmhIkroZGpKkboaGJKmboSFJ6mZoSJK6GRqSpG6GhiSpm6EhSepmaEiSuhkakqRuhoYkqZuhIUnqZmhIkroZGpKkbsdEaCTZluSlJJNJrlru8UjSx9WKD40kxwH/Bfgt4Czg0iRnLe+oJOnjacWHBrAFmKyqH1bV/wHuAi5c5jFJ0sfSsRAaa4DXh75PtZok6ShbtdwD6JARtZrRINkB7Ghf/zbJSx/6qFaGU4G/Xu5B6IgcM/csNyz3CFaMj8s9+8c9jY6F0JgCTh/6vhZ4Y7hBVd0K3Ho0B7USJJmoqs3LPQ71854de7xnMx0Lj6eeBDYkOSPJzwHjwAPLPCZJ+lha8TONqjqY5CvAw8BxwK6qen6ZhyVJH0srPjQAqupB4MHlHscK9LF7JPcR4D079njPhqSq5m8lSRLHxpqGJGmFMDRWuCS7kryV5LnFtNHRkeT0JN9N8mKS55P83mHaec9WiCR/P8kTSf5Xu2f/4TDtvGcYGseC3cC2JWijo+Mg8PtV9SvAucCVh/mzN7vxnq0U7wO/UVW/CmwEtiU5d0S73XjPDI2VrqoeBfYvto2OjqraW1Xfb/vvAi8y4i8YeM9Wjhr42/b1E237wGKv92zA0JA+JEnWA78GPL68I9F8khyX5GngLWBPVXnPDsPQkD4ESX4euA/4alX9eLnHo7lV1aGq2sjgL05sSXLOco9ppTI0jjFtofXptn15ucejD0ryCQaB8a2q+iPv2bGjqv4G+B/AF7xnox0TP+7T/1dVrzNYrNMKlCTATuDFqvoD8J6tdEnGgP9bVX+T5ATgN4Ebquq6ZR7aiuRMY4VLcifwGPDLSaaSXL6QNjpqPgt8EfiNof+pfn52I+/ZinIa8N0kzzD4W3d7quq/z27kPRvwF+GSpG7ONCRJ3QwNSVI3Q0OS1M3QkCR1MzQkSd0MDUlSN0NDktTN0JAkdft/IS041zs20YsAAAAASUVORK5CYII=\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"# Diff Freq domain\n",
"plt.bar(['1-1','2-1','3-1'],[np.mean(np.square(H1.real-H1.real)),np.mean(np.square(H1.real-H2.real)),np.mean(np.square(H1.real-H3.real))])"
]
},
{
"cell_type": "code",
"execution_count": 15,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<BarContainer object of 3 artists>"
]
},
"execution_count": 15,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAY0AAAD8CAYAAACLrvgBAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvnQurowAAFjZJREFUeJzt3X+s3fV93/Hna3Z+kEZQA5eM2c7sDrcrWG1TXMdbtCmrO+w1VcwkkBwtxdosWUWsS6dOHV6loSWyBFtVNqSBZNWeDYsAy82G1ZalFiSLJlGTm18lhjKuSoZvceOb2aF0E05N3/vjfK527uHce7+c4/heh+dD+up8z/v7+Xz8OToSL77fz/fcb6oKSZK6+CtLPQFJ0uXD0JAkdWZoSJI6MzQkSZ0ZGpKkzgwNSVJnhoYkqTNDQ5LUmaEhSeps5VJP4GK79tpra926dUs9DUm6rHzlK1/5TlVNLNbuBy401q1bx+Tk5FJPQ5IuK0n+V5d2Xp6SJHVmaEiSOls0NJIcTHImyTcH6r+c5MUkJ5P827763iRT7di2vvrNSZ5rxx5IklZ/T5LHW/1EknV9fXYlealtuy7GB5Ykja7LmcYhYHt/IcnfA3YAP1FVNwG/0eo3AjuBm1qfB5OsaN0eAvYAG9o2O+Zu4FxV3QDcD9zXxroauAf4MLAZuCfJqpE+pSTpolg0NKrqS8DZgfKdwL1Vdb61OdPqO4DHqup8Vb0MTAGbk1wPXFlVz1TvAR4PA7f29Tnc9o8CW9tZyDbgeFWdrapzwHEGwkuSdGmNuqbxo8DfaZeT/nuSn2n11cCpvnbTrba67Q/W5/SpqgvAa8A1C4wlSVoio95yuxJYBWwBfgY4kuRHgAxpWwvUGbHPHEn20Lv0xQc/+MEFJy5JGt2oZxrTwOeq51ngL4FrW31tX7s1wKutvmZInf4+SVYCV9G7HDbfWG9RVfuralNVbZqYWPS3KZKkEY0aGv8V+FmAJD8KvBv4DnAM2NnuiFpPb8H72ao6DbyeZEtbr7gDeKKNdQyYvTPqNuDptu7xeeCWJKvaAvgtrSZJWiKLXp5K8ijwUeDaJNP07mg6CBxst+F+D9jV/kN/MskR4HngAnBXVb3ZhrqT3p1YVwBPtg3gAPBIkil6Zxg7AarqbJLPAF9u7T5dVYML8pIuM+vu/t2lnsIPrG/d+7Hv+7+xaGhU1SfmOfTJedrvA/YNqU8CG4fU3wBun2esg/QCSpK0DPiLcElSZ4aGJKkzQ0OS1JmhIUnqzNCQJHVmaEiSOjM0JEmdGRqSpM4MDUlSZ4aGJKkzQ0OS1JmhIUnqzNCQJHVmaEiSOjM0JEmdGRqSpM4MDUlSZ4aGJKmzRUMjycEkZ9rzwAeP/YskleTavtreJFNJXkyyra9+c5Ln2rEHkqTV35Pk8VY/kWRdX59dSV5q265xP6wkaTxdzjQOAdsHi0nWAn8feKWvdiOwE7ip9XkwyYp2+CFgD7ChbbNj7gbOVdUNwP3AfW2sq4F7gA8Dm4F7kqx6ex9PknQxLRoaVfUl4OyQQ/cDvwZUX20H8FhVna+ql4EpYHOS64Erq+qZqirgYeDWvj6H2/5RYGs7C9kGHK+qs1V1DjjOkPCSJF06I61pJPk48CdV9Y2BQ6uBU33vp1ttddsfrM/pU1UXgNeAaxYYa9h89iSZTDI5MzMzykeSJHXwtkMjyfuAXwf+9bDDQ2q1QH3UPnOLVfuralNVbZqYmBjWRJJ0EYxypvE3gPXAN5J8C1gDfDXJX6V3NrC2r+0a4NVWXzOkTn+fJCuBq+hdDptvLEnSEnnboVFVz1XVdVW1rqrW0fuP+09X1Z8Cx4Cd7Y6o9fQWvJ+tqtPA60m2tPWKO4An2pDHgNk7o24Dnm7rHp8Hbkmyqi2A39JqkqQlsnKxBkkeBT4KXJtkGrinqg4Ma1tVJ5McAZ4HLgB3VdWb7fCd9O7EugJ4sm0AB4BHkkzRO8PY2cY6m+QzwJdbu09X1bAFeb2Drbv7d5d6Cj+wvnXvx5Z6ClqGFg2NqvrEIsfXDbzfB+wb0m4S2Dik/gZw+zxjHwQOLjZHSdKl4S/CJUmdGRqSpM4MDUlSZ4aGJKkzQ0OS1JmhIUnqzNCQJHVmaEiSOjM0JEmdGRqSpM4MDUlSZ4aGJKkzQ0OS1JmhIUnqzNCQJHVmaEiSOjM0JEmdLRoaSQ4mOZPkm321f5fkj5L8YZL/kuSH+47tTTKV5MUk2/rqNyd5rh17oD0rnPY88cdb/USSdX19diV5qW2zzxGXJC2RLmcah4DtA7XjwMaq+gngfwJ7AZLcSO8Z3ze1Pg8mWdH6PATsATa0bXbM3cC5qroBuB+4r411NXAP8GFgM3BPklVv/yNKki6WRUOjqr4EnB2o/X5VXWhv/wBY0/Z3AI9V1fmqehmYAjYnuR64sqqeqaoCHgZu7etzuO0fBba2s5BtwPGqOltV5+gF1WB4SZIuoYuxpvFPgCfb/mrgVN+x6VZb3fYH63P6tCB6DbhmgbEkSUtkrNBI8uvABeCzs6UhzWqB+qh9BuexJ8lkksmZmZmFJy1JGtnIodEWpn8B+EftkhP0zgbW9jVbA7za6muG1Of0SbISuIre5bD5xnqLqtpfVZuqatPExMSoH0mStIiRQiPJduBfAh+vqv/bd+gYsLPdEbWe3oL3s1V1Gng9yZa2XnEH8ERfn9k7o24Dnm4h9HngliSr2gL4La0mSVoiKxdrkORR4KPAtUmm6d3RtBd4D3C83Tn7B1X1S1V1MskR4Hl6l63uqqo321B30rsT6wp6ayCz6yAHgEeSTNE7w9gJUFVnk3wG+HJr9+mqmrMgL0m6tBYNjar6xJDygQXa7wP2DalPAhuH1N8Abp9nrIPAwcXmKEm6NPxFuCSpM0NDktSZoSFJ6szQkCR1ZmhIkjozNCRJnRkakqTODA1JUmeGhiSpM0NDktSZoSFJ6szQkCR1ZmhIkjozNCRJnRkakqTODA1JUmeGhiSps0VDI8nBJGeSfLOvdnWS40leaq+r+o7tTTKV5MUk2/rqNyd5rh17oD0rnPY88cdb/USSdX19drV/46Uks88RlyQtkS5nGoeA7QO1u4GnqmoD8FR7T5Ib6T3j+6bW58EkK1qfh4A9wIa2zY65GzhXVTcA9wP3tbGupvc88g8Dm4F7+sNJknTpLRoaVfUl4OxAeQdwuO0fBm7tqz9WVeer6mVgCtic5Hrgyqp6pqoKeHigz+xYR4Gt7SxkG3C8qs5W1TngOG8NL0nSJTTqmsYHquo0QHu9rtVXA6f62k232uq2P1if06eqLgCvAdcsMJYkaYlc7IXwDKnVAvVR+8z9R5M9SSaTTM7MzHSaqCTp7Rs1NL7dLjnRXs+0+jSwtq/dGuDVVl8zpD6nT5KVwFX0LofNN9ZbVNX+qtpUVZsmJiZG/EiSpMWMGhrHgNm7mXYBT/TVd7Y7otbTW/B+tl3Cej3JlrZeccdAn9mxbgOebusenwduSbKqLYDf0mqSpCWycrEGSR4FPgpcm2Sa3h1N9wJHkuwGXgFuB6iqk0mOAM8DF4C7qurNNtSd9O7EugJ4sm0AB4BHkkzRO8PY2cY6m+QzwJdbu09X1eCCvCTpElo0NKrqE/Mc2jpP+33AviH1SWDjkPobtNAZcuwgcHCxOUqSLg1/ES5J6szQkCR1ZmhIkjozNCRJnRkakqTODA1JUmeGhiSpM0NDktSZoSFJ6szQkCR1ZmhIkjozNCRJnRkakqTODA1JUmeGhiSpM0NDktSZoSFJ6mys0Ejyz5OcTPLNJI8meW+Sq5McT/JSe13V135vkqkkLybZ1le/Oclz7dgD7TnitGeNP97qJ5KsG2e+kqTxjBwaSVYD/wzYVFUbgRX0nu99N/BUVW0AnmrvSXJjO34TsB14MMmKNtxDwB5gQ9u2t/pu4FxV3QDcD9w36nwlSeMb9/LUSuCKJCuB9wGvAjuAw+34YeDWtr8DeKyqzlfVy8AUsDnJ9cCVVfVMVRXw8ECf2bGOAltnz0IkSZfeyKFRVX8C/AbwCnAaeK2qfh/4QFWdbm1OA9e1LquBU31DTLfa6rY/WJ/Tp6ouAK8B14w6Z0nSeMa5PLWK3pnAeuCvAT+U5JMLdRlSqwXqC/UZnMueJJNJJmdmZhaeuCRpZONcnvo54OWqmqmqvwA+B/xt4NvtkhPt9UxrPw2s7eu/ht7lrOm2P1if06ddArsKODs4karaX1WbqmrTxMTEGB9JkrSQcULjFWBLkve1dYatwAvAMWBXa7MLeKLtHwN2tjui1tNb8H62XcJ6PcmWNs4dA31mx7oNeLqte0iSlsDKUTtW1YkkR4GvAheArwH7gfcDR5Lsphcst7f2J5McAZ5v7e+qqjfbcHcCh4ArgCfbBnAAeCTJFL0zjJ2jzleSNL6RQwOgqu4B7hkon6d31jGs/T5g35D6JLBxSP0NWuhIkpaevwiXJHVmaEiSOjM0JEmdGRqSpM4MDUlSZ4aGJKkzQ0OS1JmhIUnqzNCQJHVmaEiSOjM0JEmdGRqSpM4MDUlSZ4aGJKkzQ0OS1JmhIUnqzNCQJHU2Vmgk+eEkR5P8UZIXkvytJFcnOZ7kpfa6qq/93iRTSV5Msq2vfnOS59qxB9qzwmnPE3+81U8kWTfOfCVJ4xn3TOM/AP+tqv4m8JPAC8DdwFNVtQF4qr0nyY30nvF9E7AdeDDJijbOQ8AeYEPbtrf6buBcVd0A3A/cN+Z8JUljGDk0klwJ/F3gAEBVfa+qvgvsAA63ZoeBW9v+DuCxqjpfVS8DU8DmJNcDV1bVM1VVwMMDfWbHOgpsnT0LkSRdeuOcafwIMAP8pyRfS/JbSX4I+EBVnQZor9e19quBU339p1ttddsfrM/pU1UXgNeAa8aYsyRpDOOExkrgp4GHqupDwP+hXYqax7AzhFqgvlCfuQMne5JMJpmcmZlZeNaSpJGNExrTwHRVnWjvj9ILkW+3S0601zN97df29V8DvNrqa4bU5/RJshK4Cjg7OJGq2l9Vm6pq08TExBgfSZK0kJFDo6r+FDiV5MdaaSvwPHAM2NVqu4An2v4xYGe7I2o9vQXvZ9slrNeTbGnrFXcM9Jkd6zbg6bbuIUlaAivH7P/LwGeTvBv4Y+Af0wuiI0l2A68AtwNU1ckkR+gFywXgrqp6s41zJ3AIuAJ4sm3QW2R/JMkUvTOMnWPOV5I0hrFCo6q+DmwacmjrPO33AfuG1CeBjUPqb9BCR5K09PxFuCSpM0NDktSZoSFJ6szQkCR1ZmhIkjozNCRJnRkakqTODA1JUmeGhiSpM0NDktSZoSFJ6szQkCR1ZmhIkjozNCRJnRkakqTODA1JUmeGhiSpM0NDktTZ2KGRZEWSryX5nfb+6iTHk7zUXlf1td2bZCrJi0m29dVvTvJcO/ZAkrT6e5I83uonkqwbd76SpNFdjDONTwEv9L2/G3iqqjYAT7X3JLkR2AncBGwHHkyyovV5CNgDbGjb9lbfDZyrqhuA+4H7LsJ8JUkjGis0kqwBPgb8Vl95B3C47R8Gbu2rP1ZV56vqZWAK2JzkeuDKqnqmqgp4eKDP7FhHga2zZyGSpEtv3DONfw/8GvCXfbUPVNVpgPZ6XauvBk71tZtutdVtf7A+p09VXQBeA64ZnESSPUkmk0zOzMyM+ZEkSfMZOTSS/AJwpqq+0rXLkFotUF+oz9xC1f6q2lRVmyYmJjpOR5L0dq0co+9HgI8n+XngvcCVSf4z8O0k11fV6Xbp6UxrPw2s7eu/Bni11dcMqff3mU6yErgKODvGnCVJYxj5TKOq9lbVmqpaR2+B++mq+iRwDNjVmu0Cnmj7x4Cd7Y6o9fQWvJ9tl7BeT7KlrVfcMdBndqzb2r/xljMNSdKlMc6ZxnzuBY4k2Q28AtwOUFUnkxwBngcuAHdV1Zutz53AIeAK4Mm2ARwAHkkyRe8MY+f3Yb6SpI4uSmhU1ReBL7b9/w1snafdPmDfkPoksHFI/Q1a6EiSlp6/CJckdWZoSJI6MzQkSZ0ZGpKkzgwNSVJnhoYkqTNDQ5LUmaEhSerM0JAkdWZoSJI6MzQkSZ0ZGpKkzgwNSVJnhoYkqTNDQ5LUmaEhSerM0JAkdTZyaCRZm+QLSV5IcjLJp1r96iTHk7zUXlf19dmbZCrJi0m29dVvTvJcO/ZAe1Y47Xnij7f6iSTrRv+okqRxjXOmcQH41ar6cWALcFeSG4G7gaeqagPwVHtPO7YTuAnYDjyYZEUb6yFgD7ChbdtbfTdwrqpuAO4H7htjvpKkMY0cGlV1uqq+2vZfB14AVgM7gMOt2WHg1ra/A3isqs5X1cvAFLA5yfXAlVX1TFUV8PBAn9mxjgJbZ89CJEmX3kVZ02iXjT4EnAA+UFWnoRcswHWt2WrgVF+36VZb3fYH63P6VNUF4DXgmosxZ0nS2zd2aCR5P/DbwK9U1Z8t1HRIrRaoL9RncA57kkwmmZyZmVlsypKkEY0VGkneRS8wPltVn2vlb7dLTrTXM60+Dazt674GeLXV1wypz+mTZCVwFXB2cB5Vtb+qNlXVpomJiXE+kiRpAePcPRXgAPBCVf1m36FjwK62vwt4oq++s90RtZ7egvez7RLW60m2tDHvGOgzO9ZtwNNt3UOStARWjtH3I8AvAs8l+Xqr/SvgXuBIkt3AK8DtAFV1MskR4Hl6d17dVVVvtn53AoeAK4An2wa9UHokyRS9M4ydY8xXkjSmkUOjqv4Hw9ccALbO02cfsG9IfRLYOKT+Bi10JElLz1+ES5I6MzQkSZ0ZGpKkzgwNSVJnhoYkqTNDQ5LUmaEhSerM0JAkdWZoSJI6MzQkSZ0ZGpKkzgwNSVJnhoYkqTNDQ5LUmaEhSerM0JAkdWZoSJI6uyxCI8n2JC8mmUpy91LPR5LeqZZ9aCRZAfxH4B8ANwKfSHLj0s5Kkt6Zln1oAJuBqar646r6HvAYsGOJ5yRJ70iXQ2isBk71vZ9uNUnSJbZyqSfQQYbUak6DZA+wp7398yQvft9ntTxcC3xnqSeht+Wy+c5y31LPYNl4p3xnf71Lo8shNKaBtX3v1wCv9jeoqv3A/ks5qeUgyWRVbVrqeag7v7PLj9/ZXJfD5akvAxuSrE/ybmAncGyJ5yRJ70jL/kyjqi4k+afA54EVwMGqOrnE05Kkd6RlHxoAVfV7wO8t9TyWoXfcJbkfAH5nlx+/sz6pqsVbSZLE5bGmIUlaJgyNZS7JwSRnknxznDa6NJKsTfKFJC8kOZnkU/O08ztbJpK8N8mzSb7RvrN/M087vzMMjcvBIWD7RWijS+MC8KtV9ePAFuCuef7szSH8zpaL88DPVtVPAj8FbE+yZUi7Q/idGRrLXVV9CTg7bhtdGlV1uqq+2vZfB15gyF8w8DtbPqrnz9vbd7XtLYu9fmc9hob0fZJkHfAh4MTSzkSLSbIiydeBM8DxqvI7m4ehIX0fJHk/8NvAr1TVny31fLSwqnqzqn6K3l+c2Jxk41LPabkyNC4zbaH16237paWej94qybvoBcZnq+pzfmeXj6r6LvBF4B/6nQ13Wfy4T/9fVZ2it1inZShJgAPAC1X1m+B3ttwlmQD+oqq+m+QK4OeA+6rqM0s8tWXJM41lLsmjwDPAjyWZTrJ7lDa6ZD4C/CLws33/p/rzg438zpaV64EvJPlDen/r7nhV/c5gI7+zHn8RLknqzDMNSVJnhoYkqTNDQ5LUmaEhSerM0JAkdWZoSJI6MzQkSZ0ZGpKkzv4faJbYRywYtNsAAAAASUVORK5CYII=\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"# Diff Phase domain\n",
"plt.bar(['1-1','2-1','3-1'],[np.mean(np.square(H1.imag-H1.imag)),np.mean(np.square(H1.imag-H2.imag)),np.mean(np.square(H1.imag-H3.imag))])"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": []
}
],
"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.7.6"
}
},
"nbformat": 4,
"nbformat_minor": 4
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment