Skip to content

Instantly share code, notes, and snippets.

@jakevdp
Created December 4, 2015 14:46
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save jakevdp/43baca1b31a97125a9cc to your computer and use it in GitHub Desktop.
Save jakevdp/43baca1b31a97125a9cc to your computer and use it in GitHub Desktop.
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Data taken from https://en.wikipedia.org/wiki/Kirkwood_gap"
]
},
{
"cell_type": "code",
"execution_count": 1,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"from __future__ import division\n",
"import warnings\n",
"warnings.simplefilter('ignore', FutureWarning)\n",
"\n",
"%matplotlib inline\n",
"import numpy as np\n",
"import pandas as pd\n",
"import matplotlib.pyplot as plt\n",
"\n",
"plt.style.use('fivethirtyeight')"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"data = pd.DataFrame({'Distance (AU)':[2.06, 2.5, 2.82, 2.95, 3.27, 1.9, 2.25, 2.33, 2.71, 3.03, 3.075, 3.47, 3.7],\n",
" 'resonance':[4/1, 3/1, 5/2, 7/3, 2/1, 9/2, 7/2, 10/3, 8/3, 9/4, 11/5, 11/6, 5/3]})\n",
"data = data.sort_values('Distance (AU)').reset_index(drop=True)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Here's the plot from Wikipedia, which is really confusing:"
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAZ4AAAD+CAYAAAAK/2PhAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3XlclWX6+PHPc3ZABAQUSBEV3E3cUpM0rTQtMzO3crKx\nbUxnbNHMflkyLbaYjZX6bbHJxjRzD1NrTDS3ccks01RccElkUzbhwFme3x/ESYSzgOxc79eL18Tz\n3Oc512GUy3u7biUjI0NFCCGEqCKa6g5ACCFE/SKJRwghRJWSxCOEEKJKSeIRQghRpSTxCCGEqFKS\neIQQQlQpSTxCCCGqlNvEM3v2bAICAop9tW3b1uVrDh8+zJAhQwgNDaV9+/a89dZbFRawEEKI2k3n\nSaPWrVuzfv16x/dardZp26ysLIYPH05MTAzx8fEcO3aMyZMn4+3tzeTJk68/YiGEELWaR4lHq9US\nHBzs0QNXrFiB2Wxm4cKFGI1G2rZtS0JCAgsWLJDEI4QQwrM5nsTERNq1a0fnzp155JFHSExMdNp2\n79699O7dG6PR6Lg2YMAAkpKSOHv27HUHLIQQonZzm3h69OjBwoULWbVqFe+99x7JyckMGjSIy5cv\nl9o+JSWFxo0bF7tW1FtKSUmpgJCFEELUZm6H2m6//fZi3/fo0YPOnTuzdOlSJk2aVKK9oigVF50Q\nQog6x6M5nqt5e3vTtm1bTp8+Xer9xo0bl+jZpKamOu45k5CQUNZQhBBC1FBRUVFO75U58ZjNZo4f\nP07fvn1LvX/TTTcxa9Ys8vPzHfM88fHxhIWFER4eXq4ghRBC1B1u53hefPFFdu7cSWJiIvv372f8\n+PHk5eUxduxYAGJjYxk2bJij/f3334+XlxdPPvkkv/32G19//TXz5s3jySefrLxPIYQQotZw2+NJ\nSkri0UcfJT09naCgIHr06MHmzZtp2rQpAMnJycVWuTVs2JA1a9YwdepU+vfvT0BAAJMnTy51PkgI\nIUT9o8gJpEIIIaqS1GoTQghRpSTxCCGEqFKSeIQQQlQpSTxCCCGqlCQeIUSNUnQUiwBVVenTpw+v\nv/56lb3nzJkzS1SsqWiSeIQQleaLL74odpZXSEgI7dq1Y8SIEXz44Yfk5OSUeI2iKGUuvZWXl8fs\n2bPZsWNHRYVeI6xcuZLExEQmTpxY6v1t27YREBBAeHg4+fn5pbbp1KkT999/f6n3zpw5Q0BAAP/6\n178c1yZPnsyvv/7Kxo0br/8DOFHmygVCiKqVlGvjYFoBB9MtHEwrACA6yEB0oJ7oIAOh3s7Px6rO\nZ19txowZtGjRAovFQkpKCtu3b2fGjBnMnz+fZcuW0aFDB0fbadOm8cwzz5Tp+bm5ubz11ltoNBpi\nYmIqJOaa4L333mPYsGFOe4BfffUV3t7e5OTksGnTpmKb+Yt4ksivvt+kSRMGDx7Me++9x+DBg6/v\nAzghiUeIGiwp18bYzeklrn97zsy358wALLs9sFwJojKffa3bbruNbt26Ob5/6qmn+OGHHxgzZgxj\nx45l7969mEwmoPD8L1eHTbqiqnVnW+LPP//Mr7/+yksvvVTqfbPZTFxcHI899hjr169n+fLlpSae\n8vxMhg8fzvjx4zl16hQtW7Ys8+vdkaE2IWqwol7I9bap6md7om/fvkybNo1z587x1VdfOa6XNsdz\n8OBBRo4cSWRkJCEhIXTq1IknnngCs9nMmTNniIyMBODNN990DOsVlek6e/YsU6dO5aabbiIsLIzm\nzZszevRojhw5Uuw9tm/fTkBAAKtWreKdd96hffv2hISEMGzYsFKLIv/000+MHj2aFi1aEBYWxs03\n38y7775brM2JEyd4+OGHadmyJSEhIfTt25d169Z59PP55ptv0Gq1Tutibty4kezsbEaPHs3o0aPZ\nvHkzGRkZHj3bnX79+qEoCt98802FPO9a0uMRogY7mG5x2+bNg9m8eTC70t5/cLhXpTwbYPTo0fzz\nn/8kPj6ehx56yHH96qGftLQ0hg8fTnBwMFOmTMHf35/z58+zadMm8vLyCA4OZu7cuTzzzDMMHTqU\noUOHAtCiRQugMEHs2rWLe++9l6ZNm5KUlMRnn33GXXfdxf/+9z+aNGlSLKZ58+ah0+n4xz/+QWZm\nJu+99x6PPfYYmzdvdrTZtm0bo0aNonHjxjzxxBOEhoZy/PhxNm3axNNPPw3AsWPHGDhwIKGhoTz1\n1FP4+PgQFxfHww8/zIcffsioUaNc/mz27t1L27Ztix2qebXly5fTsWNH2rVrh5eXF7Nnz2bNmjX8\n9a9/LcP/A6Xz8/OjRYsW/O9//+Pvf//7dT/vWpJ4hKjBKrPHURPePywsDF9f3xI9iquHh/bs2UNG\nRgZr1qwhOjracX3GjBmO/77nnnt45pln6NChAyNHjiz2rEGDBpUYgho9ejS9evXiP//5D1OnTi12\nr6CggC1btqDTFf569Pf35/nnn+e3336jXbt22O12pkyZQnBwMNu3b8ff37/Uz/b8889zww03EB8f\n70gejzzyCPfddx+xsbFuE09CQgJdunQp9d6lS5fYsmULM2fOBCAiIoKePXvy1VdfVUjiAWjevDnH\njh2rkGddS4bahBDVqkGDBqWubivi5+cHwKZNm7BarWV+ftHcERQuQrh06RK+vr60atWKgwcPlmg/\nZswYR9IB6NWrF1C4AgwK517OnDnD3/72N6dJ5/Lly2zbto1hw4aRk5NDenq642vAgAFcuHCBkydP\nuoz70qVLTp+/evVqbDZbsSQ7ZswY9uzZw9mzZ10+11P+/v6kp5ecA6wIkniEqMGigwx1/v1zcnJo\n0KCB0/sxMTEMGzaMN998k5YtWzJmzBg+//xzcnNzPXq+2WzmpZdeom3bttxwww20atWKyMhIjhw5\nQnZ2ySHKosr7RYp++RfNnxT1ztq3b+/0PU+dOoWqqrzxxhtERkYW+5o5cyaKojgOyHRGURSnCwO+\n+uorbrzxRvLy8jh16hSnTp3ixhtvRKPRsGLFCpfPdfZe11JVFY2mclKEDLUJUYNFB+odK8ycmR7t\nW655mI1n89zODUUH6sv83LL4/fffyc7Odrty6rPPPuPAgQNs2rSJrVu3MmXKFObOncvmzZsJCgpy\n+drnnnuOL774gieeeIKePXvi5+eHoijMmDEDu91eor2zFXVlWR1W9NxJkyZxxx13lNqmbdu2Lp/R\nqFGjUhcLJCYmsm/fPgC6du1a4v6KFSt49tlnHd+bTCby8vJKfY+i66XNI2VkZBAYGOgyxvKSxCNE\nDeZJj6O8vZLKfLanli9fDsCAAQPctu3atStdu3blhRdeYPPmzYwcOZLFixfz7LPPutynsnbtWsaO\nHVti9//ly5fL9Yu1aNHC4cOHncYdEREBFCaxfv36lfk9ANq0aVPsrLMiy5cvR6vV8tFHHxUbEgT4\n5ZdfeOedd/j555/p3LkzAM2aNSuxgq/I8ePHHW2udebMGZe9uushQ21C1GCh3lqW3R7I9GhfBjUz\n0cRLQxMvDYOamZge7Xtd+2wq89me2LZtG2+//TYREREuJ9ozMjJK9DZuvPFGALKysgDw8irs8V2+\nfLnE63U6XYmezcqVK7l48WK54o6OjiYiIoIPP/ywRI+kKM7g4GBuueUWFi9eTFJSUolnpKWluX2f\nnj17cuzYsRIVCVasWEGvXr247777uOeee4p9TZkyBYPBUGx5+sCBA0lOTi6xjNtqtfLpp5/i4+NT\nYtNtZmYmiYmJ9OzZ022c5SE9HiFquFBvLaHhXpWyrLkyn321zZs3c+LECaxWK6mpqfzwww9s3bqV\n8PBwli1bhsHgvGe1dOlSPvnkE4YOHUpERAR5eXl88cUX6HQ6x2o1Ly8v2rVrx+rVq4mMjCQgIICI\niAi6devG4MGD+fLLL/H19aVdu3YcOnSINWvWEBERUa7NlYqi8O677zJq1ChiYmIYN24cISEhnDp1\nir1797Jp0yYA5s6dy6BBg+jTpw/jx48nIiKC1NRU9u/fz/Hjxzlw4IDL9xkyZAizZ89m27ZtDBw4\nEIAff/yRkydPOl255uvrS58+fVi9ejWvvPIKGo2Ghx9+mCVLlvDII4/w/fff07lzZ7Kysvj66685\nePAgs2fPpmHDhsWes3XrVlRVZciQIWX++XhCEo8QotIUDYG98cYbABgMBgICAujQoQNvvPEGDz74\nID4+PiVec/XQWUxMDAcPHmTNmjWkpKTg6+tL586dmTNnTrE5jvfff5/p06fz4osvkp+fzwMPPEC3\nbt1444030Ov1rFmzhiVLltClSxdWr17Niy++WGKIztMacbfeeivffPMNb775JgsWLMBmsxEREcHo\n0aMdbSIjI4mPj+fNN9/kyy+/JD09neDgYDp27MiLL77o9j06duxI586dWbNmjSPxLF++HEVRXJay\nGTJkCFu3buWHH37g1ltvxWQysWHDBt5++22++eYbvvzyS4xGI506deLzzz937Hu62tq1a+nVqxet\nWrXy6OdRVnL0tRBC1FCrVq1iypQpHDp0qMoqdl+8eJHo6Gg+/fTTSuvxyByPEELUUCNGjCAiIoKF\nCxdW2XvOnz+fjh07VlrSAenxCCGEqGLS4xFCCFGlypR45s6dS0BAANOmTXPapuhgoWu/tmzZct3B\nlpWSmoRu+0aMH8/G+9nReD87GuPHs9Ft34iSWnKJoxBCiMrn8aq2ffv2sXjxYjp06ODRyo/Vq1fT\nsWNHx/fOag5VFiU1CZ+pY0tc1+z4Fv2ObwG4MmcZanBolcYlhBD1nUc9nszMTB5//HHmz5/vcQIJ\nCAggODjY8aXXV27pjWtpj5Ys/leeNkIIISqWR4nnqaee4t577yUmJsbjDVfjxo0jKiqKO++80+OD\njyqSJB4hhKiZ3A61LV68mMTERD755BPA/QYrX19fXn31VXr16oVWq2XDhg1MmDCBhQsXuj1/oiJJ\n4hFCiJrJZeJJSEjglVdeYdOmTY6Kraqquuz1NGrUiEmTJjm+j46O5vLly8ybN89l4klISChr7C61\nt1gp/dy+P1ks1gp/XyHK6oYbbsDb27u6wxCiTHJzc/n999+d3o+KinJ6z2Xi2bt3L+np6Y6DkABs\nNhu7d+/ms88+48KFCx7N3XTp0oUlS5a4bOMqyPLQdOoOfywicNWmot9XCCHqA29v73L//nSZeO6+\n+266devm+F5VVSZNmkRkZCTPPPOMxwsGDh06REhISLkCLC9b22jH6jWntFW74EEIIYSbxOPn5+c4\ndraIl5cXfn5+jkOMYmNjOXDggGMBwdKlSzEYDHTq1AmNRsOmTZtYtGgRsbGxlfQRSmdrG+22jW7H\nJix9B2OP7FAFEQkhhIByVKe+tnJscnJyscOKFEVhzpw5nDt3Dq1WS2RkJPPnzy92NnhVUINDuTJn\nGdqjBx1fWCwoWZdR/pijUmxWTPNeJG/W/6EGNqnS+IQQor6qd7XadLs3Y/q/V4tds4W3Iu//vQ8m\nmeAVQojKVu8SD4Bh5ScY4oovdrB2uwXz5FjQSPk6UbMoqUnFe+4UDiUXfV1P9Y3KfHZprly5UuL8\nHVH/1MvfsgX3TcDa7ZZi13Q/bsew+tNqikiI0hWVfjJ98ib6Hd+iSUtGk5aMfse3mD55E5+pY8td\nd7Aynw0we/ZsAgIC+O2333j88ceJiIjg5ptvBmDLli0MGTKEpk2b0rRpU+6//34OHTpU7PUpKSn8\n/e9/p0OHDjRp0oTWrVszcuRIjh49Wqzdp59+Su/evQkJCaFNmzY8/fTTJY6kvuuuu7jppps4evQo\nQ4cOJSwsjPbt2/Pee+8Va2exWHj99dfp378/ERERhIaGctttt/HNN9+U+HwBAQE888wzrF+/nt69\ne9OkSRN69+7N999/X6LtxYsXeeqpp2jfvj1NmjThxhtvZMqUKeTk5DjaZGZm8sILL9CxY0caN25M\ndHQ0b7/9dolju+uC+nkCqUaD+fEZeL12Ae3Zk47Lhrgl2G+IwNr79moMTog/eboR2lqOnkllPvtq\nEyZMoHnz5rz00ksUFBSwYsUKHn/8cW677TZefvllzGYzixcvZsiQIWzZssWxRHf8+PEcOXKExx9/\nnObNm5OWlsauXbs4efKkY3HT22+/zeuvv06/fv2YMGECp06d4pNPPmH//v18//33jiO1FUUhKyuL\nkSNHMnToUO677z7Wrl3Lyy+/TPv27bn99sK/81lZWSxevJgRI0bwl7/8hby8PFauXMm4ceNYsWKF\no12RoqOuH3nkEXx8fPjwww956KGH+PXXXx0HtyUnJ3Pbbbdx+fJlxo8fT7t27bhw4QLffPMNly9f\npkGDBuTl5TF06FDOnz/PhAkTCA8PZ//+/bzxxhucO3euRIKs7erlUFsRJT0Zr1l/Q5N12XFN1evJ\nm/Ee9lbtqjEyIQoZP57tfltAJbLEDCL/sRnleu3s2bN56623uOeee1i8eDFQONTWoUMHhg4dyvvv\nv+9om5GRQY8ePbj11lv5+OOPycjIoEWLFrzyyitMnjy51OenpaXRvn17+vTpw+rVqx2LnpYuXcqk\nSZN46623eOyxx4DCHs+uXbv4v//7P8fx1BaLhU6dOtGzZ09HfHa7HZvNVmyriMVioV+/fjRu3Ji1\na9c6rgcEBGAwGNizZw8REREAHD58mJiYmGLvPXHiRFasWMF3331X7Kjuq73zzju88847/PDDD0RG\nRha7/uqrr7Jv375i12u7ejnUVkQNbIJ5yquouj//kCkWC6Z5/w/lUko1RiZEoeou61QR7z9hwgTH\nf8fHx5OZmcn9999Penq648tms9GrVy+2b98OFG7bMBgMbN++vcSwWZGtW7disVj429/+Vmyl7Zgx\nY2jcuDHfffddsfbe3t6OpAOg1+vp2rVrsVW5Go3GkXQKCgq4fPkyWVlZ9O7dm4MHS/4sbrnlFkfS\nAejQoQO+vr6cOXMGKExk33zzDbfffrvTpAOwZs0aevfuTaNGjYr9XPr16wfAjh07nL62NqqfQ21X\nsUd2IH/CNEwfve64psm8hOlfL5L3/+aB0asaoxOi9mvRooXjv0+eLBzavvfee0ttW1Say2g0MmvW\nLGbOnElUVBTdu3fnjjvuYPTo0dxwww0AnDt3DihZ9USj0dCyZUvH/SKhoSWHDP38/Dh8+HCxa59/\n/jkLFizg+PHjxcqDaUpZeNS0adMS1/z9/R3JMi0tjezsbNq3b1/q5y1y8uRJDh8+TKtWrUrcUxSF\ntLQ0l6+vbep94gGw9hlIwe+nMXyzzHFNe+Y4po/fwPzky7LSTVQbW9toNNU41ObJRmx3vLz+/Mdb\n0UT5woULS00EV5s4cSJDhgxhw4YNbN26lbfffpu5c+fy5ZdfEhMTU+Y4ipLata5OLl999RVTpkxh\n8ODBPP300wQHB6PVavniiy9YsWJFuZ7pCVVV6devH08//XSp95s3b16m59V0knj+UHD/Y2h+P4Pu\n4C7HNd2+bRjWLqbgvr9WY2SiPvOk9JP50elYbxlc5mfrtm/E9Mmbbt+/IrVs2RKAwMBAxzCSK82b\nN2fixIlMnDiRCxcucMstt/DOO+8QExNDs2bNADh+/LjjuVCY3E6ePEl0dNljX7t2LS1atGDp0qXF\nri9ZssSjAzCvFRQUhK+vb4le1bVatGhBdna2Rz+TukD+KV9Eo8H8txexNW1Z7LJh3WJ0e6r+2G4h\nwLNf/OVNDpX5bGcGDBiAn58fc+fOxWKxlLifnp4OQF5eHnl5ecXuhYWFERQURFZWFgD9+/fHYDDw\n4Ycflui1pKamMmjQoDLHp9PpSlTgT0xMZP369WV+FhQOz919991s3ryZH3/80Wm74cOHc+DAAf77\n3/+WuJednU1BQUG53r+mkh7P1by8MT/1Gl6xE9Fk/zmhafz4DeyNw7C3aFuNwYn6qNTST1TMJs/K\nfLYzvr6+vPvuuzz22GP07duXESNGEBwczLlz5/j+++9p164dCxYsICEhgXvuuYfhw4fTpk0bjEYj\n3333HcePH+fVVwsrjwQGBjJt2jRee+01hg8fzpAhQxxnh3Xq1ImHHnqo+Of1YPhr8ODBxMXFMWbM\nGAYPHsyFCxf49NNPiYqKKrHPyJlr3+fll19m69at3H333Tz88MO0adOGlJQU1q9fz5IlSwgPD+cf\n//gHmzZt4oEHHmDs2LF07tyZvLw8fvvtN77++mt27drl6OHVBZJ4rqEGh2L+xz/xeuMZFJsVAMVS\nULjY4OWFqI2CqzlCUd+owaFYg0PLNZxWnc++tq5jkeHDhxMSEsLcuXP54IMPyM/PJzQ0lJ49ezpW\nwDVr1oxRo0bxww8/sHLlShRFITIykg8++IAHH3zQ8aypU6cSGBjIRx99xMyZM/H39+fBBx/k5Zdf\nLrYk2lks114fO3YsaWlpLFq0iG3bttGyZUtmz57NyZMn+fXXXz3+3Fdr0qQJmzdv5rXXXmPVqlVk\nZmYSGhrKgAEDCAwMBMBkMrF+/Xrmzp3L2rVrWb58OQ0aNCAyMpLnnnuOxo0be/TetUW93sfjiu6H\njZgWFR//trVoQ96MeWA0VVNUQghR+8kcjxPWvoMpuLP4iana08cwLnoTyrhiRQghxJ8k8bhQMPoJ\nrJ17Fbum3xOP/uv/VFNEQghR+0nicUWjxTxxJrawiGKXjas/RbtvW/XEJIQQtZwkHne8fDA//Tpq\ng4bFLps+eh1N4vFqCkoIIWovSTweUBuHkff3f6JetUtZKcjH9K8XUDLSqzEyIYSofSTxeMjeNpr8\nh4qXs9BcTsM070UoyK+mqIQQovaRxFMG1lvvpmDgiGLXtKd+w/jp27LSTQghPCSJp4wKxkzE2qlH\nsWv63ZvRr/+imiISQojaRRJPWWl1mJ98GXtoeLHLxpWfoN2/vZqCEkKI2qNMiWfu3LkEBAQwbdo0\nl+0OHz7MkCFDCA0NpX379rz11lvXFWSN492AvKdfR/XxLXbZtCAW07sz8H52NN7Pjsb48Wx02zde\n17n1QghR13hcq23fvn0sXryYDh06uCwPnpWVxfDhw4mJiSE+Pp5jx44xefJkvL29nR5hWxupTZpi\nnhyLac40FJsNAMVmRXdwt6ONZse3jpL2V+Ysq/CCi0IIURt51OPJzMzk8ccfZ/78+fj7+7tsu2LF\nCsxmMwsXLqRt27YMGzaMKVOmsGDBggoJuCaxte9K/rh/eNS2uo8wFkKImsKjxPPUU09x7733EhMT\n47a0+N69e+nduzdGo9FxbcCAASQlJXH27Nnri7YGsg4Yhu2a+Z7SSOIRQohCbhPP4sWLSUxM5MUX\nXwRKlvy+VkpKSokS3sHBwY57dZHiwT4eSTxCCFHI5RxPQkICr7zyCps2bXKcLX7t6XzXKs/xsEXv\nVVu1t1oxumljsVhr9WcUQoiyiIqKcnrPZeLZu3cv6enp9Or1Z4Vmm83G7t27+eyzz7hw4UKxw5YA\nGjduXKJnk5qa6rhXniBrOk2n7vDHIgJXbWrzZxRCiIricqjt7rvvZvfu3ezYsYMdO3awfft2unTp\nwv3338/27dtLJB2Am266id27d5Of/+fwU3x8PGFhYYSHu58LqY08OZfe3rx1FUQihBA1n8vE4+fn\nR9u2bR1f7dq1w8vLy3EdIDY2lmHDhjlec//99+Pl5cWTTz7pOC983rx5PPnkk5X7SaqRJ4lHt/Nb\nsBRUQTRCCFGzlblywbVnlCcnJ5OYmOj4vmHDhqxZs4akpCT69+/P9OnTmTx5MpMmTaqQgGsiNTiU\nK3OWYX50OpaYQdgDm6CavIq10SYex7j4XanpJoSo95SMjAz5TVgZrFZM7zyH7siBYpfzx0zEMnh0\nNQUlhBDVT2q1VRadDvOkWdib3FDssmH5/6G9qrqBEELUN5J4KlODhuQ99Tqqt4/jkqKqmBa+gub8\n6WoMTAghqo8knkqmhjXHPPElVOXPH7VizsX0rxcgO6MaIxNCiOohiacK2G7sScEDxVf1aVKT8Hr/\nZbBaqikqIYSoHpJ4qojljhFY+t1d7Jr22M8YP/+XrHQTQtQrkniqiqKQ/9AUbG07F7us3/YN+u9W\nVlNQQghR9STxVCWdnry//xN7cFixy4ZlC9H+sqeaghJCiKoliaeqNfDD/NRrqCZvxyVFtWNa8E+U\nC2eqMTAhhKgakniqgb1pC8xPXrPSLe8KXu/OgJzMaoxMCCEqn1QuqEb6TV9hXFb8ZFZr22jM0+aA\nzuNTyYUQwqmkXBsH0wo4mG7hYFphvcjoIAPRgXqigwyEemurPCZJPNVJVTF++jb6HzYUu2zpP5T8\n8c9AOc82EkIIKEw6Yzenu2yz7PbAKk8+MtRWnRSF/Ieewta6U7HL+vg49JvXVFNQQoi6oqiHc71t\nKpoknuqmN5D391ewB4UUu2xY+gHaX/dXU1BCiLrgYLr7DeqetKloknhqgob+mJ9+vdhRCordjmn+\nyyhJZ6sxMCFEbbYvJd9tG+nx1GP2pi0xT5yJetW8jpJ7Ba93X4CcrGqMTAhR25zPsfLSvkwu5dfM\nKXxJPDWILfpmCkY9UeyaJvk8pvmzwGqtnqCEELVGRr6d9w5lMz7+Ej8kue/tQOEKt6omiaeGsQwe\njaXPoGLXdEcOYFj6QTVFJISo6fJtKssSrvDg9+msPp2HrQwdnehAfeUF5oRsFqlpFIX8vz6LJvl3\ntCd+dVw2fL8W+w0tsN42rBqDE0LUJHZV5fvf8/nktxyS8+zlekZ19HhkH08NpWRdxmvW39CkJzuu\nqRoN5qlvY+vQrRojE0LUBD+lFbDwcA7HM0sfhg/x1vBYuwa089fxc7pFNpAKz2jOnsTr1Uko+WbH\nNdXHl9yXFqKGNK3GyIQQ1eVMtpUPj+SwK7n01WgN9Ap/ifJheAsvDNqauQldEk8Npz2wA9N7M1Gu\nOrPHHtqM3JkLwMe3GiMTQlSlS2Y7nx27wvqzedhL+a2tU2B4Cy/GtfbBz1Czp+/dRvfxxx/Tp08f\nwsPDCQ8PZ+DAgXz33XdO2585c4aAgIASX1u2bKnQwOsLW9cYCkY8WuyaJukcpgX/BJusdBOirjNb\nVf5zvHA3BByxAAAgAElEQVThwNdnSk86t4YZWTygEZM6+tb4pAMe9Hg2bNiA0WikVatW2O12li5d\nyrx589iyZQudOnUq0f7MmTNER0ezevVqOnbs6Lju7++PXl/1qyfqBFXF+NHr6Hf9t/hlkzdqA19s\nbaMdX2pwaDUFKYSoSDZV5btzZhYdvUKaufSFAx0D9Ezs0IAOjWrX79ZyDbW1aNGCWbNmMX78+BL3\nihJPfHw80dHRFRKkAAry8XplEtqzJ1w2uzJnmSQfIWq5/SkFLDySw8ms0kc1wry1PNHeh76hRpRa\nWEy4TMupbTYba9euJT8/n5tvvtll23HjxpGfn0+rVq2YOHEiw4bJMuDrYjBijRmEdqnrxKM9ehCr\nJB4hajRnRxXc4K3lx7QCfnZSP62hXmF8Gx/uifBCr6l9CaeIR4nn8OHDDBw4kPz8fLy8vPj3v/9N\nVFRUqW19fX159dVX6dWrF1qtlg0bNjBhwgQWLlzIqFGjKjT4+kbjprcDoP3tJ6y3DK6CaIQQ5eHs\nqIJvz5lLaV1Ir4ERLbx5sLU3vvqaP4fjjkdDbRaLhd9//53MzEzWrVvHRx99RFxcHF26dPHoTaZN\nm8auXbvYuXPndQdcn3k/OxpNWrLLNqqiwXLPOCz97kINbFJFkQkhPLXxbB5vHsz2uP3tNxh5pF2D\natlvU1nKNcdz7733EhYWxoIFC9w3BpYuXcqzzz5LUlKS0zYJCQllDaPeaf/+8xgzXR/qVERFISuy\nI2ld+5IV2Qk0decPrRC12acXvNid6b5aQGtvKyMbm4nwslVBVBXP2agYlLNkjs1mw273vDzDoUOH\nCAkJcdnGVZCikKZTd9jxrUdtFVT8ThzC78Qh7AFBWPve9UcvqHElRymEcCbdbOPnhEuA63/v+xsU\nPrwttFYuHPCE28Qza9YsBg0aRFhYGDk5OaxcuZKdO3eyatUqAGJjYzlw4ADr1q0DCns3BoOBTp06\nodFo2LRpE4sWLSI2NrZyP0k9YGsbjd7DxHM1zeU0DOsWo//6P9g698Ry61BsnXtKL0iIKnIx18aX\nJ3L55mweFg/+zW7UKnU26YAHiSclJYXHH3+clJQUGjZsSMeOHVm1ahX9+/cHIDk5mcTEREd7RVGY\nM2cO586dQ6vVEhkZyfz58xk5cmSlfYj6wtbW/fL0/BGPoNv/A9ozJYcuFdWO7uBudAd3Y28UjKXv\nXVj7DUFtJL0gISrDuRwrSxNy+e68uWwVo6uhcGdVkpI5tYySmoT26EHHF1ByA6mqojl9DH381+j+\ntwWlwPlqGVXRYOvcC0v/odhuvEl6QUJUgJOZVpYkXGHbhXzKUzN6erQvg8O93DespSTx1HV5V9Dt\n2ox+69doz5502dTeqDGWfndh7TsEtVGwZ0lOCOFw5LKFJcevOC3gCRBkVEhzczLostsD69QqtmtJ\n4qkvVBXNqaPot8Z51gtq1wXdkR9dPlKqJAgBqqpyML0w4fyYVvrGT4CmPloejPLm9qYm0sz2UjeQ\nVudRBVVJEk99lJuDbvdm9PFxaM+57gW5Yn50umxWFfWWqqrsSSlgyfFcfr3sPOG09NUyrrUP/cKM\naOvwgoGykMRTn6kqmlO/oY+PQ7dnC0qBZ2e0F7HEDCL/sRmVFJwQNZNdVdmelM+ShFwSnBzCBtDO\nX8e41j70bmJAIwmnGEk8otCVbPS7N6OLj0N7/pRHL7EHNSH3neWVHJgQVcNZ/bSi4a9gk4Ytv5v5\nIiGXMznON3V2DtTzl9Y+dAvS1+kl0ddDEo8oTlXRnDyC11vPFjv5tDT2wCbkzpXEI2o/Z/XTrhZs\nUkg1O/912bOxgXFR3nQKrNtLoStCuSoXiDpMUbBHdsDao5/7zaqWApSsy6gNA6omNiEqSVEPxxVn\nSadvqJEHo7xp41+7zsSpTrW/zKmoFJ5sVtVkXcbr5SfQnD5WBREJUXkOOjmGwBmNAnc0NfJZ/0b8\ns4efJJ0ykh6PKJUniQdAcykFr9f+Tv5fp2LtM7CSoxKicnjS4wHQKTA43MTYSB/CfOr2kufKJHM8\nwqlSN5C2bI/m7Am0F8+VaF8w8H4KxvwNtPLvGVF7nM+x8ti2S+S5KQLtrVP4rH8jGntJwrleknhE\n2VmtGJYvxPDdqpK32nfF/ORL4OtfDYEJ4bmETAtLE3I9LmszqJmJGV0aVnpc9YEkHlFuuh3fYvxs\nDoql+Pi4PagJ5n+8ir25HHUhahZVVfk53cLSE7nsTfFseK1IXa+fVpUk8Yjrojl9FNN7M9FcSi12\nXTUYyZ/wHNbet1VTZEL8ya6q7LpYwNITVzhy2fmmT1fqev20qiSJR1w3Jesypg9moT32c4l7BXeO\nomDU4zLvI6qF1a6y+byZZSdcb/ps4atlSLgX3jqFXy7Vz/ppVUkSj6gYViuGZfMxbF5T8laHboXz\nPg38qiEwUR/lWVU2nM1j+clcUvKcz+B0DNDzQJQ3vZsYpMpAFZLEIyqU7oeNGBfPRbFeM+8THIr5\nH69gD4+spshEfZBVYGfN6TxWnc4lq8D5r7ZeTQw8EOnNjVJloFpI4hEVTnPyN0zvz0RzOa3YddVg\nJP/R6Vh7DqimyERdlZJnY8XJXOLOmDE7OepTo8CAMCNjI31o5SdDv9VJEo+oFEpGOqb5s9AeP1Ti\nXsGQsRSMfFROOxVuuSvcWWBT+fJELv89b8bq5DeZQQNDwr0Y1cpbNn3WEJJ4ROWxWjB+8QH6LetK\n3urYA/PEmdBA9kWI0nlSuNMVH53C8BZe3NfCm0YmqQ5Wk0jiEZVOt3U9xv/MK2XeJwzzlFexN2tZ\nTZGJmmzj2TzePJhd5tc1MmoY1cqLoc298NFLwqmJJPGIKqE5cRjT+y+hySj+L1jVaML82PPYetxa\nPYGJMnM3/FXeZccWu0qa2U5qno1Us52vTuRyzMVBa9e6wUfL2Ehv7mhqwqiVFWo1mdvE8/HHH/PZ\nZ59x7lxhba62bdsydepUBg50XhDy8OHDTJs2jZ9++omAgAAefvhhnnvuuYqNXNQ6SkY6pvdfQnvi\ncIl7BXc/SMGICTLvU8N5MvxV2kbLPKtKqtlGWp6dVLON1Dw7qWY7aUX/nWfjsotVaK5E+el4INKb\nvnK0dK3hNvFs2LABo9FIq1atsNvtLF26lHnz5rFlyxY6depUon1WVhbdu3cnJiaG5557jmPHjjF5\n8mSmT5/O5MmTK+2DiFrCUoBxyfvot8aVuGVr2hJ70whHYrK1jXZ8qcGhVR2pKIUnw199mhjwM2r+\nSC420sx2ciyVM7ASYFBYPShI9uDUMuUaamvRogWzZs1i/PjxJe4tWrSI2NhYEhISMBqNAMyZM4dP\nP/2UI0eOXH/Eok7QxccVzvvYPBtKuTJnmSSfaqaqKi/tz2R7UtlqnFUmKdxZO5VpMbvNZmPt2rXk\n5+dz8803l9pm79699O7d25F0AAYMGMBrr73G2bNnCQ8Pv76IRZ1g7T8Ue9MWhfM+mZfcttcePYhV\nEk+VsdpVzubYSMi0kJBpJSHTyolMK1ecrVmuAAqFCwOCvDQEmzTkWlUOpLk+oC06UA5gq408SjyH\nDx9m4MCB5Ofn4+Xlxb///W+iokqvPJySkkLTpk2LXQsODnbck8QjitijOpIX+xFeMx9Bk53psq32\n6EGstwyuoshqp/JO+ufbVE5lWf9IMIWJ5lSWlQJPzgrwkFaBIJOGYC8twabCxBLspSXYS0OwSUuQ\nSUOgSYNO8+eQmSfzSdFBUnmgNvIo8bRu3ZqdO3eSmZnJunXreOSRR4iLi6NLly4l2spYqygLNSAI\nDCbAfeIRzjn7Jf3tOTPfnjMDhZP+vnqFE5lWRy8mIdPCmRwb9grsyLT01XJPhBdBpj8Ti79RQVPG\n3w2h3lqW3R5YKSvoRPXyKPHo9XoiIiIA6Ny5MwcOHODjjz9mwYIFJdo2btyYlJSUYtdSU1Md95xJ\nSEjwNGZRx7S3WjG6aWM158ufERd2ZugBb5dtHo9PJdt2ffta9IqKRXWdQPo2yKaD9RLkADmQTuFX\neUUCkQ3g/gZXXcyHnN9B/kTUXM5GxaCMczxFbDYbdnvp/fCbbrqJWbNmkZ+f75jniY+PJywszOUw\nm6sgRd2m6dQddnzrso0+N5u22cnYusZUUVS1y8qfsgCzyzZlTTqNjBqi/HRXfelRUBn7ves5uUEd\nmklPRLjkdlXbrFmzGDRoEGFhYeTk5LBy5UrmzZvHqlWr6N+/P7GxsRw4cIB16wrLomRlZdGjRw9i\nYmKYOnUqCQkJjuXUkyZNqpIPJWoX3faNmD5506O2BfeOp2DYeNDIjnSAjHw7u5LzmfdLNvnXMScT\n6q0hyk9fLNEEmkpPHpW1gVTUH257PCkpKTz++OOkpKTQsGFDOnbs6Eg6AMnJySQmJjraN2zYkDVr\n1jB16lT69+9PQEAAkydPlqQjnLK1jfa4rWHtYjRnTmB+4gXw8qnEqGquC1ds7LiYz46kfH69ZKEs\n+UYDhPtqHT2YKD8dkX46fMtQWibUW0touJccAy3KTUrmiBpBSU1Ce/Sg4wvA1iYa7DZ0//seRS3+\nx9Qe2oy8f7yKGta8OsKtUqqqkpBpdSSbU9nOT9J0pVdjA7O6+2HSyQIgUb0k8YgaT3P0IKYPZqHJ\nzih2XTV5Y/7bi9i6lL6nrDaz2lUOXbKwIymfHRfzSXZxiqanpkf7Si9F1AiSeEStoKSnYHrvRbSJ\nx0vcyx/+Vyz3/KXWz/vkWVX2pxawIymf3cn5ZHlYZiayoY7OQXpWncpz2a60GmpCVAdJPKL2KMjH\n+Nlc9DtLroCzdo3B/PiMGjfv424i3kursDu5sFezL6XAo02bGuDGQD0xoUb6hBgdyUQm/UVtIYlH\n1C6qiv6/qzAsW4ByzZJ+e1hz8qa8ihrSrJqCK86TnfcK4MlfQKMWegQbiAkx0quJEX9j7e7difpN\nEo+olbS//YRp/iyUa0rtqN4+mJ94EVt072qK7E/lPcisSEO9Qu8QIzEhRnoEG2RRgKgzJPGIWktJ\nu4jpvZlozxTfv64qCgX3TcAydBxUQwknVVU5k2PjzZ+y+C3D84PMAJp4aYgJLUw2nRrpi9UuE6Ku\nkMQjard8M8Z/z0G/e3OJW9bufTE/+jx4uS4lUxFS82wcSCvgx1QLP6YWkF6G3ZytGuqICTEQE2ok\nsqFO6h2KOk8Sj6j9VBX9tysxLF9YYt7HFhaB+alXUZs0dfLi8smx2DmYZuHHtAJ+TC3gbE759tYE\nmTSsHBhUobEJUdNJ4hF1hvbwj5gWxKLkZBW7rno3wDxxJrYbe5b72QU2lcOXC3szP6YVcOyytUwV\nA5yRg8xEfSSJR9QpSmpS4X6fsyeLXVcVhYL7H8Vy1wMk5dndLju2qyons6yFiSbVwi+XCsgvQ6fG\npIUwb63bKgOyqVPUR5J4RN2Tb8a46C30e7aUuJXdtR9Dmj1Cns7k9OW9Ghs4kmEhq8DzvxoaBdr7\n6+kWrKdrsIH2AXrSzHa3y6llU6eojyTxiLpJVdFvXI7hq49Q1OKDYqcbhPF1RH9aZp+nW+oRAH4M\nbs+PQR34Mbg9ST7Oz426WoSvlm7BBroFGegcqMenlEKbsqlTiJIk8Yg6Tfvr/sJ5nyue76cZeucH\npSafIJPGkWi6BeudHhsghHBNEo+o85SUC5jmvYj2/CmP2s/q9iTrI27FR6fQJUhfmGyCDTTz0cpS\nZyEqgCQeUedZ7Srfn8qg8edvc8uZXW7bb4roR+BTLxLlp5MNnEJUgnIdfS1EbWC1q2w+b+bz47lc\nyLVBtylsufAzDS1XXL6uZ/oR9AH6KopSiPpHEo+oc2yqyvfn8/n8+BXOX7lqObOikKP3dpt4/HLS\nsaz7HGvfIagBsrlTiIomQ22izrCpKvG/57P4+BXOOakkMGv/fO4+s82j56kaDbbom7H0H4qtY49a\nf96PEDWFJB5R69lVla0X8ll87ApnnCQcrQJDwk1MTNtO48/fLvt7BIVgufVurLcMRvUPvN6QhajX\nJPGIWsuuqvyQVJhwTjupEKBRYHAzE+Na+xDqrUVJTcJn6thyv6eq1WLr0gfLrUOxdegmvSAhykES\nj6h1VFVlx8UCPjt2hZNZpR87oFFgUFMTf2ntQ5hP8f02SmoS2qMHHV8AtrbRhV/No9Ad/RldfBza\nC4ku47AHh2LpdzfWW+6UXpAQZSCJR9Qaqqqy84+Ec8JZwgHuaGriL629adrgOtbOqCqahF/Rx8eh\n2xePYrE4b6rVYusaU9gLat9VekFCuOE28cydO5e4uDhOnjyJwWCge/fuvPzyy7Rr187pa86cOUN0\ndHSJ66tWrWLAgAHXH7Woc1yWlgnUczrbxmfHrnA803nCua2pkYda+9DsehJOaXKy0O/6Dn18HJoL\nZ1w2tTcOK5wLirkT1a+R695V22jU4NCKjVWIWsBt4hkxYgQjRoyga9eu2O12Xn/9dfbt28eePXvw\n9/cv9TVFiWf16tV07NjRcd3f3x+9XvZHiOKScm1ui2k6owADbihMOM19K3l3gKqiSTj0Ry9oq5te\nkA5rh27of9nj8pFX5iyT5CPqHbd/U1etWlXs+w8//JDw8HD27NnDoEGDXL42ICCA4ODg64tQ1HlF\nPZyy6h9mZHwbHyIqO+EUURTsrW8kv/WN5D84Gf3OP3pBSWdLNrVZ3SYdAO3Rg1gl8Yh6psx/Y7Oz\ns7Hb7U57O1cbN24c+fn5tGrViokTJzJs2LByBSnqtoPpznsOpekbauThNj60bFiN+58b+GEZNBLL\nwPvRHPsF/dY4dPu3uewFlUZ79CDWWwZXUpBC1ExlXlzw8MMPc/r0abZu3eq0YOKlS5dYtmwZvXr1\nQqvVsmHDBt555x0WLlzIqFGjKiRwUXeM/C6NVLP78zxvCTEyvo03kX41dLg2JxP9jm/Rb41Dk3TO\no5fYGwWT++6KSg5MiJqlTInnhRdeYO3atWzcuJHmzZuX6Y2mTZvGrl272LlzZ6n3ExISyvQ8UbvZ\nVTh8RcfWywZ+ydFROFvjnJ/Ozpwoz482qFaqSoOzx2n55XtoLa6HEVUgo30P0rr2I6d5a5Dq16KO\niIqKcnrP47GKGTNmsHbtWuLi4sqcdAC6dOnCkiVLnN53FaSoOy7n29lwNo+4s3lczHXfyynSK9Sb\nqKiQSoysgrVujf30IbQ7vnXZTAECjuwj4Mg+7CHNsNx6N5aYQeDrfihbiNrKo8Qzffp01q1bR1xc\nHJGRkeV6o0OHDhESUot+cYgKo6oqv1yysC4xjx8u5GMtx86x6MAaOrzmgq1tNHo3iedqmovnMH65\nEMPKT7B274ul/1DsbTpLL0jUOW4Tz9SpU/nqq69YsmQJDRs2JDk5GYAGDRrg4+MDQGxsLAcOHGDd\nunUALF26FIPBQKdOndBoNGzatIlFixYRGxtbiR9F1DQ5FjvfnTPz9Zk8Ep2UtPFUdJChgqKqOra2\nJfeyeUKxWtD/73v0//see2gzLLcOLewFNfCr4AiFqB5u53gCAgJQFAVVLd7s+eefZ/r06QA8+eST\n7Ny5k59//hmAZcuWMW/ePM6dO4dWqyUyMpKJEycycuTISvoYoiY5lmHh68Q8vv/djNlFvjFooP8N\nJoZFeBFgUAo3j5a2gTTIQKh37Txm2uUG0jad0aRdRLc1Dt3+7Si20jfHAqh6Pdbu/bD0vwd7607S\nCxK1mpTMERXCbFWJv2BmXWIeRzOc/wIFaOqj5Z4IL+5sZqKhQcrLAChZl9Ft31S4Ii7lgsu29rDm\nhXNBfQZBg4ZVFKEQFUcSj3DKZRmbP3ohZ7KtfH0mj2/PmcmxOP+jpFEgJsTIsAgvugbpnS7Fr/fs\ndrS//YQuPg7dge0oNuddRlWvx9rj1sK5oKhOKGkXpTyPqBUk8YhSeVLGpp2/jt/c9G6CTRqGNvdi\nSHMTQabaOVxWXZTMS+h2bEIfvx5NqutekK3xDWhTfnfZRsrziJpCjr4WpfKkjI2zpKMAPRobuKe5\nF72aGNBppHdTHqpfIyx3PYBl8Bi0Rw6gj/8a7U87S+0FuUs6IOV5RM0hiUeUqjz10/wMCkPCvRja\n3KvEGTjiOmg02Dp2x9axO0pGeuFc0Lb1aFKTyvQY7S97pDyPqBFkqE1gsaskZls5kWklIdPK8Uwr\nhy9Z8PQPRqdGeoZFeNE31IhBK72bKmG3oz38I/qtcWgP7ECxe7YZ1+4fhL1FG2wRrbH/8SWH2Imq\nJomnlvFkwt+VPKvKqSwrxzMtJGQWJpvT2VYsnhcRcPDSKsy/JaB6i3UKlIx0vGeMR8nNKdfr7QFB\n2CP+SEYt2hQmI79Gzt9PzhgS10kSTy3iyYT/stsDHcknq8BOwh+9mIQ/Es25HJvHPRl3BjUzMaOL\nLOetCYwfzy5TlQR37AFFPaM22CPaYG/RGrVhAEpqEj5Tx7p8rSxiEO7IP1VrEU/mXeb/WlhIMyHT\nSnJeOboxZVAby9jUVZ6U51G1WpfLs6+muZyG5nIaugN/FvW1NwpG9aCGnCxiEO5I4qlFPDm3ZsfF\n8h2qFmTS0NpPR5Sfjkg/Pf4Ghb/vzHD5mtpYxqau8qQ8T+7ri1Hy89AkHkebeBzN6WNozp/0+Awh\nzaVUuJTqtp3ux+1Yu8aAj69HzxX1jwy11SKj/5tWIb2Ypj7aPxKM7o9ko8ffWLKCwPXOJ4mqVa65\nF6sVze+n0Zw+VpiMEo+hOXcKxVq2A+1Ko3r5YA9sghrUpPB/r/7voBDUhgGgcV25QuaT6iZJPLXE\nyUwrk3Zccln77FpaBSJ8r04wOlo11OGjlzI1wgWrBc3viX8ko2NoEo9XWDK6mqrTozYKxh4UUpiU\nAhtjDwxxJCdUOz7T/+LyGTKfVDtJ4qnBVFXlQJqFL0/ksi/VsyG0Zj5aRkd6E+WnI8JXh1GWN4uK\nYLWgOX8a4+fz0J48XCVvqeLueEAwPzpd9ibVQjLHUwNZ7SrbLuTz5clcEjJdl6S51gNR3gwO96qk\nyES9pdNjj2iNpf/dbhOPvWEASl6Ox3NHznjyTybt0YOSeGohSTw1SK7VzoazZlaczC33XI5M+IvK\n5MkihryXFqAGNkHJzkBJS0ZJT0aTdrHwf9NTUNIvoklLLve+o6vp9v1AwdBxqCHNrvtZourIUFsN\nkG62seZ0HusS88h2UeG5Y4Cegc1M6DTws0z4i2pSYRP+eVcKE1B6UXJK/jMppaegZKR51OtRFQ3W\nm27FMnQc9mYty//BRJWRxFONzuZYWX4il+/Om51WDlCAPiEGxkT60LGR7JsR9Yfxw9fQ7/pvmV5j\n7dKHgqHjsLdqV0lRiYogiacaHEov4MuTuex0sedGr4E7m5kY1cqbZg1kRFTUP7rtGzF98ma5Xmvt\n0J2Ce8Zhb9NZTmutgeQ3WhWxqSq7Lhbw5YkrHL7sfMGAr17h3hZeDI/wppFJlj2L+suT+SRbSFO0\nF8+XuK47vB/d4f3YojpSMPQv2G68SRJQDSI9ngribLNlxwA9GQV2vj1n5vwV55twQrw1jGpZuCLN\nSyd/QYQAD+aTApugPbgLw9dL0J4+6vQ5tuatKbhnHLauMW43rYrKJ4mnAnhSvNOZ1n46xkR60zfU\nKAemCVFeqor28I8Y4v6D9ujPTpvZwiKwDH0Qa8/+oJUBn+oiiacCbDybx5sHs8v0mp6NDYyJ9CY6\nUI8iQwBCVBjN8V8wfL0E3aG9TtvYg8MouPsBrH0Ggl62IFQ1t4ln7ty5xMXFcfLkSQwGA927d+fl\nl1+mXTvXq0YOHz7MtGnT+OmnnwgICODhhx/mueeeq9Dga4pXf8xk8+/5btvpFLitqYnRrbzlDBsh\nKpnm9DEM679At/8Hp23sAUFYhozB0u9ulKzLUheuirhNPCNGjGDEiBF07doVu93O66+/zr59+9iz\nZw/+/qWXSM/KyqJ79+7ExMTw3HPPcezYMSZPnsz06dOZPHlypXyQqqaqKkczrGw4m8f6M2a3Z9x4\n6xQ+69+Ixl6yz0aIqqQ5fxr9+i/Q/W8Lilr6vgW7jy+aK65HLaQuXMUp81DblStXCA8PZ+nSpQwa\nNKjUNosWLSI2NpaEhASMRiMAc+bM4dNPP+XIkSPXH3U1ysi389/zZjaczeN0tucVO5t4aVh+R1Al\nRiaEcEVJ/h3DN8vQ7diEYitbKSqQunAVqczLO7Kzs7Hb7U57OwB79+6ld+/ejqQDMGDAAJKSkjh7\n9mz5Iq1GNlVlT3I+L+/L5P7v0ph/OKdMSQeklI0Q1U1tcgP5E6aS+/ZSCu4YgWowun/RVYqG38T1\nK/NEw/PPP8+NN97ITTfd5LRNSkoKTZs2LXYtODjYcS88PLysb1stLlyxseFsHpvOmUkzX985OHJa\npxA1gxrYmIJxf8dyzzj0m1ag/34tijnX7eu0v0niqShlSjwvvPACe/fuZePGjS5XYpVnlVZCQkKZ\nX1MZCuxwIFvPjgwDx3Jd/3g0qHRqYOXGBhb+c9HbZVv/7N9JSJAFhELUKF36o217E+0/mIEuP89l\nU+VSCtmfzCE9+hYsfo2qKMDaKyoqyuk9jxPPjBkzWLt2LXFxcTRv3txl28aNG5OSklLsWmpqquNe\nWYOsbKqqcizTyoYzZr7/3cwVq+sE0cxHy+BwE4OamQg0FS4WGOL2tM7SP7cQovqpPfrCjm9dtlFU\nldDt6wnZsQFb555Ybh2KrXNP0MiCobLyKPFMnz6ddevWERcXR2RkpNv2N910E7NmzSI/P98xzxMf\nH09YWFiVDrO5O7rZS6uw+XczG87kccrNnI1JC/3DTAwON9GpUcm9N6HeWkLDveQsHCFqIVvbaPRu\nEk8RRbWjO7gb3cHd2BsFY+l7F9Z+Q1AbyT8uPeV2VdvUqVP56quvWLJkCW3atHFcb9CgAT4+PgDE\nxsZy4MAB1q1bBxQup+7RowcxMTFMnTqVhIQEx3LqSZMmVeLH+ZMn1QR0Crjp3NAhQMfgcC8G3GDE\nW3s6rxoAAAVPSURBVCelNoSoi5TUJHymji3361VFg61zLyz9hxbWhZNekEtuE09AQACKoqCqxZs9\n//zzTJ8+HYAnn3ySnTt38vPPf5aqOHLkCFOnTuXAgQMEBATw17/+tUo3kJanmkARf4PCwGYmhoR7\nEeErGz2FqA9c1oWLaIM24RD6+Di0Z0+4fI69UWMs/e7C2ncIaqPgqgi91qmzJXNm/5TFt+fMHrfX\nAD2bGBgS7kXvJgapmyaEKElV0Zw+hj7+68INqQXOf8eoigZbdO/CXlCnHtILukqdTTyj/5vm0fHR\nTa9aKBBkkj8YQggP5eag2/09+q1foz170mVTe2CTP3tBAUEVd4prLVWvE08jo8KqgUFSpFMIUX6q\niubUb+i3rnffC9JosLXriu7wfpePrOvleersbLknlQJ6NDZK0hFCXB9Fwd6qPfmPPMeVeSvJ/8sU\nbE1blt7UbnebdKDuV0mou4nHg0oBUk1ACFGhvBtguX04ea8uInfmfCwxd5a5NA/U/cRTZ5dsedLj\nkfppQohKoSjYIzuQH9mB/Acmod/1X3Rb49CeP+3Ry7VHDlRygNWrzs7xgPsNpKHesphACFFFVBXN\nySN4vfUsSr7rFbcqYOvet7A6Qodude647jqdeIQQoqYxfjzb4yoJAPbgUCz97sZ6y52o/oGVGFnV\nkcQjhBBVSLd9I6ZP3izz61StFlvXmMJeUPuutboXJIlHCCGq0PWW5wGwNw7DcuvdWGPuRK2FlbIl\n8QghRBVzuYG0TTRKRir6+Dh0+7aiWCxOn6NqdVi7xmDtPxRbuy61phckiUcIIWqqnCz0O79FHx+H\nJsn16c32Jjf8ORfUMKCKAiwfSTxCCFHTqSqaY7+g3xqHbv82972gbrcU9oLaRqOkJ9e48jySeIQQ\nojbJyUS/41v0W+PQJJ1z2dQe1ARNWrLLNtVRnkcSjxBC1EaqiubYz4VzQft/QLE67wW5Yn50OtZb\nBldwcK5J4hFCiNouO+OPXtB6NBdd94KuZYkZRP5jMyopsNJJ4hFCiLpCVdEePYgu/mt0+7ej2Kxu\nX2IPakLuO8urILg/1dlabUIIUe8oCrZ2XbC160J+VgY+08eh5OZUd1Ql1I5F30IIIcqmoT/Wrn3c\nNrO1ja6CYIqTxCOEEHWUJ0lFEo8QQogKU1MTjywuEEKIOsxleZ6auoF0586dvP/++/zyyy8kJSUx\nf/58HnjgAaftz5w5Q3R0yQy6atUqBgwYcP0RCyGEqNXcrmrLzc2lY8eOjB07lokTJ6IoikcPXr16\nNR07dnR87+/vX/4ohRBC1BluE88dd9zBHXfcAcCkSZM8fnBAQADBwcHlj0wIIUSdVGmLC8aNG0dU\nVBR33nkn69atq6y3EUIIUctU+AZSX19fXn31VXr16oVWq2XDhg1MmDCBhQsXMmrUqIp+OyGEELVM\nhSeeRo0aFRuSi46O5vLly8ybN08SjxBCiKopmdOlSxeWLFnisk1CQkJVhCKEEKIKREVFOb1XJYnn\n0KFDhISEuGzjKkghhBB1h9vEc+XKFU6ePAmA3W7n3Llz/PLLLzRq1IimTZsSGxvLgQMHHAsIli5d\nisFgoFOnTmg0GjZt2sSiRYuIjY2t3E8ihBCiVnCbeA4cOMA999wDgKIozJ49m9mzZ/PAAw8wf/58\nkpOTSUxMdLRXFIU5c+Zw7tw5tFotkZGRzJ8/n5EjR1bahxBCCFF7SMkcIYQQVUqKhAohhKhSkniE\nEEJUKUk8QgghqpQkHiGEEFVKEo8QQogqJYlHCCFElZLEI4QQokpJ4hFCCFGlJPEIIcT/3yigKwAA\nz5yD+rd99sAAAAAASUVORK5CYII=\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x107e40588>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"ax = data.plot(style='-o', markersize=10, markeredgecolor='none');\n",
"plt.xlim(-1, len(data))\n",
"plt.ylim(1.4, 5.1)\n",
"plt.xticks([]);\n",
"\n",
"color = ax.lines[0].get_color()"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Here's a better view of the relationship in the data:"
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAcAAAAEtCAYAAACf/7AvAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XtcVHXeB/DPzMhlVMDhLhqiQg8ieCuvkQUmFl3QEu+v\nLppWyK5u4q1wV7InTY3WCnnUqNiU8oKIrpfMQNeohJbHVDQbc0GfQkCuogzqDM8fLrOMzAxncG5w\nPu/Xi9dr5/x+Z+Y7vz3Nx3P7HUlNTU0TiIiIREZq6wKIiIhsgQFIRESixAAkIiJRYgASEZEoMQCJ\niEiUGIBERCRKDEAiIhIlmwbg6tWroVAodP6Cg4ONrlNUVITo6Gj07NkTISEhWLt2rZWqJSKizqSL\nrQu4//778fe//137WiaTGexbV1eHSZMmITw8HLm5uTh//jzi4+PRtWtXxMfHW6NcIiLqJGwegDKZ\nDF5eXoL67ty5EyqVCqmpqXByckJwcDCUSiU2btzIACQiIpPY/BxgcXExBgwYgMGDB2POnDkoLi42\n2Dc/Px+jR4+Gk5OTdllkZCRKS0tx6dIlK1RLRESdhU0DcPjw4UhNTUVmZiY++OADlJWVYcKECaiu\nrtbbv7y8HN7e3jrLmvcey8vLLV4vERF1HjY9BPrYY4/pvB4+fDgGDx6MjIwMzJ8/v1V/iURirdKI\niKiTs/kh0Ja6du2K4OBg/Otf/9Lb7u3t3WpPr6KiQttGREQklF0FoEqlwi+//AIfHx+97SNGjMD3\n33+PxsZG7bLc3Fz4+fnB39/fWmUSEVEnYNMATExMRF5eHoqLi/Hjjz/ihRdeQENDA6ZPnw4ASEpK\nQkxMjLb/5MmTIZfLERcXh3PnzmHv3r3YsGED4uLibPUVOg2lUmnrEjoMjpVwHCvhOFbWZ9NzgKWl\npXj55ZdRWVkJT09PDB8+HEeOHEHv3r0BAGVlZTpXhbq6uiIrKwsJCQmIiIiAQqFAfHy83vOFRERE\nxkj4RHgC7vzrMygoyNZldAgcK+E4VsJxrKzPrs4BEhERWQsDkIiIRIkBSEREosQAJCIiUWIAEhGR\nKDEAiYhIlBiAREQkSgxAIiISJQYgERGJEgOQiIhEiQFIRESixAAkIiJRYgASEZEoMQCJiEiUGIBE\nRCRKDEAiIhIlBiAREYkSA5CIiESJAUhERKLEACQiIlFiABIRkSgxAImISJQYgEREJEp2E4DJyclQ\nKBRYvHixwT4lJSVQKBSt/nJycqxYKRERdQZdbF0AABQUFCA9PR0DBw6ERCJps//u3bsRGhqqfd2j\nRw9Llmf3qlRqXGnQAAB85VK4O8tsXBERkf2zeQDW1tZi3rx5SElJwZo1awSto1Ao4OXlZeHK7J/q\ntgYFFTeRcqZeJwDjQ7vjQS9HOHexmx18IiK7Y/NfyIULF2LixIkIDw9HU1OToHVmzZqFoKAgPP74\n48jOzrZwhfaroOImVhTUacMPAK40aJBYUIeCips2rIyIyP7ZdA8wPT0dxcXF+PjjjwGgzcOfLi4u\nePvttzFq1CjIZDIcOHAAs2fPRmpqKqZMmWKNku1GlUqNj87UG2xPOVOPEIUDPHg4lIhIL5sFoFKp\nxKpVq3Do0CHIZHd+pJuamozuBbq7u2P+/Pna10OGDEF1dTU2bNggugC80qBBWYs9P0PtDEAiIv0k\nNTU1wo47mtm2bdsQHx+vDT8AUKvVkEgkkMlk+P333+Hg4NDm+2RkZGDRokUoLS012EepVJqlZntS\nK/dCQqHxw5zrhznCraHCShUREVlHUFCQWd7HZnuATz31FB544AHt66amJsyfPx+BgYF4/fXXBYUf\nAJw+fRq+vr5G+5hrsOxJlUoNX3m1zvm/lnzlUvT1dIGHs7ArZJVKZaccJ0vgWAnHsRKOY2V9NgtA\nNzc3uLm56SyTy+Vwc3NDcHAwACApKQmFhYXaC10yMjLg6OiIsLAwSKVSHDp0CGlpaUhKSrJ6/bbm\n7ixDfGh3JBbU6W2PD+3Ow59EREbY/DaIliQSic6FMGVlZSguLtZpX79+PS5fvgyZTIbAwECkpKQg\nNjbWBtXa3oNejlg13NXgbRBERGSYzc4BkvlUqtTaC2J85NJ27fnx8ItwHCvhOFbCcaysz672AKl9\nPJxlPNxJRGQim98IT0REZAsMQCIiEiUGIBERiRIDkIiIRIkBSEREosQAJCIiUWIAEhGRKDEAiYhI\nlBiAREQkSgxAIiISJQYgERGJEgOQiIhEiQFIRESixAAkIiJRYgASEZEoMQCJiEiUGIBERCRKDEAi\nIhIlBiAREYkSA5CIiESJAUhERKLEACQiIlGymwBMTk6GQqHA4sWLjfYrKipCdHQ0evbsiZCQEKxd\nu9ZKFRIRUWfSxdYFAEBBQQHS09MxcOBASCQSg/3q6uowadIkhIeHIzc3F+fPn0d8fDy6du2K+Ph4\nK1YsLlUqNa40aAAAvnIp3J1lNq6IiOje2TwAa2trMW/ePKSkpGDNmjVG++7cuRMqlQqpqalwcnJC\ncHAwlEolNm7cyAC0ANVtDQoqbiLlTL1OAMaHdseDXo5w7mI3BxCIiExm81+whQsXYuLEiQgPD0dT\nU5PRvvn5+Rg9ejScnJy0yyIjI1FaWopLly5ZulTRKai4iRUFddrwA4ArDRokFtShoOKmDSsjIrp3\nNg3A9PR0FBcXIzExEQCMHv4EgPLycnh7e+ss8/Ly0raR+VSp1PjoTL3B9pQz9ahUqa1YERGRedns\nEKhSqcSqVatw6NAhyGR3zik1NTUZ3QtsKyDJfK40aFDWYs/PULsHzwcSUQdlswDMz89HZWUlRo0a\npV2mVqvx/fff47PPPsPvv/8OBwcHnXW8vb1b7elVVFRo2wxRKpVmrLzzajlON+Rebfa/cf06lFcr\nLFmS3eI2JRzHSjiOlTBBQUFmeR+bBeBTTz2FBx54QPu6qakJ8+fPR2BgIF5//fVW4QcAI0aMwMqV\nK9HY2Kg9D5ibmws/Pz/4+/sb/CxzDVZnplQqdcapSqWGr7xa5/xfS75yKfp6usDDuYe1SrQbd48V\nGcaxEo5jZX02Owfo5uaG4OBg7d+AAQMgl8u1ywEgKSkJMTEx2nUmT54MuVyOuLg4nDt3Dnv37sWG\nDRsQFxdnq6/Rabk7yxAf2t1ge3xodx7+JKIOzeZXgbYkkUh0zvOVlZWhuLhY+9rV1RVZWVkoLS1F\nREQEli5divj4eMyfP98G1XZ+D3o5YtVwV/jK/7OZ+MqleHu4Kx70crRhZURE905SU1Nj/N4DEgVj\nh18qVWrtBTE+cqno9/x4qEo4jpVwHCvrs/mN8GT/PJxlog89Iup87OoQKBERkbUwAImISJQYgERE\nJEoMQCIiEiUGIBERiZLJV4FeuHAB3377La5evYrJkycjICAAN2/eRFlZGby9vXWe1EBERGSvBAeg\nRqPBwoUL8fnnnwO4c9P68OHDERAQgMbGRowZMwZLlizBH/7wB4sVS0REZC6CD4G+99572LZtGxIT\nE/H111/rPLXBxcUFzzzzDP7+979bpEgiIiJzExyA27Ztw8yZM7Fo0SL07du3VfuAAQNw4cIFsxZH\nRERkKYID8Pfff8eDDz5osF0ul6O+3vADVImIiOyJ4AD09vbGpUuXDLb/9NNPuO+++8xSFBERkaUJ\nDsBnnnkGn3zyCS5cuNDqyexff/01vvjiC0ycONHsBRIREVmC4ABcunQpevfujbFjx2Lu3LkAgPff\nfx/jxo3DlClTEBYWhtdff91ihRIREZmT4AB0c3PDV199hUWLFqGsrAzOzs744YcfcOPGDSxfvhwH\nDhxA165dLVkrERGR2Zh0I7xcLseiRYuwaNEiS9VDRERkFYL3AOvr63H58mWD7ZcuXcL169fNUhQR\nEZGlCQ7AN998EzNmzDDYPmvWLKxYscIsRREREVma4ADMzc3Fk08+abD9qaeeQk5OjlmKIiIisjTB\nAVhWVgY/Pz+D7d7e3igtLTVLUURERJYmOAA9PDxw7tw5g+3nz5+Hm5ubWYoiIiKyNMEBGBUVhfT0\ndBQWFrZq++c//4nPPvsM48ePN2txREREliL4Nohly5bh66+/RlRUFB577DGEhIQAAIqKinDkyBH4\n+PjgzTfftFihRERE5iQ4AH19fZGTk4OVK1di//79+OqrrwDceRTS1KlT8Ze//AU+Pj4WK5SIiMic\nTLoR3sfHB6mpqdBoNLh69SoAwNPTE1Kp4COpOrZs2YLPPvtMe39hcHAwEhISEBUVpbd/SUkJhgwZ\n0mp5ZmYmIiMj21UD2acqlRpXGjQAAF+5FO7OMhtXRESdjUkB2EwqlcLb2/ueP7xXr15466230L9/\nf2g0GmRkZGDmzJnIyclBWFiYwfV2796N0NBQ7esePXrccy1kH1S3NSiouImUM/U6ARgf2h0PejnC\nuUv7/rFFRHQ3kwKwuroau3btQklJCWpqanSeCt8sJSVF8PtFR0frvE5MTERaWhoKCwuNBqBCoYCX\nl5fwwqnDKKi4iRUFdTrLrjRokFhQh1XDXfFwT2cbVUZEnY3gAPzmm2/w/PPP48aNG3BxcWm119XU\n1NTqMUmmUKvV2LNnDxobGzFmzBijfWfNmoXGxkb0798fr732GmJiYtr9uWQ/qlRqfHTG8EOVU87U\nI0ThAA8eDiUiMxAcgImJifD29sbWrVsxcOBAsxVQVFSEqKgoNDY2Qi6X49NPP0VQUJDevi4uLnj7\n7bcxatQoyGQyHDhwALNnz0ZqaiqmTJlitprINq40aFD278OextoZgERkDpKamprWxzH18PHxQVJS\nEl599VWzFnDr1i389ttvqK2tRXZ2NjZv3ox9+/Zh6NChgtZfvHgxvvvuO+Tl5Rnso1QqzVUuWVCt\n3AsJhTeN9lk/zBFuDRVWqoiI7JGhnSRTCd4D7Nevn0We9uDg4ICAgAAAwODBg1FYWIgtW7Zg48aN\ngtYfOnQotm7darSPuQarM1MqlTYfpyqVGr7yau3FL3fzlUvR19MFHs62vejJHsaqo+BYCcexsj6T\nngaRlpaG4uJiC5Zz51ygRmP4MNjdTp8+DV9fXwtWRNbi7ixDfGh3g+3xod15+JOIzEbwHmBOTg7c\n3d0xcuRIPPzww+jduzdkstY/Ru+9957gD1+5ciUmTJgAPz8/1NfXY9euXcjLy0NmZiYAICkpCYWF\nhcjOzgYAZGRkwNHREWFhYZBKpTh06BDS0tKQlJQk+DPJvj3o5YhVw10N3gZBRGQuggPw008/1f7v\nb775xmA/UwKwvLwc8+bNQ3l5OVxdXREaGorMzExEREQAuPMEipZ7nBKJBOvXr8fly5chk8kQGBiI\nlJQUxMbGCv5Msm/OXaR4uKczQhQO2gtifORS7vkRkdkJvgiGOjeefxCOYyUcx0o4jpX1cVoNIiIS\npXZNhUZkbzh3KBGZyqQAPHLkCD766COcPHkS165dg0aj0c7+0jwTTFVVlUUKJdKHc4cSUXsJ/nXY\nv38/YmNjUV5ejueeew4ajQaxsbF47rnn4OTkhLCwMCxZssSStRK10jx3aMt7B5vnDi2oMH5TPRGJ\nm+A9wOTkZAwZMgSHDx9GXV0dPvnkE8ycOROPPPIIiouL8dhjj6F///6WrJVIB+cOJaJ7IXgP8OzZ\ns5g8eTIcHBy0z/9Tq9UAgICAAMyePRvvv/++Zaok0kPo3KFERPoIDkAnJyfI5XIAQLdu3SCRSFBR\n8Z85Gf38/HDx4kXzV0hERGQBggOwb9++uHDhAgDA0dER999/v3aGlqamJhw8eJBTkpFV+cql8JUb\n3oR95VL4GGknInET/Oswfvx4ZGZm4tatWwCAuLg4HDx4EMOGDcOwYcNw+PBhzJkzx2KFEt2Nc4cS\n0b0QfBHM4sWL8corr2jn/3z++efh7OyMPXv2oEuXLliyZAmmT59usUKJ9OHcoUTUXoID0MHBAR4e\nHjrLpkyZwgfRkk1x7lAiai/OBEOdgoezjKFHRCYxeSaYzz//HMXFxaipqUFTk+482hKJBD/99JNZ\nCyQiIrIEwQH4wQcf4C9/+Qt8fHwwbNgwhISEtOrTPC0aERGRvRMcgP/zP/+DsWPHYteuXXBwcLBk\nTURERBYn+DaImpoaTJw4keFHRESdguAAfOCBB6BUKi1ZCxERkdUIDsB169Zh37592L59uyXrISIi\nsgrB5wCff/553Lp1C6+++ipef/119OzZU3tTPPCf5wGeOHHCIoUSERGZk+AA9PLygre3t9FHHvEq\nUCIi6igEB+D+/fstWQcREZFVcap8IiISJZNmgrl58ybS09Nx+PBhXL58GQDg7++PCRMm4Pnnn+ct\nEkRE1GGYdB/guHHjsGTJEpw6dQoKhQIKhQInT55EQkICxo0bh5qaGpM+fMuWLXjooYfg7+8Pf39/\nREVF4fDhw0bXKSoqQnR0NHr27ImQkBCsXbvWpM8kaq8qlRpnq2+hVu6FKpXa1uUQ0T0SvAeYlJSE\nn3/+GSkpKZg2bRqk0jvZqdFosGPHDvzxj39EUlIS3n//fcEf3qtXL7z11lvo378/NBoNMjIyMHPm\nTOTk5CAsLKxV/7q6OkyaNAnh4eHIzc3F+fPnER8fj65duyI+Pl7w5xKZQnVbg4KKm3c9cqla+8gl\n5y48k0DUEQn+L/fAgQN4+eWXMWPGDG34AYBUKsW0adPw8ssv48CBAyZ9eHR0NMaNG4eAgAD069cP\niYmJ6N69OwoLC/X237lzJ1QqFVJTUxEcHIyYmBgsWLAAGzduNOlziUxRUHETKwrqtOEHAFcaNEgs\nqENBxU0bVkZE90JwANbW1qJfv34G2wMCAkw+BNqSWq1GZmYmGhsbMWbMGL198vPzMXr0aDg5OWmX\nRUZGorS0FJcuXWr3ZxMZUqVS46Mz9QbbU87Uo5KHQ4k6JMEB2LdvX+zfv7/VI5CAOzfBHzhwwGhA\nGlJUVIRevXrBx8cHCxcuxKeffoqgoCC9fcvLy+Ht7a2zzMvLS9tGZG5XGjTaB+22p52I7Jfgc4Bz\n587FokWL8Oyzz+LVV1/VhtQvv/yCTZs24dixY0hOTja5gPvvvx95eXmora1FdnY25syZg3379mHo\n0KGt+rb3RnvOYSoMx6m1G3Kvtvtcvw7l1QorVNMxcbsSjmMljKGdJFMJDsDZs2ejsrIS69atw9Gj\nR3XaHB0d8eabb+Kll14yuQAHBwcEBAQAAAYPHozCwkJs2bJF73k9b2/vVnt6FRUV2jZDzDVYnZlS\nqeQ46VGlUsNXXq1z/q8lX7kUfT1d4OHcw8qVdQzcroTjWFmfSfcBLl68GC+99BKOHj2qvQ/wvvvu\nQ0REBDw8PMxSkFqthkaj/8dmxIgRWLlyJRobG7XnAXNzc+Hn5wd/f3+zfD5RS+7OMsSHdkdiQZ3e\n9vjQ7vBwlultIyL7ZlIAAoCnpycmT55slg9fuXIlJkyYAD8/P9TX12PXrl3Iy8tDZmYmgDu3XhQW\nFiI7OxsAMHnyZLz77ruIi4tDQkIClEolNmzYgKVLl5qlHiJ9HvRyxKrhrnfdBiHV3gZBRB2T4AD8\n+eef8euvv+LJJ5/ULvv222/x3nvvoa6uDs8++yzmz59v0oeXl5dj3rx5KC8vh6urK0JDQ5GZmYmI\niAgAQFlZGYqLi7X9XV1dkZWVhYSEBEREREChUCA+Pt7kzyUyhXMXKR7u6YwQhQPKGjS4cf36vw97\ncs+PqCOT1NTUtL6sU4/Y2FhIJBLs2LEDAPDbb79h5MiRcHJygqenJ3755Rd89NFHmDlzpkULJsvg\n+QfhOFbCcayE41hZn+DbIH766Sed+/O2b98OtVqN48eP48SJE3j88ceRlpZmkSKJiIjMzaQb4T09\nPbWvv/76a4wdOxZ+fn4AgKioKFy4cMH8FRIREVmA4AD09PRESUkJgDsTY//444949NFHte2NjY16\nb5InIiKyR4IvgomMjMTmzZvh5uaG48ePo6mpCdHR0dr28+fPo1evXhYpkoiIyNwEB+CKFStw4cIF\nrFixAo6Ojli1ahX69OkDAGhoaEBWVhamTJlisUKJiIjMSXAAent74+DBg6ipqYFcLteZkBoA9u7d\ni/vuu8/sBRIREVmCyTfC9+jResonuVyOQYMGmaUgIiIiazDpSZ5KpRJz585FcHAwPD09cezYMQBA\nZWUl4uLi8OOPP1qkSCIiInMTHICnT59GZGQkjh49iuHDh0Ot/s8z0Dw8PHDu3DneB0hERB2G4ABM\nSkqCj48PCgoK8Ne//rVV+7hx43DixAmzFkdERGQpggPwhx9+wAsvvKD3HCBw56kQpaWlZiuMiIjI\nkkw6B+js7GywraKiotWVoURERPZKcAAOGjQIX331ld62W7duITMzE8OHDzdbYUSkq0qlxtnqWzhb\nfQtVKnXbKxCRUYJvg0hISMBzzz2HP/zhD9rnAZaWluLrr79GcnIylEolkpOTLVYokVipbmtQUHHT\n4PMInbuYdCCHiP7N5KnQFi9ejK1btwIAXnvtNQCAm5sbtmzZgtGjR1umSiIRK6i4iRV3PZH+SoMG\niQV1WDXcFQ/3NHxqgogMExSA169fx5gxY/Daa6+hqKgIubm5+PXXX6HRaNCvXz9ERkbCxcXF0rUS\niU6VSo2PztQbbE85U48QhQMfzkvUDoICsFu3bqitrYWjoyO6deuGp556ytJ1ERHu7OmV/fuwp7F2\nBiCR6QSfPIiKisLhw4ctWQsREZHVCA7AP/3pTygpKcELL7yAY8eO4dKlS6ioqGj1R0Tm4yuXwldu\n+D9TX7kUPkbaicgwwRfBjBkzBgDw888/Y+/evXr7SCQSVFVVmacyIoK7swzxod2ReNdFMM3iQ7vz\n8CdROwkOwCVLlrTZRyKR3FMxRNTag16OWDXc1eBtEETUPoIDcPny5Zasg4gMcO4ixcM9nRGicNBe\nEOMjl3LPj+gemfw8QCKyDQ9nGUOPyIxsevY8OTkZERER8Pf3R2BgIKZNm4Zz584ZXaekpAQKhaLV\nX05OjpWqJiKizsCme4B5eXmYO3cuhg0bBo1Gg3feeQcTJ07EiRMnDD51otnu3bsRGhqqfd1WfyKx\nq1Kpdc4hunNvkkTOpgGYmZmp83rTpk3w9/fHiRMnMGHCBKPrKhQKeHl5WbI8ok6Bc4kS6WdXW/61\na9eg0WgE7c3NmjULQUFBePzxx5GdnW2F6og6pua5RK+0mFGmeS7RgoqbNqyMyLbsKgCXLVuGQYMG\nYcSIEQb7uLi44O2330Z6ejp27tyJsWPHYvbs2dixY4cVKyXqGITMJVrJRyuRSNnNVaBvvPEG8vPz\ncfDgQaP3E7q7u2P+/Pna10OGDEF1dTU2bNiAKVOmWKNUog6Dc4kSGWYXAbh8+XLs2bMH+/btQ58+\nfUxef+jQodpHNOmjVCrvpTzR4DgJ11HG6oa87fPkN65fh/Kq5aYx7ChjZQ84VsIEBQWZ5X1sHoBL\nly5FdnY29u3bh8DAwHa9x+nTp+Hr62uw3VyD1ZkplUqOk0AdaayqVGr4yqt1zv+15CuXoq+nCzyc\nLXMVdUcaK1vjWFmfTQMwISEBO3bswNatW+Hq6oqysjIAQPfu3dGtWzcAQFJSEgoLC7UXumRkZMDR\n0RFhYWGQSqU4dOgQ0tLSkJSUZLPvQWSvOJcokWE2DcC0tDRIJBLExMToLF+2bBmWLl0KACgrK0Nx\ncbG2TSKRYP369bh8+TJkMhkCAwORkpKC2NhYa5ZO1GFwLlEi/SQ1NTVNti6CbI+HX4TrqGNVqVJb\nfS7RjjpWtsCxsj6bnwMkIuvgXKJEuhiARNQunFqNOjoGIBGZhFOrUWfBLZWITMKp1aizYAASkWCc\nWo06EwYgEQkmdGo1oo6AAUhERKLEi2CISDBfuRS+cqnRqdV85Pr/Xc2rRsneMACJSLD2TK3Gq0bJ\nXnHLIyKTNE+t5ttiT89XLsXbw11bTa0mlUp51SjZLe4BEpFJnLtI8XBPZ4QoHNqcWq1J7trmVaMD\nFQ4AwMOjZHUMQCJqFyFTq9VoHFDWYHgvz91Ziv+tvIUtZ3l4lKyPWxcR2UT3LhKM7emEVf/k4VGy\nDQYgEVlMD+ktnXOFLY3v7YysfzUYXJc31ZOlMQCJyGIkDXWID+2ut62Hk5Q31ZNNMQCJyGI0Go3B\nq0aDXHmhC9kWL4IhIosydNWoBGj3TfVE5sCti4iswsNZhhCFA0IUDvBwlmlvqjdE3031RObEACQi\nmzHlpnoic+MhUCKyGVNuqicyNwYgEdmckJvqicyNh0CJiEiUGIBERCRKNg3A5ORkREREwN/fH4GB\ngZg2bRrOnTvX5npFRUWIjo5Gz549ERISgrVr11qhWiIi6kxsGoB5eXmYO3cuDh8+jL1796JLly6Y\nOHEiampqDK5TV1eHSZMmwdfXF7m5uVi9ejU+/PBDfPTRR1asnIg6uiqVGmerb+Fs9S1Ucco1UbLp\nRTCZmZk6rzdt2gR/f3+cOHECEyZM0LvOzp07oVKpkJqaCicnJwQHB0OpVGLjxo2Ij4+3RtlE1IHx\nAb3UzK7+n7527Ro0Gg169OhhsE9+fj5Gjx4NJycn7bLIyEiUlpbi0qVL1iiTiDowPqCXmtlVAC5b\ntgyDBg3CiBEjDPYpLy+Ht7e3zjIvLy9tGxGRIVUqdZsP6OUTKMTDbu4DfOONN5Cfn4+DBw9CIpEY\n7GesjYjImLaeMNHcznsSxcEuAnD58uXYs2cP9u3bhz59+hjt6+3t3WpPr6KiQtumj1KpNE+hnRzH\nSTiOlXD2NFY35F5t97l+HcqrFVaopjV7Git7FhQUZJb3sXkALl26FNnZ2di3bx8CAwPb7D9ixAis\nXLkSjY2N2vOAubm58PPzg7+/v951zDVYnZlSqeQ4CcSxEs7exqpKpYavvNroEyj6errAw9nwdQiW\nYm9jJQY2PQeYkJCAL774Aps3b4arqyvKyspQVlaG69eva/skJSUhJiZG+3ry5MmQy+WIi4vDuXPn\nsHfvXmzYsAFxcXG2+ApE1IHwCRTUkk0DMC0tDfX19YiJiUFwcLD2r+U9fWVlZSguLta+dnV1RVZW\nFkpLSxHmmh8ZAAARlklEQVQREYGlS5ciPj4e8+fPt8E3IKKOhk+goGY2PQRaXV3dZp+NGze2WhYS\nEoIDBw5YoiQi6uT4BApqZvNzgEREtsAnUBADkIjIiqpUap0ZaNwZwjbDACQisoK2pmAj67OrmWCI\niDqrtqZgk0r5c2xtHHEiIgsTMgWbRu5qxYoIYAASEVmckCnYajUOVqyIAAYgERGJFAOQiMjCfOVS\nnRvv9bW7SW9ZsSICGIBERBYnZAo2aUOdFSsigAFIRGQVbU3BptEYPkdIlsH7AImIrIBTsNkfBiAR\nkRVxCjb7wQAkIhIpsU/LxgAkIhKZtqZlc+4ijstDxPEtiYhIq61p2cSCAUhEJCJCpmWrVKmtWJHt\nMACJiEREyLRsxto7EwYgERGJEgOQiEhEhEzL5mOkvTMRx7ckIiIAwqZlE8t9igxAIiKRaWtaNrHg\nfYBERCLDadnuYAASEYmU2Kdls+kh0Ly8PEybNg0hISFQKBTIyMgw2r+kpAQKhaLVX05OjpUqJiKi\nzsKme4A3btxAaGgopk+fjtdeew0SiUTQert370ZoaKj2dY8ePSxVIhERdVI2DcDx48dj/PjxAID5\n8+cLXk+hUMDLy8tSZRERkQh0yKtAZ82ahaCgIDz++OPIzs62dTlERNQBdaiLYFxcXPD2229j1KhR\nkMlkOHDgAGbPno3U1FRMmTLF1uUREVEHIqmpqWmydREA0Lt3b6xbtw7Tp083ab3Fixfju+++Q15e\nnoUqIyKizqhDHgJtaejQobh48aKtyyAiog6mwwfg6dOn4evra+syiIiog7HpOcDr16/j119/BQBo\nNBpcvnwZp06dgru7O3r37o2kpCQUFhZqL3TJyMiAo6MjwsLCIJVKcejQIaSlpSEpKcmWX4OIiDog\nmwZgYWEhnnnmGQCARCLB6tWrsXr1asyYMQMpKSkoKytDcXGxtr9EIsH69etx+fJlyGQyBAYGIiUl\nBbGxsTb6BkRE1FHZzUUwRERE1tQhzwFyCjXhkpOTERERAX9/fwQGBmLatGk4d+5cm+sVFRUhOjoa\nPXv2REhICNauXWuFam2rPWMl1m1ry5YteOihh+Dv7w9/f39ERUXh8OHDRtcR4zbVzNTxEut2dbfk\n5GQoFAosXrzYaL/2blsd6j7AZpxCTbi8vDzMnTsXw4YNg0ajwTvvvIOJEyfixIkTBr9/XV0dJk2a\nhPDwcOTm5uL8+fOIj49H165dER8fb+VvYD3tGatmYtu2evXqhbfeegv9+/eHRqNBRkYGZs6ciZyc\nHISFhbXqL9Ztqpmp49VMbNtVSwUFBUhPT8fAgQON/sbfy7bVIQOQU6gJl5mZqfN606ZN8Pf3x4kT\nJzBhwgS96+zcuRMqlQqpqalwcnJCcHAwlEolNm7c2Kl/rNozVs3Etm1FR0frvE5MTERaWhoKCwv1\n/qCLdZtqZup4NRPbdtWstrYW8+bNQ0pKCtasWWO0771sWx3yEGh7cQo14Nq1a9BoNEb/JZmfn4/R\no0fDyclJuywyMhKlpaW4dOmSNcq0C0LGqpmYty21Wo3MzEw0NjZizJgxevtwm/oPIePVTKzb1cKF\nCzFx4kSEh4ejqcn4ZSr3sm11yD1AU3EKtf9YtmwZBg0ahBEjRhjsU15ejt69e+ssa/5XaHl5Ofz9\n/S1ao70QMlZi3raKiooQFRWFxsZGyOVyfPrppwgKCtLbl9uUaeMl5u0qPT0dxcXF+PjjjwGgzVNc\n97JtiSIA3d3ddQ6VDhkyBNXV1diwYUOn35haeuONN5Cfn4+DBw8a3aiEnlPtzISOlZi3rfvvvx95\neXmora1FdnY25syZg3379mHo0KGt+nKbMm28xLpdKZVKrFq1CocOHYJMdudBvU1NTUb3Au9l2xLV\nIdCWxDaF2vLly5GVlYW9e/eiT58+Rvt6e3ujvLxcZ1lFRYW2rbMzZaz0Ecu25eDggICAAAwePBh/\n/vOf8eCDD2LLli16+4p9mwJMGy99xLBd5efno7KyEqNGjYKnpyc8PT3x3XffIS0tDV5eXrh161ar\nde5l2xJtAIppCrWlS5dqf9ADAwPb7D9ixAh8//33aGxs1C7Lzc2Fn59fpz9UZepY6SOmbasltVoN\njUajt03M25QhxsZLHzFsV0899RS+//57fPvtt/j2229x/PhxDB06FJMnT8bx48fh4ODQap172bY6\nZABev34dp06dwqlTp3SmUPu///s/AEBSUhJiYmK0/TMyMrBr1y6cP38eSqUSH374IdLS0jBv3jxb\nfQWrSUhIwBdffIHNmzfD1dUVZWVlKCsrw/Xr17V97h6vyZMnQy6XIy4uDufOncPevXuxYcMGxMXF\n2eIrWE17xkqs29bKlSvx/fffo6SkBEVFRUhKSkJeXh6mTp0KgNvU3UwdL7FuV25ubggODtb+DRgw\nAHK5XLscMO+21SHPAXIKNeHS0tIgkUh0NhjgzgUeS5cuBYBW4+Xq6oqsrCwkJCQgIiICCoUC8fHx\nJt1y0hG1Z6zEum2Vl5dj3rx5KC8vh6urK0JDQ5GZmYmIiAgA3KbuZup4iXW70kcikeic5zPntsWp\n0IiISJQ65CFQIiKie8UAJCIiUWIAEhGRKDEAiYhIlBiAREQkSgxAIiISJQYgERGJEgOQyIjVq1dD\noVDYugy70NTUhIceegjvvPOO1T5zxYoVeOyxx6z2eSQuDEASjW3btkGhUGj/fH19MWDAADz33HPY\ntGkT6uvrW61z9ywUQjQ0NGD16tX49ttvzVW6Xdi1axeKi4vx2muv6W0/duwYFAoF/P39deZlbCks\nLAyTJ0/W21ZSUgKFQoG//vWv2mXx8fE4c+YMDh48eO9fgOguDEASneXLl2Pz5s1ITk7GK6+8ol02\nZswYFBUV6fRdvHgxrly5YtL737hxA2vXrkVeXp7ZarYHH3zwAWJiYgzuEe/YsQNdu3ZFfX09Dh06\npLePkH9QtGz38fHBE088gQ8++KD9hRMZwAAk0Rk3bhxiY2MxY8YMLFy4EJmZmdizZw+uXr2K6dOn\nQ6VSafvKZDI4Ojq263PaepJ1R/LTTz/hzJkzmDRpkt52lUqFffv2Ye7cuejXrx+2b9+ut197xmTS\npEn44YcfOv2jgMj6GIBEAMaOHYvFixfj8uXL2LFjh3a5vnOAJ0+eRGxsLAIDA+Hr64uwsDC88sor\nUKlUKCkp0T5G6d1339Uebm2emf7SpUtISEjAiBEj4Ofnhz59+mDq1Kk4e/aszmccP34cCoUCmZmZ\neO+99xASEgJfX1/ExMTgX//6V6v6//d//xdTp05F37594efnhzFjxuD999/X6XPhwgW8+OKL6Nev\nH3x9fTF27FhkZ2cLGp/9+/dDJpNh7NixetsPHjyIa9euYerUqZg6dSqOHDmCmpoaQe/dlkceeQQS\niQT79+83y/sRNWMAEv1b86NpcnNzdZa3PCR39epVTJo0CSUlJViwYAHWrVuHGTNm4Oeff0ZDQwO8\nvLyQnJwMAHj66aexefNmbN68GbNnzwZwJ6i+++47TJw4EWvWrEFcXBxOnTqFJ598EmVlZa1q2rBh\nA/bv348//vGP+NOf/oQff/wRc+fO1elz7NgxPP744zh79ixeeeUVrF69GhERETqHIc+fP49x48bh\n559/xsKFC/Hf//3fcHd3x4svvqgT+Ibk5+cjODgYTk5Oetu3b9+O0NBQDBgwALGxsbh9+zaysrLa\nfF8h3Nzc0LdvX/zwww9meT+iZh3ycUhEluDn5wcXF5dWe1gtD9udOHECNTU1yMrKwpAhQ7TLly9f\nrv3fzzzzDF5//XUMHDiw1eNrJkyY0OpxS1OnTsWoUaPw+eefIyEhQaft5s2byMnJQZcud/5T7dGj\nB5YtW4Zz585hwIAB0Gg0WLBgAby8vHD8+HH06NFD73dbtmwZevXqhdzcXG2IzZkzB88++yySkpIw\nZcoUo2OjVCoxdOhQvW1VVVXIycnBihUrAAABAQEYOXIkduzYgZdeesno+wrVp08fnD9/3izvRdSM\ne4BELXTv3l3v1aDN3NzcAACHDh3C7du3TX5/Z2dn7f++ceMGqqqq4OLigv79++PkyZOt+k+bNk0b\nfgAwatQoAHeumATunJsrKSnBq6++ajD8qqurcezYMcTExKC+vh6VlZXav8jISPz+++/49ddfjdZd\nVVVl8P13794NtVqtE/bTpk3DiRMncOnSJaPvK1SPHj1QWVlplvciasYAJGqhvr4e3bt3N9geHh6O\nmJgYvPvuu+jXrx+mTZuGv/3tb7hx44ag91epVPjzn/+M4OBg9OrVC/3790dgYCDOnj2La9euterf\nu3dvndfNIdR8fq15bzUkJMTgZ168eBFNTU1Ys2YNAgMDdf5WrFgBiUSCiooKo3VLJBKDF7Ds2LED\ngwYNQkNDAy5evIiLFy9i0KBBkEql2Llzp9H3NfRZd2tqaoJUyp8rMi8eAiX6t99++w3Xrl1Dv379\njPb77LPPUFhYiEOHDuHo0aNYsGABkpOTceTIEXh6ehpdd8mSJdi2bRteeeUVjBw5Em5ubpBIJFi+\nfDk0Gk2r/jKZTO/7mHI1ZfP7zp8/H+PHj9fbJzg42Oh7uLu7672opbi4GAUFBQCAYcOGtWrfuXMn\nFi1apH3t7OyMhoYGvZ/RvFzfecaamhp4eHgYrZHIVAxAon9rvnQ/MjKyzb7Dhg3DsGHD8MYbb+DI\nkSOIjY1Feno6Fi1aZPQ+tz179mD69OmtZlOprq5u1w983759AQBFRUUG6w4ICABwJ0wfeeQRkz8D\nAP7rv/4LxcXFrZZv374dMpkMmzdv1jlUCwCnTp3Ce++9h59++gmDBw8GANx3332trnht9ssvv2j7\n3K2kpMToXi5Re/CYAhHuXEm5bt06BAQEGL0gpKamptXe16BBgwAAdXV1AAC5XA7gTqjdrUuXLq32\n9Hbt2mXyzfbNhgwZgoCAAGzatKnVHlpznV5eXnj44YeRnp6O0tLSVu9x9erVNj9n5MiROH/+fKsZ\nXnbu3IlRo0bh2WefxTPPPKPzt2DBAjg6OupcZRoVFYWysrJWt1/cvn0bn3zyCbp164bw8HCdttra\nWhQXF2PkyJFt1klkCu4BkugcOXIEFy5cwO3bt1FRUYF//OMfOHr0KPz9/fHFF18YvfE9IyMDH3/8\nMZ5++mkEBASgoaEB27ZtQ5cuXbRXd8rlcgwYMAC7d+9GYGAgFAoFAgIC8MADD+CJJ57Al19+CRcX\nFwwYMACnT59GVlYWAgIC2nWTuEQiwfvvv48pU6YgPDwcs2bNgq+vLy5evIj8/HztrRDJycmYMGEC\nHnroIbzwwgsICAhARUUFfvzxR/zyyy8oLCw0+jnR0dFYvXo1jh07hqioKADAP//5T/z6668Gr/R0\ncXHBQw89hN27d2PVqlWQSqV48cUXsXXrVsyZMwfffPMNBg8ejLq6OuzduxcnT57E6tWr4erqqvM+\nR48eRVNTE6Kjo00eHyJjGIAkGs2HJtesWQMAcHR0hEKhwMCBA7FmzRrMnDkT3bp1a7VOy0Oa4eHh\nOHnyJLKyslBeXg4XFxcMHjwY69ev1zkH9uGHH2Lp0qVITExEY2MjZsyYgQceeABr1qyBg4MDsrKy\nsHXrVgwdOhS7d+9GYmJiq0OnQucgffTRR7F//368++672LhxI9RqNQICArT3NQJAYGAgcnNz8e67\n7+LLL79EZWUlvLy8EBoaisTExDY/IzQ0FIMHD0ZWVpY2ALdv3w6JRIInnnjC4HrR0dE4evQo/vGP\nf+DRRx+Fs7MzDhw4gHXr1mH//v348ssv4eTkhLCwMPztb3/D008/3eo99uzZg1GjRqF///6CxoNI\nKElNTU3nma+JiCwmMzMTCxYswOnTp632hIwrV65gyJAh+OSTT7gHSGbHc4BEJMhzzz2HgIAApKam\nWu0zU1JSEBoayvAji+AeIBERiRL3AImISJQYgEREJEoMQCIiEiUGIBERiRIDkIiIRIkBSEREosQA\nJCIiUWIAEhGRKP0/ArDycTAefTAAAAAASUVORK5CYII=\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x103794b70>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"data.plot('Distance (AU)', 'resonance', s=80, kind='scatter',\n",
" color=color);"
]
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3.4",
"language": "",
"name": "python3.4"
},
"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.4.3"
}
},
"nbformat": 4,
"nbformat_minor": 0
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment