{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Mauna Loa Time Series Example" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Gaussian processes can be used to model nonlinear time series. We consider the problem of predicting future concentrations of CO2 in the atmosphere. The data are taken from the Mauna Loa observatory in Hawaii which records the monthly average atmospheric concentration of CO2 (in parts per million) between 1958 to 2015. The purposes of testing the predictive accuracy of the Gaussian process model, we fit the GP to the historical data from 1958 to 2004 and optimise the parameters using maximum likelihood estimation. \n", "\n", "We employ a seemingly complex kernel function to model these data which follows the kernel structure given in Rasmussen and Williams (2006). The kernel comprises of simpler kernels with each kernel term accounting for a different aspect in the variation of the data. For example, the `Periodic` kernel captures the seasonal effect of CO2 absorption from plants. A detailed description of each kernel contribution is given in Chapter 5 of Rasmussen and Williams (2006).\n" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "using GaussianProcesses\n", "using DelimitedFiles\n", "\n", "data = readdlm(\"data/CO2_data.csv\",',')\n", "\n", "year = data[:,1]; co2 = data[:,2];\n", "#Split the data into training and testing data\n", "xtrain = year[year.<2004]; ytrain = co2[year.<2004];\n", "xtest = year[year.>=2004]; ytest = co2[year.>=2004];\n", "\n", "#Kernel is represented as a sum of kernels\n", "kernel = SE(4.0,4.0) + Periodic(0.0,1.0,0.0)*SE(4.0,0.0) + RQ(0.0,0.0,-1.0) + SE(-2.0,-2.0);\n", "\n", "gp = GP(xtrain,ytrain,MeanZero(),kernel,-2.0) #Fit the GP\n", "\n", "optimize!(gp) #Estimate the parameters through maximum likelihood estimation\n", "\n", "μ, Σ = predict_y(gp,xtest);" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Over the ten year predidction horizon the GP is able to accurately capture both the trend and seasonal variations of the CO2 concentrations. Arguably, the GP prediction gradually begins to underestimate the CO2 concentration. The accuracy of the fit could be further improved by extending the kernel function to include additionally terms. Recent work on automatic structure discovery (Duvenaud et al., 2013) could be used to optimise the modelling process." ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAlgAAAGQCAYAAAByNR6YAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAAPYQAAD2EBqD+naQAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAIABJREFUeJzs3Xd4VGX2wPHvTQ8JCaTQQ4L0HjpIDYuAigIubUFpgmUtv1VWcRdERNRVbGtbyyqgIi7SF1FXERCQIiUUgdBSIBDSe8+8vz+GGVKmZ9LM+TwPj2bmzr3vnSubs+c973k1pZRCCCGEEEI4jUtND0AIIYQQ4vdGAiwhqkhMTAyapjF79uyaHkqttWvXLjRNY+nSpTU9lFotLCyMsLCwMq8tXboUTdPYtWtXlVxz9uzZaJpGTExMlZxfiN87CbCEsIGmaXb9EaK2W7VqFZqmsWrVqpoeihC/S241PQAh6oLnnnuuwmtvvfUWGRkZJt8DaNmyJWfOnMHf37+qh1dn9e/fnzNnzhAUFFTTQ6lzHn30UaZNm0br1q2r5Pwvv/wyzzzzDC1btqyS8wvxe6dJkbsQjgkLCyM2Nhb5KySqmmF60JnTdatWrWLOnDmsXLlSprGFqAIyRShEFTFXgzVixAg0TaOgoIC///3vtG7dGm9vb/r06cOPP/4IQEZGBo888ggtWrTAy8uLQYMGcejQIZPXSUxM5IknnqBdu3Z4enoSFBTEH//4R06dOmXzWDMyMliyZAldunTB19cXPz8/2rVrx6xZs4iNjS1zrFKKTz/9lMGDB+Pn50eDBg3o27cvn376aYXzlq4TWrVqFb1796ZBgwaMGDECsFyDZc99nT9/njlz5tCmTRs8PT0JCAigZ8+e/OUvf7EpAC49zk8++YTu3bvj5eVFy5YteeKJJ8jKyipzfOlne+bMGSZOnEhgYGCFmqXo6GjmzZtH69at8fT0pHnz5syePbvCd2qwZcsW+vXrh7e3N02bNmX+/PmkpaVZHXN5x48fZ8aMGbRq1cp43bFjx/Lf//4X0NdXzZkzB4A5c+aYnN62VIO1cuVKBgwYgK+vL76+vgwYMMDkVGPp53v48GFuu+02GjZsiL+/PxMnTpT6LvG7JlOEQtSQqVOncvLkSe6++27y8vJYs2YN48aNY9++fTzwwAMUFhYyefJkkpKS+M9//sPYsWOJjo4uM+V48eJFRowYwZUrVxg9ejQTJkwgMTGRDRs28P3337Njxw4GDBhgcRxKKcaMGcPBgwcZPHgwY8eOxcXFhdjYWLZu3cp9991HaGio8dgZM2awdu1a2rdvz/Tp0/Hw8OCHH37g/vvv5/Tp07z22msVrrFixQp27tzJ+PHjGT16NK6urhbHZM99Xb16lf79+5OTk8Odd97J1KlTycnJ4fz587z//vu89tpruLnZ9j91b7zxBjt27GDq1Knceeed/Pjjj7z11lscOHCAn3/+GXd39zLHX7hwgYEDB9K9e3dmz55NSkoKHh4eABw8eJAxY8aQk5PDuHHjaN++PTExMaxZs4Zvv/2W/fv3c8sttxjP9dlnnzFr1iz8/Py47777aNSoEdu2bWPUqFEUFhYaz2vNhg0bmD59Okop7rrrLjp27EhiYiIHDx7kk08+4a677mLChAmkp6ezZcsWxo8fT3h4uE3nBnj88cd55513aNmyJffff7/xmnPmzOHYsWP885//rPCZX3/9lVdffZWIiAgefPBBjh07xubNmzl58iSnTp3Cy8vL5usLUWcoIYRDQkNDlaW/QtHR0QpQs2bNKvP68OHDFaCGDBmisrOzja//5z//UYBq1KiRmjx5sioqKjK+98orryhAvf7662XOdeuttypXV1f13XfflXk9KipKNWzYUHXv3t3qfZw4cUIBasKECRXey8/PV1lZWcafP/roIwWoOXPmqMLCQuPrBQUF6q677lKAOnz4sPH15557TgHKx8dHnThxosL5d+7cqQD13HPPOXxfb7/9tgLUW2+9VeH8KSkpVu+/9Dg9PDzU8ePHja/rdDo1ffp0BajXXnvN+Lrh2QJqyZIlFc5XWFiowsLCVMOGDdXRo0fLvLdnzx7l6uqqxo0bZ3wtIyND+fn5KR8fHxUVFVXmPMOGDVOACg0NNTnmnTt3Gl9LSEhQPj4+ysfHp8J1lVLq8uXLxn9fuXKlAtTKlStNfiezZs1SgIqOjja+tnv3bgWozp07q/T0dOPrqampqkOHDgpQP//8s/F1w/MF1FdffVXm/Pfdd58C1Nq1a01eX4i6TqYIhaghL774Ij4+PsafJ02ahLu7O+np6RWyLn/6058A/dSPwbFjx/jll1+YNWsWY8aMKXPuDh06MH/+fGOGwBbe3t4VXvP09MTX19f487vvvouPjw/vvfdemWyOh4cHL774IgBr166tcJ4HHniA7t272zQOR+/L1PgDAgJsuqbBzJkz6dGjh/FnTdN46aWXcHV1NTkF1qxZMxYtWlTh9W3bthETE8NTTz1Fr169yrw3ZMgQxo8fz/bt28nMzARg8+bNZGZmMnfuXDp06GA81t3d3fi92mL16tXk5OSwYMGCCtcFaNWqlc3nMnd+0E9Pls6kNm7c2LjYw9T3NGzYMKZOnVrmtblz5wL67JYQv0cyRShEDSk/LePi4kKTJk3Izc2tsDKsefPmgH46zODAgQMAXL9+3WQN09mzZ43/7Natm9lxdO7cmR49erB27VquXLnChAkTGDFiBOHh4bi43Pz/YLm5uZw8eZIWLVrwyiuvVDhPUVFRmeuW1r9/f7PXL8/e+7rrrrv429/+xiOPPMKOHTsYO3Ysw4cPLzP9ZquhQ4dWeC00NJSQkBB+++23ClN1PXv2NDl1Z7iHqKgok/eQkJCATqfj3Llz9O3b1xg4m7r+oEGDbJ7iNNTpjR492qbj7XXs2DEAYw1daREREQBERkZWeK9Pnz4VXjMEe+np6U4coRC1hwRYQtQQPz+/Cq+5ubmZfR1uBjEAqampAHzzzTd88803Zq+Tk5NjcRxubm789NNPLF26lA0bNrBgwQIAgoODefTRR1m0aBGurq6kpaWhlCI+Pp7nn3/erus1bdrU4hhKs/e+wsLCOHDgAEuXLmX79u2sW7cOgE6dOrFs2TImT55s87XNjbNp06bExMSQlZVFYGCg1eMN97BmzRqL1zPcQ0ZGBgBNmjSpcIyrq2uZa1piOE9VtVbIzMzExcWF4ODgCu81bdoUTdOMWbnSLP03XVJS4vyBClELyBShEHWU4ZfWO++8g1LK7J9Zs2ZZPVdgYCDvvPMO8fHxnD59mnfffZeAgACee+45Xn311TLX69Onj8Xr7dy5s8L57Wm+6sh9devWjfXr15Oamsr+/ftZsmQJCQkJTJ06lX379tl87evXr5t9XdM0GjZsaNN9Ge7hv//9r8V7GD58OIBxui0xMbHCuUpKSkhJSbFp/I0aNQIgPj7epuPt5efnh06nIykpqcJ7iYmJKKVMBlNC1EcSYAlRRxlW0e3fv99p59Q0jc6dO/PII4/www8/ALB161YAGjZsSOfOnTlz5kyVTutU5r7c3d0ZOHAgzz//PG+//TZKKbZt22bz5/fs2VPhtdjYWC5fvkzXrl1tXsln7z307NnT7PX3799PcXGxTecxTMX+73//s3qsYSWnPRkkQ12XqdYQhtfsWZEoxO+ZBFhC1FH9+/dnwIABrF27lv/85z8V3tfpdOzevdvqeWJiYkz2IzJkc0ovoX/88cfJzc1l/vz5JqcCo6OjK93byN77OnLkiMlpKVPjt+azzz7jxIkTxp+VUvz973+npKTErmac48ePp3Xr1rzxxhv8/PPPFd4vKipi7969ZY738/Pj008/5dy5c2WOW7x4sc3XnTVrFr6+vrz++usma6FKZ7YMCwAuX75s1/kBnn/++TLfeUZGhnHa2JaMqRD1gdRgCVGHrV27loiICKZNm8Zbb71F79698fb2Ji4ujv3795OUlER+fr7Fc0RGRnLPPffQv39/unTpQrNmzYiPj2fz5s24uLjwxBNPGI998MEHOXDgAKtXr2bfvn2MGjWKFi1acP36dc6ePcvBgwf58ssvK2xMXJX39fnnn/Phhx8ybNgw2rZti5+fH6dPn2b79u0EBAQYG2raYsyYMQwaNIhp06YRHBzMjh07OHz4MAMHDuSxxx6z+Tyenp6sX7+e22+/neHDhzNy5Ei6d++OpmnExsayZ88eAgMDjQX7/v7+vP3228yePZt+/foxbdo0/P392bZtG97e3sZFDtY0adKEzz77jGnTptG/f3/uvvtuOnbsSHJyMgcPHiQsLIzNmzcD+uJ5b29v3nrrLdLS0ox1VZYCumHDhvHYY4/xzjvv0K1bN/74xz+ilGLDhg1cuXKFxx9/nGHDhtn8PQnxeyYBlhB1WJs2bTh27BhvvPEGmzdvZuXKlbi6utK8eXOGDRvGpEmTrJ6jb9++LFy4kF27dvHNN9+Qnp5Os2bNGDVqFE899RQDBw40HmvYHPiOO+7g448/Ztu2bWRnZ9OkSRPat2/Pa6+9xqhRo6r1vv70pz+Rn5/Pvn37OHToEAUFBbRq1YqHH36Yp556yq69+p588knuvvtu3nrrLS5cuEBAQAD/93//xwsvvGDz9KBBv379OH78OCtWrGD79u3s27cPT09PWrZsyYQJE4ytNwxmzZqFv78/y5cvZ/Xq1fj7+3P33Xfz6quvmmy5YM7EiRM5ePAgL7/8Mrt372br1q0EBQURHh7O/PnzjccFBASwfv16li5dyscff0xeXh5gOcACePvtt+nVqxf/+te/+OijjwDo2rUry5YtsyuYFeL3TvYiFELUe0uXLuX5559n586dJlsQCCGEvaQGSwghhBDCySTAEkIIIYRwMgmwhBBCCCGcTGqwhBBCCCGcTDJYQgghhBBOJgGWEEIIIYSTVWuAlZuby9GjR8nNza3OywohhBBCVKtqDbDOnj1Lnz59jN2L6xvDTvfi902ec/0gz7l+kOdcfzj7WcsUYTWyZ1NVUXfJc64f5DnXD/Kc6w9nP2sJsIQQQgghnEwCLCGEEEIIJ5MASwghhBDCydxqegAGOp2Oa9eukZSUhE6nq+nhVImMjAz8/f1rehiiitXV5+zq6kpYWFidHLsQQtQ2tSLAunTpEnPnzmX37t01PRQh6jUvLy+2bNnC6NGja3ooQghRp9V4gFVQUEDPnj0JCgriyy+/pG3btri51fiwhKh3CgsLWbZsGePHjychIUEyWUIIUQk1HslERUWRnZ3Nt99+y5AhQ2p6OELUa0uWLOHbb78lJiaGnj171vRwhBCizqrxAKu4uBiABg0alHk9LluRnO/86wV5QWtfzfknFuJ3wMPDA5DeP0IIUVk1HmCZEpet6LSumLwq+N94b1c4O8XNapD13XffsXjxYgoLC2nQoAEffvih8f/RjxgxgtjYWOMUyqxZs3jiiScAePDBB/nll18IDg5m06ZN+Pv7o5Tijjvu4N1336Vt27bOvykLTp06xbhx44iJieHq1atMnTqVPXv2WPzMW2+9xbRp02jWrBkAH3zwAVlZWTz11FPVMWSH3H///ezfv5+OHTsyZ84cdu7cyZtvvlnhuNLfx+9FeHg4e/bsoWHDhjU9FCGEcFh6gULTwN/j95EEqZUBVnI+5JXAxDCNYC/nnTcpHzbF6DNjrX3NH5eWlsaMGTP4+eef6dq1K3v27GHGjBmcOnXKeMybb77JhAkTynzu1KlTnD9/npMnT7Js2TI+//xzHn30Uf79738TERHhlOCquLjY4Rq1Fi1aWA2uQB9gjRgxwhhgPfTQQw5dr7pcv36dr776iszMTFxdXQG4++67a3hUtqvMMwWIjIx04miEEKJmbItTlCiY2R40re4HWbW6D1awFzRvoDntj63B2sWLFwkMDKRr164ADB06lLi4OI4ePWrxc+7u7hQUFKDT6cjJycHDw4Nr166xdu1annzySYufHTFiBI899hj9+vWjXbt2LFiwAKWU8b3HH3+cQYMGGVd3ff/99wwZMoQ+ffrQv39/du7caTzX0qVLad++PX369OGrr74yvh4TE0OjRo2MP+/fv58hQ4bQs2dPevTowZYtW1i2bJkx0xUeHk5kZCRLly7lL3/5CwAdOnTg8OHDxnOsWrWKiRMnApCQkMCUKVPo378/3bt3Z/HixWbv95VXXqF79+707NmTgQMHGjcAX7FiBV27dqV79+7MmDHDuDfU0qVLmTp1KnfddRddunRh5MiRpKamkp6eTkREBPn5+fTp04d//OMfrFq1qkzwa+77sPQ97tq1i27duvHnP/+Znj170rVr1zL3/c0339CvXz969uxJeHg4Bw8eBODXX39l5MiR9O3bl169evH111+bvP/Zs2czd+5chg0bRrdu3ax+9sMPP6RDhw707t2bF154ocz/+GiaRnp6OgCHDx/m1ltvpUePHvTv3599+/aVefbPPfccffr0oV27dmzfvt3s8xFCiOp0IkVxIVMRnaXYn6hqejjOoarRkSNHFKCOHDli+bUkneKjQvXAz0XqucPFTvvzwM9Fio8K1ZEkncVxpqenq4CAALVv3z6llFJbtmxRgNqwYYNSSqnhw4erjh07qm7duqkpU6aoixcvGj+7aNEi1bNnTzVp0iSVk5Ojpk6dqg4fPmz1uxk+fLgaOXKkKiwsVDk5OapPnz5qzZo1xvfGjBmjCgsLlVJKXbx4UQ0cOFBlZGQopZQ6f/68atasmcrPz1fbtm1TXbp0URkZGUqn06kZM2ao0NBQpZRS0dHRyt/fXymlVEpKimrSpIn6+eeflVJKlZSUqJSUFKWUUqGhoerYsWPGsT333HPq//7v/5RSSr344ovqkUceMb43bNgwtXXrVqWUUqNHj1a7du1SSilVVFSkxowZo9atW1fhXletWqX69eun0tPTlVJKpaamquLiYrV9+3bVqVMnlZaWppRSav78+eqhhx4yjiE0NFQlJycrpZSaOnWqeumllyrcl1JKrVy5Uo0fP14ppSx+H5a+x507dypXV1d14MABpZRS//rXv9To0aOVUkpFRUWp4OBgdebMGaWUUoWFhSo9PV2lpaWp8PBwdfXqVaWUUklJSSokJERduXKlwncwa9Ys1aNHD5WZmamUUhY/e/LkSdWsWTN17do1pZRSS5YsUaX/6gIqLS1NFRQUqJCQEPXdd98ppZTas2ePatq0qcrKylLR0dEKUOvXr1dKKfXtt9+qDh06VBiXqb+PwjTD3xfx+ybPuerlFOnUK5E3f1e/cKRYJeRY/j1dFZz9rGvlFGFN8/f3Z/369fztb38jOzubQYMG0aVLF+M0zueff05ISAhKKd577z3GjRvH6dOnAVi+fDnLly8HYMuWLYSEhBAWFsacOXPIzMxkypQpTJ061eR1Z86cibu7O+7u7tx77738+OOPTJ8+HYB7770Xd3d3QF8fduHCBYYNG2b8rIuLC3FxcezYsYMpU6bg5+cH6GvC9u7dW+FahnqloUOHGj8fEBBg9buZOXMmvXr14vXXXyc+Pp5z585x++23k5OTw44dO7h+/brx2OzsbKKioiqcY9u2bTz00EPGGrbGjRsD8OOPPzJ16lRjlu3hhx9m8uTJxs+NHTuWwMBAAAYNGsTJkyetjtfS92HpewRo164dAwYMMF7vtddeA+CHH35g7NixdOrUCdBnLv39/dm+fTuXLl3i9ttvLzOGqKgoWrZsWWFskydPNtZN/fLLL2Y/e+rUKcaOHWucsp0/fz7Lli2rcL6oqChcXFwYM2YMAEOGDKFp06ZERkbSqlUrvLy8uOeee4z3c/HiRavfnxBCVLXvLytyi2/+XKxgY4yO+Z1ccHOpu1OFEmCZERERQUREBKDv1dWsWTO6dOkCQEhICKCfmnn00Uf561//SkpKivGXP0BmZiavvfYa33//PS+//DLDhw/n3nvvpWfPntx99914e3tbHUPpaSBf35tFY0opbrvtNr788ku7zuEMrVq1om/fvmzZsoXffvuNe++9Fzc3N/Lz9Us+Dxw4gJeXcwrnyo+99HldXV2NK1AdPael7zE+Pt7u6yml6Nq1K7/88otNYyn/TM19tnTtX/l7sKb0sZ6ensafXV1dZaWgEKLGXcxUHE+tOCV4PQ9+uqoY3aruBli1ugarJl27ds347y+88AIjR46kXbt2FBcXl8nSbNiwgaZNm5YJrgCeeeYZlixZQoMGDcjJyUHTNDRNo6ioiMLCQpPX/OKLLygqKiIvL48vv/ySUaNGmTxuzJgx/Pjjj5w4ccL42qFDhwAYNWoUX3/9NVlZWSil+Oijj0ye49Zbb+X8+fPGonedTkdqaioAfn5+xtonU+bMmcOnn37KZ599xty5cwF9sBAREcE//vEP43FXr17lypUrFT5/991388EHHxivkZ6eTklJCaNGjWLdunVkZmYC+rqjynYUt/R9WPoeLRkzZgzff/89Z8+eBaCoqIiMjAxuvfVWoqOj+fHHH43HRkZGmn3epVn6bEREBN9//z2JiYkAfPLJJybP0bFjR3Q6HT/88AOgz4olJCQQHh5u9fpCCFHdinWKbXHmt8bbf12RUVh367FqdQYrKR/AeV9ukh19tZYsWcKePXsoLi5m0KBBxl9qBQUF3HnnnRQUFODi4kJQUBBbt24t89l9+/aRl5fHbbfdBsAjjzzCn/70J1555RXuu+8+sx2yO3fuzODBg0lNTWX8+PFMmzbN5HHt2rXjyy+/5MEHHyQ3N5fCwkJ69erFl19+yR133MGhQ4fo3bs3fn5+FaacDBo3bsymTZtYsGABWVlZuLi48MILL3DXXXfx+OOPM3/+fBo0aMCqVasqfHb8+PE8/PDDtG/fns6dOxtfX7NmDU8++STdunVD0zR8fHz48MMPadWqVZnP33fffVy9epVbb70VNzc3fHx8+PHHH7n99ts5deoUgwYNwsXFhR49evD++++bfkA2svR9WPoeLWnXrh0rV67k3nvvpaioCFdXVz744AP69+/PN998w1//+lcWLFhAUVERrVu3ZvPmzVbH2bhxY7OfNSwYGDx4MA0bNmTs2LEm/xvy8PBg48aNPP744yxYsAAvLy/Wr1+Pr68vycnJ9n95QghRhS5mQlqB+fcVEJsFPQLNH1ObaUqpagsPjx49Sp8+fThy5Ai9e/c2+1pt6INV3UaMGMFf/vKXCq0fhADIysoy1mv985//5LvvvuPbb791+nVM/X0UpqWmptpUtyjqNnnOVWdrrI6jyZZDkL5BGuNCq2eyzdnPulZmsFr7apyd4iad3IW44ZlnnmHfvn0UFRXRokULPvzww5oekhBCOEwpxbkM6/mduByZInS61r6axWagvze7du2q6SGIWuy9996r6SEIIYTTXM2F7KKbP5/5aRN7/72c5NgogkI7MmTeYjqPnEhSHuQXK7zc6l5iRIrchRBCCFGtSmevzvy0iXVPT6ZtQTwL+4fStiCedU9P5sxPm1DA5ZyaG2dlSIAlhBBCiGoVlX4zwNr77+UMDQ1izfjePNA7jDXjezOkdRD7PtH3lIzLrpvThBJgCSGEEKLaZBQqEvJu/pwcG8WI1oHGPn2aphERGkhSjL5RdVx2TYyy8iTAEkIIIUS1KV/cHhTakV1xKcb9d5VS7IxNITisIwBXcxUlurqXxaq1Re62KCkpYc+ePVy7do3mzZszdOhQXF1da3pYQgghhDDjXEbZonafgKbsuZbM9M1HiQgNZGdsCnvjkpmy4gMAinT6oviQOrbwrc5msDZu3Ei7W9oQERHB9OnTiYiIoN0tbdi4caPTrlFYWMjChQtp164dnTt3pnv37qxevRrQr/qrjR2yd+3axXfffWf8+erVq8b9BoUQQoiaVFii2LZ5Y5mi9m4e+ir2MyW+vHIolmjvlkxZsZ7OETf7Ql6ug+0a6mQGa+PGjUyaNIk/tAnmn1MG0DHQl6iUbN45HM2kSZNYv369cVPbypg9ezYFBQUcP34cHx8fYmJiuP322ykuLqZt27ZOuBPTiouLjRtL22vXrl2kp6czduxYAFq0aGHcDkcIIYSoSZey4OePbxa1a5rG/F6hTN98lGhvfx7dcgHQZ7g+nt7H2LYh8/8Wc+ufJ9Xw6O1T5zJYJSUlLHjiL/yhTTCfjAund/NG+Hi40bt5Iz4ZF84f2gTz1yefqPRGtufPn2fz5s189NFH+Pj4ABAWFsbrr7/O888/D+gDoZkzZ9KtWzf69OlDZGSk8bODBw+mZ8+exm1OQL9n3TPPPEP//v0JDw9nypQppKWlAfpgbu7cuQwbNoxu3brx4osv8uijjxrHk52dTUBAAElJSZw8eZIhQ4bQu3dvunTpwvLl+pUWkZGRfPDBB6xZs4bw8HCWLVtGTEwMjRo1Mp7n+++/p3fv3vTo0YPhw4dz+vRpQB+YdevWjT//+c/07NmTrl27cvjwYQCSkpIYPXo03bt3p0ePHsyZM6dS360QQoj6KSpdWS1qN9W24c1HJjt1hqo61LkAa8+ePcTEXeaxvm1w0co2HnPRNB7t24bo2LhKZ22OHTtG+/btK2ziPGjQIC5fvkxSUhK//fYbs2bN4tSpUyxcuJBp06ahlOLdd99l3LhxHD9+nJMnT/Lkk08CsGLFCnx8fDh06BCRkZFlgi+AI0eO8M0333D27FlmzpzJunXrKCjQb9T09ddfExERQXBwMGFhYezYsYOjR49y5MgRNmzYwIEDBwgPD+ehhx5ixowZREZGsmTJkjJjT0xMZPr06axevZoTJ07wwAMPMGnSJGNh4dmzZ5k1axbHjx/nscceY9GiRYB+E+o2bdpw8uRJTpw4weuvv16p71YIIUT9o5TifKayWtRurm3DsmUv1OTw7VbnAqxr164B0DHQdLWb4XXDcVUpLCyMP/zhDwBMmTKFhIQELl++zLBhw/j4449ZtGgR//vf/4wZpM2bN/PFF18QHh5OeHg4a9euJTo62ni+yZMnG/ebCwkJoVevXsaNpFetWmXMHOXl5TFv3jy6d+/OwIEDiY2NNWbPLDl48CDdu3ene/fuAMyYMYOrV68SHx8P6DcxHjBg6xViAAAgAElEQVRgAKAPJC9evAjAwIED+fbbb1mwYAFbtmwxZvSEEEIIW13O0XdvHzJvMXti9UXtHx2NYfrmo+yNS2bwvGcB820boqLO1uTw7VbnAqzmzZsDEJViujGG4XXDcY7q1asX58+fJyUlpczr+/fvJyQkhODg4Aqf0TQNTdP44x//yL59++jYsaMxmwX6KP2dd94hMjKSyMhITp8+zfbt242f9/UtGzTOnTuXlStXcunSJS5cuGCsq/r73/9OUFAQx44d4/jx44wYMYL8/Mpv3Ojl5WX8d1dXV4qLiwF9sBUZGcmAAQPYuHEj/fr1q/QUrBBCiPrleIo+Y9V55ESmvPo10d4tTRa1m8twtWrbqcbG7og6F2ANHTqUsNYhvHM4Gp0qu6pApxTvHo6mTWjrSq+ca9++PXfddRcPPPAAubm5AMTExLBgwQKeffZZ4887d+4EYP369TRt2pRWrVpx/vx5mjZtysyZM3n11Vc5cOAAABMmTODNN980ni83N5fffvvN7BgmTJjAr7/+yssvv8y9995rLHxPS0ujVatWuLm5ERUVxQ8//GD8jJ+fHxkZGSbPN3DgQE6ePMmpU6cA+Oqrr2jZsiUtW7a0+F1ER0fj6+vLlClTeOeddzh37hzZ2XW085sQQohqV6xT/JZ283d255ETmbfmCH/bl828NUfKrBg0l+EaNn+xqVPXWnVuFaGrqyuvv/kWkyZN4v5tkTzat41xFeG7h6PZEZ3E+vXrndIP67PPPmPx4sV0794dDw8PXF1deeqpp5g7dy67du2ia9eurFq1iscffxwPDw/Wrl2LpmmsX7+eL774Ag8PD3Q6HR98oO/lsXDhQgoKChgwYIAx9blw4UK6du1q8vqenp5MmTKF999/nzNnzhhfX7x4Mffddx+rV6+mbdu2jBw50vjexIkT+fzzzwkPD+eee+5h5syZxveCg4NZs2YNM2fOpLi4mMaNG/P1118bx2LOrl27eOONN4xZrRUrVuDv7+/w9yqEEKJ+iUqHfBsnPgwZrn2fLOfQoSiCwzoyZcUHhAyZQE6Rwse9bmz8rCmlqq25xNGjR+nTpw9Hjhyhd+/eZl+zxcaNG1nwxF+IibtsfK1NaGtee+NNp7RoEKI+cvTvY32UmppKQEBATQ9DVDF5zs6x9oKOqIzKhxv3tXehrZ/zAqyUfEWgl/58zn7WdS6DZXDPPfcwfvx46eQuhBBC1GI5RYoLmc7J5STmKacGWHHZEOhl/ThH1NkAC/TThSNGjKjpYQghhBDCjFNpihInzZUl5lk/xlbpBYqMQgVUzZRjnStyF0IIIUTdcTzF+jG2Sqr8gnmjuCpeq1XjAZZhZZxhZZ0QouYUFhYCyFS7EMIpkvIUV3OdV+qdlK9wVul4XHbVlqDX+BRhx44d8fX15b777uOll16ibdu2Du/DJ4RwXGFhIcuWLcPLy4uwsLCaHo4Qwg5FOkVWIQR41a4VdsdTnRvEFJRARiE08qz8ueJyFF0aVd33ZTKSWblyJXPnzmXTpk1MmDCBxMREZs6cycWLF/H09OT9999n2LBhABbfs4WnpyfHjx9n7ty5TJ8+3Tl3JYRwiJeXF1u2bJE2HELUMTuvKmKyYF4nKmwjV1OUUpx0coAF+jqsygZYecWKpDygkdVDHVYhwIqJieHjjz9m4MCBxteeeeYZBg4cyHfffcevv/7KxIkTiY6Oxt3d3eJ7trrlllv46aefSEhIIDExEZ1O55y7q2UyMjLkF1c9UFefs6urK2FhYXVy7ELUZ1dzFAcSFToFu69BRIvaEWBdytJnm5wtMV/RoZKF6Zezoap7VJUJsHQ6HfPmzeOdd95hwYIFxtfXrVvHhQsXAOjXrx8tWrRg9+7djBo1yuJ79nBxcaFFixa0aNGisvdUa0k/lfpBnrMQorrolGJrrA7djWhhT4Kig7+ipU/NB1kHEqsmhHHGSsK4nKpvAVomwHrjjTcYPHgwffr0Mb6WkpJCUVERzZo1M74WFhZGXFycxfcsyc7OJjMz0/izp6cnnp5OmFAVQggh6pF9CYqd32xi77+XkxwbRVBoR+IeXsy/Hv8j7i41F2Sl5CsuOKGxqClJ+ZU/b1WvIIRSAdapU6fYsGEDP//8c5VfdPjw4WV+fvrpp1m4cGGVX7empaWl1fQQRDWQ51w/yHOuH2rzc04p0Hhr5XY2LrqXoaFB3N8/lF1x8XzyxGT8tc9YdN+dNTa2H665kptXNauRY/MhJaUQR0vNinVwIdmDEgWZWSWkpur38HH0WZubsTAGWHv27CEmJob27dsDkJCQwAMPPMDzzz+Pm5sbCQkJxkxVTEwMrVu3JjAw0Ox7luzevZvw8HDjz/UpgyVTR/WDPOf6QZ5z/VBbn/OWKB0HVq9gaGgQa8b3RtM05vcKZfrmo/znX6/z9Px7adqg+rNYecWKS7E6vL2r8CI+Pg6vmIzNUnh46Wu9/RpqBATc7FjlzGdtPOvDDz/MtWvXiImJISYmhoEDB/LRRx/x8MMPM3nyZOOGxb/++ivx8fHGLJSl98zx9fXFz8/P+Ke+BFdCCCGEMyTmKWKzFcmxUYxoHYh2I52jaRoRoYEkxUQRnVVtWw2XcTRZUVhurdqZnzbx8fQ+vDzYl4+n9+HMT5sqdY3K1GFVdf8rA5sajb7yyiv88ssvtG/fntmzZ/PFF18YVwlaek8IIYQQznfuRn1TUGhHdsWlGJtvKqXYGZtCcFjHaqkzKk+nFIeSygYwZ37axLqnJ9O2IJ6F/UNpWxDPuqcnVyrISqxEHVZcjsMftYvZjp67du0y/nvTpk353//+Z/I4S+8JIYQQwvmi0vX/HDJvMeuensz0zUeJCA1kZ2wKe+OSmbLig2pZKVfembSKrRn2/nu5yWnMfZ8sp/PIiQ5dx9EMllKKy7UpgyWEEEKI2iGnSHHlRvDUeeREprz6NdHeLXnlUCzR3i2ZsmI9nSMmkF0EaQXVG2TtN9GawdI0pqMS8xy7r8Q8yC9x+LJ2kT1phBBCiDrkfGbZJpmdR040mwmKy4bG1VTmHJ9zM/ArTT+NGc/8XqFomlZmGtNRqQVQolO42tmKojqzepLBEkIIIeqQc+m2BwnVVdANcMrMtjhD5i1mT2wy0zcf5aOjMUzffJS9cckMnvdsmeN0yvaNnEsUpBTYP8bqrEuTAEsIIYSoI0p0iot2rA68XI0ZmygzjUUtTWMa6JTi8/OKLy+aDrJMrUJ0pA6rOgNOmSIUQggh6oiYbCiwo4YoKU/fl8rbrWr7YSXmKVItZJQsTWMC7E3Q3xvAsRToHXTzPcMqxNLNVNc9PZkujb6m2/xJNo8xo1BVyd6I5kgGSwghhKgjouyYHgR9rdblapgWO2vnuEq7mqPYfU0xtBn0DID/xSuyim6er/QqxAd6h7FmfG+GtA7is7detOs61d22QgIsIYQQoo4458D+ftUxTXg23bHPFekUG2MUTRvA8OYao1tpuGnw7eWbYza3CjH+4lm7rhWbVbE/14Oje+PTwJs+vcL573//69hNmCEBlhBCCFEHXM9VpDswxVXVmZvMQsW1XMeCuB/i9dN294RpuGoaDdw0xrbSOJMOZ9ItN1MNCutIkc6Ogv9SgaZh2rFZdjwL+obim3aV2bNns3HjRofuwxQJsIQQQog6wNbsVfmC8B/+u5ESOwIRe0VlKBw5+4UMxa9JcFsrjaBS+wp2bQwd/GD7ZUV+sbK4CjHJxkL3/GJV5ljDtOMXN6Ydvxjfi6GhQby0/AUH7sQ0CbCEEEKIOiAqw/oxpralWfvUZD5d67zMTHmOTg/+L15xS0PoF1T2dU3TuLO1RmEJ/HRNWVyFGGPjqsC47LK9w0xNO45oHciZs/ZNO1oiqwiFEEKIWi6rUBFvQy2VuW1p3vjHcubP+KPTx1VQoohxYFPplHxFUj6MbKEZg5zS/Dw0+gYpTqTC7a2U2VWIlzLh1qbWr1e+waip5qe74lLo3KmT3fdijmSwhBBCiFruROrNaThTPaEMzBWER19wXmamtPMZ+qaf9jqbAW4atPUzf0xrX43s4op7G5YWl61smv4sX4dmmHacsUU/7Thjy1H2xCaz6NklNt6BdRJgCSGEELXciRtd0k1NAa57erIxyDJXEN6kEtvSWOJoe4az6Yp2fuBuYaubEF/9Py/nmD9PoQ6uWHgfoFinuFoug2WYdkxs2IrXD8eSG9CK1atXM3GiY5tPmyIBlhBCCFGLJeQqrt8o0DbXE2rfJ8sB89vSDLr/WVLynVvoXqJTnM+0/5xZRYorOdCxkeXmpw3cNAI9rbeZiLYyRXk1F4pNHNJ55EQ++P4IObl5HD56jHHjxlkduz0kwBJCCCFqseMp1ntCJcVEAZa3pbmW69xx2dtV3uBcBmhAB3/rx4b4WG+UeinL8vvVuT1OaVLkLoQQQtRSOqU4mXYzQDBVnL0zNoXgUlOA5grCE/MU+tDGOS45kL0C/fRgqK8+Q2VNK1+N46mKwhKFh6vp4+NzLL9f3R3cDSSDJYQQQtRSFzMhu+jmz5Z6QlmTmO/csVkKXMwV4heUKKKzoJOV6UGDEB99e4V4C9m3EgWxZsailOJyDWWwJMASQgghaqnS04NgeQrQmqQ85wUaxTrF1dyyndENAdXbd7U1W4hvWHXY0YbpQYBgL/B0tT5NaK4OKykf8hyYxnQGmSIUQghR7+UWK5umrKpTQYkiykT3dnNTgNakFugDIzcLK/dsFZ9zsz2DYWXj0NAg7u8fyj8PXWRISGCFXlz7PllO4+XjaeYNjTxtG4OmaYT4KK7kWJ7eNDddWVP1VyAZLCGEEPXc2XTFmgsKnaq5X8am/JYGRTrnnU+hz+g4Q2ypwKX8ysaiEh0RYUEmC/HPZdg+PWjQykfjcg7G1hOmXM+DnKKK79vS/b6qSIAlhBCi3sovVnwTpyM+R7E3oXYFWOWnB50h0ca9+6wpXX9VfmVjuwBfdsUmV+jF5RfSkUIddGpk37VCfCC/BJItBIcKiCm3mvCneB3nbdy/sSpIgCWEEKLe+iFekXWjiHz3NUVCbu0IstIKVJVMbyU5oReWUqpMb6ryzU0f7duGPXEp/GnTkTKF+H73LKKxBzTxsu96LX30k4OWGo4CXCpVh3U4ScfPNRwwS4AlhBCiXorOVBxNvvlLuETB5hidTVuvVLWDiYqqGIUzMljX88r2vyq/svHzU/EAnNX5GQvxJ736Ndc7TaBTI0zuPWiJp6tGU2/bG46eTVd8E1fzz1CK3IUQQtQ7RTrFf+N0FYKYhDx9Jmtky5oreC8oURwrNT145qdN7P33cpJjowgK7ciQeYsdKnIHQy+syoktl1kzrGzc98lyDh2KIjisI1NWfFBmZeOZdEXOJUXPQMe+11Y+FacAy0stgJOpiq2xFZ9rTZAASwghRL2z66oitcD0e3uvK7o0VjRrUDNB1tFkZcwQlV+htytO3/JgyqtfOxRkZRRisSmnLUz1v7K2sjEyWdGiATT1duy6Ib4ah5OV1dWeG6NrR3AFMkUohBCininSKQ4lmf81rFOOb2JcWUqVHZu1vQdtUbpH1UfT+7D6PxsrNUZ7a8OyChXnM6GXg9kr0Be6g/WNnc2NzFzj06okAZYQQoh65VKm9fYH1gqqq8rZdEgrlVmztvegwY54HWsu6Cq0mjBkwEo3/XxgxiQ2bnQsyEoruLkowFaRqeCqQbcAhy4JQCMP8HHDoa7spr4DQ+PTqiQBlhBCiHrlnA1L96/k1ExfrAOJZa9ZfoWeqb0HL2Qo9l6HC5lwKKns+UxlwIaGBvHS8hccGp+5LWnMUUoRmaLo0hi8KjEtqW84ejODZU9GyhlZQEdIgCWEEKLeUErZFGAVlDivZ5StruaoCgXk1vYezC1WbIlTtPWDvkHw01VFesHNc5jKgI1oHciZs2cdGqO904Ox2fri88pMDxo0b6CRkAend2y0KyNlaxbQ2STAEkIIUW9cy6XMFJelTEh1b7NiyF6VHtPefy/n1vv+anLvQaX07QiKdTC+tcaolhrervDNZWXMeJnKgO2KS6Fzp04OjdHwndiaQYpMUTT2hFBfhy5XRhNvfcPRPXZmpCxlAS11h68sCbCEEELUG6WzV9Zqc6qzDiu7SPFbmjI5pl8+f43B9y/mb/uymbfmiLH9wclUOJ0O41prNPTQ8HTVuCNE40ImnErTn9dUBmxPbDJP/f1Zu8eYU6RIzre9pim/RPFbmj57ZW/vK1OCbzQoTYk7Z1dGylwWsOesxUzaoXPqJtilSYAlhBCi3ii9ebK12pzqzGCdTlOUKNvrhTIKFdsvK7o3hq6NbwYvHRtpdGkE313RtzQw9KgqnwEbNNb+Fg+G9gy2jvFUqr55a89KFLeX1tgT3DTwadnBal1aaea+A9VrApmF4OPunPGVJ32whBBC1AuZhYpruTd/To6N4v7+oRUyIYcO6TMhGYWQWQROig8sOptu25gMdsQrPF3hjpCKmaHbQzTeO6344YpifJhmskdVUh60tnPazhBw2jJGpfTNUtv5gZ+Hc/qJuWgaQV4K7z8uYs+KqUzffJSI0EB2xqawNy6ZKSs+MPtZU9/B+ks6OjXCYl+tSo23Ss4qhBBC1DLli9ttWaEXn1v1vybzi28Wt9sypmKdIioD+gRpeJkIDnzdNQY31TiVhtltfxId2JMw7saUqS1jjMuBq7nQN9i5wUuwFxSHm85Ile4cb41Siths6BFQdc1kJYMlhBCiXjiXUfbnIfMWs+7pyRYzIfF5Vd/N/XymfirN1jHFZEGhDjo1Mn/ONg1hx1X91j8tfSq+b+8KycISxbUbG2HbMsZ9CYomXtDez77rWBPsrXEuU9EpYoLD2wWBfmVjdjH0rMIASzJYQgghfveKdIpLmab30LOUCblSDRms0l3jbRnT2Qz9yjxD0bcpzbz19UrmCvWT7MxgXcnRd7i3ZYzXc/Wd2wc3c05xe2nBXvoWGtl2NjstLzYbNKBrY6cMyyTJYAkhhPjdu5QJxcr0xsnz1hwx+7mkAq3Se/dZUqJTXDAR+JnLziiliEqH7gFYDF5cXTRa+CguZysGNql4XHYRVvf1K618wb+lMe69rmjkAd2qIHhpciOoTMyHhh6OnycuW9GsAfi4SwZLCCGEcFhUhukWCNa2TNEp6/vfVcalLIwbO9siPlc/tdWxkfXAoJWPPoNlrteTPdOEpjZ4NiWtQN+aYVBTDRcnZ68AGt1YSZiUX7nzxGQ7pzeXJRJgCSGE+F0zdG93dMuUqmzXYO+m0mfTFQ3cbm5+bEmIj0ZWkX4lpClXcmy7tk4pm4/95brC2w16Bdp0uN30KwkhsRK9q9ILFBmFEOpbtfV1EmAJIYT4XbuUpZ8Sc3TLlKpqOGrrtj2lRaVDB39syg4ZgrDLZrJP0Vm2XTMhV19Ub012keJYCgxsouHuUnXBS7CX9QyW5Q79+n/a26bCXhJgCSGE+F07nmJ7CwRTqmrj5/icstv2WJOcr0gugE7+tgUvPu4aAZ5w2Uz2KS5bmW3jUFr5/RHNOZiocNWgX5BNhzss2FsjKc/81Ke1qeDYbEWwV9X1vzKQAEsIIcTvVmGJMk7DWds4GSpmPs7t3lJlGz+ftTN7dTYd3F3gFjtaH4T4mK8hK9LZVl9mS/1VsU7xazL0CcJkby5nauIFBTrzwam1qeDYaqi/AgmwhBBC/I6dTr85vWWtvYCpzMfmxfdy5qdNNtcg2cPe+quoDEXbhtg1/dbKR+Narj7QNCU6y/oYbKlBu5ipL9bvHVT1fcOCvfX/NDdNaGkqOLtIkVJQ9fVXIG0ahBBC/I4ZpgcNLLYXKJX50DSN+b1Cmb75KPs+Wc6sqfc4dVwp+fqNk22VVaS4kgMTQu0LDEJ8QaHvqh7WsOL70Vkwwso4c4qtX+d0un7aLcir6gOXRh76lYSJedDWRDZPPxUcz/xe+u18Sk8FG7JxksESQgghHJRZqIixIUNjYCnzkVjJtgDlXSo1LksF2QbnMvSNMdv723edYC/wdDFfqH8lR1FkoQ7LlunBkhtb93S20FnemQwrCc01S7U0FRybrQjwhIZO2h/R4jir/ApCCCFEDTiRqrBnEs5SEXxSJdoCmGIIXGztzXUmTRHqa39htoum0dIHrpiZ5itREGthNaEt04PRN3p5dbGhN5clLpo+iLRFE2/zU4SWpoJjs6t+9aCBTBEKIYRwWEGJYudVxdiQ2vf/18tPD1pjaY+9nGLIKVJO6/wdeyODZWla0jCVmVGouJgFd7V27NohPnAoSR8wmur+Hp2laGdmZaItKwhPp+uzQk28HRqe0R9aaKQWwJFk69cM9tKISldm78nUVHBeseJ6HiY721eF2vc3QgghRJ3xwxXFgUTFqdSqa8bpiKs5yu5u36YyHxNeXGMsgq9s93CD9AJlbP5pS2+u4yn61YOO7psX4quRVwIpBabfN9cPK7tIkWrmMwY6dXN6sDL7DnZqpDG4mQt/aKnh7Wr9+OByKwltmWatzvorkABLCCGEg2KylDHbsP2yjpyi2hNkHXcw4Os8ciLz1hzhb/uymbfmCB2G3W18rzLdw0srXddkrTeXUorIFEXXRuDp4H6ILQ0NR83UYV3LVeQXV7w3W+qvYrMhtxg6V2J6MMDzZvF+AzeNiBbWz2VYSZiYb/s0a2y2ws9dXyRfHSoEWKNHj6ZHjx6Eh4czdOhQjh07BsB3331H37596dGjBwMHDuT48ePGz4wYMYI2bdoQHh5OeHg4b775ZvWMXgghRI0o0im2xuqMNU65xbAtrnYEWCU6+zJqBxIVH57RUWCmlYGBs3phlZ52s9abKzYb0gohvBLtD7xcNZp43azDKp/tOf3TJmLKBVM6pTiWYr0Q/3Sawt8DWjRwbGxuGky5xaVM76y+wRrNrEw3Nr6xkjApz3rfK4O4G/2vKpNps0eFAGvdunWcOHGCyMhInnzySWbPnk1aWhozZsxg9erVnDhxghUrVjBjxowyn3vzzTeJjIwkMjKSJ554oloGL4QQombsvFpx+uhMeu2YKjyVhk2tBUCfvfnhiiIhD366annszpoiLF04bq0317EUfX1Taxv2HrQkxFefwTKX7fli3YYyx2+LU5zPuBmQmfrM6R0bOVvJ6cE7WrvQrEHZz7poGne0tjzBpmkawV6QmK9smmYtLFFczYXW1dD/yqDCHTRqdHOdZUZGBpqmcfHiRQIDA+natSsAQ4cOJS4ujqNHj1bbQIUQQtQOV3P0dVembL+sI7uGpwoPJNqwcR767uObYhTB3voC60NJljdAdsYUYW5xxf5X5aclDcFVfonidBqEB2qVzrq0aKCRlA97/v2CyWzPF2/dzPb8cEXH0VKF5uYyRLv+/SLZRY5PD7ZooJltTNraV6NHgOXzBnvrM1i2bIF0OUffD6y66q/ATA3WzJkzCQkJ4dlnn+Xzzz+nffv2pKSk8MsvvwCwdetWsrKyiImJMX7mmWeeoXv37kydOpVLly5ZvGh2djaZmZnGPwUFVqrohBBC1BrfX1GYa52UWwyHzARf1SE2S3Et17Zjf7qRhbsnTOPWptC8Afw3VlFiZo+7vBIqHTzGZWNz64hTqfo2CuEBlbokcHOFX3LsOZPZnmuX9F3O9ybo2He97AjNZYjSYqPwdb+5qbS9ulgp2r+tpYalrhRNvTUS8+HW+xdZ3QIpNlvRwA2CvBwbqyNMtmn47LPPAFi9ejULFy5k+/btrF+/nr/97W9kZ2czaNAgunTpgpub/uOff/45ISEhKKV47733GDduHKdPnzZ70eHDh5f5+emnn2bhwoXOuqdaKy0traaHIKqBPOf6ob4+55xiiErysBgknL6uI9zbxjk6J/vhsht5edbXb8XlaOxPdGdEkxIaUkJBPoxuqvFZtDu7rxQzKKgEgLy8summc9eKCPN1PMj67boreXk2LJMDjia508ZX4VZSTF4l67/0Q/bEt2V7dsVdq9DlPDC0A1/+lsXF7IrfXWBr059xad6R9r4l5Oc79qybU0hqquVj/HHjipnn2cxdo0jngUf4HUxY/gUHP3uVQ4fOERjagQkv/pOwgWPIu/HFRWe608pbVRhrZlYJqan6Z+3o3+mAANMRsMU+WLNmzeKhhx4iJSWFiIgIIiIiACgoKKBZs2Z06dIFgJCQEEAf1T766KP89a9/JSUlhcDAQJPn3b17N+Hh4cafPT098fT0tP+u6iBzD0L8vshzrh/q43OOS1F4eVuegsvQoFFjF1yqqZjYIK1AEV+sw9tKgXRBieLbi/rGnUNbuuOi6ZeVhXnDoFwdvyS60TPYncAb2754lzphoWcDAgIcX4CflqjD29t6gJaYp7iWr5jSRsPb293h6xl4A37uOppOWcyef0w12evraomPye9u6PwlJvuD8fD7dA10c2h8zbyhbTPr83WdcnWkXDf9fYV5KTzjFNcKPRk2dho9x04zeVyxTnEtTzGqZcXv0q+hVuZ5OvPvdJn/StLT07l69arx582bNxMYGEhAQADXrl0zvv7CCy8wcuRI2rVrR3FxMdevXze+t2HDBpo2bWo2uALw9fXFz8/P+Ke+BFdCCFHXnbNhg+IiHTZP0znTwUTbOrfvTVDkFcP4UK1CEDiiuYafO3xz2fSZKlPoXqRTXMu1Lft1LEU/pdXBzq1xLAn2BtXLclG9KaYK8fss/hqXXhMqMT1oW/Ddysf8cS6aRqgvxFhphhqfq59qrc76KyiXwcrIyGDy5Mnk5eXh4uJCcHAw27ZtQ9M0lixZwp49eyguLmbQoEF88skngD6bdeedd1JQUICLiwtBQUFs3bq1eu9CCCFElSvRKS7auLff5WxFSwu/HJ2toKRsWwFLzqTrm3Y29qw4PncXjZEtYEOMIqNQUb5lUmVaNVzJ1v+it6agRBGZAr0CwdXFed9hsBf6pqAWNrw2p6A0c7YAACAASURBVPxnNkTraFEIHg725rI9wLL8flhDjZ+uKop1Cjcz31VsFni6QtNKdpq3V5kAKzQ0lEOHDpk88OOPPzb5uo+PD4cPH3b+yIQQQtQqMdn6PedscTkHBlbtcMo4mqxMju3MT5vY++/lJMdGERTakR6zFpESPIHbWpr/BR/WUP/PKzlwS7miaHMbDNsizsIKxdKOJkNhCQxw8pYuTbw0DiTqN3d2r0TgppQiNhu6Ozib1sQbgrxsu35DD41GHpBeaPr9sIZQrPTPyvDcyovNVrT2odqnrKWTuxBCCJucy7A9uLBlk2BnUUpxKKni9Uz1b/pu0RRcj23iFj/z5/N112jsoc/ClZdfApmFjt2bpU2VDUp0iv2Jiu4B4O/h3IDA0P28fJsIe6UV6reoCXWwp5S9m0JbmiZs5g1ervpdBUwpUYrLOY6PtTIkwBJCCGETewKsrCL9nnvVIS4b0kx0+zHVv2lw6yC8vnvJagYnxFefFTm3e0uFDuaOTBPqlDL22LK0b96JNP13N7ip8wOC4BvZuKRKrkiMvdH13dHmp7ZODxpYmibUNI0wXyp0ojdIyNXXBIaayW5VJQmwhBBCWJWYpyoEMdY22DW3952znTVTeG+qf9PI0EAKr0aZPL60Vj4aV/duYvPieyt0MP96wwarny8vIRcKdZb3zVNK8UuCooM/BHs7P8DydNUX8FdmmhMgLkvpM0eWmlSZEeQFTey8N0sZLNDXYV3J0S8iKC82W79RdvNqrr8CCbCEEELYoHz2ypYNdqtrmvCsmcyaqQ7fP8WmEFSqw7c5IT7g8s2LJjuYf/rGi3aP0bD/oKV986IyILkAhlRB9sog2LvyW/7EZkNrB1fkOdL1vXkDLDYcbdNQv3jgsoksVmy2opWPcxcL2EoCLCGEEFZFpZf92ZYNdk3VMDnb9dyKmTUDUxsp74tLZmipDt/mNPEGdd101/PLF8/aPc64G7/8Le2bt++6vhg7pArrhfT79zn++cxCRVohhDlaf2Xn9CDog6PmDcx/LtgLGrhVbNeglLqxwXP1B1cgAZYQQggrcotVhT36bNlg93qevuVAVTKXvYKK/ZsO0JKei7622PPJwEXT8GzRgV2xpve4U2a20zHHkM0zt2+ef+uOXMmBwc0cCwY6+Gv0D7b+2SZeGmkFpqfTbGGsv3IggxXgicVAyRKb6rDKLSK4nqdflFDd/a8MJMASQghh0fmMivvn2bLBrkJfKF6VzNVfGRg2Uh6/LYvixb8y/A7rwZVB2LRF7ImruMfdrfc/S4aZtgGmJOcrcm7s0GIqq7Y3Lhm3uxbRxAvaW1jdaI6/B0wM0xjZQsPH4v4slV9JGJutCPICH3f7A6U+ZjZ2toUtdVjxOVB4I6DPLVZsjlX4uEFLB4vxK0sCLCGEEBZFmQhizAUK+dkZZYreq3KaMKPQ9o2dozIUTbxMNxc1p9dt98BD/+G8V8Wu5/bUMcWVqg0y1RV9/Mtfk9BpAv2CNWNG0FauGkxu44K3m4aXm8aolpZ/rVd2JaF+ys3+z7lq0KsyAZaVa7ZpCDogLkcfXH12XpFVBDPba5Xq+VUZVmJdIYQQ9VlhieJ8ZsUgyRAo7PtkOYcOReEb2BRIpotbNiP6h7IrTl/03tT7a0Y8OqlKxmYte2VQohTnMqBfsH3nb+UD9J7IoIn30DOw7C/phFxFe3/bfnGXL/Yv3xX9SLKCOEXnRvaND2BUS41WpWqMwgPhSLJWYUrXQL+SUN1YSWh6/OWbsw6Zt5jOIyeSU6RIyoehDkxjdmms0cCBVYcG/h76FZCZRabfD/QEXzc4nab4MV7f6mJWe83uFYvOJBksIYQQZp1O1/cRMsUw/fa3fdl4+TYyWfS+6b3l6OysV7LV2fSKr5lqHRGXra/F6WTnCjZvN41ADx2XTQQrsWb6LpkSayWLdyZdv/G0vdNunRppDGpa9te4pmnc2VozEzrpWVpJaGl1aFwl6q/6ViJ7ZWBpmlDTNMIawrGU2hFcgQRYQgghLDhu4/5+5oreE6OjuF7Jxpam5BWrCoGLueBg37cb8XN3rBdSC29lso4sLltRYkOheFah+VWOoL+P6Ez72xd4uOg3qzaleQPNYr1T+ZWEpYPSrUtmMrS16dWhsdmKRh72d5hv4g2hDZ0QYFkJ7Lo11nfgrw3BFUiAJYQQwoyMQmV2C5LyLBW9V0Ud1vkMKB/fmGsdEffVS3RtjN31TQAtG+hMroYs1MFVG+q/rGW6zmXoa4c62Tk92M5Pw9vClNsfWmp4u5p+L7jUSsLyQWlxYT4jQk2vDo11sP7KGdkrsF7o3rGRxuPdXGpFcAUSYAkhhDDjRIoqs3rQUud2c0Xvg+c9S0IVZLBM1V+Zy6KVJETRK9CxX7otvfXXMZXFirYh+LTWbPVMur4Rpp+dWaHOVvpJebtp3OJn+pgmpVYSlg9KOwX6sismuUKgHBjWkYQ8+zNRHi5UqF9zVIsG+mL5ukICLCGEECadSL0ZHFiqzckvUSZXxxlW3Dmyd58lxTrFBROF9+ayaB4tOjq89UyAh8LL1fS2P9E2bN5sKcAqLFFcdGB60E2DDv7Wjwsx056g9ErC8kHp4/3bsudyCn/adLhMoNxu+mLA/gxWtwANTydFRW4umnHsdYGsIhRCiBqUX6wc2tOtqsXnqDKF0KUzHZqmMb9XKNM3H+V/Hy4nq/F45nWsuDrOoLJ735V3OVs/RQdlV7z5BDRlzzV9Fi0iNJCdsSnsjUumz7P/cvhamqYPVPSr8so+p8vZimKd4v/Ze/P4Nqsz/fs6smVLlrzK+yZntR1CcBISCEnIQqFQKLEdCE2glBnSDn2n85vf23ZYmoRpm4Vp6TYztKUFpsOS9iUQOwFK6RaShhSyk7B4yWY5sZN4X+Td1nn/OH6kZ9cjWZIVcr6fTz/F1uIjKfFz576v+7piNWwABkeprv7sVA8wShHw9uDUJGNFCxupKd978SYhK0qb8NW5ThBCcPu0TJSkJ+Jo6wAOt7iQUVSMNU8/i6aZdyOxA0iNC+ysCwyYnwZChpXg0kBkIpgmCu9gcTgcziRQVVWFsrIypNgTUFZWhqqqqsk+kgRx9wrQHr/1nGfO7W82Us1twaExpucKFY3jW33yrtrsONZmqhmz4weHXPjYnAfT13fgc3cZNxdVI99GcMENhXv7KNU3Uj3fp1be+KjpZKHJgXhzAcY7Xtk6IzVhk1BttFvb1ou7N7+MJw64sX77UZSuKPfqrwLRseXb9CNuguFK6mDxAovD4XAiTFVVFVavXo3kriY8utCJ5K4mrF69OmqKrDEPxceyAksrONmcU4wHphM09wOHW7WfM1hjSzVc46M5LVG7xZ6Mx9/rhenJw5h9SzksExxR5dmAIQ/QobINqKfD0hsPjngo6oMYD5qIcUF8rE6Gn7BJqDfaFRgeo2juAwoDzPRblhP6zmyGJfq6vVrwESGHw+FEmK2bv4+bnel4RTRuu3/3MWzbshmVlZWTfTyc7oE32kVgyfqN2PHovZLx24HGNty25ZcoSiS4Pp3ir80UJSnqa/wtgxTTDRpz6uGhvlzENlcdHl7oVHTVDh2qw/k+VhDdVTjxn5kpaJYGAYesg3KuF1ih8Ti9DcIzPcxfbFZqYGcpsutvD8rJt6l32TItBB+0UAyN6+fURrsCF/rYpmMg+qt8GzFsxBoIGUFYbUwWvIPF4XA4Eaa2rg7LZOO25YUO1NTWTvLJGGreV+JOx38ccuF95OHa77yGRbezC/MteQSWGODtRqoahByqDtbFfp/+Ss8a4ng782wqCkHQr90MxMeom3M29VHV4ORRD0WzaJQp376s6aLIsADpAXZkSgMsyAo0rA0Kxt+X8wYMU11uCmtMYOO55WHoXgEsMDqUm4TmMFZBvMDicDicCJM9Rb0wyJlSMsknA9wjFHXd6qOt0pUVePClI7D/qge5Tx1GeYVvjGSJIbijgKC+B6hRcVhvCTJcWI547KZlDXHDP27Ep51AmSPwbD81CCHItACtKuLqMSrNGhRo6mMaLb3ty0C9rwgCHylqmXM64lnheM6AR1ljgPqrfBsJSbdSDRMhii7iRCgIcOwZCLzA4nA4nAhS30Vx/UPqhcH8f9iIms7J3ZA61EIxpnOEI20sD67cSWCSXXBLUwhKkoE/nKcYlj1J64B6ZytQXH6Ck9c8/TrGysox7AGuS5vwj/OSbtGOl5HrsDyU4u+X2fe0dGKet7ZhWoCeUgV2AnuAcTrJcQSJZuX3CSEosvv0bFqMepiTfSD+Vytyw6uTCpUOy2wC8hJC8lSq8AKLw+FwIkBVVRXmzy3DnJwEvPf8Ftz05W+rCouPG4ymCQejHsqCh3Wo7aKYngSkaVzkVuQSuEeVup9hD9A9PPEzqgUnC3mI67cfRcnyVTjSSjE1EUgJcDtPj0wrQdsgVDclG0RFiodS7Dzn6wJqbV/SS3XI0/Cp0iKYMGhA2wG9KJEtJ8hd6sU097NOnFH9VaGdYJqGwWmoCNUmYaGdIEbDYiMU8AKLw+FwwoywNWjvbMKjC9iY6O8v/wiLH94oWYUHmLdSKDo9wXCyQypul2uHjv+5Cuf7WCSJFhkWaBpzanWAjNI2SNE/qn8flxto6gcWZQbuOK5HhoUVGl0qm4TN/UwsTinF7gaKT0RdSC2dWHxusaZ/lhaBjgcF8jUKuSI7s5FQG3EKuNzsvckyKC4Pl/ZKTKg6WKHQ5+nBCywOh8MJM2xrMAOvqAToyhkYYxEmk8EHLR7vf6tph954Yg1wrBozk7SfgxAyvrmmFHh/cfHcCVlR+BtnAcCByxSZFmCazhnlzEohWD1F/3IodE2O/VkpWPdQdra3GilOyOwttHRiRWs3GD8gAKedBN2R09IZpcUDiWb/VhOFdijGwao/x6YdzxNKQrVJWBSCAGo9eIHF4XA4YYZtDaapBuiqodb9CTdneqgk0kZLO2R9ZxtsfnRABTaC833Ap3+tkhRpWe6J+X35y/W73E9xugdYnGVc3J4aD9ztJChO0bcVSDQDscerceB7a1QF6282elTHq3Kd2ClLHvD1HZh/q7YtghrXT8ARPUfDcNSrw9LoYHkoHRe4G/vZZSHKHPSHI575gU2EOBMCHtEGCi+wOBwOJ8yUFBdjn4adgBrnDWx2hZoPWow5t480qxeFYgrszL19n6xIe2XVPCx1pmPbls1BndHl5335ewtFchxwjUFxewwB7pli8kYV3ZFPoGUxRQhB7NtbsbRQWXQeeGELeke0f45YJ3b9fx4BmVuumROoRkJs8PorADCbCLI0shiLEgku9gODKjqsS+OWGEb0VyYSuIVEsMSYCNLiJ/YchXblkkao4QUWh8PhhJkNm57E31TGRIvXb5LcTxin3Vdiw/y5kYvPaRukOC2zZtBybndoFIVi8hKYpUBHY72iSAvW76tnmKJLRyTfNUTxUQdwYyZBjMEL5+fyCPJEAvA0C8FNWdqPHb1Yj+VOZdGp1YlUw+WmyE0A4gIwcypzkID1WnI0dViJ2josl5sFS+ca2LSbmkiQEMFMzYnqsIoSQ3QQHXiBxeFwOGGmsrISG3+lH0ci1jw9utAJe2fk4nM+uEwVmXlq2qEDjW1Y9tVNqs8hJi6GIMsKWHJnKoq0vY3tKC0J3O9LzxUdYB04Swwwz2Hs+YqTCRZlKS+BS3MIUjQCje35M/Guy3gnUg6lFC43UBiguHp++sQLFy0dVmockKShw3K5KfJtMLRpNzstshE26RPcJCwKo/+VAI/K4XA4nAiQeVMF1s/X1t2INU+RjM8Z8VBFsDPg0w4deGELDh2qA80uRsnjv5AUhXrk24C6uzZg/3/fJ4nXea+xDVU//XXA59TTX/WPUhxrBxZlGusMEWhH6JhNBJ/PN+HVsx7Fbdc+uBH7v7tG8XrWPP2sodfQNQz0jhjXNAHAlEQCRwi25rQ6WIQQFCVSxQIBHddfLTSwjRlLgJLkCR8xIFgHK7hRenwMkBtm/RXAO1gcDocTdtoHKdpV1vvFqGmeIhGfc6bHFz0jR9AO3f/HXoxsOIybbjcuzC6wEfReW4Hyp3ZIOndbX3gdFRWBCbwBff3V4VaAUmChQSF4oZ0gUSUvUaA0lagWQdffVgE88ipq47Q7kXoIXbhAOlih6F4BQGo8gU2jpVJkJ7g4AAyM+t7jum620WpEfzUtiXh1bJFiIpuEkdBfAbyDxeFwOGFHK3pGDNM8NeGrc1l48UTGaYFQ2+U7W82earz3/Ba0ueqQ7izGkvUbUbqyArVdFLbYwLauhKy7hIUVWH+rrwNXGsQoaWCUerMM5Wdc/PAGHM8qx5w0+N1u9J7BgGB8WpJyLJloBuKvr8B1d1ViSXbgr8PlpsiywnBYs22C4nY5eTaCepU/i4IeyeUGSlKAui6K189RzEw2VmBFejwIsBFhsD2scPtfCfAOFofD4YSZ+m7/95Frnu7fdQz7XW3YsOnJsJ3LQ6n3gquXmVffDcxMNuaFJJASxwqE832yyJwgPL7Ou9mFVO2Mrz22Bt0fVOOaVONnKzVwXzX3c0IIMqxA62BwoymX27gjOsDE7aF0GtfSLaXEAclxQIOb4tNOih1nWXG1Zor/To/ZBBRHeDzIfi5BSpCbhOH2vxLgBRaHw+GEkcFR6te/CVD6JZ215uHX23cGNU4zSqMbXmd0Ld+rfc9vQfsQE4UHAiEEBTZWHIlpHQzcqV4YD2qd0fz2NsNbYbkJBMk640GBPBvrkMhJt8DbTQuEnmGKziHj+iuC0I0HBbQ27wQ/rI86gNfPUcxKBe6ZYqy4m5FEAtqIDCXBbBLGxzBfsEjAR4QcDocTRk71AB6D9UTpygqUrvQVVLPzw3vhEo8H21x1eHihU2FB8MHBOsQSYGoAzugC+XaCfRcpPJR6OyEjHib2Tg2g+yBYCOid0Wh3zejILT6GINMKNMiKqUwLwccdrEg0amYKBK6/mplMNPMeg0Uvw68okeBEB8V1acx41ej7ORnjQYEMi7HusJhI6a8A3sHicDicsKKmeTGKXkZcKBAXWFqZeTE5xShNYSOZQCmwsYLqsqxIaQmgAzTioWjup5pn3ONqR0qhMZsEwNh4UEBtTOjNJAwwuNrlpnDEA3aDOrHluaEvAvSE4XPSgLXTCFYFUFzFmYAZkzAeFMjQME/VY0oE/K8EeIHF4XA4YWJ4jEqKmECR65dCyaV+qXGnVmbe0Bc2BB2BkpPAHL7l0T+BaJia+gDBZFzbm2ujoefKsADpAXSF1LyjhCIlkCIRYMWy0+DFfWYyQU5C6Aus+BiCJLP6bSZCMDPZeMQQwM4ZTOEdKvQ6clpEwv9KgBdYHA6HEyY+7WIdHC1GPBTH2ylGNWaI7hGgcyg8RVZtF5UEMb/3/Bbc9OVvSywVpjy6A8k3lAf9r36ziSDHClyQadCaA8haFNszyHVqh0keCr+9A7NvMaZTK00J7OKq5h2VZGadm0DE+n0jFK2DxvVXy3PCVwQE0/XRIhLBznpkWNR1clrExwDZEdJfAVyDxeFwOGHjRLt+cfTnJorDrUDXELBCYyR03h2YXskov3uNBTEvdabj4YVO7G1swv6Xf4Q1P3wNpSsrMDxG8eOPKG5MQ0BdDTn5dqCuS/q9BrdxDZN8TCro1DoGKf77U4obpoR2e1CMIx6wxEi/RwhBhoWOd+GMPV/jeEFpZINwZjJBrspoMlRkWJj3WSiIlN2BFnExBElxQLfBcW0k9VcA72BxOBxOWHh5x05s+OI8PLXYjufWzUfNnmrJ7Wd7WHGVbQXeuwy0DqgXY+EYE3YNUbz1C/WNvAMvbAEA1HQxA9Jgx4MCBTYy7mDuex39o8ZGbB5KNYOva7qYg/h0g+L7lDgEPHYjhCDHqmxBMqsG48/jctNxKwT/P39ZGLtXwMQz/ASS4xByEX4w5AbwmUa6IOQFFofD4YSYqqoqPHjfPZiq4ikFMOuG3S6KKYnAPxaz7Lu3zqvbF8gF4qGgtouqOseLg4s/bKcosjMH8IkgBAVf6pd+v8GAdcWlfm2X+ZouihnJxkOTSwIcDwrkWZXnzLAQtA3C+3mJR63PrZuPP//n45Kv6/ZUG+pezUiWhk+Hg4k4oIsJJO4nnAQS2hwp/ysBXmBxOJwrlsHR8InAJ8LWzd/X7Q69c4FiaAxY5WQi4bsKCRrdwLF25XO1aHS2JkJtt/bWYEZRMTqHKBrcE+9eAaxzZFbRLDX0qt9fjFY8TvcwRVN/YEVToPorgVyVDlamhWnrOoeU5qfJnWfx95d/hGmDvuK667/ug/nDXX5/1rIg3OEDJRhhuBqRLla0MFroRdL/SoAXWBwO54pkeIzi+TqP5mhtMqmp1e4O1XRRnOgAbi/wGV5OSSQoSwP+0kThHpG+nsExoHc4dK9xaIwZn2ptDS5evwkn2iniTMY8o+xm4PM6fl1Ms6TcHGzo9W84qmVTUdvFthNnGrQISI0PLP9PTI6VKpRWQmRQY5/S/DQGwNICB7aX+4rrxYXpOL9jq+7PcdoJ8iPQFbLGEtg1NgkDIZJ2B3pkWQFrjP/7RVp/BfACi8PhXKH8tZmibRB4wxW4M3i4yZuq3h1yFBXjrUaK4mTgujTpY27NJzAR4I8XlK+lJYh4GS0a3cz4VL6RJwQXlyxfhQ87gGtS/Y/fCICKIhMWZZl0OxoZKu7nA2P+x59aDvg1XRTTEgGLwfHgXEdg9gNiLDHKsZo1liDbyopE+aj1TGc/lhelS4rrlU4HOl11uj9ndpruzSFlojqs5LiJj45DBSEEhQYK08kQ5PMCi8PhXHGcd1McamEX3/N9FB+0RE+B5aEUix5W7w5l37sRg2PAXYXKC35CLMFteQQfdwJN8vy+EHbpGmW2B+u3H8UTB9xYv/0oSleUo8HNtrKMjAdvyiKYNr6qf2cBKxDVyLAQtIo0SwLnerVfV9sgRd+o8vsjHorzffD+XH+YCCuwJoKa4agzEWhwK0et01ITsLehTXX0qnfGWUGOMINhomPCSHpJGcGIDmsyRpq8wOJwOFcUox6K3S4PxJfmPc0UHUEG8Iaa093A1GXq3aGua1ZhSqK2m/e1aUyvJNcnBROQrIXLjzv8kVbmOF6g4gElJs9GsFJkLZFhJbgxU/11ZViZWL1nRPp9PR2WS+O2pj7WgTPakZiZTJBoYHtPD7UCq8hO0D0MzP3KBkkxPUqB/efbvV+vFY1etSiyE9gMOryHgokK3aNFfyXg7zyTob8CuA8Wh8O5wth3kY0GxYx4gDcbKR6cMTHPplBwooMVevJcwYFRitdOUtxRoH0+EyHIS6C40Cf1WArUNVyLUQ9Fs47tQ8cgRU0X8IUC/ZFafIx6GPCyHIKPOih6ZYWU0DFpHWTjJQGXjh+W1njQ5WZju0yDRUIoApPVDEeFrcDY61kxfeCFLTh0qA4ZRcW46Qu34tzBP+PgoTp4sorx+a2/ROmKcs3nvyZAf66JwkaEwf+DJJDNvUiQbWV/JgbH1G8vsEVefwXwAovD4VxBXOqnOHBZ/cJwrpfiaBtwfcbkFVg9wxQ1GtE4p3vYJa3YjzC7wA4ca4Ok8AgkWkaPpj6Wo6fF+y0U1ljgOof+89yeb1LV4MTHENyWb8LOc9LNu5Q45lnVOiD1rRocAy72A7kqBYzWBqHLTeG0Gyukk+OM+2TpkWFRXsCZDouioZeiXFZMM/4De5s9ONQK3DBHr6gGSlMnfsZAmMiIMJr0VwKCDksr93OyCkI+IuRwOFcM71+m0EiVAcDGW5PJoVbt89V1U+RYgSQ/46oCG0HfqDRMeHCMFW8TRasrBLA4lw/bgRsy9PPlzCZ9Qfa1aUQxstHaJATU/bB6hqlqmPKYh+K8G4ZEzQAwLz14cbsYQtT9qZyJ+iNXlxt+i8EpiQQJsZEtWGxmgoQg2yvRpr8S0PMZm6wz8wKLw+FcEVBKcapHv8hoGWQ2BJPBiIfiqEaBN+ahON0NFBsQMgvjqPOyC3codFh6xcChVgpCgAUZ+s8xPcl/wO81Kh0ZLfdzuQ6LUop9F33vo9jE81fr5mP0qDHTzlCI28Vkq4wki+zMpb5LJS9y1ENxoc9/MTg7wuNBgWA3CaNNfyWgda44k3qHNBLwAovD4VwRXOhjESsCcvfsmj3V8FA2BpsMTrRTDGhoQBrcwJDH/3gQYKOn9Hi2HSl+jXcunouqqqqgz0cp1YzdGR5jsT3zHOzn62HEsFNNFC5sEn761yrJ5/b73VXwiDbuft9IcbTN55AuNvEsGW4Gnr0PXR9UK55fzowk4rdbGAhqIclCodegUrg297NxrFNnPBVDgBIDXmPhINgxYbTprwRyElgxJWcy/K8EFMe57bbbMGfOHJSVlWHp0qU4fvw4AOCdd97B9ddfjzlz5uDGG2/EiRMnvI9paWnB7bffjhkzZmD27Nn429/+FrlXwOFwrgrquqVdDfGFVxxFozcGCxeU6ltF1HVTJMcxU0Qj5NuB03ulrzGrtwmrV68Ousi6NAAMaRSAx9vZGFJrC1AgxqC5Z5ZVebHLsABDh6vx2mNrJJ/bb799L174XZW3uDrS5nsf5Sae28uZI/77/6Nv2gmw8WAoUStIBD8sl4rdhMvN3gO1zpfA1CTit6ANF8FsEqZEof5KwKThhzWZBaGiwNqxYwdOnjyJDz/8EN/85jfx0EMPobOzE/fffz9efPFFnDx5Ek8//TTuv/9+72Mef/xx3HjjjTh16hR+85vfYN26dRgZGZE/NYfD4QSNWMCquPCKomjOT0IH63QPFJuNApRS1HWxwsSoHqjARtBbvVXyGl9ZNQ9LnenYtmVzUGfUKjw9lOL9ForZqUCKn4tnUSKBxUBBYCIEubIuVoYViPn9VtXP7cdPbcFbsuIKgN+8RC3S4oEZBl3ejZJhgcLRHfD5YclpdFMU2qHbPYn09qCYYEaE0ZI/MdGgjQAAIABJREFUqIV8dDwzmeAGP/9oCCeKAislxdev7O7uBiEEZ86cgcPhwDXXXAMAWLp0KRobG3Hs2DEArCh75JFHAAALFixAbm4u9u3bF4nzczicq4CuISqxKtC78F7oi7yz+wctVHVkCTC38p4RoDjZ+C/6fBtAL9crXuPyQgdqamuDOqOWr9QnncxY9KYs/+cLZJwltzZIiVN/TSucDpw7VesdC4rRy0vUY0m2KeRjobgYgqQ45ffVdFgjHopGP2L8WAKUhLgIDIRgRoTRqr8SEJ/v+nSCL03zrxcMJ6oarAcffBAFBQXYtGkTXn75ZcyYMQPt7e34+9//DgB444030Nvbi4aGBrS3t2NkZATZ2dnexxcVFaGxsVHzh7rdbvT09Hj/NzQ0FOKXxeFwPkvI16/1LrxDY6HzjTJCywDFW7uqNEeWdd1AvCmwqI4MC2DKnom9Lulr3NvYjtKSkqDO2aihvzreTlFkB7IT/MfilARQJBbIOlgmQmDOUb4mvYJJnpdoxLQzJU4ZQxQq1Lo+ch3W0BjFb0+zbVK9AmpKkrFuYLhIjCOwGMjwExOt+iuB3HEd1i25BHc5Q19kB4rqouZLL70EAHjxxRfx2GOP4e2338brr7+OJ554Am63G4sWLcKsWbMQGxvcnueyZcskXz/66KN47LHHgnquK4nOzs7JPgInAvDPOfQca47FwIDv34MLH/w37Nr4ANbtOoYVTgfedbXjvcY2lG/9TwwMDODji6Mwp3p0nnHiCJ/zX5pjsP+573tHX4QQfHWuE+t2HcP+5zfDs/GLmGKnGB5SyX3RIa3yO9j/sy8pXuOL3/8xOjo6AnqujiGgtScO9ft24+CLP0B74yk4Cmeg7IHH0JBzL27NHsXAgP77lWulGHGPwOhPThgFBgakLZ+M1U9g/0/Wan5ucooW3Y7yLa/g4Es/xPsH6xGXW4zyrf+Johs/r3p/AFiSPIrurtB99uK/z5aRGAwMKKuSzHgzznRRFMaP4vVGM9qHCdYUjsAOCo1jIi1xDB0dGqK4CGEdM6NzwFgRkmSmoH0j6JikJRKj3J1JkBtPEeBfEQDB/+5OS1Ov6HUrpK985St45JFH0N7ejhUrVmDFihUAgKGhIWRnZ2PWrFlwOByIjY3FpUuXvF2shoYGFBYWaj7vvn37UFZW5v06Pj4e8fHxAb+oKxGtD4Lz2YJ/zqFjeIyizeOBVSTKve72LyEuLl7inr3m6We9btk9JoK0tPAvSSelpMLV6EF74ymsX+hUjL4OHarH8KAJi7MJrFZzQM8965bV6Bx+Fef+us37Gp/6n2fx4IOrAz6nq42i4f2d2LXxASx1pmP9Qif2Nl7EH/79AeARM+Z8pRJWP1Et1+cF9p6mAchtGUOnaEAxc+U9uDxokrwm8eemxnW3fwlzPn8fnj5JsSADuC5X+wxJZmD5VBNiQzwWEv4+T/VQfKxSiE5J9uDTTuDVxhi4R4GHZhLkJOjP4ObkmZCmsm0ZSWa4Peg06B03Ky0yf6cmykR/84byd7ekwOrq6kJ/fz9yc3MBALt27YLD4UBaWhouXryInJwcAMDmzZuxcuVKTJ8+HQBw77334tlnn8V3v/tdHD58GE1NTYoulRi73Y6kpBDY63I4nM88Z3rU3cflUTRitMZhoeZcL9vMYyPLJnx1LiuyhNFXQn4xRklwbuIFNmCkrALl6yqRPj6amhXkZpzLTSWLAeIu29F3tsG23n/RVhJEGHG+jaBTpE3KsLDXdN+XK5EYQPZe6yAwMOZfZL04m4S8uBKjpVsqshMcHN8ifWgGUbV0EDNZ2XhynHaCQwYLrGjXX0UjkgKru7sb9957LwYGBmAymZCRkYG33noLhBA8+eST2L9/P0ZHR7Fo0SK88MIL3sf94Ac/wJe//GXMmDEDcXFxeOWVV2A2B/avNQ6Hw1FDK/5Cj84hwD1CNUOVQ0XteCzOkvUbsePRexWjr+R//QWKk/17S6mRJzIcTR+/sAerLWt0U7S56vCwSpftg4P6G3kAy/1zBLF1VmADPhKNaryZhANAYgCXCJebCYbVMgEF7ObQWzPI0SqwpiUBizKBBRnEkI3BZHozidHz6FLcNwANIYchKbCcTicOHTqkesfnnntO80mysrLwpz/9KbQn43A4Vz1i9/aaPdV47/ktaHPVId1ZjCXrN6JkRbm3YJDfnrZhE/5PEOM042fzeXOVrlQG/t625Zf4U+Yq3JEW3IU0PoYgy8rMQeemB59J2DtM0TGk3mXb42pHup+NPCAwcbsYZjjqO3NqPPPSah0EpgbQ1Wt0U+Ta2CafFjdlhX9jzBJLkGiGIszabCK4Ld/4zw5k4SGc2M0E6RZtixGBJDOQFqTz+9UMD3vmcDhRS1Mf4B7xGYsudabj4YVO7G1kW3qWf34VqyoqQY/tUtz+r1+5B/n2naisrAzP2QYI3KILrXxk+YfzHtg6gel+VvHlIcJi8m1Ao8hjaWgM6B6mSA7AoVyIx1Hrsh1obMOap5/1+xyzgvRryk5g2YUj47IlEyFIt9DxQtHYc1JK4XIDc3SkMZYYtpYfCTIsBL0jExtBR9O4zWknaPNTuEfTea8kol+xxuFwrlqE8aCasejiwnSMvLkNvz9Psf/5zYrbJ2LKaYRTvdq/Psc8FB91sKIgRmcUlG4B/qnUpLkuX2Bj8TIDIhFaa4BjQsFgVOiynbPm4QeHXPgAebj2O6/pCswBIDeB+LVw0MJECHJlj82wBPYaOodZx0hPf1WSQnS7W6EkGAd0MdGivxIwYh7KC6zg4AUWh8OJSiil+KiTFQdqxqIrnQ6YLtdheAxobQitKacR6nUKrPpuJsou0wkbjiHA6ikmpMYTrMhVv5+gORLnKwY6JnSJHNxLV1Zg/fajKP99L0Y3HsayL+gXVwAwLz2gH6dArpvKsBC0DMKwGaxgkFqgM1aLpCN6sBl+AtGivxIw4m3F9VfBwQssDocTlTS64V3x1zIWzSwqxopcApo1U/X2aTOCM+X0R8sARdew9kXyw3aK3AQgU2eb7JY8gpzx7s6CDKKaWZcWz4wTW0QamUCE7gOjVPX+NZ0UmRb/wvU4E3BtkBoyAXnwc4aVjUT7xm3BtBzwBVxuimwrYNHoUFljgKkRNMAMJmJGTLTorwSS4ghSdVySkszBLThweIHF4XCilJMdvg6HnqP3wgwgsXKD5PZ147d/6V83huVswvagGu4RilM9+t2r6UkEi0QZaSZCcGehSaFKIoSJkFsHfD+vud94B+u8WywxZ4x5KOq7gVIDsTfXpBLET3D0Ju88ZY53gC4PaId2P3P3dG/BdWpvte62W0kKQUwE41AmOiKMxjy/Ip0zOfl4MGh4gcXhcKKOUQ/FJ53S0daaH76Gk7F52HrQhdOWPKx5+nWUriiHiRCsWV0JPPIqPjYzfdE5K7t96jJ1n6yJUtulfdvJDsBEgNmp6rfbYoHyIqIIfi6wE1ynUpRlWNjWnUDLANCvZgymgksl4PlsLzDkAUoN+FqFwvbAbiZIERm6p8WzrlOjzJvra/OK8MA1eQCAWbG93oKr/5n7EHN8l+bzRzowOSGWICHI9bD4GCBXx2pistArovSKL44+fIuQw+FEHXVdys26rJvK0ZuyCityCZZmS3/p5yQQ3Hh7BY7Mr8A/zyJIGfciujwQesPRnmGKixpdJEopPmynKNHxvro136Tpz3VrHkFtF5W89gwrQU0XC7AmhICC6ZJKNQo4MS638ns1XRRp8czbSo8MCyv6QkG+jaBrmL1nhBAUJbLNQLk31zNHzmFpoQPby+d7zVDX7jqGs7/bCnxRuQ1qjQnM7iFUZFiIavHqjwJbdOmvBPQ0VtGePxjN8A4Wh3OVcqjFg0GDnZBI8/zvdip0OTXdTBh+Y6b6Y1bkEJgI8LEoTqxriEXthJLaLqoYuwk097Nuk9Z4MIboBwDbzASLs5Rbd8MeoEdkCdFg4OI+4lEvBM/1Mmd5eQdNTihNO+VdG6ed4EIf4HBKtXOnO9xY7kxXLDO0u9TNUEtSJqdgCVboHq3FSmo8QXKc8vuJXH81IXiBxeFchVzup/jjBfa/aOO3O3biZ99Q6nI+/HMVpiVB00wyLoYgzwZcEMXkUPg3UQyUum7t2w63UiSZtbsqRYkEFj+u7lMSlQUWILU2aOj1/7ldcAMf/1UqID/6pyp0DfvXAcUSqI4rg0UuDHfagTEKFN8v1daZY0x4t6FNsayQoWGGGunxoECwOqxoHrep/ZmIRr3YlQQvsDicqwwPpdjtohijwPF2irM90VVkbd6s9LRaUpiOzqptmOnHUbzABpzvk1oAtIawwBr1UM3ipmeYeV/dmKndVSkxICzPSWAFjkBKHDPrDFSH9fKOnQoB+VvfWQMcq/a7dl+SQpAQRLyPFvKOT5aVmYPSueUSb65YRx7eO9+uuswgJyF2csaDQHCbhNGqvxKQe13ZYpk7Pid4eIHF4VxlvH+Z4q9vVXk7G8tvmIsdr++c7GN5OXdK6Xm1wukAvVSHmX5c0fNtBP2jPnsHgFkqhIqmPtZ5AYD6fbsl3aHf766COUbbN4rAWORMjMln3wCINgkHpZ05wR9Ki9/8RGnOuqQwHfF/2Aabn4zG+SF2RU+OY5YPAoQQOO1MIyZ4cz1xwI1/eeOMt+D6j3Ez1C8+pW6GOlnjQSC4EWG06q8ExEW3Ix54uNiEXFv0nvdKgBdYHM5VRPsgxc9fqZJ0NooGmnDfvfegqqpqso+H1gEKh4bnVVxusd/wZsHU8rzEmDN05xOEzTV7qrFr4wOS7lD9D+7DlFPVmrYGeTaCRIMRN0pzTuXr0NNheShF81n1QnXson64c04CMCUptBdWViRKn7MokemwRj3S1yEUXHP/vx4kbT6Mebeqb4JO1ngQABLjCKwa7vtaRKv+SsBhYTmLhXaC9SUmnj0YAniBxeFcJVBK8YaLYt9z6p2N734vfLEyRjnZQRWeV4Kn1Yx1/j2trLGs23O+TzwiDF0HS9jK04ru6arapvlYI+NBAfn2XoaFoHVAOvrU02Fd7Fc3Z93jakdSoX6489Ls8FwW5F0fpx0YpVKXejEut/Z2W0IsMGWSCxZ5weiPaNZfCdyaZ8KDM4jmBiwnMHiBxeFcJXzcyTowarEzK5wO1NeHL1bGCKMeimNtVJGZVxuXB3x9B5bc4T/WBWA6rAuii3bXENuomygeSr0Ceq3ong6NbTeAjbSMotbBGvIAx/7sE61vrZiP372mPtp1udUL1QONbVj8sHahmmExZkAaDHJheJaV6ZIaVKwkBkYpLg9oi6yLoiBuxp/NhZg4U3TrrwTmOAhiI2ja+lmHF1gczlXCp+PGnVqxM+nOYgyF2NIgED7q8MWniHU5zh8eRuoN5YZ1L/k2gpYBeF8LReAByWpc6geGxv2ptN5DrW23dEtgHY+kOIIks+/rDCuAY9V46ztrvGPJqUNNWLdGfbTr6lWGO5+IzYPtX3ZojtwA1r3yZ98QLHJhuMmrw1L+mRNGvFodrGgYt8lDrPUoiIKCkBN5eIHF4VwFjHkozo6PlLRGcIvXb5J0fiLNBy0exfcoZbEuxSn+fZsECmysqGoKsQ6rUVQI6L2HagTSvRIQjwlT4oDY32/F0kLpWHKpMx3btkhHu5RS71nFhWrC9w5jxjLtLmBaPDA7LeBjGkatQHbaCc67lTosVy+zu0hR8WYClBtvk0EgRV605Q9yIgMvsDicq4AGt6/7Iu9sCLEypSvKcT4Id+pQcLaHjYTkXOwHekeAYgPbdwLpFmYBIBa6h2KTUOyKXrqyAuVbXkFtHIvuqYv3vYdqGImlkSMeE5oIAb1cj+VO6VhyeaEDNbXS0W59NzAgc8H3jtx0CpPFWaawdllS46X2EwArPEYpM2gVI+iv1IrqhNjgjT5DicNCYDf7vx8QHQUhJ/LwAovDuQqo71bf1HrigBvrtx8FKMVz6+bj9mk2zJ9bFvGNwg9a1Augum4KS4wyMFgPQgjyZYajoe5gAcDMZauQ9dRhZLzQg6//9ohmcZVkBnITAv95+bIVeUvuTOx1KceSxcUl3vvUdlG8dlbZCWwcLw61Rm7JcUCZI/AzBgIhBA5ZYZSdAMSbgIZxywlKKd67RNHUzwKx1XDalTmOk4URI06zKbjPn3PlwwssDucqQF5gianZU+21bXh0gRP2ziasXr06YkVW+yDFKY3z1XUDM5KAmAAvqPk2AtffqvHcunl4arEdj39x3oReT9sg9erDBEY9wKluYFaK/gW/2M/tWuQksGgdgSnrNmJ/o3Is+Q/fZGPJjzsodpz1QM1/1OWmSNQZud2URRATAXGzmg6rcFyHRSnFn5oo/tpMsSyH4FqNcWU06K8E/Bm2Asz/KhLvLSf64AUWh/MZp3WASow35cgtB17R0PaEiw9a1LP9uobYWKtYZ7ymNaKhR6sx8ov7MG2o2ev1NZGiUc3Us6HPhGGP/627WUGMBwEgVmY4es3KCuCRV3HaIh3tTl9ejuNtFDvPeaC1LNnoZuM4tUIvzhTa3EE91CJmnIkE5/uA6gaKD1qAOwoIludoF6XRZHdgpIMVTQUhJ7LETvYBOBxOeNHrXgHMcuDhhU6FtufHR8Jv2zA4SnGiXf18H3Uyzc50jTiURDPwtRITXjzlUeQN1v92C5aMC8IJIfjqXCfu330M27ZsRmVlZcDnlI8HAaC+14R0C5Bh1b7IpsVP7AKbL7KcyLAAmFeBW9dUYqpofHaineLQmHYA9fAYRXO/dgB1cQrRzHcMNayDJT1pkR34iwf4pAtYXUQwO037LAmxgdkjhJtMKztT/6j2fbj+6uqFd7A4nM849TrhxIC65cDexnaUlpToPzAEHG2jGFZKhkApK7xmpULVGZ0AKC8yITGO4AsFyl9j7a56rDAgCDeK3EpgjFKcdpv8dq/mOiamFyoQ6bBS49nIUK4nGxyTlyxSzvex27XGWcF22IJBTZyekwDMdQDrpukXV0B06a8A9ueqUKeLZTYBeVx/ddXCCywO5zPMwCiVuJoLjHko3nR5VA0p7999DPtdbdiw6cmwn+9Dje7VhT6gfQi4TuOCe2MmwbTxLs7UJKKITUl3FisE4cEWjd3DFF3D0u+5eoHBMaK7HWgi2l0jo+SLiiKTSiahEVxuioRYtl0pJ84EzPCT7xhK0uKlujKAva67nSbv56lHNI7b9HRY+Vx/dVXDCywO5zPMqW6o6nL+doniWDvTvUxbVi6xbThrycP2HTtRUaFtSBkK2gaptxtTs6daEpz87ttVSI5Tj0PJtgKfy5NetD6fTyRhwkvWM0H4WpEgPNiisVHFabymiyLZTJGtM66anmQ8e1CL5DiWDyeglknoDz3Lg5nJkXXujjERpMUH/3gjmqdIo3emaCwIOZGDF1gczmcYNf1Vcx/F/ktsLNM3Auy7RBW2DXM/F97iCmCWAoB0i1EITj73w/uQV1utKArMJmD1FJOiK5AUR7Asx/e90pUVuGPrDnyAPPzHuCD8H3/2elBFo6vXd06hCDzxf69H+ic7dcdVoRKOi8eE6SqZhHqMeCia+rSLgMkITJZvEhrFGsPidaKN7AQW+aNGNAnyOZGHF1gczmcUD6U40yO9EI96KKpdrPNyZyHBzTkE718GLvdL76c2Vgw1tV3s/7WCk9tVgpMXZxFNUfmNmUSi8VlwWwXi/v0wFr/ei/Xbj6JwSXlQmYSNfVRRBN6AZpz50TrU7KlWfYzdDMwM0ehNPCbMtDATUbllhIC8E/j3P1RjjKqPsSI9HhRQ2yQ0gjMxuvRXAiYNHVYsAfKugPxBTvjgBRaH8xml0e1z9BYuvD9YkoiuDddjjmsXYgjBTZlMm/NmI5V0Rc6rjMVCiXuEoslPcHKnSnDyLJ2OS4yJYEWu71caIQSZVjaKBILLJOwfpWgdUBaBvy2fhyWF6TjwwhbVx5U5Qpc9lyUqKIXiRG1MqNYJ3Pvva2D+sFp18y7S40GBYDtY0Rw3o1bA5tt5cPLVDi+wOJzPKDUqI7gnbnBiEWnGHzeuQc2easSYCO4qJGjqB460+R57sZ8q8uFCSW2Xz1bAaHByWjyQqWOJACg1L+kWoEVUVAWqX2p0s8JMrQhc4XSgtUFZBAJsezBUiLtywibhpX7l/bQ6geY/bFMt9iZjPAgEH3MTzXYH8hGs3Qzcmhe95+VEBl5gcTifQcY8FB93sIJFceGVdV8K7QTzHMBfmyl6R9hj1PLhQsnLr1Z5R1kDvV2SLca1GsHJRgKTE2IJHCIRdYaFoH2IjUuBIDbwxvVXRotAgF1sHUF2adRIiiOwjGt8YsZjgNR8ubQ6gcNNyiIwzgRMn4TxIAA4LGzDMhCiVX8lkJvA9IEA88ZaX2xCno0XWFc7vMDicD6DnO7x6XSMdF8+l0dAABxqEY8Jw9PBevX1nfivf/GNsmbHMSfNmjE7th104USMenCykQILkGb4ZViAMQqvk31LgCPCxnGTT7mVxTqNIhAIjyt6uqhgc9rZZiClVKK5MhGisKbY42pHmlNZBM5Mjpy5qJxYE9GM7NEiWvVXAjEmggIbsw55uNiElPjoPSsncvACi8P5DHKyw1ccqXlCybsv1liCKYnMlFJA/N+h5Pvf26wYZS0pTIcpIRljz/Sg8n+Uwcm2WKDAoGBYIgqXaZYu9RsvGofHKC6O3790ZQXW/PA1HDXlYdtBtpVYvnW74pxJZmB2quEfYRjxWK0okWBgDPjgj1LN1ZTEWEVW4YHGNtz81Y2K59PTskWCQHVYRjL/Jptb8gjun05UjXE5Vyc8KofD+YwxOEpR1+UrJBY/vAGvPbYG63YdwwqnA++62vFeYxvWPP2s5HH5NoJ3mynGKEUMIV4Reqg5c7oOqxc4FR219w/WwRGvLmYOJDCZ2Rqws9tiAUsMK7BKAPSMAJ1DFKkGOgzn+6QeYtOXl2ModRU+n0+wMINgYEDZDgtXaLJ48y7fxkZsh/7HN/oV4oBu++37ONLSh0OX3YjNLUbWN3/BMgxFxJLJ2R4Uk2FlQd5GiWb9lQAfCXLk8A4WhxMgVVVVmHVtGWwJVsyfWxZ0gHC4+KSTaagEMhezkODaOGlIsLz7UmBjj7s8Xjf0jjAn+FAy5qHI0NAzIbsYi7LUC6niZOMXr0wrvKajhDDrhtYB3+toUAluVkPQXwlcGC+4tLoptlhgfkZ4LrLijo/ZRJCXAPReqFeMfu8tzYWHUjz2Xi/GNhzGrBVK368C++SNBwVy/CwriLHEQNfQlcOJVniBxeEEQFVVFVavXo303iZ863on7J1NWL16dVQVWeLxIADUdQGx8yvwtd8e8RqJyosrgGXCxRCpRUOgW3f+ONcL3PSwup4p9osbcF2a8jFxJmCaRuCzGiZCkCvTYYlfR0OvsaKxUTYidfVSWGOYF5UaN2SGr3CRb9457YApe6am8P7yADDkAZwqTuLRMG5TO5cWhVGWP8jhGIUXWBxOAGzd/H3cLNIPvbJqHpY607Fty+bJPhoAoGuIKjbM6roppiXB78U/1kSQkwBc6POJp4szEkLapavtol49kxDNc8aSB9PXd2DR7RWqvkHTkgL3ExLrtTKsBG2Dvk1CeXCzGmMeiguy+7ncQKFG5Ex8DLAwTN0rAEiOgyQKyJlIMPKFDZrCe5ebFctqQcPRMG6zm6XbnnrwuBnOlQovsDgcP1RVVWH+XDYSPHHyIyyTjWWWFzpQU1s7yadknOjw+UsBQN8IxYU+4yO2fBtwdp9PPP1vC5ywdzaHrEt3ptcnGheieYp/fASmeeVYkKH+GKPbg2Lkm4SjFN7A5q5hVojq0dQvHbOOetj7qBU5syCDwBIbvsKFECLZJCywAWReBeZt3OEtVMWj30Y3Rb4NisI0lrDPOBpwGiz0eNwM50qFF1gcjg7CSNDe2YxvXe+EzWzCuw1tkrHM3sZ2lJaUTPJJGSfbpYXDqR4m9zYqai6wEQzt3irZ8ntl1dyQdOl6hqnXLkFgxENxqJVibjrzsJJjIsFFzoiLCGG8JnZxb/DjVC/XXzWPF1xqoy2zCViUGf4iIF00JoyPIchNAIavK5dkSJauKAel1BvwLCea3MWNBDfHx7CsPw7nSoRvEXI4OrCRYAZeWTUXhBDkJ1rwT2+fwNrqo1hZlO7dyPvZ93812UfFBTdFu6yAqe9mnQy72egGHkAv12P5DU5Fl+7HRybWpXOpFDUn2oH+UZYjqIbTTmANojNkMxOkxQMdQ0CiGYg3MR2WYEzh6qUo03Fbl5/V5WYjOjWx9VwHgc3g+zsR1HRYH3WyIl88tmwbZO+pWgETDforASNnKbSHLnKIw4k0vIPF4ehQW1eHZYVp3gvYF2Zk457SXBxs7pSMZUbLVkmy/CaD47Lu1aiH4nQPM5VUQ3AHF5MYRxCbPVPhmxWKLp1cG0YpxfstFKUpQJqGbcKsCXhKCWNCQggyrFIX9wYdHRalVBF27XJTFNqherFflBWZAkAecu1MJOgdgaIr6HKzX+xqo8BoGrelxBMk+zEcjeb8QQ7HH7zA4nB0KCkuxr7GDkmxcalvGFnTr/WOZUApfnLvfNgSmCD8zTffjPg5+0epYnuwoRcY8QDFKiO2GAJ8ZaZJolUSyFqzQWJYef/uY9jvasOGTU9O6IxycXl9N+sw3aRRoMQSYPYEDDHFBUa6bJOwcwjoHlYvsi4NAENjkLikux69HtaTuxT3zU0ghjy1QoG8g1VoBwiU406XmyInAYiTGV7GEqkJazTgb0xoZIzI4UQrvMDicHTYsOlJ/M3Vivt3q0ekiIOUvz0uCH/ooYcibttwtJVixOP7umZPNd5YPx8x30jC7ofno2ZPteT+t+QR5CQQ3FnIInLEXLOyEqavv4qz4+LptqR8VFVVoaJC6alklIFRKtFAAcAnnRSZFqgWeQBQmhrceFCgQHRxzrQQtBw7rvpGAAAgAElEQVSoxnPr5uGpxXY8t24+Xvid+mfk6qWSz/WxhU7ciGZ8tG2N4n28Jgyu7VqkxrMiScASQ5BtlXYGP+qg+LQLmK5ia5Fnm3z/Kzl6BVScCciNEkE+hxMMvMDicHSorKzEb373Os5a1E065UHKoRKEB8KYh+Jwq+8iKxQHZWPN+M4NTkwbasaOR+/1FgdTE4lXlJ2TQHC9zF4g3wZ45lbgjl8z36yf/P7ohIorAGh0Q7LdOOqhqO8GSlO0HzNXRyNlhCyrL4B34HA1xn55H6YNNeOxhU5MG2rCtx66R7UQdrmVn+tvZQHZApGMnDERAoe8i5XIOliUUhy4RFHVQHFtKrA0J7r1VwJ6flhcf8W50uEidw7HDzNXVGD9DKUxJ8CClB9eGHpBeCB82gUc/GM13nt+C9pcdTARgqWF6dhe7otQWbfrGA68sAXzbq1AeZHUuHFlLsGnndQbDp2dwDol592s2AqF2ah8PHiulxlhlmpYMKTGA1Mm6H9kIgS5CQQuN8Wp7VuwpFAaK3P/7mPYtmUzKisrJY9rdFPVz3WF04FDh3wB2VkWY5E7oSTDQnBZ5EpfZCc42MIKq487gaXZwIocdWPOaPC/kpNuIbCbAfeI8jbuf8W50uEdLA7HD/U6mWnpKrEvkbZt+MUrOyXjrNHhQSx3OhTFQWtDHb7oNCEpTnqhtcYSfC7P96sghhDk2pjhKAC0DExcvC/fyqvpokiL94UxyylzhMa9u2C8a9PRWI8VTn3/Mkop3nJ50Deq/rnKA7JLkkQz2QiRIXu/Csdf3yedwJ0FBCtzTarvWwzxvRfRhtaYkOuvOFc6vMDicHQYHFU6o1NK8c4FDz7ppFiyXhr7EipBuFEa3RRv/Fw6zipx2LFX5tX1rqsd+dNKNEdaZQ4hJJlRYGNhx5RS9I6w9yFYhsco/vpWlVcw/ty6efhkTzVKUtRd0Qmga6EQCIK+K91ZrNiMfNfVjqLprBAe81C8dpbiSBu7Xf65yrV3AFCcNBaSMwaCOJMQYN5ht+QSrJumHPWKiUb9lYDa6DLOBORy/yvOFQ4vsDgcHc70sIBfMcfagYMtwBsuivwl5ZLYl8uJ+XjxxRcnrFkyygeX2ThLHPr7fxZOw/7z7VhbfURSHPzr45s0n4cQglvzfRfgfBuzAOgZH91MZEz4wu+q8Oq/+Tps04aaMfzz+1S38gAWjZMcF5piIGu847Nk/Ubsb2zDWlnBtOqfN2J4jGL7aYpPu3wftBDn84k5D1sPunBWpr3LtgKpfiwGwoF8kxAAlmQTTPfj1B+N+isBtU7VNakEMVFaEHI4RuEaLA5Hh/puaXXVNUTxpwsU16QyG4R3LlDcu7ICpStZQeW0E9yV0RXWM1VVVWHr5u+jtq4OqYXFSHRkY29jO746l2mGbp+WiZL0RBxtHcDhFhcyioqx9uln8cg6/aIv38ZE4SMeX5bfeTeQnMYKrGBHTD/7j83eDpugf1q76xhqX9mKJXdUKu4/UXG7mJQ49ppKV1Zgznd24OArW3H4UB0yioqx5ulnkbW4HP9bT9Hcr+zQla6swMdFq5A4AvxDsfTfotdEUNwuJi2eudvLi35/RKP+SiDTClhjgIHxhuDiLILP5UXveTkco/ACi8PRgFKKUz1U8vVuF4U1FvhiIUFdN1DdQHGqm2LGeAehuZ9iLIx+o0J0z83ODHzreif2NjZhf3Mb9gNYt+sYVjgd2ONqR21br6TjMiOZKHyR5JgIQZ6NoKGXwmYmSI2juNBPMTuNjOuwgrvouc7U4b4FUsH4SqcDPxAJxgUSYoESnc3CQGEZfsDFfuCaWypwsqgc/3e2r0PWOQR0Qv0DEyJn5jqUt81KJUB/6M5plBgTc6hvC6CjGEOk4dfRBiEEhXaC0z0UX3SaQjYe5nAmGz4i5HA0uNDHIkcEDrWylfi7CwniYwiuTQWmJgJvn6cYGW8pjHiAy4Phu0CIo3u+Nq8I21cx+4C0XCfOWfPw1EEXDiIPFT94zVtcAUBJAGHPAllWX35fsCPCMQ81JBgXmJMW+tGQEJKslkmoR/sQ0Deq7P5kWwGHZfKKALkOyx+5Cf6L68mmNJXgKzN5ccX5bMELLA5HA/F4sH2Q4i9NFAsygKlJvk7MnQVMq7Tvou++Tf3hu0jIo3uEDUF3+2U8/MoRpDzfg5IfH8acW3zjQALjgcli008WL8P+uzXITcLmfmDxw1LB+FoVwTjARl8LdYTawSIUVsK4sMVgsehys/dOHjkTSe8rNdR0WHroeU1FC2UO1sXicD5LKAqs2267DXPmzEFZWRmWLl2K48ePAwDefvttzJs3D2VlZZg9ezZefPFF72OWL1+OKVOmoKysDGVlZfjpT38auVfA4YQJcYH1ViNFYhzwuVzpRSDNQrA0m+D9yz47g+aB8P27RS26R+gGNfez0ZG8C5BnI0g0KBoXj5IyLKx4HByj6Alyk9Dlpl7B+DlrHrYddOFEjFQwLjArhSAtDJ0hoeNDCHM+b+4z9jpcvSxyJl7W/dHy7ooUmdbAfn405Q9yOFcTCg3Wjh07kJLCRBDV1dV46KGH8OGHH+KBBx7A3r17MWfOHDQ0NKCkpASVlZVITGT/PPrpT3+K8nJ1M0YOR0As0C4pLsaGTU8qjB6jge5hisvjo6SuIYoGN7C6SH3UsjgLONkBvHeJonIKCWsHa8OmJ7F69Wrcv/sYlhc68K6rHe81tmHN08/iw3aKRLPSoDMQTZPNTJAaz7RJQqekbdBnOBqo0L1x3P+qdGUFHIvK8csairunEa9mTcyS7PC8b2LvqEI78GE7K0z1fLYE/ZU8CifRrAxdjjR5AeipTMTnlcXhcCKL4p/aQnEFAN3d3ZJRRFcX247q6emBw+FAfHx8hI7J+SwgCLTtnc341vUst2/16tURz+0zwqedvi5HXTe7UM3QGLPFmghmJjPfKADoHSWaQcITpbKyEi+9qozumbFsFT7uBOakQREvUmxQfyUgjAkdMs1SoDosSqUeYnXdzN9IzaF9RjJBdkJ4Cpe0eCb0BtiWZ98oC5nWo2uYWVTILQSiYRsvNZ4gyWzsvjlXgP6Kw/msorpF+OCDD+Ldd98FwEaDhBC8+uqrqKyshM1mQ2dnJ6qqqhAX5zOCefzxx7Fp0ybMmjULTz31FKZOnar5Q91uN3p6erxfx8fH82LtKkAs0PYXVzLZnOzwFQb13RRT7MpRkZh8G8EHLRTuEYoYsM7NtWnhOdvM5RVYv13aLf6wnWJwTDkedMQH3nEpsAEfdQBmE0FqPEXrINsgbA1wk/DyADAo8uJ09VIU2llBKmdpmLpXgC/Dr2WAdXMImL6q5e++eKF0ZzGWrN/otdsQnOfl3Z9o8ZNyJhJ81OG/iC+KkvNyOFcjqgXWSy+9BAB48cUX8dhjj+GNN97Ali1bUFVVhZtvvhmHDx/G3XffjY8++gjp6el4+eWXUVBQAEopfv7zn+Ouu+7Cp59+qvlDly1bJvn60UcfxWOPPRbClxWddHZ2TvYRIs6bb76Jn/zoRzh9+jSGhofxncUzlLl9h2vR0dExySf10TpIcLadtQiGxoCG3jiszBrFwIB2NEpGDADE40zHMArMg/j04hjyEB6n7yNNsRgQ6bwoBQ5cMmOancJGRzEg2pLLTRhDR0dg57CNEAwMsNefZo7F5T5gYGAI59o96LCN+nm0j5MdJgwMsF8xYxRodMfhpvQxDAxIz5Of4IF9eBTh/CNgHfW9Z5kWM47/qRoXfrIWS53peHghs7vY8ei9KN/yCmYuW4WzXbHIjCfAyAgGRDl5yaMj6BgvbCbz73OKx/fe6t5vzHdeTnBcjb+3r1aC/azT0tT/Na37N/QrX/kKHnnkERw9ehTNzc24+eabAQALFixAfn4+jh8/jltvvRUFBQUA2AXzG9/4Br797W+jvb0dDoeKgQyAffv2oayszPv11dTB0vogPotUVVXhoYceYp5NC5z4z0Nn8G5Dm9cQ05vbV1oSVe/LsSYPrFZ2UTrbSeEBxTXpcbDqBPtarUCS2YPLI2bMTLKgy2RFWlpMyM82MErRTj2winRFp7sp2oYo7iwksFqls6MFBSakBShyTqEUSS0ejHiAbJsHH3cCVqsZ/SYE9Jq6On3vY1MfxQilmJZqhtUqtUC/fZoJaQGOMQOlaNAD1wg7y5QkD469/pTC/HTdrmM4/PLTuO72L+HCoAczkgCr6I22m4EZOdKW0GT9ub3WSrG/Sz8L0USAOfk23c4rxxjR9PuJE15C+VlLCqyuri709/cjNzcXALBr1y44HA4UFRXh4sWLqKmpQWlpKU6fPo0zZ86guLgYo6OjaG9vR1ZWFgBg586dyMrK0iyuAMButyMpKSlkL4ITnchHgvmJFvzT2yewrvooVhSlY29jO/a72lD1019P9lG9UErxyqtV+POv2OjInDMTyXdvQMq81X4fW2Bj3llwsPHYqIeqjsMmQl230sX77y0UuQnK8ZUtNjiDSRMhyE0gcLkpMixMTzY8RtEDgt5hangjsVG0rdfgZhYJ8ny5bCtUBe+hhm0SsvM47QSHLtVj+Q1Ohd3FoUN16Bmm6BxS6q+iKXw4w8JMWft1GorZVsKLKw5nEpGI3Lu7u1FeXo5rr70W1113HZ555hm89dZbyMrKwq9//WusWbMG1113HSoqKvDMM8+gsLAQQ0NDuPPOO72P+cUvfoE33nhjsl4PJ4qQezZ9YUY27inNxcHmTjx92IX+tHxUVVVFLLfPCM9ur8KL3/Ll5s3zNKP7v+5DzZ5qyf3iVRo5+XaC5n6gdu9u/GrtfCTZEjB/bllIRfy1XdLq6lI/xbleYFEmUWzFFacov2cUwftJ2MATBO4NbmOPbxukcItGa429FPk2KExEb8qKjBWf2Duq0A6QrJma5qfC5qO8YC2KIj8pwf1cj2g6L4dzNSLpYDmdThw6dEj1jmvXrsXatWsV37fZbDhy5Eh4Tse5omGeTc2SkeDl/mFkTr8Wj/z2KB6ZZfK6bEcLP35KmZu3btcxHHhhi1cADQD3TDHheJs0ILjABowdrcYbzz6Apc50LF/gxL5Gtim5c+fOCQv5RzwUZ3qkBdb7LRTJccCsVOX9r0sL/r1lm4QU6eOT+9ZBZg/g6qW41sDzunp9/+2hFK4+VgSKMZtCG4ujR7rFl+GXEEtgK9+A/f99nzdeSLC7cOTZsfuuRMRnz0TjNzZJPvNo85Ny2oFandjLaDsvh3O1wZ3cOWFjw6Yn8TdXK+7fzRy87999DPtdzMF7lAJvuKi3gxANjHgozp+pw/JCh2J01Nrgy81bkMF8nD5fQBAn+huUbQVif78VSwtZgfa1eUV4ZdVcLHWmY9uWzRM+35keFsUj0D1M8XEHcGMmUVgzOOLZplmwCH5XcTEEKXEY3yRkxqFGEN+vZYAtC8g7QtOTImchEGMiSBVJv4pXVsD+L6/inJXZXXw6ylp2s2LdeOIGJxaiGTsevdfbubTFTr7/lRy9kSUB97/icCYbXmBxwkZlZSV27tyJ/rR8/PiICy2J+RIH70Y3xZG26CmwarvgNzcvwwLcls8ubMlxRGIvEGMioJfrsdzpUGxK1tTWTvh8dbLx4MEWirgY9TDiiWa62c2ssALYa24TeWH1jfj/zMT+Vy4386Ga7MgZcYHktBO4r63Aff97BE8ccMNiT/F2Lr82rwjby1nG44EXtnjvH21kJ6iPqoXbLLHRd2YO52qCF1icsFJZWYkjx46job0f9794RBGP8pcmip4wmXIGyol2iiXrpbl560S5eTEEqJxiglmkI7opiyBdpO+x5c3EXpe0QNvb2I7SkpIJn69BVLQMjVEcbQPmpyv9uUxk4gUW4DMcTbcATe9V47l18/HUYjsWzp+rqyvrHqZ4/x3f/ff/P/OR+km1RPAfS4znI4YKsQ5L6KYJo8w2l37nMhr1TCZCUGBT/5z5eJDDmXx4gcUJO5RSvOGiGFOpoz78czXKyubClmANuSA8ENwjFGd7fbl5H5vzsPWgC2etvty8FbkEOTK38RgTwR0Fvr9Gsx/ciP2NvgJNGItu2PTkhM4nbLYJfNoFDHvUw5GnJxnPHtRDGBOOHN2Fvmfu8wr/03uadB34X/hdFXY8Kl0UaPuZdFFgWlLkN9zEej+7mcAR7xtl+utcRoODuxpaxqfRel4O52qCF1icsONyA+dVAnZr9lRjx6P3Ir+/adKjc060U6/9QenKCiRvPozpL/Xgq9uPonRFOXISgMVZ6hetaUkEs8YDgBfdXgE88ipq4pi2pzslNJuSjbLtvdpO5oqepFJIzQ1B9wrwdbBcv9uCJRJd2TxdXdkvf7hZd9wGRH48CEg7WADTKAnvq9C5XKvSuUyIVT42WlDbJOT6Kw4nOuAFFifs1HerjwDfe36L5EIcSkF4IHgoxaFW3xndIxSNbmBmiu/iNTtV3/JA0GXZzQTJN5Rj+o8O44kDbmzfc2zCxVVVVRUqls7FU4vteG7dfJz8SxXO9AKlKcrz2GJDN3rLsrKLdff5eqwIQFfmb1EghgDFER4PAsxyQvyOOe0ELeOaMteMVcAjr+JkrDTjsXRFOZz24O0uwk2ejY1bxcxNJ7By/RWHM+nwAosTduQF1vAYxYiHqupeQiUID4T/eqkKP7pnvreA+fOb1TARYJbIQqBYpZgRkxLPRk4AkGv14MJ4Z0TYvgsWISA7p4+N26YNNaH68TUYO1qNEpUiZU4aUXhNBUusiSAlfnx8ZlBX1jdC4fAzbpuWRCZFgG02ESSLNgkFXdVv6ikOtQJfWFWJb+44iicOuLF+vHMpvl80EmsiyBvvNBKwQv9uJ/+1zuFEA/xvIiesdAxStA36vh7xUDxfR/FiPUW6U2n2GCpBuFGqqqrw/z50j1cvNG2oCSe3rUFubbW3C5BugSG/LmGklmf14OK4k3vrgJ8H+YG54YvGbavYuM3yzjakqET3zE0PbeGSYSFsfNYoFf5r6cpcbuguCgCTMx4UEG8SJsexTcneEWDtNIIFKno2IDo3CMU47YAlBrh/uilixq0cDsc//G8jJ6zIu1fvNlN0DAFN/UD2GumFOFSC8ED47ve+L9ULrZqHxYXpGNi9zXufYoNRLvnjupdcK9NzNfcDLRPsYDE3fOW4bbS5TnHfPBtBZoi9mjIsTJM247FXcZj4xmeP/NfrqqNPl9u3KHDKwhYFauN847bJGg8KyLVU900l+KcSohnXY41ho9JoZlYqwVdLTJgegcghDodjHP9x7BzOBDjV4/vvRjfF+y3ArXkE7YMUn5jKcfdTO3D0f7fi0KE6zCwuQVXVryManVNfX4dHF0gz6VY6HTh8yFfA+BsPChSMu59nWijMJpZLWJQIjHlo0GM75obfJHHD3+Nqh2N83Cbm+hB3rwCh40NRsqISp6dV4DtlBLEmAhNhVhHyTUBXLysoS1dWYHBOOd5opHh4jm8kOCVxcvVB4kxCAMhO0D+LMzF69VcC/l4Dh8OZHHiBxQkbw2MUDeMX3OExil0uigIbcGMmMDRGUNtNcXZ6OdZvZxEyN2cTrMyLXFO1d5iOr+dLCxixXiiQwORMKxBnAgYICzW+0EfhoQTtQ+y2YNiw6UmsXr1aEulyoLENa55+VnK/JDMwJ3Qh8F6Ejk+GhZUl7UOso+OhbANvxng3ilKKdy5QXBKNRF1uimyr1PByMseDQODbgFo2CBwOh+MPPiLkhI0zPfB6X/2lmaJ3GFjlZLEu1liC2/IIPu6EN1/vfF9kz3eolWLxw9Ix5VqZXmhmsvEOhokQ5I7rsLKs8OqvJqLDumtVBe57+jVvpMtB5KH0iR0Kw9ZFWaETt4sRChLBTLVVVkABTFe34yzFwRbpONTlVhYoM5JCfsSAyE5Qbt3pwf2kOBxOsPACixM2BP3V2R6Kw63A5/IIHCKx+Jw0oMgO/P482yps6qPwRCibcMRDcbTNpxc6Z83DtoMuHDXlSeJ8SgyOBwWEbleGhXWuxjx0QpuEF/qAkhUVWL/9KO77Qy9GNh7G0jukI1RrDDA/DONBgGURJscB1lgCu1m6FdnQC/SPUrxUT1Eji/HpHqboGpbmIabFIyQGqBMh1uQrgv1hiWH5khwOhxMMvMDihAVKKU6Nd6b+1EThtAMLM6T3IYTgzkKC7mHgwGXmTH55glt3Rvm4A+gfZf9durICa//3COgzPVj+C1+cj9kETA2w4yJsEmZYfSO11kH9x+ghDk2u6aJIiVNe9Bdmhjc0WdigzLRAMgJs7qd4odajaiIrRNCIO1jR0g0yOvYrjGL/Kw6HE/3wAosTFpr7AfcI0DVEcXkAWJChfrFKtxDMdQAfdbCL9AWVi3U4+LiTyr4GQIBrRTqmaUlEkjtohHxvB4v9f+sg0DIQ/GsSnMYppfj/27vv8Dars3/g3yMvWZKn7Di2bMmZdrZHSFjZjEBTYsskIQmFtEDL2/L2fQuUQDP6KxmQAqWlLaRpeSmQlMbEcgIUSoEMwkqcOCEM2yHDkmMn8R7ytnV+fzySLNmSLcuPhu37c125gq11Hh2Ebs65z30X1wOpkbB7H4MkwNwxng0CLNeiVjAYjMJYig7m4y+rM/HobKF2mG0bHEAIDGOlgMxmP85f+uO5WnbBn+tfEUL8HwVYxCMs24MlDULz4Yn9rARpFAy17UKRyjKj8/uJpa2rJ/ne4nQNx+Rw+4DA1fIMtuRBDJFBHLJABnmgEGDVtsOtrc9uE7cGnJeaAWNX3+rtGTHMbsyeEGtewVIrgLZu4Iv38+16DU5oL0fuYyvsgixH+VcaPwlYkhTCv5MD8ZeAkBAyPFGARTziuwbh77MNHMkK9NvY15K3dKnZcc9Csb20R4edq3sqt3/yng5XWoE0mx5+DO63nEmQCdcQKwWqW4Um1zVubBNWtACdJuGfLzQJOUGJNicaAxhwvZP+iGKKNW9JJsqFwOT4K1v71A5LjQnDW5vvwVM3KPCX1Rmo+TzfbqUoKkQo7OkPQgLYgLWtQgKEhHhCCHEXBVhEdE0dHJdbONq7OUqNwkk8WwkyhuiQnp8jggFFoLA9WNcu9AL0FJ1Oh4d7VW7/aNNKKL7Kt5YcAIAkBYM8yL2AICFUiIpipD35V+7kYdnmX+mNQnNnic324PQo5pWgxbJFGCRhUMmAprKzdi2O/n2+EsXVTZgdK8P6ORpM6qgAdq5Cx8meFS1/q4Y+0HjUCmb3XhNCyGBRgEVE900dB4dQpsHE7St3B0mAnHEMSxN7/tVjjCFJIaxgAfDoNuHWJx1Xbg96d7vdF+qsaPe/XBNCzStYoQzV5u1Bd04SWvKvuk3C1mnvoCDDQycHe7OcIASEbULJWPsWRy8cP495SUr8IzvT7j0t/Ps263P4S4K7xUABVjLVvyKEDBEFWER0X9Za8q84xkhh1zNvaaIESinD5Ehmt7KVKGcobxGCEU9uExaX9G0wvVijRPOlnsrtQRJg+hCKdsaZK7nHSoUAs7Z98LWwOOcwmFewKlqALm6f06QIEoIdb7HkYWkUDJ23b7CrHVZcY8TC5Jg+72lVac976m8By0Dvnb8FhISQ4YcCLCKqV97Iw5NZQn5T0cOzEf1tzzZRaiRDZqxtsMWsRR8T5UK+0dVWz65gqcan9GkwbVu5HQCmRrJ+c8YGImHCNqjdScJBrmBdbRUSygEhYTxYAsTb5ARNifRuCQFLodEkBcAysjF7Y661+GlgsBSH9c7f04hgOGxM7UvyIGa9pt5CAuzfa0IIcQcFWEQ0Op0OP1rTk980FxUofnoVig7mQxEE3KGx/5KNljJrknaCTAhMypqByy0cXSbxV7E457h+gMrtAJAuwtZbolxosxMaIKxe1bQN7iSh7SlHR/lXUwdZAHWoLMGi1Jwg3pmWhfv3nMQTnxpxx5Ov4ajB+Xvqr6fxnG0TJskp/4oQMnQUYA0jOp0Os2alQRYaisz0NOh0Ol8Pyc7WLfb5Tf/IysCN6hh8+vJWZGkkDssJzItniAwWKmzHhwKXjBxdfGjFOZ251AyMWyBUbr8QqsK2Y3p8GWBfuT06RJz+c0nmIpWxUqC6TThJ6GoRVZ1Oh7WLMsynHDNQeiQfaptgQB7o/ZIHsTYV+DUKYVXNwrYa/lPH9Chg9u+pv263OQuwqP4VIUQMFGANEzqdDjk5OYhsKMej12igqCtHTk6OXwVZxcV985sWaZSo1pdgopOaUkEShlvMCe+JCuDCkXz8dU0mJihlogeRxeZ2LlMWZ+OaF07A9KdG3PPaCbu+fmlKcbbeLOUUYkJ7gkV908ArWJZ5VrVYTjlWoOulVUCh/Vart1dYYm3KGmgUDPUdQjsciymLhXY+sf/XiNTnCuzeU39tmOwsD8vfTjwSQoYnCrCGiW1bnsR8TQx2m1eHdi/PwDxNDLZv3eLroVklOMlvGjcxtd/HTY1iiAoBJIX70fynVZjQbgkiK0QNIi0BFuccn1/lmBQunPSzYLCvhTUUiiDh5F2slKHavD2odyG3zDLPvU85Fu/pOZE3Ncr7AYAiiCE0QPhnS2Bi6HU9bd0cV1rsA5TwIGEr2B9Fhgh9Fm3FSoXtakIIGSoKsPyYTqdDZnoa5LJQnDnzFRb0Wh2KkwXjzJkzkMt8t2VoO8b6+nq7/KY15lycJzZuGvB5EuUMF9/YihvVPcHF7uXpogWRVa0cNe3CP19qBirbgGt7tZgZH84QLmJdqVipsEXYxYH6DiGXig+Qh1VcUtJnnhdrlKg2n8gLDQDG+WgLyxKMWhLEbet0AcLhBA7/7D/ojG0wODWS4YFUCQIG2R6JEEIcoQDLT1m2ihR1FXhktgayIAkOlVZbv6Df/e4K9hVV4FpVFB6ZLf5qjztjnBokFLIqaFNg+zE9Loaq8MDv9+GeVdoBnytRDtQazmKRxj64WKhWoqi4eMhjLWmwb5osD+yba5Mu0oQVEVkAACAASURBVOqVRYzUpidhq9BceqDcstSUFBzp55SjL7YHLWJtTt1pFPYrWCbOcaaWQxEIuyKy/r7dplYIK5c3qRhWTpB4tGk2IWR0oQDLTwlbRbHYvTwdP85IxjNLpuGTshqszj+JXYWl+OVH3+DGJCX2ZGWKvtrj7hj3LBeS2rukEbj9rSbcv+ckVt2Z7dJzJcoZYjQpfY77HzbUYEpq/1uMrrjYBOtzFtcDKZH2p/JCA4RGymKKlQJhQUDgqXy892OhdMWCOen9BsEbNm3Gx/2ccvTF9qCFbVkDtYKhqk3oH9nezfHP8xzf1AELE+xz2Pw9YXx8GMPaiRLcOJb+U0gIERf9V8VPCVtF0dYvq9snjcWdUxJwrLwO247p0dxpwqJexR3FWu1xd4yWpHZ+pQRTo4T7uNoweWwoMP+BjThqqMZqc3Cx9kAhjuqrsWHT5iGN08Q5yszbWVdbgbqOvk2T05QMgSJvDcVKGYoP7UfXS6sws6sC6+doMKap/8MJi75nfyLvGFTI2fEmpizKgjQAGO/DgKX3SUIA+LIWeLmEw9AMrJnAkGlT4iIsCFD6af6VRbSUOT2AQQghQ0EBlp8Stopq7VZzrjR3gCfOgOTPjZBppvdJKBdrtWcoYzyor0GoKgXSAIawIEAlH+BJzAIkDDd/X4slW3LxBVR4+rgeNeGJ0Ol0yM52bRXMmcstQIe5aXJRPUdIADDOJk8ogAHXeqBpcmwo8MnftmKeOgZ7slw7nKA39pzIm7y7EQlPFWD6EuH6UyKYT/ODxticJIwIFsprfFAulKC4P6VvoOLv24OEEOJJgb4eAHFsw6bNyMnJwdoDhVioVuKQvgafGKqR+viLSEpg+ODWDTi6cxXW7C/EIo0Shw01OKqvhu75XT4d46eGasz/zUsAhCbPgyl5kCgHrr01G4fHZmGJiuHRmRIsUQ3t/wF0Oh1+9esncfG7EsRoUmC89VdIWZBtF6jMiPZM02RFEEO1vgT3zdH0WWl87oTjlUZLKQduPnWYEdNzW6qXi4v2Fh4snIw0dgo/pysZrrRyfF/NEOqgxpm/bw8SQogn0QqWn9JqtcjLy4Pe3I6kgKkQ8fNc5GizMXcMMOb6bET/z15cCFXht8f1aIkWZ7VnsGP8/av7cEEqjPHLQBWkD+Viwe1CDSRXtwctEuXCNl28DLjUzIdcbNSShB/XZKkpVY7GP66C4qv91vswADd4YPXKIn5c39IV/a00Wk7mVbcBzV09q0AM3i8u6ojtqtT8eIaV4yUOgyvA/08QEkKIJ9EKlh+bd1s2fqhZjoMVJnx6Bbg3tSdP6Htq4JW2bHxvqRazYxl+OVMCeZD3v9DC5mbj/j1Z6DJxPPcVR2aMkDweLAHGhw/uuZLM24kqOfBdg1BaYShsa4cxxvBAugar9xfi4j+3Ad8XTjZOjmB2tbDEtvp/N+LpB1dg9f5CLNb0rEQ6Wmls6uCoNZeS0BuFoMrynsSGwmElfG9TK4Bv6ga+nyIIiPHz/CtCCPEkWsHyY2cbOMqbOT65IqwWxMt6vrDUCoY0JfBRBUdzJ0dZs/fHV9HMYTCvuJxtEJoTWwp1TggffNJ4WLBQ+HGMlKGuXWiQPJSehAPVlAKAG8d6Ngi4I0uLFTtycQxC0vrFUBV+9Pt9DlcabQuR6o0cCTJYywb4Sz8/V/OqKP+KEDLaUYDlx842AIcuc8SFAjeOFX7HAASZZ+1mFQMD8GF5zyk5b/qisuc1v6zhUMl6Vi3czRdKlAvFOTl6miS7a6CaUhoFQ5KHA4FYKcPUJVqM2V6A6W804v49J5F0YxYaO/rOl2V70JJ/5Y8FO+NCAWnAwPej/CtCyGhHAZafauvi+K7RhItNwqpQgHkV5vo4hkUJwj/LAhmui2P4ug646EKfOzHYVm7/+W0ZKDqYD2Mnx3eNwCzz6pWEAZMi3Hv+RLlNcc62oTV9HqimlCdzrywsPfxipfbNnnu3zalp49ZWPnUdQFMnoLEJqvylnx9jrgWl/rLiRgghvkIBlp861wicawBMHEgxBysJMobFCQxzY5m16KNGIbRiOVXN0T2E7TRX9K7cPr6tHLmPrcCHb+dDwoDp5tpXyQrmdr5QopxBGiicVqtu46gcQh6WVqvFvc8JNaW2H9PjVIAKK5/ZhymLshAXCkz2wqm88CAgWAIkKRiutADt3cL12DZ+vtjI8bdiE5rMp/P05qKoakv+lRQ+ya9zRj1A6Q15IDya10YIIcMBBVh+xHZ1KGteOo6/r8MYqdCUNlgC5IwT6iAFSBhuTRSmLkEmrBjpjUK9J09yVrm9eM82pETAeppsKA2T42VCXapY6dBXsOraOZIXZGP130/A9KdGLHnpBKYsEk442hbE9CTGGGKkDGq5sO15yZwrZ9kOLKzm+NWuPLywSqj0/tc1mfjyQx3GhgJS8/vpb/lMmgG2K/1lO5MQQnyJAiw/0Xt1KKG5HPpnVyH623wAwG1JEruq2JMiGCZHmMsahAplDcqaPbuC5axye+flEmSYAxZpAKxV3N0hlGlgNgGW+9dk2Ya72CQEN5ZtS4a+ldw9KVYqtJmRBfYEVlVtwNt6E3a8kod//nIFJrT3lJLQP7sKYV/lWx/vbwFLggzob4HSX7YzCSHElyjA8hOOVoduUMegNm87pkUxpDtYcVmayBDAhLylsmagzOjgiUXkqHL7IX0NguJTrC1cpkcNveWMkIfFUNsGVLa6v/VpOeGoN3IoQ4Aw8zZbkoIhzAOFRZ2JDRWCUbXcvkHyyWouVHrXxGDP8gz7eddtt97P3xLGAyUMKrnz98/fAkJCCPEFCrD8hKPVocUaJeoNJfi+2vEXVrSU4doxDIlyhvoOoKje5NExCknjVVhrThpfY04an3nvRuu4HQWCg2U5SWiCcIqwpt2957HkOfU+kTdV5KbOA7GcrFQrGC41A19/pMNf1whbglfPfYWFDkpJNJUJpSSUIUJFeH+jdrJKJQ/sOaRACCGjGQVYfsLZ6lDihFRrLo4j8+MZJpu3vorqgfp2z20TarVavPDaPlwwV5c/HSBUbr91mZDXFBeKflc2XJUo7zl9V9UmlGsYLGMnR0070NIlVIS35A0xAFOjvBuwxNocSOg+mY+89SutW4LyIAkOlVb36edoKSXhr6tBzvLCNIrBtUcihJCRiiq5+4neff0Omvv6PffqX/p9XEgAQ7qSISyI4/SHOlz/v9uhP1+C1JQUbNi0GVqtVtRxqucJlduNnRy//5pjYXxPA+KhJLfbigxhGBPKIAsUgqPKNo5pGNxzW07iWf62rGCp5AzhXtweBICoECFnaawMCHx3G65TC1uCjDEkhknxk3e/xOr8k1icHGOd95XP7ATgvwFWkkIIVnuH8/7QzocQQvzBqFnBsj2hl5meBp1O5+sh2bH0HqwMS8RTx/Q4BhVydryJB9cM3FswScEQ8XU+Sp5ehdimcjwyWwNFXQVycnJEvU7OOc42CF+pBVUcEgZkmpsRBzBgZrR4wcAYqSXRnbu1gmUwJ/yXGjkig2Ft5jyUBHx3SRiDUir8za+cxUJNz5bg7ZPG4s4pCThWUYcdx/X4Aipcu/lN62lHf00YDwlgGCuz/x0DMM5PA0JCCPG2URFg9T6h54ngQwxarRYP/fMEEv/eiOTfFmBZltalelKJcgZj/jbcqBb67v04Ixm7l6djniYG27duEW18BqPQgLjTxFFQBaQre0ozTI5gotZqig0Vttaq3TxJWGqTf5Vsl3/lmwAg1pyHpUicjMN6++ryV5o7EDdxBu54pwndGwtw/W1CcBUdAq+vtg2G2mabUBoA3DVBgjFU/4oQQgCMkgCr9wk9TwQfYqhoFlre6I1CwJLiYjX0sTKgufwsFmnsk6UXqpUoKi4WbXyWSuOna4S+g9eO6fkyTRdpe9AiRiokh1sCrMGcJGzt4qhsFf6+2tqTf5UgY4gM8U0AYCkMm3HvRhw19FSXX2NTXb60iSPa5rSjv9W/6s0yvngZ8JMpEqT4KHglhBB/NGIDLNstwTNnvupzQk/s4EMMZ2o5zjUKeS2TI1zv5xcoYYgfl4JDvVZGDhtqMCU1VbTxFTdwcM7xRSXHlEggyhyshAUBE91sjeNMrHmLsPtkPnauzkSYXOby1q7BKLyHlpIIlm02X2wPWlhWsK6/LRuS/9qLb4KEgwIXQ3uqy/tr/0Fn1AqhYOt9KRLrvwuEEEIEIzLA6r0lKHNwUkvs4GOoTJzjq1ohx2lsqFBIdDBfWvf8YiM+sVkZWXugEEf11diwabMo47vawlHXDlS0ALXtwOzYnrFdE8sgEfnkWKwUqD+WD+xchckdFXj0Gte3dvU29a/Cg4DIYOH3vtoeBHpORQZKGFQ3ZiP+qQI88akR9+85iSmLsnpOO9qsWvlb/aveFEEM39dIhlz3jBBCRqIReYrQdkvQ0Umtw4YaHNVXQ/f8Ll8P1ep8I1DbzlHSANwYx5AaMbgvrVV35iD3wl4c/+d2nCgowbQpqdA9vwvZ2QMnybui2JzcXlTPIQvsWWkJCQDmxIr/BSsPYjj1962YZ3Pi7oF0DdYeKMT2rVv6PR1pqeBuWRFijGFsqFA3zFeUIUJLIxMXxnSqRgj0LauqpebTjpagyjYxnxBCyPAzIlewehfttJ7UKq/DtmN6nJSo8IdX94kWfIjhTC3Ht/VApwmYpXR9e9AiUQ7MvEmLjg0FeLrAiBOFp0S9vuJ6YXuwqF5oPm1ZscqMYf3W6RqKKr39iTtXtnZz9+XhN1kZQhHPJ2Yj+Mv9AIApXq591VuAhCHOvIqlVjA0dwkrgQBwtoHjLQNHgqwnqPL37UFCCCH96xNg3XLLLZg5cybS0tIwb948nDp1CgDw7rvvIiMjA2lpaZg+fTpeffVV62MqKyuxdOlSTJo0CdOnT8fHH3/svStwwFHRzsvNHZBqZuAH7zeh5YkCnFIvR3u3Z3v3uaqjm6O4nuPLGo5xYUCSnCFeNvDjbEWGMKSa86AKq8S9roYOjsstQGWbEBRY+vgFMOC6MZ4LBJLGp/Q5cdff1q5Op8OqFXdifJtQxPM6VoGTW1ei6GC+X5QPsJy6SzKv/pUagSOXOd44z6FRAD+YZLs96PvxEkIIcV+fLcLc3FxERgq9RPLz87Fu3TqcPn0ad999Nw4fPoyZM2eitLQUqamp0Gq1CAsLw+OPP45rr70W//73v1FQUIDs7GxcvHgRQUFBXr8gwHnRzkde3ImfZUjw6dVu5Os5lqlNyBkf4JMx2iqqF3ru6Y1AtoYhJdK9atjToiUID+rGuUaOli7uUokHV5y1bg8CIRJY+w7OUnq2p999j2zEE/etwOr9hVisUQ64tStsDQulKixbimv2F+Kz/9sK1SM5Hhunq9RyhmPgkAYwjA3leP8SR6cJWBjPMH8s7OY82U/rXxFCCHFNnxUsS3AFAA0NDXbbM/X19QCAxsZGKJVKhISEABCCsgcffBAAcM011yAhIQFHjhzx+OCdsRTtbIlOxNPm4o23bXsTT/9Ei3HhEvz3VAmMncDvvzbhcovvV7HO1AqrV8ESIDVy8NuDFoly4c+lZvfayzii0+lw1wJhy+2LhzIRV5yPAAkDA3BDnGdXWVbcmYN5/y8XX0CFp4/rcVmRCJ1O53TrU9gatt9SXKRRoqq0xFpt3pdsq5xPDBdystZMYFgQbx9QRwbDZ+UkCCGEiMNhkvs999yDQ4cOARC2Bhlj2Lt3L7RaLeRyOerq6qDT6RAcHIyamhp0dnZi7Nix1scnJyfDYDA4fVGj0YjGxkbrzyEhIdZgTSxarRZJ87Jww1vdyIgGXrk10Hra6Z7JEuzXm3DkCvBBuQn3TPLdKlZjB8e5BhO+rBXKCIQHM7dXLxLlDCo5w6EKjqutJmjChnZdltOY8zQxWDtHg0P6Cnzy7CoUjXkTOVotlB5OGo+RApm3ZONoQhZWTWC4Po4he5rza0qZnILDhnI8kK4BY8zazzF5on+cFlUEMUSHCNusCxMY5scDQQ4CP9oeJISQ4c9hgPXaa68BAF599VWsX78eb731FrZu3QqdTof58+ejoKAAd9xxB7766iu3trIWLFhg9/Njjz2G9evXuzH8/m36QgqpJBBPpLYgqM2E2rae2x4cF4BT1aF4qrATy5QNor+2I3V1dXY/v/3229jy9LMwXDwH05jJUKxbj+SU76Ghvtut5w81AVEBIejiwTh+qRkTA7uGNN7f/PrXDrfcjv5tCzavWYLaWs+v/slNwQiSBKOiqQuG4G6UVXZA7uTs69qfPoxHf3Iv1uwvxCKNEof0NfjEUI3ndv0OtbW1Hh+rRe95thXFA1He2rNw7GiGokxdqK01eWBkREz9zTMZOWieRw935zo6Otrh7/st03DvvffiwQcfxMmTJ1FRUYH58+cDELYBExMTcerUKdx8880IDAzElStXrKtYpaWlUKvVTp/3yJEjSEtLs/7siRUsE+f4vKYL8+IZslIj+gSC2ZEcB6524x/nGTpCIjFW7p0DlZaJ0Ol0WLduHeZpYrBmjgaH9RU4umUt1k3Zh+hp7ucLzRrThVwDx/l2OaKjh7aCde78OTwyW9Nny62g4CymqbxTtVOjNCFW2o367kCEhgbDGCRHkpMTgZnL7sbK38rw/s6t+OxYCeLHpWD1Mzvx8/u0Xq/V5OwDN83Ecb6j/+BppooKdw4XzuaZjCw0z6OHmHNtF1XU19ejoqLC+vP+/fuhVCqRnJyMy5cvo6ioCABw7tw5nD9/HikpKQCAFStWYOfOnQCAgoIClJeX91mlsqVQKBAeHm79I3ZwBQCfXuFo7AR+NlXicJXt7f35KPj5bEgeCsd112R4vS+hJSF7j7l34J6sDMzTxOClZ7YO6XlnRDEES4DvGoa+uiScxrQ/xXdIX4Pxk7y35RYrFf5YcsosRUQd0RuBKYuzEbm1AJNfb8T9e07iljuy/aoQ5kDNmyODQcEVIYSMAHYrWA0NDVixYgVaW1shkUgQGxuLd955B3Fxcdi1axdWrlwJiUQCk8mEP/3pT9ZVqh07duAHP/gBJk2ahODgYOzevdtnJwgt9pwzQRYI3JLY98vKkls0XxOLu+ZqcNhQjpycHOTl5fVbwFJMxSUlfVaHFqqVeO7E0Nr3qMMkiJGaUN7C0drFrc2Y3WE5jdl7y+2NN/8ypDEORowUiAtl+KaOo5tzaxHR3rpNHJeaObpMHJeagZtUln5+XhuqS5RSBkUQYOx0fLu/9x8khBDiGrsAS6PR4Pjx4w7vuHr1aqxevdrhbXFxcfjPf/4j/uiG4D/lHDfEMYerF70rvbtaIVxM4yf2TcgWo31PorxnxaeqTegX566blmVj5W/fxNsvbsXnx0owdlwK/vtPO3HXnd55jwChh59aAXRx4HILECQB2rs5QgLs57WiRSjSeqkZ6OY9ZQ6S/TBgSZIzFNU7XomjBHdCCBkZRmSrHEOTCRebgPWzHH9ZeWr1aDBueXAjfvcz12s8uUopZUiQCV/glS3cWtzSHQYjkDw/C22Ry/E9NUNmDMP31d4t/h8rBcbKhMDKYBQCyDJj3+bSPf0HAWkAMCYUCGRAop+tYAHCqlpRvePb/L3/ICGEENeMyFY5e85zSBiwarzjy3NU6d2bzZ/r2zkirs2G5tG9KGAq7DiuR11E/zWeBiM1kqHDBBQNMQ9LbxS25DiACWFC3abUyAEfJqqoECBEwpAoBww2TZxtNXRwnK7puU2tEFr5qOTMYRkEX3MW9EZQ/hUhhIwYI3IFa3+pCTOinBdrdFbp3VvNn49XcdS3c5RNzsZNL2pxWxLD/0x3nIzvjunmQxBfD7Eygd7IUWrkiDAXvtQoGORB3g0AGGNQSgG1HDheJQTDemPPGP6yJw9bt2xBZWkJYjSTcWXJBixZJmxh+lv+lcVYGRAsAXofJvTH7UxCCCHuGXErWM2dHIXVwLJ+trJsK70/d0KPUwEqzPjVm15p/tzRzVFYzXGskiMoAMiIAdJj3GuN48z0KIZABpQMYQWrtYujqhXQN/XkM031TmWGPoQ8LIbWbiGvrMKczP7bV/bhwbvvxDhz78EJ7RUwvbQKpsJ8AP6bzyRhDInyvmPT+Ol4CSGEDN6IW8F6W29CF7dvnOuIVqu1JrSv/qgLH1V4p2VOSYOwpXWyGsiMBaQBDGlKcb9Y40IZYqRAWTNHWxeH1I2ThAajEGRdaQXmxAqtcaa42cJnqGJDhdwrCYBj7+fjSu42bDWUgIFhnloodWE5rLB6fyGKdm/Dgtu11qbK/kgTBlxo6vl5cgTDVC9vvxJCCPGcEbeC9c8LHAkyICXS9UtbEM9Q1QZcbvZ89eziemGFrdMEzI0VGjtHiNwwOTZUKG9QbT5J6A69kcPQLPyzJkzIG1J4eXvQIlbKEBzAEPlNPgq3rsSE9nI8do0GXR1tWKix7z24WKNEdWkJEmT+mX9loTavYAUw4NZEhtUTmFuBMCGEEP80ogIsE+c4cpnjFtXgvqhuSxLu/16ZZ1exukxASb0JX1RyzIgW+g7OGyv+l6oskEElE4LGylb3rklvBEqbOMKCgKhgYJqPtgcBIVgEgM63tllXrH6ckYxUpQKHS6v7FEKNTU7x+9N4iQrhun6UIsF1ceLl3xFCCPEPIybA0ul0mDozHY0/DsenD80eVGV2tYIhMhg4ctmzAZa+meFUDdDYCVw3hmFcmHDSzRMmRjC0FuRj5YIMyGWhyExPc/k96ejmuNwinCBMNp/I89X2IAAoQ4QTjG0VZ+1WrH4+ZwKOltVgdf4J7Cosxer9hfjEUI0b7t/kt/lXFkEShgenSDw2/4QQQnxrRARYlsrsY5rKsWGuBvFGoTK7qwEFYwzTo4AT1Z4NsL5rkuCLSo4JYUCczDOrVxYdJ/KBnaugainHI7M1UNRVuPyelDUL+VeXW4TE6yQFQ5jI25iDESBhiA4BYjSTcUjf07pn6YQxSI0Jw8mqVjx1TI9jUOHOHW9i+uIsJMl9NlyX+VMLH0IIIeIaEQFW775+u5cLff22b93i8nNcFyfBd41Ap8kzQRbnHMdqAlDRAmTECCtX48M99wX78V+34ka17XuS7vJ7YjBylJnrX2kUwKRwjw3TZYlyhnn3b8Inhmqs2V+IXYWlWLO/EMXVTbhjy+sY/1ojknYUYNqSbCTIhJwtQgghxFdGRIBVXFKCBWpln8rsRcWuV2a/WcXQaQI+v+qZAOtSM/Dhe28jcMtsHFgWht+v9GyDacP5EizSuPeelDYBpUYOeaCwPecP221qBcOUxdkY91hPcdaLoSqsfGYfUhYuh765p4+fv+dfEUIIGflGRICVmpKCw4aaIVVmnzeWIYAB71/yTID1tzfy8O2ONbiOVWD9HA1iGge3jTlYk1NScFg/+Peky8RR3tyTfxUcwJAg88gQB8XSU3HWEi3aNxTgf4404f49JzFlURautgLt3T2FRf0hICSEEDK6jYgA67ENm3BU37N1tPZAIY7qq7Fh02aXn0MayDAuDDhW6ZkA6+Vntwhbdlnub2MOxubNm3HUUI3Vg3xPXn5DhxdXZ6Dih+GoeGI2rn6ajwA/yBWKkTLIAnuCqDKj8LfByPHGeeG0Y6JcKHswHPKvCCGEjGwjIsBqm5kFPLgXtRFCZfaWaPf6+mUoGb6uEz/AqmrluFp61u0tO3dotVos2ZKLL6DCbwtce090Oh0evPtOTGyvwIa5GkzvqsAf/3uFR7cyB0OtYIgIBsKDhDpdn17l+PtZjqhg4P4UhkAJo/wrQgghfmFEVHL/a7EJkxZm45uXVg7peebFM+Re5Khu44iRivclXVzPERQ/GYf1FXggXQPGmFcaTC9brsVHcVn46TSGF64feKothwV221RGX3ugENu3brFWvfcltQIormdQKzi+qAQ4OG6IAxYnMEgY5V8RQgjxH8N+BauqlePzSuAHk4Z+KTeayyZ8I/Iq1vEqjtbbNuCooRprD7i/jTlYM6IBBiFpvb594GsS47CAJ1mqn6dGMsgDgdUTGG5SSazBFUD5V4QQQvzDsA+wdhaZwAD819ShX8rUSAYJgAIR87AaO4Tq8iwjG9ptu60Npt3dxhyMRLkEUSFAVRuHwTjw/SdMGvphAU+KlwFBEmBaFMPDMxgmR9gHU5R/RQghxF8M2wBLp9MhMz0Nv7lOAfn22fj43fwhP2dwAEO8DPiyVrwA60IT8G09h0YBrMlehhOFp9Dc0ooThac8GlwBQKxU+FPVJuQsDeSun2/EUf3gE+O9JUDSU/neUWsZyr8ihBDiL4ZlgGWp3K6oK8ev5mqQZnK9SvlAJoYLOVNi+abWhIuNwJRIhtRwzzeTthUZAowJBapahdN2A9Es0GLqE7k4DqHOVH2E51fZBkvdzwoV5V8RQgjxF8MywBKSsWOxewiV252ZGsVwsUmEQZr9+xKHCcDcMUB8qGdb8fQmYQwTwhgaO4UVrJYu56/POUeZkaNlRhbGP1OAzZ8bcepUoV8FV4BwktAZyr8ihBDiL4ZlgCUkY0eLmoxt2XL8v5vDUL9xNnbn5g15nE0dHCeqOBLlwPVxvnmrZ8cK71GZEdD3EzheaQWMnRyXzBXRE+XML+pf9ZakEBL3e6P8K0IIIf5kWAZYYlRut9Wz5ViBx67R4DpWgR+sunPIW47fNXCcawRSIximRfkmWJkaCYQFAYZm3u82ob6Jo7wF6OKAJsx/t9tCAhjiQvv+nvKvCCGE+JNhGWD9auPQK7fb6tlyTMePM5Lxjyxxthyff10H/uRsfJwThu/fmI633357SM/njjGhEqgV5hWsfk4SlppXuEIkwNhQILmfrThfc7RN6K8BISGEkNFpWAZYwbOFyu1Xw8QpeeCpLcfdj6zA9ebeg4q6Cqxbt87rVdFjpEL9qPIWoKzZhI7uvqtYnAurW3ojh1oBBEsYVH683dY7wIqXqGKi4gAAERxJREFUARkx/hsQEkIIGX2GZYC17ZQJExZko+SrQlFKHqSmpOCIodZuy/GQfmj1n7Y8+STm2fUeTPdo70FnYqTClp+JA5eMQFlz3/tUtQFNnRxllvwrhdB2xl9ZGj8zCMVh70+RICrEf8dLCCFk9Bl2AdY3tSYcqwIenSlxWAvJHRs2bcbH+iprlfU1+wvxiWFo9Z9KSkqw0Iu9B50JkDCkRghbf4Zmx+Ua9EaOimag02TOv1J4dYiDFh7MkBzGsG6yBDepJH6ZjE8IIWR0G3YB1pOnTIgKBn6UIt7QtVot8vLyrFXWDTIV8F+5WPr9LLefM0aTgkN6/6iKHieTIEkB8zag/W0XGjkOVQi/D5YI223DodzBvZMYNMNgnIQQQkanYRNg6XQ6zJqVhn23h0G6/Rq8c2DoldttabVaa5X1Nw4XAulZQyo4Gr9yAz7xcu9BZxLlQJKcocwIGJpM6DZxcM6x9eV9mH9NOjbNVeDzhzIR9W0+QiQMiX6cf2Uh1uolIYQQ4gnDIsCylFGIaCjHhrkaTGovF61yuyOWkgonqtwLsDq6OWqmZiH+4b0wRvUk4r/66qs+KdypVjCoFUC7CahoEdr3PPxiHjbdvwLj28uxfo4GGaYKXP3dKlz9LN+v868IIYSQ4WBYBFiWMgp7PFC53ZGIYIaoYOBUjXsB1oUmDkMzMOdWLQpP9fQeXLZsmcgjdU2CDNAoAAkT8rD+ed6E3Be2Yp4mxvqe/iMrAzeqY/DBX7b5ZIyEEELISDIsAixPlFEYyPhw4Ns69wKsDy5xdJqABWP9YyUoUMKgVkiQIAPKjBzdHKjWl2Ch2j4Jf5FGCcM57ybhE0IIISPRsAiwHJVR8HTCeEoEw7lG9x778RWOYAmwROUfARYglDaQn8lHySOz8dQNCkgYw+HeSfhDLE1BCCGEEEGgrwfgig2bNiMnJwdrDxRioVqJw4YaHNVXQ/f8Lo+95sxohn9e4Ojo5oNqwdJt4viyRijYOSHcfwKsokP5KNmxCjeqY7Bojgb7iitw1CBUw1+kUeKQvgafGKqh+73n3lNCCCFktBgWK1i9yygMtXK7K2bHMpi40E9wMPTmUgjTohjkQf4TYL3y3BbM08TgH+bCp++vvg6pMWE4XtmMbcf0KAlRYcvL+3yShE8IIYSMNMNiBQsQgiytVuu115sZLQRHJ6s5pkW7/rgPLnF0mID5fpJ/ZXH2bAkema2xy7laMSUB247pkbm3EbclMtyTGuDjURJCCCEjw7BYwfKF2FAGRRBQWD24xx26zBHkZ/lXgCWPzT7n6qC+BlHqFNyeJEFIAINK5uNBEkIIISPEsFnB8oVkBfD1IE4SmjjHqRoOtRyYFOFfAZYlj8025+pTQzVWPrMTAJCkYNRyhhBCCBEJrWA5odPpUPH4bBzWhiEzPc2loqZlRo7SJiH/ShboX8GKVqvFK2/sw4VQFZ4+rsfFUBVWPrMPUxYJ7YD8vf8gIYQQMpzQCpYDlsrx8zQxWDhXgyMGoXJ8Xl5ev3lgH5ab86/i/Su4slh3Vw6qpmShuavvbcOh/yAhhBAyXNAKlgPuVo4/WMERyICb/Cz/ylaSou/YgiWAahj0HySEEEKGCwqwHHCncjznHEf+pUPQttm4Pknu8raityU5CKTUCgYJNU8mhBBCREMBlgPuVI7/2z/yUP67VZhtqsAjszVQ1FV4tCG1u9S9VrACGJCmpOCKEEIIERPlYDnQu3L8QfOJu/4qx2/fthU3qmOwJysDjDE8kK7B2gOF2L51i1frdw0kQQYEMqCLA4lyhuUahthQCrAIIYQQMdEKlgO2leN3HNfjC6jw4u7+q5yXXyjBIo3Sqw2p3REgYUgOY1iayHBfCgVXhBBCiCdQgOWEVqvFicJTOH6pGd0bC9A5K8vpfTnnCI6f3Ld5socbUrtr9QSGa+Mk1mCQEEIIIeKiAGsAM5UMEcHAkcvOC45ebeXouH0DjhqqsfZAIXYVlmLtgUIc1Vdjw6bNXhyta6igKCGEEOJZFGANgDGGqZEMp2qcB1gflnN0pmXj3ufe9GpDakIIIYT4J0pyd8G1ccAfvgZq20yIlvaNST8q5whgwKPrtJj+8J0+GCEhhBBC/AmtYLngFpUEJg68V+b49oJqjiQ5MNnP+g8SQgghxDcowHLBonggRAJ8VGHqc1tlqwkXGoEZ0QzBARRgEUIIIcTBFuEtt9yCK1euQCKRICwsDC+88ALUajWWLFlivU9LSwsuXLiAyspKREdHY+HChdDr9YiIiAAA3HvvvfjFL37hvavwsJBACSaGd+Nkdd88rHcMHK3dwAI/7T9ICCGEEO/rE2Dl5uYiMjISAJCfn49169bhyy+/xOnTp633efbZZ3HkyBFER0dbf/f8888jK8t5KYPhLjOGIfciR3u3CSEBEnSaON41cLx5Qci/ui2RAixCCCGECPpsEVqCKwBoaGhwWCvp5Zdfxn333efZkfmZRSqGtm7gyGWgqpVjV5EJp2o49EaORDkwKZICLEIIIYQIHOZg3XPPPUhKSsKmTZvw+uuv29322Wefoa6uDsuWLbP7/eOPP44ZM2Zg1apVuHDhQr8vajQa0djYaP3T3t4+xMvwvNsTGRiAPedM2FVsQlUbYOIceiMwM5ohiGpLEUIIIcTMYZmG1157DQDw6quvYv369Xj33Xett7388su45557EBjY89DXX38dSUlJ4Jzjz3/+M5YtW4Zvv/3W6YsuWLDA7ufHHnsM69evH9KFeFoggIRQOU5WcowN7EJtO8N7lwPR0sUwL7oNtbWdAz5HXV2d5wdKfI7meXSgeR4daJ5HD3fn2jZdyhbjlt4uToSGhuLSpUtQKpUwGo2Ij49HQUEBUvtpASOVSlFeXg6lUmn3+8LCQmRmZuLIkSNIS0uz/j4kJAQhISGDuR6fuOP9Lhy9wjEnluFgBUd4MHCHhuHXGQFIDht4Bau2ttbpRJCRg+Z5dKB5Hh1onkcPsefabgWrvr4eLS0tSEhIAADs378fSqXS+oJ79+7FrFmz7IKrrq4u1NTUIC4uDgCQl5eHuLi4PsGVLYVCgfDwcNEuwltuiGN428Dxn3KOuWOAJQkMskCGJLmvR0YIIYQQf2IXYDU0NGDFihVobW2FRCJBbGws3nnnHWui+8svv4wHHnjA7gna29vxve99D+3t7ZBIJIiJicFbb73lvSvwojvHMbxxHsiIYVArGIoO5uP4K1uxvbQEqSkp2LBpM7Rara+HSQghhBAfswuwNBoNjh8/7vTOn332WZ/fyeVynDhxQvyR+aHx4Qx3TZCgrRsoOpiP3MdWYL4mBj+YrcERQwVycnKQl5dHQRYhhBAyylEl90FgjCFJwRAZDHz9+lbM18Ri9/IM/DgjGbuXp2OeJgbbt27x9TAJIYQQ4mMUYA3S4gSGh6ZJUHquBAvU0dbtU8YYFqqVKCou9vEICSGEEOJrFGANgk6nw7Ib0hGhkEHCGI7oa2A5hMk5x2FDDab0c7qSEEIIIaODwzpYpC+dToecnBzM18Tikdka5BVfxseGaqzdX4iFGiUOG2pwVF8N3fO7fD1UQgghhPgYBVgu2rblSXPOVToYY3ggXYNb3/gcp6ubcbLKiCmpqdA9vwvZ2dm+HiohhBBCfIwCLBcVl5Tgkdkau5yrO1MT8NwJPZpbWn08OkIIIYT4E8rBclFqSgqOGGop54oQQgghA6IVLBdt2LQZOTk5WHugEAvVlHNFCCGEEOdoBctFWq0WeXl5aIlOxHMn9GiJToROp6OcK0IIIYT0QStYg6DVaqlKOyGEEEIGRCtYXtLe3o4dO3agvb3d10MhHkTzPDrQPI8ONM+jhyfmmnFL1rYXFBYWIjMzEydPnkRGRoa3XtYvNDY2IiIiAg0NDQgPD/f1cIiH0DyPDjTPowPN8+jhibmmFSxCCCGEEJFRgEUIIYQQIjKvJrm3tgoFOYuKirz5sn7BaDQCAE6fPg2FQuHj0RBPoXkeHWieRwea59FjqHOdmpoKmUxm9zuv5mDt2bMHd999t7dejhBCCCHE4xzllns1wKqursb777+P5ORkhIaGeutlCSGEEEI8xucrWIQQQgghowEluRNCCCGEiIwCLEIIIYQQkVGANYC2tjZkZWVh8uTJmDVrFm6++WacO3cOAFBZWYmlS5di0qRJmD59Oj7++GPr4/q7beHChRg3bhzS0tKQlpaG559/3unrv/POO0hNTcWkSZOg1WrR2NjouYsdxTwxz3PnzrXO8fTp08EYw5kzZxy+PmMMM2bMsN7/6NGjnr3gUczdud6+fTtSUlIgkUiwf/9+u+fs73G90WfaOzwxzz/84Q+tz3fDDTegoKDA6evTZ9o7PDHPon1Hc9Kv1tZW/q9//YubTCbOOed//OMf+YIFCzjnnP/whz/kv/71rznnnB8/fpyrVCre0dEx4G0LFizg+fn5A752U1MTHzNmDC8qKuKcc/6zn/2MP/rooyJeHbHwxDzbevPNN/n06dOdvj4AXldXJ94FEafcnetjx47x8+fPO/z8uvrvAX2mvccT83zgwAHe2dnJOef87bff5hqNxunr02faOzwxz2J9R1OANUgFBQXWD5VcLueXL1+23nbNNdfwDz74YMDbXJ283Nxcfuutt1p//uabb7hKpRLjMsgAxJhnW0uXLuXPP/+809ej/xj7jqtzbeHo8+vqvwf0mfYdMebZVlVVFQ8MDLQGXL3RZ9o3xJhnsb6jaYtwkP7whz9g+fLlqKmpQWdnJ8aOHWu9LTk5GQaDod/bLB5//HHMmDEDq1atwoULFxy+lsFggEajsXuOy5cvo6urywNXRmyJNc8AUFZWhiNHjgxYA27JkiWYNWsWHn74YTQ3N4t7QcQpV+a6P4N5HH2mfWeo8+zo+W6//XYEBjqv102fae8Ta57F+I6mAGsQtm/fjnPnzuGpp54a0vO8/vrrKC4uxpkzZzBv3jwsW7ZMpBESMYg1zxZ///vfsWzZMsTExDi9j16vx8mTJ/HZZ5+hqqoKv/zlL0V5bdI/seea+Cex53n37t3Izc3Frl27nN6HPtPe52/f0RRguejZZ5+FTqfDe++9B5lMBqVSicDAQFy5csV6n9LSUqjV6n5vA4CkpCQAQhLkQw89hAsXLqCmpqbPa6rVauj1ervniI+P7/f/mMjQiDnPAMA5xyuvvIL77ruv39e1PEYul+OnP/0pJcR6wWDmuj+DeRx9pr1PrHm22Lt3L37zm9/ggw8+QFxcnNP70Wfau8ScZ7G+oynAcsHvfvc7vPHGG/jggw8QGRlp/f2KFSuwc+dOAEBBQQHKy8uxYMGCfm/r6urC1atXrc+Rl5eHuLg4KJXKPq+7dOlSFBYWori4GADw4osv4q677vLYdY52Ys6zxcGDB9HV1YWbb77Z6evW1dWhpaUFAGAymbB3716kp6eLfn2khztz3R9XH0efae8Se55zc3OxceNGfPjhh/1+UdNn2rvEnGdRv6PdyCEbVcrKyjgAPn78eD5r1iw+a9YsPmfOHM4551euXOE333wznzhxIp86dSo/ePCg9XHObjMajTwzM5NPnz6dz5w5ky9evJifPn3a+rhNmzbxl156yfrzgQMHeEpKCp8wYQJfvnw5r6+v99KVjy5iz7PF6tWr+ebNm/u83ksvvcQ3bdrEOef8s88+4zNmzOAzZ87kU6dO5XfffTevqanx4NWObu7O9ZYtW7hKpeLBwcFcqVRylUrFKysrB3wcfaZ9wxPzHBgYyBMTE63PN2vWLF5dXc05p8+0r4g9z2J+R1OrHEIIIYQQkdEWISGEEEKIyCjAIoQQQggRGQVYhBBCCCEi+/9SLw3eTwt7QgAAAABJRU5ErkJggg==" }, "execution_count": 2, "metadata": {}, "output_type": "execute_result" } ], "source": [ "using Plots\n", "pyplot()\n", "\n", "plot(xtest,μ,ribbon=Σ, title=\"Time series prediction\",label=\"95% predictive confidence region\")\n", "scatter!(xtest,ytest,label=\"Observations\")" ] } ], "metadata": { "kernelspec": { "display_name": "Julia 1.0.0", "language": "julia", "name": "julia-1.0" }, "language_info": { "file_extension": ".jl", "mimetype": "application/julia", "name": "julia", "version": "1.0.0" } }, "nbformat": 4, "nbformat_minor": 2 }