{
"metadata": {
"name": "",
"signature": "sha256:67b73f8366759bb2837043e5a461a7ce4ff077967115f36e8017eaa1318f1c0f"
},
"nbformat": 3,
"nbformat_minor": 0,
"worksheets": [
{
"cells": [
{
"cell_type": "code",
"collapsed": false,
"input": [
"from __future__ import division, print_function\n",
"%matplotlib inline\n",
"import sys\n",
"sys.path.insert(0,'..') # allow us to format the book\n",
"sys.path.insert(0,'../code') \n",
"# use same formatting as rest of book so that the plots are\n",
"# consistant with that look and feel.\n",
"import book_format\n",
"book_format.load_style(directory='..')"
],
"language": "python",
"metadata": {},
"outputs": [
{
"html": [
"\n",
"\n"
],
"metadata": {},
"output_type": "pyout",
"prompt_number": 3,
"text": [
""
]
}
],
"prompt_number": 3
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"import numpy as np\n",
"import matplotlib.pyplot as plt\n",
"from gif_animate import animate\n",
"from stats import gaussian\n",
"\n",
"def plt_g (mu, variance):\n",
" xs = np.arange(2,8,0.05)\n",
" ys = [gaussian (x, mu, variance) for x in xs]\n",
" plt.plot (xs, ys)\n",
" plt.ylim((0,1))\n",
"\n",
"\n",
"mu = 2\n",
"sigma = 0.6\n",
"def ganimate(frame):\n",
" global mu, sigma\n",
" if frame < 25:\n",
" mu += .2\n",
" elif frame == 25:\n",
" mu = 5\n",
" elif frame < 37:\n",
" sigma -= 0.05\n",
" else:\n",
" sigma += 0.05\n",
" \n",
" plt.cla()\n",
" plt_g(mu,sigma)\n",
"\n",
"animate('04_gaussian_animate.gif', ganimate, frames=80, interval=50)"
],
"language": "python",
"metadata": {},
"outputs": [
{
"metadata": {},
"output_type": "display_data",
"png": "iVBORw0KGgoAAAANSUhEUgAAAcEAAAGrCAYAAAC41CzUAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XlcVXX+P/DXvey7IPsiiwKKIsgmAgoqUq5R7rZpzpQ2\nNZnTbxz61jSVU9M0X2da1DZHs2XUrNQUTVSQEFQQEWRTQRFlURHZZJF7z++PJr5duSwqcoDP6/l4\n+Hjo+XzO5X1Pt/vinPM5n49CkiQJREREAlLKXQAREZFcGIJERCQshiAREQmLIUhERMJiCBIRkbAY\ngkREJCyGIBERCavLEExOTsasWbPg7OwMpVKJzz//vMsXzcnJQWRkJIyNjeHs7Iw333yzR4olIiLq\nSV2GYENDA0aPHo333nsPRkZGUCgUnfavra3FlClT4ODggIyMDLz33nt49913sWbNmh4rmoiIqCco\n7mTGGDMzM6xduxZPPPFEh33Wr1+PuLg4VFZWwsDAAADw17/+FevXr8elS5fuvWIiIqIe0uP3BNPS\n0jB+/Pi2AASAmJgYlJWVoaSkpKd/HBER0V3T7ekXrKiowJAhQzS22dnZtbW5urq2ba+pqenpH09E\nRNTGwsKi0/YePxPs6p4hERFRX9HjIWhvb4+KigqNbZWVlW1tREREfUWPXw4dN24cVq1ahebm5rb7\nggkJCXByctK4FHq7rk5ZRZCRkQEACAoKkrkS+fFYaOLx0MTjoYnHQ9Od3Grr1iMSWVlZyMrKglqt\nRklJCbKyslBaWgoAiIuLQ3R0dFv/RYsWwdjYGIsXL0Zubi6+++47vPPOO1i5cuVdvBUiIqL7p8sQ\nTE9PR0BAAAICAtDU1ITXXnsNAQEBeO211wD8PNiluLi4rb+5uTkSEhJQVlaGoKAgPP/883jppZfw\n4osv3r93QUREdBe6vBwaFRUFtVrdYfvGjRvbbRs1ahQOHz58b5URERHdZ5w7lIiIhMUQJCIiYTEE\niYhIWAxBIiISFkOQiIiExRAkIiJhMQSJiEhYDEEiIhIWQ5CIiITFECQiImExBImISFgMQSIiEhZD\nkIiIhMUQJCIiYTEEiYhIWAxBIiISFkOQiIiExRAkIiJhMQSJiEhYDEEiIhIWQ5CIiITFECQiImEx\nBImISFgMQSIiEhZDkIiIhMUQJCIiYTEEiYhIWAxBIiISFkOQiIiExRAkIiJhMQSJiEhYDEEiIhIW\nQ5CIiITFECQiImExBImISFgMQSIiEhZDkIiIhMUQJCIiYTEEiYhIWAxBIiISFkOQiIiExRAkIiJh\nMQSJiEhYDEEiIhIWQ5CIiITFECQiImExBImISFgMQSIiEhZDkIiIhMUQJCIiYTEEiYhIWAxBIiIS\nFkOQiIiExRAkIiJhMQSJiEhYDEEiIhIWQ5CIiITFECQiImExBImISFgMQSIiEhZDkIiIhMUQJCIi\nYTEEiYhIWAxBIiISFkOQiIiExRAkIiJhMQSJiEhY3QrBdevWwd3dHUZGRggKCkJKSkqn/ePj4xEa\nGgpzc3PY2NggNjYWZ8+e7ZGCiYiIekqXIbh161asWLECr7zyCrKyshAWFoapU6eitLRUa/9z584h\nNjYWUVFRyMrKwoEDB9DU1IRp06b1ePFERET3ossQXLNmDZYsWYKlS5fC29sb77//PhwcHLB+/Xqt\n/bOysqBWq/H222/Dw8MDfn5+WLVqFYqKinD9+vUefwNERER3q9MQbGlpQWZmJmJiYjS2x8TEIDU1\nVes+4eHhMDU1xaeffgqVSoW6ujps2rQJISEhsLKy6rnKiYiI7pFCkiSpo8aysjI4OzsjOTkZERER\nbdvfeOMNfP311ygoKNC6X2pqKmJjY1FdXQ21Wo0xY8Zg7969sLGx0ehXU1PT9nfeMyQiop7g6enZ\n9ncLC4tO+/b46NDi4mLExsZiyZIlyMjIQFJSEszMzDBv3jx0krdERES9TrezRmtra+jo6KCyslJj\ne2VlJRwcHLTu8/HHH8PFxQXvvPNO27Yvv/wSLi4uSEtLQ1hYmNb9goKC7rT2AScjIwMAjwXAY3E7\nHg9NPB6aeDw0/foqY1c6PRPU19dHYGAg9u/fr7E9ISGhwzCTJAlKpebL/vJvtVrd7cKIiIjuty4v\nh65cuRKbNm3Chg0bkJ+fjxdeeAEVFRVYtmwZACAuLg7R0dFt/WfNmoXMzEy8+eabOHv2LDIzM7Fk\nyRIMGTIEgYGB9++dEBER3aFOL4cCwLx581BVVYXVq1ejvLwcvr6+iI+Ph4uLCwCgoqICxcXFbf0j\nIiKwdetW/O1vf8Pf//53GBsbY9y4cdi3bx+MjIzu3zshIiK6Q12GIAAsX74cy5cv19q2cePGdtvm\nzJmDOXPm3FtlRERE9xnnDiUiImExBImISFgMQSIiEhZDkIiIhMUQJCIiYTEEiYhIWAxBIiISFkOQ\niIiExRAkIiJhMQSJiEhYDEEiIhIWQ5CIiITFECQiImExBImISFgMQSIiEhZDkIiIhMUQJCIiYTEE\niYhIWAxBIiISFkOQiIiExRAkIiJhMQSJiEhYDEEiIhIWQ5CIiITFECQiImExBImISFgMQSIiEhZD\nkIiIhMUQJCIiYTEEiYhIWAxBIiISFkOQiIiExRAkIiJhMQSJiEhYDEEiIhIWQ5CIiITFECQiImEx\nBImISFgMQSIiEhZDkIiIhMUQJCIiYTEEiYhIWAxBIiISFkOQiIiExRAkIiJhMQSJiEhYDEEiIhIW\nQ5CIiITFECQiImExBImISFgMQSIiEhZDkIiIhMUQJCIiYTEEiYhIWAxBIiISFkOQiIiExRAkIiJh\nMQSJiEhYDEEiIhIWQ5CIiITFECQiImExBImISFgMQSIiEhZDkIiIhMUQJCIiYTEEiYhIWAxBIiIS\nVrdCcN26dXB3d4eRkRGCgoKQkpLS5T7/+te/MHz4cBgaGsLR0RFxcXH3XCwREVFP0u2qw9atW7Fi\nxQqsX78eERERWLt2LaZOnYq8vDy4uLho3WflypXYs2cP/vGPf8DX1xc1NTUoLy/v8eKJiIjuRZch\nuGbNGixZsgRLly4FALz//vvYt28f1q9fj7feeqtd/8LCQnz44YfIycmBt7d323Y/P78eLJuIiOje\ndXo5tKWlBZmZmYiJidHYHhMTg9TUVK377Ny5Ex4eHoiPj4eHhwfc3d2xePFiXL16teeqJiIi6gGd\nngleu3YNKpUKdnZ2GtttbW1RUVGhdZ/i4mKUlJRg27Zt2Lx5MwDgpZdewsyZM5GWlgaFQqF1v4yM\njLupf0Disfg/PBaaeDw08Xho4vH4maenZ7f7dnk59E6p1Wo0Nzfjiy++wLBhwwAAX3zxBby9vZGR\nkYHg4OCe/pFERER3pdMQtLa2ho6ODiorKzW2V1ZWwsHBQes+Dg4O0NXVbQtAABg2bBh0dHRw8eLF\nDkMwKCjoTmsfcH75LY7Hgsfidjwemng8NPF4aKqpqel2307vCerr6yMwMBD79+/X2J6QkICwsDCt\n+0RERKC1tRXFxcVt24qLi6FSqeDq6trtwoiIiO63Lp8TXLlyJTZt2oQNGzYgPz8fL7zwAioqKrBs\n2TIAQFxcHKKjo9v6R0dHIyAgAE899RSysrJw8uRJPPXUUwgNDeVvKURE1Kd0eU9w3rx5qKqqwurV\nq1FeXg5fX1/Ex8e3PSNYUVGhcdanUCiwe/du/P73v8eECRNgZGSEmJgYrFmz5v69CyIiorvQrYEx\ny5cvx/Lly7W2bdy4sd02e3t7bNu27d4qIyIius84dygREQmLIUhERMJiCBIRkbAYgkREJCyGIBER\nCYshSEREwmIIEhGRsBiCREQkLIYgEREJiyFIRETCYggSEZGwGIJERCQshiAREQmLIUhERMJiCBIR\nkbAYgkREJCyGIBERCYshSEREwmIIEhGRsBiCREQkLIYgEREJiyFIRETCYggSEZGwGIJERCQshiAR\nEQmLIUhERMJiCBIRkbAYgkREJCyGIBERCYshSEREwmIIEhGRsBiCREQkLIYgEREJiyFIRETCYggS\nEZGwGIJERCQshiAREQmLIUhERMJiCBIRkbAYgkREJCyGIBERCYshSEREwmIIEhGRsBiCREQkLIYg\nEREJiyFIRETCYggSEZGwGIJERCQshiAREQmLIUhERMJiCBIRkbAYgkREJCyGIBERCYshSEREwmII\nEhGRsBiCREQkLIYgEREJiyFIRETCYggSEZGwGIJERCQshiAREQmLIUhERMJiCBIRkbAYgkREJCyG\nIBERCYshSEREwupWCK5btw7u7u4wMjJCUFAQUlJSuvXiZ8+ehZmZGczMzO6pSCIiovuhyxDcunUr\nVqxYgVdeeQVZWVkICwvD1KlTUVpa2ul+LS0tWLBgASIjI6FQKHqsYCIiop7SZQiuWbMGS5YswdKl\nS+Ht7Y33338fDg4OWL9+faf7rVq1Cv7+/pg7dy4kSeqxgomIiHpKpyHY0tKCzMxMxMTEaGyPiYlB\nampqh/vt2bMHe/bswQcffMAAJCKiPku3s8Zr165BpVLBzs5OY7utrS0qKiq07lNWVoann34aO3bs\ngLGxcbcLycjI6HbfgY7H4v/wWGji8dDE46GJx+Nnnp6e3e7b46NDH3/8cSxfvhzBwcE9/dJEREQ9\nqtMzQWtra+jo6KCyslJje2VlJRwcHLTuk5iYiOTkZLz++usAAEmSoFaroaenh/Xr1+M3v/mN1v2C\ngoLupv4B5Zff4ngseCxux+OhicdDE4+Hppqamm737TQE9fX1ERgYiP3792P27Nlt2xMSEjB37lyt\n+5w+fVrj3zt27MBf//pXpKenw9HRsduFERER3W+dhiAArFy5Eo8//jhCQkIQFhaGjz76CBUVFVi2\nbBkAIC4uDunp6Thw4AAAwMfHR2P/48ePQ6lUtttOREQkty5DcN68eaiqqsLq1atRXl4OX19fxMfH\nw8XFBQBQUVGB4uLiTl+DzwkSEVFf1GUIAsDy5cuxfPlyrW0bN27sdN/Fixdj8eLFd1wYERHR/ca5\nQ4mISFgMQSIiEhZDkIiIhMUQJCIiYTEEiYhIWAxBIiISFkOQiIiExRAkIiJhMQSJiEhYDEEiIhIW\nQ5CIiITFECQiImExBImISFgMQSIiEhZDkIiIhMUQJCIiYTEEiYhIWAxBIiISFkOQiIiExRAkIiJh\nMQSJiEhYDEEiIhIWQ5CIiITFECQiImExBImISFgMQSIiEhZDkIiIhMUQJCIiYTEEiYhIWAxBIiIS\nFkOQiIiExRAkIiJhMQSJiEhYDEEiIhIWQ5CIiITFECQiImExBImISFgMQSIiEhZDkIiIhMUQJCIi\nYTEEiYhIWAxBIiISFkOQiIiEpSt3AUSikiQJ9Y01qGm4jrqbNai7eQP1jbVoarmJ5pZGNN9qxK3W\nW1CrVbhWdQ0SJORcSYRSqQMdpQ709QxhqGcEA30jGBmYwMx4EMyMLWBmPAiWptbQ1zOQ+y0S9XkM\nQaL7rO5mDcqrLqKy+hKuVF9GZfVlXK+9guq6q7jV2nJHr3Wxqvt9TYzMYWVmg8EWdrCzdIadpRPs\nrJxhb+UCPV39O3wXRAMTQ5CoBzU01eFCeSEuVJzBpSvFuHS1GDUN1+WppbEWDY21KL1SpLFdqdSB\nvaUznGzcMcRuGNzsveFs4w4dHX4dkHj4qSe6BzX113H2Ug7OXjqNosu5uHKjTO6SuqRWq1BWVYKy\nqhKkFyQBAPR09DHEbhiGOY+Ep7Mv3By8oa/Ly6k08DEEie7ArdZbKLqci7ySTORfyERl9aV7ej0D\nfSNYmlq33c8zNbKAkYExDP57r09f1wBKpQ5KLpQAANzd3aFSt0KlVqHlVhOaWhrR3NKIm811P99X\nbKxBbUM1btRXQa1Wdf99qVpQVJaHorI8/Hj8G+jq6MHDcQR83ALg4xYIO0tnKBSKe3qvRH0RQ5Co\nCzeb6pF7IQPZ544iv+QkWlqb72h/HR1dOFgNgf1gF9hZOsHW0hnWFvYYbG4LIwOTboWLoj4DABDo\nHdStn6lWq1DTcB3Xa6/gSnUZKqsvo7L6EsqrLuJ67ZUu929V3cKZ0mycKc3Gjp82wcrcFqM9xmL0\nsFB4OAyHUqnTrTqI+jqGIJEWN5vrkVN0DCfOpOBMaXa3z6oUCiUcrV3hbu8NV3tPONsMhb2Vc6/f\nb1MqdWBpZgNLMxsMdRqp0XazqR6Xrp5H6ZUiXKgoxIXywi7vW16vvYKkrB+QlPUDTI0s4DdsHAK9\nx8PDcQSUCj5pRf0XQ5Dov1pVt5B7/gTSCxKRe+EEVKrWLvdRKpQYYucJT+dR8HT2hbuDNwz0jXqh\n2rtnbGgKLxdfeLn4Avj5UY3quqs4dzkXZ0tzcPZSDq7XXe1w//rGGhzJ2YcjOftgYToYgV4RCBkx\nEY7Wbr30Doh6DkOQhHfpajGO5h7EicJkNDTVddnfwsQKI9wC4OMaAO8hfjAyMOmFKu8fhUIBK3Nb\nhJjbImTERADAleoy5JdkIu9CJs5eykGr6pbWfWvqq3AocycOZe6Es40HxvpMQtDwSJgYmvXmWyC6\nawxBElJzSyMyz6TgyOn9uFh5tsv+dlbO8BsaCl+PsRhiN2zADxKxtXSEraUjIv1noLmlEfklJ5Fd\ndAy559PR2HJT6z6Xrhbj0uFi7Ez5HP6eYQgfFQMPR58Bf6yof2MIklCuVJfhp+x4HMs7hKYOvsx/\nYWPhgADvCAR4jYfD4CG9VGHfY6BvBH/PMPh7hqFVdQsFJVnIPJOC7OJjaLnV1K5/q+oWMgoOI6Pg\nMBwGD8H40dMQPCIKBnqGMlRP1DmGIA14kiSh4GIWkk7+gPySzE77GhuaIdBrPEJGTBTijO9O6ero\nYZRHMEZ5BKPlVjNyio/hWH4iCkuyIEFq17+86iK2JX6EH45sRujIaET6z4CVua0MlRNpxxCkAetW\n6y1knklGYuYulFWVdNhPAQW8Xf0xbuQUjHIPhp6uXi9W2X/p6xkg0HsCAr0n4EZ9FY7nHUJa7gFU\n1Va269vYchOJJ3fhcNZu+HuGYVJALIbYDZOhaiJNDEEacJpaGpF6+kckZu7qdOi/mZEFQkdGI2xU\nDAZb2PVihQPPINPBiAmZi+jg2Si8eApHcn7E6eLjUEtqjX5qSY3MMynIPJMCT2dfTAmaDe8hfjzj\nJtkwBGnAaGiqw+GTu5F8ag9uNtd32M/VzhMT/KfDf1g4z/p6mFKhxAjXMRjhOgbVdVeRkr0PqbkJ\naGisbdf35+nmcuBiOxQxwXPgO3QsnzmkXscQpH6vvrEWiZk7kXxqD5q1DNQAfn6I3W9YKCaOeQju\nDt69XKGYLM1sMDP8cTw4dj4yCpORmLkTFddL2/UrvVKEDXvegeNgVzwwdh78ho1jGFKvYQhSv9XQ\nVIeDJ3Yg+dQeraMUAUBf1wDjRk1BlP9MXvKUiZ6uPsaNjEaoz2Tkl2Ti4IkdOHspp12/sqoSbIx/\nF/ZWLpgauoBhSL2CIUj9TmNzAxJP7kLSyR86fMzB2MAUE/ymY4L/dJgamfdyhaSNQqGAj1sgfNwC\ncaHiDA5kfIvsomPt+lVcL8XG+HfhZO2GaeMWYZR7MO8Z0n3DEKR+o6W1GT+d2ouEjG9xs4OZXUyN\nLDAp4CFEjJ4Kwz4+fZnI3Oy98JsZcSivuoj9x79B5pmUdo9YXL52AZ/+8BZc7b0wM+zxtmneiHoS\nQ5D6PJVahWN5h7D32BbU1GtfWt3MyAKTgx5GuO+DfCi7H3EYPARPTv0DHhw7Hz+mf4MThT9Bum1E\naUnFGXz43asY7joGM8Meg4vtUJmqpYGIIUh9liRJyCk+jl1HNqPyuvZ1+0wMzRAd9AjGj54GfT0u\nAttf2Vk544kHXsQDwXOx99hWZJ75qV2fgpKTKCg5iSDvSAwx9YWp4SAZKqWBhiFIfdK1ujKcuHAA\nlbUXtbYb6htjUsBDiBozi5c9BxA7K2csnvoHxATPxp60r5FTfLxdn4zCw8hUpGC4QzB8fIfD2MBU\nhkppoGAIUp9SXXcNP6R+gYyCw1rbdXX0MMFvGqYEzYYJB7wMWI7WbvjtzJdxvrwQPxzZjHOXczXa\n1ZIKeWVH8eamXEwNXYDwUQ/0+pqNNDDwU0N9QvOtJhzI+A6HMnfgVmtLu3YFFAgZMRFTQxfCytxG\nhgpJDu4O3nh+9mrkl2RiV8rmdtPfNTTVYXvSp/jp1F7Ejl+Mke5BMlVK/VW3H8JZt24d3N3dYWRk\nhKCgIKSkpHTYNykpCQ899BAcHR1hYmICPz8/bNy4sUcKpoFFkiSkFxzG6s2/w4/Ht2kNwBGuAfjj\non/i0ZjfMwAF9MujFX9ctAaPTnkeFqaD2/WprL6Ej3etxvodb3R4/5hIm26dCW7duhUrVqzA+vXr\nERERgbVr12Lq1KnIy8uDi4tLu/5paWnw8/PDn/70Jzg4OGDfvn14+umnYWhoiIULF/b4m6D+6WLl\nOWw//CkulBdqbR9kbItFMc9iuKt/L1dGfZFSqYOxPpMxxjMCX+5ej9OXUtGq1vylKb8kE29/dQoT\nRk/Dg6Hzeb+QuqSQJKn9+ie3GTt2LPz9/fHxxx+3bfPy8sKcOXPw1ltvdesHzZ8/HyqVCtu3b2/b\nVlNT0/Z3CwuLO6l7QMrIyAAABAUN7Es6dTdr8EPqFziWe1Dr8jtmRhYY6RiOYXb+CAkOkaHCvkeU\nz0Z3ZWRkoLGlHqU3T3f6OZoR/jjG+kwa8DPP8POh6U6ypctPRktLCzIzMxETE6OxPSYmBqmpqXdU\nlJWVVbf708CjUquQfGoPVm9+FkdzD7T74tLR0UV04CN45cn18LIPGPBfXHRvjPRNsSj6Ofy/Rf+L\noU4j27XXNdbgPwc+xD+3rkJJxVkZKqT+oMszwbKyMjg7OyM5ORkRERFt29944w18/fXXKCgo6PKH\n7N69G4888ghSU1M1flP5dVqfPcsP6UB2pbYUx4r2ovrmFa3tLlZeCHSLhrkRf1GiOydJEi5WFSDj\nwgE0NNdo7eNpNwZjXCfCUM+4l6uj3ubp6dn2967OBO/76NAjR47g0UcfxQcffMBTdQE1tjQgs+Qg\niq5ka203NxqMEPcYOFpyFhC6ewqFAq7WI+BkOQy5l9Nw+nIqVOpWjT5nK0+ipKoAAa4T4Wk3hvOR\nEoBuhKC1tTV0dHRQWam5WnRlZSUcHBw63TclJQXTp0/Hm2++iWeeeabTvgzIgXVdX61W4cjp/did\n8SUamxvatRvoGeLBsQsQ6T8dujrt1/QbSMeiJ/B4aOrseIRiHKpqH8WO5I04VXRUo62ltRFHi+JR\n3nAO8yYug4utR6/Ue7/x86Hp11cZu9LlTRd9fX0EBgZi//79GtsTEhIQFhbW4X7JycmYNm0aXn/9\ndfz+97/vdkHU/12sPIc1W1fhm8SPtQZgkHckXnlyHSYHxmoNQKJ7NdjcDktn/AnPxv4FtoMc27WX\nVJzBP7a8hO1Jn2r9jJI4unU5dOXKlXj88ccREhKCsLAwfPTRR6ioqMCyZcsAAHFxcUhPT8eBAwcA\n/Pyc4PTp0/Hcc89h4cKFqKioAADo6OjAxobPeQ1Ujc0N2J36FVKy92odrecweAjmRD0NT+dRMlRH\nIhru6o9Vj76HpJO78OPxbWhpbW5rkyQ1kk/tQdbZVDwSuRRjPMN5iVRA3QrBefPmoaqqCqtXr0Z5\neTl8fX0RHx/f9oxgRUUFiouL2/p//vnnaGpqwrvvvot33323bbubm5tGPxoYJElC5pkUfJ/8b9Te\nrG7XbqBniKmhCxDpN4NTW1Gv09PVw5Tg2Qj0noDvkjcg+7ZLpLU3q7Fp7z+QlpuAuVHPwNay/Zkj\nDVzdek7wfuFzgpr643X9qzfKsS3xIxRePKW13d8zDA+PfwqWZtZ39Lr98VjcTzwemu7leOSez8D2\npE9RVVvZrk1XRw8xwXMwOfAR6On2n0v1/HxoupNs4a/ldFdaVbdw8MQO7D/+DW6p2k91Zm1hjzlR\nT8PHLUCG6og6NtI9CJ4uvkhI344DGd9rjCJtVd1C/NH/IKMwGfMnLeelewEwBOmOFV3Ow5ZD67TO\n0fjLA+9TgmdDX5fr+1HfpK9rgOnjHkWQdyS2JX6Ms5dyNNqvVF/GB9++grE+kxEb8SRXLBnAGILU\nbQ1NddiVshlpuQla2z2dfTFv0jLYWTr1cmVEd8fOyhnPPfIGMgoP4/vkjahv1BxafyzvIE6fT8fD\n45cgeHgUB84MQAxB6pIkSThRmIzvk/+Nusb2z9+YGJnzS4L6LYVCgeDhUfBxC8QPRzYj9bTmL3kN\njbX4cv97OJ53CPMmLefAmQGGIUidulZTgW2HPkLBxSyt7aE+k/EQLxfRAGBiaIYFk3+HkBGTsOXg\nOlRcL9VoP3MpB3/76gU8EDIXkwMf5jOuAwRDkLRSqVpx6OQu7Du2Resaf3aWzpg/eTmGaZm4mKg/\n83AcgT8uWoNDJ3bgx9sGfrWqbmFP2tc4UfgTFkx+Fh6OI2SslHoCQ5DauVBxBlsOrkPZtQvt2nR0\ndBETPBfR/WwIOdGd0NXRQ0zIXIzxisC2Qx+hsFTzEaCK66X41zdxCBsVg1nhT8DYkOsW9lcMQWrT\n2HwTu1O/7HDGl2HOozB/0nIOfCFh2AxywLMP/wUZ/70nfvvAmdTT+5FTfByzI3/DGWf6KYYgQZIk\nZBcdxfakT1HTcL1du7GhGR4evxghIybxf3ISzs8DZyLh4zoGO1M+x9G8gxrtdTdvYNPef+B43iHM\nnfQMBpvbyVQp3Q2GoOCq665ie9KnyCk+rrU9eHgUYscvgZkxZ/QhsZkYmWPRlOcRPCIKWw+ux5Ub\nZRrteSWZePuL32Nq6EJEjZkJHaWOTJXSnWAICkqtViH5VDz2pH2F5ltN7dqtLewxf9JyeA/xk6E6\nor7L09kXqx79FxLSv0VCxrcaM860tDZjZ8omZBQexoJJz8LV3rOTV6K+gCEooNIrRdhycB1KrxS1\na1MqdRAd+DBiQuZyxheiDujp6mPauIUI8I7A1oPrUVSWp9F++ep5rNn6R4z3m4bp4x6FkQFXs++r\nGIICaW7NcK9GAAAWnElEQVRpxJ6j/8HhrN2QJHW7djd7byyYvByO1m69XxxRP2Rv5YLn56zG0dyD\n2JmySWNtQgkSkk/twamio5gT+RuMHhrKe+p9EENQENlFx7A96RPcqK9q12aob4yZ4Y8j3PcBKBVd\nrrNMRL+iVCgRNmoKRrkH4/uf/o0Thcka7TX1Vdiw5x2Mcg/GnKinYWXONVX7EobgAHe99iq2H/4U\npzsY+OLvGYbZkb+BhYlVL1dGNLCYmwzCkw+uRMiIidh26KN2SzWdPp+OM6XZPw+c8efamn0F/ysM\nUCq1CoezfkD80S1o0TLwxcrMBnMnPoOR7lx/jKgnjXAdg7jH3sePx7fhYOYOqNWqtrZfBs6kFyRh\n/qRlcHcYLmOlBDAEB6TisgJsS/xI64wvSoUSEwMewoNj58NAz7D3iyMSgL6eAWaGP45A7wnYemg9\nzpcXaLSXXbuAf22Lw7hRUzAz/HGYGJrJVCkxBAeQhsZa7DryRYdLHbnZe2P+pGVwsnHv5cqIxORo\n7YoX5r6FtNMJ2HVkc7uBM6mn9+NU0VHERixGyIiJHDgjA4bgAKCW1DiWdwi7Uj5HQ1Ndu3YjAxPM\nCn8C40ZN4cAXol6mVCgR7vsAfD3GYkfKRmQUHNZob2isxVcJ7+No7gHMnfgMHK1dZapUTAzBfu7S\n1WJsS/wYF8oLtbYHDY9EbMQSmJsM6uXKiOjXzE0G4YkHXkSoz2RsO/RRuxlnisry8PevX0TUmJl4\ncOwCGOobyVSpWBiC/dTN5nrsPboFyafitT7zZ2fpjLkTn4GXi68M1RFRR7xcRmPVo+/hUOb32H98\nu8ZSTWpJjUOZO3Gi8Cc8POEpTsrdCxiC/YwkSUgvSMLOnzZpXeVdT1cfDwTPxaTAWC76SdRH6enq\n4YGQeQj0noDtiZ8gryRTo72m4To27f0HUnN+xOyop+Ew2EWmSgc+hmA/UnqlCN8mfYbi8nyt7aM8\nQjA7cilnsSfqJ6wt7PHMQ68iu+gYvjv8Garrr2m0n7mUg3e+XoEo/xl4IGQ+p1+7DxiC/UB9Yy32\npH6F1NP7ta7zN9jcDo9ELoWvR4gM1RHRvVAoFPAbForhrv748dg2JJ7cpTEpt1qtwqHMncgoSMas\niCcQNDySA9x6EEOwD1OpVTiSsw/xaf/Bzeb6du26OnqYEjQbk4Me5mTXRP2cgZ4hZkU8gbEjJ2N7\n0icovKi5mn3tzWp8uf89pGTvw+zIpXC195Kp0oGFIdhHFV48he+SN6C86qLW9lHuwXgkcimsLex7\nuTIiup/sLJ3wbOxfcOpcGr5P/ne7S6QXKgrxv1v/iLEjJmFG+GOc8vAeMQT7mNrG6/hs99vILjqm\ntd1mkCNmRy6Fj1tgL1dGRL1FoVDA3zMMPm6BSMj4FgdPfI9W1S2NPsfyDyHrXCpigufCXHKCjpJf\n53eDR62PuNlcj4zzB1BQfhxqLY88GOgZIiZkHqL8Z0JPl6M+iUSgr2eA6eMWYazPJHyf/G/k3DYR\nfvOtJvyQ+gVMDSwQ6BaNQCmQj1TcIYagzFSqVhw5/SP2Ht2idbYXAAgZMREzwx/nZQ8iQVlb2OO3\nM19GQUkWvkvegIrrpRrt9c01OFz4LUrr8vDw+CW8X3gHGIIykSQJp8+nY2fK57hSfVlrHzd7bzwS\nuRRu/EATEYDhrv5YteifSMnZh71Ht7QbMFdclo//3fpHBHiNx8zwx/i4VDcwBGVwoeIMdv60CUVl\neVrbLU2tMSviCQR4jeelDSLSoKOji0j/GQgaHol9x7bip1Px7W6hZJ75CaeK0hDpNx1TgudwlYpO\nMAR70dUb5fgh9QtknU3V2q6r1MMo5zA8NuNZ6OvxkQci6piJoRlmR/4GEb4P4vPd7+FS9VmNdpWq\nFYcydyIt9wBigudgvN80PkqlBUOwF9TUX8e+49uQlpugscDmLxQKJUJ9JsPJ2AfG+mYMQCLqNjsr\nZ0zymY/yG+eRV3kEl29bR7SxuQE7Uz7H4azdmDp2AUJ8JkFHqSNPsX0QQ/A+amiqw8ETO3A46wfc\nam3R2sfHNQCzIp6Ao7UbMjIyerlCIhooHAa5Y/qkR5BekIQ9aV/jRn2VRvuN+ir85+BaHMzcgenj\nHoXfsFDOPAOG4H3R1NKIpJO7cChzJ5pabmrt42zrgdiIxfByGd3L1RHRQKVU6mCsz2SM8YrA4aw9\nSEjf3u476Er1ZWyM/zucbTwwfdwi+LiJ/VgFQ7AHNd9qQkr2Xhw48T0aGmu19rG2sMf0cYswxiuC\nv4UR0X2hr2uAKUGPIGxkNBIyvkXyqfh2D9tfulqMj3ethrvDcEwLXQgvl9FChiFDsAe03GpGSs4+\nHMz4TuvyRgBgZjwID46dj3Ejo7nEERH1ChMjc8SOX4IJfjOw99gWHM9PbLf+6PnyAqz9/jUMdfTB\n1NCFwq1ByhC8Bz+f+e3DocwdqLt5Q2sfYwNTTA58GBP8p8NAz7CXKyQiAqzMbfDolOcxOTAW8Uf/\no3WEelFZHj787lUMdRqJB0PmCXNmyBC8C43NN5GSvReHTu7s8LKngZ4hJo55CBMDZsHIwKSXKyQi\nas/eygVPTfsjSq8UIz7ta+ReaD8Yr+hyLtZ+/xrcHLzxQPDcAX/PkCF4B+pu1uBw1m78dGoPGjsY\n8KKvZ4hIv+mYFPAQTIzMe7lCIqKuudh64JmHXsGFijPYe3QL8m9b2R4ALpQX4uNdq+Fk446Y4Dnw\nGxoK5QB8tIIh2A1VtZVIzNyFo7kH0NLarLWPvp4hxo+eikkBsTAztujlComI7pybvReWx/4Z58sL\nsPfoFhRczGrX5/LV89gY/y5sBjlicmAsgodHQU9XX4Zq7w+GYCcuXS3GwYzvcfLsEa0rOwCAgb4R\nIv1mIGrMTJjyzI+I+iF3h+F49uG/4ELFGfx4bJvWy6RXb5Rhy8F12JP2NSL9ZyDC90EYG5rKUG3P\nYgjeRi2pkXf+BBJP7sLZSzkd9jM2NEOU/wxM8Js+ID4IRERu9l545qFXUHqlCPvTtyP73FFIkDT6\n1N28gd2pX2J/+naE+kxGpP8M2AxykKnie8cQ/K/mlkYcz0/E4VN7OlzVAQAsTKwwKSAWYaOmwEDf\nqBcrJCLqHS62Q7F0+ipUXC/FgYzvkFGY3G7Kx5ZbTUg+tQc/nYqH79AQRPrPwDCnUf1uEI3wIXj1\nRjl+yt6LY7kHOhzsAvw8P1904MMI9J7A5/yISAj2Vi54LOYFTAtdhKSsH5B6ej9abjVp9JEgIbvo\nGLKLjsFxsCsm+M9AkPeEfjMHspAhqFarkHvhBFKy92kdFfVrw5xGYmLAQxjpHsQZXohISFbmNnhk\nwlN4MGQeUrL3Ijk7HrUN1e36lVWVYMvBtdiZsgljfSYjwvcB2Fo6yVBx9wkVgjX113E07wDSTifg\net3VDvsplToI8IxA1JiZGGI3rBcrJCLqu4wNTRETMheTAmOReSYFiZk7261aAfy8ckXSyV1IOrkL\nXi6jEe77AHw9QvrkVbQBH4IqtQr5FzKRlpuA3PMZHY7yBABTIwuE+8Yg3PdBDDId3ItVEhH1H7o6\neggZMRHBw6Nw7nIukrN2I7v4eLsp2QDgTGk2zpRmw9TIAmN9JiLUJxp2Vs4yVK3dgA3ByurLOJZ7\nEMfzE1F7s/1p+6+52nkiYvRUBHhFDKjnX4iI7ieFQgFP51HwdB6F67VXkJK9D2l5B7TOpFXfWIOD\nJ3bg4Ikd8HAYgbEjJ8N/WBiMDIxlqPz/DKgQbGisReaZFBzPT0RJ5dlO++rp6CPAKwLj/abxkicR\n0T2yMrfFrIgnMDV0IbLOHcFP2XtxobxQa9/i8nwUl+dje9InGD00FCEjJsLbZbQsM9L0+xBsudWM\n0+fTkVFwGPklJ6FSt3ba32HwEISNikHw8Cg+30dE1MP0dPUQPDwKwcOjcPnqBaTl7kd6fpLW0fe3\nWltwojAZJwqTYW5siQCvCAR6T8AQu2G99qhFvwzBW60tyC85iZNnj+B08XE03zZk93YGeoYI8BqP\n0JGT4Wbv3e+eYyEi6o+cbNwwJ+ppzAp/EifPHsGxvIM4dzlXa9/am9VIyvoBSVk/wMbCAWO8IjDG\nMxyO1q739Tu734RgS2szCkpOIutcGk4Xp3e4YvuvDXMehbEjJsLfM5zLGBERyURfzwBjfSZhrM8k\nXL1RjmN5B5Gen4Tq+mta+1+tKcf+9G+wP/0b2Fk6w99zHPyGjYOTtXuPB2KfDsGbTfXIvXACOUXH\nkHfhRIeTV/+a7SBHBI+YiKDhEzDY3K4XqiQiou6yGeSAGWGPYdq4RTh3KRfH8w/h1Lm0Dq/oVVZf\nwo/Hv8GPx7/BYAs7+A0Nha9HCNwdhvfIPcQ+F4JXqi8j9/wJnD6fjqLLuZ0+0vALCxOrtmvJLrZD\nebmTiKiPUyqU8HLxhZeLL+ZNWobc8xk4UZiM3AsnoFJpH9tRVVOJQ5k7cShzJ0wMzTDSPQgj3YMx\nfIj/XY8y7TMh+O3hz5B3/gSu1pR3q7+ZkQX8PMMQ4BUBD8cRnM2FiKif0tc1wBjPcIzxDMfN5nrk\nFB3HyTMpKCg91W7O0l80NNXheH4ijucnQqnUgZezL5bF/vmOs6DPhODhrN1d9rEwHQy/oWMxeug4\nDHPyGZALPBIRiczYwLTt/mFDUx1OF6fj1Lk0FFzMQqvqltZ91GoVWtWtd3Uy1GdCsCP2Vi7w9QjB\nKI8QuNp78oyPiEgQJoZmbYHY1NKIvAsncLo4HbkXMtDY3KDRd6Rb0F39jD4XgjpKXQxzGgkf90CM\ndAuCraWj3CUREZHMDPWNEOAVgQCvCKhUrSgqy0Pu+QzkXjiBK9WXMdI98K5et8+EYLjvgxjh6g8v\nFz8Ycp0+IiLqgI6OLrxcRsPLZTQenvAUqmoqYWVue1ev1WdCcP6kZXKXQERE/dBgi7t/HI432IiI\nSFgMQSIiEhZDkIiIhMUQJCIiYTEEiYhIWAxBIiISVrdCcN26dXB3d4eRkRGCgoKQkpLSaf+cnBxE\nRkbC2NgYzs7OePPNN3ukWCIiop7UZQhu3boVK1aswCuvvIKsrCyEhYVh6tSpKC0t1dq/trYWU6ZM\ngYODAzIyMvDee+/h3XffxZo1a3q8eCIionvRZQiuWbMGS5YswdKlS+Ht7Y33338fDg4OWL9+vdb+\nX331FZqamvD555/Dx8cHs2fPxqpVqxiCRETU53Qagi0tLcjMzERMTIzG9piYGKSmpmrdJy0tDePH\nj4eBgYFG/7KyMpSUlPRAyURERD2j02nTrl27BpVKBTs7zSlpbG1tUVFRoXWfiooKDBkyRGPbL/tX\nVFTA1dVV6341NTXdLnqg8vT0BMBjAfBY3I7HQxOPhyYej7vX46NDuao7ERH1F52GoLW1NXR0dFBZ\nWamxvbKyEg4ODlr3sbe3b3eW+Mv+9vb291IrERFRj+r0cqi+vj4CAwOxf/9+zJ49u217QkIC5s6d\nq3WfcePGYdWqVWhubm67L5iQkAAnJ6d2l0ItLCzutX4iIqK71uXl0JUrV2LTpk3YsGED8vPz8cIL\nL6CiogLLlv289FFcXByio6Pb+i9atAjGxsZYvHgxcnNz8d133+Gdd97BypUr79+7ICIiugtdric4\nb948VFVVYfXq1SgvL4evry/i4+Ph4uIC4OfBLsXFxW39zc3NkZCQgN/97ncICgqClZUVXnrpJbz4\n4ov3710QERHdBYUkSZLcRRAREcmh1+cOffvttxEcHAwLCwvY2tpi1qxZyM3N7e0y+oy1a9fCz88P\nFhYWsLCwQFhYGOLj4+Uuq894++23oVQq8fzzz8tdiiz+8pe/QKlUavxxdHSUuyxZlZeX48knn4St\nrS2MjIwwcuRIJCcny12WLNzc3Np9PpRKJWbMmCF3ab2utbUVL7/8Mjw8PGBkZAQPDw+8+uqrUKlU\nne7X5eXQnnb48GE899xzCA4Ohlqtxp///GdER0cjLy8PlpaWvV2O7FxcXPD3v/8dnp6eUKvV2LRp\nE2JjY5Geng4/Pz+5y5PV0aNH8emnn2L06NFCP3ozfPhwJCUltf1bR0dHvmJkduPGDYSHh2PChAmI\nj4+HjY0NiouLYWtrK3dpsjhx4oTGl3xZWRkCAwMxf/58GauSx1tvvYWPP/4Ymzdvhq+vL06dOoXF\nixfDwMAAr7zySsc7SjKrr6+XdHR0pN27d8tdSp9hZWUlffLJJ3KXIasbN25IQ4cOlZKSkqSoqCjp\n+eefl7skWbz22mvSqFGj5C6jz4iLi5MiIiLkLqPPWr16tWRpaSk1NTXJXUqvmzFjhrR48WKNbU88\n8YQ0c+bMTveTfSml2tpaqNVqIc8Cb6dSqbBlyxY0NTVhwoQJcpcjq6effhpz585FZGQkJMFvWxcX\nF8PJyQkeHh5YuHAhzp8/L3dJstmxYwdCQkIwf/582NnZYcyYMVi7dq3cZfUJkiRhw4YNeOyxxzSm\nrRTF1KlTcejQIRQWFgIA8vLykJiYiGnTpnW+4/3L5e6ZO3euFBAQIKnVarlLkU12drZkYmIi6erq\nSmZmZsKfFX/yySdSUFCQ1NraKkmSJPSZ4N69e6VvvvlGysnJkQ4cOCBFRUVJ9vb2UlVVldylycLA\nwEAyNDSUXn75ZSkrK0vauHGjZGpqKn344Ydylya7H3/8UVIoFFJ2drbcpcgmLi5OUigUkp6enqRQ\nKKRXX321y31kDcEXX3xRcnJyks6fPy9nGbJraWmRioqKpMzMTCkuLk4yNTWV0tPT5S5LFgUFBZKN\njY1UWFjYti0yMlJ67rnnZKyq72hoaJBsbW2lNWvWyF2KLPT09KTw8HCNbS+//LI0YsQImSrqO+bM\nmSONHTtW7jJk895770n29vbS1q1bpdOnT0tffPGFZGVlJW3YsKHT/WQLwRUrVkiOjo4aX3b0s+jo\n6HbXtkWxceNGSaFQSLq6um1/FAqFpFQqJT09PamlpUXuEmU3ceJE6dlnn5W7DFm4urpKv/3tbzW2\nbd68WTIxMZGpor6hsrJS0tfXlz777DO5S5GNra2t9P7772tsW716tTRs2LBO9+v10aEA8MILL+Cb\nb75BYmIivLy85CihT1OpVFCr1XKXIYuHH34YISEhbf+WJAlLliyBl5cXXn75Zejp6clYnfyampqQ\nn5+PSZMmyV2KLMLDw1FQUKCx7cyZM3Bzc5OnoD5i06ZNMDQ0xMKFC+UuRTaSJEGp1BzmolQquxxT\n0Osh+Lvf/Q5ffvklduzYAQsLi7bJts3MzGBiYtLb5cjuT3/6E2bMmAFnZ2fU1dXh66+/xuHDh7Fv\n3z65S5PFL89L/pqxsTEsLS3h4+MjU1XyeemllzBr1iy4uLjgypUrePPNN9HY2Ignn3xS7tJk8eKL\nLyIsLAxvvfUW5s2bh5MnT+KDDz7A22+/LXdpspEkCZ999hkWLFgAY2NjucuRTWxsLP72t7/B3d0d\nPj4+OHnyJP75z392/f/K/To17cgvl7YUCoXGn9dff723S+kTFi9eLLm6ukoGBgaSra2tNGXKFGn/\n/v1yl9WniDwwZsGCBZKjo6Okr68vOTk5SXPmzJHy8/PlLktWe/bskfz8/CRDQ0PJ29tb+uCDD+Qu\nSVaHDh2SlEqlsOMIflFfXy/94Q9/kNzc3CQjIyPJw8ND+p//+R+pubm50/04bRoREQlL9ucEiYiI\n5MIQJCIiYTEEiYhIWAxBIiISFkOQiIiExRAkIiJhMQSJiEhYDEEiIhLW/weD5xVc2/tVUQAAAABJ\nRU5ErkJggg==\n",
"text": [
""
]
}
],
"prompt_number": 4
},
{
"cell_type": "code",
"collapsed": false,
"input": [],
"language": "python",
"metadata": {},
"outputs": [],
"prompt_number": 2
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
""
]
}
],
"metadata": {}
}
]
}