Interactive-3D-Solar-System.ipynb

{
 "metadata": {
  "name": ""
 },
 "nbformat": 3,
 "nbformat_minor": 0,
 "worksheets": [
  {
   "cells": [
    {
     "cell_type": "markdown",
     "metadata": {},
     "source": [
      "# Interactive 3D Solar System\n",
      "\n",
      "This notebook can be used in two ways.\n",
      "\n",
      "To enable a live Mayavi scene\n",
      "that you can rotate with the mouse\n",
      "and scale with the `+` and `-` keys,\n",
      "set the following value to `False`.\n",
      "\n",
      "But if you are a contributor to the project\n",
      "who wants to re-render the inline images\n",
      "for a fresh presentation of this notebook online,\n",
      "set the following value to `True`."
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "take_screenshots = True"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 1
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "%pylab wx\n",
      "if take_screenshots:\n",
      "    %pylab inline"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "Populating the interactive namespace from numpy and matplotlib\n",
        "Populating the interactive namespace from numpy and matplotlib\n"
       ]
      }
     ],
     "prompt_number": 2
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "from mayavi import mlab\n",
      "from mayavi.sources.builtin_surface import BuiltinSurface\n",
      "from sgp4.earth_gravity import wgs72\n",
      "from skyfield.constants import AU_KM\n",
      "from skyfield.api import JulianDate, nine_planets, earth\n",
      "\n",
      "pylab.rcParams['figure.figsize'] = (8.0, 8.0)\n",
      "\n",
      "mlab.figure(bgcolor=(0.0, 0.0, 0.0))\n",
      "planet_color = (0.7, 1.0, 1.0)\n",
      "yellow = (1.0, 1.0, 0.0)"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 3
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "# One Earth year\n",
      "\n",
      "day_array = arange(1, 366.3, 0.25)\n",
      "jd = JulianDate(utc=(2014, 6, day_array, 0, 0, 0))\n",
      "print 'Start:', jd.tt[0]\n",
      "print 'End:  ', jd.tt[-1]"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "Start: 2456809.50078\n",
        "End:   2457174.75078\n"
       ]
      }
     ],
     "prompt_number": 4
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "def reposition_camera(distance, focalpoint=(0.0, 0.0, 0.0)):\n",
      "    \"\"\"Move the camera above the solar system.\"\"\"\n",
      "    mlab.view(azimuth=90.0, elevation=360.0 - 23.4,\n",
      "              distance=distance, focalpoint=focalpoint)\n",
      "    if take_screenshots:\n",
      "        imshow(mlab.screenshot(antialiased=True))\n",
      "        axis('off')"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 5
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "# Get ready to draw!\n",
      "\n",
      "mlab.clf()"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 6
    },
    {
     "cell_type": "markdown",
     "metadata": {},
     "source": [
      "Our first step will be drawing the orbits of the planets\n",
      "over the period of a single Earth year.\n",
      "Note that we are plotting positions in astronomical units (AU),\n",
      "so that the radius of the Earth\u2019s orbit in the plot\n",
      "should be roughly 1.0.\n",
      "Positioning our camera at a distance of 10\u00a0AU\n",
      "shows us the four inner planets, out to the orbit of Mars:"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "for planet in nine_planets:\n",
      "    x, y, z = planet(jd).position.AU\n",
      "    mlab.plot3d(x, y, z, color=planet_color,\n",
      "                tube_radius=None, line_width=2.0)\n",
      "    print planet.jplname, '- plotted', len(x), 'points'\n",
      "\n",
      "reposition_camera(10.0)"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "mercury - plotted 1462 points\n",
        "venus - plotted 1462 points\n",
        "earth"
       ]
      },
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        " - plotted 1462 points\n",
        "mars - plotted 1462 points\n",
        "jupiter"
       ]
      },
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        " - plotted 1462 points\n",
        "saturn - plotted 1462 points\n",
        "uranus"
       ]
      },
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        " - plotted 1462 points\n",
        "neptune - plotted 1462 points\n",
        "pluto"
       ]
      },
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        " - plotted 1462 points\n"
       ]
      },
      {
       "metadata": {},
       "output_type": "display_data",
       "png": "iVBORw0KGgoAAAANSUhEUgAAAeIAAAGFCAYAAAAhGPBBAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3Xd0nNWd//H3FI00M5qRRr03W8VyEy64GzfANraxg8Em\nQAhZWoAAYUOW/Ai7YZPNZjdLzwZYuk0J3aa54oKLZLkXWb33PpJmNNJImuf3h4032dgUzYxGlr+v\nc57DiQ++9ytF6DPP89z7vSpFURSEEEII4RNqXxcghBBCXMokiIUQQggfkiAWQgghfEiCWAghhPAh\nCWIhhBDChySIhRBCCB+SIBZCCCF8SIJYCCGE8CEJYiGEEMKHtL6YVKVS+WJaIYQQwmcu1MhS7oiF\nEEIIH5IgFkIIIXxIglgIIYTwIQliIYQQwockiIUQQggfkiAWQgghfEiCWAghhPAhCWIhhBDChySI\nhRBCCB+SIBZCCCF8SIJYCCGE8CEJYiGEEMKHJIiFEEIIH5IgFkIIIXxIglgIIYTwIQliIYQQwock\niIUQQggfkiAWQgghfEiCWAghhPAhCWIhhBDChySIhRBCCB+SIBZCCCF8SIJYCCGE8CEJYiGEEMKH\nJIiFEEIIH5IgFkIIIXxIglgIIYTwIQliIYQQwockiIUQQggfkiAWQgghfEiCWAghhPAhCWIhhBDC\nhySIhRBCCB+SIBZCCCF8SIJYCCGE8CEJYiGEEMKHJIiFEEIIH5IgFkIIIXxIglgIIYTwIQliIYQQ\nwockiIUQQggfkiAWQgghfEiCWAghhPAhCWIhhBDChySIhRBCCB+SIBZCCCF8SIJYCCGE8CEJYiGE\nEMKHJIiFEEIIH5IgFkIIIXxIglgIIYTwIQliIYQQwockiIUQQggfkiAWQgghfEiCWAghhPAhCWIh\nhBDChySIhRBCCB+SIBZCCCF8SIJYCCGE8CEJYiGEEMKHJIiFEEIIH5IgFkIIIXxIglgIIYTwIQli\nIYQQwockiIUQQggfkiAWQgghfEiCWAghhPAhCWIhhBDChySIhRBCCB+SIBZCCCF8SIJYCCGE8CEJ\nYiGEEMKHJIiFEEIIH5IgFkIIIXxIglgIIYTwIQliIYQQwockiIUQQggfkiAWQgghfEiCWAghhPAh\nCWIhhBDChySIhRBCCB+SIBZCCCF8SIJYCCGE8CGtrwsQQnw/KrWaoPBwgiIiUGs0ABhMJoIiItCb\nTN/4d3tsNhw2G329vdja28/9eWdzMx1NTfT39Xm1diHE35MgFmIYCrRYMIWEEJ6YSKDFAoApJITk\nrCxMISH4GwzoDAZUKhUAfjod/gYDWn//bxy33+mk3+nENTCAs6fn3J87HQ56u7txDQxgbWykrriY\n5qoqulpbz1xtbXS1taG4XN77ooW4RKkURVGGfNKzvzyEuNT5BQRgDArCHBpKSGwsarWasXPnEmix\noD97lxsQGAiAPjCQyJQUAoxG7FYrdqv1XDD2OhzYrVacDsc3zqcLCECn16Px80N/dlwAQ1AQxuBg\nNFotdquV1tpaOpub6e7qwvH1ZbNht1opOXiQ5qqqMzV0dNDb3e29b5AQI8iF4laCWIghFhQejiU6\nGktUFJHJyUQkJmIODycsLg6VSsX4+fPp7+vD0dVFZ3MzPWeDrqeri6aKChw2G7azQewaGACgt7sb\nu9X6raHor9ej0+vR+vkR8FePsY3BwQSeDWJDcDAhUVGYw8IIMJnQBwae+2d3ZydFBw7QVFFx7sOA\ntbGR+rIy2uvraa+vx2GzyZ2zEOchQSyED6g1GvQmEyExMQDEjB5NTFoasWlpxKanE5mSgsFspufs\n3aaiKNQVF9PR3ExXayuN5eXYz77L7WptpfTo0b95t+sNlqgoYtLSiEhMxBQaiik0FHNYGKbQUAIM\nBiKSkggMCTkT6gYDPV1dFOXmUltURF1hIS01NdSVlGBtaMDR1cVAf79X6xXiYiFBLMQQ0Wi16E0m\n/I1GjEFBjJ48mcw5cwCYs2YNarUah82Go6uL9sZGmquqaK2upq6oCJfLxe533qHvr97fDifmsDCm\nLltGwrhxBEdEEBQRQVB4OIEhIegDA9GbTNitVg5s3EhhTg7NlZW01tXR3tBAr90ui8HEJU2CWAgv\n0wUEoPHzIyQmhsxZs0gYO5aIpCSmLl9+7t1tj81GxcmTlB46RMnhw5QeOUJTZaWPK3ePJTqay66+\nmtGTJzNq0iTCExIwmM1nvh9aLad272bbq69SdeoUbXV1OGw2nA7HBX8pCTFSSRAL4S0qFSpg6vLl\nxGdkMHvNGkZNmgSceXdbXVDAwU8+AeDTZ5/1+qNlXwuNjeXqO+9k0uLFxKanYzSbQaXC2dND9enT\nbHzqKfZ98AF9vb0gYSwuIRLEQniYn05HRHIyly9fztTly4lNS8PfYMDpcFB29ChVp0/TWFbG3vff\n/987Yrt9xC9kUms0BBiN+Pn7o9FqmbBgAYt+8hMSxo0j0GKh5+wd8ZHNm9n+2muUHT1Kj93u67KF\n8DoJYiE8JCgigqTx48mYOZNZ11+PKSQEU2goRQcO0FxVxZHNmyk+dIheu52+nh46W1t9XbJPGYOC\nCI6KItBiISY1lbk33sjERYvosdmwNjbSWlND7mefkfPxxzRXVfm6XCG8RoJYCDfFZ2YSNWoUY2bM\nYMKCBZhCQ9HqdOTv20fB/v2UHTuG3WqlpaYGu9Xq63KHHbVajd5sJnr0aMLi4kibNo3JS5cSkZhI\nZ3MzDWVllB45wtaXX6ahtNTX5QrhcRLEQgxS8sSJBIWHM+3aa0maOBGD2YziclF56hRHtmyhvriY\nuuJiujs6cI3wx86eFJmSQvKECWc+3MyaReasWbhcLk7t2kVTZSXbXn2VxrIy2f4kRowLxa20uBTi\nPLRnVz+HxMay8Mc/JiwuDnNYGCqViuKDBzm1axe1hYUUHzrk61IvWo1lZTSWleFvMFCUm0tBdjbJ\nEycyceFCjEFBBFosVJ46xandu6kvKZEOXmLEkjtiIf4PvclESlYWk5csISY1lYSxY3G5XBzetIny\nY8eoOHGCipMnfV3miBSdmsqCH/2IuPR0MufMQa3RcGDjRk7s2EHx2daa/U6nr8sUYlDk0bQQ38IY\nHIwlKorEceOYe+ONjJ8/n+7OTna//TZN5eUc276dxvJyX5d5SdCbTPzw8ccZPXkyMampOGw2Dmzc\nyKEvvqC2oIC2+npflyjE9yZBLMQF+Ov1aP39mTB/Pgt//GNGT51Kv9PJyV27KD96lO2vvUZ3Z6ev\ny7wkWSIjWXb//UxeupSolBTa6urY/dZbfPnGG3Q0Nw/bDmRCnI8EsRDnodZoyLrySpInTmTS4sWM\nmTWLuqIiPnn6aba+/LKvyxNnxY0ZwzX33svM1asJCgujKi+PF+69l6LcXHlULS4aEsRCnMfCW29l\n+QMPEJOWRsnBg3z1zjuc2LmTlurqvzmvV/iWWqNBp9eTPHEiV91+O/NuvpnOlhby9uzhuZ/8BIfN\n5usShfhWEsRCnOVvMDBx0SIsUVEsveceujs6aKqo4NAXX3B0yxYcNptsmRmmdAEBmMLCyJgxgxse\nfZSolBSaKit589FHOblrl7xCEMOaBLEQnDkLeNbq1Vx9113o9Hqq8/LYuX49NQUFWJua6LrEu2Bd\nLAxmM6MmTyZz1ixW/+pX1BUXU3jgAFtfeokS2VImhikJYnHJC4mJYfrKlSy67Tb6nU5s7e188fzz\nnN67l+6ODl+XJ74ntVZLSFQUWVdeyZW3305YfDxVp06x/dVXObxpk/SvFsOOBLG4ZIXExJBy2WWM\nmzuX1KlTcdhsbH/1Vbo7OynKzaVH3i9e9KavXMm4K65g6rJlWJua2LluHUe3bpXtZmJYkSAWl6SI\npCSmXnMNM1evJiolhbKjR/n8T3/i2Pbtvi5NeFhYXBwzV69m7o034ufvT+6nn3Jk82by9+3zdWlC\nABLE4hIUEhPDvJtvZsaqVRjMZopyc9nzl79wZMsWX5cmvGjmddex9J57iExJofLkSd7653+m4vjx\nC/4SFGKoSBCLS4opJIQrb7+dOTfcQGtdHTvXr+f49u3Y2tt9XZoYAmmXX86Se+4hc9YsmioqeOWh\nh6gvLcXpcEggC5+RIBYjn0qFRqNB4+fHsvvuY/Fdd1GYk8Onzz5LUW6ur6sbUiq1GrVa/Vd/oPrG\n/+6Us6dGKS4XiqKMiLAKDAlh2vLl3P3nP+NvMPAvV19N/r599HR3wwj4+sTFR4JYjHgBRiNp06cz\nY9Uq5qxZw+EvvuCTZ5+l/NgxXAMDvi5vyPj5+xMaG0tITMy5PwuKiEAXEHDBv9NaV4erv5/OlhZa\na2txdHUNRalepVKpMAYHM2nxYh5av56+3l7+68YbObZ9u5zkJHxCgliMaMGRkVy+YgVX3X47sWlp\nHNi4kU+fe46qU6fo6+31dXleERwZiSUqivjMTEwhIViioki+7DJiRo9G4+eHRvu/p5xqtNpv/O+u\nv78fFAXXwAADfX24XC66Wlspys0l76uvUFwuyo4do62+noG+vqH48jxCpVYTYDQSEhPDL95+m0CL\nhTcfe4wDGzbI9iYx5CSIxYgVGhvL1OXLmbNmDRqtlsKcHHa99Ra1BQUjpk2lKSSE+MxMIpOT0ep0\njLviCsLj4/E3GtGbTGh1OnQBARiCgtD6+dHR3ExnS8u5v29rb//GnszmsDDUGg3G4GCCwsPxNxjO\n7bXubG0FRcHe0UFfTw/1paUUZGdTW1hIc1UVrTU1Q/EtcNvoKVP4yX/9F8agIL58/XX2vPce7XKK\nkxhCEsRiRAqLj2fqsmVMWryYXrudAxs3UnrkCI0VFRfVndv5hMXFET16NMagIOIyM7nsyisJtFhQ\nazREJCXR53Rit1pprqyku7MTW1sb9SUltNTU0GO3/83jV6fD8Y2P5/0NBlRqNX4BAQQYDGh1OvQm\nE3FjxpA0fjwqlYro0aMxhYbi6OqipbqazpYW7FYrlXl5HP7iCzpbWmgoLR3W75fHz5/P6kcewRwe\nzp533mHXm2/KkYpiyEgQixHHHBZG1qJFTF66FID9H37IyZ07L+p+wyq1msxZs9D4+Z05FWrCBAJM\nJgKMRvz8/bE2NjLQ10flqVM0V1djt1qxNjTQ291Nd2cnrbW1f3Mn7A5/o5GIhASiRo1CpVIRGheH\nMSiIkJgY4jIyCI2LIyg8HNfAAAXZ2XR3dNBaW8vhTZuoLymhq7WV/mH4YWjS4sVcdfvtGIKC2LV+\nPTvWrfN1SeISIUEsRhSdXs+oSZMYO2cOerOZkoMHObFjB/aLsFWlTq8nODKS0NhYkrOymLhwIX46\nHVGjRqG4XPR0d1NfUsLJnTtpra1loK+P0iNHsLW14Tq72nkohcbGkjh+PBGJiViio0nIzCQ5Kwu9\nyYQlMpL9H31ETUEBzVVVlB8/Tk1+/rD7cDRu3jwW3HILhqAg9rz7Lqd276ajqcnXZYkRToJYjBjB\nkZGkZGWRNHEidquV4oMHaSwru+hC2BgcTKDFQlRKCuPnzSNpwgSmLl9OQ1kZroEB8vfto+TQIWzt\n7dSXlFB88KCvSz6vyORkZv7gBwRFRBCXkUHC2LGYw8MZ6O/nxJdfcnTbNmry86ktKhpWYTd27lyu\nue8+okeN4vP//m8ObNiArb19WD9aFxc3CWIxIgRaLMxcvZqZ111HV0sL219/neMXWbtKg9mMTq8n\nY+ZMMmbMIGn8eNJnzMDR1UV7fT273nyTvp4eDm/eTHNVla/L/d5WPPgg6dOnE5mcTFhcHKbQUGqL\nitjy4ouc+uorGsvL6bXbh0XgpU2bxqpf/IKolBQ+/9Of2Pf++/QMk9rEyCNBLC5uKhV+Oh1X/PCH\nXHXnnbRUVbFz/Xryvvpq2D32vBC1RoOfvz9TrrmGxLFjueyqq0jKyqK9ro7T+/aRv3cvW156yddl\neowpNJR5N93EFTfdRExqKgFGI12trfz3XXdRlJuLvaODfqfzXDMRXxk1eTIPvv46UcnJ/OuyZRQf\nPEhvd7eEsfC4C/5MKT4AyCXX97r8AgKU8fPnK08dOaL8x759yvSVKxWD2ezzur7rpVarlfD4eGXF\ngw8q/1Naqmzo71eeyM1V7nruOWXOmjWKRqtV1BqNz+v0+Net0ShaPz8lc84c5fHNm5VPXC7lo95e\n5cWSEmXujTcqgRaLz2tUqVRKTGqq8nJlpfJhT4+StWiRogsI8Hldco2860LkjlgMe1qdjpjRo1ly\nzz1MXLSIve++y4433qCxosLnd1PfhTksjKtuv53Fd92F3mym7MgRGsrKyP3kEwpycujr6RnxzSX8\n/P0xWixEJiVxxzPPkDxhAt1dXZQeOcKHf/gDJ3ft8ml9Gq2WyJQUnszNpdfh4I9r13Jq926f1iRG\nngvFrQSxGPYsUVEs/PGPWf7AA2faVj7zDLWFhcO+Y1ZYfDyjJk0ifswYlvz0p2j8/MjbvZstL71E\nY3k5nS0tI6KV5Peh1emIHzOGqFGjWHzXXSSNH09XWxs5Gzbw1TvvUJWX57Pa1BoNU5ct467nnqOx\nvJy3f/MbTu7c6bN6xMgjQSwuSmFxccxavZr5P/oRdcXFbHzyScqPHx/2HbPix4xhztq1TF66FBVQ\nW1jIvg8+oKmyktqiInpH+B3wt/E3GEiaOJH0adNYcOutmMPCKNi/n+yPP2bPX/7i07pm/OAHrHzo\nIayNjWx8+mmOyrGZwkMuFLfa8/6pEMNAeGIiM1etYsYPfkBdcTGfPPUUJUeODPuOWfGZmcy76SYm\nLFjAgNNJfWkpnz33HCWHD18Uj9KHQm93N4XZ2dSXlNBcVcXlK1aQOmUKEYmJhMXFsfWll3yyHa23\nu5ucDRswh4Wx4Ec/4qrbb8fpcJD31VdDXou4dEgQi2EpODKSKUuWMH3VKro7O/niz3+mICfH12V9\no6iUFJInTmTiokUkZGZSlZfHkc2baaurG7Z7gH2ts7mZ7I8+or64mKnLlzNl6VJW/uM/4qfTcWLn\nTkqPHBnyVxA9Nht733sPS1QUk5cs4YqbbsLa2EhtYeGQ1iEuHRLEYtjx8/cnfcYMpq1ciU6vZ/OL\nLw7rhTMqlYqY1FTm/vCHTF+1Cn1gIEc2b2b7a69RcuiQr8u7KFScPElLTQ3V+fksuftuVv7iFyRn\nZbFj3TqfbFFrq6tj/4cfEp6YyNg5c2iuqGDjU08N+1ci4uKk/vZ/RYihFZGUxJQlSwhPSODkzp0c\nGcbv6DRaLeEJCVxz333Mu+kmtFotORs2sO3llyWEvydbezs5H3/Mn+++m9N795I+fTo3/+53jJ8/\nH4PZPOT1FB88yKHPP8c1MMC0a69l1OTJQ16DuER4f9fw32MY7OeSa3heAUajsuoXv1CePXZMeeSD\nD5T06dN9XtOFLrVGo0QkJip3PP20sr6xUXngtdeUzNmzlQCj0ee1XeyXv8GgrHjgAeXF4mLl6SNH\nlJnXXeeT76vBbFZWPPCA8kZ9vfL41q2KOSzM598buS7e64KZOIT5+7+TDoNviFzD81p0223KcydO\nKI9v2aJMuvpqn9fzTVd4QoJyx9NPKxv6+5Uf/8d/KMGRkT6v6TtdKpXva/iOV/q0acpzJ08q6xsb\nlUU/+YlPagiJjlZW/+pXyrudncqvPvxQUavVPv++yHVxXhci74jFsDF56VKW3nsvdquVLS+95PMm\nD98kNj2dpT/9KbOuv54NTzzBhiefpLO52ddl4W8wYImOJjgiAlQqQqKjCYmJQevn93f/bn9fH+0N\nDQC019fTVldHe309vQ7HUJd9QaVHj/LS/ffzw9/+ljWPPUaA0chnzz03pDVYm5oozM6moayMjJkz\nSZ8xg5JDh4b9PnZx8ZAgFj6nUqsJjoggfdo0tH5+FOfl0VBaOmx/0UUmJzN37VomLFzIrrfe4pNn\nnqGjudknW5Ni09JIHD/+3ClOo6dOJTgiAq1OB5xpoKHV6c67d19RFPqdTgD6nc5zl81qpfLkSUoP\nHz6zgOr0aTp89CGj3+mk8MABXn/4YVY9/DA/+OUvCQoP561//uchq8E1MEDlyZNsev557nz2WW76\n7W/549q1w+okKXFxk4Yewuf8/P2Zs2YN1/785zSUl7Pt5ZfJ27sXxzA8zCElK4vpq1aRkpVFTWEh\nW158kfrS0iGbPyAwkLRp0zCHhnL58uWExMZiCgnBz98fvcmEOSwM18AADWVl5/6Ora2NztbWc6EL\n4KfTYQgKQqVWYw4LA86cbGUOC0OlUtHV1kZnSws9Nht2qxWHzcb2V1+lra6OhtJSHDbbkH3NcOZn\nJG3aNJbffz8RSUm8//vfk/3RR0M2v1qjITwhgavvvJNl993H1ldeYcMTT9BSXT1kNYiL34XiVoJY\n+JS/wUDatGnc/NvfotFq+fSZZziydStdra2+Lu3vpGRlseDWW0kYO5aC7Gz2/OUvVOfne31eU0gI\nkcnJJIwdy5jZs4lKSUEfGEhyVhYdzc00lpfTY7Nha2+nsbyc5qqqv3lM3mO302OzMdDff+7PNFot\nOoMBlUpFQGAgKsDfaEQfGIglOpq4MWOIGT0ao8VCYHAwBrOZvD17sFutdDQ1cXLXLopyc2mvrx+y\nPtl+/v5kzJzJ/JtvJnr0aPZ9+CGfPfvskMwNoPXzIzYjg0fefx9/o5EX77uP419+Sc8QfygRF68L\nxa08mhY+FRAYyGVXXknatGlsev558vfvH5YhnDR+PAtuvZWUrCyKDx1i3wcfeD2EIxITCYqIIHXq\nVDJnzyZuzBgSMjNpraujrbaW3W+9RWVeHo1lZWfuXDs6aKqsdPuRqSk0lNi0NCJTUjAGBREYEkJ4\nfDypU6eSMHYsxuBgkrOyKD1yhOrTpynKzaXy1CmvB1Jfby8F2dmgKFxz330s/9nPcHZ3s3P9+iF5\njdHf10djWRlbXnqJNb/+NbOuv56GsjIqT570+txiZJMgFj6jN5kYNWkSmXPmUH36NCd27sRmtfq6\nrL8Tl5HBop/8hNTLL6dg/352rlvntV++KpWK8IQEAkwmpq1YQcLYscSkpREUFkZ7QwP7PviAihMn\nqC8uPnNH2tjo8ZafXa2tFGRnnwk9zjyWNYeGMnnpUmLT0wlPSCA2LY3ZN9xAX08PeXv2cHz7dsqO\nHaMmP58eu91rZ/n29fRQkJ1Nv9PJmsce48bHH6fX4eDgp58OSdOPXoeDHW+8wegpU8icPZua/Hy6\nOzporqry+txi5JJH08In1BoNyRMnsubXvyZz7lzWPfIIe959d1idRqRSqwkKD+e6X/6SSUuWcGTz\nZra98opXTgjS6nSYQkPxNxiYd/PNhMXFkT5tGsbgYKry8sjft4/SI0coPHDA508MAi0Wxs+fz5Sl\nS0kcPx5LVBT+BgOn9+5l73vvUZCdTVttLf1e7Amu1elIHDuW+197jaDwcJ6/5x5O7NgxZD8/8ZmZ\nPPzOO+jNZt589FF2v/32kMwrLm4XjFtv7xk+H4bBfi65fHsFWizKkrvvVj5wOJR//+orxc/f3+c1\n/fWlUquVQItFWX7//cprNTXKL95+W0nJyvLKXFqdTonPzFRW/+pXyl1/+pOyvqlJebutTfnF228r\nqx95REmbNs3n348LXZHJycq1Dz2k/DE7W1nf1KR87HQqt/z+90ratGmKIShIUXlxz61fQICSOXu2\n8nxBgfLU4cPKlGuuUfwNhiH72n/10UfK+3a78sPf/EYJtFh8/v+FXMP/umAmDmH+/u+kw+AbIpfv\nLpVaraROnar84q23lNdra5Vrf/5zRaPV+ryuv770gYHKtGuvVV4sKVF+v2uXMnHhQo//klep1YpG\nq1VSsrKUn7/xhvKJy6V82NOjPPjGG8oPf/MbJXrUKJ9/H77rFRgSolx9553Kn/LylA8cDmWjy6Ws\nePBBxRwWpqg1Gq/+LKXPmKH8T2mp8uyJE8qUa64ZsoYllshI5bHPPlPeam1Vrv9//09RXUSNUuTy\nzXUh0mtaDLn4MWNY+OMfkzFzJrvffptNL7zwNyt6hwNDcDAzVq0iPCGBo1u3Up2fj9PDjS6SJ05k\n0W23ce1DDzF7zRrsHR1seOIJXnv4YT78z/+ksaLCo/N5k729nR3r1vFPs2ax7eWX6Wpt5ebf/Y6H\n33mHy5cv99q8istFycGDNJSVEZWcTGx6OgaTyWvz/TVrczPlx47hdDiIHzOG6NTUIZlXjEBDeCN8\nDsPgk4lcvrsmL16s/O7LL5X/zstT5qxd6/N6/u8VnpCg/PA3v1HeqK9X7n3xRSV61CiP37GPnjxZ\neeD115XXamqUJw8dUm7+3e+UCQsWKMGRkV59nDsUV0hMjJIwdqzy8F/+oqxvalL+dOqUsvZf/kUx\nh4Z6bc60adOUJw8eVP779Gll3s03D9nXGjdmjPJP77+vvFpdrax57DGff+/lGt7XhWh+85vf/IYh\n9vjjjw/1lGKYiMvI4Iof/pDE8eM5unUru958c8j2oX4XITExzPzBD5h/yy1UnDjB+//2bzSUl3v0\njj3ryitZ9rOfMemqq7C3t/PVO+/w5euvU1tQgL2jA4Z+/aRHObq66GhupiY/n+bKSsISErhs0SJC\nYmMpPHDA408WADpbWrC3t5M0YQJJEybQ3dlJTUGBx+f5v7o7O/HX64nLyCA0Lo7W2loa/6qZihB/\n7UJxK9uXxJCKy8ggeeJEnA4HhQcOYG1s9HVJ5/jr9aROncq0lSsZ6O9n2yuvUFtYiMuDrSsnLFjA\nkp/+lJCYGHI2buTEjh0U5uSMyO0vVXl5dDQ10VRRwawbbiDrqqtQFIVPn3mGhrIyXAMDHpurr7eX\no1u3kjZtGtNXrmTS1VdTnJtLa22tx+Y4n36nk4L9+8mYOZMp11zDpKuv5vj27V6dU4w8EsRiyEQm\nJzN27lz0ZjOn9+wZduf1hiUkMHbuXCISEzn0+efk7dnjsRDW+PmRMX0619x7L6aQEHI+/phDX3xB\nfUmJV+4Qh4uO5maObd9OW309tvZ2Zq1ejZ9Ox97336dg/36PPg3pamvj0BdfEJ+ZScaMGUxavJht\nr7zisfGR0OfKAAAgAElEQVQvpKmqioL9+0m7/HLSp08ndepUig8e9Pq8YuSQIBZDZvTkyYyZNYse\nu52Tu3bRMIQ9mr9NUEQEExcuJHXKFKry8sj+6CO6Ozo8MrZOr2fUpEmsePBBwuLj2bluHdkffURr\nXZ1Hxh/u+np7KT1yBJvVilanY9qKFYTExGBra6Py1CmPdsUqyM5m/4cfcvUddzBn7VpqCgrI37fP\nY+Ofj9PhoDAnh9N79zLzuuuYs2aNBLH4XmTVtBgSQeHhpF1+Of4GAxUnTgzJ+7vvI2nCBC5fvpzA\nkBByP/3UY0cw+hsMpGRlcdXtt5M4fjw7169nz3vvXTIh/Ncay8p4+7HH2PPuu0SmpDB12TJSLrsM\ng9nssTl6bDYObNzI4c2bSbv8cpb97GceG/ubNFVWUpidTa/DwbgrriAqJQW1Wn69iu9GflLEkBg3\nfz6Tliyhs6WFw59/TvXp074u6ZwAo5HEceMIjoqiobSUylOnPDKuLiCAlKwsrr7zTmIzMvjs2WfZ\n8cYbl/TxeZ2trbzxT//Ehv/6LyYsXMjy++8ndepU/AICPDaHtbGR6rw8emw2EsePJyolBbzcza+3\nu5vigwc59NlnRI8ezeK770br7+/VOcXIIUEsvE6j1TIqKwuD2UxdcTFNw2hhkkqtZuzcucxYtYqB\nvj72ffABhWd7LLtDo9UyatIklv3sZ4yfP58977zDF88/PyT9kIe7gf5+tr36Kke3bCFj5kwW3Hor\nSePHe7T1be5nn7H+0UeJSk7mn95/H7+z5zN7U3N1NXl79hAQGMjlK1ag9fPz+pxiZJAgFl6lUqkY\nNXky01aupLuzk8ObNlHthV7NgxUcEcHcG28kND6eo1u2cGzbNo+Mmz5jBtc+9BApWVlsfPJJNj3/\nvEdXCY8EH/3xj5zes4epy5cza/VqAi0Wj43dY7dTV1xMdUEBEYmJpE+fjp+X71B7HQ46mpvRaLWE\nxcWhkkfT4juSnxThVSqVitQpUzAGBVFbWEh7XZ1HtwO5a+LChSRnZVFfXExBdjadLS1ujzl+/nxW\n3H8/senp7HnvPXasW+fVAxAuVv1OJ6/98pcc37aNq26/nWt//nPPDa4oVJ8+zafPPIPebOZHf/iD\nR4P+QnOWHzvGUz/6EbqAAB7dsIGIpCTvzilGBAli4TVqjQZLdDSLbrsNgJwNG6gaLu+GVSqCwsOZ\nvHQpAYGB5O/bR8WJE27ftY6aNIlFt91GwrhxHNu2jW2vvIJ9GB7tOFy019ez4YknOLl7N7PXrmXl\nL37hsbG7OzooOXSI03v3kjxxIhMXLcIYHOyx8c/H3tFB2dGjqDUa0qZNw1+v9+p8YmSQIBZeowsI\nYOKiRUSnplJ29ChVp04Nm2MOVSoViePGkTRhAp1NTVScOIHVzUVUxuBgpq9cydi5c7E2NXFq9+4R\n2ajD08qOHWPbyy9Tk5/PvJtuOvfBzV0ul4vWmhp2vP46fjod82+5BUtUlEfGvhDF5aKrtZXcTz9F\no9EwcdEigiIivDqnuPhJEAuv8fP3Z/TkyWi0WqpOnz7TvnGYUKlUjJ46FXNYGHUlJTRXVdHvdLo1\n5qhJkxg3bx7dHR0c+uwzSg4f9lC1I1tfby/5+/ZxdMsWjEFBXHHTTYQnJHhkpbPDZqNg/34URSF9\n+nSPbpW6kO7OTnauW0ef08n0VasIjYnx+pzi4iZBLLxCq9MRlpBA6tSpdLW2Unb0KLZh8ohWrdEQ\nGhvLxIUL6evpoWD/flpqagY9nkqlInr0aOauXUuA0ciBjRvJ2bCBVjfGvNTYOzooyM4mf98+EsaO\nZc7atR5ZRe0aGMDa3Mzut95CrdEw5ZprCE9M9EDFF+bs6aEgOxtXfz+jLrsMo7ffTYuLnnTWEl6h\n0+uJSU0latQoCrKzKczJwd7e7uuygDO1TViwgIwZMzi1ezf5+/YNem+vSqUiwGhkztq1XH7tteTv\n3cvBzz+nrrjYw1W7LyAwkECLBX+DAQCNRvM370ztHR302O24BgYY6O+nx2Yb0lcJ1fn57H77beIy\nMrjyH/6BguxsinNz3e685ejq4rWHH2bCwoUs+9nPKDl0iObKSg9V/fdcAwN0trZSmJND5pw5JE+c\neK7vthDnI0EsvMI1MEBHczPZH39M5YkTw+rdcGBwMDOvuw5UKo5u3Up7ff2gx9P6+RGfmcm8m27C\nNTBA8aFDbo3nSSqViqCICFRqNaaQEGLS0kgaP56Qs49K/Q0G4jMzz/37tQUFNFdX43Q46LHbaaqo\noK6oCFt7OzarlV67HcWLJ0M5HQ4K9u9n41NP8eDrr3P7U0/x2+XL3f5+Ki4XHc3NtNfXkzxxIpao\nKAKMRq+e+jXgdPLmY4/xy3ff5ao77qCuuJhDn3/utfnExU2CWHhFj83GqV27OOWhVpGeogsIIDY9\nnSnXXENdSQnHtm0b/CItlQq92cycNWsIi4tj68svs/fdd326QEulVqP180Or02Ewm1l8113o9Hpm\nXncd4fHxDAwM4BoYOG+gxmVkoFKpUGs0Zy61ms7WVnI/+YTcTz6h5PBh7FYrTofDa1vQ7B0dVJw8\nSXdXF6MnT8YSGUlnSwsDHtj+tf3VV7nukUeYsnQpVadPe7UHtaIotNXW0t/XR1B4OAFGo9fmEhc/\nCWJxSTGHhTHlmmtQXC5yN250a2uRn05HXHo6yx94gNa6Og59/jntDQ0erPb7M4eGkj59OuPnz2fZ\nffedaSqhUqFWqagtLqYmP5/Wmhq62trO//fDwgiNiyMsLo7Q+HiCwsNZdNttLLj1VlAU3vqXf2HH\nunW0efF4wbqiIv516VL+Y98+fvXxx/z7dddRduSI2+Me/PxzrrrjDtJnzCAmNdXrQdxaX8/A2SD+\n+nWAEOcjQSwuKVp/f0JiYs78ojx7xzJY5rAwJi9dCsDBTz6hprDQp0cajpo0iavvvJPpK1ei0+tR\nqdUc+uILWqqr2ff++9SVlDDQ18dAf/8F90urNRo0Wu25f4YnJDDvpptImTSJUZMm8YOHH2bJT3/K\nV++8w5evv05Nfr7Hv45+p5OWmhoURSEsLg6dhzpitdXX09fbizksDL3J5JExv5Gi8NrDD3PrH/7A\n4rvvxt7RQc7HH3t/XnHRkSAWlwxdQAAxqamMnzePtvp6Dm/aNOijDnUBASSOG8f8W26hra6O7I8/\npqu11avvUC8kOSuL1ClTmLhoERMXLsQ1MMDOdes4uXMnVfn5OLu7sTY24uzp+d5jWxsbaampwWAy\noTebueqOO5iydCkLb72VxHHj2P/BB+Rs2IDNgwvxFEWhs7WVP65Zwz0vvMDN//ZvrH/0Ubd7gLv6\n+/ny9dcJDA7msquvpiovjxM7dnio6vOryc/H0dVFVEoKlshIr84lLl4SxOKSYbRYSBg7Fr3JROH2\n7bQ3NDDQ3z+osYKjosicMwejxcLed9+l8uTJQQWdu8bPm8fCH/+Y0VOm4OjqYu9771GYk0NRbi4t\nZxdeufPhoK+3l6aKinP/u7uzk5yPPmLy0qWMnTOHax96iLiMDDa98AKN5eUe+IrO6O/t5fDmzfQ6\nHGRMn+6xRhyn9+xh9g03kDhuHPGZmV4P4vazH4ACQ0LQD8EeZnFxkiAWlwxzaCgJY8fS19tLYU7O\noBcA+QUEEJeRwfj58+nv7aXs6FG6u7pQhriHduacOSy+6y4iU1Loam3lxM6d7Hv/fZoqKry2Irjy\n5EkqT56koayMmoICJi9ZwrSVKwkIDGTve++R99VXHnkqoCgKvd3d7Fy/niv/4R+YsmQJDaWlVJw4\n4da4zVVVOLq6CI6MHJKOV06Hg+LcXOIzMogfM4b4MWOo9sLjfHFxkyAWlwStTkdkcjKjp0yho7mZ\nY9u2DbqTlikkhOSsLGLT0miqqKDwwAGPrOr9rtQaDTFpaSy9+25i09M5sWMHRQcOUH7iBLWFhUNS\nQ+WpU7TW1VGTn8+ctWuZvmoVQeHhNFVU0FJT47GTpr565x2mLlvGhEWLOLFrl9tB3GO301heTm93\nN6GxsYTGxtLqxYVnAIc3bSJ9+nRGT5nCmFmzJIjF35HOWuKSoAsIIDwhgZjUVKyNjZQcOjTox9Ih\nMTEkjh0LQEFOjltjfV+6gABSsrJY8tOfkpyVxeHNm9n99tsc/OyzIQvhr9na2jiyZQsbnniCE19+\nScK4ccxes4aMmTMJCAz0yBzt9fXUFhai9fMjZvRogj3wnvXkrl1Unz5N4rhxjJ071wNVfrPaoiI6\nW1oIjYsjXE5jEuchQSwuCf56PYagIPqdTrpaWwc9jkarJTIpibgxY+hsbqYwJ8eDVX4zrU5HbHo6\nKx54gDnXX8/Bzz5j0/PPU3rkCL0+XK1deeoUrz78MNkffcTiu+5i7WOPET1qFGoPnMfb1dbGnr/8\nhY6mJrKuvJKMmTPdHrMwO5u64mJi09NJmzbN7fG+Tb/TiWtgAI1Gg1YrDyHF35MgFpeEqFGjGD15\nMra2No5v3z7ocYxBQaRcdhmx6elUFxS4Ndb3odZoiExOZt4ttzBt5Ur2fvABH/7nf9JSXT0k838b\na0MDG554gpaaGsbOnUv69OmYwsLO7GN2U0F2Nu0NDcRnZhKXkeH2eK6BARSXC7Va7ZEPC9/G1t5O\nV1sbWp0Os+wpFuchQSwuCXqTCVNoKH29vbS50TLRGByMMTgYZ3c3rTU1Q9bAIzgighmrVrH8Zz+j\nvrSUT59+etBbr7zF1t7O/7viCtrq6rjnhReYff316D30iJqzC8A8cRBEV1sbPTbbuZ8JT5zy9G2O\nbtlC0YEDJE+YwLRrr/X6fOLiIkEsLgmWmBhi0tLo7upyq0tT8sSJJGdl0dnaSuWpUx6s8JtFjRpF\n5uzZ2K1Wcj76iObq6iF7L/19/X7lSiqOH+fahx5i6rJlbo/X0dTE8S+/pKOpibTLLyd9+nS3xyw7\nepTy48cJT0hg/BVXuD3etyk+eJDq06eJGjWKMbNne30+cXGRIBYjnsFkIiw2Fn1gIK21tbTW1Q16\nrIiUFCKTk+m124fsbnjUpEks+NGPiE1PZ/+HH7LlpZfc6gjmbdUFBbz9+ON0NDWx+O67zxyw4QaX\ny0XB/v3Ul5YSN2aMR97rdra20tXWhj4w8NwhGN7U73TS39+PRqvFT6fz+nzi4iJBLEY8tVaL1t8f\nlVp95hfiILctGcxmwuPi8DcYqCkspOjAAQ9X+vf0JhMZM2Ywbt48+vv6KMjJOfMBwAcdvL6rfqeT\nkzt2YG1sZNSkSSSNH+/2mC01NdhaWzGFhhIaG+v+eNXVtNbUEHj2VCpv6+7qotduR6fXYwwK8vp8\n4uIiQSxGPH+DgQCDgX6nk+7OzkGPYwoJwRQaikqlwtbejrWx0YNVnl9cejoZM2bg5+/P6T17OL1n\nj9fn9AR7Rwc5GzZQV1RE1lVXMXP1arfG63c6Gfj6jtIDvac7mprobGnBYDYTFh/v9njfpt/ppLG8\nnJbqaoKjokgcN87rc4qLhwSxGPEsUVEER0XR3dFBY1nZoMcJi4/HFBqKo7OTrpYWD1Z4YfFjxxKf\nmUlnSwt5X33l0TaS3nb4iy8oO3qUlMsu47KrrnJrrN7ubvqcTvz8/dHp9W7X5uztpa+nB42f35Ct\nYu5obqa9oQF9YCARsp9Y/BUJYjHiWaKjCYmOxt7RQUNp6aDHCU9IICg8nLb6euqKiz1Y4fkZg4NJ\nHDeOoMhIGsrKKDt2zOtzepK1qYny48exNjYSkZhIwtkmKIPRY7fTVltLr92OJTJySN7repqtrY3O\nlhb8jUZC4+J8XY4YRiSIxYgXYDQSEBhIX28vdje2/OhNJnR6PV1tbV5viwgQl5FB8sSJKC4XJYcO\nUZWX5/U5Pe3Ejh2c3LmT+MxM5t54o1tj1RQWUldUhCU6mtFTprg11tfHQao1miFbPFVfUkJNfj5B\n4eGkD0EjEXHxkCAW4jvS6fVodToUl4sBD/VS/iaxaWmYw8KoKyqi8uRJr8/nDS3V1bRUV2MKCSHW\nzWYcXWdXOgcYjYRER7s1lqOrC4fNhi4gAFNIiFtjfec5bTa6Ozvx8/fHFBY2JHOKi4MEsRjxVCoV\nKpUKRVFwDfKEJJVKRdq0acSkpuLo6vL6O2KVSkVQRAQBgYHYrFY63WjL6Uv2jg7sVit+AQEEh4ej\ncaPFY1NlJc1VVZjCwojPzHS7NkVRzpwUpVJ5pAOYEIMlP31ixAtPTCQiKYmulhaqBtmEIyg8HP+z\ni4Ta6uq8/o44ODKSrCuvxBgcTGF2NqWHD3t1Pm+qKSykKCeHsLg4Ft1226DHsVut2Ds68DcYCAoP\nd7uuttpa6oqLMZjNpGRluT2eEIMlQSxGvECLhUCLhd7u7kFvOVKp1aBS0d3Zic1q9dgxfxecT6VC\n6+eHSq1moL9/2HbR+i5qCwspPHAAc0QEU9zotKW4XCgu15knHB64g+3u7KSzpQW/gIAhWfzl6OrC\n3tGBLiDAIx8kxMghQSxGPLVGg1qjOfNu181AG+jv93oIA5hCQ9F6YL/scGBra6Otrg6NVothGDWz\ncA0M4OrvR6VSofHzG5L5miorqTh5ElNICGPnzPFI72xx8ZMgFmIY8vP3R63R+LoMjxjo7x90N7OR\nptdmo7OlBa2fHxY3F5yJkUOCWAghhPAhCWIhhBgi/X199PX0oNJoCDAafV2OGCYkiIUYhvp6e4fk\nXbQYWvaODlpra9H5+xOVkjIkZyGL4U+CWIx4X5+4pNZq0QUEuDWW5uzCL2/ram2lv7fX6/NcTDRa\nrccW3cHZRXxaLYqiDNmq9K7WVpoqK9Hp9cRnZiIxLECCWFwCmquqaKqsxBQSMuhGEIrLBYqCISgI\nY1CQrHb9Hr5uqOIug9mMISiIXoeDTg80VDGYzZhDQ+nr6aHNjTOqhXCXBLEY8ewdHf/bCCIiYlBj\ndDQ30+twAGd6TptCQz1Z4og2XL9fITExxKSm0t3ZSfnRo74uR1zCJIjFyKcoZy437swURaE4N5e6\n4mIsUVFEjx7t4SL/VkdzM329vQQGBxOXkXFRH5tniY4mevRoHF1dVBw/PuhxwhMTCY+Pp6ulhZr8\nfLdqMgYFYQwOxtnTQ0dz86Bbn35fPXY7Xa2taHU6gqOj5R2xACSIxSXA0dWFo6vL7Qb/vd3d9Dud\nGIOCCB7knfV3NdDfz7Ft22goKyMiMfHMwp6LlDEoiKDwcLo7Oylz487TFBpKYEjImSMR6+vdqkml\n0aBWq4f0/TCcecXh6u8HlcqtvttiZJEgFiNea20trbW1GIODiU5NHfQ4NQUFtNbWEp6YSOKECR6s\n8PyOf/kldcXFxKSlkTxxotfn84YAo5H4zEwSx4/H4WYQG4ODMQYH09vdTUdTkwerFMK3JIjFiNfR\n1IS1sRGD2UxEYuKgx2koK8Pa0IAxOJjgyEgPVnh+dSUldDY3ExoTQ7gbdfuSJTqauDFjCDAaqS8t\npaGsbFDjBBiNRCYlEWix0FpbS2N5uVt1+ev16PR6+p1Oemw2t8YSwl0SxGLEc/b04OzpQePnR0Bg\n4KDH6e3upm8ItxT12u1UnDhBW3098WPGkD59+pDN7SkRSUlEJCbS2dJCxfHjOLq6BjWOKTSU8IQE\n1BoNLdXVgz6842vBUVEER0bS3dFBU0WFW2MJ4S4JYjHi9ff1MeB0otFo8HPjIIV+p5NehwNFUQi0\nWLz+nhjg2PbtFGRnk5KVxewbbvD6fJ7kp9ORkJlJRFISDeXl5O/fP+ixQqKjMVos2Nrb3X4/DBAa\nG0tIdDS29nbqS0vdHu+70mi1aP39UVwunA7HmUWE4pInQSxGvB6bjY6WFvr7+giKiMAYHDyocbo7\nO6kvKcHW3k78mDGMmT3bw5X+vebKSmry83H29hKRmOjWo/WhFhIbS9q0aZhCQijYv58TO3YMeqzI\nlBSCIyKwNjbS5OZjaQBzaCim0FAcNtuQ7iHW6fUYg4MZ6OujywN7ocXIIEEsLgkNpaUUHzqEKSSE\nrCuvHPQ4p/fuJX/vXsyhoSSOG+fBCs/P3tFB9scfc2TzZsbPn8+Njz/u9Tk9QaPVMnbOHGJSU6kp\nKKDowIFBj6VWq7FERaE3mbC1tWF1c6GWSq0mIimJ0Lg4Opub3d4KJYS7JIjFJcHR1UVXayt+/v6E\nuHH8XFNlJU0VFZjCwkgYgiAGqD59mtJDhwi0WEgegtXanpAwbhwLbr2VqFGjKDxwgMKcnEGPFRQR\nwYQFCwiKiKAoN9etUAcwhYQQFheHSq2mubqajuZmt8YTwl0SxOKSUFdURFFODuawMCYvXTrocbqt\n1nNduobiHTGcaSZy7Msv+fA//oOYtDR+/emnQ3KQvTsmzJ9PeEICFcePk79vH11tbYMeK2HsWMxh\nYWc+BFVWorj5XtUcGoreZKKnq8utugYjODKS2LQ0eux2Sg8fRt4QC5AgFpeI/r4+nD09qDUa/A2G\nQY/jcrmoOHGCU7t2EZGUxNJ77vFglRfWUl3NjnXr+OzZZxk7eza3P/30sGwbCTD7+uuZu3YtnS0t\n7Fy/nryvvnLrJKnEceMwhYbSVF5Oc2Wl2/VFjR5NcGQk7Q0NNA5yO9Vg6U0mgiMj6Xc6aa6ulsVa\nApAgFpeIXoeD7q4uNFqt2wFWffo0eXv2uP2++fvodzqpLynh8KZN9Pf1MXv1asLj44fVnbHWz4+x\nc+ZwzX330d/Xx7ZXXuHwpk3YrdZBj2kKCWH8ggUEBAZSfPAgVadPu12nJSoKY3AwtvZ2rA0Nbo/3\nfWg0GjR+figu15BuhRPDmwSxuCQM9PVRX1zMyV27CAoLY/YNNwx6K5O1qYnG8nL8DQZi09OJSkkZ\nkqMR+51Oqk6f5t3f/Q6b1crKf/xHksaNQ6vTeX3ub+NvMJAxaxbXPfIIQRER5H76KUe3bqXdzaBL\nnz6duPR0mioqqDhxApubj5L1JhMpWVkEWizUFhVRXVDg1nhCeIIEsbhkNJSVkfvpp+hNJhbddtug\nzybudzqpysvj8KZNhERHs/C229zan/x9dLW2svWll9j45JMER0ay8Cc/YeqyZW4tQHNXoMXCxEWL\nWHL33ajVara98grZH3/sVgirVCoMZjNTly9Hp9eT99VX1BUXu11rZHIyCWPH0t/XR1VeHq01NW6P\nKYS7JIjFJcNutdJQWopWpyNx/HjUbjTdrzhxgs0vvuh2qA+G0+Fg84svsvGpp8icPZtrf/5zZt9w\nA9GjRg1ZDV8LtFgYd8UVLL7zThLHjWPv+++z6YUXqCsqot/pHPS4aq2WmLQ0MufMwWa1UnjgAC0e\nCE1Pv2/+vr5+R9zndNJSXT3k84vhSYJYXDL6enuxNjZibWoiKCKCsLg4tIN8x9rX20tbXR0NpaUE\nR0SQnJWFTq/3cMXf7PCmTRz6/HMCQ0JYeu+9LL77bpLGjx90w5LvIzAkhLiMDKYuW8aV//APRKak\nsOfdd7E2Nrq1MOtrWj8/0qdPJzQ6mtrCQhrLy890onJnTJ2O0ZMnY7RYqCksHHTf68HSaLWEJyYS\nn5lJV2srBdnZbq8AFyODBLG4pLTV1fHla6+hVqu5+o47MIeHD3qsuuJi3njkETRaLT9fv57I5GQP\nVvrdvPnrX/PBv/87rbW1LL7rLh5ct45pK1YQFhdHgNE46POXL0RvMmEMCmLqNddwxzPPcNsf/0hE\nYiI7Xn+dbS+/zKEvvnA7MNVqNcbgYDJnz8blclF65AidHuhCFZGYyIQFCwAoyM6mZojfDwcEBmIw\nm+nr6aGztXVI5xbDmwSxuKR0trRweNMmVGo1U1escOvu0dHVRdmxY7hcLkJjYvDz9/d48H0Xu958\nk+duv50d69aRNGECD7z+Ok8dPsyM665Dp9ejVqvP1DXI2lQqFSq1GrVGw/xbbuHahx7i2gcfJOvK\nK2muquL93/+eD//zP8/0gPbAHZ7RYmHW9dcza/VqTuzYwf4PPvDIu9ypy5ZhCg2lMDvbI++b3SJ3\nwuKvSBCLS0q/00l7QwMqlYrQmJhBP5r+WltdHb+aOxe71crP160jfcYMD1X6/TRVVLD+0Uf53fLl\nnNy5E3NYGHc+8wz/U1rKY59/ztwbbyQuLQ2D2fy9xjWYzYybN49V//iP/NuuXfzw8ce59qGHcPb0\nsPXll9n45JPkbNzosa9DFxBA4rhxXHXHHbTW1bH5xRfd7gWt1mgIi4vj8hUrGBgYOHfO81AzBgUR\naLHgdDjcPj1KjCyDX60ixEVIURSsjY28cO+93PqHP3DDr3/NX/71X6k8eXJQ4/U7nZQfP07JkSOk\nTp7MnBtuwNbWNuSPPV0DA3R3dHBq927qiouJSEpixYMPkrVoEePmzCE+I4Meu51+p5PG8nKKcnOx\nNjb+zSNfc2joua1QwZGRJE6YQMzo0QQEBhIQGEigxUJhTg7NlZVkf/wxpUeO4OjsdPtR9F8Ljopi\n0tVXExYby+633qIqL8/t8TVaLaOnTiUmNZXawkJq8vPp7ujwUMXfXWhcHBFJSditVqry8oZ8fjF8\nSRCLS05fTw+n9+5loK+PUZMmEWixDHosRVHodTjY+OST3P2nPzFh4UKObd8+5EH8tR67nbriYpqr\nquhqa2PzCy+QPn06Y2bNImnCBAItFmJSU0m57DKcDgfOnp5zf1cXEIBKfeYhmU6vxxQSglqtpqaw\nkOKDB8nft4/y48exW6201tRg93CYhcTEcPmKFUxdtoxTu3fz+Z//TFdbm9sLmr5e+KU3mag6fdpj\nC8q+L1NICOawMHrsdtk2Jf6GBLG45AwMDNBSXc2ed99l6vLlTFm6lPb6+kE/rlRcLvJ272bziy8y\n/5ZbmH3DDXQ0N7t9OIE7+np7KT18mNLDh6k8dYqTu3YRFheHOSyM8MREYtPS0JvNqFSqcx9Eujs7\nsbe3n9t21FZXR2FuLtaGhjPfn5ISHJ2duFwur9QcGhtL5qxZmMLC2PrKK9QUFKC4OZff/2/vzqOj\nKuUJnyQAABlfSURBVA8+jv9mskwyM8lkJsuQPZksJIRAAlkMO7IEoQLWjVoV6Gs9trWt2HPU19b2\ntUWPx1dsbY8eRSmKSBVUdgsiJmGTkLAlIQnZSMhkmWQyS/b9vn9AeRFBxUzyTMLvc85zPIc/7n0y\nB/zmzr33eRSKS093L1kCm8mEgoMHHbKf8Q/h4uYGV3d3DA4MfO0XICKGmG450uAgOux2ZH/wAZIW\nLMCUO+5A8eHDQ7pv2NXejqz334chORkTZs6EubYWLUYjWurqHDjzH8Z04QJMl/fwdffwgF9oKMZF\nRcFTrQZkMqg0GkAmQ1dbGzrt9ishtplMqCkqGpE5+oeFYdLcudBHROBcTg4Ks7Ic8kCT0tsbSQsX\nIiQ+HlmbNuHC2bPoamtzwIxvjrunJwLCwqDx90d9eTkaKipGfA7kvBhiuiVJg4MwlpSguqAAMWlp\niJg0CeX5+UN6iKalrg7HPvkE+shIJGdmoqGiAof+9S+nuvrp7e5GfXm5+KeGr6Lx90fKkiWYungx\nbE1N2L9+PaoLCoZ8XDeFAuOiopC6ZAnsTU04sWsX2gS9NuShVEIfFQUfvR6F2dnCbl2Qc+JT03TL\narNYkLNlCzqsVkxZtAgxqalDPuaRrVvx1aefwkevx6yf/ASG5GQHzHTsUmu1SF60COnLlqG/txcH\n330XBVlZDjm2l68vEufMgWHKFBTl5ODc4cPobG11yLFvlkKphLtCgb6eHnR3dAiZAzkvhphuaecO\nHcLF4mLoDQbEpKYO6cGt//j3m2/i87ffRkh8PJY9+SR8g4MdMNOxx8XVFfHTpyPz5z+HSqPBFxs3\n4tgnnzjs2PrISEyeNw8dNhuObtuGns5Ohxz7hxgXHQ3fkBDYTCbU8WqYrsEQ0y3N3tyMj/7yF1Sd\nPo15q1Zh4c9/PuRjdthsOPjuu8h+/30kL1yIX2/YMGKbQowmuqAgRE+dCrmLC6oLCx265KSPXo/E\nOXMQm56O4sOHcXzHDqEhDoyKgm9ICKyNjfxamr6BIaZbnrG0FDaTCb4hIQ7bOKHFaMT+y7skjU9P\nx683bHCqvYOdwdyHH8bchx9GfXk5DmzYgMqTJx127HFRUZgwYwbarVbk7t4tfE3ncdHR8A0Ohs1k\ngvH8eaFzIefDENMtb3BgAAc3bsSJXbswcc4cLF2zZujHHByEubYWBzZswPrf/AbacePwq7fegmYI\na1uPFS5ubrj32WfRWFmJtx5/HNteeAHVZ89ioL/fIccPiYtDxl13ITA6Gsd37EDerl1Dfg1qKLx8\nfeEfGgqZXA6z0YjW5mZhcyHnxBATAZdWxzp5EmqtFomzZ8NLpxvyMQf6+2FpaMCpfftwYMMGpCxe\njJ88/7yQzSGcgVqrRcrixVi2Zg3m/+xnaKyqwrlDh9BYVeWwJ8uV3t6YdPvtSJg1C/UVFcj54AOH\nLzxys4JiYuATEABbYyOaqqsd9gsHjR0MMREurUh1ev9+nDlwAJFJSch89FGHHHdwYABtLS0oOHgQ\nX2zciKT583H/c88hftq0Ed3DWDRvPz8kLViARb/4BeavXo3SY8fQUleH7o4Oh65yFRofj4RZs+Dq\n7o6i7OwRew/6u+akCQhAc23tiG+9SKMD3yMmuqy6sBBZ778PtVaLOQ89BFN1NY5u2zbkUAwODqLV\nbMa/33wTbgoF0pYuhZdOhy/few+F2dlot1od9BM4J7/QUKQsXozUH/0I7kolaoqK8Nkbb6DNAVsb\nXm2cwYDb7roLAeHhOHfoEPL37hX6gBZwadvIuGnToNbpcObAAa4xTdfFEBNd1tfdjdKvvoI2MBD/\ntW4dlj7xBE5//jk6bLYh32McHBxEc00N9vzjHxjo60P68uVY9Nhj8A0NRfHhw6g6fdpBP4XzUPn4\nwEevR9L8+Zh+zz2QABz6179QXVCAivx8h35Fq1AqMXXx4ktbM9bU4OjHHzvF1bCPXo+wCRMAAHVl\nZVxjmq6LISa6SqfdjgunT6M8Lw9xGRmYumgRcnfudNgiDKYLF/DxSy/B2tiIWQ88gOVPPomYlBTs\ne+stXCwqEn4/01F8g4MxPiMD49PTEZWcjJ6uLuRs2YJT//43Wodhdavg2Fik/uhHcHFzw9mDB1GR\nn+/wc/wQAeHhUCiVsJlMaLNYRE+HnBTvERNdo7akBFvXrkWr2YxVL7+MyKQkh97PbbdasfOvf8WO\nV15BfVkZpmRm4rf//Cem3nEH9JGRcPPwgEwmc9j5RopMJoPKxweagADMfeghLP3tbzHj/vvR2daG\nna++iuzNmx0eYZlMBo2/P9KXL4cuMBDnjx/H+ePHhawnfT0xqalQ63SoO38ezTU1oqdDzkoSAAAH\nh1MPD5VKmnrHHdKWlhbplRMnpPG33Sa5uLo6/Dy6oCDpnmeekT6026WtbW3SC1lZUty0aZJaq5Vc\n3NwkmUwm/LP4riGTyyU3hUJSa7XSoscekx5cu1ZaX1kprcvLkx79+9+lxLlzh+e8MpnkoVJJP/nj\nH6V36+qkx99+W4pITBT+efxnuLi6Sn/54gtpa3u7tGzNGknp7S18Thxixw2bOIL9/f+TOsEHwsHx\nXcPdw0NKWrBA2mg0Sv/71VfSxNmzHX4OmUwmubq7S7Hp6dIrJ05IO/r7pY87O6WnPvpImjR3ruSh\nVgv/HL5r+AYHS/NWr5YeWrtW2tbZKX3a0yM9sXGjlLxwoeTt5yfJXVyG5byeXl7SvFWrpC0Wi/SH\nXbskQ3KyJJPLhX8eACS5XC6Nz8iQ1ldWSn8vKJBSFi+WMAp+qeIY3nEjvEdMdAN9PT0oPnIEebt3\nX9m3uKWuzqFb2EmShP7eXlw4cwav/vSniJs2Dfc++yySFixA/LRpKDl2DPmffYZzOTkwVVc77LyO\nEBofj6ipUxGTmoqZ998PSZJwet8+WBsasP/tt1F3/jx6e3qGZTENbz8/pC1dijt+8QsU5eRg35tv\norGyUujCHVeTubggNjUVHioVKvLzYWtqgiO2daSxiSEmugFJktDb1YXt69ZBHxWFjLvvhtloxMF3\n33X4Pci+nh40Vlai1WyGsaQEk+fPR8aPf4yEWbNgSE6GeeVKVOTnI2/PHhhLSmA3m4X8jz04NhZh\nCQkYn5GB6JQU+AUHo6+3F+V5eTDX1mLv66+jt6sLLfX16Bum7R/dPTwQPnEi5q1eDaW3Nz585x2U\n5eWhu719WM53s+QuLtAFBiL1zjvR19uLwqws3h+mb8UQE32HhooKbF27Fvc88wzmPPggBvr7kb15\ns8NjPDg4iHarFWUnTsBSX4+yEycQl5GBpAULYEhKQkhcHGLT0mBvaoLZaMSp/fthqa9HU00NuoZh\nez83hQJqnQ4aPz8ExcQgJj0dIXFx0I0bB9/QUHS3t6OpuhrVBQU4+vHH6GprQ21x8bCv6xwQHo7b\nli9HoMGA7M2bcf7ECYe8YuYobgoFDElJiJg0CcbSUlQXFIz5d8VpaBhiou/h3KFDcFMosPCRRzD9\nnnsASUL2Bx8M29O5ZqMRZqMRxpISVJ46hbCEBASPH4/YtDRETJoEmVyO8MREtLW0wGYywVxbe+U9\n1YH+fnTa7eiw29Fpt39rGGUyGZQaDVQaDWRyOXz0enioVFBqNAiMjkZgVBRUGg18Q0IQEh+P7vZ2\ntBiNOPnZZ6g8dQr15eWw1NXhYnHxsHwO1/IPD8eURYsQmpCAvD178OWmTWi3WJwmwgDg7umJ2Ntu\ng9LbGzUFBbA2NnJZS/pWDDHR93TmwAH0dXdj9gMPYEpmJvp7e3H4o4+GdaP3lro6tNTV4ewXX0Af\nGYn4GTOgj4iAX2goIidPhiEpCSqNBu02G2pLStBcU4OBvj6022zosFrRbrN9LcT9PT2QJAlul1/H\nkslkUPv4QKXVQi6XQxsYCE8vL6h9fKA3GKDW6WBrbITNZEJjZSXy9+7FxeJi1BQWoqmmBr1dXcP2\ns18rKCYGSQsWQG8wIG/PHhTl5DjdSlVyFxdo/P0xed489HV349yRI2gbhvemaWxhiIluwrnDh2Fr\nasKiRx/F0jVr0NPZidydO9Hb3T2sX8n29fTAWFoKY2kpZHI5vH19kbZ0KcITE6ELDISXTgdNQACC\noqPh6u4OmVwOdw8PuHt6fu2d5O6ODgwODkLp5fW14/d2daG3u/vSf7u6MNDfD3NtLQqzsmAsLUV9\nWRkA4OS+feiw2Ybt57weH70e3n5+mH7vvYhJTcXpzz/H7tdeG9E5fF+eajUMycmITknB+ePHUZ6b\nO2YWaaHhwxAT3aS68+fx1fbtiE5NxaqXX0ZXWxsqTp6E3WzG4Ah8BSkNDsLe3IwDGzZc+TPfkBDE\nZWQgasoUqLVayF1d4aPXwycgAHK5HC5ubnB1c4PcxQUAYMWle9I9l6/mbU1NsJlMsJlMsJtM6O7o\nQENlJQq+/HLYf55vo9ZqkbJkCRJmzkT4xIloqKhATWGh0DndiFwuhy4oCFMXLUJvZydyt29H5zDc\nu6exRyYN95MV1zvpKFw1iOhakZMn44lNmxAaF4ct//M/OLJ1K8wXL6K/r0/01L7B288PuqAgeKhU\nV/6su6MD1QUFAmd1YzKZDAqVCjPuuw+3r1yJkPHjkbtzJ/a+8Qaqz54VPb3r8vTywpTMTKxetw6d\ndjtevu8+NFRWYsAJ/z6QGDfM7XAv3nE9cIIXqzk4hjxkMikgPFx6p7pa2t7XJ72SmyvFpqeLn9e3\nzFd21XDmBSY8VCrp/j/8QdpisUh/3LtX+vFTT0lhCQlOPWe9wSA9/OKL0od2u/Sr9eslhVIpfE4c\nzjVuhGtNE/1QkgSz0Yj/njUL2Zs3Qx8RgZUvvYRpd98Npbe36Nl906WV9K4MZ11gQuPvj6VPPIFl\nTz6JkqNH8dkbb+DL995DXVmZ084ZADR+fghLSEBfdzfKcnP5pDR9b/xqmmiIZDIZ9JGRmDBjBpY8\n/jhUPj7I3rwZBzduRHNtrejpjRp+oaGInz4dsWlpSLvzTpw7fBifvf466svL0d3ejkEnekXpWnqD\nAfNXrcLsBx5A0aFDeO/pp9FqNg/7O9U0utzo7wMf1iIaIkmS0FhVhXarFXIXF9y+ciVm3HcfgseP\nR/7evcjZskX0FJ1eWEICZq1YgdQ774SrmxuKjxzBzr/+FXWlpU55z/1a/qGhMCQnAwBKjh2DvblZ\n8IxoNGGIiRyk3WrFid27Yamvx8wVKzBx9mwExcQgICICOVu2oLmmhldI1wiKjb30pHFmJibOnYt2\nqxXVBQXI2rQJF4uKRsXn5R8WhomzZ8M/LAwlx46hKDtb9JRolGGIiRyo1WzGqf370WaxwGw0InH2\nbCxbswaagACUfvUVirKz0Wo2Y3BgQPRUhYucPBkzV6xA2IQJUHp7o7qgAHl79qCxstLpFur4NiFx\ncYibNg0yuRyF2dkO3RSEbg0MMdEwKM/LQ315OSry8pD56KOYtWIFkhcsQLbBgLITJ1BbXAy72XzL\nvdri6u6O4NhY+IWFIWP5csSkpcHdwwNnDhzAFxs3ovLUKdFTvCmagABEp6bCJyAANYWFqDx5UvSU\naBTiw1pEw8w3OBj3P/cc4jIyoDcYUF1QgOPbt6MoJwemqip0tLaO+SArlEpoAwMREB6OuQ89hEm3\n3w6ZTIb8vXtRU1SEkqNHUXXmjOhp3rSpixfj7qeegodajd2vvYas998XPSVyYnxYi0iQlro6vPX4\n44iaMgX3/v73iEpOxk///Gecz83Fka1bUXzkCMy1tehqaxtTX1m7urnBVaGAu6cnIidPxryVK5F4\n++3wUCrR09mJLzZuxOdvv42mUbpFoIdajchJk6DW6VBbXAzj+fOip0SjFK+IiUZYcmYmFv/yl0iY\nORMqHx80VVcje/NmHHjnHZjr6i695+vEr+p8F5lMBshkCJswAXEZGUjOzETCrFnw8vVFW0sLdv3t\nbzi6bRuaqqtHxRPRN7LwkUewbM0aNFVXY8err+LswYOip0RO7ka5ZYiJRpibQgGFUomkBQsw58EH\nEZueDg+lEp1tbWiqqcG+N99E7s6dI765gqOkLFmCcZGRmLliBcInToSrQoHGykp8uWkTTu/bh4aq\nKvR1d2NgYMCpF+j4NrrAQDy1dSsCY2Lw4fPP48CGDejv7RU9LXJyDDGRk1F6e8Pb3x/B48cjfelS\npC9fDrWPD2wmE9osFpQeO4YzBw6gLDcXloYG0dO9IaVGg4jEREQmJWHy/PkIjY+Hh0qFDpsNDRUV\nKD12DMVHj6K+rAwdNtuovgr+j3krV2LFn/4EY2kpPnn5Zb6yRN8LQ0zkpBRKJfxCQzHOYIB/WBjS\nly3DhJkz0dvVBUtdHezNzWipq8OFggKUHT+O+vJytLa0CLuaVGk00BsM8AsJQVBsLKJTUqCPjISX\nTgff4GBUnT4Na2MjTuzejQtnz8Le1IRWs3lMXDG6ursjfvp0PPTCCxgcGMBnr7+OU/v3o91qFT01\nGgUYYiInJ5PL4alWIzIpCUExMYicPBmGKVMQGhcHdw8PmI1GmC5cgO3yNoWmCxdw4cwZWBoa0GG1\not1mQ09np8Pm4+ruDrVWCy9fX+gjIgAA4ZMmQR8RAV1QENRaLXz0eqi0WtgaG9FYWYnKU6dQnp+P\ndosFdWVlaLdYHDYfZ+CpVuPHTz+NO3/zGxzasgW7XnsNxtJS0dOiUYIhJhplAqOjEZ6YiKCYGPiF\nhCAgIgJBMTHQBQVB6e2N5tpaXCwqgs1kQqfdjs7WVnTY7ehqa0OHzQZLfT0AoKu9/crT2H3d3ejr\n6QEAyF1c4Onl9Y0NKtRaLXyDg6EJCIBKo4Faq4VfWBgAIDQ+HgqlEm0WC9paWmCpr0dtSQmMJSVo\nqqlBdUEBOuz2Uf2w2Y14qNVImDULDzz/PPp7erDj1VdR8OWXo/ZePo08hphoFFP5+CB0wgRETZmC\ncVFR8NLpoPH3h39oKFRaLRRKJRQqFSBJaLNY0NrcjMbKSgCXlt4cuBzins7OK1fNLi4uUOt08PL1\nxdX/IjV6PQKjo+Gl06GnowNd7e1oNZsBAObaWphra9FYVQVLfT0aq6pQW1yMztbWEf08RppCqURM\nairu+t3vEBIfj20vvojcnTvRNsau+Gl4McREY0xIfDwSZsxAQGQkNP7+0Pj7Q+XjA7fL7+7+50rX\nTaG48m/OxdUVcjc3AIA0OIi+nh70dXd/7bi9XV3osNvRbrXC3twMa0MDjCUlAC5taGCurb1yVX2r\n0BsMWPLLX+KOxx7D/nfewYfPP8/7wnTTuKAH0RhjvPyV8NXcPT2hHTcOfiEhiJw8GcClK1xX10v/\n1D29veF5OdADfX2wNjTA2tj4tQe/mi5eRNXp0zBzC0cAgJuHB/Th4YiYNAk2kwnnDh/+xi8vREPB\nK2IiohuQyeWITUvDsieewPiMDOxYtw57/vGPUbErFDmfG/29kY/wPIiIRo3whATMW70aMWlpyPng\nA+xbv54RJodjiImIriMgIgLzVq1C4pw5OLlvH3b97W/8SpqGBUNMRHQdUcnJiElLQ4fNhsKsLNia\nmkRPicYohpiI6Brjb7sNGXffDZlcjuPbt6Pk6FHRU6IxjCEmIrpKTGoqFj7yCPyCg3F6/34c+/TT\nK4ujEA0HhpiI6CoxKSnwDQ5GQVYWDn34IerLykRPicY4vr5ERHSV1CVLoPLxQenx41dWJyNyBK6s\nRUREJBDfIyYiInJCDDEREZFADDEREZFADDEREZFADDEREZFADDEREZFADDEREZFADDEREZFADDER\nEZFADDEREZFADDEREZFADDEREZFADDEREZFADDEREZFADDEREZFADDEREZFADDEREZFADDEREZFA\nDDEREZFADDEREZFADDEREZFADDEREZFADDEREZFADDEREZFADDEREZFADDEREZFADDEREZFADDER\nEZFADDEREZFADDEREZFADDEREZFADDEREZFADDEREZFADDEREZFADDEREZFADDEREZFADDEREZFA\nDDEREZFADDEREZFADDEREZFADDEREZFADDEREZFADDEREZFADDEREZFADDEREZFADDEREZFADDER\nEZFADDEREZFADDEREZFADDEREZFADDEREZFADDEREZFADDEREZFADDEREZFADDEREZFADDEREZFA\nDDEREZFADDEREZFADDEREZFADDEREZFADDEREZFADDEREZFADDEREZFADDEREZFADDEREZFADDER\nEZFADDEREZFADDEREZFADDEREZFADDEREZFADDEREZFADDEREZFADDEREZFADDEREZFADDEREZFA\nDDEREZFADDEREZFADDEREZFADDEREZFADDEREZFADDEREZFADDEREZFADDEREZFADDEREZFADDER\nEZFADDEREZFADDEREZFADDEREZFADDEREZFADDEREZFADDEREZFADDEREZFADDEREZFADDEREZFA\nDDEREZFADDEREZFADDEREZFADDEREZFADDEREZFADDEREZFADDEREZFADDEREZFADDEREZFAriJO\nKkmSiNMSERE5HV4RExERCcQQExERCcQQExERCcQQExERCcQQExERCcQQExERCcQQExERCcQQExER\nCcQQExERCcQQExERCcQQExERCcQQExERCcQQExERCcQQExERCcQQExERCcQQExERCcQQExERCcQQ\nExERCfR/m/i+vthJCsgAAAAASUVORK5CYII=\n",
       "text": [
        "<matplotlib.figure.Figure at 0xcde3ecc>"
       ]
      }
     ],
     "prompt_number": 7
    },
    {
     "cell_type": "markdown",
     "metadata": {},
     "source": [
      "We have to step much farther back from the scene,\n",
      "to a distance of over 100\u00a0AU,\n",
      "to bring the orbits of the outer planets and Pluto into view.\n",
      "\n",
      "Jupiter and Saturn are the planets which, over one Earth year,\n",
      "have travelled the blue arcs shown close in against the bottom edge\n",
      "of the inner Solar System.\n",
      "Uranus and then Neptune have moved only small fractions\n",
      "of their huge orbits up near the top of the diagram,\n",
      "while minor planet Pluto sits at the right size of the diagram."
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "reposition_camera(110.0)"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "metadata": {},
       "output_type": "display_data",
       "png": "iVBORw0KGgoAAAANSUhEUgAAAeIAAAGFCAYAAAAhGPBBAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAEHlJREFUeJzt3VuMloWdx/HfOwyDIwwMg8NhBgYPIwpUZRXLtrLdHohb\nrUl7UbVGuzfbxO629aJXTZMmmzTdZi/a2JM92KQ1aWJN2qZrD1Z7sMUzdj0WRUQExBEGOQwDc2CA\n2YtmTZuWm1bmDzOfzyXvO3l+TwJ8877PvM/bGB8fHw8AUKKpegAATGVCDACFhBgACgkxABQSYgAo\nJMQAUEiIAaCQEANAISEGgELNFQdtNBoVhwWAMie6kaVXxABQSIgBoJAQA0AhIQaAQkIMAIWEGAAK\nCTEAFBJiACgkxMCUtmzNmvz3gw/m/NWrM6255B5HTHH+1gFT2ozW1vSuXp2P3357hg4ezCvPP5+7\nb701Ozdtqp7GFNEYP9E9t07mQd3iEjhFtC9YkLXXXpurP/axdJ1/fg4fOJAXH388T953X/73nnsE\nmTfNiXLrFTEwpR3YvTu/vuOOHNy7N/O6u9PZ05N//MAHcubs2dm1dasQc9IJMTDlDQ8OZv2ddyZJ\nOnt60r99e+YuWpThwcHiZUwF3poG+Cs6urpybGwsA3v2VE9hkjhRboUYACaAr0EEgFOQEANAISEG\ngEJCDACFhBgACgkxABQSYgAoJMQAUEiIAaCQEANAISEGgEJCDACFhBgACgkxABQSYgAoJMQAUEiI\nAaCQEANAISEGgEJCDACFhBgACgkxABQSYgAoJMQAUEiIAaCQEANAISEGgEJCDACFhBgACgkxABQS\nYgAoJMQAUEiIAaCQEANAISEGgEJCDACFhBgACgkxABQSYgAoJMQAUEiIAaCQEANAoebqATBVzT7r\nrMxsb8/w4GAO7N5dPQcoIsQwARqNRlrb2tJ72WVpO+ustC9YkLMWL86xo0ezcf36PHnffdUTgSJC\nDCdRa1tbFvX2Zk5nZ7ovuCCXXXVVOrq7s+Dss9NoasqLGzZkx3PPVc8ECgkxnERzFy7Muz784Vyy\nbl26enuzr68vo0NDeWXjxuzeti2bN2zIKxs3Vs8ECgkxnERHx8YyOjSUWe3tefH3v8+vvvOd7Nqy\nJTs2bszgvn3V84BTQGN8fHx8wg/aaEz0IaHMGTNnZulFF2Xzhg0ZP368eg5Q5ES5FWI4yRqNRpqa\nm3NsbKx6ClBIiAGg0Ily64YeAFBIiAGgkBADQCEhBoBCQgwAhYQYAAoJMQAUEmIAKCTEAFBIiAGg\nkBADQCEhBoBCQgwAhYQYAAoJMQAUEmIAKCTEAFBIiAGgkBADQCEhBoBCQgwAhYQYAAoJMQAUEmIA\nKCTEAFBIiAGgkBADQCEhBoBCQgwAhYQYAAoJMQAUEmIAKCTEAFBIiAGgkBADQCEhBoBCQgwAhYQY\nAAoJMQAUEmIAKCTEAFBIiAGgkBADQCEhBoBCQgwAhYQYAAoJMQAUEmIAKCTEAFBIiAGgkBADQCEh\nBoBCQgwAhYQYAAoJMQAUEmIAKCTEAFBIiAGgkBADQCEhBoBCQgwAhYQYAAoJMQAUEmIAKCTEAExJ\n//r5z+ffv/71nLNqVekOIQZgSlq2Zk3WvP/9WbJ8eVpaW8t2CDEAU9KrmzZl+owZ6V62LLPnzSvb\nIcQATEkvPfFERg4dSvcFF6R94cKyHUIMwJR0eGAgA/39SZLm6dPLdggxAFPS2OhoBvbsycG9ezM2\nOlq2Q4gBmJJe37Ejs9rbc/7q1eno6irb0Vx2ZAAotPWpp3Lv7beno6sr+197rWxHY3x8fHzCD9po\nTPQhAeAvTJ8xI01NTTk6NpZjR4+e1GOdKLdCDAAT4ES5dY0YAAoJMQAUEmIAKCTEAFBIiAGgkBAD\nQCEhBoBCQgwAhYQYAAoJMQAUEmIAKCTEAFBIiAGgkBADQCEhBoBCQgwAhYQYAAoJMQAUEmIAKCTE\nAFBIiAGgkBADQCEhBoBCQgwAhYQYAAoJMQAUEmIAKCTEAFBIiAGgUHP1AOAkaTQyq709K9/xjnR0\ndaV5+vQ3HhrYsyfPP/RQ9uzYUTgQSJLG+Pj4+IQftNGY6EPClNE0bVrmzJ+fi9/97vSsWJHlV1yR\nuQsXZtqfhPjg669n84YN2fXSS3n0xz/O3p07c+zo0cLVMPmdKLdeEcMkMv2MMzK/pyeXrFuXaz7x\niXR0deW1l15Kc0tL5i5cmJbW1hwZHs4ZM2dm0XnnpaW1Ne0LFmT9nXfm1c2bMzYyUn0KMOUIMUwi\ns+fNy6orr8w7b7opo8PD2fTII9m5aVN6Vq7MmbNnp6W1NaNDQ9n32mtpNDWlZcaMvPfmm9Pc0pK7\nb701r7/ySvUpwJQjxDAJNDU15Yy2tizq7c3iCy/Mwddfz12f/WzOWrIk//aFL6Szp+eN57bNm5e2\nefOyr68vP7vttpy7alUuWbcuf/jd7zJ6+HCGDh70NjVMINeIYRKYPW9e1l5/fVa/733Z9swzeezH\nP84rmzblc/ffn6UrV6a5peUvf2h8PIP79+cLN96Y1VdfnbXXXZf7vv3t3Hv77dmzffvEnwRMcifK\nrY8vwSRwzqpVueyqq3LW4sUZ6O/PvMWL81+//W1GDh3K0bGxv/ozhw4cyFP33ZdrbrklfVu2ZHhw\nMB/45Cdz+TXXTPB6mNq8NQ2nuYXnnpsrrr02iy+8MI/dfXdeePTRXPvpT6dnxYocHRvLrq1bs/3Z\nZzNr7tycc8kl6ejqypGRkQzu3585Cxbk/NWrc87FF+e5Bx9MGo1c9M53ZusTT2TTI49UnxpMCUIM\np7kD/f35xTe/mQfvuitDg4PpXrYsK9auTXNLS5pbWrKvry/NLS3p7OnJrLlzkyRjIyMZPngw565a\nlda2trS2tWXJ8uWZOWdOupctS9f55wsxTBAhhtPcyKFD2frkk0mSOfPn58K3vS0z29vfeHxsZCTd\nF1yQfX19efiHP8zRsbFcsGZNLliz5o0wJ0lzS0uapk3L7M7OzO7snPDzgKlKiGGS29fXl56VKzM6\nPJxnfvObDO7dm84lS/4itgN79qS1rS1tHR1/FnLg5BJimESODA9ncO/eP/uzsSNH8urmzZm7YEGu\n+OAHc2j//izq7c2eHTv+7GNNA/396Vi0KEfHxnLsBL/gBbz5hBgmkeHBwezaujWHDxzIzPb2HBkZ\nyfmXX56dmzZl/tKled/HP54kObB7d15++ukcP3Ys888+O0MDA9n98stZeO65Sf4YZWBi+PgSTDKD\n+/bliXvvzUB/fwb37cuyyy/PirVrM2f+/Deec+acOeno7s7jP/1pDuzenT+sX58dGzf+8eNOR47k\n6JEjhWcAU4sQwyTTt3lzbvvoR/OTL385rbNmpam5Od3LlmXOn1wTbjnjjCxYujSLly/Pz7761fz8\na19L+/z5mdXRkde2bMnul18uPAOYWoQYJqGhgYHcfeut2fjAAxkdGvqrz2luaUnbvHl5+Ic/zOjQ\nUP7phhuSRiPrv//9/GH9+gleDFOXEMMkND4+ntGhodx28835/c9/nsMDA3/2+P5du/Lb730vd332\ns9n98su55D3vSaOpKY/86Ed58fHHc/zYsaLlMPW41zRMckuWL8/szs4sOOectHV0ZHDfvry2ZUsG\n+vsz0N+fwwMDWdTbm0W9vdm9dWv6t2/P2Oho9WyYdE6UWyGGKeL/vwbxyPBwhg4e/IvHW1pbc2xs\nzDcvwUkixABQyLcvAcApSIgBoJAQA0AhIYYijaamNJr8E4Spzv8CUGTmnDmZOWdO9QygmBBDgUW9\nvfnw5z6X/7znnlxzyy3Vc4BCQgwF9vX15YVHH83rO3dm7sKFaevoqJ4EFBFiKDA6NJTnH3oozz3w\nQNrmzs0/XHll9SSgiBBDkddffTU7N23K9NbWXLJuXTp7eqonAQWEGIqMjYxk5wsvZNvTT6d94cKs\nve666klAASGGQv3btmXDT36Sbc88k7XXX5/LrroqzS0t1bOACSTEUOy1LVuy/s47Mzo0lP/4xjdy\n9kUXZfqMGdWzgAkixHAK6N+2Lf/zxS+mpbU113/mM1lw7rlpmjatehYwAYQYTgHDg4N55te/zh2f\n+lQuete7csW116ajq6t6FjABmqsHAH80cvhwHvrBD7L87W/Pv3zkI2lubs6vvvvd7N66tXoacBL5\nPmI4hTQajZy1ZEmu/fSnc+Hb355Go5GffuUrufdb36qeBvydfB8xnAbGx8ezZ8eO3P2lL/3xZh8d\nHXnnjTdmzfvfXz0NOEm8NQ2noJ3PP597b789mx55JBe+7W1Z1NtbPQk4Sbw1Dae48y69NDNmzsxz\nDzxQPeWU02g00tHdnaampgz09+fIyEj1JDihE+VWiIHTUtO0aWlfsCDv+NCHcqC/P0//6lfZv2tX\n9Sw4IdeIgUmj0dSUme3tufjd784/33RTOnt6Mq3ZlTZOT0IMnHZmtLamZ8WKrFq3Lvv6+vLUL3+Z\ngT17qmfB30SIgdPOnPnzc9nVV+f48ePZ/Nhj6d+2LWOjo9Wz4G8ixMBpZcaZZ2bJihW57Kqrcuzo\n0dxz2205uHdv9Sz4m7moApw2zpg1K0tXrsx5l16azRs25N5vfSuHBwYyfvx49TT4mwkxcNpY9ta3\n5h033JAjIyO5/447smPjxhw7dqx6FvxdvDUNnBZ6Vq7Mpe99b+Z0dmbTww/npaeeypHh4WTiP4EJ\nbyohBk55888+O5dfc02WvuUteeX55/Ps/fdn9PDh6lnwphBi4JQ3f+nSzOvuTt+LL+aJX/zCjTuY\nVFwjBk55h/bvz7P335/+7dvz0hNPVM+BN5VbXALABHCLSwA4BQkxABQSYgAoJMQAUEiIAaCQEANA\nISEGgEJCDACFhBgACgkxABQSYgAoJMQAUEiIAaCQEANAISEGgEJCDACFhBgACgkxABQSYgAoJMQA\nUEiIAaCQEANAISEGgEJCDACFhBgACgkxABQSYgAoJMQAUEiIAaCQEANAISEGgEJCDACFhBgACgkx\nABQSYgAoJMQAUEiIAaCQEANAISEGgEJCDACFhBgACgkxABQSYgAoJMQAUEiIAaCQEANAISEGgEJC\nDACFhBgACgkxABQSYgAoJMQAUEiIAaCQEANAISEGgEJCDACFhBgACgkxABQSYgAoJMQAUEiIAaCQ\nEANAISEGgEJCDACFhBgACgkxABQSYgAoJMQAUEiIAaCQEANAISEGgEJCDACFhBgACgkxABQSYgAo\nJMQAUEiIAaCQEANAISEGgEJCDACFhBgACgkxABQSYgAoJMQAUEiIAaCQEANAISEGgEJCDACFhBgA\nCgkxABQSYgAoJMQAUEiIAaCQEANAISEGgEJCDACFhBgACgkxABQSYgAoJMQAUEiIAaCQEANAISEG\ngEJCDACFhBgACgkxABQSYgAoJMQAUEiIAaCQEANAISEGgEJCDACFhBgACgkxABQSYgAoJMQAUEiI\nAaCQEANAISEGgEJCDACFhBgACgkxABQSYgAoJMQAUEiIAaCQEANAISEGgEJCDACFhBgACgkxABQS\nYgAoJMQAUEiIAaCQEANAISEGgEJCDACFhBgACgkxABQSYgAoJMQAUKi54qDj4+MVhwWAU45XxABQ\nSIgBoJAQA0AhIQaAQkIMAIWEGAAKCTEAFBJiACgkxABQSIgBoJAQA0AhIQaAQkIMAIWEGAAKCTEA\nFBJiACgkxABQSIgBoND/ASnowDwRUIxLAAAAAElFTkSuQmCC\n",
       "text": [
        "<matplotlib.figure.Figure at 0xc5432ac>"
       ]
      }
     ],
     "prompt_number": 8
    },
    {
     "cell_type": "markdown",
     "metadata": {},
     "source": [
      "The scale of the Solar System is many times larger\n",
      "than the diameters of any of the objects within it.\n",
      "Even the Sun, which is by far the largest object in the system,\n",
      "cannot easily be seen\n",
      "unless we move in close to the orbit of Mercury:"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "# Draw a sphere for the Sun\n",
      "\n",
      "sun_r = 695500.0 / AU_KM  # km\n",
      "sphere = mlab.points3d(0, 0, 0, name='Sun',\n",
      "      scale_mode='none', scale_factor=sun_r * 2.0,\n",
      "      color=yellow, resolution=50)\n",
      "\n",
      "reposition_camera(distance=1.3)"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "metadata": {},
       "output_type": "display_data",
       "png": "iVBORw0KGgoAAAANSUhEUgAAAeIAAAGFCAYAAAAhGPBBAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3XdglfXB9vHvWdl7kEkYIRD23lNkC4IDFdS6bavV2j52\nvl3a9fRtq6it1i32URQRlT1kyF5hhiQkhEyyyZ4nJznvH9C89qkoKDn3SXJ9/tHDuc99X4eEc517\n/X4mp9PpRERERAxhNjqAiIhIV6YiFhERMZCKWERExEAqYhEREQOpiEVERAykIhYRETGQilhERMRA\nKmIREREDqYhFREQMZDVioyaTyYjNinR4FpuN8Lg4rr/nHkKio+kxeDC+QUFYrFa8/fzw9vfHYrNh\nsf77P+3qCxdobmwkPy2NlN27KTp3jvrqanJTUgCoLC6m1eEAwNHcTEtzs8vf29dlsdn41yeKh48P\nNg8PALz8/OjWowce3t70GjqU4KgoBl93Hf4hIfgGBuLl5weA0+nEYbdjb2igpbmZ+poaygsKsDc0\nkJeayrHNmzl3/Dgtzc1UlZYa9C6lM7jcQJYmI4a4VBGLXBmfgAD6jh3L4OuuIyo+Hr+QEDx9fC4W\njJcXPoGB/1a69sZGqkpLKc3N5eDHHwNwPj2d+qoqWltbaaipobqsjKa6OlpbWqivrjbqrbU7s8WC\n1WbDZDbjGxSEh5cXAeHhWD088A8JwcffH9/gYKL69CE0JoYBkycTGBaGyWym2W7H2dJCfU0NlUVF\n1FVV4Whqojg7m9O7dnF61y7KCwpwdKAvLGI8FbGImzNbLPhdKoY+I0cS068f/qGhhHXvTmTv3viH\nhODh7d22fG1FBa0tLZw9coRzx45RXlBAXVUVjXV11FVWkndpb7e6rIyWS3u7cpHN0xOLzYaHt3fb\nHnK3nj3x8vXFNyiIbj17EhwRQXTfvkQnJODl53fxi0x1NcVZWRRnZ9NQXU1tZSXn09LIOHyY82fO\n0FBba/RbEzemIhZxQzZPT6ITEhg4ZQqhMTF4+/sTFhtLbP/+bYdVAZrtduqrqqgoKiL9wAFK8/Ko\nr67G2dpKXmoq58+cobq0FHtjo8HvqOPz9vcnOCoK/5AQwuPiiOrTh+CICPxCQujWowdxgwbhGxiI\nyWymsb6e0uxsclNSKMzIoCw/n6Jz5zizf//Fn4/m1JHPURGLuAlPHx96DB5MdJ8+BISFETdwICPm\nzCGse/d/W66poYG806c5n55OdVkZtRUVXMjP5+SOHRSfO6cPeRexeXnhFxREYHg4UX36ED9qFOGx\nsXgHBBAQFkZ4XBxBEREAVJWUkH/mDCe3baO2spKq4mIqi4spycmhqrQUh91u8LsRI6mIRQxiMpvx\n8vWl+4AB+AUHExAayrhFixg0bRoBYWEAtLa20lhbS3FWFrUVFdgbGqgpL+fI+vWc3LaNypISg9+F\nfF5AaCgB4eFExseTMHo0PQYNIjgyksBu3Qjr3h0PLy8ACtLTyU9LI/3gQbJOnuRCfj6VJSWUFxQY\n/A7ECCpiERcyWyz4BAYSGB6Oj78/4T16cNOTT9J7+HBsnp4AtLa0UFtZSXVpKfXV1RSdO8dn77xD\n9okTlBcVdagrl7s6Lz8/Inr2pP/EiYy/5RbCYmKwenjg6eODT2AgHt7eVBQWkn3yJKd37eLwunXU\nVVZSW1lJU12djm50ESpiEReweXq2XQA0cu5c5j3yCN0HDGh7vrW1FXt9Pc12OzXl5Rz8+GO2vPoq\n59PTDUwt15LFaiUkOpqQ6Gj6jBzJyHnz6DlkCL5BQVhtNqw2G3WVlRz45BMOfvIJ6YcO0VBTQ3NT\nE60Oh0q5E1MRi7jA5NtvZ8H3v0/8iBGYzGYsVuu//b5fOH+eD/7wBw6tXUtVSQmtra368O2ETCYT\nmEyYzWZMZjOx/frRb9w4hs6YweBp0wgMD7/4s29poTQvj8Pr1rFv1SqyTpygoabG6PjSTlTEIu3A\nbLEQGhPDbb/4BYnjxhEYEYFvYGDb4Wd7YyM5p06x+s9/5tyxY7Q0N1NbUUHjpft4pWuw2mzYvLzw\n9PEhJCqKPqNGkTBqFCPmzSOoWzeaGxtprKujsriYjCNHSNqwgZQ9e6guKzM6ulxDKmKRayiiVy8m\nLV7MqBtuwOblRVR8PL5BQZgtFgBS9uxh14oVZJ88SUNtLcWXRrISsdhs+AQE4BMQQGC3bgyfNYuJ\nt95KdN++WKxW6quqKM3NJTs5mfQDB9i3ahVVpaU6atIJqIhFviFvf38GT5vG4GnTiO7Xj+6JiUTG\nx7c9X1tRwaG1a0k/dIicU6fIS0nRHo18pfC4OGL69aNbjx70HDKEvmPHEp2QAEBlURH5aWmU5udz\nbPNm0vbto7aiwuDE8nWpiEW+Jp+AAIbOmEHC6NEkjh9Pz6FD8QsKAqChpoad77xDaU4OddXVpB84\nwPkzZ2isqzM4tXQ0nj4+baUcFR9PdN++9J84kR4DB1JfU0Pm0aOkHzrEuWPHyD55kpLsbJrq642O\nLVdBRSxyFcxmM/6hoUy67Tb8QkIYMWcOPQcPxtvfH4fdTklODmePHKEoM5Ntb71FcXa2zvnKNWPz\n8CAsLo7B113HgEmTiO7Th6iEBCxWK0WZmaQdOEDavn1knThBaU6Ovvh1ECpikSvg4eVFWFwcvYcN\nI6x7d+757/9um1ShoaaGvNRUCjIyOHfsGPtXr6Y4K8vgxNKZmS0WAsPD6T9hAokTJ9J72DBiExPx\n8PGhJCuL1H37SD9wgOxTpyjIyMDe0KBzyW5MRSzyJTy8vAiJjiY2MZGRc+dyw/e+1/ZcRVERFUVF\nnE9LY8trr3Fi2zYDk0pX1nfMGMbfcgt9x4whNCaGoIgInE4nx7ds4bN33yUvNZXygoKLs2u1thod\nV/4XFbHIF7DabPgGBRGTmMj8Rx9l0u23A9DicFBXUUGz3c7W119n8yuvcOH8eYPTilzkExDAwMmT\nmX7PPQyaOhVPHx+8/Pw4vG4dW19/nbNJSVSXldHc1IRThew2VMQiX6DPyJE8/MIL9Bs37t9+Lwsz\nM1l2zz2cOXDg4p6FDveJG/L08aF7//5MWbKEhT/8ISaTCafTyakdO9j08ssc37pVV1m7ERWxyCUW\nq5Ww2Fh+sXYtgd264RcUhNXDA4AT27ax6r//m9zkZGrKyzVbjrg1k8mExWbD28+P6IQERs+fz8wH\nHsAvJITGujqyjh9n13vvcfDjj6kqLTU6bpenIpYuz8PLiz4jR3LzT35CePfuxA0a1HYh1mfvvsve\nVavIS0mhLC9Pt4VIx2IyYfPwwD8khLhBgxg4ZQrjb7qJiF69qCot5fyZMxxet449K1dSVVKiC7oM\noiKWLi1x/HimLFlC/IgR9B4xAk9vbwA2vvQS6YcOkXn0KIVnz6qApcPz8PIiKCKCnkOHkjhuHEOm\nT6fn0KGU5eaSn5bGnpUrObR2rUZ6M4CKWLqk/hMnMnDyZAZOnkzfcePwDwnB3tDAyR07OLVjB0kb\nN1KUmYm9sdHoqCLXlM3Li5CoKHoPG8aAyZPpM2oU/SdM4NyxY5zcsYOzSUmk7tmjixBdSEUsXYq3\nnx9jFy1i/E030X/CBIIiI2mqryc/LY3D69ZxetcuTu/erXPA0ulZPTwIjoykz6hRjJo7lz6jRxMS\nFUVRVhYnt23j0Nq1ZB0/ri+jLqAili7BYrMxYOJEYvr146Yf/Yio+Hgcdjv5aWmk7d9P2v79bF++\n3OiYIi73r8FBRt1wA0OmT2fg5Ml4+flxfOtW9q9eTfapUxRnZWFvaDA6aqelIpZOzcvXl249exIe\nF8fSp56i55Ah2Dw9Kc3N5Xx6Ojv++U92vP220TFF3EJAWBjzH3+cYTNnEtmrF421tRzbsoX9q1eT\nl5bGhfx8oyN2Sipi6ZRMZjM+AQEkjh/Pwh/+kGEzZgBQX1VFXVUVa5YtY9d771FRWGhwUhH302fk\nSCYvWcLIOXMIioigvLCQw2vXsumVV6gsKsJht+sK62tIRSydjtliwS84mNkPP8zdv/89cPEXvaW5\nmfV/+xtrnnuO0txcg1OKuL8+I0dy4xNPMPn223EC5QUF/GXJErJOnND41deQilg6nbnf/S6zHnqI\nmL598fL1xdHcTH5aGk/NnUtdVRX2hgbNiCRyBSxWK/4hIcSPHMnshx9mzIIF1FdXs/2tt9jw0ksU\nZGQYHbFTUBFLpxEQFsbdv/89o264gcDwcKweHuQkJ/PBH/9I5tGjFJw5o2/wIlfJZDbj6e1NWPfu\nTFy8mOvuvhtvPz+Kzp3j8Pr1HPj4Y/JTU42O2aGpiKXDM1ss3PKTn9Bv7FgSJ0wgICyMxro6jm7a\nxJply8hJTqaustLomCIdmtliISQqioQxYxh6/fUMnz0bi8XC0c2b2b58OWn79xsdscNSEUuH5h8a\nyrf+8AeGzZxJWGwsFpuN41u3cmzLFk7v3k36wYNGRxTpVP51J0LiuHFMvuMOInr3Jjc5mZM7d3Jq\nxw6yjh83OmKHoyKWDikgPJx+Y8cyeNo0Fv3XfwHQUFPDwU8+Yec773DmwAHtBYu0Iy9fXwZOmcLY\nhQtJnDABL19fTu7YwbY33yT94EFaHA6jI3YYKmLpcKITEhhz442MXbiQgZMnt03vlpuSwoYXX9T5\nKhEXik5IYNQNNzBmwQLiBg0iZdcu9n30Eae2b6eqtFQXRl4BFbF0GFabjag+fZj54IPMfOABfPz9\nqa+u5szBg7z/9NNkHjum0X9EDODh5UX8yJHc+IMfED98OJ4+Pnz05z9zaudOCjIyaKipMTqiW1MR\nS4fgHxpKZK9e3PzjHzNx8WLsDQ1cKCjgzIEDPHPXXUbHExHANyiI6++9l1kPPUREz54c27KFLa++\nStr+/dRXV2vv+DJUxOL2PH19ufnJJ7n+3nvp1rMnDrud1L17+eCPf+T41q1GxxORzzGZTITGxPDo\nK6+QOGECFYWFbHn1VXa/9x4VJSW06tzxf1ARi9v71YYNDJgwAW9/f0xmMxtefJG3fvxjmpuadEGI\niBsymc1E9urF9ffey+Q77sDLz49jW7aw9fXXOb1rl9Hx3I6KWNyWb1AQv9+xg5i+ffHw9qa2spK1\ny5ax6eWXqSwuNjqeiHwJi9WKb1AQQ6ZPZ/a3v03CqFFkJiWx4cUX2btqldHx3IqKWNyOt78/iRMm\nsPinP2XglCmYzGZ2rVjBwU8+IUUTlot0KL5BQfQYNIipS5cy5sYbKcvLY/vy5Wz8xz+MjuY2VMTi\nVmL69mXCrbcycu5cBkyaBMCq//5vDn78MbkpKbr6UqQDsnp4ENWnDyNmz2bE7NmExMSQtGkTG//+\nd8ry87v8KabL1a3VxTlE6DduHNO/9S2Gz55NZO/eNNXX89Ff/sLW11+nvKCgy/9jFemoHHY7eSkp\nVBUXU1FYyLxHH2X2gw9iNps5sn49mUePagCeL6A9YnGpgVOmcMOjjzJq3jxsXl5UFBayZ+VKlv/s\nZ7Q0NxsdT0SuEd+gIIbNmMHYRYtIGD2a3NOn2f3ee5zauZOqkhKj4xlCh6bFcAmjR3Pn008zYs4c\nasrLyU1O5tTOnbz7618bHU1E2oHZbMY/LIylv/kNg6ZNoywvj93vvcfRTZsoLyw0Op7LqYjFMFab\njcj4eB59+WUGTplCdWkpRzZsYONLL3FGkzWIdHpB3boxcfFipt15Jyazmd0rV7L3gw+oKCzsUqei\nVMTieiYTVpuN2MREfrJyJTH9+tFQW8u655/no7/8hdqKCqMTiogLjVu4kAXf/z49hgzh8Lp1rPzt\nbynOyqK1tdXoaC6hi7XE5WyeniSOH8+P3nuPwPBwAN7+6U/Z9PLLtGgIPJEu59DatbS2tHDjE08w\n/qab8PD25vUnnuiSh6k/T3vE0i5CoqKY+eCDzP3OdwiKjMRsNvO3hx5iz8qV1FdXGx1PRAzi6eND\n9/79mXH//Uy+4w5S9+1jy6uvkrRhQ6c/TK1D0+IysYmJzP3Od5hw662ExsQA8MKDD3Lgo4+oKS83\nOJ2IGM3m6Ulk796MW7SImQ89RGNNDe/++tcc//RTGmtrjY7XblTE4hI9Bg9m5v33M+m22wiKiKCq\ntJR3f/1r9n7wgc4Ji0gbi9VKeFwcE2+7jYVPPEFJTg5bX3uNA5980mlvb1IRS7vrOXQo0+++m4m3\n3UZgeDh5qalsfOkldv7P/2j+YBH5DxarleDISMYuWsTUpUvB6WT78uUcWru2U543VhFLuxo0dSpT\nly5lxOzZePn5kbZvH3s//JDty5cbHU1E3JzFZmPcokXMuO8+rDYbez/4gEPr1lFeUGB0tGtKRSzt\npu/YsSz++c8ZNmMG9sZGTm7fzsaXXuLk9u1GRxORDmTknDnc8L3vYbHZ2Pfhhxxas4aKoiKjY10z\nun1J2kX3/v258+mnGTh5Mg67nVM7dqiEReRrSdq0CUwmbnj0USbffjtms5n9H33U6adD1R6xfC0m\nkwkPb2+e2rKFhFGjaHE4SNqwgQ0vvsipnTuNjiciHdiI2bOZ/9hjhPfowa4VK1j/t791itsetUcs\n15TVw4PxN99MdEICNk9PPn3zTdb/7W/knj5tdDQR6eCObt6Mo7mZRT/8IXO+/W0sVivvPf200bHa\njYpYrlp4XBwPPfccAyZNwi8khD0rV7J9+XLOp6cbHU1EOonUvXsxmUzMe+QRpt15JxWFhWx+9VWj\nY7ULHZqWqxIVH8/8xx9nxn334e3vz56VK1n7/POcO3aMpvp6o+OJSCfi5etLv3HjmPPtb9N7+HCO\nb93K2z//eYed01iHpuUbi+jViylLlzLx1luxeXqy78MPVcIi0m4a6+pIP3QITCbmP/ooE269lab6\nej559lkunD9vdLxrRkUsVyQqPp5Jt93G1KVL8Q0KInXvXtYsW8bZpCSaGxuNjicinVRDTQ1n9u/H\nZDJh9fRk/C23UFVaymfvvENZfr7R8a4JFbF8pZDoaMYsXMi0u+4iqk8fCtLTWfv886Tt30+rZlES\nkXbWWFfH6V27wOlk0X/9FxNuuYX6qioOfPIJFZ1gBC4VsXwpm5cXw2bOZModd9B9wABKcnLY/vbb\nHPj4Y6OjiUgX0tzUxPFPP8Xq4cH8xx5j3KJFOOx2jmzY0OEH/TAbHUDcW8KoUUy76y4SRo+moqiI\nvR98wId/+pPRsUSkizqyYQPr//53nMDkO+5g5Lx5BISGGh3rG1ERy2WFxcay9KmnGDZjBrUVFWx7\n6y3e/NGPjI4lIl3c4XXrWP/CCwDMvP9+xixciLe/v8Gpvj4VsXwhq4cHv1i7liHTp9PicLBm2TLe\n/tnPjI4lIgLA4fXr+fTNN7HYbCx4/HEm3norZnPHrDSdI5b/YLFaee7YMaL79gVg+U9/ysZ//MPg\nVCIi/27/Rx/R2tLC/MceY94jj1BRXEzShg1Gx7pqHfPrg7SbwG7d+MWaNUTGx2OxWlnx1FPsXbmS\npro6o6OJiPyb5sZGUvbs4fC6dQRHRbH0N79h5Ny5WKwdax9TRSxtwuPiuONXv2LQtGlYPTw4umkT\nRzdt6pQTdItI51BVUsKBjz5i21tv0X3AAOZ85zuExsRgtliMjnbFVMQCQEhUFONuuomJt96Kp7c3\ntZWVbHvrLfJTU2lxOIyOJyLyhVocDoqysvjsnXfY9I9/MGDSJG756U8JCAszOtoV61j779Juovv2\nZfrddxMUEUFraytrn3uOkzt2UFdVZXQ0EZEv1dLcTEFGBptffZXgqCgm3HILVSUlbHntNcry8oyO\n95VUxELcgAFMuu024keOpLG2lt3vv8+aZ5/tFPN/ikjX0OJwUJyVxeo//5mInj2Z/fDD1FdXs/u9\n99x+XGodmu7iovv25YZLVxzWVVZy4OOPeeHBB1XCItLhOOx2ck+fZu3zz1NZVMTMBx5gxJw5+AUF\nGR3tS6mIuzC/4GAWPvEEc7/zHZqbmkjdt49n7r7b6FgiIl9bS3Mzu99/n23Ll2O12Zh8220MnDIF\nD29vo6Ndlg5Nd2FTli5l0LRpABSkp7P7vfeMDSQico2sWbaMHoMHM+Hmm7m+vp6aigpS9+y57JzA\nRlIRd2H9xo4lpl8/AM6np7NLRSwinch7Tz1FQGgoo+fPp7qsjPzUVKrLyoyO9R9UxF3UD5YvZ+zC\nhZjNZna8/TZv//zntDQ3Gx1LROSaKS8o4OO//hWzxcLIefOoKCrinV/+0uhY/0HniLugJ956i9EL\nFuATEMDmV19l5R/+4PZXFYqIXK0Wh4PMo0fZs3IlpdnZjL/5ZhY8/rjRsf6DirgLsdpsjJ4/n+Gz\nZuEXHMypHTvY9+GHnD9zxuhoIiLtorGujhOffsru99/H0dTElCVLGDF7ttGx/o2KuAuxeXmx4LHH\n8L10Kf/RzZs5d/SowalERNpXeWEhRzdv5sS2bYTHxTHv0UfpPWIEVpvN6GiAirjL8AsJYdZDDzFo\n2jQ8vL05uGYNJ7dvp6q01OhoIiLtrjQnh8Pr1nE2KYkh06cz47778PT1NToWoCLuEnyDghg1dy4P\n/PWvWD08yDx6lBW/+Q0Zhw8bHU1ExCXsjY2cO3aMT994g/KCAibddhsDJk3Cy8/P6Ggq4q4gNjGR\nW3/2s7bH7//2txSePWtgIhER16uvriZ13z62vPoqFouFO3/7W+JHjMDm5WVoLhVxJ+fh5UV4XBxx\nAwfidDqpLisj68QJGmpqjI4mIuJyVSUlbH71VXa99x7RCQnMfvBBYhISMJlMhmVSEXdyk2+/nf96\n5x0A6ior+XZCAiXZ2caGEhExUGNtLZ888wyNtbVMufNO+o4da+ghahVxJzbuppuY893vYrFaKcvL\n428PPURDTY1bDvEmIuIqLQ4Hpbm5rPjNbyg/f55FTz7JsFmzDMujIu7Eovv0oefgwQA01NZyfOtW\nWltbDU4lImK8FoeDvR9+yKG1a/Hy9WXM/PnEjxhhSBYNcdlJzXzgAabdfTeePj5knzzJ8p/+VFMb\nioh8TnVpKbtWrKD38OEMnzWLmgsXaKyrc/kgR9oj7oQGTJrE5Ntvp+fgweSnpbHltddI2rjR6Fgi\nIm7n3LFjfPbuu9SWlzNizhwGX5qRzpVUxJ3QuEWL6DV0KAA5ycns+/BDgxOJiLinxro69q1aRdKm\nTXj7+zNk+nRiExNdmkFF3MkMvf56Bk2dSmC3buSlpnJqxw7KCwqMjiUi4rYqioo4tnkzeSkp9B0z\nhuvuuouQ6GiXbV9F3Il4+/tz/1//Sp9Ro6itqODTN95gw4svGh1LRMTtpR86xPGtW7F6eDDpjjsY\ns2CBy7atIu4kTCYTg6ZOxScgAIAT27aRtm+fwalERDqG+upqdr//PptefpmAsDCmLFmCt78/uGCg\nDxVxJ2G2WnnstdeI6NULnE5OfPopqSpiEZErduH8eU5s20b6gQPEDRzIDd/7HjYPj3bfroq4E/AL\nDub+P//54rc3YNWf/sTxrVsNTiUi0vHkpaSwd9Uq/EJCmP/YY3Tr2ROLtX3v9FURdwJevr5MvfNO\nPC4NXH7u2DHK8vIMTiUi0vHUV1WRfugQ+1evJjA8nIeff57gyMh23aaKuIML696dm370IwLCwjCZ\nzbz5ox9xevduHM3NRkcTEelwWltbKcrMZOvrr1NbUUH/iRMZfN11+AUHt9s2VcQdXGC3bky67TYA\naisqOLZlCxWFhQanEhHpuBrr6shMSmLbW29htdmYvGQJoTEx7bY9FXEHFtOvHzPuu4/gyEicra1s\nePFFKouLjY4lItLh1VVWsv3tt8k6cYKEUaMYOHkygeHh7bItFXEH1mPwYG549FEcdjtZJ07wybPP\nqohFRK4BR3MzBenpbH/7bXA6mbJkCb2HD8dqs13zbamIO6jQmBjiBgwAoKGmhjXLltFUX29wKhGR\nzqOluZm9K1eSl5pKn9GjGTZjBqGxsdd8OyriDshiszFlyRKWPvUUjuZmSrKz2f7229gbGoyOJiLS\naTidTipLStj62muU5uYy/uabGTZzJrZLd6hcKyriDiiiZ08ievcGoDQnh/d//3uDE4mIdF47/ud/\nKMnOJqJ3b4bPnk388OHXdP0q4g5o2l13MeO++wCouXCBpA0bDE4kItK57Vu1iryUFAZMmsSQ6dOv\n6bpVxB1MRK9eRPTsiYeXF7mnT7Px5ZdpbmoyOpaISKe2d9UqUvfuxebpyeDrrmP4rFnXbN0q4g5m\n+j33MHTGDABKsrM5sn69wYlERDq/2ooK9n74IRmHDtF9wAAGTplyzdatIu5AwmJjiR8xgpDoaPJS\nU0nauJGqkhKjY4mIdAkZhw6RefQoHt7eJI4fz4BJk67JelXEHciUpUvbblk6d+yY9oZFRFyorrKS\nlD17yDp+nJ6DBzNlyZJrsl4VcQfh6ePDmPnziYyPpzgri/SDBynOzjY6lohIl3Lm4EGObd6M2Wpl\n4JQpxA0c+I3XqSLuIIbNnElAWBgAJ7Zt4+jmzQYnEhHpeqpLS0m/dIg6ODKS2Q899I3XqSLuIO77\n85+J7d+f5qYm8lJSOH/mjNGRRES6pJzkZI5s2ICXnx8TFy8mMDwcs/nr16mK2M2ZTCa8/f3bfsh5\nqamU5uYanEpEpOuqKikh/eBBso4fx9vfn4U/+AFefn5fe30qYjfn7e/P77ZvJ6x7dwDWPf88Bz/5\nxOBUIiJdW2FGBp+9+y7efn4s/OEP8fL1/drrUhG7ObPFQkxCApZLM3402+20OBwGpxIR6dqqL1wg\ndd8+zp04gc3Dgxn33UdQt25fa10qYjcWEBbG/Mcew8PHB5PJxBtPPsnJbduMjiUi0uW1trRQlpvL\nzn/+E0wmZj/8MEGRkV9rXSpiN+YXEsL199yDxWoFIDMpiYqiIoNTiYgIXLyv+MS2bVSXlRHWvTvD\nZ80iOCrqqtejInZTPgEB9Bo6lIjevTGZTBzdvJmq0lKjY4mIyCWO5mZKc3M5uGYN9sZGxixYQLe4\nuKtej4rYTYX36MG0u+5qe7z19dcpy8szMJGIiPxvjXV1bHntNZrq6+k5dCg9hgzBNyjoqtahInZT\nIdHRjL1/+JpQAAAgAElEQVTxRpytrVzIz+fskSPUV1cbHUtERD7HYbeTfuAA+Wlp2Dw9GT5rFrGJ\niVe1DhWxG/Ly82u7+q7Zbmf7229TX1VlcCoREfkiTqeTT994g6qSEoZMn078yJF4+vhc8etVxG4o\nYdQoZj7wAAD2hgZW/v731JSXG5xKREQu5+imTVQWF+MbGMiASZPoMWjQFb9WReyGInr1YuDkyUbH\nEBGRK1ReVMSRDRsoLyggfvhwegwefMWvVRG7mdjERPqMGoXJbKaiqIh3fvUrHM3NRscSEZEv43Sy\n6913KTx7lsjevQm/iqunVcRuZvB11zHh5puBi/eo7Xn/fVo1kpaIiNsrzs4m6+RJ6qqqGDhlCqPn\nz7+i16mI3Yx/aChBkZE47HaqSkupKi3F6XQaHUtERL6Cw27nzIEDFGdl0XvYMAZNnXpFr1MRu5Hu\n/fsTFR8PXPxmte2tt4wNJCIiVyVt3z4Kz57FNyiI2P79ienb9ytfoyJ2I0Ouv57+EycCcCE/n32r\nVhmcSERErkZpbi55KSlUlZQQ3acPI+bO/crXqIjdSNyAAUQnJNBYW8uF8+c1gIeISAd05sABMo8e\nJbxHDwZNnYqHt/eXLq8idhPBUVH4BAYCcD49naSNGw1OJCIiX0f6oUNkHDmCxWolvHt3uvfv/6XL\nq4jdxHV3302/sWMByDp+nF0rVhicSEREvo766mpKsrMpLywkrHt3rrv77i9dXkXsJgZMnEjkpQu1\nRESkY0s/dIijmzYRFBHBiNmzv3RZFbGbyUtJIePwYaNjiIjIN1B+/jxFZ88C4B0QwMTFiy+7rIrY\nDYxbtIjwHj0AyDh8mOOffmpwIhER+Sbqqqo4d+IEZw4cwCcggEm33XbZZVXEbmDaXXe1HZauKCqi\nNDfX4EQiIvJNtLa0cCE/n6wTJ/D08WHwtGmXXVZF7AYievbE28+Pxro66qurcdjtRkcSEZFvqCw/\nn/SDBzFbLASEhV12ORWxkUwmYhMT2+4xO3vkCNmnThkcSkREroX6qiryUlMpyMj40uVUxAYymUxM\nvuMO/ENDAUjeuZP0AwcMTiUiItdKZXExR79iXAgVsYFMJhPjb7oJ/5AQACqKi6kqLTU4lYiIXCvV\nZWWc3LnzS5dRERsssFs3LDYbDTU1NDc2Gh1HRESuoYaaGvJSUmioqbnsMipio5hMWKxWTCYTAAc/\n+YSsEycMDiUiItdafXU1u99//7LPq4gN4uPvz8InnsDTxwe4WMTZJ08anEpERK612gsX+OSZZy77\nvIrYIN7+/tz4gx/g6esLQLPdTovDYXAqERG51hzNzRRmZl72easLs8jnmC0WgiIijI4hIiLtzOl0\nfun4ENojdgOnd+2isqjI6BgiImIAFbEBfAIC2qY8BPjsnXcoOnfOwEQiImIUFbEBgqOimPXQQ22P\nM5KSqC4rMzCRiIgYRUVsAP+QEIbNnGl0DBERcQMqYoNVX7igq6VFRLowFbGLmcxmLNaLF6s7W1vZ\n8sorVBUXG5xKRESMoiJ2sbDYWIZMnw5cvKT94Jo11FZUGJxKRESMoiJ2saiEBCbdcUfb4wv5+Zp/\nWESkC1MRu5hPQACRvXq1PW5xOHA6nQYmEhERI2lkLRczWyzYPD1paW7mfHq6LtQSEenitEdskGa7\nnTMHDuBobjY6ioiIGEhF7EJ+wcGERkcD0FRfz/7VqzUHsYhIF6dD0y4UlZBAv/HjAbA3NHBkwwaD\nE4mIiNG0R+xCkb1703f0aKNjiIiIG1ERu1BQRASR8fFGxxARETeiIhYRETGQitgANRcukKTzwyIi\ngorYZUKiowmNiQGgvqaG9EOHDE4kIiLuQEXsIgGhoQSEhQFQW17O8U8/NTiRiIi4AxWxiwRHRbXd\nQ9zc1ERZXp7BiURExB2oiF0kpl8/Yvv3NzqGiIi4GRWxi4TFxhIeF2d0DBERcTMqYhdraW6msa7O\n6BgiIuImVMQuVlNRQV5qqtExRETETaiIXcxeX09VcbHRMURExE2oiF3M3tBAVUmJ0TFERMRNqIhd\noPfw4W1jTNeUl5ObkmJwIhERcRcqYhfw9vPDw9sbuFTEyckGJxIREXehInYBLz8/PC8VsaOpibqq\nKoMTiYiIu1ARu0BYbCz+l4a3FBER+TwVsQuExcURGB5udAwREXFDKmIX8A8Jwdvf3+gYIiLihlTE\nLuAbGIiXr6/RMURExA2piF3MaXQAERFxKypiF2qoqaG6rMzoGCIi4kasRgfoShpqa6m5cMHoGNJF\n+fiAry80NUF1tdFpRORfVMQu1FhToyIWQ4wYAbNmXfxvSQmsWgU7dxqdSkRARexSLQ4HDrvd6BjS\nBQ0aBNOmwahR0NoKLS0qYhF3oXPELtTc1KS5iMXlPD0hMBBsNnA4Lv5/XJzRqUTkX7RH7ELNdjtN\nKmJxMZMJSkuhuBhCQuDMGdi+3ehUIvIvKuJ2ZvP0xGyxANDqcNCsQ9PiYo2NcOIEBAfDuXNw9Chs\n3Gh0KhH5FxVxO/MNCsLDywuAxvp6TfgghkhNhZyci+eHGxuNTiMin6ciFuki6uuNTiAiX0QXa4mI\niBhIRdzOAsLC8PDxMTqGiIi4KRVxO7NYrZjM+msWEZEvpoZwIXtDAw0aW1BERD5HRexCZrMZs1XX\nx4mIyP+nInYhq6cnnjpfLCIin6MibmfNTU20OhxGxxARETelIm5ntRUVNGsEBRERuQwVsYiIiIFU\nxCIiIgZSEYuIiBhIRdzOasvLsescsYiIXIaKuJ05mptpbWkBwD8khLDYWIMTiYiIO1ERu5Cntze+\ngYFGxxARETeiIhYRETGQitiFLDYbNk9Po2OIiIgbURG7kKe3N94BAUbHEBERN6IidiGT2YzZYjE6\nhoiIuBEVsQt5+fnhHxpqdAwREXEjKmIXqC4vp76mBm9/fwJUxCIi8jkqYhcoysykorAQk8mE2ay/\nchER+f/UCi7Q3NCAw243OoaIiLghFbELVBQXU1tRYXQMERFxQypiF6gtL6ehpsboGCIi4oZUxC5Q\nV1VFU309AJ5+foTGxBicSERE3IWK2AVKc3KoKi4GICAsjJ5DhxqcSERE3IWK2AXqqqporKsDwGyx\n4KFhLkVE5BIVsYv5h4QQN2iQ0TFERMRNqIhdzNPbm8Bu3YyOISIibkJF7CL11dXUVVXh4eOjIhYR\nkTYqYhe5cP48ZXl5eHp7E6QiFhGRS1TELlKSk0NJdjZmi0VzEouISBsVsYtUlZRQUVQEgMVqxUfz\nEouICCpil6mrrKSushIAD29vouLjDU4kIiLuQEXsIvVVVW1F7O3nR++RIw1OJCIi7kBF7CINtbU0\n1NYC4OXnR8/Bgw1OJCIi7kBF7EKtLS20OBzYPD013rSIiAAqYpcqyMggbd8+vP38iNehaRERQUXs\nUhUFBeSfOQMmE2az/upFRERF7FL11dVUlZQAFyd/CIqIUCGLiHRxagEXaqitpaasDACrzUavYcOw\n2GwGpxIRESOpiF2osaaG6vJyAGyeniSOG6ciFhHp4lTELtRst1NZXExJTg5WT0/6jR+PVUUsItKl\nqYhdrKqkhMykJKw2Gz2HDNEesYhIF6cidrHG2lou5OdfvFirWzd8AgIwWyxGxxIREYOoiF2spryc\nnOTkiw9MJnoNG4anj4+xoURExDAqYheruXChrYhNQEh0NFYPD2NDiYiIYVTERjKZGDN/Pn7BwUYn\nERERg6iIDVBbUcGpHTswmUz0GzcOLz8/oyOJiIhBVMQGqK+uJnX/fgC8/f11sZaISBemIjZAQ00N\n55KS2h7HDRiAb1CQgYlERMQoKmID2BsaKMjIaHvcb/x4Art1MzCRiIgYRUVsgBaHg+oLFyjOzqa1\npYWeQ4bogi0RkS5KRWyQpvp69n3wAfbGRmL79cM3MNDoSCIiYgAVsUHs9fXs//hjmhsbCQgLwzco\nSMNdioh0QSpigzTb7aQfPIjDbgeg/8SJRMbHG5xKRERcTUVsJKeTjKQkGmtr6T9xIlEqYhGRLkdF\nbCAnkHnkCI11dcQNGEBwZKTRkURExMVUxAZyOp0c3bSJuspKPLy96TFoEJG9exsdS0REXEhFbCSn\nk5zTp2lqaAAgJjFRRSwi0sWoiA3WWFvL+TNnqK+uJrJXL8Li4oyOJCIiLqQidgPHt26lvKCAkJgY\ngiIijI4jIiIupCJ2A6d37aKyuBhvPz9Co6MJCAszOpKIiLiIitgNFGRkUFFYiMNuJ2H0aIZef73R\nkURExEVUxG4iZe9eCjMziUlMJH7kSKPjiIiIi6iI3UTqnj0UZWbiExCAf0iI0XFERMRFVMRuIi8t\njYqiIkwmE93792f0/PlGRxIRERdQEbuJ5sZGck6dovDsWSLj4xkyfbrRkURExAVUxG7kbFISOcnJ\n+IeEENGrl9FxRETkGvANCmLGffdd9nkVsRvJT0ujKDMTi81GaEwMCaNHGx1JRES+IZ/AQK771rcu\n+7yK2I3UXLhAYWYmlcXFhMXGMuHWW42OJCIi35DN05Megwdf9nkVsZvJPnmS1L178QsNZcCkSYRE\nRWEymYyOJSIiX4OXnx/d+/cnIDT0ssuoiN1M9okTnN61C7PZTETPnoycNw+zxWJ0LBER+Rqi4uNZ\n8PjjX7qMitjNNNTWUpqXR3lBAZ6+vvQbNw6zWT8mEZGOKCA8/CvvgtEnvBvKTU5m/+rV+AYEMG7R\nIsxWq9GRRESknaiI3VBlSQm5p0+DyYSHlxdjbrwRLz8/o2OJiMhViB8xglkPPPCVy6mI3VBDTQ25\np0+TuncvNi8v5n73u/gGBRkdS0RErkJUnz4MnzXrK5dTEbuh1pYWSnNzObp5M2aLhT4jRxISGYnV\nZjM6moiIXAFPb28Cw8Pxu4K5A1TEbqrmwgVS9+yhvqoKLx8fhs2ahb/mKRYR6RAGTpnC+JtvBi5+\nnn/4pz9ddlkVsZuyNzZScPYsJ3fupLW1lUm33UZIVJTRsURE5ArEjxzZdrV0zYULrP6///eyy6qI\n3VhteTnb3nyTFoeDXkOHEjdoED4BAUbHEhGRLxEUEUHApSOY9sZGyvLzqSkvv+zyKmI31lhXx+nP\nPqMsP58Wh4NxixYR3bev0bFERORLDJ0+nQGTJgFQXlDAvg8//NLlVcRuzt7UxIa//526ykr6T5xI\nZO/eWHRfsYiI2+o/cWLbpD3F586x4cUXv3R5FbGba7Hb2b96NY21tQR168bg664jpl8/o2OJiMg1\noiJ2c62trZTl53No7VqqSkvpN24csYmJRscSEZEv8K0//pHxt9wCQMaRI2x8+eWvfI2KuANobWlh\n+/LlXMjPJ6ZfP/qOGUNobKzRsURE5H+Jio8nMDwcgOLMTE58+ulXvkZF3EHkpqSQsmcPDdXVDJg0\niX7jxhkdSUREPmfCrbcS3bcvZouF0txczh0/Tl1l5Ve+TkXcQdgbGjixbRulubl0HziQhFGjNP60\niIgbmbR4MZG9ewOQtm8fh9etu6LXqYg7kIzDh8lJTsZsNtNr6FB6Dx1qdCQREQFi+val+4ABePv7\nA1B07hw5yclX9FoVcQdSXlBA8mefUZCRQY/Bgxkxdy4eXl5GxxIR6dpMJq6/916CunUDIPPoUXJO\nnbril6uIO5jknTvJOHwY/9BQBkycqFuZREQMZgKm33MPgZeKeOsbb7Drvfeu+PUq4g6mJCeHtH37\nKMjIIKpPHybdfjtmi8XoWCIiXZLZYqFbz55tn8O1lZU01ddf3TraI5i0r9O7d3Nqxw5CY2IYOWcO\nodHRRkcSEemSgqOiePXcOYIiIgD46M9/5tjmzVe1DhVxB1Sam8uZAwfIT08nOCqKWQ89ZHQkEREB\nzqenU1lUdFWvURF3QK0tLWQmJXF040YCwsIYt2gRwZGRmM36cYqIuEqPQYP49YYNbY83v/IKucnJ\ntLa2XtV69MndQZXm5nL800/JOn6csO7due5b38LD29voWCIiXYanry89Bw9ue3xo7VpKcnKuej0q\n4g7K3thI1okT7P/oI7x8fZn+rW/RfcAAbJ6eRkcTEen0eg0bxs0/+lHb472rVpF7+jT2hoarXpeK\nuAOrKi4m+bPPKMzMpHv//ky45Rb8goONjiUi0ulF9u7NhEuTOzidTra89hqVxcVfa10q4g7M0dxM\n0blz7Fm5EpPJxPibb6bnkCF4+vgYHU1EpNOKTUxk0NSpADjsdlL37iUzKemqb1v6FxVxB1dTVsb+\n1aspzs5u+4YW1r270bFERDqtoTNmsODxxwGor67m/d/+loaamq+9PhVxB+dobqY4K4v1f/87TfX1\njF20iNjERKweHkZHExHpdALCwtqGsnQ0N1Oam8uxLVtobmr62utUEXcCTfX17F+9mobqagLDwxk6\nYwaR8fFGxxIR6XQW/fCH3P7LXwIX7155/3e/+8brVBF3Ai0OB6V5efzzF7+gLD+fcYsW0WvIEKNj\niYh0arUXLnDkCqc6/DIq4k6i1eFg/0cfkbJnDx5eXoxesIC+Y8YYHUtEpNO4+w9/YNaDDwJwNimJ\n955+Gkdz8zder4q4E6mvqmLTSy9RdO4cw2fNYtjMmfgEBBgdS0Skw4vo1YsegwYREB4OQEVhIad3\n774m61YRdzJnk5JI2rCBpro6hs+ezZDp042OJCLS4c168EHiR4wA4OT27axZtoz66uprsm4VcSfT\nVF/P0U2byE9Lo+fgwYyeP18XbomIfANDZ8xg9Pz5hMbEAJCXksKJbduu2fpVxJ1QbkoKp3fvpray\nksTx4xk+a5bRkUREOqyZ999PeI8eAGQcPkzqvn3XdP0q4k6ooaaG/atXc2jtWnyDghg5dy7hcXFG\nxxIR6VAsVivDZ81iwOTJ+AYGUpaXx6aXX2bXihXXdDsq4k4qPy2Ngx9/TNaJE/QeNowZ992Hb1AQ\nmExGRxMR6RA8fXz46apVhMXGArDvww85vWvXNd+OirgTO3vkCMe2bCGwWzdmP/ww42++GYvFYnQs\nERG3Z7ZY8AsJwXRp56WhtpYzBw9SkJFx7bd1zdcobqO+poaU3bvZv3o1AWFhzHvkESw2m9GxRETc\nXnSfPryYmoqXnx8Ar3zveyRt2NAu21IRd2ZOJznJyWx86SWSP/uM6IQE5n33u/gGBhqdTETEbUX1\n6cOiJ5/Ew8ur7c/KCwtpqK1tl+2piDu55qYmcpKT+ezdd/H282PB44/Tb/z4tm95IiLy70Kiohi3\naFHb42fvuYeMw4dxtra2y/ZUxF1AfXU1afv3c2T9ekJjY5n3yCOaKlFE5AsMmDyZO379awLCwnDY\n7Xz65psc37KF2oqKdtumtd3WLG6jtaWFsrw81r7wAv5hYRfvLZ45k7qKCiqKioyOJyLiFiJ69mTM\n/PkMvf56ABx2O5tefrldSxi0R9xlNNXXc3rXLj594w2cTidTli5lyPTp+IeEGB1NRMQtJIwdy4g5\ncwCorahg40svkXX8+Deaa/hKqIi7EIfdzoGPPyZ55066JyYy+Y47iBs0yOhYIiKG6zV0KGMWLKDn\nkCE01NRwcvt23vzxj9u9hEGHprsUp9NJdVkZK556Cp/AQPqMHMngadPIT0ujqqTE6HgiIobwDwnh\nxieeYNqddwJQlJnJuhdecNn2tUfcBeWePs3W11+nvKCA6++9l+n33IOHt7fRsUREDDHmxhvpPWwY\ncPE0Xv6ZMyR/9pnLtq8i7qL2fvABqfv24RcSwoSbbmL0DTcYHUlExBBjbryRXpeKeO+qVfxlyRKX\nbl9F3EU5nU7WPv88h9eupfeIEUxcvJjASxNei4h0Fb9Ys6btAq0T27axa8UKnE6nSzOoiLuw0pwc\nDq9fz7ljxxg4eTI3//jHRkcSEXGZBd//PvHDh+N56dRcTnIyqXv3ujyHirgLa3E4OLVjB5tffZXC\ns2cZNnMmo3SIWkS6iClLlhBw6Ujg1tdfZ+trr9FQU+PyHCriLq6yuJikjRvZu2oVNk9P5n/vewyf\nNQurh4fR0URE2oVPYCAPPPMMPQYOxObpyZ6VK9n8yivkJCcbkke3LwkVhYUc3byZmH79mHH//QBc\nOH+egowMHHa7welERK6d4MhI5j3yCAufeAJMJrJPnmT78uVknThhWCbtEQsAJdnZfPbuu5w7epTh\ns2Yx+fbbCYmOxqz5i0Wkk/AJCGDwtGnc/stfwqV5hnetWME5F4ye9WVUxAJcnKUp6/hxVjz1FGX5\n+cx88EGGzZiBn4bAFJFOIqpPH6YsXdr2OOPwYQ5+8gnlBQUGplIRy+c01tWRtn8/a5Ytw9nSwrxH\nH6XPiBHYPD2NjiYi8o34BASQMHo0YxYsaJsI53c33kheaqrR0VTE8u8aampY89xzfPzMMwSGh3P9\nvffSa+hQTJcO44iIdDRmi4Xr772XR/7xD3A6qSwu5qeTJ1NVWmp0NEBFLF/A2drKppdf5rN33iFh\nzBgW/uAHJE6YYHQsEZGv5Zaf/IS7fvc7ACqKitj40kuU5efT2tJicLKLdNW0fKGm+np2/s//EJOY\nyMCpU2moqaGxttbQKwtFRK7WjPvuY+LixXj7+wMXi3jzK6+4TQmD9ojlS5zPyGDtsmXkp6QwePp0\nRs6bh29goNGxRESuyJDrrmPirbcS268fAGcOHODV73+fSjebbU5FLJfV3NhI2oEDbH3jDcry8phw\nyy1Mu/tulbGIuD3fwECmLFlCwtixeHh7k/zZZ6x46ilO795tdLT/oCKWL2VvaODYli3sef99Wlta\nmHn//Uy6/XajY4mIfKnZDz/M8NmzCQgN5dSOHXz0179ydNMmo2N9IRWxfKXqsjIOr1/PvlWr8PD2\nZs53vsOQ6dMxm/XrIyLuxWKzMWXJEhb+4AeEx8WRc+oUm15+mcNr1xod7bL0SSpXpCwvj/0ffcTe\nDz7Ax9+fO59+mh6DB2tMahFxG54+PoyYPZvHXnuN4KgoKgoL2fL665zcscPoaF9KRSxXrPDsWba8\n9hq7VqwgMj6epU89RVSfPlhsNqOjiUgX5+HtTf+JE/nl2rV4+vjQ4nCw54MPOLZpE1VudnHW/6bb\nl+SqlObmsu6FFwC46cknKS8o4JNnn6UwMxNna6vB6USkK7JYrSSOH88Tb73V9mfF2dnsWrGC/DNn\njAt2hVTEctWqy8pY89xzhHXvzsTFi8FkYt0LL5CXkmJ0NBHpgkbMncuDy5YRHBnZ9md/WbKErOPH\nDUx15VTEctWcTie1FRX88//8H1qamxk2YwYAm19+mXMd5BdfRDqHQVOnMvP+++kWF4fJbKaxro5f\nTJ9OzqlTtDgcRse7Iipi+Vqcra2UFxSw+i9/oaWlhcTx48HpZOsbb3D2yBGj44lIF5AwejTX3X03\ng6ZOxWK1UpaXx3P33Ufm0aMdpoRBRSzfUEF6Ohtfeon6qir6jRvHzPvvB1AZi0i7Shg9mpkPPMDI\nuXPxCw4mLzWVd375S05s22Z0tKumIpZvLPvkSRzNzVhsNgZMnIjT6cTZ2krm0aNGRxORTqjPyJHM\nfOABxsyfT0h0NFnHj/PxM8+w78MPjY72taiI5ZrIT01l14oVeHh50XfsWAAcdjs5yckGJxORziSo\nWzem33MPYxcuJDgyknPHj7Pu+efZ8c9/Gh3ta1MRyzVz9sgRWltaMJnNDJo6FYA1zz5LUVaWW810\nIiIdj8lsxi8oiEm33864RYvwCwkh68QJ1ixbxvbly42O942oiOWaOnfsGOv//ncsVitjbrwRnE4+\nfuYZSnNzO9TFEyLiPixWK4HdujFizhyWPvUUfsHBZB49yru/+Y1bD115pVTEcs3lJiez6o9/xGK1\nMv2ee2htbWXNsmWU5ORoz1hErorZYiEkOpppd9/N3b/9LZhM1FdV8c6vfsWR9euNjndNqIilXRRn\nZfHWj39MeUEBN3zve5jMZtY+9xwFGRlGRxORDiQqPp4F3/8+sx9+GEwmAJ6eP5+0/fsNTnbtqIil\nXTidTmouXOCTZ5/FZDIxcfFi/ENCWPPcc6QfPGh0PBHpAHoOGcKcb3+baXfdhcV6sa5+MmkSZ5OS\nOtXRNZPT6XS6fKOXvtVIF2AyERYby7xHHmHUDTdQkpPDjuXLObZlC/XV1UanExE3NWjKFGY99BDD\nZ80isFs3qsvK+NPixaTt309zU5PR8b6Wy9Wt9oilfTmdlOXl8ekbbwAwcu5c5j/+OP6hoexasUJl\nLCL/YcyCBcy47z4GTJ5MQFgYRefO8eaTT3Jq506jo7ULFbG4REFGBjv++U+a6usZu3Ah1997L62t\nrex5/32VsYi0+X/t3VtM1GcexvEvjOPIQQEBgxaxW8ED1LNoabFQO0WBGrpdjSLaak/B0BhjojZp\nbZvWxLS9aWPVZK2tGENS6NrY9QAVj2EWzwIqbj2N4qqgIrg4DA6W2YtuvNjE7Al5meH5JL+buXqu\n/k/ed97/+58yZw4vvfEGwyZPJrhfP+rOnGHr559T+eOPpqM9Nipi6TJXa2vZ73bTdu8ek3NySM/L\nA6+X0wcOcOPChUdu24hIz5Cel8eMJUt4ctQorDYbl0+dYvuaNezdvNl0tMdK/xFLlwuLjmaM3c7z\nc+cyYMgQTpSWsm/zZq6dO8cDj8d0PBHpYpZevZiQmcmideuIjI2l3ePhb7W17Fi3jp83bDAdr9M8\nqm5VxGKELSiIkampzHr/fZ5ISOBEaSk71q7lam0tnrY20/FEpAsEBAYSFBrK0PHjWbl9O31CQnjQ\n3o6zupo/ffYZf/nhB9MRO5WKWLqd3n36EDtiBNPeeYfnc3OpraigZPVqnFVVeNxubVWL+DGL1Uq/\nyEiSpkxheXExAG0uFw1OJ4UrVnBs507DCTvfI59pXgMAjebh2IKDvbM/+MBb1Njo3XT9ute+cKE3\nuF8/47k0Gs3jm9+NHu1dsmmTd1tHh/cnr9f7k9frzf/6a2/M0KHGsz2ueRStiMW4gIAAgvr2Zcbi\nxbMwkPEAAAWOSURBVGS8/TbtHg+VW7dysKgIZ3W16Xgi0skm5+Tw+2XLGDp2LLaQEDxuN+vy86kq\nL+fuzZt+ey/9o+pWRSzdRkRMDE+npTH19dd5Yvhwzh0+zM/ffEPN3r2mo4lIJ8lZupTM/HyiYmPp\nHRRE47VrrF+0iFqHg9bmZjo6OkxHfGxUxOITQsLDSUhOJm3uXBKSk6m/eJEDRUVU7d5Ny507puOJ\nyP8oNDycrIICXlywgIHx8QBcOXWK75Yv58zBg9x3u8HPz4WoiMVn9A4KIi4pieTsbEZPnUqgxcLx\n0lJ2rl2Lq7lZh7hEfMyo9HQmZmfz3KxZDBgyBE9bG78cOsT2NWuo3LrVdLwu86hnly70kG7H43Zz\n4dgx7t66xd9v32bKnDlkvPUWvW02qvfs4fzRo7S5XKZjisi/YQsOJi0vj+TsbEalpxMcFkaD00lt\nRQWOkhKO+MG3hDuDVsTSrfWLjmas3Y59wQKGp6Rwav9+9m/ZwlmHg+aGBr891CHiywItFiJiYkh+\n+WXyPv2UvpGRBAYGcv7oUQ5v28bhbdu4cvq06ZhdTlvT4rN6Wa0MGjaM11avJi4xkbbWVvYWFnKy\nrIwGp1OrY5FuJCQ8nEHx8Yyx23lt9eqHvzurqij66COqysu539pqMKE5KmLxaYEWC/0HDSI5K4us\nd98lJCyMqvJyDhYVce7IEdpcLr/6PqmILwqNiGBCVhbZBQWMSEnB6/XS5nJxr7GRVa+8grOqynRE\no1TE4jcmZmeTVVDAyJQUbl6+zJ7CQiqKi2mqr8frx68+iHRXAQEBBFoszFu1ij+sWAH8VjrulhYc\nJSV8t2wZ95qaDKc0T0UsfsNqsxE1eDDPzZxJ6uzZRMTEcNbhoGzDBk6WlZmOJ9LjDE5MZM6HHzJp\nxgxswcEA3LhwgT8uXsxZhwN3S4vedkCnpsWPtN+/z83Ll9m9cSO3r17lpTffZFxGBtFxcUQOHMjJ\n3btpqq/XVrVIF8j9+GMmTJ9O7MiR2IKDaXO5OFlWxpaVK7l15YrOcPwHtCIWn9Y3MpK4pCQmZmfz\nwvz5uJqacNbUUPH995zavx9Xc7PpiCJ+x2qzEZeUxMz33mP4M88QHhNDL6uVsw4HjpISavbt43JN\njemY3Y5WxOKXWhob+aWykjvXrnG7ro5JM2YwITOT6Lg4BiUkcLy0lBvnz+vTiiKdJGrwYF6YP5+R\nzz7LuIwMLFYrAHsLCznwz/vhmxsaDKf0LVoRi9/oExpKYmoqEzIzSXn1VdwtLZw5eJDTBw5w7vBh\n6i9dMh1RxGeFRUcTn5zMxKwsJufkEBUbC8CF48cf3gt/6eRJwym7Nx3Wkh4hICCA0P79yczPZ4zd\nTsxTT9Hc0EB1eTlHd+zg4okTeNxu0zFFfEbvPn2InziRp9PSGGO3Myo9HYBfHzygoriYE6WlnCgt\n5e6tW2aD+gAVsfQ446dN4/ncXBImTSI0PBxnTQ271q/n6tmz3K6rw3P/vt9fMi/y/4gdMYJBw4aR\nXVDAuIwMANpcLuovXuTWlSt8kZurw1j/BRWx9EihERGMtdtJmzePxNRUfm1v5/iuXezeuJHr58/T\n0tioazJF/oXFamXAkCHkffIJ46dPJzQigo6ODlqbm/lrZSXbvvyS6vJy0zF9jopYerTBiYm8uHAh\n6fPmERYVxc26OhzFxfz5q69oqq83HU+k2wi0WIiKjeWLQ4eIiIl5+Lvr7l3Kv/2WjUuXGkzn27pV\nEYuIiMhvAk0HEBER6clUxCIiIgapiEVERAxSEYuIiBikIhYRETFIRSwiImKQilhERMQgFbGIiIhB\nKmIRERGDVMQiIiIGqYhFREQMUhGLiIgYpCIWERExSEUsIiJikIpYRETEIBWxiIiIQSpiERERg1TE\nIiIiBv0D84mKoCS3S64AAAAASUVORK5CYII=\n",
       "text": [
        "<matplotlib.figure.Figure at 0xd39be8c>"
       ]
      }
     ],
     "prompt_number": 9
    },
    {
     "cell_type": "markdown",
     "metadata": {},
     "source": [
      "The Earth by contrast is tiny.\n",
      "To even see it,\n",
      "the camera must zoom in very close to the Earth\u2019s orbit \u2014\n",
      "close enough that you can see the gap\n",
      "between where the Earth was sitting at the end of 2013\n",
      "and where it was sitting at the end of 2014.\n",
      "The gap between what ought to be \u201cthe same position\u201d each year\n",
      "is due to the different positions of the other planets,\n",
      "whose gravity tugs and distorts the Earth\u2019s orbit in a different direction\n",
      "every year:"
     ]
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "# Add the Earth\n",
      "\n",
      "x, y, z = earth(JulianDate(tt=jd.tt[-1])).position.AU\n",
      "r = wgs72.radiusearthkm / AU_KM\n",
      "ocean_blue = (0.4, 0.5, 1.0)\n",
      "sphere = mlab.points3d(x, y, z, name='Globe',\n",
      "    scale_mode='none', scale_factor=r * 2.0,\n",
      "    color=ocean_blue, resolution=50)\n",
      "#sphere.actor.property.specular = 0.20\n",
      "#sphere.actor.property.specular_power = 10\n",
      "\n",
      "reposition_camera(distance=0.025, focalpoint=(x, y, z))"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "metadata": {},
       "output_type": "display_data",
       "png": "iVBORw0KGgoAAAANSUhEUgAAAeIAAAGFCAYAAAAhGPBBAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3flz3Hdh//HXHlppdayu1bm6b1mWLdmWndhxQhw7h3N9\niQmhMKSQlClDp0y/Q5kWhk6PDNMCQ6cHhQHaJoEEEhJISUMgwXZiO3Ziy7cly7pvyVrdt7S6vj+E\nLxNCQnxI+97j+fgH8hr/8sx79fm8P5aVlZUVAQAAI6ymBwAAEM4IMQAABhFiAAAMIsQAABhEiAEA\nMIgQAwBgECEGAMAgQgwAgEGEGAAAg+wm/qMWi8XEfxYAwprD6dRf/uhHqtqzR88+9phe++EPNdLX\nZ3pW2Hi/iyw5EQNAmLBarSqorpYjKkp9TU2aHBoyPQkydCIGAPhXXFKS9n7uc0pIS9OBJ55QZ12d\nFnw+07MgTsQAEBZiEhJ0+2c+owiHQ4efeUZD3d2mJ+E3CDEAhDirzaaY+HilZGdrbmZGwz098s3N\nmZ6F3yDEABDiYhMTVXLDDVpaWlL94cOam542PQnvQIgBIMSl5uVp7+c+p0WfTz/7+tc1MThoehLe\ngRADQAizWK2KTUhQTkWFFubn1Xb2LD9LBxhCDAAhzOV2K3vdOvnm5tR04oSWFxdNT8K7EGIACGFZ\nZWXa/pGPaGpkRC//x3/wylIAIsQAEMIS09NVWF0t3+ys6g4f5kQcgLjQAwBClDsrS7mVlVpaXFTj\n8eOaGR83PQnvgRMxAIQoT1mZSrZt0/jgoE798pem5+B9EGIACFHpBQXK37BBk8PDunjkiOk5eB+E\nGABCUHphoQqrq2WxWNR0/LgGu7pMT8L7IMQAEIK23XefNt11l1rPnNGvvvc903PwBxBiAAgxVqtV\nsUlJina5NDc5yU1aAY4QA0CIiUtOljsrSxarVYM9PRonxAGNEANAiLntU5/Sup07df7gQR16+mnT\nc/ABCDEAhJi8DRuUkp2twc5OddXXm56DD0CIASCEpOXnKyEtTZPDwxrs6pJvdtb0JHwAQgwAIWTz\nXXcpJSdHzbW1aj192vQcXAFCDAAhZMOuXUrOylLb2bPqvHDB9BxcAUIMACHClZys1Lw8LS0uaqi7\nW5MjI6Yn4QoQYgAIEes/9CHFJCSo7cwZeTs7Tc/BFSLEABAibn34YSWkpenMK6+o4/x503NwhQgx\nAIQAm92uwk2bFOl0qrexUWMDA6Yn4QoRYgAIctEulx755jcVk5CgF//lX9TOaTioEGIACHI2u12e\nkhLZ7HYNdnVpbnLS9CRcBUIMAEHMZrcrMSNDhZs3a2ZiQp11dZoaGzM9C1eBEANAEItPTdVtn/qU\nXG63Dj39tHqbmrTo85mehatAiAEgiMUlJ+umj35UFotFR59/nk8eBiFCDABBKiIyUkmZmUrOytLo\n5cvydnRoYX7e9CxcJUIMAEEqIS1N5du3a2FuTqdeflm+uTnTk3ANCDEABKnUvDxt37dP8zMz+tX3\nvqe5qSnTk3ANCDEABCGLxSJXcrJy1q2Tb25OzbW1PKQVpAgxAAShhPR0eUpLNTs9rZZTp0zPwXUg\nxAAQhPIqK1Vzzz0aHxjQr//zP7WysmJ6Eq4RIQaAIOTOylL+xo2amZzUxSNHJEIctAgxAASZos2b\nVXnrrZocHtaxn/5U0+PjpifhOhBiAAgyntJS5VZWanRg4O3TMIIaIQaAIOKIilJyVpbikpM17vVq\noL3d9CRcJ0IMAEEkq7xcRVu2aGlhQQ3Hjmmwq8v0JFwnQgwAQWTrvfdq3U03qeXUKR39yU9Mz8Eq\nIMQAECRsdrvSCwqUkJqq0f5+9be2mp6EVUCIASBIZJWXKz41VYPd3fxtOIQQYgAIEjd++MPKXrdO\nl44d0/mDB03PwSohxAAQJPKrqpSUkaGBtjb1NjWZnoNVQogBIAj8/78NTwwNaainR77ZWdOTsEoI\nMQAEgZp77lFiRoYuvfmmOuvqTM/BKiLEABAENt15pxLT09VcW6vuixdNz8EqspseAAB4fxarVZnF\nxcosLtbE0JD6mpo0NTpqehZWESdiAAhgNrtdt37yk4pNTNSpX/6Sh7RCECEGgABmtdlUvHWrHE6n\nui9e1LjXa3oSVhkhBoAAZrXZVFJTo0inUz2XLml8cND0JKwyQgwAASopI0P/9wc/UGRMjP7rC19Q\n+7lzpidhDRBiAAhQkTExqtq9Wza7XQ1Hj2pyeNj0JKwBQgwAAcgeESGX262o2FhNDA9renxcS4uL\npmdhDRBiAAhASZmZuunBB7W0uKiDTz7JK0shjBADQABKzMjQ9t+E+PWnntL02JjpSVgjhBgAAowj\nKkpJGRlKzszU2MCABtrbtejzmZ6FNUKIASDAJGVmqnjrVs1NT+vsq69qaWHB9CSsIUIMAAEms6RE\n2+6/X9NjY3rl+9/Xwvy86UlYQ4QYAAKIxWJRfEqKssvL5ZudVcvJk1peWjI9C2uIEANAAEnKzFRG\nUZFmxsfVeuaM6TnwA0IMAAGkcPNmVd1+u4Z6evT6D39oeg78gBADQABJy8tT7vr1mhodVcOxY6bn\nwA8IMQAEiMpbb337NNzdrUM/+hHvDocJQgwAAaK4pkaF1dUaaGvT2VdfNT0HfkKIASAAOKKiFJ+S\nIkd0tCZHR/ncYRghxAAQALIrKuQpLdXk8LDaz57VzMSE6UnwE0IMAAHgpgcfVMnWrWo4elRvvfCC\n6TnwI0IMAIZZbTZlFhcrPiVFI319utzWZnoS/IgQA4BhORUViktO1kB7u7ydnabnwM8IMQAYdtND\nDymjqEh1hw6p4ehR03PgZ4QYAAwr3rJF8Skp6m1qUn9Li+k58DNCDAAGpRcWKj41VeODgxru7ZVv\ndtb0JPgZIQYAg7bdf79cbrcuHjmi3qYm03NgACEGAIO23Xff2yE+elR9hDgs2U0PAIBwZLPblfOb\nSzyGurvV09DA3dJhihMxABgQFROjD//lXyo6Lk6vP/UUP0uHMUIMAP5mscjmcCglN1cWq1Wj/f2a\nn5kxvQqGEGIA8DN7RIRSc3JUduONmhgeVsvp05ocHjY9C4YQYgDwM3dWlj72t38rq82m5776VXk7\nOkxPgkGEGAD8zOlyaf0tt8hisaj+yBEe0gpzhBgA/MgZF6eMoiJFRkeru6FBkyMjWl5aMj0LBhFi\nAPCjZI9Hm++6S77ZWR16+mnNT0+bngTDCDEA+FGSx6PqPXu0MDeno889x9PSIMQA4C+OqCglZWQo\nMSNDY16v+ltbtbS4aHoWDCPEAOAn7uxsFVRXa2ZiQhdef10rKyumJyEAEGIA8JPc9eu1Ze9ejXu9\nOvD444QYkggxAPhNQnq6ssrKNDc9rZZTpyRCDBFiAPCLzXfdpR379qmzrk7/88//bHoOAgghBgA/\niE9NVWJGhuampjTQ1mZ6DgIIIQaANeaMi1NmUZFcKSnydnaq59Il05MQQAgxAKyx9bfcovKbbtJQ\nV5dqX3qJKy3xOwgxAKyxdTt3Kn/jRvVcuqSz+/ebnoMAQ4gBYA3FJCQoJSdH9ogIDff2auzyZdOT\nEGAIMQCsoeKaGrk9HvU1N6urvt70HAQgQgwAa6jmnnvkKStT0/Hjqjt82PQcBCBCDABrxGK1Kru8\nXC63W0M9PfJ2dJiehABEiAFgjeRv2KDo+Hj1t7RoqKfH9BwEKEIMAGvklk98Qik5OTq3f7+aa2tN\nz0GAIsQAsEbKtm9XXHKyui5e1EB7u+k5CFB20wMAINRYbTZV7d4td1aWWk+dUs+lS/LNzpqehQDF\niRgAVpnVZtPW++5TtMuli2+8wUNa+IMIMQCsMqvVquKaGjmcTnU3NGh8cND0JAQwQgwAqygiKkpF\nW7Yoq6xM/c3N6q6v18z4uOlZCGCEGABWUbzbrU994xtyOJ167p/+ST2NjaYnIcARYgBYRXaHQ8ke\njyRp7PJlLczNGV6EQEeIAWCVRERGKjkrS8mZmRrp79dQT498hBgfgBADwCpJLyzUA1/8opYWF/XM\n3/2dRvnSEq4AIQaAVRKXlKSKnTu1sryscwcPam5qyvQkBAFCDACrINrlUlp+vhxOp3oaGzU5PKzl\npSXTsxAECDEArILUvDxt3L1bs5OTOvLss1pcWDA9CUGCEAPAKkjNzdWGXbs0OzWlY88/ryWfz/Qk\nBAlCDADXKTYpSZ6SEsUlJam3sVH9LS1aXl42PQtBghADwHXK37BBlbt2abS/X288+6zpOQgyhBgA\nrlOSx6PMoiLNTEyo9cwZ03MQZAgxAFyn5MxMZRQXa2ZiQu1nz5qegyBDiAHgOux6+GHd8vGP69Kx\nY3rmH/5BKysrpichyBBiALgO2evWyVNWpuG+PjUdP256DoIQIQaAa5SclaWkzEzNTU/L296uuelp\n05MQhAgxAFyjip07lVVWpu76el08etT0HAQpQgwA16h8+3ZlFherp7GRn6VxzQgxAFyDmIQEpeTm\nyma3a7inR2MDA6YnIUgRYgC4BqXbtikxLU3dly6pt6nJ9BwEMUIMANfgxn37lF5YqItHjujSsWOm\n5yCIEWIAuFoWi3LXr1dccrK8HR0a7OoyvQhBjBADwFX62N/8jdLy8rT/8cdVf/iw6TkIcoQYAK5S\nRnGxHE6nhnt7NTkyYnoOgpzd9AAACBYWq1Wu5GTlb9yoRZ9PPQ0NGh8cND0LQY4TMQBcocjoaN3+\nJ3+ilJwcnXz5ZbWfOyff7KzpWQhyhBgArpDD6dRtn/60omJiVPvSS/J2dJiehBBAiAHgCtgiIpSY\nlqb0ggJNj43J29nJ3dJYFYQYAK5AtMulDbt2yWKx6MJrr2l6dNT0JIQIQgwAVyA+JUW7/viPZbFY\n9NpTT2mUKy2xSggxAHwQi0XOuDgVbtqk5ZUVtZ89q7mpKdOrECJ4fQkAPkBsYqJyKyu16POpubZW\nC/PzpichhHAiBoAPkJaXp50PPaT5mRn94lvf0uzkpOlJCCGEGAA+QHxKikq3bdOiz6dzBw7INzdn\nehJCCD9NA8AfkL9xo3Z89KPyzc3p5W9/WxPDw1pZXjY9CyGEEzEA/AEZxcWq/NCHNDkyomM//alW\nVlZMT0KIIcQA8D7sDodcycmKSUzU9NiYRvr6JEKMVUaIAeB9pOXnK6+yUks+nxqPH9cUl3hgDRBi\nAHgflR/6kLbcfbd6m5r06ve/b3oOQhQhBoD3YLVa5c7OVkpOjiZHRtR98aLpSQhRhBgA3kN6YaHc\nWVkaGxhQz6VLpucghBFiAHgPVXv2qHjrVnVcuKATL75oeg5CGCEGgPeQXV6utPx8Dff2quP8edNz\nEMIIMQC8S9n27coqL9dAe7uaTpzgAw9YU4QYAN6lfPt2pebmqvviRbWfPWt6DkIcIQaAd8kuL1d8\naqqGurt1ua3N9ByEOEIMAO9QUFWlnIoKjXu9aj93TuNer+lJCHGEGADeYd+XvqSs8nId+tGPdPqV\nV0zPQRggxADwG1abTbkVFXLGxmqgo0Oj/f2mJyEMEGIAkGSxWJSUkSGb3a7xwUHNT0+bnoQwQYgB\nQG9/aekTjz2mxPR0vfr976vpxAnTkxAmCDEA6O27pStvvVVRsbFqPX367U8eAn5AiAGEPbvDoazy\ncsUmJsrb0aExr1eLPp/pWQgThBhA2HPGxWnnxz4mu8OhEy+9pBEe0oIfEWIAYc8ZF6ft+/bJ7nDo\nzCuvaOzyZdOTEEYIMYCw5oyLU+GmTUrLy1P72bPqbWrS/MyM6VkII4QYQFhzZ2dr96c/rcWFBb38\n7W9rYmjI9CSEGUIMIGxZLBZFRUcr2ePR0uKiBru6tLSwYHoWwgwhBhC2IqOjlVFUpNz16+Xt7FTT\niRP8LA2/I8QAwlbh5s3a+2d/pvGhIT35V38lHxGGAYQYQNhKTE9X4aZN8s3MqO7QIS0tLpqehDBk\nNz0AAExwud1Ky8vTysqKehobNTs5aXoSwhQnYgBhKaeiQut27tTYwIDeeuEF03MQxggxgLDkKS1V\nydatGh8c5LvDMIoQAwg7zrg4ubOy5HS5NNLXp6HubtOTEMYIMYCwk11eLk9pqUb6+tR0/LjpOQhz\nhBhA2CnfsUOlN96ovuZmHf/5z03PQZgjxADCis1uV3pBgZIzMzV2+bJ6Ll0yPQlhjhADCCv/5wtf\n0Oa9e3X85z/Xi//6r6bnAIQYQHhxOJ2KcDi04PPJNzdneg5AiAGEj9jERGUWF8seGamB9naelkZA\n4GYtIEx48ncpNt6jEe85DfadNz3HiBsfeEAFVVVqOXlS5w8e1NzUlOlJACdiIBwkpKaqsOpWbbv3\n4yq94UbTc4zZfNddSsnNVXNtrVpPnTI9B5BEiIGwULZ9u6rv3qKMcqdsUbOm5xgRn5Ki1NxcLfp8\nGuru1tToqOlJgCRCDISFqj17tPXDWzU+Wq+6QwdMzzFiw65diklMVOvp0xpobzc9B/gtQgyEOIfT\nqYyiIkXFxGiop0fDvb2mJxlxyyc+ocS0NJ3dv1+ddXWm5wC/xcNaQIgrrK5WbFKS+pqbNdjVZXqO\n/1ksskdEKH/jRkVERqq/uVljAwOmVwG/xYkYCHG3Pvyw0vLydPbVV9X45pum5/idIzJS2/ftU2R0\ntM4dPMgrSwg4hBgIcQXV1YpNTFR/S4sGwzBCEVFR2vHgg4qMjtb5gwc11NNjehLwOwgxEKJsdrs+\n8dhjSi8o0K//+79Vf+SIlhYWTM/yO3tEhCp27pQjKkrt585pYnDQ9CTgd/A3YiBEWW027di3TzEJ\nCbrw2mvqb201PcnvEtLTde+f/7nikpL08ne+o56GBi2G4f+MILBxIgZCkNVqVWxiohLS0jQ/M6Px\nwUH5ZsPv/eF4t1u7Hn5YFqtVh55+moe0EJAIMRCCImNitPW++2R3OHTuwIGwDJDNbldMQoKSPB7N\nTk5qqKdHC/PzpmcBv4cQAyEoJiFB937+83I4nXr9qafkDcMLLFxut0q2bdOiz6dzBw5ogS8tIUAR\nYiDEWG02xSYmKqeiQr7ZWfU2NWk2DD9ukFlcrD2PPqq56Wk994//qOnxcdOTgPdEiIEQEx0fr+Ka\nGi0vL6vh2LGw/MKQ1WpVnNutrLIyzc/MqO30aS36fKZnAe+JEAMhJtnj0Z5HH9Xy4qJe/va3NTE0\nZHqS3yWkpSmrtFRz09NqOXlSKysrpicB74sQAyEm2uVS4aZNWlleVnNtreanp01P8ru8DRu07f77\nNe716lff/a5WlpdNTwLeFyEGQkj2unXa9fDDmp+Z0fNf+5qmx8bC8jSY5PEob8MGzU5OquHo0bD8\nN0DwIMRACHG53cqtrNTC3Jwa33orLG/SKt66VVW7d2vc69WRZ57R7OSk6UnAH0SIgRARERWlZI9H\nqbm5GhscVOeFC1paWjI9y+8yi4qUvW6dxgcH1XDsmOk5wAcixECI8JSUaMOuXbJYLDrx4osa7u0N\nu7+NOpxOubOz5XK7NTYwoIG2NtOTgA9EiIEQUVhdrU133qmRvj4dfPJJ03OMyK2oUHFNjRbm5lR/\n5AhfWkJQIMRACHBERSnJ41FCWprGBwfV39JiepIRm+++W+U7dqjtzBkde/5503OAK0KIgRCQXlSk\njKIijXu9aj1zxvQcI+wOh9Lz8+VyuzVy+bIu87M0ggQhBkLA+ptv1vqbb1ZvY6PeePZZ03OMyC4v\nV3xqqga7uogwggohBkJAal6eUvPyNDE0pK6LF03PMeLGBx5Qzvr1unTsmM7t3296DnDFCDEQ5LZ/\n5CNad9NNajpxQm8891xYvjssvX2bVlJ6ugY6OtTX3Gx6DnDFCDEQ5DIKC5WYnq7R/v6wfUgro6hI\n8SkpGvN6NdzTwycPEVQIMRDEol0uZZaUKDImRt7OTg12dpqeZMTWe+9VUmamGo4dU8f586bnAFeF\nEANBbOPu3SqoqlJfU5PqjxzR1Oio6UlGVN9+uxLS09Vy8qS6GxpMzwGuCiEGgtiNDzygjKIiNRw9\nqoajR03P8TuL1arirVuVWVIib3u7uuvrNT02ZnoWcFUIMRCkol0uZRQWymqzabCrS+Ner+lJfmez\n23XHZz6juORkvfXzn4ftE+MIboQYCFJlN96omPh4ddXVabC72/Qc/7NYZHc4lJqXJ6vVqtG+Ps1O\nTZleBVw1QgwEqdseeURJHo9O/vKXajt92vQcv7Pb7UrPz1dxTY0mhofVdvasJgYHTc8CrhohBoJU\nYXW1nHFx6m9u1kh/v+k5fudyu/XHX/uaIqOj9bOvf109jY2mJwHXxG56AICr43A6teuTn1RcUpKO\n/8//qLexMew+dyi9/e9QcfPNstntajl5UlMjI6YnAdeEEzEQZOwREdpw221yREWp+eRJjQ4MmJ7k\nd1ExMcoqL1ek06nuixc1NTqq5aUl07OAa0KIgSBji4hQybZtioiKUld9vSaHhkxP8ruE9HTd+MAD\nWl5a0uEf/5hXlhDUCDEQRGKTknTLxz+uZI9Hb73wgrovXpQvDK9zdLnd2nTHHVpeWlLtSy9pZnzc\n9CTgmvE3YiCIJKSl6d7Pf15Wm00vf+c7GurpMT3JiNmJCV184w0t+nzqb23VYph+6AKhgRADQcIe\nEaH4lBSlFxZqenRU/a2t8s3Omp5lRH9Li77/+c9rZWVF8zMzpucA14UQA0Eizu1W6bZtWvT5dO7g\nwbCNsCQtLixoLAxvEkNo4m/EQJBIy8vTzX/0R/LNzuoX3/qWZiYmTE8CsAoIMRAkYhMTlbdxo5YW\nFtR04gTf3AVCBCEGgkB8Soo8paVamJtTy+nTWlpcND0JwCohxEAQyF63TpvvuksTQ0M68MQTWuHy\nCiBkEGIgCLizs1W0ZYtmJyd1/sABLa+smJ4EYJXw1DQQ4Ao3bVL17bdrdnJSR559VuN8YQgIKZyI\ngQBXvmOHqu+4Q5fb2nTwySdNzwGwyggxEMDsDocS0tIUl5SkqZERDXZ1mZ4EYJURYiCApeXlKdnj\n0bjXq96mJtNzAKwBQgwEsBs+/GFtvO02NdfW8rM0EKIIMRDA0gsL5c7O1pjXq57GRtNzAKwBQgwE\nKE9pqRJSUzXU06PLbW0SrywBIYkQAwGq5p57lL1unZpra1X3+uum5wBYI4QYCFBFmzfLnZ2t/pYW\nddXXm54DYI0QYiAAbbrzTmWVl6vj/Hk1HT+u2clJ05MArBFCDASg6j17lJierpZTpzgNAyGOEAMB\nJjI6WmmFhbLZ7Rru7tbk8LDpSQDWECEGAojFalVxTY3yKis12N2ttrNnuVsaCHGEGAggNptNH/nr\nv1Zierpe/+EP1XLqlOlJANYYIQYCiNVmU9GWLXI4neq5dEnjXq/pSQDWGCEGAkREVJRKtm2TLSJC\nHefPa3pszPQkAH5AiIEAEe1y6a7PflaOqCgd/vGP5e3sND0JgB8QYiBARDqd2rhnj+wREWp86y1N\nDA2ZngTADwgxEADiU1O1fd8+xSYk6Ojzz2ugvV2LPp/pWQD8gBADAcCVnKzqO+7QyvKyan/xC02P\nj5ueBMBPCDFgmM1uV1xystLy8zU7NaW+piYtzM2ZngXATwgxYJgrJUWlN9ygZI9H9UeOqPviRS3M\nz5ueBcBPCDFgWO769dr9yCOaGh3Vs489pvmZGdOTAPgRIQYMslgscrndyi4rk292Vi0nT2ppcdH0\nLAB+RIgBgxIzMpRZXKyZyUm1nj5teg4AAwgxYFDhpk3adOedGunr08EnnzQ9B4ABhBgwKCUnR3mV\nlZoeHVXD0aOm5wAwgBADhuRUVKho82ZNj4/rwqFDvDsMhClCDBiSV1mp3MpKDff0qP7QIdNzABhC\niAFDUvPylJqXp4mhIXXW15ueA8AQQgwYkFFUpLwNG7S8uKi2M2c01N1tehIAQwgxYMDuRx7Rhl27\ndO7gQe1//HHTcwAYRIgBP7PabMosKpIrJUUjfX263NZmehIAgwgx4Gc5FRWKS06Wt6ND3s5O03MA\nGEaIAT/b8eCDyiwpUf3hw7p4+LDpOQAMI8SAnxVt3qyE1FT1NTerr6XF9BwAhhFiwI/SCwsVn5qq\n8cFBjfT2yjc7a3oSAMMIMeBH2+6/X/EpKbr4xhvqvnTJ9BwAAYAQA35Uc/fdik9J0aU331RvY6Pp\nOQACgN30ACAc2CIiVFxTo6yyMvU1N6vzwgVNj42ZngUgAHAiBvwgKiZGj3zzm4pLTtb//tu/qe3M\nGdOTAAQIQgysNYtFEZGRKqyuls1uV29jo6ZGR02vAhAgCDGwxuwREUrKyJDFYtHI5ctamJ83PQlA\nACHEwBpLyszUJx57TFa7Xc999asa4EpLAO9AiIE15oyN1fpbbpHValXD0aOa4iEtAO9AiIE1FBUb\nq8ziYkVGR6u7oUFTIyNaXloyPQtAACHEwBpK9ni05Z57tDA3p0NPP63ZqSnTkwAEGEIMrKHEjAxV\n79mjBZ9Px372M81NT5ueBCDAEGJgjSRlZqrshhsUm5ioukOHdLmlRUsLC6ZnAQgwhBhYIwXV1br5\n4x/X6MCAXvjGN7TE34YBvAdCDKwBi9Uqu8Mhm92upYUF+Xh3GMD7IMTAGoiJj1d2WZlSc3PV39qq\n1lOntLK8bHoWgABEiIE1sOmOO3TTQw+pu6FBL3zjG6bnAAhghBhYA2n5+coqLdXk8LCaa2tNzwEQ\nwAgxsMqSMjOVnJWl+dlZ9be2an5mxvQkAAGMEAOrrGz7duVXVamvqUnn9u83PQdAgCPEwCor2rRJ\n2eXlGmhvV8PRo6bnAAhwhBhYRdnr1imnokKzU1NqO3tWYwMDpicBCHCEGFhFW++9VzkVFWo6flzn\nDhwwPQdAECDEwCqx2e1KzMhQZEyMJoeHNTk0ZHoSgCBAiIFVkuzxqKCqSpLUWVengY4Os4MABAVC\nDKySD3/xi8pet05v/OQnOv3KK6bnAAgShBhYJUWbNysuKUl9TU0a5DQM4ArZTQ8Agp3FYlFWebli\nExM12NWlkf5+LfK5QwBXiBMxcJ0sVqtu/tjHFJuUpLP796uvqcn0JABBhBAD18lqtarmnnsUEx+v\nS8eOyctaTS+3AAAHMUlEQVTP0gCuAj9NA9fBERWl8h07lF5QoI7z59XT0KDZqSnTswAEEU7EwHVw\nuly680//VBFRUTryk5/I29VlehKAIEOIgWtksVgU6XQqraBAy8vLGuru5ktLAK4aIQaukcPplKe0\nVIWbNsnb2am2M2c0Mz5uehaAIEOIgWuUlp+vj37lK1peWtLTX/mKRvr6TE8CEIQIMXCNouPjVVxT\no5XlZTUdP6656WnTkwAEIUIMXIOYhAR5SkpktVrVfu6cZqemtLK8bHoWgCBEiIFrkJafry13363Z\nyUkdeOIJLczPm54EIEgRYuAaJHs8qti5U/MzMzrxv/+rRUIM4BoRYuAqpeTkqHTbNkVERqru0CEN\ndXdrmZ+lAVwjQgxcpbzKSm247TaNXr6sQz/+sek5AIIcIQaugs1uV5zbrfiUFM1OTmqQm7QAXCdC\nDFyF+NRU5VVWKs7tVmddnbrq6kxPAhDkCDFwFbbdf7+23X+/2s+c0Uv//u+m5wAIAYQYuAqpublK\nzcvTxPCwuurrTc8BEAIIMXCFUvPylJSZqcmREfW3tGjR5zM9CUAIIMTAFdpw663Kq6xU14ULuvDa\na6bnAAgRhBi4QkVbtii9sFD9bW1qOXnS9BwAIYIQA1egZNs25VRUaLi3Vy21tZocGTE9CUCIIMTA\nFdiyd6/SCwrUdOKEGk+cMD0HQAghxMAHsDscSsnJUURkpMYGBjQ5NGR6EoAQQoiBD+ApLVV+VZVm\nJyfVfu6chnt7TU8CEEIIMfABHvzyl5WWn6/XnnpKdYcOmZ4DIMQQYuADFFZXK9rlUn9zs0b7+03P\nARBi7KYHAIHKZrcrt7JSUbGx6m1q0vjgoJaXlkzPAhBiOBED7yMiKkp7Hn1Uzrg41b74ovpbW01P\nAhCCCDHwPuwOh7bs3auo6Gg1vPmmRvr6TE8CEIIIMfAeol0ubb7zTrmzsnTh0CH1NTbKNztrehaA\nEESIgfcQm5ionQ89JIvNpjd/9jONeb2mJwEIUYQYeBer1Sqny6W0ggItzM9roL2d0zCANUOIgXdx\nulwqqKpSVlmZOs6fV/u5c5qbnjY9C0CIIsTAu2SVluq+v/gLzU1P64df+pKm+MADgDVEiIF3iUlM\nVP7GjVpaWFDTiRPyzc+bngQghBFi4B1cbrc8JSVamJ9X25kzWlxYkFZWTM8CEMIIMfAOWeXl2rx3\nryaHh3XgySe5SQvAmiPEwDukZGerpKZGM5OTOvfrX2t5edn0JAAhjrumgd8oqKpS1Z49mp+Z0bHn\nn+fdYQB+wYkY+I2imhpV3nqrvF1dOvLMM6bnAAgThBiQZI+IUFxyspwul2bGxzU2MGB6EoAwQYgB\nSSm5ufIUF2txfl6ddXWa5N1hAH5CiAFJNz7wgDbu3q3m2lodeOIJ03MAhBFCDEhKLyiQOztb416v\nehoaTM8BEEYIMcKep6RECWlpGunr0+W2Nq1wgQcAPyLECHub9+5V9rp1ajl5Uhdef930HABhhhAj\n7BVv2aKU7Gz1tbSoq67O9BwAYYYQI6xt37dPeRs2qLm2VvWHDmlmYsL0JABhhhAjrG3ft08pOTmq\nP3xYzSdPmp4DIAwRYoQth9OpxIwMyWLR+OCgZjkNAzCAECMsWaxWFW3aJHdWlgY6OtTf2qq56WnT\nswCEIUKMsGS12fSRL39ZiRkZeu0HP1ALP0sDMIQQIyxZLBYV19QoKiZGPZcucbc0AGP4DCLCTkRk\npIq3bpU9IkId589ranTU9CQAYYwTMcJOtMuluz77WTmionTk2Wc12NlpehKAMEaIEXYcTqeq9+yR\n3eHQpTff1PjgoOlJAMIYIUZYScrI0G2f+pRik5K0//HH1dfcrEWfz/QsAGGMECOsJHk82vPoo1pe\nXtavvvtdTXi9picBCHOEGGHDZrcrNiFB7uxszU1NaaC9XQuchgEYRogRNuKSk5VbWamFuTldevNN\nLS0smJ4EAIQY4SO3slJ7Hn1UU6Oj+vHf/73mZmZMTwIAQowwYbHI5XYrp6JCvrk5NZ84wYkYQEAg\nxAgLSRkZyiwu1szEhFpPnzY9BwB+ixAjLORXVWnTHXdotL9fB3/wA9NzAOC3CDHCQmpurvI2bNDU\n2Jga3njD9BwA+C1CjJBXcfPN2nTHHRrp79ehp5/mbmkAAYUQI+St27FD63buVF9Tk449/7zpOQDw\nOwgxQlpMfLzc2dmKiIzUcG+vRvr7TU8CgN9BiBHS8quqlJafr4H2drWdPWt6DgD8HkKMkLZx927l\nV1Wp4/x5nfv1r03PAYDfQ4gR0jIKC5WYnq6xy5fV39pqeg4A/B7LysrKiukRAACEK07EAAAYRIgB\nADCIEAMAYBAhBgDAIEIMAIBBhBgAAIMIMQAABhFiAAAMIsQAABhEiAEAMIgQAwBgECEGAMAgQgwA\ngEGEGAAAgwgxAAAGEWIAAAwixAAAGESIAQAw6P8BJaaEWEAYYz4AAAAASUVORK5CYII=\n",
       "text": [
        "<matplotlib.figure.Figure at 0xd640f6c>"
       ]
      }
     ],
     "prompt_number": 10
    },
    {
     "cell_type": "markdown",
     "metadata": {},
     "source": [
      "The tiny dark-blue dot in the center,\n",
      "at the top of the orbit segment that enters from the bottom,\n",
      "is our home.\n",
      "\n",
      "If you view this live in a Mayavi window,\n",
      "be sure to first press `+` (or the same key un-shifted, `=`)\n",
      "several times to zoom in on the Earth\n",
      "and get a feeling for its size.\n",
      "Then press `-` until you have zoomed back out\n",
      "to the size of the Solar System,\n",
      "to get a feel for how different the scales are!"
     ]
    }
   ],
   "metadata": {}
  }
 ]
}