Skip to content

Instantly share code, notes, and snippets.

@hannorein
Created April 2, 2024 20:05
Show Gist options
  • Select an option

  • Save hannorein/db28bc21a502cef01f554e3b9d9de89e to your computer and use it in GitHub Desktop.

Select an option

Save hannorein/db28bc21a502cef01f554e3b9d9de89e to your computer and use it in GitHub Desktop.
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 3,
"id": "ad75f25a-b3f1-4eab-8b07-53edaf479ebc",
"metadata": {},
"outputs": [],
"source": [
"import rebound\n",
"import numpy as np\n",
"import matplotlib.pyplot as plt"
]
},
{
"cell_type": "code",
"execution_count": 14,
"id": "c2d50e6b-ecb0-4a21-b527-8242fb803097",
"metadata": {},
"outputs": [],
"source": [
"sim = rebound.Simulation()\n",
"sim.add(m=1)\n",
"sim.add(a=1)\n",
"times = np.linspace(0,1e6,10000)\n",
"semias = np.zeros(len(times))\n",
"mass = 1.0 - 0.5*times/times[-1]\n",
"for i in range(len(times)):\n",
" sim.integrate(times[i], exact_finish_time=False)\n",
" semias[i] = sim.particles[1].a\n",
" sim.particles[0].m = mass["
]
},
{
"cell_type": "code",
"execution_count": 19,
"id": "d08504e2-f3fd-4497-82c3-983b6225d29e",
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"[<matplotlib.lines.Line2D at 0x12f411520>]"
]
},
"execution_count": 19,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAjcAAAGwCAYAAABVdURTAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy89olMNAAAACXBIWXMAAA9hAAAPYQGoP6dpAABr/klEQVR4nO3deVzUdf4H8NfMAAPILTLM4AheILc3oplnXmC51uav3KSsdts1K+3Y3MrWarPttMOyttKObS0rLcEob/PIGwaUQwQFnYNLGA4ZYOb7+wMdJREZYBhmeD0fDx678/ke855v6rz4fD/fz0ckCIIAIiIiIgchtnUBRERERJ2J4YaIiIgcCsMNERERORSGGyIiInIoDDdERETkUBhuiIiIyKEw3BAREZFDcbJ1AV3NZDJBrVbD09MTIpHI1uUQERFRGwiCgKqqKigUCojFrffN9Lhwo1aroVQqbV0GERERtUNRURH69u3b6j49Ltx4enoCaLo4Xl5eNq6GiIiI2kKv10OpVJq/x1vT48LN5VtRXl5eDDdERER2pi1DSjigmIiIiByKTcPNypUrMWrUKHh6eiIgIABz5sxBTk7ODY/bsGEDhgwZAldXV0RHR2PLli1dUC0RERHZA5uGm927d2PRokX47bffsHXrVjQ0NGDatGmoqam57jH79+/HXXfdhfvvvx/Hjx/HnDlzMGfOHGRmZnZh5URERNRdiQRBEGxdxGUlJSUICAjA7t27cfPNN7e4z7x581BTU4Pk5GRz25gxYzB06FCsWbPmhu+h1+vh7e2NyspKjrkhIiLqZHtPlUIiFiF+YO9OPa8l39/dasxNZWUlAMDPz++6+xw4cABTp05t1jZ9+nQcOHCgxf0NBgP0en2zHyIiIupc1YZGPLspA3/65CAe/yYN+roGm9XSbcKNyWTCY489hnHjxiEqKuq6+2m1WshksmZtMpkMWq22xf1XrlwJb29v8w/nuCEiIupcu3NLMP2tPfjyt0IAwJRwGSQ2nCi32zwKvmjRImRmZmLv3r2det5ly5Zh6dKl5teXn5MnIiKijqmsbcBLKSex4eg5AIDSzw2vzI3BuEH+Nq2rW4Sbhx9+GMnJydizZ88NZx0MDAyETqdr1qbT6RAYGNji/lKpFFKptNNqJSIiIuDnE1o8uykTJVUGiETAvWND8OT0MLi72D5a2PS2lCAIePjhh7Fx40bs2LED/fv3v+Ex8fHx2L59e7O2rVu3Ij4+3lplEhER0SVl1QY8/NUx/OWLoyipMmBAn17Y8Jd4PD87slsEG8DGPTeLFi3CV199hR9++AGenp7mcTPe3t5wc3MDACxYsABBQUFYuXIlAODRRx/FhAkT8MYbbyAhIQHr16/HkSNH8NFHH9nscxARETk6QRDwY7oa//zxBC7UNkAiFuHPNw/Ao1MGw9VZYuvymrFpuPnggw8AABMnTmzWvnbtWtx7770AgMLCwmarf44dOxZfffUVnn32WfzjH//A4MGDsWnTplYHIRMREVH76fR1eGZjJrZlNQ0LGRLoidfuiEV0X28bV9aybjXPTVfgPDdERERtIwgCNhw5hxdTTqKqrhHOEhEWTx6MhyYMhItT145sseT7u3vcHCMiIqJupai8Fv/YmIFfT5UCAGL7euPVO2IRFnjjVbltjeGGiIiIzBqNJnx+4Cxe/yUHtfVGSJ3EeHxaKBaO6w8nSbeZHq9VDDdEREQEAEgrqsAzGzNwQt00m//oED/8+44Y9PfvZePKLMNwQ0RE1MPp6xrw+s85+OK3sxAEwMvVCU/PDMf/jVJCLLbdTMPtxXBDRETUQwmCgM0qDV5MPomSKgMAYO6wIPwjIRz+HvY7AS7DDRERUQ90prQGz/2QaR4wPMC/F16aE4WxNl46oTMw3BAREfUghkYjPtqdj3d35qG+0QQXJzEenjQIf5kwAFKn7jUZX3sx3BAREfUQ+0+X4tlNmcgvqQEA3DTIHy/OibK7AcM3wnBDRETk4EqrDXh5Sxa+P3YeAODvIcVzieG4NVYBkcj+BgzfCMMNERGRgzKaBKw/XIhXU3NQebEBIhHwp7hgPDE9DN5uzrYuz2oYboiIiBzQ0bPleP7HE8g83zRnTYTcC//6QxSG9fO1cWXWx3BDRETkQIr1dVj5UzY2Hm+6BeXp6oQlU0OxID7YbmYY7iiGGyIiIgdQ32jC2n0FeGf7KdTUGyESAXeOUOLJGWF2PWdNezDcEBER2bldOcV4YfNJ5Jc2PQU1VOmDFbdGIlbpY9vCbIThhoiIyE6dLavBi8lZ2JalA9D0FNTTM4dg7rAgu1w2obMw3BAREdmZ2vpGvL/zND76NR/1jSY4iUW4b1wIFk8ZDC9Xx30Kqq0sDjepqanw8PDATTfdBABYvXo1/vOf/yAiIgKrV6+Gr6/jj8ImIiKyBaNJwHdHz+H1X3JQfGktqPGD/fH87AgMCvC0cXXdh8XDpp988kno9U2PlWVkZODxxx/HrFmzUFBQgKVLl3Z6gURERATsPVWKhHd+xVPfqVBcZYDSzw1r/jQCny8czWDzOxb33BQUFCAiIgIA8N133yExMREvv/wyjh07hlmzZnV6gURERD1ZXnEVXt6SjR3ZxQCaHu1+ZPJgLBgb7DBrQXU2i8ONi4sLamtrAQDbtm3DggULAAB+fn7mHh0iIiLqmNJqA1Zty8X/DhXBaBLgJBbhT2OC8eiUwfDt5WLr8ro1i29L3XTTTVi6dClefPFFHDp0CAkJCQCA3Nxc9O3b16Jz7dmzB7Nnz4ZC0bS2xaZNm254zH//+1/ExsbC3d0dcrkcCxcuRFlZmaUfg4iIqFuqazDig12nMfG1Xfjyt0IYTQKmRcjwy5Kb8c9bIxls2sDicPPee+/ByckJ3377LT744AMEBQUBAH766SfMmDHDonPV1NQgNjYWq1evbtP++/btw4IFC3D//ffjxIkT2LBhAw4dOoQHH3zQ0o9BRETUrQiCgB/SzmPKG7vx79RsVBsaERXkhfV/HoOPFozEgD4eti7RbogEQRBsXQQAiEQibNy4EXPmzLnuPq+//jo++OADnD592tz27rvv4t///jfOnTvX4jEGgwEGg8H8Wq/XQ6lUorKyEl5eXp1WPxERUXsdOVOOF1OykF5UAQCQe7viyelhmDO0Z89XczW9Xg9vb+82fX+3qefm6rE0er2+1R9rio+PR1FREbZs2QJBEKDT6fDtt9+2OpB55cqV8Pb2Nv8olUqr1khERNRW+SXV+Nt/j+KONQeQXlQBdxcJHr8lFDsen4i5w/sy2LRTm3puJBIJNBoNAgICIBaLIRJde7EFQYBIJILRaGxfIW3ouQGADRs2YOHChairq0NjYyNmz56N7777Ds7OLU9axJ4bIiLqbrSVdXh7+yl8c6RpsLBYBMwbpcSSW0IR4Olq6/K6JUt6btr0tNSOHTvg5+dn/v8thZuucPLkSTz66KNYvnw5pk+fDo1GgyeffBIPPfQQPvnkkxaPkUqlkEp71oJhRETUPVXWNuCD3aexdl8BDI0mAMCUIQF4ckYYhgTyF+7OYldjbu655x7U1dVhw4YN5ra9e/di/PjxUKvVkMvlN3wfS5IfERFRZ7hYb8Ta/QVYs+s09HWNAICRwb74+8whGBXiZ+Pq7EOnj7m52j//+U+YTKZr2isrK3HXXXdZejqL1NbWQixuXrJE0jSBUTfJaERERGYNRhO+/O0sJry2E6+m5kBf14ghgZ749N6R2PBQPIONlVg8id8nn3yCX375BV9++SUGDBgAANi1axcWLFiAwMBAi85VXV2NvLw88+uCggKkpaXBz88P/fr1w7Jly3D+/Hl8/vnnAIDZs2fjwQcfxAcffGC+LfXYY49h9OjRUCgUln4UIiIiqzCZBCRnaPDGLzk4W9Y08W1fXzc8Pi0Ut8YGQcKBwlZlcbhRqVT4y1/+gqFDh+KNN95Abm4u3n77bTz55JNYsWKFRec6cuQIJk2aZH59eW2qpKQkrFu3DhqNBoWFhebt9957L6qqqvDee+/h8ccfh4+PDyZPnox///vfln4MIiKiTicIAvacKsWrqdk4oW56gtjfwwWLJw/GXaP7wcXJ4hsm1A7tHnPzj3/8A6+88gqcnJzw008/YcqUKZ1dm1VwzA0REVnDb/lleGtrLg4WlAMAPKRO+MvNA7Dwpv7oJbW4L4F+p9Oflvq9d999F2+//TbuuusuHD16FI888gi++uorxMbGtqtgIiIie3XkTDne2paLfXlNSwG5OImRFB+Mv04cBD8ulWATFoebGTNm4MiRI/jss89wxx134OLFi1i6dCnGjBmDFStW4KmnnrJGnURERN3KscILeGtrLn49VQoAcJaIMG+UEosmDYLc283G1fVsFocbo9EIlUplHsDr5uaGDz74AImJiXjggQcYboiIyKGpzlXgra252JlTAgBwEovwx5F9sWjSIPT1dbdxdQR08jw3paWl8Pf376zTWQXH3BARUXtknq/Eqm2nsC1LBwCQiEW4fXgQFk8eDKUfQ421WX3MzfV092BDRERkqSyNHm9vO4XUE1oAgFgEzBkWhEcmD0aIfy8bV0ctaddtqbfeegvffPMNCgsLUV9f32x7eXl5pxVHRERkK2lFFXhvR565p0YkAm6NVeCRKYMxsI+Hjauj1lgcblasWIGPP/4Yjz/+OJ599lk888wzOHPmDDZt2oTly5dbo0YiIqIuc6igHO/uOGUeKCwSAbOi5XhsymAMlnnauDpqC4vH3AwcOBDvvPMOEhIS4OnpibS0NHPbb7/9hq+++spatXYKjrkhIqLfEwQBe/NK8e6OPBy6NE+NRCzCnKFB+OvEgRgUwJ4aW7PqmButVovo6GgAgIeHByorKwEAiYmJeO6559pRLhERkW0IgoBtWcV4b2ce0osqAAAuEjHuGNkXf50wkAOF7ZTF4aZv377QaDTo168fBg4ciF9++QXDhw/H4cOHIZVKrVEjERFRpzKaBKRmavHujlPI1lYBAFydxbhrdD/8+eYBnKfGzlkcbv7whz9g+/btiIuLw+LFi/GnP/0Jn3zyCQoLC7FkyRJr1EhERNQpGo0m/JiuxuqdeThdUgMA6OUiwT3xIbj/pv7o48lf0h1Bh+e5+e2337B//34MHjwYs2fP7qy6rIZjboiIep66BiM2HD2Hj/acRlH5RQCAl6sT7hvXH/eNC4GPO5dJ6O66dJ6bMWPGYMyYMR09DRERUaerqK3HFwfOYt3+MyiraZq6pHcvF9w/vj/uGRMMT1dnG1dI1sBlSomIyOGoKy7i418LsP5wIWrrjQCAIB83PDi+P+4cpYS7C7/+HBn/6xIRkcPI0Vbhwz2n8WOaGo2mplEX4XIvPDRhABKi5XCSiG1cIXUFhhsiIrJrgiDg8JkLWLP7NHZkF5vb4wf0xkMTB+Lmwf4QiUQ2rJC6mkXhxmg0Yt++fYiJiYGPj4+VSiIiIrqxRqMJP2Vq8fHeAvMcNSIRMDMqEH+5eSBilT42rY9sx6JwI5FIMG3aNGRlZTHcEBGRTVRebMDXhwvx2f6zOF/R9OSTi5MYd4zoiz+PH8DFLMny21JRUVHIz89H//79rVEPERFRiwrLarF2fwG+OVyEmkuDhP09XHDPmBDMH9MP/h6co4aaWDyy6qWXXsITTzyB5ORkaDQa6PX6Zj+W2LNnD2bPng2FQgGRSIRNmzbd8BiDwYBnnnkGwcHBkEqlCAkJwaeffmrpxyAiIjsgCAKOnCnHQ18cxcTXd2LtvjOoqTciVOaBV2+Pwd6/T8ajUwcz2FAzFvfczJo1CwBw6623NhugJQgCRCIRjEZjm89VU1OD2NhYLFy4EHPnzm3TMXfeeSd0Oh0++eQTDBo0CBqNBiaTybIPQURE3VrDpfE0n1w1ngYAJoT2wQPj++OmQRwkTNdncbjZuXNnp735zJkzMXPmzDbvn5qait27dyM/Px9+fn4AgJCQkE6rh4iIbOvyeJp1+85AXVkHoGk8zdxhQVh4U3+EyjxtXCHZA4vDzYQJE6xRR5v8+OOPGDlyJF599VV88cUX6NWrF2699Va8+OKLcHNreZEzg8EAg8Fgfm3prTMiIrK+HG0VPjtwBhuPncfFBo6noY5p1zw3FRUV+OSTT5CVlQUAiIyMxMKFC+Ht7d2pxf1efn4+9u7dC1dXV2zcuBGlpaX429/+hrKyMqxdu7bFY1auXIkVK1ZYtS4iIrJco9GErSd1+OzAGfyWX25uD5N54v6b+uPWoQq4OktsWCHZK4sXzjxy5AimT58ONzc3jB49GgBw+PBhXLx4Eb/88guGDx/evkJEImzcuBFz5sy57j7Tpk3Dr7/+Cq1Waw5S33//Pe644w7U1NS02HvTUs+NUqnkwplERDZSVm3A+sNF+PK3s9BcuvUkEYswLUKGpLEhiOvvx/E0dA2rLpy5ZMkS3HrrrfjPf/4DJ6emwxsbG/HAAw/gsccew549e9pXdRvI5XIEBQU16yEKDw+HIAg4d+4cBg8efM0xUqkUUim7M4mIbC29qAKfHTiD5HQN6o1ND4L49XLBXaOVmB8XDIVPy8MLiCxlcbg5cuRIs2ADAE5OTnjqqacwcuTITi3u98aNG4cNGzaguroaHh4eAIDc3FyIxWL07dvXqu9NRESWMzQasSVDg8/2n0XaVU89xfb1RtLYEMyKlvPWE3U6i8ONl5cXCgsLMWTIkGbtRUVF8PS0bBR7dXU18vLyzK8LCgqQlpYGPz8/9OvXD8uWLcP58+fx+eefAwDuvvtuvPjii7jvvvuwYsUKlJaW4sknn8TChQuvO6CYiIi6nrriIv53qBD/O1SI0up6AICLRIyEGDmSxoZgKJdGICuyONzMmzcP999/P15//XWMHTsWALBv3z48+eSTuOuuuyw615EjRzBp0iTz66VLlwIAkpKSsG7dOmg0GhQWFpq3e3h4YOvWrVi8eDFGjhyJ3r17484778RLL71k6ccgIqJOZjQJ2J1bjK8OFmJHdjEuLcqNQC9XzI/rh/8b3Q99PDlMgKzP4gHF9fX1ePLJJ7FmzRo0NjYCAJydnfHXv/4Vr7zySrcf32LJgCQiIrqxYn0dvj5chPWHi8xrPQHAmAF+uGdMCKZFyuAssXhCfKJmLPn+tjjcXFZbW4vTp08DAAYOHAh3d/f2nKbLMdwQEXWcySRg3+lS/Pe3QmzL0qHxUjeNt5sz7hjRF3eN7odBAR42rpIciVWflrrM3d0d0dHR7T2ciIjsUFm1ARuOnsP/DhXibFmtuX1ksC/ujuvHAcLULbQp3MydOxfr1q2Dl5fXDdeA8vDwQGRkJB566CGrT+pHRETWJwgCfssvx1eHCpGaqUGDsamXxlPqhLnDg3B3XDDCArksAnUfbQo33t7e5gmVbhRYDAYD1qxZg3379uHHH3/seIVERGQTOn0dvj16DhuOFOHMVb00sX29MT8uGImxcri7tPsGAJHVtHvMTWtOnjyJUaNGoaamprNP3WEcc0NEdH0NRhN2ZBfjm8NF2JVbAuOlsTS9XCS4dagC8+OCERXEXnnqel0y5qY1YWFh2L9/vzVOTUREVnC6pBrfHC7Cd8fOo7T6ypI1I4N9cecoJRKi5eglZS8N2Yd2/Uk9cuQIvvnmGxQWFqK+vr7Ztu+//x4SiQSxsbGdUiAREVlHbX0jUlQafHOkCIfPXDC3+3u44PbhffHHkUo+8UR2yeJws379eixYsADTp0/HL7/8gmnTpiE3Nxc6nQ5/+MMfrFEjERF1EkEQkFZUgW+OFGFzugbVhqb5ysQiYGJYAOaNUmLykADOS0N2zeJw8/LLL+Ott97CokWL4Onpibfffhv9+/fHX/7yF8jlcmvUSEREHaSuuIiNx8/j+2PncLrkynjI4N7uuHOkErcP74tAb1cbVkjUeSwON6dPn0ZCQgIAwMXFBTU1NRCJRFiyZAkmT56MFStWdHqRRERkudr6RqRmavHdsXPYf7oMlx8fcXUWY2aUHHeOVCKuvx/EYpFtCyXqZBaHG19fX1RVVQEAgoKCkJmZiejoaFRUVKC2tvYGRxMRkTWZTAJ+KyjDd0fP46dMDWrrjeZtcf39cPvwvpgZHQhPV2cbVklkXRaHm5tvvhlbt25FdHQ0/vjHP+LRRx/Fjh07sHXrVkyZMsUaNRIR0Q3kl1Tj+2PnsfH4+WbrOwX3dsfcYX0xd3gQlH72sUwOUUdZHG7ee+891NXVAQCeeeYZODs7Y//+/bj99tvx7LPPdnqBRETUssraBmxWqfH9sXM4VlhhbveUOiExVo7bh/fFiGBf8ySsRD2FVSbx6844iR8R2bO6BiN25ZTgx/Tz2HayGPVGE4Cmp51uDu2D24f3xS0RMq7vRA6nSybxKy4uRnFxMUwmU7P2mJiY9p6SiIhaYDQJOJhfhh/S1NiSqUFVXaN525BAT9w+vC9uG6pAgBefdiIC2hFujh49iqSkJGRlZeH3nT4ikQhGo/E6RxIRUVsJgoDM83psSjuPzelqFFddmTU40MsVtw5V4NZYBSIVXrztRPQ7FoebhQsXIjQ0FJ988glkMhn/UhERdaKC0hr8kHYeP6apkV96ZT4abzdnzIqW47ahCowO4ePbRK2xONzk5+fju+++w6BBg6xRDxFRj1Osr8NmlQY/pp1H+rlKc7ursxhTw2W4bWgQbg71h9SJ42iI2sLicDNlyhSkp6cz3BARdUBFbT1+OaHDj+lq7D9dikuLb0MiFuGmQf64bagC0yID4cHFKoksZvHfmo8//hhJSUnIzMxEVFQUnJ2bTwR16623dlpxRESOpLK2AT+f1GJLhgZ7T5Wi0XRl3OKIYF/cNlSBWdFy+HtIbVglkf2zONwcOHAA+/btw08//XTNNksHFO/ZswevvfYajh49Co1Gg40bN2LOnDltOnbfvn2YMGECoqKikJaW1ub3JCLqSpUXG7D1pA4pKjX25pWiwXgl0ITLvZAYI8etsQpOsEfUiSwON4sXL8af/vQnPPfcc5DJZB1685qaGsTGxmLhwoWYO3dum4+rqKjAggULMGXKFOh0ug7VQETU2fR1Ddh2UocUlQZ7TpU0CzRDAj2REC3HrBg5BvbxsGGVRI7L4nBTVlaGJUuWdDjYAMDMmTMxc+ZMi4976KGHcPfdd0MikWDTpk2t7mswGGAwXHmEUq/XW/x+REQ3UlXXgO1ZxUhWabAnt8Q8uR4AhMo8kBCtQEJMIAYFeNqwSqKeweJwM3fuXOzcuRMDBw60Rj03tHbtWuTn5+PLL7/ESy+9dMP9V65cyZXKicgqqg2N2J7V1EOzK7cE9Y1XAs2gAA8kxsiREC3HYBkDDVFXsjjchIaGYtmyZdi7dy+io6OvGVD8yCOPdFpxv3fq1Ck8/fTT+PXXX+Hk1LbSly1bhqVLl5pf6/V6KJVKa5VIRA6uvKYe27J0+DlTi1/zSpsFmgF9eiExRoHEGDlCGWiIbKZdT0t5eHhg9+7d2L17d7NtIpHIauHGaDTi7rvvxooVKxAaGtrm46RSKaRSPnlARO2nqbyIX07okJqpxcGCMlz1kBMG+PdCQowcCTFyhMk8ObEpUTdgcbgpKCiwRh03VFVVhSNHjuD48eN4+OGHAQAmkwmCIMDJyQm//PILJk+ebJPaiMjx5JdU4+cTOqSe0CK9qKLZtkiFF2ZEBmJGVCAGBXgw0BB1Mx2aHWrfvn0YOXJkl/SMeHl5ISMjo1nb+++/jx07duDbb79F//79rV4DETkuQRBwUqPHz5lapJ7QIldXbd4mEgEjg30xPTIQ0yMD+dg2UTfXoXAzc+ZMpKWlYcCAAe06vrq6Gnl5eebXBQUFSEtLg5+fH/r164dly5bh/Pnz+PzzzyEWixEVFdXs+ICAALi6ul7TTkTUFiaTgGOFF/DziaZAU1R+0bzNSSzC2EH+mB4pwy0RMgR4csVtInvRoXDz+1XBLXXkyBFMmjTJ/PrywN+kpCSsW7cOGo0GhYWFHXoPIqKr1TUYsS+vFNuydNiWVYySq1bbdnUWY0JoH8yICsTkMBm83Z1bORMRdVcioQMJxdPTE+np6e3uubEFvV4Pb29vVFZWwsvLy9blEFEXKKkyYEd2U5j59VQJ6hquPOHk6eqEqeEyTI8MxITQPnBz4eKURN2RJd/fHeq5+fDDDztlMj8ios4kCAJOFVdj60kdtmXpkFZUgat/jQvyccPU8ABMCZdhzIDecHES265YIup0HQo3d999d2fVQUTUIQ1GEw6fKce2k8XYlqVDYXlts+0xfb0xNVyGqeEyhMv5yDaRI2tTuJk7dy7WrVsHLy+vG64B9f3333dKYUREN1JZ24Ddp0qwPUuHndnF0Nc1mre5OIlx0yB/TA2XYUp4AGReHBBM1FO0Kdx4e3ubf8vx9va2akFERNcjCAJydFXYkV2MXdklOFp4AcarZtTr3csFk4cEYGqEDOMH+8PdpUOd00Rkpzo0oNgecUAxkX2pMTRi/+mypkCTUwxNZV2z7YMDPDAlXIZbIgIwVOkLiZi3m4gcUZcNKCYisoaC0hrszC7GzpxiHMwvb7bCtquzGGMH+mPSkABMDO3DCfWI6BoWh5uysjIsX74cO3fuRHFxMUwmU7Pt5eXlnVYcEfUMhkYjDhWUX+qdKUFBaU2z7Uo/N0wOC8DEIQGIH9Abrs58XJuIrs/icHPPPfcgLy8P999/P2QyGZ84IKJ2KSqvxZ5TJdiVU4J9eaWorTeatzmJRRjd3w+ThwRgYlgABvbpxX9riKjNLA43v/76K/bu3YvY2Fhr1ENEDqrG0IiDBWXYk1uKPbklyP9d70yApxSTwgIwaUgfjBvkD09Xzg5MRO1jcbgZMmQILl68eOMdiahHu7wQ5eUwc+RsORqMV55fkIhFGN7PBzcP7oNJQwIQqfBi7wwRdQqLw83777+Pp59+GsuXL0dUVBScnZv/dsUnkIh6rtJqA/aeagoze06VorTa0Gx7X1833BzaBzcP7oOxg3rDi70zRGQFFocbHx8f6PV6TJ48uVm7IAgQiUQwGo3XOZKIHE19ownHCi9cCjMlyDyvb7bdzVmC+IG9cfNgf9wc2gf9/Tl2hoisz+JwM3/+fDg7O+Orr77igGKiHubymk17T5Vi/+lSHDhdhpr65r/QRMi9LvXO+GNEiC+kTnyyiYi6lsXhJjMzE8ePH0dYWJg16iGibuZ8xUXsyyvF/rxS7DtdhpKq5reaevdywfhLPTM3DfZHgCeXOSAi27I43IwcORJFRUUMN0QOqrK2AQfyS7E3rxT78squmXPG1VmMUSF+GDfIHzcN8keE3AtizgpMRN2IxeFm8eLFePTRR/Hkk08iOjr6mgHFMTExnVYcEVlfXYMRR85cwN68pltNGecrcfWiLGIREKv0wbiB/hg3yB/Dg314q4mIujWL15YSi8XXnkQkspsBxVxbino6o0lAxvlK7Msrxb68Uhw5ewH1jc1nGh8c4IFxg5rCTNwAPz7VREQ2Z9W1pQoKCtpdGBF1PaNJwEm1Hr/ll+G3/DIcKihHlaGx2T6BXq6XwkxvjBvkD5kXx80Qkf2yONwEBwdbow4i6iRGk4AszZUwc7CgHFV1zcOMp6sTxg7sbe6dGcBHtInIgbRrVfAvvvgCa9asQUFBAQ4cOIDg4GCsWrUK/fv3x2233dbm8+zZswevvfYajh49Co1Gg40bN2LOnDnX3f/777/HBx98gLS0NBgMBkRGRuKf//wnpk+f3p6PQeQQmoeZchwqKIP+92FG6oTR/f0QP7A3xgzojXC5FyQcBExEDsricPPBBx9g+fLleOyxx/Cvf/3LPMbGx8cHq1atsijc1NTUIDY2FgsXLsTcuXNvuP+ePXtwyy234OWXX4aPjw/Wrl2L2bNn4+DBgxg2bJilH4XILplMArK0evyWX26+zVR5saHZPh6Xw8yApjAToWCYIaKew+IBxREREXj55ZcxZ84ceHp6Ij09HQMGDEBmZiYmTpyI0tLS9hUiEt2w56YlkZGRmDdvHpYvX96m/TmgmOxNo9GELE0VDp0px8FLt5laCjOjQnwxZkBvxA/sjQi5F5wk1w7+JyKyV1YfUNxSL4lUKkVNTU0LR1iPyWRCVVUV/Pz8rruPwWCAwXBl0jG9Xn/dfYm6g7oGI9KKKnC4oByHzpTj2NkL18wC3MtFglH9/ZrCzIDeiFQwzBARXWZxuOnfvz/S0tKuGVicmpqK8PDwTiusLV5//XVUV1fjzjvvvO4+K1euxIoVK7qwKiLLVF5swLGzF3DoTDkOFZQj41wl6o3NH832dHXCyGBfjO7fG2MG+CE6yJthhojoOiwON0uXLsWiRYtQV1cHQRBw6NAh/O9//8PKlSvx8ccfW6PGFn311VdYsWIFfvjhBwQEBFx3v2XLlmHp0qXm13q9HkqlsitKJGpRsb4Oh86UX+qZuYBsrR6/vzkc4CnFqP5+GB3ih1EhfggL9OSYGSKiNrI43DzwwANwc3PDs88+i9raWtx9991QKBR4++238X//93/WqPEa69evxwMPPIANGzZg6tSpre4rlUohlUq7pC6i3xMEAWfLas1h5vCZcpwpq71mv5De7hgV4ofR/Zt++vm589FsIqJ2atej4PPnz8f8+fNRW1uL6urqVntOOtv//vc/LFy4EOvXr0dCQkKXvS9RW9Q3mnBSo8fRsxfMt5p+v9CkSASEB3phdP+mXplRIb4I4KR5RESdxuJwc/HiRQiCAHd3d7i7u6OkpASrVq1CREQEpk2bZtG5qqurkZeXZ35dUFCAtLQ0+Pn5oV+/fli2bBnOnz+Pzz//HEDTraikpCS8/fbbiIuLg1arBQC4ubnB29vb0o9C1GFl1QYcK6wwh5n0cxUw/G4pAxeJGDF9vZtuM/X3w4hgXy5nQERkRRY/Cj5t2jTMnTsXDz30ECoqKhAWFgYXFxeUlpbizTffxF//+tc2n2vXrl2YNGnSNe1JSUlYt24d7r33Xpw5cwa7du0CAEycOBG7d+++7v5twUfBqb1MJgGniqtx9OyFpjBTeOGaFbMBwNfdGSOCfTE82Bcj+vkiVukDV2cuNElE1BGWfH9bHG78/f2xe/duREZG4uOPP8a7776L48eP47vvvsPy5cuRlZXVoeKtjeGG2qra0Ii0S70yRwsv4HjhhWuWMQCaFpk0h5lgXy5lQERkBVad56a2thaenp4AgF9++QVz586FWCzGmDFjcPbs2fZVTGRjlwf+Hi+6cKlnpgI5Wj1Mv4v+7i4SDFX6mMPMcKUvvN15i4mIqDuxONwMGjQImzZtwh/+8Af8/PPPWLJkCQCguLiYPSFkNypq65FWVGH+SS+qwIXahmv26+vrhhGXemSG9/PFkEBPzi9DRNTNWRxuli9fjrvvvhtLlizBlClTEB8fD6CpF4frO1F3VN9oQrZWj7SiChwvbAozLY2VcZGIEaHwwsjLYSbYFzI+xUREZHcsHnMDAFqtFhqNBrGxsRCLm36LPXToELy8vDBkyJBOL7IzccyNYxMEAecuXMTxogqkFVYgregCMtV61P/uCSagaW6ZoUqfpp9+vgiXe0LqxIG/RETdkVXH3ABAYGAgAgMDm7WNHj26Paci6hB9XQPSzUGm6aespv6a/bzdnK8KMj4Y2tcHvr1cbFAxERFZW7vCDZEt1DUYkaXRI+N8JdKLKpFWdAGnS669veQsESFc7oWhSh8M6+eDoUpfhPTmjL9ERD0Fww11S/WNJuTqqqA6V4mM8xVIL6pErq4Kjb9/fAmA0s8NQ5W+5p6ZSIUX55UhIurBGG7I5owmAXnF1VCdq4DqXCVU5yuRpWl5nEzvXi6I6euN6L4+iO3rjVilD/w9uHYYERFdwXBDXcpkElBQVoOMc5VNQeZcBU6o9bjYYLxmXy9XJ8T09UFMX29zoFF4u/L2EhERtYrhhqzGZBJwtrwWJ9SV5jCTeb4SVYZrZ/nt5SJBVJB3s14ZroxNRETtwXBDnaLRaMLpkhpknq9EproSJ9R6nFTrUd1CkJE6iRGp8GrWKzPA3wNiMYMMERF1HMMNWayuwYhcXRUyz+txQl2JTLUe2Rr9NathA4CLkxjhci9EKrwQ29cbMX19MDjAg7P8EhGR1TDcUKuqDY3I0uhx4nxTiMk8X4m84uoWn1rykDohQu6FyCAvRCm8ERnkhYF9PODMIENERF2I4YbMLtTU46RGf+nWUlOvTEFpDVqaw9rX3RlRQd6IVHgjUuGFqCBvBPu589YSERHZHMNND2Q0CThTVoMsTdO4mCyNHlmaKmj1dS3uH+jliqggr2ZBRs6nloiIqJtiuHFwVXUNyNZWXQowepzUVCFHq0ddw7XjYwCgn587ooO8EXEpxEQqvDiPDBER2RWGGwdxecHIk5dDjFqPLK0eReUXW9zf1VmMsEAvRMg9ES73QoTcC2GBnvB0de7iyomIiDoXw40dulhvRI7uSm9MlkaPbE1Vi/PHAIDc2xXhci+EXwoy4XIvhPTuBQnHxxARkQNiuOnGjCYBheW1yNHqka2tQq6uCtnaKpwprUELDyvBRSLGoAAPc5CJUHghPNCLq18TEVGPwnDTDQiCgOIqA3K0VcjRVpmDzKniquuOjfH3cDH3wlzukeFj10RERDYON3v27MFrr72Go0ePQqPRYOPGjZgzZ06rx+zatQtLly7FiRMnoFQq8eyzz+Lee+/tkno7g76uAbnaKuToqsxhJkdXhYrahhb3lzqJESrzRFigJ8Iu/e+QQE8EeLl2ceVERET2wabhpqamBrGxsVi4cCHmzp17w/0LCgqQkJCAhx56CP/973+xfft2PPDAA5DL5Zg+fXoXVNx2hkYjThfXmG8l5Wj1yNVV43xFywN8xSIgxL8XhgR6IkzmhbBAD4QFeqGfnzvHxhAREVnApuFm5syZmDlzZpv3X7NmDfr374833ngDABAeHo69e/firbfeum64MRgMMBgM5td6vb5jRV9HabUBXxw4i7ySauRoq1BQWgNjSwNj0DTAN1TW1ANzuVdmUIAHXJ0lVqmNiIioJ7GrMTcHDhzA1KlTm7VNnz4djz322HWPWblyJVasWGHlyoBGo4C3t59q1ubl6oQhgV4IvdQLEyZrurXk7c7HrYmIiKzFrsKNVquFTCZr1iaTyaDX63Hx4kW4ubldc8yyZcuwdOlS82u9Xg+lUtnptcm8pJgf1w/Bvd3NvTGBXpzFl4iIqKvZVbhpD6lUCqnU+jPsikQi/OsP0VZ/HyIiImqdXT03HBgYCJ1O16xNp9PBy8urxV4bIiIi6nnsKtzEx8dj+/btzdq2bt2K+Ph4G1VERERE3Y1Nw011dTXS0tKQlpYGoOlR77S0NBQWFgJoGi+zYMEC8/4PPfQQ8vPz8dRTTyE7Oxvvv/8+vvnmGyxZssQW5RMREVE3ZNNwc+TIEQwbNgzDhg0DACxduhTDhg3D8uXLAQAajcYcdACgf//+SElJwdatWxEbG4s33ngDH3/8cbeb44aIiIhsRyQIQsuTsTioyspK+Pj4oKioCF5eXrYuh4iIiNrg8tPOFRUV8Pb2bnVfh39a6veqqqoAwCqPgxMREZF1VVVV3TDc9LieG5PJBLVaDU9Pz06fg+ZyqmSvkHXxOncNXueuwevcdXitu4a1rrMgCKiqqoJCoYBY3Pqomh7XcyMWi9G3b1+rvoeXlxf/4nQBXueuwevcNXiduw6vddewxnW+UY/NZXb1KDgRERHRjTDcEBERkUNhuOlEUqkUzz//fJcs99CT8Tp3DV7nrsHr3HV4rbtGd7jOPW5AMRERETk29twQERGRQ2G4ISIiIofCcENEREQOheGGiIiIHArDDRERETkUhhsiIiJyKAw3RERE5FAYboiIiMihMNwQERGRQ2G4ISIiIofCcENEREQOheGGiIiIHArDDRERETkUhhsiIiJyKAw3RERE5FAYboiIiMihMNwQERGRQ2G4ISIiIofCcENEREQOheGGiIiIHArDDRERETkUhhsiIiJyKAw3RERE5FAYboiIiMihMNwQERGRQ2G4ISIiIofiZOsCuprJZIJarYanpydEIpGtyyEiIqI2EAQBVVVVUCgUEItb75vpceFGrVZDqVTaugwiIiJqh6KiIvTt27fVfXpcuPH09ATQdHG8vLxsXA0RERG1hV6vh1KpNH+Pt6bHhZvLt6K8vLwYboiIiOxMW4aUcEAxERERORSbhps9e/Zg9uzZUCgUEIlE2LRp0w2P2bVrF4YPHw6pVIpBgwZh3bp1Vq+TiIiI7IdNw01NTQ1iY2OxevXqNu1fUFCAhIQETJo0CWlpaXjsscfwwAMP4Oeff7ZypURERGQvbDrmZubMmZg5c2ab91+zZg369++PN954AwAQHh6OvXv34q233sL06dNbPMZgMMBgMJhf6/X6jhV9HYIg4P8++g2RCm8kxsoxTOnDR82JiIhswK7G3Bw4cABTp05t1jZ9+nQcOHDgusesXLkS3t7e5h9rPQaeeV6PgwXl+HRfAea+vx83/XsnVm7JgupcBQRBsMp7EhER0bXsKtxotVrIZLJmbTKZDHq9HhcvXmzxmGXLlqGystL8U1RUZJXaBss88NE9I3DbUAV6uUhwvuIiPtyTj1vf24eJr+/Cq6nZOKnWM+gQERFZmcM/Ci6VSiGVSq3+Pq7OEkyLDMS0yEDUNRixM7sYyRkabM/S4WxZLd7fdRrv7zqNAX16ITFajsRYBUJlN35Wn4iIiCxjV+EmMDAQOp2uWZtOp4OXlxfc3NxsVNW1XJ0lmBktx8xoOWrrG7E9qxjJKjV25pQgv6QG7+zIwzs78hAq80BCtAKJsXIM7ONh67KJiIgcgl2Fm/j4eGzZsqVZ29atWxEfH2+jim7M3cUJs2MVmB2rQLWhEdtO6pCsUmNPbilyddXI1eXirW25CJd7ITFGjsQYOYJ797J12URERHZLJNhwEEh1dTXy8vIAAMOGDcObb76JSZMmwc/PD/369cOyZctw/vx5fP755wCaHgWPiorCokWLsHDhQuzYsQOPPPIIUlJSrvu01O/p9Xp4e3ujsrLSpjMUV15swNZLQWfvqVI0mq78Z4gO8kZCjBwJ0XIo/dxtViMREVF3Ycn3t03Dza5duzBp0qRr2pOSkrBu3Trce++9OHPmDHbt2tXsmCVLluDkyZPo27cvnnvuOdx7771tfs/uEm6uVlFbj59PaJGs0mD/6TIYrwo6Q5U+SIyRY1a0HAqf7nPrjYiIqCvZTbixhe4Ybq5WVm1A6gktktM1+K2gDFf/1xkZ7GsOOgFerrYrkoiIqIsx3LSiu4ebqxVX1eGnDC1SVBocPltuDjoiETA6xA+JsQrMjAqEv4f1nwYjIiKyJYabVthTuLmatrIOKRkapKjUOFZYYW4Xi4D4gb2RGKPAjMhA+PZysV2RREREVsJw0wp7DTdXO3ehFlsyNEhRaZB+rtLcLhGLMG6QPxJj5JgeEQhvd2cbVklERNR5GG5a4Qjh5mqFZbVIzlAjRaXBCfWVdbOcJSKMH9wHiTFy3BIhg6crgw4REdkvhptWOFq4uVp+STVSVBokqzTI0VWZ212cxJgQ2hR0pobL0EtqV9MbERERMdy0xpHDzdVO6aqQrNIgWaXG6ZIac7ursxiThwQgMUaBSWEBcHOR2LBKIiKitmG4aUVPCTeXCYKAbG3VpR4dNc6U1Zq3ubtIMCVchsQYOSaE9oGrM4MOERF1Tww3rehp4eZqgiDghFpv7tE5d+HKSuoeUifcEtEUdG4a7A+pE4MOERF1Hww3rejJ4eZqgiAg/VwlUlRNg5HVlXXmbZ6uTpgeGYjEGDnGDfKHs0Rsw0qJiIgYblrFcHMtk0nA8aIL2JyuwZYMDYqrDOZtPu7OmBEZiMQYBcYM8IMTgw4REdkAw00rGG5aZzIJOHymHMkqDX7K1KC0ut68rXcvF8yIago6o/v7QSIW2bBSIiLqSRhuWsFw03aNRhMOFZRjs0qD1EwNLtQ2mLf18ZQiIVqOhBg5RvTzhZhBh4iIrIjhphUMN+3TYDThwOkyJKvUSM3UQl/XaN4m93bFrEtBZ5jSByIRgw4REXUuhptWMNx0XH2jCfvySrFZpcbWEzpUGa4EnSAfNyTGNAWd6CBvBh0iIuoUDDetYLjpXHUNRuzJLUFKhgbbTupQU280bwvu7W6+dRUh92LQISKidmO4aQXDjfXUNRixM7sYyRkabM/Soa7BZN42wL8XEmPkSIxVIFTmacMqiYjIHjHctILhpmvU1jdie1YxklVq7MwpQX3jlaATKvNAQrQCibFyDOzjYcMqiYjIXjDctILhputVGxqx7aQOySo19uSWot54JeiEy72aenRi5Aju3cuGVRIRUXfGcNMKhhvbqrzYgK2Xgs7eU6VoNF354xcd5I2EGDkSouVQ+rnbsEoiIupuLPn+7hbTza5evRohISFwdXVFXFwcDh06dN19Gxoa8MILL2DgwIFwdXVFbGwsUlNTu7Ba6ghvN2fcMaIv1t03GkeenYp/3x6N8YP9IRGLkHG+Eq/8lI3xr+7EnNX78PGv+VBXXLzxSYmIiK5i856br7/+GgsWLMCaNWsQFxeHVatWYcOGDcjJyUFAQMA1+//973/Hl19+if/85z8YMmQIfv75ZyxduhT79+/HsGHDbvh+7LnpnsqqDUg9oUVyuga/FZTh6j+VI4N9kRgjx6xoOQK8XG1XJBER2Yxd3ZaKi4vDqFGj8N577wEATCYTlEolFi9ejKeffvqa/RUKBZ555hksWrTI3Hb77bfDzc0NX3755Q3fj+Gm+yuuqsNPGVqkqDQ4fLbcHHREImB0iB8SYxWYGRUIfw+pbQslIqIuY8n3t1MX1dSi+vp6HD16FMuWLTO3icViTJ06FQcOHGjxGIPBAFfX5r+9u7m5Ye/evdfd32C4shCkXq/vhMrJmgI8XZE0NgRJY0OgraxDSoYGKSo1jhVW4GBBOQ4WlOP5HzIRP7A3EmMUmBEZCN9eLrYum4iIugmbhpvS0lIYjUbIZLJm7TKZDNnZ2S0eM336dLz55pu4+eabMXDgQGzfvh3ff/89jEZji/uvXLkSK1as6PTaqWsEervi/pv64/6b+uPchVpsydAgRaVB+rlK7Msrw768Mjy7KRPjBvkjMUaO6RGB8HZ3tnXZRERkQza9LaVWqxEUFIT9+/cjPj7e3P7UU09h9+7dOHjw4DXHlJSU4MEHH8TmzZshEokwcOBATJ06FZ9++ikuXrx28GlLPTdKpZK3pexcYVktkjPUSFFpcEJ9pTfOWSLC+MF9kBgjxy0RMni6MugQETkCu7kt5e/vD4lEAp1O16xdp9MhMDCwxWP69OmDTZs2oa6uDmVlZVAoFHj66acxYMCAFveXSqWQSjk2w9H06+2Ov00chL9NHIT8kmqkqDRIVmmQo6vCjuxi7MguhouTGBNCm4LO1HAZeklt+sediIi6iE3/tXdxccGIESOwfft2zJkzB0DTgOLt27fj4YcfbvVYV1dXBAUFoaGhAd999x3uvPPOLqiYuqMBfTyweMpgLJ4yGKd0VUhWaZCsUuN0SQ22ntRh60kdXJ3FmDwkAAnRCkweEgA3F4mtyyYiIiux+dNSX3/9NZKSkvDhhx9i9OjRWLVqFb755htkZ2dDJpNhwYIFCAoKwsqVKwEABw8exPnz5zF06FCcP38e//znP1FQUIBjx47Bx8fnhu/Hp6V6BkEQkK2tutSjo8aZslrzNncXCaaEy5AYI8eE0D5wdWbQISLq7uzmthQAzJs3DyUlJVi+fDm0Wi2GDh2K1NRU8yDjwsJCiMVX5hqsq6vDs88+i/z8fHh4eGDWrFn44osv2hRsqOcQiUQIl3shXO6Fx6eF4oRab+7ROXfhIjanq7E5XQ0PqRNuiWgKOjcN9ofUiUGHiMjedbjnRq/XY8eOHQgLC0N4eHhn1WU17Lnp2QRBQPq5SqSomgYjqyvrzNs8XZ0wPTIQiTFyjBvkD2dJt5jAm4iIYOVJ/O68807cfPPNePjhh3Hx4kXExsbizJkzEAQB69evx+23396h4q2N4YYuM5kEHC+6gM3pGmzJ0KC46spTdT7uzpgRGYjEGAXGDPCDE4MOEZFNWTXcBAYG4ueff0ZsbCy++uorPP/880hPT8dnn32Gjz76CMePH+9Q8dbGcEMtMZkEHD5TjmSVBj9lalBaXW/e1ruXC2ZENQWd0f39IBGLbFgpEVHPZNVw4+bmhtzcXCiVSixYsAAKhQKvvPIKCgsLERERgerq6g4Vb20MN3QjjUYTDhWUY7NKg9RMDS7UNpi39fGUIiFajoQYOUb084WYQYeIqEtYdUCxUqnEgQMH4Ofnh9TUVKxfvx4AcOHChWuWRSCyR04SMcYO8sfYQf544bZIHDhdhmSVGqmZWpRUGbBu/xms238Gcm9XzLoUdIYpfSASMegQEXUHFvfcvP/++3j00Ufh4eGB4OBgHDt2DGKxGO+++y6+//577Ny501q1dgr23FB71TeasC+vFJtVamw9oUOVodG8LcjHDYkxTUEnOsibQYeIqJNZfVXwI0eOoKioCLfccgs8PDwAACkpKfDx8cG4cePaV3UXYbihzlDXYMSe3BKkZGiw7aQONfVX1jYL7u1uvnUVIfdi0CEi6gRWDzdXMxqNyMjIQHBwMHx9fTtyqi7BcEOdra7BiJ3ZxUjO0GB7lg51DSbztgH+vZAYI0dirAKhMk8bVklEZN+sGm4ee+wxREdH4/7774fRaMSECROwf/9+uLu7Izk5GRMnTuxI7VbHcEPWVFvfiO1ZxUhWqbEzpwT1jVeCTqjMAwnRCiTGyjGwj4cNqyQisj9WDTd9+/bFpk2bMHLkSGzatAmLFi3Czp078cUXX2DHjh3Yt29fh4q3NoYb6irVhkZsO6lDskqNPbmlqDdeCTrhcq+mHp0YOYJ797JhlURE9sGq4cbV1RV5eXno27cv/vznP8Pd3R2rVq1CQUEBYmNjodfrO1S8tTHckC1UXmzA1ktBZ++pUjSarvy1iw7yRkKMHAnRcij93G1YJRFR92XVR8FlMhlOnjwJuVyO1NRUfPDBBwCA2tpaSCRcl4eoJd5uzrhjRF/cMaIvKmrr8fMJLZJVGuw/XYaM85XIOF+JV37KxlClDxJj5JgVLYfCx83WZRMR2SWLw819992HO++8E3K5HCKRCFOnTgXQtFr3kCFDOr1AIkfj4+6CeaP6Yd6ofiirNiD1hBbJ6Rr8VlCGtKIKpBVV4KWULIwM9jX36AR4cQ4pIqK2atfTUt9++y2Kiorwxz/+EX379gUAfPbZZ/Dx8cFtt93W6UV2Jt6Wou6quKoOP2VokaLS4PDZclz+mykSAaND/JAYq8DMqED4e0htWygRkQ106aPg9obhhuyBtrIOKRkapKjUOFZYYW4Xi4D4gb2RGKPAjMhA+PZysV2RRERdyOrhpqamBrt370ZhYSHq6+ubbXvkkUcsPV2XYrghe3PuQi22ZGiQotIg/VyluV0iFmHcIH8kxsgxPSIQ3u7ONqySiMi6rBpujh8/jlmzZqG2thY1NTXw8/NDaWkp3N3dERAQgPz8/A4Vb20MN2TPCstqkZyhRopKgxPqK08mOktEGD+4DxJj5LglQgZPVwYdInIsVg03EydORGhoKNasWQNvb2+kp6fD2dkZf/rTn/Doo49i7ty5HSre2hhuyFHkl1QjRaVBskqDHF2Vud3FSYwJoU1BZ2q4DL2kFj83QETU7Vg13Pj4+ODgwYMICwuDj48PDhw4gPDwcBw8eBBJSUnIzs7uUPHWxnBDjuiUrgrJKg2SVWqcLqkxt7s6izF5SAASohWYPCQAbi6croGI7JNV57lxdnaGWCwGAAQEBKCwsBDh4eHw9vZGUVFR+yomog4ZLPPEkls88djUwcjWVl3q0VHjTFkttmRosSVDC3cXCaaEy5AYI8eE0D5wdWbQISLHZHG4GTZsGA4fPozBgwdjwoQJWL58OUpLS/HFF18gKirKGjUSURuJRCKEy70QLvfC49NCcUKtN/fonLtwEZvT1dicroaH1Am3RDQFnZsG+0PqxKBDRI5DbOkBL7/8MuRyOQDgX//6F3x9ffHXv/4VJSUl+Oijj9pVxOrVqxESEgJXV1fExcXh0KFDre6/atUqhIWFwc3NDUqlEkuWLEFdXV273pvIUYlEIkQFeePpmUPw61OTsGnRODw4vj8U3q6oNjRi4/HzuP+zIxj50jY8sSEdu3KK0XDV+ldERPbK5vPcfP3111iwYAHWrFmDuLg4rFq1Chs2bEBOTg4CAgKu2f+rr77CwoUL8emnn2Ls2LHIzc3Fvffei//7v//Dm2++ecP345gb6ulMJgHHiy5gc7oGWzI0KK4ymLf5uDtjRmQgEmMUGDPAD04Si3//ISKyCruaxC8uLg6jRo3Ce++9BwAwmUxQKpVYvHgxnn766Wv2f/jhh5GVlYXt27eb2x5//HEcPHgQe/fuvWZ/g8EAg+HKP956vR5KpZLhhghNQefwmXIkqzT4KVOD0uor81b17uWCGVFNQWd0fz9IxCIbVkpEPZ0l4cbiX8vKysqwaNEiREREwN/fH35+fs1+LFFfX4+jR4+a16cCALFYjKlTp+LAgQMtHjN27FgcPXrUfOsqPz8fW7ZswaxZs1rcf+XKlfD29jb/KJVKi2okcmRisQhxA3rjxTlR+G3ZFHz1QBzuGt0Pvu7OKKupx38PFuKu//yGMSu34/kfMnH4TDlMph41qTkR2SGLe25mzZqFvLw83H///ZDJZBCJmv82l5SU1OZzqdVqBAUFYf/+/YiPjze3P/XUU9i9ezcOHjzY4nHvvPMOnnjiCQiCgMbGRjz00EPm1cl/jz03RJZrMJpw4HQZklVqpGZqoa9rNG+Te7tiVrQcCTFyDFP6XPNvABGRNVj1UfBff/0Ve/fuRWxsbLsL7Ihdu3bh5Zdfxvvvv4+4uDjk5eXh0UcfxYsvvojnnnvumv2lUimkUi40SGQJZ4kYN4f2wc2hffDSnGjsyyvFZpUaW0/ooKmswyd7C/DJ3gIE+bghMaYp6EQHeTPoEFG3YHG4GTJkCC5evNgpb+7v7w+JRAKdTtesXafTITAwsMVjnnvuOdxzzz144IEHAADR0dGoqanBn//8ZzzzzDPmOXiIqHO4OIkxaUgAJg0JQF2DEXtyS5CSocG2kzqcr7iID/fk48M9+Qju7Y6ESz06EXIvBh0ishmLk8D777+PZ555Brt370ZZWRn0en2zH0u4uLhgxIgRzQYHm0wmbN++vdltqqvV1tZeE2AkkqY5OnrYAudEXc7VWYJpkYF4+/+G4ehzt+CD+cORECOHq7MYZ8tq8f6u00h4Zy+mvLEbb/6Sgxxt1Y1PSkTUySzuufHx8YFer8fkyZObtQuCAJFIBKPRaNH5li5diqSkJIwcORKjR4/GqlWrUFNTg/vuuw8AsGDBAgQFBWHlypUAgNmzZ+PNN9/EsGHDzLelnnvuOcyePdsccojI+lydJZgZLcfMaDlq6xuxPasYySo1duaUIL+0Bu/syMM7O/IwOMADiTEKJMbKMbCPh63LJqIewOJwM3/+fDg7O+Orr75qcUCxpebNm4eSkhIsX74cWq0WQ4cORWpqKmQyGQCgsLCwWU/Ns88+C5FIhGeffRbnz59Hnz59MHv2bPzrX//qUB1E1H7uLk6YHavA7FgFqg2N2HZSh2SVGntyS3GquBpvbcvFW9tyES73QmKMHIkxcgT37mXrsonIQVn8tJS7uzuOHz+OsLAwa9VkVZzEj6jrVF5swNZLQWfvqVI0XvUYeXSQNxJi5EiIlkPp527DKonIHlj1aamRI0eiqKjIbsMNEXUdbzdn3DGiL+4Y0RcVtfX4+YQWySoN9p8uQ8b5SmScr8QrP2VjqNIHiTFyzIqWQ+HjZuuyicjOWdxzs2HDBvzzn//Ek08+iejoaDg7OzfbHhMT06kFdjb23BDZXlm1AakntEhO1+C3gjJc/a/QyGBfc49OgJer7Yokom7FqssvtPSotUgkaveA4q7GcEPUvRRX1eGnDC1SVBocPltuDjoiETA6xA+JsQrMjAqEvwfnqyLqyawabs6ePdvq9uDgYEtO1+UYboi6L21lHVIyNEhRqXGssMLcLhYB8QN7IzFGgRmRgfDt5WK7IonIJuxq4cyuxnBDZB/OXajFlgwNUlQapJ+rNLdLxCKMG+SPxBg5pkcEwtvduZWzEJGjYLhpBcMNkf0pLKtFcoYaKSoNTqivTBbqLBFh/OA+SIyR45YIGTxdGXSIHBXDTSsYbojsW35JNVJUGiSrNMjRXZkB2cVJjAmhTUFnargMvaQWPwxKRN0Yw00rGG6IHMcpXRWSVRokq9Q4XVJjbpc6iTF5SAASYxSYPCQAbi6cvZzI3lkt3BiNRuzbtw8xMTHw8fHpaJ02wXBD5HgEQUC2tupSj44aZ8pqzdvcXSSYEi5DQrQcE8P6wNWZQYfIHlm158bV1RVZWVno379/h4q0FYYbIscmCAJOqPXmHp1zFy6at3lInXBLhAyJMXLcNNgfUicGHSJ7YdUZiqOiopCfn2+34YaIHJtIJEJUkDeigrzx9xlhSD9XiRRV02BkdWUdNh4/j43Hz8PT1QnTIwORGCPHuEH+cJZcO4cXEdkni3tuUlNTsWzZMrz44osYMWIEevVqvvhdd+8NYc8NUc9kMgk4XnQBm9M12JKhQXGVwbzNx90ZMyIDkRijwJgBfnBi0CHqdrpshuKrVwTnDMVEZC9MJgGHz5QjWaXBT5kalFbXm7f17uWCGVFNQWd0fz9IxKJWzkREXcWq4Wb37t2tbp8wYYIlp+tyDDdEdLVGowmHCsqxWaVBaqYGF2obzNv6eEoxKyoQibEKjOjnCzGDDpHN8FHwVjDcENH1NBhNOHC6DMkqNVIztdDXNZq3yb1dMStajoQYOYYpfZr1XBOR9XVJuKmtrUVhYSHq6+ubtXNVcCJyBPWNJuzLK8VmlRpbT+hQZbgSdIJ83JAY0xR0ooO8GXSIuoBVw01JSQnuu+8+/PTTTy1u55gbInI0dQ1G7MktQUqGBttO6lBTf+XfueDe7ki41KMTIfdi0CGyEquGm/nz5+Ps2bNYtWoVJk6ciI0bN0Kn0+Gll17CG2+8gYSEhA4Vb20MN0TUEXUNRuzMLkZyhgbbs3SoazCZtw3w73WpR0eBsEBPG1ZJ5HisGm7kcjl++OEHjB49Gl5eXjhy5AhCQ0Px448/4tVXX8XevXs7VLy1MdwQUWeprW/E9qxiJKvU2JlTgvrGK0FncIAHEmMUSIyVY2AfDxtWSeQYLPn+tngyh5qaGgQEBAAAfH19UVJSAgCIjo7GsWPH2lEusHr1aoSEhMDV1RVxcXE4dOjQdfedOHEiRCLRNT/dvceIiByPu4sTZscq8OE9I3HsuVuwat5QTA0PgItEjFPF1XhrWy6mvLEbM1btweqdeThbVnPjkxJRh1k8Q3FYWBhycnIQEhKC2NhYfPjhhwgJCcGaNWsgl8stLuDrr7/G0qVLsWbNGsTFxWHVqlWYPn06cnJyzCHqat9//32zQcxlZWWIjY3FH//4R4vfm4ios3hInTBnWBDmDAtC5cUGbD2pQ7JKjb2nSpGtrUK2Ngev/ZyD6CBvJMTIkRAth9LP3dZlEzkki29Lffnll2hsbMS9996Lo0ePYsaMGSgvL4eLiwvWrVuHefPmWVRAXFwcRo0ahffeew8AYDKZoFQqsXjxYjz99NM3PH7VqlVYvnw5NBrNNbMlt4S3pYioK1XU1uPnE1okqzTYf7oMRtOVf3KHKn2QGCPHrGg5FD5uNqySqPvr0nluamtrkZ2djX79+sHf39+iY+vr6+Hu7o5vv/0Wc+bMMbcnJSWhoqICP/zwww3PER0djfj4eHz00UctbjcYDDAYrkyzrtfroVQqGW6IqMuVVRuQekKL5HQNfisow9X/+o4M9jX36AR4udquSKJuyqoLZ/6eu7s7hg8f3q5jS0tLYTQaIZPJmrXLZDJkZ2ff8PhDhw4hMzMTn3zyyXX3WblyJVasWNGu+oiIOlNvDynmxwVjflwwiqvq8FOGFikqDQ6fLceRsxdw5OwFvJB8EqND/JAYq8DMqED4e0htXTaR3WlTz83SpUvbfMI333yzzfuq1WoEBQVh//79iI+PN7c/9dRT2L17Nw4ePNjq8X/5y19w4MABqFSq6+7Dnhsi6u60lXVIydAgRaXGscIKc7tYBMQP7I3EGAVmRAbCt5eL7YoksrFO77k5fvx4m97Y0smr/P39IZFIoNPpmrXrdDoEBga2emxNTQ3Wr1+PF154odX9pFIppFL+5kNE3Vegtyvuv6k/7r+pP85dqMWWDA1SVBqkn6vEvrwy7Msrw7ObMjFukD8SY+SYHhEIb3dnW5dN1G3ZfG2puLg4jB49Gu+++y6ApgHF/fr1w8MPP9zqgOJ169bhoYcewvnz59G7d+82vx8HFBORvSgsq0VyhhopKg1OqPXmdmeJCOMH90FijBy3RMjg6cqgQ47PrhbO/Prrr5GUlIQPP/wQo0ePxqpVq/DNN98gOzsbMpkMCxYsQFBQEFauXNnsuPHjxyMoKAjr16+36P0YbojIHuWXVCNFpUGySoMcXZW53cVJjAmhTUFnargMvaQdHkpJ1C11+m2puXPntvnNv//++zbvCwDz5s1DSUkJli9fDq1Wi6FDhyI1NdU8yLiwsBBicfO5BnNycrB371788ssvFr0XEZG9GtDHA4unDMbiKYNxSleFZJUGySo1TpfUYOtJHbae1EHqJMbkIQFIjFFg8pAAuLlIbF02kU20qefmvvvua/MJ165d26GCrI09N0TkKARBQLa26lKPjhpnymrN29xdJJgSLkNCtBwTw/rA1ZlBh+ybXd2W6moMN0TkiARBwAm13tyjc+7CRfM2D6kTboloCjrjQ/0hdWLQIftj9XDT2NiIXbt24fTp07j77rvh6ekJtVoNLy8veHh07wXiGG6IyNEJgoD0c5VIUTUNRlZX1pm3ebo6YXpkIBJj5Bg3yB/OEouXGCSyCauGm7Nnz2LGjBkoLCyEwWBAbm4uBgwYgEcffRQGgwFr1qzpUPHWxnBDRD2JySTgeNEFbE7XYEuGBsVVV+b98nF3xozIQCTGKDBmgB+cGHSoG7NquJkzZw48PT3xySefoHfv3khPT8eAAQOwa9cuPPjggzh16lSHirc2hhsi6qlMJgGHz5QjWaXBT5kalFZfWYS4dy8XzIhqCjqj+/tBIrZs3jIia7NquOnduzf279+PsLAweHp6msPNmTNnEBERgdra2hufxIYYboiIgEajCYcKyrFZpUFqpgYXahvM2/p4SjErKhCJsQqM6OcLMYMOdQNWXVvKZDLBaDRe037u3Dl4enpaejoiIrIBJ4kYYwf5Y+wgf7xwWyQOnC5DskqN1EwtSqoM+OzAWXx24CwCvVwxK1qOxFg5hil9LJ6JnsgWLO65mTdvHry9vfHRRx/B09MTKpUKffr0wW233YZ+/frxUXAiIjtW32jCvrxSbFapsfWEDlWGRvO2IB83JMbIkRAjR3SQN4MOdSmr3pY6d+4cpk+fDkEQcOrUKYwcORKnTp2Cv78/9uzZg4CAgA4Vb20MN0REbVPXYMSe3BKkZGiw7aQONfVXeu2De7sjIbop6ETIvRh0yOq65FHwr7/+Gunp6aiursbw4cMxf/58uLm5tbvorsJwQ0RkuboGI3ZmFyM5Q4PtWTrUNZjM2wb497rUo6NAWCCHJ5B1WDXc7NmzB2PHjoWTU/PhOo2Njdi/fz9uvvlmyyvuQgw3REQdU1vfiO1ZxUhWqbEzpwT1jVeCzuAADyTGKJAYK8fAPt173jOyL1YNNxKJBBqN5prbT2VlZQgICGhxsHF3wnBDRNR5qg2N2HZSh2SVGntyS1FvvBJ0hgR6YnasAgnRcoT497JhleQIrBpuxGIxdDod+vTp06w9NzcXI0eOhF6vt7ziLsRwQ0RkHZUXG7D1UtDZe6oUjaYrXy9RQV5IjGkKOko/dxtWSfbKKuHm8srgP/zwA2bMmAGpVGreZjQaoVKpEBYWhtTU1A6Ubn0MN0RE1ldRW4+fT2iRrNJg/+kyGK8KOkOVPkiMkWNWtBwKn+4/VpO6B6vMc+Pt7Q2gac0ST0/PZoOHXVxcMGbMGDz44IPtLJmIiByJj7sL5o3qh3mj+qGs2oDUE1okp2vwW0EZ0ooqkFZUgZdSsjAy2BcJMXIkRMsR4OVq67LJQVh8W2rFihV44okn0KuXfd4/Zc8NEZHtFFfV4acMLVJUGhw+W47L30AiETA6xA+JsQrMjAqEv4e09RNRj2PVMTcXL16EIAhwd2+6Z3r27Fls3LgRERERmDZtWvur7iIMN0RE3YO2sg4pGRqkqNQ4VlhhbheLgPiBvZEYo8D0yED49XKxXZHUbVg13EybNg1z587FQw89hIqKCoSFhcHFxQWlpaV488038de//rVDxVsbww0RUfdz7kIttmRokKLSIP1cpbldIhZh3CB/JMbIMT0iEN7uzjaskmzJquHG398fu3fvRmRkJD7++GO8++67OH78OL777jssX74cWVlZHSre2hhuiIi6t8KyWiRnqJGi0uCE+soTuM4SEcYP7oPEGDluiZDB05VBpyex6sKZtbW15gUyf/nlF8ydOxdisRhjxozB2bNn21cxERHRJf16u+NvEwfhbxMHIb+kGikqDZJVGuToqrAjuxg7sovh4iTGhNCmoDM1XIZeUou/zsiBWdxzExMTgwceeAB/+MMfEBUVhdTUVMTHx+Po0aNISEiAVqu1Vq2dgj03RET26ZSuCskqDZJVapwuqTG3S53EmDwkAIkxCkweEgA3F4kNqyRrseT7W2zpyZcvX44nnngCISEhiIuLQ3x8PICmXpxhw4a1q+DVq1cjJCQErq6uiIuLw6FDh1rdv6KiAosWLYJcLodUKkVoaCi2bNnSrvcmIiL7MFjmiSW3hGLb0gn46dHxeHjSIIT0doeh0YSfMrVY9NUxjHhpKxb/7zhSM7Woa+jeM+aT9bRr4UytVguNRoPY2FiIxU356NChQ/Dy8sKQIUMsOtfXX3+NBQsWYM2aNYiLi8OqVauwYcMG5OTktLjCeH19PcaNG4eAgAD84x//QFBQEM6ePQsfHx/Exsbe8P3Yc0NE5DgEQcAJtd7co3PuwkXzNg+pE26JkCEhWo7xof6QOrFHx55ZfVXwzhQXF4dRo0bhvffeAwCYTCYolUosXrwYTz/99DX7r1mzBq+99hqys7Ph7Gz5YDKGGyIixyQIAtLPVSJF1TQYWV1ZZ97m6eqE6ZGBSIyRY9wgfzhLLL5xQTZmN+Gmvr4e7u7u+PbbbzFnzhxze1JSEioqKvDDDz9cc8ysWbPg5+cHd3d3/PDDD+jTpw/uvvtu/P3vf4dEcm0qNxgMMBgM5td6vR5KpZLhhojIgZlMAo4XXcDmdA22ZGhQXHXle8DH3RkzIgORGKPAmAF+cGLQsQtWfVqqM5WWlsJoNEImkzVrl8lkyM7ObvGY/Px87NixA/Pnz8eWLVuQl5eHv/3tb2hoaMDzzz9/zf4rV67EihUrrFI/ERF1T2KxCCOC/TAi2A/LEyNw+Ew5klUa/JSpQWl1PdYfLsL6w0Xo3csFM6Kags7o/n6QiEW2Lp06gU17btRqNYKCgrB//37zwGQAeOqpp7B7924cPHjwmmNCQ0NRV1eHgoICc0/Nm2++iddeew0ajeaa/dlzQ0RElzUaTThUUI7NKg1SMzW4UNtg3tbHU4pZUYFIjFVgRD9fiBl0uhW76bnx9/eHRCKBTqdr1q7T6RAYGNjiMXK5HM7Ozs1uQYWHh0Or1aK+vh4uLs2n6ZZKpc1WMCciop7LSSLG2EH+GDvIHy/cFokDp8uQrFIjNVOLkioDPjtwFp8dOItAL1fMipYjMVaOYUofiEQMOvbEpjcaXVxcMGLECGzfvt3cZjKZsH379mY9OVcbN24c8vLyYDKZzG25ubmQy+XXBBsiIqLrcZaIcXNoH7x6RyyOPHsL1t47CnOHB8FT6gStvg6f7ivA3Pf346Z/78TLW7KgOlcBGz+DQ21k86elvv76ayQlJeHDDz/E6NGjsWrVKnzzzTfIzs6GTCbDggULEBQUhJUrVwIAioqKEBkZiaSkJCxevBinTp3CwoUL8cgjj+CZZ5654fvxaSkiImpNXYMRe3JLkJKhwbaTOtTUX5kvJ7i3OxKi5UiIkSNC7sUenS5kN7elAGDevHkoKSnB8uXLodVqMXToUKSmppoHGRcWFprn0gEApVKJn3/+GUuWLEFMTAyCgoLw6KOP4u9//7utPgIRETkQV2cJpkUGYlpkIOoajNiZXYzkDA22Z+lwtqwW7+86jfd3ncYA/15IjJEjIUaBsEBPW5dNV7F5z01XY88NERG1R219I7ZnFSNZpcbOnBLUN14ZHjE4wAOJMQokxsoxsI+HDat0XHYzz40tMNwQEVFHVRsase2kDskqNfbklqLeeCXoDAn0xOxYBRKi5Qjx72XDKh0Lw00rGG6IiKgzVV5swNZLQWfvqVI0mq58rUYFeSExpinoKP3cbVil/WO4aQXDDRERWUtFbT1+PqFFskqD/afLYLwq6AxV+iAxRo5Z0XIofNxsWKV9YrhpBcMNERF1hbJqA1JPaJGcrsFvBWW4+tt2ZLAvEmLkSIiWI8DL1XZF2hGGm1Yw3BARUVcrrqrDTxlapKg0OHy23Bx0RCJgdIgfEmMVmBkVCH8PTjp7PQw3rWC4ISIiW9JW1iElQ4MUlRrHCivM7WIRED+wNxJjFJgeGQi/XpyY9moMN61guCEiou7i3IVabMnQIEWlQfq5SnO7RCzCuEH+SIyRY3pEILzdnW1YZffAcNMKhhsiIuqOCstqkZyhRopKgxNqvbndWSLC+MF9kBgjxy0RMni69sygw3DTCoYbIiLq7vJLqpGi0iBZpUGOrsrc7uIkxoTQpqAzNVyGXlKbLzTQZRhuWsFwQ0RE9uSUrgrJKg2SVWqcLqkxt0udxJg8JACJMQpMHhIANxeJDau0PoabVjDcEBGRPRIEAdnaqks9OmqcKas1b3N3kWBKuAwJ0XJMDOsDV2fHCzoMN61guCEiInsnCAJOqPXmHp1zFy6at3lInXBLRFPQGR/qD6mTYwQdhptWMNwQEZEjEQQB6ecqkaJqGoysrqwzb/N0dcL0yEAkxsgxbpA/nCViG1baMQw3rWC4ISIiR2UyCThedAGb0zXYkqFBcZXBvM3H3RkzIgORGKPAmAF+cLKzoMNw0wqGGyIi6glMJgGHz5QjWaXBT5kalFbXm7f17uWCGVFNQWd0fz9IxCIbVto2DDetYLghIqKeptFowqGCcmxWaZCaqcGF2gbztj6eUsyKCkRirAIj+vlC3E2DDsNNKxhuiIioJ2swmnDgdBmSVWqkZmqhr2s0bwv0csWsaDkSY+UYpvSBSNR9gg7DTSsYboiIiJrUN5qwL68Um1VqbD2hQ5XhStAJ8nFDQowciTFyRAd52zzoMNy0guGGiIjoWnUNRuzJLUFKhgbbTupQU280b+vn547EGDkSYuSIkHvZJOhY8v3dLYZKr169GiEhIXB1dUVcXBwOHTp03X3XrVsHkUjU7MfV1bULqyUiInI8rs4STIsMxNv/NwxHn7sFH8wfjoQYOVydxSgsr8X7u04j4Z29mPLGbrz5Sw5ytFU3PqmN2HxRiq+//hpLly7FmjVrEBcXh1WrVmH69OnIyclBQEBAi8d4eXkhJyfH/NrWXWVERESOxNVZgpnRcsyMlqO2vhHbs4qRrFJjZ04J8ktr8M6OPLyzIw+DAzyQGKNAYqwcA/t42LpsM5vfloqLi8OoUaPw3nvvAQBMJhOUSiUWL16Mp59++pr9161bh8ceewwVFRXtej/eliIiImqfakMjtp3UIVmlxp7cUtQbTeZtQwI9MTtWgYRoOUL8e3X6e1vy/W3Tnpv6+nocPXoUy5YtM7eJxWJMnToVBw4cuO5x1dXVCA4OhslkwvDhw/Hyyy8jMjKyxX0NBgMMhiuTGOn1+hb3IyIiotZ5SJ0wZ1gQ5gwLQuXFBmy9FHT2nipFtrYK2docvPZzDqKCvPC/B8fA09XZJnXadMxNaWkpjEYjZDJZs3aZTAatVtviMWFhYfj000/xww8/4Msvv4TJZMLYsWNx7ty5FvdfuXIlvL29zT9KpbLTPwcREVFP4+3mjDtG9MW6+0bjyLNT8e/bozF+sD8kYhEajYLNgg1g49tSarUaQUFB2L9/P+Lj483tTz31FHbv3o2DBw/e8BwNDQ0IDw/HXXfdhRdffPGa7S313CiVSt6WIiIisoKyagM0lXWICvLu1PPazW0pf39/SCQS6HS6Zu06nQ6BgYFtOoezszOGDRuGvLy8FrdLpVJIpdIO10pEREQ31ttDit4etv3eteltKRcXF4wYMQLbt283t5lMJmzfvr1ZT05rjEYjMjIyIJfLrVUmERER2RGbPwq+dOlSJCUlYeTIkRg9ejRWrVqFmpoa3HfffQCABQsWICgoCCtXrgQAvPDCCxgzZgwGDRqEiooKvPbaazh79iweeOABW34MIiIi6iZsHm7mzZuHkpISLF++HFqtFkOHDkVqaqp5kHFhYSHE4isdTBcuXMCDDz4IrVYLX19fjBgxAvv370dERIStPgIRERF1Izaf56arcZ4bIiIi+2M3A4pt4XKW43w3RERE9uPy93Zb+mR6XLipqmpaC4Pz3RAREdmfqqoqeHu3/ph5j7stZTKZoFar4enp2elrUl2eQ6eoqIi3vKyI17lr8Dp3DV7nrsNr3TWsdZ0FQUBVVRUUCkWzsbgt6XE9N2KxGH379rXqe3h5efEvThfgde4avM5dg9e56/Badw1rXOcb9dhcZtN5boiIiIg6G8MNERERORSGm04klUrx/PPPc7kHK+N17hq8zl2D17nr8Fp3je5wnXvcgGIiIiJybOy5ISIiIofCcENEREQOheGGiIiIHArDDRERETkUhhsLrV69GiEhIXB1dUVcXBwOHTrU6v4bNmzAkCFD4OrqiujoaGzZsqWLKrVvllzn//znPxg/fjx8fX3h6+uLqVOn3vC/CzWx9M/zZevXr4dIJMKcOXOsW6CDsPQ6V1RUYNGiRZDL5ZBKpQgNDeW/HW1k6bVetWoVwsLC4ObmBqVSiSVLlqCurq6LqrU/e/bswezZs6FQKCASibBp06YbHrNr1y4MHz4cUqkUgwYNwrp166xeJwRqs/Xr1wsuLi7Cp59+Kpw4cUJ48MEHBR8fH0Gn07W4/759+wSJRCK8+uqrwsmTJ4Vnn31WcHZ2FjIyMrq4cvti6XW+++67hdWrVwvHjx8XsrKyhHvvvVfw9vYWzp0718WV2xdLr/NlBQUFQlBQkDB+/Hjhtttu65pi7Zil19lgMAgjR44UZs2aJezdu1coKCgQdu3aJaSlpXVx5fbH0mv93//+V5BKpcJ///tfoaCgQPj5558FuVwuLFmypIsrtx9btmwRnnnmGeH7778XAAgbN25sdf/8/HzB3d1dWLp0qXDy5Enh3XffFSQSiZCammrVOhluLDB69Ghh0aJF5tdGo1FQKBTCypUrW9z/zjvvFBISEpq1xcXFCX/5y1+sWqe9s/Q6/15jY6Pg6ekpfPbZZ9Yq0SG05zo3NjYKY8eOFT7++GMhKSmJ4aYNLL3OH3zwgTBgwAChvr6+q0p0GJZe60WLFgmTJ09u1rZ06VJh3LhxVq3TUbQl3Dz11FNCZGRks7Z58+YJ06dPt2JlgsDbUm1UX1+Po0ePYurUqeY2sViMqVOn4sCBAy0ec+DAgWb7A8D06dOvuz+17zr/Xm1tLRoaGuDn52etMu1ee6/zCy+8gICAANx///1dUabda891/vHHHxEfH49FixZBJpMhKioKL7/8MoxGY1eVbZfac63Hjh2Lo0ePmm9d5efnY8uWLZg1a1aX1NwT2Op7sMctnNlepaWlMBqNkMlkzdplMhmys7NbPEar1ba4v1artVqd9q491/n3/v73v0OhUFzzF4quaM913rt3Lz755BOkpaV1QYWOoT3XOT8/Hzt27MD8+fOxZcsW5OXl4W9/+xsaGhrw/PPPd0XZdqk91/ruu+9GaWkpbrrpJgiCgMbGRjz00EP4xz/+0RUl9wjX+x7U6/W4ePEi3NzcrPK+7Lkhh/LKK69g/fr12LhxI1xdXW1djsOoqqrCPffcg//85z/w9/e3dTkOzWQyISAgAB999BFGjBiBefPm4ZlnnsGaNWtsXZrD2bVrF15++WW8//77OHbsGL7//nukpKTgxRdftHVp1EHsuWkjf39/SCQS6HS6Zu06nQ6BgYEtHhMYGGjR/tS+63zZ66+/jldeeQXbtm1DTEyMNcu0e5Ze59OnT+PMmTOYPXu2uc1kMgEAnJyckJOTg4EDB1q3aDvUnj/Pcrkczs7OkEgk5rbw8HBotVrU19fDxcXFqjXbq/Zc6+eeew733HMPHnjgAQBAdHQ0ampq8Oc//xnPPPMMxGL+/t9R1/se9PLyslqvDcCemzZzcXHBiBEjsH37dnObyWTC9u3bER8f3+Ix8fHxzfYHgK1bt153f2rfdQaAV199FS+++CJSU1MxcuTIrijVrll6nYcMGYKMjAykpaWZf2699VZMmjQJaWlpUCqVXVm+3WjPn+dx48YhLy/PHB4BIDc3F3K5nMGmFe251rW1tdcEmMuhUuCyi53CZt+DVh2u7GDWr18vSKVSYd26dcLJkyeFP//5z4KPj4+g1WoFQRCEe+65R3j66afN++/bt09wcnISXn/9dSErK0t4/vnn+Sh4G1h6nV955RXBxcVF+PbbbwWNRmP+qaqqstVHsAuWXuff49NSbWPpdS4sLBQ8PT2Fhx9+WMjJyRGSk5OFgIAA4aWXXrLVR7Abll7r559/XvD09BT+97//Cfn5+cIvv/wiDBw4ULjzzjtt9RG6vaqqKuH48ePC8ePHBQDCm2++KRw/flw4e/asIAiC8PTTTwv33HOPef/Lj4I/+eSTQlZWlrB69Wo+Ct4dvfvuu0K/fv0EFxcXYfTo0cJvv/1m3jZhwgQhKSmp2f7ffPONEBoaKri4uAiRkZFCSkpKF1dsnyy5zsHBwQKAa36ef/75ri/czlj65/lqDDdtZ+l13r9/vxAXFydIpVJhwIABwr/+9S+hsbGxi6u2T5Zc64aGBuGf//ynMHDgQMHV1VVQKpXC3/72N+HChQtdX7id2LlzZ4v/3l6+rklJScKECROuOWbo0KGCi4uLMGDAAGHt2rVWr1MkCOx7IyIiIsfBMTdERETkUBhuiIiIyKEw3BAREZFDYbghIiIih8JwQ0RERA6F4YaIiIgcCsMNERERORSGGyIiIuoUe/bswezZs6FQKCASibBp0yaLzyEIAl5//XWEhoZCKpUiKCgI//rXvyw6B8MNEdmFXbt2QSQSoaKiwtalENF11NTUIDY2FqtXr273OR599FF8/PHHeP3115GdnY0ff/wRo0ePtugcnKGYiLqliRMnYujQoVi1ahUAoL6+HuXl5ZDJZBCJRLYtjohuSCQSYePGjZgzZ465zWAw4JlnnsH//vc/VFRUICoqCv/+978xceJEAEBWVhZiYmKQmZmJsLCwdr83e26IyC64uLggMDCQwYbIjj388MM4cOAA1q9fD5VKhT/+8Y+YMWMGTp06BQDYvHkzBgwYgOTkZPTv3x8hISF44IEHUF5ebtH7MNwQUbdz7733Yvfu3Xj77bchEokgEomwbt26Zrel1q1bBx8fHyQnJyMsLAzu7u644447UFtbi88++wwhISHw9fXFI488AqPRaD63wWDAE088gaCgIPTq1QtxcXHYtWuXbT4oUQ9SWFiItWvXYsOGDRg/fjwGDhyIJ554AjfddBPWrl0LAMjPz8fZs2exYcMGfP7551i3bh2OHj2KO+64w6L3crLGByAi6oi3334bubm5iIqKwgsvvAAAOHHixDX71dbW4p133sH69etRVVWFuXPn4g9/+AN8fHywZcsW5Ofn4/bbb8e4ceMwb948AE2/OZ48eRLr16+HQqHAxo0bMWPGDGRkZGDw4MFd+jmJepKMjAwYjUaEhoY2azcYDOjduzcAwGQywWAw4PPPPzfv98knn2DEiBHIyclp860qhhsi6na8vb3h4uICd3d3BAYGAgCys7Ov2a+hoQEffPABBg4cCAC444478MUXX0Cn08HDwwMRERGYNGkSdu7ciXnz5pl/cywsLIRCoQAAPPHEE0hNTcXatWvx8ssvd92HJOphqqurIZFIcPToUUgkkmbbPDw8AAByuRxOTk7NAlB4eDiApp4fhhsicnju7u7mYAMAMpkMISEh5n8oL7cVFxcDaNtvjkRkHcOGDYPRaERxcTHGjx/f4j7jxo1DY2MjTp8+bf67nZubCwAIDg5u83sx3BCR3XJ2dm72WiQStdhmMpkAtO03RyJqv+rqauTl5ZlfFxQUIC0tDX5+fggNDcX8+fOxYMECvPHGGxg2bBhKSkqwfft2xMTEICEhAVOnTsXw4cOxcOFCrFq1CiaTCYsWLcItt9xyzS8lrWG4IaJuycXFpdlA4M7Qlt8ciaj9jhw5gkmTJplfL126FACQlJSEdevWYe3atXjppZfw+OOP4/z58/D398eYMWOQmJgIABCLxdi8eTMWL16Mm2++Gb169cLMmTPxxhtvWFQHww0RdUshISE4ePAgzpw5Aw8PD3PvS0e05TdHImq/iRMnorXp85ydnbFixQqsWLHiuvsoFAp89913HaqDj4ITUbf0xBNPQCKRICIiAn369EFhYWGnnHft2rVYsGABHn/8cYSFhWHOnDk4fPgw+vXr1ynnJyLb4wzFRERE5FDYc0NEREQOheGGiIiIHArDDRERETkUhhsiIiJyKAw3RERE5FAYboiIiMihMNwQERGRQ2G4ISIiIofCcENEREQOheGGiIiIHArDDRERETmU/weoOq9dNqTxrgAAAABJRU5ErkJggg==\n",
"text/plain": [
"<Figure size 640x480 with 2 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"fig, axs = plt.subplots(2,1, sharex=True)\n",
"axs[1].set_xlabel(\"time\")\n",
"axs[0].set_ylabel(\"semi-major axis\")\n",
"axs[1].set_ylabel(\"stellar mass\")\n",
"axs[1].plot(times, mass)\n",
"axs[0].plot(times, semias)"
]
},
{
"cell_type": "code",
"execution_count": null,
"id": "599dc8c0-0821-46c2-9526-0da9f2abbde6",
"metadata": {},
"outputs": [],
"source": []
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3 (ipykernel)",
"language": "python",
"name": "python3"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.9.6"
}
},
"nbformat": 4,
"nbformat_minor": 5
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment