{ "metadata": { "name": "", "signature": "sha256:825f01b55aa6ccaeb57d4bd83256817c03351e3bf9c414914854dcc5ffaadf7e" }, "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": 1 }, { "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.028), feh=(-0.27,0.08))" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 2 }, { "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": 3 }, { "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": "iVBORw0KGgoAAAANSUhEUgAAAYwAAAEXCAYAAAC+mHPKAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XucFNWZ//HPw4DADIqjCF4QXDQqgQ2iKEjUzUIiCnHJ\nLzExatRsjC4aE3/GTaLRZTFuxEuibMSgJmqU+POH10gQzUUhiiiKGMgiURG5ODgIglwGhuuzf1T3\n9GVuNUN3V3X39/161YupU9XVT5dlP33qnDrH3B0REZHWdIg6ABERKQ5KGCIiEooShoiIhKKEISIi\noShhiIhIKEoYIiISihKGiIiEUpCEYWb7mNl9ZrbczDaZ2Ztmdkba9pFm9nczqzOzF8ysTyHiEhGR\n8ApVw+gIrAROc/f9gOuBR82sj5n1AJ4ErgOqgfnAtALFJSIiIVlUT3qb2ULgBqAHcKG7n5IorwTW\nAce5+zuRBCciIo1E0oZhZr2Ao4H/AQYAC5Pb3H0rsBQYGEVsIiLStIInDDPrBDwM/CZRg6gCNmXt\ntgnoVujYRESkeR0L+WZm1gGYCtQDVySKtwD7Ze3aHdic9VqNkigi0k7ubnt7jILVMMzMgPuAg4Cv\nuPvuxKbFwKC0/aqAIxPlIiISEwVr9DazuwkSw+fdvS6tvAdBm8W3gJnAT4BT3H141usdQMOxQ5B7\ndS5A5yJdbM9FIi4KGFdsz0UE0s7FXtcwCpIwzKwv8D7BrajdaZsudfdHzGwkMBnoC7wKfNPdV2Yd\nQwkjQf8zpOhcpMT2XChhRKroEkYuKGGk6H+GFJ2LlNieCyWMSOUyYRS00VtyQ/8TpOhcpOhcpOhc\n5IfGkhIRkVCUMEREJBTdkhKR/NLtoZKhGoaIiISihCEiIqEoYYiISChKGCIiEooShoiIhKJeUlL0\n5s6Fr3yl+e1/+hMM1Owq0YngSW/JDyUMKXo7dkBtbfPbd+0qXCwipUy3pEREJBTVMKSkDBsGTz4J\nn/88vPVWqqyiIrXPN74B99wDI0ZAXV3Tx7nkEvj2t/Mfr0gxUcKQktK5MxxyCHTqlCrbvj1znx07\ngn/feAM2ZU8OnDB6dH7iEylmuiUlZWf9evjb32D37tb3FZEUJQwpSa+8Alu2pJYpU1Lbpk+Hz3wm\n83bUn/8M3/pW4eMsC+7qIVUilDCkJHXtClVVqaVz55b3HzIEevcuTGwixUptGFIWqqtbfhajg346\nibRKCUPKwpe+FCwi0n76XSUiIqGohiFFZ/ZsWLUqtb5kSWShiJQVJQwpOpMmwdNPRx2FhFbkY0nN\nnz+fqVOnMmTIEF5++WV+8IMfcOSRR7Z533nz5vHyyy+zadMm5s6dy/XXX89pp51WyI+y15QwRKSs\nTZkyhcsuu6zJbdu3b+fss89m3rx59OrVi/79+3Puuefy2muvtWnfrVu38rvf/Y6JEycC8Pjjj3Pm\nmWfy7rvvcuihh+b18+WS2jCkqJ12WjDUR3L5wheijkiKzUcffdTsthdffJFu3brRq1cvAE444QSW\nLFnC8uXL27Tv0qVLueWWW1i2bBkAp59+Otu2bWPu3Lm5/0B5pBqGFLWrrlLvJ8mf5cuXc+CBBzas\nmxnV1dUsXryYI444IvS+Y8aMYe7cufTr1w+ADz74AIBPfepT+f8QOaSEISJ58/bbb/NjoC9gV1/N\nunXruO222+jZs2err122bBk33HADXbt2paKigrq6Om655ZaGX/DZtm7dyrBhw1i0aFHO4l+3bh2V\nlZUZZV26dGHz5s1t3nfYsGEN5RMnTuTqq69m0KBBOYu1EJQwRCQvNm7cyMiRI7kVOA/g5z9n4sSJ\njBw5kgULFtApfYTILO+//z4nnXQSv/71r/lSogr58MMPM2rUKN544w0q0ocfBl5//XUuv/xyFi9e\n3GJMH330EZMmTcLTGuDnzJlDfX19w3q3bt247rrrANh///0z9gXYsmULPXr0aHTssPved999HHbY\nYdx8880txhpL7l4UC+BBuFLuxo5NDk7k/tRTuTvu+PGp406YkLvjlqsf//jH3qtXL9+9e3dD2fr1\n671Tp04+ZcqUFl87duxYP/DAAzPKtm7d6l26dPHf/OY3DWVLlizxMWPG+De/+U0fNmyYd+jQoc1x\nTmjhP/YLL7zggwYNaljfuXOnd+7c2d9555127Ttjxgy///773d29vr7ely9f3uZ42yrtu3Ovv4fV\n6C0iefHYY48xdOhQOqSNu1JdXc2xxx7LY4891uzrduzYwTPPPMNRRx2VUd61a1f69u3L448/3lB2\n7LHHMmPGDB544AGOOeaYRr/w99app57K2rVrWZV48Ocvf/kLAwYMaGh7eP7551m4cGGoff/yl7+w\nZs0aRo8eTW1tLTNnzuTDDz/Mabz5pltSIk2orYXE9wAQTMCkecHD27x5M0uXLmXUqFGNth1yyCHM\nmzev2deuX7+e3bt306VLl0bb9ttvvya7tOZLx44dmTp1Kj/96U8ZPnw4s2bNYtq0aQ3b77rrLo4/\n/ngGDRrU4r7Lli3jrLPOYsuWLQ2vNTM2btxYsM+SC0oYIk24++5gSaquDubRkHBWrFgBBF/w2aqq\nqti0aRM7d+5ssh3joIMOorKyMqNdIWn16tWsW7eOPXv2ZNRc9kbXrl1b3D5ixAhGjBgBwIUXXpix\n7cknnwy1b79+/djU3GxdRUS3pEQk55Jfjvvss0+jbVVVVQB88sknTb62oqKCc845h/feey+jvLa2\nltWrVwNBLSRXfvjDH+bsWKVOCUMkoVevYGKl5PLpT0cdUfFK9mKy5LAgaXbu3AnA7hamPLzllls4\n4IADuPfeewHYtWsXd911F4MHD8bdG/WSksJQwhBJuPzyoN0iubz0UtQRFa+DDjootWKWGk8KqKur\nw8zYd999m319jx49mDdvHrW1tVx99dXccMMNXHLJJZgZXbt2pbq6Op/hSzPUhiEiOderVy/MjA0b\nNjTaVldXx/77799wa6o5+++/P+PHj88oW7NmDcOHD89prBKeahgiknNVVVUMHjyYlStXNtq2dOlS\njjvuuDYfc+3atdTU1PDVr341FyE2a/78+Vx55ZVMnTqVcePGNWpLCbvvvHnzuP3225kwYQKnn346\nL774Yl7jLohcPMxRiAU9uCcJ+XpwL9vHH6fep7o6f+9TqsaPH++HHXZY6iS6+9KlS93M/Je//GXG\nvu+8845v27atYf2OO+7wnj17ek1NTUPZbbfd5ocffrhv3769yfe76KKL2vXgXrr6+nrv27ev19bW\nurv766+/7ieeeGKb962rq/NrrrmmYd/HHnvMKysrMz5PoaAH90Qk7i677DLq6ur4bVrZnXfeyYAB\nA7jkkksayl588UWOOeaYhiFAILhtVVVV1dC4/eabb/Kzn/2MRx99tMmeVwD19fW4O9u2bWt3zBqd\ntmVqwxCRvDj44IOZPXs21x13HAuAzZdcwoYNG3juuefo2LFjxn49e/bMeLL76quvZuPGjVx77bVs\n3ryZjRs38vTTTzN06NCM91i7di0XXHABNTU1LF68GDOjb9++DBw4kG9/+9ucd955bYpZo9O2TAlD\nRPJm0KBBzGhluI6jjz6a2trajLIuXbpw6623tnr8gw46iOeee26vYkyn0WlbpoQhIiVv8uTJLTZe\nDxkyhPPPP5/u3bs3Go+qrEenzaKEISIl74orrgi1X//+/RseFoTggcH169fTt2/fRvsee+yx3HPP\nPS3u+8wzz9ChQwduvvlmtm/fTm1tbZPHKhZq9JZY27MH6uoyl127oo5KSpVGp22ZahgSa6tWQVZb\no0jeaHTalln2Pbi4MrPgYYwiiVdyY8WKlhPGU0/lb07v9esh2QmmY0c444zM7VddBYmBSUViKzme\nl7s3HtirjVTDkKKS1SmFQo1Bt2sXzJiRWZbnB45LR3IcKf3YK3oFa8MwsyvMbL6Z1ZvZA2nlR5jZ\nHjPbnLZcV6i4pHj06dO4PeOss6KOSqR8FLKGUQPcCIwCmpqxZD/X/SaJkW7dYPr0zLLbb4fZsyMJ\nRyRyBUsY7v4UgJkNAXo3sUsHoPkB8kUKbJ99GtdgWpiKWqTkRdGttrmGlxVmtsrM7jezA5vZR0RE\nIhJFwsi+7bQWGAL0AU4A9gUebu7FZtawiIhIY/n6noyil1TGJ3D3OmBBYvUjM7sC+NDMqhLbRKSY\nqWmyZESRMMJePU3WftQuLiLSsvTvyVzWMgqWMMysAuiUeM8KM+tM0Mh9PLAReBeoBn4BzHL3xsND\niohIZApZw/gPIH2C3m8AE4B3gJuAnsAm4I/AuQWMS6Rd3noLZs1KrXfrBieeGF08IvkWamgQM/t3\nd/9ZE+Xfd/fb8xJZ4/fS0CBl4L33oKYmtV5bC+ecE/zdp08wVEiULrwQpk5tetvgwbBgQdPbRKIS\nxdAg/wk0ShgEtYaCJAwpD5Mnw6RJUUchIk1pMWGY2QiCXk0Vib/THUlwC0mkbPTvD5/7XGp982Z4\n443IwikOGkuqZLRWw7ifoFdTZ+C+tHIH1gDfzVNcIvzDP0DvtDEBevWKLpaka68NlqQFC+CEE6KL\nR6SQWkwY7n4EgJlNdfcLChKRSMJ3vxsMIS4i8RCqDSM9WZhZh6xte3IdlIiIxE+ooUHM7AQze8XM\ntgK70pad+QxORETiI2wvqQeB6cDFwNb8hSMiInEVNmH0Aa7TfBUi0mb62igZYUerfYpg4iMRESlT\nYWsYXYGnzOwlgu60Se7uF+Y+LBERiZuwCeOtxJJNdU0RkTIRtlvthDzHISIiMRcqYTQxLEgDd38h\nd+GIiEhchb0llRwiJOkgguFCVgH9ch2UiJQQjSVVMsLekjoifT0xGdL1wJY8xCRl5PLLM9fnzIkm\nDhFpXbsmUHL33WZ2E/AB8PPchiTlZMqUqCMQkbDCPofRlC8QTLEqIiJlIGyj96qsokqgC3B5E7uL\ntMvkyanb3QCf/Wx0sYhIY2FvSWUPbV4HvOPuG3Mcj5Sxyy/PTBgiEi9hG71nQ8PQ5r2ANRrWXERC\nUe+okhF2ePP9zOwhoB6oAerN7CEz657X6EREJDbCNnrfCVQBAwnaL5L/3pmnuEREJGbCtmGcAfRz\n97rE+jtm9k1gWV6iEhGR2Albw9hG8HR3uh4Et6hERKQMhK1h/Br4k5n9HFgBHAFcBfwqT3GJiEjM\nhE0YNwGrgfOBQxJ/30IwxpSIAH/7G/TunVn2xBMwdGg08cSGxpIqGWG71e4hSA5KECLN2LULamoy\ny3bsiCYWkXwI2632TjMbnlU23Mwm5ScsERGJm7CN3ucCb2SVLSC4RSVStgYOhFWrMpcTT4w6KpH8\nCNuGsYfGyaUDoIEcpKzts0/jdovOnaOJRSTfwtYw5gD/lRgaJDkfxg3AS/kKTERE4iVsDeNKYAZQ\na2YrgD7Ah8BZ+QpMREqEekeVjLC9pFaZ2fHAScDhBFOzztMAhNIW27dDvR71FClaoWfcc/fdwCuJ\nRaTN7rkHrrwy6ihEpL32ZsY9EREpI+2a01tkb3XqBF27Rh2FiLSFahgSiXHjYOPGzEWz7YnEW7M1\njGQX2tao4VtEWqSxpEpGS7ekdmWtO40f1HOgIqcRiYhILLWUMPql/T0GOJtg1NqVBM9hXAM8kb/Q\nREQkTppNGO6+PPm3mX0fGOLuGxJFb5vZfGA+8Mu8RihSxF56CdavT60fcgicdFJ08YjsjbC9pPYj\nmMN7Q1pZJdA95xGJlJDrrstc//KXgzkyRIpR2ITxIPBnM7uD4CnvPsD3gIfyFZiIiMRL2ITxQ2Ap\n8HWCGfc+BO5EU7SKNHLKKXDAAan11ath/vzo4omcekeVjLbMuHd3YhGRFkycmLn+xBNw9tnRxCKS\nS2Fn3OtgZpea2Qtm9rdE2Wlm9rX8hiciInER9knvG4CLCW5B9UmU1RB0rQ3FzK4ws/lmVm9mD2Rt\nG2lmfzezukRS6tPccUREJBphE8a/Al9090cIZt8DeJ/MZzVaUwPcCNyfXmhmPQie57gOqCboqjut\nDccVEZECCNvo3QHYklVWBWwO+0bu/hSAmQ0B0ie1/DLwP+7+RGL7BGCdmR3t7u+EPb6IiORX2BrG\ns8DtZtYFGsaZuhH4fTveM3t4kQHAwuSKu28l6JE1sB3HFpG4MdPIkiUibML4PnAw8AnBQ3xbgCNo\nQxtGmuw+dlXApqyyTUC3pl5sZg2LxFt9PWzbllp27Ig6IpHykK/vybDdajcC/8fMegF9gVXu/mE7\n3zP7E2whSELputOG210ST5/+NLz/ftRRiEiuhO1W+4CZjXL3Ne7+WjJZmFl7xpHKrmEsBgalvVcV\ncGSivPGL3RsWERFpLF/fk2FvSX0DuN/MfpBVfkHYNzKzikQbSEegwsw6m1kF8BQw0My+nNj+n8Bf\n1eBdWvbZB7p0SS2dOkUdkYi0VdheUvXAUOBpMxsEXOzu29v4Xv8BjE9b/wYwwd1/YmZfASYDvwVe\nJRiCRErIkiXQry2dsEUkdkJP0eruHwCnJF7zkpn1buUl2a+f4O4dspafJLY97+793b3S3Ue4+8o2\nfQoRiS93jSdVIto0p7e7b3P384AngXnAPnmJSkREYifsLakb01fc/WYzWwhoLCkRkTIRtlvtrU2U\nPUvwQJ+IiJSBZhOGmf3B3Ucl/n6pmd3c3U/LS2QiIhIrLdUw0mfTu6+ZfdSSJSJSJppNGO7+cNrf\nvylINCJSepLDU6inVNFr6ZbUxYSoQbj7/a3tIyIixa+lW1IXEO6WkxKGiEgZaOmW1OcKGIeIiMRc\n2OcwGlgwXm7DiLPuvqeF3UUkzcyZ0CdrAuJXX4VDD40mHpG2CJUwzOwwgrGe/olg6PFkwnCgIj+h\niZSe+npYtSqzbPfuaGIRaauwQ4PcDewERhDMX3E88DRwWZ7iEpFSobGkSkbYW1KfBfq4+xYzw93/\nmuhFNRe4N3/hSbFYuxYWLcos27Ytmlji5swzYcWKzLKhQ6G2Npp4RNorbMLYlVgANphZT2AjcFhe\nopKi88orMHZs1FHEU2Vl43aLjm1uPRSJXthbUq8BZyb+/gMwjWDio/n5CEpEROIn7O+cb5BKLlcB\nVwPdgEn5CEqK24EHwqBBmWVdukQTi4jkTtjRaj9J+3srWcOdi6QbPhymT486ChHJtbDdajsB5wKD\nCWoWSe7ul+YjMBEpERpLqmSEvSU1FfhHgvkvkn07DI1WKyJSNsImjDMIutVuymcwIiISX2F7SS0B\nDshnICIiEm9t6SV1v5k9C6xJlBlBG8ZDzb9MRERKRdiEcRHB0977AdnP7yphiIiUgbAJ4/8Cg939\nrXwGIyIlSL2jSkbYNow1wMp8BiIiIvEWtoZxOzDVzG4l1YYBgLsvy3lUIiISO2ETxl2Jf7OHl9N8\nGCJ76ayzoHPn1Pr558P3vhddPCLNaTVhJGbY+xSw0t135j8kkfKycGHm+qmnRhOHSGvCtmEsAjQv\nmIhIGWu1huHubmZ/BY4heIBPRPbS9OmwY0dq/be/hcmTo4snrzSWVMkI24YxC3jWzH4DrCJou0g+\nuHd/nmITKVmDB2euv/RSNHGItEXYhHEKsBz4pya2KWGUoRkzYPXq1Hr2fXgRKT1h58P4XJ7jkCJz\nxx3wwgtRRyEihRR6ZmEzqwb+BTgUqAFmuPv6fAUmIiLxEnYCpZOBZ4C/AyuAs4BJZvZFd5+bx/ik\nCIweDYcdllr/zGeii0VE8idsDeO/gcvd/f8nC8zsnET5ifkITIrH978PI0dGHYXElnpHlYywz2Ec\nDTyaVfYEwQN9IiJSBsImjHcJ5vRO91VgaW7DERGRuAp7S+pK4Bkz+y7BqLV9CWodX8xXYCIiEi9h\nu9XONbMjgTEEvaSmA8+6+8f5DE5EROIjdLfaRBfaqXmMRUSAefPg5ptT62bwox9FF49IUosJw8xm\ntfJ6d/cROYxHpOzNmRMsSRUVRZ4wNJZUyWithvFwE2UOHEbQrlGZ84gkdmprYdKkzLKl6u4gUnZa\nTBju/uv0dTPrAVwDXAJMA36Sv9AkLtauhVtuiTqK0nbSSZm1iD174LbbootHpClhn/TuDvw78F1g\nBnC8u7+Xz8BEyslppwVL0q5dShgSP621YVQS3Hr6d2A28Fl3X1yAuCSmevaEq67KLDvyyGhiEZHC\naq2G8T7Bw323AvOBXmbWK30Hd8/JmKVmNhsYCuxKFH3g7v1zcWzJnZ494Zproo5CRKLQWsLYlvh3\nXAv7/EOOYnHgO5qQSaTEqHdUyWit0fuIAsWRZAV+PxERCSnsWFKFMtHM1prZHDNranY/kbKzezdU\nVmYud98ddVRSjkI/6V0APwIWAzsIBjr8vZkd5+7L0ncyS1VCXFVdKRPbtmWu79rV9H4ikPk9mUux\nqWG4+2vuXufuO939IeBlYHTUcZWj3buDL6T0RUQkNgkjLHdvWCQ/Lr0UOnVKLccfH3VE5aeiAurq\nMpdLL406KikW+fqejEXCMLPuZjbKzLqYWUczOx84FXgu6thEomDWuN2iU6eoo2ons9R4UlLU4tKG\n0Qm4ETgW2A0sAca6u0YsipgZdEj7WVFREV0sIhKtWCQMd18HnBR1HNLYr34FF18cdRQiEgexuCUl\nIiLxp4QhIiKhKGGIiEgosWjDEJESpi7wJUM1DBERCUUJQ0REQlHCEBGRUJQwREQkFCUMEREJRb2k\nRIrQlCnw7LOp9UMPDZ7Kj6XkOFLqLVX0lDDK3O7dsGhRZtnHH0cTi4T31lvBknTUUdHFIuVDCaPM\nbdum4ctFJBwlDJEiMW4cnHlmar2mBv7t36KLR8qPEoY0MINBgzLLDjggmliksYEDgyXp3Xeji0XK\nkxKGNKishDffjDoKEYkrJQyRErBlC0yfnlk2dCj06hVNPBnUO6pkKGGIlIDaWhg7NrNsxgwYMyaa\neKQ06cE9EREJRTUMkSJVVQVnnZVZ9uqrsHZtNPFI6VPCEClShx7auN1izBiYOTOaeKT06ZaUiIiE\nooQhIvlllhpPSoqaEoaIiISihCEiIqGo0VukRI0fD5Mnp9YHD4abboouHil+ShgiJWrBgsz1Xbui\niUNKhxJGmXngAbj33tT67t3RxSIixUUJo8zU1AQPd0lpuuEGuOKK1PqCBXD99dHFA2gsqRKihCFS\nQoYMyVzvqP/DJYd0OZWxiy7KnICng/rMiUgLlDDKWO/ecPLJUUchIsVCvylFRCQU1TBK3He+E0yu\nk7RwYXSxSLSefx66ds0sW7AA+vePJh4pPkoYJW7aNPj446ijkDhwh/r6xmV5lxxHSr2lip5uSYmI\nSCiqYZSRO++EffdNrX/mM9HFIoUxYgRs3ZpZNngwvP12NPFIcVPCKCNf/zr06BF1FFJIFRWN2y3U\nfVraS5eOSBkbMCA1XYUZXHpp1BFJnKmGUULmzIFf/jKzbPPmaGIRkdKjhFFCVqyARx6JOgqRLOod\nVTKUMETKzOLFmev33gvjxgV/P/ssjB6duf3JJ6FLl8LEJvGmhFHExo6FRYtS6+kP6A0fHjy0l65b\nt8LEJfGWPb12+voHHwRLutGjobIytX7hhfC1r+UvPokvJYwisXo1PPpoZtn06c3vf8QRcN55eQ1J\nysSsWZnrn/1sNHFI9JQwisT778NVV0UdhZSiM86AZ57JLBszpvn9//jHzGc7qqrgmmvyE5vEi3mR\nNEiZmQMUS7y59vLLcMopzW+fNg1OPDG1XlUFPXvmPy4pTTNnwp49qfUHH4THH2963549Yc2awsQl\nbWeJe47ubq3s2qrY1DDM7ADgPuALwDrgWndXn58mHHwwnHNOZtmwYdCnTzTxSOnJbvh++eXm9/3o\no+CJ8nTXXQcjRyZWNJZUyYhNwgDuAuqBnsBg4BkzW+jub0UbVvz06weTJkUdhZSTL3whs+F7zRq4\n667UenY7x4ABqSfMh+c/PCmQWNySMrMqYD0wwN2XJsoeBFa7+7WJ9ZK9JfXxx003UO/aBStXwo9/\nHNxiStUqGqqYBYsxrtKq2xFHEr1Cnos1a4KabhhO8zWMLVuCGkpz+vRp3zSzui5SSvGW1NHArmSy\nSFgIfC6acHKre3fYsaPpbT/9aTDG0x//2Pzrv/Wt/MQl0l7V1cH8GunOPbflL/+ZM4N/t24NngU5\n+WT4wx/g9tubf01tbfAMyCefNL3dLEgqv/gFPPdc4+2jR8Phh8M997T8eSScuNQwTgUedfdD0sou\nAc5z939OrBdtDaNz5+YTRnL79u1tOaJ+PSXpl2RK1Odi/Hj4859T66++GlQqkjUMIz9xVVVBXV12\nafLHtNO/P7xVxje2S7GGsQXYL6usO9BoJCTLfuqoBLQtWaSU4rloL52LlLidC2vir1xqnCwy333J\nksYPK0r7xGW02neAjmZ2VFrZIOB/IopHRESyxOKWFICZPQI48G3geGAGcLK7L4k0MBERAeJTwwC4\nHOgKfAT8FhinZCEiEh+xqWGIiEi8xamGISIiMRZpwjCzA8zsKTPbYmbLzezcFvb9DzNbZWafmNks\nM/t0e44TVzk8F7PNbJuZbU4sRXdbL+y5MLPOZnaHmdWY2Xozu8vMOrb1OHGWw3NR1NeFmV1hZvPN\nrN7MHmhl36vM7EMz22hm95nZPmnbSuGayNW5aPM1EXUNI304kPOBKelffklm9i/AOOBU4ADgFWBq\nW48Tc7k6Fw58x933TSz98x557oX973kNQQeJAQQPfx4PXN+O48RZrs5FsV8XNcCNwP0t7WRmo4Af\nASOAvkA/4Ia0XUrhmsjVuWj7NeHukSxAFbAdOCqt7EFgYhP7XgtMS1sfAGxr63HiuuTqXCTWZwEX\nR/2ZCnQuXgfOTls/F1hZptdFs+eiFK6LtM9xI/BAC9v/H/Bfaev/DHxYKtdErs5Fe6+JKGsYzQ0H\nMqCJfZ8HTjazT5lZJ+Ai4Nl2HCeucnUukiaa2Vozm2Nm/5SfkPOmrf890x/J6gD0NrN923GcOMrV\nuUgq5usiqbVH8D5NcI6SFgG9zKya0rgm0u3NuUhq0zUR5ZPe3YBNWWWbgX2zd3T31ywYjPBtYDew\nEkgOnhz6UeDHAAAEMklEQVT6ODGWq3MBQRV0MbCD4Ffm783sOHdflo/A86At/z2fA640s1kE1/L3\nCKrZlW08Tlzl6lxspvivi6TWunV2AzamrSfP376UxjWRbm/OxQbacU1EWcNoy3AgVxB8KfYGOgM/\nAV4ws65tOU6M5epc4O6vuXudu+9094eAl4HR2ceJsbb89/wp8CbwV2AO8BTBL8g1bTxOXOXqXJTC\ndZHU2q/q7HPWPfHv5ia2JbcX0zWRbm/ORbuuiSgTRluGAzkDeMTdV7v7Hnd/EKgG+rfxOHGVq3NR\nCkKfC3evd/fvuntvdz+KYIj8+W09Tozl6lyUktZ+VS8GjktbHwSscfcNlMY1kW5vzkU73zHaRptH\nCBpmKoFTgE+A/k3sdxPwEkHPhg7ABQRZcr+2HCfOSy7OBcEviFFAF4LbEucT/Mo4qtCfp0Dn4tDE\nYsAwgttzn2/rceK85OJclMJ1AVQk4p8IPERQu65oYr9RwIcEP6CqgdnATSV2Tez1uWjvNRH1B68m\nqDpvAZYDX0+U90l8CfZOrFcCvwZqCe7JzQdOb+04xbTk4lwAPYDXCO5VbgDmAiOj/mx5PBenAu8D\ndcAS4NwwxymmJRfnohSuC2ACsCdrGZ99HhL7XpX2/8d9QKcSuyb2+ly095rQ0CAiIhJK1A/uiYhI\nkVDCEBGRUJQwREQkFCUMEREJRQlDRERCUcIQEZFQlDBERCQUJQyRdjKzI8xsj5l1SKzPNLMLoo5L\nJF+iHK1WpKS4ezEO5icSmmoYIgnpU5qKSGNKGFLWEvM6/9DMFgFbzOw6M1tqZpvMbLGZfSlt3w5m\n9rPEhDPvAWOyjjXbzC5O/D3BzKambcu+ffVNM3sv8T7LzOy8wnxikfbTLyoR+DpwJvAx8EXgFHev\nNbOvAb81syM9mFfiUoIkcRywFXiSzCGmPW292UHazKwK+G9giLu/a2a9gANz/JlEck41DCl3DvzC\n3Ws8mFPicXevBXD3R4F3gZMS+34NuCOx7waCoeabm8Smtclt9gD/aGZd3X2Nu7+19x9FJL+UMERg\nVfIPM7vQzN40sw1mtgEYSDAUNMAh6fsSzDnRZu5eB5wDjANWm9kMMzumfaGLFI4Shkji9pGZ9QXu\nBb4DHODu1QSzsSVrCx8SzDmQlP53ti0Ec5ckHZzxhu5/dPfTE+V/B361Nx9ApBCUMERSqgiSxzqg\ng5n9K0ENI+lR4HtmdpiZVQPXtHCsvwKnmdnhZtYduDa5wcx6mtnYRFvGToJJj3bn+LOI5JwShkhC\noh3h58ArBLOUDQTmpO3yK+APwEKCmQ6foJnGbXf/MzANWAS8Dvw+bd8OBDOh1RA0tJ8KXJbbTyOS\ne5pxT0REQlENQ0REQlHCEBGRUJQwREQkFCUMEREJRQlDRERCUcIQEZFQlDBERCQUJQwREQlFCUNE\nREL5X6Osmd8YaAQJAAAAAElFTkSuQmCC\n", "text": [ "" ] }, { "metadata": {}, "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAYwAAAEXCAYAAAC+mHPKAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XmcFOW97/HPj1XZMbIpEURcEBNESdyCeMAjuOKWHBc0\ner1G8WISEw2IGwaDaI5RD6LGuBH0GsWAURQ9HlkU8IKEiMYFZFcCKLLMsAzr7/5RPTPdPVvNTHdX\nzfT3/XrVa7qeqq76dVH0r5+qp57H3B0REZGqNIg6ABERqRuUMEREJBQlDBERCUUJQ0REQlHCEBGR\nUJQwREQkFCUMEREJJScJw8yGmdkCMysys2fSlp1vZp+YWUHi7+BcxCQiItVjuXhwz8wuAPYBA4H9\n3f3qRHl7YAVwobu/ZWZnAZOALu6+IeuBiYhIaDmpYbj7FHf/G/Bt2qLuwFZ3fyux3hvANuCwXMQl\nIiLh5foehqXNLwL2mNk5ZtbQzM4HioCPchyXiIhUoVGO95dy/cvdt5nZdcCLQBNgF3Cxu+/IcVwi\nIlKFXCeMlBqGmR0HPAH0dfeFZtYHeNXMznT3RUnrqYdEEZEacvf0qzs1kutLUulf/AOA/+fuCwHc\nfQEwDzg9x3GJiEgVclLDMLOGQOPE/hqaWVNgL8E9jN+YWS93X2RmvYG+wPjytqOu2MEs+KGgY6Fj\nkUzHolROj0ViX8T0uBcfi0zJVQ3jDmA7MBwYAuwARrr7fwP3A5PNrBB4Gfidu/9PjuISEZGQcvIc\nRm0V38OoC7Fmm35JltKxKKVjUUo1jFJJxyIjVY1c3/SWWtIXQikdi1I6FqV0LLJHfUmJiEgoShgi\nIhKKLklJpEaNgmnTyl/WvDlMn57TcESqJ88ufylhSKSWLYP588tf1qpVbmMRkcopYUjO7NgBq1en\nlm3ZEk0sIlJ9ShiSMx9/DCecUPHyO++EU0+F0/Wcv0gs6aa3xEb37tCnT9RRiEhFVMOQSDRtCl26\npJbpnoVIvClhSCS+//3yb3brnobUKTF/0jvTdElKRERCUcIQEZFQlDBERCQU3cOQGpk5E/buLX9Z\nx47Qs2dOwxGRHFDCkBoZPBgKCspfNmQITJyY23hEJPtycknKzIaZ2QIzKzKzZ9KWNTOzR83sGzPb\nbGazchGTiEituedNCynIXQ1jDTAaGAjsn7bsCYLEdRSwETg2RzFJhpx6KnzzDXz2WdSRiEg25SRh\nuPsUADPrA3QuLjezo4BzgYPdfWui+B+5iEky59VXg+nKK6OORESyKdetpNKHCfwhsAr4beKS1Edm\ndmGOYxIRkRByfdM7/WJfZ+AY4GWgE3Ay8LqZferun6e/uXh8WtAwjLlUWAjbtqWW7dsXTSwiUrnk\n78lMy3XCSP8kO4DdwD3uvg9418xmAGcAZRKGRGPsWBgzpvrvu+wy+Oc/S+e3b89cTCKSe1HXMD5K\n/E1PJOVWH1SrqFuWLg26NBepyxYsWMDEiRPp06cPc+bM4ZZbbuGwww4LFqb1JVXZuvPmzWPOnDkU\nFBQwd+5cbr/9dk499dSMx5v8PZnp2kZOEoaZNQQaJ/bX0MyaAnuAWcBq4FYzGwucAJwG3JyLuKT6\nmjeHFi1Sy9LPyWXLYMKEoOWUSNw99thjDB06tNxlO3fu5OKLL2bevHl06NCBHj16cOmllzK/nJ4z\nK1t3+/btvPLKK9x7770AvPzyy5x55pl88cUXHHTQQVn9fJmUq5vedwDbgeHAEIJLUbe5+x5gMHAW\nsBn4I3CFuy/JUVxSTbfeCuvWpU7p3ZK//z5cdRWsXFla9txzsGhR6fT887mMWqRiX3/9dYXL3n33\nXVq0aEGHDh0AOP744/nss89YmXxyh1h36dKl3HfffSxfvhyAM844gx07djB37tzMf6AsylWz2lHA\nqAqWfUpws1vqscMPD7o0F6lLVq5cyXe+852SeTOjbdu2fPLJJ3Tt2jX0umeffTZz586lW7duAHz1\n1VcAHH744dn/EBmkrkEkI7p1q/w5jAMPzF0sEh+LFy9m5MiRdOnSBTNjw4YN/P73v6d9+/ZVvver\nr75i7NixmBk7d+5k+/btjBgxgmOOOabS9/3nf/4ne/bsYcSIEbWOf8OGDTRr1iylbL/99qOwsLDa\n65544okl5ffeey+//vWv6dWrV61jzCUlDMmIU04JJpFiW7ZsYcCAAdx///1cdtllQPBFOWDAABYu\nXEjjxo0rfO+3337LjTfeyOOPP15yiWflypUMGjSIqVOn0r1793Lft2rVKu666y6GDx9e7vKvv/6a\nhx56KOXG8OzZsykqKiqZb9GiBbfddhsAbdq0KdPYZuvWrRxYzi+gsOs+9dRTHHzwwYwdO7aijx9f\n7h77iaDVlEs0Ro4s7jDH/Z57sruvzZtL99WqVXb3Jdk1cuRI79Chg+/du7ekbOPGjd64cWN/7LHH\nKn3v/fff7+PGjStTPmbMGB8+fHiF77v22mvdzPzuu+8OHeeoUaMqXDZ9+nTv1atXyfzu3bu9adOm\nvmTJkhqtO3XqVH/66afd3b2oqMhXrlwZOs6aSPruzMh3scbDEJGsmDRpEieccAINGpR+zbRt25aj\njjqKSZMmVfre5cuXM3369DLlTZs2Zc+ePeW+Z/LkyQwYMKB2Qafp27cv33zzDV9++SUAs2bNomfP\nniX3Ht555x0WLVoUat1Zs2axfv16zjrrLNatW8cbb7zB2rVrMxpvtumSlIhkXGFhIUuXLmXgwIFl\nlnXq1Il58+ZV+v7evXtz/fXXc/nllzN+/HjatGnDzp07ef7555kwYUKZ9bdu3cq0adP405/+xKWX\nXpqxz9GoUSMmTpzI7373O04++WRmzJjBiy++WLJ8/PjxHHfccfTq1avSdZcvX865557L1q1bS95r\nZmypY4PYK2GISMatWrUKgFbpba6B5s2bU1BQwO7duyu8j3HVVVfx5JNP8sILLzBjxgzuv/9+3nrr\nLcaMGVPuTe+xY8cycuTIGsW6//7pHWin6t+/P/379wfgyrSWHZMnTw61brdu3SioaACZOkSXpEQk\n44q/HJs0aVJmWfPmzQHYvHlzhe9v0qQJ06dPZ+DAgaxbt44rr7ySzZs307t37zLrfvjhh7Ro0YJD\nDz20RrH+5je/qdH78pEShohkXMOGDYHyu6bYvXs3AHsrGuM34ZFHHuGQQw5h6tSpHHroobz++usc\ne+yx/DOpg7J9+/bx0EMPcfPN6hwiF3RJSmJvzRpIPOdUrt69oZwfshKhdu3aVbhs27ZtmBktW7as\ncJ0HH3yQOXPm8NprrwFw2mmncfvtt/Pwww8zZMgQPvzwQwCeeOIJrr76aho1iuirLK0vqfpOCUNi\n79ln4fbbK16+di107JizcCSEDh06YGZs2rSpzLJt27bRpk2bkktT6dyd0aNHM3v27JKyZs2a8Yc/\n/IEjjzySoUOHsnjxYlq3bs2nn37K9ddfX+42JPOUMCS2CgvhiCPgiy+ijkSqq3nz5vTu3ZvVq1eX\nWbZ06VKOPbbikZi/+eYbNm/eXG63Gddddx3Dhw+nqKiIBQsW8Pnnn3PBBReULC++3PWXv/yFDz/8\nkCuvvDJleRiV9k6bvi4wEegzcWJkvdPmVKYe6MjmhB7ci1RUD+5VNP3wh+6NGpXOr12b3ZikZu68\n804/+OCDU8qWLl3qZuaPPvpoSvmSJUt8x44d7u6+b98+79ixo8+ePbvMNgsKCrxz586+a9eucve5\ncuXKaj+4l6yoqMi7dOni69atc3f3Dz74wH/wgx9UvC74usR3U/K627Zt8xEjRpSsO2nSJG/WrJmv\nWbOmRnHVFHpwT/LZyJEwbx4k9fEmMTV06FC2bdvGc889V1I2btw4evbsybXXXltS9u6773LkkUdy\n/vnnA8GN8ocffpjrrruupHkuBC2vrrnmGh5++OEKm+Pu2rUr5W91Vbt3WqBDYr4+9k6bTpekJFZa\ntIDFiytefsABuYtFaqdjx47MnDmT2267jYULF1JYWMimTZt48803U25Sd+zYkfbt26f0D/WTn/yE\nTp068ctf/pKmTZtiZuzdu5dhw4aVe1mnsLCQwYMHs3jxYsyMBx54gOnTpzN8+HAGDx4cOuZq906b\nNF8fe6dNp4QhsdKwYXDfQuqHXr16MXXq1ErXOeKII1i3bl2Z8r59+9K3b99Q+2nZsmW5XYlUV7V7\npx00CKZNK3fd+tA7bbpcjbg3DLgKOAZ4wd2vLmedOwnGzDjd3Wv/Ly8ikvDII4+wbNmyCpf36dOH\nyy+/nNatW5dpYZXXvdOmyVUNYw0wGhgIlHkO38wOAy4G/pWjeKQC69cHY1skq+HlYJHYGDZsWKj1\nevTowRNPPFEyv2fPHjZu3EiXLl3KrHvUUUfxxz/+sdJ1X3/9dRo0aMDYsWPZuXMn69atK3dbdUVO\nbnq7+xR3/xvwbQWrPEIwfOvuXMQjFXOH7dtTpwo6BxWpd9Q7beVyfQ+jTD8BZvZjoMjdp5XXjYCI\nSK6od9rKWfo1uKzuzGw00Ln4HoaZtQT+TnDfYrWZrQCuSb+HYWZlgsxl3Plk3Tro1Cl43a4drFiR\nurxx43h0w9GxY3D5DPSkt0iy8n54u3tGfo1HXcMYBUx099WVrCMRadAAKui9QUQg7/qSyvWDe+lH\ntT/wczNba2Zrge8CL5nZLeW+OfXpbxERSZPN78lcNattCDRO7K+hmTUF9gIDkmIw4APgJuDNXMQl\nIiLh5aqGcQewnaAl1BBgBzDS3Te6+9eJaT1BEtnk7ttyFJeIiISUkxqGu48iuF9R1Xo1GzJLRESy\nTl2D5LnhwyG5R4XddfBJmCVLIHnYhZYtoXPn6OIRqa+UMPLcsmWwYEHUUdROv36p8xddBC+/HE0s\nkmfyrAGOujcXEZFQVMOQEvfeCwMGlM5HNUxyGIcfDm3bls4XFMC/1BOZSFbF+CtBcq17d/jBD6KO\nIpz33kudf/ll+PGPo4lFJF+EuiRlZjdXUP6rzIYjIiJxFfYexl0VlN+RqUBERCTeKr0kZWb9CZ7A\nbph4newwoCBbgYmIxF6e9SVV1T2Mpwn6f2oKPJVU7sB64MYsxSUiIjFTacJw964AZjbR3a/ISUQi\ntbR+Pbz9dmpZnz6prapEpPpCtZJKThZm1iBt2b5MByVSG7NnwxlnpJbNnFn2AT8RqZ6wraSON7P3\nzWw7sCdpqoMdSYiISE2EfQ5jAvAqcA1Br7MisdK+PZx+emrZBx9AHR8RUyRWwiaMQ4DbXCMXSUyd\nemrZ+xb9+sG770YTj+SJPPtKDPscxhRgYDYDERGReAubMPYHppjZf5vZxKTpz2F3ZGbDzGyBmRWZ\n2TNJ5Sea2dtm9q2ZfW1mL5lZx+p+EBERya6wl6Q+TUzpqlMfWwOMJqip7J9U3gZ4HHiLYMS9R4Bn\ngDOrsW0REcmysM1qR9V2R+4+BcDM+gCdk8pTxu82s/HAzNruT0REMitUwiinW5AS7j69omUVba6K\n5acC/6zmNkVEJMvCXpIq7iKkWDuC7kK+BLpVc58VXsYys+8TdGh4XgXLSzeSZ60TRCSGYtiXVPL3\nZKaFvSTVNXnezBoCtwNba7DPcj+NmXUH3gB+7u5zarBdERHJohoN0erue4ExwG9q8vb0AjPrArwN\n/Nbdn69kvyWTiIiUlc3vydqMuPfvBK2aQknUShon9tnQzJoSdC/SEZgOPOLuT9QiHhERyaKwN72/\nTCtqBuwH3FCNfd0B3Jk0PwS4m6DGcSgwysxGJZa5u7eqxrYlhDlz4JxzUsu21uSioojkpbA1jPSu\nzbcBS9w9dE89iaa5oypY/Nuw25Ga27MHNm+OOgoRqavC3vSeCSVdm3cA1qtbcxHJe3l2PzVs9+at\nEt2AFBE8sV1kZn82s9ZZjU6y5qSTYOPG1Om8chszi4gEwraSGgc0B44huH9R/HdcluKSLGvcOBiB\nLnlq0iTqqEQkzsLewxgEdHP3bYn5JWZ2FbA8K1GJiEjshK1h7CB4ujvZgQSXqEREJA+ErWE8Cbxt\nZg8Aq4CuwE3An7IUl0hGDRoEjZLO9qOPhnnzootHpC4KmzDGAP8CLgc6JV7fR9DHlEjsFaXVhbdr\noGHJhBj2JZVNYZvV7iNIDkoQIiJ5Kmyz2nFmdnJa2clm9lB2whKpvWnToKCgdHr//agjEqnbwt70\nvhT4e1rZQoJLVCKx1KwZtGxZOjVvHnVEInVb2ISxr5x1G1D1YEgiIlJPhE0Ys4F7El2DFPc8ezfw\nXrYCExGReAnbSuoXwFRgnZmtAg4B1gLnZiswEZHYy5PWUcXCtpL60syOA34IfJdgaNZ56oBQRCR/\nhB5xz933uvv77v5S4m/oZGFmw8xsgZkVmdkzacsGmNnnZrbNzKab2SHV+QAiIpIbNRqitQbWAKNJ\ne47DzA4E/grcBrQFFgAv5igmERGphtoM0Rqau08BMLM+QOekRRcC/3T3vyaWjwI2mNkR7r4kF7HV\nZ8cdB598Ujq/TxcQRaQWclXDKJbeDLcnsKh4xt23A0sJuk+XWtq9G3btKp327Ik6IhGpyyqsYRQ3\noa1KNW98pzcpaA58k1ZWALSoxjZFRKKhvqRKpP8edcrWEBxoWI39pb9/K9Aqraw1UFjum6307Z4n\n/0CZsmABHJNUbzM9cilSL1kW/3NXljC6Jb0+G7iYoNfa1QTPYYwguGFdHenf8p8APy2eMbPmwGGJ\ncsmgJk2gadOooxCRuqzChOHuK4tfm9mvgD7uvilRtNjMFhC0anq0qp0kngxvnNhfQzNrSlCDmQL8\n3swuBN4A7gI+rOiGt2oVkinr1sGIEallF1wAJ5wQTTwimZL8PZnp2kbYVlKtCMbw3pRU1ozg8lEY\ndwB3Js0PAUa5+2/N7CLgEeA54P8Bl4TcpkiNbdgA992XWtatmxKGSGXCJowJwP+Y2YMET3kfAvwc\n+HOYN7v7KGBUBcveAXqEjENERCISNmH8hqC56yUEI+6tBcahIVqlDunQAe69N7Vs8mT44INo4pF6\nIM8uk1dnxL3HE5NIndS+fdn7FsuXK2GIhBV2xL0GZvazRF9PHyfKTjWzn2Q3PBERiYuwT3rfDVxD\ncAmquHPANQRNa0VEJA+ETRhXA+e4+wsEo+8BrCD1WQ0REanHwiaMBgRPZSdrTgVPZIuISP0TNmFM\nA/5gZvtBST9To4HXshWYiEjsmeVVPzthE8avgI7AZoKH+LYCXdE9DBGRvBG2We0W4AIz6wB0Ab50\n97VZjUxERGIlbLPaZ8xsoLuvd/f5xcnCzKrsR0pEROqHsJekhgBPm9ktaeVXZDgeERGJqbAJowg4\nAbjEzJ5L9DYrIiJ5JPSY3u7+lZn9CHgKeC/RJbnExD/+AW++mVr29dfRxCKSN9SXVMXcfQdwmZmN\nAOYBTbISlVTb/PkwcmTUUYhIfRY2YYxOnnH3sWa2CFBfUiIieSJss9r7yymbRvBAX62ZWWfgMeBk\nYBfwMvBLd9+bie3nm1694MwzU8vatYsmFhGpPypMGGb2lrsPTLx+r4LV3N1PzUAc/wVsIBhroy3w\nNnADwZgbUk0nnFB23AcRkdqqrIaRPJreUxWsk6k7Pj2BX7j7LmC9mb2ZKBMRkZioMGG4+/NJr5/N\nchxvEdxMnwUcAJwJ3J7lfYqkeOEF+Pjj0vmGDeGhh6KLR+qA4n6k8qS1VGWXpK4hRA3C3Z/OQByj\ngP8BCoCGwLPu/rcMbFcktJkzg6lY48ZKGCLJKrskdQXhLjnVKmGYmRHUMCYRPBzYkuCp8vvcfXja\nuiWvPU8yuohIdVgWe8+t7JLUaVnba6oDgeOB/u6+G9hoZs8SNOUdXtkbRWrrkkvge98rnd+zB371\nq+jiEYmzaj24ByU1gpIU5u77Klk9jA3AWmComT1AUMP4KbAofUXVKiTT+vcPpmK7dilhSN2W/D2Z\n6dpG2N5qDzazKWa2EdiTNO2ubQAefLoLgXMJkscXwE7gptpuW0REMidsDeNxYAfQH5gF9APuIkMP\n7rn7PKBvJrYlIpIzeXbVI2zCOAU4xN23mhnu/mGiFdVc4InshSciInERtnvz4ktQAJvMrD2wDTg4\nK1GJiEjshE0Y8wkepoOgCeyLwBRgQTaCEhGR+Al7SWoIpcnlJuDXQAtAjzWJiOSJsL3Vbk56vZ20\n7s5FRKT+C5UwzKwxcCnQm6BmUczd/WfZCExEJPbUl1S5JgLfI2hGuy5RZmSut1oREYm5sAljEEGz\n2oJsBiMiIvEVNmF8RtDtuBJGTFx3HexNGo/ws8+ii0VE8kN1Wkk9bWbTgPWJMiO4h/Hnit8m2fLM\nM7C71h2ziIiEFzZh/JTgae9WBF2EJFPCEBHJA2ETxi+B3u7+aTaDkZp57DFolPQvedRR0cUiklfy\npHVUMQvTZbiZfUGQMLZmP6Ry9++g7s2TNWlSeklq585gXmpv1y5o2jR43bhxMC9SVxV3b+7uGenn\nPGwN4w/ARDO7n9J7GMWBLM9EICJxs3s3HH98atnll2u8DMlfYRPG+MTfwWnlTjAGt0i9tHBh6ny/\nftHEIRIHVXY+mBhh73Cgqbs3SJsylizM7BIz+8zMtprZUjP7Uaa2LSIitRe2hvERwdCpWWFm/w6M\nBX7i7vPNrBNJw8CK5EqjRrAgrQ/m556Dh9TNpkjVCcPd3cw+BI4keIAvG+4G7nb3+Yl9rs3SfkQq\n1aBB2fsWs2ZFE4vUAepLqlwzgGlm9izwJcG9i+IH956uTQBm1hA4HvhbojXWfsArwC3uXlSbbYuI\nSOaETRg/AlYSjOWdrlYJA+gANAYuSuxnD/A34PbEVKK4iRioia2ISHmSvyczLex4GKdlLYLSJ8fH\nuft6ADP7A+UkjHw1Ywb8/e+pZcn9SImI5ELYGgZm1hY4DzgIWANMdfeNtQ3A3TeZ2Vch163t7uqk\n116DBx+MOgoRqQuSvyczXdsINaa3mZ0ELAOuA74PXA8sNbOTMxTHM8CNZtYukZhuAl7L0LZFRCQD\nwtYwHgZucPe/FBeY2X8kyn+QgThGAwcCS4Ai4EXgdxnYbr3Trx/06ZNa1iBU2heRjMuzqx5hE8YR\nwEtpZX8F/piJINx9D/B/EpNU4rzz1DWFiEQj7G/TLwjG9E72Y2BpZsMREZG4ClvD+AXwupndCKwG\nuhDUOs7JVmAiIhIvYZvVzjWzw4CzCVpJvQpMc/dvsxmciIjER+hmtYkmtBOzGIuIiMRYpQnDzGZU\n8X539/4ZjEdEpO5QX1Ipni+nzIGDCe5rNMt4RCIx9uijMGFCatnChdClSzTxiORSpQnD3Z9Mnjez\nA4ERwLUEz0r8NnuhicTPzp3BlGzfvmhiEcm1sE96tzaz0QTNaDsCx7n7z9w9VJceIiJS91WaMMys\nmZndCiwHjgZOcfch7r4sJ9GJxMDQobBhQ+r03e9GHZVI7lV1D2MFQVK5H1gAdDCzDskruPv0LMUm\nEgv77x9MyRpqJHvJQ1UljOKux6+vZJ1DMxSLiEjdkieto4pVddO7a47iEBGRmFM/pyIiEooShoiI\nhKKEISIiocQqYZjZ4WZWZGbqs0pEJGZilTCA8cB8gu5HRETizay0P6k8EJuEYWaXAJuAd4D8+RcQ\nEakjYpEwzKwVcDdwE0oWIiKxFHo8jCwbDTzp7v8yswovR1lS1c/z7IEZEZEwLIuXyCJPGGZ2LDAA\n6F1cFGE4IiJSgcgTBtAP6AqsTmTGFkBDM+vh7n2SV8yHWsWGDXDbballc+dGE4uI1D3J35OZrm3E\nIWE8AbyQeG3AzQQJpLL+q+qtwkJ44omooxCRUPLgR2yyyBOGu++gtJNDzGwrsMPdv40uKhERSRd5\nwkjn7ndHHUNcHHAAjBmTWnbSSdHEIiISu4QhpVq1guuuizoKqcqRR0KDpAbqgwbBK69EF49Itihh\niNTS7t2Vz4vUF7F4cE9EROJPCUOkBj7/HHbsKJ0mT446IolEnvUlpUtSIjXQtGnqfJMm0cQhkkuq\nYYiISCiqYYhkWGEhfPJJatkhh0DLltHEI5IpShgiGfbee3DMMall06YFzW1F6jJdkhIRkVBUw4jY\nlCmwYkXp/LfqEKVOatECjj46tWzVKti2LZp4JEfUl5Tk0pNPwhtvRB2F1Fa/fmXvWwwaBG+9FU08\nItmgS1IiIhKKahgxMngwdOtWOn/AAdHFIiKSTgkjRq69Fs4+O+ooRETKp0tSIiISSixqGGbWBHiM\nYGzvA4BlwK3u/makgWXY7Nlwww2pZcktpESkjinuRypPWkvFImEQxLEaONXdV5vZ2cBLZvY9d18V\ncWwZs3UrfPxx1FFIFP7jP1L7n+raFebPjywckRqJRcJw9+3A3Unzr5vZCuA4oN4kDMlfBQWp861a\nRROHSG3EImGkM7MOwBHAJ1WtW1edfDI89lhqWdeukYQiIhJK7BKGmTUGngeedfclactKXnsdv2bY\nsiV8//tRRyHZ9Je/wK5dpfMrVsCJJ0YXj+QHy+L4HLFKGGbWAJgIFAHDIg5HpFbatEmdT78sJVLX\nxCZhWJAWnwLaAWe5+970depareKee2DDhtL5lSsjC0ViZs0a6N8/teyGG+Dii6OJR2ooht9Jyd+T\nma5txCZhEDSrPQo43d13Rh1MJjz7LCxbFnUUEkdFRTBjRmrZeedFE4tIWLF4cM/MugA/A3oB68ys\nMDFdGnFoIiKSEIsaRuJZi1gkr2wZORLatSud79IlulgkGgcdBO+8k1o2fjxMnhxNPCLVFYuEkQ+u\nvhq6d486ColSs2Zl71u8+mo0sYjURL3+VS8iIpmjGkaGTJgAY8aklq3SM+pSDf/1X2UvT82YAQ0b\nRhOPhKC+pKQmNm2CJUuqXk+kIitWqDNKiTddkhIRkVBUw8iCK66A225LLTvkkGhikXi78Ua48MLU\nstNOy5srHFLHKGFkwQEHwJFHRh2F1AWHHRZMycyUMCSedElKRERCUQ2jBvbtg2FpXSMuWhRNLFK/\n9e1b2hAH4Pzz4ZZbootH0uRZVVAJowbcy45lIZIN77+fOq8u8SVKShghFBXBnj2l83vL9KMrIlL/\nKWGEcNNcli3qAAAIKUlEQVRN8PjjFS8fPz51Xr8CpTbefTd1fsoUeOCBaGIRSaaEUUsNGgTjGIhk\nyimnpM5/9FE0cYikU8KopiZNoHHj0vkGamcmEXjrLSgsLH9Zu3bQr19u45H8EIuEYWYHEIy29+/A\nBuBWd38h2qjK99BDMHRo1FFIvlq2DF5+GX7848rXmzQpdf6kk+Dgg7MXV95SX1KRGE8wjnd7oDfw\nupktcvdPs73jKVOC4TLL07YtXH55tiMQCe/tt4OpKukJZfJkuOCC7MQk+SPyhGFmzYELgZ7uvh2Y\nY2Z/A64Abs32/seNKztUZrJrroHdu7MdRXjFY/TWtfHNs0HHItCjB3z2WfHDGuUfiwsvTO2epkMH\nmD8/+7FFQedF9kSeMIAjgD3uvjSpbBFwWqZ3tGwZvPFGalllyQJgZ70YXVzqsm7d4KKLKl7+4IOl\nySB5vb/+NXW91atTX597bvDaHaZNC5JKcc1kyxZYty6olRTfsysogBYtoFOn1O22aAGN4vBNIlln\nUWdhM+sLvOTunZLKrgUuc/d/S8w7VP6LYfx4uO++1LIvvwz+NmkS/N21q/JYevSAjh0rTyKPPhrt\nPQz9eiqlY1GqvGNxwQXwyivZ33dxAlm7Fr7znWBkwXQPP5y7S2I5PS9ifg8j6VhYFauG217U/9nM\nrDcw292bJ5XdDJzq7ucl5uP5ryEiUgdkKmHEoVHoEqCRmSWPeN0L+GdE8YiISDkir2EAmNkLBHfr\n/jdwHDAVOMndP4s0MBERKRGHGgbADcD+wNfAc8D1ShYiIvESixqGiIjEX1xqGCIiEnORJgwzO8DM\nppjZVjNbaWaXVrDe42ZWmDQVmVlB0vKZZrYjaXmdu5wV9lgk1r3DzL40s81mNsPMjq7JduIqg8ci\nb84LM2tqZg+a2Roz22hm482sUXW3E2cZPBZ1+rwws2FmtiDxPfhMFeveZGZrzWyLmT1lZk2SllX/\nnHD3yCbghcTUDDgF2AwcHeJ9zwBPJs3PAP5XlJ8lV8cCOA9YA3QlSPhjgL/X9pjGacrgscin8+Iu\nYBbQBjgQeB8YlafnRVXHok6fF8AFwGDgUeCZStYbCKwDeiSOxQzg3tqcE1F+6ObATqB7UtmE5A9U\nyfsKgL5pJ8A1Uf9D5uJYEHSX8mLSfE9gR22OaZymTB2LPDwvPgAuTpq/FFidp+dFhceiPpwXSZ9j\ndBUJ4/8C9yTN/xuwtjbnRJSXpCrqEqRnFe+7CPja3d9LK7/XzL4xs9lmVtc6d67OsXgHOMnMDjez\nxsBPgWk12E5cZepYFMuX8wIg+eGsBkBnM2tZg+3EUaaORbG6fF4Uq+phvKMJjlGxj4AOZtaWGp4T\nUfYA04KgppCsEGhZzrrJfgr8Oa1sOPAJsIvg18RrZnasuy/PRKA5EPpYuPt8M5sALAb2AquBAdXd\nToxl6lhAHp0XwJvAL8xsBsH/658TPNvUrJrbiatMHYtC6v55UayqJq4tgC1J88XHryU1PCeirGFs\nBVqllbUmCLpcZnYI0I+0hOHu8919m7vvdvc/A3OAszIcbzaFPhZmNozgS7Ez0BT4LTDdzPavznZi\nLFPHIq/OC+B3wD+AD4HZwBSCX5Drq7mduMrUsagP50WxqmoY6cesdeJvYTnLipdXek5EmTBq0iXI\nFQT9Tq3MZmARqM6xGAS84O7/cvd97j4BaEtwY6s+dLOSqWNRH4Q+Fu5e5O43untnd+8ObAQWVHc7\nMZapY1GfVFXD+AQ4Nmm+F7De3TdR03Mi4ps2LxDcmGkG/IjgLn2PStZfDFyVVtaaoDXAfgTVz8sJ\nsmf3bMYe1bEgaAn0HsFgUw0Ikmgh0KomxzSOUyaORR6eFwclJgNOJLg8d3p1txPnKRPHoj6cF0DD\nRPz3ElxtaQo0LGe9gcBagh9QbYGZwJjanBNRf/C2BNXFrcBK4JJE+SGJ//idk9Y9KVHWPG0bBwLz\nCa7HbQLmAgOi/kfN1rFI/OM+SdBcbgvBL6czqtpOXZoycSzy8LzoC6wAtgGfAZeG2U5dmjJxLOrD\neQGMAvalTXdW8L15U9L/j6eAxrU5J9Q1iIiIhKKuQUREJBQlDBERCUUJQ0REQlHCEBGRUJQwREQk\nFCUMEREJRQlDRERCUcIQEZFQlDBERCQUJQwREQlFCUMkITGu8c1m9lFirOenzKyDmU1LjIn8tpm1\nSaw7KTFW8mYzm5U2lvhZZvaJmRWY2Vdm9utE+YFmNtXMNpnZt2b2rplV1UW1SGwoYYiUcuBCgjE2\njgTOIRjBbwSlPeL+PLHuG0B3oB2wEHg+aTtPAT9z91YEI5hNT5T/GviSoAO89sCtrs7cpA5RwhBJ\nNc7dv3H3fxF0nf6+uy9y950EPXv2BnD3ZzwxCA9wN9AraQjQXUBPM2vl7lvc/R9J5Z2Aru6+193n\n5PSTidSSEoZIqvVJr3ekzRcBLcysgZmNNbOlZraFoCttJ6g5QDDu/FnASjObaWYnJsp/DywF/tvM\nlpnZ8Kx+EpEMU8IQqVx59xguA84jGEehNXBoYj0DcPcF7n4+weWqV4CXEuVb3f1mdz8s8f5fmVn/\nHHwGkYxQwhCpvpbATmCjmTUnGPkPADNrbGaXm1lrd99LMKDN3sSyc8yse+JGd0GifG/uwxepGSUM\nkcp52msnGBZzFbCGYAzk99PWGwKsSFyu+hnBMKAQ3CR/myCJzAXGu/usrEYvkkEacU9EREJRDUNE\nREJRwhARkVCUMEREJBQlDBERCUUJQ0REQlHCEBGRUJQwREQkFCUMEREJRQlDRERC+f9Az6HWL0pD\nmQAAAABJRU5ErkJggg==\n", "text": [ "" ] }, { "metadata": {}, "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAYUAAAEXCAYAAABCjVgAAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xe8VNW5//HPQ28CR6UI6MFG+aEi6i+iiSWisWDsUdGL\n5UZ9qVExJna9gEZNNInGckk0KgG9tp/yUzEQvRbEaESCUSlXBOmKoCDlUKQ89489ZximnT3n7Nkz\n55zv+/XaL2fWLutZzmGeWbusZe6OiIgIQJNSByAiIuVDSUFERJKUFEREJElJQUREkpQUREQkSUlB\nRESSlBRERCQptqRgZmvNbE3KstnM7o+rfhERqVmzuCpy93bVr82sLbAUeDau+kVEpGalOn10BvCV\nu79TovpFRCSLUiWF84ExJapbRERysLjHPjKzSmAusKe7L4i1chERySu2awophgKTsyUEM9PofCIi\nteDuFsVxSnH66DzgLyWoV0REahBrT8HMDgW6Ac/l264hDudtFiTxhtg2UPvqO7Wv/qpuW1Ti7imc\nBzzv7lUx1ysiIiHEfqE5n+prCuUUU1Qa8i8VUPvqO7Wv/kppWyRdhlJcaG6UGuIfYyq1r35T+6Sa\nxj4SEZEkJQUREUlSUhARkSQlBRERSVJSEBGRJCUFERFJUlIQEZEkJQUREUlSUhARkSQlBRERSdIw\nFyJSZxs2wHff5d+mTRtopm+csqePSKSRGzcOqvKMW9yqFZxxRv5jDB8Od9+df5u33oIjjig4PImZ\nkoJII3f11bBwYe71XbvWnBSk4VBSEJG8li6FJ57ILN911+y//Fu0gJYtg9dVVbB1a3Hjk2gpKYhI\n0sknQ7t28OWX8MYb28qHDs3c9pRTsieF226D668PXh9xBLz9dnFileJQUhCRpD/8ASorg6TQrVup\no5FSUFIQkQytWsG552aWL1wIkyfHH4/ER0lBRDJUVGS/jjBunJJCQ6ekICKxePVVWLQos3zIEGja\nNP54JDslBRGJxZ13Zi//yU+UFMqJkoKI1MqECdsuRq9eXdpY4jJ16lTGjh3LQQcdxN///neuvfZa\n9txzz1ptO2bMGBYtWkTXrl1ZtWoV11xzTVzNyCv2pGBmZwPDgV2BpcAF7v5O3HGISN1s3BjcpZTP\nMccEzzOk+6//AvfixFUXo0aN4rLLLsu6buPGjZxxxhm8//77dOnShb59+zJkyBCmTJlS8LaPPfYY\nc+bM4c4772TBggX06tWLCy64gB133LGo7Qsj1qRgZscAvwbOdPcpZrYLYHHGICLxueWW7OXPPVfz\nWEmlsGzZspzr3n77bdq1a0eXLl0AOPDAA5k1axbz58+nZ8+eobft1q0b119/PdOmTQOgsrKSTz/9\ntCwSAsTfUxgJjHT3KQDuXsPvDBEpJ8cdB0uW5N9mhx3iiSVu8+fPZ6eddkq+NzMqKiqYMWNGRlLI\nt+38+fP55ptvmDdvHu+99x5Tp07lmGOOyThGqcSWFMysKXAg8KKZfQa0Av4/cK27b4grDpHGavJk\nWLAgs3zt2vDHaN06WLJZt24dAwcO5OOPP866/tNPP+Wmm26isrKSzZsN+Bq4B+hcY72p+5oZX3/9\nNffccw+dO2fft6ZYauPrr7+mTZs225W1atWKNWvWFLTtqlWrAGjWrBlnnnkmJ554It27d2f69Ol0\n7949snhrK86eQhegOXA68ANgM/AicEtiSTLbdkbJy/HEo0g99OCD8OyzxTn2Bx98wOWXX86MGTOy\nrl+1ahWDBg3i7rvv5pxzzuGhh+C77+4CBrFp0zRatmye89jp+wLcddddDBo0iGnTptG8+fb71hRL\ntWXLlnHfffdt9x3zzjvvsGHDtt+o7dq14+abbwagY8eOGd9Ha9euZeedd844dr5tN27cCASnlADa\ntGlD69atefnll7n00kvzxlwt9Tsycu4eywJUAFuBoSllpwHTUt57+iIi0TjzTPfg8m7uZf78wo45\na9YsHzx4sF9wwQU+cOBAb9KkSdbtbrrpJu/SpYtv2bLF3d1btHCHFQ7N/YEHRuWtI31fd/cVK1Z4\n8+bNfdSobfuGjSWfESNG5Fz3xhtveP/+/ZPvN23a5C1btvTZs2cXtO3ixYvdzLyqqiq5vlu3bv7g\ngw+GjjPHd2Uk39Wxzbzm7iuBxSG3TU0mIhKx738/GMYifWnbtrDj9OnTh/Hjx/P444/Tu3fvnP9m\nn3vuOQ4++GCaNEn9yqkA+vD888/lrSPbvhUVFfTp04fnntu2b9hYauuwww5j+fLlLEo8gTdp0iT6\n9evH3nvvDcDrr7/ORx99VOO23bt358gjj+Sdd4KbLpcvX87atWs59dRTQ8dSzO/IuC80Pw5caWYT\nCU4f/Rx4OeYYRBq9K6+Es86Kp641a9YwZ84cjj322Cxrd+HDD9+v1b677LIL77+fe9+oNWvWjLFj\nx3LHHXdw6KGH8uabb/LMM88k1z/00EMccMAB9O/fv8Ztx4wZw8iRI5k5cyazZs3ipZdeoluZjEAY\nd1K4HdgZmA1sAJ4B7og5BhGJ0YLE1e327dtnWduW1atXs2nTpoxrAzXt27Zt/n1ro3Wuq+gJRx11\nFEcddRQA55133nbrXnjhhdDb9ujRg0ceeaSu4RZFbKePANx9s7v/zN0r3H0Xd7/a3cvwbmURicrq\nxOPOLVq0yLI2OF/17bffFrxv27b5962N6667LrJj1VexJgURaXyaJgY2yn7HzCYAtmzZUvC+mzbl\n31dqR0lBRIqqU6dOedZWYWbskOOJt3z7VlXl31dqR0lBRIqqS5cumBkrV67MsraKjh07Jk8FFbJv\nVVX+faV2lBREpKjatm3LgAEDWLhwYZa1c9hvv/1rte+cOXPYf//c+0Zh6tSpDBs2jLFjx3LppZcy\nd+7cWm3bo0cPKioqtlt+9rOfFTX22tLQ2SISmVxP2g4ePJhHH300rXQusJjTTrtxu9LPPvuMXXfd\nlVatWuXcd+7cuSxevJgbb9x+3zCxhBXVqKhLly7luuuu46STTkrG9Nvf/pbbb7+9TvEVi3oKIhKJ\nDRs24O6sX78+Y91ll11GVVUVT2w3x+cDQD/+/d8vTpa8/fbb9O7dm1NOOSXvvg888AD9+vXj4osv\nJpt8sYSVb6TTsNvOmzcPM2Po0KH07NmTyspK3n//fc455xw6duxY69iKST0FEam15cuXM3ToUJYs\nWcKMGTMwMyorK9lnn3246KKLkmMVde3albfeeoubb76ZadOmsXnzGmAlMJFmzbZ9DXXt2pXOnTuz\n1157bVeWuu+aNWtYuXIlEyduv2/YWMKKYlTUmTNnMnjw4GT5kiVLmDlzJmeeeWZBscRJSUFEaq1T\np05MnDgx1Lb9+/dn/PjxAIkB8TK36dWrF0uXLs27bxSxhBHVqKipbr75ZkaMGBFZjMWgpCAiJfXU\nU5D+QHJFBZxwQnHqe/DBB/NeMD7ooIM499xz6dChQySjolZbtmwZb775JqNHj65bA4pMSUFESurC\nCzPL9t+/eEnhiiuuCLVd3759efjhh5PvN2/ezIoVK6isrMzYtk+fPvzpT3/Ku+2ECRO2O8VUrnSh\nWUQki6hGRa02ffr0jFNM5Ug9BRGJ3dlnw+bN25etWAERXhKosyhHRQXo0KEDvXv3jrsZBbNymrPA\nzIKZdsooJpGG4qyzts289vTT8Q2dHdaHH8IBBwSv998/eC81q372wd0jmY5Np49ERCRJSUFERJKU\nFEREJElJQUREkpQUREQkSUlBRESSlBRERCRJSUFERJKUFEREJCnWpGBmb5nZejNbk1hmxVm/iIjk\nF2rsIzP7pbv/Nkv5Ne7++wLqc+Bn7v5YAfuISBZr18Knn+bfpnNn2HXXeOKRhiHsgHjDgYykANwK\nFJIUACIZn0OksfvkEzj00PzbXH013HtvPPFIw5A3KZjZUQRf4k0Tr1PtCayuRZ13mdmvgU+Bm919\nUi2OISIiRVBTT+ExglM+LYFHU8od+Aq4ssD6rgdmAN8BQ4CXzWx/d/88daPqUf9AI6aKhNGmDfTp\nE7z+6itYsmTb62nTgtcrV5YmNole6ndk1PImBXfvmQhgrLsPrWtl7j4l5e0YMxsCnAA8WNdjizRm\n++0H770XvL73XrjmmuD1U08Fi0hYoa4ppCYEM2uStm5r1EGpdyAiklvqd2TUvYawdx8dSPBrvj/Q\nKjU2oGnIY3QABgKTgM3AWcBhFH4KSkTy6NwZBgzIv03HjvHEIvVP2LuP/gK8BPwUWFfLupoDtwN9\ngC3ALOBkd59Ty+OJSBbnnhssIrURNinsRnCnUK3P67j718D3aru/iIgUX9gnmscBxxYzEBERKb2w\nPYXWwDgzm0xwK2o1d/fzog9LRERKIWxSmJlY0uk2IRGJ3Lx5cPbZmeVHHw0XXRR/PI1J2FtSRxQ5\nDhGRpFWr4JlnMssrKpQUii3sLanpQ1wkufsb0YUjIiKlFPb0UfVwF9U6EQx9sQjYI+qgRKTx6dkz\n+9PXr70Gj2lc5diEPX3UM/W9mTUFbgHWFiEmEWmEKiqyX0dYuVJJIU61mmTH3bcAdwLXRRuOiIiU\nUl1mXjuG4MlkERFpIMJeaF6UVtSGYAykyyOPSERESibsheb0YbOrgNnuvirieEREpITCXmh+C5LD\nZncBvirGkNkiIlJaoa4pmFl7MxsDbACWABvMbExiOGwREWkgwl5ofgBoC+xDcD2h+r8PFCkuEREp\ngbDXFI4D9nD3qsT72WZ2AfB57l1ERKS+CdtTWE/wFHOqnQlOJ4mISAMRtqfwZ+A1M/sdsADoCfwc\neKRIcYmISAmETQp3Al8A5wK7JF7/hmBMJBERaSDC3pK6lSABKAmIiDRgYW9JfcDMDk0rO9TM7itO\nWCIiUgphTx8NAX6ZVjYNeBG4OtKIRGQ77jByZGb5ovTBZ0QiEDYpbCWzV9EEsGjDEZFssiUFkWII\nmxTeAX5lZte6+9bEfAojgcmFVmhmewOfAM+5e/qYSiIiOX32GTz9dGb5wQfD7rvHH09DFDYpDAPG\nA0vNbAGwG/Al8ONa1PkQMIXtZ3ITkZCGD88s69Ej/jhK4fXXgyXd6NFKClEJe/fRIjM7APgesCvB\nNJzvFzoonpmdDawEZgJ7FRiriAAjRpQ6AmnIwvYUqmdbey+xFMzM2hOccvohcEltjiEijc9ee8FZ\nZ2WW/+MfsGBB/PE0dKGTQgRuB/7s7l+YWd5TR2bbrl+76yyTSGN2zDHBku7882HMmPjjKQep35FR\niyUpmNn+wCBgQHVRHPWKiEhh4uopHEEwXtLCRIZrBzQ1s77uflD6xuodiIjklvodGXWvIWdSSMyy\nVqOQF5sfBp6qPjTBg3A9gUvD1CEiIvHI11PYnPbeyTzt40DTmipx9/UEw28DYGZrgfXu/k3IOEVE\nJAb5ksIeKa8HA2cQjJa6kOA5hRuA52tTqbvr+UwRkTKUMym4+/zq12Z2DXCQu69MFH1qZlOBqcB/\nFjVCERGJTdiZ19oTzMmcqg3QIdpwRESklMLeffQX4L/N7F6Cp5l3A64CGuldwiLFceutMHfu9mW6\nGU/iFDYpXAfMAc4mmHntS+ABNB2nSKRefRWmTCl1FNKYFTLz2h8Ti4iINFChkkLimYWLCHoKndx9\nXzM7HOjq7s8WM0CRxmrkyGDcH5E4hT19NBL4EXAf23oLSxLvlRREiuDYY4N5AkTiFPbuowuBE939\nKYJZ2ADmsf2zDCIiUs+FTQpNgLVpZW2BNdGGIyIipRQ2KUwAfm9mrSB5jeF24OViBSYiIvELmxSu\nAboC3xI8yLaWYEC7G4oTloiIlELYW1JXAaeaWRegEljk7l8WNTIREYldqJ6CmT1uZse6+1fuPqU6\nIZiZxj0SEWlAwp4++jfgMTO7Nq18aMTxiIhICYVNChuAg4GzzewJM2tZxJhERKREwiYF3H0x8IPE\nPpPNrEfRohIRkZIoaI7mxAxq55jZDcD7QIuiRCUiUoDhw+H++zPLn38eevaMPZx6LWxSuD31jbv/\n2sw+As6MPiQRkcIsWBAs6TZujD+W+i7sLal3ZymbQPBQm4iINBA5k4KZ/c3dj028npxjM3f3w4sS\nmYhIHiNGwFVXZZafdhosXBh7OA1Gvp5C6qxqj+bYRnNCidTCPffA009nls+aFX8s9dXuuwdLulat\n4o+lIcmZFNz9yZTXo2OJRqSRWLwYpk0rdRQimfKdPvopIXoC7v5Y2MrM7AlgEMEIq18Dj7r7HWH3\nFxGR4sp3+mgo4U4PhU4KwF3ARe6+wcx6A5PM7J/uPrGAY4g0KL/4BQwZklnep0/8sYjkO310ZNSV\nufuMtKLNwLKo6xGpT3bdFQ48sNRRiARCP9FczQJNqpda7P+fZlYFzAB+5e46syoiUibCjpLa3czG\nmdkKgl/31cumQit098uBdsDRwK/M7HtZ6ksuIiKyvWJ+R4b9pf9HggRwFMEEOwcALwKX1aZSD7wF\nPAdkOZsqIiKlEHaYi+8Du7n7WjPD3f+VuDvpXeDhOtTfHPgmvdBdjz+IiOSS+h0ZdW8hbE+h+nQR\nwEoz6wxUAd3DVmRmnczsbDNra2ZNzexY4CcEPQ4RESkDYXsKU4DjgXHA34BngPXA1ALqcuBSYBRg\nwGxgqLt/UMAxRESkiMImhX9jW6/i58AvCC4W3xe2Inf/GjiykOBERCReYUdJ/Tbl9TrShtIWEZGG\nIVRSMLPmBHcJDSDoIVRzd7+kGIGJiEj8wp4+GgvsSzB/wtJEmaFRUkVEGpSwSeE4gltSVxczGBER\nKa2wt6TOAnYsZiAiIlJ6hdx99JiZTQC+SpQZwTWFMbl3E2nc/vlPGJPlX8ikSfHHIhJG2KRwPsFT\nze0Jnk9IpaQgksPs2XD//aWOQiS8sEnhamCAu88sZjAiIlJaYZPCV4CmwhapgwED4PzzM8sPPzz+\nWBqL776DjRszy1u2jD+W+iJsUvg9MNbM7mbbNQUA3P3zyKMSaYB694Zhw0odReOy336ZZR07wsqV\n8cdSX4RNCg8l/ntyWrkDTaMLR0RESqnGpGDBuKx7AwvdveBJdURE4tS8ObRosX2ZO2zSt1coYXsK\nHwM7FDMQEZEoTJ+eWbZyJeyoJ61CqfHhNQ9mc/gX0Lv44YiISCmF7Sm8CUwws9HAIoJrCdUPrz1W\npNhERCRmYZPCD4D5wBFZ1ikpiIg0EGHnUziyyHGIiEgZCNtTwMwqgJOAbsASYLy7ryhWYCIiEr+w\nk+wcArwC/A+wAPgxcJ+Zneju7xYxPpF6YelSuPDCzPIvvog/FpG6CNtT+ANwubs/XV1gZmclyv9v\nMQITqU/Wr4eJE0sdhUjdhZ1PoRfwbFrZ8wQPtYmISAMRtqfwGcEczU+mlP0EmBN5RCL1XOfOMHp0\nZvkuu8QeikjBwiaFYcArZnYlwWiplQS9hxPDVmRmLYBRwCCCWdzmAje6uzrd0qC0aQPHH1/qKERq\nJ+wtqe+a2Z7AYIK7j14CJrj7NwXWtRA43N0Xmtlg4Fkz29fdFxQauIiIRC/0LamJ20/H1rYid18H\njEx5/4qZzQMOILijSURESixvUjCzN2vY3939qNpUbGZdCE5BzajN/iIiEr2aegpPZilzoDvBdYY2\ntanUzJonjj3a3WdnWb+tMvfaVCEi0mClfkdGLW9ScPc/pwWyM3ADcDHwDHBboRWaWROC01AbgCsK\n3V9ERIon7BPNHYBfAlcC44ED3H1uoZUlJux5FOgEnODuW7Jtp96BiEhuqd+RUfcaarqm0IbgNNEv\ngbeA77t7Xa4BjAL6AEe7e5bptEVEpJRq6inMI3jq+W5gKtAlcYE4yd3fCFORmVUClxCcNlqakt0u\ncfenCglapJSGDIFvv92+bN260sQihduyBT75JLO8XTvYfff44yk3lu9UjZnNT7zMuZG7R/a/0cw8\nccyoDikSuS5dYNmy3Ot79oR582ILR0IIMx3nccfBhAnxxBOl6h/Y7h7JeaSaLjT3jKISERGpH0I/\nvCYimZ58Eioqti9r3bo0sUhuTZrAPvtklq9ZAwv06Ox2lBRE6mDQoOB0kpS3Dh2yX0eYMAFOOCH+\neMpZ2KGzRUSkEVBSEBGRJCUFERFJUlIQEZEkJQUREUlSUhARkSQlBRERSVJSEBGRJCUFERFJUlIQ\nEZEkJQUREUnS2EciIsCKFTBpUmZ5ZWUwHHpjoaQgksO778KXX2aWb9gQfyxSfFOmwJFHZpaPGAHD\nh8cdTekoKYjk8JvfwEsvlToKkXgpKYhIo1VRAYcfnlk+fz4sXBh7OGVBSUEkhIEDoVu3zPKWLeOP\nRaIzcGD26wjDh8Ntt8UfTzlQUhAJ4YYb4OSTSx2FSPHpllQREUlSUhARkaTYkoKZXWFmU81sg5k9\nHle9IiISXpzXFJYAtwPHAq1jrFdEREKKLSm4+zgAMzsI6BFXvSIiEl4p7j6yEtQpktOHH8IHH2SW\nz58feygiJVeKpOA1bWC2LW+417i5SJ1MnAg33VTqKETCS/2OjFop7j5ST0FEpEyVZU9BvQMplf79\n4XvfyyyvrIw/FikPTzwRDJaX7rbb4MAD448Htv+OjLrXEFtSMLOmQPNEnU3NrCWw2d23xBWDSE2O\nPx7uuqvUUUg5mTMnWNJddVX8scQhzp7CrcB/pLz/N2AE0EhHGJG4LV8Okydnlk+fHn8sIuUqzltS\nRxAkAZGSmD4dTj+91FFIfXDuudlPI956a3C3WkOmAfFERNL06hUs6e6/P/5Y4qakII3STjtlH0e/\nX7/4YxEpJ0oK0ijtuy+88EKpoxApPxolVUREkpQUREQkSaePpF57/XVYtSr3+iZN4JRT4otHGo+V\nK+GrrzLLO3UK/u7qKyUFqdeuuQY+/jj3+latYP36+OKRxmPIkOzlq1ZB+/bxxhIlJQVp0Ny3/Zpb\nsaK0sYjUB0oK0mD88IfQsSNs2QIvvRSUbdwIXbuWNi5pODp2hM6dM8uXLYs/lmJRUpAG4957gwHt\n1q+HNm1KHY00RM88k728fXtYsybeWIpFSUEapGy/5lJVVMQTh0h9o6QgDU7r1tnvChGRmtXjG6dE\nRCRqSgoiIpKkpCAiIkm6piBlZfx4uOOO7Ou2bg2mRdxjDzj11KDsiy/ii00kjOXLYcOG7cuaNg1G\n5q0PlBSkrCxfDv/4R/5tPv8cfve7eOIRKdRee2WW9e0LM2fGH0tt6PSRiIgkqacgJbFoEVx2Wfby\naiecALfcErx+9VWYOxf22y/3MXfZJdoYRcLaeedgnK1UmzcHg+bVN0oKUnT33hsMN5Fq1ix45ZX8\n+3XuDIccEryu/q9IOfr888yymTO3zeS3cCGceGLmNkcfDVdfXdzYCqWkIEU3YgSsXl3qKERKp6oq\n+4+gmp68L4VYk4KZ7Qg8ChwDfA3c6O5PxRmDlKeXX84s69Ej/jhE4jR/Powbl1m+//6w++6xhwOA\nuXt8lZlVJ4CfAgOAV4BD3X1mYr0DxBmTROfbb+G77zLL99wT1q4NXg8blnnutVMn+MUvih+fSJzW\nrIFJkzLLX3sN7r8//75/+hNcckm4eswMAHe3AkPMKraegpm1BU4D+rn7OuDvZvYiMBS4Ma44SiXl\ng4ulvtWrM++VTrfjjtCsWTDU9OzZubdbtCiYl2C33YL3s2bBZ59Br17btjnttOq/x/ztGzkSOnSo\nOf5yE/fnFze1L3o77JD9OkK5D7Md5+mjXsBmd5+TUvYRcGSMMdRb7sGS7t134YUXwCxYqoW5j/+u\nu6B37+AOiZ/+NLpYRSS3ysrsU8R++CEsWBB/POniTArtgPTLjWuAHWKMIadhw+Cvf80sn5NIYdke\nSKleV/1r4NVXg9MnlZXQvHn2esxg4MDgdfVDWtmGca6+lW3AAGjZMviDSb+Dp65uLHL/rFOnzDKL\npIMrUn8NGhQs6S65BB55JP540sV2TcHMBgDvuHvblLJfAoe7+0mJ9w2z7yoiUmRRXVOI84nm2UAz\nM0v9zd0fmB5jDCIikkcp7j5y4CLgAGA8cIi7z4otCBERySnusY8uB1oDy4AngEuVEEREykesPQUR\nESlvGiVVRESSYk8KZtbXzN4ws2/N7DMzy3LHbnLbPcxsvJmtNrPlZvabOGOtjbDtM7M/mtmalGWD\nmZX1CEEFfna3mtmixLZvmtn/iTPW2ijgs2tpZvea2RIzW2FmD5lZ2Y0jZmZXmNnUxN/W42nrBpnZ\n/5hZVaLNu+U5zo5mNs7M1prZfDMbUvzo84uwbTmPU0pRtM/MWpjZo4nPbLWZfWhmx9VUd6xJIfEP\n50XgJaACuAR4wsz2zrJtC+A14L+BLkB3gusQZauQ9rn7pe6+Q/UCPAU8G2vABSjwszsJuBQ4DNgR\neA8YG1+0hSukfcANBDdK9CN4KPMA4JaYQi3EEuB24LHUQjPbGXgeuJmgrVOBZ/Ic5yFgA9AZOBcY\nVQZJPqq2ZT1OGYiifc2AhQS3/bcn+Bt91swq89bs7rEtwD7AmrSyvwG3Zdn2EmBSnPHF2b60bdoS\nPNh3WKnbENFndyPwTMr7fsD6UrchwvZ9AJyR8n4IsLDUbcjTttuBx1PeX0LwzFD1+zbAOqBXln3b\nAhuBvVLK/gLcVep21bVt+Y5TLktU7UvZ/iPg1HzblMM1hSYE/yDTDQQWmNlfE6eO3jSzbNuVu1zt\nS3U6sMzdJ8cQT5Ryte114BAz29vMmgPnAxNijSwa+T47S9uuh5mVxdP5WaQ/1NSP4MsBAA/GIptD\n9rbmGp6mX9RB1lJd2pbvOOUiqvZhZl0IPs8Z+baLOyl8Ciwzs2vNrLmZ/Qg4nOA21XQ9gLOBPwC7\nEIyo+mLiS6ZcFdK+VOcDY4oeXd2Ebpu7TyH4Nfkpwa+Y04Fr4gy2Fgr57CYCw8xsZzPrClxF8PxN\nm/jCLUj6LYbVPdNUqwmGoklX1sPTULe25TtOuYikfYnvzSeB0e6eZ/jLmJOCu28CTgEGA18CPyc4\nj744y+brgMnu/jd33+zuvwV2AvrEFW+hCmwfAImLREdQ5kmhkLaZ2RXAIILE3hK4DXjDzGpKjiVT\n4Gd3B/DLWNYcAAADz0lEQVQh8C/gHWAcwa/pr+KJtmDpvzbXAu3TyjoQfNmnK2TbUqhL2/Idp1zU\nuX1m1oTgmt4G4IqaKoz99JG7f+LuR7r7zu5+PLAnMCXLph+nvjGrH0OpFdC+akMJzhHOjyXAOiig\nbccBT7n7F+6+1d3/QnBRrG+c8RYqbPvcfYO7X+nuPdx9L2AFwQW/cpX+a3MGwRAzQHJY+z3Jflqh\n3IenqUvb8h2nXNSpfYnvzUeBTsDp7r6lpgpLcUvqvmbWyszaWDAgXhdgdJZNnwAGJm6/agpcDSwH\nyvoJ6ALaV+28GtaXjQLa9jFwppl1NrMmZjaU4E6IOVm2LRth22dm3RKLmdlAgrs6hsccbo3MrKmZ\ntSL4f9/UgltpmxL0bPYxs9MS64cD/8p2WsHdq4AXgNsS/19+APyYEt9NFkXbajhOSUXVPmAUwdmV\nk9w93DjLJbiafjfBL6s1BNcJ9kiU75Yo65Gy7anAZ8Aq4A2gbynvBChC+w5JlLUtddxRto3g3Pqf\ngaWJz24q8KNSxx9h+w4D5gFVBD9ShpQ69hztGQFsTVv+I7FuUCL2dYl/W7ul7HcT8NeU9xUEX0Zr\ngfnA2Q2obTmPU9/bB1Qm9luX+PutXvL+vWqYCxERSSqHW1JFRKRMKCmIiEiSkoKIiCQpKYiISJKS\ngoiIJCkpiIhIkpKCiIgkKSmIiEiSkoKIiCQpKYiISJKSgjR6ZnaDmc1JzGM7wxJzMycGJftdYpKn\nzxPz5m5NDEWMmXVIzIH7hZktNrPbq9eJ1FdlN9m4SAnMAX7g7kvN7EyCuZn3Iphf4TiCoYrXAf+P\n7YcyHk0w6N+eBJOcjAcWAQ/HF7pItDQgnkgaM/uQYEjiYcDT7v5IonwQ8BrBj6lOwAKgo7tvSKwf\nAlzs7keVJHCRCKinII2emZ1HMNNaz0RRO2BnoBvBL/9qqbOwVQLNgS9T5n9qAiwsZqwixaakII2a\nmVUSnO45CnjP3T3RUzCCaTl3Tdk89fUiYCOwk7tvjStekWLTRTFp7NoSXCf4GmhiZhcC+yTWPQsM\nS8yy1hG4PrEt7v4l8CrwezPbITHD3J5mdnj8TRCJjpKCNGruPhP4HfAewUXjfYB3CL78HyH44v8Y\n+CfBbGxbUnoG5wEtgJkEM7Y9B3SNM36RqOlCs0hIZnY8MMrde5Y6FpFiUU9BJAcza2VmJ5hZMzPr\nTnBH0guljkukmNRTEMnBzFoDk4A+wHqC5xCGufvakgYmUkRKCiIikqTTRyIikqSkICIiSUoKIiKS\npKQgIiJJSgoiIpKkpCAiIkn/C7Fot5UxisfWAAAAAElFTkSuQmCC\n", "text": [ "" ] } ], "prompt_number": 4 }, { "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": {} } ] }