{ "metadata": { "name": "", "signature": "sha256:92507752e7bb847826e05f21eb9f820b3add651e9aeed3bb7600ac570198f0fc" }, "nbformat": 3, "nbformat_minor": 0, "worksheets": [ { "cells": [ { "cell_type": "code", "collapsed": false, "input": [ "%matplotlib inline\n", "\n", "import numpy as np\n", "import matplotlib.pyplot as plt\n", "\n", "from isochrones.dartmouth import Dartmouth_Isochrone\n", "dar = Dartmouth_Isochrone()\n", "\n", "from isochrones.starmodel import StarModel" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 2 }, { "cell_type": "markdown", "metadata": {}, "source": [ "Here, I create a stellar model object with observations of $T_{eff}$, $\\log g$, and $[Fe/H]$ corresponding to properties in the Huber stellar properties catalog paper, which agree closely with the Borucki (2011) values, if not exactly. " ] }, { "cell_type": "code", "collapsed": false, "input": [ "mod = StarModel(dar, Teff=(5642,50), logg=(4.44,0.06), feh=(-0.27,0.08))" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 3 }, { "cell_type": "markdown", "metadata": {}, "source": [ " I then run an MCMC chain using the Dartmouth isochrone grids, which can predict $T_{eff}$ and $\\log g$ (and $[Fe/H]$) as a function of mass, age, and $[Fe/H]$, where the likelihood surface is determined by the observations above. This gives chains of mass, age, and $[Fe/H]$ that can be used to read off distributions of other stellar properties, such as radius." ] }, { "cell_type": "code", "collapsed": false, "input": [ "mod.fit_mcmc()" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 4 }, { "cell_type": "markdown", "metadata": {}, "source": [ "We can look at the result of these chains, which support the idea of an older, smaller star than the literature values (plotted with red vertical lines) seem to indicate. This discrepancy is pretty severe, and I'm not sure how to explain it." ] }, { "cell_type": "code", "collapsed": false, "input": [ "mod.plot_samples('radius')\n", "plt.axvline(0.979, color='r', ls='--')\n", "mod.plot_samples('mass')\n", "plt.axvline(0.970, color='r', ls='--')\n", "mod.plot_samples('age')" ], "language": "python", "metadata": {}, "outputs": [ { "metadata": {}, "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAYsAAAEXCAYAAABcRGizAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XmcXFWZ//HPN90J2QgJZBEJCYPsizGAMCJiTxgngoiO\nC4iI4jg44KCMw4jDmgRUBGXQQYiDoAwBFX5icAQBZUiDogOEVTESAgkJCQFiIBsEO83z++Pe7lRX\nqvre7q6tO9/361WvusupOk9v9fQ959xzFBGYmZl1Z1C9AzAzs8bnZGFmZpmcLMzMLJOThZmZZXKy\nMDOzTE4WZmaWycnCzMwy1TRZSLpe0vOS1kp6RtI5ZcqdJKld0rqCx+G1jNXMzDZrrnF9FwH/GBEb\nJe0J3CPpoYi4o0TZ+yLCCcLMrAHUNFlExBNFhzYBL5YpriqHY2ZmOdW8z0LSlZI2AE8AX4mIh0sU\nC2CqpJckPSnpXElNtY3UzMw6qB5zQ0kS8G7gJ8BREfFA0fm/At6IiGcl7QfcCMyJiK/XPFgzM6tP\nsuisXJoNbIyIL2aUOw74UkQcVHDMMyCamfVCRPS4mb/eQ2cHAxtylnUfhplZndQsWUgaJ+ljkkZI\napI0Hfgo8LMSZY+UNCHd3gs4F7il1PtGRGUeJB0lFXu/9DFjxoyKv2c1Ho7TcTbyw3FW5tEXtbyy\nCOAU4Dngz8CFwIkR8aCkSem9FBPTstOAxyStB24Dbga+VsNYzcysQM2GzkbEKqClzLmlwLYF+18C\nvlSbyMzMLEu9+ywGvJaWlnqHkIvjrCzHWVmOs/7qOhqqLzpGQ1UsfqX95/30+2FmlkXp51z0w9FQ\nZmbWD9R6bqjG5SsKM7OyfGVhZmaZnCzMzCyTk4WZmWVysjAzs0xOFmZmlsnJooO0+V4LMzPrwsnC\nzMwyOVmYmVkmJwszM8vkZGFmZpm2quk+li+Ht72t9LmX0udx4+D556F5q/rOmJl1b6v6SHzjDVi1\nqvQ5kc4NVea8mdnWzM1QZmaWaatNFjvuCC++uPkxaKv9TpiZZduqmqEKNTUl/RMdfD+emVl5/n/a\nzMwyOVmYmVkmJ4vUpnYRuC3KzLo3f/58Tj/9dObMmcMpp5zC008/3auy//AP/8DgwYMZNmwYBx98\nMA8//HAtwu+1mvVZSLoeOAIYQTJA9ZqI+GqZsl8EzgSGAz8BTo2Iv9QqVjPbes2ePZtTTz215LnX\nX3+dj3zkI9x///1MmDCBvffem+OPP54HHnigx2UnT57M8uXLiQgmTJhQ1a+pEmp5ZXER8FcRMQo4\nEvi8pPcWF5I0HfgyMA2YDOwKzKphnGa2FXvxxRfLnrv33nsZOXJk54f7gQceyIIFC1iyZEmvyo4f\nP75fJAqoYbKIiCciYmPBoU1AqZ/Kp4CrI2JBRLwCXACcVIMQzcy6tWTJEnbYYYfOfUmMGTOGJ554\nosdlX3vtNa666ip+/OMfc/LJJ7NgwYLqfwF9UNOhs5KuJEkG2wCnRUSpRrp9gLkF+48DEySNiYiX\naxCmmVXJk08+ydlnn83kyZORxKpVq/jGN77B+PHjM1/7zDPPMGvWLIYNG0ZTUxMbNmzg4osv3uI/\n877UkWXVqlUMHz68y7GhQ4eybt26Hpfdf//9OfbYYxk8eDDjx4/ngx/8IH/6059Qg47jr2myiIjP\nSfpn4N3ATyQ9HBHFjX0jgTUF+2vT520BJwuzfmrNmjUcccQRXHLJJXz84x8H4KKLLuKII47g4Ycf\nZvDgwWVfu3jxYg4++GCuvvpqPvjBDwJwww03MH36dB566CGampp6VceLL77It771LSKi89hvfvMb\nNm7c3AgycuRIzjnnHABGjx7dpSzA+vXrGTt27BYxZ5U97rjjaE4nodt111156qmn+P3vf89b3/rW\nst+HuoqIujyA2cBlJY4/CnykYH8s8AYwpqhcADFjxozOx7x586I7S5dGQPKYOLHruaamzefa2rp9\nGzPrhbPPPjsmTJgQ7e3tncdWr14dgwcPjtmzZ3f72g984AOxww47dDn26quvxtChQ+Paa6+tSB0d\nZs6cWfbc3XffHVOmTOncb2tri2222SYWLlzYo7K/+93vYuTIkbFx48aIiHj88cdDUixYsCBXjD0x\nb968zs/Ijs/N6M1ndm9eVIkHcDXwlRLHbyg8TjKC6vkS5Tq+6NycLMzqZ/fdd49jjjlmi+P7779/\nTJs2rezrXn/99Whubo5DDjlki3N77rlnHH300X2uo1B3yaKtrS3e/OY3x9KlSyMi4q677ooDDjig\n8/xdd90Vjz76aGbZ5557LmbNmtX5utmzZ8dhhx2WK76+6EuyqEkzlKRx6Yf+z4GNwN8CH02fi10H\nXCvpBmAlcB7wg1rEaWbVsW7dOhYtWsT06dO3OLfjjjty//33l33t6tWraW9vZ+jQoVucGzVqVOdQ\n1L7UkVdzczNz5szhq1/9Koceeijz5s3jxhtv7Dx/xRVXcMABBzBlypRuy+60004ccMABXHrppbS3\nt/PUU08xd+7cctU2hFr1WQRwCknTk4CFwIkR8aCkScATwN4R8VxE3CnpEmAeMIzkPosZNYrTzKrg\n2WefBZIP92IjRoxg7dq1tLW1ley3GDduHMOHD+/Sj9BhxYoVrFq1ijfeeKNPdRQaNmxYt+enTZvG\ntGnTAPjkJz/Z5dxPf/rT3GWPPvpojj766G7raiQ1GTobEasioiUixkTE6Ig4OCL+Jz23NCK2jYjn\nCspfFhFviojtIuIzEdFWizjNrDrWrk3GqQwZMmSLcyNGjADglVdeKfnapqYmjjvuuC3ulF65ciUr\nVqwAkquPvtRR6Mwzz8wsszXydB9mVnUdo5VKDQtta0v+F2xvby/7+osvvpjtt9+eq666CoBNmzZx\nxRVXMHXq1M7372sd1j0ni5TnhjKrnnGF6wEU2bBhA5LYdttty5YZO3Ys999/PytXruSMM85g1qxZ\nnHzyyUhi6NChjBkzps91WPe22vUszKx2JkyYgCRefnnLW6U2bNjA6NGjO5uKyhk9ejTnn39+l2Mv\nvPAChx56aMXqsPJ8ZWFmVTdixAimTp3K0qVLtzi3aNEi3va2t/X4PV966SWWL1/ORz/60arVAZWb\nZbbQGWecwX333dereOqmN+NtG+FBhe+z6Djh+yzMquP888+PnXbaqcuxRYsWhaS48soruxxfuHBh\nvPbaa537l112WYwfPz6WL1/eeewb3/hG7LzzzvH666/3qo48Nm7cGJMnT46VK1dGRMSDDz4Yb3/7\n2/tUtrW1NSZMmBD33HNPj+PpK/pwn4WvLMysJk499VQ2bNjA9ddf33ns8ssvZ9999+Xkk0/uPHbv\nvfey5557dk7rAUkz0ogRIzo7sR955BG++c1vctNNN3UZ/ZS3jrwqPcvsmjVrePzxx9l77717HEu9\nuc/CzGriTW96E62trZxzzjk8/PDDrFu3jpdffpk77rijc46kjnLjx49nt9126zx2xhlnsGbNGs46\n6yzWrVvHmjVr+NnPfsYhhxzSqzry6m7m2F122aXHZb/3ve/xhS98YYv7MfoDJ4tUc1PgUXVm1TVl\nyhRuvfXWbsvssccerFy5ssuxoUOHcskll1SsjrwqOcvsz3/+c4466qiS94H0B04WZrbV+c53vtNt\nR/VBBx3ECSecwHbbbdfRR9qpN7PMrlixgldeeYV99tmn81xx2UbnZGFmW53TTjstV7m9996780ZA\nSG4GXL16NZMnT96i7F577cV//dd/lSz7y1/+kpUrV3LxxRcDsHDhQn74wx+yYcMGjjrqqD5+NbXh\nZGFmVsa73vUuXnrpJZYtW8bOO+/MPffcw7777svuu+8OwP/+7/8yduxYpkyZ0m3ZjvIdvvvd73LC\nCSdw+OGH1+PL6hUnCzOzMio1y2yH5cuX8+1vf5uVK1dy6aWXsm7dOt73vvfV+svqFfW3drMOkpKb\nLXoQ/7JlMGlSsj1xYrLfobmZzg7utrZk38xsIOmYNysiejy3ke+zSHluKDOz8pwszMwsk5OFmZll\ncrIwM7NMThZmZpbJycLMzDI5WaSamwLRP4cRm5lVm5OFmZllcrIwM7NMNUsWkoZIukbSEklrJT0i\n6b1lyp4kqV3SuoJH/5lExcxsgMmVLCT9W5nj/9qDupqBpcDhETEKOBe4SdKW0zcm7ouIbQse9/ag\nLjMzq6C8VxYzyhw/L29FEfFqRMyKiKXp/m3AYuCAMi/x3BtmZg2i2+nyJE0j+dBuSrcLvQVY29uK\nJU0A9gCeKHE6gKmSXgJWA3OAiyKiamvZbWpPcpNHRJmZbSlrbtXvk3xwbwNcU3A8gBeAz/emUkmD\ngRuAayNiYYki9wL7RsSzkvYDbgQ2AV/vTX1mZtY33SaLiNgFQNKciDixEhVKGkRypbARKLlcVUQs\nLtj+g6QLgC9RIlnMnDmzc7ulpYWWlpZKhGlmNiC0trbS2tra5/fp8XoW6Yd9p4h4owevFcnVyiTg\nqIh4PefrjgPOjIgDC45VdD0LtLkZyutZmNlAVPX1LCQdKOl3kl4laQ7qeLT1sL7ZwF7AMd0lCklH\npn0aSNqLZOTULT2sy8zMKiTv/8//DfwP8Bng1d5UlA6R/SxJ89PKjgyXHruPpKN774h4DpgG/EDS\nSJK+kTnA13pTr5mZ9V2uZihJa4HtooHWYPWyqmZmPVOLZVXnAtN7+uZmZjYw5P3/eRgwV9KvSZqF\nOkREfLLyYZmZWSPJmyz+mD6KNUyzlJmZVU+uZBERM6sch5mZNbBcyaLEVB+dIuLuyoVjZmaNKG8z\nVMe0Hx3GkUwBsgzYtdJB1YPnhjIzKy9vM9QuhfuSmkhulFtfhZjMzKzB9Grxo3T2168BZ1Y2HDMz\na0R9WSnvPUDVpgw3M7PGkbeDe1nRoeHAUOBzFY/IzMwaTt4O7uLpyTcACyNiTYXjMTOzBpS3g7sV\nOqcnnwC80JOpyfuD5qbonBvKzMy6yjtF+ShJ15HMGLsc2CjpOknbVTU6MzNrCHk7uC8HRgD7kfRX\ndDxfXqW4zMysgeTts3gvsGtEbEj3F0o6CXimKlGZmVlDyXtl8RrJXduFxpI0S5mZ2QCX98riauBX\nki4FngV2Ab4IfK9KcZmZWQPJmyy+BqwATgB2TLcvJpkzakDw3FBmZuXlHTr7BkliGDDJwczM8ss7\ndPZySYcWHTtU0reqE5aZmTWSvB3cxwMPFR17mKRZyszMBri8yeKNEmUHAapsOGZm1ojyJovfAF9J\np/voWM9iFvDrPC+WNETSNZKWSFor6RFJ7+2m/BclPS9pTfq6ITnjNDOzKsibLE4H/hZYKelBktFQ\n7wG+kPP1zcBS4PCIGEWycNJNkiYXF5Q0HfgyMA2YTLIS36yc9fRac1N4JJSZWRmKyPcBmV5NHAzs\nTLKc6v19mUxQ0mPAzIiYW3T8h8AzEXFuuv83wA8jYseicgGQN36AZctg0qRke+LEZL9DczOdEwm2\ntSX7ZmYDiZT0HEREj7sQcn8kpqvj/S599ImkCcAewBMlTu8DFCaQx4EJksZExMt9rdvMzHquLyvl\n9YqkwcANwLURsbBEkZFA4ToZa9Pnbasdm5mZlVbTxpa0g3wOyZxSp5Upth4YVbDfMQ36ulKFZ86c\n2bnd0tJCS0tLX8M0MxswWltbaW1t7fP75O6z6HNFSWPZ94FJwFER8XqZcjcAiwv6LI4ArnefhZlZ\n3/Slz6JsM5SkQXkePahrNrAXcEy5RJG6DviMpL0ljQHOA37Qg3p6ZVO7CN82YmZWUncf9puKHm1l\njmVKh8h+FphCMvx2Xfo4XtKkdHsiQETcCVwCzAOWAE8DM3rxtZmZWYV019iya8H2+4CPkMw+u5Sk\nKenfgZvzVBIRz9J9YurSeR0RlwGX5XlvMzOrvrLJIiKWdGxL+lfgoIKhq09Kmg/MB66saoRmZlZ3\nefscRpGsuV1oOJtHKpmZ2QCWd8zPfwN3SbqM5O7tSSRTfVxXrcDMzKxx5E0WZwKLgI+RrJT3PHA5\nA2hZ1eam6Bw6a2ZmXfVkpbzvpg8zM9vK5F0pb5Ckz0q6W9Lv02OHSzq2uuGZmVkjyNvBPQv4DEmz\nU3oPNMtJhs+amdkAlzdZfBo4OiJ+RLJqHsBiut6LYWZmA1TeZDGIZIK/QiMoM7mfmZkNLHmTxe3A\nf0gaCp2zx14I/LxagdWa54YyMysvb7L4V+BNwCskN+itB3bBfRZmZluFvENn1wB/n65wNxlYFhHP\nVzUyMzNrGHmHzv5A0vSIeCEiHuhIFJI8L5SZ2VYgbzPUJ4DvS/pS0fETKxyPmZk1oLzJYiNwCPAx\nSddL2qaKMZmZWYPJvdJdRDwHHJa+5tcdixUNFM1NgajNErNmZv1NT5ZFJSJei4iPAz8F7geGVCUq\nMzNrKHlnnb2wcCcivi7pMcBzQ5mZbQXyDp29pMSx20lu1jMzswGubLKQdGdETE+3f12mWETE4VWJ\nzMzMGkZ3VxaFq+BdU6aMe4TNzLYCZZNFRNxQsH1tTaKpo03tybxQHhFlZral7pqhPkOOK4eI+H6e\niiSdBpwE7Af8KCI+XabcSSRXMq8WHH5fRNybpx4zM6u87pqhTiRfM1OuZEGyWNKFwHRgWEbZ+9wX\nYmbWOLprhmqpZEURMRdA0kFA1g19nivczKyB9OimPAAlBnU8elFnViIIYKqklyQ9KelcSU29qMfM\nzCok76yzO0maK2k1sKng0daLOrOatu4F9o2IccCHgeOB4gkMzcyshvLewf1d4DVgGnAP8G5gBr27\nKa/bK4uIWFyw/QdJF5Aki6+XKj9z5szO7ZaWFlpaWnoRUjI3VHt7r15qZtawWltbaW1t7fP7KCK7\nDzu9opgUEeslrYmI7SRtD/w2IvbqUYXShcDEcqOhSpQ/DjgzIg4sOh4AeeLvsGwZTJqUbE+cmOx3\naG6mM1m0tSX7ZmYDiZT8rx4RPe4Xztvn0NHsBPCypPHABmCnvBVJakrX8G4GmiRtU6ovQtKR6Yp8\nSNoLOBe4JW89ZmZWeXmTxQPAken2ncCNwFxgfg/qOo/k3okvkyym9BpwjqRJktYVTHk+DXhM0nrg\nNuBm4Gs9qMfMzCosbzPUaGBQRKyWNBw4AxgJfKtea3G7GcrMrGf60gyVd9bZVwq2X6VoynIzMxvY\nciULSYNJhrBOJbmi6BAR8dlqBFZrnhvKzKy8vI0tc4D9SYbKrkyPCc86a2a2VcibLN5LMnR2bTWD\nMTOzxpR3NNQCYPtqBmJmZo0r75XFJ4DvS7odeCE9JpI+i+vKv8zMzAaCvMniU8A7gVEk90cUcrIw\nMxvg8iaLfwGmRsQfqxlMPXluKDOz8vL2WbwALK1mIGZm1rjyXln8BzBH0iVs7rMAICKeqXhUZmbW\nUPImiyvS5w8UHQ/ACxOZmQ1wmc1QSiYT2R3YJiIGFT2cKMzMtgJ5+yweB9z9a2a2lcpMFpFM6/oo\nsGf1w6mfTe0iMpcHNzPbOuXts5gH3C7pWmAZSV9Fx015369SbGZm1iDyJovDgCUka28Xc7IwMxvg\n8q5n0VLlOMzMrIHlXg9O0hjgGODNwHLg1ohYXa3AzMysceQaDSXpHcDTwD8BbwVOARZJOrSKsZmZ\nWYPIe2XxbeBzEfHjjgOSjkuPv70agdWa54YyMysv730WewA3FR27meRmPTMzG+DyJounSNbgLvRR\nYFFlwzEzs0aUN1mcDnxH0v9JuknS/cCV6fFcJJ0mab6kjZJ+kFH2i5Kel7RG0jWShuStx8zMKi9X\nsoiI3wJvIZlQ8CHgcmC3iLivB3UtBy4k474MSdOBLwPTgMnArsCsHtRjZmYVlnvobDpMdk5vK4qI\nuQCSDgImdlP0U8DVEbEgLX8B8EPgrN7WbWZmfdNtspA0L+P1ERHTelhn1gRM+wBzC/YfByZIGhMR\nL/ewrtw2tSdhiahWFWZm/VbWlcUNJY4FsBNJf8XwXtSZ9Wk8ElhTsL82fd4WqFqyMDOz8rpNFhFx\ndeG+pLHAvwMnAzcCF/Sizqwri/XAqIL97dLndaUKz5w5s3O7paWFlpaWXoRkZjYwtba20tra2uf3\nUTIDeUYhaTvg34DPA7cCMyLi6V5VKF0ITIyIT5c5fwOwOCLOTfePAK6PiB2LygVAnvg7LFsGkyYl\n2xMnJvsFb5g8EbS1QXPu3hwzs/5B6edcRPR4PYZuR0NJGi7pLOAZkr6Ed0bEJ3qTKCQ1SRpKcjXT\nJGkbSaVW2rsO+IykvdP5qM4Duh1qa2Zm1ZX1//NikoRyCTCfpKN5QmGBiLg7Z13nAecX7H8CmJmu\nkfEEsHdEPBcRd0q6hGQNjWHAT4AZOeswM7Mq6LYZStKSdLNsoYj4qwrHlEulm6Gam+mcG8rNUGY2\nEPWlGSqrg3uX3oVkZmYDSd7pPszMbCvmZGFmZpmcLMzMLJOThZmZZXKySG1qF5F5c7mZ2dbJycLM\nzDI5WZiZWSYnCzMzy+RkYWZmmZwszMwsk5NFqrkpvEqemVkZThZmZpbJycLMzDI5WZiZWSYnCzMz\ny+RkYWZmmZwsUp4bysysPCcLMzPL5GRhZmaZnCzMzCyTk4WZmWWqWbKQtL2kuZLWS1oi6fgy5U6S\n1C5pXcHj8FrFaWZmW2quYV1XABuB8cBU4DZJj0XEH0uUvS8iapogmpuC9vZa1mhm1n/U5MpC0gjg\nQ8B5EfFqRNwH/Aw4sdxLahGXmZnlU6tmqD2ATRGxqODYY8C+JcoGMFXSS5KelHSupKaaRGlmZiXV\nqhlqJLC26Ng6YNsSZe8F9o2IZyXtB9wIbAK+Xt0QzcysnFoli/XAqKJj25EkjC4iYnHB9h8kXQB8\niTLJYubMmZ3bLS0ttLS09D1aM7MBorW1ldbW1j6/jyKqv+BP2mexmuSKYVF6bA6wLCLOznjtccCZ\nEXFg0fEA6En8y5bBpEnJ9sSJyX6H5mY6O7jb2pJ9M7OBREq6gyOix/3CNemziIgNwE+BCyQNl3QY\n8H5gTnFZSUdKmpBu7wWcC9xS7Rg9N5SZWXm1vCnvc8Aw4EXgeuCUiFggaVJ6L8XEtNw04DFJ64Hb\ngJuBr9UwTjMzK1KTZqhqqHQzFOnlmQg3Q5nZgNTwzVBmZta/OVmYmVkmJwszM8vkZJFqbgpE/+y/\nMTOrNicLMzPL5GRhZmaZnCzMzCyTk4WZmWVysjAzs0xOFinPDWVmVp6ThZmZZXKyMDOzTE4WZmaW\nycnCzMwyOVmYmVkmJ4uU54YyMyvPycLMzDI5WZiZWSYnCzMzy+RkYWZmmZrrHUAjmj8fmkt8Z/bb\nD4YOrX08Zmb1VrNkIWl74BrgPcAq4KyI+FGZsl8EzgSGAz8BTo2Iv1Qzvk3tybxQInjHO0qXWbgQ\ndt+9mlGYmTWmWjZDXQFsBMYDJwCzJe1TXEjSdODLwDRgMrArMKuGcZqZWZGaJAtJI4APAedFxKsR\ncR/wM+DEEsU/BVwdEQsi4hXgAuCkWsQJcOCBXR9DhvTt/VpbWysSV7U5zsqqR5yrVsHtt3f/uOee\n+sfZG46z/mrVDLUHsCkiFhUcewxoKVF2H2Buwf7jwARJYyLi5eqFmJg/v+v+7rvDokWly5aydi1s\n3Lh5/7bbWtlnn5YuZYYPh5Ejs9/noYe6LzN4MBx2WP7YutPa2kpLS0tl3qwPbrkF2tpKn2tvhxtv\nbGXPPVsYM2bL84MHQ1NTZeNZuhROP737MkOGwI03bt6/6CK4/fZWDj+8pexr/umfYOedk+0VK2Dd\nutLlIuCVV2DsWNhhh9J1jxiRbD/8MBx1VPex7rYbPPXU5v1G+bln6U2cf//3sHx592WuuALe/vbe\nx1Ws0t/PY4+FJUu6L3PZZfDOd1asyrJqlSxGAmuLjq0Dti1Tdk3BfsfrtgWqniy6M21a8oEEsHhx\n8jxuXNcyL7205eu++c0tj735zV33V6xIno8+Onm+//7S71Xsr/968/b//R+MHg177LFluTVr4Mkn\nYZ99YNddk2O33po8H3VU0h/z4IPJ/i9+kTz/zd+UrnPevOT5yCO7Hr/99q5fw6JF8Kc/Jb/I25b4\nSS9YAM8+u7l8YUzdueWW8uc63qvjfcrFuNtuoILlSzo+QIcM6Xo1uX59djwAd9yxeXtt+hv761+X\nL//Vr8JOOyXbWR9oWTo+KO67L7vsokXw/vdv3n/yyex/Sqrhl7+Ev/wFjjgCBpVo3/jVr5Lnjp9n\nb+LM87v07ncnMVRKpb+feb6GNWuyy1SCIqo/xYWkqcBvImJEwbF/Aw6PiGOKyj4KfCUifpLujwVe\nBHYovLKQ5Lk5zMx6ISJ6vNJbrTq4FwLNknYrODYF+EOJsk8Abysq90ItmqDMzKy0mlxZAEj6ERDA\nPwIHALcC74iIBUXlpgPXkoyGWknSf/HbiDi7JoGamdkWajl09nPAMJImpeuBUyJigaRJktZJmggQ\nEXcClwDzgCXA08CMGsZpZmZFanZlYWZm/ZfnhjIzs0z9MllI2l7SXEnrJS2RdHy9YyomaYika9L4\n1kp6RNJ76x1XOZJ2l7RR0px6x1KOpI9JWpD+3BdJqtBdJpUjaaKkn0v6s6TnJV0uqcJ3f/Q4ptMk\nzU9/vj8oOneEpD9J2iDpbkmTGi1OSX8t6Vfp9/RFSTdJelOjxVlU5nxJb0iaVuv4CmLo7uc+XNKV\nkl6S9Iqke8q9T4d+mSzIOXVInTUDS0mGB48CzgVukjS5vmGVdQXwADTmcoGS3gN8HfhURIwE3gU8\nU9+oSvpPkrnPdiQZ1fdukv66eloOXAh8v/BgOiz9ZuAcYAwwH7hxi1fXTsk4gdHAd0mm/5lMco9W\nyQ/pGikXJwCS3gJ8BFhRy6BK6C7Oq0i+r3uR/Oz/JevN+t2sswVTh+wbEa8C90nqmDrkrLoGVyCN\nbVbB/m2SFpOMBHu2boGVIOljJDc8/hHYLaN4vcwCZkXEAwAR8Xyd4ylnX+D0dOLLFyTdkR6rm4iY\nCyDpIGBiwakPAX+IiJvT8zOBVZL2iIiFjRJnRNxRWE7SFUBrTYMr0M33s8N3SOa3u7KWcRUrF6ek\nvYD3Azs0fpeBAAAFZElEQVRFRMdtp49kvV9/vLIoN3VIXf8gs0iaQBL7E/WOpZCkUSQfxF8Eenyj\nTi2kzTgHAuMlPSVpWdq804gTxt8JfFzSMEk7AUcCt9c5pg7FP999Sf52gM5/cBYB+9UyqBKyfg8P\np/Q9WrW2RZySPgpsjIhG+ZnDlnEeTPIP6wVpM9Tjkj6U9Sb9MVn0ZOqQhiBpMHADcG09/mPLcCHJ\nxI0raNAmKGACMBj4MHAYSfPOVJKmvUYzk+TDdi2wDHgwIn5W14g2K/75jmDLv6W1JH9j9VT291DS\nW4HzgC/VLpyyusQpaVvgq0DGbGI1V/z9nEjyO/oKSXPpacB/p1ccZfXHZLEeGFV0bDuShNFwJA0C\n5pD0sZxW53C6kPQ24AjgWx2H6hhOd15Lny+PiBci4s/AfwAZ0+bVliSRXFn8P5K1WMYC20u6uK6B\nbVb8823Uv6WSv4fpDBC/AL6Qzlxdb8VxzgTmRMTSbsrUQ3EMrwFtJNMqbYqIe0nua/u77t6kPyaL\nnkwdUlfph8c1wDjgwxHRXueQir0b2AVYKul54Azgw5Lmd/uqGkunenmu3nHkMJakuew7EdEWEatJ\nZiNolKRW/B/mEyR/O0Bnf+BbqH9T6RZXFunAkF8BF0TEDbUPqaTiOKcBX0hHwT0P7EwyqKXeV0HF\ncT6ePhcnkW5bFvpdsoiIDcBPSdrbhqfDJ99P8t97o5lNMtrgmIh4vd7BlHAVyeJSU0iadr4L3AZM\nr2dQZfwA+LykcZLGkPSx/LzOMRVbBTwPnCqpSdJokvVZHuv+ZdWVxjKUZEBLk6Rt0n6gucB+kj6U\nnp8BPFqvptJycaZ9P3eTJOGr6hFboTJxNpNcpe/L5r+nFcBnqVNHdzc/93tIRmqeJalZ0jtJlou4\ns9s3jIh+9yAZ6jWX5DJ6CfCxesdUIsbJwBvAqySX9R2P4+sdWzcxzwCuq3ccZWJrJhne+zLJB/K3\ngCH1jqtEnIcAv07jfAn4MTCuzjHNTH8XCx/np+eOABakv6d3A5MaLM4ZwPnpduHf0doGi/P8EuUW\nA9MaMU6SdYN+m36G/gH4QNb7eboPMzPL1O+aoczMrPacLMzMLJOThZmZZXKyMDOzTE4WZmaWycnC\nzMwyOVmYmVkmJwuzCpG0S7rgzaB0/xeSTqx3XGaV0O/WszDrLyKiUeaEMuszX1mYlZHO92NmOFmY\ndaFkzfQzJT0OrJd0jpL1vtdKekLSBwvKDpL0zXQBmaeB9xW9V6ukz6TbM1WwvnmJJquTJD2d1vOM\npI/X5is2y8f/OZlt6WMkK9z9GTgaOCwiVko6Frhe0lsi4gWSGUXfRzLD6KsksyEXTrYWBfvdLegz\nAvg2cFBEPJWuqrhDhb8msz7xlYVZVwH8Z0Qsj4iNEfGTiFgJEBE3AU+RLEsJcCxwWVr2ZeBrlF/s\nJmsRnDeA/SUNi2SBpz/2/UsxqxwnC7MtLevYkPRJSY9IelnSyyTLUY5NT+9YWJZkjYAei2SNluOA\nU4AVkm6VtGfvQjerDicLsy0FdK7OdhXwz8D2ETGGZO7/jquE54FJBa8r3C62nmSp1Q5v6lJhxC8j\n4u/S438CvteXL8Cs0pwszMobQZI4VgGDJH2a5Mqiw00ky2julK7e9+/dvNejwOGSdpa0HXBWxwlJ\n4yV9IO27aAM2AI22BK9t5ZwszMpI+w0uBX4HrCRJFL8pKPI9kqUoHwPmAzdTpiM7Iu4CbiRZ//hB\nkiVhO8oOIlkmdjlJp/q7gFMr+9WY9Y1XyjMzs0y+sjAzs0xOFmZmlsnJwszMMjlZmJlZJicLMzPL\n5GRhZmaZnCzMzCyTk4WZmWVysjAzs0z/HwfFCXnu/4CbAAAAAElFTkSuQmCC\n", "text": [ "" ] }, { "metadata": {}, "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAYgAAAEXCAYAAAC3c9OwAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XmUFOW9//H3F4RBQBbDMBIVXMGICYsYjVfIJNx70Nxc\nRZPrSlyOv4jmYiRmwbgOISYG89MQRI0JEkWPMcSABqPGK7YIcUNRoxAJ6ogb+74K4/f+UT1DT09N\nT81Md1fNzOd1Tp2peurp6m8XRX+76nnqKXN3REREsrWLOwAREUkmJQgREQmlBCEiIqGUIEREJJQS\nhIiIhFKCEBGRUEoQIiISqmgJwszGmdkiM9tpZjMyyk8wsyfNbJ2ZrTazP5rZAcWKS0REwhXzDOJD\nYBJwd1Z5D+BOoF962gLMQEREYmXFvpPazCYBB7n7RfWsHwqk3L1bUQMTEZFa4miDsAbWjwDeKEYg\nIiJSv31ieM96T1nM7AvAdcCpxQtHRETCxJEgQs8gzOwI4K/Ad919YdY6jSgoItJE7t7QlZtQcVxi\nqvNlb2b9gCeBn7j7/cUPSUREshWzm2t7M+tEcNbS3sxK0mUHAvOA29z9rlzbcPe2NxFk1MyyG264\nIf64EjJpX2hfaF/UPzVXMS8xXQdcn7E8BphI8P13KFBhZhXpde7qxSQiEquiJQh3rwAq6ln9k2LF\nISIi0WiojRaovLw87hASQ/tiL+2LvbQv8qPoN8o1RXUvppYQa95ZuvNBW/zsItIslv7+8BbUi0lE\nRFqAOO6DkMbQmYOIxERnECIiEkoJQkREQilBiIhIKCUIEREJpQQhIiKhlCCSzmzvvRAiIkWkBJFn\no0dDu3a5p1deiTtKEZGG6T6IPHPXrQsi0jroDEJEREIpQRTQn/8MVVXBNHhw3NGIiDSOLjEVUHWb\nA6idWURaHiWIpFODhojERJeYREQklBKEiIiEUoIQEZFQShAiIhJKCUJEREKpF1PS6ZnUIrFatGgR\nM2fOZNiwYSxcuJAf/vCHHH744Y2um0ql+Pjjj9m+fTtPP/00F110ESNHjizmR2k8d0/8BHgQavKd\nemr1YBvuc+bsLR8yZG/5yy83YoPVLxKRgrj99tvrXbdz507v16+fr1y50t3dX3rpJT/uuOOaVLdX\nr15+zz33uLv7n/70J+/cubNv3bo1Xx8jVMZ3Z5O+e3WJSUTatNWrV9e7bv78+XTt2pWysjIAjj32\nWJYuXUplZWWj6z7zzDN885vfBODTTz9lz549+f0gBaAEISJSj8rKSj7zmc/ULJsZPXv25M0332x0\n3aOPPprOnTsDMHv2bCoqKujSpUuBP0HzqA1CRArmrbfe4uqrr6Zfv36YGWvXruXmm2+md+/eDb72\ngw8+4KabbsLM2LVrF9u3b+eqq67imGOOqVN32bJlTJw4kX333ZeSkhL23XdfKioq6Nq1a7PiX7t2\nbc2XerVOnTqxZcuWJtV99dVXeeqpp+jatSvjx49vVmzFoAQhIgWxadMmRo4cyeTJkzn33HMB+PnP\nf87IkSN55ZVX6NChQ72vXbduHZdffjl33nlnzSWbyspKTj75ZObOncsRRxxRU3fJkiWccsopPPjg\ng5xwwgmsWrWKE044gf32248bbrih1nZXr17Nr371q+q2TQAWLFjAzp07a5a7du3KNddcA0CPHj1q\n1QXYunUrvXr1qhNzlLqDBw9m8ODB/Pa3v2X48OE888wzyT6LaGrjRWMmYBywCNgJzMhaNxL4J7AN\nmAf0DXl9222kFmmhrr76ai8rK/OqqqqasvXr13uHDh38jjvuyPnayZMn+9SpU+uU/+xnP/MJEybU\nLO/evdv79+/vU6ZMqSlbsWKFl5aW+q9//etIcVZUVNS7bt68eT5o0KBa71dSUuLLli1rVN3nnnvO\ny8rKvLKy0t3dly5d6mbmDz30UKQYm4oW0kj9ITAJuDuz0Mx6AQ8B1wA9CZLIg0WKSUQKaNasWRx/\n/PG0a7f3a6Znz54cddRRzJo1K+dr33nnHebNm1envKSkpFbj7r333ktlZSUXXnhhTdnBBx/M6tWr\nufzyy5v9GYYPH86aNWt4//33gaCheeDAgRx55JEAPPXUU7z22msN1m3fvj3HHHMMffr0qfl8HTt2\nZHDCnwNQlEtM7j4bwMyGAQdlrDoDeMPdH0qvrwDWmll/d19WjNhEJP+2bNnC8uXLGTVqVJ11ffr0\n4YUXXsj5+iFDhnDppZdy3nnnMW3aNHr06MGuXbu4//77ueeee2rqPfDAAxx55JF069Yt758BYJ99\n9mHmzJnceOONnHjiiTz99NM8+ODe37DTpk1j6NChDBo0KGfd4447josuuojbbruNdu3asWDBAubO\nncthhx1WkLjzpqmnHk2ZgJ+ScYkJmAJMy6rzOnBGVpkuMYm0IP/4xz/czPzqq6+us+700093M/NP\nPvmk3tfv2rXLjzvuODcz79Onj8+cOdPHjBnjjz/+eE2dTz/91Dt16uTl5eU+f/58v/baa338+PE+\nevRoX7x4ceRYf/GLXzTuw7UgtJBLTDX5KGu5C7A5q2wz0LyuByISq82bg//WHTt2rLOuulF248aN\n9b6+Y8eOzJs3j1GjRrFy5UrOP/98Nm7cyJAhQ2rqrFu3jl27drFq1SqWLl3KpEmTuPXWWxk/fjwj\nRoxgyZIlkWL90Y9+1JiP1qYUO0FkP1dtK5B9btgdqNuHTERajPbt2wPBvQDZdu/eDUBVVVXObdx2\n22307duXuXPncuihh/Loo48yePBg3njjjVqvr6ys5IILLqh53Ze//GW6d+/OhAkT8vJZ2rJid3PN\nPoN4E6j5lzWzLsDh6fI6KioqaubLy8spLy/Pe4CJo7GYpAUqLS2td922bdswM/bbb79669x6660s\nXLiQv/zlL0Dw//3aa69lypQpjBkzhldffZX9998fgCOOOIKSkpJarz/wwAN58skn2b17d87utK1R\nKpUilUrlZVtFSRBm1h7okH6/9mZWAuwBZgM3m9kZwF+BG4BXvZ4G6swEISLJVVZWhpmxYcOGOuu2\nbdtGjx496u3/7+5MmjSJBQsW1JR17tyZW265hQEDBnDZZZfx1ltvMWDAAEpLS+nZs2edbZSUlLB7\n927Wr19fcx9FW5H543nixInN2laxLjFdB2wHJgBjgB3ANe6+FvgGcCOwHhgGnF2kmESkQLp06cKQ\nIUNYsWJFnXXLly/P2b1zzZo1bNy4saYraaaxY8fSrVu3mhvbTjjhBNasWVOn3q5duygpKcl5JlOf\nRYsWccUVVzBz5kwuvfRS3n777SbVTaVSPPDAA0yfPp0xY8bw1FNPNTqW2DW1dbuYE225F5NGc5UW\n6vrrr/cDDzywVtny5cvdzOqMoLps2TLfsWOHuwe9kw444ABfsGBBnW1u3rzZDzrooJoeUA888IB3\n7ty55rXVr+/evbufddZZjY65pY/emo0W1otJRNqIyy67jG3btnHffffVlE2dOpWBAwfy7W9/u6Zs\n/vz5DBgwgNGjRwNBw/aUKVMYO3Ys7733Xk29zZs3c/HFFzNlypSadoWzzjqLQYMGMW3atJp6s2bN\noqSkhMmTJzc65rY+ems2jcUkIgVxwAEHkEqluOaaa3jllVfYsmULGzZs4PHHH2efffapVa937961\nxlc688wz6dOnD+PHj6ekpAQzo6qqinHjxjFixIiaembGo48+ypVXXsmZZ55Zkzief/55+vbt2+iY\nc43IesghhzSq7tFHH12zrqWM3ppNCSLp1HtJWrBBgwYxd+7cnHX69+/PypUr65QPHz6c4cOHN/ge\nPXv2ZMaMGU2OMVNbH701mxKEiLR6t912W87G5mHDhnHeeefRvXv36nbPGm1q9NYsShAi0uqNGzcu\nUr3Pfe5z3HXXXTXLe/bsYf369fTr169O3aOOOorf/OY3oXWff/55Ro8ezQsvvEC/fv0YPnw4Y8eO\n5YknnuCMM85o/gcqEjVSi4iktfXRW7PpDEJEJK3Nj96axbKvoSWRmQU3Q7SAWE87DR55JJifMydY\nBhg6FBYvDuZffjlYFhEppOqxsNy97qBYEegSU9KZ7R2PSUSkiJQgREQklBKEiIiEUoIQEZFQShAi\nIhJKCUJERELpPoikawFde0WkddIZhIiIhFKCEBGRUEoQIiISKlKCMLMf1FN+ZX7DERGRpIh6BnFD\nPeXX5SsQERFJlpy9mMzsq4AB7dPzmQ4HNhcqMEmrHodJvZlEpMga6uZ6N+BACTA9o9yBVcDlBYpL\nRERiljNBuPshAGY2092/VZSIREQkESLdKJeZHMysXda6T/MdlIiIxC9qL6Zjzew5M9sO7MmYdhcy\nOBERiU/UoTbuAR4BLga2Fy4cERFJiqgJoi9wjbeEZ362NtrlIhKTqPdBzAZGFTIQERFJlqhnEPsC\ns83sWYLurdXc3c9vbhBmdhBwB3Ai8AnwJ2C8u1c1d9siItI0URPEkvSULV/XP34NrAX6AD2BJ4Hv\nAFPztH0REWmkqN1cKwocx0DgCnf/BFhlZo+ny0REJCaREkTIMBs13H1eHuJ4AjjXzJ4B9gdOAa7N\nw3ZFRKSJol5iqh5yo1opwfAb7wOH5SGOCuB/CcZ2ag/83t0fzsN2E2nRItgcMorViSdCx45ZhRqL\nSURiEvUS0yGZy2bWnuAX/tbmBmBmRnAGMQs4HtgPuNvMfuHuEzLrVlRU1MyXl5dTXl7e3LePxdix\n4eWrV0NpaXFjEZHWJZVKkUql8rIta+qtDWbWAfjA3cuaFYBZKUHPqO7uviVdNhqY5O6fTy87QEu4\nDeO00+CRR4L5OXOCZYChQ2Hx4tyvDU0QOoMQkSay9PeHu1tTXh/1ElOY/wDy0Q11LfAxcJmZ/X+C\nM4gLgNfysO3EOPZY6NatbvnChbBnT/HjERFpSNRG6vezijoDnQi6ojaLu7uZnQH8EvgxwRhPTwHf\na+62k+S3vw0vLy2FtWuLG4uISBRRzyCyh/reBixz9035CMLdXwCG52NbIiKSH1EbqVNQM9R3GbBK\nw3wXidoeRCQmUYf77mZm9wI7gQ+BnWZ2r5l1L2h0IiISm6iD9U0FugDHELQ/VP/VUBgiIq1U1DaI\nk4HD3H1benmZmV0IvFOQqEREJHZRzyB2ENw9nakXwSUnERFphaKeQfwOeDJ9n8J7wCEE3VDr6bwp\nIiItXdQE8TPgI+A8giG5PwJ+QTBGkxSS7qQWkZhE7eb6KUEyUEIQEWkjonZznWpmJ2aVnWhmvypM\nWCIiEreojdTnAC9nlb1CcMlJRERaoagJ4tOQuu2AJo0QKCIiyRc1QSwAfpoeaqP6eRATgWcLFZiI\niMQrai+mK4C5wEozew/oSzBE938VKjBJU+8lEYlJ1F5M75vZUOCLwMEEjxp9QQP2iYi0XpEfGOTu\nVcBz6UlERFq5qG0QIiLSxihBiIhIKCUIEREJVW8bRHWX1oaoobrANBaTiMQkVyP1nqxlp+6NcQ60\nz2tEIiKSCLkSxGEZ8/8JfJNgVNcVBPdBXAU8VLjQREQkTvUmCHevrJ43syuBYe6+IV30lpktAhYB\ntxc0QhERiUXURupuBM+gztQZ6J7fcEREJCmi3ih3D/C/ZnYrwV3UfYHvAvcWKjAREYlX1ATxI2A5\ncDbBE+U+BqaiR44WnnoviUhMGvNEuTvTk4iItAFRnyjXzswuMbN5ZvaPdNkIMzuzsOGJiEhcojZS\nTwQuJrik1Ddd9iFBV9e8MLOzzWypmW01s+VmdlK+ti0iIo0XtQ3iImCIu68xs+pure9S+16JJjOz\n/wBuAs509xfNrA96Wp2ISKyiJoh2wNassi7AljzFMRGY6O4vArj7x3naroiINFHUS0yPAbeYWSeo\nGadpEvCX5gaQfnzpsUBvM/uXmb1vZlOr36vNM9s7HpOISBFFTRBXAgcAGwlumtsKHEJ+2iDKgA7A\nN4CTgMHAEODaPGxbRESaKGo3103A6WZWBvQD3s/jZaAd6b9T3X0VgJndQpAgaiWJioqKmvny8nLK\ny8vzFIKISOuQSqVIpVJ52VakBGFmM4A/uPsTwKqM8tvd/TvNCcDdN5jZB1HqZiYIERGpK/PH88SJ\nE5u1raiXmMYAd5vZD7PKv9Wsd99rBnC5mZWaWU/ge+ShfUNERJouaoLYCRwPnG1m95lZSZ7jmAS8\nBCwDlgAvAzfm+T1ERKQRIj9y1N0/IGhEbgc8a2YH5SsId9/j7v/j7j3dvY+7j3f3T/K1/RbNXeMx\niUgsGvVManff4e7nAn8GXgA6FiQqERGJXdQb5SZlLrj7TWb2GqCxmEREWqmo3Vwnh5Q9RnADnYiI\ntEL1Jggze8LdR6Xnn62nmrv7iIJEJiIiscp1BpH5tLjp9dRR66mISCtVb4Jw9/sz5n9flGikrupx\nmNSTSUSKLNclpouJcIbg7nfnNSIREUmEXJeYvkW0S0hKECIirVCuS0zlRYxDREQSJup9EDXMzMh4\n2pu7f5rXiEREJBEi3UltZgea2WwzWw/syZh2FzI4ERGJT9ShNu4kSAZfJXhY0FDgYeCyAsUl1TQW\nk4jEJOolpn8D+rr7VjPD3V9N93L6O3BX4cITEZG4RD2DqL6kBLDBzHoD24ADCxKViIjELmqCeBE4\nJT3/BPAgMBtYVIigREQkflEvMY1hbzL5HvB9oCvwq0IEJSIi8Ys6muvGjPntZA3/LSIirU+kBGFm\nHYBzgCEEZw7V3N0vKURgkqaxmEQkJlEvMc0EPk/w/IeV6TJDo7mKiLRaURPEyQTdXDcXMhgREUmO\nqL2YlgL7FzIQERFJlsb0YrrbzB4DVqXLjKAN4t76XyYiIi1V1ARxAcHd1N2AHVnrlCBERFqhqAli\nPDDE3ZcUMhgJod5LIhKTqG0Qq4AVhQxERESSJeoZxC3ATDObzN42CADc/Z28RyUiIrGLmiCmpf+e\nllXuQPv8hSMiIknR4CWm9BPkjgRK3L1d1pTX5GBmR5rZTjObmc/tiohI40Vtg3gdqCpkIGnTCEaO\nVcusiEjMGkwQ7u7Aq8CAQgZiZmcDG4CnyHjmdZtntnc8JhGRIoraBvE08JiZ/R54n+AXfvWNcnc3\nNwgz6wZMBL4CaPA/EZEEiJogTgIqgS+HrGt2giAYPvx37v6RmenykohIAkR9HkR5oQIws8HASIKh\nxCHH5aWKioqa+fLycsrLCxaWiEiLlEqlSKVSedmWecQ7dc2sJ3Aq8FngQ2Cuu69vdgBmVwA3AlvS\nRV0Jus4ucfdh6ToOEDXWOJ12GjzySDA/Z06wnEtpKaxdG8yvXh0s16LnQYhIE1n6+8Pdm9SQGakX\nk5l9CXgbGAt8AbgUWG5mJzblTbPcBRwGDAIGA3cCjwKj8rBtERFpoqhtEFOA77j7H6oLzOysdPlx\nzQnA3XeQMQCgmW0Fdrj7uuZst9XQmYOIxCRqgugP/DGr7CHgN/kNB9x9Yr63KSIijRf1Rrl/ETyT\nOtN/A8vzG46IiCRF1DOIK4BHzexyglFd+xGcVXy9UIGJiEi8onZz/buZHQ78J0EvpkeAx9ROICLS\nekU9gyDdpVWD6ImItBE5E4SZPd3A693dv5rHeCSb7oMQkZg0dAZxf0iZAwcStEt0zntEIiKSCDkT\nhLv/LnPZzHoBVwHfBh4EflK40EREJE5R76TubmaTCLq1HgAMdfdL3P2DgkYnIiKxyZkgzKyzmf0Y\neAc4Gvg3dx/j7m8XJToREYlNQ20Q7xIkkcnAIqDMzMoyK7j7vALFJiIiMWooQVSPkXRpjjqH5ikW\nCaPeSyISk4YaqQ8pUhwiIpIwUcdiEhGRNkYJQkREQilBiIhIKCUIEREJpQSRdGZ7x2MSESkiJQgR\nEQmlBCEiIqGUIEREJJQShIiIhFKCEBGRUJEfOSox0VhMIhITnUGIiEgoJQgREQmlBCEiIqGUIERE\nJFQiEoSZdTSz6WZWaWabzWyxmZ0cd1wiIm1ZIhIEQW+qFcAId+8GXAv80cz6xRtWAmgsJhGJSSK6\nubr7dmBixvKjZvYuMBR4L7bARETasKScQdRiZmVAf+DNuGMREWmrEpcgzKwDcD/we3dfFnc8IiJt\nVaIShJm1A2YCO4FxMYcjItKmJaINAsDMDJgOlAJfc/eq7DoVFRU18+Xl5ZSXlxcrvDpWr4ZNm+qW\nb91a/FhERKqlUilSqVRetmWekLF+zOxOYBDw7+6+LWudAyQlVoBx42DatNx15syB007LXae0FNau\nDeZXrw6WRUTywdI9IN29SV0hE3GJKd2d9RKCBLHSzLakp3NiDk1EpM1KxCUmd3+PhCSrpujVC3r0\nqFvepUvxYxERyZdEJIiW7oYbgktOIiKtSYv91S4iIoWlBCEiIqGUIJJOYzGJSEyUIEREJJQShIiI\nhFKCEBGRUEoQIiISSglCRERC6Ua5pEvQ+FMi0rboDEJEREIpQYiISCglCBERCaUEISIioZQgREQk\nlBJE0mksJhGJiRKEiIiEUoIQEZFQulEuQb7yFejQoXbZ4nhCERFRgkiSN9+MOwIRkb10iUlERELp\nDCJm8+bBnj21y3btgi99KZjvVOLs3Fn8uERElCBi9vnP1y1TQhCRJNAlJhERCaUEISIioZQgREQk\nlBKEiIiESkSCMLP9zWy2mW01s0ozOyfumJJi5y6NxSTSVu3YAT175p6GDSvc+yelF9M0YCfQGxgC\nPGpmr7n7knjDEpE4zZoFq1fnrnPhhdClS1HCicXGjbnXb95cuPeOPUGYWRfgDGCgu28HFprZw8C3\ngB/HGlxCpVIpysvL4w4jEbQv9krSvpg5E375y9x1vv99OP/8YP766+GZZ+rWmT+/4fc6/fS6CSJJ\n+6Iliz1BAP2BPe6+PKPsNaA8nnBq++c/4fXX65b/61/FjWP9+uD+iCVL4A9/SNGtW3lova5doX//\n3NvatQs++aRu+apVsGVLsI36fPazDf9aO+kk2Lq1/vVdusDChbm3EVX1F8H27Q3fP9KlC5SU5Od9\nk+jSS1Ns2VIeuq76+LnsMjj44LrrJ0yAdnm84Lx2bfj/m0wXXAA/Tv8E/Oij5r+nOyxdGszPmpWi\nd+/yOnU6dIAjjwzmN22CN97Ivc1evWDAgObHlg8lJfDxx8H88uXwxS8W/j2TkCC6AtknSVuA/WKI\npY6HH4arroo7CvjMZ2ovT59ef92+fYO/K1YEf9u1g/bt967fvbt5sXTqtHe++ku5UycYODCYf/nl\nhrcxalTt5b//PUgqp5yS+3WPPRb8HTIk+FwffQRz50Z7T4Cvfz33+rlzg7/VnyXTBx8EXyojRgSf\nd8MGeOmlIGmWldWtv3hx7fd87jlYt67uoIx/+1vwd/jwus1N7vDss9Fif+ut3OsB7rgjvHz69NrH\nyLJl0d6zPtX7sSGNSQzHHgvHHx/Mz5gRXJ8H+MIXgtirqoL9W+3228O3U/3F+uKL0d63qfsgH6qq\n9s6bBW0OAN27F+f9zd2L8071BWA2BFjg7l0yyn4AjHD3U9PL8QYpItKCuXuTerokoRfTMmAfMzsi\no2wQ0MDJn4iIFFLsZxAAZvYA4MD/A4YCc4EvufvSWAMTEWnDknAGAfAdYF9gNXAfcKmSg4hIvBJx\nBiEiIsmTlDMIERFJmEQliKhDbpjZhWZWZWZbMqYRxY63UMxsnJktMrOdZjajgbrfM7OPzWyTmU03\ns47FirMYou6L1n5MAJhZx/S/caWZbTazxWZ2co76rfbYaMy+aO3Hhpndl/533mxm75jZNTnqNuqY\nSFSCoPaQG+cBd5jZ0fXUXeju+2VMEe65bDE+BCYBd+eqZGajgAnAV4F+wGHAxIJHV1yR9kVaaz4m\nILhvaQVBF/BuwLXAH82sX3bFNnBsRN4Xaa352Pg5cGh6P5wCXB6WLJtyTCQmQWQMuXGdu29394VA\n9ZAboS8pWnBF5u6z3f1hYF0DVS8AfufuS919I/AT4MJCx1dMjdgX0IqPCYD0/4uJ7r4ivfwo8C5B\nz79srfrYaOS+gFZ8bLj7m+6eOY7AHoIOP9kafUwkJkFQ/5AbIfe04sAQM1tjZm+Z2bVm1j6kXkvX\n0EF9NME+qvY6UGZmPQsXUmwa2hdt5ZioYWZlBP9v3gxZ3ZaOjYb2Ras/NszsdjPbRvD5f+rur4RU\na/QxkaQE0ZghN+YTDO5XCnwDOAf4YWHDi0VDXcy6Apsylqv3XyKGKcmzhvZFWzkmADCzDsD9wO/d\nfVlIlTZzbETYF63+2HD37xD8m/878FMzCxupqdHHRJISxFagW1ZZd4IkUYu7v+vu76Xn3yA4Vfpm\nwSMsvoZ+NWfvs+oRWurss1Yg575oQ8cEZtYOmEnQXjeunmpt4tiIsi/ayrHhgRQwiyAJZmv0MZGk\nBNHcITda4zXGhn41vwkMzlgeBKxy9w2FCyk2Tblhp9UdE2ZmwHSgFPiGu1fVU7XVHxuN2BehLy9M\nVInQAdgWUt7oYyIxCcLdtwF/Bn5iZp3N7CTgvwh+HdRiZqekrzliZkcR9GCYU8x4C8nM2ptZJ4Ke\nGu3NrKSea6b3Aheb2efS1xGvA3J2i21pou6L1n5MZLgDOAo41d135ajX6o8NIu6L1nxsmFmpmZ1t\nZl3S/1dGAf9N0MEnW+OPCXdPzAT0BGYTnApVAmeny/sSnAYdlF6+GViZrvc2UAG0jzv+PO6HCuDT\nrOn67P2Qrvu99L7YRPBrqkPc8cexL1r7MZH+jP3Sn397+rNXT+e0tWOjMfuiNR8bQC8gBWwANgIv\nEiTMOt+bTTkmNNSGiIiESswlJhERSRYlCBERCaUEISIioZQgREQklBKEiIiEUoIQEZFQShAiIhJK\nCUJEREIpQYiISCglCBERCaUEIZKWfr7xD8zs9fRzi6ebWZmZPZZ+hu+TZtYjXXdW+tm+G83smcxH\n45rZ18zszfQzgj8ws++ny3uZ2Vwz22Bm68xsfnpEUpFEUoIQ2csJHns7EhgAfB14DLiK4Dnp7YDv\npuv+FTiCYKjpVwgeWFNtOnCJB88IHgjMS5d/H3ifYIC13sCPXYOhSYIpQYjUNtXd17j7R8CzwHPu\n/poHw0lDKUtbAAABQUlEQVTPBoYAuPsMd9/m7rsJHvw+yMyqn8z1CTDQzLq5+yZ3X5xR3gc4xN2r\nPHjuukhiKUGI1LYqY35H1vJOoKuZtTOzm8xsuZltAt4lOPvola73DeBrQKWZpczshHT5zcBy4G9m\n9raZTSjoJxFpJiUIkdzC2gjOBU4FRrp7d+DQdD0DcPdF7j6a4PLTHOCP6fKt7v4Ddz88/forzeyr\nRfgMIk2iBCHSePsBu4D1ZtYF+Fn1CjPrYGbnmVl3Dx6BuQWoSq/7upkdkW6Y3pwub8xjMkWKSglC\nJDfPmneCRze+B3xI8Mz057LqjQHeTV9+ugQ4L11+BPAkQdL4OzDN3Z8paPQizaAnyomISCidQYiI\nSCglCBERCaUEISIioZQgREQklBKEiIiEUoIQEZFQShAiIhJKCUJEREIpQYiISKj/A2iKwzUXGrCH\nAAAAAElFTkSuQmCC\n", "text": [ "" ] }, { "metadata": {}, "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAZAAAAEXCAYAAACDChKsAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XmYVNW19/HvoplHQQYNkwOiIoo4XjVii8+VaN6oMRo1\nSjAxxiHkKhpnDY0aBZMregmaKCgRxWBUNEqirzfSopCI84AYFGQSGkGUGZR23T9OdVNVXVV9qrrG\n7t/neerpOvvsqr36UPSqc/Y+e5u7IyIikq5mhQ5ARERKkxKIiIhkRAlEREQyogQiIiIZUQIREZGM\nKIGIiEhGlEBERCQjeU0gZvawma0ysw1mttjMbkhS73wzqzazjVGPIfmMVUREUmue5/ZuB37m7tvM\nbF/gJTN7w92fS1B3jrsraYiIFKm8JhB3nx9XtAP4LEl1y3E4IiLSAHnvAzGze8xsMzAfuNXd30xQ\nzYHBZrbGzP5tZjeaWVl+IxURkVSsEHNhmZkBxwGPAye7+7y4/XsC37j7UjMbCEwHprr72LwHKyIi\nCRUkgdQ2bnYvsM3dR9VT7yzgKnc/LKpMs0CKiGTA3bPSRVDoYbwtgM0h66pPRESkiOQtgZhZNzM7\n28zamVmZmQ0DzgSeTlD3JDPrEXm+H3Aj8FSi93V3PdwZPXp0wWMoloeOhY6DjkXiR7bl8wzEgYuB\nFcDnwC3AcHd/zcz6RO716BWpOxR4x8w2ATOBJ4Db8hiriIjUI2/DeN19LVCeZN8yoEPU9lXAVfmJ\nTEREMlHoPhDJkvLy8kKHUDR0LAI6DjvpWORGQUdhNUTNKKxSjV9EJN+COygazygsEREpUUogIiKS\nESUQERHJiBKIiIhkRAlEREQyogQiIiIZUQIREZGMKIGIiEhGlEBERCQjSiAiIpIRJRAREcmIEoiI\niGRECURERDKiBCIiIhlRAhERkYwogYiISEaUQEREJCNKICIikhElEBERyYgSiIhIGl5//XUuu+wy\npk6dysUXX8yiRYtS1t+4cSNnnHEGy5cvT2tfKchbAjGzh81slZltMLPFZnZDirqjInXXm9lkM2uZ\nrzhFpGm79957k+7bvn07Z5xxBtdffz3Dhw/nZz/7Geecc07S+pMnT+bOO+/kySefxN1D7ysVlq/A\nzewAYJG7bzOzfYGXgPPd/bm4esOAPwHHA6uAGcC/3P26uHoOlOyBF5HsuuYa+OKLxPu6d4dbbw33\nPmPGjGH06NEJ973wwguMGjWK999/Hwj+/nTs2JH33nuPPfbYI+l7NmvWjCVLltCnT5+09mWbmQHg\n7paN92uejTcJw93nxxXtAD5LUHUEMMndFwCY2c3ANOC6BHVFRACYNg1WrEi8r1+/8AkklSVLlrDr\nrrvWbpsZnTt3Zv78+SkTSGOVtwQCYGb3ECSIVsBId38zQbUBBGcdNd4FephZZ3dP8v1CRErB4sWL\nGTNmDG3atKGsrIzNmzczbtw4evToUe9rV6xYwdixYzEztm/fzpYtW7j22msZOHBgfCvAGKANUAZs\nZseOcUD9bdRn7dq1tG3bNqasdevWbNy4scHvXYrymkDc/VIz+wVwHPC4mb3p7vPiqrUH1kdtb4j8\n7AAogYiUqE8++YQjjjiCSZMmcdpppwHwyCOPMGzYMN544w3KysqSvvbzzz/nl7/8JX/4wx/o0aMH\nEybAzJlLOOqo73D00c/Svn0/1q4F+AQ4gvPOm0T//qfx618DPMLKlcOorq7bxmeffcZdd90Vcyn8\nlVdeYdu2bbXb7du354Ybgi7bXXbZpc5l802bNtG1a9eGHJrS5e4FeQD3AuMTlL8NnBG13RX4Bugc\nV88BHz16dO1j1qxZLiLF6dRTT/Vdd901pmzLli3eunVrnzJlStLXPfWU+2673eG77z7BDzzQ/cAD\n3aHmcZvDNVHbpzrs6suWuS9cWFO2xc1StxGtoqIi6b4XX3zRBw0aVLv99ddfe6tWrXzhwoUp39PM\nfOnSpWnvy4ZZs2bV/o2s+bvpWfo7ntczkDgtgM8TlM8HDgYej2wPAlZ7kstXFRUVOQlORLLnq6++\nYubMmRx66KEx5W3atKFv3748/vjjjBgxIuFrv/gCqqoWA6tZtWpk3N5WBN2pAF8BM4FD4+q0oXnz\n1G2Edeyxx7JmzRqWL19O7969eemllzjggAPYZ599APjHP/5B165dGTRoUJ3XeooBP6n2NVR5eTnl\n5eVAMEAgm/IyjNfMupnZ2WbWzszKIiOtzgSeTlD9IeACM9vfzDoDNwEP5iNOEcmNdevWUV1dTevW\nrevs69ixI/PmxV/JjjcYeAo4F/gyUrYdeITjjjufxx+H++9fh1k1Awa0JqqfG4BmzcK0Ub/mzZsz\ndepUfvOb3/DQQw/x8MMPM3369Nr9EydO5JlnnqndnjZtGpdeeilmxrXXXsvEiRND7SsZ2TqVSfUg\nuAxVSdCH8SUwDzglsq8PsBHoFVV/FFBF0BcyGWiR4D1rTsVEpMjt2LHD27Vr50ceeWSdfT179vRm\nzZp5dXV1wtc++KA7bHc43MG8W7fd/Te/merf/e55fs89z/nnnyduY+clLPeystRtRBs3blzGv2ex\nI8uXsPJyBuLua9293N07u/su7n6Eu/81sm+Zu3dw9xVR9ce7+27u3sndL3D3r/MRp4jkRllZGWed\ndVadu7arqqpYuXIlEJylJNcSeJGePYexdm0VN974Y8y+5Ac/GEyXLqnbgCqqq8O0Ebj66qtD/lai\nqUxEJC/GjRtHly5duO+++wDYsWMHEydOZPDgwbh7ylFYgd/Trl0fnn32Wfbcc09mzpzJwQcfXHtT\nX6I2gv6RibRqFbYNSYcSiIjkRdeuXXn11VepqqriyiuvZMyYMVx44YWYGW3atKFz584pXj0emMNR\nR/2Rk08+mffee4/LL7+c1atXc9555yVs4/bbryS4H+RCIEwbkq68TWWSbZrKRKRx6N27N/vttx8v\nvPBCwv0PPuj89Ke7Aq8wYsQApkzZue+Pf/wjl1xyCQsWLGDfffeNed1HH0H//sHz5s17U16evI2m\nIttTmegMREQKZs2aNXz66aeceeaZSets3LiGYOzNPnX2XXTRRXTs2DHmxr8ErbBjR+o2Ukln9t1X\nX32VO++8k4qKCk488URmz54dal/JylZvfL4faBSWSEkZP368d+/e3T/99NPast/+9rfeu3dv3759\ne0zdhQsX+tatW93d/YEHvnHYzeEVHzEi9j03bNjgvXr18q+++qpOGztHYf3Wmzev20YY27Zt8759\n+3pVVZW7u7/22mt++OGHJ6y7efNmv/baa2u3//KXv3jbtm195cqVSfdFH4t8oBRHYYmIbN68mXbt\n2tV2ZL/11lv87ne/47HHHqNly50rNsyePZt99923drqT4LLL3cBFbNq0tLbehg0buOCCC7j77rtp\n0aJFwjbgLeB37LZbbBthzZ49m/bt29fO1XXooYeyYMEClixZUqfuxx9/zLhx41i8eDEAJ554Ilu3\nbmXOnDlJ982dOzftmIpJIe9EF5Em5Morr2T9+vVcd911bNy4kfXr1/P0009z5JFHxtTbbbfd6N69\nO/369Ysq/SGwO/PmXc7ZZ7fCzKiurmbkyJEMGTIkYRurVm0kuJXsaVq3jm0jrHRm3z3ooIOYO3cu\ne+21FxBM/giwzz77pNxXypRARCQvWrduzR133FFvvf79+1NVVZVgz7EMHXpsTCd6qjaiO9Ezle7s\nu//xH/9R+/z222/nyiuvrJ3WJNW+UqUEIiJNzu9///uUneGHHXYY5557Lp06daoz0jPM7LuTJ0+m\nZ8+ejB07Nq19pUYJRESanJEj4ydlTGz//fePuikxuPlx3bp19O3bN+lrZs6cSbNmzRg7dizbt2+n\nqqqqtn6qfaVInegiIklEz74LJJx995133qmt/9JLL7F69WpOPvlkqqqq+Nvf/saqVavq3VeqdAYi\nIpJE9Oy7Rx99NLNmzaoz++4hhxzCoEGDWLx4Md/73vfYtGlT7X4zY/369Sn3lTLdiS4iRW3KFPjJ\nT4LnI0aQshM9WnQner9+wXZTpzvRRUSkKCiBiIhIRpRAREQkI0ogIiKSEY3CEpGisX493HVXbNlb\nbzX8fbdvh6jRtgB06QK9ezf8vZsyJRARKRrr10NFRfbfd/lyOPjg2LILL4SoewQlA7qEJSIiGdEZ\niIgUpQ4dYNSo2LL4s4hUWraEgw6KLfv8c/j004bHJgElEBEpSrvsAmPGZP76vn3r9nvcdx9cdFHD\n4pKd8nYJy8xamtlkM1tiZhvM7C0z+06SuuebWbWZbYx6DElUV0RECiNUAjGzXyUpvyKNtpoDy4Ah\n7t4RuBF4zMySTUU5x907RD0awQLCIiKNR9gzkNFJym8K25C7b3H3Me6+LLI9E/gEOCTJS7IyV4uI\niORGyj4QMxtK8Ie8LPI82t7AhkwbNrMeQH9gfoLdDgw2szXAOmAqcLu7V2fanoiIZFd9negPEPwx\nbwVMjip3YDXwy0waNbMWwCPAFHdfmKDKbOAAd19qZgOB6cAOoPSX8BIRaSRSJhB33wPAzKa6+/Bs\nNGhmzQjOKLYBCZcFc/dPop6/b2Y3A1eRIIFURN11VF5eTnl5eTbCFBFpFCorK6msrMzJe4caxhud\nPCIJIHrfN2Ebs2Ay+slAN+DkNC9JJewTqcjFbasiIo1E9BfrMQ0ZF51A2FFYh5rZP81sC8GlpJrH\n12m2dy+wH3CKu29P0d5JkT4SzGw/ghFbT6XZloiI5FDYGwn/BPwVuADYkklDkeG6Pye4dFVVszJW\npGwOQWf6/u6+AhgKPGhm7Qn6WqYCt2XSroiI5EbYBNIHuMEbsH6suy8l9RlPh6i6VxH0eYiISJEK\nex/IDGBYLgMREZHSEvYMpA0ww8xeJrikVMPd/cfZD0tERIpd2ATyQeQRL+NLWiIiUtrCDuOtyHEc\nIiJSYkIlkATTmNRy9xezF46IiJSKsJewaqY0qdGNYHqT5cBe2Q5KRESKX9hLWHtEb5tZGcHNfZty\nEJOIiJSAjBaUikxBchtwdXbDERGRUtGQFQn/E9D06iIiTVTYTvTlcUVtgdbApVmPSERESkLYTvT4\nqdw3AwvdfX2W4xERkRIRthO9Emqncu8BrE5nGncREWl8wk7n3tHMHiKYSfdTYJuZPWRmnXIanYiI\nFK2wnegTgHbAQIL+j5qfE3IUl4iIFLmwfSDfAfZy982R7YVmdj6wOCdRiYhI0Qt7BrKV4O7zaF0J\nLmmJiEgTFPYMZBLwgpn9N7AU2AMYBdyfo7hERKTIhU0gtwErgXOB3SPPxxHMkSUiIk1Q2GG83xAk\nCyUMEREBwt+JPgF41N3nRpUdDfzQ3S/PVXAi0rgNHw4rV+7c3qZe1ZIS9hLWOcCv4sreBJ4GlEBE\nJCP//CcsWlToKCRTYUdhfZOgbjPAshuOiIiUirBnIK8At5rZVe7+TWQ9kDHAy2FebGYtgXuBE4Au\nwCLgOnd/Lkn9UQRTxbcFHgcucfevQsYqIiXoT3+Cnj13brdqVbhYJJywCeQy4FmgysyWAn2AVcD3\n0mhnGTDE3ZeZ2XeBx8zsQHdfGl3RzIYB1wDHR9qYQZCsrgvZloiUoKOPhn79Ch2FpCPsKKzlZnYI\ncATQm2Ap21fDTqjo7lsIkkDN9kwz+wQ4hOC+kmgjgEnuvgDAzG4GpqEEIiJSVMKegdSsQvjPyKNB\nzKwH0B+Yn2D3AIKzjhrvAj3MrLO7f9HQtkVEJDsasiJhRsysBfAIMMXdFyao0h6IXmdkQ+Rnh1zH\nJiIi4YU+A8mGyHoiUwnm0BqZpNomoGPUds2U8RsTVa6oqKh9Xl5eTnl5eUPDFBFpNCorK6msrMzJ\ne+ctgZiZAZMJJmU8OXJJLJH5wMEEo68ABhEsYJXw8lV0AhERkVjRX6zHjBmTunKakl7CMrNmYR5p\ntHUvsB9wirtvT1HvIeACM9vfzDoDNwEPptGOiIjkQaoEsCPu8XWSsnqZWV/g5wRnE1VmtjHyOMfM\n+kSe9wJw9+eBO4BZwBKCe0ZGZ/C7iYhIDqW6hLVX1PPvAmcQzMq7jOA+kGuBJ8I0ErnXI1Wyiukg\nd/fxwPgw7y0iIoWRNIG4+5Ka52Z2BXBYVD/Ev83sdeB14J6cRigiIkUpbB9GR4JpRaK1ZecIKRER\naWLCjsL6E/C/Zjae4C70PsB/EXR4i4hIExQ2gVwNfAycTbAi4SpgAlrSVkSkyUpnRcI/RB4iIiLh\n+kAi93z83MxeNLP3ImVDzOyHuQ1PRESKVdhO9DHABQSXrPpEyj4lGMorIiJNUNgE8hPg/7n7owSr\nEwJ8Quy9IiIi0oSETSDNCCY5jNaOJBMciohI4xc2gfwduNPMWkPtrLq3AM/kKjARESluYRPIFcBu\nwJcENxVuAvZAfSAiIk1W2GG864HvR1YS7Assd/dVOY1MRESKWqgEYmYPAn+OzJS7Oqr8Hne/NFfB\niUjjMW0azJ0bW7ZmTWFikewIeyf6ecCJZnaXu/82qnw4oAQiIvWaNQsmTSp0FJJNYftAtgFHAmeb\n2cNm1iqHMYmISAkIvaStu68ws28TLEv7spmdnruwRKQxO/tsOOaY2LKuXQsTi2QurTXR3X0r8CMz\nuxZ4FWiZk6hEpFEbOhQuvLDQUUhDhU0gt0RvuPtYM3sH0FxYIiJNVNhhvHckKPs7wQ2GIiLSBCVN\nIGb2vLsPizx/OUk1d/chOYlMRESKWqozkOjVBicnqeNZjEVEREpI0gTi7o9EPZ+Sl2hERKRkpLqE\ndQEhzjDc/YEwDZnZSOB8YCDwqLv/JEm98wnOeLZEFX/X3WeHaUdEJIyvv4ZNcXOMt2oFLVoUJp5S\nlOoS1nDCXaIKlUAIFqC6BRgGtKmn7hz1rYhILk2ZEjyiPfpocI+KhJPqElZ5Nhty9xkAZnYY0Kue\n6pbNtkVEJPvCTmVSywLNah4ZtFlfcnBgsJmtMbN/m9mNZlaWQTsiIjGaN4d27WIfZfrrkrFQCcDM\neprZDDNbB+yIenydQZv1XRabDRzg7t2AHwDnAFdl0I6ISIyf/jTo94h+nHFGoaMqXWHvRP8DsBUY\nCrwEHAeMJrMbCVOegbj7J1HP3zezmwkSyNhE9SsqKmqfl5eXU15enkFIIiKNU2VlJZWVlTl577AJ\n5Bigj7tvMjPc/e3IKK25wH1ptpnJvSNJk050AhERkVjRX6zHjBmT1fcO24dRc8kK4Asz6w5sBnqG\nbcjMyiJrqjcHysysVaK+DTM7KbLyIWa2H3Aj8FTYdkREJD/CJpB5wEmR588D04EZwOtptHUTwb0d\n1xAsULUVuMHM+pjZRjOrGZk1FHjHzDYBM4EngNvSaEdERPIgnRUJa5LNKOBKoD1wV9iG3L0CqEiy\nu0NUvatQp7mISNELOxvvl1HPtxA3vbuIiDQ9oRKImbUgGE47mODMo4a7+89zEZiIiBS3sJewpgIH\nEgzbrYqUGZqNV0SkyQqbQL5DMIx3Qy6DERGR0hE2gSwAugBKICJSr5dfhvHjY8vefLMwsUjupDMK\n6wEz+zuwOlJmBH0gDyV/mYg0RStWwIwZhY5Cci1sAhlBcDd6R4L7N6IpgYiINEFhE8jlwGB3/yCX\nwYhI43PMMXDFFbFlgwcXJhbJrrAJZDWwLJeBiEjj1Ls3nH56oaOQXAibQO4EpprZHezsAwHA3Rdn\nPSoRESl6YRPIxMjPU+PKHdByLCIiTVC9kymamQH7AK3cvVncQ8lDRKSJCjsb77tAdS4DERGR0lJv\nAnF3B94G9s19OCIiUirC9oHMAv5uZlOA5QR9HzU3Ej6Qo9hERKSIhU0g3waWEKyFHk8JRESkCQq7\nHkh5juMQEZESE/YMBDPrDJwCfAv4FHjW3dflKjARESluoUZhmdlRwCLgIuAg4GLgYzM7OoexiYhI\nEQt7BnI3cKm7/7mmwMzOipQfnovARESkuIW9D6Q/8Fhc2RMENxiKiEgTFDaBfESwJnq0M4GPsxuO\niIiUirAJ5DLg92b2LzN7zMxeBe6JlIdiZiPN7HUz22ZmD9ZTd5SZrTKz9WY22cxahm1HRETyI1QC\ncfe5wN4Ekyq+AUwA+rn7nDTa+hS4hXruGzGzYcA1wFCgL7AXMCaNdkREJA9CD+ONDNmdmmlD7j4D\nwMwOA3qlqDoCmOTuCyL1bwamAddl2raIiGRfygRiZrPqeb27+9A027R69g8AoldTfhfoYWad3f2L\nNNsSEZEcqe8M5JEEZQ70JOj/aJtBm17P/vbA+qjtDZGfHQAlEBGRIpEygbj7pOhtM+sKXAtcCEwH\nbs6gzfrOQDYBHaO2O0V+bkxUuaKiovZ5eXk55eXlGYQkItI4VVZWUllZmZP3DtUHYmadgF8BvwSe\nBQ5x90UZtlnfGch84GDg8cj2IGB1sstX0QlERERiRX+xHjMmu+ORUo7CMrO2ZnYdsJigb+IYdz8v\nk+RhZmVm1pogaZWZWSszS7Si4UPABWa2f2T+rZuAlMN+RUQk/+o7A/mEIMncAbxO0JndI7qCu78Y\nsq2bgF9HbZ8HVETWGJkP7O/uK9z9eTO7g2ANkjYEZyKjQ7YhIiJ5Ul8C2Rr5eXGKOnuGacjdK4CK\nJLs7xNUdD4wP874iIlIY9XWi75GnOESkRH31FcybF1v24YeFiUXyK/SNhCIiiXzxBRx7bKGjkEII\nOxeWiIhIDJ2BiEjWNG8ORx4ZW7bvvoWJRXJPCUREsqZLF3jllUJHIfmiS1giIpIRJRAREcmIEoiI\niGRECURERDKiBCIiIhlRAhERkYxoGK+ISMQHH0D80hlHHw0tWxYknKJn7vUtz1GczMwBSjV+kVJ1\n0EGwYMHObXeorg6ed+8Oq1cXJq5MnX02TJ+efP/q1cHv1RiYBev5uXt9C/uFojMQEUnLjh3BQ0QJ\nRESatAED4LjjYsvmzFGSDEOXsEQkLQMG7LyE9fbbwXa0Fi3yH1O2de8Oa9YEz3UJKzmdgYhIxlq0\naBwJQzKjYbwiIpIRnYGISFIVFfDww7FlS5cWJBQpQkogIpLU2rWwaFGho5BipUtYIiKSEZ2BiEgo\nv/41DB8eW9anT2FikeKQtzMQM+tiZjPMbJOZLTGzc5LUO9/Mqs1sY9RjSL7iFJHEunWDfv1iH5ri\no2nL5xnIRGAb0B0YDMw0s3fc/YMEdee4u5KGiEgRy8sZiJm1A04HbnL3Le4+B3gaGJ7sJfmIS0RE\nMpevS1j9gR3u/nFU2TvAAQnqOjDYzNaY2b/N7EYzK8tLlCIiElq+LmG1BzbElW0EOiSoOxs4wN2X\nmtlAYDqwAxib2xBFRCQd+Uogm4COcWWdCJJIDHf/JOr5+2Z2M3AVSRJIRUVF7fPy8nLKy8sbHq2I\nSCNRWVlJZfwiJ1mSl8kUI30g6wjOLD6OlE0Flrv79fW89izganc/NK5ckymKZNGbb8If/xhb9sor\nwSJLABMmwMiR+Y+rEDSZYjh5OQNx981m9iRws5n9DDgE+B5wVHxdMzsJeNPdV5vZfsCNwGP5iFOk\nKVuyBO67r9BRSCnJ553olwJtgM+Ah4GL3X2BmfWJ3OvRK1JvKPCOmW0CZgJPALflMU4REQkhb/eB\nuPsXwPcTlC8jqjPd3a8i6PMQkQI5+GC4+OLYsmOOKUwsUrw0lYmI1LHnnnDRRYWOQoqdJlMUEZGM\nKIGIiEhGdAlLpAn66it47rnYsnnzChOLlC4lEJEmaPNmOPXUQkchpU4JRKQJ+Oab1NsimVACkZy4\n91644ork+x99FE47LX/xNHWHHgpvv514X/PmcNJJsWWHH577mKT0KYFITlRXw7ZtyffrG3DxaN8e\n/vrXQkchpUijsEQameuvh112iX3En32YxT5EMqEzEMm5iy6C8ePh7LP1TTcftm2D9euT73/zTRg8\nOH/xlLof/hBatdq53bev5gyroQQiOde8ObRpA2VaFkxK0EsvxW4PGFCYOIqREoik7YUXgvsIEmnf\nHo47Lr/xSHK33gq/+EVsWfv2hYlFGh8lEEnbeefBZ58l3jdgAMyfn994JLk2bYI+EEnP9OmxX5KW\nLKk7uaQogYiUjB074M9/Trzvww+hbVs4/nhYuTK/cTVGxx8fu60vRYkpgZQQ9+CRSrM8j6s74YTg\nm9rLLwfb69bBPffA7Nn5jaMxWr4ctm/fub1lCwwfXrh4ROIpgZSQ5cuDESCpXHhh7Pb558PRR+cs\nJB55BNauhYEDg+2qqrrX3CUzZ5yh+amkuCmBNDL33x+7fcwxuU0gkh3PP193tE99yePcc4M1y5cu\nhUMOgZYt69bZfffsxSgSTwmkhJnVf0krHzp3hksuSb7/2GPzF0upqqyEsWNT19l7753PW7eGhx/O\naUgi9VICKVG9e8OyZbBwYfDHJ9qDD8K//pW/WL71raDfI1cefRQmTky+/447Suss64UXYNiw2LL6\nvgj8619w5JG5i0kkE0ogJa5//+ARbc6c/CaQXFu5Mvidkvnii/zFkq4dO+qegdX3b9O+PVx3XWxZ\nr17ZjUsy9/XXdUe6tWsHnToVJp5CUgIRyaL/+R9YvHjn9o4d6Sfz66+vm0CkeHz0EfTsGVs2ciRM\nmFCYeAopbwnEzLoAk4H/BNYC17n7o0nqjgKuBtoCjwOXuHuSe5+lKTnrrOA/61VX7fzDXFFR9xLa\nE08E/QSZuvvu5JeVevaEM8+E//1f+Mc/YvfV148Rb+jQ4JJWNE1uKKUin2cgE4FtQHdgMDDTzN5x\n9w+iK5nZMOAa4HhgFTADGAPoO1mObdgAn3+efH+3brmfBuPDD+FHP4otW7165/OePeHb34YuXXaW\nvf563fdJNl38c8/BM8/EltUkn+i5uqqrk8fYuzccdhhMnQoPPZS8XiJz58Zud+qU/3t3JH3Nm9cd\n0bZpE2zcWJh4ikVeEoiZtQNOBw5w9y3AHDN7GhhO3cQwApjk7gsir70ZmJagnkSprKykvLw8dP3V\nq+FXv4otq29Uz7RpcM456ceWSvy3/C1b4K23GvqulSxaVM5HHwU33vXqtfNs5N13k78qVdKItnw5\n7LVX/fW5FYvrAAAIG0lEQVS+/30YMmTnthkcdVS4NrIh3c9EY9bQY7HvvnX7PSZMgP/6r4bFVery\ndQbSH9jh7h9Hlb0DlCeoO4DgrKPGu0APM+vs7kXcXZp9F10U3JhXY8uW5HWT/Qe56y548sng+bp1\n8OqrMGlSMN13usNAf/QjuOAC2Lo1vdclMn58cJlp2rTMXj96NFx6aWzZ6afXzF9UyUEHldeWL1yY\naZSBkSOD+y2SregHwT0YFRWxZeeeC336NKzthlAC2UnHIjfylUDaAxviyjYCHZLUjV7NoOZ1HYBG\nm0Buuw2eeiq27LXXGv6+b79d9w/fiBHhXrvHHsEkctGykTwg+IMcxhtvxG537x78POKIunXTnS6+\nY8fguENwJlRdDaecElyiin/fefOC9UySuewyuPzy9NoXKXXmebgTzcwGA6+4e7uosl8BQ9z9lLi6\nbwO3uvvjke2uwGfArtFnIGZWBLfQiYiUHnfPylCNfHXfLQSam1m/qLJBwPsJ6s4HDo6rt7qpXb4S\nESl2eTkDATCzRwEHfgYcAjwLHFXTWR5VbxgwBRgKVBH0h8x19+vzEqiIiISSzwGElwJtCC5HPQxc\n7O4LzKyPmW00s14A7v48cAcwC1gCLAJG5zFOEREJIW9nICIi0rjoFiYREclIUScQM+tlZs+Y2edm\ntsrMJphZwsGaZjYqUme9mU02swSrI5SusMfCzM43s+rIZcGax5BE71mqzGx/M3vRzL40s4/M7LQU\ndRv75yLUsWhsnwszG2lmr5vZNjN7MG7fCWb2oZltjhybpHfjmFkXM5thZpvMbImZZflW2dzL4rGo\nNLOtUZ+PBcnq1ijqBAL8D8G8WbsTjMw6jqAvJUbU9CdDgb7AXgTTnzQmoY5FxBx37xD1aDQLzJpZ\nc+Bp4K9AZ+DnwMNmtk+Cuo36c5HOsYhoTJ+LT4FbgAeiCyPD/p8AbiA4Jq8D01O8T/QUS+cC95rZ\ngFwEnEPZOhYO/CLq87F/fQ0XewI5AJju7l+5+2rguUhZvNrpT9z9S+Bm4Pz8hZkXYY8FQGOejm8/\nYHd3v8sDs4A5BNPixGvsn4t0jgU0os+Fu89w96eB+NnbTgfed/cnIhOwVgCDzKx//HtETbF0k7tv\ncfc5BAm5pFaez8axiJLWZ6TYE8jzwI/MrI2Z9QROAv6eoN4AgqlRatROf5KHGPMl7LFwYLCZrTGz\nf5vZjcku+zUizYCBCcqbwuciXrJj0Vg/F/F/8A4g6t88MvfexyQ+JsmmWEr2xazYNeRY1Lg98hl5\nxcyOq6/BYk8gFQS/7AZgOfBaJNPGSzX9SWNRQbhjMZtg0spuwA+Ac4Cr8hVkHvwb+MzMrjKzFmZ2\nIjCEYIh4vMb+uUjnWDTWz0X8MNJ21J02aQPBZyFeOlMslYKGHAsILvfuCXwLuA94xsxSThtatAnE\nzIzgW/dfCNYF6Qp0MbNxCapvAjpGbdesDdYoJltO51i4+yfuvjTy/H2CyzZn5DHcnHL3r4HTgO8S\nTPc/CngMWJGgeqP+XKRzLBrx5yL+W3f8vzkE/+6J/s3TqVsKGnIscPd57r7Z3b9294cILoeenKrB\nok0gBH8kDwV+H/mF1hHcoZ7oF2rs05+kcywSaTTXvgHc/T13L3f3ru5+ErA3MC9B1cb+uUjnWCTS\nGD4X8d+65xP8OwO1/Rx7R8rjpTPFUiloyLHISDEnkLUE36ouMbMyM9uFoFP0nQR1HwIuiAxp7Azc\nBDyYoF6pCn0szOwkM+sReb4fcCPwVHy9UmZmB5pZazNra8GknD0IEmq8xv65CH0sGtvnIvL/oDXB\njOJlZtYq0qczAxhoZqdH9o8G3nb3OpP6u/tm4Eng5sjx+zbwPWBq/n6ThsvGsTCzTmY2LPJZam5m\n5wLHEgzWSc7di/YBHAm8TDCN+xrgz0A3oA/BaVivqLqjCObOWk+wdG6LQsdfiGMB/DZyHDYRTANT\nAZQVOv4sH4s7gHWR33smsFekvCl+LkIdi8b2uYjE/03c49eRfScAC4AtwItAn6jXXQ/8LWq7M8Ef\n2k0EUyedXejfrRDHIvK3ZB5BH8kXwFzghPra1lQmIiKSkWK+hCUiIkVMCURERDKiBCIiIhlRAhER\nkYwogYiISEaUQEREJCNKICIikhElEBERyYgSiIiIZEQJREREMqIEIpIGM7vWzD42sw1mNr9mDfLI\nhHb/HVmMZ3FknepvzKxZZH8nC9ZkX2lmK8zslpp9IqWqeaEDECkxHwPfdvcqM/shwRrk/QjW5fgO\nwfTZW4DHiZ1eewrBZIZ7Eyzo8yzBwmD35S90kezSZIoiDWBmbxFMk30Z8Gd3vz9SfgLwAsGXtG7A\nUmAXd98W2X8OcKG7Dy1I4CJZoDMQkTSY2Y8JpojfI1LUnmDBr28RnFHUiF4VsC/QAlgVLC4JBJeP\nl+UyVpFcUwIRCcnM+hJcchoK/NPdPXIGYgQLfvWOqh79fDmwHdjV3b/JV7wiuaZOPJHw2hH0a6wF\nmpnZT4CBkX2PAZeZ2bciK0ZeE6mLu68C/j9wp5l1MLNmZra3mQ3J/68gkj1KICIhufsHwH8D/yTo\nEB8IvEKQKO4nSBLvAm8QrA5YHXXG8WOgJfABwQqCfwF2y2f8ItmmTnSRHDCzk4B73X2PQscikis6\nAxHJAjNrbWYnm1lzM+tJMDLryULHJZJLOgMRyQIzawO8BOwHbCW4z+Myd99U0MBEckgJREREMqJL\nWCIikhElEBERyYgSiIiIZEQJREREMqIEIiIiGVECERGRjPwfyPLgTvrv/tUAAAAASUVORK5CYII=\n", "text": [ "" ] } ], "prompt_number": 5 }, { "cell_type": "markdown", "metadata": {}, "source": [ "Just for a sanity check, let's look directly at the Dartmouth isochrone grid (the MASTERDF dataframe), query within the literature range of mass and logg, and show the resulting effective temperatures:" ] }, { "cell_type": "code", "collapsed": false, "input": [ "from isochrones.dartmouth import MASTERDF\n", "subdf = MASTERDF.query('feh > -0.31 and feh < -0.19 and M > 0.95 and M < 1 and logg > 4.42 and logg < 4.46')\n", "plt.hist(np.array(10**subdf['logTeff']));" ], "language": "python", "metadata": {}, "outputs": [ { "metadata": {}, "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAX0AAAEFCAYAAAAPCDf9AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAD85JREFUeJzt3X+sZGV9x/H3h92F1V1AdP2BP5YfGoJlW6wS0kaKTbeN\nhqYlgTbRqG1tLZiW+E/TkIgISq3a2qZVq4kENVo0kippbKupKdBUoTGkSgsuQaoIxV9gKXfZBYPw\n7R9zLg6zd/feuXPmx93n/UpO7s555pzzvc/O/cyZZ848k6pCktSGI+ZdgCRpdgx9SWqIoS9JDTH0\nJakhhr4kNcTQl6SGGPqS1JBDhn6SI5NcleSuJEtJvprkVV3biUkeT7J3aLlkNmVLktZj8xra7wbO\nrqq7k/wqcE2SXUP3Oab8hJckbQgZN6+T3AJcDnwV+Cawpaoe6780SVLfxhrTT/Js4BTgtqHV305y\nT5KPJHlGr9VJknq15tBPsgW4GvhYVd0B3AecAewEXgYc3bVLkhbUmoZ3khwBfBLYDpy70nBO9yrg\nu8DRVbVvaL3j/ZK0TlWVPve32hu5JAlwFfBM4Jw1jN97GagkLahVQx/4EHAq8MtV9aPllUnOBB4E\nvgEcB7wPuL6q9q60k416gc/gOW9j1r+Rawfrnzfrn6/l+vt2yNBPcgJwAfAI8L2hIi4EHgf+FHgW\nsAT8M/CaqVQpSerF2Jdsjn2Abkx/oz/bbsT6N3LtYP3zZv3zNVR/r6f8hr4kLaBphb5vukpSQwx9\nSWqIoS9JDTH0Jakhhr4kNcTQl6SGGPqS1BBDX5IaYuhLUkMMfUlqiKEvSQ0x9CWpIWuZT1/q3bTm\nCl8rJwBUqwx9zdG8gne+TzjSPDm8I0kNMfQlqSGGviQ1xNCXpIYY+pLUEENfkhpi6EtSQwx9SWqI\noS9JDTH0Jakhhr4kNcTQl6SGGPqS1BBDX5IaYuhLUkMMfUlqiKEvSQ0x9CWpIYa+JDXE0Jekhhwy\n9JMcmeSqJHclWUry1SSvGmrfneT2JPuSXJdk5/RLliSt12pn+puBu4Gzq+oY4K3ANUl2JtkBfBa4\nBDgOuBn49DSLlSRNJlU13gbJLcDbgR3Ab1XVWd36pwL3Ay+pqjuG7l8A4x5Hh7ckwLweE/HxqIU3\n+BuBqkqf+x1rTD/Js4FTgFuB04Bbltuqaj9wJ7CrzwIlSf1Zc+gn2QJcDXysO5PfBiyN3G0J2H6Q\n7Z9YJEkHmkVOrin0kxwBfAJ4BLioW/0QcMzIXY8F9vZWnSSpV6uGfgZPOVcBzwTOr6rHuqbbgNOH\n7rcNeGG3/gBV9cQiSTrQLHJyLWf6HwJOBX69qn40tP5aYFeS85JsBS4Dvjb8Jq4kabEc8uqdJCcA\n32IwrPPYUNMFVfWpJLuBDwAnAP8O/E5V3T2yD6/e0QG8ekc6tGldvTP2JZtjH8DQ1woMfenQFuKS\nTUnSxmboS1JDDH1JaoihL0kNMfQlqSGGviQ1xNCXpIYY+pLUEENfkhpi6EtSQwx9SWqIoS9JDTH0\nJakhhr4kNcTQl6SGGPqS1BBDX5IaYuhLUkMMfUlqiKEvSQ0x9CWpIYa+JDXE0Jekhhj6ktQQQ1+S\nGmLoS1JDDH1JaoihL0kNMfQlqSGGviQ1xNCXpIYY+pLUEENfkhpi6EtSQ1YN/SQXJbk5ySNJPjq0\n/sQkjyfZO7RcMt1yJUmT2LyG+9wLXAG8EnjKCu3HVFX1WpUkaSpWDf2quhYgyRnA81e4yxHAYz3X\nJUmagnHG9HOQ9d9Ock+SjyR5Rh9FSZKmY5zQHx3CuQ84A9gJvAw4Gri6p7okSVOwljH9ZU8606+q\nfcB/dDd/kOQi4LtJtnVtT944Gd52HaVK0uFtOCenZZIz/T72KUmaoVXP9JNsArZ0992U5CgGb9y+\nFHgQ+AZwHPA+4Pqq2rvSfjy7l6RDG87JaZ31r+Ws/FJgP3Ax8DrgYeAtwMnA54El4L+69a+ZSpWS\npF5k2mfgSQo809eTDc5i5vWYiI9HLbzlM/2q6vWU3/F3SWqIoS9JDTH0Jakhhr4kNcTQl6SGGPqS\n1BBDX5IaYuhLUkMMfUlqiKEvSQ0x9CWpIYa+JDXE0Jekhhj6ktSQcb4uUTpszOJr6VbilM6aN0Nf\njZpH+M7niUYa5vCOJDXE0Jekhhj6ktQQQ1+SGmLoS1JDDH1JaoihL0kN8Tr9BTGvDwvp8Dfvx5Yf\nSFsshv5CmdcHhvyg0uFvXsHr//OicXhHkhpi6EtSQwx9SWqIoS9JDTH0Jakhhr4kNcTQl6SGGPqS\n1BBDX5IaYuhLUkNWDf0kFyW5OckjST460rY7ye1J9iW5LsnO6ZUqSZrUWs707wWuAD4yvDLJDuAz\nwCXAccDNwKf7LlCS1J9VJ1yrqmsBkpwBPH+o6Tzg1qr6TNd+OXB/klOq6o4p1CpJmtA4Y/qj0+Wd\nBtyyfKOq9gN3Art6qEuSNAXjhP7o3KzbgKWRdUvA9okqkiRNzTjz6Y+e6T8EHDOy7lhg74obD32R\ng1+qIEkHmsUX3owT+qNJfRvw28s3kmwDXtitlyRgft/c5cnlytZyyeamJFsZPEFsSnJUkk3AtcCu\nJOd17ZcBXzvYm7hV9cQiqSU1h2VjmkVOrmVM/1JgP3Ax8DrgYeCSqrofOB94J/C/wBnAq6dUpySp\nB5n2mXeSAl9qrWbwEri178id5/e2zud3nsffwfweW9BaX/dpeVisqnodH3MaBklqiKEvSQ0x9CWp\nIYa+JDXE0Jekhhj6ktQQQ1+SGmLoS1JDDH1JaoihL0kNMfQlqSGGviQ1ZJz59GfiwQcf5IEHHpjL\nsTdt2sQLXvCCuRxbkmZh4UL/yiuv5C1veQdHHvn0mR636sfs33/vTI8pSbO2cKEPUHUB+/a9d8ZH\n/Q7wPOY7Ba0kTZdj+pLUEENfkhpi6EtSQwx9SWqIoS9JDTH0Jakhhr4kNWQhr9OXDleJn8fQfBn6\n0kzN48N/PtHoJxzekaSGGPqS1BBDX5IaYuhLUkMMfUlqiKEvSQ0x9CWpIYa+JDXE0Jekhhj6ktQQ\nQ1+SGjJx6Ce5IcnDSfZ2y54+CpMk9a+PM/0C/rCqju6WF/ewT0nSFPQ1vOM0fpK0AfQV+u9Kcl+S\nLyV5RU/7lCT1rI/Qvxg4CXgu8GHgc0lOHr1TkicWSdKBZpGTE4d+VX2lqvZV1aNV9XHgy8A5k5cm\nSes3HKCzXhbZzL45q2oe3xgkqV3zypz1h/5wTk7ryWOiM/0kxyZ5ZZKtSTYneS3wC8AX+ilPktSn\nSc/0twBXAKcCjwF7gHOr6s5JC5Mk9W+i0K+q+4Eze6pFkjRlTsMgSQ0x9CWpIYa+JDXE0Jekhhj6\nktQQQ1+SGmLoS1JDDH1JaoihL0kNMfQlqSGGviQ1xNCXpIYY+pLUEENfkhpi6EtSQwx9SWqIoS9J\nDTH0Jakhhr4kNcTQl6SGGPqS1BBDX5IaYuhLUkMMfUlqiKEvSQ0x9CWpIYa+JDXE0Jekhhj6ktQQ\nQ1+SGmLoS1JDDH1JaoihL0kNMfQlqSETh36Spye5NslDSe5K8po+CpMk9W9zD/v4G+AR4FnAzwL/\nmOSWqvp6D/uWJPVoojP9JNuA84BLq2p/VX0Z+Hvg9X0UtxjSLRvRRq4drH/erP9wNOnwzinAj6vq\nzqF1twCnTbhfSdIUTBr624GlkXV7gaMn3K8kaQomHdN/CDhmZN2xDIL/SZJxX2b9xXprmtDB6pzF\ny8RpHWO1/c7rJfBajzuN+mb5O+cg/56lSY47ac3zHmKZ/fHHz7vZmfRM/w5gc5IXDa07Hbh1wv1K\nkqYgVTXZDpJPAQW8EXgp8A/Az1fVnsnLkyT1qY8PZ/0B8BTgB8DfAm8y8CVpMU18pi9J2jichkGS\nGjJW6Ce5IcnDSfZ2y56htjcm+Ua3/vNJjh/Z9j1J7u+Wd4+0nZjk+iT7kuxJsnuyX6vf+pNcnuTR\noe2Wkpw46/q7Y726O8ZDSe5Mcla3fneS27sarkuyc2S7uff/eutf8P5/eZItSf4uybeSPJ7kFSts\nN/f+X0/tC973ZyX5uSRfTPLDJD9Ick2S54xsN/e+X2/9U+n/qlrzAlwP/O4K638R+D7wYmAL8EHg\nhqH2C4Hbged2y23AhUPtNwHvBY5i8AnfB4Ad49Q25fovAz5+iP3Oqv5fAe4CzuxuH9/15w7gQeB8\n4Ejgz4CbFrD/11v/ovf/FuDNwMuB7wBnj2w39/6foPZF7/tXdY+b7QzeW7wK+Pwi9f2E9ffe/+MW\nfj3weyusfy/wgaHbxwOPAyd1t28E3jjU/oblP2oGn+p9BNg21P6vw/8xPXb8euu/HPjEQfY5y/pv\nBN6wwvoLgC8N3X4qsB84ZcH6f731L3T/j9znHg4Mzrn3/wS1b5i+7+73UmBpkfp+wvp77//1jOm/\nK8l9Sb409FKwePInIJb3u6v7+VMMpmdY9p/8ZKqG04BvVtW+ofZpTuWwnvoL+LXuJditSd40dN+Z\n1J9kE/Ay4FkZDEPdk+T9SbZ2x3qif6tqP3DnUA1z7/8J61/0/l/NXPt/wto3Wt+fzZM/J7Toj/1R\no/X33v/jhv7FwEkMXpZ8GPhckpOBLwC/meSnkzwFeFtX7FO77bYzePm+bKlbt1Lbcvs0pnJYb/3X\nAKcyGIb4feBtSV494/qfzeCl+PnAWcBLGMxq+lZgGwdOhzFcwyL0/yT1L3r/r2be/T9J7Rum75P8\nDHAp8MdDq+fd9zBZ/b33/1ihX1Vfqap9VfVoVX0c+DJwTlX9C4OXIZ8BvtUte4H/6TYdna7h2G7d\nSm0AT+PAEJjYeuuvqj1V9b0auAn4a+A3Zlz/w93P91fV96vqh8BfAuccpIbh6TAWof/XXf8G6P/V\nzLv/1137Run7DGYF+CfgzTWY7XfZvPseJqh/Gv3f2yWbVfXBqjqlqp4DfJbBvD7LL1NuY/Dstuz0\nkbaTk2wfab+tr9rWYpX6D2Um9VfVA/zkSXSlGk5fvpHBlNcvHKph7v0/Yf2Hsgj1r2au/T9h7Yey\nEH2f5ATgi8A7qurqFWpc5Mf+avUfyvrqH+ONiGOBVwJbGQTiaxk807yIwTvHuxiMi+8EbgD+ZGjb\nC4GvMxhWeV5X1AVD7TcBf97te/kd6Ges902TKdR/LnBc134mcC/w+lnW3x3n7cBXgGd29fxbt24H\n8H/dsbcyuPrlxkXq/wnrX+j+79qO6o5/D4MrNbYuUv9PUPtC933Xn/8N/NFBtpt7309Yf+/9P07R\nO7qil7od3wjs7tqexuANhIeA7wLvpPu079D27wF+2C3vHmk7gcGVNfuBPcAvTaHT110/8EngfgbD\nDXuAi2Zdf3eczQy+qeyBrs6/Ao7s2nZ3x94PXAfsXKT+n6T+DdL/dzG44uuxoZ/Dv8O8H//rqn3B\n+/4oBpc0Pt7Vt7wsjWy7qI/9VeufRv87DYMkNcRpGCSpIYa+JDXE0Jekhhj6ktQQQ1+SGmLoS1JD\nDH1JaoihL0kNMfQlqSH/D6G8az34av6PAAAAAElFTkSuQmCC\n", "text": [ "" ] } ], "prompt_number": 5 }, { "cell_type": "markdown", "metadata": {}, "source": [ "Nowhere do we see anything close to 5640 K. What is going on here?" ] }, { "cell_type": "code", "collapsed": false, "input": [ "mod.save_hdf('kep22.h5')" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 6 }, { "cell_type": "code", "collapsed": false, "input": [ "mod2 = StarModel.load_hdf('kep22.h5')" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 7 }, { "cell_type": "code", "collapsed": false, "input": [ "mod2.plot_samples('radius')" ], "language": "python", "metadata": {}, "outputs": [ { "metadata": {}, "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAYwAAAEXCAYAAAC+mHPKAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XucFNWZ//HPw4DADIooggqCi0YlsEEUBYm6WUhEIS75\nJSZGjZpdo4vGxJ9xk6i4LMaN16yyEYOaqFH05w+vK8FLslGIIooiBrJIVEQuDg6CIJeB4TI8+0d1\nT1/mVjN0d1V3f9+vV72Yc6q6+umipp85darOMXdHRESkNR2iDkBERIqDEoaIiISihCEiIqEoYYiI\nSChKGCIiEooShoiIhKKEISIioRQkYZjZPmZ2n5mtMLPNZva2mZ2etn60mf3VzGrN7CUz61eIuERE\nJLxCtTA6AquAU919P+A64DEz62dmPYGngIlAD2ABMKNAcYmISEgW1ZPeZrYIuB7oCVzg7icn6iuB\n9cCx7v5eJMGJiEgjkfRhmFlv4Cjgf4BBwKLkOnffBiwDBkcRm4iINK3gCcPMOgGPAL9NtCCqgM1Z\nm20GuhU6NhERaV7HQr6ZmXUApgN1wOWJ6q3Aflmbdge2ZL1WoySKiLSTu9ve7qNgLQwzM+A+4CDg\nG+5en1i1BBiStl0VcESiXkREYqJgnd5mdjdBYviyu9em1fck6LP4J+A54GfAye4+Muv1DqDh2CHI\nvToWoGORTsciRcciJe1Y7HULoyAJw8z6Ax8SXIqqT1t1ibs/amajgalAf+B14LvuviprH0oYCfpl\nSNGxSNGxSNGxSCm6hJELShgp+mVI0bFI0bFI0bFIyWXCKGint+SGfglSdCxSdCxSdCzyQ2NJiYhI\nKEoYIiISihKGiIiEooQhIiKhKGGIiEgoShgiIhKKEoaIiISihCEiIqHowT0pevPmwTe+0fz6//5v\nGKzZVUT2mhKGFL2dO6Gmpvn1u3cXLhaRUqZLUiIiEopaGFJSRoyAp56CL38Z3nknVVdRkdrmO9+B\ne+6BUaOgtrbp/Vx8MXzve/mPV6SYKGFISencGQ45BDp1StXt2JG5zc6dwb9vvQWbsycHThg7Nj/x\niRQzXZKSsrNhA/zlL1Bf3/q2IpKihCEl6bXXYOvW1DJtWmrdzJnwhS9kXo764x/hn/6p8HGKFBMl\nDClJXbtCVVVq6dy55e2HDYO+fQsTm0ixUh+GlIUePVp+FqOD/nQSaZUShpSFr30tWESk/fR3lYiI\nhKIWhhSdOXNg9epUeenSyEIRKStKGFJ0pkyBZ56JOgopFwsWLGD69OkMGzaMV199lR//+MccccQR\nbd52/vz5vPrqq2zevJl58+Zx3XXXceqppxbyo+w1JQwRKWvTpk3j0ksvbXLdjh07OOuss5g/fz69\ne/dm4MCBnHPOObzxxhtt2nbbtm3813/9FzfddBMATzzxBGeccQbvv/8+hx56aF4/Xy6pD0OK2qmn\nBkN9JJevfCXqiKTYfPLJJ82ue/nll+nWrRu9e/cG4Pjjj2fp0qWsWLGiTdsuW7aMW265heXLlwNw\n2mmnsX37dubNm5f7D5RHamFIUbvySt39JPmzYsUKDjzwwIaymdGjRw+WLFnC4YcfHnrbcePGMW/e\nPAYMGADARx99BMDnPve5/H+IHFLCEJG8effdd7n22mvp378/Zsb69eu57bbb6NWrV6uvXb58Oddf\nfz1du3aloqKC2tpabrnlloa/4LNt27aNESNGsHjx4pzFv379eiorKzPqunTpwpYtW9q87YgRIxrq\nb7rpJq666iqGDBmSs1gLQQlDRPJi06ZNjB49mltvvZVzzz0XCL4oR48ezcKFC+mUPkJklg8//JAT\nTzyR3/zmN3wt0YR85JFHGDNmDG+99RYV6cMPA2+++SaXXXYZS5YsaTGmTz75hClTpuDuDXVz586l\nrq6uodytWzcmTpwIwP7775+xLcDWrVvp2bNno32H3fa+++6jT58+3HzzzS3GGkvuXhQL4EG4Uu7G\nj3eHYHn66dztd9Kk1H4nT87dfsvVtdde67179/b6+vqGug0bNninTp182rRpLb52/PjxfuCBB2bU\nbdu2zbt06eK//e1vG+qWLl3q48aN8+9+97s+YsQI79ChQ5vjnNzCf/ZLL73kQ4YMaSjv2rXLO3fu\n7O+99167tp01a5bff//97u5eV1fnK1asaHO8bZX23bnX38Pq9BaRvHj88ccZPnw4HdLGXenRowfH\nHHMMjz/+eLOv27lzJ88++yxHHnlkRn3Xrl3p378/TzzxREPdMcccw6xZs3jggQc4+uijG/2Fv7dO\nOeUU1q1bx+rEgz9/+tOfGDRoUEPfw4svvsiiRYtCbfunP/2JtWvXMnbsWGpqanjuuef4+OOPcxpv\nvumSlEgTamog8T0ABBMwaV7w8LZs2cKyZcsYM2ZMo3WHHHII8+fPb/a1GzZsoL6+ni5dujRat99+\n+zV5S2u+dOzYkenTp/Pzn/+ckSNHMnv2bGbMmNGw/q677uK4445jyJAhLW67fPlyzjzzTLZu3drw\nWjNj06ZNBfssuaCEIdKEu+8OlqQePYJ5NCSclStXAsEXfLaqqio2b97Mrl27muzHOOigg6isrMzo\nV0has2YN69evZ8+ePRktl73RtWvXFtePGjWKUaNGAXDBBRdkrHvqqadCbTtgwAA2NzdbVxHRJSkR\nybnkl+M+++zTaF1VVRUAn332WZOvraio4Oyzz+aDDz7IqK+pqWHNmjVA0ArJlZ/85Cc521epU8IQ\nSejdO5hYKbl8/vNRR1S8kncxmVmjdbt27QKgvoUpD2+55RYOOOAA7r33XgB2797NXXfdxdChQ3H3\nRndJSWEoYYgkXHZZ0G+RXF55JeqIitdBBx3U7Lra2lrMjH333bfZbXr27Mn8+fOpqanhqquu4vrr\nr+fiiy/GzOjatSs9evTIR9jSCvVhiEjO9e7dGzNj48aNjdbV1tay//77N1yaas7+++/PpEmTMurW\nrl3LyJEjcxqrhKcWhojkXFVVFUOHDmXVqlWN1i1btoxjjz22zftct24d1dXVfPOb38xFiM1asGAB\nV1xxBdOnT2fChAmN+lLCbjt//nxuv/12Jk+ezGmnncbLL7+c17gLIhcPcxRiQQ/uSUK+HtzL9umn\nqffp0SN/71OqJk2a5H369MmoW7ZsmZuZ/+pXv8qof++993z79u0N5TvuuMN79erl1dXVDXW33Xab\nH3bYYb5jx44m3+/CCy9s14N76erq6rx///5eU1Pj7u5vvvmmn3DCCW3etra21q+++uqGbR9//HGv\nrKzM+DyFgh7cE5G4u/TSS6mtreXhhx9uqLvzzjsZNGgQF198cUPdyy+/zNFHH90wBAgEl62qqqoa\nOrfffvttfvGLX/DYY481eecVQF1dHe7O9u3b2x2zRqdtmfowRCQvDj74YObMmcPEiRNZuHAhW7Zs\nYePGjbzwwgt07NgxY7tevXplPNl91VVXsWnTJq655hq2bNnCpk2beOaZZxg+fHjGe6xbt47zzz+f\n6upqlixZgpnRv39/Bg8ezPe+972GMazC0ui0LVPCEJG8GTJkCLNmzWpxm6OOOoqampqMui5dunDr\nrbe2uv+DDjqIF154Ya9iTKfRaVumhCEiJW/q1Kktdl4PGzaM8847j+7duzcaj6qsR6fNooQhIiXv\n8ssvD7XdwIEDGx4WhOCBwQ0bNtC/f/9G2x5zzDHcc889LW777LPP0qFDB26++WZ27NhBTU1Nk/sq\nFur0lljbswdqazOX3bujjkpKlUanbZlaGBJrq1dDVl+jSN5odNqWWfY1uLgys+BhjCKJV3Jj5cqW\nE8bTT+dvTu8NGyB5E0zHjnD66Znrr7wSEgOTisRWcjwvd288sFcbqYUhRSXrphQKNQbd7t2QfbNP\nnh84FomdgvVhmNnlZrbAzOrM7IG0+sPNbI+ZbUlbJhYqLike/fo17s8488yooxIpH4VsYVQDNwBj\ngKZmLNnPdb1JYqRbN5g5M7Pu9tthzpxIwhGJXMEShrs/DWBmw4C+TWzSAWh+gHyRAttnn8YtmBam\nohYpeVHcVttcx8tKM1ttZveb2YHNbCMiIhGJImFkX3ZaBwwD+gHHA/sCjzT3YjNrWEREpLF8fU9G\ncZdUxidw91pgYaL4iZldDnxsZlWJdSIiEgNRJIywHdtNtn7ULy4i0rL078lctjIKljDMrALolHjP\nCjPrTNDJfRywCXgf6AH8Epjt7o2HhxQRkcgUsoXxr0D6BL3fASYD7wE3Ar2AzcAfgHMKGJdIu7zz\nDsyenSp36wYnnBBdPCL5FmpoEDP7F3f/RRP1P3L32/MSWeP30tAgZeCDD6C6OlWuqYGzzw5+7tcv\nGCokShdcANOnN71u6FBYuLDpdSJRiWJokH8DGiUMglZDQRKGlIepU2HKlKijEJGmtJgwzGwUwV1N\nFYmf0x1BcAlJpGwMHAhf+lKqvGULvPVWZOGIFFRrLYz7Ce5q6gzcl1bvwFrgB3mKS4S/+RvomzYm\nQO/e0cWSdM01wZK0cCEcf3x08YgUUosJw90PBzCz6e5+fkEiEkn4wQ+CIcRFJB5C9WGkJwsz65C1\nbk+ugxIRkfgJNTSImR1vZq+Z2TZgd9qyK5/BiYhIfIS9S+pBYCZwEbAtf+GIiEhchU0Y/YCJmq9C\nRKR8hR2t9mmCiY9ERKRMhW1hdAWeNrNXCG6nTXJ3vyD3YYmISNyETRjvJJZsukQlIlImwt5WOznP\ncYiISMyFShhNDAvSwN1fyl04IiISV2EvSSWHCEk6iGC4kNXAgFwHJSIi8RP2ktTh6eXEZEjXAVvz\nEJOUkcsuyyzPnRtNHCLSunZNoOTu9WZ2I/AR8B+5DUnKybRpUUcgImGFfQ6jKV8hmGJVRETKQNhO\n79VZVZVAF+CyJjYXaZepUyF9vvovfjG6WESksbCXpLKHNq8F3nP3TTmOR8rYZZdlJgwRiZewnd5z\noGFo897AWg1rLiJSXsIOb76fmT0E1AHVQJ2ZPWRm3fManYiIxEbYTu87gSpgMEH/RfLfO/MUl4iI\nxEzYPozTgQHuXpsov2dm3wWW5yUqERGJnbAtjO0ET3en60lwiUpERMpA2BbGb4D/NrP/AFYChwNX\nAr/OU1wiIhIzYRPGjcAa4DzgkMTPtxCMMSUiwF/+An37ZtY9+SQMHx5NPCK5Fva22j0EyUEJQqQZ\nu3dDdXVm3c6d0cQikg9hb6u908xGZtWNNLMp+QlLRETiJmyn9znAW1l1CwkuUYmUrcGDYfXqzOWE\nE6KOSiQ/wvZh7KFxcukAaCAHKWv77NO436Jz52hiEcm3sC2MucC/J4YGSc6HcT3wSr4CExGReAnb\nwrgCmAXUmNlKoB/wMXBmvgITEZF4CXuX1GozOw44ETiMYGrW+RqAUNpixw6o06OeIkUr9Ix77l4P\nvJZYRNrsnnvgiiuijkJE2mtvZtwTEZEy0q45vUX2VqdO0LVr1FGISFuohSGRmDABNm3KXDTbnki8\nNdvCSN5C2xp1fIuIlIeWLkntzio7jR/Uc6AipxGJiEgstZQwBqT9PA44i2DU2lUEz2FcDTyZv9BE\nRCROmk0Y7r4i+bOZ/QgY5u4bE1XvmtkCYAHwq7xGKFLEXnkFNmxIlQ85BE48Mbp4RPZG2Luk9iOY\nw3tjWl0l0D3nEYmUkIkTM8tf/3owR4ZIMQqbMB4E/mhmdxA85d0P+CHwUL4CExGReAmbMH4CLAO+\nTTDj3sfAnWiKVpFGTj4ZDjggVV6zBhYsiC4ekVxpy4x7dycWEWnBTTdllp98Es46K5pYRHIp7Ix7\nHczsEjN7ycz+kqg71cy+ld/wREQkLsI+6X09cBHBJah+ibpqgltrQzGzy81sgZnVmdkDWetGm9lf\nzaw2kZT6NbcfERGJRtiE8Y/AV939UYLZ9wA+JPNZjdZUAzcA96dXmllPguc5JgI9CG7VndGG/YqI\nSAGE7fTuAGzNqqsCtoR9I3d/GsDMhgHpk1p+Hfgfd38ysX4ysN7MjnL398LuX0RE8itsC+N54HYz\n6wIN40zdAPyuHe+ZPbzIIGBRsuDu2wjuyBrcjn2LiEiehE0YPwIOBj4jeIhvK3A4bejDSONZ5Spg\nc1bdZqBbUy82s4ZF4q2uDrZvTy07d0YdkUh5yNf3ZNjbajcB/8fMegP9gdXu/nE73zP7E2wlSELp\nutOGy10ST5//PHz4YdRRiEiuhL2t9gEzG+Pua939jWSyMLP2jCOV3cJYAgxJe68q4IhEfeMXuzcs\nIiLSWL6+J8NekvoOcL+Z/Tir/vywb2RmFYk+kI5AhZl1NrMK4GlgsJl9PbH+34A/q8O7tOyzD3Tp\nklo6dYo6IhFpq7B3SdUBw4FnzGwIcJG772jje/0rMCmt/B1gsrv/zMy+AUwFHgZeJxiCRErI0qUw\noC03YYtI7ISeotXdPwJOTrzmFTPr28pLsl8/2d07ZC0/S6x70d0Hunulu49y91Vt+hQiIpJ3bZrT\n2923u/u5wFPAfGCfvEQlIiKxE/aS1A3pBXe/2cwWARpLSkSkTIS9rfbWJuqeJ3igT0REykCzCcPM\nfu/uYxI/v9LMZu7up+YlMhERiZWWWhjps+nd18w2ehhCRKRMNJsw3P2RtJ9/W5BoREQktlq6JHUR\nIVoQ7n5/a9uIiEjxa+mS1PmEu+SkhCEiUgZauiT1pQLGISIiMRf2OYwGFoyX2zDirLvvaWFzEUnz\n3HPQL2sC4tdfh0MPjSYekbYIlTDMrA/BWE9/RzD0eDJhOFCRn9BESk9dHaxenVlXXx9NLCJtFXZo\nkLuBXcAogvkrjgOeAS7NU1wiIhIzYS9JfRHo5+5bzQx3/3PiLqp5wL35C0+Kxbp1sHhxZt327dHE\nEjdnnAErV2bWDR8ONTXRxCPSXmETxu7EArDRzHoBm4A+eYlKis5rr8H48VFHEU+VlY37LTq2ufdQ\nJHphL0m9AZyR+Pn3wAyCiY8W5CMoERGJn7B/53yHVHK5ErgK6AZMyUdQUtwOPBCGDMms69IlmlhE\nJHfCjlb7WdrP28ga7lwk3ciRMHNm1FGISK6Fva22E3AOMJSgZZHk7n5JPgITEZF4CXtJajrwtwTz\nXyTv7TA0Wq2ISNkImzBOJ7itdnM+gxERkfgKe5fUUuCAfAYiIiLx1pa7pO43s+eBtYk6I+jDeKj5\nl4mISKkImzAuJHjaez8g+/ldJQwRkTIQNmH8X2Cou7+Tz2BERCS+wvZhrAVW5TMQERGJt7AtjNuB\n6WZ2K6k+DADcfXnOoxIRkdgJmzDuSvybPbyc5sMQ2UtnngmdO6fK550HP/xhdPGINKfVhJGYYe9z\nwCp335X/kETKy6JFmeVTTokmDpHWhO3DWAxoXjARkTLWagvD3d3M/gwcTfAAn4jspZkzYefOVPnh\nh2Hq1OjiEQkjbB/GbOB5M/stsJqg7yL54N79eYpNpGQNHZpZfuWVaOIQaYuwCeNkYAXwd02sU8Io\nQ7NmwZo1qXL2dXgRKT1h58P4Up7jkCJzxx3w0ktRRyEihRR6ZmEz6wH8A3AoUA3McvcN+QpMRETi\nJewESicBzwJ/BVYCZwJTzOyr7j4vj/FJERg7Fvr0SZW/8IXoYhGR/AnbwvhP4DJ3///JCjM7O1F/\nQj4Ck+Lxox/B6NFRRyEi+Rb2OYyjgMey6p4keKBPRETKQNiE8T7BnN7pvgksy204IiISV2EvSV0B\nPGtmPyAYtbY/Qavjq/kKTERE4iXsbbXzzOwIYBzBXVIzgefd/dN8BiciIvER+rbaxC200/MYi4gA\n8+fDzTenymbw059GF49IUosJw8xmt/J6d/dROYxHpOzNnRssSRUVShgSD621MB5pos6BPgT9GpU5\nj0hip6YGpkzJrFum2x1Eyk6LCcPdf5NeNrOewNXAxcAM4Gf5C03iYt06uOWWqKMobSeemNmK2LMH\nbrstunhEmhL2Se/uwL8APwBmAce5+wf5DEyknJx6arAk7d6thCHx01ofRiXBpad/AeYAX3T3JQWI\nS2KqVy+48srMuiOOiCYWESms1loYHxI83HcrsADobWa90zdw95yMWWpmc4DhwO5E1UfuPjAX+5bc\n6dULrr466ihEJAqtJYztiX8ntLDN3+QoFge+rwmZRETiqbVO78MLFEeSFfj9REQkpLBjSRXKTWa2\nzszmmllTs/uJlJ36eqiszFzuvjvqqKQchX7SuwB+CiwBdhIMdPg7MzvW3Zenb2SWaoS4e0EDFInK\n9u2Z5d27m95OBDK/J3MpNi0Md3/D3WvdfZe7PwS8CoyNOq5yVF8ffCGlLyIisUkYYbl7wyL5cckl\n0KlTajnuuKgjKj8VFVBbm7lccknUUUmxyNf3ZCwShpl1N7MxZtbFzDqa2XnAKcALUccmEgWzxv0W\nnTpFHZWUu7j0YXQCbgCOAeqBpcB4d9eIRREzgw5pf1ZUVEQXi4hEKxYJw93XAydGHYc09utfw0UX\nRR2FiMRBLC5JiYhI/ClhiIhIKEoYIiISihKGiIiEooQhIiKhKGGIiEgoShgiIhKKEoaIiISihCEi\nIqHE4klvEWmbadPg+edT5UMPDZ7KF8knJYwyV18Pixdn1n36aTSxSHjvvBMsSUceGV0sUj6UMMrc\n9u0avlxEwlHCECkSEybAGWekytXV8M//HF08Un6UMKSBGQwZkll3wAHRxCKNDR4cLEnvvx9dLFKe\nlDCkQWUlvP121FGISFwpYYiUgK1bYebMzLrhw6F372jikdKkhCFSAmpqYPz4zLpZs2DcuGjikdKk\nB/dERCQUtTBEilRVFZx5Zmbd66/DunXRxCOlTwlDpEgdemjjfotx4+C556KJR0qfLkmJiEgoShgi\nIhKKEoaIiISihCEiIqGo01ukRE2aBFOnpspDh8KNN0YXjxQ/JQyRErVwYWZ59+5o4pDSoYRRZh54\nAO69N1Wur48uFhEpLkoYZaa6Oni4S0rT9dfD5ZenygsXwnXXRRePlBYlDJESMmxYZrmjfsMlh3Q6\nlbELL8ycgKeD7pkTkRYoYZSxvn3hpJOijkJEioX+phQRkVDUwihx3/9+MLlO0qJF0cUi0XrxReja\nNbNu4UIYODCaeKT4KGGUuBkz4NNPo45C4sAd6uoa14mEpUtSIiISiloYZeTOO2HffVPlL3whulik\nMEaNgm3bMuuGDoV3340mHiluShhl5Nvfhp49o45CCqmionG/hW6flvbSqSNSxgYNArPUcsklUUck\ncaYWRgmZOxd+9avMui1boolFREqPEkYJWbkSHn006ihEpFQpYYiUmSVLMsv33gsTJgQ/P/88jB2b\nuf6pp6BLl8LEJvGmhFHExo+HxYtT5fQH9EaODB7aS9etW2Hikngza7780UfBkm7sWKisTJUvuAC+\n9a38xSfxpYRRJNasgccey6ybObP57Q8/HM49N68hSZmYPTuz/MUvRhOHRE8Jo0h8+CFceWXUUUgp\nOv10ePbZzLpx45rf/g9/yHy2o6oKrr46P7FJvJgXydgAZuYAxRJvrr36Kpx8cvPrZ8yAE05Ilauq\noFev/Mclpem552DPnlT5wQfhiSea3rZXL1i7tjBxSdtZ4pqju1srm7YqNi0MMzsAuA/4CrAeuMbd\ndc9PEw4+GM4+O7NuxAjo1y+aeKT0ZHd8v/pq89t+8knwRHm6iRNh9OjcxyXRik3CAO4C6oBewFDg\nWTNb5O7vRBtW/AwYAFOmRB2FlJOvfCWz43vtWrjrrlQ5u59j0KDMJ8y7ddNQNKUgFpekzKwK2AAM\ncvdliboHgTXufk2iXLKXpD79tOkO6t27YdUquPba4BJTqlXR0MQsWIxxldbcjjiS6BXyWKxdG7R0\nwxo+vOm55LduDVoozenXr33TzOq8SMnlJam4JIyhwFx3r0qr+xHwJXf/h0S5aBNG9+6wc2fT637+\n82CMpz592rJH/TIk6YshpZDHYufOYGSBdOec0/KXf7Jjfdu24FmQk06C3/8ebr+9+dfU1ATPgHz2\nWdPrzYKk8stfwgsvpOqffz44Fmec4Rx2GNxzT4gPVaJKMWGcAjzm7oek1V0MnOvuf58oF23C6Ny5\n+YSRXL9jR1v2qC/JJCWMlKiPxaRJ8Mc/psqvv16Y+TaqqqC2Nrs2+d3oDBwI75Txhe1S7PTeCuyX\nVdcdaDQSkmU/dVQC2pYsUkrxWLSXjkVKuR2LxskinbF0aeOHFaV94jJa7XtARzM7Mq1uCPA/EcUj\nIiJZYnFJCsDMHgUc+B5wHDALOMndl0YamIiIAPFpYQBcBnQFPgEeBiYoWYiIxEdsWhgiIhJvcWph\niIhIjEWaMMzsADN72sy2mtkKMzunhW3/1cxWm9lnZjbbzD7fnv3EVQ6PxRwz225mWxJL0V3WC3ss\nzKyzmd1hZtVmtsHM7jKzjm3dT5zl8FgU9XlhZpeb2QIzqzOzB1rZ9koz+9jMNpnZfWa2T9q6Ujgn\ncnUs2nxORN3CSB8O5DxgWvqXX5KZ/QMwATgFOAB4DZje1v3EXK6OhQPfd/d9E8vAvEeee2H/P68m\nuEFiEHBU4ufr2rGfOMvVsSj286IauAG4v6WNzGwM8FNgFNAfGABcn7ZJKZwTuToWbT8n3D2SBagC\ndgBHptU9CNzUxLbXADPSyoOA7W3dT1yXXB2LRHk2cFHUn6lAx+JN4Ky08jnAqjI9L5o9FqVwXqR9\njhuAB1pY//+Af08r/z3wcamcE7k6Fu09J6JsYRwF7PbE2FEJiwi+ALO9CJxkZp8zs07AhcDz7dhP\nXOXqWCTdZGbrzGyumf1dfkLOm7b+f6Y/ktUB6Gtm+7ZjP3GUq2ORVMznRVJrj+B9nuAYJS0GeptZ\nD0rjnEi3N8ciqU3nRJRPencDNmfVbQH2zd7Q3d+wYDDCd4F6YBWQHDw59H5iLFfHAoIm6BJgJ8Ff\nmb8zs2PdfXk+As+Dtvx/vgBcYWazCc7lHxI0syvbuJ+4ytWx2ELxnxdJrd3W2Q3YlFZOHr99KY1z\nIt3eHIuNtOOciLKF0ZbhQC4n+FLsC3QGfga8ZGZd27KfGMvVscDd33D3Wnff5e4PAa8CY7P3E2Nt\n+f/8OfDh210zAAADrUlEQVQ28GdgLvA0wV+Qa9u4n7jK1bEohfMiqbW/qrOPWffEv1uaWJdcX0zn\nRLq9ORbtOieiTBhtGQ7kdOBRd1/j7nvc/UGgBzCwjfuJq1wdi1IQ+li4e527/8Dd+7r7kQRD5C9o\n635iLFfHopS09lf1EuDYtPIQYK27b6Q0zol0e3Ms2vmO0XbaPErQMVMJnAx8BgxsYrsbgVcI7mzo\nAJxPkCX3a8t+4rzk4lgQ/AUxBuhCcFniPIK/Mo4s9Ocp0LE4NLEYMILg8tyX27qfOC+5OBalcF4A\nFYn4bwIeImhdVzSx3RjgY4I/oHoAc4AbS+yc2Otj0d5zIuoP3oOg6bwVWAF8O1HfL/El2DdRrgR+\nA9QQXJNbAJzW2n6KacnFsQB6Am8QXKvcCMwDRkf92fJ4LE4BPgRqgaXAOWH2U0xLLo5FKZwXwGRg\nT9YyKfs4JLa9Mu334z6gU4mdE3t9LNp7TmhoEBERCSXqB/dERKRIKGGIiEgoShgiIhKKEoaIiISi\nhCEiIqEoYYiISChKGCIiEooShkg7mdnhZrbHzDokys+Z2flRxyWSL1GOVitSUty9GAfzEwlNLQyR\nhPQpTUWkMSUMKWuJeZ1/YmaLga1mNtHMlpnZZjNbYmZfS9u2g5n9IjHhzAfAuKx9zTGzixI/Tzaz\n6Wnrsi9ffdfMPki8z3IzO7cwn1ik/fQXlQh8GzgD+BT4KnCyu9eY2beAh83sCA/mlbiEIEkcC2wD\nniJziGlPKzc7SJuZVQH/CQxz9/fNrDdwYI4/k0jOqYUh5c6BX7p7tQdzSjzh7jUA7v4Y8D5wYmLb\nbwF3JLbdSDDUfHOT2LQ2uc0e4G/NrKu7r3X3d/b+o4jklxKGCKxO/mBmF5jZ22a20cw2AoMJhoIG\nOCR9W4I5J9rM3WuBs4EJwBozm2VmR7cvdJHCUcIQSVw+MrP+wL3A94ED3L0HwWxsydbCxwRzDiSl\n/5xtK8HcJUkHZ7yh+x/c/bRE/V+BX+/NBxApBCUMkZQqguSxHuhgZv9I0MJIegz4oZn1MbMewNUt\n7OvPwKlmdpiZdQeuSa4ws15mNj7Rl7GLYNKj+hx/FpGcU8IQSUj0I/wH8BrBLGWDgblpm/wa+D2w\niGCmwydppnPb3f8IzAAWA28Cv0vbtgPBTGjVBB3tpwCX5vbTiOSeZtwTEZFQ1MIQEZFQlDBERCQU\nJQwREQlFCUNEREJRwhARkVCUMEREJBQlDBERCUUJQ0REQlHCEBGRUP4XTctvVZTjLtkAAAAASUVO\nRK5CYII=\n", "text": [ "" ] } ], "prompt_number": 8 }, { "cell_type": "code", "collapsed": false, "input": [], "language": "python", "metadata": {}, "outputs": [] } ], "metadata": {} } ] }