Code.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Vaccination data: https://www.kaggle.com/gpreda/covid-world-vaccination-progress\n",
    "# New covid cases in Lithuania: https://open-data-ls-osp-sdg.hub.arcgis.com/datasets/ba35de03e111430f88a86f7d1f351de6_0/data\n",
    "# New covid cases worldwide: https://github.com/CSSEGISandData/COVID-19/tree/master/csse_covid_19_data/csse_covid_19_daily_reports"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 29,
   "metadata": {},
   "outputs": [],
   "source": [
    "import pandas as pd\n",
    "import numpy as np"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Prepare confirmed cases data"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 30,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Country/Region</th>\n",
       "      <th>1/22/20</th>\n",
       "      <th>1/23/20</th>\n",
       "      <th>1/24/20</th>\n",
       "      <th>1/25/20</th>\n",
       "      <th>1/26/20</th>\n",
       "      <th>1/27/20</th>\n",
       "      <th>1/28/20</th>\n",
       "      <th>1/29/20</th>\n",
       "      <th>1/30/20</th>\n",
       "      <th>...</th>\n",
       "      <th>4/11/21</th>\n",
       "      <th>4/12/21</th>\n",
       "      <th>4/13/21</th>\n",
       "      <th>4/14/21</th>\n",
       "      <th>4/15/21</th>\n",
       "      <th>4/16/21</th>\n",
       "      <th>4/17/21</th>\n",
       "      <th>4/18/21</th>\n",
       "      <th>4/19/21</th>\n",
       "      <th>4/20/21</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>Afghanistan</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>...</td>\n",
       "      <td>57160</td>\n",
       "      <td>57242</td>\n",
       "      <td>57364</td>\n",
       "      <td>57492</td>\n",
       "      <td>57534</td>\n",
       "      <td>57612</td>\n",
       "      <td>57721</td>\n",
       "      <td>57793</td>\n",
       "      <td>57898</td>\n",
       "      <td>58037</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>Albania</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>...</td>\n",
       "      <td>128393</td>\n",
       "      <td>128518</td>\n",
       "      <td>128752</td>\n",
       "      <td>128959</td>\n",
       "      <td>129128</td>\n",
       "      <td>129307</td>\n",
       "      <td>129456</td>\n",
       "      <td>129594</td>\n",
       "      <td>129694</td>\n",
       "      <td>129842</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>Algeria</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>...</td>\n",
       "      <td>118516</td>\n",
       "      <td>118645</td>\n",
       "      <td>118799</td>\n",
       "      <td>118975</td>\n",
       "      <td>119142</td>\n",
       "      <td>119323</td>\n",
       "      <td>119486</td>\n",
       "      <td>119642</td>\n",
       "      <td>119805</td>\n",
       "      <td>119992</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>Andorra</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>...</td>\n",
       "      <td>12545</td>\n",
       "      <td>12581</td>\n",
       "      <td>12614</td>\n",
       "      <td>12641</td>\n",
       "      <td>12641</td>\n",
       "      <td>12712</td>\n",
       "      <td>12771</td>\n",
       "      <td>12805</td>\n",
       "      <td>12805</td>\n",
       "      <td>12874</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>Angola</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>...</td>\n",
       "      <td>23457</td>\n",
       "      <td>23549</td>\n",
       "      <td>23697</td>\n",
       "      <td>23841</td>\n",
       "      <td>23951</td>\n",
       "      <td>24122</td>\n",
       "      <td>24300</td>\n",
       "      <td>24389</td>\n",
       "      <td>24518</td>\n",
       "      <td>24661</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>269</th>\n",
       "      <td>Vietnam</td>\n",
       "      <td>0</td>\n",
       "      <td>2</td>\n",
       "      <td>2</td>\n",
       "      <td>2</td>\n",
       "      <td>2</td>\n",
       "      <td>2</td>\n",
       "      <td>2</td>\n",
       "      <td>2</td>\n",
       "      <td>2</td>\n",
       "      <td>...</td>\n",
       "      <td>2693</td>\n",
       "      <td>2705</td>\n",
       "      <td>2714</td>\n",
       "      <td>2733</td>\n",
       "      <td>2758</td>\n",
       "      <td>2772</td>\n",
       "      <td>2781</td>\n",
       "      <td>2785</td>\n",
       "      <td>2791</td>\n",
       "      <td>2800</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>270</th>\n",
       "      <td>West Bank and Gaza</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>...</td>\n",
       "      <td>268132</td>\n",
       "      <td>270856</td>\n",
       "      <td>272767</td>\n",
       "      <td>274690</td>\n",
       "      <td>276407</td>\n",
       "      <td>278135</td>\n",
       "      <td>279753</td>\n",
       "      <td>280741</td>\n",
       "      <td>282270</td>\n",
       "      <td>284280</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>271</th>\n",
       "      <td>Yemen</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>...</td>\n",
       "      <td>5357</td>\n",
       "      <td>5446</td>\n",
       "      <td>5507</td>\n",
       "      <td>5582</td>\n",
       "      <td>5657</td>\n",
       "      <td>5715</td>\n",
       "      <td>5770</td>\n",
       "      <td>5812</td>\n",
       "      <td>5858</td>\n",
       "      <td>5918</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>272</th>\n",
       "      <td>Zambia</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>...</td>\n",
       "      <td>90029</td>\n",
       "      <td>90064</td>\n",
       "      <td>90218</td>\n",
       "      <td>90389</td>\n",
       "      <td>90532</td>\n",
       "      <td>90750</td>\n",
       "      <td>90844</td>\n",
       "      <td>90918</td>\n",
       "      <td>90942</td>\n",
       "      <td>91042</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>273</th>\n",
       "      <td>Zimbabwe</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>0</td>\n",
       "      <td>...</td>\n",
       "      <td>37288</td>\n",
       "      <td>37307</td>\n",
       "      <td>37330</td>\n",
       "      <td>37369</td>\n",
       "      <td>37422</td>\n",
       "      <td>37534</td>\n",
       "      <td>37699</td>\n",
       "      <td>37751</td>\n",
       "      <td>37859</td>\n",
       "      <td>37875</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>274 rows × 456 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "         Country/Region  1/22/20  1/23/20  1/24/20  1/25/20  1/26/20  1/27/20  \\\n",
       "0           Afghanistan        0        0        0        0        0        0   \n",
       "1               Albania        0        0        0        0        0        0   \n",
       "2               Algeria        0        0        0        0        0        0   \n",
       "3               Andorra        0        0        0        0        0        0   \n",
       "4                Angola        0        0        0        0        0        0   \n",
       "..                  ...      ...      ...      ...      ...      ...      ...   \n",
       "269             Vietnam        0        2        2        2        2        2   \n",
       "270  West Bank and Gaza        0        0        0        0        0        0   \n",
       "271               Yemen        0        0        0        0        0        0   \n",
       "272              Zambia        0        0        0        0        0        0   \n",
       "273            Zimbabwe        0        0        0        0        0        0   \n",
       "\n",
       "     1/28/20  1/29/20  1/30/20  ...  4/11/21  4/12/21  4/13/21  4/14/21  \\\n",
       "0          0        0        0  ...    57160    57242    57364    57492   \n",
       "1          0        0        0  ...   128393   128518   128752   128959   \n",
       "2          0        0        0  ...   118516   118645   118799   118975   \n",
       "3          0        0        0  ...    12545    12581    12614    12641   \n",
       "4          0        0        0  ...    23457    23549    23697    23841   \n",
       "..       ...      ...      ...  ...      ...      ...      ...      ...   \n",
       "269        2        2        2  ...     2693     2705     2714     2733   \n",
       "270        0        0        0  ...   268132   270856   272767   274690   \n",
       "271        0        0        0  ...     5357     5446     5507     5582   \n",
       "272        0        0        0  ...    90029    90064    90218    90389   \n",
       "273        0        0        0  ...    37288    37307    37330    37369   \n",
       "\n",
       "     4/15/21  4/16/21  4/17/21  4/18/21  4/19/21  4/20/21  \n",
       "0      57534    57612    57721    57793    57898    58037  \n",
       "1     129128   129307   129456   129594   129694   129842  \n",
       "2     119142   119323   119486   119642   119805   119992  \n",
       "3      12641    12712    12771    12805    12805    12874  \n",
       "4      23951    24122    24300    24389    24518    24661  \n",
       "..       ...      ...      ...      ...      ...      ...  \n",
       "269     2758     2772     2781     2785     2791     2800  \n",
       "270   276407   278135   279753   280741   282270   284280  \n",
       "271     5657     5715     5770     5812     5858     5918  \n",
       "272    90532    90750    90844    90918    90942    91042  \n",
       "273    37422    37534    37699    37751    37859    37875  \n",
       "\n",
       "[274 rows x 456 columns]"
      ]
     },
     "execution_count": 30,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "url = \"https://raw.githubusercontent.com/CSSEGISandData/COVID-19/master/csse_covid_19_data/csse_covid_19_time_series/time_series_covid19_confirmed_global.csv\"\n",
    "dfCases = pd.read_csv(url, index_col=0)\n",
    "dfCases = dfCases.drop(columns=['Lat', 'Long'])\n",
    "dfCases = dfCases.reset_index(drop = True)\n",
    "dfCases"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 75,
   "metadata": {},
   "outputs": [],
   "source": [
    "def prepareCasesData(country, countryPop):\n",
    "    data = dfCases[dfCases['Country/Region'] == country]\n",
    "    data = data.melt(id_vars=[\"Country/Region\"], \n",
    "        var_name=\"Date\", \n",
    "        value_name=\"Confirmed Cases\")\n",
    "    data['New Cases'] = data['Confirmed Cases'].diff() # Calculating new cases by substracting nth row value from nth-1 row value\n",
    "    data['New Cases per pop.'] = data['New Cases'] / countryPop\n",
    "    data['MA7'] = data['New Cases'].rolling(7).mean() # 7 day Moving Average for new cases \n",
    "    data['MA7'].fillna(0, inplace = True)\n",
    "    data['New Cases'].fillna(0, inplace = True)\n",
    "    data['Date'] = pd.to_datetime(data.Date)\n",
    "    \n",
    "    return data"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 76,
   "metadata": {},
   "outputs": [],
   "source": [
    "casesLT = prepareCasesData('Lithuania', 2722289)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 77,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Country/Region</th>\n",
       "      <th>Date</th>\n",
       "      <th>Confirmed Cases</th>\n",
       "      <th>New Cases</th>\n",
       "      <th>New Cases per pop.</th>\n",
       "      <th>MA7</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>Lithuania</td>\n",
       "      <td>2020-01-22</td>\n",
       "      <td>0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>Lithuania</td>\n",
       "      <td>2020-01-23</td>\n",
       "      <td>0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>Lithuania</td>\n",
       "      <td>2020-01-24</td>\n",
       "      <td>0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>Lithuania</td>\n",
       "      <td>2020-01-25</td>\n",
       "      <td>0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>Lithuania</td>\n",
       "      <td>2020-01-26</td>\n",
       "      <td>0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>450</th>\n",
       "      <td>Lithuania</td>\n",
       "      <td>2021-04-16</td>\n",
       "      <td>231601</td>\n",
       "      <td>1139.0</td>\n",
       "      <td>0.001139</td>\n",
       "      <td>1041.714286</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>451</th>\n",
       "      <td>Lithuania</td>\n",
       "      <td>2021-04-17</td>\n",
       "      <td>232829</td>\n",
       "      <td>1228.0</td>\n",
       "      <td>0.001228</td>\n",
       "      <td>1074.428571</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>452</th>\n",
       "      <td>Lithuania</td>\n",
       "      <td>2021-04-18</td>\n",
       "      <td>233631</td>\n",
       "      <td>802.0</td>\n",
       "      <td>0.000802</td>\n",
       "      <td>1065.428571</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>453</th>\n",
       "      <td>Lithuania</td>\n",
       "      <td>2021-04-19</td>\n",
       "      <td>234232</td>\n",
       "      <td>601.0</td>\n",
       "      <td>0.000601</td>\n",
       "      <td>1064.142857</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>454</th>\n",
       "      <td>Lithuania</td>\n",
       "      <td>2021-04-20</td>\n",
       "      <td>235383</td>\n",
       "      <td>1151.0</td>\n",
       "      <td>0.001151</td>\n",
       "      <td>1070.571429</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>455 rows × 6 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "    Country/Region       Date  Confirmed Cases  New Cases  New Cases per pop.  \\\n",
       "0        Lithuania 2020-01-22                0        0.0                 NaN   \n",
       "1        Lithuania 2020-01-23                0        0.0            0.000000   \n",
       "2        Lithuania 2020-01-24                0        0.0            0.000000   \n",
       "3        Lithuania 2020-01-25                0        0.0            0.000000   \n",
       "4        Lithuania 2020-01-26                0        0.0            0.000000   \n",
       "..             ...        ...              ...        ...                 ...   \n",
       "450      Lithuania 2021-04-16           231601     1139.0            0.001139   \n",
       "451      Lithuania 2021-04-17           232829     1228.0            0.001228   \n",
       "452      Lithuania 2021-04-18           233631      802.0            0.000802   \n",
       "453      Lithuania 2021-04-19           234232      601.0            0.000601   \n",
       "454      Lithuania 2021-04-20           235383     1151.0            0.001151   \n",
       "\n",
       "             MA7  \n",
       "0       0.000000  \n",
       "1       0.000000  \n",
       "2       0.000000  \n",
       "3       0.000000  \n",
       "4       0.000000  \n",
       "..           ...  \n",
       "450  1041.714286  \n",
       "451  1074.428571  \n",
       "452  1065.428571  \n",
       "453  1064.142857  \n",
       "454  1070.571429  \n",
       "\n",
       "[455 rows x 6 columns]"
      ]
     },
     "execution_count": 77,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "casesLT"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Prepare vaccines data"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 34,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>country</th>\n",
       "      <th>iso_code</th>\n",
       "      <th>date</th>\n",
       "      <th>total_vaccinations</th>\n",
       "      <th>people_vaccinated</th>\n",
       "      <th>people_fully_vaccinated</th>\n",
       "      <th>daily_vaccinations_raw</th>\n",
       "      <th>daily_vaccinations</th>\n",
       "      <th>total_vaccinations_per_hundred</th>\n",
       "      <th>people_vaccinated_per_hundred</th>\n",
       "      <th>people_fully_vaccinated_per_hundred</th>\n",
       "      <th>daily_vaccinations_per_million</th>\n",
       "      <th>vaccines</th>\n",
       "      <th>source_name</th>\n",
       "      <th>source_website</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>Afghanistan</td>\n",
       "      <td>AFG</td>\n",
       "      <td>2021-02-22</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.00</td>\n",
       "      <td>0.00</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>Oxford/AstraZeneca</td>\n",
       "      <td>Government of Afghanistan</td>\n",
       "      <td>http://www.xinhuanet.com/english/asiapacific/2...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>Afghanistan</td>\n",
       "      <td>AFG</td>\n",
       "      <td>2021-02-23</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>1367.0</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>35.0</td>\n",
       "      <td>Oxford/AstraZeneca</td>\n",
       "      <td>Government of Afghanistan</td>\n",
       "      <td>http://www.xinhuanet.com/english/asiapacific/2...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>Afghanistan</td>\n",
       "      <td>AFG</td>\n",
       "      <td>2021-02-24</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>1367.0</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>35.0</td>\n",
       "      <td>Oxford/AstraZeneca</td>\n",
       "      <td>Government of Afghanistan</td>\n",
       "      <td>http://www.xinhuanet.com/english/asiapacific/2...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>Afghanistan</td>\n",
       "      <td>AFG</td>\n",
       "      <td>2021-02-25</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>1367.0</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>35.0</td>\n",
       "      <td>Oxford/AstraZeneca</td>\n",
       "      <td>Government of Afghanistan</td>\n",
       "      <td>http://www.xinhuanet.com/english/asiapacific/2...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>Afghanistan</td>\n",
       "      <td>AFG</td>\n",
       "      <td>2021-02-26</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>1367.0</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>NaN</td>\n",
       "      <td>35.0</td>\n",
       "      <td>Oxford/AstraZeneca</td>\n",
       "      <td>Government of Afghanistan</td>\n",
       "      <td>http://www.xinhuanet.com/english/asiapacific/2...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>12509</th>\n",
       "      <td>Zimbabwe</td>\n",
       "      <td>ZWE</td>\n",
       "      <td>2021-04-12</td>\n",
       "      <td>234579.0</td>\n",
       "      <td>205275.0</td>\n",
       "      <td>29304.0</td>\n",
       "      <td>11087.0</td>\n",
       "      <td>12724.0</td>\n",
       "      <td>1.58</td>\n",
       "      <td>1.38</td>\n",
       "      <td>0.20</td>\n",
       "      <td>856.0</td>\n",
       "      <td>Sinopharm/Beijing</td>\n",
       "      <td>Ministry of Health</td>\n",
       "      <td>https://twitter.com/MoHCCZim/status/1383174384...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>12510</th>\n",
       "      <td>Zimbabwe</td>\n",
       "      <td>ZWE</td>\n",
       "      <td>2021-04-13</td>\n",
       "      <td>248355.0</td>\n",
       "      <td>218516.0</td>\n",
       "      <td>29839.0</td>\n",
       "      <td>13776.0</td>\n",
       "      <td>12246.0</td>\n",
       "      <td>1.67</td>\n",
       "      <td>1.47</td>\n",
       "      <td>0.20</td>\n",
       "      <td>824.0</td>\n",
       "      <td>Sinopharm/Beijing</td>\n",
       "      <td>Ministry of Health</td>\n",
       "      <td>https://twitter.com/MoHCCZim/status/1383174384...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>12511</th>\n",
       "      <td>Zimbabwe</td>\n",
       "      <td>ZWE</td>\n",
       "      <td>2021-04-14</td>\n",
       "      <td>262128.0</td>\n",
       "      <td>231632.0</td>\n",
       "      <td>30496.0</td>\n",
       "      <td>13773.0</td>\n",
       "      <td>11816.0</td>\n",
       "      <td>1.76</td>\n",
       "      <td>1.56</td>\n",
       "      <td>0.21</td>\n",
       "      <td>795.0</td>\n",
       "      <td>Sinopharm/Beijing</td>\n",
       "      <td>Ministry of Health</td>\n",
       "      <td>https://twitter.com/MoHCCZim/status/1383174384...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>12512</th>\n",
       "      <td>Zimbabwe</td>\n",
       "      <td>ZWE</td>\n",
       "      <td>2021-04-15</td>\n",
       "      <td>285881.0</td>\n",
       "      <td>252932.0</td>\n",
       "      <td>32949.0</td>\n",
       "      <td>23753.0</td>\n",
       "      <td>13172.0</td>\n",
       "      <td>1.92</td>\n",
       "      <td>1.70</td>\n",
       "      <td>0.22</td>\n",
       "      <td>886.0</td>\n",
       "      <td>Sinopharm/Beijing</td>\n",
       "      <td>Ministry of Health</td>\n",
       "      <td>https://twitter.com/MoHCCZim/status/1383174384...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>12513</th>\n",
       "      <td>Zimbabwe</td>\n",
       "      <td>ZWE</td>\n",
       "      <td>2021-04-16</td>\n",
       "      <td>304701.0</td>\n",
       "      <td>269732.0</td>\n",
       "      <td>34969.0</td>\n",
       "      <td>18820.0</td>\n",
       "      <td>14071.0</td>\n",
       "      <td>2.05</td>\n",
       "      <td>1.81</td>\n",
       "      <td>0.24</td>\n",
       "      <td>947.0</td>\n",
       "      <td>Sinopharm/Beijing</td>\n",
       "      <td>Ministry of Health</td>\n",
       "      <td>https://twitter.com/MoHCCZim/status/1383174384...</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>12514 rows × 15 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "           country iso_code        date  total_vaccinations  \\\n",
       "0      Afghanistan      AFG  2021-02-22                 0.0   \n",
       "1      Afghanistan      AFG  2021-02-23                 NaN   \n",
       "2      Afghanistan      AFG  2021-02-24                 NaN   \n",
       "3      Afghanistan      AFG  2021-02-25                 NaN   \n",
       "4      Afghanistan      AFG  2021-02-26                 NaN   \n",
       "...            ...      ...         ...                 ...   \n",
       "12509     Zimbabwe      ZWE  2021-04-12            234579.0   \n",
       "12510     Zimbabwe      ZWE  2021-04-13            248355.0   \n",
       "12511     Zimbabwe      ZWE  2021-04-14            262128.0   \n",
       "12512     Zimbabwe      ZWE  2021-04-15            285881.0   \n",
       "12513     Zimbabwe      ZWE  2021-04-16            304701.0   \n",
       "\n",
       "       people_vaccinated  people_fully_vaccinated  daily_vaccinations_raw  \\\n",
       "0                    0.0                      NaN                     NaN   \n",
       "1                    NaN                      NaN                     NaN   \n",
       "2                    NaN                      NaN                     NaN   \n",
       "3                    NaN                      NaN                     NaN   \n",
       "4                    NaN                      NaN                     NaN   \n",
       "...                  ...                      ...                     ...   \n",
       "12509           205275.0                  29304.0                 11087.0   \n",
       "12510           218516.0                  29839.0                 13776.0   \n",
       "12511           231632.0                  30496.0                 13773.0   \n",
       "12512           252932.0                  32949.0                 23753.0   \n",
       "12513           269732.0                  34969.0                 18820.0   \n",
       "\n",
       "       daily_vaccinations  total_vaccinations_per_hundred  \\\n",
       "0                     NaN                            0.00   \n",
       "1                  1367.0                             NaN   \n",
       "2                  1367.0                             NaN   \n",
       "3                  1367.0                             NaN   \n",
       "4                  1367.0                             NaN   \n",
       "...                   ...                             ...   \n",
       "12509             12724.0                            1.58   \n",
       "12510             12246.0                            1.67   \n",
       "12511             11816.0                            1.76   \n",
       "12512             13172.0                            1.92   \n",
       "12513             14071.0                            2.05   \n",
       "\n",
       "       people_vaccinated_per_hundred  people_fully_vaccinated_per_hundred  \\\n",
       "0                               0.00                                  NaN   \n",
       "1                                NaN                                  NaN   \n",
       "2                                NaN                                  NaN   \n",
       "3                                NaN                                  NaN   \n",
       "4                                NaN                                  NaN   \n",
       "...                              ...                                  ...   \n",
       "12509                           1.38                                 0.20   \n",
       "12510                           1.47                                 0.20   \n",
       "12511                           1.56                                 0.21   \n",
       "12512                           1.70                                 0.22   \n",
       "12513                           1.81                                 0.24   \n",
       "\n",
       "       daily_vaccinations_per_million            vaccines  \\\n",
       "0                                 NaN  Oxford/AstraZeneca   \n",
       "1                                35.0  Oxford/AstraZeneca   \n",
       "2                                35.0  Oxford/AstraZeneca   \n",
       "3                                35.0  Oxford/AstraZeneca   \n",
       "4                                35.0  Oxford/AstraZeneca   \n",
       "...                               ...                 ...   \n",
       "12509                           856.0   Sinopharm/Beijing   \n",
       "12510                           824.0   Sinopharm/Beijing   \n",
       "12511                           795.0   Sinopharm/Beijing   \n",
       "12512                           886.0   Sinopharm/Beijing   \n",
       "12513                           947.0   Sinopharm/Beijing   \n",
       "\n",
       "                     source_name  \\\n",
       "0      Government of Afghanistan   \n",
       "1      Government of Afghanistan   \n",
       "2      Government of Afghanistan   \n",
       "3      Government of Afghanistan   \n",
       "4      Government of Afghanistan   \n",
       "...                          ...   \n",
       "12509         Ministry of Health   \n",
       "12510         Ministry of Health   \n",
       "12511         Ministry of Health   \n",
       "12512         Ministry of Health   \n",
       "12513         Ministry of Health   \n",
       "\n",
       "                                          source_website  \n",
       "0      http://www.xinhuanet.com/english/asiapacific/2...  \n",
       "1      http://www.xinhuanet.com/english/asiapacific/2...  \n",
       "2      http://www.xinhuanet.com/english/asiapacific/2...  \n",
       "3      http://www.xinhuanet.com/english/asiapacific/2...  \n",
       "4      http://www.xinhuanet.com/english/asiapacific/2...  \n",
       "...                                                  ...  \n",
       "12509  https://twitter.com/MoHCCZim/status/1383174384...  \n",
       "12510  https://twitter.com/MoHCCZim/status/1383174384...  \n",
       "12511  https://twitter.com/MoHCCZim/status/1383174384...  \n",
       "12512  https://twitter.com/MoHCCZim/status/1383174384...  \n",
       "12513  https://twitter.com/MoHCCZim/status/1383174384...  \n",
       "\n",
       "[12514 rows x 15 columns]"
      ]
     },
     "execution_count": 34,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "file = 'country_vaccinations.csv'\n",
    "dfVaccines = pd.read_csv(file)\n",
    "dfVaccines"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 35,
   "metadata": {},
   "outputs": [],
   "source": [
    "def prepareVaccinesData(country, countryPop):\n",
    "    \n",
    "    data = dfVaccines[dfVaccines['country'] == country]\n",
    "    data = data[['country', 'date', 'total_vaccinations', 'people_vaccinated', 'people_fully_vaccinated']]\n",
    "    data['ppl_vac_percent'] = data['people_vaccinated'].apply(lambda x: (x / countryPop) * 100) # Share of people who received at least one dose of vaccine\n",
    "    data['Date'] = pd.to_datetime(data.date)\n",
    "    data.drop(columns = ['date'], inplace = True)\n",
    "    \n",
    "    return data"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 36,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>country</th>\n",
       "      <th>total_vaccinations</th>\n",
       "      <th>people_vaccinated</th>\n",
       "      <th>people_fully_vaccinated</th>\n",
       "      <th>ppl_vac_percent</th>\n",
       "      <th>Date</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>6676</th>\n",
       "      <td>Lithuania</td>\n",
       "      <td>2420.0</td>\n",
       "      <td>2420.0</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.088896</td>\n",
       "      <td>2020-12-27</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6677</th>\n",
       "      <td>Lithuania</td>\n",
       "      <td>5207.0</td>\n",
       "      <td>5207.0</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.191273</td>\n",
       "      <td>2020-12-28</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6678</th>\n",
       "      <td>Lithuania</td>\n",
       "      <td>8194.0</td>\n",
       "      <td>8194.0</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.300997</td>\n",
       "      <td>2020-12-29</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6679</th>\n",
       "      <td>Lithuania</td>\n",
       "      <td>9726.0</td>\n",
       "      <td>9726.0</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.357273</td>\n",
       "      <td>2020-12-30</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6680</th>\n",
       "      <td>Lithuania</td>\n",
       "      <td>9751.0</td>\n",
       "      <td>9751.0</td>\n",
       "      <td>NaN</td>\n",
       "      <td>0.358191</td>\n",
       "      <td>2020-12-31</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6784</th>\n",
       "      <td>Lithuania</td>\n",
       "      <td>729457.0</td>\n",
       "      <td>529676.0</td>\n",
       "      <td>199781.0</td>\n",
       "      <td>19.457008</td>\n",
       "      <td>2021-04-14</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6785</th>\n",
       "      <td>Lithuania</td>\n",
       "      <td>751976.0</td>\n",
       "      <td>547173.0</td>\n",
       "      <td>204803.0</td>\n",
       "      <td>20.099740</td>\n",
       "      <td>2021-04-15</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6786</th>\n",
       "      <td>Lithuania</td>\n",
       "      <td>771584.0</td>\n",
       "      <td>562166.0</td>\n",
       "      <td>209418.0</td>\n",
       "      <td>20.650489</td>\n",
       "      <td>2021-04-16</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6787</th>\n",
       "      <td>Lithuania</td>\n",
       "      <td>773239.0</td>\n",
       "      <td>563640.0</td>\n",
       "      <td>209599.0</td>\n",
       "      <td>20.704635</td>\n",
       "      <td>2021-04-17</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6788</th>\n",
       "      <td>Lithuania</td>\n",
       "      <td>773324.0</td>\n",
       "      <td>563724.0</td>\n",
       "      <td>209600.0</td>\n",
       "      <td>20.707721</td>\n",
       "      <td>2021-04-18</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>113 rows × 6 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "        country  total_vaccinations  people_vaccinated  \\\n",
       "6676  Lithuania              2420.0             2420.0   \n",
       "6677  Lithuania              5207.0             5207.0   \n",
       "6678  Lithuania              8194.0             8194.0   \n",
       "6679  Lithuania              9726.0             9726.0   \n",
       "6680  Lithuania              9751.0             9751.0   \n",
       "...         ...                 ...                ...   \n",
       "6784  Lithuania            729457.0           529676.0   \n",
       "6785  Lithuania            751976.0           547173.0   \n",
       "6786  Lithuania            771584.0           562166.0   \n",
       "6787  Lithuania            773239.0           563640.0   \n",
       "6788  Lithuania            773324.0           563724.0   \n",
       "\n",
       "      people_fully_vaccinated  ppl_vac_percent       Date  \n",
       "6676                      NaN         0.088896 2020-12-27  \n",
       "6677                      NaN         0.191273 2020-12-28  \n",
       "6678                      NaN         0.300997 2020-12-29  \n",
       "6679                      NaN         0.357273 2020-12-30  \n",
       "6680                      NaN         0.358191 2020-12-31  \n",
       "...                       ...              ...        ...  \n",
       "6784                 199781.0        19.457008 2021-04-14  \n",
       "6785                 204803.0        20.099740 2021-04-15  \n",
       "6786                 209418.0        20.650489 2021-04-16  \n",
       "6787                 209599.0        20.704635 2021-04-17  \n",
       "6788                 209600.0        20.707721 2021-04-18  \n",
       "\n",
       "[113 rows x 6 columns]"
      ]
     },
     "execution_count": 36,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "vaccLT = prepareVaccinesData('Lithuania', 2722289)\n",
    "vaccLT"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Merge vaccines and new cases data"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 78,
   "metadata": {},
   "outputs": [],
   "source": [
    "def mergeData(country, countryPop):\n",
    "    cases = prepareCasesData(country, countryPop)\n",
    "    vacc = prepareVaccinesData(country, countryPop)\n",
    "    \n",
    "    data = cases.merge(vacc, how = 'inner', on = 'Date')\n",
    "    \n",
    "    data.fillna(method='ffill', inplace=True)    \n",
    "    return data"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 79,
   "metadata": {},
   "outputs": [],
   "source": [
    "LT = mergeData('Lithuania', 2722289)\n",
    "df = LT[['Date', 'MA7', 'ppl_vac_percent']]"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### In order to have smooth animated plot, I create extra 10 rows for every row. In this way, plot doesn't jump from one point to another and gives a smoother look"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 80,
   "metadata": {},
   "outputs": [],
   "source": [
    "rows = 10\n",
    "df1 = pd.concat([df.iloc[0:1], pd.DataFrame([[np.nan]*df.shape[1]]*rows, columns=df.columns), df.iloc[1:2]], ignore_index=True)\n",
    "df1['MA7'] = df1['MA7'].interpolate()\n",
    "df1['ppl_vac_percent'] = df1['ppl_vac_percent'].interpolate()\n",
    "df1['Date'] = df1['Date'].ffill()\n",
    "\n",
    "for i in range(1, len(df)):\n",
    "    df1 = pd.concat([df1, pd.DataFrame([[np.nan]*df.shape[1]]*rows, columns=df.columns), df.iloc[i+1:i+2]], ignore_index=True)\n",
    "    df1['MA7'] = df1['MA7'].interpolate()\n",
    "    df1['ppl_vac_percent'] = df1['ppl_vac_percent'].interpolate()\n",
    "    df1['Date'] = df1['Date'].ffill()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 81,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>Date</th>\n",
       "      <th>MA7</th>\n",
       "      <th>ppl_vac_percent</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>2020-12-27</td>\n",
       "      <td>2669.857143</td>\n",
       "      <td>0.088896</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>2020-12-27</td>\n",
       "      <td>2665.428571</td>\n",
       "      <td>0.098203</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>2020-12-27</td>\n",
       "      <td>2661.000000</td>\n",
       "      <td>0.107510</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>2020-12-27</td>\n",
       "      <td>2656.571429</td>\n",
       "      <td>0.116817</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>2020-12-27</td>\n",
       "      <td>2652.142857</td>\n",
       "      <td>0.126124</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1228</th>\n",
       "      <td>2021-04-17</td>\n",
       "      <td>1068.701299</td>\n",
       "      <td>20.706599</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1229</th>\n",
       "      <td>2021-04-17</td>\n",
       "      <td>1067.883117</td>\n",
       "      <td>20.706879</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1230</th>\n",
       "      <td>2021-04-17</td>\n",
       "      <td>1067.064935</td>\n",
       "      <td>20.707160</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1231</th>\n",
       "      <td>2021-04-17</td>\n",
       "      <td>1066.246753</td>\n",
       "      <td>20.707440</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1232</th>\n",
       "      <td>2021-04-18</td>\n",
       "      <td>1065.428571</td>\n",
       "      <td>20.707721</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>1233 rows × 3 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "           Date          MA7  ppl_vac_percent\n",
       "0    2020-12-27  2669.857143         0.088896\n",
       "1    2020-12-27  2665.428571         0.098203\n",
       "2    2020-12-27  2661.000000         0.107510\n",
       "3    2020-12-27  2656.571429         0.116817\n",
       "4    2020-12-27  2652.142857         0.126124\n",
       "...         ...          ...              ...\n",
       "1228 2021-04-17  1068.701299        20.706599\n",
       "1229 2021-04-17  1067.883117        20.706879\n",
       "1230 2021-04-17  1067.064935        20.707160\n",
       "1231 2021-04-17  1066.246753        20.707440\n",
       "1232 2021-04-18  1065.428571        20.707721\n",
       "\n",
       "[1233 rows x 3 columns]"
      ]
     },
     "execution_count": 81,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "#Another way:\n",
    "rows = 10\n",
    "(df.assign(enum=df.groupby('Date').cumcount())\n",
    "   .set_index(['enum', 'Date'])\n",
    "   .unstack('Date')\n",
    "   .reindex(range(rows + 1))\n",
    "   .stack('Date', dropna=False)\n",
    "   .sort_index(level='Date')\n",
    "   .reset_index('Date')\n",
    "   .interpolate()\n",
    "   .reset_index(drop=True)\n",
    ").iloc[:-rows]"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Plotting"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 83,
   "metadata": {},
   "outputs": [],
   "source": [
    "%matplotlib notebook"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "application/javascript": [
       "/* Put everything inside the global mpl namespace */\n",
       "/* global mpl */\n",
       "window.mpl = {};\n",
       "\n",
       "mpl.get_websocket_type = function () {\n",
       "    if (typeof WebSocket !== 'undefined') {\n",
       "        return WebSocket;\n",
       "    } else if (typeof MozWebSocket !== 'undefined') {\n",
       "        return MozWebSocket;\n",
       "    } else {\n",
       "        alert(\n",
       "            'Your browser does not have WebSocket support. ' +\n",
       "                'Please try Chrome, Safari or Firefox ≥ 6. ' +\n",
       "                'Firefox 4 and 5 are also supported but you ' +\n",
       "                'have to enable WebSockets in about:config.'\n",
       "        );\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure = function (figure_id, websocket, ondownload, parent_element) {\n",
       "    this.id = figure_id;\n",
       "\n",
       "    this.ws = websocket;\n",
       "\n",
       "    this.supports_binary = this.ws.binaryType !== undefined;\n",
       "\n",
       "    if (!this.supports_binary) {\n",
       "        var warnings = document.getElementById('mpl-warnings');\n",
       "        if (warnings) {\n",
       "            warnings.style.display = 'block';\n",
       "            warnings.textContent =\n",
       "                'This browser does not support binary websocket messages. ' +\n",
       "                'Performance may be slow.';\n",
       "        }\n",
       "    }\n",
       "\n",
       "    this.imageObj = new Image();\n",
       "\n",
       "    this.context = undefined;\n",
       "    this.message = undefined;\n",
       "    this.canvas = undefined;\n",
       "    this.rubberband_canvas = undefined;\n",
       "    this.rubberband_context = undefined;\n",
       "    this.format_dropdown = undefined;\n",
       "\n",
       "    this.image_mode = 'full';\n",
       "\n",
       "    this.root = document.createElement('div');\n",
       "    this.root.setAttribute('style', 'display: inline-block');\n",
       "    this._root_extra_style(this.root);\n",
       "\n",
       "    parent_element.appendChild(this.root);\n",
       "\n",
       "    this._init_header(this);\n",
       "    this._init_canvas(this);\n",
       "    this._init_toolbar(this);\n",
       "\n",
       "    var fig = this;\n",
       "\n",
       "    this.waiting = false;\n",
       "\n",
       "    this.ws.onopen = function () {\n",
       "        fig.send_message('supports_binary', { value: fig.supports_binary });\n",
       "        fig.send_message('send_image_mode', {});\n",
       "        if (fig.ratio !== 1) {\n",
       "            fig.send_message('set_dpi_ratio', { dpi_ratio: fig.ratio });\n",
       "        }\n",
       "        fig.send_message('refresh', {});\n",
       "    };\n",
       "\n",
       "    this.imageObj.onload = function () {\n",
       "        if (fig.image_mode === 'full') {\n",
       "            // Full images could contain transparency (where diff images\n",
       "            // almost always do), so we need to clear the canvas so that\n",
       "            // there is no ghosting.\n",
       "            fig.context.clearRect(0, 0, fig.canvas.width, fig.canvas.height);\n",
       "        }\n",
       "        fig.context.drawImage(fig.imageObj, 0, 0);\n",
       "    };\n",
       "\n",
       "    this.imageObj.onunload = function () {\n",
       "        fig.ws.close();\n",
       "    };\n",
       "\n",
       "    this.ws.onmessage = this._make_on_message_function(this);\n",
       "\n",
       "    this.ondownload = ondownload;\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._init_header = function () {\n",
       "    var titlebar = document.createElement('div');\n",
       "    titlebar.classList =\n",
       "        'ui-dialog-titlebar ui-widget-header ui-corner-all ui-helper-clearfix';\n",
       "    var titletext = document.createElement('div');\n",
       "    titletext.classList = 'ui-dialog-title';\n",
       "    titletext.setAttribute(\n",
       "        'style',\n",
       "        'width: 100%; text-align: center; padding: 3px;'\n",
       "    );\n",
       "    titlebar.appendChild(titletext);\n",
       "    this.root.appendChild(titlebar);\n",
       "    this.header = titletext;\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._canvas_extra_style = function (_canvas_div) {};\n",
       "\n",
       "mpl.figure.prototype._root_extra_style = function (_canvas_div) {};\n",
       "\n",
       "mpl.figure.prototype._init_canvas = function () {\n",
       "    var fig = this;\n",
       "\n",
       "    var canvas_div = (this.canvas_div = document.createElement('div'));\n",
       "    canvas_div.setAttribute(\n",
       "        'style',\n",
       "        'border: 1px solid #ddd;' +\n",
       "            'box-sizing: content-box;' +\n",
       "            'clear: both;' +\n",
       "            'min-height: 1px;' +\n",
       "            'min-width: 1px;' +\n",
       "            'outline: 0;' +\n",
       "            'overflow: hidden;' +\n",
       "            'position: relative;' +\n",
       "            'resize: both;'\n",
       "    );\n",
       "\n",
       "    function on_keyboard_event_closure(name) {\n",
       "        return function (event) {\n",
       "            return fig.key_event(event, name);\n",
       "        };\n",
       "    }\n",
       "\n",
       "    canvas_div.addEventListener(\n",
       "        'keydown',\n",
       "        on_keyboard_event_closure('key_press')\n",
       "    );\n",
       "    canvas_div.addEventListener(\n",
       "        'keyup',\n",
       "        on_keyboard_event_closure('key_release')\n",
       "    );\n",
       "\n",
       "    this._canvas_extra_style(canvas_div);\n",
       "    this.root.appendChild(canvas_div);\n",
       "\n",
       "    var canvas = (this.canvas = document.createElement('canvas'));\n",
       "    canvas.classList.add('mpl-canvas');\n",
       "    canvas.setAttribute('style', 'box-sizing: content-box;');\n",
       "\n",
       "    this.context = canvas.getContext('2d');\n",
       "\n",
       "    var backingStore =\n",
       "        this.context.backingStorePixelRatio ||\n",
       "        this.context.webkitBackingStorePixelRatio ||\n",
       "        this.context.mozBackingStorePixelRatio ||\n",
       "        this.context.msBackingStorePixelRatio ||\n",
       "        this.context.oBackingStorePixelRatio ||\n",
       "        this.context.backingStorePixelRatio ||\n",
       "        1;\n",
       "\n",
       "    this.ratio = (window.devicePixelRatio || 1) / backingStore;\n",
       "    if (this.ratio !== 1) {\n",
       "        fig.send_message('set_dpi_ratio', { dpi_ratio: this.ratio });\n",
       "    }\n",
       "\n",
       "    var rubberband_canvas = (this.rubberband_canvas = document.createElement(\n",
       "        'canvas'\n",
       "    ));\n",
       "    rubberband_canvas.setAttribute(\n",
       "        'style',\n",
       "        'box-sizing: content-box; position: absolute; left: 0; top: 0; z-index: 1;'\n",
       "    );\n",
       "\n",
       "    // Apply a ponyfill if ResizeObserver is not implemented by browser.\n",
       "    if (this.ResizeObserver === undefined) {\n",
       "        if (window.ResizeObserver !== undefined) {\n",
       "            this.ResizeObserver = window.ResizeObserver;\n",
       "        } else {\n",
       "            var obs = _JSXTOOLS_RESIZE_OBSERVER({});\n",
       "            this.ResizeObserver = obs.ResizeObserver;\n",
       "        }\n",
       "    }\n",
       "\n",
       "    this.resizeObserverInstance = new this.ResizeObserver(function (entries) {\n",
       "        var nentries = entries.length;\n",
       "        for (var i = 0; i < nentries; i++) {\n",
       "            var entry = entries[i];\n",
       "            var width, height;\n",
       "            if (entry.contentBoxSize) {\n",
       "                if (entry.contentBoxSize instanceof Array) {\n",
       "                    // Chrome 84 implements new version of spec.\n",
       "                    width = entry.contentBoxSize[0].inlineSize;\n",
       "                    height = entry.contentBoxSize[0].blockSize;\n",
       "                } else {\n",
       "                    // Firefox implements old version of spec.\n",
       "                    width = entry.contentBoxSize.inlineSize;\n",
       "                    height = entry.contentBoxSize.blockSize;\n",
       "                }\n",
       "            } else {\n",
       "                // Chrome <84 implements even older version of spec.\n",
       "                width = entry.contentRect.width;\n",
       "                height = entry.contentRect.height;\n",
       "            }\n",
       "\n",
       "            // Keep the size of the canvas and rubber band canvas in sync with\n",
       "            // the canvas container.\n",
       "            if (entry.devicePixelContentBoxSize) {\n",
       "                // Chrome 84 implements new version of spec.\n",
       "                canvas.setAttribute(\n",
       "                    'width',\n",
       "                    entry.devicePixelContentBoxSize[0].inlineSize\n",
       "                );\n",
       "                canvas.setAttribute(\n",
       "                    'height',\n",
       "                    entry.devicePixelContentBoxSize[0].blockSize\n",
       "                );\n",
       "            } else {\n",
       "                canvas.setAttribute('width', width * fig.ratio);\n",
       "                canvas.setAttribute('height', height * fig.ratio);\n",
       "            }\n",
       "            canvas.setAttribute(\n",
       "                'style',\n",
       "                'width: ' + width + 'px; height: ' + height + 'px;'\n",
       "            );\n",
       "\n",
       "            rubberband_canvas.setAttribute('width', width);\n",
       "            rubberband_canvas.setAttribute('height', height);\n",
       "\n",
       "            // And update the size in Python. We ignore the initial 0/0 size\n",
       "            // that occurs as the element is placed into the DOM, which should\n",
       "            // otherwise not happen due to the minimum size styling.\n",
       "            if (width != 0 && height != 0) {\n",
       "                fig.request_resize(width, height);\n",
       "            }\n",
       "        }\n",
       "    });\n",
       "    this.resizeObserverInstance.observe(canvas_div);\n",
       "\n",
       "    function on_mouse_event_closure(name) {\n",
       "        return function (event) {\n",
       "            return fig.mouse_event(event, name);\n",
       "        };\n",
       "    }\n",
       "\n",
       "    rubberband_canvas.addEventListener(\n",
       "        'mousedown',\n",
       "        on_mouse_event_closure('button_press')\n",
       "    );\n",
       "    rubberband_canvas.addEventListener(\n",
       "        'mouseup',\n",
       "        on_mouse_event_closure('button_release')\n",
       "    );\n",
       "    // Throttle sequential mouse events to 1 every 20ms.\n",
       "    rubberband_canvas.addEventListener(\n",
       "        'mousemove',\n",
       "        on_mouse_event_closure('motion_notify')\n",
       "    );\n",
       "\n",
       "    rubberband_canvas.addEventListener(\n",
       "        'mouseenter',\n",
       "        on_mouse_event_closure('figure_enter')\n",
       "    );\n",
       "    rubberband_canvas.addEventListener(\n",
       "        'mouseleave',\n",
       "        on_mouse_event_closure('figure_leave')\n",
       "    );\n",
       "\n",
       "    canvas_div.addEventListener('wheel', function (event) {\n",
       "        if (event.deltaY < 0) {\n",
       "            event.step = 1;\n",
       "        } else {\n",
       "            event.step = -1;\n",
       "        }\n",
       "        on_mouse_event_closure('scroll')(event);\n",
       "    });\n",
       "\n",
       "    canvas_div.appendChild(canvas);\n",
       "    canvas_div.appendChild(rubberband_canvas);\n",
       "\n",
       "    this.rubberband_context = rubberband_canvas.getContext('2d');\n",
       "    this.rubberband_context.strokeStyle = '#000000';\n",
       "\n",
       "    this._resize_canvas = function (width, height, forward) {\n",
       "        if (forward) {\n",
       "            canvas_div.style.width = width + 'px';\n",
       "            canvas_div.style.height = height + 'px';\n",
       "        }\n",
       "    };\n",
       "\n",
       "    // Disable right mouse context menu.\n",
       "    this.rubberband_canvas.addEventListener('contextmenu', function (_e) {\n",
       "        event.preventDefault();\n",
       "        return false;\n",
       "    });\n",
       "\n",
       "    function set_focus() {\n",
       "        canvas.focus();\n",
       "        canvas_div.focus();\n",
       "    }\n",
       "\n",
       "    window.setTimeout(set_focus, 100);\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._init_toolbar = function () {\n",
       "    var fig = this;\n",
       "\n",
       "    var toolbar = document.createElement('div');\n",
       "    toolbar.classList = 'mpl-toolbar';\n",
       "    this.root.appendChild(toolbar);\n",
       "\n",
       "    function on_click_closure(name) {\n",
       "        return function (_event) {\n",
       "            return fig.toolbar_button_onclick(name);\n",
       "        };\n",
       "    }\n",
       "\n",
       "    function on_mouseover_closure(tooltip) {\n",
       "        return function (event) {\n",
       "            if (!event.currentTarget.disabled) {\n",
       "                return fig.toolbar_button_onmouseover(tooltip);\n",
       "            }\n",
       "        };\n",
       "    }\n",
       "\n",
       "    fig.buttons = {};\n",
       "    var buttonGroup = document.createElement('div');\n",
       "    buttonGroup.classList = 'mpl-button-group';\n",
       "    for (var toolbar_ind in mpl.toolbar_items) {\n",
       "        var name = mpl.toolbar_items[toolbar_ind][0];\n",
       "        var tooltip = mpl.toolbar_items[toolbar_ind][1];\n",
       "        var image = mpl.toolbar_items[toolbar_ind][2];\n",
       "        var method_name = mpl.toolbar_items[toolbar_ind][3];\n",
       "\n",
       "        if (!name) {\n",
       "            /* Instead of a spacer, we start a new button group. */\n",
       "            if (buttonGroup.hasChildNodes()) {\n",
       "                toolbar.appendChild(buttonGroup);\n",
       "            }\n",
       "            buttonGroup = document.createElement('div');\n",
       "            buttonGroup.classList = 'mpl-button-group';\n",
       "            continue;\n",
       "        }\n",
       "\n",
       "        var button = (fig.buttons[name] = document.createElement('button'));\n",
       "        button.classList = 'mpl-widget';\n",
       "        button.setAttribute('role', 'button');\n",
       "        button.setAttribute('aria-disabled', 'false');\n",
       "        button.addEventListener('click', on_click_closure(method_name));\n",
       "        button.addEventListener('mouseover', on_mouseover_closure(tooltip));\n",
       "\n",
       "        var icon_img = document.createElement('img');\n",
       "        icon_img.src = '_images/' + image + '.png';\n",
       "        icon_img.srcset = '_images/' + image + '_large.png 2x';\n",
       "        icon_img.alt = tooltip;\n",
       "        button.appendChild(icon_img);\n",
       "\n",
       "        buttonGroup.appendChild(button);\n",
       "    }\n",
       "\n",
       "    if (buttonGroup.hasChildNodes()) {\n",
       "        toolbar.appendChild(buttonGroup);\n",
       "    }\n",
       "\n",
       "    var fmt_picker = document.createElement('select');\n",
       "    fmt_picker.classList = 'mpl-widget';\n",
       "    toolbar.appendChild(fmt_picker);\n",
       "    this.format_dropdown = fmt_picker;\n",
       "\n",
       "    for (var ind in mpl.extensions) {\n",
       "        var fmt = mpl.extensions[ind];\n",
       "        var option = document.createElement('option');\n",
       "        option.selected = fmt === mpl.default_extension;\n",
       "        option.innerHTML = fmt;\n",
       "        fmt_picker.appendChild(option);\n",
       "    }\n",
       "\n",
       "    var status_bar = document.createElement('span');\n",
       "    status_bar.classList = 'mpl-message';\n",
       "    toolbar.appendChild(status_bar);\n",
       "    this.message = status_bar;\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.request_resize = function (x_pixels, y_pixels) {\n",
       "    // Request matplotlib to resize the figure. Matplotlib will then trigger a resize in the client,\n",
       "    // which will in turn request a refresh of the image.\n",
       "    this.send_message('resize', { width: x_pixels, height: y_pixels });\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.send_message = function (type, properties) {\n",
       "    properties['type'] = type;\n",
       "    properties['figure_id'] = this.id;\n",
       "    this.ws.send(JSON.stringify(properties));\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.send_draw_message = function () {\n",
       "    if (!this.waiting) {\n",
       "        this.waiting = true;\n",
       "        this.ws.send(JSON.stringify({ type: 'draw', figure_id: this.id }));\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_save = function (fig, _msg) {\n",
       "    var format_dropdown = fig.format_dropdown;\n",
       "    var format = format_dropdown.options[format_dropdown.selectedIndex].value;\n",
       "    fig.ondownload(fig, format);\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_resize = function (fig, msg) {\n",
       "    var size = msg['size'];\n",
       "    if (size[0] !== fig.canvas.width || size[1] !== fig.canvas.height) {\n",
       "        fig._resize_canvas(size[0], size[1], msg['forward']);\n",
       "        fig.send_message('refresh', {});\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_rubberband = function (fig, msg) {\n",
       "    var x0 = msg['x0'] / fig.ratio;\n",
       "    var y0 = (fig.canvas.height - msg['y0']) / fig.ratio;\n",
       "    var x1 = msg['x1'] / fig.ratio;\n",
       "    var y1 = (fig.canvas.height - msg['y1']) / fig.ratio;\n",
       "    x0 = Math.floor(x0) + 0.5;\n",
       "    y0 = Math.floor(y0) + 0.5;\n",
       "    x1 = Math.floor(x1) + 0.5;\n",
       "    y1 = Math.floor(y1) + 0.5;\n",
       "    var min_x = Math.min(x0, x1);\n",
       "    var min_y = Math.min(y0, y1);\n",
       "    var width = Math.abs(x1 - x0);\n",
       "    var height = Math.abs(y1 - y0);\n",
       "\n",
       "    fig.rubberband_context.clearRect(\n",
       "        0,\n",
       "        0,\n",
       "        fig.canvas.width / fig.ratio,\n",
       "        fig.canvas.height / fig.ratio\n",
       "    );\n",
       "\n",
       "    fig.rubberband_context.strokeRect(min_x, min_y, width, height);\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_figure_label = function (fig, msg) {\n",
       "    // Updates the figure title.\n",
       "    fig.header.textContent = msg['label'];\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_cursor = function (fig, msg) {\n",
       "    var cursor = msg['cursor'];\n",
       "    switch (cursor) {\n",
       "        case 0:\n",
       "            cursor = 'pointer';\n",
       "            break;\n",
       "        case 1:\n",
       "            cursor = 'default';\n",
       "            break;\n",
       "        case 2:\n",
       "            cursor = 'crosshair';\n",
       "            break;\n",
       "        case 3:\n",
       "            cursor = 'move';\n",
       "            break;\n",
       "    }\n",
       "    fig.rubberband_canvas.style.cursor = cursor;\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_message = function (fig, msg) {\n",
       "    fig.message.textContent = msg['message'];\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_draw = function (fig, _msg) {\n",
       "    // Request the server to send over a new figure.\n",
       "    fig.send_draw_message();\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_image_mode = function (fig, msg) {\n",
       "    fig.image_mode = msg['mode'];\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_history_buttons = function (fig, msg) {\n",
       "    for (var key in msg) {\n",
       "        if (!(key in fig.buttons)) {\n",
       "            continue;\n",
       "        }\n",
       "        fig.buttons[key].disabled = !msg[key];\n",
       "        fig.buttons[key].setAttribute('aria-disabled', !msg[key]);\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_navigate_mode = function (fig, msg) {\n",
       "    if (msg['mode'] === 'PAN') {\n",
       "        fig.buttons['Pan'].classList.add('active');\n",
       "        fig.buttons['Zoom'].classList.remove('active');\n",
       "    } else if (msg['mode'] === 'ZOOM') {\n",
       "        fig.buttons['Pan'].classList.remove('active');\n",
       "        fig.buttons['Zoom'].classList.add('active');\n",
       "    } else {\n",
       "        fig.buttons['Pan'].classList.remove('active');\n",
       "        fig.buttons['Zoom'].classList.remove('active');\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.updated_canvas_event = function () {\n",
       "    // Called whenever the canvas gets updated.\n",
       "    this.send_message('ack', {});\n",
       "};\n",
       "\n",
       "// A function to construct a web socket function for onmessage handling.\n",
       "// Called in the figure constructor.\n",
       "mpl.figure.prototype._make_on_message_function = function (fig) {\n",
       "    return function socket_on_message(evt) {\n",
       "        if (evt.data instanceof Blob) {\n",
       "            /* FIXME: We get \"Resource interpreted as Image but\n",
       "             * transferred with MIME type text/plain:\" errors on\n",
       "             * Chrome.  But how to set the MIME type?  It doesn't seem\n",
       "             * to be part of the websocket stream */\n",
       "            evt.data.type = 'image/png';\n",
       "\n",
       "            /* Free the memory for the previous frames */\n",
       "            if (fig.imageObj.src) {\n",
       "                (window.URL || window.webkitURL).revokeObjectURL(\n",
       "                    fig.imageObj.src\n",
       "                );\n",
       "            }\n",
       "\n",
       "            fig.imageObj.src = (window.URL || window.webkitURL).createObjectURL(\n",
       "                evt.data\n",
       "            );\n",
       "            fig.updated_canvas_event();\n",
       "            fig.waiting = false;\n",
       "            return;\n",
       "        } else if (\n",
       "            typeof evt.data === 'string' &&\n",
       "            evt.data.slice(0, 21) === 'data:image/png;base64'\n",
       "        ) {\n",
       "            fig.imageObj.src = evt.data;\n",
       "            fig.updated_canvas_event();\n",
       "            fig.waiting = false;\n",
       "            return;\n",
       "        }\n",
       "\n",
       "        var msg = JSON.parse(evt.data);\n",
       "        var msg_type = msg['type'];\n",
       "\n",
       "        // Call the  \"handle_{type}\" callback, which takes\n",
       "        // the figure and JSON message as its only arguments.\n",
       "        try {\n",
       "            var callback = fig['handle_' + msg_type];\n",
       "        } catch (e) {\n",
       "            console.log(\n",
       "                \"No handler for the '\" + msg_type + \"' message type: \",\n",
       "                msg\n",
       "            );\n",
       "            return;\n",
       "        }\n",
       "\n",
       "        if (callback) {\n",
       "            try {\n",
       "                // console.log(\"Handling '\" + msg_type + \"' message: \", msg);\n",
       "                callback(fig, msg);\n",
       "            } catch (e) {\n",
       "                console.log(\n",
       "                    \"Exception inside the 'handler_\" + msg_type + \"' callback:\",\n",
       "                    e,\n",
       "                    e.stack,\n",
       "                    msg\n",
       "                );\n",
       "            }\n",
       "        }\n",
       "    };\n",
       "};\n",
       "\n",
       "// from http://stackoverflow.com/questions/1114465/getting-mouse-location-in-canvas\n",
       "mpl.findpos = function (e) {\n",
       "    //this section is from http://www.quirksmode.org/js/events_properties.html\n",
       "    var targ;\n",
       "    if (!e) {\n",
       "        e = window.event;\n",
       "    }\n",
       "    if (e.target) {\n",
       "        targ = e.target;\n",
       "    } else if (e.srcElement) {\n",
       "        targ = e.srcElement;\n",
       "    }\n",
       "    if (targ.nodeType === 3) {\n",
       "        // defeat Safari bug\n",
       "        targ = targ.parentNode;\n",
       "    }\n",
       "\n",
       "    // pageX,Y are the mouse positions relative to the document\n",
       "    var boundingRect = targ.getBoundingClientRect();\n",
       "    var x = e.pageX - (boundingRect.left + document.body.scrollLeft);\n",
       "    var y = e.pageY - (boundingRect.top + document.body.scrollTop);\n",
       "\n",
       "    return { x: x, y: y };\n",
       "};\n",
       "\n",
       "/*\n",
       " * return a copy of an object with only non-object keys\n",
       " * we need this to avoid circular references\n",
       " * http://stackoverflow.com/a/24161582/3208463\n",
       " */\n",
       "function simpleKeys(original) {\n",
       "    return Object.keys(original).reduce(function (obj, key) {\n",
       "        if (typeof original[key] !== 'object') {\n",
       "            obj[key] = original[key];\n",
       "        }\n",
       "        return obj;\n",
       "    }, {});\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.mouse_event = function (event, name) {\n",
       "    var canvas_pos = mpl.findpos(event);\n",
       "\n",
       "    if (name === 'button_press') {\n",
       "        this.canvas.focus();\n",
       "        this.canvas_div.focus();\n",
       "    }\n",
       "\n",
       "    var x = canvas_pos.x * this.ratio;\n",
       "    var y = canvas_pos.y * this.ratio;\n",
       "\n",
       "    this.send_message(name, {\n",
       "        x: x,\n",
       "        y: y,\n",
       "        button: event.button,\n",
       "        step: event.step,\n",
       "        guiEvent: simpleKeys(event),\n",
       "    });\n",
       "\n",
       "    /* This prevents the web browser from automatically changing to\n",
       "     * the text insertion cursor when the button is pressed.  We want\n",
       "     * to control all of the cursor setting manually through the\n",
       "     * 'cursor' event from matplotlib */\n",
       "    event.preventDefault();\n",
       "    return false;\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._key_event_extra = function (_event, _name) {\n",
       "    // Handle any extra behaviour associated with a key event\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.key_event = function (event, name) {\n",
       "    // Prevent repeat events\n",
       "    if (name === 'key_press') {\n",
       "        if (event.which === this._key) {\n",
       "            return;\n",
       "        } else {\n",
       "            this._key = event.which;\n",
       "        }\n",
       "    }\n",
       "    if (name === 'key_release') {\n",
       "        this._key = null;\n",
       "    }\n",
       "\n",
       "    var value = '';\n",
       "    if (event.ctrlKey && event.which !== 17) {\n",
       "        value += 'ctrl+';\n",
       "    }\n",
       "    if (event.altKey && event.which !== 18) {\n",
       "        value += 'alt+';\n",
       "    }\n",
       "    if (event.shiftKey && event.which !== 16) {\n",
       "        value += 'shift+';\n",
       "    }\n",
       "\n",
       "    value += 'k';\n",
       "    value += event.which.toString();\n",
       "\n",
       "    this._key_event_extra(event, name);\n",
       "\n",
       "    this.send_message(name, { key: value, guiEvent: simpleKeys(event) });\n",
       "    return false;\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.toolbar_button_onclick = function (name) {\n",
       "    if (name === 'download') {\n",
       "        this.handle_save(this, null);\n",
       "    } else {\n",
       "        this.send_message('toolbar_button', { name: name });\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.toolbar_button_onmouseover = function (tooltip) {\n",
       "    this.message.textContent = tooltip;\n",
       "};\n",
       "\n",
       "///////////////// REMAINING CONTENT GENERATED BY embed_js.py /////////////////\n",
       "// prettier-ignore\n",
       "var _JSXTOOLS_RESIZE_OBSERVER=function(A){var t,i=new WeakMap,n=new WeakMap,a=new WeakMap,r=new WeakMap,o=new Set;function s(e){if(!(this instanceof s))throw new TypeError(\"Constructor requires 'new' operator\");i.set(this,e)}function h(){throw new TypeError(\"Function is not a constructor\")}function c(e,t,i,n){e=0 in arguments?Number(arguments[0]):0,t=1 in arguments?Number(arguments[1]):0,i=2 in arguments?Number(arguments[2]):0,n=3 in arguments?Number(arguments[3]):0,this.right=(this.x=this.left=e)+(this.width=i),this.bottom=(this.y=this.top=t)+(this.height=n),Object.freeze(this)}function d(){t=requestAnimationFrame(d);var s=new WeakMap,p=new Set;o.forEach((function(t){r.get(t).forEach((function(i){var r=t instanceof window.SVGElement,o=a.get(t),d=r?0:parseFloat(o.paddingTop),f=r?0:parseFloat(o.paddingRight),l=r?0:parseFloat(o.paddingBottom),u=r?0:parseFloat(o.paddingLeft),g=r?0:parseFloat(o.borderTopWidth),m=r?0:parseFloat(o.borderRightWidth),w=r?0:parseFloat(o.borderBottomWidth),b=u+f,F=d+l,v=(r?0:parseFloat(o.borderLeftWidth))+m,W=g+w,y=r?0:t.offsetHeight-W-t.clientHeight,E=r?0:t.offsetWidth-v-t.clientWidth,R=b+v,z=F+W,M=r?t.width:parseFloat(o.width)-R-E,O=r?t.height:parseFloat(o.height)-z-y;if(n.has(t)){var k=n.get(t);if(k[0]===M&&k[1]===O)return}n.set(t,[M,O]);var S=Object.create(h.prototype);S.target=t,S.contentRect=new c(u,d,M,O),s.has(i)||(s.set(i,[]),p.add(i)),s.get(i).push(S)}))})),p.forEach((function(e){i.get(e).call(e,s.get(e),e)}))}return s.prototype.observe=function(i){if(i instanceof window.Element){r.has(i)||(r.set(i,new Set),o.add(i),a.set(i,window.getComputedStyle(i)));var n=r.get(i);n.has(this)||n.add(this),cancelAnimationFrame(t),t=requestAnimationFrame(d)}},s.prototype.unobserve=function(i){if(i instanceof window.Element&&r.has(i)){var n=r.get(i);n.has(this)&&(n.delete(this),n.size||(r.delete(i),o.delete(i))),n.size||r.delete(i),o.size||cancelAnimationFrame(t)}},A.DOMRectReadOnly=c,A.ResizeObserver=s,A.ResizeObserverEntry=h,A}; // eslint-disable-line\n",
       "mpl.toolbar_items = [[\"Home\", \"Reset original view\", \"fa fa-home icon-home\", \"home\"], [\"Back\", \"Back to previous view\", \"fa fa-arrow-left icon-arrow-left\", \"back\"], [\"Forward\", \"Forward to next view\", \"fa fa-arrow-right icon-arrow-right\", \"forward\"], [\"\", \"\", \"\", \"\"], [\"Pan\", \"Left button pans, Right button zooms\\nx/y fixes axis, CTRL fixes aspect\", \"fa fa-arrows icon-move\", \"pan\"], [\"Zoom\", \"Zoom to rectangle\\nx/y fixes axis, CTRL fixes aspect\", \"fa fa-square-o icon-check-empty\", \"zoom\"], [\"\", \"\", \"\", \"\"], [\"Download\", \"Download plot\", \"fa fa-floppy-o icon-save\", \"download\"]];\n",
       "\n",
       "mpl.extensions = [\"eps\", \"jpeg\", \"pdf\", \"png\", \"ps\", \"raw\", \"svg\", \"tif\"];\n",
       "\n",
       "mpl.default_extension = \"png\";/* global mpl */\n",
       "\n",
       "var comm_websocket_adapter = function (comm) {\n",
       "    // Create a \"websocket\"-like object which calls the given IPython comm\n",
       "    // object with the appropriate methods. Currently this is a non binary\n",
       "    // socket, so there is still some room for performance tuning.\n",
       "    var ws = {};\n",
       "\n",
       "    ws.close = function () {\n",
       "        comm.close();\n",
       "    };\n",
       "    ws.send = function (m) {\n",
       "        //console.log('sending', m);\n",
       "        comm.send(m);\n",
       "    };\n",
       "    // Register the callback with on_msg.\n",
       "    comm.on_msg(function (msg) {\n",
       "        //console.log('receiving', msg['content']['data'], msg);\n",
       "        // Pass the mpl event to the overridden (by mpl) onmessage function.\n",
       "        ws.onmessage(msg['content']['data']);\n",
       "    });\n",
       "    return ws;\n",
       "};\n",
       "\n",
       "mpl.mpl_figure_comm = function (comm, msg) {\n",
       "    // This is the function which gets called when the mpl process\n",
       "    // starts-up an IPython Comm through the \"matplotlib\" channel.\n",
       "\n",
       "    var id = msg.content.data.id;\n",
       "    // Get hold of the div created by the display call when the Comm\n",
       "    // socket was opened in Python.\n",
       "    var element = document.getElementById(id);\n",
       "    var ws_proxy = comm_websocket_adapter(comm);\n",
       "\n",
       "    function ondownload(figure, _format) {\n",
       "        window.open(figure.canvas.toDataURL());\n",
       "    }\n",
       "\n",
       "    var fig = new mpl.figure(id, ws_proxy, ondownload, element);\n",
       "\n",
       "    // Call onopen now - mpl needs it, as it is assuming we've passed it a real\n",
       "    // web socket which is closed, not our websocket->open comm proxy.\n",
       "    ws_proxy.onopen();\n",
       "\n",
       "    fig.parent_element = element;\n",
       "    fig.cell_info = mpl.find_output_cell(\"<div id='\" + id + \"'></div>\");\n",
       "    if (!fig.cell_info) {\n",
       "        console.error('Failed to find cell for figure', id, fig);\n",
       "        return;\n",
       "    }\n",
       "    fig.cell_info[0].output_area.element.on(\n",
       "        'cleared',\n",
       "        { fig: fig },\n",
       "        fig._remove_fig_handler\n",
       "    );\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_close = function (fig, msg) {\n",
       "    var width = fig.canvas.width / fig.ratio;\n",
       "    fig.cell_info[0].output_area.element.off(\n",
       "        'cleared',\n",
       "        fig._remove_fig_handler\n",
       "    );\n",
       "    fig.resizeObserverInstance.unobserve(fig.canvas_div);\n",
       "\n",
       "    // Update the output cell to use the data from the current canvas.\n",
       "    fig.push_to_output();\n",
       "    var dataURL = fig.canvas.toDataURL();\n",
       "    // Re-enable the keyboard manager in IPython - without this line, in FF,\n",
       "    // the notebook keyboard shortcuts fail.\n",
       "    IPython.keyboard_manager.enable();\n",
       "    fig.parent_element.innerHTML =\n",
       "        '<img src=\"' + dataURL + '\" width=\"' + width + '\">';\n",
       "    fig.close_ws(fig, msg);\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.close_ws = function (fig, msg) {\n",
       "    fig.send_message('closing', msg);\n",
       "    // fig.ws.close()\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.push_to_output = function (_remove_interactive) {\n",
       "    // Turn the data on the canvas into data in the output cell.\n",
       "    var width = this.canvas.width / this.ratio;\n",
       "    var dataURL = this.canvas.toDataURL();\n",
       "    this.cell_info[1]['text/html'] =\n",
       "        '<img src=\"' + dataURL + '\" width=\"' + width + '\">';\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.updated_canvas_event = function () {\n",
       "    // Tell IPython that the notebook contents must change.\n",
       "    IPython.notebook.set_dirty(true);\n",
       "    this.send_message('ack', {});\n",
       "    var fig = this;\n",
       "    // Wait a second, then push the new image to the DOM so\n",
       "    // that it is saved nicely (might be nice to debounce this).\n",
       "    setTimeout(function () {\n",
       "        fig.push_to_output();\n",
       "    }, 1000);\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._init_toolbar = function () {\n",
       "    var fig = this;\n",
       "\n",
       "    var toolbar = document.createElement('div');\n",
       "    toolbar.classList = 'btn-toolbar';\n",
       "    this.root.appendChild(toolbar);\n",
       "\n",
       "    function on_click_closure(name) {\n",
       "        return function (_event) {\n",
       "            return fig.toolbar_button_onclick(name);\n",
       "        };\n",
       "    }\n",
       "\n",
       "    function on_mouseover_closure(tooltip) {\n",
       "        return function (event) {\n",
       "            if (!event.currentTarget.disabled) {\n",
       "                return fig.toolbar_button_onmouseover(tooltip);\n",
       "            }\n",
       "        };\n",
       "    }\n",
       "\n",
       "    fig.buttons = {};\n",
       "    var buttonGroup = document.createElement('div');\n",
       "    buttonGroup.classList = 'btn-group';\n",
       "    var button;\n",
       "    for (var toolbar_ind in mpl.toolbar_items) {\n",
       "        var name = mpl.toolbar_items[toolbar_ind][0];\n",
       "        var tooltip = mpl.toolbar_items[toolbar_ind][1];\n",
       "        var image = mpl.toolbar_items[toolbar_ind][2];\n",
       "        var method_name = mpl.toolbar_items[toolbar_ind][3];\n",
       "\n",
       "        if (!name) {\n",
       "            /* Instead of a spacer, we start a new button group. */\n",
       "            if (buttonGroup.hasChildNodes()) {\n",
       "                toolbar.appendChild(buttonGroup);\n",
       "            }\n",
       "            buttonGroup = document.createElement('div');\n",
       "            buttonGroup.classList = 'btn-group';\n",
       "            continue;\n",
       "        }\n",
       "\n",
       "        button = fig.buttons[name] = document.createElement('button');\n",
       "        button.classList = 'btn btn-default';\n",
       "        button.href = '#';\n",
       "        button.title = name;\n",
       "        button.innerHTML = '<i class=\"fa ' + image + ' fa-lg\"></i>';\n",
       "        button.addEventListener('click', on_click_closure(method_name));\n",
       "        button.addEventListener('mouseover', on_mouseover_closure(tooltip));\n",
       "        buttonGroup.appendChild(button);\n",
       "    }\n",
       "\n",
       "    if (buttonGroup.hasChildNodes()) {\n",
       "        toolbar.appendChild(buttonGroup);\n",
       "    }\n",
       "\n",
       "    // Add the status bar.\n",
       "    var status_bar = document.createElement('span');\n",
       "    status_bar.classList = 'mpl-message pull-right';\n",
       "    toolbar.appendChild(status_bar);\n",
       "    this.message = status_bar;\n",
       "\n",
       "    // Add the close button to the window.\n",
       "    var buttongrp = document.createElement('div');\n",
       "    buttongrp.classList = 'btn-group inline pull-right';\n",
       "    button = document.createElement('button');\n",
       "    button.classList = 'btn btn-mini btn-primary';\n",
       "    button.href = '#';\n",
       "    button.title = 'Stop Interaction';\n",
       "    button.innerHTML = '<i class=\"fa fa-power-off icon-remove icon-large\"></i>';\n",
       "    button.addEventListener('click', function (_evt) {\n",
       "        fig.handle_close(fig, {});\n",
       "    });\n",
       "    button.addEventListener(\n",
       "        'mouseover',\n",
       "        on_mouseover_closure('Stop Interaction')\n",
       "    );\n",
       "    buttongrp.appendChild(button);\n",
       "    var titlebar = this.root.querySelector('.ui-dialog-titlebar');\n",
       "    titlebar.insertBefore(buttongrp, titlebar.firstChild);\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._remove_fig_handler = function (event) {\n",
       "    var fig = event.data.fig;\n",
       "    if (event.target !== this) {\n",
       "        // Ignore bubbled events from children.\n",
       "        return;\n",
       "    }\n",
       "    fig.close_ws(fig, {});\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._root_extra_style = function (el) {\n",
       "    el.style.boxSizing = 'content-box'; // override notebook setting of border-box.\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._canvas_extra_style = function (el) {\n",
       "    // this is important to make the div 'focusable\n",
       "    el.setAttribute('tabindex', 0);\n",
       "    // reach out to IPython and tell the keyboard manager to turn it's self\n",
       "    // off when our div gets focus\n",
       "\n",
       "    // location in version 3\n",
       "    if (IPython.notebook.keyboard_manager) {\n",
       "        IPython.notebook.keyboard_manager.register_events(el);\n",
       "    } else {\n",
       "        // location in version 2\n",
       "        IPython.keyboard_manager.register_events(el);\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._key_event_extra = function (event, _name) {\n",
       "    var manager = IPython.notebook.keyboard_manager;\n",
       "    if (!manager) {\n",
       "        manager = IPython.keyboard_manager;\n",
       "    }\n",
       "\n",
       "    // Check for shift+enter\n",
       "    if (event.shiftKey && event.which === 13) {\n",
       "        this.canvas_div.blur();\n",
       "        // select the cell after this one\n",
       "        var index = IPython.notebook.find_cell_index(this.cell_info[0]);\n",
       "        IPython.notebook.select(index + 1);\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_save = function (fig, _msg) {\n",
       "    fig.ondownload(fig, null);\n",
       "};\n",
       "\n",
       "mpl.find_output_cell = function (html_output) {\n",
       "    // Return the cell and output element which can be found *uniquely* in the notebook.\n",
       "    // Note - this is a bit hacky, but it is done because the \"notebook_saving.Notebook\"\n",
       "    // IPython event is triggered only after the cells have been serialised, which for\n",
       "    // our purposes (turning an active figure into a static one), is too late.\n",
       "    var cells = IPython.notebook.get_cells();\n",
       "    var ncells = cells.length;\n",
       "    for (var i = 0; i < ncells; i++) {\n",
       "        var cell = cells[i];\n",
       "        if (cell.cell_type === 'code') {\n",
       "            for (var j = 0; j < cell.output_area.outputs.length; j++) {\n",
       "                var data = cell.output_area.outputs[j];\n",
       "                if (data.data) {\n",
       "                    // IPython >= 3 moved mimebundle to data attribute of output\n",
       "                    data = data.data;\n",
       "                }\n",
       "                if (data['text/html'] === html_output) {\n",
       "                    return [cell, data, j];\n",
       "                }\n",
       "            }\n",
       "        }\n",
       "    }\n",
       "};\n",
       "\n",
       "// Register the function which deals with the matplotlib target/channel.\n",
       "// The kernel may be null if the page has been refreshed.\n",
       "if (IPython.notebook.kernel !== null) {\n",
       "    IPython.notebook.kernel.comm_manager.register_target(\n",
       "        'matplotlib',\n",
       "        mpl.mpl_figure_comm\n",
       "    );\n",
       "}\n"
      ],
      "text/plain": [
       "<IPython.core.display.Javascript object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/html": [
       "<img src=\"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAoAAAAHgCAYAAAA10dzkAAAgAElEQVR4XuxdB5hURda9hBUVd3R+0UHECOKsmNMKKhhQUVkxjmAgqMgKRlSCIskcMIIJXQVcFQyYc86uCQMmzGRUEAZEJP3fKayhp6fD67536Jo3p76PT6Rf1bt16oZTt8KrIyIrhIUIEAEiQASIABEgAkSg1iBQhwSw1ow1O0oEiAARIAJEgAgQAYcACSAVgQgQASJABIgAESACtQwBEsBaNuDsLhEgAkSACBABIkAESACpA0SACBABIkAEiAARqGUIkADWsgFnd4kAESACRIAIEAEiQAJIHSACRIAIEAEiQASIQC1DgASwlg14obs7ceJE2WGHHWTx4sXSpEkTmTNnTqFFSvv++vXrywknnCAdO3aUnXfeWTbYYANZsWKFzJw5Uz788EOZMGGCPPDAA7JkyZKC9OH777+XzTff3P358ccfI8tw1113Sbdu3dyf0aNHR65Xt25dOeKII2SXXXap+LP++uvL0qVL5W9/+1vWdg488EA566yzZPfdd5e///3vMmPGDHnmmWfkkksukWnTpmWtH+cHXn75Zdlnn33cn1dffTXOXQ2qb7BnlDp1EAqjlXzqRGs5t6dg9z169JCDDjpIWrRoIQ0aNHB+4KGHHpJLL71UFixYkFuDfLrWIUACWOuGvHAd3nXXXeW9996rEABk4MYbbyycQBnevNNOO8mDDz4oW265pSxfvlw+/vhj+fbbb93f4XhBCEEQf/jhB9lmm21k0aJFq70fq5sArrvuuvLbb79V6WcUAjhs2DC56KKLXN0PPvhAIPt2220nW2+9tcydO1f2228/weSgthYSwMKMfD5kLp861dG7P//80028vvrqK+efYJ/77ruvrLHGGu7/99hjD/njjz+q49VsMyYIkADGZCBrQjduueUW+fe//y1Tp06Vpk2byieffOKygaEVkL/XX39dGjZsKI8//riceeaZjugllkaNGsk555wj5557rpSUlMi8efNWezdAThEAQExBwqKWfDOAa6+9ttx2223y0UcfuQwosrcINNkI4MEHHyxPPfWULFu2TI499liXofBl0KBBMnToUPnuu+/kH//4hyCo1cayySabCPD96aefCjKZqI2Yo8/5kDlMWlBAvApZkCnu37+/vP322xVitGzZUt566y0pKiqS3r17y80331xIEfnuwBEgAQx8gOIi3lprreWW/Pws9bHHHnPLgLvttpu8//77wXQTWb0vv/xSmjVr5pZ4jzrqqIogkUpIZDU/++yzGjXTzpcAJvd/s802c8Q4GwHEMi+Wqe6++27p3r17FRgx/lhWPvnkk+U///lPMLpAQeKPQD4EMHRUMME68sgj5dZbb5XTTjstdHEpXwERIAEsIPi16dVdunRx+80+/fRT2X777WXUqFFyyimnZHRSictiyB717dvXLWtg39lJJ51UsX9tzTXXlF69eskxxxwjpaWlgv/HXphHH31Urrzyypz2GXo5sUcRBGfWrFk5DxPILrKGZWVlstVWW0m9evXckicI5fDhwysto2Jf3LPPPitffPGFW0pOVVAfWdPGjRu7jCkypyiZloCLi4sF2TXs2UM99AOkG/927bXX5rUHMF8C+Ouvv8r//d//uf2U//3vf6t0EWOEsX366aflkEMOyYq3BjNMOI4++mi3127TTTd1cmEJ+n//+5/ccMMN8uKLL6Z9P8by7LPPlv3331+QsQPxxbi88sorLtMyadKkSnXXW289OeOMM+Swww5zeoClOUyCQHhBwkGMfUm3BJxI1pHxwVL6AQccIBhfvPu+++6Tiy++OG3mFFsV+vTpI3vvvbfLVC9cuNBtw0BfgXdyga4gq4SsLfDBlgeM39dff+2eh/5mK5jkoc78+fOdrXqShXqw0fHjx7smMNaJMgAfjAX2mqJ/icuX+P/zzjvP7cfdYostXDYZMo0bN05uuummKhOwtm3bunHBH7wHfcK7YdOzZ892baCkI4CQYcSIEY5AwWehDeCdqU6iPaL9AQMGuAku/BHs+/rrr5exY8emhA96OHjwYDn88MPdOGGfMfzFkCFDXL2oe3axlxH+Ydttt3Xvv+KKK7INF3+vxQiQANbiwV+dXUfwatOmjVs2hUNr1aqVW6rAnrKNNtooZQbNB8WRI0e6pWNk5uDc4CyRTbr//vtdXQRSkEoEHSxNlpeXuz162KsHp4xgj6W1KMXPnkGWEGxyLQhUIBFYRsayMAIQDokgIOEQCZY6sd/NH9qAw8bfQShAbt99990qr+zQoYNbisbeOWQcfUlHADfccEO3hI2N4VimBY4IaMABGIGooG+5HgLJlwCCTCO4g1Qkkh7f3vnnny9XXXWVC8wIftmKBrPnn3/e7ZMCBlOmTHGECNleZCBR0u1L7dy5s8tO+skFSBQwxTI8dA/EDEvZvuDfnnzySbfVATr+5ptvOr0Eqdpxxx0d4YQcvmQjgLAZZE9BkKAjsIE999zTLRuDKCDjk1wwCQHZxwQCy/bffPONmwzgEA4ODGAyAPLoC7CHjm288cZOJ1EHJAyHtbC0iHag31EKZMR78Cdx3y+2EJx66qmuCciGLRS+AI+XXnrJ/QHJ9gVkCv8Ge4aOvPbaa27rA57HUidkbteuXaWJlSeA77zzjhsnTK5Qz5NSTCJQUhFAbP0ASQXpe+6559yEAWPnSzrS6O0RujBw4EDniyZPnuzkhr9DwQQC5DuxYExgr82bN3f2CZ8BmdEH9BfkERO5KPYKgg7CDwKJPba//PJLlOHiM7UUARLAWjrwq7PbyH5gto79XQgu3il9/vnnbt/XiSeeKPfcc08VkXxQxA/I8GEPYXKB49xrr73kjjvucOTSn3xDsMLsF1mD5ICSqe8IfAjSOLCA06m5FmRkOnXqJAg8hx56aEX2MTGogAxAZl8QhBEw0i3ZeFJ6+umnC8iwL+kIIIIXsh0IeP/6179c0ENB8MZePBBNlCgBJVP/oy4Bg2iBCKXbk+T3huJdwOn333/PCnu+mLVv394dNkGATCzABOQU2VsQjunTp1f8jMkE9llhewAIIsYgMasFfcGeUAR8FJAyBG38O7Le6DeIpi8gLcgMJWYbsxFA1IU+IkuErBwKSBn0bJ111nEEA3/3BQQH2TUQCmxjgJ34guwQ9ACTDpAM6AkKdB7kBSQNE67Egr5jAgdbilJwCvWCCy5wfy6//PKKKtivCvKJP8A4cQ9wujrAHuODjP5xxx1XoR/AHGMG8o7MMjLMvngCiP/HPlVsQUiVzU8mcyC7TzzxhJvAgfD37Nmzyv7abAQQfg6EHBMAX7p27eomrakmvN6+oQOYmHmyiUwqJn7I3kaxV0wSoJ8YcxDi2nyoKoqO8hl+C5g6sBoQQADAEgxO1YKY+AJydvXVV7sZb2I2xP/ugyICJRxacoFTRwBApgIBFctCiQWZIjh/zIQR9JKX6FJ1HeQDJACO//bbb88JHQRUkDK8F1keLB0lFgQXZGHQfuvWrSs2byOLhMCI4IBsADJmviDI4YoUBH3URwbIl1QEEEQL+/IgA/oMIpJYEHB9YFhdBBAkCEvrIEg+0+ZlwjIp+o6MFgoyusnkLNUg5ItZpgH1BCR5svHwww+7DAxOrCPAZivIvCHLA71ExtYTtkz1shFALBtDx5MLlp6xTJk8YfGkCeQP8icXZLVwhVGiTWLJE2QVfX3kkUeydTPj7z6bh34h440CYo0MOIgQSDK2SCDriCwXSqqsIbKcb7zxhiPQGHP/rH85yDkygLB9TEj8dUKJBBAECm2kKolkDn4CpA12nGkCmI0AIgsH35Zc/IQ3UR5MEmDHaBOEPvlgCf4Nqx7ICGayV3/KHraDLQLYl8xCBLIhwAxgNoT4uwoBZOKQAUJgT97zg6VK7KvBM8gSIjgkFh8UEeCQHUsuCMjYY4XsGYJ3quKDGpadsO8wW9EQwOOPP95lMpOXahPfieU67PO58MIL5bLLLqv4CSQYQQtLjVja9gWEAzN77HVCZjGxpCKAyJAgG5KOMKC+v4txdRFALL9BHhBf9B/BFSQVARf7t7BciiVilERCkG2s8sEMbYJsIjsLgoysqL/DEDoIwgHyhqU6FAReLOUjyxZ1EoHsGpa7kzO2mfqTjQBiL1jiErNvC/oPO0jMHmPfHYgSlm9BsFPdU4mJxc8//+yycMjKo/j9ryAh2D+G5c/EzGW28Uj8PXE/H/DGNUm4sw6TKtgJMr3+7/fee687HJZq3yDsBJlPZMlAWlMVEG1MuNAu2kLxBBBZP0yq0hVP5jCZBCFGZhKHkVLtVfVtZCOA6e5y9Jk+2DHsGcX7DCyTY7k8U//S2StIJHwnSDAmWCR/uWhq7X6WBLB2j3+19x4b4LF0A6KHGXpyNgSZBix7gMCByCUWHxTT7R3DUg0CeZSSTLjS1dEsAeMgAw40JGc6E991zTXXuH1PyNwg2+KLD74IughGvniyhqVLHBZJLKkIoJcBwQwZllTFk9DEgIIxADFNLlhax5J1qhJ1CRh1MYYIqsl7yBD0sawJog7dwB67qBdr54MZDh5dd911jtClKzh4gUNGKNi36bNOIIpRrttBphmkF1nrTIdKUul6Mnnwh0AgD/6eXPzSYuIJ6+T7NjPZB7D25BtZ4zFjxlQspaKvyFohewadhj3mUqCvWIqGPkOvQXpA4jAZBAEEafFYI+uITCX+IGvpi89wwm6wVzRVgWyoA93HigKKJ4DIKvotD6nqejIHHDC+iSQyXV+zEUBkKmGbySXV6ft+/fq5rSq52mti27Bb2HSmSV8u48Znaw8CJIC1Z6wL0lOQP5BA7PvDIY7kguwDloawdONPHfpnsl2O6zMt2N+EZcRMBUQTsmQrfrkvn0MgGgKI7BiWbxAY/VIW9iFh2RTkGdgk7jtDPywJIEgYskzJJVOWMBcCiHaR5QEBQMYPe8pALpDtxDIeCAL0A3tCo5ZcMcN7kGlBpgQTAuyvwuEgv+fQZ6gSyVQoBDDdOKQigMgkgfhgL1ninYvpcE2+mgdjgIkVsqH4AwxQYBMgalGWtPG8P9zjyRsyjrBzZOtQYLMgXdBt7KvE0nvy8ruWAKbbXuKx8DaFiQ4mByClOICSfO9nInbZCGC6L/NkIoDYt4s7MlMV75PS6YDPImbra1S74nO1BwESwNoz1qu9p1h6wfIvgn2UgqADUheVAPoThdhvE+V6iigy+ICKfXgITMl7jjK14R1xqr1uvp535qkykv5qHP8bliKxnyxVdjQdAbRYUoqCE57JlQCma9fvX8JyMPqbS8kFM78XNfn0qX8fTiKDtCQSQCwB4xANiHnUJWAcvkDG1nIJOBcCiAwblnYhN5aAkycOueCLZ7GHD0urWJ4HWQQ+UYqfwCCLjQwm7CJxf5y3X1zdBHKJU+v4g5Ozvvgl4OTMYOL70S7elWoJOBspSiRzyAxj6R8TLmRv0130bEkAcXAFV8PgVPg///nPlLD6/qXTAWQ4sWcZEyisQLAQgagIkABGRYrP5YwADn4g6OJ0or8GIVUjWALBUkiyk8+WAfRLH9mWeXIRHGQVJ5aRlYQ8yFhlCqDIKiGThf1WiYdAEJD8fX3+/SDEyDAkHwLxv/urcfB+7I9DtgR7tbA3DYdHkkuqDCBkQPYi3aZyZN+wZyrbpvIomFkQQJyIReBClgmyJR9aySZHLph5wpHqqhfs/cI4Yvku+cJqv2SeuDcwk1x+3yYCNw5uRMmYZdsDmAsBhGx+60DyvttseKb7HaQZp+yjYuDbQdYPewBRHxO1xK0M/k5ALNuCeCMbC51KLP4QCLK0sMnkCRmyidBnZHWRefN39SXeA5jqgJl/RzKZ85MRvAfL1zhEllwsCSD6C5+ANpF5TSS/eC/+DXv6LOw137FnvfgiQAIY37EteM9AZEBe0h3i8AJivxT2TSVfE5ONAGLPEsgllrwQtBFEku+9QgYESytY4kk+JZwOIOyhwtUYIGrITCCgJy8JYS8bAiLeCWLnPwWX7hoYnHrEcieuZUm+BiZRDhAgZERw8APZCCxv4/qNVCXdNTB+T1TytRLAAvsmEVRRVtchELwLRCjxPjj8G5b/sS8QwRoHYpDtyadExQx4IssDYoZ3+iuDQP6w1Oj3/SUTQOgD7qxEEMahi+TriJKvgUG2EDKBjOMqEdRJvNoGX8CBzuZyDUyuBNDfHYlDEFjaxLgnF8iAZXncjYiC65hgh/46G/889kvCHjCpwR1zwDBqwb4/7EXFIRB/wbP/bjaIIQgisu2wNWCFAxjJxZ9oxjYOHHLy9XHYBSsG6Ee6a2ByyQD69/rbCXDiHgQ68XodPGNJANGe3wcNfcCk1uslJkfwP9DVTPaKOphoY8KY6raEqGPF52ofAiSAtW/MV0uPQVpw+TMyY1iSwhUnmYr/HBgOSCBbgJKNAOIZtI2rGxCc4DgxY0cmARvbkc1BJg1ZPRwuSLxeJRsICPrYmI2sAjI4yDJgpo6/Y9aO39Eu9jFhadB/tQBBDY4cmQn0GX3AZno4cZx6Tr4IOlkOv4/Q/3umJbd0BBBLddi4n3ixLMgyMiHai6CxVwtZTxQQJ+COkhgkMR7JdygimOJEKbIZIOm40gZX4aANZOYwSch3qTIqZiA7yIxhTCEDyDUmBbiWwxMQkMRUn6wDOcIkAnqFyUDiRdC4Wif5ImiMP8gJ9BN9B+mHfoIUArNcL4LOlQBiTEA8seSKfXbILGFJExMVZFshM/QE2Xec+EXxmU4QCeAEuTHRwYQBkwdca4Qx8wQlmw3hd7+vEn9PdR8nTsx7fQK5wwQquSReBA1Cm3gRNMY000XQ+RBAvB9LqpgUgLhjD3PiARhrAggdgX6gn9BL+E3YK3wG7BWTCRzSSoeP37YCvfRfOIkyNnyGCJAAUgeqBQF/91um022JL/Z3p2EZDndfRSWAeA5BGQESmT4sIyIo4AsY2AcFx4pZtM9y5NJZBE4EfjhfBCksx8L543NeCDrYYI8/ySdD/afgIA/2NCHz4T8Fhw3xmchw4r5JBFr8f7qrODJ9Cs5/Wgqb9kE8saSFLBCuYIEMUT8tlYxX4uXc6bBMRaBwwARLaiCl/soPZHZwfUk+Y5P47lwwQ9YI16ngZCruTESABTHBv+Fybsie7pvFWI5DBgx74kBgQfpBllAfWcHk5Wu8C9lj6A8mI9AD6A7IIw4E4OCLL9ZLwL5d2BKIIMg/yCcmMDhsBEIIog79hUwo6D8uMAbJ85/Jgx1Bz7AHEDJHuaQ7cWz83X/4t+RLofFv/jOAkAtEKN2eW/8pOGS70CaeB6FFhhHX4CR+Ng7t5rsEnCg7CBd0AfaN5Wp/sbM1AcQ74VtwEAu6AntN/HQjfCgOpsB+UtkKCWAuXp3PJiJAAkh9IAJEgAgQASIQIAKYLGHVAP9FxhYTFhYiYIUACaAVkmyHCBABIkAEiEAeCKTaI4usILYdICuIK6xS3dOZx6tYhQhUIEACSGUgAkSACBABIlBABLCsjCuzsI0AWT4ckMJeURwYwuX0WJ73J5wLKCZfHTMESABjNqDsDhEgAkSACNQsBHCICPv8mjVr5g7e4EYEHDDDvl0cisNeTBYiYI0ACaA1omyPCBABIkAEiAARIAKBI0ACGPgAUTwiQASIABEgAkSACFgjQAJojSjbIwJEgAgQASJABIhA4AiQAAY+QBSPCBABIkAEiAARIALWCJAAWiPK9ogAESACRIAIEAEiEDgCJICBDxDFIwJEgAgQASJABIiANQIkgNaIsj0iQASIABEgAkSACASOAAlg4ANE8YgAESACRIAIEAEiYI0ACaA1omyPCBABIkAEiAARIAKBI0ACGPgAUTwiQASIABEgAkSACFgjQAKYAtHvv/9ezj77bPcB7hBL8+bN5d5775XS0lK5/fbb5bzzzgtGzK5duzrs8B3L6iyhj1GmvpeXl0urVq3ks88+M4foo48+kuuvv15Gjx5t3jYafPnll+WRRx6RG264oVrarymNXnfddbLeeutJ9+7dgxL5tttuk6OPPtp9SmyjjTYKQrbq1PcgOkghiEANRaBWEsAWLVrINddc44LwGmusIdOnT5e77rpLrrrqKjeMoZOLUaNGSZ06deSUU04JTu1IAAs7JBoCCBv47bff5JxzzknbidVBANu2betIJr6JGmoJkQDuueee8thjj8kWW2wh8+fPDxW6apXrrLPOkgsuuMD1/6STTpLXX3/dvW/dddeVN998U/bZZx/55ZdfqlUGNk4EagoCtZIATp48We6//3657LLLZPHixS6Tts0228iDDz5oTgDr168vS5cuNdWHF154QR5//PEgszBxIIDVMWamCpChMRLA1YN0iATwhBNOkHPPPbfas++rB+Hc31JSUuKy6tttt53stttucumll8r222/vGrr11lvlnXfekbvvvjv3hlmDCMQUgVpHANdff303A9xkk01k6tSpKYcVGcCbb75ZjjzySGnZsqV8+OGHAufqn7/yyivl2GOPlf/7v/+TKVOmyODBgyvIo89eDBgwQPBn1qxZsvvuu8v+++/vCCeyj9OmTXO/gcSlKiAgw4YNk+OPP17WWmsteemll+T00093cr/77ruyyy67OFK5ZMkSOfzww+XFF1+s1Azk2XXXXd27jznmGPffvn37uqyKL2eccYb06tVLGjduLBMnTpTTTjtNvvzyS/fzhhtuKDfddJPsu+++smjRIhk7dqzr47Jly8T376KLLnIzbWQisew0ZMgQVzeZADZs2FCuuOIKOeyww2TNNdeUZ555RvDuVBkKYD5v3jyHDcrs2bNd3zp37uz+//3333cYPvzwwy5Lm2mMmjVrJiNHjnSBYO7cua4/6ZYtvcxY8u/Zs6fLFGAZDWOMPm666aaCSQOyC2+//baT5W9/+5sMHDjQjRHw+uGHH1zfQcAwfsAHv2GZEO39+9//lhkzZri6K1askB133NFlnqFTW221lfz000/uN2SkZ86cKYcccogLWFtuuaVb0t1jjz3k999/F2R/gQHaQOndu7f069dP1l57bRfkDj300LRLwND5O++8070bMr711luu/o8//ujGZPjw4a5dLB/i37bddtsq6pmcAcRSP+rtsMMOMmfOHIFt3HHHHa4e3gPcMbmC7mDiAj3GcyjHHXec0xvoIPQB8mNMgQn0fsGCBe65gw8+WN54440qsgDfCy+80C11IvCfeeaZDn8UyImx2nnnnaV169Zu/DA+ftk9F71Ee3vvvbfTJ2TXnnvuOadT9erVq1gChk1Cv+AvMK4XX3yxm2SiACP0CzgAW8gFe0DZYIMNBGRyv/32c9iPHz/ejSeeS1UOOOAAhzH04ttvv3V2DRvB+GEFA3oJm8VkNnl5Gvb/wQcfONl8gVx169Z1+gk8+/fvL5tttpnrH7YRDBo0qOJZECysnMCXYXw++eQTOeigg+SPP/6QTL95ff/444+dHwFW0C/4VIz7+eef7/rtSybflBKUv/4Rfha2gvFu0KCB0zOMM/4ffYbcLESACKxCoNYRQHT9iy++cEEG++dAqHzw9bCAXMAxdezY0QVtEA6QEe9QEbief/55+fXXXx3BwqzyH//4hyMBIEhwyAjUffr0cU1iz96rr74qRx11lLzyyivOIT355JOOGH799ddV9BHkAe0imMOJIaCCbMLZ+uCWaR8WnCzICYI7Aj4C6AMPPOCC03fffefIHogOSA76CiKI4IkABVKJQA0SgmdAmJ966in573//K5dffnlF/8aMGePaBzkCFngf/i2ZAI4bN86RVbwDbaMv+P8uXbpU6Tf6jAwGyA5m7sAdxKZJkyaOSGEMEGgQnDKNEQIzAj2WwxDAQLpBPLFX8r777qvyXsgMuYYOHeqCK8gRloqgHwjUIMgg2hhTtIUxAelp06aNI6fffPON+3cEQugS2kCQQ4CDjoCwYayhGyiJARGTABACPIMCHcHft956axdkP//8cxfUEKhBlDAWIAz/+c9/HEGfMGGCtG/f3gV2jDvIA7YGpNoDiMAOPQU5AtGEbgDXAw880L071yVgjMWkSZOcPj300EOubZCjE0880U1aMIZ///vfnY1Bf6GDX331lZx66qluXLHcjKCMZTos0YEIg+RHWQIGIYMNwUaAH3QRZB1twHbRR9gdfoeMwA9jBMxQctFLYAR9A7benkCwoEvwCZAdOgD9AYn19g17BcnGBACywn5A0P75z39WLE1CdvwOm8d4o12Q3UTi5RUWk5pPP/3UETXoNnQSNge79hOQTPtvE3ULbUIW+DfoOOSEHoFUgiyD0D/77LPOh2G/MSZ6GEdgiS0C2NcHO33vvfecXaf7DUQ2mQBinGAb0Bn0BcQaNg7Cn803ZQre2DIAe8HkF6QbPmmvvfZyEynYKfrFQgSIQC0ngAhcmHXC4WH5F0EJ2R0QHxQ4e2StkNlCAeHDzNgvJyQrELIOV199tXOUCF4geQgayGahjBgxws3oPSHEv91zzz0u43bJJZdU0UeQQjgvPytGhgNZBThJOOxs+7BABMrKylxg8AXEAYEGyyIgR3DCCCK+gBB36tTJEURkKIERCBcKnCcyNSAlvn/Iev3888/ud2QhkJnAn0QC2KhRI0ck8V8EexQEZQQRBLvly5dX6juyIXg3SOfJJ5/sMjsgr8jE4d0IisjoZBsjBGD0F+0hOKEgqwhS50l04oshMwgdnveZtSeeeMKRmRtvvLHiUQRmBHiMHYIVZPN7jBLbQ3DEfixkSFCQjVi4cKFsvvnmbuKRGBBBepHtBXlCAbEHCYJegKBjnHyf8TvIHcapXbt2jrRiCwPIDwqIK7LE0OUoh0AQ5BEcQcYgU64EEIQaWCNT7gvkBlFNtT8VEyrYCYgY3gn9AZkAkQJmvkQhgCDnfmLh68GeQMLQHmwEffPZZMiJSUBRUZHTx1z0EmQF7STbEzLrIIDwD7BXTKB88YG7vjIAACAASURBVL4Dkyj4A/gYjDP02xcQFciUqHcYV+gY7CS5QBdAaJAd9gUkDe2DXGbbfgHSD/+B+iBsIJAYD5DmVAUTDWTQQNgxgcHEFrJiopNYMv2G55IJIPwu9l/7gvYwPlhpyeSbUmWBk+WGzYCoYxIAO4DOgbxiMonsJfwOMrWJqyEpO89/JAK1AIFamQFMHFfMGrGMBEeNbJbPLiWeAoYTQRYGyz8o+A0BrmnTps65rbPOOi5zBbKA4AVihayALyATWOJBsPYFwRpLq8iMJRcs9YGs/O9//6v4Ccs6yDhhxh2FACLLkBgobrnlFicr3gfyAvKFZTlfEBxAukAAEVAQoH1BW8joIBigf08//XSl30HQQDoRABODEAIc5PXkz7eHpWBkM0BqkwsyHHDgWJJC1gaEDTN3EHWQaH/iOfmgTuIYgfyCsCYGZMgIApkYxP27ITPIOQiRLwhEIGyeQOLfkTHBUhKybyDHIKggEokF5ALEBuTfk0n8DhKIbBcyPokBEf+ONpCF8+QbZBdLZOgrsoEgj75guQ7bDrA8C5KLsfKHl/AMstuYvKQigJANwQ/ZM6+fIET4OwJmrgQQExvYAXTTF2RfQYqRecMYg1hjGR42AtmBJ7KBKLAJ2A3IMsYdWTD0JwoBTNV36CX0FKQm2UYwtsjkggzkqpeY4MAek+0JegwCmOp36DDsFThguRb2AV2GfwFuyHqBrGCZOJH8Qj5giMxpcoE9oCT6DNg1bBmTgGwEEHUTnwcpAiYgpijQQcgJgg5dh24CUxB8TFSg+7DD5JLpNzybTACxNeCII46oaAaYgIxilSSTb/JL6lUESPMPINHI9CJT+dprr7lxgp5hYobJfLJfitounyMCcUGg1hNADCSCE5wwlu0wC81ELhCssJyD4IXMH5wb/otlYATXVMELThfOxmcjsilPcgYQ2TiQBE0GEDJjmQcZQJAEkFhkD5IL3pGcAQR5QpBIlwFENhXBIzkDCLnRFoJZIknI1H+QaBBTZF1AwIAzlouRIUNQRT9QMo2RzwCC8PgDOKiLttJlAJOXzhD4kCXwmZxkmRGokMlIzkoggHtdQtYnVUkMiPgd2SxkQfA8yKtfKgbukCsxW5LYXq4ZQCxhA5MePXq4TKEnRT5bjeVNEMGop4CBKYK536OZ3FdsDYAuI3MFQgySDjtJPt2LyRBIDXQMvyHLhUlUplPAqTKA0Gu04TOAidskEglgrnqZKgMIPfTbQlJlAGHzILyYWCYW+A+sNOC/IFlYwofNRSmpMoDQU5CbKBlAvAO2AWz9sjEmSbAlyAIihnEA0cJkCxlA2CDIms/yQX8SJ7JoM9NvuRLATL4pCkaJz0D/MBHEygcm1dBz9At+ENlBTE5ZiEBtRqDWEUA4AWQdsIyHzBJmuSAw2AOHDCAcRSZygWU/1PWb+EFO/F186QggnsVSD5bu4KwR8LCsB1LoD14kKiEyVZh1d+jQwTlltI+lF79XK0oGEEtS2E+DrA5ID/bb4HQc9iphIz72aeEPAjQIGvZGIXsCYoOlHhA3ZOGwHItg509N+z2OCORoB5hhqRSOFv+WnIXAniYEf8y+sR8OwReEJt0SDIINZMYyMYIkMlTY3wSSDll8tiTTGAFf1EefIReWuEB2IQOW6ZNLqswJsEfWF4QMkwIsHSF4YryADYIj/h/BH/umEvcAYqkJmRLghyVfZLywtOeX9JMJIAgPsjHYPwjiACKGgowrMhbIaCHriOwZshrIPCJbgowi+gi9wCQEGTRMMtLtAfT73tBf4AoCCbw9AUS2EcQARC1dSdQ9EBe8F3rmtxOgPsgElrGxzIjnIdPGG2/sdAi/g9hhCwHwAxmCznXr1s0t0UHPMV7Yy4V++m0GyfIgu4b9k7BHvAcywG5QF/qWKQOItnLRS8gLfUNGNtGe0B9kAIEf7Ar4w1ah38hQQjaQD9gZ9A+EEf1HZh99x9giIwz7wb5R4AB7AimDv0gufg8g/AhsEnv34IuQDYZ8UTKAaBN+D/LCppANRsHf4Y8w9n5/MvAFWYKOYGIDwoSDHMiWQ9bEPYDpfku1BzBTBjCbb0KGEtlYv5cznZ7CLwNj2CAKMvqYsMCWkQHEOPgtLrWZALDvtRuBWkcAsbSJJRgQGQQhZF7gFECY/IwwE7mAI8QeHWRnMBPGMi6WlDCTT0cAoWJwWFhCQSYLSzZYekFAgUNNLgig2EsFcoFlJgQzOEbvsKIQwMRTwKiHPYwgC74gYGLZCCdDQaqQyQJxgGMHSQNGcLTI3OEACIIbsmnJp4CR5UDQg2MGsUkOQggs2JeFIAICh31TICJoL1VBsEV2CmTEP4O9XGg7MROWaYzQLogA+uBPAWPJ7dprr035znSBE0t0IC9YwsNYI3ADMyzBYskc5BLZL/TLB2CMK8YPZBPtYj8ciC9Itd8Xl0wAIRRIJAiPPxHrBcW7scQLMgxdwHMghMAQBScm8S4Q1GyngEFKsTQMggFiCjwwdp4A4l0gqfgv+pi4JO7lSdY9BHOQF0xooAvI4ICIYTIBmZFBxeETTDRAVvAbxhj9RKYO9VEPv2MihiVgFNTDJAhkHmQcRCq5IMhjfPwpYGR1cBgGJRsBzFUvoffQJ2TEkFnCQaDEU8DQM0wY/ClgZNphNyjAHCQd74T+A3e/nAvCC/yQPQcpxyQAfce7UhVknZHtwxhhywBIDQgkSlQCiDGATWJvH7LCviBbid8gJ8YB2xDgH/xyLXCG7PBlmDhD10Fy4UMz/ZbLEjBkyeSbMDnCPkb463QF9gj5MbHy+7AxfphEoW/ww+nwrd10gL2vbQjUOgJYGwYYZCx5lm3V7yj7s6zexXaIABEgAokIIHuHiam/SojoEAEikD8CJID5YxdsTRLAYIeGghEBIkAEiAARCAIBEsAghsFWCBJAWzzZGhEgAkSACBCBuCFAAhi3EWV/iAARIAJEgAgQASKQBQESQKoIESACRIAIEAEiQARqGQIkgLVswNldIkAEiAARIAJEgAjElgDi2hNcI4GrL3CVCe6zwpUJid/exVUROFGWWHCVBq4hYCECRIAIEAEiQASIQFwRiC0BxA35uKgVd/vhLjHcK4cLU3EHGi57RgEBBCFM/PA6fkv8NFNcB579IgJEgAgQASJABGovArElgMlD6r/Rii8I4FulngDiMtNMn76qvarBnhMBIkAEiAARIAJxRaDWEEB8RgmfP0IWEJ8J8wQQN/fj6x741i4+fYRb4qN+tzauSsF+EQEiQASIABEgAvFGoFYQQBA8fKsUn7zy377EsPbo0cN97mj69Omy/fbbu08y4XNfRx11VMpRx+e/8AmkxILvvPJW+ngbCXtHBIgAESAC8UPg73//u4v/tbXUCgKI727im5X4NuS0adPSjjW+cYlvmCJbiO9sJhdcsIzvv7IQASJABIgAESACNR+BjTfeuNaSwNgTwJtuukk6duwo2Pv3ww8/ZNTWtddeWxYuXCgHHXRQxQfWEyskZwAxewChbN68ucydOzdvS2jcuLFbgs63aOvjvdo2tPX/9re/SefOneW+++6TJUuW5AWFVoYQcNDKYIGjVgaL+hZtaPWBWK4ywxCw1MoQgk5pZaBO2ulkcXGx2xZWVFRUaw9+xpoAgvwdccQR7qoXDHS20rp1a3nzzTfdcvCnn36a7XEBAZw/f76sv/76qmXgjTbaSGbMmJH1feke0NZHu9o2tPXh2Lp16yZ333133gRQK0MIOGhlsMBRK4NFfYs2tPpALFd5nBCw1MoQgk5pZaBO2ukktm/9+uuvJIB5M4+AK44cOVKOO+44l/376quvKiSdN2+e/PHHH7Llllu635966imnBCB91113nUydOrXK3YDpukkCaGeMFo6NAULEAkdtkLKob9GGVh+IJe072fdrdUqr19RJO50kARSJbQZwxYoVKXkbskyjR4+Wpk2byj333ONOBTds2FCmTJkiEyZMkEsuuSRyOpgE0M4YLRxboZ2z1rl7NDX9sMDRoh+aPljgYNEHYkn7JgFMnf6Ig32TAMaYAK6O5CQJIAMEA0R8AwQJIO2b9h1f+yYBJAFU8UQSQAYIBoj4BggSQNo37Tu+9k0CSAJIAvgXAtqUvra+RbDVymCxbFhoGSxwjAMOFn0gliSAJIAkgCqSEHjl2O4BXB24MwPIAMEAEd8AQQJI+6Z9x9e+mQFkBlDFE0kAGSAYIOIbIEgAad+07/jaNwkgCSAJIJeAK+mAdglXW1+7dEnSEhZp0Y6nRX2LNrR6baGXWhlCwMHLgIv/Tz31VOnVq5e7kmz27Nny6KOPCr42VV5e7pQYX6TC1WT4fOnSpUvlgQcekIEDB7pPlfr7UvHMZZdd5r5yhdsscH8t2sAXrVKVDh06uG/eN2nSJNLds2gTn0jFO9dZZx157bXX5JxzznEyprq7Fs9AhldeeUW6d++eMb4mj+c111wju+yyi+CLXIkFH2jAF7iOOeYYd+fuZ5995nD44osvIvUhnRAkgCSAJIAkgCSAKaxAG2y19UMI1hakJYR+hCCDBZZx0Ck/Fl26dHFXjl199dXy4osvSosWLeTiiy+WDz/8UA488EBZd9115ZNPPnFfiLr00ktlww03lKuuusp9q/6hhx5yBBArUHgG99iiLXyUAN+3P/zww2W//fZzZC2xtG3b1pE/1ItKAEFK99hjD+nbt69rH+Rygw02kP3331++/PLLKp5j1KhRcsoppzj5ciGAffr0keHDhzvimEwAx44d6z7ocMEFF8jnn3/u+oj/Rz9xj2++hQSQBDBf3XH1uAQcVrYlTgEi3y/DWATaEAhDCDIQS9p3coCw8DEgX8hi3XvvvXL66adXvKKsrEzGjRsnu+66qxxwwAEuy7XZZps5gofSvn17efrpp11G7qKLLnJ1r7jiCtliiy0qvmVbt25d+fjjj91nT//1r3+5esjK9e/fX/r16yf4EAKyaFEIIIjf22+/LQcffLA888wzrq1GjRrJ999/LyNGjJABAwZUggfP3X///bJ8+XJ55JFHIhHABg0aOOJ32GGHyYIFC2TixImVCOCaa67pso2XX365DBo0yL2vXr16Toa33npLOnXqlHcMJwEkAcxbeUgAK0OndYwWwVYrQwikQyuDBY5aGSzqW7Sh1QdiSQJYHQRw6623dku/IHsgMb7ga1Qgb8cee6ycdtpp7otVIFW+gNzhm/NvvPGGy37hz84771yFiD388MNSWloq22yzjat61llnuewZSCNsAsupUQggsn3nnXeey0aC1PnyxBNPuIzk7rvvXvFv6623niO1yGKCbEZdAr755pvdUnfv3r1l2LBhrr3EDCC+04s+g7xiidgXLDNPnz5dDjrooLxjOAkgCWBeygPjhcLCIGFoLVu2dEqabykpKZFZs2blW1209fFibRva+gi2mAGPHz8+728Ba2UIAQetDBY4amWwqG/RhlYfiOUqlxQClloZQtCpTDJgGRSECwQIfvCxxx5zpC2xvPDCC24/ILJ7S5YsqRIzQNZAKvFNe+wxREEWEXsMFy1aJOeee677s9NOO2WNObfeeqs0b95c2rVrV+k9IGrYE4i458tNN93kSCHI67vvvutkwF7BTAXjCeLoP9X64IMPusePPvroStWuvfZaadOmjfTs2VMmT54snTt3dkvRWJbG17zyLcXFxTJp0iR+CzhfAGt7PS4Br9KAELItWhnQG20b2vpaGZi1CksnteNpUd+iDa1eW+ilVoYQcEgnA7JpyJo9//zz7vv1IGvIeCUTQGT/QFx23HHHKgSwTp06jjgiM9i6dWt57733qoRoEKeoGUAsN2P5GIdQEguyfCBfWL5FwfvGjBnjPqv6008/ueXZqBnAxK0uL7/8smsveQ8giOKTTz7pDoj4AlzuvPNOHgJRkjDeA6gAkAQwrGAb5wARVU0tAm3IgTIqDhZ9IJa072R9qw4fA7KGZVWQIZCtOXPmyOLFi92hj1QEEFk+LP0mZgDr16/vDl4cf/zxbnUKS6upSioCCOKI1azEsmzZMnn22WcFJ3BTEUBkKtdaay23JxBZNOzPu+2221wTyQQwXfvJWKYigDhw8v7777u+QvZp06a5vZDIluJwzNChQ3NxCZWe5RIwl4DzVh5UJAFkgFgdASIXJSVpCUsnLYhodZCOXHTKog8WehkHHJKxxLYXELevv/7aERuc+kXBki0OiZx99tmVhgqnfrE3cM8996wggCCE2PeHU77Y7zdy5Mi0w5uKAPp/S6zks4k4nYxsY2K54YYb3DIslnxxNQ0OleDk8ooVK9xj33zzjbz66qty8sknC4hkuvajEEAcNME1N1tttZVr1xcchEFfsZcRhDmfQgJIApiP3lTUIQEMK9jGMUDkqqAWgdYi4MdhLIgl7bs6J3jYi4csH5ZLca0Jrlnx5fXXX3cndnFvny/I0v3222/ueezBQ1Zs4403dsvGOAl84oknit9Hl85vpCKAsFUQqcTywQcfuAMdZ555ptun58kdnsFyLP4NJDTx31O9c5999nHkNlX7UQgg9iFiiblx48aVmj/00ENd1vSf//ynuxonn0ICSAKYj96QAKZATRvwLYKtVoY4EB8LHOOAg0UfiCUJYHURQBziwJIprk3BnYDJBzqw9It9diB2v/zyixPDXwODq1+w5IorUrDPD+QIRBH7A7OVXPYAYukXdwmmugYGhz5wsjhxX55/Nw6vgEBieRYHPHC9S6oShQAiA4prYv7xj384IukLsoLABwQ43wOUJIAkgNnsJePvzAAyQFRXgOA9gIU/kEMCSPuuDvvGdS+4Xw/EBVk7nOpNLN9++637X1x6jD1vIFJYZkW2ECdsJ0yY4JaNsQyKk7Ygddivl1iwhxB36iWXXAgg6uKLIjvssIMjW7iPEAdIIAvuKUx1ETTqWB4CQYz96KOPBEvSePfUqVPdu7EHcfTo0e5S6HwLCSAJYL664+qRADJAVEeA0GQySVrC0klIoxlPi/oWbWj7YKGXWhlCwAEygLThapN0pVu3bo7c4JqV66+/3p3oxWXIuFwZe+JwTQoIIPbEbbrppimbwUXQyB5qCSCWeiErlmGxBI3rZTJ9Cs6aAKI9ZDiR9TzkkEPcoRRcBYN9jliKzneijHZJAEkASQD/QkDrXLX1GSBsiIsFjqEESq1OaesTSxudRCsWWGrHMw56bYFjHHCw6AMJIAkgCSAJYCUd0AYZbX2tY2OACIu0aMfTor5FG1q9ttBLrQwh4KCVwQJHrQwW9S3a0OoDCSAJIAkgCSAJYAor0DpXbf0QAgSDbVhkOg46pdVr6qSdTpIAkgCSAJIAkgCSAKb0Awy2dsHWAksSQJuldC0Jtahv0YZWH0gASQBJAEkASQBJAEkAs3hCbbAlAbQh0xY4hkC+QpCBBJAEkASQBJAEkASQBJAEMFIs0BJhLfEhAbQh0miFBJAEMJLRp3uI18DYGaOFYyu0c9Y6d4+mph8WOFr0Q9MHCxws+kAsad/Jvr/Qek2dtNNJEkASwLwIYK9evdwHt3EvUmlpqbuvae7cuXm1hUolJSV532ZuUd+iDW0f4NjwXczx48dXuRU/KrBaGULAQSuDBY5aGSzqW7Sh1QdiucryQsBSK0MIOqWVgTppp5PFxcUyadIkKSoqcvcs1sZSR0RWfsGZJWcEmAG0m41ZzGwLPTu3yDpp27DAUSuDRX2LNrT6QCxp38lBQatTWr2mTtrpJDOAzADmTPoSK5AA2hmjhWMrtHPWOnePpqYfFjha9EPTBwscLPpALGnfJICpQ2Qc7JsEkASQBPAvBLQGra1vEWy1MliQhkLLYIFjHHCw6AOxJAEkASQBVJGEwCtzCVgxQMwAMkAwQMQ3QJAA0r5p3/G1b2YAmQFU0D8REkAGCAaI+AYIEkDaN+07vvZNAkgCSALIJeBKOqBdwtXW1y5dkrSERVq042lR36INrV5b6KVWhhBw0MpggaNWBov6Fm1o9YEEkASQBJAEkAQwhRVonau2fggBgsE2LDIdB53S6jV10k4nSQBJAEkASQBJAEkAU/oBBlu7YGuBJQkgvwWcaKhafSABJAEkASQBJAEkASQBzOIJtcGWBNCGTFvgqM1CWtS3aEOrkySAJIAkgCSAJIAkgCSAJICRYoGWdGiJDwmgDZFGKySAJICRjD7dQzwFbGeMFo6t0M5Z69w9mpp+WOBo0Q9NHyxwsOgDsaR9J/v+Qus1ddJOJ0kASQBJAJkBZAaQGUBmAJkBjBQLSABXwlRoHCxkIAEkAYxk9MwAZodJ6xAsZrZaGSycSqFlsMAxDjhY9IFY2mVbLLAstG1Z6JS2DQsctTJY1LdoQ6sPJIAkgNmZTYYnuATMAMElotQGonXOIQQIBlvaN+07vvZNAkgCmBcB7NWrl/Tu3Vvq1q0rpaWl0rJlS5k7d25ebaFSSUmJzJo1q2D1Q5ABwbasrEzGjx8vS5YsyQsLLY4h4KCVwQJHrQwW9S3a0OoDsVxlhiFgqZUhBJ3SykCdtNPJ4uJimTRpkhQVFUl5eXleMaemV+K3gBUjyAwgMwTMEMQ3Q8AMIO2b9h1f+2YGkBlABf3jt4ATwdMu+VkEW60MISw7amWwwFErg0V9iza0+kAsSQBJAEkAVSQh8MrMACoGiBlABggGiPgGCBJA2jftO772zQwgM4AK+scMIDOAVdVHm3XSZr5IWsIiLdrxtKhv0YZWry30UitDCDhoZbDAUSuDRX2LNrT6QAJIAkgC+BcCWmPS1rdwbFoZQnBKWhkscNTKYFHfog2tPhDLsMi0djxD0CmtDNRJO50kASQBJAEkAaykA9ogo63PAGHn4LVjwWBrNxYWWGrHU2tbFvW1bVjgqJXBor5FG1p9IAEkASQBJAEkAUxhBVrnqq0fQoBgsCUBTDaNQus1ddJOJ0kASQBJAEkASQBJAFP6AQZbu2BrgWWhyZfFpETbhgWOWhks6lu0odUHEkASQBJAEkASQBJAEsAsnlAbbC2Ii1aGEEiHVgYLHLUyWNS3aEOrDySAJIAkgCSAJIAkgCSAJICRYoGWdGiJDwmgXVaaBJAEMJLRp3uI9wDaGaOFYyu0c9Y6d4+mph8WOFr0Q9MHCxws+kAsad/Jvr/Qek2dtNNJEkASQBJAZgCZAWQGkBlAZgAjxQISwJUwFRoHCxlIAEkAIxk9M4DZYdI6BIuZrVYGC6dSaBlywrG+iDQQkcUisrTyGGv7oa0fwljkhGUGE9Fioa0fFyzjgIN2LKiTzABmj8bRn+Cn4KJjVeVJLgHbGaOFY2OAEImE46YisoeIlIpIXRFZLiJfisjbIjLFZoYfh7GIhGUE/6HFQltfSzos6ltgGQcctFha4KiVwaK+RRtafWAGkBnACO676iO9evWS3r17S926daW0tFRatmwpc+fOzastVCopKZFZs2YVrH4IMsCxlZWVyfjx42XJkiV5YaHFMQQctDJkxXEHEWn3F+kD+fMFJBD//4KIfEydBCxZsYyopVq91NbX6pRFfQss44CDFksLHLUyWNS3aEOrD8XFxTJp0iQpKiqS8vLyiNYcr8eYAVSMJzOAq8DTzsYsZrZaGUKYlWplyIgjMn/dZeW0L11ZISL/Edlo6UYyY8aMvK0jDmNhoZPa8bSob9GGdjwtsNTKEAIOWhkscNTKYFHfog2tPjADyAxg3gEOFUkASQCTFUjrlLSOMWOAKBORrUWkXga1XyYiX4ls9DoJIIMt7btG2XcO0Uzrp7T1tX7Ooj4JIAlgDiZT9VESQAaIGhMgcODjApHmjZpL06KmsknRJtJ9R6QDV5bL3rhMXvgO678rl4c3umsjmTGFGcBu3brJ3Xffnfe2BItAFYdga0Gm44CDVh8scNTKYFHfog2tPpAAkgCSAP6FgNaYtPUtHJtWhhCcklaGtDg2FDlvwnlySPNDUur8XRPvkrGfjK34baOxG8mMb0kASQBXqoTWtmjfNpNlCxwtxlOrDyHIQAJIAkgCSAJYSQe0jk1bX+sY0waI+iLXvXad7NB4B5kyf4p8NvszKf+zXGYvnC3fzvlWPp718SocmAF0WDDY2pAWKywLbVta2/RoavpBnbTTSRJAEkASQBLA2kEARaTff/vJQa0Okjs+ukPu/fTe1LrPPYAVuDDY2gVbCyw1xMmCfJEA2umDBZZafSABJAEkASQBrDUEsGufrtL19K7yxOQn5Nq3r02t+zwFTAKYQjO0wZYE0IY8WeAYAvkKQQYSQBJAEkASwFpDAA866CDpd0k/eX/x+9L32b6VTwMj84d7AJ8Ukff1e760hCGEAMFga0NauARcOcxobIM6aaeTJIAkgCSAJIC1hgDusMMOct1118nUuVOly7gu/BJIFutnsLULthZYaogTl4DtSKjF5MyiDa0+kADGmAD2799fjjzySPeljkWLFslbb70l/fr1k6+//rrCEho0aCDDhw+XTp06Cf7+7LPPCr7yMXv27EjEkNfAMEAkK4rWKWkdY6ZAi5vz77vvPlm6dKkceOCBIvwWcEY7tyAt2vG0qG/RhlavLbDUyhACDloZLHDUymBR36INrT6QAMaYAD799NNy//33y3vvvSf169eXyy67TLbddlvZZptt5Pfff3eO/+abb5ZDDz1UcNXDvHnzZMSIEbJ8+XLZa6+9SAAjIUACWJMIID5d+Nxzz7lPGB599NEyZ86ctKOsda7a+iEECAZb2ndNsu9cXLbWPrX1Q7BvEsAYE8BkY2jUqJH8/PPP0qZNG3n99dfd9//w/8cdd5w89NBD7vGtt95avvzyS9ljjz3k3XffzWpPzAAyQNS0AIEMIDKBZ5xxhvsOZrqidfDa+iEECBJA2ndNs++sQeuvB7T2qa0fgn2TANYiAtisWTP55ptvXBYQgW/fffeVl156SdZbbz2X/fPlhx9+kOuvv979yVZIABkgouqg3gAAIABJREFUalqAePTRR90nDHv27CmTJ08mAcxg5CSAtO+aZt/ZYpb/XUvgtPVJAKOOVPU+h8/C4+KHWJc6derIY4895sje3nvv7frauXNnueuuu2TNNdes1Hdk/l5++WXBHsLkssYaa7i9gr4gkE6bNk0aN26ccTktG7haY9LWD8EYEWy7dOkiY8aMyfuzW3HAQTsWmXDcfPPN5fbbb5c//vhDDj/8cLfdIV3RYqmtr8XBor6FTlrIQSxXamkccND2gzppNylBBnDmzJluNbC8vDxbmI7l77WCAGKv38EHH+z29oGw5UsABw8eLEOGDKmiCOecc44sXLgwbwVZd911K2Uhc21IWx/v07ahrV+vXj1p3bq1O6yzbBnuJMm9aGUIAQetDJlwxCngdu3ayY8//igPPvhgRoC1WGrra3GwqG+hkxZyEMuVqhoHHLT9oE6ucltafWjYsKG7FYEEMPdYW2Nq3HTTTdKxY0e39w/Lu77kswTMDGD6YdfOzi1mtloZtLNzi/raNjLhOGDAALf1YezYse5PpqLFUltfi4NFfQudtJCDWK7U1DjgoO0HdXKV19LqAzOAMd8DCPJ3xBFHyD777OP2/yUWfwgES8EPP/yw+6lFixby1Vdf8RBIHvRea4wW+620Mmids0V9bRuZcMSp+A033FDOPfdc+eijj0gAs+i5hU5qx9OivkUbWtuywFIrQwg4aGWwwFErg0V9iza0+sBDIDEmgCNHjnQnfJH9A6nzBQc+sAcKBUvDhxxyiLsGZv78+QLCiLLnnntGokA8BGI3G7NwbFqHEIJT0sqQDscNNthAxo0b5/b9dejQocIG0im6FkttfS0OFvUtdNJCDmK5UkvjgIO2H9RJu5hDAhhjArhiReqzLSB7o0ePdlrkL4JGFjDxIuhZs2aRAEZCwM4YLRwbA4RIOhz3228/GThwoLvmCJedZytaLLX1tYHSor6FTlrIQSxJAL29UiftYg4JYIwJYLYAZ/E7M4B2xmjh2OIQKLWEIR2OZ599thx22GHywAMPyC233JJV/bVYautrcbCob6GTFnIQSxJAEsCqLktrFySAJIBZA2GmB0gASQCT9UPrlLSEIR1pufPOO2WLLbaQQYMGyRtvvJFV77X90NbX4mBRnwSQ9l1T7DurQSc9oLVPbX0L+9TKQAJIApir3VR6ngSQAaImBAgQPxBAFByKSrz4PJ0BaJ2rtn4IAYIEkPZdE+w7nyCmtU9t/RDsmwSQBDAf26moQwLIAFETAsSVV14pu+22m7z66qsydOjQSDqvdfDa+iEECBJA2ndNsO9IBs0MYBWYSABJAPOxHRLAFKhpA75FsNXKEALp0MqQjOOuu+4qV111lSxdulS6du0qM2bMiKTzWiy19bU4WNS30EkLOYjlSpWNAw7aflAn7SYlJIAkgJGCYbqHmAG0M0YLx8YAUfkUMEjfqFGjZMstt3Rf/sC1R1GLFkttfW2gtKhvoZMWchBLEkBvt9RJu5hDAkgCGDUepnyOBNDOGC0cWxwCpZYwJOK4//77S9++fWXBggVy/PHH5/S9Sy2W2vpaHCzqW+ikhRzEkgSQBLBqCNbaBQkgCSAJ4F8IaI1JW98i2GpliEOw9jjed999ctddd8n666/vrn3B9S+5FC2W2vohjIWFTobQjxBksMAyDjqlHQsLHLUyWNS3aEOrDySAJIC5xMQqzzIDyAxgslJonZLWMfoAsWTJEunSpYvMnDnT/RfLwbkUbT+09bU4WNRnsKV9h2rfd999t8DG8y1a+9TWt7BPrQwkgCSAedkPvqTQu3dvqVu3rpSWlkrLli1l7ty5ebWFSiUlJRL16yOpXqKtH4IMCLZlZWUyfvz4vB1bHHDQjgVwxGEPfPkDf7/++uvlzTffzFk3tVhq62txsKhvoZMWchDLleobBxy0/aBOrnJlWn0oLi6WSZMmSVFRUU7bY3J2pgFXqCMiqb+ZFrDQoYjGDCAzBCFmCO655x7ZcMMN5Ysvvoj02bdU9qSdXWvrh5AhYAaQ9h2ifeNzpswA6k+FMwPIDKCKS5IAMkCEFiCaNWsmjz/+uMson3nmmfLpp5/mpeNaAqetTwJoZ1txwTIOOqUdC05K7OyCBJAEMK/g6CuRANoZo4Vjq+0Bol69eoJLnzt06CAPPfSQXHjhhXnrtxZLbX1toLSob6GTFnIQy5VqHAcctP2gTtrFHBJAEsC8AyQqkgDaGaOFY6vNAWLttdeWYcOGyS677OKWf0ECv//++7z1W4ultr42UFrUt9BJCzmIJQmgN2TqpF3MIQEkAcw7QJIAVoZOG6QsHJtWhpoarDfYYAOX+dt8881l8eLF8tFHH8ngwYPzPkxTU3FINmatPljoJLFcOSoWWGrHMw5jYYFjHHCw6AMJIAkgCeBfCGidq7a+hWPTymDhVFa3DM2bN5fLL7/c3ff366+/yqBBg2SvvfbiJnGDJUMLnayJOpXKKWr12gJLrQxxGAsLHOOAg0UfSABJAEkASQAr6YA2yGjr5+LYdt99dxkyZIisueaabrm3f//+8ttvvwlPCdosGTLY2i23WWC5Om0rXWAotAwWOObiY0LFwaIPJIAkgCSAJIA1kgBij9/ZZ5/t7qL84IMP3JLv77//brLUZuFcCx0oLfrAYEsCaL2tQKuX1Ek7nSQBJAEkASQBrHEEsEePHtK5c2cn9zPPPCPDhw+XZcuWuf9ngLALEMQyLCzjMKkgAbTTKa0+kACSAJIAkgDWGAIIQoJl3n333dfJjG/9jh07tpL8JC12AYZYhoWlNuBryZdFfW0b1Ek7nSQBJAEkASQBrBEEEFcOXXLJJbLddtu5bN9VV10lzz//fBX9ZYCwCxDEMiwsSQCZ4U90eFp9IAEkASQBJAGsEQRw6NChsvfee8vChQvloosukokTJ6bUXZKWsEiLNuNjUd+iDW2wtdBLrQwh4KCVwQJHrQwW9S3a0OoDCSAJIAkgCWCNIID4vm+TJk1kwIAB8u6776bVWwYIEsBUyqENltr6FnqplSEE0qGVwQJHrQwW9S3a0OoDCSAJIAkgCWCNIIDjxo0TXPjcs2dPmTx5MglgBMvVBggG27DItHY8QyAdWhmok3Y6SQJIAhghjFR9pFevXtK7d293BUdpaam0bNlS5s6dm1dbqFRSUiKzZs0qWP0QZIBjKysrk/Hjx+f9BQstjiHgkE6GO++8U4qKiqRPnz4yZcqUjARQi2PIOORiJFp9sNBJYrlyxCyw1I5nHMbCAsc44GDRh+LiYpk0aZLzq+Xl5bm4ltg8W0dEVsSmN6u5I/wWsN1szGJmG8cMAS553mmnndw9f2ussYaceOKJMm3aNGYAI9i6Vh8sdFKb8bGob9FGCFhqZQgBB60M1Em7mMMMIDOAEcJI+kdIAO2M0cKxxSlAICM8cOBA91m3+vXrO6BXrFghHTt2lAULFpAARrBcrT5Y6KQ24FvUt2gjBCy1MoSAg1YG6qRdzCEBJAGMEEZIAKOApHXOFo5NK4PWOVvU923Uq1dPxowZ46CfMWOG/O9//5MXXnjBLVlkKhY4WvQjDmNBLO2CrQWWcdAprW1Z4KiVwaK+RRtafSABJAGMwm3SPsMMIANEsnJonZJ3jMuXL5f77rtP/vzzT2nfvn1kPWWACEsnQwh0IchgoZdWtoUJVb6l0DJY4BiCPoQgAwkgCWC+fsDVIwEMK9gW2jlbODXfBr7rO2HCBAfwfvvtF1lPGSDC0kkLnYiDXlvoZRxw0OqDBY5aGSzqW7Sh1QcSQBLAyIE11YMkgGEFW61DCMEpeRlwqvypp55yAB988MGyePHiSLrKABGWToakU4XMfFnoZZzsO9+xsMCROrnSR5AAkgBGCqrpHiIBDCvYxilAzJw5U1588UUH8BFHHCHz5s2LpKsMEGHpJIPtyvGw0Ms42TcJoIh2PLX1SQBJACMFVRLA7DBpjZEBoipxefbZZ13g7NSpk8yePTv7IBgFWpIWO9JCLO2w1PqYOIyFhZ+MAw4WfSABJAGMFFRJALPDpHXOFo5NK4OFU7GU4dFHH3X7TLt16yY//fRT9kEgAayEkXYsLHQyNJ2KpEQpHgoBS60McRgL6qRdhp8EkAQwX3/o6nEJ2M4YLRxb3AIEvorSqFGjrJ9/S1RiCxzjECgt+kAsad/JAaLQPoY6aaeTJIAkgCSAfyGgdWza+haOTSuDBWmwkgFLvg8//LCbZOCzg1988UUkXbXAMSQc8t0rZdEHYmkXbC2wtLKtQuqUVi8tcNTKYFHfog2tPpAAkgBGCqpcAs4Ok9YYLRybVoYQnJKXYYcddpC+ffvKnDlzpHPnzpG/j2yBY0g4FDJYE0sSQGYAU/v+OPhaEkASwOzMJsMTXAJmgKiuAHHFFVfIJptsIrfddpuMGzcusp6StISlkyTTK8fDQi/jQDq0+mCBo1YGi/oWbWj1gQSQBDByYE31IAlgWMFW6xBCcEqQAd/7Peuss9w3f4899lhZtGhRZD1lgAhLJ0PRKa1taOtb6KVWhjiMhQWOccDBog8kgCSAkQNr4oO9evVy+7Lq1q0rpaWl0rJlS8HFvfmWkpISmTVrVr7VRVsfL9a2oa0Px1ZWViY4+LBkyZK8sNDKEAIOkOGGG26QJk2auD2A+BxcLsUCx1Bw0I6ntj6xXKV5IWCplSEOek2dtNPJ4uJi9231oqIiKS8vz8XNxubZOiKyIja9Wc0dYQYwrGxLHDIEO+64o4wYMaIi+xf1Amg/EswQhKWTFpmKOOi1hV7GAQetPljgqJXBor5FG1p9YAaQGUAVZSQBDCvYah1CCE7p6quvllatWrl9fzfddFPO+skAEZZOhqBTIchgoZdxsG/tWFjgqJXBor5FG1p9IAEkAcw5wCZWIAEMK9hqHUKhnVKLFi3k1ltvlfr168sxxxyT17YABoiwdLLQOuXR0NqGtr6FXmpliMNYWOAYBxws+kACSAJIAvgXAlrnqq1v4di0Mlg4FY0MQ4YMkTZt2sjbb78tF154YV66aYFjoXGIE2khlitH00IvNbYVik5p9cECR60MFvUt2tDqAwkgCWBeQdZXYgYwrGyL1iEU0ik1bdpURo8eLXXq1JF+/frJe++9l5duMkCEpZOF1KlEBdLahra+hV5qZYjDWFjgGAccLPpAAkgCmFeQJQGsCpvWOVs4Nq0MFk4lXxnOP/98Ofjgg+Wtt96SkSNHSr4XIFvgWEgc4kZaiCUzgMneMl8fYZVJpU6uHBESQBJAEsC/ENA4JQuHYkFctH2w6Ec+MuB7v/fee6/b+3f66ae7K4VIAEXywdKSQFroZKF0ypJ0WPTBAkutPlj0o9AyWOAYBxws+kACSAJIAkgCWEkHtA4+n/qnnXaaO/TxySefyNlnn60iPgwQq4Yzn7FIVAZiGRaW2vG0IA2FloE6aaeTJIAkgCSAJIAFJYDrrLOOu/x6zTXXrNj7pwkyDBB2AYJYhoWlxi58T7RtaOtrSSh10k4nSQBJAEkASQALSgA7d+4sPXr0kO+++05OOeUUJ4smyDBA2AUIYhkWlhq7IAGsHOq0WGrra/2cRX0SQBJAEkASwIIRwHr16sn9998v66+/vlx++eXy/PPPkwAmjIY2yGjrkwCSACYHCK1OaYkLddJOJ0kASQBJAEkAC0YADzjgABkwYID8+uuv0qlTJ1m2bBkJIAlgFZ9UaNKhJS2ob0Fc4oCDFksLHLUyWNS3aEOrDySAJIAkgCSABSOAo0aNkmbNmskdd9zhTgFbLFMxQNhlCIhlWFhqA34IpEMrA3XSTidJAEkASQBJAAtCAHfaaScZPny4LF682J0AXrBgAQlgkjVqA762PoOtXbC1wFI7nlryZVFf24YFjloZLOpbtKHVBxJAEsC8CGCvXr2kd+/eUrduXSktLZWWLVu6u9vyLSUlJXl999W/T1sf7Wjb0NaHYysrK3MnYpcsWZIXlFoZVicOF1xwgYAEPvPMM3LnnXdW6q+mHxY4rk4cMg20BgeLPhDLVaOjHQsLLLUyWOhEoWWwwDEOOFj0obi4WCZNmiRFRUVSXl6eV8yp6ZXqiMiKmt6JQsnPT8ExQ5Cse1FmpZtuuqncfffdsmLFCjnxxBNl+vTpZllIZgjC0skQMh0hyGChl1FsK1ss0Lahra8dCwsctTJY1LdoQzsWzAAyA5jNX2T8nQQwrGCrdQiryyn16dNHOnToIG+88YYMGjSoio5p+sEAEZZOri6dyubINDpl0QcLvdT2waIfhZbBAsc44GDRBxJAEsBsfpMEMCJCWsdo4di0Mlg4lWwyrLvuum6ZG/0988wz5bPPPiMBTKNj2bDMppra+hY6uTp0KhsOIchggaV2PEPAQSuDBY5aGSzqW7Sh1QcSQBLAKL4z7TPMAIaVbdE6hNXhlLp27Sr48+WXXwr2kqYqmn4wQISlk6tDp6I4MY1OWfTBQi+1fbDoR6FlsMAxDjhY9IEEkAQwiu8kAYyAktYxWjg2rQwWTiWTDGussYa7+Hm99daTYcOGySuvvEICmEG3tOOprW+hk9WtUxFM0z2ixUJb3wJLrQwh4KCVwQJHrQwW9S3a0OoDCSAJYFT/mfI5ZgDDyrZoHUJ1O6VDDz1Uzj33XJk9e7Ycd9xxsnz5chJAEsCsPih0vc7aAV4EXQkizXiSANrFHBJAEsAovosZwAgoaZwamrdwbFoZqpsA4uQvTgDffPPN8uCDD6ZFVdMPCxyrG4cI6hSbrBWxXDnaFnqpsQuvc9o2tPW1+mCBo1YGi/oWbWjHggSQBDBqLGIGMAtSWmO0cGxaGarTKe2+++5yxRVXyO+//+4ufl60aBEJYDXrlFYfLHSyOnUqF+elxUJb3wJLrQxxGAsLHOOAg0UfSABJAHPxoVWe5RKwXTrewrGFGiDWWWcdl/Vr2rSpPPDAA3LLLbdk1DtNPyxwtHCumj6Ekq0hlrTvZEMttF5TJ+10kgQwxgRw7733lvPPP1922WUXadKkiRx++OHy6KOPVmjPXXfdJd26datk3/gqw8EHHxyZFJIA2hmjhWMrtHNORZzq1KnjMn+77bab2/t36qmnyvz580kAI1iZdjy19S10kmR65UBbYKkdzziMhQWOccDBog8kgDEmgO3bt5c999xTPvjgA5kwYUJKAojP+nTv3r0iFOG7rL/99luE0LTyERJAEsBsGQIQvk6dOrlv/p5++uny7bffZtUvTaBjgAhLJy0ClUYf4pRNjQMOWn2gfdvZNwlgjAlgYpTFJ7dSZQBxHccRRxyRNSCne4AE0M4YLRxbaAGiXbt2gm/+ogwdOlReffXVSLqm6YcFjtogZVHfog0NjlZZqxD6EYIMFnqpHc8QcNDKYIGjVgaL+hZtaPWBBLCWE0CQwj///FPmzp0rL730kgwcOFDmzJmTNkjjHrcGDRpU/A4COG3aNGncuHHGetmivlaRtfVDMEY4ti5dusiYMWNkyZIl2SBL+XtIOLRo0UKuu+46t/R13333CbYcRC2afljgGII+hCADsbSd4MXJvqPacqrnaN8rUdHgYFEfBHDmzJlSVFQk5eXlmiGtsXXriMiKGit9RMFTZQCPPfZYdyLz+++/l2bNmslll10mCxYskFatWqW9n23w4MEyZMiQKm8955xzZOHChRGlqfoYPg82b968gtXHiwstQ7169aR169by1ltvybJly/LCQtsHKxyWLl0qJ5xwguDwB5Z8H3nkkZz6o+mHBY5WOGh0OgQZiOUqtdXoJFqxwFIrQwg6pZXBAketDBb1LdrQ6kPDhg3dJJ0EMKfwVPMeTkUAk3uxxRZbyHfffSf777+/ywamKswAph977WzOItuilcFiVgkZsK90n332kSlTpsgZZ5zhJhq5FE0/LHC0wmHGjBm5dLvKsxocLPpALFcNiXYsLLDUymChE4WWwQLHOOBg0QdmAGvxEnCqyIRTmlgGvv322yMFLu4BtA0QOJWNy5Jr+hIw9pUeffTR8uOPP8opp5ySc0ZTE2S4RygsnbQIVBp98Gho29DWt9BLrQxxGAsLHOOAg0UfuAeQBLAiWmy88cby008/ucMijz/+OAlgJATCCrahBAjsJxk7dqxbVr/ttttk3LhxOaGp6QcDRFg6aRGoNPpAAljZ9LRYautr9YH2bWffJIAxJoBY32/evLnTlokTJwr26b388svusAb+YD/fQw895DaBYg/gVVdd5a512W677dzBkCiFGUA7Y7RwbIV2zonO/aCDDpJ+/frJH3/84Q63/PLLL1FUyj2j6YcFjloZLOpbtKHBEe8nlrTvZKPV6pRWr6mTdjpJAhhjAti2bVt55ZVXqgRdLDGedtppbmP+TjvtJLgKZvr06fLcc8/JRRdd5C7rjVpIAO2M0cKxFdo5Jzv3G2+8Ubbddlt3/QuugYlaNP2wwFEbpCzqW7ShwZEE0DZzZqGX2vEMQae0MljgqJXBor5FG1p9IAGMMQGMGmw1z5EAkgBmyhDgYNGoUaOkbt267qs0uJQ8StE4NgaIsHQyhEAXggwWeqmxC68V2ja09bVjYYGjVgaL+hZtaMeCBJAEMEo8TvsMCWBYwVbrEKrDKfXv318OPPBAmTp1qlsKjlI0/WCACEsnq0OnouhQpolJIepb6KXGLkgAbTO6cRgLEkASwHx8YUUdEsCwgm1ITmnHHXeUzp07u28Ao+AqmKOOOsp9Ei5b0fTDItCStKwcIWJJ+7Ym0lrbok7a6SQJIAlgtlic8XcSQDtjtHBsGuKkzRBssMEGAtLXpk0b2XLLLd1BDhTcQYnDR/fcc4/88MMPkfRN0w8LHLVByqK+RRsaHEkAbTNGFnqpHc8QdEorgwWOWhks6lu0odUHEkASwEgBOd1DJIC1mwDic2+HHXaYI35NmjRxYOCycJwix5+nn37aXQGDk+a5FI1jY4AISydDCHQhyGChlxq70E7wrOprx8ICR60MFvUt2tDqAwkgCWAucbnKsySAYQVbrUPIxSkdeeSR7jQ5Ps2EgkzfV1995T4tiFO/n376qSxatCgv/dL0gwEiLJ3MRafSKYtGH6yIi1YGC73UyhCHsbDAMQ44WPSBBJAEMK8A7SuRAIYVbFdHgFhrrbXc/X5Y6kV5/fXX5cknn6wgfKtDhkxKywARlk5aBKpC65RFHyz0Mg44aLG0wFErg0V9iza0+kACSAJIAvgXAlpj0ta3cGxaGbI5JVzrgvv8mjZtKkuXLpVbbrlFJkyYUEmHqluGbAprgWM2HLLJYFHfog3tWBDLsMi0djxD0CmtDNRJO50kASQBjBLL0j7DDKCdMVo4tuoMEPiyB74mgz1+uCx8yJAh8uWXX1bRjeqUIYqyWuCoDVIW9S3a0I4Fsaw99h3Ftix0UtsGddJOJ0kASQCj2n2l53r16iW9e/d2F/yWlpZKy5YtZe7cuXm1hUolJSUya9asgtUPQQY4trKyMhk/frwsWbIkLyy0OKbDoVWrVtKnTx8nEz4reMMNN8iCBQtSylhdMkQFxALHEPQhBBmI5Sqt0+q1BZZaGULQKa0MFjhqZbCob9GGVh+Ki4tl0qRJUlRUJPh+e20sdbB/vTZ23KLPzADazcYsZrbajE+62fmxxx4rPXv2lPfff1/69u2bUXWqS4ao+mqBozZLYVHfog3tWBDL2mHfUW3LQie1bVAn7XSSGUBmAHOx/SrPkgDaGaOFY9MG/HTOuWPHjnLWWWe5b0sPGzaMBDCC1VTXWER4dcUjWhksdFIb8C3qW7QRApZaGULAQSsDddIu5pAAkgDmEk9IADOgpXXOFo5NK0M659y+fXuX+Xv33XdlwIABJIARrKa6xiLCq0kAU4CkHQ9t/ZDte3XqFAmgHYHT6iQJIAlgLrZPAlhLCeA+++wjgwYNcvv//F7AdFBonRIDRDgBwoK0aMfTor5FG1q9tsBSK0MIOGhlsMBRK4NFfYs2tPpAAkgCSAL4FwJaY9LWt3BsWhnSOaU99thDLrvsMnfqFweAMpXqkiGqolrgGIJzDkEGYmlHxi2wLLRtWeiktg0LHLUyWNS3aEOrDySAJIBR42rK57gHsHYECHzq7dprr3Xf8j3ppJNIACNYjdY5hxAgGGxrh31HUOeKRwqt19RJO50kASQBzMX2qzxLAmhnjBaOrbqc8zbbbCMjRoyQGTNmyPHHH08CGMFqqmssIrzaLFhb6GQIRDYEGSywjINOacfCAketDBb1LdrQ6gMJIAlgLvGEBDADWlpjtHBsWhnSOaV//etf7hLojz/+2P03U6kuGaIqqgWOITjnEGQglrVjghfVtix0UtsGddJOJ0kASQBzsX0SwFpKAAcPHixt27aVu+66S8aOHUsCGMFqCk2EtYEW9Rls7YKtBZZx0CmtXlrgqJXBor5FG1p9IAEkAYwQytI/wiXg2hEgHn30UcFYn3766fL555+TAEawGq1zDiFAMNjWDvuOoM4VjxRar6mTdjpJAkgCmIvtMwNYCzOAzZs3l9tvv10WLVokhx12mCxbtowEMILVFDpQWhBIBlu7YGuBZRx0SquXFjhqZbCob9GGVh9IAEkAI4QyZgCjgKQ1RgvHppUhlVPyn4F755135IILLsgKRXXIkPWlCQ9Y4BiCcw5BBmJJAphse7TvlYgUGgcLGUgASQBzia3MANbCDOCVV14pu+22m9x8883y4IMPZtWXQjtGkpawSItFoCq0Tln0wUIv44CDFksLHLUyWNS3aEOrDySAJIBZA3qmB7gHMKxgq3UIyU6pQYMG8thjj7nDAN27d5cff/wxq75Yy5D1hUkPMECEpZMhBLoQZLDQy0LblgWO2jYscNTKYFHfog2tPpAAkgDmGl/d8/gaRO/evaVu3bpSWloqLVu2lLlz5+bVFiqVlJTIrFmzClY/BBng2MrKymT8+PGyZMmSvLDQ4piMw6677ir9+vWTn3/+Oes79FtpAAAgAElEQVQXQLzA1jLkCoQFjiHoQwgyEMtV2qfVawsstTKEoFNaGSxw1MpgUd+iDa0+FBcXy6RJk6SoqEjKy8tzdbWxeL6OiKyIRU8K0AlmAMPKtmhnhMmzUnz3t0OHDvLII4/IjTfeGEnDrGWI9NKEh5ghCEsnQ8h0hCCDhV4W2rYscNS2YYGjVgaL+hZtaPWBGUBmAHONr5WeJwEMK9hqHUKyU0I2slGjRi4L+N5770XSFWsZIr2UBDAlTNqxYLCNt33nalsWpEXbBnXSTidJAEkA8/EBFXVIAO2M0cKxaQN+onNu1qyZjBo1ShYvXuyuf4m6LG0pQz7KaYGjNkhZ1LdoQzsWxDK+9p2PbVnopLYN6qSdTpIAkgDm6wdcPRJAO2O0cGzagJ/onE844QQ56aST5M0335SLLroosp5YyhD5pcwAMgOYRVm0eqmtH5p952NbWvLm36nB0gJHi35o+mCBg0UfSABJAPP1AySASchpHYKFY9PKkOhURowYIdtss40MHz5cnnzyych6YilD5JeSAJIAkgBGMhetfWrra4mLhZ/UymBR36IN7ViQAJIARnIa6R5iBjCeGcDff/9dHn74YalTp44cc8wx8uuvv0bWE61T0jpGBoiwdFI7nhb1LdrQ6rWFXmplCAEHrQwWOGplsKhv0YZWH0gASQAjB/ZUD5IAhhVstQ7BO6XttttO+vfvL5MnT5aePXvmpCNWMsyYMSOn9/qHGSDC0skQAl0IMljoZaFtywJHbRsWOGplsKhv0YZWH0gASQDzCrK+EglgWMFW6xC8U8J9hB07dpT777/ffQc4l2IlAwmg/nNT2rFgsI2nfedrWxakRdsGddJOJ0kASQBzie1VniUBtDNGC8emDfjeOYP8gQSOGzdObrvttpx0xEqGfIOUBY7aIGVR36IN7VgQy3jad762ZaGT2jaok3Y6SQJIAphTcE9+mATQzhgtHJs24HvnvO+++8opp5wiTz/9tFx99dU56YiVDPkGKQsctUHKor5FG9qxIJbxtO98bctCJ7VtUCftdJIEkAQwp+BOApgerhCCrVYG75x33313Oeuss3K+Akbr3D26mn4wQNgFCGIZFpYau7CwLdq3nT5YYKnVBxJAEkASwL8Q0BqTtr5FsNXKACg222wzGTJkiPtvLp+AswowWsdogaNWBov6Fm1o9YFY2gV8Cyy14xmCTmllsMBRK4NFfYs2tPpAAkgCSAJIAlhJB3DpM5aAcfXLySefLPPnz89JR7ROSesYGSDCIi3a8bSob9GGVq8t9FIrQwg4aGWwwFErg0V9iza0+kACSAKYU3BPfph7AMMKtrk6hLXWWkvatm0re+21l+DvcK4777yz/Pnnn3L++efLBx98kLN+5CpDqhdo2mCACEsnQwh0IchgoZcau/BaoW1DW187FhY4amWwqG/RhnYsSABJAHMO8KjQq1cv6d27t9StW1dKS0ulZcuWMnfu3LzaQqWSkhKZNWtWweqHIAMcG07ejh8/PvJ3d5MBi4pjixYtpF27dtK6dWtp0KBBpWbWWGMNdwn06NGj8xqPqDJkalzThgWOIehDCDIQy1VaqtFJtGKBpVaGEHRKK4MFjloZLOpbtKHVh+LiYpk0aZIUFRVJeXl5Xv6+pleqIyIranonCiU/M4BhZVuyzQhB7s444ww59NBDKwSfOnWqPPPMM4KTgcuWLZO1115bnn32WVmxIj+zyCZDFF3VtMEMQVg6GUKmIwQZLPRSYxdeK7RtaOtrx8ICR60MFvUt2tCOBTOAzABGicdpnyEBDCvYZnIIm2yyiQwePFi23HJLR+6ee+45943fzz77rNL4ap2Ktr7WMTJAhKWT2vG0qG/RhlavLfRSK0MIOGhlsMBRK4NFfYs2tPpAAkgCSAL4FwJaY9LWt3Bs6WTYf//95dxzz5U111zTLdVfcskl8tFHH6Uce20/tPW1jtECR60MFvUt2tCOBbEMi0xrxzMEndLKQJ2000kSQBJAEsAYE0As+Z555plyyCGHuF6C9IH8ZdqvqQ0y2voMEHYOXjsWDLZ2Y2GBpXY8tbZlUV/bhgWOWhks6lu0odUHEkASQBLAmBJALPniPr8tttjCLfniUMfYsWOz7u3TOhVtfa1jZIAIi7Rox9OivkUbWr220EutDCHgoJXBAketDBb1LdrQ6gMJIAkgCWAMCSBO9w4cOLBiyffiiy+WiRMnRhprrVPR1tc6RgYIEsBUiq7VS219C73UyqC1LYv62jYscNTKYFHfog2tPpAAkgBGIgXpHuIhkLCCLRxC48aN5corr5T69etHWvJNHlutU9HW1zpGBoiwdFI7nhb1LdrQ6rWFXmplCAEHrQwWOGplsKhv0YZWH0gASQBJAGOUAcSFzgMGDHCXOr/66qsybNiwrEu+JICpTUDrXLX1QwgQDLZhkek46JRWr6mTdjpJAkgCSAIYEwLYtGlTGTVqlLvY+cMPP5T+/fvL0qVLcx5fbZDR1meAsHPw2rFgsLUbCwssteOptS2L+to2LHDUymBR36INrT6QAJIA5kwQEitwCTiMANGoUSMZMWKEgAR++umn0qdPH1m0aFFeY6t1Ktr6WsfIABGGTiYqn1YntPW1OmVR30Iv44CDFksLHLUyWNS3aEOrDySAJIB5kQRfiQSw8MEWX+4YOXKkbLbZZvLLL79Ijx49ZN68eXmPq9apaOtrHSMDROF1Mln5tDqhra/VKYv6FnoZBxy0WFrgqJXBor5FG1p9IAEkAcybKKAiCWDhg+3QoUNl7733duQPe/6Sv+yR6wBrnYq2vtYxMkAUXidJAKtanYVeFtq2tLbpUdH0wwJHi35o+mCBg0UfSABJAHPlB5WeJwEsbLA9+uijpVevXm6vH77xO3/+fPdNX03ROjZtfa1jY4AorE6m0j2tTmjra3XKor6FXsYBBy2WFjhqZbCob9GGVh9IAEkANVyBGcAE9LTGmKtj22abbeSGG26QevXquf8++uijopUhBKeklSFXHNMZgBZLbX0tDhb1iWVYZDoOOqXVS+qknU6SAJIA5kUAkXXq3bu31K1bV0pLS6Vly5YZPy+W7SUlJSUya9asbI+l/V1bHw1r29DWh2MrKyuT8ePHy5IlSzJiUVRUJFdffbXAgN988025/vrr3fNaGSzaKLQMueCYCWRtP7T1QxgLYrlKQ7TjaYGlVoYQdEorgwWOWhks6lu0odWH4uJimTRpkiCelJeX5x1/a3LFOiKyoiZ3oJCycwnYbjYWdWZbp04dd9HzrrvuKlOmTJF///vfFSd+mSEQiYpjNrvRYqmtr82UWNQnlqvfvjPpZRx0SquX1Ek7nWQGkBnAbHEw4+8kgHbGGNWxnXjiidK9e3dZvHixnHbaafLDDz9UCMEAQQKYaLBafYiqk9mciFYObX0t6bCob4FlHHDQYmmBo1YGi/oWbWj1gQSQBDCb7yYBjIiQ1hijOLZWrVrJJZdcIsgCXnHFFfLcc89Vkk4rQwhOSStDFByjDKkWS219LQ4W9Ynl6p/gMQOY2Tqpk3Y6SQJIAhglFqZ9hhlAO2PM5NhA+Lp06eL+4O9PPfWUXHPNNVXGJQ6kQ0tcGCBWj07m4ji0eqmtr9Upi/oWehkHHLRYWuColcGivkUbWn0gASQBzMWPV3mWBLD6g+26664rAwcOlF122cW97IknnpAbb7wx5WfetA4hBKeklYEBovp1MlenodVLbX2tTlnUt9DLOOCgxdICR60MFvUt2tDqAwkgCWCuvrzS8ySA1Rtscbp68ODBgk+9Yc/ftddeK88//3zaMdM6hBCcklYGBojq1cl8HIZWL7X1tTplUd9CL+OAgxZLCxy1MljUt2hDqw8kgCSA+fjzijokgNUXbNu1ayf9+vVz9/zhtO+gQYPkxx9/zDheWocQglPSysAAUX06ma+z0Oqltr5WpyzqW+hlHHDQYmmBo1YGi/oWbWj1gQSQBDBfn+7qkQBWT7DFPr9x48YJln9ffvllt99v0aJFWcdK6xBCcEpaGRggqkcns91NmUk5tXqpra/VKYv6FnoZBxy0WFrgqJXBor5FG1p9IAEkAcxKKjI9QAJYPcH2wAMPlHPPPdddjn388cfL8uXLI42T1iGE4JS0MjBAVI9OkgBupPrMooVe0r55zVNiINDqAwlgjAng3nvvLeeff747PNCkSRM5/PDD3efCEsvQoUOlR48est5667kvSuBeuW+++SYS2WAGsDJMWmNMDBD/+c9/ZOONN5bbb79d7r///sjjoZVBS74s6mvbsAi0Whks6lu0odUHYhkWmdaOZwg6pZWBOmmnkySAMSaA7du3lz333FM++OADmTBhQhUC2LdvXxkwYIB07dpVvv/+e7n44otlu+22E3xjFgcOohRmAO2M0Ts2fJoHd/35gqtfpk6dGmU4+C1gYYbAMkPAYGtv33fffXfWTz2mM3YSQNq3pX2TAMaYACYqyooVK6oQwOnTp8vw4cPdHxR8DxBLjt26dXP7z6IUEkCbAIE9f7vvvruAlMMoMV6+YDx++umnKMNBAkgCWElPtISBBNDGvtGKBZba8dRm3yzqa9uwwFErg0V9iza0+kACWEsJ4BZbbCHfffed7LjjjvLxxx9XeLlXXnlFJk6cKGeffXZKwrHGGmtIgwYNKn4DAZw2bZo0btxY5syZE4mkpHpIq8ja+oU2xosuukjatGnjcJw5c6Yj4jj88cwzzzh8o5aajoPvp6YfCBDImo4ZMybvTEuh9cECB4s+EEtbAqjVS41dhKJTWr2kTtrpJAgg4g2SP+Xl5VHDTKyeqyMiq9Itseraqs4kZwDxSbG33nrLZYygAL4g84dnO3XqlBIJ3Ek3ZMiQKr+dc845snDhwrzRw2nXefPmFaw+XlxIGY499lhp2rSp6z/GIOqSbzJg2j4UGgffH00/cG1O69atnX4vW7asYDql6YMFDhZjSSxXqY92PC2w1MpgoROFlsECxzjgYNGHhg0bynXXXUcCmHeUqCEVrQggM4DpBzyf2TmcWVlZmYAAwhhx1ctRRx2Vd+YqHxmSe6RtQ1ufGQK7Gb52LJhtsRsLCyy146m1LYv62jYscNTKYFHfog2tPjADyCXgnJaAk8kC9wDmHyA23XRTufDCC2WrrbZyjUyePFl++OEHd+dfvlduaB1CCE5JKwP3COWvk8n2TSzDwpL2bbOXUutjLOpbtKHVB+4BrKUEEMqHQyAgG/i8GArI3OzZs3kIJM+sbi7GeMQRR0jPnj0FGVXsvcC3fV977TWHPU8JiuowC0lLWKQlhEAXggwWepmLj0nnxrRtaOtrx8ICR60MFvUt2tCOBQlgjAkglhSbN2/u/AAOdmCfHg4W4LAGPi2GE6f9+/evdA3M9ttvz2tgqpEAwuCA+a677ure8t5778mVV17pxsTCsWkdQghOSSuDBY5aGSzqW7Sh1QdiGRaZ1o5nCDqllYE6aaeTJIAxJoBt27YVnOpNLsgwde/e3f0zLoI+9dRT3UXQb7zxhvTq1cstRUYtXAKOboxrrbWWjB49Who1aiR//vmn3HLLLZUu5rZwbAwQXCJKtF2tPljopDbgW9S3aCMELLUyhICDVgbqZPSYky2OkwDGmABmG3yL30kAoxvjXnvtJcOGDZNff/3VfeYt+W4/C8fGAEECSAKY2rNpbUNbn/Yd3Vdmik0WOGpJqEV9iza0OkkCSAKo4oEkgNGdGu5WPOyww+Thhx+WESNGVMHdwrFpHUIITkkrgwWOWhks6lu0odUHYhndvrM5UgssteMZgk5pZbDAUSuDRX2LNrT6QAJIApjNb2X8nQQwe4CoW7eu+xbzVVdd5S56vuCCC+Sdd94hAcygWRrHxgCRXSejGj2xDAtLjV34nmjb0NbXEh/qpJ1OkgCSAEaNBSmfIwFMb4w77bST9O7dW3DdS/369d2DuJgYWUDc95dcLBxboZ2z1rlbBCkLHC36EYexIJZ2wdYCyzjolNa2LHDUymBR36INrT6QAJIAkgD+hYDWmJLr9+nTRzp06OBa/+OPP9yev6effrrSwY9E8C0cm7YPITglrQwWOGplsKhv0YZWH4glCWBygNDqlFavqZN2OkkCSAJIAlhNBPCAAw6QAQMGyI8//lhx6joT2BaOrdDOWevcmQGsrCHa8dTWt9BJC53Q9iMEGSywjAMO2rGwwFErg0V9iza0+kACSAJIAlhNBBAf2J4wYYLUqVPHfert559/zoi1hWPTOoQQnJJWBgsctTJY1LdoQ6sPxNIu22KBpXY8Q9AprQwWOGplsKhv0YZWH0gASQBJAKuJAKJZfOFj2223leuvv14ee+wxEsCI2qZxbAwQYZGWEAJdCDJY6KXGLiyy6xY4atuwwFErg0V9iza0+kACSAIYMSSnfoyHQDIH286dO0uPHj3k3XffdcvBmYqFY9M6hBCcklYGCxy1MljUt2hDqw/EMiwyrR3PEHRKKwN10k4nSQBJAPMigPhiCE644oqT0tJSadmypcydOzevtlCppKREZs2aVbD61SEDyDE+87bBBhvI0qVL5cQTT3T/TVfg2MrKymT8+PGyZMmSvLDQ4lgdOOTTEU0/LHCMAw4WfSCWq7RXo5NoxQJLrQwWOlFoGSxwjAMOFn0oLi6WSZMmCbYr4Zv0tbHUEZEVtbHjFn1mBjD9bKxVq1Zy6aWXugfw7eWuXbtmhNxiZssMAb8EkqhkWn2w0EltxseivkUbIWCplSEEHLQyUCeZAbTgLr4NEkAFmiSA6Y2xRYsWcuutt7oHcCIYdwBmKhaOjQGCBJAEMLWVaW1DW5/2bUNcLHDUklCL+hZtaHWSS8BcAlbQPxESwPRObZ111qk4+NG+fXv5888/SQAjapvGsTFA2ARatEIsw8JSYxe+J9o2tPW1xIc6aaeTJIAkgBFDcurHSAAzG+Pjjz8uDRs2lG7durmLoJkBjKZumiDDAGEXIIhlWFhq7IIEsLLv0WKpra8lwhb1SQBJAKNF5DRPkQBmDhC33XabbLXVVmm//5sIq0WwjYNT0jo2Cxy1MljUt2hDqw/EkgQw2fVrdUqr19RJO50kASQBJAH8CwGtY0tVf8iQIdKmTRu56aab3KXQzABGUzfNWDBA2AUIYhkWlhq7YAaQGcBk70sCSAIYLSIzA5gVp1TO+dRTT5VOnTrJQw89JCNHjiQBzIriygc0gY6kJSzSoh1Pi/oWbWh0Eu+30EutDCHgoJXBAketDBb1LdrQ6sP/t3ce0FZV1xpegkNFGJYAagSxC4pCfHbz9PnsUZ+VEjEadGB5ENGYaGJ5aExRnyWWWGMhiA1bRJMomigxtvjUWMAaFRQRlCIgqCC+8a24L+ceTtlnz3m5m33/NcYZ98Ldc525/jXXnP+eq4kAigCmDMmVH9MUcO1gO2jQoDBkyJAwbty4cP7554sAprQ2i2NTgBABrGRmFpvyCNYedmltg0c7WlsHDxyLgINHG0QARQBThmQRwHpAVXKM3P7BETA33HBDuOWWW0QA64H49d8tQUYBQgRQBLD6QLOMLQ/SYa1D49tvfIsAigCmDMkigPWAquRYr7/++rDRRhuFM888Mzz11FMigPVAFAFshlBrB2sFW79g64Gl1R6s5MtD3lqHB45WHTzkPeqw2oMIoAhgypAsAlgPqPLB2L59+/CnP/0prLjiinEd4PTp00UA64EoAigCWMFGrIEuD8HWg7gUAQdrX3jgaNXBQ96jDqs9iACKAKYMySKA9YAqH4wbbLBBuPHGG8P8+fPDAQccUE9ci8RLELI4NgWIfGWt8hDo8qCDh11axkViFdY6rPLWvvDA0aqDh7xHHda+EAEUAaxLTGo9oE0g1YPt7rvvHs4666zwyiuvhOHDh9fF2cOxWR1CHpySVQcPHK06eMh71GG1B2GZLzJt7c882JRVB9mkn02KAIoA1iUmIoDpICp3znvssUdc+zdhwoRw4okn1q3Ew7EpQPgct2ENUh7yHnVY7cHDJvPQjjzo4IGltT/zgINVBw8crTp4yHvUYbUHEUARwLrEpNIDQ4cODcOGDQvt2rULvXr1Cr179w6zZs3KVBdCa6+9dpg2bVqrybeEDt26dQuXXnppvAP4yCOPDIsXL67ZPhzbgAEDwpgxY8LChQszYWHFsSVwyNIQSzs8cCwCDh5tEJZLrNdik9TigaVVBw+baG0dPHAsAg4ebVhzzTVjgmK11VYLc+fOzeKql3uZFUIIXy33rWilBmgKeAnwld7G/vCHP4QOHTqEY445Jrz77rt1CSB3Bo8cOTIzAbS+EebhrdSqgzIEtW2yEVchLPOFpca3Mvyl49dqD8oAKgPYSDxY6lkRwNoBggxgnz59wquvvhpmzJhRE2t2DfPsSy+9FL788stM/dKpU6cwb968pWTJPn766adNnzlz5oTx48eHmTNnLvWs1alY5UUA/UiHtS9EAP36wgNLa39ax5aHvLUODxytOnjIe9RhtQcRQBHATEQjERIBrB0gjj766Dj9m6assMIKTVegffVVtqT0SiutFKec0xQIKZtTpk6d2uxxq1OxylsdowJEvkiLtT895D3qsNq1h11adcgDDlYdPHC06uAh71GH1R5EAEUA03CFqs+IANYOtiuvvHLYZZddAj/rFc4LZOPIn//857Bo0aJ6j1f8OwO6UlaPujt27BjIEPJz6623Dt27d4/rLiGBH330UVN9Vqdilbc6RgUIEcBKg8Nql1Z5D7u06mAdWx7y1jo8cLTq4CHvUYfVHkQARQAzEY1ESAQwX8E2rUNg8e9ll10WSeD7778fTjrppKZNPGnrqGY4VnmrY1SAyJdNWvvTQ96jDqtde9ilVYc84GDVwQNHqw4e8h51WO1BBFAEUATwawSsg8kq7+HYGtGha9eu4fLLL487sN95551w8sknx51gjdTREpkWq2P0wNGqg4e8Rx3WvhSW+SLT1v7Mg01ZdZBN+tmkCKAIoAhgGyWANHvdddeNmcDOnTuH119/PZxyyilhjTXWWGpdYCNG0tpBSgHCL0AIy3xh2dpjy0reEjQt7ZBN+tmkCKAIYCOxfalnNQXsNxg9HFsWx9qjR4+YCeQsqN/97ndh3LhxIoAhLPeZUI9g7WGTHnpksetyZ2WtwyrvgaVVhyL0hQeORcDBow0igCKAIoBtOAOYdD4HUJ9wwglxAwr3F5fvDG7ESFo7SClA5OulxCNQtbZNebTBwy6LgIMVSw8crTp4yHvUYbUHEUARwEZiuzKANdCyDkYPx5ZVhz333DOcccYZ4YUXXggXXXSRCKAygNHSPWwyD4EuDzp4YJl1fJe6LWsdVnlrX3jgaNXBQ96jDmtfiACKAIoAKgMYttlmm3DhhRfG20q4v1gZQE0BiwA2d43WYOtBXKw65IF0WHXwwNGqg4e8Rx1WexABFAEUARQBDBtuuGG44YYbwuzZs+MdzyKAIoAigCKAlYKDlXRYiY8I4JJesfaFCKAIoAigCGDc+XvPPfcEbiA56qijwpQpUzLbhdUpKUD4OXhrXyjY+vWFB5bW/rSOLQ95ax0eOFp18JD3qMNqDyKAIoCZAz2C2gVcjADBNXT33XdfvCnk2muvDXfccUdmu7A6JatjVIDIl01a+9ND3qMOq1172KVVhzzgYNXBA0erDh7yHnVY7UEEUAQwU6AfOnRonCps165d6NWrV+jdu3fTTRJZKuQwYq4ly1qs8nyvtQ6rPI6N3bhjxowJCxcuzASFRYdDDjkkDBo0KMybNy8cd9xxraJD0mhLOzxwzIM95EEHYblkGFpsklo8sLTqkAebsurggaNVBw95jzqs9sCNUBMmTIhHgHEJQFssK4QQvmqLDfdoszKA+cq2WN4IV1pppTB69Oh4ODTnAkJEsxSLDsn3WepQhiBfNpmHTEcedPCwS8u48BhbHjha6/DA0aqDh7xHHVZ7UAZQGcAsMb5JRgQwX8HW6hD23XffcNZZZ4WZM2c2ZQMbNRCrDlbHqACRL5u09qeHvEcdVrv2sEurDnnAwaqDB45WHTzkPeqw2oMIoAhgo/G92fMigPkKtlaHwFrAW2+9NU6H33777eG6665r2D6sOlgdowJEvmzS2p8e8h51WO3awy6tOuQBB6sOHjhadfCQ96jDag8igCKADQf4UgERwHwFW6tDoDUHHnhgOPnkk8OcOXPCYYcdFr788suGbMRDB0sdChD5ssk8BLo86OBhl5ZxkViFtQ6rvLUvPHC06uAh71GHtS9EAEUAGwru5Q+LAOYr2FodAq3p1q1buPTSS0Pnzp3D2WefHR5//PGGbMRDB0sdChD5ssk8BLo86OBhl5ZxIQLY3I1ZsbTK58EmRQBFABsK7iKA1eGyOoQ8BYj9998/rgF85plnwumnn96QjVhxsDpGDxytOnjIe9Rh7QthmS8ybe3PPNiUVQfZpJ9NigCKADYU3EUA2wYB5Hifm2++OR4M3b9//7gpJG1p7SClAOEXIIRlvrBs7bFlJW8eWUjZpJ9NigCKAKaN6xWf0xSw32D0cGyeAeKyyy4LW221Vbj66qvDnXfemdpOPHVI/aUlD3rg6BHoWhsHjzYIy+KO7yxjy8OmrHXIJv1sUgRQBDCrH4hyIoB+g9HDsXmSjoEDB4bjjz8+PPbYY+Hcc89NbSeeOqT+UhHAilBZ+8LDJq0B30Peo448YGnVIQ84WHWQTfrFHBFAEcAsMbZJRgTQbzB6ODbPALHNNtuECy+8MN4LfOSRR6a2E08dUn+pCKAIYB1jsdqlVT5v4zvL2LKSt+Q7LVh64OjRDksbPHDwaIMIoAhgVj+gDGAZclaH4OHYrDqUOpXVV1893HvvvbGVbApZsGBBKlvx1CHVF5Y95IGjh3NtbRw82iAsi/uCl2VsediUtQ7ZpJ9NigCKAGb1AyKAy4AAbr755oEP17O1b98+cFBz8u7AoigAACAASURBVGGjBr+X/uROR+7yrfQ31K30/+X1kdWdP39+fLZnz57xO4YPHx5eeeWVVLbS2sRHAcIvQAjLfGHZ2mPLSt48Ml+yST+bFAEUAUwV1Ks9pClgv8FY7tiSA5kb6SDu8/3iiy8aEVnq2Up1DBs2LLz66qup6m3tIKUA0XI2mcoAKjxktQmrvAdxsergYZdWHfKAg1UHDxytOnjIe9RhtQcRQBHATD596NChAVJAlqhXr16hd+/eYdasWZnqQoirx6ZNm9Zq8nnQAcc2YMCAMGbMmLBw4cLQt2/feC8v5Z133gnPPfdcWLx4cTyaJfnJ78mH/2PaliNbkv9DtlyGf1NK6yl9JnEKyd9nzJgRJk2alLpvrH1p7YtyHFMrXvagtR1WeSsOHvLCcolRWPvTA0urDh420do6eOBYBBw82rDmmmuGCRMmBGaO5s6dm9VVLtdyKxALl+sWtKLyygAuAd/6Nlbpzfakk04KBx10UJg9e3YYMmRI3fP4rDrk4a3UqoMyBC1rk1ncjdUurfJWm/KQ97DLIuBgxdIDR6sOHvIedVjtQRlAZQCz+PMmGRHAlg22OLurrroqbLzxxuGFF14IP/7xj2PmrlqxOoQ8OCWrDgoQLWuTWRyG1S6t8lab8pD3sMsi4GDF0gNHqw4e8h51WO1BBFAEMIs/FwGsgJp1MFZzbN27dw/XXXddWGWVVcJNN90Ub+gQAaxutgoQIoCVrMM6Pq3yHnZp1SEPpMOqgweOVh085D3qsNqDCKAIoAjg1whYB5NVvpZj69evX2DdJWv1OJLl888/r9hvVh3y4JSsOihAiACKAFZ361YfYZXX+PYbn9a+EAEUARQBzDkB7NChQ7yOrUePHnEa+Ec/+pEygDWsVgTQL8AIy3xhaQ34VvLlIW+tQzbpZ5MigCKAIoA5J4DnnHNO2HXXXcPHH38cjjvuuLghRFPAmgJOM3CthEHB1i/YemBp7U8r+fKQt9bhgaNVBw95jzqs9iACKAKYJo5UfUabQFo2QHAszAknnBAWLVoU2BFc7yw+q0PIg1Oy6qAA0bI2mcVhWO3SKm+1KQ95D7ssAg5WLD1wtOrgIe9Rh9UeRABFALP48yYZEcCWC7bf+ta3wsUXXxxv4rj00kvD2LFj6/aV1SHkwSlZdVCAaDmbrGuAVR6w2qVV3mpTHvIedlkEHKxYeuBo1cFD3qMOqz2IAIoAZvXpUU4EsGWCLQc6s/N3jTXWCOPGjQvnn39+qn6yOoQ8OCWrDgoQLWOTHE6etVjt0ipvtSkPeQ+7LAIOViw9cLTq4CHvUYfVHkQARQCz+nQRwDLkrIOx1LGx0WPvvfcOb731VvjBD36Q+no3qw55cEpWHRQgRAArOTXr2LDKe9ilVQfr2PKQt9bhgaNVBw95jzqs9iAC2MYJ4Nlnnx3YZFBaXnvttbD55punIobKAPoH29tvvz1eB7fyyitH8jdx4sRUfeHhUDzqsDolqw4KEP42OXLkyHg9YdZitQmrvNWmPOQ97LIIOFix9MDRqoOHvEcdVnsQARQBDJwxt+eeezb5djYccP9rmiIC6B9s33///XDqqacGfh511FFpuqHpGatDyINTsuqgAOFvkyKAIVjHloddWnWwji0PeWsdHjhadfCQ96jDag8igCKA4eCDDw5bb711Q0QjeVgE0D/YbrHFFqFv377hxhtvDKNHj26oX6wOIQ9OyaqDAoS/TYoAigAmVtXaPkbj2298iwCKAMZs0yeffBI+++yz8NRTT4XTTz89vPfeexWJx0orrRSnJksJ4JQpU8I666wTZs6c2RBZKX3Y6lSs8lbS4SGPYxs2bFi86YP7fo888sgwbdq0hjAtAg5WLMGRzOmoUaM0bfnNb4apU6c2ZEOlDwtLv2DrgaXGdwgeOFp9jIe8Rx1We4AAfvjhh2G11VYLc+fOzewnlmfBFUIIXy3PDbDovu+++4ZOnTqF119/PU5xsCawW7duYcsttwzz5s1bqupKawZ56Ic//GH49NNPM6vCrldIaNZiled7rXVY5du3bx+OOOKI0LVr10jAWQfYaLHqkAccrDqA48477xyefPLJ8OWXXzYKYdPzViyt8lYcPOSF5RLzsfanB5ZWHTxsorV18MCxCDh4tKFjx47h17/+tQhg5ihRMEEG96RJk8Ipp5wSpyDLizKA1Tvc+ja24447hmuuuSZ89NFH4cILL4zHvzRarDrk4a3UqoMyBEusxmoPwjJfWFr70zq2POStdcgm/WxSGcA2PgVciWD8/e9/D4888kg444wz6vIPrQH0GYzbbbdd+OUvfxnv+73rrrvCiBEjwuLFi+viX/6AAsS/pogGDx4ctG4tH+vWrAHfQ96jDuvY8rBLqw55wMGqgweOVh085D3qsNqD1gCKADbjEKSEJ0+eHI+GueKKK+oSEBFAOwGE/P3iF78IZFeZdu/fv39cj5mlWB1CHpySVQcFCLtNJjUIy3xhqfGtF7zSuGC1BxHANk4AmWq8//7747TvuuuuG372s58FriBjJ+rHH39cl4OIANoCBFhfcMEFMWvFmjXW/jH1nvXMNatDsJIvD3lrHSItNpssHfTCMl9YanyLAIoA1qUlDT3QpjeB3HbbbWHXXXcNnTt3jmvP/va3v4UzzzwzvP3226lAFAHMHiDYOc2aP3ZggfuvfvWruPPXMnWpAKEA4RkgRACzj+9yB+qBpca3xrfn+FYGsI1nAFOxvBoPiQBmCxDgduWVV4bu3bvHHdgnnXRSPPrFunZNAUIBwjNAeJAWa0bXQ96jDuvY8sDSqkMecLDq4IGjVQcPeY86rPYgAigCaOKAIoDZCOB5550Xdthhh3hMyaBBg2L21cOxWR1CHpySVQcPHK06eMh71GG1B2GZbXxXcqoeWFr7Mw82ZdXBA0erDh7yHnVY7UEEUARQBPBrBKyDKY08Z1iddtppYa+99orfykabe++9N/7u4djS6FCvw611WOWtjtEDR6sOHvIedVj7QliKAJb7C6tNWe1aNulnkyKAIoD1+EDNvysDmH4wcoMKu6vJ/HHEC5s/Hn744aYKPBxbaztnq3NPwLC0wwNHj3ZY2uCBg0cbhGX68V3PkXpgWQSbstqlB45WHTzkPeqw2oMIoAhgPb8lApgSoVqDkdtWmPbt3bt3+PzzzyMRfOaZZ5rV7OHYrA4hD07JqoMHjlYdPOQ96rDag7AUAVQGsHIAsI6tPIxvEUARwJT0pvJjygCmCxCXX355vF6PwkHb9913X3j11VfD7NmzlQGsYFoW5yrSks4m0wx8YZkvLC3jIi9ZZSvxkU362aQIoAhgmjhQ9RkRwHSDkfsW+/btuxSO06ZNi0Rw4sSJ4a233go77bRTuP7663UO4De/GaZOnZrJNhUg0tlkGnCFZb6wFAH0WSttJaEe8h51WO1BBFAEME0cEAFMgVKtwbjCCiuEDTfcMGy++ebxkG1+rr/++oH/Twq/r7322vFMwAkTJkRSCDl8//33U3z7vx6xOgSPOlpbB5GWfJGWItiURxs87LK1x5YHDtY6PHC06uAh71GH1R5EAEUAU5OLSg8qA5g92Hbo0CH07NmziRCyPhBySOaLMwGTwvVwr732Wnj55ZfD3XffHebPn1+1z6wOIQ9OyaqDAkR2myw3LGGZLyw1vpUBLB2jVnsQARQBzEQAhw4dGoYNGxbatWsXevXqFTc3zJo1K1NdCJH5Yjo0a7HK50EHgu0xxxwTM38bbbRR2HTTTeNP/j8pt956a9OxMZWwKgIO1r4ArwEDBoQxY8Zknkq36uAh71GH1R6E5ZJRlgcsrTrkwaasOsgm/WxyzTXXjLNN3EY1d+7crOF3uZZr01fBWXtOGcCWzRBwbiBTxwMHDgx77LFHGD9+fLyvuVqxvhFas28e8tY6lLVqWZvM4jOsdmmVt9qUh7yHXRYBByuWHjhadfCQ96jDag/KACoDmMWfN8mIAC6bYLvbbruFESNGxC97/PHHw7nnnhtvESkvVoeQB6dk1UEBYtnYZCOOw2qXVnmrTXnIe9hlEXCwYumBo1UHD3mPOqz2IAIoAtiIH1/qWRHAZRNs11hjjXDVVVeFddZZJ34hO4WZDhYBXNp8FSCWjU024jisgcoqn4dg62GXRcDB2hceOFp18JD3qMNqDyKAIoCN+HERwBpoWQdjGse29957h5/+9Kdh0aJF4dlnn43rLjlX8K9//WvUzKqDRx2trUMaHNMYvbUdVvk89IWwzBeZLoJNWe1aNulnkyKAIoBpYmHVZ5QB9BuMaR0bawB32WWXpi9mx3C/fv0iGVSA0C7B0sFqtYe0NlnPiVj1sMpbSYeHvAeWRcDBiqUHjlYdPOQ96rDagwigCGA9313z7yKAy54Acl7gtttuG7p27Ro3h6y33nrhkksuCQ888IAIYBABFAGs7LKswdIq70FcrDrkgXRYdfDA0aqDh7xHHVZ7EAEUARQB/BoB62CyymdxbIcffng49thj43TwT37yExFAEcBm47k1bLKSQ7HqYZXPQ7DNMr7LsSwCDta+8MDRqoOHvEcdVnsQARQBFAFcTgkgR8R8//vfD9/73vfimsD9998/dOnSJfMVaokhWJ2KVd7qGBUgln1Wup4TsdqEVd5qUx7yHnZZBBysWHrgaNXBQ96jDqs9iACKANbz3ZoCTomQdTA24tg4FobMH99Jef311+PB3BwUm/UOXRHA5h1t7U+rfB4CRCM2WWuYWLGwyhcFyyLgYO0L2aTfC54IoAhgSnpT+TGtAfQbjGkcW9++fcPxxx8fb1+hsPFj5MiR4Y9//GM8F1ABQmsAS0eq1R7S2GQaB2LVwypvJR0e8h5YFgEHK5YeOFp18JD3qMNqDyKAIoBp/HfVZ0QAlw0B5DaQ4447Luywww7xCxcsWBBuv/32cOedd4bPPvusSQmrQ8iDU7LqoACxbGyyEcdhtUurvNWmPOQ97LIIOFix9MDRqoOHvEcdVnsQARQBbMSPL/WsCGDLB9t99903nHrqqYHdv2T57r///jBq1Kgwe/bspfrD6hDy4JSsOihAtLxNNuo0rHZplbfalIe8h10WAQcrlh44WnXwkPeow2oPIoAigI368mbPiwC2fLAdPHhwOOqoo+IXQfp++9vfhoceeigsXrxYBLCC9SpAtLxNNuo0rIHKKp+HYOthl0XAwdoXHjhadfCQ96jDag8igCKAjfpyEcAqiFkHYzXHRuZvr732Csccc0xYa6214rdPmjQp/P73v48/33vvvTBjxoz4/1YdPOpobR0UIEQAKw1Rq11a5T3s0qqDxrff2ChCX4gAigBmIoBDhw6Nu07btWsXNyT07t07bkjIWti9Om3atKzicferRZ4vttZhlSdADBgwIIwZMyYsXLhwKSz4+z777BNv/ejYsWOzv/P8J598Enhm/vz5MTvIDSH8ZNqYn3y++OKLMGfOnJhJpL+QmT59epgyZUpT/1nbYZW39kU9HNMambUdVnkrDh7ywnKJtVj70wNLqw4eNtHaOnjgWAQcPNqw5pprhgkTJoTVVlstzJ07N61rLNRzK4QQvipUi5ZhYzQF7PdGmTZDsOqqq4ZDDjkkku7u3buHddddNxJxykorrRRJXpYCcZw8eXIkgm+++WYkiRBEPsnvkMdK5LT0+1r7zTgtjvUwsrbDKp+HbI2wXPbju5ZdFsGmrHYtm/SzSWUAlQGsFwdr/l0E0G8wZnVsHAjNtXBrrLFGzGLOnDkzEsJKnw4dOgTe+pIPDqBbt27xw1RzGhLJDuRyYsi/ycC+8MILcVdya55FmBXHckO3BlurvDVQesgLy9Yf33l6ufKwKWsdskk/mxQBFAEUAfwaAWvAtsp7OLasOqy44oqRBPbo0SNw1iD/Xn311Zs+kEv+nWQaaxkNBPD555+PaxPff//9mFX86KOPwueffx6zk3xaMovogaM1SHnIe9SR1R6S/hWWfsHWA0trf+bBpqw6eOBo1cFD3qMOqz2IAIoAigCKADazgVpOhbWHCTFMSGHy74022iiSR7LCaaahIYE8B1F88MEHwyOPPBLXL1odowJEvkiLtT895D3qsAZbD7u06pAHHKw6eOBo1cFD3qMOqz2IAIoAigCKAKYmgPWMhenoXXfdNZLA9dZbL65R5Gfnzp3j+sRkmrlSPWQIH3300TB+/Pjw8ccfh3/+85/1vq7i3xUgRAArGYY1WFrlPezSqkMeSIdVBw8crTp4yHvUYbUHEUARwExBNhHSGsB8BVurQ2hppwRBhAgmH9YkbrfdduGAAw4I66+/fhOYOPnXXnstvPjii+Gll16Kn0oHX1cyXgWIfNlkS9tUWgdmHRtWeQ+7tOpQhL7wwLEIOHi0QQRQBDCt/6z4nAhgvoLt8hwgtthii8CtJ1tvvXXg6rvyaWTWE0IEufd44sSJVe1WASJfNukRqJZnu056w8Mui4CD1R48cLTq4CHvUYfVHkQARQBFAL9GwDqYrPIejs2qQx6cEjpABtnZ3KdPn7iuEEKYTB9zFM3BBx8sAphi5FrtwcMm82JTViys8h5YWnUoQl944FgEHDzaIAIoApgijFR/RBnAfGVbihogOnXqFK688sq4nvDxxx8PZ599tghgipFrtQcFW43vcjOz2pSVuMgm/WxSBFAEMEUYEQFMA5LVMXo4NqsOVufsIV+pjnXWWSfceuutsRvYIMJGkSeeeCKuEeTGk9LigaNHO4rQF8LSL9h6YFkEm7KOLQ8crTp4yHvUYbUHEUARwDTcpuozygAqQCyLDAGbR4YPHx7vRF5llVWafSVX3VG46o6jZfjJmYbsImZnMVPGHI7NJ7n+LjnIuvRn+ZpDq3O1yuchQCjYanwvi/HdSBCSTfrZpAigCGAjY2+pZ0UA/Qajh2MrAumoRXzYPTxixIiw8847V7Vb1gomOJRnB2sZOwSwlBBCJNl4khDIp59+Oh5onbYUoS88bDIPRDYPOnhgWQSbsvaFB45WHTzkPeqw2oMIoAhg2nhW8TkRQBHA1sgQrLXWWoE7kefNm9c0BcwtJdxgwv8PHDgwjB07Nm4c4aJzHB0frsDj36W3nPA7cuWl/F5lyCQ3nIwbNy5eecdUdK1idc55CBAKthrfrTG+a40r2aSfTYoAigCKAH6NgDVgW+U9HJtVhzyQDqsOWXDkPMKEFCYEcYMNNojTyfw/ZxSyI7m0TJ8+Pbz88stNn3fffbfZesQi9EUWLCs5FCsWVnmrTXnIe2BZBBysWHrgaNXBQ96jDqs9iACKAGYigEOHDg3Dhg2Ld8P26tUr9O7dO66vylrWXnvtMG3atKziwSrPF1vrsMrj2AYMGBDGjBlT967cakBZdcgDDlYdPHCspANZx9122y1su+22AXJYfqsJ19i98cYbAWJItpDs+Ny5c2NXVZqKLv8/1jJyjzLrFpOfENMPP/ww/h+fBQsWNP2e3K1ca9BY7aGlsGx0oFvbYbUpD3kPLIuAgxVLDxytOnjIe9RhtQdmRCZMmBBnRRJf1ejYXN6fXwH/vLw3orX01xTwEuStb2Meb7ZWHfLwVmrVwQPHejpAzDbffPOw5ZZbhq222iq+AJVvTimfQs4yRuvVAYmcMWNGJJ6vv/5604e1jPXakEafZYFlGj2KYNceWBYBB6tdeuBo1cFD3qMOqz0oA6gMYBr/W/UZEUARwHLjsDolq2NsjQBBJnyjjTaKZJC3aQpv15WuryvNHJb+zjpEiCVEMvl06dIlkBlM/p38nTbWKmQhIYRk1dm8AkFkvWSjpTWwrKRja9uU1SaR98CyCDhYsfTA0aqDh7xHHVZ7EAEUAWw0JjR7XgRQBFAEsPIQsjrnWgECwrnyyitHwrjuuuuGnj17Nn04LDsppRnEDz74oFmWEFLIlHKtomCr8a3xvezHd9qgbPUxIoAigGltreJzIoAKEAoQ+QoQ7ILedNNNIyFkvSIEkU+l8t577zUjhW+++WZcY5gUEUCNb43vfI3vUm1EAE30JQprDaABQxFABQgFiPwHCK7S22yzzZplCllAXl5YUzhp0qQmUvj222+Hb3/72+H666/PvDGpViYzreuxBro86OBBpouAg7UvPHC06uAh71GH1R6UARQBTOuDlQGsg5R1MHo4NqsOeXBKVh08cLTq4CHvUUcte+B4m9KpY37v3LlzMytnjSJE8amnngoTJ05sRgyTG1jSOBCrXVrlWxrLNBh42GURcLD2hQeOVh085D3qsNqDCKAIYBrfVfUZZQCVAVQGMP8ZwLSDnIBQSgrZ6cwxT1OnTm12lM2iRYviVXulO4/JHFYjhdZAZZXPQ7D1IC5FwMHaFx44WnXwkPeow2oPIoAigGljgzKAygCmshWrU7I6RgUI35eSE088Md6AsvHGGzeRQ178ygtX6b311lvNSOHkyZMjcbTahFXealMe8h52WQQcrFh64GjVwUPeow6rPYgAigCmCurVHlIG0DfYDh48OIwcOTLzeiurQ8iDU7LqoADR8ja5zjrrLDV9zOaT8sKB1mws4TDrZ599NkAIp0yZEj799NOG/E4R7NrDLouAg8a33/i02oMIoAhgQ464/GERQL/BrADhg6UHjtYg5SHvUYc1QDSCZbdu3ZqRQjadJIdjlx9oPWfOnEgEkw9H1CS/J4dYl/oaazuWNyyrOeUi4GDti0ZsslZws2Jplbfi4CEvAigCKAL4NQLWAW2V93BsVh08nEpr6+CBYxFw8GiDBUs2kHAmIWsKt99++8BVepBEgk6twpV6peSQ35le/sc//hFmzpyZ2V9Z7dIqb8EyabRVBw+baG0dPHAsAg4ebRABFAHM7FARVAbQJ2tFLR6OrbWds4dTstbhgaNVBw95jzqs9tASWJIV5FxCyGDySf4NSaxUkgwiZxRWyxx+9NFHFe9c9iJPecDSqkMebMqqQ0vYZJYgWIS+EAEUAcxi+00yIoAigOUG1NqOUQEiXzbZSMCH6LG+sJQc8vsmm2wSOLam9Oq8crtjZ/KsWbMCU8hcwcfP5MO/uWqPncvJ3+fOnRsWL16c2v9Z7drDLq06NNIX1YBpbR08cCwCDh5tEAEUAUztAEsfHDp0aBg2bFjgSiqOiejdu3d0vlkLZ41xb2nWYpXne611WOVxbAMGDAhjxozJvAnEqkMecLDq4IGjVQcPeY86rPaQJyw//vjjOI0MQeQDEeEnbeTTvn37mu6jfB0iu5OZbmY9Ih+IIaQw+Xfp7/wfmUtuTslaPLC09mcebMqqgweOVh085D3qsNoD95VPmDAh3l+OvbfFoptADL2uDGC+si2t/Xbu8VZqrUMZgnzZpLU/08iTGezatWsgoJEpTD5rrLFG/J2fTDNjG/y70jE29dxgOYFMnh8+fHh45ZVX6olriUcJQhY/pfHtN76VAVQGsK7jqvWACKDfYPRwbBbHmrTEWodVPk3Ar2WTHjhadfCQ96jD2hdFxZKZC3xXKUHk95NPPrmqaVUjgAiMGDEiZhCZBWG6ed68eUvV44GltT/zYFNWHTxwtOrgIe9Rh9UeRABFAEUAv0bAOpis8h6OzapDHpySVQcPHK06eMh71GG1h7aG5W677RbJXKVSiwCWP8+NKMlaQwghH6bY+vTpE8aOHRtmzJjRRBb5W9pzEa39mQebsurQ1myyVoC22oMIoAigCKAIYDMbsDoVq7wChF9W2doXCrYhri1k2phzDTmOhkxh+SeZeuZnx44dK/pUpqmT/mD9YWlhA0tCGMkiTp8+PYwePToeoF1arP1pHVse8tY6ZJN+/kEEUARQBFAE0DXItHaQUoDwCxDCsnEs2W2crDtMiCLEkGC71157xevykuln/l7pBhW+9e233w7HH398s/uVW3tsWclbgqalHbLJxm2yWpAXARQBFAEUARQBrDAKLEEqD4HSQwcFW79gWw1L/j9ZiwhR5HdOWGBn5j333BOeeeaZuDOZdYVkFyGQjRxfU27ay7tdyyb9bFIEUARQBFAEUARQBLCiH1Cw9Qu2jWC57777htNOO22pPknWIXJ8DWsKkw/kEJLIv/m99G+lv7PWkKNzpk6dmtnvWwmk9cWkERxrNdLaDqu8FQcPeRFAEcDMjgBB7QJunQBRrdOK4JSsjk0BIl82ae1PD3mPOqxjq1G7HDJkSNhiiy2ij+XTqVOnmB1kHWLWwtpDNqhwa0opUVywYEH47LPP4ofbVir95P8WLlwYunTpEuXL1zGW/jv5vdoznOeYnPta7ZlK9dFucDz00EPD3Xff3QyLevWU/50z9JI1lpX05XvYrAPZrlSs9pAHmxQBFAHM6kuinAhgvoJtEZyS1TE2GmhbikwXoS+EZb7Gd/fu3WN2DzKYEMNaP5PnmE5eeeWVY2Ma2c1caWxY5a061NpM00gwS9uOQYMGLbUZx+qjEj2tPsIqLwIoAtjImFnqWRHAfAUIq0PwcGytrYNIS75ssgg25dEGD7u0jC12M+OvuVaPrFYpceSWEwhivZ+Qpg4dOsQMISW5mq/8ir7Sf1d6Jk0dleqnLs5w7NGjR5g8eXJTFrKaHtXq4P9pL9lOZNGnWjnkkEPiDu3yYukLEUAT7XAVbvM3gXCt26mnnhrXhrz44ovhxBNPDM8++2wqkEUA8xVsi+CUrMHWI9BadfCQ96jDag/CUuM7b8TH2ya32mqrcNFFF8WpZYjto48+Gp5++unwxhtv1Lye1Dq28jC+lQFs4xlA7p4dNWpUOOGEE+JuM07C79+/f+jZs2dc51GviAAqQBQ9QNQbA9X+XoQA4R1sheXgMHLkyMx3fRfBpqzEx9smTzrppHDQQQc1M00O52a3demHe6BL1xEWoS9EANs4AeRNh2wfWb8kXY6hX3HFFeGCCy6o669FAEUARQArD5MiBAjvYFvXoVR5QFj+C5gi4GBth7dNsk5yv/32C5tuummcHmd6uXxKG53vc2f/VwAAG4FJREFUuOOOcO211zZZaBH6QgSwDRNABhJrQfr16xfuu+++JsPmDZXdZgcffPBS7pg1IMliYv4IAZwyZUocOJxgn7Uw/Vx+6n0jdVnl+S5rHVZ5+uPwww8Pt912W+YMgVWHPOBg1cEDR6sOHvIedVjtQVgu8UJ5wNKqQx5syqpDS9skMW799dcPG264YdNngw02CNdcc0147LHHmgyiCH3BuZNkOdkkxOaitlja7BpA3mA++OCDsNNOO8U1D0kh8/cf//EfYccdd1zKHs4+++xwzjnntEU7UZuFgBAQAkJACBQOAQjupEmTCteuNA0SAWyAAFbLAHbr1s30BvH3v/89bL/99mn6q+IzVnkqtdZhlU+yqRYsrTrkAQerDh44WnXwkPeow2oPwnKJu8kDllYd8mBTVh1kk/42qQxgZuqx/ApmmQIub22yBtBqQBMmTAi9e/fODKZVni+21mGV98DSqkMecLDq4IGjVQcPeY86rPYgLJe4pDxgadUhDzZl1UE2mS+bzBy0cyLYZjOA4M/UL2+Vw4cPj93B4lfOV/rNb37T0CYQKwHkKJqrrroqs0lY5fliax1WeQ/HZtUhDzhYdfDA0aqDh7xHHVZ7EJZLXFIesLTqkAebsuogm8yXTWYO2jkRbNMEkGNgfve734Xjjz8+EkGOgeH/evXqFaZPn163i7wGY90vagMPCEufThaOPjhSi7AUln4I+NQkm/TBUeP7Xzi2aQIIAMOGDWs6CPof//hHzAZCBtMU1gSefvrp4bzzzjPdUZnmu4r+jLD06WHh6IMjtQhLYemHgE9NskkfHDW+RQD9LEk1CQEhIASEgBAQAkJgOUKgzWcAl6O+kqpCQAgIASEgBISAEHBBQATQBUZVIgSEgBAQAkJACAiB5QcBEcDlp6+kqRAQAkJACAgBISAEXBAQAXSBUZUIASEgBISAEBACQmD5QUAEMGNfcSbVqaeeGu/QffHFF8OJJ54Ynn322Yy1FVOs0tV5r732Wth8881jg7lX+eKLLw7f/e534+8PPfRQPI+w9Aie9dZbL1x99dXhP//zP8O8efPisT3svP7yyy+LCVoIYZdddom2tc0224R111033ktdel81Df/Zz34Wjj322Hhv9RNPPBH++7//O95rmRTuubziiivCf/3Xf4XFixeHu+++O5x00knh008/bXpmq622CldeeWXYbrvtwkcffRSfv/DCCwuFaz0sb7rppjB48OBmbX7wwQfDd77zHWFZgspPf/rTcOihh8YjshYsWBCefPLJ8JOf/CS88cYbTU95jWeu4rzkkkvi4fjvvfde+MUvfhHHfVFKGiwfffTRsNtuuzVrMvfxMs6TksY3Fh3LE044IWLCdW4UDgs/99xzA2PYM8YUFUcRwAxehbMCR40aFTC+Z555Jp4f2L9//9CzZ88YSFX+hQAEsF+/fmHPPfdsgmTRokVhxowZ8d8cfr3//vvHAPzJJ5/EA7ghK//+7/8e/96uXbvA0TwffvhhJETc3wzuv/3tb8OZZ55ZWJj33Xff8O1vfzs899xz4d57712KAJ522mmRBH//+98P77zzTvj5z38eIHNbbLFF+PzzzyMuf/zjHyNenHHJrTcQHV5QjjjiiPh3zhMjeD/yyCPxGCPkb7zxxmjL4FuUUg9LcFl77bXD0Ucf3dRkMJw9e3bTv4VlCH/605/C7bffHm1oxRVXDL/61a/ClltuGW1u/vz5buOZQP7KK68EyM71118f9thjj3DppZdGPzFu3LhCmGUaLCGAjM8RI0Y0tRmc586dm9o3tgUsDzjggJgMePPNN+NFDvhEYsXWW28dJk6c6BJjioyjCGAGl8INIjhCsn4UDI83VTIoF1xwQYYaiykCASR7xWAsL9yeAlkeNGhQzE5RINBkCHfcccdIrAneDzzwQMyCJVlBCA0Yd+3aNSxcuLCYwJW06quvvlqKAH7wwQcxc8qHApbTpk2LRPqOO+6IWZpXX301bLvttpFEUvbZZ59ICrt37x6mTp0aX15++ctfxgx2giNEkP5KMrRFA7cSlhBAsqiHHHJIxeYKy8pW0KVLlzh+d9111/D4449HG/QYz+eff34ke7yQJOW2226LfVSalS2SbZZjSdsggLz8/vCHP6zY1DS+sS1iCVgkGCCBd911l2yyzkARAWzQk3jcIdzgVy63j0MAGYhk9z777LPw1FNPxcwVZJkp3b/85S/RsfP3pLz77rvxjZ8P05wHHnhgMwLJ2xhZL0glDrLopZy0bLjhhuHtt98O3/rWt+LSg6Q89thjEQ8yeGSzIIff+MY3mv7evn372Adkqn//+9/HKTWCdinxYcqJwMP0cWkGrCgYVyOAkN4vvvgizJo1K9rkWWedFWbOnBmbLSwr9/7GG28clxyQBWTazWs8jx8/Pjz//PPNiA8vNvgDfEURSzmWCQFkCpzkAjMg999/f8z0M/1OSeMb2xqWzBjh3/BtxAdebj1iTJFxFAFs0KMwrUYGZqeddop3CSeFrBTrBMheqfwLAd5SO3XqFF5//fU4HQkh7NatWwwarE0j+7LKKqs0g4vMHySEdTLXXnttWH/99WM9SenQoUOcciIbkKzzKDLe5aQFu2P9FXgSGJJC5o9nWU+ZTA+TvSotZAnpA6bXWG8JkSYTmBQyf0yb8JNMbNFKJQI4cODAaE9gQSBmapO1puDMcgRhubQVQErGjh0bCRlrLCmHH364y3jGV+AXyF4lhbFO9pqxz0tMkUolLGkf63snTZoUY02fPn3irAc3VB122GGx+Wl8Y1vBknhCcoFYwthlVolpdtlk/ZEiAlgfo2ZPiAA2CFjJ46uvvnp0aqecckp8kxUBrI+lCGB9jNI+UYkAlssmGVbWnpE9EAFcGl3W7kLKWKs7ZcoUEcC0BljhuUpYVqouybDyksIMgAjgEpSYlevRo0cgvrDmfMiQITEZwyyJR4wpMpEWAWxw8GoKuEHAyh7nLZaNBw8//LBLet6mTf6lNQXs10dpCCDfxnpTpoGvu+46TQGXwc8654MOOiiu/WO5RlI0Bdy4nVbDslJNq666atzBz1peNsNoCrg63sSWf/7zn3E9tKaAa9ulCGDj4zZO/UJkhg8fHqVJ40+ePDnuYtUmkOqAduzYMeJ0zjnnxHUaLBonTX/PPfdEoc022yxOF5dvAiHrmuyuZmqEo0rWWmutuG6r6KXaJpCLLrooHpVBYUcvpKV8EwjHyLCeirLXXnvFKfPyTSDsgGVnNoVNIRz10ZY2gZTbD0sUsFHWBbLuKtkEIixD3OTGmlHWipYeOQSGySYQ63hm6ne//faL055JueWWW+J61iJtAqmFZSWftvPOO8fjnsDl5ZdfbtogV8s3thUsy/H685//HMcwx155xJgi4ygCmIFBcAwMBIYdqRBBFt7zfwSL0jPsMlRdKBGIGkGUaV928vLWSlqeoyM+/vjjuEUfZw9xmTNnTgwwFI5AoSTHwLAOhqNPWNR78803x+MhinwMDER5k002iRgkOwFZF8nGBDbQgAVrJEuPgSEwlB8DA7ljjV9yDMz//d//NR0DQ8CGbJNN4KWFdTQcA8OuwyIdA1MLS/BkTSS70FlPyfTa//7v/0ZCzS7U5AWD9WdtHUvOi2RtFdk/7CYpyQYv/u0xnpMjN/g+7HH33XcPl19+eaGOgamH5UYbbRSxxu7Y0crY/vWvfx3ef//9prMB0/jGtoAla3ZZ7wfhY9yCG+dTkillpkk2WZtSiABmpFzDhg1rOgiaIE02EDKosgQBjm9gqqhz587xTexvf/tbJG6sYaEkB8eSNSg9CJrNCklhbQcHQZN1YAoE4g35KfJB0KxfYVdveRk5cmTTeXWQ6eOOOy4uxAdXDtDmLKyksJOXjHTpQdDYaLWDoCHkEHAIUJFKLSw5QJYd0ewYBEdeNCDE//M//9PsRU5YhrjBqFLh5S05pNlrPNNnEB5eaCA97H4t0kHQ9bAkSz969Oj4UsYLDC99nAfKgdjJOYD0RRrfWHQsk7MiyYTyMvLSSy/FF1rIn2eMKSqOIoBFinZqixAQAkJACAgBISAEUiAgApgCJD0iBISAEBACQkAICIEiISACWKTeVFuEgBAQAkJACAgBIZACARHAFCDpESEgBISAEBACQkAIFAkBEcAi9abaIgSEgBAQAkJACAiBFAiIAKYASY8IASEgBISAEBACQqBICIgAFqk31RYhIASEgBAQAkJACKRAQAQwBUh6RAgIARsC3KnL4coc1ErhLDQOBX/xxRdtFRdAmvPzODOT69TeeOON8KMf/Sj+e7311mu11nGLArehoNOyLtwxXKv9ffv2jQekcwOTij8Crdn3/q1RjbUQEAGUfQgBIVATAW4h2WmnneLNGIsXL44H0z700EOBK5I4QDpLEQFcgtoRRxwRzjjjjHggdV6uN0xDAg444IB4GH6iN/evcnvHtddem8UkUsvUI4DcfHPrrbeGbbfdNnAbBkT2vvvua1Y/7eNgdK6U5AB/DlVH/+W98NIwceLEZs1YZZVV4q0i3OKSpqTp+zT16Jn8IyACmP8+koZCoFURgAByY8Zll10W9eCu4BEjRgTuJ91uu+0yXX8oArikS88666zAXcPcs5umtG/fvsVvwqlHArhikLujeW7s2LHxKsd/+7d/i9c9cvtMS5Y0BJCbcZ577rlIBE888cRmBPC73/1uvPGG+7EhfVwHeNhhh8WbN3jBKVKBDHPDDWSXrGqaUq/v09TR2s8sizHS2m30+H4RQA8UVYcQKDAC5QSQpuJgmYYjs8Ddm1xZdcstt8RMIVOaTO0SeLmaiUKQZco3ITkJAZwyZUrMKEIq33333fgs8lOnTo0X3pdfr7j++uvH54488shYZ5cuXSI5PfbYY8OiRYuiPBmpiy++OEAUuO+Xq6G4MgqdZ82aFUkr99mSweKuar6HjCYE4PHHH49XF1YiAhAc7remHVxFePvtt8fASoFM8D3c4wqp4L5mLqWn3HTTTWHhwoXxrtL9998/BmTuER8/fny46KKLIk7c7frZZ59FvRO8uYKOwr/Bge/lnmwIDNf2kWk79NBD45Vpf/3rX+M9z1wXxjQ7Vy9yV/RTTz0V61hxxRXjFXc8w7V3TzzxRLwnGpwp1HHDDTeE3r17B+6MfvbZZ8P2229fcQq4U6dOgX6DZNHn1Uo1TGgH+nK/8oIFC6I493xzZ/iGG24YNt1009inSftXX331cN1114W999473tnM/a7cD5xmCvidd96Jd7WXZgDvuOOO8Nprr0X7SbDhikTqp08qFTDnCkr0w4aQ5Xo62kJfgB82wRVkP/jBD6LdUeiLc845J7YPknzNNdfEPqplp8nfaCf1khWmHw888MCGvUz//v3jd3IX++eff15Rvl7fY9cDBw4M3/jGN+JYpe133XVXrOv555+PL4alV/VxNy82W+laSe4a5wpG8Jg+fXq88o+7kSn0OaSdq/+SAgaMDWy1ng3jU8CeZ7EhfAN+odr38R277LJL/H4yxQ8//HDsW/zE0UcfHVVgPF966aVhxx13DPPnz4/3pHP/cLXr/BruoFYWEAFs5Q7Q1wuBvCNQiQCiM46aII9zhNx85zvfCQ888EAMhAQNiFWvXr1i86oRQIjinXfeGV555ZWYPaIQbMmKQcjKS0IAyexAovjeZ555JpIbghCkYsKECdHp33333ZFYcr8vhPEvf/lL+MMf/hBJH0HxkksuidkqniO4Q+Z23333OGVYXgigTK2xlpFpToISWTvuYd54443Dyy+/HMkV2TDkR40aFckUZBUCCFEjgEMwqWPIkCGRTFTChntHSwkQ+BOkIawQM6b0Xn311TB79uyoP8QCQrfSSivFdkB2yNBCtiHBFPoDfb/3ve+FGTNmxCAGweO7CHgQYggtfcBz4AR5r7QGkD6H+K+66qqR2FYq9TABr/POOy9m6Cise9xnn30iCStvP/0KYRgwYED8TjBmejcrAcTewA+MKBALgjs6QObKC7hzDzdkiju6IRbdunWLL0B9+vRpskEIEnWDJVPK6Eof7bHHHrHfIbIQEwh2PTulP+kDMCKLt8MOO8Q6Gi0PPvhgXFeavKiUy6fpe0gs5Ai7AQOwSF7Yhg0bFvr169dkJ4wT7nqHUEHWywvjgJcZ7njmfnfsaM899wxPPvlkzMJilz179oxitJsXFMYNf69lwzwPKeM59EFXbJMxUO37eBHiBeGUU06J4xXbu+eee+I4gAB26NAhjnkIIEQUG0RfSCs+oAhFBLAIvag2CIEWRKAaAeRNG+e52WabLfXtBDuCH4GSjFctAghR5C0c0kAhaBFwyIaVl4QAQiwJtBSyQ2TPCHI//vGP49Q0gSYpZFxw3pAu/k4GkIwGATzJwBFgIV2QRLJL5YWMHgSFYF5eWL/HxoX99tuv6U9kFCELBHAIIKTt8MMPj38nSJJBI9NIpqgcm0oEEF3JniSFwEXdtJ1CcCSbQdspBGhINd9LIJw7d27MHiYZWbKsZL0I1Hwg7hCbJItKwKOOSgQQQkDffPOb36xqdfUwYe0gZJuXBkrSF6NHj25GACHakDPaBvmlQATJ4mUlgIMHD44vL/TnW2+9FUkv/QshZFq7vBD0ycCVZqaqNZx1dmSJGRMQQDKx9BvTr/RBUurZKbaDfZ977rnRVrKUHj16RDJG5jrp9/J6sNtG+h75F154IbYR8g6JQr/kZYcXEOws7TKAe++9N/YrxI8XGAgf44iXOl6k+B5IM6WWDUMoIYCV1nuWtrn0+3gZYvaCzWlJAQv6DAIIkcSOwS8p+BBeUCGtRSgigEXoRbVBCLQgAmkygBANSAHOm0wIU6gEBzJQBJ9aBJBAPnny5JidevPNN+MUKkSPKaLykhBA6v7kk0/in3kj59847d/85jeR6CVTi/ydLAfZE6ZfyW4R0AlYZB9ZNE/QIetAsNx1110jcSovEFTqgfSWF8gShSnRpFx99dURAzIkEEDIcELgEnIM8WLaMw0BJOtFO5NSPrVZXkcpTmRSCGrgVTp1BQmE0IIB05RkGZNCYISYVyKAZErITtXKANbDBBJMG9CTjRhkvMiKQfZKCTB/mzZtWvxbYg+Q9aeffjoSQAgMU45JISNcWipNAfN3Ajv2Ql+Q0YKI0sdkhssLGWXIH5mh8sJLC3bPSwVT4xBWCDdjgALJJbMIKSLrSaYaclfPTpl6pE8hqUxL8nwyVZp2qCNP9hLdqhXIdL2+ZwqdMdW9e/doP7STNiUvShBByCpEmswq2EK0KhVeHpDF9sEKG6JdvEhSSscN2TheDCDBvJzUsmEIOrpB1iCoSan1fZB+spClL27YLZk/bAOSDjHlRSkp6Mw0eKXZibT9kqfnRADz1BvSRQjkEIFqawBxtJAp3vrPPPPM+NbPNA4ZgYTkJEe91CKANJkAS4AhkDClXGkalufqEUCIC9+ZZNvK4cSBMz3EuiWyNAQIpoAJzkwbQzQqFYIFxKfSm3+lbBekhHVuSQbQSgBLN+GgXyMEkClisieQ3yRrWtrGSlkggjKEsBIBTNYAglcyhVuOWT1MeJ6pebK9ZBIheqxZpJQSQPoLMo+OSQaQacgxY8ZkzgCW65qsbYMole+g5VlsnKnFZO1eqTyZaqZYaS8EmwwghDJZv5g8yzQzLwiQGf6GPdWy09LvgDyytpCfrLlLUyDHSZa41q7sen3Pd0L2IbKMd0gWP2ljsimMlwgy0bzAsfwAcp9kkkt15UWDlyxeLCDBLBWBKLJMInk5IoPNy06SUcQGaQftqWXDfE/5xrJ631cvA8gsAeSXdc1FLSKARe1ZtUsIOCFQTgDJlpHJYFou2QXMFCSZGd6mcdYsACfgpSWArIfjbZ8sDxkCgkClUo8AEnwIUJCTpA6CCVkw1l5RCFKQDKYhWdRNlo6pP6ZtcfqVClPZZDeQIaNHNrB8DSBTQwRL1iwxlUmWgODlkQG0EECICVPdTJuTwWS6DNIDmYVIQU4gMTfffHMkOWyigfRAiqudA0g9PMu0O2sqCc5keyE4kKBkDWA1TMCYAAyeZHeOOuqopk0z5VPg6MUz1EXGiP7D7mpNATOdyN8hvGRy0JHMHFlZXk4g+rQZ8knWiSwPazgrFdrDJiJebsgkl64BZKqS8cG6TmyELCH2BsmD1EJoIG/z5s0LTD3TD127do3fW8tOWbOKPTIeqI91bNRF1hp7oiQbFSrpnKxnYzzwAlCt1Ot7MqPYMuOYpRz0E2QPTBMCSN0QO+wMLJJsXvl3sqQAm8K+yLJDBNlMQn2lyxuYBWBqnhcNfExSatlwJQJY7/voI8Ynu55pI+MBQsryArBlYxv6Mg3Nmj/sZ5NNNol9V22zkJPLXWbViAAuM6j1RUJg+USg/BxAMnxkuCB9TMtQCKhkg9hYwNmAEEQCd1oCSB3smiXYEUjJDmQhgMjwnejGdBAZpGTBP+v7KLzVM52KM2e6GWLEM5CaWtkSSAdBiIX/7MxkXRd1UQhmZPuYuiMYkokkw0XJAwGEAJN1IsvGekiyoODN1B4FsgrJAX8ybZBl2lvrIGgyvsk5gOwwBUvqgFTXw4S/Q+bYKABxYL1aMj1dTgCZ3qdONp8wXZ9mFzCBnWnG0gIBY0MJLxEQXH5CXNnJTAav1hmMEDLwQ4Z1m4l9kyHDZvh/CCVEgrWEkAtwxkawR+yQv4MX2a96doqekDhIEFPgbFhKptXpN+oF62oFEkPmlDbXK7X6HhLNtDgvRvQxY5oNOBClUgIIHpB/xgakqVphmpgXQ16geEFjGQI2UEoAwY/n2EiT2BL11bPhSkdL1fs+poDZ+IOtMF4h6uCWLPVgPPMySz+zzAUbhxCCbxGKCGARelFtEAIFQIBjSAiuBEkVISAElkaAzCZraiFtlaZZWwszCDLZ4FrrDVtLt0a+lyUJLN1g7V9bKCKAbaGX1UYhkHMEeNNmSoxpVaZ/VISAEFg+EGCqlKwmU+nL2/EoZJVZV8kaXXb9kuEki8k5kW2hiAC2hV5WG4VAjhFgionNGFwt11bevHPcHVJNCKRGgHWcjF/WObJGstrSjdQVLuMH2cDGmkWWI7BsgB3RbAprK0UEsK30tNopBISAEBACQkAICIGvERABlCkIASEgBISAEBACQqCNISAC2MY6XM0VAkJACAgBISAEhIAIoGxACAgBISAEhIAQEAJtDAERwDbW4WquEBACQkAICAEhIAREAGUDQkAICAEhIASEgBBoYwiIALaxDldzhYAQEAJCQAgIASEgAigbEAJCQAgIASEgBIRAG0NABLCNdbiaKwSEgBAQAkJACAgBEUDZgBAQAkJACAgBISAE2hgCIoBtrMPVXCEgBISAEBACQkAIiADKBoSAEBACQkAICAEh0MYQ+H86sfV7nicOJgAAAABJRU5ErkJggg==\" width=\"640\">"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "MovieWriter ffmpeg unavailable; using Pillow instead.\n"
     ]
    }
   ],
   "source": [
    "import numpy as np\n",
    "from matplotlib import rc\n",
    "import matplotlib.pyplot as plt\n",
    "import matplotlib.animation as animation\n",
    "\n",
    "x = df1['MA7']\n",
    "y = df1['ppl_vac_percent']\n",
    "date = df1['Date']\n",
    "\n",
    "fig, ax = plt.subplots()\n",
    "plt.style.use('dark_background')\n",
    "plt.rcParams[\"font.family\"] = \"arial\"\n",
    "\n",
    "\n",
    "\n",
    "line, = ax.plot(x, y, color='white', alpha = 0.75)\n",
    "sctr  = ax.scatter([],[], s=50, color='green')\n",
    "time_text = ax.text(0.8, 0.94, '', transform=ax.transAxes, fontsize = 'large')\n",
    "\n",
    "\n",
    "\n",
    "def update(num, x, y, line):\n",
    "    line.set_data(x[:num], y[:num])\n",
    "    ax.set_xlim([0, 3000]) # fix the x axis\n",
    "    ax.set_ylim([0, 25]) # fix the y axis\n",
    "    sctr.set_offsets((x[num],y[num]))\n",
    "    time_text.set_text(date[num].strftime('%Y-%m-%d'))\n",
    "    return line, sctr\n",
    "\n",
    "ani = animation.FuncAnimation(fig, update, len(x), fargs=[x, y, line],\n",
    "                              interval=50, blit=True)\n",
    "\n",
    "# Grid and ticks\n",
    "x_major_ticks = np.arange(0, 3001, 500)\n",
    "x_minor_ticks = np.arange(0, 3001, 100)\n",
    "\n",
    "y_major_ticks = np.arange(0, 26, 5)\n",
    "y_minor_ticks = np.arange(0, 26, 1)\n",
    "\n",
    "ax.set_xticks(x_major_ticks)\n",
    "ax.set_xticks(x_minor_ticks, minor=True)\n",
    "ax.set_yticks(y_major_ticks)\n",
    "ax.set_yticks(y_minor_ticks, minor=True)\n",
    "\n",
    "ax.grid(which='minor', alpha=0.1)\n",
    "ax.grid(which='major', alpha=0.4)\n",
    "\n",
    "# Labels\n",
    "ax.set_xlabel('Daily new confirmed Covid-19 cases, 7 day average', fontsize = 9)\n",
    "ax.xaxis.set_label_coords(0.7, -0.075)\n",
    "ax.set_ylabel('Share of people who received at least one dose of vaccine, %', fontsize = 9, rotation = 'horizontal')\n",
    "ax.yaxis.set_label_coords(0.35, 1.02)\n",
    "\n",
    "# Title\n",
    "ax.set_title('Are Covid-19 vaccines working?', y = 1.07, fontsize = 16)\n",
    "\n",
    "# Animate\n",
    "ani = animation.FuncAnimation(fig, update, len(x), fargs=[x, y, line],\n",
    "                              interval=50, blit=True, repeat=False)\n",
    "\n",
    "ani.save('Final.gif')\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.6.4"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}