Skip to content

Instantly share code, notes, and snippets.

@TTY6335
Created September 19, 2020 13:31
Show Gist options
  • Save TTY6335/a4f42faf7bd92894a7366ada82cb76f0 to your computer and use it in GitHub Desktop.
Save TTY6335/a4f42faf7bd92894a7366ada82cb76f0 to your computer and use it in GitHub Desktop.
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "markdown",
"metadata": {},
"source": [
"http://data.astropy.org/ より計算に必要なデータを取ってくる。 \n",
"astropyを初回使用するときのみでOK"
]
},
{
"cell_type": "code",
"execution_count": 1,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"'http://data.astropy.org/'"
]
},
"execution_count": 1,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"from astropy.utils.data import conf\n",
"conf.dataurl"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {},
"outputs": [],
"source": [
"#conf.remote_timeout = 10000"
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {},
"outputs": [],
"source": [
"from astropy import units as u\n",
"from astropy.time import Time, TimeDelta\n",
"from astropy.coordinates import solar_system_ephemeris\n",
"solar_system_ephemeris.set(\"jpl\")\n",
"\n",
"#惑星がclassオブジェクトとして定義されている\n",
"from poliastro.bodies import Sun, Earth, Mars\n",
"from poliastro.ephem import Ephem\n",
"from poliastro.frames import Planes\n",
"from poliastro.twobody import Orbit\n",
"from poliastro.plotting import StaticOrbitPlotter\n",
"from poliastro.plotting.misc import plot_solar_system\n",
"from poliastro.util import time_range, norm"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"datetime形式で計算したい時刻を指定。 "
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {},
"outputs": [],
"source": [
"EPOCH = Time(\"2020-10-06 00:00:00\", scale=\"tdb\")"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"描写"
]
},
{
"cell_type": "code",
"execution_count": 5,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"([<matplotlib.lines.Line2D at 0x23c52906a00>],\n",
" <matplotlib.lines.Line2D at 0x23c517aa670>)"
]
},
"execution_count": 5,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAjUAAAGBCAYAAACepkBmAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nOzdd3hUZfo38O+ZPslMJj2TOmmTSgKa0EsQBIL0ruAqikhXBFZl2R9FxWVVEGTFRUCKq4IiIAJGERBRFEmQOum99zqZTD3vHxBehISWcjLJ/bkuL+HMKd9JQs49z3kKw7IsCCGEEEKsHY/rAIQQQgghrYGKGkIIIYR0ClTUEEIIIaRToKKGEEIIIZ0CFTWEEEII6RSoqCGEEEJIp2D1RQ3DMJ8wDFPCMMzV+9jXh2GYUwzD/MkwzGWGYZ5oj4yEEEIIaXtWX9QA2AUg9j73/SeAL1mWfQTAkwC2tFUoQgghhLQvqy9qWJb9GUDFrdsYhglgGCaOYZgEhmHOMAwT0rg7ALsbf1YAKGjHqIQQQghpQwKuA7SRjwHMZVk2lWGY3rjeIjMEwGoAPzAMswiALYDHuYtICCGEkNbU6YoahmFkAPoB+IphmMbN4hv/fwrALpZl1zMM0xfApwzDdGNZ1sJBVEIIIYS0ok5X1OD6I7UqlmV7NPHaLNzof8Oy7G8Mw0gAOAMoacd8hBBCCGkDVt+n5nYsy9YAyGQYZgoAMNd1v/FyDoChN7aHApAAKOUkKCGEEEJaFWPtq3QzDPMFgMG43uJSDGAVgJMAPgLgDkAIYC/Lsm8wDBMGYBsAGa53Gn6VZdkfuMhNCCGEkNZl9UUNIYQQQgjQCR8/EUIIIaRroqKGEEIIIZ2CVY9+cnZ2Zn19fbmOQQghHUZCQkIZy7IurXxOV4FAsB1AN9CHYcIdC4CrJpPphaioqCZHLVt1UePr64v4+HiuYxBCSIfBMEx2a59TIBBsVyqVoS4uLpU8Ho86YhJOWCwWprS0NKyoqGg7gLFN7UMVNyGEkHvp5uLiUkMFDeESj8djXVxcqnG9xbDpfdoxDyGEEOvEo4KGdAQ3fg6brV2oqCGEEGI1GIaJmj17tlfj31euXOm2ZMkSDy4ztdSRI0fkjz32WGB7X3fJkiUeK1eudGvv67YlKmoIIYRYDZFIxB47dsyhsLDQqvuEkrZBRQ0hhBCrwefz2Weeeab07bffvqOF4fPPP1dERkaGhIaGhvXr1y8oNzdXAFxvkZg4caJv//791Z6enhG7d++2nzt3rldQUFDYwIED1Xq9ngGAM2fO2PTs2TM4PDw8dMCAAers7GwhALz11luuAQEB4UFBQWGjR4/2v/26ycnJoqioqOCwsLDQsLCw0OPHj9sC11tgevXqFRwbG+vv5+cXPnbsWD+L5fr6yfv377fz8/MLj4qKCt6/f799U+/VZDJhzpw5Xt26dQsNCgoKe/fdd50bzxsdHR08bNiwgICAgPDp06f7mM1mAMDWrVsdg4KCwtRqdfi8efM8G8+1f/9+u7CwsNDg4OCwvn37BjVuT0xMlPbq1SvYy8sr4q233nIFgJqaGt7gwYMDg4ODw9Rqdfi2bdscHu671f6oqCGEEGJV/v73v5ccOHDAsby8nH/r9mHDhtVdvHgxKTExUTN58uSKN954Q9n4WnZ2tvjkyZNp+/fvT5s7d67fkCFDalJSUjQSicTy5ZdfKvR6PfPSSy/5fPPNN+nXrl1LfPbZZ8uWLVvmCQAffPCB8urVq5qUlBTNrl277hhd5uHhYTpz5kyKRqNJ3LdvX8Yrr7zi0/haYmKi9MMPP8xNS0u7lpOTIz5+/Lisvr6eWbhwoe/hw4fTzp8/n1xSUiJs6n1u3LjRWaFQmK9evZp46dKlxN27d7skJSWJAODKlSu2mzZtyk1OTr6WlZUl3rNnj0NWVpZw9erVnj/99FOKRqO59ueff9p++umn9gUFBYKFCxf6HjhwID05OVlz6NCh9MZrpKWlSU6fPp1y/vz5xPfee89Dr9czBw4csFMqlcbk5GRNamrqtYkTJ9a0/LvWPqj5jhBCiFVxdHS0TJkypXzdunWuUqnU0rg9MzNTNH78eK/S0lKhwWDgeXt76xtfe/zxx6vFYjHbq1cvndlsZiZPnlwDAOHh4brMzEzR5cuXxampqdIhQ4YEAYDFYoGLi4sRAIKDg3UTJkzwGzt2bNWMGTOqbs9jMBiYWbNmqTQajZTH4yE7O1vc+FpERIQ2ICDAeONa9enp6SK5XG728vLSR0RE6AFgxowZ5du3b79jbqEff/zRLikpyebw4cMOAFBbW8vXaDQSkUjERkREaMPCwgwAMHXq1IozZ87IhEIh26dPn1oPDw8TAEybNq3i9OnTMj6fz/bq1as2JCTEAABubm7mxmsMHz68SiqVslKp1OTo6GjMy8sTPProo7oVK1Z4z5s3z3PcuHHVsbGxdS35frUnaqkhhBBidZYvX178+eefO2u12pv3sYULF/rMnz+/JCUlRfOf//wnW6/X33xNLBazAMDn8yEQCFge7/pLPB4PJpOJYVmWCQwM1CUlJWmSkpI0KSkpml9//TUVAE6dOpW6YMGC0oSEBNvu3buHGY3Gv2RZu3atm6urqzExMVFz5coVjdFovOO6jdc2mUwMADAMc8/3yLIss379+pzGTPn5+VcaW01uP55hGDS3liPLss1er6l8kZGR+gsXLmgiIiJ0K1as8Fy2bJn7PcN2EFTUEEIIsTpubm7mMWPGVH7++efOjdtqa2v5Pj4+RgDYtWuX04OcLzIysqGiokLw448/2gKAXq9n4uPjJWazGenp6aIxY8bUbtmyJa+2tpZfXV39l8de1dXVfHd3dyOfz8eWLVucGvu3NKdHjx4NeXl5omvXrokBYO/evY5N7Tds2LDqjz76yKWxz8/ly5fFNTU1POD646ekpCSR2WzG/v37HQcOHFg7aNAg7blz5+SFhYUCk8mEr776ynHw4MF1jz32mPbcuXPyxkdXxcXF/Kau1ygrK0sol8st8+fPr1i8eHHxxYsXbe7zy8g5evxECCHEKq1YsaJo9+7dLrf8veCpp54KcHNzM0RHR2tzcnLEdzv+VhKJhN27d2/6Sy+95FNbW8s3m83MvHnziiMiIvTTp0/3q62t5bMsy8yZM6fY2dn5L1XL4sWLSyZNmhRw6NAhhwEDBtTe+kisKTY2NuzmzZuzR48eHejo6Gjq3bt3XWJiovT2/V555ZWyrKwscURERCjLsoyjo6Px2LFj6QDQo0ePuqVLl3olJSVJe/fuXfu3v/2tis/nY+XKlfkxMTFBLMsyQ4cOrX766aerAOCDDz7ImjBhQqDFYoGTk5Px7Nmzqc3lS0hIkC5fvtyLx+NBIBCwW7ZsafVZqtsK01xzlTWIjo5maZkEQgj5/xiGSWBZNro1z3np0qWs7t27l7XmOcnDO3LkiHz9+vVup06dSuM6CxcuXbrk3L17d9+mXqOWGkLugTWbYKqvg8XQANagh8VkgMVgAGs0wGIyAAAYhgcwDMDjXX92zeODL5KAJ5GCJ5aCL5aCJ5aA4d211ZcQQkgLUFFDuiSWZWGuq4ahshSGihIYqkphrCiBobIUptoqmLTVMGlrYdbWwKzTttp1eWIJBHIHCOX2EMgdILCzh1DuAKG9M8ROSoiclRA5uoEvsZpH2ISQdjZ69Oja0aNH13KdoyOiooZ0ahajHg1FOWgoykVDcS70RbloKM5BQ3EeLA31f9mX4QsgtHeG0M4BArkjJEoVBDI78G3tILCRgyeWgicUgScUgWn8v0AEMABYFrBYro8+YC1gzWZYDA0w63Ww6Btg0etg1utgrq+7XjTVVsFQUYz67GQYa6sAy187FgpkCoiclJAovSFx94XUwxcSd1+IXdyptYcQQppBRQ3pNMwN9ajPTYMuJxX1OSnX/1yQ9f8LBoaByNEVEjcfOPXrBrGLB0QOrhA5ukJk7wKBnQMYXvsPCGQtFphqK6EvL4KhrAiG8qKbf65NuYyKcz/e3JcRCCFxV8FGFQRbVQhsfUMg9fQHI6B/yoQQQr8JidUyVJSgLu0K6lIvoy7tCnQFmddbTAAI7Bxg4xMERWRfSL0CIFH6QOLqBZ7ovgdDtBuGx4NQ4QShwgnwD7/jdbNOi4aiHOgKstBQmAVdXgaq/vwF5b8cu368QAipVyBk/mGQBXWHXB0JgbzJWdcJIaRTo6KGWA1jdTlqNPGo0ZxHXeoVGCqKAQA8sRSygG6wj4qBrSoYNt5qCO0faIqKDo0vtYWtXyhs/UJvbmNZFoayQmizk1GflQRtVhJKzxxBycmvAQASDz/Ig7pDHtQD8tBHIbC14yo+IYS0GypqSIdlMRpQl34VNdfOo+baH9DlXV+uRCB3gDy4B9yGT4MsMAJSTz8w/K71o8wwDMQuHhC7eMAx+jEAgMVkRH12MmpTLqEu+SLKf4tD6U+HAIYHW/8wKLr1hiKiD6Tegfc1mykhhFgbmlGYdCjmhnpUJvyEjO1v4tLS8UjdsAQlP34Fga0dPCe+iND/24bId7+G/4ur4DpkImx81F2uoGkOTyCELKAb3EfOgHrxu+jx/hEEv/Yh3Ef9DazJiIJvdiDxrdm48upkZO15B9VXz4E1mbiOTTohs4VF3JUqxXtxhe5xV6oUZkvL5kNLS0sT9u7dO8jf3z88MDAw/M0333RtfK24uJjfr18/tUql6tavXz91aWkpHwAOHjxoFx4eHhoUFBQWHh4eevjwYXnjMWfOnLEJCgoK8/Hx6TZz5kzvxpWzb7do0SJPpVIZaWNj88it23U6HTNq1Ch/Hx+fbpGRkSHJycmipo6/236pqami/v37q/39/cMDAgLCmzpHc+8NAJYvX6708fHp5uvr2+3rr79usim2pccDQGxsrL9GoxEBgKenZ0RQUFBYSEhIWEhISNjMmTO9mzuuKUeOHJE3rmAOAJMmTfLduXPnHSuAFxQUCAYOHKh+kHM3oqKGcM5cX4fy335A2ocrcGnpeGRsXY3axAQ49hyCgAVvo/v7hxG09H0oY6fDxlvNSWdea8QIBJAFhMNj7HMIXbEVke8dgO/M1yFTR6Iy/iekffAaLi2bgKxd/0b1lXOwmIz3Pikh92C2sJi9KzNg6b6cwJ2/lHos3ZcTOHtXZkBLChuhUIj169fnZWRkXDt//nzijh07XBMSEiQAsGrVKvfBgwfXZmdnXx08eHDtypUrlQDg6upqPHr0aNqNlbUzX3jhBb/G882fP1+1ZcuW7KysrKsZGRmS/fv3N3lTHz9+fNW5c+cSb9++adMmZ4VCYcrJybm6cOHC4iVLlng1dfzd9psxY4bfsmXLijMyMq5duHAhsXERyls1994SEhIkBw4ccExOTr4WFxeXsnjxYh9TEx9QWnr8jWUimMaFMwHg9OnTKY1rUe3atSu3qffdFKPRiJMnT8rPnDkju9e+Hh4eJjc3N+MPP/xge699b0cfcQknWLMJNZp4lP/+A6ou/gLWaIDQwQUuA8fA/pGBkKkjaOhyKxPaOcKpXyyc+sXCYtSj5lo8Ki+cRuWFn1F+9jvwbWRwiH4MTv1iYesXRo+oyEM5fq1acS6j7i891c9l1Nkfv1atiI2wr36Yc6pUKqNKpTICgIODgyUgIECXk5MjioqKaoiLi7M/ffp0MgDMmTOnPCYmJhhAfv/+/XWNx0dFRTUYDAaeTqdjSkpKBHV1dbzHH39cC1xfIfvQoUMOU6dOrbn9ukOHDm1ykqojR47Yr169ugAAnnvuucrXXnvNx2KxgHfbB67m9vvzzz8lZrMZEyZMqAEAhULRZFNRc+9t//799hMnTqyQSqVsSEiIQaVS6X/66SfbxvfUWsfv2rXLacyYMXesSn679evXO+/cudPFaDQyvr6++v3792fK5XLLpEmTfB0cHExXrlyxsbe3NyckJMh4PB775ZdfOm3cuDEHAE6fPi374IMP3EpLS4Vvvvlm3nPPPVcJXC8o9+zZ4zR8+PAHmiiMihrSrupz01D+2/eo+ONHmGoqwbe1g3P/J+DYZxjdSNsRTyiGfY/+sO/RHxajATWJ8ag8fwrlv/2Asp+/hdjNG879YuHYZzhEDi73PiEhN1zN1zU5c+S1Ap30YYuaWyUnJ4s0Go1NTExMHQCUl5cLGgselUplrKiouOO+tnv3boewsLB6qVTKZmdnC93d3W82S6pUKkNhYaHwQTIUFxeL/Pz8DMD1ViSZTGYuLi4WuLu7m+5nP41GI7GzszMPHz48IDc3Vzxo0KCaDz/8ME9w29QMzb23/Px8UZ8+feoa9/Pw8DDk5uaKAGinTZumWrBgQemgQYPqH+b4W69/7tw52TPPPFNx67aYmJigxuLtqaeeKlu1alXJjBkzKpcuXVoGAC+99JLHBx984LxixYoSAEhPT5f8+uuvKQKBAEuWLPGQyWTmN954oxgAtm3b5lxcXCyMj49PunjxomTChAmBjUVN//79tW+88YbHg3xfACpqSDuwGPWojP8JpT8dgjYzEQxfAEVkPzj1HQ67br3BEzzQ7xPSynhCEewj+8E+sh98pi9GZcJplP8Wh/yD25B/aAfswnvB9bHxsAvvRY/+yD1185TWN7U93EOqa2r7g6iuruZNnDgxYN26dbmOjo53XTSyUXx8vGTlypWecXFxqcD1kYO3e9APU82c446Nze1nMpmY+Ph42blz5zRqtdowevTogM2bNzu/8sor97W+1t2uv2/fvnsuPnm/+UtLS4VKpfIvz6VPnz6dcnvxlpCQIF25cqVnbW0tX6vV8mNiYm4WrxMnTqy8vVi71dixY6v4fD6ioqIaysvLb94MPDw8TCUlJU32VbobKmpIm9GX5KP058Mo+/U7mLU1kCh94D1tERz7DKMhxh0UX2oL5wFPwHnAE2goyUP5b9+j7MxRpG1+HWIXD7jEjINT/5H0/SPNGhauqO7tL6u69RFUb39Z1bBwRYtaafR6PTNq1KiAKVOmVDz77LM3H4k4OTmZsrOzhSqVypidnS10dHS8ecNNT08XTp48OXDHjh2Z4eHhegDw9fU13toyk52dLVIqlUaTyYRu3bqFAUBsbGzVxo0bC5rLolQqDZmZmaKAgACj0WhEXV0d39XV1bxo0SLP48ePKwAgKSlJ09x+Pj4+htDQUF1jX5WxY8dW/v7773f0NWnuvXl5eTW2rAAACgoKRF5eXnd0imvp8WKx2KLT6e75SebFF1/0279/f1rfvn11H3zwgdPp06dvdsqWyWR3LT4lEsnNYurWYqu+vp4Ri8X3Vbjeij52kVbFsixqUy4hbfPruPrPGSj+cT/kwY8gaMn7CFuzG65DJ9EN0UpIXL3gOW4WItbtg9/slRAqnJC3/yNcfnUKsve8i4bCe34gJF0Qn8dg20y/9PXTfNKeH+iSv36aT9q2mX7pfN7DP1q2WCx48sknVUFBQQ2rV68uvvW1ESNGVG3dutUJALZu3eoUGxtbBQBlZWX8J554Qr169eq8W/tlqFQqo62treXEiRO2FosFn332mdO4ceOqBAIBGjvA3q2gAYBRo0ZVffLJJ04AsHPnToe+ffvW8ng8bN68Ob/xHHfbLyYmRltdXc0vKCgQAMCpU6fswsLC7mjJau69TZo0qerAgQOOOp2OSUpKEmVlZUkGDx58R9+Tlh6vVqsbEhMT7zljaX19Pc/Hx8eo1+uZvXv3Oja3n1wuN9fW1t5XZ8mrV69KgoKCHrh1j1pqSKtgLWZU/fkLir7fi/qsRAhkCriPfhbOg8ZAZO/MdTzSAjyBEI49h8Cx5xDU56ah9KdDKP/9B5T9egz2PQZAGTv9LxMDEsLnMYiNsK9ujT40AHD8+HHZoUOHnNRqtS4kJCQMANasWZM/bdq06jVr1hROmDAhQKVSOXt4eBgOHTqUDgDvvPOOa05OjnjdunUe69at8wCAEydOpHh6epq2bNmSPWvWLL+Ghgbmscceq5kyZUqTOefOnet18OBBx4aGBp6bm1vkjBkzyjZs2FDw8ssvl02aNMnPx8enm0KhMO/bty+9qeOb208gEGDdunV5gwcPDgKAiIiI+qYePTX33qKjoxvGjx9fERQUFM7n87Fhw4bsxkc8t/apeZjjbzVy5MiqkydPysePH39z8cxb+9SEhobWHzx4MOv1118v6NWrV6inp6chNDS0vq6ursnCZdKkSVWTJ08O+O677+wbOwo35/jx4/LY2NgH/vlhmnq2Zi2io6PZ+Ph4rmN0aRajAeVn41B8fB/0JfkQOXtAOXwanPqOAE8s4ToeaSPGmkqUnDqA0lMHYa6vgyyoB5Qjp8MurCd19uYYwzAJLMtGt+Y5L126lNW9e/f76u9BOo+6ujqmf//+wQkJCUl36xfTFqKjo4O/++67NBcXF/Ptr126dMm5e/fuvk0dRy015KGwJhPKf4tD4dFPYagoho0qGP4vrob9owNpKHYXILRzgOe4WVCOeAplZ46g+PiXSNv0Kmz9w+AxbhbkIY9ScUOIlZPJZOzKlSsLMjMzRWq12nDvI1pHQUGB4OWXXy5uqqC5F2qpIQ+ENZtQce5HFBzZA0NZAWx8Q+Ax7nn6hN7FWUxGlJ+NQ+HRPTBWlkIW1AOe456HTB3JdbQuh1pqSGdHLTWkxViWRdWFn5F/aDv0xbmQeqsRsPBtKCL6UjFDwBMI4TJoDJz6DkfZmaMoPPYpkt99CXbhPeE5cQ5svAO5jkgI6QKoqCH3pM1KQt6XH6Iu7Qok7r7wn/sG7B8ZSMUMuQNPKIbrkIlw7v8ESn46hKK4z5H41mw49R8Jz3GzIFR0ntXTCSEdDxU1pFmGihLkH9yGinPHIZA7wOfppXDuP5IWkCT3xBNLoBzxJJwHjELh0T0oPXUQledPQTlyBtwenwKe6J6jRAkh5IHR3YncwWLUoyjuCxTFfQGwFihHzoAydjr40gdeW4x0cQJbObynLoBLzDjkf/1fFBzajrKfv4XXtIWw7zGAWvsIIa2KJt8jf1GjiYdmzfMo/HYX7Lv3Q/gbe+A5YTYVNKRFJG5eCJj/FoKWvg++jQwZH/0f0j/8B/RlhVxHI22AtVhQcf6UIm//f90rzp9SsJYHnhj2L9LS0oS9e/cO8vf3Dw8MDAx/8803XRtfKy4u5vfr10+tUqm69evXT11aWsoHgIMHD9qFh4eHBgUFhYWHh4cePnz45iy3Z86csQkKCgrz8fHpNnPmTG9LM/kWLVrkqVQqI21sbB65dbtOp2NGjRrl7+Pj0y0yMjIkOTm5yen877ZfamqqqH///mp/f//wgICA8KbO0dx7A4Dly5crfXx8uvn6+nb7+uuvm5zRtKXHA0BsbKy/RqMRAYCnp2dEVFRU8K2vh4SEhKnV6vDmjn8Qo0eP9r9y5UqLmnGpqCEAAGN1OTK2v4nUjcsAMFAvfg/+L66C2Nmd62ikE5EHP4LQFVvhNXkeapP/xLVVM1EU9zlYk+neBxOrwFosSN24LCBz25rA4h/2emRuWxOYunFZQEsKG6FQiPXr1+dlZGRcO3/+fOKOHTtcExISJACwatUq98GDB9dmZ2dfHTx4cO3KlSuVAODq6mo8evRoWkpKimbXrl2ZL7zwgl/j+ebPn6/asmVLdlZW1tWMjAzJ/v37m7ypjx8/vurcuXOJt2/ftGmTs0KhMOXk5FxduHBh8ZIlS7yaOv5u+82YMcNv2bJlxRkZGdcuXLiQ6OHhccc/gubeW0JCguTAgQOOycnJ1+Li4lIWL17sY2ri31BLj4+Pj5eYzWamcTkHANBqtfy0tDQhAFy4cOGBJyNr6jqN5s2bV7J27Vrlg57zVlTUdHGsxYLSn77BtZXPoOrCz3Af/SzCVu2AXVirjggl5CaGL4Db8GkIX70bivCeyD/wMTRvvQBtVhLX0UgrqEw4rahNumB/67bapAv2lQmnFQ97TpVKZRwwYEA9ADg4OFgCAgJ0OTk5IgCIi4uznzNnTjkAzJkzp/y7775zAID+/fvrfH19jQAQFRXVYDAYeDqdjsnOzhbW1dXxHn/8cS2Px8OMGTPKDx065NDUdYcOHaptXOX6VkeOHLF//vnnywHgueeeqzx79qy8qdae5vZLSEiQmM1mTJgwoQYAFAqFRS6X33GC5t7b/v377SdOnFghlUrZkJAQg0ql0v/00093NKe39Phdu3Y5jRkzpurWbePHj6/Ys2ePIwDs2bPHcdKkSTdX8U5OThZFRUUFh4WFhYaFhYUeP37c9sbXQd67d++gMWPG+AUHB4fX1NTwBg8eHBgcHBymVqvDt23b5gAAsbGxdWfOnLEzGu/4kt83Kmq6MH15EVI3LkPO5+/DRhWMsFWfwGPsc+AJqRMnaXsiJzcEzH8LAQvehrm+Dknr5iP/4HZYjO02xxdpA/XZyTZNbs9JlrbG+ZOTk0UajcYmJiamDgDKy8sFjYWHSqUyVlRU3NFXdPfu3Q5hYWH1UqmUzc7OFrq7u9+8a6pUKsOtC1zej+LiYpGfn58BuN6KJJPJzMXFxXdct7n9NBqNxM7Ozjx8+PCA0NDQsDlz5ng11YLR3HvLz88XeXt73/yH4uHhcXOBymnTpql+/vlnm4c9/lbnzp2T9enT5y+rrk+fPr3y22+/dQCA77//3n7ixIlVt5zHdObMmRSNRpO4b9++jFdeecWn8bXLly/bvvvuu/np6enXDhw4YKdUKo3Jycma1NTUaxMnTqwBAD6fD5VK1fD77783+TN0P6ijcBfEsizKfzmK3K+2ACx7fVTTwNHUaZNwwr57P8jUEcj7cguKvvsfqi79Cr/nlsNGFcR1NPIQbFTB9U1u9wl+4MUJb1ddXc2bOHFiwLp163IdHR3v63lWfHy8ZOXKlZ5xcXGpwF9Xgm70oL/7mjnHHRub289kMjHx8fGyc+fOadRqtWH06NEBmzdvdm5q/acHvf6+ffvuudLs/eYvLS0VKpXKvzSbuLi4mBUKhenjjz92CAwM1N26CrfBYGBmzZql0mg0Uh6Ph+zs7JufkCMjI7UhISEGAHj00Ud1K1as8J43b57nuHHjqmNjY+sa93N2djbl5uY+UJF5qw7TUsMwjDfDMKcYhklkGBQjNMQAACAASURBVOYawzAvc52pMzJUliJt8+vI/vQ92PgEIWzVJ3AZNIYKGsIpgY0cvjNfQ+DCf8GsrUHiv+ai4PBOsGbqa2NtHKJiquUhj/7lkYU85NEqh6iYFi1uqdfrmVGjRgVMmTKl4tlnn715ficnJ1N2drYQALKzs4WOjo43f2jS09OFkydPDtyxY0dmeHi4HgB8fX2Nt7bMZGdni5RKpdFkMiEkJCQsJCQkbPHixR53y6JUKg2ZmZkiADAajairq+O7urqaFy1a5Nl4jrvt5+PjYwgNDdWFhYUZhEIhxo4dW3nhwoU7Wieae29eXl5/aVkpKCgQeXl53fHMpqXHi8Vii06nu6NOmDx5cuWrr76qeuqppypu3b527Vo3V1dXY2JioubKlSsao9F481gbG5ubxU9kZKT+woULmoiICN2KFSs8ly1bdrPzpl6v592674PqMEUNABOApSzLhgLoA2ABwzBhHGfqVKr+PAPNmudRm3wR3k++hKAlG6gjMOlQFJF9EbZ6Fxx7DkHhkd1IWf8KDOXFXMciD4Dh8aBe/F663+xVaW4jnsz3m70qTb34vXSG9/C3G4vFgieffFIVFBTUsHr16r/8QIwYMaJq69atTgCwdetWp9jY2CoAKCsr4z/xxBPq1atX5w0fPlzbuL9KpTLa2tpaTpw4YWuxWPDZZ585jRs3rkogECApKUmTlJSk2bhxY8Hd8owaNarqk08+cQKAnTt3OvTt27eWx+Nh8+bN+Y3nuNt+MTEx2urqan5BQYEAAE6dOmUXFhZ2R0tWc+9t0qRJVQcOHHDU6XRMUlKSKCsrSzJ48GBtax+vVqsbEhMT7+iPMGPGjMoFCxYUNT42alRdXc13d3c38vl8bNmyxclsbnrppqysLKFcLrfMnz+/YvHixcUXL168WdBlZmaKH3nkkYZmvvT31GEeP7EsWwig8MafaxmGSQTgCUDDabBOwGLUI2//f1F66iBsVEHwe2ElJG5NdtYnhHMCWzn8Zv0Tdt16I+d/G6B58wWonn0VDo8M5DoauU8MjwfHno9VO/Z8rEWtM42OHz8uO3TokJNardY1toKsWbMmf9q0adVr1qwpnDBhQoBKpXL28PAwHDp0KB0A3nnnHdecnBzxunXrPNatW+cBACdOnEjx9PQ0bdmyJXvWrFl+DQ0NzGOPPVYzZcqUJnPOnTvX6+DBg44NDQ08Nze3yBkzZpRt2LCh4OWXXy6bNGmSn4+PTzeFQmHet29felPHN7efQCDAunXr8gYPHhwEABEREfVNPXpq7r1FR0c3jB8/viIoKCicz+djw4YN2Y2raE+bNk21YMGC0kGDBtU/zPG3GjlyZNXJkyfl48ePr711u4ODg2Xt2rVFt++/ePHikkmTJgUcOnTIYcCAAbVSqbTJFpeEhATp8uXLvXg8HgQCAbtly5ZsAMjNzRWIxWK2qc7Z96tDLmjJMIwvgJ8BdGNZtua2114E8CIA+Pj4RGVn3/PxYZfWUJyLjI/XQJebBtfHp8Bz4ovgCR76cSUh7aqhJA+Z295AfXYKXGLGwWvKfJqN+B5oQUvSWurq6pj+/fsHJyQkJDVV9LS2NWvWuNrZ2Vnu1bfobgtadqTHTwAAhmFkAL4GsPj2ggYAWJb9mGXZaJZlo11cXNo/oBUpP3cciW/NhqGiBAEL34b31AVU0BCrInH1QvBrH8Jt2FSUnv4Gye8sgr78jg+IhJA2IJPJ2JUrVxY09gtqa/b29uaFCxe2qHjuUEUNwzBCXC9oPmNZ9gDXeawVazIhd99mZO1Ye70z8MrtsI/sx3UsQh4KTyCE15T5CFj4NvSl+Uh860XUaOK5jkVIlzBp0qQatVrdLvMsvPzyy+VCYcs+eHeYooa5PvxmB4BElmU3cJ3HWhlrq5CycRlKTnwN16GTEbRkA0QOrvc+kJAOzj6yH0L+sRVChSNSN72KorgvmhyaStqExWKx0BBJwrkbP4fNjo7qMEUNgP4A/gZgCMMwF2/89wTXoaxJfU4qktbOgTbjGnyfWw7vaQtpRW3SqUjcvBDy+hY4RMUg/8BWZH68Bhb9Qw+UIPfvamlpqYIKG8Ili8XClJaWKgBcbW6fDnPHY1n2FwD0D+YhVfxxAlm734FAZofgVzfD1jeE60iEtAm+xAZ+s1fCRhWM/ANboS8vRuCCtyBUOHEdrdMymUwvFBUVbS8qKuqGjvVhmHQtFgBXTSbTC83t0CFHP92v6OhoNj6+az9bZ1kWRXGfo+DgNsgCI+A/dw2Edo5cxyKkXVRd/AWZ29+CQKZA4KJ/Qerpz3UkzrXF6CdCrAVV3FaMNZuQ87/1KDi4DQ69hkL9ynoqaEiXYt9jAIL+vgms2YSkfy9C9bU/uI5ECOEQFTVWytxQj7QPV6DszBEoR86A3/MrwBO2y6g7QjoUW1UwQpZ/BLGzEmmbX0fZ2e+4jkQI4QgVNVbIWF2O5HdfRs218/CZsQSeE2ajJVOQE2LtRI6uCH51M+TBjyJ7179RfPxLriMRQjhAd0Iroy8vQvK7L0FfnIvAhW/DJWYs15EI6RD4EhsELnwb9o/GIO+rLcg/uJ2GfBPSxVBRY0UaivOQ/M5LMNVWQf3Keigi+nAdiZAOhScUwf/FlXAeOBpF3/0POZ9tAGtpelE9Qkjn02GGdJO70+VlIGXjUsBiQdDSjbDxUXMdiZAOieHx4fP0Ughs7VAU9zkseh18n1sOhsfnOhohpI1RUWMFtFlJSN30KnhCEdRLN0LqruI6EiEdGsMw1xdvldig4NB2AKDChpAugIqaDk6blYSUDUshkNkh6JX1ELt4cB2JEKvh/sTTAMui4JsdAKiwIaSzo6KmA6vPTUPqxr9DYCtH8LJNEDnSGk6EPCj3UX8DACpsCOkCqKjpoHQFWUjduAw8sQRBS9+ngoaQFri1sGH4QqiefRXX19AlhHQmVNR0QA3FeUh9fynA8BC0ZAPEzu5cRyLE6rmP+htYswmFR3ZDYCuH5+R5VNgQ0slQUdPBGCpKkLJhCVizCUHLNkHi5s11JEI6DfcxM2HS1qD4+JfgyxRwHzmD60iEkFZERU0HYtLWInXTq7A0aBG0bBOkHr5cRyKkU2EYBt7TFsGsrUXBwW0Q2MhpAktCOhEqajoIi1GP9C0roC/Nh/rld2DjHch1JEI6JYbHg+/M12Gqr0PO5+9DIFPAISqG61iEkFZAMwp3AKzFjMwda1GXehm+z/8D8uBHuI5ESKfGCAQImLMatn6hyPxkLbSZiVxHIoS0AipqOMayLHL3/QdVF36G19QFcIx+jOtIhHQJPLEEAfPXQqhwQtqH/4C+vIjrSISQFqKihmMlJ79G6amDcBs2FW6PT+E6DiFditDOAYGL1oE1GpC2eTnM9XVcRyKEtAAVNRyquXYeeV9ugX2PAfCcNJfrOIR0SVJ3FfznvoGGohxkbFsD1mziOhIh5CFRUcORhuI8ZGxbA6mHL3yf/wcYHn0rCOGKXWgUVE8vQc2188g/uJ3rOISQh0Sjnzhg1mmRvmUFwPAQsGAt+BIbriMR0uU5DxiF+pxUFP+wFza+wdS/jRArRM0D7Yy1mJG5/U00FOchYO4ami2YkA7Ea+oC2AZ0Q/buf0OXn8F1HELIA6Kipp0VHvsfqq/8Du8nF9HQbUI6GJ5ACP85q8ET2yD9o/+Dqb6W60iEkAdARU07qkm6gMJvd8Gxz3C4xIzjOg4hpAkie2cEzFkNfVkRsnauA8uyXEcihNwnKmraibG6HJnb34TEzRs+0xfTQnqEdGAydSS8Js1F9aVfUXrqINdxCCH3iYqadtA4Y7BZVw//OaupYzAhVsD18cmw69YHefs/Qn1uGtdxCCH3gUY/tYPCo5+iNukCVM++BqmnP9dxuhyDyYLCKiMKqw0oqDKisMqAcq0ZVfUmVGpNqNKZUa+3oMFogd5kgd7EwmxhwTAMeAzAYwAhn4GtiA8bMQ+2Yh7kEj6cZQI4yQRwlgnhLBfA21EEb0cxHGz41BLXCTAMA9/nXodmzfPI3PYGQlZsBV8s5ToWIeQuqKhpY3VpV1B4ZA8c+wyHU79YruN0ahYLi8wyPTQFOqSVNCC9RI/00gbkVRhguaVbBMMA9lI+7G0EsLfhw9NeBJmYB7GQB4mQgUjAg4DHgGVZWFjAwgJGswVa/Y3/DGbU6MzIKtOjrM4Eg+mvfS5sxTz4OIqgdpNA7SZBkFKKYKUELnJhO39FSEsJ5fbwm7UCqRuXIXffZvg+8yrXkQghd0FFTRsyN9Qja+e/IHJyo340baCuwYwL2Vqcz9LiSm49NIU6aPUWAICAz8DXSYQQpRRPRNjD21EED3sR3O2FcLMTQiRonSevLMuiTm9BcbUReZUG5FYYkFuhR3a5AecytDh8sermvkqFED28bRDpbYNHfGwR5iGFgE8/Ex2dXWgUlLHTUfTdZ7CP7A/7Hv25jkQIaQYVNW0o7+v/Ql9WiKBlm6gfTSswmllczNHi55Ra/JFRB02BDhb2egET6i7B2B4OCPeUIsLTBr7O4nYpGBiGgVzCh1zCR6Cb5I7Xq+pNSClqQGKhDlfy6nExpx5xV6sBADIxDz39ZOgbKEO/ABl8ncVU+HZQ7mNmovrK78j+33rIArtBIFNwHYkQ0gTGmocrRkdHs/Hx8VzHaFL11XNI++A1uA2fBq/J87iOY7W0ejNOJdXgp6Qa/Jpah5oGMwR8Bt29bNDTzxY9/WTo7m0Dqch6+ryX1BjxZ44Wv6fX4be0OuRWGgAAKicRhoQqMDTMDt29bMDjUYHTkdTnpiJx7Vw4Rj8Gvxf+yXWcZjEMk8CybDTXOQjhAhU1bcCkrYFmzfPg28gQumIreEIx15Gsis5gwenkGnx3pQpnUmqhN7FwshVgULAcMcF26Bsgg0zC5zpmq8mt0OPX1DqcTKrGuQwtTGYWzjIBYiPsMbaHPcI8pNSC00EUHNmNwsM74T/3DTg8OojrOE2iooZ0ZVTUtIHMT95GxR8nEPqP/8LGR811HKvAsiyu5OtwIL4Cx65UQau3wFkmwPBuCsR2s8cjPl2j5aK2wYwzKbX4/moVTifXwmhm4e8ixpjuDpgQ5UCdjTnGmkxI/Nc8GKtKEb56FwRye64j3YGKGtKVUVHTymo08UjduAzuo56Bx7jnuY7T4Wn1ZnzzZyW+Ol+BlOIGSIQMRnSzx7hHHBDtawt+FyhkmlOtM+GHq9U4cqkK8VlaCHjAkFAFpvVyQm9/W2q94YguLwOat2bDqe8I+D7b8UZDUVFDujLqKNyKLAY9cj5/H2JXLyifmMF1nA6toNKAz8+VYX98BWobLAj3kGLVWE+MjLSHvBM9WmoJhVSAKT2dMKWnE7LK9PjqfDkOXqjED9eq4e8ixsz+LhjTw77VRnKR+yP18ofbsKko/v4LOPWLhVwdyXUkQsgN1FLTivK/2YGio59CvWQD7EIe5TpOh5Re0oCPT5fg2OUqMAwwLEyBv/VzRndvG2p5uA8NRgu+v1qNT8+WIrGwAS5yAWb0cca0Xk6wk1Ix2F7Meh00q2aCJ5Ei9J/bwBN0nMeC1FJDujIqalqJriALiW++AIdeQ+H33HKu43Q4yUU6bP2pBD9cq4ZEwGBqLyc83dcZHvYirqNZJZZl8Xt6HT45U4qz6XWwk/Axc4Aznu7rDFsxFTftoerSWaR/+A94TnwRytjpXMe5iYoa0pXR46dWwLIscv63HjyJDQ3fvk1ehQGbTxThyKUq2Ip5mD3IFc/0c4aDLf3otQTDMOgbKEffQDk0BfX48GQxPvixGHvOlmHWQBdM7+MMiZAeS7Ul++79YN9jAAq+3Q2HnkMgdlJyHYmQLo9+67WCyvMnUZd2BV4TX4SwA46G4EKl1oR/HS3AqE3J+FFTjdmDXHB8WQheHqakgqaVhXnY4MOn/fDF3ECEe0qx/vsijN6YjKOXKmHNLbHWwHvaIgAs8g9u4zoKIQRU1LSYxaBH3tdbIfVWw6n/SK7jcM5sYbH3XDme2JiMz38vw7geDjj2SggWD3eHQkrFTFuK9LLBx8/6Y+fz/rC34ePVr3Ix4+N0XM6t5zpapyVycoPbsKmo/OME6jKucR2HkC6PipoWKj7+JYyVJfCeugAMr2v3ZbiYo8W0j1Lx5rf5CFFKcGBhEN6Y4AU3u47TibIr6OUvw755arw1wQv5lQZM/zgNbxzOQ43OzHW0TkkZOx0CO0fkfbmFWsYI4RgVNS1grCpHUdxnsH9kIOTBPbiOwxmt3oy3vs3HjI/TUaE1471pPvjkeX+om1gLibQPPo/BhChHHHslGH/r64yvzldg9KZkHLtcRTfeVsaX2MBz/AvQZlxD5fmTXMchpEujoqYF8r/ZAdZkguekOVxH4czZtFqM35yCvX+U4+m+zvj25SCMjLCn4dkdhK2Yj9ee8MC+eYFQ2gnx9y9zsPCzLJTVGbmO1qk49RsBqXcg8g98DItRz3UcQrosKmoeki4vA+Vnv4PrkImQuHpxHafd6QwWrPkmD7N3ZUIs4OHT2QFYPsqDhhN3UGEeNvhibiD+PtIdZ9PqMP6DFPxwrZrrWJ0Gw+PDa9JcGCqKUXbmKNdxCOmyqKh5SAXf7gRPbAPlE09zHaXdJRboMGVLKr6Kr8BzA1zw9QI1HvGx5ToWuQc+j8HM/i74ar4aHvYivPJFNpbvz4VWT31tWoM8NAqyoO4oPPY/WPQNXMchpEuiouYhaLOTUfXnGbgNmwKBrR3XcdoNy7LY/WspntyaBq3ejO0z/bEs1h1img/FqgS6SvDZnEDMe8wVRy5VYtpHaUgp0nEdy+oxDAOPcbNgqqlAyelDXMchpEuiu9FDKPhmJ/i2dnAbOpnrKO1Gqzdjyd4cvPNdIQYFyXFwURD6BMi4jkUekpDPYOFQJXY854/aBjOe/G8avjpfTp2IW0iujoRdWE8Ux30BcwMNpSekvVFR84Dq0q+h5urvUA6fBr5N17ipp5c04Mn/puFHTTWWjlDig+kq2NvQnDOdQS9/Gb5eqMajKlus/iYfq7/Jh8Fk4TqWVfMY9zxMddUoOfE111EI6XKoqHlABYc/gUBuD5fHJnAdpV2cSanBU1vTUFVvwvbn/PH8QFca2dTJOMuE2PqsH2YPcsH++ArM2pmB8joT17Gslq1fKBQRfVByYj/1rSGknVFR8wDqMq6hNjEByhFPgS+x4TpOm/vs9zLM/zQL3o4ifDlPjd7+XaNlqivi8xgsHu6Od6f6QFOgw9SPUpFM/WwemjJ2Okx11Sj79RjXUQjpUqioeQDF3+8D30YO50FjuI7SpswWFm8fycfbRwoQE2yHPS8EwJ1W0+4Snoi0x6ezA8CywDPb0vF7eh3XkaySTB0J24BuKD6+D6yJWr0IaS9U1NynhuI8VF08A5eYcZ26lcZgsuDvX+bgs9/L8Uw/Z2yarqK5Z7qYMA8bfPZiAJQKIebsycSxy1VcR7JKytjpMJQXoyLhFNdRCOkyqKi5T8XH94HhC+A6pPP2pdHqzVjwaRa+v1qNv8e647UnPMDnUf+ZrsjdXoQ9swPQw9vmepH7WxnXkayOIqIPJO6+KI77gkaVEdJOqKi5D8aaCpSfjYNT3xEQKpy4jtMmanRmzN6Vid8z6vDWBC/MHODCdSTCMYVUgI+f9cPQUDu8fbQAn5wp4TqSVWF4PLiNmAZdfgZqky5wHYeQLoGKmvtQcuogWLMJbsOmch2lTdQ2mDF7VwauFejw/lMqTIhy5DoS6SDEQh7WP6nCyAgF1n9fhI9OFXMdyao49hwCgUyB0p9oMj5C2gNNNnIPFqMBZacPQxHZDxKlD9dxWl1tgxkv7spEUlED3n/SB0NCFVxHIh2MkM/g31N8IOTn4T8nisGywPwhblzHsgo8oRhOA55A8ff7YKgogcjRletIhHRqHaqlhmGYTxiGKWEY5irXWRpV/XkGprpquAwex3WUVqfVmzF3dyY0BfVU0JC74vMYrJ3ohfGPOODDk8XY/Wsp15GshsugsQBYlP58mOsohHR6HaqoAbALQCzXIW5VevowRM4esAuN5jpKqzKYLHhlbzYu59XjvWkqKmjIPfF4DNaM98LwcAXe+a4QX50v5zqSVRA7u0MR2Q9lZ47CYjRwHYeQTq1DFTUsy/4MoILrHI10BZmoS70El0FjwPA61JeqRSwWFisO5OHX1DqsGe+FYeFU0JD7I+Az+PcUbwwMkmPN4Xx8f5WGe98P18HjYaqtRNWfZ7iOQkinZnV3aoZhXmQYJp5hmPjS0rZtAi/9+VswAiGc+o9s0+u0t3fiCnHschVeGa7EROoUTB6QSMDDxqdU6O5tg+X7c3ExR8t1pA5PHhoFkaMbyn+L4zoKIZ2a1RU1LMt+zLJsNMuy0S4ubTfs2KzXoeK37+HwaAyEcvs2u05723uuHJ+eLcPTfZ0xayAN2yYPRyLk4T8zfOGmEGLh/7KQU67nOlKHxvB4cOo7AjWaBBgqqT8SIW3F6oqa9lJ18ReYdVo4DxzNdZRWczatFm8fzUdMsByvjnSnhSlJizjYCvDfv/mBBTDv00zU6MxcR+rQnPqOAFgLKn7/gesohHRaVNQ0o+LcjxA5ukGmjuQ6SqvIKtNjyd5sBLhI8O5UH5opmLQKlbMYH0z3RV6lEa99lQOLhWbObY7Y1RMydSTKzsbRDMOEtJEOVdQwDPMFgN8ABDMMk8cwzCwuchhrKlCjOQ/H3o93ig7C9QYLXv48CwIegw+f9qW1nEirivK1xT9GeeDnlFp8eJIm57sbp74joC/OhTZTw3UUQjqlDnXHZln2KZZl3VmWFbIs68Wy7A4uclSePwVYLHDsPYyLy7cqlmWx6lAe0kv1eHeqDzwcaLVt0vqm9nTExEcd8N+fSnAysZrrOB2WQ9RgMALh9d8xhJBW16GKmo6i/NyPkHoHQurhy3WUFvviXDmOXa7CoqFu6Bso5zoO6aQYhsE/x3gi3FOKFV/noaCS5mNpCl9qC7vwXqi8cBqsxcJ1HEI6HSpqbtNQnIv6rEQ4dYJWmuQiHd75rhCDguSYPYimZydtSyzk4b2pPjCzLF79KgcmM/UbaYrDozEwVpZCm5XIdRRCOh0qam5TmXAaAODQcwjHSVqmwWjBa1/lwk7Kx9pJ3uBRx2DSDnycxFg51hN/5tTT4pfNsO/eDwxfcPN3DSGk9VBRc5uqi7/AxjcUIgfrnsPl/R8KkVrcgLUTveFoS+uWkvYzursDxj/igK2nS/AnTcx3B76NDHZh0ahKOE2joAhpZVTU3MJQVYb6rCTY9+jPdZQWOZdRh//9Vo4ZfZwwMIj60ZD294/RHlDaCfF/B/LQYKS+I7ezj4qBoaIY9TkpXEchpFOhouYW1Zd+BQDY9xjAcZKHpzNYsOpgHnycRFgywp3rOKSLshXz8eYEL2SW6fGfE/QY6naKbn0AADVXznGchJDOhYqaW1Rd/BViV09I3FVcR3lo/zlRhNxKA94Y7wWJkL69hDt9A+WYEu2I3b+W4nJuPddxOhShnQNsfENQfZWKGkJaE931bjDrtKhNugD77v2tdvmAq/n12HO2DFN7OqKnn4zrOIRgWaw7XORCvHE4D2aabfgvFN16Q5uZCFMdzetDSGuhouaGmsQEsGYTFN2tsz+NxcLirW/z4WgroMdOpMOQSfh4daQ7EgsbsPePcq7jdCiKbr0B1oIaTTzXUQjpNKiouaE2MQE8sRQy/3CuozyUb/6sxJU8HZaOcIdcQssgkI5jRDcF+gbIsPnHIpTVGbmO02HY+AaDb2tHj6AIaUVU1NxQk3QBMnV3MALrG/5cozNjww+F6OFjgzE97LmOQ8hfMAyDFaM9oDOy2PRDEddxOgyGx4ddaBRqky7Q0G5CWgkVNQAMlSXQF+fCLvRRrqM8lI9Pl6Cy3owVoz2stj8Q6dz8XCSY0ccJh/6sREqRjus4HYY8qAeMVWUwlBZwHYWQToGKGgC1iRcAAPIQ6ytqiqoN+Pz3Mozpbo8wDxuu4xDSrBdjXCET8/E+tdbcJAvqDgCoTbnEcRJCOgcqanD90ZNApoDU05/rKA9sy8liWFhg4VAl11EIuSt7GwFmx7jg55Ra/J5ex3WcDkHiroJAbo/aVCpqCGkNXb6oYVkWtUkXIA95FAzPur4c6SUNOHihEk/2doKng4jrOITc04w+zlAqhNj8YxH1I8H1/kYydSTqUi5yHYWQTsG67uJtwFhZCmNVGWSBEVxHeWAfny6BWMDgxRhagZtYB7GQh9mDXHExtx7nMqi1BgDk6u4wlBfDUE4zLxPSUl2+qNFmaAAAtv5hHCd5MDnlehy7XIVpvZxowUpiVSZGOcDNToiPTpVwHaVDaPzdo81K4jgJIdaPippMDRiBEFKvAK6jPJDtP5dCwGcwc4B1ryZOuh6RgIdZA10Qn6XF+UxqrZF6BYDhC1Cfncx1FEKsHhU1mYmw8VGDJxByHeW+ldQY8c3FSkyMcoSL3HpyE9JoUrQjnGwF2PlLKddROMcTiiD19KeWGkJaQZcualiTCdrsFNj6Wdejp71/lMNsYTGzvzPXUQh5KBIhD9N6OeJ0ci2yyvRcx+GcjSoY9dnJ1HmakBbq0kWNLj8DrFEPW79QrqPcN73Rgq/Ol2NwsB28HcVcxyHkoU3r7QQhn8Gnv5VxHYVztr7BMOu00Jfmcx2FEKvWpYua+txUANfXYLEWxy5XoUJrxtN9nbiOQkiLOMuEGN3dHt9cqECNzsx1HE7ZqK7/DqrPSeU4CSHWrUsXNbqCLDBCMcTOHlxHuW97/yhHoKsYvf1lXEchpMWm93GCzsji6OVKrqNwSqL0ARgGDQVZXEchxKp17aImPxNSD5XVTLqXUqTD1XwdJkc70RpPpFMI87BBqLsEX8dXcB2FUzyRGGJndzQUZnMdhRCrZh138zbSUJgFqYcf1zHu24ELlRDwGYzuTitxKw5dOwAAIABJREFUk85jUrQjEgsboCmo5zoKpyTuvtBRUUNIi3TZosakrYGxqgwSKylqDCYLjlysxNBQOzjQZHukExkV6QCxgMHXCV38EZS7CvriXLBmE9dRCLFaXbao0d14dm0tLTVn0+pQWW/GuEccuI5CSKuyk/IxJNQOP1ytgsncdYc0S91VYM0m6EsLuI5CiNXqskVNY4c8iYeK2yD3Ke5KFeykfPQNoA7CpPMZGWGPCq0Zf3ThGYbFSm8AQENxHsdJCLFeXbao0ZcXguELIHLo+ItB6o0WnEyqweOhdhAJuuy3jHRiA9RyyMQ8xF2p4joKZ8ROSgCAobyI4ySEWK8ue4c0lBVB5ORmFSOffk2rhVZvwYgI6iBMOiexkIchoQocv1YDYxd9BCWQO4ARCGGooNW6CXlYHf+O3kYM5UUQ3fhk1NGdSqyBXMKjuWlIpzY0zA41DWZczNFyHYUTDI8HkZMShnIqagh5WF22qNGXF91s7u3ILBYWP6fUon+gHEI+zU1DOq++ATII+QxOJ9dwHYUzIkc36OnxEyEPrUsWNRaDHqaaSqtoqUks1KGszoRBwXZcRyGkTdmK+ejpa4vTybVcR+GM2JlaaghpiS5Z1DR2xBM5d/yi5ueUWjDM9Y6UhHR2g4LtkFGqR25F11y5W6hwgqm2kuaqIeQhPXBRwzCMLcMw/LYI014MVddXBRbZu3Cc5N7+yKhDiFICJxlNuEc6v76B1/uNnc/smv1qBHbX56EyabvuIzhCWuKeRQ3DMDyGYaYzDHOUYZgSAEkAChmGucYwzLsMw6jbPmbrMtVVAwAE8o49mshgsuBSbj16+lEHYdI1BLiI4WjLx/kuOl+NUHb9d5KppusObSekJe6npeYU8P/au+/wtqr7j+Pvo2l57xXbWXZ2CBmQQELZEPYoUMIoo0CBltLSlrbQPWlpadmUQtk/CpSwGnYIJEAZGZA97Ti2470t2Zrn94ccGohJYsfSlXS/r+fxQ2LLV5+I5Pqjc8+5h7HAT4BCrXWp1jofOAL4ALhFKXVRBDMOu0B3+IRhS80wOMnerantxRvQHDI6xegoQkSFUopDRqfyUZUbrc23tNuWFj4n+Xuk1AgxFPtzTeM4rbX/i5/UWrcBzwLPKqXsw54sggI9naAUttTYnny7fHsPSsGMkVJqhHnMGpXCa2s7qWv3U5LtMDpOVNnS+i8/dUupEWIo9llqdhUapdQs4GZgZP/3qfCX9UEDlZ5YFujpxJqchrLE9tSg1bUeRuc6yUyW+TTCPA4uSwZgTZ3HhKWm//KTlBohhmQwPy2fAH4IrAFCkYkTHYHujpi/9ASwrq5X9noSplOen4TdqlhX5+Ekk91F25oULnTBPo/BSYSIT4MpNc1a6xcjliSKAj2dMV9qmrr8NHcHmDwi2egoQkSVw2ZhQlESa2t7jY4SdcpmB4uVkNd8f3YhhsNgSs0vlFIPAIuBz24iobVeOOypIizo7cOWEtv3fVlbF36nNqnYZXASIaJv8ohkXlzVTiiksVjMcydtpRRWZxIhb5/RUYSIS4MpNZcBEwA7/7v8pIG4KzXa78XiyDU6xl5tawr3xvGFSQYnESL6xhck4fGFaOj0U5xlrnk1FqeLoE9GaoQYisGUmmla66kRSxJFIb8XZY/tE+W2pj4K0+2kOGN7MrMQkTAmzwnAtmav+UqNQ0ZqhBiqwdxR+AOl1KSIJYmikM+Lxe40OsZeVTZ7GZMf2xlNIeQNf4ioGpMfHqGsbDbfD3eLXH4SYsgGM1IzD7hEKVVFeE7NZ0u6I5IsgkJ+H5YYHqnRWlPZ4uWrM7ONjmIu2g8dz0PXYvCsgN41oPsLjSUZXNMgeSZkzIf0k0CZcuu0qMhKtpLhslLZbL5CqSwWdCiuF5gKYZjBlJr5EUsRZSGfF4sjdkdB2j1Ben0hSk12jw7DBLuh8c/QfD8EGgZ+TMgD7v+GP5rvAscoyLsW8q8Hi/x/Gm5KKcpyHNS1+4yOYgDzTIwWYrgN5q1mrta6evcPIO5GaQB0wB9eOhmj6jvCJ/LizNjNmDC6FsP6qVD/6y8vNAPxbYe6G2HDDHB/HLF4ZlaU4aChM67u6zmMzLdFhBDDYTCl5h9Kqc8mCiulFgA/Hf5IUaBUTJ8zdnaET+RFGTICEFE7fwlbjgNf9dCP0bcONh4GzfcOWywRVpRpp77TZ749oJQCs/2ZhRgmgyk15wCPKKUmKqWuBK4FTohMrMhSFgtaB42O8aXqO8MjNUUyUhM5dT+B+l8N08GCsONaaLpzmI4nIFzq+/yazt7Y/bcaGYqYftclRAzb71Kjta4Ezie8ieU5wAla685IBYsoZYFQ7J402noC2CyQ4ZLl3BHR8hA03DL8x625HjpfGf7jmlReWnjKX3N3wOAkQoh4sc+JwkqpNXz+bUM2YAU+VEoRj6uflEWhdeyuLmj3BMlItqGUTBgcdr4aqPluhA6uofpKmLwOrLG9DUc8yEwOl/pOj7lKjQ74UDa59CzEUOzP6qdTI54i2pQVYnjJZIcnQFayjNJExI5vQagrcsf310Htj2DkfZF7DpPI6N+dvsNkl59ifXWmELFsf0rNDr2PmXpKKbWvx8SSmB+pcQfJTB7ManuxX/o2QudLkX+e1odhxO/AlhP550pguy6/dnpMVmr8sX9zUCFi1f7MqVmilLpOKVW2+yeVUg6l1DFKqUeASyITLzKUzYH2x+79Lzy+IClOubHbsIvWCiXthZYHo/NcCSw9KVxquvtMVmpkpEaIIdufn5zzgSDwpFJqp1JqvVKqEtgCLAD+qrV+OIIZh53VlUKw1210jC/lDWiS7FJqhl37M1F8rqej91wJymELzynzBeJmEHhY6Bi/47kQsWyf1zi01n3APcA9Sik7kAv0aq07hjuMUmo+cDvhicgPaK0jsEQFrEnJBPtidxfcPn8Ip01KzbDy7QR/ffSer3cNhHxyt+EDYLcqlAJvIHYvFQ83HQqFLz85koyOIkRcGtRPTq21X2tdH6FCYwXuBk4CJgELIrWBptWVQqgvhkdq/Joku6x8GlaeFdF9Pu2DvrXRfc4Eo5TCYVV4TTRSE/R0g9ZYU9ONjiJEXIql4YBDga1a60qttQ/4F3BGJJ7I4nQR7PNE4tDDQqNlOfdw8++M/nP6DHjOBGO3KvxB85SagDu8Ms+WIqVGiKGIpVIzAqjZ7fe1/Z/7HKXUVUqp5Uqp5c3NzUN6olifU6NQ5rs1fKRpA/YQMuI5E0xIgzWWzlIRFugJ38/Ulir3ORJiKPb7dKGUelMpNS2CWQYamtjjJ7vW+n6t9Syt9ay8vLwhPZE1KTm2S41s/TL8LC5zPGeCCYY0FhONWgZ6+kdqpNQIMSSDeQ90I/BXpdRDSqmiCGSpBUp3+30JEJHxe3t6NkFPN6EYXdZtURCSVjO8ksab4zkTTFBGaoQQgzCYvZ9Waq2PAf4DvKqU+oVSajjfin4MVCilRiulHIT3mXpxGI//GXtm+KZo/q62SBz+gDltFlNNjowK13SierXVmg3O0dF7vgSktSYY0lgtJhqp6W4HpNQIMVSDOsur8OzVTcC9wHXAFqXUxcMRRGsdAL4NvAZsAJ7WWq8bjmN/kT2jv9R0tkbi8AfM5bTg8ZlnGWtUWFPANTV6z5cyJ3rPlaB6/RqtMdWNKH1tTVhdKViTko2OIkRcGsycmneBOuCvhCfwXgocBRyqlLp/OMJorV/WWo/TWo/VWv9uOI45kM9Gajpis9QkOyx4vFJqhl3u5Yn5XAnK7Q3fSTjVaZ590HxtTTiyC4yOIUTcGswGQ1cD6wbY4+k6pdSGYcwUcfaMXAD8HS0GJxlYitNCl8k28YuKnEug7iYIRXiSuH0EZEbkbgSmsmt7hFQzjdS0N+HIzjc6hhBxazBzatbuZdPKU4YpT1TYUjPAYsXfGZtzatKTrHSYbBO/qLBmQMEPIv88RT8HJRuSHqievvBoZWqSjNQIIfbPsLwF0lpXDsdxokVZLDgyc/G1NRodZUA5qXZaewJGx0hMRTeDK4J3Jkg7HvKuitzxTaTVHf43kJ1ijoIY7PMQdHfJSI0QB8A847pf4Mwfgbe5zugYA8pNteHxheiVycLDT9lh9GNgSRv+Y9sKYNQDw39ck2ruCt+8MD/dbnCS6PC1NgDgyCk0OIkQ8cvUpaavsdboGAPKSQ2/M5XRmghxTYXyF8GSMnzHtOVCxevgKBu+Y5pcc7cfpSDHJCM1fQ07AEgqlL9DQgyVaUtNUn4JQXfXZ3utxJKC/nemDZ2xeXPAhJB2FIxbDI6RB36spEkwbikkH3TgxxKfaeoOkJ1iw2Y1x31q+urDpcZZUGJwEiHil2lLjbMgfPNibwyO1pRkOwCoaZdSE1Eps2HSGsj95hAPYIXCH8PEleCaOKzRBOxs91GUYY5LTxAeqXHkFGB1yvYaQgyVaUtNUn54r8xYvARVlGHHoqC2TUpNxFnTYOR9MPFTyL0aLKn78T3ZkH8DTN4II/4AFmfkc5pQdauXkTnmeW37GnbIpSchDpA5LlYPwJFXDMpCX1PNvh8cZQ6bhYJ0O7UyUhM9yQfByHuh5FZwfwSe5dC7GoKdgAJrFiRPh+SZkHIIWJKMTpzQfIEQ9Z1+Ts9xGB0lKnQoRF/DDnIr4uruGELEHNOWGovNHp4sXFdldJQBjc51UtXsNTqG+VhTIf2Y8IcwTG27j5DGNCM1vpZ6Qr4+koplvzAhDoRpLz8BJJeW46nZYnSMAZUXJLG1qY9gSDa2FOazrSlc6EflmqPUeHaEz0PJZRUGJxEivpm71JSNw9faGJMroMYVJOENaJlXI0xpY30vFgUVBea4zOfesQllteGSkRohDojJS034XdGud0mxpLz/ZL65sc/gJEJE38aGXkbnOUmym+MU5anejGvEGCx2c8whEiJSzHHG+BLJpeUAMXkJqqIgCZsF1tV5jI4iRNRtrO9jQqE5ljZrrfFUbyZ55DijowgR90xdamxpmdiz8undsdXoKHtIslsYX+Ti0xopNcJcmrv9NHT6mVRsjlLja20g6OkmuUxKjRAHytSlBsKXoNzVm4yOMaCDSpJZU9crk4WFqazaES7yB5clG5wkOtzb1gGQMmqCwUmEiH+mLzWpYyfjbazB39VudJQ9TCtNptcXYovMqxEmsqrajcOmTDNS073lUyxJKbhKxxodRYi4J6WmYhoAPVvXGJxkTzNHhjdcXL7dbXASIaJn1Q43U0a4cNjMcXrq2bKG1PIpKIvV6ChCxD1znDX2InnkOJTdSc+WT42OsofiLAelWQ4+2NZjdBQhoqKnL8iGnb2fFfpE5+/uoK9+O2kVshmqEMPB9KXGYrOTOmYSPVtib6QGYPbYVD6u6iEQlHk1IvF9WNlDIASHl6cZHSUqdo0Qp0qpEWJYmL7UQPiE4qnZSrA39i7zzBmbSo83xLqdsgpKJL73tvaQ7LCYZpJwz+ZPUHYHySPHGx1FiIQgpYb+d0k6RM/WtUZH2cNhY1OxKHhnU7fRUYSIuPe3dHPomBTTzKfpWr+c1PKpctM9IYaJOc4c+5AyZhLKZqd74wqjo+whM9nGjJEpLNkYe1s5CDGctjX1UdPuY16FOS49eVsb6KuvJmPKbKOjCJEwpNQAVqeL1IppdK790OgoAzp6YjqbG/pkHyiR0N5Y1wnAsRMzDE4SHV1rPwIgXUqNEMNGSk2/jKmz6auvxttSb3SUPRw9IR2AtzZ2GpxEiMh5Y30nB5cmk59uNzpKVHSu/RBHdgFJhWVGRxEiYUip6ZcxdQ5ATI7WjMxxMqEoiVdWS6kRiammzcvG+j6On2yOUZpQwE/3xpWkTzkUpZTRcYRIGFJq+iUVlOLMH0HXmg+MjjKgUw7KZHWthx2tXqOjCDHsXl7dAWCaUtOzZTUhb6/MpxFimEmp2U3GlDl0bVxFyBd7xeHkgzJR6n8nfyEShdaaF1a1M2tUCiOyzLEKqGPlUiyOJNInzjI6ihAJRUrNbtKnzkb7vXRvXGl0lD0UZjiYNTKFFz9pR2u5EZ9IHJ/WeKhu9XHm9Cyjo0SFDgVpX7WU9CmzsTiTjI4jREKRUrObtHEHY3Wl0L7ibaOjDOjsmdlUt/r4uCr2bhIoxFC9sKodl11xwhSTXHraupZAVztZM480OooQCUdKzW4sdgeZ04+g45N3Cfljb/n0CVMySE+y8szyNqOjCDEs3N4gi1Z3cPzkTFKc5tjQsX3lOyi747PFCUKI4SOl5guyZh1NsNdN17qPjY6yhyS7hdOnZ/LGuk7a3QGj4whxwF78pB23N8SCOTlGR4kKHQrRsXIpGZMPxZpkjq0ghIgmKTVfkD5hJtaUdNqXLzE6yoDOnZWDP6h55mMZrRHxTWvN/33QypQRLg4qMccP+J4tq/F3tJA18yijowiRkKTUfIGy2ciafgQdn74Xk6ugyguSOHxsKk980IIvEDI6jhBD9mFlD5XNXhbMNscoDUDrf1/D4nSRcfBco6MIkZCk1Awga9bRhLy9dMboPWsumZdHS09AlneLuPbgsmZyUmycNDXT6ChREezz0L5iCVmzjsbqdBkdR4iEJKVmAGnjD8aWnk3rf18zOsqA5panUp7v5JH3WmR5t4hL6+o8vL+1h4sPz8VpN8dpqGPVMkLePnIOn290FCESljnOJoOkrDZyD59P55oP8LU3Gx1nD0opLj8in82NfbJ7t4hLDyxtJtVp4XwzXXp6/1UcucWklk81OooQCUtKzZfImXcy6BCt779qdJQBnXJQJqXZDu5+q1FGa0Rc2dbUxxvrO1kwO4e0JHMs4/a21NO9aRW5h8+XvZ6EiCApNV8iKb+EtPHTaXnvZXQo9ibk2qyKq4/KZ2O9jNaI+HLX4kZcdgtfn5tndJSoaVm2CJSF7MNOMDqKEAlNSs1e5B5xKr7+d1ix6NRpWZTlOLhzcSPBkIzWiNi3rs7D6+s6uXRuLtkpNqPjREXI76Vl2UtkHHQYzpxCo+MIkdCk1OxF5vR5WFPSaVn2H6OjDMhmVXznuEI2N/Tx0iftRscRYp/+9kYDmclWLjHRKE378ncI9HSSf/RZRkcRIuFJqdkLi91Jzpzj6Vi1DH9Hq9FxBjR/SgZTS1zc8WYjvb7Yu0wmxC7vb+3m/a09XHlkPqkmmUsD0LTkOZwFpaRNnGl0FCESnpSafcg7+ix0KEjTO88bHWVASil+OL+Yxi4/j74feyu1hADwBzV/WLSTshwHF5hoxZN7+0Y82zeQf/RZMkFYiCiQUrMPSfklZE6bS/PbLxDy9hkdZ0AzR6Vw3KR0/vFOEzs7Ym8jTiGe/LCFymYvN55UjMNmntNO01sLsThd5Bx2otFRhDAF85xdDkD+8ecRdHfR+sHrRkf5Uj86qRgN/PHlnUZHEeJzWnsC3PNWI3MrUjlqfJrRcaLG29pA20eLyZ13ClZXitFxhDAFKTX7IbV8Kskjx9P45jMxubwboDjLwdVHFfDm+i6WbpIl3iJ23PLyTvr8mh+fXGyqSzCNrz+FUoqCE84zOooQpiGlZj8opSg4/jy8jTV0ro3N/aAALpmby5g8J799qQ63N2h0HCF4Z1MXL6/u4Koj8xmTl2R0nKjxd7XT8u4isuecgCMr3+g4QpiGlJr9lDXjSOxZeTS+9i+jo3wph83CL88oYWenn9tebzA6jjA5tzfIr1+oozzfyRVfMc8SboCmt55FB/wUnni+0VGEMBUpNftJ2WwUnriAni2r6d4Ymzfjg/Ck4YsOy+VfH7byYWWP0XGEid36Sj2N3X5+fVaJqSYHB3vdNL/9PJnTjyCpsMzoOEKYinnONMMg94hTsGfmsvOlh2J6v6XrjyukLMfBzxbW0NMnl6FE9C3Z2MUzy9u4dG4e00rNNUm26a1nCXp6KJx/gdFRhDAdKTWDYLE7KZx/QXi0Jka3TgBwOSz8/qul1Hf6+c1LdTFdwETiaenx8/PnaphQlMR3jiswOk5UBdxdNLz2FJkHzyNl1ASj4whhOlJqBmnXaE39i7E9WjO9LIVrjyngP5928OInHUbHESYRCml+trAWtzfEH88tM9VlJ4CG154k5PVQfMblRkcRwpTMdcYZBha7k8KTLqRn6xq6N6wwOs5eXXVkPoeMSuE3L9WxvcVrdBxhAg8ua2bp5m5+ML+I8nzzrHYC8He20rR4IdmHHotrxBij4whhSlJqhiB33snhuTUv/jOmR2usFsUt55bitCmu/7/tssxbRNSHlT3c8WYDJ03NYIGJtkLYpf7lx9FBP0WnXWp0FCFMS0rNEFjsTopPuxR35Xo6VrxjdJy9Ksxw8Ofzyqhs9vLThbUxXcJE/Grq8vPDp3YwKtfJr84sMdVN9gC8zTtpWfoSuXNPISm/xOg4QpiWlJohypl7Eq4RY6hd+HdC/ti+tHNYeRrfO6GQ19d18uAy2fRSDK9eX4hvP7Edjz/EXxeMJMVpnh24d6n9970om42iU79udBQhTC0mSo1S6lyl1DqlVEgpNcvoPPtDWayUnHstvpZ6mt56zug4+3TZvDzmT8ngb2808NaGTqPjiAQRCml+urCG9Tt7+dO5paabRwPQtXElHauWUXjSRTiyzHWTQSFiTUyUGmAtcDaw1Oggg5E+aRYZU+dQv+gx/N2xvcJIKcVvzy5lcrGLG5/ewdo6j9GRRAK4Z0kjr67t5IYTCjlmYobRcaJOBwPUPnUXjpxCCo4/1+g4QpheTJQarfUGrfUmo3MMxYhzriHk66X+pYeNjrJPLoeFuy8eRXaKjWsf287Odp/RkUQcW7iijXuXNHHWjCwum2fOEYqWdxfRW1dJyTnXYLE7jY4jhOnFRKmJZ66ikeR95XSal75Ib12l0XH2KTfVzr1fH40/oLnqkSpaewJGRxJxaPH6Tn7xfC1zK1L5+ekjTDcxGCDg7qbu+QdJHTeNzBlfMTqOEIIolhql1JtKqbUDfJwxyONcpZRarpRa3twcG5Nei0+/DFtyGtWP34YOhYyOs09j85O466JR1Hf6uOqRSrp6Zam32H/Lq3r4wdM7mDIimb8tGGW6G+ztUvfsfQR7eyj92nWmLHVCxKKonY201sdpracM8PHCII9zv9Z6ltZ6Vl5ebAx521IzKDn3Wtzb1tLy7iKj4+yXmaNSuP2CUWxt8nLNY1V4fLFfxoTxPq1x863Ht1OS5eDer48i2WHOQtO96RNa3l1EwfFfI7m03Og4Qoh+5jwjRUD2nBNIGz+duoV/x9/ZanSc/TKvIo1bzytjdY2Hax+tkpvzib1aXePhqoeryEqx8Y9LR5OZbDM6kiFCfh/Vj/8FR24RxadeYnQcIcRuYqLUKKXOUkrVAocBi5RSrxmdabCUUpRd+D1CPi81z9xjdJz9dsLkDG45p5SVO9xcLcVGfInVtR6ufLiSrBQbD39jDIUZDqMjGabhlSfwNtZQduH3sDjNt4RdiFgWE6VGa/2c1rpEa+3UWhdorU80OtNQJBWWUXjShbR/tJjOdR8ZHWe/nTIt67MRmysfrpI5NuJzPqrs4YqHKslKtvHQ5eYuNL311TS88gTZs48jY/KhRscRQnxBTJSaRFI4/wKcBaXsePwvBHvdRsfZbydOyeQv549k3c5eLnlgG41dfqMjiRiweH0n33y0isIMOw9fMYaiTPMWGh0IsP2hP2BJSqbk3G8ZHUcIMQApNcPMYncw6tIf42trpuapu4yOMyjHTcrgvotHUdvu48K/b6Wyuc/oSMJAC1e08d0nq5lQ6OLRK8aaeoQGoP6Vx/Fs38jIi27Anp5ldBwhxACk1ERA6tjJFM6/gNb3X6Hjk3eNjjMoh5Wn8cgVY/AFNBfdv40V2+NntEkMj1BIc/sbDfzsuVrmjE3lgcvMOyl4F3fVBuoXPUr27OPJmnmU0XGEEF9CSk2EFJ12Ca7SCqof+zP+rnaj4wzKpOJkHr9qLJnJVi5/qJKFK9qMjiSipNcX4vtP7eD+d5r46sxs7rl4tCk3qNxdyNtH1T9/jz0jh9IF3zE6jhBiL6TURIjFZmf0N24i2Oum+rE/o7U2OtKglOU4efLqcg4ZlcLPnqvlT6/sJBiKrz+DGJz6Dh+XPLiNN9Z38sP5RfzqzBHYrXJTudqFf8fbWMOoS3+MLTnN6DhCiL2QUhNBruLRjDjrSjo/fY/W9142Os6gZbhs3Pf10Vw4J4dH3mvhqoeraOmRCcSJ6L9buzn3ni1UNXu544KRXDovT+6SC3R88i7NS54j/9ivkj5xptFxhBD7IKUmwvKPPYe08dOp+dedcbE31BfZrIqbTh3Bb84qYdUON+fcvUXm2SSQUEjz97cbufKRKnJSbTx9Tbkpd9seiLelnu0P30LyyHGMOPubRscRQuwHKTURpiwWRn/jp1iSXGy77xcE+zxGRxqSs2dm8+TV5SQ7LFz2z23c/3aTXI6Kc01dfr75SBV3vNnIyVMzefLqCkbnyc3kIHzX4Mr7fwVaM+aqX2Kxm3vllxDxQkpNFNgzcxhz5c/xNtXF5fyaXcYXunj6mgpOmJzB7W82cOmD26ht8xkdSwzBG+s6OfPOzaza4eYXZ4zgj+eWmnYfp4HUPntfePn2pT/GmVdsdBwhxH6Ss1iUpI2fTvEZl9P+8Vs0vzOoPTxjSmqSlVvPK+OWc0rZ3NDHWXdtZuGKtrgtambT2Rvg5mdr+O6T1ZRkOfj3t8Zx3iE5Mn9mN+0r3qb5rYXkH3sOWdOPMDqOEGIQpNREUeH8C0ifMofap+/GvX2j0XGGTCnFaQdn8dx145g8wsXPnqvlqoerqGnzGh1N7MXr6zo57fbNvPRpO1cdmc/jV41lVK7T6Fgxpbe2ku0P30LK6ImM+KrMoxEi3qh4foc9a9YsvXz5cqNjDEqgp5MNv70KgAk33Ys9PdvgRAcmFNI8/XEbt71eTzCkufaYAr5+eJ4sBY4h9R0+bnnztBzoAAAZZElEQVR5J2+u72JiURK/OauUicUuo2PFHH93Bxt/fzU64GfCTffhyMozOtKQKKVWaK1nGZ1DCCPISE2U2VIzGHPNr/F3d7Dtnp8R8sf36IbFojh/dg4vXT+eueVp3PZaA1+9ezP/3dptdDTT8/pD/P3tRk67fRPLNnfzvRMK+dfVFVJoBhAK+Km87xf4O1sZe+1v47bQCGF2UmoMkDJyPKMv/wnuynVUP3prQsxHKUi3c8eFo7jzwpH4/JorHq7i249vp7o1vktbPNJas3h9J2fcsZk73mxk3rh0Xrp+PFd8JR+bjKDtQWtNzZN30LPlU0ZeciMpoycaHUkIMUTm3tDFQFkzj6L4jG+w84UHSSocSdEpFxsdaVgcMzGDueVpPPp+C39/p4nT79jMubOy+eZR+eSl2Y2Ol/A+rurhb6838EmNh7H5Th64dDSHlctdcPemeclztCx7iYL5C8iZfbzRcYQQB0BKjYEKT76Ivobq/mJTRtbMI42ONCycdgtXHpnPmTOyuHdJI8983MpzK9u4cE4ulx+RZ/rNESNhbZ2HO99s5N0t3RSk2/nVmSM4c3q2jMzsQ/uqZdQ8dRcZBx3OiDOvMDqOEOIAyURhg4X8Xjb/5QY8NVsZ//2/kjJmktGRht2OVi/3vNXIf1Z3kGS3cN4h2Xz98FwKM+SGZgdCa83HVW7+8U4T72/rId1l5aoj81kwO4cku1xZ3peeLavZ/Nfvk1xaTsUNt2F1JsZcI5koLMxMSk0M8He1sfGWbxPs7WH8jXfiKhppdKSI2NrYxwNLm3h5TQdKKU6flsnX5+ZRUSB3sR2MQFCzeEMnj77fwic7POSk2rhkbi5fOySH1CRz76i9v3p3VrHpT9dhS8tiwo13YkvLNDrSsJFSI8xMSk2M8DbvZOMfv42yWBj/o7tw5hQaHSli6tp9PPRuMwtXtOENaA4ZncKC2TkcMzFDloLvRZs7wL8/buOpj1pp6PIzItPOZUfkcdaMbBmZGQRfexMbb/kWOhhkwo/vxplbZHSkYSWlRpiZlJoY4qndxuZbv4MtPZvxN96JPYHePQ6k3R1g4co2nvqwlboOP/lpNk6fnsUZ07MYI3sQAeFRmfe3dvPcynaWbOzCH9QcNjaVC+bkcOT4dKwWKYGDEejuYNNfvouvrYnxP7yd5NIKoyMNOyk1wsyk1MSYni2r2fy3H+IqHsm4G/6K1ZVidKSIC4Y0Szd38/RHrby3tZtgCKaWuDj94CyOn5xhulVTWmvW1vXy2toO/vNpB83dAbKSrZwyLYtzD8mmPF8K31AE3F1svu0G+hp2UHHdH0mbMN3oSBEhpUaYmZSaGNS55gO23n0zqeVTKL/uloSZwLg/mrv9LPq0gxdWtbO5sQ+lYHpZMsdPzuC4iRkUZyXm5OJgSLO61sOb6zp5fV0nOzv82CwwtyKNs2Zkc+T4NBw2ucQ0VAFPN1v++n1666oY+63fkTH5UKMjRYyUGmFmUmpiVNtHi6l68HekVhxE+XV/MFWxgfBoxbYmL6+v6+SN9Z1sbugDYEyek7kVacwtT2PWqBRccbyzdFOXn/e3dvPulm7e29pDV28Qm1Vx+NhUTpiSwdET0mX5+zAI9rrZcvsP8VRvZsw1vybzoMONjhRRUmqEmUmpiWFtH75J1T9//78Rm6RkoyMZprrFyzubu3h3SzcfV7nxBTQ2q2JysYsZI5OZMTKFaaUp5KTGZgkIhTTbW718ssPDiu1uVmx3U9PuAyA31ca8ijTmjQuXtXSXrGAaLsE+D1vv+BE9lesZ+81fkmmCXbel1Agzk1IT49o+fouqB39LyuhJVFz/J1MXm116fSGWb3fzcVUPK6vdrK3rxR8M/z3OT7MxsdjFhCIX4wuTGJnjpCzHSfJeRnRCIc27W7pZv7OXScUu5lWkYRniBFytNa3uANUtPqpa+tjU0MeGnb1sbOij1xcCIDPZysyRKcwYlcLs0amML0wa8vOJLxdwd7P1zh/jrtrAmCt/Ttaso4yOFBVSaoSZSamJA+3L36bygV+TMmpiuNiYYPLwYHj9IdbW9bKmzsPGnb1sqO+lstlLaLe/2gXpdkqyHOSn2yhIt5Ofbic31UZqkpUHlzaxotrz2WPnVaTy56+V4bBaCIQ0/mD/R0DT7Q3S3Rukqy9IZ2+Q5q4Ajd1+mrr8NHb62dHmpbsv9NmxXA4LE4uSmFjkYmKxi6klyYzNc6KUlJhI8ne1seVvP6SvYQejr/w5WSYYodlFSo0wMyk1caJ95VIq7/8VyWUVVFx3S0LdLCwS+vwhqlq87Gj1sr3FS3Wrj7p2H41d4QLiDQzf3/t0l7W/KNkozXYyOtfJyBwno3IdjMh0yChMlPlaG9n8tx/gb29m7LW/JX2SuX6+S6kRZialJo50fPo+lff/Ekd2ARXX/ynhbhoWLVprOnuDtPQEeOz9Fv69vG2Px8wtT+WQ0anYrAq7VWGzhP+bmmQhw2UlLclKustKbqo9ricrJ5q+xho23/Z9Qn1uyq+7hdTyqUZHijopNcLMpNTEmZ6ta9h6100om52K6/9Ecmm50ZHi2tJNXVzz2PY9Pn/vxaP4yvj06AcSQ+auXM/Wu28GNBXX30pyWeLdWG9/SKkRZiZvMeNMavlUxt94J8pqZdOt19O9cZXRkeLavIo0jp7w+fJy9IR05lWkGZRIDEX7yqVs+st3sThdjP/hHaYtNEKYnYzUxClfexNbbr8Rb1Mdoy6/iexZRxsdKW7tWv20ob6XiUUHtvpJRJfWmqY3nqb22ftIGT2Rsd/6fcJvL7IvMlIjzExKTRwLuLvZdvdN9GxdQ9Hpl1F08sUoiwy+CXPQwQA1/7qT5ndeIHPGkYy+/CYsDqfRsQwnpUaYmfwEjGO2lDQqvvdnsg87kfoXH6Ly/l8S9PYaHUuIiAu4u9l61000v/MCBScuYMxVv5BCI4QgNm+/Kvabxe5k1KU/JrlkLLX/vo9NTd9m7Ld+hzOn0OhoQkSEp3Yb2+75Gf72Jsou/gF5R5xqdCQhRIyQkZoEoJSi4PjzKL/uD/haG9j4u2/SvflTo2MJMezaPlrMxj9ciw74GPeD26XQCCE+R0pNAsmYMpsJP7kXa0o6m2+7gcY3niGe50wJsYsOBqh5+m6qHvgNKSPHMfHmv5M6drLRsYQQMUZKTYJJKixjwk/uIWPqHGqfuZttd99MwN1ldCwhhszX1sTm275P05vPkHf0WVTccBv2jByjYwkhYpCUmgRkS05j7LW/peRr36Zr3Ues//UV9Gxba3QsIQatfdUy1v/6G3h2bGLUZTdRtuB6LDa70bGEEDFKSk2CUkpRcOw5jP/RXSirjU23foeGV/8PHQrt+5uFMFjI20f1E7dRee/PcOYVMfGnD5Bz2AlGxxJCxDhZ/ZTgUkZNYNJP76f6sT9Tt/B+ujasZNQlN+LIzjc6mhAD6q2tpPIfv6avfjsFJ3yN4jOvkNEZIcR+kZvvmYTWmpZl/6H2mbtBWSn92rfJOXw+Ssmdc0Vs0MEADa8/Rf1LD2NNTmX0ZTeRPvkQo2PFHbn5njAzGakxCaUUeV85jfSJM9n+8B+pfuSPdKx8h7KLf4AjM9foeMLkemsr2f7ILXiqN5M540jKLvgu9vQso2MJIeKMjNSYkA6FaFqykLqF/8Bit1N6/nfInn28jNqIqNOBAA2v/h/1ix7FmpxK2QXfJWvmUUbHimsyUiPMTEZqTEhZLBQcew4ZU2az/aFb2P7P39P20VuULfgOzrxio+MJk3Bv30j1Y3+mt2YrWYccQ9n538Fm8s0ohRAHRkZqTE6HgjS9tZCdL/wTHQxQePJFFJ54Pha77KMjIiPQ3UHd8w/Q8u4ibGlZjLzwe2ROP8LoWAlDRmqEmUmpEQD42pupfeYe2pcvwZk/gtIF15Mx+VCjY4kEokNBWt59mbrn/kGwt4f8Y8+h+NRLsLpSjI6WUKTUCDOTUiM+p2v9cnY8eTvexhoyZxxJyTlX48wtMjqWiHPuqg3sePJ2PNs3kloxjbILrsc1YozRsRKSlBphZlJqxB5Cfh+Nrz9F/cuPgw6Rf/RZFJ58EbaUdKOjiTjjbaqj7vkHaF++BFt6NqXnXkvWocfKpPQIklIjzExKjfhSvvYmdr7wEK3/fRWrK4XCky4i/5izZL6N2Cd/Vzv1ix6l+Z0Xsdjs5B9/LoUnnC+XmqJASo0wMyk1Yp96ayupXXg/XWs/wJFdQPEZl5M9+ziUxWp0NBFjgn0emt78Nw2vPUnI7yV33ikUn3op9kzZgDJapNQIM5NSI/Zb18aV1P37Pjw7NuMsKKXo5IvIPvRYlFXuDGB2QU8PTW8/T+MbTxN0d5E5/QhGnHkFSUUjjY5mOlJqhJlJqRGDokMhOlYupf7lx+it3YYjt5iiky4g+7ATZX8eEwq4u2ha/CxNbz1L0NND+pQ5FJ16MaljJhsdzbSk1Agzk1IjhkSHQnSufp/6RY/hqd6EPSufwvkLyJ17MhaHzLlJdP7O1nCZeft5Qn0eMg+eR+EpF5MycrzR0UxPSo0wMyk14oBorela9xH1ix7DvW0t1pR08r5yGnlHnYEjS3YCTzTu7RtpWvws7cuXoENBsmYeRdHJF+MqkeXZsUJKjTAzKTViWGit6dn8KU1vPUvHJ++BgqwZR5J/zNmkjJ0iS3jjmA4GaF+1jKbFz+LethaL00XO3JPIP/pskgpKjI4nvkBKjTAzmeEphoVSirTxB5M2/mC8LfU0v/08Le8uon35EpLLxpF31JlkzToKa1Ky0VHFfvK2NtD639doWbYIf3sTjtxiSs77FrmHn4Q1OdXoeEIIsYeYGKlRSt0KnAb4gG3AZVrrjn19n4zUxLagt5e2D96g6a2F9NVvx+JMImvGkeQcPp/Uimkoi8XoiOILQn4vHavepeW9V+jeuAK0Jm3iTPKPPpuMg+bIMv44ICM1wsxipdScALyltQ4opf4IoLX+0b6+T0pNfNBa465cR+v7r9L28VuE+jw4covIOexEcg47UbZhMJgOhXBXbaDtwzdo++hNgp4eHDkF5Bx+kvz/iUNSaoSZxUSp2Z1S6izgHK31hft6rJSa+BPy9tH+yTJa33+V7o0rQWuSR00ka+aRZM34Cs68YqMjmsKuItO+fAntK9/B396MstnDI2lzTyJt/HQZSYtTUmqEmcViqXkJeEpr/fi+HiulJr75Whtp++hN2le+g6d6MwCu0nKyZnyFzBlH4pIbtw0rHQjQU7mOjlXLPldk0icdQtbMo8icdrjMlUkAUmqEmUWt1Cil3gQKB/jSzVrrF/ofczMwCzhbf0kwpdRVwFUAZWVlM6urqyOUWESTt6W+/4ftUtzb1gLgzC8hfdIs0icdQtqE6TLJeAh87c10rfuIzrUf0r1hBcFe9/+KzKyjyDxIikyikVIjzCxmRmqUUpcAVwPHaq09+/M9MlKTmHwdLXSsWkbX2g/p3vQJIV8fWKykjp3SX3JmkVxWIdszDCDg7sa9bS3dW1bTte4jemu3AWDPzCV9yqFkTJ5N+sSZUmQSmJQaYWYxUWqUUvOB24AjtdbN+/t9UmoSX8jvw71tHV3rP6Zr/XI8O8KXqSyOJJJHTSB17GRSxk4mdfQkbGmZBqeNPl97Ez1b1tCzdQ09W1bTu7MKtA6XwPIpZEyZTfqU2bhGjJF7BZmElBphZrFSarYCTqC1/1MfaK2v3tf3SakxH39XO92bVuHeto6ebWvx1GyBUAgAZ0EpKaMn4ioZg2vEGJJHjMWWkZ0QP8x1KISvpR5PzRY8O7biqdlKb80W/J3hfzIWp4vUsVNILZ9KasVUUkZNxOJMMji1MIKUGmFmMVFqhkpKjQh5+3BXb8K9bS09levxbN/42Q96AGtKerjglIzFWViKM6cIZ14RjpwCLPYv36NKh0J0rfsIz47NJJeNI33yoVFZDRTs8+BtrsPbWEdfcx3exlr6Gmvorasi1OcOP8hixVU8CldpOSll40gtn4qrZIxcjhOAlBphbnIWFHHN4kwibdw00sZN++xzgZ5Oeusq6a2tDP+3rpKWdxeF5+bsxp6ZizO3CEduEfb0bGxpmdjTMrGmptP42r/o2bL6s8dmTJvL2Gt+M6RiE/L7CPV5CPZ5CPa68Xe27vbRhr+jBX9nK97WRgJdbZ/PmJGDM7+EnDnHk1xWgau0HFfxqL0WMiGEMCsZqRGmoEMh/F1t+Frq8TbX422t/+zXvtYG/F1t6IB/r8ewuFKwuVKxOJxYHEkoqxWtNegQOhT+b/i5goR8XoK9bkLe3r0e15qchj0jB3tGNo7sfJz5JSTll+AsGIEzb4Ss+BKDJiM1wsxkpEaYgrJYcGTm4sjMJbV86h5f11oT8vYS6O6g4dX/o2XZf/Z4TFJBKUlFIwl5+9B+L6FAIDxyoywoiwr/VylQCosjCWtSMhZXMtakFKxJyf2/T8Geno09I/whIy5CCDF8pNQIQXhDzl3FI/PgeQOWmuLTLiVj6hwD0gkhhNgfch90Ib4gffKhZEyb+7nPZUybS/rkQw1KJIQQYn/ISI0QX6AsFsZe85v+1U9bSC6riNrqJyGEEEMnpUaIASiLhYypc+RykxBCxBF56ymEEEKIhCClRgghhBAJQUqNEEIIIRKClBohhBBCJAQpNUIIIYRICFJqhBBCCJEQpNQIIYQQIiFIqRFCCCFEQpBSI4QQQoiEIKVGCCGEEAlBaa2NzjBkSqlmoDqCT5ELtETw+PFAXgN5DUBeA4if12Ck1jrP6BBCGCGuS02kKaWWa61nGZ3DSPIayGsA8hqAvAZCxAO5/CSEEEKIhCClRgghhBAJQUrN3t1vdIAYIK+BvAYgrwHIayBEzJM5NUIIIYRICDJSI4QQQoiEIKVmL5RStyqlNiqlViulnlNKZRqdyQhKqXOVUuuUUiGllGlWfyil5iulNimltiqlfmx0HiMopf6plGpSSq01OotRlFKlSqklSqkN/f8Orjc6kxBiYFJq9u4NYIrW+iBgM/ATg/MYZS1wNrDU6CDRopSyAncDJwGTgAVKqUnGpjLEw8B8o0MYLAB8X2s9EZgDfMukfxeEiHlSavZCa/261jrQ/9sPgBIj8xhFa71Ba73J6BxRdiiwVWtdqbX2Af8CzjA4U9RprZcCbUbnMJLWul5rvbL/193ABmCEsamEEAORUrP/LgdeMTqEiJoRQM1uv69FfpCZnlJqFDAd+NDYJEKIgdiMDmA0pdSbQOEAX7pZa/1C/2NuJjwE/UQ0s0XT/rwOJqMG+JwsFTQxpVQq8CzwXa11l9F5hBB7Mn2p0Voft7evK6UuAU4FjtUJvP59X6+DCdUCpbv9vgTYaVAWYTCllJ1woXlCa73Q6DxCiIHJ5ae9UErNB34EnK619hidR0TVx0CFUmq0UsoBnA+8aHAmYQCllAIeBDZorW8zOo8Q4stJqdm7u4A04A2l1CdKqfuMDmQEpdRZSqla4DBgkVLqNaMzRVr/BPFvA68Rnhj6tNZ6nbGpok8p9STwX2C8UqpWKfUNozMZYC5wMXBM/3ngE6XUyUaHEkLsSe4oLIQQQoiEICM1QgghhEgIUmqEEEIIkRCk1AghhBAiIUipEUIIIURCkFIjhBAxYDCbhyqlyvo32VzVv+GurMYSAik1QggRKx5m/zcP/Snh2wxMJ3wPpXsiFUqIeCKlRgghYsBAm4cqpcYqpV5VSq1QSi1TSk3Y9XAgvf/XGcjdroUApNQIMShKKZdS6h2llFUpNWp/LhV8yXEcSqmlSinTb1Ui9up+4Dqt9UzgB/xvROaXwEX9N8V8GbjOmHhCxBYpNUIMzuXAQq118EAOorX2AYuBrw1LKpFw+jfQPBx4Rin1CfB3oKj/ywuAh7XWJcDJwGNKKTmfC9OTfwRC9OufeHl8/69/q5S6Y4CHXQjssWu5UmpM/6TNQ/pHcDYqpR5QSq1VSj2hlDpOKfWeUmqLUurQ/m97vv94QgzEAnRorQ/e7WNi/9e+ATwNoLX+L5AE5BqUU4iYIaVGiP/5BXCzUupCYDrwvd2/2L+x5Rit9fYvfH484R2cL9Naf9z/6XLgduAgYAJwATCP8CWEm/ofsxY4JCJ/EhH3tNZdQJVS6lwIb6yplJrW/+UdwLH9n59IuNQ0GxJUiBgipUaIfv0TNRVwA3D+AJeYcoGOL3wuj/DIzUVa6092+3yV1nqN1joErAMW6/BGa2uAUf3PFwR8Sqm0Yf/DiLjzJZuHXgh8Qyn1KeG/R2f0P/z7wJX9n38SuFTLRn5CIJMUheinlJpKeM5Ci9a6e4CH9BJ+R7y7TqCG8E7Ou+/i7d3t16Hdfh/i8//unEDfAcQWCUJrveBLvrTHMm+t9XrCf+eEELuRkRohAKVUEfAE4XfCbqXUiV98jNa6HbAqpXYvNj7gTODrSqkLBvmcOUCz1to/9ORCCCF2kVIjTE8plQwsBL6vtd4A/IbwktmBvE54bsxntNZu4FTge0qpMwb8roEdTXg5rhBCiGGg5DKsEPtPKTUduEFrffEwHGsh8BOt9aYDTyaEEEJGaoQYBK31KmCJUsp6IMfpX0n1vBQaIYQYPjJSI4QQQoiEICM1QgghhEgIUmqEEEIIkRCk1AghhBAiIUipEUIIIURCkFIjhBBCiIQgpUYIIYQQCUFKjRBCCCESwv8DD0oUXFGKcEUAAAAASUVORK5CYII=\n",
"text/plain": [
"<Figure size 1008x432 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"plotter = StaticOrbitPlotter(plane=Planes.EARTH_ECLIPTIC)\n",
"plotter.plot_body_orbit(Earth, EPOCH, label=\"Earth\")\n",
"plotter.plot_body_orbit(Mars, EPOCH, label=\"Mars\")"
]
},
{
"cell_type": "code",
"execution_count": 6,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"Ephemerides at 1 epochs from 2020-10-06 00:00:00.000 (TDB) to 2020-10-06 00:00:00.000 (TDB)"
]
},
"execution_count": 6,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"earth = Ephem.from_body(Earth, EPOCH, plane=Planes.EARTH_ECLIPTIC)\n",
"earth"
]
},
{
"cell_type": "code",
"execution_count": 7,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"Ephemerides at 1 epochs from 2020-10-06 00:00:00.000 (TDB) to 2020-10-06 00:00:00.000 (TDB)"
]
},
"execution_count": 7,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"mars = Ephem.from_body(Mars, EPOCH, plane=Planes.EARTH_ECLIPTIC)\n",
"mars"
]
},
{
"cell_type": "code",
"execution_count": 8,
"metadata": {},
"outputs": [
{
"data": {
"text/latex": [
"$62074610 \\; \\mathrm{km}$"
],
"text/plain": [
"<Quantity 62074610.23119402 km>"
]
},
"execution_count": 8,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"distance = norm(mars.rv(EPOCH)[0] - earth.rv(EPOCH)[0])\n",
"distance.to(u.km)"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": []
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3",
"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.8.5"
}
},
"nbformat": 4,
"nbformat_minor": 4
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment