Skip to content

Instantly share code, notes, and snippets.

@Gijs-Koot
Created February 21, 2020 10:31
Show Gist options
  • Save Gijs-Koot/44efc106f229442f1401bab23f3dc733 to your computer and use it in GitHub Desktop.
Save Gijs-Koot/44efc106f229442f1401bab23f3dc733 to your computer and use it in GitHub Desktop.
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 2,
"metadata": {},
"outputs": [],
"source": [
"import numpy as np\n",
"import matplotlib.pyplot as plt"
]
},
{
"cell_type": "code",
"execution_count": 47,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.image.AxesImage at 0x7faf81dfb8d0>"
]
},
"execution_count": 47,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAPsAAAD6CAYAAABnLjEDAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAgAElEQVR4nO2deZyNdfvHP1+zMwszZsZYx75kX6JoI1kS9UuJrCmhDSVSiTYeW1QeQpZItqyhSApZk53s21gGs5lhZsyM7++POZ6X43PNkwcz6L7er5cX83HdZ+77Puc69zmf+/pel7HWQlGUfz65bvcOKIqSM2iyK4pD0GRXFIegya4oDkGTXVEcgia7ojiEm0p2Y0xjY8xeY8wBY0zfW7VTiqLcesyN3mc3xngA2AegIYAoAJsAtLbW7s5qm6BgTxtWyMtNO3YhWIzNdZHfh4qFnRFjj54JE3ZQ3gfvM8mkFaiQJMYmXfYlLTYmgDT/4Ivi9iGe/Lh5stivA7uDSLO+XkIkkO7L58Z6yo97WXiI3LlT5dijHqSlhLHm5ZMubp+WzL/MJ/6yvGMF+TEK+cSLoUdP8fPrFXJJjPXOxY+bHJ2btPCCseL2Fy97k5YLco7EnePXwmV/+Xi9zvITf7lAhhibfkk45+flF05anmu2jYtFRtIFMTiLl8h1cS+AA9baQwBgjJkBoAWALJM9rJAXRi0o6aa9sq6NGJt7ux9pY7t9KcZ2/+JV0qQXOQAUHrODtD4Lfxdj11woQ9qsKfVJq/fcn+L2z4esJa2Oj7xfT1RvTFpa6YJibFw5Pjcp+eUXw8WC/IKqXu2gHNs1hLQ9r/GbUMHIc+L20bs4KYvPl99YMIAfY3CJ78XQlz95g7QCbY+IsUXzxJG2dXhV0t78cLq4/eYLkaT5e8jHMHtcA9JSHkoUYwuO5TeRC28miLHRx/ORVmiZ/CH8VF335/3k8JFiHHBzH+MLATh+1c9RLk1RlDuQbDfojDFdjDF/GGP+SIiVP7YoipL93EyynwBQ5KqfC7s0N6y146y1Na21NYOC+buIoig5w80YdJ7INOgaIDPJNwFoY63dldU2xSoG2D5zarhpD+Q+IMZmCA7b0sRKYuySfo+Q1nLwT2Ls7KjqpCXNLyDGep/nc/P1J5+RFpRL/sTSbu/zpLUq9IcYu6BWJGmH+lUWY0vM5u96B5/j79YAkBbE+1Zgtfwe32fgNNJ6rn6ONL/D/P0TAC6VZ6MyX6BsXl5Yn5+01LJsngLApLqTSOu+lc8tAHit4PMgPY9+7U6J2x89xvtVeIl8kVr6+SjSKi9gfwEAxjaaSNqw9rJfdaALW2n+O2Szp8isY24/rz09HQmp0bfWoLPWphtjXgXwEwAPABP/W6IrinJ7uRk3HtbaJQCW3KJ9URQlG9EKOkVxCJrsiuIQNNkVxSHcsBt/I+QpE2Hv+byjmxbySpoYu/u9UNK2NvpCjK01tRdpBTbIDnlUAzYqyw87KcYe6FKYNKlscXq3EeL2z4/m/Sre7JAYu2NXUdKCt8ku8IWCvA/9Ws0SY0f8+1nSEqrKpaaVSkaRtvMIV/H57eUyYgDIfZpfS8Vf2CfGHhtTmrSQtbJDnmcqlx1/WWyhGPt861dIs558TTtbRT4G3zg+hpIv/yXGbl9UnrRLgXI+pRfmKjwTI9/V+LbFaNLaLOLjAoDCK93Lc7euHIWkuCjRjdcru6I4BE12RXEImuyK4hA02RXFIeSoQRfoX8jeW6Wbm7a/k2xSlJjJ64IrDdkmxvYL+420utPfkh/3nY2kld8kLw/dOKQmaTGVOLbo8hRx+9LD9pBWzDdGjK2dm5edlvU6L8Y2HvE2aeEbLoixwcOOk/Zh4UVibIplQ/CVXlz+6b9Efh7qbOD9nbmfy5MBIPIlWkYBkzdQjD3QmU3Cd1rKy2EvCQv757bhcup6UzaL28/+mpcwd3lZPl+j5jcjLT1AXs9eciYbdLkGykuFD55mc3rm/V+JsZ2+6OG+7bQRSD59XA06RXEymuyK4hA02RXFIWiyK4pD0GRXFIdwU0tc/1dsLoMMv2t+ZRY3A1JCuGNkw6CdYmyDz3uTNqO73HivTRq7y2cmym58Ql0uufWO4/fHo03k0sukEdzoMFevLWLs+O31SPM4LD9uxH4uMd7fXm5uUL4DN3P417RGYuyBjyuQVmMgu9b7N0eI23sZbvYYNNdfjP1rQFnSwsucFWMjH+c7KMMSW4qxhZdJjT3Y5Z8bvF3cft6j3DCkU6DcoPNz4bUbulG+fh5rxE1CQ8ZzOTYA/DJkGGkduvQUY1M6upcS23lZdPOFXtkVxTFosiuKQ9BkVxSHcFPf2Y0xRwAkAsgAkG6t5ZIzRVHuCG6FQfeItVau+7uG9LDLOPu6e7fRUgHy6KXhjWaT9uFxLk8EAE+hgWnvl7uLsREePB7I74Q8xcN68GiqhPK8/ZzG8qSarhXbklbMTz5VBx7hDqo1f+0mRAK+P7HJV369XGqaMYtHH+0bco8Y+8nnXJI56KnWpBkfeT38mnrhpAUmbhBjg2azAXusr3ytyFeKtdBHudwWAFK281Sa8I1sWtU7LHeBTc/NZu3Qomy0AkDukxwbsiFajA08xs/P+5P4OQeAIdE8aabmp3JX4osj67j9fCoh6w/r+jFeURzCzSa7BbDMGLPZGNPlVuyQoijZw81+jK9nrT1hjAkDsNwY85e1dtXVAa43gS4A4BUqDzJQFCX7uakru7X2hOvvMwDmIXOy67Ux/xn/5BHI3x8VRckZbjjZjTF5jDEBV/4N4DEAcomboii3nZv5GB8OYJ4x5srjTLfW/vjfNsjtlYbKYe4dRH/fW1KMbT+fO7NmhVcTdrh9f5bLBi+W4tnXMHK5bGIkP4ZvNJ+ydl/3IA0AHv+/daSV9ZE7qDZo15m0rybwLDEA+Ok1nnm3vLfsGEfm2Uvato55xNg5sbVIi6mWl7RzNeRzm1eYz55/mzzr7Xx/vgtTKrfceTf1Ey6jPXqqihgb/nosaQlr+C5ByG65+/Alf77+bWgglwe3+ZXnCU5rSB9uAQCFOh/lfcglz7bb90Y50oZNl5tXlH7H3f0f9Ee8GAfc3Ky3QwDkM64oyh2H3npTFIegya4oDkGTXVEcQo6uZw/xTESn8NVu2tp9skGX9yCv2c6zXR7TtKcmGzABxeT3MZ+z3Am27MT9YmzioBqk+Z1lU+W7aXK57L0L2WTcNaqEGPvlMh5t9eS0N8XYPJXZhGo5dLUQCUz+gbulVqx3QIw9ciGEtDj2ihBWUu6QGzyGz3nwhDNirH86l8seTRDMUwCPbj5N2pkJ8lr/UW+xkeVRjheer75YRt5+/aOktejNY7EAYP6HDUkLeEE+3j1D+XkfGc2/C8js+3AtXXZz6TUABL/pHht/SB63BeiVXVEcgya7ojgETXZFcQia7IriEDTZFcUh5Kgb72cyUMHLvfvn05X/FGOfe4ibHrzVTR5Iv7k+l5U+NUNuThBXnrudnkuVO6CuGMUu+1MNnyftz1QuKQWAsHXsqh76WF4M9MTvfGy5Y+Uy3paRW0lrFijPX1uy+WHStuSPFGMrDGInOb0fN+vwHsdNPQDg8kHehynFuGQYAEqveJG05vfIHV//jC1CWp7TcrnrwAbPkHakFc+Ki/zmiLi916t8l2DlWdm5P1mfy4YfD5Gbavg3O0ba8UrcYRcAesyeQdrOZD4HADC258NuP6d8zPP6rqBXdkVxCJrsiuIQNNkVxSFosiuKQ8hRg+5gSgie2d3eTascIpfAPreOW9qVWCabeZMSKpL2yODfxdhpSx8irUuBX8XY5w81Ie1MXS4pHd6WO7ACQPmRu0hbtV7u7OqZyGbcxZryWvAVPXlUVMuJ8rnxjeGy49r3HBFjo2qUJu23RsNJeyhD7jXQcSB36W1e7ykxtlwMl+z6r0oVY4+d4zLaSw3kuWF+QmzZJlwOnTpTXtP/8/NDSXtoudyvAF5s0C3Zxa9FACifn0tuX5s9R4ztt/tJ0jy/l03RwHzur5szybKpC+iVXVEcgya7ojgETXZFcQia7IriEIy1WQxIvxJgzEQAzQCcsdZWdGnBAGYCiARwBMCz1loezn0NxSv62w/mujdL/OIAr7cGgPhdbISV+Upu1vjGssWk9fz6JTE2LZCPN72gbAwhiaupyvbeQVrhlbIpsncQm3FRT8qVX9VLckPCThFrxNjR1bmp4aUawowkAMmh3qQF7ZKfqtFLviatyZTepOWrJa/ZTvid+wosfmmIGNv63bdIe6rPCjG2mt8R0l6ZyxV4ABC5iPsVeJ3mme3xNbg5JgAEvMRGWqfCstkb6nmetK5z5FkpRVawUZrhI19rT7bm8VpbHhwrxj7bwH2d+7ojk5GQfEp8QV7PlX0ygMbXaH0BrLDWlgawwvWzoih3MH+b7K4JL9e2RmkBYIrr31MA8L0CRVHuKG70O3u4tfbKZ+rTyOwhL2KM6WKM+cMY80diHH+UURQlZ7hpg85mfunP8ov/1eOfAvLxd2BFUXKGG032aGNMBAC4/pYdG0VR7hj+1o0HAGNMJIAfrnLjhwKIsdYONsb0BRBsrX377x7nnsredsYP7i5ou8Fy6eXsd7hssfF37AwDQDBXpSKuvLwPGb58vJ0e/VWMnTaP7xTk38Fu+sf/Gi9uP6T2IywGy2vfjw3mbqkRI+RPQjGV/Ejr8cZsMXZK9xaknY9khx4AQldxF1eTxuvZ/3qjkLh9xFo+t3Gl5fXVF8vwHRAvP/lrnsdu7jfQsPkmMdYnF+9vmDe75gcuym78RxE/kzYhrroY++OpCqQFvCgfw1+DQkn7+QHuKAwAjYTX+SP1uYeBxPz2i3F2d8yNufHGmO8ArANQ1hgTZYzpDGAwgIbGmP0AHnX9rCjKHczfLoSx1sqrPIAGt3hfFEXJRrSCTlEcgia7ojiE6zLobhW+hYrYIt17ummL2g8TY9v349FHIz8aLcZ2+PZV0kpNkEf2hM3k+dVbp/K8cwD46i1uZNlhMjey9Koml5/+VJ2Nu7orXxdjy78bTdqFibJBt6DCd6SNjOFRVQDwQehuUZdoVJBnvO8by6W55cbybHUAODGAX0u5ZweJsQMGTiLtje9eEGNLTRLKpC/LM+InrJpO2sgYXv+/bPz94vb/6sXP2cgW/yfGPjSDewh8teUBMdbruA9pDR7bIsYW9OHXaLMAuaHoa73cX0/bfhmFpLjjN1wuqyjKPwBNdkVxCJrsiuIQNNkVxSFosiuKQ8hRNz7IJ9zeX9B9fNLJz+Uun37Tuaz09APyvpbts5O0w5NKiLGXD3HpZXDls2JsqbznSNsTw2WWS6uyswwA85O4W+sXk+TVwKnVLpDmuVM+N888+xtpYV5cEgoAX+1jJ7pQf/k8Hn6aO5gKxjBqtpLHNEXdn0xaWn12+AEgrgy700++/KsYeyyZ92vtkspibPHPuXa6/Qbe3ylPNxK3Ty4SQFpcGfmuyMXa3P3Xe4c83kuo2EXB7w+Ksbs/LEpa2Gq5/u3hXu7jtaa2WYHTu+S5YXplVxSHoMmuKA5Bk11RHIImu6I4hBw16AKCCtvq97/mprX6bKkYu/Cp+0h7ZfEPYmz/oZ1IC35WLpf9pjSXmvY7yWOeAODgIF4UnxzC67N9zsulmwH7uKvp6QfkMT5Pv/wLaRNX8agqAMi7i9+jI344LsYWn8vm40/75cX+rSpsJm36Gi4rLbpEPl7veO6KumzOFCESiErnktsHfuopRAIBe9ggC24sjw3ze4qPt8AKPl9rf5bHNA1txft76JK89n3Ra7zw02fbYTH22PgI0oJmsBkIAAEHeYzWgecCxVjfc+5e3OHJI5B8SstlFcXRaLIrikPQZFcUh6DJrigO4Xp60E00xpwxxuy8ShtgjDlhjNnq+tM0e3dTUZSb5W970CFz/NOXAL65Rv/MWit3nsiCAkVj0Gf0VDet+9rnxdjihdhQlGacAUDxxftJe70QdwkFgN5RzUg7MVCekxaw/RhpCa24DDdGbj6Kgi+yOz2qzL/F2I5ruHFD09pyR9HlQrfTyzHXDu3J5MGgv0g7MFw+3tUF+Q5Ixbd4Bt3FOQXF7Q9252vHI53kmWx+x/hORb5/ySW/cRncAKNWIJcyA8CeZtyI5HKd9aSlj5DvQjXPwyWwdQZcO/0sE+8+3HDE84PCYmzR93gGXWwWr5vYiuy8V79vnxj7ZbGFbj83WiKfF+DGxz8pinKXcTPf2V81xmx3fczPd8v2SFGUbOFGk30MgJIAqgI4BWB4VoFXz3pLiJXHFSuKkv3cULJba6OttRnW2ssAxgOQv0zDfdZbULA8HURRlOznRsc/RVyZ4mqM6QmgtrX2ub97HP/gIrZiox5uWnIbYcE0gClVJpNWIgs78RdhvfOIw4+JsSfO8Tr5gFU8TgkAUvKzSTig/bekTW4h/y7PsVz2mHiJ13EDwMn1bHoVWi2PEur7by7p/OZMXTH23CPCGvO6cqnoM1/+SNrclxqSVm6UMG8LwLp/1yRt/PsjxdinlnNH4IAwuWvtkIpzSVsSL69nbxPMZtyHtfgYTG553fnJJ3gteXw1+XkI2s5lvFvfkQ3Y4j91Jq38p3JX4o+WzSDtvSc7iLETF7l3w23a9By2bU8Ty2X/1o13jX96GEB+Y0wUgA8APGyMqYrM6a1HALz8d4+jKMrt5UbHP32dDfuiKEo2ohV0iuIQNNkVxSFosiuKQ7iectlbRnpu4FxVd6Pwm0rsbgPAqz14JlpiJ7mcMmASlxeW7ivPOLv0dQHSkoTSXABILs7lrifTuH7o0EBfcfu840NICzzMZZMA4H8Pa6l55afnrTEvkRZ0SK5hOD+Tz1m3Muy6A8DME+ym+yalkrZqSi1x+8vN2V3u95g88TtX/3TSgr/izr8A0LvL06SZNXxXBQA2NihG2sKtfPei9mK5UUaxRey8x1cTQxG6ne909ImWu+lGzuTr6olm/FoEgFf2tCHNu7jc6GJtivtdnCQr5wigV3ZFcQya7IriEDTZFcUhaLIrikPI0e6ylSp72bmL87tpzce+LcYW+pXHIZ2tJo9DKv/8HtLK+vNaYwBYV5NNoNp/8O8CgE2PFSJtb19ez1773r3i9k+G/kla/xlsvgBA6xY80mltVbm09oeoP0ibdL6IGDu13xOkRT0mP+d1q/Ca6dhULivdGxUubv9KNT6Gbz+Xxyzln7CRtEOD5CUW4Zu4m+2MEfLaqxef4mLOvS/y62bgI1yCCwBfDn6GtIvhsoHrG8PnMXgXr4cHAK/Dp0nLiJNLxaNfrEFawcVyt+SQ79wfY2nHBYjZc067yyqKk9FkVxSHoMmuKA5Bk11RHIImu6I4hBx140tUymM/nuteFzq2Z0sx9tiz7MDaFLnTTZMa20nbmyDP58KgUJLOVJdd78tCtWrXtotJG7mlvrj9g6UOkObnITdCWL6/HGklwuVOoQXzcGfWjwrKM/MeWNKLtMj58nMeX5KbMRRYIzRY2M8dZwHAZymXdNYN4XMAAL+dK0Parj3yHQXfaH4i0n3lY8glnN7AQ6yN7y831TiSziXOHw5vJ8aO7zOKtPcbthJjTzXmWW8F1shu/NIl00m7v1dXMdYj1f08bPtlFJLidNabojgaTXZFcQia7IriEK5n/FMRY8xKY8xuY8wuY8wbLj3YGLPcGLPf9bf2jleUO5i/NeiMMREAIqy1fxpjAgBsBvAkgI4AYq21g40xfQHks9b2+W+PFeQXYe8r6T7m6ExdNkQAIOxXLnc90Eku00zLy2u5H6jGY48AoEPY76RtS+Y10ABQ0fc4aW/u4HLKgJm8nh4AghZzF9YRO34SYxcmViFt6SlhkTuA1Aw2KqMP5RcigbxF2ATKtYi78QLA3PeHkta1Io/xS6tWUtx+5BTurNrsxzfE2ND1fAz5154RY59byGW4qxPKirGN8+0gbfBgHjEWPIm70AJAox1sfi47U16M7VGUR4x91lY26A69ztdVe0rug9DiYS4lfjFkjRjbq5y7Obw+ZQkSLsfcmEFnrT1lrf3T9e9EAHsAFALQAsCVrgBTkPkGoCjKHcr/9J3d1T++GoANAMKv9I4HcBqAfNlVFOWO4LqT3RjjD+B7AD2sde99YzO/C4jfB64e/3QpQ14RpChK9nNdyW6M8UJmon9rrb2yNjDa9X3+yvd68QvX1eOfvD3kKRyKomQ/12PQGWR+J4+11va4Sh8KIOYqgy7YWisvTncREFTY1rjvNTet2iBe8w0Anrm4gs7LyE0V6/nzOuyRrZ8VY3MP53XFKRlcOQYA8V/xKKDYCux9FJ8jj/EpO4mrx3a8zUYcAJy+l6v4IqezQQgA09bMJG1tqmy69R/SibQCP58SIoHo+lzlFbqRj80kcyNOALg4mp+z6N/lWe6pIfxcRi7kJpQAENWZ9dJvx4ixu9/jY8grjGlKkf1MpIbyfuXbIV8TY6tw7KAGs8XYIXt5RFiBd+R9OPEJ/77gMXIvh5iK3m4/H5w6Asmn5Qq66+kuWxdAOwA7jDFbXVo/AIMBzDLGdAZwFICcXYqi3BFcz/inNQDkVh1Ag1u7O4qiZBdaQacoDkGTXVEcgia7ojiEHB3/lBGWgfhXE920pYcryMFbuQS15uM7xdBua3lQfccJcnnhfXn2k9Z7p7ymPvEBvlNReBk7zmkh8i3Fk8lBpOXqJ5eE+s4sTFqrn9aJsS33Pkda7AV5H1IfSyRtZv8ZYmyVld1IG9h3Pmmjn+COtQDQr8Qi3q9i8kinb2pXZtFTfjnmWcOlsa+tlI/hg4EvkJYqrNpo1JxLUgFgwWae9WTS5Wti+dFcWvtekFxIWkK4U2G95N4GKbv4dWOsHBvc6KTbz0cXyHGAXtkVxTFosiuKQ9BkVxSHoMmuKA4hRw06E+cJr9nuZZ1N35SNNN8ybDQkZsjrfz2TeG302UvyPOvPOvB69MsfiKHwOcOP679LKDWNl2di7/yRjSUr98xEiYUHSRvhm0VRolDhbBvIJbu/1x5H2rP15QaKXeeuIu3TA7yePeDQMXH7z6rVIS29irz2PaoHG4rhm2Rz6YFOm0ir4i2Xy8Y341Fel2L4dbPvEdnQLPog149F1c+iueX/sen2aWl5rNS7T/LYL1NUHjtWYii/npKKy0Znyiz38uCMOLn0G9Aru6I4Bk12RXEImuyK4hA02RXFIWiyK4pDyFE3Pj03cLaWu4MZnSp3Zm0Xyl1gl52vJMaWGc2D6ks9zt1pAWBlUx50v6HWCDG28TTujPrpci7TfKWP3EH1+xeHkdZy80tibGq5QqR5JckusGETGB3KyHc1dqexE93xB+6KCgDjO/8faa9N4HLZgF3J4vaja91HWtSjsustjW+6ECa/HKv787ipkxneQiTweqWVpE07ei9ptqjcVKPuR9x19s86fmJslfWppM0/W12MzfDlJ81nt+ywpwXx46b7yqvMJ7zrPsaq7Xr5dQ/olV1RHIMmu6I4BE12RXEINzP+aYAx5oQxZqvrD5daKYpyx3Az45+eBZBkrWUXKgtKV/KzIxaUctNeX99ajPXwYEPD00vuLostbPIVn3xEDP1wDRtO7cf3ECKBwg24LPTwWe7iGviz3PnzQiE2VXzPiqF4uusvpP3WnctPAeDwk2y6hcvLs3EhnN/Pc6XLz/nbPdh8HLSnMWmFOp4kDQCOTeQ1+QU+5665ADB96hek1fud19MDgD3K5zfvPXK57Kbqs0j79ByXLadlUbe88fHipN23hEuZAeDb7+uTlhIhd8jd0exz0lrVeVqMPf4lv55X15woxj7ew/21u/3nUUiKvcHusq6pL6dc/040xlwZ/6Qoyl3EzYx/AoBXjTHbjTETdYqrotzZ3Mz4pzEASgKoiswr//AstvvP+KeE2Cw+hiuKku3c8Pgna220tTbDWnsZwHgAXLkA9/FPQcFZrO9UFCXbuR433gD4GsAea+2Iq/SrF9I+BUDuBqkoyh3B9bjx9QCsBrADwBWLvB+A1sj8CG8BHAHw8lUjnEV8ihW2Bd51Ly0tN4a7nwLAB/OmkVbLRy4ZfHT3U6Rd+FYuh5QaP+T/nee/AUCl2YdJO3whhLRNe0qI25f/jLuPjvtRdlV3XuLHfXdPCzHWczbHXoiQz03F5n+RlthKLv881roYae+/9C1pJb3kWwr9Wr9I2sk+sjvdtexq0kYubCbGlp7Evy+hijys7WIYX7+CmvPdA19PuVHGj+UWk9akhHxX5Ozz3Il21Lujxdj3u3YhzfOCfG6SivDdltNN5fl6pca4fzXeuG0MzieduGE3PqvxT0v+bltFUe4ctIJOURyCJruiOARNdkVxCH9r0N1K/EoVtKVGuJs4AVPl9ewxlfg2XZGf5G6c+7uw9VCm81YhEvj6yG+ktXrrLTH2kj9bFdK64iV9h4jbP9PzTdKsh2yk+cawYZTuJ9+qrDJgC2nr/l1TjL0Yzr8vJVRYEA8gsGQ8aeEtD5HWZgdrADBgGZd/vvAwn28AWB/LZandCvNadADoM4FHOhWdI3vBKSXYvKw/jNf6l/eVS34HjGtLWoPWci1y1TxcTj3riXpi7P6B/DpPvyB3gv2t0WekNf7qbTG2SH33fdjQdTrO740WX2R6ZVcUh6DJrigOQZNdURyCJruiOARNdkVxCDnaXTbIOwWPFXUv39y9hUs0ASC5I3cPPZiXG0cAwKMVdpB2MrKIGPvSPU1IK/LjPjE2tjc/xsFu/P7YdLDslHr5852Osw/JZZp+gdxRtNjHsmu+9it23mPqyo+bdzOfx2ITZTf9XBOey5bQkrulDthYUdy+3MD9pM3b+YgYa4XLTPV+58TY5HA+D+mh8l0c303caGL2IS5rTd0u71eqMGMwv1eSGJsLwvPjLTvs4d9zE4/A3w6IsQGN+eQEHJPvmh1dU9Tt50tJctddQK/siuIYNNkVxSFosiuKQ9BkVxSHkKMGXXx8Hvww331E0IKVQ8XYvkefJK3AUNkoKbwkjrR1z1QRY997gcs3D6SGi7G/ebF5GLiO14JH/CyP3DldP4y0MuPkdcm5trNpdraNfAzhS3gc0roBi8TYUw15VNND1eVuunoqho8AAA3/SURBVB5+KaTZs2wslfuY1+kDwJ4hvK7/cJMxYuz9vbqS1r7da2Js6Q+Ok3Y0XjZ2l8/k89D4j5dJSyvBxwoAH9daSNrUimxcAsDBKbzOPaO3XOIcOY3Xrk/ewr8LAGYnlSbtzH1yS7e8u91TOJe8RD7z/7L+L0VR/klosiuKQ9BkVxSHcD0NJ32NMRuNMdtc458GuvTixpgNxpgDxpiZxpis7+YrinLbuZ6GkwZAHmttkqul9BoAbwDoBWCutXaGMWYsgG3WWtmNceFbspAtNsTdLElNliuOynbbS1rygxXE2GNN+D2r/2PzxNhZdcqTFrxEXmPevxA3H2z3Hq99T4yUt5/YiUccffhEGzE2tirP2Mh9VnZbfKPOk3bweV7HDQA7O35JWpUxshE2uMNk0nrP6ECaf1V59FLiDt6HivXkKrEHQ7ja7gl/uUFxl05vkHakufy6CSrBZu2a6lNJe+QdfkwAyBPNFXQ2l/z8JudnMy54uVydiMucZ2ll5cFK94/eRNrq3vcJkUCFT9yrR+e0W4ozu2NubD27zeSKDe7l+mMB1Acwx6VPQeb8N0VR7lCud0iEhzFmK4AzAJYDOAgg3lp75dITBZ3/pih3NNeV7K7JL1UBFEbm5Jdy1/sLrh7/lHFebiulKEr28z+58dbaeAArAdwHIK8x5sod/cIATmSxzX/GP3kEyqONFUXJfq7HjQ81xuR1/dsPQEMAe5CZ9C1dYR0ALMiunVQU5ea5Hje+MjINOA9kvjnMstZ+aIwpAWAGgGAAWwC0tdbyouyrCCoXbuuOa+WmeWdR35cizKk/3I3LCAEg+C9eV5x3k9x9dN8n7HobI58Djz38SaTYp9xp9ERPcaYlfB/i9dlPFd0mxo5f9yBvHyyXdLYry/swb2R9MTa9BbvTCUfyirG5i/AorvsKHSHt563yXZFS0/m59N4rd3E9PpZ7E1QOk5+zKoFcLruw/6NirP+P3Nvg8GQud710Sv6U+dyDa0mbubOGGDuk9vekjXi3tRj7/qDJpH3cp6MYO3HECNIa/9BLjH3hAffy769arcaJXfE3PP5pOzJnsl+rH0IWk1sVRbnz0Ao6RXEImuyK4hA02RXFIeTo+Keg3AVtnXIvuWmJg2UTys+LyxYPHJHXnR9oPI600vO6ibGhG/j9LXgbm1gA4PE5r9s+OzGStJB18nr2vd2F9eyT5bXgNhfvl7ksN5yMGsD+S8BMuQGjNMJqwrsjxViJeee54eSSUWwmAkCtbjyWqkfYCjH2yQm9SZPKiwHg5VFc3rv4TXnk1tvHnyAt5k1uHNp/+mRx+7YreY56yHq5NNdDeOn6dxDvQKNUIM+Y98jCGC7ow2O4Jq14WIz9toX7PPgXmp/AX9tTdfyTojgZTXZFcQia7IriEDTZFcUhaLIrikPIUTfep2gRG9HHvWlAwEG5G2e+fezGRzWQY8N4rT+s3G8AiUX5/a1i87+ESGBisR9JqzGBO7MWXCVXCR/tzG66927uTgsA/lH8PPi2OS3GRp3mkt8Sk8VQxJfk7rCJcmNWhOzkfej78TekNcnNZbUA8MTe5qTZ+rI7HdOZmzH4n8yidDqYn/eYx+W7OKV7cnnuia+4qYbPPLlk+Fx1PgcFuIIWANBxAHeHXdCQR00BQIMf95CWcll2+dfUL0za5WIFxNhcie7nYd2RKUhIOaVuvKI4GU12RXEImuyK4hA02RXFIeTo+CeJgOPyWJvca3lmuk+Ve8TYvNu4FPFA/9xibOFJbIocii4rxj79B5e7Rh7l9dJxM+Uy3qAFoawdlM28tAB+Kpbc850YW/s3Xtuc4S2PlUoOZa9mY6fhYmzTnmw+9lrYnrQvqsqmW9x0NpbSespuoG8MG2HdR80SYyc9+zhpHpcCxNhzE4NIW1dlCmkVz/P4KQDIs9OXtHwb5eMdN7wFadW+3y7Gzn+X198fbyyGAh+zVP4d7sYLAHs+de/xkPKpbGIDemVXFMegya4oDkGTXVEcws2Mf5psjDlsjNnq+lM1+3dXUZQb5XoMulQA9a8e/2SMWer6v97W2jn/ZVtFUe4Q/qdyWWNMbmTOeuvm+vPD/5LsYRVCbMupTdy0LqG/ibF9H+cZY+l55VJTz3NJpHmPk0s664bw7LHvj8kljiGd+XG/XM+HG+whO6B70/i9tP038owxv6qxpMWfkhtSmDR22Mv2leekRXWvQpr/CbkpxrphY0X9Wjanys5/t4F8bCFbuBEDAMRU43LVBLl5MEpO4rLhMjO54ywA7G/Jg4msl3BN85SfM+vFujktz7ZLLxlBWkoYlycDQNtBP5A2aF1TMTZ0Fd8xyr+IZx8CwNNrdrv9/MnTW3FkZ+KNl8teO/7JWrvhymMbY7YbYz4zxshHqSjKHcENjX8yxlQE8A4yx0DVQmbv+D7StlePf0qOkxcvKIqS/dzo+KfG1tpTrgmvqQAmIYse8lePf/LLxwULiqLkDDc6/ukvY0yESzPIHNcsf2lUFOWO4GbGP/0CIBSAAbAVQNer5riLBOYKsXV83A26e9byunUA+LTABtJqDZfNrYjHj5FWMI/cxfXUa5Gk5RoiGzCX3+Z10HYwd6L1eFlel3yyKRs4P7wld0V9YAUf2+oGo8TYOYkVSVv+qDxY92JlLmFtOHS1GLv4JJcjB77lTdrerlySCgC+BXhKb8pJeczSiw//StrX2+4XY3vXWEbazF6yueXdm0dIxaWwsTumwrfi9gfTuMR5ZH95pFNcGb5W+tZgoxUALqbweWxXnsd4AcD282wyHh1XRoyNqeT+84mRnyH1+PFbPv5JHi6mKModiVbQKYpD0GRXFIegya4oDkGTXVEcQo52l424J5/tMN3d1/v5pNw44mw0O77l32HXHQCOvlCKtNRg+bhKzeQy2hPvyuWjFw9zuWrQPjY6g57mjqYA4PsEz4CL7iSX5oZ/w00PgpbJRYlRI7iu1K+7vA/J/y5IWsP+shvfJHAbafPja5C2ObaouD0aRJHkGSF3RY2pH0naxQLytSdDOA3FFsh3UPZ14s67ZYcd5t9Vhee/AUBiEaFUtY38utt7lI/NXJQ9b5+wi6T5/iY34AjdzHc1TBZ5GlXf3+3nI1+PQPIp2Y3XK7uiOARNdkVxCJrsiuIQNNkVxSHkaHfZ2CR/zPjdfexPuY/YPAGAt1ZzieQ7feWyRc8kNi/KfLJbiASMH5dOLq0hL8nv0oFLMhfsWUna/Vvk/eq0ldfOfzxPNujSu/K6830bZeNw6QjuDpuYxSihlo1fIa2sL5eUAsCwE9zuNKEeG2EXfuTSTwAILMtG6TMLVomxn84tTlpahNx51+c4/75co8+LsUuKTyStRWJv0ooOlGc6xU5mQ9J/KJucABAWyumTkk+eO/Zmg6WkBVaVV4EOHtiWtHQf+XGvNSpPxssjtAC9siuKY9BkVxSHoMmuKA5Bk11RHIImu6I4hBx1431iL6PkTPfOpGnT5ZLQMa8/S5pfFfm9qchSbhiQXEde7H+sLbuVf6byTDcA2DugAmlNnmfHOaWm3G5r6uBmpPUat1CMHb6oOWk2WO7i2qMK3yWwGbJzHzydm3iMGPScGHuuOj9G2TzsGCculUtg4z7l3zXo+6fF2Nyn2V0OrX1GjPUazK+Ro4ns5gNA6Ov8uP7H+W7N2KNrxO0bLKxNmu85LnUFgOjaXO6aGirPLlydwGXhJ5rKr/3AcnzOM3zkbrj73nW/u5TSL+vrt17ZFcUhaLIrikPQZFcUh6DJrigOIUfXsxtjzgI46voxP4BzOfbLcw49rruPf9KxFbPWcotc5HCyu/1iY/6w1ta8Lb88G9Hjuvv4Jx/b1ejHeEVxCJrsiuIQbmeyj7uNvzs70eO6+/gnH9t/uG3f2RVFyVn0Y7yiOIQcT3ZjTGNjzF5jzAFjTN+c/v23EmPMRGPMGWPMzqu0YGPMcmPMftff3Nv4DscYU8QYs9IYs9sYs8sY84ZLv6uPzRjja4zZaIzZ5jqugS69uDFmg+s1OdMYI7fiucvJ0WQ3xngAGA2gCYAKAFobY3i1yd3DZADX9nLqC2CFtbY0gBWun+820gG8aa2tAKAOgFdcz9PdfmypAOpba6sAqAqgsTGmDoB/AfjMWlsKQByAzrdxH7ONnL6y3wvggLX2kLX2EoAZAFrk8D7cMqy1qwBcu+SuBTJHXMP195M5ulO3AGvtKWvtn65/JwLYA6AQ7vJjs5lcGSvu5fpjAdQHcKUR4V13XNdLTid7IQDHr/o5yqX9kwi31l7p6HgaQPjt3JmbxRgTicyR3RvwDzg2Y4yHMWYrgDMAlgM4CCDeWntl7fM/8TUJQA26bMVm3uq4a293GGP8AXwPoIe11q2d6916bNbaDGttVQCFkflJs9xt3qUcI6eT/QSAq4dsFXZp/ySijTERAOD6W+7IcIdjjPFCZqJ/a62d65L/EccGANbaeAArAdwHIK8x5kojl3/iaxJAzif7JgClXe6nN4DnAMitW+5eFgLo4Pp3BwALbuO+3BDGGAPgawB7rLUjrvqvu/rYjDGhxpi8rn/7AWiITD9iJYCWrrC77riulxwvqjHGNAUwEoAHgInW2k9ydAduIcaY7wA8jMxVU9EAPgAwH8AsAEWRucLvWWst9826gzHG1AOwGsAOAFd6VfVD5vf2u/bYjDGVkWnAeSDzQjfLWvuhMaYEMs3iYABbALS11soTK+5itIJOURyCGnSK4hA02RXFIWiyK4pD0GRXFIegya4oDkGTXVEcgia7ojgETXZFcQj/Dw5auMbL7pB6AAAAAElFTkSuQmCC\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"x = np.random.random((40, 40)) * .2\n",
"\n",
"plt.imshow(x)"
]
},
{
"cell_type": "code",
"execution_count": 48,
"metadata": {},
"outputs": [],
"source": [
"for i in range(40):\n",
" x[i, :] += np.arange(0, 40) * .04 + 5 "
]
},
{
"cell_type": "code",
"execution_count": 49,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.colorbar.Colorbar at 0x7faf81d94c90>"
]
},
"execution_count": 49,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAS4AAAD6CAYAAAAfmKrOAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAgAElEQVR4nO2de4xkd3Xnv6duPfsx3dPd0z09Y4/Hb8dAsMGPkLAYjAjG2ZVBQigkSyALchIlEdFmV2YjJSTZjUQECWSVhewECI6WhCBjAkJAsFgngAjGNjF+Gxt77HnPeGb6VdX1unX2j6peuup8f123Z6q7q9rnI1111+nf/d3fffSp3z3nd84RVYXjOM4gkdrqATiO46wXV1yO4wwcrrgcxxk4XHE5jjNwuOJyHGfgcMXlOM7AcV6KS0RuEZGnROQZEflArwblOI6zFnKu67hEJALwIwBvAnAYwP0A3qmqj4f2yUYFLWTG2oXrOL5m0nwsDduHivC25HhxPko8BjZejfixlHwtMBkApOrJ21I5H8K6xiAN0pZdmtT5r/1Lpe3BohQZAIC4YQeciWLer5D7S054KKom3j9iFybQb1bIjQQQkznCkFRoW/bohv5NctJ+HY4cjnHmTCPwNCTjzW8Y1tNn+PXt5MGHK/+kqrecz/HOBa4JknEDgGdU9VkAEJHPAbgNQFBxFTJj+NkL3tUujJNdIACo7Zmg8mjRPgBayNC2UrUP1uJlO2hbppCiqn2Il3dyxVcfsvvHedoUhVP2yawXeFvWbz3QLztefYT/I6aLtt/qTttWC8nvGWL+PzQ0VTKy8eFl2na+ZC/E7NgCbTuaKdv9q3b/a3YepvsXopqRjaXtWAHgbG3YyPbnX6Rt52M7hlcVDtK2GbHXt0a/QYBL0vNtn9/6C/z46+H0mRjf/6d9idpGs09PnfcBz4HzUVx7ARxa9fkwgBvPbziO42w1CqAB/uXWL5yP4kqEiNwO4HYAyKdHN/pwjuOcJwpFTdcxq94CzkdxHQFw4arPF7RkbajqAQAHAGAsv9sDIx1nANjOM677AVwuIhejqbB+EcAvrbmHwlgZdWSINm3k7dBSdX4xpWbtVpXd1v7Q7CRnRPkz1q4BAI2MNagu7bG2swY3pyGqWj3dSIes6FZUz/O2KTbcgD0MxNicXuLW+cq0vY5CbFRS4fYWHbL7p4kMACoVe9EWAkZ/Zs+qxXwMJ2t2Vp9L2zEsxfY5AIAKuZknKtwG+sZxa849Vedtry88Z2RPVmZp21fkDxlZlti9AOD5juNV9Sxttx4UirjPky+cs+JS1bqI/BaAfwIQAfi0qj7Ws5E5jrNlNNg3aR9xXjYuVf0qgK/2aCyO4/QBCiDezorLcZztybaecTmOs/1QALXtauNyHGd7olB/VWxDG0C5fZW7nuQrfVN7ZoysPs09NlLIGlnuFF/tXNtp3W9RmXu+6kPEAzlnvTvFae7hSi+T8JMs9xRWx6w8FVhKw1bU10b4g5Zetv3WhnlbqZGxrSealXgbpcA9tnHVPnrFRb78f7lgV8NPDy3StocWdhoZC+M5VR6h++cj+yy8fPQobfvo8gVGNpXh4zpeHzOy3Zl50hLYE9lndz7guu5cZc/Odd0oEPe33vIZl+M47TRXzvc3rrgcx+lAEIei9vsEV1yO47TRNM73t+LyRIKO47TRXMclibYkiMi4iNwlIk+KyBMi8hrS5vUi8pCIPCYi/9Ktz02dcWk9RvzimTZZtHuaN14sGlFqjIcH1XZao260zA3umTPW8Fm8hBv90yX7ph/nrK4fPsmt6MUZa6wOZCeh6WcCUSlgWVbyp/lDVN5FHATD3IKhGZKLapEY3AOGW91rjegNkksLAJQZ8nP8Oi6W7YUYz/MUOJfutM6ew4vjRnbB0Bzdf65mPR95GmPFU81EAW/c81Wb/eXGoR/TtidICpwrMzx319GOkKxUj7yBjd7OuP4CwNdV9e0ikgXQ9o8sIuMAPg7gFlV9QUQCSuEn+Kui4zhtrMy4eoGIjAF4HYD3AICqVgF0ZnH8JQB3q+oLrTYnu/Xrr4qO47ShEMRIJdoScDGAUwD+RkT+TUQ+KSKdGRCuALBTRP5ZRB4UkV/p1qkrLsdxDA2VRBuAKRF5YNV2e0dXaQCvAvAJVb0WQBFAZ32KNIBXA/gFAG8G8PsicsVa4/NXRcdx2lAIqiFjrOVFVb1ujb8fBnBYVe9rfb4LVnEdBnBaVYsAiiLyLQCvRLOmBcVnXI7jtNFcgJpKtHXtS/U4gEMicmVL9EbYuhRfAvBaEUmLyBCaKeCfWKvfzZ9xpRIa/SLikUvzC5U5bd1s5Vke0pHO2H6z89wDWRtJdnnqeT6u7IL18CxP8fOPiJMsxR1JFJI/rzmGs/Z4jUxgDAv2PKpT9tpEC4HrwirslHjbzIitsjM8xE94ccl62U5m+P091rAXYu8Om4hwJOLHYoUqhgI3YpG4gllRDACYzVgv5sEarzNx89CzRnZ/xYYyAcArsu2JA9M98ir2eAHqbwP4bMuj+CyAXxWRXwcAVf0rVX1CRL4O4GE0F+1/UlUfXatDf1V0HKcNVaGl1869P30IQOfr5F91tPkwgA8n7dMVl+M4hoaH/DiOM0g0jfP9rRrOa3QichDAIoAYQL2Ld8FxnAFgxTjfz/RCrb5BVROVz5V0hGiiw8iY4UOId9oqPVILJKgiESz5Y0u8ad7mNdKY36RUzhry4xyp9jyavHJPaAZOiiIjy1M7QerJO14mFeEyi/x8y7utIT5z1t6fRpobgPWkNVYLCSMCgBqpEl9O81CkVMrKR3MBg3nFhgctVq3syUWb7w0AxrI2bGkxx/OElRo2D1woVIYZ/d88wu3Piw373F2d5dV7TsTtz3OtR694cZ8HWff3fNBxnE1nZeV8P3O+o1MA32gt0+9cMes4zoDS0FSibas43xnXa1X1SCua+x4ReVJVv7W6QUuh3Q4A+YivvXEcp39oBllv4xmXqh5p/TwJ4IsAbiBtDqjqdap6XTYVKrfsOE6/oBDUNEq0bRXnrLhEZFhERld+B/DzANZc7eo4Tv+jCsSaSrRtFefzqjgD4IsistLP36nq19fcQwRIdxyyxsNtopO2AormrRcHAOJJ+woaLa4jXiZAI209K8Q5REN7AKAybvcP3euhE7aP4h7u2dHIytM272IT0kV9KFDlp07Cg4gHMR5ZRymFQNbBwoSNccpnecK+hapNIFmL+bf99LD1JjNPX7XBH/1KbOUPzl9E29488aSRRYEyEy9UJ40sVJHnkcoeI3vLMK80tCtqf6Dy0gtvoGzfBaiq+iyaEdyO42wjFNjS2VQSfDmE4ziGfjfOu+JyHKcNhfQ653zPccXlOE4bzfJk/a0aNnd0IjbEJ2CcR80aaqXBDZ9Cqv8wA3azD5Ija5qX00kvWUu8kPCg0jS/jCzlU3aBn8PiPttvoJo76iQ8qDrODb1R2V6HWsi4TsS08nsUMO4v23NIT9sQGgBIpYjRP1ARaM+UzWVVrPKS9PtGbWjMjowdQzbFn7s6Cbe5qMAj2p4q7TayVwwfpm2vyB83skUSMhTixZiHux3pyP9V7MlMyQvCOo4zYCiwpavik9Dfo3McZ0vY7IKwrXbXi0hdRN7erU+fcTmO04aq9HrGtWZBWAAQkQjAnwL4RpIOXXE5jtNG0zjfm3CehAVhgWZe+i8AuD5Jv/6q6DhOB9LLkJ+uBWFFZC+AtwH4RNIRbvqMSztCFJDmHr3G9JiRRae5m622w/aRPc29MA1S5SdUGKW023quCqesN0oCTrqoYjuujHG7QNoWKgrCPH0h73WK9JsioT0AkC6SKj8T9jqmioFvY3Id83kexlOpkGs7zNuypHaNgAfy8RdtgsDdo/a5YaE9AHDJ6GnblrpWgcmMjbPKCz+HR5cvMLK52HrDAeCWYVuZayhgTtqfbg9xyiGQbHMdNI3zib2KUyLywKrPB1T1wKrPKwVhf1tV7xORv0CzruLvr2rzMQB3qGpDEoYs+aui4ziGdayc70VB2OsAfK6ltKYA3CoidVX9x1Cnrrgcx2mjlyvnVfW4iBwSkStV9SmQgrCqevHK7yLyGQBfWUtpAa64HMch9LhYxpoFYc+lQ1dcjuO0oQrUAjbEc+uve0HYVW3fk6TPzVVcDYWU2sMvdIhXUEktklCRMs+xlapag2R1ihs+U3VrSWd5twAgXbbW5kbW3tChk9wgWpyxRuxUyHZKcjPVAxbZHCn4UgollyUG80ag8g4L+dFhO2Ct8oc6GrWG6cVjo/xYJM+XjNgcXQDPvTWUYx51oFq3j/REznoojpV20P0vJuE9SzF3II2l7XhDBndW5efnCs/QtnMkFOhQg/+f7IraHQS9CNVpvir294IDn3E5jmPwWEXHcQaKdS6H2BJccTmO04G/KjqOM4AMfM55Efk0gH8P4KSqvrwlmwDwDwD2AzgI4B2qymuEt3cGZNtXIUs9YK1mK2ijwIptknsrs8CNt/GQPeV0iS99V2K0T5fseBf28bxKOVJEI7RyPs5beZ3bY5Ev2X7zpwL9Erty9iz/Nl2+wEYFpBbs9WqM8FxWcdG2HZ7hVTxKS9zgzRjNWadMqcZXszOeW5gwsskCD1U4W7fG9b3MGwJgPLJ9hIzzQyQ524NlXoTjZbkjRnYdcTAAwA+q7V6Zeg+WMTS9iltXeiwJSc7yMwBu6ZB9AMA3VfVyAN+EXQnrOM6AsrIANcm2VXRVXK3K1Gc6xLcBuLP1+50A3trjcTmOs4U0WiXKum1bxbnauGZU9Vjr9+No1likiMjtAG4HgHw6sKbHcZy+4SXhVVRVFQlUtmz+/QCAAwAwltsdbOc4Tv+wXb2KJ0RkVlWPicgsgJO9HJTjOFuHqqC+TRXXlwG8G8CHWj+/lGgvIfm45m3JdABo7Js2slSJe5Iyp6znKt7BXXJKvJVxPhDCUrXeRtY2lAsrqtoJZnaBt41zdlyFF/kEtTZs25YnedvMIsllleVtcyeIB5G0beQCDzXptrTA70OmYMODyjV+IStEPruDX8iRjPXejZEqP4WI582aytjnMZQNlHkQD5WtBxMArh95zsiuzh0jLYFTsS3jNNfg4VCvyraf21AoOdw6GfhXRRH5ewCvRzNh2GEAH0RTYX1eRN4L4HkA79jIQTqOs3lsCxuXqr4z8Kc39ngsjuP0CQOvuBzHeWnRy0SCG4UrLsdxDAMf8tNT6jFwuqOc+h5rhAeA6LgNs6jt28Xblmx4T7TAS79Xx20RjtwcN9QWZ20oT3rZWqAjfigs7SH5uGoBw/iclZd284enOkpyWQVssuVp+wcNGufteOsjdv+oxI3zjUl7HTVQmGMob+8ZM8IDQL1mx3UoHqdtX73nkB0X+Sc8uszzce3K2sIaTxf5M/q6nT8ysssLJ2jbYsOGOC2SvFsAsCeyYxgWfs2fqrXLyz2YKakC9R4mEhSRcQCfBPByNE1o/0lV/3XV338ZwB0ABMAigN9Q1R+u1afPuBzHMfT4VbFbQdjnANykqmdF5C1orvu8ca0OXXE5jtNGL21cSQrCqup3V338HgBby62D/l5l5jjOlqAqibYEdC0I28F7AXytW6euuBzHMawjyHpKRB5Ytd3e0dVKQdhPqOq1AIoIZJMRkTegqbju6DY+f1V0HKcN1XXZuHpREBYi8tNoGvDfoqq2nHgHm6q4NI4Rz7V7FdPj3LvDSAc8hVK04RDVC3joRXrZJsGrjfDLwDyALDSH5J5rjiu2+6e4AxOVcdtvyFuZIkkHGznuKUyXSL+n+US7toN4K2skRGonTySYOmtDsliVIABYWLRliSZ38vCv7IjtI+T1+uGJvUZ2w+wLRrZY46FIjy3OGtl0jo9rPg6VVrIMpawX9f7lS2jb6wvPGtmJmD+jnQkG8+F8B+tAEPfIq5ikIKyI7ANwN4B3qap11RJ8xuU4jiGh/Sop3QrC/gGASQAfl2Yscb3LLM4Vl+M47fQ6VrFbQVhVfR+A962nT1dcjuO0o007Vz/jistxHIOH/KxGBJJtD3PQbKBay6I1iMajU7Qpy5akaW5crI4lP+XcWWuEjgv2aFGV32SWi60yyscVk+iPVKAAEkkZBVngYyiRpLMNEjIEAPEwiRsiY2CVfwBAM7bf0Wlu2M5EtuPlKn8WlmHl06O8X9ZHhVSsmczz6kM/NXLcyEYDXpKD5Ukjywi/aTeO/tjISiQMCAAO1mxo28uyR2nbQ/X2e0ZSwK0b7aFxfqPwGZfjOAZ/VXQcZ+DosVex57jichynDVVXXI7jDCD9nkiwqwVORD4tIidF5NFVsj8UkSMi8lBru3Vjh+k4zmaimmzbKpLMuD4D4C8B/G2H/KOq+pH1HEyyGUR7O0IqyrYqS3Nkdmjp53kVtNrFth5tHKhEE1Xs1U4vBeJwIhLeM0KS7ZEQHACoDVl5bZS3jchlqAdi6JnjKsWjcMAK1GQDHshU3V6z6qQ9mCwHvKjDdhCLJ0f4wCJ7H3ZMck9fmngg04HMiUM5G1rz5Bn7fFw2/iLfn9yIp0q83vFPDdsqPSerPITtmfJuI9uTtckyAWA0ZUPYJgJVifal269vIZBwcD0oBI0+9yp2HZ2qfgvAmU0Yi+M4fYIm3LaK81GrvyUiD7deJXf2bESO42wt2tN8XBvCuSquTwC4FMA1AI4B+LNQQxG5fSVXTzUuhZo5jtNP9PmU65wUl6qeUNVYVRsA/hrADWu0PaCq16nqddkokP/FcZy+ot9nXOe0HEJEZlV1xTL5NgCPrtX+/1Oro3G83cCemuGVe+K9NrynkechIZq2FzBd5qEXDRIKpBmuv+skvKeyw8ryZ/mxasOkyk/AFyHERZOdIw0BlGaZ0Z9//RVOJM8fVt1JjOBzxBkxYw3gAICKbTs8zQ3u1Sova88YL9iQGwnknbp0zOage2bOPksnlkfp/o+kbD6v2fwCbcsM8a8bfYq2PR1bJ8WZOndc3ETycT1S5eFu+9KBpG3ngQJoNPp7OURXxSUifw/g9WimaD0M4IMAXi8i16B5jgcB/NoGjtFxnM1EAfT5Oq6uiktV30nEn9qAsTiO0yf0e6xify/WcBxna+ihcV5ExkXkLhF5UkSeEJHXdPxdROR/isgzrZUKr+rWp4f8OI7TQc8N790Kwr4FwOWt7UY0Vy2sWRDWZ1yO41h6NONaVRD2U0CzIKyqdrqdbgPwt9rkewDGRcRWLVnF5s64UgIptFdX0QJPphYdtd4hmebrXKVqQ02qu7nXKF20oRP1Ie6tbGTItw65WZUxrv8zSyS8KMXvdm2YVO4h3lIAyBInl6Z42+pOe7z6MB+D1Ox5sF5lKfDYjNr7UDzJ45ZyEzaspVjiz0Iha+9ZKNFdOmVDga7bdcjI5mq8Qk85DiS2JCzFdrxzMXfZHq2NG9mFGR6QcqZhs0rmhYf81LTdE6y9WFylgPbOq7i6IOwrATwI4P2qutrdvBfA6pt0uCWzMVUtfMblOA5BEm69Kwi7HtzG5TiOJfnErRcFYY8AuHDV5wtasiA+43Icx9IjG5eqHgdwSESubIlMQVgAXwbwKy3v4s8AmF+1wJ3iMy7Hcdrp/QLUbgVhvwrgVgDPACgB+NVuHW6u4ooiYKLdSCln5mlTJSvgpMTjZeq7rCE+ZKxOLVsjZ203L8eeO2vbMoN5NVC5p7bDtmX5wACgStrmz/C2pT1EGHjO0kX7B5ajCwDqpMpPbZzkvQqE20RZGzIUx3xglTl7zQvEYA8ACyXb9qrpE7Tt8/PWgcNyd106coruf6Jiw3iWAwb7AsmR9XDpQtISyJGEaYcwQduORvY61JT/q5a0vdpRo0eRz71cgJqgIKwC+M319OkzLsdxLIMeq+g4zkuPwKS6b3DF5ThOO1ud3jQBrrgcx+lABj87hOM4L0F8xrUKBaTeEaJQDiRC2zVpZdlAOIbYb4dUnVeB0aw95cxCKBGgbVsbsh7EVJ3fZVZ5p04q/wBARHLzlafWUxGIjyHOW3kUqtJDlvWxlhpY/dcYsiec38E9wSkSmkNuIwAgIm2PFXk1nR15e7xi3YbQsNAgALh82FaSOhGo3POywnNG9kjAq7g3wyv6MF6VtRWInqoFqgfV2l3E5V7NlPjl6Rt8xuU4TjvbIZGg4zgvPdyr6DjO4NHniqtrrKKIXCgi94rI4yLymIi8vyWfEJF7ROTp1k+vreg4zqaQZMZVB/C7qvoDERkF8KCI3APgPQC+qaofEpEPoBnxfcdaHWm1ivrzh9tk0YTNUwQAOEtCgXbxEIlUmeTY2sFzOy3P2DxMNO8WgHqehOxU7VdRyDifKVr5/MU83iZk8GYI8SWkl/g5xAU2Nj7e+rjtOHvajjfOBkJ+Mnb/avn8J/X7dtu8VdkUd6jszNvanafLNifYfaf20/1vmnnayPbkeLmlBvnez7CbA2AyWjKyfIrn2DrTsNf82hyvlnRfud1oH6/nQVqDfn9V7HqWqnpMVX/Q+n0RwBNoJvm6DcCdrWZ3AnjrRg3ScZxNRNEM+UmybRHr+joUkf0ArgVwH4CZVaknjgOY6enIHMfZOvp8xpVYcYnICIAvAPgdVV2QVYtuVFUlUKGzlRHxdgDImxz5juP0IwP/qggAIpJBU2l9VlXvbolPrCS0b/20K/cAqOoBVb1OVa/LgNudHMfpM3pYnmwjSFLJWtCs0PGEqv75qj99GcC7AXyo9fNLXY8mAsl0HDJgnNe8XSWvOb5yPh4iq+FPWSMtAMRZUvZcuP4unLaG1tqINZymavwOLu2xbYdO8iXJ1RG2+p82xeJFVhbKsZUuJSv4AQCZOdsJ65ektwIA1KvkcQrYQdI5e3KNQAGMYtWufD9ZtjIAAHmcUuSEp4cW6e6nqrzICoP1++rhg4n3P1Ql0SEA9qTtKvv7yvxt5cJ0uxMr5BxYN30+40ryqvhzAN4F4BEReagl+z00FdbnReS9AJ4H8I6NGaLjOJuJaG9fFUXkIIBFADGAemeO+lYJs/8DYB+aOukjqvo3a/XZVXGp6ncQzK+JN3YftuM4A0fvPYZvUFUbhNnkNwE8rqr/QUR2AXhKRD6rqiSCt4mvnHccx7DJxnkFMNoyS40AOIPm+tEgXuXHcRxLb43zCuAbIvIgqbsIAH8J4KcAHAXwCJoFY9fMT+EzLsdx2lmfjWtKRB5Y9fmAqh7oaPNaVT0iItMA7hGRJ1X1W6v+/mYADwG4GcClrTbfVlVSs73JpiouSUeIJjpCGqs87EGIPJ7kHp9oyb4KV2d46fdGhuScagRyWeVs26hivwhCIUPsG6kyxie5zDO5PB3Ix0VSmGkgVRkL+amN8S8zqdrjNXbYGXs0xx+bNAn5YTIAWD5rQ6+EVAkCgOWqPblLp07TtqeXrfdtJGOfj+mcDcEBgEbQnGs5XLHhuRNp3u98TMbFbiSAxQavOsXYFbXf39CjuG56VxAWqnqk9fOkiHwRwA0AViuuXwXwoVa1n2dE5DkAVwH4fqhPf1V0HMcgjWRb135EhlsxzhCRYQA/D+DRjmYvoOXoE5EZAFeiWX8xiL8qOo6zkcwA+GIr0iYN4O9U9esdBWH/O4DPiMgjaK5guGMNDyRWOnIcx2mnR15FVX0WwCuJfHVB2KNozsQS44rLcZx2erwAdSPYXMUVN9BYaA+1kL27aVMW8oMUN8lJjRh1Axe+cMzmNSrP8HAKiW0nTBYyFKZJRfmY7A8ApJo7MjwqBQ1yaWqBO1knxvn8cR4fVJ4l17HGnBmB3F8nrMG9VuAG92iUO2UYO4fshUwHQlvGctbgzXJ3PTbHn7tdBWtcv2jI5gMDgNmszRm3K81vGsu9tTtNcs4BOF63cUu70zwn2NG4/V6SdHHnhisux3EGDldcjuMMEoJkHsOtxBWX4zjtuI3LcZyBxBWX4zgDhyuuVaQEUugIZ1B+hTQKZMZj3S6QpIGT1sMFAHHBuuRqw9wvmDtrw12WLrAJ7HILgeSAY9b7xrySAJ+ah+wMrG0jUHmH9VHexTvOnLXXoTpjr0G825a5BwA5Y6+N5ALhRUyW4m0jIm8EqtmMZOzYUuSCFes8EeEVIzaR74+WpmnbS/O2bTkQe5UiNe0XAqE9EWkbBTTJyzLt51GQ3sT8+Kui4ziDhysux3EGCnWvouM4g0ifz7i6ZocQkQtF5F4ReVxEHhOR97fkfygiR0TkodZ268YP13GczWAl73y3batIMuOqA/hdVf1BKz3FgyJyT+tvH1XVjyQ+WhRBdnTk1Ark40oVrZFRC9yg2hi1hvjcUZ6DrHiZzaGUP8OzxFbGSfWgkr1boarn2XnSNuBzWJ6y55td5E9GfZi0nQ+E4eRIjq0M77cyxUJ+7P7RWV5mLp62ea+0wi/O+JQNjanU+cWZK3FHC2P/qA3PoVV+xnhozsPze43s300+Tdseq9nQnBj8ubsm/7yR3b98CW37s0P2eFcH8po9VG2/vuTxPDf6fMaVpFjGMQDHWr8visgTAOzddRxne7DFNROTsK5EgiKyH8C1AO5riX5LRB4WkU+LiJ3KOI4zcAj6/1UxseISkRE0q1n/TisX9CfQzA99DZozsj8L7He7iDwgIg9UY16k1XGc/qLfFVcir6KIZNBUWp9V1bsBQFVPrPr7XwP4Ctu3lTj/AACM5Xf3+QTUcRwAPX1V7FYQttXm9QA+BiCDZh77m9bqs6viatU6+xSAJ1T1z1fJZ1v2LwB4G2weacdxBpXeTzGCBWFFZBzAxwHcoqovtKoBrUmSGdfPAXgXgEdE5KGW7PcAvFNErkHzFA8C+LWuPcUN6GJ7Ij8t80onetk+u/twIJyiblfLSSDpYO5FHq5CxzBpj9cgV0wDYRaZZRLmsS/gViTDre7g/ZKcdKgHCsPUh+0TmF7k/aYX7djqE9bj2pjl9wzz1usbTfDrXSUexHrAq1iuJzfFPh3vMrJ9o2eNLE2SCwLcg3j/3H7a9uJhW2nowgyvPvTd0uVGNpri1/Fo3ZqLh+UUbVvW9geyJ/pm818DfwnA3ar6AtCsBtRthyRexe+Ah5Z9dd3DcxxnMOit4lopCKsA/jepu3gFgIyI/DOAUQB/oap/u1aHvnLecVuprEkAABEPSURBVBzDOkJ+elEQNg3g1WiWKCsA+FcR+Z6q/ih0UFdcjuMY1vGq2IuCsIcBnFbVIoCiiHwLzcpAQcXlBWEdx2lH17F1IWFB2C8BeK2IpEVkCMCNAJ5Yq9/NnXEJIOkOA2yaDyE1Z6utRC9yg2o8Y0MvUkVu+CztmzSy3BwPO8rOk1xUeavrizPcqMzsv4UX+d1mRv+YR9agTMKDUjbaBgAQlWzb+khgDDkrz560A6tdwO+Dkoo+6TRvW4/tdRwbJmWRABQr1ujPjPsAkOl8vgDkInsfRyLuNHh8aY+RvXrsBdp2ilT0+fbilbTtq4cPGtloxM/3FVnrfDseeBiuzrQ/5/leWdV7Z+PqWhBWVZ8Qka8DeBhAA8AnVXXNVQr+qug4ThsrK+d7QZKCsK3PHwbw4aT9uuJyHMcgjf5eK+6Ky3GcdgYgyNoVl+M4Bs857zjO4OGKaxUiQLY9jEYXeEI35EnSwEBoTSNrPUkp4l0CgOGD9nilfaOkJZAuWo9Ydcx6w/Jn+Wq9Rpp49AI58VjCv6jMn56MdbiiNsL7Zd+cmUDIDx0XuQ1aDnhRl608FGCVytlruzPgVUyTKj+q/BwqNftIL8c2dOtIbD3RADCVLRpZah2rMVlFIQA4WiNZn5gMwCuJV7HU4F7FU412d3K9RwrHZ1yO4wwerrgcxxkovMqP4ziDRi/XcW0Urrgcx7EEKsz3C5uruOIGdL69CorsvyDx7jLHDfk6MWRkjWFuzCxP28RVzIgeIqqQ3F+BaXV11PabWeIPRJpktQ7m4yIW79oM71dJRR8hlXsAHvKjkZVFJG8XAMQkd9fohDV2A8DSgvVSHD89RtvOTM4bWbbBQ4mu2GlTOZWJcX44zWOkdmXtMxbKx3X9+EEjuzjH82YtNuxzd+sIj2opEcfDRMTTntc6Skz1St34jMtxnMHCF6A6jjOIuHHecZyBo98VV9d8XCKSF5Hvi8gPReQxEfmjlvxiEblPRJ4RkX8QEV5m2nGcwULRNM4n2baIJDOuCoCbVXWpVabsOyLyNQD/GcBHVfVzIvJXAN6LZq3FMCmB5DuMlEvc6KgLZHn4GF/hHhVtPq3qJK8eUThs+y3u5/2WdlujbnbBfhXVhrn+L+9MvsqeFbuQOGBwj6zxNrvADe7lKWJc53Zt1PN2bKmKPYd4LNBBbMdQLvMCJ6M77Cr50TxfZ//iwrCRTe7gRv8zFdt2//AZI3uhxFetF+v2+3c0zcd1srrDyI4rdzDsSNv8cN8uXUbb3jRkC3Ycr/Nn9Gfz7c6EbHI/05r0u3G+64xLm6z8t2damwK4GcBdLfmdAN66ISN0HGfz6VEG1I0iUepmEYlapclOArgHwI8BzKnqiv/7MIC9GzNEx3E2k5UFqL2qZC0iB0XkERF5qKOwRme760WkLiJv79ZnIuO8qsYArmkVbvwigKuSDRkQkdsB3A4A+SgQCew4Tv+guhGJBIMFYYHm5AjAnwL4RpLO1lUsQ1XnANwL4DUAxkVkRfFdAOBIYJ8Dqnqdql6XTQVSIziO019s/qvibwP4AppvdV1J4lXc1ZppQUQKAN6EZgWOewGsTOnejWalDsdxtgG9fFXETwrCPth6A2s/lsheAG9DN+feKpK8Ks4CuLM1lUsB+LyqfkVEHgfwORH5HwD+DcCnunclpqqPZpIvJdMcX3GRKtvwjfxRXrmnMmu9M6kavwO5Oes9YyE/9SGu/6OK7Xd5irdl62aYpxEA6tZxhogXNUJ62bqZNPB1xTyIDeJpZN5DAJC6ldcCJxERb+dywN25e9yG4cwv836XyjbUa4iE94wEPIWXDNu3mUYg99cVheNG9nxlira9aeRJI/txdZq2HSb5x3aTikIAkOqYewgtOr9OFEDyV8VeFIT9GIA7VLUhgZx7nXTVGqr6MIBrifxZNAs7Oo6z3djcgrDXoTkJAoApALeKSF1V/zHUp6+cdxzH0Kt1XK0isClVXVxVEPaPV7dR1YtXtf8MgK+spbQAV1yO4xB66FXsWhD2XDp1xeU4Tjs99BgmLQi7Sv6eJP1uquLSeg3x8RNtsmgfz8els7usLMWtytVpa62Olm1uKACISlYeqISO0h5rAGbhNullHsZTK9i8VYUXeduYxGqkeMooNDK2bWhqHxHjfHmaG8E1SzpJkXxe5YCDYcIOeGiInwQLBcqP8nt2umjzrV0+yfNeNYjn4XjRhubcuOsg3X+uRnK7BQzej5XsmuuXDdFVQSiqdSxdlTvK2zaIkyTgIMhJ+79wLyJ+mgtQ+zvmx2dcjuNY+jw7hCsux3EMPuNyHGew8AyojuMMHhsSq9hTXHE5jmPxV8WfIFGE1Hh7ojUt8kSCUiKuvgmepC13bMHINMcT2NVHrHenupOHEtVzxHtHHHIhe4CSYjh14pUMta2OJ/cRsXGFyCxyryBzXNWmbehUapJ7ChskFKh0jGcEkR22jxNnrPcPAHI5O4YzZRL3BB4KdOWUjdt9epGH28wU7LPEPJUAcNXIMSN7pMi95G8ee8TIMoGb9t3lS4xsd2aOtl2otMd6LQW8j+vCC8I6jjOQ+IzLcZyBo7/1lisux3Es0ujvd0VXXI7jtKPwBahdqfK8WXSmWg9YoMm3Q32SG2+jiu0js8hDTTKk0FCDGNeLe7hxP10mJe0DubCYIyC2qaUA8BClkE2WGf1rY/ypTC/YwUnRPiJBP0CdnNwQv7ZswNM7rWEcAE7P23sZyse1d2zeyKaytiJQMeL3rNqw51uJ+b/Jj8q7jexVI8/Tto+WLzSyy3I2nxcAXE7kp+rccfGKbHv+sFz47iRGoL4A1XGcAcQVl+M4A4crLsdxBortYOMSkTyaaVZzrfZ3qeoHW5kKbwKwYlR4j6o+tFEDdRxn89gOXsUKgJtVdUlEMgC+IyJfa/3tv6rqXWvs6zjOwKE9fVUUkYMAFtH069Q7c9SLyC8DuAPNVGCLAH5DVX+4Vp9JimUogBX/Wqa1ndtZpVKQQnttRd3BvX84YautSIWHmrBKQZrm7rvyDhIKFPDIZZasR2zpIuvqC3n/GmQMEfE0AkBthCQotM6wZr/kFEJVftidyp4mrkYAld3Ew5u137ypNP82bjCvYoUfi3kb54uBupvEAzkzyqveHF2w3rezZdvvUIZ7s4dJRaB8mrct1u2NZ55GAHgdqfJTVh6WNhfbZIYhZqL2MaRlXaVSOYqNsHGtVRD2OQA3qepZEXkLgAMAblyrs0RnKSKRiDyEZrHGe1T1vtaf/kREHhaRj4pI4N/XcZyBo5Fw6wGq+l1VPdv6+D00C0yvSSLFpaqxql7T6vAGEXk5gP8G4CoA1wOYQHOqZxCR20XkARF5oBoHciQ7jtNXiGqiDa26iqs2U/AVXQrCdvBeAF/r0mZ9XkVVnRORewHcoqofaYkrIvI3AP5LYJ8DaE79MJab6W8fq+M4TZK/Knatq4juBWEBACLyBjQV12u7HbTrjEtEdonIeOv3AoA3AXhSRGZbMgHwVgCPduvLcZwBQBWIG8m2RN39pCAsgJWCsG2IyE8D+CSA21T1dLc+k8y4ZgHcKSIRmoru86r6FRH5vyKyC03T9kMAfr3rCVRrqB863CaLrriUtm3snzWy1AJ/1axPkpxPAYN7/qTtY3mWG0M1ZTuJKvabKFS5pzJmvxdKM/y7IkXsvzUe5UHDhiJeUZ72W97Jxytla0hPn7WPSG2ch5VI3sol4scaHrXeBA3ELc2MWUP8MWKEB4CLdp41siFicL9o6Azdf4iUVjpc3knbTmZtTNi+LP+fOx3bZ3RfOjQGkn8sYLA/Wm+/8TXtmeGpJ90kKQgrIvsA3A3gXar6oyT9JvEqPgzgWiK/OckBHMcZQHrnVUxSEPYPAEwC+HirnVky0YmvnHccpx0F0KOc80kKwqrq+wC8bz39uuJyHKcDBXr1yrlBuOJyHKcdRWLD+1bhistxHItnh/gJks0gvacjoVqNJ5pLzdvqP7LEKwJFeRs6IcoTxdFxxYEwnGF7eTIl+020PMnDWnILtm3+DD8WK/hSH+JeNhbeU5mgTYPeRobm7HhreTLeKHAOZ+19aOT5NzerRpPN89Ca00XrUWs0uHe2WLP3/VTJhpUtVHkiwomcfcauGuEJ/1iVnlKDB5CMkuyPP67xSkMRWZLO9geAcoeLuRFyp68XV1yO4wwWvQ2y3ghccTmO046CpkPvJ1xxOY5j8RmX4ziDhbpXsY1GA1psN35qKZAx4or9VkbybgEAhITmzPN+NWcNyLVhbugdOm7DP+Yut0bdVKCQTXG37ZeF4DTHZWURTz+GyqT9NmwEfBEsPChUaUhq9jpm5q3jocrydgFQYrQv7OIOlVqVOD4ygVAiIpsd4xWBZgo2POhMxRr3S8SIDwCXjNiQnedKU7TtBKkeNMVKQwGYjKx8PMWf0UN162kZTfGEa6ca7edWT5bwZW0UUF/H5TjOwNGjlfMbhSsux3EsbuNyHGegUHWvouM4A4jPuBzHGSwUGnNHSb+wqYpL4wYa8+1en+jCPbzxcevd0R0kYSAAqVovlw7zkI7qhK34IoFZcXXMeiCHj1kXYm2Uh/ykiSNoaTZ5IsGQ92/omJXFWR7qsUyyZWfneMf1YXIhSLcs4SAAYIe9NuUi994J8UBmIv7PslSyLtf5DL+/s0PzRlZr2PHeuOsg3X+xbvutBW7EEHH7skSEAK/o8y/Fi2jbPRmbDLEUqEXzyo7EhXn0QOH0MK3NRuEzLsdxLL4cwnGcQUIBqM+4HMcZKNQTCTqOM4D0u3FedBPdniJyCsDzrY9TAEIluQcZP6/BYzud20Wquut8OhCRr6N5TZLwoqrecj7HOxc2VXG1HVjkgQSFJAcOP6/BYzuf23alBxGZjuM4m4srLsdxBo6tVFwHtvDYG4mf1+Cxnc9tW7JlNi7HcZxzxV8VHccZODZdcYnILSLylIg8IyIf2Ozj9xIR+bSInBSRR1fJJkTkHhF5uvVz51aO8VwQkQtF5F4ReVxEHhOR97fkA31uIpIXke+LyA9b5/VHLfnFInJf65n8BxFJXtvO2RI2VXGJSATgfwF4C4CrAbxTRK7ezDH0mM8A6FzD8gEA31TVywF8s/V50KgD+F1VvRrAzwD4zdZ9GvRzqwC4WVVfCeAaALeIyM8A+FMAH1XVywCcBfDeLRyjk4DNnnHdAOAZVX1WVasAPgfgtk0eQ89Q1W8BONMhvg3Ana3f7wTw1k0dVA9Q1WOq+oPW74sAngCwFwN+btpkJfl7prUpgJsB3NWSD9x5vRTZbMW1F8ChVZ8Pt2TbiRlVXUk8cxzAzFYO5nwRkf0ArgVwH7bBuYlIJCIPATgJ4B4APwYwp6orOXm24zO57XDj/AaiTZftwLptRWQEwBcA/I6qtpXVGdRzU9VYVa8BcAGabwBXbfGQnHNgsxXXEQAXrvp8QUu2nTghIrMA0Pp5covHc06ISAZNpfVZVb27Jd4W5wYAqjoH4F4ArwEwLiIrCQe24zO57dhsxXU/gMtbXpwsgF8E8OVNHsNG82UA7279/m4AX9rCsZwTIiIAPgXgCVX981V/GuhzE5FdIjLe+r0A4E1o2u/uBfD2VrOBO6+XIpu+AFVEbgXwMQARgE+r6p9s6gB6iIj8PYDXoxlJfwLABwH8I4DPA9iHZiaMd6hqpwG/rxGR1wL4NoBHAKwkZvo9NO1cA3tuIvLTaBrfIzS/tD+vqn8sIpeg6SiaAPBvAP6jqla2bqRON3zlvOM4A4cb5x3HGThccTmOM3C44nIcZ+BwxeU4zsDhistxnIHDFZfjOAOHKy7HcQYOV1yO4wwc/w8yGIVfHVMGBwAAAABJRU5ErkJggg==\n",
"text/plain": [
"<Figure size 432x288 with 2 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"plt.imshow(x)\n",
"plt.colorbar()"
]
},
{
"cell_type": "code",
"execution_count": 50,
"metadata": {},
"outputs": [],
"source": [
"x, y = np.gradient(x)"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": []
},
{
"cell_type": "code",
"execution_count": 60,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.colorbar.Colorbar at 0x7faf814ec810>"
]
},
"execution_count": 60,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAATwAAAD8CAYAAADqmhgGAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAgAElEQVR4nO2dd5hU5fXHP8cVsP5AagARCBJQLBg2YCQRBBHFKCYasaNBJfbeG2LD2GKMGrFiQbFExR5EsUZkUSygCAgKiCBNERUEz++PuUt29pzLzhYXdud8nmceZr7zvTP3zgzvvvc995wjqkoQBEE+sMG63oEgCILqIga8IAjyhhjwgiDIG2LAC4Igb4gBLwiCvCEGvCAI8oYY8IIg+FkRkT1FZKqITBeRc53nTxeRKSLyvoiMFZHWJZ4bKCLTktvASu9LXIcXBMHPhYgUAJ8AfYA5wATgYFWdUsKzGzBeVb8TkeOAnqo6QEQaAkVAIaDARKCLqi6p6P7EDC8Igp+TrsB0Vf1UVVcCDwH9SxpU9WVV/S55+BawZXK/LzBGVRcng9wYYM/K7MyGldlYRPYEbgQKgDtUddja/Bs33lQ3b9MwS9tqxRzfvNhKnzZv41pXOYfxq6+n+6+70EoTN+ziewus1GWDiUZb3mJjd/NPaWu0uvzoejfhO6M15wvX+wkdjNaEr1yv99l8sXBLxwld5tlj+3L7pkabTzN3+1Z8brSGi792ve/U39FoutL/OXZZZPdr5patHScsXtLYaC22sL+xxWzhbu99Xqs+q+fv1092v2a0td85QLupM4028f/8313L5rONthE/uN4ZS9pnC1/NQr9ZKK45R7YWUftr9JkHkyFr54ar6vASj1sCJQ9oDtBtLS85CHhuLdu2zHHXXCo84CVT1ZspMVUVkdElp6ql2bxNQw4oOi1Lu2XmGb75QSsNOH+Ia11II6ONfXYf/3XvsJI0KvK9DaxUtIn9Lb116a/czQ/lAaO1SBnECrH7cBFDXW8v7jbaYP7lehdhB4CL7rzW9RZdYY/tqqJDjHYDpxkN4ApONtrB9z/pejfe5yWj/TCnoeOEorvtfh1y7cWu98HH/mK04/Y/02ijGOBu7/2Wvjzml/5+/WD364D7/O/s0Z6HG036+L+7Ey6wn28Hprre/R97Nls4u9D1lYfvgME5eofAD6pa+TcFROQwMqevPari9Twqc0pb5lQ1CIKah5CZCeVyy4G5QKsSj7dMtOz3FNkduADYV1VXlGfb8lCZAS+n6aaIHCsiRSJS9P1XyyvxdkEQVAcbABvneMuBCUB7EWkrInWBg4DRJQ0ishNwG5nBbkGJp14A9hCRLURkC2CPRKswlVrDy4XkfH44QNPCVhESDoL1HAHqVNFrqeoqETmRzEBVANylqpNFZChQpKqjgWuAzYBHRATgc1XdV1UXi8hlZAZNgKGq6qzu505lBrxyTzdXUpfZWZtA07afud5h55vLdejirHMBbMhqox3eb7jjhEv7OWs//jITv+7/ut2es43WCrvIDLB0tV0EbFNgF68BjucWoz3BH13v0c5C5G6Mc71XcIHRDh50l+tdfpid8K92Ijc78a67vbe/ux821vV6n9m0F/01vHeu3cZoB/Co6x31O7s29wUtrLbaagAFBfa3dOLtf3O9V3Gq0Tbme9c7Ydx2RjsY/3vYjZeNlhYo2mP/7B/vW1ctdX3lofiUtqpQ1WeBZ0tpF5e4v/tatr0LUj6oClCZ41ozVSUz0B0E2BXuIAhqFFU5w1vfqPCAlzZVrbI9C4JgnVDVM7z1iUodlzdVDYKgZhMzvCAI8obiKG1tJAa8IAiyiBleFfH191vw9OQ/Z2nbdZrgege9MdJoZ3S/3PVe+/FFRuvS8TXX227yPKN92L+d672XI4z2L+ca9C9f8K/Ev67v8UZLi7a1f9amP9Xt96LrfdvJzDmG213vYG4z2mELHnO9VzW1UccX6Gu0RU42AsCRTgbI1ZzjejdxoplN9repaQB/4t9Gu4azXO+nzWxq14VcYbTLC2z0GmAsNmA4bMV5rrelc1HC8fVstB2g6/gPjDavm5PKA9zm/Ma8KD5APVZkPT6Z+a6vvNTWmVBtPa4gCCpIzPCCIMgbIkobBEHeEEGLIAjyhtp8SlutFY/rFm6vzYqyF6DH0Mf1/tUpd7QtfuWpAlYZ7aYv/MXyk1pcbbTd8dOfprCt0bxUK2+RGaCfc4liG2a5Xo/NWebqL9LbaAtSgiHeaxzOfa73LGwK1fWcbrSbnDJQAA9ykNF68Krr9YIOfVf7eeH131pptDodv3G9Jzf6h9F2xwZ/en/9irv9+Pq/NlraZzvVqUv4Adu73rqlggvgl+4C6MZ4o73oBFMAhpGdgnlU4Sd8VPRdperhdRJRpzqby44wsarKQ1UHMcMLgiCL2jzDiwEvCIIsImgRBEHeEDO8IAjyBiGitEEQ5AkC1Ml1ZLDxwvWaah3wOjCVf9MrSxtHT9c77uO9jHZaxytd79/n2qYnn7S0ETTw05/2m+pHB6/rYFPDVjlR2tkf+018ju54k9G+S/nb2ZhFRjtphp8u9ut2bxhtJya53gbYgpCHv+cXzzx7x0uN5qVapRXfbDN2gRW39ru0FbW2gb2mBc72QNfuNqLajhmu97rPbHS+TetZRjumvv/Zer+P79nE9S51ujwV4XciG8Aoo/VNqVa+EtslrQ1+4djS6Y+LuMH1lQcR2DAGvCAI8gERqOO0KK0NRCPuIAiyKJ7h5XLL7fVkTxGZKiLTRcT0bhCRXUXkHRFZJSIHlHputYhMSm6jS29bXmKGFwRBFiJQx+89XoHXyql/9efAkYBtIAzfq2rnqtmbGPCCIChN1V6It6Z/NYCIFPevXjPgqeqs5LmfquxdU6jUYYnILGAZsBpYVVaKSQGrzSK611EK4PqOxxnt7hVHud4nWtr+3/s96S8I1+1v05T6d/ATabzadQVOh7S0YMqlXGK0x1M6kXnpaZ+028r1tljxhdH61Bvjet9cYNPQlu24ues9k2uNtgPvGe20lIXxa3p/arS9UzoAdGCq0d5kF9e7h7O4v3VK0OL61jYVzqvJN+eN9u72z3bvZ7R3bv6d6/3vCXbikZa2dxy3Gu2yBX6AQ9616Z76Sz9b7N/ts4N7/8FPuSsXVTvgef2rbUHHdDYSkSIy4ZFhqvpEZXamKg5rN1VdWAWvEwTB+kLuI0PjZEAqZnjSi7qqaK2qc0Xkl8BLIvKBqvp/7XIgTmmDIMhGwLn6Ko2FZZzZlbt/dUlUdW7y76ciMg7YCVKm9zlQ2SitAv8RkYkicqxnEJFjRaRIRIoWfVV9lVmCIKggxae0udzKZk3/ahGpS6Z/dU7RVhHZQkTqJfcbA90hpWRSjlR2hve7ZLrZFBgjIh+ralY9oGR6Oxygc2FBjHhBsL4j4Fz7XCHS+leLyFCgSFVHi8hvgMeBLYB9RORSVe0EbAPclgQzNiCzhrfuBrwS080FIvI4mYiMXwAtCIKaQRWXS/H6V6vqxSXuTyBzqlt6uzchpcBgBalwAVAR2RTYQFWXJffHAENV9fm0bbYurK/XFe2cpT3Ofq53xEQbpWUz/3Uf7rCP0d5P+Zy+xUYob1hwvusd3PTvRpuEjczthx84mopNORsx3jku4MRutvhmBz5xvScsvtNo4jf8ol9n26HsRk5xvRdiu8J53bK25313+3uwUfTTP7TRScCdQTzS/g+u9c+znzZat1bjXO/4mT1zeq/PWjRxtx/Hbkabjt/VzivgOWy1ua4WgCkFtpjszZzgekfOP8RoRzWzKW8Ay0r9nscUDmFx0cxKFQAt3Ey0KMdhRt7KnwKgzYDHRaT4dUaubbALgqAGUUtTyyo84CUXEu5YhfsSBMH6QC2uAFpLDysIggoTA14QBHlDFUZp1zeqdcD7CTH14NJqgl3fxaYInYLtSAV+zTevvhz4Xbx+2NS1cj5XGO1gHjLaAEcD2NnpPtWrm12ABzjCSUnqtfwl1zu/YVOjXbboKtdbhF1PvtbN0YaH3xtotAE72jpuf+YRd3uvbtzpi/ygxfE9rjNaWs03Dy/FD6Du/31ttDMb2ZS5GSmBCO/3OHSx/9nu3vApoz1QYAMOmffb2mj3P3OM69Vf2ZjD0c1sbUWwqZk/VkVx9pjhBUGQN8SAFwRB3lC+1LIaRQx4QRBkEzO8IAjyhghaBEGQN8QMr2pYTCNGcmiW5nVzAvgE23WsAUtcb+n0GoCjV9zhet+sZ4tM3r2pjU4C7M5Yox3OvUZr9+Q8d/uH+9uUt7RuW3tjo7fXb2oj1QBXcIHRWrS3RUHBTwNrynzXu+WO04zWmXeN1iqrnuP/uGe27fImv0hJXbR1RRm+4+GudUCre4w2mNtc738/7mW0T7rbFL/H3jjM3f7M7jaiu2dDP3XwLKdgalqBV0+/Ze8jXe80m1bKu05KI8DFXJb1+IyU77ZcxIAXBEFeUUtHhlp6WEEQVJiI0gZBkDfEKW0QBHlDLY7SVrgeXkVoUriV/qnojCxtBXVdr1eXbHdedL0b853RuvG2613pvN+L7O56X+P3RuvndOFKO4a09DaPdkw32r7Y1CWALV9wXjelVuCF3W2tv1EMcL2f/McujA/ew9YEPDklxW+F879kJH6qlZfylvb9tmGW0eqxwvXOzmqfkMELGKQFPbzX9YJiAL9yOq9d6KQjAhyNDaItopHr9fb3CCdYBjaANafwQFYUTa5cPbwWokV+1ptBhuZPPbwgCGojcUobBEHeUIsHvMp2LQuCoDZSkOMtB0RkTxGZKiLTRcTUwBeRXUXkHRFZJSIHlHpuoIhMS27+BbPloJaO40EQVJgqnOGJSAFwM9AHmANMEJHRpbqPfQ4cCdl1y0SkIXAJUEimJezEZFs/AyEHyjwsEbkL+AOwQFW3K7Ejo4A2wCzgwFx2oj5fm0X/s7DNawBGYAfzFvjZBL9ZMcFoXkYF+IEAr0Ye+FkRdZ1F7ds+PtXdXpbbgNA2Xd5xvedwtdF2w6+Hd0vfI432ATu43vJwyB62OdCDD/zFaJsf6n9ety0fbLQnN+3ves90shRO5J+u9wOnIVNahs7xy239vYJNbe28w76wzY0ALmpxntEKmeh6J7GT0aZgm/WAH5C522l6BNDMyZb4Pa+53oc4KOvxEXzq+spF1UZpuwLTk5YQiMhDQH9K9JdV1VnJcz+V2rYvMEZVFyfPjwH2BB6s6M7kckp7T/ImJTkXGKuq7YGxyeMgCGoD5WvE3VhEikrcji31ai0hKxdxTqLlQmW2dSlzhqeqr4pIm1Jyf6Bncn8EMA44pzI7EgTBekL5TmkX1qTLUioatGimqsUZ81+SadnoIiLHFo/+33y1soJvFwRBtVGcWlY1QYu5kHVx5JaJ9nNv61LpKK1mrlxOvXpZVYeraqGqFv5fE/8C3SAI1iPKd0pbFhOA9iLSVkTqAgcBo3PckxeAPURkCxHZAtgj0SpMRQe8+SLSHCD5d0FldiIIgvUIATbK8VYGqroKOJHMQPUR8LCqThaRoSKyL4CI/EZE5gB/Bm4TkcnJtouBy8gMmhOAocUBjIpS0eDzaGAgMCz598lcNppLS84juwPUPzjF9d6Cra026o0jXW+D7rZrWT380+duTiexmznB9V6PX4+uNEd39DtK6TSb4XMRNgoIcORUG3Xcv8P9rve4B0cY7Yd9/X3byCldd9mnfhcu+chO1K87w34PQ1dc7G4/fFObj9R76puu99s29nxoej3b2QvgvWd2NtqUvf1o6J8+es5opzQbbrS/tTrR3b4Dnxhtv1f8ScVxPa432gUpqWVvYq8aOIq7Xe/La5bH/4f3/wFsett3+LUZy0UVV0tR1Wch+/IMVb24xP0J4BQBzDx3F3BXVe1LLpelPEgmQNE4GYUvITPQPSwig4DPgAOraoeCIFjH1OJMi1yitAenPNW7ivclCIL1hXwd8IIgyDOiAGgQBHlDPp/SViU/LN6Ejx74dZZWcKhN+wF4eLJNLdu++weu16sld1BK6tHLTp29BtigB/j1yho5Ne7SGvPgxJO64y/in93hUqN5ARaAaQfb9d07ONr1HtHR1lEb2tEPOozr181o37Ox0UbU83O4vYX1tMDLoVjdSzcDaLG3TSksHfwqxkvnu6zVmUa7aMY17vbMsYGm1j0+9r0O2+P/Rq/F7sP+d9raigBXDzrJaI9ygOO0dfZ+rIr/0rW4AGgtHceDIKgwMcMLgiBviAEvCIK8IQa8IAjyiojSBkGQF8QMr2qo33AJvz/0kSwtLfp0YCebPpXWPaqt09Xq8OV+dPBxpyBlWjHJLu9NMdqgHW2RSq8DF8D5c24w2h+6PeI44amPbbLKIR1tQU6A/Wc8Y8Uv/UZVU7rbFKwZtHO9PcbbTm/ndBtiNK9AJfjFL2dP+5Xr/Vt7m9qV1sWrMQuN9tKKXq53s0k26r+0RwOj/b3dX93t57ezhX+unHmZ65UNbPriyF/4XdpurGdTKIsG+b+bcybbVMVfdprseksXxa3Lj66vXESUNgiCvCFmeEEQ5A0x4AVBkDfEgBcEQT6hEaWtPFt/+ylPvZG9OH9Sd9utC6CpszB+L4e73nkz7SL88217uF6vQ9mf+Lfr9T6dmbQx2lLsojjAwfvbMl7XcJbrPafjEKMdz82ud2o7Gwh458Pfud4jsKll7ztdwAAO6Haf0U7CLqD3WTTG3X5hg8ZGu6i9X//P6/iV1j3uRXY3Wt16fr3Dh0/Zx2gHjrCd6v47sLO7/c5F7xnt+MLrXO91Ti1Hb1/BD86l1WG8u5PtZpb2uykdMJP04uM5oxvAyhyKe9ZEYoYXBEEWKrCqINdi6KU7K67fxIAXBEEWKsLqDXMdGmpWY65KN/EJgqD2sbqgIKdbLojIniIyVUSmi4jpYS0i9URkVPL8+OK2sCLSRkS+F5FJye1flT2umOEFQZCFIqyuotwyESkAbgb6kGmkPUFERqtqyav6BwFLVHVrETkIuBoYkDw3Q1X9BdcKUOYMT0TuEpEFIvJhCW2IiMwtMfL2q6odCoJg3aIIqyjI6ZYDXYHpqvqpqq4EHgJKpzv1B4pTqx4FeouInzpUSXKZ4d0D/BNMuO8GVfUrNqbw5WZNuap7durNPFq43kdfsRHZtj1mud6v2m5mtFfZ1fU+ix2b67LC9d7TaYDR9nFaap7CP9zt/7XCpi8V1fPTiV50WoR4na4AnmcvozX5zbeuV16xUbtterzjeg9hpNHO5m9G+7KR33d9FPbzSotgT8VGmv/Kba7XS/3z3gvg9NG3WnE7K/0V/+xoxjY24t81pRBrIUVGe/b+/V3vSYfZqxEe54+u99ki+xpdCl9zvWPok/X4hpTfcnlQhJW555Y1FpGSH8RwVS3ZJq4lULJ33hygdKXZNR5VXSUiX8OaPMO2IvIu8A1woar6H0SO5NLE59Xic+ogCGo/5TylXaiq/l/xyjMP2EpVF4lIF+AJEemkqt9U9AUrE7Q4UUTeT055t6jE6wRBsJ6xmoKcbjkwF2hV4vGWieZ6RGRDoD6wSFVXqOoiAFWdCMwA59SgHFR0wLsVaAd0JjMK+1dmAiJyrIgUiUjR8q++r+DbBUFQXVTxGt4EoL2ItBWRusBBYNaFRgPFjVIOAF5SVRWRJknQAxH5JdAe+LQyx1ahKK2qrkmDEJHbgafX4h0ODAfYsrBZ5S8DD4LgZyVzSls1F3Aka3InAi+QKSt6l6pOFpGhQJGqjgbuBO4TkelkWl8dlGy+KzBURH4kc4XzX1XVaY2VOxU6KhFprqrzkod/BD5cm7+YhTTmbo7M0m7nWNd7eY8zjHbmcn8i2XZT2zVsd8a6Xq/m2qVc4np78rKzvU2f8havAXau95bRPnzjN653r+42vS2tzt5QLjLari38tdw/tLD1956aaWvvAXzctrV9r0W2w9kTjfzF9oXOZ/M9m7jeZ/iD0fbkeddbgK1xdy9HuN6e+44zmt5jA36Xd7G/L4CL3rJxuJN72/Q6gGO43WhH9XvQ9W7tdLbz6jgCvF74a6PNZivXeyFXZD3+AluDsbxkghZ1K/06a15P9Vng2VLaxSXu/wD82dnuMeCxKtsRchjwRORBoCeZaMwc4BKgp4h0BhSYBQyuyp0KgmDdoZDr6WqNI5co7cGO7JfiDYKgFlB1p7TrG7XzqIIgqDBVmWmxvhEDXhAEhhjwgiDIC2KGV0Vs9+MUir7IzgPu18IPwjw39k9Gu6jNNa53YDubJlTIRNc7cr7tKvX7Zn6Es7+TRuYx4evfuvom9b8z2uPd/Qjn+aWibQDXcqbr9SJ+X6Sk6HkFQF9p29X1jsR+Np83stHBE7Gd2wBudApibjX/c9c7ajMnbW9TW6gT4NTP7Pvd0/ow19tlrO00N+PI5ka7sMiP+C/tbVPh0jrreZ/tBi/7V1412d9+DgfwqOs9d9Ewo53QyC8GW7og7UrquL7yoAgramnbspjhBUGQRczwgiDIG2LAC4Igr8jb6/CCIMgvqjK1bH2jWo/quzob8U6Ltlnac+/Z4ATA2N62FlxaDTQvNexqznG9xzezi79p3dC25wOjtWGm0S6sP8Td/mjuMNohK2zNOYAX69luV1dwoetdxuZG23mm7bYFMLjt3402/LPjXW/X1m8a7R+cbLTZWcUv/of3/az+xtYqBLi12UCjecEYgJWtbZqTlyII0L+3Te1qN3We0d4u9Du3XbvApu3t19T/zu772KZFXre1/51dik3Re4p9Xe+YRn2MNo6ervfJz0oFVFaWq0SlS5zSBkGQN2SitFWXS7s+EQNeEARZxCltEAR5RZzSBkGQF8QaXhAEeUMMeFXER0s60eWx7GKZx+1/vet9nP2M5kVNAf6K7VSVVlzxptE2envpvme73vteslG4W3vZ6OKVXOBu/7dRNno8a0BT11tEF6P1n/0f17tLK6e4aVsrAfyTE402v7W/D153sEMmP2G0bTr5Xc+8Qp33t/e7eHkR2bPwUwffWNHdaAPrjXCc8MTHNj1uv47W1/UB/7c071CbWnYCt7jexR03Mtp5TjQW4I/Yz7EvL7her2jqJS/Z7nEAQ+qV6oa2svLdDSO1LAiCvKE2z/Aq07UsCIJaShV2LUNE9hSRqSIyXUTOdZ6vJyKjkufHl2wLKyLnJfpUEelb2eOKGV4QBFkUdy2rCpKuYzcDfcg04Z4gIqNVtWRZm0HAElXdWkQOAq4GBojItmQa+nQCWgAvisivVNWuneRIzPCCIMii+Dq8XG450BWYrqqfqupK4CGgfylPf6B4UfZRoLeISKI/lPSnnQlMT16vwuTSxKcVcC/QjEx/j+GqeqOINARGAW3INPI5UFWXrPXFNga2y5YuT1nwfwE7e90Ev6+t16HsvOVXud5z9r3aaLtgU6oAjuxlF6snpnQS8xgywAZIOn79sevdu/6zRpOfVrrekbbBE29iU/EAruR8ozVmkeudxE5GO6+TXYRvgP81eylnd3C06+3AVKNtgq0fCDCqnk1Zm0oH13tjRxto8mojvnmoDYQA7O10HO2LHzw6HRtw+2CS///x+c49jPYqu7reF9jDaLv28ms23lKqC+CwzWa5vvJSjjW8xiJSMhI5PGnNWkxLYHaJx3OAbqVeY40naev4NdAo0d8qtW3LXHfMI5chehVwhqq+IyKbAxNFZAxwJDBWVYcl5+XnQkoCaxAENYZytmlcqKq5zwLWMWWe0qrqPFV9J7m/DPiIzChbcho6ApzrSIIgqHEUr+HlcsuBuZA19d8y0VyPiGwI1AcW5bhtuSjXGl4SPdkJGA80K9GM+0syp7zeNseKSJGIFLHkq0rsahAE1UEVr+FNANqLSFsRqUsmCFG6d8JooPgC1wOAl1RVE/2gJIrbFmgPvF2ZY8s5Sisim5HpAn6qqn6TWVPMoKoqIm4x/+R8fjiAbFfoF/wPgmC9oqquw0vW5E4EXgAKgLtUdbKIDAWKVHU0mT7X94nIdGAxmUGRxPcwMIXM0toJlYnQAkhmIC3DJFIHeBp4QVWvT7SpQE9VnScizYFxquqvJCfsWLihvlCUXR9tEp1db13sgn3auoK3+Ls5y1yvV0tuBu1cr7cw7u1DT8a5298y8wyj9WvrNy161/kchnGe651VqnELQKOUQIQXSLh6wRDX27npf4026SXboOikXjbwA1DP+c76YYMxAFc4wZQdUjJp7l19hNEuL/CDXacs+ofRLm401GhpTXx6F9pGQgv8kxc+HPsbq/X2f0v7Og2h0oI/3uf43zd6uV5tm51ZUbgXFL2nlUq3aFrYSv9cdFpO3lvkjIm1ag0vCQ/fCXxUPNgllJyGDgSerPrdC4KguqniNbz1ilxOabsDhwMfiMikRDsfGAY8LCKDgM+AA3+eXQyCoDrJRGnzNJdWVV8H0qbIvat2d4IgWNfU5lzaSC0LgsAQA14QBHlBVebSrm9U64C3mIaMJLvL0hJs/TGADZ3aamk1up5wrnl+72tbXw5gox++NZrXyQz8emW3MdhoBaxyt+/S1qYDnYsf4VzEmUZrxELXe+RcWwtur5Z+NPSZxbYe3UVN/ejvk04Xra962a5jaddfXTfV1hX8oIPfHWxbphitAUtd7zUFZxntL8/a7mQAy/rZKHwzFtj3L7TpZgBNHe9p3OB6V/S2EfuHUjrrdWO80dLSKnd2vJd1t7+PzL5dmfV4dh3/t1weoqdFEAR5QzlTy2oUMeAFQZBFnNIGQZBXxCltEAR5QVyWUkV8z8a8T/YitreYC37AYFQmxc7wndP05Mr6fmOepvXnG+0A/HSvoVxktPGzexrtb61soxyA0/i70V7l9673V3xitIOXP+R6t2tpU7AOYaTrfaGhrcM2jxaudy+eN9o+2FSrWybZlDmAWyefbrTpHbZ2vdfjePG9nXnXaJf38/dhMLcZzUtf/Gjqr93tL+/Qz2g34KdZvfaFrVt3e4vDXK9Xs/GTlJp+OznHm1YrsHTw56mUlMryEANeEAR5RazhBUGQF/zEBvmbWhYEQf4Rp7RBEOQFsYYXBEHeoMQaXpWwlX7OzSuyI5qH1rvf9Zana9k+TnHFIRP9FC51aiPuVwz8PnoAABfMSURBVOhHOE/mJqPNbdXQaF6aFMB405wp/fqm+z603bYGbOdf7d6Tl412rZOaBvDeXBuhfKzln1xvG2YZzSuYemlnPwL+7Tb2P8mwlL5OXsFS770ABnKv0dJmIF503/Pq934BoMMZbrRLucT17tPiYaMtopHrfdC5wmDHFe+73vfq7WC01tP89gintc9OLVuG/5rlI1LLgiDIE2rzKW004g6CIAtFWEHdnG6VQUQaisgYEZmW/LtFim9g4pkmIgNL6ONEZKqITEpuTct6zxjwgiDIooq7lq2Nc8n0tm4PjE0eZyEiDYFLyDTv7gpcUmpgPFRVOyc3W+qmFDHgBUFgWE1BTrdKkktv677AGFVdrKpLgDHAnhV9wzKHaBFpBdxLpu+sAsNV9UYRGQIcAxSvpp6vqn5RtoQvpRnX1huYpT352QGu98LWdqG4g5N+BXDAcpsaNrLzoa53cIFN93pyxMGu918DjzPaifzTaEdzh7v95U9eabRx/W0gA4BNrXQ8t7hWL6CTtuB/cEsbkNn/Gf9rGrm3/b1NYVujFTi1CgHurneU0dK+M6/zmvdeAK+84vy+f+Fa2beD7SXlfT+7bP+mu30/njHay/R0vUXYmotXcKHrPYwHjNatnp9WOcqpqbdDe7+j29/HlqptuMxPkywP5VzDaywiRSUeD09as+ZCLr2tWwKzSzyek2jF3C0iq8m0kL1cy2jDmMucdBVwhqq+IyKbAxNFZEzy3A2qem0OrxEEQQ1BEVb/lPOAt3BtbRpF5EX8P09Z1U/X1tt6LRyqqnOTcekxMs3GbEi/BLk08ZkHzEvuLxORj8geYYMgqEXoT8KKH6omtUxVd097TkTmi0jzEr2tvTW4uZA1xd4SMo2gVXVu8u8yERlJZo1vrQNeudbwRKQNsBOsKXFyooi8LyJ3rSXCcqyIFIlI0Xdf+dfRBUGw/qAqrF5VkNOtkuTS2/oFYA8R2SIZY/YAXhCRDUWkMYCI1AH+AHxY1hvmPOCJyGZkpo2nquo3wK1AO6AzmRmg28pdVYeraqGqFm7SZONc3y4IgnWFUl0D3jCgj4hMA3ZPHiMihSJyB4CqLgYuAyYkt6GJVo/MwPc+MInMTPD2st4wp7hyMoI+Bjygqv9OdmR+iedvB57O8SCDIFiPURVW/fjzX3isqotwelurahFwdInHdwF3lfIsBydqVAa5RGkFuBP4SFWvL6E3LxFh+SM5TCeXs5lJt/pta9vZC2CGUwzyHmwUEOD7Te3M8RAnKgZ+B6vDBvp/GPowxmgt+MJofRe/4m5/YH/bXcyLwAH0rGtf40r87mK78qrR0gpEehHV3fY+3PWOxEa2n55ho+j6pX9icHR3m4p369cnud4/1H/CaCfxD9d7Yo+/Ge2Or492nPDlM7802uq97c/8XwW2+xzATkwy2vdOgVmA552rI/7E467X+35GcITr9X77+37tne3Bh73bZT0+cPM5rq98CD+trp1JWLkcVXcy0Y8PRKT413A+cLCIdCZzqcoscPoXBkFQ81Cg8qer6yW5RGlfB7xM67VecxcEQQ3lJ4Ef8neGFwRBvuH3lq/xxIAXBEE2mYJ4tZJqHfCa8JVJ83mTXVzvB6W6mwGsSKmzv4QGRptNK9frpS+9uHgf17vBspVG69zaLmr/ueF97vZex6+rU+rD9Wr5otGuXe3XuDu5wC7ut8rKvvkfXsctb7/SXmODzexie7fu49ztt8emP/Wq/x/XO9Xp2OVtD36gaWV9/7cwZW/7/a50qnqMxb8e9hZOMNqs8du43q+6bWa0T1/o5Hov62u/yyuzkw3WcMdMG+jZpa2fCjezVIreSvy6eeUiBrwgCPIGBX5c1zvx8xADXhAE2SiwYl3vxM9DDHhBEGQTp7RBEOQNMeAFQZA3xIBXNSykEfeSndb0xOxDXO/zrXoY7XH+6Hq3ZobRdstUkDHc7aSnLWxoo20AAxvawpEtmGe0qx4Y6m6/y6E2svZ7/FS6pU6k+YKCK1zvGRNtYVBt7nfh2qeFjcjegZ+W9T02Re/TZm2N1uYVv5L2Pj3se6Ud71Autq870X/d/3T5vdEOwe801+PDt402YLt7jJYW8X+J3Ywm3/pl2vZwint80te/OqC/4/U6mQEMbmuL1O6bElm/rNTn+AW2+125iQEvCIK8Iga8IAjygp+AH9b1Tvw8xIAXBEE2cUobBEHeUIsHPCmjyU+V0qiwjfYtyl5kvTwlvabdZ58brVfrF1xvPWwK2LOT9ne9F3X2a8x5XDbpKqP17Pyc0bxOVwAnL7cdzjpt6pcNbMxCo+3Po67XS8fz0q8AbltwqtFeb/pr17vDapva9WiBrYe3kEbu9n2x38/KlOBAU+Yb7Tb+6npfwwYtLk/pDtZtuQ1ajN+0q9HSAjef8CujFc33e9QUbGhrDQ5u9C/XO48WRuuZEljzvt8GLHW9B5T6jRxX+D5Ti771I1g5Iq0LlQuKyjYCDJaJa2vis74RM7wgCCy1dIYXjbiDIMim+JQ2l1slEJGGIjJGRKYl/6Y1AnteRJaKyNOl9LYiMl5EpovIKBGxVSJKEQNeEATZ/AR8n+OtcpwLjFXV9sDY5LHHNYDXl+BqMr2xtwaWAIPKesMY8IIgyEaB1TneKkd/oLjxywhgP3d3VMcCy0pqSa+dXrBmETN1+5Lk0sRnI+BVMm3RNgQeVdVLRKQt8BDQCJgIHK6qNnpQgoYsNs11Tud63/xtHSNdyiWuda/lzxvt9539Omx1nTIQm2d/lmv4d+e9jOZlI+yQUsftzU1/a7QZC7ZzvXObNjSaVxMQ4OwPbTAkrWn7lE62Ply3NW2FsxlXYLMMvBp1adsf6jROeu/GnV3va6fYhlNnco3r9Ron3cLxrrfnHXbfzj7lUqOlfeddnWMb07iP6321wAZTrk6ZpLw2ew+jyQL/O7uuiz22hTR2vaU/mzpVVdcp99PVxiJSMsIxXFWH57htsxKNwL4EmuX8rplxZ6mqFu/pHKBlWRvlErRYAfRS1W+Tdo2vi8hzwOlkppMPici/yEwnby3HDgdBsD5SvstSFq4tSisiLwK/cJ7KujxDVVXS/mpXIbk08VHg2+RhneSmZKaTxYmwI4AhxIAXBDWfKrwOT1X90tKAiMwvbvcqIs0h5doqn0VAAxHZMJnlbUmmGfdayWkNT0QKkhaNC4AxwAxynE6KyLEiUiQiRV9/tdYz3iAI1geKU8tyuVWO0cDA5P5AcCospJBMxF4Gii8UzWn7nAY8VV2tqp3JjKJdgY7l2LHhqlqoqoX1m5QZNQ6CYH2gGi5LAYYBfURkGrB78hgRKRSRNaWKROQ14BGgt4jMEZG+yVPnAKeLyHQya3p3lvWG5brwWFWXisjLwG+pwHQyCIIaQDWllqnqIqC3oxfB/1JhVNVGhzL6p2QmYDmTS5S2CfBjMthtDPQhc/1L8XTyIXKcTs6iNcdwe5Z2Lle73ie2sHXyXkk5tvc23cFoB/OQ6z2Ta42231x/11e3tB/PHk76lFfbDfwad/s39dPFZmHrzjVgievdfjsbOT1wTXQ/m4cnDzTaQ538OmxercBGTsrbX7jL3f4azjLawlP8NLRjS/0OAPql9HbflilGS4vYdzvFRlnPeMzWD/x+Tz/7qmhTm3Z3d4H9XMBPsXt1sY3GAu4A8u12Ba51mNPZbldedb2l0+6+TfkMy0WeN/FpDowQkQIyp8APq+rTIjIFeEhELgfeJYfpZBAENYDi6/BqIblEad8HdnL0ck8ngyCoAdTiailRPCAIgmyUqkgbWy+JAS8Igmxq8SlttdbD26Lwl9q7KLvhzUzauN4TsAvN81MyT75z0r068InrfcJJt+vNi673AB4z2gKaGu0SbOoSQHMnJWos/nWYbZhltLQGOF59OC/9CmAKuaeWveq87lHcY7TdUz4vL13LW4AHmIi9OH/KCruvaTxZr7+re5/NamxwIC1tz0vh8rYH2AXbpGkG7VxvK2Yb7SL85k+N7rTTqx0HveV635tYKnXv8EJ0SlHl6uE1KFR2zbEe3lNRDy8IgppMrOEFQZA35PllKUEQ5Bu1dA0vBrwgCLKJNo1BEOQNcUpbNfyCL01q1xEpKVGDpo402mEdbDoSwP0jjjHadQP9ApFNnQo0m6RcdPSL+V8arWezl432CH92tz+Rm412JHe7Xq/z2ikv+XUUx/ey13unFbQs3dUKYEVKJ7F3Jv/OaJt0sp9NWpR3T2wh1vlf26g2wKH17fe7Sz0b9QSYSgej9XzP34cTd/yb0W56xUaKu/UY527fjhlGW8bmrncTvjOaV0AU4MoVtjvf9Hpbu95eg5422ifOZwCwXZcJ2a+5yXLXVy5q8WUpMcMLgsASUdogCPKCuCwlCIK8IYIWQRDkDTHDqxp+pI5JD3sLv6vVeR2GGS0tHejggbY+W1+nbh34aWR3/K/WYBY/PbGp0e4dbOvLNSzy/xweXXiH0bzUJ4BuvG20w3v5QYu2ThpaWtqdl96W9jlqgc1I+qfT6vNet0Wov4h/ff3TXe9stjLaTryb4m1ltHt2HOB6Zzmpig12nme0tN+Hl4rXhpmud8hkW8vxzk62jiPAafVuMNqb7OJ6T+Ymoz1LP9d724enZj0urKqk/xjwgiDIC2rxZSnRiDsIgmyqqRG3iDQUkTEiMi35d4sU3/MislREni6l3yMiM0VkUnLrXNZ7xoAXBEE2xWt4P38Tn3OBsaraHhibPPa4BlLWUeAsVe2c3CaV9YYx4AVBkM1PZAqA5nKrHP1hTebBCHBqtwGqOhZSrqwvJ2UOeCKykYi8LSLvichkEbk00cs9nQyCoIZQDae0QDNVLY4ofQkpkbe1c4WIvC8iN4iIn0JUglyCFiuAXqr6rYjUAV4XkeeS585SVb8Nl8NnTtey97EdxwBuu/9Uo51z2BDXe/EKW0hxZD0/WuZF4bwilwCNBy8yWssvFhtt4/bWB35Htks+tKlPAC9t91ujpRX1PIQHjHat0zEM4G26Ga2uk8YG8E7HbYzmRbA7MNXd/kZOMdr0lIKYJzhpd4c6xwUwbuZeRjuz7WWutwFLjbZtPdv17HinwCxAKz432mvs6np36zTOaGmR16OxEfsu+EU2l2KXstI6up25XfbnMGfjW11fucm9LnBjESl5IMNVdc3lBSLyIvALZ7usXDtVVREpbzXi88gMlHWB4WT61PpVVRNyaeKjwLfJwzrJrfrKJAdBsD6zcG0Vj1XVL/ENiMh8EWmuqvNEpDk4ie5rocTscIWI3A2cWdY2Oa3hiUiBiExKdmiMqhZnSJc5nRSRY0WkSESKfvrKzo6CIMhbRpPpaQ059rYuSTJIIiJCZv3vw7K2yWnAU9XVqtoZ2BLoKiLbkZlOdgR+AzQEv3mBqg5X1UJVLdygScOcDiQIgrxgGNBHRKYBuyePEZFCEVmzBiAirwGPAL1FZI6I9E2eekBEPgA+ABoDl5f1huW68FhVl4rIy8Ceqlpc5ynn6WQQBDWB4jDtz4uqLgJ6O3oR/G/xWFXd9CRV7VXe9yxzwBORJsCPyWC3MdAHuLrEuXfO08kmLGBwqcXqgSn18G48zC6AL6KR632tnv08HuePrve5yX8yWkEnP9x0lFO7TtT+EHrUH+duP+ROG7T4YlAL1zv8M1u/777WB7ves7nGaC984Ub0ObKFXZwfxnmudygXG81Lq/qOTXLevnNKutgs2hpte953vWPa2jp910290PVO7GCDUktpYH10cbf30tiOWPyw6323oQ3ytGO6670CWw+vIOVCts3XLJn/j7QahKV/50t40PWVj9qbapHLDK85MEJECsicAj+sqk+LyEvJYCjAJOCvP+N+BkFQbdTe6gG5RGnfB3Zy9HJPJ4MgqAnk9wwvCIK8Iga8IAjyBqU6ghbrghjwgiAoRR6v4VUlTVhoUmyGjLeRTIBru802Wmf8Ygg3cJrR/jO+v+s9rtv1RtsFv1uW97rPtdzTaGmdqpYNst2ufpOSTnRga9vF6zyucr3nOClraTHyo1rcY7S7Ocr1tsJ+5t+xsdG6p3xe5662RVsvL/CjqadzndF2Z6zrPYUbjTarg98NbSL2on+v69l1N/v71eME23ltYkM/orvEif56UV6AE5xUtrS0u1s4wWhbp0R/G7Ak6/GXVdJuLE5pgyDIG2KGFwRB3hAzvCAI8oaY4QVBkDdUT2rZuqBaB7z3V+7Ilp+9laVd1s1PwX2KfYx2OPe53sedQql7dPMLL3h12LxFbYCmTrWaBzjUaF7aEPipaROcRXWA7510LS9VC+BFm37ISS1tvTWAiU79Py8YA9CYhUb7aOyvrdbbvibAwq+3NFpRw+1c7wBGGW2nlKCUV2MurfOal8rmhWie6n2gu/397G+0Y772u8cNqG+Pwes+BzAfG2R5in1d7+Hca7SClGBE6Q5nQ8pXYSmFOKUNgiCviFPaIAjygpjhBUGQN8SAFwRB3hBR2iAI8oaI0lYNqwWW1smSLjzLphgBXDTzWqPdNOEk1/vSjL2N9lo7Pxra6bpPjTbsDL//732PHGs0OXC50Tb71moAgwptupg+KK53RGcbNdyNl13vsXKEo45zvVPURlS/OmYr17v57bb155W9bUT3/M/8zmvyhNW6bv2B6+2/ty1U2YFPXO/rk/tYbSO/N8yV7U43mhfl/Wab54wG8F+nuK4XjQU4ANuwr15KR7jbGGy018UWEAXYXHNvwfrsg9lR5ZuqpG1MnNIGQZA31N5T2pya+ARBkE8Uz/ByuVUcEWkoImNEZFryr2nIKyKdReS/IjI56ZA4oMRzbUVkvIhMF5FRIlK3rPeMAS8IglIUz/ByuVWKc4GxqtoeGJs8Ls13wBGq2gnYE/i7iBSXqbkauEFVtwaWAIPKesMY8IIgKEVx0CKXW6XoD2u6eI0AmzKlqp+o6rTk/hdkemM3SZqH9YI1C6nu9qURVa3sTueMiHwFfJY8bAxOLlPNJ46r5lGbjq21qjapzAuIyPNkPpNc2Aj4ocTj4arq5+LZ91mqqg2S+wIsKX6c4u9KZmDrRKYX9lvJ7A4RaQU8p6p+LmNCtQYtSn4RIlKkqn4otQYTx1XzqM3HVhFU1Va5rSAi8iLwC+eprAR0VVURSZ19iUhz4D5goKr+lBkfy09EaYMg+NlQVf/6IUBE5pfob90c/MoHIvJ/wDPABapaXH1kEdBARDZU1VXAlsDcsvYn1vCCIFhXjAYGJvcHAqbEURJ5fRy4V1XXXPiombW4l4ED1rZ9adblgJfTeX4NJI6r5lGbj219ZhjQR0SmAbsnjxGRQhEprnd2ILArcKSITEpunZPnzgFOF5HpQCPgzrLesFqDFkEQBOuSOKUNgiBviAEvCIK8odoHPBHZU0SmJukgftZ+DUFE7hKRBSLyYQmtzHSZ9R0RaSUiL4vIlCSl55REr9HHJiIbicjbIvJeclyXJnq5U5SCmkm1DngiUgDcDOwFbAscLCJ+g4SawT1k0l1Kkku6zPrOKuAMVd0W2Bk4IfmeavqxrQB6qeqOQGdgTxHZmQqkKAU1k+qe4XUFpqvqp6q6EniITHpJjURVXwVKF+QpM11mfUdV56nqO8n9ZcBHQEtq+LFphm+Th3WSm1KBFKWgZlLdA15LYHaJx3MSrTbRTFXnJfe/BJqty52pLCLSBtgJGE8tODYRKRCRSWQuch0DzACWJhevQu38TQYJEbT4GUkujqyx1/2IyGbAY8CpqvpNyedq6rGp6mpV7UzmyvyuQMd1vEtBNVLdA95coFWJxzmlg9Qw5idpMsX5f1XRKLTaEZE6ZAa7B1T134lcK44NQFWXkrlS/7ckKUrJU7XxNxkkVPeANwFon0TF6gIHkUkvqU2UmS6zvpNUrrgT+EhVry/xVI0+NhFpUlxLTUQ2BvqQWZ8sd4pSUDOp9kwLEekH/B0oAO5S1SuqdQeqEBF5EOhJppTOfOAS4AngYWArMqWwDlTVKuk0UF2IyO+A14APyBRHAzifzDpejT02EdmBTFCigMwf+4dVdaiI/JJMAK0h8C5wmKquWHd7GvxcRGpZEAR5QwQtgiDIG2LAC4Igb4gBLwiCvCEGvCAI8oYY8IIgyBtiwAuCIG+IAS8Igrzh/wHKQsRJG3jzpgAAAABJRU5ErkJggg==\n",
"text/plain": [
"<Figure size 432x288 with 2 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"plt.imshow(x, cmap=\"jet\", norm=None, vmax=.2, vmin=-.2)\n",
"plt.colorbar()"
]
},
{
"cell_type": "code",
"execution_count": 59,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.colorbar.Colorbar at 0x7faf815aef90>"
]
},
"execution_count": 59,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAATwAAAD8CAYAAADqmhgGAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAgAElEQVR4nO2de7zVc/b/n8vRhTTlKBRNF6IhU3RGxq0mxM+gqSnKRBlhhmZyv85IuYwI04xcmiY1LiFqJn0ZUuT21dgRkSiV6aJSKSlKWb8/9uf4nn3W2rVP5zg6Z6/n47Efnf3ar8/en8+5vHt/3uu91hJVJQiCIB/Y6fs+gSAIgsoiBrwgCPKGGPCCIMgbYsALgiBviAEvCIK8IQa8IAjyhhjwgiD4ThGRk0TkAxGZJyJXO69fKiKzReQdEZkiIk1LvNZHROYmjz7lPpfYhxcEwXeFiBQAHwInAIuBN4Beqjq7hOdnwHRV3SAivwU6quoZIlIIpIAiQIEZQDtV/Wx7zydmeEEQfJccDsxT1fmqugl4FOhS0qCqL6jqhuTp68C+ydcnApNVdXUyyE0GTirPyexcnoNF5CRgGFAAjFTVW7fmr9Wgru7SrGGGVoPNrrfpV4uMtr52bde7koZGq8VG17v3mhVW/Ny18vEPmxhtD1YZbbfVG4wGwA+sNGvng12rYGfarb+e7ThBnf+mPA1gJ+fbsLT23q63Nl8arXD1WqN9Ubire/xqCo3W5JvFrleWOue1715ZzsteROG6Na53cd3GRtuL5Uar8dkW9/h1u9trW0t917sr6422k/NzBKjNV1b7cpPrXbJLI6PVw/4cAHZbk/m7t3AFrPxcxTXnyP4imuU32vAJvAcZFzdCVUeUeL4PUPKPeTHQfitveS7wzFaO3SfHU3PZ7gEvmaoOp8RUVUQmlpyqlmaXZg3pmLo5Q9sLZwAC7p9zsdHeaLW/7+UCo+3HR673mgl/tuIU10q/u680Wl8eMNrRY990j9cTrda88AnXW9P5o04tbet6v65jtU21/RGvzoJvjPbHVue43oOwP7peD/3LaK/0buUe/wi/Mtqd6y9zvbUHW23QEH+JxvtZ9p76pOu9vNNvjXYJdxltn3Gr3eOn9WhttKc41fW2I2W0WviDmPe9bfXux6732tb9jHYy/+N6j56Y+btXdKlrKxMbwPmL8rkBvlLVovJ/KohIb9K3rx0q4v08ynNLu82pahAEVQ8hPRPK5ZEDS4CSt0r7JlrmZ4ocD1wHnKaqG8tybFkoz4CX03RTRM4XkZSIpDZ9uq4cHxcEQWWwE7BLjo8ceANoKSLNRaQm0BOYWNIgIocC95Me7Ere8j0LdBaR3UVkd6Bzom035VrDy4Xkfn4EQP2iFhESDoIdHAFqVNB7qepmEelPeqAqAEap6nsiMhhIqepE4HZgN2CciAD8V1VPU9XVInIj6UETYLCq+msRObLd21JE5KfADarplSoRuQZAVf+U7ZjD2om+/FrmpHJlrT1c7/scZLRbMVt4ABjG743WZupc13thpzuMtgcrXe+NI51Lcb7do67s5R7fGLsyfxeXuN4jec1o/Rjpeldiv2dtpvvXO7a9XWXo9a5dlwOgnpUuaXJLzud10Or5Rnui8BTX6wV/Ok38X9c757SmRruZ61zvIbxjtCvn3G2NTkAJ4KrGNxhtyDSrAdDSSlMb/9S1TnfW6bMF1n639l6jDahnf28BVrBnxvMpRdfzWWp+uYIWTUXU/0uzXAgzKmoNrzIozwzv26kq6fvqnsCZFXJWQRB8b1TkDG9HY7sHvGxT1Qo7syAIvheKgxbVkXJdl6o+DTxdQecSBMEOQMzwgiDIG4qjtNWRGPCCIMggZngVxFtftGO31zN3pz/ToaPrPen0aUYreNxPB/ot9xltVKdfu97jed5oXsQQ8Pd7O8HQtsx0D5/tRJoforfrvYcLczoe4IRxrxhtag8/OrgCJ13Lz+bj7SY27NiMhUZ7jSPd4w9eYaO0TQptiiDAYK43Wqfdful6W53uZCQ87lrdbI+Lmt9jtDW1/HSxIakbjLakg02ZA5hAV6OdveUfrrdxgY3Yt5rjZ1q83cr+HD7gANfbhMzvbwH+30hZqa4zoep6XUEQbCcxwwuCIG+IKG0QBHlDBC2CIMgbqvMtbaVWPC7aWzRVas3+vaEtXK+XpiQ2ywmAYUPPN9rVa/0Mty9XO6lstfz3va1xf6MVMcNoneb6KVHeDsVlA5z8LfzF5iu4zfXOwGbyzJp5uH8Oo620fohfM+K0WjblbOaWQ422auy+Rku/sZXaXfCyaz2EWUYb/aoN3AAwztH8bD6Oaf+c0W7nCntea992j/+wnv193A2/6MWTdDdazSzloc7ZaMuK1UnZ0l0AX9sKVTxWzw/onPV2qXJjZxah76XKlVp2sIiOzdHbJo9Sy4IgqIZU5xleDHhBEGQQQYsgCPKGmOEFQZA3CBGlDYIgTxCgRq4jQ5asnR2Vyo3S7iOaKtUdZP71fgetFu8uM9qy1n6Ec5VTEDNbtGwgg4w2iIGu12usswjbycxN3wK6LXjGill+Qe5saZvPbMDvDnYmDxstWxraiWunGm1Vvdy/j/OwjZP2Z557/H7rbWT9/jo2gg4wIDXCaBcUOQ2WgLa8ZbTu+E181rGb0VqcZX+XPh/t37TNK7DXmy2tayS22c7t2MZPAFsoMNpPJrzreqd1tRH3bOl81/TP/J4VPQ6pFeXrWnboTqLTsuxcKE29ryJKGwRBFUYEatjxuVoQjbiDIMhABHbeObdHbu8nJ4nIByIyT0RM9XgROVZE3hSRzSLSvdRrW0RkZvKYWPrYshIzvCAIMhCBGjne0m77vXLqX/1foC9wufMWX6qq36B5O4gBLwiCTCp2I963/asBRKS4f/W3A56qLkxe81NPKpByXZaILATWAVuAzdtavPy6UQHLrs9cVG6xwC4oA7zX2qb4eDXUAOo6qT+/YILr9Tpu1WeN62346hdGa7rgU6Mt6+0HASY3P9poJ7xra9mBHxy4f7nf/73+XvZ8+w/+u+v96jLn+PVrXe/MOvY/0j9wk9HuxG9vv6iODeh43dgANw3tWm52rU2X2u/5ksZ+jbr7+Y3RLnnwLv8cHA57+n2jbTq5puvd4vz5LKSZ6+32GyeAleWvxQu8LKWR6605KPNnufkVr4hjGanYAc/rX21buGWntoikSIf7blXVf5bnZCrisn6mqn6fwyAIqia5jwwNkgGpmBFJL+qKoqmqLhGRFsBUEZmlqh9t75vFLW0QBJkIOLtosrFyG3d2SyBjL9e+iZYTqrok+Xe+iLwIHAps94BX3iitAs+JyAwRcTdcicj5IpISkdSqT7/zW/QgCMpL8S1tLo9t823/ahGpSbp/dU7RVhHZXURqJV83AI6ixNrf9lDeGd7RyXRzT2CyiMxR1ZdKGpLp7QiANkU7V94u5yAItg8ha8m0spKtf7WIDAZSqjpRRH4CTAB2B04VkUGqejDwI+D+JJixE+k1vO9vwCsx3VwhIhNIR2Re2vpRQRDs0FRwuRSvf7WqXl/i6zdI3+qWPu414JCKO5NypJaJSB1gJ1Vdl3w9GRisqv/Odsy+RXvpRakzM7RN+BGw/Zzb9N4L/HSiUc1tNcgPOND1no3tKnXwUpsSBeBkljG1ue0O1undLAVAbaMqzug82rV66W2thvtdrY65yBa5vIk/uN4OU/9jtPc6+UVXJ3Ka0Rpg41HZIq/e8dfM8dPFSFnpxN5+AO5CbNexLrfZ7wHA+Cv/n/NRdonpSF51j19HXaP1GmsLowJupLltP/93YeYC+3uzrLkf3V9KY6MddouNHgOMvzbzeq8sepV5qbXlSi0r2k00leMwI6/nT2rZXsAEESl+n0e2NtgFQVCFqKapZds94CUbCdtU4LkEQbAjUI0rgFbTywqCYLuJAS8IgryhAqO0OxqVOuCtoy4v0jFD25xlsWC5U2Ou93o/aHEqTxnNS9UCOHi0DVAs6eunKQ1gmNGeGHiWNfZwD6dJ5w+NtuhVv7Yazvr13Red61pfntnZaGPbdnG9kzs56W0r/PS21/a0AYofOduevOAEwBk8asVprpXJF9jzuhvbJQ7gaU422rIr/QX/bqNtClfdvjb1cHqW7KYTedaeVy//59AkI2MqjRecAJjf3NZ9fJCzXa8XFJp2rd+V7pdvl2qN92UFxA9ihhcEQd4QA14QBHlD2VLLqhQx4AVBkEnM8IIgyBsiaBEEQd4QM7yKoS7r6MiLGdrPeMH1/oXfG21O66au1+vuNYN2/knUsdI+g1e71o3X2//mJg+y0cX2W6a7x1/FECvarCHAT4+7YK1f1BOn2dXTbX/uWr0OY+v2tOlTAO2YYbRmLDCa160LoMXlTjHX1q7VTZ86hFmud5aTTvkyx7jeEX3PM5pX4PUMHnOPbzlssdH2GJB7uce+zW0aHMAqGhjtqdNP99/EKWh+27V+BLtPm3sznk/axRZLLTMx4AVBkFdU05Ghml5WEATbTURpgyDIG+KWNgiCvCGitBVDQz7loi3DM7RUgZ8K8wE2BWu/tX59uBpXWa3jfS+63o97NDRa0xX+Qu+ubDDaA5xjtBNe8VO1+u9hgw7DWruV8N0ua6Pr9Xa957V9yGi/5y+u16sL2Ja3XG+TtTboUONz6+vaxK9bd/nQG432PMe73pkDnRQs2/QMgHb9bPG8hTR3vbuu/8poa+rUN1rRaifyA+kau6XwfuYA3XnCfj5fut7Rc50Ahc0QBOCMfqON1jFLcO92rsx4PtP5nS0zMcMLgiBviAEvCIK8oRoPeOXtWhYEQXWkIMdHDojISSLygYjME5GrndePFZE3RWSziHQv9VofEZmbPPqU65qotuN4EATbTQXO8ESkABgOnAAsBt4QkYmluo/9F+gLXF7q2EJgIFBEuiXsjOTYz7b3fLZ5WSIyCjgFWKGqrUucyGNAM2AhcHouJ7Gcvbi9IHMBeGe2uN6DsE1LaszM8saXWWkNdqEa0tkepZm7p2mYBEBNNhntH6tt0OFrZ2c8wBP1bI26AROzNGV3Fssfbvkr1/pwa6u/sNo2rwHYUGizUFrc4mREgFvX7+6WthZce/zMkqGL/mi0SU38xfZJgzoZ7ZT7p7per7HOaxzpeofUudJofXnAaDLOPdzl6o1Oxgxway0bLZvNQa53VEubSfPr1WNd72NX9TXahUPucL0XHjI6U5hXQfXwKi5KezgwL2kJgYg8CnShRH9ZVV2YvFa6cfWJwGRVXZ28Phk4CfC/cTmQyy3t6ORDSnI1MEVVWwJTkudBEFQHytaIu4GIpEo8Ss8I9oGMSqmLEy0XynOsyzZneKr6kog0KyV3gW9LF48BXgSczSFBEFQ5ynZLu7IqtWnc3qDFXqr6SfL1MnDqsSeIyPnFo//6Tytgj1AQBN8txallFRO0WELmDst9E+27Ptal3FFaTXfyztrNW1VHqGqRqhbVaWjXk4Ig2MEo2y3ttngDaCkizUWkJtATmJjjmTwLdBaR3UVkd9JbtW3TkTKwvQPechFpBJD8u6I8JxEEwQ6EALVzfGwDVd0M9Cc9UL0PPK6q74nIYBE5DUBEfiIii0mHze4XkfeSY1cDN5IeNN8ABhcHMLb70tITtG2Y0mt4k0pEaW8HVqnqrcm+mkJVteGx0u9zQJHy18w0oZkn+l28GvOJ0fowxvUOYqDRDlnvpw7VtuXdmNba7wjl1X178DkbpX2o8y/d44c4y5qzRvqfNb+f7Wr1NH6Nu/632JS1f13r5ylN4BdG+xWPuN5FTm7Xr6fbgNjnRTXc4700wU5z/9f1fuXUBXy2jn8NP9tiI70/mPW1653T1tZM9KK8zVjoHr+KPYyWLT3uQWwHu5v4g+s94UybfvjvRzq43s+cHQa9FvzL9b7ZPDO837toAbNTX4przpGi/UVTt+fmlW7MqEpreLlsSxlLOkDRIBmFBwK3Ao+LyLnAx0CWSoZBEFQ5qnGmRS5RWruBKM1xFXwuQRDsKOTrgBcEQZ4RBUCDIMgbqvEtbU5Bi4qiaDfRVKk0rPGv+ClRjVlqtAP4wPUewetG+3BBlnwvJ2Xm8738RfiVBXYB2wsk9B/rN9uZ1ssGKHpmyYpJ8ROj7dPND0gtG1/PaN5iO0CzjQuN9oM1TpE74Pd72Zp613KL0TZR0z1+F6cWW+EKW58OYPWeNsTXgZdcr9dw5w9n+qlWtLTSLwbZIM1w/KY4XoCiScZm//+j03M2IHNTZyfPEbhk411GqzO3dCZVmtWt7femP8MdJ4y949eZwp+L0EWp8gUtWommsvSPKo0cXc2CFkEQ5BnVeIZXTS8rCILtJga8IAjyhhjwgiDIKyJKGwRBXlCNZ3iVG6WtKZoq1TRs/Xw/nfeFWh2N9glOPhLQ2cknbjrV70Q2t5Mt9tly0WLX+0aT1kZru9amrG2q7V/DebX+ZrTjed71ekU191s/3/XWtk28GN/Bj3bfxSVGuzBLxK/napu+JKus76aWfiTS65D24Lt+lzZedTT77QZg0FE2azFbutfLU2162vxONm0vW9Tz6RVOmuB//fN6r6iF0bIVJvXSxS7aeI/r3XW9jd7KZv8c6Jb5tGgmpL7Q8kVpDxFN5ZjeLy0iShsEQVWmGs/wqullBUGw3cSAFwRB3hADXhAE+YRGlLb8bPoaFpfKGHu31jGud38+MtqJa/2uVl53sKaz/Pph2GZZHNlkims9kteMtms9mz71Doe4x//J6W3U9BY/mPL5VTa97Yo6fresWh1sN7Wh02zHMIBuK54x2ns97GI7wAuFPzXayEJbE/AQZrnHP3ixE6DIsth++d03Gs2rxwcwjAFG8zrKQZag1DQblBrV4ddGA9wAxXFFT7nWFU5ng1lT/XqHX7ez2oB6fnrcsK9sUKhGlu/juFdOyXj+WdHLvrEM6E6wKYfinlWRmOEFQZCBCmwuyLUYup8PvKMSA14QBBmoCFt2znVo8GfaOyrlbuITBEH1Y0tBQU6PXBCRk0TkAxGZl7SEKP16LRF5LHl9enFbWBFpJiJfisjM5HFfea8rZnhBEGSgCFsqKLdMRAqA4cAJpBtpvyEiE1V1dgnbucBnqrq/iPQEhgBnJK99pKpZar2VnW3O8ERklIisEJF3S2g3iMiSEiPvyRV1QkEQfL8owmYKcnrkwOHAPFWdr6qbgEeB0lHGLvBth64ngONEpFzZItnYZmqZiBwLfAH8o0TXshuAL1R1aFk+rH5RC+2YujlDW8PurrcZtr3Y6LEX+m/sBNHefsSpBAkMwBa5LMgSShxHD6MVPucUtGzunxZO/c7J7Y92rRuwPXuzndeNXG+06dM6ut5lHWyx0POwKW8AT422vZhe6XuY0dY4aVIAzZ1OYPdxgev961Tb0S3b38+SDoVGm8iprtfrNJdabXcCZPsv+qHXbWrZCvZ0vV5B2kf4leu9k0uNVnfjOtdb51UnEOCkEwImHa9oAKTmli+17MdFNXRSyn7PPZrKio+BlSWkEao6oviJiHQHTlLVfsnzs4D2qtq/hOfdxLM4ef4R0B7YDXgP+BD4HPiDqpYrDJ1LE5+Xiu+pgyCo/pTxlnbld5hL+wnwQ1VdJSLtgH+KyMGq6pfszoHyBC36i8g7yS2vP00LgqBKsoWCnB45sAQyNljum2iuR0R2BuqR7nu9UVVXAajqDOAjwG9knSPbO+DdC+wHtCU9CmdpMAAicr6IpEQktelTfwofBMGOQwWv4b0BtBSR5iJSE+gJlK7FMhHok3zdHZiqqioiDZOgByLSgnTHEr+EUI5sV5RWVZcXfy0ifwMmbcU7AhgB6TW87fm8IAgqj/QtbcVs4FDVzSLSH3iW9CrtKFV9T0QGAylVnQj8HXhQROaRXvnumRx+LDBYRL4mvcP5N6rqd7bKkZzq4SVreJNKBC0aqeonydeXkF6E7Jn9HdIU7S2a6l1K9NeD4TQrvdLKLqADzMAuIRSwxfWuo67R9mee6+0x1o7jw3rZ9KlsXa26jbRpXff26+M4YTYHGc1d2AdYbqUmvT50rR150WgPnuPXqHvjAVuQzgumHMRso4EfZNl9td+1bHjhuUbr/67fKuvT1rsZLVsgwQvILHdSwD5akKX43vtWevtkPwDWZupco33V3n/bH9b52GgrZjZ1ve+1tal/2X6fH5XM9x0BLNXyBS0OKtpFH05li8Rlcpi8X73q4YnIWKAj0EBEFgMDgY4i0hZQYCFkCcUFQVDlUMj1drXKkUuUtpcj59i1MgiCqkfF3dLuaFTPqwqCYLupyEyLHY0Y8IIgMMSAFwRBXhAzvIpiI2YXzZihNp0JYAO7GK09/3G93XnCaNnSn56gu9F6LMqyq2YPK53M00a7mWvdw1/rZztYncMDrtd73zmd/CieF2letDTLfkwn7W7OA/77XsqdRjvVbJmCDhP8n8PfupYOwcOPCv2Irvczm9PaP689vbB0Fl4beZzRJvWzVV/bN3/RPX5E8/OM5n2/AbyU9o/q+MVV+zHSilk2WCygmdGyFb8tnUg33n/LMqEIG6lVAe+04xEzvCAIMogZXhAEeUMMeEEQ5BV5uw8vCIL8oiJTy3Y0KvWqZjRuh9yYWdirDa+7Xq9T1WEjnbwfgB9aqUfnca71bvobbU4Tf7H82SYnGm3ALSOMNrqzX6dvWJFN4VpJA9f7AOcYbdDGga63Tj9bL639gy+63jMueMxol0641/X+z2k/N9rzBccbbVlXW2MvG79mlKu/iq0L2GqFTb8C8H5F/nOa3+HsH/3ONtqQgTcYbdUg/+cw2vk59M0SaHqg0O7Jz3Yr2J7pVvQz1tzUvRpzfG/puIefgFY24pY2CIK8IR2lrfl9n8Z3Qgx4QRBkELe0QRDkFXFLGwRBXhBreEEQ5A0x4FUUGzDdl8452I+AdZhr05e+8gpVARvq1DZaM6eDFsDL2A5WF6y3kVeAAdOt/vlVNYyWKvDrH56zxV6bF/UEONXJAdtYy184rnOHLap5Ife43nuwEeQmXf2CpXuwymjz2M9oXmFTgPNSDxmt3X0zXO9j37YdLUGWYrCHnDbLaI1Y6nrdaKhTlLPP6Mfd45/ua9uZ1WKTf2IOL/AzV7+Om63onwJ1L/vCaO3a+826hmnmboYtRVnCuWUgUsuCIMgbqvMMrzxdy4IgqKZUYNcyROQkEflAROaJyNXO67VE5LHk9ekl28KKyDWJ/oGI2I2xZSRmeEEQZFDctawiSLqODQdOABYDb4jIRFUtubv6XOAzVd1fRHoCQ4AzROQg0g19DgYaA8+LyAGqut37q2OGFwRBBsX78HJ55MDhwDxVna+qm4BHgS6lPF2AMcnXTwDHiYgk+qNJf9oFwLzk/babXJr4NAH+AexFur/HCFUdJiKFwGNAM9KNfE5X1c+29l7tas4g1SyzodIb+N2jPm1pO1U1nGYXcwFqt7CL+I+ssl2xAOa33dtos+v8yPV+0MnWmPMW23fOktDzxB1nGa3bf/0F/7/dbWvJZavpt2VP+79vtk5iv+OvRuuxwK//d3dz+z3zrvfKp+92j2ewlQYyyLV65zvkoRv893XK4b19mZ+X9Rq2BuGAk4cZbdE4v37gYwP7Gm3uoH1d70tOAOxIXnO9B8+x7VTX9/fnGw2nOr/ntqQfAMc0LBUUWlMxDcTKsIbXQERKhiJHJK1Zi9kHMtr6LcaGkb71JG0d15KuRrkPmYmFixNtu8lliN4MXKaqb4pIXWCGiEwG+gJTVPXW5L78aiBLX8EgCKoKirAp99SylVWpTeM2b2lV9RNVfTP5eh3pzp37kDkNHQP84rs6ySAIKo/iNbxcHjmwBChZ6WHfRHM9IrIzUA9YleOxZaJMa3hJ9ORQYDqwV3EzbmAZON2O08ecLyIpEUl9urYcZxoEQaVQwWt4bwAtRaS5iNQkHYQo3TdgIlDcob47MFVVNdF7JlHc5qTry/j9BXIk5yitiOwGPAlcrKqfp9cU06iqioh6xyX38yMAig70PUEQ7FhU1D68ZE2uP/AsUACMUtX3RGQwkFLViaT7XD8oIvNIV7zqmRz7nog8DswmvbR2UXkitACSHki3YRKpAUwCnlXVOxPtA6Cjqn4iIo2AF1X1wK29T1EN0VRhKdHZgA7wdj+7KF2fNa73t9xnNLdpCtDtIRs0uKT3La73AD4wWmM+MVqXOc+5x89vZQMkD2LrtYGfIXDS6Gmu16sPd+p9/rb9hU5DmFH82vXuygajpbDLM30eypIi4KzkjGt1imvdwK5G68gLrncWhxitLn4Aq8NSOwH4d+MORjtyix9cmF5g0zKyNRFqM3au0Zb18msFduFfRruX37rePVhptOfwt6CdNyczu6WoO6TeVXHNObJnURPtkbokJ+89ctmMarWGl4SH/w68XzzYJZSchvYB5ycaBEGVo4LX8HYocrmlPQo4C5glIjMT7VrgVuBxETkX+Bjw+y0GQVClSEdp8zSXVlVfAbJNkW0T0CAIqjTVOZc2UsuCIDDEgBcEQV5Qkbm0OxqVOuB90WZXXkm1ytA+Yn/XeyLPGi3b/zpeh7OWMxe73vG9/5/RvDpwAIucdmjrqGuNH/pR2hbvLzPa5q7+NRzu1XFr61p5qO8vjXY5Q11vhzk2avl6qzaudymNjZainT2t3m+5x7e53EYtrxh6u+utyUaj9ZnrR3+b3jLVaNMe8FMq+zW2qXQjJ/7OaOtP9ON1TxfYeng1s9TDa9PhBqNlSwf0Us4OW+R34ZvbxKay9Vttaw0CcFGp5//1bWUheloEQZA3lDG1rEoRA14QBBnELW0QBHlF3NIGQZAXxLaUCmK3ORs4+og3M7SjL3vT9b7Zw9aoez7Ltr+f8aLR/tS2u+v16rD9YcEdrrdvc9sYp5az2I5fWg0v+ylbQ5hZ/NhoHer5edK9Fz1pxdLp2Al9L7LX0JV/ut4uo23wZUBfW0suW6DpoaG/MtrCV/1ag3OPcmrMZcuStLEUVrGHa/UCUGee9nejuc1+gLueu9aK2fbg1rHSrMY2DQ78NMXbmvR3vd7s6qJ6w13vhimZKXpfF/kpd2UhBrwgCPKKWMMLgiAv+Iad8je1LOhnaMMAABhtSURBVAiC/CNuaYMgyAtiDS8IgrxBiTW8imFvbJsfJwIHfjHKkfRzvT/naaNdM+fPrnd1q9pGG9/cppuBX0S0LuuM9nYrv4NWm5k21eopTnW93vUe2NxG9gD2HufUys/yk7zCSTk7j7+53iZ9Fxlt9kCbWnbAoJlGA7iF64w2dPMfXe8HTmi75Rd+OqD3O9IEe67gR2+9TmZ3cql7/L8722Khe7LC9R52nU0NW17kdjpwMxeydbs7kA+N9oPrv3a9n9+S+btQMSXFI7UsCII8oTrf0kYj7iAIMlCEjdTM6VEeRKRQRCaLyNzk392z+Poknrki0qeE/qKIfCAiM5PHntv6zBjwgiDIoIK7lm2Nq0n3tm4JTEmeZyAihcBA0s27DwcGlhoYf6WqbZOHv/ZQghjwgiAwbKEgp0c5yaW39YnAZFVdraqfAZOBk7b3A7c5RItIE+AfpPvOKjBCVYeJyA3AecCnifVaVbXRgxJ8VL853bsOztA6OmlhAMfzvNEGMdD13mr/Y+DB9893vYWbvzLautZOjTvgxI22Jt/QWpcbzUtXA9x6dtNTHV3r+CIbOHmEM13vL3s8YTSvs1e2c3ttQpbK/M75egGKh+ntHt6MhVb01/A5ZZqtcXdhBz/F7/dFfzHaT8a963r36zHPaLs5gSavixj4ncQOm+nXrbv35j5G6/+QTWMDvzNeI5a63i7dbIrfoPFXut4bvhmSKYwufwOxMq7hNRCRVInnI5LWrLmQS2/rfSAjQrU40Yp5QES2kG4he5Nuow1jLnPSzcBlqvqmiNQFZojI5OS1u1TVrzwZBEGVRBG2fJPzgLdya20aReR50vszSpMR0t9ab+ut8CtVXZKMS0+Sbjb2j60dkEsTn08g3YxVVdeJyPtkjrBBEFQj9Bth41cVk1qmqsdne01ElotIoxK9rb01uCVAxxLP94X0baGqLkn+XScij5Be49vqgFemNTwRaQYcCt+WmugvIu+IyKitRFjOF5GUiKQ2fvp5WT4uCILvAVVhy+aCnB7lJJfe1s8CnUVk92SM6Qw8KyI7i0gDABGpAZwC+OscJch5wBOR3UhPGy9W1c+Be4H9SK/8fAK4CzCqOkJVi1S1qFbDH+T6cUEQfF8olTXg3QqcICJzgeOT54hIkYiMBFDV1cCNwBvJY3Ci1SI98L0DzCQ9E/R31Jcgp7hyMoI+CTysquOTE1le4vW/AZNyvMggCHZgVIXNX3/3G49VdRVOb2tVTcH/pVWp6ihgVCnPenA6TG2DXKK0AvwdeF9V7yyhNyoRYelKDtPJr6jNbA7K0J5IneWb77TS9Efau9alNDLaoK5+VGs/PjJan4l+tywWWKnJAJvS9HCWaOo/X3X0cf5HLSpqYrTp+NfrpbedN9XvanVnJxt13NjVX5+5pr9Nx/vwKBu6PaaX36XtDB4zWv/xftRy/LU2Kj0zS5u2FHZNvFWHj11v4TQbhe+1h71T6lXLj9Ie1/IpK2bpHjfldJsm+O/HbWoa+D/LMdgoLwCXWSlbwdKWf3o74/l/J3/pv2eZEL7ZUj2TsHK5qqNIRz9miUjxHoVrgV4i0pb0VpWFwAXfyRkGQVC5KFD+29UdklyitK8A4ry01T13QRBUUb4R+Cp/Z3hBEOQbm7/vE/huiAEvCIJM0gXxqiWVOuAJSs1SXbseKvql693rkeVGOxGb6gX+In636c/4J2FjA5x8mtMFDLibi4z26zvG2re8zK/N1vEoew43HWVrxoG/KD1gtJ+hM7+v3bh+UydnpRu4ZONdRtt1/Teu98K77c6iYWvt+748obN7/CtdDzPaVwNcK215y2jdsSlzAENMEUUo2NP/i+w13AlGON3FrrryBvf4qxlitBOefsX1Hve4DXD0dAI3AK8956TzZantMemoTkb7JEvhyA9HZ0ZUimzTtrITA14QBHmDAn690SpPDHhBEGSi4LVfrg7EgBcEQSZxSxsEQd4QA14QBHlDDHgVw5fL6vD2kCMytPpXfeZ6T3jVRsaWHFXoer3Uo6HtbaFOgNfusNGyMZf5KT4vc4zRGlxsC5POYz/3+MlrbfpUjSxFqK9qeYPRpvf1U8sGTLPR21M7THS9V9S63Whn1/Ir6HgR2RrrrW9U117u8V7RyKNXvOl6W0xcZrS+Ax5wvQ9wjtH+kyXtrld7J0prG5m5BWYBTpjpRGT9pnRuifNnOdH1HtDZdiLrcN1/XO/Pf2iLo55X+FfXe16HUimFu7m2shEDXhAEeUUMeEEQ5AXfALYGQ7UgBrwgCDKJW9ogCPKGGPAqhnafzSD1RKnCK6f73iZH2UXeaRzrem+Zc6PR5rbyF8A5wko9seliAD9mltHqF6wxWnv8xeeF9fY12lP1bA01gKvSxV4zWJGt5derVmozfK5r3fXxDUZ7Hr/NwOZ6NujwSL1fGe332C5iAK262Rp1d48/1/X+ZcDvjPbhNL/w3KwfHe7qHl855eiG1rHBmOtW+x3Szmxr6/cdw0uu98XRNij1x77XuN4O053fET8GxxWF9ve5dEpmMX9qfnHG82W1HvHftCzEgBcEQV5RTQe8aMQdBEEmxTO8XB7lQEQKRWSyiMxN/s3WCOzfIrJGRCaV0puLyHQRmScij4lIzW19Zgx4QRBk8g3wZY6P8nE1MEVVWwJTkucet5Ouul6aIaR7Y+8PfAb46ycliAEvCIJMFNiS46N8dAHGJF+PAX7hno7qFMisAZf02ukE39YUy3p8SXJp4lMbeIl0W7SdgSdUdaCINAceJb2PfQZwlqr6K6sJ635Uh6lv/DhDuz9LK4wLucdo3u58wG2285dWdlEc4OdH2cr0U7r5gYTLx9vF40echj0jh/ufNeyi8432Aj9zvb90asHN4hDX26qf08Cmv2t1a8ytob7r3R0vIGPr9G3Ev3OYP97W6euf8pv4XNDS6qd28JspHcgHRhv63B9db23nNqvgZPuXeWjh/7rH38mlRmvm/YIB4/vaoEVzFrreZe3rGW3vcWtd7znYgNs66rreI/bJbOLz5KeurezkfrvaQERSJZ6PUFW/kKNlrxKNwJZBtiidyx7AGlUtPtPFwD7bOiiXoMVGoJOqfpG0a3xFRJ4BLiU9nXxURO4jPZ28twwnHATBjkjZorQrVdXmdiaIyPOA/Z8QMirhqqqKiOb8qdtJLk18FPgieVojeSjp6WTxdGcMcAMx4AVB1acCt6Woqr8HChCR5cXtXkWkEZAl09xlFVBfRHZOZnn7km7GvVVyWsMTkYKkReMKYDLwETlOJ0XkfBFJiUhqzafVtIxqEFQnilPLcnmUj4nwbXPePoDfLNghmYi9AHQvy/E5DXiqukVV25IeRQ8HWpXhxEaoapGqFtVvWCPXw4Ig+D6phG0pwK3ACSIyFzg+eY6IFInIyGKTiLxMuoX9cSKyWESKS9JcBVwqIvNIr+n5C8YlKNPGY1VdIyIvAD9lO6aTQRBUASop00JVVwGmXpuqpoB+JZ7bOm1pfT7pCVjO5BKlbQh8nQx2uwAnkN7/UjydfJQcp5M78zV7kdmNLFv0yetE9rWTNgQw9LSLjfbXObbTFYCTLcb48TbaBrCQZka7nKHWmGVq79WtO7CDjTgCvMaRRuu11P+WDmtso7+NH1/qek/eaKPSs2sd5HoPXjHfaMv3tIGzIVm2S+3BSqP9dbX/c6gx3GqPXtbT9Xo1/Q7qPMP1zp7ZzmjXTPyz0X5/op8et7xWllZiDt7vx3FZ6uwtdbqOrRnqR8sPXmp/DpMa205mANxc6vkg31Ym8ryJTyNgjIgUkL4FflxVJ4nIbOBREbkJeIscppNBEFQBivfhVUNyidK+Axzq6GWeTgZBUAWI4gFBEOQNSkWkje2QxIAXBEEm1fiWVtLbWSqHotqiqR+WEq/1vZP62kXaUxbY5iYArzdvY7Srnfpy4Nc2u3Han1zv50fbbTRnF4wx2qPrbboZwEd1WhhtAl1d78nY4MJhE993vZ+eZju1POjmVsOldzh7wbM0epl0gf2ej/y/YNm3XM0Q9/gjhr1ttNsG+Dlvx/Ky0TZSy/WWDnQBtEo56XXAIUW27tyrW44y2g9e91flzzhqtNE68oLr/e1M+7uAHzti/sk22WBTlut9gY5GO42nXO8+56zOeF70FKRWqrjmHJH6RcqxqW0bAZ6SGVvLtNjRiBleEASZxBpeEAR5Q55vSwmCIN+opmt4MeAFQZBJtGkMgiBviFvaimFF6wYMS3XL0AqyzJ29wpVzmjd1vfPY32iP4qcp7f2cU3Txv66VBzqcY7SPnM/aUKe2e/wBa22KUP16tsgmwGFLnYhsltjXQCd/6FDe8s1O8PaVPQ9zradMtVHwU6ZYbdzNp7jHH9XTtlO7nsGut91aG9GdUc9G28GPyI4p8tvdjaOH0X6wyv71TjrKT9W6hLuMdsQce64AY9t2Mdr+bT9yvT+56l2jtR3iFyH1otLHM8X1Tn3gwIzn64recX1lohpvS4kZXhAElojSBkGQF8S2lCAI8oYIWgRBkDfEDK9i+JJdTCeuXdjgemuyMScNoPv6J422vE5D1zuss60ld4hXJC/L580a6xSIOcI9nEua32K0U7OkCHGblVb/2Q+GTHeK1PRjpOP033fV0Aau9cJOdxjtHi7z39fh9r2uMNoxTgoZQA0nS3BpV1szDqB70eVGe2Kcn0rHTCtNvfmnRnuQs93D92CV0Ta2esz19rrN1iuceqX9LIAxQ2yQZeYC33thc/tzyFY3slPHzMBHXb/cYtmJAS8IgrygGm9LiUbcQRBkUkmNuEWkUEQmi8jc5N/ds/j+LSJrRGRSKX20iCwQkZnJo+22PjMGvCAIMilew/vum/hcDUxR1ZbAlOS5x+24O0oBuEJV2yYPZ0EjkxjwgiDI5BvSBUBzeZSPLqR7WpP8+wvPpKpTwGlysx1sc8ATkdoi8h8ReVtE3hORQYle5ulkEARVhEq4pQX2UtVPkq+XAbZj1La5WUTeEZG7RMQvMFiCXIIWG4FOqvqFiNQAXhGRZ5LXrlBVmwOWhaYbFzFy7u8ytJNb2ggrwAlMNtrfsBFWgBYTltnPmvup6203yBY2HMAw1ztjkdMdbqyVXu/lp0R9hu1K9QEHuN79/jzPaKUj2sUUYTt2LaWR6z3sSpuy1mX4c663y1KrexHO4VzkHr+KPYzWkRddLz+y0gZ2da1PzLV3M2/2cN4AuKaHLeZ6J5caLVvK20psBLuB040NYP0AO1/odIufLkYPq/dteY9r3dkZSQ4b6xeDvWFa5vMs9UfLTu51gRuISMk/qhGq+m27PhF5HrDVT+G6jI9TVREpazXia0gPlDWBEaT71Po/2IRcmvgo8EXytEbyqLwyyUEQ7Mis3FrFY1U9PttrIrJcRBqp6ici0ghYUZYPLjE73CgiDwB2/1IpclrDE5ECEZmZnNBkVZ2evLTN6aSInC8iKRFJfbracwRBkKdMJN3TGnLsbV2SZJBERIT0+p+t0FCKnAY8Vd2iqm2BfYHDRaQ16elkK+AnQCHp6aR37AhVLVLVooaFOV1HEAT5wa3ACSIyFzg+eY6IFInItzvpReRlYBxwnIgsFpETk5ceFpFZwCygAXDTtj6wTBuPVXWNiLwAnKSqQxM55+lkEARVgeIw7XeLqq4CjnP0FPxf9yhVdRbTQVX9Gl9bYZtdy0SkIfB1MtjtAjwHDAFmJPfeAtwFfKWq2fbRANBUxBh+28H3XvviH43WOMuSbP/U3412VtEIxwkXcJ/Rdi3DD9dL8fFSvQCuPOtuo31lPx6A2XXsIrxbIw9o1NjWXGvPdMcJ/7zK6ahmfsXSdOz8jNF+z1+M1m2C9QG82dVew4Eb/Vyn5bX2NNqXWYIWXqe3bOmAXabawMvX7axvS5b/6rvVsUG0g5jteofOtL+j63/k3zRt2bnAaBsK/Ov16h0W4XcROy/1UKbvbEjNLmfXMjlUYdq2jQDUq3ZdyxoBY0SkgPQt8OOqOklEpiaDoZDOYPzNd3ieQRBUGtW3ekAuUdp3gEMdvczTySAIqgLVN5k2igcEQVCKGPCCIMgblMoIWnwfxIAXBEEp8ngNr0JptwekMjtejWKTa/3QScHKFqlivZWO5SXX6kXcCt/161kPa21T2fpttIU2Z9RywoDAnx682Ggn8z+u97DpNiKrLV0rE+hmtJpZvo9LhtjNj/vc5u8Av72zLeA5m4OMdlBXm9oGsB82evzUUr+7WIvNNh1wWct6rtfr4nUIfneuaZ1sxHy2k8f226VjjAbw9JRfGk2Pdq2Ma2u7t2WLHntpgq9xpOvtyAtG67Uiy57cRaWe+78GZSRuaYMgyBtihhcEQd4QM7wgCPKGmOEFQZA3VE5q2ffBNlPLKpKd27XR3V57NkNbM9ev4+YVirmq0w2udRM1jVafNa534Gm2jdfrE/16dnWdIqsHD59vtEkX+XuwvTS0bAvVv+F++1kP2c8CoImV/tjhGtd6PFOMVrTRD/60qjXHaGPpabTmLHSPr7vFfr9+8L5/a3R363ONli09rpnzeQ0nfmGNwN2n2fe9Zr2tkTe/Tgv3+Ppr7fuOrNfHcUJ3bBpaw9H+eVHHSmf2sCmRAAVOPbyuTHC9pYMh9xeNYWlqWTlTy1oqTkqhz8nVLrUsCIK8I25pgyDICyJoEQRB3hADXhAEeUNEaYMgyBsiSlshFO0qmjowU5vzVlPX+ywnGi1bIcY1TnewHs9NcpykK/eVYvKVfu7QzZmNlQC4Ol2FOoM9s/QeOew5p4BnlnSxe5vbSOBH7O96VzrdwUYPv9B/49Os9K8mnV3rIif828TkLkGX+/2uZx9f0NBoT3Oy6/3NapvaNbKwt+v1upkNONMv8Dr1EdtlzYu2H7Leb3/wkzr/MdowBrjeTjOdDmW1XSu/azXEaCuydCW81yktWTjNT38s3UWv6ElIfVreAqBNNXtP7NJcGFHaIAiqMtX3ljanJj5BEOQTxUGLXB7bj4gUishkEZmb/Lu742krIv8rIu8lHRLPKPFacxGZLiLzROQxEbEbcksRA14QBKUonuHl8igXVwNTVLUlMAX/PnoDcLaqHgycBPxZRIrXsIYAd6nq/sBngN11XooY8IIgKEVx0CKXR7noAhQv5o4h3Vs2A1X9UFXnJl8vJZ2D1TBpHtYJeGJrx5emUoMWIvIp8HHytAGwstI+vPKI66p6VKdra6qqNnpUBkTk36S/J7lQGygZURmhqn5EyX7OGlWtn3wtwGfFz7P4Dyc9sB1Muhf268nsDhFpAjyjqq239pmVGrQo+YMQkVRViu7kSlxX1aM6X9v2oKonVdR7icjzwN7OSxlbIFRVRSTr7EtEGgEPAn1U9Zv0+Fh2IkobBMF3hqoen+01EVkuIo2S/taNcEuGgIj8APgf4DpVfT2RVwH1RWRnVd0M7Ass2db5xBpeEATfFxOB4g2ofQBTxz6JvE4A/qGqxet1aHot7gWg+9aOL833OeDldJ9fBYnrqnpU52vbkbkVOEFE5gLHJ88RkSIRKW4eczpwLNBXRGYmj7bJa1cBl4rIPGAPwK+3VYJKDVoEQRB8n8QtbRAEeUMMeEEQ5A2VPuCJyEki8kGSDpJrhvIOiYiMEpEVIvJuCW2b6TI7OiLSREReEJHZSUrPgESv0tcmIrVF5D8i8nZyXYMSvcwpSkHVpFIHPBEpAIYD/w84COglIrbTc9VhNOl0l5Lkki6zo7MZuExVDwKOAC5Kfk5V/do2Ap1UtQ3QFjhJRI5gO1KUgqpJZc/wDgfmqep8Vd0EPEo6vaRKoqovAatLydtMl9nRUdVPVPXN5Ot1wPvAPlTxa9M0xV12aiQPZTtSlIKqSWUPePtARoG1xYlWndhLVT9Jvl4GWYqeVRFEpBlwKDCdanBtIlIgIjNJb3KdDHwErEk2r0L1/J0MEiJo8R2SbI6ssvt+RGQ34EngYlX9vORrVfXaVHWLqrYlvTP/cKDV93xKQSVS2QPeEjK7quaUDlLFWJ6kyRTn//nlkHdwRKQG6cHuYVUdn8jV4toAVHUN6Z36PyVJUUpeqo6/k0FCZQ94bwAtk6hYTaAn6fSS6sQ202V2dJLKFX8H3lfVO0u8VKWvTUQaFtdSE5FdgBNIr0+WOUUpqJpUeqaFiJwM/BkoAEap6s2VegIViIiMBTqSLqWzHBgI/BN4HPgh6VJYp6tq6cDGDo2IHA28DMwiXRwN4FrS63hV9tpE5MekgxIFpP+zf1xVB4tIC9IBtELgLaC3qm78/s40+K6I1LIgCPKGCFoEQZA3xIAXBEHeEANeEAR5Qwx4QRDkDTHgBUGQN8SAFwRB3hADXhAEecP/B8eDMyB22tthAAAAAElFTkSuQmCC\n",
"text/plain": [
"<Figure size 432x288 with 2 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"plt.imshow(y, cmap=\"jet\", vmax=.2, vmin=-.2)\n",
"plt.colorbar()"
]
},
{
"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.5"
}
},
"nbformat": 4,
"nbformat_minor": 4
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment