{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Effective potential for null geodesics in Schwarzschild spactime\n", "\n", "This Jupyter/SageMath notebook is relative to the lectures\n", "[Geometry and physics of black holes](https://relativite.obspm.fr/blackholes/).\n", "\n", "Click [here](https://raw.githubusercontent.com/egourgoulhon/BHLectures/master/sage/ges_effective_potential_null.ipynb) to download the notebook file (ipynb format). To run it, you must start SageMath with `sage -n jupyter`" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "'SageMath version 9.0.beta9, Release Date: 2019-12-08'" ] }, "execution_count": 1, "metadata": {}, "output_type": "execute_result" } ], "source": [ "version()" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "%display latex" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "text/html": [ "" ], "text/plain": [ "-(2/r - 1)/r^2" ] }, "execution_count": 3, "metadata": {}, "output_type": "execute_result" } ], "source": [ "var('r')\n", "U(r) = (1 - 2/r)/r^2\n", "U(r)" ] }, { "cell_type": "code", "execution_count": 4, "metadata": { "scrolled": false }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAnYAAAHWCAYAAAD6oMSKAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAAPYQAAD2EBqD+naQAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi40LCBodHRwOi8vbWF0cGxvdGxpYi5vcmcv7US4rQAAIABJREFUeJzs3Xd4FOX6xvFvGgECAYRIEwiIlNCblIOoB6WoPztgA1FQOGJBFBXFI0oLxRCRACLSq4AgQgSCUoIEDL0EpAkBIZBQEhIggWR/f8whRSkJzO7sbu7PdeVyZjM780Ru4uO7M+/rYbPZbIiIiIiIy/O0ugARERERMYcaOxERERE3ocZORERExE2osRMRERFxE2rsRERERNyEGjsRERERN6HGTkRERMRNqLETERERcRNq7ERERETchBo7ERERETfhlo3d2LFjqVy5MgULFqRRo0ZERkbe8PgFCxYQFBSEr68vQUFBLFy4MMf3u3btioeHR46vZs2a2fNHEBEREckzt2vs5s6dS+/evfnkk0/YunUr9913H+3btyc2Nvaax0dFRdGpUyc6d+7M9u3b6dy5Mx07dmTjxo05jmvXrh0nTpzI/AoPD3fEjyMiIiKSax42m81mdRFmatq0KQ0bNmTcuHGZr9WsWZMnn3ySoUOH/uP4Tp06kZSUxM8//5z5Wrt27ShRogSzZ88GjBG7c+fOsWjRIvv/ACIiIiK3yNvqAsyUlpbG5s2b+eijj3K83qZNG9avX3/N90RFRfHuu+/meK1t27aEhobmeG316tXceeedFC9enPvvv5/Bgwdz5513XreW1NRUUlNTM/czMjI4c+YMJUuWxMPDI68/moiIiLghm83G+fPnKVeuHJ6et/9Bqls1dgkJCaSnp1O6dOkcr5cuXZq4uLhrvicuLu6mx7dv354OHTpQqVIl/vzzTz799FP+/e9/s3nzZnx9fa953qFDh/L555/f5k8kIiIi+cHRo0e56667bvs8btXYXfX3ETGbzXbDUbKbHd+pU6fM7dq1a9O4cWMqVarE0qVLefrpp695zn79+tGnT5/M/cTERCpWrMjRo0fx9/fP08/jbFJSUihXrhwAx48fx8/Pz+KKXMvBgwe5++67rS5DLKYcCCgHAklJSVSoUIGiRYuacj63auxKlSqFl5fXP0bnTp069Y9RuavKlCmTp+MBypYtS6VKldi/f/91j/H19b3maJ6/v7/LN3ZeXl6Z2/7+/mrs8mjlypU0aNDA6jLEYsqBgHIgWcy6TcutnootUKAAjRo1IiIiIsfrERERtGjR4prvad68+T+OX7FixXWPBzh9+jRHjx6lbNmyt1+0C/L19WXJkiUsWbLkuh9Fy/W99dZbVpcgTkA5EFAOxHxu1dgB9OnTh4kTJzJp0iT27NnDu+++S2xsLD179gSgS5cu9OvXL/P4d955hxUrVjBs2DD27t3LsGHDWLlyJb179wYgOTmZ999/n6ioKA4fPszq1av5v//7P0qVKsVTTz1lyc9oNW9vbx599FEeffRRvL3datDXIUaOHGl1CeIElAMB5UDM53bTnYAxQfHw4cM5ceIEtWvXZtSoUbRq1QqABx54gMDAQKZMmZJ5/Pz58+nfvz+HDh3i7rvvZvDgwZn3zl28eJEnn3ySrVu3cu7cOcqWLcuDDz7IwIEDqVChQq5rSkpKolixYiQmJrr8R7EiIiJiDrP7A7ds7JyROzV2ly9fZubMmQC8+OKL+Pj4WFyRaxk0aBD9+/e3ugyxmHIgoByIGjuX5U6NXUpKCkWKFAGMj6r18ETexMXFUaZMGavLEIspBwLKgZjfH7jdPXYizu6nn36yugRxAsqBgHIg5lNjJ+Jg9evXt7oEcQLKgYByIOZTYyfiYKdPn7a6BHECyoGAciDmU2Mn4mAXL160ugRxAsqBgHIg5lNjJ+Jg99xzj9UliBNQDgSUAzGfGjs7CwsLIygoiCZNmlhdijiJVatWWV2COAHlQEA5EPNpuhMHcafpTq5cucLChQsBeOqpp7T6RB4lJiZSrFgxq8sQiykHAsqBaLoTcQLe3t506NCBDh06qKm7BV9//bXVJYgTUA4ElAMxn0bsHMSdRuxERETEHBqxE8tduXKFefPmMW/ePK5cuWJ1OS5n0KBBVpcgTkA5EFAOxHwasXMQdxqx05Jit+fcuXMUL17c6jLEYsqBgHIgGrETcXkzZsywugRxAsqBgHIg5lNjJ+Jg999/v9UliBNQDgSUAzGfGjsRBzt48KDVJYgTUA4ElAMxnxo7EQcrUKCA1SWIE1AOBJQDMZ8aOxEHK126tNUliBNQDgSUAzGfGjsRB9u8ebPVJYgTUA4ElAMxn5YNkDwrUKAAkydPztyWvHniiSesLkGcgHIgoByI+TRiZ2dhYWEEBQXRpEkTq0sxjY+PD127dqVr1674+PhYXY7L+fbbb60uQZyAciCgHIj5NEGxg7jTBMUiIiJiDk1QLJa7cuUKS5cuZenSpVpS7BZoCSEB5UAMyoGYTSN2DuJOI3ZaUuz2XLx4kUKFClldhlhMORBQDkQjdiIub+zYsVaXIE5AORBQDsR8auxEHOzxxx+3ugRxAsqBgHIg5lNjJ+JgW7ZssboEcQLKgYByIOZTYyfiYAEBAVaXIE5AORBQDsR8auxEHMzX19fqEsQJKAcCyoGYT42diIMdPHjQ6hLECSgHAsqBmE9LikmeFShQgDFjxmRuS9488MADVpcgTkA5EFAOxHwasZM88/HxoVevXvTq1UtLit2C6dOnW12COAHlQEA5EPNpgmIHcacJiuX2ZGRk4Omp/6fK75QDAeVANEGxywkLCyMoKIgmTZpYXYpp0tPTWb16NatXryY9Pd3qclzOkCFDrC5BnIByIKAciPk0Yucg7jRipyXFbo/+D11AORCDciAasRNxcUOHDrW6BHECyoGAciDmU2Mn4mCdO3e2ugRxAsqBgHIg5lNjJ+Jgq1evtroEcQLKgYByIOZTYyfiYHfffbfVJYgTUA4ElAMxnxo7EQdLTU21ugRxAsqBgHIg5lNjJ+Jg8fHxVpcgTkA5EFAOxHxaUkzyzMfHh+HDh2duS940bNjQ6hLECSgHAsqBmE8jdpJnBQoUoG/fvvTt21drxd6CxYsXW12COAHlQEA5EPNpgmIHcacJiuX2XLx4kUKFClldhlhMORBQDkQTFIsTSE9PJzo6mujoaC0pdgu+/PJLq0sQJ6AcCCgHYj6N2DmIO43YaUkxERERc2jEzsWEhYURFBREkyZNrC5FnMSgQYOsLkGcgHIgoByI+TRi5yAasZOrTp48SenSpa0uQyymHAgoB6IROxGX9+OPP1pdgjgB5UBAORDzqbETcbBGjRpZXYI4AeVAQDkQ86mxE3GwkydPWl2COAHlQEA5EPOpsRNxsLS0NKtLECegHAgoB2I+LSkmeebj48Nnn32WuS15c/fdd1tdgjgB5UBAORDzacRO8qxAgQIMGDCAAQMGaEmxW7BmzRqrSxAnoBwIKAdiPk134iDuNN2J3J5z585RvHhxq8sQiykHAsqBaLqTXBk7diyVK1emYMGCNGrUiMjIyBsev2DBAoKCgvD19SUoKIiFCxde99gePXrg4eFBaGio2WW7jIyMDHbv3s3u3bvJyMiwuhyXM2bMGKtLECegHAgoB2I+t2vs5s6dS+/evfnkk0/YunUr9913H+3btyc2Nvaax0dFRdGpUyc6d+7M9u3b6dy5Mx07dmTjxo3/OHbRokVs3LiRcuXK2fvHcGoXL16kdu3a1K5dm4sXL1pdjsvp37+/1SWIE1AOBJQDMZ/bNXYhISF069aN7t27U7NmTUJDQ6lQoQLjxo275vGhoaE8/PDD9OvXjxo1atCvXz9at279jxG5v/76izfffJOZM2fqgQG5LVpCSEA5EINyIGZzq8YuLS2NzZs306ZNmxyvt2nThvXr11/zPVFRUf84vm3btjmOz8jIoHPnzvTt25datWrlqpbU1FSSkpJyfIkAvPXWW1aXIE5AORBQDsR8XgMGDBhgdRFmOXnyJCNHjuTtt9+mQoUKma9v3bqV33777Zp/gT788ENeeukl6tSpk/nan3/+ycyZM/n0008BCA4OJjY2ljFjxmTeX9esWTOaNWt23VoGDRpEu3btCA4OJjg4mFGjRgHQqVMnpk2bRuPGjQkODqZVq1YMGjSIqlWr8v3335Oens6uXbvYsWMHNpuN2bNnU7NmTUaMGJF5bO3atZk0aRIFCxZk48aNHDhwgOTkZBYvXkylSpUYNWpU5rGNGzdm9OjRlCxZkl9//ZW4uDji4uL45ZdfKFmyJF9//XXmsf/6178YNmwYFStW5KeffuL8+fMcOnSIDRs24Ovry6RJk2jUqBHBwcGsXbsWgG7dujF//nwyMjLYuXMnO3fuJD09nblz51KjRg1GjhyZef66desyceJEChUqRFRUFIcOHSIpKYklS5Zw1113ERoamnlskyZNGD16NAEBAaxcuZJTp05x/PhxVq9eTYkSJRgzZgwtW7Zk8ODBtGzZkiFDhhAYGMiPP/5IcnIyBw4c4Pfff8fX15cpU6bQoEEDhg0blnn+6tWrM2vWLAC2bdvG7t27uXz5MvPmzaNatWo56q5Xrx7ffvstfn5+/Pbbbxw5coSzZ88SHh5O+fLlc9R97733MmrUKO68804iIiJISEjg2LFjrFmzhuLFizNu3Dh2797NypUrM99TuXJlFi5cyIULF9i3bx+bNm3C29ubqVOnUr9+/Rx116xZkxkzZuDp6cmWLVvYu3cvqampzJ8/n6pVq/Lll19mHlu/fn2++eYbihYtSmRkJLGxsZw9e5Zly5ZRpkwZvvrqq8xjmzZtSkhICGXKlGH58uWcPn2ao0ePEhkZib+/P+PHj6dZs2YMGTIk8z1VqlRhwYIFpKamsmfPHrZs2YKXlxfTp0+nTp06DB8+PPPYoKAgpk2bhre3N5s2beKPP/7g4sWLLFy4kCpVqhASEpJ5bIMGDRg3bhz+/v6sXbuWY8eOkZCQwIoVKyhdunSOups3b87IkSMpV64c4eHhnDt3jiNHjrBu3TqKFCnChAkTaNq0KUOHDs3xd23evHmkpaURExPDtm3bAJg5cya1atXK8XetVq1aTJ06FR8fH6Kjo9m/fz8XLlxg0aJFBAYG5qi7YcOGhIWFUaJECVatWsXx48eJj48nIiKCgIAARo8enXlsixYteO6552jevDlLliwhMTGRP//8k6ioKAoXLszEiRNd/nfE1fdUq1aN2bNn63fEdX5HTJs2jaJFi+b4HdGiRYscf9f0O8K9f0csX76c8PBw+vXrh6+v7y33QFe51VOxx48fp3z58qxfv57mzZtnvj548GCmT5/O3r17//GeAgUKMHXqVJ5//vnM12bOnEm3bt24dOkSmzdv5tFHH2XLli2Z99YFBgbSu3dvevfufd1aUlNTSU1NzdxPSkqiQoUKbvFUbEpKCkWKFAEgOTkZPz8/iytyLbt27aJ27dpWlyEWUw4ElAPRU7E3VKpUKby8vIiLi8vx+qlTpyhduvQ131OmTJkbHh8ZGcmpU6eoWLEi3t7eeHt7c+TIEd577z0CAwOvW4uvry/+/v45vkQA9u/fb3UJ4gSUAwHlQMznVo1dgQIFaNSoERERETlej4iIoEWLFtd8T/Pmzf9x/IoVKzKP79y5Mzt27GDbtm2ZX+XKlaNv374sX77cPj+IuLVChQpZXYI4AeVAQDkQ87ndkmJ9+vShc+fONG7cmObNmzNhwgRiY2Pp2bMnAF26dKF8+fIMHToUgHfeeYdWrVoxbNgwnnjiCX788UdWrlzJunXrAChZsiQlS5bMcQ0fHx/KlClD9erVHfvDOQkfHx/ef//9zG3Jm7/nSfIn5UBAORDzuV1j16lTJ06fPs0XX3zBiRMnqF27NuHh4VSqVAmA2NhYPD2zBipbtGjBnDlz6N+/P59++il33303c+fOpWnTplb9CE6vQIECjBgxwuoyXNa2bdto0qSJ1WWIxZQDAeVAzOdWD084My0pJlfFxcVRpkwZq8sQiykHAsqB6OEJcQIZGRkcPnyYw4cPa0mxWzBx4kSrSxAnoBwIKAdiPo3YOYg7jdhpuhMRERFzaMROxMVpCSEB5UAMyoGYTSN2DqIRO7nq0qVLFCxY0OoyxGLKgYByIBqxE3F5X3/9tdUliBNQDgSUAzGfGjsRB3vyySetLkGcgHIgoByI+dTY2VlYWBhBQUGap0gybd682eoSxAkoBwLKgZhPjZ2d9erVi5iYGKKjo60uRZyE5qwSUA7EoByI2dxu5QmxP29vb954443Mbckb/TsTUA7EoByI2TRiJ3nm6+tLWFgYYWFh+Pr6Wl2Oy/nzzz+tLkGcgHIgoByI+dTYiThYq1atrC5BnIByIKAciPnU2Eme2Ww24uPjiY+Pxx7TIF6+DPPmwSOPgL8/VK4MDz0EY8dCaqrpl3O4WbNmWV2COAHlQEA5EPNpgmIH0QTFuXPuHDz2GPz227W/X7EiTJoErVubdkmHS09Px8vLy+oyxGLKgYByIJqgWNxYfDw8+GDOpq50aShRIms/NhbatzdG9FzV0KFDrS5BnIByIKAciPk0YucgGrG7sfR04+PW1auN/YAAmD4dHn4YPD1h82bo2xdWrTK+7+FhfDTbs+dtX1pERMQyGrETtzR5clZTV7YsREZC27ZGUwfQqBGsWAGvvmrs22zwn/9AcLAl5d4WLfotoByIQTkQs6mxE8ulp+ds0GbOhOrV/3mctzdMnAgffJD1Wr9+8P339q/RTC+//LLVJYgTUA4ElAMxnxo7sdz8+XDwoLHdurVxn931eHjAsGEweHDWa6+8Ajt22LdGM/3yyy9WlyBOQDkQUA7EfGrsxFI2G2S/d/ijj3L3vn79oEsXY/vCBXjqKUhKMr8+e6hWrZrVJYgTUA4ElAMxnxo7OwsLCyMoKIgmTZpYXYppvL29efnll3n55Zdvezmc5cth+3Zju3Hj3E9j4uEB48cb7wE4dAh6976tUhzm4sWLVpcgTkA5EFAOxHxq7OysV69exMTEEB0dbXUppvH19WXKlClMmTLltpcUyz5a16+f0bDlVqFCxv11RYsa+5Mnww8/3FY5DnH69GmrSxAnoBwIKAdiPjV2Ypk9e2DtWmO7enV48sm8n6NyZRg9Omu/Rw9w9t+T9erVs7oEcQLKgYByIOZTYyd5ZrPZSElJISUl5baWFFu0KGu7R4+sqU3y6uWXjXvsABIS4L33brkkh1iyZInVJYgTUA4ElAMxnyYodhBNUPxPzZrBxo3G9sGDUKXKrdd04gTUrAmJicZ+RIQx4bEzunDhAoULF7a6DLGYciCgHIgmKBY3ceJEVlNXu/btNXVgTGo8fHjWfs+ecOnS7Z3TXkJCQqwuQZyAciCgHIj51NiJJX76KWv7iSfMOWf37nDffcb2wYMQGmrOec3Wv39/q0sQJ6AcCCgHYj41dmKJH3/M2jarsfP0hDFjsu7VGzTIGBl0NlpCSEA5EINyIGZTYycOl5oKv/5qbJcrZ6wDa5a6deH1143tlBT4+GPzzm2WHj16WF2COAHlQEA5EPOpsROH+/33rPvfHnro1p+GvZ4vvoDixY3tKVPA2aYQ/MEVJtsTu1MOBJQDMZ8aO3G41auzth94wPzzBwTAgAFZ+717G0uXOQt3WoVEbp1yIKAciPnU2EmeeXl58eyzz/Lss8/i5eWV5/evW5e1ff/9JhaWzRtvQI0axvb69TBnjn2ucytOOOONf+JwyoGAciDmU2MneVawYEHmzZvHvHnzKFiwYJ7ea7NlfTR6553GyhH24OMDo0Zl7X/wATjLkozp6elWlyBOQDkQUA7EfGrs7CwsLIygoCANt//PwYNw9qyx3aRJ3taGzat27eCRR4ztY8dg7Fj7XSsvAgMDrS5BnIByIKAciPnU2NlZr169iImJIdrZ7uC3SPZ/DY7odYODs5rHoUOzVqawUmRkpNUliBNQDgSUAzGfGjvJs5SUFDw8PPDw8CAlJSVP73V0Y1enDrz4orF9+jR8+aX9r3kzzz//vNUliBNQDgSUAzGfGjtxKEc3dgCff27ccwcQEgInTzrmutcz1lk+ExZLKQcCyoGYz8Nmc6aJINyX2Yv8WiklJYUiRYoAkJycjJ+fX67ed+UKFCsGFy5AYCD8+acdi/ybt94yVqUAePtt+Oorx11bRETkeszuDzRiJw4TE2M0deC40bqr+veHwoWN7XHj4PBhx14/Oy0hJKAciEE5ELOpsROH2bUra7thQ8deu3RpePddY/vy5ZwTGDva22+/bd3FxWkoBwLKgZhPjZ04zB9/ZG3XrOn467//PpQoYWxPmwa7dzu+BoApU6ZYc2FxKsqBgHIg5lNjJw6zb1/WdrVqjr9+8eLQr5+xbbMZH89a4aGHHrLmwuJUlAMB5UDMp8ZO8szLy4tHHnmERx55JE9Lil0dsfP0hLvvtlNxN/Hmm1CunLG9aBFs2OD4Gvbu3ev4i4rTUQ4ElAMxnxo7ybOCBQuydOlSli5dmuslxWy2rBG7ypWhQAE7FngDhQrBZ59l7ffrZ9TmSLl9iljcm3IgoByI+dTYiUMcPw5X5zKuXt3aWl55Be65x9hevRoiIhx7/TvuuMOxFxSnpBwIKAdiPjV24hBW31+XnY8PDByYtf/JJ44dtduxY4fjLiZOSzkQUA7EfGrs7CwsLIygoCCaOHriNjtKSUnBz88PPz+/XC8plv2JWKtH7AA6dID69Y3tTZuM++0cpX379o67mDgt5UBAORDzqbGzs169ehETE0N09rW03MCFCxe4cHW24VxwphE7MB7gyD4v6KefQnq6Y649adIkx1xInJpyIKAciPm0pJiD5PclxR57DJYuNbaPHYPy5e1ZYe7YbPCvf0FUlLE/Ywa8+KK1NYmISP6iJcXEJV39KNbPL2u6Eat5eMDgwVn7n31mrEphb1pCSEA5EINyIGbTiJ2D5OcRu7Q0Y53W9HRo0AC2bHFElbn30EPwyy/G9oQJ8Npr9r1eamoqvr6+9r2IOD3lQEA5EI3YiQs6dCjr/jVnuL/u77KP2n3xBVy6ZN/rffXVV/a9gLgE5UBAORDzqbETu8v+4IQzPBH7d02bwv/9n7F97BiMH2/f6z3zzDP2vYC4BOVAQDkQ86mxkzzz9PTk/vvv5/7778fT8+YRyj7ViTOO2EHOee2GDIHkZPtd6/fff7ffycVlKAcCyoGYzy0bu7Fjx1K5cmUKFixIo0aNiIyMvOHxCxYsICgoCF9fX4KCgli4cGGO7w8YMIAaNWrg5+dHiRIleOihh9i4caM9fwSnVqhQIVavXs3q1aspVKjQTY93tqlOrqVePejUydiOj4fRo+13rbJly9rv5OIylAMB5UDM53aN3dy5c+nduzeffPIJW7du5b777qN9+/bExsZe8/ioqCg6depE586d2b59O507d6Zjx445Grdq1aoxZswYdu7cybp16wgMDKRNmzbEx8c76sdyaa4wYgfw+efG/HYAI0bAuXP2uY6Xl5d9TiwuRTkQUA7EfG7X2IWEhNCtWze6d+9OzZo1CQ0NpUKFCowbN+6ax4eGhvLwww/Tr18/atSoQb9+/WjdujWhoaGZx7zwwgs89NBDVKlShVq1ahESEkJSUpKWgsmlqyN2pUtDsWLW1nIj1avDyy8b2+fOwciR9rnOkSNH7HNicSnKgYByIOZzq8YuLS2NzZs306ZNmxyvt2nThvXr11/zPVFRUf84vm3bttc9Pi0tjQkTJlCsWDHq1atnTuEuJiUlhYCAAAICAm66pFhiIpw8aWw744MTf/ff/xpryQKEhsKpU+Zfo2XLluafVFyOciCgHIj53KqxS0hIID09ndKlS+d4vXTp0sTFxV3zPXFxcbk6fsmSJRQpUoSCBQsyatQoIiIiKFWq1HVrSU1NJSkpKceXO0lISCAhIeGmxx08mLV9zz12LMgkgYHw+uvGdkoKBAebf405c+aYf1JxOcqBgHIg5vMaMGDAAKuLMMv58+cJCQmhe/fuVKhQIfP1tWvXsmPHDt58881/vGfgwIE89dRT1KlTJ/O1HTt2sHjxYvr375/5Wvny5XnppZd4/vnnSUpKYvjw4bzwwgvXnZx30KBBtGvXjuDgYIKDgxk1ahQAnTp1Ytq0aTRu3Jjg4GBatWrFoEGDqFq1Kt9//z3p6ens2rWLHTt2YLPZmD17NjVr1mTEiBGZx9auXZtJkyZRsGBBNm7cyIEDB0hOTmbx4sVUqlSJUaNGZR7buHFjRo8eTcmSJfn111+Ji4sjLi6OX375hZIlS/L1119nHvuvf/2LYcOGUbFiRX766SfOnz/PoUOH2LBhA76+vkyaNIlGjRoRHBzM2rVrAejWrRvz588nIyODnTt3snPnTtLT05k7dy41atTgk08Ws2tXLQDat7/Inj3fUKhQIaKiojh06BBJSUksWbKEu+66i9DQ0MxamjRpwujRowkICGDlypWcOnWK48ePs3r1akqUKMGYMWNo2bIlgwcPpmXLlgwZMoTAwEB+/PFHkpOTOXDgAL///ju+vr5MmTKFBg0aMGzYsMzzV69enVmzZgGwbds2du/ezeXLl5k3bx7PPVedCRO8yMjwYtOmKzz7bDJz5kzAz8+P3377jSNHjnD27FnCw8MpX758jrrvvfdeRo0axZ133klERAQJCQkcO3aMNWvWULx4ccaNG0efPn1y/NlXrlyZhQsXcuHCBfbt28emTZvw9vZm6tSp1K9fP0fdNWvWZMaMGXh6erJlyxb27t1Lamoq8+fPp2rVqnz55ZeZx9avX59vvvmGokWLEhkZSWxsLGfPnmXZsmWUKVOGr776KvPYpk2bEhISQpkyZVi+fDmnT5/m6NGjREZG4u/vz/jx42nWrBlDhgzJfE+VKlVYsGABqamp7Nmzhy1btuDl5cX06dOpU6cOw4cPzzw2KCiIadOm4e3tzaZNm/jjjz+4ePEiCxcupEqVKoSEhGQe26BBA8aNG4e/vz9r167l2LFjJCQksGLFCkqXLp2j7ubNmzNy5EjKlStHeHg4586d48jD0mfwAAAgAElEQVSRI6xbt44iRYowYcIEmjZtytChQ3P8XZs3bx5paWnExMSwbds2AGbOnEmtWrVy/F2rVasWU6dOxcfHh+joaPbv38+FCxdYtGgRgYGBOepu2LAhYWFhlChRglWrVnH8+HHi4+OJiIggICCA0aNHZx7bokUL1q1bR4UKFViyZAmJiYn8+eefREVFUbhwYSZOnOjyvyOuvqdatWrMnj37ur8jRo4cmXls3bp1mThxotP/jqhWrVqOuuvVq8e33357S78j7rnnHo4fP57jd0SLFi1y/F3T7wj3/h2xfPlywsPD6devnzmTVdvcSGpqqs3Ly8v2ww8/5Hj97bfftrVq1eqa76lQoYItJCQkx2shISG2ihUr3vBaVatWtQ0ZMuS637906ZItMTEx8+vo0aM2wJaYmJjLn8Z5JScn2wAbYEtOTr7hsaNH22zGqqw22+TJjqnPDH37ZtXdo4e55x44cKC5JxSXpByIzaYciM2WmJhoan/gVh/FFihQgEaNGhEREZHj9YiICFq0aHHN9zRv3vwfx69YseK6x19ls9lITU297vd9fX3x9/fP8ZUfHT2atV2xonV15NWHH0LRosb2d98Zq2eYJftIsORfyoGAciDmc6vGDqBPnz5MnDiRSZMmsWfPHt59911iY2Pp2bMnAF26dKFfv36Zx7/zzjusWLGCYcOGsXfvXoYNG8bKlSvp3bs3YDwo8PHHH7NhwwaOHDnCli1b6N69O8eOHaNDhw6W/IyuJPssM9k+HXd6JUvCe+8Z21eugJk3LGjRbwHlQAzKgZjN7Rq7Tp06ERoayhdffEH9+vVZu3Yt4eHhVKpUCYDY2FhOnDiReXyLFi2YM2cOkydPpm7dukyZMoW5c+fStGlTwJhjaO/evTzzzDNUq1aNxx57jPj4eCIjI6lVq5YlP6MryT5id9dd1tVxK95912jwAGbMgN27zTlv165dzTmRuDTlQEA5EPN52Gw2m9VF5AdJSUkUK1aMxMREl/9Y9uLFi7Rq1QowHky50eoTFSsazV1AgH2mDrG3ESPggw+M7aefhgULbv+ckydP5pVXXrn9E4lLUw4ElAMxvz9wuxE7sb9ChQoRHR1NdHT0DZu6K1fg+HFj25U+hs2uVy8oU8bY/uEH2Lz59s9Zs2bN2z+JuDzlQEA5EPOpsRO7OXEC0tONbVd6cCK7woUh+73NZtznfP78+ds/ibg85UBAORDzqbETu8l+f52rjtgBvPYa/O8WTZYtg8jI2zvfOXstQisuRTkQUA7EfGrsJM8uXLhAYGAggYGBXLhw4brHuUtjV6AAfPZZ1v6HHxoz3N2q7JNhS/6lHAgoB2I+NXaSZzabjSNHjnDkyBFu9OyNuzR2AJ07w9VbYaKiYOHCWz/Xzz//bE5R4tKUAwHlQMynxk7sJvtyu+XKWVeHGby9YdiwrP2PPoLLl2/tXK+99po5RYlLUw4ElAMxnxo7sZts0wVmPlnqyh57DO6/39jevx8mTLi184SGhppXlLgs5UBAORDzaR47OwsLCyMsLIz09HT27dvnFvPYpaSkUKRIEQCSk5Px8/O75nGtW8OvvxrbSUlZS3S5st9/h//NXU1AABw4AC7+xykiIhbSPHYuplevXsTExBAdHW11KQ53dcSucGH4Xx/o8u69Fzp1Mrbj440JjPNKSwgJKAdiUA7EbBqxcxB3WnkityN2d9wBZ89ClSpw8KAjK7SvQ4egRg3jHrtChYyPZcuXz/37ExISKFWqlP0KFJegHAgoB6IRO3ECHh4eBAUFERQUhIeHxzWPuXTJaOoAypZ1YHEOUKWKsSIFwMWLOadCyY358+ebX5S4HOVAQDkQ86mxkzwrXLgwu3fvZvfu3RQuXPiax5w8mbXtDg9O/F3//lCsmLE9eTLs2pX79zZr1sw+RYlLUQ4ElAMxnxo7sYvsT8S624gdQMmS8PHHxnZGhjFpcW4dO3bMPkWJS1EOBJQDMZ8aO7GL7HPYueOIHcBbb2VNvBwenvUE8M3otlYB5UAMyoGYTY2d5NmFCxeoVasWtWrVuu6SYu4+YgfGgxPZH2j74ANj9O5mKrj6MhxiCuVAQDkQ86mxkzyz2WzExMQQExNz3f/bdLfJia/npZegXj1je/NmmDnz5u+Jioqyb1HiEpQDAeVAzKfGTuwi+0ex7jpiB+DpCSNHZu1/+CGcP3/j93Ts2NG+RYlLUA4ElAMxnxo7sYv8MmIH8NBD8PjjxvaJEzBkyI2PHzdunP2LEqenHAgoB2I+TVDsIPltguImTWDTJmNEKy0NvLwcXaVjHTwIQUHGz1qgAOzeDVWrWl2ViIg4O01QLC7h6ohdQID7N3UAd98NffoY22lp8N571z9WSwgJKAdiUA7EbBqxs7OwsDDCwsJIT09n3759+WLEzmYzRq2uXDEeLNi2zYoqHe/8eahePaupXb4c2rS51nHnKVq0qGOLE6ejHAgoB6IRO5fTq1cvYmJiiI6OtroU03h4eFCpUiUqVap0zSXFEhONpg7gzjsdXJyFihaFYcOy9nv3NtaT/bvvvvvOcUWJ01IOBJQDMZ+3mSeLj49n+/btHD58mMTERFJTUylSpAjly5enevXq1K5d28zLiUUKFy7M4cOHr/v9hISs7fy2tvWLL8LYsbBhA+zZA2FhRoOXXZtrDeNJvqMcCCgHYr7bbuwOHTrEd999x8KFC/njjz+w2Wx4eXlRrFgxfH19OXPmDKmpqQAUK1aMRx55hC5dutC2bdvbLl6cU35u7Dw9YfRouPdeY3/AAKPZCwjIOmbPnj0EBQVZUp84D+VAQDkQ893yR7Hx8fG8/PLLNG7cmD/++IOePXuyefNmzpw5w+XLl0lISOCvv/7i4sWLJCcnc/DgQcLCwihatChvvPEG9erVY82aNWb+LOIk8nNjB8YTwa+8YmwnJsInn+T8vu6nEVAOxKAciNluacRu+fLlvPfeezz//PMcOnSI4sWL3/D4woULU7lyZSpXrswLL7yAzWZj3rx59OnTh9atWzNs2LBr3qslzunixYu0atUKgLVr11KoUKEc3z99Oms7PzZ2YMxlN3++8UDFxInQowc0amR8r1ixYtYWJ05BORBQDsR8eR6xW7p0KV9//TWRkZF88sknN23qrsXDw4OOHTuyefNm7rjjDrp165bnc4h1MjIy2LRpE5s2bSLjGouj5vcROzAmZf7vf41tmw169oT0dGN/9+7d1hUmTkM5EFAOxHx5buz27NnD4sWLKVGihCkFfPTRRzz33HPs2rXLlPOJ9dTYGd5+25i0GIzJmseONbbbtWtnXVHiNJQDAeVAzJfnxu7999/H09PcWVLatGmjJ2bdiBo7Q4EC8M03WfuffAJ//QWTJk2yrihxGsqBgHIg5jO1Q9uxYwfvvPMOv/zyi5mnFRejxi5Ly5bQvbuxff68MfVJ//79rS1KnIJyIKAciPlMbewGDBjA+PHj6dSpk5mnFReTvbErWdK6OpzFsGFZ053Mnw9dusy2tiBxClpKSkA5EPOZOkFxmTJlWL58OaezPxYp+c7Vxq5oUfD1tbYWZ3DHHfDll9Cli7EfGfkcFy5A4cLW1iXW+vDDD60uQZyAciBmM3XEztfXl+rVq/PMM8+YeVpxQqVKlaLUdT5nvdrY5fePYbN76SX497+N7cOHPfjiC2vrEeuFhIRYXYI4AeVAzGZqY/fZZ5/xzjvvsH37djNP69LCwsIICgqiSZMmVpdiGj8/P+Lj44mPj8fPzy/H99LT4cwZY1uNXRYPDxg3znigAowRvJ07ra1JrNWhQwerSxAnoByI2Uxt7H777TfWrFlDkyZN+Pe//80XX3zB2rVrSUtLM/MyLqVXr17ExMQQHR1tdSkOce4cXJ3aTvfX5VStGnz8sbF95YoxafE1pgGUfCIqKsrqEsQJKAdiNlMbu+DgYN5++226d+/O6dOn+fzzz3nwwQcpXrw4//73vxk4cCB79+4185LiZPRE7I199BFUqHABgKgomDDB4oLEMhUqVLC6BHECyoGYzdTGrnLlyvTr14+xY8eyfft2Tp06xZw5c+jatStxcXF89tln/PvqjUbisi5evMgDDzzAAw88wMWLF3N8T43djfn6Qp8++zP3+/aFI0csLEhERNyKqU/Fvv/++7z22muUL1+eZ555hnr16tGhQ4fMewhOnjzJyZMnzbykWCAjI4M1a9Zkbmenxu7mAgJ20b17PSZOhORkeO01WL7cuA9P8o+jR49aXYI4AeVAzGbqiF3dunX57rvv6Nq1KxcuXPjH90uXLk3dunXNvKQ4mewz3aixu7bmzZszciTcdZexHxEB331nbU3ieM2bN7e6BHECyoGYzdy1wf6nSpUqCms+pRG7m5s3bx7FisG332a91qcPxMZaV5M43rx586wuQZyAciBmy3Njd/LkSVJSUkwv5I8//jD9nOJ4auxurk+fPgC0awevvGK8dv48vP462GwWFiYOdTUHkr8pB2K2PDd2np6evPrqq5w4ccK0IqZPn87o0aNNO59YR43dzQ0bNixzOyQEypUztpcvh8mTLSpKHC57DiT/Ug7EbB42W97HCA4fPkznzp155ZVX6Nq1K56et/aJ7pEjR/jiiy/IyMhg4sSJeHl53dJ5XEFSUhLFihUjMTERf39/q8u5LSkpKRQpUgSA5OTkHJMUP/44/PSTsR0XB6VLW1Gha1m6FB57zNj294fdu7PuvxMREfdmdn9wSx1ZYGAgP//8Mxs2bKBGjRoMHTqU7du3k5seMSkpifDwcDp37kyTJk1o2rQpkydPduumzh0VLlyYwtdY7DT7iN0ddziwIBfy90W/H30UXn7Z2E5Kgq5dNXFxfqDF3wWUAzHfLY3YZbd161ZGjBjBokWL8Pb2pkmTJpQvX57ixYtTrFgx0tLSOHPmDGfOnOHPP/9k586dBAQE0K1bN955553rrjfqbtxpxO5GataEvXuNkafERKurcU7Hjx+n3NXPX//n7FmoWxeOHTP2Q0Lg3XctKE4c5lo5kPxHORCnGLHLrkGDBsyaNYuTJ08yZcoUGjZsyJkzZ4iMjGTOnDn8+OOP7Ny5Ex8fH5599llWrVrFX3/9xcCBA/NNU5efnDtn/LN4cWvrcGbLli37x2slSsCUKVn7/fppLVl3d60cSP6jHIjZbnmC4tatW/PLL79k7hctWpSnn36ap59+2pTC3EVYWBhhYWGkp6dbXYpDqLG7uVq1al3z9datjWlPQkIgNRVefBGio43VKsT9XC8Hkr8oB2K2Wx6x27JlS75pVm5Hr169iImJITo62upSTHPp0iUeffRRHn30US5dupTtdeML1NjdSOINPqMePBjq1DG2d+6E/v0dVJQ43I1yIPmHciBmu+XGLjExkf/85z9aDiUfSk9PJzw8nPDw8BzNffbfT2rsru/8+fPX/V7BgjBjBhQoYOx/+SX8+quDChOHulEOJP9QDsRst3WPXUZGBr1796Z58+a88cYbLFiwgLNnz5pVm7iYqx/Dghq7G6lZs+YNv1+3LgwZYmzbbMZHsqdOOaAwcaib5UDyB+VAzHbLjV3FihWZOHEiCxYsYP369XTr1o2DBw/ywgsv0KpVKz766CNWrlzJlStXzKxXnJgau9xZsWLFTY959114+GFjOy4OunTRFCjuJjc5EPenHIjZbrmxGz9+fOa2h4cHjRo14oMPPuDnn39m5cqVtGvXjtWrV9O8eXNGjBiRqznuxLWpscudbt263fQYT0+YPj1rgufly2HECDsXJg6VmxyI+1MOxGy33Ni1a9fuut8rUKAADzzwAIMGDWLDhg1cvnz5Hzfai/tRY5c7X331Va6OK10aZs4EDw9j/5NPYP16OxYmDpXbHIh7Uw7EbLc9j93NeHl58fHHH/Paa6/xytUVz8UtqbHLnf55eNS1deusJ2PT0+G55+DMGTsVJg6VlxyI+1IOxGy33Nht3bo1T8c/9dRTBAYGsnjx4lu9ZK6NHTuWypUrU7BgQRo1akRkZOQNj1+wYAFBQUH4+voSFBTEwoULM793+fJlPvzwQ+rUqYOfnx/lypWjS5cuHD9+3N4/hsvJ/tyMGrvry+sSQv/9L7RqZWwfPQqvvGI8VCGuTUtJCSgHYr5bbuxeffXVPL+nT58+fPPNN7d6yVyZO3cuvXv35pNPPmHr1q3cd999tG/fntjY2GseHxUVRadOnejcuTPbt2+nc+fOdOzYkY0bNwJw4cIFtmzZwqeffsqWLVv44Ycf2LdvH48//rhdfw5n5ufnh81mw2az4efnl/m6Ruxy5z//+U+ejvf2hlmzoGRJY3/xYmMSY3Ftec2BuCflQExnu0UVKlS4pfc1bdr0Vi+ZK/fee6+tZ8+eOV6rUaOG7aOPPrrm8R07drS1a9cux2tt27a1Pffcc9e9xu+//24DbEeOHMl1XYmJiTbAlpiYmOv3uJoePWw2YyzJZtuyxepqnNfYsWNv6X1Ll2b9+/XystlWrTK3LnGsW82BuBflQMzuD255xO7YsWM8/PDDhIaGsmvXrly/z56rVaSlpbF582batGmT4/U2bdqw/jp3nUdFRf3j+LZt2173eDAmZ/bw8KC4hqVy0Ihd7jRv3vyW3vfII8YDFGDcb9epExw7ZmJh4lC3mgNxL8qBmO22Hp745ZdfeO+996hXrx5ly5blxRdfZNKkSdf92HPXrl342nHhy4SEBNLT0yl9dY6I/yldujRxcXHXfE9cXFyejr906RIfffQRL7zwAv7+/tetJTU1laSkpBxf7uLSpUt06NCBDh065HjSWY1d7tzOai2ffw5X/z/k1Cno0MFYV1Zcj1btEVAOxHxeAwYMGHArbxw/fjzh4eFUrlyZjIwMDhw4wLZt2/jpp5/46quvmDFjBjExMZw5c4b4+HiWL19Ojx496N+/P7Vr1zb5xzCcP3+ekJAQunfvToUKFTJfX7t2LTt27ODNN9/8x3sGDhzIU089RZ2rC3QCO3bsYPHixf94Wuny5ct06tSJc+fOMW/evBs2qYMGDaJdu3YEBwcTHBzMqFGjAOjUqRPTpk2jcePGBAcH06pVKwYNGkTVqlX5/vvvSU9PZ9euXezYsQObzcbs2bOpWbMmI0aMyDy2du3aTJo0iYIFC7Jx40YOHDhAcnIyixcvplKlSowaNSrz2MaNGzN69GhKlizJr7/+SlxcHHFxcfzyyy+ULFmSr7/+OvPYf/3rXwwbNoyKFSvy008/cf78eQ4dOsSGDRvw9fVl0qRJNGrUiCFDhjB+/HhiYmJ49dVXmT9/PhkZGUye7Mfp0wUBKFt2DDVrVmfkyJGZ569bty4TJ06kUKFCREVFcejQIZKSkliyZAl33XUXoaGhmcc2adKE0aNHExAQwMqVKzl16hTHjx9n9erVlChRgjFjxtCyZUsGDx5My5YtGTJkCIGBgfz4448kJydz4MABfv/9d3x9fZkyZQoNGjRg2LBhmeevXr06s2bNAmDbtm3s3r2by5cvM2/ePKpVq5aj7nr16vHtt9/i5+fHb7/9xpEjRzh79izh4eGUL18+R9333nsvo0aN4s477yQiIoKEhASOHTvGmjVrKF68OOPGjaNUqVLMnz8/8z2VK1dm4cKFXLhwgX379rFp0ya8vb2ZOnUq9evXz1F3rVo1SU6ex86dNUlO9ubYMdi3L4FTpyZTtWpVvvzyy8xj69evzzfffEPRokWJjIwkNjaWs2fPsmzZMsqUKcNXX32VeWzTpk0JCQmhTJkyLF++nNOnT3P06FEiIyPx9/dn/PjxNGvWjCFDhmS+p0qVKixYsIDU1FT27NnDli1b8PLyYvr06dSpU4fhw4dnHhsUFMS0adPw9vZm06ZN/PHHH1y8eJGFCxdSpUoVQkJCMo9t0KAB48aNw9/fn7Vr13Ls2DESEhJYsWIFpUuXzlF38+bNGTlyJOXKlSM8PJxz585x5MgR1q1bR5EiRZgwYQJNmzZl6NChOf6uzZs3j7S0NGJiYti2bRsAM2fOpFatWjn+rtWqVYupU6fi4+NDdHQ0+/fv58KFCyxatIjAwMAcdTds2JCwsDBKlCjBqlWrOH78OPHx8URERBAQEMDo0aMzj23RogXBwcE0aNCAJUuWkJiYyJ9//klUVBSFCxdm4sSJLv07Invd1apVY/bs2WRkZLBz50527txJeno6c+fOpUaNGvn+d8S+ffvw8PDI8TuiRYsWOf6u5eV3RM2aNZkxYwaenp5s2bKFvXv3kpqayvz58/U7wkl/Ryxfvpzw8HD69etnzuDXrX6G27Bhwxz7ly5dsq1cudL28ccf2+69916bt7e3zcPDw+bp6Zn51bdv39v85PjGUlNTbV5eXrYffvghx+tvv/22rVWrVtd8T4UKFWwhISE5XgsJCbFVrFgxx2tpaWm2J5980la3bl1bQkLCTWu5dOmSLTExMfPr6NGjbnOPXXJysg2wAbbk5OTM16tXN+7/8ve3sDgXsHXr1ts+x6ZNNpuvb9Y9dxMmmFCYOJQZORDXpxyI09xj98MPP+TY9/X1pXXr1gwePJiNGzdy5swZlixZwsiRIwkNDWXnzp0MHz78tprQmylQoACNGjUiIiIix+sRERG0aNHimu9p3rz5P45fsWJFjuMvX75Mx44d2b9/PytXrqTk1ccTb8DX1xd/f/8cX+7u6kex+hj2xjZs2HDb52jUCMaOzdp/4w1Yvfq2TysOZEYOxPUpB2I6U9pDJzJnzhybj4+P7bvvvrPFxMTYevfubfPz87MdPnzYZrPZbJ07d87xhOxvv/1m8/LysgUHB9v27NljCw4Otnl7e9s2bNhgs9lstsuXL9sef/xx21133WXbtm2b7cSJE5lfqampua7LnZ6Kvd6I3dURpHr1LCzOBcTHx5t2rnfeyRq1u+MOm+3AAdNOLXZmZg7EdSkH4jQjds6qU6dOhIaG8sUXX1C/fn3Wrl1LeHg4lSpVAiA2NpYTJ05kHt+iRQvmzJnD5MmTqVu3LlOmTGHu3Lk0bdoUMJ7+Xbx4MceOHaN+/fqULVs28+tGT87mN5cuZd3ErxG7G8u+zvLtGjkS2rY1ts+cgf/7P0hMNO30Ykdm5kBcl3IgZvOw2TSHvSMkJSVRrFgxEhMTXf5j2ZSUFIoUKQJAcnIyfn5+xMVB2bLG9594AhYtsrDAfCYxEZo1g717jf327eGnn8DLy9q6RETk5szuD9xuxE6soalOcs/sJYSKFTMauTvuMPZ//hk++MDUS4gdaCkpAeVAzKfGTvKscOHCJCcnk5ycTOHChQE1dnnRu3dv089ZtSrMn28sPwbGkmMTJ5p+GTGRPXIgrkc5ELOpsZM88/DwwM/PDz8/Pzw8PAA1dnnx7bff2uW8Dz4IY8Zk7ffsaYzeiXOyVw7EtSgHYjY1dmIKNXa51759e7udu0cPeOcdYzs93ViZYtMmu11OboM9cyCuQzkQs6mxkzxLTU2la9eudO3aldT/PQqrxi73du7cadfzf/klPPussZ2SAo8+CgcP2vWScgvsnQNxDcqBmE2NneTZlStXmDp1KlOnTuXKlSuAGru8KG7nf0FeXjB9Otx3n7F/6hS0awfx8Xa9rOSRvXMgrkE5ELPlurEbPnw4lSpVws/Pj/bt2/PHH39kfm/ixIn06tWLrl272qNGcQFq7HKvaNGidr9GwYLw448QFGTsHzgAjz1mjOCJc3BEDsT5KQditlw1dl9++SXDhw/ngQceoGfPnly8eJGGDRuydOlSALp3785jjz3G9OnT7VqsKwoLCyMoKIgmTZpYXYpdZW/sihWzrg5XsGfPHodcp0QJWLYMypc39n//3bjnLi3NIZeXm3BUDsS5KQdiNu/cHPTzzz+zc+dOyl6dgRaIjo7mjTfeoGDBgrRu3RpPT32qey29evWiV69emRMQuqukpKxtF59/2e4efvhhh12rQgXjydiWLY0/o59/hhdfhNmzs6ZGEWs4MgfivJQDMVuuurHGjRvnaOoAmjRpwpo1a5g2bRq//vqrXYoT16HGLvemTJni0OvVqWNMYFyokLE/fz507w4ZGQ4tQ/7G0TkQ56QciNly1dh5eXlx6dIl9u7dy+LFizNfL1y4MFOnTuW3335j2bJlditSnN/581nbauxurH///g6/ZqtWsHAh+PgY+1OnwltvgRYUtI4VORDnoxyI2XLV2L311lu8+eabdOjQgddff/0f3//000+pX78+Plf/qyH5ztUROx8f8PW1thZnZ9USQm3bwty5WWvIjh0LH32k5s4qWkpKQDkQ83nYbLn/tX7gwAH8/Pz+8bHsVfv37+eee+4xrTh3YvYiv1ay2WwkJCQAUKpUKTw8PKha1ZgrrWRJ+N+35DquXLmCt4U3uM2YAV26ZDV0X3wBn35qWTn5ltU5EOegHIjZ/UGennioWrXqdZs6QE1dPuHh4UFAQAABAQGZS4pdHbHTk/s3N3LkSEuv/9JLMG5c1v5//wuDB1tXT35ldQ7EOSgHYjaHPMo6adIkBg4cyNmzZx1xObHA1XvsXHww0iGee+45q0ugRw/I/t+T/v1h4EDr6smPnCEHYj3lQMxml8YuPj6ejGyP3L366qs89dRTvP3226xatcoelxQHSk1NzZzGJTU1lbQ0uHTJ+J4au5tbt26d1SUA8N57MGJE1v5//wuff25dPfmNs+RArKUciNlMbexWrVpFqVKlKFOmDP7+/rRv355JkyaRmJhI7dq1mTZtGuPHjzfzkmKBK1euMHbsWMaOHcuVK1dyPBGrj2JvrlKlSlaXkOn99421Za8aMMBo8PRAhf05Uw7EOsqBmM3UOzb/+9//0qtXL0qXLk1MTAwRERF0796dN998kyeeeIJWrVpl3nQv7kNTneRNenq61SXk0KcPeHrCu+8a+wMHGn+mISHwv1soxQ6cLQdiDeVAzGZqY1enTh0+/3hKy4cAACAASURBVNtnOTt37mTWrFksXLiQiIgIvv76azMvKU4g++TEGrG7uRMnTlhdwj/07m1Mg/L228Z+aKjx5zphQtb0KGIuZ8yBOJ5yIGaz+8MTderUYejQoezdu5eEhASef/55e1/SqeSHtWI1Ypc39957r9UlXNNbb8GkScboHRjbzz2ntWXtxVlzII6lHIjZTG3s7rvvPhYtWmTmKV1er169iImJITo62upS7EYjdnmzYMECq0u4rldeMSYxvjrX+Pz58MQTcOGCtXW5I2fOgTiOciBmM7Wxe/7559m1axcTJ04087Ti5DRilzfvvPOO1SXc0LPP5lxbdtkyaN1aE0+bzdlzII6hHIjZTG3sfv/9dyZPnkyPHj2oWrUq3bt3Z8aMGRw7dszMy4iTyT5ip8bu5kZkn2PESbVtCytWZP15btgALVoYq4uIOVwhB2J/yoGYLU9Lit3MfffdR61atfD392fHjh2sX7+e5ORkPDw8qFKlCg8++CAvvfQSrVq1MuuSLsOdlhTLyMggNjYWgIoVKxIa6sl77xnfmzMHOnWysDgx1fbt8MgjcPy4sR8QAEuXghvfMioi4lCWLil2M7Vr12b8+PEMHz6cZcuWcfbsWdavX8+gQYOoUqUKs2fPppP+q+/yPD09CQwMJDAwEE9PT30Um0eutOh3vXoQFQVBQcZ+fDw88AAsWWJpWW7BlXIg9qMciNlMbewKFCiQY9/Ly4tmzZrRr18/li9fztmzZ9m5c6eZlxQnoIcn8ubVV1+1uoQ8qVgR1q2D++839i9cMB6oyL7erOSdq+VA7EM5ELOZ2thVr16dNWvWXPf73t7elCpVysxLigXS0tLo27cvffv2JS0tTSN2efTzzz9bXUKelSgBy5cb058AZGTAG28Y895duWJtba7KFXMg5lMOxGymNnb/+c9/+OGHH4iIiDDztOJkLl++zMiRIxk5ciSXL1/WwxN5VLduXatLuCW+vjBzJvTtm/Xa119D+/Zw9qx1dbkqV82BmEs5ELOZ2tgtW7aMOXPm0K5dOx555BFGjRrF1q1bzbyEOKHsI3ZFilhXh6s4c+aM1SXcMk9PGD7cmLz46lx3K1dC06bwxx/W1uZqXDkHYh7lQMxm6pJiAwcOpEuXLvz111+sW7eOZcuW4eHhQfHixWnVqhUPPvggTz31FBUqVDDzsmKx5OSsbd1jd3MpKSlWl3DbXnkF7rkHnnrKmN9u/35o1gxmz4Z27ayuzjW4Qw7k9ikHYjZTR+zuueceRowYwaxZs4iNjWX//v2MGzeOtm3bsnHjRnr37k2LFi3MvKQ4gauNnbc3/O35GbmGGjVqWF2CKVq2hOhoqFPH2D93zpga5fPPjXvw5MbcJQdye5QDMZupjd3rr7/Oyy+/zFdffUVsbCx33303r7/+OrNmzeL48ePExMQwZ84cMy8pTuBqY1ekCHh4WFuLK1i5cqXVJZgmMBB++w2efNLYt9lgwAB49FE4fdrKypyfO+VAbp1yIGYzdYJiAJvNRnR0NGlpabRs2dLMU7uksLAwwsLCSE9PZ9++fW4xQXFKSgpF/nczXXJyMtWq+XH8ONx1Fxw9anFxLiApKcnlM/B3GRkwYgR8/HHWaF3FisZas5rM+NrcMQeSd8qBOPUExQAeHh7ce++9aur+p1evXsTExBAdHW11KXaTfcRObm706NFWl2A6T0/48EPjQYo77zRei401Pq4dP94YyZOc3DEHknfKgZjN9BE7uTZ3W1Jsz549ANSoURNfX0/S06FxY+OeK8nf/voLOnaE9euzXuvcGcaOVfMvIvJ3Tj9iJ+7P09OTWrVqUatWLS5fNpo60H+0c8vdlxAqXx5Wr4bevbNemz4dGjaETZssK8vpuHsOJHeUAzGbGju5LdmnOvHzs64OV/LGG29YXYLd+fjAqFEwd25Ww79/PzRvDsOG6alZyB85kJtTDsRsauwkz9LS0hgwYAADBgzgzJm0zNc1Ypc7s2fPtroEh+nYEbZsMT6mB2P5sY8+gocegmPHrK3NavkpB3J9yoGYTY2d5Nnly5f5/PPP+fzzzzl37nLm62rscue+++6zugSHuuce4367fv2ypsNZtQrq1oUffrC2NivltxzItSkHYjY1dnJbsk+arsYudw4fPmx1CQ7n4wNDhsCvvxrT4oCxvuwzz0DXrsbkxvlNfsyB/JNyIGZTYye3Jfs9dmrscsfLy8vqEizzwAOwfTs8+2zWa1OnQq1aEB5uWVmWyM85kCzKgZhNjZ3cFo3Y5V3ZsmWtLsFSd9wB338PkyfD1Sf7jx83Vqt45ZX8M3qX33MgBuVAzKbGTm6LRuzyzp0nq84tDw/jI9hdu6Bt26zXp0yB2rXhp5+sqsxxlAMB5UDMp8ZObotG7PLu6aeftroEp1GhAvz8M0ycmDV699df8Pjjxse1f/1lbX32pBwIKAdiPjV28v/t3XtcVHX+P/AXDjAggqAiiHEzL4FXREUJL20baalZW9m3Xbvs1mYrtWS7tWb9tF0DMe+Cprner5lSZqxKiYqihXgjb2iaqDhyUUBABxjO748TAyNgM3CYz8zwej4e5yEczsy8cV+effc553w+TcLGznRLly4VXYJFsbMD/vKXuqN3W7cCQUFAfDz0k2DbEuaAAOaAlMclxZpZQkICEhISoNPpkJWVZRNLiul0Ohw9ehQA8M03/fGf/8g3//7vf8DIkSIrI2snScDGjcA77wC5uTX7Bw4Eli0D+vUTVxsRUXPgkmJWZtKkSTh9+rRN3UehUqkwcOBADBw4EGVlNU90ccTOOFxCqGF2dsCLLwJnzgCvv16zPz1dnuR48mSgqEhcfUpiDghgDkh5bOyoSfjwhOkmT54sugSL166dPEKXmgoEB8v7dDp5mbJu3eR78qz98ixzQABzQMpjY0cmKy8vx6effopPP/0URUVcUsxUS5YsEV2C1YiIAI4dA/7zH8DJSd6XlyeP5g0aBBw4ILa+pmAOCGAOSHls7MhkFRUVeO+99/Dee+/h9m0uKWaq0aNHiy7Bqjg6Ah9+CJw9K689W+3oUWDoUPnS7ZUr4uprLOaAAOaAlMfGjpqk9lOxLi7i6rAmJ06cEF2CVfL3BzZvBvbuBfr2rdm/cSPQo4fc/FnT/XfMAQHMASmPjR01SVlZzdetW4urw5q0b99edAlWbfhwICMD+OwzoPqv8s4d4JNPgAcflO/D02rF1mgM5oAA5oCUx8aOmqS6sVOrAS55aBxnZ2fRJVg9lQp44w3g/HkgOhpwcJD3FxTIT8527y6vQWvJD1gwBwQwB6Q8NnbUJNWNHUfrjJeVlSW6BJvh4SGP0J07B/zxj/J0KQCQnS0vWdavH7Bjhzw/nqVhDghgDkh5NtnYLV68GIGBgXByckJoaChSU1Pve/zWrVsRHBwMtVqN4OBgJCYmGvx827ZtePzxx9GhQwfY2dnh+PHjzVm+VWFjZ7pHH31UdAk2JzAQWLdOfqCi9iTZP/0EjBkjP127e7dlNXjMAQHMASnP5hq7zZs3Izo6GlOnTsWxY8cwdOhQjBo1CtnZ2fUef+jQIYwfPx4TJkzAiRMnMGHCBDz//PP44Ycf9MeUlpbi4YcfxsyZM831a1iNO3fkP/nghPFWr14tugSb1a+fvALKnj3ydCjV0tLk5crCw4GdOy2jwWMOCGAOSHk2t6RYWFgY+vfvbzA3UFBQEMaNG4fY2Ng6x48fPx7FxcX43//+p983cuRIeHh4YOPGjQbH/vLLLwgMDMSxY8fQz8S1jZReMkQknU6nHwWNjByKigoV+vWT5xsjshSSBGzbBnz0kbySRW2DBgH/7/8BTzxRc/mWiEgELil2H+Xl5cjIyEBkZKTB/sjISKSlpdX7mkOHDtU5/vHHH2/weGNptVoUFxcbbLZCpVJhxIgRiIgYgYoK+YkJXoo1HpcQMg87O+APfwAyM+VpUnr2rPnZjz8Co0fLa9B+/TVQVWX++pgDApgDUp5q+vTp00UXoZQbN25g9uzZePvtt+Hr66vff+zYMRw8eBBvvfVWnde8//77+NOf/oTevXvr9126dAnr16/HRx99ZHBsYWEhFixYgIkTJ8Lb2/u+tcyYMQMjR47EzJkzMXPmTMybNw+APEK4Zs0aDBgwADNnzsSwYcMwY8YMdO3aFV988QV0Oh1++uknnDx5EpIkYePGjQgKCsKnn36qP7ZXr15YsWIFnJyc8MMPP+DChQsoKSnB9u3b4e/vj3nz5umPHTBgABYuXIj27dtjz5490Gg00Gg0+P7779G+fXssWrRIf+zDDz+MuLg4+Pn54ZtvvsHt27dx8eJFHD58GGq1GitWrEBoaKi+7mnTZmH//ocBAJ06FcHXdy8yMzOh0+mwefNmPPTQQ5g9e7b+/fv06YPly5fD2dkZhw4dwsWLF1FcXIwdO3bggQcewPz58/XHDhw4EAsXLoSnpye+++475ObmIicnB3v37oWHhwfi4+MRERGBTz75BBEREYiJiUFAQAC+/vprlJSU4MKFC/jxxx+hVquxatUqhISEIC4uTv/+PXr0wIYNGwAAx48fx6lTp1BRUYEtW7age/fuBnX37dsXn3/+OVxcXHDw4EFcvnwZt27dQlJSEjp37mxQ96BBgzBv3jx07NgRycnJyM/Px9WrV7Fv3z64u7tjyZIl+Oc//4nY2Fj9awIDA5GYmIiysjJkZWXhyJEjsLe3x+rVq9GvXz+DuoOCgrBu3Tq0atUKR48exdmzZ6HVavHll1+ia9eumDNnjv7Yfv36YenSpXB1dUVqaiqys7Nx69Yt7Ny5E97e3liwYIH+2LCwMMydOxfe3t7YtWsXCgoKcOXKFaSmpsLNzQ2fffYZBg8ejJiYGP1runTpgq1bt0Kr1eLMmTM4evQoVCoV1q5di969e2PWrFn6Y4ODg7FmzRrY29vjyJEjOHfuHO7cuYPExER06dIFc+fO1R8bEhKCJUuWwM3NDfv378fVq1eRn5+P3bt3w8vLy6DuIUOGYPbs2fDx8UFSUhIKCwtx+fJlHDhwAG3atMGyZcsweHAYtm+PxYYNw5CV9SV0um7Iz5f/g+T6dWDTJmDdOi0OHUrBiBEdMXduzb+1nj17YvXq1XBwcEB6ejrOnz+PsrIyfPXVVwgICDCou3///khISICHhwdSUlKQk5ODvLw8JCcnw9PTEwsXLtQfGx4ejrS0NPj6+mLHjh0oKirCpUuXcOjQIbRu3RrLly+3iXPEjBkz0L17d2zcuBFVVVXIzMzkOeKec0SPHj2Qk5NjcI4IDw83+LfGc0TzniPCwsIMzsldu3bFli1bUF5ejtOnT+vvq1+/fj169uxp8G9NiXPErl27kJSUhClTpkCtVpvU99RLsiHXrl2TAEhpaWkG+2fMmCH16NGj3tc4ODhIGzZsMNi3bt06Sa1W1zn20qVLEgDp2LFjv1nL3bt3paKiIv125coVCYBUVFRkwm9kmcrLy6X4+HgpJiZeAsolQJLGjhVdlfWIiYkRXUKLptNJUmKiJIWESJJ8wbZm8/aWpNhYSbp1q/nrYA5IkpgDkqSioiJF+wObuhTboUMHqFQqaDQag/25ubnw8vKq9zXe3t4mHW8stVoNNzc3g81WlJeXIyoqCh98EAVAXiuWl2KN9+KLL4ouoUVr1QoYN06e5Pibb+RlyappNMCUKYCvrzwfXgPPXCmCOSCAOSDl2VRj5+joiNDQUCQnJxvsT05ORnh4eL2vGTJkSJ3jd+/e3eDxVD82dsbbv3+/6BII8j14o0cD+/cDhw/L9+NVP0hRUiLPj9eli7w+bWqq8k/SMgcEMAekPJtq7ABg8uTJWL58OVasWIEzZ87gnXfeQXZ2NiZOnAgAeOmllzBlyhT98X//+9+xe/duxMXF4ezZs4iLi8N3332H6Oho/TE3b97E8ePHcfr0aQDAuXPncPz48TojfS0ZGzvjBQYGii6B7hEWBnz5JZCVBbz5JuDkJO/X6YAtW4Bhw4CQEGD5csNl9JqCOSCAOSDl2VxjN378eMyfPx///ve/0a9fP+zfvx9JSUnw9/cHAGRnZ+P69ev648PDw7Fp0yasXLkSffr0wapVq7B582aEhYXpj9m+fTtCQkLw5JNPAgBeeOEFhISE4LPPPjPvL2fB2NgZr7KyUnQJ1ICuXYHFi+VLsNOmAbXvyDhxAnj9dfky7XvvARcvNu2zmAMCmANSns3NY2epbGkeu9LSUrRp0+bX70oAuGD6dPn/COm3bdq0CS+88ILoMsgIWq08krdoEVBrznK9xx4D/vpXYOxYwNHRtPdmDghgDojz2JGF4soTxgsNDRVdAhlJrZbXoD18WJ77bsIEwwYuORl47jl5FO/994ELF4x/b+aAAOaAlMfGjhTBS7HG++qrr0SXQI0wcCCwZo18mTYmRn6wolpuLjBrFtCtG/C73wFr18oPYNwPc0AAc0DK46VYM7GlS7GVlZW/TqgILF78OAB7rFwJvPKK6Mqsw927d+FUfXc+Wa2qKnlN2mXLgMRE4N5bpVxc5CdtX3oJGDECUKkMf84cEMAcEC/FkgWwt7fHk08+iQcffBKAPQCO2Jli9uzZoksgBbRqBfz+98AXXwBXrwJxcfLDF9VKS+URvt//HggIkOfHq71mLXNAAHNAyuOInZnY0ohdtRkz5AXWAXmi19GjxdZDJJokAQcPyg3dF18ARUV1jxkwQB7Fe/55w6duiahl4ogdCVdRUYFVq1bh8OFVACoAcMTOFFz023bZ2QEREfLlWY1Gbu5Gjza8DHvkCPD220CnTlV49FFg6VIgL09czSQWzwekNI7YNbOEhAQkJCRAp9MhKyvLJkbs6pvu5NAhYPBgkVVZD41GA29vb9FlkBnl5gIbN8ojeUeP1v25SiU/dDF+PPD000C7duavkcTg+YA4YmdlJk2ahNOnTyM9PV10Kc2KI3bG++abb0SXQGbWsSPw97/L69NmZgIffAB4etZcp9Xp5KlTXntNvjz7+OPAkiXAtWsCiyaz4PmAlMbGjhTBxs54/fr1E10CCdSrF/DJJ8COHVnIyJBXsQgIqPl5ZSWwezfwt78BDzwADBokT69y+rTy69WSeDwfkNLY2JEi2NgZr6CgQHQJZAFu3ixA//7y07QXL8orW0yeDPj5GR6Xng5MnQr07Al07w7885/yAxo6nZi6SVk8H5DS2NiRIrjyhPHu3LkjugSyALVzYGcnj8zNmQP88ot8H960aUDfvoavuXABmD1bfkCjUyfgz38GtmwBCgvNWzsph+cDUhobO1IER+yM161bN9ElkAVoKAd2dkBICDB9OnD8OHDpEjB/PvDII4ZP1+blAStXytOmdOgADB0qX7I9doyXbK0JzwekNDZ21GT29oCDg+gqrEdKSoroEsgCGJuDgAD5wYs9e4AbN4DVq4FnnjH8jymdDjhwQL5k278/4OMDvPoqsGmT/EQuWS6eD0hpnO7ETGxpguLKykokJiYiOhrIyXkabm729U7ESvUrKipC27ZtRZdBgjU1B3fvAvv3A//7n7ydO9fwsX36AI8+Km/DhgGuro3+WFIYzwfE6U5IOHt7ezz33HOwt38OgD0vw5po0aJFoksgC9DUHDg5AZGRwLx5wNmz8gMYCQnyhMj3/ps8eVI+bvRoeY68iAj5Hr7UVKC8vEllUBPxfEBK44idmdjSiF01T08gPx/o0gX4+WfR1RBRtbt35abt++/lLSOj4fvuXFyAhx+WR/KGDpUf4uCa9ETmwxE7Eq6yshJbtmxBcfEWAJUcsTMRlxAioHlz4OQEPPYYMHOmPF1Kfj6wdSvw5pvylCm1lZbK8+Z9+CEwfDjQtq3c5E2dCuzaBdy+3WxlEng+IOVxxM5MbGnE7t4lxQYNcsEPPwgtyaoUFhbC3d1ddBkkmMgcXLlSM5r3/ffA9esNH9uqlfyU7tCh8iXcwYOBzp3NV6ut4/mAOGJHFocjdqZZt26d6BLIAojMga8v8MorwNq18rJl588DK1bI+x580PDYqir5Uu78+cCzz8qrYfj5ydOszJ0LpKXJl36pcXg+IKXZiy7A1iUkJCAhIQE6G54mnpMTm2b48OGiSyALYCk5sLMDunaVt1dflffl5Mj36O3fL/+ZmWn4mitX5G3LFvl7Bwd5VG/w4JotIEB+b7o/S8kB2Q5eijUTW74U+9xzLvjiC6ElWZWvvvoK48aNE10GCWZNObh5U17G7PBhefvxR6Ck5P6vadcOCA013Njs1WVNOaDmoXR/wBE7ajJeijWNo6Oj6BLIAlhTDtq1A8aMkTdAnhD51KmaRu/QIXnKldpu3gSSk+Wt9vv072/Y7AUGtuxmz5pyQNaBjR01GRs703h5eYkugSyANedApZInPe7TB/jrX+V9t27JI3nVI3oZGfJKGbXdvAl89528VfPwkJu9/v3l9+vbF+jRA2gp/Y4154AsExs7ajI2dqbJyMhAaGio6DJIMFvLgYcH8Pjj8gbI8+bl5ABHjshNXvV2b7N361bN07nVHByAoKCa5rG64fPysr3RPVvLAYnHxo5M5ujoiHffXYk5cwDAkY2diZ566inRJZAFsPUc2NnJ06J07gxU/6rVzV7tRi8jA9BoDF9bUSGvlnHypOF+T8+aRq9nTyA4WG4ArXm2EFvPAZkfGzsymYODAyIiXvm1seOInak+//xzfPjhh6LLIMFaYg5qN3tjx9bsz8kBTpyoaeZOnpTv2ausNHx9Xl7d0T0A6NRJbvLu3Tp0aP7fqalaYg6oefGpWDOxpadiAWDDBuCPf5S/XrAAePttsfUQkW3RaoEzZwybvRMngNxc49/D01Me0atu9IKC5JU3HnhAnniZyBLwqVgSrrKyEocP7/r1u8fRujVjZIoZM2bwv9CJOfgNajXQr5+81XbjhtzknTkDnD4t/3nqFFBQUPc98vLkbf9+w/1OTvK8fd27y1u3bjV/duxo3vv4mANSGkfszMSWRuzuncdu/XoXvPii0JKsyp07d+Ds7Cy6DBKMOVBWXp7c6FU3e9Vf32+5tPq4udVt9rp2Bbp0kS/tKt30MQfEETuyOFx5wjSLFy/Gu+++K7oMEow5UJanJzB8uLzVdutWTaN37hyQlSUvoXbhgvyQxr2Ki+UneY8cqfuzNm3kBi8wUP6zegsMlCdfbkx/xhyQ0tjYUZPx4QnTjK191zi1WMyBeXh4AOHh8labTgdkZ8uNXnWzV/315cvyGrn3Kimp/2ndaj4+DTd9nToB9vX8Py5zQEpjY0dNxsbONEePHkW3bt1El0GCMQdiqVRy0xUYWDP3XjWtFrh4sabRu3gRuHRJ/vOXX+of6QPkp3tzcuTl1+r7vAceAPz8DLdffsnBhAnd4OcHuLoq/mtSC8TGrpklJCQgISEBOp1OdCnNhreHmMbT01N0CWQBmAPLpVbLT9AGBdX9mU4HXLtW0+jVbvouXqw7AXPt112+LG+GhiMuTv7K3b2m4fP3l/+snh6mc2d5RJC3vtBvYWPXzCZNmoRJkybpb460RU5OoiuwLmq1WnQJZAGYA+ukUtU0X/fezwcApaXyqF51o3fxonzJt3rLz2/4vQsL5a2hS70A0LatYbNX39axI6dzacnY2FGTsbEzzc8//4yHH35YdBkkGHNgm1xc5FUxevas/+elpcCVKzWN3jffnIS7ex/991euNHypFwCKiuTt9OmGj7G3l+/pq93sdeoEeHsbbp6ecqNKtoWNHZnM0dERffvG48QJAHBkY2eiESNGiC6BLABz0DK5uAAPPSRvABAZ6Q4/v5qf63Ty5dzqRu/atfo3rbbhz6islBvEK1fuX0urVnJz5+VVt+m7d3N3t711em0VGzsymYODAzp3nvRrYyffj0LGW7t2LaZOnSq6DBKMOSCgbg5UKvleOh8fYPDg+l8jScDNmzVNXk5O/c1fXt79P7uqSm4iqyd9vh9HR7kB7NhR3jw9G946dpSnhmEjKAYnKDYTW5qgGAAefRTYs0f+uqSEN/SaoqqqCq14A0yLxxwQ0Lw50GrlCZqvXQM0mprtxg3D7zWa+1/+bQy1uuGmr/b3HToA7dvLI4It9bIwJygm4XQ6HW7cSP31u6FQq1vov8ZGiomJ4RJCxBwQgObNgVotz6EXEHD/4yRJnsj53mbv3mbw+nV56bb65vi7l1YLXL0qb8aws5Obu3bt5Eavffuar+/9s/bXbm4cGbwXR+zMxJZG7GovKdaqVQl0Og7XmYIjNQQwBySzthzodHITmJtbsxZvQ1turvwUcHPO9qVSyU3evU2fh4fcKFb/Wb3V/t5SLhdzxI4sCh+cMF1sbCzvrSLmgABYXw5UKvnyaYcOxh1fVSVP4dJQ85efL98vWFBQ82dhoTyKaAydrua9GvO7tG17/+bvfvucnS2jMbwXR+zMxFZH7Nq1K0FBAUfsTJGdnQ2/2o/BUYvEHBDAHNRHp5Obu9rNXu2vG9p3+7Z567S3lxtDN7earfb39/tZ7e+12mK4u3PEjiyEo6PoCqzP3r178dJLL4kugwRjDghgDuqjUtVcVjVFeXlNo1dUJF8yrp70ubCw7vf37jPm3sHaKitrGsymqG8N4Sa9n7JvRy0NL8Wa7sEHHxRdAlkA5oAA5kBJjo418+6ZSpLkGR6MaQCrvy8urtmKiuTGsjEqKxv3uoawsaMmYWNnOu39ZhalFoM5IIA5sBR2doCrq7z5+jbuPbRaw0bv3savoe9v3gQyM5X7XdjYNbOEhAQkJCRA15yPBQnEyYlNl9eYu3zJ5jAHBDAHtqT23H2mKC6W77dTivU8Y22lJk2ahNOnTyM9PV10KYpp1coBwCwAs+Dk5CC6HKvTv39/0SWQBWAOCGAOSHls7MhkkuQI4J8A/onWrfn0hKm2b98uugSyAMwBAcwBKY/TnZiJLU13cvNmzdNKo0YBSUli67E2d+7cmUIwGwAAGqpJREFUgbOzs+gySDDmgADmgJTvDzhiRyYrK9MBSAeQDkdH27x3sDnNmTNHdAlkAZgDApgDUh5H7MzElkbsTp0qRa9e8gTFzz5bgi1bOEExERFRY3DEjoS7e7fmaz4Va7oZM2aILoEsAHNAAHNAymNjRyarPQkjGzvTvf7666JLIAvAHBDAHJDy2NiRyWqP2HGCYtN9/fXXoksgC8AcEMAckPLY2Jlg8eLFCAwMhJOTE0JDQ5Gamiq6JCF4KbZpQkNDRZdAFoA5IIA5IOWxsTPS5s2bER0djalTp+LYsWMYOnQoRo0ahezsbNGlmV3tS7EcsTPdjRs3RJdAFoA5IIA5IOWxsTPS3Llz8Ze//AWvvfYagoKCMH/+fPj6+mLJkiWiSzM7jtg1TXljV4omm8IcEMAckPK4VqwRysvLkZGRgX/9618G+yMjI5GWllbva7RarcHizsXFxc1aozlVVjoAmAYAaN2aS4qZ6sEHHxRdAlkA5oAA5oCUp5o+ffp00UVYuhs3bmD27Nl4++234evrq99/7NgxHDx4EG+99Vad18yYMQMjR47EzJkzMXPmTMybNw8AMH78eKxZswYDBgzAzJkzMWzYMMyYMQNdu3bFF198AZ1Oh59++gknT56EJEnYuHEjgoKC8Omnn+qP7dWrF1asWAEnJyf88MMPuHDhAkpKSrB9+3b4+/tj3rx5+mMHDBiAhQsXon379tizZw80Gg00Gg2+//57tG/fHosWLdIf+/DDDyMuLg5+fn745ptvcPv2bVy8eBGHDx+GWq3GihUrEBoaik8+2YkzZyYCGIFHHy3BmTNrUVVVhczMTGRmZkKn02Hz5s146KGHMHv2bP379+nTB8uXL4ezszMOHTqEixcvori4GDt27MADDzyA+fPn648dOHAgFi5cCE9PT3z33XfIzc1FTk4O9u7dCw8PD8THxyMiIgKffPIJIiIiEBMTg4CAAHz99dcoKSnBhQsX8OOPP0KtVmPVqlUICQlBXFyc/v179OiBDRs2AACOHz+OU6dOoaKiAlu2bEH37t0N6u7bty8+//xzuLi44ODBg7h8+TJu3bqFpKQkdO7c2aDuQYMGYd68eejYsSOSk5ORn5+Pq1evYt++fXB3d8eSJUtw7do1JCcn618TGBiIxMRElJWVISsrC0eOHIG9vT1Wr16Nfv36GdQdFBSEdevWoVWrVjh69CjOnj0LrVaLL7/8El27dsWcOXP0x/br1w9Lly6Fq6srUlNTkZ2djVu3bmHnzp3w9vbGggUL9MeGhYVh7ty58Pb2xq5du1BQUIArV64gNTUVbm5u+OyzzzB48GDExMToX9OlSxds3boVWq0WZ86cwdGjR6FSqbB27Vr07t0bs2bN0h8bHByMNWvWwN7eHkeOHMG5c+dw584dJCYmokuXLpg7d67+2JCQECxZsgRubm7Yv38/rl69ivz8fOzevRteXl4GdQ8ZMgSzZ8+Gj48PkpKSUFhYiMuXL+PAgQNo06YNli1bhrCwMMTGxhr8W9uyZQvKy8tx+vRpHD9+HACwfv169OzZ0+DfWs+ePbF69Wo4ODggPT0d58+fR1lZGb766isEBAQY1N2/f38kJCTAw8MDKSkpyMnJQV5eHpKTk+Hp6YmFCxfqjw0PD8cbb7yBIUOGYMeOHSgqKsKlS5dw6NAhtG7dGsuXL7fqc0Tturt3746NGzfyHNHAOWL79u1wdXU1OEeEh4cb/FvjOcK2zxG7du1CUlISpkyZArUCl8E4QbERcnJy0LlzZ6SlpWHIkCH6/Z988gnWrl2Ls2fP1nlNfSN2vr6+NjFB8bJlwBtvyF//97/An/8sth5rU1hYCHd3d9FlkGDMAQHMAXGCYiE6dOgAlUoFjUZjsD83NxdeXl71vkatVsPNzc1gsxVlZVUATgE4BUfHKtHlWJ34+HjRJZAFYA4IYA5IeWzsjODo6IjQ0FAkJycb7E9OTkZ4eLigqsQpKbkDoBeAXmjV6o7ocqzOhx9+KLoEsgDMAQHMASmPjZ2RJk+ejOXLl2PFihU4c+YM3nnnHWRnZ2PixImiSzO7WleY+VRsI3AJIQKYA5IxB6Q0PhVrpPHjx6OgoAD//ve/cf36dfTq1QtJSUnw9/cXXZrZcbqTpqnvYRtqeZgDApgDUh5H7Ezwt7/9Db/88gu0Wi0yMjIwbNgw0SUJUXvEjhMUm27NmjWiSyALwBwQwByQ8tjYkcnY2DXNI488IroEsgDMAQHMASmPjR2ZrPalWEdHcXVYq/Pnz4sugSwAc0AAc0DKY2NHJuOIXdM4OzuLLoEsAHNAAHNAyuPDE2SyigoHAP8AALRpwyXFTNW+fXvRJZAFYA4IYA5IeRyxI5OVlzsC+BTAp3B15bVYU1UvT0MtG3NAAHNAymNjRyarfY8dL8WabsyYMaJLIAvAHBDAHJDy2Ng1s4SEBAQHB2PgwIGiS1HM3btVAH4B8AuXFGuE5cuXiy6BLABzQABzQMqzkyRJEl1ES6D0Ir8ihYaW4ujRNgCA27dL0KaNi+CKiIiIrJPS/QFH7MhktZ+KtbMTV4e14hJCBDAHJGMOSGkcsTMTWxqx69q1FD//LI/YlZSUwMWFI3amuHv3Lpx4c2KLxxwQwBwQR+zIAtQesSPTLVq0SHQJZAGYAwKYA1IeGzsyGRu7phk3bpzoEsgCMAcEMAekPDZ2ZLLa052Q6TIyMkSXQBaAOSCAOSDlsbEjk3HErmm8vb1Fl0AWgDkggDkg5XFJMTKJJAHl5fYA/oaOHQF7e0bIVPw7I4A5IBlzQErjiB2ZRB6tUwNIQFBQAtRqteCKrM+lS5dEl0AWgDkggDkg5bGxI5PUvgzLJ/QbZ9iwYaJLIAvAHBDAHJDy2NiRSeQHJyQAeQDywGkQTbdhwwbRJZAFYA4IYA5IeZyg2ExsZYLiy5eBgIBSAJyguLF0Oh1UKpXoMkgw5oAA5oA4QbHVSUhIQHBwMAYOHCi6FEXwidimi42NFV0CWQDmgADmgJTHETszsZURu5Mngb59OWJHRESkBI7YkVCcnLjpuOg3AcwByZgDUhobOzIJG7ume/nll0WXQBaAOSCAOSDlsbEjk/Aeu6b7/vvvRZdAFoA5IIA5IOWxsSOTcMSu6bp37y66BLIAzAEBzAEpj2uZkEnkxs4ewMsYMIDL4TTGnTt3RJdAFoA5IIA5IOXx/5XJJDVLiq3CK68AXFHMdAUFBaJLIAvAHBDAHJDyeCmWTFL7Uiybusbp27ev6BLIAjAHBDAHpDw2dmSSmiXFSgGUckmxRtixY4foEsgCMAcEMAekPE5QbCa2MkHxnDnAP/7BCYqboqysDK1btxZdBgnGHBDAHBAnKCbB+FRs082dO1d0CWQBmAMCmANSHkfszMRWRuw++giYMYMjdkRERErgiJ2VSUhIQHBwMAYOHCi6FEVwguKm4xJCBDAHJGMOSGkcsTMTWxmxe/ttYNEijtg1RV5eHjw9PUWXQYIxBwQwB8QROxKM99g13bZt20SXQBaAOSCAOSDlsbEjk7CxazpbuSxPTcMcEMAckPK48gSZRL7HTgXgWTzxBKBSqQRXZH2uX78uugSyAMwBAcwBKY+NHZlEHrFzArAFK1YATk6CC7JCOp1OdAlkAZgDApgDUh4vxZJJal+KZVPXOAEBAaJLIAvAHBDAHJDy2NiRSWpPd8LGrnFSU1NFl0AWgDkggDkg5bGxI5PII3alAOzg5GSH0tJSwRVZn//7v/8TXQJZAOaAAOaAlMfGjkzCp2KbbvHixaJLIAvAHBDAHJDyOEGxmdjKBMVBQcDZs5ygmIiISAmcoJiE4ohd03EJIQKYA5IxB6Q0jtiZia2M2HXqBGg0HLFriuLiYqvOACmDOSCAOSCO2FmdhIQEBAcH28zs4rWfiqXGWbVqlegSyAIwBwQwB6Q8NnbNbNKkSTh9+jTS09NFl6IIXoptut///veiSyALwBwQwByQ8tjYkdEkqbqxU6Ft2yfwxBNPcEmxRjh79qzoEsgCMAcEMAekPC4pRkarqJCbO8AJffp8i2+/FV2RdeI9iQQwByRjDkhpHLEjo3HVCWW0a9dOdAlkAZgDApgDUh4bOzIa14lVxsmTJ0WXQBaAOSCAOSDlsbEjo9U0dqVISnKBi4sLlxRrhFGjRokugSwAc0AAc0DKY2NHRqt9KVanK0NZWZm4YqzYihUrRJdAFoA5IIA5IOXZVGMnSRKmT58OHx8fODs7Y8SIETh16tRvvm7x4sUIDAyEk5MTQkNDkZqaavDzZcuWYcSIEXBzc4OdnR0KCwub61ewaJzqRBkffvih6BLIAjAHBDAHpDybauxmzZqFuXPnIj4+Hunp6fD29sZjjz2G27dvN/iazZs3Izo6GlOnTsWxY8cwdOhQjBo1CtnZ2fpjysrKMHLkSHzwwQfm+DUsFhs7ZXAJIQKYA5IxB6Q0m1lSTJIk+Pj4IDo6Gu+//z4AQKvVwsvLC3FxcXjjjTfqfV1YWBj69++PJUuW6PcFBQVh3LhxiI2NNTh27969eOSRR3Dr1i24u7ubVJ8tLCl28CAQEQEAXFKsKbRaLdRqtegySDDmgADmgLikWIMuXboEjUaDyMhI/T61Wo3hw4cjLS2t3teUl5cjIyPD4DUAEBkZ2eBrWjKO2CljwYIFoksgC8AcEMAckPJsZoJijUYDAPDy8jLY7+XlhcuXL9f7mvz8fOh0unpfU/1+jaXVaqGt9bRBUVERALkzt1Y3b1Z/VfMkbHFxMXQ6nZB6rFVkZKRV54CUwRwQwBxQTV+g1AVUq23s1q9fb3B59dtfl0Gws7MzOE6SpDr77tWY1/yW2NhYfPzxx3X2+/r6Nul9LY2Pj4/oEoiIiKxeQUEB2rZt2+T3sdrGbuzYsQgLC9N/Xz06ptFo0KlTJ/3+3NzcOiNy1Tp06ACVSlVndO5+rzHWlClTMHnyZP33hYWF8Pf3R3Z2tiL/w93PwIEDkZ6e3qyfwc9pnOLiYvj6+uLKlSvNfq+lLf292drnMAf8HIA54OfIioqK4Ofnp9gqJFbb2Lm6usLV1VX/vSRJ8Pb2RnJyMkJCQgDI99Dt27cPcXFx9b6Ho6MjQkNDkZycjKefflq/Pzk5GU899VST6lOr1fXeENu2bdtm/wesUqnM8oAGP6fx3NzcmIMW/jkAc8DPkTEH/BwAaNVKmccerLaxu5ednR2io6MRExODbt26oVu3boiJiUHr1q3x4osv6o979NFH8fTTTyMqKgoAMHnyZEyYMAEDBgzAkCFDsGzZMmRnZ2PixIn612g0Gmg0Gly4cAEAkJmZCVdXV0U7bCVNmjSJn2PBn2Mutvb3ZmufYy629vdma59jLrb292Zrn6MoyYZUVVVJ06ZNk7y9vSW1Wi0NGzZMyszMNDjG399fmjZtmsG+hIQEyd/fX3J0dJT69+8v7du3z+Dn06ZNkwDU2VauXGl0bUVFRRIAqaioqLG/HtkA5oAkiTkgGXNAkqR8DmxmHjtLp9VqERsbiylTpnDOohaMOSCAOSAZc0CA8jlgY0dERERkI2xmgmIiIiKilo6NHREREZGNYGNHREREZCPY2JnB4sWLERgYCCcnJ4SGhiI1NVV0SWRm06dPh52dncHm7e0tuixqZvv378eYMWPg4+MDOzs7fPXVVwY/lyQJ06dPh4+PD5ydnTFixAicOnVKULXUXH4rB6+88kqd88PgwYMFVUvNJTY2FgMHDoSrqys6duyIcePG4dy5cwbHaLVavPXWW+jQoQNcXFwwduxYXL161aTPYWPXzDZv3ozo6GhMnToVx44dw9ChQzFq1ChkZ2eLLo3MrGfPnrh+/bp+y8zMFF0SNbPS0lL07dsX8fHx9f581qxZmDt3LuLj45Geng5vb2889thjuH37tpkrpeb0WzkAgJEjRxqcH5KSksxYIZnDvn37MGnSJBw+fBjJycmorKxEZGQkSktr1l+Pjo5GYmIiNm3ahAMHDqCkpASjR482bU12RSZNoQYNGjRImjhxosG+hx56SPrXv/4lqCISYdq0aVLfvn1Fl0ECAZASExP131dVVUne3t7SzJkz9fvu3r0rtW3bVvrss89ElEhmcG8OJEmSXn75Zempp54SVBGJkpubKwHQz51bWFgoOTg4SJs2bdIfc+3aNalVq1bSzp07jX5fjtg1o/LycmRkZCAyMtJgf2RkJNLS0gRVRaKcP38ePj4+CAwMxAsvvICLFy+KLokEunTpEjQajcH5Qa1WY/jw4Tw/tEB79+5Fx44d0b17d7z++uvIzc0VXRI1s6KiIgDQr2CVkZGBiooKg3OCj48PevXqZdI5gY1dM8rPz4dOp4OXl5fBfi8vL2g0GkFVkQhhYWFYs2YNdu3ahc8//xwajQbh4eEoKCgQXRoJUn0O4PmBRo0ahfXr12PPnj2YM2cO0tPT8bvf/Q5arVZ0adRMJEnC5MmTERERgV69egGQzwmOjo7w8PAwONbUc4LNrBVryezs7Ay+lySpzj6ybaNGjdJ/3bt3bwwZMgQPPvggVq9ejcmTJwusjETj+YHGjx+v/7pXr14YMGAA/P398e233+KZZ54RWBk1l6ioKJw8eRIHDhz4zWNNPSdwxK4ZdejQASqVqk6nnZubW+e/0qllcXFxQe/evXH+/HnRpZAg1U9F8/xA9+rUqRP8/f15frBRb731FrZv346UlBQ88MAD+v3e3t4oLy/HrVu3DI439ZzAxq4ZOTo6IjQ0FMnJyQb7k5OTER4eLqgqsgRarRZnzpxBp06dRJdCggQGBsLb29vg/FBeXo59+/bx/NDCFRQU4MqVKzw/2BhJkhAVFYVt27Zhz549CAwMNPh5aGgoHBwcDM4J169fx08//WTSOUE1ffr06UoVTXW5ubnho48+QufOneHk5ISYmBikpKRg5cqVcHd3F10emck//vEPqNVqSJKErKwsREVFISsrC0uXLmUObFhJSQlOnz4NjUaDpUuXIiwsDM7OzigvL4e7uzt0Oh1iY2PRo0cP6HQ6vPvuu7h27RqWLVvGReFtyP1yoFKp8MEHH8DV1RU6nQ7Hjx/Ha6+9hoqKCsTHxzMHNmTSpElYv349vvzyS/j4+KCkpAQlJSVQqVRwcHCAk5MTcnJyEB8fj759+6KoqAgTJ06Eq6sr4uLi0KqVkWNxCj65Sw1ISEiQ/P39JUdHR6l///76R5up5Rg/frzUqVMnycHBQfLx8ZGeeeYZ6dSpU6LLomaWkpIiAaizvfzyy5IkyVOeTJs2TfL29pbUarU0bNgwKTMzU2zRpLj75aCsrEyKjIyUPD09JQcHB8nPz096+eWXpezsbNFlk8LqywAAaeXKlfpj7ty5I0VFRUnt2rWTnJ2dpdGjR5ucBbtfP4yIiIiIrBzvsSMiIiKyEWzsiIiIiGwEGzsiIiIiG8HGjoiIiMhGsLEjIiIishFs7IiIiIhsBBs7IiIiIhvBxo6IiIjIRrCxIyIiIrIRbOyIiIiIbAQbOyIiM9uyZYvoEojIRrGxIyIyo7S0NKxevVp0GURko9jYERGZ0fr16/HHP/5RdBlEZKPsJEmSRBdBRNQSVFZWolu3bjh16hRat24tuhwiskEcsSMiMpOdO3ciIiKCTR0RNRs2dkREZsLLsETU3HgplojIDEpKStC7d2+cP38e9vb2osshIhvFETsiIiOdOHECEyZMQEREBHbs2IFbt24hOjoab775Jp5++mkcP368wddu27YNY8aMqdPUHThwAC+99BKGDBmCb7/9FlVVVUhISMBbb72FiRMnIiIiAocPHwYAbNiwAW+++SbefPNNREREYM+ePc36+xKRFZKIiMgor776qlRRUSH95z//kTw9PaU//OEP0tWrV6U9e/ZIDg4OUlRUVIOvjYyMlA4fPmywT6fTSRMmTJCqqqqkd999V/L09JSioqKkPXv26I/561//KgUGBkoff/yxtHXrVv3+jz76SGrXrp2k0+mU/0WJyGpxxI6IyAg///wzfHx8YG9vj5ycHNy8eRMffPABOnfujIKCArRp0wZjxoyp97U3btzA5cuXERYWZrA/PT0dISEhsLOzQ05ODvLy8vDkk0/ikUce0R/j5uaGS5cuwcvLC88884x+v5eXF27evIm8vLzm+YWJyCrxHjsiIiOkpaXB2dkZISEh6NOnD9q1a4e9e/ca9doFCxbg5s2b+Pjjjw3279+/H35+fggICICvry969uyJnTt3Ghzz+OOP4/z58/j5559hZ2en3x8VFYUVK1aguLiY9+wRkR5H7IiIjBAeHo6QkBDk5+fjp59+wogRI4x+7fr16/Hiiy/W2T9s2DAEBATgwoULuHr1Kh577DGDn1dWViItLQ0jRowwaOoAYPfu3Rg6dCibOiIywMaOiMgEKSkpkCTJ6MYuKysLkiShR48eDR5T/RDEve+Znp6OkpKSOvszMzNx/vx5PPvss6aUTkQtABs7IiITpKSkQK1WY/DgwUYdb8zcdSkpKXB3d0dISIjB/upLvfc2dhs2bIBardY3dv/973+NK56IbB7H8ImITJCSkoLBgwfDycnJqOM3b978m/fi7d27F8OGDUOrVob/rZ2SkoIuXbrAz8/PYP+2bdswduxYeHh4IDMzE1evXjXpdyAi28UROyIiI2k0Gpw9e9boy7CHDx+Gn58fvL29GzzmzJkz0Gg0dd6zoqJCf3/dvfLz8zF8+HBUVVUhLi4OUVFRJvwWRGTL2NgRERnpxo0b8PLywvPPP2/U8cZchi0oKICPj4/BVCYAUFhYCBcXF/zpT3+q85qFCxfiyy+/xPjx4/Hqq6+iffv2xv8SRGTTON0JEVEzqKysRNeuXZGZmQlXV1fR5RBRC8EROyKiZrB7926EhYWxqSMis2JjR0TUDIy5DEtEpDReiiUiUlhpaSmCg4Nx4cIFODg4iC6HiFoQjtgRESksMTERTzzxBJs6IjI7NnZERArbt28fXnrpJdFlEFELxEuxRERERDaCI3ZERERENoKNHREREZGNYGNHREREZCPY2BERERHZCDZ2RERERDaCjR0RERGRjWBjR0RERGQj2NgRERER2Qg2dkREREQ2go0dERERkY1gY0dERERkI/4/Vd/g0dEOw/IAAAAASUVORK5CYII=\n", "text/plain": [ "Graphics object consisting of 2 graphics primitives" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "rmax = 20\n", "ymin = -0.01\n", "ymax = 0.05\n", "graph = plot(U(r), (r, 1, rmax), plot_points=300, thickness=2,\n", " xmin=0, xmax=rmax, ymin=ymin, axes_labels=[r'$r/m$', r'$m^2 \\; U(r)$'],\n", " axes=False, frame=True, gridlines=True, axes_pad=0)\n", "graph += line([(2, -0.1), (2, ymax)], color='black', thickness=1.5, linestyle='--',\n", " ymin=ymin)\n", "show(graph, xmin=0)" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAkwAAAGFCAYAAAAPa6wiAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAAPYQAAD2EBqD+naQAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi40LCBodHRwOi8vbWF0cGxvdGxpYi5vcmcv7US4rQAAIABJREFUeJzs3XdYFdfWwOHfERQVDXrtXWMndkUFa+wlzSRKzA1qYprRJGg+WzT3eqNGsTfUqLGlKPbeOyoaYlessWBDwQJSpM73x4QmIAdOmcOw3ufhYc6wZ2YdlsJiz569DYqiKAghhBBCiAzl0ToAIYQQQghbJwWTEEIIIUQmpGASQgghhMiEFExCCCGEEJmQgkkIIYQQIhNSMAkhhBBCZEIKJiGEEEKITEjBJIQQQgiRCSmYhBBCCCEyIQWTEEIIIUQmdFkwzZ07lypVqpA/f34aN26Mr6/vS9uvXbsWZ2dnHBwccHZ2Zv369am+3q9fPwwGQ6qP5s2bW/ItCCGEEMKG6K5g8vHxwdPTk1GjRnHq1ClatWpF165dCQwMTLe9n58f7u7ueHh4cObMGTw8POjVqxfHjx9P1a5Lly7cv38/6WPbtm3WeDtCCCGEsAEGvS2+26xZMxo1asS8efOS9tWuXZt33nmHCRMmpGnv7u5OWFgY27dvT9rXpUsXihYtyooVKwC1h+np06ds2LDB8m9ACCGEEDbHXusAzCkmJoYTJ04wYsSIVPs7derE0aNH0z3Gz8+PwYMHp9rXuXNnZsyYkWrfgQMHKFmyJEWKFKFNmzaMHz+ekiVLZhhLdHQ00dHRSa8TEhJ4/PgxxYoVw2AwZPWtCSGEECKLFEXh2bNnlC1bljx5TLuppquCKSQkhPj4eEqVKpVqf6lSpQgKCkr3mKCgoEzbd+3alZ49e1KpUiVu3LjBDz/8QLt27Thx4gQODg7pnnfChAn873//M/EdCSGEEMJUt2/fpnz58iadQ1cFU6IXe3AURXlpr05m7d3d3ZO269SpQ5MmTahUqRJbt27l3XffTfecI0eOZMiQIUmvQ0NDqVixIrdv3+aVV17J0vsROUNgYCAVK1bUOgxhAZJbfZP86ldYWBgVKlSgcOHCJp9LVwVT8eLFsbOzS9Ob9PDhwzS9SIlKly6dpfYAZcqUoVKlSly9ejXDNg4ODun2Pr3yyis5rmBSFIWQkBBA/R7LLcX0bd68mZEjR2odhrAAya2+SX71zxy/t3T1lFy+fPlo3Lgxu3fvTrV/9+7duLm5pXuMq6trmva7du3KsD3Ao0ePuH37NmXKlDE96BwgMjKSkiVLUrJkSSIjI7UOx2YNGzZM6xCEhUhu9U3yK4yhq4IJYMiQISxatIjFixdz8eJFBg8eTGBgIF9++SUAffr0SfWXxLfffsuuXbvw8vLi0qVLeHl5sWfPHjw9PQEIDw/n//7v//Dz8+PmzZscOHCAN998k+LFi9OjRw9N3qOwTek9hSn0QXKrb5JfYQzdTSsA6sSVkyZN4v79+9SpU4fp06fTunVrANq2bUvlypVZunRpUvs1a9YwevRorl+/TtWqVRk/fnzS2KSoqCjeeecdTp06xdOnTylTpgyvv/46Y8eOpUKFCkbHFBYWhpOTE6GhoTnullxERASFChUC1ALS0dFR44iEEEJY1NV56kf4TfW102tQ9z9QtqumYWWVOX/36rJgskVSMOnfuHHjGD16tNZhCAuQ3Oqb5DcddzaDwQ4KV1Nf31gGFydDl1NQ5DVtY8sCc/7u1dWgbyG01LdvX61DEBYiudU3yW86yr+Z+nX98WqP06NjOapgMifdjWESQit79+7VOgRhIZJbfZP8ZiIhHm6uhLgIKO6qdTSakR4mIcykRo0aWocgLERyq2+S3ww8PQe7XCH+OdgXglbrwclZ66g0IwWThXl7e+Pt7U18fLzWoWSbvb19Upe1vb38k8lIVFSU1iEIC5Hc6pvkNwOFa0LX0xDzFG6vhWN9ocPBXFs0yW8/Cxs4cCADBw5MGniWEzk4OKR6qlCk79GjR1qHICxEcqtvkt8M2OVLHvRdrAk88ofLM6Hpz9rGpREZwySEmdSvX1/rEISFSG71TfJrLAXiozNvplNSMIlMKYpCREQEERERyCwUGduyZYvWIQgLkdzqm+Q3Hae/h4e+6jxMT8/BmVHw8ABU/rfWkWlG5mGyEpmHyYyenIbQAKjYE/Lk1TaWFCIjIylYsKDWYQgLkNyaiaLA3wuhRCtwqq11NEkkv+k41h8e7IWo+5DXCYrUA+fhUKaj1pFliTl/90oPk8h5Lk6Fo/+GTa/C379AQqzWEQEwbdo0rUMQFiK5NZPIO/DnF7D1NTjcG0Ivah0RIPlNV/Nf4O2b8EE0vPcQ2u/JccWSuUkPk5VID5MZHfk33FrxzwsFCpaHumOgSh+b6nESQrwgIhA2VlK3DfagxENFd3XJDRvqcRL6IT1MQhjyAP/U+pF34finmvc4jRs3TpPrCsuT3FqAEgcocHuN5j1Okl9hDCmYhA7YRuH0xRdfWPV6wnoktxZkA4WT5FcYQ+ZhEtmjKBAfqc21E+LU66eRWDjdUQuns/+B10ZB5Q+tcqtu07oV9O/f3+LXEdYnuTWTuJf8zFDi1M+BqyFwJZR/H14bAU61LB6W5FfH4iLMdioZw2QluhvDFBcBqwppHJkQQgiRsbBIcPoMs/zulR4mC9PD0ih2dna8//77SdtCCCFEbiM9TFaSk3uY0qXlLTm/T9TxDiSk80UDoEDhGlDnByj/9j8DxC1v69atdO/e3SrXEtYluTWTiNuw9WVPw+UBDPBqX6g9FBwrWCUsya9+hYWF4VSsrPQwCQ0ZDGCv0fQCeezV66cq9fMACfBKTag3Fiq8a7VCKVGFKrW0+54Ii5Lcmol9BpNDGuwAA1TtD699D44VrRqW5FfH7M13d0eekhM53D//hF+pAS1XQ/cLUPF9qxdLAL6+vla/prAOya2FGOzU+Ziqfgpv/Q1N51u9WALJrzCO9DCJTNncxJWgTngHaqGkUY/Si3r37q3p9YXlSG4twGCvWY/SiyS/whjSwyRynqL1oEgDzXuUXjR37lytQxAWIrk1k3z/gsLVodoXmvYovUjyK4whg76tJCcP+rbJHiYhhBAiE7I0ihA2SJZX0C/Jrb5JfoUxpIfJSqSHSf/CwsJyXG6FcSS3+ib51S/pYRLCBi1dulTrEISFSG71TfIrjCEFk4V5e3vj7OyMi4uL1qEIC+vQoYPWIQgLkdzqm+RXGEMKJgsbOHAgAQEB+Pv7ax1KttnZ2dGtWze6desmS6O8xKVLl7QOQViI5FbfJL/CGDIPk8hU/vz52bp1q9Zh2DwZ26Vfklt9k/wKY0gPkxBm8q9//UvrEISFSG71TfIrjCEFkxBmcvbsWa1DEBYiudU3yW9aFxYuZEevXqxycWFtq1Yc+vprwm7c0DosTUnBJDIVERGBo6Mjjo6OREREaB2OzeratavWIQgLkdzqm+Q3rYf+/tTo3ZtOK1bQbuFCEuLj2ffZZ8RFRmodmmakYBJGiYyMJDIX/0cxxuLFi7UOQViI5FbfJL9pvb5gAa/26EGRatUoWqsWzceNI/L+fR4HBGgdmmZk4korkYkrhRBC5FTPbt1ic7dudNuwgSLVq2sdjtFk4kohbJAsr6Bfklt9k/y+nKIonJw0iRKNGuWoYsncpIfJSqSHSf+io6NxcHDQOgxhAZJbfZP8vpz/2LHcO3SIjr/+SsHSpbUOJ0ukhykHkZm+c4+ZM2dqHYKwEMmtvkl+M/bX+PHcPXCA9kuW5Lhiydykh8lKpIdJ//7++2+qVq2qdRjCAiS3+ib5TUtRFP4aP547e/fSfulSXqlUSeuQskV6mIRV5cmThzZt2tCmTRvy5JF/Mhn5888/tQ5BWIjkVt8kv2n9NXYsN7dswW3SJPIWLEhUcDBRwcHEPX+udWiakaVRRKYKFCjAgQMHtA7D5pUpU0brEISFSG71TfKb1lUfHwD29uuXan/zceN4tUcPDSLSnhRMOhUXFcW1Vat4eu0azX78EYPBoHVIuicLE+tXbspt+O3bBCxZQvG6dXPNL8bclF9jfXjhgtYh2By5v6IzcVFRXFq2jA0dOnBy0iSur1tHQkyM1mHlCrdu3dI6BGEhuSG34bdvc+yHH9jcrRvXfHy4vWeP1iFZTW7IrzCdLgumuXPnUqVKFfLnz0/jxo3x9fV9afu1a9fi7OyMg4MDzs7OrF+/PsO2X3zxBQaDgRkzZpg7bJOkKpQmTybm6VOznTsiIoISJUpQokQJWRrlJVq2bKl1CMJC9JzblIXSjY0bURIStA7J6vScX2E+uiuYfHx88PT0ZNSoUZw6dYpWrVrRtWtXAgMD023v5+eHu7s7Hh4enDlzBg8PD3r16sXx48fTtN2wYQPHjx+nbNmyln4bRku3ULLAg48hISGEhISY/bx6snLlSq1DEBaix9ymVygp8fFah6UJPeZXmJ/uphVo1qwZjRo1Yt68eUn7ateuzTvvvMOECRPStHd3dycsLIzt27cn7evSpQtFixZlxYoVSfvu3r1Ls2bN2LlzJ927d8fT0xNPT0+j40p8tPHR/ftmmVYg7vlz/l67loAlS4gNDX1p2/eOHMEuX75sXysiMpJSJUsC8ODhQxwLFsz2ufQsLi4Oe3sZFqhHespt+N27BCxezM3Nm8FggAx6lMq0akWradOsHJ029JRfkVpYWBjFypQxy7QCuvoXEhMTw4kTJxgxYkSq/Z06deLo0aPpHuPn58fgwYNT7evcuXOqW24JCQl4eHgwdOhQXnvtNaNiiY6OJjo6Oul1WFgYAOvatqWglQcYrm3RwuRzLK5dG4CtbdqYfC4hhI14yd/L9319WSUT7oocLtKMvaa6uiUXEhJCfHw8pUqVSrW/VKlSBAUFpXtMUFBQpu29vLywt7fnm2++MTqWCRMm4OTklPRRoUKFLLwTIYQQQtgSXfUwJXrxEXpFUV76WP3L2p84cYKZM2dy8uTJLD2aP3LkSIYMGZL0OiwsjAoVKvDugQMmdQtGBAXx9+rV3Ny+nefBwUYdI7fkrMNr0iSGDxumdRjCAvSQ22d37nBh4UICd+5Ub8NlMhqjRKNGNBw+3KhzFyhWDAcnJ3OEqQk95FekLywsjM/MNM+Wrgqm4sWLY2dnl6Y36eHDh2l6kRKVLl36pe19fX15+PAhFStWTPp6fHw83333HTNmzODmzZvpntfBwSHdxRztCxbE3oSCw+nVV2k0fDgNhgzhxqZNnPP2JvLBA3UsQgY/AO0LFMDOhIUl7RWF6H/ObV+ggEnx61nfTz+V741O6SG3RWvUoOXkyUQNG0bA4sVcXbECJT4+w6fi8r3yCsXr1LFylNrQQ35F+uzj4sx2Ll3dksuXLx+NGzdm9+7dqfbv3r0bNze3dI9xdXVN037Xrl1J7T08PDh79iynT59O+ihbtixDhw5l586dlnkjRsiTNy9V33uPt3bupNmPP1Lwnx4gLDBBZZ48eWjSpAlNmjSRpVFe4sV/R0I/9JTbAiVK0Hj4cN7evZsaH31Enrx5MeTy/9d6yq+wHF31MAEMGTIEDw8PmjRpgqurKwsWLCAwMJAvv/wSgD59+lCuXLmkJ+a+/fZbWrdujZeXF2+//TYbN25kz549HD58GIBixYpRrFixVNfImzcvpUuXpmbNmtZ9c+lILJyqvPWW0T1OWVWgQAH8/f3Nci49q/3PwHihP3rMbWLh5PzJJ0b1OOmZHvMrzE93f1a4u7szY8YMfvzxRxo0aMChQ4fYtm0blf5ZaTkwMJD79+8ntXdzc2PlypUsWbKEevXqsXTpUnx8fGjWrJlWbyFbMuxxElbz7NkzrUMQFqLn3KbX45Tb6Dm/wnx0Nw+TrUqch8kcc0EYIyE2lhubN/P02jUaDR0qa8lZwerVq+nZs6fWYQgLyE25jQoO5tKvv/IvZ2cqdemidThWkZvym9uY83ev7m7JCVWevHmp+u67ZjlXZGQkzs7OAAQEBFBQBkemq27dulqHICwkN+W2QIkSNEzxhG9ukJvyK7JPd7fkhPkpisKtW7e4desW0iGZsZSzxQt9kdzqm+RXGENuyVmJtW/JmVNERASFChUCIDw8HEdHR40jsk3h4eFJ3yehL5JbfZP86pc5f/dKD5MwK19fX/73v/+xY8cOrUOxupTL6Qh9kdzqm+RXGEN6mCzM29sbb29v4uPjuXLliq57mO7evcuaNWtwcHCgcuXKdMklA0aFEELYJulhykEGDhxIQECA7ucxiomJYd26dbz55pvkz59f63A0MW7cOK1DEBYiudU3ya8whhRMwiy2bdtG9erVefXVV7UORTOJk6MK/ZHc6pvkVxhDCiaRKYPBgLOzM87OzunO53T+/Hnu379Phw4dsnReRYHTp+GHH+Ddd+HLL2HKFLh3z1yRW9eaNWu0DkFYiORW3yS/whgyD5PIVMGCBblw4UK6XwsNDWXHjh189NFH2Nsb/88pLAz69IGNG9N+7fvvoV8/mDYNctKDK82bN9c6BGEhklt9k/wKY0jBJExy//59IiIiWLBgQdK+xHmb/vzzT0aPHp1mwd7bt6FzZ7h4Mf1zxsbCwoVw5AisXw81aljyHZjPnTt3aNCggdZhCAuQ3Oqb5FcYQwomYZIqVaowYMCAVPs2btxI8eLFadGiRZpi6flzePPN5GKpaFEYOxY6dYKnT9Uep1mz4NkzCAgAFxdYvhzeftta7yj75IFT/ZLc6pvkVxhDCiaRqcjISFxcXADw9/dPtTSKg4MDJV9Y6Ddv3rwUKFAgzX6A2bPhzBl1u2pV2LlT/ZzIxQX69oUePeDCBfXWXY8e8Pvv0Lu3+d+bOVWoUEHrEISFSG71TfIrjCGDvnVMURRiYmLMcp6AgAACAgJM+kvs0SMYP17dNhhg7drUxVKi6tXh2DHo1Svx+up4p61bs31pq/Dz89M6BGEhklt9k/wKY8jElVZizaVRFEXh0qVL7N+/n9DQUIYPH57m1lhWJMYOcOLECRo2bJju03KZ+e47dSA3wMcfw+LFL2+vKOqTc4nDowoUAF9faNw4y5e2ikePHlGsWDGtwxAWILnVN8mvfsnElTmIt7c3zs7OSbe0LElRFC5evMi8efNYtWoVwcHBxMTEkJCQYNJ5Ux6/du1a5s+fz6VLl7LU23TjBsyZo24XKAA//pj5MQYDzJ0L7u7q66goeOcdePAgK9Fbz7x587QOQViI5FbfJL/CGNLDZCWW7GFK2aMUHByMwWBIVcyMGjUqS4/8v+jp06cULVoUgO+//x4HBwcURaFkyZK8/vrr1KxZM9Mep/79k3uUvv8++dacMaKjoW1b9TYdQJs2sHcv2Nll480IIYTINaSHSQBpe5RCQkKS9lv6ugDBwcH4+Phk2uMUEQGrVqnbTk4wbFjWrufgAOvWQdmy6uuDB8HLK7vRW44sr6Bfklt9k/wKY8hTclYWExNj8kBsRVG4cuUKhw4dSiqSEvdndE1TbsvFxsaSN2/edOMAePjwIT4+PpQoUYLWrVtTvXr1VD1OGzaovUR580LPnnlxcsr6+KcyZcDHR+1dSkiA//4XOnSApk2z/bbM7ttvv9U6BGEhklt9k/wKY8gtOStJ7BYcMWJErl2cFqBVq5G0a5cv28f/8AMk/jFYrRqcOmU7s4HPmDEDT09PrcMQFiC51TfJr37JLTmRY7m5mXb8f/6T3Kt07RrY0h+GnTp10joEYSGSW32T/ApjyC05K/vuu+9MrnKjoqL4888/+fPPP4mLi8u0/dChQ00a9B0bG8uUKVMybWdvb0+zZs1wcXGhQIECAPzyS3JR89134OCQ9tZeVuTNq05i2bAhhIerA8m7dYP33jPptGZx8eJFnJ2dtQ5DWIDkVt8kv8IYUjBZWb58+ciXL/u3pBLP0bFjR1q2bImfnx/Hjh0jLi4uwzFM+fLlM6lgetkTcAaDAXt7e1xdXWnevHlSoZRoxw51bThQi5psTN+URrVq6ozhH3+svh4wQB3bVLy46ec2ReHChbUNQFiM5FbfJL/CGFIw5WAFChSgXbt2uLq6GlU4ZVdUVFTS4roff/wxefPmzbRQAnWg95496napUtCokfli6tsXNm1SF+cNDgZPT/jtN/OdPzsSJ/cU+iO51TfJrzCGjGHSgcTCafDgwbRs2TKpoDGXhIQE7t27x7179wB1rbhWrVoxePBgXn/99XSLJYC//oLISHW7UycwYbLxNAwG8PaGIkXU17//rv3SKRcuXNA2AGExklt9k/wKY0jBZGHWnOk7vcIpT548JhdPKY93c3PLtFBKdOhQ8nabNiaFkK4yZWD69OTXAwcmF2ha6NKli3YXFxYludU3ya8whhRMFjZw4EACAgLw9/e32jVTFk79+/fHzsQpsVOOf2rTpk2mhVIiX9/k7datTQohQ337Qrt26vatWzB2rGWuY4zFmS2OJ3Isya2+SX6FMWQeJiux5uK75hYREUGhfyY7Cg8Px9HRMdNj4uOhaFF49kwdv3T/vnkGfKfn8mWoVw9iYsDeHk6fhtdes8y1hBBC5BwyD5OweWfOqMUSqL1LliqWAGrWhOHD1e24OPWpOS3+DJDlFfRLcqtvkl9hDCmYhEWkHL9kqdtxKY0cCVWrqtu+vrBsmeWv+aLhiVWb0B3Jrb5JfoUxpGASRilevDjFszDR0cmTydumzu5tjAIF1KfmEg0dCo8eWf66KU2bNs26FxRWI7nVN8mvMIYUTCJTjo6OBAcHExwcbNT4JYBz59TPdnZgrQl0O3eGXr3U7ZAQGDHCOtdN1LNnT+teUFiN5FbfJL/CGFIwCbOLi4OLF9Xt6tXBmmsNT58OiZP2LloER45Y79p+fn7Wu5iwKsmtvkl+hTGkYBJmd+2aOss3QN261r122bKQcvzmgAHJS7NYWoUKFaxzIWF1klt9k/wKY0jBJDIVFRVF27Ztadu2LVFRUZm2T7wdB1CnjgUDy8BXXyUvw3LuHMycaf0YhBBC6IsUTBZmzZm+LSUhIYGDBw9y8OBBEhISMm2fsmCydg8TqHMxzZ+fPJXBmDFw547lr3v79m3LX0RoQnKrb5JfYQwpmCxMi5m+tXb+fPK2FgUTgIsLfPGFuh0RAd99Z/lrurq6Wv4iQhOSW32T/ApjSMEkzC6xh6lAAahSRbs4xo+HxJkQVq2CPXsse73Vq1db9gJCM5JbfZP8CmPI0ihWkluWRomIUJ9SUxRo0gS07lhbvBj691e3a9SAs2fBwcEy14qNjSVv3ryWObnQlORW3yS/+iVLowibdfFi8rIkWt2OS6lfP0jsbb9yBSw5P52Xl5flTi40JbnVN8mvMIb0MFlJbulhWrIEPvlE3Z46FYYMsUaEL3f6NDRuDAkJ6m3CixehUiWtoxJCCGFp0sOUiblz51KlShXy589P48aN8fX1fWn7tWvX4uzsjIODA87Ozqxfvz7V18eMGUOtWrVwdHSkaNGidOjQgePHj1vyLdicggULUrBgwUzbaf2EXHoaNICBA9XtqCgYPNgy15EFPPVLcqtvkl9hDN0VTD4+Pnh6ejJq1ChOnTpFq1at6Nq1K4GBgem29/Pzw93dHQ8PD86cOYOHhwe9evVKVRDVqFGDOXPmcO7cOQ4fPkzlypXp1KkTwcHB1npbmnJ0dCQiIoKIiIhMl0axhSfk0jN2LJQqpW6vXw/bt5v/Gp8kdq0J3ZHc6pvkVxhDd7fkmjVrRqNGjZg3b17Svtq1a/POO+8wYcKENO3d3d0JCwtje4rfoF26dKFo0aKsWLEi3WskdvHt2bOH9u3bGxVXTr4llxVlykBQkPp02sOHyXMh2YLffgMPD3W7alW1uDPnsi2LFy+WH7w6JbnVN8mvfsktuQzExMRw4sQJOnXqlGp/p06dOHr0aLrH+Pn5pWnfuXPnDNvHxMSwYMECnJycqF+/vnkC14mQELVYAnWGb1sqlgD+/W9o3Vrd/vtvmDTJvOd/7bXXzHtCYTMkt/om+RXG0FXBFBISQnx8PKUS7738o1SpUgQl/iZ/QVBQkFHtt2zZQqFChcifPz/Tp09n9+7dFE+c5Ccd0dHRhIWFpfrIqZ4/f0737t3p3r07z58/z7DdpUvJ27b488dgAG9vsLNTX0+YANevm+/8oaGh5juZsCmSW32T/Apj2I0ZM2aM1kGYy7Nnz5g2bRqffvppqsUUDx06xNmzZxk0aFCaY8aOHUuPHj2om2LAzdmzZ9m0aROjR49O2leuXDk++ugjevfuTVhYGJMmTeLDDz/McEzPuHHj6NKlCxMnTmTixIlMnz4dUG8BLl++nCZNmjBx4kRat27NuHHjqFatGqtWrSI+Pp7z589z9uxZFEVhxYoV1K5dm8mTJye1rVOnDosXLyZ//vwcP36ca9euER4ezqZNm6hUqRLTp09PatukSRNmzZpFsWLF2LdvH0FBQQQFBbF3716KFSvG7Nmzk9q2aNECLy8vKlasyObNm3n27BnXr1/n8OHDTJkyhatXrzJ48GCmTZuWdEyNGjVYsWIFCQkJrF37CD8/tfgsX/4QXboUYcqUKUlt69Wrx6JFiyhQoAB+fn5cv36dsLAwtmzZQvny5ZkxY0ZSWxcXF2bNmkWJEiXYs2cPDx8+5N69exw4cICiRYsyZ84cWrZsyfjx42nZsiU//fQTlStXZuPGjYSHh3Pt2jX+/PNPHBwcWLp0KQ0bNsTLy4v332/Nrl1+3L5dgbg4OHnyKaVK7SMg4AKxsbGsXr2aGjVqpIq7fv36LFy4EEdHR44cOcKtW7d48uQJ27Zto1y5cklxjx49mjfeeIPp06dTsmRJdu/eTUhICHfu3OHgwYMUKVKEefPm4ebmxk8//ZR0/ipVqrB+/XoiIyO5cuUKf/31F/b29ixbtowGDRrg5eWV1LZ27dr89ttv5MmTh5MnT3Lp0iWio6NZs2YN1apVY+rUqUltGzRowM8//0zhwoXx9fUlMDCQJ0+esGPHDkqXLs3MmTOT2jZr1oxp06ZRunRpdu7cyaNHj7h9+za+vr688sorzJ8/n+bNm6eK+9VXX2Xt2rVER0dz8eJFTp48iZ0Ce0MmAAAgAElEQVSdHb/++it169Zl0qRJSW2dnZ1Zvnw59vb2/PXXX1y+fJmoqCjWr1/Pq6++murfVMOGDZk3bx6vvPIKhw4d4s6dO4SEhLBr1y5KlSqVKm5XV1emTJlC2bJl2bZtG0+fPuXWrVscPnyYQoUKsWDBApo1a8aECRNS/V9bvXo1MTExBAQEcPr0aQB+//13XnvttVT/11577TWWLVvGnTt3ePToEVevXiUyMpINGzZQuXLlVHE3atQIb29vihYtyv79+7l37x7BwcHs3r2bEiVKMGvWrKS2bm5uTJo0ifLly7NlyxZCQ0O5ceMGfn5+FCxYkEWLFuWonxHHjh3DwcGBxYsX07hx41Rxp/wZce7cOc6dO0d8fDw+Pj7UqlXLJn5GeHt706VLl1T/12rWrMkff/wBwOnTp7lwwfSfEePGjaNp06byM8KKPyMaNWrEpEmTGDlyJA6mTsKn6Eh0dLRiZ2enrFu3LtX+b775RmndunW6x1SoUEGZNm1aqn3Tpk1TKlas+NJrVatWTfnpp58y/Prz58+V0NDQpI/bt28rgBIaGmrku7Ed4eHhCqAASnh4eIbt/vtfRVFnYVKUTZusF19WhYUpStmyybFu2GCe8164cME8JxI2R3Krb5Jf/QoNDTXb715d3ZLLly8fjRs3Zvfu3an27969Gzc3t3SPcXV1TdN+165dGbZPpCgK0dHRGX7dwcGBV155JdWH3qW8vfXqq9rFkZnChVNPYPnttxAZafp5d+3aZfpJhE2S3Oqb5FcYxfT6zbasXLlSyZs3r/LLL78oAQEBiqenp+Lo6KjcvHlTURRF8fDwUEaMGJHU/siRI4qdnZ0yceJE5eLFi8rEiRMVe3t75dixY4qiqL0rI0eOVPz8/JSbN28qJ06cUPr37684ODgo58+fNzouc1a51mZsD1OLFsm9NhERVgwwGxISFKV9++R4R40y/ZxhYWGmn0TYJMmtvkl+9Ut6mF7C3d2dGTNm8OOPP9KgQQMOHTrEtm3bqPTP1M6BgYHcv38/qb2bmxsrV65kyZIl1KtXj6VLl+Lj40OzZs0AsLOz49KlS7z33nvUqFGDN954g+DgYHx9feXJihf8/bf6uXRpMGKOS00ZDDBnDiQuHzV5srp0iilmzpxpemDCJklu9U3yK4yhu3mYbFVOnofJmKVRIiMhcbebGxw5Ys0Is2/kSJg4Ud3u1Al27LC96RCEEEJkj8zDJGzOjRvJ27Y8fulFo0dD4gOVu3bBmjXZP5csr6Bfklt9k/wKY0jBJDLl6OiIoigoipLhNAo5ZcD3ixwdYcaM5NfffANPnmTvXAMGDDBPUMLmSG71TfIrjCEFkzCLW7eStytX1iyMbOnRA954Q90OCoLhw7N3nlWrVpkvKGFTJLf6JvkVxpCCSZjFnTvJ2ynmDM0RDAaYOxf+GabFwoVw8GDWz+Pq6mrewITNkNzqm+RXGEMKJgvz9vbG2dkZFxcXrUPJtufPn9OzZ0969uyZ4dIot28nb+e0ggnUmH/6Kfn155/DS1aBSdftlN8EoSuSW32T/ApjSMFkYQMHDiQgIAB/f3+tQ8m2+Ph41qxZw5o1a4iPj0+3TcoepvLlrRSYmX31FTRvrm5fuQLjx2fteIM8Xqdbklt9k/wKY0jBJMwi8Q+0okWTpxfIaezs1NtxiXMzTZwI588bf3z5nFopikxJbvVN8iuMIQWTMFlCAty9q27n9J87deokD/qOi4PPPoMMOtXSOHbsmOUCE5qS3Oqb5FcYw6wTVwYHB3PmzBlu3rxJaGgo0dHRFCpUiHLlylGzZk3q1KljrkvlOHqeuPLBA3V2b4Bu3WDrVmtHaF7Pn0ODBnD5svp69mwYNCjz40JCQihevLhlgxOakNzqm+RXv8z5u9fe1GCuX7/OL7/8wvr167l8+TKKomBnZ4eTkxMODg48fvw4aZFaJycnunXrRp8+fejcubOplxY2Qg/jl1LKnx8WLIA2bdTXI0fC229nPph9/vz5jB492vIBCquT3Oqb5FcYI9u35IKDg+nbty9NmjTh8uXLfPnll5w4cYLHjx8TGxtLSEgId+/eJSoqivDwcP7++2+8vb0pXLgwX331FfXr1+dgdp7dFjYnJ08pkJHWrdXbcQDh4TBggLpM78vID1z9ktzqm+RXGCNbBdPOnTt5/fXXqVGjBtevX2fNmjV88803NGjQgCJFiqRpX7BgQapUqcKHH37I/PnzuXbtGqNGjWLIkCEMGzYMWc4uZ0v5RK4eepgSTZqUfKtx61ZYtuzl7WV5Bf2S3Oqb5FcYI8tjmLZu3cq8efP49ddfKVq0qMkBTJw4kStXrrB48WKTz2XLcvIYJkVRiIyMBNTi98VHcEeMAC8vdXvPHmjf3toRWs6mTertOIBXXlGfmsuoFy08PDxprJfQF8mtvkl+9UvTxXcvXrzIpk2bzFIsAYwYMYIPPviA81l5fltYlcFgwNHREUdHx3TnK9FrDxPAW29B377qdlgY9O+f8a25hQsXWi8wYVWSW32T/ApjmPUpOZGWt7c33t7exMfHc+XKlRzZw5SZNm3g0CF1Ozw8587DlJGnT9XpBhKnTpg3D778Mm27S5cuUatWLesGJ6xCcqtvkl/90rSH6WXOnj3Lt99+y969e8152hxNDzN9R0dH069fP/r165f0xGNK9++rn195RX/FEkCRIvDLL8mv/+//4Pr1tO3OnTtnvaCEVUlu9U3yK4xh1oJpzJgxzJ8/H3d3d3OeVmgsLi6OZcuWsWzZMuLi4tJ8/cED9XPiAGk96txZXV8OICICPv5YnbAzpfQeeBD6ILnVN8mvMIZZC6bSpUuzc+dOfv75Z3OeVtiw58/VsT0ApUppG4ulTZkClSur24cOwaxZqb9euHBhq8ckrENyq2+SX2EMsxZMDg4O1KxZk/fee8+cpxU2LLF3CfRfMBUuDCkf5hw5Ei5dSn598eJF6wclrEJyq2+SX2EMsxZM//3vf/n22285c+aMOU8rbFhQUPK23gsmgNdfh6+/VrefP4fevSFxWFfHjh21C0xYlORW3yS/whhmLZiOHDnCwYMHcXFxoV27dvz4448cOnSImJgYc15G2JDc1MOUaOJEqF1b3T59Wu1pAli6dKlmMQnLktzqm+RXGMOs0wq0atWKLl26cPfuXY4cOZI0t5KDgwPNmzfn9ddfp2fPnrny8c2cPHHlyxbfXbgweTD0zz8nb+vdmTPQtCkk/i2wbRt07aptTEIIIVKz2WkFqlSpwsiRI5k7dy5nzpzh4cOHrFy5kn79+hEUFMR///tf2rVrZ85LCo3lxh4mgPr1YfLk5Nd9+8Lw4dO1C0hYlCydoW+SX2EMs/YwnT17lpkzZ1KuXDnee+896tevn+rrDx484MGDB9SrV89cl8wxcnIPk6IohISEAFC8ePFUs31//TXMmaNuHzsGzZppEaE2FAXefFNdZw6gY8cEduzIQx6z/hkibEFcXBz29vZahyEsRPKrXzbbw1SvXj1++eUX+vXrl7T2WEqlSpXKdcWSt7c3zs7OuLi4aB1KthkMBkqUKEGJEiXSLI2SW3uYAAwGWLIkef6p3bvzMF06mXRpypQpWocgLEjyK4whS6NYSU7uYXqZ1q3B11fdjoyEAgW0jUcLu3dDp07qdt684OcHjRtrG5Mwr5s3b1I5cRIuoTuSX/3StIfpwYMHREREmHTR9Fy+fNns5xTmER0dzcCBAxk4cGCapVESe5gKF86dxRJAx44wdKi6HRurTjWQOJmn0IfDhw9rHYKwIMmvMEaWC6Y8efLwySefcD9xATEz+PXXX5n14rTJwmbExcUxd+5c5s6dm2ZplMSCKbfdjnvRuHFQq9YzAK5ehU8+Ucc4CX2oVKmS1iEIC5L8CmNkuWAqUaIEXl5e9OrVi8WLF5Pw4oJaWXDr1i369+/Pvn37pGDKgZ4/h9BQdTu3F0z58sEPP1zAyUl9vXYtzJihbUzCfOLj47UOQViQ5FcYI1uDvitXrsz27ds5duwYtWrVYsKECZw5cwZjhkOFhYWxbds2PDw8cHFxoVmzZixZsgQ7O7vshCI09PBh8nZuL5gADIYbLF+e/HrYMDhyRLt4hPmYs0dd2B7JbzoOHVIfAy5bVn3CZcMGrSPSXLafoyxUqBALFizg1KlTTJ48mbFjx2Jvb4+LiwvlypWjSJEiODk5ERMTw+PHj3n8+DE3btzg3LlzlChRgv79+xMQEEDx4sXN+X6EFaVcFiXxSbHcrGnTplStCiNGqLOBx8VBr15w8qQUlDld06ZNtQ5BWJDkNx0REeqEcx9/DLI+LGCGaQUaNmzIH3/8wYMHD1i6dCmNGjXi8ePH+Pr6snLlSjZu3Mi5c+fImzcv77//Pvv37+fu3buMHTtWiqUcLjdPKZCetWvXAjB2rLrmHMC9e+og8BeGfokcJjG3Qp8kv+no2lUdnPnuu1pHYjOyPa1A+/bt2bt3r7nj0a2cPK1ARkujLFoEn32mtpk/H774QqsIbUN0dDQODg6AWkw2bAiJPf3/93+pZwYXOUvK3Ar9kfxmwmCA9evhnXe0jiTLbGLiypMnT8pAuVwuODh5u2RJ7eKwFZNTVESlSsGqVZA4efCUKfDbbxoFJkw2WapdXZP8CmNku2AKDQ1lwIAB3L5925zxCBtUoEABbty4wY0bNyiQYrKlR4+S2xQrpkFgNmb06NGpXrdsCTNnJr/+9FP46y8rByXM4sXcCn2R/ApjmDSGKSEhAU9PT1xdXfnqq69Yu3YtT548MVdsuqCHpVHy5MlD5cqVqVy5MnlSLJQmBVNq6S3gOWBA8m3L6Gi1R1seyMl5ZHFWfZP8CmNkewxT5cqVuXnzJqAuznry5En27t3L/v37iYiIwM3NjQ4dOtC2bVtZ1JCcPYYpI2+9BZs3q9v378uTcvfu3aNs2bJp9sfEQLt2yVMMuLrC/v0gQyZyjoxyK/RB8psJGcMEmNDDNH/+/KRtg8FA48aNGTZsGNu3b2fPnj106dKFAwcO4OrqyuTJk42ao0nYppiYGIYOHcrQoUOJiYlJ2v/4cXKbf/1Lg8BszPbt29Pdny+fOpFlhQrqaz8/mQk8p8kot0IfJL/pCA+H06fVD4AbN9TtwEBt49KQxRffjY+Px8vLi8OHD7Nu3Try589vycvZrJzcw5TRU3K1a8OlS1CoEDx7pmWEtsHf3/+lt15PnoRWrdRFigFGj1anIBBWdvMmFC1K0rTsRsgst7pz967aBZpLpn7Jdfk1xoEDyfOjpNS3Lyxdau1oss0mepiMZWdnx/fff89nn33Gxx9/bOnLCStKHMMk45dUj1N2uaWjUSP44w+1dxvUKU5y0M8d/WjRQu3u+/FHePrUqEMyy61uXL8O/ftDpUrwzTdaR2M1uSa/WdG2rdoN/uJHLv6hle2C6dSpU1lq36NHDypXrsymTZuye0lhQxQl+Zac3I5TRUREZNrm7bdh+vTk1599Bvv2WTAokdazZ+rH//4HFSsaVTgZk9scLbFQqlEDli+H+Hj1lkwuofv8CrPIdsH0ySefZPmYIUOG8PPPP2f3kkabO3cuVapUIX/+/DRu3BhfX9+Xtl+7di3Ozs44ODjg7OzM+vXrk74WGxvL8OHDqVu3Lo6OjpQtW5Y+ffpw7949S78NmxYWpv5MBelhSlSrVi2j2n37LXz9tbodF6dOpBsQYMHARPoSEowunIzNbY6TXqGUC6el121+hVllu2B6lPKZciOVKFEiW8dlhY+PD56enowaNYpTp07RqlUrunbtSmAGA9X8/Pxwd3fHw8ODM2fO4OHhQa9evTh+/DgAkZGRnDx5kh9++IGTJ0+ybt06rly5wltvvWXR92HrZEqBtPbs2WN02+nT4Y031O3QUOjSBWRKM40YUThlJbc5ghRKqeguv8Iisj3oO0+ePLRv357u3bvToUMH6tSpY9RxLi4u+Pv7Z+eSRmnWrBmNGjVi3rx5Sftq167NO++8w4QJE9K0d3d3JywsLNVTEl26dKFo0aKsWLEi3Wv4+/vTtGlTbt26RcWKFY2KS2+Dvv/6CxLHSA4YAHPnahigjQgLC8tSbsPDoU0bOHlSoSCR1KwBu3blmnG22ild+uW3mwwGcHRUuwIHDIAiRbKcW5t14wZMmgS//66+z5et1tC1K6xebb3YNKSb/Io0wsLCcCpb1iy/e02aIGnv3r3s+2cARsmSJWnXrh3t27enQ4cO6RYS58+ft+h6PTExMZw4cYIRI0ak2t+pUyeOHj2a7jF+fn4MHjw41b7OnTszY8aMDK8TGhqKwWCgSJEiGbaJjo4mOjo66XVYWJgxbyHHkB6mtGbNmpWlGYMLFYJt26CjWyRnrxeCK0Bli4UnjKUoakE1frz6AeTKX6Xbt6v/SHOBXJlfkWXZviVXqlQpDh48yP/+9z9at27N06dPWbFiBZ999hlVqlShevXqfPnllyxfvpw9e/Ywb948unTpwqBBg8wZfyohISHEx8dTqlSpNLEGBQWle0xQUFCW2j9//pwRI0bw4YcfvrRanTBhAk5OTkkfFf6ZhOf69etMmjSJ6OjopNllx40bx7179/jll1/w9/dn586drFu3joCAAGbNmkVYWFiqto8fP8bb25uzZ8+yadMmtm7dysmTJ/n5558JDg5O1TYyMpKpU6dy+fJlVq1axd69ezl69ChLly7l9u3bqdrGx8czYcIEbt26xa+//srhw4c5cOAAGzduZOvWrQwZMgSDwcC4ceNSFUxXrvjh7+/Pjh07WL9+PefPn2f27NmEhoamOv/Tp0+ZM2cO586dY8OGDWzbto0TJ06wYMECHjx4kKptVFQUU6dO5erVq/j4+LBv3z6OHDnC8uXLCQwMZPz48SQkJDBu3DgSEhIYP348gYGBLF++nCNHjrBv3z58fHy4evUqU6dOJSoqKtX5Hzx4wIIFCzhx4gTbtm1jw4YNnDt3jjlz5vD06dNUbUNDQ5k9ezbnz59n/fr17NixA39/fxYuXEhQUFCqWYKfP3/O5MmTuXr1KitXruTAgQMcPnyYX3/9lVu3bjFhwgTi4+OTjlm4cBwLFsjU30IIYeuyfUuucePGnDhxIul1dHQ0hw8fZt++fezZsydpcV5D4jPUwHfffcekSZNMjzoD9+7do1y5chw9ehRXV9ek/ePHj+fXX3/l0qVLaY7Jly8fy5Yto3fv3kn7fv/9d/r378/z589TtY2NjaVnz54EBgZy4MCBlxZM6fUwVahQIUfekkvPnDnJA5eXLYM+fbSNxxaMGzcue2tSKQrXVvyJ40fv8rayjgvU4Z231aEldnbmjzPXy+yWHKjdpkOHquN8ChTAy8uL4cOHWyc+S3rwAKZNgwUL1NtxCQkZt81Ft+R0k1+Rhk3cklu3bl2q1w4ODrRv35727dszfvx4nj17hq+vL5cvX8be3p727dvj7OxsUrCZKV68OHZ2dml6hx4+fJimFylR6dKljWofGxtLr169uHHjBvv27cv0G+/g4GDR249ak1tyaX311VfZO9BgoNqHzfCrcpcLHdSJLf/YCAW/U3+vpfibQ5hDRt9Qg0H9xzxqFHzxBaRYaPozT091XFNO9+qr6l87o0eDl5c6+DA+Pv2xTPb2+njPRtBNfkVaLxunl0XZviVXqVKll369cOHCdOvWjcGDB/P1119bvFgCtbeocePG7N69O9X+3bt34+bmlu4xrq6uadrv2rUrVfvEYunq1avs2bOHYrmsQoiJiWHMmDGMGTMmaWkUKZjSyughAWO5uqpLqCQuvbhokdqLJ0uoWJjBoI60nzZNXfbB0zNVsQSm59bmlC6tPqp565b6jyxfvlzdnam7/ArLUHRm5cqVSt68eZVffvlFCQgIUDw9PRVHR0fl5s2biqIoioeHhzJixIik9keOHFHs7OyUiRMnKhcvXlQmTpyo2NvbK8eOHVMURVFiY2OVt956Sylfvrxy+vRp5f79+0kf0dHRRscVGhqqAEpoaKh537AVhIeHK4ACKOHh4YqiKMq//5089evlyxoHaCPOnDmT/YMvXFCUunUV5cIFZcUKRTEYkr+/33yjKAkJ5osz1ytcWP3GGgyKUry4okyfriiRkS89xKTc5gT37yuKp6ei5MunKHZ26vfnzTe1jspqdJ/fXMycv3stvjSKtbm7uzNjxgx+/PFHGjRowKFDh9i2bVtSj1hgYCD37ycPsnVzc2PlypUsWbKEevXqsXTpUnx8fGjWrBkAd+7cYdOmTdy5c4cGDRpQpkyZpI+MnrzLDaSHKa2bN29m/+DChdWlCAoX5oMP1NUHEu8czZoF330nPU1mU6RIpj1KLzIptzlBej1OL3kKWG90n19hFhZffFeo9DYPU9Om4O+v/lKPjc3VvflJtm7dSvfu3c12viVL1DHHif9D/+//1Cl0ZEyTiR4+VAvUTIqklMydW5v35Il6b7hwYa0jsYpcl99cxJy/e02ah0nkXonryBUpIsVSojJlymT/4OfP4c4dKF8e8ucH4OOP1YeYPv1UbTJlivpZiiYTlSyZ5UNMym1OVLSo1hFYVa7Lr8gW3d2SE9aReEtObsclM2kG+4AAqF49zaJy/ftDyuUXp0yBQYNe/jS4MD9Lrk4gtCf5FcaQgklkWVxc8lJbUjAle/fdd7N/cPXqsH+/+vkFn3+uFk2JvUpz58Inn+Tqpb+szqTcCpsn+RXGkILJwry9vXF2dsYlceE1HUi5Lum//qVdHLbm55RdQVmVYtB3ej7/PPVElsuWwYcfwj+zPAgLMym3wuZJfoUxZNC3leTkQd/x8fGcPHkSgEaNGnH9uh01aqhf+/BDdR1PYaIHD9SKqE8fyGCSVYB16+CDD9SB9gCdO8OaNblmyS8hhMgSc/7uNbqHadKkSVSqVAlHR0e6du3K5cuXk762aNEiBg4cSL9+/UwKRtgmOzs7XFxccHFxwc7OLlUPUy4bG/pSKdeUy7KgIJgwQf38Eu++C5s2JY0LZ+dOtWPqwYPsX1pkzqTcCpsn+RXGMKqHaerUqUyYMIHu3btTvHhxTpw4gb+/P6tWrUp6FHP79u288cYbxJtxGnI9yck9TC/aswc6dlS3v/8+aUH3XC8yMpKCBQta5Vq+vvDWW8m3R6tUUYundIZACTOwZm6F9Ul+9cvqPUzbt2/n3LlzLFu2jKlTp3LgwAEOHDjAmDFj2Lt3r3qiPDIcSq9iYmKYPHkykydPJiYmJlUPUy6a2y5T8+bNs9q1WrWCw4ehQgX19Y0b4OYGx49bLYRcxZq5FdYn+RXGMKrKadKkSZp5KlxcXDh48CDLly9n3759FglO2IbY2FiGDRvGsGHDiI2NlYIpA2+88Ub2D758WV1MLsWt7sy89hr4+UHduurrkBB4/XXYsiX7YYj0mZRbYfMkv8IYRhVMdnZ2PH/+nEuXLrFp06ak/QULFmTZsmUcOXKEHTt2WCxIYVukYErfmTNnsn9w/vxqBZQ4OMlI5cqpt+fatlVfR0XB22+Dt7cspWJOJuVW2DzJrzCGUQXT119/zaBBg+jZsyeff/55mq//8MMPNGjQgLx585o9QGF7QkOTt52ctIvD1hQzZVKqSpVg0SL1cxY5OcGOHeDurr5OSFAnt/zqq+Sn6YRpTMqtsHmSX2EMo5ZGKV26NIsWLeLatWs4Ojqm26Zv3764ubmZNThhm6SHKX0FsrA2WRqxseo9teLFIRt/eDg4wB9/QOXK4OWl7ps/Hy5dUqcdkN8HpjEpt8LmSX6FMbI0UrtatWovXXOnujyikytIwZS+K1euZP/gc+egbFn1czblyQMTJ6rTOeXLp+47cACaNoULF7IfmjAxt8LmSX6FMazyaNvixYsZO3YsT548scblbIreZ/qWgilZ+/bts39w1arqaO2qVU2Ow8NDLZQS57+8fl0dT75hg8mnzrVMyq2weZJfYQyLFEzBwcEkpFgd9JNPPqFHjx5888037N+/3xKXtFkDBw4kICBAV4s7yhim9C1btiz7Bzs5QffuZvuGurqCvz80bKi+fvYMevSAESNkDbrsMCm3wuZJfoUxzLo0yv79++nZsydPnjyhQIECtGrVip49e/Lee+/h5OSEoih88MEH+Pj4mOuSOUZOnrgyPj4eX19fAFq1akXDhnacO6c+0BUVpXFwehEcDKtWQa9eUKKE2U4bEQH9+0PK/3Lt2sGKFVCypNkuI4QQNkmTpVGM8Z///IeBAwcye/Zs+vXrx/Xr1/n0008pU6YMvXv3Zv78+YSEhJjzksIK7OzsaNu2LW3btk21NIrcjkvNpOUV7tyBIUPUz2bk6KgWR9Ong/0/j3js2weNGqnTEQjjyNIZ+ib5FcYwaw/TV199xdy5c1PtO3fuHH/88Qfr168nJCSE2bNn07t3b3NdMsfIyT1ML3JygrAwqFULLl7UOhrbER8fj52dndZhZOjwYbUD6/599XWePPCf/8CoUcnFlEifredWmEbyq18228OUnrp16zJhwgQuXbpESEhIriyWcrrY2Fi8vb3x9vbm+fNYwsLU/TJ+KbVJkyZpHcJLtWwJJ08mT3KZkABjxqi36G7f1jIy22fruRWmkfwKY5i1YGrVqhUb5FEc3YmJiWHQoEEMGjSIkJCYpP1ySy61Dz/8MPsHX70KHTqony2odGl18eSxYyHxD2pfX6hfH9avt+ilczSTcitsnuRXGMOsBVPv3r05f/48ixYtMudphQ1J7F0C6WF60aFDh7J/sL29OtjbCvfG7Oxg9Gg4eBAqVlT3PXkC774LAwZAZKTFQ8hxTMqtsHmSX2EMsxZMf/75J0uWLOGLL76gWrVqfPrpp/z222/cMfNAVqEdKZgyVqVKFVMOVkdnm3KOLGrRAk6fhvffT943fz40aABHj1otjBzBpNwKmyf5FcYwa8H03Xff0bFjR7777juqVavGque3hvoAACAASURBVFWr6NOnD5UqVaJ69ep8/vnnUsnncM+eJW/n8LHrZhdnygRH8fFqNRofb76AjFC0qDqbwYIFkLg6xNWr6ninoUPh+XOrhmOzTMqtsHmSX2EMs/b/16lTh3nz5iW9jo+Px9/fn/3793PgwAFWrFjB5s2buZ/4mI7IcaRgylhQUFD2Dz5zBho3hhMn1Gf+rchggM8+g9atoV8/OHYMFAWmTFEnH1+2TF1eJTczKbfC5kl+hTHM2sOUL3EBq3/Y2dnRvHlzRo4cyc6dO3ny5AnnTFgrKyfS29IoKW/JScGUWuPGjbN/cJUqalePhrcGatZUpx7w8kpei+7SJXXW8OHDc/fYJpNyK2ye5FcYw6wFU82aNTl48GCGX7e3t6d48eLmvKTN09vSKNLDlDGTnhAtWhR69lQ/a8jODoYNU6cfSPwdkpAAkyZB3bqwe7em4WlGnv7VN8mvMIZZC6YBAwawbt06dufWn6o65eDgwJYtW9iyZQuRkQ5J+6VgSu3rr7/O/sGPHsHSpepnG/Daa+DnB+PGJfc2Xb8OnTpBnz7qSi65iUm5FTZP8iuMYdaCaceOHaxcuZIuXbrQrVs3pk+fzqlTp8x5CaEBe3t7unfvTvfu3YmISB72JgVTalOmTMn+wbduwccfq59tRN686izgZ85Aq1bJ+3/9Vb195+2dexbyNSm3wuZJfoUxzLo0ipubGy1atODu3bscPnyYO3fuYDAYKFKkCK1bt+b111+nR48eVKhQwVyXzDH0sjSKpyfMnKlu+/lB8+baxqMbiqI+IWdnp47CtjEJCbB4sfrkXOJaggD16sHs2eqAcSGEsDU2uzRK9erVmTx5Mn/88QeBgYFcvXqVefPm0blzZ44fP46npydubm7mvKSwgtjYWJYuXcrSpUt5+jQ2aX/hwhoGZYNMWsDTYFAnrbTBYgnUdec+/VRdO7BPn+T9Z89CmzbQu7fZ1w22KbI4q75JfoUxzNrDdOTIERYsWECjRo3o0aMHFROnEf7HpUuXePToES1atDDXJXOMnNzDFBERQaFChQB4++1wNm50BCAwEHJhZ2GGgoKCKF26dPYO/vtvGDwYpk+HqlXNG5gFHD0KgwZByjvujo7qLbwhQ8DBIeNjcyKTcitsnuRXv2y2h6lFixYsXboUV1dXAgMD03y9Vq1aubJY0hN5Si5jmzdv1joEq3FzA39/+PlnKFZM3RcRAd9/D7Vqwe+/q7fx9CI35TY3kvwKY5i1YAIwGAw0bdqUli1bmvvUwgakLJj+6XQS/2jQoEH2D65aFTZtyhG9S4ns7ODzz+HKFRg4UL1tB3DzJnz0ETRpoi70qwcm5VbYPMmvMIbZCyahb4kFU6FCyavdC9UjU6YEUBT1kTPz3SG3mn/9C+bMUedu6tQpef+pU9CxI3TurK5Zl5OZlFth8yS/whhSMFmYXmf6lttxaUVFRWX/4FOn1Of4c/A0HPXrw86dsGuXuoBvol271NVePvhAHTSeE5mUW2HzJL/CGFIwWZheZ/qWgimt6tWrZ//gSpVgyRL1cw7XsaO6JN5vvyW/HUUBHx91Qsx//xsuX9Y2xqwyKbfC5kl+hTGkYBJZIgVTxvbv35/9g4sVU1e+TRxBncPlyZNcGE2bBiVKqPsVBf74A5yd1ekJrl7VNk5jmZRbYfMkv8IYUjCJTDk4OLBq1SqWLVsFqM+LS8GUVp+UExRl1ZMnsHq1+llHHBzU2RJu3FAX9U1cSjIhQZ0xvFYt+PBDdTZxW2ZSboXNk/wKY0jBJDJlb29Pz5496dChJ6AujSKTVqY1e/bs7B984wb06qV+1iFHR3VR3xs3YMIEdaA4qIXTihXqmKdu3eDQIdsc925SboXNk/wKY5h14kqRsZw8cWWiK1fUNcQAPDxg+XJt49GV+Hh1IiNHx1zx+OGzZ+qSKtOnQ0hI6q+5usKIEfDGG8lTFQghRHbY7MSVtmLu3LlUqVKF/Pnz07hxY3x9fV/afu3atTg7O+Pg4ICzszPr169P9fV169bRuXNnihcvjsFg4HROf0Y6i+Li4li9ejXr168G1NVWZQ6mtExaXsHOTr3PmQuKJVB7KL//Xl1reM6c1GPd/fzg7behdm31aynn/tKKLJ2hb5JfYQzdFUw+Pj54enoyatQoTp06RatWrejatWu6M48D+Pn54e7ujoeHB2fOnMHDw4NevXpx/PjxpDYRERG0aNGCiRMnWutt2JTo6Gh69erFiBG9gGhACqb0DBo0KPsH37ihLsim01tyGSlYUJ308upV9am6OnWSv3blCnz9NZQvry63cv26dnGalFth8yS/whi6K5imTZtG//79+fTTT6lduzYzZsygQoUKzJs3L932M2bMoGPHjowcOZJatWoxcuRI2rdvz4wZM5LaeHh48J///IcOHTpY623YPEdHrSOwPb/99lv2D46Lg+Bg9XMulDev+lTd2bOweTO0a5f8tbAw9dZdtWrwzjuwb5/1xzmZlFth8yS/whi6KphiYmI4ceIEnVJONwx06tSJo0ePpnuMn59fmvadO3fOsL2xoqOjCQsLS/WhJ9LDlFabNm2yf3D16uo6Irl8PhiDQR27tHev+uRc//6QP7/6NUWBjRuhfXv16bqpU9OOf7IUk3IrbJ7kVxjDbsyYMWO0DsJcHjx4wJQpU/jmm2+oUKFC0v5Tp05x5MgRvv766zTHDB8+nI8++oi6desm7btx4wa///47P/zwQ6q2T58+ZebMmXz55ZeZrmw9btw4unTpwsSJE5k4cSLTp08HwN3dneXLl9OkSRMmTpxI69atGTduHNWqVWPVqlXEx8dz/vx5zp49i6IorFixgtq1azN58uSktnXq1GHx4sXkz5+f48ePc+3aNcLDw9m0aROVKlVi+vTpSW2bNGnCrFmzKFasGPv27SMoKIigoCD27t1LsWLFmD17dlLbFi1a4OXlRcWKFdm8eTPPnj3j+vXrHD16lI0bN/7zzr4H8vHuu7Bjxzhq1KjBihUrSEhI4Ny5c5w7d474+Hh8fHyoVasWU6ZMSTp/vXr1WLRoEQUKFMDPz4/r168TFhbGli1bKF++PDNmzEhq6+LiwqxZsyhRogR79uzh4cOH3Lt3jwMHDlC0aFHmzJlDy5YtGT9+PC1btuSnn36icuXKbNy4kfDwcK5du8aff/6Jg4MDS5cupWHDhnh5eSWdv2bNmvzxxx8AnD59mgsXLhAbG8vq1aupUaNGqrjr16/PwoULcXR05MiRI9y6dYsnT56wbds2ypUrlxT3kCFD6N69O9OnT6dkyZLs3r2bkJAQ7ty5w8GDBylSpAjz5s3Dzc2Nn376Ken8VapUYf369URGRnLlyhX++usv7O3tWbZsGQ0aNEgVd+3atfntt9/IkycPJ0+e5NKlS0RHR7NmzRqqVavG1KlTk9o2aNCAn3/+mcKFC+Pr60tgYCBPnjxhx44dlC5dmpkzZya1bdasGdOmTaN06dLs3Pn/7d17XJRV/gfwD9cBEcQLgqgg3vOSAiZKeKkt0lrNrF9UW7ltVrZia1a/sstP20yEytIVLTOvq2XeqhW3pETSzM0QExUVlQTUQUEBAR1wOL8/znIZuT3MMzDPDJ/36/W8nplnnjOc4eszfjnnPOd8h/z8fGRnZ2PPnj3w8vLCxx9/jBEjRpjUu2fPntiyZQsMBgPS09Nx8OBBODk5Yd26dRg8eDDi4uKqzh0wYADWrl0LZ2dn/Prrrzhx4gSuXbuGbdu2oWfPnli4cGHVucHBwVi2bBn69PGCt/ePmDQpFx07OuP0aWdcv+4KAMjPlzOIL1oksGVLOjw9y3H48DcoKCjA2bNnsXfvXrRt2xbLly9HWFgYYmJiTK61TZs2oaysDMeOHasak7h+/XoMHDjQ5FobOHAg1qxZg5MnTyIvLw8ZGRkoLS3FV199hR49epjUOyQkBPHx8Wjfvj2SkpJw/vx5XLp0CYmJifDx8cHixYurzg0PD0dcXBy6deuG7du3o7CwEJmZmfj555/Rpk0brFixwqa+I/bv3w+dToeVK1ciNDTUpN628B0RFxeHe+65p9m/I+bNm4fhw4fzO8JC3xFeXl748ccfkZOTg7y8POzcuRO+vr4m9Q4JCUFcXBxmz54NnU6nPKGoi7Aj586dEwDEvn37TI7PmzdP9OvXr84yLi4uYsOGDSbH/vnPfwqdTlfr3MzMTAFApKamNlqX69evi8LCwqotOztbABCFhYVN+ETaUFxcLAD8dysWgBDr1lm7VtqTkJBgfuGDB4VwdZV7qlNZmRBffCHE2LFCyPYm0613byFiY4XQ6y3/s1XFljSP8bVfhYWFFvu/16665Dp16gQnJyfo9XqT4xcvXoSvr2+dZfz8/Jp0vlI6nQ5eXl4mmz1hl1xtqv7NdOsmp8Tu1s1yFbIzLi5AVBSQlCRnEH/lleqJMAHg1Cng1VeBrl2BCRPkPKDXr1vmZ6v9PiBtY3xJCbtKmFxdXREaGorExEST44mJiQgPD6+zzMiRI2udv3PnznrPJ4kJU20pKSnmF/bxkbeLVa4hQg3q2xeIiwNycuQadTUHiRuNwPbtch5QPz/gueeAn35SN1BcVWxJ8xhfUsKuEiYAmDVrFlasWIGVK1ciPT0dL774IrKysjBt2jQAcgr82bNnV53/t7/9DTt37kRsbCyOHz+O2NhYfP/995g5c2bVOZcvX8ahQ4dw7NgxAMCJEydw6NChWi1T9srV1RWrVq3C3XevAiDHkDBhqu3+++83v3BhIZCQIPekmE4nE6MffpBTE7z2mmkjXWEhsHw5EBEhx9O/8YYcTN7U5ElVbEnzGF9SxAJdhJoTHx8vAgMDhaurqwgJCRHJyclVr40ZM0ZMmTLF5PxNmzaJfv36CRcXF9G/f3+xZcsWk9dXrVpVYwxP9TZnzhzFdbJkP6q1TJ1aPV7k8GFr10Z73nnnHfMLp6TIX2xKiuUq1ErduCHE998L8eSTQnh41D3eqV8/Id56S4i0NGXvqSq2pHmMr/2y5P+9XBqlhdjD0iiPPgp88YV8fOYMEBRk3frYlfJyeY98p05ysA5ZREkJsG0bsGaNnL+poqL2OQMGyFaqqCg5XQER2Q8ujUIt6saNG0hISMDp0wng0ij1U7W8gosL0KULkyUL8/AAHn8cSEwEzp8Hli4FxoyR8z1VOnYMmDtXLsUyeDDw5pvAgQOm3XZcOsO+Mb6kBFuYWogttzCVlJSgbVWGVAzAA6WlgLu7NWulPdeuXYO7ub+Us2eBd94B3nrLdGE1ahYXLgCbNwNffgns3Vv3OV27AhMnynXtwsKuwdub/+DtlaprlzSNLUxkVQ4O1bMvU7WlS5eaX/j6deDoUcvdB08N6tJFrlO3Zw+QnS2XXhkxwvScc+eAZcuAceMAPz9H/M//ACtXypYqsi+qrl1qNdjC1Mzi4+MRHx8Po9GIkydP2kULk6enB+xspReLyMjIQJ9WvrSJrbtwQa5l9/XX8s47g6Hu84YMkYnU+PFAeDh7Um0dr137ZckWJiZMLcSeuuS6dPHgX9l12LhxI6KioqxdDbKQ4mLgu+9k8rR1qwElJXUvq+DpCdx1l0yexo0DaqzKRDaC1679YpccWZWHh7VroE0+aiad/O03oEMHuSdNaNsWePBBYO1aYNu2ffjpJznEbNgw0/OuXpV34j37LBAQICfVnDZNjo+6dMk6daemUXXtUqvhbO0KkO3hHXJ1U7Wwo58fMHu23JPmtGnjivBw2f32978DFy/KxX///W/ZCpWfX31uRobcPvlEPh88WM5EfuedwOjRgLe3dT4D1U/1oqzUKrCFiZqMCVPdTp8+bX5hX1+5OBrXtNKkm2PbubOcrmD9eiA3F/jPf+TUBBERgPNNf4ampQGLFsm77Tp2BIYPB15+GfjqK7ZAaYWqa5daDY5haiG2PIapvLwcH320HP/7vwDwLMaNc8G//23tWmlPVlYWAgICzCt89SqQkgKEhspBMaQpTYltcbFcu27XLrmlpDS8FEu/fjLRqtx69TKdJ4qan6prlzSNY5ioRbm4uOCRR6YDmA7AhWOY6rFu3TrzC2dkAHfcIfekOU2Jbdu2wD33ALGxcgLM/HzZmvTCC8CgQbXPP3EC+Owz4Kmn5Hp3XboADz0ELFwoE6/SUgt+EKqTqmuXWg22MLUQW25hAoDjx+VMyAAwZQqwerVVq6NJFRUVcHQ082+Q69eBnBy5ciwnudIcVbG9yaVLwL59csLMvXtlC1R5ef3nOzkBt94qu/KGDwfCwuQSLk5OFqkOwbLxJW1hCxO1KKPRiOTk3QB2AzByDFM95s+fb35hNzegd28mSxqlKrY38fGR45neew/4+WegoADYvRt49105NcHN3+lGI5CaKgeRP/20bKXy9pYNkq++CmzdCvz+e8PdftQwS8aX7BdbmFqILbcw3TwP06uvemDBAqtWSZNU/ZWanS3/B33lFU7ko0Et2QJhNAJHjsiB5L/8IvdHjzaeELVvDwwdCgQHV2/9+tUehE61sYXJfrGFiayKY5jqFhMTY37hq1dlM8PVqxarD1mOqtg2kZOTnEn82WeBFSvkXXaFhfKfR2ysnBuqW7fa5a5cAZKS5NinJ56QLVGenrILb9o04OOPZfJ17VqLfRSb0ZLxJdvFFqZmZo9Loyxc6IEXX7RqlTSJd9rYLy3G9vx5Oaj8l19kl11qKqDXN17O0VH2/g4aVL0NHCgHnLfWJV60GF+yDC6NYoPsqUtu+XIPPPOMVaukSWvXrsWTTz5p7WpQM7CV2Or1wKFD1QlUaipw6pSysi4ucjB5zURq0CCgRw+ZZNkzW4kvNZ0l/+9l7zY1Gbvk6tarVy/zC6elyRG///63nBqaNEVVbFuQn59cz27cuOpjRUVyxZ3UVJlMHToEpKfLGzNrKi+X/wzT0kyPu7vL1qf+/eWYqMp9v372M4mtrcSXrIsJEzWZvXxJWpqhvqXtlejUCZg6Ve5Jc1TF1sq8vIBRo+RWyWgEzpyRg8mPHKneTpwAbtwwLX/tGnD4sNxu1rVr7SSqf385xsqWpj2w5fhSy2HCRE3GhKlul9Ssc9Gli1xbgzRJVWw1yMlJthr16QNMmlR9vKwMOHnSNIk6fhw4fbp2IgUA587Jbdcu0+OurkBQkJy1vFcvOWaq8nFQEKC1pdvsLb7UPJgwUaNcXFxw++1x+OknAHBhwlSPkJAQ8wuXlADHjgEDBrDPU4NUxdaGuLpWj12qqbxctkidOCETqJr7mgsPVyork6+dOFH7NQcH2QJ1cyLVq5ccL9W+fcsvDdNa4kvqMGGiRrm6uqJv31f+mzDx//P6fPPNN3jppZfMK3zihJzGOSUF4Je35qiKrR1wcanucps40fS1vLzq5KgyiTp9WiZYdU1hIIScdiw7W06VcLO2bYHAQJk8BQZWb5XPfX0tn1C19viSMrxLroXY8l1yAPDww8CmTfLx77/LLy4yde3aNbi7u5tbWK4j16ePHGVLmqIqtq2UEMCFCzJ5qtxOnap+fPmyee+r05kmUpVbt25yTFXXrk3/o47xtV+8S45alNFoxLlzB//7LARt29rQaM4W9MEHH+DNN980r7C7u1wwjDRJVWxbKQcHwN9fbjUHnFcqKDBNpk6fBs6elVtWFlDfOGyDQY6zOnmy/p/t7V2dQNW379ChuqWK8SUl2MLUQmy5henmeZiuX/fQ3KBNm3fuHLBkCRAdLb/NiVqxigrg4kXZml2ZRJ09a/q8uFjdz3BzM02gunaV9174+VVvXboA7dq1/Jgqshy2MNmQmjN92wMHBzkwlGqbN2+e+X+lFhTIPs/HH2fCpEGqYktN5uhYnbSMGFH7dSFkl15l8pSdDeTkyL87cnKqHzc0W8D169UtWw3R6WonUTWfVx7z9dXe3X9kWWxhaiH20sLUpk0xSko46rsuubm58PX1tXY1qBkwtrZHCHkHX2USdXMyVbkvLLTcz/T2Bnx8gM6d5b6xjQlW82MLE1kNx0XW7+uvv8azzz5r7WpQM2BsbY+Dg5wHtlMnuZhxfa5eBT78cCMiIqKg16PO7cIFeTdgYwoK5JaRoayOXl6NJ1UdOgAdO8q9t7f9L1OjZUyYqEnc3KxdA+0KDQ01v/DRo8DkycDWrXIlVNIUVbElTfP0BO67rzcaC3F5uRxXdXMiVfN5bi5w6ZLyVquiIrk11i1YycFBzlPVoYNpItXY43btmGhZAhMmahK2MNUvNzfX/MLt2skJbtq1s1yFyGJUxZY0T0l8XVyqB4c3xmCQLVKXLtXeLl6sfezKFWX1rBy71dQpGSoTrY4d5d7bu3pr1870+c1bu3ZAmzYc+A4wYaImYgtT/crKyswv3K0b8N57lqsMWZSq2JLmWTq+Op3y5AqQrVf5+bUTqry86gTp8mV5TuW+oEB5fcxNtCo5OzecUFU+9vKSLXZeXrUfe3jYftLFhIka5ejoAmAOAKBNGxfrVkbDVK14fu2anBq5Z08242kQV7O3b9aOr4tL9R13ShmNMmmqTKJqJlT1Pc7PN2+Q+40bMnlTMo6rPg4O1QlUfUlVY489PeXm7m6d5IsJEzWqosIVwFwAXBalIcnJyRg8eLB5hdPTgdBQLo2iUapiS5pni/F1cpJdbB07Nq2c0SgHulcOUC8srH5881bXa4WFcp6sphKiesyWWg4OcgkdJZuzBbMcTivQQmx5WoH8fHmnCQDcey+QkGDd+mhVQUEBvL29zStcXAykpQGDB4OrG2uPqtiS5jG+ylVUyK+rupKpq1dlQlS5b+ix2olHlSsCwGkFqIWUlFQASAcAuLndAoC3W9RlyZIl5k9u2LYtMHKkZStEFqMqtqR5jK9yjo7VXWXdu5v/PpWJl9IEq/JxSYksd/Wq3Fc+vnHDcp+xPmxhamY1Z/o+efKkTbYw/fZbCYYOla0ejzxSjM8/Z7+cxV24AHz2GfD003LaYCIiUqysrDqBqrnl5hbhsccs08LEpoJmNn36dBw7dgwHDhywdlXMdu1a9WPOTFu/efPmmV/40iUgPl7uSXNUxZY0j/G1fa6uct6pgABgwABg+HDgzjuB++6z3M9glxw16vr16se8gat+M2bMML/wrbfKVibSJFWxJc1jfEkJtjBRo2q2MHEepvqtXbvW2lWgZsLY2jfGl5RgwkSNYguTMnfccYf5hdPT5XQC6emWqxBZjKrYkuYxvqQEEyZqVM0WJiZM9ctQuuJmXTw85F1ynOhKk1TFljSP8SUlmDA1wdKlSxEUFAQ3NzeEhoZiz5491q5Si2CXnDLuarLJgAA56DsgwHIVIotRFVvSPMaXlGDCpNDGjRsxc+ZMvPHGG0hNTcWoUaMwfvx4ZGVlWbtqza683AXAywBeRtu2XBqlPh2bOuVuTQYD8Pvvck+aoyq2pHmMLynBhEmhhQsX4umnn8bUqVNxyy234KOPPkL37t2xbNkya1et2ZWXuwJ4D8B78PR0tXZ1NOvQoUPmFz56FAgKknvSHFWxJc1jfEkJJkwKlJWVISUlBZGRkSbHIyMjsW/fvjrLGAwGFBUVmWy2imOYlJkwYYL5hXv3BhIT5Z40R1VsSfMYX1LCae7cuXOtXQmty83Nxfvvv48XXngB3WvMBZ+amoqffvqpzjk85s2bh3HjxmHBggVYsGABPvzwQwBAVFQU1q5di2HDhmHBggUYPXo05s2bh969e+PLL7+E0WjEkSNHcPjwYQgh8Pnnn+OWW27Be++9V3XuoEGDsHLlSri5ueE///kPTp06heLiYnzzzTcIDAzEhx9+WHXusGHDsHjxYnTs2BG7du2CXq+HXq/HDz/8gI4dO+If//hH1bm33347YmNjERAQgH/961+4evUqzpw5g61bC5Gefh1AAaKi3LBlS3W9+/bti88//xwVFRVIS0tDWloajEYjNm7ciP79++P999+vOvfWW2/FihUr4O7ujp9//hlnzpxBUVERtm/fjm7duuGjjz6qOve2227D4sWL4ePjg++//x4XL17E+fPnsXv3brRv3x5LlixBREQE3n33XURERGD+/Pno0aMHvv76axQXF+PUqVP45ZdfoNPpsHr1agQHByM2Nrbq/fv164cNGzYAkH9dHj16FOXl5di0aRP69u1rUu8hQ4bg008/hYeHB3766SecPXsWV65cwY4dO9C1a9eqek+aNAkPP/wwPvzwQ3Tu3BmJiYnIy8tDTk4OkpOT4e3tjWXLliE8PBzz58+vev+goCBs27EDpX5+OHn2LH799Vc4OztjzZo1GDp0qEm9b7nlFvzzn/+Eo6MjDh48iOPHj8NgMGDz5s3o3bs3Pvjgg6pzhw4dik8++QSenp7Ys2cPsrKycOXKFXz77bfw8/PDokWLqs4NCwvDwoUL4efnh++++w75+fnIzs7Gnj174OXlhY8//hgjRowwqXfPnj2xZcsWGAwGpKen4+DBg3BycsK6deswePBgxMXFVZ07YMAArF27Fs7Ozvj1119x4sQJXLt2Ddu2bUPPnj2xcOHCqnODg4OxbNkyeHl54ccff0ROTg7y8vKwc+dO+Pr6mtR75MiReP/99+Hv748dO3agoKAAZ8+exd69e9G2bVssX74cYWFhiImJMbnWNm3ahLKyMhw7dqyqdWH9+vUYOHCgybU2cOBArFmzBps3b0ZFRQUyMjJQWlqKr776Cj169DCpd0hICOLj49G+fXskJSXh/PnzuHTpEhITE+Hj44PFixdXnRseHo64uDh069YN27dvR2FhITIzM/Hzzz+jTZs2WLFihU19R+zfvx86nQ4rV65EaGioSb1t4Tti2rRpePDBB5v9O2LevHkYPny4ed8R27ahtLQUJ0+e5HdEE74jQkJCEBcXh9mzZ0OncuZlLo2iwPnz59G1a1fs27cPI2us9/Xuu+9i3bp1OH78eK0yBoMBhhrjUYqKWAOguwAAFztJREFUitC9e3ebXBpl1qwSfPihXBolIaEY997LO7ksLjcXWL8e+NOfAF9fa9eGiMguWHLhe3bJKdCpUyc4OTlBr9ebHL948SJ86/nPTafTwcvLy2SzVZyHSRlVyytcuADMncvZvjWKS2fYN8aXlGDCpICrqytCQ0ORmJhocjwxMRHh4eFWqlXL4RgmZV5++WXzCw8dKpfiHjrUchUii1EVW9I8xpeUYMKk0KxZs7BixQqsXLkS6enpePHFF5GVlYVp06ZZu2rNrmYLE+dhqt8//vEPa1eBmglja98YX1KCi+8qFBUVhfz8fPz973/HhQsXMGjQIOzYsQOBgYHWrlqzY5ecMpMmTTK/8IkTwNNPA599BvTrZ7lKkUWoii1pHuNLSrCFqQn++te/4vfff4fBYEBKSgpGjx5t7Sq1CM70rUxKSor5hXU6OaWAyrs4qHmoii1pHuNLSrCFiRrFFiZl/Pz8zC/cowewerWlqkIWpiq2pHmMLynBhIkadf26M4C/AgA8PflPpj7Ozip+N+XlQEEB4O0NuHD5Ga1RFVvSPMaXlGCXHDXKYNABiIerazzc3dllVJ/MzEzzC6elAZ07yz1pjqrYkuYxvqQEEyZqVOUYJnbHNUzVmLaePYGvv5Z70pzWMl6xtWJ8SQkmTM0sPj4eAwYMwG233WbtqpittFQAuARX10vgxPD1q1xGwSze3sDEiXJPmqMqtqR5jC8pwaVRWoglp2dvaT4+JcjLk0ujFBcXw8ODS6PUxWg0wsnJybzCly4BW7cCkycDPj6WrRippiq2pHmMr/3i0ijUompOK0D1i4mJMb9wdjYwfbrck+aoii1pHuNLSrCFqYXYcguTk1MJKirYwkRERLaFLUzUYsrLgYoKa9fCNnABT/vF2No3xpeUYMJEDWJ3nHJTpkwxv3BGBnDPPXJPmqMqtqR5jC8pwYSJGsSESbkffvjB/MJOToCXl9yT5qiKLWke40tKcHpTahATJuX69u1rfuGePYFNmyxXGbIoVbElzWN8SQkmTNQgmTA5A5iCXr24hEBDrqnJLo1GuWifmxtbmTRIVWxJ8xhfUoJdctQgufCuDsBq3H33auh0XBqlPvn5+eYX/u03oG1buSfNURVb0jzGl5RgwtTMbH2m75p/eHFplIYNGTLE/MI9egAbNsg9aY6q2JLmMb6kBBOmZjZ9+nQcO3YMBw4csHZVzCITJgGgBE5OJVwapQHbt283v3CHDsCjj8o9aY6q2JLmMb6kBCeubCG2OnHl9u3AhAklADhxZWNKS0vRpk0b8wpfvgzs2AHcey+TJg1SFVvSPMbXfnHiSmoxHAup3MKFC80v/PvvwBNPyD1pjqrYkuYxvqQEW5haiK22MK1dC0yZwhamZldRIadVd3EBHPl3DBGRJbCFiVoMW5iUU7W8gqMjoNMxWdIoLp1h3xhfUoLfztQgJkzKPffcc+YXPnMGmDxZ7klzVMWWNI/xJSWYMFGDmDApt3XrVvMLV1QABgNXOtYoVbElzWN8SQlO20wNYsKknKq5tnr3BhISLFcZsihbnUeNlGF8SQkmTNQgmTA5AXgIY8cCTly2o14XLlywdhWomTC29o3xJSXYJdfM7GOmbzcAm/Dee5vg5uZm5Rppl9FoNL/wwYOAg4Pck+aoii1pHuNLSjBhamb2MdO3xFypYT3ULGsSEAB8+qnck+aoii1pHuNLSjBhogZxLTnl9uzZY37hTp2AqVPlnjRHVWxJ8xhfUoIJEzVIJkwlABzQu7cDSkpKrFwj7Xr00UfNL3zlCrB1q9yT5qiKLWke40tKMGGiBvEuOeWWLl1qfuHMTODBB+WeNEdVbEnzGF9SgkujtBBbXRpl7FggOZlLozS7GzeAoiLAywtw5s2rRESWwKVRqMWwhUk5VcsrODsDHTowWdIoLp1h3xhfUoIJEzWICZNyL7zwgvmFMzOBxx9nl5xGqYotaR7jS0owYaIGMWFSbvXq1eYXLi8HcnLknjRHVWxJ8xhfUoLt/9QgJkzK3XXXXeYX7tsX2L3bYnUhy1IVW9I8xpeUYAtTM7OPmb6d0KbNvbj33nu5NEoDjh8/bu0qUDNhbO0b40tKMGFqZvYx07cbevVKQEJCApdGaYCquwcPHQI8POSeNId3hto3xpeUYMJE9RKiukuOs3w3rkOHDuYX9vcHYmLknjRHVWxJ8xhfUoIJE9XLYKh+zISpcYcPHza/cOfOwAsvyD1pjqrYkuYxvqQEEyaqV/WA7xLs3esBDw8PLo3SgPHjx5tfuKgI+O47uSfNURVb0jzGl5RgwkT1qnmHnNFYitLSUutVxgasXLnS/MKnTgHjxsk9aY6q2JLmMb6khF0lTEIIzJ07F/7+/nB3d8fYsWNx9OjRRsstXboUQUFBcHNzQ2hoaK2Vq5cvX46xY8fCy8sLDg4OKCgoaK6PoCmcUqBp3nzzTfMLDxoEZGfLPWmOqtiS5jG+pIRdJUxxcXFYuHAhlixZggMHDsDPzw933303rl69Wm+ZjRs3YubMmXjjjTeQmpqKUaNGYfz48cjKyqo6p7S0FOPGjcPrr7/eEh9DM5gwNY2q5RVcXYFu3eSeNIdLZ9g3xpeUsJvFd4UQ8Pf3x8yZM/Hqq68CAAwGA3x9fREbG4vnnnuuznJhYWEICQnBsmXLqo7dcsstmDRpEmJiYkzO3b17N+644w5cuXIF3t7eTaqfLS6+e+AAMHw4AHDxXSUMBgN0Op15hbOygHffBd54AwgIsGzFSDVVsSXNY3ztFxffrUNmZib0ej0iIyOrjul0OowZMwb79u2rs0xZWRlSUlJMygBAZGRkvWVaE7YwNc2iRYvML1xaChw8KPekOapiS5rH+JISdrM0il6vBwD4+vqaHPf19cXZs2frLJOXlwej0Vhnmcr3M5fBYIChxn35hYWFAGS2ayvy8ysfVd8ZV1RUBKPRaJX6aF1kZKT58fX3B374QT62oX8jrYWq2JLmMb72qzKuluhMs9mEaf369SbdbAkJCQAABwcHk/OEELWO3cycMo2JiYnB22+/Xet49+7dVb2vtflzYkUiIrIx+fn5aNeunar3sNmEaeLEiQgLC6t6Xtmao9fr0aVLl6rjFy9erNWCVKlTp05wcnKq1ZrUUBmlZs+ejVmzZlU9LygoQGBgILKyslQHrT633XZbsy7BwvevX1FREbp3747s7OxmG6Nmy7+f5n5/xpbvby7G177fv7CwEAEBARaZzd1mEyZPT094enpWPRdCwM/PD4mJiQgODgYgxyglJycjNja2zvdwdXVFaGgoEhMT8cADD1QdT0xMxP3336+qfjqdrs5BhO3atWu2i9LJyalZB5Tz/Rvn5eXF+Frh/Rlbvr9ajK/9vj8AODqqH7JtswnTzRwcHDBz5kzMnz8fffr0QZ8+fTB//ny0adMGjz32WNV5f/jDH/DAAw8gOjoaADBr1iw88cQTGDZsGEaOHInly5cjKysL06ZNqyqj1+uh1+tx6r+TCqalpcHT09NiWaulTJ8+ne9vxfdvbrb++2nO92ds+f5aZuu/H1t/f4sRdqSiokLMmTNH+Pn5CZ1OJ0aPHi3S0tJMzgkMDBRz5swxORYfHy8CAwOFq6urCAkJEcnJySavz5kzRwCota1atUpx3QoLCwUAUVhYaO7HIw1jfO0XY2vfGF/7Zsn42s08TFpnMBgQExOD2bNnc74PO8T42i/G1r4xvvbNkvFlwkRERETUCLuZuJKIiIiouTBhIiIiImoEEyYiIiKiRjBhIiIiImoEE6YWsHTpUgQFBcHNzQ2hoaHYs2ePtatEFjB37lw4ODiYbH5+ftauFpnpxx9/xIQJE+Dv7w8HBwd89dVXJq8LITB37lz4+/vD3d0dY8eOxdGjR61UW2qqxuL75z//udb1PGLECCvVlpoiJiYGt912Gzw9PdG5c2dMmjQJJ06cMDnHYDBgxowZ6NSpEzw8PDBx4kTk5OQ06ecwYWpmGzduxMyZM/HGG28gNTUVo0aNwvjx45GVlWXtqpEFDBw4EBcuXKja0tLSrF0lMlNJSQmGDBmCJUuW1Pl6XFwcFi5ciCVLluDAgQPw8/PD3XffjatXr7ZwTckcjcUXAMaNG2dyPe/YsaMFa0jmSk5OxvTp07F//34kJibixo0biIyMRElJ9cLxM2fOxLZt2/DFF19g7969KC4uxh//+MemLSaveiYnatDw4cPFtGnTTI71799fvPbaa1aqEVnKnDlzxJAhQ5pW6NAhIR5/XIjbbxfiX/8S4vJlIf72NyGmTRNi0iQhUlObp7LUJADEtm3bqp5XVFQIPz8/sWDBgqpj169fF+3atRMff/yxNapIKtwcXyGEmDJlirj//vutVCOypIsXLwoAVZNQFxQUCBcXF/HFF19UnXPu3Dnh6Ogovv32W8XvyxamZlRWVoaUlBRERkaaHI+MjMS+ffusVCuypIyMDPj7+yMoKAiPPPIIzpw503CBRYuAVauAceOAv/wFeOYZ4JVXgIcfBhISgM8+a5mKU5NkZmZCr9ebXMs6nQ5jxozhtWxHdu/ejc6dO6Nv37545plncPHiRWtXicxQWFgIAFVLl6WkpKC8vNzk+vX398egQYOadP0yYWpGeXl5MBqN8PX1NTnu6+sLvV5vpVqRpYSFhWHt2rX47rvv8Omnn0Kv1yM8PBz5+fl1Fzh9GvD3B5ydgfPngcuXgddfB7p2BfLzgbZtgQkTWvZDkCKV1yuvZfs1fvx4rF+/Hrt27cIHH3yAAwcO4M4774TBYLB21agJhBCYNWsWIiIiMGjQIADy+nV1dUX79u1Nzm3q9Ws3i+9qmYODg8lzIUStY2R7xo8fX/V48ODBGDlyJHr16oU1a9Zg1qxZtQvk5gIPPigf790LREQAISHy+UMPyY00jdey/YqKiqp6PGjQIAwbNgyBgYFISEjA5MmTrVgzaoro6GgcPnwYe/fubfTcpl6/bGFqRp06dYKTk1OtDPbixYu1/lIl2+fh4YHBgwcjIyOj7hPCw4HgYCAvDzhyBBg7tkXrR+arvPuR13Lr0aVLFwQGBtZ/PZPmzJgxA9988w2SkpLQrVu3quN+fn4oKyvDlStXTM5v6vXLhKkZubq6IjQ0FImJiSbHExMTER4ebqVaUXMxGAxIT09Hly5dGj4xKQkQggmTDQkKCoKfn5/JtVxWVobk5GRey3YqPz8f2dnZjV/PZHVCCERHR2Pr1q3YtWsXgoKCTF4PDQ2Fi4uLyfV74cIFHDlypEnXr9PcuXPnWqrSVJuXlxfeeustdO3aFW5ubpg/fz6SkpKwatUqeHt7W7t6pMLLL78MnU4HIQROnjyJ6OhonDx5Ep988knDsY2PB9LSgCVL5Himm33xBfDoo8DXXwNXrgApKcCyZcDhw8CoUYCDgyz/2mvAqVPA0qVAjx5Aly7Axo1AVBSwZ4/sAnzpJcDRERg6tNl+D/aiuLgYx44dg16vxyeffIKwsDC4u7ujrKwM3t7eMBqNiImJQb9+/WA0GvHSSy/h3LlzWL58OVe5twENxdfJyQmvv/46PD09YTQacejQIUydOhXl5eVYsmQJ46tx06dPx/r167F582b4+/ujuLgYxcXFcHJygouLC9zc3HD+/HksWbIEQ4YMQWFhIaZNmwZPT0/ExsbC0VFh25EF7+SjesTHx4vAwEDh6uoqQkJCqm51JNsWFRUlunTpIlxcXIS/v7+YPHmyOHr0aOMF+/cXYsyYhs9ZsUKIwEA57YAQQhiNQgwZIsTcuULk5grRqZMQaWnytWPHhOjWTQiDQT6Pj5evHzkixKefCrF1q7kfsVVJSkoSAGptU6ZMEULIqQXmzJkj/Pz8hE6nE6NHjxZplTEgzWsovqWlpSIyMlL4+PgIFxcXERAQIKZMmSKysrKsXW1SoK64AhCrVq2qOufatWsiOjpadOjQQbi7u4s//vGPTY6vw39/GBG1BL1etgTNmQM01Li7erXcdu+uPvbJJ7JV6fXX5dQEx45Vv9anD/Dpp7Kbb/VqICYGuGmmWyIiMh/vkiNqSbm5gK+vnHepqXr2BAoKZKIkhEyMKo0eDdTsNuASLUREFsWEiaglDRkiW5nMkZkJdOgADBoEHDwI/PnP1a/VfAzIcU5ERGQxvEuOSKsyM4GiIvm4vBxYvlx2x/3lL3IweM2FX7//HkhNtU49iYhaAbYwEWmVj4/sdvPwAH75BZgyBZgxQ762ezcwfz4QECBnCO/aFbjrLuDLL4HFi+X4palTgf/7P3kOERGpwkHfRFpU16BvIiKyGnbJERERETWCCROR1mzcCHz0kZywcvp0a9eGiIjALjkiIiKiRrGFiYiIiKgRTJiIiBqwadMma1eBiDSACRMRUT327duHNWvWWLsaRKQBTJiIiOqxfv16/OlPf7J2NYhIAzjom4ioDjdu3ECfPn1w9OhRtGnTxtrVISIrYwsTEVEdvv32W0RERDBZIiIATJiIiOrE7jgiqoldckRENykuLsbgwYORkZEBZ2cuuUlEbGEiIjv122+/4YknnkBERAS2b9+OK1euYObMmXj++efxwAMP4NChQ/WW3bp1KyZMmFArWdq7dy+efPJJjBw5EgkJCaioqEB8fDxmzJiBadOmISIiAvv37wcAbNiwAc8//zyef/55REREYNeuXc36eYmomQkiIjv01FNPifLycvHOO+8IHx8f8eCDD4qcnByxa9cu4eLiIqKjo+stGxkZKfbv329yzGg0iieeeEJUVFSIl156Sfj4+Ijo6Gixa9euqnOeffZZERQUJN5++22xZcuWquNvvfWW6NChgzAajZb/oETUItjCRER25/Tp0/D394ezszPOnz+Py5cv4/XXX0fXrl2Rn5+Ptm3bYsKECXWWzc3NxdmzZxEWFmZy/MCBAwgODoaDgwPOnz+PS5cu4b777sMdd9xRdY6XlxcyMzPh6+uLyZMnVx339fXF5cuXcenSpeb5wETU7DiGiYjszr59++Du7o7g4GDceuut6NChA3bv3q2o7KJFi3D58mW8/fbbJsd//PFHBAQEoEePHujevTsGDhyIb7/91uSce+65BxkZGTh9+jQcHByqjkdHR2PlypUoKirimCgiG8UWJiKyO+Hh4QgODkZeXh6OHDmCsWPHKi67fv16PPbYY7WOjx49Gj169MCpU6eQk5ODu+++2+T1GzduYN++fRg7dqxJsgQAO3fuxKhRo5gsEdkwJkxEZLeSkpIghFCcMJ08eRJCCPTr16/ecyoHb9/8ngcOHEBxcXGt42lpacjIyMBDDz3UlKoTkcYwYSIiu5WUlASdTocRI0YoOl/J3EtJSUnw9vZGcHCwyfHKLr+bE6YNGzZAp9NVJUyfffaZssoTkaawfZiI7FZSUhJGjBgBNzc3Redv3Lix0bFOu3fvxujRo+HoaPr3ZlJSEnr27ImAgACT41u3bsXEiRPRvn17pKWlIScnp0mfgYi0gS1MRGSX9Ho9jh8/rrg7bv/+/QgICICfn1+956Snp0Ov19d6z/Ly8qrxSzfLy8vDmDFjUFFRgdjYWERHRzfhUxCRVjBhIiK7lJubC19fXzz88MOKzlfSHZefnw9/f3+TKQMAoKCgAB4eHnj88cdrlVm8eDE2b96MqKgoPPXUU+jYsaPyD0FEmsFpBYio1btx4wZ69+6NtLQ0eHp6Wrs6RKRBbGEiolZv586dCAsLY7JERPViwkRErZ6S7jgiat3YJUdErVpJSQkGDBiAU6dOwcXFxdrVISKNYgsTEbVq27Ztw7333stkiYgaxISJiFq15ORkPPnkk9auBhFpHLvkiIiIiBrBFiYiIiKiRjBhIiIiImoEEyYiIiKiRjBhIiIiImoEEyYiIiKiRjBhIiIiImoEEyYiIiKiRjBhIiIiImoEEyYiIiKiRjBhIiIiImoEEyYiIiKiRvw/4jnULdnnlvMAAAAASUVORK5CYII=\n", "text/plain": [ "Graphics object consisting of 20 graphics primitives" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "U1 = 0.02\n", "r1 = find_root(U(r) == U1, 3, 10)\n", "graph += line([(r1, U1), (rmax, U1)], color='red', thickness=1.5) \\\n", " + arrow((9, U1), (10, U1), color='red', arrowsize=4, width=1) \\\n", " + arrow((14, U1), (12, U1), color='red', arrowsize=4, width=1) \\\n", " + text('1', (0.9*rmax, U1 + 0.002), color='red' )\n", "graph += line([(r1, U1), (r1, -0.02)], color='red', linestyle=':') \\\n", " + text(r'$r_{\\rm per}$', (6, -0.014), color='red', fontsize=16)\n", "U2 = 0.04\n", "graph += line([(0, U2), (rmax, U2)], color='brown', thickness=1.5) \\\n", " + arrow((13, U2), (12, U2), color='brown', arrowsize=4, width=1) \\\n", " + arrow((0.8, U2), (0.7, U2), color='brown', arrowsize=4, width=1) \\\n", " + text('2', (0.9*rmax, U2 + 0.002), color='brown' ) \n", "U3 = 0.046\n", "graph += line([(2.2, U3), (rmax, U3)], color='orange', thickness=1.5) \\\n", " + arrow((3, U3), (4, U3), color='orange', arrowsize=4, width=1) \\\n", " + arrow((13, U3), (14, U3), color='orange', arrowsize=4, width=1) \\\n", " + text('3', (0.9*rmax, U3 + 0.002), color='orange' )\n", "U4 = 0.035\n", "r4 = find_root(U(r) == U4, 2, 3)\n", "graph += line([(0, U4), (r4, U4)], color='grey', thickness=1.5) \\\n", " + text('4', (0.9*r4, U4 + 0.002), color='grey' ) \\\n", " + arrow((0.9*r4, U4), (0.95*r4, U4), color='grey', arrowsize=4, width=1) \\\n", " + arrow((0.8, U4), (0.7, U4), color='grey', arrowsize=4, width=1) \n", "show(graph, xmin=0)" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [], "source": [ "graph.save(\"ges_eff_pot_null.pdf\", xmin=0)" ] } ], "metadata": { "kernelspec": { "display_name": "SageMath 9.0.beta9", "language": "sage", "name": "sagemath" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.7.3" } }, "nbformat": 4, "nbformat_minor": 2 }