bamford/cumhist.ipynb

## cumhist.ipynb
{
 "worksheets": [
  {
   "cells": [
    {
     "metadata": {},
     "cell_type": "heading",
     "source": "How to avoid nasty vertical lines at the right edge of cumulative histograms",
     "level": 1
    },
    {
     "metadata": {},
     "cell_type": "code",
     "input": "%matplotlib inline",
     "prompt_number": 1,
     "outputs": [],
     "language": "python",
     "trusted": true,
     "collapsed": false
    },
    {
     "metadata": {},
     "cell_type": "code",
     "input": "from matplotlib import pyplot as plt\nimport numpy as np",
     "prompt_number": 2,
     "outputs": [],
     "language": "python",
     "trusted": true,
     "collapsed": false
    },
    {
     "metadata": {},
     "cell_type": "code",
     "input": "x = np.random.uniform(1, 12, 100)\ny = np.random.uniform(1, 15, 90)",
     "prompt_number": 3,
     "outputs": [],
     "language": "python",
     "trusted": true,
     "collapsed": false
    },
    {
     "metadata": {},
     "cell_type": "code",
     "input": "high = max(x.max(), y.max())\nrange=(0,high*1.10)",
     "prompt_number": 4,
     "outputs": [],
     "language": "python",
     "trusted": true,
     "collapsed": false
    },
    {
     "metadata": {},
     "cell_type": "code",
     "input": "plt.hist(x, cumulative=True, normed=True, histtype='step', bins=100, range=range)\nplt.hist(y, cumulative=True, normed=True, histtype='step', bins=100, range=range)\n_ = plt.axis(xmax=high*1.05)",
     "prompt_number": 5,
     "outputs": [
      {
       "png": "iVBORw0KGgoAAAANSUhEUgAAAXAAAAEDCAYAAAAsr19QAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xl4VPWhxvHvZCMJSYBIEsK+CQpIKLIJiAFrRcVHarBU\ne1FU9JZNQYSoyKaFErgiYqu40aLSokaieHGpFCK3QFllMUFZTYpAhpiQhKyQnPvHkQHMMAlkJieT\neT/PM49zcobzey3w9njmnN/PZhiGgYiIeB0/qwOIiMiVUYGLiHgpFbiIiJdSgYuIeCkVuIiIl1KB\ni4h4qYCqPnDixAmeffZZ9uzZw9atWyvtT0pKIisri9jYWLZv385zzz1H586dPRJWRETOq/IMfOPG\njQwfPpxL3S5eWFjIokWLmDp1KgkJCUydOtXtIUVEpLIqCzwhIYGwsLBL7n/uuecc78vLywkPD3dP\nMhERcanKSyjVVVZWxttvv80rr7zirkOKiIgLbvkSs6ysjHHjxjFv3jzatWvnjkOKiEhVjGpYv369\n0atXL8d2Tk6OkZ+fbxiGYRQWFhqjR4820tLSDMMwjOTkZKfH6NChgwHopZdeeul1Ga+4uLhLdnOV\nBf7VV18ZDz/8sNGyZUtj7ty5RnFxsTFt2jQjKSnJMAzD+PWvf2107tzZiI+PN+Lj440+ffo4H4hq\n/X9FnTBr1iyrI1SbsnqGt2T1lpyGoaxXylV3VnkNfNCgQQwaNOiinyUlJTner1q1qqpDiIiIB+hB\nHhERL6UCdyI+Pt7qCNWmrJ7hLVm9JScoqyfYfrrG4vmBbLZLPgwkIiLOuepOnYGLiHgpFbiIiJdS\ngYuIeCkVuIiIl1KBi4h4KRW4iIiXUoGLiHgpFbiIiJdSgYuIeCkVuIiIl1KBi4h4KRW4iIiXUoGL\niHgpFbiIiJdSgYuIeCkVuIiIl1KBi4h4KRW4iIiXUoGLiHgpFbiIiJdSgYuIeCkVuIiIl1KBi4h4\nKRW4iIiXUoGLiHgpFbiIiJeqssBPnDjBmDFj6NOnj9P9JSUlTJgwgfnz5/Pwww9z4MABt4cUEZHK\nqizwjRs3Mnz4cAzDcLp/8eLFtG3blqeeeorJkyfz8MMPuz2kiIhUVmWBJyQkEBYWdsn9n376KTfc\ncAMA3bp1Y/fu3Zw+fdp9CUVExKkaXwO32+2Eh4c7tiMiIrDb7TU9rIiIVCGgpgeIjo4mPz/fsZ2f\nn090dHRNDysi9cSGDbBwodUpas93zWeRF7qzVsa6ogLPzc0lICCA8PBw7rjjDjZv3szAgQPZu3cv\nPXr0uOQll9mzZzvex8fHEx8ffyXDi4gX2b4d/PxgzBirk9SOiekp3NV0DDFB7a7o13+/by8Z+/Y6\ntl1dz7AZl/p28icbNmzg7bff5osvvmDs2LE88cQTzJo1i8jISBITEykpKeHJJ58kNjaWgwcPMn36\ndDp27Fh5IJvtkl+Eikj9csstcPSo+f7HH+Ghh2D+fGszecrBnIM8v+F5R7+t/m41Gx7cQPeY7m45\nvqvurLLA3UUFLuI7wsIgNRUaNjS3W7c+/76+SU5PJmljEhP7TAQg0C+QEV1GEOgf6Jbjq8BFpFaF\nhcGJE+Y/66PXtr9G+sl0AA7kHCA0MJTk3yR7ZCxX3VnjLzFFxDdlZ0NhofN9FRW1m6W2Ldm6hDs7\n3UlsWCztm7SnTwvnDzp6mgpcRK5Ihw7QqBHYbJX3tW8PQUG1n8lTfsj/gTd3vomBeSZsL7Qzqvso\nukZ3tTSXClxEHHJy4Ouvq/fZoiLIyoLgYM9mqgv+lfkv3k9/n3u63APApL6TaN+kvcWpdA1cRC4w\ndy688YZ5dl2ViAhITgZ/f8/nssL/7v9fDuUcAmDniZ2UnC3hvRHv1XoOXQMXkWqpqIBRo+D5561O\nYr0p/5hC3xZ9aRLchMYNGvPLa39pdaRKVOAi4rNSv08lMy/T6b68kjym3zidzk0713Kq6lOBi4jP\nenj1w3SP6U54UHilfcM6DaN5eHMLUlWfClxEfNr/3PI/dIisxkX/Okgr8oiIeCkVuIiIl9IlFBEf\ns3IlzJvnfF9WFowfX7t55MqpwEV8zL59MGAAjB3rfH917gH3VmcrzvJd9neO7bLyMgvT1JwKXMQH\nNWsG3d0z26lXeT/tfcauGUvLiJYANA1tSpOQJhanunIqcBEfkJsLxcXm+4ICcw4TX3Sm/AzDrxnO\n8uHLrY7iFipwkXquogJiYyEy8vzPLnUNXLyLClyknjMMOHMGjh2zOom4m24jFBHxUjoDF6mHjhwx\nlzSD+r+4gi9TgYvUQ2++CZ99BnFx5vYTT1ibRzxDBS5STyUkwPTpVqcQT1KBi0i98sXBL3jw4wed\n7is6U8TIriNrOZHnqMBF6oGKCvj97yEvz9zes8dcmMEXHc0/yqA2g1h06yKn+68KuaqWE3mOClyk\nHjhzBpYtgxUrzO2774ZBg6zNZKWGgQ3r/Fze7qACF6kn/P1hZP25OiDVoPvARUS8lApcxEstWWLO\nadKoEURFQViY1YmktukSiogXKSiA0lLz/YEDMHmy+QIICrIul9We+OIJDuQcACAzL5M+zftYnKh2\nqMBFvERpKTRtCuEXrL/70ku+O7Pghd5Le4+5Q+Y67jDpHuMbc+WqwEXqsNxc+O6n9QeKiyEgALKz\nrc1UF+w+sZsXNr/g2M4pzuGW9rfQIqKFhalqX5UFvnbtWlJSUoiOjsZmszFz5syL9h85coSpU6fS\np08fdu/ezejRo7nllls8FljEl8yfD3//OzT/6Y64oUOtzVNX/Pvov/lP/n94sIf5wM6wTsN84rbB\nn7MZhmFcamdRURFxcXGkp6cTGBjIiBEjGDduHEOGDHF8ZuzYsVxzzTU8/vjj7Nq1i/vuu4/09PTK\nA9lsuBhKRJyYMsUs7ylTrE5S+9Lsaby+43Wn+/ba93J15NW8dudrtZyq9rnqTpdn4Js3b6ZNmzYE\nBgYCMGDAANasWXNRgTdr1gy73Q6A3W4nNjbWXblFxId9lfEVO47v4J4u91Ta175Je25qe5MFqeoW\nlwVut9sJv+Abk4iICEdZnzN58mTuvvtupkyZwvbt23nppZc8k1TEB1RUwKJFcPq0ub1pE4wYYW2m\n2vTmzjc5mn8UgG3HttE9pjuP93vc4lR1l8sCj4mJoaCgwLGdl5dHTEzMRZ8ZPXo0jzzyCCNHjiQ7\nO5sePXqwf/9+QkNDKx1v9uzZjvfx8fHEx8fXLL1IPVNQYM4g+PTT5vatt8Idd1ibqTY9tfYpHvrF\nQ4QGhtK7eW9u7XCr1ZFqXWpqKqnnJnOvQrWugaelpREUFMSIESMYP348PXr0ICAggPDwcPr06cPC\nhQu56aabOHv2LFFRUezfv5+oqKiLB9I1cJEq5eVB69bnJ6Wq7zLzMknZl+LYfmbdM2RMyqBpaFML\nU9UtrrrTZYGDeRdKcnIyUVFRBAUFMWPGDBITE4mMjCQxMZGNGzeyePFievbsyaFDh+jevTuPPfbY\nZYUQEZOvFfiSLUt46+u3iG8TD0BoYCjPD3meAD/d4XxOjQq8NkKIiMkXC/xgzkGW3LbE6ih1lqvu\n1FwoIiJeSv+dImKxRYtg3jzzfUUFOPn+v17599F/c7LwJADf2L8hOCDY4kTeSwUuYrHMTJg0yVxR\nByAkxNo8nnbbitvo26Ivgf7m8yUPxD1gcSLvpQIXqQPCwsyJquqj02Wn2XVil2P7TPkZVo5YSePg\nxhamqh9U4CLiUW/tfIukjUm0b9IegAGtBxAaWM+vE9USFbiIuN3R/KPklZi30vxQ8AO/7fbbSy4y\nLFdOBS4ibnfdq9cR0zAGfz9/ACb3m2xxovpJBS4iNVZeUU5WYZZju/RsKdsf3U5YkNZ58yQVuIjU\n2CvbXuGpfz5Fowbm8kCtGrUiyN+H13irJSpwEbkiZ8rPcLbiLAD5pfmM7z2eBbcssDiVb1GBi9Sy\ns2fN+74LC83tLVvg0UetzXQlWixqQX5pPjabDYA/DP6DxYl8j+ZCEallOTnQqhX86U/nfzZ0KHjb\nWijBfwjm1FOn9CSlh13xijwi4hkNGsCDD1qd4vK8ufNNJn420bHt7+ePv83fwkSiAheRarEX2pnQ\newLPDX4OMAv83OPwYg0VuEgteOEF8wXmhFXeMt9Jj6U9sBeayygWlBUwtf9UQgK9JLwPUIGL1IIj\nR2DcOHjoIXO7YUNr81RX2sk0Dk486DjTjm4YbXEiuZAKXKSWNG4MzZtbncK199PeZ85Xcy76WYuI\nFlohp47S74qIOHyX/R03tr6RiX3MLyvDgsJU3nWYfmdEfNxdK+/ix6IfAXOR4Yd+8RBdo7tanEqq\nQwUu4uNWf7eaDaM3OB7I6RbdzeJEUl0qcBEf8/nBz/nztj87tm3YGNh6oKPAxXuowEU8oKICXnkF\niorM7Z074ZprrM10zub/bCY8KJx7u90LwLT+01TeXkqP0ot4gN0ObdvChAnnf/bQQ9aUuGEYTPp8\nErkluQDsOrGLhGsTmBU/q/bDyGVz1Z0qcBE3+ewzSEkx3xcXwxdfmEVutfKKcgKeD2D58OWOnw1u\nO5hWjVpZmEqqS3OhiNSCjz82Zxi86SZze+RIa/NcyM/mx/1x91sdQ9xMBS7iRjfe6J1Tw4p3UoGL\n1EO7Tuzio28/AqDCqLA4jXiKn9UBRLxVfj507gxt2pivd981p4mtCz5I+4ANGRsA8/LJS0NfsjiR\neILOwEUug2GYK+oA5Oaar23bzu9v5ebvBdPsaaw9vPayf93WY1u5ud3NTB803b2BpE6pssDXrl1L\nSkoK0dHR2Gw2Zs6cWekzixYtws/Pj2PHjvHjjz/y1ltveSSsiNWmToUXXwS/n/7btVs38+zbU5Z9\nvYytx7bSs1nPy/p1XZp24daOt3ooldQVLgu8qKiIsWPHkp6eTmBgICNGjGDdunUMGTLE8Zl33nkH\nf39/Hn/8cQD27t3r2cQiHlZQAEePOt935Ai89hqMGVN7eYZ3Hs6U/lNqb0DxGi4LfPPmzbRp04bA\nQHMu4AEDBrBmzZqLCnzFihXceOONLFmyhKysLB46N+GxiJeaMgU++cSc/tWZceNqN4/IpbgscLvd\nTnh4uGM7IiIC+8+eTMjMzOTkyZNMnz6dQ4cOMXToUPbv369Hc8VrlZXBH/8Io0dbnUTENZcFHhMT\nQ0FBgWM7Ly+PmJiYiz4TERFBv379AOjQoQOFhYX85z//oXXr1pWON3v2bMf7+Ph44uPjaxBdRKT+\nSU1NJTU1tVqfdVng/fr1IyMjg7KyMoKCgti0aRPjx48nNzeXgIAAwsPDufnmmzl8+DAAp0+fpry8\nnGbNmjk93oUFLiLOfX38awrPFALwQ8EPNA+v48v4iFv9/OR2zpw5l/ysywIPDQ3l1Vdf5bHHHiMq\nKoq4uDgGDx5MYmIikZGRJCYmkpiYyOTJk5k7dy6ZmZksX76coKAgt/3LiPiSgtICer/Rm34t+zl+\n9tAv9L2SOKfJrER+ZvRoiI/37DXww7mHKT5TXOnnBWUF3PrureQ9lee5wcWraDIrkTokvzSfTi93\nonPTzk7392/Vv5YTibdSgYvPKyiA//ovKC01t/fuhQvulHWLsvIysouyAcgtziWiQQRp49LcO4j4\nHBW4+LyTJ2HLFvjrX8//bOBA947xzD+f4Y2db9AwsCEAXaK6uHcA8UkqcPFJ8+bBsmXm+zNnzId2\nhg517xhl5WWUV5QDcLrsNHOHzGVCnwlV/CqR6lOBi0/atw9+/3sYPtzcvtRTl1cqrySP6P+Jxs92\nfsLPC1fEEXEHFbj4rJgY6NjRM8cuLS+lUYNG2KfWgTXVpN5SgYvPyMw0H5MHcy5vd5u1fhYLNi0A\nzIWEW0a0dP8gIhdQgYtPyMgwz7bPTf1qs4GT2R5qJKswi6RfJvFIz0cACPDTXy/xLP0JE59QVgZt\n28KBA54dJ8g/iJDAEM8OIvITLakmIuKldAYu9da2beb93QB2fZco9ZAKXOqtRYvgxx+hUydz+6dF\no0TqDRW41BuFhTBnjvlgDsDOnTB7Ntx7r4fGKyvkrpV3UXK2BICDOQe5PvZ6zwwm4oRmI5R6Y98+\nuOkmePppc9tmg/vug+hoz4x3rOAY1716HR//9mPHz3rG9iQ0MNQzA4pP0myE4jOuugomT/bc8Vfs\nWcHKtJUAlJwtISQghIGt3Txxikg1qcBFLsP679fTtlFbftXhVwA0C3O++pRIbVCBi1ymuGZx3Nn5\nTqtjiOg+cBERb6UzcPFq48fD66+b7w0Devd27/FLz5YyY/0MSs+aqz1sPrr5ovUqRaykM3DxOseO\nmXObZGTAkSPw9ttQVATFxbBxo3vHyirM4vUdr9O+SXvaN2nPoz0f5farb3fvICJXSGfg4lX274eu\nXaF5c3PbZjMnqAoMdN8YGzI2sO7IOsCc1zuiQQSP99NTQFL3qMDFq5SWwjXXmOtWesqyr5dxquQU\nPZr1oFFwI54f/LznBhOpARW4+LyzFWd5c+ebjuvc6SfTGdd7HKN7jLY2mEgVdA1cfF5mXibTvpzG\n4dzDHM49zA0tb2BAqwFWxxKpks7ARYCmoU156baXrI4hcllU4FLnFBTAH/8IZ89W3nfyZO3nEamr\nVOBS5xw8CH/9K0yaVHlf06bnV5KvyuHcw6TZ06r83InTJy4voEgdoQKXOikmBqZNq9kxnl33LN9m\nf0uLiBZVfvaeLvfUbDARC6jApd4yMJjafyr3XuehCcFFLFZlga9du5aUlBSio6Ox2WzMnDnT6edW\nrFjBqFGjOH36NKGhmg9ZLs/s2bBqlfm+uBiaNLE0johXcFngRUVFjB07lvT0dAIDAxkxYgTr1q1j\nyJAhF31u37597Nu3z6NBpX7bsQPGjIH4eHM7KsrSOCJeweV94Js3b6ZNmzYE/vSc8oABA1izZs1F\nnykqKmLhwoXMmjXLcynFJ7RrB927m6/YWKvTiNR9Ls/A7XY74eHhju2IiAjsP1vee/r06cyaNctR\n8lo2TarLbj9/q2BJiXuOebLwJGcqzEUxi84UueegInWUywKPiYmhoKDAsZ2Xl0dMTIxj++jRo5w6\ndYqVK1c6fvbiiy9y2223cf31lRd3nT17tuN9fHw88ef+e1l8zt690LPn+Uslfn7QrIaL2xzJPcLV\nL19NdENzEUybzcYT/Z6oYVKR2pWamkpqamq1PutyUeOioiLi4uJIS0sjKCiIESNGMH78eHr06EFA\nQMBFZ+cAfn5+l/wSU4say4W2bYNx48x/1kSFUeGYw+Tb7G+5b9V97Buv72Ok/rjiRY1DQ0N59dVX\neeyxx4iKiiIuLo7BgweTmJhIZGQkiYmJAGRnZ7N06VJsNhsLFizg0Ucfpfm5+T5FPOiR1Y/wzp53\n8PfzB6B/q/4WJxKpPS7PwN06kM7Afd727eengT1yBD77rOZn4Pd8cA+/6fIb7umqB3GkfrriM3AR\nd5o5E8rKoFUrc/uBBy7/GBmnMrju1escX1SWlZcx5hdj3JhSxHuowMVj8vJg0aLzd5p89x28/DLc\nXoMVyfJL82nVqBXbH9kOmGcnwQHBbkgr4n1U4OIx+/bBsmXw+9+b22PGQN++l3+cqf+Yyoq9KwA4\nU3GGto3bEhIY4sakIt5J18DFrV59Fc7dAZWdbS42vHlzzY756/d+zbCrh3Hb1bcBENEggrCgsJod\nVMRL6Bq41JqUFOjVC+LizO2OHd1z3MiQSJqH684mkQupwMXt4uPhV7+yOoVI/acClzrnxOkT/Db5\nt5ytML/9TD+ZzoM9HrQ4lUjdowKXGklPh4EDobzc3C4shJrOa3a84Dg/FPzAX+76CwA2bPRu0buG\nSUXqHxW4XLbSUnPdSoDDh6FTJ/jiC3Pbzw9+NsNCtfx565/5/NDnAOSV5BHRIIKBrQe6KbFI/aQC\nl8uWkAAbNkBQkLl9++3QqFHNjvn5oc/pEdODPi36ANCmcZsaphSp/1Tg4tTRo+bLmePH4aOP4Gfr\nelyW7KJspn05zXGd++vjX/Noz0e5s/OdV35QER+jAhen7r8fTpyAiIjK+4KDoXXrmh3/SO4R1n+/\nnjnxcwD4VYdfEd82vmYHFfExKnAfUlgIH38MFRVVf/aHH+CNN2DQoCsfL6c4h7kb5jrOsi+UVZhF\n09Cm3B93/5UPIOLjVOA+5F//gilT4Oabq/5sv35w9dU1G2//j/tJ3pfsdFGF9k3a89AvHqrZACI+\nTo/S13OrVsGePeb7gwfNx9s//9x9xz9VcoqXt7xMuVFead/R/KN8Y/+Gf4/5t/sGFPExrrrT5aLG\n4v2Sks5/GdmxIzz2mHuP/439G97Y+YbTfS0jWjJtwDT3DigiDrqE4gMeeeTKZgGsrtaNWjM7frbn\nBhARp3QGLiLipXQGXg888QRs2uR8X1oaBLj5d3nbD9vYl20uHLz/x/3uPbiIVJsKvB746iuYOBGu\nuabyvoAA+MUv3Dte4tpE/P38iQ2LBeC+6+5z7wAiUi0qcC9x993m4+vOnDplXuO+9tray/P0wKcZ\n0q4Gj2KKSI2pwL3EoUOQnAzdulXeFxAAjRvXfiYRsZYK3Is0aQJNm7rveN+f+p6j+ZeY8MSFUyWn\n3BdCRK6YCtxCZ87A3r1QneebiorcP/6olFGOqVsvR3iDcNo00myBIlZTgVto9Wpzpfb27av+7FVX\nQXS0e8cvryhn6bCl9G/V370HFpFaoQKvZZs3Q2am+X7TJnPtyPfeszaTiHgnFXgte/BBaNfu/DSt\nw4Z5drzCskLySvOc7isrL/Ps4CLiUSrwWmYYsHgxdO5cO+Pd/rfbSbOnEeQfVGmfv58/TUPd+K2o\niNQqFXg9VHK2xDF7WdGZItbct4a+LT04GYqIWKJaBb527VpSUlKIjo7GZrMxc+bMi/YnJSWRlZVF\nbGws27dv57nnnqNzbZ1iykU2ZGxg8PLBjjPuAL8AGgfrJnGR+qjKAi8qKmLs2LGkp6cTGBjIiBEj\nWLduHUMuWBCxsLCQRYsWAfD+++8zdepUVq9e7bnUckkFpQUM7TiUNfetsTqKiHhYlbMRbt68mTZt\n2hAYGAjAgAEDWLPm4nJ47rnnHO/Ly8sJDw93c0xxZdjfhhEyN4SQuSH8+r1fX/Z93SLinao8A7fb\n7RcVckREBHa73elny8rKePvtt3nllVfcl1CqdPz0cdaOWkvP2J4ATr+wFJH6p8oCj4mJoaCgwLGd\nl5dHTExMpc+VlZUxbtw45s2bR7t27dybUqoUHBBMSGCI1TFEpBZVWeD9+vUjIyODsrIygoKC2LRp\nE+PHjyc3N5eAgADCw8MpKipi/PjxTJ06lS5duvDhhx+SkJBQ6VizZ892vI+Pjyc+Pt6d/y4iIl4v\nNTWV1NTUan22Wosar127luTkZKKioggKCmLGjBkkJiZy1VVXMW3aNO6++27S09OJjTXnhy4qKmLL\nli0XD6RFjQHz/u/Vq2t2H/jO4zu5P+V+DMz/PQ/lHGLrI1vpHtPdTSlFpK5w1Z1ald7DzpyBb789\nP2HVnXfCP/5RswL/+NuPWbxlMX+67U+Aeatgp6s6YbPZ3JBYROoSV92pB3k8LCUFHn0UWrc2t6Oi\nzImpqrL7xG4mfDbB6W/cj8U/cm3Ta+ka3dXNaUXEm6jAPWDrVjj60zTbmzbBHXfAihWXd4yDOQfx\ns/kx9+a5Tve3a6wvikV8nQrcA373O3PCqrAwc3v48Or9uhnrZrAraxcAxwuO065JOwa2HuihlCLi\n7VTgbnDsGLz77vntnBz47DPo2PHyjrPq21WM7TXWsVjCtVG1uMiliHgdFbgbrF0Ly5ebl0oAxo6F\nli2v7FiD2w7WtW0RqRYVuAtJSeZiwlXZvx+uvx4WLPB8JhGRc1TgLrz8MkycaC4m7EqvXtBXs7WK\nSC1TgVfhd7+78sshIiKepAKvZX/b+zf2/7jf6T57ofNJwkREnKlyOllf8vnn0KbN+deJE9CggXvH\nmPPVHLKLsp3um9R3Eu2bVGOJehERdAZOeTlUVJjvDx+G/v1h/nxzu0ED88nJmrAX2ln5zUrHE5W5\nxblM7DORzk21YpGI1IzPF3jXrnDwIJybRmTaNPPs212+PPQlf972Z4Z2GArAA3EP0KpRK/cNICI+\ny+cL/NQp87H3Zs08N0av5r146baXPDeAiPgkXQMXEfFSPn8G7gk7j+/kh/wfHO9FRDxBBe4BI5NH\n0iqiFQ2DGgJw9zV3W5xIROojFbgHGIbB63e+TsfIy5zNSkTkMugauIiIl/K5M/C0NHj11fPb+fnV\n+3Vnys9c8gnKnystL72CZCIil8fnCvyrr2D3bhg50tx++WWIjq761638ZiUTP5tIi4gWVX62aWhT\nIkMia5hURMQ1nytwgOuugwkTqv5cfmk+p8tOA+YTlXdfezfL7lrm4XQiItXjkwVeXT1f60lBWQH+\nNn8AxvUeZ3EiEZHzfL7Az5Sf4WzFWaf7is8Ws+PRHbSM0HyyIlL3+HyBd32lK5l5mdjOTYZygeCA\nYBoGNrQglYhI1XyiwDt1ggMHzm9Pn37+fX5pPt9P+p5mYR6cDEVExAN84j7wnBw4eRIMAz5IS+aF\nkBBC5pqv3JJcggOCrY4oInLZfOIM/EL2Qjujuo/ipaHm7ID+fv4E+QdZnEpE5PL5XIEDBPgFEBIY\nYnUMEZEa8YlLKCIi9VGVZ+Br164lJSWF6OhobDYbM2fOvGh/SUkJTz75JC1btuTAgQM89dRTXH31\n1R4LfCkffgjHj5vv08o+5aOiqY59OffCwL+BfwDkFOfwmy6/qfV8IiLuZjPOLdboRFFREXFxcaSn\npxMYGMiIESMYN24cQ4YMcXxm/vz5BAQE8OSTT/LNN98wbtw4NmzYUHkgmw0XQ9VYkyaQkADBwbA7\nZAlZAVu4vugZAEKC4fFJEGA+j0PrRq0JbxB+yWOlpqYSHx/vsazupKye4S1ZvSUnKOuVctWdLi+h\nbN68mTZt2hAYGAjAgAEDWLNmzUWf+fTTT7nhhhsA6NatG7t37+b06dPuyO3S3r0wefL5V/4vf8c3\nfQayq+c6NDSIAAAGTklEQVRAfmi5mDtuiubvL3Xl7y91ZVlSV+Jiu9I12ny5Km8wf/O8hbJ6hrdk\n9ZacoKye4PISit1uJzz8fNlFRERgt9ur9ZmwsDA3R73Y2rWwcycMH25uBzf5lBnxK2gUHAHANU2v\n8ej4IiJWc1ngMTExFBQUOLbz8vKIiYm56DPR0dHkXzAna35+PtGXmN4vZvKdNcl6kaIiiL4R1v30\n/E3Z6dMMbD2ARsGN3DaGiEidZrhQWFhodOzY0SgtLTUMwzASEhKMdevWGTk5OUZ+fr5hGIYxf/58\nY8GCBYZhGMaePXuMQYMGOT1Whw4dDEAvvfTSS6/LeMXFxV2yo11+iQnmXSjJyclERUURFBTEjBkz\nSExMJDIyksTERMddKLGxsRw8eJDp06fTsaOWEhMR8bQqC1xEROomPcgjIuKlVOAiIl5KBS4i4qVU\n4CIiXkoFLiLipTw+nWxVk2HVJYcOHWLGjBn07NmTo0ePctVVVzFjxgyrY11ScXExffv25dZbb2Xh\nwoVWx7mk7777juXLl9O4cWP+7//+j5kzZ9K7d2+rYzm1cOFCMjIyiIqK4vvvv+f11193TCVhtRMn\nTvDss8+yZ88etm7dCkBOTg5PP/007du358CBA8ybN++SD9LVJmdZJ0+eTMOGDQkLC2P37t0sXry4\n0oOBVnCW9Zy5c+eyePFiTp48aVG6Krh6kKemzj0IVFZWZhiG+SDQP//5T08OWSPbtm0zVq9e7dju\n0qWLsWPHDgsTufbEE08YDzzwgPHkk09aHeWSzp49a9x4442OPwMnTpwwTp48aXEq544fP25ERkY6\ntu+66y7jvffeszDRxZKTk41PPvnE6NWrl+Nn//3f/2188MEHhmEYxieffGKMGjXKqngXcZb12Wef\ndbxPSkoyJk6caEW0SpxlNQzDWL9+vTFlyhQjKirKomRV8+gllOpMhlWX9OrVizvvPP+4f0VFhcfn\ndLlS7777LgMHDqRdu3ZWR3Fp27ZtGIbB0qVLWbBgASkpKTRt2tTqWE6FhoYSFBREXl4eANnZ2TRr\nVnfWSk1ISKj05/HCyeT69+9fZ/5+Ocv6/PPPO96Xl5dfNIeSlZxlzcrK4v3332fixIkenUW1pjx6\nCaU6k2HVVSkpKQwdOpROnTpZHaWS9PR09u3bx9y5c9m9e7fVcVzKyMhg+/btJCcnExMTw+jRo2nQ\noAEPPvig1dEqiYiIYOHChYwcOZLmzZvTs2dPBg0aZHUsly78OxYREUFubi4VFRX4+dXdr7dOnTrF\nl19+yapVq6yO4lRFRQXTp0/nhRdeIDc31+o4Lnn0d7k6k2HVRevXr+err77ixRdftDqKUx999BHB\nwcEkJSWxceNGtm7dypIlS6yO5VRERAStW7d2/L7379/f6XzxdcGuXbtYuHAhn376KcuWLSMoKIik\npCSrY7kUHR3t+DuWn59PkyZN6nR55+XlMWHCBP7yl7/QuHFjq+M49fXXXxMYGMjSpUtZunQpxcXF\nLFiwgIMHD1odrRKPnoH369ePjIwMysrKCAoKYtOmTYwfP96TQ9bYmjVr+Ne//sXixYs5duwYmZmZ\n9OvXz+pYF3nmmWcc74uLiyksLOSxxx6zMNGl9e3bl4KCAkpLS2nQoAEZGRl07tzZ6lhOHTt2jMjI\nSEcBNmvWjCNHjlicyrU77riDTZs2cc8997Bx40aGDRtmdaRLys7OZtKkSSxYsIDmzZvz4YcfkpCQ\nYHWsSq6//nquv/56AL7//nveeustpk2bZnEq5/xnz54921MHDwwMpEuXLrzwwgts2bKFFi1a8MAD\nD3hquBrbsWMHw4YNwzAM/vrXv/LOO+/QoUMH4uLirI7m1KpVq0hOTiYzM5OQkBC6detmdaRKQkJC\n6NChAy+//DK7d+8mKyuLOXPmEBBQ99bT7tChAzt27GDdunVs3bqVHTt28Ic//IGIiAirowGwYcMG\n3nnnHXbv3k1JSQm9e/fmpptu4rXXXmPPnj1s2rSJhQsX0rBhQ6ujOs06ePBgsrKy+PTTT1m+fDnf\nfvst9957r9VRnWYNCAjg0KFDvPLKK2zZsgUwi72u3JF0jiazEhHxUnX3YpmIiLikAhcR8VIqcBER\nL6UCFxHxUipwEREvpQIXEfFSKnARES+lAhcR8VL/D02xEVmWZuiZAAAAAElFTkSuQmCC\n",
       "text": "<matplotlib.figure.Figure at 0x107dd0a50>",
       "output_type": "display_data",
       "metadata": {}
      }
     ],
     "language": "python",
     "trusted": true,
     "collapsed": false
    }
   ],
   "metadata": {}
  }
 ],
 "metadata": {
  "name": "",
  "signature": "sha256:0795d55207b5da4a51d307bb6088633cdbc813297e7a9f5ee8bba7ecb1765be8",
  "gist_id": "29457bec83da6865f5c2"
 },
 "nbformat": 3
}
	{
	"worksheets": [
	{
	"cells": [
	{
	"metadata": {},
	"cell_type": "heading",
	"source": "How to avoid nasty vertical lines at the right edge of cumulative histograms",
	"level": 1
	},
	{
	"metadata": {},
	"cell_type": "code",
	"input": "%matplotlib inline",
	"prompt_number": 1,
	"outputs": [],
	"language": "python",
	"trusted": true,
	"collapsed": false
	},
	{
	"metadata": {},
	"cell_type": "code",
	"input": "from matplotlib import pyplot as plt\nimport numpy as np",
	"prompt_number": 2,
	"outputs": [],
	"language": "python",
	"trusted": true,
	"collapsed": false
	},
	{
	"metadata": {},
	"cell_type": "code",
	"input": "x = np.random.uniform(1, 12, 100)\ny = np.random.uniform(1, 15, 90)",
	"prompt_number": 3,
	"outputs": [],
	"language": "python",
	"trusted": true,
	"collapsed": false
	},
	{
	"metadata": {},
	"cell_type": "code",
	"input": "high = max(x.max(), y.max())\nrange=(0,high*1.10)",
	"prompt_number": 4,
	"outputs": [],
	"language": "python",
	"trusted": true,
	"collapsed": false
	},
	{
	"metadata": {},
	"cell_type": "code",
	"input": "plt.hist(x, cumulative=True, normed=True, histtype='step', bins=100, range=range)\nplt.hist(y, cumulative=True, normed=True, histtype='step', bins=100, range=range)\n_ = plt.axis(xmax=high*1.05)",
	"prompt_number": 5,
	"outputs": [
	{
	"png": "iVBORw0KGgoAAAANSUhEUgAAAXAAAAEDCAYAAAAsr19QAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xl4VPWhxvHvZCMJSYBIEsK+CQpIKLIJiAFrRcVHarBU\ne1FU9JZNQYSoyKaFErgiYqu40aLSokaieHGpFCK3QFllMUFZTYpAhpiQhKyQnPvHkQHMMAlkJieT\neT/PM49zcobzey3w9njmnN/PZhiGgYiIeB0/qwOIiMiVUYGLiHgpFbiIiJdSgYuIeCkVuIiIl1KB\ni4h4qYCqPnDixAmeffZZ9uzZw9atWyvtT0pKIisri9jYWLZv385zzz1H586dPRJWRETOq/IMfOPG\njQwfPpxL3S5eWFjIokWLmDp1KgkJCUydOtXtIUVEpLIqCzwhIYGwsLBL7n/uuecc78vLywkPD3dP\nMhERcanKSyjVVVZWxttvv80rr7zirkOKiIgLbvkSs6ysjHHjxjFv3jzatWvnjkOKiEhVjGpYv369\n0atXL8d2Tk6OkZ+fbxiGYRQWFhqjR4820tLSDMMwjOTkZKfH6NChgwHopZdeeul1Ga+4uLhLdnOV\nBf7VV18ZDz/8sNGyZUtj7ty5RnFxsTFt2jQjKSnJMAzD+PWvf2107tzZiI+PN+Lj440+ffo4H4hq\n/X9FnTBr1iyrI1SbsnqGt2T1lpyGoaxXylV3VnkNfNCgQQwaNOiinyUlJTner1q1qqpDiIiIB+hB\nHhERL6UCdyI+Pt7qCNWmrJ7hLVm9JScoqyfYfrrG4vmBbLZLPgwkIiLOuepOnYGLiHgpFbiIiJdS\ngYuIeCkVuIiIl1KBi4h4KRW4iIiXUoGLiHgpFbiIiJdSgYuIeCkVuIiIl1KBi4h4KRW4iIiXUoGL\niHgpFbiIiJdSgYuIeCkVuIiIl1KBi4h4KRW4iIiXUoGLiHgpFbiIiJdSgYuIeCkVuIiIl1KBi4h4\nKRW4iIiXUoGLiHgpFbiIiJeqssBPnDjBmDFj6NOnj9P9JSUlTJgwgfnz5/Pwww9z4MABt4cUEZHK\nqizwjRs3Mnz4cAzDcLp/8eLFtG3blqeeeorJkyfz8MMPuz2kiIhUVmWBJyQkEBYWdsn9n376KTfc\ncAMA3bp1Y/fu3Zw+fdp9CUVExKkaXwO32+2Eh4c7tiMiIrDb7TU9rIiIVCGgpgeIjo4mPz/fsZ2f\nn090dHRNDysi9cSGDbBwodUpas93zWeRF7qzVsa6ogLPzc0lICCA8PBw7rjjDjZv3szAgQPZu3cv\nPXr0uOQll9mzZzvex8fHEx8ffyXDi4gX2b4d/PxgzBirk9SOiekp3NV0DDFB7a7o13+/by8Z+/Y6\ntl1dz7AZl/p28icbNmzg7bff5osvvmDs2LE88cQTzJo1i8jISBITEykpKeHJJ58kNjaWgwcPMn36\ndDp27Fh5IJvtkl+Eikj9csstcPSo+f7HH+Ghh2D+fGszecrBnIM8v+F5R7+t/m41Gx7cQPeY7m45\nvqvurLLA3UUFLuI7wsIgNRUaNjS3W7c+/76+SU5PJmljEhP7TAQg0C+QEV1GEOgf6Jbjq8BFpFaF\nhcGJE+Y/66PXtr9G+sl0AA7kHCA0MJTk3yR7ZCxX3VnjLzFFxDdlZ0NhofN9FRW1m6W2Ldm6hDs7\n3UlsWCztm7SnTwvnDzp6mgpcRK5Ihw7QqBHYbJX3tW8PQUG1n8lTfsj/gTd3vomBeSZsL7Qzqvso\nukZ3tTSXClxEHHJy4Ouvq/fZoiLIyoLgYM9mqgv+lfkv3k9/n3u63APApL6TaN+kvcWpdA1cRC4w\ndy688YZ5dl2ViAhITgZ/f8/nssL/7v9fDuUcAmDniZ2UnC3hvRHv1XoOXQMXkWqpqIBRo+D5561O\nYr0p/5hC3xZ9aRLchMYNGvPLa39pdaRKVOAi4rNSv08lMy/T6b68kjym3zidzk0713Kq6lOBi4jP\nenj1w3SP6U54UHilfcM6DaN5eHMLUlWfClxEfNr/3PI/dIisxkX/Okgr8oiIeCkVuIiIl9IlFBEf\ns3IlzJvnfF9WFowfX7t55MqpwEV8zL59MGAAjB3rfH917gH3VmcrzvJd9neO7bLyMgvT1JwKXMQH\nNWsG3d0z26lXeT/tfcauGUvLiJYANA1tSpOQJhanunIqcBEfkJsLxcXm+4ICcw4TX3Sm/AzDrxnO\n8uHLrY7iFipwkXquogJiYyEy8vzPLnUNXLyLClyknjMMOHMGjh2zOom4m24jFBHxUjoDF6mHjhwx\nlzSD+r+4gi9TgYvUQ2++CZ99BnFx5vYTT1ibRzxDBS5STyUkwPTpVqcQT1KBi0i98sXBL3jw4wed\n7is6U8TIriNrOZHnqMBF6oGKCvj97yEvz9zes8dcmMEXHc0/yqA2g1h06yKn+68KuaqWE3mOClyk\nHjhzBpYtgxUrzO2774ZBg6zNZKWGgQ3r/Fze7qACF6kn/P1hZP25OiDVoPvARUS8lApcxEstWWLO\nadKoEURFQViY1YmktukSiogXKSiA0lLz/YEDMHmy+QIICrIul9We+OIJDuQcACAzL5M+zftYnKh2\nqMBFvERpKTRtCuEXrL/70ku+O7Pghd5Le4+5Q+Y67jDpHuMbc+WqwEXqsNxc+O6n9QeKiyEgALKz\nrc1UF+w+sZsXNr/g2M4pzuGW9rfQIqKFhalqX5UFvnbtWlJSUoiOjsZmszFz5syL9h85coSpU6fS\np08fdu/ezejRo7nllls8FljEl8yfD3//OzT/6Y64oUOtzVNX/Pvov/lP/n94sIf5wM6wTsN84rbB\nn7MZhmFcamdRURFxcXGkp6cTGBjIiBEjGDduHEOGDHF8ZuzYsVxzzTU8/vjj7Nq1i/vuu4/09PTK\nA9lsuBhKRJyYMsUs7ylTrE5S+9Lsaby+43Wn+/ba93J15NW8dudrtZyq9rnqTpdn4Js3b6ZNmzYE\nBgYCMGDAANasWXNRgTdr1gy73Q6A3W4nNjbWXblFxId9lfEVO47v4J4u91Ta175Je25qe5MFqeoW\nlwVut9sJv+Abk4iICEdZnzN58mTuvvtupkyZwvbt23nppZc8k1TEB1RUwKJFcPq0ub1pE4wYYW2m\n2vTmzjc5mn8UgG3HttE9pjuP93vc4lR1l8sCj4mJoaCgwLGdl5dHTEzMRZ8ZPXo0jzzyCCNHjiQ7\nO5sePXqwf/9+QkNDKx1v9uzZjvfx8fHEx8fXLL1IPVNQYM4g+PTT5vatt8Idd1ibqTY9tfYpHvrF\nQ4QGhtK7eW9u7XCr1ZFqXWpqKqnnJnOvQrWugaelpREUFMSIESMYP348PXr0ICAggPDwcPr06cPC\nhQu56aabOHv2LFFRUezfv5+oqKiLB9I1cJEq5eVB69bnJ6Wq7zLzMknZl+LYfmbdM2RMyqBpaFML\nU9UtrrrTZYGDeRdKcnIyUVFRBAUFMWPGDBITE4mMjCQxMZGNGzeyePFievbsyaFDh+jevTuPPfbY\nZYUQEZOvFfiSLUt46+u3iG8TD0BoYCjPD3meAD/d4XxOjQq8NkKIiMkXC/xgzkGW3LbE6ih1lqvu\n1FwoIiJeSv+dImKxRYtg3jzzfUUFOPn+v17599F/c7LwJADf2L8hOCDY4kTeSwUuYrHMTJg0yVxR\nByAkxNo8nnbbitvo26Ivgf7m8yUPxD1gcSLvpQIXqQPCwsyJquqj02Wn2XVil2P7TPkZVo5YSePg\nxhamqh9U4CLiUW/tfIukjUm0b9IegAGtBxAaWM+vE9USFbiIuN3R/KPklZi30vxQ8AO/7fbbSy4y\nLFdOBS4ibnfdq9cR0zAGfz9/ACb3m2xxovpJBS4iNVZeUU5WYZZju/RsKdsf3U5YkNZ58yQVuIjU\n2CvbXuGpfz5Fowbm8kCtGrUiyN+H13irJSpwEbkiZ8rPcLbiLAD5pfmM7z2eBbcssDiVb1GBi9Sy\ns2fN+74LC83tLVvg0UetzXQlWixqQX5pPjabDYA/DP6DxYl8j+ZCEallOTnQqhX86U/nfzZ0KHjb\nWijBfwjm1FOn9CSlh13xijwi4hkNGsCDD1qd4vK8ufNNJn420bHt7+ePv83fwkSiAheRarEX2pnQ\newLPDX4OMAv83OPwYg0VuEgteOEF8wXmhFXeMt9Jj6U9sBeayygWlBUwtf9UQgK9JLwPUIGL1IIj\nR2DcOHjoIXO7YUNr81RX2sk0Dk486DjTjm4YbXEiuZAKXKSWNG4MzZtbncK199PeZ85Xcy76WYuI\nFlohp47S74qIOHyX/R03tr6RiX3MLyvDgsJU3nWYfmdEfNxdK+/ix6IfAXOR4Yd+8RBdo7tanEqq\nQwUu4uNWf7eaDaM3OB7I6RbdzeJEUl0qcBEf8/nBz/nztj87tm3YGNh6oKPAxXuowEU8oKICXnkF\niorM7Z074ZprrM10zub/bCY8KJx7u90LwLT+01TeXkqP0ot4gN0ObdvChAnnf/bQQ9aUuGEYTPp8\nErkluQDsOrGLhGsTmBU/q/bDyGVz1Z0qcBE3+ewzSEkx3xcXwxdfmEVutfKKcgKeD2D58OWOnw1u\nO5hWjVpZmEqqS3OhiNSCjz82Zxi86SZze+RIa/NcyM/mx/1x91sdQ9xMBS7iRjfe6J1Tw4p3UoGL\n1EO7Tuzio28/AqDCqLA4jXiKn9UBRLxVfj507gxt2pivd981p4mtCz5I+4ANGRsA8/LJS0NfsjiR\neILOwEUug2GYK+oA5Oaar23bzu9v5ebvBdPsaaw9vPayf93WY1u5ud3NTB803b2BpE6pssDXrl1L\nSkoK0dHR2Gw2Zs6cWekzixYtws/Pj2PHjvHjjz/y1ltveSSsiNWmToUXXwS/n/7btVs38+zbU5Z9\nvYytx7bSs1nPy/p1XZp24daOt3ooldQVLgu8qKiIsWPHkp6eTmBgICNGjGDdunUMGTLE8Zl33nkH\nf39/Hn/8cQD27t3r2cQiHlZQAEePOt935Ai89hqMGVN7eYZ3Hs6U/lNqb0DxGi4LfPPmzbRp04bA\nQHMu4AEDBrBmzZqLCnzFihXceOONLFmyhKysLB46N+GxiJeaMgU++cSc/tWZceNqN4/IpbgscLvd\nTnh4uGM7IiIC+8+eTMjMzOTkyZNMnz6dQ4cOMXToUPbv369Hc8VrlZXBH/8Io0dbnUTENZcFHhMT\nQ0FBgWM7Ly+PmJiYiz4TERFBv379AOjQoQOFhYX85z//oXXr1pWON3v2bMf7+Ph44uPjaxBdRKT+\nSU1NJTU1tVqfdVng/fr1IyMjg7KyMoKCgti0aRPjx48nNzeXgIAAwsPDufnmmzl8+DAAp0+fpry8\nnGbNmjk93oUFLiLOfX38awrPFALwQ8EPNA+v48v4iFv9/OR2zpw5l/ysywIPDQ3l1Vdf5bHHHiMq\nKoq4uDgGDx5MYmIikZGRJCYmkpiYyOTJk5k7dy6ZmZksX76coKAgt/3LiPiSgtICer/Rm34t+zl+\n9tAv9L2SOKfJrER+ZvRoiI/37DXww7mHKT5TXOnnBWUF3PrureQ9lee5wcWraDIrkTokvzSfTi93\nonPTzk7392/Vv5YTibdSgYvPKyiA//ovKC01t/fuhQvulHWLsvIysouyAcgtziWiQQRp49LcO4j4\nHBW4+LyTJ2HLFvjrX8//bOBA947xzD+f4Y2db9AwsCEAXaK6uHcA8UkqcPFJ8+bBsmXm+zNnzId2\nhg517xhl5WWUV5QDcLrsNHOHzGVCnwlV/CqR6lOBi0/atw9+/3sYPtzcvtRTl1cqrySP6P+Jxs92\nfsLPC1fEEXEHFbj4rJgY6NjRM8cuLS+lUYNG2KfWgTXVpN5SgYvPyMw0H5MHcy5vd5u1fhYLNi0A\nzIWEW0a0dP8gIhdQgYtPyMgwz7bPTf1qs4GT2R5qJKswi6RfJvFIz0cACPDTXy/xLP0JE59QVgZt\n28KBA54dJ8g/iJDAEM8OIvITLakmIuKldAYu9da2beb93QB2fZco9ZAKXOqtRYvgxx+hUydz+6dF\no0TqDRW41BuFhTBnjvlgDsDOnTB7Ntx7r4fGKyvkrpV3UXK2BICDOQe5PvZ6zwwm4oRmI5R6Y98+\nuOkmePppc9tmg/vug+hoz4x3rOAY1716HR//9mPHz3rG9iQ0MNQzA4pP0myE4jOuugomT/bc8Vfs\nWcHKtJUAlJwtISQghIGt3Txxikg1qcBFLsP679fTtlFbftXhVwA0C3O++pRIbVCBi1ymuGZx3Nn5\nTqtjiOg+cBERb6UzcPFq48fD66+b7w0Devd27/FLz5YyY/0MSs+aqz1sPrr5ovUqRaykM3DxOseO\nmXObZGTAkSPw9ttQVATFxbBxo3vHyirM4vUdr9O+SXvaN2nPoz0f5farb3fvICJXSGfg4lX274eu\nXaF5c3PbZjMnqAoMdN8YGzI2sO7IOsCc1zuiQQSP99NTQFL3qMDFq5SWwjXXmOtWesqyr5dxquQU\nPZr1oFFwI54f/LznBhOpARW4+LyzFWd5c+ebjuvc6SfTGdd7HKN7jLY2mEgVdA1cfF5mXibTvpzG\n4dzDHM49zA0tb2BAqwFWxxKpks7ARYCmoU156baXrI4hcllU4FLnFBTAH/8IZ89W3nfyZO3nEamr\nVOBS5xw8CH/9K0yaVHlf06bnV5KvyuHcw6TZ06r83InTJy4voEgdoQKXOikmBqZNq9kxnl33LN9m\nf0uLiBZVfvaeLvfUbDARC6jApd4yMJjafyr3XuehCcFFLFZlga9du5aUlBSio6Ox2WzMnDnT6edW\nrFjBqFGjOH36NKGhmg9ZLs/s2bBqlfm+uBiaNLE0johXcFngRUVFjB07lvT0dAIDAxkxYgTr1q1j\nyJAhF31u37597Nu3z6NBpX7bsQPGjIH4eHM7KsrSOCJeweV94Js3b6ZNmzYE/vSc8oABA1izZs1F\nnykqKmLhwoXMmjXLcynFJ7RrB927m6/YWKvTiNR9Ls/A7XY74eHhju2IiAjsP1vee/r06cyaNctR\n8lo2TarLbj9/q2BJiXuOebLwJGcqzEUxi84UueegInWUywKPiYmhoKDAsZ2Xl0dMTIxj++jRo5w6\ndYqVK1c6fvbiiy9y2223cf31lRd3nT17tuN9fHw88ef+e1l8zt690LPn+Uslfn7QrIaL2xzJPcLV\nL19NdENzEUybzcYT/Z6oYVKR2pWamkpqamq1PutyUeOioiLi4uJIS0sjKCiIESNGMH78eHr06EFA\nQMBFZ+cAfn5+l/wSU4say4W2bYNx48x/1kSFUeGYw+Tb7G+5b9V97Buv72Ok/rjiRY1DQ0N59dVX\neeyxx4iKiiIuLo7BgweTmJhIZGQkiYmJAGRnZ7N06VJsNhsLFizg0Ucfpfm5+T5FPOiR1Y/wzp53\n8PfzB6B/q/4WJxKpPS7PwN06kM7Afd727eengT1yBD77rOZn4Pd8cA+/6fIb7umqB3GkfrriM3AR\nd5o5E8rKoFUrc/uBBy7/GBmnMrju1escX1SWlZcx5hdj3JhSxHuowMVj8vJg0aLzd5p89x28/DLc\nXoMVyfJL82nVqBXbH9kOmGcnwQHBbkgr4n1U4OIx+/bBsmXw+9+b22PGQN++l3+cqf+Yyoq9KwA4\nU3GGto3bEhIY4sakIt5J18DFrV59Fc7dAZWdbS42vHlzzY756/d+zbCrh3Hb1bcBENEggrCgsJod\nVMRL6Bq41JqUFOjVC+LizO2OHd1z3MiQSJqH684mkQupwMXt4uPhV7+yOoVI/acClzrnxOkT/Db5\nt5ytML/9TD+ZzoM9HrQ4lUjdowKXGklPh4EDobzc3C4shJrOa3a84Dg/FPzAX+76CwA2bPRu0buG\nSUXqHxW4XLbSUnPdSoDDh6FTJ/jiC3Pbzw9+NsNCtfx565/5/NDnAOSV5BHRIIKBrQe6KbFI/aQC\nl8uWkAAbNkBQkLl9++3QqFHNjvn5oc/pEdODPi36ANCmcZsaphSp/1Tg4tTRo+bLmePH4aOP4Gfr\nelyW7KJspn05zXGd++vjX/Noz0e5s/OdV35QER+jAhen7r8fTpyAiIjK+4KDoXXrmh3/SO4R1n+/\nnjnxcwD4VYdfEd82vmYHFfExKnAfUlgIH38MFRVVf/aHH+CNN2DQoCsfL6c4h7kb5jrOsi+UVZhF\n09Cm3B93/5UPIOLjVOA+5F//gilT4Oabq/5sv35w9dU1G2//j/tJ3pfsdFGF9k3a89AvHqrZACI+\nTo/S13OrVsGePeb7gwfNx9s//9x9xz9VcoqXt7xMuVFead/R/KN8Y/+Gf4/5t/sGFPExrrrT5aLG\n4v2Sks5/GdmxIzz2mHuP/439G97Y+YbTfS0jWjJtwDT3DigiDrqE4gMeeeTKZgGsrtaNWjM7frbn\nBhARp3QGLiLipXQGXg888QRs2uR8X1oaBLj5d3nbD9vYl20uHLz/x/3uPbiIVJsKvB746iuYOBGu\nuabyvoAA+MUv3Dte4tpE/P38iQ2LBeC+6+5z7wAiUi0qcC9x993m4+vOnDplXuO+9tray/P0wKcZ\n0q4Gj2KKSI2pwL3EoUOQnAzdulXeFxAAjRvXfiYRsZYK3Is0aQJNm7rveN+f+p6j+ZeY8MSFUyWn\n3BdCRK6YCtxCZ87A3r1QneebiorcP/6olFGOqVsvR3iDcNo00myBIlZTgVto9Wpzpfb27av+7FVX\nQXS0e8cvryhn6bCl9G/V370HFpFaoQKvZZs3Q2am+X7TJnPtyPfeszaTiHgnFXgte/BBaNfu/DSt\nw4Z5drzCskLySvOc7isrL/Ps4CLiUSrwWmYYsHgxdO5cO+Pd/rfbSbOnEeQfVGmfv58/TUPd+K2o\niNQqFXg9VHK2xDF7WdGZItbct4a+LT04GYqIWKJaBb527VpSUlKIjo7GZrMxc+bMi/YnJSWRlZVF\nbGws27dv57nnnqNzbZ1iykU2ZGxg8PLBjjPuAL8AGgfrJnGR+qjKAi8qKmLs2LGkp6cTGBjIiBEj\nWLduHUMuWBCxsLCQRYsWAfD+++8zdepUVq9e7bnUckkFpQUM7TiUNfetsTqKiHhYlbMRbt68mTZt\n2hAYGAjAgAEDWLPm4nJ47rnnHO/Ly8sJDw93c0xxZdjfhhEyN4SQuSH8+r1fX/Z93SLinao8A7fb\n7RcVckREBHa73elny8rKePvtt3nllVfcl1CqdPz0cdaOWkvP2J4ATr+wFJH6p8oCj4mJoaCgwLGd\nl5dHTExMpc+VlZUxbtw45s2bR7t27dybUqoUHBBMSGCI1TFEpBZVWeD9+vUjIyODsrIygoKC2LRp\nE+PHjyc3N5eAgADCw8MpKipi/PjxTJ06lS5duvDhhx+SkJBQ6VizZ892vI+Pjyc+Pt6d/y4iIl4v\nNTWV1NTUan22Wosar127luTkZKKioggKCmLGjBkkJiZy1VVXMW3aNO6++27S09OJjTXnhy4qKmLL\nli0XD6RFjQHz/u/Vq2t2H/jO4zu5P+V+DMz/PQ/lHGLrI1vpHtPdTSlFpK5w1Z1ald7DzpyBb789\nP2HVnXfCP/5RswL/+NuPWbxlMX+67U+Aeatgp6s6YbPZ3JBYROoSV92pB3k8LCUFHn0UWrc2t6Oi\nzImpqrL7xG4mfDbB6W/cj8U/cm3Ta+ka3dXNaUXEm6jAPWDrVjj60zTbmzbBHXfAihWXd4yDOQfx\ns/kx9+a5Tve3a6wvikV8nQrcA373O3PCqrAwc3v48Or9uhnrZrAraxcAxwuO065JOwa2HuihlCLi\n7VTgbnDsGLz77vntnBz47DPo2PHyjrPq21WM7TXWsVjCtVG1uMiliHgdFbgbrF0Ly5ebl0oAxo6F\nli2v7FiD2w7WtW0RqRYVuAtJSeZiwlXZvx+uvx4WLPB8JhGRc1TgLrz8MkycaC4m7EqvXtBXs7WK\nSC1TgVfhd7+78sshIiKepAKvZX/b+zf2/7jf6T57ofNJwkREnKlyOllf8vnn0KbN+deJE9CggXvH\nmPPVHLKLsp3um9R3Eu2bVGOJehERdAZOeTlUVJjvDx+G/v1h/nxzu0ED88nJmrAX2ln5zUrHE5W5\nxblM7DORzk21YpGI1IzPF3jXrnDwIJybRmTaNPPs212+PPQlf972Z4Z2GArAA3EP0KpRK/cNICI+\ny+cL/NQp87H3Zs08N0av5r146baXPDeAiPgkXQMXEfFSPn8G7gk7j+/kh/wfHO9FRDxBBe4BI5NH\n0iqiFQ2DGgJw9zV3W5xIROojFbgHGIbB63e+TsfIy5zNSkTkMugauIiIl/K5M/C0NHj11fPb+fnV\n+3Vnys9c8gnKnystL72CZCIil8fnCvyrr2D3bhg50tx++WWIjq761638ZiUTP5tIi4gWVX62aWhT\nIkMia5hURMQ1nytwgOuugwkTqv5cfmk+p8tOA+YTlXdfezfL7lrm4XQiItXjkwVeXT1f60lBWQH+\nNn8AxvUeZ3EiEZHzfL7Az5Sf4WzFWaf7is8Ws+PRHbSM0HyyIlL3+HyBd32lK5l5mdjOTYZygeCA\nYBoGNrQglYhI1XyiwDt1ggMHzm9Pn37+fX5pPt9P+p5mYR6cDEVExAN84j7wnBw4eRIMAz5IS+aF\nkBBC5pqv3JJcggOCrY4oInLZfOIM/EL2Qjujuo/ipaHm7ID+fv4E+QdZnEpE5PL5XIEDBPgFEBIY\nYnUMEZEa8YlLKCIi9VGVZ+Br164lJSWF6OhobDYbM2fOvGh/SUkJTz75JC1btuTAgQM89dRTXH31\n1R4LfCkffgjHj5vv08o+5aOiqY59OffCwL+BfwDkFOfwmy6/qfV8IiLuZjPOLdboRFFREXFxcaSn\npxMYGMiIESMYN24cQ4YMcXxm/vz5BAQE8OSTT/LNN98wbtw4NmzYUHkgmw0XQ9VYkyaQkADBwbA7\nZAlZAVu4vugZAEKC4fFJEGA+j0PrRq0JbxB+yWOlpqYSHx/vsazupKye4S1ZvSUnKOuVctWdLi+h\nbN68mTZt2hAYGAjAgAEDWLNmzUWf+fTTT7nhhhsA6NatG7t37+b06dPuyO3S3r0wefL5V/4vf8c3\nfQayq+c6NDSIAAAGTklEQVRAfmi5mDtuiubvL3Xl7y91ZVlSV+Jiu9I12ny5Km8wf/O8hbJ6hrdk\n9ZacoKye4PISit1uJzz8fNlFRERgt9ur9ZmwsDA3R73Y2rWwcycMH25uBzf5lBnxK2gUHAHANU2v\n8ej4IiJWc1ngMTExFBQUOLbz8vKIiYm56DPR0dHkXzAna35+PtGXmN4vZvKdNcl6kaIiiL4R1v30\n/E3Z6dMMbD2ARsGN3DaGiEidZrhQWFhodOzY0SgtLTUMwzASEhKMdevWGTk5OUZ+fr5hGIYxf/58\nY8GCBYZhGMaePXuMQYMGOT1Whw4dDEAvvfTSS6/LeMXFxV2yo11+iQnmXSjJyclERUURFBTEjBkz\nSExMJDIyksTERMddKLGxsRw8eJDp06fTsaOWEhMR8bQqC1xEROomPcgjIuKlVOAiIl5KBS4i4qVU\n4CIiXkoFLiLipTw+nWxVk2HVJYcOHWLGjBn07NmTo0ePctVVVzFjxgyrY11ScXExffv25dZbb2Xh\nwoVWx7mk7777juXLl9O4cWP+7//+j5kzZ9K7d2+rYzm1cOFCMjIyiIqK4vvvv+f11193TCVhtRMn\nTvDss8+yZ88etm7dCkBOTg5PP/007du358CBA8ybN++SD9LVJmdZJ0+eTMOGDQkLC2P37t0sXry4\n0oOBVnCW9Zy5c+eyePFiTp48aVG6Krh6kKemzj0IVFZWZhiG+SDQP//5T08OWSPbtm0zVq9e7dju\n0qWLsWPHDgsTufbEE08YDzzwgPHkk09aHeWSzp49a9x4442OPwMnTpwwTp48aXEq544fP25ERkY6\ntu+66y7jvffeszDRxZKTk41PPvnE6NWrl+Nn//3f/2188MEHhmEYxieffGKMGjXKqngXcZb12Wef\ndbxPSkoyJk6caEW0SpxlNQzDWL9+vTFlyhQjKirKomRV8+gllOpMhlWX9OrVizvvPP+4f0VFhcfn\ndLlS7777LgMHDqRdu3ZWR3Fp27ZtGIbB0qVLWbBgASkpKTRt2tTqWE6FhoYSFBREXl4eANnZ2TRr\nVnfWSk1ISKj05/HCyeT69+9fZ/5+Ocv6/PPPO96Xl5dfNIeSlZxlzcrK4v3332fixIkenUW1pjx6\nCaU6k2HVVSkpKQwdOpROnTpZHaWS9PR09u3bx9y5c9m9e7fVcVzKyMhg+/btJCcnExMTw+jRo2nQ\noAEPPvig1dEqiYiIYOHChYwcOZLmzZvTs2dPBg0aZHUsly78OxYREUFubi4VFRX4+dXdr7dOnTrF\nl19+yapVq6yO4lRFRQXTp0/nhRdeIDc31+o4Lnn0d7k6k2HVRevXr+err77ixRdftDqKUx999BHB\nwcEkJSWxceNGtm7dypIlS6yO5VRERAStW7d2/L7379/f6XzxdcGuXbtYuHAhn376KcuWLSMoKIik\npCSrY7kUHR3t+DuWn59PkyZN6nR55+XlMWHCBP7yl7/QuHFjq+M49fXXXxMYGMjSpUtZunQpxcXF\nLFiwgIMHD1odrRKPnoH369ePjIwMysrKCAoKYtOmTYwfP96TQ9bYmjVr+Ne//sXixYs5duwYmZmZ\n9OvXz+pYF3nmmWcc74uLiyksLOSxxx6zMNGl9e3bl4KCAkpLS2nQoAEZGRl07tzZ6lhOHTt2jMjI\nSEcBNmvWjCNHjlicyrU77riDTZs2cc8997Bx40aGDRtmdaRLys7OZtKkSSxYsIDmzZvz4YcfkpCQ\nYHWsSq6//nquv/56AL7//nveeustpk2bZnEq5/xnz54921MHDwwMpEuXLrzwwgts2bKFFi1a8MAD\nD3hquBrbsWMHw4YNwzAM/vrXv/LOO+/QoUMH4uLirI7m1KpVq0hOTiYzM5OQkBC6detmdaRKQkJC\n6NChAy+//DK7d+8mKyuLOXPmEBBQ99bT7tChAzt27GDdunVs3bqVHTt28Ic//IGIiAirowGwYcMG\n3nnnHXbv3k1JSQm9e/fmpptu4rXXXmPPnj1s2rSJhQsX0rBhQ6ujOs06ePBgsrKy+PTTT1m+fDnf\nfvst9957r9VRnWYNCAjg0KFDvPLKK2zZsgUwi72u3JF0jiazEhHxUnX3YpmIiLikAhcR8VIqcBER\nL6UCFxHxUipwEREvpQIXEfFSKnARES+lAhcR8VL/D02xEVmWZuiZAAAAAElFTkSuQmCC\n",
	"text": "<matplotlib.figure.Figure at 0x107dd0a50>",
	"output_type": "display_data",
	"metadata": {}
	}
	],
	"language": "python",
	"trusted": true,
	"collapsed": false
	}
	],
	"metadata": {}
	}
	],
	"metadata": {
	"name": "",
	"signature": "sha256:0795d55207b5da4a51d307bb6088633cdbc813297e7a9f5ee8bba7ecb1765be8",
	"gist_id": "29457bec83da6865f5c2"
	},
	"nbformat": 3
	}