{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "### Example 6.5.1" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Populating the interactive namespace from numpy and matplotlib\n" ] } ], "source": [ "%pylab inline" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [], "source": [ "x = linspace(-1.5,1.5,1000)\n", "y = linspace(-1.5,1.5,1000)" ] }, { "cell_type": "code", "execution_count": 10, "metadata": {}, "outputs": [], "source": [ "### Create mesh\n", "X, Y = meshgrid(x,y)" ] }, { "cell_type": "code", "execution_count": 11, "metadata": {}, "outputs": [], "source": [ "### Define xdot and ydot\n", "\n", "def x_dot(x, y, b):\n", " return y\n", "\n", "def y_dot(x,y,b):\n", " return -x**3+x-b*y" ] }, { "cell_type": "code", "execution_count": 12, "metadata": {}, "outputs": [], "source": [ "### Euler and Runge-Kutta integrators\n", "def intgrator_euler(xn, yn, dt, b):\n", " \n", " xnp = xn + dt*x_dot(xn, yn, b)\n", " ynp = yn + dt*y_dot(xn, yn, b)\n", " return xnp, ynp\n", "\n", "def intgrator_rk4(xn, yn, dt, b):\n", " \n", " k1x = x_dot(xn, yn, b)*dt\n", " k1y = y_dot(xn, yn, b)*dt\n", " \n", " k2x = x_dot(xn+0.5*k1x, yn+0.5*k1y, b)*dt\n", " k2y = y_dot(xn+0.5*k1x, yn+0.5*k1y, b)*dt\n", " \n", " k3x = x_dot(xn+0.5*k2x, yn+0.5*k2y, b)*dt\n", " k3y = y_dot(xn+0.5*k2x, yn+0.5*k2y, b)*dt\n", " \n", " k4x = x_dot(xn+k3x,yn+k3y, b)*dt\n", " k4y = y_dot(xn+k3x,yn+k3y, b)*dt\n", " \n", " xnp = xn + 1./6*(k1x+2*k2x+2*k3x+k4x)\n", " ynp = yn + 1./6*(k1y+2*k2y+2*k3y+k4y)\n", " \n", " return xnp, ynp" ] }, { "cell_type": "code", "execution_count": 20, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "Text(0,0.5,'$y$')" ] }, "execution_count": 20, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAV4AAAFFCAYAAACkOoVCAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzsnXdYFNf79u9RsGtMbImSaDR2TX6JxsTevoqiAgrBFhUVBdHE3iMqYiEiduzGAioaRRRFQKSpdKUoIAKCdJAOC2yZ5/2DsK/I7syQCGoyn+va64Lde59zpuyzZ2fOfR6GiCAiIiIiUnfUe9cdEBEREfmvISZeERERkTpGTLwiIiIidYyYeEVERETqGDHxioiIiNQxYuIVERERqWPExCsiIiJSx4iJV0RERKSOeS8TL8MwpxmGyWIY5oma10cwDFPAMEzYXw+Luu6jiIiIyN9F4113QA1nABwCcI5D40dEE2sStHXr1tSpU6d/0C0RERGR6oSGhr4iojZC9e9l4iUiX4ZhOr3tuJ06dUJISMjbDisiIvIfh2GYpJro38tLDQIZyDBMOMMwrgzD9FYnYhhmIcMwIQzDhGRnZ9dl/0RERERU8qEm3kcAOhLRNwAOAriuTkhEx4moPxH1b9NG8C8BERERkVrjg0y8RFRIRMV//X0bgCbDMK3fcbdEREREBPFBJl6GYT5lGIb56+8BqNiOnHfbKxERERFhvJeJl2GYiwD8AXRnGCaFYZj5DMOYMQxj9pfEEMAThmHCARwAMI1qeWHh8vJyeHh41GYTIiIi/xHe11kN03leP4SK6WZ1QkJCAoyMjLBx48a6alJERKQOICL89eO5TnkvR7zvE0SEHTt2IDExETo6Orx6lmXh5OSE3NzcOuidiIjIP2Ht2rW4f/9+nbcrJl4ekpKSkJWVhaCgIDRs2JBT6+bmhl69esHV1RWffPKJoPgsy76NboqIiPxFQUEBEhISBGm3bNmCAwcO4ObNm7Xcq6qIiZcDmUyGJUuW4PDhw+jcuTOnVi6X49q1a2jevDl27tzJqU1KSsKSJUswePBgODk5vc0u1zoSiQSlpaWCtBkZGRBy6Z1lWcjl8n/aNZF/OeXl5YJ0paWl2LJlC5YvX45Xr15xaps0aQIHBwfcuHEDZ86ceQu9FAgR/Wce/fr1o5qwdu1acnZ25tWVlJSQkZERXb9+nUpLSzm1SUlJZGJiQp9//jnZ29tzal+8eEHDhw+nzp07k7GxMWVnZ6vVpqen065duyg0NJQUCgUREfn7+9O+ffvo+vXrFBYWRnl5eUREFBkZSbm5uVXev2jRIjIyMqKdO3eSl5cXFRcXK+O+ztOnT8nIyIgmT55cpb2SkhKSyWRVtHZ2dqSrq0tz5swhe3t7yszMJCKioqKiKrq0tDSaMmUKGRoa0u+//05BQUEkk8movLy8mvbYsWM0ffp02rlzJ7m5uSn3SWJiIrEsq9QpFAratm0bnT9/nvz8/Cg5OZnkcjlJJBIKCQmpEvPVq1cUFxdH6enpVFhYSHK5nIiIPD09KSwsrEpcVeTl5ZGFhQW5u7tTSUmJWl1paSkdOnSITE1N6cCBAxQTE6NSx7IsxcTE0LBhw0hbW5u2bdum9rwqLCwkS0tL6tChAw0YMIB8fHzUxoyOjqbevXtTixYtaN68ecrz4U0kEgkdO3aMWrZsSV26dKErV66o3Qfp6ek0bdo06tevH82ePZtSUlI4258xYwZZWVnR7du3SaFQUHFxMSUnJ9OrV6+opKREee4+fPiQfH19lf8TEVlYWNCkSZNo6dKl5OTkpDyH3zyXKwkKCiI9PT3avn0753Gp7N+aNWvo999/59SpA0AI1SAXMfQfqjLcv39/EmoZvnPnDtzc3LB3715OXU5ODubMmYP169dj8ODBanWvXr3Czp07kZGRgU2bNuGTTz5B27Ztla8TEeLj4+Hj4wM/Pz/k5+ejY8eOkMvl6NOnDxYuXIj69esDAIyMjKqNJFmWhbOzM4YPHw5zc3NMmTIFeXl5ePLkCRITE5WP/Px8vHjxAs+fP4exsTHWrFmDL774AkDFiOLx48cICAhAcHAwJBIJgoKCoKuri61bt1bpr0wmQ2BgIDw8PBAWFgYiwsuXL3Hs2DH88MMP1bb97t278PDwQHZ2NlJTUzFq1Chs2bIFTZs2VeqkUimCg4Ph7e2NR48eAQCCgoKwe/duTJ06VXkTRCKRIDIyEqGhoXj06BFycnIQExOj/LUxevRosCyLiIgI5Xa/ePECaWlpkMvlcHZ2xpAhQ7B69WpMnDgRDx8+hLOzM0pKSlBSUgKJRAKFQoGYmBhERUVh6tSpsLGxgZaWFszMzJCdnV3lhoympiacnZ0xePBgzJkzB9OnT8eLFy9w5swZJCYmKn8hNGzYEBKJBNnZ2Vi+fDn09fWhoaGBFy9eICgoCKGhoYiLiwMAfPnll4iKioKuri7mz5+PBg0aAACKi4vx4MEDeHl5ITo6Gs2bN8d3332H8PBwWFtb49NPP1X2Ky8vD56ennB3d0dGRgZ69OiB+vXr4/vvv8fkyZOV28CyLMLCwuDh4YGAgABoampiyJAhCA8Px+bNm5Xnh1QqRVhYGAICAhAUFITi4mK0a9cOLMuiXbt2WLt2LZo3bw4A2LNnD6Kjo5GbmwuGYcAwDNq1awcXFxfMnDkT5ubm0NLSQnh4OP744w+UlpairKwMZWVlUCgUSEpKQkhICLS1tbF//350795d+TlJTEyEl5cXfH19kZeXh2fPnmHs2LGwtLREy5Ytq5x7RARXV1ccPHgQBgYG6Nq1K4YNG6b2hpqNjQ2ysrKwa9cu1Ksn/IIAwzChRNRf8BtqkqU/9IfQEW9aWhqNGzeOysrKOHVJSUk0duxYevLkiVpNcXExbdu2jfT09CggIECt7uTJk7RixQpydnau8g3ON9qqpLCwUDmi5MPNzY3c3NwoJSWFM35BQQGtWbOGrKys6MKFC8qRoCquXr1KP//8My1YsIAeP36sVieVSsnMzIxmzJhBW7ZsIYlEolbr7+9PpqamtHLlSrp3755aHcuyZGdnR/v27aNLly6pHQEREb18+ZIuXLhA4eHhvMc3JCRE0D6VSqXVfhm8evWKnj59Wm2kpWrkdfr0abp48SLFxsZWGeG9/jdRxah+xowZZGdnR1FRUcpjp+oYHjx4kGbNmkVnzpyhtLQ0tX0vLS0lPT09srCwID8/P5JKpWq1N27coO3bt5OnpycVFhYqn1fVflRUFOXk5FT7FcI38qzkyZMn1fapKsrKymjp0qW0ZMkSsrW1VXtM5XI52draEgBauXIl53l/5swZmj9/PsXFxQnqK1HNR7zvPBnW5UNI4pXL5aSvr0+xsbGcusjISBo7diy9fPlSrSY1NZXGjRtHrq6ughOoiIhI7ZCVlUXXrl2jdevW0e3btzm1Bw4coIYNG9L9+/cFxa5p4hUvNbyBlZUVOnfujBkzZqjV+Pn5wdraGufOneOcvVC5k2vyk0VEROTdk5qaChsbGwQGBsLX1xcaGtyWh5peahAzwmv4+fkhOTlZZdKVyWTIzMwEAERHR+Py5cu8U8YYhhGTrojIB0iHDh2wd+9eXL9+HTk5b381AjEr/EVOTg6srKxga2ur8rXx48cr/1+4cCGaNGlSl90TERF5B7Rt2xbt2rV763HFxPsXMTExsLGxqXKXHQCKioowfvx4lJaWCj4AiYmJtdBDERGRt01oaCi8vLzqvF0x8f7F4MGD0bdv32rPFxUVoW3btlizZg1vjNjYWEydOvWDM0WIiPzbUCgUgnS9e/eGnZ1dnSff93KRnPcFIsLy5ctha2uLbt26cWpZlsW8efOQmpqKc+e4SsVVJPMHDx7A398f8+fPV86TFBEReTssWLAA3333HczMzDhvjDVq1Ajnz5/HrFmzAAAjR46sk/6JI14Ozp49i0GDBvEmXaBiNsSkSZPg7OzMuaaDQqGAmZkZdHR00KNHD86km5mZibVr12LTpk24ePEiuGagPH36FHl5eVWeIyJERESguLi4yvOFhYXV3u/i4gJra2sEBARAJpMp+6oKY2NjWFhYICgoSLnWhKq+PXjwAEZGRjh+/DjS0tKq9OtNNm7cCCsrKwQFBSnbVbWORWpqKvbt2wd/f3+UlZUpn6/s8+skJiYiLS2tSnuqtl0VCoUCycnJgrQxMTF4+fKloJiOjo4ICgpCQUEBpzY8PBz79u2Dvb09Z51AIsKhQ4ewbt06bN++nfNGUHJyMszMzLBw4UI4ODhwnk/29vZYsmQJVq5cyXnprLCwENu3b8e+fftw+fJlzpFmYGAgPD09qx0TVahaa+HevXswNTXF1atXkZ+fD+D/zxx6k1OnTqFly5aYNGkS7yI4lcm3Lke+4nQyNaSkpMDc3BzXr1/nnJlARNiyZQuaN2+OVatWccZ8+PAhrKysYGBggI8//hhTpkyppiksLIS3tzfc3d3x8uVLBAUFYf78+di4caPyhp6joyO8vb2VTh+ZTIa4uDhERERAV1cXtra26Ny5M6RSKXbu3ImEhARl8m3atCnu37+PCRMmYPny5co1KMrLyxEaGgpfX1+EhoZCoVBAS0sLSUlJsLW1RZcuXapsc1hYGFxcXBASEoLWrVuDiPC///0P06dPr+IKkkgk8PT0xM2bN5GRkYH+/fsjKipK6QSrRKFQIDQ0FO7u7ggJCUHTpk2hqamJ7t27Y+XKlUrnVllZGQICAhAYGIiwsDCUlZVBS0sLAQEBWLt2LQwMDJTt37x5E25ubso1I+rVq4ecnBzI5XKsWrUKEydORL169SCRSDBnzhzltgFAvXr14Orqim+//RZz5szBnDlzoKGhAR8fH9y+fRvFxcUoKipCaWkpUlNT4e/vj6FDh+LQoUP4+uuvAVR8GT5//hxxcXGIj49HZmYm/P39Ua9ePSxbtgyLFy9WHlOZTIaIiAilKywvLw8eHh5YunQp1q9fj48++ki5r1iWRXh4ONzc3BAYGIjExER89NFHOHr0KHr06FHlfCovL4ePjw9u3ryJxMREREZGYtWqVTA3N692XiclJeHWrVvw8vJCVlYWSkpKcPr0aeX2VJKVlQU/Pz/4+voiMTERAQEBGDt2LKytrdG+fXsAFet0REREICEhAQkJCUhKSkJKSgoCAgJgYGAACwsL9OnTBwAwe/ZsSCSSKm14e3ujffv2mD17NpYsWYJGjRoBANLS0uDu7g5PT08UFBTg22+/xaNHj3Dw4EGoqiBeWFgIS0tL5ObmYvv27crzWhVlZWWYNWsWzM3NazzyFZ1r/9BAQVThxDEwMOB1rrAsS7/99hvZ2tpy6lJTU2nu3Lm0YsUKtf745ORk0tPTo5kzZ9KJEycoMTFR+d43yc3NpaysLCoqKlKuj+Dl5VXN266K7Oxsmj9/PllZWZGTk5Nap5JMJiNbW1saOHAg6erqUnBwsNqYaWlp9MMPP1DHjh3p119/rba+QiUKhYJOnz5NHTp0oJ49e3JOYi8oKKApU6ZQnz59aPbs2ZSVlaVSx7IsBQQE0MiRI+mnn34ie3t7tWYVuVxOe/fupRUrVtCJEycoKSlJbfsSiYQsLS3Jw8NDuW4FUcUk/Li4OMrMzKSSkhJiWZYePXpEXl5e1Zx9v//+O507d44ePHhAmZmZxLIs+fr6Unl5eRWdj48PGRoakqWlJbm7u1NeXh6VlJQoz4HXWb9+Penq6pKFhQU9ePCAZDIZpaenV9vmSpfblClT6MCBA/TixQsiqnDUvcmJEydo4sSJZG5uTrdv3yaJREK5ubnV1t4oKCigSZMmkYmJCZ09e5ZevHhBLMtSfHx8tZj37t2jo0ePkru7O8XFxZFUKqXY2FhKSEhQu88rUSgU9Pvvv5Obmxun000ul9OxY8eoe/fu9OOPP5Knp6da7dOnT0lbW5s6d+6sdj0Jogonn6GhId27d6/KcecDooFCPUJHvCdOnIBMJoO5ublaDRFh48aNaNeuHZYuXapWd+DAAXh7e8PKygq9evVSq6v8WS3O+xXhQiqVKkf/XJSWlkIqlVYZKasjLy8PLVu2fCcLgtclN2/ehJOTE7Kzs3H27Fm18/DLy8thZGSE8vJyuLq6Ctov4oj3H454ExMTSV9fn3PkWLmS0cGDB3njxcbGinZhEZH3DK41KViWpaNHj1Lbtm3J0dFRUDzUcMQrDq9eg2VZLFu2DPv27as28nx9vdg7d+6gU6dOWLJkCW/Mrl27/utHEiIiHxqamppqX2MYBqampoiPj4dCoQDVwlUB8VLDa9jZ2UFTUxMLFiyo8nxRURGsrKxgbW1d210UERH5ABHXavibVK5Fa2JiUuX5zMxMjBw5skYVEoRO3hYREXm3EBGioqLqvF0x8f5Fw4YNceLEiWqXBVJTU5Gbm4vRo0fzxigtLcXmzZsRFBRUW90UERF5ixBVFLOta+eamHh52L9/P27duoUxY8Zw6jIzM9G/f3/cvXsXP/74I29ciUQCb29vsdiliEgtcP36dVy/fp1XV69ePZw8eRL79u3D48eP66Bnf7VbZy19gHh6eqJ9+/bo2bMn58V4ALh//z66du2K3377jfNmmlQqxdy5c9GmTRvk5+dzTh978eIFIiMj8eTJE6Snp6vVsSz7j28AyOVyeHt7Cyoo6OHhgZs3b1ZxjqmCiGBtbY2IiAje/kVFRcHd3V1Q+8+fP4dUKuXV/ROE7k8iEnxpSSKRCP6iTUlJQV5eHu8lroKCAoSHhyMsLAyRkZGc/Y6Ojoafnx88PDzw4sULtbrS0lIEBwfjwYMH8PHxUekKrCQjIwNJSUnIzMys5pB8HZZlBRdJVbcNRIQzZ85UcUGqQ0dHB5GRkZg9ezays7M5tY0aNcK5c+ewfv16wdWJ/yniWg1qKC8vx++//45r167xar28vHD16lVcvXqVM5GWlpbC0tISKSkpsLS0hL6+vtq2fX19sXv3bnh4eGDevHnYs2eP8vWCggI8fPgQpaWlKC0tRUlJCXbv3o3+/ftDW1sbM2bMUM713LdvH3x9fdG0aVN069YNDRo0QEJCAtasWVPFjSaVSvHo0SMcPnwYRISBAwdCJpNhyJAhGDJkSJX+ff3113BycsLMmTPRtGlT6Onp4eXLl1i0aJHSYVTJ2LFjcfHiRaxfvx7fffcdevXqhU8//bSaM+iTTz6Bp6cnjh07hkaNGkFbWxtpaWkwMTFB69atq2jv3r0LT09PKBQK9O3bF61btwbDMDAzM6vyBRkTE4ONGzeiXr166NixI7p374579+7ByMgIenp6VY7VmTNnEBsbi+bNm6N58+ZISEhAfHw8fv75Z0ycOBGNGzdW7qfQ0FAUFhYqH7a2tujYsSNGjx4NMzOzKivcVdYPe/bsGVxdXXH8+HH07t0bO3fuxNixY6sd+/T0dPj5+WHr1q2IiorCrFmzYGNjU6XeXSXPnj3D5cuXYWFhgZ49e+LMmTMqv/RlMhm8vLywY8cO+Pn5YePGjRgxYkQ1XWX7V69exZo1a9CuXTucO3dO5aCjrKwMDx8+xOHDh3Ht2jWMGDECR44cqeKcIyJkZ2fj+fPniI2Nxbp169C4cWNMnz4dmzZtUrr2bt26haCgIDRt2hRNmjRBVlYW7t+/j+nTp2PKlClo1aoVgIrk3apVK2zZsgUZGRkYMGAA2rZti27dulXbngYNGmDTpk148uQJjI2NMWvWLOjq6qJx48Yq99FHH32EP/74A8bGxjh//rzK/f02EWc1qGHHjh3o3bs39PT0OHWPHj2CpaUlLl26VC3pvI6fnx+srKywYsUKjBo1ChoaGtVOgCtXrsDJyQnl5eUYPnw4vv32W5SWllb7gGZlZeHmzZto3LgxGjdujEaNGuHUqVMYPXo0DAwMVJ40xcXFiI2NxdmzZ+Hu7o4ePXpg3bp11QpTAhXJxd/fHz///DNSU1NhZmaGXbt2oUWLFtW0r169wtGjR7Fp0yb07t0b9vb2+L//+79qOiLCo0ePsGjRIoSEhGDVqlXYtm2bynUtCgoKcPnyZZiamqJ169aws7ODoaFhNZ1CocCTJ0+wfv16eHt7Y9CgQbCzs6u2toZCocDLly8RHh6OFStWQEtLC7q6uli6dKkyqeTk5CAnJwdFRUUoKiqCh4cHwsPDMXz4cEyYMEFpfikpKcHRo0fRokULtGjRAs2bN8eVK1fQv39/GBkZoU2bNsp27e3t4ezsrEz6n332Ge7fv49ly5ZVKUpJRDA1NUVmZiY+/fRTDBs2DCzLokuXLhg0aFCVbYmNjcWJEycQHR2N7t27Y9KkSUhJSYGhoWG18+/evXtwdHRERkYGRo4cicGDB0MqlVYrypqTk4Pz58/Dy8sLbdq0wZQpUyCVSjFq1Khqx9ze3h7Xrl1D/fr1MXjwYHz77bdITk7GzJkzq5zP2dnZMDU1Rdu2bdG1a1d069YNAQEBGDVqFEaNGlVF++rVK2RkZEAikUAikSA0NBQeHh4YNWoUxo4dq/J8YlkWwcHBMDExwZMnT7B48WLs2rULzZo1q6aVy+WwtbXFlStXYGhoiLVr11bTVPL8+XMsXboUjo6OaNasmeCpoKKB4h8aKIiIEhISyMjIiNf4EBsbS+PGjaOCggK1moKCAlq8eDEtXbpUrZW2ksDAQM5YbwOhZg6WZam0tJTXgkxUUdAzJSWFEhISeC2h6enplJCQQFFRUSrtq5WUlZVRbm4u5eTkVCua+CZcRThf5302sqizkr9JQkICRUdHC9qWK1eucNYErCQlJYXc3d05TQWVJCYmVrM8vytYlqWcnBzKzs6mrKwsTntxSUkJmZqaUuvWrenatWuccYODg0lbW5uOHz8uuC8QLcPqETrinTp1Knbt2oUvv/xSrSYtLQ1z587l/Vni7+8PIqo2chEREal7FAoFIiMj8c0336gdzRYUFGDy5MlITExEbGwsb701oOYjXjHxvsGNGzcQGRmJjRs3qtWUl5dDX18fR44cUbkikoiIyIcNy7L4448/oKWlBW1tbV69mHg54Eu8JSUlmDJlCm7cuFHl2mNxcTEePnxY5Vprfn4+WrZsWav9FRERebcQUa0skiNOJ3uNXbt2YfXq1VWSbmXNtTenTolJV0Tk3wHXdMDaWmdFTLyvsXz5cvzvf/+r8tzRo0dx//59fP/997zvJyJcuHBB0FxUERGR9wNbW1vcvHmzTtsUE+9rqFqfMyYmBnfu3MFnn33G+d6UlBSMHz8e9+7d4yz98zpCJ5SLiIjUDKlUiocPHwrSLlu2DKdOnUJ0dHQt9+r/IyZeDsLCwtC4cWNBF9f9/Pzw8OFDzsXTgYo5k8eOHcOYMWPq1KIoBKELAcnlciQlJQnSJiYmVivrog6uemEib5cP9d7Oq1evBLn/WJbFpUuX8Ntvv/Ge15qamjh58iSWLVumrOVW69Rk7lldPQCcBpAF4Ima1xkABwDEAYgA8J2QuELn8VZiYGBAGRkZvLqUlBQaN24cJScn82ovX75MGhoatHv3bk5dUVERjRs3jvr27UurV6+m0tLSahqZTEYFBQV0584dsrCwIG9vbyorK6uiiYqKotWrV9OZM2coODiYHBwc6OnTpyrb3LVrF+nq6pK1tTVFRUXR9evXVc7jzc3NJVNTU9LX16eTJ09SVFQUxcbGqozp4eFBkydPpgULFpCfn5+yXM2byOVymj9/Punr65OdnR09f/6c/P39VcZ0dXUlPT092rp1Kz148ICcnZ3VzoPdvHkz2drakpubGwUGBpK9vb3K9l+9ekXe3t70+PFjSkhIIFtbW/Lw8OCcI1xaWkr+/v60YsUKcnJyUlueqHL7Nm/eTIsWLSJra2tKS0tTq01ISKAxY8aQjo4OHT16lHMutY2NDfXt25dGjhyp9rgSEeXn59PgwYOpffv2tG7dOs65uJcuXaLu3bvTN998Q76+vmp15eXlNH36dBoyZAj98ssvnHPQ/fz8aNasWbR582a6e/dulddKSkrI39+fIiMjKSEhgY4cOUJXr15Vec5fvHiRJkyYQKtWraKQkBC6e/cu5zG6evUqTZw4kZKSknjnpD9+/JgMDAwEzwt/HdRwHu87T7IqOwUMA/AdR+LVAeD6VwL+EUCgkLg1Sbzu7u60ZcsWXp1EIiEdHR1lTSt1SKVSWrduHf36668UGhqqcgJ8WVkZ3bx5k4yNjUlfX58mTpxIx48fV6lduXIlTZs2jebNm0dTpkyhtm3b0uLFi6vVRpPL5fT8+XNycnIiKysr0tLSIoZhyNzcXGWykslk5OPjQ4sXLyYANGLECLW150pLS+nKlSv0ww8/ULNmzejixYtqtz8xMZEsLS2pSZMmNGLECEpPT1epKysroxs3btCYMWOofv36tHXrVrWJOjg4mHbu3EnNmzendu3aqay3lpOTQ76+vnTkyBHS1dUlANS3b1968OBBFV1ycjIdOHCAtm3bRqtWrSItLS3q0KEDrVq1qlpC2bFjBxkYGNCMGTNo9uzZ9Mknn5C5uTkFBARU0WVnZ9OsWbNo8uTJZGhoSMOGDSM9PT0KDQ2ttj0RERG0efNm0tXVpRkzZtCAAQPo4sWLKo/9s2fPyNLSkiZMmEALFiygn376iXJzc6vpKuvRLVq0iPT09MjQ0JCuXLmicr/n5ubSsWPHSFdXlxYtWkQTJkxQWe8vPz+fHBwcaNasWaSvr0/jxo0jS0tLlYncxsaGJk+eTJMnT6a5c+dSp06d6MiRI9UGBzk5OWRjY0NbtmyhVatWUc+ePenTTz8lU1NTyszMVNnfJ0+e0MaNG0lTU5OGDRvGaRRJTEykiRMnkqGhocpteh1HR0dau3Ytp0YV/4rEW7Ed6MSReI8BmP7a/88AfMYXU2jiVSgUNH78eF6nGcuyZGxsTN7e3py6uLg40tHRIWdnZ06dt7c3HT58WJmUVCUcVeTl5QnSsixLUVFRVFhYyKtNT08nHx8f8vHxoeDgYE6n1IMHD+jSpUt07NgxzgKS+fn55OTkRKdOnSIHBwfOEUhERAS5uLjQlStXKCoqSq2urKyMAgICKCwsjGJiYqp9qF8nJSWF0tLSBLnxkpOTBbnDioqK1LbsS9LoAAAgAElEQVQpk8mq7Guu4on+/v4UEhKi7Ju6mFKplFavXk3+/v7EsizJ5XK1/XR2dqbNmzcri1FybfeKFSvo2rVrynbVaWNiYujChQvKRM/ldisuLlb2raysTJAzjqiiwKuQfZ+Tk0Ourq509epVcnFx4XzP48ePqW3btjRs2DDez8qGDRvo119/VZv0VfFfSbwuAIa89r8ngP5qtAsBhAAI+eKLLwTtxPPnz9PRo0d5ddbW1nTkyBFOTXl5Oc2ePZuzsqmIiEjtI5PJyNXVlbNiNlGFZVhDQ4N27twpOHZNE++HenNN1eQ6lXcLiOg4EfUnov6vL2CijrKyMtjb22P+/Pmculu3buHly5cwMzPj1DVo0ABnz55Fhw4deNsWERGpPTQ0NDBu3Dj078/tc+jfvz/8/f3h5+dXa9VkPtTEmwLg89f+1wLAv0inAA4fPgwzM7Mq/uz8/Pxqd9wbN26MvXv3vo0mRURE3jP69++PP//8k3fN6b/Lh5p4bwCYzVTwI4ACIlK/UrhA8vLy4OPjU2UpyNLSUujp6SnXY61k1KhRvIuji4iIfLg0bty4ytrKb5P3MvEyDHMRgD+A7gzDpDAMM59hGDOGYSp/198GkICK6WQnAHBPnhVISUkJtm/frrQJyuVyzJw5E0lJScpFm/l4/vz52+iKiIhIHZGcnIwzZ87UaZvvZQUKIprO8zoBWPy229XS0oKWlpbyfw0NDTRr1gxr1qzhfa9UKsWqVavQs2dPdO3a9W13TUREpIaUlpZW+6WqCi0tLbi5ueH7779H796966Bn7+mI933h2bNnaNy4Ma8bDQDWr1+vtlLCm7x48QLbtm2rO5fMO4T+/+wSEZE6Zfv27Th8+DCvjmEYHDx4ECtXrqy1a7pvIiZeDnbv3o3Vq1fz6srKyvD06VNcunQJfDMnjhw5gi5duqB169acK5zFxsbCzs4OVlZWuHPnDmdMb29vzmKYfn5+SE5OBsuyCAwMVKsLCQnBihUrEBgYiOzsbE4L74YNG7B69WrEx8cjKSlJbXJNTU3FxIkTcfr0aZSWliI5OVltzMOHD2PNmjWIiYlBYmKi2jvKUqkUa9asgZeXF2QyGZ48eaI2Znx8PB4/fgypVIqwsDBB62Pk5+cjMDBQ0BfGw4cPERcXx6vNz8/HjRs3EBwcjJSUFE7ba1BQEBwcHGBvb4+nT5+q1clkMpw5cwbW1tawsbHhTBrR0dHYvXs3li5dCnd3d7U6IoKjoyO2bduGDRs2IC8vT6325cuXOHPmDI4cOQIXFxe1OqDiHPXz88OLFy84F5GSSqXw8vJSe+zj4uIwdepUXL9+HWVlZZzW9W3btqG0tBRLlizhLNgJAK1bt8bKlSuxYcMGTt1boyZzzz70R02ca4mJiTRnzhxB2mXLlpGbmxuvLigoiMaOHUurV69WO0FdLpeTm5sbTZs2jRo3bkxbtmxRaWG8ffs2GRgYkIGBAXXu3JkYhiEzMzPKycmppnVwcKAlS5aQnp4eaWpqko6ODj169Ehl+9HR0bR161YaNGgQtW/fnjw8PNRuT1hYGJmZmVGnTp1o7ty5aif9l5eX04ULF0hHR4fatWtH7u7unDGXLl1K3bp1o//9739qbbgxMTG0e/dumjBhgnI/qXJPPXnyhLZu3UpGRkbUtWtX6tChAx07dkzlZH4zMzMyMDAgfX19AkA9evSga9euVZuY/+jRI5o2bRoZGBhQ7969CQAZGRmptIz7+fnRjh07lMezT58+5OzsrHKy/4sXL+jAgQM0YsQIatiwIf3+++9q7auJiYm0ZcsW+uKLL2jgwIFq54mzLEv379+nadOmUfPmzenw4cMqdUQVBpcDBw5Q3759qWfPnmpt4JUOy1mzZlHDhg3JxMREpdkoJSWFDh8+TAsXLqRevXpR48aNac2aNZSdnV1N+9tvvynP5wYNGlCHDh3o5MmTKj8neXl5ZGdnR9ra2tSqVSvez97NmzdJX1+fDhw4wKkjqnCE3r59m1f3Jvi3GChq41GTxLt48WKKjIzk1bm4uNDq1at5dQ4ODmRkZEQFBQVqHTZSqZT09fXJ1taWMjIyKDw8XG2812O4ubkJMmgkJyfT+fPnydPTk549e8bp9HFxcSFra2uytLRUOp9UUVxcTDt27FCuQcDlCnrw4AGtXr2aZs2aRUFBQWp1MpmMDh48SOvWrSNLS0vOWloRERF08OBBOnjwIPn4+KjVsSxLly9fphs3blBISAinKzE1NZWuX7+uNulLpVJl4vbz8+OsM+fm5kZBQUEkkUjI3d2dM5EuXryY3NzcKCMjg6Kjo9XG3LdvHy1cuJD8/f0pOjpa7boLJSUlNGHCBNq8eTOlpaVRTEyM2piOjo5kaGhIN27coKSkJLXrLrAsS/PmzSM7OzvKyMjgNCMkJiaSu7s7ZWZm0qNHj9TaxF+noKCALl++LEjr5eVFFhYWtGTJEoqIiODUbty4kQDwJtWysjLS1tbmHByoQky8byHxpqWl0bRp0wTpxo8fz1v8z9bWljZs2CDIqioi8jYRajv/t5+bLMuSn58frVixgtO6TUS0detWAiBo0atKapp4xdI/Kli1ahWmTZvG6XBhWRaGhoawtrbmncVQVFSE5s2b17i/IiIidY9EIoG5uTm+/vprrFixQtB7xNI//5CcnBy8ePGiStJVdUMmJiYGhoaGgqaOiUlXROTDoUmTJvjjjz8wdOjQWmtDTLxvsH//fixdulT5f0lJCbZs2VJN16tXL8yYMaMOeyYiIlJXMAwjqNzX30VMvG9Qr149DBs2DEDF9e+5c+ciMzNT8PvFemsiIh8e9+/fr9P2xMT7Bq+PbrOyshAdHY3u3bvzvo+IsGPHDtEyLCLyAfLnn3/iwYMHddaemHg5KCwsxNChQ3lNFHK5HKampjh06JAgy6FcLq/zb1gRkf8agYGBOHnypCCtpaUltm3bxmu0eFuIiZeDI0eOYNGiRVWWiFQFESE+Ph7z5s1TLrCjCoVCAXt7e/Tu3VvQ5YuSkhJERUUhKytLraa8vJy3+F+lCyg3N5dXe/PmTZSUlPDamf39/XHu3DnIZDIUFRVxanfv3o3w8HCUlJRwtp+cnIzz589DKpWisLCQM+aTJ09QUFCA/Px8wdsvtJhmQUGBIJ1cLq+19VqJCGVlZbyJQC6XIzExkdO9VxkvPT0dnp6enOcTUOGIe/bsGfz8/MA364llWeTk5CA3N5dTJ5SCggLeNo8fP46XL1+isLCQUztgwAA8e/YMBw4c4G23RYsWmD9/fp0t9SomXjVIJBLEx8ejb9++vNpLly5BT08PlpaWnLry8nLY2dlBS0sLU6ZMUauTyWSYMmUKmjVrBjs7O5Vl5+VyOY4dO4YdO3agTZs20NPTw6VLl1SeiC4uLtDX18fkyZPx/fffw9vbW2W7RISSkhLMmjULQ4YMwcaNG9Umln79+oFlWUyaNAlDhgxBeHi42u2ZPn06Tp8+jfHjx8PY2FhtMmnbti1YloW+vj5+/PFHuLm5qY2Znp4OU1NTTJgwAePHj0dGRoZKnUQigZGRESZPnoyvv/4aO3bsUGsbvnr1KoyNjTFx4kR8//33OHLkiNovICcnJxw6dAjt2rXD6NGjYWdnp7KaLREhISEBJiYm+PTTTzFy5EiEhYWp3a60tDT06NEDDRs2hIWFBWdicXJyQqNGjTB+/HjO1fMUCgX09fXRvn17BAQEcNraw8PD8dFHH2HUqFH4+OOP1Q4kZDIZ5syZA01NTRgbG6vV5eXl4cqVK2jTpg26deuGAwcOqDynfH19sWDBAkyfPh29evXCzp071drL+/fvj82bN2PMmDFYtGiR2irCDMPg999/R05ODnbs2AEfHx+12w0AhoaGCA4ORmJiIqfubSDO41XDqVOn0KxZM0ydOpVTl5OTg9mzZ+PGjRuoX7++Wl1hYSFmzZqF1atXo2vXrmjXrp1KXaVHf9CgQWjZsiV++eUXlTqWZeHp6Yn69evDwcEB06dPx8iRIzn7cP78eTx+/BitWrXCtGnT0KVLF5W60tJSbNiwAQkJCfjuu++wYcMGtWsPBwYGYvv27cjIyMDevXsxePBglToigpWVFVxdXdGlSxccP35c7cpRcXFx+O2335CSkoKVK1di8uTJarfpwIEDCA4ORrNmzbB161a0bdtWpS45ORk7duxAq1at8O2332LKlCnVkgURobi4GH/88Qfy8/MxcOBADBgwAB999FG1eDdv3kTTpk1hb28PQ0NDaGtrq9z3fn5+uHbtGjp06IDs7Gxs2LBBZTwAWLlyJfLy8tCyZUtoa2tDW1tbpS4pKQkrV65E+/btoampid9++w0ff/yxSq2TkxOOHz+OYcOG4eOPP1ZbMaWwsBBr165FUVEROnToAHNzc3Ts2FGl1tHREfb29vjhhx/AMAzWr1+PevWqj+GICEuWLEHv3r0RGhoKU1NTDBgwQGVMhUKBwsJCXL9+HfHx8Rg6dCgGDhyIFi1aqNVbWFggICAAPXr0wO7du9V++UilUnz//feQSCR4+vQpGjRooFIHVCxgtXbtWtjZ2aF169ZqdW9S03m879xNVpcPoc41lmVJR0eH15FGRLRgwQK16x5U8urVK9LR0VFZXfZNMjMzBRX6ex8R6n5iWVawo+q/hNDjXlJSorKqsCqePHkiqMikXC7ntBS/jqqy6+87qamptHjxYs61KogqnH6jR4+mYcOG1Sg+auhcey/X433XBAQEoF+/fpzfjEDFiksfffQRvv32W06dQqGAjY0Nevbsydu2uhHbh4CqUY8qGIbhvW7+X4Tr/sDrNGnSRPDC/ELXl61fv76g2TsA0KhRI0G694n27dvj0KFDvNfvNTQ0sH79evz000/Iy8tT+0vinyJe41XB8ePHYWpqyqkhIhw+fFilueJN2rZtKyjpioiI1C7qLvO8zujRo+Ht7Y1nz57VWj/EYccbZGVlQSqVKqsC5+bmIjg4uNr1NoZhcOnSJc5rqiIiIh8mX3/9da3GF0e8b3Dx4kWYmJgAqBjVmpiYqJ0qIyZdEZF/BzKZrNamBqpCTLxvsGTJEowYMQIA4OXlBWdnZzRr1kzQe+3t7WuxZyIiIrWFo6MjLly4UGftiYn3DerXr6+8ydGmTRv88ssvnMtDVrJnzx44OjoKbofPdCAiIvLPUCgUcHV1FaSdNm0aHBwcIJVKa7lXFYiJl4MLFy5g7ty5+Oyzzzh1L1++hLW1Nfr168cbMyMjAzNnzkR8fPzb6qZghPyUKikpAREJsk5GRkaCZVneuElJSXj16pXaie6vk5KSAgC87bMsCyIS/EGpy5+R7wISOB+/vLxcsFbI8XoXZGZmCnL11a9fHx4eHrh06RJvTA0NDcydO1ewxfifIiZeNbAsi4iICEEX2VNSUjB37lwsWLCAU5ednY0BAwYgJycH//d//6dWV1xcDEtLS+jr68PKykqtJTYvLw8HDx7Ehg0bcPv2bbUuq6KiIhgbG0NHRwenT5/mtNi6urpCV1cX2tranF8OCoUCDg4O+Omnn7BgwQLOmCUlJZg3bx50dHTg5+enVgcAp0+fhpGREQwNDTmLNz579gwTJkyAtrY2bt26xRnT1tYWJiYmmD17NqdVu6ysDA4ODtDT08Mff/yBkpIStVqZTIYrV65g9erVcHR05Cy6yLIsrK2tYWFhgb1796K4uFitNi4uDhs2bICxsTGuXbvGuV3Xrl3D3LlzMXnyZM5ip8XFxbCwsMDw4cNx4MABzmlrfn5+MDExwYgRIxAXF6dWV1ZWhl27dmHOnDnYtGkTZxK8d+8eNmzYgO3btyM0NFStzsnJCTNmzMDBgwc57d0JCQnQ19fH+PHj8ejRI7U6ALCxsYGLiwtu3LjBuY8A4KeffoKLiwskEgmn7q1Qk0m/H/qjJjXX/Pz8yMrKilfHsizp6uqqLDL5JseOHSNzc3Py9fVVqykoKKDNmzfTl19+SYsXL+Y0JZw9e5bGjh1L3333He3bt48yMzPValNTU2nEiBE0bNgw2rZtG6c55OHDh9SvXz/66quvKDAwUK2OZVnasmULtW7dmszNzTkNABkZGTRw4ED65JNPKCQkRK2OiMjOzo6+/PJL0tbW5qy3lpWVRaNHj6Zu3brRvn371LYvk8lo06ZNNGjQIDI0NKRnz56p1JWWltKFCxeoW7duNGXKFDpy5Iha84GjoyPp6elR165daffu3SoLOBJV7KOpU6dSjx49SEdHh7M+2+PHj8nQ0JDatm1LR44cUasjqih2OmbMGPrqq6/o5cuXanXl5eW0Z88e+uqrr2jFihWcx+jp06c0efJk+uKLLygsLEytrqCggNauXUtdu3almTNnqi1ySkR0+vRp0tbWps8++4wuXrzI2f7Fixepd+/eNGnSJNqzZw/nsU9ISKABAwZQu3btOOv3ERHFx8dT06ZNad26dZw6oopag7t27aqxwQdizbW3k3gXLVrE+SGp5Pbt27Rjxw5enbOzMy1cuJDXnZSamkqenp6UmpoqyMlUWFjIq/m7sCwrKD7LshQbG0t5eXm82sLCQnr8+LGg9vPy8tQmNFV9eBuOv5rEKCgoEOTWY1mWUlJSeGOnpaVRXl4eZyKt5NKlS1RcXEypqamcuoiICLpy5QplZGTwtn/27FlKTk7mLTQpkUgoMDBQ0DYRVewnoedRTY+hVCrl/ZyyLEsXL16k7t27qy3iWUlOTg598803tH379hr1o6aJV1yrQQUymQyGhoZwdnbm1LEsiwkTJuDPP/9E06ZN1eqICDt37sSaNWtEx5aIyDsiPz8fubm56Ny5s1qNQqHArFmzkJeXJ/jGHCDWXHsruLu7Y+zYsbw6R0dHTJo0iTPpAhVmiw0bNohJV0TkHdKyZUvOpAtU3JA7d+4cunTpgtoclIqZQAWXL1/G7t27AQDPnz9Hx44dVa7bMGDAAHz++ed13T0REZFaRENDA3v37q24JCBw/YyaIo5436CkpAQSiUS5NuyCBQvUTqvp0qUL70I6IiIiHx6ampqCF336O4iJ9w0YhsGOHTsAACdPnoSPjw9vhYNKYmJiarNrIiIitURRURGsra3rrD0x8b5BkyZN0LVrVwBA06ZNoaenJ+jnxqlTp3Dz5s3a7p6IiEgN4CshVUnz5s0RGBgoWP9PERMvB+Hh4di5cyfvzTMfHx8sWrQIX3zxBW9MhUKBPXv2cJoD3geEjPIr3WNCUCgUgrVCf2EAwh1bIv8+hJwnJ06cwMWLFwXFmzt3Ls6cOfMPeyUMMfFyEB0djR49evDqPv30U4wbNw59+vTh1GVmZkJbWxvh4eG8i0lnZWXBwcEBDx8+VKtRKBRwdnZW2my58PHxwYIFC3iLWGZlZcHExASLFy/mTWoHDx6EkZERr3OsrKwMurq6mDt3LqdrC6hwY02bNk3QB2Djxo2YNWsW77oXKSkp+PXXX+Hr68sbMywsDGfOnFFbl60SIsLdu3eRmJjIu58kEgnS09ORlJTEW8+LiJCSkoLbt2/j5cuXnFqpVAovLy8cP36cNwmlpqZi//798Pf359QpFArcv38f1tbWvIOD/Px8uLm5cboRWZZFamoqkpOTkZSUxOkKS05OxrFjx3hHnSUlJZg8eTJMTEx4LePLli2Do6MjHj16xGuB1tHRwe3bt2v0xf+3qcmk3w/9URMDRVJSEpmbmwvSLlmyhKKionh1fn5+1LZtW7pz5w6n7tatWwSA9PX1SS6Xq9U9ffqUWrduTc2bN6dt27ZxOoj2799PzZs3p65du9K9e/c429+4cSM1atSIli9fzjmhvaioiAYPHkwtWrTgLRtz69YtatasGU2dOpUzZnl5OU2aNIkaNGhAd+/e5Yx58eJFatWqFf3www+czsGUlBTq2LEjNWvWjDZs2MC5n1atWkUA6Pvvvyd/f3+1uoyMDOrQoYNyP3FNzD979iwBoB49elBwcLBaHcuyZGxsTABo/fr1nPuppKSEunXrRh9//DHvvr937x4xDEMzZszgjMmyLE2aNIk0NDQ43ZVERK6ursQwDA0fPpzXXailpUWampq0Y8cOTkfY3r17CQD16tWLXF1dOds/e/YsNW3alObOnctrunBzc6NGjRqRm5sbp46IyNbWlm7dusWrexP8G5xrAMYBeAYgDsA6Fa8bA8gGEPbXw0RI3Jok3qNHj5KLiwuvLi8vjwwMDHh15eXlpKOjQy9evOC1I65Zs4YsLCw4T+hKXF1dBelqSnl5OQUGBgpyjkVHR9PVq1d5PwASiYScnJwoOjqaN2Z6ejo5OztzfvEQVSSLtLQ0QS5DIaSlpfG6myq5c+eOoP2TkJBANjY2JJFIOHXl5eV0+PBh2rVrF2/My5cv08aNG8nHx4dTV5nMt27dSsXFxZxaPz8/2rRpEx0/fpy3/f3799OmTZt4HWmlpaV08uRJioiI4I2ZmZlJGRkZvLpKysrKyMPDgxITEzl1hYWFNHXqVDI2NuaNmZiYSIMHD+YdnLzJB594AdQHEA+gM4AGAMIB9HpDYwzgUE1j1yTxGhkZCUpoe/bsodu3b/Pq7Ozs/tY3qYiIyD+HZVk6f/487xd5fHw8tWrVivbv31+j+DVNvO+jgWIAgDgiSgAAhmEuAdADEFVXHSgvL4dCoeAtKKhQKODh4YFly5bxxjQzM3trk7HL5GWwj7DH6cen4Z+i+pqdZj1N6HbXhVl/M4z6chTqMeLlfJH/LgzD4Oeff+bVde7cGdeuXYObm1ut9ud9TLwdACS/9n8KgB9U6AwYhhkGIBbAciJKVqEBwzALASwEIGjWAQD4+vpi+PDhACqmo7Ro0UKlLjc3F7/88ougidb/NOnKWTmW3VmGw8GHBellrAxXo6/iavRV5XO63XVxcPxBfPGRsP0gIvJfZNiwYbzW4n/Ke7dIDsMwPwHQJiKTv/6fBWAAEf3ymqYVgGIiKmcYxgyAERGN4ostdJGc5cuX45dffoGWlhbWrl2LvXv3/u3t+adIFVIM/WMoglKDqjyv30Mfpv1MMabzGNSvV732W35ZPq48vYKjoUfxKL36mqXbRm7DxqEba80SKSLyX+LfsEhOCoDXF0DQApD2uoCIcoio/K9/TwDgL/1QA/r27YvOnTvj8OHDCAoK4n9DRZ84F87+O+z024mGVg2VSVevux5KNpSANhOcpjph3FfjVCZdAGjZqCUW9FuA0IWhoM0E6W9S7B6zW/n6Jq9NqGdZD8bXjaFg/93VGUREhHD79u06a+t9TLzBALoyDPMlwzANAEwDcON1AcMwr9fi0QUQ/TY7MG/ePLAsC39/f17zRCVHjx7Fixcv3kr7ZfIyMFsZbLi3AQAwrc80sBYsrk+7jiaa3Ned1aFZXxOrBq0CbSZkr87GDx0qrt6cDT8LjW0aWOm28q30XUTkQ8Xb2xtPnz6tk7beu8RLRHIASwC4oSKhXiaipwzDWDIMo/uX7FeGYZ4yDBMO4FdUzHJ4q9SrVw+tWrUSNJE/OjoaK1asQOPGjQXF5rrckZSfhMbb/3+czFWZuGhw8a1eEmjdpDUCTAJQtrEM474aBwCwDbAFs5WBe7w7anL5qTa079vlL5H3EyHnyfPnz3HgwAFB8YyMjHD58uV/2i1BvHeJFwCI6DYRdSOiLkS0/a/nLIjoxl9/ryei3kT0DRGNJKJaWZ0mPT0d7du359U9fPgQn3zyCa8bjYiwbds2tU6viMwIdNrfCQAw6stRyDLPQkN5wxr3WxUymQx3796t8lxDjYZwnemKwnWF+KTxJwAAbXttfGL1Cc45nOONmZCQgJ9//llQ4U47Ozts2rSJV1dUVIR58+bh8ePHvNro6GhYWVnx6gDAy8tLUC2t0tLS96ICdFlZmaACnXK5XHDh1KioKEFrESgUCoSHhwuKKZVKUVBQIEgrpF2hsc6cOYNt27bx6r766ivcu3dPkLuzX79+CAkJqZMv/vcy8b4PFBYWonnz5oK0ubm5cHBw4K1G7OzsDAsLCwwZMqTaa89znuObo98AAKZ3nI52d9ph/vz5amdUVDJx4kTo6urCy8uLU3f9+nWMGTMGFhYW1SyRzRs2R86aHAQvCAYA5LP5mBM3B6e9TnPGZFkWHh4eguzFhYWF2L59O1xcXDh1RUVF8PHxwZw5c3jtoGfPnoWFhQXvrxKpVIrp06djyJAhvDbcp0+f4vPPP8e6deuQl5fHqV20aBF0dHRw7do1zu2Pi4vDzz//jKlTp3JahhUKBU6dOgVdXV1YWFhwzpbJzs7Gtm3b0KtXL94EHRISgilTpmDz5s2c51NpaSns7OzQvXt3Xmv5gwcPYGJiAh0dHc6lUYkIFhYWmDRpEi5cuMAZMz09HZ06dYK5uTlnYUoiwqNHj7Blyxbe4qkMw0BPTw/9+vVDamoqr7ZPnz51c7mhJpN+P/RHTQwUHh4edPDgQUHaiRMn8k7MJqqwI1pYWFBRUVGV53MluYQtIGwBWXpbko2NDbVv356Sk5N5Y164cEGQXVkul1NGRgY9fvyYXr16pVbHsiyNPjta2Z8FNxZwxi0uLqaDBw9SQEAAbx8iIyNp06ZNnIU2iSrqaF25coXXWk1UYd29cOFCtX36JoWFhfTs2TN6+vQppy49PZ0CAwMpPz+ft21HR0d68OABry4tLY3Gjx/P67AiqnAtamtr87obMzIyqGXLlrR3717emNu3b6e2bdtSUlISp668vJy6detGc+fO5Y158eJF+uyzzyg+Pp5Xu2rVKjp//jyvrqCggLy8vAQVjiWqOJ9sbGx4P3vPnj0jLS0tunLlCm/M4OBgmjt3Lm8tuzfBh+5cq81HTRKvpaUlhYaG8uri4uJoyZIlvLrs7GzS1dWtZqtlWVaZ5IyuGBFRRZIUknRrk3sJ95T9whaQVK662q6IMIR8MRNVrMEgJOkTESUnJwsqtros2L8AACAASURBVFlYWEhpaWmCYubn5wsuMMqXyCsR0sfaJiEhQdBAyt7engDwVsJ+k5om3vduHm9tInQeLwAYGhri0qVLvHXSdu7ciaFDh6q8fPA6crkcr169wqefflrl+VFnR8Er0Qv1mHpQWLxf07ryy/LxsfXHyv+zV2ejdZPW77BHIiJ/H7lczvt5JiIYGBhg5cqVGDx4sODY/4Z5vO8clmWhUCigoaHBu0RcVlYWBg0axBtTQ0OjWtJ1j3eHV2LFtdnSjdzLEL4LWjZqCdpM6NG6YmnMNrvbIDr7rc7cExGpM4QUm2UYBidOnEDLli1rtS9i4lVBTEwMevbsiYyMDNy5c4dTu3fv3r9Vm4mIoG2vDQDwn++PBvXf39pt0YujMfub2QCAXna94J/MvaariMiHTKtWrdC7d+9abUNMvCqIi4vDsGHDsGfPHsF3OGt6yWbmtZkAgM4fd8aPWj/WuI91zVn9s0rn26DTg8TkK/KvQyKRvHX3qTrExKsCXV1dDBgwAPb29oLdaA4ODoLjF5UX4eKTinIkUeZ1tujaP2bVoFU4MK5iMvqg04MQkBLwjnskIvL28PLy4p3y9rYQE68aPv74YwwdOlTQpP/4+HhYWFgIjj34VMVF++U/LkdDjbdjkKgrfvnhF2XyHXhqIBLzE99th0REeLh27Zog3eDBg/HgwYNa7k0FYuJVA8MwkMvlvKYIALC2tkZ2draguCfPn0RkdiQAYM/YPZzat137SYgTCqiYSB8bG6v29V9++AWbhlV8IX25/0s8jeO/HJORkYE//vhDUPvR0dGIiRFmRuSbFF+J0G3/EBG6bTUpsMpnXnnX5OfnV3NiqsPZ2ZnTkFFJy5YteY0jbwsx8apBKpVCU1NTkLZ///7Ys4c7iQJAaGgoTL1MAVSMdrnWXwgMDMTOnTt5Y0ZFRcHV1RW5ubmcOoVCAR0dHTx79ow35vbt26Gnp8dZHNBypCVGfz4aANDHoQ9Y4v6SOHfuHExMTHgTpUQiwdSpU2FjY8Pbz9DQUPTo0YPXjQYAx44dE7Q/AWDXrl1ISkri1SUnJ+PWrVu820REiImJgaOjI2QyGa/2zp07uHHjBqcOqEg+S5cu5T32AHD37l3s2rWLV1dWVoZVq1YJOk+ioqKwa9cu3vsbOTk5cHFx4fwyr+TYsWOIiuK//GZnZ4cpU6YIshgPHjwYOjo6gq7famlpITlZ5dLebxUx8aohLi4OXbt2FaS9f/8+jIyMeHXR0dFgO1YkKJux6hOLp6cnBg4ciLFjx/LGzM7O5lysvRKiinUiIiMjeT8oW7duxerVq3Hv3j1O3d15d9FMsxkAQMtWi1O7evVqnD17lvdnX5MmTeDt7Y3GjRsjMzOTU9uvXz+4ubkJGvno6uriyy+/RFpaGqeOiNC9e3ekpKTw7qdXr14hMTERbdq04Y25du1afPbZZ7xf5g4ODliwYAFGjeJeXrq4uBhDhgxBaWkpb/uXLl3CmDFjMH78eE5deXk5Jk6cCH9/f/Tt25dTGxcXh6+//hoDBw7kXcDp/v37cHFxEfR50tLSQmZmJu8ofv369Th06BCuX7/OG7NTp06IiIjgPZ8AYOjQobh9+3btj/hr4rb40B81ca79+eefZG9vz6urrMwqhOGrhhO2gHoc6sGpk8lk5OTkJChmbcJXvJKISMEqlO627b7befU1KczJZ5v9kJBKhTn/5HK5YEdYaGgob4VhooqCrLa2toKO5+PHj8nd3V1Q+0LdXQqFotbca0LPkV27dtHDhw95dRs3bqR69epxVqJWBUTnmnpq4lyzsrKCjo4OvvvuO05dYmIi9u7di/379/PGZLZWjAxeLnuJzz/6nEf94ZBSmILP91ZsT+SiSPRp2+cd90hEpCpEhJycHLRuze28fPnyJb755hveBZLeRHSuvSViYmLQvXt3Xp2vry+GDRvGqyuT/z/2zjyupvyN45/ToiIxdtkGjWUGY98GY4xtSlLZyjZiYsi+M8NkEBnbYBgMKltEGHuyJpKUpagopX3fu7fuvc/vjzs10r3nfC9F/ea8X6/7ktun7/mec+99zrnf83ye598bG/9PQRcAGhs1xmErZTpd+13tNc5pFhEpbziOEwy6gLIvI+u9gPdBDLxqyMvLQ0hIiGAQGTVqFMzMzATHW3ldmW426ctJZTK/ioZte9viOg4Wxyw+8mxERN4de3v7ct+GGHhVoFAowHEc1qxZI3jH2sDAQLAAOgBs9FW6vrYO3Vomc6yIxM5XHqu/w/7Gk8QnH3k2IiLvxruUANB4G+W+hUrKsmXLcO7cOQQEBDDp+Sr2v3nVXFO/fItvfEyqaFfBqdHKrIUOuzt85NmIiGgGEWHhwoUfZFti4FWBlpYW6tati6+//hqNG/OnSQHA48ePsWPHDrW/Pxem7LrQqUGn//v1T8u2lsU/HwhkM0yIiJQnycnJgvnTgHIduKwa1gohBl41NGvWDJ988gm6dBHuHL9hwwbe5PAlV5codQM3wMXFhelN8LFhPUEQEfLzS5a0jJyjfPPanbUrMU5OTg7z9vnMGyIiAJg/RwcPHoS/vz+TluO4MneMqkIMvGpgKZoMKANP3bp10atXL7WaZynKGrbdanfD0qVLmZopSiQS+Pr6sk+YAX9/fybbaHh4OGxtbQV1CoUCY8eOxebNm0s8/2nNT2FSywQAsPb2WgBKl5OxsTFSUlIEx92wYQNWrFghqCssLMSsWbOYAnpubi5vd+c3YWmK+S6wNtHMyclhas4IAF5eXky69PR0XLlyhUn7999/M534FAoFHj58yDRmYWEhkylBIpEwnfTDwsLQtm1bpnm2aNECmzdvZrIDN2rUiPnYvw9i4FVDUlIS6tWrJ6jjOA46Ojr44YcfBLVubm5ITEwU/ABKJBJYWloy5RJKJBLs37+f6Sydm5sLHx8fQV3z5s1hZGQkaEXV0tLC8uXLVVp2/ab6AQB+vq6s6dC2bVtMnjyZyYo6fvx4QYcZAOjq6qJXr16CDQ8BZf2JO3fuMB2n+/fvM82TiLBjxw6mKy93d3dmG/CAAQNQrVo1Qa2rqyt2794tqMvPz8fw4cOZgtS1a9fw008/CV50yOVy2NnZMdmAc3JyYG9vz3QhExISwtRh+rPPPsN3333HVNMjOTkZJ0+e5G3K+ea44eHh5b8kqInborI/NHGuPXjwgNauFXZiERGNHDlSrYMmJCmE8Auo3sZ6dP36dVqzZg1lZWXxjldQUEB79+6lhIQE5vmWNQqFgtllFhMTo/L5WhtqEX4B7Xmwh4iUDqbo6GimMVn7fhGxu8I+JgUFBcx913x8fJh0t27dor///ltQl5aWRosWLaKIiAimbW/atIl5+6zvUaEGp++CQqFgep/IZDLq0qULk3NvypQppEmcKAKic009mjjXzp07h6SkJNjZ2Qlqrays1NYgmHtpLrb5bcPKfivRPrU9WrVqhQ4d/ht3/N90tNGq/877TKTiERQUhI4dOwrq/vjjDxw6dEjjZT5NnWvC1/7/UeLj45kyGuRyOW/e37GnxwAAkzpOQotPWpTZ/CoDjY3+PX7x2fFoWF24xKaISHnAEnQBYPLkyXj16lX5TgbiGq9a4uPjoacnXKRcJpNhzpw5an+fmKusiPRfC7pFbB6svPFmc9LmI89EREQYAwMDphu774sYeNUQHx+Pffv2Cer09PTQt2/fDzCjysmcnsqT0s2omx95JiIibNSoUaPctyEGXjW0bNkSp06dYkq/kkgk5ZaCVNnR4v59i6XkCaeSiYh8TDw8PD5I9w0x8KqhZcuWaNSoEV68eCGovXDhgsrF+Dcrkv2Xmd9zPgBgvY9wBwQRkY/J1atXkZWVVe7bEQOvGszNzdG1a1e0aydcW/bw4cO4ebP0V+ngJGUvsg71/81iYK398P/Egt4LAACb7gq3RxIRKWtCQ0Oxfft2Jq2hoaFGDst3RQy8apBIJExVx6RSKWJjY1WaDYISggAAHRso76gWFhbCwcGBaftyuRypqakazJgNlvRBmUwGDw8PJmOAl5cXb50KADCublz8c1BQkOCYL1++ZDIGAMoPFetxqiypk5o0pWRpZwMozTOsWpYedkWwuvHKGplMhq1btzL1R9PS0sLjx4+Z+rNVr15dDLwfE4lEAgMDA0Gdnp4eBg8ejJkzZ5b63ZMkZWnE9vWU/avOnTuHe/fuISkpSXBcd3d3ZovrhQsXmHTPnz9nal+to6MDV1dXJqdR9erVeSuzvY1Qfy4AqFatGnMH2StXriAiIkJQR0T466+/mMaMj49nOkEA7Mdek664W7duLVX/QhVZWVlYtWoV05h79uzBkyfCpTrlcjnmzZvHNOadO3dw8eJFJu2ZM2eYdHl5ebhx44agTkdHBz4+PkwlHD08PLBv3z4m556hoeEHOZmIebxqYL3iBZRn/erVq5d6/nWW8mzcrEYzAMqme/369WO6Effq1Sum1vIAYGpqyqRr06YNcyvwHTt2MAXJnj17Mr35q1epjuyCbBjUFT6ZNWjQAI6OjkzznD59OlP3WI7jMHnyZKYxGzZsWObHvmrVqujUqROT1tjYGIWFhYIn/ufPn6Nq1aqC4xER7t+/j+bNmwtqQ0NDkZmZqXRXCbz+ycnJzMfpq6++YtJVrVoV/fv3Z9KuXr0aDRo0ENRNmjQJLi4uqF27Nq9OoVDgyZMnSElJ4a29UhaIgVcNmgTenJwcGBoalno+OTcZAFCvmrLmQ/369dGhQwd8+umngmMuX76cfbIaoK2tzaRr2rQp85jdu3cX1Aw1GYoTISfgn+KPVsatBPVffPEF07Z1dXVRsyZbjWPWfS8PqlSpIvjBL2LixIlMuu7duzMde47jcPToUaYxP//8c+Yr8xEjRjDpADC13dGUzz//nElnbGzMdCLX0tJCWlraB6keWCGXGjiOG8pxXCjHcS84jluq4vd6HMe5//N7P47jPi3rOeTn5zMtNQDqr3iT85SBt241Zfvthg0bYsuWLWU3ySIOHwY+/RTQ0lL+e/hw2W/jPWldW9m/Ljwt/CPP5B2oBMe3UvMBju/o0aOZdLa2tvj666/LfPtvU+GueDmO0wawE8AgADEA/DmOO0tEIW/IpgBIJyITjuPGAtgAYExZzkMikUBXV5dJ27lzZ5XroWn5yhtutQxqAfi3klmZcvgwYG8PFC1fREUp/w8A48aV7bbegyLnXkS68HpshaKSHN9Kywc6vizLZoAymyk5ObnMtquOinjF2x3ACyKKIKICAMcAvN090QKAyz8/ewD4lmM9soxEREQwV6NfvHixyuflCuV6qjZXjl9xV6z4901bRF6e8vkKRFHLo0yp8J3lCkUlOb6Vlgp2fKtVq8a0FPi+VMTA2wjAmzkiMf88p1JDRDIAmQBULqBxHGfPcdwDjuMeaHIme/LkCVOtUQBqb1jpaiuvmAsV5bhmpC71R4OUoA9BkZnEQIdt+abCUEmOb6Wlgh3fjIwMJCQklPt2KmLgVXXl+nYCJotG+STRHiLqSkRd69atyzyJmJgYvHz5kimncs+ePSpzRHW1/gm88nIMvOpugmlwc+xDkFuozDww0K1kgbeSHN9KSwU7vnfv3mVOEXwfBAPvP2uoH5IYAE3e+H9jAG+3IyjWcBynA6AGAP52CRqyYMECfP/990xV6w8cOIB79+6Vel5fR5kVkS/7NyezzNeP1q4F3k4pqlpV+XwF4nmKslNAq1rCGQ0VikpyfCstH+j4Hjx4kEmnLkOprGG54nXlOO4ax3FsuRvvjz+AzziOa85xXBUAYwG83TPlLIBJ//w8EsA1KmNbUvXq1SGVSgVzVLOysvDkyRN4enqW+l1R/dn47HgAyuRw1vxU5maP48YBe/YAzZoBHKf8d8+e974xkZaWxuT0SkxMxL179wRTcALilVbpBhDOu5RIJHj+/Dlz369ydaSV0/EtK1hTnzQp4sTqnCuT4/4ex/fZs2dMbjQiwrlz55imk5OTw9R26X1hCbxdAOgCCOQ47jeO48r1dPDPmq0DgMsAngE4TkTBHMet5jhu+D+yvwDU5jjuBYD5AEqlnL0v+vr6TG9APT09zJo1S2Wnik9rfgoAeJXxCoCyToOrqytTwv+FCxeYCvQAwIPWrYFXrwCFQvmvmjdtdnY2Ll++zDTmtGnT8PTpU0FdYmIipk2bJpitcSvqFgAg65lwARItLS0MHDiQKaicOHEC+/fvF9QBSgcTS7BIS0sr6YYbN07t8WV1FwJg6g8GKB2OrAF127ZtTLqjR48yv5+cnJyYdIGBgXj8+DGTNjycJ43wjeMrDQ3F0y+/ZBrTycmJKed4165d8PX1RXBwsKC2wlzxEtETIuoLwB7AeAChHMeVa1VrIrpARK2IqCURrf3nuZVEdPafnyVENIqITIioOxGVeY6SJoHXyMhIpfbTGp8C+DfwpqWlwdjYmOkD+Pr1a+zZs4dprqxmBx0dHebawTNmzED79u0FdS1atMCKFSsE03UUpGwyOX3SdMExq1Spgq1btzLlUbdr1w7jGK8+v/nmG2bLMmvCvyZGE5YGngAQFRXF1L03IyMD27dvZzqRh4aGMgVpIsL27dv5A+Ub2z9+/LigDgBSUlKYvsHo6uoyu+HMzc1haWkpqDMxMUH16tUFDRcSiQRRUVGQSqVM238vNGnQBuVa6k4AMgDXAXyhyd9/7IcmTezy8/Np/PjxTNotW7bQ7du3Sz3vEuRC+AVk42FT/Jy1tTXTmIWFhSSXy9km+5FhmSd+AeEXMI/J0pjwv05hYSFzo8/k5GTKzMwU1MlkMkpMTGQa82O/RqzbT09PJ2dnZ0GdXC6nWrVq0Y8//qjxXKBhs0uNshqIKJOIZgLoBqAOlMsPmziOK23bquTo6ekxn/nUlZLr0rALAMA/zr/4uUWLFjGNqaOjw1QDoSIgNM9HCcoiOrUN2CyzAHvC+38ZHR0dZpNPnTp1YGRkJKjT1tZGvXr1mMb82K8R6/Zr1qypsojV22hpaeHrr7+Gra3t+05NeFssIo7jdDmO685x3GyO444AOAngCyidbzMBPH9j/fX/Ao7jmG8eGBkZqbx50bZuWwDAi7R/19Z69OhRNhOsRKy4pkyGX/ftuo88E5H/KizFhABlnYzevXuX82wYLMMcx/kC6ASgCgAFgEcA/gbgA+AOgBwAqwB4cBw3m4jYCqlWcFgLhQDqfeBvtr35L3M+/DwAYHJHtupgIiIfCwsLiw9yJc9SOCAHgBOUQfYeEalayV/AcVwigOUA/i8C75EjR5jq5orwkyn5N92nyMknIlJR+VDLJyxZDYOJaDUReasJukXcgtLs8H9BQEAAWJ1ut2/fFtQU3dX/rzHnkrLLcFHfNRGRig5rWc73oSy/Cz9C6WI2lRIigqWlJfr06SOoVSgUmD59utpUmc4NOwMAfKJ9Pkidz4qGyyNlLaO134pOL5GPx8uXL5la+qSlpTG5Vd+XMgu8RJRPRH+X1XgfE47jMGDAAKasgpiYGISEhMDDw0Pl78d8oaxWeezpMfj5+cHHx6dM51pesHaqAKC2qIjXS6/in/V19JGSksLsoGLJTQX+TYcU4Yf1GFWmY5mbmwuFgu2b5Nq1a5ny8l++fAkTE5P3nZog4t0fNTRo0ADx8fGCOi0tLZiZmak1G4xrr0zu3/twL7KysjB79mymoHb+/PkyDz5EBHd3dyats7MzU8J9ZmYmGjdurNLlNvjQYACA/w/KdLodO3Zg+PDhTPvfuXNn3Lp1S1D39OlTzJ07l2mfvLy8kJKSIqgD1J9M3oY1UKWnpzP1UQOUhVpYxo2OjmbqoQcALi4uTO8nqVQKNzc3pjHv3bvH3J+MtdEmAObXaMuWLbCzsxM8VpcuXUJISAjTST88PFwMvB+TBg0aMH34GjdujK+++krti9rISFnRUqaQISsrC69evVJZUOdtsrKyYGlpKZhLTEQ4duwY09eo8PBw1K9fn+kGQk5ODiZNmiSoCw0NhaOjI9q1a1fi+cD4wOKfuxp3BaB0Ou3atUuwBU96ejoGDRrE5LJ7+vQppk6dKrhPRITk5GSmtKLExESEhoYK6gClvZelbnPVqlUxc+ZMppOOj48P1q0TTr1r1KgRrK2tERAQIKhNT0/HmDFjBN1j+vr6WLNmDf7+m+3Lq5mZGVO++5EjR5gaiGZnZ8Pf319QByhP+mvWrBF87WNjY/HixQuVXWLeRCaT4cWLFx8k8H50N9mHfGjiXFMoFGRlZcWkPXHiBB0+fFjt74tcW+mZ6TRx4kRKT08XHDM6Oppu3brFPN+yhtURlZaWptJBVLTPl8IvFT+XkpLCNGZubi6za+9ju6c0gXWfkpKSmFxmRERHjx6l/Px8Qd39+/fJ3d2d6XU9ePAgRURECOoKCgro/v37TPMsj9cpOzubSXfixAnatWuXoO7Vq1dUo0YN+v333zWeCzR0rlW41j8VBU3SSkxMTHD27NsF1P5lUItB8IrwwtnIsxg7dixSUlIEGzQ2adIETZo04dWUJ6yOqE8++aTUc94R3sU/DzEZUvwza7NH1mR34OO7pzSB1YmoSd3osWPZqrZ269YN3bp1Y9KyfNMBlO8R1jHL43ViLWbTtWtXps9SkyZNUFhYiMGDB7/v1AQRA28ZYGJiwlv5ae2AtfCK8MK8y/OQujj1A87s4zDQbSAAwG+q30eeiYgImFv5aGlpYeHChWjdunX5TgjiGq9aIiMjoa2tzVRRydDQEPPnq89T7dZIeVVQ1Pzy/5lpf08DANStWhfdGwm3HhcRqUgsXVrmFWZVIgZeNcyePRt169Zl7hjRsWNHJt3LtJfvM60KTUxWDPY8VJayjJzD1ihURKQiwVKKtCwQA68KCgoKcOHCBbRq1UpwHSk5OZkpP3Cn6U4AwKTTbOtnlZEmW5TraJsHb0a1KuVfxV9EpDxwc3NjTml7V8TAq4LU1FR069YN9erVE0xBiYiIwKZNmwTH/LHrjwCAO6/Z8i4rG2M9/r3JM6/XvI84ExGR0hARsyno1KlTKm8alyVi4FVBw4YN8f3336NWrVqC2szMTKxbtw4xMTG8Oo7joMMp72XejhKu7fCxISJmV9DZ4LNwD3YHAOQtV5+kLpPJmJ1rBQUFTDqR/y5ExPw+uX37Nn7//XdBXUFBAbS1tQVzzd8XMfCqoW3btnj27JmgTl9fHy1atEBGRoagdqreVABAv4P9eHVEhLt377JNVEOuXr3KdOb38vJiaswZnx0PCw9liQ7vid687dsvXLiAXr16MbXAcXBwwJ9//inoSsrLy2M2JkRFRTG9poDScFCWsF5tSaVS5hZBgYGBTMYZImWzR6FjCSgNKa9fv2bafkhICJOO9QReRE5ODtPfXL58GdbW1oLHICIiAvv27WNKvXv48CG6dOnCPNd3RQy8amANvP369YO5uTnTmL2N/y2wnJ6v/oPNcRwuXrwIa2trxMbG8o6ZkZGBGTNmYNeuXYLb9/T0xMaNGwUdeUQEJycn6Ojo8LqS5Ao5jDcbAwAmfzYZA5oP4B13//79+PbbbwW/Sbx69QrXrl1D06ZNBfM/nZ2dkZmZKbgml5iYiHnz5sHPTzjF7fr167C0tBS8Os/Pz8eyZcuYTlB3797FN998g0ePHvHq9PT0sHTpUjg4OAhuv0aNGujQoQMuXrzIq+M4DpcvX4aNjY3gCaVZs2b46quvcPjwYcFA7ePjg6FDhwr2EMzLy8P48eOZmmgGBgZi+PDhTCefHTt2oH379oKNVs+cOQNPT0+mi6Pbt28z9yV8LzRxW1T2h6bONUtLSybtxYsX6Y8//hDUvXr1iros7EL4BdRtTzdebWJiIvXq1Yup/5VCoaCcnBxBHatzKi0tjeLj4wV1Re60xpsaC2rz8/Pp2rVrTNt//PgxkytJoVBQTEwM05jl1b9OoVAwORGJiG7fvk25ubmCuujoaFq5ciXTmFOmTKEDBw4IOsOCgoKoefPmFBoaKjimg4MDLVy4UNDlVlBQQGPGjKG0tDTBMTMyMujhw4eCOlZyc3MpKCiISbt06VJydHQU1D148IC+/PJL8vf313g+0NC59tGD4Yd8aBJ4iYg58GZmZjI3xgwODi4OWDlS/mDJatv9GDT8rWHxflQm2+7/G5q8R1gttiwnhyIqw2t/+vRpphOvn58f6ejoUEhIiMbb0DTwiksNPNSuXRupqcJOMyMjI2RlZSnPZAJ8/vnnxRkOPfbx919jte1+aHr/1RvxOcrKbQU/FVQq2+7/G5q8R1gttv9vlm0LCwsmu3aNGjUwc+ZMtG3bttznJAZeFeTn50OhUKBNmzbMN2NGjRrFXPavKKc3ODkYsVn8a7gVjc5/dsbdGOWNv6ylWWI7H5H/Gxo0aIBVq1Z9kG2JgVcFL168wJEjR9C2bVvBGwdFjB8/nvlKgeO44uDbeEvl6ZbU8veWCExQlnvMWJKB6nr8Oc4iIpWJGjVqlHv+bhFi4FVBamoqli1bhlatWgnegWdxraliRrcZxT+v91n/TmN8KGQKGThHDhHpEQCA7GXZqKFf4yPPSkSkfMjMzGSqHfw+iIFXBZmZmZDJZIiPj4eVlRWvdteuXUydKlSRsECZ1rXMexmiMqLeaYzyJi0/Dbq//rucIFkhgWEVtrVCEZHKyOnTpxEWFlau2xADrwrMzMxgbm6OVq1aCWpzcnIwefJkpoTv4ODgEv+vb1gfGwZuAAB8uu1TyBXsfc4+BFcjrqK2s7KGbqPqjaBYqYCejt5HnpWIiOYQEX7++Wcm7fnz52FmZlau8xEDrwp0dHTQu3dvJvdYjRo1EB4ejoiICEHtypUrce3atRLPLf5qMT6r9Zlyu78qE8HT09PfeQlDiNDQUKZE8pFHRmKQ2yAAwLQu0xAzP0btHezNmzfDxcVFcMzY2FjMnDmT6WriwIEDOHDggKCOiLBpBTn0lQAAIABJREFU0yamnl6FhYVMbXIApduKteddWcJShhRQ7reQuaaI+Ph4ZideeHg4U3ZOfn4+81zfBZYb1YmJiZg5c6bgN06ZTAZ3d3fUq1dPcMyMjAzo6uqiWrVyLvKkSe5ZZX9oksf7/PlzWrx4saBOKpWSnZ0dk+HA0dGRatSoQc+fPy/xvEKhKM6JretclyQSCfXv358WLlxIL1684B3z3r17ZGdnx9Sm6Ny5c2RkZESrV69Wq3md+bp4LvgFtNdrL++Yvr6+BIDmzJkjmNM5duxYMjQ0FMyTjIqKIiMjI/r2229JKpXyards2UIAaOvWrby65ORk6t+/PzVq1IgKCwvV6iQSCe3du5dat25Nx44d4x0zODiYpk2bRqamprzzlMlktGHDBpo8eTIdOnSId8zw8HAaMmQIOTo6lnqfvM3GjRvJzMyMLl68yJunKpFIqEuXLrR06VJBQ46npyd169aNPDw8SCaTqdUVFBSQtbU1OTg4kK+vL++YDx48oPHjx9OCBQsoKytLrU4mk9Hhw4epR48e5OzszDsmEZG1tTVVrVqVHj16xKvbsWMHtWvXjuLi4gTHdHFxEXzdVQHRQFE2gVehUNCIESOYtOfPn6edO3cK6gICAsjW1pYCAgJK/U4qkxYHO8N1hhQQEEAGBgZMjjh/f3/asWMHr0YikZCXlxc9fvyYUlNTVWpmnJtRIujGp/KfTGQyGS1fvpxu3rwpGHQfPnxIK1asoKSkJP6dIaLff/9dMJgQKd19e/fupejoaMExHz58SPfu3aPo6GjegKJQKCguLo6uXbvG5PJ6+vQp/fzzz4L7X1BQQDNnzhQ8QRAReXt7k5GRET1+/JhXJ5PJqEePHtSnTx/efSJSuit1dXXJ3d2dV6dQKMjc3JxatmxJeXl5vNqYmBhq0KAB/frrr7w6IqLr16/TgAEDeE96RErX5I0bNwRdbvfv36eNGzeqfS+/yejRo2no0KGC+/Ps2TOytLRkNpq8iRh4yyjwEhGNHDlS8IqLSHnVyxqkk5KSyMLCQuUHNa8gr0TgC3gaQBKJRKM5vwtnnp8psd1t97Yx/Z0mrqXK4HD6ELA6zSIjI5l04eHhzE6zZ8+eMTXRjIqKYjpBEim/GbK+tkJBt7xwcHBg+hz98ccfpK+vz2xFfhNNA6+4xstDx44dBYuaAECVKlVQu3ZtpuyGunXrYv161eljBroGkK/89wZbF48uOPrsKPuENeTGqxvgHDlYHFNWF/tE/xNkLc3C7B6zmf5eE9dSZXA4fQhYnWasfcJMTEyY88fbtGkDIyMjQV3Tpk2ZG262bt2a+bUVKmZTHhAp7wHo6QnfFJZKpVi4cCG+/PLLcp9XhQq8HMfV4jjOi+O48H/+VZnNzHGcnOO4oH8e6tv7viPKExjQq1cv5vKMo0aNwu3bbHV227Rpo/bNqsVpgVZRsa148pnJ4Bw5BCcFq9RrChFh5/2d4Bw5fOPyTfHzT398irQlaaIpQuT/Co7jUKVKFSZt+/btsXLlynKekRKuKMhUBDiOcwaQRkTrOY5bCuATIlqiQpdDRBonk3bt2pUePHggqHN3d4elpSUKCgrw448/ws3NTdNNlQkv0l7gs+2flXjur+F/YXLHyRpdQRIRzoSewfenv0emNLPE7+5OuYuejXuWyXxFRP6rcBwXQERdWfUVrb27BYD+//zsAuAGgFKBt7wJCgpCaGgoVq5cyVRco7wwqWUCWkU4G3q2eDlgytkpmHJ2CgCg/6f9Yd7KHH2b9kUDwwaoZVAL+bJ8pOen49jTYzgQdACRGaWbTnZs0BF/2/yNxkaVx64sIvIhyc3NLdeUsop2xZtBRDXf+H86EZVabuA4TgYgCIAMwHoiOs0yPusV7+LFi+Hq6oqAgAA0atSIV+vm5oYJEyawbP69iUiPgKW7JR4nPtb4bzs16IRtQ7ehb7MPUORZRKQSI5FIMHr0aJw9y76KWeGveDmOuwqggYpfrdBgmKZEFMdxXAsA1ziOe0JEKvumcxxnD8AeUN40YGH27Nl49uwZjI2NBbV79uyBgYEBRo4cyavLyspC9erV3+smU4tPWuDRdOXNPqlMCq8IL5wPOw+/WD8k5SYhXZKO9d+uR039msgrzEOvJr3QoX6Hd96eiMj/E/fv30fjxo0FP9fu7u4YPXp0+U5GkxSI8n4ACAXQ8J+fGwIIZfibgwBGsoyvSTrZ/PnzmfI4bW1tqW3btpSRkcGrS0hIoNGjRwvqiMo39UoikfAmsRchl8spNjaWSZuQkECPHj0SzLstLCykuLg4pnSduLg4io2NZUrni42NZU5VSklJYdIRlV/XirJAk/dIfn4+k66goKDCpP2xzEMmk1FcXBxFRUUJalNTU2nw4MGCqXcKhYJMTU01TuNEJU8nOwtg0j8/TwJw5m0Bx3GfcByn98/PdQB8BYCt654GmJqa4sKFC4K6DRs2oHnz5tDX1+fV1a9fHxkZGejTpw+Sk5N5tZcvX8aoUaOwd+9e3v5oMpkMfn5+2LZtG+zs7ARtoY8ePUL37t3h4ODAqyMiODg4oG/fvrhzh78dfWhoKLp3745NmzYJdnxdsWIF+vTpI5h2l5qaigEDBjDt05kzZ9ClSxfBOqr5+fmYO3cuvvrqKwQGBqrVERHu3LkDGxsbLFq0SK1OoVDA19cXu3btwowZM3jtu7GxsXB2dsbmzZvh6enJO08XFxcsWLAArq6uSEpKUqtLS0uDjY0NVq1aBX9/f94xL126hO+++w4HDhxAZmYmr3bSpEmYO3cubt++XZzdowp/f39MmjQJO3fuRHh4uFqdQqGAm5sbfvvtN2zdulXw9Xz48CGmT5+OSZMm8eoAYOnSpejdu7dgyueNGzfg4OCADRs2CKbe3b17F127dmVKP3svNInS5f0AUBuAN4Dwf/6t9c/zXQHs++fn3gCeAHj0z79TWMfX5IpXIpGQtbU1k9bV1ZUOHjwoqPPw8KBhw4YxtRZZt24d6ejo0L59+9RqFAoF+fj4kJ2dHTVt2pSePXumVhsaGkr29vY0depUOnHiBO+Ya9asoaFDh9Kvv/7K6+LJzMykHj16kK2tLd26dYt3f06ePElt27alVatW8SbxSyQSMjc3p5EjR5KHhwfvmF5eXjRkyBCaO3cuPX36VK1OKpWSk5MT2dnZ0YIFC3ivAMPCwsjR0ZEsLS1p+fLlvNt//vw5OTg4UN26dQX3/8aNG9S0aVPq3bs3r8tMLpeTvb096enpCbrMnj17RrVq1SIzMzNBt9X8+fPJwMCAvL29eXWJiYnUuHFj6tu3r+C3M1dXV9LV1aUFCxbw6lJTU2ncuHFUp04dunLlilpdbGwsOTs706hRo2jixIm8Y+7bt4969OhBjo6OgqaQAQMGULt27QR1eXl5NHHiRIqNjeXVqQKic61sAi8R0ZgxY5iaSEokEjI1NWX6epSYmEhDhw4V/PqnUCjo4cOHgh/oInJycsrka6JCoWD+ip2cnMz8lezVq1dM88vJyanQveZUkZGRIWjZLdK9fv1aUFf02ickJAhqHz16RHfu3BFcEiooKKCAgADBYE6ktEH7+flRWFiYoDYwMJC5R9m79DJThUKhYGqwSaRcsurZsye9evVKULt06VJq06YNs2vvTTQNvBUqq6G8Yc1qKGLPnj0wNjbGsGHDBLWOjo749ttv0adPH0Ht9evXUb9+fXz++efMcxEREdGc169fo169ekxLB1ZWVvj000+xefNmjbejaVZDRVvjrVB89913uHjxIpN22rRp2L17N5P2m2++EYOuiMgHoEmTJszrte3bt8fGjRvLeUZKxMCrgocPHwJQvmivX7/mvclQRIMGDbBu3brynpqIiEg58dNPP0FbW/uDbEsMvCqYO3du8Z13TRpesuYJi4iIVDxYCxiVBWLgVUF+fj4mTJgAhUKBefPmqU24joyMRFpa2geenYiISHmzfft23jS590UMvCqYNWsWmjdvDiJCgwYNUKOG6o66BgYGGDJkCFMrnf/STUwRkYoIETG1iMrMzMSVK1dgYmJSbnMRA68KJk6ciOTkZMjl/M0nGzRogNjYWPzwww+CzS7d3Nzg5eVVltOsUBCRoIGiCE36yYknLBEhWBrNAoCTkxPOnCnlySrFpk2bsGDBgnKtIV3RqpNVGIYOHYrLly/D3NycV7dr1y4cOHBA8EUaOnQoPv30U9ja2uKPP/5QWyNUJpNh69atMDAwQMeOHdGtWzfeeqK5ubnw8/ODTCbD4MGDeeeQlJSE33//HR06dEDPnj1516T37duHtLQ0REZGYteuXWp1ycnJ2LlzJx4/fgx7e3sMHTpUrdbV1RVxcXG4e/cuPDw81K6pSSQSODk5ITk5GZ06dcIPP/ygdszHjx/Dw8MDsbGx2LhxI2rVqqVSR0RwdXVFYWEhjI2NYWpqqnbMpKQkuLm5YeTIkahfv75aV2Jubi6CgoKQmpqK4cOHqx3v6dOnSEtLg0QiQatWrdQWOc/Ozsbhw4dRp04dNGrUCD169FBbHc/b2xu3bt1C586d0bNnT9SvX1+lTiqVwtHRETVq1EC/fv3QvXt3tTeQgoKCcODAAfTs2RN9+/ZF48aqq9cpFAq4uLhALpfDxMQEXbt2haGh6iqtqampCA8PR25uLnR0dPD111+r1BWRkZGBY8eO4dtvv0Xjxo1hYGCgVuvi4oK4uDgEBATg+PHjao9VVlYWzp49i8zMTCxfvpx3+8nJyQgODsbq1at5de+NJkm/lf2hiYEiOTmZJkyYwKR1dHTkdeQU8dNPP9GgQYMEE+Pj4+Ppiy++oDp16tC9e/fU6qKiomjcuHGkp6dHHTt2pA0bNqjUKRQKcnV1pVq1ahEA2rx5M6Wnp6vVrly5kgCQqakpr5khLS2NOnbsSDo6OuTv78+7T6dOnSItLS0aOHAgr0EjNzeXBg0aRADI1dWVd0wPDw+qWrUq1apVizfZ/8WLFzRgwAACQAMHDqQnT56o1CkUCjp8+DA1b96cANCUKVPUtuBJSkoiBwcHql69OrVu3Zq3MWpkZCSZm5sTx3E0atQo3uaMd+7cofr161OLFi3o8uXLak0nCoWCli5dSgBo9erVvK9TXFwctWzZklq0aEERERFqdUREe/bsIY7jaM6cObx1MlJSUqh3796kp6dHR48eVWsgKSgooN9++40MDQ2pTZs2vL0BL168SO3atSMAZGlpSQ8ePFCr/euvv4jjOBo4cKCgeWXo0KH07bffMplc5s+fT35+foK6t4HoXCubwEtENGrUKKbGd1lZWUzONblcTmFhYfTdd98JjpucnEze3t60YMECwUCdnp4u2BiRSPlhzc/P5/1AJSQk0O3btykhIUFwfy5cuECXL18WtGJmZmbSnj17KCgoSPDN7+fnR97e3oJNDAsKCigwMJDi4+OZPlBFf8PiRCRS2ldZChplZ2cLnnSKuHDhAiUnJ1NycjKvLjo6mm7evEm7d+/mfQ0UCgV5e3uTq6srubi48I4ZGRlJt27donHjxgm6vm7dukXHjx+ndevW8ery8vLo8uXL9Mcffwi63GJiYiggIIDJQZaYmMjbDTktLY3c3d0pLCxM0GXp6+tLtWvXpgMHDgi+n7ds2UJjxox5JweopoFXdK7xcOTIERARxo0bJ6jdunUrmjRpAmtra0Ht/fv3ce7cufL/OiMi8h/H09MTffv2RZ06dQS1zZo1Q9euXXHs2DGNU8s0da6JgZeHnJwc2NnZ4fjx44JaiUSCESNG4Ny5c0xN/YhIbAApIlJBSE1NxeTJk+Hh4cHco+1NRMvwe5Kamlp8l9TQ0BD6+vqCZRwBQF9fHwsWLBAse1eEGHRFRCoWJ06ceKeg+y6IgfctoqKisGbNmuL/jx49Gh4eHkx/O2jQIOa22CIiIhWH2rVrl38N3jcQA+9bNGrUCKtWrSouWD148GC1KUqpqakfcmoiIiIfiJycHMGi9e+DGHjfom7duujatStevlS2cKtSpQrGjBmjUnvw4EG4u7t/yOmJiIi8B6z3tBwdHfHJJ6X67JYZYuB9Cy0tLVy9ehX37t0T1H799dcYO3Ys/vzzT16dTCbD4cOHIZVKy2qaIiIiGpKYmIixY8cKOlIfP36MtLQ09O/fv9zmIgZeFdSoUQN169bFixcveHWdOnWCra2t4M03HR0dJCYm4rPPPsOhQ4d4ta9evcLp06cRGBjIfKMuOTmZqVhPYWEhnj9/jsLCQt5ebgDw4sULZGZmYteuXcVX/6qQSqWIj4/HmTNneHuUAcqi1ElJSViyZAm8vb3V6ogIr169QkBAAObNm8drCU1JSUFycjLu3LkDX19f3u0/ffoURCRY6lMulzNXpMvMzGTy/0skEsTFxUEmk6nVvHr1CiEhIYiOjkZWVpZanUwmw8mTJ3Hz5k1ER0fzBpLnz59j79698PX15a0pQkTw8PDA2bNnERYWxjvPpKQk3Lp1CyEhIUw3nhUKBdLS0pjez0SEkBC2FoqvXr1Camoq1q1bh7t37/JqY2NjMXHiRKxZs4a39KNCocCKFSvg5OTENId3RpOk38r+0MRAERQUJNhLqggbGxvBBHKpVEp9+/al0aNH87bLKep5hn+cY+qS3RUKBXl6elL79u0JAI0ZM4Zmz56t0iCgUCjo7Nmz1KZNG2rdujWNHDmSDh8+rHbcHTt2UI0aNcjKyoo8PT3VuqIyMzOL+1lt2bKF15zg5eVF9evXJysrK15nkEQiofHjx1PTpk1p8eLFFB0drVbr7e1NxsbG1Lt3b1qzZo3a5PyYmBiaNGkSGRgYkIWFBTk4OKjsNpuamko//fQTNW7cmFq2bEnW1ta0fft2lWNev36dunTpQhzHkYWFBS1fvpzi4uJK6XJzc2nJkiWkra1NtWrVon379ql9ryQnJ9OwYcMIAE2ePJni4+PV7vvDhw/J2NiYjIyMeI+nQqGgdevWEQCaMGECb8up2NhY6tChA+nq6tLly5d5x3R2diYtLS3q0qULvXz5Uq32ypUr1KJFCwJA9vb2dOXKFZVzyM3NpbVr11LLli2pYcOGZG1tTb/88ovacT09Pal27dpkaWlJ3t7evA6/DRs20ODBg9W6EN9kz549tHfvXkHd20A0UKhH0zxeCwsLHDt2jNcvDgAxMTGYNWsWTp06xZsmlpeXhzt37uDPP/+Eq6srb8fT48ePQ1dXF0eOHMHw4cMxYcIElbr4+HgcO3YMI0aMgK6uLho2bFjqjJ6eno7bt28jNTUVHMdh0qRJaud55coVXLx4EfXq1YOFhYXaThl5eXlYvnw5cnNz0alTJ0ybNk3tlYSfnx+cnJzQsmVLjB07Ft26dVOpk8vl2LhxI6Kjo9GiRQv8+OOPqFatmkptZGQkDh48CAMDA7Rv3x5mZmYqdUTKrsHR0dFIS0uDhYUFmjRpolILKL89+Pn5ITU1VbDTbWZmJo4dO1Zcq6Bp06Zq3yvPnz+Hl5cXmjVrVpwzqm6+O3bsQIsWLeDh4YEpU6aobScVExODK1euIDg4GImJiXBzc1P7uh45cgS6urrYv38/Zs2apbZWRUZGBjw9PREWFobnz5/j0KFDal+Da9euIT4+Hj4+PmjTpg3mzJmjUieRSLB9+3b06tULjx49wqhRo1CvXr1SupycHPj5+SEsLAzTp09Xuy937tzB7t270axZM1hZWaFz584qdYCyC/WIESPw+++/Y9asWWp1AODj44MNGzbgzJkzaus+qEPTPN6PfhX6IR+aWoaPHDlC+/fvZ9Ju27aN3NzcmLQ+Pj40fPhwwYaXRRQWFjLpRP67sDYIlcvlaut0vI1QA82KjlQqpSlTptDNmzcFtRkZGaStrU1Lly59p21BvOJVj6ZXvFKpFNbW1jh37pygVi6Xw9zcHG5ubqhdu7ag/vXr12jcuLFopBARKSfkcjm0tLSYPmOnT5/Gtm3b4OnpiZo1a2q8LdG59p68eaNCT08PHTt2hL+/v+DfaWtrY82aNcxmiyZNmohBV0SkHNHW1mb+jNWsWROXLl16p6D7LoiB9y3++OMPpKSkFP/f3t5eMF2siM6dO2PatGnlNTUREZFyon///qJz7WNjZmaGnJwcAMoGloaGhigsLCyhUSgU+C8t04iI/NeIiIgot8+4GHjfok2bNggNDS1u8Q4oSz6+XSaOiLB06VLk5uZ+6CmKiIiUM9HR0ZgzZw5zWyFNEQPvW/Tp0wcuLi4Qugmnra0NqVSKHj16IDQ0lFcbHR2NzMzMspymiIjIO3Du3Dn4+fnxauRyOWbNmoXt27fzmi3eBzHwvoWBgQHMzc3h5eUl2JTRxsYGEokE8fHxvLpatWqhb9++mD59umCQjo+PR25urkZfcTTRsp7Bi8bMzs7mnXORrqCgAMePH+d1xBVpIyIiBB18RWk39+/fR2RkpOCYRISYmBjeMTW5eimv418WpKWlITk5mdddBijnFRQUhOjoaEFteno6QkNDkZWVJbg/5XUVWATr8SzSJSYm8johi7S//fYbLl26xJv3CwDOzs4YO3as2t54ZYEYeFWgpaWFSZMmwcXFhVfXvXv3YnNA0ZqwKgwNDbF9+3YcO3YMV65c4R0zPDwcJiYm0NfX562OlJ2djb1792LcuHHo168fPDw8cOnSJURFRanUP378GMOHD8fSpUsxYcIEWFlZqT2xnD17FhMnToStrS2mTp2KV69eqdRJpVJMmzYNK1euhKWlJbKysmBkZKRSGxgYiBkzZmDs2LFYv3492rVrp1Inl8uxcuVKLFy4EGZmZjh37pzaZpPBwcEYMWIE5s+fj+HDh2PPnj0qdTk5OVi7di1GjhyJOXPmwNLSUmW32by8PBw5cgRDhgzBokWLsGHDBixbtqxUoJFKpdizZw/GjRuHTp064dSpU7h16xZev36tco6WlpZo164dpk+fjqioKJWBRS6XY+vWrWjQoAGMjY3h6+urNgBpa2vD3t4eurq6sLGxUVsDhOM4KBQK9O7dGwYGBrhx44ZKHQAYGRlh//79xU0x+ZbQrl+/jubNm8PIyAg7duxQG4jz8/OxatUq9OrVCx06dICvry/y8vJUaiUSCTw8PDBixAhMnz4dW7ZswerVq9XajH19fTFlyhTY2Nhg8eLFauvoEhEiIiIwbdo06OvrY/v27bzdJfz9/REeHg4bGxu1mjJBk6Tfyv7QxEAhk8loyJAhTOYFX19fmjp1qqAuLi6OHB0daf78+bzjxsfHk4ODAy1ZsoTMzMxo+/btKi2Rz549o5UrV5KVlRUdPXqU9u7dq7JBoK+vL40aNYp69+5Nv/32m9rEeIlEQsuWLaOGDRtS586deXu9xcXFUc+ePcnAwIB+/fVX3j5VJ0+epKpVq1KjRo14my0mJyfT4MGDycjIiKZOncrbS+348eP0xRdf0GeffUaXLl1Sq0tPT6eFCxfSwIEDydTUlKKiotTO9eHDh7R27VoyNTUlJycn8vPzo7CwMJX6hIQE2rx5M5mampKrqytt27aNLl68qHLchIQEmj17Ni1btozmzp1LlpaWFBMTo1IbEBBAdnZ2tHjxYho2bBht2rRJpa7Itrt27VoyNzen2bNnq7VXx8XF0ZIlS8jBwYGsrKzo77//VqkjIjp8+DA5OjqShYUFzZkzR+3rlZKSQlOmTKF169aRmZkZzZkzR6WJQyaT0aFDh8jGxoacnJzI1taW/vrrr1K68PBw+v3332nkyJG0ePFiun37Nj169KhUf0CFQkHbtm2jqlWrUsuWLSkqKkrtvhApTVBVq1alEydO8OqIlL3cBg8eLNhDUBUQm12WTeAlUnYyFep0W8Tq1avp6NGjTNqTJ0+SlZUVk4OosLCQfH19mcb9kCQmJlJ6ejpTg8/w8HCVtRHeJjc3l7cRZ2WHtYnimzpWl1lISIjKWhFvk5GRQcePH2ca8+nTp7wdkd9EqHvxuzSQVIVMJqPMzEym8RISEqht27a0YsUKio2N5dU+ffqUmjZtSmfPnn2neWkaeEXnGg8FBQUYPnw4Lly4IOjdlsvlsLKywrZt25jWhh49eoTc3Fz07t2beT4iIiLsREVFoX79+mqXqt5k5cqVOHHiBE6ePKm2PgkfonOtDKlSpQrMzMyYKtFra2tj+/btmDt3LlhOZl9++aUYdEVEypFmzZoxBd2iq9D79++/U9B9F8Qr3jeQyWS4fPlyiUpX+fn5GDlyJM6dO8dkP8zIyPhgtkMREZH3pygGvo+Fv1Jf8XIcN4rjuGCO4xQcx6ndCY7jhnIcF8px3AuO45aW1fZ1dHSwdu1anDx5svg5AwMD9OvXD5cvX2YaQwy6IiKVC47jVAbd8nSnVqjAC+ApACsAt9QJOI7TBrATwHcAPgdgw3FcmX0/6NSpE+bOnYukpKTi52bNmoV+/fqV0AUGBiI4OLisNisiIlLBWL16NW7fvl0uY1eowEtEz4iI32EAdAfwgogiiKgAwDEAFmU1h3nz5mHIkCEl8gerVq1aqmh569atMXDgQDg7Owv2cBJ7rYmIVAzkcrmgcw1QlonMzs4udcFVVlSowMtIIwBvZqrH/PNcmWBiYoKVK1di9erVvLqqVavC1tYWu3btKlHNTBUHDx6EjY0Nzp49i4KCgrKa6geD1akkk8kQFxfHpI2OjmY+FqLdWkSIjIwMwWWB2NhYWFlZqTUZFREcHAwXFxds2LChLKdYgg8eeDmOu8px3FMVD9arVlUr4GqPOMdx9hzHPeA47gFLYz5AWZGsUaNGgs0Tly9fjv3792P58uW8L7q9vT2kUikmTZrE20AzJycHM2fORLdu3WBjY1OqItqbyOVy+Pv7w9nZGefPn0d+fj7vHJKTk7F//36Eh4fD09MTa9aswbVr11Rqg4KCEBwcDHd3d0yaNAmjRo1SOReFQgFfX18EBgZi3rx5sLS0VFu7+PXr1whsWTvNAAAgAElEQVQLC8OhQ4dgZWWFNWvWqGzoSES4fv06Xr58iXXr1sHMzAx//fWXyjGTkpLg7e0NHx8f/PLLL7CwsEBgYKDKeXp4eODhw4e4efMmdu7cibVr16rcdlBQENavX4/4+Hg8f/4cfn5+pU4mBQUFuHv3LjZu3IhDhw6pPYEUbffnn3/G999/r9LZVkRcXBzs7e0xfPhwrFy5kve1PHHiBLp06YKePXsiICBArS4jIwMTJ06EiYkJxo0bx3sCvXHjBjp27IgOHTrg7NmzanVSqRQ///wzBg4cCGtra16H2/Pnz7F582bMmjULBw4cEAyMISEh+PXXXxETE4OIiAg8fvxYrQ385cuXePbsGY4ePYoxY8ZgxowZak/QRITz58/D3t4eW7ZswejRo9XOIT09HfPnz8e+ffugo6PDO9/3QpOk3w/1AHADQFc1v+sF4PIb/18GYBnLuJoYKFJTU8nc3JwpUXvHjh3022+/8WrS09PJx8eHLCws6MiRI2p1MpmMnJycyNTUlMzNzcnBwYFu3LhRysUVFBREM2fOpHr16tGQIUNo/PjxZG1tTTNmzCihKygoIGdnZ6pWrRrp6enRokWLyMXFhR48eFDK1CCVSunnn38mHR0d+uqrr8jV1ZVSUlJUzjMpKYm+++470tfXp6VLl1JISIjafbp48SLVrl2bOnfuTIcOHVJrpkhPTycbGxvS0tKiqVOn0u3bt0kul5fSyeVy2r17N9WsWZOMjY3J2dmZHj58qFJ779496tixIwGgoUOH0tatW+nq1aulGklKpVJydnamZs2aUbVq1cjBwYFWrFhBGzZsoMePH5fQRkZG0i+//EItWrSgvn37ko2NDVlbW9OYMWNU7tMvv/xCzZo1I1tbW7K0tKRRo0bR6dOnS2kTEhJo3LhxNHToUBo+fDhNnjyZDh06RHl5eaW0fn5+9PXXX5OtrS2NHj2a3NzcVDqu5HI5rVu3jkaPHk2mpqa0atUqtWaHyMhIGjRoENnZ2ZGVlRW5ubmpbGAqk8loy5Yt1K9fPxoxYgSNGzeO3NzcSjVmlcvldPr0aerSpQsNGTKELC0tydLSkn788ccSnyu5XE779u2jdu3akba2Nk2bNo0WLVpEjo6Opdr2KBQKOnjwIBkaGlLXrl3J1dWV12l29OhRGjlyJC1cuFDQnCOTycjKyorZNPIm+H9wrgkEXh0AEQCaA6gC4BGAL1jG1dS5tnHjRpUfkLdRKBQ0ffp0tZbRN5HJZLRmzRqaOnUqZWdnq9UV2S9fvnxJ69evJ39/f5W6wsJCSk5OFtxuTk4OBQcHq/wQFxEZGUlXr16lq1evUmBgoFqdQqEgd3d32rdvH504cUJtcCZS2lW3bNlCBw4coGvXrvHO8dq1a3T06FE6c+YMb3fh7OxsevToEfn5+akNzm8ik8koOjpa0F5aRHR0NJNVXKFQMB17IirhtpJKpWq7RxNRcS++pKQkOnr0qNr3SdF7JCsriw4fPkxbt25VO6ZMJiOFQkH37t2jH3/8UW1wKdrvnJwcOnToEC1cuFDtmEXdstPT0+nYsWNq36MKhaJEcORz4yUmJvK6HKOjo8nV1ZU8PDzozp07anVERKGhoWRoaEgDBgwQdE7evn2bpk2bxmQtVkWlDrwALKFcs5UCSCy6sgVgDODCGzpTAGEAXgJYwTq+poE3Pz+fhg4dyvQhlEqlNGzYMHr+/DnT2L6+vjRv3jyN5iMiIsKGQqGgI0eOMF+9WltbU+PGjXnrk/ChaeAVDRQCuLq6QiqV4ocffhDUJiUlYcKECXB3dxfzeUVEKglxcXGws7PDnj170LRp03cao1IbKD42kZGRpeqWjhs3DqdOnWLqNFGvXj1s3LhRvAsvIlKJ0NPTw4ULF9456L4LYuB9g/T0dIwcORL5+fnFz2lra2POnDnYsmUL0xgdOnRAs2bNymuKIiIiZUzt2rUFi2CVNWLgfYOOHTvi1q1bsLOzK2GKGDJkCL788ssS2pCQkA89PRERkQ9Ieno6Nm/eXC5ji4H3DbS0tLBjxw5kZWWV8G5zHAdzc/MS2hs3bsDS0lJtdwYREZGKyZvfaPk0kyZNwuDBg8tlDmLgfQtbW1tYWlpi3759vDo7Ozv4+flh5syZvInp2dnZWLZsGTw8PFQaBv7L/Jdu7IpUDC5duoQRI0YgIyNDrUYmk2Hy5MlYtGiR2hZV74sYeFVgZ2eHs2fPgs/ppq+vj8uXL6NZs2a89XqrV68OCwsLTJgwASNGjOAN0g8fPsS8efMwb948XL16VXCeubm5uHz5MlMtiMLCwmIXkEwmQ3h4OO7evatSW3SCyM3NxcWLFzF//ny1veIkEgmIlLVMlyxZgpEjR6rcx6I5ZmZmws3NDWPGjIGrq6va/QKU7qQtW7bAysoKaWlpKsfMyckBESEkJAR//vmn2r5rRd9MiAhJSUnw8fFROWZmZibTMlJOTg68vLyQmpqqViOVSvHs2TOcO3euRMU7Vdy4cQM7d+7EunXreLefnZ0NR0dHzJ49G+vWreN9P929exfjx4+HnZ0dbt1SW3cKBQUFcHJygp2dHRYsWMB7RRgWFgYnJyesWrVKZd+6IogIycnJ8PX1xcGDB3k/S0Xk5uYiKChIUAeguHdbdnY2PDw88P3338Pd3V2tPj09HT/88ANu3bqFM2fOqM06IiLMmjULY8eORd++fZnm8i6I6WRqCAwMxO7du/Hnn3/y6uRyOSZOnIipU6fim2++Uas7ffo00tLS4OnpiY0bN6JNmzYqdZcuXYK9vT3atWuH+vXrw9TUFIMHD0aNGjWKNTKZDLt378ZPP/2EwsJCmJqalhjjl19+wRdffAFA+ZVp//79+O2331CnTh00adIEOjo6aN68Obp06VLCPqlQKLBnzx5s2rQJ7dq1g56eHvr164dBgwbBxMSkxPJLdnY25s2bh+joaFSpUgXdunXDyJEj8fnnn5cqsXfv3j3MmjULLVu2hJaWFiwsLGBmZgZDQ8MSuvz8fPz888/w8fFBnTp1YGJiAnNzc/Tt27dEM0Miwrlz5zB//ny0atUKurq6aNOmDfr27YuvvvqqxIfq3r17cHJywp07d9CvXz9wHIe6deuidevWsLW1Rf369QEoM1pWrFiBU6dO4YsvvkCLFi2Kx2jevDmcnZ2Lt3306FEsWrQI8fHxsLS0LC4rWK1aNUyfPh09e/YEoAziu3btwubNm1GnTh20adMGDRo0QJs2bdC2bVsMGjSoeBu5ubnYsGEDtmzZggEDBqBu3bro06cP+vXrh+bNm5c4pomJiZgxYwbCwsJgYmKCfv36wdzcHCYmJqXeTzdv3sTkyZPRrl07NGzYEKNGjUL//v1L2WElEglWrFiBc+fOoXXr1ujTpw+sra3RsmXLUmNeuHABDg4OaNq0KYyNjfH1119j4MCBaNGiRfE8iQgXLlzA+vXr4e/vj++++w4GBgb47LPP0KpVKwwYMAANGzYs3p+ffvoJR44cQdOmTUtcZRoZGZWyjF+4cAFLlixB69atoauriyFDhsDMzAx169YtNVc/Pz8kJSVh586dWLNmDbp25c/4WrlyJRo3bgx7e3te3dtomk720U0TH/KhqYFizpw5gu4YIqXRwsLCggICAgS1cXFxZGNjQ2vXrlXZHJBI6XLKzc2lpKQkcnFxIVtbW1qzZk0pXXZ2Nnl4eJQweCgUihJOroSEBLp48SJt376d15WjUCjo0KFD9MMPP9DEiRN5XV5JSUk0efJkGjRoEC1btozXVn3+/Hnq378/DRw4kO7du6dWl5OTQ0uWLKERI0bQjz/+yNvo8sGDB7Rs2TKaOXMmr1uwsLCQ7ty5Qy4uLvTrr78KupeysrLo0qVLdOXKFV4dkdI1dv78+WJLrUKhoOzs7GLX2Zvk5ubSo0ePSKFQUHx8PF2/fp1cXFxUjhsZGUkymYySkpLo1KlTNHfuXJUWc4VCQYmJiVRQUEDe3t40Z84cMjMzU2mfTU9PJ7lcTlFRUbRp0yYaNmwYbdu2TeX2U1JSSCaT0c2bN2nWrFlkZmZG4eHhKvcpOTmZsrOz6fz58zRv3jxatGiRyjHv3btHcrmc8vLy6MmTJ3Ty5El68eJFCY1UKiUfHx9Bp6inpycNHz6cRowYIWiOCAwMJENDQxo2bJigXTguLo527txJjo6OvDp1oDI718r7oWngzcjIYO40nJmZSUOHDqWwsDCmsT08PGjo0KGlagaIiLwrmjSUVHWCUIVcLuc9CVZUUlJSaOLEibRp0ybBz6RCoaCePXtS7969mT7rqtA08IpLDQIcOXIESUlJmDt3rqA2MTEREydOxIEDB2BsbCyoT09Ph5GREbS1tTWak4iICD9ExNzK58qVK3BwcMDOnTtLLP9oguhcew9U3VSwsbHBjRs3mOrM1q9fH7t378bRo0eZtvfJJ5+IQVdEpBzQpH9alSpV8OTJk3cOuu+CGHjfYN26daXqm3IcBycnJyxbtoxpjObNm2PBggXlMT0REZFyoH///tDT0/ug2xQD7xt07twZffr0gYeHR4nn27Ztiy+++KI4JUYikTB3ZRAREam83L17l6lOi6aIgfcNBg8ejO7du5dodFnE4sWLi9NVOI7D6NGj4eXlhf/SGrmIyH+JK1euYOvWrdDV1S3zscXA+wbVqlXD9evX4eXlhcjISLU6PT09WFpaYvDgwTh8+DDvmAEBAbh69Wq5nDVFREQ0QyaT4ffff8exY8d4dV5eXvjrr7/g5uZWIoe8rBAD71toaWlh69atmDt3Lu9ygo2NDbZt24ZDhw4hISFBra59+/bYuHEjatasyev0kcvl8Pb2xpkzZ3D8+HFeS+ObZGdnM+mA0hZdIhJcMlEoFHj06BH+/vtvwbGDg4Ph7OwMHx8fXm1OTg48PT0xb948wW8MSUlJOHLkCLy8vHh1gNJx5+3trbLrs6rtqNp3IipVGlQdLK8RESEhIQHx8fFqNUWuwqCgIHh7e/NuPyQkBKdPn8bx48d5+63JZDIcPXoUe/fuxV9//YWcnBy12oiICLi6usLFxUVtH76ifXn27Bl8fX1x6dIlJCYmqtUWIZVKER4eLqgrgq/PIN82rl27JljQ5tatWzA3N0fdunUxZswYtbryDrqAso2OyFs0a9YM5ubm+OOPP+Dg4KBSo6WlhdmzZ2PYsGH4/vvvsX//fpUpZFWqVMGpU6fw888/w83NDSEhIZg/f36pxXxtbW3Ur18fK1euxOPHjzF69Gh8++23GDhwIOrVq1dCGxsbi40bN+LKlSuoXr06zM3Noa+vDwMDA+jr66NVq1bFdkciwpUrV7Bp0yaYmpoiNjYWUVFRxcFp/fr1+Oyzz4rHvnDhAh4+fAh9fX3cv38fhYWFaN++PQYOHFhiDv9r786jmrq2P4B/jyMq9FVE6lAQpE6ooIBWAQewQJ9VUUTk4VB9PqeqtE9ba9HWtlafdRbqVEUtVSzq4jk8UYsFRBc4IDKIA6CgjApEQIYAIfv3R4YfMWTAQqB6PmtlEZLDzb43N5ube885u7q6Gtu2bcOQIUNw7tw5ZGZmYuDAgZgwYYJ85JZMWloaLl26BMYY/vjjD7Ru3RouLi748ssvFa4+19bWYvfu3TAyMkJiYiLu3r0LY2NjuLi4wMbGRmGZ8fHx2L9/P2xsbHD9+nUUFhbCwMAAw4YNg4ODg7y3SFZWFnbu3IkXL17AxMQEDx8+lA+JdnZ2xpIlSwBIhsIGBwfj2LFj8Pb2BiD5ByG7rVu3Dr1794ZAIMCGDRtw7tw5iMVizJgxAwYGBtDX14e+vj5MTEzg6OgIAIiMjMSSJUtw7949eHh4wMHBARYWFnjvvffQu3dvdOjQQb4+p06dwrfffgvGGGbPno0RI0ZgxIgRePfddxW2kbm5OYKDg/Hjjz9i1KhRGD9+PNzc3DBo0CCFdm3atMHIkSPxr3/9C4mJiYiJicFHH30EV1dXpRGD5ubmiIqKwvLly9GvXz989NFHmDBhAoYOHao0WVR5eTk+//xz3LhxAz4+Phg7dixGjx4NCwsLhbaFhYVYtWoVDh06BHNzc8yYMQMWFhbo3bs3evfujW7dusmnYkxLS0NISAhCQkLg7OyMLl26oKKiAhUVFaisrMTixYsV3v/4+HhERkaiXbt2iI6OBhHBwcEBH374IV5WUVGB4uJirF69Gl27dsXx48dhYGCg1E7m0qVLOHDgQJMmXQB8AIUqYrGYpk6dWu+onZc9evSIXF1dKTs7W+MyT58+Ta6urnT27Nl6O7zX1tZSdHQ0PX36lIKDg2nevHk0ceJE+uKLLxRqb+Xm5lJAQABt2rSJbt++TbGxsRQZGUnnz5+nuLg4+euFhobS9OnTydramvbt20fp6en1jpgTCAQ0b9486tixI1lYWNDly5dVdrKPi4sjKysrAkC+vr6UkZFRbzuRSETbtm0jPT09MjQ0pDNnzqis+Xb79m2ys7MjAOTl5UUpKSn1bh+RSEQ7d+4kMzMzevvtt+m3335TWffs8ePHtHTpUhozZgx5e3vTrVu3qLi4WKmdWCymxMRE2rlzJ02ZMoVOnTpFycnJlJGRQQUFBVRZWakQS1lZGR09epR++OEHSkhIoKtXr9L58+fpxIkTSiPphEIhBQYGytsdPnyY1qxZQz4+PkqjIouLi+nixYv08OFDOnr0KC1btowmT55M8+fPV4r5zp079PTpU7p8+TL5+fnRxIkTac6cOUqjA8ViMcXHx1N+fj4FBgbS9OnTydPTk3bv3i2vmSaTk5NDT548oRs3btA333xDEydOpPnz51N0dLTSMsPCwqiwsJBOnz5NK1asIHd3d/Lx8aGcnByFttevX6ejR49SUlISnTp1irZt20ZLly6lb7/9Vt4mMzOTDh48SDNnzqRffvmFkpKSKC0tjXJyckggECjUgVu+fDm1atWK3nnnHbpw4YLaGoJxcXFkY2NDU6ZMUVuMlUjyuQsPDycvLy+l7aIN8JFrjZN4iYiysrJo4sSJGospEv1/hdasrCyNbSsrK+mHH34gT09PevDggcb2YrGY7t27p7MRRDU1NWpHQT1//pwKCgro2bNnVFpaqrKdUCiknJwcys3NpdzcXI2jgkQiEVVUVNRb2VaVhozW+qvSdjTVs2fP1BYJlamoqKBz585ptV/n5uZSQkKCVq//4sULlcPgG4tsHykpKVEbf2pqKpmZmdEHH3ygMekKhULy8vJ65aRLxBNvoyZeIqLDhw/T9u3btWqbmZlJLi4uWu38RJKKqZGRkQ2OieM49TIzMzXOzyCzY8cO6tix4ytXGCZqeOLlF9c0mD17NmJjY5Gamqqxba9evRAYGIj58+drNdmyiYkJxo4d2whRchxXV69evbQ6R1tZWQmBQICHDx/C09NTB5FJ8LkapAQCAY4fP46FCxcqDTfMzc3FwoULcerUKa2G+AqFQujp6TVKzBzHtXx8roZXZGhoiNDQUHh7eytViujRowe8vb1VTrL9Mp50Oe71QEQN6rKpLZ5465gxYwauXr2KrKwsped8fHwwZ84cAJJ+fkFBQa/U55DjuL8GkUiETz75RKs+5A3FE28d06ZNQ3h4ONauXavUuZ4xJu936eTkhO3bt2Pw4MFqO8aLxWI8fPiQz+vAcS3IixcvcOjQIbVthEIhZs+eDVdXV3h4eDR6DDzx1tGxY0dYWlpiwoQJ2Lp1q8p2bdq0we7du9G3b1/s27cPqs6Tt2rVCrGxsejatSumTJmidvQQIBm1IxAItJqCUtdKS0tRXV2tsd2zZ89U1nKrSyQSISYmpt5vFy+rrq5GcnKyVnECqHfk2l8BEUEoFKptU1NTg+LiYuTk5MjrjtWnpKQEGRkZSE1NRXp6utrXLCgoQFZWFtLT07UajScWi/H8+XOV+31zICKkpaXVW0dPpqamBrt27YK3tzd69eqlsl1JSQmmT5+OBQsWYMqUKU0RLh+5Vp+PP/4Y8+bNw/Xr1/H+++/X22bkyJE4ffo09u3bhyVLliAgIKDeC28zZ84EAAQGBmLu3LlYs2YNrK2tldrl5eXh3//+N0JCQjBkyBCsXbsWo0ePhqGhoVLbrKwshIeH49KlS5g8eTL69u2L9u3bo3379tDT00P79u1haGgov0iYlpaG+Ph42NvbIzMzU+Fma2uLTz75RL7s5ORk9OzZE0+fPkVsbCyuXbuGgoICGBgYYOPGjfLReUSE+/fvw8zMDFeuXEF4eDgePHgAY2NjTJ48WSFeWUFIWXHOyMhIVFZWYvjw4Zg9e7ZC24SEBFhaWuLmzZuIiopCfHw8WrVqBUdHRwwePFhhe+Xn58PY2Bi3bt1CfHw8UlJSIBKJMG/ePEyYMAGA5Kr1kSNH8Pe//x2FhYUK656bm4vg4GC0adMGFRUV+N///ofU1FR4eHigvLxcfquoqECvXr0wfPhwpKenIyIiApGRkbC3t4eDgwP09fXRqVMn+a1t27YQCoXYsWMHIiIikJ2djU2bNqFPnz7o1auX0jWAlJQULF26FNevX4erqys+++wzDB06VKHOHiD5h3LgwAGsXr0apaWl2LJlC9zc3NC/f/9655/dvHkz9uzZg2HDhuGrr76Cs7Oz0jIZY0hKSsLixYuRlpYGPz8/TJ06FUOGDJGPLJMRCARYuXKlfDTa5s2bMWLECHnttLouX76M/fv3486dO5gwYQLc3d1hZmYGIyMjpVgrKytx4cIFxMTEYO7cuaisrIRQKERlZSUqKyvRpUsX2Nvby9tnZ2dDT08PJSUliIqKQnR0NIqLi9GnTx8sW7ZM6TNDRDh58iQCAwPx8ccf4+zZs0rrJvP06VPMmTMHGzZswNChQ+tt0xh4rwYVXrx4AU9PT4SEhKisSCoTGhqKkydP4uDBgyovrFVXV6OoqAjr169HeXk5vvnmG5ibmyu1i4qKQnZ2NmpqahAdHQ2BQABTU1OMHTsWY8eORZcuXZCUlIQTJ07gxIkTmDRpEt566y1UVVVBKBTKfwYEBKCqqgpr165FQECAfHy6mZmZwk22bvn5+fj6669x8OBBjBs3Dra2tvJhq7KCkDKpqalYvHgxUlNTMXLkSIwaNQqurq7o27evwoeKiPDrr79ixYoVGDhwIAYOHAg3Nzc4OTkpDdt8+PAhPvvsM1y5cgUffPABhg8fDicnJwwdOlShMKNAIMCmTZvg7+8Pc3NzODk5wdbWFjY2NrC0tJTPJEVEOHPmDHx9ffHkyRNMnz4dgwYNUlj37t27o3Xr1igpKcHGjRtx5swZCIVCzJgxQ55EO3bsiE6dOqF///6wtrZGfn4+Tpw4geDgYNjY2KBnz54oKyuTJ+mysjJs3boVPXv2RFVVFY4dO4awsDC4uLggIyMDjx8/lh/V2traws/PD4Dkq+3evXuhr6+P8vJy3L59G6WlpejYsSOsra3h4+ODnj17ApAMxw0ODkbfvn0RFRWF+/fvQ19fH46OjnByckK/fv3k2ysmJgYlJSXIy8tDREQEysrKYGdnBzc3N9ja2soTkFAoxKFDh2BlZYXw8HAkJCRAX18fTk5OcHFxgampqXyZkZGRSExMROfOnXHt2jXk5eWhU6dOGDZsGNzc3DBgwAAAkkKTO3fuxJAhQ1BbW4vHjx+jsLBQfqS8Z88eGBgYwN/fH//973+RmZmJBQsWyIe+d+jQAR06dICZmRnGjBmD6upq7NixA99//z3s7OwwZMgQODk5YfTo0ejcubPSZykxMRERERH4448/4OLigkWLFqmdd/fRo0dYvHgxdu3aVW/hUHUa2quBJ1414uPjsX37dgQFBWmc0f7y5cvYunUrgoKCNCbqtLQ0rFu3DoaGhvIx5OpkZmbi8uXLMDExgbOzs/xxIsmELpqmrSsvL0d6ejqsrKxUrsejR4/w5MkTFBUVoU+fPrCysqq3XVVVFU6cOIGSkhIQEebPn69yZ7579y4iIiJQXV0NKysrpfkeZMRisbxUellZGTw9Pes90gckY+9zcnKQl5cHsVissR90bW0tMjMzYWBgoDTnRX0EAoHK136ZSCRSqtbbENr8vazkuYWFBbp166ayXWlpKa5evYr79+9j+fLlKtvV1tYiLi4OFy9ehJ2dnVKF6rpkR5SRkZHYsmWL2ljLysoQFxeHtm3bwsHBQeE5IuUyPLK8U/fxkpISvPXWWyr30bt37yIxMRHl5eUYMGCA0uvUlZSUBGdnZ7Rt2xaJiYlq3/uSkhKcPXsWR48excGDB+s9gteEVxlu5JFrO3bsoL1792rVNjExkdzc3DTO2SBz8+ZNmj9/vlZDNzmO005NTQ0dOnSIIiMj1Q5pl/n000+pdevWWg+Nrg/4kOHGTbxisZimT59OSUlJWrXPyMggV1dXunfvXoNfi+M43UpKSiI3NzeKiIj4U/N+NDTx8lMNUjdv3oSlpSU6deqk9FxRURF8fHwQGhpa7/MvKywsRFFRkcK5No7jWp6SkhKlC46vgo9ce0V6enoYOnQorl27pvRcly5dsGbNGrXnzuoyMjLiSZfj/gIaI+m+Cp54pQYPHgxjY2MsWbIEVVVVSs+PGjUKs2bNglgsxpkzZ9T2oeQ47vXw+PFjLF68GI19ZqBFJV7G2DTGWApjTMwYU3nYzhjLZIwlM8YSGGPad1PQ4JdffoGLiwvOnDlT7/OOjo5o1aoVqqqqYGJigu+///4v21mf495ksh406sTGxmLRokVYvXq1xl5NDdWiEi+AOwA8AERr0daJiIY05LyKJhYWFli/fj1+++03tSOlpk2bBhcXF1y5ckVtPanc3FxMnToVM2fO1Fhc73UiEom0HiatzfSZMm/S9Qju1WjaR4qKirB69WosWLAAhYWFKtsFBwcjICAAJ0+exLvvvtvYYbbMXg0AogDYqXk+E4BRQ5erba+GoqIicnFxoaKiIpVtKioq6NmzZ8KQ0TwAAA04SURBVDRx4kS6cOGCynbZ2dlkaWlJTk5OdPjwYZXVBMLCwmj06NE0YMAA2rdvn8qyO0SSMiVxcXF0+PBhrdbnxYsXlJeXJ/+9pKSEkpOTKSwsjPbv369wNbduKZXCwkK6ePEirV+/nry8vGj//v0Ky5VVxKiurqbY2FjauHEjTZ06lTw8PCg/P1+pHRHR06dPKSQkhBYuXEiTJk2izZs3KyxTVsGgtraWUlJS6Oeff6Y5c+aQu7s7xcfHy9tVVlbKJ7quqqqixMRECgoKos8//5yWLVsmbycWiyk9PV1h+RkZGXT58mU6cuSIQm+VlJQUSk1NVYhHLBZTVVWVvGRQTU0N3bhxgwICAtRWhnjy5AkdO3aM/Pz8KDAwUGWXwdraWjp48CDNmjWLnJ2d6caNGyqvrufk5NCMGTNo0KBB5OHhQbm5uSpfPzQ0lAYMGEBdu3alPXv2UHl5eb3tiouLadGiRWRubk62trZ0//59la9/9epV8vT0JAcHB/L391e5jwqFQoqMjKSdO3fSp59+qlVFkQcPHtDt27eVHheLxSQUChViqlvloqCggE6dOkVffPEFTZo0iTZt2lTv8gsLC8nPz488PDzo6tWrKuOora2lr7/+mlatWtWgbp54HXo1MMaiAHxORPWeRmCMZQB4DoAA7CMilfM1MsYWAFgAAKampraPHz/WKobk5GR89913CAkJUTsHb3V1NZYtWwZLS0v4+vrW+5WkqKgIHTp0kBf08/LywqxZs5QGPty9excrV66Eo6MjEhISIBQKYWFhAXt7e9jb26N79+4oLy/Hhg0b8PPPP6Nt27YKQyllr01EmDt3Lvr3749du3YhMDAQNjY28mHEBgYGMDU1hYmJCUxMTODi4oKKigps27YNoaGhsLS0RGVlJQwNDWFnZwc7OzsMHjxYPiovPz8fq1atQqtWreRTaNra2mLMmDGws7NTmID6999/x6ZNmzBgwAA8fvwYxsbGGDt2LJycnOQjsQCgoKAAa9euRU5ODtq3bw+RSIR+/frBwcEBI0eORJcuXQBIBnAEBgbiP//5D8aMGYOKigq0a9cO/fv3h5WVFaytrWFubo7q6mocO3YM/v7+EAgEsLW1BWMMbdu2Rc+ePeXr//777yM3Nxe+vr64du2afHRg3e3Zrl076OvrIyAgAAEBAdi7dy8KCwsxbtw4hfdP1n7ixImYNm0agoODsXPnThgZGcHQ0FDygWMM+vr6MDU1hampKby9vaGnp4dff/0Vu3btwujRo5GZmQkiQo8ePeTrNHjwYHTq1AmVlZXYunUrsrOz5e/F22+/jZEjR8LBwQGWlpby0WgZGRn4+uuvMXz4cMTExKCqqgrW1tZwcnLCiBEj5ANfxGIxjhw5gpiYGOjr6yMtLQ0GBgYYNWoUxowZg379+snXLSEhAX5+frC0tEROTg6EQiGMjIzk+8mgQYPQrl07xMTEICAgALGxsbC3t5fP89GqVSt069YNJiYm8hGMvr6+CA8Px/Dhw2FiYqL0+Wnfvj12794NAwMDBAcHY/fu3Rg4cCAKCgpgZGQEe3t7ODo6ok+fPgqfv0ePHuFvf/sbtm3bhnv37mH58uXyQqT1qaiowMKFCzFu3Dj5TITaavEDKABcguSUwss39zptoqD+iLeH9KcxgEQAo7V57Yb24w0JCaGvvvpKYzuxWExbt26lxYsXa6w5VVNTQ0FBQeTm5kb79u2rtzyJ7L+77GgtKCiIFi1aRAsWLFBYjqyopSoCgYCio6Ppp59+oitXrqhsV11dTT/99BPNnz+fvLy8qLCwUGXb1NRU8vLyIgcHB/Lz81O5vmKxmA4cOECOjo5ka2tL169fV7nMvLw8WrJkCY0fP56WLFmitu5VeHg4rV27lv75z3/SjRs31K5TcnIyBQcHk7+/v8Y+mlVVVXTlyhVKTk5W2062bomJiVr1+xSLxSQQCBQeKy0tpZSUFDp//rzCN4y621IsFlNOTg6FhYXRxo0byd/fX2m5Ms+fP6ewsDDy8/OjyZMnq+xzLhKJ6NatW7RlyxaaNm0a+fr6qlxmcXExnTt3jlauXEnu7u70+++/Ky1L5tmzZ3T+/Hlat24dzZ49W2E5L39rFIlElJOTQ9euXaP79+8TkeQoMyEhgW7evFlv3LI2e/bsIXd3d3J2dpb/rSq3b9+mbt260Ycffqh235eta25uLo0fP16psKe28DoMoNCUeF9q+y0kR8eNnniJiL788ks6fvy4Vm3PnTtH7u7uahOXjEgkoqNHj5Kbm5vGBMpxTeF1LRT66NEj8vHxoZUrV1J4eLjG9tOnTydHR0eFU1IN9donXgCdABjUuR8D4ENtlvsqiVckEpGHh4fWI9dSUlLI1dWVUlJStF7+q1Y25TjuzwkNDaXu3bvT+vXr/1SF5IYm3hbVq4ExNoUxlg1gJIBzjLGL0sd7MMbCpM3eAXCVMZYI4AaAc0R0oaliat26Nfbv348VK1aonetTxtLSEkePHsXp06e1Xr66GZM4jmsaRITOnTsjMzMTfn5+Giebakwt8uJaU2no7GR1aXuxjeO4Nw8fMtxEBg8ejH/84x+4d+9ec4fCcdxfHK9A0QBTp05t7hA4jnsN8CNejuM4HeOJl+M4Tsd44uU4jtMxnng5juN0jCdejuM4HXuj+vEyxgoAaDdLjoQRANVzx7U8PN6mxeNtWn/leHsRkfpy4XW8UYm3oRhjcQ3pFN3ceLxNi8fbtN6kePmpBo7jOB3jiZfjOE7HeOJVT+UE6y0Uj7dp8Xib1hsTLz/Hy3Ecp2P8iJfjOE7HeOKto7nLyzdUA+L9kDH2gDGWzhhbpcsYX4rDkDEWzhhLk/7srKJdrXTbJjDGzjRDnGq3F2OsPWMsRPr8dcaYma5jfCkeTfHOYYwV1Nmm/2qOOKWxHGSMPWOM3VHxPGOM+UvXJYkxZqPrGF+KR1O8YxljJXW27TdaLbghs6a/7jcAAwD0QxNVOW6OeAG0BvAQQG8A7SCpUWfZTPFuArBKen8VgB9VtCtrxm2qcXsB+ATAXul9bwAhLTzeOQB+aq4YX4plNAAbAHdUPD8ewHkADMAIANdbeLxjAfyvocvlR7x1ENE9InrQ3HFoS8t4hwNIJ6JHRFQN4DcA7k0fXb3cAfwivf8LgMnNFIc62myvuutxEsA4Vl95ad1oSe+vRkQUDUBdKRd3AEEkcQ3A24yx7rqJTpkW8b4SnnhfDQH4nTF2S1o+viXrCSCrzu/Z0seawztElAcA0p/GKtrpMcbiGGPXGGO6Ts7abC95GyISASgB0EUn0SnT9v2dKv3qfpIxplxDveVoSfurtkYyxhIZY+cZYwO1+YM3biJ0xtglAN3qeWo1EWlXKA1wIKJcxpgxgHDG2H3pf8ZG1wjx1nck1mRdWdTF24DFmEq3b28AEYyxZCJ62DgRaqTN9tLpNtVAm1jOAjhGRFWMsUWQHK07N3lkr6YlbVttxEMyXLiMMTYewCkAfTT90RuXeInog0ZYRq705zPG2H8h+brXJIm3EeLNBlD3COddALl/cpkqqYuXMfaUMdadiPKkXx+fqViGbPs+YoxFARgKyXlMXdBme8naZDPG2gD4G5rg66iWNMZLREV1ft0P4EcdxPWqdLq//llEVFrnfhhjbDdjzIiI1M45wU81NBBjrBNjzEB2H4ArgHqveLYQNwH0YYyZM8baQXIxSOc9BaTOAPhYev9jAEpH7Iyxzoyx9tL7RgAcANzVWYTaba+66+EJIIKkV1qagcZ4XzpHOglASy4ceAbAbGnvhhEASmSnp1oixlg32fl9xthwSHJqkfq/Au/V8NIVyimQ/MetAvAUwEXp4z0AhEnv94bkynEigBRIvvK32Hilv48HkArJUWNzxtsFwB8A0qQ/DaWP2wE4IL1vDyBZun2TAcxrhjiVtheA7wFMkt7XA3ACQDqAGwB6N/N+qyne/0j31UQAkQD6N2OsxwDkAaiR7rvzACwCsEj6PAOwS7ouyVDTu6iFxLu0zra9BsBem+XykWscx3E6xk81cBzH6RhPvBzHcTrGEy/HcZyO8cTLcRynYzzxchzH6RhPvBzHcTrGEy/HcZyO8cTLcRynYzzxchwAxth7jLEaxth3Lz2+hzH2Qt1E8xzXUDzxchwAIkoHcADAv6VzREBaTeCfAKYQUbNVGuFeP3zIMMdJMca6QTJHwG4A9yGpIvsPIjrerIFxr503blpIjlOFiPIZYzsArIDks+HLky7XFPipBo5TlAagPYBYItrV3MFwryeeeDlOijHmDGAfgFgADowx62YOiXtN8cTLcQCkZcRPQXKBbSyAJwA2NGdM3OuLJ17ujccYew+SkuK/A1hGkmq93wEYzxgb3azBca8l3quBe6NJezLEQHKE60ZEVdLHW0NS0uk5Edk3Y4jca4gnXo7jOB3jpxo4juN0jCdejuM4HeOJl+M4Tsd44uU4jtMxnng5juN0jCdejuM4HeOJl+M4Tsd44uU4jtMxnng5juN07P8A+tQvlpvB7pYAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "### Set the values of xdot and ydot all over the mesh\n", "xdot = x_dot(X,Y, b=0.)\n", "ydot = y_dot(X,Y, b=0.)\n", "\n", "### Create figure\n", "fig = figure(figsize=(5,5))\n", "ax = subplot(1,1,1)\n", "\n", "### Plot vector field (one vector every 40 points)\n", "ax.quiver(X[::40, ::40], Y[::40, ::40], xdot[::40, ::40], ydot[::40, ::40],\n", " pivot='mid', units='inches')\n", "\n", "### Calculate one trayectory starting at the position (xmap, ymap).\n", "xmap = -1\n", "ymap = -0.5\n", "\n", "### Time-step size\n", "dt = 0.01\n", "\n", "### Damping rate (b=0 means no damping)\n", "b = 0.\n", "\n", "### Create lists to store the solutions x(t) and y(t)\n", "sol_x = []\n", "sol_y = []\n", "\n", "### Fill the list with the solutions using one of the integrators (euler or rk4)\n", "npoints = 1000\n", "\n", "for n in range(npoints):\n", " \n", " ### Store initial condition\n", " sol_x.append(xmap)\n", " sol_y.append(ymap)\n", " \n", " ### Integrate\n", " xmap1, ymap1 = intgrator_rk4(xmap, ymap, dt, b)\n", " \n", " ### Update for the next integration\n", " xmap = xmap1\n", " ymap = ymap1\n", "\n", "### Plot the solutions x(t) and y(t) in the yx-plane.\n", "ax.plot(sol_x,sol_y,ls='-', color='g')\n", "\n", "### Plot fixed points\n", "ax.plot(-1,0,'ro')\n", "ax.plot(0,0 ,'ro')\n", "ax.plot(1,0 ,'ro')\n", "\n", "ax.set_xlabel(r'$x$', fontsize=16)\n", "ax.set_ylabel(r'$y$', fontsize=16)" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 2", "language": "python", "name": "python2" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 2 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython2", "version": "2.7.15" } }, "nbformat": 4, "nbformat_minor": 2 }