notebooks/pca-introduction.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "%matplotlib inline\n",
    "plt.rcParams[\"figure.dpi\"] = 150\n",
    "plt.rcParams[\"figure.facecolor\"] = \"white\""
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [],
   "source": [
    "# create some random data\n",
    "x = 10 * np.random.randn(25)\n",
    "y = np.random.randn(25)\n",
    "y2 = np.zeros(25)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAABCQAAAG4CAYAAACdCjaWAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/Z1A+gAAAACXBIWXMAABcSAAAXEgFnn9JSAAA/GUlEQVR4nO3de7hVdb0v/vdC7guUmxcUMcVLRiJgIoiWtslLPqaQO4/lNk093bxt62R7J+pW22W6K/t1yh51ezt7l79oY9soE7M6KpdEQPqhIgGiorsQlwoLuen8/dFZHJHFWmsu5hpzXV6v5/F5ZNzmZ8wx4TPHe47xHTWlUqkUAAAAgAJ1q3YBAAAAQNcjkAAAAAAKJ5AAAAAACieQAAAAAAonkAAAAAAKJ5AAAAAACieQAAAAAAonkAAAAAAKJ5AAAAAACieQAAAAAAonkAAAAAAKJ5AAAAAACieQAAAAAAonkAAAAAAK12UDifXr1+e+++7L+eefn0MOOSS9e/dObW1tDj/88Fx77bVZt27dDte98847M27cuPTr1y+DBg3KRz/60cyaNavA6gGAltDvAaD9qimVSqVqF1ENt912Wy688MIkyaGHHpr3v//9eeONNzJr1qysXbs2733ve/P73/8+e+yxxzbrXXbZZbn55pvTp0+fnHDCCdmwYUN+85vfpFQqZdq0aTn99NOrsDcAQGP0ewBov7psIHHXXXdl1qxZueyyy3LooYdunf7yyy/nlFNOyYIFC3LWWWfl3//937fOe+ihh/KRj3wkgwcPzuzZs3PQQQclSWbPnp3jjjsuffv2zYoVKzJgwICidwcAaIR+DwDtV5cNJJoye/bsHH300enVq1feeOON9OzZM0ny0Y9+NL/61a/yne98J5dddtk261x66aX53ve+l5tuuilf+tKXqlA1AFAO/R4AqqvLjiHRlMMPPzxJsnHjxqxZsyZJ8uabb+bhhx9OkpxxxhnbrdMw7f777y+oSgBgZ+j3AFBdAolGLF++PEnSo0ePDBo0KEmyZMmSbNy4MbvvvnuGDRu23Tpjx45NkixatKi4QgGAVtPvAaC6ule7gPbo5ptvTpKcdNJJ6dWrV5Lk+eefT5JGv5wkSW1tbQYMGJC6urqsXbs2/fv3b/Z1Ro4c2ej0JUuWpE+fPhk+fHhrygeAinv++edTW1ub//qv/6p2KRVTRL/X6wHoKKrR610h8S6//OUvc/vtt6dHjx657rrrtk5veCxY3759d7hubW1tkmTt2rU7VUOpVMrmzZt3ahtA8vqKFXl9xYpqlwGdwubNm1NfX1/tMiqm2v1er4fK0e+hMqrR610h8Q7PPPNMzj777JRKpdx4441b7y1tK4sXL250esOvKTuaD7TMv40fnyT51Jw5Va4EOr4d/dLfERXZ7/V6aHv6PVRGNXq9KyT+j1WrVuWkk05KXV1dLr/88lx66aXbzO/Xr1+SZP369TvcRkOa1JLbNQCA4un3ANB+CCSSvPrqqznhhBOycuXKnHfeebnpppu2W6bhHs8XX3yx0W3U19fntddey8CBA31BAYB2SL8HgPalywcS69aty8knn5ynnnoqU6ZMya233pqamprtljvkkEPSq1evrF69OqtWrdpu/vz585Mko0aNavOaAYDy6PcA0P506UBi48aNOe200/KHP/whJ554Yn784x9nl112aXTZPn365MMf/nCS5Kc//el286dNm5YkOfXUU9uuYACgbPo9ALRPXTaQeOutt3LWWWfl4YcfzrHHHpv/+I//SM+ePZtc5/LLL0+SXH/99Vm6dOnW6bNnz86PfvSjDBgwIOeff36b1g0AtJx+DwDtV5d9ysb3v//9TJ8+PUkyZMiQfOELX2h0uZtuuilDhgxJkkyaNCmXXnppbr755owePTof+chHsmnTpsycOTOlUil33HFHBgwYUNQuAADN0O8BoP3qsoFEXV3d1v9v+KLSmGuuuWbrF5Qk+e53v5vRo0fn+9//fmbOnJmePXtm0qRJmTp1ao4++ug2rRkAKI9+DwDtV02pVCpVuwi25dnkUBmeSw6VozdVlvcTKke/h8qoRm/qsmNIAAAAANUjkAAAAAAKJ5AAAAAACieQAAAAAAonkAAAAAAKJ5AAAAAACieQAAAAAAonkAAAAAAKJ5AAAAAACieQAAAAAAonkAAAAAAKJ5AAAAAACieQAAAAAAonkAAAAAAKJ5AAAAAACieQAAAAAAonkAAAAAAKJ5AAAAAACieQAAAAAAonkAAAAAAKJ5AAAAAACieQAAAAAAonkAAAAAAKJ5AAAAAACieQAAAAAArXpQOJJ554It/85jczZcqUDBs2LDU1Nampqdnh8tdcc83WZRr776tf/WqB1QMAzdHrAaD96l7tAqrpuuuuy89//vOy15s4cWIOPPDA7aYfccQRlSgLAKgQvR4A2q8uHUhMmDAho0aNypFHHpkjjzwy73nPe7Jx48Zm17vgggty7rnntn2BAMBO0esBoP3q0oHEFVdcUe0SAIA2pNcDQPvVpceQAAAAAKqjS18h0VoPP/xwFi5cmA0bNmTYsGE5+eST3VMKAJ2IXg8AbU8g0Qr33HPPNn+eOnVqPv7xj+fOO+9Mv379qlQVAFApej0AtD2BRBkOPPDA3HTTTTn55JOz3377pa6uLv/7f//vfOUrX8nPfvazvPXWW5k+fXqLtzdy5MhGpy9btiwjRoyoVNkAQAvp9QBQHIFEGc4+++xt/lxbW5tPfvKTOf7443PYYYflvvvuy5w5czJ+/PgqVQgA7Ay9HgCKI5CogKFDh+a8887LTTfdlAceeKDFX1IWL17c6PQd/ZoCAFSHXg8AlecpGxVy0EEHJUlefvnlKlcCALQFvR4AKksgUSF1dXVJ/nppJwDQ+ej1AFBZAokKKJVKWwe4Gjt2bJWrAQAqTa8HgMoTSLTQ6tWr8z//5//M2rVrt5m+bt26fP7zn8/cuXOz1157ZcqUKVWqEADYGXo9ABSrSw9qOWPGjFx33XVb/7xp06Yk2WagqqlTp+aUU05JfX19Lrroonz1q1/NkUcemaFDh2b16tWZP39+1qxZkwEDBmTatGnp27dv4fsBADROrweA9qtLBxKrV6/O3Llzt5v+zmmrV69OkgwePDhXXHFF5syZk2effTazZs3KLrvskv333z/nnntu/v7v/z777LNPYbUDAM3T6wGg/aoplUqlahfBthoeBbajR4UBLfNv/+cX0E/NmVPlSqDj05sqy/sJlaPfQ2VUozcZQwIAAAAonEACAAAAKJxAAgAAACicQAIAAAAonEACAAAAKJxAAgAAACicQAIAAAAonEACAAAAKJxAAgAAACicQAIAAAAonEACAAAAKJxAAgAAACicQAIAAAAonEACAAAAKJxAAgAAACicQAIAAAAonEACAAAAKJxAAgAAACicQAIAAAAonEACAAAAKJxAAgAAACicQAIAAAAonEACAAAAKJxAAgAAACicQAIAAAAonEACAAAAKFyXDiSeeOKJfPOb38yUKVMybNiw1NTUpKamptn17rzzzowbNy79+vXLoEGD8tGPfjSzZs0qoGIAoBx6PQC0X92rXUA1XXfddfn5z39e1jqXXXZZbr755vTp0ycnnHBCNmzYkJkzZ+bBBx/MtGnTcvrpp7dNsQBA2fR6AGi/unQgMWHChIwaNSpHHnlkjjzyyLznPe/Jxo0bd7j8Qw89lJtvvjmDBw/O7Nmzc9BBByVJZs+eneOOOy7nnXdejjvuuAwYMKCgPQAAmqLXA0D71aUDiSuuuKKs5b/97W8nSa688sqtX1CSv37Z+dznPpfvfe97uf322/OlL32ponUCAK2j1wNA+9Wlx5Aox5tvvpmHH344SXLGGWdsN79h2v33319oXQBAZej1AFAsgUQLLVmyJBs3bszuu++eYcOGbTd/7NixSZJFixYVXRoAUAF6PQAUq0vfslGO559/Pkka/YKSJLW1tRkwYEDq6uqydu3a9O/fv9ltjhw5stHpy5Yty4gRI1pfLABQNr0eAIrlCokWWrduXZKkb9++O1ymtrY2SbJ27dpCagIAKkevB4BiuUKiihYvXtzo9B39mgIAdCx6PQDsmCskWqhfv35JkvXr1+9wmfr6+iRp0SWcAED7otcDQLEEEi00fPjwJMmLL77Y6Pz6+vq89tprGThwoC8pANAB6fUAUCyBRAsdcsgh6dWrV1avXp1Vq1ZtN3/+/PlJklGjRhVdGgBQAXo9ABRLINFCffr0yYc//OEkyU9/+tPt5k+bNi1JcuqppxZaFwBQGXo9ABRLIFGGyy+/PEly/fXXZ+nSpVunz549Oz/60Y8yYMCAnH/++dUqDwDYSXo9ABSnSz9lY8aMGbnuuuu2/nnTpk1JkvHjx2+dNnXq1JxyyilJkkmTJuXSSy/NzTffnNGjR+cjH/lINm3alJkzZ6ZUKuWOO+7IgAEDCt0HAGDH9HoAaL+6dCCxevXqzJ07d7vp75y2evXqbeZ997vfzejRo/P9738/M2fOTM+ePTNp0qRMnTo1Rx99dJvXDAC0nF4PAO1XTalUKlW7CLbV8GzyHT27HGiZf/s/v4B+as6cKlcCHZ/eVFneT6gc/R4qoxq9yRgSAAAAQOEEEgAAAEDhBBIAAABA4QQSAAAAQOEEEgAAAEDhBBIAAABA4QQSAAAAQOEEEgAAAEDhBBIAAABA4QQSAAAAQOEEEgAAAEDhBBIAAABA4QQSAAAAQOEEEgAAAEDhBBIAAABA4QQSAAAAQOEEEgAAAEDhBBIAAABA4QQSAAAAQOEEEgAAAEDhBBIAAABA4QQSAAAAQOEEEgAAAEDhBBIAAABA4QQSAAAAQOEEEq1w3HHHpaamZof/PfDAA9UuEQDYCXo9ALS97tUuoCP7+Mc/nn79+m03fZ999qlCNQBApen1ANB2BBI74aabbsp73vOeapcBALQRvR4A2o5bNgAAAIDCCSQAAACAwrllYyfcfvvtWbNmTbp165aDDz44p59+eoYPH17tsgCACtHrAaDtCCR2wvXXX7/Nn7/85S9n6tSpmTp1apUqAgAqSa8HgLYjkGiFD37wg7ngggty9NFHZ+jQoXnhhRcybdq0XH/99bnqqquy66675tJLL212OyNHjmx0+rJlyzJixIhKlw0AtJBeDwBtr6ZUKpWqXURn8eCDD+bEE0/MgAED8tJLL6VPnz5NLt/cl5TFixe3RZnQZfzb+PFJkk/NmVPlSqDja+hZXb036fXQ/uj3UBnV6PWukKigE044IR/4wAcyb968zJ07N8cdd1yTy+/oQO/oywsAUF16PQBUjqdsVNhBBx2UJHn55ZerXAkA0Bb0egCoDIFEhdXV1SVJamtrq1wJANAW9HoAqAyBRAWtXr06jzzySJJk7NixVa4GAKg0vR4AKkcgUaZZs2blvvvuy1tvvbXN9Oeeey6TJ09OfX19Pvaxj2XYsGFVqhAA2Bl6PQAUw6CWZXr22Wdz3nnnZa+99srYsWMzYMCArFy5Mk888UQ2bNiQkSNH5tZbb612mQBAK+n1AFAMgUSZjjrqqHz+85/P3Llz8/jjj6euri61tbUZPXp0/vZv/zaf//znm30EGADQfun1AFAMgUSZDj300PzgBz+odhkAQBvR6wGgGMaQAAAAAAonkAAAAAAKJ5AAAAAACieQAAAAAAonkAAAAAAKJ5AAAAAACieQAAAAAAonkAAAAAAKJ5AAAAAACieQAAAAAAonkAAAAAAKJ5AAAAAACieQAAAAAAonkAAAAAAKJ5AAAAAACte92gUAxSmVSpm9fE1+vuClrKnfmMG1vXLamL0z4YDBqampqXZ5AABAFyKQgC7i9fWb89n/NS9zlr+6zfR7572Q8QcMyo/O/kB269ujStUBAABdjVs2oAsolUqNhhEN5ix/NZ/9X/NSKpUKrgwAAOiqBBLQBcxevmaHYUSDOctfbXYZAACAShFIQBfw8wUvtWi5+xasauNKAAAA/soYEtAFrKnfWNHlKq0jDbbZkWoFAID2TCAB7VilTn4H1/aq6HKV1JEG2+xItQIAQHsnkIB2qpInv6eN2Tv3znuh2eVOH7NPq2ptrZYOtvnjC8dX/eqDjlQrO+YKFwCA9kMgAe1QpU9+JxwwOOMPGNTkoJXjDxiU8QcManXNrVHOYJsTRgwuqKrGdaRaO6KmgoIkFQkRXOECANC+CCSgBYr+VbWlJ7+3P7oi5x+zf7M11NTU5Ednf2CHIUfDyVjRvxCXM9hmtU/y75vfsgE/20OtHU1TQcEH9huQpCbzVtZtN6+cEKEtr3Bx1QUAQOsIJKAZ1fhVtaUn6tfPeDoPPf3nFtWwW98e+fGF4zNn+au5b8GqrSdOp4/ZJ+MPGFSVE6f2Pthmg9fXb87Mp//comWrXWtH01xQMG/laztct5wQoa2ucHHVBQBA63nsZyu8+eabueqqq3LwwQend+/e2XvvvfOZz3wmq1Z5ZGJn09JfVUulUkVft5yT2nJqqKmpyYQRg3PDGaNy26ePzA1njMqEEdX7Fbc1g22WSqXMWvZKrpi2KBfc9XiumLYos5a9UvFj8M7X++z/mpe69ZvLrpXmtSQoaEpDiNCctnj0bbX+faA4+j0AtC2BRJk2bNiQD3/4w7nuuuuybt26nHbaadl3331zxx13ZMyYMVm+fHm1S6SCyvlVtZLKPaltixqKcNqYvVu0XMNgm6+v35yzbp2TT946N/fOeyEPPf2X3DvvhXzy1rk569Y5eb2FoUE5yj1hLnpg0I6upUFBU1oSIrTF1TjV+veBYuj3AND23LJRpuuvvz5z5szJhAkT8uCDD6Zfv35Jkm9/+9v50pe+lM985jP53e9+V90iqZhqjXHQ0qdiVLKGatwHX85gm9V6ykU5J8zVGBi0o6vELS4t2UZbPPq2I42BQvn0+2JVsge1dlvVGg/mna/7yroNeevtJDXJLt1qMqQVNRQxSHBTr3vfglX501/WpW795gzs2yMH7t4vp4/dJ+P3H5Q5K17d4fx311CN49HwmtPnv5iFL7yWl1/fmF26JYP69swnxw/P+RP3T7dujf+eW069nWXsoXL3o6jPZlu8v+/e5sC+PdKze7fM+tOarN24Jf1798gnj9q3yc8IjaspVfBa0n/6p3/K+eefn2HDhlVqk+3Kpk2bsscee+T111/P/PnzM2bMmG3mH3744Vm0aFHmzZuXI444otWvM3LkyCTJ4sWLd6pedt4Fdz2eh57+S7PLTTp0j9z26SMr9rqlUiln3TqnrF9Wd6aGHd0Hn6TN74Nv6WvPWvZKPnnr3Ga39+MLx289+fu38eOTJJ+aM6fV9bX0MzCwb4/87svHGy+gTFdMW1R2+PZuZ35g39xwxqgml2nN56c5Lf1s7L1b79z0icM73BfNd+tKvamIft+V3s/mVLIHtXZb1eqDTb1ua2poansf2G9gklKjY/Ps7D62ZD927d09b2zYssP576yh3ONRiX7f0n2YcfGx2Xdw3xav++56q/mdq5LK3Y+mP5sD0tgA1jvaViXr2tltvtuOPiMdRTV6U0Xjm3/6p3/K/vvvn1NPPTX/+Z//mbfffruSm6+6xx57LK+//npGjBix3ZeTJDnjjDOSJPfff3/RpdFG2uJX1ZZoeCpGOb+2t7aGat8H3zDY5o8vHJ8zP7BvJh26R878wL5bpzU0jbYYA6AlWvq+Tjp0jw7xBaK9aeltO01pyW0yDVfjNKXcK1xa+tl46fUNbXpbEZWn3xenkj2otduqVh9s7nXLraH5QYLrdjhQ8M7sY0v3o6kw4p01vP3224Ufj3L24ZT/55FtznHK+fxU+ztXpZS7Hy0ZwLqxMKKxbVWyrpYo5+9p0vhnhKZVNJC4/vrrM3z48MyYMSOTJ0/Ovvvum6lTp+a5556r5MtUzZNPPpkkGTt2bKPzG6YvWrSosJpoW+WOcVBJDSfqV55yaJvW0B7ug28YbPObHz8snzlm/yTJ7Y8uz1d/9setA1ZW64kcLf0MLH5prZPNVmhJUNCUloYIzYV8rXn0bblhSkf5ool+X6RK9qDWbqtafbDcMYqaq6GoQYIr/brvruFfH3uu8ONRzj68sWFL7njsubLWbai3PXznqoRy96Ooz2ZbvL+tqf3dnxGaVtFbNho89NBDue2223Lfffdl06ZN6datWyZNmpQLL7wwp512Wrp375hDV1x++eX5zne+k7//+7/Pt7/97e3mP/nkkxk9enTGjh2bJ554otntNVwS827Lli3LkG7dcsOopi9B7krWbtiSuvWbsuWtUrrvUpOBfXumf+/WfY7K3dbSv6zLuiZS/X69u+egPfq1qpaWassann91fdas29TscoP79czwQW13+dlbb5ey/JX6RvezX+/u6blLt7xaX16dr/x//1+SZMj7379TtTX3/r+zzrb+LHRGTR372l67JEnqN7613bx+vbvngCG12aVbebdBVPLfk5Z+Nt7pwD36tfr1qumKRYuy2/77d4lbDCrZ7/X6plWyB7V2W9Xqgy193ZbW0JrtlbP9HanE675Trx7dsnFz878uV7Lfl7sPvXp0y/uG7lrWuoP79UySdvGda2eV+3emqM9mW/xdbm3t7/yMdCTV6PVt8o1o0qRJmTRpUtasWZO77rort99+ex588MHMnDkzu+++e84999ycf/75Oeigg9ri5dvMunXrkiR9+zb+Aa6trU2SrF27trCaOrsdnaisWbep7JOR1m7rgCG1TZ4oHzCktow9ap22rGHLWy3LJFu6XGvtaP+SZN2GLenTc5cWbWdg356VLCvJX9//pX9Zlzc3bX9S/E7rNmzJ2g1bOuTJZjXt0q0mB+3Rr8mgoJIhQv/e3St2jJr6u7kjdes3+Yy0c/p9cSrZg1q7rWr1wdZsr6l1KlFfpWtqjbfeLv54lLutd9bYFp+ftv7OtbPK3eeiPpvt6Vi09HNMGz9lY/Dgwbn88stz+eWXZ9asWfnBD36Qf//3f8+NN96YG2+8Mccdd1wuuuiiTJ48uS3LaLd2lDw1/JqyMwPzdAYtGdhx/AGDWvRUhZ3dVqlUypzlfx0ZumG03tPH7JPxBwwqbJC6tqrhimmLMq0Fgwqe+YF986lmBg5srVnLXsnXWjDg4PuG7pqnXn5jh/PHHzAoN77jGFZikKsGX5n2ZKbNe7HZ5VoywCKdS8PfzS/9vwvz0usbml1+0qF75B8qOAhuUf55B7/00zS9vmmV7EGt3Va1+mBLX7elNbRme+Vsf0cq8brvtP+Q2qx4pb7Z5d5Z6872+3L3Yf8htfnml48ra90zP7BvklT9O1cllPt3pqjPZlv8XW5t7e/8jHQk1ej1hTyTZNmyZbn//vvzm9/8Zuu0YcOG5be//W3OOOOMjBs3Li+8ULl/yNpKwyO/1q9f3+j8+vq//uPZv3//wmrqzNrDfaUNGsY4uOGMUbnt00fmhjNGZcKIYkfMb6saihono1QqZdayV3LFtEW54K7Hc8W0RVvHh2jpgJWH7Nm/omMAlKMlt4sklR/DoqNq6nh3Ng1/N489aPcWLV/pQXCpPP2+OJXsQa3dVrXGi2rNoL5N1VDUIMFt8brv9KmjhrdouUoej3L34ex31FjO56eaY5NVUrn7UdRnsy3e39bWfnYLP8e0YSCxefPm/OQnP8nf/M3f5OCDD84NN9yQLVu25PLLL88zzzyTlStX5rHHHsvJJ5+cefPm5aKLLmqrUipm+PC/frBefLHxX0kbpu+3336F1dSZVfKpCtV6QkNH0BZPH3i319dvzlm3zsknb52be+e9kIee/kvunffC1icP/Ncbzf+qnCRrN25u0RM52kJRT1zpDCfyzR3vzjr4Z2f5ool+X6RK9qDWbquIPtiYcgf1ba6GogYJrvTrvruGz0x8T+HHo5x92LV395w38T1lrdtQb7U+a5VW7n4U9dlsi/e3NbW/+zNC0yoeSDz99NO5/PLLs/fee+dTn/pUfvvb32bChAm5++678+KLL+amm27aOnbEhAkT8otf/CLjxo3L73//+0qXUnGHH354kmT+/PmNzm+YPqqLD1BVKZV8qkK1ntDQEbTF0wfeqVQq5b/f83iTj2Bq6jaMdxpc26tqV6sUcbLZGU7kO8sjzVpj/P6DsmszY0Ps2rt7jtp/YEEV0Vr6fXEq2YNau6227oOtrbfcGprb3gf2G5gP7Deg1dtv7es2aO7fx4YaunXrVvjxKGcfZlx8bLp1+7+nUOV8fqr1Wau0cvej+c/mgHxgv8Z7YxH/BuzMNt+tsc8ITavoGBLHHHNMZs+enVKplF133TWf//zn87nPfS7vb2bE25EjR+bxxx+vZCltYuLEidltt92ybNmyLFy4MKNHj95m/rRp05Ikp556ahWq63wq+Yt0Ub9ud1QNjxhtizEqZj7958xd0fizpRusXtuyIKiavyw3JOTNjUPS2l8AWnoi35IxU6qpnNujJowYXFBVxZiz4tW80czglm9s2JK5K+o63b53Nvp9sSrZg1q7rbbsg+XU+8q6DXn77SQ1SbduNRlSZg3N7UeSNtnHd77u9AUv5k9/qU/d+k0Z2LdHDty9XyaPHZaj9h+YuSvqdjj/nTVU43i88zX/Y/4LWfjC63n59Q3ZpVsyqG/PnD1+v5w38T2NnmiWU2+1PmuVVu5+FPXZbIv3t7FtDuzbI7177JJHl76StRu3pH/vHjn7qOE7/IywYxV97Ge3bt0yduzYfO5zn8snP/nJHY5O/W6zZ8/Os88+m09/+tOVKqXNXHnllfn617+eo48+Og8++ODWkba//e1v50tf+lI+9KEP5Xe/+91OvUbDQFdd4dFqTZm17JV8sgUDHf74wvHNfrmv5LZouVKplKP++Tf5SwsChz3692pyuZYOYPpOlRzUMvnrFQw7Cg0aUvd33zZSKpUye/ma/HzBS1ub4mlj9s6EA7a9qqOzfEavmLYo97ZwQKnONvhnZ9/3rtab2rrfd7X3E9pSpfs9dFXV6E0VvULi8ccfzxFHHFH2ehMmTMiECRMqWUqbufLKK/PQQw9l1qxZOeigg3Lsscdm5cqVmTt3bnbffff867/+a7VL7DQq+Yt0W/+6TeNmL1/TojAiSQ4d2j8H7F7b5Ml+tX81KDd131GAce+8F7YLMMoZ56Q9BxJd+faorrzvnZF+DwBtr6LXk7QmjOhoevfund/+9reZOnVq+vbtm/vuuy8rV67Mueeem/nz5+eAAw6odomdRnu4r5Sd09KT7CTZa9c+VRuwshwtHcOi3LEUOsvJbFe+Paor73tnpN8DQNur6BUSXUWfPn1y7bXX5tprr612KZ1ee7ivlNYr5+T59DH7bD3Zb89XALRUuWMpdJaT2dPG7N2i2xY645MmuvK+d1b6PQC0LYEE7V4lT1I70wlvR9DSk+fd+/fqdLfLlHsLRmc5me3Kt0d15X0HAGgNQ4ACbaalj8r8xuTDOt0VKuXegtFZnk3elW+P6sr7DgDQGq6QANpMS38x/ptD9yiwqmKUewtGw8lsc0/x6Agns1359qiuvO8AAOUSSABtpjOdZJerNbdgdKaT2a58e1RX3ncAgHIIJIA21ZlOssvR2vEEnMwCANBVCCSANtcVT7K78tUhAADQEgIJgDbSVa8OAQCAlhBIALShrnh1CAAAtITHfgIAAACFE0gAAAAAhRNIAAAAAIUTSAAAAACFE0gAAAAAhRNIAAAAAIUTSAAAAACFE0gAAAAAhRNIAAAAAIUTSAAAAACFE0gAAAAAhRNIAAAAAIUTSAAAAACFE0gAAAAAhRNIAAAAAIUTSAAAAACFE0gAAAAAhRNIlOF3v/tdampqdvjf+PHjq10iALAT9HoAKE73ahfQEY0YMSLHHHNMo9MBgI5PrweAtieQaIVjjjkmd955Z7XLAADaiF4PAG3PLRsAAABA4QQSAAAAQOHcstEKS5cuzT/8wz9kzZo1GTJkSI455picdNJJ6dZNvgMAnYFeDwBtTyDRCrNmzcqsWbO2mXbYYYflZz/7WQ466KAqVQUAVIpeDwBtTyBRht122y3/43/8j3z84x/f+mVk4cKF+drXvpY5c+bkhBNOyMKFC7Pbbru1aHsjR45sdPqyZcuM4g0AVaDXA0BxulQgMXny5Dz99NNlrXP33Xdn3LhxSZIxY8ZkzJgx28z/8Ic/nEcffTTHH398HnnkkfzgBz/IP/zDP1SsZgCg5fR6AOg4ulQgsWLFiixZsqSsddavX9/sMrvsskuuuOKKPPLII/n1r3/d4i8pixcvbnT6jn5NAQCaptcDQMfRpQKJhQsXttm2Gy7rfPnll9vsNQCApun1ANBxGCq6Qurq6pIktbW1Va4EAGgLej0AVJZAokJ+9rOfJUnGjh1b5UoAgLag1wNAZQkkyvDd7343L7zwwjbTSqVSfvSjH+U73/lOampq8vnPf75K1QEAO0uvB4DidKkxJHbWd7/73Xz5y1/O2LFjs//++2fDhg354x//mBUrVqRbt2753ve+lyOOOKLaZQIAraTXA0BxBBJl+NKXvpQHH3wwixcvzlNPPZXNmzdn6NChOfvss3PJJZfkyCOPrHaJAMBO0OsBoDgCiTJcfPHFufjii6tdBgDQRvR6ACiOMSQAAACAwgkkAAAAgMIJJAAAAIDCCSQAAACAwgkkAAAAgMIJJAAAAIDCCSQAAACAwgkkAAAAgMIJJAAAAIDCCSQAAACAwgkkAAAAgMIJJAAAAIDCCSQAAACAwgkkAAAAgMIJJAAAAIDCCSQAAACAwgkkAAAAgMIJJAAAAIDCCSQAAACAwgkkAAAAgMIJJAAAAIDCCSQAAACAwgkkAAAAgMIJJAAAAIDCCSQAAACAwgkkAAAAgMJ12UCivr4+99xzTy6++OIcddRR6dWrV2pqanLNNdc0u+6LL76Y8847L3vvvXd69+6dgw8+OFdffXU2bNjQ9oUDAC2i1wNA+9a92gVUy9KlS3POOeeUvd6f/vSnTJgwIa+88kre//7359hjj828efNy7bXX5je/+U1+85vfpFevXm1QMQBQDr0eANq3LnuFRP/+/XP++efnlltuyRNPPJFrr722Reude+65eeWVV3LJJZfkj3/8Y+69994sWbIkkydPzmOPPZZvfOMbbVw5ANASej0AtG9dNpAYMWJEbrvttnz2s5/N2LFj06NHj2bX+cMf/pDHHnsse+yxR771rW9tnd69e/f88Ic/TI8ePfK9730vW7ZsacvSAYAW0OsBoH3rsoFEa8yYMSNJcuqpp253qeaee+6ZY489NnV1dXn00UerUR4AsJP0egAojkCiDE8++WSSZOzYsY3Ob5i+aNGiwmoCACpHrweA4nTZQS1b4/nnn0+SDBs2rNH5DdNXrlzZou2NHDmy0enLli3LiBEjWlEhALAz9HoAKI4rJMqwbt26JEnfvn0bnV9bW5skWbt2bWE1AQCVo9cDQHE67BUSkydPztNPP13WOnfffXfGjRvXRhWVb/HixY1O39GvKQDQlej1ANC5ddhAYsWKFVmyZElZ66xfv36nXrNfv35Nbqe+vj7JXx8zBgDsHL0eADq3DhtILFy4sPDXHD58eBYsWJAXX3yx0fkN0/fbb78iywKATkmvB4DOzRgSZTj88MOTJPPnz290fsP0UaNGFVYTAFA5ej0AFEcgUYZTTjklSXL//fdn48aN28z785//nEceeSQDBw7MxIkTq1EeALCT9HoAKI5Aogzjxo3LxIkT85e//CVXXHHF1ulbtmzJF77whWzevDmXXHJJevToUcUqAYDW0usBoDg1pVKpVO0iqmXy5Ml5+eWXkyQvvfRSXnjhheyzzz5bnzE+dOjQTJ8+fZt1li5dmgkTJmTNmjU57LDD8r73vS+PP/54li9fnqOPPjoPP/xwevXqtVN1NYy8vaORuYGW+bfx45Mkn5ozp8qVQMfXUXuTXg+dn34PlVGN3tRhB7WshAULFmTlypXbTFu1alVWrVqVpPEBqw466KAsWLAgV111VR544IFMnz49w4cPz9SpU/OP//iPO/0FBQCoHL0eANqvLh1IPPfcc61ab999980dd9xR2WIAgIrT6wGg/TKGBAAAAFA4gQQAAABQOIEEAAAAUDiBBAAAAFA4gQQAAABQOIEEAAAAUDiBBAAAAFA4gQQAAABQOIEEAAAAUDiBBAAAAFA4gQQAAABQOIEEAAAAUDiBBAAAAFA4gQQAAABQOIEEAAAAUDiBBAAAAFA4gQQAAABQOIEEAAAAUDiBBAAAAFA4gQQAAABQOIEEAAAAUDiBBAAAAFA4gQQAAABQOIEEAAAAUDiBBAAAAFC4LhtI1NfX55577snFF1+co446Kr169UpNTU2uueaaJterqalp8r8NGzYUswMAQJP0egBo37pXu4BqWbp0ac4555xWrVtbW5szzjij0Xm77LLLzpQFAFSIXg8A7VuXDST69++f888/P0ceeWSOPPLIzJgxI1dddVWL1h0yZEjuvPPOti0QANgpej0AtG9dNpAYMWJEbrvttq1/fvDBB6tYDQBQaXo9ALRvXXYMCQAAAKB6uuwVEjujvr4+X//61/P888+nb9++GTNmTKZMmZJ+/fpVuzQAoAL0egBoewKJVnjllVdy5ZVXbjPt8ssvz1133ZVTTjmlSlUBAJWi1wNA23PLRpnOOeecPPDAA1m1alXWrVuXBQsW5O/+7u+yZs2aTJkyJY8//niLtzVy5MhG/1u2bFkb7gEA0BS9HgCK0WGvkJg8eXKefvrpsta5++67M27cuJ163bvuumubP48ePTp333139t133/zzP/9zrrzyyvz617/eqdcAAPR6AOjsOmwgsWLFiixZsqSsddavX99G1SRf+cpXcsMNN+R3v/tdNm3alJ49eza7zuLFixudPnLkyEqXBwAdjl4PAJ1bhw0kFi5cWO0StrHbbrtljz32yMsvv5w1a9Zk6NCh1S4JADo0vR4AOjdjSFTI22+/nTfeeCNJUltbW+VqAIBK0+sBoLIEEhXywAMPpL6+PiNGjMiuu+5a7XIAgArT6wGgsgQSZfjJT37S6Mjav//973PhhRcmSb74xS8WXRYAUCF6PQAUp8OOIVEJkydPzssvv5wkeemll5Ikt912Wx544IEkydChQzN9+vStyz/wwAO56667cvDBB2fkyJHp0aNHnn322a33uP63//bfcumllxa7EwDADun1ANB+delAYsGCBVm5cuU201atWpVVq1YlSfbbb79t5p155pnZsmVLnnjiifz2t7/NunXrMmjQoJx88sn5zGc+kzPOOKOw2gGA5un1ANB+delA4rnnnitr+ZNPPjknn3xy2xQDAFScXg8A7ZcxJAAAAIDCCSQAAACAwgkkAAAAgMIJJAAAAIDCCSQAAACAwgkkAAAAgMIJJAAAAIDCCSQAAACAwgkkAAAAgMIJJAAAAIDCCSQAAACAwgkkAAAAgMIJJAAAAIDCCSQAAACAwgkkAAAAgMIJJAAAAIDCCSQAAACAwgkkAAAAgMIJJAAAAIDCCSQAAACAwgkkAAAAgMIJJAAAAIDCCSQAAACAwgkkAAAAgMIJJAAAAIDCCSQAAACAwnXZQOKZZ57JDTfckOOPPz5DhgxJjx49stdee2XKlCl55JFHmlz3xRdfzHnnnZe99947vXv3zsEHH5yrr746GzZsKKh6AKA5ej0AtG9dNpCYNGlSvvrVr2bevHkZM2ZMpkyZkt133z3Tp0/Phz70oXz3u99tdL0//elPGTNmTO68884MHjw4p512Wt56661ce+21mTRpUjZu3FjsjgAAjdLrAaB967KBxHvf+97cfffdWb16dWbOnJl77703f/zjH3PLLbekVCrly1/+cp566qnt1jv33HPzyiuv5JJLLskf//jH3HvvvVmyZEkmT56cxx57LN/4xjeqsDcAwLvp9QDQvnXZQOKhhx7K3/3d36V3797bTP/sZz+bE044IW+99VZ++tOfbjPvD3/4Qx577LHsscce+da3vrV1evfu3fPDH/4wPXr0yPe+971s2bKlkH0AAHZMrweA9q3LBhJNOfzww5MkL7300jbTZ8yYkSQ59dRT06tXr23m7bnnnjn22GNTV1eXRx99tJhCAYBW0esBoPoEEo1Yvnx5kmSvvfbaZvqTTz6ZJBk7dmyj6zVMX7RoURtWBwDsLL0eAKqve7ULaG+WLVuWX/ziF0mSj33sY9vMe/7555Mkw4YNa3TdhukrV65s0WuNHDmy0enPPPNMevToscP5QMu8vmJFkuSf/V2CnbZs2bL06NGj2mVUhF4PnYt+D5VRjV4vkHiHLVu25Nxzz83GjRtz5pln5ogjjthm/rp165Ikffv2bXT92traJMnatWt3qo633347mzdv3qlt0PaWLVuWJBkxYkSVK2FHXnn77STJblWug6b5u9QxbN68uVOMm6DXUw7/PnUM+n375+9Sx1CNXt9hA4nJkyfn6aefLmudu+++O+PGjdvh/EsuuSSPPvpoDjjggPzgBz/Y2RKbtXjx4kanN/xasqP5tA+OU/vnGHUMjlPHUI1f8vV6qs1x6hgcp/bPMeoYqtHrO2wgsWLFiixZsqSsddavX7/DeV//+tfzwx/+MHvuuWd+/etfZ9CgQdst069fvya3U19fnyTp379/WXUBANvT6wGgc+uwgcTChQsrtq1bbrklV155ZXbbbbc88MADOfDAAxtdbvjw4VmwYEFefPHFRuc3TN9vv/0qVhsAdFV6PQB0bl3+KRs/+clP8sUvfjF9+/bNjBkzMnr06B0u2/CIsPnz5zc6v2H6qFGjKl4nANA6ej0AtE9dOpD45S9/mXPOOSfdu3fP9OnTM3HixCaXP+WUU5Ik999/fzZu3LjNvD//+c955JFHMnDgwGa3AwAUQ68HgParywYSjz32WM4444yUSqXce++9OeGEE5pdZ9y4cZk4cWL+8pe/5Iorrtg6fcuWLfnCF76QzZs355JLLuk0j0UDgI5MrweA9q2mVCqVql1ENQwcODCvvfZa9t9//3zwgx9sdJljjjkmF1xwwTbTli5dmgkTJmTNmjU57LDD8r73vS+PP/54li9fnqOPPjoPP/xwevXqVcQuAABN0OsBoH3rsoFETU1Ns8t8+tOfzp133rnd9BdeeCFXXXVVHnjggbz66qsZPnx4zjrrrPzjP/5jevfu3QbVAgDl0usBoH3rsoEEAAAAUD1ddgwJAAAAoHoEEgAAAEDhBBIAAABA4QQSAAAAQOEEEgAAAEDhBBLtxDPPPJMbbrghxx9/fIYMGZIePXpkr732ypQpU/LII480ue6LL76Y8847L3vvvXd69+6dgw8+OFdffXU2bNhQUPVdR319fe65555cfPHFOeqoo9KrV6/U1NTkmmuuaXZdx6k4b775Zq666qocfPDB6d27d/bee+985jOfyapVq6pdWpfyxBNP5Jvf/GamTJmSYcOGpaampkWPYbzzzjszbty49OvXL4MGDcpHP/rRzJo1q4CKu57169fnvvvuy/nnn59DDjkkvXv3Tm1tbQ4//PBce+21Wbdu3Q7XdZzKp9d3DHp9x6DXtw96ffvX7nt9iXZhn332KSUp9evXrzRp0qTSJz7xidL73//+UpJSTU1N6Tvf+U6j6y1durQ0ZMiQUpLS+9///tInPvGJ0gEHHFBKUpo4cWJpw4YNxe5IJ7dgwYJSku3+u/rqq5tcz3EqzptvvlkaP358KUlp6NChpU984hOlcePGlZKUdt9999KyZcuqXWKXcdpppzX696Upl156aSlJqU+fPqXTTjutdOKJJ5a6d+9e2mWXXUrTp08vpvAu5NZbb916XA499NDS3/7t35ZOPPHEUv/+/UtJSu9973tLf/7zn7dbz3FqHb2+Y9Dr2z+9vv3Q69u/9t7rBRLtxN/8zd+U7r777tKbb765zfRbbrmllKS0yy67lBYvXrzdehMnTiwlKV1yySVbp23evLk0efLkFjVPyvOnP/2pdP7555duueWW0hNPPFG69tprW/Q+O07F+drXvlZKUpowYUJp7dq1W6f/y7/8SylJ6UMf+lD1iutivvnNb5amTp1a+s///M/Syy+/XOrVq1eTX1JmzpxZSlIaPHhw6dlnn906fdasWaWePXuWBgwYUKqrqyug8q7jzjvvLP33//7fS0899dQ201966aXSmDFjSklKZ5111jbzHKfW0+s7Br2+/dPr2w+9vv1r771eINEBnHDCCaUkpWuuuWab6XPnzi0lKe2xxx7bpe7/9V//VerRo0dp4MCBpc2bNxdZbpfyjW98o9kvGY5TcTZu3FjabbfdSklK8+fP327+qFGjSklK8+bNq0J1NPcl5eSTTy4lafRX4ksuuaSUpHTTTTe1YYW806xZs0pJSr169Spt3Lhx63THqW3o9e2XXt++6PXtm17fsbSHXm8MiQ7g8MMPT5K89NJL20yfMWNGkuTUU09Nr169tpm355575thjj01dXV0effTRYgqlUY5TcR577LG8/vrrGTFiRMaMGbPd/DPOOCNJcv/99xddGs1488038/DDDyf5v8fpnRy74jX0no0bN2bNmjVJHKe2pNd3bI5TcfT6jksPaX/aQ68XSHQAy5cvT5Lstdde20x/8sknkyRjx45tdL2G6YsWLWrD6miO41Qc73XHtWTJkmzcuDG77757hg0btt18x654Db2nR48eGTRoUBLHqS3p9R2b41Qc73XHpYe0P+2h1wsk2rlly5blF7/4RZLkYx/72Dbznn/++SRp9IPyzukrV65swwppjuNUHO91x9Xcsautrc2AAQNSV1eXtWvXFllal3XzzTcnSU466aStv/g6Tm1Dr+/4HKfieK87Lj2k/WkPvV4g0Y5t2bIl5557bjZu3JgzzzwzRxxxxDbzGx7R0rdv30bXr62tTRJ/oavMcSqO97rjau7YJY5fkX75y1/m9ttvT48ePXLddddtne44VZ5e3zk4TsXxXndcekj70l56ffed3gJJksmTJ+fpp58ua527774748aN2+H8Sy65JI8++mgOOOCA/OAHP9jZEknbHCeAzuSZZ57J2WefnVKplBtvvHHr/aXo9R2FXg/QtPbU6wUSFbJixYosWbKkrHXWr1+/w3lf//rX88Mf/jB77rlnfv3rX2+9p+ed+vXr1+R26uvrkyT9+/cvq67OrNLHqSUcp+J4rzuu5o5d4vgVYdWqVTnppJNSV1eXyy+/PJdeeuk287v6cdLrOwa9vnPzXndcXb2HtBftrdcLJCpk4cKFFdvWLbfckiuvvDK77bZbHnjggRx44IGNLjd8+PAsWLAgL774YqPzG6bvt99+Fauto6vkcWopx6k4w4cPTxLvdQfU3LGrr6/Pa6+9loEDB/qS0kZeffXVnHDCCVm5cmXOO++83HTTTdst09WPk17fMej1nZte33F19R7SHrTHXm8MiXbmJz/5Sb74xS+mb9++mTFjRkaPHr3DZRsurZk/f36j8xumjxo1quJ10nKOU3G81x3XIYcckl69emX16tVZtWrVdvMdu7a1bt26nHzyyXnqqacyZcqU3HrrrampqdluOcepMvT6zsdxKo73uuPSQ6qrvfZ6gUQ78stf/jLnnHNOunfvnunTp2fixIlNLn/KKack+eszYDdu3LjNvD//+c955JFHMnDgwGa3Q9tynIozceLE7Lbbblm2bFmjv5BNmzYtyV+fE0/70qdPn3z4wx9Okvz0pz/dbr5j13Y2btyY0047LX/4wx9y4okn5sc//nF22WWXRpd1nHaeXt85OU7F0es7Lj2ketp1ry/RLjz66KOlPn36lLp3716aPn16i9ebOHFiKUnp0ksv3Tpt8+bNpSlTppSSlK6++uqK18r/9Y1vfKNF77PjVJyvfe1rpSSlo48+urRu3bqt0//lX/6llKT0oQ99qHrFdXG9evUqNdV2Zs6cWUpSGjx4cOnZZ5/dOn3WrFmlXr16lQYMGFCqq6sroNKuY8uWLaXJkyeXkpSOPfbYUn19fbPrOE6tp9d3THp9+6PXt196ffvT3nt9TalUKlUm2mBnDBw4MK+99lr233//fPCDH2x0mWOOOSYXXHDBNtOWLl2aCRMmZM2aNTnssMPyvve9L48//niWL1+eo48+Og8//PDWZ8pSGZMnT87LL7+cJHnppZfywgsvZJ999tn6rN6hQ4dm+vTp26zjOBVnw4YNOe644zJ37twMHTo0xx57bFauXJm5c+dm9913z5w5c3LAAQdUu8wuYcaMGds8RuoPf/hDSqVSjjrqqK3Tpk6duvWXxSS57LLLcvPNN6dv3775yEc+kk2bNmXmzJkplUqZNm1aTj/99CJ3odO7+eabc9lllyX5679tu+66a6PL3XTTTRkyZMjWPztOraPXdxx6ffum17cfen371+57fUViDXZakmb/+/SnP93ous8//3zp3HPPLe21116lnj17lg488MDS1KlTS2+++WaxO9FF7Lfffk0ep/3226/R9Ryn4qxfv740derU0ogRI0o9e/Ys7bXXXqVzzz239MILL1S7tC7ljjvuaPbftTvuuKPR9Y444ohS3759SwMGDCiddNJJpccee6z4HegCrr766hb1nxUrVmy3ruNUPr2+49Dr2z+9vn3Q69u/9t7rXSEBAAAAFM6glgAAAEDhBBIAAABA4QQSAAAAQOEEEgAAAEDhBBIAAABA4QQSAAAAQOEEEgAAAEDhBBIAAABA4QQSAAAAQOEEEgAAAEDhBBIAAABA4QQSAAAAQOEEEgAAAEDhBBIAAABA4QQSAAAAQOEEEgAAAEDhBBJAp3LmmWempqYmX/nKV7ab9+yzz6Zfv37p169fli5dWoXqAICdpddD51FTKpVK1S4CoFLq6uoyatSovPTSS3nooYdy/PHHJ0k2b96co48+OvPmzcutt96aCy64oMqVAgCtoddD5+EKCaBTGThwYO6+++4kyTnnnJO6urokyTXXXJN58+bl9NNP9wUFADowvR46D1dIAJ3SV77yldx44435xCc+kYsuuijHHXdc9txzzyxatChDhgypdnkAwE7S66HjE0gAndKmTZty1FFHZeHChdl1112zdu3a/OpXv8qJJ55Y7dIAgArQ66Hjc8sG0Cn17Nkzd911V5LkjTfeyOc+9zlfUACgE9HroeMTSACd1r333rv1/xcuXJi33nqritUAAJWm10PHJpAAOqVHH300N9xwQ/baa69MmjQps2fPzte//vVqlwUAVIheDx2fMSSATueNN97I4Ycfnueeey6/+tWvMmbMmBx22GGpq6vLo48+mqOOOqraJQIAO0Gvh87BFRJAp3PRRRflueeey0UXXZSTTjope+65Z2677bZs2bIlZ599durr66tdIgCwE/R66BwEEkCn8tOf/jT33HNP3ve+9+Vb3/rW1ukf+9jHcuGFF+ZPf/pTLr300ipWCADsDL0eOg+3bACdxqpVq3LYYYelvr4+c+fOzejRo7eZX19fnzFjxmTp0qX5j//4j0yePLk6hQIAraLXQ+cikAAAAAAK55YNAAAAoHACCQAAAKBwAgkAAACgcAIJAAAAoHACCQAAAKBwAgkAAACgcAIJAAAAoHACCQAAAKBwAgkAAACgcAIJAAAAoHACCQAAAKBwAgkAAACgcAIJAAAAoHACCQAAAKBwAgkAAACgcAIJAAAAoHACCQAAAKBw/z9YRUQIhwlngAAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 1200x450 with 2 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.figure(figsize = (8, 3))\n",
    "\n",
    "plt.subplot(121)\n",
    "plt.plot(x, y, \"o\", markersize=4)\n",
    "plt.axvline(c=\"darkred\", lw=0.75)\n",
    "plt.axhline(c=\"darkred\", lw=0.75)\n",
    "plt.xlabel(\"x\")\n",
    "plt.ylabel(\"y\")\n",
    "plt.xlim(-20, 20)\n",
    "plt.ylim(-20, 20)\n",
    "\n",
    "plt.subplot(122)\n",
    "plt.plot(x, y2, \"o\", markersize=4)\n",
    "plt.axvline(c=\"darkred\", lw=0.75)\n",
    "plt.axhline(c=\"darkred\", lw=0.75)\n",
    "plt.xlabel(\"x\")\n",
    "plt.xlim(-20, 20)\n",
    "plt.ylim(-20, 20)\n",
    "\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import numpy as np\n",
    "from sklearn import datasets\n",
    "import matplotlib.pyplot as plt\n",
    "%matplotlib inline\n",
    "plt.rcParams[\"figure.dpi\"] = 150\n",
    "plt.rcParams[\"figure.facecolor\"] = \"white\"\n",
    "\n",
    "iris = datasets.load_iris()\n",
    "data = iris.data\n",
    "\n",
    "# data[:4, :]\n",
    "# iris.feature_names\n",
    "\n",
    "sepal_length = data[:, 0]\n",
    "petal_length = data[:, 2]\n",
    "\n",
    "sepal_length_centered = sepal_length - np.mean(sepal_length)\n",
    "petal_length_centered = petal_length - np.mean(petal_length)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAA/gAAAGwCAYAAAAdT86/AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/Z1A+gAAAACXBIWXMAABcSAAAXEgFnn9JSAAB9KklEQVR4nO3deXwU9f3H8fcmJATCIQkghCskoEIqEFCSiEjEgngRQKvlkENErHc96k+riCC1VdRapYqKoCCIVU5RqlZEBBIEE7BclYRLQDnCGUggyfz+oJsSc80kM5ndzev5eOTxgNnvzvczE8g7392Zz3oMwzAEAAAAAAD8WpDbBQAAAAAAgKpjgQ8AAAAAQABggQ8AAAAAQABggQ8AAAAAQABggQ8AAAAAQABggQ8AAAAAQABggQ8AAAAAQABggQ8AAAAAQABggQ8AAAAAQABggQ8AAAAAQABggQ8AAAAAQABggQ8AAAAAQACo5XYBvqpZs2bKyclR69at3S4FAABJ0q5duxQeHq6ffvrJ7VICBnkPAPAlVc163sEvQ05Ojs6cOeN2GQAAFDlz5oxycnLcLiOgkPcIZEe3b9fR7dvdLgOABVXNet7BL4P3lfyNGze6XAkAAGfFxcW5XULAIe8RyN5LTJQkDU1NdbkSAGZVNet5Bx8AAAAAgADglwv8r776Sh6Pp8KvCRMmuF0qAAAAAADVwi8v0W/WrJlGjBhR6mMFBQWaNWuWJKlnz57VWRYAAAAAAK7xywX+RRddpBkzZpT62KeffqpZs2apVatWSk5Orta6AAAAAABwi19eol8e77v3Q4cOlcfjcbkaAAAAAACqR0At8HNycrRw4UJJ0q233upyNQAAAAAAVJ+AWuDPmzdPOTk5io+PV8eOHd0uBwAAAACAahNQC3zv5fm8ew8AAAAAqGn8ssleafbt26d//etfCg4O1uDBg00/Ly4urtTtmZmZio2Ntas8AAAAAAAcFTAL/Dlz5qigoED9+vVTs2bN3C4HAOAnDMPQ6qxDWpi+V4dy8hQZXlsp8VFKiomkWSsAAAGipuR9wCzwK3t5/saNG0vdXtY7+wCAwHH05BmNnbVWqVnZxbbPXbtbiTERmjrsEjWsG+JSdQAAwA41Ke8D4h78zZs3Kz09XfXq1dOAAQPcLgcA4AcMwyg17L1Ss7I1dtZaGYZRzZUBAAC71LS8D4gF/syZMyVJgwYNUt26dV2uBgDgD1ZnHSoz7L1Ss7IrHAMAAHxXTct7v1/gG4ah2bNnS6J7PgDAvIXpe02NW5C+x+FKAACAU2pa3vv9An/FihXauXOnWrRood69e7tdDgDATxzKybN1HOz34osvatCgQWrfvr0aNmyo2rVrq02bNho+fLi+//57t8sDAPiBmpb3fr/A9zbXGzJkiIKC/P5wAADVJDK8tq3jYL8//elP+vTTTxUREaGrrrpK1113ncLCwjRz5kx169ZNH3/8sdslAgB8XE3Le79eEefl5enDDz+UJA0bNszlagAA/iQlPsrUuAHxLRyuBGVZuHChDh8+rLS0NM2bN0/z5s3T1q1bNWXKFJ05c0a333678vPz3S4TAODDalre+/UCv3bt2srOzpZhGOrUqZPb5QAA/EhSTKQSYyLKHZMYE1HhGDinR48eCgsLK7H9rrvuUmxsrH7++Wdt2rTJhcoAAP6ipuW9Xy/wAQCoLI/Ho6nDLikz0L2fi+vxeKq5MpgREnL284pDQ0NdrgQA4MtqWt6zwAcA1FgN64bovdHdNSyhterVrqXQYI/q1a6lYQmt9d7o7mpYN8TtElGKmTNnauvWrWrfvr3at2/vdjkAAB9Xk/K+ltsFAADglt2HTuq6V1boWO7/7uM+XZCvWWm7tGj9Xi25t6daRdZ1sUJI0vPPP6+NGzcqJydHmzdv1saNGxUVFaU5c+YoODjY1D7i4uJK3Z6ZmanY2Fg7ywUA+JialPe8gw8AqJEKCwtLhP25juXm67pXVqiwsLCaK8Mv/fOf/9Q777yjDz/8UBs3blSbNm00Z84cdevWze3SAAA+rqblPQt8AECNNG3l9jLD3utYbr6mr9xRPQWhTF988YUMw9Dhw4f19ddfq3379urVq5cmTZpkeh8bN24s9Yt37wEgsNW0vGeBDwCokWan7TY1blbaLocrgVnnnXeeevbsqU8++UTdunXTk08+qW+//dbtsgAAPqym5T334AMAApJhGFqddUgL0/fqUE6eIsNrKyU+SkkxkfJ4PDqee8bUfsyOQ/UJCQnRLbfconXr1mnx4sW69NJL3S4JAOCCirJeMp/jgZL3LPABAAHn6MkzGjtrrVKzsottn7t2d9HH4dQPC9HBE6cr3Ff9sMDprBtIGjduLEk6cOCAy5UAANxgJusb1g2pcXnPJfoAgIBiGEapge+VmpWtsbPWakj3Vqb2NyyhtZ3lwSbLly+XJO6hB4AayGzWG4ahIQk1K+9Z4AMAAsrqrENlBr5Xala2OjZvoAZh5V/I1iCslkb1iLaxOpi1cuVKLV26tERX4zNnzuiVV17RzJkzVadOHd1yyy0uVQgAcIvZrE/NytboHm1rVN6zwAcABJSF6XtNjVu0fp+W3NuzzNBvEFZLS+7tqaAgotINP/zwg6655hqdf/756tevn4YOHaqrr75abdq00X333afQ0FDNmDFDrVqZe2cGABA4zGb9gvQ9CgoKqlF5zz34AICAcignz/S4VpF1lTGuj6av3KFZabt0PPeM6oeFaFhCa43qER0wYe+PevXqpccff1zLly/Xhg0bdPDgQYWGhio6Olo33XST7rvvPrVr187tMgEALrCS9ZJqVN6zwAcABJTI8NqWxgUFBWl0zxiN7hnjZFmwqG3btpY+5x4AUHNYzXqp5uR94LxUAQCApJT4KFPjBsS3cLgSAADgBLK+bCzwAQABJSkmUokxEeWOSYyJqHAMAADwTWR92VjgAwACisfj0dRhl5QZ6t7PxvV4PNVcGQAAsANZXzbuwQcABJyGdUM0Z0yiUrOytSB9jw7l5CkyvLYGxLdQYkxEjQx8AAACCVlfOhb4AICA5PF4lBQbqaTYSFfmNwxDq7MOaWH63qJfOlLio5QUE1ljf+kAAMBOZH1JLPABALDZ0ZNnNHbWWqVmZRfbPnft7qLLBhvWDXGpOgAAUFW+mvXcgw8AgI0Mwyg18L1Ss7I1dtZaGYZRzZUBAAA7+HLWs8AHAMBGq7MOlRn4XqlZ2RWOAQAAvsmXs54FPgAANlqYvtfUuAXpexyuBAAAOMGXs5578AEAjvDFxjPV4VBOnq3jAADwVWS9PePsxAIfAGA7X208Ux0iw2vbOg4AAF9E1ts3zk5cog8AsJUvN56pDinxUabGDYhv4XAlAAA4g6z33axngQ8AsJUvN56pDkkxkUqMiSh3TGJMRIVjAADwVWS972Y9C3wAgK18ufFMdfB4PJo67JIyQ9172WIg35sIAAhsZL3vZj334AMAbOXLjWeqS8O6IZozJlGpWdlakL6nqPHQgPgWSoyJYHEPAPBrZL3vZj0LfACArXy58Ux18ng8SoqNVFJspNulAABgK7L+LF/Mei7RBwDYypcbzwAAgKoj630XC3wAgK18ufEMAACoOrLed7HABwDYypcbzwAAgKoj630X9+ADAGznq41nAACAPch638QCHwBgmmEYWp11SAvT9xYFeUp8lJJiIksEudnGM1b26QS35wcAwJc4kfVW9+sEt+evLizwAQCmHD15RmNnrVVqVnax7XPX7i66FK9h3RDX9+lP8wMA4EucykW389bt+asT9+ADACpkGEapweiVmpWtsbPWyjAMV/dphdvzAwDgS5zKRbfz1u35q5vfL/APHDighx9+WBdeeKHq1KmjiIgIde3aVY888ojbpQFAwFiddajMYPRKzcqucIzT+7TC7fkBAPAlTuWi23nr9vzVza8X+OvWrVOHDh30wgsvKCQkRCkpKUpMTFR2drZeeuklt8sDgICxMH2vqXEL0ve4uk8r3J4fAABf4lQuup23bs9f3fz2HvwDBw6oX79+OnXqlBYuXKj+/fsXe3zNmjUuVQbASTWlQYqvOZSTZ+s4p/ZphdvzAwBKR9a7w6lcdDtv3Z6/uvntAv+pp57SwYMHNWXKlBKLe0nq3r27C1UBcFJNapDiayLDa9s6zql9WuH2/ACAksh69ziVi27nrdvzVze/vET/1KlTmjVrlsLDwzVq1Ci3ywFQDWpagxRf079Lc1PjUrpEmd5nSry5sQPiW5jepxVuzw8AKI6sd5cTWS+5n7duz1/d/HKBv3btWh0/flzx8fGqU6eOPv30Uz344IO666679Ne//lV795q7zwKA/6hpDVJ8jUfmLom0cuVkUkykEmMiyh2TGBNR4ZjKcnt+AEBxZL27nMh6yf28dXv+6uaXC/xNmzZJkpo2baoBAwbo2muv1UsvvaTXXntNv//979WuXTvNmTPH1L7i4uJK/crMzHTyEABYVNMapPiahRlmz7/5F1g9Ho+mDrukzED1Xorp1P2Wbs8PACiOrHeXE1kvuZ+3bs9f3fzyHvzDhw9LkhYtWqTg4GBNmTJFv/nNb3Ty5Em9+uqrmjx5skaMGKEOHTqoS5cu7hYLwBaB3CDF7WZCZuZ36vw3rBuiOWMSlZqVrQXpe4rmHxDfQokxEY4fv9vzAwD+J5CzXvL9vHfy/Ludt27PX538coFfWFgoScrPz9ekSZN01113FT32/PPPa+fOnfrHP/6h559/Xu+99165+9q4cWOp2+Pi4uwrGECVBWqDFLebCZmd38nz7/F4lBQbqaTYSMvPtYPb8wMAzgrUrJf8I++dPv9u563b81cXv7xEv169ekV/Lq3Jnnfb8uXLq60mAM4KxAYpbjcTsjJ/IJ5/AIBvCdSs8Ze8N9s8z9/Of03jlwv8Nm3aSJLq1q2rJk2alHg8OjpakrR///7qLAuAgwKxQYrbzYSszB+I5x8A4FsCNWv8Je8lBeT5r2n8coEfHx8v6ezH5eXllbwHJDv77D/Qc9/pB+DfArFBitvNhKzMH4jnHwDgWwI1a/wl7xdm7A3I81/T+OU9+K1bt1bnzp21fv16LV++XH379i32uPfSfO8LAQACg1MNUpxqeuNmMxs79+sd17BuiGbfnqC3V+7Qe2m7dDz3jOqHhWhoQmvd1iNaQUF++ZoxfNTJkyf12WefafHixfrmm2+0c+dOBQcHq127drrxxhv14IMP8kI+EICcbIbmRN672ajWLCvzk/X+zy8X+JL0hz/8QUOHDtXDDz+sf/7zn2revLkkKSMjQy+88IIk6c4773SzRAAOsLtBilNNb3yhmY1d+/WOK+2YDp44rWeWbNYXm392vEEQapbZs2drzJgxkqQOHTqof//+OnbsmFatWqWnnnpKc+bM0fLly9W0aVOXKwVgNyeaoTmR977QqNbO/UaG1ybrA4DfvgQzZMgQjRgxQt9//706duyo6667Tr1791ZiYqKys7M1ZswY/eY3v3G7TAA+zKmmN/7SzMZKMyO3GwSh5gkJCdEdd9yhTZs2adOmTfrggw+0dOlSbd26VfHx8dqyZYseeOABt8sE4AecyDB/alRrdv6ULlFkfQDw2wW+JE2fPl1vvPGGYmNj9dVXX2nNmjXq2rWrZsyYoTfeeMPt8gD4OKea3vhLMxsrzYzcbhCEmmfEiBGaOnWqOnToUGx78+bNNWXKFEnSvHnzdPr0aTfKA+BHnMgwf2pUa3Z+SWR9APDrBb7H49GYMWO0du1a5eTk6MSJE1q1apVGjBjhdmkA/IBTTW/8pZmNlWZGbjcIAs7VuXNnSVJeXp4OHTrkcjUAfJ0TGeZPjWrNzr8wg6wPBH57Dz4AVJVTTW+sNrNxqpmQGWbnd7tBEHCurKwsSWcv44+I4OOaAJTPiQyrTKNaX897sj4wsMAHUGM51fTG6n6tNBNyogOwmfndbhAEnOvll1+WJPXr10+1a/NvDkD5nMiwyuzTbN479ek+Fc1P1gcGFvgAaqyU+CjNXbu7wnFWm944tV+nOv6b4dQxAVZ98sknmjZtmkJCQjRx4kTTz4uLiyt1e2ZmpmJjY+0qD4APciLDyHr4Kr++Bx8AqsKppjcJ0Y0UHFT+K+zBQR51jz7P9D7d7mLvdoMgQJK2bNmiYcOGyTAMPf/880X34gNAeZzIMLIevooFPoAay6mmN2+v2qGCwvLDt6DQ0IxVO03v0+0u9m43CAL27Nmjfv366fDhw3rwwQd1//33W3r+xo0bS/3i3Xsg8DmRYWQ9fBWX6AOo0ZxoejM7reLL2yRpVtouje4ZY2qslW69Zu7lrwy3GwSh5srOzlbfvn21c+dOjRo1SpMnT3a7JAB+xu4MI+vhq1jgA3CME01iCgsLNW3lds1O263juWdUPyxEQxJaaXSPtgoKCqrUWCtN7sw4nnvG1nGS73Sxt/tcARU5ceKErrnmGm3atEmDBg3Sm2++yS+YgA/xl6yX7M0wsh6+yrYF/k8//aSlS5dq+fLlysjI0IEDB3TkyBGdd955atKkibp06aJevXqpX79+atasmV3TAvBRTjSJ2X3opK57ZYWO5eYXbTt44rQmLdmiV/61TUvu7alWkXUtj7Vb/bAQHTxx2tQ4s+hsi5ooLy9PKSkpWrNmja6++mrNmTNHwcHBbpcF4L/IerIevqfK9+D/85//1IABA9S6dWuNHj1a77zzjtavX6+9e/fq5MmT2rt3r9avX6933nlHo0ePVuvWrTVw4EB99tlndtQPwAc50SSmsLCwRIif61huvq57ZYUKCwstjXXC4O6tTI0b2r216X2mxEeZGkdnWwSKgoICDR48WF9++aV69uypefPmKTQ01O2yAPwXWU/WwzdV+h38tWvX6uGHH9aKFSskSZdffrmSk5PVvXt3XXTRRYqIiFCDBg109OhRHT58WJs3b9aaNWu0bNkyLVq0SIsWLdIVV1yhyZMnq1u3brYdEAD3WWkSY/byr2krt5cZ4l7HcvM1feUOFcowPdbsfXFWxDVvYGpcx6j6pvfp7Wxb3nmlsy0Cyauvvqr58+dLkho3bqy77rqr1HGTJ09W48aNq7M0ACLryXr4qkov8Lt3767zzz9ff/rTnzR06FC1bNmy1HERERGKiIhQbGysrr/+eknS7t27NWvWLP3tb39T9+7dVVBQUNkyAPggJ5rEWGlmY5aVxjdWLFq/z9S4hRn7dFm7JqbGejvblvVuCZ1tEWgOHz5c9GfvQr8048ePZ4EPuICsJ+vhmyq9wH/ppZc0duxYhYWFWX5uq1at9Nhjj+n3v/+9pk6dWtkSAPgoJ5rEONHMxjvW7gZBVo/f7Px0tkVNMn78eI0fP97tMgCUwd+yXrI37ytz/GbmJ+tRVZVe4Fv9/NnShIWF2bIfAL7FiSYxVpvZmB3rRIMgK8dvdX462wIAfIE/Zb1kf0NAq8dvZX6yHlVR5SZ7APBLTjSJGZJgrpnNsITWGty99FuGfmlo91a2NwiSrBx/lCPzAwDgNH/KeicaAlo5fifmB8rCAh+A7bxNYspjtUnM6B5t1SCs/IuOGoTV0qge0Ypr3tDUPgtlmG4QZIXZ4y80nJkfAACn+UvWd4xqYKkhoFlWjt+J+YGy2LrAP3DggJ5++mn17t1bHTp0UExMTKlfsbGxdk4LwMd4m8SUFXyVaRITFBSkJff2LDP4G4TV0pJ7eyooKMh045s5a340NW5B+h7TdUrmj39Rhrk6rc4PAIDT/CXrF2bss9QQ0Cwrx+/E/EBZKn0P/i99//336t27t7Kzs7m8BIClJjFmm960iqyrjHF9NH3lDs1K26XjuWdUPyxEwxJaa1SPaAUFnX3N0mzjG7NNeqw0CPIyc/xONegBAKA6WG0IZybD7M56KxluNe/NHr9TzXeB0ti2wL///vt16NAh3XrrrXr44YcVExOj8PBwu3YPwA+ZaRJjtelNUFCQRveMKfcjb8w2vjHbzMdKg6BzVXT8TjboAQCgOphtCGclw+zMeisZXpm8N3P8TjbfBX7Jtkv009LS1KlTJ73zzju6+OKLWdwDqJBTTWfMNr4ZmtDa1DgrDYKsoEEPAKAmcLvJnRMNAa0wO39KF5rvoupsW+DXq1dPHTt2tGt3AGoAp5rOmG18c1uPaNsbBFlBgx4AQE3gdpM7JxoCWmF2fklkParMtgV+7969tX79ert2B6AGcKrpjNnGN0FBQbY3CHKiThr0AAD8mdtN7pxoCGiF2fkXZpD1qDrb7sF/5plnlJSUpCeeeEJPP/20goOD7do1gADlRIMcL7ONb6w2CLKbUw16AADwFU5lmJUM94e8J+thB9sW+LGxsVq1apVSUlI0d+5cJScnq0WLFkWdLs/l8Xj05JNP2jU1AD/lRIOccxmGoX/vPaI1O7KLuvC2Oz9cCW0bFQtysw2CnGJ3gx4AAHyJkxlmNusl3897sh52sG2Bf+bMGU2aNElbtmyRYRjKzMwscywLfADS2aYzc9furnBcZZre7D50Ute9skLHcvOLth08cVqTlmzRK//apiX39lSryLqW9+sWJ88VAABOcirDyHqgJNsW+E888YTeeecdnX/++RoyZIhiYmJUr149u3YPIAB5m86U1yymMk1vCgsLSwT+uY7l5uu6V1YoY1yfUq8y8kVOnSsAAJzmRIaR9UDpbFvgv/fee2rSpInWr1+vpk2b2rVbAAHM23SmrI+EqWzTm2krt5cZ+F7HcvM1feWOcj9j15c4da4AAHCaExlG1gOls22Bf/jwYV1zzTUs7gFY4kTTm9lpFV/eJkmz0nb5TehL7jcIAgCgsuzOMLIeKJ1tC/y4uDgdP37crt0BqEHsbnpzPPeMpXGGYWh11iEtTN9bFKQp8VFKiomsUpA6sV+3GwQBAFBZdmaY1ayXnMllsh6+xrYF/kMPPaRbb71V6enpio+Pt2u3AGBZ/bAQHTxx2tS4oyfPlHop3Ny1u4suhWtYN8RyDU7tFwAAWMt6yZlcJuvhi2zrOJGUlKR77rlHycnJGj9+vFauXKkdO3Zo165dpX4BgFOGJLQyNW5o91Zl3ucmSalZ2Ro7a60Mw7A0v2EYjuwXAACcZTbrhyW0diSXyXr4KtsW+NHR0Xr55Zd1/PhxTZw4UVdccYViY2PVtm3bEl8xMf5zHwwA/zO6R1s1CCv/AqUGYbXUoXmDcjvVSmcDuqIxv7Q665Aj+wUAAGeZzfpRPaIdyWWyHr7Ktkv0r7jiCpo+APAJQUFBWnJvzzI/PqdBWC0tubenXl22zdT+FqTvsXQf3ML0vY7sFwAAnGU264OCghzJZbIevsq2Bf5XX31l164A2KSgoEDjFm3Uwox9Op1foNBawUrp0lwT+scpODi4aJyVBjFONaQzy+z8rSLrKmNcH01fuUOz0nbpeO4Z1Q8L0bCE1hrVI1pBQUE6lJNnak6z46yOt7pfyf3zDwDwLeVl/S+ZzRC3s8bOrJecyWWyHr7KtgU+AN+yac9R3TBlpQoK/3fv1+mCfL2Xtlvvf/ujFt/dQx1bNLTUIMbtZjJW5w8KCtLonjFlfjxOZHhtU/OaHef0ft0+/wAA31JR1j92ukB1Qs++oG82Q9zOGruzXnIml8l6+Crb7sEH4DsKCgpKBH6xxwsN3TBlpfLz8003iHG7mYwT86d0iTI1bkC8uXFF+zU5fkB8C9P7dPv8AwB8i5ms3/rzcUsZXlhYGHBZLzmT92Q9fJXt7+AvX75cy5cv1759+5SXV/olKR6PR9OmTavSPMnJyVq+fHmZj3/66afq169fleYA/NW4RRvLDHyvgkJDv3sv3XSDGEOG6bFO3GtmpZmN2fkNmQtIqzmaFBOpxJiIcutNjIlQYkyE6X06cfwAAP9lJusNQ/rxyCnTGfL2yh0Bl/WSM3lP1sNX2bbAz87O1sCBA/XNN99U+KqSHQt8rxtvvFH16tUrsb1FC/OvlgGBZmHGPlPjlm09YGrcgvQ9pud2qpmME81sFpk8Twsz9uqydo1NjZXO/oybOuySMl+F915iZ+U+Opr5AADOZTbrD+ecMZ0h76WZ+yhrf8p6yZm8J+vhq2xb4P/+97/XihUrFBcXpzvuuEMxMTGlLrztNnnyZEVHRzs+D+ArzDReOZ1fYGpfFb3y7+VE0xmr/K1BTsO6IZp9e4LeXrlD753T+GdoQmvddk7jH1+oFQDgW+zM+kLD0GGT2XA894ypcf6U9U7ul6yHL7Jtgb948WK1bNlSq1evrpaFPVATmW28ElorWKcLSn5kzC8FB3mUb2KR70TTGav8qUGOVPr36uCJ03pmyWZ9sflny01ynKwVAOA77M76II/HdDbUDwvRwROnKxznT1nv5H7Jevgi25rsFRQUKCkpicU94BArjVdSujQ3tc8rL2xiatyA+Bbqb3KfZhvZWOXE/E40yJEcagjoUK0AAN/hRNY3Cg8xnSFDurcyNc6fsl7yn4Z4ZD3sYNsCv1u3bvrpp5/s2p1p06ZN01133aV77rlHf/vb37Rrl7l7hwB/Y6XxyoT+cQoOKv+er+Agj14bGl9h8xdvgxiPzN1D5tTHszoxv7dBTnmsNsiRrH2vzHKqVgCA77A76z0eqeV5dUxnSMfmDUzV6U9ZLzmToWQ9fJVtC/wnnnhCq1ev1tKlS+3apSnPPPOMXnvtNU2ZMkX333+/2rVrp4kTJ5p+flxcXKlfmZmZDlYNWGel8UpwcLAW392jzOAPDvJo8d09VKtWLU0ddkmZQXFug5iFGWbnNzfOKifm9zbIMXP8Vlj5XpnlVK0AAN9hd9ZfeH59eTwe0xmyaL25ZnT+lPWSMxlK1sNX2XYPfu/evTV79mwNHz5c1157rfr06aMWLVqU2VziiiuuqNJ8V1xxhW6//XZddtllat68uXbv3q0PP/xQzzzzjMaNG6cGDRro/vvvr9IcQHUx00zHauOVji0a6j8Tr9b4xZs0P32vTucXKLRWsAbGR2n8DR0VHBws6WyDmDljEpWala0F6XuK5h8Q3+LsO/eVnN/scZnlZIMcM8fvK7Xa2cwHAFC9KspFu7P+/S8nFD3HTN4FatZL9uc9WQ9fZdsCX5JOnDihkJAQzZw5UzNnzix3bEGBuc6fZZkwYUKxv19wwQV6/PHHdckll+jqq6/W+PHjdccdd6hOnTrl7mfjxo2lbo+Li6tSfYBZZpvpVKbxSnBwsCYOuFgTB1xc7nM8Ho+SYiPL/cgVq/ObPS6znGw8Y+b4rfCXZj4AgOpjJhedzHqp4rwL5KyX7M17sh6+yraXgGbMmKHRo0dr37596tKliwYNGqThw4eX+eWUvn376pJLLtGRI0eUlpbm2DyAHSw103G58YqV+Z1oPNO/s9nGO+bGOclfmvkAAKqH2Z/hZpvMkfVkPVAW297Bf+6551S7dm198sknSk5Otmu3ldK+fXutXbtW+/aZu48IcIuVBi3exivljXey8YqV+S0dl8lX0TfuO2pq3Ka9x3RZO3OfDuAUJ75XTpxTAED1MPsz/P6rPGS9CWQ9WY+y2fYO/o4dO9SrVy/XF/eSdPjwYUlSeHi4y5UA5bPSoMXtxitW5nei8cycNT+aGvfemt2m9+kUf2nmAwCoHuZ/hu8l600g68l6lM22d/BbtGihunXr2rW7Sjtw4IBWrFghSeratavL1QDls9qgxYmGcFY42ZCvIsdzz9g6zmn+0swHAOA8Kz/DyXr7xjmNrIcvsm2BP3z4cL300kvKzs5WRISzn824atUq7d+/XzfccENRJ3Dp7FUEw4YNU05Ojvr376+WLVs6WgdQVZVp0GJ3QzjJegdcQ0a5f7d6XGbmrx8WooMnTle4z/phvtN4xuz3yszxO914CPBF69at0+eff641a9ZozZo12rPn7LtW3H8Kf2P1Z7gTWS9Zy3u7s97M/IGc9VLFx0/Www62LfAff/xxZWRk6Morr9TLL7+sXr16OfYK43/+8x+NGjVKzZo1U9euXXXeeedp586dWrdunXJzcxUXF6c333zTkbkBO6XER2nu2oovM3OqmY5krQOu2bFWjsvsPocktNKkJVsq3OewhNYVjvElTpxTIFBMnDhRCxcudLsMoMp84We42bxxKpfM7DdQs14yd/y+8O8E/s+2e/Dbt2+v9PR0/fvf/9ZVV12lOnXqKDo6WjExMSW+YmNjqzRXQkKCfve73ykqKkrffvutPvjgA/373/9Wly5d9MILL+jbb79V06ZNbToywDneBi3lcbKZjpVurVbGmj2uhLaNTO9zdI+2ahBW/muSDcJqaVSP6HLH+BInzqlT/1YANyQlJenJJ5/UokWLtG/fPtWuzbtW8E9u/ww3mzeFhYWO5JLZ+W+7LDrgsl4yf/4T21b8b4CsR0VsbbK3Y8eOooXA6dOntWvXrqLt535t3769SnN16NBBf//737Vu3Trt379fZ86c0ZEjR7R69Wo9+OCDqlOnjk1HBTjL7cZ5Vrq1Whlr9rhSt2eb3mdQUJCW3NuzzOBvEFZLS+7tqaAg236sOc6Jc+r0vZlAdXr00Uc1YcIE3XDDDWrWrJnb5QCV5vbPcLN58/bKHY7kktn51+w4EnBZL5k//2nbD5P1qDLbLtEvLCy0a1dAjeJmMx0nurUuSN+jpNhIU8dlZf6k2Ei1iqyrjHF9NH3lDs1K26XjuWdUPyxEwxJaa1SPaL8LfKvH73bjJQBA5flD3r+XtsvUOKu5ZCXv/nJTp4DKesl63pP1qArbFvgAKs+pZjoVcaJb67ljDcPQv/ce0Zod2UUB3e78cCW0bSSPx1Op+T0ejzpENVD36Iii0OsQ1cAnA6+iZjqVPX67mvkAAKqXr+e92e70VrLeyvzecYGU9VLljp+sR2WxwAdqMCe6tXrH7j50Ute9skLHcvOLHjt44rQmLdmiV/61TUvu7Wl5fisNAd1mplYnu+X607kCqkNcXFyp2zMzM6vcGwjwdWZzxGwXeytZ3yqyrqW886f8MlurU3nvT+cK1afS17hs3LjRlgLs2g8A61Lio0yNGxDfwtLYwsLCEoF/rmO5+brulRXq37m56X1aaUjnNrO19u9i/vidmN8XzhUAwHlmM3yoye70VrK+sLDQwu8QUX6TX1ay1srvUE7Mj5ql0gv8Tp06afDgwdqwYUOlnp+enq6bb75ZnTt3rmwJAKrISgdcK2OnrdxeZuB7HcvN1+Z9x0zv00pDOreZrTXI43GkW64/nSugumzcuLHUL969R01gNsNv6xFte9ZPX7nD9PyFhuE3+WUla534FAWyHmWp9AL/qaee0pIlSxQfH68uXbroL3/5i1JTU5WXV/o9Jrm5uVq9erWeffZZXXzxxbrkkku0dOlSPfXUU5UuHkDVWOmAa2Xs7LSKP8NVkt5bs9v0Pp1oCOgU87XudaRbrj+dKwCA88xmeFBQkO1ZPyttl+n5F2XsM7VPX8gvK1nrxKcokPUoS6XvwR83bpzuvPNOTZo0Se+++64ee+wxeTwe1apVS61atVKjRo1Uv359HT9+XNnZ2dq9e7cKCgpkGIYaNmyo+++/X4899piaNGli5/EAsMhKV1+zY8026Tmee8b0Pp1oCOhlpUGN3c10nOiq7OS5AgD4J7N540TWm92v0/llNu+daJxnd96T9ShLlZrsNW3aVC+//LL+/Oc/64MPPtDHH3+sb775RllZWSXGNmvWTD179tR1112nm2++WWFhYVWZGoCNrHRrNTPWbJOe+mEhpvfpCw1qnGqmY3dXZSeb9wEA/JfZvHEi683s1xeazzrZOM/OvCfrURZbPkiyTp06GjFihP7xj39o3759+vnnn/X999/rm2++0ffff6+ff/5Ze/fu1dy5czV8+HAW90CAG9y9palxQ7u3Mr1PtxvUuN1Mxwq35wcABD5/yXrJfN4XFhaS9fB7tizwf6lJkyaKi4vTZZddpri4OC7DB2qYuOYNTY3rGNXA9D7dblDjdjMdK9yeHwAQ+Pwl6yXzef/2yh1kPfyeIwt8ADXbovXmmuQsNNlMR7LWENAsKw1q3G6mY4Xb8wN2WbJkiRITE4u+Tp8+eznwuduWLFnicpVAzeQvWS+Zz/v30naZGkfWw5dV6R58oCJWmpf5CyeOye4mb04qLCzUtJXbNTttt47nnlH9sBANSWil0T3aKijo7GuGVhu/mD0mqw1qKtqvEw1qnGqmY5Xb8wN2OHDggNLS0kpsP3fbgQMHqrMklMLtXHKCU8dkZ5M3JzmR9ZK547I7663UarZxIFkPX+YxDMNwuwhfFBcXJ+ns5+aicspqUiL971VFb/Myf+HEMVnZp9vndPehk7rulRWlfu5tg7BaWnJvT7WKrKtHP9yguWsr/vicWy5ppcev7eDIMZk5V3/6ZLPpOiWZHvuXmzpZrhcwg2yyH+e0atzOJSc4dUxm92vn/O8lJkqShqammq7Tiaz/y02dXP0dymytbRuHa/vBnArHkfVwUlVziUv04QgrDcn8hRPH5FSTNycUFhaWGfiSdCw3X9e9skKFhYUWGr9EOXJMZs9VShfzDWpoZgMAxbmdS05w6picaPLmBGeyvoXrv0OZrXVoQmtT48h6+DIW+HCElYZk/sKJY3KqyZsTpq3cXmbgex3Lzdf0lTtMN34pNAxHjsnsufLWUVGdiTERNLMBgF9wO5ec4NQxOdHkzQlOZH1iTITrv0OZrfW2HtFkPfweC3w4wkpDMn/hxDE51eTNCbPTKr60TZJmpe0y3fhlkcnGO1aPyey5Wpix13SDGprZAEBxbueSE5w6JieavDnBiaz3eDyu/w5lttagoCCyHn6PJntwhBPNy9zmZEM2N/ZpldnGM95xZhq/OPXvxMp+rTSooZkNAPwPWW//fq02ebObE1lvpV4nf4cyWytZD3/HAh+OiAyvbes4X1CZY6qos6sT58mpc1o/LEQHT5w2Nc7L4/EoKTZSSbGRpY514pxWZr8V1flLhoxy/26lVgDwV2T9WXbmktms9aesl6yfVyey3mytXmQ9/JXtC/yvvvpKX3/9tfbt26e8vNJfWfN4PJo2bZrdU8OHpMRHmepW6k9NSqweU1mdXeeu3V10iZeVfRoyXD2nQxJaadKSLRWOG2ayQY3kzDltWDfEsX9/Zue3OhYA/BFZb38uDU1orWeWbDY9v92cyHrJ2nkl64Gqse0e/KNHjyo5OVlXXXWVxo8fr6lTp2rGjBllfiGwBWJDMivHZLaza2Lbis+BrzR5G92jrRqElf+aYIOwWhrVI9r0Pp04p4ZhOHKu/OkTDwCgOridS05wO5fcbvLmRNZL5s9rQttGZD1QRbYt8B999FF9/fXXio2N1eTJk7Vo0SItW7as1K8vv/zSrmnhowKxIZmVYzLb2TVt+2G/afIWFBSkJff2LDP4vZ+NGxRk/seKE+c0NSvbkXPlT594AADVwe1ccoLbueR2kzcnsl4yf15Tt1ecjWQ9UD7bLtFfuHChzj//fKWmpioiwn9eqYVzArFJidljstLZNSk20m+avLWKrKuMcX00feUOzUrbpeO5Z1Q/LETDElprVI9oy4EvOXdO7T5XTnQA9tYKAP7K7Vxygtu55PY5dSLrJXPHRdYDVWfbAv/o0aO69tprWdyjGKvNy/yBmWOy2tnVynly+5wGBQVpdM8Yje4ZY+t+K2pmU9kOvGab5Fjdb1XHWR0LAL7K7VxyghNZb3a/VsY5xamsl8rPZbIeqDrbFvjt27dXTk6OXbsD/FogdhZ2itkGNVbPqd2Nb/zpEw8AAM4j660xk8tkPVB1tt2Df++99+qrr77Stm3b7Nol4LdS4qNMjfOnzsJOsNKgxso5daLxjZX5+f4DQODjZ715ZnO5f5fmpvZH1gNls22Bf/vtt+u+++5Tr169NH36dP3444927RrwO4HYWdgJVhrUWDmnTjS+sTI/338ACHz8rDfPbC4HeTxkPVBFlV7gBwcHl/iaPHmy9u3bp9tvv11t2rQpdUxwcLBq1bLtzgDAJwViZ2EnWGlmY+WcOtEkx8r8fP8BIPDxs94887m8l6wHqqjSK+1WrVrxDxYoh5XOroZhaHXWIS1M31s0LiU+SkkxkSX+nxUWFmrayu2anba7qLPtkIRWGt2jbaU721qZ305Wm9mYPadONMmxMr/VsQAA/2T1Z73ZvHUi663MbzcruUzWA1VT6QX+jh07bCwDCExmuuBaaRCz+9BJXffKCh3LzS8ad/DEaU1askWv/GubltzbU60i61qq0e4GNVZUppmNmXPqZOMjf/rEAwCA88z+rDebt05kvZX5nWA1l8l6oPJsuwcfgHVWGsQUFhaWCPxzHcvN13WvrFBhYaEj8zvBqQY1NL4BAPgSs3lbUFBge9Zbmd+f8p6sB0pn2wK/d+/eeu655yocN3nyZPXu3duuaQG/ZqVBzLSV28sMfK9jufmavnKHI/M7wakGNTS+AQD4ErN5+9SiTbZnvZX5/SnvyXqgdLYt8L/66itt2bKlwnFbt27V8uXL7ZoW8GtWGsTMTtttauystF2OzO8EpxrU0PgGAOBLTOdthrlxVrLe0vx+lPdkPVC6am9nn5ubSxd94L+sNIg5nnvG1Nhzx1XUTMepBjVWONWghsY3AABfYTZHT+cXmBpnNuutzu9veU/WAyVV60r72LFjWrVqlZo3b16d0wI+y0qDmPphITp44nSFY+uHnW2QY6aZjpMNaqxwqkENjW8AAL7AbI6G1grW6YLyL9GXzGf9tYWGgoM8AZ33ZD1QXJUu0Y+JiSn6kqQPP/yw2LZzv1q3bq0mTZpox44dSklJsaV4wN+ldDHbICZKg7u3NDV2aPdWppvp9O9i7sU2GtQAAFB5ZvM+pXMzU+PMZn3WwZyz+6UhHVBjVOkd/HM/Ks/j8ejEiRM6ceJEqWNDQkIUFRWl/v3769lnn63KtCUcOnRIHTp00IEDBxQbG6tt27bZun/AKYbMdas1DCmueUNTYztGNTDdTOf+qzxKjIkodywNagAAqBqzed8mMtzUOLNZfyI3X8dz84sa0pH3QOCr0gL/3I/oCAoK0siRI/X2229XuSirHnroIR08eLDa5wWqalHGPlPjFppsunN2rLl9StKC9L2aOuySMt8BoEENAABVZzbv53z7o6lxVrL+8MnTRQ3pyHsg8Nl2D/706dPVrl07u3Zn2r/+9S+98847uuOOO/TGG29U+/xAVTjR9MbqWBrUAADgLLPZbLahrpWszy84e/UAeQ/UDLYt8EeMGGHXrkw7deqUxo4dq44dO+rhhx9mgQ+fU1hYqGkrt2t22m4dzz2j+mEhGpLQSqN7tFVQUJAjTW8qM9ZKg5qKOvNXdiwAAP6ooqyXzGez2Ya6VrK+VvD/8tZs3pP1gP+ybYH/9ddfmxoXGhqqyMhItWvXrsr/6Z9++mllZWVp+fLlCgkJqdK+ALvtPnRS172yQsdy/9cN9+CJ05q0ZIte+dc2Lbm3p1LiozR3bcWfbz8gvoUMGY6MtcJMZ/6Gdc138feOBQDAH5nJ+laRdU3n/dCE1npmyeYKx1nJ+kZ1Qysccy6yHvBvti3wk5OTLS3Y69Wrp8GDB2vSpEmKjLT+sRYbNmzQCy+8oFGjRqlnz57FGv4BbissLCwR+Oc6lpuv615ZofQnf22p6Y3ZsYZhqEFYrTLnl6QGYbWU0LaRySOS6c78c8YkSpLpsby6DwDwR2azPmNcH9NN7m7rEa0vNv9sW9YHB3lUL8z8r/tkPeD/qvQxeecaPny4brjhBhmGIY/Ho/j4eKWkpCglJUVdu3Yt+o99/fXX66qrrlLt2rX1xhtvKCkpSYcPH7Y0V2FhoW6//Xadd955eu655+w6BMA201ZuLzdwpbPBP2PVTk0ddkmZXWvPbXrjbZBjZmzq9mxT86dtN/9/z2xn/tSsbEtjAQDwR2azfvrKHaYzPCgoyNasLyg0dKKCMeci6wH/Z9s7+M8//7wSEhJ09dVX629/+5vat29f7PFt27bpvvvu0/fff6/U1FSFh4drzJgxmjt3rp5//nn96U9/Mj3XK6+8om+//VbTp0+v1Lv/54qLiyt1e2ZmpmJjY6u0b9Rcs9MqvmROkmal7dLonjGmm96YbZCzMN1c1/0F6XtM3XdvdZ9mWZkfAABfYjXrzWa43Vl/+GTF9/R7kfWA/7Ntgf/YY48pLy9P8+fPV1hYWInH27Vrp48++kjt2rXTY489pmnTpun111/XZ599pkWLFple4O/atUtPPPGEevXqpZEjR9pVPgKUW41fzHbBPXfcLz8jt6zPzDXTIMfN7vxOdfwHAKA0/pT1krm8tzPrvV30zSDrAf9n2wJ/yZIlSk5OLnVx71WnTh317NlTn3zyiSSpQYMGio+P16pVq0zPc/fdd+v06dN6/fXXq1yzJG3cuLHU7WW9sw//4WbjF7NdcOuHhThSp5vd+Z3q+A8AwC/5S9Y7UavZDD23i75d+yTrAd9l2z34R48e1dGjRyscd+zYsWLjGjdubGmejz/+WHXr1tWdd96p5OTkoq/f/va3kqQ9e/YUbfvpp5+sHQQChtkmMYZh/lVtK4YktDI1bmj31o7UmdIlytS4AfHmxklSismxA+JbWBoLAEBl+EvWD0to7UitZrPeShd9sh7wf7Yt8Nu3b69ly5Zp8+ayP9pj8+bNWrZsmS644IKibfv27bO8yD9y5IiWL19e7CstLU2SlJubW7QtNze3cgcDv+d245fRPdqqQQVdaxuE1VKH5vUdqbOsy/tLjLPwO4+3A3B5vJ19rYwFAKAy/CXrR/WIdqRWs1lvBVkP+D/bFvh333238vLydMUVV+jZZ5/VDz/8oNzcXOXm5uqHH37Qn//8Z/Xq1UunT5/W3XffLUk6deqU1q1bp27dupmexzCMUr+2b98uSYqNjS3aFh0dbdfhwc840STGiqCgIC25t2eZwd8grJaW3NtTi9bvM7U/q3UuyjC334UZ5s6TJEtd/K2MBQCgMvwl64OCghyp1WzWW2myR9YD/s+2e/DvuOMObd68WS+//LKeeOIJPfHEEyXGGIahBx54QGPGjJEk7dy5U8OGDdONN95oVxmAJGeaxFjVKrKu0p/8tZ5avEkL0vfqdH6BQmsFa0B8lJ6+oaOCg4Mdq9Op/Zrt7Gt1LACU5dSpU3r22Wf1/vvva9euXYqIiFC/fv00ceJEtWjBpb81mb9kvVO1OtFkTyLrAX9n2wJfkl566SXddNNNeu2117Rq1Srt23f2lcXmzZurR48eGjt2rC6//PKi8RdddJFee+01O0sAJDnTJMaq0prpnC7I16zUXdq2/4SmDrvEsTqdPH4znX0rMxYAfik3N1e9e/dWamqqmjdvrpSUFO3YsUPTp0/Xxx9/rNTUVMXExLhdJlziL1nfsG6Iq83rrDTZ8yLrAf9l2yX6Xj169NCsWbOUlZWlU6dO6dSpU8rKytLMmTOLLe4BJ7nd+MVsM53+XZqb2p/VOt0+fgCwwzPPPKPU1FQlJSXpP//5j+bOnau0tDS98MILOnDggG677Ta3S4SL3M46K43znKjV7D6tNNkD4P9sX+C7JTo6WoZhaNu2bW6XAh/gduMXs810gjweR+p0+/gBoKpOnz6tV199VZI0ZcoU1atXr+ixBx98UJ06ddLy5cu1bt06t0qEy9zOOiuN85yo1cw+64XVUv0KGgECCCyOLPDz8/P1888/a9euXWV+AU5yu/GL+WY6ex2p0+3jB4CqWrlypY4eParY2FjFx8eXePymm26SJC1evLi6S4OPcDvrrDTOc6JWM/uMaRxuen8AAoOtL+l98cUXRZfTnTlzpsxxHo9H+fn5dk4NlOBm4xcrzXScqpPGNwD82fr16yVJXbt2LfVx7/YNGzZUW03wPf6S9U7VWtE+Z79N1gM1jW0L/I8//lgDBw5UQUGBGjVqpLZt26p+/fp27R6oFLcav1htpuNUnTS+AeCvvFf7tWzZstTHvdt37txpan9xcXGlbs/MzFTjoCC9l5hYiSrhSzqd8+es/345KTb7pM47UfFH0EXWC9V7r9Utts2JWkvb58F//1uS+PcN+JGj27erYdu2lX6+bQv8p59+WoWFhXrppZd0zz33FH0sCFAT9e/SXHPX7q5wXEoXcw1yAKCmOXHihCSpbt26pT4eHn720uPjx49XW03AuRrVDdUhEwt8mtwBqE62LfA3btyopKQk3X///XbtEvBbHpm7JI6r5AGgemzcuLHU7d539oemplZnOQgAq7Yd1J/fSqtw3JwxCUqKbVwNFZXkfeeef9+A//hTGVecmWVbk7169eqpdevWdu0O8GsLM8w32QMAlOTtmn/y5MlSH8/JyZEkbgeEa8h6AL7ItgX+r3/9a61du9au3QF+zWrjHQBAcd43DX788cdSH/dub9OmTbXVBJyLrAfgi2y7RP8vf/mLLr30Uj366KOaNGmSatXiMzdRc1ltsmcYhlZnHdLC9L1FHXBT4qOUFBNJt3sANVLnzp0lSd99912pj3u3d+rUqdTHAadZzXqJvAfgPNtW4dOnT9c111yjyZMn66OPPlJycrJatmypoKCSFwl4PB49+eSTdk0N+JyU+ChTTfYGxLfQ0ZNnNHbWWqVmZRd7bO7a3UWfi9uwbohTpQKAT+rRo4caNmyozMxMZWRkqEuXLsUe//DDDyVJN9xwgwvVAdayXhJ5D6Ba2LbAHz9+vDwejwzDUFZWlrKyyv7ADxb4CHRJMZFKjIkoEeLnSoyJUELbRhryVlqZ41KzsjV21lrNGZPIK/sAapTQ0FDdc889mjRpku6++2599tlnRZ3zX3zxRW3YsEG9evVSt27dXK4UNZXZrE+MiZBhGKUu7r3IewB2sfUdfABneTweTR12SZlh7n2lPnV7drm/GEhnQz81K5vPsgdQ4zzxxBP64osvtGrVKrVv3149e/bUzp07lZaWpiZNmujtt992u0TUYGaz3uPxaFXmQfIeQLWwbYE/YsQIu3YFBISGdUM0Z0yiUrOytSB9T9G9dgPiWygxJkIej0cLTXbWXZC+h8AHUOOEhYVp2bJlevbZZzV79mwtWLBAERERGjlypCZOnKiWLVu6XSJqODNZL4m8B1Bt6IQHnxCoTWc8Ho+SYiPLDOvKdOAN1HMFAKWpU6eOJkyYoAkTJrhdCmwQiBlWUdZL1vM+EM8TgOph+wI/Pz9fS5Ys0Zo1a3Tw4EElJCTotttukyTt3btXBw8eVMeOHemyjyI1uemM1Q68NflcAQD8W03OMCt5X5PPE4CqK9nivgq++eYbtWvXToMGDdKzzz6rt956S998803R46tXr1Z8fLwWLVpk57TwY2abzhiGUc2VVY+U+ChT4wbEt6jx5woA4L9qeoaZzfuULlE1+jwBqDrbFvibNm1Sv379tG/fPt1777364IMPSvzwueGGG1S3bl199NFHdk0LP7c665DppjOByNuBtzzeDrw1/VwBAPxXTc8ws3kvqUafJwBVZ9sCf+LEicrNzdXixYv117/+VTfddFOJMaGhoeratavS09PtmhZ+zkrTmUDk7cBbVuif24G3pp8rAID/qukZZjbvF2bU7PMEoOpsuxF+2bJl6t69u/r27VvuuBYtWmj9+vV2TQs/V5kmc4HGbAdezhUAwF+RYebynvMEoKpsW+AfOXJErVq1qnBcTk6Ozpw5Y9e08HNWm8wFKjMdeDlXAAB/RYadVVHec54AVJVtl+g3bdpU27Ztq3Dc5s2bTb0QgJrBSpO5mo5zBQDwV2SYOZwnAFVl2wK/d+/eysjI0LJly8ocM3/+fG3btk19+vSxa1r4OStN5mo6zhUAwF+RYeZwngBUlW0L/P/7v/9TaGioBgwYoNdee00//fRT0WOHDx/W22+/rdGjRys8PFwPPvigXdPCz1lpMlfTca4AAP6KDDOH8wSgqjyGjR+kuWDBAt166606efJkqY+HhYVpzpw56t+/v11TOiYuLk6StHHjRpcr8W+GYWh11iEtTN9b1EwmJT5KSTGRxcKpsLBQb6/coffSdul47hnVDwvR0ITWuq1HtIKCbHsdKiAYhlFhQz4AgYlssh/ntOrMZr1E3ptlV9a/l5goSRqamupUqQBsVtVcsnWBL0k7d+7USy+9pM8//1w7duxQYWGhWrZsqT59+uihhx5SbGysndM5hsCvuqMnz2jsrLWlflar9xXohnVDTI8DgJqObLIf57RqrGQ4eV/9WOAD/qequWT7S6Vt2rTRX//6V23cuFE5OTk6deqUfvjhB/3973/3m8U9qs4wjDJDXJJSs7I1dtZaFRYWmhpn8+tQAACgisxmvWEYlsYCACqPa6HgiNVZh8oMca/UrGy9vXKHqXEVjQEAANXLbNanZmVbGgsAqDwW+HDEwvS9psa9l7bL1LgF6XuqUg4AALCZ2axfkL7H0lgAQOXVquwTY2JiKj2px+NRZmZmpZ8P33coJ8/UuOO5Z2zdHwAAqB5ms9lKhpP3AFA1lV7g79ixw8YyEGgiw2ubGlc/LEQHT5y2bX8AAKB6mM1mKxlO3gNA1VT6Ev3CwsIqfSGwpcRHmRo3NKG1qXED4ltUpRwAAGAzs1k/IL6FpbEAgMrjHnw4IikmUokxEeWOSYyJ0G09ok2Nq2gMAACoXmazPjEmwtJYAEDlscCHIzwej6YOu6TMoPZ+3m1QUJCpcR6Px8lyAQCARWaz3uPxWBoLAKi8St+DD1SkYd0QzRmTqNSsbC1I36NDOXmKDK+tAfEtlBgTURTiZscBAADfYiXDyXsAcB4LfDjK4/EoKTZSSbGRtowDAAC+xUqGk/cA4Cwu0QcAAAAAIACwwAcAAAAAIAD47QL/xRdf1KBBg9S+fXs1bNhQtWvXVps2bTR8+HB9//33bpcHAAAAAEC18tsF/p/+9Cd9+umnioiI0FVXXaXrrrtOYWFhmjlzprp166aPP/7Y7RIBAAAAAKg2fttkb+HCherWrZvCwsKKbf/73/+uu+++W7fffrt+/PFH1arlt4eIMhiGodVZh7QwfW9RB96U+CglxUTSgRcAgABA1gNA5VR69fvuu+9WaeLhw4dX6fk9evQodftdd92lF198UZmZmdq0aZM6depUpXngW46ePKOxs9YqNSu72Pa5a3cXfYZuw7ohLlUHAACqiqwHgMqr9AJ/5MiRlXoF1TAMeTyeKi/wyxMScvaHfmhoqGNzoPoZhlFq4HulZmVr7Ky1mjMmkVf3AQDwQ2Q9AFRNpRf448aN88kfrDNnztTWrVvVvn17tW/f3u1yYKPVWYfKDHyv1KxspWZl8/m6AAD4IbIeAKqm0gv88ePH21hG5T3//PPauHGjcnJytHnzZm3cuFFRUVGaM2eOgoODK3x+XFxcqdszMzMVGxtrd7mogoXpe02NW5C+h9AHAMAPkfUAUDV+34Hun//8p/71r38V/b1NmzZ699131a1bNxerghMO5eTZOg4AAPgWsh4AqsbvF/hffPGFJOnIkSP6/vvvNWHCBPXq1UvPPPOM/vjHP1b4/I0bN5a6vax39uGeyPDato4DAAC+hawHgKqxfYG/a9cuLV68WD/88IOOHz8uwzBKjPF4PJo2bZqt85533nnq2bOnPvnkEyUlJenJJ59U3759demll9o6D9yTEh+luWt3VzhuQHyLaqgGAADYjawHgKqxdYE/YcIETZw4UYWFhUXbvAt8b0M+bxd9uxf4XiEhIbrlllu0bt06LV68mAV+AEmKiVRiTES5zXcSYyKUGBNRjVUBAAC7kPUAUDVBdu1o7ty5Gj9+vFq1aqU33nhDffr0kXT2HvnXXntNvXr1kmEYevDBB/Xll1/aNW2pGjduLEk6cOCAo/Ogenk8Hk0ddkmZoe79bFxf/HQHAABQMbIeAKrGtnfw//73vys0NFTLli1TmzZt9M0330hS0UJ/7Nixeumll/SHP/xBAwYMsGvaUi1fvlyS6IIfgBrWDdGcMYlKzcrWgvQ9OpSTp8jw2hoQ30KJMREEPgAEgJycHM2bN09r1qzRmjVrlJGRodOnT+upp57ymU/xgXPIegCoPNsW+Bs2bNBll12mNm3aSCp5Sb4k/f73v9e0adP0zDPPaOnSpZWea+XKlTp+/Lj69u2roKD/XYRw5swZvf7665o5c6bq1KmjW265pQpHBF/l8XiUFBvJx+MAQID64YcfNHz4cLfLgIvIegCoHNsW+Hl5eWrWrFnR38PCwiSd7W7fqFGjou2dO3eu0uJeOhv8o0aNUuPGjdWtWzdFRkbq4MGD+v7777Vv3z6FhYVpxowZatWqVZXmAQAA1a9+/foaPXq0Lr30Ul166aVasmSJxo0b53ZZAAD4PNsW+M2bN9f+/fuL/t6ixdnuphs3btTll19etP3HH39UQUFBlebq1auXHn/8cS1fvlwbNmzQwYMHFRoaqujoaN10002677771K5duyrNAQAA3BEbG6u33nqr6O+fffaZi9UAAOA/bFvgX3zxxVq3bl3R35OTk2UYhp566iktWrRI4eHh+uCDD7RixQolJSVVaa62bdtq0qRJVS0ZAAAAAICAYVsX/RtuuEF79uwp6pDfo0cPXXnllVq2bJkaNWqkxo0ba/DgwfJ4PHryySftmhYAAAAAAMjGBf6wYcO0efNmdenSpWjb/PnzdccddygiIkLHjx9Xx44dNXPmTPXr18+uaQEAAAAAgGy8RL927dq68MILi21r0KCBXn/9db3++ut2TQMAAGBZXFxcqdszMzP5WF0AQMCwbYH/9ddfq1mzZrrgggvKHffDDz9o3759uuKKK+yaGgAA+JCBAwdq8+bNlp7z7rvvqnv37g5VBABAzWDbAj85OVmjRo3StGnTyh333HPP6e23365yJ30AAOCbtm/frq1bt1p6zsmTJx2q5qyNGzeWur2sd/YBAPBHti3wJckwDFvGAAAA/5WRkeF2CQAA1Ei2Ndkza+/evapXr151TwsAAAAAQECr0jv47777brG/b9u2rcQ2r/z8fG3dulVffPGFEhMTqzItAAAAAAD4hSot8EeOHCmPxyNJ8ng8WrlypVauXFnmeMMwFBYWpnHjxlVlWgAAAAAA8AtVWuCPGzdOHo9HhmFowoQJ6tKli1JSUkodGxoaqqioKPXt21fNmzevyrQAACDADRw4UPv27ZN09vY+SXrrrbe0dOlSSVLz5s01f/581+oDAMAXVWmBP378+KI/z5gxQ7/+9a/11FNPVbUmAABQw6Wnp2vnzp3Ftu3Zs0d79uyRJLVp08aNsgAA8Gm2ddHfsWOHXbsCAAA1HL9XAABgna0fk+e1fv16rVmzRgcPHlRcXJz69+8vScrLy1NeXp4aNGjgxLQAAAAAANRYtn5M3tatW3XZZZepa9euuvPOO/XEE09owYIFRY/Pnj1bjRo1Krp/DgAAAAAA2MO2Bf7u3bt1xRVXKDU1VTfccIOee+45GYZRbMzNN9+s0NBQffTRR3ZNCwAAAAAAZOMCf8KECTp48KDeeustLViwQA899FCJMeHh4erSpYvS0tLsmhYAAAAAAMjGBf7SpUvVqVMn3XbbbeWOi46OLuqACwAAAAAA7GHbAn///v268MILKxx35swZnTx50q5pAQAAAACAbFzgR0ZGateuXRWO+89//qPmzZvbNS0AAAAAAJCNC/wePXro22+/VUZGRpljli9frn//+99KTk62a1oAAAAAACAbF/gPP/ywDMNQSkqKPv30UxUUFBR7/Msvv9Stt96qWrVq6YEHHrBrWgAAAAAAIBsX+AkJCfrb3/6mvXv36vrrr9d5550nj8ejjz76SI0aNVKfPn20d+9evfrqq+rUqZNd0wIAAAAAANm4wJeku+66SytWrNANN9wgj8cjwzB0/Phx5eXl6eqrr9by5ct1xx132DklAAAAAACQVMvuHSYmJmrBggUyDEMHDx5UYWGhGjdurODgYLunAgAAAAAA/1XlBf4nn3yiBQsWaPfu3apdu7Y6deqkUaNGqW3btmrSpIkdNQIAAAAAgApUaYE/dOhQvf/++5IkwzAkSYsXL9bkyZP1/vvvq3///lWvEAAAAAAAVKjSC/xp06Zpzpw5qlWrlm699VbFx8fr+PHj+vjjj7V69WoNHz5cO3fuVMOGDe2sFwAAAAAAlKLSC/x33nlHQUFB+vTTT3XVVVcVbX/sscc0atQovfvuu5o3b55GjRplS6EAAAAAAKBsle6i//333ysxMbHY4t7r8ccfl2EY+v7776tUHAAAAAAAMKfSC/xjx44pNja21Me8248dO1bZ3QMAAAAAAAsqvcA3DKPMj74LCjq728LCwsruHgAAAAAAWFDpBT4AAAAAAPAdVVrgv/POOwoODi71y+PxlPl4rVpV+nQ+AAAAAADwC1VaaRuGUa3PAwAAAAAApav0Ap/76wEAAAAA8B1+eQ/+yZMntWDBAo0ePVoXXnihwsLCFB4ers6dO2vChAk6ceKE2yUCAAAAAFCt/HKBP3v2bA0cOFBvv/22goOD1b9/f/Xs2VPbt2/XU089pUsvvVT79+93u0wAAAAAAKqNXy7wQ0JCdMcdd2jTpk3atGmTPvjgAy1dulRbt25VfHy8tmzZogceeMDtMgEAAAAAqDZ+ucAfMWKEpk6dqg4dOhTb3rx5c02ZMkWSNG/ePJ0+fdqN8gAAAAAAqHZ+ucAvT+fOnSVJeXl5OnTokMvVAAAAAABQPQJugZ+VlSXp7GX8ERERLlcDAAAAAED1qPTH5Pmql19+WZLUr18/1a5du8LxcXFxpW7PzMxUbGysrbUBAAAAAOCUgHoH/5NPPtG0adMUEhKiiRMnul0OAAAAAADVJmDewd+yZYuGDRsmwzD0/PPPF92LX5GNGzeWur2sd/YhGYah1VmHtDB9rw7l5CkyvLZS4qOUFBMpj8cT8PMDABDo3M5at+cHAH8VEAv8PXv2qF+/fjp8+LAefPBB3X///W6XFLCOnjyjsbPWKjUru9j2uWt3KzEmQlOHXaKGdUMCdn4AAAKd21nr9vwA4M/8/hL97Oxs9e3bVzt37tSoUaM0efJkt0sKWIZhlBq4XqlZ2Ro7a60MwwjI+QEACHRuZ63b8wOAv/PrBf6JEyd0zTXXaNOmTRo0aJDefPNNLtty0OqsQ2UGrldqVnaFY/x1fgAAAp3bWev2/ADg7/x2gZ+Xl6eUlBStWbNGV199tebMmaPg4GC3ywpoC9P3mhq3IH1PQM4PAECgcztr3Z4fAPydXy7wCwoKNHjwYH355Zfq2bOn5s2bp9DQULfLCniHcvJsHedv8wMAqseWLVv0l7/8RVdeeaUaN26skJAQNWvWTIMGDdKKFSvcLi+guZ21bs8PAP7OL5vsvfrqq5o/f74kqXHjxrrrrrtKHTd58mQ1bty4OksLaJHhtW0d52/zAwCqx69//Wvt2bNH9erVU2JioiIiIrRp0ybNnz9fCxYs0IsvvqgHHnjA7TIDkttZ6/b8AODv/HKBf/jw4aI/exf6pRk/fjwLfBulxEdp7trdFY4bEN8iIOcHAFSPiy66SM8++6x+85vfKCwsrGj71KlTdeedd+rhhx9W37591bFjRxerDExuZ63b8wOAv/PLS/THjx8vwzAq/IqOjna71ICSFBOpxJiIcsckxkRUOMZf5wcAVI8vvvhCt956a7HFvSSNHTtWffv2VUFBgf7xj3+4VF1gcztr3Z4fAPydXy7w4Q6Px6Opwy4pM1S9n03r1CcZuD0/AMB9nTt3liTt3WuuGRuscTtr3Z4fAPydX16iD/c0rBuiOWMSlZqVrQXpe3QoJ0+R4bU1IL6FEmMiHA/chnVDNPv2BL29cofeS9ul47lnVD8sREMTWuu2HtEKCnL+NavCwkJNW7lds9N2F80/JKGVRvdoWy3zG4ah1VmHtDB9b9H5T4mPUlJMJL/wAAh4WVlZkqRmzZq5XEngIuvJegD+y2MYhuF2Eb4oLi5OkrRx40aXK8G5jp48o7Gz1pb6+bfeV/Ub1g1xbP7dh07quldW6FhufonHGoTV0pJ7e6pVZF3H5nf7+AG4q6ZnU2ZmpuLi4pSXl6e1a9eqW7dupp/rPXel7TM2NrbGnlNf5HbWBVLWv5eYKEkamppqa40AnFPVrOcSffgNwzDKDDxJSs3K1thZa+XUa1aFhYVlBr4kHcvN13WvrFBhYaEj87t9/ADgpvz8fI0cOVJ5eXm65ZZbLC3u4T/czjqyHoC/4xJ9+I3VWYfKDDyv1KxspWZlKyk20vb5p63cXmbgex3Lzdf0lTs0umeM7fO7ffwAYNbAgQO1efNmS89599131b179zIfv++++/TNN98oJiZGf//73y3XVNY7IWW9sw93uJ11ZD0Af8cCH35jYbq5hkoL0vc4Enqz0yr+2B5JmpW2y5HQd/v4AcCs7du3a+vWrZaec/LkyTIfmzRpkl577TWdf/75+uc//6mICDqoByq3s46sB+DvWOBDkrVmLk40nikoKNC4RRu1MGOfTucXKLRWsFK6NNeE/nEKDg6WJB3KyTO1L7PjrDqee8bWcVa5ffwAYFZGRoZt+3r99df1xBNPqGHDhlq6dKnatWtn275rGrezXqo4793OOrIegL9jgY8ym7nMXbu7RDOX0hrPHDxxWpOWbNEr/9pWqcYzm/Yc1Q1TVqqg8H/3k50uyNd7abv1/rc/avHdPdSxRUNFhtc2tT+z46yqHxaigydOmxrnBLePHwCq2/vvv6+7775bdevW1ZIlS9SlSxe3S/Jbbme9ZC7v3c46sh6Av6PJXg1npZmLE41nCgoKSoR9sccLjbOPFxSof5fmpvaZ0iXK9PxWDEloZWrcsITWjsyfEm/uuAbEt3BkfgCoTp988omGDx+uWrVqaf78+erRo4fbJfktt7NeMp/3113c1NT+yHqyHkDpWODXcFaauVhpPGPWuEUbywx7r4JCQ+MXb5JH5j731amPhx3do60ahJV/0UuDsFoa1SPakfmTYiKVGFP+faeJMREVjgEAX7dy5UrddNNNMgxDc+fOVd++fd0uya+5nfWS+byfmWruHniynqwHUDoW+DWclWYuVhrPmJ4/Y5+pcfPT92phhtlazY2zKigoSEvu7Vlm8Hs/G7cq9yaWx+PxaOqwS8oMde8llr+8jxIA/M3111+vU6dOqVWrVlqwYIFGjhxZ4uutt95yu0y/4XbWS+bzftnWA6bGkfVkPYDScQ9+gKuooY6VZi6VaTxTUZOe0/kFpvZ5Or+gUo1nzDYJMtt4qFVkXX33xFW68710fbX1gAoKDQUHeZR8YRO9PjRetWo5+1+qYd0Qzb49QW+v3KH30nYVHdPQhNa6rUe0Y79wAEB1OnLkiKSz3fi3b99e5rjbb7+9mirybb6e9ZJM531F7/KfW6uV+b3M5D1ZD8CfscAPYGYa6lhp5mK18YyZJj2htYJ1uqD8SwElKbRWsOXGM2abBFW18VB+oaEvNu9X12f+VenGQ2aVVuvBE6f1zJLN+mLzz8VqBQB/ZRjmFnnwj6xvFVnXdN4HB3mUb2KRbzXrJfONBsl6AP6MlwADlNmGOmab1AyIb2Gp8YzZJj39O5lrnDcwPko3dG5mamz/zs1Mz19QUOB64yGzrDRJAgAEPn/J+sLCQqWYbJTb6wJzn+1uJesLCwtNn6uCggKyHoBfY4EfoMw21JFkupmLlcYzZpv0SOYC6uqO5+uT783dv/fpv38yPf/4RZtcbzxklpUmSQCAwOcvWT995Q71izu/3HFehQXmFs5Wsn76yh2mz9VTizeS9QD8Ggv8AGW2oc7CjL2mm7lYaTxjtknPog0/mRq3eMNPWrTe3NgFGftMzz/fdOM+5xoPmWWlSRIAIPD5S9bPStuljzf8bGrsiszDpsZZyfpZabssZKi5NxPIegC+invwA5SVhjoN64ZozphEpWZla0H6nqKmMwPiWygxJqJEk7mMcX00feUOzTqn8cuwhNYadU7jF7NNesw23TmUk2epIZ8T81em8ZCdKtNkEAAQuPwl64/nnjFdq9kme1ay3sr8Zn8vIOsB+CoW+D7CbBd3s6w2pPN4PEqKjVRSrLl73wp/cWn9L/9utkmP2aY7keG1LTXkc2J+q42HJHu/r1a/pwAA31JTs75+WIjpWs022bOS9VbmN/t7AVkPwFexwPcBVrq4m5USH6W5ayu+dG1AfAtL+zXbrXZIQitNWrLFxPxRmpVa8WVuA+JbqFaw9J6Jy/EGxkepdWRdU/MP7BKlmSYusxsQ30Ltzg83tc9hCa0l2f99dep7CgBwXk3O+mEJrdUhqoGpWpMvaKIvtuyvcJyVrLcy/4D45pqVWvE4sh6Ar+IefJc51S01IbqRgoPKf9U4OMij7tHnmd6nlW61o5LamNrnk9dcYLrO8dd3MLXPcdddpNsuiza136du6OBI4yEnvq9JMZGmawUA+I6anvUjklqbrvXV315sap9Wsn7kZW1MZ+jTN8SR9QD8Ggt8lznVLfXtVTsqvI+toNDQjFU7Te/TSrfa6avN7ffu9zeYrnOGiXf6Jend1N1K23HY1H6/3XnUkcZDTnxfPR6P6VoBAL6jpmf9O6t3ma71nrnfm9qnlaxfs+OI6QwNDg4m6wH4NS7Rd5mVbqlm75mTZKmz7OieMbbv06yvth6wfZ+z0nape7S5V7a959XuxkNOfV+tNEkCAPiGmp71Tv1eYDXrzWYoWQ/An7HAd1lluqWaaeZiteO7E/s0w2y3XCv7tNIt1zvOSuMhj8ejDlEN1D06ouhcdYhqUCxwneyCa7VJEgDAXTU96536vaAy59VshpL1APwVC3yXWe2WaraZi5XOsk7sU5KpsWa75VrZp5VuuVa70Jo9V3TBBQB41fSsd+r3ArIeAEriHnyX3dC5malx/Ts3t9TMZUhCK1P7HZbQ2pF9mh3b64LGpsYN7W5+n8MSWislPsrUWCtdaK2c//5dmpvaZ0oXc3UCAPxXTc/6YQmtNbi7ubHJJn8vIOsBoHQs8F22ed9xU+O2/HTcUjMXsx3fL2xWz/Z9juoRbXrsyMvMdeDtGFXf0vxOdKG1cv49Mnd/HLfRAUDgq+lZP6pHtOKaNyh3nNfIy9qQ9QBQBSzwXTZ7jflmNlaauZjt+L54/U+27zMoKMj02I83/Gxq/oUZ+yzN70QXWivnf2GG2bHmxgEA/FdNz/qgoCAtWr/PVA2LN/xM1gNAFXAPvsusNLOx2szFTBdYJ/bp5fb8dnehdaKZTmUa71hhpqESAMBZNT3rz923mRrIemvKy3oANQ8LfJdZaWZTmWYuQUFBGt0zpsyPx3Fin+dye347u9A60UzHycY7ZpsEAQCcVdOzvjI1kPXmVJT11xYaCg7iBX2gJuESfZe53TjOiX1aYbbxzACTdTrJyrly+7xaaRIEAHBWTc96yX/yPtCyPutgjiNzA/BdLPBd5nbjOCf2aYUhcwtMX1iHWjlXbp9XK02CAADOqulZL/lP3gda1p/Izdfx3HxH5gfgm1jgu8ztxnFO7NOKRRnmmu6YbWTjJCvnyu3zaqVJEADAWTU96yX/yftAzPrDJyu+PQRA4OAefIeZaXLmZjMZp/ZplhPNbJxk5VxxXgGgZiDrK+ZPuRRoWZ9f4AOXQQKoNizwHWSlyZlbzWSc3KcZTjSzcZqVc8V5BYDARtab42+5FEhZXyuYJntATcIl+g6hyZk5bjeoCVScVwBwHllvHrlkP7PntFHdUIcrAeBL/HaBv27dOv35z3/WoEGD1LJly6J7oXwFTc7McbtBTaDivAKA88h688gl+5k5p/XCaql+BQ0eAQQWv13gT5w4UY899pjmz5+vPXt8r1EYTc7McbtBTaDivAKA88h688gl+5k5pzGNw6u5KgBu89uX9JKSktSpUyddeumluvTSSxUdHa28PPcbs3j5UzMZt7nd+CdQcV4BwFlkvTXkkv0qOqez3+acAjWN3y7wH330UbdLKJevNJMx09nXF7jZ+CeQcV4BwDlkvXXkkv04pwDO5bcLfF+XEh+luWt3VzjOyWYyVjr7AgAAa8h6AICv8dt78H2d281k6OwLAICzyHoAgK9hge8Qt5vJ0NkXAABnkfUAAF9T4y/Rj4uLK3V7ZmamYmNjq7RvN5vJWOnsyz1bAABUDlkPAPAlNX6B7zS3Gp/Q2RcAgOpB1gMAfEWNX+Bv3Lix1O1lvbPvL3ylsy8AAHAGWQ8A+CXuwQ9QKfFRpsY52dkXAAA4h6wHAPwSC/wA5XZnXwAA4CyyHgDwSyzwA5TbnX0BAICzyHoAwC/V+HvwA5mbnX0BAKisDRs26I033tDatWu1a9cuHTp0SGFhYerYsaOGDBmiO++8UyEhIW6X6RPIegDAufx2gb9kyRJNnDix6O+nT5+WJCUmJhZte/LJJ3XddddVe22+xK3OvgAAVNbXX3+tKVOmqE2bNurYsaOaNGmiAwcOaOXKlUpNTdVHH32kzz77TKGhoW6X6hPIegCAl98u8A8cOKC0tLQS28/dduDAgeosCQAA2ODaa6/Vtddeq5iYmGLbf/75Z/3617/W8uXL9cYbb+iee+5xqUIAAHyT396DP3LkSBmGUe7XyJEj3S4TAABYFBMTU2JxL0nnn3++Hn30UUnSl19+Wd1lAQDg8/x2gQ8AAGoe7733XJ4PAEBJLPABAIBfOHz4sF544QVJqvE9dgAAKI3f3oMPAAAC2w8//KBJkyapsLBQP//8s1atWqUTJ07ozjvv1NChQy3tKy4urtTtmZmZio2NtaNcAABcxwIfAAD4pJ9//lnvvPNOsW333XefJk6cqKAgLkIEAOCXWOADAABbDRw4UJs3b7b0nHfffVfdu3cvtu3yyy+XYRgqKCjQrl27NH/+fD399NP69NNP9dlnnyk6Otr0/jdu3Fjq9rLe2QcAwB95DMMw3C7CF9WvX19nzpzhsj0AgM/IzMxUSEiIjh8/7nYp5erSpYvWr19v6TnLli1TcnJyhePmzZunG2+8Uddff70WL15cyQr/h7xHIDu6fbskqWHbti5XAsCsqmY97+CXITw8XDk5OW6XUeNkZmZKEr9o+QG+V/6D75V/MPN9CgkJUXh4eHWVVGkZGRmO7XvgwIGqV6+eli5dqtOnT1e5mz55bx0/U/zHwcJCSVJDl+tA+fg/5T+q43tV1aznHXz4FO+lkmVdSgnfwffKf/C98g98n8xr06aNdu3apZ9++knnn3++2+XUOPxb9R98r/wD3yf/4Q/fKzrUAAAAv5GVlaXdu3erQYMGaty4sdvlAADgU1jgAwAAn/LKK6/op59+KrF969atGjJkiAzD0PDhwxUcHOxCdQAA+C7uwQcAAD7lhRde0AMPPKDOnTurXbt2MgxDO3fu1Lp161RYWKgrrrhCzz77rNtlAgDgc1jgAwAAnzJp0iR98sknWrt2rf75z3/q1KlTioiIUJ8+fTR48GDdeuutCgriIkQAAH6JBT4AAPApQ4cO1dChQ90uAwAAv0MXfQAAAAAAAgDXtwEAAAAAEABY4AMAAAAAEABY4AMAAAAAEABY4AMAAAAAEABY4AMAAAAAEABY4AMAAAAAEABY4AMAAAAAEABY4MN1ycnJ8ng8ZX4tXbrU7RLxCwcOHNDDDz+sCy+8UHXq1FFERIS6du2qRx55xO3SIOmrr74q9/+U92vChAlulwpJ3377rW6++WZFRUUpJCRE5513nnr27Knp06fLMAy3ywMqlJOTo5kzZ+ree+9VQkKCateuLY/Ho/Hjx7tdWo106tQpjRs3ThdccIHCwsIUFRWl2267TXv27HG7NJxj3bp1+vOf/6xBgwapZcuWRdkM33Ly5EktWLBAo0eP1oUXXqiwsDCFh4erc+fOmjBhgk6cOOF2iSV4DH57gMuSk5O1fPly3XjjjapXr16Jxx966CFdfPHFLlSG0qxbt05XX321Dh06pLi4OP3qV7/SsWPHtGnTJv3444/Kz893u8Qab8uWLfrzn/9c6mMFBQWaNWuWJOnLL7/UlVdeWZ2l4Rc++ugj3XLLLSooKFDXrl3Vrl07HThwQCtWrFB+fr6GDBmi9957z+0ygXJlZGQoPj6+xPannnqKRX41y83N1ZVXXqnU1FQ1b95cPXv21I4dO7RmzRo1adJEqampiomJcbtMSBowYIAWLlxYYjtLM9/y1ltvacyYMZKkDh06FP3eu2rVKh0/flwXXXSRli9frqZNm7pc6f/UcrsAwGvy5MmKjo52uwyU48CBA+rXr59OnTqlhQsXqn///sUeX7NmjUuV4VwXXXSRZsyYUepjn376qWbNmqVWrVopOTm5WutCcfn5+brrrrtUUFCg9957T0OGDCl6bPPmzbr88ss1e/Zs3X777bwQA59Wv359jR49WpdeeqkuvfRSLVmyROPGjXO7rBrpmWeeUWpqqpKSkvTZZ58VvXHy4osv6qGHHtJtt92mr776yt0iIUlKSkpSp06div7fREdHKy8vz+2y8AshISG644479MADD6hDhw5F2/ft26frrrtO6enpeuCBBzR79mwXqyyOBT4A05566ikdPHhQU6ZMKbG4l6Tu3bu7UBWs8L57P3ToUC4FdNmWLVu0f/9+XXjhhcUW99LZdwmGDRumv/3tb/r2229Z4MOnxcbG6q233ir6+2effeZiNTXX6dOn9eqrr0qSpkyZUuyqyAcffFDvvPOOli9frnXr1qlbt25ulYn/evTRR90uASaMGDFCI0aMKLG9efPmmjJlii677DLNmzdPp0+fVmhoqAsVlsQ9+ABMOXXqlGbNmqXw8HCNGjXK7XJQCTk5OUWXA956660uV4PatWubGhcZGelwJQACwcqVK3X06FHFxsaWesvETTfdJElavHhxdZcGBKTOnTtLkvLy8nTo0CGXq/kf3sGHz5g2bZoOHTqkoKAgXXDBBRowYIBat27tdln4r7Vr1+r48eO6/PLLVadOHX366af6/PPPlZubqwsuuKCoSRh817x585STk6P4+Hh17NjR7XJqvJiYGMXGxmrr1q2aPXt2iUv0Z82apUaNGmngwIEuVgnAX6xfv16S1LVr11If927fsGFDtdUEBLKsrCxJZy/jj4iIcLma/2GBD5/xzDPPFPv7ww8/rCeffFJPPvmkSxXhXJs2bZIkNW3atNTGMI8//rimTZumwYMHu1EeTPBens+7974hODhY77zzjq6//noNHTpUL7zwgtq3b6/9+/drxYoV6tixo2bMmOFTvzQA8F27du2SJLVs2bLUx73bd+7cWW01AYHs5ZdfliT169fP9FV51YFL9OG6K664QjNnzlRmZqZOnjyprVu3atKkSapVq5bGjRtX9J8H7jp8+LAkadGiRVq6dKmmTJmi/fv3a8eOHXr44Yd16tQpjRgxQhkZGe4WilLt27dP//rXvxQcHMyLMD6kR48eWr58uWJiYvTdd99p7ty5WrZsmYKCgtSnTx+6XQMwzftxXXXr1i318fDwcEnS8ePHq60mIFB98sknmjZtmkJCQjRx4kS3yymGd/Dhul9+FvcFF1ygxx9/XJdccomuvvpqjR8/XnfccYfq1KnjUoWQpMLCQklnO39PmjRJd911V9Fjzz//vHbu3Kl//OMfev755/lYLx80Z84cFRQUqF+/fmrWrJnb5eC/5syZo1GjRikxMVFz5sxRXFyc9u7dq8mTJ+uFF17QsmXLtGrVKp96ZwCBZ+DAgdq8ebOl57z77rs0VgVQI23ZskXDhg2TYRh6/vnni+7F9xUs8OGz+vbtq0suuURr165VWloaH+nlsnO78ZbWZG/UqFH6xz/+oeXLl1dnWTCJy/N9zw8//KARI0aoadOm+vjjj4v+j7Vv315Tp07V3r179fHHH+vtt9/W7373O5erRSDbvn27tm7dauk5J0+edKgaVJb3Z0hZ35ucnBxJZz/WEEDl7NmzR/369dPhw4f14IMP6v7773e7pBK4RB8+rX379pLOXl4Md7Vp00bS2Uv/mjRpUuLx6OhoSdL+/fursyyYsHnzZqWnp6tevXoaMGCA2+Xgv95//32dOXNG/fr1K/YCmtfNN98sSfr666+ruzTUMBkZGTIMw9IXL7r7Hm9j4h9//LHUx73bvXkOwJrs7Gz17dtXO3fu1KhRozR58mS3SyoVC3z4NO993977xuAe70funDp1Snl5eSUez87OlqRSFypw18yZMyVJgwYNKvPeTFQ/7y/bDRs2LPVx73bvz0EAKI/3MuHvvvuu1Me92zt16lRtNQGB4sSJE7rmmmu0adMmDRo0SG+++aY8Ho/bZZWKBT581oEDB7RixQpJZX/kC6pP69at1blzZxmGUepl+N5tpX32LtxjGIZmz54ticvzfY23F8LatWtLffzbb7+V9L+rYwCgPD169FDDhg2VmZlZasPbDz/8UJJ0ww03VHNlgH/Ly8tTSkqK1qxZo6uvvlpz5sxRcHCw22WViQU+XLVq1SotWLBABQUFxbbv2LFDAwcOVE5Ojvr371/mR76gev3hD3+QdPYjDM+9bSIjI0MvvPCCJOnOO+90pTaUbsWKFdq5c6datGih3r17u10OzpGSkiLp7CX4r732WrHHUlNT9dJLL0mSbrrppmqvDYD/CQ0N1T333CNJuvvuu4vuuZekF198URs2bFCvXr3UrVs3t0oE/E5BQYEGDx6sL7/8Uj179tS8efMUGhrqdlnl8hiGYbhdBGquGTNmaNSoUWrWrJm6du2q8847Tzt37tS6deuUm5uruLg4ffnll2ratKnbpeK/Ro4cqXfeeUfnnXeeLrvsMp06dUqrVq1SXl6exowZozfeeMPtEnGOO+64Q2+++aYeeeQRPffcc26Xg1945JFHiu7hi4uLU8eOHbV3716tXr1ahYWFuuOOOzR16lSXqwQqNnDgwKIXfvfu3avdu3erRYsWRS/QN2/eXPPnz3ezxBohNzdXycnJSktLU/PmzdWzZ0/t3LlTaWlpatKkiVJTU/n4TR+xZMmSYh+vtmbNGhmGoYSEhKJtTz75pK677jo3ysN/vfzyy3rggQcknf0516BBg1LHTZ48WY0bN67GysrGAh+u2rx5s1555RWlpaVp9+7dOnz4sMLDw9WhQwf95je/0e9+9zs+Hs/HGIaht956S1OnTtXmzZvl8XjUqVMnjR07ViNGjHC7PJwjLy9PzZs31+HDh7V+/Xruu/RR8+fP1+uvv65169bp6NGjql+/vrp06aIxY8Zo8ODBbpcHmBIdHa2dO3eW+XibNm20Y8eO6iuoBjt16pSeffZZzZ49W7t371ZERIT69euniRMnckWkD/G+yVWe6dOna+TIkdVTEEo1fvx4Pf300xWO2759u8/cUscCHwAAAACAAMA9+AAAAAAABAAW+AAAAAAABAAW+AAAAAAABAAW+AAAAAAABAAW+AAAAAAABAAW+AAAAAAABAAW+AAAAAAABAAW+AAAAAAABAAW+AAAAAAABAAW+AAAAAAABAAW+AAAAAAABAAW+AAAAAAABAAW+KhRli1bphtvvFEtWrRQaGioGjVqpAsvvFC/+c1v9Oqrr+ro0aNul2hadHS0PB6P6fEzZsyQx+PR+PHjnSuqmn311VfyeDwaOXKkrfvNz8/XhRdeqO7du9u6X6cZhqH4+HhdfPHFKiwsdLscAHANeU/em0HeIxCxwEeNMWHCBPXu3Vvz5s1Tw4YNdf3116tv376qU6eO5s2bp3vvvVebN292u0ycw61fUqZOnar//Oc/fvfLkcfj0bhx4/Tvf/9bM2bMcLscAHAFee9/yHtryHuUp5bbBQDVYd26dRo/frxCQkL0wQcfaMCAAcUe/+mnnzRr1iydd955rtQH35GXl6cJEyaoU6dOuvbaa90ux7IBAwbooosu0rhx4zR8+HDVqsWPeQA1B3kPs8h7BCrewUeNMG/ePBmGoZtvvrlE2EtSs2bN9PDDD+uiiy6q/uLgUz788EPt379fw4cPd7uUSvF4PBo6dKj27NmjRYsWuV0OAFQr8h5mkfcIVCzwUSMcOHBAktSkSRPLzz158qSeffZZxcfHq169eqpXr54SExP1zjvvlDre4/EoOjpap0+f1lNPPaXY2FiFhYUpJiZG48aNU25ubonnbNu2TePHj1dSUpKaNWum0NBQtWzZUsOHD9d//vMfyzVbZRiG5syZo969e6tRo0YKCwtThw4dNH78eJ08ebLE+OTkZHk8Hu3YsUMLFixQYmKiwsPDFRERocGDB+vHH38sdZ6DBw/qd7/7naKiolSnTh396le/0pQpU2QYRtF5O3eOUaNGSZKefvppeTyeoq/SLkfLzs7W7373OzVv3ly1a9fWr371K7399tuWz8Vbb70lj8ej3/72tyUeO/cSwszMTN18881q3LixGjRooGuuuUabNm2SdPaevj/96U+64IILFBYWpnbt2mnKlCkl9rdjxw55PB4lJycrJydHDz74oFq1aqU6deqoa9euWrx4cdHYf/zjH0pISFB4eLjOP/983XfffTp16lSpxzBkyBBJ0ptvvmn5+AHAn5H35SPv/4e8R8AygBpgwoQJhiSjVatWxs8//2z6eT///LPRqVMnQ5LRrFkz49prrzWuueYao2HDhoYk45577inxHElG69atjeuvv96oU6eOcf311xuDBg0qes5VV11l5OfnF3vOo48+ang8HuPiiy82rr/+euPGG280OnToYEgyGjRoYKxfv77EPG3atDGs/BeePn26Icl46qmnim0vKCgwBg8ebEgy6tWrZyQnJxsDBw40WrVqZUgyunfvbpw8ebLYc3r16mVIMh555BEjODjYSE5ONm666aai57Rv377Ecw4cOGC0b9/ekGRERUUZN998s9GnTx8jJCTEuP/++w1JRps2bYrGP/vss0aPHj0MSUbnzp2NESNGFH2tWLHCMAzDWLZsmSHJSElJMS644AIjKirK+M1vfmNceeWVRnBwsCHJePPNN02fo6NHjxrBwcFG+/btyz2Hw4cPNyIiIowOHToYt9xyi3HxxRcbkowmTZoY+/btM1JSUoyGDRsaAwYMMK6++mojNDTUkGS88cYbxfa3fft2Q5KRlJRkJCQkGE2bNjVuuukmIzk52QgKCjKCg4ONzz//3HjxxReNWrVqGVdddZUxcOBAIzIy0pBkDBkypMxjadWqlREaGlri+wAAgYy8J+/NIO8RyFjgo0bIzMw06tSpY0gy6tevb4wYMcJ48803je+++65E+J7r2muvNSQZ999/v5Gbm1u0/aeffjIuueQSQ5Lx6aefFnuOJEOS0bJlSyMzM7No+/79+41f/epXhiTjpZdeKvac1atXG1lZWSXmf/vttw1JxpVXXlniMbsC/7nnnjMkGcnJyca+ffuKtufl5RmjR482JBmPPvposed4A79u3brGqlWrirbn5OQYl112mSHJmDZtWrHnePfVv39/49SpU0Xb161bV/TL0LmBX17NXt7Al2T89re/LfY9mj9/ftEvX2Z9+umn5Qaptx5Jxv/93/8ZhYWFhmEYRmFhoTFy5EhDktGxY0fjV7/6lbF///6i533xxRelHp838CUZvXv3Nk6cOFFirnbt2hmNGjUyvv3226LH9uzZYzRt2tSQVOzf2LluvPFGQ5Lx5Zdfmj5+APB35D15bwZ5j0DGAh81xhdffFH0ivO5X+edd57xu9/9zti7d2+x8enp6YYk49JLLzUKCgpK7O+7774rCrBzeff7y1dvDeN/gRIbG2u67h49ehgej8c4cuRIse12BP6ZM2eMxo0bG+Hh4cZPP/1U4jknT540mjVrZjRq1KjYOfAG/h//+McSz/nwww8NScaIESOKth0/ftwICwszgoODje3bt5d4zh//+McqBX6DBg2MgwcPlnjc+wtWaXOW5i9/+YshyXj66adLfdxbT0xMjHH69Olij61fv77oe//FF1+UeG58fHyJWryBHxQUZGzdurXY+IKCAqNx48aGJOOJJ54osb/f//73hiRj+vTppdbqPacvvvhiBUcNAIGFvCfvK0LeI5BxDz5qjKuuukrbtm3TvHnzdOedd6pr166qVauWjhw5otdee01dunTR1q1bi8Z/9tlnks52KQ0KKvlfxXuP3po1a0qdr7R7uvr166dGjRopMzNT+/btK/bYiRMnNGfOHD366KMaM2aMRo4cqZEjR2rfvn0yDEOZmZlVOfxSfffddzp48KAuu+wynX/++SUer1Onjrp166bDhw/rhx9+KPF43759S2y74IILJKnY8a1bt065ubm69NJLi91353XLLbdU4Sikbt26KTIy0lQt5dm/f78kqVGjRuWOS05OVkhISLFtMTExkqSQkBAlJyeXeI738dJqiY6OLqrVKygoSG3atJFU+nkub3+SFBERIel/96MCQE1B3pdE3hdH3iOQ8XkKqFFCQ0M1cOBADRw4UJJ05MgRvf/++3r88ce1f/9+3XPPPfr8888lnW2IIkl//OMf9cc//rHMfZbWRKdRo0aqX79+qePbtGmjw4cPa+/evWrevLkk6csvv9Rvf/vbcn84Hz9+3NQxWuE9xs8//1wej6fcsQcPHtSFF15YbFvLli1LjPMed15eXtE2byi1atWq1H23bt3adM2lKa2Osmopz9GjR4s9rywtWrQosa1evXqSznZoDg4OLvPx0mopbX/nPqe8+co6tgYNGkg6+28cAGoa8r448r448h6BjAU+arTzzjtPd955p6KiopSSkqJly5bp5MmTqlu3rgoLCyVJl19+uWJjYx2r4cSJE7r55puVnZ2tcePG6be//a3atGmjOnXqyOPxaMiQIZozZ44Mw7B9bu8xtmvXTj169Ch3bGmvmJf2Tocb7KqjYcOGkir+5aq8+SpTS0XPqcw+vb+88FnPAEDek/fFkfcIZCzwAUm9e/eWJBUUFOjIkSOqW7du0avEAwYM0EMPPWRpf4cPH9bx48dLfWV4165dkqSoqChJ0ooVK3To0CHddNNNevrpp0uMz8rKsjS3Fd5jvOiii0r9KBq7eN+52L17d6mPl7W9ujVt2lTS2Y/g8XeHDx+WVLmPigKAQEXek/cSeY/A5hsvxwEOq+jV8G3btkk6e0lf48aNJUl9+vSRJM2fP79Sc37wwQcltn322WfKzs5WTExMUQh6fzCXdtnZtm3b9N1331VqfjMuvfRSNWzYUMuXL3c05Lp166awsDCtXbu26Beec5V2rqSz3w/p7OfMVofOnTtLUrF7M/3V5s2bJUldunRxtxAAqEbkfenI++LIewQyFvioEZ588kk98sgjpTau2bNnj8aOHStJ6t+/f1HIJCQkqE+fPlq5cqXuvvtuHTt2rMRz169fr6VLl5Y659NPP110z5t09p62Rx55RJJ09913F233NluZN29esXvyjhw5otGjR+vMmTMWj9a82rVr6w9/+IOOHz+uQYMGlfruwZ49ezRz5swqzVOvXj0NHTpU+fn5uv/++4vdR7Z+/Xq98sorpT7P+65HdQXwZZddpuDgYH377bfVMp+T1qxZo9DQUCUmJrpdCgBUG/K+dOR9ceQ9AhmX6KNGOHHihF5++WVNnjxZF1xwgTp27KiwsDD9+OOPSktL05kzZ9SuXTv99a9/Lfa8WbNmqV+/fvr73/+u2bNnq0uXLoqKitLRo0e1YcMG7d69W/fff7/69etX7HmtW7dWp06dFBcXp6uuukohISH68ssvdeTIEV155ZW67777isZecskl6tOnjz7//HNdcMEFRR1Zv/rqKzVu3FgpKSlauHChY+fm//7v/7RlyxbNnDlTHTp0UHx8vNq2bavTp09r69at2rRpkzp16qRbb721SvP8+c9/1vLly7VgwQLFxsbq8ssv15EjR/Tll19q7NixevXVV4t+2fJKTExU06ZN9eGHHyo5OVkxMTEKCgrSbbfdpssuu6xK9ZSmfv366tmzp7766iv9+OOPZTbz8XWZmZn68ccf1a9fP9WpU8ftcgCg2pD3ZSPv/4e8RyDjHXzUCE888YRmzpypYcOGqXbt2lqxYoU+/PBDbdq0Sd27d9dzzz2njIyMEt1LmzZtqlWrVulvf/ubOnbsqPT0dH344YfasGGDYmJi9Pzzz+vhhx8uMZ/H49GHH36oBx54QN9//70+/vhjNWzYUH/84x+1ZMkS1apV/LW1hQsX6o9//KOaNGmiTz/9VOvWrdNvf/tbpaamOt40JSgoSO+++64WLlyoPn36aPv27froo4/0zTffKCwsTI888ojefvvtKs/TuHFjrVq1SmPHjlVBQYEWLFigXbt26fnnn9cf/vAHSSUb+4SFhWnJkiXq06ePMjIyNGPGDE2bNk3/+c9/qlxPWcaMGSNJmjNnjmNzOG327NmS/ncsAFBTkPdlI++LI+8RqDyGE606gRrM4/GoTZs2xS7XQ/nef/99DR48WHfeeadee+01V2vJy8tTmzZt1LRpU23YsMHVWirDMAx16NBBJ06c0I4dO0r8cgkAsAd5bx15bx/yHmXhHXwA1WbdunUltmVkZBTdqzhs2LDqLqmE2rVra9y4cUXvxPibBQsWaOvWrZowYQJhDwBwBXnvPPIeZeEdfMBmvKJftrCwMDVr1kwdOnRQgwYNtH37dq1bt06FhYW65557ymy+U93y8/MVFxenBg0a+FUDHsMw1LVrV+Xn52v9+vU+87nFABCIyPuykffOIu9RHhb4gM0I/LI9/fTT+uSTT5SVlaUjR46oXr166tKli26//XYNHTrU7fIAADCNvC8beQ+4hwU+AAAAAAABgOs5AAAAAAAIACzwAQAAAAAIACzwAQAAAAAIACzwAQAAAAAIACzwAQAAAAAIACzwAQAAAAAIACzwAQAAAAAIACzwAQAAAAAIACzwAQAAAAAIACzwAQAAAAAIACzwAQAAAAAIACzwAQAAAAAIACzwAQAAAAAIAP8Pn0qcp/IGowYAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 1200x450 with 2 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.figure(figsize = (8, 3))\n",
    "\n",
    "plt.subplot(121)\n",
    "plt.plot(sepal_length, petal_length, \"o\", markersize=4)\n",
    "plt.xlabel(\"Sepal length (mm)\")\n",
    "plt.ylabel(\"Petal length (mm)\")\n",
    "\n",
    "plt.subplot(122)\n",
    "plt.plot(sepal_length_centered, petal_length_centered, \"o\", markersize=4)\n",
    "plt.axvline(c=\"darkred\", lw=0.75)\n",
    "plt.axhline(c=\"darkred\", lw=0.75)\n",
    "plt.xlabel(\"Sepal length (mm)\")\n",
    "\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.7.3"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}