{ "nbformat": 4, "nbformat_minor": 0, "metadata": { "colab": { "provenance": [] }, "kernelspec": { "name": "python3", "display_name": "Python 3" }, "language_info": { "name": "python" } }, "cells": [ { "cell_type": "markdown", "source": [ "###Super simple Ridge Regression###" ], "metadata": { "id": "389Ux-WmmlYW" } }, { "cell_type": "code", "execution_count": 363, "metadata": { "colab": { "base_uri": "https://localhost:8080/" }, "id": "p_ziA-X42KY7", "outputId": "c62b4374-3d4d-4e17-945f-6868582adc1b" }, "outputs": [ { "output_type": "stream", "name": "stdout", "text": [ " Size in feet^2 (x) Price ($)(y)\n", "0 2104 399900\n", "1 1600 329900\n", "2 2400 369000\n", "3 1416 232000\n", "4 3000 539900\n" ] } ], "source": [ "\n", "#Importing stuff\n", "from sklearn.linear_model import Ridge\n", "import numpy as np\n", "import matplotlib.pyplot as plt\n", "import pandas as pd\n", "\n", "# Load dataset downloaded from https://www.kaggle.com/datasets/kennethjohn/housingprice\n", "# The dataset contains some information about housing in Portland, such as size, number of bedrooms and price\n", "dataset = pd.read_csv('./Housing in portland dataset.csv', header=None)\n", "dataset = dataset[[0,2]]\n", "# For this example i'll only consider size and price columns\n", "size_label = 'Size in feet^2 (x)'\n", "price_label = 'Price ($)(y)'\n", "dataset.columns = ([size_label,price_label])\n", "print(dataset.head())" ] }, { "cell_type": "code", "source": [ "def plot_df(df):\n", " df.plot(\n", " x=size_label,\n", " y=price_label,\n", " title='Ridge Regression',\n", " xlabel = 'Size',\n", " ylabel = 'Price',\n", " style = '.',\n", " label='Price'\n", " )\n", "plot_df(dataset)" ], "metadata": { "colab": { "base_uri": "https://localhost:8080/", "height": 295 }, "id": "p-pVhfR0QTB9", "outputId": "33196d2d-ad04-4648-a473-dd7d5ea0e1dd" }, "execution_count": 364, "outputs": [ { "output_type": "display_data", "data": { "text/plain": [ "
" ], "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZcAAAEWCAYAAACqitpwAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nO3de7hV1Xnv8e9Prt5BpFbdBrASDRolsqukiW0Sq0GTA56eaE1tpZaEtmKbHNvT6Olpzcml1Z7YJrZGa71BSmKISSonjwSJl/bk8YG4SbxfyhblcRMUAniLN7a85485Fix21l577c2c6/r7PM969lxjzTXHOycP+91zjDHHUERgZmaWp30aHYCZmbUfJxczM8udk4uZmeXOycXMzHLn5GJmZrlzcjEzs9w5uVhHknS9pL+q8nlIOqaeMTU7Se+Q9KqkUY2OxZqf/JyLtSNJzwKHAW8DrwLfBy6JiFdr/H4A0yOiN+e47gNmA/3AG8B/AIsiYlOe9Zg1mu9crJ39l4g4AJgJvAe4vMHxlFyS4joGOAD4Ut4VSBqd9zHNhsPJxdpeRDwPrCRLMgBIulXSF8re/w9JmyT9VNIflH9f0iRJ/1fSy5IekPQFST8s+/w4SaskbZP0lKTzaozrReDfBsQ16LFqiCMkLZK0DliXyj4q6UFJL0q6X9KJZft/RtJGSa+kuk5P5adI6kn1vCDp71P51FTH6PT+CEnLU6y9kj5ZduzPSlomaUk6/mOSumu5LtYenFys7UnqAs4CKjZxSZoD/DlwBjAd+M0Bu1wL/Bz4ZWB+epW+uz+wCvg68EvA+cBXJc2oIa5JwG+V4qrhWIPGUeYc4FRghqT3ADcDfwhMAv4ZWC5pnKRjgUuAX42IA4EPA8+mY3wF+EpEHAT8CrBskFO4DegDjgA+BvyNpA+VfT437TMBWA7801DXxNqHk4u1s3+T9ArwHLAZuGKQ/c4DbomIRyPi58BnSx+kzuv/BlwREa9FxOPA4rLvfhR4NiJuiYj+iPgJ8G3g3CpxXSPpJeBnwKHAnwx1rBriKPnbiNgWEa8DC4F/jog1EfF2RCwG3iTr83kbGEeWhMZExLMR8XQ6xg7gGEmHRsSrEbF6YCWSjgLeB3wmIt6IiAeBG4ELy3b7YUTcGRFvA18DTqpyTazNOLlYOzsn/VX+AeA4sl/klRxBloBKNpRtTwZGD/i8fHsKcGpqdnpR0ovABWR3F4P504g4GDgRmAh01XCsoeIYLLY/G3C8o4Aj0kCFT5Ml0s2SbpN0RPreAuCdwJOp+e2jFeo5AtgWEa+UlW0Ajix7/3zZ9mvAePcFdQ4nF2t7EfHvwK0M3nG+ieyXbsk7yra3kI3s6iorK9/3OeDfI2JC2euAiPjjGuJ6BPgCcK0kDXGsoeLYddgBsX1xwPH2i4hvpPq/HhHvJ0tCAVyVytdFxMfJmuauAm5PTXblfgocIunAsrJ3ABuHOm/rDE4u1im+DJwhqVLTzDLg9yXNkLQfZc1nqUnnO8BnJe0n6Tj2bPr5HvBOSb8naUx6/aqkd9UY12KyIdNzqx2rhjgq+RfgjySdqsz+kj4i6UBJx0r6kKRxZEOiXwd2Akj6XUmTI2In8GI61s7yA0fEc8D9wN9KGp8GCiwA/rXG87Y25+RiHSEitgBLgL+u8NkKsuRzD1nn+j0DdrkEOJismedrwDfI+i5IzUJnknW+/zTtcxVZf0Ytcb1F1oH+VzUca9A4Bjl2D/BJso707encfj99PA64kqzf53myu5TSUO05wGOSXk2xnZ/6cAb6ODA1xfpdsv6gH9Ry3tb+/BCl2TBJugr45YioNFqr4+Iwq8R3LmZDSM+enJialk4ha/75bqfGYVYLj9wwG9qBZE1QRwAvAFcDd3RwHGZDcrOYmZnlzs1iZmaWOzeLJYceemhMnTq10WGYmbWUtWvX/iwiJg8sd3JJpk6dSk9PT6PDMDNrKZI2VCp3s5iZmeXOycXMzHLn5GJmZrlzn0sVO3bsoK+vjzfeeKPRoRRu/PjxdHV1MWbMmEaHYmZtwMmlir6+Pg488ECmTp1KNmlte4oItm7dSl9fH9OmTWt0OGbWBgprFkuzrj5Y9npZ0qclHaJsGdd16efEtL8kXZOWS31Y0sllx5qf9l8nqXwVwFmSHknfuSZNW85gdQzXG2+8waRJk9o6sQBIYtKkSR1xh2Zm9VFYcomIpyJiZkTMBGaRLRb0XeAy4O6ImA7cnd5Dtgzt9PRaCFwHWaIgmwL9VOAU4IqyZHEd2ayvpe/NSeWD1TFs7Z5YSjrlPM1st7UbtnPtvb2s3bA992PXq0P/dODpiNgAzGP38qyLydb8JpUvicxqYIKkw8nW9l6Vlm7dTrbG+Jz02UERsTqyOWyWDDhWpTrMzIwssVxw42quvuspLrhxde4Jpl7J5XyyCfcADouITWn7ebKFkiBbHrV8ida+VFatvK9CebU69iBpoaQeST1btmwZ9knVw6hRo5g5cyYnnHAC5557Lq+99lrF/X7t136tzpGZWStbvX4rb/XvZGfAjv6drF6/NdfjF55cJI0lW2XvWwM/S3cchc6cWa2OiLghIrojonvy5F+YvaAp7Lvvvjz44IM8+uijjB07luuvv36Pz/v7+wG4//77GxGembWo2UdPYuzofRglGDN6H2YfPSnX49fjzuUs4McR8UJ6/0Jq0iL93JzKN7LnmuBdqaxaeVeF8mp1FK7INszTTjuN3t5e7rvvPk477TTmzp3LjBkzADjggAN27XfVVVfx7ne/m5NOOonLLsu6m55++mnmzJnDrFmzOO2003jyySdzj8/MWsesKRNZ+onZXHrmsSz9xGxmTRnRuKdB1WMo8sfZ3SQGsByYT7bE6nx2r0exHLhE0m1knfcvRcQmSSuBvynrxD8TuDwitqURaLOBNWTrif/jEHUUqtSG+Vb/TsaO3ifXf7D+/n5WrFjBnDnZmIUf//jHPProo78wdHjFihXccccdrFmzhv32249t27YBsHDhQq6//nqmT5/OmjVruPjii7nnnoGr+ZpZJ5k1ZWLuSaWk0OQiaX/gDOAPy4qvBJZJWgBsAM5L5XcCZ5Ot8/0acBFASiKfBx5I+30uIral7YuBW4F9gRXpVa2OQlVqw9zbf7jXX3+dmTNnAtmdy4IFC7j//vs55ZRTKj6T8oMf/ICLLrqI/fbbD4BDDjmEV199lfvvv59zzz13135vvjno0utmZnut0OQSET8HJg0o20o2emzgvgEsGuQ4NwM3VyjvAU6oUF6xjqKV2jB39O/MrQ2z1Ocy0P7771/zMXbu3MmECRMqHsfMrAieWyxHRbdh1uKMM87glltu2TWqbNu2bRx00EFMmzaNb30rG1MRETz00EN1j83MOoeTS85mTZnIog8e05DEAjBnzhzmzp1Ld3c3M2fO5Etf+hIAS5cu5aabbuKkk07i+OOP5447vPS6mRVHWWuUdXd3x8DFwp544gne9a53NSii+uu08zWzvSdpbUR0Dyz3nYuZmeXOycXMzHLn5DKETmk27JTzNLP6cHKpYvz48WzdurXtf/GW1nMZP358o0MxszbhxcKq6Orqoq+vj2ad1DJPpZUozczy4ORSxZgxY7wyo5nZCLhZzMzMcufkYmZmuXNyMTOz3Dm5mJlZ7pxczMwsd04uZmaWOycXMzPLnZOLmZnlzsnFzMxy5+RiZma5c3IxM7PcObmYmVnunFzMzCx3Ti5mZpY7JxczM8udk4uZmeXOycXMzHLn5GJmZrlzcjEzs9wVmlwkTZB0u6QnJT0h6b2SDpG0StK69HNi2leSrpHUK+lhSSeXHWd+2n+dpPll5bMkPZK+c40kpfKKdZiZWX0UfefyFeD7EXEccBLwBHAZcHdETAfuTu8BzgKmp9dC4DrIEgVwBXAqcApwRVmyuA74ZNn35qTyweowM7M6KCy5SDoY+HXgJoCIeCsiXgTmAYvTbouBc9L2PGBJZFYDEyQdDnwYWBUR2yJiO7AKmJM+OygiVkdEAEsGHKtSHWZWo7UbtnPtvb2s3bC90aFYCxpd4LGnAVuAWySdBKwFPgUcFhGb0j7PA4el7SOB58q+35fKqpX3VSinSh1mVoO1G7ZzwY2reat/J2NH78PST8xm1hS3LlvtimwWGw2cDFwXEe8Bfs6A5ql0xxEFxlC1DkkLJfVI6tmyZUuRYZi1lNXrt/JW/052Buzo38nq9VsbHZK1mCKTSx/QFxFr0vvbyZLNC6lJi/Rzc/p8I3BU2fe7Ulm18q4K5VSpYw8RcUNEdEdE9+TJk0d0kmbtaPbRkxg7eh9GCcaM3ofZR09qdEjWYgpLLhHxPPCcpGNT0enA48ByoDTiaz5wR9peDlyYRo3NBl5KTVsrgTMlTUwd+WcCK9NnL0uanUaJXTjgWJXqMLMazJoykaWfmM2lZx7rJjEbkSL7XAD+BFgqaSywHriILKEtk7QA2ACcl/a9Ezgb6AVeS/sSEdskfR54IO33uYjYlrYvBm4F9gVWpBfAlYPUYWY1mjVlopOKjZiyLgnr7u6Onp6eRodhZtZSJK2NiO6B5X5C38xGzMOVbTBFN4uZWZvycGWrxncuZjYinT5c2Xdt1fnOxcxGpDRceUf/zo4bruy7tqE5uZjZiJSGK69ev5XZR0/qqF+ule7aOun8a+HkYmYj1qnDlTv5rq1WTi5mZsPUyXdttXJyMTMbgU69a6uVR4uZmVnunFzMzCx3Ti5mZpY7JxczM8udk4uZWZNrxdkAPFrMzKyJtepsAL5zMTNrYq06h5uTi5nt0orNL+2uVZecdrOYmQGt2/zS7lp1NgAnFzMDPBljM2vF2QDcLGZmQOs2v1hz8p2LmQGt2/xizcnJxcx2acXmF2tObhYzM7PcObmYmVnunFzMzCx3Ti5mZpY7JxczM8udk4uZmeXOycXMzHLn5GJmZrkrNLlIelbSI5IelNSTyg6RtErSuvRzYiqXpGsk9Up6WNLJZceZn/ZfJ2l+WfmsdPze9F1Vq8PMzOqjHncuH4yImRHRnd5fBtwdEdOBu9N7gLOA6em1ELgOskQBXAGcCpwCXFGWLK4DPln2vTlD1GHW1DzlvbWLRjSLzQMWp+3FwDll5UsisxqYIOlw4MPAqojYFhHbgVXAnPTZQRGxOiICWDLgWJXqMGtapSnvr77rKS64cbUTjLW0opNLAHdJWitpYSo7LCI2pe3ngcPS9pHAc2Xf7Utl1cr7KpRXq2MPkhZK6pHUs2XLlmGfnFmeWnXFQbNKip648v0RsVHSLwGrJD1Z/mFEhKQoMoBqdUTEDcANAN3d3YXGYTaU0pT3O/p3esp7a3mFJpeI2Jh+bpb0XbI+kxckHR4Rm1LT1ua0+0bgqLKvd6WyjcAHBpTfl8q7KuxPlTrMmpanvLd2UlizmKT9JR1Y2gbOBB4FlgOlEV/zgTvS9nLgwjRqbDbwUmraWgmcKWli6sg/E1iZPntZ0uw0SuzCAceqVIdZU5s1ZSKLPniME4u1vCLvXA4DvptGB48Gvh4R35f0ALBM0gJgA3Be2v9O4GygF3gNuAggIrZJ+jzwQNrvcxGxLW1fDNwK7AusSC+AKwepw8zM6kDZQCvr7u6Onp6eRodhZtZSJK0te9RkFz+hb2ZmuXNysY7nBxfN8lf0UGSzplZ6cPGt/p2MHb0PSz8x253pZjnwnYt1tEY8uOg7JesEvnOxjlbvBxd9p2SdwsnFdlm7YXvHPcBX7wcXK90pdcq1ts7i5GJAZ/9FPWvKxLqdq6d4sU7h5GKA/6KuF0/xYp3CycUA/0VdT/W8UzJrlJqSi6R3ki3MdVhEnCDpRGBuRHyh0OisbvwXtZnlqdahyP8CXA7sAIiIh4HziwrKGsOTJg7Ow4fNhqfWZrH9IuJHaRLKkv4C4jFrOp082MFspGq9c/mZpF8hW1kSSR8DNlX/ill78AqRZsNX653LIrIVG4+TtBF4BvjdwqIyayIe7GA2fDUll4hYD/xmWvRrn4h4pdiwzJqHBzuYDV+to8X+Bvi7iHgxvZ8I/FlE/K8igzNrFh4+bDY8tfa5nFVKLAARsZ1s1UgzM7NfUGtyGSVpXOmNpH2BcVX2NzOzDlZrh/5S4G5Jt6T3FwGLiwnJzMxaXa0d+ldJehg4PRV9PiJWFheWmZm1sprnFouIFcCKAmMxM7M2UTW5SPphRLxf0iukByhLHwEREQcVGp2ZmbWkqsklIt6ffh5Yn3DMzKwdDDlaTNIoSU/WIxgzM2sPQyaXiHgbeErSO+oQj5mZtYFaO/QnAo9J+hHw81JhRMwtJCozM2tptSaXvyo0CjMzaytDjRYbD/wRcAzwCHBTRHgdFzMzq2qoPpfFQDdZYjkLuHq4FaQBAT+R9L30fpqkNZJ6JX1T0thUPi69702fTy07xuWp/ClJHy4rn5PKeiVdVlZesQ4zM6uPoZLLjIj43Yj4Z+BjwGkjqONTwBNl768C/iEijgG2AwtS+QJgeyr/h7QfkmaQLal8PDAH+GpKWKOAa8mS3gzg42nfanVYB/ISxWb1N1Ry2VHaGElzmKQu4CPAjem9gA8Bt6ddFgPnpO157J6v7Hbg9LT/POC2iHgzIp4BeoFT0qs3ItZHxFvAbcC8IeqwDlNaovjqu57ightXO8GY1clQyeUkSS+n1yvAiaVtSS/XcPwvA38B7EzvJwEvliWqPuDItH0k8BzsSmQvpf13lQ/4zmDl1erYg6SFknok9WzZsqWG07FW4yWKzRqjanKJiFERcVB6HRgRo8u2q079IumjwOaIWJtrxDmKiBsiojsiuidPntzocKwApSWKRwkvUWxWRzVPXDkC7wPmSjobGA8cBHwFmCBpdLqz6AI2pv03AkcBfZJGAwcDW8vKS8q/U6l8a5U6rMN4iWKzxqh1sbBhi4jLI6IrIqaSdcjfExEXAPeSDQ4AmA/ckbaXp/ekz++JiEjl56fRZNOA6cCPgAeA6Wlk2NhUx/L0ncHqsA40a8pEFn3wGCcWszoqLLlU8RngUkm9ZP0jN6Xym4BJqfxS4DKAiHgMWAY8DnwfWBQRb6e7kkuAlWSj0ZalfavVYWZmdaDsD33r7u6Onp6eRodhZtZSJK2NiO6B5Y24czEzszbn5GJmZrlzcjEzs9w5uZiZWe6cXCxXnTqPV6eet9lginyI0jpMaR6vt/p3Mnb0Piz9xOyOeLakU8/brBrfuVhuOnUer049b7NqnFwsN506j1ennrdZNX6IMvFDlPlYu2F7R87j1annbTbYQ5Tuc7FczZoysSN/uXbqeZsNxs1iZmaWOycXaxoezmvWPtwsZk1hpMN53ddh1pycXKwpVBrOO1SyGG5CciIyqx8nF2sKpeG8O/p31jycdzgJyQ86mtWXk4s1hZEsRzychDSSOyMzGzknF2sawx3OO5yENJI7IzMbOT9EmfghyvbnPhez/PkhSstd0b+s8z5+LXdGTkBm+XBysREpuoO8ER3w7vQ3y48forQRKXom4EbMNOzZjc3y4+RiI1L0TMCNmGnYsxub5ccd+ok79Iev1fpcmrVOs1Y2WIe+k0vi5GJmNnyDJRc3i9mQWmlCyVaKFVovXrNaebSYVdVKI6haKVZovXjNhsN3LlZVK42gaqVYofXiNRsOJxerqpVGULVSrNB68ZoNR2Ed+pLGA/8BjCNrfrs9Iq6QNA24DZgErAV+LyLekjQOWALMArYCvx0Rz6ZjXQ4sAN4G/jQiVqbyOcBXgFHAjRFxZSqvWEe1eN2hP7hWGkHVSrFC68VrNlDdR4tJErB/RLwqaQzwQ+BTwKXAdyLiNknXAw9FxHWSLgZOjIg/knQ+8F8j4rclzQC+AZwCHAH8AHhnquY/gTOAPuAB4OMR8bikZZXqqBZvOyaXVvzFNZyYW/H8atXO52btpe5zi0WWtV5Nb8ekVwAfAn4nlS8GPgtcB8xL2wC3A/+UEtQ84LaIeBN4RlIvWaIB6I2I9QCSbgPmSXqiSh1toZZfPK3YWTycmFvx/GrVzudmnaPQPhdJoyQ9CGwGVgFPAy9GRH/apQ84Mm0fCTwHkD5/iaxZa1f5gO8MVj6pSh0D41soqUdSz5YtW/bmVOum9Ivn6rue4oIbVw86hLUVO4uHE3Mrnl+t2vncrHMUmlwi4u2ImAl0kd1tHFdkfcMVETdERHdEdE+ePLnR4dSk1l88rdhZPJyYW/H8atXO52adoy7PuUTEi5LuBd4LTJA0Ot1ZdAEb024bgaOAPkmjgYPJOvZL5SXl36lUvrVKHS2v1kWvRrKyY6MNJ+ZWPL9atfO5WecoskN/MrAjJZZ9gbuAq4D5wLfLOtsfjoivSloEvLusQ/+3IuI8SccDX2d3h/7dwHRAZB36p5MljweA34mIxyR9q1Id1eJtpQ59d/buPV9Ds3w0YrGww4HFkkaRNb8ti4jvSXocuE3SF4CfADel/W8CvpY67LcB5wOkZLEMeBzoBxZFxNvppC4BVpINRb45Ih5Lx/rMIHW0heEuB9zs6v2L3h3mZsUrcrTYw8B7KpSvZ/dor/LyN4BzBznWF4EvVii/E7iz1jqs+TTiF32lfisnF7N8+Ql9a6hGjIxyh7lZ8TxxpTVUrQMU8uQOc7PieT2XpJU69NuNO9fNWlcjOvTNatJuAxTMzH0uZmZWACcXMzPLnZOLmZnlzsmljXTKeuydcp5mrcwd+m2iU54675TzNGt1vnNpE+04TXulO5R2PE+zduQ7lzbRiIcRR2pvFjtrpfM062ROLk1oJA8V7s1T5/V8iLHWZq3B5v/y0/VmrcHJpcnsTZ/CSB5GrHcfRq2TRla7Q/FDl2bNz8mlydR7xt5619fOi52Z2W5OLk2m3n0Kg9VXa79Ikc13vkMxa12euDJppokri/rFXmt9tTSVNWJIsCe4NGs+nriyhQz1F3vev9gH1ldLU1m9m9P8fItZa/FzLi2o6Gc9allMq94Lbvn5FrPW4juXFpRXv8xgzUy19IvUu8Pdz7eYtRb3uSTN1OdSi/LEAAz7l3wrNjO5z8Ws+bjPpc2U+klGmiTq3WeSB48eM2sd7nNpcSPti6h3n4mZdRbfubS4kfZF+CFFMyuS+1ySVutzKee+CDNrFPe5tDH3RZhZs3GfS4fxKo5mVg++c+kgrTj82Mxak+9cOoifcjezenFy6SAefmxm9VJYcpF0lKR7JT0u6TFJn0rlh0haJWld+jkxlUvSNZJ6JT0s6eSyY81P+6+TNL+sfJakR9J3rpGkanV0utLw40vPPNZNYmZWqCLvXPqBP4uIGcBsYJGkGcBlwN0RMR24O70HOAuYnl4LgesgSxTAFcCpwCnAFWXJ4jrgk2Xfm5PKB6uj7Qy3g37WlIks+uAxuxKLO/jNrAiFdehHxCZgU9p+RdITwJHAPOADabfFwH3AZ1L5ksgevFktaYKkw9O+qyJiG4CkVcAcSfcBB0XE6lS+BDgHWFGljraytx307uA3s6LUpc9F0lTgPcAa4LCUeACeBw5L20cCz5V9rS+VVSvvq1BOlToGxrVQUo+kni1btgz/xBpsbzvo3cFvZkUpPLlIOgD4NvDpiHi5/LN0l1LoFAHV6oiIGyKiOyK6J0+eXGQYhdjbDnp38JtZUQp9zkXSGLLEsjQivpOKX5B0eERsSs1em1P5RuCosq93pbKN7G7iKpXfl8q7KuxfrY62srfzg3l+MTMrSpGjxQTcBDwREX9f9tFyoDTiaz5wR1n5hWnU2GzgpdS0tRI4U9LE1JF/JrAyffaypNmprgsHHKtSHW1nYAd9vb9vZlZJkXcu7wN+D3hE0oOp7H8CVwLLJC0ANgDnpc/uBM4GeoHXgIsAImKbpM8DD6T9Plfq3AcuBm4F9iXryF+Rygerw8zM6sCzIicjnRXZMxKbWSfzrMgF8FBeM7PKPP3LXvBQXjOzypxc9oKH8pqZVeZmsb0wkqG87qMxs07g5LKXhrMKpPtozKxTuFmsjtxHY2adwsmljtxHY2adws1ideTpVsysUzi51Nlw+mjMzFqVm8XMzCx3Ti5mZpY7JxczM8udk4uZmeXOycXMzHLn5GJmZrnzei6JpC1kC4vVw6HAz+pU195wnPlynPlynPkaaZxTImLywEInlwaQ1FNpcZ1m4zjz5Tjz5TjzlXecbhYzM7PcObmYmVnunFwa44ZGB1Ajx5kvx5kvx5mvXON0n4uZmeXOdy5mZpY7JxczM8udk0tOJN0sabOkR8vKDpG0StK69HNiKpekayT1SnpY0sll35mf9l8naX6d4vyspI2SHkyvs8s+uzzF+ZSkD5eVz0llvZIuKyDOoyTdK+lxSY9J+lQqb6prWiXOprqmksZL+pGkh1Kc/zuVT5O0JtX5TUljU/m49L43fT51qPgLjPFWSc+UXcuZqbxh/49SHaMk/UTS99L7prmWQ8RZn+sZEX7l8AJ+HTgZeLSs7O+Ay9L2ZcBVaftsYAUgYDawJpUfAqxPPyem7Yl1iPOzwJ9X2HcG8BAwDpgGPA2MSq+ngaOBsWmfGTnHeThwcto+EPjPFE9TXdMqcTbVNU3X5YC0PQZYk67TMuD8VH498Mdp+2Lg+rR9PvDNavEXHOOtwMcq7N+w/0epnkuBrwPfS++b5loOEWddrqfvXHISEf8BbBtQPA9YnLYXA+eUlS+JzGpggqTDgQ8DqyJiW0RsB1YBc+oQ52DmAbdFxJsR8QzQC5ySXr0RsT4i3gJuS/vmGeemiPhx2n4FeAI4kia7plXiHExDrmm6Lq+mt2PSK4APAben8oHXs3SdbwdOl6Qq8RcZ42Aa9v9IUhfwEeDG9F400bUcLM4h5Ho9nVyKdVhEbErbzwOHpe0jgefK9utLZYOV18Ml6Vb45lJTU5V46hpnakZ4D9lfsk17TQfECU12TVPzyIPAZrJfEE8DL0ZEf4U6d8WTPn8JmFR0nANjjIjStfxiupb/IGncwBgHxFKPf/MvA38B7EzvJ9Fk13KQOEsKv55OLnUS2f1ls477vg74FWAmsAm4urHh7CbpAODbwKcj4uXyz5rpmlaIs+muaUS8HREzgS6yv5CPa3BIv2BgjJJOAC4ni/VXyZpmPtPAEJH0UWBzRKxtZOHWslUAAANOSURBVBxDqRJnXa6nk0uxXki3laSfm1P5RuCosv26Utlg5YWKiBfSf+qdwL+w+9a8oXFKGkP2C3tpRHwnFTfdNa0UZ7Ne0xTbi8C9wHvJmj5GV6hzVzzp84OBrfWKsyzGOanpMSLiTeAWGn8t3wfMlfQsWfPlh4Cv0HzX8hfilPSvdbuee9NR5NcvdIhNZc+O8v/Dnp3Pf5e2P8KeHWc/it0dZ8+QdZpNTNuH1CHOw8u2/ztZOzDA8ezZ4bierON5dNqexu7O5+NzjlHAEuDLA8qb6ppWibOprikwGZiQtvcF/h/wUeBb7NkJfXHaXsSendDLqsVfcIyHl13rLwNXNsP/o1TXB9jdUd4013KIOOtyPXM/iU59Ad8ga/7YQdYmuYCsXfVuYB3wg9I/SPrHu5aszfsRoLvsOH9A1rHXC1xUpzi/luJ4GFjOnr8Y/zLF+RRwVln52WQjo54G/rKAON9P1uT1MPBgep3dbNe0SpxNdU2BE4GfpHgeBf46lR8N/Chdm28B41L5+PS+N31+9FDxFxjjPelaPgr8K7tHlDXs/1FZPR9g9y/tprmWQ8RZl+vp6V/MzCx37nMxM7PcObmYmVnunFzMzCx3Ti5mZpY7JxczM8udk4tZg0n6yzQL8MNpltpTJd0oaUajYzMbKQ9FNmsgSe8F/h74QES8KelQYGxE/LTBoZntFd+5mDXW4cDPIpuKg4j4WUT8VNJ9krolzS1bd+MpSc8ASJol6d8lrZW0sjQljlmzcHIxa6y7gKMk/aekr0r6jfIPI2J5RMyMbDLHh4AvpbnM/pFsTY5ZwM3AF+seuVkVo4fexcyKEhGvSpoFnAZ8EPimKqxCKekvgNcj4to0U/AJwKpsWRBGkU3pY9Y0nFzMGiwi3gbuA+6T9Agwv/xzSb8JnEu2iihkc0A9FhHvrWecZsPhZjGzBpJ0rKTpZUUzgQ1ln08hm0zw3Ih4PRU/BUxOgwGQNEbS8fWK2awWvnMxa6wDgH+UNAHoJ5t1diG7l8v9fbKZoP8tNYH9NCLOlvQx4BpJB5P9P/4y8FidYzcblIcim5lZ7twsZmZmuXNyMTOz3Dm5mJlZ7pxczMwsd04uZmaWOycXMzPLnZOLmZnl7v8DbVVtIre77MwAAAAASUVORK5CYII=\n" }, "metadata": { "needs_background": "light" } } ] }, { "cell_type": "code", "source": [ "def fit_model(df):\n", " X = df[size_label].values\n", " X = np.reshape(X, (len(X),1)) #Reshape array X in order to make it compatible with the sklearn.Ridge.fit method\n", " Y = df[price_label].values\n", "\n", " model = Ridge(alpha=1.0)\n", " model.fit(X,Y)\n", " # Returns the predicted y values (using linear model)\n", " return X,model.predict(X)" ], "metadata": { "id": "nFBtTlaR-D-j" }, "execution_count": 365, "outputs": [] }, { "cell_type": "code", "source": [ "def plot_ridge_function(df):\n", " X,y_pred = fit_model(df)\n", " plot_df(df)\n", " plt.plot(X,y_pred,label='Ridge Function')\n", " plt.legend()\n", " plt.show()\n", "\n", "plot_ridge_function(dataset)" ], "metadata": { "colab": { "base_uri": "https://localhost:8080/", "height": 295 }, "id": "4471FVu_RMYp", "outputId": "6263d510-bc26-4cd0-9d61-e2671ecaacbf" }, "execution_count": 366, "outputs": [ { "output_type": "display_data", "data": { "text/plain": [ "
" ], "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZcAAAEWCAYAAACqitpwAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nO3de5xVVfnH8c8jV1HkJhmCAiYqSIAMAqbkFUQyvKRmaaGhpGJa/fylZile6ielppiXzBsWiqaZZpKCSmkGOijhBQ1UUBAFuaiIt2Ge3x97DXPOzJkzZ2bOPtfv+/Wa15zznH32es5GzzN7rbXXNndHREQkm7bKdwIiIlJ6VFxERCTrVFxERCTrVFxERCTrVFxERCTrVFxERCTrVFykLJnZjWb28zSvu5ntmsucCp2Z7WxmG82sVb5zkcJnus5FSpGZLQN2ADYDG4G/A2e6+8YM3+9AP3dfmuW85gIjgSrgE+CfwGR3X5XNdkTyTWcuUsq+7u7bAkOAvYDz85xPjTNDXrsC2wJXZLsBM2ud7X2KNIWKi5Q8d38HeISoyABgZreb2WUJz//XzFaZ2dtm9r3E95tZNzP7q5l9YGbPmtllZvZUwut7mNlsM1tnZq+a2XEZ5rUB+EudvBrcVwZ5uJlNNrMlwJIQO9zMFprZBjN72swGJWx/rpmtNLMPQ1sHh/hwM6sM7bxrZleFeJ/QRuvwfEczezDkutTMTk3Y9xQzu8fM7gj7f8nMhmVyXKQ0qLhIyTOzXsBhQMouLjMbC5wDjAb6AYfU2eQ64CPgi8CE8FPz3m2A2cCdwBeA44HrzWxABnl1A46uySuDfTWYR4IjgRHAADPbC7gV+D7QDfgd8KCZtTOz3YEzgb3dvSNwKLAs7OMa4Bp33w74EnBPAx9hJrAC2BE4BvilmR2U8Pr4sE1n4EHgt40dEykdKi5Syv5iZh8CbwGrgYsa2O444DZ3f9HdPwKm1LwQBq+/AVzk7pvc/WVgesJ7DweWuftt7l7l7s8D9wHHpslrmpm9D7wHbA/8oLF9ZZBHjf9z93Xu/jEwCfidu893983uPh34lGjMZzPQjqgItXH3Ze7+WtjH58CuZra9u29093l1GzGznYB9gXPd/RN3XwjcDHw3YbOn3P1hd98M/AEYnOaYSIlRcZFSdmT4q/wAYA+iL/JUdiQqQDWWJzzuDrSu83ri497AiNDttMHMNgAnEJ1dNOQsd+8EDAK6AL0y2FdjeTSU2//U2d9OwI5hosIPiQrpajObaWY7hvdNBHYDXgndb4enaGdHYJ27f5gQWw70THj+TsLjTUB7jQWVDxUXKXnu/g/gdhoeOF9F9KVbY+eEx2uIZnb1SoglbvsW8A9375zws627n55BXi8AlwHXmZk1sq/G8tiy2zq5/aLO/jq4+12h/TvdfT+iIuTA1BBf4u7fIuqamwrcG7rsEr0NdDWzjgmxnYGVjX1uKQ8qLlIurgZGm1mqrpl7gJPMbICZdSCh+yx06fwZmGJmHcxsD5K7fh4CdjOz75hZm/Czt5n1zzCv6URTpsen21cGeaTye+A0MxthkW3M7Gtm1tHMdjezg8ysHdGU6I+BagAzO9HMurt7NbAh7Ks6ccfu/hbwNPB/ZtY+TBSYCPwxw88tJU7FRcqCu68B7gAuTPHaLKLi8zjR4PrjdTY5E+hE1M3zB+AuorELQrfQGKLB97fDNlOJxjMyyeszogH0n2ewrwbzaGDflcCpRAPp68NnOym83A64nGjc5x2is5SaqdpjgZfMbGPI7fgwhlPXt4A+Idf7icaD5mTyuaX06SJKkSYys6nAF9091WytsstDJBWduYg0Ilx7Mih0LQ0n6v65v1zzEMmEZm6INK4jURfUjsC7wJXAA2Wch0ij1C0mIiJZp24xERHJOnWLBdtvv7336dMn32mIiBSVBQsWvOfu3evGVVyCPn36UFlZme80RESKipktTxVXt5iIiGSdiouIiGSdiouIiGSdxlzS+Pzzz1mxYgWffPJJvlORDLVv355evXrRpk2bfKciUtZUXNJYsWIFHTt2pE+fPkSL1kohc3fWrl3LihUr6Nu3b77TESlrsXWLhVVXFyb8fGBmPzSzrhbdxnVJ+N0lbG9mNi3cLnWRmQ1N2NeEsP0SM0u8C2CFmb0Q3jMtLFtOQ2001SeffEK3bt1UWIqEmdGtWzedaYoUgNiKi7u/6u5D3H0IUEF0s6D7gfOAx9y9H/BYeA7RbWj7hZ9JwA0QFQqiJdBHAMOBixKKxQ1Eq77WvG9siDfURpOpsBQX/XuJZG7B8vVc98RSFixfn/V952pA/2DgNXdfDhxB7e1ZpxPd85sQv8Mj84DOZtaD6N7es8OtW9cT3WN8bHhtO3ef59EaNnfU2VeqNkREhKiwnHDzPK589FVOuHle1gtMrorL8UQL7gHs4O6rwuN3iG6UBNHtURNv0boixNLFV6SIp2sjiZlNMrNKM6tcs2ZNkz9ULrRq1YohQ4YwcOBAjj32WDZt2pRyu6985Ss5zkxEitm819fyWVU11Q6fV1Uz7/W1Wd1/7MXFzNoS3WXvT3VfC2ccsa6cma4Nd7/J3Ye5+7Du3eutXlAQtt56axYuXMiLL75I27ZtufHGG5Ner6qqAuDpp5/OR3oiUqRG7tKNtq23opVBm9ZbMXKXblndfy7OXA4DnnP3d8Pzd0OXFuH36hBfSfI9wXuFWLp4rxTxdG3ELs4+zFGjRrF06VLmzp3LqFGjGD9+PAMGDABg22233bLd1KlT+fKXv8zgwYM577xouOm1115j7NixVFRUMGrUKF555ZWs5ycixaOidxdmnDKSH4/ZnRmnjKSid7PmPTUoF1ORv0VtlxjAg8AEolusTqD2fhQPAmea2Uyiwfv33X2VmT0C/DJhEH8McL67rwsz0EYC84nuJ35tI23EqqYP87Oqatq23iqr/2BVVVXMmjWLsWOjOQvPPfccL774Yr0pt7NmzeKBBx5g/vz5dOjQgXXr1gEwadIkbrzxRvr168f8+fM544wzePzxunfzFZFyUtG7S9aLSo1Yi4uZbQOMBr6fEL4cuMfMJgLLgeNC/GFgHNF9vjcBJwOEInIp8GzY7hJ3XxcenwHcDmwNzAo/6dqIVao+zJb+w3388ccMGTIEiM5cJk6cyNNPP83w4cNTXssxZ84cTj75ZDp06ABA165d2bhxI08//TTHHnvslu0+/bTBW6+LiLRYrMXF3T8CutWJrSWaPVZ3WwcmN7CfW4FbU8QrgYEp4inbiFtNH+bnVdVZ68OsGXOpa5tttsl4H9XV1XTu3DnlfkRE4qC1xbIo7j7MTIwePZrbbrtty6yydevWsd1229G3b1/+9KdoToW785///CfnuYlI+VBxybKK3l2YfOCueSksAGPHjmX8+PEMGzaMIUOGcMUVVwAwY8YMbrnlFgYPHsyee+7JAw/o1usiEh+LeqNk2LBhXvdmYYsXL6Z///55ykiaS/9uIrljZgvcfVjduM5cRETKkTs8MBmmdILPUl+c3RJaFVlEpNwsmQMzvlH73DdnvQkVFxGRcvHRWvj1LrXPO+0MZz4DbbbOelMqLiIipc4d7v8+LLq7Nvb9f0KPwbE1qeIiIlLKXv073PXN2ucH/Ry+ek7szaq4iIiUoo1r4Ipda593/RKc8W9o3S4nzWu2WIFLXHL/61//Ohs2bADg7bff5phjjkn5ngMOOIC606qbY+7cuXTq1IkhQ4YwZMgQDjnkkBbvM9HVV1+ddAuBcePGbfl8ItJM7vCnk5ILy2n/grOey1lhARWXgpe45H7Xrl257rrrANhxxx259957Y29/1KhRLFy4kIULFzJnzpys7rtucXn44Yfp3LlzVtsQKSuL/woXd4aX7o+ej74UprwPX6y3SlbsVFyKyD777MPKldFdBZYtW8bAgdF/MB9//DHHH388/fv356ijjuLjjz/e8p5bbrmF3XbbjeHDh3Pqqady5plnArBmzRq+8Y1vsPfee7P33nvzr3/9K6McEtsFuOKKK5gyZQoQnTGde+65DB8+nN12240nn3wSgM2bN3POOecwcOBABg0axLXXXsu0adN4++23OfDAAznwwAMB6NOnD++99x4AV111FQMHDmTgwIFcffXVW9ru378/p556KnvuuSdjxoxJ+qwiZevDd6LrVe4+MXrefQ/42RrY96y8paQxl0zNOg/eeSG7+/zil+GwyzPadPPmzTz22GNMnDix3ms33HADHTp0YPHixSxatIihQ4cCUdfZpZdeynPPPUfHjh056KCDGDw4mh1y9tln86Mf/Yj99tuPN998k0MPPZTFixfX2/eTTz65ZVXmY489lhNOOCFtnlVVVTzzzDM8/PDDXHzxxcyZM4ebbrqJZcuWsXDhQlq3bs26devo2rUrV111FU888QTbb7990j4WLFjAbbfdxvz583F3RowYwf7770+XLl1YsmQJd911F7///e857rjjuO+++zjxxBMzOoYiJcc9KiivPFQbO2M+fGGP/OUUqLgUuJol91euXEn//v0ZPXp0vW3++c9/ctZZ0V8ogwYNYtCgQQA888wz7L///nTt2hWIisN///tfIFqa/+WXX96yjw8++ICNGzcm3XQMom6xhx6q/Q932bJlafM9+uijAaioqNiy7Zw5czjttNNo3Tr6z60mn4Y89dRTHHXUUVtWfj766KN58sknGT9+PH379t1S7BLbECk7L90fja3UGHs5jDw9b+nUpeKSqQzPMLKtZsxl06ZNHHrooVx33XVbCklLVFdXM2/ePNq3b9+k97Vu3Zrq6uotzz/55JOk19u1iwYMW7VqteUWzNlUs/+aNtQtJmXng7fhqoS18774ZTj1CWjVJn85paAxlyLRoUMHpk2bxpVXXlnvS/urX/0qd955JwAvvvgiixYtAmDvvffmH//4B+vXr6eqqor77rtvy3vGjBnDtddeu+V5pvd62WGHHVi9ejVr167l008/TTqracjo0aP53e9+tyXvmrtjduzYkQ8//LDe9qNGjeIvf/kLmzZt4qOPPuL+++9n1KhRGeUnUrKqq2HGscmFZfKzcNpTBVdYQMWlqOy1114MGjSIu+66Kyl++umns3HjRvr378+FF15IRUUFAD179uSnP/0pw4cPZ99996VPnz506tQJgGnTplFZWcmgQYMYMGAAN954Y0Y5tGnThgsvvJDhw4czevRo9tij8b7dU045hZ133plBgwYxePDgLYVw0qRJjB07dsuAfo2hQ4dy0kknMXz4cEaMGMEpp5zCXnvtlVF+IiXprm/BJV1gyaPR83FXRLPAuu+W37zS0JL7QakuuV8zjlJVVcVRRx3F9773PY466qh8pxWrUvh3EwHgzXlw66G1z3sMgVMeg1aFM6LR0JL7hZOhxGLKlCnMmTOHTz75hDFjxnDkkUfmOyURaczmKri0zm3ST7gP+mX3QuY4qbiUuJo7UYpIkZg+Ht74R+3zHYfCpCfyl08zqbg0wt0xs3ynIRlSN68UrTeehOmHJ8d+tgZat81PPi2k4pJG+/btWbt2Ld26dVOBKQLuztq1a5s8vVokr6o+g8u6J8e++wDsckA+sskaFZc0evXqxYoVK1izZk2+U5EMtW/fnl69euU7DZHM3DIG3ppf+7z3vnDyw/nLJ4tUXNJo06YNffv2zXcaIlJqlj4Gfzw6Ofbz9wryepXmUnEREcmVqk/hsi8kx06eBb2/kp98YqTiIiKSCzful7z47ZcOhu/8OX/5xEzFRUQkTv99BO48Ljn287UFdSFkHEr704mI5MvnH8MvvpgcmzgHdto7P/nkWKxri5lZZzO718xeMbPFZraPmXU1s9lmtiT87hK2NTObZmZLzWyRmQ1N2M+EsP0SM5uQEK8wsxfCe6ZZmC/cUBsiIjlx7bDkwrLH4dFaYGVSWCD+hSuvAf7u7nsAg4HFwHnAY+7eD3gsPAc4DOgXfiYBN0BUKICLgBHAcOCihGJxA3BqwvvGhnhDbYiIxGfxX6M7Qq5dUhu7cB0cPyN/OeVJbN1iZtYJ+CpwEoC7fwZ8ZmZHAAeEzaYDc4FzgSOAOzy6xHpeOOvpEbad7e7rwn5nA2PNbC6wnbvPC/E7gCOBWWFfqdoQkQwtWL6eea+vZeQu3ajorZP/tD77CH65Y3Js0lzYsXxX845zzKUvsAa4zcwGAwuAs4Ed3H1V2OYdYIfwuCfwVsL7V4RYuviKFHHStCEiGViwfD0n3DyPz6qqadt6K2acMlIFpiG/GQjvJ3xFDTwGjrklf/kUiDi7xVoDQ4Eb3H0v4CPqdE+Fs5RYF4NK14aZTTKzSjOr1FX4IrXmvb6Wz6qqqXb4vKqaea+vzXdKhefZW6IusMTCcuF6FZYgzuKyAljh7jVrG9xLVGzeDd1dhN+rw+srgZ0S3t8rxNLFe6WIk6aNJO5+k7sPc/dh3bt3T7WJSFkauUs32rbeilYGbVpvxchdujX+pnLx0dqoqPztx7Wx056KBuy30v0Xa8R2JNz9HeAtM9s9hA4GXgYeBGpmfE0AHgiPHwS+G2aNjQTeD11bjwBjzKxLGMgfAzwSXvvAzEaGWWLfrbOvVG2ISAYqendhxikj+fGY3dUllmhKJ/j1LrXPO+0cFZUvfjl/ORWouK9z+QEww8zaAq8DJxMVtHvMbCKwHKi5uuhhYBywFNgUtsXd15nZpcCzYbtLagb3gTOA24GtiQbyZ4X45Q20ISIZqujdRUWlxiMXwL9/mxy7aANotfQG6TbHQarbHItImfvwXbiyzn3qJzwEfUflJ58CpNsci0jWlfR05Smdkp932hl+9ELqbaUeFRcRaZaSna780I+hss6ML3WBNZmKi4g0S6rpykVdXN5fCb8ZkBz73qOw84iUm5f0WVsWqLiISLPUTFf+vKq6+Kcr1+0C674HTJ6feltK+Kwti1RcRKRZaqYrF/Vf7/efDv+5MzmWQRdYyZ21xUDFRUSarWinK69fDtcMSo5NmpvxWmAlddYWExUXESkvdbvAelbAqY83aRclcdYWMxUXESkP90yAl/+SHJvyfrN3V7RnbTmi4iIipW3ta3Dt0OTYaf+CLw7MTz5lQsVFREpX3S6wPqPgpIfyk0uZUXERkdIz41hY8mhyrAVdYNJ0Ki4iUjpWvwLX17nocfIz0H331NtLbFRcRKQ01O0C2+0w+PbM/OSSZcW4GoCKi4gUt9vGwfJ/JcdKqAusWFcDUHERkeK0ahH8rs7S9z94Drp9KT/5xKRYVwNQcRGRLYqm+6VuF9jAb8Axt+Ynl5gV62oAKi4iAhRJ98tNB8DbzyfHSqgLLJViXQ1AxUVEgALvflm5AH5/UHLshy9A553zk0+OFeNqACouIgIUaPeLO1zcOTm214lwxHX5yUcypuIiIkABdr/8dm9477/JsRLvAislKi4iskVBdL+8OQ9uPTQ59uPFsN2O+clHmkXFRUQKQ6ousL1Pha9dkZ98pEVUXEQk/67sDx++nRxTF1hRU3ERkfx5458w/evJsXOWwLZfyE8+kjUqLiKSe9XVcEmdsZ2vnAVjLs1PPpJ1Ki4iklv/tzN8WqfLS11gJUfFRURyY8kcmPGN5NhP3oAOXfOTj8RKxUVE4lW9GS6pU0D2PxcO/Gl+8pGcUHERkfhc3AW8OjmmLrCysFWcOzezZWb2gpktNLPKEOtqZrPNbEn43SXEzcymmdlSM1tkZkMT9jMhbL/EzCYkxCvC/peG91q6NkQkR175W7RycWJhOXe5CksZibW4BAe6+xB3Hxaenwc85u79gMfCc4DDgH7hZxJwA0SFArgIGAEMBy5KKBY3AKcmvG9sI22IFLQFy9dz3RNLWbB8fb5TaZ7NVVFRmfnt2tghU6KisnXnht4lJSgf3WJHAAeEx9OBucC5IX6Huzswz8w6m1mPsO1sd18HYGazgbFmNhfYzt3nhfgdwJHArDRtiBSsoljyPp2691gBnamUsbjPXBx41MwWmNmkENvB3VeFx+8AO4THPYG3Et67IsTSxVekiKdrI4mZTTKzSjOrXLNmTZM/nEg2pVryvii8eF/9wnL+ChWWMhf3mct+7r7SzL4AzDazVxJfdHc3M48zgXRtuPtNwE0Aw4YNizUPkcYU5JL36VR9CpfVuZJ+7OUw8vT85CMFJdbi4u4rw+/VZnY/0ZjJu2bWw91XhW6v1WHzlcBOCW/vFWIrqe3iqonPDfFeKbYnTRsiBavglrxPR11g0ojYusXMbBsz61jzGBgDvAg8CNTM+JoAPBAePwh8N8waGwm8H7q2HgHGmFmXMJA/BngkvPaBmY0Ms8S+W2dfqdoQKWgVvbsw+cBdC7ewzJ2qLjDJSJxnLjsA94fZwa2BO93972b2LHCPmU0ElgPHhe0fBsYBS4FNwMkA7r7OzC4Fng3bXVIzuA+cAdwObE00kD8rxC9voA0RaY7PNsEveyTH9j0bRl+Sn3yk4Fk0OUuGDRvmlZWV+U5DpPCoC0zSMLMFCZeabKEr9EUktUd/Bk9fmxz76Spo2yE/+UhRUXGRsrdg+friGETPlY/Xw9Q+ybEDfgoH6FIxyZyKi5S1or9wMdvUBSZZouIiZS3VhYtxF5eCPFO6ZgisfyM5dsE70Gbr/OQjRU/FRcpari9cLLgzpQ/fgSt3T471/zp884/5yUdKhoqLbFGQf1HHLNcXLubjTKlB6gKTGKm4CFCAf1HnUEXvLjn7rAWxxMvlveGTDcmxC96FNu1zn4uULBUXAQrsL+oSltclXja8CVd/OTk25AQ48vrc5SBlQ8VFgAL5i7pM5PJMaQt1gUmOZVRczGw3ohtz7eDuA81sEDDe3S+LNTvJmaJaNFEyl6qo/Pw9aNUm97lIWcl04crfA+cDnwO4+yLg+LiSkvwo+EUT86jo7hC59rX6hWXEadHZigqL5ECm3WId3P2ZsAhljaoY8hEpOEU32UFdYFIAMi0u75nZl4juLImZHQOsSv8WkdJQNJMdUhWVC9fBVq1yn4uUvUyLy2SiOzbuYWYrgTeAE2PLSqSAFPxkh3dfhhv2SY599X/hoJ/lJx8RMiwu7v46cEi46ddW7v5hvGmJFI6CnuygLjApUJnOFvsl8Ct33xCedwH+x931p5GUhbxMH04nVVG5aAMkj4uK5E2ms8UOqyksAO6+nuiukSKSSyufq19YDpkSna2osEgByXTMpZWZtXP3TwHMbGugXXxpiUg96gKTIpJpcZkBPGZmt4XnJwPT40lJRJKoC0yKUKYD+lPNbBFwcAhd6u6PxJeWiLD8abjtsOTYuCtg+Kn5yUekCTJeW8zdZwGzYsxFRGqoC0yKXNriYmZPuft+ZvYh4QLKmpcAd/ftYs1OpNyoqEiJSFtc3H2/8LtjbtIRKVMvPwj3fCc5NnYqjDwtP/mItFCj3WJm1gp4yd33yEE+IuVHZytSghotLu6+2cxeNbOd3f3NXCQlUhZUVKSEZTqg3wV4ycyeAT6qCbr7+FiyEillz/weHj4nOTb+tzD0O6m3FylCmRaXn8eahUi50NmKlInGZou1B04DdgVeAG5xd93HRaSpVFSkzDS2tth0YBhRYTkMuLKpDZhZKzN73sweCs/7mtl8M1tqZnebWdsQbxeeLw2v90nYx/kh/qqZHZoQHxtiS83svIR4yjZEcm72hfULy9ipKixS8horLgPc/UR3/x1wDDCqGW2cDSxOeD4V+I277wqsByaG+ERgfYj/JmyHmQ0guqXynsBY4PpQsFoB1xEVvQHAt8K26dqQMpS3WxRP6QT/uqZO7H1NL5ay0Fhx+bzmQXO6w8ysF/A14Obw3ICDgHvDJtOBI8PjI6hdr+xe4OCw/RHATHf/1N3fAJYCw8PPUnd/3d0/A2YCRzTShpSZmlsUX/noq5xw87zcFJgpneqfrUx5X2crUlYaKy6DzeyD8PMhMKjmsZl9kMH+rwZ+AlSH592ADQmFagXQMzzuCbwFWwrZ+2H7LfE672konq6NJGY2ycwqzaxyzZo1GXwcKTapblEcm7+eXb+oHHG9ioqUpcau0G/2zbfN7HBgtbsvMLMDmrufOLn7TUS3b2bYsGHeyOZShHJyi2J3uLhz/biKipSxjBeubIZ9gfFmNg5oD2wHXAN0NrPW4cyiF7AybL8S2AlYYWatgU7A2oR4jcT3pIqvTdOGlJnYb1GsWWAiKWV6J8omc/fz3b2Xu/chGpB/3N1PAJ4gmhwAMAF4IDx+MDwnvP64u3uIHx9mk/UF+gHPAM8C/cLMsLahjQfDexpqQ8pQRe8uTD5w1+wWltsPr19YjrpJhUUkiPPMpSHnAjPN7DLgeeCWEL8F+IOZLQXWERUL3P0lM7sHeBmoAia7+2YAMzsTeARoBdzq7i810oZIy6gLTCQjFv2hL8OGDfPKysp8pyGFTF1gIvWY2QJ3H1Y3no8zF5Hi8ut+8NHq5Njxd8Ee4/KTj0gRUHERacjmKrg0xewyna2INErFRSQVdYGJtIiKi0iiVEXlhPug3yG5z0WkiKm4SFYtWL4+vmtK4lT1KVz2hfrxDM9WivZzi8RExUWypmYdr8+qqmnbeitmnDKyOL5oW9gFVrSfWyRGsV1EKeUnp+t4ZUOqBSZP/nuTx1aK7nOL5IDOXCRrcrKOVzZ8tgl+2aN+vJkD9kXzuUVySBdRBrqIMjsKfuwhpllgBf+5RWKiiyglJyp6dynML9dUReX0f8MOA+rHm6FgP7dInqi4SGn7eANM7V0/rmtWRGKl4iIFI+tdS7oQUiRvVFykIDR3Om/KgpSqqPzgOej2pSxnLSINUXGRgpBqOm9jxaVuQbr7hC8xeObe9TcMZysadBfJHRUXKQjNmc6bWJBeaXU8zKyzQUIXmC50FMktFRcpCM25HfHIXbrxertv13/hRy9Dp55JoeacGYlI86m4SMFo0nTedW9QcduQ+vEGBux1oaNIbqm4SPFpxiyw5pwZiUjzqbhIs8U9QF5v/6mKyjlLYdvuGe0vkzMjDfqLZIeKizRL3APkifsf0HolD7X63/obZfmaFQ36i2SPios0S9wD5DX7TzlgH9OFkBr0F8keFRdplrgHyCf/o4LJ7eoEz10OW3fOagCvyQMAABA/SURBVDuJNOgvkj1aFTnQqshNF8v4xIpKuPng+vEcLduiMReRptGqyJJ1WV8JuADWAtPqxiLZoeIijYr9r/lUReWnq6BthybvqtjOPIotX5FMqbhIWrHOoFoyG2YcUz/ezLOVYpvtVWz5ijSFioukFdsMqhi6wIpttlex5SvSFCouklbWZ1ClKio/WwOt27ZsvxTfbK9iy1ekKWKbLWZm7YF/Au2Iiti97n6RmfUlWr+2G7AA+I67f2Zm7YA7gApgLfBNd18W9nU+MBHYDJzl7o+E+FjgGqAVcLO7Xx7iKdtIl69mizUsK+MCi/8Kd59YPx7DhZDFNIZRbPmK1NXQbLE4i4sB27j7RjNrAzwFnA38GPizu880sxuB/7j7DWZ2BjDI3U8zs+OBo9z9m2Y2ALgLGA7sCMwBdgvN/BcYDawAngW+5e4vm9k9qdpIl28pFpeC+eJqQhdYU3IumM8Xg1L+bFJacj4V2aOqtTE8bRN+HDgIqLnsejowBbgBOCI8BrgX+G0oUEcAM939U+ANM1tKVGgAlrr76wBmNhM4wswWp2mjJGTyxVMQg8WpisqF62CrVik3b0rOBfH5YlLKn03Kx1Zx7tzMWpnZQmA1MBt4Ddjg7lVhkxVAzY03egJvAYTX3yfq1toSr/OehuLd0rRRN79JZlZpZpVr1qxpyUfNmZovnisffZUTbp7HguXrU26XarA4Z56fUb+wbNM9OltpoLBA03LO6+eLWSl/NikfsQ7ou/tmYIiZdQbuB/aIs72mcvebgJsg6hbLczoZyXSGUd4Gi1swC6wpOZfyYHgpfzYpHzmZLebuG8zsCWAfoLOZtQ5nFr2AlWGzlcBOwAozaw10IhrYr4nXSHxPqvjaNG0UvUy/eHJ+/5JUReWiDWCW8S6aknMp35+llD+blI84B/S7A5+HwrI18CgwFZgA3Jcw2L7I3a83s8nAlxMG9I929+PMbE/gTmoH9B8D+gFGNKB/MFHxeBb4tru/ZGZ/StVGunyLaUC/oAZ7K2+Fh36UHOszCk56KD/5ZKigjqFIEcvH2mI9gOlm1opobOced3/IzF4GZprZZcDzwC1h+1uAP4QB+3XA8QChWNwDvAxUAZNDdxtmdibwCNFU5Fvd/aWwr3MbaKMkFMz6V1m6EDLXX/QaMBeJX5yzxRYBe6WIv07tbK/E+CfAsQ3s6xfAL1LEHwYezrQNyZIsXl2fjy96XRkvEr9YZ4tJiUk1C2zEaS26EDIfM6Nqxq1aGRowF4mJln+RzMS0HH4+ZkZpwFwkfrpZWFBMA/o5lYN7rGhwXaR46WZh0jTzboS/n5scO/BnsP//Zr2pgpmgICJZo+Ii9RXAHSFFpLipuEgtFRURyRLNFhOYO7V+YTns1yosItJsOnMpIc0aGC/CsxVNABApfCouJaLJFyMWYVEBXV0vUizULVYiMr4Y8d/X1y8sR95QkIVlwfL1XPfE0qTbCmg5epHioDOXEtHoxYjucHHn+m/MQ1Fpyc3OtBy9SHFQcSlAzRlTSHvVeSNdYLkcw8i0W6uh9b90db1IcVBxKTAtGVOodzHik1fCY5ckb3TS36DPfllprzmycbMzXXQpUvhUXApMVlbsbUIXWK5XCC7Ym52JSFapuBSYFo8pNHEWWEPtZTouktXuuxTbqqiIFCctXBkU0sKVzfpif/q38OgFyRud+jj0rGhye5l0leVjSrCubxEpPFq4sog09hd74hd7u9awuNW362/UhFlgddvLpKss191pur5FpLiouBShmi/219u1rKg0JJOuuVxPCdbdI0WKi4pLETq0VSWT230/OXjmAth+1ybtp6FupkzGRXI94K7rW0SKi8ZcgkIac2lQdTVckvwl/l6/b7J8v6lN/pIvxm4mjbmIFB6NuRS7X/SAzzclx6a8z/JmFoli7GbS7DGR4qG1xQrdq3+PphcnFpafvLFlbKW5a23VdDO1MtTNJCJZpzOXQlW9GS7pmhw78ALY/ydJoeaORegiRRGJk8ZcgoIac7mkG1RXJcfSzALTWISI5IvGXIrBikq4+eDk2HlvQfvt0r5NYxEiUmhUXApBqi6w8dfC0O9mvSmd5YhILqi45Nsfj4Gls2uff2FPOOPpWJoqxunHIlKcVFzyZfm/4baxybEL3oU27WNrshinH4tIcVJxybVUXWAn3Af9Dom9aV3lLiK5Ett1Lma2k5k9YWYvm9lLZnZ2iHc1s9lmtiT87hLiZmbTzGypmS0ys6EJ+5oQtl9iZhMS4hVm9kJ4zzQzs3Rt5N3zM5ILS8+KaBZYDgoL1E4//vGY3dUlJiKxim0qspn1AHq4+3Nm1hFYABwJnASsc/fLzew8oIu7n2tm44AfAOOAEcA17j7CzLoClcAwwMN+Ktx9vZk9A5wFzAceBqa5+ywz+1WqNtLlG+tU5PXL4JrBtc/bbhtdCNm6bYt33dIBeg3wi0hL5HwqsruvAlaFxx+a2WKgJ3AEcEDYbDowFzg3xO/wqNrNM7POoUAdAMx293Xhg8wGxprZXGA7d58X4ncQFa9ZadrIrerNcPvX4M1/18bOWghd+2Zl9y0doNcAv4jEJSfLv5hZH2AvojOMHULhAXgH2CE87gm8lfC2FSGWLr4iRZw0bdTNa5KZVZpZ5Zo1a5r+wdJZMD3qAqspLEdcH3WBZamwQPOXfsnW+0VEGhL7gL6ZbQvcB/zQ3T8IwyIAuLubWaxLBKRrw91vAm6CqFssKw2ufQ2uHVr7vO9X4Tt/ga1aZWX3iVo6QK8BfhGJS6zFxczaEBWWGe7+5xB+18x6uPuq0O21OsRXAjslvL1XiK2ktourJj43xHul2D5dG/HZXAW3joGVC2pjP3wBOu8cW5MtXR9M64uJSFxiKy5h5tYtwGJ3vyrhpQeBCcDl4fcDCfEzzWwm0YD++6E4PAL8MmHG1xjgfHdfZ2YfmNlIou627wLXNtJGPJ69Gf72P7XPj7oJBn8z1iZrtHTpFy0dIyJxiPPMZV/gO8ALZrYwxH5K9IV/j5lNBJYDx4XXHiaaKbYU2AScDBCKyKXAs2G7S2oG94EzgNuBrYkG8meFeENtZN+jP4enp0WPdz0Evv0n2Ep3MhCR8hbnbLGnAGvg5YPrBsIssckN7OtW4NYU8UpgYIr42lRtxGFp+4Fsv/XOvPn1uxk0YEAumhQRKXj6E7sFFixfz+Gzt2Pohss57q7lLFi+Pt8piYgUBBWXFtBUXhGR1FRcWkC3ChYRSU0LV7ZAc6byarkVESkHKi4t1JSpvFpuRUTKhbrFckhjNCJSLlRcckhjNCJSLtQtlkNabkVEyoWKS45puRURKQfqFhMRkaxTcRERkaxTcRERkaxTcRERkaxTcRERkaxTcRERkayz6DYqYmZriG4slgvbA+/lqK2WUJ7ZpTyzS3lmV3Pz7O3u3esGVVzywMwq3X1YvvNojPLMLuWZXcozu7Kdp7rFREQk61RcREQk61Rc8uOmfCeQIeWZXcozu5RndmU1T425iIhI1unMRUREsk7FRUREsk7FJUvM7FYzW21mLybEuprZbDNbEn53CXEzs2lmttTMFpnZ0IT3TAjbLzGzCTnKc4qZrTSzheFnXMJr54c8XzWzQxPiY0NsqZmdF0OeO5nZE2b2spm9ZGZnh3hBHdM0eRbUMTWz9mb2jJn9J+R5cYj3NbP5oc27zaxtiLcLz5eG1/s0ln+MOd5uZm8kHMshIZ63/49CG63M7Hkzeyg8L5hj2UieuTme7q6fLPwAXwWGAi8mxH4FnBcenwdMDY/HAbMAA0YC80O8K/B6+N0lPO6SgzynAOek2HYA8B+gHdAXeA1oFX5eA3YB2oZtBmQ5zx7A0PC4I/DfkE9BHdM0eRbUMQ3HZdvwuA0wPxyne4DjQ/xG4PTw+AzgxvD4eODudPnHnOPtwDEpts/b/0ehnR8DdwIPhecFcywbyTMnx1NnLlni7v8E1tUJHwFMD4+nA0cmxO/wyDygs5n1AA4FZrv7OndfD8wGxuYgz4YcAcx090/d/Q1gKTA8/Cx199fd/TNgZtg2m3mucvfnwuMPgcVATwrsmKbJsyF5OabhuGwMT9uEHwcOAu4N8brHs+Y43wscbGaWJv84c2xI3v4/MrNewNeAm8Nzo4COZUN5NiKrx1PFJV47uPuq8PgdYIfwuCfwVsJ2K0KsoXgunBlOhW+t6WpKk09O8wzdCHsR/SVbsMe0Tp5QYMc0dI8sBFYTfUG8Bmxw96oUbW7JJ7z+PtAt7jzr5ujuNcfyF+FY/sbM2tXNsU4uufg3vxr4CVAdnnejwI5lA3nWiP14qrjkiEfnl4U67/sG4EvAEGAVcGV+06llZtsC9wE/dPcPEl8rpGOaIs+CO6buvtndhwC9iP5C3iPPKdVTN0czGwicT5Tr3kRdM+fmMUXM7HBgtbsvyGcejUmTZ06Op4pLvN4Np5WE36tDfCWwU8J2vUKsoXis3P3d8D91NfB7ak/N85qnmbUh+sKe4e5/DuGCO6ap8izUYxpy2wA8AexD1PXROkWbW/IJr3cC1uYqz4Qcx4auR3f3T4HbyP+x3BcYb2bLiLovDwKuofCOZb08zeyPOTueLRko0k+9AbE+JA+U/5rkwedfhcdfI3ng7BmvHTh7g2jQrEt43DUHefZIePwjon5ggD1JHnB8nWjguXV43Jfawec9s5yjAXcAV9eJF9QxTZNnQR1ToDvQOTzeGngSOBz4E8mD0GeEx5NJHoS+J13+MefYI+FYXw1cXgj/H4W2DqB2oLxgjmUjeebkeGb9Q5TrD3AXUffH50R9khOJ+lUfA5YAc2r+QcI/3nVEfd4vAMMS9vM9ooG9pcDJOcrzDyGPRcCDJH8xXhDyfBU4LCE+jmhm1GvABTHkuR9Rl9ciYGH4GVdoxzRNngV1TIFBwPMhnxeBC0N8F+CZcGz+BLQL8fbh+dLw+i6N5R9jjo+HY/ki8EdqZ5Tl7f+jhHYOoPZLu2COZSN55uR4avkXERHJOo25iIhI1qm4iIhI1qm4iIhI1qm4iIhI1qm4iIhI1qm4iOSZmV0QVgFeFFapHWFmN5vZgHznJtJcmooskkdmtg9wFXCAu39qZtsDbd397TynJtIiOnMRya8ewHseLcWBu7/n7m+b2VwzG2Zm4xPuu/Gqmb0BYGYVZvYPM1tgZo/ULIkjUihUXETy61FgJzP7r5ldb2b7J77o7g+6+xCPFnP8D3BFWMvsWqJ7clQAtwK/yHnmImm0bnwTEYmLu280swpgFHAgcLeluAulmf0E+NjdrwsrBQ8EZke3BaEV0ZI+IgVDxUUkz9x9MzAXmGtmLwATEl83s0OAY4nuIgrRGlAvufs+ucxTpCnULSaSR2a2u5n1SwgNAZYnvN6baDHBY9394xB+FegeJgNgZm3MbM9c5SySCZ25iOTXtsC1ZtYZqCJadXYStbfLPYloJei/hC6wt919nJkdA0wzs05E/x9fDbyU49xFGqSpyCIiknXqFhMRkaxTcRERkaxTcRERkaxTcRERkaxTcRERkaxTcRERkaxTcRERkaz7fwS1bdssA5XKAAAAAElFTkSuQmCC\n" }, "metadata": { "needs_background": "light" } } ] }, { "cell_type": "markdown", "source": [ "###Let's try to use another Dataset###" ], "metadata": { "id": "E8pmFfJXpX6o" } }, { "cell_type": "code", "source": [ "n_samples = 50\n", "X = np.sort(np.random.rand(n_samples,1),axis=0)\n", "f = lambda x: np.sin(2*np.pi*x)\n", "# Adding Gaussian noise\n", "noise = np.random.normal(0,0.2,size=(n_samples,1))\n", "Y = f(X) + noise\n", "plt.plot(X,Y,'ro')\n", "plt.show()" ], "metadata": { "colab": { "base_uri": "https://localhost:8080/", "height": 265 }, "id": "bSk5udgypeyV", "outputId": "68f77fb3-6216-4443-a31b-4f30b6f62b5d" }, "execution_count": 371, "outputs": [ { "output_type": "display_data", "data": { "text/plain": [ "
" ], "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXwAAAD4CAYAAADvsV2wAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAUZ0lEQVR4nO3df4xlZ13H8c9nd1vIBAS6O2JtuztVi6GCCp00EBN/pMWU/aMr8iPFSyimuCkNhogxNplEDWYjaFRC0loWaCh0pEUSZAxrGhFIE2VJb4PUtlgd193tFrTTXWxCNrCW+frHOdPeuXvvnTu95/fzfiWTuffc03uec2f7uec8z/c8xxEhAED37ai7AQCAahD4AJAIAh8AEkHgA0AiCHwASMSuuhswzp49e2JhYaHuZgBAqzz44INPRcT8qNcaG/gLCwvq9/t1NwMAWsX2iXGv0aUDAIkg8AEgEQQ+ACSCwAeARBD4AJAIAr8LlpelhQVpx47s9/Jy3S0C0ECNLcvElJaXpYMHpbNns+cnTmTPJanXq69dABqHI/y2W1p6Luw3nD2bLQeAAQR+2508OXr5ibHXXgBIFIHfdnv3jl5u05cPYBMCv+0OHcrCfVgE3ToANiHw267Xy8J9lHHdPQCSROB3wb59o5eP6+4BkCQCvwsOHZLm5jYvm5vLlgNAjsDvgl5POnw4O9K3s9+HD1OHD2ATLrzqil6PgAcwEUf4AJAIAh8AEkHgA0AiCHwASASBDwCJIPABIBEEPgAkgsAHgEQQ+ACQCAIfABJRSODbvtP2k7YfHvO6bX/E9qrth2y/tojtAgCmV9QR/iclXTfh9TdKuiL/OSjprwraLqaxvCwtLEg7dmS/uRMWkKRCAj8i7pd0ZsIqByR9KjJHJb3U9sVFbBtbWF6WDh7M7nEbkf0+eJDQBxJUVR/+JZIeH3h+Kl+2ie2Dtvu2+2traxU1reOWlqSzZzcvO3uW2x8CCWrUoG1EHI6IxYhYnJ+fr7s53TDuNofc/hBITlWB/4SkywaeX5ovQ9nG3eaQ2x8Cyakq8FckvTOv1nmdpKcj4jsVbTtt3P4QQK6osszPSPqapJ+2fcr2TbZvtn1zvsoRScckrUr6mKRbithuq1VVOdPk2x9SPQRUyhFRdxtGWlxcjH6/X3czyrFROTM4mDo315wgrgKfAVAK2w9GxOLI1wj8GiwsZOWRw/btk44fr7o19eAzAEoxKfAbVaWTjDIqZ9rWPUL1EFA5Ar8ORVfOtPHiKqqHgMoR+HUounKmqourijyLoHoIqByBX4eiK2eq6B4p+iyiydVDQEcxaNsFVQyAjtvGxnYOHSo+rJeXs7OUkyezrp4ytgF0DIO2XVdF98iks4UyxgzaOC4BNByB3wVVdI9sNZha9JgBk74BhaNLB9MZdaHUMFtaXy9mezt2ZEf2ZW4D6CC6dDC7wbOIcYosqaRsEygcgY/p9XrZIPDdd5c/ZkDZJlA4Ar/Jmnr1bBVjBpRtAoWjD7+pmFwMwPNAH34bUaUCoGAEflMxuRiAghH4TVVmlUpTxwYAlIrAb6qyqlS4ghVIFoHfVGVVqTA2ACSLKp3UcAUr0GlU6eA5XMEKJIvATw1XsALJIvBTwxWsQLJ21d0A1KDXI+CBBHGEDwCJIPABIBEEPtqBq4OBmdGHj+Ybnjl04+pgibEIYBs4wkfzcXUwUAgCH83HzKFAIQh8NB9XBwOFIPDRfFwdDBSCwEfzTXt1MJU8wESFBL7t62w/ZnvV9q0jXn+X7TXb/5L/vLuI7SIhvZ50/Hg2o+fx46PDnnn+gYlmDnzbOyXdJumNkq6U9HbbV45Y9d6I+Pn85+OzbhfYhEoeYEtFHOFfLWk1Io5FxDlJ90g6UMD7AtOjkgfYUhGBf4mkxween8qXDXuz7Ydsf872ZaPeyPZB233b/bW1tQKahmRQyQNsqapB27+TtBARPyvpHyTdNWqliDgcEYsRsTg/P19R09AJVPIAWyoi8J+QNHjEfmm+7FkRcToifpA//bikqwrYbjGKquygQmT7ivzMmOcf2FpEzPSjbD6eY5Iul3ShpG9K+pmhdS4eePwmSUe3et+rrroqSnf33RFzcxFZXUf2MzeXLa/jfVLCZwaUQlI/xuRqITcxt71f0ocl7ZR0Z0Qcsv2BfMMrtv9E0vWSnpF0RtJ7IuLfJr1nJTcxX1jIyveG7duXlf5V/T4p4TMDSjHpJuaFBH4ZKgn8HTuyY8thdlbvXfX7pITPDCjFpMBP+0rboio7qBDZPj4zoHJpB/6oyo4LLpC+973tDSRSIbJ9fGZA5dIO/OHKjt27s9+nT2/v8nwqRLavzs+MiiokKu0+/GEMJHbf8N2zpOzMgi9odAR9+NPi8vzuY84dJIzAH8RAYvfxpY6EEfiDGEjsPr7UkTACfxCDr+20nUFYvtSRsF11N6Bxej0Cvk2GB2E3Kquk0X/HjWVLS1k3zt69WdjzN0cCqNJBu1FZBWxClQ66a9IgLPX2wCYEPtpt3GDrRRdxj1tgCIGPdhs3CCtRbw8MIfDRbuMqq86cGb0+9fZIGIGP9uv1sgHa9fXsd69XXL094wDoEAIf3VREvf1GySfjAOgIAh/dNM1FdFsdvTPvDjqGOnykaZpZM7krF1qIOnxg2DRH78y7g44h8JGmaWbNZN4ddExagU/FBTZMc/TOZHromHQCn4oLDDp0KLt/8aALLjj/6H1UySfQUukEPhUXGGZPfg50TDqBz52OMGhpSTp3bvOyc+c4AEC9Su52TifwqbjAIA4A0DQVdDunE/hUXGAQBwBomgq6ndMJfCouMIgDADRNBWed6QS+RMUFnsMBAJqmgrPOtAIfGMQBAJqkgrNOAh8AmqCCs85dhb0TAGA2vV6pZ5oc4QNAIgoJfNvX2X7M9qrtW0e8/gLb9+avf932QhHbBSrHfExosZkD3/ZOSbdJeqOkKyW93faVQ6vdJOm7EfFTkv5S0odm3S5QOeZjQssVcYR/taTViDgWEeck3SPpwNA6ByTdlT/+nKRrbCYuQcswHxNarojAv0TS4wPPT+XLRq4TEc9IelrS7uE3sn3Qdt92f21trYCmAQViOga0XKMGbSPicEQsRsTi/Px83c0BNmM6Bgxr2ZhOEYH/hKTLBp5fmi8buY7tXZJeIul0AdsGqsN0DBjUwjGdIgL/AUlX2L7c9oWSbpC0MrTOiqQb88dvkfTlaOrd04FxmI4Bg1o4pjNz4Od98u+VdJ+kb0n6bEQ8YvsDtq/PV/uEpN22VyW9X9J5pZtA7aY5PWc6Bmxo4ZhOIVfaRsQRSUeGlv3BwOPvS3prEdsCSrFxer5xxLZxei4R6hht797s38mo5Q3VqEFboDYtPD1HzVo4pkPgA1IrT89RsxaO6TB5GiC18vQcDVDyZGdF4wgfkFp5eg5sF4EPSK08PQe2iy4dYEPLTs+B7eIIHwASQeADQCIIfABIBIEPAIkg8AFgO1o2JfIgAh8AplX2lMglf5kQ+AAwrTLnXKpgfn0CHwCmVeacSxVM4EfgA8C0yrzNZQUT+BH4ADCtMudcquCeyQQ+AEyrzDmXKpjAj7l0AGA7yppzaeM9l5aybpy9e7OwL3BbBD4ANEXJE/jRpQMA47T4IqtRCHwAkM4P91tuKb0uvmqOiLrbMNLi4mL0+/26mwEgBRsXPQ3WwdtZ0A/bt086fryypm2X7QcjYnHUaxzhA8Coi57GHQy3+Mb2BD4AbCfEt1MX37AxAAIfAMaFuL35+Xbq4iuYG2e7uhf4DftGBdAC4y56uvnm53+RVQVz42xXtwK/gd+oAFpg3BW0t9+eDdCur2e/R4X9uIPMCubG2a5uVeksLGQhP6zho+oAWmpUdc/cXPZl8b73SadPn//f7N4tPfVUaU1Kp0qngd+oADqsgd02k3Qr8CuYbQ4AnjXpIPPMmdGvjVtegW4FfgWzzQHAsyYdZDbwALRbgV/m1KUAMGzSQWYDD0BnGrS1fZGkeyUtSDou6W0R8d0R6/1Q0r/mT09GxPVbvTdTKwBoheXl56Y0vuiibNmZM9mR/P790pEjpU13PMqkQdtZA/9PJZ2JiA/avlXSyyLi90es972IeNF23pvAB9Aqkyp2KuxlKLNK54Cku/LHd0n6tRnfDwDaqQUVO7MG/ssj4jv54/+W9PIx673Qdt/2UdtjvxRsH8zX66+trc3YNACo0LiKnVHXBtVky8C3/SXbD4/4OTC4XmR9Q+P6h/blpxi/IenDtn9y1EoRcTgiFiNicX5+frv7AgD1mTQfT0Ou9t8y8CPi2oh41YifL0j6H9sXS1L++8kx7/FE/vuYpK9Kek1hewAATXDo0PmTrUnZNC8N6daZtUtnRdKN+eMbJX1heAXbL7P9gvzxHkm/IOnRGbcLAM3S6zV+Dv1ZA/+Dkt5g+z8kXZs/l+1F2x/P13mlpL7tb0r6iqQPRgSBD6B79u0bvbwhV/vPFPgRcToiromIK/KunzP58n5EvDt//M8R8eqI+Ln89yeKaDgAbNKEqdEbeLHVoG5daQsgTU2ZGr3hV/t3a3pkAGliavRnpTM9MpCyJnRp1IWp0adC4ANd0JQujbo0cGbKJiLwgS5owWX9pWr4YGlTEPhAF6TepbHdwdJEu7921d0AAAXYu3f0oGVKXRq93nTVMLfcIt1xx3MXSW10f228R4dxhA90AV0a01le3hz2GxLp/iLwgS5oeP13YywtNX76gzIR+EBbDfdDS1nN+fp69puwP/8zmjRVcQLdX/ThA200fHelhPqhpzbqM7JHH+HbSXR/cYQPtFHqZZjTGPUZRZw/hbEt3XxzEl+UBD7QRqmXYU5j3GcRsXms49Oflm6/vdq21YQuHaCNKMPc2rjPKMH5dTZwhA+0EWWYW+MzOg+BD7QRZZhb4zM6D106QFtNe2VpyjY+n6WlrE9/Y1A70c+NwAfQXZSvbkKXDoDuonx1EwIfQHdRvroJgQ90RaJT/k7EjVE2IfCBLkj9jlfjUJq5CYEPdAF91aNRmrmJY9xUoTVbXFyMfr9fdzOAdtixY/ykYOvr1bcHtbH9YEQsjnqNI3ygC+irxhQIfKAL6KvGFAh8oAsG+6olaefO5/rwUx+4xbO40hboio2BSK4sxRgc4QNdQrUOJiDwgS7hylJMQOADXUK1DiYg8IEuoVoHE8wU+LbfavsR2+u2Rxb65+tdZ/sx26u2b51lmwAm4MpSTDBrlc7Dkn5d0kfHrWB7p6TbJL1B0ilJD9heiYhHZ9w2gFG4MQrGmCnwI+JbkmR70mpXS1qNiGP5uvdIOiCJwAeAClXRh3+JpMcHnp/Kl53H9kHbfdv9tbW1CpoGAOnY8gjf9pck/diIl5Yi4gtFNiYiDks6LGWTpxX53gCQui0DPyKunXEbT0i6bOD5pfkyAECFqujSeUDSFbYvt32hpBskrVSwXQDAgFnLMt9k+5Sk10v6ou378uU/bvuIJEXEM5LeK+k+Sd+S9NmIeGS2ZgMAtmvWKp3PS/r8iOXflrR/4PkRSUdm2RYAYDZcaQsAiSDwAcxmeVlaWMhus7iwwPz7DUbgA3j+ob28nM23f+JEdk/djfn3Cf1GIvCB1M0S2sy/3yoEPpC6WUK7ivn36TIqDIEPpG6W0C57/n26jApF4AOpmyW0y55/ny6jQhH4QOpmCe2y59/nlo2FIvCB1M0a2r2edPy4tL6e/S5yLn5u2VgoAh9AuaE9C27ZWCgCH0BzccvGQs16i0MAKBe3bCwMR/gAkAgCHwASQeADQCIIfABIBIEPAIlwRNTdhpFsr0k6MeXqeyQ9VWJzmor9Tgv7nZbnu9/7ImJ+1AuNDfztsN2PiMW621E19jst7HdaythvunQAIBEEPgAkoiuBf7juBtSE/U4L+52Wwve7E334AICtdeUIHwCwBQIfABLRqsC3fZ3tx2yv2r51xOsvsH1v/vrXbS9U38riTbHf77f9qO2HbP+j7X11tLNoW+33wHpvth22W1+6N80+235b/vd+xPZfV93GMkzxb3yv7a/Y/kb+73x/He0smu07bT9p++Exr9v2R/LP5SHbr51pgxHRih9JOyX9p6SfkHShpG9KunJonVsk3ZE/vkHSvXW3u6L9/hVJc/nj96Sy3/l6L5Z0v6SjkhbrbncFf+srJH1D0svy5z9ad7sr2u/Dkt6TP75S0vG6213Qvv+ipNdKenjM6/sl/b0kS3qdpK/Psr02HeFfLWk1Io5FxDlJ90g6MLTOAUl35Y8/J+ka266wjWXYcr8j4isRsXGn56OSLq24jWWY5u8tSX8s6UOSvl9l40oyzT7/lqTbIuK7khQRT1bcxjJMs98h6Ufyxy+R9O0K21eaiLhf0pkJqxyQ9KnIHJX0UtsXP9/ttSnwL5H0+MDzU/myketExDOSnpa0u5LWlWea/R50k7Ijgrbbcr/z09vLIuKLVTasRNP8rV8h6RW2/8n2UdvXVda68kyz338k6R22T0k6Ium3q2la7bb7//9E3PGqQ2y/Q9KipF+quy1ls71D0l9IelfNTanaLmXdOr+s7Ezuftuvjoj/rbVV5Xu7pE9GxJ/bfr2kT9t+VUSs192wNmnTEf4Tki4beH5pvmzkOrZ3KTv1O11J68ozzX7L9rWSliRdHxE/qKhtZdpqv18s6VWSvmr7uLL+zZWWD9xO87c+JWklIv4vIv5L0r8r+wJos2n2+yZJn5WkiPiapBcqm1ys66b6/39abQr8ByRdYfty2xcqG5RdGVpnRdKN+eO3SPpy5CMfLbblftt+jaSPKgv7LvTpSlvsd0Q8HRF7ImIhIhaUjV1cHxH9eppbiGn+jf+tsqN72d6jrIvnWJWNLME0+31S0jWSZPuVygJ/rdJW1mNF0jvzap3XSXo6Ir7zfN+sNV06EfGM7fdKuk/ZqP6dEfGI7Q9I6kfEiqRPKDvVW1U2EHJDfS0uxpT7/WeSXiTpb/Ix6pMRcX1tjS7AlPvdKVPu832SftX2o5J+KOn3IqLVZ7FT7vfvSvqY7d9RNoD7rg4czMn2Z5R9ge/Jxyf+UNIFkhQRdygbr9gvaVXSWUm/OdP2OvCZAQCm0KYuHQDADAh8AEgEgQ8AiSDwASARBD4AJILAB4BEEPgAkIj/B51Bow8v/0a8AAAAAElFTkSuQmCC\n" }, "metadata": { "needs_background": "light" } } ] }, { "cell_type": "code", "source": [ "# I didn't find a way to perform polynomial ridge regression using sklearn.Ridge package, so i found here https://www.kaggle.com/code/elyas19/implementing-polynomial-ridge-regression\n", "# a code example that implements it from scratch\n", "\n", "def ridge_regression_from_scratch(N, x_train, y_train, poly_deg,alpha):\n", " a = alpha\n", " Q_train = np.zeros(shape = (N,poly_deg+1))\n", " Q_train[:,0] = 1\n", " for i in range(1,poly_deg+1):\n", " Q_train[:,i] = np.power(x_train,i).reshape((N,)) \n", "\n", " W = np.linalg.pinv((Q_train.T.dot(Q_train) + a*np.eye(poly_deg+1))).dot(Q_train.T).dot(y_train)\n", " plt.plot(x_train,Q_train.dot(W),'b',label='poly fit alpha =' + str(alpha))\n", " plt.plot(X,Y,'ro')\n", " plt.legend()\n", " plt.show()\n", "\n", "# Trying different values of alpha\n", "ridge_regression_from_scratch(n_samples,X,Y,3,0.001)\n", "ridge_regression_from_scratch(n_samples,X,Y,3,0.01)\n", "ridge_regression_from_scratch(n_samples,X,Y,3,1)\n", "ridge_regression_from_scratch(n_samples,X,Y,3,10)" ], "metadata": { "colab": { "base_uri": "https://localhost:8080/", "height": 1000 }, "id": "Iq3jcLStsIun", "outputId": "e997b9b4-f9b9-493d-a946-3e475ad2868d" }, "execution_count": 372, "outputs": [ { "output_type": "display_data", "data": { "text/plain": [ "
" ], "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXwAAAD4CAYAAADvsV2wAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nO3dfZzM9fr48de163Yj93Vk2eVEFEVRiU4kp5KbOqk4GwpJfVXKzVEbyfnp6FTq6FCR3C5KJSodUpRKtN0pOkVarOjILoXcrH3//njPam9mZmd3PnP7uZ6Pxzx2Z+Yz87lmrGve87653mKMQSmlVPxLiHQASimlwkMTvlJKuYQmfKWUcglN+Eop5RKa8JVSyiUqRDoAX+rWrWtSU1MjHYZSSsWUTz/99GdjTD1v90Vtwk9NTSUzMzPSYSilVEwRke2+7tMuHaWUcglN+Eop5RKa8JVSyiWitg9fKbc5fvw42dnZHDlyJNKhqBhQpUoVkpOTqVixYsCP0YSvVJTIzs6mevXqpKamIiKRDkdFMWMM+/btIzs7m8aNGwf8OO3SiQcZGZCaCgkJ9mdGRqQjUuVw5MgR6tSpo8lelUpEqFOnTpm/DWoLP9ZlZMCQIXD4sL2+fbu9DpCWFrm4VLlosleBKs/firbwY116+u/JvsDhw/Z2pZQqRBN+rNuxw/vt232uvVDKMZ06dSrTAsnFixfTokULOnfuTGZmJnfffTcAa9as4aOPPirTudesWUP37t2DPiZYc+bMoWnTpjRt2pQ5c+Z4PSYnJ4euXbvStGlTunbtSm5uLmD74u+++27OPPNMzj33XD777LOTj7nqqquoWbOmo/Frwo91jRp5v11E+/JV1Jk5cyYzZsxg9erVtG3blilTpgDlS/jRICcnh4cffpj169ezYcMGHn744ZPJvLBJkybRpUsXtmzZQpcuXZg0aRIAb731Flu2bGHLli1Mnz6dO+644+RjRo0axbx58xyNVxN+rJs40Sb34ozRbh1VJllZWTRv3py0tDRatGhB7969OezpLnznnXdo06YNrVq1YuDAgRw9erTIY1944QWGDx9+8vqMGTO49957ixwzYcIEPvjgAwYNGsSoUaNOtr6zsrJ49tlnefLJJ2ndujVr164t8rgNGzbQvn172rRpwyWXXMK3335bIvbx48fTr18/2rdvT9OmTZkxY8bJ+w4ePEjv3r1PvraCXf4mTJhAu3btaNmyJUOGDKE8u/+tWLGCrl27Urt2bWrVqkXXrl35z3/+U+K4pUuXMmDAAAAGDBjAa6+9dvL2/v37IyJcfPHF7N+/n927dwPQpUsXqlevXuaY/NFB21iXlgY33+z9Pl/dPSrqDR8OX3zh7HO2bg1PPeX/mG+//ZaZM2fSoUMHBg4cyLRp0xg2bBi33HIL77zzDs2aNaN///4888wzRRL8jTfeyMSJE3nssceoWLEis2bN4rnnnivy3OPGjePdd9/l8ccfp23btqxZswawdbOGDh1KtWrVGDlyZImYmjdvztq1a6lQoQKrVq3igQce4JVXXilx3MaNG/n44485dOgQbdq04ZprrgHg888/Z9OmTZxxxhl06NCBDz/8kI4dOzJs2DDGjRsHQL9+/XjjjTfo0aNHkefMyMjgscceK3GuM888k5dffpldu3bRsGHDk7cnJyeza9euEsf/9NNP1K9fH4A//OEP/PTTTwA+H19wrNM04ceDlBTvffa+unuU8qFhw4Z06NABgJtvvpkpU6bQtWtXGjduTLNmzQDbQp06dWqRhF+tWjUuv/xy3njjDVq0aMHx48dp1aqVIzEdOHCAAQMGsGXLFkSE48ePez2uV69eVK1alapVq9K5c2c2bNhAzZo1ufDCC0lOTgagdevWZGVl0bFjR1avXs0///lPDh8+TE5ODuecc06JhJ+Wlkaaw7PdRCRis7E04ceDiROLTs0ESEqyt6uYVFpLPFSKJ6KyJKbBgwfzyCOP0Lx5c2699VbHYho7diydO3dmyZIlZGVl0alTJ6/H+Yq9cuXKJ29LTEwkLy+PI0eOcOedd5KZmUnDhg0ZP3681zntpbXwGzRocPKbCtjFc97iO/3009m9ezf169dn9+7dnHbaaQA0aNCAnTt3Fnl8gwYNfL4XwdI+/HiQlgbTp9uWvoj9OX26zsNXZbZjxw7WrVsHwIIFC+jYsSNnnXUWWVlZbN26FYB58+Zx2WWXlXjsRRddxM6dO1mwYAF9+/Yt03mrV6/Or7/+6vW+AwcOnEyCs2fP9vkcS5cu5ciRI+zbt481a9bQrl07n8cWJPe6dety8OBBXn75Za/HpaWl8cUXX5S4FBx/5ZVXsnLlSnJzc8nNzWXlypVceeWVJZ6nZ8+eJ2fwzJkzh169ep28fe7cuRhj+Pjjj6lRo0bIunNAE378SEuDrCzIz7c/NdmrcjjrrLOYOnUqLVq0IDc3lzvuuIMqVaowa9YsbrjhBlq1akVCQgJDhw71+vgbb7yRDh06UKtWrTKdt0ePHixZssTroO3o0aO5//77adOmDXl5eT6f49xzz6Vz585cfPHFjB07ljPOOMPnsTVr1uS2226jZcuWXHnllX4/HPypXbs2Y8eOpV27drRr145x48ZRu3ZtwH7jKZiyOmbMGN5++22aNm3KqlWrGDNmDADdunWjSZMmnHnmmdx2221Mmzbt5HNfeuml3HDDDbzzzjskJyezYsWKcsVYmJRnZDoc2rZta3QDFOUm33zzDS1atIjY+bOysujevTtff/11uZ+je/fu3HvvvXTp0sXByEo3fvx4n4O+8czb34yIfGqMaevteG3hK6WCtn//fpo1a0bVqlXDnuxV4HTQVikF2OmR5W3d16xZk++++87hiAI3fvz4iJ07lmgLX6koEq1drCr6lOdvRRO+UlGiSpUq7Nu3T5O+KlVBPfwqVaqU6XHapaNUlEhOTiY7O5u9e/dGOhQVAwp2vCoLTfhKRYmKFSuWafcipcpKu3SUUsolNOErpZRLaMJXSimXcCThi8gLIvI/EfE6iVesKSKyVUQ2isj5TpxXKaVU4Jxq4c8GrvJz/9VAU89lCPCMQ+dVgcjIgNRUSEiwP3UnLKVcyZGEb4x5H8jxc0gvYK6xPgZqikjoSsKp32Vk2NLJ27fbXbC2b7fXNekr5Trh6sNvAOwsdD3bc1sRIjJERDJFJFPnIjskPb1onXyw13X7Q6VcJ6oGbY0x040xbY0xbevVqxfpcOKDr20OdftDpVwnXAl/F9Cw0PVkz20q1Hxtc6jbHyrlOuFK+MuA/p7ZOhcDB4wxu8N0bnebONFud1iYbn+olCs5NS1zIbAOOEtEskVkkIgMFZGCbXGWA9uArcAM4E4nzhvTwjVzJpq3P9TZQ0qFle54FQkFM2eKbzoeLYk4HPQ9UCok/O14pQk/ElJT7fTI4lJS7H60bqDvgVIhoVscRptQzJyJte4RnT2kVNhpwo8Ep2fOxOLiKp09pFTYacKPBIdmzuTkwIoVsP//vC+uOjAsnVmzYPly+PRT+OmnION28luEzh5SKux0A5RIKBiUTE+3XRiNGtlE52ew8tAh+Owz+OST3y/ff2/vO4H3bpDq+3cwcGDR2xo1go4doUMH6NQJzj47wJiLD7IWfIso/HrKohzvgVIqODpoG8W++gr+/W/46CPYvBny8+3tDRtCu3a/XzrdkkpidskB0PyGKWStyeKnn2zr/ocfYN06+OAD2O1ZBdGxIwwfDr16QQV/H/++BlnBDrSGIllnZOgHglJl5G/QFmNMVF4uuOAC41YbNxrTu7cxYEy1asZcfbUx48YZ8/rrxuzZ4+UB8+cbk5RkH1BwSUqyt3uRn2/Mtm3GTJ5sTOPG9vCUFGMef9yY3FwfQYkUff7iFz/nK5cyviallAVkGh95NeKJ3dfFjQn/yy+Nuf56+69y6qnGjB1rzL59AT54/nybtUXszwATY16eMUuWGHPZZfa8p5xizLBhxnz3XbEDU1L8J/yCTw2n+Dqfk+dQKg75S/japRMFvvwSJkyAV1+FU0+1XSzDh0OtWuGN44sv4F//ggUL4NgxuOYaG0eXLiALvCyUKk7k936nYCUk2BQfynMoFYd0Hn6U+vZbuP56aN0aVq2CcePsmqOHHw5/sgcbx6xZtst8/Hg7MNy1K7RqBTMOp3Hs354SDb44OaVSp20q5ThN+BGyciVceGF0JPriTj8dHnrIJv7Zs6FiRdu4P2NUGulpWRyYNj/0Uyp12qZSzvPV1xPpSzz34T/7rDGJicace64xO3b4ObCc/fJOy8835r33jLnuOhtKrVrGrBo43+Q3CnFsUfL6lYol6KBtdDhxwpgRI+y7fvXVxvzyi5+Do3SWyubNxnTqZMNp396YL76IaDhKqWL8JXzt0gmTw4ehd2944gn4v/+DZcugenU/D4jSrQlbtIB334U5c2DLFmjbFqZO9T6+qpSKLprww2DPHruq9bXX4Kmn4OmnS1nkBFFdXEwE+ve3g85XXQXDhsGtt8Jvv0U6MqWUP5rwQ+zrr+Gii2DTJpvw77nHJsxShXKWikM1cWrXhqVL7QDvnDlw6aVR8XmklPJBE34IrVgBl1wCeXmwdi307FmGB4dqlorDlTUTEuwUzqVLbRfPBRfA6tXBhaiUCg1N+CHy7LN24VKTJrB+PZx/fhmfIFRbE4ZobKBnT9iwAerVs3P3J0/Wfn2loo2utHXYiRPwt7/Zwdlu3WDRolIGZ8MtxCtYf/0VBgyAJUugb194/vmSX1SUUqGjK23D5NCh32fiDBtmuzmiKtlDyFewVq8OL79se54WLbJdWj/84MhTK6WCpAnfISdO2G6NpUttPZqAZuJEQhhWsCYkwAMPwJtv2iGCtm3tymKlVGRpwnfII4/Y+ekzZsDdd0c6Gj9CNTbgxdVXQ2YmNGhgf3/0Ue3XVyqStA/fAWvX2nn2ffvCvHkBTrt0kUOHYNAgePFF2+U1axZUqxbpqJSKT9qHH0I5ObZx3LgxPPOMJntvTjkFFi6Exx6zJaAvvthO4VRKhZcm/CAYY1uue/ZE4WycKCMCI0fatQl79titGd98M9JRKeUumvCDMG2aXT07aZIdmFSlu+IK26/fuDH06GE3fgloNqhDq4OVcjNN+OX05ZcwYoQdjBw+PNLRxJbUVPjwQ9sV9tBDdt5+Xp6fBzi8Olgpt9KEXw6HDkGfPnazktmzbaNTlU1SEsydC3//O8yfDzffDMeP+zg4SiuHKhVronGmeNS75x5bKfLtt+G00yIdTewSgQcfhCpVYNQom/AXLoRKlYodGMWVQ5WKJdo2LaNFi2DmTLj/fru5twreyJF2sdqrr9ppm0ePFjtA97dVyhGa8Mtg2zbbddy+va0QqZxz9912EPz11+Haa4vV1tf9bZVyhCb8AB07ZhdWJSbCggV2Y2/lrDvusMXWVqywM3gOHfLcEejqYJ3Jo5RfjiR8EblKRL4Vka0iMsbL/beIyF4R+cJzGezEecPpwQdt+d/nn7e5RIXGoEF2M5XVq2210V9/9dyRlgZZWXYOZ1aW92SvM3mU8ivohC8iicBU4GrgbKCviJzt5dAXjTGtPZfngz1vOK1YYVeJDh0K118f6WjiX79+Nk9/+KHdQvGXXwJ4kM7kUapUTrTwLwS2GmO2GWOOAYuAXg48b1TYs8fu39qypd3UQ4VHnz629s6GDbalf7J7xxedyaNUqZxI+A2AnYWuZ3tuK+56EdkoIi+LSENvTyQiQ0QkU0Qy9+7d60BowcnPt63NX3+1s3OqVo10RO5y/fV2vGTdOjuQe+SIn4N1Jo9SpQrXoO3rQKox5lzgbWCOt4OMMdONMW2NMW3r1asXptB8e+wxWLXKThk855xIR+NON9xgq2uuWgU33uhncZbO5FGqVE4k/F1A4RZ7sue2k4wx+4wxBbOrnwcucOC8zvAxs2P9ejtQe8MNMDiQIWadIVJ2Ab5n/fv/PmUzLc1uNlNCGOv8KxWzjDFBXbCrdbcBjYFKwJfAOcWOqV/o9+uAj0t73gsuuMCE3Pz5xiQlGWPnddhLUpI5OH2+SU01JiXFmNzc8j+PmT8/1K8gdpXjPXv8cXvYkCHG5OeHMValYgiQaXzla193lOUCdAO+A74H0j23TQB6en7/B7DJ82GwGmhe2nOGJeGnpBRNOJ7L/5JSTGKiMevWBfc8JiUldLHHunK+Zw88YA978MGwRKlUzPGX8N2941VCgtc99/IRJk3M54EHgnseRAKs/etC5XzPjLHT659/HqZMgbvuCmGMSsUg3fHKFx8zOHZXaMSoUcE/j84Q8aOc75mI3Vns2mttEbtFi0IQm1Jxyt0J38vMjiNUpF7SQSpWLsPgq84QKbsg3rMKFWxVzT/9yQ7orlwZohiVijPuTviFZnYYEX6mDokiVPplX9mW5+sMkbIL8j2rUgWWLoWzz4a//MUu0AqYzqhSLuXuPvxCbr8dHpieSgrbS96ZkmLrt6ios2cPdOgABw7ABx9A8+alPKCg5k7hMgxJSfoBreKGvz58TfjY7QrPPx/y8hMQdPA11nz/PVxyiW31r1sHZ5zh5+DUVPvNrTj9UFdxQgdtSzFqFNSsCfnJOvgai/74R3jrLdi3D665plCFTW+05o5yMdcn/Pfes1sVpqdD4iQdfI1V558PixfDV1/Z1dE+SzDojCrlYq5O+MbA2LFQv77dfEMHX2OUZxD26msSyK2RSu0VGdx5p/dp/jqjSrmZqxP+22/D2rW2dX+yEmZpG22o6FJs45PqOduZXWEIh57P4NFHvRyvH+rKxVw7aGsMXHyxneXx3XdQuXLITqVCyccg7N6kFE47nMXChba2vlJuoYO2Xrzxhp27PW6cJvuY5mOwte5vO3j4rAza903F6Hx7pQCXtvDz8+0g36FDsHmzbkge03xNs6xTB3P4N+Q3nW+v3EVb+MW8+qqde//QQ5rsY56vQVgomuxB97hVrue6hH/ihO3GadEC+vaNdDQqaL4GYXNyvB5udL69cjHXJfyFC+Gbb2DCBEhMjHQ0yhHeZlb5mFefc0oj79M1fdG6OyqOuCrhHz8ODz8MrVvbglsqjnnp6jlWIYm7Dk7k6acDfI5iUz4DLqanVJRyVcJfvBi2brVJP8FVr9yFvHT1VJg1nd+uTePee2H5ckpvvaenFy2yBjoOoGKaa2bpGANt29r/r5s2acJ3q0OH4NJLofXmDJ6XISQc8TOLR3cyUzFIZ+lga+Z89hncd58mezc75RRYtgwm5KUXTfZQsvWudXdUnHFN6nviCahXD/r1i3QkKtKSk6FBfgBVM7Xujoozrkj4//2vXVn7zKUZVGmeqjMuFBJI613r7qg444qE/+ST0D8xg7+8pTMulMfEiSVX3VWsWLL1rsX0VByJ+4S/dy/MnQuTq6bryktVlEiRq/lGfByoVHyI+4Q/bRocOQK1D+lOR6qQ9HQ4dqzITQl5xzg+WhsAKoJCvNAvrhP+b7/B1Kl227uA+myVe/j4oE/8cQe5uWGORSkIy0K/uE748+fbLp0RI9AZF6ooHx/0O2lEnz6QlxfmeJQKw0K/uE34+fkweTK0aQOdOqEzLlRRPhoA2wZPZOVK+NvfIhOWcjFf3csOdjtXcOyZosxbb9npmPPnFxqbS0vTBK+sgr+D9HT7H6pRI5g4kc5padxV1TYWzj0XBgyIbJjKRRo18r63g4PdznFbWuHyy2HLFti2TWveq7LJy4OrrrL7Hb/3nt0KU6mQK+jDPxzcpj2uK63w+eewejXcfbcme1V2FSrAiy/aFbnXXQe7dkU6IuUKYeh2jssW/s03w9KlsHMn1KzpcGDKNTZtsq37Fi1sS79q1UhHpNxg/36oUaPEMpGAuaqFn51tW2eDB2uyV8E55xw7BvTJJ/abdpS2jVScGTwYOnYMzd+bIwlfRK4SkW9FZKuIjPFyf2URedFz/3oRSXXivN7UqWPn3t9zT6jOoNykVy/4+99t4p88Gd0BS4XUnj22d6J9+/K38P0JepaOiCQCU4GuQDbwiYgsM8ZsLnTYICDXGHOmiPQBHgVuCvbc3lStaltjSjklPd1uev/ZyAxOVB5C4lHPoFrBwhjQ2V/KEdOm2X23Q5XDnGjhXwhsNcZsM8YcAxYBvYod0wuY4/n9ZaCLSCg+v5RyngjMng2PVUz/PdkX0HpMyiGHD9uE36MHNGsWmnM4kfAbADsLXc/23Ob1GGNMHnAAqFP8iURkiIhkikjm3r17HQhNKWeccgrUz9N6TCp05syBfftg5MjQnSOqBm2NMdONMW2NMW3r1asX6XCUKkLrMakSHBrTyc+3ZdzbtbMDtqHiRMLfBTQsdD3Zc5vXY0SkAlAD2OfAuZUKHy/lGIzWY3IvB4udvf66XSg6cmRoBmsLOJHwPwGaikhjEakE9AGWFTtmGVCwSL038K6J1gUASvlSaGGMQcgiheW9tB6TazlY7Ozxx+06q7/8xaHYfAg64Xv65IcBK4BvgJeMMZtEZIKI9PQcNhOoIyJbgfuAElM3lYq4QL6ee3bAMifyGdU7i+4L03jllXAHqqKCQ8XO1q+HDz6Ae++1q7xDKS5X2ipVZuWoY/Lbb9Cly++lPLTmjsukpnovdpaSYrfDDNCNN8LKlbYyQPXqwYflqpW2SpVLOb6eV61qF8k0aAA9e8L334c4RhVdHNhj44cf4JVXYOhQZ5J9aTThKwXl/nperx4sX24Xy3TrBjk5IYhNRScHip099ZTtQbzrrhDGWYgmfKXA99TKAKZcNmtmW/pZWXDttXD0qLOhqSjmGdMhP9/+LEOyz82FmTPhr3+13xLDQRO+UhD01/OOHe3CmbVr4dZb7f9/pfx57jk4dAjuuy9859SErxQ48vW8Tx/4xz9g4UIYNy6EsaqYd+wYTJkCXbvCeeeF77xxu8WhUmXmwBaYf/ub3WVt4kRo3BgGDXIoNhVXFi6E3bth1qzwnlcTvlIOErHluXfsgNtvh4YN4c9/jnRUKpoYA088AS1bhv9vQ7t0lHJYxYrw0kt2A5XeveGrryIdkYomb79t/yZGjAhtGQVvNOErFQKnngpvvmnnVnfrBj/+GOmIVLR44gmoXx/69g3/uTXhKxUiyck26e/fD9dcA7/+GumIVKRt3GhX1d51F1SuHP7za8JXKoRat7bdO199ZWfx5OVFOiIVtCBKIk+ebPdWuP32kEXnlyZ8pULs6qvtQO7y5bZlF6Xlq1QggiiJ/OOPsGABDBwItWv7ef4Q7pmsCV+pMLj9dhg9Gp591vbhqhgVREnkp5+2JTiGD/dxgIP19X3RaplKhUl+vu3WWbzYdvPccEOkI1JllpDg/SuaiN/l1QcP2im6V1xh//29cqj6plbLVCoKJCTY8guXXAL9+sG6dZGOSJVZOWsuvfCCHbwfMcLPQQ7V1/dHE75SYVRQUrlhQy2pHJPKUXMpL89WxezQoZQ9E8KwZ7ImfKXCrG5dO4BrjJ2jv093d44d5ai5tGSJrXvvt3UPjtTXL4324SsVIR9+aHfMatfOrr6sUiXSESmnGWNb9Tk58N//QmJiKQ/IyLADwDt22Jb9xIllru+kffhKRaEOHWyf/gcfaEnleLV8OWzYACNHBpDsIaj6+oHQhK9UBN10E0yaBIsWwdixkY5GlRDEvPj8fNtY/+Mf7dz7aKDVMpWKsNGj7eDtI4/YksqDB0c6Ipcq3p3SrZv9ClYw775gXjwE1PJ+8UX48kv7tBUrhjDuMtA+fKWiwPHj0KMHrFpluwG0pHKYFSx6KryoSsT7nPsA5sUfPw4tWtgyCp9/br8ghIv24SsV5YqXVN64MdIRuYy3FbS+GsMBzIufPv33b23hTPaliaJQlHK3wiWVr7kGdu2KdEQuUpbFTaXMi//lF3j4YejUCbrtD21tnLLShK9UFClcUrl7dy2pHDa+knjxHUoCmBf/6KOwdy/M7JKBhLg2TlnFX8IPcbU5pUKtcEnlm27Sksph4WvR09ChZVpktW2bLY6XlgZNni9/obWQMcZE5eWCCy4wZTZ/vjFJScbYz1N7SUqytysVY557zv4JDx1qTH5+pKNxgfnzjUlJMUbE/gw0bxR63E9VU8wtleab7Gxjn6dwLiq4iITuNRhjgEzjI6/G1ywdh6rNKRUtxoyxXQT//CeMGhXpaFQJXmb3HKuYRKVZ0+Gee7zXzahTB37+OWQh+ZulE1/z8MNQbU6pcHrkEVuHZfRo257RkspRxsvsnkrHI9xt40d89eGHodqcUuFUvKTyRx9FOiJVhL9GZk6O9/t83R4G8ZXww1BtTqlwq1KlaEnlrVsjHZE6yV8jMwoboPGV8MtRulSpWFBQUhm0pHJUmTiRo4k+GplR2AANKuGLSG0ReVtEtnh+1vJx3AkR+cJzWRbMOUsV4mpzSkVK06a2pb9jB1x7LRw5EumI1Js107j1xHT21/A0MuvUsbvc9Otn+/EHDIiqBmiwLfwxwDvGmKbAO57r3vxmjGntufQM8pxKuZaWVI4ev/xip+lvPCeNqj9lwbx58Ntv9utXwUKrOXNsiz5KGqDBJvxewBzP73OAa4N8PqVUKQqXVH7wwUhH415jxtjyFzNnQuXKeK/HE+mFVsUEm/BPN8bs9vy+Bzjdx3FVRCRTRD4WEZ8fCiIyxHNc5t69e4MMTan4NXq0nf79j3/AjBmRjsZ93nsPnnkGhg+Hiy7y3Ohrxo63tUERUurCKxFZBfzBy13pwBxjTM1Cx+YaY0r044tIA2PMLhFpArwLdDHG+N2+WcsjK+VfXp4tqfz223b9z003RToid8jNhfPOg0qVbL37U07x3OFr4aeI7e4JU3dOUOWRjTFXGGNaerksBX4Skfqek9QH/ufjOXZ5fm4D1gBtyvlalFIeFSrA4sV2jv5f/woLF0Y6ovhnDNx2G+zebbvUTiZ7sH31xYutFTwoSrp1gu3SWQYM8Pw+AFha/AARqSUilT2/1wU6AJuDPK9SCqhWzU7X7NgRbr4ZFiyIdETxbcYMeOUVuwK6bfE2dFpaUDX0wyHYhD8J6CoiW4ArPNcRkbYi8rznmBZApoh8CawGJhljNOEr5ZCCpP+nP9nZgPPmRTCNg5QAAAviSURBVDqi+LR5s+2z79oVRozwcVBKivfbo2S1f1AJ3xizzxjTxRjT1NP1k+O5PdMYM9jz+0fGmFbGmPM8P2c6EbhS6nennGLr6HfqZKd+u3IgN4Sl0Y8cgT597Ifr3Ll+drGKwsVWhcXXSlulXCwpCd54A66+2s7gefxx3z0McaegamWINhsZPdruTzB7NvzB2xSWAlG+2j++yiMrpTh2zPbnL15sc96//233zI1rISyN/vrrtobR8OHw5JNBPVVY6CbmSrmBp0ujUpUEXlyfyqKeGUyfblv8ubmRDi7EQlQa/bvvbBdZ69Z2sVus04SvVDwo1qUhO7Zz06ohrBmSwfvvQ/v28L3flS8xLgSVKX/+2W4mn5hoZ+ZUrlzup4oamvCVigc+lvVftiKdVavsptoXXQTvvx+Z8ELO4cHSo0fhuutg505YtgyaNHEgxiigCV+peOCnS+NPf4L1622J5SuusPW84k5ZB0v9zOgxBgYOtAXq5s61347ihq/NbiN9Kdcm5kq5VUqK9w2zU1JOHpKTY0yXLvbm++835sSJiEUbWXfcUXKD8aSkk5uWjxtnb3rkkQjHWU742cRcW/hKxYMAujRq1YK33vq96NqNN5bsBYp7GRnw7LMl56t6qlrOnQsTJtgW/hhfxd5jmCZ8peJBgF0aFSvafPfEE/Dqq3DZZbYujGukp/tcnGB27GDwYLj8cvseeSuLE+s04SsVq4r3Q0NAu72JwH332d2zvvkGLrwQvvgiTDGHW/H3yE+p4h2mEa1a2Rk58bpuQRO+UrHIgZWlPXrAhx/a3zt2tLNR4oq398hHsz0f4fnGE1m1CmrW9HpIXNCEr1Qscmh3pfPOgw0b4Oyz7T65Dz0UR3vlenuPjCmR9PMRXj1tKCM/T6OW112544cmfKVikYMrS+vXhzVrbA/QhAlwzjm2nECUVl0JnK/3whhIScEgbCeFCU3n8ect06hRI7zhRYImfKVikcMrS5OSbFnlVavsitKePaF7d9iyJYgYI83He2EapfD3QVkkkM/Ay7MY8Wkap54a5tgiRBO+UrEoRGV4u3Sx2/ZNngxr10LLlvDAA3DwYFBPGxle3iNTNYmnTpvIuHHQv7+dplq9eoTiiwBN+ErFohCW4a1YEe691xYO69PHztlv3hxefDHGunmKvUfH6qcw4tTpjPo8jSeftKWOK1WKdJDhpeWRlVJ+ffQRDBsGn39uN1h5+mnb8o8ZGRkcvCedpH072JXQiEPpE2k+ITrq04eClkdWSpXbJZfAJ5/AM8/Axo22VPDw4bB/f6QjK92eyRkcGTCEavu2k4ChYf52mj/h3MYosUYTvlKqVImJMHSo7ea57TaYMgXOOgtmzbLrvKLN4cMwdiwcHZFOlRPBT1+NF5rwlVIBq1PHtvQzM+GPf7Q1Zy65xF6PBsbYkhEtWsD/+3/QiNBsjBKrNOErFS9CuIl3ceefb8sHz5ljqzhceCH07Wv30V2yxHb9hHNmz4EDtpTx5ZfD9ddDjRrw3nsgKc5vjBLLKkQ6AKWUAwrKCBSsLC0otQAh20A7IcFObezVyy7Ymj0bFi0qeszpp9tvAk2a2J+FL6edFlyBskOH7KbtixbZ6ZVHj9o8/q9/wZ13QoUK2KmZhd8XcGT6aqzSWTpKxYMQbuJdFvv3260UvV2ys4tO66xWzfsHQZMmNuwKXpqjR47Af/5jk/zrr9s8Xr++LfXcp4/d1avEh0hGhu2z37HDfiJMnBiyD8Fo4G+WjiZ8peJBQoL3SfIiUTOqevSo/ezx9mGwbZu9v0Biok36hT8Evv4aXnsNfvnF7t7Vu7dN8h072uOV5S/ha5eOUvGgUSPvLfwo6quuXNnO7DnrrJL35efDjz8W/QAo+P2llyAnx1ax7N0bbrrJ9tV7+wag/NO3TKl4EON91QkJkJxsL5ddVvL+AwegalX3rYx1ms7SUSoeFC4jALaPo2C+eRwsMqpRQ5O9E7SFr1S8KBiIDPNsHRU7tIWvVDxxaGMUFZ804SsVTxzcGEXFH034SsUThzdGUfFFE75S8SREG6Oo+BBUwheRG0Rkk4jki4jXif6e464SkW9FZKuIjAnmnEopP0K4MYqKfcHO0vka+AvwnK8DRCQRmAp0BbKBT0RkmTFmc5DnVkp5k5amCV55FVTCN8Z8AyD+KyBdCGw1xmzzHLsI6AVowldKqTAKRx9+A2BnoevZnttKEJEhIpIpIpl79+4NQ2hKKeUepbbwRWQV8Acvd6UbY5Y6GYwxZjowHWzxNCefWyml3K7UhG+MuSLIc+wCGha6nuy5TSmlVBiFo0vnE6CpiDQWkUpAH2BZGM6rlFKqkGCnZV4nItlAe+BNEVnhuf0MEVkOYIzJA4YBK4BvgJeMMZuCC1sppVRZBTtLZwmwxMvtPwLdCl1fDiwP5lxKKaWCoyttlVLKJTThK6WCk5Fh99RNSLA/46D+frzShK+UKn/Szsiw9fa3b7d76hbU39ekH5U04SvldsEkba2/H1M04SvldsEk7XDU39cuI8dowlfK7YJJ2qGuv69dRo7ShK+U2wWTtENdf1+7jBylCV8ptwsmaYe6/r5u2egoTfhKuV2wSTstDbKyID/f/nSyFr9u2egoTfhKqdAm7WDolo2O0oSvlIpeumWjo4Ld4lAppUJLt2x0jLbwlVLKJTThK6WUS2jCV0opl9CEr5RSLqEJXymlXEKMMZGOwSsR2QtsD/DwusDPIQwnWunrdhd93e5S3tedYoyp5+2OqE34ZSEimcaYtpGOI9z0dbuLvm53CcXr1i4dpZRyCU34SinlEvGS8KdHOoAI0dftLvq63cXx1x0XffhKKaVKFy8tfKWUUqXQhK+UUi4RUwlfRK4SkW9FZKuIjPFyf2URedFz/3oRSQ1/lM4L4HXfJyKbRWSjiLwjIimRiNNppb3uQsddLyJGRGJ+6l4gr1lEbvT8e28SkQXhjjEUAvgbbyQiq0Xkc8/febdIxOk0EXlBRP4nIl/7uF9EZIrnfdkoIucHdUJjTExcgETge6AJUAn4Eji72DF3As96fu8DvBjpuMP0ujsDSZ7f73DL6/YcVx14H/gYaBvpuMPwb90U+Byo5bl+WqTjDtPrng7c4fn9bCAr0nE79Nr/BJwPfO3j/m7AW4AAFwPrgzlfLLXwLwS2GmO2GWOOAYuAXsWO6QXM8fz+MtBFRCSMMYZCqa/bGLPaGFOw0/PHQHKYYwyFQP69Af4OPAocCWdwIRLIa74NmGqMyQUwxvwvzDGGQiCv2wCnen6vAfwYxvhCxhjzPpDj55BewFxjfQzUFJH65T1fLCX8BsDOQtezPbd5PcYYkwccAOqEJbrQCeR1FzYI2yKIdaW+bs/X24bGmDfDGVgIBfJv3QxoJiIfisjHInJV2KILnUBe93jgZhHJBpYDd4UntIgr6/9/v3THqzgiIjcDbYHLIh1LqIlIAjAZuCXCoYRbBWy3TifsN7n3RaSVMWZ/RKMKvb7AbGPMEyLSHpgnIi2NMfmRDiyWxFILfxfQsND1ZM9tXo8RkQrYr377whJd6ATyuhGRK4B0oKcx5miYYgul0l53daAlsEZEsrD9m8tifOA2kH/rbGCZMea4MeYH4DvsB0AsC+R1DwJeAjDGrAOqYIuLxbuA/v8HKpYS/idAUxFpLCKVsIOyy4odswwY4Pm9N/Cu8Yx8xLBSX7eItAGewyb7eOjThVJetzHmgDGmrjEm1RiTih276GmMyYxMuI4I5G/8NWzrHhGpi+3i2RbOIEMgkNe9A+gCICItsAl/b1ijjIxlQH/PbJ2LgQPGmN3lfbKY6dIxxuSJyDBgBXZU/wVjzCYRmQBkGmOWATOxX/W2YgdC+kQuYmcE+LofA6oBiz1j1DuMMT0jFrQDAnzdcSXA17wC+LOIbAZOAKOMMTH9LTbA1z0CmCEi92IHcG+Jg8YcIrIQ+wFe1zM+8RBQEcAY8yx2vKIbsBU4DNwa1Pni4D1TSikVgFjq0lFKKRUETfhKKeUSmvCVUsolNOErpZRLaMJXSimX0ISvlFIuoQlfKaVc4v8DXSK+YuPJSqYAAAAASUVORK5CYII=\n" }, "metadata": { "needs_background": "light" } }, { "output_type": "display_data", "data": { "text/plain": [ "
" ], "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXwAAAD4CAYAAADvsV2wAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nO3deXxV1bn/8c8TBiFCZbZIIEENg4qCBGRqBbkUqwLWqpVfVBwRlXq1t1osV0QqFqdqaR2AolCJOF6FqnUWxQE0ilIckBTDJCoSiiKiYNbvj3UCSThJTjjz2d/363VeOcPO2WsfwnPWXuvZzzLnHCIikvmykt0AERFJDAV8EZGAUMAXEQkIBXwRkYBQwBcRCYiGyW5ATdq0aePy8vKS3QwRkbTy9ttvf+mcaxvutZQN+Hl5eRQXFye7GSIiacXM1tT0moZ0REQCQgFfRCQgFPBFRAIiZcfwRcTbuXMn69evZ8eOHcluiqSQJk2akJOTQ6NGjSL+HQV8kRS3fv16mjdvTl5eHmaW7OZICnDOsXnzZtavX0/nzp0j/j0N6WSCoiLIy4OsLP+zqCjZLZIY2rFjB61bt1awl93MjNatW9f7rE89/HRXVARjx8L27f7xmjX+MUBhYfLaJTGlYC/V7cvfhHr46W7ixD3BvsL27f55EZFKFPDT3dq14Z9fU+O1FyJxN3jw4HpdOPnwww/TvXt3hgwZQnFxMZdddhkAixYt4vXXX6/XvhctWsRJJ50U9TbRmjt3Lvn5+eTn5zN37tyw25SVlTFs2DDy8/MZNmwYW7ZsAeCjjz6if//+7Lffftxyyy0xa5MCfrrr1Cn882Yay5e0MXv2bGbNmsVLL71EQUEB06dPB/Yt4KeCsrIyrrvuOpYuXcqbb77JddddtzuYVzZt2jSGDh3KqlWrGDp0KNOmTQOgVatWTJ8+nd/+9rcxbZcCfrqbOtUH9+qc07COxERpaSndunWjsLCQ7t27c+qpp7I9NIz4wgsv0KtXL3r06MF5553Hd999V+V377nnHi6//PLdj2fNmsUVV1xRZZspU6bw6quvcv7553PllVfu7n2XlpZy9913c9ttt9GzZ08WL15c5ffefPNN+vfvT69evRgwYAArV67cq+2TJ0/mrLPOon///uTn5zNr1qzdr23bto1TTz1197FVrP43ZcoU+vTpwxFHHMHYsWPZl1UBn3nmGYYNG0arVq1o2bIlw4YN4+mnn95ruwULFjBmzBgAxowZw+OPPw5Au3bt6NOnT71SLiOhSdt0V1gIZ54Z/rWahnskbV1+Obz7bmzfs2dPuP322rdZuXIls2fPZuDAgZx33nnceeedjB8/nnPOOYcXXniBLl26cPbZZ3PXXXdVCfCnn346U6dO5eabb6ZRo0bce++9zJgxo8p7T5o0iRdffJFbbrmFgoICFi1aBPh6WuPGjaNZs2Zhe7rdunVj8eLFNGzYkOeff57f//73PProo3ttt3z5cpYsWcI333xDr169OPHEEwFYtmwZ77//PgcddBADBw7ktddeY9CgQYwfP55JkyYBcNZZZ/HEE08wYsSIKu9ZVFTEzTffvNe+Dj30UB555BE2bNhAx44ddz+fk5PDhg0b9tr+888/p3379gD8+Mc/5vPPPw/7+ceKAn4myM0NP2Zf03CPSD117NiRgQMHAnDmmWcyffp0hg0bRufOnenSpQvge6h33HFHlYDfrFkzjjvuOJ544gm6d+/Ozp076dGjR0zatHXrVsaMGcOqVaswM3bu3Bl2u1GjRtG0aVOaNm3KkCFDePPNN2nRogV9+/YlJycHgJ49e1JaWsqgQYN46aWXuOmmm9i+fTtlZWUcfvjhewX8wsJCCmOcBWdmcc/GUsDPBFOnVk3NBMjO9s9LRqmrJx4v1QNRfQLTBRdcwA033EC3bt0499xzY9ama665hiFDhvDYY49RWlrK4MGDw25XU9v322+/3c81aNCAXbt2sWPHDi655BKKi4vp2LEjkydPDpvrXlcPv0OHDrvPVMBfPBeufQceeCAbN26kffv2bNy4kXbt2kVw5PtOY/iZoLAQZs70PX0z/3PmTOXhS8ysXbuWN954A4D777+fQYMG0bVrV0pLSykpKQHgvvvu49hjj93rd4855hjWrVvH/fffz+jRo+u13+bNm/P111+HfW3r1q106NABgDlz5tT4HgsWLGDHjh1s3ryZRYsW0adPnxq3rQjubdq0Ydu2bTzyyCNhtyssLOTdd9/d61ax/fDhw3n22WfZsmULW7Zs4dlnn2X48OF7vc/IkSN3Z/DMnTuXUaNG1di2WFDAzxSFhVBaCuXl/qeCvcRQ165dueOOO+jevTtbtmzh4osvpkmTJtx7772cdtpp9OjRg6ysLMaNGxf2908//XQGDhxIy5Yt67XfESNG8Nhjj4WdtL3qqqu4+uqr6dWrF7t27arxPY488kiGDBlCv379uOaaazjooINq3LZFixZceOGFHHHEEQwfPrzWL4fatGrVimuuuYY+ffrQp08fJk2aRKtWrQB/xlORsjphwgSee+458vPzef7555kwYQIAn332GTk5OfzpT3/i+uuvJycnh6+++mqf2lKZ7csMdCIUFBQ4LYAiAh9++CHdu3dP2v5LS0s56aSTWLFixT6/x0knncQVV1zB0KFDY9iyuk2ePLnGSd9MEO5vw8zeds4VhNtePXwRiZv//Oc/dOnShaZNmyY82MveNGkrIrXKy8vb5959ixYt+Pjjj2PcoshNnjw5aftORerhi6SBVB16leTZl78JBXyRFNekSRM2b96soC+7VdTDb9KkSb1+T0M6IikuJyeH9evXs2nTpmQ3RVJIxYpX9aGAL5LiGjVqVK9VjURqoiEdEZGAUMAXEQkIBXwRkYCIScA3s3vM7AszC5usa950Mysxs+VmdnQs9isiIpGLVQ9/DnB8La//HMgP3cYCd8VovxKJoiLIy4OsLP9TK2GJBFJMAr5z7hWgrJZNRgF/d94SoIWZtY/FvqUORUW+dPKaNX4VrDVr/GMFfZHASdQYfgdgXaXH60PPVWFmY82s2MyKlXMcIxMnVq2TD/6xlj8UCZyUmrR1zs10zhU45wratm2b7OZkhpqWOdTyhyKBk6iAvwHoWOlxTug5ibealjnU8ocigZOogL8QODuUrdMP2Oqc25igfQfb1Kl+ucPKtPyhSCDFKi1zPvAG0NXM1pvZ+WY2zswqlr95ClgNlACzgEtisd+0lqjMmVRe/lDZQyIJpRWvkqEic6b6ouOpEogTQZ+BSFzUtuKVAn4y5OX59MjqcnP9erRBoM9AJC60xGGqiUfmTLoNjyh7SCThFPCTIdaZM+l4cZWyh0QSTgE/GWKdOZOoi6tieRah7CGRhFPAT4ZYZ84kYngk1mcRqZw9JJKhNGmbCRIxAVrTPir2M3Vq7IN1UZE/S1m71g/1xGMfIhlGk7aZLhHDI7WdLcRjziAd5yVEUpwCfiZIxPBIXZOpsZ4zUNE3kZjTkI5EJtyFUtWZQXl5bPaXleV79vHch0gG0pCORK/yWURNYplSqbRNkZhTwJfIFRb6SeB58+I/Z6C0TZGYU8BPZal69Wwi5gyUtikScxrDT1UqLiYi+0Bj+OlIWSoiEmMK+KlKxcVEJMYU8FNVPLNUUnVuQETiqmGyG5ApvvsOli+HjRvhhx98qvgPP9R9Pzvbx9zOnSEnBxpW/ItMnRp+DD/aLJXqcwMVV7CC5gZEMpwmbffBzp3wwQfw1ltQXOxvy5f756PRoIEP+nl5/jbi6yKGvzyR/cvWQsdO2A0xqCWjhUdEMppWvIpCeTmsXOmDekWAX7YMduzwrx9wABQU+FufPj6eNmjgb1lZdd//6isffz/5xMfbitsnn8Cnn+652LRNGxg0yN9+8hPo1QsaNdqHA9IVrCIZrbaAryGdar75Bp5+Gt54wwf3t9+Gbdv8a/vvD0cfDRdfvCfAH3KIj6H7qm1b/x7hfPcdrF4Nr78Or74KixfD44/717KzoV+/PV8A/fpBs2YR7LBTp/A9fF3BKpLx1MPHD8U88wzMnw8LFvig37gx9Ozpg3pFcO/WzffKk2njRh/8K74A3nvPd8wbNvTBf9QoOPVUPzQUlvL7RTKahnRq8NlnMGMG3H23v9+qlQ+WZ5wBAwf6oJ/qvvrKn40sWgRPPAErVvjRmZ/8BEaP9sfTpk21X1KdeZGMpYBfiXOwdCn85S/w8MO+d3/CCTBuHAwfnh5BvjarVvkzlfnz4aOPfM//Zz/zwX/UKGjePNktFJF4UsDHB/qnn4Zrr/WTrz/6EZx7Llx6KeTnx2w3KcM5P9wzfz488IDvzDdtCr/6lR/R6dfPnwmISGYJfGmFpUthyBDfk//yS/jrX2H9erj99swM9uCDec+ecOONPuPn1VfhrLPgkUdgwAA46ii4447ay9uLSGbJ6ID/0Ufwy1/63uyHH/phnI8+8r36IA1tZGX5OYkZM3yq54wZfuhq/Hiffn/DDfCf/yS7lXXQ1cEiUcvIgL9hA1x4IRx+ODz7LFx3HZSU+ACX7mP00Wre3A/pFBf7LJ8+ffz8bW4u/P738MUXyW5hGFrfViQmMm4M/+OP/XDFDz/4fPmJE6Fduzg0MIMsWwZ//KMf7tlvP/9l+dvfplBqvq4OFolYoMbw8/N9kP/oI/jznxXsI9GrFzz0kB/2Gj0a7rrLXwx23nn+CzTpVDlUJCYyLuCbwf/+Lxx8cLJbkn66doV77oF//9ufHc2f7y82+9Wv4N13k9gwrW8rEhMZF/Alep06wfTpfhRlwgSfztqrF5x4Irz2WhIapPVtRWJCAV9q1K6dz+BZswauvx7efNOXbxg2DN55J4ENiXR9W2XyiNQqJgHfzI43s5VmVmJmE8K8fo6ZbTKzd0O3C2KxX0mMFi38vEhpKdx6q5/k7d3bx9tPPklQIwoLfQPKy/3PcMFemTwitYo64JtZA+AO4OfAYcBoMzsszKYPOud6hm5/i3a/knj77w+/+Y0f47/6avi///Nj/JMmwbffJrlxWgNYpE6x6OH3BUqcc6udc98DDwCjYvC+kqIOOMAP9axaBaedBn/4g7/m4YknktgoZfKI1CkWAb8DsK7S4/Wh56r7pZktN7NHzKxjuDcys7FmVmxmxZs2bYpB0ySecnJg3jx48UVo0gRGjICTTw6fMh93yuQRqVOiJm3/AeQ5544EngPmhtvIOTfTOVfgnCto27Ztgpom0RoyxKdtTpsGzz0H3bv7C7m+/z6BjVAmj0idYhHwNwCVe+w5oed2c85tds59F3r4N6B3DPYbG7HK7Ah4hkjjxvC73/mLt44/3pdpOOoo3/uvUSw/s0gzeUSCzDkX1Q2/TOJqoDPQGHgPOLzaNu0r3f8FsKSu9+3du7eLu3nznMvOds7ndfhbdrZ/Phnvk0GefNK5gw/2H8Xo0c59+mm1DfSZicQFUOxqiKsxqaVjZicAtwMNgHucc1PNbEpoxwvN7I/ASGAXUAZc7Jz7qLb3TMgSh7Gq0aJaL2F9+60vzzxtmj8DmDLFF7Br2BB9ZiJxogVQapKV5fuW1Zn5fO9Ev0+GKimBX//aX7F75JG+Vs+AQfrMROIhUMXT6iVWmR3KEKnVoYfCU0/Bo49CWZmvzf9ltj4zkUQLdsAPl9nRqBFs21a/iURliNTJDE45xU/qXnUVXLF9KtvRZyaSSMEO+NUzO1q39j83b67f5fnKEIlYs2Z+XH/Cvwq5tdtMSsmlHGPnQQn8zAKeUSXBFewx/Oo0kZhQzsHcuXD55T5n/8Yb/fKTWfHshlTU3KlchiE7W1/QkjE0hh8pXZ6fUGZwzjnw/vsweDBcdhkcdxysXh3HnarmjgSYAn5lmnxNig4d4MknYfZsX4nzyCPhzjvjlKyjL3UJMAX8yjT5mjRmfknFFSt8Fs+ll/q6+zEfSdOXugSYAn5lmnxNuo4dfb7+rFnw1lvQo4e/X+tUU30mYfWlLkFW0yW4yb4lpLSCpLQ1a5wbOtRXXTjpJOc++yzMRvtSomHePOdyc50z8z9VzkEyCPEurRAPScnSkZRTXg5/+YtfW7d5c3/CdfLJlTZQZpVIFcrSkbSVlQX//d/w9tu+/v4vfgEXXuivjQNqn4RVvr1IFQr4khYOOwyWLPE9/dmz4eij/Rh/jZOtrVppjVuRahTwJW00buwXVnnpJdixAwYMgAV9p+LCTcKC8u1FqlHAl7Rz7LHw3nu+Ns/JDxfyh44z2dWhWmZVWVn4X1a+vQSYAr6kpZYt4YEHYM4cuHlDIW22lfLg/HI/UVtYGLt8e80DSAZRwJe0ZQZjxvj1dLt1gzPO8BdvbdtGbPLtK+ruaB5AMoQCvqS9Qw6BxYv98PycOdC7Nyw7LIKL6OrqvavujmQY5eFLRlm0yMf0L7+Em27yBdnMwmwYSdVMrWQmaUh5+BIYgwf7Cd3hw33Z5REjYNOmMBtG0ntX3R3JMAr4knHatIEFC2D6dHjuOejZ0/f8q4ikaqbq7kiGCVbAV8ZFYJj5hdOXLvWrbB13HFx7LfzwQ2iDSHrvKqYnGSY4AV8ZF4HUs6cvy3DmmTBlig/8Gzbge+mNGlXduFGjvXvvhYU+1bO8UsqnSJoKTsBXxkVgNWsGf/+7X07x7bfhqKP8Qit7zeaGnd0VyRzBCfha6Sjwzj7bB/wOHaDlrRP9QrqVff+9OgCSXHEedg5OwFfGhQBdu/px/VzUAZAUk4Bh5+AEfGVcSEiTJmC56gBIiknAsHNwAr4yLqSyMB2A7xpm8/1kdQAkSRIw7BycgA/KuJA9KnUAnBllzXM5d9dM+k0vpKQk2Y2TQErAsHOwAr5IZaEOgJWX0+qrUkYvLKS01C+u8sADyW6cBE4Chp0V8EVCRozwlTd79IDRo+Gii+Dbb5PdKgmMBAw7q3iaSDU7d8I118CNN/rg/9BDvvyySDpQ8TSRemjUCKZNg3/+EzZuhIICuO++ZLdKJHoxCfhmdryZrTSzEjObEOb1/czswdDrS80sLxb7FYmn44/3Qzy9e/uLts49F3bco3pMkr6iDvhm1gC4A/g5cBgw2swOq7bZ+cAW59yhwG3AjdHuVyQROnSAF17wQzzfzynCXaB6TJK+YtHD7wuUOOdWO+e+Bx4ARlXbZhQwN3T/EWComQqXSHpo2NAXXvtbu4k0darHJOkrFgG/A7Cu0uP1oefCbuOc2wVsBVpXfyMzG2tmxWZWvCnsqhUiydN0k8oxSHpLqUlb59xM51yBc66gbdu2yW6OSFU1XACzvY3KMQRWmq2xEYuAvwHoWOlxTui5sNuYWUPgAGBzDPYtkjhhLoz51rK5YNNULroIvvkmSe2S5EjDNTZiEfDfAvLNrLOZNQbOABZW22YhMCZ0/1TgRZeqFwCI1CTMhTEN751JzpWFzJrlF1tZujTZjZSEScM1NqIO+KEx+fHAM8CHwEPOuffNbIqZjQxtNhtobWYlwG+AvVI3RZIuktPzavWYGo0p5Kab4KWXfDn9gQNh8mTYtSuxTZckSMM1NnSlrQjsOT2v3GPLzq7Xpe1bt8L48TBvHvTt63/m58epvZJ8eXl+GKe63FzfKUgSXWkrUpcYnJ4fcIC/IvfBB2HVKj/EM3OmH96VDJSGa2wo4ItATE/PTz8d/vUvGDDAF2AbORI+/zzK9knqScM1NhTwRSDmtcg7dIBnnoE//xmee84XYfvHP6Jon6SmNFtjQwFfBOJyep6VBZddtmfh9JEj/TTBtm1RtlVkHyngi0BcT88PP9yna/7ud/C3v0GvXrBkSQzaLFJPCvgiFeJ4et64sS+5vGiRr7c/aBBce62/L5IoCvgiCfTTn8J77/nvkilTfN7+xx8nu1USFAr4Igl2wAEwdy48/DD8+99+iOfuu5W+KfGngC+SJKee6tM3Bw2Ciy/2a+p+9lmyWyWZTAFfJIkOOsgvpTh9ul9opUcPWLAg2a2STKWAL5JkWVnw61/DO+9Ax45w8slwwQXw9dfJbpmElWYlkStTwBdJEd27+3TNq6+Ge+/1pRneeCPZrZIq4l0SOc5fJgr4IimkcWO44QZ4+WWfHTpokF9PV+mbKSKeJZETUF9fAV8kBQ0a5NM3zz4brr/e1+VZuTLZrZK4lkROQH19BXyRFPWjH/mhnUcfhU8+8embd96p9M2kinHNpSoSUF9fAV8kxZ1yik/fPPZYuPRSOPFEpW8mTTxLIsfzyyREAV8kDbRvD089BX/9q19d64gj4LHHkt2qAIpnSeQE1NdXwBdJE2a+h79smU/gOOUUOO88pW8mXLxqLiWgvr6WOBRJQzt3+lo8N9zg48Jdd8Hw4clulaQCLXEokmEaNYI//AEWL/apnMcfD2ecARs3JrtlGSaNL7IKRwFfJI0NGODTN6dMgccfh27d/Dj/Dz8ku2VpqHpwv+SSuOfFh7N1K6xbF5/31pCOSIYoKfEx6rnnfArnX//qvxAkAhUXPVXOgzcLnwObm+vH7uOgpMSX1mjY0K+U1qBB/d9DQzoiAXDooX4d3Ycegi++8LX2TzvNBxGpQ7iLnmrqDMcwL76yBx+Eo4+GTz+FW2/dt2BfFwV8kQxi5oP8ypVw3XW+Eudhh8Hll8PmzcluXQqrTxCvT158BHMA334L48b5OZjDD/dZWEOHRr6L+lDAF8lA++8PkybBqlVwzjnwl7/AIYfALbfAjh3Jbl0KqimIm1V9XJ+8+Ahq46xcCf36wYwZcNVV8MorfsQoXjIv4GfYrLpINNq396nc773nx/OvvNJX5XzgAZVoqKKmi57Gjdv3vPg6auMUFUHv3rBhAzz5JNx4o8++iivnXEreevfu7ept3jznsrOd83/L/pad7Z8XEffss84ddZT/r9G3r3OvvJLsFqWQefOcy811zsz/jDRu1PR7ZlVjUehWbubOP98/HDTIuXXrYnsYQLGrIa4mPbDXdNungJ+bG/YDdrm59X8vkQy1a5dzc+Y416GD/+9x8snOrVyZ7Falqdo6ma1bh41HZQ1aO3Du6qud27kz9k2qLeBnVlpmVlb481Qzfxm0iOy2fTvcdhtMm+bH9ceN87X327VLdsvSSF6eH5uvLjcXtm0LO1O+mdYUP/1l3K6MDk5aZgKqzYlkiuxsP5xcUuKXVLzzTv9f5YIL4P33k926NFFbSeOysrAvtbKypJXByKyAn4BqcyKZ5sADfS2eDz7wGT333++rcQ4fDk8/rcndWtXWyazhNUtiBzSzAn4Cqs2JZKquXeHuu/1l/VOn+hr8P/+5D/6zZvl8cammhk7mliuncmfOVL4htTqgUY3hm1kr4EEgDygFTnfObQmz3Q/Av0IP1zrnRtb13iqtIJJc33/vr/687TZ/MVCbNn6c/9JL4cc/TnbrUkhRkR8bW7uW8pat2LEDmmwvYx2d2Nj7BPp+8RRZ69f6Hv/UqXHvgNY2hh9twL8JKHPOTTOzCUBL59zvwmy3zTnXrD7vrYAvkhqc84uq33Yb/OMfvs7L6NFwxRXQs2eyW5cadu6Ely8qYsC9Y8mmUu59dnbCRxniOWk7Cpgbuj8XODnK9xORFGMGgwfDggX+ytCLLvLr7PbqBccd578EgpoE9/33cM89vkrpofdOrBrsIeaLkEcr2oB/oHOuogL3Z8CBNWzXxMyKzWyJmdX4pWBmY0PbFW/atCnKpolIrOXn+zIN69b5K0NXrYKRI33Au/12n90ThEnerVvhT3+Cgw+G88+HFi0g12rI2AmXtpkkdQZ8M3vezFaEuY2qvF0o4b+mf+rc0CnG/wNuN7NDwm3knJvpnCtwzhW0bdu2vsciIgnSsqWv/bJ6Ncyf7x9fcYWf4G3XDn75S5g+Hd59N3N6/87BG2/4ZSUPOgj+53/8F+A//wnFxbVk35ilTImXaMfwVwKDnXMbzaw9sMg517WO35kDPOGce6S27TSGL5I+nPPB/+WXfQGwl1/eUzK+RQv4yU/gpz+FY4/1Q0ENGya1ufVSVgbz5vlMpRUroFkzP4cxdiwUVB4pLyqCs85KeA396uI5aXszsLnSpG0r59xV1bZpCWx3zn1nZm2AN4BRzrkPantvBXyR9LZ2bdUvgFWr/PPNm/ta/RVfAAUFfpnGVFJaCgsX+tvLL8OuXb6dY8f6MsbNm9fwi9Wra1Z+PkGnOvEM+K2Bh4BOwBp8WmaZmRUA45xzF5jZAGAGUI4fQrrdOTe7rvdWwBfJLJ9+6oN/xRfAB6EuX9OmvpJnxRfAMcdAkyaJbdvXX8PSpb5dCxfC8uX++e7d/RzFr37lz0zqVFuphXTv4ceTAr5IZtu0yS/C/vLL/rZ8uR8N2W8/H/QrvgD69/f1/etUKR++tpz3rVv92cby5fDOO/Daa/5+ebkvxzVokA/yI0f6Mfp6CbdUYoJTMxXwRSTlbdkCr7665wvgnXd8EG7Y0A+ndOniA3929t4/uy8ros+ssTT8bk+g3dkom4UjZrK4YyGbN/s5hlWr/BdNhWbN/JfLwIH+dswxcMABUR5IhF888aKALyJp56uv4PXX98wDbNgA33zjO8/V1xX5hDzy2HsopZRcjvpRKa1a+VGV/Pw9t8MP96uAxWPt2GRSwBcJgiT3LBOpvNyXdK74AujUOQsLE8ucGZYpeaERqi3gp1FylIjUqPrYccX6qZCRQT8ryw/l7K5b1qlT2MnSZFamTEWZVS1TJKjqWD8146k0ekQU8EUyQW0LcQRBfUujFxX5FMqsLP8zRa6EjTcN6YhkghqGNAK12lthYWTDV5dc4gv/V4z5Z/jwV2Xq4YtkAg1pRKaoqGqwrxCQ4S8FfJFMoNXeIjNxYs3lPAMw/KWAL5Kuqo9Dg798v7zc/1Sw3/szqq1UcQCGvzSGL5KOApaGuU/CfUZm4Xv4ZoEY/lIPXyQdBT0NMxLhPiPn9q5oaeYX6w3AF6UCvkg6CnoaZiRq+iycqzrXcd99cOediW1bkmhIRyQdKQ2zbjV9RgksVZxq1MMXSUdKw6ybPqO9KOCLpCOlYdZNn9FeNKQjkq4ivbI0yCo+n4oqohWT2gH93BTwRSRzKX21Cg3piEjmUvpqFQr4IpK5lL5ahQK+SKYIaMnfWtWUpj77Q3gAAAZbSURBVBrQ9FUFfJFMUDFWvWaNv7CoYqw66EFfqZlVKOCLZAKNVYen1MwqtIi5SCbIyqq5KFjAFvEOutoWMVcPXyQTaKxaIqCAL5IJNFYtEVDAF8kElceqARo02DOGH/SJW9lNV9qKZIqKiUhdWSo1UA9fJJMoW0dqoYAvkkl0ZanUQgFfJJMoW0dqoYAvkkmUrSO1iCrgm9lpZva+mZWbWdhE/9B2x5vZSjMrMbMJ0exTRGqhK0ulFtFm6awATgFm1LSBmTUA7gCGAeuBt8xsoXPugyj3LSLhaGEUqUFUAd859yGAmdW2WV+gxDm3OrTtA8AoQAFfRCSBEjGG3wFYV+nx+tBzezGzsWZWbGbFmzZtSkDTRESCo84evpk9D/w4zEsTnXMLYtkY59xMYCb44mmxfG8RkaCrM+A75/4ryn1sADpWepwTek5ERBIoEUM6bwH5ZtbZzBoDZwALE7BfERGpJNq0zF+Y2XqgP/CkmT0Tev4gM3sKwDm3CxgPPAN8CDzknHs/umaLiEh9RZul8xjwWJjnPwVOqPT4KeCpaPYlIiLR0ZW2IiIBoYAvItEpKoK8PL/MYl6e6u+nMAV8Edn3oF1U5Ovtr1nj19StqL+voJ+SFPBFgi6aoK36+2lFAV8k6KIJ2omov68ho5hRwBcJumiCdrzr72vIKKYU8EWCLpqgHe/6+xoyiikFfJGgiyZox7v+vpZsjCkFfJGgizZoFxZCaSmUl/ufsazFryUbY0oBX0TiG7SjoSUbY0oBX0RSl5ZsjKlolzgUEYkvLdkYM+rhi4gEhAK+iEhAKOCLiASEAr6ISEAo4IuIBIQ555LdhrDMbBOwJsLN2wBfxrE5qUrHHSw67mDZ1+POdc61DfdCygb8+jCzYudcQbLbkWg67mDRcQdLPI5bQzoiIgGhgC8iEhCZEvBnJrsBSaLjDhYdd7DE/LgzYgxfRETqlik9fBERqYMCvohIQKRVwDez481spZmVmNmEMK/vZ2YPhl5famZ5iW9l7EVw3L8xsw/MbLmZvWBmucloZ6zVddyVtvulmTkzS/vUvUiO2cxOD/17v29m9ye6jfEQwd94JzN7ycyWhf7OT0hGO2PNzO4xsy/MbEUNr5uZTQ99LsvN7OioduicS4sb0AD4N3Aw0Bh4Dzis2jaXAHeH7p8BPJjsdifouIcA2aH7FwfluEPbNQdeAZYABcludwL+rfOBZUDL0ON2yW53go57JnBx6P5hQGmy2x2jY/8pcDSwoobXTwD+CRjQD1gazf7SqYffFyhxzq12zn0PPACMqrbNKGBu6P4jwFAzswS2MR7qPG7n3EvOuYqVnpcAOQluYzxE8u8N8AfgRmBHIhsXJ5Ec84XAHc65LQDOuS8S3MZ4iOS4HfCj0P0DgE8T2L64cc69ApTVssko4O/OWwK0MLP2+7q/dAr4HYB1lR6vDz0Xdhvn3C5gK9A6Ia2Ln0iOu7Lz8T2CdFfncYdObzs6555MZMPiKJJ/6y5AFzN7zcyWmNnxCWtd/ERy3JOBM81sPfAU8OvENC3p6vv/v1Za8SqDmNmZQAFwbLLbEm9mlgX8CTgnyU1JtIb4YZ3B+DO5V8ysh3PuP0ltVfyNBuY45241s/7AfWZ2hHOuPNkNSyfp1MPfAHSs9Dgn9FzYbcysIf7Ub3NCWhc/kRw3ZvZfwERgpHPuuwS1LZ7qOu7mwBHAIjMrxY9vLkzzidtI/q3XAwudczudc58AH+O/ANJZJMd9PvAQgHPuDaAJvrhYpovo/3+k0ingvwXkm1lnM2uMn5RdWG2bhcCY0P1TgRddaOYjjdV53GbWC5iBD/aZMKYLdRy3c26rc66Ncy7POZeHn7sY6ZwrTk5zYyKSv/HH8b17zKwNfohndSIbGQeRHPdaYCiAmXXHB/xNCW1lciwEzg5l6/QDtjrnNu7rm6XNkI5zbpeZjQeewc/q3+Oce9/MpgDFzrmFwGz8qV4JfiLkjOS1ODYiPO6bgWbAw6E56rXOuZFJa3QMRHjcGSXCY34G+JmZfQD8AFzpnEvrs9gIj/t/gFlmdgV+AvecDOjMYWbz8V/gbULzE9cCjQCcc3fj5ytOAEqA7cC5Ue0vAz4zERGJQDoN6YiISBQU8EVEAkIBX0QkIBTwRUQCQgFfRCQgFPBFRAJCAV9EJCD+P806wnEUV8BRAAAAAElFTkSuQmCC\n" }, "metadata": { "needs_background": "light" } }, { "output_type": "display_data", "data": { "text/plain": [ "
" ], "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXwAAAD4CAYAAADvsV2wAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nO3deXyU1b3H8c8vAQQERZZWMUBwQVBQ0IAg0IK4oKB0Ua8aEbdyEbGVqvdWqZRLS+vSVusVRShqK0GttlaqtLSi1BU1XgS3QlEDxBWDpSKiYM7940zMwkwyyTzzzPJ836/XvJKZeTLPeQb9zZnfOed3zDmHiIjkv4JMN0BERMKhgC8iEhEK+CIiEaGALyISEQr4IiIR0SrTDUika9eurri4ONPNEBHJKS+++OKHzrlu8Z7L2oBfXFxMeXl5ppshIpJTzGxDoueU0hERiQgFfBGRiFDAFxGJiKzN4YtIMHbu3EllZSU7duzIdFMkQG3btqWoqIjWrVsn/TcK+CJ5rrKyko4dO1JcXIyZZbo5EgDnHFVVVVRWVtK7d++k/04pnXxQVgbFxVBQ4H+WlWW6RZJFduzYQZcuXRTs84iZ0aVLl2Z/a1MPP9eVlcHkybB9u7+/YYO/D1Bamrl2SVZRsM8/Lfk3VQ8/182YURvsa2zf7h8XEalDAT/XbdwY//ENCddeiGS9UaNGNWvh5f3330+/fv0YPXo05eXlfPe73wVgxYoVPPPMM80694oVKxg/fnzKx6Rq7NixdOrUKdDzKKWT63r2jB/czXy6R2kdiYCFCxeyYMECRowYAUBJSQngA3OHDh045phjMtm8FrnyyivZvn07t99+e2CvqR5+rpszxwf3hpxTWkeyQkVFBX379qW0tJR+/fpx2mmnsT2Whly+fDmDBg1iwIABXHDBBXz22Wf1/vaOO+7gsssu+/L+ggULmD59er1jZs+ezVNPPcWFF17IlVde+WXvu6Kignnz5nHjjTcycOBAnnzyyXp/9/zzzzNs2DAGDRrEMcccw9q1a3dr+6xZs5g4cSLDhg3j4IMPZsGCBV8+t23bNk477bQvr61m98DZs2czePBg+vfvz+TJk2nproJjxoyhY8eOLfrbRNTDz3WlpXDOOfGfS5Tukci67DJ46aVgX3PgQLjppsaPWbt2LQsXLmT48OFccMEF3HrrrUybNo3zzjuP5cuX06dPH84991xuu+22egH+jDPOYM6cOdxwww20bt2aO++8c7ce78yZM3nsscf4+c9/TklJCStWrAB8Pa4pU6bQoUMHrrjiit3a1LdvX5588klatWrFo48+ytVXX83vf//73Y5bs2YNK1eu5JNPPmHQoEGMGzcOgFWrVvHqq6/SvXt3hg8fztNPP82IESOYNm0aM2fOBGDixIk8/PDDnHLKKfVes6ysjBtuuGG3cx100EE88MADjb+ZKVDAzwe9esVP6/TsGX5bROLo0aMHw4cPB+Ccc87h5ptv5vjjj6d379706dMHgEmTJjF37tx6Ab9Dhw4ce+yxPPzww/Tr14+dO3cyYMCAQNq0detWJk2axD//+U/MjJ07d8Y9bsKECbRr14527doxevRonn/+eTp16sSQIUMoKioCYODAgVRUVDBixAgef/xxrr/+erZv386WLVs47LDDdgv4paWllGYg3aqAnw/mzKk/NROgfXv/uEgdTfXE06XhFMLmTCm86KKL+OlPf0rfvn05//zzA2vTNddcw+jRo3nwwQepqKhg1KhRcY9L1PY99tjjy8cKCwvZtWsXO3bsYOrUqZSXl9OjRw9mzZoVd658pnr4yuHng9JSmD/f9/TN/M/58zVgK1lj48aNPPvsswAsXryYESNGcMghh1BRUcH69esBuPvuu/n617++298effTRbNq0icWLF3PWWWc167wdO3bk448/jvvc1q1b2X///QG46667Er7GQw89xI4dO6iqqmLFihUMHjw44bE1wb1r165s27YtYfAuLS3lpZde2u2WzmAPCvj5o7QUKiqgutr/VLCXLHLIIYcwd+5c+vXrx0cffcTFF19M27ZtufPOOzn99NMZMGAABQUFTJkyJe7fn3HGGQwfPpx99tmnWec95ZRTePDBB+MO2v7Xf/0XV111FYMGDWLXrl0JX+Pwww9n9OjRDB06lGuuuYbu3bsnPLZTp0585zvfoX///px44omNfjg0ZeTIkZx++uksX76coqIili1b1uLXqmEtHUFOt5KSEqcNUERS9/rrr9OvX7+Mnb+iooLx48fzyiuvtPg1xo8fz/Tp0xkzZkyALWvarFmzEg76ZoN4/7Zm9qJzriTe8erhi0jW+te//kWfPn1o165d6ME+H2nQVkTSqri4uMW9+06dOrFu3bqAW5S8WbNmZezc6aAevkgEZGvqVlquJf+mCvgiea5t27ZUVVUp6OeRmnr4bdu2bdbfKaUjkueKioqorKxk8+bNmW6KBKhmx6vmUMAXyXOtW7du1q5Ikr+U0hERiQgFfBGRiFDAFxGJiEACvpndYWYfmFncybbm3Wxm681sjZkdGcR5RUQkeUH18O8Cxjby/EnAwbHbZOC2gM4rySgrg+JiKCjwP8vKMt0iEcmAQAK+c+4JYEsjh0wAfuu8lUAnM9sviHNLE8rKfOnkDRv8LlgbNvj7CvoikRNWDn9/YFOd+5Wxx+oxs8lmVm5m5ZozHJAZM+rXyQd/X9sfikROVg3aOufmO+dKnHMl3bp1y3Rz8kOibQ61/aFI5IQV8N8GetS5XxR7TNIt0TaH2v5QJHLCCvhLgHNjs3WGAludc++GdO5omzPHb3dYl7Y/FImkoKZl3gM8CxxiZpVmdqGZTTGzmu1rlgJvAuuBBcDUIM6b08KaOZPN2x9q9pBIqLTjVSbUzJxpuOl4tgTiMOg9EEmLxna8UsDPhOJiPz2yoV69/H60UaD3QCQttMVhtknHzJlcS49o9pBI6BTwMyHomTO5uLhKs4dEQqeAnwlBz5wJa3FVkN8iNHtIJHQK+JkQ9MyZMNIjQX+LyObZQyJ5SoO2+SCMAdBE56g5z5w5wQfrsjL/LWXjRp/qScc5RPKMBm3zXRjpkca+LaRjzCAXxyVEspwCfj4IIz3S1GBq0GMGKvomEjildCQ58RZKNWQG1dXBnK+gwPfs03kOkTyklI6kru63iESCnFKpaZsigVPAl+SVlvpB4EWL0j9moGmbIoFTwM9m2bp6NowxA03bFAmccvjZSsXFRKQFlMPPRZqlIiIBU8DPViouJiIBU8DPVumcpZKtYwMiklYK+NkqXbNUtIJVJLIU8LNVumapaGxAJLI0SydqtIJVJK9plo7U0gpWkchSwI8arWAViSwF/KjRClaRyGqV6QZIBpSWKsCLRJB6+CIiEaGALyISEQr4khu0OlgkZcrhS/ZrWDm0ZnUwaCxCpBnyrodfXQ2//CX85S9QWRl/jZHkGK0OFglE3vXwKyvh8str7++9Nxx2GPTvX//WrVvm2ijNpMqhIoHIu4Dfsyd8+CG8+iq88krtzwce8NPNa+y7LwwYAIcfXvuzXz9o2zZzbZcEevb0aZx4j4tI0vIu4AN06QJf+5q/1XAO3nvPB/+XX/a3NWtg7lzYscMfU1gIffv64D9wIAwa5G9du2bmOiRmzpz4u39pdbBIs+RlwI/HDPbbz9+OP7728V27YP16H/zXrIHVq+Gpp+Cee2qPKSqqDf41t549/WtKCGoGZmfM8Gmcnj19sG84YFtW1vQxIhEWSLVMMxsL/AooBH7tnLu2wfPnATcAb8ceusU59+vGXjPT1TKrquCll2DVqtrb2rW1BSU7d67/LWDQIDjkEP8tQTJAewCLAI1Xy0w54JtZIbAOOB6oBF4AznLOvVbnmPOAEufctGRfN9MBP57t2/23gFWraj8M1qyBzz7zz7drB0ccASUlMHiw/6kPgZAUF8fP8/fqBRUVYbdGJGMaC/hBpHSGAOudc2/GTnYvMAF4rdG/ykHt28PQof5WY9cu+Mc/ar8F/N//wV13wS23+Of33BOOPLL2A6CkBA46SOmgwGkmj0iTggj4+wOb6tyvBI6Oc9y3zexr+G8D051zmxoeYGaTgckAPXNkBkarVrVTPSdO9I998QWsWwcvvADl5f526621g8OdOsFRR9V+Ezj6aD9OICnQTB6RJoU1aPsn4B7n3Gdm9p/Ab4BjGx7knJsPzAef0gmpbYErLPRTPPv1g3PP9Y/t3AmvvVb/Q+CXv/SPA3Tv7gP/0KH+51FHQYcOmbuGnKOZPCJNCmKl7dtAjzr3i6gdnAXAOVflnItluvk1cFQA5w1GUDVamnid1q19fv+ii2DePB/w//1vWLkSfvUrGDXKjwf893/73/fe2x8/eTLceacfMM67VcNB1sdRnX+RpjnnUrrhvyW8CfQG2gCrgcMaHLNfnd+/Caxs6nWPOuool3aLFjnXvr1zPpb6W/v2/vFMvI5zbvNm5x55xLmZM5074QTnOnWqfcnOnZ0bN865n/zEuccec+7jj5v98tkjwPdMRGoB5S5BXA1qWubJwE34aZl3OOfmmNns2ImXmNnPgFOBXcAW4GLn3D8ae81QZukENbMjjTNEqqt97/6ZZ+DZZ/3P11/3zxUU+G8Bw4bBMcfA8OG1Hdysp1k1ImmR1mmZ6RJKwC8oiJ8nMaudcB/m6yRpyxZ47rnaD4HnnoNt2/xz3bvDiBH+NnKkLxuRldNCQ37PRKIi3dMyc1dQMztCniHSuTOcdJK/gZ8V9PLL8PTT/vbkk/C73/nn9trL9/5rPgSOPjpL6gVpVo1I6PKuPHKzzJnjZ3LU1bq17y43ZyAx3uuEOEOksNCv+r3kEli8GDZt8lmRRYvgrLP8/R/+sHYweORIX4Fg2bLabwahy/B7JhJJiZL7mb6FMmjrnB8k7NXLOTPnunRxrk2blg0k1n2dXr2ybvCxqsq5hx5y7oornBsyxLnCQn95rVo5N2yYc1dd5dyyZc5t2xZiozL1nmX5v5VIKkj3oG06ZKS0QoQGErdt82MAf/87PP64Xx+wa5dfSDZkiP82MGqUTwftuWemWxsg1dyRPKdB22RFeCCx5gPg8cdhxQr/AfDFFz7D1fADoGEmJqdE6ENdokkBP1kKBl/6+OP6HwDl5bUfAEOHwrHHwpgxfhC4TZtMt7YZIvyhLtHQWMCP9qBtQxpI/FLHjnDiiXDttX418EcfwZ//DNOnw6efwuzZfoOZffaBcePgppt86Ygs7T/USjQLSLODJAIU8OvS8vyEOnaEsWPhuut8uqeqCh58EM4/328gM3263zu4Rw//2OLF8MEHITWuOSUa9KEuEaaUjgRiwwb429/gr3+F5cv94jDw00VPOMHvMjZiRBrWALRkEFY7Y0keUw5fQvXFF35vgL/+1X8IPP20rwratq1PA9V8AAwYEEAZCI27iNSjHL6EqrDQ1/q/+mo/6LtlCzzyCPznf0JlJVxxha8BtN9+fg+BRYtg8+YWnqyxjU+CrMYpkgfUw5fQVVbCo4/WfgP48EPf0y8pqS0ZMXhwkjWAEvXwu3Txo8uaby8Ro5SOZK3qanjxRT8D6C9/8YXgqqt9vD7hBB/8x46Fbt0SvECiHH67dn5kuSGleiTPKeBLzqiq8r3+mg+ADz7wvf8hQ/z0z3HjYNCgBrn/eIOwEydqvr1EkgK+5KTqaj/4+8gjsHQpPP+8j+Hdu/vAP348HHdcgpW/QQ3makaP5BgN2kpOKijwe/vOnOkXf733Htx1l9/o5d57YcIEn/oZNw5uu63B+G0Q8+1r0kUbNvhPmg0b/H0N/kqOUsCXnPGVr8CkSb7W/4cf+vn+F18M69bB1Km+8z5woP+AKD+klOrbm1hE19Qsnhkz6o8NgL8/Y0a6L1UkLZTSkbywdi386U+wZImf919d7VM/p5zivwkceyzssUedP0hmwZbq7kgOUg5fIqWqyuf8lyzxA7/btkGHDn7Gzze+ASefDJ0GFjed49eiLslByuFLpHTp4ifp3H+/X9C1dCmcfTY88YTvvHfrBtUbGlmwVUN1dyTPRCvga+Vl5LRt63v2t98O77zjSz5ffjm82yqJqpkqpid5JjopHe10JHWVlVF93vkU7Nr55UM7aM2Piu6k7QWlnHYa9O8fQK0fkZAppQOacSG7KSioH81bFRp77gk/+QkcfjgccghcdZXf/CVL+0UizRKdgN9YkS2Jnhkz4PPP6z3U6ovPmbljBu+8A/Pm+azfDTf4uj69e/uib889p+AvaZTmtHN0UjqacSF1JTnlcssWeOgheOABX/Jh506/ycsZZ/iB4N3KPIi0VEBpZ6V0QDMupL4ktzrs3Nnv4PXII76uz29/6xd33XyzXwV86KHw4x/DG2+E0GbJbyGknaMT8DXjQupqQQegUyc/3XPJEl/mYf58+OpX/cregw7ym7vffDO8/36a2y75KYS0c3RSOiINBVQYbdMmX9tn8WJ46SWfLTruOP9S3/gG7LVXGtou+SegtLNW2oqE5LXXfOBfvBjeesuvAzj1VJ/vP+kkaNMm0y2UrKUcvkhuOfRQP63zjTf8Iq8LL4THHvM9/X339f8///3vKsUjcYSQdlYPXyTNdu70WzouXgwPPgiffAJFRXDWWb7nf8QRmukjwVEPXySDWrf26Zy77/YDuosX+5k+N97op3X27++HD956K9MtlXwXSMA3s7FmttbM1pvZD+I8v4eZ3Rd7/jkzKw7ivCKhS3FhzJ57+p79n/7kZ/rcdpsv9vbDH8IBB8Axx8Att/gpoCJBSzngm1khMBc4CTgUOMvMDm1w2IXAR865g4AbgetSPa9I6ALeAatLF5gyxVfxrKiAa6/1pZwvvdTX8j/pJFi0yD8mEoSUc/hmNgyY5Zw7MXb/KgDn3M/qHLMsdsyzZtYKeA/o5ho5uXL4knVCWq398su1M302boR27fwmLmefDSeeqJk+0rh05/D3BzbVuV8ZeyzuMc65XcBWoEuchk42s3IzK9+8eXMATRMJUEj1mAYMgJ/9zOf0n3wSzjvPl3U49VTYbz//reDJJzXTR5ovqwZtnXPznXMlzrmSbt26Zbo5IvUlWY4hKAUFMGIE3Hqrr+X/8MO+h3/33fC1r/mCbj/4gf9GIBmSY3tsBBHw3wZ61LlfFHss7jGxlM7eQFUA5xYJTwbrMbVpA+PG+TTP++/73P5hh8HPf+5LOR95JPzv//rN3SUkAY/phCGIgP8CcLCZ9TazNsCZwJIGxywBJsV+Pw14rLH8vUhWypJ6TB06+FMuXQrvvuvr95jBd7/rB3u//W1f7G3XrlCbFT05uMdGygE/lpOfBiwDXgd+55x71cxmm9mpscMWAl3MbD3wfWC3qZsiGZfM1/PSUj9AW13tf2a4+F63bn5Wz4svwurVMG2az++PH+8zTTNmqJJn2uTgHhtaaSsCebUF5uef+x7+woXw5z/7z6bRo+Gii+Bb3/L1fSQAWbrHhlbaijQlB7+eJ9KmDXzzm36Qd8MGX6+/5stI9+4+9bN6daZbmQdycI8NBXwRyMmv58koKvKreNev9/V8TjwRbr/dl3YYPNj//u9/Z7qVOSpLxnSaQykdEcjar+fpUFXlM1gLFsArr/hO6emn+5TP8OEq5JbrlNIRaUoOfj1vqS5dfFpnzRq/Kfs558Dvfw8jR/ryzr/6FXz0UaZbKemggC8COfn1PFVmMGSIT+u8+64f5N17b7jsMth/f7jgAnj++fh7vUtuUkpHROpZtQrmzfNpn08+8Yu6pkzxVT47dMh066QpSumISNIGDfK9/nfe8WUddu70M1a7d4dLLlEph1ymgC8ice21F1x8sZ/C+fTTfpvGhQt9KYcRI3x5hx07Mt1KaQ4FfBFplJnfmOW3v4XKSl+/5/33YeJEP+3zyivhn//MdCslGQr4IpK0rl3h8sth7Vo/r3/0aLjpJujTB44/3u/kpbLN2UsBX0SaraAAxoyB++/3a9Nmz4bXX/c1+/v08R8CW7dmupVpkmMlketSwBeRlOy3H1xzjd+w5b774KtfhenTfbrn0kth3bpMtzBA6S6JnOYPE03LFJHAlZf7ss333utn+Zx0Enzve3DCCTm+kjedK7IDKuDX2LRMBXwRSZv33vNTPG+7zQ/09u3re/3nnpujc/oLCuKvRDNLffAioA8TzcMXkYzYd1/40Y98nv/uu32Qv+QSn+65/HKfBsop6dzmMoQCfgr4IpJ2bdr4mj3PPw/PPANjx/qaPQce6Of3P/54jpRwSGfNpRD2TFbAF5HQmMGwYT63X1EBV10FTz0Fxx4LRxwBv/41fPppplvZiHTWXAqhgJ9y+CKSUZ9+Cvfc43v8a9ZA585+7HLqVOjRI9OtC1lZmd90Z+NG37OfM6fZHyYatBWRrOccPPGEn93zxz/6DvS3vuVn9xxzTI7P7gmRBm1FJOuZwde/7mvzv/EGfP/78Le/+bo9gwf72j07d4bcqBxeZBWPAr6IZJ3iYrj+el+7Z948PzV94kQ44ABfyyctq3gbBvepU9O7yCoDlNIRkaxXXQ3Llvlg/9hj0LGjj73f+15Aef54i57M4k8dyvJtL5XSEZGcVlDgV+suXw4vvginnOLr9RxwgJ/uuWpViieYMaN+sIfE80RzeGN7BXwRySlHHuk75G++6ffmfegh/9hxx8Ff/tLC+fzNCeLNmRefZWMACvgikpN69oRf/AI2bfL5/n/8w38LOOKIFgzwJgriDacGNWdefLoLrbVA/gX8LPtEFZH06tTJb8Ly5pvwm9/4fP/EiXDQQX5u/7ZtSbxIokVPU6a0fJFVvDTR9u3+8UxxzmXl7aijjnLNtmiRc+3bO+c/T/2tfXv/uIhEwhdfOPfww86NHOlDwD77OHfNNc598EETf7hokXO9ejln5n8mGzcS/Z1Z/VhUczNr8bUlAyh3CeJqfs3SSWfpUhHJOc8+C9dd5/P87drBd77ji7YFVp6msZLG3/seVFXt/jddusCHHwbUgN1FZ5ZOCNXmRCR3DBvmV+2+9hr8x3/Arbf6gm2TJvnHUpaNaZtG5FfAD6HanIjknn794M47fZ5/2jR44AE47DD45jfhhRdSeOHGOplbtsR/LtHjIcivgB9CtTkRyV09esCNN/rM78yZsGIFDBnid+JasaIFUzob62RmYQc0vwJ+OkuXikje6NoV/ud/fEf8+ut9lc7Ro2H4cHjkkWYE/sY6mVnYAU1p0NbMOgP3AcVABXCGc+6jOMd9Abwcu7vROXdqU6+t0goiEpZPP/Upn+uu8x8CgwbBD3/oN2cpaKpbXLekcefO/rEtW3xP/uSTYenSlModN1fayiOb2fXAFufctWb2A2Af59x/xzlum3OuWTtYKuCLSNh27vSLtn76U1i/Hg49FK6+2g/4tmrVxB8HtAl5qtI5S2cC8JvY778BvpHi64mIZEzr1nD++X7V7uLFvnd/zjl+0Peuu5pYvZsDM3ZSDfhfdc69G/v9PeCrCY5ra2blZrbSzBJ+KJjZ5Nhx5Zs3b06xaSIiLVNYCGedBatXwx/+AHvt5T8I+vSBBQvg88/j/FGiGTvx1gZlSJMB38weNbNX4twm1D0utsIrUX6oV+wrxtnATWZ2YLyDnHPznXMlzrmSbt26NfdaREQCVVDgp26Wl8PDD8NXvuKzNgcdBLfdBp99VufgxurxZEmJlyYDvnPuOOdc/zi3h4D3zWw/gNjPDxK8xtuxn28CK4BBgV2BiEiamcG4cbBypa/IWVTk90c58ECYOxd27MAPyMbbh9G5rEnrpJrSWQJMiv0+CXio4QFmto+Z7RH7vSswHAhijZuISKjM4MQT4emn4dFHfT3+adN84L/lo1ISToLJktX+qQb8a4HjzeyfwHGx+5hZiZn9OnZMP6DczFYDjwPXOucU8EUkZ5nBmDHw97/7TVkOPBAuvRTeLuwV/w+yZLV/SgHfOVflnBvjnDs4lvrZEnu83Dl3Uez3Z5xzA5xzR8R+Lgyi4SIi9WSgNLoZHHtsbeBfeOAcPiG7FlvVlV8rbUUkmjK82UhN4J/5j1Leumo+7+3Ri2qMysJeLD9zPp+fnh2r/fOrPLKIRFOWlUZ3Dv72N/jRj/xAb+/e/vfS0iQWcKUoOuWRRaIsyru9ZVlpdDNfkO2ZZ3xlhc6d4bzzoH9/uO8+vytXJijgi+SDLNw/NVRZWJkSfOA/6SRfgvkPf/C9+zPP9LV6lixp4YbrKVDAF8kHObCsP62ysDJlXWZ+Adfq1f4zePt2mDABhg71qZ+wAr8Cvkg+yLKURuiaWxo9Q+mvwkI4+2x4/XVYuBDee8+nfkaNgqeeSv/5FfBF8kGWpjRCVVrqB2irq/3PRMF+6lSYODGj6a9WreCCC2DdOrjlFv9z5EgYOzbFHbiaoIAvkg+yPKWRNcrKYN683XMoGUp/7bEHXHIJvPEG3HCDr9kzZIjP86cjzaOAL5IPtNtbcmbMSBxJM5j+at8errgC3noLZs/2JRvileVJlQK+SK5qmIeG5FIaUdLwPWqsVHEWpL86doRrrvEbsKRDmpcAiEhaNNxdqSYPDQr0NeK9R2bxe/hmkUh/qYcvkouiPg0zGfHeI+d2z5WYwZQpkfigVMAXyUVRn4aZjETvhXP1xzruvhtuvTXctmWIUjoiuahnz/j56CzIQ2eNRO9RhurrZAP18EVykaZhNk3v0W4U8EVykaZhNk3v0W6U0hHJVaWlkQ5eSal5f2bM8Dn9mkHtiL5vCvgikr80fbUepXREJH9p+mo9Cvgikr80fbUeBXyRfBHlHa8SURXRehTwRfJB1He8SkRTM+tRwBfJB8pVx6epmfWYC3tTxSSVlJS48vLyTDdDJDcUFCQuCpapHbMlI8zsRedcSbzn1MMXyQfKVUsSFPBF8oFy1ZIEBXyRfFA3Vw1+t+yaHH7UB27lS1ppK5IvagYitbJUElAPXySfaLaONEIBXySfaGWpNEIBXySfaLaONEIBXySfaLaONCKlgG9mp5vZq2ZWbWZxJ/rHjhtrZmvNbL2Z/SCVc4pII7SyVBqR6iydV4BvAbcnOsDMCoG5wPFAJfCCmS1xzr2W4rlFJB5tjCIJpESWABUAAAV+SURBVBTwnXOvA5hZY4cNAdY7596MHXsvMAFQwBcRCVEYOfz9gU117lfGHtuNmU02s3IzK9+8eXMITRMRiY4me/hm9iiwb5ynZjjnHgqyMc65+cB88MXTgnxtEZGoazLgO+eOS/EcbwM96twvij0mIiIhCiOl8wJwsJn1NrM2wJnAkhDOKyIidaQ6LfObZlYJDAMeMbNlsce7m9lSAOfcLmAasAx4Hfidc+7V1JotIiLNleosnQeBB+M8/g5wcp37S4GlqZxLRERSo5W2IiIRoYAvIqkpK4PiYr/NYnGx6u9nMQV8EWl50C4r8/X2N2zwe+rW1N9X0M9KCvgiUZdK0Fb9/ZyigC8SdakE7TDq7ytlFBgFfJGoSyVop7v+vlJGgVLAF4m6VIJ2uuvvK2UUKAV8kahLJWinu/6+tmwMlAK+SNSlGrRLS6GiAqqr/c8ga/Fry8ZAKeCLSHqDdiq0ZWOgFPBFJHtpy8ZApbrFoYhIemnLxsCohy8iEhEK+CIiEaGALyISEQr4IiIRoYAvIhIR5pzLdBviMrPNwIYkD+8KfJjG5mQrXXe06LqjpaXX3cs51y3eE1kb8JvDzMqdcyWZbkfYdN3RouuOlnRct1I6IiIRoYAvIhIR+RLw52e6ARmi644WXXe0BH7deZHDFxGRpuVLD19ERJqggC8iEhE5FfDNbKyZrTWz9Wb2gzjP72Fm98Wef87MisNvZfCSuO7vm9lrZrbGzJabWa9MtDNoTV13neO+bWbOzHJ+6l4y12xmZ8T+vV81s8VhtzEdkvhvvKeZPW5mq2L/nZ+ciXYGzczuMLMPzOyVBM+bmd0ce1/WmNmRKZ3QOZcTN6AQeAM4AGgDrAYObXDMVGBe7Pczgfsy3e6Qrns00D72+8VRue7YcR2BJ4CVQEmm2x3Cv/XBwCpgn9j9r2S63SFd93zg4tjvhwIVmW53QNf+NeBI4JUEz58M/BkwYCjwXCrny6Ue/hBgvXPuTefc58C9wIQGx0wAfhP7/QFgjJlZiG1Mhyav2zn3uHOuZqfnlUBRyG1Mh2T+vQF+DFwH7AizcWmSzDV/B5jrnPsIwDn3QchtTIdkrtsBe8V+3xt4J8T2pY1z7glgSyOHTAB+67yVQCcz26+l58ulgL8/sKnO/crYY3GPcc7tArYCXUJpXfokc911XYjvEeS6Jq879vW2h3PukTAblkbJ/Fv3AfqY2dNmttLMxobWuvRJ5rpnAeeYWSWwFLg0nKZlXHP//2+UdrzKI2Z2DlACfD3TbUk3MysAfgmcl+GmhK0VPq0zCv9N7gkzG+Cc+1dGW5V+ZwF3Oed+YWbDgLvNrL9zrjrTDcsludTDfxvoUed+UeyxuMeYWSv8V7+qUFqXPslcN2Z2HDADONU591lIbUunpq67I9AfWGFmFfj85pIcH7hN5t+6EljinNvpnHsLWIf/AMhlyVz3hcDvAJxzzwJt8cXF8l1S//8nK5cC/gvAwWbW28za4AdllzQ4ZgkwKfb7acBjLjbykcOavG4zGwTcjg/2+ZDThSau2zm31TnX1TlX7Jwrxo9dnOqcK89McwORzH/jf8T37jGzrvgUz5thNjINkrnujcAYADPrhw/4m0NtZWYsAc6NzdYZCmx1zr3b0hfLmZSOc26XmU0DluFH9e9wzr1qZrOBcufcEmAh/qveevxAyJmZa3EwkrzuG4AOwP2xMeqNzrlTM9boACR53XklyWteBpxgZq8BXwBXOudy+ltsktd9ObDAzKbjB3DPy4POHGZ2D/4DvGtsfOJHQGsA59w8/HjFycB6YDtwfkrny4P3TEREkpBLKR0REUmBAr6ISEQo4IuIRIQCvohIRCjgi4hEhAK+iEhEKOCLiETE/wM/pzgNfIYDLAAAAABJRU5ErkJggg==\n" }, "metadata": { "needs_background": "light" } }, { "output_type": "display_data", "data": { "text/plain": [ "
" ], "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXwAAAD4CAYAAADvsV2wAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nO3de3xU1b338c8PiFwqEkiiAgGC5WJELJd4odBHKVotRekFbW28YzlobY+0tTceLS+fQ1trH9vHc6wKorUSq6d6VEptaYuiVkGNl1q8gFQDBLWEoFSkKJH1/LFmmFxmJpPMnuv+vl+v/ZrMzM7stSfw22uv9VtrmXMOEREpfj1yXQAREckOBXwRkZBQwBcRCQkFfBGRkFDAFxEJiV65LkAi5eXlrqqqKtfFEBEpKM8888wO51xFvPfyNuBXVVVRX1+f62KIiBQUM9uc6D016YiIhIQCvohISCjgi4iERN624YtI1+zbt4/Gxkb27t2b66JIFvTp04fKykpKSkpS/h0FfJEi0djYSP/+/amqqsLMcl0cySDnHM3NzTQ2NjJy5MiUf09NOsWgrg6qqqBHD/9YV5frEkkO7N27l7KyMgX7EDAzysrKunw3pxp+oaurg3nzYM8e/3zzZv8coLY2d+WSnFCwD4/u/K1Vwy90CxfGgn3Unj3+dRGRVhTwC92WLfFf35xw7IVI3jjppJO6NMDyN7/5DdXV1UyfPp36+nq+/vWvA7BmzRqeeOKJLh17zZo1zJo1K+190nXaaadRWlra4Tivv/46xx9/PKNGjeKLX/wiH3zwQdrHUsAvdMOHx3/dTG35UnSWLVvG0qVLefjhh6mpqeH6668Huhfw88UVV1zBHXfc0eH173znOyxYsIBNmzYxcOBAli1blvaxFPAL3eLFPri355yadSSrGhoaOPLII6mtraW6upo5c+awJ9LcuHr1aiZOnMj48eO56KKLeP/999v87q233srll19+4PnSpUtZsGBBm32uvvpq/vKXvzB37lyuuOKKA7XvhoYGbrrpJn72s58xYcIEHnvssTa/99RTTzFlyhQmTpzIxz/+cTZs2NCh7IsWLeLcc89lypQpjB49mqVLlx54b/fu3cyZM+fAuUVXCbz66qs59thjOfroo5k3bx7dXT1wxowZ9O/fv81rzjkeeugh5syZA8D555/P/fff363Pb02dtoWuthbOOSf+e4mae6ToXX45PP98sJ85YQL8/OfJ99mwYQPLli1j6tSpXHTRRfziF7/gsssu44ILLmD16tWMGTOG8847jxtvvLFNgD/rrLNYvHgx1157LSUlJdx2223cfPPNbT77qquu4qGHHuKnP/0pNTU1rFmzBvDzbs2fP5+DDz6Yb33rWx3KdOSRR/LYY4/Rq1cv/vznP/P973+fe++9t8N+L7zwAuvWreO9995j4sSJfOYznwHgueee48UXX2TIkCFMnTqVxx9/nGnTpnHZZZdx1VVXAXDuueeycuVKTj/99DafWVdXx7XXXtvhWKNGjeKee+5J+D02NzdTWlpKr14+RFdWVrJt27aE+6dKAb8YjBgRv80+UXOPSIYMGzaMqVOnAnDOOedw/fXXc8oppzBy5EjGjBkD+NrqDTfc0CbgH3zwwXzyk59k5cqVVFdXs2/fPsaPHx9ImXbt2sX555/Pq6++ipmxb9++uPvNnj2bvn370rdvX6ZPn85TTz1FaWkpxx13HJWVlQBMmDCBhoYGpk2bxsMPP8xPfvIT9uzZw86dOxk3blyHgF9bW0ttHmXLKeAXg8WL26ZmAvTr51+XUOqsJp4p7VMFu5I6ePHFF/PDH/6QI488kgsvvDCwMl155ZVMnz6d++67j4aGBk466aS4+yUqe+/evQ+81rNnT1paWti7dy+XXnop9fX1DBs2jEWLFsXNie9uDb+srIx33nmHlpYWevXqRWNjI0OHDk3ldJNSG34xqK2FJUt8Td/MPy5Zojx8ybotW7awdu1aAO68806mTZvG2LFjaWhoYNOmTQDccccdnHjiiR1+9/jjj2fr1q3ceeednH322V06bv/+/Xn33Xfjvrdr164DwfKXv/xlws944IEH2Lt3L83NzaxZs4Zjjz024b7R4F5eXs7u3bsTBu/a2lqef/75DluyYA/+YjN9+vQD+91+++3Mnj076e+kQgG/WNTWQkMD7N/vHxXsJQfGjh3LDTfcQHV1NW+//TaXXHIJffr04bbbbuPMM89k/Pjx9OjRg/nz58f9/bPOOoupU6cycODALh339NNP57777ovbafvtb3+b733ve0ycOJGWlpaEn3HMMccwffp0TjjhBK688kqGDBmScN/S0lK+8pWvcPTRR3PqqacmvTh05hOf+ARnnnkmq1evprKyklWrVgFwzTXXcN111zFq1Ciam5uZO3dut48RZd3tWc60mpoapwVQRFL38ssvU11dnbPjNzQ0MGvWLNavX9/tz5g1axYLFixgxowZAZasc4sWLUrY6ZvP4v3NzewZ51xNvP1VwxeRnHvnnXcYM2YMffv2zXqwDxN12opIIKqqqrpduy8tLWXjxo0Blyh1ixYtytmxs0k1fJEikq9NtBK87vytFfBFikSfPn1obm5W0A+B6Hz4ffr06dLvqUlHpEhUVlbS2NhIU1NTrosiWRBd8aorFPBFikRJSUmXVj+S8FGTjohISCjgi4iEhAK+iEhIBBLwzexWM9tuZnGTcM273sw2mdkLZjYpiOOKiEjqgqrh/xI4Lcn7nwZGR7Z5wI0BHVdSUVcHVVXQo4d/1EpYIqEUSMB3zj0K7Eyyy2zgV85bB5Sa2eAgji2dqKvzUydv3uxXwdq82T9X0BcJnWy14Q8FtrZ63hh5rQ0zm2dm9WZWr1zigCxc2HaefPDPtfyhSOjkVaetc26Jc67GOVdTUVGR6+IUh0TLHGr5Q5HQyVbA3wYMa/W8MvKaZFqiZQ61/KFI6GQr4K8Azotk65wA7HLOvZmlY4fb4sV+ucPWtPyhSCgFlZb5a2AtMNbMGs1srpnNN7PosjYPAq8Bm4ClwKVBHLegZStzJp+XP1T2kEhWacWrXIhmzrRfdDxfAnE26DsQyYhkK14p4OdCVZVPj2xvxAi/Hm0Y6DsQyQgtcZhvMpE5U2jNI8oeEsk6BfxcCDpzphAHVyl7SCTrFPBzIejMmWwNrgryLkLZQyJZp4CfC0FnzmSjeSTou4h8zh4SKVLqtC0G2egATXSM6HEWLw4+WNfV+buULVt8U08mjiFSZNRpW+yy0TyS7G4hE30GhdgvIZLnFPCLQTaaRzrrTA26z0CTvokETk06kpp4A6XaM4P9+4M5Xo8evmafyWOIFCE16Uj6Wt9FJBJkSqXSNkUCp4Avqaut9Z3Ay5dnvs9AaZsigVPAz2f5Ono2G30GStsUCZza8POVJhcTkW5QG34hUpaKiARMAT9faXIxEQmYAn6+ymSWSr72DYhIRing56tMZaloBKtIaCng56tMZamob0AktJSlEzYawSpS1JSlIzEawSoSWgr4YaMRrCKhpYAfNhrBKhJavXJdAMmB2loFeJEQUg1fRCQkFPBFREJCAV8Kg0YHi6RNbfiS/9rPHBodHQzqixDpAtXwJf9pdLBIIBTwJf9p5lCRQCjgS/7T6GCRQCjgS/7T6GCRQCjgS/5LdXSwMnlEkgok4JvZaWa2wcw2mdl347x/gZk1mdnzke3iII4rIVJbCw0NfkbPhob4wV7z/IsklXbAN7OewA3Ap4GjgLPN7Kg4u97tnJsQ2W5J97gibSiTR6RTQdTwjwM2Oedec859ANwFzA7gc0VSp0wekU4FEfCHAltbPW+MvNbeF8zsBTO7x8yGxfsgM5tnZvVmVt/U1BRA0SQ0lMkj0qlsddr+Fqhyzh0D/Am4Pd5Ozrklzrka51xNRUVFloomRUGZPCKdCiLgbwNa19grI68d4Jxrds69H3l6CzA5gOMGI6jMDmWIdF2Q35nm+RfpnHMurQ0/H89rwEjgIOCvwLh2+wxu9fPngHWdfe7kyZNdxi1f7ly/fs75vA6/9evnX8/F54SJvjORjADqXYK4Gsgi5mY2E/g50BO41Tm32Myujhx4hZn9CDgDaAF2Apc4515J9plZWcS8qsqn77U3YoRP/cv254SJvjORjEi2iHkgAT8TshLwe/Twdcv2zHy+d7Y/J0z0nYlkRLKAH+6RtkFldihDpOv0nYlkXbgDfrzMjpIS2L27ax2JyhDpOn1nIlkX7oDfPrOjrMw/Njd3bXi+MkS6LpffmTKqJKTC3YbfnjoSi1/71bPA31noAi1FQm34qdLw/OKnOXckxBTwW1NHYvHTRV1CTAG/NXUkFj9d1CXEFPBbU+drYepKJ6wu6hJivXJdgLxTW6sAX0jad8JGM6sg/t8x+trChb4ZZ/hwH+z1N5cQUJaOFDZlVom0oSwdKV7JOmGVby/ShgK+FLZEna2DBmmNW5F2FPClsCXqhAXl24u0U3Sdtu++C6edBqWlMGCAf0y0DRwY+7mkJNcll25J1Al77rnx91e+vYRY0QX8Dz7wFbzt22HjRnjnHb+1tCT/vYMP9heAQYNij2VlyR8HDYLevbNzXpJEvMyqhQvjd+Z2Nd++rk4ZPVI0ii7gl5XBn/7U9jXn/N38rl3w9tv+AvD22x1/jm47d8Irr/jH5mbYty/x8T7ykbYXgrIyKC+Hior4j+XlcNBBmf0OBB+Y482Z05V8+66mfIrkOaVldsI5eO+9WPBv/dj+tei2Y4d/LdFXe8ghsYtARQUcemjH7bDD/GN5OfQqustylnRWO+/sfaV8SgHSilc50NLi7xaamvy2Y0f8x+3bY4+Jmp3Ky+Hww2PbYYfFf15W5jMQJQWpzJqpVbmkACngFwDnfPPS9u1++8c/Yo/R7a23YtvevR0/o1cvH/yHDPHb4MFtH6NbebkuDCnV3lXDlwKULOCrsSBPmPnO4oEDYezY5Ps657OR3nordiF4802/RX9+7TV4/HF/J9FeSYm/CFRWtt2GDYv9PHgw9OyZmXPNC6nMmhlEP4BIHglXwC+SjAsz3w9wyCEwZkzyfT/4wF8E3njDb9u2+e2NN6CxEZ57Dn77W/jXv9r+Xs+esYtC6wvBsGF+Gz7cNyMV7J3C8OGdZ/Fo3h0pMuFp0tFKRwk55/sbtm71F4HGxo4/b93a8aJQUuKDf1VVbBs5MvZzXt8l1NXBhRe2TcEqKYHbbgv9vwcpbGrSgeQrHYX8P7hZbFzBxz4Wf59oH0M0+G/Z4reGBl9R/v3vfVNSayUlvlLc+iLQehs8OMd3CGbJn4sUmfDU8JVxkXH/+lfsIhDdXn899vM//tF2/4MO8v2fRxwBH/2o36I/H3GEH+OQMeqQlXwUQLOzaviQWputpKVvX9/hnKjTec8e/++49UXg9dd9B/OTT/o7iNYOP9wH/1Gj/DZ6dOzxkEPSLKyWOpR8k4WBfuGp4asNP+/t3Al//3v8bdu2tvsedpgP/KNH+47r6OOoUf7C0ynV8CXfBPRvUnn4UUWSpRNG773n7wRefTW2bdzoH996K7ZfdGXK6J3G2LFw5JFQXe3vGA4006sCIPkmoGZnBXwpav/8J2za5C8AGza03d57L7bfgAE++Ee3U5rqOOauhfR6cwumCoDkmmr4CvjSfc75pqANG+Dll/32yiv+sXVGUe/e/k5g3Dg46ij/OG6c7z/I27RSKT4B3XWq01ZCySw2YGzGjLbv7drlg/9LL8W2tWvh17+O7dO7t28KGjcOjj4axo/327BhyuCUDMjCQD/V8EVa2b3b3wG8+KLf1q/3W2NjbJ9DDoFjjoFJk/w2ebJvItKsppIP1KQjkqZdu2LB/29/g+ef91u0j6BPHz9oLXoRmDTJ3xlogRzJtow36ZjZacD/A3oCtzjnftzu/d7Ar4DJQDPwRedcQxDHFsmGAQNg6lSY2lAHK/0ttxs2nDd+sJhHhtbyzDPw7LO+GfbGG/3vlJT4pqDoXcCkSf7OIKW0UZEMSLuGb2Y9gY3AKUAj8DRwtnPupVb7XAoc45ybb2ZfAj7nnPtiss9VDV/yTgqdavv3+8Fkzz4b2555xi+MA74TuLq67Z3AhAnQv38OzkeKUkabdMxsCrDIOXdq5Pn3AJxzP2q1z6rIPmvNrBfwFlDhkhxcAV/yTjfT5pzzfQDRu4DoFs0UMvMDx1pfBCZN8lNli3RVppt0hgJbWz1vBI5PtI9zrsXMdgFlQJvZ2s1sHjAPYLimPJB8083pGMxi00p/9rOx1998009PHb0ArF0Ld90Ve3/kyI4XgUMPDeA8JLTyKq/AObcEWAK+hp/j4oi0FfB8TIMH+23mzNhrO3a0vQg8+yzce2/s/aFD22YHTZrkVzFTmmiOFNjo/SAC/jZgWKvnlZHX4u3TGGnSGYDvvBUpHFlYAau8HE45xW9Ru3b5jKBof8Czz8LKlbFR+IceCiecENtqatQnkBVZmOwsaEG04ffCd9rOwAf2p4EvO+debLXPV4HxrTptP++cOyvZ56oNX/JSntTo3nsP/vpXH/yfftrPNrphg3+vRw+fHTRlSiSzaKpvHtJdQMDydAK+jOfhm9lM4Of4tMxbnXOLzexqoN45t8LM+gB3ABOBncCXnHOvJftMBXzJujwJ5t21c6cP/E8+6fsD1q71ax+DnzguGvynToWJE33aqKQhT9fY0MArkc4U4eyZH37oB4o9/nhsi1ZI+/aF446LXQCmTFFWUJeFtYafCQr4klV5+p83aNu2tb0APP+8vzCAHxn8iU/ASSfBiSf6uwJJIk8rCQr4Ip3J09vzTNu9G556qu1FYPdu/97YsT7wRy8AQ4bktKj5KQ+bARXwRToTkhp+Z1pafFromjXwyCPw2GN+vQHwg8Oiwf/EE/0spJJ/FPBFOpOnt+e59uGHvtknegF49FGfJgp+vYDoBeCkk/zAMsk9BXyRVOTh7Xm++fBDnw76yCP+IvDoo7HF5484Aj71Kb998pN+wjnJPgV8EcmIDz/000U/8gisXg0PP+z7AHr2hOOP98H/1FP9YDCtF5AdCvgikhUffADr1sEf/+i3+nrfF15a6lcdi94BVFXluqTFSwFfRHKiudnX/P/4R1i1KrZy2OjRseA/fbqmggiSAr6I5JxzfvqHaPBfs8b3kffqBR//eOwCMGmSFo9PhwK+iOSd99+HJ56INf88+6x/fdAgOPnk2AVA2T9dkyzg98h2YUREwK/3O306/OhHfhbQ7dvhzjvhjDPgL3+Biy/2yVJHHQWXXw5/+APs3ZvrUuOzuaqq/GC9qir/vEAo4ItIXqiogLPPhttu823969fDddf5oH/zzfDpT0NZmV9E5pZb4I03clDI6HiNzZt9G1V0SuSggn6GLyZq0hGRvLd3r2/zX7kSfvvb2CJjkyfDrFl+mzTJx8mMyuSI7IAG/6kNX0SKhnO+9r9ypd/WrvWvHX44fOYzcPrpfgGZfv0ycPBMzrkU0MVEAV9EitaOHfD73/vg/4c/+Ll/+vb1QX/2bF/7D2wt4EzW8AO6mKjTVkSKVnk5nHsu3H03NDXBn/4Ec+f6OYDmzo0t/nLttbBpU5oHW7y4461DUMtcJlobuZtrJsejgC8iReOgg3xK53/+p69wP/cc/OAHvln829/2A76OOQauuspfELrcwFFb69vUR4zwNe8RI4KbYC+TF5MINemISCg0NMD998N99/m0z/37/Vq/n/88fOELfu6fjHf6diaACfzUhi8i0kpTE6xYAf/zP74JaN8+GDoUPvc5mDMHpk0r3NG+asMXEWmlosK37//udz74L1/u1/i95RY/t39lJXz1q/DS/67DjagqyEFW8WjCUhEJtQEDfKtJLXXsr1+INW5hxz+Hc+9NMxmx/3aMSF785s24efMwKNh1ElTDFxGJDHrqsXUz5hwVezbzb+4mPsKeNrvZnj3svGQhjz9emEsdK+CLiCxc2HaEK2AJ+jdL393CtGk+Qedb3/LzAOVpV2gHCvgiItG5GlIxbDjLl8OECXD99X41r7FjfarnK6+02zfPJlpTwBcRSTS4yazt83796PGjxdTW+jl93noLli71Uzj/x39AdbWf0+enP4Wd/5Xhida6ofgCfp5dUUWkACQa9DR/ftJBVoMG+WmcV6/2M3z+7Gd+QZcrroB/fq1jMxF79vjmo1xxzuXlNnnyZNdly5c716+fc/566rd+/fzrIiLJLF/u3IgRzpn5x1TjRpzf27jRuf2+F6DDtt8sc+fgnAPqXYK4WlwDrzI5sZGISHvJpjT+93/3i/q202xl/PibOzjvPBg/PvgihWekbSanLhURaS9ZJXP37rgBf1dJGeVuBy0tvuP3/PPhy18ObkbP8Iy0zcJscyIiByTK7tmyBXbujPvWgJadvPGGz/Dp2RMWLIAhQ/w8/vfck9llHIsr4GdhtjkRkQOSVTKTvFdRAV/7GtTX+8VcvvlNv4j7mWfC4MHwjW9kprjFFfAzOXWpiEh7ySqZKVZAx42Da67xNwWrVsHMmfD225kpblpt+GY2CLgbqAIagLOccx2KamYfAn+LPN3inDujs8/WbJkiUhBaT2k8aJB/bedOX8OfORMefLDL0x0713EIQKoy1mlrZj8Bdjrnfmxm3wUGOue+E2e/3c65g7vy2Qr4IlJQAlqEPF2Z7LSdDdwe+fl24LNpfp6ISGGKMx9PzgdatZNuwD/MOfdm5Oe3gMMS7NfHzOrNbJ2ZJbwomNm8yH71TU1NaRZNRCSLEmXsxEvbzJFOA76Z/dnM1sfZZrfeLzLCK1H70IjILcaXgZ+b2Ufj7eScW+Kcq3HO1VRUVHT1XEREcifZfDx5MsVLpwHfOXeyc+7oONsDwD/MbDBA5HF7gs/YFnl8DVgDTAzsDERE8sHixfF7Wp3Lm2addJt0VgDnR34+H3ig/Q5mNtDMekd+LgemAi+leVwRkfxSW5t4YvyuTL+cQekG/B8Dp5jZq8DJkeeYWY2Z3RLZpxqoN7O/Ag8DP3bOKeCLSPEZMSL+63ky2j+tgO+ca3bOzXDOjY40/eyMvF7vnLs48vMTzrnxzrmPRR6XBVFwEZE28mFq9Dwf7V9cI21FJJzq8mSxkTwf7V9cs2WKSDhpavQDwjNbpkiY5UOTRq4km7VSDlDAFykG+dKkkSuaGj0lCvgixaAAhvVnVJ53luYLBXyRYhD2Jo2udpaGtPmrV64LICIBGD48fqdlmJo0amtTy4a59FK46abYIKlo81f0M4qYavgixUBNGqmpq2sb7KNC0vylgC9SDPI8/ztvLFyY99MfZJICvkihat8ODT7nfP9+/6hg3/E7SjZVcQiav9SGL1KI2q+uFKJ26JTF+47M4tfwzULR/KUavkghCnsaZirifUfxFos1g/nzQ3GhVMAXKURhT8NMRaLvwrm2fR133AG/+EV2y5YjatIRKURKw+xcou8ohPPrRKmGL1KIlIbZOX1HHSjgixQipWF2Tt9RB2rSESlUqY4sDbPo97NwoW/Tj3Zqh/R7U8AXkeKl9NU21KQjIsVL6attKOCLSPFS+mobCvgixSKkU/4mpYVR2lDAFykGYV/xKhGlZrahgC9SDNRWHZ9SM9swl2iq0Byrqalx9fX1uS6GSGHo0SPxpGD792e/PJIzZvaMc64m3nuq4YsUA7VVSwoU8EWKgdqqJQUK+CLFoHVbNUDPnrE2/LB33MoBGmkrUiyiHZEaWSoJqIYvUkyUrSNJKOCLFBONLJUkFPBFiomydSQJBXyRYqJsHUkirYBvZmea2Ytmtt/M4ib6R/Y7zcw2mNkmM/tuOscUkSQ0slSSSDdLZz3weeDmRDuYWU/gBuAUoBF42sxWOOdeSvPYIhKPFkaRBNIK+M65lwHMLNluxwGbnHOvRfa9C5gNKOCLiGRRNtrwhwJbWz1vjLzWgZnNM7N6M6tvamrKQtFERMKj0xq+mf0ZODzOWwudcw8EWRjn3BJgCfjJ04L8bBGRsOs04DvnTk7zGNuAYa2eV0ZeExGRLMpGk87TwGgzG2lmBwFfAlZk4bgiItJKummZnzOzRmAK8DszWxV5fYiZPQjgnGsBLgNWAS8D/+2cezG9YouISFelm6VzH3BfnNffAGa2ev4g8GA6xxIRkfRopK2ISEgo4ItIeurqoKrKL7NYVaX59/OYAr6IdD9o19X5+fY3b/Zr6kbn31fQz0sK+CJhl07Q1vz7BUUBXyTs0gna2Zh/X01GgVHAFwm7dIJ2puffV5NRoBTwRcIunaCd6fn31WQUKAV8kbBLJ2hnev59LdkYKAV8kbBLN2jX1kJDA+zf7x+DnItfSzYGSgFfRDIbtNOhJRsDpYAvIvlLSzYGKt0lDkVEMktLNgZGNXwRkZBQwBcRCQkFfBGRkFDAFxEJCQV8EZGQMOdcrssQl5k1AZtT3L0c2JHB4uQrnXe46LzDpbvnPcI5VxHvjbwN+F1hZvXOuZpclyPbdN7hovMOl0yct5p0RERCQgFfRCQkiiXgL8l1AXJE5x0uOu9wCfy8i6INX0REOlcsNXwREemEAr6ISEgUVMA3s9PMbIOZbTKz78Z5v7eZ3R15/0kzq8p+KYOXwnl/w8xeMrMXzGy1mY3IRTmD1tl5t9rvC2bmzKzgU/dSOWczOyvy937RzO7MdhkzIYV/48PN7GEzey7y73xmLsoZNDO71cy2m9n6BO+bmV0f+V5eMLNJaR3QOVcQG9AT+DtwBHAQ8FfgqHb7XArcFPn5S8DduS53ls57OtAv8vMlYTnvyH79gUeBdUBNrsudhb/1aOA5YGDk+aG5LneWznsJcEnk56OAhlyXO6Bz/1/AJGB9gvdnAr8HDDgBeDKd4xVSDf84YJNz7jXn3AfAXcDsdvvMBm6P/HwPMMPMLItlzIROz9s597BzLrrS8zqgMstlzIRU/t4A/we4BtibzcJlSCrn/BXgBufc2wDOue1ZLmMmpHLeDjgk8vMA4I0sli9jnHOPAjuT7DIb+JXz1gGlZja4u8crpIA/FNja6nlj5LW4+zjnWoBdQFlWSpc5qZx3a3PxNYJC1+l5R25vhznnfpfNgmVQKn/rMcAYM3vczNaZ2WlZKzkIH5YAAAHtSURBVF3mpHLei4BzzKwReBD4WnaKlnNd/f+flFa8KiJmdg5QA5yY67Jkmpn1AK4DLshxUbKtF75Z5yT8ndyjZjbeOfdOTkuVeWcDv3TO/V8zmwLcYWZHO+f257pghaSQavjbgGGtnldGXou7j5n1wt/6NWeldJmTynljZicDC4EznHPvZ6lsmdTZefcHjgbWmFkDvn1zRYF33Kbyt24EVjjn9jnnXgc24i8AhSyV854L/DeAc24t0Ac/uVixS+n/f6oKKeA/DYw2s5FmdhC+U3ZFu31WAOdHfp4DPOQiPR8FrNPzNrOJwM34YF8MbbrQyXk753Y558qdc1XOuSp838UZzrn63BQ3EKn8G78fX7vHzMrxTTyvZbOQGZDKeW8BZgCYWTU+4DdltZS5sQI4L5KtcwKwyzn3Znc/rGCadJxzLWZ2GbAK36t/q3PuRTO7Gqh3zq0AluFv9TbhO0K+lLsSByPF874WOBj4TaSPeotz7oycFToAKZ53UUnxnFcBnzKzl4APgSuccwV9F5vieX8TWGpmC/AduBcUQWUOM/s1/gJeHumf+AFQAuCcuwnfXzET2ATsAS5M63hF8J2JiEgKCqlJR0RE0qCALyISEgr4IiIhoYAvIhISCvgiIiGhgC8iEhIK+CIiIfH/AZksqJQeDMN0AAAAAElFTkSuQmCC\n" }, "metadata": { "needs_background": "light" } } ] } ] }