From 2a3815f19efe9575e01b54777c03e928191ee496 Mon Sep 17 00:00:00 2001 From: 3eb2734abf9c03158722854f9624a07d <3eb2734abf9c03158722854f9624a07d@app-learninglab.inria.fr> Date: Thu, 8 Dec 2022 14:12:48 +0000 Subject: [PATCH] prediction --- module3/exo3/Correction.ipynb | 64 +++++++++++++++++++++++++++++++++++ 1 file changed, 64 insertions(+) diff --git a/module3/exo3/Correction.ipynb b/module3/exo3/Correction.ipynb index 9498488..45bea42 100644 --- a/module3/exo3/Correction.ipynb +++ b/module3/exo3/Correction.ipynb @@ -205,6 +205,70 @@ "plt.show()" ] }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "## Prédiction" + ] + }, + { + "cell_type": "code", + "execution_count": 10, + "metadata": {}, + "outputs": [], + "source": [ + "slope, intercept, r_value, p_value, std_err = stats.linregress(data_year['year'].dt.year, data_year['concentration'])\n", + "\n", + "def predict(x):\n", + " return slope * x + intercept" + ] + }, + { + "cell_type": "code", + "execution_count": 11, + "metadata": {}, + "outputs": [ + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAA6gAAAJCCAYAAADN6ep4AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzs3Xd4VvX9//HnSSBhhb03ypI9AqiIioATK1IXinu1atVaa1VEBHHVWlddaHGiuK11lA0O9t57y96EDDLO74/k+6t1MnLnvpM8H9fFZXLn3OfzQjBXXp5z3p8gDEMkSZIkSYq2uGgHkCRJkiQJLKiSJEmSpBhhQZUkSZIkxQQLqiRJkiQpJlhQJUmSJEkxwYIqSZIkSYoJFlRJkiRJUkywoEqSJEmSYoIFVZIkSZIUE0pEOwBA1apVw4YNG0Y7hiRJkiQpAmbNmrUjDMNqv3ZcTBTUhg0bMnPmzGjHkCRJkiRFQBAE6w7lOG/xlSRJkiTFBAuqJEmSJCkmWFAlSZIkSTHBgipJkiRJigkWVEmSJElSTLCgSpIkSZJiggVVkiRJkhQTLKiSJEmSpJhgQZUkSZIkxQQLqiRJkiQpJlhQJUmSJEkxwYIqSZIkSYoJh1xQgyCID4JgThAEn+V9/ngQBEuDIJgfBMHHQRBU/N6x9wRBsDIIgmVBEJwRieCSJEmSpKLlcK6g3gYs+d7nY4BWYRi2AZYD9wAEQdACuARoCZwJPB8EQXz+xJUkSZIkFVWHVFCDIKgLnAO88n+vhWE4OgzDrLxPpwJ18z4+DxgZhmFGGIZrgJVA5/yLLEmSJEkqig71CupTwF1Azs98/Rrgy7yP6wAbvve1jXmv/Y8gCG4IgmBmEAQzt2/ffogxJEmSJElF1a8W1CAIegPbwjCc9TNfHwBkASP+76WfOCz80QthOCwMw+QwDJOrVat2GJElSZIkSUVRiUM4pivwmyAIzgZKAeWDIHgrDMP+QRBcCfQGeoRh+H8ldCNQ73vvrwtsys/QkiRJkqSi51evoIZheE8YhnXDMGxI7vCj8Xnl9EzgL8BvwjBM/d5bPgUuCYIgMQiCRkATYHoEskuSJEmSipBDuYL6c/4BJAJjgiAAmBqG4e/CMFwUBMF7wGJyb/29OQzD7KOPKkmSJEkqyg6roIZhOBGYmPdx41847iHgoaMJJkmSJEkqXg5nH1RJkiRJUow5mPVzm60UPhZUSZIkSSqkJq/cQfe/TWT+xj3RjpIvLKiSJEmSVMjk5IQ8M24F/f85jdIJ8ZQuGR/tSPniaIYkSZIkSZIK2M6UDG5/dy5fr9hBn3a1eej81pRNLBrVrmj8LiRJkiSpGJi5dhe3vD2HXakHefj81vTrXI+8XVWKBAuqJEmSJMW4MAx5+evVPPafZdStVJqPfn8irepUiHasfGdBlSRJkqQYtif1IHe+P4+xS7ZxVquaPHZBG8qXKhntWBFhQZUkSZKkGDVvwx5uGjGbbfvTGXRuC646sWGRuqX3hyyokiRJkhRjwjDkjSnrGPr5YqonleK9G0+gff1K0Y4VcRZUSZIkSYoh+9MzufvDBXy+YDM9mlfniYvaUrFMQrRjFQgLqiRJkiTFiMWb9nHTiFls2J3G3Wc154ZuxxAXV3Rv6f0hC6okSZIkRVkYhrw1dR0Pfr6ESmVK8s71x9O5UeVoxypwFlRJkiRJiqK9qZn85cP5/GfRFk5pWo0nLmpL1XKJ0Y4VFRZUSZIkSYqSWet2ces7c9m6L517z27OdScVr1t6f8iCKkmSJEkFLCcn5IVJq/j7mOXUqViaD35/Iu3qVYx2rKizoEqSJElSAdq2P5073p3HNyt30LtNLR7u25rypUpGO1ZMsKBKkiRJUgGZtHw7f3pvLikZWTz229ZclFyPICi+t/T+kAVVkiRJkiIsMzuHv41exkuTVtOsRhLvXH88TWokRTtWzLGgSpIkSVIEbdiVyh/emcPcDXu4rEt9BvZuQamS8dGOFZMsqJIkSZIUIZ/P38zdH86HAJ6/rANnt64V7UgxzYIqSZIkSfksOydk6OeLefXbtbSvX5FnLmlPvcploh0r5llQJUmSJCkfpWdmc/vIufxn0Rau7tqQe88+jpLxcdGOVShYUCVJkiQpn+w+cJDr3pjJ7PW7ub93C645qVG0IxUqFlRJkiRJygcbdqVy5fDpbNyTxnOX+rzpkbCgSpIkSdJRWrBxL1e/NoPM7BxGXNeFTg0rRztSoWRBlSRJkqSjMGHZNm4eMZtKZRIYeUMXGld3f9MjZUGVJEmSpCP07oz13PvxQprXTOLVqzpRvXypaEcq1CyokiRJknSYwjDkqbEreHrcCk5uWo3nL+tAuUTr1dHy36AkSZIkHYbM7BwGfLyA92Zu5MKOdXm4b2u3kcknFlRJkiRJOkQHMrK4acRsJi3fzq09mvDHnk0IgiDasYoMC6okSZIkHYJt+9O55rUZLNm8n0f7tuaSzvWjHanIsaBKkiRJ0q9Yt/MAl70yjZ0pB3nlimS6N68e7UhFkgVVkiRJkn7Bup0HuGTYVNIzsxl5w/G0rVcx2pGKLAuqJEmSJP2M75fTEdcdT4va5aMdqUhz1JQkSZIk/YR1Ow/Qz3JaoLyCKkmSJEk/sH5nKv2GTSU1M5u3LacFxiuokiRJkvQ963emcsmwKZbTKLCgSpIkSVKe9TtT6fdy7pXTEdd1sZwWMAuqJEmSJAEbduWW0wMHsxhxXRda1q4Q7UjFjgVVkiRJUrG3YVcqlwybSkpGFm9dazmNFguqJEmSpGLt++V0xHVdaFXHchotFlRJkiRJxVaRKKdhGO0E+caCKkmSJKlYKhLldNV4eLk77FwV7ST5woIqSZIkqdgp9OV02xJ46wJ483xI3Qkp26KdKF+UiHYASZIkSSpIk1ft4I5355GWt5VMoSqnKdtgwsMw+3VISIJeQ6DzjVCyVLST5QsLqiRJkqRiISMrmydGL+flr1fTsEpZhl/VqfDsc5qZBlOeg2+egqw06HQdnHI3lK0S7WT5yoIqSZIkqchbtmU/t42cw9It+7msS30GnHMcZRIKQR3KyYEF78O4IbBvIzQ7O/eqadUm0U4WEYXgT0SSJEmSjkxOTsirk9fy2H+WUr5UCf55ZTI9jqsR7ViHZu23MHoAbJoDtdrC+S9Co27RThVRFlRJkiRJRdKWvenc+f48vlm5g57HVefR37aharnEaMf6dTtXwZj7YelnkFQb+rwIbS6GuKI/49aCKkmSJKnI+WLBZu75aAEHs3J4+PzW9OtcjyAIoh3rl6Xugkl/hRkvQ3winHYfHH8zJJSJdrICY0GVJEmSVGTsT89k0KeL+Gj2d7StW4EnL27HMdXKRTvWL8vKgOnD4KvHIWM/tL8cug+ApEJyK3I+sqBKkiRJKhKmr9nFHe/NZdOeNG7t0YQ/nNaYkvExfFtsGMKij2HsA7BnHTTuCb0ehBotop0saiyokiRJkgq1g1k5PD1uOS9MXEXdSmV4/3cn0rFBpWjH+mXrp+UOQNo4A6q3hP4fQeMe0U4VdRZUSZIkSYXW1n3p/O6tWcxZv4eLk+sx8NwWlEuM4ZqzazWMHQyLP4FyNeA3z0K7yyAuPtrJYkIM/8lJkiRJ0s+bsXYXv39rNqkHs3ju0g6c06ZWtCP9vLTd8NXfYNpLEF8STrkbTvwDJMb487EFzIIqSZIkqVAJw5C3pq5j8L8XU7dSad6+vgtNayRFO9ZPyzoIM16BSY9B+l5ofxl0vw/Kx3CZjiILqiRJkqRCIz0zm4GfLOT9WRvp3qwaT13SngqlS0Y71o+FISz5FMYMgt1r4JhT4fShULN1tJPFNAuqJEmSpEJh0540fvfWLOZv3MutPZpwe48mxMXF4N6m382CUQNg/RSo1hwu+yB3Qm+s78MaAyyokiRJkmLelFU7ueXt2WRk5TDs8o6c3rJmtCP92J71MG4ILHgfylaD3k9C+ysg3tp1qPw3JUmSJClmhWHIq9+u5aEvltCgShmGXZ5M4+oxNlgofS988yRMeT73Kmm3O+Gk2yExRp+LjWEWVEmSJEkxKe1gNvd8NJ9P5m6iV4sa/P2itiSViqHnTbOzYPZrMOERSN0BbS6GHvdDhbrRTlZoWVAlSZIkxZwNu1K58c1ZLNmyjzt6NeWW7o1j53nTMIQVo2H0QNixDBp0hdPfhzodop2s0LOgSpIkSYop36zYwR/emU1WTsg/r0zmtOY1oh3pv7YsyB2AtGYSVD4WLnkbmp3tAKR8YkGVJEmSFDM+m7+JW9+Zw7HVyjHsimQaVS0b7Ui59m2G8UNh7ggoXRHOfAySr4ESCdFOVqRYUCVJkiTFhFnrdnHHe/PoUL8Sr13TmXKJMVBXDh6Ab5+Byc9AThaccDOcfCeUrhTtZEVSDPyJS5IkSSru1u08wPVvzKJ2hVIMuyI5+uU0Jxvmvp171TRlC7Q8H3oMgsqNopuriLOgSpIkSYqqPakHufq1GeSEIa9e3ZnKZaN82+yqCTD6Pti6EOp2govegPpdopupmLCgSpIkSYqag1k53PjmLDbuSuOt67pE95nTbUthzMDcCb0V68MFw6FlXwcgFSALqiRJkqSoCMOQuz+cz7Q1u3jq4nZ0blQ5OkFStsPEh2HW65BQDnoNgc43QslS0clTjFlQJUmSJEXFM+NW8tGc7/hjz6b0aV+n4ANkpsHU5+HrJyEzFTpdC6fcDWWrFHwWARZUSZIkSVHw8ZyNPDl2OX071OHWHo0LdvGcHFj4AYwbAns35O5j2msIVG1SsDn0IxZUSZIkSQVq2uqd/OWDBRx/TGUe7duGoCCf8Vw3GUYNgE2zoVZb6PMCNOpWcOvrF1lQJUmSJBWY1dtTuPGtWdStXJoX+3ckoURcwSy8cxWMuR+WfgZJtaHPi9DmYogroPV1SCyokiRJkgrErgMHuea1GcQFAa9e1YmKZQpgO5nUXTDprzDjZYhPhO73wQk3Q0KZyK+tw2ZBlSRJkhRx6ZnZ3PDGTDbtTeed67vQoEqEt5PJyoDpL8NXf4WM/dD+cug+AJJqRHZdHRULqiRJkqSIyskJ+fMH85m5bjf/uLQ9HRtEcDuZMITF/4Kxg2D3Wji2B5w+FGq0iNyayjcWVEmSJEkR9eTY5fx73ibuOrMZvdvUjtxCG2fmDkDaMBWqt4D+H0LjnpFbT/nOgipJkiQpYt6fuYFnx6/k4uR6/P6UYyOzyO51MG4wLPwQylaHc5+Gdv0h3rpT2PgnJkmSJCnfZeeEPD1uBc+OX0HXxlUYen6r/N9OJm0PfP0ETHsRgng4+c/Q9TZITMrfdVRgDrmgBkEQD8wEvgvDsHcQBJWBd4GGwFrgojAMd+cdew9wLZAN3BqG4ah8zi1JkiQpRm3bl85tI+cyZfVO+naow9A+rSgZn4/buWRnwszhMPFRSNsNbfvBafdBhTr5t4ai4nCuoN4GLAHK531+NzAuDMNHgyC4O+/zvwRB0AK4BGgJ1AbGBkHQNAzD7HzMLUmSJCkGfb1iO398dy4pGVk8fkEbLkyul38nD0NY9kXufqY7V0LDbnDGQ1Crbf6toag6pP+NEQRBXeAc4JXvvXwe8Hrex68Dfb73+sgwDDPCMFwDrAQ6509cSZIkSbEoKzuHJ0Yv44rh06lUJoF/33JS/pbT72bDa71h5KUQxEG/d+HKf1tOi5hDvYL6FHAX8P2buWuEYbgZIAzDzUEQVM97vQ4w9XvHbcx77X8EQXADcANA/fr1DzO2JEmSpFixdV86f3hnDtPX7OLCjnUZfF5LyiTk07ibvRth3BCY/y6UqQrnPAEdroT4kvlzfsWUX/1bEwRBb2BbGIazgiA49RDO+VNPPoc/eiEMhwHDAJKTk3/0dUmSJEmxb9Ly7dzx7lxSD2bz94va0rdD3fw5cfo++OZJmPp87q29J90BJ90OpSrkz/kVkw7lf2t0BX4TBMHZQCmgfBAEbwFbgyColXf1tBawLe/4jcD3r+XXBTblZ2hJkiRJ0ZWVncPfxyzn+YmraFYjiecua0/j6vkwPTc7C2a/DhMehtQd0Poi6HE/VMzH24UVs361oIZheA9wD0DeFdQ7wzDsHwTB48CVwKN5//xX3ls+Bd4OguDv5A5JagJMz//okiRJkqJh8940bn1nDjPW7uaSTvUYdG5LSifEH91JwxBWjIbRA2HHMmjQFU5/H+p0yJ/QKhSO5sbwR4H3giC4FlgPXAgQhuGiIAjeAxYDWcDNTvCVJEmSioYJy7Zxx7tzycjK4amL29GnfT5s7bJlAYwaAGsmQeVj4eIR0PwcyO99UxXzgjCM/uOfycnJ4cyZM6MdQ5IkSdLPOJiVO6X3pa9W07xmEs9d1oFjq5U7upPu2wzjh8LcEVC6IpxyNyRfAyUS8ie0YkYQBLPCMEz+tePyabSWJEmSpKJq1fYUbhs5h4Xf7ePSLvW5v3cLSpU8ilt6Dx6Ab5+Byc9AThaceAt0+xOUrpR/oVUoWVAlSZIk/aQwDBk5YwND/r2YxJJxvHR5R85oWfPIT5iTDXPfzr1qmrIFWp4PPQZB5Ub5F1qFmgVVkiRJ0o/sPnCQuz+az6hFW+nauAp/v6gdNcqXOvITrpqQOwBp6wKo2wkuegPqd8m/wCoSLKiSJEmS/sfklTu447157DyQwb1nN+e6k44hLu4IBxZtWwpjBuZO6K3YAC54NffKqQOQ9BMsqJIkSZKAvEFIY5Yx7KvVNKpalleu7EqrOhWO7GQp22HiwzDrdUgoB70ehC43QonE/A2tIsWCKkmSJInV21O4beRcFny3l36d6zOw93GUSTiCupCZBlOfh6+fhKw06HRt7nTeslXyP7SKHAuqJEmSVIyFYci7MzYwOG8Q0ov9O3JmqyMYhJSTAws/gHFDYO8GaHYO9BoMVZvkf2gVWRZUSZIkqZjak3qQez5awJcLt9C1cRWeuLAdNSscwSCkdZNh1ADYNBtqtYU+L0CjbvkfWEWeBVWSJEkqhuZt2MONb85i54EM7jmrOdd3O4JBSDtXwZj7YelnUL4OnP8StL4I4uIiE1pFngVVkiRJKmbGL93KzSPmUDUpgY9vOoJBSKm7YNJfYcbLEJ8Ip90Hx98MCWUiE1jFhgVVkiRJKkZGTl/PgE8W0qJWeYZf1YlqSYcxVTcrA6a/DF/9FTL2Q/vLofsASKoRucAqViyokiRJUjEQhiFPjV3B0+NWcGqzajx3aQfKJh5iHQhDWPwvGDsIdq+FY3vA6UOhRouIZlbxY0GVJEmSirjM7BwGfLyA92Zu5MKOdXm4b2tKxh/ic6IbZ+YOQNowFaq3gMs+hCY9IxtYxZYFVZIkSSrCDmRkcfPbs5m4bDu39mjCH3s2IQgOYRjS7nUwbjAs/BDKVodzn4Z2/SHeCqHI8W+XJEmSVERt35/BNa/NYNGmvTzStzX9Otf/9Tel74Wvn4CpL0IQByf/GbreBolJkQ+sYs+CKkmSJBVBa3Yc4Mrh09m2P52Xr0imx3G/MsgoOxNmvQYTH4HUndC2H5w2ECrUKZC8ElhQJUmSpCJnzvrdXPv6TADeuf542tev9PMHhyEsHwVjBsKO5dCwW+4ApNrtCiit9F8WVEmSJKkIGbt4K7e8M5sa5Uvx2tWdaVS17M8fvHkejL4P1nwFVRrDJe9As7PgUJ5RlSLAgipJkiQVESOmrWPgJwtpVacCw6/qRNVyP7PH6b5NMH4ozH0bSleCsx6H5KshvmTBBpZ+wIIqSZIkFXJhGPL3Mct5dvxKujerxj9+bo/TjP3w7TMw+VkIs+HEP0C3P0HpigUfWvoJFlRJkiSpEEs7mM2dH8zj8/mbuTi5Hg+d34oSP9zjNDsL5r4F4x+CA9ug1W+hxyCo1CA6oaWfYUGVJEmSCqnNe9O4/o2ZLNq0j3vPbs713Y758R6nK8bmDkDathjqHQ/93oG6ydEJLP0KC6okSZJUCM1Zv5sb3pxF2sFsXvmpbWS2LMwtpqvGQ6VGcNEbcNxvHICkmGZBlSRJkgqZT+Z8x10fzqdG+URGXNeFpjWS/vvFfZthwkMwdwQkloczHoFO10GJhOgFlg6RBVWSJEkqJHJyQh4fvYwXJq6iS6PKvNC/I5XL5hXPgwdyhx99+zRkZ8LxN+UOQCpTObqhpcNgQZUkSZIKgZSMLG4fOZexS7bSr3M9Bv+mFQkl4iAnG+a9k7ttzP7N0KIP9BwElY+JdmTpsFlQJUmSpBi3YVcq178xk+Vb9/PAuS248sSGucOQVk2A0QNh6wKo2wkufB3qd4l2XOmIWVAlSZKkGDZ9zS5+99YssrJzeP2aznRrUg22LYUx98OKUVCxPlwwHFr2dQCSCj0LqiRJkhSj3p2xnvs+WUi9SmV45cpkjimdBp/9EWa9DgnloNcQ6HwjlCwV7ahSvrCgSpIkSTEmKzuHR75cyj+/WUO3JlX5xwXHUWH+S/D1k5CZCp2uhVPuhrJVoh1VylcWVEmSJCmGrN6ewj0fLWDaml1cfUJ97muwiPjh18LeDdDs7NyrplWbRDumFBEWVEmSJCkGZGRl8+LE1Tw3YSWJJeN4tXsm3dfdCnNmQ8020Od5aHRytGNKEWVBlSRJkqJs6uqd3PvxAlZvP8A1x+VwV/zblJryBSTVhj4vQpuLIS4u2jGliLOgSpIkSVGy+8BBHv5iCe/P2kiLSll802YcdVeOgPhE6D4ATrgFEspEO6ZUYCyokiRJUgELw5CP53zH0M+XkJ6WyvBms+i+9TWCFfuh/eXQ/V5IqhntmFKBs6BKkiRJBWjNjgMM+HgBk1ft4Kbqi7i97FskrFsPx/aA0x+EGi2jHVGKGguqJEmSVAAOZuXw0qRVPDthJcklVjG95gdU3zMHqh0Hl30ITXpGO6IUdRZUSZIkKcKmr9nFvR8vIG3bGkZU+xed9o+HzOpw7tPQrj/E+2O5BBZUSZIkKWLCMOSpsSsYPm4ed5X9nEvLfE58Wjyc/GfoehskJkU7ohRTLKiSJElSBGRm53Dfh3MoOe9NppT9mHLZe6BtPzhtIFSoE+14UkyyoEqSJEn5LCU9k2H/fIHrt75A45KbCOt2hTMegtrtox1NimkWVEmSJCkf7Vo5k/Uj/8gdWfPZl9QQzn2HoNlZEATRjibFPAuqJEmSlB/2bWLfF4OouPR9grAcy5Pvp+nZt0J8yWgnkwoNC6okSZJ0NDJSYPIzZH/7DKWyMnkr+A0dLx9Ky2PrRzuZVOhYUCVJkqQjkZMNc0fA+KGQspX/5JzAiKSreezac6lXuUy000mFkgVVkiRJOlyrxsPogbB1IVsrtOF3B28ivl5nXr4imUplE6KdTiq0LKiSJEnSodq2JLeYrhxDWLEBHx4zlDsXN+KsVrV48uJ2lCoZH+2EUqFmQZUkSZJ+Tco2mPAQzH4DEpPI7DGEO9d24V+LdnJ114bcd04L4uOc0isdLQuqJEmS9HMOpsKU5+DbpyArHTrfyN7Ot3Pd+6uZsXYn951zHNd1OybaKaUiw4IqSZIk/VBODswfCeMehP2boHlv6DWEhelVufW1OWzclcY/Lm1P7za1o51UKlIsqJIkSdL3rZ4EowfAlgVQuwNc8E8O1jmef0xYyXMTvqVquQTevLYzXY6pEu2kUpFjQZUkSZIAti+DMffD8v9AhXrQ9xVo9VsWb0nhzue+ZfHmffTtUIdBvVtSoUzJaKeViiQLqiRJkoq3lO0w8RGY9RoklIWeD0CX35MZl8CLE1bxzPgVVCidwMtXJNOrRY0oh5WKNguqJEmSiqfMNJj6PHz9JGSmQvI1cOrdULYqy7fu50/vzWTBd3v5TdvaDP5NS/c3lQqABVWSJEnFS04OLPwAxg2BvRug6VnQawhUa0pWdg4vT1zFk2OWk1SqBC9c1oGzWteKdmKp2LCgSpIkqfhYNxlGDYBNs6FmG+jzPDQ6GYCV21K48/15zN2wh7Nb1+TB81pRpVxilANLxYsFVZIkSUXfzlW5A5CWfgZJtaHPC9DmEoiLIzsnZPg3a3h89DLKJMTzbL/29G5TiyAIop1aKnYsqJIkSSq6UnfBV4/D9JchPgG6D4ATboGEMgCs2XGAP78/j5nrdtOrRQ0eOr8V1ZNKRTm0VHxZUCVJklT0ZB2EGS/DpL9Cxj5o3z+3nCbVBCAzO4eXv17N02NXkFgijicvbkufdnW8aipFmQVVkiRJRUcYwpJPYcwg2L0Gjj0NTh8KNVr+/0PmbtjD3R/OZ+mW/ZzZsiaDz2tJjfJeNZVigQVVkiRJRcPGWTB6AKyfAtWOg8s+hCY9//+XD2Rk8bfRy3ht8lpqJJXipcs7ckbLmlEMLOmHLKiSJEkq3Pash7GDc7eOKVsNej8F7S+H+P/+qDth6Tbu+2Qhm/amcfnxDfjzGc1IKlUyiqEl/RQLqiRJkgqn9L3w9d9h6gsQBNDtTjjpdkhM+v+HbN+fweB/L+Kz+ZtpUr0cH/zuBDo2qBzF0JJ+iQVVkiRJhUt2Jsx6DSY+Aqk7oW0/OG0gVKjz/w8Jw5D3Zm7goc+XkJ6Zwx29mvK7U44loURc9HJL+lUWVEmSJBUOYQjLR8GYgbBjOTTsljsAqXa7/zls9fYU7v14AVNX76Jzw8o83Lc1jauXi1JoSYfDgipJkqTYt3kejL4P1nwFVRrDJe9As7Nyb+3NczArh2FfreKZ8StJLBHHI31bc3FyPeLi3DpGKiwsqJIkSYpd+zbBuAdh3jtQuhKc9TgkXw3x/zvgaNa6Xdz70UKWbd3P2a1r8sC5Lanu1jFSoWNBlSRJUuzJSIFvn4bJz0KYDSf+Abr9CUpX/J/D9qZl8tf/LGXEtPXUrlCKV65IpmeLGlEKLeloWVAlSZIUO3KyYc5bMOEhSNkKLftCz0FQqeH/HBaGIZ/N38zgfy9m14EMrj2pEXf0akrZRH+8lQoz/wuWJElSbFg5DkYPhG2LoG5nuHgE1Ov0o8M27Epl4L8WMnHZdlrXqcBrV3eiVZ0KUQgsKb9ZUCVJkhRd25bkDkBaORYqNoALX4MWff5nABJAZnYOw79Zw5NjlxMXBAzs3YIrT2hAiXi3jpGKCguqJEmSoiNlW+6tvLPfgMSk3C1jOt9wHDvGAAAgAElEQVQAJRJ/dOjcDXu456MFLNm8j57H1WDIeS2pXbF0FEJLiiQLqiRJkgpWZhpM+Qd88xRkpeeW0lP+AmUq/+jQ/emZ/G3UMt6Yuo4aSaV4sX9HzmhZgyBw6xipKLKgSpIkqWDk5MCC92DcENj3HTTvDT0HQ9XGP3n4fxZuYdCnC9m2P4MrT2jIn05vSlKpkj95rKSiwYIqSZKkyFv7DYwaAJvnQq120HcYNDzpJw9Nz8xmyGeLeXvaeo6rVZ6XLk+mXb2KP3mspKLFgipJkqTI2bECxgyCZZ9D+Tpw/jBofSHE/fRgow27Urn57dnM37iXG085hjtPb0ZJhyBJxYYFVZIkSfnvwE6Y9CjMHA4lSsFpA+GEm6Hkzw82mrBsG398dy7ZOSHDLu/I6S1rFmBgSbHAgipJkqT8k5kO016Er5+AgynQ8So49R4oV/1n35KdE/L02OU8O2ElzWok8WL/jjSsWrbgMkuKGRZUSZIkHb0whIUfwtjBsHc9NDkDeg2B6s1/8W27DhzktpFz+HrFDi7oWJcHz2tF6YT4AgotKdb8akENgqAU8BWQmHf8B2EYDgqCoB3wIlAKyAJuCsNwet577gGuBbKBW8MwHBWh/JIkSYq2dVNg9AD4bhbUaA3n/QuOOfVX3zZ7/W5uHjGbnQcO8mjf1lzcqZ7bx0jF3KFcQc0ATgvDMCUIgpLAN0EQfAkMAQaHYfhlEARnA38FTg2CoAVwCdASqA2MDYKgaRiG2RH6PUiSJCkadq6CsYNgyb8hqRac9zy0vQTifvkKaBiGvDl1HQ9+tpga5Uvx0e9PpFWdCgUUWlIs+9WCGoZhCKTkfVoy71eY96t83usVgE15H58HjAzDMANYEwTBSqAzMCUfc0uSJClaUnfBpL/CjFcgPgG6D8gdgJTw68+NHsjI4p6PFvDpvE30aF6dv1/Ujgpl3NtUUq5DegY1CIJ4YBbQGHguDMNpQRDcDowKguBvQBxwYt7hdYCp33v7xrzXfnjOG4AbAOrXr3/EvwFJkiQVkKwMmP4yfPVXyNgP7S+H7vdC0qFN2125LYXfvzWLVdtT+PMZzfj9KccSF+ctvZL+65AKat7tue2CIKgIfBwEQStyy+UfwzD8MAiCi4B/Aj2Bn/ouE/7EOYcBwwCSk5N/9HVJkiTFiDCExf/KvZ1391o49jQ4fSjUaHnIp/h8/mbu+mAepUrG8+a1XejauGrk8koqtA5rim8YhnuCIJgInAlcCdyW96X3gVfyPt4I1Pve2+ry39t/JUmSVJhsnAmjBsCGqVC9BfT/EBr3POS3Z2bn8MgXSxn+7Ro61K/I85d1pGaFUhEMLKkwO5QpvtWAzLxyWprcq6SPkVs6TwEmAqcBK/Le8inwdhAEfyd3SFITYHr+R5ckSVLE7F4H4wbnbh1Ttjqc+zS06w/xh359Y8vedG55ezYz1+3m6q4Nuees40goERfB0JIKu0P5DlMLeD3vOdQ44L0wDD8LgmAP8HQQBCWAdPKeJw3DcFEQBO8Bi8ndfuZmJ/hKkiQVEul74esnYOqLEMTByX+GrrdBYtJhnWbyqh3c+s4cUg9m82y/9pzbtnaEAksqSoLcIb3RlZycHM6cOTPaMSRJkoqv7EyY9RpMfARSd0LbfnDaQKjwo1mXvygMQ16ctJrHRy2lUdWyvNi/I01qHF65lVT0BEEwKwzD5F877rCeQZUkSVIRE4awfBSMGQg7lkPDbrkDkGq3O+xT7UvP5M735jF68VbOaVOLx37bhnKJ/rgp6dD5HUOSJKm42jwPRt8Ha76CKo3hkneg2VkQHP7WL0s27+P3b81i4+407u/dgqu7NiQ4gvNIKt4sqJIkScXNvk0w7kGY9w6UrgRnPQ7JV0N8ySM63YezNjLgkwVUKF2SkTccT3LDyvkcWFJxYUGVJEkqLjJS4NunYfKzEGbDiX+Abn+C0hWP7HRZ2Qz592JGTFvP8cdU5tl+HaiWlJjPoSUVJxZUSZKkoi4nG+aOgPFDIWUrtOwLPQdBpYZHfMqNu1O5acRs5m/cy+9OOZY7T29KiXi3kJF0dCyokiRJRdmq8TB6IGxdCHU7w8UjoF6nIz5dGIZ8vmAz932ykOzskJcu78gZLWvmY2BJxZkFVZIkqSjatiR3ANLKsVCxAVz4GrToc0QDkP7P0i37eODTRUxdvYuWtcvz3KUdaFi1bP5lllTsWVAlSZKKkpRtMOFhmP06JCTlbhnT+QYoceTPhu5NzeTJsct5c+o6kkqVYGifVvTrXJ/4OKf0SspfFlRJkqSiIDMNpjwH3zwJWem5pfSUv0CZI5+om50T8u6MDTw+ail70zK5rEsD/nR6UyqWScjH4JL0XxZUSZKkwiwnBxa8B+OGwL7voHlv6DkYqjY+qtPOXLuLQZ8uYtGmfXRuVJkHzm1Ji9rl8ym0JP00C6okSVJhtfYbGDUANs+FWu2g7zBoeNJRnXLL3nQe/XIJn8zdRK0KpXi2X3t6t6lFcBTPrkrSobKgSpIkFTY7VsLYQbD0MyhfB84fBq0vhLgj3+YlIyub4d+s5dnxK8jKDrmle2Nu6n4sZRL8cVFSwfE7jiRJUmGRugsmPgoz/wklSsFpA+GEm6Fk6aM67YSl2xj870Ws3ZlKrxY1uO+c42hQxem8kgqeBVWSJCnWZWXAtJfgq7/Bwf3Q4QroPgDKVT+q04ZhyJNjV/DMuBUcU60sr1/TmVOaVsun0JJ0+CyokiRJsSoMYdHHMPYB2LMOGveEXg9CjRZHfeqs7BwG/msR70xfz0XJdRnapzUJJY78FmFJyg8WVEmSpFi0YXruAKSN06F6S+j/ETTukS+nTs/M5tZ35jB68VZu7n4sd57ezCFIkmKCBVWSJCmW7FoD4wbnXjktVwPOfQba94e4+Hw5/d60TK5/fSYz1u3igXNbcFXXRvlyXknKDxZUSZKkWJC2O/cZ0+nDIIiHU/4CJ94KieXybYkte9O5cvh0Vu9Iyds+pna+nVuS8oMFVZIkKZqyDuZO5Z30GKTtgXaXwWkDoHz+lseV21K4cvh09qZl8trVnenauGq+nl+S8oMFVZIkKRrCMHcf0zH3w67V0OgUOH0o1GqT70vNWb+ba16bQXxcwMgbjqdVnQr5voYk5QcLqiRJUkH7bhaMug/WT4ZqzeHS96FJL4jAoKIJy7Zx01uzqV4+kTeu6ez+ppJimgVVkiSpoOxZD+OGwIL3oWw16P0ktL8C4iPzI9mHszZy14fzOa5WEq9e1ZlqSYkRWUeS8osFVZIkKdLS98LXf4epL+ReJe12J5x0OyQmRWzJYV+t4uEvltK1cRVeujyZcon+2Ccp9vmdSpIkKVKys2DWqzDxEUjdCW0ugR4DoULdiC2ZkxPy8BdLeOWbNfRuU4snLmpLYon82aJGkiLNgipJkpTfwhCWj4IxA2HHcmhwEpwxFGq3j+iy+9MzueejBXw2fzNXndiQ+3u3IC4u/59rlaRIsaBKkiTlp83zYfQAWPMVVGkMl7wNzc6OyACk75u1bje3vzuH73ancfdZzbnx5GMIIrymJOU3C6okSVJ+2LcJxg+FuW9D6Upw1uOQfDXEl4zoslnZOTw3YRXPjF9BrQqleP93J9CxQeWIrilJkWJBlSRJOhoZKTD5GZj8LORkwYl/gG5/gtIVI770hl2p/PHducxct5s+7WozpE8rypeKbCGWpEiyoEqSJB2JnGyYOyL3qmnKVmjZF3oOgkoNC2T5f839jvs+XkgIPHVxO/q0r1Mg60pSJFlQJUmSDteq8TB6IGxdCHU7w8UjoF6nAll6f3omg/61iI/mfEeH+hV5+pL21KtcpkDWlqRIs6BKkiQdqm1LcovpyjFQsQFc8Cq0PD/iA5D+z+z1u7l95Fw27k7l9p5NuKV7Y0rExxXI2pJUECyokiRJvyZlG0x4GGa/DglJcPpQ6HwDlEgskOWzc0Kem7CSp8c5CElS0WZBlSRJ+jmZaTDlOfjmSchKzy2lp/wFyhRcOdy4O3cQ0oy1DkKSVPRZUCVJkn4oJwcWvA/jhsC+jdC8N/QcDFUbF1iErOwc3p25gUe/XEoYwpMXt+X89nULbH1JigYLqiRJ0vet/RZGD4BNc6BWW+j7EjQ8qcCWD8OQUYu28NdRy1i9/QCdG1Xmbxe0pX4VByFJKvosqJIkSQA7VsLYQbD0MyhfB85/CVpfBHEFN4Ro+ppdPPLlEuas30Pj6uV4+Ypkeh5XnaCAhjBJUrRZUCVJUvGWugsmPQYzXoESpeC0++D4myGh4K5YLt+6n8e+XMq4pduoUT6Rx37bmt92qOuEXknFjgVVkiQVT1kZMH0YfPU4ZOyHDlfAqfdCUo0Ci7BpTxpPjlnOh7M3UjaxBHed2YyrT2xE6YT4AssgSbHEgipJkoqXMITFn8DYB2D3WmjcE3o9CDVaFFiEvamZPD9pJa99u5YwhGtPasRNpzamUtmEAssgSbHIgipJkoqPDTNyByBtmAbVW0L/j6BxjwJbPj0zmzemrOW5CavYl57J+e3rcEevptSt5AAkSQILqiRJKg52r4Wxg2HRR1CuBpz7DLTvD3EFcyvtrgMHeW/mBl6fvJbNe9M5tVk1/nJmc46rVb5A1pekwsKCKkmSiq60PfD1EzDtRQji4eS7oOttkFgu4kuHYcjs9Xt4a+o6Pp+/mYPZOXRpVJknLmrLicdWjfj6klQYWVAlSVLRk50JM1+FiY9A2m5o2y93Om+FOhFfOvVgFv+au4k3p6xj8eZ9lEssQb/O9bjs+AY0rZEU8fUlqTCzoEqSpKIjDGH5f2D0QNi5Ahp2gzMeglptI770ym0pvDV1HR/O2sj+jCya10ziofNb0addHcom+iOXJB0Kv1tKkqSiYdOc3GK69muo0gT6jYSmZ0IQRGzJzOwcxizeyptT1jFl9U4S4uM4u3VNLj+hAR3qVyKI4NqSVBRZUCVJUuG2dyOMexDmj4QyVeDsv0HHqyC+ZMSWTDuYzfBv1/D65LVs259BnYqluevMZlyUXI+q5RIjtq4kFXUWVEmSVDhl7IdvnoQpz+Xe2nvSH3N/laoQsSXDMOSLBVt4+IslfLcnjVOaVuORvg04tVl14uO8WipJR8uCKkmSCpfsLJj9eu4ApAPbofVF0GMgVKwf0WWXbN7HA58uYtqaXRxXqzx/v6gtXY6pEtE1Jam4saBKkqTCIQxhxRgYMxC2L4X6J8Kl70KdjhFddveBgzwxZhlvT1tPhdIlGdqnFf061/eKqSRFgAVVkiTFvs3zYfR9sGYSVD4WLh4Bzc+J6ACkrOwc3p6+nidGLyclI4srTmjI7T2bULFMQsTWlKTizoIqSZJi175NMP4hmDsCSleCs/4KHa+GEpEtiZNX7WDwp4tZtnU/Jx5bhUHntqRZTfcwlaRIs6BKkqTYk5ECk5+Byc9CThaceAt0uxNKV4zosht2pfLwF0v4cuEW6lYqzYv9O3BGy5puFyNJBcSCKkmSYkdOdu7V0vEPQcoWaHk+9HwAKjWM6LJpB7N5YdIqXpq0iiCAO3o15YaTj6FUyfiIritJ+l8WVEmSFBtWjYfRA2HrQqjbGS5+E+p1jviy+9IzuXL4dOas38O5bWtzz1nNqV2xdMTXlST9mAVVkiRF17YlucV05Rio2AAufA1a9InoAKT/sy89kyv+OZ2F3+3lhcs6cFbrWhFfU5L08yyokiQpOlK2wYSHc/c0TUiCXg9ClxuhRGKBLL83LZMrhk9n8aa9PH9ZB05vWbNA1pUk/TwLqiRJKliZaTDlOfjmKchKg07Xwyl/gbJVCizC3tRMLh8+jSWb9/H8ZR3p1aJGga0tSfp5FlRJklQwcnJgwfswbgjs2wjNzoFeQ6Bq4wKNsTc1k/7/nMayLft5sX9HehxnOZWkWGFBlSRJkbf2Wxg9ADbNgVptoe9L0PCkAo+xJ/Ug/f85jeVbUnjx8g6c1txyKkmxxIIqSZIiZ+cqGHM/LP0MyteB81+C1hdBXFyBR9l94CCXvTKNldtTeOmKjnRvVr3AM0iSfpkFVZIk5b/UXTDpMZjxCpQoBafdB8ffDAllohJnV145XbU9hWGXd+RUy6kkxSQLqiRJyj9ZGTB9GHz1OGTshw5XwKn3QlL0bqXddeAgl748lTU7DvDKFcmc3LRa1LJIkn6ZBVWSJB29MITFn8DYB2D3WmjcM3fbmBotohprZ0oGl70yLbecXplMtyaWU0mKZRZUSZJ0dDbMyB2AtGEaVG8J/T+Cxj2inYodKRlc+vJU1u1M5Z9XduKkJlWjHUmS9CssqJIk6cjsXgfjBsPCD6FcDTj3GWjfH+Lio52M7ftzy+mG3akMv6oTXRtbTiWpMLCgSpKkw5O2B75+Aqa9CEE8nHwXdL0NEstFOxnZOSFfr9jO0M+XsDGvnJ54rOVUkgoLC6okSTo02Zkw81WY+Aik7Ya2/XKn81aoE+1kbN6bxnszNvLezA18tyeNKmUTePWqzpxwbJVoR5MkHQYLqiRJ+mVhCMu+zN3PdOcKaNgNTh8KtdtFNVZmdg7jl25j5PT1TFq+nZwQTmpclXvObk6vFjVILBH9W40lSYfHgipJkn7e5nkwagCs/RqqNIF+I6HpmRAEUYu0bucB3p2xgfdnbWT7/gyqJyVy06mNubhTPepVjs4+q5Kk/GFBlSRJP7b3Oxg/FOa9A2Uqw9l/g45XQXzJqMTJyMpm1KKtvDtjPd+u3ElcAN2bVeeSzvXp3qwaJeLjopJLkpS/LKiSJOm/MlLg26dg8j8gzIaut0K3P0GpClGJk5KRxXMTVjJy+np2p2ZSp2Jp7ujVlAuT61KrQumoZJIkRY4FVZIkQU42zHkLJjwEKVuh1W+hxyCo1CAqccIw5MuFWxj870Vs25/BmS1r0q9zfU5qXJW4uOjdXixJiiwLqiRJxd3KsTB6IGxbDPW6wCVvQ93kqMVZvzOV+z9dyMRl22lRqzwv9u9I+/qVopZHklRwLKiSJBVXWxflFtNV46BSQ7jwdWhxXtQGIGVkZfPyV6t5dvxKSsQF3N+7BVec0MDnSyWpGLGgSpJU3Ozfknsr75y3ILE8nPEwdLoOSiRGLdKUVTu575MFrNp+gLNb1+T+3i2pWaFU1PJIkqLDgipJUnFx8EDu8KNvn4bsg9Dl93DynblTeqNkR0oGD3+xhI9mf0e9yqV59apOdG9ePWp5JEnRZUGVJKmoy8nO3S5m/FDYvzn3Nt6eD0DlY6IXKSdk5IwNPPafpaQezOKW7o25uXtjSifERy2TJCn6LKiSJBVlqybkPme6dQHUSYYLX4P6x0c10uJN+7jvkwXMXr+HLo0q89D5rWhcPSmqmSRJscGCKklSUbRtKYwZCCtGQ8X6cMFwaNk3agOQAA5m5fDU2OW89NVqKpQuyRMXtqVvhzoEUcwkSYotv1pQgyAoBXwFJOYd/0EYhoPyvvYH4BYgC/g8DMO78l6/B7gWyAZuDcNwVGTiS5Kk/5GyDSY+ArNeh4Ry0OtB6HwDlIzuwKEVW/dz+7tzWbRpHxd2rMuAc46jYpmEqGaSJMWeQ7mCmgGcFoZhShAEJYFvgiD4EigNnAe0CcMwIwiC6gBBELQALgFaArWBsUEQNA3DMDsyvwVJkkRmGkx5Dr55CrLScqfynvIXKFslqrFyckLemLKWR75cStnEErx8RTK9WtSIaiZJUuz61YIahmEIpOR9WjLvVwj8Hng0DMOMvOO25R1zHjAy7/U1QRCsBDoDU/I5uyRJysmBBe/DuCGwbyM0Owd6DYGqjaOdjK370rnz/Xl8vWIH3ZtV47EL2lA9ya1jJEk/75CeQQ2CIB6YBTQGngvDcFoQBE2BbkEQPASkA3eGYTgDqANM/d7bN+a99sNz3gDcAP+PvfsOz7I++z/+vhNCBoEQRlhhh72HIKIis27rrFqtdtrW1q7nZ6sguKjaql3ax9FaV+tordbHxVK2TBWRsAKEESAhIYEEsnP9/riptVUrwh3ujPfrODgkyXV/z/NWDpKP3+s6v9ClS5fjehOSJDVK2Utg9lTY/S50GAIXPgTdT4t2VwC8tnYPN7+4lrLKau784kC+PLqLz5pKkj7TUQXUI7fnDg2FQi2BF0Oh0MAjr00FTgZOAp4PhUI9gE/67hN8wpqPAI8AjBw58mNflyRJn6JgC8yZDhtegRad4MKHYdBlEBMT7c4oLqvk1pczeeGdXQxOT+FXXxpKz7bJ0W5LklRPfK4pvkEQFIVCofnAmYR3Rv9+5BbgFaFQqAZoc+TznT/ysnRgd2TalSSpETu8HxbcAyv/AE0SYMI0OPl6aJoU7c4AWLFtPz9+/j12F5Vyw4QMvj+xF3Gx0Q/NkqT642im+LYFKo+E00RgEnAP4edSJwDzj9zu2xTIB14G/hIKhe4nPCSpF7CilvqXJKnhqyqHFY/Awl9CeTEM/wqccTM0rxvDhv55fMz/LthC59Qk/vrtMYzo2irabUmS6qGj2UHtADxx5DnUGOD5IAheCYVCTYHHQqHQB0AFcM2R3dR1oVDoeSCT8PEz1zvBV5KkYxAEkPkSzL0VCrMhY1L42Jh2/aPd2Yey8sLHx3yQc5DLRqYz/bwBJMd7zLok6diEwpkyukaOHBmsWrUq2m1IklR37FwZHoC0czmk9Ycpd4QDah1RXRPwpyXb+OWsjSQ1jeWuiwZz5sD20W5LklRHhUKh1UEQjPys6/xfnJIk1SWF2TD3Nlj3d0huB+f9FoZdBTGx0e7sQxv3FnPjC++zZmcRE/umcddFg0hr4fExkqTjZ0CVJKkuKC2CRffB8ocgFAun3whjfwDxdWcCbkVVDb+fn8WDb2XRPCGO314xjPMGd/D4GElSxBhQJUmKpupKWPUnmH8XlBbCkCvC03lTPnaEeFS9t7OIn/7tfTbmFnPB0I7MOG8ArZo1jXZbkqQGxoAqSVI0BAFsfD18nmnBZuh2Gky5EzoOjXZn/6a0opr752zkj4u3kdY8gT9eM5KJ/erG9GBJUsNjQJUk6UTbswZmTYXsRdC6F1zxLPQ+E+rYrbJLt+TzsxfWsmP/Ya4c3YWfndWXFglx0W5LktSAGVAlSTpRDuTAm3fCmmcgqRWcfS+MuBZi61boO1hWyV2vbeCZFTvo2jqJZ755MmN6to52W5KkRsCAKklSbSsvgSW/hqUPQFANY2+AU38MiS2j3dnHzM3MZepLa9lXXM63Tu/Bjyb1JrFp3ZkgLElq2AyokiTVlppqePfp8K7poTwYeDFMnAGpXaPd2YdqagI255WwMns/8zfmMXd9Hn3bN+eRq0cypHPdC9CSpIbNgCpJUm3Imgezb4G8ddB5NFzxDKR/5vnkta6sspq1OQdYmb2fVdmFrMrez8GyKgDaNo/nx5N78+1xPWnaJCbKnUqSGiMDqiRJkZSbCXNugay5kNoNLn0C+l8QtQFIBw5XsnrHflZmF7Jy237e33WAiuoaAHq2bcbZgzowslsrTuqWSpdWSZ5pKkmKKgOqJEmRUJwLb82Ed5+C+OYwZSaM+iY0iT+hbRworWTZ1gKWZOWzfOt+NuYWA9AkJsSg9BSuHduNkV1TGdE1ldbJJ7Y3SZI+iwFVkqTjUXEYlj0Ii38NVWUw6joYd2N4Su8JUF5VzTvbi1iSlc/irHze31VETQCJcbGM7JbKuYM7cFL3VgxJb+mwI0lSnWdAlSTpWNTUwPvPwZt3wMEc6HsuTL4dWves5bIB6/cePBJIC1ixrYCyyhpiY0IMSU/he+MzGJvRhmFdUn2OVJJU7xhQJUn6vLYtgtlTYc8a6DgMLnoUuo2ttXLFZZW8tnYPizbns3RLAfsPVQCQkZbM5Sd1YWxGG0b3aEWLhLp1nqokSZ+XAVWSpKO1bxPMmQ6bXoeUznDRH8JHx8TUzk7lgdJKHl+SzWNLtnGgtJK05vGc0bstYzPaMDajDe1TEmqlriRJ0WJAlSTpsxzKh/l3wao/QdNmMOlWGP1tiEuslXJFhyt4bPE2/rQkm+LyKib1a8f143sytHNLp+xKkho0A6okSZ+msgyW/R4W3Q+Vh2HkV+GMm6BZm1opV1BSzh8Wb+PJpdkcqqjmrIHt+d6EDAZ0TKmVepIk1TUGVEmS/lNNDXzwAsy7DQ7shN5nweTboG2fWimXV1zGowu38vSyHZRVVXPOoA58f0Iv+rRvXiv1JEmqqwyokiR91PalMGsq7H4H2g+GCx6EHuNqpVTuwTIeWrCFvyzfQWV1DRcM7cT143uSkWYwlSQ1TgZUSZIACraEByBteAWad4QvPgSDv1QrA5B2F5Xy0IItPLtyJ9U1ARcO68T14zPo3qZZxGtJklSfGFAlSY3b4f2w4B5Y+QeIjYfx02DM9dA0KeKl1u46wB8Xb+WV9/cQCsElI9L5zrgMurSOfC1JkuojA6okqXGqKocVj8LCX0B5MQy7GsZPhebtIlqmuiZgTmYujy3exors/STHN+ErY7rxtVO7kZ5qMJUk6aMMqJKkxiUIIPMfMHcGFGZDz4kw5U5o1z+iZUrKq3h+5U4eX5rNjv2H6dQykWnn9OOykzrTIiEuorUkSWooDKiSpMZj16rwAKSdyyCtP1z1AmRMimyJwsM8sTSbZ1fspLi8ihFdU/nZWX2Z0r8dTWIj/zyrJEkNiQFVktTwFW4PHxnzwQvQLA3O+034lt6Y2IiVWL29kMcWb+ONdXsBOHtQB75+aneGdm4ZsRqSJDV0BlRJUsNVdgAW3QfLHoJQDJx+I4y9AeIjc4xLEATMWpfLwwu38O6OIponNOEbp3bnK6d0o1PLxIjUkCSpMTGgSpIanupKWP04zL8rPKV3yBUwYRqkdIrI8kEQ8NbGPO6bvYl1uw/StXUSt50/gEtGpNMs3m+tkiQdK7+LSpIajiCATW/A7FugYDN0Oy08AKnj0AgtH7A4K5/7Zm/ivaTrY3wAACAASURBVJ1FdGmVxH2XDuGCoR19vlSSpAgwoEqSGoY9a8IDkLIXQetecMWz0PtMCIUisvzyrQXcN3sTK7L30zElgbsuGsQlI9KJM5hKkhQxBlRJUv12cDfMuwPWPANJreDse2HEtRAbmaNc3tlRyP2zN7E4K5+05vHcdv4ALh/VmfgmkRuwJEmSwgyokqT6qbwElvwGlv4Ogurw8KPTfgIJKRFZ/oOcA9w/ZxNvbsijdbOmTDunH1ed3JWEOIOpJEm1xYAqSapfaqrh3afhrZlQkgsDL4aJMyC1a0SW37D3IL+as4lZ63JJSYzjxjP7cM2Ybg4/kiTpBPC7rSSp/siaFx6AlLcOOo+Gy/8C6SMjsnRNTcD9czbx4Pwskps24YeTevG1U7vTIiEytwpLkqTPZkCVJNV9eeth9jTImgup3eDSJ6D/BREbgHSovIofPfceszNzuWxkOjef3Y+WSU0jsrYkSTp6BlRJUt1Vkhe+lfedJyG+OUyZCaO+CU3iI1ZiV+FhvvHEKjblFjPjvP5ce0o3QhEKvpIk6fMxoEqS6p6Kw7DsQVj8a6gqg1HXwbgbw1N6I2hV9n6ue2o1FdU1PP7VUZzeu21E15ckSZ+PAVWSVHfU1MDa52He7XAwB/qeC5Nvh9Y9I17qr6t2cvOLa0lPTeIP14ykZ9vkiNeQJEmfjwFVklQ3ZC+GWVNhz3vQYShc9Ch0GxvxMtU1AXe/vp5HF23j1Iw2PHjlcFKSHIQkSVJdYECVJEVXfhbMmQ4bX4UW6XDhIzDoUoiJiXipg2WV3PDMu8zfuI9rT+nGtHP60SQ28nUkSdKxMaBKkqLjUAEsuAdW/RGaJMCEW2DM9RCXWCvlsvMP8Y0nV5Gdf4iZFw7ky6Mjc26qJEmKHAOqJOnEqiqH5Q/DwnuhohiGXwPjb4bktForuXRLPt/98zsAPPX10Yzp2brWakmSpGNnQJUknRhBAOv+DnNvhaIdkDEZptwBaf1qtezTy7Zz68vr6N6mGX+4ZiRdWzer1XqSJOnYGVAlSbVvx3KYdTPkrIJ2A+HqF6HnhFotWVpRzczXMnl62Q7G92nLb68YRvMEhyFJklSXGVAlSbVn/9bwjmnmPyC5PVzwIAy5AmJia61kEATMzszl9v/LJKeolG+d3oOfntmX2JhQrdWUJEmRYUCVJEXe4f3hZ0xXPAKxcXDGTXDK96Fp7d5eu73gEDNeXsf8jfvo0645z183hlHdW9VqTUmSFDkGVElS5FRVwMpHYcEvoOwADLsKJkyD5u1rtWxZZTW/n7+FhxZsoWlsDNPO6cc1p3QjziNkJEmqVwyokqTjFwSw/mWYMwMKt4WfL518B7QfWOul563P5db/W8fO/aWcP6QjU8/pR7sWCbVeV5IkRZ4BVZJ0fHatgllTYecySOsPV70AGZNqvezO/Ye57f8ymbs+l4y0ZP7yzdGc0rNNrdeVJEm1x4AqSTo2hdth3m3wwQvQLA3O+w0MvQpia/dbS1llNY8s3MqDb2URGxPiZ2f15Wtju9O0ibfzSpJU3xlQJUmfT9kBWHQfLHsIQjFw+o0w9gcQn1zrpRds2seMf3xAdsFhzh7Unmnn9Kdjy8RarytJkk4MA6ok6ehUV8Lqx2H+XeEpvUOuCA9ASulUq2VragIWbN7H40uyWbBpHz3aNOPJr43i9N5ta7WuJEk68QyokqT/Lghg0xsw+xYo2AzdToMpd0LHobVa9sDhSv66eidPLdvO9oLDtEmO58Yz+/D1U7sT36T2zlGVJEnRY0CVJH26PWvCA5CyF0HrXnDFs9D7TAiFaq3kut0HeOrt7bz0Xg5llTWM7JrKjyf35qyBHXzOVJKkBs6AKkn6uIO7Yd4dsOYZSGoFZ98LI66F2LhaKVdRVcPrH+zhybe3s3p7IQlxMXxxaCeuHtOVAR1TaqWmJEmqewyokqR/KS+BJb+Bpb+DoBrG3gCn/QQSaick7jlQyl+W7+CZFTvJLymna+skpp3Tj0tHdCYlqXbCsCRJqrsMqJIkqKmGd5+Gt2ZCSS4MvBgmzoDUrrVS7oOcAzz4VhazM3OpCQLG90njK2O6cnqvtsTE1N7tw5IkqW4zoEpSY5c1LzwAKW8ddB4Nl/8F0kfWSqmcolLunbWRF9/NISUxjq+f2p2rRnelS+ukWqknSZLqFwOqJDVWeeth9jTImgup3eDSJ6D/BbUyAOlgWSW/f2sLjy3ZBsC3x/XkO2f0JCXR23glSdK/GFAlqbEpyQvfyvvOkxDfHKbMhFHfhCbxES9VUVXDn5dv57fzNlN4uJILh3XiJ1N6k57qjqkkSfo4A6okNRYVh2HZg7D411BVBqOug3E3hqf0RlgQBLzxwV7ueWMD2QWHOaVna24+ux8DOzmRV5IkfToDqiQ1dDU1sPZ5mHc7HMyBvufC5Nuhdc9aKbd6eyE/f209q7cX0rtdMn+69iTO6NOWUC2enSpJkhoGA6okNWTZi2HWVNjzHnQcBhc9Ct3G1k6p/EP8YtYGXlu7l7bN47n7okFcMiKdJrExtVJPkiQ1PAZUSWqI8rNgznTY+Cq0SA8H04GXQEzkw2JBSTm/ezOLPy/fTlxsDD+c1ItvntaDZvF+i5EkSZ+PPz1IUkNyqAAW3AOr/ghNEmHidDj5uxCXGPFSJeVV/GHRVh5duJXSymq+dFIXfjSpF2ktEiJeS5IkNQ4GVElqCKrKYfnDsPBeqCiGEdfCGTdBclrES5VVVvPn5Tt48K0s9h+q4OxB7fnx5D5kpCVHvJYkSWpcDKiSVJ8FAax7EebOgKIdkDEZptwBaf0iXqqquoa/v5vDb+ZuJqeolFMz2vD/vtCHIZ1bRryWJElqnAyoklRf7VwBs26GXSuh3UC4+kXoOSHiZYIgYNa6XO6dvZGsvBKGpKfwi0sGMzajTcRrSZKkxs2AKkn1zf5tMPdWyHwJktvD+Q/A0CshJjbipZZuyeeeNzayZmcRPds246GrhvOFAe09MkaSJNUKA6ok1RelheFnTFc8AqFYGPdTOOUGiI/8s59rdx3gF7M2sGhzPh1TEvjFxYO5aHgnj4yRJEm1yoAqSXVdVQWsegwW3A2lReHd0gnToEXHiJfaVXiYu1/fwCvv7yE1KY5p5/TjqpO7khAX+d1ZSZKk/2RAlaS6KghgwyswZwbs3wLdx8GUO6HD4IiXOlRexUMLtvDIwq2EQnDDhAy+eXoPmifERbyWJEnSpzGgSlJdlLMaZk2DHUuhTR+48q/QazJE+NnPmpqAF9/N4RezNpB7sJwLhnbkp2f2pWPLyJ+bKkmS9FkMqJJUlxTtgHm3w9q/QlIbOOd+GH4NxEb+r+vV2wu5/ZVM1uwsYkh6Cr//8ghGdE2NeB1JkqSjZUCVpLqg7AAs/hW8/fvwLulpP4GxP4SEFhEvtbuolLtf38DLa3bTrkU89182hC8O7URMjJN5JUlSdBlQJSmaqith9eMw/y44XACDL4eJt0BKesRLlVZU89CCLTy8cAtBAN+fkMG3x/WkWbzfCiRJUt3gTyWSFA1BAJvegDnTIX8TdD0VvnAndBxWC6UC/vHebu55YwN7DpRxzuAO3HRWX9JTkyJeS5Ik6XgYUCXpRNuzBmZNhexF0DoDLn8G+pxVKwOQFmfl86u5m3h3RxGDOqXw2yuGcVK3VhGtI0mSFCkGVEk6UQ7uhnl3wJpnIDEVzvoljPwqxEb2KJcDhyv56+qd/Hn5DrblH6Jt83h+cclgLhme7nOmkiSpTjOgSlJtKy+BJb+Bpb+DoBrG3hAegpSQEtEy63Yf4Km3t/PSezmUVdYwomsqP5jYi7MGtSe+SWxEa0mSJNWGzwyooVAoAVgIxB+5/m9BEMz4yNf/B/gl0DYIgvwjn7sJ+DpQDdwQBMGsWuhdkuq2mmp492l4ayaU5MKAi2DSDEjtFrES5VXVvL52L0++nc07O4pIiIvhi0M7cdXJXRnYKbIBWJIkqbYdzQ5qOTAhCIKSUCgUBywOhUKvB0GwLBQKdQYmAzv+eXEoFOoPXA4MADoCc0OhUO8gCKproX9Jqpuy5sHsWyBvHXQeDV/6M3Q+KWLL5xSV8pfl23lu5U7ySyro3qYZt5zbn0uGp5OSFNlbhiVJkk6UzwyoQRAEQMmRD+OO/AqOfPwr4EbgHx95yQXAs0EQlAPbQqFQFjAKeDtSTUtSnZW3HmZPg6y54Z3SS5+A/hdEZABSEAQsySrgybezmbs+F4AJfdvxlTFdOTWjjc+XSpKkeu+onkENhUKxwGogA3gwCILloVDofCAnCII1oX//wasTsOwjH+868rn/XPNbwLcAunTpcmzdS1JdUZIXvpX3nSchvjlMuRNGfQuaxB//0uVVvLB6F0++nc2WfYdo1awp143ryZdHd/GoGEmS1KAcVUA9cnvu0FAo1BJ4MRQKDQamAlM+4fJP+l/4wcc+EQSPAI8AjBw58mNfl6R6oeIwLHsQFv8aqspg1HUw7kZIOv6jXLLySnjq7WxeeCeHkvIqhnRuyf2XDeGcwR0ceiRJkhqkzzXFNwiColAoNJ/wbbzdgX/unqYD74RCoVGEd0w7f+Rl6cDuiHQrSXVFTQ2sfR7m3Q4Hc6DvuTDpNmiTcVzLVtcEvLkhjyffzmbR5nyaxsZw7uAOfOWUbgzt3DIyvUuSJNVRRzPFty1QeSScJgKTgHuCIEj7yDXZwMggCPJDodDLwF9CodD9hIck9QJW1Er3khQN2Yth1lTY8x50GAoXPQrdxh7XkoWHKnh+1U6eWradXYWldEhJ4P99oQ9fOqkzbZKP/zZhSZKk+uBodlA7AE8ceQ41Bng+CIJXPu3iIAjWhUKh54FMoAq43gm+khqE/CyYMx02vgot0sPBdOAlEBNzzEuu232AJ5eGzy4tr6rh5B6tmHp2Pyb3b0eT2GNfV5IkqT46mim+7wPDPuOabv/x8Uxg5nF1Jkl1xeH9MP9uWPVHaJIIE6fDyd+FuMRjXnLLvhJm/GMdi7PySYyL5eIR6XxlTFf6tm8RwcYlSZLql8/1DKokNSpV5bD8YVh4L1QUw4hr4YybIDntM1/6aSqqanhowRYeeDOLhLgYpp7dj8tGdvbsUkmSJAyokvRxQQCZL8GcGVC0HXpNgcl3QFrf41p29fZCbvr7+2zKLeHcwR2Yfl5/0ponRKhpSZKk+s+AKkkftXNFeADSrhXQbiBc/RL0HH9cSxaXVfLLWRt5atl2OrRI4LFrRzKhb7sINSxJktRwGFAlCaAwG+beCutehOR2cP7vYOiXIeb4zhudvW4v0/+xjtziMq4Z043/+UIfkuP9q1eSJOmT+FOSpMattAgW3QfLH4JQLIz7KZxyA8QnH9eyeQfLmPHyOl7/YC992zfnoatHeI6pJEnSZzCgSmqcqith1WPh6bylhTD0SpgwDVp0PK5la2oCnl25k7teX095VQ3/7wt9+NbpPYjzyBhJkqTPZECV1LgEAWx8LXyeaUEWdD8dpsyEDoOPe+msvBJu/vtaVmTvZ0yP1vz8okF0b9MsAk1LkiQ1DgZUSY3H7vdg9jTIXgRtesOVz4cn9IZCx7dsUSl/XLyNp97eTmLTWH5xyWAuHZFO6DjXlSRJamwMqJIavgM58OYdsOZZSGoFZ98bPtM09vjOHt2UW8xDC7bw8nu7CYAvDu3Ez87qS9vm8RFpW5IkqbExoEpquMqLYfGv4e0Hwrf2jv0BnPZjSEg55iWDIGBldiEPL9jCvA15JMbFctXJXfnGad1JT02KYPOSJEmNjwFVUsNTXQXvPglv/RwO7YNBl8KEWyC16zEvWVMTMGd9Lg8v2MI7O4po1awpP5rUm6+M6Upqs6YRbF6SJKnxMqBKajiCADbPgTm3wL4N0GUMXPEcpI845iXLq6r5x7u7eXjhFrbsO0R6aiK3XzCAS0d0JrHp8Z2RKkmSpH9nQJXUMOxdGx6AtHU+tOoBX3oa+p57zAOQissq+cvyHTy2ZBu5B8vp36EFv71iGGcPbE8Tj4yRJEmqFQZUSfXbwT3w5p3w3p8hsSWceQ+M/Bo0Obbbbssqq3l8aTYPvpVFcVkVYzNa88tLhnBarzZO5ZUkSaplBlRJ9VN5CSz9LSz9HdRUwSnfg9N+Aompx7RcEAS8tnYvd7+xnp37S5nQN40fTurF4PSWEW5ckiRJn8aAKql+qakO75a+ORNK9sKAC2HiDGjV/ZiXfG9nEXe8ksnq7YX0bd+cp78+mlN7tYlg05IkSToaBlRJ9UfWPJh9C+Stg/RR8KWnoPOoY14up6iUX7yxgX+8t5s2yfHcfdEgLh3ZmdgYb+WVJEmKBgOqpLovb314AFLWXGjZFS59HPp/8ZgHIJWUV/HQ/C08umgrAN8bn8G3z+hJcrx/JUqSJEWTP41JqrtK8uCtmfDOkxDfHKbcCaO+BU3ij2m56pqAv63eyb2zN7GvuJwLhnbkxjP70qllYoQblyRJ0rEwoEqqeyoOw7IHYfGvoaosHErH/RSSWh3zkkuy8rnjlUw27C1mRNdUHrl6BMO6HNtAJUmSJNUOA6qkuqOmBtY+D/Nuh4M54XNMJ90GbTKOecldhYe59eVM5q7PJT01kQevHM7Zg9p7ZIwkSVIdZECVVDdkL4ZZU2HPe9BhKFz0CHQ79ZiXq6yu4bHF2/j13M2EQvDTM/vy1bHdSIiLjWDTkiRJiiQDqqToys+COdNh46vQohNc+AgMuhRiYo55ydXbC5n64lo27C1mUr923HbBAJ8zlSRJqgcMqJKi41ABLLgHVv0RmiTAhFtgzPUQd+xB8sDhSu6ZtYFnVuygfYsEHr56BF8Y0D6CTUuSJKk2GVAlnVhV5bD8YVh4L1QUw/CvwPipkJx2zEsGQcDLa3ZzxyuZ7D9UwdfGdudHk3t7bIwkSVI9409vkk6MIIB1L8LcW6FoO2RMhil3QFq/41o2O/8Qt/zjAxZtzmdIegqPf3UUAzulRKZnSZIknVAGVEm1b+eK8ACkXSsgbQBc/SL0nHBcS5ZXVfPwgq088FYWTWNjuP2CAXx5dFdiY5zOK0mSVF8ZUCXVnsLs8I7puhchuR2c/zsY+mWIOb5Jusu2FjD1xbVs2XeIcwZ1YPp5/WnXIiEiLUuSJCl6DKiSIq+0CBbdB8sfglAsjPspnHIDxCcf85IHDlcyK3Mvr7y/h4Wb9pGemsifvnoS4/sc+7OrkiRJqlsMqJIip7oSVj0G8++G0kIYeiVMmAYtOh7TcgdKK5mTmcur7+9mcVY+ldUBnVsl8sNJvbju9J4kNvVMU0mSpIbEgCrp+AUBbHw9fJ5pwWbofjpMuRM6DPncSxWXVTJ3fS6vrNnDws37qKwO6NQyka+N7c45gzswqFMKoZDPmUqSJDVEBlRJx2f3ezB7GmQvgta94IpnofeZ8DlCZEl5FfPW5/LK+3tYsGkfFVU1dExJ4Jox3Th3SEeGpBtKJUmSGgMDqqRjcyAH3rwD1jwLSa3g7HthxLUQG3dULy86XMGbG/KYtW4v8zfuo7yqhvYtErhqdFfOGdyBYZ1bEuNEXkmSpEbFgCrp8ykvhiW/gaUPQFADY38Ap/0YEj777NGd+w8zJzOXOZm5rMjeT3VNQLsW8VwxqgvnDu7A8C6phlJJkqRGzIAq6ehUV8G7T8FbP4dDeTDwEpg4HVK7fupLgiBg3e6DzD4SStfvOQhA73bJfHtcD6b0b8+gTimGUkmSJAEGVElHI2suzL4F8jKh88lwxTOQPvITL62srmHFtv0f7pTmFJUSE4KRXVsx9ex+TO7fjm5tmp3gNyBJkqT6wIAq6dPlrgsH0y3zILU7XPYk9Dv/3wYgHSitZF3OAdbmHOD9XQdYtHkfB8uqiG8Sw2m92vKDSb2Y2DeN1snxUXwjkiRJqg8MqJI+rjgX3poZvqU3vgV84edw0jc5UBHigy0FrD0SSD/IOcD2gsMfvqxTy0Qm92/PlAHtOK1XG5Ka+leMJEmSjp4/PUr6l4pD8PaDsPjXBNUV7O5zDbNaf4XVW2HtoiXs2P/vYXRQpxQuG9mZgZ1SGNQphVbNmkaxeUmSJNV3BlRJUFMdPi7mzTugeA8lPc7mu3nns/C9FsAe0lPDYfRLJ3VmUKcUBhpGJUmSVAsMqFJjt3U+zJ4Ge9dCpxEsHvpLvrUgjsS4WB68ciCn9GxNqmFUkiRJJ4ABVWqs8jbAnOmweRakdKHywke5c1s/npizg5FdW/DAlcNpn5IQ7S4lSZLUiBhQpcamJA/m3wWrn4CmzWDSbeT0vYbvPpfJmp07+Map3fnpWX2Ji42JdqeSJElqZAyoUmNRWfrhACSqSuGkr8O4n7Egp4Yf/n4lldUB//vl4Zw1qEO0O5UkSVIjZUCVGrqaGlj7V5h3OxzcBX3Ogcm3Ud0qg9/M28zv3txMn3bN+d+rRtC9TbNodytJkqRGzIAqNWTZi2HWVNjzHnQYAhc+BN1PY/+hCn7wpxUs2pzPxcPTufOLA0lsGhvtbiVJktTIGVClhig/KzwAaeOr0KITXPgwDLoMYmJ4Z0ch1//5HQoOVXDXRYO4/KTOhEKhaHcsSZIkGVClBuXwfph/N6z6IzRJgAnT4OTroWkSQRDw+JJtzHx1PR1aJvD375zCwE4p0e5YkiRJ+pABVWoIqsph+cOw8F6oKIbhX4HxUyE5DYDdRaXMfG09r76/h0n90rjv0qGkJMVFuWlJkiTp3xlQpfosCGDdizD3VijaDhmTYcodkNaPquoa5mfm8syKHby1MQ+AG8/sw7dP70lMjLf0SpIkqe4xoEr11c4V4QFIu1ZA2gC4+kXoOYGcolKem7OJ51fuZO/BMto2j+c7Z/TkSyO70KV1UrS7liRJkj6VAVWqbwqzwzum616E5HZw/u+oGnQFb23ezzOPr2T+xjwC4PRebbn1/AFM7JdGXGxMlJuWJEmSPpsBVaovSotg0b3hZ01DsXD6jeQM+BbPrdnPc79cQO7BctKax3P9+AwuG9mZzq3cLZUkSVL9YkCV6rrqSlj1WHg6b2khpQO+xJIu3+bpzEoWzFkBwBm923LHBV2Y0DeNJu6WSpIkqZ4yoEp1VRBQmfkq1bNvIeHAVjYmDeeuplOZv7oDrN5N+xYJfH9CLy4bmU56qrulkiRJqv8MqFIdknewjHd2FLF7/duM2nw/AyveZ3tNR35e9T+sbzqG4V1bMa1LS4Z3TWVwpxR3SyVJktSgGFClKCopr+KND/ayYNM+3tleSE3RLv4n7jm+FruYA6EWvNblf2D4Nczs3pYOKYnRbleSJEmqVQZU6QSrrglYnJXPi+/s4o11eymrrKF78xqmJb/O5Iq/ERMKqBr9A1JO/wlnJ6REu11JkiTphDGgSifIhr0H+fs7Obz0bg55xeW0SGjCJcM68M3kpXRZ8ytChXkw8BKYOJ2Y1K7RbleSJEk64QyoUi3KKy7j5fd288I7Oazfc5AmMSHG903jomGdmNh0LU3nfQfyMqHzyXDFM5A+MtotS5IkSVFjQJUirLSimtmZe3nx3RwWbtpHTQBDOrfktvMHcN6QjrQq2Qyzr4ct8yC1O1z2JPQ7H0KhaLcuSZIkRZUBVYqQ4rJKHl24lceWZFNSXkWnlol854yeXDgsnYy0ZCjOhXk/gXefhvjm8IWfw0nfgCbx0W5dkiRJqhMMqNJxqqiq4S/Lt/O7N7MoOFTBOYM6cNXJXRndvRUxMSGoOAwLfgmLfwXV5TDqOhh3IyS1inbrkiRJUp1iQJWOUU1NwKtr93Dv7I1sLzjMmB6t+dlZfRnSueU/L4D3noV5t0Pxbuh7Lky+HVr3jG7jkiRJUh1lQJWOwdKsfO56fQNrcw7Qt31z/vTVkzijd1tC/3yOdNsimD0V9qyBjsPg4j9At7HRbVqSJEmq4wyo0uewfs9B7n59Aws27aNjSgL3XTqELw7rRGzMkWCavxnmTIeNr0GLdLjo0fDRMTEx0W1ckiRJqgcMqNJR2FV4mPvnbOLFd3NokRDHzWf35StjupEQFxu+4FA+zL8bVj0GcUkwcTqc/F2IS4xu45IkSVI9YkCV/ouiwxU8+FYWT7y9HYBvndaD756RQUpSXPiCyjJY/r+w6H6oOAQjroUzboLkttFrWpIkSaqnDKjSfyg8VMGbG/KYk5nLws37KK2s5uLh6fxocm86tTyyIxoE8MELMPc2OLADep8ZHoDUtk90m5ckSZLqMQOqBGTnH2Lu+lxmZ+ayKns/NQG0axHPhcM6cdXJXenXocW/Lt7+dngAUs5qaD8ILvgH9DgjWq1LkiRJDYYBVY1STU3Aml1FzMnMZU5mLpvzSgDo27453xufwaT+7RjUKeVfU3kBCrbA3Bmw/v+geUf44v/C4MsdgCRJkiRFiAFVjUZZZTVLt+QzJzOXuevz2FdcTmxMiNHdW3Hl6C5M6teOzq2SPv7Cw/thwS9g5R8gtimMnwZjroemn3CtJEmSpGNmQFWDFQQBW/YdYuGmfSzYtI9lWwsor6ohOb4J4/q0ZXK/dozvk/avgUf/qaocVjwKC38B5cUw7GoYPxWatzuxb0SSJElqJAyoalCKyypZklXAgk37WLhpHzlFpQD0aNuMK0d34Yw+aZzcoxXxTWI/fZEggMyXYO6tUJgNGZNg8h3Qrv8JeQ+SJElSY2VAVb1WUxOQuecgCzbtY8HGfbyzo5CqmoDk+CaMzWjNd8f35PRebT/51t1PsnNleADSzuWQNgCu+jtkTKzdNyFJkiQJMKCqHgqCgGVb9/P8qp0s3LSPgkMVAAzs1ILrxvXg9F5tGd41lbjYzzG8qDA7fGTMur9Dcjs477cw7CqI+S87rZIkSZIiyoCqeqO6JmDWur08vGALa3YdoGVSHGf0bsu4Pm05NaMtbZvHf/5FS4tg0X2w/CEIxcLpN8LYH0B8cuTfgCRJkqT/yoCqOq+sspq/rd7Fo4u23sx95AAAFI1JREFUsr3gMN1aJzHzwoFcPDydhLhj3OGsroRVf4L5d0FpIQy9MjwAKaVTZJuXJEmSdNQMqKqzig5X8PSy7Ty+NJv8kgqGpKfwsy8PZ8qA9sTGhD57gU8SBLDxdZgzHQo2Q/fTYcqd0GFIZJuXJEmS9Ll9ZkANhUIJwEIg/sj1fwuCYEYoFPolcB5QAWwBvhoEQdGR19wEfB2oBm4IgmBWLfWvBiinqJQ/LtrGsyt3cLiimjP6tOXb43oyunsrQqFjDKYAu9+D2dMgexG07gVXPAu9z4TjWVOSJElSxBzNDmo5MCEIgpJQKBQHLA6FQq8Dc4CbgiCoCoVC9wA3AT8NhUL9gcuBAUBHYG4oFOodBEF1Lb0HNRDr9xzkkYVbeXnNbkLA+UM68s3Te9CvQ4vjW/hADrx5B6x5FpJawdn3wohrIfZTzj+VJEmSFBWfGVCDIAiAkiMfxh35FQRBMPsjly0DLjny+wuAZ4MgKAe2hUKhLGAU8HbEulaDsnr7fn73ZhbzN+4jqWks157Sja+d2p1OLROPb+HyYljyG1j6AAQ1MPYGOO0nkJASmcYlSZIkRdRRPYMaCoVigdVABvBgEATL/+OSrwHPHfl9J8KB9Z92Hfncf675LeBbAF26dPl8XatBWLa1gN/O28zSLQW0ataUn0zuzdVjutIyqenxLVxdBe89DW/OhEN5MPBimDgDUrtGpnFJkiRJteKoAuqR23OHhkKhlsCLoVBoYBAEHwCEQqGpQBXw5yOXf9IDfcEnrPkI8AjAyJEjP/Z1NUxBELAkKxxMV2Tvp01yPNPO6ceVo7uQ1DQCM7uy5sLsWyAvEzqPhiuegfSRx7+uJEmSpFr3uRJBEARFoVBoPnAm8EEoFLoGOBeYeORWYAjvmHb+yMvSgd0R6FX1WBAEzN+0j9/O28y7O4po3yKBW8/rz+Wjuhz7UTEflZsZHoC0ZR6kdoNLn4D+FzgASZIkSapHjmaKb1ug8kg4TQQmAfeEQqEzgZ8C44IgOPyRl7wM/CUUCt1PeEhSL2BF5FtXfRAEAXPX5/G7Nzfz/q4DdGqZyJ1fHMilI9OJbxKBYFqcC2/NhHefgvjm8IWfw0nfgCbxx7+2JEmSpBPqaHZQOwBPHHkONQZ4PgiCV44MP4oH5hw5+mNZEATfDoJgXSgUeh7IJHzr7/VO8K3bKqpq2HOglF2FpewqPExO4ZHfF5WSU1jKwdJK0lrE07FlIp1aJtLxw18JdGqZSPuUhI+FzZqagFnr9vLbN7NYv+cgXVolcc/Fg7hwWDpNm8REoOnD8PaDsPhXUF0Oo66DcTeGp/RKkiRJqpdC/7ozN3pGjhwZrFq1KtptNHj5JeW8tSGP7IJD7Cos/TCI5haX8dE/BjEhaN8igfTUJNJTE2mRGEfuwTJ2F5WSU1RGfkn5x9ZukxxPp5YJdDwSWJdk5bMpt4QebZpx/fgMLhjakSaxEQimNTXw/nMw73Yo3g19z4XJt0Prnse/tiRJkqRaEQqFVgdB8JnDYSIwlUZ1WUFJOW+s28ur7+9h2dYCagKIjQnRISWB9NRETu3Vhk4tE0lPTfwwkLZPSSDuv4TJsspq9h4IB9bd//xnUSk5RaVsyi1m/sZ9dGmVxG8uH8q5gzsSGxOh50C3LYLZU2HPGug4DC7+A3QbG5m1JUmSJEWdAbUBKigpZ9a6XF5du5u3t4RDafc2zfjuGRmcNag9fdo1P67dzIS4WLq1aUa3Ns0+8etBEBCK5HCi/M3hybybXocW6XDRozDwEoiJwI6sJEmSpDrDgNpA7D9UwawjO6Vvby2guiagW+skvnNGT84Z1JF+HZpHNjT+FxGrc6gAFtwNqx6DJokwcTqc/F2IS4zM+pIkSZLqFANqPVZQUs6czFxeXbuHpVv+FUq/Pa4HZw/qQP8OLU5YKI2oyjJY8TAsvA8qimHEtXDGzZDcNtqdSZIkSapFBtR65MDhSpZtK+DtLeFfG3OLAejaOonrTg+H0gEd62koBQgCWPd3mHsrFO2AXlNg8h2Q1jfanUmSJEk6AQyodVhxWSUrtu0PB9KtBWTuOUgQQEJcDCO7tuL8oR0Z17tt/Q6l/7RjeXgA0q6V0G4gXP0S9Bwf7a4kSZIknUAG1DrkcEUVK7MLPwyka3cVURNA0yYxDO/Skh9O7M2Ynq0Z0jnlY+eO1lv7t8Lc2yDzJUhuD+c/AEOvhJgG8v4kSZIkHTUDapT9c7jRa2v38PaWAqpqAprEhBjauSXXj89gTI/WDO+aSkJcAwtspYWw8F5Y/jDExsG4n8Ep34f45Gh3JkmSJClKDKhRkF9S/mEoXbZ1P9U1AV1bJ/H1U7szNqMNI7ulktS0gf6nqaqAlX+ABfdA2QEY9mUYPw1adIh2Z5IkSZKirIGmoLpnX3E5b6zby2vv72H5tn+dTVrvJ+4erSCA9S/DnBlQuA16jIcpd0L7gdHuTJIkSVIdYUCtRXnFZbzxQXindMW2/dQE0KNtM64fn8FZAzuc0LNJo2rXKpg1FXYug7Z94ct/g4xJ0BjeuyRJkqSjZkCNkIKScjbnlbA5r4Ss3GI+2H2Qd3YUEgSQkZbM9yb04uxB7enTrpGEUoDC7TDvNvjgBWjWFs79NQy7GmL9YydJkiTp40wKn0MQBOwrPhJEc4v/FUjzSth/qOLD65o1jaVXu+bcMKEX5wzuQO92zaPYdRSUFsGi+2D5QxCKhdP/H4z9AcQ3sn8PkiRJkj4XA+pRuP3/Mlmzq4jNucUcLKv68PMtEprQu11zvjCgHRlpzclIS6ZXWjIdUhIazy7pR1VXwqo/wfy7wlN6h1wBE6ZBSqdodyZJkiSpHjCgHoWt+SU0iQlx/tCO9EprTq+0ZDLaJdM2Ob5xBtH/FASw8XWYMx0KNkO308IDkDoOjXZnkiRJkuoRA+pRePyro6LdQt21+z2YPQ2yF0HrXnDFs9D7TAcgSZIkSfrcDKg6Ngdy4M07YM2zkNQKzr4XRlwLsXHR7kySJElSPWVA1edTXgxLfgNLH4CgGsbeAKf9BBJSot2ZJEmSpHrOgKqjU1MN7z4Fb86EQ3kw8GKYOANSu0a7M0mSJEkNhAFVny1rLsy+BfIyofNouOIZSB8Z7a4kSZIkNTAGVH263MzwAKQt8yC1G1z6BPS/wAFIkiRJkmqFAVUfV5wLb80M39Ib3xymzIRR34Qm8dHuTJIkSVIDZkDVv1QchrcfhCW/hqoyGHUdjLsxPKVXkiRJkmqZAVVQUwPvPxc+NuZgDvQ9FybfDq17RrszSZIkSY2IAbWx27YIZk+FPWug4zC46FHoNjbaXUmSJElqhAyojVX+ZpgzHTa+Bi3Sw8F04CUQExPtziRJkiQ1UgbUxuZQASy4G1Y9Bk0SYeJ0OPm7EJcY7c4kSZIkNXIG1MaiqhyWPwwL74WKYhhxLZxxEySnRbszSZIkSQIMqA1fEMC6F2HuDCjaAb2mhAcgpfWLdmeSJEmS9G8MqA3ZzhUw62bYtRLaDYSrX4Ke46PdlSRJkiR9IgNqQ7R/G8y9FTJfguT2cP4DMPRKiImNdmeSJEmS9KkMqA1JaWH4GdMVj0BMExj3Mzjl+xCfHO3OJEmSJOkzGVAbgupKWPnH8HTe0iIY+mWYMBVadIx2Z5IkSZJ01Ayo9VkQwIZXw+eZ7t8C3cfBlDuhw+BodyZJkiRJn5sBtb7KeQdm3wLbF0ObPnDl8+EJvaFQtDuTJEmSpGNiQK1vinbCvNth7fOQ1AbOuQ+GXwux/qeUJEmSVL+ZauqLsoOw+Few7PfhW3tP/TGc+iNIaBHtziRJkiQpIgyodV11FbzzOLx1FxzOh8Ffggm3QMvO0e5MkiRJkiLKgFpXBQFsmgVzboH8TdB1LEz5K3QaHu3OJEmSJKlWGFDroj1rYPY02LYQWmfA5X+BPmc7AEmSJElSg2ZArUsO5MCbd8KaZyAxFc76JYz8KsTGRbszSZIkSap1BtS6oLwYlvwGlj4AQTWMvSE8BCmxZbQ7kyRJkqQTxoAaTTXV8O5T8OZMOJQHAy+GiTMgtWu0O5MkSZKkE86AGi1Zc2H2LZCXCZ1HwxXPQPrIaHclSZIkSVFjQD3RcjPDA5C2zIPUbnDpE9D/AgcgSZIkSWr0DKgnSnEuvDUzfEtvfHOYMhNGfROaxEe7M0mSJEmqEwyota3iMLz9ICz5NVSVwajrYNyNkNQq2p1JkiRJUp1iQK0tNTXw/nPw5h1wMAf6nguTb4fWPaPdmSRJkiTVSQbU2rBtEcyeCnvWQMdhcNGj0G1stLuSJEmSpDrNgBpJ+ZthznTY+Bq0SA8H04GXQExMtDuTJEmSpDrPgBoJhwpgwd2w6jFokggTp8PJ34W4xGh3JkmSJEn1hgH1eFSWwYqHYeF9UFEMI66FM26C5LRodyZJkiRJ9Y4B9VgEAaz7O8y9FYp2QK8pMPkOSOsb7c4kSZIkqd4yoH5eO1fArJth10poNxCufgl6jo92V5IkSZJU7xlQj9b+beEd08yXILk9nP8ADL0SYmKj3ZkkSZIkNQgG1M9SXQVzZ8CKRyAUC+N+CqfcAPHJ0e5MkiRJkhoUA+pniYmFvPUw6DKYMBVadIx2R5IkSZLUIBlQP0soBFc+D7H+q5IkSZKk2hQT7QbqBcOpJEmSJNU6A6okSZIkqU4woEqSJEmS6gQDqiRJkiSpTjCgSpIkSZLqBAOqJEmSJKlOMKBKkiRJkuoEA6okSZIkqU4woEqSJEmS6gQDqiRJkiSpTjCgSpIkSZLqBAOqJEmSJKlOMKBKkiRJkuoEA6okSZIkqU4woEqSJEmS6gQDqiRJkiSpTjCgSpIkSZLqBAOqJEmSJKlOMKBK0v9v515DLSvrOI5/fzVeakwqbcIbaaiRQzalWFGkGGNB1FgyNCUVFGFg0MUKJxp6US8KUnrRjShnpESRUjKUBpFK0qi8OxdLK7PJockumErm6L8X69mxZzhnOvvMnnP2Ovv7gWGv/awLC75nNvs5a50lSZKkieAEVZIkSZI0EZygSpIkSZImghNUSZIkSdJEcIIqSZIkSZoIqarFPgeS/BX442Kfh8bqSOCRxT4JLQrbTy/bTye7Ty/bTy/bT6/9af+SqnrR/9toIiaoWnqS3FZVpy/2eWjh2X562X462X162X562X56LUR7b/GVJEmSJE0EJ6iSJEmSpIngBFUHyrcW+wS0aGw/vWw/new+vWw/vWw/vQ54e/8GVZIkSZI0EbyCKkmSJEmaCE5QNSdJLkuyK8mWobFXJvlFknuT/CjJ4UPrTm3rtrb1h+51vOuGj6XJNUr7JOcnuWvo3zNJVu11PNv3xIjtD0pyeRvfnmT9DMezfU+M2P7gJBvb+N1JzprheLbvgSTHJflJ+z+8NclH2/gLk9yY5P72+oKhfdYneSDJb5K8eYZj2r4HRm2f5Ii2/WNJvjrLMW3fA/NovzrJ7e0z//YkZ89wzP1q7wRVc7UJeMteY98GLq6qVwDXAp8CSLIM+B7w4apaCZwFPDXYKck7gccO/ClrTDYxx/ZVdUVVraqqVcB7gQer6q7BTrbvnU3MsT2wFjikjZ8GXJDk+MFOtu+dTcy9/YcA2vhq4JIk//t+Yfte2Q1cVFUvB14LXJjkFOBi4KaqOgm4qb2nrVsHrKT7efl6kmcPDmb7XhmpPfBvYAPwyZkOZvteGbX9I8Db2mf++4HvDh9sHO2doGpOqupm4O97Db8MuLkt3wic15bPAe6pqrvbvn+rqqcBkhwGfAL4wgE/aY3FiO2HvRu4cvDG9v0zYvsClrdfUD0H+A/wKNi+j0ZsfwrdlxeqahfwT+B0sH3fVNXOqrqjLf8L2A4cA6wBLm+bXQ6c25bXAFdV1ZNV9QfgAeAMsH3fjNq+qh6vqp/TTVT3YPt+mUf7O6vq4Ta+FTg0ySEwvvZOULU/tgBvb8trgePa8slAJdmc5I4knx7a5/PAJcATC3eaOgBmaz/sXQxNULH9UjFb++8DjwM7gYeAL1fVYIJj+6VhtvZ3A2uSLEtyAt0V9ME62/dUuwPiVcAvgRdX1U7ovswCK9pmxwB/GtptRxsD2/fWHNvvi+17ah7tzwPurKon2/uxtHeCqv3xAbrbAG4Hnkd3xQRgGfAG4Pz2+o4kb2p/i3hiVV27KGercZqtPQBJXgM8UVVb2nvbLx2ztT8DeBo4GjgBuCjJS22/pMzW/jK6icltwFeAW4Hdtu+vdhXkB8DHqurRfW06w1jZvr9GaD/b/rbvqVHbJ1kJfAm4oL0fW/tl+3sATa+quo/udl6SnAy8ta3aAfysqh5p624AXk13P/ppSR6k+9lbkeSnVXXWAp+69tM+2g+sY8+rp6/D9kvCPtq/B/hxVT0F7EpyC91tnkdg+yVhtvZVtRv4+GC7JLcC9wNnYvveSXIQ3ZfUK6rqmjb8lyRHVdXOJEcBu9r4Dva8g+ZY4GH8zO+lEdvPxvY9NGr7JMfSPYvgfVX1uzY8tvZeQdW8JVnRXp8FfBb4Zlu1GTg1yXPb36OdCWyrqm9U1dFVdTzdldXf+oHVT/toPxhbC1w1GLP90rGP9g8BZ6eznO5BC/fZfumYrX37rF/ellcDu6vKz/weShLgO8D2qrp0aNV1dA9Dob3+cGh8XZJD2u3dJwG/sn3/zKP9jGzfP6O2T/J84HpgfVXdMth4nO29gqo5SXIl3dN4j0yyA/gccFiSC9sm1wAbAarqH0kuBX5N9+CUG6rq+oU/a43DKO2bNwI7qur3C3qiGrsR23+tLW+hu+1vY1Xds7BnrHEZsf0KYHOSZ4A/0z3BW/30erp+9yYZPIH9M8AXgauTfJDul1FrAapqa5KrgW10TwK9cPBQRPXOSO0B2pWyw4GDk5wLnFNV2xb0rDUOo7b/CHAisCHJhjZ2TntI3likqsZ1LEmSJEmS5s1bfCVJkiRJE8EJqiRJkiRpIjhBlSRJkiRNBCeokiRJkqSJ4ARVkiRJkjQRnKBKkiRJkiaCE1RJkiRJ0kRwgipJkiRJmgj/BTVUiQTRzs9zAAAAAElFTkSuQmCC\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], + "source": [ + "plt.figure(figsize=(16,10))\n", + "plt.plot(data_year['year'], data_year['concentration'])\n", + "plt.plot(data_year['year'], [predict(year) for year in data_year['year'].dt.year])\n", + "plt.show()" + ] + }, + { + "cell_type": "code", + "execution_count": 12, + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "413.75316542940527" + ] + }, + "execution_count": 12, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "predict(2025)" + ] + }, { "cell_type": "code", "execution_count": null, -- 2.18.1