cdecker/Untitled0 (1).ipynb

## Untitled0 (1).ipynb
{
 "metadata": {
  "name": "Untitled0"
 },
 "nbformat": 3,
 "nbformat_minor": 0,
 "worksheets": [
  {
   "cells": [
    {
     "cell_type": "code",
     "collapsed": false,
     "input": "from sympy.mpmath import *\nlibmp.BACKEND",
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "pyout",
       "prompt_number": 1,
       "text": "'gmpy'"
      }
     ],
     "prompt_number": 1
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": "mp.dps = 50\n\nimport operator\ndef mult(s):\n    return reduce(operator.mul, s, 1)\n\nfrom sympy.core.cache import cacheit\n\n@cacheit\ndef Pbt(t, p, ht):\n    s = mpf(0)\n    for i in range(0, ceil(ht-1)):\n        s += binomial(t,i)*(1-p)**(t-i)*p**i\n    return s\n\ndef Pst(t, p, t0=3):\n    num = [Pbt(a, p, mpf(1)/3.0 * a) for a in arange(t-1, t0, -2)]\n    denom = [Pbt(a, p, mpf(1)/3.0 * a) for a in arange(t, t0, -2)]\n    return (1-Pbt(t+1,p, mpf(1)/3.0 * t+1)) * (mult(num) / mult(denom))",
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 32
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": "pyplot.plot([Pst(x, 0.25, 4) for x in arange(3,300,1)])",
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "pyout",
       "prompt_number": 34,
       "text": "[<matplotlib.lines.Line2D at 0x54cfcd0>]"
      },
      {
       "output_type": "display_data",
       "png": "iVBORw0KGgoAAAANSUhEUgAAAXgAAAD9CAYAAAC2l2x5AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XtYVXW6B/DvBjYQFwEvXALME5iB4obSKMvErkcttLGe\n0dPFTGd4fLqfpqbOZdQ55WQ1pwdzpunqaNNkPc5pZE5oTpOY6UGmxLJsRixRQCUVUEEQ2Lznj3cW\n+8Lmomz23iy/n+fZz2bvvVhr/Vj6rne/v9/6LYuICIiIyHSC/L0DREQ0MBjgiYhMigGeiMikGOCJ\niEyKAZ6IyKQY4ImITKrHAF9VVYWpU6di7NixGDduHFasWOFxuYceegijR4+GzWZDeXn5gOwoERGd\nnZCePrRarXjxxReRnZ2NxsZGXH755bjxxhuRkZHRuUxxcTH27duHiooK7NixA4sWLUJpaemA7zgR\nEfWsxww+MTER2dnZAICoqChkZGTg0KFDLssUFRVh3rx5AIDc3Fw0NDSgtrZ2gHaXiIj6qs81+MrK\nSpSXlyM3N9fl/ZqaGqSmpna+TklJQXV1tff2kIiIzkmPJRpDY2Mjbr/9dhQWFiIqKqrL5+6zHVgs\nli7LeHqPiIh6d64zyvSawbe1tWH27Nm46667MGvWrC6fJycno6qqqvN1dXU1kpOTu91Jsz4WL17s\n931g29g+ts98j/7oMcCLCBYsWIDMzEw88sgjHpfJz8/HmjVrAAClpaWIjY1FQkJCv3aKiIj6r8cS\nzbZt2/C73/0O48ePR05ODgBg2bJlOHjwIACgoKAA06dPR3FxMdLT0xEZGYlVq1Z1u77164EjR4CC\nAi+2gIiIPOoxwF9zzTXo6OjodSUrV67s08b27QOcqjmmkpeX5+9dGDBmbhvA9g12Zm9ff1ikv0We\nvm7IYsFzzwmqqoBurpciIiI3FovlnGvxPp2qQASw2325RSKi85dPA3xHhz6IiGjgMcATEZmUz0s0\nDPBERL7BDJ6IyKSYwRMRmRQzeCIik2KAJyIyKZZoiIhMihk8EZFJMYMnIjIpn2fwnKqAiMg3WKIh\nIjIplmiIiEyKGTwRkUkxgyciMilm8EREJsUAT0RkUizREBGZFDN4IiKTYgZPRGRSzOCJiEyKAZ6I\nyKR8XqLhXDRERL7BDJ6IyKQY4ImITIqjaIiITIoZPBGRSTGDJyIyKWbwREQmxQBPRGRSLNEQEZkU\nM3giIpNiBk9EZFI+z+A5VQERkW+wRENEZFIs0RARmRQzeCIik2IGT0RkUszgiYhMigGeiMikWKIh\nIjIpZvBERCbFDJ6IyKSYwRMRmRQDPBGRSfUa4O+77z4kJCQgKyvL4+clJSWIiYlBTk4OcnJy8PTT\nT3e7LhHORUNE5CshvS0wf/58PPjgg7jnnnu6XWbKlCkoKirqdWPM4ImIfKfXDH7y5MmIi4vrcRkR\n6dPG2MlKROQ7vWbwvbFYLNi+fTtsNhuSk5PxwgsvIDMz0+OyBw4sQXs7sGQJkJeXh7y8vP5unojI\nVEpKSlBSUuKVdVmkD+l3ZWUlbr31VuzevbvLZ6dOnUJwcDAiIiKwYcMGPPzww9i7d2/XDVksuPJK\nQVkZ6/BERH1lsVj6XCVx1+9RNNHR0YiIiAAATJs2DW1tbairq/O4LEs0RES+0+8AX1tb23l2KSsr\ng4hg6NChHpc1gvs5noyIiOgs9FqDnzt3LrZs2YJjx44hNTUVS5cuRVtbGwCgoKAA69atw8svv4yQ\nkBBERERg7dq13a7LCPAdHUBwsHcaQEREnvWpBu+VDVksuOwywc6dQGsrYLX6YqtERIObX2vwZ8M5\ngyciooHl88nGAAZ4IiJf8EsGz2GSREQDjyUaIiKTYomGiMikmMETEZkUM3giIpNiBk9EZFIM8ERE\nJsUSDRGRSTGDJyIyKWbwREQm5fMM3mJhgCci8gWfB/iQEAZ4IiJf8HmJJjiYc9EQEfkCM3giIpPy\neQbPAE9E5Bs+z+CDgxngiYh8gSUaIiKTYomGiMikWKIhIjIpZvBERCbFGjwRkUkxwBMRmRRLNERE\nJsVOViIik/JLiYZz0RARDTy/TDbGDJ6IaOCxk5WIyKTYyUpEZFLsZCUiMimfBniAAZ6IyFcY4ImI\nTMqnAT4oSB8M8EREA48BnojIpHwa4C0WlmiIiHyFGTwRkUn5PIMPCuJUBUREvsAMnojIpBjgiYhM\nyi8lGgZ4IqKBxwyeiMikmMETEZkUM3giIpNigCciMimWaIiITIoZPBGRSfUa4O+77z4kJCQgKyur\n22UeeughjB49GjabDeXl5d0uxwyeiMh3eg3w8+fPx8aNG7v9vLi4GPv27UNFRQVeffVVLFq0qPuN\nBXGyMSIiX+k1wE+ePBlxcXHdfl5UVIR58+YBAHJzc9HQ0IDa2lrPGwviXDRERL4S0t8V1NTUIDU1\ntfN1SkoKqqurkZCQ0GXZkyeXoKwM2LsXGD8+D3l5ef3dPBGRqZSUlKCkpMQr6+p3gAcAEXF5bbFY\nPC4XF7cEkyYBKSkAYzsRUVd5ea7J79KlS895Xf0eRZOcnIyqqqrO19XV1UhOTva8MQ+jaAoKWJMn\nIhoI/Q7w+fn5WLNmDQCgtLQUsbGxHsszQNdRNB0dwKuvAq2t/d0LIiJy12uJZu7cudiyZQuOHTuG\n1NRULF26FG1tbQCAgoICTJ8+HcXFxUhPT0dkZCRWrVrV7brcM/gzZ/SZna5ERN7Xa4B/5513el3J\nypUr+7Qx9wy+pUWf29v79OtERHQW/HolqxHgmcETEXlfQAR4ZvBERN7n18nGWIMnIho4zOCJiEzK\n5xl8cLAjY2eAJyIaOH7N4FmiISIaOCzREBGZlF87WTlMkoho4AREiYYZPBGR9wVEBs8AT0TkfQFR\ng2eJhojI+1iiISIyqYAo0djtwI9/DHz1lS/3hojI3AKiRNPeDuzZA9TU+HJviIjMLSAy+PZ2oK3N\nUbIhIqL+C4gavN2uQZ53diIi8h6/BHhPc9Ewgyci8i6/TDbmqZOVGTwRkXcFRImGGTwRkff5vZPV\nYtEA397OAE9E5E1+HyYZEcESDRHRQPB7iSYqiiUaIqKB4PcSTWQkM3giooHg8wzeanUEciPAM4Mn\nIvI+nwf46Gjg1Cl97VyiYScrEZF3+bxEExXlCPDOJZq2NpZoiIi8ya8ZfEsLM3giooHi8wzePcBH\nRmrmLsIMnojIm/xeg4+MdExZwAyeiMh7/BLgGxs1YzdKNAzwRETe5/MSTUgIEBoKnD6tJZmICEeA\nZ4mGiMh7fJ7BA5rFHz+uY+KtVs8Z/BdfAI8/7su9IyIyF59n8IAG+KNHgfBwnT7YUwb/3XfAzp2+\n3DsiInPxWwZvBPiQEKC5Wd93zuCbm3VsPBERnRu/Bfj9+4FhwzTAe8rgm5tZkyci6g+/lWj27gWS\nklxLNMzgiYi8x28ZfEWFBngjgw8P7xrgmcETEZ07v2bwiYmODN64otXAAE9E1D9+yeCjonSUjHMG\nHxnJEg0RkTf5rURjt7vW4JnBExF5l99KNAAzeCKigeS3DB7QGrxzgDdmlQSYwRMR9ZffM3ijRBMa\nqsHeyNpbWpjBExH1h98y+PBwICbGkcGHhABhYY4yjZHBt7QA117ry70kIjIHvwX4xETN6I0M3mrV\nAG+UZYwa/KlTwKefOko3RETUN34p0aSnA7Nm6c8hIY7n0FDXDB5wnTueiIj6zi8Z/MiRwIsv6s9G\ngDcyePcAb9z9qanJd/tJRGQGvQb4jRs34tJLL8Xo0aOxfPnyLp+XlJQgJiYGOTk5yMnJwdNPP93t\nuowM3llwsD4bNXjnEg3AAE9EdK5CevrQbrfjgQcewEcffYTk5GRMnDgR+fn5yMjIcFluypQpKCoq\n6nVjQR5OJ84ZfHclGoABnojobPWYwZeVlSE9PR2jRo2C1WrFnDlzsH79+i7LSR97QD0F+J4y+OBg\nRwZ/+rTjd9auBbZv79MmiYjOWz1m8DU1NUhNTe18nZKSgh07drgsY7FYsH37dthsNiQnJ+OFF15A\nZmamx/Vt27YES5boz3l5ecjLy+sxgx8yBDh5Ul87Z/AffghkZwOTJvW9oUREg0FJSQlKSkq8sq4e\nA7zFU9HczWWXXYaqqipERERgw4YNmDVrFvbu3etx2WuvdQR4Q08ZfHKy5xJNUxNH1RCRORnJr2Hp\n0qXnvK4eSzTJycmoqqrqfF1VVYWUlBSXZaKjoxEREQEAmDZtGtra2lBXV+d5Y73U4I1RNO3tQEeH\nTl/gqZO1sZEBnoioNz0G+AkTJqCiogKVlZVobW3Fu+++i/z8fJdlamtrO2vwZWVlEBEMHTrU4/o8\nfSHwNA6+uRm44AJ97SnANzW5TkxGRERd9ViiCQkJwcqVK3HzzTfDbrdjwYIFyMjIwCuvvAIAKCgo\nwLp16/Dyyy8jJCQEERERWLt2bbfr66mT1flK1r4EeGbwREQ96zHAA1p2mTZtmst7BQUFnT/ff//9\nuP/++/u0sd4y+LAwDe5GgLdaGeCJiM6VX65kdeacwY8dC3z2GUs0RETe4PcA75zBT58ObNjADJ6I\nyBv8MtmYM+dRNOPH6wVNX36p0wkbGXxsrOcA39gIHDjgm30nIhps/J7BO4+Dt1iAadOA//kf1xLN\n0KGOAG+3a3nmzBng/feBJ57w3f4TEQ0mAZXBA1qm+fBDR4mmsdE1wBvPRgZfX++6Pt7mj4hIBVQG\nDwA33KBZunMGP2yY5wDf1AScOOG6vnHjgO+/H5j9JyIaTPwe4N0z+CFDgMmTXTtZnTN4Y+qCM2e0\nXt/Q4Lq+Q4e6ZvVEROejgCnRhDiNyJ85U4N6aKhm8+4lmqAgzeBPn3bN4O12/dx55kkiovNVrxc6\neZNRjvH0npHBA8ADD+hcND/5ib52L9HExXkO8EZ2b8wlT0R0PvNpBn/PPV3f85TBBwdrwDeCvnsG\nP3Soo0TT0uK46MmYWpgZPBGRjwO8pznI3DtZnYWGOn7POcAPG+boZAUcWbynm4MQEZ2vfBrgPbFY\ntKbuXKIxGO+5l2iMAG8EcqOjlRk8EZGD3wM8oFl8Txl8dLTW5Ddt0iDuXKIBPGfw27cDR44M/L4T\nEQWqgAjwISGeM3gjwIeHA/HxQH6+XgQ1dKhe0NTYqJ95yuCXL9d5bQwHDwI/+tHAtoOIKJAETID3\nlMEbQT8sTOecufdeoLQUiIrS4F9fDyQlec7gT50Cjh1zrGvfPuCTTwa0GUREASUgArwxasadkcGH\nhennxlWqkZGaudfVARde6DmDdw/wDQ2OYZREROeDgAjwfcngASArS5/dA7xzBh8SogH+5Eng6FHH\nuhoaHBk+EdH5ICACfG8ZfHi4Po8bp8+RkRr07XYt0Thn8PHxeqGTewZfX68Z/D9uH0tEZHoBEeC7\ny+CdSzSADo9MSnJk8EDXDD4x0ZHBu5doRHiVKxGdPwImwPc0Dt4I8ABwxRVAQoIG+KAg/Xn3buBf\n/1WDekKCZupNTV0DPKAngZde0mGXhk2bdHQOEZGZBESAnzNHM2937hk8oDcDmTpV34uI0Ls9lZQA\nhYVac09M1I5Yi6VriQbQbP+hh4DaWsdnGzcCH33k9WYREflVQAT4Z5/VaYLdWa2a3TtPUmZMORwe\nrgE+MRFITgYuukhv2J2YqBc4JSRott7erssbGXx1tT4fOuRY57FjXeeVJyIa7AIiwHcnNNQ1e3dm\nBPgrrwT+/nfgsst0VE1CgmbnMTE662RdnS7f0KBZfVWVvj582LGu48cZ4InIfAI6wFutjs5Ud2Fh\n2tkK6HNmpv6cmKjlmCFDgOHDHUMljYuiDh7U184Bnhk8EZlRQAf4vmTwhrFj9TkhQZ+jozXAr1ql\nj4YGIDXVkcF7KtGcOAFs2eK6nbvvBiorvdIcIiKfMk2ANzJ4I8AbGfx//zewdq1m8M4B3lOJpqTE\ncZMRw+bNwLffeqU5REQ+5dM7Op2tqCjPna+AYxSN4ZJLtA4/fLi+jo7Wk0BKCrBzJ9DWpsH/m2+A\nESMcAb6tTYN7ZKSWcw4ccKxTREfkGHV8IqLBJKAz+Esu0Qzak/BwRw0e0ID/+eeO94YMAW67Dfjt\nb3Vq4dhYDfpVVZrtGyWaujq9wfeJExrgjx51TENcX68ngO4CvN3ulWYSEQ2IgA7wgI6E8cS9ROP8\nPqDBfNo04LrrdIqDuDj9RnDypNbrDx8GFi4E/vIXHWLZ0uKYP97oiDXGytfV6beA775zbOfwYUdZ\niIgoEAV8gO+Oe4nGEBSkGblzaScrSzP4qCh9femlGqDfeEPr88OH6/L79unnRpnGOcCvWAGsXu1Y\n53ff6YNz2xBRoBq0Ad69ROMsIkIzeMO4cY4SDaC1+KFDgTFjtGN1+HAdN79vnw6lNAK8kdHX1WlJ\np6LCsc6aGr2IyriAiogo0AR0J2tPbrtNa+ueRES4ZvCzZgEjRzqWj4vTuWsmTgRuvFEnMYuJAfbs\nAW6+WZ8ffhi4+GK9Svb4cc34jekOAA3wgHbC1tRo2cdi0ffeew/429+An/3M++0mIuqrQRvgjbnh\nPXHP4FNT9VFcrK/j4oCnntLb/lmtjgy+rQ2YMAF47jntaP3hD4GMDM3gDx/Wz0U0kBsB/uhR4Ac/\n0AnLsrP1vc8/B/buHZh2ExH11aAt0fTEPYM3GEHf6LgNDdXyjZHBAxrgT5/W7H3DBg3wR45o52xw\nsGbsTU0a4C0WvQjq6FHgq68c26msdL2QytmuXbyzFBH5hikD/KxZWl93Z3SyOo/MueUWLa/ExGhd\n/6qrdLbJggIN6hkZeqFTfLwO21ywQGe/rKkB0tM1WwdcA/yBA5rxf/898G//5roPixa53gyciGig\nmDLA/+xneoGTOyODN7J1APj5z4Hp0/W9ESO087WwEMjN1c8zMnTu+KQkDfAffKA37z54EMjJ0Rks\ng4N1TnpjRE1lpQb4sjK9kratzbG9igrHMEx3xkyXRETeYMoA352YGC3HOE8/7PzZiBGO15dfrkMu\nR43SzD4pScs5t9+uHa8HDmiALy8Hrr5aM/jbbtOhlA0NejLZsUM7dr/+WtdZX68dtlVVwLZtrhdx\n1dcDaWnddxwTEZ2tQdvJei5GjNALljxxD/BRUXpzkZEjNatPSgIefVSz9Acf1Az9oou0Hj9limby\nRUUa+FNT9aSwdaueJD77zHUmy4MHgbfe0hLO1Kn63tdfa6dvZWXX8lJdnZaJJk70+p+EiEzsvMrg\nAQ3Yngwb1vWuUjNnaoAeOlTv/Wq1asfstddqFh8fr8uNGqXB9xe/0E7Uiy7SgL5jh15Ju2aNLv/B\nB7rswYM6FPPTTx1lnT179Pnbb4Ef/1jnuDe88w7w+OPe/CsQ0fngvMrge3LnnVpi8cS42bdh5kwd\nWmkE+JEjdZhkaKjOfTNqlF4E1dKiHbILF+rv//rXwK23An/6k9b17Xb9lrB/v5ZtQkK0Rv/73+so\nnief1PWXl+u4ehGgtFQ7gg2LF+vrf/7ngfirENFgdt5l8N0JD9dM3ZPrrtOavCEqSi+IMko6I0c6\n7h87Y4aWWJKSdBjl7NnAY48Bv/mNDqe85hqd2KyjQ0fw3Hkn8B//Afzf/+k3g02btOzz4Yea1VdV\n6beC2lqdq37yZNfJz957D/j4Yw3+xu0JDS+95NrBS0TnF4uIb2ZTsVgs8NGmfKa9XU8KtbU6/w2g\nWbnFAqxcCTz/vGP++ZYWzfo3bgTuvVenS3jgAb0ZyQUXaP3+xReBJ57Quvz27dpRe/XVjtLOxIla\n7vntb3W8flKSfru46Sb99rFli5ZzAO0juPBCrf9fdpnjKltA71VbWqpX8RJRYOtP7GQG3w8hIVpP\nN4I7oCN0goJ0mObo0Y73w8P1StrcXO2EzcwE5s7VTH3uXP1WkJenGfc11wCTJulsmP/7v1rTv/xy\nYN06zfJ/8Qtg/Hjgpz/VkTdffqnr/sMfNLB/9pnOsQPoUM077tDPDO++C8yfr1m/e4a/ezdvcEJk\nFqzB91NsrOf3Z8zQLNvZtdfq88UX65BLQDPrWbP0xJCWpu/l5OhcORERmvmL6AyYp08DTz8N5OcD\n99+vJZgnngBeew346CNd/1VXaZC/4QY9iXz8sX47qKrSjL6jQ08GNTXAX/+q1wB89pl+QwB0uyNH\nan/BoUPAP/2TY/9ffFG/fXQ3hTMRBRaWaPygoUEzek83FM/O1jKOMaLnwAEN7N98o8G1vl6DdFCQ\nXnhVWAi88ILeV/Y3vwGefVaD8sqVwJtvagdvXp5ObVxfr8NBT5zQk0hrqw4bnT1bTyT33ac/Dx+u\nJ5AVK7TT98QJ/TwtDXj9dS0LjRyppR9AO4AfecRxha5zOaihQcf2G7dSJKKz05/YyQA/SHz/vV5E\n5Txc8uRJrdM/8ohm+StW6PtHj+qUCG+/rQH9uef0ZBAVpSeEAweAe+7R8fyvvqrj+ydN0hr+rbdq\nRn/qlJaSWls1yM+YoaWgtDQt48TEAP/5nzqvzhdf6Mll40bd7urV2tEcF6c3LT9xQkcG/eUveqGY\n4d//HbjrLl3WYnGcGER0ZNHFF/vu70sUqPoVO6UXGzZskDFjxkh6ero8++yzHpd58MEHJT09XcaP\nHy87d+70uEwfNjWobd682W/bPn5cpK7O82c/+YlIVZXj9dGjIl98IfLllyJDhog0N4scOybS0SGy\nYIHIxo0iP/qRSEaGSGWlyKOPijz11GYBRD7+WMRqFbn5ZpHZs0VGjxZJTRWJihJZtEhk6FCR7GyR\nhASRiAiR664TGTFCZNgwkVtuEQkKEnnqKZE77hB5/XV9PXeuyL/8i0h+vsjXX4ts2SKybp1+9tVX\nIosXi3zzjciZMyJ2u8jvfy/y2GMi7e0in3/u2tZNmxxt7ehw/ay+vvu/nz+PnS+wfYNbf2Jnj7/Z\n3t4uaWlpsn//fmltbRWbzSZ79uxxWeaDDz6QadOmiYhIaWmp5Obmen0nB4PFixf7exfOWm2t5/e/\n+UakpMTxevHixbJ7twbNRYtEPvlEA+6ZMyK7dokUFop8+61IWJhIRYXI229rsJ4xQ+T55/UkMXSo\nrnPKFJEnnhCxWEReflkkNlZPFNOm6c/JyXriueMOkQsvFLn4YpH4eP3siitEhg8XSUzUdQcFidx1\nl8jUqSL/9V8iF1wgMnmyyL336rKrVoksWyayZo1IcLDIH/4gcv/9Im+9pW3YvVvkzTdFLr10sRw9\nKvLaa3qiPHBATyCvv64njfp6kb/9Tf8WHR36KCoSOX1apKVFlzW0t3f/d/WXwfhv82yYvX39iZ09\ndrKWlZUhPT0do/7RAzdnzhysX78eGRkZncsUFRVh3rx5AIDc3Fw0NDSgtrYWCSy6BjzjQi13l16q\nD2dGp/Cvf+36vs2mD0DLSEOG6CybgI4AAnT+nblzdZinMbpn4UJdLjpaSzQ2m9bqGxt1mZkzdT7+\nwkLtX4iL0wvEhg3TO3ktXar9EqtX60VeK1fqNQHPP6/9DXPm6MihmBhg+XJdzw9/qMvu2KGjh44c\n0auTL7hARyrZbFq2CgvT92Ji9JqDsDCguVmHxB45ovtbWakd0wcOaOlr5EjtF2lp0QvT7rlHr2mY\nNEn7UNLTtSM8NBS4/nqdhfSqq/Rvk5am+37TTbqdjg5dX1ub7sfmzdqeXbt0Wbtdr8H49lv9e11z\njZbNMjO1bJeSonMdZWRo+/7+dy27nTmjbfrrX7VD/sABx9xMERFaljt6VPtn9u3TfTh9WrdVUaHL\nRkdrv0p8vK4vOlrvfTB6tP5+eLi20WLRYcT19dq5b5QTRXR7jY36HB6uy1mtjn9Tzc2u/VPOfTp2\nu+e5pMizHgN8TU0NUlNTO1+npKRgx44dvS5TXV3NAH8e8jQHP6CBYdgw1/eMIaR33ul4z2rVYDln\njr4uKtJnYyTPo486lv3BD/Q//jPPuK7nuut0+GpoqGP59nZ9b+pUDc5GgGhq0mC6fLleQHbTTdpf\nEBOjQTo9XYNXXZ0G0W++0cC2bZuOfHr5Zb1OQUSDXlOT/v7VV+tnb7+tQTkuTi9ae/JJDXi7d+u0\n0zt36oll5069g9if/qSd0ZGROjdRUJCeNK++Wk9ckybpNRHh4TpSasQIPfmsW6cjtt57T0d17d+v\nr9es0bbW1+sJNSxMA3h2tvbljB6t+22xaMANDdXjtHu3/p1qajTwf/+9HoOjRzWox8bqz+HhGtRH\njdKZUMPC9KRkDL0NCdG2nD6t22ht1WMREaF/d2OZ5mY9kZ05o8f/1Ck9GbS06O+FhTlODMYJo65O\nl7Xbdb0rVugysbGOk0lbm27r+HF9327X7YaH67bOnNHX8fHa/o4O3Z8LLtD9bG3VfwsnT+q/zbAw\nfT55UtdhnNBDQ/V94+8YE+MY/Wa16sNYn/F/xG7XR3u7PsLDddvXXqt9Zl7TU3q/bt06WbhwYefr\nt956Sx544AGXZW655Rb59NNPO19ff/318rl7cfQfXzP44IMPPvg4+8eAlGiSk5NRZVyKCaCqqgop\nbhOtuy9TXV2N5OTkLusSjqAhIvKpHq9knTBhAioqKlBZWYnW1la8++67yM/Pd1kmPz8fa9asAQCU\nlpYiNjaW5RkiogDQYwYfEhKClStX4uabb4bdbseCBQuQkZGBV155BQBQUFCA6dOno7i4GOnp6YiM\njMSqVat8suNERNSLcy7u9FFfxtEPNhdddJFkZWVJdna2TJw4UUREjh8/LjfccIOMHj1abrzxRqnv\naeB1gJk/f77Ex8fLuHHjOt/rqT3Lli2T9PR0GTNmjHz44Yf+2OWz4ql9ixcvluTkZMnOzpbs7Gwp\nLi7u/Gwwte/gwYOSl5cnmZmZMnbsWCksLBQR8xy/7tpnluPX3NwsV1xxhdhsNsnIyJAnn3xSRLx3\n/AY0wPdlHP1gNGrUKDl+/LjLe48//rgsX75cRESeffZZ+elPf+qPXTsnn3zyiezcudMlAHbXnq+/\n/lpsNpu0trbK/v37JS0tTex2u1/2u688tW/JkiXyy1/+ssuyg619hw8flvLychEROXXqlFxyySWy\nZ88e0xzc+6RGAAADWUlEQVS/7tpnluMnItLU1CQiIm1tbZKbmytbt2712vEb0NkkncfRW63WznH0\nZiBuncbO1wPMmzcPf/zjH/2xW+dk8uTJiHObQay79qxfvx5z586F1WrFqFGjkJ6ejrKyMp/v89nw\n1D7Ac8f/YGtfYmIisrOzAQBRUVHIyMhATU2NaY5fd+0DzHH8ACAiIgIA0NraCrvdjri4OK8dvwEN\n8J7GyBsHZzCzWCy44YYbMGHCBLz22msA4HJxV0JCAmpra/25i/3WXXsOHTrkMpJqMB/Tl156CTab\nDQsWLEBDQwOAwd2+yspKlJeXIzc315THz2jflVdeCcA8x6+jowPZ2dlISEjA1KlTMXbsWK8dvwEN\n8BbnS9BMZNu2bSgvL8eGDRvwq1/9Clu3bnX53GKxmKrtvbVnMLZ10aJF2L9/P3bt2oWkpCQ89thj\n3S47GNrX2NiI2bNno7CwENHR0S6fmeH4NTY24vbbb0dhYSGioqJMdfyCgoKwa9cuVFdX45NPPsHm\nzZtdPu/P8RvQAN+XcfSDUdI/btA6YsQI3HbbbSgrK0NCQgKOHDkCADh8+DDiu5sHYJDorj19ve4h\n0MXHx3f+x1m4cGHn19zB2L62tjbMnj0bd999N2bNmgXAXMfPaN9dd93V2T4zHT9DTEwMZsyYgc8/\n/9xrx29AA3xfxtEPNqdPn8apU6cAAE1NTdi0aROysrKQn5+P1atXAwBWr17d+Q9xsOquPfn5+Vi7\ndi1aW1uxf/9+VFRU4IorrvDnrp6Tw4cPd/78/vvvIysrC8Dga5+IYMGCBcjMzMQjjzzS+b5Zjl93\n7TPL8Tt27Fhneam5uRl//vOfkZOT473jN6DdwyJSXFwsl1xyiaSlpcmyZcsGenMD7rvvvhObzSY2\nm03Gjh3b2abjx4/L9ddfPyiHSc6ZM0eSkpLEarVKSkqKvPnmmz2255lnnpG0tDQZM2aMbNy40Y97\n3jfu7XvjjTfk7rvvlqysLBk/frzMnDlTjhw50rn8YGrf1q1bxWKxiM1m6xwyuGHDBtMcP0/tKy4u\nNs3x+/LLLyUnJ0dsNptkZWXJc889JyI9x5OzaZ/PbvhBRES+xZtuExGZFAM8EZFJMcATEZkUAzwR\nkUkxwBMRmRQDPBGRSf0/A0/kLjsF/p8AAAAASUVORK5CYII=\n"
      }
     ],
     "prompt_number": 34
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": "",
     "language": "python",
     "metadata": {},
     "outputs": []
    }
   ],
   "metadata": {}
  }
 ]
}
	{
	"metadata": {
	"name": "Untitled0"
	},
	"nbformat": 3,
	"nbformat_minor": 0,
	"worksheets": [
	{
	"cells": [
	{
	"cell_type": "code",
	"collapsed": false,
	"input": "from sympy.mpmath import *\nlibmp.BACKEND",
	"language": "python",
	"metadata": {},
	"outputs": [
	{
	"output_type": "pyout",
	"prompt_number": 1,
	"text": "'gmpy'"
	}
	],
	"prompt_number": 1
	},
	{
	"cell_type": "code",
	"collapsed": false,
	"input": "mp.dps = 50\n\nimport operator\ndef mult(s):\n return reduce(operator.mul, s, 1)\n\nfrom sympy.core.cache import cacheit\n\n@cacheit\ndef Pbt(t, p, ht):\n s = mpf(0)\n for i in range(0, ceil(ht-1)):\n s += binomial(t,i)(1-p)(t-i)p*i\n return s\n\ndef Pst(t, p, t0=3):\n num = [Pbt(a, p, mpf(1)/3.0 a) for a in arange(t-1, t0, -2)]\n denom = [Pbt(a, p, mpf(1)/3.0 * a) for a in arange(t, t0, -2)]\n return (1-Pbt(t+1,p, mpf(1)/3.0 * t+1)) * (mult(num) / mult(denom))",
	"language": "python",
	"metadata": {},
	"outputs": [],
	"prompt_number": 32
	},
	{
	"cell_type": "code",
	"collapsed": false,
	"input": "pyplot.plot([Pst(x, 0.25, 4) for x in arange(3,300,1)])",
	"language": "python",
	"metadata": {},
	"outputs": [
	{
	"output_type": "pyout",
	"prompt_number": 34,
	"text": "[<matplotlib.lines.Line2D at 0x54cfcd0>]"
	},
	{
	"output_type": "display_data",
	"png": "iVBORw0KGgoAAAANSUhEUgAAAXgAAAD9CAYAAAC2l2x5AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XtYVXW6B/DvBjYQFwEvXALME5iB4obSKMvErkcttLGe\n0dPFTGd4fLqfpqbOZdQ55WQ1pwdzpunqaNNkPc5pZE5oTpOY6UGmxLJsRixRQCUVUEEQ2Lznj3cW\n+8Lmomz23iy/n+fZz2bvvVhr/Vj6rne/v9/6LYuICIiIyHSC/L0DREQ0MBjgiYhMigGeiMikGOCJ\niEyKAZ6IyKQY4ImITKrHAF9VVYWpU6di7NixGDduHFasWOFxuYceegijR4+GzWZDeXn5gOwoERGd\nnZCePrRarXjxxReRnZ2NxsZGXH755bjxxhuRkZHRuUxxcTH27duHiooK7NixA4sWLUJpaemA7zgR\nEfWsxww+MTER2dnZAICoqChkZGTg0KFDLssUFRVh3rx5AIDc3Fw0NDSgtrZ2gHaXiIj6qs81+MrK\nSpSXlyM3N9fl/ZqaGqSmpna+TklJQXV1tff2kIiIzkmPJRpDY2Mjbr/9dhQWFiIqKqrL5+6zHVgs\nli7LeHqPiIh6d64zyvSawbe1tWH27Nm46667MGvWrC6fJycno6qqqvN1dXU1kpOTu91Jsz4WL17s\n931g29g+ts98j/7oMcCLCBYsWIDMzEw88sgjHpfJz8/HmjVrAAClpaWIjY1FQkJCv3aKiIj6r8cS\nzbZt2/C73/0O48ePR05ODgBg2bJlOHjwIACgoKAA06dPR3FxMdLT0xEZGYlVq1Z1u77164EjR4CC\nAi+2gIiIPOoxwF9zzTXo6OjodSUrV67s08b27QOcqjmmkpeX5+9dGDBmbhvA9g12Zm9ff1ikv0We\nvm7IYsFzzwmqqoBurpciIiI3FovlnGvxPp2qQASw2325RSKi85dPA3xHhz6IiGjgMcATEZmUz0s0\nDPBERL7BDJ6IyKSYwRMRmRQzeCIik2KAJyIyKZZoiIhMihk8EZFJMYMnIjIpn2fwnKqAiMg3WKIh\nIjIplmiIiEyKGTwRkUkxgyciMilm8EREJsUAT0RkUizREBGZFDN4IiKTYgZPRGRSzOCJiEyKAZ6I\nyKR8XqLhXDRERL7BDJ6IyKQY4ImITIqjaIiITIoZPBGRSTGDJyIyKWbwREQmxQBPRGRSLNEQEZkU\nM3giIpNiBk9EZFI+z+A5VQERkW+wRENEZFIs0RARmRQzeCIik2IGT0RkUszgiYhMigGeiMikWKIh\nIjIpZvBERCbFDJ6IyKSYwRMRmRQDPBGRSfUa4O+77z4kJCQgKyvL4+clJSWIiYlBTk4OcnJy8PTT\nT3e7LhHORUNE5CshvS0wf/58PPjgg7jnnnu6XWbKlCkoKirqdWPM4ImIfKfXDH7y5MmIi4vrcRkR\n6dPG2MlKROQ7vWbwvbFYLNi+fTtsNhuSk5PxwgsvIDMz0+OyBw4sQXs7sGQJkJeXh7y8vP5unojI\nVEpKSlBSUuKVdVmkD+l3ZWUlbr31VuzevbvLZ6dOnUJwcDAiIiKwYcMGPPzww9i7d2/XDVksuPJK\nQVkZ6/BERH1lsVj6XCVx1+9RNNHR0YiIiAAATJs2DW1tbairq/O4LEs0RES+0+8AX1tb23l2KSsr\ng4hg6NChHpc1gvs5noyIiOgs9FqDnzt3LrZs2YJjx44hNTUVS5cuRVtbGwCgoKAA69atw8svv4yQ\nkBBERERg7dq13a7LCPAdHUBwsHcaQEREnvWpBu+VDVksuOwywc6dQGsrYLX6YqtERIObX2vwZ8M5\ngyciooHl88nGAAZ4IiJf8EsGz2GSREQDjyUaIiKTYomGiMikmMETEZkUM3giIpNiBk9EZFIM8ERE\nJsUSDRGRSTGDJyIyKWbwREQm5fMM3mJhgCci8gWfB/iQEAZ4IiJf8HmJJjiYc9EQEfkCM3giIpPy\neQbPAE9E5Bs+z+CDgxngiYh8gSUaIiKTYomGiMikWKIhIjIpZvBERCbFGjwRkUkxwBMRmRRLNERE\nJsVOViIik/JLiYZz0RARDTy/TDbGDJ6IaOCxk5WIyKTYyUpEZFLsZCUiMimfBniAAZ6IyFcY4ImI\nTMqnAT4oSB8M8EREA48BnojIpHwa4C0WlmiIiHyFGTwRkUn5PIMPCuJUBUREvsAMnojIpBjgiYhM\nyi8lGgZ4IqKBxwyeiMikmMETEZkUM3giIpNigCciMimWaIiITIoZPBGRSfUa4O+77z4kJCQgKyur\n22UeeughjB49GjabDeXl5d0uxwyeiMh3eg3w8+fPx8aNG7v9vLi4GPv27UNFRQVeffVVLFq0qPuN\nBXGyMSIiX+k1wE+ePBlxcXHdfl5UVIR58+YBAHJzc9HQ0IDa2lrPGwviXDRERL4S0t8V1NTUIDU1\ntfN1SkoKqqurkZCQ0GXZkyeXoKwM2LsXGD8+D3l5ef3dPBGRqZSUlKCkpMQr6+p3gAcAEXF5bbFY\nPC4XF7cEkyYBKSkAYzsRUVd5ea7J79KlS895Xf0eRZOcnIyqqqrO19XV1UhOTva8MQ+jaAoKWJMn\nIhoI/Q7w+fn5WLNmDQCgtLQUsbGxHsszQNdRNB0dwKuvAq2t/d0LIiJy12uJZu7cudiyZQuOHTuG\n1NRULF26FG1tbQCAgoICTJ8+HcXFxUhPT0dkZCRWrVrV7brcM/gzZ/SZna5ERN7Xa4B/5513el3J\nypUr+7Qx9wy+pUWf29v79OtERHQW/HolqxHgmcETEXlfQAR4ZvBERN7n18nGWIMnIho4zOCJiEzK\n5xl8cLAjY2eAJyIaOH7N4FmiISIaOCzREBGZlF87WTlMkoho4AREiYYZPBGR9wVEBs8AT0TkfQFR\ng2eJhojI+1iiISIyqYAo0djtwI9/DHz1lS/3hojI3AKiRNPeDuzZA9TU+HJviIjMLSAy+PZ2oK3N\nUbIhIqL+C4gavN2uQZ53diIi8h6/BHhPc9Ewgyci8i6/TDbmqZOVGTwRkXcFRImGGTwRkff5vZPV\nYtEA397OAE9E5E1+HyYZEcESDRHRQPB7iSYqiiUaIqKB4PcSTWQkM3giooHg8wzeanUEciPAM4Mn\nIvI+nwf46Gjg1Cl97VyiYScrEZF3+bxEExXlCPDOJZq2NpZoiIi8ya8ZfEsLM3giooHi8wzePcBH\nRmrmLsIMnojIm/xeg4+MdExZwAyeiMh7/BLgGxs1YzdKNAzwRETe5/MSTUgIEBoKnD6tJZmICEeA\nZ4mGiMh7fJ7BA5rFHz+uY+KtVs8Z/BdfAI8/7su9IyIyF59n8IAG+KNHgfBwnT7YUwb/3XfAzp2+\n3DsiInPxWwZvBPiQEKC5Wd93zuCbm3VsPBERnRu/Bfj9+4FhwzTAe8rgm5tZkyci6g+/lWj27gWS\nklxLNMzgiYi8x28ZfEWFBngjgw8P7xrgmcETEZ07v2bwiYmODN64otXAAE9E1D9+yeCjonSUjHMG\nHxnJEg0RkTf5rURjt7vW4JnBExF5l99KNAAzeCKigeS3DB7QGrxzgDdmlQSYwRMR9ZffM3ijRBMa\nqsHeyNpbWpjBExH1h98y+PBwICbGkcGHhABhYY4yjZHBt7QA117ry70kIjIHvwX4xETN6I0M3mrV\nAG+UZYwa/KlTwKefOko3RETUN34p0aSnA7Nm6c8hIY7n0FDXDB5wnTueiIj6zi8Z/MiRwIsv6s9G\ngDcyePcAb9z9qanJd/tJRGQGvQb4jRs34tJLL8Xo0aOxfPnyLp+XlJQgJiYGOTk5yMnJwdNPP93t\nuowM3llwsD4bNXjnEg3AAE9EdK5CevrQbrfjgQcewEcffYTk5GRMnDgR+fn5yMjIcFluypQpKCoq\n6nVjQR5OJ84ZfHclGoABnojobPWYwZeVlSE9PR2jRo2C1WrFnDlzsH79+i7LSR97QD0F+J4y+OBg\nRwZ/+rTjd9auBbZv79MmiYjOWz1m8DU1NUhNTe18nZKSgh07drgsY7FYsH37dthsNiQnJ+OFF15A\nZmamx/Vt27YES5boz3l5ecjLy+sxgx8yBDh5Ul87Z/AffghkZwOTJvW9oUREg0FJSQlKSkq8sq4e\nA7zFU9HczWWXXYaqqipERERgw4YNmDVrFvbu3etx2WuvdQR4Q08ZfHKy5xJNUxNH1RCRORnJr2Hp\n0qXnvK4eSzTJycmoqqrqfF1VVYWUlBSXZaKjoxEREQEAmDZtGtra2lBXV+d5Y73U4I1RNO3tQEeH\nTl/gqZO1sZEBnoioNz0G+AkTJqCiogKVlZVobW3Fu+++i/z8fJdlamtrO2vwZWVlEBEMHTrU4/o8\nfSHwNA6+uRm44AJ97SnANzW5TkxGRERd9ViiCQkJwcqVK3HzzTfDbrdjwYIFyMjIwCuvvAIAKCgo\nwLp16/Dyyy8jJCQEERERWLt2bbfr66mT1flK1r4EeGbwREQ96zHAA1p2mTZtmst7BQUFnT/ff//9\nuP/++/u0sd4y+LAwDe5GgLdaGeCJiM6VX65kdeacwY8dC3z2GUs0RETe4PcA75zBT58ObNjADJ6I\nyBv8MtmYM+dRNOPH6wVNX36p0wkbGXxsrOcA39gIHDjgm30nIhps/J7BO4+Dt1iAadOA//kf1xLN\n0KGOAG+3a3nmzBng/feBJ57w3f4TEQ0mAZXBA1qm+fBDR4mmsdE1wBvPRgZfX++6Pt7mj4hIBVQG\nDwA33KBZunMGP2yY5wDf1AScOOG6vnHjgO+/H5j9JyIaTPwe4N0z+CFDgMmTXTtZnTN4Y+qCM2e0\nXt/Q4Lq+Q4e6ZvVEROejgCnRhDiNyJ85U4N6aKhm8+4lmqAgzeBPn3bN4O12/dx55kkiovNVrxc6\neZNRjvH0npHBA8ADD+hcND/5ib52L9HExXkO8EZ2b8wlT0R0PvNpBn/PPV3f85TBBwdrwDeCvnsG\nP3Soo0TT0uK46MmYWpgZPBGRjwO8pznI3DtZnYWGOn7POcAPG+boZAUcWbynm4MQEZ2vfBrgPbFY\ntKbuXKIxGO+5l2iMAG8EcqOjlRk8EZGD3wM8oFl8Txl8dLTW5Ddt0iDuXKIBPGfw27cDR44M/L4T\nEQWqgAjwISGeM3gjwIeHA/HxQH6+XgQ1dKhe0NTYqJ95yuCXL9d5bQwHDwI/+tHAtoOIKJAETID3\nlMEbQT8sTOecufdeoLQUiIrS4F9fDyQlec7gT50Cjh1zrGvfPuCTTwa0GUREASUgArwxasadkcGH\nhennxlWqkZGaudfVARde6DmDdw/wDQ2OYZREROeDgAjwfcngASArS5/dA7xzBh8SogH+5Eng6FHH\nuhoaHBk+EdH5ICACfG8ZfHi4Po8bp8+RkRr07XYt0Thn8PHxeqGTewZfX68Z/D9uH0tEZHoBEeC7\ny+CdSzSADo9MSnJk8EDXDD4x0ZHBu5doRHiVKxGdPwImwPc0Dt4I8ABwxRVAQoIG+KAg/Xn3buBf\n/1WDekKCZupNTV0DPKAngZde0mGXhk2bdHQOEZGZBESAnzNHM2937hk8oDcDmTpV34uI0Ls9lZQA\nhYVac09M1I5Yi6VriQbQbP+hh4DaWsdnGzcCH33k9WYREflVQAT4Z5/VaYLdWa2a3TtPUmZMORwe\nrgE+MRFITgYuukhv2J2YqBc4JSRott7erssbGXx1tT4fOuRY57FjXeeVJyIa7AIiwHcnNNQ1e3dm\nBPgrrwT+/nfgsst0VE1CgmbnMTE662RdnS7f0KBZfVWVvj582LGu48cZ4InIfAI6wFutjs5Ud2Fh\n2tkK6HNmpv6cmKjlmCFDgOHDHUMljYuiDh7U184Bnhk8EZlRQAf4vmTwhrFj9TkhQZ+jozXAr1ql\nj4YGIDXVkcF7KtGcOAFs2eK6nbvvBiorvdIcIiKfMk2ANzJ4I8AbGfx//zewdq1m8M4B3lOJpqTE\ncZMRw+bNwLffeqU5REQ+5dM7Op2tqCjPna+AYxSN4ZJLtA4/fLi+jo7Wk0BKCrBzJ9DWpsH/m2+A\nESMcAb6tTYN7ZKSWcw4ccKxTREfkGHV8IqLBJKAz+Esu0Qzak/BwRw0e0ID/+eeO94YMAW67Dfjt\nb3Vq4dhYDfpVVZrtGyWaujq9wfeJExrgjx51TENcX68ngO4CvN3ulWYSEQ2IgA7wgI6E8cS9ROP8\nPqDBfNo04LrrdIqDuDj9RnDypNbrDx8GFi4E/vIXHWLZ0uKYP97oiDXGytfV6beA775zbOfwYUdZ\niIgoEAV8gO+Oe4nGEBSkGblzaScrSzP4qCh9femlGqDfeEPr88OH6/L79unnRpnGOcCvWAGsXu1Y\n53ff6YNz2xBRoBq0Ad69ROMsIkIzeMO4cY4SDaC1+KFDgTFjtGN1+HAdN79vnw6lNAK8kdHX1WlJ\np6LCsc6aGr2IyriAiogo0AR0J2tPbrtNa+ueRES4ZvCzZgEjRzqWj4vTuWsmTgRuvFEnMYuJAfbs\nAW6+WZ8ffhi4+GK9Svb4cc34jekOAA3wgHbC1tRo2cdi0ffeew/429+An/3M++0mIuqrQRvgjbnh\nPXHP4FNT9VFcrK/j4oCnntLb/lmtjgy+rQ2YMAF47jntaP3hD4GMDM3gDx/Wz0U0kBsB/uhR4Ac/\n0AnLsrP1vc8/B/buHZh2ExH11aAt0fTEPYM3GEHf6LgNDdXyjZHBAxrgT5/W7H3DBg3wR45o52xw\nsGbsTU0a4C0WvQjq6FHgq68c26msdL2QytmuXbyzFBH5hikD/KxZWl93Z3SyOo/MueUWLa/ExGhd\n/6qrdLbJggIN6hkZeqFTfLwO21ywQGe/rKkB0tM1WwdcA/yBA5rxf/898G//5roPixa53gyciGig\nmDLA/+xneoGTOyODN7J1APj5z4Hp0/W9ESO087WwEMjN1c8zMnTu+KQkDfAffKA37z54EMjJ0Rks\ng4N1TnpjRE1lpQb4sjK9kratzbG9igrHMEx3xkyXRETeYMoA352YGC3HOE8/7PzZiBGO15dfrkMu\nR43SzD4pScs5t9+uHa8HDmiALy8Hrr5aM/jbbtOhlA0NejLZsUM7dr/+WtdZX68dtlVVwLZtrhdx\n1dcDaWnddxwTEZ2tQdvJei5GjNALljxxD/BRUXpzkZEjNatPSgIefVSz9Acf1Az9oou0Hj9limby\nRUUa+FNT9aSwdaueJD77zHUmy4MHgbfe0hLO1Kn63tdfa6dvZWXX8lJdnZaJJk70+p+EiEzsvMrg\nAQ3Yngwb1vWuUjNnaoAeOlTv/Wq1asfstddqFh8fr8uNGqXB9xe/0E7Uiy7SgL5jh15Ju2aNLv/B\nB7rswYM6FPPTTx1lnT179Pnbb4Ef/1jnuDe88w7w+OPe/CsQ0fngvMrge3LnnVpi8cS42bdh5kwd\nWmkE+JEjdZhkaKjOfTNqlF4E1dKiHbILF+rv//rXwK23An/6k9b17Xb9lrB/v5ZtQkK0Rv/73+so\nnief1PWXl+u4ehGgtFQ7gg2LF+vrf/7ngfirENFgdt5l8N0JD9dM3ZPrrtOavCEqSi+IMko6I0c6\n7h87Y4aWWJKSdBjl7NnAY48Bv/mNDqe85hqd2KyjQ0fw3Hkn8B//Afzf/+k3g02btOzz4Yea1VdV\n6beC2lqdq37yZNfJz957D/j4Yw3+xu0JDS+95NrBS0TnF4uIb2ZTsVgs8NGmfKa9XU8KtbU6/w2g\nWbnFAqxcCTz/vGP++ZYWzfo3bgTuvVenS3jgAb0ZyQUXaP3+xReBJ57Quvz27dpRe/XVjtLOxIla\n7vntb3W8flKSfru46Sb99rFli5ZzAO0juPBCrf9fdpnjKltA71VbWqpX8RJRYOtP7GQG3w8hIVpP\nN4I7oCN0goJ0mObo0Y73w8P1StrcXO2EzcwE5s7VTH3uXP1WkJenGfc11wCTJulsmP/7v1rTv/xy\nYN06zfJ/8Qtg/Hjgpz/VkTdffqnr/sMfNLB/9pnOsQPoUM077tDPDO++C8yfr1m/e4a/ezdvcEJk\nFqzB91NsrOf3Z8zQLNvZtdfq88UX65BLQDPrWbP0xJCWpu/l5OhcORERmvmL6AyYp08DTz8N5OcD\n99+vJZgnngBeew346CNd/1VXaZC/4QY9iXz8sX47qKrSjL6jQ08GNTXAX/+q1wB89pl+QwB0uyNH\nan/BoUPAP/2TY/9ffFG/fXQ3hTMRBRaWaPygoUEzek83FM/O1jKOMaLnwAEN7N98o8G1vl6DdFCQ\nXnhVWAi88ILeV/Y3vwGefVaD8sqVwJtvagdvXp5ObVxfr8NBT5zQk0hrqw4bnT1bTyT33ac/Dx+u\nJ5AVK7TT98QJ/TwtDXj9dS0LjRyppR9AO4AfecRxha5zOaihQcf2G7dSJKKz05/YyQA/SHz/vV5E\n5Txc8uRJrdM/8ohm+StW6PtHj+qUCG+/rQH9uef0ZBAVpSeEAweAe+7R8fyvvqrj+ydN0hr+rbdq\nRn/qlJaSWls1yM+YoaWgtDQt48TEAP/5nzqvzhdf6Mll40bd7urV2tEcF6c3LT9xQkcG/eUveqGY\n4d//HbjrLl3WYnGcGER0ZNHFF/vu70sUqPoVO6UXGzZskDFjxkh6ero8++yzHpd58MEHJT09XcaP\nHy87d+70uEwfNjWobd682W/bPn5cpK7O82c/+YlIVZXj9dGjIl98IfLllyJDhog0N4scOybS0SGy\nYIHIxo0iP/qRSEaGSGWlyKOPijz11GYBRD7+WMRqFbn5ZpHZs0VGjxZJTRWJihJZtEhk6FCR7GyR\nhASRiAiR664TGTFCZNgwkVtuEQkKEnnqKZE77hB5/XV9PXeuyL/8i0h+vsjXX4ts2SKybp1+9tVX\nIosXi3zzjciZMyJ2u8jvfy/y2GMi7e0in3/u2tZNmxxt7ehw/ay+vvu/nz+PnS+wfYNbf2Jnj7/Z\n3t4uaWlpsn//fmltbRWbzSZ79uxxWeaDDz6QadOmiYhIaWmp5Obmen0nB4PFixf7exfOWm2t5/e/\n+UakpMTxevHixbJ7twbNRYtEPvlEA+6ZMyK7dokUFop8+61IWJhIRYXI229rsJ4xQ+T55/UkMXSo\nrnPKFJEnnhCxWEReflkkNlZPFNOm6c/JyXriueMOkQsvFLn4YpH4eP3siitEhg8XSUzUdQcFidx1\nl8jUqSL/9V8iF1wgMnmyyL336rKrVoksWyayZo1IcLDIH/4gcv/9Im+9pW3YvVvkzTdFLr10sRw9\nKvLaa3qiPHBATyCvv64njfp6kb/9Tf8WHR36KCoSOX1apKVFlzW0t3f/d/WXwfhv82yYvX39iZ09\ndrKWlZUhPT0do/7RAzdnzhysX78eGRkZncsUFRVh3rx5AIDc3Fw0NDSgtrYWCSy6BjzjQi13l16q\nD2dGp/Cvf+36vs2mD0DLSEOG6CybgI4AAnT+nblzdZinMbpn4UJdLjpaSzQ2m9bqGxt1mZkzdT7+\nwkLtX4iL0wvEhg3TO3ktXar9EqtX60VeK1fqNQHPP6/9DXPm6MihmBhg+XJdzw9/qMvu2KGjh44c\n0auTL7hARyrZbFq2CgvT92Ji9JqDsDCguVmHxB45ovtbWakd0wcOaOlr5EjtF2lp0QvT7rlHr2mY\nNEn7UNLTtSM8NBS4/nqdhfSqq/Rvk5am+37TTbqdjg5dX1ub7sfmzdqeXbt0Wbtdr8H49lv9e11z\njZbNMjO1bJeSonMdZWRo+/7+dy27nTmjbfrrX7VD/sABx9xMERFaljt6VPtn9u3TfTh9WrdVUaHL\nRkdrv0p8vK4vOlrvfTB6tP5+eLi20WLRYcT19dq5b5QTRXR7jY36HB6uy1mtjn9Tzc2u/VPOfTp2\nu+e5pMizHgN8TU0NUlNTO1+npKRgx44dvS5TXV3NAH8e8jQHP6CBYdgw1/eMIaR33ul4z2rVYDln\njr4uKtJnYyTPo486lv3BD/Q//jPPuK7nuut0+GpoqGP59nZ9b+pUDc5GgGhq0mC6fLleQHbTTdpf\nEBOjQTo9XYNXXZ0G0W++0cC2bZuOfHr5Zb1OQUSDXlOT/v7VV+tnb7+tQTkuTi9ae/JJDXi7d+u0\n0zt36oll5069g9if/qSd0ZGROjdRUJCeNK++Wk9ckybpNRHh4TpSasQIPfmsW6cjtt57T0d17d+v\nr9es0bbW1+sJNSxMA3h2tvbljB6t+22xaMANDdXjtHu3/p1qajTwf/+9HoOjRzWox8bqz+HhGtRH\njdKZUMPC9KRkDL0NCdG2nD6t22ht1WMREaF/d2OZ5mY9kZ05o8f/1Ck9GbS06O+FhTlODMYJo65O\nl7Xbdb0rVugysbGOk0lbm27r+HF9327X7YaH67bOnNHX8fHa/o4O3Z8LLtD9bG3VfwsnT+q/zbAw\nfT55UtdhnNBDQ/V94+8YE+MY/Wa16sNYn/F/xG7XR3u7PsLDddvXXqt9Zl7TU3q/bt06WbhwYefr\nt956Sx544AGXZW655Rb59NNPO19ff/318rl7cfQfXzP44IMPPvg4+8eAlGiSk5NRZVyKCaCqqgop\nbhOtuy9TXV2N5OTkLusSjqAhIvKpHq9knTBhAioqKlBZWYnW1la8++67yM/Pd1kmPz8fa9asAQCU\nlpYiNjaW5RkiogDQYwYfEhKClStX4uabb4bdbseCBQuQkZGBV155BQBQUFCA6dOno7i4GOnp6YiM\njMSqVat8suNERNSLcy7u9FFfxtEPNhdddJFkZWVJdna2TJw4UUREjh8/LjfccIOMHj1abrzxRqnv\naeB1gJk/f77Ex8fLuHHjOt/rqT3Lli2T9PR0GTNmjHz44Yf+2OWz4ql9ixcvluTkZMnOzpbs7Gwp\nLi7u/Gwwte/gwYOSl5cnmZmZMnbsWCksLBQR8xy/7tpnluPX3NwsV1xxhdhsNsnIyJAnn3xSRLx3\n/AY0wPdlHP1gNGrUKDl+/LjLe48//rgsX75cRESeffZZ+elPf+qPXTsnn3zyiezcudMlAHbXnq+/\n/lpsNpu0trbK/v37JS0tTex2u1/2u688tW/JkiXyy1/+ssuyg619hw8flvLychEROXXqlFxyySWy\nZ88e0xzc+6RGAAADWUlEQVS/7tpnluMnItLU1CQiIm1tbZKbmytbt2712vEb0NkkncfRW63WznH0\nZiBuncbO1wPMmzcPf/zjH/2xW+dk8uTJiHObQay79qxfvx5z586F1WrFqFGjkJ6ejrKyMp/v89nw\n1D7Ac8f/YGtfYmIisrOzAQBRUVHIyMhATU2NaY5fd+0DzHH8ACAiIgIA0NraCrvdjri4OK8dvwEN\n8J7GyBsHZzCzWCy44YYbMGHCBLz22msA4HJxV0JCAmpra/25i/3WXXsOHTrkMpJqMB/Tl156CTab\nDQsWLEBDQwOAwd2+yspKlJeXIzc315THz2jflVdeCcA8x6+jowPZ2dlISEjA1KlTMXbsWK8dvwEN\n8BbnS9BMZNu2bSgvL8eGDRvwq1/9Clu3bnX53GKxmKrtvbVnMLZ10aJF2L9/P3bt2oWkpCQ89thj\n3S47GNrX2NiI2bNno7CwENHR0S6fmeH4NTY24vbbb0dhYSGioqJMdfyCgoKwa9cuVFdX45NPPsHm\nzZtdPu/P8RvQAN+XcfSDUdI/btA6YsQI3HbbbSgrK0NCQgKOHDkCADh8+DDiu5sHYJDorj19ve4h\n0MXHx3f+x1m4cGHn19zB2L62tjbMnj0bd999N2bNmgXAXMfPaN9dd93V2T4zHT9DTEwMZsyYgc8/\n/9xrx29AA3xfxtEPNqdPn8apU6cAAE1NTdi0aROysrKQn5+P1atXAwBWr17d+Q9xsOquPfn5+Vi7\ndi1aW1uxf/9+VFRU4IorrvDnrp6Tw4cPd/78/vvvIysrC8Dga5+IYMGCBcjMzMQjjzzS+b5Zjl93\n7TPL8Tt27Fhneam5uRl//vOfkZOT473jN6DdwyJSXFwsl1xyiaSlpcmyZcsGenMD7rvvvhObzSY2\nm03Gjh3b2abjx4/L9ddfPyiHSc6ZM0eSkpLEarVKSkqKvPnmmz2255lnnpG0tDQZM2aMbNy40Y97\n3jfu7XvjjTfk7rvvlqysLBk/frzMnDlTjhw50rn8YGrf1q1bxWKxiM1m6xwyuGHDBtMcP0/tKy4u\nNs3x+/LLLyUnJ0dsNptkZWXJc889JyI9x5OzaZ/PbvhBRES+xZtuExGZFAM8EZFJMcATEZkUAzwR\nkUkxwBMRmRQDPBGRSf0/A0/kLjsF/p8AAAAASUVORK5CYII=\n"
	}
	],
	"prompt_number": 34
	},
	{
	"cell_type": "code",
	"collapsed": false,
	"input": "",
	"language": "python",
	"metadata": {},
	"outputs": []
	}
	],
	"metadata": {}
	}
	]
	}