{ "cells": [ { "cell_type": "code", "execution_count": 1, "id": "5c32b417", "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "import matplotlib.pyplot as plt\n" ] }, { "cell_type": "code", "execution_count": 14, "id": "4acf77bc", "metadata": {}, "outputs": [], "source": [ "a = 2\n", "b = 1\n", "x = np.linspace(0,10,10)\n", "y = a*x + b" ] }, { "cell_type": "code", "execution_count": 35, "id": "d967b4dc", "metadata": {}, "outputs": [], "source": [ "y_bruit = y + np.random.randn((10))*2" ] }, { "cell_type": "code", "execution_count": 36, "id": "b4a2ef7d", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 36, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAh8AAAGdCAYAAACyzRGfAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAA9hAAAPYQGoP6dpAABFwElEQVR4nO3dd3hUZeL28e+k90AIpEACoYZAggULwS4KSJFmZcVVV1fFgqiAihQBQ3Ata1/97br4rmVZuqggIogiKKKShBBq6JBQMymkzZz3j2g0EoHAZM6U+3NduTTnnMzcGULm5pzzPI/FMAwDERERESfxMTuAiIiIeBeVDxEREXEqlQ8RERFxKpUPERERcSqVDxEREXEqlQ8RERFxKpUPERERcSqVDxEREXEqP7MD/J7dbmffvn2Eh4djsVjMjiMiIiKnwTAMiouLiY+Px8fn5Oc2XK587Nu3j4SEBLNjiIiIyBnYvXs3rVq1OukxLlc+wsPDgZrwERERJqcRERGR02G1WklISKh9Hz8Zlysfv1xqiYiIUPkQERFxM6dzy4RuOBURERGnUvkQERERp1L5EBEREadyuXs+TodhGFRXV2Oz2cyOIk7m6+uLn5+fhmGLiLgxtysflZWV7N+/n7KyMrOjiElCQkKIi4sjICDA7CgiInIG3Kp82O128vPz8fX1JT4+noCAAP0L2IsYhkFlZSUHDx4kPz+fDh06nHIiGxERcT1uVT4qKyux2+0kJCQQEhJidhwxQXBwMP7+/uzcuZPKykqCgoLMjiQiIg3klv9s1L92vZv+/EVE3JtbnfkQERExhd0GO7+BkgIIi4HW6eDja3Yqt6Xy4aLmz5/P8ePHueWWWxz+2JWVlfztb39j8ODBdO7c2eGPLyLiUXIXwuKxYN3367aIeOiTCSkDzcvlxnT+2gWtWbOGhx56iB49epzy2BUrVmCxWDh27NhpP/6jjz5KdnY2ycnJZ5HyzP373/+mSZMmpjy3iEiD5C6EWSPqFg8A6/6a7bkLzcnl5lQ+XMzhw4e56667mD9/Pm3atDnl8enp6ezfv5/IyEjg1G/ss2bNYsOGDcycOdO0kUI33XQTmzdvNuW5RUROm91Wc8YDo56dP29bPK7mOHdht0H+V5A9u+a/JmXXZRcX06xZMzZs2HBax1ZVVREQEEBsbOxpP/6NN97IjTfeeKbxTqqysvK05t4IDg4mODi4UTKIiDjMzm9OPONRhwHWvTXHJV3qtFhnzIUuH7n9mQ/DMCirrDblwzDqa8Mneuutt4iPj8dut9fZfv3113PnnXfWfr5gwQLOO+88goKCaNu2LZMnT6a6urp2v8Vi4Y033mDgwIGEhoYybdq0OpddVqxYwR133EFRUREWiwWLxcKkSZMAqKio4LHHHqNly5aEhoZy0UUXsWLFij/MfOutt3LTTTfV2VZVVUV0dDTvvvsuAFdccQUPPPAAo0aNIjo6mt69ewPwwgsvkJqaSmhoKAkJCdx///2UlJTUPk59Z2dO9b2LiDhdSYFjjzOTi10+cvszH8erbKRMWGLKc+c+05uQgFO/hDfccAMPPvggy5cv5+qrrwbgyJEjLF68mE8++QSAr776ihEjRvDyyy9z6aWXsm3bNu655x4AJk6cWPtYkyZNYvr06bz00kv4+fmxffv22n3p6em89NJLTJgwgU2bNgEQFhYGwAMPPEBubi4ffvgh8fHxzJs3jz59+pCdnU2HDh1OyDx8+HBuuOEGSkpKah9jyZIllJWVMXjw4NrjZs6cyX333ceqVatqt/n4+PDyyy+TlJTE9u3buf/++xkzZgyvv/56va/P6X7vIiJOFRbj2OPMcsrLR5aay0fJ/Zw2gsftz3y4g6ZNm9K3b1/ef//92m2zZ88mOjqaK6+8EoDJkyczbtw4br/9dtq2bcs111zDlClT+Mc//lHnsW699VbuuOMO2rZtS2JiYp19AQEBREZGYrFYiI2NJTY2lrCwMHbt2sU777zD//73Py699FLatWvHY489xiWXXMI777xTb+bevXsTGhrKvHnzare9//77DBw4kPDw8NptHTp0YMaMGXTq1IlOnToBMGrUKK688kratGnDVVddxdSpU5k1a9Yfvj6n+72LiDhV6/SayxL80f1xFohoWXOcK2vI5SMncfszH8H+vuQ+09u05z5dw4cP5+677+b1118nMDCQ9957j5tvvrl2wqz169ezatUqpk2bVvs1NpuN8vJyysrKamd07d69e4NzZmdnY7PZ6NixY53tFRUVNGvWrN6v8fPz48Ybb+S9997jtttuo7S0lAULFvDhhx/WOe78888/4Ws///xzMjIyyMvLw2q1Ul1dfcL38Vun+72LiDiVj2/N/RCzRlBTQH575uDnQtJnuuvP9+GCl4/cvnxYLJbTuvRhtgEDBmAYBh9//DEXXHABX331FS+++GLt/pKSEiZPnsyQIUNO+NrfTiEeGhra4OcuKSnB19eXdevW4etb9y/JL5dU6jN8+HAuv/xyCgsLWbp0KcHBwfTp06fOMb/Ps2PHDvr37899993HtGnTiIqK4uuvv+auu+6isrKy3iJxut+7iIjTpQyEG9/9gxs1p7vHPB8uePnI9d+1PURQUBBDhgzhvffeY+vWrXTq1Inzzjuvdv95553Hpk2baN++/Vk9T0BAADZb3aFT5557LjabjcLCQi699PTvyE5PTychIYH//ve/fPrpp9xwww34+/uf9GvWrVuH3W7n+eefrz2rc7JLLuC4711EpFGkDKy5H8JdZzj95fKRdT/13/dhqdnvxMtHKh9ONHz4cPr378+GDRv405/+VGffhAkT6N+/P4mJiQwbNgwfHx/Wr19PTk4OU6dOPe3naNOmDSUlJSxbtoxu3boREhJCx44dGT58OCNGjOD555/n3HPP5eDBgyxbtoy0tDT69ev3h49366238uabb7J582aWL19+yudv3749VVVVvPLKKwwYMIBVq1bx5ptvnvRrHPW9i4g0Gh9f9xhOWx8XvHykG06d6KqrriIqKopNmzZx66231tnXu3dvFi1axGeffcYFF1zAxRdfzIsvvkjr1q0b9Bzp6ence++93HTTTTRv3pwZM2YA8M477zBixAgeffRROnXqxKBBg1i7du0JN63+3vDhw8nNzaVly5b07NnzlM/frVs3XnjhBTIzM+natSvvvfceGRkZJ/0aR33vIiLyB365fBQRV3d7RHzNdidfPrIYpztZhZNYrVYiIyMpKioiIiKizr7y8nLy8/NJSkrSvQBeTD8HIiJnqBEXyDvZ+/fv6bKLiIiIt3CRy0e67CIiIiJOpfIhIiIiTqXyISIiIk6l8iEiIiJOpfIhIiIiTqXyISIiIk6l8iEiIiJOpfLhAlasWIHFYuHYsWNmRxEREWl03jvJWCPO8iYiIiJ/zDvLR+7CP1geOdM9lkcWERFxY9532SV3Yc3Kfr8tHlCz1PCsETX7G0FFRQUPPfQQLVq0ICgoiEsuuYS1a9fWOWbVqlWkpaURFBTExRdfTE5OTu2+nTt3MmDAAJo2bUpoaChdunThk08+aZSsIiIijcm7yofdVnPGg/rW0vt52+JxNcc52JgxY5gzZw4zZ87khx9+oH379vTu3ZsjR47UHvP444/z/PPPs3btWpo3b86AAQOoqqoCYOTIkVRUVLBy5Uqys7PJzMwkLCzM4TlFREQam3eVj53fnHjGow4DrHtrjnOg0tJS3njjDZ577jn69u1LSkoKb7/9NsHBwfzzn/+sPW7ixIlcc801pKamMnPmTAoKCpg3bx4Au3btomfPnqSmptK2bVv69+/PZZdd5tCcIiIizuBd5aOkwLHHnaZt27ZRVVVFz549a7f5+/tz4YUXsnHjxtptPXr0qP3/qKgoOnXqVLv/oYceYurUqfTs2ZOJEyeSlZXl0IwiIiLO4l3lIyzGscc50V/+8he2b9/ObbfdRnZ2Nt27d+eVV14xO5aIiEiDeVf5aJ1eM6oFyx8cYIGIljXHOVC7du0ICAhg1apVtduqqqpYu3YtKSkptdvWrFlT+/9Hjx5l8+bNdO7cuXZbQkIC9957L3PnzuXRRx/l7bffdmhOERERZ/CuobY+vjXDaWeNoKaA/PbG058LSZ/pDp/vIzQ0lPvuu4/HH3+cqKgoEhMTmTFjBmVlZdx1112sX78egGeeeYZmzZoRExPDU089RXR0NIMGDQJg1KhR9O3bl44dO3L06FGWL19ep5iIiIicjqLjVdjtBk1DA0zL4F3lA2rm8bjx3T+Y52N6o83zMX36dOx2O7fddhvFxcV0796dJUuW0LRp0zrHPPzww2zZsoVzzjmHjz76iICAmh8Om83GyJEj2bNnDxEREfTp04cXX3yxUbKKiIhnWppbwPj52VyY1IxXbjnXtBwWwzDqG3dqGqvVSmRkJEVFRURERNTZV15eTn5+PklJSQQFBZ3dE2mGU7fl0J8DEREvcLikgkkf5fLR+pp/dCdFhzL//p5Ehvg77DlO9v79e9535uMXPr6QdKnZKURERBqNYRgsXL+PyR/lcqS0Eh8L3H1ZWx7p1ZEgf/P+we295UNERMSD7S86zvh5OSzLKwQgOTacGcPSSGvVxNxgqHyIiIh4FMMw+OC73WR8spHiimr8fS08eFUH7r28HQF+rjHIVeVDRETEQ+w8XMq4Odms3n4YgHMSmjBjWBodY8JNTlaXyoeIiIibs9kN3lmVz98+20R5lZ0gfx8eu7YTd/RMwtfnj+a2Mo9blg8XG6AjTqY/fxGRX20uKGbM7Cx+2n0MgPR2zZg+JI3EZiHmBjsJtyof/v41Q4LKysoIDg42OY2YpaysDPj150FExBtVVtt5Y8U2Xl2+hSqbQXigH0/268zNFyRgsbje2Y7fcqvy4evrS5MmTSgsrLlzNyQkxOVfYHEcwzAoKyujsLCQJk2a4OureVlExDtl7TnGmNlZ5B0oBuDq5BZMHdyVuEj3+Ie5W5UPgNjYWIDaAiLep0mTJrU/ByIi3qS8ysaLSzfz9lfbsRsQFRrAxAEpDOwW71b/GHe78mGxWIiLi6NFixZUVVWZHUeczN/fX2c8RMQrfbv9MOPmZpN/qBSAgd3imTgghWZhgSYna7gGlY+MjAzmzp1LXl4ewcHBpKenk5mZSadOnWqPKS8v59FHH+XDDz+koqKC3r178/rrrxMT49hl6n19ffUmJCIiHq+4vIrMxXn8Z80uAGIiApk2KJVeKY59X3WmBs028uWXXzJy5EjWrFnD0qVLqaqq4tprr6W0tLT2mEceeYSPPvqI//3vf3z55Zfs27ePIUOGODy4iIiIp1u+qZDeL66sLR63XJjAZ49c7tbFA85yYbmDBw/SokULvvzySy677DKKiopo3rw577//PsOGDQMgLy+Pzp07s3r1ai6++OJTPmZDFqYRERHxREdLK5myKJe5P+4FIDEqhOlDUklvH21ysj/mtIXlioqKAIiKigJg3bp1VFVV0atXr9pjkpOTSUxM/MPyUVFRQUVFRZ3wIiIi3uqT7P1MWJDDoZJKLBa4s2cSj17bkZAAt7tN8w+d8Xdit9sZNWoUPXv2pGvXrgAcOHCAgIAAmjRpUufYmJgYDhw4UO/jZGRkMHny5DONISIi4hEKreVMWLCBxRtq3i87tAgjc1ga5yU2NTmZ451x+Rg5ciQ5OTl8/fXXZxXgiSeeYPTo0bWfW61WEhISzuoxRURE3IVhGPxv3R6mLsrFWl6Nn4+F+69ox8ir2hPo55kDK86ofDzwwAMsWrSIlStX0qpVq9rtsbGxVFZWcuzYsTpnPwoKCv5wXobAwEACA91vmJCIiMjZ2n2kjCfnZfPVlkMApLaMJHNoGinxnn3PY4PKh2EYPPjgg8ybN48VK1aQlJRUZ//555+Pv78/y5YtY+jQoQBs2rSJXbt20aNHD8elFhERcWN2u8H/W7OTzMV5lFXaCPTz4ZFrOvKXS5Lw83WNZe8bU4PKx8iRI3n//fdZsGAB4eHhtfdxREZGEhwcTGRkJHfddRejR48mKiqKiIgIHnzwQXr06HFaI11EREQ83baDJYydncX3O48CcGGbKKYPTaVt8zCTkzlPg4ba/tHUre+88w5//vOfgV8nGfvggw/qTDJ2utNha6itiIh4oiqbnbdWbufvy7ZQWW0nNMCXcX2TGX5Ra3xccNn7hmrI+/dZzfPRGFQ+RETE02zYV8SY2Vls2FczncRlHZvz7OCutGrqusveN5TT5vkQERGRP1ZeZePVL7by5pfbqLYbRAb7M6F/CkPOa+lWC8E5msqHiIhII1i38whjZmex7WDNEiTXpcYyeWBXmodrhKfKh4iIiAOVVlTz3JJNzFy9A8OA6LBApg7qQp+ucWZHcxkqHyIiIg7y9ZZDjJubxZ6jxwEYdn4rxvfrTJOQAJOTuRaVDxERkbNUdLyKaR/nMuv7PQC0bBLMs0NSubxjc5OTuSaVDxERkbPw2YYDjJ+fQ2FxzSKpt/dozeN9kgkL1FvsH9ErIyIicgYOlVQwceEGPs7aD0Db6FAyh6VxQZsok5O5PpUPERGRBjAMgwU/7WPyRxs4WlaFr4+Fey5ry8NXdyDI3zMXgnM0lQ8REZHTtO/YccbPz+GLvEIAOsdFMGNoGqmtIk1O5l5UPkRERE7Bbjf4YO0uMj7Jo6SimgBfHx66uj1/vbwd/l6wEJyjqXyIiIicxI5DpYybm8Wa7UcAODexCTOGptEhJtzkZO5L5UNERKQeNrvBP7/ezvOfbaai2k6wvy+P9e7En9Pb4OsBC8GZSeVDRETkdzYdKGbM7PWs31MEQM/2zcgYnEZiM89ZCM5MKh8iIiI/q6y28/qKrby2fCtVNoPwID/G9+vMjd0TvHohOEdT+RAREQHW7z7GmNlZbCooBqBX5ximDe5KTESQyck8j8qHiIh4teOVNl78fDP/99V27AY0Cw1g0sAu9E+L09mORqLyISIiXmvN9sOMm5PFjsNlAFx/TjwTB3QhKlQLwTUmlQ8REfE6xeVVTP80j/e+3QVAbEQQ0wZ35erOMSYn8w4qHyIi4lWW5xXy5Lxs9heVA3DrRYmM65tMRJC/ycm8h8qHiIh4haOllTyzKJd5P+4FoHWzEDKGpJLeLtrkZN5H5UNERDyaYRh8nL2fiQs2cLi0Eh8L3NkziUev7URwgBaCM4PKh4iIeKxCaznj5+fwWW4BAB1jwsgcmsa5iU1NTubdVD5ERMTjGIbB/77fw5SPcykur8bPx8LIK9sz8sr2BPhpITizqXyIiIhH2X2kjCfnZfPVlkMApLWKZMawNJJjI0xOJr9Q+RAREY9gtxvMXL2D55ZsoqzSRqCfD49e25E7eybhp2XvXYrKh4iIuL2thSWMnZPFup1HAbgwKYrMoWkkRYeanEzqo/IhIiJuq8pm562V2/n751uotNkJC/RjXN9kbr0wER8te++yVD5ERMQt5ewtYszsLHL3WwG4olNznh2cSnyTYJOTyamofIiIiFspr7Lx8rIt/GPldmx2gyYh/kwckMKgc1pqITg3ofIhIiJuY93OI4yZncW2g6UA9EuNY9LALjQPDzQ5mTSEyoeIiLi80opqnluyiZmrd2AY0Dw8kCnXd6VP11izo8kZUPkQERGX9tWWgzwxN5s9R48DcMP5rRjfL4XIEC0E565UPkRExCUVlVUx7ZNcZn2/B4CWTYKZPjSVSzs0NzmZnC2VDxERcTlLNhxg/PwcDhZXYLHA7T3a8HjvToQG6m3LE+hPUUREXMbB4gomLdzAx9n7AWjbPJQZQ9Po3ibK5GTiSCofIiJiOsMwmP/TXiZ/lMuxsip8fSzce3lbHryqA0H+Wvbe06h8iIiIqfYdO85T87JZvukgAClxEcwYlkbXlpEmJ5PGovIhIiKmsNsN3v9uF9M/zaOkopoAXx8e7tWBey5ri78WgvNoKh8iIuJ0+YdKGTcni2/zjwBwfuumZA5No32LMJOTiTOofIiIiNNU2+z8a1U+z3+2mYpqOyEBvozp3YnberTBVwvBeQ2VDxERaTx2G+z8BkoK2FUZzsPfBPHj3hIALmkfTcaQVBKiQkwOKc6m8iEiIo0jdyEsHgvWfQAkAq8ZUcwIuoP0fndwQ/dWWgjOS6l8iIiI4+UuhFkjMDD4bb2IsxzhRV7AEnoOWBLMSicm0+3EIiLiWHYb9k/HnlA8ACw/f7B4XM0lGfFKKh8iIuJQG1Yvxqd43wnF41cGWPfW3AsiXkmXXURExCGKy6vI+DSPkrXf8HLAaXxBSUGjZxLXpPIhIiJnbXleIU/Oy2Z/UTkX+zQ5vS8Ki2nUTOK6dNlFRETO2JHSSkZ9+CN3/Hst+4vKad0shIfvuB0i4uEPL7xYIKIltE53ZlRxITrzISIiDWYYBouy9jNp4QYOl1biY4G/XNqWR3p1JDjAF/pkwqwR1BQQ4zdf+XMh6TMdfLRgnLdS+RARkQYpsJYzfn4OS3Nr7tnoFBPOjGFpdEto8utBKQPhxnfrzPMB1JwR6TO9Zr94LZUPERE5LYZhMOv73Uz9eCPF5dX4+1oYeWV77r+iPQF+9VzFTxkIyf1qZzglLKbmUovOeHg9lQ8RETml3UfKeGJuNl9vPQRAt1aRzBjWjU6x4Sf/Qh9fSLrUCQnFnah8iIjIH7LZDWZ+s4PnlmzieJWNQD8fHru2E3dekqSF4OSMqXyIiEi9thYWM2Z2Fj/sOgbARUlRZA5No010qLnBxO2pfIiISB1VNjv/+HIbLy/bSqXNTligH09cl8wtFyTio7Md4gAqHyIiUitnbxFjZmeRu98KwJWdmjNtcCrxTYJNTiaeROVDREQor7Lx92VbeGvldmx2g6Yh/kwc0IXrz4nXsvficCofIiJebu2OI4ydncX2Q6UA9EuLY/LALkSHBZqcTDyVyoeIiJcqrahmxuI83l2zE8OAFuGBTBnUld5dYs2OJh5O5UNExAut3HyQJ+Zms/fYcQBu6p7Ak9d1JjLE3+Rk4g1UPkREvEhRWRVTPs5l9ro9ALRqGsz0IWlc0iHa5GTiTVQ+RES8xOKc/Ty9YAMHiyuwWODP6W147NpOhAbqrUCcq57J+E9u5cqVDBgwgPj4mjug58+fX2f/n//8ZywWS52PPn36OCqviIg00MHiCu5/bx33/ucHDhZX0K55KLPv7cHEAV1UPMQUDf6pKy0tpVu3btx5550MGTKk3mP69OnDO++8U/t5YKDumBYRcTbDMJj3416eWZTLsbIqfH0s3Ht5Wx68qgNB/lrcTczT4PLRt29f+vbte9JjAgMDiY3V3dIiImbZe+w4T83LZsWmgwCkxEUwY1gaXVtGmpxMpJHu+VixYgUtWrSgadOmXHXVVUydOpVmzZrVe2xFRQUVFRW1n1ut1saIJCLiFex2g/e+28X0TzZSWmkjwM+Hh6/uwD2XtcXft8FX2kUahcPLR58+fRgyZAhJSUls27aNJ598kr59+7J69Wp8fU88zZeRkcHkyZMdHUNExOvkHypl7Jwsvss/AkD31k2ZPjSN9i3CTE4mUpfFMAzjjL/YYmHevHkMGjToD4/Zvn077dq14/PPP+fqq68+YX99Zz4SEhIoKioiIiLiTKOJiHiNapudf36dzwtLN1NRbSckwJcxvTsxokcbLQQnTmO1WomMjDyt9+9Gv825bdu2REdHs3Xr1nrLR2BgoG5IFRE5Qxv3Wxk7J4usPUUAXNohmmcHp5IQFWJyMpE/1ujlY8+ePRw+fJi4uLjGfioREa9RUW3jteXbeH35VqrtBhFBfozvn8IN57fSQnDi8hpcPkpKSti6dWvt5/n5+fz0009ERUURFRXF5MmTGTp0KLGxsWzbto0xY8bQvn17evfu7dDgIiLe6sddRxkzO4sthSUA9O4Sw5Tru9IiIsjkZCKnp8Hl4/vvv+fKK6+s/Xz06NEA3H777bzxxhtkZWUxc+ZMjh07Rnx8PNdeey1TpkzRpRURkbN0vNLG3z7bxL9W5WMYEB0WwOSBXbkuNVZnO8StnNUNp42hITesiIh4i2+2HWLcnGx2HSkDYMi5LXm6fwpNQwNMTiZSw6VuOBURkTNnLa8i45M8PvhuFwBxkUE8OziVK5NbmJxM5MypfIiIuKhlGwt4al4OB6zlAPzp4kTG9kkmPEjL3ot7U/kQEXExh0sqeGZRLgt+2gdAm2YhTB+axsVt658pWsTdqHyIiLgIwzD4KGs/kxZu4EhpJT4W+MulbXmkV0eCA7QQnHgOlQ8RERdwoKic8fNz+HxjAQDJseFkDk2jW0ITc4OJNAKVDxHxHHYb7PwGSgogLAZap4OPa58xMAyD/67dzbRPNlJcXo2/r4UHruzAfVe0I8BPC8GJZ1L5EBHPkLsQFo8F675ft0XEQ59MSBloXq6T2HW4jHFzs/hm22EAuiU0YcbQNDq1CIGdq9yqRIk0hMqHiLi/3IUwawTwu2mLrPtrtt/4rksVEJvd4N/f7OBvSzZxvMpGkL8Pj13biTt6JuGb9xG8714lSqShVD5ExL3ZbTVnPH5fPODnbRZYPA6S+7nE2YMtBcWMmZPFj7uOAXBx2yimD0mjTXSo25UokTOl8iEi7m3nN3XPEpzAAOvemuOSLnVarN+rstl5c8U2XvliK5U2O2GBfjx5XWduviChZtl7NytRImdD5UNE3FtJgWOPawTZe4p4fPZ68g4UA3BVcgumDe5KXGTwrwe5SYkScQSVDxFxb2Exjj3OgcqrbLz0+Rbe/mo7NrtB0xB/Jg3swsBu8ScuBOcGJUrEUVQ+RMS9tU6vuSHTup/6L1lYava3TndqrO/yjzBuThbbD5UCMKBbPBMHpBAd9gcrfLtwiRJxNJUPEXFvPr41I0FmjQAs1C0gP59d6DPdafdJlFRUM2NxHu+u3glAi/BApg1O5ZqUU5QGFy1RIo1BM9iIiPtLGVgzEiQiru72iHinjhD5cvNBer+4srZ43NQ9gaWjLz918YBfSxRQW5pqOb9EiTQmi2EY9VVs01itViIjIykqKiIiIsLsOCLiTkya4fRYWSVTFm1kzg97AEiICmb6kDR6to9u+IPVO1lay5rioWG24sIa8v6tyy4i4jl8fJ0+EuTT7P08vWADh0oqsFjgjvQkHuvdkZCAM/z1mjKwZjitm00TL9IQKh8iImegsLiciQs28GnOAQDatwgjc2ga57duevYPbkKJEnEmlQ8RkQYwDIM5P+xlyqJcio5X4edj4b4r2vHAVe0J9NPZCZHTofIhInKa9hwt48l5OazcfBCAri0jyByaRpf4SJOTibgXlQ8RkVOw2w3+8+1OMj/No7TSRoCfD4/06sjdlybh56tBgyINpfIhInIS2w+WMHZOFmt3HAXggjZNmT40jXbNw0xOJuK+VD5EROpRbbPz9lf5vPj5Ziqr7YQE+DKubzJ/uqh1zUJwInLGVD5ERH4nd5+VMXPWk7PXCsClHaLJGJJKq6YhJicT8QwqHyIiP6uotvHqF1t5Y8U2qu0GEUF+PN0/hWHntzpxITgROWMqHyIiwA+7jjJmdhZbC0sA6NMllmcGdaFFeJDJyUQ8j8qHiHi1sspq/rZkM+98k49hQHRYIFOu70Lf1LhTf7GInBGVDxHxWqu2HmLc3Cx2HzkOwJDzWjKhfwpNQgJMTibi2VQ+RMTrFB2vIuOTjXy4djcA8ZFBTBuSypWdWpicTMQ7qHyIiFdZmlvA+PnZFFgrABjRozVj+iQTFqhfhyLOor9tIuIVDpdUMOmjXD5aX7NUfVJ0KJlD07gwKcrkZCLeR+VDRDyaYRgsXL+PSQs3cLSsCh8L3H1ZWx7p1ZEgfy0EJ2IGlQ8R8Vj7i44zfl4Oy/IKAUiODWfGsDTSWjUxN5iIl1P5EBGPYxgGH3y3m4xPNlJcUY2/r4UHr+rAvZe3I8BPC8GJmE3lQ0Q8ys7DpYybk83q7YcBOCehCTOGpdExJtzkZCLyC5UPEfEINrvBO6vy+dtnmyivshPk78Nj13bijp5J+GohOBGXovIhIm5vc0Exj8/OYv3uYwD0aNuM6UNTad0s1NxgIlIvlQ8RcVuV1XbeWLGNV5dvocpmEB7ox5P9OnPzBQlaCE7Ehal8iIhbytpzjDGzs8g7UAxAr84tmDooldhILQQn4upUPkTErZRX2Xhx6Wbe/mo7dgOiQgOYNLALA9LidLZDxE2ofIiI2/h2+2HGzc0m/1ApAAO7xTNxQArNwgJNTiYiDaHyISIur7i8iszFefxnzS4AYiICmTYolV4pMSYnE5EzofIhIi5t+aZCnpqbzb6icgBuuTCBJ67rTESQv8nJRORMqXyIiEs6WlrJlEW5zP1xLwCJUSFMH5JKevtok5OJyNlS+RARl2IYBp/mHGDCghwOlVRiscCdPZN49NqOhAToV5aIJ9DfZBFxGYXWcp5ekMOSDQUAdGgRRuawNM5LbGpyMhFxJJUPETGdYRj8b90epi7KxVpejZ+PhfuvaMfIq9oT6Kdl70U8jcqHiJhq95EynpyXzVdbDgGQ2jKSzKFppMRHmJxMRBqLyoeImMJuN3h39Q5mLNlEWaWNQD8fHrmmI3+5JAk/Xy17L+LJVD5ExOm2HSxh7Owsvt95FIAL20QxfWgqbZuHmZxMRJxB5UNEnKbKZuetldv5+7ItVFbbCQ3wZVzfZIZf1BofLXsv4jVUPkTEKTbsK2LM7Cw27LMCcHnH5kwb3JVWTUNMTiYizqbyISKNqrzKxqtfbOXNL7dRbTeIDPZnQv8UhpzXUgvBiXgplQ8RaTTrdh5hzOwsth2sWQjuutRYJg/sSvNwLQQn4s1UPkTE4UorqnluySZmrt6BYUB0WCBTB3WhT9c4s6OJiAtQ+RARh/p6yyHGzc1iz9HjAAw7vxXj+3WmSUiAyclExFWofIiIQxQdr2Lax7nM+n4PAC2bBPPskFQu79jc5GQi4mpUPkTkrH224QDj5+dQWFwBwO09WvN4n2TCAvUrRkROpN8MInLGDpVUMHHhBj7O2g9A2+hQMoelcUGbKJOTiYgrU/kQkQYzDIMFP+1j8kcbOFpWha+PhXsua8vDV3cgyF8LwYnIyal8iEiD7Dt2nPHzc/girxCAznERzBiaRmqrSJOTiYi7UPkQkdNitxt8sHYXGZ/kUVJRTYCvDw9d3Z6/Xt4Ofy0EJyIN0ODfGCtXrmTAgAHEx8djsViYP39+nf2GYTBhwgTi4uIIDg6mV69ebNmyxVF5RcQEOw6Vcuv/reGpeTmUVFRzbmITPn7oEh64qoOKh4g0WIN/a5SWltKtWzdee+21evfPmDGDl19+mTfffJNvv/2W0NBQevfuTXl5+VmHFRHnstkN3lq5jd4vrWTN9iME+/syoX8Ks+9Np0NMuNnxRMRNNfiyS9++fenbt2+9+wzD4KWXXmL8+PFcf/31ALz77rvExMQwf/58br755rNLKyJOs+lAMWNmr2f9niIAerZvRsbgNBKbaSE4ETk7Dr3nIz8/nwMHDtCrV6/abZGRkVx00UWsXr1a5UPEDVRW23l9xVZeW76VKptBeJAf4/t15sbuCVoITkQcwqHl48CBAwDExMTU2R4TE1O77/cqKiqoqKio/dxqtToykog0wPrdxxgzO4tNBcUA9Oocw7TBXYmJCDI5mYh4EtNHu2RkZDB58mSzY4h4teOVNl78fDP/99V27AY0Cw1g0sAu9E+L09kOEXE4h96mHhsbC0BBQUGd7QUFBbX7fu+JJ56gqKio9mP37t2OjCQip7Bm+2H6/n0lb62sKR7XnxPP0tGXM6BbvIqHiDQKh575SEpKIjY2lmXLlnHOOecANZdRvv32W+677756vyYwMJDAwEBHxhCR01BcXsX0T/N479tdAMRGBDFtcFeu7hxziq8UETk7DS4fJSUlbN26tfbz/Px8fvrpJ6KiokhMTGTUqFFMnTqVDh06kJSUxNNPP018fDyDBg1yZG4ROQvL8wp5cl42+4tqhsDfelEi4/omExHkb3IyEfEGDS4f33//PVdeeWXt56NHjwbg9ttv59///jdjxoyhtLSUe+65h2PHjnHJJZewePFigoJ0w5qI2Y6UVjJlUS7zftwLQOtmIWQMSSW9XbTJyUTEm1gMwzDMDvFbVquVyMhIioqKiIiIMDuOiEcwDIOPs/czccEGDpdW4mOBO3sm8ei1nQgO0EJwInL2GvL+bfpoFxFpXAXWcp6en8NnuTU3gneMCSNzaBrnJjY1OZmIeCuVDxEPZRgG//t+D1M+zqW4vBo/Hwsjr2zP/Ve2I9BPZztExDwqHyIeaPeRMp6cl81XWw4BkNYqkhnD0kiO1aVMETGfyoeIB7HZDd5dvYPnlmyirNJGoJ8Pj17bkTt7JuGn1WdFxEWofIh4iK2FxYydk826nUcBuDApisyhaSRFh5qcTESkLpUPETdXZbPz1srt/P3zLVTa7IQG+DLuus4MvzARHx/NUCoirkflQ8SN5ewtYszsLHL31yzIeEWn5kwbnErLJsEmJxMR+WMqHyJuqLzKxsvLtvCPldux2Q2ahPgzcUAKg85pqfVYRMTlqXyIuJnvdxxhzJwsth8sBaBfahyTBnahebjWSBIR96DyIeaw22DnN1BSAGEx0DodfDT3xMmUVlTz3JJNzFy9A8OA5uGBTLm+K3261r9itIiIq1L5EOfLXQiLx4J136/bIuKhTyakDDQvlwv7astBxs3JZu+x4wDccH4rxvdLITJEC8GJiPtR+RDnyl0Is0YAv1tSyLq/ZvuN76qA/EZRWRVTP87lf+v2ANCySTAZQ1K5rGNzk5OJiJw5lQ9xHrut5ozH74sH/LzNAovHQXI/XYIBlmw4wPj5ORwsrsBigdt7tOHx3p0IDdRfWxFxb/otJs6z85u6l1pOYIB1b81xSZc6LZarOVhcwaSFG/g4ez8AbZuHMmNoGt3bRJmcTETEMVQ+xHlKChx7nIcxDIN5P+7lmUW5HCurwtfHwl8va8tDV3cgyF9ngkTEc6h8iPOExTj2OA+y79hxnpyXzYpNBwFIiYtgxrA0uraMNDmZiIjjqXyI87ROrxnVYt1P/fd9WGr2t053djLT2O0G73+3i+mf5lFSUU2Arw8P9+rAPZe1xV8LwYmIh1L5EOfx8a0ZTjtrBGChbgH5eVbOPtO95mbT/EOljJuTxbf5RwA4v3VTMoem0b5FmMnJREQal8qHOFfKwJrhtPXO8zHdK4bZVtvs/GtVPs9/tpmKajvB/r6M6dOJET3a4GvWQnCa9E1EnEjlQ5wvZWDNcFovfLPLO2BlzOwssvYUAXBJ+2gyhqSSEBViXihN+iYiTqbyIebw8fWq4bSV1XZeXb6V15dvpdpuEB7kx9P9UriheytzF4LTpG8iYgKVD5FG9tPuY4yZvZ7NBSUAXJMSw9RBXYmJCDI3mCZ9ExGTqHyINJLjlTZeWLqJf36dj92AZqEBTL6+C/1S41xj2XtN+iYiJlH5EGkEq7cdZtzcLHYeLgNg8Lktebp/ClGhASYn+w1N+iYiJlH5EHEga3kVGZ/k8cF3uwCIiwxi2uCuXJXsghOnadI3ETGJyofI2fjNENUfjgTwwNdB7CuuAuDWixJ5om8y4UEuuuy9Jn0TEZOofIicqd8NUT0PmG1E8Xrk3fS78a/0aNfM3HynoknfRMQkmr9Z5EzkLsSYNQLjdzdsxlmOMKViBj0qVpkUrIF+mfQtIq7u9oh4DbMVkUajMx8iDWW3YftkDD4Y/H7MSu3n7jRE1YsnfRMRc6h8iDSAYRgs/2w+V5XsP9lR7jdE1csmfRMRc6l8iJym3UfKGDc3i2bbf+Cq0xkxqyGqIiL1UvkQOQWb3WDmNzt4bskmjlfZuNS/6el9oYaoiojUS+VD5CS2FhYzZnYWP+w6BsBFSVFMGXwP/OefGqIqInKGVD5E6lFls/OPL7fx8rKtVNrshAX68cR1ydxyQSI+PhYNURUROQsqHyK/k7O3iMdnZ7FxvxWAKzs1Z9rgVOKbBP960C9DVOtdin66hqiKiJyEyofIz8qrbPx92RbeWrkdm92gaYg/Ewd04fpz4utfCE5DVEVEzojKhwiwdscRxs7OYvuhUgD6pcUxeWAXosMCT/6FGqIqItJgKh/i1UoqqnlucR7vrtmJYUCL8ECmDOpK7y6xZkcTEfFYKh/itVZuPsgTc7PZe+w4ADd2b8VT16UQGeKiC8GJiHgIlQ/xOkVlVUz5OJfZ6/YA0KppMNOHpHFJh2iTk4mIeAeVD/Eqi3P28/SCDRwsrsBigT+nt+GxazsRGqi/CiIizqLfuOIVCovLmbRwA59kHwCgXfNQZgxL4/zWUSYnExHxPiof4tEMw2DuD3t5ZlEuRcer8PWxcO/lbXnwqg4E+WtIrIiIGVQ+xGPtPXacJ+dm8+XmgwCkxEUwY1gaXVtGmpxMRMS7qXyIx7HbDd77difTP82jtNJGgJ8PD1/dgXsua4u/r4/Z8UREvJ7Kh3iU/EOljJ2TxXf5RwDo3rop04em0b5FmMnJRETkFyof4hGqbXb+7+t8Xly6mYpqOyEBvozp3YkRPdrULAQnIiIuQ+VD3N7G/VbGzM4ie28RAJd2iObZwakkRIWYnExEROqj8iFuq6LaxmtfbOX1FduothtEBPnxdP8Uhp3fqv6F4ERExCWofIhb+nHXUcbMzmJLYQkAvbvEMOX6rrSICDI5mYiInIrKh7iVsspqnv9sM/9alY9hQHRYAJMHduW61Fid7RARcRMqH+I2vtl6iHFzs9l1pAyAIee25On+KTQNDTA5mYiINITKh7g8a3kVGZ9s5IPvdgMQFxnEs4NTuTK5hcnJRETkTKh8iEtbtrGAp+blcMBaDsCfLk5kbJ9kwoO07L2IiLtS+RCXdLikgskf5bJw/T4A2jQLIXNoGhe1bWZyMhEROVsqH+JSDMNg4fp9TP4olyOllfhY4O5L2/LINR21EJyIiIdQ+RCXcaConPHzs/l8YyEAybHhZA5No1tCE3ODiYiIQ6l8iOkMw+DDtbt59uONFFdU4+9r4YErO3DfFe0I8NNCcCIinkblQ0y163AZ4+Zm8c22wwB0S2jCc8PS6BgTbnIyERFpLCofYgqb3eDf3+zgb0s2cbzKRpC/D49d24k7eibhq4XgREQ8msqHON2WgmLGzMnix13HALi4bRTTh6TRJjrU3GAiIuIUKh/iNFU2O2+s2MarX2yl0mYnPNCPJ67rzM0XJGjZexERL6LyIU6RvaeIx2evJ+9AMQBXJ7dg6uCuxEUGm5xMRESczeFDCSZNmoTFYqnzkZyc7OinETdRXmVj+qd5DHp9FXkHimka4s/fbz6H/7u9u4qHiIiXapQzH126dOHzzz//9Un8dILFG32Xf4Sxc7LIP1QKwIBu8UwckEJ0WKDJyURExEyN0gr8/PyIjY1tjIcWN1BSUU3mp3n8vzU7AWgRHsi0walckxJjcjIREXEFjVI+tmzZQnx8PEFBQfTo0YOMjAwSExPrPbaiooKKioraz61Wa2NE8hx2G+z8BkoKICwGWqeDj+tMO/7l5oM8OTebvceOA3DzBQk8cV1nIoO1EJyIiNSwGIZhOPIBP/30U0pKSujUqRP79+9n8uTJ7N27l5ycHMLDT5w4atKkSUyePPmE7UVFRURERDgymvvLXQiLx4J136/bIuKhTyakDDQvF3CsrJJnFuUy94e9ACREBTN9SBo920ebmktERJzDarUSGRl5Wu/fDi8fv3fs2DFat27NCy+8wF133XXC/vrOfCQkJKh8/F7uQpg1Avj9H9fPQ1RvfNe0AvJJ9n4mLMjhUEklFgvckZ7EY707EhKge31ERLxFQ8pHo787NGnShI4dO7J169Z69wcGBhIYqBsQT8puqznjcULx4OdtFlg8DpL7OfUSTGFxORPmb2DxhgMAtG8RRubQNM5v3dRpGURExP00+qpdJSUlbNu2jbi4uMZ+Ks+185u6l1pOYIB1b81xTmAYBrPX7eGaF1ayeMMB/HwsPHhVez5+6BIVDxEROSWHn/l47LHHGDBgAK1bt2bfvn1MnDgRX19fbrnlFkc/lfcoKXDscWdhz9EynpyXw8rNBwHo2jKCGUO7kRKvS2QiInJ6HF4+9uzZwy233MLhw4dp3rw5l1xyCWvWrKF58+aOfirvEXaaQ1RP97gzYLcb/OfbnWR+mkdppY0APx8e6dWRuy9Nws9Xy96LiMjpc3j5+PDDDx39kNI6vWZUi3U/9d/3YanZ3zq9UZ5+28ESxs3JYu2OowBc0KYp04em0a55WKM8n4iIeDYNR3AHPr41w2lnjaBmdMtvC8jPo136THf4zabVNjtvf5XPi59vprLaTkiAL+P6JvOni1prITgRETljKh/uImVgzXDaeuf5mO7wYba5+6yMmbOenL01k75d1rE5zw7uSqumIQ59HhER8T4qH+4kZWDNcNpGnOG0otrGq19s5Y0V26i2G0QG+/N0/xSGntcSi0VnO0RE5OypfLgbH19IurRRHnrdzqOMnZPF1sISAPp0ieWZQV1oER7UKM8nIiLeSeVDKKus5rklm/j3NzswDIgOC2TK9V3om6q5WURExPFUPrzcqq2HGDc3i91HahaCG3peK57u35kmIQEmJxMREU+l8uGlio5XkfHJRj5cuxuAlk2CmTa4K1d0amFyMhER8XQqH15oaW4B4+dnU2CtWdBvRI/WjOmTTFigfhxERKTx6d3GixwuqWDiwg0sytoPQFJ0KJlD07gwKcrkZCIi4k1UPryAYRgsXL+PSQs3cLSsCl8fC3df2pZRvToQ5O+8VXBFRERA5cPj7S86zvh5OSzLKwQgOTac54Z1I7VVpMnJRETEW6l8eCi73eDDtbvJ+GQjxRXVBPj68OBV7fnr5e0I8NNCcCIiYh6VDw+083ApY+dksWb7EQDOSWjCc8PS6BATbnIyERERlQ+PYrMbvLMqn799tonyKjtB/j48dm0n7uiZhK8WghMREReh8uEhNhcU8/jsLNbvPgZAertmTB+SRmIzLQQnIiKuReXDzVVW23ljxTZeXb6FKptBeKAfT/XrzE0XJGghOBERcUkqH25s/e5jjJ2TRd6BYgB6dW7B1EGpxEZqITgREXFdKh9uqLzKxotLN/P2V9uxGxAVGsCkgV0YkBansx0iIuLyVD7czLfbDzNubjb5h0oBuP6ceCb0T6FZWKDJyURERE6PyoebKC6vInNxHv9ZswuAmIhApg1KpVdKjMnJREREGkblww0s31TIU3Oz2VdUDsAtFybwxHWdiQjyNzmZiIhIw6l8uLCjpZVMWZTL3B/3ApAYFcL0Iamkt482OZmIiMiZU/lwQYZh8En2ASYuzOFQSSUWC9zZM4lHr+1ISID+yERExL3pnczFFFrLeXpBDks2FADQoUUYmcPSOC+xqcnJREREHEPlw0UYhsH/1u1h6qJcrOXV+PlYuP+Kdoy8qj2Bflr2XkREPIfKhwvYfaSMJ+dl89WWQwCktoxkxrA0OsdFmJxMRETE8VQ+TGS3G7y7egczlmyirNJGoJ8Pj1zTkb9ckoSfr5a9FxERz6TyYZJtB0sYOzuL73ceBeDCNlFMH5pK2+ZhJicTERFpXCofTlZls/PWyu38fdkWKqvthAb4Mq5vMsMvao2Plr0XEREvoPLhRBv2FTFmdhYb9lkBuLxjc54dkkrLJsEmJxMREXEelQ8nKK+y8eoXW3nzy21U2w0ig/2Z0D+FIee11EJwIiLidVQ+Gtm6nUcYMzuLbQdrFoK7LjWWyQO70jxcC8GJiIh3UvloJKUV1Ty3ZBMzV+/AMCA6LJCpg7rQp2uc2dFERERMpfLRCL7ecohxc7PYc/Q4AMPOb8X4fp1pEhJgcjIRERHzqXw4UNHxKqZ9nMus7/cA0LJJMM8OSeXyjs1NTiYiIuI6VD4c5LMNBxg/P4fC4goAbu/Rmsf7JBMWqJdYRETkt/TOeJYOlVQwaeEGFmXtB6BtdCiZw9K4oE2UyclERERck8rHGTIMgwU/7WPyRxs4WlaFr4+Fey5ry8NXdyDIXwvBiYiI/BGVjzOw79hxxs/P4Yu8QgA6x0UwY2gaqa0iTU4mIiLi+lQ+GsBuN/hg7S4yPsmjpKKaAF8fHrq6PX+9vB3+WghORETktKh8nKYdh0oZNzeLNduPAHBuYhNmDE2jQ0y4yclERETci8rHKdjsBv/6Op/nl26ivMpOsL8vj/fuxO3pbfDVQnAiIiINpvJxEpsOFDNm9nrW7ykCoGf7ZmQMTiOxWYjJyURERNyXykc9KqvtvL5iK68t30qVzSA80I+n+nXmpgsStBCciIjIWVL5+J31u48xZnYWmwqKAejVOYapg7oSGxlkcjIRERHPoPLxs+OVNl78fDP/99V27AY0Cw1g0sAu9E+L09kOERERB1L5ANZsP8y4OVnsOFwGwPXnxDNxQBeiQrUQnIiIiKN5dfkoLq9i+qd5vPftLgBiI4KYNrgrV3eOMTmZiIiI5/La8rE8r5An52Wzv6gcgFsuTOSJ65KJCPI3OZmIiIhn87rycbS0kmcW5TLvx70AtG4WQsaQVNLbRZucTERExDt4TfkwDIOPs/czccEGDpdW4mOBO3sm8ei1nQgO0EJwIiIizuI15WPJhgIeeP9HADrGhJE5NI1zE5uanEpERMT7eE356NW5Bee3bkrP9tGMvLIdgX462yEiImIGrykffr4+zPprD63HIiIiYjKvWgdexUNERMR8XlU+RERExHwqHyIiIuJUKh8iIiLiVCofIiIi4lQqHyIiIuJUKh8iIiLiVF4zzwd2G+z8BkoKICwGWqeDjyYaExERcbZGO/Px2muv0aZNG4KCgrjooov47rvvGuupTi13IbzUFWb2hzl31fz3pa4120VERMSpGqV8/Pe//2X06NFMnDiRH374gW7dutG7d28KCwsb4+lOLnchzBoB1n11t1v312xXAREREXGqRikfL7zwAnfffTd33HEHKSkpvPnmm4SEhPCvf/2rMZ7uj9ltsHgsYNSz8+dti8fVHCciIiJO4fDyUVlZybp16+jVq9evT+LjQ69evVi9evUJx1dUVGC1Wut8OMzOb04841GHAda9NceJiIiIUzi8fBw6dAibzUZMTEyd7TExMRw4cOCE4zMyMoiMjKz9SEhIcFyYkgLHHiciIiJnzfShtk888QRFRUW1H7t373bcg4fFnPqYhhwnIiIiZ83hQ22jo6Px9fWloKDu2YSCggJiY2NPOD4wMJDAwEBHx6jROh0i4mtuLq33vg9Lzf7W6Y3z/CIiInICh5/5CAgI4Pzzz2fZsmW12+x2O8uWLaNHjx6OfrqT8/GFPpk/f2L53c6fP+8zXfN9iIiIOFGjXHYZPXo0b7/9NjNnzmTjxo3cd999lJaWcscddzTG051cykC48V2IiKu7PSK+ZnvKQOdnEhER8WKNMsPpTTfdxMGDB5kwYQIHDhzgnHPOYfHixSfchOo0KQMhuZ9mOBUREXEBFsMw6rsZwjRWq5XIyEiKioqIiIgwO46IiIichoa8f5s+2kVERES8i8qHiIiIOJXKh4iIiDiVyoeIiIg4lcqHiIiIOJXKh4iIiDiVyoeIiIg4lcqHiIiIOJXKh4iIiDhVo0yvfjZ+mXDVarWanERERERO1y/v26czcbrLlY/i4mIAEhISTE4iIiIiDVVcXExkZORJj3G5tV3sdjv79u0jPDwci8Xi0Me2Wq0kJCSwe/durRvTiPQ6O4deZ+fQ6+w8eq2do7FeZ8MwKC4uJj4+Hh+fk9/V4XJnPnx8fGjVqlWjPkdERIR+sJ1Ar7Nz6HV2Dr3OzqPX2jka43U+1RmPX+iGUxEREXEqlQ8RERFxKq8qH4GBgUycOJHAwECzo3g0vc7OodfZOfQ6O49ea+dwhdfZ5W44FREREc/mVWc+RERExHwqHyIiIuJUKh8iIiLiVCofIiIi4lReUz5ee+012rRpQ1BQEBdddBHfffed2ZE8TkZGBhdccAHh4eG0aNGCQYMGsWnTJrNjebzp06djsVgYNWqU2VE8zt69e/nTn/5Es2bNCA4OJjU1le+//97sWB7FZrPx9NNPk5SURHBwMO3atWPKlCmntT6InNzKlSsZMGAA8fHxWCwW5s+fX2e/YRhMmDCBuLg4goOD6dWrF1u2bHFKNq8oH//9738ZPXo0EydO5IcffqBbt2707t2bwsJCs6N5lC+//JKRI0eyZs0ali5dSlVVFddeey2lpaVmR/NYa9eu5R//+AdpaWlmR/E4R48epWfPnvj7+/Ppp5+Sm5vL888/T9OmTc2O5lEyMzN54403ePXVV9m4cSOZmZnMmDGDV155xexobq+0tJRu3brx2muv1bt/xowZvPzyy7z55pt8++23hIaG0rt3b8rLyxs/nOEFLrzwQmPkyJG1n9tsNiM+Pt7IyMgwMZXnKywsNADjyy+/NDuKRyouLjY6dOhgLF261Lj88suNhx9+2OxIHmXs2LHGJZdcYnYMj9evXz/jzjvvrLNtyJAhxvDhw01K5JkAY968ebWf2+12IzY21njuuedqtx07dswIDAw0Pvjgg0bP4/FnPiorK1m3bh29evWq3ebj40OvXr1YvXq1ick8X1FREQBRUVEmJ/FMI0eOpF+/fnV+tsVxFi5cSPfu3bnhhhto0aIF5557Lm+//bbZsTxOeno6y5YtY/PmzQCsX7+er7/+mr59+5qczLPl5+dz4MCBOr8/IiMjueiii5zy3uhyC8s52qFDh7DZbMTExNTZHhMTQ15enkmpPJ/dbmfUqFH07NmTrl27mh3H43z44Yf88MMPrF271uwoHmv79u288cYbjB49mieffJK1a9fy0EMPERAQwO233252PI8xbtw4rFYrycnJ+Pr6YrPZmDZtGsOHDzc7mkc7cOAAQL3vjb/sa0weXz7EHCNHjiQnJ4evv/7a7CgeZ/fu3Tz88MMsXbqUoKAgs+N4LLvdTvfu3Xn22WcBOPfcc8nJyeHNN99U+XCgWbNm8d577/H+++/TpUsXfvrpJ0aNGkV8fLxeZw/m8ZddoqOj8fX1paCgoM72goICYmNjTUrl2R544AEWLVrE8uXLadWqldlxPM66desoLCzkvPPOw8/PDz8/P7788ktefvll/Pz8sNlsZkf0CHFxcaSkpNTZ1rlzZ3bt2mVSIs/0+OOPM27cOG6++WZSU1O57bbbeOSRR8jIyDA7mkf75f3PrPdGjy8fAQEBnH/++Sxbtqx2m91uZ9myZfTo0cPEZJ7HMAweeOAB5s2bxxdffEFSUpLZkTzS1VdfTXZ2Nj/99FPtR/fu3Rk+fDg//fQTvr6+Zkf0CD179jxhqPjmzZtp3bq1SYk8U1lZGT4+dd+KfH19sdvtJiXyDklJScTGxtZ5b7RarXz77bdOeW/0issuo0eP5vbbb6d79+5ceOGFvPTSS5SWlnLHHXeYHc2jjBw5kvfff58FCxYQHh5ee90wMjKS4OBgk9N5jvDw8BPuowkNDaVZs2a6v8aBHnnkEdLT03n22We58cYb+e6773jrrbd46623zI7mUQYMGMC0adNITEykS5cu/Pjjj7zwwgvceeedZkdzeyUlJWzdurX28/z8fH766SeioqJITExk1KhRTJ06lQ4dOpCUlMTTTz9NfHw8gwYNavxwjT6exkW88sorRmJiohEQEGBceOGFxpo1a8yO5HGAej/eeecds6N5PA21bRwfffSR0bVrVyMwMNBITk423nrrLbMjeRyr1Wo8/PDDRmJiohEUFGS0bdvWeOqpp4yKigqzo7m95cuX1/s7+fbbbzcMo2a47dNPP23ExMQYgYGBxtVXX21s2rTJKdkshqFp5ERERMR5PP6eDxEREXEtKh8iIiLiVCofIiIi4lQqHyIiIuJUKh8iIiLiVCofIiIi4lQqHyIiIuJUKh8iIiLiVCofIiIi4lQqHyIiIuJUKh8iIiLiVCofIiIi4lT/H2Yo/H/K1mqjAAAAAElFTkSuQmCC", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "plt.plot(x,y,label=\"vérité vraie\")\n", "plt.plot(x,y_bruit,'o',label=\"obs\")\n", "plt.legend()\n" ] }, { "cell_type": "code", "execution_count": 37, "id": "a86af8ed", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "2.155853967974699 1.079475435345719\n" ] } ], "source": [ "n=10\n", "xbar = np.mean(x)\n", "xvar = np.var(x)\n", "ybar = np.mean(y_bruit)\n", "cov = np.sum((x-xbar)*(y_bruit-ybar))/n\n", "a_star = cov/xvar\n", "b_star = ybar - a_star*xbar\n", "print(a_star,b_star)" ] }, { "cell_type": "code", "execution_count": 38, "id": "5cdbc3db", "metadata": {}, "outputs": [], "source": [ "def f(x,a_star,b_star):\n", " return a_star*x + b_star" ] }, { "cell_type": "code", "execution_count": 39, "id": "7b6371a3", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "[1.079475435345719, 3.4748687330953847, 5.87026203084505, 8.265655328594717, 10.661048626344382, 13.056441924094047, 15.451835221843714, 17.84722851959338, 20.242621817343043, 22.638015115092713]\n", "[ 1. 3.22222222 5.44444444 7.66666667 9.88888889 12.11111111\n", " 14.33333333 16.55555556 18.77777778 21. ]\n", "[-9.80348061e-03 4.71610318e+00 8.63473352e+00 7.54105339e+00\n", " 8.02789833e+00 1.38513245e+01 1.37912889e+01 1.64466151e+01\n", " 2.28816547e+01 2.27065846e+01]\n" ] } ], "source": [ "y_pred = [f(xi,a_star,b_star) for xi in x]\n", "print(y_pred)\n", "print(y)\n", "print(y_bruit)" ] }, { "cell_type": "code", "execution_count": 41, "id": "223fa0bc", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 41, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAh8AAAGdCAYAAACyzRGfAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAA9hAAAPYQGoP6dpAABR40lEQVR4nO3dd3hUZf7+8ffMpFcIJJMEkhB6DRIEBERBooAUabpFF9vqFrBFheiiwCoScEXXXcvqdxd1f5ZlKYKoqCDFgoJkJYFI6KGm0FJJmzm/PyLRQCiByUxmcr+uK5fOOc/M+WQCmZvnPMVkGIaBiIiIiJOYXV2AiIiINC0KHyIiIuJUCh8iIiLiVAofIiIi4lQKHyIiIuJUCh8iIiLiVAofIiIi4lQKHyIiIuJUXq4u4Ex2u53Dhw8THByMyWRydTkiIiJyEQzDoKioiOjoaMzm8/dtNLrwcfjwYWJiYlxdhoiIiFyCAwcO0Lp16/O2aXThIzg4GKguPiQkxMXViIiIyMUoLCwkJiam5nP8fBpd+Dh9qyUkJEThQ0RExM1czJAJDTgVERERp1L4EBEREadS+BARERGnanRjPi6GYRhUVVVhs9lcXYpcgLe3NxaLxdVliIhII+J24aOiooIjR45QWlrq6lLkIphMJlq3bk1QUJCrSxERkUbCrcKH3W5n7969WCwWoqOj8fHx0UJkjZhhGOTn53Pw4EE6dOigHhAREQHcLHxUVFRgt9uJiYkhICDA1eXIRQgPD2ffvn1UVlYqfIiICOCmA04vtGyrNB7qmRIRkTO5Vc+HiIiIS9htkP01FOdCkBXiBoBZvbmXSuGjiamoqOAvf/kL48aNo0uXLq4uR0Sk8ctcDiunQeHhn46FRMPwudB1jOvqcmO6f+HG1q5di8lk4uTJkxf9nIcffpiMjAw6d+7ccIWJiHiKzOWwcFLt4AFQeKT6eOZy19Tl5hQ+3NiAAQM4cuQIoaGhALzxxhs0a9bsnO0XLlzItm3bePPNNzUWQ0TkQuy26h4PjDpO/nhsZUp1O3dht8HeLyBjUfV/XVS7bru4qcrKSnx8fIiMjLzo59xyyy3ccsstDViViIgHyf767B6PWgwoPFTdLn6Q08q6ZJnLsa2cRlrFMfItFsJtNhJ9WmBxwe0jt+/5MAyD0ooql3wZRl1p+GyvvfYa0dHR2O32Wsdvuukm7rrrLgCWLVtGYmIifn5+tG3bllmzZlFVVVXT1mQy8corrzBmzBgCAwOZPXt2rdsua9eu5c4776SgoACTyYTJZGLmzJkAlJeX88gjj9CqVSsCAwPp168fa9eudcj7LyLisYpzHdvOlTKXs2rF7xgWCndFWZkW0ZK7oqwMC4VVK37n9NtHbt/zcarSRtcnP3HJtTP/PIwAnwu/hTfffDP33Xcfa9asYejQoQAcP36clStX8tFHH/HFF18wadIkXnzxRQYNGsTu3bu59957AZgxY0bN68ycOZPU1FReeOEFvLy82LNnT825AQMG8MILL/Dkk0+SlZUFULOq6JQpU8jMzOS9994jOjqapUuXMnz4cDIyMujQoYPD3g8REY8SZHVsO1ex21j1eQrJES3OuoGUZ7GQHNGC+Z+nkNR5pNNm8Lh9z4c7aN68OSNGjOCdd96pObZo0SJatmzJkCFDmDVrFikpKdx+++20bduW66+/nqeeeop//OMftV7n17/+NXfeeSdt27YlNja21jkfHx9CQ0MxmUxERkYSGRlJUFAQ+/fvZ8GCBfz3v/9l0KBBtGvXjkceeYSrr76aBQsWOOX7FxFxS3EDqme1cK4xciYIaVXdrhGz7fuSVH+jOnicMd7P+PHxXH8D274vnVaT2/d8+HtbyPzzMJdd+2Ldeuut3HPPPbz88sv4+vry9ttv88tf/hKz2cyWLVv46quvmD17dk17m81GWVkZpaWlNau5XnnllfWuMSMjA5vNRseOHWsdLy8vp0WLFvV+PRGRJsNsqZ5Ou3AS1QHk5/0GP36ID09t9Ot9pOVsItfr3B/3hslEjpcXaTmb6NP2WqfU5Pbhw2QyXdStD1cbPXo0hmHw4Ycf0qdPH7744guef/55AIqLi5k1axbjx48/63l+fn41/x8YGFjv6xYXF2OxWNi8efNZy5trszcRkQvoOgZueesc63ykusU6H/mWi7vJcbHtHKHxf2p7CD8/P8aPH8/bb7/Nrl276NSpE4mJiQAkJiaSlZVF+/btL+saPj4+2Gy1p0316tULm81GXl4egwa5wWhsEZHGpusY6DzSbVc4DY8ZANsvfJs9PMZ5t48UPpzo1ltvZdSoUWzbto3bbrut5viTTz7JqFGjiI2NZeLEiTW3YrZu3crTTz990a/fpk0biouLWb16NT179iQgIICOHTty6623MmnSJJ577jl69epFfn4+q1evJiEhgZEjRzbEtyoi4lnMFveYTluHxMg+WL1DyKsoqBnj8XMmw8Dq24zEyD5Oq0kDTp3ouuuuIywsjKysLH7961/XHB82bBgrVqzg008/pU+fPlx11VU8//zzxMXF1ev1BwwYwO9//3t+8YtfEB4ezrx58wBYsGABkyZN4uGHH6ZTp06MHTuWTZs2nTVoVUREPI/FbCFl4CwwmTCdsUSEyTDAZGLagJlYnNiTYzIudrEKJyksLCQ0NJSCggJCQkJqnSsrK2Pv3r3Ex8fXGgshjZd+ZiIijcOq7FWkbkwlt/SndUkiA6xM65tCUlzSZb/++T6/z6TbLiIiIk1AUlwSQ2KGkJaXRn5pPuEB4SRGJDq1x+M0hQ8REZEmwmK20MeJYzvORWM+RERExKkUPkRERMSpFD5ERETEqRQ+RERExKkUPkRERMSpFD6asIqKCp555hl++OEHV5ciIiJNiMKHCw0ePJgHH3zQZdd/+OGHycjIoHPnzi6rQUREmp6mu86H3ea2mwQ5wsKFC9m2bRsrV67EVMda/yIiIg2laYaPzOXn2B55rltsj+wIt9xyC7fccouryxARkSao6d12yVwOCyfVDh4AhUeqj2cub5DLlpSUMGnSJIKCgoiKiuK5556rdf7EiRNMmjSJ5s2bExAQwIgRI9i5c2fN+TfeeINmzZrxySef0KVLF4KCghg+fDhHjhypaXPHHXcwduxY/vKXvxAVFUWLFi2YPHkylZWVNW3Ky8t55JFHaNWqFYGBgfTr14+1a9fWquXLL79k0KBB+Pv7ExMTw/33309JSUnN+ZdffpkOHTrg5+eH1Wpl4sSJDn63RETEkzWt8GG3Vfd4UNdeej8eW5lS3c7BHn30UdatW8eyZcv49NNPWbt2LWlpaTXn77jjDr777juWL1/Ohg0bMAyDG2+8sVZwKC0t5S9/+Qv//ve/Wb9+Pfv37+eRRx6pdZ01a9awe/du1qxZw5tvvskbb7zBG2+8UXN+ypQpbNiwgffee4/09HRuvvlmhg8fXhN0du/ezfDhw5kwYQLp6en85z//4csvv2TKlCkAfPfdd9x///38+c9/Jisri5UrV3LNNdc4/P0SEREPZjQyBQUFBmAUFBScde7UqVNGZmamcerUqUt78T3rDWNGyIW/9qy/zO+itqKiIsPHx8dYuHBhzbFjx44Z/v7+xgMPPGDs2LHDAIyvvvqq5vzRo0cNf3//mucsWLDAAIxdu3bVtHnppZcMq9Va8/j222834uLijKqqqppjN998s/GLX/zCMAzDyM7ONiwWi3Ho0KFa9Q0dOtR47LHHDMMwjLvvvtu49957a53/4osvDLPZbJw6dcpYvHixERISYhQWFl7U937ZPzMREXEL5/v8PlPTGvNRnHvhNvVpd5F2795NRUUF/fr1qzkWFhZGp06dAPjhhx/w8vKqdb5FixZ06tSp1jTYgIAA2rVrV/M4KiqKvLy8Wtfq1q0bFoulVpuMjAwAMjIysNlsdOzYsdZzysvLadGiBQBbtmwhPT2dt99+u+a8YRjY7Xb27t3L9ddfT1xcHG3btmX48OEMHz6ccePGERAQcMnvj4iINC1NK3wEWR3bzsm8vb1rPTaZTBiGccE2drsdgOLiYiwWC5s3b64VUACCgoJq2vzud7/j/vvvP+v6sbGx+Pj4kJaWxtq1a/n000958sknmTlzJps2baJZs2aX+y2KiEgT0LTGfMQNqJ7VwrmmlpogpFV1Owdq164d3t7efPvttzXHTpw4wY4dOwDo0qULVVVVtc4fO3aMrKwsunbt6rA6evXqhc1mIy8vj/bt29f6ioyMBCAxMZHMzMyzzrdv3x4fHx8AvLy8SEpKYt68eaSnp7Nv3z4+//xzh9UpIiKerWmFD7OlejotcHYA+fHx8FSHr/cRFBTE3XffzaOPPsrnn3/O1q1bueOOOzCbq9/+Dh06cNNNN3HPPffw5ZdfsmXLFm677TZatWrFTTfd5LA6OnbsyK233sqkSZNYsmQJe/fuZePGjcyZM4cPP/wQgGnTpvH1118zZcoUvv/+e3bu3MmyZctqBpyuWLGCF198ke+//57s7Gzeeust7HZ7zS0kERFp3ApOVXKipMKlNTSt8AHV63jc8haERNU+HhJdfbyB1vl49tlnGTRoEKNHjyYpKYmrr76a3r1715xfsGABvXv3ZtSoUfTv3x/DMPjoo4/Ouo1yuRYsWMCkSZN4+OGH6dSpE2PHjmXTpk3ExsYCkJCQwLp169ixYweDBg2iV69ePPnkk0RHRwPQrFkzlixZwnXXXUeXLl149dVXeffdd+nWrZtD6xQREcf7LDOXG55fx5PLt7m0DpNx5qABFyssLCQ0NJSCggJCQkJqnSsrK2Pv3r3Ex8fj5+d3eRdq4iucOotDf2YiInJJjhWXM/ODTD7YUr3GVXzLQN7/40BCAxz3D9zzfX6fqWkNOP05swXiB7m6ChERkQZjGAbLtxxm1geZHC+pwGyCe65py0NJHfHzdt0/uJtu+BAREfFgRwpOMX3pVlZvr16SoXNkMPMmJpDQuplrC0PhQ0RExKMYhsG7Gw8w56MfKCqvwtti4r7rOvD7a9vh49U4hnoqfIiIiHiI7GMlpCzOYMOeYwBcEdOMeRMT6GgNdnFltSl8iIiIuDmb3WDBV3v5y6dZlFXa8fM288gNnbhzYDwW87nWtnIdhQ8REZELsNltpOWlkV+aT3hAOIkRiVgayQzJHblFTF2UzvcHTgIwoF0LUscnENui8W57ofAhIiJyHquyV5G6MZXc0p/2/bIGWEnpm0JSXJLL6qqosvPK2t38fc1OKm0Gwb5ePD6yC7/sE4PJ1Ph6O35O4UNEROQcVmWvInltMga1l8TKK80jeW0y8wfPd0kAST94kqmL0tmeUwTA0M4RPD2uO1Gh/k6v5VIofIiIiNTBZreRujH1rOABYGBgwsTcjXMZEjPEabdgyiptPP/ZDl7/Yg92A8ICfZgxuitjekY3+t6On2scc27EYQYPHsyDDz540e3feOMN7UYrIlKHtLy0WrdazmRgkFOaQ1pemlPq+XbPMUb89Qv+sb46eIzpGc1nD13DTVe0cqvgAfUMH3PmzKFPnz4EBwcTERHB2LFjycrKqtWmrKyMyZMn06JFC4KCgpgwYQK5uef+4bmKzW5jU84mPtrzEZtyNmGz21xdkoiINCL5pfkObXepisoqmf5+Br947Rv2Hi3BGuLL/026khd/1YsWQb4Neu2GUq/bLuvWrWPy5Mn06dOHqqoqHn/8cW644QYyMzMJDAwE4KGHHuLDDz/kv//9L6GhoUyZMoXx48fz1VdfNcg3cCka6+AhERFpPMIDwh3a7lKsycrjT0syOFxQBsCv+saQMqILof6O3XTU2erV87Fy5UruuOMOunXrRs+ePXnjjTfYv38/mzdvBqCgoIB//vOfzJ8/n+uuu47evXuzYMECvv76a7755psG+Qbq6/TgoTO70k4PHlqVvapBrjt48GDuu+8+HnzwQZo3b47VauX111+npKSEO++8k+DgYNq3b8/HH39c85x169bRt29ffH19iYqKIiUlhaqqqprzJSUlTJo0iaCgIKKionjuuefOum55eTmPPPIIrVq1IjAwkH79+rF27drz1rps2TISExPx8/Ojbdu2zJo1q9Z1RUSagsSIRKwBVkzUfUvDhInIgEgSIxIdfu0TJRUk/+d77lywicMFZcSGBfDOb/sxZ3yC2wcPuMwxHwUFBQCEhYUBsHnzZiorK0lK+qn3oHPnzsTGxrJhw4Y6X6O8vJzCwsJaXw3lQoOHAOZunNtgt2DefPNNWrZsycaNG7nvvvv4wx/+wM0338yAAQNIS0vjhhtu4De/+Q2lpaUcOnSIG2+8kT59+rBlyxZeeeUV/vnPf/L000/XvN6jjz7KunXrWLZsGZ9++ilr164lLa32vccpU6awYcMG3nvvPdLT07n55psZPnw4O3furLPGL774gkmTJvHAAw+QmZnJP/7xD9544w1mz57dIO+JiEhjZTFbSOmbAnBWADn9eFrfaQ4fbPpRxhGuf34dS/53CJMJ7r46npUPDmJA+5YOvY5LGZfIZrMZI0eONAYOHFhz7O233zZ8fHzOatunTx9j6tSpdb7OjBkzDOCsr4KCgrPanjp1ysjMzDROnTp1STVvPLLR6P5G9wt+bTyy8ZJe/3yuvfZa4+qrr655XFVVZQQGBhq/+c1vao4dOXLEAIwNGzYYjz/+uNGpUyfDbrfXnH/ppZeMoKAgw2azGUVFRYaPj4+xcOHCmvPHjh0z/P39jQceeMAwDMPIzs42LBaLcejQoVq1DB061HjssccMwzCMBQsWGKGhobXOPfPMM7Xa//vf/zaioqIu6fu+3J+ZiIirfbbvM2PowqG1PieSFiYZn+37zKHXyS04Zfzure+MuGkrjLhpK4yk59Yam7OPO/QaDamgoOCcn99nuuSptpMnT2br1q18+eWXlxV+HnvsMZKTk2seFxYWEhMTc1mveS6uHjyUkJBQ8/8Wi4UWLVrQo0ePmmNWqxWAvLw8fvjhB/r3719rBPPAgQMpLi7m4MGDnDhxgoqKCvr161dzPiwsjE6dOtU8zsjIwGaz0bFjx1p1lJeX06JFizpr3LJlC1999VWtng6bzUZZWRmlpaUEBDTeFfNERBpCUlwSQ2KGNNgKp4Zh8N/NB3l6RSaFZVV4mU38cXA7Jl/XHl+vxrGKqqNdUviYMmUKK1asYP369bRu3brmeGRkJBUVFZw8ebLW9M3c3FwiIyPrfC1fX198fZ0zWtfVg4e8vWvfpzOZTLWOnQ4adrvdIdcrLi7GYrGwefNmLJbaf4CDgoLO+ZxZs2Yxfvz4s875+fk5pC4REXdjMVvoE9nH4a974Hgpjy/N4IudRwHo0SqUuRMS6Bod4vBrNSb1Ch+GYXDfffexdOlS1q5dS3x8fK3zvXv3xtvbm9WrVzNhwgQAsrKy2L9/P/3793dc1Zfo9OChvNK8Osd9mDBhDbA2yOCh+urSpQuLFy/GMIyaUPLVV18RHBxM69atCQsLw9vbm2+//ZbY2FgATpw4wY4dO7j22msB6NWrFzabjby8PAYNGnRR101MTCQrK4v27ds3zDcmIiLY7Qb//iabuSu3U1phw9fLzEPXd+S3V8fjZfH8JbjqFT4mT57MO++8w7JlywgODiYnJweA0NBQ/P39CQ0N5e677yY5OZmwsDBCQkK477776N+/P1dddVWDfAP1cXrwUPLaZEyYagWQhhw8dCn++Mc/8sILL3DfffcxZcoUsrKymDFjBsnJyZjNZoKCgrj77rt59NFHadGiBREREfzpT3/CbP7pD23Hjh259dZbmTRpEs899xy9evUiPz+f1atXk5CQwMiRI8+67pNPPsmoUaOIjY1l4sSJmM1mtmzZwtatW2sNdhURkUuzO7+YaYvS+S77BAB924SROqEHbcPr7pH2RPUKH6+88gpQPW305xYsWMAdd9wBwPPPP4/ZbGbChAmUl5czbNgwXn75ZYcU6whJcUnMHzy/znU+pvWd1mjW+WjVqhUfffQRjz76KD179iQsLIy7776b6dOn17R59tlnKS4uZvTo0QQHB/Pwww/XzEA6bcGCBTz99NM8/PDDHDp0iJYtW3LVVVcxatSoOq87bNgwVqxYwZ///Gfmzp2Lt7c3nTt35re//W2Dfr8iIp6u0mbntfV7+OvqnVRU2Qn0sZAyojO39ovD3Ai3vW9IJsMwzr7/4EKFhYWEhoZSUFBASEjte15lZWXs3buX+Pj4yx5/0Ji3R/YkjvyZiYi4q22HC5i6KJ1th6uXk7imYzjPjOtO6+aeM4j/fJ/fZ2qyG8s11OAhERGR08oqbfz98128um43VXaDUH9vnhzVlfGJ7rcfiyM12fAhIiLSkDZnH2fqonR255cAcGOPSGaN6U54sHvux+JICh8iIiIOVFJexbOfZPHmhn0YBrQM8uXpsd0Y3j3K1aU1GgofIiIiDvLlzqOkLEnn4IlTAEzs3ZrpI7vQLMDHxZU1LgofIiIil6ngVCWzP8xk4XcHAWjVzJ9nxvfg2o4Nt+OtO3PL8NHIJujIeehnJSKe7tNtOUx/fyt5ReUA3N4/jkeHdybI1y0/Yp3Crd6Z00uRl5aW4u/v7+Jq5GJUVFQAnLW8u4iIuztaXM6M5dv4MP0IAG1bBjJ3YgJ92oS5uLLGz63Ch8VioVmzZuTl5QEQEBDQpKcqNXZ2u538/HwCAgLw8nKrP2oiIudkGAbLvj/MrA+2caK0EovZxL3XtOWBoR3w89Y/tC6G230inN6g7nQAkcbNbDYTGxurkCgiHuHwyVNMf38rn2+v/gzqEhXCvAkJ9Ggd6uLK3IvbhQ+TyURUVBQRERFUVla6uhy5AB8fn1r7zYiIuCO73eDdTfuZ89F2isur8LGYuX9oe353bTu8m8BGcI7mduHjNIvFonEEIiLS4PYdLSFlSTrf7DkOQK/YZsybkEAHa7CLK3Nfbhs+REREGpLNbvDPL/fw3Kc7KK+y4+9t4ZFhnbhjQBssTWwjOEdT+BARETlDVk4RUxdtYcvB6p3CB7ZvwZxxCcS28JyN4FxJ4UNERORHFVV2Xl67i5fW7KLSZhDs58X0kV245coYDZx3IIUPERERYMuBk0xdlE5WbhEASV2szB7XHWuIn4sr8zwKHyIi0qSdqrDx/Kod/N8Xe7Ab0CLQh5ljujEqIUq9HQ1E4UNERJqsb/YcI2VxOvuOlQJw0xXRzBjdjbBAbQTXkBQ+RESkySkqqyT14+28/e1+ACJD/Jg9rjtDu1hdXFnToPAhIiJNyprteTy+NIMjBWUA/LpfLCkjOhPi5+3iypoOhQ8REWkwNruNtLw08kvzCQ8IJzEiEYvZNQtEniip4M8rMln6v0MAxLUIYM74Hgxo19Il9TRlCh8iItIgVmWvInVjKrmluTXHrAFWUvqmkBSX5LQ6DMPgw4wjzFi2jWMlFZhNcNfAeB6+oRP+Plop2xUUPkRExOFWZa8ieW0yBkat43mleSSvTWb+4PlOCSB5hWVMf38rn2ZWB6CO1iDmTkigV2zzBr+2nJvCh4iIOJTNbiN1Y+pZwQPAwMCEibkb5zIkZkiD3YIxDIP/fneQpz7MpKisCi+ziclD2jN5SHt8vLQRnKspfIiIiEOl5aXVutVyJgODnNIc0vLS6BPZx+HXP3C8lMeXZvDFzqMAJLQOZd7EBDpHhjj8WnJpFD5ERMSh8kvzHdruYtntBm9u2Mezn2RRWmHD18vMwzd05K6B8Xhp2/tGReFDREQcKjwg3KHtLsauvGKmLU5nc/YJAPrGhzF3QgLxLQMddg1xHIUPERFxqMSIRKwBVvJK8+oc92HChDXASmJE4mVfq9Jm57X1e/jrqp1U2OwE+XqRMqIzv+4bi1nb3jda6ocSERGHspgtpPRNAaqDxs+dfjyt77TLHmy69VABN/39K579JIsKm53BncL59KFruO2qOAWPRk7hQ0REHC4pLon5g+cTERBR67g1wHrZ02zLKm3MW7mdm176iswjhTQL8Ob5X/RkwR19iG7mf7mlixPotouIiDSIpLgkhsQMcegKp5uzjzN1UTq780sAGNkjipljuhEe7OuossUJFD5ERKTBWMwWh0ynLSmv4tlPsnhzwz4MA8KDfXnqpu4M7x7pgCrF2RQ+RESkUftiZz6PLcng4IlTANzcuzXTR3YlNEAbwbkrhQ8REWmUCkormf1RJgu/OwhAq2b+pE7owaAOjpuiK66h8CEiIo3OJ9tymP7+VvKLyjGZ4Pb+bXh0WCcCffWx5Qn0UxQRkUYjv6icmcu38WHGEQDahgcyb0ICV7YJc3Fl4kgKHyIi4nKGYfD+94eY9UEmJ0srsZhN/P7attx3XQf8vLXtvadR+BAREZc6fPIUf1qawZqs6r1eukaFMG9iAt1bhbq4MmkoCh8iIuISdrvBOxv3k/rxdorLq/CxmHkgqQP3XtMWb20E59EUPkRExOn2Hi0hZXE63+49DkDvuObMnZBA+4ggF1cmzqDwISIiTlNls/Ovr/by3Kc7KK+yE+BjYeqwTvymfxss2o+lyVD4EBGRhmO3QfbXUJzL/opgHvjaj/8dKgbg6vYtmTO+BzFhAS4uUpxN4UNERBpG5nJYOQ0KDwMQC7xkhDHP704GjLyTm69sjcmk3o6mSOFDREQcL3M5LJyEgcHP40WU6TjPMx9T4BVginFVdeJiGk4sIiKOZbdh/3jaWcEDwPTjFytTqm/JSJOk8CEiIg61bcNKzEWHzwoePzGg8FD1WBBpknTbRUREHKKorJI5H2+neNPXvOhzEU8ozm3wmqRxUvgQEZHLtmZ7Ho8vzeBIQRlXmZtd3JOCrA1akzReuu0iIiKX7HhJBQ++9z/ufGMTRwrKiGsRwAN33g4h0XDOGy8mCGkFcQOcWao0Iur5EBGRejMMgxXpR5i5fBvHSiowm+C3g9ryUFJH/H0sMHwuLJxEdQAxfvbMHwPJ8FQwa8O4pkrhQ0RE6iW3sIzp72/ls8zqMRudrMHMm5hAz5hmPzXqOgZueavWOh9AdY/I8NTq89JkKXyIiMhFMQyDhd8d4OkPf6CorApvi4nJQ9rzx8Ht8fGq4y5+1zHQeWTNCqcEWatvtajHo8lT+BARkQs6cLyUx5Zk8OWuowD0bB3KvIk96RQZfP4nmi0QP8gJFYo7UfgQEZFzstkN3vx6H89+ksWpShu+XmYeuaETd10dr43g5JIpfIiISJ125RUxdVE6aftPAtAvPoy5ExJo0zLQtYWJ21P4EBGRWiptdv6xbjcvrt5Fhc1OkK8Xj93YmV/1icWs3g5xAIUPEZFGxma3kZaXRn5pPuEB4SRGJGJx0iDNrYcKmLooncwjhQAM6RTO7HE9iG7m75TrS9Og8CEi0oisyl5F6sZUckt/WnrcGmAlpW8KSXFJDXbdskobf129k9fW78FmN2ge4M2M0d246YpobXsvDqfwISLSSKzKXkXy2mSMWotyQV5pHslrk5k/eH6DBJBN+44zbVE6e46WADAyIYpZY7rRMsjX4dcSAYUPEZFGwWa3kbox9azgAfy4Nb2JuRvnMiRmiMNuwZSUVzFv5Xbe+iYbw4CIYF+eGtudYd0iHfL6Iuei8CEi0gik5aXVutVyJgODnNIc0vLS6BPZ57Kvt35HPo8tyeDQyVMA/OLKGB6/sQuhAd6X/doiF6LwISLSCOSX5ju03bkUlFby1IeZLNp8EIDWzf1JHZ/A1R1aXtbritSHwoeISCMQHhDu0HZ1Wbn1CE8s20Z+UTkmE9wxoA2P3NCJQF99FIhz1bEY//mtX7+e0aNHEx1dPQL6/fffr3X+jjvuwGQy1foaPny4o+oVEfFIiRGJWAOsmM6xDb0JE5EBkSRGJNb7tfOLyvnj25v5/f9LI7+onHbhgSz6fX9mjO6m4CEuUe/wUVJSQs+ePXnppZfO2Wb48OEcOXKk5uvdd9+9rCJFRDydxWwhpW8KwFkB5PTjaX2n1WuwqWEYLEk7yPXPr+OjjBwsZhOTh7Tjw/sH0TsuzHHFi9RTvSPviBEjGDFixHnb+Pr6Ehmp0dIiIvWRFJfE/MHz61znY1rfafWaZnvo5Cn+tDSDtVnVY0S6RoUwb2IC3VuFOrxukfpqkP62tWvXEhERQfPmzbnuuut4+umnadGiRZ1ty8vLKS8vr3lcWFjYECWJiLiFpLgkhsQMueQVTu12g7c37if1ox8oqbDh42XmgaEduPeatnhb6t3ZLdIgHB4+hg8fzvjx44mPj2f37t08/vjjjBgxgg0bNmCxnP2XZ86cOcyaNcvRZYiIuC2L2XJJ02n3Hi1h2uJ0Nu49DsCVcc1JnZBA+4ggR5cocllMhmGcvaLNxT7ZZGLp0qWMHTv2nG327NlDu3btWLVqFUOHDj3rfF09HzExMRQUFBASEnKppYmINBlVNjv//HIv8z/bQXmVnQAfC1OHdWJS/zbaCE6cprCwkNDQ0Iv6/G7wYc5t27alZcuW7Nq1q87w4evri6+vlvAVEbkUPxwpZNridNIPFgAwqENLnhnXg5iwABdXJnJuDR4+Dh48yLFjx4iKimroS4mINBnlVTZeWrObl9fsospuEOLnxfRRXbm5d2ttBCeNXr3DR3FxMbt27ap5vHfvXr7//nvCwsIICwtj1qxZTJgwgcjISHbv3s3UqVNp3749w4YNc2jhIiJN1f/2n2DqonR25hUDMKybladu6k5EiJ+LKxO5OPUOH9999x1DhgypeZycnAzA7bffziuvvEJ6ejpvvvkmJ0+eJDo6mhtuuIGnnnpKt1ZERC7TqQobf/k0i399tRfDgJZBPswa050be0Sqt0PcymUNOG0I9RmwIiLSVHy9+ygpizPYf7wUgPG9WvHEqK40D/RxcWUi1RrVgFMREbl0hWWVzPloO+9u3A9AVKgfz4zrwZDOES6uTOTSKXyIiDRSq3/I5U9Lt5JTWAbAbVfFMm14Z4L9tO29uDeFDxGRRuZYcTl/XpHJsu8PA9CmRQCpExK4qm3dK0WLuBuFDxGRRsIwDD5IP8LM5ds4XlKB2QS/HdSWh5I64u9z8RvKiTR2Ch8iIo1ATkEZ09/fyqofqjeU6xwZzNwJCfSMaebawkQagMKHiHgOuw2yv4biXAiyQtwAqMcW9K5gGAb/2XSA2R/9QFFZFd4WE1OGdOAPg9vh46WN4MQzKXyIiGfIXA4rp0Hh4Z+OhUTD8LnQdYzr6jqP/cdKSVmSzte7jwHQM6YZ8yYk0CkiALK/cqsQJVIfCh8i4v4yl8PCScAZyxYVHqk+fstbjSqA2OwGb3y9j798ksWpSht+3mYeuaETdw6Mx7L9A3jHvUKUSH0pfIiIe7Pbqns8zgwe8OMxE6xMgc4jG0Xvwc7cIqYuTud/+08CcFXbMFLHJ9CmZaDbhSiRS6XwISLuLfvr2r0EZzGg8FB1u/hBTivrTJU2O6+u3c3fPt9Fhc1OkK8Xj9/YhV/2iane9t7NQpTI5VD4EBH3Vpzr2HYNIONgAY8u2sL2nCIAruscwexx3YkK9f+pkZuEKBFHUPgQEfcWZHVsOwcqq7TxwqqdvP7FHmx2g+YB3swc040xPaPP3gjODUKUiKMofIiIe4sbUD0gs/AIdd+yMFWfjxvg1LI27j1OyuJ09hwtAWB0z2hmjO5Ky6Bz7PDdiEOUiKMpfIiIezNbqmeCLJwEmKgdQH7sXRie6rRxEsXlVcxbuZ23NmQDEBHsy+xxPbi+6wVCQyMNUSINQSvYiIj76zqmeiZISFTt4yHRTp0hsm5HPsOeX18TPH5xZQyfJV974eABP4UooCY01XB+iBJpSCbDMOqK2C5TWFhIaGgoBQUFhISEuLocEXEnLlrh9GRpBU+t+IHFaQcBiAnzJ3V8AgPbt6z/i9W5WFqr6uChabbSiNXn81u3XUTEc5gtTp8J8nHGEZ5Yto2jxeWYTHDngHgeGdaRAJ9L/PXadUz1dFo3WyZepD4UPkTEY9jsNtLy0sgvzSc8IJzEiEQsDfShnVdUxoxl2/h4aw4A7SOCmDshgd5xzS//xV0QokScSeFDRDzCquxVpG5MJbf0p6mo1gArKX1TSIpLcth1DMNgcdohnlqRScGpSrzMJv4wuB1TrmuPr5d6J0QuhsKHiLi9VdmrSF6bjHHGLJG80jyS1yYzf/B8hwSQgydKeXzpVtbvyAege6sQ5k5IoFt06GW/tkhTovAhIm7NZreRujH1rOABYGBgwsTcjXMZEjPkkm/B2O0G/+/bbOZ+vJ2SChs+XmYeSurIPYPi8bJo0qBIfSl8iIhbS8tLq3Wr5UwGBjmlOaTlpdEnsk+9X39PfjHTFqezad8JAPq0aU7qhATahQddcs0iTZ3Ch4i4tfzSfIe2O63KZuf1L/by/KodVFTZCfCxkDKiM7f1i6veCE5ELpnCh4i4tfCAcIe2A8g8XMjUxVvYeqgQgEEdWjJnfA9aNw+4pBpFpDaFDxFxa4kRiVgDrOSV5tU57sOECWuAlcSIxAu+VnmVjb9/votX1u6mym4Q4ufFE6O6MrF367M3ghORS6aRUiLi1ixmCyl9U4DqoPFzpx9P6zvtgoNN0/afYOSLX/K3z3dRZTcY3i2SVQ9fy81Xxih4iDiYwoeIuL2kuCTmD55PREBErePWAOsFp9mWVlTx5w8ymfDK1+zKK6ZlkC+v3JrIq7/pTUSwX0OXLtIk6baLiHiEpLgkhsQMqdcKp1/tOkrKknQOHD8FwPjEVjw5qivNAnycVbZIk6TwISIew2K2XNR02oJTlcz56Afe23QAgOhQP2aP78GQThEXeKaIOILCh4g0KZ9l5jL9/QxyC8sBmNQ/jqnDOxPkq1+HIs6iv20i0iQcKy5n5geZfLCleqv6+JaBzJ2QQN/4MBdXJtL0KHyIiEczDIPlWw4zc/k2TpRWYjbBPde05aGkjvh5ayM4EVdQ+BARj3Wk4BTTl25l9fY8ADpHBjNvYgIJrZu5tjCRJk7hQ0Q8jmEYvLvxAHM++oGi8iq8LSbuu64Dv7+2HT5eWmFAxNUUPkTEo2QfKyFlcQYb9hwD4IqYZsybmEBHa7CLKxOR0xQ+RMQj2OwGC77ay18+zaKs0o6ft5lHbujEnQPjsWgjOJFGReFDRNzejtwiHl2UzpYDJwHo37YFqRN6ENci0LWFiUidFD5ExG1VVNl5Ze1u/r5mJ5U2g2BfLx4f2YVf9tF+LCKNmcKHiLil9IMnmboone05RQAkdYng6bE9iAzVfiwijZ3Ch4i4lbJKG89/toPXv9iD3YCwQB9mjunG6IQo9XaIuAmFDxFxG9/uOUbKkgz2Hi0BYEzPaGaM7kqLIF8XVyYi9aHwISKNXlFZJXNXbuf/fbMfAGuIL7PH9iCpq9XFlYnIpVD4EJFGbU1WHn9aksHhgjIAftU3hsdu7EKIn7eLKxORS6XwISKN0omSCp5akcmS/x0CIDYsgNTxPRjQvqWLKxORy6XwISKNimEYfLw1hyeXbeVocQUmE9w1MJ6Hb+hIgI9+ZYl4Av1NFpFGI6+wjCeWbeWTbbkAdIgIYu7EBBJjm7u4MhFxJIUPEXE5wzD47+aDPL0ik8KyKrzMJv44uB2Tr2uPr5e2vRfxNAofIuJSB46X8vjSDL7YeRSAHq1CmTshga7RIS6uTEQaisKHiLiE3W7w1oZ9zPski9IKG75eZh66viO/vToeL4u2vRfxZAofIuJ0u/OLmbYone+yTwDQt00YqRN60DY8yMWViYgzKHyIiNNU2uy8tn4Pf129k4oqO4E+FlJGdObWfnGYte29SJOh8CEiTrHtcAFTF6Wz7XAhANd2DGf2uO60bh7g4spExNkUPkQEm91GWl4a+aX5hAeEkxiRiMXsmFkmZZU2/v75Ll5dt5squ0GovzdPjurK+MRW2ghOpIlS+BBp4lZlryJ1Yyq5pbk1x6wBVlL6ppAUl3RZr705+zhTF6WzO796I7gbe0Qya0x3woO1EZxIU6bwIdKErcpeRfLaZAyMWsfzSvNIXpvM/MHzLymAlJRX8ewnWby5YR+GAS2DfHl6bDeGd49yVOki4sYUPkSaKJvdRurG1LOCB4CBgQkTczfOZUjMkHrdgvly51FSlqRz8MQpACb2bs30kV1oFuDjsNpFxL0pfIg0UWl5abVutZzJwCCnNIe0vDT6RPa54OsVnKpk9oeZLPzuIACtmvnzzPgeXNsx3GE1i4hnUPgQaaLyS/Md1u7TbTlMf38reUXlANzeP45Hh3cmyFe/YkTkbPrNINJEhQdcXI/E+dodLS5nxvJtfJh+BIC2LQOZOzGBPm3CHFKjiHgmhQ+RJioxIhFrgJW80rw6x32YMGENsJIYkXjWOcMwWPb9YWZ9sI0TpZVYzCbuvaYtDwztgJ+3NoITkfNT+BBpoixmCyl9U0hem4wJU60AYqJ6/Y1pfaedNdj08MlTTH9/K59vzwOgS1QI8yYk0KN1qPOKFxG3pt2bRJqwpLgk5g+eT0RARK3j1gDrWdNs7XaDt7/N5obn1/P59jx8LGYeuaEjy6cMVPAQkXqpd8/H+vXrefbZZ9m8eTNHjhxh6dKljB07tua8YRjMmDGD119/nZMnTzJw4EBeeeUVOnTo4Mi6RcRBkuKSGBIz5LwrnO47WkLKknS+2XMcgF6xzZg3IYEO1mBXlS0ibqzePR8lJSX07NmTl156qc7z8+bN48UXX+TVV1/l22+/JTAwkGHDhlFWVnbZxYpIw7CYLfSJ7MONbW+kT2SfmuBhsxu8tn43w15Yzzd7juPvbeHJUV1Z9PsBCh4icsnq3fMxYsQIRowYUec5wzB44YUXmD59OjfddBMAb731Flarlffff59f/vKXl1etiDhNVk4RUxdtYcvBAgAGtm/BnHEJxLbQRnAicnkcOuB079695OTkkJT0033i0NBQ+vXrx4YNGxQ+RNxARZWdl9fu4qU1u6i0GQT7eTF9ZBduuTJGG8GJiEM4NHzk5OQAYLVaax23Wq01585UXl5OeXl5zePCwkJHliQi9bDlwEmmLkonK7cIgKQuVmaP6441xM/FlYmIJ3H5VNs5c+Ywa9YsV5ch0qSdqrDx/Kod/N8Xe7Ab0CLQh5ljujEqIUq9HSLicA6dahsZGQlAbm7t/SJyc3Nrzp3pscceo6CgoObrwIEDjixJRC7gmz3HGPHX9by2vjp43HRFNJ8lX8vontEKHiLSIBza8xEfH09kZCSrV6/miiuuAKpvo3z77bf84Q9/qPM5vr6++Pr6OrIMEbkIRWWVpH68nbe/3Q9AZIgfs8d1Z2gX6wWeKSJyeeodPoqLi9m1a1fN47179/L9998TFhZGbGwsDz74IE8//TQdOnQgPj6eJ554gujo6FprgYiIa63ZnsfjSzM4UlA9Bf7X/WJJGdGZED9vF1cmIk1BvcPHd999x5AhQ2oeJycnA3D77bfzxhtvMHXqVEpKSrj33ns5efIkV199NStXrsTPTwPWRFzteEkFT63IZOn/DgEQ1yKAOeN7MKBdSxdXJiJNickwjLN3lHKhwsJCQkNDKSgoICQkxNXliHgEwzD4MOMIM5Zt41hJBWYT3DUwnodv6IS/jzaCE5HLV5/Pb5fPdhGRhpVbWMYT72/l08zqgeAdrUHMnZBAr9jmLq5MRJoqhQ8RD2UYBv/97iBPfZhJUVkVXmYTk4e0549D2uHrpd4OEXEdhQ8RD3TgeCmPL83gi51HAUhoHcq8iQl0jtStTBFxPYUPEQ9isxu8tWEfz36SRWmFDV8vMw/f0JG7BsbjZXHosj4iIpdM4UPEQ+zKK2La4gw2Z58AoG98GHMnJBDfMtDFlYmI1KbwIeLmKm12Xlu/h7+u2kmFzU6gj4WUG7twa99YzGatUCoijY/Ch4gb23qogKmL0sk8Ur0h4+BO4cwe14NWzfxdXJmIyLkpfIi4obJKGy+u3sk/1u/BZjdoFuDNjNFdGXtFK+3HIiKNnsKHiJv5bt9xpi5OZ09+CQAje0Qxc0w3woO1R5KIuAeFD3ENuw2yv4biXAiyQtwAMGvtifMpKa/i2U+yeHPDPgwDwoN9eeqm7gzvXveO0SIijZXChzhf5nJYOQ0KD/90LCQahs+FrmNcV1cj9sXOfFIWZ3Do5CkAbu7dmukjuxIaoI3gRMT9KHyIc2Uuh4WTgDO2FCo8Un38lrcUQH6moLSSpz/M5L+bDwLQqpk/c8b34JqO4S6uTETk0il8iPPYbdU9HmcGD/jxmAlWpkDnkboFA3yyLYfp728lv6gckwlu79+GR4d1ItBXf21FxL3pt5g4T/bXtW+1nMWAwkPV7eIHOa2sxia/qJyZy7fxYcYRANqGBzJvQgJXtglzcWUiIo6h8CHOU5zr2HaNgM1uIy0vjfzSfMIDwkmMSMRyib02hmGw9H+H+POKTE6WVmIxm/jdNW25f2gH/LzVEyQinkPhQ5wnyOrYdi62KnsVqRtTyS39KSxZA6yk9E0hKS6pXq91+OQpHl+awdqsfAC6RoUwb2IC3VuFOrRmEZHGQDtNifPEDaie1cK5FsEyQUir6naN3KrsVSSvTa4VPADySvNIXpvMquxVF/U6drvB//smmxueX8/arHx8LGYeHdaJZVMGKniIiMdS+BDnMVuqp9MCZweQHx8PT230g01tdhupG1Mx6hg4e/rY3I1zsdlt532dvUdL+NXr3zD9/a0Ul1fRO645Hz0wiMlD2uOtHWhFxIPptos4V9cx1dNp61znI9Utptmm5aWd1ePxcwYGOaU5pOWl0Seyz1nnq2x2/vXVXp77dAflVXb8vS1MHd6JSf3bYHHVRnBa9E1EnEjhQ5yv65jq6bRu+mGXX5p/ye225xQydVE66QcLALi6fUvmjO9BTFiAQ2usFy36JiJOpvAhrmG2uO102vCAi1vg6+ftKqrs/H3NLl5es4squ0GwnxdPjOzKzVe2du1GcFr0TURcQOFDpJ4SIxKxBljJK82rc9yHCRPWACuJEYkAfH/gJFMXbWFHbjEA13e18vTY7lhD/Jxa91m06JuIuIhGtYnUk8VsIaVvClAdNH7u9ONpfadRUQWzP8xk/MtfsSO3mBaBPvz917147Te9XR88oH6LvomIOJDCh8glSIpLYv7g+UQERNQ6bg2wMn/wfAKrejH8r+t5/Yu92A0Y16sVnyVfy6iEaNfeZvk5D1z0TUTcg267iFyipLgkhsQMqbXCafuQHsxbuZN3N34DQFSoH7PHdee6zo1w4TQPW/RNRNyHwofIZbAAfU6VQUkpaQeOMPLLExwuqgTg1/1ieWxEZ4L9Gum296cXfSs8Qt3jPkzV591g0TcRcS8KHyKX6owpqonAIiOMl0PvYeQtv6N/uxaure9CTi/6tnAS1Yu8/TyAuM+ibyLifjTmQ+RSZC7HWDgJ44wBm1Gm4zxVPo/+5V+5qLB6Or3oW0hU7eMh0ZpmKyINRj0fIvVlt2H7aCpmjHMtEu9eU1TdfNE3EXE/Ch8i9WAYBms+fZ/rio+cr9VPU1TdZSE1N170TUTcj8KHyEU6cLyUlCXptNiTxnU+F/EETVEVEamTwofIBdjsBm9+vY9nP8niVKWNQd7NL+6JmqIqIlInhQ+R89iVV8TURemk7T8JQL/4MJ4ady/8v39qiqqIyCVS+BCpQ6XNzj/W7ebF1buosNkJ8vXisRs786s+sZjNJk1RFRG5DAofImfYeqiARxel88ORQgCGdApn9rgeRDfz/6nR6SmqdW5Fn6opqiIi56HwIfKjskobf129k9fW78FmN2ge4M2M0d246Ypz7MeiKaoiIpdE4UME2LTvONMWpbPnaAkAIxOimDWmGy2DfM//RE1RFRGpN4UPadKKy6t4duV23vomG8OAiGBfnhrbnWHdIl1dmoiIx1L4kCZr/Y58HluSwaGTpwC45crW/OnGroQGNNKN4EREPITChzQ5BaWVPPVhJos2HwSgdXN/UscncHWHli6uTESkaVD4kCZl5dYjPLFsG/lF5ZhMcMeANjxyQycCffVXQUTEWfQbV5qEvKIyZi7fxkcZOQC0Cw9k3sQEeseFubgyEZGmR+FDPJphGCxJO8SfV2RScKoSi9nE769ty33XdcDPW1NiRURcQeFDPNahk6d4fEkG63bkA9A1KoR5ExPo3irUxZWJiDRtCh/icex2g7e/zSb14+2UVNjw8TLzwNAO3HtNW7wtZleXJyLS5Cl8iEfZe7SEaYvT2bj3OABXxjUndUIC7SOCXFyZiIicpvAhHqHKZuf/vtzL85/toLzKToCPhanDOjGpf5vqjeBERKTRUPgQt/fDkUKmLkon41ABAIM6tOSZcT2ICQtwcWUiIlIXhQ9xW+VVNl76fBcvr91Nld0gxM+LJ0Z1ZWLv1nVvBCciIo2Cwoe4hM1uIy0vjfzSfMIDwkmMSMRSj91g/7f/BFMXpbMzrxiAYd2sPHVTdyJC/BqqZBERcRCFD3G6VdmrSN2YSm5pbs0xa4CVlL4pJMUlnfe5pRVVPPfpDv711V4MA1oG+TBrTHdu7BGp3g4RETeh8CFOtSp7FclrkzEwah3PK80jeW0y8wfPP2cA+XrXUVKWZLD/eCkA43u14olRXWke6NPgdYuIiOMofIjT2Ow2UjemnhU8AAwMTJiYu3EuQ2KG1LoFU1hWyZyPfuDdjQcAiAr145lxPRjSOcJptYuIiOMofIjTpOWl1brVciYDg5zSHNLy0ugT2QeA1T/k8qelW8kpLAPgtqtimTa8M8F+2vZeRMRdKXyI0+SX5l90u2PF5cz6IJPlWw4D0KZFAHMnJNCvbYuGLFFERJxA4UOcJjwg/KLaZR2G6W+v53hJBWYT3DOoLQ9d31EbwYmIeAiFD3GaxIhErAFW8krz6hz3AeBDGH9dYQMMOkcGM3dCAj1jmjm1ThERaVjaZUucxmK2kNI3BQATdUyLNaDg4I14Wyw8lNSR5VOuVvAQEfFACh/iVElxScwfPJ+IgNozVeyVoZw6dBvdml3Nh/cP4oGkDvh46Y+niIgn0m0XcbqkuCSuaTWYZ9Z8xMK0bZSXB+Jd2Y6UG7pw58B4LNoITkTEoyl8iNPtzC1i6uJ0/rffAiRwVdswUscn0KZloKtLExERJ1D4EKeptNl5Ze1u/v75LipsdoJ9vXjsxi78sk+Mtr0XEWlCFD7EKTIOFvDooi1szykCYGjnCJ4e152oUH8XVyYiIs7m8BF9M2fOxGQy1frq3Lmzoy8jbqKs0kbqx9sZ+/JXbM8ponmAN3/95RX83+1XKniIiDRRDdLz0a1bN1atWvXTRbzUwdIUbdx7nGmL09l7tASA0T2jmTG6Ky2DfF1cmYiIuFKDpAIvLy8iIyMb4qXFDRSXVzH34+38+5tsACKCfZk9rgfXd7W6uDIREWkMGiR87Ny5k+joaPz8/Ojfvz9z5swhNja2zrbl5eWUl5fXPC4sLGyIkjyH3QbZX0NxLgRZIW4AmBvPsuPrduTz+JIMDp08BcAv+8Tw2I1dCPXXRnAiIlLNZBhG3etcX6KPP/6Y4uJiOnXqxJEjR5g1axaHDh1i69atBAcHn9V+5syZzJo166zjBQUFhISEOLI095e5HFZOg8LDPx0LiYbhc6HrGNfVBZwsreDPKzJZknYIgJgwf1LHJzCwfUuX1iUiIs5RWFhIaGjoRX1+Ozx8nOnkyZPExcUxf/587r777rPO19XzERMTo/BxpszlsHASnLUnyo9TVG95y2UB5KOMIzy5bCtHiyswmeDOAfE8MqwjAT4a6yMi0lTUJ3w0+KdDs2bN6NixI7t27arzvK+vL76+GoB4XnZbdY9HnZuxGYAJVqZA55FOvQWTV1TGk+9vY+W2HADaRwQxd0ICveOaO60GERFxPw2+eUZxcTG7d+8mKiqqoS/lubK/rn2r5SwGFB6qbucEhmGwaPNBrp+/npXbcvAym7jvuvZ8eP/VCh4iInJBDu/5eOSRRxg9ejRxcXEcPnyYGTNmYLFY+NWvfuXoSzUdxbmObXcZDp4o5fGlW1m/Ix+A7q1CmDehJ12jdYtMREQujsPDx8GDB/nVr37FsWPHCA8P5+qrr+abb74hPDzc0ZdqOoIucorqxba7BHa7wf/7Npu5H2+npMKGj5eZh5I6cs+geLws2n1WREQunsPDx3vvvefol5S4AdWzWgqPUPe4D1P1+bgBDXL53fnFpCxOZ9O+EwD0adOc1AkJtAsPapDriYiIZ9N0BHdgtlRPp104ierZLT8PID/Odhme6vDBplU2O69/sZfnV+2gospOgI+FlBGdua1fnDaCExGRS6bw4S66jqmeTlvnOh+pDp9mm3m4kKmLt7D1UPWib9d0DOeZcd1p3TzAodcREZGmR+HDnXQdUz2dtgFXOC2vsvH3z3fxytrdVNkNQv29eWJUVyYktsJkUm+HiIhcPoUPd2O2QPygBnnpzdknmLY4nV15xQAM7xbJn8d2IyLYr0GuJyIiTZPCh1BaUcWzn2Txxtf7MAxoGeTLUzd1Y0QPrc0iIiKOp/DRxH216ygpS9I5cLx6I7gJia15YlQXmgX4uLgyERHxVAofTVTBqUrmfPQD7206AECrZv7MHtedwZ0iXFyZiIh4OoWPJuizzFymv59BbmH1hn6T+scxdXhngnz1x0FERBqePm2akGPF5cxYvo0V6UcAiG8ZyNwJCfSND3NxZSIi0pQofDQBhmGwfMthZi7fxonSSixmE/cMasuDSR3w83beLrgiIiKg8OF2bHYbaXlp5JfmEx4QTmJEIpbzrPNxpOAU05duZfX2PAA6Rwbz7MSe9Ggd6qySRUREalH4cCOrsleRujGV3NKfdq+1BlhJ6ZtCUlxSrbZ2u8F7mw4w56MfKCqvwsdi5r7r2vO7a9vh46WN4ERExHUUPtzEquxVJK9NxjhjY7m80jyS1yYzf/D8mgCSfayEaYvT+WbPcQCuiGnGsxMT6GANdnrdIiIiZ1L4cAM2u43UjalnBQ8AAwMTJuZunMs1rQbz1ob9/OXTLMoq7fh5m3nkhk7cOTAeizaCExGRRkLhww2k5aXVutVyJgODnNIcRv3fW+zYZwVgQLsWpI5PILaFNoITEZHGReHDDeSX5l9Uuz3HjxDs24o/jezCL/rEaCM4ERFplBQ+3EB4QPhFtUtsFcuLY68lMlQbwYmISOOlaQ9uIDEiEWuAFRPn6MkwINQ7nHdu+5WCh4iINHoKH27AYraQ0jfl3A1MMHPg43hZ1JElIiKNn8KHm+hnvZYrAx7AXll7cbDIgEieH/z8Wet8iIiINFb6p7IbWJOVx5+WZHC4IBKYxvWJJQzvGUhsaOQFVzgVERFpbBQ+GrETJRU8tSKTJf87BEBsWACp43swoH1LF1cmIiJy6RQ+GiHDMPgoI4cZy7dytLgCkwnuGhjPwzd0JMBHPzIREXFv+iRrZPIKy3hi2VY+2Va9qFiHiCDmTkwgMba5iysTERFxDIWPRsIwDP67+SBPr8iksKwKL7OJPw5ux+Tr2uPrpTEdIiLiORQ+GoEDx0t5fGkGX+w8CkCPVqHMm5hAl6gQF1cmIiLieAofLmS3G7y1YR/zPsmitMKGr5eZh67vyG+vjsfLolnQIiLimRQ+XGR3fjHTFqXzXfYJAPq2CSN1Qg/ahge5uDIREZGGpfDhZJU2O6+t38NfV++kospOoI+FlBGdubVfHGZtey8iIk2AwocTbTtcwNRF6Ww7XAjAtR3DeWZ8D1o183dxZSIiIs6j8OEEZZU2/v75Ll5dt5squ0GovzdPjurK+MRW2vZeRESaHIWPBrY5+zhTF6WzO78EgBt7RDJrTHfCg31dXJmIiIhrKHw0kJLyKp79JIs3N+zDMKBlkC9Pj+3G8O5Rri5NRETEpRQ+GsCXO4+SsiSdgydOATCxd2umj+xCswAfF1cmIiLiegofDlRwqpLZH2ay8LuDALRq5s8z43twbcdwF1cmIiLSeCh8OMin23KY/v5W8orKAbi9fxyPDu9MkK/eYhERkZ/TJ+NlOlpczszl21iRfgSAti0DmTsxgT5twlxcmYiISOOk8HGJDMNg2feHmfXBNk6UVmIxm7j3mrY8MLQDft7aCE5ERORcFD4uweGTp5j+/lY+354HQJeoEOZNSKBH61AXVyYiItL4KXzUg91u8O6m/cz5aDvF5VX4WMzcP7Q9v7u2Hd7aCE5EROSiKHxcpH1HS0hZks43e44D0Cu2GfMmJNDBGuziykRERNyLwscF2OwG//pyL899lkVZpR1/bwuPDuvE7QPaYNFGcCIiIvWm8HEeWTlFTF20hS0HCwAY2L4Fc8YlENsiwMWViYiIuC+FjzpUVNl5ee0uXlqzi0qbQbCvF38a2YVf9InRRnAiIiKXSeHjDFsOnGTqonSycosASOpi5emx3YkM9XNxZSIiIp6hyYQPm91GWl4a+aX5hAeEkxiRiMX803ocpypsPL9qB//3xR7sBrQI9GHmmG6MSohSb4eIiIgDNYnwsSp7FakbU8ktza05Zg2wktI3haS4JL7Zc4yUxensO1YKwE1XRDNjdDfCArURnIiIiKN5fPhYlb2K5LXJGBi1jueV5vHQ2mSu9H+ANWmRAESG+DF7XHeGdrG6olQREZEmwaNXxrLZbaRuTD0reADVxwyDjYULADu/6hvLp8nXKHiIiIg0MI/u+UjLS6t1q+UsJjB7FzDjZn/u7N3DeYWJiIg0YR7d85Ffmn9R7azNKxq4EhERETnNo8NHeEC4Q9uJiIjI5fPo8JEYkYg1wIqJuqfKmjARGRBJYkSikysTERFpujw6fFjMFlL6pgCcFUBOP57Wd1qt9T5ERESkYXl0+ABIikti/uD5RARE1DpuDbAyf/B8kuKSXFSZiIhI0+TRs11OS4pLYkjMkPOucCoiIiLO0STCB1TfgukT2cfVZYiIiDR5Hn/bRURERBoXhQ8RERFxKoUPERERcSqFDxEREXGqJjPgFLsNsr+G4lwIskLcANBsFxEREadrsJ6Pl156iTZt2uDn50e/fv3YuHFjQ13qwjKXwwvd4c1RsPju6v++0L36uIiIiDhVg4SP//znPyQnJzNjxgzS0tLo2bMnw4YNIy8vryEud36Zy2HhJCg8XPt44ZHq4wogIiIiTtUg4WP+/Pncc8893HnnnXTt2pVXX32VgIAA/vWvfzXE5c7NboOV0wCjjpM/HluZUt1OREREnMLh4aOiooLNmzeTlPTTsuVms5mkpCQ2bNhwVvvy8nIKCwtrfTlM9tdn93jUYkDhoep2IiIi4hQODx9Hjx7FZrNhtVprHbdareTk5JzVfs6cOYSGhtZ8xcTEOK6Y4lzHthMREZHL5vKpto899hgFBQU1XwcOHHDciwdZL9ymPu1ERETksjl8qm3Lli2xWCzk5tbuTcjNzSUyMvKs9r6+vvj6+jq6jGpxAyAkunpwaZ3jPkzV5+MGNMz1RURE5CwO7/nw8fGhd+/erF69uuaY3W5n9erV9O/f39GXOz+zBYbP/fGB6YyTPz4enqr1PkRERJyoQW67JCcn8/rrr/Pmm2/yww8/8Ic//IGSkhLuvPPOhrjc+XUdA7e8BSFRtY+HRFcf7zrG+TWJiIg0YQ2ywukvfvEL8vPzefLJJ8nJyeGKK65g5cqVZw1CdZquY6DzSK1wKiIi0giYDMOoazCEyxQWFhIaGkpBQQEhISGuLkdEREQuQn0+v10+20VERESaFoUPERERcSqFDxEREXEqhQ8RERFxKoUPERERcSqFDxEREXEqhQ8RERFxKoUPERERcSqFDxEREXGqBlle/XKcXnC1sLDQxZWIiIjIxTr9uX0xC6c3uvBRVFQEQExMjIsrERERkfoqKioiNDT0vG0a3d4udrudw4cPExwcjMlkcuhrFxYWEhMTw4EDB7RvTAPS++wcep+dQ++z8+i9do6Gep8Nw6CoqIjo6GjM5vOP6mh0PR9ms5nWrVs36DVCQkL0B9sJ9D47h95n59D77Dx6r52jId7nC/V4nKYBpyIiIuJUCh8iIiLiVE0qfPj6+jJjxgx8fX1dXYpH0/vsHHqfnUPvs/PovXaOxvA+N7oBpyIiIuLZmlTPh4iIiLiewoeIiIg4lcKHiIiIOJXCh4iIiDhVkwkfL730Em3atMHPz49+/fqxceNGV5fkcebMmUOfPn0IDg4mIiKCsWPHkpWV5eqyPF5qaiomk4kHH3zQ1aV4nEOHDnHbbbfRokUL/P396dGjB999952ry/IoNpuNJ554gvj4ePz9/WnXrh1PPfXURe0PIue3fv16Ro8eTXR0NCaTiffff7/WecMwePLJJ4mKisLf35+kpCR27tzplNqaRPj4z3/+Q3JyMjNmzCAtLY2ePXsybNgw8vLyXF2aR1m3bh2TJ0/mm2++4bPPPqOyspIbbriBkpISV5fmsTZt2sQ//vEPEhISXF2Kxzlx4gQDBw7E29ubjz/+mMzMTJ577jmaN2/u6tI8yty5c3nllVf4+9//zg8//MDcuXOZN28ef/vb31xdmtsrKSmhZ8+evPTSS3WenzdvHi+++CKvvvoq3377LYGBgQwbNoyysrKGL85oAvr27WtMnjy55rHNZjOio6ONOXPmuLAqz5eXl2cAxrp161xdikcqKioyOnToYHz22WfGtddeazzwwAOuLsmjTJs2zbj66qtdXYbHGzlypHHXXXfVOjZ+/Hjj1ltvdVFFngkwli5dWvPYbrcbkZGRxrPPPltz7OTJk4avr6/x7rvvNng9Ht/zUVFRwebNm0lKSqo5ZjabSUpKYsOGDS6szPMVFBQAEBYW5uJKPNPkyZMZOXJkrT/b4jjLly/nyiuv5OabbyYiIoJevXrx+uuvu7osjzNgwABWr17Njh07ANiyZQtffvklI0aMcHFlnm3v3r3k5OTU+v0RGhpKv379nPLZ2Og2lnO0o0ePYrPZsFqttY5brVa2b9/uoqo8n91u58EHH2TgwIF0797d1eV4nPfee4+0tDQ2bdrk6lI81p49e3jllVdITk7m8ccfZ9OmTdx///34+Phw++23u7o8j5GSkkJhYSGdO3fGYrFgs9mYPXs2t956q6tL82g5OTkAdX42nj7XkDw+fIhrTJ48ma1bt/Lll1+6uhSPc+DAAR544AE+++wz/Pz8XF2Ox7Lb7Vx55ZU888wzAPTq1YutW7fy6quvKnw40MKFC3n77bd555136NatG99//z0PPvgg0dHRep89mMffdmnZsiUWi4Xc3Nxax3Nzc4mMjHRRVZ5typQprFixgjVr1tC6dWtXl+NxNm/eTF5eHomJiXh5eeHl5cW6det48cUX8fLywmazubpEjxAVFUXXrl1rHevSpQv79+93UUWe6dFHHyUlJYVf/vKX9OjRg9/85jc89NBDzJkzx9WlebTTn3+u+mz0+PDh4+ND7969Wb16dc0xu93O6tWr6d+/vwsr8zyGYTBlyhSWLl3K559/Tnx8vKtL8khDhw4lIyOD77//vubryiuv5NZbb+X777/HYrG4ukSPMHDgwLOmiu/YsYO4uDgXVeSZSktLMZtrfxRZLBbsdruLKmoa4uPjiYyMrPXZWFhYyLfffuuUz8YmcdslOTmZ22+/nSuvvJK+ffvywgsvUFJSwp133unq0jzK5MmTeeedd1i2bBnBwcE19w1DQ0Px9/d3cXWeIzg4+KxxNIGBgbRo0ULjaxzooYceYsCAATzzzDPccsstbNy4kddee43XXnvN1aV5lNGjRzN79mxiY2Pp1q0b//vf/5g/fz533XWXq0tze8XFxezatavm8d69e/n+++8JCwsjNjaWBx98kKeffpoOHToQHx/PE088QXR0NGPHjm344hp8Pk0j8be//c2IjY01fHx8jL59+xrffPONq0vyOECdXwsWLHB1aR5PU20bxgcffGB0797d8PX1NTp37my89tprri7J4xQWFhoPPPCAERsba/j5+Rlt27Y1/vSnPxnl5eWuLs3trVmzps7fybfffrthGNXTbZ944gnDarUavr6+xtChQ42srCyn1GYyDC0jJyIiIs7j8WM+REREpHFR+BARERGnUvgQERERp1L4EBEREadS+BARERGnUvgQERERp1L4EBEREadS+BARERGnUvgQERERp1L4EBEREadS+BARERGnUvgQERERp/r/vUiELRVrJ3AAAAAASUVORK5CYII=", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "plt.plot(x,y,label='verité')\n", "plt.plot(x,y_bruit,'o',label='données')\n", "plt.plot(x,y_pred,'o',label='modele')\n", "plt.legend()\n" ] }, { "cell_type": "code", "execution_count": null, "id": "e7c8e687", "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "id": "35c33046", "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "mlenv", "language": "python", "name": "mlenv" }, "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.10.6" } }, "nbformat": 4, "nbformat_minor": 5 }