From 6d33769e012b8c4f50d5aa5dd90ecc6cd7219766 Mon Sep 17 00:00:00 2001 From: 0e57e1b8e1a23d7aaeee34c5821213f6 <0e57e1b8e1a23d7aaeee34c5821213f6@app-learninglab.inria.fr> Date: Sun, 23 May 2021 14:58:06 +0000 Subject: [PATCH] versiondef --- module3/exo2/exercice.ipynb | 263 ++++++++++++++++++++++++++++++++++-- 1 file changed, 254 insertions(+), 9 deletions(-) diff --git a/module3/exo2/exercice.ipynb b/module3/exo2/exercice.ipynb index 0932720..5bc1361 100644 --- a/module3/exo2/exercice.ipynb +++ b/module3/exo2/exercice.ipynb @@ -2250,9 +2250,32 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 51, "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 51, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYcAAAEKCAYAAAD5MJl4AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzsfXmYHUW59+8958xM9n0lCUkggRC2QEKIgAoEJAoKsnjD/S6i8l28iIrX64LeCy6XKFw/hYteUBQEuSogqKAIyr5IIAkkEMKWnWwkk32dycyc+v7oru7q6reqq8/pmcnM6d/znOecU11dVV1dVW+9a5EQAjly5MiRI4eKQmc3IEeOHDlyHHjIiUOOHDly5IghJw45cuTIkSOGnDjkyJEjR44YcuKQI0eOHDliyIlDjhw5cuSIIScOOXLkyJEjhpw45MiRI0eOGHLikCNHjhw5Yih1dgMqxZAhQ8S4ceM6uxk5cuTI0aXw8ssvbxZCDE3K12WJw7hx47BgwYLObkaOHDlydCkQ0WqXfLlYKUeOHDlyxJAThxw5cuTIEUNOHHLkyJEjRww5cciRI0eOHDHkxCFHjhw5csSQE4ccOXLkyBFDThxy5MiRI0cMOXHI0S2xp7kVf1y4rrObkSNHl0WXdYLLkcOGa/74On6/cB3GDu6F4w4e2NnNyZGjyyHnHHJ0S7y3swkAsHd/Wye3JEeOromcOOTo1hCis1uQI0fXRE4ccnRLEHV2C3Lk6NpIJA5E1IOI5hHRq0S0hIi+46d/m4jWEdEi//MR5Z5vENEyInqbiM5S0qcS0WL/2s1E3hQmogYiutdPf4mIxmX/qDly5MiRwxUunEMzgNOFEMcCmAJgFhHN8K/dKISY4n/+AgBENBnAbABHApgF4BYiKvr5bwVwOYCJ/meWn34ZgG1CiAkAbgRwQ/WPliNHjhw5KkUicRAedvt/6/yPTZJ7LoB7hBDNQoiVAJYBmE5EIwH0E0LMFUIIAL8CcJ5yz13+7/sBzJRcRY4cOXLk6Hg46RyIqEhEiwBsAvCYEOIl/9Lnieg1IrqDiKS94CgAa5Tb1/ppo/zfenrkHiFEK4AdAAYz7biciBYQ0YLGxkanB8yRI0eOHOnhRByEEG1CiCkARsPjAo6CJyI6FJ6oaQOAH/rZuR2/sKTb7tHbcZsQYpoQYtrQoYkHGeXIkRr79rfhr0ve6+xm5MjR6UhlrSSE2A7gaQCzhBAbfaJRBvBzANP9bGsBjFFuGw1gvZ8+mkmP3ENEJQD9AWxN9SQ5cmSAax58HZ+9+2W8vm5HZzclR45OhYu10lAiGuD/7gngDABv+ToEiY8DeN3//RCA2b4F0nh4iud5QogNAHYR0Qxfn/BJAA8q91zq/74QwJO+XiJHjg7Fu1v3AgB2N7d2ckty5OhcuITPGAngLt/iqADgPiHEn4nobiKaAk/8swrAZwFACLGEiO4D8AaAVgBXCiGkm+oVAO4E0BPAI/4HAG4HcDcRLYPHMczO4Nly5ICw2k5Y7su3JjlqHInEQQjxGoDjmPRLLPfMATCHSV8A4CgmvQnARUltyZGjvZGbyOXI4SH3kM7RLUHKMr+icTdm/vBpbN2zP/k+/7ZKOY4cOboLcuKQo9vjtmdXYHnjHicrpICo5LQhR40jJw45cjDIaUOOWkdOHHLkUJD75efI4SEnDjlydBLunf8u7n5xdWc3I0cOFvlJcDm6PSoxS+0IU9avP7AYAHDJjLHtX1mOHCmRcw45cijIxUo5cnjIiUOOHAxyU9YctY6cOOTo1lDFQy5MAeVucDlyAMiJQ45uilw8lCNHdciJQ44cDPLYSjlqHTlx6CJ4b0cT3tm4q7Ob0eUgkE5/kHMcOXJ4yE1ZuwhmfP8JAMCq68/u5JbUBnLGIUetI+cccuTIkSNHDDlxyFEzyEVGOXK4IycOOWoGaZTM+UGE3QsPLlqHlZv3dHYzuhRy4pAjhwLy2YucNHQf7NjbgqvuWYRL75jX2U3pUsiJQ45uDZUDcBErBVly6tBt0FouAwjPB8/hhkTiQEQ9iGgeEb1KREuI6Dt++iAieoyIlvrfA5V7vkFEy4jobSI6S0mfSkSL/Ws3k79NI6IGIrrXT3+JiMZl/6g5aglUoYKhu+sl3tywE+OufhgvrdjS2U3pMOR0vjK4cA7NAE4XQhwLYAqAWUQ0A8DVAJ4QQkwE8IT/H0Q0GcBsAEcCmAXgFiIq+mXdCuByABP9zyw//TIA24QQEwDcCOCGDJ4tU2zc2YQ/LlzX2c3IkaMq/H3ZZgDAX5ds7OSWdBxy9VFlSCQOwsNu/2+d/xEAzgVwl59+F4Dz/N/nArhHCNEshFgJYBmA6UQ0EkA/IcRc4fH6v9LukWXdD2AmVbr1aydccvtL+NK9i7CzqaWzm5KjA9DdA+8dWLOrfdHd32V7wUnnQERFIloEYBOAx4QQLwEYLoTYAAD+9zA/+ygAa5Tb1/ppo/zfenrkHiFEK4AdAAZX8kDthfd2NAEARLmTG5KjXdHd18x8F53DFU7EQQjRJoSYAmA0PC7gKEt2bn4JS7rtnmjBRJcT0QIiWtDY2JjU7Bw5vPAZB+hhP52J7k4EI+jm77K9kMpaSQixHcDT8HQFG31REfzvTX62tQDGKLeNBrDeTx/NpEfuIaISgP4AtjL13yaEmCaEmDZ06NA0Tc9xAGNPcyu27tnfPoWnXBgOMGlmarSVBf64cB3KZf7Ba1HEUntPnA1crJWGEtEA/3dPAGcAeAvAQwAu9bNdCuBB//dDAGb7Fkjj4Sme5/mip11ENMPXJ3xSu0eWdSGAJ8UB6oX08OINeHBRrpjOEjN/+AyO/8/HnPMvWb8DF976Appa2hLzVroYHpijLxl3vbAKX7p3Ee5bsIa9Lp+ri9PAVOiq77Kz4RJ4bySAu3yLowKA+4QQfyaiuQDuI6LLALwL4CIAEEIsIaL7ALwBoBXAlUIIOYuvAHAngJ4AHvE/AHA7gLuJaBk8jmF2Fg/XHvjmH7xzf8+dMiohpxsWvrsNx44egEKhhmarhvd2NqXK/50/vYEFq7dh0ZrtmHGIXTUVPewnuY+76ltYvWUPCBRwYI27mq35uzqHlAa1yC1lgUTiIIR4DcBxTPoWADMN98wBMIdJXwAgpq8QQjTBJy61hL8v24z/84uX8B9nH4H/+/5DOrs53RIvr94W/O7Oi8QHf/A0AOCqmRMBAK1GsZKH2iENOedQKXIP6U7Euu37AABvvZef05A1pFTylqeXV3Z/lo3pQJR8DlR6BesQNUgduuq77GzkxKETIednOd/aZI59++P6CCexUhdfNEtFb0qbOAeJ/KzsHEnIiUMnIpD75rShItho6l6GOKQru2u+FMk5tLW5WSu9u2Uvrn7gNbS2dV8Hnq76LjsbOXHoRBRy2tBu2OdgycSja++oi4FYKT6qtu7Zj627PYW13Jdcde9C3DN/DV5du73D2tjRyGlDZciPCe1EyAmai5Wyx7C+DVi5eQ8KVFvEt1Q06xw4c2FJQwpdXZ6WI3PknEMnQk7IBPFwzhZrcFnGjh0zAABw2Snj05Xdxbm5UsHXORjEShKyD+XY6s7EIZ8+lSEnDp0ICoiDffTu3NfaEc3pVqiUoIaLZnZt6UhInUNLEnHQuNZuTRy6LKnvXOTEoRMR+L3lYzdzyMU9ssinWv+65kuRzpSuxFFKn7oxbchRIXLi0ImQ5oS5zqEy2HaEQvt2RVdfJF2NHPSx1605h3x6VYScOHQiAvl2Pngzh+zT9jLR/NhPnsf/vWt+JO25pY3YvNsetqK9UXAUVepjr9CNV4J8elWGbjwkDnyEuzz78M1lpukh++yuuauxrYKIr0kE+7W1O/D4m5uC/0IIXHL7PMy+7cXUdbng5dXb8PtX1rLXLrj1heB3qEtwK7c2OId8/lSC3JS1U+FmreQ60WsFLuuYGrL6qbc3WXJqZVfo5yDXn2WbdtszVghJAM4/fnTsmho/6v6X1/rtcbNWColDBo08QJFPn8qQcw6diIBzSJjItaiTeGnFFvzp1fXJGQ1oUYiDjbgKIfDUW5vQVhYQQuDRJe956RXX3LmYt9I7BsV1yMi+6c5RWmtw+mSCnDh0Ilz9HGqROPzDbS/iC79dWPH9pvAROh5/cxM+fed8/Py5FZi7fEvF9VXzhnY1teArv3s1k/PJXR0r97cJtLaVa2Rs1cIzZo+cOHQipBIwiXPoLvP3lXe34e4XV2dXoKVfWgxRSXVs3eMpkJdv2o21fpTcippSxUu64/lVuP/ltfj5sysqLkPC1QLup88sxyd+NjfI113GWI7skOscOhHkqHNo6yZKh/Nv8eTml8wY2+51ufZZfcmj0M2tZezYG+7cO2OxzFKw49L+V95V4yl1jzHGQe2LJ97ciIZSEadMHNJ5DeoiyIlDZ8LRJr02WP8KYFlNk8JHSDSUigCA/a1l7K/C7LWaNxRYo2Uo90+7n+jOQ0x9tMvuWgAAWHX92Z3TmC6EXKzUiZA6h1oRK2UOS7+YDrvRUe+ff7C/rRyxcEprPlzNO1rvi7OyIA2y3WnFXN15iOXzpzIkEgciGkNETxHRm0S0hIiu8tO/TUTriGiR//mIcs83iGgZEb1NRGcp6VOJaLF/7WbyTSSIqIGI7vXTXyKicdk/arbIwnbaNY5PzjmkB8c5cIuvFCvtby2jrYp+rsYX5b4FvP9CJZD0LW1rshpi72zcxR601JnI/YQqgwvn0Arg34QQRwCYAeBKIprsX7tRCDHF//wFAPxrswEcCWAWgFuIqOjnvxXA5QAm+p9ZfvplALYJISYAuBHADdU/Wvsii8nk6s3aTVQOHYoWx06T76C5tS3KOXSGzsHAOqTaiPhZ024oslhAm1ra8KEbn8Xnf/NK1WVliXxvVRkSiYMQYoMQ4hX/9y4AbwIYZbnlXAD3CCGahRArASwDMJ2IRgLoJ4SYK7zR/isA5yn33OX/vh/ATDrADa+zGG8FR7ND9Xp39vbM8tnaHMVKklvY31quighn0XSTA14lBglp25NF+1t8nc2LK8wmwXe/uBqfvXtB9ZWlQDeeMu2KVDoHX9xzHICX/KTPE9FrRHQHEQ3000YBWKPcttZPG+X/1tMj9wghWgHsADA4Tds6GpksZI6xlYTo3B1tRyHLZ3NRSL+8ehvm+05jzZpY6UDq5jS0QXIAqTmHLIaz1KFZ8lzzx9fx1yUbq68sR7vDmTgQUR8ADwD4khBiJzwR0aEApgDYAOCHMitzu7Ck2+7R23A5ES0gogWNjY2uTW8XZLF4uNqkq4vDgbBoPbe0sSqHMRNcn80lxIXLYnfBrS/gJ08tAyA5h87l0Ey8cpqFXmZ9bulmvLx6awatcseBehZGrnOoDE7EgYjq4BGGXwshfg8AQoiNQog2IUQZwM8BTPezrwUwRrl9NID1fvpoJj1yDxGVAPQHEBvZQojbhBDThBDThg4d6vaE7YQsJ0DSwSydvWjpuOT2ebj459kHmMtS8Z52QVixeQ+aW6owZc2OkYwhjVhJ7cMLbp3rfF8WC6gIvpPLunf+u9i0s6nqOl1wAEyZLgkXayUCcDuAN4UQP1LSRyrZPg7gdf/3QwBm+xZI4+EpnucJITYA2EVEM/wyPwngQeWeS/3fFwJ4UhwIq6AFWe5G9rfaFyVVfN6RndLU0oZP/XIeHnujY8QAWb7xSsq684VVFd+fxXgwcQ7VWFG5Iosq5JR1KevrDyzGlR2guP75syvwqV/OT86YIwYXJ7iTAVwCYDERLfLTvgngYiKaAm+9WgXgswAghFhCRPcBeAOepdOVQghp23YFgDsB9ATwiP8BPOJzNxEtg8cxzK7usdKjcVczlm7ahZMOdfOczGQy+QtKS4LzVbmTdA7v7WjC0283Yt22fThz8vB2ry+9hU1l19zK7gyxEk8dRAqGpjN3VAHn4NiIXU3tf/ztnL+82e51dFckEgchxPPgOd6/WO6ZA2AOk74AwFFMehOAi5La0p645PaX8NZ7u7D8ex9BsYPjFyd55kaIQwdOf1nTpl0dc4BNtpxDdYV1hrWPCWk4h0rbkQ3n4H87jtHeDfblZ9/+Nqzbvg8ThvWptmk5KkDuIe1j1ZY9AICd+9wiY2ayGPhltCSJlVSFdAduDcOgbMmVPvNOIzZWKUPOkvBxJaWxjk4vVmo/dERsrUz6XhIHx6J61Ret17/w24U440fPoKnlwHKqqxXkxMFH3x51AIBte91ODctSgbc/lUK66mqdkaauS++Yh4/95Pmq6nNdA13W+Gr7Kb2HcfvpHDrCQz6LKoLNhGP+pNPn5i7fDABozb1AOwU5cfAhY+w0K7v4M3/0DE6Y8zibP1trJTvnIDpLrJRysm/cWZ34KUsbhKp1Dp1hypqhE1xaZMgIZ9Z3knh0l6jEXQ15VFYfnKXFUsuRj5lMJilWStQ5xO/pCAT1ZlCny6luB9IakJpzaJdWeOgYziEDTjjlZiIpX8HX/ZUPpIFRQ8g5Bx9pbLSB6ifTvJVbg2icyaasKufQMdjd3Iqzbno2s/JcTnVLHUnUaq7UsXKlbDyM+XTHSCBVwaX5jyzegEaLcUJaa6Wk9y3tQlwPbsqRLXLOwYftRKz3djTFLIqqXQs+8bPQQSlJphrlHDqGPKz2FfSA/VnXbtuLEf16ZFJnptZKVd/fGWIlHh3h55CEnU0tuOLXr+DIg/rh4S++n82TdTOl1WAuVuoc5JyDD9v4m/H9J9Ckec92rGK44zkHU/0qGnc145QbnsJ1D2djS56ph3S1jEMnyJWMTnAdoXNIqEKeyb3OcpRqlgR1V1MLNu/2jENscbL+7b5X8bNnlmdWb44QOXHwEdhoZyUwNWB/axnn3/L3VPd0hs5BrcdUpbTseuaddHGu1m7by9fpeH+StdLb7+3C4nU7UrWp0raE+TOwVjLwDh1znoerODX9tT3NrdjuaAUocc6PQ8s3G2f9wCtr8f1H3kpVdg435MTBR6hMc5wkFS4Gq7fs0c7uTUbZZaVuR5gmvWvIcR0X/ZSP+eNaji1buSyMuhITTeGcrA6k8xz+7KDMrxZZOsHpOPX/PY0p330s8f7Zt83Fid/zrANXbwk3EK7h13Nki5w4+JDj2pWDr3Qy7a3glCx10TywToWrTCZsVGpm8GiVWJgVmVW5I48JTcLNTy5rv8J9mJq/fe/+iBOajWsz9Znpfet99uKKraw5dFJgSgA45YYnE/NE6z6Q5tGBiZw4+EjjDQxUvo5VSxw6RedgqFUuFFnNsyxE65VEPuEWvAPJQ7ojYHreKd99DP9w24uJ+YD0sZJcCbDL5mPttn2pFvxcx52MnDj4COPCuOavbHS5HnyvIhKVtYN2PBFJlqFKuaZmpTDNgiuyxcUyXeHuSW9Wm4HOoRMPP7S1/9U1bmLQL92zKDlTBUjyA5JI40mdW0AlIycOPmymrBz+8vp77diaKDqDc3DZ1cnFLCtTyywU0pUETSxxxCF1KdXDpeVb9+xvF6ewLPp+lWL+7FSnqwjXsTyX0/8kuI3IzqYWfP+RNxP9jmoFOXHwEY4VtwFWqZKwknW0szykgzoN6eHJXyl32Yb0LBa9SohDgeUc0pWRxWvhxVvRko//z8fwo8feyaA2vR779bfe25WYr70i2bqOrzRcOUccfvjXt/GzZ1bgj4vWOZeThLXb9nbZwIE1Txweff09HPudvwUvsL0X30qK76zYSkqlhrZ4OJBY9IqIA6eQztJbuwqoBxBJ/O2NjuNaJbiT/15evQ1v+0QDaD9jCddi04xDLq90dE3DgSThlBuewj//akFm5XUkap44XPfwG9ixryWQV7b3MlfJBIqM407xcxBsunyWtLTBtHynP+wnnr+S9Ym3Vup4cP24YNW2Dqm7kk3HBbe+EDEbbi/fENdyXayaJGxMxp9eXY9HX9/gXFYSnlu6ObOyOhI1Txz0daHdg3xVJFbqWJ1Dc2sb/vuJpey1CJ2SxEHps/UWD1quDBVZdH0lRfAK6bT1prvhjfU7mTqjZTS3tuHhxfFFyuQsVxUcm281ZW0nzsF1TqbiHCxtnbtiC/7lf6s/wrSrm8vmxEGbaO39OivZobme59DU0oZvP7QE465+GPfOf9fcBiHw1nvxxUniVy+sxpNvbQr+qzsydcDLubirOTRh3LI7nSes3i6XtEphWthY4pC28JQ3nH9r3Etef9RfPLcybSsqhruVnvlae+2r3DmH6nQObiYB7ujitCEnDvqC0d4vtBJnT06Uw+GOv68MZNS/mbfGmO+e+Wsw66bn8Kwh7MU+iwJNNRfk2uJijWkWK8XTuMe17ZwrISa8Qrr9BsLc5VtisbqAeH/+4K9vs/e3h8VrFo+rtv8pZXNRbZ2mMa9zFFv3uG9MOG4k637t4rQhmTgQ0RgieoqI3iSiJUR0lZ8+iIgeI6Kl/vdA5Z5vENEyInqbiM5S0qcS0WL/2s3k20ISUQMR3eunv0RE47J/VMPzaf9dd/aVDqRKBoyrWEldcGzNkyKNNKaHF976AgDg/ykLVqVRDczPEL+SWg9RQQdzpqw2rNkajw3lWm1zaxur3AU61zEri3Gv9v0bG7wx5hLiOxGGjLpo6NN3znctsUMi3R5Y0QzSw4VzaAXwb0KIIwDMAHAlEU0GcDWAJ4QQEwE84f+Hf202gCMBzAJwCxHJw2JvBXA5gIn+Z5affhmAbUKICQBuBHBDBs/mhJjjkavstUIWtJIdqWvIbnWNq2YXxN26YLWnGF3WGIanyHrws5xDpjXwaCjFp4Ht0c73CaVrfhU2ubj7Ap096+BuVupe5qI1240nKaaBqcv0vkzDOXDvIete7fbEQQixQQjxiv97F4A3AYwCcC6Au/xsdwE4z/99LoB7hBDNQoiVAJYBmE5EIwH0E0LMFd4K9yvtHlnW/QBmUnvMAAZxzqF9Ucnu0FXnoO6A26vzkjynXWz1Xcq2pVWa34S6IkMcLCOB2w27HxLldq2jlZlZ11YgwtsWvZatUl13YOrbNDoGHR1ygFLXpg3pdA6+uOc4AC8BGC6E2AB4BATAMD/bKACqwHutnzbK/62nR+4RQrQC2AFgMFP/5US0gIgWNDamCxPtivZ/oRVwDo4URZWdtxdtTQoCyHFU1ciWO2L3xUmVbNX2baj8jCzb06gE4V1GdCWR5s26E+ZsjQFchp9p0f/h36JOfqZmLG9M55GtolJ9WRrUDHEgoj4AHgDwJSGEbUvAdbGwpNvuiSYIcZsQYpoQYtrQoUOTmuwG3ZS1nd9oe3pIq/b6HcE5cH3F+h7EyjAoGDPo+6ycBK2lVBGoz7bgqu+ZU1hXgmpCVFTzOgpUuej1p9rhPaZx0VyF53Guc0iGE3Egojp4hOHXQojf+8kbfVER/G9pnrAWwBjl9tEA1vvpo5n0yD1EVALQH8DWtA9TCSoVKyXtMsZd/TD+889vxNKrFitZWliJdzAH17DM3LNwstysxUquSlFXcFyWrRynnYwBtnzqe7bF90mzw01aoK49Z3LF99rgQhiqja1UjQKftVbKeEvV7YmDL/u/HcCbQogfKZceAnCp//tSAA8q6bN9C6Tx8BTP83zR0y4imuGX+UntHlnWhQCeFO0odL3g1hfwP095MfL1hSGLamUZtz+/Mn6tgp1tJGSF5XY1DEQlLLIQIvH5VVktJ/NliYNz/Uxa2nMVUuX2wC/27TP8XHUO+pnlKtK826SnCPYTTMZqFl8X3ZP7uOBzpn1H3z33yOA3xzm0pymrGmakq8CFczgZwCUATieiRf7nIwCuB3AmES0FcKb/H0KIJQDuA/AGgEcBXCmEkPzfFQB+AU9JvRzAI3767QAGE9EyAF+Gb/nUHhBC4OXV24w25FksCbbQwZXQHnXBtd2uDm7TLujFFVtw94ur2Wsf+MFTeN/37YemqLuhbz+0xHpdQk8ydQ832W0L1PVZHQ+ZUkzEcxrhDZt2NmHc1Q/jgZfXxvLZXqBahk1kkmaHm7R7lXoqvu+r0TnE2/jEm8m+DxxMzUjbvHrF8KBDzuVW6PtZNz2LXU0t7V5nlkjUrAkhnodZhD3TcM8cAHOY9AUAjmLSmwBclNSWLLBzX/RAkn364TsJY2bq2IF4ebU53o0QAq9Yrlcy4VRi4yqz5t5Yc2sbZt/G29gDwJqtDqEvlDreYnZD3IbXdYfHO8GZiY20pTflP+OI4Xj8zY2J9fJKdHObkw4HenSJFxjvt/PexQVTR0fzWfpCfX6bI2IaJA03yW1mZfkVlhtP25byHOmkdqQmDorJckdIfPS53tVCgdech7Qe1len5rbJe9jwPvjk+8YCACaP7Bekl8sCrf6qeMffV0VOznKFbTGKeiW7lcFR82bHwWmzdEoibrzOwalag4Lbnm/c1Q9jwaqtbP7+Pevw37OnBP9TiWNS6hxUXPugx1FxYjerWElpve3EwDTPkWT7byMO1XAOXKRbPTCeqwjX6CGdsn11CZyDa7fe//JanPuT5xPz6TV88Z6FjjUcGKg54qC/MP2/zf65QIQPTR4BABjStyFI//xvX8GEf/ckZFxAtUh9FeyCogPZtvNUiAMz0vVQxJWIWJOmYzXsOrt7Zd6HXseFP53LltHc6rb7ZjkB7f9jb2zECXMeR3Nrm7OZMBcl1K6QDn/HOFoFad7bBYzDngq5m+baVc3mmuui/Y7vw7UdadtXVwwbVY210ld+9ypeXbsjMZ9OvP6+bEvFdXYGao84JFAH25ApEAWDXn3xf1kcxtdPEqGYdju2u1ojge9sZYe/OVFJNU5DYR0JnIODzsGEPy5chz8sXBtZGNOH5Q4v7tjnJuNlFdJaHd/50xI07mrGpp3NTvkB/vAZ225ZvbR3v+U85hSsw4YdTdbr0juca1d11kpx6Er2ahXSadtXLITLHWetlFaUt2W3OTQIUAPWSt0N+mIToxWWF1oskJUNZwuMlc+n2wZSm7LI2IpXFbT8zi0L4pBwnbVWcpskd7+4Gv9676u47uHQBJjrFtdd366mVqddPs85uNUhCRmXv7VN4K4XVuErv3tVKdcMdezZrJWyhJVzqKIJRBSjEHGxkltZaeaabf6qEQQ4DnfttmSdm4qp1yWEBunatKG6xC3dAAAgAElEQVT2iEOMU0hhXlegUNFm8lpOGg/Gca5d+MFf38KVv/FiykcV0vz9O/ZGd8kscchiwUkydWV3oOmq2Lgz3JGl9ZqOmIO2lp1EMC5e3ep/tW837DAvKPvbyvjWQ0twv2K1ZOf83DjELC0uQ87B3h4bZv7wadyhmW1nuTkxjR+ufTaxphpBoGOc4Nq9inZFzRGHJJ2D1Y9A4RyM5pgVLJ5c+v88tRwPv+Yd9BI1ZeXv/9BNz1jrBTqGc+Aeb/7KdP6MqqULr5C21K/8TnOmsK2cBxetw7rIIUZhA3vWezElWbESq3NwI2z6M86cNAztgdCCp3Kx0vLGPfiu5vDJKaT1A6ScxUqGnNVwlVy2rJxIJXKxUhdD3OY+iVyEMOkc3O62Z7DF0nHhHNTdNsDvhrMgDpXoVO6xHDzEQV1Y2B2tox9JS5uI7GC/9WDcLwPwdrnvOyQayksS+T3NrbjqnkXG+mQ/cy1iiZOjzkjv58NG9I20l8NHjh6Bgb3qzBUw4CLScu1JC6dltoqFHDBwlUqX6xs1NT9rrZSxF1zXJg21SBy0V6bPX9t4LRKBfAJhVpKlq1/icssh5K2KOMh1M8KNc9eQxtWEp0i70+eg6A2rUpS2tJUjRHJnE6/kJQKu0cJISIL8LcbRb7OiiLQRS44Y21puEyupBNH0egiE/559HHvt1TXb2XRp3sn7OTCcj2PfMyqHipFmI6ZyDjHTbcHnk9AZh2qjJbT7kcPtjNojDtr70geJ7X3KRbNAVPGuyhib3tHPIU3IgKaWNlz/yFuB5UtjgnVFUIfj7pa/t/oJoe7gXIgNRcRQ4UV366z4MiZ3lht3Rq19YmeO+9Vxz727OU6MbN2z3eIkpu50TTtcAYEPHDY0xgUBYB03P33yOKNCuqWtzId/cd6ctH/Efa7P1X7aafFhclm4q/WidvUrOlBReezhLgr9desDwLb4hsTBvJgnLY6my/17msUBbQ5iJR1EhP99cTV++sxy1BcJowf2wi9fWOV2swWVPF+yl250wd+mcDguO1qTGEoXK6WBJCxJXIpsiyvHZBtff1y0Hpe8byymjh0Ue8Y0z+F6JvZ/nD0ZyzZ5hzfpj3nXC6vws2dXxO5x5dpcmvvq2h0QQiQSEiPnkCBy1KMhqFICF2dNE23wJAf8NYn5q7biIsX/piuiBjmHyhZvIAwvUSoU2MH1pXsWJoYIMC0O/XrwxGHz7mYnnYOO+iIF1kn72wS+9sBreJMJN5EWSYsDP6Hs95QK0WGoLm4u1koRBbZKHFytlSj+XuT71d+zHmIkLaOUlH+J70Sp92PB4SAnWTa31haZtAKFefXnN4ngXB/XFl9MhYsFnVnnYK+3SfNbiJ5FwtQTK78yzmFnUwvu/Puqiu49kFCDxMF+/Q8L1xmvSeVdQ10hNvAAb+eX5AVpmjMm4jDtuseDBQNwFyuVCgXrYmGDieXetKvJqjgH+PYl9bm+0y0kiJX0BdsUiK6lXHZ6dtLaOKRPQ2CTr+uk9POf5X3ONvva/8+cPD7y3+RHU0zxEjlLIY6bIAp7Lm6XESZcNXNi8Nt1wRzQs65q0dIPLzrWr5O/njZYoEjIp28cTWIlNVurRtz27W/DMd/+Gx5evMHYjq6CmiMOKiZd80gs7cm3zJEjpfKuR6mI5koPYzEM3kOG9jbesmxTGODOdRGqs1ihSNz85DI23SQymz7nicRDaLj5lGSPri9c6n9b4L0ARp2DowhEW8PqihRM+iTTyPDZ4vlGDegZS9Of59qPRhXhRUOU1GLklD9rk9igdwWDmWbIOZhxxMi+gSltEkNwyBDzOOZg697JB/Xz8/CZks4TifuqCDafCS6E8EM3PRshEB3lvNgRqDnioL7vtKdtyYnUUFdAU8ZxYnr59vKV3q9DjSNjWksadzXjCSZqaTVWFuxizuS76fHwKEgr5+BAbNS7f/ZMKCf3FgAHD2ktT6lIRrGSjoA0MNnSBt4DwoVdzxc5qyPhmVjOwUhRJKdiblipUMDw/j38hlmrhplUGvJbMpYCQmm6177pkAT2Ox87Eg9/8ZQId+4SWt62UZRY0bgHP39upVJI4i1dBrVHHFK8vRH9ekTv9W9tKBUq5hyMntWWZqly1N8tWGPOqKCu4PZq75q7OpZWjZGGq0JaDVCob2rVMNsujk7qYvjQq+sj11wlG2qRpUIBLX4nJOtYzNc50WMSpChGLzVCQBOeiXvmJM7BhlIxJEdJ/WEUUxlgKy/gotIopJXEW5/2jhsdM6gnjjyoP9536GD825mHATAdShVN+4W66FuwdY/do7+rovaIQ4p3F4/D5P3vUVfMnHOwNUtt869fijuUSYsTFaUiBZMqrei3mtACLmdI67DJplm5skaXq3Vs1RXSpQLFZMkmBDoH5hqnlE3Uv1DIOhAB8745E3O/cTorojKVzfWniXNw6bq6YiHR+TMoz6DgNsFWnjRUSOUE5ye9umY7HnndC4gpSRsR4dwpoyL5VOjFVRKosvuQhlokDiny6gNIDp4CUcU20KaBXs2G44wfxUNnqLHr7e3hFt/KG8M6BSeImmyLu4tYaU+V4a3VPP16lFAqFoKF3dmU1ZFjSoLcLZeF165h/XpgZP+eOP/4UWx756+KhybhCEESI5kk3tE9wU27eTIo1E2wDbViMSlUDVeel7hF2c2rHSb7gQ0QqSWVOBOvhHbknEMXRpIp69Gj+it5tXv972p2qmZvzzDdFq65znnAVj5I/6SJZtKAq5V75qg81/xMroH3khYrG4gIPes8nc/7Dh2MYiEUO3DxkSL1Wq7xhxe577zVtpueg7Ol1wmBEII9tU8t1xatuKRwDknWWYHOhL8ch404BHHM+Exq+uHDvfAictFXb1F7ThJfjjvW63HdYNnK6MpIfHoiuoOINhHR60rat4lonXamtLz2DSJaRkRvE9FZSvpUIlrsX7uZ/FFJRA1EdK+f/hIRjcv2EaOwvbp+PUroUad2STS3VJARUeaDQC3OFuaiT4Ob32JZhGX+z1PLHeoPG7A+4QwAe73pRSnWcB1M2sGDejH1xvNNGTPAOYTDxOF98eOLj8MPLjoWBQpFckkconxefXE96dDBrIesq1mvEJYwGYYLsg06Ibl3/pqIoj5SlkO7GkqFWFtM2dOcbw0kiJWKPHHYsrsZjbuagzY8+qX34wszJ/h54+WoOilJcHidQxRjmXH2VJKSuvvQBifO4U4As5j0G4UQU/zPXwCAiCYDmA3gSP+eW4hImuHcCuByABP9jyzzMgDbhBATANwI4IYKn8UJtklQLFBw/bI752Pz7ugifeiwPgC83ZHMlySb1ne0rxlOkFJz2dp46UnjrPWF5blJfV3s9E+fNAxjB8cniqHi1LAtJ1y7jh87MJamLyBXf3gSfvPPJ6Zqx0ePPQj9etT5xN9LS7RWMvRfbwMRT+qewM8BZiKQ1lppieV0Qp0j4NCzvhgjOKZFPSxPOJEJG3GQRhW6SfLU6x7HCXMeD+ZW3x51wcYtINYq5xARK1m4ES3pj4vWx0Lhf/rO+bHbnlAIRhcPpxRBInEQQjwLwDXm8rkA7hFCNAshVgJYBmA6EY0E0E8IMVd4b/RXAM5T7rnL/30/gJmkj8RMYX57xQIFV59gdgg9Sh6dI4Scwy1P23fl+sR828DeuwppXVndtIPUlr2kEM3ketNzDjaw4iIHPcSIfj3Qq77k5gSn5fHCefhiJUfOYZfmUVwyyB6557nhgqOVulXOIZ2FUaCQ1tJtBgamqLLqLb3r44TOTBzsHa6LRW3dKzkH0wYs1AGG9fJnQ4d1BiH3HQ+l+uULK80N9LFy857gd02JlSz4PBG95oud5FZuFADV1nKtnzbK/62nR+4RQrQC2AEgHjksIyRzDtEM15wzGf/+kSMAhANbja2yShkYLmgplzFz0jB8ysIB2NroqiwWQkRMQo35/AmRZFLoan3C5RrY2x5G2vV8Bmu9Mf2QFLEk36svwp5Yyfvtyjl84mdR2b/pbACutH844eDYfZ7OwVp1DJKQ6VXbxoy60zfB4xyiaabsSWIqPWqs1b9CEgdD+2UygQJxEVec2vZArMTtOZg0VzFuUEaq3Ac2KiUOtwI4FMAUABsA/NBP54azsKTb7omBiC4nogVEtKCxsTFdi20F+ygSxa4P7l0fRK4sMcShoc7uvKaz+EJ4IoeJw/sY22U9MtRxZyKEWYSl51O/ORSq5BzGDbZ7zbqeqxykOeRLs4GzcQ6HK+couLYPsHEO9raoTnBGsZIhXSr5dZ8Gl92snkOto2ddMUZAjcQhwZRVJ5o22huKlQycA0JiGFghMQ2LipX8fFrFbWWBBUzk2r49eOIwffwgNr2rh+lWURFxEEJsFEK0CSHKAH4OYLp/aS2AMUrW0QDW++mjmfTIPURUAtAfBjGWEOI2IcQ0IcS0oUOHVtL01Itgz/piMDilSIcQ7qSjCmyuzOj/shAeG2yZbDZRhuvYS8ve2jiDIrkTBy6fewiKdNdc8qVVkAJRg4PRA+3+BaZ+MzmdJe0tVccvo1gp4ZmOPKh/5L9LGHq9WWp3chZyxudWxWJMM3UzW9u4KxQIxQIZicMfZRw0RawUGgiEUPvLZK1099xVbB09lM2f2o4j/dAeOrqRVKky4uDrECQ+DkBaMj0EYLZvgTQenuJ5nhBiA4BdRDTD1yd8EsCDyj2X+r8vBPCkyOJQAAOsi6Cvc1A9W3vWFTGgVz0A4KD+3kJRKIQTrqFk5xz0J/GIA8UPFlHaZRJleIpw18XSKVtYf0pxm7mceL6k3ZTtKss5OJghBn+dxEpRqCHEkx7b9Ggmp7Ok8uTlSjgHic+cPC4itrQRWM6U9fV1O/Dgq+sieWKiqiTOwdQvKTgHwJsLJmu7+au8nb43n8yK5gjnYNBN6MYnXHv3Kv409YbYZWkiMBzoSBSoEdFvAZwKYAgRrQXwLQCnEtEUeGN5FYDPAoAQYgkR3QfgDQCtAK4UQsgevQKe5VNPAI/4HwC4HcDdRLQMHscwO4sHM8G6CPryolk3PRuk9aov4vzjRqGuSDjnmIMA+JyD42KpD8JyGcFpcqZ2mXZKaZzvqtnpc/W6Ehsum82gq6mlDds1i5BoeZURwzRTVH8X3vstO9VvGgcmByqZ++uzJmHmEcy50CL8MtEALnaSCiLCURZ/nUheJu2cHz/PlBn9n6SQNrU/JvJyHFjNrW3GjRhB9Ynw0tT3otYYcjbRetVmPf/103DKDU9F8uv3NBgMQ7qRVCmZOAghLmaSb7fknwNgDpO+AMBRTHoTgIuS2pEVEsVKAFZtCcNS96grolAI3e4BGW5BlmcfDfpCL6RYKcZehzByDoUUi3RK4mXnHFIs0kwD2zi3aR/3zrfHiuI5By5N5xy8/y5CJf1dFAqA3NJUzDkk6BwOHtQLhw2P6zNUU0yT5Y9JYmU6MGreymRjw7S8uil/qDPhWZ86rfGu9e5pthAHhbPhxp9KkKQuaL+mkVb7evTA0GzbdFaIiXPIrZW6MNLK1nsy0VJVO/ikwRDjHESUDQ7a5aBzSCdWcsvXUk4+8azA9IsJXDYuuF+QP6Fgd0U4345KrKILis4huXo+R0srT6ySrKgCAwGLtZKJIH3ZDyoHRIniuu374pllvkCBbEd8vBo4B3ndUI6+qO7Yx3ONf/jcSZH/r68zG1c0lAqhKSsnVlJ+FwqEvg0l7NTqNYvwFM5BSTeZlLejRLzDUXvEIZFziGbgQmkXKCwoaSevL/RlIVAoAG+/Z3ZM+vQv5/HtSyNWcsoVhoew5SeLWEmP36/nS4pM+uiS96zXF2uLQlNLm3OIDiB9bCUg+ryVcg5Pvh31k9GtwpLU1UKY85ieNXJanCNNDPwckgaMoylrQaE2XBt04vC1B15jyznu4Kij46otZpPxHnXFiGd5Evr1rIsRB5OoTk1V+71kJA7J9UfzC9y3YI1TBN9KovxWg5ojDjZwzl49GVNVAqriHIgoftyksuSZjmgsppD9u+ZrDcRK5hu8eRO9/sHDeGsxvZyv/O5Va/0vrrCLPL75h8XB7827mzHpmkfZg+9j7yGV0iH6N8qhVcbZ6App3XPXzDmE78PE9aQ1n7UhbRRVCdO4P93XowgItp36jnvDDjNXo6LBIMapK3oWTYFYieH49Gb07VHC7uboHHPpOfV5TF3NzbsrTzvUWObT7zTia/e/husfecta90+fWY5J1zyK7QnHEGeJmiMOaU1ZezDEoUAhh5G0UzDpHNTJ2FAqOC1mRO7iItd80vvUllsNFyJx9tG+wVrCjvLFFfZjU5Ogyu43bLfEfNLqlc9fqRNcEudwenAyGp9B1zmEpCYQeEWuf/Wsw7V88ba/9M2ZGNS73uxDQJVwDn59CcNF31lzi+Bphw/Feb5uzlSezjkkKdclTEYNcn7qHtKuVmFBOxycFtWNz1sb+EgH3Hg4ZEgfJqcH6Vm/eXezMQ8A/PfjSwEgRtTaE7VHHKw6h/ig4XYsRGFo6pcZxxkVelRPacqqTi5i6uXgKaSrEP4zkMHh7BYt8XpNkyl2SluVkVBUzs327oymlQ57wueWRh0qVSc4U78kLaom6x4T53Da4f6OW1VIa2UO79cDE4b2YcOiAzbfiig+ddI4LLzmTL8h3ldKqRL7Lvr0qEvUYdQXdeKQULEPadSwojF6dokkDrpYKcl8V7+cpAMCos803nCsL1ft+cePwrFjBrAOda46Cnn8qMW2I3PUHnGwvAvdnv+8KQexixspntS6TFzHHi38tlRIX3B86BPoahpbJHfi4Jpvn5Rj2ogDxXcsusjVFKq52oN4TBY4OmJ+Dinq0ONfRXQOhpKSwk4khZuI6zmi+QQEu6u2cY+uEaYbSgUM7F3vtyPBMcEALjsp5Zn8NHTOQc6vpACWkiPQIxZLJ1RdrGTTzYXGtmpakhYofObPnzYB//z+Q9jcvJ8FYerBAxPmmH2iBCHkO5A61B5xsFzTJ6PxaEW4U/wn3owqJsvCs0KZddQIAN6uQg3HYcLUsQOtimEdrsShyXfsSbJW0iNj6n0lT+3Sy9m4084uJ+EfTzw4ORNTb5JsX8Vgf6GUUHUO5m7hlhj1Km+NFrYrel3fcZcNi6vNcsz1nGnVSKJyayUuT7IOQ9c5SFHY/gTiEB6+FE3vVVeKtK8tyBdf1CW4+WbaxKj5VK/5JFNlHSo3yuV33UNVeshYJUgXVaobwLaoFzWdg0kemsa086FX1+Pmi8NgY0KE5b5z3YdRKhD++vp77FS69/IZOGRoHzTuasaYQT0x84fPpAi859a+vT7nkHbIFbRVpVggoM2t3iNG9ouJB0xQJ6Gt7Lg+WsTuN+G686LuNy6mrPLxjZZDOufglyQ5ML1VatiJ8Dsd5+BylCgQXWB08ZjJIsbFCY6IEsVtcZ2D972fOfuCq0+fvz3qpc5B5jO3L2xn/L0axUrMb9uG4665q4x1skeTJpg260iKEpwlao84WK7p0UdN60oaxbCOtrIIypUThZOBAsDQvg3BR7Yvqd5e9UXs3d/mzGH09idXsrVSFPqiWwpkvskV920opfBfqEyMJv+6WPDoHKI6kYXwwn+/tzOqDA/PDwB+/MTSWJnx8w+876894FlvrdesdCjIJ9DSVsZv573L6rsKFA8OKTFtXBgMzrbYqKIJvZ3f+dMb5hsVcO9l4bvbEnUYepymQoGwesseq5c8ECqk9cWxpy9WKmrjT2VEeA94OycsEdE5BEYO8bzSuuz+l9dG0p/+yqlB+Tad2YHIOdSeWClB7qdet8XT14s56dB4lPFTDx8a85Pw/By0xcjQHp0Fdwljcc05k3HkQf2Mi/SMQ6LRJI8dMwAAsHrrXi57UG88zf/hfxc0haANLkROwp2IaPch2i4b4vL/cPEQECz3IdOWrN+BHz72TmKZsjxpwmxywhIIo+lyJ8mZNibD+zXE0kxgT0Hzy1y7jR8HMY9+5r3s2d8W0WFw86dOi0R50ICe+OAPnsa5//P3hDZ7faFzNlIhHcZW8tJVDpsVK0FPSzawsImATIv28H49gvJZziHlWt+RnEPNEQcb70DaVdOh7OriIfGbf54RW0SG9W1Avx5RhaoqVoq2imPT4/+TxEqtbWXr7vKKUydE/svJdv4tLxjL5CaDPpkkIXMZu+mIQ2Wcg4QT56DrmtRNgkH2L8s1BdiLi2Gi/+OTPNz5Jul/uD4umQYrA9WCLhADyW9D1XFiF8/Tv2ddonhEJ9YzJzHxpRhITkBfhCV3JYuVHtL2A46SDQQkImIlSRyYh9R1ckG5FH5zY9lWpsQvnguPeLWFoskaNUccuDEzpI+nkKwvFbQXaFNIx9P1gauzkmu37eUVbwaFtD5guEVVt97Y2dTKKr/U0NPvOyTkcrjdabwdbJMjqC9KBW1Yr0mO3Ku+mJkz3+Uf8KxGYv0n3HUO+vNFTFmZ61653tQxOWfF6hVh2UDcxFmtw0YPTSFUdHGNbbFhFdKGOq/+8CQ2ndvMqGaqrvtb/f1ePH0MnvnqqbF8bUGYl2i6jLekB9Szil+YzZNZIa1wDv5dXN6kY1M5XyH1PhtRve7hN4Pf+rhpT9QecWDS5nz8aLx67Yc84qCkm16YbWeuQldCff43CwEASzdGlbGu8kZu1/jdPy2J/G8oFQAmnzQJFULg388+IkjnFJBfOmNirF4Vk0b0DRcfv566Upxz4MzufvSJYzHZj4UfiXKZEMhMCIGn3o4f3SpP6ooF3vO/TTvqvyzeEPyOEwfVCc4TjwzsFeUAkzgSnWOUzyHDLrRofaMqpJM4B27h08M5mFrXt0cJn/3gIUo+SdThf0fLvtw32XRxglPrdRWX6PV9aPIIjGUOhwo5gWj++oBzkGKlcLyYQMx1nZj+/JPTYvcFJ88xndtmWORl/6pGDiqEli8Juc6hHcGNmYZSAf171XmvR7levUI6qsMwW4Lwfg56/URxdlnf+V960jh2d6nuZAcqppsc5/C5UyfghatPVx8j3mb/W+5CSwU/DLlf75bdzexJdOcfPzomHwaAH31iSrwSJc99C9bgpsfjit/Qvj2aLh/fZPv/uV+/Ej4LY1aqcw5fnBklmMVi/BlU6J71MpuMShrjHCCfwy5WKhYMxMHRoeT5r52OSSOUg2oCzkEuqtH8UgwU99sw7ZTDHfyzS5NPa9SLMUU7lQx3nHOIEgdO/BQnBInNwqG+k5t6q3SW5BZyP8I7JgyNekPL10JM29O2CQD+8RcvdViMpRokDtzEUqyGlHQTNXcxZW0oFfyB4aZL4Iob4SuzovV6Oddu24vGXc145PVo4Lq6YsH3S4gu+nLSEhFGDeiJx7/8AZx//Ch2oJUKhIMUs0irzbxySp6cAOWywNTrHsfs215k75ETRg2tYTpRTz7vOkPoDPlctzy9jL2v6CCL159Ofb9CeNf1BTm0Voq/uVMmDMH1FxwdSZP5DvZ3xfoJc6p4xza26ooFVimpi7GMiw0z9iR2NbXgheVu4U4SOQcAv39lHZ9JLUcryEwceE/+QKykHRNqk75wYmH9/eqHILW0lXHtg0v8a0z7/ALlmS96OSrRjCBhHeE2Arrkob1Qe8SBSZMTS2c3jZwD7JzDjz5xLB7/8gdjYiV5j+4Nqg/WEf164BPTRsd1DkSB+/wpNzyFE+Y8bmyfvjPVn2XCsL7o01BiOQddaRjnYMLDiqQOpVSkQMcyNyGeknyumx4PrXxM3ZnEocl2PLhofSQ9FCuZCHy8Peq1KOcQ3wwULcThmnMmB+bHQXv8bGcf7Tk/6s59qngniXPgvIn157RZ2kXzhe37viX4m14aq1xF3GooCXo2k3gx5Bw0Il0MxTZqu2yGG8SYlcbCvsj2+clqYDxOnxOOF74cqRPSdY6Bn4OhrdzBX9c+9DqTM3vUHnFgxowcYPoO3qTUM/klSJw+aRjGDOoVs6eW41W3bNDraS2X2ZDAruKsAhFatMnBWdWYRBRcvTpkHH4pPikVCgExbG61s73cAmKyLklqnh6rR0IWZzJlVXfaepaIE5wQ7GZALsacfUGxQPGzkjWxjX495ByE9ZlLxdBbfdmmcAd5MCOn56D3hrpD3mMJ6hYTK5ny+a+DIx5XnBqPTqr3a9Lxm/HYXd53TKyU0lopFjtM0QEBwIJVYfRgjtuWxEh/d7J9Uje2W4u4HFor8W3lDFgWvrudz5wxao84MMM6yjkoeS0WCEIIz+lHwfsnDgGg2l5HJ9Ewfzf5sSlR1lNvV0ubiJ2Y5ZXnZgJKFOdOeF8FN09v/V5CKPL6gB+6u75YCIhmUpl6HBzAOxntp/80NXbIS9Lzql7BqpJ5hm+RZXqHtlAT0dhKABhdT9EQLsS73xyVVZYb61MK89meua5QCAj6GT96BoBnIHD9+VExlmmxielXlN92mXj0PmN8J00xrEIPU8LVqRP7n/yjF13A1CWyPlWsVC6LyHnPOrz5G02Lb9jC8oDo++QshrjjSb1yvPv6+QYhu3TiIPMZeIcWB2vC9kIicSCiO4hoExG9rqQNIqLHiGip/z1QufYNIlpGRG8T0VlK+lQiWuxfu5n8XiOiBiK6109/iYjGZfuIGjjOIVC62b0YJeSi/78vvhtJv+2SafjT50+JhBFW2dthfRswakBPfGLamMh9OlFqKwuWcygmHBP66Jfe79fLiJWYN83tiNXF9l8+eCh++k/Hs0HiTpowBM9+9TR866OTcdrhQ3Hdx49yjjlVDHbd0byzjhqB8drhQUnFjegf6mVUJXNoEWVvA2DQAQVKB36nOcJ3OjNxXqbF1GS6qMqkZd1cOIxikWJWYKMH9kTvhuqCHSRZSXH5Oeg7eIA/MAtAxIBBokFT5J9zzEEY0Ksu1n8S8j2qYqXv/vkN3Kx4retNJcTnucnsOhRPqma6XuohSmTWNkP7JALOwcCdmYi5yX+iI+DCOR97TogAACAASURBVNwJYJaWdjWAJ4QQEwE84f8HEU0GMBvAkf49txCRfNu3ArgcwET/I8u8DMA2IcQEADcCuKHSh3FBss7BntfL5+3g92oRV3vWF3H06P6RNLWMppYyejdwx45G87W0lVlZuSl4FwBMGTMgsELxxErRwa6HFwD48ys+ffK44PfVH56EWUeNNI7cgwf3QkOpiF9+ejoOG943MPFNWmP02PtQ7ulVH13kkpz+uPM2VHDHvAL2sw9Um3TvuE6KtePkCR6XyHIOTH2hgtsLvGjawQsRhmX+6T9NjZVTV6CYQnobE3rC9NwcoQf88Wfp6rgpazyzECKmGAaAjx3rccqDfM7h6a+cinnfnMmaZnNiQt1AINr+KHFoKwv85qXopi0GlnPQDTjkQ3lfJcWPhBuSJrGShBSX6Tq+pLnC6Rw6ConEQQjxLAD9uK5zAdzl/74LwHlK+j1CiGYhxEoAywBMJ6KRAPoJIeYKb3X6lXaPLOt+ADNJnzkZgnsZcmAl6RIkJFuadPBGQVn1hRBo3N3Mniynm7y2lkVkMIb1mnUEKi1ZsHobVjTu0a5z5Zll6aay5X0ciLxJktSFgTcr8yy6zDlJ55BEHIb0acD3Pn50LH2X8u74w35CHQEBuGDq6MhOXg5RjngN7uNxFWccMRxnHBE9FEiGbNfBiZW4ft7fVsb2vS34m3K8Kvc6TjWc1BdXSIeydWvsH13nkMA5qP3yf99/CG6++Dh8/DjvIKBxQ3pjWL8eLOfK6RzUfHp3B2KlQAwEs3ZXKU9vvs456NZKETEh8/AmzkZCEj19sZcch251qJfbGahU5zBcCLEBAPxv6QM/CsAaJd9aP22U/1tPj9wjhGgFsANAPFBRRuAmgDoZIx6RhvciZdJJ3sXq4nvnC6vw8uptaNwVD2HtjUO5GAlPrMTIgfSosSrUBccW4VKlu5zOgfModnXQCTkH+4CWbVV3wLo4SUL23xom9tPCa85kFxM9ftTRo/rH8qiIcw7Rk+CIvBhAf1d8P0hdjBS8fd2swOHwF5dOw1lHjgjKkc/D6chln+xqajHqJQDgDf8EspufNItNvPYRjhkdf+64fkWWwR/rKaE79Z3z4+cBAOcfNwq3XRJyOKHOIcxbII97iFuFxZ1JOc5BNTGPK479b8V6LGm0kkYdduxtCaITP3DF+7w8/jVZnbppYjkHhrNRHSelk6gu7m3zicWOfXzgwQ70eYsha4U0y1Fb0m33xAsnupyIFhDRgsbGZAcbDtwEkGy8vqMw7aSkk1nSIqiW97clGwEAmzjioLRLdSrj6nU5XP5TJ42ztstWHmfd48rHyR2eqVfOOGJ4UC8Q3V2qYpCxg3sFv2Uf/2Fh3GZ+YO961grrsOF9I/+TeBlOOVwWAht27MOjS95jiWMgxtBl5qUoJxPqErz/3jkN8fIkUb716eXBM3N6orU+kRw9IOwjo96DSzO8yySdwz+c4OnJTp4Q3bdddcbEQFxUIFKIZliWsVx2/MWzyfnWVhaxxVXuxHXrIhWHaKe2qToHIQSO/e7f8JuX3sWwvg2YOtbbWETEbYiOkWAzqZQZnCOhvIvPnx46TtYxnMOfXl2PX724Ot5gptzOQKXEYaMvKoL/LeMarAWgaltHA1jvp49m0iP3EFEJQH/ExVgAACHEbUKIaUKIaUOH8mxzEriujnIOYbpJ3FcqENr8wQoAJ44fxOaTMv0fPfZOYPs/uE/cYkOtVw5+3pTVLFZS28AdR1hS9CpBeYhbP3GLretpbvI5uAk6aURf/OLSaX6+KOdg29knzQ1uMXFhxQ8e1Mt4TXJUF/tOfJxZpOySbXvsB74XgkUmXIy4/hzii6Kmjh1o5Rxu+T/He/n7huPIvPgmvzg1i97XfRQld7FAOO7gAbE2DenToDh6IfCBiRIHvm4u3hA3/jxxJfCJn83FN/+wOHJtj2+VpAfek1h1/dkY1jfqTKrON1Xhq3KhqrhN3iOh2CooaXGxl/oknJ/DF367MCb+1eEaeLI9UClxeAjApf7vSwE8qKTP9i2QxsNTPM/zRU+7iGiGr0/4pHaPLOtCAE+KduwRruhhvuWJbsVgUoYWCoS2tpA4mBord9Kq5cRwzetZ1ishFcl6IDXAmzimnvnSGYeF5TET7IYLjsFnTh4fKFKBuKkt4PloxNoXU57yi04YBiTeSLXd+qldA7S4RWrpZSGwz2KWyC0m+u6S67PDR/RVrsdFFWUhglPsuPrlM9wzf03smgpd/CTPEOcwakBPDOhVHyysHBE58ZDBGNKnPvJMxqBvlvaEeeQiGBcr3XP5jMj/IrM5CSMBhGXrx9ma9WTxsDGcWFM6V3LntcvzSPTYSjaoBiBNik+O6oCnE3UVsg5VB5VkrWTSOSSB25hIDry94WLK+lsAcwEcTkRriegyANcDOJOIlgI40/8PIcQSAPcBeAPAowCuFELI3r8CwC/gKamXA3jET78dwGAiWgbgy/Atn9oLalcfMrQ35v/7GRjZ33vJOudgip0ecA6JZjnxZdLotAWBTbua8J2H3gjq0FEoWOzLlfzcgjmsXw9c+9HJmgmnR2xUn4hhDPFyhSQ2rOhOtZKSu7yy5NjMO9yyAO5/Ob4AXzJjLAB+MVm9JaqfSLvTkAppGdJDt0oD3EVtMc9dwXMEssxyWSicrJkIRxwIDesNd7tR5yDiDlcTh2txghinyZLv3wJEdUnrlXAnNuIVdxrjdF5m7uNzp02I1C11RDZI59Rlm3bhl8+vUtqpVRpJCwuVj/Pji48LfJvumbfGv2ZYMyokDnyE7o7hJhKNo4UQFxsuzTTknwNgDpO+AMBRTHoTgIuS2pEZlH4tEkXCHOgmpabY6cWCF98meM9GtjlOHUwyfSGAf//D63jsDU83UTSY9EmCdFD/Hli/w5uAQ/rEz0B2gZxQSef3mhy2dJC/qHLd0RYhDlGxkq29psl2nm/5wvXntr12UU+sDu2/XHx71hWxDS2sQ1VaezqVczDdW/Q3HfKRjUQE0X4xcbimsyb0sgCvD559J9TjfeToETH9SZHiPhZQ2qIu7KqOyGZE4WYdaM6nOpwC3oYjyYBCTsuP/PfzkbG/cnMo4lEPLZL3SMgRPqBXPaaNHYTnlm7GnS+swrc/dqSRiElJQEtruoXddu50e6MGjwmNL1Ih4ialHEr+Diq0bbbtjJJl+gRvsKo7eLOHtHyOEFdpEUNdTj/zyvO+k87vdd8lmwfuj5VztCmYyF69tuKTzmh2WQCTpJT6ZakA5WzT//PcI3H4iH5Wbidalrr8enWZFv2i/35tYiVZpovClyP6MbGSkQAxYs0CoZlZ3IJQ1nxzzfODUoiBEvLZYl1x4DglvU6At5pRl4VQIS6vhRfVrpUSg6SNmA5OJNdR5q21RxxU2TcbydJBVuoTB24XpZenl8B6Kvs7I3XR4BTSaihudYDokUdNiw/XPsDBJDfhv3qlLKJ9fNslU3Hi+MHor+gVdJ2Drb2mnZhcDDixkj53otwgf+ynCm/R4tt1yfvGAQA2aGdA28oCdJ2Dpd5yGFvJ1C8FivaLabHgQqabntxlvSkUKIh2Oqh3PaaNlYER4rtrFTYvchf1osv5KYHfiYtYyaU8/1s2728+R+/VEd4d6ggLeGfjLvxaccBTmyGtlbigiUA8Si9X15yPH4V756/pIKFSLcZWUolDTDmnWyuZOQcg5CxML4tTuBkPMvc9cSU4hbR60ItarE5H3K2LfLFSAnFwJTYyRLnKnfWqL0UIAxASyFDnYC7TzDlE5dwqbKar+xxi4QfvzdKuNL4fgO7nYOAcJEfKiDIidVM0Mmsaa0ejZ7beZ0zdRQpFWAUiDPFFskliMKu1kkObbSbcsXwJcZUA/9GY8pZ856wwj6YrUhGRLvjUsq5YwIdufDaSTzU+kX4OpnAYpujB6nPXFQpWEVvWqD3ioPzmYuDru0wO8r7mFhlj3iZWiqYZvWO1HQ9rteHHVlq1eU/EX0Iv0+VoTPW+JFY3Ln0ziCLIU6Cp3cFl1XUONhGNEHwhsgwnzkH5L5XL6l2xMMu+eMfWi84E2P9WPXxNzyvFRYGfgylfAdiuOE2lsYU3vcqkMCWAx6HKutrKZeVMC9l+/j41/pUK90CS9phiYfvI6EwWLY8nSmoMqKgwEBgzKNzZq/Ndhgb50JFRC6JffuoEzDpqRPDfFLJbwkQ01OcuFYkVVbcXao84RGSC+i4qutO36RyAMDS18VUxk5v1QPYHq3qJ281LsdJV9yzU0nUil07nIIlctZAmh2p/cC2R7ZPiLNtCa3b6k3XGr9mmzt7mhHeW0J4ATJ6/fPH98Wwa52Dyc/Dq9RbLe33zWLNYibB5d7g5SEUcGJ1DQ6kQEy2eEIiMQhQLYV2tingudAqLt/f75x/NBhCUbXFpOqdz6F1fjEXwJXIkDuDFaOq8Ua24vHvCa+q9R/oBHsdpIdNPmzQsKgko2K2VTO9QTZei5r3724ziqSxRezoH5TcXMygN5yAXcBMh56a2TUSjXlu1Oe4cI8VKsfgy2oO0u7WSIZ/klCK7UJZz0O+zcCKGvlUj6erQF5IpYwYEv9lnNXB3tnVLb/MvP31CEAk2ki/QOficQ9ksLioWCE+8tUkR05jqjoZ+/vKZh/EZuXuZyjnicCnjZS+tqQA/cnAQzFFpmAbp3Gdqi4vOgVvMX/zmTPTVQnq4iqlcoi+TNgbUblPHpCnCsI5CgVAqxE9olDClR8VKhK179uPdrXtxz/w1+CffnLu9UIOcQ/j7FMUhDHDXOcgBIVlBm84hfm88n+RY1Pzci/fESvEzhnXFuqtYKVBIJ8jhXa2VpDxUjR9lCz0hoYsdjh7tLeZ9GkpGzkE3s1Sh36L2Bxf6QK+BC++hg9NXcdDTWw1xswDEouTaxE9S1v2rz0zHBVNHs/lc0VBXjB3QxNVdoDA6bWubCAwhbKeZJYnmXLgeVdeW2L4KiU0sT8A5eBlVncBHjx0ZaYdudTWN4boATy/BiY8G9a436klUm5diISRqrnrAalBzxEEuBbddMjXiVQzEdzImayQ5UFLL6mHWOQggMpO4+PzSLjy2+FUoVpL5kmLGL2+MnllrIj6S/d+8JyQOnEhBv/2rZx0e+f9fFxyD33/uJIzo18O4szQdJ5mEJK92wK1fYiJJy0IOhIuHd8qfKW/0v0lJSaQcz+rKJlrQo67gJFpUOYdWRedgU0jbjvAmxMOPcygQBYHxwjSmfUROh+N44zQhj/8dGhIAHz32IKy6/mwcMjTqHKh7jsuDfXSUihRIG1Sv+/4967CvpY3djKhEZ7oSIieD156ImiMOsq/HDOrFKzOV36aBGzu03qKQ1sESB8Tt37kFpECeM1n8wBOtfSl1Dvvb7JyDDBoIAJNH9sNN/zDFUJ5nIqjqSw4eHI9hpC+ketjtnvVFHH/wQH+nyvejlXMwXuE9TuN+DpI4mBcaZ85Bk123tgnjgq6/t6KRiITWSq4+LTY0lIqJ5syyfdK/pyzCMRr4ZRjMtE1w5RyIuHPX+U2Wmwdy3JRVnjgXliXFSqEYzfA6fBNfETjUbtrVxOarKxaCDecR1z4apMtYaE3M8bqSGP/+cydhQK96KyHOGrVHHPxvrm+l1ZCE0fNUpw2GutiJy9ZLMYV0HReyO7Bm0erRHkYt5/rzj8Ydn5rGt08qhhN2jSrL+62PTsYYQ9C6ppY2zFu5NbVprG1hNcaqsXAONtHCT59dHqtTlz9Li6bG3fEIukHbDGEodOj6i5a2cmDzHsurjRduDMgypVjJVYRoQ12xkHjut2xfuRyGjZFE7pjRA9CnoYQv+FFI/+PsIyJtNYHIjXMg5XnDNL59Lk5mnIJbj1cUSAf8sdxWFubzyH1xm9SvbGcOXwLi/ikSMhz6nmbz2dRy4xAoyHPOIXtw1gcSpUJ0sLpyDtValkn7bHUi8WE2vJ206ajEsLzw/5SDB+D0SXygrpBz8CaAGpdfxRdOnxD8NolEAC8c+brt+xLDV+hro21hVT2GVdiIgw0Pv7YBe/e3WrkL6chke6+kP4OBxMlF9555Xpmmg5wA732or9ImvpPvLAviUCxEFy2TyK5I0WjEch7071mH179zVnBut8oZ2VpHBLQ5HINJiItwWX0ekVN4Cp4TjT5zj7oihvSpx5qtnrNjWQgjR+75pwADfHHSToPFFOf3JOsCTN7QYR16ne2N2iMOFm/OksL2AclOcHqZOljuhMm3vHEPHl68IZElls5ASQpp9b9NJq2blKpxplR85OhQAefCziY6IbnK6wveDu/aB5fErqmTQ1cAJhHrydf+Fc8tDeMI6fn19zCCjaSr/Td0y9Y93kIho7d6R8AaFt9C1J7frHMIZdyuIkQbdJn5pJFxqyvAF5+U7WeOANG4YEnWeYnBK+GNA13/x3P+aTiH+L06+vWow579rXh93Q5s2NFkXJDlvOzX0xMPyWCA8XzEijUb/ACPnI5TdxRN8oHJErVHHALOIY66AkWUkGbOwbtbHvn5Xxccy+ZL+wIXr4uHO9Dr5QLb6QuEOobtbL33vdU/k8Ak7lB36S47lqYE66d+zHkTHEzWJ8989dTI/8NGRA/3cbFYaVJEaXpufVPwsSkHsW1TYeoV/W21tgnW+50r0yTGKFA4NtPsIPswRg6ynkhMIMP9RT8q8COLN1jrVomGrXnOOgdQTFRpMuzYtJOX9+vlqe9FHlYUa5/fL/LUO6NYyfdsJxAOH94X//LBQ/nyDGLSBn/ecYTjdX9N0OdmR4iVatbPwcQ5qEjiHNqEwGmHD2Xt2wE7S80haVGTnqLJOgeVczDTf5nv+kfe8vIaFq16x52gRFOCDoM704KDyTN2rOZwlDWHrcu3v6ZZUwG8MxkHXZ/jWfmYdQkuiIgfU6wSpn7SOQdTkaVCAa1tZXz1/te8/4bxooeFN8Fd54DgbI0gjcnHHcFrKk+dQ0MNvhh6vxj7rxCG8Le9DjJwSpJz4K795KllAMI1R+bIOYd2QCjzY3QO/mCXB4h80HBIu9xB7G8txy2XFGT9/qRuIuloRXWHYzMl1Ae7iXOoq5BzuFc7LEbCZmmkt8/Nbl0T86XUAelio0kjQ05k8sh+bBDEuFiCL1u3QGlpM+sc3MOehL9tOqBrz5nsVL5+ToOpxALpYi9+vER0DtbFMhoWf6QhzMZra3fEPJ9drLRM/gaB6XhCG6UYTcJ2aFG57IU9sb1DeQb8Ro27CQNRmjdVsq9za6UOANe30jpkz/42fPioERGrCxXq4LfJ9LkXqJttqkgI8hqIWfTdqL44RhYPK/GK3meyjlHFILbBL0+Ra24t48Txg3DiIYPZfGk8uF0W+mo5h8M1sdQNFxwTlm3oPr1Zpib007x4W8sWa6WU/ilJ94zVzIhN704/uc0ENXwGYDvvxK19qtUVkO05BdPHD8KPNfNUCT1MjjkQYnRzst+g7A4DJiY9r1feid97IpLeHFhEGW8NzJpDJzhz3qxQc8TBpnNQd2FHjerP7hgB7dQ1y86NwzVnTzZeW7fdHgZaKix1O+od+6LWQVGrJ3N5+kCuK/HPohIN2+A/c7JnFdXU0ma1JkoT5dU1MJuKak+ZVc/gNil8e9cXMWFY6AxlEp9cPP1gAGGAtuaWstEayHXCq1XZiLWr6Kuo6drMBgIUiWq7ggnxAkQ3JLZnIkSJTZYB5b4+a1JwwiNXsVqTaY7oYiWTwYhUrNvCsct8HPPRrJjLmiDPdwk4h9xaKXuE1kqcWCnsDhtHoA5+7lAeCX0R/OBhQ2Phq9NA2v3rY0iPX+Oqc9BbbgvrIGFbjOSlppY2qwezK0dMjjtavbxql5jI4maxGPqvC49R/vNlFQuEsYN7haFKWs3EIW00XcBuraSPcVtY6EVrtof3GcrbpMn9TUHuXHUOurVShoyDffxplZmj30bHn4k4FH3/D1vEXcCfvwwBeL8fxsdGHGzm6u2F2lNIJ1grSVgXQWXcvWexjtCtUioN+RDUSxQZXFPHDsSNn5gS80KO2MrbJqfWHJMVjQqXxaippWzlHNJEjXWJ2qlPlGo3oGnEIhK2J1LFY82tbUbRokqIDhnSm83jlae01fLObIYKKrZpTlumR9bNRM8/jo/plE7noIqVsiMPNvGtftiPkaPyxV5yHNqIQ5uQvkrmNhUL8c3Oc187DUs37QLAK6Ql5MY1VEib68kKVa1WRLSKiBYT0SIiWuCnDSKix4hoqf89UMn/DSJaRkRvE9FZSvpUv5xlRHQzua4eFcDmYahyDjbioO4sX1yx1ZxPKW/6uEH43vlHp2hpHN4xjYoJphBseApXsZe+WJjEaJF7HMRUTa1tEQsnHa5vl3MaOu3wuJGAXpxu2poWER2LpbHqNesiiFBkYuMc1Pdx02w+RAmgRQW1EWu9vSbOQQ9qZyB16qI/68gROGXiEDafOuYSPaTbSeeQxDlEdQ58PhkWI3onk49CQxE7Zx0nDmMG9VIU0hbioIuVuohC+jQhxBQhhIzRcDWAJ4QQEwE84f8HEU0GMBvAkQBmAbiFiCR5vxXA5QAm+p9ZGbSLRWirFO9cdT2ziZVc2X+1jG+efYQ1fLELiKK+F6ahVKnYoadltxWU5yBWEsLuwZzGZFOdL186YyLu+NQJsXxqULY7PjUN/2NQRLoiwjnYAsdFHsMuThAAnl+6Gdv3thjjGKl03NZHL6/eFvy2WSvZwqqocHFE89rntulI4+fwxoadwX+TCNHkuW+DNA/lQBQNqmjq65IiLpLt5SAX/bIQCWKl6Hg+0Q+kJzebbmIlsxNv1mgPncO5AO7yf98F4Dwl/R4hRLMQYiWAZQCmE9FIAP2EEHOFR85/pdyTOeSOgevcpPAVEq6RMFWrlCw8WdUY/oBZ5OKqsNQvuRA9q5hKudanwaxbce0JqWOROHPycHbyzVCsok6eMCQW5z8tVM7wnY27jfmKjougjPZ785NLAQDLNvFlqmPONlyOVc6nMFmYcW0yjenYomzZSUtYOaoKxHIAMNDgjFZJqAirWAlh/CyvHXy+osY5JB/vmkwMVS5NGj7ox+ZykGuJbE5XCJ8hAPyNiF4mosv9tOFCiA0A4H8P89NHAVij3LvWTxvl/9bTYyCiy4loAREtaGxs5LI4NdgEdeHJhHNQ2esUPf0p5qAVAPiNcni5l28smy8iVrJZs6R207MTTXXujB/CB+cD0plsusT8cdWxuEKtR3qP8/Wqi7mdCAsRLgwmU1a17bY++tDkMFaWjXPQX69pTDvSBmcTbpsRRKQepYje9UXcfdmJbL5KRCh9Dd7gXr0UCfFiO3FPNRs3x8QKj3e1tVXXOQRe7v49NuOLGKHvAmKlk4UQxwP4MIAriegDlrzc0whLejxRiNuEENOEENOGDuUd1BJh0TmoLyDpJbug5GgCesK40Fnn2x+djG9/7MjEsm+44Gh83KAQjHrQmsuoZPPhyjk0WBWCbnUVPOEwAO9ciAnDeF1CRIxmeCj9zAjAW0Ae/uIp8fIcO0Zd4213yMVDKjTrHHQOtvGi7opNhIYrw1SmvmM9YdwgNl/Ev8Jx85Skc5C4aNoY43GilZhtWsU70BXhfL5igSJir2s/ypuhS86hrWwnDvpmR7ZBEh3OW7xvQwmfOXl88Dzht7GazFAVcRBCrPe/NwH4A4DpADb6oiL435v87GsBjFFuHw1gvZ8+mklvF9hMWZPOUwiuOe7MXZ3HfvKPxwe/j1FEBjq+OHNi8NsqtiH1t9vk/Ncz3I6atE1Utf9slk9pFNJywlxoOe3MZQd/LhMf6aozJuLIg/q7NYYBd+awCWUB7Pfl3PWmHaijmKpnhDi4i20+c/J4Np9qhXPVzIn4V8Oxo2pxds7BjWNWS0gj/qwWhUI0bIfxnHKt/4b15T24C761UptImh/RuqQyPjgQiiEO+9vKEf+jw4Z7vjUu529Ui4qJAxH1JqK+8jeADwF4HcBDAC71s10K4EH/90MAZhNRAxGNh6d4nueLnnYR0QzfSumTyj2Zw2bKWgnnYJucJceYRGp5Niuf6cqObnAfXj4LpD/sB4h705pg9whX89msRVzbF565a18Ek8viFp++jgEAjWVGTFnti7QQwJv+LtQoVlI5ESvnUHDKpz7ywmvOxCdOGMPmUx3bJo3oa40+GvzOgHNQi7BFJM5CVKhCN5YwB9RzK69I3sLe2lZO9HtSiUMY+pwi/yVuf34lmlvLmi7PPxgoIfJxFqhmdgwH8Ad/cJYA/EYI8SgRzQdwHxFdBuBdABcBgBBiCRHdB+ANAK0ArhRCyCe8AsCdAHoCeMT/tAsCayX2HaoyVfPIiBIHcz5Xv4kiuZWnljHYoLwD3NnwQoqdr4Tduc1NJq1eGp9gzy+jcdrMbF3kr9wC4xrjyVimKoNP4JRU00mX8Bm2V+hicizr5X7rUM9USDLFDPJZClT7IkncJmE7ICprb2B9A2ZTSLtAipVaLXGzAK/PVEdCGaI7IA4aB/Off34DQPRI0Z5+3Ld9CZGPs0DFxEEIsQJALFa1EGILgJmGe+YAmMOkLwBwVKVtSQPbYT8RxaZl/qmDxrZYRjkHS3lFt0UmQhwsZrGuU0md364KLvsi7ZZPJSKcuEdtU0AcbIuWw1rJLTDVOiW6ilkKFHW6OmY0L8py0Z0k1RVtX7K4DfDOK9/V7FnvuOowbKbKUSc4c71quJhzjrGPgywR4xwM5du4eP3+NiHQUi4njHtgmXIeexLnIKFyVVKkmBT5OAvUcPiM+LXo5LQNfiV8hmUwlBydgdTJZBuQ6oIxwHCIOeAWI19vUxamca46G7UqebQkm68QeuXa+tnFZ5LlHCy28C5wtgqjqLLxslN42b+rOKaSMBu2cfUzxY/Aav2kwEZY1bljq3fDDi+6wEH9exgd6rzynJrkDL1NprH1zDtuFpGL1mzHvJVbsaJxD55523yPHkhS8tdbnAAAEKdJREFUt1ZS561qaqtyVecf7+nepo/nI85miTx8hgKVHtimiOskVm3QXdl1V87Bxmo70gaNc3C7xwa1/2z29+RIlFRTQtf+M7ctnqe+WJ1YyTmGFRH2K6G7jYHtlGTXnbkNai5bearfhGs/2/wI1EXO5ThXm2WbXm8WUNs0sn8PXH8BH7lADytigqoc3t3casynDxGdc4jGcQp/q1GDp48fhFXXn+3UrmpRg5yDD2a8uS5a6rVrzuHDegOVcQ5WTsRxUXCNUVPpoTEmuPafKwpECufgprMxlxVP61lffayroA02nQOik90EeaQokLQz9+qalBAmxJWzUWEXK4W/be1Tr7nE67LFkfLqjZYxzHCcrStU4jDn40eZo7dmDP0kPp04qKFEWv1xf96Ug4ycZnuj5oiDZB14nUOY1qvevJtRJ9qso0Ya89U5OsG5KrhdJ7hrOIQ0BOH68482OudJRBXrWRCHUKxkP9EuuSyu744bY2bNZxsse0z1JingbQpXicff3Bj8tu24bWNEhauISIVt3Bccx+mYQaHlmwvn8K2P2v169HdX7T5GfcaOcCaT0E9AlGKlILYS4yB3wvhBzoEqs0bNEQebtZI6Bm2htV138OqC5hrjyLpDdqy3ErHSu1v3WvPOnn5wonNee+gw5HyxWwO5cA7RPGdOHm4VzV2vHPhjQslRbEhEaG5Nti5RT25z0T0lPbfrOFVhsuUHvNAkEq5nL7gQMi54pApXZz5XHKRwCrb3dqzF50iFycBAh246La2V5NhW/RwCE+40oRUyRu0RB6ufQ5g6aQR/LjRQYfgMxwGdBedQyaEpNlmpK6Jih+rk+YD7TrUSzuGbHzGLAyVOGDfQ6hyobiCSOAcXsdIHlYizNqsXWVfSY7uGsVDRr6dZDakem+s6Fm1ErreFS4nka4jmM82lM5WwIjYMUN6b7Tl+8clpxmsqvj5rklM+XcclTYjlxlE1WpAipkq4v6xQgwpps4e0KyoJvOdqq22N6eTYZs7TkoOqODz7aLN4zBVqn/ao0hIISCG2cehbve9c+vJ3/3JSYh6JJNNdKVaSkTg5uO4SZT8nLdCVLCyubXD1tbC9m/n/cYaTZZ2uEzC9utsumYrx3/hLYnmR0PyWcTDI4kukQiVW9xjOTQfi1nGBWElTSL/93i6s3OyZvHZEgD0Tao84+N9cl7tuuEvFAhb8xxnobzEnBdzDbKiwES3XMsb5Cj5byAkglHF+4LChmDi8ujMQgKhIzGbN4grVCsTVNNaYJ2O5tQ4b8SKEz/KPJx5szOd65Gyz7wBl0w94bUpPoLPYxLiiV73b8qPrLeTRqzpcN3zqOLVN+TROcBI2ZblU1HthNDxRLRA3ZT3rpmeVtnaecKf2iIMl8F4acYzL2QzqpiiLDYAtZIaKQ4f2wSvXnImBCUeSlgM766qbBiC6I8uCOCx8Nzy6MslUNC2yJg5JJqDSlDULa7Q9vsdsb0vkUaB9RRIfPmqE9fq3PjrZGtG2Uiz/3keqnkuqqMs1hMrNF5vPCFFf6aiBZssnSeTqigUs+c5ZwZiRmwKOi8qCCFeK2iMO/jdnreRq5eOKIcpinoVVRK/6Ek47fCjGD+mTmNeFJW7TrCWqhTqQq/U+BqKTxbbQuYrRVGRtpWI3OAh1DrbJ7rpLPO7gAagvFvC5Uw+15mtPZeaAXvbx9WlDkL9qkYWYRe1n18CLNmIYiUZs0bXJa2UhWNEWRxxyzqEDISwa6SzPsAWiO9osFksA+OWnp2dSDhBySlnJNesy5hzUyWKbJFyo4yRkzTkkRat18fR2fQ1D+jTgnTkfTsznKqaqNZhCplvvsXJ8buVJzkEnAqbYSuq1zkDNWStJsGKldgxXkiT6+K8Lj8EDV7grQLOADNmSHXHwyhncuz6TMtXJYttxu4YLufK0cKfdkfbtRPHY/Xy+bNvUmSKJAxlpfXDebwntAbjF9gLCDaI+XOV7amUs2rLYZFWKGuQcvG9ueFRiApqEf/3/7d17jFTlGcfx77O7XIIgt7KKIJdUJSCKwMZLwGqsN4wRG9tGY4SKjbWVaJs0UXuxaa2JNq2plzaWKERrW5q2GrG2NbTxElurQL0iIkgvSqnUiLpgBVee/nHekWFnzszZ2TNzznF/n2SyZ99958x55p2d97zvec/7nnoEL/3nnbr5PttV/6artJW+fNOa9bJ0dpXWl1ySmUwhecth8ph9d+K28mtz//tYWnc+Vnqt8sWkBoK4RYNKkk6oB/DydxfUPdFJeiIU13vQ0d5Ge5tVvVEyjVF/jRp4lUONxX4a6J2o68pT4yeWy9reXhN/9deHF9hS+jyXtwhq/QMmvebgZWNT0qrAnrjmlLozZJa/Ur2z1jmTRvGpObVHmSXV3mY8eMV8poytPT3FR0n5hd44wwZ3MGxwO188qfY1G0h2h3fS/59a+xrS0cbung8quraTvH6zDLzKoUbLoVQw4/o5d0tRlL4s0+5+aGRt6mp69rvmEL/PpC2HI1IYrttbknl5yt/eet1Z935pXn8PaT/9WemuiOqN4ILoC3f9t89I7QQhacu71sXqwR1t7O7Zy+btO/dLT+t/qREDr3IIP6sPZY1+DpRm+IKZ43lyy5t8tcr6yo0o9acfOiadicz2m92zxgyqR02IvgDHj4yf+gFg9qTRHy7MsqfGymNpe7hsGufu9/p/J3qW5kwaxZY3dmV9GP2W5vWd5C3v+JOYIR1t7OnZy8bXu/dLz/Ky0YCrHBbMPJgjDhpetRbfm8Ld00UydFB7ojmEkuocMZRbLpjNvI+PTWV/5RfoOg+Mb80dPHJo4mmMz5l1CPc9vTW10WN9lXQ51lb78YVz2Lrjf3Xzpd2ySeqi4yczvJ/LujZLGoMvdux6n5VrXmXlmlf3S6+1UmKz5fPdbqLJYw9gckwf7PTxUbfD6QnnaJFK58yKX9Grr8rP7tMatXHDeUfx+ROnJrqJsRniPntZOyuF6VOa6bpzW7JQZEOS3kg3Z9Jorj17BmfPqnyvq7Vk71jUlXiakmbITeVgZmcCNwPtwB3ufkOrj+GwzhG8dN2ZmQ4fGyjuXNxV9yz6OwuP5Nr716f6ukM62jPrh681dYYUV2kanXozEpgZSxKszTCo3dh0/VmpHFt/5OI+BzNrB34ELABmABeY2Yzaz2oOVQyt8cnpB3FYZ+0LxItOmAIkm6qkCK46I9nsnVIsZsY9lxzHqqXzU9nfydM6U9lPf+Wl5XAssNndtwCY2UpgIfBipkclmVu1dF7LVupqlmGD23l3zwf9XrNa8qvWGthJdE0ezdp/7gBg+vj45QJaKS+VwwSg/ErMa8BxGR2L5MjRE5MtuJJnq5bOY/WL29UqlVi/DrMjbHq9+8NZlbOWl8qh6m0HFZnMLgUuBZg0Sf23UgyHdY6o24UmAqQydX5a8tLOfQ0onz9iIvDv3pncfZm7d7l717hx43r/WUREUpKXymENcLiZTTWzwcD5wKqMj0lEZMDKRbeSu/eY2VLgIaKhrMvdPd0xjCIiklguKgcAd/8dUH8BWBERabq8dCuJiEiOqHIQEZEKqhxERKSCKgcREalgvVceKgoz6wY2VvnTJOBfCXYxEng7x/kUR77yKY585VMcjeeb5u7177Zz90I+gLUx6f9N+PxlOc+nOPKVT3HkK5/iaDBf3Hdn78dHsVvprYT5Hsh5PsWRr3yKI1/5FEf/8tVV5G6lte7elTS9aBRHviiOfFEczX/NIrcclvUxvWgUR74ojnxRHE1+zcK2HEREpHmK3HIQEZEmyX3lYGbLzWy7mb1QljbLzJ4ws+fN7AEzOzCkDzazFSH9WTM7uew5c0P6ZjO7xcyqrSFRhDgeMbONZvZMeLR0TUEzO9TMHjazDWa23syuDOljzGy1mW0KP0eXPeea8L5vNLMzytIzK5OU48isTPoah5mNDfl3mtltvfZVmPKoE0eRyuM0M1sX3vd1ZnZK2b4y/c6qO5wp6wfwCWAO8EJZ2hrgpLC9BLgubF8OrAjbncA6oC38/hRwAtHCQr8HFhQ0jkeArgzLYzwwJ2yPAF4mWvf7e8DVIf1q4MawPQN4FhgCTAVeAdqzLpOU48isTBqI4wBgPnAZcFuvfRWpPGrFUaTymA0cErZnAlvzUB7uBRjK6u6PAW/2Sp4GPBa2VwPnhe0ZwJ/C87YTDRPrMrPxwIHu/oRH7/rdwLnNPvZyacTRgsOsy923ufvfwnY3sIFomdeFwF0h213se38XAivdfbe7/x3YDBybdZmkFUerjjdOX+Nw913u/jjwXvl+ilYecXFkrYE4nnb30sJm64GhZjYk6/KAAnQrxXgBOCdsf4Z9q8g9Cyw0sw4zmwrMDX+bQLTaXMlrIS1rfY2jZEVoLn+z5U3NMmY2hejM50ngIHffBtE/CFGLB6qvDz6BHJVJP+MoybxMEsYRp2jlUU8Ry+M84Gl3300OyqOolcMS4HIzW0fUdNsT0pcTvYlrgR8CfwF6SLhGdQb6GgfAhe5+FHBieFzU0iMOzGw48Bvgy+7+Tq2sVdK8RnpLpRAH5KBM+hBH7C6qpOW5PGopXHmY2ZHAjcAXSklVsrW0PApZObj7S+5+urvPBX5B1P+Lu/e4+1fc/Rh3XwiMAjYRfdFOLNtF1TWqW62BOHD3reFnN/BzMujaMLNBRB/8n7n7vSH59dAULnVRbA/pceuDZ14mKcWReZn0MY44RSuPWEUrDzObCNwHLHL3V0Jy5uVRyMqhNPrAzNqAbwC3h9+HmdkBYfs0oMfdXwzNuG4zOz40MRcB92dz9Pv0NY7QzfSxkD4IOJuoa6qVx2zAncAGd7+p7E+rgMVhezH73t9VwPmhH3UqcDjwVNZlklYcWZdJA3FUVcDyiNtPocrDzEYBDwLXuPufS5mzLo/SQeT6QXRGvQ14n6g2vQS4kmgUwMvADey7mW8K0UytG4A/ApPL9tNF9CF5Bbit9JwixUE0QmMd8BzRxaubCSNmWhjHfKLm7XPAM+FxFjCW6CL6pvBzTNlzvh7e942UjbjIskzSiiPrMmkwjn8QDY7YGT6LMwpaHhVxFK08iE4Kd5XlfQbozLo83F13SIuISKVCdiuJiEhzqXIQEZEKqhxERKSCKgcREamgykFERCqochBpAjO7zMwW9SH/FCubsVckax1ZH4DIR42Zdbj77Vkfh0h/qHIQqSJMmvYHoknTZhPdqLgImA7cBAwH3gA+5+7bzOwRojmw5gGrzGwEsNPdv29mxxDd/T6M6IamJe6+w8zmEs2j9S7weOuiE6lP3Uoi8aYBy9z9aOAdonU2bgU+7dF8WMuB68vyj3L3k9z9B732czdwVdjP88C3QvoK4Ap3P6GZQYg0Qi0HkXiv+r75bu4Bvka0IMvqMAt0O9GUKCW/7L0DMxtJVGk8GpLuAn5VJf2nwIL0QxBpjCoHkXi955bpBtbXONPf1Yd9W5X9i+SGupVE4k0ys1JFcAHwV2BcKc3MBoV5+GO5+9vADjM7MSRdBDzq7m8Bb5vZ/JB+YfqHL9I4tRxE4m0AFpvZT4hm07wVeAi4JXQLdRAtxrS+zn4WA7eb2TBgC3BxSL8YWG5m74b9iuSGZmUVqSKMVvqtu8/M+FBEMqFuJRERqaCWg4iIVFDLQUREKqhyEBGRCqocRESkgioHERGpoMpBREQqqHIQEZEK/wdLy/Gu+VIpkAAAAABJRU5ErkJggg==\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], "source": [ "sorted_raw_data['inc'].plot()" ] @@ -2266,9 +2289,32 @@ }, { "cell_type": "code", - "execution_count": null, + "execution_count": 52, "metadata": {}, - "outputs": [], + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 52, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYcAAAEKCAYAAAD5MJl4AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzsvXmYI1d57/89tWnvvXtmerrHs9vM2OPxglcwNsbYkBBzEwj2cxOcewlOuJBwQ5afneS5Idz4EpIQLhA2h3BjnLCFJTgBY4xtwAYvzHgbj2dfe5vpTVJrrVKVzu+PqlMqSSW1ultLSTqf5+mne6oldammdN7zfVdCKQWHw+FwOE6EVp8Ah8PhcLwHNw4cDofDKYMbBw6Hw+GUwY0Dh8PhcMrgxoHD4XA4ZXDjwOFwOJwyuHHgcDgcThncOHA4HA6nDG4cOBwOh1OG1OoTWC1DQ0N08+bNrT4NDofDaSv2798/TykdXu5xbWscNm/ejH379rX6NDgcDqetIIScqeVx3K3E4XA4nDK4ceBwOBxOGdw4cDgcDqcMbhw4HA6HUwY3DhwOh8MpgxsHDofD4ZTBjQOHw+FwyuDGgcNxYf+ZRRyaWWr1aXA4LYMbBw7Hhf/13YP4+A+Ptvo0OJyWwY0Dh+NCRjOQUvVWnwaH0zK4ceBwXFD1PDI5o9WnweG0DG4cOBwXVD2PLDcOnC6GGwcOxwVVN7hx4HQ13DhwOC5wtxKn2+HGgcMpgVIKTc8jo3HjwOleuHHgcEpQ9TwAIGt953C6EW4cOJwSmHHQ9DyMPG3x2XA4rYEbBw6nBFUvuJN4UJrTrSxrHAgh44SQJwghhwghBwkhH7SODxBCHiWEHLO+9zuecy8h5Dgh5Agh5FbH8SsIIQes332KEEKs4z5CyNet488SQjbX/61yOLWh5gruJB6U5nQrtSgHHcAfUkpfA+AaAO8nhOwCcA+AxyilOwA8Zv0b1u/uALAbwG0APksIEa3X+hyAuwHssL5us46/B0CUUrodwCcAfKwO743DWRWqI9bAlQOnW1nWOFBKZyilz1s/JwAcArARwO0AHrAe9gCAt1s/3w7ga5RSlVJ6CsBxAFcRQjYA6KGUPk0ppQC+XPIc9lrfBHAzUxUcTrPhbiUOZ4UxB8vdcxmAZwGso5TOAKYBATBiPWwjgAnH0yatYxutn0uPFz2HUqoDiAMYXMm5cTj1wqkcMhrPWOJ0JzUbB0JIGMC3APxPSmm1XsZuO35a5Xi155Sew92EkH2EkH1zc3PLnTKHsyo0p1tJ58qB053UZBwIITJMw/CvlNJvW4fPW64iWN9nreOTAMYdTx8DMG0dH3M5XvQcQogEoBfAYul5UErvp5ReSSm9cnh4uJZT53BWTLFy4MaB053Ukq1EAPwTgEOU0r93/OohAHdZP98F4LuO43dYGUhbYAaen7NcTwlCyDXWa7675Dnstd4B4HErLsHhNB3VEWfg2UqcbkWq4THXA/hNAAcIIS9ax/4UwF8D+AYh5D0AzgJ4JwBQSg8SQr4B4FWYmU7vp5SyT9j7APwzgACAh60vwDQ+DxJCjsNUDHes8X1xOKuGZytxODUYB0rpU3CPCQDAzRWecx+A+1yO7wNwscvxLCzjwuG0Gm4cOBxeIc3hlOFMZa015qAbfP4Dp7PgxoHDKaG4Qrq2VNbP/+QEfvnTTzXqlDicpsONA4dTQlG2Uo1q4ORcClPRTKNOicNpOtw4cDglMLeSIgpFmUvViGVy0AxeMMfpHLhx4HBKUPU8ZJEg5BNrVg7xTA5Gnha1+D45l0RK1Rt1mhxOQ+HGoUP4jS8+iwefOdPq0+gI1FwePkmEXxZrDkjH0hqAQnX1UjaHt37qSXzpqVMNO08Op5Fw49AhPHd6ES+cjbb6NDoCzTDgkwQEZLHmaXDxTM58rvX4nx6dQzaXx3Q827Dz5HAaCTcOHUDOyEPT81iyFijO2jCVg1CzcqCU2sZBNczH/+jV8wCAeEZr3IlyOA2EG4cOIK2aC1KcG4e6oOp5+GQRAUWsqXYhrRnIGWasQc3loRt5PHHEbAwZTfH/E057wo1DB5DSzKCnm3F44vAsnufuphWh6oalHISaAtLO664Zeew7E0U8k4NfFhBNc+XAaU+4cegAWEZMLF1sHBLZHD7wlefxyR8da8VptS2qni/EHGowDs7rrul5PHNyAYQAN1+0jqs5TtvCjUMHkFTdlcN3XphCSjP47nWFFGUrrVQ56HkkszoCsoiN/QF+7TltCzcOHUDaCpqqeqG/D6UUDz5tprbyBWplqLoBxQpIZ2sISDuDzpqRt5VHX1BGNsd7LnHaE24cOoCko9CKZSztPxPFsdkkhsIKYjwouiKK3Eo1pLKWKodszoBfFtEXUACUu/s4nHaAG4cOIK0VjANbqE7NpwAAN144goSqF42+5FTHzFYSEFBqS2UtjTkw49IflAFw5cZpT7hx6ACSamEBY8YhkTUNxgUDQQBAjOfb14yZrVSIOSw3lNCpHFTdKCiHoKkcuHHgtCPcOHQAaYdbie1imatpnBkH7tqoGU3P26msQHGXVifn4lkcmllCrMg4FMccAH7tOe1JLWNCOR7H2dytoBzMPPuhsA8AEE3x3WutOGMOAGwlUMrf/OAwfnpsDldvGYQiCtCsSvVszoBPFtEf5DEHTvvClUMH4OZWSqo6In7Z3r1G+QJVM2rOqpC2DEKldNYzi2nMJzW8cDaK4YhphEuzlQDuVuK0J9w4dABpTcdQ2NylMuOwlNUR8UkYCHG/d608eWwOR88n7ArpgGIZhwpB6emYOdxnOp4tGAdHtpLfMjAxfu05bQh3K3UATCWouXxBOWR1RPyS7drgxgF47tQivvfyNP7y9otdf3/vtw/gwnUR5CngkwT4JOZWKo855Iw8zi8VOq6OOIwDi1kAQF9Q5m4lTlvClUMHkNYMhHwieoOyXeeQyOYQ9ksIKCJ8ksAXKABPHJnFA0+fgV5hYttSJofD5xIAAJ8kIuQzjUPSZWDPuXgWeQoIxPz3SI9pHFSHcgCAvqDCXXqctoQbhw4gqeoIKRJ6A7KdOZNUdUR8ps+7P6jwgDSAnJV1lHJxE1FKkVR1TFmuIp8sYENvAAAwFUuXPZ65lK7bNgQAGAj5IJDiOgcA6A/K3K3EaUu4cegA0pqOkM80Ds46h7Df9Br2BeWO2r0upjR86OsvIr7C95SzFIOzaJCRzeXhmPAJnyRgrN80DmcXMmWPn46bx37l0lEAQF9AhiKZGUvFykEuSnXlcNoFbhw6gJRqlBmHZFZH2Gcah4GQYu9e83mKH7wyU9G10g68cDaKb78whUdePbei52nWzIWUWq4cEmrxAs6K4Nb1+DARdVMOZrzhLZesx7uvvQA3XTRiprOWKIe+oMKVA6ct4cahAzDdSqJtHPJ5iqSmo8dSDv1BBYvWAvXEkVn87r88j5+fWGjlKa8JVpT2s+PzK3peNeVQajAUa3Ef7w9iYrFgHP79hSl8/icnMBXLYDCkIOKX8ZHbL8aWoRAUSURa06Hnqa0c+q2A9HJV1hyO1+DGoQNIq8VupZSmg1IUuZVYQPqnR80JZZXmDHz3xSl86OsvNufEV4mqmwv5z08srGjRZcbBTTmkSoLObOe/aaDYOHzr+Ul8/IdH8Or0Ekb7AmXPYW1L2PN7AzL0PHWNc3A4XoYbhzYnby08IZ+EnoAMTc9jPmmqhIi/EJCOpTXk8xRPHjN326WLIeO5U4v4oTX/2KuwJoJzCRXHZpM1P69gHMrfO1vUGSyNdWwgiJmlrP03E1kdOYPixYkYRvv8Jc8RsJQ1jS5TDuy7ytt2c9oMbhzaHFa9G1IK7RrOLJgdWVnMoS8oI0+BV2eWcNLq1uqWngmYCyjbmXsVZ6+jlbiWNN2KObi6lcxjGy014JMLyoHSQnZSIltQXKXKQXEoB9aXiSkIrY1jPJzuhBuHNoctaiGfhA3WTvboeTNXP+IvBKQB02VUeJ67AdANipxBYeS96yNXraK0obAPjx2arfl5hZhD+XtnxvLijT0ACov6OMtYslxLSVW3axs2uhgHVmfClAeLXaguhXQcjpfhxqHNYb7skE/EmLVYsUIuZhzW9ZhG4x+fPIXRXj98kuAalAUKO1wvz39gyubd116Ap47P46ljtamHam4lZhz2jPUBgN1XiXW1ZRlLiayON71mHcI+yX4sQxHdlIPlVvLw9eRw3ODGoc2xlYMiYaO1yz08w4yDGXO4ZusgPv8bl+PuG7bif71tF8I+qcitdHIuiedOLQIoLKBedi2peh4CAe6+YSs2DQTxkf88WFNqbi3K4V2vHcd9/+VibB8JAzANqyIKOLuYhm7kkdYM7Brtwct/8WZctWWg6DUUR8yBGQWmQLx8PTkcN7hxaHPYohb2SQgqZqO941aQlsUcRIHgtos34E/f+hrcdvEGhHxS0e75M0+cwJ988yUAQM6qBfDyTtesIzDrEP7o1gtx9HwSz1rGrRp2nUOFmINAgMGQgv969QUgxPQdiQLBxv4AJhcz9rWO+GUIzLfkQJEE+/r5uHLgtDncOLQxPz4yi9NWgDloGYKx/oDtGmJupVKCiliUWpnWdHs3bSsHD/vI1ZxhL747rB3+Ug1VyKx9Rtol3pK00oGZUXCysS+A6XjGdhlVuq6KWPg42cpB5jEHTnvCu7K2KUae4rcf2AdJNBezsNUkbmNfAC9PxgGYriY3wiXKQdXzZbEGL7tBVD1vL8RBq622m6uolKoxB0dFeSkDIQWT0XTBOFR4nM8xEKg0W8nL15PDcYMrhzYlmTUrcVk76ZBDOQCmAXBzfbDHOhfIbM6wd9WFmIN3d7qqnrd35CxwnK6hjsA2Dm5uJas/lRsDIQWLKc1OY2WxnFJclYP13csBfg7HDW4cPM5cQsXrPvZ4WT4/q3CWLAMQtFQCS6+s5PoAUBaQdiqHQszBuztdcxiPueiygTzZmpSD+d7cVEaiinLoDypYyur2TIyKbiWp8HFiysFOZeXGgdNmdJ1xeGUqjgd+frptGs99c/8kJqMZHJiKFx1nxuGPb70Q977lIvQGzN3sWL+ZellpoQNMV4xzgVR1AzmDglJaFHOIp3N4eTJW1/dTD9RcobGdrRxqMA5aFbdSSq3iVrKm7LFah0rGwSe5KQfuVuK0J8saB0LIlwghs4SQVxzHPkwImSKEvGh9vdXxu3sJIccJIUcIIbc6jl9BCDlg/e5TxIr8EUJ8hJCvW8efJYRsru9bLOZnx+fxFw8dbIudXD5P8bVfnAWAsnkMLGXy0vE+/M4bttnHWTprNeUQKlEOzDWlGQUFoep5/PPPT+Mdn3/aNhheQTMKxkESBSiigHTOvW7DSbVU1pRqVDYOwVLjUMGt5KIcfFw5cNqUWpTDPwO4zeX4Jyile62v7wMAIWQXgDsA7Lae81lCCIvSfQ7A3QB2WF/sNd8DIEop3Q7gEwA+tsr3UhOS5RfWDe9WADOePrmAMwvmgrRQYhyYcmCKgcGMQ7jCAgYUAtKsaR3b1eYMWlTnEMto0PQ8zsWzFV+rFZjKoRD8DShibW4lvXLMgWUrudEfMq8l+79YWbYS663EjQOnvVjWOFBKfwpg+SRyk9sBfI1SqlJKTwE4DuAqQsgGAD2U0qepuSJ9GcDbHc95wPr5mwBuJm75hHVCtrJ7cnnvf1j/bd8EegMytg2HsFiqHCzj0FNiHHr8Mnr8UlXlEPSJyNPCbpYtXJqeR04v1DlkrSCv2zyDVqLqhVRWoNxNVgk75lAhlZVlfJUyGDJHgJ5ZSEMWSZH7yAlTDgIp3GfcrcRpV9YSc/gAIeRly+3Ubx3bCGDC8ZhJ69hG6+fS40XPoZTqAOIABt3+ICHkbkLIPkLIvrm5uVWdtGzt7rzmKilF0/N47NAsbt29DqN9gTLjUEk5AMDv37wD77h8rOJrM/cJcy0xI5Az8kUxB7bgTkXLJ6G1EucwHcCMOyyXrUQptV1mpU0H2YjQcAWDypTDVCyDiF92rYUACobAJ4n2YySBQCDcrcRpP1ZrHD4HYBuAvQBmAHzcOu72qaFVjld7TvlBSu+nlF5JKb1yeHh4ZWdswbJ7vO5W+vmJeSRUHbfuXm+nUjpZyuYgCgQhpXy3+9uv34qbLhqp+Nqs/oEFZtnCpenOmIOBDDMOMS8ah5W5lXSrkaAiClD1fFFCgqrnYeRpZbeSFXMw8rSqImPKwe9QNYQQc3woNw6cNmNVxoFSep5SalBK8wD+EcBV1q8mAYw7HjoGYNo6PuZyvOg5hBAJQC9qd2OtmHZRDo8cPI+QIuL67UPoDyplAel4Jocev3tF73KEHMqBUlowDk7loOftduCTXlMOOaNIOdTiVmLvqzdoqgCn0mDFbZUC0rIo2FP1ajEOTsPF/s2VA6fdWJVxsGIIjP8CgGUyPQTgDisDaQvMwPNzlNIZAAlCyDVWPOHdAL7reM5d1s/vAPA4beBMRWYcdA+3pDbyFI++eg43XjQCvyxiMKQgoepFfut4Rnd1KdVCyFdI/3QuWpqeL+qtlPGyW0l2ZgYt71ZisZQ+65o54w4ptbpxAAptz6s9hgWkncoBMN1NPObAaTeWbZ9BCPkqgBsBDBFCJgH8BYAbCSF7Ybp/TgP4HQCglB4khHwDwKsAdADvp5SyT8X7YGY+BQA8bH0BwD8BeJAQchymYrijHm+sEqzdhJdl/umFFOaTGm7cabrOWJ59NJXD+l5zYV/K5MqC0bXiVA5O45DNGfYcBzVn2MrBi24lRSzszoOKiNkltepzmLuMuYicGUtJx0yMSvSHFJxeSFdMYwWqKAdZ4NlKnLZjWeNAKb3T5fA/VXn8fQDuczm+D8DFLsezAN653HnUC5ZF4mXlwHbsTBmwPPvFlIb1veZshngmt2rlwHa/qRI1UlwYV1AO07EMjDyFWKEdR7Mpz1aSlq1zKHMrOZSD3W21inEYtJRDNbeSzyXmYB7nbiVO+9F1FdK2W8nDMQe2kMnWYsNcGs6g9FI2h54qu9hqhJzGwbGjLW3Gl8kZEIhpSGcT3qh1MPLmpDqfVOxWytQYc2BuJadySNWiHCwDXe2aV445cLcSp/3oOuMgCd6f6ctcXsyHPWi5lRZSBdfJmtxKVoZTSjWKFi3ngqnqBtKagU3WJDSvxB3YtXEuwEFlBcbBUg5OQ1iLW2mgJuVQ3Ka7cFzgyoHTdnSdcbDdSh5OZWVBYbYTZbtWlrFEKcXSmgLSBeWQLVIODrdSzlQO20ciALyTsVQwDiXZSjkD1fIYNBaQtq7lvjNR/J/vHwKl1K4Z6QmszThUUg4KNw6cNqTrjIPdPsPDFdK2W8k6176gAkIKbqVszqxHqLaYVUMWBSiSgKRWHHNw7qbTmgFNz9vjMr0SlGbnW5qtRGn1QjM75mAZ1C/85ATu/+lJnF9ScS6ehSQQuxLajX47W6mKW6lithKPOXDaj64zDnb7DA8rB7aQsHMVBYL+oGL3V6pWHV0rYZ+EtGoUxxwcrplYxvxbAyGzHcfskjdiDmoFtxKAqq6lUrcSy0eYjmdwLp7Fuh5/1YA7SwpYXZ2DALWGeRMcjpfoQuPg/SI4dm5O10l/ULbnCbCOrKsNSAPWqNCSVFancoilzb8RkEUMhX2YT2plr9EKbOVQ4lYCqg/80eyAtLnIs13+dCyD6XgGG6wssEpsHjJjL2yYkhtuFdKA2XzPy6nTHI4bXWcc2qF9hqYXu5UAs/nbQrK+yiGp6nZfJaCCcVAkyzhUryNoFixGUpqtBAAZl26rDKYUA4qIa7cO4t63XgTANA7n4lk7RbgS20ciePZPb8Zlm/orPoYZHPdspc4wDqfmU/jEo0erxnc47rTbNes649BOysE5H2AgpBSUQx2MQ8gnIaWVKAfLLUNIwa0UkEUMRRTPGAfbrSQ73Uqmqyejlf+fnl1I44kjs3a7bkUU8NW7r8F/u34LIn4JU9EMZuJZjPZVVgSMdT3VDUjlOofOSWX9/oEZfPKxY1jKLD8/g2NCKcVHHz6E6/768bZSkF1sHLxrxUsD0oCZznpqPoUPP3QQJ+aSAMrbda8E061kuAakw4pk79CDimiqllQbuJVclMPf/fAIPvjVFxy1I4W4wmhvAAenl6DqeaxfZuGvhdK50c7jnVIhbXfy7RBj1wz+z/cP4Qs/OYmZeLYoHd3rrC7dpY2R7App735YVRe30ntetwVLWR1fefasfe5rUQ49fhnTsUxJKqtlHPwSEtbPfivmEEvnkDPyRefUCti1UVzcSm4xh32nF5FQdTvm4Dz/0T4/fnZiwf55rQQUEQMhBZsGi1WIIglQPaxUV0KqpM07pzrZnIEvPnUKG/sCmIplEE3lsKF3eZXqBbpPOQjtoBzMc3PujrcOh/HpOy/Dg++5yl7gqmXOLEdPQEI8o7sWwTmbywUV0S7CK20b3gpUl5gDUw6lbbunYhlMx7OgtOCKU4qMQ8CW+evr8IFVJAFP3/tGvH3vxqLjPqtld7v5nN0ozADpDGPXaCYW06AUuGHnEICCu7Yd6D7jILGAtHdvbreANOPqrYP4f//ttfiDN+1c0y6+JyBjKZuzF1tRIHa/IefQm4BiKgcAmEu0XhIX3Erlqaylbbv3nS50fo9aAXa5xDgwlstWqhXnoB/7mNw5c6SZcshw5VATbLTspWN9AIC4dR+2A93nVhLaIyAtCqRi3v1124Zw3bahNf2N3oAMTc8jnjGHBgVl0d4VOpVDQBYxHDGVgxeC0qpLhXSggltp3+mo/XPMNg6Fa7rRMg6SQGwD2AiYIVP1vO0Ca1dKpwdyqnN6IQUAuHTcNA7RNjIO3acc2qAIzvTtN7YDKotXnE+o8EkCZEmwd96REuXAKocXPFDrUMhWchiHCm6lX5xetA1szMr0kqVy5bBcAdxa6aQ50klLXXLjUBtnF9OI+CVsGQoB4G4lT0MIgSQQzwekGx34ZQV0s0tZ+GURiii4xhzMVFbTOHhBObg13rOVg8M4JLI5HDmfwBVWXQJLA3bGHJgrqV4upUrYxqED/PQpHnNYEWcW0rhgMAi/LMIvC7aCbQe6zjgAZsaS15WD023SCJhymEsy5UDA4qXO/kEBWURIEeGThJams+pGHs+eXHBNZZVEAYooFM10OHIuAUqBa7YOAHCPOazv9YMQLFsAt1YUqfNiDp2ggprBmYUULhgwVUNfQLEVbDvQlcZBFgRPxxy0JigH2zgsmcbBuaNmAWmfJEAQCAgxffLzLQxIP3Z4Fu+6/xm8PBG3z81JQBGL3EpHzicAAJddYCqHWFqDQFDkPpJFAb9y6SjeeNFIQ8+dqZx2KoCqBI851I5u5DEZzWDToNl6pS8ot1XMoesC0oDpd/Zy+4xm1BMw45BQdWws6RfEJqKxLCAAGIr4MNdCtxLL8jgwFYciCWUZQUFFLHIrHT2XQEgRscPqKhtN51yv6SfvuKyBZ21SyFZq7wWVUsrdSitgJp6FnqfY7DAO7ZSt1JXKwesxh5xBi4q8GoGzutonCUU7caYcAo7MmqGQ0tKANKvInYplXF1uAVksylY6cj6BnesjiFgusqVsrkgdNRNfh7iVsrm83c2WK4flYZlKmxxupSh3K3kbWRTs4S9epDkB6YJo9Eli0d9jw4D8TuXQ4uZ7zsWotD0FUO5WOno+iQvXRRDymY+ltDhTqZk4U1nbmaSjMSNXDsvDahwusJRDf0hGLMOVg6eRRK8rh3zDlYMkCnZWkk8W7L8niwR+6+dit5KCxZSGfL41RtW5GLkpB6dbaT6pYjGlYee6CCSxoIoanR5ciUK2Unvvtp1de3lvperMJ1X845MnMRT22X27eq2AdLtUynelcZBFb8ccND0PpQkLGYs7OJWDLAp2x1OnW6k/qEDPU7vnUrPJFCmH8tvW73ArHT1nBqMvXG+OOGVGsFV9ofwdUiHtVA7LzezudJ4+sYC/feQwHjl4rqzbAqUUd395H84vZfGF37wCgpUE0R+UkTNoWSW/V+lK4yAJxG7E5kWa1eCOFbsVK4fCTjugOHsssbbYrbmxnW4lN1UVUiR7ngPLVNq5zjQOQcu11KqYgyJ2nlup3YPra+WzPz6OzzxxAr/z4H48cvB80e+WMjqePxvD+2/cjisuKMz/YFMI2yXu0JXGwVQO3v2gNsOtBDiVQyGVtcg4yM44ROW22M2gyK3k0oIi6DNbkAPAsdkk+oMyhqyGgSGltcqBZSu1eyprisccbKZiGbx+h9nC5sxiquh3SeszMtJT3JKlzxo12y6FcF1qHAj0FvnOa6EZAWmgYBz8slgUc2AB1KBSXCkNlDe3axZqzrCNlptbKaRItuGKpjQMhX12uqvtVpJaHHNo8902Uw6iQLo6W4lSiplYFjvXRdAbkDETK56vnsyyTgPFLfX7rM9bvE2C0l1pHCTR20VwOSPfFBdIj0M5sGCtGXNgE80KO3SWwZRqUcwhqxsYHwgiaFVrlxLySbZySKq6fb7sd0ALlUOHZCux6zsQUrraOMQzOWRyBjb0+rGh14+ZeKbo98yIMrXN6A+ZyoG7lTyM7LH2GT89Oocnj83Z/9aa7lYqVQ7MreTSFrtFi0JGMxCQRewd78Ooy+yFkCJCM/LQ9DySql7UH4p9SFtlHOz2GW3uimEbg8GQ0tVupWlLKYz2BTDaF7D/zWDGoXTeClMObm6lD33jRXz0+4cacbqrpisrpCVBgG54YwZuUtXxe199AUFFxM/+vzdCEAhyOm1K2qUz5qDqzpgDcyuVz2lmMx+aTTaXh18W8E93vda1gypTB2lNR0rVMRIp+HtZzKFVAWlRIJBF0jFupcGw0jYZN42AKYUNvX6s7/XjhbPRot+nbOVQvLz2BplxKFcO+05H7RiZV+hS5SC0RDmkNb0sKPm1584inslhJp7FCxPmTdbsgLQz5qBIBbdSQClXDqlWBaR1A35ZREAR3bOVfOz8DKRUo4JbqTUxB8A0TO3vVtIRVEQEZKm7lUPcoRx6/Yimc0VZfIWYQ7Fx8EkigoqIeZdOA9GUZrenoZQiltYwFcu0tCaiS41Da4rg7rz/GXz04YJ01PQ8vvjkKVw63gdFFPC9l8/Zx5vhAukJFBoWXxW5AAAgAElEQVTsObOVIj4Jv//G7bh193r7sWyBbV0qa/VBOQVlo5e5lVpd5wAAg2EfJhbTLfv79SClmbEcvyy0fUHfWpiJZewBUWwetDPu4DY0i3HVlgF854WpIvWQM/JIqDrmE+axe799AHs/8iiu/+vH8cJErJFvpSpdaRykFimHU/MpvOT4z37u1CLOLWXxgZu244adQ3j4lRnk89SMOTQxW8knO42D2YX1Q2++ENutpnWAB5RDzqhqHNgHMamabiWncmB1Dq1qnwEA120bxNMnFjydCLEcSdVA2CfBL4tdHZCeiWftAVEb+vz2MUayglsJAO55y0VIZHP49OPH7WMsQJ3JGUipOvafKbiYWjmatyuNgxmQbu6HNG9VF5+aL+REs6lQmweDeOslGzATz+KV6XgL6hxEe+GstLv2SQIE0sqYg2G39XCDGa9oWoOep67KoVUxBwC4YecwEqqOF1u4E1wrptE1h9Z08wzp6VgGo5ZRYMkR07GCckipupUBWH6/XbS+B79+5Ti+/PRpu0NrNFUIUM8nVZyLZ+3q/lZuJrrTOAjNb5+RUHVQaraOZpLS9k36JftmmFjMIE+b4wLZ2BeELBKM9vnthbPSAkoIsWoJWmgcqigHtkubXTJ3WiFHvKRQBNe6mMP124YgEDMzrV1JZnWEFAl+SezqmMNMPGu7k9igKKdySKh6WaaSkzfsHEbOoJiyDIoztfXUfAoJVbc7ubaycLIrjYPUAuWw5Ch8YerBKT9Za+lFNuu4CcZhfa8f+/7sFly7dXBZ5QCYAepWVkg7A+Sl2MbBkuFeqnMAzEyVveN97W0crFiOXxaR1Y22aSBXT/J5ipl4xnYn+WURgyGlyDiUujVLGQybmXSL1mTFqGPC4itT5jAr1smVG4cmI7egCM5ZFcn6vCcs5RBSJHuGwqKVydAMtxJgLlqEEPhYzKHK3w35JKRaoBwopcgs41ZiSuH8kvkh9VpAGjBdSy9PxdumQrYUZ0CaUni6P1mjmE+pyBm0qNZmQ19xIVwyq7sGoxmDVjxhIWVuZBYdyuHlScs4DJjGgbuVmkwr2mcUKYc50zikVB0hRYQoEDsVk0nMZnRldcJaS1RzvQQV0W5u10xYCqhbTyVGNeXAAtKNnsu9HDtGIqC02D/dTrAdMXPvdaNr6XzcvL+cc8c39AaKWmiUVuiXMmhVSrOUVjflwEaLtjL9uSuNg9SClt1st0gIcMoaApJUdce8ZjN/n0nNZikHBuscWi1oG1QKze2aCassDlQxDux3s5ZyCHlQOQxbhXmtHJq0Fky3kmgbh25MZ01kzc9xr2OS4lDYhwXHAp9UdXvUrhs9fhmSQLBg3QfRdA4hRcRASMF0PAtCgLF+phw8XOdACPkSIWSWEPKK49gAIeRRQsgx63u/43f3EkKOE0KOEEJudRy/ghBywPrdp4jVFY0Q4iOEfN06/iwhZHN932I5stWyu5k+0yXrpto2HMap+SQAM3DlXMQiPsk2Ds1eyJy9lSoRdDS3ayZssEy1gLQgEIQU0VYOYY/FHICCcWhleuJqSao6srk8BkK+rlYOzK0acjSl7PFL9ucbWD7mIAgEA46xu9GUhv6QYqevDoV99v3r9ZjDPwO4reTYPQAeo5TuAPCY9W8QQnYBuAPAbus5nyWEsE/05wDcDWCH9cVe8z0AopTS7QA+AeBjq30ztcIWCaOJriWmHPaO9+H0fNoe1u7cYYT9rTMOSg0B6ZBPbEm2Eiu888vVr0nQJ9kLr7PpWdj6IEstzFYCYH/421E5nLN86qN9fvv/4fhcAh/5j1eb+jlqNWxzFHTcXz0BGZqet2s/nB6BSgw61MZiWkN/ULE3Dxt6zRoKgXg85kAp/SmAxZLDtwN4wPr5AQBvdxz/GqVUpZSeAnAcwFWEkA0AeiilT1Nzu/7lkuew1/omgJuZqmgUkrUANjPusJTRIRDg4tEeJFUdc0nVDFz5i90frXMrsYB05UsfkFuTylqLcgDM68f+T53KoScg4XffsA237FrXuJOsgbAVzG1H5cCay23oDcBv9d761v4pfOlnp3C2zSu/VwJzqzqVA0tbZQkmpRX6bgyGFDsgHU3nLOVgGgc2VlSRhJYG/Ve7Aq2jlM4AgPV9xDq+EcCE43GT1rGN1s+lx4ueQynVAcQBDK7yvGqCuVCaeeHjmRx6ArLtSzwXz5qBK6XEONgB6dYoh2p/N+QTW1Ihna0h5gAUNwp0ynpCCO55y0X2ZLhWQYjZcqEdjYOz2Rwz0ofPLQEopGR2A27KoWAcctCNPLK5/PLGIVzsVhoIyrZxGO0zM6FkUfC8W2kluG07aZXj1Z5T/uKE3E0I2UcI2Tc3t/p8ceY6aWZQOp7JoTcg2z3dF1JamfyM+CX7Zmi2cnDOkK5EsEVFcEyu+5ZxK9ndVytUp3qB4YjPtfFaq5hLqHYxVjVmrEDpup6CW+m0lVjRXcbBvBeDjo1Kj98MTieyekFZLKscfIWAtBVzYG4llgnla1PlcN5yFcH6PmsdnwQw7njcGIBp6/iYy/Gi5xBCJAC9KHdjAQAopfdTSq+klF45PDy8ylMv+J6bOSp0KZtDj1+209iilnGIuGTVAN6MOQQVEZqeb7oflLVqWM6txOIMy+3aWonXlMOH/+MgPvCV55d93Ewsi6GwD4ok2P8PLNYQ7SLjkNLM1hiS43MSsYzDUjaHhGrGFqtlKwGmckhpBhLZHBKqjv5gwa20wTIOsigg14bK4SEAd1k/3wXgu47jd1gZSFtgBp6fs1xPCULINVY84d0lz2Gv9Q4Aj9MGpxHJgvm2c00OSDuVw2JKM9sRlASk7XNsdp2Do/FeJeyBP01WDyxlkvm6KxG0rqWXjYOpHLxjHCajmZrOZzqesRetUiO92CaTzepBuqQdPFDobrwy5WCuAyesmqf+kIJNVuHb5kGzdUarYw7LfooIIV8FcCOAIULIJIC/APDXAL5BCHkPgLMA3gkAlNKDhJBvAHgVgA7g/ZRStpK8D2bmUwDAw9YXAPwTgAcJIcdhKoY76vLOqsCUQzOt8lImh9HeAHr8EmSRYDqWNRvEFQWkC7nTzXYr+RzzHCrhbNvtzPNuNHbMoUr7DKCQlbTcB7OVDIV9WExr0I180e6zVSwk1ZqM/bl4FluHzUWrNGus25RDsOQ+tJVDJoekpRxqyVYCgGPnEwCAgaCC127ux/d///XYNdoDoDWdHJws+ymilN5Z4Vc3V3j8fQDuczm+D8DFLsezsIxLs7BjDk2c6RDP6OgJSCCEoD+o2BkeTvnpbNbV7IB0T0CGXxYwEvFXfEyr2nZnc7WmsjK3UnUj0kqGIz5QairHkZ7K17pZLCQ1GDUI9Zl4FtdvHwJQrhwWusg4pFWjKIkEKM5WSlrKYbl7kLXQOD5r1jz1W21smGEAzDWgkwLSbQFznTSr+pBSiiUrWwkwB7SfXTTlZGkqa+Ecm/tf0xuQ8fQ9N+PNVdI92UCdo+cSuO97ryLfJLdcpka3UqgNlMOwtWOc9UDcIa3pyOSMZeNIiWwOSVW321Q7jcP6Hn/3KYeShT+sSCDEuk72FLjqynooZN4HvzhthleZknCiSK2dHtiVxkFiMYcmSTZVz0Mz8rYrxjQOpnIoTWVlNNutBJh+T8FlPjODKYcvP30G//jkKUxGm9MjiLmVlg9It4FxiHinEG7BkTXlNqdD1Q3c9Hc/xid/dAwAsN5qNscaIIZ9EraPhLsr5qCVKwdBIAj7JCxldcf86NqUw/NnY3jNhh7scAzWYigtdit1pXFgnUebpRxYdTRLeesPKfaCV6Qc/K1TDrXAjMPz1kB1Z8uARmKnsi5jMO1sJcXDxiFs7r69kLE05zBQSRdX4cHpJZyaT+GLT50CAIxaAWlJFCAJBJsGghgIKV2VysrmaJfS45etbCXzOkaWUQ5BRbTv5w/evN11U6ZIrXUrefdT1EBkobmprKwjK1MOLFMBKL6JIi1WDsvBduRM6rKK0EaT1Q1zEl0VVQO0h1tpyFIOcx5TDmzH6+SFs8VT6zb0FdpU+2URFwx2n3FIa+XZSoAZd0isQDkQQrCux4+gIuLNu9a7PsacWNm61iTe/RQ1kGa3z4iXGIcBh3Eo6gHUwoB0LZRWKDdNOWjVp8AxQm0QkA4qEkKKaA+TbyULDgPlZhxenIhhtNePyy/ox6OvnsdIpOAXf/tlo7h26xBOzCWRyOrIGXlPqt16k3bJVgIs5ZAxYzN+WagpE+2Td+zFQBVXLlcOLUBqYvsMVTfsitgeF+NQOSDd2iZxbpTumJqmHHL5ZVtnAIWAuZeVAwCM9PjtoUStxJll5NaK/YWzUezd1IeP/doevO/GVNHi/1dvvwQA8OAzZwCYc0iqZbp1CimXOgfArHWYjmVxLp61i9mW47JN/VV/3+pU1s439S4oTWyf8dZPPon3WxWobsrB6VZihkIg8EQOfCmlO6ZEnZTDMycX8GffOVDWQj2RzeH7B2aQ1Y1l01iB9ghIA8BYfwCT0dY3q3MGxZMlymEuoWIymsFl4/0I+STsHu11fY2BYKGos9Mx8uZEQjflELFiDqcXUtgyFKrL3+PZSi2gWe0zKKU4u5jGxr4Artk6YFeYMuMgkOLcfWYovCrPfZIAgcC++Zcy9VEOTxyexb8+e7as59C/vzCF//Gvz2Pf6WhNbqUey7hWG+7uBcb6g5hoUqZXNRaSGiTLpVHqVnpxwow37N3UV/U1+kPW7PMuMA4spbo0WwkoxBxOzdXROHDl0HxYKmuj3UqqnkfOoHjXa8fxtbuvtRc4ZhzCPrMojuGXBYgC8WS8ATCDaBG/jL3jfQgpYt2UA/vQHbWqRRlsaPtULFN1RChj+0gYH/3VS1remns5xgcCZvsUFz9/M5lPqtjYbwaZS4c4PX82ClEguLiCYmAMWvn60VR7zsVeCWm1vCMro8cvI54xs5XqqRzasfFeW9MstxL78JfuZJlxYGX3DELMfGkvZioxPn3nZfjQLTttGV0P2DCfI+eKjcP5pYLbI1CDW4kQgjuv2mTHHrwK66Ez0eI5CAtJzT6XZEnM4bFD53HFBf3LtiyxlUMX1Dq4TYFjOD/j9TIO7dp4r62x3UoNbp9RqJYsvpn6LT+tW7pb2Cd51q0EADfsHMb4QBA9AaluAel0BeUwmygEbWtxK7UL4/0eMQ4pFWP9QQik2K10ci6Jo+eTeMvF7imWTti9vOihNuSNgl0j12wlR68xrhzaGKlJ7TOYcig1DrIooMcvuXYPjfi9rRwY9VQOWaYcSo3DkmqnTy7XOqOdGGfKoYVxByNPsZjSMBxWEPJJRS6uRw6eBwC8effyxkEWBUT8EqJdoBxYg8JKdQ6AmWW40VEPshbMmANt6qx7J95fhRoAcys1OtjDdtZuHRoHwz6E/eVVlKZy8F4aayksAFcP7JjDuUTRB+F8Ios37VqH/qBst0XuBPqDMsI+qaXKIZbWkKfmfRhSpKKYww8OnsOesd6aF7mRiK+mgUHtDms4WSlbCTANf70yDdkm0ake4ukc3vbpp/DwgZm6/I1qdM4nbgVIzY45uJTS/9Z1m+3ZDk56ArK9WHqZHr+MMwv1WdzYjiylGZiKZTDWH0Q2ZyCWzmG0149/+e2rbfdFJ0AIwVh/oKXGgdU4DIYVc/yrFXOIpjS8NBHDH96ys+bX2j3aazeQ62TSVWY1sEy5rXVyKQHOTSwF+5NJTceBqXhTClC70zhY6Xu5BsccmI/STTncdd1m1+d86JadrtWqXiPil+y2IGslmzMwFFYwn9Rw9HwCY/1Bu/fQSI+/Yo59OzM+EMSZhVTL/j5rnTEQUhB2uJXOW3GercPljeAqsWesFw+9NI3ZRLajC+FqUQ5sUE89sGfd63nAqqtjyRvNiMF1pVupWTOkEzX2WXFy8cZeXL11sFGnVDd6AjISWb0u/tBMzsCeMTOf/sWJOIBCMNrZsqGTGO8PYmIx0zJ/MosRDIZ8CCqSvSFhgWWWhVQLl46b/3cHJuN1PktvwVJZ3bKVhiM++GUBe8ar14WsBMWKszlbaLAmlM3IyOtK4yAKBAJpfMyBZSst16GxHYn4JWhGvi4VnGnNwLoeH27YOYz/99QpzCaydhrrOg8MxGkEmwYCyOSMssK/ZsGK1vqDMkI+yU7TXHQYjVrZPdoDgQAvdbhxYNfIrc6hNyDjmXtvxtv2bKjb3yvMnSl8xpgLtpZ2MmulK40DYMYdGp+tlIMokJpaP7QbzqHqayWrGQjIEj78tl1Q9Tzu+94hu/dQpxoH9r4a2br7z//9AP7F6n1UChvQ0xdUEPaJBeWQKribaiWoSNgxEsHLk7HlH9zGpDUdUpUi1b6gUlTUulZYQNq5AWPxyIDS+DWl81atGpEF0vD2GcmsXlYF3SmwAFw9WmhkcgYCioCtw2G894Yt+O6L03jm5AJkkaA/2HmqCyjMw25k8sGPXp3Fj4/Muv5uMa0hYhVcBn0FtxKLRaz0uu8Z68WByXjL3GTNIKWafZWa9Xl2y6rkMYcm0B9S8NzpxYYaiISqu9YydAJscNFaW2hoeh56ntoy+d3XboZAgB++eh4jEX9HGlag4BbINtA4qLpRcb5zNKXZ2XLOgPRiSkNfUF5xOuaesV4spDRMx1vfbbZRpDW9qU0d7VRWl5gDdys1kHvechFenozjcz8+0bC/kczqnm8Ct1pY3cFaax0KMtl8vXU9frx+xzAoBUZ6OjMYDRQCisyH3AiyuXzFhniL6ZxtHEKKBFXPQzfyWExrdqfVlbDBGiG64IEhRo0ipbl3ZG0Uspty4AHpxvPLe0bxtktH8cnHjtnDeOpNSutc5VCvmIPbTugdV4wB6NxMJaDgM26UW4lSClU3Kra1iKY023XEsulSmvn4lcQbGCxI6zYXolNIq61XDjwg3SRuvmgEep42bNh7Mtvcm6mZMEW0VuVg3+yOANstu9ZhKKysKNe+3WBKKeMyu7ke5AyKPDVdm6pevmAvpgoKgd2jKVU3j6/COIRsJeT9Gp3VksjqrmmsjcKtQpptpvw8IN1Y2AKXXOMCpxt5/MYXn8XPjs8XHU+oumsBXCfAYg5rLYTL2DuhwnXyyyIe+Z834IM371jTa3uZoLXzyzTIrZR1GAQ311IsXYg5OI3DQkrDYHgVxsGhPjqVxZSGgVVcm9XCAtJO5ZDRDAikOWOEu9o4MJfPWvvqL6Y1PHV8viwzJJnVEelQ5RBURIgCqWPMoVgmD4Z9HdWJtRT2ftMNciupucKCslDiWsrmDKQ0wzFXxDyXhKojml6lcrDu83QbVPevlsW0hsFVXJvVwpSDM+U+kzMQkJuTMdXdxsF2jaxt9xtPm88/W9IrJ9nB2Urm4B9pzTGHTBN9qF7CJwkgpIHKIVesHI6cS+D0vNmuI2bdr6xfFQtunotnYeTpqvpYsdfoNOXw5/9+AP/3R0ehG3nE0rmm9vhiAWnNKFxTM+27OWtKVxuHQjrm2nY7LKB9drHQmdLIU6Q1o2PdSkB9OrMWsi+6yzgQQhCQxYYZB2fh1GJKwwe/9gJ+7XM/x8Ri2lHoZt7/bAPDGgGuxq3E/v/aoS/YSnji8Bx+enQOUcugrubarBZbOegO5aAZTSmAA7rcONTLrcR2YpOLabsIqNIsh04i5OjJs1qYcehkF1IlgorYMLeSUznMJVScmEtiIaXhvV/eh2mrvXaftQtm1drPnTI7qw6soHUGQxYFKJJgN6frBCilmEuomE2odi+q1bjcVgtrn6GWFME1S2V3t3GoU8YNUw4JVbcNRaURoZ1EyCetOU+fZessN46yE/HLoj3oaDXkjDxiFYbsOJXDS5Mx5AyKW3evw+FzCXxz/ySAwkI3HPFh92gPfnx0DgBW7VcPKaLd1roTiGdy0Iw8ZpdUO6NxNTUgq8UnljfeYzGHZtDVxkEWBfhlYe3KwZGxMxE1pXnK7sjaucYhqIhrvnbMrRLsVuWwBuPwwM9P440f/wmMfHnLCtWhHPadjgIA/vv1WxDxS3jcSpxw+s9vunDEfp3V7o7NBn6doxxmrb5XmpHH6Xnzc93MbCVZKm+8l8kZTVPZXW0cACDsk+umHIBCUDpRYX50JxH2SWvOa89YWTXdqBwCsrimIrjTCykspjTXOh1nKus5q4nhjnURXL9tyN6J9jn6J9100bD986qNgyJ1lHKYXSpc18PnlgA0163klsqazTWvSrvrjYMZVF1rtpIGnxU8YsahG9xK5hyAtbuVCIF9/bqJgLK2gHQ0Zd630y4jOlkqK6uC7gvKGAgpeP3OIQBm40TZkSu/d7wffUEZQUVc9c406BM7TDkU+kQdnjHnmzczW0kUCAgpb9ndrI1U930iS4j4pTW7RuKZHNb1+DEYUjBhZSwlbeXQmV1FATM/fq2LQTPztr3GWpUDC5LOuDS7Y8ph1JoDzcZX3rDDVAilO2BRILh11/o1TTKrR4KCl5h1tFM/dG6pzKA2GkLM9uClRXDNcit17ra2RsI+ac0V0vFMDr0Bc2c2YSsHc1fXyamszlbPq6WZATavEVQkpLXVz5FmKanVlMNoXwAHp5fsViTjA0FsHgy6zi//y9t3F7VqWClBRWxYK5pWMLukmouzkUciq2PzYLDp56BIQln7DB6QbhJh39pz9WOZHPqCMjYNBMtjDk3sxdJsQoqInEGLdjYrJd3EnZDXCCgisrnVXzumHKZjWTz00jSuuu9Hdgor+z7aa6apbh0uKIKP/doe3HPbRWWv55dFu/ZnNYQ7LiCdxWif3+5y0Mx4A6NMOfCYQ/OI+OW6uJV6AjIuGAxiKpZBNmfgxFwSfUHZbm3didgtE9awIDQzwOY1ArK46mtHKbULs2biGTx5dA6zCdVWrlm9oBwAYOtQoYnh1VsHGzKnPOjrrFTW2YSKkYgfw1br+NXUf6wVRRLsmAOllKeyNpP6BKRz6AvI2DPWByNP8cpUHC9NxHHJxt6O9qWzDpVrMa6ZJgbYvMZaUlnTmmHvKKfjWRyysmmYcmVupd2jvZBFgt2jPXU44+qElM5SDnMJFcM9Prt1PKsobyayQzmoeh6UAn6uHJoDm4K12vGGlFI75nDZpj4AwNMnFnD0fAJ7xnrreaqeo6AcVr9b7Ga3kl8Woep55F3qFJaDuZRkkWByMY2j55MACsYhqxuQRYLX7RjCvj+/BeMDjfeXBxUJ2Vy+4eN3m8XsUhYjEZ9dQd465WDeH83uQ7Ym40AIOU0IOUAIeZEQss86NkAIeZQQcsz63u94/L2EkOOEkCOEkFsdx6+wXuc4IeRTpInb7YhfQp6ufoFLaQb0PEVfUMZQ2IdNA0F87RcT0PMUl2zsq/PZegs24GUtyqGb3UrBNcyRZmmsO9dFsJDS7N2lUzn4JPP1ewPN2fGytt2NagnSTFKqjpRmYCTit5VDMzuyMmRRsKvdMy6DsRpJPZTDTZTSvZTSK61/3wPgMUrpDgCPWf8GIWQXgDsA7AZwG4DPEkLYu/wcgLsB7LC+bqvDedUEyyZa7QLHCuDYB/CyTX2YsrJHOl45sAEva/Azd3O2UmAtxsFSDrs2FNxFfllwxBwM+OXmOgaCdbgfvAJLY3UqB7cMr0bjjDlUam/fKBpx99wO4AHr5wcAvN1x/GuUUpVSegrAcQBXEUI2AOihlD5NTd/Olx3PaTisgnm1GUustw0zDpdvMoXSUFjBBitTpFMJ1UE5pJvYSMxrBNYw8IcZBxZLUEQB124ddFUOzaIw8Kf94w6zVlX5uh4/hluoHBSR2KqwrdxKACiAHxJC9hNC7raOraOUzgCA9X3EOr4RwITjuZPWsY3Wz6XHyyCE3E0I2UcI2Tc3N7fGUzcptO2uPSg9m8ji+KxZMVlQDuaNw4xDpwejgfqMhszmujcgvRblwGocdo2a6nT7SBhbh8M4a3UGzuoGfE1WDvVQkl7hPFMOPT7sHe/DlqEQLtoQafp5tFI5rDXP8npK6TQhZATAo4SQw1Ue67ZS0irHyw9Sej+A+wHgyiuvXF0EuYTVuJX+9gdHsP9MFI//0Y32oB+mHC7aEMFwxIfrtw/V4/Q8jT1ecg0B6Wa2IPYaLOawmnhXNJ0DIbAXrNds6MGmgSCyuTzmkirUnNF05VCPGJRXYO65sf4AgoqEJ/7oxpachywKtlej2cphTcaBUjptfZ8lhHwHwFUAzhNCNlBKZyyXEZudOQlg3PH0MQDT1vExl+NNwZ7psAK30kw8a/skmXJgTcxkUcCTf3JTU2a8thrbjbDKxYBSinQXKwf/WtxKKQ29ARk9fhm//botuGXXOtvITCymoer5pscc6qEkvcLEYhpDYZ8dR2kVQ2EffnZ8Hn/98GFcstFUiZ7vykoICRFCIuxnAG8G8AqAhwDcZT3sLgDftX5+CMAdhBAfIWQLzMDzc5brKUEIucbKUnq34zkNJ7KKmQ4LKQ1JVTf76ZcEpAHzP08QOtulBMDqibT6ucF23nbXKgfz3svkVn79omnNbgL357+8C1dvHbTTVc8uppHNGU1vZliIObS/W+nsYhqbBgKtPg3c85aL8LZLR/H5n5zAV587C6B5UxPXYhbXAfiO5VeXAHyFUvoDQsgvAHyDEPIeAGcBvBMAKKUHCSHfAPAqAB3A+yml7C56H4B/BhAA8LD11RQiVmO8xAoWuMVUQTXEMznIIunKdExCCEKKhOQqfcyFcZXND/R5ATuVVVt5XYBpHIpTVMf6zcXs7EIGqp5v+nW16146wa0UTdvxw1YyFPbh4++8FI8dmsVzp81JfZ6POVBKTwK41OX4AoCbKzznPgD3uRzfB+Di1Z7LWrAzbmpUDpRSe1GLZ3KIpTX0BpSODz5XwqzyXd1i4EwX7EaY73g11y+aypVlw/llEcMRH6atFi7+ZscclLXHoJpNPJ3Dg8+cxmxCxa9ePoa9433IGXlMx7J4+97mN9pzg8M6E5cAABbgSURBVBCCyzb14cdHzCScdslWanskUUBQEWvOVkqoul2xGEvnsJDUWpLi5hXMZmurWwxYuuBIpLNTfivBdoDZVdY5uOXdD4V9WEipUPV807OV7AB7GymHv/req/i7Hx7Fg8+cwWeeOA4AmIllYeRpU6rKa2XveKGglrfsbiIr6cy6mCzM7F3K5LCQ0jDYxNGBXiPoE1cdkJ51pAt2IwXlULtx+Ld9E3jg6dOYT6plbiXArK+ZS2otUQ6yKECRBCTbJCA9n1Tx3Ren8RvXbEI2l8djh84jn6d2rch4v3eMw2WWi6uZg7G6XjkAZqHLdLy8J74bC6mCcYhlNCymtK71mQNsGtzqjQMhrSku8gKrMQ6Pvnoer0wtIWdQbOgtD5gOhX1YSKrI5pqvHACzjXu71Dn86zNnoRl5/NZ1W/Dazf2IpnM4OZ+058BvasH8hkrsHTOVQ7CJg7G4coDZ654NYV+ORadxSOewkFS7dnEDTNXlHKe4EuYSWQyGfJC6IO3XDUEg8EnCitxKJ+dTuGXXOvzeG7dj57ryoqyhsIL5pNqyLLBQG8x0OHIugb995AieOj6HN+wcxvaRMFhy4XOnopiIpiGLBOt7vOPu7A3K2DYcKppX32i681NZwrbhMKZiGWQ0A3/x3Vfw6ceOVXwsy1QCTFm6lNUxGO5OtwhgBaRVA6pufi1HPk/xiUePYjqWweySarcm6FZW0rZbN/I4s5DCtuEw9oz1uS7+Q2Efsrm8GXNowVzudhgV+q3nJ/HEkVm8bc8o/urtZh7MlqEQhsIK9p1exNnFNDb2BSB6LB399TuG7fkczYArBxSmZB2fTeLr+yYgiwLee8NW1w8fcyspkoDT86b87Ga3Emt5/ptffA6Hzy3h168cxx/fdmHF6tyT80l88rFjkARiDVPpbuOwkjnSE9EMcgYtmupWinOj0grl0BuQEUs3b3e7GuaTKtb3+PG37ywkWxJCcOUFA/jZiXkokrCmWdqN4k/f+hoYq2jvvlq4coCpHADgewdmkM2Z82JZ2lgpi0kNAVnESMSHE3NmD/1udisFFQmLKQ3PnV5E2Cfhi0+dwsuT8YqPn1g0YzuHzyUwm8hy46CINVdIn7Tut21VjMOQIzmiFcphIKQUuV69yHzSPYnkmq0DOL+kYnZJxa2717fgzKqjSEJTuwlw5QBTUhICfOcFs/9fUBHx0EtTuO3i8huEBaB7AzKOz1rGoYvdSiGfCN3azfzezTtw77cPYKmKX5QF+w5OxzGf1Lo2U4kx2hfAixMx5Iw85GViLyfnUgCKR36WMuS4F30tUA6DYQW/OO1t47CQVO023E7uuGoTNg+FcOXmAbutTjfDlQNM+b2xL4DzSyqGwgreecUYfnRo1rX2gaWu9gVlewhHN7uVWFXscMSH124eAFC9FQlraHZ6IQ0jT7u2xoHxW9dtxlQsg/94afl2YifnkxgIKVXnCjhjOP4WKIfBkILFtNZU98dKqVSb5JdF3HjhCDcMFtw4WGy1XEt7x/vxpl3roOl5HHBxjzDl0Bco3Fzd7FYKWTL3DTuH0cP6VFUJSDK3EqPb3Uo3XTiCC9dF8PmfnFh2XOiJ2RS2DlX3hTs3Kq1RDj5QWphz4jUopVhIqRjq8vuuFrhxsGB+3Ms29WHMKn45t1SeosmMQ4/VaE8USNPGMHoRphxuvHAYkRpmY0xE00U+8253KwkCwXtv2Iqj55N4aTJW9bEn55NVg9GAWYjGOgS3Qjkw47Tg0bjDUsbscNDNG7pa4cbBggWlL9vUZ+c3z8TLjcNCyqxrYB/A/qDSFR1YK3HN1kG844ox3HThCPyyAEkgy7qVrt02aEv3bncrAYXWCKcXUq6/n4ym8eGHDmI+qdkKtxos7tAS5cCMQ9KbxmHeSkUf6uI4Ya1w55rFL+/ZgGhKw1WbByCJAnoDMs6VGIe0piOby2Mg5AMrUuz2HchoXwB/50gJjPilik0M45kclrI6Ng0EceH6CPafiXZ9nQNQ6KZa6nJj/P2jR/HQi9N46yXrccdrx10f42QwpOA4WhRzsBbdBUc9kJeYT3DjUCvcOFj0BRX83s077H9v6PWXuZXYbmgwpCBPTf9wNwej3Qj7pYpupQlHz5rLN/VhOpbp2lkOTvxWavSklclVynxSw+6Nvfjsf72iptdj/vRWXFv2efBqOitzd3VzP7Ra4cahAut7/WXKgaWujg8EEc/wm8yNiE+u6FZii9/4QBA3XjiC33791maemqcZ6w9UVA7RlFZUv7Acw7ZbqfnKoT8ogxDToHmRhaSpHPjndnl4zKECG3r9ZTGHFyZiEAiwZ6zXDkh3u1uplIhfqpitxBa/8f4gAorommverYwPBO0akFIWU+7tuSvBDEmzu7ICZgv8voBc1GamUeTzFD85OgdKa0+bnUtqIAQYCPLP7XJw41CBdT1+zCdVaHphSteLEzHsXBdByCfZqawDIe67dBLxV25/fvhcAhG/hF6XVtPdznh/EDPxLHSjfCpcNK2taDHbOhyGIhWylppNs6qknzw+j7u+9Bz2n6mtaSZgKof+oNK1zR5XAr9CFWBTtljHUUopXpqI2ZklQxEFhAAb+vju10nEL7vGHL65fxLfen4Sv3TJhhaclfcZHwjAyNMytZrNGUhrxoqUw2271+OpP7kJfS3aHQ+GfU1xK52xsrtOL7grLjcWkitz0XUzPOZQgfVWr/xjs0k88PPTuOnCEcQzOds4jET8+M7/uB67NvS08jQ9R8RvNuJzcmYhhXu+9TJet30IH7m9JdNgPQ8bLDMRTRdNIIumVz5nWxAIRlroshsMKThmxecayWTUdFNOx2qbxQKYTfcGudqvCa4cKsBqHf7h8eP4xydP4e4H9wMA9m4qjOvbO94HpQXpgl6GTdVz+oG/f+Ac9DzFx96xh1+vCrDCy8mSoDRzz/S3kY98MKzYgd9GMmUZh6loBvNJFb/+hafLMr50I4+PPnzIzpTr9smNK4F/Uiuw3nIr7T8TRW9ARlLVEVJE7BgpH7DCKRDxyzDytKgN9Q9emcGlY73Y2MRe9O3Ghj4/BIKyoHQ0Zbro2illeiDkQyyTc42f1JNJSzFMxTJ49uQinjtlfjl5aTKGL/zkJB6yelfNJ1Ve41Aj3K1UgR6/ZA9i+cBN2xFNa9D0vOcGgHiNiNVfKZnVEVQkTMUyeGkyjj+57cIWn5m3kUUBG3oD9g6XsWi7ldoniD8UVkApEE3nGlrkaCuHWAbHZhMAyrsaPHPSNBYTi2lkcwYSWZ3HHGqEG4cKEEKwvtePU/Mp/MreUZ52WSPMOCxldYz0AD88eA6AGSTlVGd8IGAPt2dE29Ct5CyEa5RxyOYMzCdVSALBVCyDo+dN41Aaf3jWUhJnFtI4NW8GsC/w4CAfL8LdSlXYO96H23av54ZhBTDjwDKWnjo2j63DoZp6AnU7W4bCODlf3F9pMWXm5bdTc0fWL2smXnugeKVMWUbg4o290PQ8nrUUglM55Iw89p02j59dTNtFrNv4vVgTXDlU4e9/fe+KCmw4sDuzsoyls4tpbOcfxprYNhxCLJ2zO/8CZrZSb0Buq7z8zYNmcP3MClJMVwpzKV29ZQAvTsTsthhO5fDKVBxpzcCOkTCOzyVxaGYJhGDZzrYck/a541oEITzGsBIKysHMWJqKZexMHE512I6WjQMFrBbxbeRSAsyBQ0FFrNhlth4w5XDVlgH7WG9ALlIOzKX0zivHQCnw4yNzGO8P8n5eNcKNA6eusFbciWwO0XQOac3Axn6epVQLbEfLxoECpnJYSQGcFyCE4ILBUEOVw2Q0DVEguHxTv33sdduHEM/kkNZM1frsyQVsGw7Zj3l1ZgnbR7iKrRVuHDh1pTDwR7el/xg3DjUx1h+EIgo4Me9UDjn0t2G7kc2DwYYoh2hKw98+chj7z0SxvseP/pCCiE8CIcDrdgwBAKZjWRh5in2no7h66yA2DRaUKzcOtcNjDpy6UlAOul2QxOsbakMUCC4YDOLkXAo/PTqHaFpDNKXh4tH2q8LfPBTCjw6dh27k6xov+cpzZ/GZJ04AKLiUNvYHkNYMe4TqTDyDjGYgoeq4essAhsM+Oy19G4831Aw3Dpy6IgoEIUU0lUOs0IWVUxtbh0M4NJPAH/7bS4ilNVDaXgVwjM2DQeQMs1eUsx3IWvnRofN4zYYe/NIl63Hxxl4AwDuuGEOeUoxam5CZWBZHzpmprddsHQQhBJsGgjh8LsGVwwrgxoFTdyJ+GUk1h8koRdgnoSfAb7Na2TocxiMHzwMABALkKdou5gAUaglOL6TqZhzmEipenIjhD960Ex94Y2EwF5sLoul5EGIGqw9OL2HzYNBOQ79g0DIOw7zDQa3wmAOn7rC23ZPRDMb6AzzjawUw18jFG3tw9w3bALTn7IHNtnFYfVDayFO871/248FnzgAAnjg8C0qBm18z4vp4RRIwFPZhKpbBc6cWcPWWQft3120bwhUX9PN28SuAb+k4dSfsl7CY0hDP5Hi8YYVcMma6Sj5w03Zcu20Is0tZXLttcJlneY91PT74ZQFnrKK+p47NI5rW8LZLRzETz2AuoWLPWF/V1/jGvgk8/Mo5PHZoFtdsGcB/HpjBaK+/aifk0V4/vv38JPIUeP3OIfv4Xddtxl3Xba7Le+sWuHHg1J2rtwziCz89AVkQcLUjD52zPBet78Gzf3qz7Q75+3ftbfEZrQ5CCDYPhuyWFf/7P1/FdCyDW3evx5995xU8dWwe33n/ddg92uv6/KSq4+M/PII9Y704s5DG2/7hKWRzefz+zTuqKtHLNvVjJp7FH9yyk88OWSPcrcSpO79zw1aEFAmakec1DqugU9q1XDrWh2dPLeL0fApHzieQUHU8dug8njo2D83I4/e+8kLZ7A/GvzxzBvNJDR+5/WJ85PbdiPhl/M2v7cEfvGmH6+MZH/6V3Xjuz96EO6/axN2Za4QbB07d6Q8peK8VJOTV0d3LL+3ZgKSq4yP/+SoAM8D+V987BM3I449vvRCnFlL4h8ePI5HN4V1feBpffvq0/dx/f2EKV1zQj73jfbh970b84s/ehF9/7Thf8JsIdytxGsJ7b9gCSSR4w87hVp8Kp0Vcu20Q/UEZjx+exQYrVvDY4VkMR3x43xu24fC5BL789GksJFU8e2oRz55ahJGnuHbbIA6fS+Ajt+9u9VvoajyjHAghtxFCjhBCjhNC7mn1+XDWRlCR8P6btiPk4/uPbkUWBdx2sen3v/HCYdyyax0A4M271kEQCD548w5kcwb+bf8kfvXyjbh19zr85X+8ig99/SWIAsFbecygpXjCOBBCRACfAfAWALsA3EkI2dXas+JwOGvl9r2jAIA3vWYd3rx7PS4d68WdV20CYLay+LXLx9AXlHHvW16DT995Od68ax1enVnC67YP8YltLYZ4oSU1IeRaAB+mlN5q/fteAKCUfrTSc6688kq6b9++Jp0hh8NZLSfmktg6FHKNFxh5imRWt+sPND2PzzxxHG/eva5iJhNnbRBC9lNKr1zucV7R/BsBTDj+PQng6hadC4fDqSPVhuuIAikqTFMkAX9wy85mnBZnGTzhVgLgloJQJmkIIXcTQvYRQvbNzc014bQ4HA6nO/GKcZgEMO749xiA6dIHUUrvp5ReSSm9cniYZ8FwOBxOo/CKcfgFgB2EkC2EEAXAHQAeavE5cTgcTtfiiZgDpVQnhHwAwCMARABfopQebPFpcTgcTtfiCeMAAJTS7wP4fqvPg8PhcDjecStxOBwOx0Nw48DhcDicMrhx4HA4HE4ZnqiQXg2EkASAIzU+vBdAvI5/vt6vBwBDAObr8Drt8F4b8br1un6A96+hl68d4P336+Xr18j3ys7zAkrp8rUAlNK2/AKwbwWPvb/Of7uur7fS99Pu77VB51mX69cO19DL165N3q9nr18j3+tKz7Nb3Er/4fHXqyft8l75NfTO69Ubr79fL18/z7zXdnYr7aM1NI9qFzrt/TQbfv1WD792a6Ndrt9Kz7OdlcP9rT6BOtNp76fZ8Ou3evi1Wxvtcv1WdJ5tqxw4HA6H0zjaWTlwOBwOp0Fw49AgCCHjhJAnCCGHCCEHCSEftI4PEEIeJYQcs773W8cHrccnCSH/UPJadxJCDhBCXiaE/IAQMtSK99RM6nz93mVdu4OEkL9pxftpJqu4drcQQvZb99h+QsgbHa91hXX8OCHkU8RtYk+HUefrdx8hZIIQkmzV+1k19Uyb4l9FKWQbAFxu/RwBcBTmCNS/AXCPdfweAB+zfg4BeB2A3wXwD47XkQDM/v/t3V+IVGUYx/HvQyuFmVmmkVQs3VQmoRSk/aGL6MJuCuyiiHazm6wgugsjqJsukpLILrZIQyvCwiIrKkpIsrIuUvyTUG1IrSxFZLquFEW/Lt53aNgz6zq7Z/bMjL8PHM7sO+e8vO/D7HnOOXPmfYHz8t9rSLPmVd7HDonfXOAnYF7+eyNwU9X9a7PYLQEW5NeLgEN1dX0NLCPNufIBsLzq/nVY/Jbm+o5V3a9mF185tIikYUnf5NcjwAHSjHe3kg5Q5PVteZtRSTuAP8dUFXk5M5+1zabBXBfdpsT4XQJ8J6k2O9QnwIoWN79Sk4jdLkm1z9R+4IyIOD0iLgBmS/pS6Ui3qbZPNysrfvm9nZKGp7P9ZXFymAYR0Us6u/gKOL/2Ycnr+SfaV9LfwP3AXlJSWAisb2Fz285U4gf8AFwWEb0R0UP6h75ogn26xiRitwLYJekv0gFxqO69oVx2yphi/Dqak0OLRcQsYAvwsKSjk9h/Bik5LAEWAHuA1aU2so1NNX6SDpPitxn4DDgI/FNmG9tVs7GLiCuAp4D7akUNNjtlHm8sIX4dzcmhhfKBfQvwmqS3cvEv+XKdvP51gmoWA0gazJf2bwDXtqjJbaWk+CHpXUnXSFpGGo/r+1a1uV00G7uIuBB4G+iTNJiLh0hT9tY0nL63G5UUv47m5NAi+fuB9cABSWvr3toK9OfX/cA7E1R1CFgYEbWBsm4m3QPtaiXGj4iYn9fnAA8AL5Xb2vbSbOwiYg7wPrBa0ue1jfOtk5GIWJrr7OMk4t3pyopfx6v6G/FuXUhPzoh0G2h3Xm4hPT2zjXT2ug04t26fg8DvwDHSWdvCXL6KlBD2kMZKmVt1/zosfq8D3+bljqr71m6xAx4DRuu23Q3Mz+9dDewDBoHnyT+c7eal5PityZ/Ff/P6iar7d7KLfyFtZmYFvq1kZmYFTg5mZlbg5GBmZgVODmZmVuDkYGZmBU4OZi0QEasioq+J7XsjYl8r22TWjJ6qG2DWbSKiR9JA1e0wmwonB7MG8oBrH5IGXFtCGra5D7gcWAvMAn4D7pE0HBGfAl8A1wFbI+Is0jDNT0fEYmAAmEn6Mdm9kg5HxFXABuA4sGP6emc2Md9WMhvfpcCLkq4EjgIPAuuA2yXVDuxP1m0/R9KNkp4ZU88m4JFcz17g8Vz+MvCQ0phPZm3FVw5m4/tZ/4+V8yrwKGkyl4/zhGinAfVj9W8eW0FEnE1KGttz0UbgzQblrwDLy++C2eQ4OZiNb+zYMiPA/hOc6Y82UXc0qN+sbfi2ktn4Lo6IWiK4E9gJzKuVRcSMPIb/uCQdAQ5HxA256G5gu6Q/gCMRcX0uv6v85ptNnq8czMZ3AOiPiBdII3GuAz4Cnsu3hXqAZ0lTQ55IPzAQETOBH4GVuXwlsCEijud6zdqGR2U1ayA/rfSepEUVN8WsEr6tZGZmBb5yMDOzAl85mJlZgZODmZkVODmYmVmBk4OZmRU4OZiZWYGTg5mZFfwHRNvqBP6YBHsAAAAASUVORK5CYII=\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], "source": [ "sorted_raw_data['inc'][-200:].plot()" ] @@ -2287,7 +2333,7 @@ "Etant donné que le pic de l'épidémie se situe en hiver, à cheval entre deux années civiles, nous définissons la période de référence entre deux minima de l'incidence, du 1er septembre de l'année $N$ au 1er septembre de l'année $N+1$.\n", "Notre tâche est un peu compliquée par le fait que l'année ne comporte pas un nombre entier de semaines. Nous modifions donc un peu nos périodes de référence: à la place du 1er septembre de chaque année, nous utilisons le premier jour de la semaine qui contient le 1er septembre.\n", "Comme l'incidence de syndrome de varicelle est très faible en été, cette modification ne risque pas de fausser nos conclusions.\n", - "Encore un petit détail: les données commencent en decembre 1990, ce qui rend la première année incomplète. Nous commençons donc l'analyse en 1992." + "Encore un petit détail: les données commencent en decembre 1990, ce qui rend la première année incomplète. Nous commençons donc l'analyse en 1992.L'année 2020 est incomplete nous arrêtons l'analye à 2019" ] }, { @@ -3283,18 +3329,18 @@ }, { "cell_type": "code", - "execution_count": 26, + "execution_count": 87, "metadata": {}, "outputs": [], "source": [ "first_septembre_week =[pd.Period(pd.Timestamp(y, 9, 1), 'W')\n", " for y in range(1991,\n", - " sorted_raw_data.index[-1].year)]" + " 2020)]" ] }, { "cell_type": "code", - "execution_count": 30, + "execution_count": 73, "metadata": {}, "outputs": [ { @@ -3332,7 +3378,7 @@ " Period('2020-08-31/2020-09-06', 'W-SUN')]" ] }, - "execution_count": 30, + "execution_count": 73, "metadata": {}, "output_type": "execute_result" } @@ -3341,6 +3387,205 @@ "first_septembre_week" ] }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "En partant de cette liste des semaines qui contiennent un 1er septembre, nous obtenons nos intervalles d'environ un an comme les périodes entre deux semaines adjacentes dans cette liste. Nous calculons les sommes des incidences hebdomadaires pour toutes ces périodes.\n", + "Nous vérifions également que ces périodes contiennent entre 51 et 52 semaines, pour nous protéger contre des éventuelles erreurs dans notre code." + ] + }, + { + "cell_type": "code", + "execution_count": 88, + "metadata": {}, + "outputs": [], + "source": [ + "year = []\n", + "yearly_incidence = []\n", + "for week1, week2 in zip(first_septembre_week[:-1],\n", + " first_septembre_week[1:]):\n", + " one_year = sorted_raw_data['inc'][week1:week2-1]\n", + " assert abs(len(one_year)-52) < 2\n", + " yearly_incidence.append(one_year.sum())\n", + " year.append(week2.year)\n", + " \n", + "yearly_incidence = pd.Series(data=yearly_incidence, index=year) \n", + "\n", + " \n", + "\n", + " \n", + " \n", + " " + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "Visualisation des incidences annuelles puis de la repartition" + ] + }, + { + "cell_type": "code", + "execution_count": 89, + "metadata": {}, + "outputs": [ + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAaEAAAEWCAYAAADPZygPAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzt3X2cVdV97/HPF4eAGjGggaiomGo0mlTUCZjmNh1qA5qkoo2mRKM0scEYNTGvpkYbe7ECubVNasO1esVKfIpPIbHaRqIEOwltCQhKoqgEEgiiRDSDwphAQH73j70ObMZ5ODPOYZ+H7/v1Oq/Z53f22rPWbJjfrLXXXlsRgZmZWREGFF0BMzNrXE5CZmZWGCchMzMrjJOQmZkVxknIzMwK4yRkZmaFcRIySyT9oaQV3Xx+q6Tpe7JOtUhSq6S/TNt/Iem/iq6TVS8nIasLktZI+pM3c4yIWBARR/dXncysZ05CZmZWGCchqzulISBJX5O0UdJqSaflPh8m6ZuSXkif/1uKt0hal9vvBEmPS9os6V5gcIfv81FJyyS9Iul/JP1+7rM1kr4k6aeSXpV0r6TBuc8nprKbJP1c0qkpvr+kWyStl/S8pOmS9uqinWMkLUzff72k6yW9Jfd5SPqspJWpnf8iSWX+jHbrWUq6WtKdufcnpza/IuknklrKPDfHSJonqU3SCkkfL6ec1S8nIatXY4EVwIHAPwC3lH4BA3cA+wDHAcOB6zoWTr/M/y3tOwz4NvCx3OcnArOBC4EDgJuAByUNyh3m48CpwBHA7wN/kcqOAW4H/hp4G/BBYE0qcxuwHTgSOAEYD/xlF218HfhiauP7gVOAz3XY56PA+4DjU30mlPkz6pKkQ4DvAdPJfjZfAr4j6e09lNsXmAfcRfZz/wRwg6TjevqeVr+chKxe/TIibo6I18l+sR8EjJB0EHAa8NmI2BgR2yLih52UPxkYCPxz2mcO8Fju888AN0XEooh4PSJuA7amciUzI+KFiGgD/h0YneIXALMjYl5E7IiI5yPiWUkjUt0ui4jXImIDWYKc1FkDI2JpRPw4IrZHxBqyRPhHHXb7+4h4JSLWAv+Zq0OXP6OufqA5nwQeioiHUv3nAUuAD/dQ7qPAmoj4Zqrz48B3gLPK+J5Wp5qKroBZhfyqtBERv0l/4L+V7C/3tojY2EP5g4HnY/cVfn+Z2z4cmCzp0lzsLancG+oA/Cb32aHAQ518z8PJEt/6XIdkAPBcZxWU9C7gn4Bmsp5dE7C0w24d6/DWzj7r8DPqyeHA2ZL+NBcbSJbkeio3VtIruVgTWW/TGpSTkDWa54Bhkt4WEa90s9964BBJyiWiw4Cf544zIyJm9LEOv9dFfCtwYERsL+M4NwJPAJ+IiM2SLqP/ehWvkSW2knd0qOcdEfGZXh7zOeCHEfGhN1s5qx8ejrOGEhHrgblk1yKGShoo6YOd7LqQ7NrM5yU1SfozYEzu85uBz0oaq8y+kj4iab8yqnEL8ClJp0gaIOkQScekuj0CfF3SkPTZ70nqOMRWsh+wCWiXdAxwUXk/hbIsAyaln08zuye3O4E/lTRB0l6SBqdJHSN7OOZ/AO+SdF467kBJ75P07n6st9UYJyFrROcB24BngQ3AZR13iIjfAX9GNplgI/DnwHdzny8huy50ffp8Vdq3RxGxGPgU2fWeV4Efkg1VAZxPNqz3dDruHLJrNZ35EnAOsJksKd5bzvcv09+S9dY2An9HNpmgVP/ngInA3wAvkfVw/poefp9ExGayiRaTgBfIhgOvBQZ1V87qm/xQOzMzK4p7QmZmVhgnITMzK4yTkJmZFcZJyMzMCuP7hJIDDzwwRo0aVXQ13pTXXnuNfffdt+hqVFS9t7He2wf138ZGa9/SpUtfjohul2zqjpNQMmrUKJYsWVJ0Nd6U1tZWWlpaiq5GRdV7G+u9fVD/bWy09kn6Zdd796yiw3GSvihpuaSnJN2dbmq7Oq0OvCy9Ppzb/0pJq9LquhNy8ZMkPZk+m5lbCXhQWp14laRFkkblykxOqwevlDS5ku00M7O+qVgSSivtfh5ojoj3AHuxayHG6yJidHo9lPY/Nn1+HNnKwzfklrC/EZgCHJVep6b4BcDGiDiS7Ma/a9OxhgFTyVYJHgNMlTS0Um01M7O+qfTEhCZgb0lNZOtQvdDNvhOBeyJia0SsJrsDfUxa9XhIRCxMa3jdDpyRK3Nb2p4DnJJ6SROAeRFRWqhyHrsSl5mZVYmKJaGIeB74GrCWbDHIVyPikfTxJcoe9jU710M5hN1XC16XYoek7Y7x3cqkBR9fJXu2S1fHMjOzKlKxiQkpuUwke6DXK8C3JX2SbGhtGhDp69eBTwOdPUwruonTxzL5Ok4hG+ZjxIgRtLa2dt2gGtDe3l7zbehJvbex3tsH9d9Gt693Kjk77k+A1RHxEoCk7wJ/EBH5RwTfTLayLmS9lUNz5UeSDd+tS9sd4/ky69KQ3/5AW4q3dCjT2rGCETELmAXQ3NwctT6jpd5n5UD9t7HW2rdh0xYuufsJrj/nBIbvN7jnAtReG3vL7eudSl4TWgucLGmfdJ3mFOCZdI2n5EzgqbT9INnS8YMkHUE2AWFxWt5+c3qmvchWGX4gV6Y08+0s4NF03ehhYHxaqn8o2cq9D1euqWaNaeb8lTy2po2ZP1hZdFWsRlWsJxQRiyTNAR4ney7LE2S9jn+VNJpseGwNcGHaf7mk+8iWsN8OXJweOwzZc1JuBfYmexbM3BS/BbhD0iqyHtCkdKw2SdPY9Tjma9Ijls2sHxx91Vy2bt+x8/2di9Zy56K1DGoawIrppxVYM6s1Fb1ZNSKmkk2Vzjuvm/1nAG94UmV6dst7OolvAc7u4lizgdm9qa+ZlWfB5eOY/tAzPLL8V2zZtoPBAwcw4bh38JWP+Pl01jteO87Mem34kMHsN6iJrdt3MKhpAFu372C/QU1lXxcyK/GyPWbWJy+3b+XcsYdzzpjDuGvxWl7avKXoKlkNchLqB32ZIWRW6246r3nn9vQz3jBablYWD8f1A88QMjPrG/eE3gTPEDLrnQ2btvDVRb/l2JO2eNTAAPeE3pQFl4/j9NEHM3hg9mMcPHAAE0cfzIIvjyu4ZmbVaeb8lazcuMOjBraTe0JvgmcImZXHowbWFfeE3qTSDKH7P/cBzh17OC+1by26SmZVx6MG1hX3hN4kzxAy61l+1GDgADxqYDs5CZnZHlEaNXjXgBf52Y4Rvq/IACchM9tDSqMGra0vc36LRw0s42tCZmZWGCchMzMrjJOQmZkVxknIzMwK4yRkZmaFcRIyM7PCOAmZmVlhnITMzKwwTkJmZlYYJyEzMyuMk5CZmRWmoklI0hclLZf0lKS7JQ2WNEzSPEkr09ehuf2vlLRK0gpJE3LxkyQ9mT6bKUkpPkjSvSm+SNKoXJnJ6XuslDS5ku00M7O+qVgSknQI8HmgOSLeA+wFTAKuAOZHxFHA/PQeScemz48DTgVukLRXOtyNwBTgqPQ6NcUvADZGxJHAdcC16VjDgKnAWGAMMDWf7MzMrDpUejiuCdhbUhOwD/ACMBG4LX1+G3BG2p4I3BMRWyNiNbAKGCPpIGBIRCyMiABu71CmdKw5wCmplzQBmBcRbRGxEZjHrsRlZmZVomKPcoiI5yV9DVgL/BZ4JCIekTQiItanfdZLGp6KHAL8OHeIdSm2LW13jJfKPJeOtV3Sq8AB+XgnZXaSNIWsh8WIESNobW3te4OrQHt7e823oSf13sZ6bx/Ufxvdvt6pWBJKw18TgSOAV4BvS/pkd0U6iUU38b6W2RWImAXMAmhubo6WlpZuqlf9WltbqfU29KTe21jv7YP6b6Pb1zuVHI77E2B1RLwUEduA7wJ/ALyYhthIXzek/dcBh+bKjyQbvluXtjvGdyuThvz2B9q6OZaZmVWRSiahtcDJkvZJ12lOAZ4BHgRKs9UmAw+k7QeBSWnG2xFkExAWp6G7zZJOTsc5v0OZ0rHOAh5N140eBsZLGpp6ZONTzMzMqkglrwktkjQHeBzYDjxBNvT1VuA+SReQJaqz0/7LJd0HPJ32vzgiXk+Huwi4FdgbmJteALcAd0haRdYDmpSO1SZpGvBY2u+aiGirVFvNzKxvKpaEACJiKtlU6bytZL2izvafAczoJL4EeMND6SNiCymJdfLZbGB2L6tsZmZ7kFdMMDOzwjgJmZlZYZyEzMwKsGHTFj5+00I2bN5SdFUK5SRkZlaAmfNX8tiaNmb+YGXRVSlURScmmJnZ7o6+ai5bt+/Y+f7ORWu5c9FaBjUNYMX00wqsWTHcEzKrIA+5WEcLLh/H6aMPZvDA7Nfv4IEDmDj6YBZ8eVzBNSuGk5BZBXnIxToaPmQw+w1qYuv2HQxqGsDW7TvYb1ATw/cbXHTVCuHhOLMK8JCLdefl9q2cO/ZwzhlzGHctXstLDdxTdhIyq4AFl49j+kPP8MjyX7Fl2w4GDxzAhOPewVc+8u6iq2ZV4KbzmnduTz/jDffhNxQPx5lVgIdczMrjnpBZhXjIxaxnTkJmFeIhF7OeeTjOzMwK4yRkZmaFcRIyM7PCOAmZmVlhnITMeslL8Zj1HychK1Qt/kL3Ujxm/cdTtK1Q+V/o0898b9HV6ZaX4jHrf05CVoha/IXupXjM+p+H46wQtbicvZfiMet/FUtCko6WtCz32iTpMklXS3o+F/9wrsyVklZJWiFpQi5+kqQn02czJSnFB0m6N8UXSRqVKzNZ0sr0mlypdlrf1Oov9NJSPPd/7gOcO/ZwXmrfWnSVzGpaxYbjImIFMBpA0l7A88D9wKeA6yLia/n9JR0LTAKOAw4GfiDpXRHxOnAjMAX4MfAQcCowF7gA2BgRR0qaBFwL/LmkYcBUoBkIYKmkByNiY6Xaa71Xi2urNcJSPBs2beGSu5/g+nNOqPo/Cqz27alrQqcAP4+IX6ZOTGcmAvdExFZgtaRVwBhJa4AhEbEQQNLtwBlkSWgicHUqPwe4PvWSJgDzIqItlZlHlrjurkDbrI8a4Rd6LaqlySJW+/ZUEprE7gngEknnA0uAv0o9lEPIejol61JsW9ruGCd9fQ4gIrZLehU4IB/vpMxOkqaQ9bAYMWIEra2tfWxedWhvb6/5NvSk3ttYZPs+88hrbNs1V2TnZJGBA+Dm8fv22/fxOaxt/d2+iichSW8BTgeuTKEbgWlkw2TTgK8DnwY66yJFN3H6WGZXIGIWMAugubk5WlpaumpGTWhtbaXW29CTem9jke377xO3dDn7rz+H5XwOa1t/t29PzI47DXg8Il4EiIgXI+L1iNgB3AyMSfutAw7NlRsJvJDiIzuJ71ZGUhOwP9DWzbHMrAu1OlnEatueSEKfIDcUJ+mg3GdnAk+l7QeBSWnG2xHAUcDiiFgPbJZ0crrecz7wQK5MaebbWcCjERHAw8B4SUMlDQXGp5iZdcOz/2xPq+hwnKR9gA8BF+bC/yBpNNnw2JrSZxGxXNJ9wNPAduDiNDMO4CLgVmBvsgkJc1P8FuCONImhjezaExHRJmka8Fja75rSJAUz65oni9ieVtEkFBG/IZsokI+d183+M4AZncSXAG/4HxERW4CzuzjWbGB2L6tsZmZ7kFdMMDOzwjgJmZlZYZyEzMysME5CZmZWGCchMzMrjJOQmZkVxknIzMwK4yRkZmaFcRIyM7PCOAmZmVlhnITMzKwwTkJmDWDDpi18/KaFbKiBR6hbY3ESMmsA+Ud2m1WTPfV4bzMrwNFXzWXr9l3P7C49sntQ0wBWTD+twJqZZdwTMqtjCy4fx+mjD2bwwOy/+uCBA5g4+mAWfHlcwTUzyzgJmdUxP7Lbql1Zw3HpsdrnAu+MiGskHQa8IyIWV7R2ZvamlR7Zfc6Yw7hr8Vpe8uQEqyLlXhO6AdgB/DFwDbAZ+A7wvgrVy8z6iR/ZbdWs3CQ0NiJOlPQEQERslPSWCtbLzMwaQLnXhLZJ2gsIAElvJ+sZmZmZ9Vm5SWgmcD8wXNIM4L+Ar1asVmZm1hDKSkIR8S3gcuD/AOuBMyLi292VkXS0pGW51yZJl0kaJmmepJXp69BcmSslrZK0QtKEXPwkSU+mz2amiRJIGiTp3hRfJGlUrszk9D1WSprcmx+KmZntGWUlIUknA89HxL9ExPXAOkljuysTESsiYnREjAZOAn5D1pu6ApgfEUcB89N7JB0LTAKOA04FbkhDgAA3AlOAo9Lr1BS/ANgYEUcC1wHXpmMNA6YCY4ExwNR8sjMzs+pQ7nDcjUB77v1rKVauU4CfR8QvgYnAbSl+G3BG2p4I3BMRWyNiNbAKGCPpIGBIRCyMiABu71CmdKw5wCmplzQBmBcRbRGxEZjHrsRlZg3M6+hVl3KTkFICACAidtC7JX8mAXen7RERsT4dZz0wPMUPAZ7LlVmXYoek7Y7x3cpExHbgVeCAbo5lZg3O6+hVl3ITyS8kfZ5dvZ/PAb8op2Cayn06cGVPu3YSi27ifS2Tr9sUsmE+RowYQWtraw9VrG7t7e0134ae1Hsb6719UFwbP/PIa2zLzektraM3cADcPH7ffvs+9X4O+7t95Sahz5LNkLuK7Jf5fNIv7zKcBjweES+m9y9KOigi1qehtg0pvg44NFduJPBCio/sJJ4vs05SE7A/0JbiLR3KtHasWETMAmYBNDc3R0tLS8ddakprayu13oae1Hsb6719UFwb//vELUx/6BkeWf4rtmzbweCBA5hw3Dv4ykfe3a/LGNX7Oezv9pU7O25DREyKiOERMSIizomIDT2XBOAT7BqKA3gQKM1Wmww8kItPSjPejiCbgLA4DdltlnRyut5zfocypWOdBTyahg0fBsZLGpomJIxPMTNrUF5HrzqVu3bc24HPAKPyZSLi0z2U2wf4EHBhLvz3wH2SLgDWAmenYy2XdB/wNLAduDgiXk9lLgJuBfYG5qYXwC3AHZJWkfWAJqVjtUmaBjyW9rsmItrKaauZ1S+vo1d9yh2OewBYAPwAeL2HfXeKiN+QTRTIx35NNluus/1nADM6iS8B3rDoVURsISWxTj6bDcwut65mVv+8jl71KTcJ7RMRX65oTczMqtCGTVu45O4nuP6cEzx0VwHlTtH+D0kfrmhNzMyqkKd0V1a5PaEvAH8j6XfA78imQEdEDKlYzaxq+C9Ba0R+NPqeUe7suP0iYkBEDI6IIem9E1CD8F+C1oj8aPQ9o7dPVj0iIqZJOhQ4yE9WrW/+S9Aamad07xnlXhO6AXg/cE563w78S0VqZFWjL38Jel0uqyelKd33f+4DnDv2cF5q31p0leqOn6xqXerLX4L5obvpZ753D9bWrP95SnfllZuE/GTVBlXuzX0eujOzvig3CXV8supZZOvIWZ0r9y/BBZeP63JdLjOzrpSVhCLiW5KWkq10ILInqz5T0ZpZTfFFXDPri26TUHpCackGcguRShrm9dgsz+tymVlv9dQTWsqu5/McBmxM228jW3z0iIrWzmqKL+KaWW91O0U7Io6IiHeSPQbhTyPiwIg4APgo8N09UUEzM6tf5d4n9L6IeKj0JiLmAn9UmSqZmVmjKHd23MuSrgLuJBue+yTw64rVyszMGkK5PaFPAG8nm6b9b8DwFLMK8uoDZlbvyl3AtC0ivhARJ6TXFzwzrvK8cKiZ1buepmj/c0RcJunfSasl5EXE6RWrWQPz6gNm1ih6uiZ0R/r6tUpXxHbx6gNm1ii6TUIRsTRtLgF+GxE7ANI6coMqXLeG5dUHzKxRlDsxYT6wT+793sAP+r86VtLbJeQ3bNrCVxf91pMYzKymlDtFe3BEtJfeRES7pH26K2BvTm9XH5g5fyUrN+7wIxTMrKaU2xN6TdKJpTeSTgJ+21MhSW+TNEfSs5KekfR+SVdLel7SsvT6cG7/KyWtkrRC0oT895P0ZPpsZnrSK5IGSbo3xRdJGpUrM1nSyvSaXGY7a87RV81l1BXf485FawmySQyjrvgeR181t+iqmZn1qNwkdBnwbUkLJC0A7gUuKaPcN4DvR8QxwPFAaeXt6yJidHo9BCDpWGAScBxwKnBDuvYEcCMwBTgqvU5N8QuAjRFxJHAdcG061jBgKjAWGANMlTS0zLbWlL48/dTMrFqU+yiHxyQdAxxNtoDpsxGxrbsykoYAHwT+Ih3jd8DvUiemMxOBeyJiK7Ba0ipgjKQ1wJCIWJiOeztwBjA3lbk6lZ8DXJ96SROAeaV7mSTNI0tcd1Nn8pMYBg7AkxjMrKaUe00I4H3AqFTmBElExO3d7P9O4CXgm5KOJ1uR+wvps0sknU826+6vImIjcAjw41z5dSm2LW13jJO+PgcQEdslvQockI93UmYnSVPIeliMGDGC1tbWbppTvZ5ZvYVxI5sYc8A2Fv96IE+vfp7W1vpcVam9vb1mz1M56r19UP9tdPt6p6wkJOkO4PeAZcDrKRxAd0moCTgRuDQiFkn6BnAFcD0wLZWfBnwd+DRZD6uj6CZOH8vsCkTMAmYBNDc3R0tLSzfN6R8bNm3hkruf4PpzTui33kqp2q2trXx2Uku/HLNatba2sifOU1HqvX1Q/210+3qn3GtCzcAHIuJzEXFpen2+hzLrgHURsSi9nwOcGBEvRsTr6Z6jm8mu2ZT2PzRXfiTwQoqP7CS+WxlJTcD+QFs3xyqcl+IxK4/XTmwM5Sahp4B39ObAEfEr4DlJR6fQKcDTkg7K7XZmOjbAg8CkNOPtCLIJCIsjYj2wWdLJ6XrP+cADuTKlmW9nAY9GRJA9/2i8pKFpQsL4FCvMbrPYwrPYzHriP9gaQ7nXhA4kSyCLgZ13TZaxdtylwLckvQX4BfApYKak0WTDY2uAC9Oxlku6D3ga2A5cHBGlob+LgFvJbpKdm14AtwB3pEkMbWSz64iINknTgMfSftcUveCql+IxK4/XTmws5Sahq/ty8IhYRjaUl3deN/vPAGZ0El8CvOGOzYjYApzdxbFmA7N7U99K8lI8ZuXxH2yNpdwp2j+sdEUaQWkpnnPGHMZdi9fykse6zd7Af7A1lp4e5bCZTmaVkc0+i4gYUpFa1aneLsVj1pVKzLKsJv6DrXH0tIr2fnuqImZWvvxF+3pcK9B/sDWO3tysamYF80V7qzflTtE2swor53EcXivQ6o2TkFmVyD+Ooyu+aG/1xsNxZgXr7RCbL9q/Ub1P1KhnTkJmBevtfTG+aP9G9T5Ro545CZkVzI/j6DtP1Kh9TkJmVaA0xPauAS/ysx0jPMRWJq+uUPuchKwu1do1gtIQW2vry5zf4iG2cnmiRu3z7DirGeVMYS7xCsyNo9SLvP9zH+DcsYfzUvvWngtZ1XBPyGpGfgpzVxeffY2g8XiiRm1zErKq15vE4msEZrXFw3FW9XqzSoCvEZjVFveErOr1dgqzb+Y0qx1OQlYTejOF2dcIzGqHk5DVBE9hNqtPviZkZmaFcRIyM+snvbmXzTJOQmZm/aScx3HY7iqahCS9TdIcSc9KekbS+yUNkzRP0sr0dWhu/yslrZK0QtKEXPwkSU+mz2ZKUooPknRvii+SNCpXZnL6HislTa5kO82ssR191VxGXfE97ly0liC7l23UFd/j6KvmFl21qlfpntA3gO9HxDHA8cAzwBXA/Ig4Cpif3iPpWGAScBxwKnCDpL3ScW4EpgBHpdepKX4BsDEijgSuA65NxxoGTAXGAmOAqflkZ2bWn/zE276rWBKSNAT4IHALQET8LiJeASYCt6XdbgPOSNsTgXsiYmtErAZWAWMkHQQMiYiFERHA7R3KlI41Bzgl9ZImAPMioi0iNgLz2JW4zMz6lR/H0XeVnKL9TuAl4JuSjgeWAl8ARkTEeoCIWC9peNr/EODHufLrUmxb2u4YL5V5Lh1ru6RXgQPy8U7K7CRpClkPixEjRtDa2trXtlaF9vb2mm9DT+q9jfXePqjfNj6zegvjRjYx5oBtLP71QJ5e/Tytrb8uulr9rr/PXyWTUBNwInBpRCyS9A3S0FsX1Eksuon3tcyuQMQsYBZAc3NztLS0dFO96tfa2kqtt6En9d7Gem8f1G8bS01qbW3ls5NaiqxKRfX3+avkNaF1wLqIWJTezyFLSi+mITbS1w25/Q/NlR8JvJDiIzuJ71ZGUhOwP9DWzbHMzKyKVCwJRcSvgOckHZ1CpwBPAw8Cpdlqk4EH0vaDwKQ04+0IsgkIi9PQ3WZJJ6frPed3KFM61lnAo+m60cPAeElD04SE8SlmZmZVpNLL9lwKfEvSW4BfAJ8iS3z3SboAWAucDRARyyXdR5aotgMXR8Tr6TgXAbcCewNz0wuySQ93SFpF1gOalI7VJmka8Fja75qIaKtkQ83MrPcqmoQiYhnQ3MlHp3Sx/wxgRifxJcAbFgyLiC2kJNbJZ7OB2b2pr5mZ7VleMcHMzArjJGRm1sA2bNrCx29aWNh6d05CZmYNbOb8lTy2pq2w9e78PCEzswZ09FVz2bp9x873dy5ay52L1jKoaQArpp+2x+rhnpCZWQOqlvXunITMzBpQfr27QU0DClvvzsNxZmYN6uX2rZw79nDOGXMYdy1ey0sFTE5wEjIza1A3nbfrNs7pZ7zhVsw9wsNxZmZWGCchMzMrjJOQmZkVxknIzMwK4yRkZmaFcRJqUEWvF2VmBk5CDavo9aLMzMD3CTWcalkvyswM3BNqONWyXpSZGTgJNZxqWS/KzAw8HNeQqmG9KDMzcBJqSNWwXpSZGXg4zszMCuQkZGZW5er5vr6KJiFJayQ9KWmZpCUpdrWk51NsmaQP5/a/UtIqSSskTcjFT0rHWSVppiSl+CBJ96b4IkmjcmUmS1qZXpMr2U4zs0qq5/v69sQ1oXER8XKH2HUR8bV8QNKxwCTgOOBg4AeS3hURrwM3AlOAHwMPAacCc4ELgI0RcaSkScC1wJ9LGgZMBZqBAJZKejAiNlaslWZm/awR7uurpuG4icA9EbE1IlYDq4Axkg4ChkTEwogI4HbgjFyZ29L2HOCU1EuaAMyLiLaUeOaRJS4zs5rRCPf1VbonFMAjkgK4KSJmpfglks4HlgB/lRLFIWQ9nZJ1KbYtbXeMk74+BxAR2yW9ChyQj3dSZidJU8h6WIzPuJM4AAAKTUlEQVQYMYLW1ta+t7QKtLe313wbelLvbaz39kH9t7G/27f511vZum0HAwfA1m072PTrDTy99FWe7rfv0Dv93b5KJ6EPRMQLkoYD8yQ9Sza0No0sQU0Dvg58GlAn5aObOH0ssyuQJcVZAM3NzdHS0tJtY6pda2srtd6GntR7G+u9fVD/bezv9t393BLOPXnwbvf1tbQ0d1tmw6YtXHL3E1x/zgn9fiN6f7evokkoIl5IXzdIuh8YExE/Kn0u6WbgP9LbdcChueIjgRdSfGQn8XyZdZKagP2BthRv6VCmtV8aZWa2B/Xlvr78RIbpZ763UlXrFxVLQpL2BQZExOa0PR64RtJBEbE+7XYm8FTafhC4S9I/kU1MOApYHBGvS9os6WRgEXA+8H9zZSYDC4GzgEcjIiQ9DHxV0tC033jgykq11cysGtTiRIZK9oRGAPen2dRNwF0R8X1Jd0gaTTY8tga4ECAilku6D3ga2A5cnGbGAVwE3ArsTTYrbm6K3wLcIWkVWQ9oUjpWm6RpwGNpv2sioq2CbTUzK9yCy8cx/aFneGT5r9iybQeDBw5gwnHv4CsfeXfRVetSxZJQRPwCOL6T+HndlJkBzOgkvgR4Qz80IrYAZ3dxrNnA7F5U2cysptXiAsVeO87MrI7U2gLFTkJmZnWk1hYorqabVc3MrME4CZmZWWGchMzMrDBOQmZmVhgnITMzK4yTkJmZFcZJyMzMCuMkZGZmhXESMjOzwjgJWcPbsGkLH79pIRuqfHkTs3rkJGQNL//sFTPbs7x2nDWsWnz2ilm9cU/IGtaCy8dx+uiDGTww+28weOAAJo4+mAVfHldwzcwah5OQNaxafPaKWb3xcJw1tFp79opZvXESsoZWa89eMas3Ho4zM7PCOAmZmVlhnITMzKwwFU1CktZIelLSMklLUmyYpHmSVqavQ3P7XylplaQVkibk4iel46ySNFOSUnyQpHtTfJGkUbkyk9P3WClpciXbaWZmfbMnekLjImJ0RJSuAF8BzI+Io4D56T2SjgUmAccBpwI3SNorlbkRmAIclV6npvgFwMaIOBK4Drg2HWsYMBUYC4wBpuaTnZmZVYcihuMmArel7duAM3LxeyJia0SsBlYBYyQdBAyJiIUREcDtHcqUjjUHOCX1kiYA8yKiLSI2AvPYlbjMzKxKVHqKdgCPSArgpoiYBYyIiPUAEbFe0vC07yHAj3Nl16XYtrTdMV4q81w61nZJrwIH5OOdlNlJ0hSyHhZAu6QVfW1olTgQeLnoSlRYvbex3tsH9d/GRmvf4W/mYJVOQh+IiBdSopkn6dlu9lUnsegm3tcyuwJZUpzVTZ1qiqQluWHPulTvbaz39kH9t9Ht652KDsdFxAvp6wbgfrLrMy+mITbS1w1p93XAobniI4EXUnxkJ/HdykhqAvYH2ro5lpmZVZGKJSFJ+0rar7QNjAeeAh4ESrPVJgMPpO0HgUlpxtsRZBMQFqehu82STk7Xe87vUKZ0rLOAR9N1o4eB8ZKGpgkJ41PMzMyqSCWH40YA96fZ1E3AXRHxfUmPAfdJugBYC5wNEBHLJd0HPA1sBy6OiNfTsS4CbgX2BuamF8AtwB2SVpH1gCalY7VJmgY8lva7JiLaKtjWalE3Q4vdqPc21nv7oP7b6Pb1grKOg5mZ2Z7nFRPMzKwwTkJmZlYYJ6EqJ2m2pA2SnsrFjpe0MC1l9O+ShqT4WyR9M8V/IqklV6Y1LYe0LL2Gd/Lt9jhJh0r6T0nPSFou6Qsp3m/LOxWpn9tXF+dQ0gFp/3ZJ13c4Vs2fwx7aV3XnsA/t+5Ckpek8LZX0x7lj9f78RYRfVfwCPgicCDyViz0G/FHa/jQwLW1fDHwzbQ8HlgID0vtWoLno9nTSvoOAE9P2fsDPgGOBfwCuSPErgGvT9rHAT4BBwBHAz4G90meLgfeT3Sc2FzitztpXL+dwX+B/AZ8Fru9wrHo4h921r+rOYR/adwJwcNp+D/D8mzl/7glVuYj4EdnMv7yjgR+l7XnAx9L2sWTr8RHZvVmvAFV901xErI+Ix9P2ZuAZstUt+nN5p8L0V/v2bK17p7dtjIjXIuK/gN0eY1sv57Cr9lWrPrTviUj3gALLgcHKbq3p0/lzEqpNTwGnp+2z2XVj7k+AiZKalN1rdRK737T7zTQE8LfVMMzRkbJV0E8AFtFheSeynh10vSTTIXS9vFNVeJPtK6mHc9iVejmHPanac9iH9n0MeCIittLH8+ckVJs+DVwsaSlZ9/l3KT6b7MQvAf4Z+B+ye64Azo2I9wJ/mF7n7dEa90DSW4HvAJdFxKbudu0kVvZSTUXph/ZB/ZzDLg/RSawWz2F3qvYc9rZ9ko4je3LBhaVQJ7v1eP6chGpQRDwbEeMj4iTgbrLrBkTE9oj4YmSPzpgIvA1YmT57Pn3dDNxFFQ3xSBpI9o//WxHx3RTuz+WdCtVP7aunc9iVejmHXarWc9jb9kkaSbYU2/kR8fMU7tP5cxKqQaUZNZIGAFcB/y+930fZEklI+hCwPSKeTsNzB6b4QOCjZEN6hUvDEbcAz0TEP+U+6s/lnQrTX+2rs3PYqTo6h10dpyrPYW/bJ+ltwPeAKyPiv0s79/n8FTUjw6+yZ67cDaxn1yMtLgC+QDaD5WfA37Nr5YtRwAqyC4s/AA5P8X3JZsr9lOxC4jdIM66KfpHNIopUt2Xp9WGyR3LMJ+vJzQeG5cp8haz3t4Lc7BuySRhPpc+uL/1c6qF9dXgO15BNuGlP/66PrbNz+Ib2Ves57G37yP7wfS237zJgeF/Pn5ftMTOzwng4zszMCuMkZGZmhXESMjOzwjgJmZlZYZyEzGqMpAuVW/DUrJY5CZlVmKQzJYWkY/rhWP8baIuIjf1QNbPCeYq2WYUpe2z9QcD8iLi64OqYVRX3hMwqKK3H9QGym4wnpVhLeq7MHEnPSvpWaSFLSWsk/Z2kx9NzWY5J8X2VPVvqMUlPSJqY4ntJ+scU/6mkC1P8IEk/SgtlPiXpDwv5AZj1wEnIrLLOAL4fET8D2iSdmOInAJeR3Un/TrJEVfJyRJwI3Ah8KcW+AjwaEe8DxgH/mJZougB4NcXfB3wmLfdzDvBwRIwGjie7q92s6jgJmVXWJ4B70vY96T1k692ti4gdZAliVK5MaQHJpbn4eOAKScvIHow2GDgsxc9P8UVkS60cRfbgw09Juhp4b2QLZppVnaaiK2BWryQdAPwx8B5JAexFtkbXQ8DW3K6vs/v/xa2dxAV8LCJWdPgeAi6NiIc7+f4fBD4C3CHpHyPi9jffKrP+5Z6QWeWcBdweEYdHxKiIOBRYTbZgZG89DFyau3Z0Qi5+UVqVGUnvStePDgc2RMTNZCskn9jZQc2K5iRkVjmfIHvmSt53yK7X9NY0YCDwU0lPpfcA/wo8DTye4jeR9Z5agGWSniB7+uU3+vA9zSrOU7TNzKww7gmZmVlhnITMzKwwTkJmZlYYJyEzMyuMk5CZmRXGScjMzArjJGRmZoX5/2qjjwmixuEUAAAAAElFTkSuQmCC\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], + "source": [ + "yearly_incidence.plot(style='*')\n", + "plt.title('Incidence annuelle') # Ajout d'un titre\n", + "plt.xlabel('Années') # Nom de la grandeur en abscisse\n", + "plt.ylabel('Incidence') # Nom de la grandeur en ordonnée\n", + "plt.grid()" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "Une liste triée permet de plus facilement répérer les valeurs les plus élevées (à la fin)." + ] + }, + { + "cell_type": "code", + "execution_count": 90, + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "2002 516689\n", + "2018 542312\n", + "2017 551041\n", + "1996 564901\n", + "2019 584066\n", + "2015 604382\n", + "2000 617597\n", + "2001 619041\n", + "2012 624573\n", + "2005 628464\n", + "2006 632833\n", + "2011 642368\n", + "1993 643387\n", + "1995 652478\n", + "1994 661409\n", + "1998 677775\n", + "1997 683434\n", + "2014 685769\n", + "2013 698332\n", + "2007 717352\n", + "2008 749478\n", + "1999 756456\n", + "2003 758363\n", + "2004 777388\n", + "2016 782114\n", + "2010 829911\n", + "1992 832939\n", + "2009 842373\n", + "dtype: int64" + ] + }, + "execution_count": 90, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "yearly_incidence.sort_values()" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "Enfin, un histogramme montre bien que les épidémies fortes, qui touchent moins de 1,5% de la population française, sont assez rares: il y en eu 4 au cours des 20 dernières années." + ] + }, + { + "cell_type": "code", + "execution_count": 91, + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "" + ] + }, + "execution_count": 91, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXwAAAEICAYAAABcVE8dAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAFGtJREFUeJzt3XuQJWV9h/Hnxy4oMLggDKgLukbIRGFVZL2gVTqDxhKXaGm8IZpotNYqFTFlyqCoaLytF7zESyobJZB4nShUkDUaDI6oUZQV4oI4amBVlighKjCI4Oovf3QvGadmdmZO95lzOu/zqdraPrfu73mn53v69Ok+E5mJJOn/v70GHUCStDIsfEkqhIUvSYWw8CWpEBa+JBXCwpekQlj4klQIC1+SCmHhS1IhVq/kwg488MA88sgjV3KRrbn11lvZf//9Bx1j2bqaG8w+KF3N3tXcsHj2bdu23ZiZo02Xs6KFf9hhh3HZZZet5CJbMzU1xfj4+KBjLFtXc4PZB6Wr2buaGxbPHhE/bGM57tKRpEJY+JJUCAtfkgph4UtSISx8SSpEo8KPiLGIuGLWv5sj4uVthZMktafRYZmZOQ08GCAiVgE7gfNbyCVJalmbu3QeC/xnZrZyvKgkqV3R1t+0jYizgW9l5vvnXL8J2AQwOjp63OTkZCvLW2kzMzOMjIwMOsay9Zp7+86b+pBmcevXrrlzuqtjDmYfhK7mhsWzT0xMbMvMDU2X00rhR8Q+wPXA0Zn504XuNzY2ltPT042XNwhdPYuv19zrTt/afpgl2LF5453TXR1zMPsgdDU3LOlM21YKv61dOidSbd0vWPaSpMFqq/BPBj7e0rwkSX3QuPAjYj/gD4HzmseRJPVL42/LzMxfAge3kEWS1EeeaStJhbDwJakQFr4kFcLCl6RCWPiSVAgLX5IKYeFLUiEsfEkqhIUvSYWw8CWpEBa+JBXCwpekQlj4klQIC1+SCmHhS1IhLHxJKoSFL0mFsPAlqRAWviQVwsKXpEJY+JJUiMaFHxEHRsSnIuK7EXF1RBzfRjBJUrtWtzCP9wKfy8ynRcQ+wH4tzFOS1LJGhR8RdwMeDTwPIDPvAO5oHkuS1LbIzN4fHPFgYAvwHeBBwDbgtMy8ddZ9NgGbAEZHR4+bnJxsFHhQZmZmGBkZGXSMZes19/adN/UhzeLWr11z53RXxxzMPghdzQ2LZ5+YmNiWmRuaLqdp4W8Avg48KjMvjYj3Ajdn5mvnu//Y2FhOT0/3vLxBmpqaYnx8fNAxlq3X3OtO39p+mCXYsXnjndNdHXMw+yB0NTcsnj0iWin8ph/aXgdcl5mX1pc/BTyk4TwlSX3QqPAz8yfAjyNirL7qsVS7dyRJQ6aNo3ROBT5aH6FzDfD8FuYpSWpZ48LPzCuAxvuWJEn95Zm2klQIC1+SCmHhS1IhLHxJKoSFL0mFsPAlqRAWviQVwsKXpEJY+JJUCAtfkgph4UtSISx8SSqEhS9JhbDwJakQFr4kFcLCl6RCWPiSVAgLX5IKYeFLUiEsfEkqROM/Yh4RO4BbgN8AuzLTP2guSUOoceHXJjLzxpbmJUnqA3fpSFIhIjObzSDiWuDnQAJ/m5lb5ty+CdgEMDo6etzk5GSj5Q3KzMwMIyMjg46xbL3m3r7zpj6kWdz6tWvunO7qmIPZB6GruWHx7BMTE9va2F3eRuHfKzOvj4hDgYuAUzPzkvnuOzY2ltPT042WNyhTU1OMj48POsay9Zp73elb2w+zBDs2b7xzuqtjDmYfhK7mhsWzR0Qrhd94l05mXl//fwNwPvCwpvOUJLWvUeFHxP4RccDuaeDxwJVtBJMktavpUTqHAedHxO55fSwzP9c4lSSpdY0KPzOvAR7UUhZJUh95WKYkFcLCl6RCWPiSVAgLX5IKYeFLUiEsfEkqhIUvSYWw8CWpEBa+JBXCwpekQlj4klQIC1+SCmHhS1IhLHxJKoSFL0mFsPAlqRAWviQVwsKXpEJY+JJUCAtfkgph4UtSIRoXfkSsiojLI+LCNgJJkvqjjS3804CrW5iPJKmPGhV+RBwObAQ+1E4cSVK/RGb2/uCITwFvBQ4A/iIzT5rnPpuATQCjo6PHTU5O9ry8QZqZmWFkZGTQMZat19zbd97UhzSLW792zZ3TXR1zMPsgdDU3LJ59YmJiW2ZuaLqc1b0+MCJOAm7IzG0RMb7Q/TJzC7AFYGxsLMfHF7zrUJuamqKL2XvN/bzTt7YfZgl2nDJ+53RXxxzMPghdzQ0rl73JLp1HAU+KiB3AJ4ATIuIjraSSJLWu58LPzFdl5uGZuQ54FnBxZj6ntWSSpFZ5HL4kFaLnffizZeYUMNXGvCRJ/eEWviQVwsKXpEJY+JJUCAtfkgph4UtSISx8SSqEhS9JhbDwJakQFr4kFcLCl6RCWPiSVAgLX5IKYeFLUiEsfEkqhIUvSYWw8CWpEBa+JBXCwpekQlj4klQIC1+SCtGo8CPirhHxjYj4j4i4KiLe0FYwSVK7Vjd8/O3ACZk5ExF7A1+JiH/JzK+3kE2S1KJGhZ+ZCczUF/eu/2XTUJKk9jXehx8RqyLiCuAG4KLMvLR5LElS26LaSG9hRhEHAucDp2bmlbOu3wRsAhgdHT1ucnKyleWttJmZGUZGRgYdY9l6zb195019SLM8h+0LP71t0Cl6s5zs69eu6W+YPZjv57wS496P57yUdX2Q6/WenvNi2ScmJrZl5oamGVorfICIOBO4NTPfOd/tY2NjOT093dryVtLU1BTj4+ODjrFsveZed/rW9sMs0yvW7+Ks7U0/ZhqM5WTfsXljn9MsbL6f80qMez+e81LW9UGu13t6zotlj4hWCr/pUTqj9ZY9EbEv8Djgu01DSZLa1/Rl/J7AuRGxiurFYzIzL2weS5LUtqZH6XwbOLalLJKkPvJMW0kqhIUvSYWw8CWpEBa+JBXCwpekQlj4klQIC1+SCmHhS1IhLHxJKoSFL0mFsPAlqRAWviQVwsKXpEJY+JJUCAtfkgph4UtSISx8SSqEhS9JhbDwJakQFr4kFcLCl6RCNCr8iDgiIr4YEVdHxFURcVpbwSRJ7Vrd8PG7gFdk5rci4gBgW0RclJnfaSGbJKlFjbbwM/O/MvNb9fQtwNXA2jaCSZLaFZnZzowi1gGXAMdk5s2zrt8EbAIYHR09bnJysqf5b995U/OQDRy2L/z0toFG6ElXc0M52devXdPfMHsw3+9VV8d92HPv6ec8MzPDyMjIgrdPTExsy8wNTTO0UvgRMQJ8CXhzZp630P3GxsZyenq6p2WsO31rj+na8Yr1uzhre9M9YCuvq7mhnOw7Nm/sc5qFzfd71dVxH/bce/o5T01NMT4+vuDtEdFK4Tc+Sici9gY+DXx0T2UvSRqspkfpBPBh4OrMfFc7kSRJ/dB0C/9RwHOBEyLiivrfE1vIJUlqWaMdXpn5FSBayiJJ6iPPtJWkQlj4klQIC1+SCmHhS1IhLHxJKoSFL0mFsPAlqRAWviQVwsKXpEJY+JJUCAtfkgph4UtSISx8SSqEhS9JhbDwJakQFr4kFcLCl6RCWPiSVAgLX5IKYeFLUiEsfEkqRKPCj4izI+KGiLiyrUCSpP5ouoV/DvCEFnJIkvqsUeFn5iXAz1rKIknqo8jMZjOIWAdcmJnHLHD7JmATwOjo6HGTk5M9LWf7zpt6TNiOw/aFn9420Ag96WpuMPugdDX7sOdev3bNgrfNzMwwMjKy4O0TExPbMnND0wx9L/zZxsbGcnp6uqflrDt9a0+Pa8sr1u/irO2rB5qhF13NDWYflK5mH/bcOzZvXPC2qakpxsfHF7w9IlopfI/SkaRCWPiSVIimh2V+HPgaMBYR10XEC9qJJUlqW6MdXpl5cltBJEn95S4dSSqEhS9JhbDwJakQFr4kFcLCl6RCWPiSVAgLX5IKYeFLUiEsfEkqhIUvSYWw8CWpEBa+JBXCwpekQlj4klQIC1+SCmHhS1IhLHxJKoSFL0mFsPAlqRAWviQVonHhR8QTImI6In4QEae3EUqS1L5GhR8Rq4APACcCDwBOjogHtBFMktSuplv4DwN+kJnXZOYdwCeAJzePJUlqW2Rm7w+OeBrwhMx8YX35ucDDM/Ols+6zCdhUXzwGuLL3uAN1CHDjoEP0oKu5weyD0tXsXc0Ni2e/T2aONl3I6oaPj3mu+51XkMzcAmwBiIjLMnNDw2UORFezdzU3mH1Qupq9q7lh5bI33aVzHXDErMuHA9c3nKckqQ+aFv43gaMi4r4RsQ/wLOCC5rEkSW1rtEsnM3dFxEuBzwOrgLMz86o9PGRLk+UNWFezdzU3mH1Qupq9q7lhhbI3+tBWktQdnmkrSYWw8CWpEBa+JBVi6As/Ik6IiPsOOsdydTU3mH1Qupq9q7mhvOxD+6Ft/Z08nwB+AfwWODMzvzTYVIvram4w+6B0NXtXc0O52YdmCz8iDo+Iu8266pnApzPz0VRP7uSIOH4w6RbW1dxg9kHpavau5gaz7zbwwo+I+0fEZ4GvAH8VEbu/fO1XwH719CTV90w8IiLm+zqHFdfV3GD2Qelq9q7mBrPPNZDCj4j9Z118MHBdZq4DLgbeWV//M+D2iDggM38GfA84DFi3glF/R1dzg9kHpavZu5obzL4nK1b4EXFQRJwTEd8ENkfEaP2K9EDgqxERmXkB8IuI2Ej1JA4A1tez+D7VN8rdsVKZu5zb7GYvJbfZl559JbfwHw3sAp5I9S2brwbuVme4R/7fp8fnAs8GvgHcQvXHVcjMrwEnADevYGbobm4wu9mXp6u5wexLy56Zrf6j+k6dFwFfovoe/EPq6yeBl9XT9wU217c/lGof1ar6thHgv+v5rAWuBl4K/D3wQWC/tjN3ObfZzV5KbrM3z96PLfyTgCcBbwCOB95eX38R8Mh6+sfAl4ETM/ObVK9qEwCZOQNcCjw0M3cCz6XaP/UT4DWZ+cs2w876oOOPupR7jk6NOTjug8jumLuu9/RtmfU+pYyIh1K9xfgysDUzbwd+H7gmMy+OiGuBd0TE44FtwFMi4pDMvDEivg/cGhH3Bt4HPCciDqX6Tv3/oXrbQmZeBlzWS8495N8AvJDqbdE7gBuA3xv23HX2To55nd1xd11fSuZO90v9HIZy3Je9hR8Rd6t/GOPA2VSHCD0OeGt9l98C34uIfTPzWqq3HQ+sn/j1VMeQAvyG6q3JXsCnqb4e9BTgOGBLZv52udkWyT0SEXeNiHPr5V0LvDczb4iIvaheSYcud5394Pr/RwLn0JExrzMfFNVxxFuA8+jWuB8aEQfXv7zn0pFxj4hDIuIhEfF+4Hy6Neb71P3yGKpdFZ0Y8zr73hGxb0Scw7B2zBL3Pe0H/Cnwb1QH/AP8OfCSevog4NvAsXXozcC6+raT6sCH1NPbgTVUHzJ8Fthn1nL2arKPbA+5LwY+WV/3NuBFs+6zuv7/pcCbhiT37jOgn061D28KGANe3oExn539q8AXqFbqoR/3ep77A8+jept9E7CxI+O+O/cXqLYCn9qhMd8beDHwGeBvgCOB04Z9zOdkv5Bqa/zoOt9QjvuiW/gRsTdwFfA04B2Z+cf1Tcfuvk9m/hz4Z+BlVL8ohwL3r2++BHgMcEdmXgh8GPgU8AGqraZfz5pPa6+6c3K/PTN3v3puBx4QEW+tX4n/LCLuDnwOuMegc9fzy4hYAzwDeE9mjmfmNNUr/O77DN2Yz5P93Zn5uMz8Nh0Y94g4kmr3wWOB1wA7gR9RfZC2ql7m0I37nNxnUO3T/WGdZ/0wj3ntJVQl9x6qv4n91Pryb+plDt2Yz5P93VTryFOo1pmxiNg8dOO+xFex84BT5lz3TODSWZfvBVxfT7+E6pTfg4B9qV657z3rvoe0/Uq7jNz3rrN9EjiZaiX73DDlrpf1YuCN9fTureanDvuYz81eX94XuGedb3JYx53qF/Yusy6fTbXB8ORhHvd5cn+Iamv/0I6s658B/qSefgFwat0v3xjWMV8g+/OBV9Ud88lhHPelPqmTqA72P4tq98LrgPtQnfF16Kz7XQQ8vJ5+E/CvVG8v/3KlfgB7yP1F4PVUH3ysnXWfvevndEJ9+c2Dzl3neEad6xTgW1Sv+E+m2s1wyKz7DdWYz5N9G/ARqiMUZq/YQznus/IdTLUP+XFU+1J/RnVM9NCO+5zcj68vd2FdfzHV7uJ/ojpK5UNUu0Z+DowO85jPyf4jqo2E+wNHDOO4L+eJfZ7qhIAjqF65TgO+RlX+Ady9/kHde9aTPAa466B+GPPk/ijV6clHzbr9HlQfhD5wyHIfVY/vX1NtDTwbeBfVcbivrEtoWMd8bvanU23Znzjs4z7neVwOPLWe/gjwtnp6KMd9Tu4n1dN7zbp+aMec6oiWs6n2Z78eOBOYrn93h3Zdnyf7GcDfAY8YxnFfzlE6T8nMt2Tmj4G3UB0e9T6qtyUXUL3KrcrMHwFk5q8z88rM/NUyltEPs3O/lWrwD4yIB0XEGVS7fX6Z1X7mYcr9I6pTpVdntQ/zwvq6C6k+jB7mMZ+bfSvwHeDYiPiDIR936iMqoPrA/Mh6+s3VTfEZhnTcF8i9V33EzqsZ4jEHHgBMZeaNVO9QVgEfZ/j7BX43+7nAbcD9IuLoYVvXl3wcfv7uQf2/oCr812bmxyLiOcBVmXl52wGbmpP7FqoTFa6meiXem+pIgGHMfXtEbKZ66wrVVs7xwObMvHzIx3xu9n2o3ua+keqoqX0Y0nGH6gOyiNj9bYQ/qK+7GnhlRJwCfGcYs8/JvaO+bldEPItq/RnKMY+IVcB1VF8x8A/ArVRnmb44M68d5nV9nuy/onqH+3aq3T13YYjGfcl/ACUi7gI8gersrqOpDp/6YGbu6l+85ubJvSUz3z3YVEsXEW+mOmHjWKr9fmfWW81Db072LwCvzP6ehdmqiJgGXpeZn9x9MtCgMy3F7NyDzrJUEXE/qkMU76BaZ84H3pTV2aVDbZ7sFwBnZOaKfxHbYpb1F68i4kVUJz7845C8lVqSrubeLSLGgB+afWXMOtPzwVSH9u7qQtl3Nfdu9ZmlRwH/npm3DTrPcnQl+9D+iUNJUrsG/hevJEkrw8KXpEJY+JJUCAtfkgph4UtSISx8SSqEhS9JhfhfHCg3qEW8+c0AAAAASUVORK5CYII=\n", + "text/plain": [ + "
" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + } + ], + "source": [ + "yearly_incidence.hist(xrot=20)" + ] + }, + { + "cell_type": "code", + "execution_count": 95, + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "842373" + ] + }, + "execution_count": 95, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "max(yearly_incidence)" + ] + }, + { + "cell_type": "code", + "execution_count": 93, + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "516689" + ] + }, + "execution_count": 93, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "min(yearly_incidence)" + ] + }, { "cell_type": "code", "execution_count": null, -- 2.18.1