PageRank.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# PageRank"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "<h3>Part 1 - Worksheet</h3><br>\n",
    "<b>Introduction</b> <br>\n",
    "PageRank (developed by Larry Page and Sergey Brin) revolutionized web search by generating a ranked list of web pages based on the underlying connectivity of the web. The PageRank algorithm is based on an ideal random web surfer who, when reaching a page, goes to the next page by clicking on a link. The surfer has equal probability of clicking any link on the page and, when reaching a page with no links, has equal probability of moving to any other page by typing in its URL. In addition, the surfer may occasionally choose to type in a random URL instead of following the links on a page. The PageRank is the ranked order of the pages from the most to the least probable page the surfer will be viewing."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Populating the interactive namespace from numpy and matplotlib\n"
     ]
    }
   ],
   "source": [
    "# Before we begin, let's load the libraries.\n",
    "%pylab notebook \n",
    "#%pylab is a magic function which imports numpy and matplotlib\n",
    "import numpy as np\n",
    "import numpy.linalg as la\n",
    "from readonly.PageRankFunctions import *\n",
    "np.set_printoptions(suppress=True)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "<h3>PageRank as a linear algebra problem</h3>"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Let's imagine a micro-internet, with just 6 websites (Avocado, Bullseye, CatBabel, Dromeda, eTings, and FaceSpace). Each website links to some of the others, and this forms a network as shown,\n",
    "\n",
    "![A Micro-Internet](readonly/internet.png \"A Micro-Internet\")\n",
    "\n",
    "The design principle of PageRank is that important websites will be linked to by important websites.\n",
    "This somewhat recursive principle will form the basis of our thinking.\n",
    "\n",
    "Imagine we have 100 *Procrastinating Pat*s on our micro-internet, each viewing a single website at a time.\n",
    "Each minute the Pats follow a link on their website to another site on the micro-internet.\n",
    "After a while, the websites that are most linked to will have more Pats visiting them, and in the long run, each minute for every Pat that leaves a website, another will enter keeping the total numbers of Pats on each website constant.\n",
    "The PageRank is simply the ranking of websites by how many Pats they have on them at the end of this process.\n",
    "\n",
    "We represent the number of Pats on each website with the vector,\n",
    "$$\\mathbf{r} = \\begin{bmatrix} r_A \\\\ r_B \\\\ r_C \\\\ r_D \\\\ r_E \\\\ r_F \\end{bmatrix}$$\n",
    "And say that the number of Pats on each website in minute $i+1$ is related to those at minute $i$ by the matrix transformation\n",
    "\n",
    "$$ \\mathbf{r}^{(i+1)} = L \\,\\mathbf{r}^{(i)}$$\n",
    "with the matrix $L$ taking the form,\n",
    "$$ L = \\begin{bmatrix}\n",
    "L_{A→A} & L_{B→A} & L_{C→A} & L_{D→A} & L_{E→A} & L_{F→A} \\\\\n",
    "L_{A→B} & L_{B→B} & L_{C→B} & L_{D→B} & L_{E→B} & L_{F→B} \\\\\n",
    "L_{A→C} & L_{B→C} & L_{C→C} & L_{D→C} & L_{E→C} & L_{F→C} \\\\\n",
    "L_{A→D} & L_{B→D} & L_{C→D} & L_{D→D} & L_{E→D} & L_{F→D} \\\\\n",
    "L_{A→E} & L_{B→E} & L_{C→E} & L_{D→E} & L_{E→E} & L_{F→E} \\\\\n",
    "L_{A→F} & L_{B→F} & L_{C→F} & L_{D→F} & L_{E→F} & L_{F→F} \\\\\n",
    "\\end{bmatrix}\n",
    "$$\n",
    "where the columns represent the probability of leaving a website for any other website, and sum to one.\n",
    "The rows determine how likely you are to enter a website from any other, though these need not add to one.\n",
    "The long time behaviour of this system is when $ \\mathbf{r}^{(i+1)} = \\mathbf{r}^{(i)}$, so we'll drop the superscripts here, and that allows us to write,\n",
    "$$ L \\,\\mathbf{r} = \\mathbf{r}$$\n",
    "\n",
    "which is an eigenvalue equation for the matrix $L$, with eigenvalue 1 (this is guaranteed by the probabalistic structure of the matrix $L$).\n",
    "\n",
    "Complete the matrix $L$ below, we've left out the column for which websites the *FaceSpace* website (F) links to.\n",
    "Remember, this is the probability to click on another website from this one, so each column should add to one (by scaling by the number of links).\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Replace the ??? here with the probability of clicking a link to each website when leaving Website F (FaceSpace).\n",
    "L = np.array([[0,   1/2, 1/3, 0, 0,   0 ],\n",
    "              [1/3, 0,   0,   0, 1/2, 0 ],\n",
    "              [1/3, 1/2, 0,   1, 0,   1/2 ],\n",
    "              [1/3, 0,   1/3, 0, 1/2, 1/2 ],\n",
    "              [0,   0,   0,   0, 0,   0 ],\n",
    "              [0,   0,   1/3, 0, 0,   0 ]])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([16.        ,  5.33333333, 40.        , 25.33333333,  0.        ,\n",
       "       13.33333333])"
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "eVals, eVecs = la.eig(L) # Gets the eigenvalues and vectors\n",
    "order = np.absolute(eVals).argsort()[::-1] # Orders them by their eigenvalues\n",
    "eVals = eVals[order]\n",
    "eVecs = eVecs[:,order]\n",
    "\n",
    "r = eVecs[:, 0] # Sets r to be the principal eigenvector\n",
    "100 * np.real(r / np.sum(r)) # Make this eigenvector sum to one, then multiply by 100 Procrastinating Pats"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We can see from this list, the number of Procrastinating Pats that we expect to find on each website after long times.\n",
    "Putting them in order of *popularity* (based on this metric), the PageRank of this micro-internet is:\n",
    "\n",
    "**C**atBabel, **D**romeda, **A**vocado, **F**aceSpace, **B**ullseye, **e**Tings\n",
    "\n",
    "Referring back to the micro-internet diagram, is this what you would have expected?\n",
    "Convince yourself that based on which pages seem important given which others link to them, that this is a sensible ranking.\n",
    "\n",
    "Let's now try to get the same result using the Power-Iteration method that was covered in the video.\n",
    "This method will be much better at dealing with large systems.\n",
    "\n",
    "First let's set up our initial vector, $\\mathbf{r}^{(0)}$, so that we have our 100 Procrastinating Pats equally distributed on each of our 6 websites."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([16.66666667, 16.66666667, 16.66666667, 16.66666667, 16.66666667,\n",
       "       16.66666667])"
      ]
     },
     "execution_count": 5,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "r = 100 * np.ones(6) / 6 # Sets up this vector (6 entries of 1/6 × 100 each)\n",
    "r # Shows it's value"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 55,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([16.        ,  5.33333333, 40.        , 25.33333333,  0.        ,\n",
       "       13.33333333])"
      ]
     },
     "execution_count": 55,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "r = L @ r # Apply matrix L to r\n",
    "r # Show it's value\n",
    "# Re-run this cell multiple times to converge to the correct answer."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We can automate applying this matrix multiple times as follows,"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 56,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([16.        ,  5.33333333, 40.        , 25.33333333,  0.        ,\n",
       "       13.33333333])"
      ]
     },
     "execution_count": 56,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "r = 100 * np.ones(6) / 6 # Sets up this vector (6 entries of 1/6 × 100 each)\n",
    "for i in np.arange(100) : # Repeat 100 times\n",
    "    r = L @ r\n",
    "r"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Or even better, we can keep running until we get to the required tolerance."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 57,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "18 iterations to convergence.\n"
     ]
    },
    {
     "data": {
      "text/plain": [
       "array([16.00149917,  5.33252025, 39.99916911, 25.3324738 ,  0.        ,\n",
       "       13.33433767])"
      ]
     },
     "execution_count": 57,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "r = 100 * np.ones(6) / 6 # Sets up this vector (6 entries of 1/6 × 100 each)\n",
    "lastR = r\n",
    "r = L @ r\n",
    "i = 0\n",
    "while la.norm(lastR - r) > 0.01 :\n",
    "    lastR = r\n",
    "    r = L @ r\n",
    "    i += 1\n",
    "print(str(i) + \" iterations to convergence.\")\n",
    "r"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "See how the PageRank order is established fairly quickly, and the vector converges on the value we calculated earlier after a few tens of repeats.\n",
    "\n",
    "Congratulations! You've just calculated your first PageRank!"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Damping Parameter\n",
    "The system we just studied converged fairly quickly to the correct answer.\n",
    "Let's consider an extension to our micro-internet where things start to go wrong.\n",
    "\n",
    "Say a new website is added to the micro-internet: *Geoff's* Website.\n",
    "This website is linked to by *FaceSpace* and only links to itself.\n",
    "![An Expanded Micro-Internet](readonly/internet2.png \"An Expanded Micro-Internet\")\n",
    "\n",
    "Intuitively, only *FaceSpace*, which is in the bottom half of the page rank, links to this website amongst the two others it links to,\n",
    "so we might expect *Geoff's* site to have a correspondingly low PageRank score.\n",
    "\n",
    "Build the new $L$ matrix for the expanded micro-internet, and use Power-Iteration on the Procrastinating Pat vector.\n",
    "See what happens…"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 58,
   "metadata": {},
   "outputs": [],
   "source": [
    " # We'll call this one L2, to distinguish it from the previous L.\n",
    "L2 = np.array([[0,   1/2, 1/3, 0, 0,   0, 0 ],\n",
    "               [1/3, 0,   0,   0, 1/2, 0, 0 ],\n",
    "               [1/3, 1/2, 0,   1, 0,   1/3, 0 ],\n",
    "               [1/3, 0,   1/3, 0, 1/2, 1/3, 0 ],\n",
    "               [0,   0,   0,   0, 0,   0, 0 ],\n",
    "               [0,   0,   1/3, 0, 0,   0, 0 ],\n",
    "               [0,   0,   0,   0, 0,   1/3, 1 ]])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 59,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "131 iterations to convergence.\n"
     ]
    },
    {
     "data": {
      "text/plain": [
       "array([ 0.03046998,  0.01064323,  0.07126612,  0.04423198,  0.        ,\n",
       "        0.02489342, 99.81849527])"
      ]
     },
     "execution_count": 59,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "r = 100 * np.ones(7) / 7 # Sets up this vector (6 entries of 1/6 × 100 each)\n",
    "lastR = r\n",
    "r = L2 @ r\n",
    "i = 0\n",
    "while la.norm(lastR - r) > 0.01 :\n",
    "    lastR = r\n",
    "    r = L2 @ r\n",
    "    i += 1\n",
    "print(str(i) + \" iterations to convergence.\")\n",
    "r"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "That's no good! *Geoff* seems to be taking all the traffic on the micro-internet, and somehow coming at the top of the PageRank.\n",
    "This behaviour can be understood, because once a Pat get's to *Geoff's* Website, they can't leave, as all links head back to Geoff.\n",
    "\n",
    "To combat this, we can add a small probability that the Procrastinating Pats don't follow any link on a webpage, but instead visit a website on the micro-internet at random.\n",
    "We'll say the probability of them following a link is $d$ and the probability of choosing a random website is therefore $1-d$.\n",
    "We can use a new matrix to work out where the Pat's visit each minute.\n",
    "$$ M = d \\, L + \\frac{1-d}{n} \\, J $$\n",
    "where $J$ is an $n\\times n$ matrix where every element is one.\n",
    "\n",
    "If $d$ is one, we have the case we had previously, whereas if $d$ is zero, we will always visit a random webpage and therefore all webpages will be equally likely and equally ranked.\n",
    "For this extension to work best, $1-d$ should be somewhat small - though we won't go into a discussion about exactly how small.\n",
    "\n",
    "Let's retry this PageRank with this extension."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 71,
   "metadata": {},
   "outputs": [],
   "source": [
    "d = 0.5 # Feel free to play with this parameter after running the code once.\n",
    "M = d * L2 + (1-d)/7 * np.ones([7, 7]) # np.ones() is the J matrix, with ones for each entry."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 72,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "8 iterations to convergence.\n"
     ]
    },
    {
     "data": {
      "text/plain": [
       "array([13.68217054, 11.20902965, 22.41964343, 16.7593433 ,  7.14285714,\n",
       "       10.87976354, 17.90719239])"
      ]
     },
     "execution_count": 72,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "r = 100 * np.ones(7) / 7 # Sets up this vector (6 entries of 1/6 × 100 each)\n",
    "lastR = r\n",
    "r = M @ r\n",
    "i = 0\n",
    "while la.norm(lastR - r) > 0.01 :\n",
    "    lastR = r\n",
    "    r = M @ r\n",
    "    i += 1\n",
    "print(str(i) + \" iterations to convergence.\")\n",
    "r"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "This is certainly better, the PageRank gives sensible numbers for the Procrastinating Pats that end up on each webpage.\n",
    "This method still predicts Geoff has a high ranking webpage however.\n",
    "This could be seen as a consequence of using a small network. We could also get around the problem by not counting self-links when producing the L matrix (an if a website has no outgoing links, make it link to all websites equally).\n",
    "We won't look further down this route, as this is in the realm of improvements to PageRank, rather than eigenproblems.\n",
    "\n",
    "You are now in a good position, having gained an understanding of PageRank, to produce your own code to calculate the PageRank of a website with thousands of entries.\n",
    "\n",
    "Good Luck!"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Assessment"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 62,
   "metadata": {},
   "outputs": [],
   "source": [
    "# PACKAGE\n",
    "# Here are the imports again, just in case you need them.\n",
    "# There is no need to edit or submit this cell.\n",
    "import numpy as np\n",
    "import numpy.linalg as la\n",
    "from readonly.PageRankFunctions import *\n",
    "np.set_printoptions(suppress=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 63,
   "metadata": {},
   "outputs": [],
   "source": [
    "# GRADED FUNCTION\n",
    "# Complete this function to provide the PageRank for an arbitrarily sized internet.\n",
    "# I.e. the principal eigenvector of the damped system, using the power iteration method.\n",
    "# (Normalisation doesn't matter here)\n",
    "# The functions inputs are the linkMatrix, and d the damping parameter - as defined in this worksheet.\n",
    "# (The damping parameter, d, will be set by the function - no need to set this yourself.)\n",
    "def pageRank(linkMatrix, d) :\n",
    "    n = linkMatrix.shape[0]\n",
    "    M = d * linkMatrix + (1-d)/n * np.ones([n, n]) # np.ones() is the J matrix, with ones for each entry.\n",
    "    r = 100 * np.ones(n) / n # Sets up this vector (n entries of 1/n × 100 each)\n",
    "    lastR = r\n",
    "    r = M @ r\n",
    "    i = 0\n",
    "    while la.norm(lastR - r) > 0.01 :\n",
    "        lastR = r\n",
    "        r = M @ r\n",
    "        i += 1\n",
    "    print(str(i) + \" iterations to convergence.\")\n",
    "    return r\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 64,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([[0.2, 0. , 0. , 0. , 0.2],\n",
       "       [0.2, 0.5, 0.5, 0. , 0.2],\n",
       "       [0.2, 0.5, 0. , 0. , 0.2],\n",
       "       [0.2, 0. , 0.5, 0.5, 0.2],\n",
       "       [0.2, 0. , 0. , 0.5, 0.2]])"
      ]
     },
     "execution_count": 64,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# Use the following function to generate internets of different sizes.\n",
    "generate_internet(5)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 65,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Test your PageRank method against the built in \"eig\" method.\n",
    "# You should see yours is a lot faster for large internets\n",
    "L = generate_internet(10)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 66,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "13 iterations to convergence.\n"
     ]
    },
    {
     "data": {
      "text/plain": [
       "array([ 3.70374804, 18.52134933,  3.70374804,  7.40856259,  3.70374804,\n",
       "        3.70374804,  3.70374804, 11.11256037, 18.51520347, 25.92358404])"
      ]
     },
     "execution_count": 66,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "pageRank(L, 1)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 67,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([ 3.70370371, 18.51851852,  3.70370371,  7.40740741,  3.70370371,\n",
       "        3.70370371,  3.70370371, 11.11111111, 18.51851851, 25.92592591])"
      ]
     },
     "execution_count": 67,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# Do note, this is calculating the eigenvalues of the link matrix, L,\n",
    "# without any damping. It may give different results that your pageRank function.\n",
    "# If you wish, you could modify this cell to include damping.\n",
    "# (There is no credit for this though)\n",
    "eVals, eVecs = la.eig(L) # Gets the eigenvalues and vectors\n",
    "order = np.absolute(eVals).argsort()[::-1] # Orders them by their eigenvalues\n",
    "eVals = eVals[order]\n",
    "eVecs = eVecs[:,order]\n",
    "\n",
    "r = eVecs[:, 0]\n",
    "100 * np.real(r / np.sum(r))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 68,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Populating the interactive namespace from numpy and matplotlib\n",
      "36 iterations to convergence.\n"
     ]
    },
    {
     "data": {
      "application/javascript": [
       "/* Put everything inside the global mpl namespace */\n",
       "window.mpl = {};\n",
       "\n",
       "\n",
       "mpl.get_websocket_type = function() {\n",
       "    if (typeof(WebSocket) !== 'undefined') {\n",
       "        return WebSocket;\n",
       "    } else if (typeof(MozWebSocket) !== 'undefined') {\n",
       "        return MozWebSocket;\n",
       "    } else {\n",
       "        alert('Your browser does not have WebSocket support.' +\n",
       "              'Please try Chrome, Safari or Firefox ≥ 6. ' +\n",
       "              'Firefox 4 and 5 are also supported but you ' +\n",
       "              'have to enable WebSockets in about:config.');\n",
       "    };\n",
       "}\n",
       "\n",
       "mpl.figure = function(figure_id, websocket, ondownload, parent_element) {\n",
       "    this.id = figure_id;\n",
       "\n",
       "    this.ws = websocket;\n",
       "\n",
       "    this.supports_binary = (this.ws.binaryType != undefined);\n",
       "\n",
       "    if (!this.supports_binary) {\n",
       "        var warnings = document.getElementById(\"mpl-warnings\");\n",
       "        if (warnings) {\n",
       "            warnings.style.display = 'block';\n",
       "            warnings.textContent = (\n",
       "                \"This browser does not support binary websocket messages. \" +\n",
       "                    \"Performance may be slow.\");\n",
       "        }\n",
       "    }\n",
       "\n",
       "    this.imageObj = new Image();\n",
       "\n",
       "    this.context = undefined;\n",
       "    this.message = undefined;\n",
       "    this.canvas = undefined;\n",
       "    this.rubberband_canvas = undefined;\n",
       "    this.rubberband_context = undefined;\n",
       "    this.format_dropdown = undefined;\n",
       "\n",
       "    this.image_mode = 'full';\n",
       "\n",
       "    this.root = $('<div/>');\n",
       "    this._root_extra_style(this.root)\n",
       "    this.root.attr('style', 'display: inline-block');\n",
       "\n",
       "    $(parent_element).append(this.root);\n",
       "\n",
       "    this._init_header(this);\n",
       "    this._init_canvas(this);\n",
       "    this._init_toolbar(this);\n",
       "\n",
       "    var fig = this;\n",
       "\n",
       "    this.waiting = false;\n",
       "\n",
       "    this.ws.onopen =  function () {\n",
       "            fig.send_message(\"supports_binary\", {value: fig.supports_binary});\n",
       "            fig.send_message(\"send_image_mode\", {});\n",
       "            if (mpl.ratio != 1) {\n",
       "                fig.send_message(\"set_dpi_ratio\", {'dpi_ratio': mpl.ratio});\n",
       "            }\n",
       "            fig.send_message(\"refresh\", {});\n",
       "        }\n",
       "\n",
       "    this.imageObj.onload = function() {\n",
       "            if (fig.image_mode == 'full') {\n",
       "                // Full images could contain transparency (where diff images\n",
       "                // almost always do), so we need to clear the canvas so that\n",
       "                // there is no ghosting.\n",
       "                fig.context.clearRect(0, 0, fig.canvas.width, fig.canvas.height);\n",
       "            }\n",
       "            fig.context.drawImage(fig.imageObj, 0, 0);\n",
       "        };\n",
       "\n",
       "    this.imageObj.onunload = function() {\n",
       "        fig.ws.close();\n",
       "    }\n",
       "\n",
       "    this.ws.onmessage = this._make_on_message_function(this);\n",
       "\n",
       "    this.ondownload = ondownload;\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._init_header = function() {\n",
       "    var titlebar = $(\n",
       "        '<div class=\"ui-dialog-titlebar ui-widget-header ui-corner-all ' +\n",
       "        'ui-helper-clearfix\"/>');\n",
       "    var titletext = $(\n",
       "        '<div class=\"ui-dialog-title\" style=\"width: 100%; ' +\n",
       "        'text-align: center; padding: 3px;\"/>');\n",
       "    titlebar.append(titletext)\n",
       "    this.root.append(titlebar);\n",
       "    this.header = titletext[0];\n",
       "}\n",
       "\n",
       "\n",
       "\n",
       "mpl.figure.prototype._canvas_extra_style = function(canvas_div) {\n",
       "\n",
       "}\n",
       "\n",
       "\n",
       "mpl.figure.prototype._root_extra_style = function(canvas_div) {\n",
       "\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._init_canvas = function() {\n",
       "    var fig = this;\n",
       "\n",
       "    var canvas_div = $('<div/>');\n",
       "\n",
       "    canvas_div.attr('style', 'position: relative; clear: both; outline: 0');\n",
       "\n",
       "    function canvas_keyboard_event(event) {\n",
       "        return fig.key_event(event, event['data']);\n",
       "    }\n",
       "\n",
       "    canvas_div.keydown('key_press', canvas_keyboard_event);\n",
       "    canvas_div.keyup('key_release', canvas_keyboard_event);\n",
       "    this.canvas_div = canvas_div\n",
       "    this._canvas_extra_style(canvas_div)\n",
       "    this.root.append(canvas_div);\n",
       "\n",
       "    var canvas = $('<canvas/>');\n",
       "    canvas.addClass('mpl-canvas');\n",
       "    canvas.attr('style', \"left: 0; top: 0; z-index: 0; outline: 0\")\n",
       "\n",
       "    this.canvas = canvas[0];\n",
       "    this.context = canvas[0].getContext(\"2d\");\n",
       "\n",
       "    var backingStore = this.context.backingStorePixelRatio ||\n",
       "\tthis.context.webkitBackingStorePixelRatio ||\n",
       "\tthis.context.mozBackingStorePixelRatio ||\n",
       "\tthis.context.msBackingStorePixelRatio ||\n",
       "\tthis.context.oBackingStorePixelRatio ||\n",
       "\tthis.context.backingStorePixelRatio || 1;\n",
       "\n",
       "    mpl.ratio = (window.devicePixelRatio || 1) / backingStore;\n",
       "\n",
       "    var rubberband = $('<canvas/>');\n",
       "    rubberband.attr('style', \"position: absolute; left: 0; top: 0; z-index: 1;\")\n",
       "\n",
       "    var pass_mouse_events = true;\n",
       "\n",
       "    canvas_div.resizable({\n",
       "        start: function(event, ui) {\n",
       "            pass_mouse_events = false;\n",
       "        },\n",
       "        resize: function(event, ui) {\n",
       "            fig.request_resize(ui.size.width, ui.size.height);\n",
       "        },\n",
       "        stop: function(event, ui) {\n",
       "            pass_mouse_events = true;\n",
       "            fig.request_resize(ui.size.width, ui.size.height);\n",
       "        },\n",
       "    });\n",
       "\n",
       "    function mouse_event_fn(event) {\n",
       "        if (pass_mouse_events)\n",
       "            return fig.mouse_event(event, event['data']);\n",
       "    }\n",
       "\n",
       "    rubberband.mousedown('button_press', mouse_event_fn);\n",
       "    rubberband.mouseup('button_release', mouse_event_fn);\n",
       "    // Throttle sequential mouse events to 1 every 20ms.\n",
       "    rubberband.mousemove('motion_notify', mouse_event_fn);\n",
       "\n",
       "    rubberband.mouseenter('figure_enter', mouse_event_fn);\n",
       "    rubberband.mouseleave('figure_leave', mouse_event_fn);\n",
       "\n",
       "    canvas_div.on(\"wheel\", function (event) {\n",
       "        event = event.originalEvent;\n",
       "        event['data'] = 'scroll'\n",
       "        if (event.deltaY < 0) {\n",
       "            event.step = 1;\n",
       "        } else {\n",
       "            event.step = -1;\n",
       "        }\n",
       "        mouse_event_fn(event);\n",
       "    });\n",
       "\n",
       "    canvas_div.append(canvas);\n",
       "    canvas_div.append(rubberband);\n",
       "\n",
       "    this.rubberband = rubberband;\n",
       "    this.rubberband_canvas = rubberband[0];\n",
       "    this.rubberband_context = rubberband[0].getContext(\"2d\");\n",
       "    this.rubberband_context.strokeStyle = \"#000000\";\n",
       "\n",
       "    this._resize_canvas = function(width, height) {\n",
       "        // Keep the size of the canvas, canvas container, and rubber band\n",
       "        // canvas in synch.\n",
       "        canvas_div.css('width', width)\n",
       "        canvas_div.css('height', height)\n",
       "\n",
       "        canvas.attr('width', width * mpl.ratio);\n",
       "        canvas.attr('height', height * mpl.ratio);\n",
       "        canvas.attr('style', 'width: ' + width + 'px; height: ' + height + 'px;');\n",
       "\n",
       "        rubberband.attr('width', width);\n",
       "        rubberband.attr('height', height);\n",
       "    }\n",
       "\n",
       "    // Set the figure to an initial 600x600px, this will subsequently be updated\n",
       "    // upon first draw.\n",
       "    this._resize_canvas(600, 600);\n",
       "\n",
       "    // Disable right mouse context menu.\n",
       "    $(this.rubberband_canvas).bind(\"contextmenu\",function(e){\n",
       "        return false;\n",
       "    });\n",
       "\n",
       "    function set_focus () {\n",
       "        canvas.focus();\n",
       "        canvas_div.focus();\n",
       "    }\n",
       "\n",
       "    window.setTimeout(set_focus, 100);\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._init_toolbar = function() {\n",
       "    var fig = this;\n",
       "\n",
       "    var nav_element = $('<div/>')\n",
       "    nav_element.attr('style', 'width: 100%');\n",
       "    this.root.append(nav_element);\n",
       "\n",
       "    // Define a callback function for later on.\n",
       "    function toolbar_event(event) {\n",
       "        return fig.toolbar_button_onclick(event['data']);\n",
       "    }\n",
       "    function toolbar_mouse_event(event) {\n",
       "        return fig.toolbar_button_onmouseover(event['data']);\n",
       "    }\n",
       "\n",
       "    for(var toolbar_ind in mpl.toolbar_items) {\n",
       "        var name = mpl.toolbar_items[toolbar_ind][0];\n",
       "        var tooltip = mpl.toolbar_items[toolbar_ind][1];\n",
       "        var image = mpl.toolbar_items[toolbar_ind][2];\n",
       "        var method_name = mpl.toolbar_items[toolbar_ind][3];\n",
       "\n",
       "        if (!name) {\n",
       "            // put a spacer in here.\n",
       "            continue;\n",
       "        }\n",
       "        var button = $('<button/>');\n",
       "        button.addClass('ui-button ui-widget ui-state-default ui-corner-all ' +\n",
       "                        'ui-button-icon-only');\n",
       "        button.attr('role', 'button');\n",
       "        button.attr('aria-disabled', 'false');\n",
       "        button.click(method_name, toolbar_event);\n",
       "        button.mouseover(tooltip, toolbar_mouse_event);\n",
       "\n",
       "        var icon_img = $('<span/>');\n",
       "        icon_img.addClass('ui-button-icon-primary ui-icon');\n",
       "        icon_img.addClass(image);\n",
       "        icon_img.addClass('ui-corner-all');\n",
       "\n",
       "        var tooltip_span = $('<span/>');\n",
       "        tooltip_span.addClass('ui-button-text');\n",
       "        tooltip_span.html(tooltip);\n",
       "\n",
       "        button.append(icon_img);\n",
       "        button.append(tooltip_span);\n",
       "\n",
       "        nav_element.append(button);\n",
       "    }\n",
       "\n",
       "    var fmt_picker_span = $('<span/>');\n",
       "\n",
       "    var fmt_picker = $('<select/>');\n",
       "    fmt_picker.addClass('mpl-toolbar-option ui-widget ui-widget-content');\n",
       "    fmt_picker_span.append(fmt_picker);\n",
       "    nav_element.append(fmt_picker_span);\n",
       "    this.format_dropdown = fmt_picker[0];\n",
       "\n",
       "    for (var ind in mpl.extensions) {\n",
       "        var fmt = mpl.extensions[ind];\n",
       "        var option = $(\n",
       "            '<option/>', {selected: fmt === mpl.default_extension}).html(fmt);\n",
       "        fmt_picker.append(option)\n",
       "    }\n",
       "\n",
       "    // Add hover states to the ui-buttons\n",
       "    $( \".ui-button\" ).hover(\n",
       "        function() { $(this).addClass(\"ui-state-hover\");},\n",
       "        function() { $(this).removeClass(\"ui-state-hover\");}\n",
       "    );\n",
       "\n",
       "    var status_bar = $('<span class=\"mpl-message\"/>');\n",
       "    nav_element.append(status_bar);\n",
       "    this.message = status_bar[0];\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.request_resize = function(x_pixels, y_pixels) {\n",
       "    // Request matplotlib to resize the figure. Matplotlib will then trigger a resize in the client,\n",
       "    // which will in turn request a refresh of the image.\n",
       "    this.send_message('resize', {'width': x_pixels, 'height': y_pixels});\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.send_message = function(type, properties) {\n",
       "    properties['type'] = type;\n",
       "    properties['figure_id'] = this.id;\n",
       "    this.ws.send(JSON.stringify(properties));\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.send_draw_message = function() {\n",
       "    if (!this.waiting) {\n",
       "        this.waiting = true;\n",
       "        this.ws.send(JSON.stringify({type: \"draw\", figure_id: this.id}));\n",
       "    }\n",
       "}\n",
       "\n",
       "\n",
       "mpl.figure.prototype.handle_save = function(fig, msg) {\n",
       "    var format_dropdown = fig.format_dropdown;\n",
       "    var format = format_dropdown.options[format_dropdown.selectedIndex].value;\n",
       "    fig.ondownload(fig, format);\n",
       "}\n",
       "\n",
       "\n",
       "mpl.figure.prototype.handle_resize = function(fig, msg) {\n",
       "    var size = msg['size'];\n",
       "    if (size[0] != fig.canvas.width || size[1] != fig.canvas.height) {\n",
       "        fig._resize_canvas(size[0], size[1]);\n",
       "        fig.send_message(\"refresh\", {});\n",
       "    };\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_rubberband = function(fig, msg) {\n",
       "    var x0 = msg['x0'] / mpl.ratio;\n",
       "    var y0 = (fig.canvas.height - msg['y0']) / mpl.ratio;\n",
       "    var x1 = msg['x1'] / mpl.ratio;\n",
       "    var y1 = (fig.canvas.height - msg['y1']) / mpl.ratio;\n",
       "    x0 = Math.floor(x0) + 0.5;\n",
       "    y0 = Math.floor(y0) + 0.5;\n",
       "    x1 = Math.floor(x1) + 0.5;\n",
       "    y1 = Math.floor(y1) + 0.5;\n",
       "    var min_x = Math.min(x0, x1);\n",
       "    var min_y = Math.min(y0, y1);\n",
       "    var width = Math.abs(x1 - x0);\n",
       "    var height = Math.abs(y1 - y0);\n",
       "\n",
       "    fig.rubberband_context.clearRect(\n",
       "        0, 0, fig.canvas.width, fig.canvas.height);\n",
       "\n",
       "    fig.rubberband_context.strokeRect(min_x, min_y, width, height);\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_figure_label = function(fig, msg) {\n",
       "    // Updates the figure title.\n",
       "    fig.header.textContent = msg['label'];\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_cursor = function(fig, msg) {\n",
       "    var cursor = msg['cursor'];\n",
       "    switch(cursor)\n",
       "    {\n",
       "    case 0:\n",
       "        cursor = 'pointer';\n",
       "        break;\n",
       "    case 1:\n",
       "        cursor = 'default';\n",
       "        break;\n",
       "    case 2:\n",
       "        cursor = 'crosshair';\n",
       "        break;\n",
       "    case 3:\n",
       "        cursor = 'move';\n",
       "        break;\n",
       "    }\n",
       "    fig.rubberband_canvas.style.cursor = cursor;\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_message = function(fig, msg) {\n",
       "    fig.message.textContent = msg['message'];\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_draw = function(fig, msg) {\n",
       "    // Request the server to send over a new figure.\n",
       "    fig.send_draw_message();\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_image_mode = function(fig, msg) {\n",
       "    fig.image_mode = msg['mode'];\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.updated_canvas_event = function() {\n",
       "    // Called whenever the canvas gets updated.\n",
       "    this.send_message(\"ack\", {});\n",
       "}\n",
       "\n",
       "// A function to construct a web socket function for onmessage handling.\n",
       "// Called in the figure constructor.\n",
       "mpl.figure.prototype._make_on_message_function = function(fig) {\n",
       "    return function socket_on_message(evt) {\n",
       "        if (evt.data instanceof Blob) {\n",
       "            /* FIXME: We get \"Resource interpreted as Image but\n",
       "             * transferred with MIME type text/plain:\" errors on\n",
       "             * Chrome.  But how to set the MIME type?  It doesn't seem\n",
       "             * to be part of the websocket stream */\n",
       "            evt.data.type = \"image/png\";\n",
       "\n",
       "            /* Free the memory for the previous frames */\n",
       "            if (fig.imageObj.src) {\n",
       "                (window.URL || window.webkitURL).revokeObjectURL(\n",
       "                    fig.imageObj.src);\n",
       "            }\n",
       "\n",
       "            fig.imageObj.src = (window.URL || window.webkitURL).createObjectURL(\n",
       "                evt.data);\n",
       "            fig.updated_canvas_event();\n",
       "            fig.waiting = false;\n",
       "            return;\n",
       "        }\n",
       "        else if (typeof evt.data === 'string' && evt.data.slice(0, 21) == \"data:image/png;base64\") {\n",
       "            fig.imageObj.src = evt.data;\n",
       "            fig.updated_canvas_event();\n",
       "            fig.waiting = false;\n",
       "            return;\n",
       "        }\n",
       "\n",
       "        var msg = JSON.parse(evt.data);\n",
       "        var msg_type = msg['type'];\n",
       "\n",
       "        // Call the  \"handle_{type}\" callback, which takes\n",
       "        // the figure and JSON message as its only arguments.\n",
       "        try {\n",
       "            var callback = fig[\"handle_\" + msg_type];\n",
       "        } catch (e) {\n",
       "            console.log(\"No handler for the '\" + msg_type + \"' message type: \", msg);\n",
       "            return;\n",
       "        }\n",
       "\n",
       "        if (callback) {\n",
       "            try {\n",
       "                // console.log(\"Handling '\" + msg_type + \"' message: \", msg);\n",
       "                callback(fig, msg);\n",
       "            } catch (e) {\n",
       "                console.log(\"Exception inside the 'handler_\" + msg_type + \"' callback:\", e, e.stack, msg);\n",
       "            }\n",
       "        }\n",
       "    };\n",
       "}\n",
       "\n",
       "// from http://stackoverflow.com/questions/1114465/getting-mouse-location-in-canvas\n",
       "mpl.findpos = function(e) {\n",
       "    //this section is from http://www.quirksmode.org/js/events_properties.html\n",
       "    var targ;\n",
       "    if (!e)\n",
       "        e = window.event;\n",
       "    if (e.target)\n",
       "        targ = e.target;\n",
       "    else if (e.srcElement)\n",
       "        targ = e.srcElement;\n",
       "    if (targ.nodeType == 3) // defeat Safari bug\n",
       "        targ = targ.parentNode;\n",
       "\n",
       "    // jQuery normalizes the pageX and pageY\n",
       "    // pageX,Y are the mouse positions relative to the document\n",
       "    // offset() returns the position of the element relative to the document\n",
       "    var x = e.pageX - $(targ).offset().left;\n",
       "    var y = e.pageY - $(targ).offset().top;\n",
       "\n",
       "    return {\"x\": x, \"y\": y};\n",
       "};\n",
       "\n",
       "/*\n",
       " * return a copy of an object with only non-object keys\n",
       " * we need this to avoid circular references\n",
       " * http://stackoverflow.com/a/24161582/3208463\n",
       " */\n",
       "function simpleKeys (original) {\n",
       "  return Object.keys(original).reduce(function (obj, key) {\n",
       "    if (typeof original[key] !== 'object')\n",
       "        obj[key] = original[key]\n",
       "    return obj;\n",
       "  }, {});\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.mouse_event = function(event, name) {\n",
       "    var canvas_pos = mpl.findpos(event)\n",
       "\n",
       "    if (name === 'button_press')\n",
       "    {\n",
       "        this.canvas.focus();\n",
       "        this.canvas_div.focus();\n",
       "    }\n",
       "\n",
       "    var x = canvas_pos.x * mpl.ratio;\n",
       "    var y = canvas_pos.y * mpl.ratio;\n",
       "\n",
       "    this.send_message(name, {x: x, y: y, button: event.button,\n",
       "                             step: event.step,\n",
       "                             guiEvent: simpleKeys(event)});\n",
       "\n",
       "    /* This prevents the web browser from automatically changing to\n",
       "     * the text insertion cursor when the button is pressed.  We want\n",
       "     * to control all of the cursor setting manually through the\n",
       "     * 'cursor' event from matplotlib */\n",
       "    event.preventDefault();\n",
       "    return false;\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._key_event_extra = function(event, name) {\n",
       "    // Handle any extra behaviour associated with a key event\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.key_event = function(event, name) {\n",
       "\n",
       "    // Prevent repeat events\n",
       "    if (name == 'key_press')\n",
       "    {\n",
       "        if (event.which === this._key)\n",
       "            return;\n",
       "        else\n",
       "            this._key = event.which;\n",
       "    }\n",
       "    if (name == 'key_release')\n",
       "        this._key = null;\n",
       "\n",
       "    var value = '';\n",
       "    if (event.ctrlKey && event.which != 17)\n",
       "        value += \"ctrl+\";\n",
       "    if (event.altKey && event.which != 18)\n",
       "        value += \"alt+\";\n",
       "    if (event.shiftKey && event.which != 16)\n",
       "        value += \"shift+\";\n",
       "\n",
       "    value += 'k';\n",
       "    value += event.which.toString();\n",
       "\n",
       "    this._key_event_extra(event, name);\n",
       "\n",
       "    this.send_message(name, {key: value,\n",
       "                             guiEvent: simpleKeys(event)});\n",
       "    return false;\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.toolbar_button_onclick = function(name) {\n",
       "    if (name == 'download') {\n",
       "        this.handle_save(this, null);\n",
       "    } else {\n",
       "        this.send_message(\"toolbar_button\", {name: name});\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.toolbar_button_onmouseover = function(tooltip) {\n",
       "    this.message.textContent = tooltip;\n",
       "};\n",
       "mpl.toolbar_items = [[\"Home\", \"Reset original view\", \"fa fa-home icon-home\", \"home\"], [\"Back\", \"Back to  previous view\", \"fa fa-arrow-left icon-arrow-left\", \"back\"], [\"Forward\", \"Forward to next view\", \"fa fa-arrow-right icon-arrow-right\", \"forward\"], [\"\", \"\", \"\", \"\"], [\"Pan\", \"Pan axes with left mouse, zoom with right\", \"fa fa-arrows icon-move\", \"pan\"], [\"Zoom\", \"Zoom to rectangle\", \"fa fa-square-o icon-check-empty\", \"zoom\"], [\"\", \"\", \"\", \"\"], [\"Download\", \"Download plot\", \"fa fa-floppy-o icon-save\", \"download\"]];\n",
       "\n",
       "mpl.extensions = [\"eps\", \"jpeg\", \"pdf\", \"png\", \"ps\", \"raw\", \"svg\", \"tif\"];\n",
       "\n",
       "mpl.default_extension = \"png\";var comm_websocket_adapter = function(comm) {\n",
       "    // Create a \"websocket\"-like object which calls the given IPython comm\n",
       "    // object with the appropriate methods. Currently this is a non binary\n",
       "    // socket, so there is still some room for performance tuning.\n",
       "    var ws = {};\n",
       "\n",
       "    ws.close = function() {\n",
       "        comm.close()\n",
       "    };\n",
       "    ws.send = function(m) {\n",
       "        //console.log('sending', m);\n",
       "        comm.send(m);\n",
       "    };\n",
       "    // Register the callback with on_msg.\n",
       "    comm.on_msg(function(msg) {\n",
       "        //console.log('receiving', msg['content']['data'], msg);\n",
       "        // Pass the mpl event to the overridden (by mpl) onmessage function.\n",
       "        ws.onmessage(msg['content']['data'])\n",
       "    });\n",
       "    return ws;\n",
       "}\n",
       "\n",
       "mpl.mpl_figure_comm = function(comm, msg) {\n",
       "    // This is the function which gets called when the mpl process\n",
       "    // starts-up an IPython Comm through the \"matplotlib\" channel.\n",
       "\n",
       "    var id = msg.content.data.id;\n",
       "    // Get hold of the div created by the display call when the Comm\n",
       "    // socket was opened in Python.\n",
       "    var element = $(\"#\" + id);\n",
       "    var ws_proxy = comm_websocket_adapter(comm)\n",
       "\n",
       "    function ondownload(figure, format) {\n",
       "        window.open(figure.imageObj.src);\n",
       "    }\n",
       "\n",
       "    var fig = new mpl.figure(id, ws_proxy,\n",
       "                           ondownload,\n",
       "                           element.get(0));\n",
       "\n",
       "    // Call onopen now - mpl needs it, as it is assuming we've passed it a real\n",
       "    // web socket which is closed, not our websocket->open comm proxy.\n",
       "    ws_proxy.onopen();\n",
       "\n",
       "    fig.parent_element = element.get(0);\n",
       "    fig.cell_info = mpl.find_output_cell(\"<div id='\" + id + \"'></div>\");\n",
       "    if (!fig.cell_info) {\n",
       "        console.error(\"Failed to find cell for figure\", id, fig);\n",
       "        return;\n",
       "    }\n",
       "\n",
       "    var output_index = fig.cell_info[2]\n",
       "    var cell = fig.cell_info[0];\n",
       "\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_close = function(fig, msg) {\n",
       "    var width = fig.canvas.width/mpl.ratio\n",
       "    fig.root.unbind('remove')\n",
       "\n",
       "    // Update the output cell to use the data from the current canvas.\n",
       "    fig.push_to_output();\n",
       "    var dataURL = fig.canvas.toDataURL();\n",
       "    // Re-enable the keyboard manager in IPython - without this line, in FF,\n",
       "    // the notebook keyboard shortcuts fail.\n",
       "    IPython.keyboard_manager.enable()\n",
       "    $(fig.parent_element).html('<img src=\"' + dataURL + '\" width=\"' + width + '\">');\n",
       "    fig.close_ws(fig, msg);\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.close_ws = function(fig, msg){\n",
       "    fig.send_message('closing', msg);\n",
       "    // fig.ws.close()\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.push_to_output = function(remove_interactive) {\n",
       "    // Turn the data on the canvas into data in the output cell.\n",
       "    var width = this.canvas.width/mpl.ratio\n",
       "    var dataURL = this.canvas.toDataURL();\n",
       "    this.cell_info[1]['text/html'] = '<img src=\"' + dataURL + '\" width=\"' + width + '\">';\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.updated_canvas_event = function() {\n",
       "    // Tell IPython that the notebook contents must change.\n",
       "    IPython.notebook.set_dirty(true);\n",
       "    this.send_message(\"ack\", {});\n",
       "    var fig = this;\n",
       "    // Wait a second, then push the new image to the DOM so\n",
       "    // that it is saved nicely (might be nice to debounce this).\n",
       "    setTimeout(function () { fig.push_to_output() }, 1000);\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._init_toolbar = function() {\n",
       "    var fig = this;\n",
       "\n",
       "    var nav_element = $('<div/>')\n",
       "    nav_element.attr('style', 'width: 100%');\n",
       "    this.root.append(nav_element);\n",
       "\n",
       "    // Define a callback function for later on.\n",
       "    function toolbar_event(event) {\n",
       "        return fig.toolbar_button_onclick(event['data']);\n",
       "    }\n",
       "    function toolbar_mouse_event(event) {\n",
       "        return fig.toolbar_button_onmouseover(event['data']);\n",
       "    }\n",
       "\n",
       "    for(var toolbar_ind in mpl.toolbar_items){\n",
       "        var name = mpl.toolbar_items[toolbar_ind][0];\n",
       "        var tooltip = mpl.toolbar_items[toolbar_ind][1];\n",
       "        var image = mpl.toolbar_items[toolbar_ind][2];\n",
       "        var method_name = mpl.toolbar_items[toolbar_ind][3];\n",
       "\n",
       "        if (!name) { continue; };\n",
       "\n",
       "        var button = $('<button class=\"btn btn-default\" href=\"#\" title=\"' + name + '\"><i class=\"fa ' + image + ' fa-lg\"></i></button>');\n",
       "        button.click(method_name, toolbar_event);\n",
       "        button.mouseover(tooltip, toolbar_mouse_event);\n",
       "        nav_element.append(button);\n",
       "    }\n",
       "\n",
       "    // Add the status bar.\n",
       "    var status_bar = $('<span class=\"mpl-message\" style=\"text-align:right; float: right;\"/>');\n",
       "    nav_element.append(status_bar);\n",
       "    this.message = status_bar[0];\n",
       "\n",
       "    // Add the close button to the window.\n",
       "    var buttongrp = $('<div class=\"btn-group inline pull-right\"></div>');\n",
       "    var button = $('<button class=\"btn btn-mini btn-primary\" href=\"#\" title=\"Stop Interaction\"><i class=\"fa fa-power-off icon-remove icon-large\"></i></button>');\n",
       "    button.click(function (evt) { fig.handle_close(fig, {}); } );\n",
       "    button.mouseover('Stop Interaction', toolbar_mouse_event);\n",
       "    buttongrp.append(button);\n",
       "    var titlebar = this.root.find($('.ui-dialog-titlebar'));\n",
       "    titlebar.prepend(buttongrp);\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._root_extra_style = function(el){\n",
       "    var fig = this\n",
       "    el.on(\"remove\", function(){\n",
       "\tfig.close_ws(fig, {});\n",
       "    });\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._canvas_extra_style = function(el){\n",
       "    // this is important to make the div 'focusable\n",
       "    el.attr('tabindex', 0)\n",
       "    // reach out to IPython and tell the keyboard manager to turn it's self\n",
       "    // off when our div gets focus\n",
       "\n",
       "    // location in version 3\n",
       "    if (IPython.notebook.keyboard_manager) {\n",
       "        IPython.notebook.keyboard_manager.register_events(el);\n",
       "    }\n",
       "    else {\n",
       "        // location in version 2\n",
       "        IPython.keyboard_manager.register_events(el);\n",
       "    }\n",
       "\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._key_event_extra = function(event, name) {\n",
       "    var manager = IPython.notebook.keyboard_manager;\n",
       "    if (!manager)\n",
       "        manager = IPython.keyboard_manager;\n",
       "\n",
       "    // Check for shift+enter\n",
       "    if (event.shiftKey && event.which == 13) {\n",
       "        this.canvas_div.blur();\n",
       "        event.shiftKey = false;\n",
       "        // Send a \"J\" for go to next cell\n",
       "        event.which = 74;\n",
       "        event.keyCode = 74;\n",
       "        manager.command_mode();\n",
       "        manager.handle_keydown(event);\n",
       "    }\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_save = function(fig, msg) {\n",
       "    fig.ondownload(fig, null);\n",
       "}\n",
       "\n",
       "\n",
       "mpl.find_output_cell = function(html_output) {\n",
       "    // Return the cell and output element which can be found *uniquely* in the notebook.\n",
       "    // Note - this is a bit hacky, but it is done because the \"notebook_saving.Notebook\"\n",
       "    // IPython event is triggered only after the cells have been serialised, which for\n",
       "    // our purposes (turning an active figure into a static one), is too late.\n",
       "    var cells = IPython.notebook.get_cells();\n",
       "    var ncells = cells.length;\n",
       "    for (var i=0; i<ncells; i++) {\n",
       "        var cell = cells[i];\n",
       "        if (cell.cell_type === 'code'){\n",
       "            for (var j=0; j<cell.output_area.outputs.length; j++) {\n",
       "                var data = cell.output_area.outputs[j];\n",
       "                if (data.data) {\n",
       "                    // IPython >= 3 moved mimebundle to data attribute of output\n",
       "                    data = data.data;\n",
       "                }\n",
       "                if (data['text/html'] == html_output) {\n",
       "                    return [cell, data, j];\n",
       "                }\n",
       "            }\n",
       "        }\n",
       "    }\n",
       "}\n",
       "\n",
       "// Register the function which deals with the matplotlib target/channel.\n",
       "// The kernel may be null if the page has been refreshed.\n",
       "if (IPython.notebook.kernel != null) {\n",
       "    IPython.notebook.kernel.comm_manager.register_target('matplotlib', mpl.mpl_figure_comm);\n",
       "}\n"
      ],
      "text/plain": [
       "<IPython.core.display.Javascript object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/html": [
       "<img src=\"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAoAAAAHgCAYAAAA10dzkAAAgAElEQVR4Xu2dD6ylRXn/50KDqPWuYsQUuPwRRQVSBbEqRgPRFS21QIwUMRWt0UaMitQKYstvN6WwuomStLpdxYppS2ltAQsUJbZlQbCKVhqEYsEGXAQrbeF2aXWJcH+Zl97LvWfPve8zzzzzzpzzfE5C0rrPPDPzme/M8z3znnPuzMLCwkLgBQEIQAACEIAABCDghsAMBtDNWjNRCEAAAhCAAAQg0BHAACIECEAAAhCAAAQg4IwABtDZgjNdCEAAAhCAAAQggAFEAxCAAAQgAAEIQMAZAQygswVnuhCAAAQgAAEIQAADiAYgAAEIQAACEICAMwIYQGcLznQhAAEIQAACEIAABhANQAACEIAABCAAAWcEMIDOFpzpQgACEIAABCAAAQwgGoAABCAAAQhAAALOCGAAnS0404UABCAAAQhAAAIYQDQAAQhAAAIQgAAEnBHAADpbcKYLAQhAAAIQgAAEMIBoAAIQgAAEIAABCDgjgAF0tuBMFwIQgAAEIAABCGAA0QAEIAABCEAAAhBwRgAD6GzBmS4EIAABCEAAAhDAAKIBCEAAAhCAAAQg4IwABtDZgjNdCEAAAhCAAAQggAFEAxCAAAQgAAEIQMAZAQygswVnuhCAAAQgAAEIQAADiAYgAAEIQAACEICAMwIYQGcLznQhAAEIQAACEIAABhANQAACEIAABCAAAWcEMIDOFpzpQgACEIAABCAAAQwgGoAABCAAAQhAAALOCGAAnS0404UABCAAAQhAAAIYQDQAAQhAAAIQgAAEnBHAADpbcKYLAQhAAAIQgAAEMIBoAAIQgAAEIAABCDgjgAF0tuBMFwIQgAAEIAABCGAA0QAEIAABCEAAAhBwRgAD6GzBmS4EIAABCEAAAhDAAKIBCEAAAhCAAAQg4IwABtDZgjNdCEAAAhCAAAQggAFEAxCAAAQgAAEIQMAZAQygswVnuhCAAAQgAAEIQAADiAYgAAEIQAACEICAMwIYQGcLznQhAAEIQAACEIAABhANQAACEIAABCAAAWcEMIDOFpzpQgACEIAABCAAAQwgGoAABCAAAQhAAALOCGAAnS0404UABCAAAQhAAAIYQDQAAQhAAAIQgAAEnBHAADpbcKYLAQhAAAIQgAAEMIBoAAIQgAAEIAABCDgjgAF0tuBMFwIQgAAEIAABCGAA0QAEIAABCEAAAhBwRgAD6GzBmS4EIAABCEAAAhDAAKIBCEAAAhCAAAQg4IwABtDZgjNdCEAAAhCAAAQggAFEAxCAAAQgAAEIQMAZAQygswVnuhCAAAQgAAEIQAADiAYgAAEIQAACEICAMwIYQGcLznQhAAEIQAACEIAABhANQAACEIAABCAAAWcEMIAZC/7YY4+F++67LzztaU8LMzMzGZloCgEIQAACEIDAUAQWFhbCjh07wj777BN22223obptqh8MYMZy3HvvvWFubi4jA00hAAEIQAACEKhFYPv27WG//far1X3VfjGAGfjn5+fD05/+9BAFNDs7m5GJphCAAAQgAAEIDEXgv//7v7sLnIceeiisW7duqG6b6gcDmLEcUUBRONEIYgAzQNIUAhCAAAQgMCAB6ncIGMAMwSGgDHg0hQAEIAABCFQiQP3GAGZJDwFl4aMxBCAAAQhAoAoB6jcGMEt4CCgLH40hAAEIQAACVQhQvzGAWcJDQFn4aAwBCEAAAhCoQoD6jQHMEh4CysJHYwhAAAIQgEAVAtRvDGCW8BBQFj4aQwACEIAABKoQoH5jALOEh4Cy8NEYAhCAAAQgUIUA9RsDmCU8BJSFj8YQgAAEIACBKgSo3xjALOEhoCx8NIYABCAAAQhUIUD9dm4ADzzwwHDPPffsIr7TTz89fOpTn+oVJQLqRUQABCAAAQhAoDkC1G/nBvCBBx4Ijz766JIwv/vd74b169eHf/iHfwjHHHNMr2ARUC8iAiAAAQhAAALNEaB+OzeAo4o844wzwlVXXRXuvPPOMDMz0ytYBNSLiAAIQAACEIBAcwSo3xjAJVE+8sgjYZ999glnnnlmOOecc0RiRUAiTARBAAIQgAAEmiJA/cYALgnyL//yL8Opp54afvCDH3RGcNxr586dIf63+IoCmpubC/Pz82F2drYpcTMYCEAAAhCAAATGE8AAYgCXlHHccceFPfbYI1x55ZWr7pcNGzaEjRs37vLvXgzggWdfvTT3uzcdz7kCAQhAAAIQmEgCGEAMYCfc+E3g5zznOeGyyy4LJ5xwwqpi9n4DiAGcyHOOQUMAAhCAwAgBDCAGsJNEvNnbunVr2L59e/i5n/s58UbxJiAMoFgaBEIAAhCAQMMEvNXvcUsxs7CwsNDwGhUf2mOPPRYOOuig8Ja3vCVs2rQpqT9vAsIAJsmDYAhAAAIQaJSAt/qNARxD4Nprrw3x83/f+973wiGHHJIkVW8CwgAmyYNgCEAAAhBolIC3+o0BNBaiNwFhAI0FRDoIQAACEKhCwFv9xgAay8ybgDCAxgIiHQQgAAEIVCHgrX5jAI1l5k1AGEBjAZEOAhCAAASqEPBWvzGAxjLzJiAMoLGASAcBCEAAAlUIeKvfGEBjmXkTEAbQWECkgwAEIACBKgS81W8MoLHMvAkIA2gsINJBAAIQgEAVAt7qNwbQWGbeBIQBNBYQ6SAAAQhAoAoBb/UbA2gsM28CwgAaC4h0EIAABCBQhYC3+o0BNJaZNwFhAI0FRDoIQAACEKhCwFv9xgAay8ybgDCAxgIiHQQgAAEIVCHgrX5jAI1l5k1AGEBjAZEOAhCAAASqEPBWvzGAxjLzJiAMoLGASAcBCEAAAlUIeKvfGEBjmXkTEAbQWECkgwAEIACBKgS81W8MoLHMvAkIA2gsINJBAAIQgEAVAt7qNwbQWGbeBIQBNBYQ6SAAAQhAoAoBb/UbA2gsM28CwgAaC4h0EIAABCBQhYC3+o0BNJaZNwFhAI0FRDoIQAACEKhCwFv9xgAay8ybgDCAxgIiHQQgAAEIVCHgrX5jAI1l5k1AGEBjAZEOAhCAAASqEPBWvzGAxjLzJiAMoLGASAcBCEAAAlUIeKvfGEBjmXkTEAbQWECkgwAEIACBKgS81W8MoLHMvAkIA2gsINJBAAIQgEAVAt7qNwbQWGbeBIQBNBYQ6SAAAQhAoAoBb/UbA2gsM28CwgAaC4h0EIAABCBQhYC3+o0BNJaZNwFhAI0FRDoIQAACEKhCwFv9xgAay8ybgDCAxgIiHQQgAAEIVCHgrX5jAI1l5k1AGEBjAZEOAhCAAASqEPBWvzGAxjLzJiAMoLGASAcBCEAAAlUIeKvfGEBjmXkTEAbQWECkgwAEIACBKgS81W8MoLHMvAkIA2gsINJBAAIQgEAVAt7qNwbQWGbeBIQBNBYQ6SAAAQhAoAoBb/UbA2gsM28CwgAaC4h0EIAABCBQhYC3+o0BNJaZNwFhAI0FRDoIQAACEKhCwFv9xgAay8ybgDCAxgIiHQQgAAEIVCHgrX5jAI1l5k1AGEBjAZEOAhCAAASqEPBWvzGAxjLzJiAMoLGASAcBCEAAAlUIeKvfGEBjmXkTEAbQWECkgwAEIACBKgS81W8MoLHMvAkIA2gsINJBAAIQgEAVAt7qNwbQWGbeBIQBNBYQ6SAAAQhAoAoBb/UbA2gsM28CwgAaC4h0EIAABCBQhYC3+o0BHCHwwx/+MJx11lnhmmuuCT/5yU/CIYccEj73uc+Fl7zkJSJBehMQBlAkC4IgAAEIQKBxAt7qNwZwGYEHH3wwHHHEEeHYY48N73nPe8Lee+8dvv/974cDDzwwHHzwwSLpehMQBlAkC4IgAAEIQKBxAt7qNwZwGYGzzz473HjjjeGGG25Qy9SbgDCAaqnQEAIQgAAEGiLgrX5jAJcROPTQQ8Nxxx0X7r333rBt27aw7777htNPPz28613vEkvUm4AwgGJpEAgBCEAAAg0T8Fa/MYDLCOy5557d/3fmmWeGN7/5zeGb3/xmOOOMM8LWrVvD2972trGy3blzZ4j/Lb6igObm5sL8/HyYnZ1tWOo2Q8MA2nAkCwQgAAEI1CWAAQxhZmFhYaHuMtTpfY899ghHHXVUuOmmm5YG8P73vz/cfPPN4etf//rYQW3YsCFs3Lhxl3/DANZZQ3qFAAQgAAEIaAhgAB0bwAMOOCCsX78+XHTRRUva2bJlSzjvvPNC/HbwuBc3gFcvYbl70/GaPUcbCEAAAhCAQHUCGEDHBvDUU08N27dvX/ElkA9+8IPhG9/4xopbwbVU6k1APAKufmYxAAhAAAIQMCDgrX6PQ+b2EXB81Hv00Ud3j3RPPvnk7jOA8Qsgn/nMZ8Jb3/pWkby8CQgDKJIFQRCAAAQg0DgBb/UbAzhC4Kqrrgof+chHwp133hkOOuig7gshfAt49V2LAWz8RGN4EIAABCAgIoABdPwIWKSQniBvAsIAWqiGHBCAAAQgUJuAt/rNDaCx4rwJCANoLCDSQQACEIBAFQLe6jcG0Fhm3gSEATQWEOkgAAEIQKAKAW/1GwNoLDNvAsIAGguIdBCAAAQgUIWAt/qNATSWmTcBYQCNBUQ6CEAAAhCoQsBb/cYAGsvMm4AwgMYCIh0EIAABCFQh4K1+YwCNZeZNQBhAYwGRDgIQgAAEqhDwVr8xgMYy8yYgDKCxgEgHAQhAAAJVCHir3xhAY5l5ExAG0FhApIMABCAAgSoEvNVvDKCxzLwJCANoLCDSQQACEIBAFQLe6jcG0Fhm3gSEATQWEOkgAAEIQKAKAW/1GwNoLDNvAsIAGguIdBCAAAQgUIWAt/qNATSWmTcBYQCNBUQ6CEAAAhCoQsBb/cYAGsvMm4AwgMYCIh0EIAABCFQh4K1+YwCNZeZNQBhAYwGRDgIQgAAEqhDwVr8xgMYy8yYgDKCxgEgHAQhAAAJVCHir3xhAY5l5ExAG0FhApIMABCAAgSoEvNVvDKCxzLwJCANoLCDSQQACEIBAFQLe6jcG0Fhm3gSEATQWEOkgAAEIQKAKAW/1GwNoLDNvAsIAGguIdBCAAAQgUIWAt/qNATSWmTcBYQCNBUQ6CEAAAhCoQsBb/cYAGsvMm4AwgMYCIh0EIAABCFQh4K1+YwCNZeZNQBhAYwGRDgIQgAAEqhDwVr8xgMYy8yYgDKCxgEgHAQhAAAJVCHir3xhAY5l5ExAG0FhApIMABCAAgSoEvNVvDKCxzLwJCANoLCDSQQACEIBAFQLe6jcG0Fhm3gSEATQWEOkgAAEIQKAKAW/1GwNoLDNvAsIAGguIdBCAAAQgUIWAt/qNATSWmTcBYQCNBUQ6CEAAAhCoQsBb/cYAGsvMm4AwgMYCIh0EIAABCFQh4K1+YwCNZeZNQBhAYwGRDgIQgAAEqhDwVr8xgMYy8yYgDKCxgEgHAQhAAAJVCHir3xhAY5l5ExAG0FhApIMABCAAgSoEvNVvDKCxzLwJCANoLCDSQQACEIBAFQLe6jcG0Fhm3gSEATQWEOkgAAEIQKAKAW/1GwNoLDNvAsIAGguIdBCAAAQgUIWAt/qNATSWmTcBYQCNBUQ6CEAAAhCoQsBb/cYAGsvMm4AwgMYCIh0EIAABCFQh4K1+YwCNZeZNQBhAYwGRDgIQgAAEqhDwVr8xgMYy8yYgDKCxgEgHAQhAAAJVCHir3xhAY5l5ExAG0FhApIMABCAAgSoEvNVvDOAyAhs2bAgbN25cweTZz352+NGPfiQWozcBYQDF0iAQAhCAAAQaJuCtfmMARwzgX/3VX4WvfvWrS//r7rvvHp71rGeJJetNQBhAsTQIhAAEIACBhgl4q98YwBEDeMUVV4RbbrlFLVFvAsIAqqVCQwhAAAIQaIiAt/qNARwxgJs3bw7r1q0LT3rSk8LLXvaycP7554fnPOc5q0p0586dIf63+IoCmpubC/Pz82F2drYhaZcZCgawDFeyQgACEIDAsAQwgCHMLCwsLAyLvY3errnmmvC///u/4ZBDDgn//u//Hs4777xwxx13hNtuuy0885nPHDvIcZ8bjIEYwDbWlFFAAAIQgAAEJAQwgI4N4KhA/ud//iccfPDB4cMf/nA488wzx+qHG8Crl7jcvel4yR4jBgIQgAAEINAcAQwgBnCFKNevXx+e+9znhi1btojE6k1APAIWyYIgCEAAAhBonIC3+j1uOdw+Ah6FEW/34g3gu9/97nDuueeKpOtNQBhAkSwIggAEIACBxgl4q98YwGUEPvShD4U3vvGNYf/99w8//vGPu88Abtu2Ldx6663hgAMOEEnXm4AwgCJZEAQBCEAAAo0T8Fa/MYDLCJxyyinh+uuvD//xH//R/fbfy1/+8vB7v/d74dBDDxXL1puAMIBiaRAIAQhAAAINE/BWvzGAxmL0JiAMoLGASAcBCEAAAlUIeKvfGEBjmXkTEAbQWECkgwAEIACBKgS81W8MoLHMvAkIA2gsINJBAAIQgEAVAt7qNwbQWGbeBIQBNBYQ6SAAAQhAoAoBb/UbA2gsM28CwgAaC4h0EIAABCBQhYC3+o0BNJaZNwFhAI0FRDoIQAACEKhCwFv9xgAay8ybgDCAxgIiHQQgAAEIVCHgrX5jAI1l5k1AGEBjAZEOAhCAAASqEPBWvzGAxjLzJiAMoLGASAcBCEAAAlUIeKvfGEBjmXkTEAbQWECkgwAEIACBKgS81W8MoLHMvAkIA2gsINJBAAIQgEAVAt7qNwbQWGbeBIQBNBYQ6SAAAQhAoAoBb/UbA2gsM28CwgAaC4h0EIAABCBQhYC3+o0BNJaZNwFhAI0FRDoIQAACEKhCwFv9xgAay8ybgDCAxgIiHQQgAAEIVCHgrX5jAI1l5k1AGEBjAZEOAhCAAASqEPBWvzGAxjLzJiAMoLGASAcBCEAAAlUIeKvfGEBjmXkTEAbQWECkgwAEIACBKgS81W8MoLHMvAkIA2gsINJBAAIQgEAVAt7qNwbQWGbeBIQBNBYQ6SAAAQhAoAoBb/UbA2gsM28CwgAaC4h0EIAABCBQhYC3+o0BNJaZNwFhAI0FRDoIQAACEKhCwFv9xgAay8ybgDCAxgIiHQQgAAEIVCHgrX5jAI1l5k1AGEBjAZEOAhCAAASqEPBWvzGAxjLzJiAMoLGASAcBCEAAAlUIeKvfGEBjmXkTEAbQWECkgwAEIACBKgS81W8MoLHMvAkIA2gsINJBAAIQgEAVAt7qNwbQWGbeBIQBNBYQ6SAAAQhAoAoBb/UbA2gsM28CwgAaC4h0EIAABCBQhYC3+o0BNJaZNwFhAI0FRDoIQAACEKhCwFv9xgAay8ybgDCAxgIiHQQgAAEIVCHgrX5jAI1l5k1AGEBjAZEOAhCAAASqEPBWvzGAxjLzJiAMoLGASAcBCEAAAlUIeKvfGEBjmXkTEAbQWECkgwAEIACBKgS81W8MoLHMvAkIA2gsINJBAAIQgEAVAt7qNwbQWGbeBIQBNBYQ6SAAAQhAoAoBb/UbA2gsM28CwgAaC4h0EIAABCBQhYC3+o0BNJaZNwFhAI0FRDoIQAACEKhCwFv9xgAay8ybgDCAxgIiHQQgAAEIVCHgrX5jAI1l5k1AGEBjAZEOAhCAAASqEPBWvzGAa8jsggsuCOecc074wAc+EC688EKRIL0JCAMokgVBEIAABCDQOAFv9RsDuIogb7755nDyySeH2dnZcOyxx2IAV+GEAWz8RGN4EIAABCAgIoABDGFmYWFhQURrSoMefvjhcOSRR4ZPf/rT4bzzzgsvfvGLMYAYwClVO9OCAAQgAIFIAAOIAQynnXZa2GuvvcInP/nJcMwxx2AA1zgbuAHk4IQABCAAgWkggAF0bgAvvfTS7tbvW9/6Vthzzz17DeDOnTtD/G/xFQU0NzcX5ufnu8fH0/7CAE77CjM/CEAAAj4IYAAdG8Dt27eHo446Klx77bXhRS96Uaf4vhvADRs2hI0bN+6yOzCAPg4MZgkBCEAAAtNBAAPo2ABeccUV4aSTTgq77777kpofffTRMDMzE3bbbbfupm/5v8UgbgCvXmJ196bjp+MUYBYQgAAEIOCOAAbQsQHcsWNHuOeee1aI/h3veEd4wQteEM4666xw+OGH924IbwLiEXCvJAiAAAQgAIEJIOCtfo9bEvffAl4Ope8R8ChAbwLCAE7AqcYQIQABCECgl4C3+o0B7JEEBnBtQBjA3jOFAAhAAAIQmAACGEDHj4At9OlNQBhAC9WQAwIQgAAEahPwVr+5ATRWnDcBYQCNBUQ6CEAAAhCoQsBb/cYAGsvMm4AwgMYCIh0EIAABCFQh4K1+YwCNZeZNQBhAYwGRDgIQgAAEqhDwVr8xgMYy8yYgDKCxgEgHAQhAAAJVCHir3xhAY5l5ExAG0FhApIMABCAAgSoEvNVvDKCxzLwJCANoLCDSQQACEIBAFQLe6jcG0Fhm3gSEATQWEOkgAAEIQKAKAW/1GwNoLDNvAsIAGguIdBCAAAQgUIWAt/qNATSWmTcBYQCNBUQ6CEAAAhCoQsBb/cYAGsvMm4AwgMYCIh0EIAABCFQh4K1+YwCNZeZNQBhAYwGRDgIQgAAEqhDwVr8xgMYy8yYgDKCxgEgHAQhAAAJVCHir3xhAY5l5ExAG0FhApIMABCAAgSoEvNVvDKCxzLwJCANoLCDSQQACEIBAFQLe6jcG0Fhm3gSEATQWEOkgAAEIQKAKAW/1GwNoLDNvAsIAGguIdBCAAAQgUIWAt/qNATSWmTcBYQCNBUQ6CEAAAhCoQsBb/cYAGsvMm4AwgMYCIh0EIAABCFQh4K1+YwCNZeZNQBhAYwGRDgIQgAAEqhDwVr8xgMYy8yYgDKCxgEgHAQhAAAJVCHir3xhAY5l5ExAG0FhApIMABCAAgSoEvNVvDKCxzLwJCANoLCDSQQACEIBAFQLe6jcG0Fhm3gSEATQWEOkgAAEIQKAKAW/1GwNoLDNvAsIAGguIdBCAAAQgUIWAt/qNATSWmTcBYQCNBUQ6CEAAAhCoQsBb/cYAGsvMm4AwgMYCIh0EIAABCFQh4K1+YwCNZeZNQBhAYwGRDgIQgAAEqhDwVr8xgMYy8yYgDKCxgEgHAQhAAAJVCHir3xhAY5l5ExAG0FhApIMABCAAgSoEvNVvDKCxzLwJCANoLCDSQQACEIBAFQLe6jcG0Fhm3gSEATQWEOkgAAEIQKAKAW/1GwNoLDNvAsIAGguIdBCAAAQgUIWAt/qNATSWmTcBYQCNBUQ6CEAAAhCoQsBb/cYAGsvMm4AwgMYCIh0EIAABCFQh4K1+YwCNZeZNQBhAYwGRDgIQgAAEqhDwVr8xgMYy8yYgDKCxgEgHAQhAAAJVCHir3xhAY5l5ExAG0FhApIMABCAAgSoEvNVvDKCxzLwJCANoLCDSQQACEIBAFQLe6jcG0Fhm3gSEATQWEOkgAAEIQKAKAW/1GwNoLDNvAsIAGguIdBCAAAQgUIWAt/qNARwhsGXLlhD/u/vuu7t/Oeyww8K5554b3vCGN4gE6U1AGECRLAiCAAQgAIHGCXir3xjAEQJXXnll2H333cNzn/vc7l++8IUvhM2bN4fvfOc7nRnse3kTEAawTxH8OwQgAAEITAIBb/UbAyhQ5V577dWZwHe+85290d4EhAHslQQBEHBHgHPB3ZJPxYS91W8M4BqyffTRR8MXv/jFcNppp3U3gIceeugu0Tt37gzxv8VXFNDc3FyYn58Ps7OzU7Ep1poEB/3ULzEThEAyAc6FZGQ0aIAABjCEmYWFhYUG1qLaEG699dbwile8Ivz0pz8NP//zPx8uueSS8Mu//Mtjx7Nhw4awcePGXf4NA1ht+egYAhCoTAADWHkB6F5FAAOIAQyPPPJI+MEPfhAeeuih8Nd//dfhoosuCtu2beMGcMyW4qBXnTM0gsBUE+BcmOrlndrJYQAxgLuI+7WvfW04+OCDw9atW3uF701AHPS9kiAAAu4IcC64W/KpmLC3+j1u0dw/Ah6F8prXvKb7XN/FF1/cK3JvAuKg75UEARBwR4Bzwd2ST8WEvdVvDOAIgXPOOaf7zb9o+Hbs2BEuvfTSsGnTpvDlL385rF+/vlfk3gTEQd8rCQIg4I4A54K7JZ+KCXur3xjAEQLxp17+7u/+Ltx///1h3bp14Rd/8RfDWWedJTJ/MZU3AXHQT8W5xyQgYEqAc8EUJ8kGIuCtfmMAjYXlTUAc9MYCIh0EpoAA58IULOKAU2hFL97qNwbQWOTeBNTKxjVeRtJBAAIZBDgXMuA5bNqKXrzVbwyg8WbzJqBWNq7xMpIOAhDIIMC5kAHPYdNW9OKtfmMAjTebNwG1snGNl5F0EIBABgHOhQx4Dpu2ohdv9RsDaLzZvAmo1MZdnjcu0d2bjjdeKdJBAAKlCJQ6F0qNl7x1CbSiF2/1GwNorHtvAiq1cTGAxsIkHQQGJFDqXBhwCnQ1IIFW9OKtfmMAjUXuTUClNi4G0FiYpIPAgARKnQsDToGuBiTQil681W8MoLHIvQmo1MbFABoLk3QQGJBAqXNhwCnQ1YAEWtGLt/qNATQWuTcBldq4GEBjYZIOAgMSKHUuDDgFuhqQQCt68Va/MYDGIvcmoFIbFwNoLEzSQWBAAqXOhQGnQFcDEmhFL97qNwbQWOTeBFRq42IAjYVJOggMSKDUuTDgFOhqQAKt6MVb/cYAGovcm4BKbVwMoLEwSQeBAQmUOhcGnAJdDUigFb14q98YQGORexNQqY2LATQWJukgMCCBUufCgFOgqwEJtKIXb/UbA2gscm8CKrVxMYDGwiQdBAYkUOpcGHAKdDUggVb04q1+YwCNRe5NQKU2LgbQWJikg8CABEqdCwNOga4GJNCKXrzVbwygsci9CajUxsUAGguTdBAYkECpc2HAKdDVgARa0Yu3+o0BNBa5NwGV2rgYwH5hlmLf3zMREFibANpEISkEWtGLt/qNAUxRqSDWm4BKbVwMYL/YSrHv7+lLVb8AACAASURBVJkICGAA0YAdgVbOMm/1GwNop+EukzcBldq4GMB+YZZi398zERDAAKIBOwKtnGXe6jcG0E7DGMBNx5vRxAD2o2zl0OwfKRHeCKBNbyueN99W9IIBDGFmYWFhIW85/bb2JqBSGxcD2L+HSrHv75kICHADiAbsCLRylnmr39wA2mmYG0BuAI3VRJEdFCidmRFopaCbTYhERQm0ohcMIDeAWUL3JqBSG5cbwH4ZlmLf3zMREODNCRqwI9DKWeatfnMDaKdhbgC5ATRWE0V2UKCrdMabkfRVaKWgp4+cFjUItKIXDCA3gFn69yagUhuXotsvw1Ls+3v2FYEW09cbbaYz89yiFb14q9/cABrvOm8CKrVxKbr9wizFvr9nXxFoMX290WY6M88tWtGLt/qNATTedd4EVGrjUnT7hVmKfX/PviLQYvp6o810Zp5btKIXb/UbA2i864YWUO3iVGrj1p6XsSyKpCvFvshgJzgpWkxfPLSZzsxzi1b0MnT9bnHN+R3AjFUZWkC1i1OpjVt7XhkSGKxpKfaDTWBCOkKL6QuFNtOZeW7Ril6Grt8trjkGMGNVhhZQ7eJUauPWnleGBAZrWor9YBOYkI7QYvpCoc10Zp5btKKXoet3i2uOAcxYlaEFVLs4ldq4teeVIYHBmpZiP9gEJqQjtJi+UGgznZnnFq3oZej63eKaYwAzVqW0gEY3Su3iVGrj1p5XhgQGa1qK/WATmJCO0GL6QqHNdGaeW7Sil9L1exLWGAOYsUqlBYQBzFicKWvayqE5ZVh3mQ4GMH2Fp0GbrHv6umtbtKKX0vVby2fIdhjADNqlBYQBzFicKWvayqE5ZVgxgAYLOg3axAAaCEGYohW9lK7fQhxVwzCAGfhLCwgDmLE4U9a0lUNzyrCqDKCVWbDKU3tNpkGb07IWk6CFVvRSun7XXgtJ/xhACaVVYkoLCAOYsThT1rSVQ3PKsGIADRZ0GrSJATQQQghBogVJjM1o1s5Sun4PMYfcPjCAGQRLCwgDmLE4U9a0lUNzyrBiAA0WdBq0iQE0EAIG0AbigFkwgBmwMYAZ8JY15fDt5zgNRbZ/lvUjJFqUxEhmYpVH0lfJmGnQ5rSsRcl1luSWaEESI+krN6Z0/c4d3xDtMYAZlEsLiBvAjMWZsqatHJpThpUbQIMFnQZtYgANhMANoA3EAbNgADNgYwB18FoztrpZDNtqGorssMR0vUmMgCRG0rtVHklfJWOmQZvTshYl11mSW6IFSYykr9yY0vU7d3xDtHdtAC+44IJw2WWXhTvuuCM8+clPDkcffXT42Mc+Fp7//OeL2JcWUGtGyWrjtjYv0WJXDrJiX3kazXcvMQKSGMlErfJI+hoXY6UpqzzaeVi0q70WFnNoIYdEC5KYIeZSun4PMYfcPlwbwNe//vXhlFNOCS996UvDz372s/DRj3403HrrreH2228PT33qU3vZlhZQa0bJauO2Nq/ehW4gwIp9A1NpeggSIyCJkUzSKo+kLwzg2pRqr4V2DVtrJzmnJDFDzKt0/R5iDrl9uDaAo/AeeOCBsPfee4dt27aFV7/61b1sSwuoNaNktXFbm1fvQjcQYMW+gak0PQSJEZDESCZplUfSFwawLQM4rftZMi9JjFbTKe1K1++UsdSKxQAuI3/XXXeF5z3ved0t4OGHH77LmuzcuTPE/xZfUUBzc3Nhfn4+zM7Omq9ha0bJauO2Ni/zhSuQ0Ip9gaFNVUqJKZPESKBY5ZH0hQHEAGp1ktJOck5JYlL61MZiAEPAAP6fehYWFsIJJ5wQHnzwwXDDDTeM1dSGDRvCxo0bd/k3DGDaFsQApvGK0a0cmukjn6wWElMmiZHM2iqPpC8MIAZQq5OUdpJzShKT0qc2FgOIAVzSznvf+95w9dVXh6997Wthv/32G6spbgCvXuJy96bjtftuFzNTuxCqJzJgw1YOzQGnXKUriRYlMZLBW+WR9IUBxABqdZLSTnJOSWJS+tTGYgAxgJ123ve+94UrrrgiXH/99eGggw4S66m0gFq7KbPauK3NS7zgFQOt2FecwkR0LTFlkhjJZK3ySPrCAGIAtTpJaSc5pyQxKX1qY0vXb+24hmzn+hFwfOwbzd/ll18errvuuu7zfymv0gJqzShZbdzW5pWy5rVirdjXGv+k9CsxZZIYyXyt8kj6wgBiALU6SWknOackMSl9amNL12/tuIZs59oAnn766eGSSy4JX/rSl1b89t+6deu63wXse5UWUGtGyWrjtjavvnVu4d+t2Lcwl5bHIDFlkhjJHK3ySPrCAGIAtTpJaSc5pyQxKX1qY0vXb+24hmzn2gDOzMyMZf35z38+vP3tb+9dh9ICas0oWW3c1ubVu9ANBFixb2AqTQ9BYsokMZJJWuWR9IUBxABqdZLSTnJOSWJS+tTGlq7f2nEN2c61AcwFXVpArRklq43b2rxydTBEeyv2Q4x1kvuQmDJJjISBVR5JXxhADKBWJyntJOeUJCalT21s6fqtHdeQ7TCAGbRLC6g1o2S1cVubV4YEBmtqxX6wAU9oRxJTJomRTN8qj6QvDCAGUKuTlHaSc0oSk9KnNrZ0/daOa8h2GMAM2qUF1JpRstq4rc0rQwKDNbViP9iAJ7QjiSmTxEimb5VH0hcGcDoNYGvngmQ8khitplPala7fKWOpFYsBzCBfWkCtGSWrjdvavDIkMFhTK/aDDXhCO5KYMkmMZPpWeSR9YQAxgFqdpLSTnFOSmJQ+tbGl67d2XEO2wwBm0C4toNaMktXGbW1eGRIYrKkV+8EGPKEdSUyZJEYyfas8kr4wgBhArU5S2knOKUlMSp/a2NL1WzuuIdthADNolxZQa0bJauO2Nq8MCQzW1Ir9YAOe0I4kpkwSI5m+VR5JXxhADKBWJyntJOeUJCalT21s6fqtHdeQ7TCAGbRLC6g1o2S1cVubV4YEBmtqxV4z4NpGRTNmbRvJXCUxkv6t8kj6wgBiALU6SWknOackMSl9amNL12/tuIZshwHMoF1aQK0ZJauN29q8MiQwWFMr9poB1zYqmjFr20jmKomR9G+VR9IXBhADqNVJSjvJOSWJSelTG1u6fmvHNWQ7DGAG7dICas0oWW3c1uaVIYHBmlqx1wy4tlHRjFnbRjJXSYykf6s8kr4wgCsJjLIf5XP3puO1WEXtrPazVR7RoAVBkvFIYgRdZYeUrt/ZAxwgAQYwA3JpAbVmlKw2bmvzypDAYE2t2GsGXNuoaMasbSOZqyRG0r9VHklfGEAMoFYnKe0k55QkJqVPbWzp+q0d15DtMIAZtEsLqDWjZLVxW5tXhgQGa2rFXjPg2kZFM2ZtG8lcJTGS/q3ySPrCAGIAtTpJaSc5pyQxKX1qY0vXb+24hmyHAcygXVpArRklq43b2rwyJDBYUyv2mgHXNiqaMWvbSOYqiZH0b5VH0hcGEAOo1UlKO8k5JYlJ6VMbW7p+a8c1ZDsMYAbt0gJqzShZbdzW5pUhgcGaWrHXDLi2UdGMWdtGMldJjKR/qzySvjCAGECtTlLaSc4pSUxKn9rY0vVbO64h22EAM2iXFlBrRslq47Y2rwwJDNbUir1mwLWNimbM2jaSuUpiJP1b5ZH0hQHEAGp1ktJOck5JYlL61MaWrt/acQ3ZDgOYQbu0gFozSlYbt7V5ZUhgsKZW7DUDrm1UNGPWtpHMVRIj6d8qj6QvDCAGUKuT1dqN06/knJLEWI91XL7S9XuIOeT2gQHMIFhaQK0ZJauN29q8MiQwWFMr9poB1zYqmjFr20jmKomR9G+VR9IXBhADqNUJBtCaXDv5MIAZa4EB1MHDAKZzwwCmM9O0kJgySYykb6s8kr4wgBhArU4wgNbk2smHAcxYCwygDh4GMJ0bBjCdmaaFxJRJYiR9W+WR9IUBxABqdYIBtCbXTj4MYMZaYAB18DCA6dxaM4A1x5NOT95CYsokMZIerfJI+sIAYgC1OsEAWpNrJx8GMGMtMIA6eBjAdG41DZf2w97ps6zfQmLKJDGSmVjlkfSFAcQAanWCAbQm104+DGDGWmAAdfAwgOncMIDpzDQtJKZMEiPp2yqPpC8M4HAGULJXJTGSdbXKI+lrNEb7xrDmmJfPoXT91jAdug0GMIN4aQG1ZpSsNm5r88qQwGBNrdhrBqw96DV91W4jMWWSGMk8rPJI+sIAYgC1OuEG0JpcO/kwgBlrgQHUwcMApnPDAKYz07SQmDJJjKRvqzySvjCAGECtTjCA1uTayYcBzFgLDKAOHgYwnRsGMJ2ZpoXElEliJH1b5ZH0hQEsYwC1t+NW+9kqz6g+JNqsPXet7hfbla7fueMboj0GMINyaQG1ZpSsDpvW5pUhgcGaWrHXDFh70Gv6qt1GW/g045b0pckrbWOlKas80nFbxI2yH81596bjRd1o94YVM6s8GMBZ0XpPWxAGMGNFMYA6eBjAdG6lDnrJSLRFTpK7tRiJKZPESOZllUfSFzeA3ACm6ESiTe25UPMsW86gdP1O4V0rFgOYQb60gFozSlYbt7V5ZUhgsKZW7DUD1h70mr5qt9EWPs24JX1p8krbWGnKKo903BZx3ACuTVGiTe250IpeStdvC52WzoEBzCBcWkCtGSWrjdvavDIkMFhTK/aaAWsPek1ftdtoC59m3JK+NHmlbaw0ZZVHOm6LuEkwgBKukhgNL4k2tedCqTGnzrN0/U4dT414DGAG9dICas0oaTau5JCQHDYZy1StqeW8NOytJi5ZQ6u+aueRrJkkRjIPqzySvsbFaDQ1LVrAAHIDWLp+a/flkO0wgBm0SwsIA5ixOA00tSzwmmJthWBair6Eh2TNJDFWfUnyaGM0mpoWLWAAMYCl67d2Xw7ZDgOYQbu0gDCAGYvTQFMroxCnoinWVgimpehLeEjWTBJj1ZckjzZGo6lp0QIGEANYun5r9+WQ7TCAGbRLC8izAdQUp4ylLNLUyihgAIssz9ikkjWTxEhGbJVH0te4GM0ewwCuJKnlIWFvFaPRh0SbJeeuGXNqm9L1O3U8NeIxgBnUSwsIA/j44kh/kytjKYs0lRyi0o4lxUCaKzVOe9Cn9tNCvGTNJDGSuVjlkfSFAVzbuI3ykZ452r0h2c9WMRp9SLRZcu6aMae2KV2/U8dTIx4DmEG9tIAwgBjARXlKikGGlNdsqj3oS42nZF5t4dOMqTZXjaZqj1nDuc/8jvt3DODVK7CM46HVgkZ3Vuu+PE/p+l1izNY5MYAZREsLCAOIAcQAZmxQRdNJMIBWBVSTR1v0FUtRtMnoPLgBXPuGFANYVI7VkmMAM9BjAPvhSQqGJKa/p/YiJGZCOmpNsZbm7oub1vWR3AxJC18fQ2lfknWWxEjGo8kzLVrAAK6tEMnZpdWCRncSPafGlK7fqeOpEY8BzKBeWkDcAHIDyA1gxgZVNNUWPkVXK77ZHdtHsykpjpIYyXg0ebRFXzIeqxjNGnIDyA2glf4mKQ8GMGO1MID98CQFQxLT31N7EZJCJB21plhLc/fFTev6SG/lRuOs1lXL1UoLmjzaMfdpzPLfJevDDSA3gKXrt6WmS+XCAGaQLS0gbgC5AeQGMGODKppqzIP0CwMSIykxZZIYydQ1eTCA/TdlEq5Dxki0INGmJMZqXpoxp7YpXb9Tx1MjHgOYQb20gDCAGEAMYMYGVTTFAKbfDEmKvmIp1E00azjamdTUaw2xhJlVjAakhuHQH2HQzGt5m9L1O3d8Q7R3bQCvv/76sHnz5vDtb3873H///eHyyy8PJ554oph7aQFhADGAGEDxdjQJ1BY+TeclzYNkPBKDYXXrIxmPVYxmDWsbwNpakKyzJEaiKUmMlRbWylO6fg8xh9w+XBvAa665Jtx4443hyCOPDG9605swgD1q0mxcycEmickVeo32kkIkHZeGvTR3X9y0rs+4eUvWTBLTxzT+u5arlRY0ebRjlvCwipGsz2gMBnAlAQ1DbgCtFDxcHtcGcDnmmZkZDCAG0HTnSQ5RaYeaYi3N3Rc3CUW/bw7Sf5esmSRG0p+Wq5UWNHm0Y5bwsIqRrA8GcG3aGoYYQCsFD5cHA/h/rDGA/aIrVTAmoaj009k1QnKISvNq2Etz98VN6/pwA/jEX3so/Zm3Po1Z/rtk30kMoCaP1gRp91ipc2HIuVuufUouHgGHgAFMMIA7d+4M8b/FVxTQ3NxcmJ+fD7OzsynaE8WObm7JphQlVgZpDhvJwSaJUQ65ajPL9dKwt5r8tK4PBhADuNoeGTVyMU76o+CSvSo52zV5Su750dzac0EyL6t5rJUHA4gBXNKH5AZww4YNYePGjbtoCgO4+jaTHBKSmCEOBOs+MIDWRMvnk6yZJEYyUq3urQqoJo92zBIeVjGS9eEGcG3aGoba20+rdU/NgwHEACYZQG4Ay9wYTEJRST1cYrzkEJXm1RRrae6+uGldH24Ay+znPj2V/nfJvsMAYgAxgBjAJAM4umVKC0jymKD0Ybo8v8aESMyDJGbIeVr1JSlE0r407KW5++KmdX2sDOBontKfp7PSgibPkFqQ7J++8cS1kTy6HbeGFv2vpgXJ2S5ZH0lM3/622gfcAGpI123j+jOADz/8cLjrrru6FTjiiCPCJz7xiXDssceGvfbaK+y///69K4MB7EUk+qmLvkNcWlD7RzNsRN+8VitOfQfy0Dz65jH0eEquoqboYwDTbxIla6hZC+1n9zCAKwlYsJeaX4kWSsSUrt8lxmyd07UBvO666zrDN/o67bTTwsUXX9zLurSAJO8SewdpGKB5tykxD5IYw2kMlqpvXhjAwZZC3JGm8GEAMYCLGtDegvWdFUObKc0+0M5dvDmNA0vXb+PhFknn2gDmEi0tIAzg4ytU+oZJcthptNJ3qKfMTWO+NWMe16ZvHqXXRzMP7ZpK2o3GYAAxgBjAx89pyTklidHs+dQ2pet36nhqxGMAM6iXFhAGEAO4KM+ah6ZnAzjO3GEAn6BSsuhrzDiPgDMK2rKmFuyHvrVMnXnp+p06nhrxGMAM6qUFhAHEAE6TAZQUlYztuKKpti+JuZPESOahNdZWbwY0eYYcs2QN+8az2i27ZA0t+peaoL55SPNIdCeJGXLukvGUiCldv0uM2TonBjCDaGkBYQAxgBhA3QaVFLBxmTXGYNwtoWTUtYs+BnD1VdLeJGpvRK20oNX9KAlJHqsxS/ZKiZjS9bvEmK1zYgAziJYWEAYQA4gB1G1QSQEb2gBK9rPElEliJNQ0eYYs+pI17BuPtxtACbNUbUgZas2vZDwlYkrX7xJjts6JAcwgWlpAkoKRMfzkpqUKRt8hXvpLBlaHpuRdtLYvDfvkBV6lgdX6aOeumYe2r5I3gJL9LFlnSYyEmSaPVgtWfVntMc06S35PUGuCrLhqdS/hKomRrLMkRqLf3JjS9Tt3fEO0xwBmUC4tIEnByBh+clPNxpUcbJKY5MEmNLA6NFMPyNXeWffdTJU2xKnzkI6nFOc+XlrO4/KOezQ4Gif9vJZW95p92McoZw0l45HESHQniZHoDAO49gGpYag1vwlHtWlo6fptOthCyTCAGWBLCwgD+PjiSIuTdiklh50md1+BT5mbpoBqxtxnFBbHrBlPKc7SMUt4aIwBBrD/Z2AketGcd9o9pllnbgBXKr2PvfSNkGRflogpXb9LjNk6JwYwg2hpAWkOxIzp9DaVHOKp79BX+7C1VV+9kzL+m73L++s7IDGAx0uWRxWjNZsaY4ABxAAuakB7C9Z3VkjPSa3uJee2JEZybktiVJs+sVHp+p04nCrhGMAM7KUFhAF8fHGk7yS1h5+2XZ90+g51DCAGcLmGtOahT4er/bumEPdpWrpXx41Jc971jWe1PaYx+rVvAPvmKv3mskQvkjNRMh7JOkvGUyKmdP0uMWbrnBjADKKlBaQ5EDOm09tUUjAkY9bESA4SiwMpxZT1AesbT0pfEvZ949H+e988pI/oJUVFO0bJ7YQkt8YYjOaVmiAtVystaPKUHLPkXJCss0RnmnWeFgOo4VNy7pJ9WSKmdP0uMWbrnBjADKKlBaQ5EDOm09tUUjAkY9bEYAD7H7P1LqAyQFv0JcVaOaTeZpIi16epcf/Ol0Ce0OHiGxjNudDHXnqb1afN1d5kYQCfWAGpuZPsZyst9G5wg4DS9dtgiMVTYAAzEJcWkMQoZQw/ualmc/cd0NLPtmgLRuqYVysYybBW+WyhhIdkrprxaNtoxywpGNox9bXDAPYRCqK/2ypZw9Q9lnNDmjoeDOCuH7OQ7A1tjJUW+tWbH1G6fuePsHwGDGAG49ICGtIAWm1cyZg1MRJTpDUqksNOI5O+8ax2g9LXTvrIVTPmPs6rjbmvnfRGRzJmyXpJYvrGPO7fvd0AltqrfeyleunbK5YGcHTMQ7957ZtrDjOJsZbEWOlFcg7kxpSu37njG6I9BjCDcmkBSTZTxvBXNMUAPoHDymD1Hdg1DKBGU33zsLzRkehZYu76xlzaGGgNjtU+TOWYs4ZWYy6hzdLr3KczS66SvrR7Q2LuJDGSNZToRaLf3JjS9Tt3fEO0xwBmUC4tIMlmyhg+BvD/CEgOTQ3nvgMbA6j7FrBkvaTs+4qa5NZnXAwGcCUBSdHXnHe117mvfwyg7jOjmvM2tU3p+p06nhrxGMAM6qUFpDkQtdPRHNBWRW7IQ7RvzKvdGGi49s0LA4gBXK4r7c/AWBhiS6PSt8dy+uoz7NrHoFqj37fHc+YqOf8lMRJmVjGS8Uhqjea8TW1Tun6njqdGPAYwg3ppAUk2U8bwuQHkBnDFFwFWM78li5xGvxaGRzpXrTGQmKAhuZbsS3JOSYq+JI/EqGj0oV3n1rhq5l7yW8BaPppzIbVN6fqdOp4a8RjADOqlBaQ5ECXT0W5Kq0NcMq9SffUVZm4AZd8atFofiV4lRV8SoymOWmPQp7PSXyBI5TH0TZWWT+q8xvUj/SLPqF4kWrC8xdWckxqNYwDnw+zsrOYomvg2GMCMJcQA7grP6tCyMhipeTCAGMC1jgSpedAanFS9ah951jYqWj4YwCcIjFtDDKC8oJeu3/KR1IvEAGawLy0giZnSDJ8bwJXUJIdmCc6LZlOyzhJjIBmjpK/UIptze1RrzKsZfc2tz7ibIa3Bkaxza2soGU+teXEDKHtDJ9nzkhgrLUjOhdyY0vU7d3xDtMcAZlAuLSDJZtIMHwM4nQawVpHFAK7UUw6PSVxDyTlVa14YQAzgajWydP3W1Oah22AAM4iXFpDkYNUMHwOIARx3+yi9FbN8fCjRr2Yf9GlcOtdxt3uSW0JuAJ8gUFIvknXWrqFknfv6t3wzIOlL8jSjZIxkr0reDEjOhdyY0vU7d3xDtMcAZlAuLSDJZtIMX3KQSAqYJEbSlyTGqq++PKsZgxKcx5mw0l8O0Giq5PpIuErGnBqDATxe9KfgUrmW1u+oXvq0yQ0gN4DcAK5+ymIAJRVolZgWDaD23Z3kXZkmpu+ALl0wUseMAZQVjFSu0i8r9Bl0qV76dDcpBrBvHlKukjxDsc+5FcMArn2zqj3/JVwlMZo3DBklOKtp6fqdNbiBGmMAM0CXFpBkM0k2pSRGU9AnoWCkzgsDiAFc60gY+lvAEuNmFTMJ+zn1LOMGULafJVwlMZKaJTmTM8qyuGnp+i0eSMVADGAG/NICkmwmyaaUxEg2pWQ8pWKsilNfHgygrGBY6UWy/UpoatpvADXM+vaG1e0rN4Bhl8fvViaeG0DJifJ4TOn6LR9JvUgMYAb70gKSHOIScyeJsSrokjFrYoYqThhADCA3gE8Q0OxViZnR7ufUs4wbQNl+lnCVxFjpJaMsi5uWrt/igVQMxABmwC8toFKbqeQBbTXmWoYUAygrGFbrI9l+JTTFDWDdL4FIziDtbdZou1GNSR/ja/JYfuPZSvcS42YVU2rMknMiNaZ0/U4dT414DGAG9dICKrWZJIev9h261ZitDEbqeHIMoFXBqL0+kmJgtT6S7Ze6hpJHlRhADKDFTW/tvSrZG5L9bBUjGY8kRnIu5MaUrt+54xuiPQYwg3JpAUk2iiRGsrmtCrpkPJqYoQyp1AD2HfzjxisxJi3EDKkXyfbT6EWyPpK/gaq9PRpKr5Z6mYQxp2pTug+169yns5zPO1rpXsLMKqbUmCXnRGpM6fqdOp4a8RjADOqlBVRqM03CoWU199Q8GEAeAVvcDE2CmUrdG5ZmU3IGjYuRGBXNTfw0G0CrdU5ln6OXjLIsblq6fosHUjEQA5gBv7SAam5crwVsNQNYqqhI/6C7lRYkeTQHvVYvku0nGXNqTOl11vJInUdOkZ3EvlK1yQ3gro/6hzLfOdqUnAu5MaXrd+74hmiPAcygXFpAVgd06qHZwmMLq7mn5iltDLSHb+o8LA9fyZi1hkejTcl4+mJKr7OWx5DrPIl9peoFA2hnAIfUS0ZZFjctXb/FA6kYiAHMgF9aQFYbLvXQxAD2PwYdZar9ZmGLN4Aa3WkNj0abfeZOYn5LG0CLMUrm4S0mVS8YQAzgaiW+dP3OsBaDNcUAZqAuLSBNIbb43AwGEAO4uC2kBgMD+AQBKTNM4urMrP7EHQYQA4gBXN3kYAAxgKHF366yMr+pebgZGrZgpN7oWJkr1nnYdU7dhxjAyVyfUuucUaZXbVr6AqfEmK1zYgAziJYWkNVmIs/jiywxDy3cGHAzlH4zlKpx1tmHwWCdp2OdM8o0BnANeBjADGVhAK9eQW8aPs9GwZiOgtFnolln1nmto9/7Z3r79o/kzbRlTEaZxgBiAEvIp/wfk0691bDccK0dAEONB2OAMcAYjL8xn8Tzhf08Hfu5RAUvfYFTYszWObkBzCBaWkAYQPmjW6viRMGYjoLR94aBdWadMfqTY/QzyjQ3gNwArk7g05/+dNi8eXO4//77w2GHHRYuvPDC8KpXvUqkNwwgj4AXhcIjo8nSAgYQrmKVBQAADSFJREFUA4gBxACuW7cuzM/Ph9nZWVHNn7Yg1zeAf/EXfxF+/dd/PUQT+MpXvjJs3bo1XHTRReH2228P+++/f+9aYwAnq+hLbgkxBhgDjMHkGANuep84gyXn26TG9BZjRUDp+q0Y0uBNXBvAl73sZeHII48MW7ZsWQL/whe+MJx44onhggsu6F2M0gLiETCPgCOBafhyTWuFB6OP0cfoT47R7y3GioDS9VsxpMGbuDWAjzzySHjKU54SvvjFL4aTTjppCfwHPvCBcMstt4Rt27btshg7d+4M8b/FV7w6jjeF27dvL3KFfPj/+8pSX9/deFxY/v/Hfxj934h5gpeWz7gdOI7raJw2hjWso3HWeeV5Mq1nB+s8HetcwhlFAzg3NxceeuihEB8Fe3y5NYD33Xdf2HfffcONN94Yjj766KW1P//888MXvvCF8L3vfW8XPWzYsCFs3LjRo06YMwQgAAEIQGDqCMQLnP3222/q5iWZkHsDeNNNN4VXvOIVS6x+//d/P/zJn/xJuOOOO3pvAB977LHwX//1X+GZz3xmmJmZkfBOill8h1LqhjFpMFMeDOthFhjOcB6GwHC9oOlhWFtzXlhYCDt27Aj77LNP2G233YaZRGO9uDWAmkfAQ68dn1EYjjish2ENZzgPQ2C4XtD0MKzhbM/ZrQGMKOOXQF7ykpd03wJefB166KHhhBNOEH0JxH45VmZE8KUJP5Ef1sOwhjOchyEwXC9oehjWcLbn7NoALv4MzB/90R91j4E/85nPhM9+9rPhtttuCwcccIA97cSMCD4RWEY4rDPgJTSFcwKsjFA4Z8BLbArrRGDKcDgrwa3RzLUBjFzi7d/HP/7x7oegDz/88PDJT34yvPrVr7YnrcgYv3Ecf47mIx/5SHjSk56kyEATKQFYS0nlxcE5j5+0NZylpPLjYJ3PUJIBzhJKaTHuDWAaLqIhAAEIQAACEIDA5BPAAE7+GjIDCEAAAhCAAAQgkEQAA5iEi2AIQAACEIAABCAw+QQwgJO/hswAAhCAAAQgAAEIJBHAACbhIhgCEIAABCAAAQhMPgEMYMNrGL+hvHnz5u4byocddli48MILw6te9aqGR9z20OI3qi+77LLur7w8+clP7v4E4Mc+9rHw/Oc/f2ng8ZtmH/rQh8Kf//mfh5/85CfhNa95TfdNca9/KshiRSP3c845J8S/sx01HF9wtiD7eI4f/vCH4ayzzgrXXHNNp9lDDjkkfO5zn+t+4zS+4l88iH/CMv7M1YMPPtj9/umnPvWp7kzhJSPws5/9LMQ/Bfpnf/Zn4Uc/+lH4hV/4hfD2t789/M7v/M7SX5GAs4zlaNT111/f1blvf/vbXa27/PLLw4knnrgUJuEadf3+978//M3f/E3X7ld/9VfDH/zBH4SnP/3pukE5aYUBbHShF3+jMJqPV77ylWHr1q3hoosuCrfffnvYf//9Gx1128N6/etfH0455ZTw0pe+NMQD/aMf/Wi49dZbO6ZPfepTu8G/5z3vCVdeeWW4+OKLuz/x91u/9Vvdn/uLh9Puu+/e9gQbHN3NN98cTj755DA7OxuOPfbYJQMIZ5vFioXviCOO6NhGpnvvvXf4/ve/Hw488MBw8MEHd53ENznxT1xGTUdzeN5554VYdOPfO3/a055mM5ApzxL5xZ8Ii38nPhrnb33rW+Ed73hHxzK+sYGzXgDxjcuNN94YjjzyyPCmN71pFwMo0e8b3vCGcO+993ZvcuLr3e9+d7cH4lnOa3UCGMBG1RHfpccNsWXLlqURvvCFL+zeGcUbFV75BB544IGuYG7btq377cf5+fnwrGc9q/tb0L/2a7/WdXDfffeFubm58Ld/+7fhuOOOy+/UUYaHH36403B8ExML5Ytf/OLOAMLZTgRnn312VzxvuOGGsUnj7Un8W6dnnHFGd0sYX/H29dnPfnZnDH/zN3/TbjBTnOlXfuVXOmbxZnXxFc3KU57ylO68gLPN4s/MzKwwgBKu//Iv/xLiX/D6x3/8x+52O77i/x3/uEN82rP8CY/NKKcnCwawwbWchL9T3CC25CHddddd4XnPe153Cxh/BPzv//7vu0e+8cbvGc94xlK+F73oRZ3xjo/ReMkJnHbaaWGvvfbqbk6OOeaYJQMIZznDvshY+OIbk3j7Ed/I7LvvvuH0008P73rXu7qm//Zv/9bdBP7TP/1Td1O4+Ip/7jI+Hos3Wrz6CWzatCnEvxh17bXXdreo//zP/xxe97rXdW9o3vKWt8C5H6EoYtQASvT7x3/8x+HMM88MDz300Io+or7j2RNvanmNJ4ABbFAZ8dYpHuTxnX38nNri6/zzz+8O7PjohlcegfjOMhbB+Aht8fbkkksu6Q6LeEOy/BUP+oMOOqh7DM9LRuDSSy/tbv3io7I999xzhQGEs4yhJCqyja9YAN/85jeHb37zm91tX9Tq2972tnDTTTd1HyGJnxOMN4GLr/iI7J577glf+cpXJN24j4nnRfwca7w1jR8FefTRR7vH6vGvNMUXnG0kMmoAJVxjXYwfb/jXf/3XFYOIRj2e54trZDPC6cqCAWxwPRcNYBR/vMZefMUDJz5uiNfavPIIvPe97w1XX311+NrXvrb0BY/VjMn69eu7W5R4A8Crn8D27dvDUUcd1d2WxNvT+Fp+AwjnfobSiD322KNjHc+KxVf8MHz87OXXv/71JWMSz5T4xYXFV7whjOv05S9/WdqV67j4hua3f/u3uy8rxM8A3nLLLZ3R/sQnPhHiTfeiUYFznkxWM4BrcV3tYiQ+3XnnO98Z4sckeI0ngAFsUBk8Ai67KO973/vCFVdc0X0QPt7sLb54NGnDPbI96aSTVnxpJt6YxMN9t912626dXvva1/Ko3QD3AQccEOIblPgFscVX/NxwvH2Nt36SR2gGw5j6FPFzwNFIxDeOi6/I+E//9E+7N+RwtpEAj4BtOEqzYAClpAaOix9mjT/jED9Av/iKn/eJjy35EohuMeJjnGj+4s8MXHfddd3n/5a/Fr+cEA/1+M3V+Io/SxB/AoYvgciZ79ixo3u8uPwVH8W84AUv6L6IEItp/LINnOVMV4s89dRTu5u85V8C+eAHPxi+8Y1vdLdSix+ij//bhz/84S5NfIMZv/zEl0Dk/OMvAkTDF79pvfiK5/DnP//57tEjnOUs14pc7Usga+l38UsgUfO/9Eu/1KWP//fLX/5yvgTSsywYQBvdmmdZ/BmY+NgxPgaOX2//7Gc/G2677bYQ3/XzSicQPxwfHz9+6UtfWvHNsHXr1nW/Cxhf8YC/6qqrus+UxC8wxN8E/M///E9+BiYd94oWyx8BwzkT5rLm8VFv/Jxw/IJSfNMSPwMYH+/G8+Ktb31rFxmN3qJZiW964iOz+AaIn4GRr0P8zb+vfvWr3Wcr4yPg73znO91PjfzGb/xGxxfOcpajkfHXAuIX8uIrflEpPlaPP2sUz9/4k2cS/cafgYmPiRc/px3XJtZJfgZm7XXBAOp1W7xlvP37+Mc/3t1CxW+pxm80xZ8r4aUjEN9djnvFd/HxgI+vn/70p91nfaJRXP5D0PHWipeewKgBhLOe5WjL+IYlftD9zjvv7D7SEL8Qsvgt4Bi7+EO6sTgu/yHoeKbwkhGIt9q/+7u/2z09+PGPf9x9oSZ++/fcc88N8XOYcJZxHBcV34xEwzf6ip+tjG/EJfqNv9ww+kPQf/iHf8gPQfcsCwZQr1taQgACEIAABCAAgYkkgAGcyGVj0BCAAAQgAAEIQEBPAAOoZ0dLCEAAAhCAAAQgMJEEMIATuWwMGgIQgAAEIAABCOgJYAD17GgJAQhAAAIQgAAEJpIABnAil41BQwACEIAABCAAAT0BDKCeHS0hAAEIQAACEIDARBLAAE7ksjFoCEAAAhCAAAQgoCeAAdSzoyUEIAABCEAAAhCYSAIYwIlcNgYNAQhAAAIQgAAE9AQwgHp2tIQABCAAAQhAAAITSQADOJHLxqAhAAEIQAACEICAngAGUM+OlhCAAAQgAAEIQGAiCWAAJ3LZGDQEIAABCEAAAhDQE8AA6tnREgIQgAAEIAABCEwkAQzgRC4bg4YABCAAAQhAAAJ6AhhAPTtaQgACEIAABCAAgYkkgAGcyGVj0BCAAAQgAAEIQEBPAAOoZ0dLCEAAAhCAAAQgMJEEMIATuWwMGgIQgAAEIAABCOgJYAD17GgJAQhAAAIQgAAEJpIABnAil41BQwACEIAABCAAAT0BDKCeHS0hAAEIQAACEIDARBLAAE7ksjFoCEAAAhCAAAQgoCeAAdSzoyUEIAABCEAAAhCYSAIYwIlcNgYNAQhAAAIQgAAE9AQwgHp2tIQABCAAAQhAAAITSQADOJHLxqAhAAEIQAACEICAngAGUM+OlhCAAAQgAAEIQGAiCWAAJ3LZGDQEIAABCEAAAhDQE8AA6tnREgIQgAAEIAABCEwkAQzgRC4bg4YABCAAAQhAAAJ6Av8fs0hViQo7h3IAAAAASUVORK5CYII=\" width=\"640\">"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# You may wish to view the PageRank graphically.\n",
    "# This code will draw a bar chart, for each (numbered) website on the generated internet,\n",
    "# The height of each bar will be the score in the PageRank.\n",
    "# Run this code to see the PageRank for each internet you generate.\n",
    "# Hopefully you should see what you might expect\n",
    "# - there are a few clusters of important websites, but most on the internet are rubbish!\n",
    "%pylab notebook\n",
    "r = pageRank(generate_internet(100), 0.9)\n",
    "plt.bar(arange(r.shape[0]), r);"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.6.5"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}