{
"cells": [
{
"cell_type": "markdown",
"metadata": {},
"source": [
"# Simulating a model with basico\n",
"First some jupyter magic for plotting and convenience"
]
},
{
"cell_type": "code",
"execution_count": 1,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Using matplotlib backend: Qt5Agg\n",
"Populating the interactive namespace from numpy and matplotlib\n"
]
}
],
"source": [
"%pylab\n",
"%matplotlib inline\n",
"import sys\n",
"if not '../..' in sys.path:\n",
" sys.path.append('../..')"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Now import basico"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {},
"outputs": [],
"source": [
"from basico import *"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"now we are ready to load a model, just adjust the file_name variable, to match yours. \n",
"The file can be a COPASI or SBML file. For this example, we use the brusselator model, that is distributed with the package. "
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {},
"outputs": [],
"source": [
"file_name = get_examples('brusselator')[0]"
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {},
"outputs": [],
"source": [
"model = load_model(file_name)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"now we are ready to simulate. Calling `run_time_course` will run the simulation as specified in the COPASI file and return a pandas dataframe for it. "
]
},
{
"cell_type": "code",
"execution_count": 5,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"
"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"df = run_time_course()\n",
"df.plot()"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## The run_time_course command\n",
"you can change different options for the time course by adding named parameters into the `run_time_course_command`. Supported are: \n",
"\n",
"* `model`: incase you want to use another model than the one last loaded\n",
"* `scheduled`: to mark the model as scheduled\n",
"* `update_model`: to update the initial state after the simulation is run\n",
"* `duration`: to specify how long the simulation is run\n",
"* `automatic`: in case you would like automatic step size being used\n",
"* `output_event`: in case you would like to have the event values before and after the event hit listed\n",
"* `start_time`: to change the start time\n",
"* `step_number` or `intervals`: to overwrite the number of steps being used\n",
"* `method`: a method name to use for the simulation.\n"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"so lets run two simulations that will be different slightly, as we will use the `update_model` flag:"
]
},
{
"cell_type": "code",
"execution_count": 7,
"metadata": {},
"outputs": [],
"source": [
"df1 = run_time_course(update_model=True)\n",
"df2 = run_time_course(update_model=True)"
]
},
{
"cell_type": "code",
"execution_count": 8,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"(, )"
]
},
"execution_count": 8,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXQAAAEGCAYAAAB1iW6ZAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/d3fzzAAAACXBIWXMAAAsTAAALEwEAmpwYAABNiUlEQVR4nO2dd5hcV3nwf+dO266VtLuyerOae5FsywU3DAYDDgRjQxwMDjgxIQESkpjwfR+kEQiE0Gxsy5VebIohLhh3W12yZDVLVu/a1a6k1bZp93x/nHtn7szO7M5qZ+beO3N+z6Nndmdn7n11nnPf+963CiklGo1Go/E/htsCaDQajaY4aIWu0Wg0FYJW6BqNRlMhaIWu0Wg0FYJW6BqNRlMhBEtx0JaWFjljxoxSHFqj0WgqkjVr1hyVUraO5hglUegzZsxg9erVpTi0RqPRVCRCiD2jPYZ2uWg0Gk2FoBW6RqPRVAhaoWs0Gk2FUBIfei7i8Tj79+9nYGCgXKccMTU1NUyZMoVQKOS2KBqNRjNiyqbQ9+/fT2NjIzNmzEAIUa7TFoyUks7OTvbv38/MmTPdFkej0WhGTNlcLgMDA4wfP96TyhxACMH48eM9/QSh0Wg0Q1FWH7pXlbmN1+XTaDSaodBB0VJy8gi88Us4+pbbkviPvi5Y/TAMnHBbEv9hmtCxDTb+Cno73ZbGf3TugNUPQTLutiQjpmoU+r59+5g5cyZdXV0AHDt2jJkzZ7Jnz6hz+QcjJTz9BfjvufCrT8DvPlP8c1Qye5fDvVfA7z8L914O+1a6LZF/iPbAD94Hdy+Cxz4OL33VbYn8xeqH1J77/efgJzdD9KTbEo2IqlHoU6dO5c477+Suu+4C4K677uKOO+5g+vTpxT2RlPDMF2H5PXDBbXDBR2HPUuhpL+55KpV9K+HR90IgBDfeo9bzlx93Wyp/EOuFn3xI7bfr/g2mXwbbn3NbKv+wcolS5FMvhnf8O+x8EX77abelGhFly3LxAp/73Oe48MIL+da3vsWrr77Kd7/73eKfZNUDsPxuuPiv4PqvwpFNsPYH8Ob/wkKtmIbkxH742Z9B02T45PNQNw4634Kl33NbMu8jJTzxt7B3GXxgCZz9QQiE4el/gmO7YewMtyX0Ntufg6f+CebdADf/CAxD3Rg7d7gt2YhwRaH/y+82sflgd1GPecakJr703jOH/EwoFOLrX/86119/PX/4wx8Ih8NFlYFDb8Az/wxz3gHv/E8QAiacCWNnwpbfaYU+FMkE/PJjEO+H236nlDlAqA7MuPJnBnR9QF7W/gA2PgbX/B+lzAFmX6NedzwPC293Tzavc/IwPP4X0LYAPnC/UuYAwQgk/JX1VjUuF5unnnqKiRMnsnHjxuIeON4Pj90OtePgT76f3hRCwBnvg10vQf+x4p6zknjtW7B/Fbz3W9A2P/1+qE69xvvckMofHN0OT/0jzLoKLv+79Pstc6BpilLomtxIqdwq8QG46RGINKT/FohAMuqaaKeCKxb6cJZ0qVi3bh3PPvssy5cv5/LLL+eWW25h4sSJxTn4S19T7oE//zXUt2T+bd4N8Nq3YdcrSrlrMjm8AV78Kpz5/rR1aROqVa+xPqgZU37ZvI5pwhN/o6zJ998HRiD9NyHg9Gtg02/VE1CgqjyshbH2Udj+LLzr6+oG6CQYhkTMHblOkaqx0KWU3HnnnXzrW99i2rRp/MM//AOf//zni3Pwg+vgte/A+bemH3OdNFgtjmO9xTlfJWEmle+3thlu+Obgv4fr1au20HOz5iHYuxTe+RVoPG3w32deCdET0L65/LJ5nZNH4Nn/BzOugEWfGPx3H1roVaPQlyxZwrRp07juuusA+NSnPsWbb77JSy+9NLoDm6aKjNe3qMh4LoRlNUlzdOeqRFY/BAfXKoVk+82daJdLfk4egT/+C8y6Gs77s9yfqW1Wr/H+sonlG575Z7Uu7/mftIvUSTDiOwu9ap7B7rjjDu64447U74FAgDVr1oz+wOt/ohTS+++H2rG5PyOszSKToz9fJdHTAc/9m7Iiz74p92dSCl0rpEE8+/9U0O6G/1bulVxoYyI3O19SQeSrvjDY1WITCGsLvaoY6FYW0pSL4JwP5f9cSqHriyqDF/4d4r3w7m/kV0hhS6Frd1Ume5bCGz+DS/8Gxs/O/zltTAzGTKrCv+ZpcNln838uGAEzoZ7CfULVWOgl4eX/gt52+MjP8yskSAeqtEJPc3ijSrW76C+hdW7+z9lBUW2hpzFNePoula9/xd8P/VltTAxm7Q+gfRPc9CiEavJ/LhhRr8koGLXlkW2UaAv9VDm6HZbfC+fdCpMvGPqz9kVlaisJsKpp/xkiTXDlPw792ZAOig5iwy/h0Hq49kvpoHE+tDGRycAJeP7fYdqlcMaNQ382YCn0hH/cLlqhnyrPfAGCNXDt/xv+sykrSZZWJr+w7WmVl3/1P+cOhDpJpS1qlwugnlSe+1eYeF7+uIMTbUxk8vI3oK8Trv/K0E/VoNIWAZL+CYxqhX4q7HgB3voDXPkP0Dhh+M/rx940iZjqddMyt7DqxVTaona5AKpHUPd+lVGVKzMjG21MpOnaCSvuhfM+ApPOH/7zlWqhCyE+J4TYJITYKIT4qRBiCMdThWOa8Mcvw5ipyv9bCDowlWb1g9C1QymkQkr5U1ku2kKnpwNe+R+Y926YeUVh39FZLmn++GUwQnDN/y3s88EKVOhCiMnA3wILpZRnAQHgllILVkyklFx++eU89dRTqfd+8YtfcP3114/8YJt/A4fWKXfBUAEVJ9pCV0RPwstfh5lvU/1uCiEYAYS20AFe/E8VS7juXwv/ju1WqHZj4sAa2PxbuPTT0FRgdXjAdrn4R6EXmuUSBGqFEHGgDjhYOpGKjxCCe++9l5tuuomrr76aZDLJF7/4RZ5++umRHSgZh+f/DdrOgHNuLvx7OjClWHa38l9e++Xh/Zc2Qii3S6zKg6KdO2DNI6rBW7686Vzovad47l9Vn6XFI2iH60MLfViFLqU8IIT4BrAX6Af+IKX8Q/bnhBB3AHcATJs2rdhyjpqzzjqL9773vXzta1+jt7eXj370o8yePUT+bi7WPqr8cB/+eWbPjOHQgSk1OWfp92D+e2DKhSP7bqhWZ7m88BWlYN42TFZQNnrvqb7mO19U1cg1TYV/L+C/oOiwCl0IMRa4EZgJHAd+KYS4VUr5I+fnpJT3A/cDLFy4cOgIzFN3qYZMxeS0s+FdQ09n+dKXvsQFF1xAOBxm9erVIzt+rBde+i+YthjmvnNk39UuF3j1m8oPXqj/0kmorroV+uGNsPFxuPxzhQXhnVT73pNSWedNU2DhX4zsu5VooQNvB3ZJKTsAhBC/Ai4FfjTktzxIfX09N998Mw0NDUQikZF9efn3oecIfOiHhbsLbFKBqSrNNDixX02DOffDma1xCyVcX91piy/8h8rZv+xvR/7dag+Kvvl75T9/3/cKj3nZ2FkulWSho1wtlwgh6lAul2uBEZq3WQxjSZcSwzAwCkn3cjJwApZ+F+a8E6ZdPPKTVnuWy4tfBSRcddepfT9UW71B0X2rYOuTanBFvl5BQ1HNFrqZVL2CWuYqY2Kk2HnoPrLQh9VsUsoVwGPAWmCD9Z37SyyXt1hxHwwcP3WFlMo0qMKLqnMHrPuJyjlvPsXYSjW7XJ7/N6hrgYvvPLXvV7NC3/g4HN0KV3/x1HrBBxyl/z6hoP+llPJLwJdKLIs36T8Oy74Hc981fIl/PoRQF1Y1XlSv/LfKN7/8c6d+jFAd9Bwunkx+YdfLqqL2+q9mTtIZCUaVKnQzqYbOTDgLFpziUJmUhV5ZLpeK4stf/vLIvrDiXuVyufoLozuxMKov06BzB6z/GVz8l7mHLxRKuK460xZf/Co0ToQLRzGLtlqzXDY+Dp3bVcxrpC5Wm6Dlc/eRha5L/4ei/zgsu0el2k08d3THqkYL3bbOL/vM6I4Tqqs+H/quV2DPa+rJZqTBPCfV6HJxWufz33Pqx0mV/vvHQtcKfSiW36PGd52q79yJCFTXRWVb5wtvH511DpZCr7Isl5e+Bg2nwQW3je441ZjlYlvnV/7TqVvn4GjOpS30nEiPp+1lyNfXpVIVF7xX5biPlmqz0ItlnUP1uVx2vwa7X4HLPzs66xyqL8OqWNY5OCz0gdHLVSbKptBramro7Oz0rFKXUtLZ2UlNjXUBLbsbot1wZRGsc6guhV5M6xyUhZ6MVo8f+KWvQsMEuPBjoz9WtZX+b3isONY5pJvH+cjlUrag6JQpU9i/fz8dHR3lOuWIqampYcqUKco6X3EvnPEncNpZxTm4UUVB0WJa55A5KDrSWJxjepU9y1R2yzu/ku4FPxpSQdEqUOjJhJoiVgzrHFR2WiDiK5dL2RR6KBRi5syZ5Trd6Hj5e6oysRi+c5tqsdCLldnixDmGrtIV+ktfg/rW0WW2OKmmoGgxMluyCUZ8ZaHroGg2/cdVmfoZN0LbguIdt1qCosW2ziE95KLSy//3rYSdL8Clf5sejj1aqkWhF9s6twmEfWWha4Wezaolync+3PDdkSKMyg9Mde1S1vmFHy+edQ7VMyj6xa+qqtBFI2wiNRTVEhTd/BvLd/6PxbPOQVvovibWqzJb5rwDJp5T3GNXg8tl6XdUEO5UmkgNRTUMit6/GnY8B5f+zfCDn0dCNVjoUsIr34TW+TD/vcU9trbQfczaH6gBDMW2zkEpukq+qE4ehtd/pOY1Nk0q7rFt90Mlu1xe/oYawLDoE8U9bjVkuWx7Bto3qSKsYlrnYFnoWqH7j0QMXvsOTL8cpl1S/OMLg/5YnJ+v2lv8Y3uBZXeDmSiu79ym0l0uhzfCtqfgkjtPvWdLPiq99F9KeOUbMGYanPWnxT9+IOyr9rlaodus/ymcPAhX/F1pji8EB7p6+afHN3C8zz8bpCD6j8Hqh+DMD8C4WcU/fsrlUqEW+qv/A+EGuOiTxT92yuXizfqPUbP7Vdi/Srn5Chk6PlK0he5Dkgl1UU08D2ZfU5pziADSspKiiQp7/F25BGI9o+uoOBSVbKF37YRNv1JFWKfS73w4Kr30/9VvQn0bnH9raY4frNEWuu/Y/Bs4tkv5zkc6jahQHEHRWCUpdDuQPPf64hVhZZNKW6zAoOhr3wYjBIv/ujTHT/Xir0CXy4G1sON5WPyp4hRh5SIQ1qX/vsKOkLfMK27+ajZOhZ6sIIW+5lHo7ypNINkmZaFXmELvPqSGf5x/a3HTPJ04evF7te3GKfPqNyEyZuSzQkeCTlv0GXaE/Iq/K36E3Ikjy6ViLPREVI3mm345TL2odOcJVqhCX/Y9FawsdppnNsJg7Z5ObvjOq6U9Tznp2Apbfg8X3wE1TaU7j05b9BF2hLy5RBFyJ8JAWI+9FaPQ3/h5aQPJNoah+rlUUtpiXxesfhjO/iCMnVHacwmDk31R9nRW0Pq9+i3l3774r0p7Hh0U9RG7X7Ei5J8tTYTcicPlEq8El4uZVBfVxHNLF0h2UmmDolfcp7J2ShVIdiICgFk5rr7je2HDL1Q3yvqW0p5Lpy36iFf+W7UpPe/PSn8uYaRSxyrCQt/8W+jaUdpAspNQfeW4XKInVTfPeTcUt19QPoSBMJPEkxLTrAA/+tLvAgIu/XTpz6UtdJ+wfw3sfBEWf3r0QwQKwdHLJep3SykVSJ5b/FLrfARCvrKUhmTNIzBwvPSuKhthpAKivrfSe9pVRfe5N8OYKaU/XyDiq31XvQr9tf+BmmZYWKQ2pcMhDITtcvG7hf7Ws3Bkg3JVlTKQ7MQIVEa1YyKm5tTOuAKmLCzPOQ3H3vO7Ql/+fWUxX1YGVxWoMXTaQvc4nTtUhHzRJ8rXX9sIpF0ufr+oXv0mjJkK53yofOcUgcrIpd74mAokX/bZ8p3T8XToa3ffwAlY9YBqbd1yennOGYiAGffNgJDqVOhLv6uCHRf/ZfnOWSlZLvtWwt5lylVV6kCyEyPgm4sqL1Kqvdd2Jpx+bfnOWyk1EGseUa2tL/9s+c6ZGhTtD7dL9Sn0ng5VzHHuLdDQVr7zVkql6GvfVq6qUpVa56MS+slv/yO0b1YtcssRSLYRAYTfA/KJGCy/F2a+DSadX77z2oOifZKLXn0KfeX96m576d+U97wi4H8/ZucOePN/LVdVkbsCDkcl+NBf+zY0TS59zUM2leBy2fi4clVdWuIirGyClkL3SbVodSn0WK+aSDT/BmiZU95zC5Gy0H3bnGvZ95Sb5aI7yn9uv/vQD6xVdQ+X3Jl+jC8XjiplX+69lKvqDDj97eU9d0qh+6OfS3Up9Nd/pFq9lvsuD2AEEPjYj+l0VTVOKP/5/W6hL/0ORJrggtvKf25hIPBxQH77c6o9R7ldVeBwuWgL3VskE8rCnHoxTLu4/OfPSFv0YXHHqiXKSllcZleVjZ+HbHftUoVYC28vbd+RfAiBMH3scln6bWicBGd9sPzntp+mfJK6WD0KfctvVclwKSbqFIJDoceSPrM0Y32q5/m8d0PrXHdk8LOFvuxudUMqdd+RfFil/+DD+M3BdbDrZbjkr8rvqgIdFPUkUqqA1Pg5MPdd7sjguKh8ZyWt+7FqkVvuQLITv2a59HYqV985N0PTRHdkcBoTftt7S78D4UbVt8UNUha6drl4h10vw6H1qvdDuSobs/HrRWUmlatq8kKYttg9Ofxqoa9aAol+D9wMfbj3ju2BTb+BC2+DmjHuyKAtdA+y9DtqTNU5t7gnQ4bLxUc+9C2/g2O7Vc/ucgeknPgxyyXWp9Jk514PbfPdk8MIOPaejxT68nvUnrvkTvdk0GmLHuPwRlXQcfFflqcJVz4MR6aBX6wkKdXNcNys0k5zKgQ/WujrfwJ9ne5kVTkRRirDyjdpi31dsPaHKhBajiZc+QjYlaLaQvcGS7+rWq8uKuGYqkJwlv77xUrasxQOrFHzLo2Au7L4LcvFTMJSy1U1/VJ3ZfGjy2X1Q6pfvJuuKnBY6Fqhu8+J/aoZ0gUfLc1E9ZEgjFT5tW+6LS79DtSNh3M/4rYkKvbhJwt9y+/U4HG3XVXgv/hNfEANAJl9bekGjxdKQPdy8Q4r7lOWiZs+OBvhs8Kijq2w7WlY9EkI17ktjf986MvuVqPl3HZVQVb8xgd7742fQ2976WetFkIlWuhCiGYhxGNCiDeFEFuEEC6mOxRItEdNpD/jRhg73W1p/GclLf2umtl40SfdlkThJx/6vlWwfyVc8in3XVWQWaXs9b1nmiqr6rSzYeaVbkuTznLxiUIPFvi5bwNPSyk/KIQIAx4w2YZh3Y8hegIu+Wu3JVE4AlOet5J6OuCNX8B5Hyn9zMZC8ZOFvvxuiIwpz2jDQhAGQiYAHxQWbf8jHN0G77/ffVcVONrn+kOhD2uhCyGagLcBDwJIKWNSyuMllmt0mEk12WTKIpi6yG1pFEbAP1kuqx9UG3ixR26G4J9+6Mf3wuYnVO50uTtS5kMYCHxS+r/8bmicCGe+321JFEErM84nFnohLpdZQAfwsBDidSHEA0KI+uwPCSHuEEKsFkKs7ujoKLqgI2LrUyog5SWFJASGH1wu8QE1FWbOO8rfkXIo/GKhr7hPvZZzeMpwOPqhezpt8cgmNef3ok+6U+afiwoMigaBC4DvSynPB3qBu7I/JKW8X0q5UEq5sLW1tchijpDl98CYaeUbYFwIfnG5bHwMejuU/9dL+CHLJXpSDTA+40Z3c6ezEQYGPkiZXX4PBGvhwjLN+S0EIZQx4fW9Z1GIQt8P7JdSrrB+fwyl4L3Jwddhz2vKQgoUGiIoAyKQstA968eUUg0wbjsTZl3ltjSZ+MFCf/3HakTa4k+7LUkmQnh/YlFPB7zxSzjvw1A3zm1pMjECYCbclqIghlXoUsrDwD4hxDzrrWuBzSWVajQsuwfCDXDBn7stSSZOC92rF9Wul1Tf6Uvu9EZAyonXs1zMpLIwp14MUy50W5pM/JDlYsdtvPZkCGAEvW9MWBRqwv4N8GMrw2Un4KFnIgfdB2HTr1TutFvNfPIhfFD6v+weqG+Fs29yW5LBeL3b4tYn4fgeeMe/uS3JYLxuTHg1bmPjI5dLQQpdSrkOWFhaUYrAyiWqkMhLASkbI4DhZR/60bfgrWfgyrvc7XmTD+HxLJdld0PzNG8UEmXjHBLtxb234ZfejNvY+CF+Y1E5laKxXtX/Yf4NMG6m29IMxuuFRcu/ryL6bve8yYfhYR/6gTWwd5kaYOGFQqJsHBa65+I3Uqq958W4jY0RrBwfum9Y/1MYOO6dQqJshIHhmOsopYda6PZ1qfU750PQ0Oa2NLkRHraSlt2jhjCc77G4jY0wUk+Hnktb3Pmiitss/pT34jY2fgjIW1SGQjdNdVFNugCmXeK2NLlxXFRSQsL0kEJf8wjE+7z7yAvetdBPHIDNv1EN4NyYF1oIRsC7NRDLrbiNG/NCC0Vb6GXmrT9A1w5VSOTZu3z6sRc89OibiKkhDLOugglnui1NfqzA1NIdR72zdqDWzqtxGxshvBmQ79imrt1Fn/Bm3MbGL1XKVIpCX343NE1WBR1exQikXC7goQtr82/g5CHvuqpsjAAg+ciS5fxx8xG3pVFEe2DNw7Dgvd5oAJcPx9Ohp4KiK76vml8t9GjcxqaS8tA9z6E31MzQi+6AQMhtafLjuKjAIwpdSpWdMX4OnP52t6UZGqGCjQFMTvTHXRbGYv1PYeCE9wqJshEe7CPU1wXrfgrn3AQNLleWD4f2oZeR5d+HUJ1qhuRlrKCoEB5KH9u7DA6tU4VEbg3PLhRLvgCmNwJ7pqn8v5MXwtSL3JZmaISBIT3WnGvNw2p4tpfjNjbah14mTh5WOazn3+r+RKLhsCzM2qBack9cWMvuVut27ofdlmR4rPUzMBmIe8Ba2vY0dO1U2RkeR2ZlWLlOIqZqRrwet7ExAvRHY/zp95dyrNfbTbr8rdBXPaDunBf/lduSDI9QS90QUkFb1y+srl3w5v+qRkhemEg0HEba5eIJC33Z3TBmKizwcNzGwqnQPTH+0C9xGxsjwMn+KGv2HGN7R4/b0gyJfxV6fEAVEs17F4yf7bY0w2Nl39QGLYXu9oW14j6lJL0ykWg4HD501y30Q2/AnletuI2HGsDlwURgCCsP3W1DQko1kahlrvfjNjYigEwql0s07oEb4hD4V6FvfAz6Or0xL7QQLAuzNqwUuqupdwMn4PUfwpkfgKZJ7skxEoy0y8V1C33lfSpu47UGcPkQKsMqEjSIJVwuatu7DA6t90fcxsYIIq2iNteNiWHwyYpmISWsuBfazoAZV7gtTWFYLpd6y+XiqlJ6/ccQ6/GF/zeFVyz03qOqzeu5t3g/bmNhIjAwqQ2rNYwnXVToy7+v1u2cW9yTYaQYAaQVFB1IaIVefPYshcMbVDGHVwuJsrEUuusuF9NUxTBTL4FJ57sjw6lgWXMG0t2b4ZpHVJvXizxcSJSFRKXM1oaUQnctfnN8H7z5e7jgNn/EbWyMYKrthHa5lIIV96q7/NkfcluSwrEszLqQWnLXrKTtz6rxfF6ubMyFF7JcknFY9SDMuhra5rsjwylgB0VTCt2tG+KqB9SrVxvA5UMYqbRFbaEXG7/e5VMWuvrVtYtqxb3QOElVN/oJL2S5bPkdnDzoj6wqB8pCl0TcVOjxflj7qOqG2jyt/OcfDUYQaZX+D2gLvcisWqJeF33CXTlGiqXQbQs9lnThTt+xFXY8D4tu93ZVbS5sC124qNBX3AtjZ6pBDD5CCkEAk9qQizUQG34J/cd8dzME1MQn20LXQdEiEuuDNY+qIQLNU92WZkRI20IP2PnALrhcVt6vemd4aQhvoRguB0UPrIV9K5Sryi/ZGRYmqvTfDoqW3YcupUqTnXAWTL+svOcuBo4RdK5nWA2Dv3bmhl9YPc99kqrowEQFQ2vsLJdyX1QDJ1TvjLM/CPUt5T13MRAuu1xW3q9m1Z73kfKfe5SY2Ba6Sy6XPa/BkY3+SmJwIgKIVFBUW+jFwb7Ln3Y2TFvstjQjxlborj32vv5jiPeqYhg/kspyMct/UfW0w8bHlTL32qzaApBCpS3WuJXlsuI+K4nBg7NqC8HRi1+7XIrFrpehfbPywfnwLi8tC7PGcKHrnZlUxTBTL4FJ55XvvMXETQt99cOQjPn2ZphyubhhoTuTGEK15TtvMTECCGlVimqXS5FYcR/Ujff2ZJMhyHa5lLVS9K1n4dhu/6UqOnHLh56IweoH4fTrvDmRvgCk7XIJu6DQU6mKPkticGIEU/OAtYVeDI7thq1PqmCelyebDIFtoYcNMESZL6qV9/kzVdGJcKn0f/NvoOeIP7MzLMzsPPRyZVjF+qxURf8lMWQgAoiUy0Vb6KNn5RKV9ue3ggQHplSWeQBJOGiUz4+ZSlX8C/+lKjpxy0JfcS+MPx1mX1O+cxYZKQ0MoXq5QBmNCT+nKjoxgimFHtWFRaMk2gNrf6jGy/mlkVQObJdLwJCEAkb5LqpUquLHynO+UuGGD33/ajiwRikkn6UqOjHtTp/W/bws6yel2nsTzoLpl5b+fKXESA8I0Rb6aHnjZxA94ctURSemlYcesCylsljofk9VdOLIckmasjwxiBX3QqRJNeLyMSYu9BHye6qiEyOIYfvQtYU+CuxUxUnnw5RFbkszKkxpj1BLEi6Xhe73VEUnKQtdZQmV3MrsPgSbfq2mYUUaS3uuEmM/HdYFy9hHKNVvyaepik5EAAPdnGv07Hgejm7zbaqiE/uxNwDKh15qhWSnKk5b7N9URSdGuvQfypBtsPohtYZ+GQAyBBK77YT6PVZqK/P4Xmsa1sf8m6roRFvoRWLFfVDfBme+321JRo1toRtC+dBL7jKwUxUrwTqHDB86lNhCT0SVQp97PYybVbrzlImkdZnXlGv84aoHAQEL/ZvEkIGhLfTR07kD3voDLLwdghG3pRk1qaAoZnksdL92VcyHkanQS2qhb3wc+o76O2/fQapKOVAGH3oqVfEGf6cqOjECGOXYd0XAuwp91QPqIl7ow0ZSOUgpdFGGtMWOrbDzBf+nKjpx5KFDCS0lO27TMk9Npa8AbJdLRC1haRV6paQqOhEBAro51yiI9aqA3hl/Ao2nuS1NUZCOLJeSpy1WSqqiE8MOKpfYl3lgDRxap3znPo/b2NguF9uYKFljuFRXxbP9n6roxAiW58mwCHhTob/xC5WqWAEBKZtkKstFpS2W7E7ff9xKVbzJ/6mKTrJ96KWy0FcugXCj71MVncjU06FZ2gyr3a9C+6bKSFV0YApVmBU0JAlTkijFDdEszo3CewpdSnVRnXY2TL3YbWmKRrpSVHW9K9mdfp2VqnhxhQRDbSwfeo3VT74kFnpPB2z6FZz3Yd+nKjpxGhMljd+svA9qx6m6hwoiidp7YyNqHUtijG17uiiH8Z5C37NU3eUvuqPi7vKgfMB14QD9pVDoZlK5W6YthonnFv/4bpI1k7UkFvraR1VXxUWV82QIjpRZIQmXKsOq0lIVHdguq+Za9VoSY2zl/UU5jPcU+qolUNPs266K+UhaFrrdJKk/VoJNUQldFfNhWej1qfL1Iq9fMqHa5M66ClrnFvfYLmM69l4kZJSmfH3VA4Dwdb+lfMStJ5zmGrUHB4ptoXdsg50vFuVQ3lLo3QfVIN7zb/XXAOgCSDqyXGrDJVLoK+9XqYrz31P8Y7tNaiarNfGp2Epp21PQvb/irHPIDIrWhkrwdBjvh7U/UKmKY6YU99gewDbGxtbYT4dFXr9VD4BRnGw0byn0NY8ot0EF3uXtfhqGNQqs6BfV0e2w4zmV5lkpqYpOLAu9Nqh+LboPfeX9MGaqKiaqMExHULQ2XIL4zcZfqVTFCkpicBK3fOhjatQ6FvUJJ3oS1v+0aMWTBSt0IURACPG6EOL3RTlzNomYeuSd846KqM7Lxtk+tzYUIFHsBlOrH1R3+QtuK94xvYTtQ7cUelEt9PY31USshbdDIFi843oEW6EbUsVv+or9dLhqCbTOhxlXFPe4HiFhqvUbE7Z86MU0Jt74OUS7i3YzHImF/hlgS1HOmostT0Bve8Xe5VMWukhPXy/ahZXK238fNE4ozjG9RraFXkwrc9UDEAjDBR8t3jE9RCrLRZQgfrN/DRx8XU0kqqAkBieJlIVe5KColLDyAZXAUKTmgwUpdCHEFOAG4IGinDUXK5fA2Jkw+9qSncJNbB+64RgFVrSNseGXKm+/Av2/KeyJTwE18aloqWMD3eqR96w/ray8fQdmxt4LFtfdt2oJhBvgnJuLd0yPYT8MNkWs+E2x9t7uV6Fji7pui3QzLNRC/xbwj0De/4kQ4g4hxGohxOqOjo6RSXHoDdi3XFnnPh4kMBTJLJcLUBxLyb7LTzgLpl0y+uN5FctCD4ki5/Gv/xnEeir2yRCyAvIhg75YojgH7u1U/vNzb4GapuIc04PEpdp7TWE7IF+kvbdqidViuHgZfcNqTyHEe4B2KeWaoT4npbxfSrlQSrmwtbV1ZFKsWgLBWjjvIyP7no9IOlwudcV0uexbAUc2VPQjL5DKcgmKIlbaSqn23qQLYPKFoz+eR3EWFtWFg8Vzubz+A0hG/T0AugDiljHWGCliUPTEAdjye5XRV8S8/ULM4cuA9wkhdgM/A64RQvyoaBL0dcEbv4RzPqTuVhVKOhdYWZhAcR59Vy6ByBi1fpWMZaEHhSyehb7rJdVvv1JaDOfBdrkIy91XlH1nJmHVQyoQ2rZg9MfzMDFTqcnGVKVoEdZvzSMgzaK3GB5WoUspvyClnCKlnAHcAjwvpby1aBKs+zEk+iv6kRccPnRpUhdWkb1RW0o97bD5t+rJJlw/WhE9jd3cLCTM4lnoK5dA3fiK6Lc/FNlB0XiyCBlWb/0BTuyt+OsW0hZ6Q6hIFnoiphT6nHfAuJmjlC4Tdx3WpqkyDKYtVr1bKpiUlSQcPvTRWkprHgUzXvGPvJCplIpioR/fB1ufVGmeoZoiSOhdUj50mXb3jXrvrVyiitjm3TBa8TxPykIP2wp9lGuXyugr/pPhiBS6lPJFKWXxyhC3/9GaqlP5d/lEyo9pUmvls47qokomYM3DMOtqaDm9GCJ6mgS2D71IFvrqh9TrwttHKZn3sd19QqQzrEb1dNi5w1HEVnl5+9nErTz0hnCRslxW3q9qbWZfM1rRBuGuhb7yfmg4DeZXyFSdIXD206hNuVxGkW2w9UnoPlAVN0OAmHVRBTGJjNZCjw+oRlzz3l05U3WGIFX6X6wMq1UVXsSWRdTaezUBSdAQo9t7h9arRIZFnyhJRp97Cr1zh7LQL/wYBMOuiVEuklml/zDai2pJxZaq58Ku1gsUI8tl06+hr7NqboZ2yqywOn3CKDKsYn2w7keVXcSWhW1MhA3b3TeKvbeytBl97in01Q+pzIVKmqozBEnVxjtToZ/qxujYmi5Vt7I/Kp24KUlKQbAYeeirlkDLXJh5ZfEE9DDJrF78AP3xU3w63PBLGKjwIrYsbIUekEnLmDjFvdfXBRseK2lGnzsKPdYLr/8QFrwPmia6IkK5cbbPrbF6ep+yy6XCS9VzEU+aJDEIoHzopzykYf8aNWauiNV5Xif9dCgdGVansH728JlKL2LLwna5CGmOzkIvQ0afOwrdvstXySMvpC8qgUQIceodF6Mn1Yi5Mz9QsaXquYgnJSYGgdFa6HapegWNmBsO59Nh2uVyCsZEtRSxZRFNWv9XM6H6yZ+KhV6mjL7yK/SMUvXFZT+9W6QtdLUZTrnr3Rs/h9jJqroZAiRSFvoofOi9R2Hj4xVfqp6NvfeQoyxqq5YitixslwsySSQYOLXSfzujr8QpxuVX6HuXq7t8BU1VL4REqrBImUs1p2Khp7qznVfRpeq5iDlcLjWhwKkp9LU/qMgRc8NhPx0i0xb6iAPyVVTEls1AykJPUhM6RWNi1RJomKDczCWk/Ap9lXWXP/umsp/aTezCGGH1NzulQQN7XlPd2arsZgiQcLhcIkFj5GtnJlW//RlXQNv80gjpUVIWupkOyI/46bCKitiyiToVevAUrttju9V4yAtuK3lGX3kVek87bH6iKu/ydqW1gbLQT8nlstLqznbWnxZZOu+TGRRVA0ISIylff+tZVapehQrJLmpDOgqLRqKUqqyILZuBbB/6SIOiqx9WzeXKkNFXXoX++g/VXb4KqvOySToaJIHlchmJQu8+BG8WvzubX0gFRR1ZQiN69F39oFXEVvml6tnYbSeQ6unGECN0uVRZEVs2A/Y2k8pCH1HaYiKq9N68d8GYySWRz0n5FLqZVA1pZlxRcVPVC8FZWARQO9JMDXveapG7s/kF20I3SI68OCb1yPvRypy3OgwJh0KyM6xG9HRYZUVs2WT70EdkoW9+QhWxlcmILZ9C3/4cHN9bldY5pCsd7aDoiFwuybjVne26ondn8wsJM53lUh9RudS90QJT79Y8qmIOF1ZHqXo2CUdQFBjZ1KJUEdvHq6aILZtMhT5CQ2z1g2oS26yrSyNcFuVT6KsfhPo2mF+83l5+ItvlMqI89C2/g57DVZed4SSWkJhSWei2Qu8pRKEnYuqRd+71MGZKiaX0Jk6XC0Bt2Ci8qC1VxFadN0OAAXupzMTIstOObIK9y6yK7vKo2vKc5fhe2PaMeuStgr4tuUjIdLUZqCyXgv2Yqx6EsTPg9LeXSDrvY1voBib14RFY6FuegN6OqnVVQTogj2nVQIQKtNBjvWpE35nvr6oitmz67aWSSRoiQXqjCaT1pD0kqx+CQETFvcpEeRR66pH3Y2U5nRexXS4jttA7tsGeV+HCj1fsvNVCSPvQTeoj6tG/txArc/VD6mZYglalfsGeiYmlhGoLdfdtfByi3VXrJrVxZrnUR4KYsoAhF9EeWP9zdTOsG1d6IS1KryESMVXQMecdVdGqNB9mqvw6fVH1x5PD3+nXPqpalZ73ZyWW0NvYWS6GNGlI+dCHUUrtW1TufpXfDFOqR6r1qi00w2rNI9A6H6ZeXCrRfEF/Iu1Db7CMiZPR+NBf2vALVdFd5pth6Xf51v9V0zmq+JEXVD+NpBSpx97acAAph0m9iw/Aup/A/HdDwwgHb1cYmRZ6gS6X1Q8r/28ZH3m9SMLM9KHXFTJX9NAbqonZhR+vuiK2bPoTltFlJh17b4j1k1LNW51wFky9qAwSpim9Ql/1IDRPg9OvLfmpvIwpwcRIB6YK6Yn+5u+hv6uqXVU2dqWoU6EPGRSN9cL6n8IZf1LV/l9It51gJPGbNY9AsKbq+rbkos+20GWyMGNi/2rV3mTh7WW/GZZWoXdsg92vKIVUpSlPNqaUuRX6UJbSmkegeTrMvKrk8nkdZaELhExSb+WhD2klbXhM+X8XVfeTIUDSTLsMoID4TawX3viFuhmW0f/rRaSUDNhLZSZS7r4hjYnVD6qOni7cDEur0Nc8rPy/51dP3+58JE2p0sdk2uUCQxTHHN1u3Qxvq2r/r008KZXLRZoEAwaRoDF0UHT1Q9B2RtX7fyE9td7pchkyKLrxV8r/q58MiSdlunWCWYCF3n9Mrd85N0OksUxSpimdpkjGVcrT/Buq3v8LSqFLIdKZBpaFnrdIYe0jYAThvOr2/9rEkyYmRipLyE4fy8mRzXBoncqdrnL/L2S2zwWoGc7lsuZhaJlXVUMs8hFNJNPdKh1B0bwW+qZfQzLq2vCZ0in07X9U/t9zP1yyU/gJUyoLM20lqTt9TkspEVXB0Hnvqpq5jcNh90MX1hNO/VAK/Y2fqZvh2R8so4TexY7ppfZeKEgsaeZubnZ4gwqGLtTBUIBYwnS0Hy4gKLrup9C6ACaeWyYJMymdQl//M6gbX/XBUJukKZEYjiwXawxdLgt9x/Oq/4N2VaWIJSWSAMJUSqg+EqQn10VlJpX/9/Trqj4YapNI5aGnXS6QZ+9tfBxEAM7WwVBQWWgSA4lI5aFDHpdL5w7Yv1INUHHpZlgahS6TsPUp1ea1Cpsh5UIFRUX6sXeoLJctv1M942ddVUYJvU0iaSKFkYpB1IcDuS+qXS/DyUNw7s1lltC7ZFvoNfmGXEiphljMvALqx5dRQu9ipxVLEVA+9PAQQdH1P1Vtcs9xb++VRqH3n1B+pHOqZ27jcKigaFohpYb1Zk9fTyZUu9J511dtm4RcxJMmpkg/4dRHgrmDohsfh0gTzH1XmSX0LomsLJe6fBlWRzZB104448ZyiudpYhkKPUHAENTlMyY2PwEzLnd18H1pFHq0GxonwuQLSnJ4P5I0US6XQXnoWX7MPa+pSPmC95ZbRE8TN6W6qKwbYs6gqJSw80WYdSWEasovpEfJznLJm2G15QlAVG0DvVzYvc+lEUytX05j4sQBOLpVufpcpEQKvUf1PddBlRRSSuWHy1Log6avb/kdBGthto49OIknLJdLykIPDA5Mde2EE/u0qyqLZMrlkpkyO8hC3/wETL8UGtrKKJ23SVVyCwNMda025Irf7HxRvc4uT5vcfJRGoZtxmPm2khzarySltFwGmVbSID/mtmdUIDlcV24RPU3ClCpYN1SWi31RzbyqnKJ5nkR2Hnqu+M3Jw2pe7TztqnKScrkYwSxjInvvvQD1rdB2ZrlFzKB0WS4zryjZof1IKsvFuqjCQYNw0MgMrvR0qLmX0xa7JKV3idlBUdORhx7LamO66yVomgzjZ7skpTfJ7uWS0+Vy8HX1OmVROUXzPKlxc5YPHaA+HMy8bk3TcvVd5XoRYGnOHgirlqWaFNlZLgBNNUG6Bxwb49B69epSDquXUVkugYygsikdbgPThF2vwMwrtasvi1SWi5ldA+HYewdfV26F084us3TeJmq3yTUy4zc9zuu2fZPque+BFs2lUegulLx6nWRWUA+gqSZE94CjDeehdep14jnlFc4HxJMylToGpCr2Un70IxtUIdusK90S0bPEzcwRdI01SqGfHMhS6K3zIVxfbvE8TcwqvhJGIH+G1a5X1KsHYjclUugNJTmsn0maDLLQG2uCmRfVoXUwbhbUjCm/gB4nnjSVBenwoYOjwMN2GejeLYNQnT7Te6+pRtWGpIwJKdX6TTrfLRE9S9pCD2YqdKfLpX0z1LVA0yQXJMykRAq9qSSH9TOmzPShAzTWhDiZYaGv1+6WPCiFnmklgaPA4+hbqt1r8zS3RPQsyVSGlVq7mpBBKCDSxkT3AeUy0Ap9ECkfuhFwZLkEMn3oR9+C1nkuSDeY0ih0I1iSw/oZ5XJJp90BNNUG6e63FHpfl5q9OvE8dwT0OImktPyY6aAoOCz0zu0wbnbVt2nOhWlKTJFeOyEEjTWh9N6zn260Qh+EnbYoAsGMp8OBuKMXztFt0DLHLREz0H1Zy4SZlYcO0BgJpa0kHRAdkrgdFLWrHe2e6DGHhd5yulvieZpkjr2XEZA/+Loywia4m3LnRVIKPcNCt4yJWBJ6O1XspmWuazI60Qq9TJjSstAdaXYZPnSt0IcknrLQ05kGYAVFEzE4thvGe8NK8hrZjeEAmmod7r6D61SHwFCtOwJ6mLRCD2akzIL1dHh0q/pgSyW7XDSDyO7lAuqi6o8nlX+4fbPKoa7yCTH5iCdNRFYvF7AuqmO71Lp6xEryGqaZ25hIuVy6dnrGZeA1ookk4aChLPRcAfmj29QHPbJ+wyp0IcRUIcQLQogtQohNQojPlEOwSsM0sS6qzCwXsNLHju1WGS6anCTMTAs9Iyh69C31Ie1yyUkyFZDPTplNqBvkiX26biQPsYRJJGhYWS6ZLpfU3gvWwJipboqZohALPQH8vZRyAXAJ8NdCiDNKK1blkcyT5QKoR99je2DsdLfE8zz9sSRGIJB67M2YK9ppKXTtcslJMocx0WRnWHUfUIpKK/ScRBMmkWAgZ4ZVbzSpLPTxc1yvELUZVgop5SEp5Vrr55PAFmByqQWrNHJmuVgWek/PSeg5DM0zXJLO+/TGEgQcmQYZc0WPboeGCVCj02VzkTtlNkh3v/VkCNqYyEM0PthCr3eOoTu6DVq94+ob0W1FCDEDOB9YkeNvdwghVgshVnd0dBRJvMrBlBIGFRYpCz3WsUu9oa2knEgp6Y8llUJ33BBTLXSPbtP+8yFIzbM1B8dvEp167w1FLGkr9MEps/19PerJ2kN7r2CFLoRoAB4HPiul7M7+u5TyfinlQinlwtZWPRQ6m2RWLjCkfejmsT3qDX1R5SSWNEmYkkAwmOEHTlXsdb4F47X/PB9mqu1EZtoiQPzoLuVOaJrilnieJhpXQVFnYZHtcgkc3wVIT+29ghS6ECKEUuY/llL+qrQiVSZJU1ql6+mLakytstCN47vVG/qxNyd2m9dgloVeHwli9h1XA0F0h8W85MpDt58Ok527YMwUCOhiwFxE7aBoRh8hS6F371MfGjvTLfEGUUiWiwAeBLZIKb9ZepEqEynJcVGpjRHq3guhOtVPWTOIXluhZ1noDZEA9f0H1C/N+maYj1T8xmmhp4yJPfrJcAhidlDU4UOPBA0ChiDcY+89b2S4QGEW+mXAnwPXCCHWWf/eXWK5Ko6kzOwWCOk7fW3vfnVR6bavOem3qkEDQWvguGl3DQxR139Qveehi8prpIvaBrtcgif36ifDIYgmkkRCRkbKrBCC+nCAur4DKmXRQ4bYsM9ZUspXAa1pRolpSqWwHRdVMGBQHw7Q0L8fpnonsOI17Ba5oaC1XWUSMGiuC9E0cEi9py30vCRNCTlSZusYIDzQqS30IYgmTMbW2T70tDHWWBOiof+gyj/3kCHmjeTJKiCZw0oCaIwEaY4e1BfVEPQ5XS6QurDG1oUZGz8E4QaoHeuWeJ5GSokpydkYbqpoV7/ovZeXWMJUFnrW03VzXYim6CHPdffUCr1M5BpwATA10kvE7NcX1RD0x5XLJRSyXC7SVughTjPbMcdM8ZSV5CVMq9o/lw99qrDSi3X9Q16iOXzooIyJcfHDWqFXK7lcLgCzw53qB+3HzMsgl4tlKTXXhZksjhJv0Cl3+UjaGt0xHASgIRxkmqEt9OGIJpKEA5k+dIAJtUnGyG7PxW60Qi8T6dJ/mfH+dH1RDYudthgKOX3oMK4+zBTRQV+dLlzOh2ntt2wL3TAEs4MdRI063RBuCFIuFyMzZXZGwDLEPBa70Qq9TCRNBvnhAKZgKXSPbQwvYfc8D4cys1zGB/sZI/rojkx0SzTPk7LQcxoTR+kMTdTuqiFI56FnxiCmiKMAJBq9ZUxohV4mpLTKr7NcLhPNIxylGcJ17gjmA+ygaCiY6UNvTaqbYWfoNFfk8gNJp4WeZUxM5QiHA3rthiKfD/00eQSA7hr354g60Qq9TCTl4EpRgNbEIfaarcgs60mTpi+WIGCIwVkuscMAHDHa3BLN85gpH3pm6T9ScprZzkGh1y4fiaRJ0pTp0n+HD70l2U5UBukS3squ0gq9TKjS/8FZLmOjB9kj21KTUTSD6YslqQsF1JABSE8tGlBFRYfwTmGH10gHRbOeDnvaiRBlr6kVej5i1szQ7G6LAM3RQxyU4znWn8j3dVfQCr1MmDl6uZCM0xA9wj7ZSrc9DkwziL5okrpIID0A2rLQQ9376ZdhDsbqXZTO22S4XJzGhNU2d2dS3wzzEY07FLpI9+IHqO8/yH7ZSldvzC3xcqIVepnI6XI5sQ8Dk32yLT1bVDOIvniSunBQrR+kFdPxPRw22jxnJXmJlA7KdrlYCv2tmM5wyYf91BwOBjK6LQJEevayX7ZyvE8r9KokNYLOHHxR7TUncMxjd3ov0RdNUBe2psZAeg27dtERnOi5i8pL2Bb6IGPiuGrZvDU6Lu1n12QQSzhdLg53aV8Xgf4udshJdPV668laK/QykdNCt/qg75OtdJyMuiSZ9+mLJZVCd/rQzSR0bqc9Mp1jWqHnxXQWFmUZE72RVgZkWD8d5iGaUAo8nYdurZM1w3afMclze08r9DKRTA05zryopBHiMONo1wo9L32xRKbLxUwqCzMZ5Xj9DI73ectK8hJ2UFQa2cbEbgbqVZVjR8+AG6J5npTLJWCkXVZSpmbYdtXO8NyTtVboZcDMKO7ICkw1T0MYAdpP6osqHzktdMtK6m2a7TkryUvYLhcxKCi6B9MqZjvSrY2JXNgKPRKy8tBBGRNH3wIjRF/dFM/tPa3Qy4Bdfk0OK0mMnU5LQ1i7XIZAKfSgw4duTVsHEs2zOdEf137gPNjrkjGCLhGF7gMEx80A4Ei3NiZykXK5BA117ULamBg/m+aGGo557OlQK/QykA5MBQYHpsbOoLUxol0uQ6BcLtkW+jaoa6F2TBumRKd95iFnUPT4XkBSO0GN7dN7LzfpLBfDYaGnZ9g214W1y6UaSaeOOQJTvUfVLMxxs2hrrKFdP/bmpS9m5aE7s1yOvgUtcxlbr9oBeC0f2CtkdFu0S/8PbwCgZvLZ1IcD2kLPQ0aWi733ElHo2gktcxhXF6ZLu1yqj5Qf03BU6x1ar15PO4e2xggdPVqh5yJpSqIJk7pQMPOxt2MrtM6luS4M4LlHX69g2w/CGZA/vEFZnG0LmNCkjYl8pHzoQcfTYddOZaW3zGVsXYgT/XFHAzT30Qq9DCQz/JiWlWQr9IlKoXf2RD21MbxCn9Vpsd5pofe0Q3+XdVEpha5z0XOTdrmIdLfFw29A63wIRmhriuiAfB6icacP3XK5tG9Rr+PnMLY+jJRwot87xoRW6GXADkwJpx/z0HrVMrd2LK2NEUwJndpKH4TdabHW6UPveFO9WlYSaAs9HzkHXBzeAKedDUBbY43OcslDRi8XO2U2tfdOZ1y9/XToHWNCK/QykDMwdfgNmHguAK2NNYAOTuXCVuj1ziyX9s3qtWVOyuWiLfTcpDOsLJfLySPQcwROOweACU0RjnQP6G6fOUj3cnGkLR5aDw0ToHZs2t3nofiNVuhlIJVSZ08OHzihfHET1UXV1hQB0KmLOeiNKpdLhoW+byXUNMOYaTTVBKkLB9h/rN89IT1MMvvp0AqI2ntvQlMN0YRJt64WHUQ6D91I7739q2DKIgDGWQq9Uyv06iLdwTT7ojoPgLZGpdC1L3Mw/fEcFnr3AZi2GAwDIQQzW+rZ3dnropTexczOcjlsxW4mnAVAq733dKbLIGLOSlHbQk8MqL0HTGpWT9ZeMia0Qi8DSWdhEdIREFUul5YG+6LSFno2mRa6Y7tOX5z6cUZLPbuOaoWei2S2y+XwBit20wwoCx10tWguookkoYDAMETamACYdgmgZto214XY0dHjkoSD0Qq9DGRMjQE4+Do0ToQGNVygJhRgTG1I+9BzYA+IzshyAZh2aerHWS317OvqS1lUmjSZLpek2nuWuwXSCl0/HQ4mNX4O0i6XYG0q/iCEYFZLPTu1Qq8u4la0XNgW5s6XUn44m7bGiPah56DXUugqD91xUVlPNwAzW+oxJew71ueGiJ5mR4d6cqmNhKH7kOofNONtqb/b7j5toQ8mmkiqKlFI773JF0IwnPrM7NaG1Bp7Aa3Qy8Ab+08AMLZeWUP0tsO8d2V8RpX/ayspm34rDz2jUnTKwoyLakaLmli0y0MXlldYuv0oM8bXURcJgWmlds59R+rv9ZEgDZGgrhbNQSxhqpRFSPvQLXeLzey2BjpORj2Ti64Vehl4/s12WhrCnNZsj0oTMOcdGZ9pa4xoKykHKQvdmeUy/dKMz8yyFbr2o2eQSJqs2NXF4tkt6Tzq1vkwdkbG59qa9NNhLqJOhR5UTzJ2QNTG3ntecbtohV5iEkmTl7Z1cOXcNgzb5TL1Iqhvyfjc/IlNHDjery+sLPpiSYSAmmBABfMWvA/OuTnjM811Kji1S2e6ZLDhwAl6ogkuO318+ukmy5AAmDSmVmcJ5SAaN9Mulxlvg/ffD7OvyfjM7LYGAHZ65OlQK/QS8/q+45zoj3PN/La0lZTlbgG4eKaa7bhyV1c5xfM8x/ti1IYCKtMgVAM3/xDGzx70uZkt9drlksXSHZ0ALJ41Pr335r5z0OcWzhjL5kPdnNDVtimklBw43k9t2HK1BMNw7s2ZmVbAtHF1BA3hmUwXrdBLzPNvthM0BFfMbYGg5UOfO1ihnzV5DHXhACt3dZZZQu+yr6uPx9bs57LTW4b97EydujiIpTuOMv+0RsY3RKBuvKpwnHrxoM9ddnoLUsKynXrv2Tyx/iAbDpzgQwunDPm5UMBg2vi6ylbodlZHtbOns5dfrz3AwhljaaoJwdkfhA//DNrmD/psKGBw4fSxrNAWOqAspC/+ZiMC+PL7zhz28zPH13O4eyDVzKvaef7NI6zc1cXl9s3wbZ+HO5dCIDTos+dOaaYuHGDpjqNlltKb9EQTfOXJLZw9eQy3LJo27OdntzZUtsuls8c7pbBukDQlT288xPvvWUo0keQL71qg/lA7Nqe7xeaiGeN48/BJT/WGcINdR3u5/ZFVvLytg8+/cx6Tm2uH/c7pli/z9b3HSyydt+mNJrj/5R188gdrWDCxiU9dfbr6QzAyKG5jEw4aLJoxjte2V7dCl1KydMdRbvzeq7SfjPIvN55JwBDDfm92awO7O3tTNRNuEizFQTt7Y5wciNNYM9gaqETiSZM9nb1sPdzDil2dPLelnQPH+5nVUs+DH1vEzJb64Q8CXDxrPAArd3fxzjNPK6XInkFKSfvJKJsOnmDTgW5efquDVbuP0RAJ8n9uWMBti2cUdJyr57fR0hDhnhe3F+SiqRT6Ygm2t/fw5qGTvLbjKM+/2c7JgQRXzWvlex+5gIZIYZf4ZaeP5ytPdnD4xACnjakpsdTewDQle7r62HjgBBsOnOCPW46ws6OXKWNrefhji7hg2tiCjnPtgjbufWkHjy7bzV9dOTi+U05KotBNKfnZyn188m2zSnH4shBPmpwcSNDdH+fkQIKTA3GO9cVpPzlA+8ko7d1R2k8OcPB4P3s6+0hYFXl14QCXzBrP/33PAt6+YALBQOEPQedOHUMkaPDrtQe4bsEEFQj0IaYp6Ysn6Y2qdevqjdPVG+NYX4yu3hhHe6IcOj7AwRP97D/WnzFtaP5pjfzddXO5ZdFU2poKVyw1oQB/+bZZ/MeTW1izp4sLp48rxX+tLAzEk3QPxDP2X1dvTO297qjafycH2H+sP6OPSEtDmOvOmMCtl0wvWBnZXDpb3QSf3HCI2y+fWdT/TzlJJE16o0l6YmrtjvXG6OqLcaw3lrp+Dxzr5+DxAfYf60ulxYYCggunj+X2y2bypxdMUa0mCmTRjHFcPa+Ve17YzocXTWNMnXuGrChF28xx0xfIWZ/8Ll95/1m8fYSKSUpJ0pTEk5K4aRJPmCRMScx6jSdN658kkTSJJU0SSZl6L540SZgm8YQkmkjSH08yEDetV/XzgPVzf+rVJGr93h9LcnIgkWoKlYtQQNDaEKG1McKEphpOb2tg7oTG1Gsq1ekU+PYf3+J//riND180jS+8e77yvY9w7RL2P3udTLVG2e/ZaxhN2P+SxOyf48nU+zHrb+p9x88Jk95ogt5YQl1E0QR90UTqIslHfTjApOZaJjbXMrm5hrkTGjlz0hgWTGwc1VNdXyzB5V97gSlja/nK+8/mrMljRvR9KWV6bRLW/rP2V+Y+y7P/TEk8of4+EE8ykDBT+ywaN+mPJRlIJK1XkwHr9wHH3uvuT6T6cOeiJmTQ1ljDhKYIE8fUMqetgTkTGpgzoZGZ4+tP2QgwTckt9y9n9Z4u/vXGs7hp4ZR02XuBa6f2V3q/xa3XpGlflzK1nvbP2XstmlBrFUtm77Wk9Xv6vd5ogp6oMhzUz4lUh8R8NNeFmNxcy6TmWiY317Jgotp7o71uNx/s5t3feYX3nTuJu941n0kFuAmzEUKskVIuPGUhKJFCP/Oc8+W4W/+bfV391IYCjG8IEwoYJE2JKSWmKTGlsuRN6yJyXjilIGgIakMBIqEAtWGDmmCA2nCAmmCASMigNhSgJhSgNhSgqTZIY02IppogTbWh1M9j6kK0NdbQXBsqmfUspeTrz2zlnhd3AKpfddAwCBiCgCGUwraURzLjAlFKphQEDUE4aBAJGkSs9bJ/rgsHqI8ErYrDAPXhIHX2zxG1juPqwoytDzGuPszYujA1ocIVxUj5zesH+Odfb6AvlmRMbYjmuhAC1aTKNNN7Lmn9HHcoZPspq9iEg/b+Su+zmqzfa0MBasIBmmpCNFr7rqkmmPq9uS5MW1OExkgQIUqz9/piCe780Vpe2tZB0BBMaKrBMMAQIrWGaq2cRkJaQZeCcMBw7D2DSCiQ+rku7Nh31h6sDwepjwRosPdefVjtu/oQY+uUHioVX3lyC0te2QnA+PoIjTVBpFS6LqX7HHsvYe25hHUt7/zPG8qj0IUQ1wPfBgLAA1LKrw71+YULF8rlK1byvxsO8cb+E3T1xkiakoAhEEJtkIAQGIZqcBMOGAQNQTBgEA6o11DAIBQQhAIGQes19bthEA4Kgsbgz4UDBkHreDUhW2kbI3J9uI2UkmU7O1m75xh7OvtSmyEpISCw1kcpeLUG1ppZaxgwhHrP8Te1vo73jPS61oQMwoFMRa0uHiO1nn6ieyDOr9ceYEdHT6ok2xDC+oe1DwUBA2svqfXJ3Et59qIhCAUNQo61DWft0WAgrcAjwUBBgTWvEE+aPLv5CBsPnOBw9wAyZXipG7tzHwUda+Lcf6nPON5z7rmg470ah4JO7TvLaAgHDN+5Hfd19fHbdQc4cHyAnmgCASm9FxBi0N4LOtbsH66fX3qFLoQIANuA64D9wCrgw1LKzfm+s3DhQrl69erRyKXRaDRVRTFcLoWYXhcB26WUO6WUMeBnwI2jOalGo9Foik8hCn0ysM/x+37rvQyEEHcIIVYLIVZ3dHQUSz6NRqPRFEghCj2XE2uQn0ZKeb+UcqGUcmFra+voJdNoNBrNiChEoe8Hpjp+nwIcLI04Go1GozlVClHoq4A5QoiZQogwcAvwRGnF0mg0Gs1IGbZSVEqZEEJ8GngGlbb4kJRyU8kl02g0Gs2IKKj0X0r5JPBkiWXRaDQazSjwV8WIRqPRaPJSktJ/IUQHsKfoB/Y2LUB19x/NRK/HYPSaZKLXI5N5UsrG0RygJN0WpZRVl7cohFg92iqvSkKvx2D0mmSi1yMTIcSoy+u1y0Wj0WgqBK3QNRqNpkLQCr143O+2AB5Dr8dg9Jpkotcjk1GvR0mCohqNRqMpP9pC12g0mgpBK3SNRqOpELRCHyFCiKlCiBeEEFuEEJuEEJ+x3h8nhHhWCPGW9TqyKb0+RwgREEK8LoT4vfV7ta9HsxDiMSHEm9ZeWVzNayKE+Jx1vWwUQvxUCFFTbeshhHhICNEuhNjoeC/vGgghviCE2C6E2CqEeGch59AKfeQkgL+XUi4ALgH+WghxBnAX8JyUcg7wnPV7NfEZYIvj92pfj28DT0sp5wPnotamKtdECDEZ+FtgoZTyLFRPqFuovvV4BLg+672ca2DplFuAM63v3GNNjxsaKaX+N4p/wG9R4/m2AhOt9yYCW92WrYxrMMXajNcAv7feq+b1aAJ2YSUdON6vyjUhPSRnHKqY8ffAO6pxPYAZwMbh9gTwBeALjs89Aywe7vjaQh8FQogZwPnACmCClPIQgPXa5qJo5eZbwD8CpuO9al6PWUAH8LDlhnpACFFPla6JlPIA8A1gL3AIOCGl/ANVuh5Z5FuDgibFZaMV+ikihGgAHgc+K6XsdlsetxBCvAdol1KucVsWDxEELgC+L6U8H+il8t0JebH8wjcCM4FJQL0Q4lZ3pfI8BU2Ky0Yr9FNACBFCKfMfSyl/Zb19RAgx0fr7RKDdLfnKzGXA+4QQu1EDxK8RQvyI6l0PUNbUfinlCuv3x1AKvlrX5O3ALillh5QyDvwKuJTqXQ8n+dbglCbFaYU+QoQQAngQ2CKl/KbjT08At1k/34byrVc8UsovSCmnSClnoII4z0spb6VK1wNASnkY2CeEmGe9dS2wmepdk73AJUKIOuv6uRYVJK7W9XCSbw2eAG4RQkSEEDOBOcDK4Q6mK0VHiBDicuAVYANpn/E/o/zovwCmoTbwTVLKLleEdAkhxFXA56WU7xFCjKeK10MIcR7wABAGdgIfRxlQVbkmQoh/AW5GZYm9DnwCaKCK1kMI8VPgKlTb4CPAl4DfkGcNhBBfBG5HrdlnpZRPDXsOrdA1Go2mMtAuF41Go6kQtELXaDSaCkErdI1Go6kQtELXaDSaCkErdI1Go6kQtELX+BohxHghxDrr32EhxAHr5x4hxD1uy6fRlBOdtqipGIQQXwZ6pJTfcFsWjcYNtIWuqUiEEFc5erN/WQjxqBDiD0KI3UKIDwgh/ksIsUEI8bTVygEhxIVCiJeEEGuEEM/YJdkajV/QCl1TLcwGbkA1ifoR8IKU8mygH7jBUurfBT4opbwQeAj4D7eE1WhOhaDbAmg0ZeIpKWVcCLEBNWDhaev9Dage1fOAs4BnVbsRAqhWrxqNb9AKXVMtRAGklKYQIi7TwSMTdR0IYJOUcrFbAmo0o0W7XDQaxVagVQixGFSLZCHEmS7LpNGMCK3QNRpAShkDPgh8TQixHliH6tmt0fgGnbao0Wg0FYK20DUajaZC0Apdo9FoKgSt0DUajaZC0Apdo9FoKgSt0DUajaZC0Apdo9FoKgSt0DUajaZC+P+fvz9J3k7a3gAAAABJRU5ErkJggg==\n",
"text/plain": [
"