Example-2.2.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "b16ccb14",
   "metadata": {},
   "outputs": [],
   "source": [
    "import numpy as np\n",
    "from utils import *\n",
    "%matplotlib notebook"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "6a06026c",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "r = k/h**2 = 1\n",
      "Number of discretization on x: 11\n",
      "Number of discretization on t: 21\n"
     ]
    }
   ],
   "source": [
    "L = 1\n",
    "T = .2\n",
    "\n",
    "h = .1\n",
    "k = 0.01\n",
    "\n",
    "r = k / h**2\n",
    "print('r = k/h**2 = %g' % r)\n",
    "\n",
    "Nx = int(np.ceil(L / h)) + 1\n",
    "Nt = int(np.ceil(T / k)) + 1\n",
    "\n",
    "print('Number of discretization on x:', Nx)\n",
    "print('Number of discretization on t:', Nt)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "3e255551",
   "metadata": {},
   "outputs": [],
   "source": [
    "@np.vectorize\n",
    "def initial(x):\n",
    "    if 0 <= x <= 1 / 2:\n",
    "        return 2 * x\n",
    "    elif 1 / 2 <= x <= 1:\n",
    "        return 2 * (1 - x)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "1093e809",
   "metadata": {},
   "outputs": [],
   "source": [
    "x = np.array([i * h for i in range(Nx)])\n",
    "t = np.array([i * k for i in range(Nt)])\n",
    "\n",
    "u0 = initial(x)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "2e15648b",
   "metadata": {},
   "outputs": [],
   "source": [
    "u = np.zeros((Nt, Nx))\n",
    "u[0, :] = u0\n",
    "N = Nx - 2\n",
    "\n",
    "A = (2 + 2 * r) * np.eye(N) - r * np.eye(N, k=1) - r * np.eye(N, k=-1)\n",
    "B = (2 - 2 * r) * np.eye(N) + r * np.eye(N, k=1) + r * np.eye(N, k=-1)\n",
    "\n",
    "for j in range(Nt - 1):\n",
    "    u[j + 1, 1:-1] = np.linalg.solve(A, B @ u[j, 1:-1])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "be837689",
   "metadata": {
    "scrolled": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "╒══════════════╤═══════╤════════╤════════╤════════╤════════╤════════╤════════╤════════╤════════╤════════╤═══════╕\n",
      "│      x       │  0.0  │  0.1   │  0.2   │  0.3   │  0.4   │  0.5   │  0.6   │  0.7   │  0.8   │  0.9   │  1.0  │\n",
      "╞══════════════╪═══════╪════════╪════════╪════════╪════════╪════════╪════════╪════════╪════════╪════════╪═══════╡\n",
      "│   j=0(t=0)   │   0   │  0.2   │  0.4   │  0.6   │  0.8   │   1    │  0.8   │  0.6   │  0.4   │  0.2   │   0   │\n",
      "├──────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=1(t=0.01)  │   0   │ 0.1989 │ 0.3956 │ 0.5834 │ 0.7381 │ 0.7691 │ 0.7381 │ 0.5834 │ 0.3956 │ 0.1989 │   0   │\n",
      "├──────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=2(t=0.02)  │   0   │ 0.1936 │ 0.3789 │ 0.5397 │ 0.6461 │ 0.6921 │ 0.6461 │ 0.5397 │ 0.3789 │ 0.1936 │   0   │\n",
      "├──────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=3(t=0.03)  │   0   │ 0.1826 │ 0.3515 │ 0.4902 │ 0.5843 │ 0.6152 │ 0.5843 │ 0.4902 │ 0.3515 │ 0.1826 │   0   │\n",
      "├──────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=4(t=0.04)  │   0   │ 0.1683 │ 0.3218 │ 0.4461 │ 0.5267 │ 0.5555 │ 0.5267 │ 0.4461 │ 0.3218 │ 0.1683 │   0   │\n",
      "├──────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=5(t=0.05)  │   0   │ 0.1538 │ 0.2932 │ 0.4047 │ 0.477  │ 0.5019 │ 0.477  │ 0.4047 │ 0.2932 │ 0.1538 │   0   │\n",
      "├──────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=6(t=0.06)  │   0   │ 0.1399 │ 0.2664 │ 0.3672 │ 0.4321 │ 0.4546 │ 0.4321 │ 0.3672 │ 0.2664 │ 0.1399 │   0   │\n",
      "├──────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=7(t=0.07)  │   0   │ 0.127  │ 0.2418 │ 0.333  │ 0.3916 │ 0.4119 │ 0.3916 │ 0.333  │ 0.2418 │ 0.127  │   0   │\n",
      "├──────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=8(t=0.08)  │   0   │ 0.1153 │ 0.2193 │ 0.3019 │ 0.355  │ 0.3733 │ 0.355  │ 0.3019 │ 0.2193 │ 0.1153 │   0   │\n",
      "├──────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=9(t=0.09)  │   0   │ 0.1045 │ 0.1989 │ 0.2738 │ 0.3219 │ 0.3385 │ 0.3219 │ 0.2738 │ 0.1989 │ 0.1045 │   0   │\n",
      "├──────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=10(t=0.1)  │   0   │ 0.0948 │ 0.1803 │ 0.2482 │ 0.2918 │ 0.3069 │ 0.2918 │ 0.2482 │ 0.1803 │ 0.0948 │   0   │\n",
      "├──────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=11(t=0.11) │   0   │ 0.086  │ 0.1635 │ 0.2251 │ 0.2646 │ 0.2782 │ 0.2646 │ 0.2251 │ 0.1635 │ 0.086  │   0   │\n",
      "├──────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=12(t=0.12) │   0   │ 0.0779 │ 0.1483 │ 0.2041 │ 0.2399 │ 0.2522 │ 0.2399 │ 0.2041 │ 0.1483 │ 0.0779 │   0   │\n",
      "├──────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=13(t=0.13) │   0   │ 0.0707 │ 0.1344 │ 0.185  │ 0.2175 │ 0.2287 │ 0.2175 │ 0.185  │ 0.1344 │ 0.0707 │   0   │\n",
      "├──────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=14(t=0.14) │   0   │ 0.0641 │ 0.1219 │ 0.1678 │ 0.1972 │ 0.2074 │ 0.1972 │ 0.1678 │ 0.1219 │ 0.0641 │   0   │\n",
      "├──────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=15(t=0.15) │   0   │ 0.0581 │ 0.1105 │ 0.1521 │ 0.1788 │ 0.188  │ 0.1788 │ 0.1521 │ 0.1105 │ 0.0581 │   0   │\n",
      "├──────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=16(t=0.16) │   0   │ 0.0527 │ 0.1002 │ 0.1379 │ 0.1621 │ 0.1705 │ 0.1621 │ 0.1379 │ 0.1002 │ 0.0527 │   0   │\n",
      "├──────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=17(t=0.17) │   0   │ 0.0478 │ 0.0908 │ 0.125  │ 0.147  │ 0.1546 │ 0.147  │ 0.125  │ 0.0908 │ 0.0478 │   0   │\n",
      "├──────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=18(t=0.18) │   0   │ 0.0433 │ 0.0824 │ 0.1134 │ 0.1333 │ 0.1401 │ 0.1333 │ 0.1134 │ 0.0824 │ 0.0433 │   0   │\n",
      "├──────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=19(t=0.19) │   0   │ 0.0393 │ 0.0747 │ 0.1028 │ 0.1208 │ 0.1271 │ 0.1208 │ 0.1028 │ 0.0747 │ 0.0393 │   0   │\n",
      "├──────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│ j=20(t=0.2)  │   0   │ 0.0356 │ 0.0677 │ 0.0932 │ 0.1096 │ 0.1152 │ 0.1096 │ 0.0932 │ 0.0677 │ 0.0356 │   0   │\n",
      "├──────────────┼───────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼────────┼───────┤\n",
      "│      x       │   0   │  0.1   │  0.2   │  0.3   │  0.4   │  0.5   │  0.6   │  0.7   │  0.8   │  0.9   │   1   │\n",
      "╘══════════════╧═══════╧════════╧════════╧════════╧════════╧════════╧════════╧════════╧════════╧════════╧═══════╛\n"
     ]
    }
   ],
   "source": [
    "print_table(u, x, t)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "id": "66ec12c2",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "application/javascript": [
       "/* Put everything inside the global mpl namespace */\n",
       "/* global mpl */\n",
       "window.mpl = {};\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(\n",
       "            '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",
       "\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 = document.createElement('div');\n",
       "    this.root.setAttribute('style', 'display: inline-block');\n",
       "    this._root_extra_style(this.root);\n",
       "\n",
       "    parent_element.appendChild(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 (fig.ratio !== 1) {\n",
       "            fig.send_message('set_dpi_ratio', { dpi_ratio: fig.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 = document.createElement('div');\n",
       "    titlebar.classList =\n",
       "        'ui-dialog-titlebar ui-widget-header ui-corner-all ui-helper-clearfix';\n",
       "    var titletext = document.createElement('div');\n",
       "    titletext.classList = 'ui-dialog-title';\n",
       "    titletext.setAttribute(\n",
       "        'style',\n",
       "        'width: 100%; text-align: center; padding: 3px;'\n",
       "    );\n",
       "    titlebar.appendChild(titletext);\n",
       "    this.root.appendChild(titlebar);\n",
       "    this.header = titletext;\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._canvas_extra_style = function (_canvas_div) {};\n",
       "\n",
       "mpl.figure.prototype._root_extra_style = function (_canvas_div) {};\n",
       "\n",
       "mpl.figure.prototype._init_canvas = function () {\n",
       "    var fig = this;\n",
       "\n",
       "    var canvas_div = (this.canvas_div = document.createElement('div'));\n",
       "    canvas_div.setAttribute(\n",
       "        'style',\n",
       "        'border: 1px solid #ddd;' +\n",
       "            'box-sizing: content-box;' +\n",
       "            'clear: both;' +\n",
       "            'min-height: 1px;' +\n",
       "            'min-width: 1px;' +\n",
       "            'outline: 0;' +\n",
       "            'overflow: hidden;' +\n",
       "            'position: relative;' +\n",
       "            'resize: both;'\n",
       "    );\n",
       "\n",
       "    function on_keyboard_event_closure(name) {\n",
       "        return function (event) {\n",
       "            return fig.key_event(event, name);\n",
       "        };\n",
       "    }\n",
       "\n",
       "    canvas_div.addEventListener(\n",
       "        'keydown',\n",
       "        on_keyboard_event_closure('key_press')\n",
       "    );\n",
       "    canvas_div.addEventListener(\n",
       "        'keyup',\n",
       "        on_keyboard_event_closure('key_release')\n",
       "    );\n",
       "\n",
       "    this._canvas_extra_style(canvas_div);\n",
       "    this.root.appendChild(canvas_div);\n",
       "\n",
       "    var canvas = (this.canvas = document.createElement('canvas'));\n",
       "    canvas.classList.add('mpl-canvas');\n",
       "    canvas.setAttribute('style', 'box-sizing: content-box;');\n",
       "\n",
       "    this.context = canvas.getContext('2d');\n",
       "\n",
       "    var backingStore =\n",
       "        this.context.backingStorePixelRatio ||\n",
       "        this.context.webkitBackingStorePixelRatio ||\n",
       "        this.context.mozBackingStorePixelRatio ||\n",
       "        this.context.msBackingStorePixelRatio ||\n",
       "        this.context.oBackingStorePixelRatio ||\n",
       "        this.context.backingStorePixelRatio ||\n",
       "        1;\n",
       "\n",
       "    this.ratio = (window.devicePixelRatio || 1) / backingStore;\n",
       "\n",
       "    var rubberband_canvas = (this.rubberband_canvas = document.createElement(\n",
       "        'canvas'\n",
       "    ));\n",
       "    rubberband_canvas.setAttribute(\n",
       "        'style',\n",
       "        'box-sizing: content-box; position: absolute; left: 0; top: 0; z-index: 1;'\n",
       "    );\n",
       "\n",
       "    // Apply a ponyfill if ResizeObserver is not implemented by browser.\n",
       "    if (this.ResizeObserver === undefined) {\n",
       "        if (window.ResizeObserver !== undefined) {\n",
       "            this.ResizeObserver = window.ResizeObserver;\n",
       "        } else {\n",
       "            var obs = _JSXTOOLS_RESIZE_OBSERVER({});\n",
       "            this.ResizeObserver = obs.ResizeObserver;\n",
       "        }\n",
       "    }\n",
       "\n",
       "    this.resizeObserverInstance = new this.ResizeObserver(function (entries) {\n",
       "        var nentries = entries.length;\n",
       "        for (var i = 0; i < nentries; i++) {\n",
       "            var entry = entries[i];\n",
       "            var width, height;\n",
       "            if (entry.contentBoxSize) {\n",
       "                if (entry.contentBoxSize instanceof Array) {\n",
       "                    // Chrome 84 implements new version of spec.\n",
       "                    width = entry.contentBoxSize[0].inlineSize;\n",
       "                    height = entry.contentBoxSize[0].blockSize;\n",
       "                } else {\n",
       "                    // Firefox implements old version of spec.\n",
       "                    width = entry.contentBoxSize.inlineSize;\n",
       "                    height = entry.contentBoxSize.blockSize;\n",
       "                }\n",
       "            } else {\n",
       "                // Chrome <84 implements even older version of spec.\n",
       "                width = entry.contentRect.width;\n",
       "                height = entry.contentRect.height;\n",
       "            }\n",
       "\n",
       "            // Keep the size of the canvas and rubber band canvas in sync with\n",
       "            // the canvas container.\n",
       "            if (entry.devicePixelContentBoxSize) {\n",
       "                // Chrome 84 implements new version of spec.\n",
       "                canvas.setAttribute(\n",
       "                    'width',\n",
       "                    entry.devicePixelContentBoxSize[0].inlineSize\n",
       "                );\n",
       "                canvas.setAttribute(\n",
       "                    'height',\n",
       "                    entry.devicePixelContentBoxSize[0].blockSize\n",
       "                );\n",
       "            } else {\n",
       "                canvas.setAttribute('width', width * fig.ratio);\n",
       "                canvas.setAttribute('height', height * fig.ratio);\n",
       "            }\n",
       "            canvas.setAttribute(\n",
       "                'style',\n",
       "                'width: ' + width + 'px; height: ' + height + 'px;'\n",
       "            );\n",
       "\n",
       "            rubberband_canvas.setAttribute('width', width);\n",
       "            rubberband_canvas.setAttribute('height', height);\n",
       "\n",
       "            // And update the size in Python. We ignore the initial 0/0 size\n",
       "            // that occurs as the element is placed into the DOM, which should\n",
       "            // otherwise not happen due to the minimum size styling.\n",
       "            if (fig.ws.readyState == 1 && width != 0 && height != 0) {\n",
       "                fig.request_resize(width, height);\n",
       "            }\n",
       "        }\n",
       "    });\n",
       "    this.resizeObserverInstance.observe(canvas_div);\n",
       "\n",
       "    function on_mouse_event_closure(name) {\n",
       "        return function (event) {\n",
       "            return fig.mouse_event(event, name);\n",
       "        };\n",
       "    }\n",
       "\n",
       "    rubberband_canvas.addEventListener(\n",
       "        'mousedown',\n",
       "        on_mouse_event_closure('button_press')\n",
       "    );\n",
       "    rubberband_canvas.addEventListener(\n",
       "        'mouseup',\n",
       "        on_mouse_event_closure('button_release')\n",
       "    );\n",
       "    rubberband_canvas.addEventListener(\n",
       "        'dblclick',\n",
       "        on_mouse_event_closure('dblclick')\n",
       "    );\n",
       "    // Throttle sequential mouse events to 1 every 20ms.\n",
       "    rubberband_canvas.addEventListener(\n",
       "        'mousemove',\n",
       "        on_mouse_event_closure('motion_notify')\n",
       "    );\n",
       "\n",
       "    rubberband_canvas.addEventListener(\n",
       "        'mouseenter',\n",
       "        on_mouse_event_closure('figure_enter')\n",
       "    );\n",
       "    rubberband_canvas.addEventListener(\n",
       "        'mouseleave',\n",
       "        on_mouse_event_closure('figure_leave')\n",
       "    );\n",
       "\n",
       "    canvas_div.addEventListener('wheel', function (event) {\n",
       "        if (event.deltaY < 0) {\n",
       "            event.step = 1;\n",
       "        } else {\n",
       "            event.step = -1;\n",
       "        }\n",
       "        on_mouse_event_closure('scroll')(event);\n",
       "    });\n",
       "\n",
       "    canvas_div.appendChild(canvas);\n",
       "    canvas_div.appendChild(rubberband_canvas);\n",
       "\n",
       "    this.rubberband_context = rubberband_canvas.getContext('2d');\n",
       "    this.rubberband_context.strokeStyle = '#000000';\n",
       "\n",
       "    this._resize_canvas = function (width, height, forward) {\n",
       "        if (forward) {\n",
       "            canvas_div.style.width = width + 'px';\n",
       "            canvas_div.style.height = height + 'px';\n",
       "        }\n",
       "    };\n",
       "\n",
       "    // Disable right mouse context menu.\n",
       "    this.rubberband_canvas.addEventListener('contextmenu', function (_e) {\n",
       "        event.preventDefault();\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 toolbar = document.createElement('div');\n",
       "    toolbar.classList = 'mpl-toolbar';\n",
       "    this.root.appendChild(toolbar);\n",
       "\n",
       "    function on_click_closure(name) {\n",
       "        return function (_event) {\n",
       "            return fig.toolbar_button_onclick(name);\n",
       "        };\n",
       "    }\n",
       "\n",
       "    function on_mouseover_closure(tooltip) {\n",
       "        return function (event) {\n",
       "            if (!event.currentTarget.disabled) {\n",
       "                return fig.toolbar_button_onmouseover(tooltip);\n",
       "            }\n",
       "        };\n",
       "    }\n",
       "\n",
       "    fig.buttons = {};\n",
       "    var buttonGroup = document.createElement('div');\n",
       "    buttonGroup.classList = 'mpl-button-group';\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",
       "            /* Instead of a spacer, we start a new button group. */\n",
       "            if (buttonGroup.hasChildNodes()) {\n",
       "                toolbar.appendChild(buttonGroup);\n",
       "            }\n",
       "            buttonGroup = document.createElement('div');\n",
       "            buttonGroup.classList = 'mpl-button-group';\n",
       "            continue;\n",
       "        }\n",
       "\n",
       "        var button = (fig.buttons[name] = document.createElement('button'));\n",
       "        button.classList = 'mpl-widget';\n",
       "        button.setAttribute('role', 'button');\n",
       "        button.setAttribute('aria-disabled', 'false');\n",
       "        button.addEventListener('click', on_click_closure(method_name));\n",
       "        button.addEventListener('mouseover', on_mouseover_closure(tooltip));\n",
       "\n",
       "        var icon_img = document.createElement('img');\n",
       "        icon_img.src = '_images/' + image + '.png';\n",
       "        icon_img.srcset = '_images/' + image + '_large.png 2x';\n",
       "        icon_img.alt = tooltip;\n",
       "        button.appendChild(icon_img);\n",
       "\n",
       "        buttonGroup.appendChild(button);\n",
       "    }\n",
       "\n",
       "    if (buttonGroup.hasChildNodes()) {\n",
       "        toolbar.appendChild(buttonGroup);\n",
       "    }\n",
       "\n",
       "    var fmt_picker = document.createElement('select');\n",
       "    fmt_picker.classList = 'mpl-widget';\n",
       "    toolbar.appendChild(fmt_picker);\n",
       "    this.format_dropdown = fmt_picker;\n",
       "\n",
       "    for (var ind in mpl.extensions) {\n",
       "        var fmt = mpl.extensions[ind];\n",
       "        var option = document.createElement('option');\n",
       "        option.selected = fmt === mpl.default_extension;\n",
       "        option.innerHTML = fmt;\n",
       "        fmt_picker.appendChild(option);\n",
       "    }\n",
       "\n",
       "    var status_bar = document.createElement('span');\n",
       "    status_bar.classList = 'mpl-message';\n",
       "    toolbar.appendChild(status_bar);\n",
       "    this.message = status_bar;\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",
       "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",
       "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], msg['forward']);\n",
       "        fig.send_message('refresh', {});\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_rubberband = function (fig, msg) {\n",
       "    var x0 = msg['x0'] / fig.ratio;\n",
       "    var y0 = (fig.canvas.height - msg['y0']) / fig.ratio;\n",
       "    var x1 = msg['x1'] / fig.ratio;\n",
       "    var y1 = (fig.canvas.height - msg['y1']) / fig.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,\n",
       "        0,\n",
       "        fig.canvas.width / fig.ratio,\n",
       "        fig.canvas.height / fig.ratio\n",
       "    );\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",
       "        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.handle_history_buttons = function (fig, msg) {\n",
       "    for (var key in msg) {\n",
       "        if (!(key in fig.buttons)) {\n",
       "            continue;\n",
       "        }\n",
       "        fig.buttons[key].disabled = !msg[key];\n",
       "        fig.buttons[key].setAttribute('aria-disabled', !msg[key]);\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_navigate_mode = function (fig, msg) {\n",
       "    if (msg['mode'] === 'PAN') {\n",
       "        fig.buttons['Pan'].classList.add('active');\n",
       "        fig.buttons['Zoom'].classList.remove('active');\n",
       "    } else if (msg['mode'] === 'ZOOM') {\n",
       "        fig.buttons['Pan'].classList.remove('active');\n",
       "        fig.buttons['Zoom'].classList.add('active');\n",
       "    } else {\n",
       "        fig.buttons['Pan'].classList.remove('active');\n",
       "        fig.buttons['Zoom'].classList.remove('active');\n",
       "    }\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",
       "            var img = evt.data;\n",
       "            if (img.type !== 'image/png') {\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",
       "                img.type = 'image/png';\n",
       "            }\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",
       "\n",
       "            fig.imageObj.src = (window.URL || window.webkitURL).createObjectURL(\n",
       "                img\n",
       "            );\n",
       "            fig.updated_canvas_event();\n",
       "            fig.waiting = false;\n",
       "            return;\n",
       "        } else if (\n",
       "            typeof evt.data === 'string' &&\n",
       "            evt.data.slice(0, 21) === 'data:image/png;base64'\n",
       "        ) {\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(\n",
       "                \"No handler for the '\" + msg_type + \"' message type: \",\n",
       "                msg\n",
       "            );\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(\n",
       "                    \"Exception inside the 'handler_\" + msg_type + \"' callback:\",\n",
       "                    e,\n",
       "                    e.stack,\n",
       "                    msg\n",
       "                );\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",
       "    }\n",
       "    if (e.target) {\n",
       "        targ = e.target;\n",
       "    } else if (e.srcElement) {\n",
       "        targ = e.srcElement;\n",
       "    }\n",
       "    if (targ.nodeType === 3) {\n",
       "        // defeat Safari bug\n",
       "        targ = targ.parentNode;\n",
       "    }\n",
       "\n",
       "    // pageX,Y are the mouse positions relative to the document\n",
       "    var boundingRect = targ.getBoundingClientRect();\n",
       "    var x = e.pageX - (boundingRect.left + document.body.scrollLeft);\n",
       "    var y = e.pageY - (boundingRect.top + document.body.scrollTop);\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",
       "        }\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",
       "        this.canvas.focus();\n",
       "        this.canvas_div.focus();\n",
       "    }\n",
       "\n",
       "    var x = canvas_pos.x * this.ratio;\n",
       "    var y = canvas_pos.y * this.ratio;\n",
       "\n",
       "    this.send_message(name, {\n",
       "        x: x,\n",
       "        y: y,\n",
       "        button: event.button,\n",
       "        step: event.step,\n",
       "        guiEvent: simpleKeys(event),\n",
       "    });\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",
       "    // Prevent repeat events\n",
       "    if (name === 'key_press') {\n",
       "        if (event.key === this._key) {\n",
       "            return;\n",
       "        } else {\n",
       "            this._key = event.key;\n",
       "        }\n",
       "    }\n",
       "    if (name === 'key_release') {\n",
       "        this._key = null;\n",
       "    }\n",
       "\n",
       "    var value = '';\n",
       "    if (event.ctrlKey && event.key !== 'Control') {\n",
       "        value += 'ctrl+';\n",
       "    }\n",
       "    else if (event.altKey && event.key !== 'Alt') {\n",
       "        value += 'alt+';\n",
       "    }\n",
       "    else if (event.shiftKey && event.key !== 'Shift') {\n",
       "        value += 'shift+';\n",
       "    }\n",
       "\n",
       "    value += 'k' + event.key;\n",
       "\n",
       "    this._key_event_extra(event, name);\n",
       "\n",
       "    this.send_message(name, { key: value, 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",
       "\n",
       "///////////////// REMAINING CONTENT GENERATED BY embed_js.py /////////////////\n",
       "// prettier-ignore\n",
       "var _JSXTOOLS_RESIZE_OBSERVER=function(A){var t,i=new WeakMap,n=new WeakMap,a=new WeakMap,r=new WeakMap,o=new Set;function s(e){if(!(this instanceof s))throw new TypeError(\"Constructor requires 'new' operator\");i.set(this,e)}function h(){throw new TypeError(\"Function is not a constructor\")}function c(e,t,i,n){e=0 in arguments?Number(arguments[0]):0,t=1 in arguments?Number(arguments[1]):0,i=2 in arguments?Number(arguments[2]):0,n=3 in arguments?Number(arguments[3]):0,this.right=(this.x=this.left=e)+(this.width=i),this.bottom=(this.y=this.top=t)+(this.height=n),Object.freeze(this)}function d(){t=requestAnimationFrame(d);var s=new WeakMap,p=new Set;o.forEach((function(t){r.get(t).forEach((function(i){var r=t instanceof window.SVGElement,o=a.get(t),d=r?0:parseFloat(o.paddingTop),f=r?0:parseFloat(o.paddingRight),l=r?0:parseFloat(o.paddingBottom),u=r?0:parseFloat(o.paddingLeft),g=r?0:parseFloat(o.borderTopWidth),m=r?0:parseFloat(o.borderRightWidth),w=r?0:parseFloat(o.borderBottomWidth),b=u+f,F=d+l,v=(r?0:parseFloat(o.borderLeftWidth))+m,W=g+w,y=r?0:t.offsetHeight-W-t.clientHeight,E=r?0:t.offsetWidth-v-t.clientWidth,R=b+v,z=F+W,M=r?t.width:parseFloat(o.width)-R-E,O=r?t.height:parseFloat(o.height)-z-y;if(n.has(t)){var k=n.get(t);if(k[0]===M&&k[1]===O)return}n.set(t,[M,O]);var S=Object.create(h.prototype);S.target=t,S.contentRect=new c(u,d,M,O),s.has(i)||(s.set(i,[]),p.add(i)),s.get(i).push(S)}))})),p.forEach((function(e){i.get(e).call(e,s.get(e),e)}))}return s.prototype.observe=function(i){if(i instanceof window.Element){r.has(i)||(r.set(i,new Set),o.add(i),a.set(i,window.getComputedStyle(i)));var n=r.get(i);n.has(this)||n.add(this),cancelAnimationFrame(t),t=requestAnimationFrame(d)}},s.prototype.unobserve=function(i){if(i instanceof window.Element&&r.has(i)){var n=r.get(i);n.has(this)&&(n.delete(this),n.size||(r.delete(i),o.delete(i))),n.size||r.delete(i),o.size||cancelAnimationFrame(t)}},A.DOMRectReadOnly=c,A.ResizeObserver=s,A.ResizeObserverEntry=h,A}; // eslint-disable-line\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\", \"Left button pans, Right button zooms\\nx/y fixes axis, CTRL fixes aspect\", \"fa fa-arrows icon-move\", \"pan\"], [\"Zoom\", \"Zoom to rectangle\\nx/y fixes axis, CTRL fixes aspect\", \"fa fa-square-o icon-check-empty\", \"zoom\"], [\"\", \"\", \"\", \"\"], [\"Download\", \"Download plot\", \"fa fa-floppy-o icon-save\", \"download\"]];\n",
       "\n",
       "mpl.extensions = [\"eps\", \"jpeg\", \"pgf\", \"pdf\", \"png\", \"ps\", \"raw\", \"svg\", \"tif\"];\n",
       "\n",
       "mpl.default_extension = \"png\";/* global mpl */\n",
       "\n",
       "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.binaryType = comm.kernel.ws.binaryType;\n",
       "    ws.readyState = comm.kernel.ws.readyState;\n",
       "    function updateReadyState(_event) {\n",
       "        if (comm.kernel.ws) {\n",
       "            ws.readyState = comm.kernel.ws.readyState;\n",
       "        } else {\n",
       "            ws.readyState = 3; // Closed state.\n",
       "        }\n",
       "    }\n",
       "    comm.kernel.ws.addEventListener('open', updateReadyState);\n",
       "    comm.kernel.ws.addEventListener('close', updateReadyState);\n",
       "    comm.kernel.ws.addEventListener('error', updateReadyState);\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",
       "        var data = msg['content']['data'];\n",
       "        if (data['blob'] !== undefined) {\n",
       "            data = {\n",
       "                data: new Blob(msg['buffers'], { type: data['blob'] }),\n",
       "            };\n",
       "        }\n",
       "        // Pass the mpl event to the overridden (by mpl) onmessage function.\n",
       "        ws.onmessage(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 = document.getElementById(id);\n",
       "    var ws_proxy = comm_websocket_adapter(comm);\n",
       "\n",
       "    function ondownload(figure, _format) {\n",
       "        window.open(figure.canvas.toDataURL());\n",
       "    }\n",
       "\n",
       "    var fig = new mpl.figure(id, ws_proxy, ondownload, element);\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;\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",
       "    fig.cell_info[0].output_area.element.on(\n",
       "        'cleared',\n",
       "        { fig: fig },\n",
       "        fig._remove_fig_handler\n",
       "    );\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_close = function (fig, msg) {\n",
       "    var width = fig.canvas.width / fig.ratio;\n",
       "    fig.cell_info[0].output_area.element.off(\n",
       "        'cleared',\n",
       "        fig._remove_fig_handler\n",
       "    );\n",
       "    fig.resizeObserverInstance.unobserve(fig.canvas_div);\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.innerHTML =\n",
       "        '<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 / this.ratio;\n",
       "    var dataURL = this.canvas.toDataURL();\n",
       "    this.cell_info[1]['text/html'] =\n",
       "        '<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 () {\n",
       "        fig.push_to_output();\n",
       "    }, 1000);\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._init_toolbar = function () {\n",
       "    var fig = this;\n",
       "\n",
       "    var toolbar = document.createElement('div');\n",
       "    toolbar.classList = 'btn-toolbar';\n",
       "    this.root.appendChild(toolbar);\n",
       "\n",
       "    function on_click_closure(name) {\n",
       "        return function (_event) {\n",
       "            return fig.toolbar_button_onclick(name);\n",
       "        };\n",
       "    }\n",
       "\n",
       "    function on_mouseover_closure(tooltip) {\n",
       "        return function (event) {\n",
       "            if (!event.currentTarget.disabled) {\n",
       "                return fig.toolbar_button_onmouseover(tooltip);\n",
       "            }\n",
       "        };\n",
       "    }\n",
       "\n",
       "    fig.buttons = {};\n",
       "    var buttonGroup = document.createElement('div');\n",
       "    buttonGroup.classList = 'btn-group';\n",
       "    var button;\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",
       "            /* Instead of a spacer, we start a new button group. */\n",
       "            if (buttonGroup.hasChildNodes()) {\n",
       "                toolbar.appendChild(buttonGroup);\n",
       "            }\n",
       "            buttonGroup = document.createElement('div');\n",
       "            buttonGroup.classList = 'btn-group';\n",
       "            continue;\n",
       "        }\n",
       "\n",
       "        button = fig.buttons[name] = document.createElement('button');\n",
       "        button.classList = 'btn btn-default';\n",
       "        button.href = '#';\n",
       "        button.title = name;\n",
       "        button.innerHTML = '<i class=\"fa ' + image + ' fa-lg\"></i>';\n",
       "        button.addEventListener('click', on_click_closure(method_name));\n",
       "        button.addEventListener('mouseover', on_mouseover_closure(tooltip));\n",
       "        buttonGroup.appendChild(button);\n",
       "    }\n",
       "\n",
       "    if (buttonGroup.hasChildNodes()) {\n",
       "        toolbar.appendChild(buttonGroup);\n",
       "    }\n",
       "\n",
       "    // Add the status bar.\n",
       "    var status_bar = document.createElement('span');\n",
       "    status_bar.classList = 'mpl-message pull-right';\n",
       "    toolbar.appendChild(status_bar);\n",
       "    this.message = status_bar;\n",
       "\n",
       "    // Add the close button to the window.\n",
       "    var buttongrp = document.createElement('div');\n",
       "    buttongrp.classList = 'btn-group inline pull-right';\n",
       "    button = document.createElement('button');\n",
       "    button.classList = 'btn btn-mini btn-primary';\n",
       "    button.href = '#';\n",
       "    button.title = 'Stop Interaction';\n",
       "    button.innerHTML = '<i class=\"fa fa-power-off icon-remove icon-large\"></i>';\n",
       "    button.addEventListener('click', function (_evt) {\n",
       "        fig.handle_close(fig, {});\n",
       "    });\n",
       "    button.addEventListener(\n",
       "        'mouseover',\n",
       "        on_mouseover_closure('Stop Interaction')\n",
       "    );\n",
       "    buttongrp.appendChild(button);\n",
       "    var titlebar = this.root.querySelector('.ui-dialog-titlebar');\n",
       "    titlebar.insertBefore(buttongrp, titlebar.firstChild);\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._remove_fig_handler = function (event) {\n",
       "    var fig = event.data.fig;\n",
       "    if (event.target !== this) {\n",
       "        // Ignore bubbled events from children.\n",
       "        return;\n",
       "    }\n",
       "    fig.close_ws(fig, {});\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._root_extra_style = function (el) {\n",
       "    el.style.boxSizing = 'content-box'; // override notebook setting of border-box.\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._canvas_extra_style = function (el) {\n",
       "    // this is important to make the div 'focusable\n",
       "    el.setAttribute('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",
       "    } else {\n",
       "        // location in version 2\n",
       "        IPython.keyboard_manager.register_events(el);\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",
       "\n",
       "    // Check for shift+enter\n",
       "    if (event.shiftKey && event.which === 13) {\n",
       "        this.canvas_div.blur();\n",
       "        // select the cell after this one\n",
       "        var index = IPython.notebook.find_cell_index(this.cell_info[0]);\n",
       "        IPython.notebook.select(index + 1);\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_save = function (fig, _msg) {\n",
       "    fig.ondownload(fig, null);\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(\n",
       "        'matplotlib',\n",
       "        mpl.mpl_figure_comm\n",
       "    );\n",
       "}\n"
      ],
      "text/plain": [
       "<IPython.core.display.Javascript object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/html": [
       "<img src=\"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAoAAAAHgCAYAAAA10dzkAAAAAXNSR0IArs4c6QAAIABJREFUeF7t3Qm0lXW5P/CfmGJUOGOFIGVaN01wAEXN8Dan3Uuz1k3MhmsGpg1oliF3aWZ2zXAIGmm4luhKXQ2WWVfL7JJ5M0dIMpVEyjIIszSV/3re/3+f/+FwDmzOs895z9n78651Fgr72cPn/Z79PPud9mZr165dWywECBAgQIAAAQIdI7CZAbBj1rUXSoAAAQIECBCoBAyAgkCAAAECBAgQ6DABA2CHrXAvlwABAgQIECBgAJQBAgQIECBAgECHCRgAO2yFe7kECBAgQIAAAQOgDBAgQIAAAQIEOkzAANhhK9zLJUCAAAECBAgYAGWAAAECBAgQINBhAgbADlvhXi4BAgQIECBAwAAoAwQIECBAgACBDhMwAHbYCvdyCRAgQIAAAQIGQBkgQIAAAQIECHSYgAGww1a4l0uAAAECBAgQMADKAAECBAgQIECgwwQMgB22wr1cAgQIECBAgIABUAYIECBAgAABAh0mYADssBXu5RIgQIAAAQIEDIAyQIAAAQIECBDoMAEDYIetcC+XAAECBAgQIGAAlAECBAgQIECAQIcJGAA7bIV7uQQIECBAgAABA6AMECBAgAABAgQ6TMAA2GEr3MslQIAAAQIECBgAZYAAAQIECBAg0GECBsAOW+FeLgECBAgQIEDAACgDBAgQIECAAIEOEzAAdtgK93IJECBAgAABAgZAGSBAgAABAgQIdJiAAbDDVriXS4AAAQIECBAwAMoAAQIECBAgQKDDBAyAHbbCvVwCBAgQIECAgAFQBggQIECAAAECHSZgAOywFe7lEiBAgAABAgQMgDJAgAABAgQIEOgwAQNgh61wL5cAAQIECBAgYACUAQIEahe4/vrry1VXXVVOOOGEss022wz681m1alWZPXt2ueyyy8rDDz9cpkyZUv7zP/+z7LPPPht8Lk888UT5yle+Ur75zW+WX/7yl+XBBx8sz3rWs8oRRxxRPvCBD5StttpqvfovfOEL5ZOf/GT57W9/W8aNG1eOP/74MmvWrPVud99995UTTzyxconHOfTQQ8unPvWp8uxnP3vQfTwgAQLtJ2AAbL916hURGHYCMRB98IMfrIaiCRMmDOrzj+HqhS98YfnVr35VPYcddtihXHjhhWX58uXlxhtvLLvttlufz+ehhx4qT3va08oBBxxQDj/88DJmzJjys5/9rHz5y18uhxxySPnRj35UNttss676BQsWlGOPPba87nWvKy9/+cvLT37yk/LVr361fPzjHy8nnXRS1+3ifmP4XL16dXn/+99ftthii2r4W7t2bbnpppvK9ttvP6hGHowAgfYTMAC23zr1iggMO4E6B8BFixaVN73pTeWSSy4pr3/96yu7Bx54oOy+++7lla98Zbnooov69Hz00UfLL37xi3LggQeuc5v/+I//KHPmzCk/+MEPykte8pLq3/72t79VW/xiWPz2t7/ddft/+7d/K5dffnk1cG677bbV33/iE5+oBsKf//znZfLkydXfLVmypOy5557VlsqPfexjw24de8IECAwtAQPg0Fofng2BjhM47bTTyty5c9d73YO1NfCNb3xj+fGPf1xWrFhRRowY0fU8/v3f/7187Wtfq3brjhw5cpPWyy233FL22muvMm/evK7du9/97nfLYYcdVr7zne+UV73qVV33F1sMY4CMLYExDMYSu6BjiQGw+xJbDX/zm9+UZcuWbdLzcWMCBAj0FDAAygQBArUK3HzzzdUu0K9//evVbs7YBRvLa17zmvKUpzyl1+cWx+nFz8aWzTffvGurWl+3jV288RMDWvcljtV7xzveUeL5veAFL9jYQ63z77Hl72Uve1m19fDII4+s/u2MM84oH/nIR8rvf//7aldxY4mtiE9+8pOr4x/juMPYJT1q1KhyzDHHVLuiuy+nnnpqOf3008tf/vKXatezhQABAv0VMAD2V04dAQItE9jUXcB9bTXs+YR22WWXcvfdd2/weT71qU+tdgHHwNd9aWyx+973vlcdr7cpy0tf+tJq690999zTdVLLzJkzy/z588tjjz223l3FQPjiF7+4GoL/+Mc/lh133LHEbuQY+LovMRC+5z3vqXYHP/e5z92Up+S2BAgQWEfAACgQBAjULrCpA+Bdd91V4mdjS2xZO+iggzZ4s9hKGLt7e25tixM4YiiLM4OnT5++sYfq+vc4Pu/DH/5wdX/vfve7u/7+7W9/ezXg9bblcvz48dVJH41jAeP/zzrrrOp4v+7LF7/4xRL3E2ccT5o0qenn5IYECBDoKWAAlAkCBGoX2NQBsJVPuJVbAC+++OJql2/svv385z+/ztO0BbCVa819ESCQFTAAZgXVEyCQFtjUATAukxI/G1ti617sTt3Q0qpjAOO4v7gUTBz7F1sNn/SkJ63zsI4B3Nja8u8ECAymgAFwMLU9FgECvQrEyQ9x4eRmz/xt5TGAb3jDG6rr8fU8C/hd73pX+a//+q+mzgJevHhxtbt44sSJ5eqrr65O6ui5xNm/MSD2PAs4LoIdu6njgtJvfetbq7K49EtcP7DnWcAxXMZZwPFjIUCAQEbAAJjRU0uAQEsE4uSIOF6u2WPbWnkMYOy2jW/u6H4dwDgRI7YMxskf3/jGN7peY2Pw2nXXXbv+7o477qguJP30pz+9GiQb1/LrCRPXAdx5552rS75861vf6vrnGPrim0TiOoDbbbdd9fdx/N/JJ59cbrjhhrLffvtVf7d06dKyxx57VINynDVtIUCAQEbAAJjRU0uAQEsEYtCJa9/F9fFiGItvvnj1q1/d52VgWvKg/+9OHn/88XLwwQeXW2+9dZ1vArn33nurAaz72baNbylpnFm8Zs2aaiiLr22Lkz/Gjh27zlOLQXHq1Kldf9c4izcuON34JpDY8he7h0855ZSu28X97r333iX+jIEvPM4555wSzzW+CWRju7Vb6eO+CBBoTwEDYHuuV6+KwLATiOvbxZbA+++/v7oWXrO7g1vxQv/85z9Xw1+chRtb6mIXbByX2Nj61niMngNgDILx3b99LTNmzCgLFy5c558/97nPVdf7a3wXcJwc8t73vnedr4yLgt/97nfrfBfwtGnTquskPuc5z2nFS3YfBAh0uIABsMMD4OUTIECAAAECnSdgAOy8de4VEyBAgAABAh0uYADs8AB4+QQIECBAgEDnCRgAO2+de8UECBAgQIBAhwu0zQD44x//uJx99tnlxhtvrA4ib+brm6655pryvve9r9x2221l3Lhx1Re1H3300UMiEo888kg588wzy4c+9KEycuTIIfGc2vFJcB6ctcqZ8+AIDM6jyDPnwREY2EdpmwHwyiuvLD/96U/LvvvuW1772tdudACMM/D23HPPcuyxx5Z3vOMd5Yc//GE54YQTqou0buoXvw/EKvrLX/5Stt5667J69eoyevTogXgI91lK4Tw4MeDMeXAEBudR5Jnz4AgM7KO0zQDYnSmuoL+xLYAnnXRSNezFtb8aS1x/bNWqVeV73/vewKo3ce/eYJpAasFNOLcAsYm74NwEUgtuwrkFiE3cBecmkFpwE84tQNzAXXTsAHjIIYeUffbZp5x77rldPF/60peqrYCx1a3uRfAHZw1w5jw4AoPzKPLMeXAEBudR5HlgnTt2ANx9993L2972tuoYu8by3e9+txx22GHl4Ycf7vW7POO4j/hpLHGx2gcffLBsv/32613ENbvaIvhxXGJ8PZRdwFnNvus5D5xt93vmzHlwBAbnUeR5+DuvXbu2/OlPf6q+fnHEiBFdLyiOue+U4+4NgJswADb7BfSD86vhUQgQIECAAIFWCsyZM6dEr++EpWMHwP7sAu65BTB2FY8fP95Wuk74TfEaCRAgQKBtBBpbceM7v+OEy8ZiC+AwX8XNngQSu3xvueWWrlf75je/udql2+xJII5PGOZB8fQJECBAoCMF9O9S2mYL4EMPPVSWLVtWBXnvvfcu55xzTjn00EOr/fuxlS6O9bvvvvvKV77yleo2jcvAvOc97ynHHHNM+dGPflSOP/74TboMjAB15PuGF02AAAECw1xA/26jATAu6hwDX89lxowZZeHChdUFnu++++4St2ss8d8nnnhiuf3228vOO+9cTj311E26ELQADfN3AE+fAAECBDpSQP9uowGwjgQLUB3qHpMAAQIECOQE9G8DYCpBApTiU0yAAAECBGoR0L8NgKngCVCKTzEBAgQIEKhFQP82AKaCJ0ApPsUECBAgQKAWAf3bAJgKngCl+BQTIECAAIFaBPRvA2AqeAKU4lNMgAABAgRqEdC/DYCp4AlQik8xAQIECBCoRUD/NgCmgidAKT7FBAgQIECgFgH92wCYCp4ApfgUEyBAgACBWgT0bwNgKngClOJTTIAAAQIEahHQvw2AqeAJUIpPMQECBAgQqEVA/zYApoInQCk+xQQIECBAoBYB/dsAmAqeAKX4FBMgQIAAgVoE9G8DYCp4ApTiU0yAAAECBGoR0L8NgKngCVCKTzEBAgQIEKhFQP82AKaCJ0ApPsUECBAgQKAWAf3bAJgKngCl+BQTIECAAIFaBPRvA2AqeAKU4lNMgAABAgRqEdC/DYCp4AlQik8xAQIECBCoRUD/NgCmgidAKT7FBAgQIECgFgH92wCYCp4ApfgUEyBAgACBWgT0bwNgKngClOJTTIAAAQIEahHQvw2AqeAJUIpPMQECBAgQqEVA/zYApoInQCk+xQQIECBAoBYB/dsAmAqeAKX4FBMgQIAAgVoE9G8DYCp4ApTiU0yAAAECBGoR0L8NgKngCVCKTzEBAgQIEKhFQP82AKaCJ0ApPsUECBAgQKAWAf3bAJgKngCl+BQTIECAAIFaBPRvA2AqeAKU4lNMgAABAgRqEdC/DYCp4AlQik8xAQIECBCoRUD/NgCmgidAKT7FBAgQIECgFgH92wCYCp4ApfgUEyBAgACBWgT0bwNgKngClOJTTIAAAQIEahHQvw2AqeAJUIpPMQECBAgQqEVA/zYApoInQCk+xQQIECBAoBYB/dsAmAqeAKX4FBMgQIAAgVoE9G8DYCp4ApTiU0yAAAECBGoR0L8NgKngCVCKTzEBAgQIEKhFQP82AKaCJ0ApPsUECBAgQKAWAf3bAJgKngCl+BQTIECAAIFaBPRvA2AqeAKU4lNMgAABAgRqEdC/DYCp4AlQik8xAQIECBCoRUD/NgCmgidAKT7FBAgQIECgFgH92wCYCp4ApfgUEyBAgACBWgT0bwNgKngClOJTTIAAAQIEahHQvw2AqeAJUIpPMQECBAgQqEVA/zYApoInQCk+xQQIECBAoBYB/dsAmAqeAKX4FBMgQIAAgVoE9G8DYCp4ApTiU0yAAAECBGoR0L8NgKngCVCKTzEBAgQIEKhFQP82AKaCJ0ApPsUECBAgQKAWAf3bAJgKngCl+BQTIECAAIFaBPRvA2AqeAKU4lNMgAABAgRqEdC/DYCp4AlQik8xAQIECBCoRUD/NgCmgidAKT7FBAgQIECgFgH92wCYCp4ApfgUEyBAgACBWgT0bwNgKngClOJTTIAAAQIEahHQvw2AqeAJUIpPMQECBAgQqEVA/zYApoInQCk+xQQIECBAoBYB/dsAmAqeAKX4FBMgQIAAgVoE9G8DYCp4ApTiU0yAAAECBGoR0L8NgKngCVCKTzEBAgQIEKhFQP82AKaCJ0ApPsUECBAgQKAWAf27zQbACy64oJx99tll5cqVZeLEieW8884rU6ZM6TNc5557bvnMZz5T7r333rLDDjuU17/+9eXMM88sW221VVOBFKCmmNyIAAECBAgMKQH9u40GwIsvvrgcddRRZf78+WX//fcvMdxdcsklZenSpWXMmDHrBe+iiy4qxxxzTPniF79YDjzwwPLrX/+6HH300eWII44o55xzTlNBFaCmmNyIAAECBAgMKQH9u40GwBj6Jk+eXM4///wqZE888UQZN25cmTVrVjn55JPXC97MmTPLHXfcUX74wx92/dv73//+snjx4nLdddc1FVQBaorJjQgQIECAwJAS0L/bZAB89NFHy6hRo8qll15apk+f3hWyGTNmlFWrVpUrrrii1y2Axx13XLnqqquq3cR33XVXOeyww8pb3/rWcsoppzQVVAFqismNCBAgQIDAkBLQv9tkAFyxYkUZO3Zsuf7668vUqVO7QjZ79uxy7bXXVlv1elvmzZtXPvCBD5S1a9eWxx57rBx77LHVMYF9LY888kiJn8YSAYqtjKtXry6jR48eUuH2ZAgQIECAAIHeBRoD4PLly9fp3yNHjizx0wnLZmtj+hnmS38GwGuuuaY63u/000+vjhlctmxZee9731ve+c53llNPPbVXkdNOO63MnTt3vX8zAA7zAHn6BAgQINBRAo0BsOeLnjNnTole3wlLWwyA/dkF/MIXvrAccMAB1VnDjeVrX/taede73lUeeuihMmLEiPXWvy2AnfAr4TUSIECAQLsL2ALYJruAI6ixFS+O5YtLv8QSJ4GMHz++xMkevZ0Esu+++5aXvOQl5ayzzurK+de//vXy9re/vaxZs6ZsvvnmG82/Ywg2SuQGBAgQIEBgyAno3200AMZlYOKkjwULFlSDYFwGZtGiRWXJkiVlp512qi4RE8cJxnX+YolNvHG5l89+9rNdu4Df/e53lxgM476aWQSoGSW3IUCAAAECQ0tA/26jATCiFZeAaVwIetKkSSVO8ogtg7FMmzatTJgwoSxcuLD6/zjp44wzzihf/epXy3333Vd23HHH8upXv7r6u2222aappApQU0xuRIAAAQIEhpSA/t1mA+Bgp0uABlvc4xEgQIAAgbyA/m0ATKVIgFJ8igkQIECAQC0C+rcBMBU8AUrxKSZAgAABArUI6N8GwFTwBCjFp5gAAQIECNQioH8bAFPBE6AUn2ICBAgQIFCLgP5tAEwFT4BSfIoJECBAgEAtAvq3ATAVPAFK8SkmQIAAAQK1COjfBsBU8AQoxaeYAAECBAjUIqB/GwBTwROgFJ9iAgQIECBQi4D+bQBMBU+AUnyKCRAgQIBALQL6twEwFTwBSvEpJkCAAAECtQjo3wbAVPAEKMWnmAABAgQI1CKgfxsAU8EToBSfYgIECBAgUIuA/m0ATAVPgFJ8igkQIECAQC0C+rcBMBU8AUrxKSZAgAABArUI6N8GwFTwBCjFp5gAAQIECNQioH8bAFPBE6AUn2ICBAgQIFCLgP5tAEwFT4BSfIoJECBAgEAtAvq3ATAVPAFK8SkmQIAAAQK1COjfBsBU8AQoxaeYAAECBAjUIqB/GwBTwROgFJ9iAgQIECBQi4D+bQBMBU+AUnyKCRAgQIBALQL6twEwFTwBSvEpJkCAAAECtQjo3wbAVPAEKMWnmAABAgQI1CKgfxsAU8EToBSfYgIECBAgUIuA/m0ATAVPgFJ8igkQIECAQC0C+rcBMBU8AUrxKSZAgAABArUI6N8GwFTwBCjFp5gAAQIECNQioH8bAFPBE6AUn2ICBAgQIFCLgP5tAEwFT4BSfIoJECBAgEAtAvq3ATAVPAFK8SkmQIAAAQK1COjfBsBU8AQoxaeYAAECBAjUIqB/GwBTwROgFJ9iAgQIECBQi4D+bQBMBU+AUnyKCRAgQIBALQL6twEwFTwBSvEpJkCAAAECtQjo3wbAVPAEKMWnmAABAgQI1CKgfxsAU8EToBSfYgIECBAgUIuA/m0ATAVPgFJ8igkQIECAQC0C+rcBMBU8AUrxKSZAgAABArUI6N8GwFTwBCjFp5gAAQIECNQioH8bAFPBE6AUn2ICBAgQIFCLgP5tAEwFT4BSfIoJECBAgEAtAvq3ATAVPAFK8SkmQIAAAQK1COjfBsBU8AQoxaeYAAECBAjUIqB/GwBTwROgFJ9iAgQIECBQi4D+bQBMBU+AUnyKCRAgQIBALQL6twEwFTwBSvEpJkCAAAECtQjo3wbAVPAEKMWnmAABAgQI1CKgfxsAU8EToBSfYgIECBAgUIuA/m0ATAVPgFJ8igkQIECAQC0C+rcBMBU8AUrxKSZAgAABArUI6N8GwFTwBCjFp5gAAQIECNQioH8bAFPBE6AUn2ICBAgQIFCLgP5tAEwFT4BSfIoJECBAgEAtAvq3ATAVPAFK8SkmQIAAAQK1COjfBsBU8AQoxaeYAAECBAjUIqB/GwBTwROgFJ9iAgQIECBQi4D+bQBMBU+AUnyKCRAgQIBALQL6twEwFTwBSvEpJkCAAAECtQjo3wbAVPAEKMWnmAABAgQI1CKgfxsAU8EToBSfYgIECBAgUIuA/m0ATAVPgFJ8igkQIECAQC0C+rcBMBU8AUrxKSZAgAABArUI6N8GwFTwBCjFp5gAAQIECNQioH+32QB4wQUXlLPPPrusXLmyTJw4sZx33nllypQpfYZr1apV5cMf/nD55je/WR588MGyyy67lHPPPbe86lWvaiqQAtQUkxsRIECAAIEhJaB/t9EAePHFF5ejjjqqzJ8/v+y///7VIHfJJZeUpUuXljFjxqwXvEcffbQcdNBB1b+dcsopZezYseWee+4p22yzTTU8NrMIUDNKbkOAAAECBIaWgP7dRgNgDH2TJ08u559/fpWyJ554oowbN67MmjWrnHzyyeslLwbF2Fq4ZMmSssUWW/QrmQLULzZFBAgQIECgVgH9u00GwNiaN2rUqHLppZeW6dOnd4VqxowZJXbzXnHFFesFLXbzbrfddlVd/PuOO+5Y3vzmN5eTTjqpbL755r0G85FHHinx01giQDFkrl69uowePbrWMHtwAgQIECBAoDmBxgC4fPnydfr3yJEjS/x0wrLZ2rVr1w73F7pixYpqF+71119fpk6d2vVyZs+eXa699tqyePHi9V7i8573vHL33XeXt7zlLeW4444ry5Ytq/48/vjjy5w5c3olOe2008rcuXPX+zcD4HBPkOdPgAABAp0k0BgAe77m6P/R6zth6dgBcPfddy9///vfy29/+9uuLX7nnHNOtVv4/vvvtwWwE9LvNRIgQIBARwrYAtjBu4Bf9KIXVcf+XX311V3hv/LKK6szgGM375ZbbrnRXwrHEGyUyA0IECBAgMCQE9C/22QAjGTFSSBxyZe49EsscRLI+PHjy8yZM3s9CSTO/L3ooovKXXfdVUaMGFHVfPrTny5nnXVWiV3KzSwC1IyS2xAgQIAAgaEloH+30QAYl4GJkz4WLFhQDYJxGZhFixZVZ/nutNNO1SVi4jjBM888s0phHPi5xx57VDVxpvCdd95ZjjnmmOoYwLg2YDOLADWj5DYECBAgQGBoCejfbTQARrTiEjCNC0FPmjSpzJs3r9oyGMu0adPKhAkTysKFC7tS+LOf/ayceOKJ5aabbqqGw7e//e0bPAu4Z3wFaGj9Qns2BAgQIECgGQH9u80GwGZWeitvI0Ct1HRfBAgQIEBgcAT0bwNgKmkClOJTTIAAAQIEahHQvw2AqeAJUIpPMQECBAgQqEVA/zYApoInQCk+xQQIECBAoBYB/dsAmAqeAKX4FBMgQIAAgVoE9G8DYCp4ApTiU0yAAAECBGoR0L8NgKngCVCKTzEBAgQIEKhFQP82AKaCJ0ApPsUECBAgQKAWAf3bAJgKngCl+BQTIECAAIFaBPRvA2AqeAKU4lNMgAABAgRqEdC/DYCp4AlQik8xAQIECBCoRUD/NgCmgidAKT7FBAgQIECgFgH92wCYCp4ApfgUEyBAgACBWgT0bwNgKngClOJTTIAAAQIEahHQvw2AqeAJUIpPMQECBAgQqEVA/zYApoInQCk+xQQIECBAoBYB/dsAmAqeAKX4FBMgQIAAgVoE9G8DYCp4ApTiU0yAAAECBGoR0L8NgKngCVCKTzEBAgQIEKhFQP82AKaCJ0ApPsUECBAgQKAWAf3bAJgKngCl+BQTIECAAIFaBPRvA2AqeAKU4lNMgAABAgRqEdC/DYCp4AlQik8xAQIECBCoRUD/NgCmgidAKT7FBAgQIECgFgH92wCYCp4ApfgUEyBAgACBWgT0bwNgKngClOJTTIAAAQIEahHQvw2AqeAJUIpPMQECBAgQqEVA/zYApoInQCk+xQQIECBAoBYB/dsAmAqeAKX4FBMgQIAAgVoE9G8DYCp4ApTiU0yAAAECBGoR0L8NgKngCVCKTzEBAgQIEKhFQP82AKaCJ0ApPsUECBAgQKAWAf3bAJgKngCl+BQTIECAAIFaBPRvA2AqeAKU4lNMgAABAgRqEdC/DYCp4AlQik8xAQIECBCoRUD/NgCmgidAKT7FBAgQIECgFgH92wCYCp4ApfgUEyBAgACBWgT0bwNgKngClOJTTIAAAQIEahHQvw2AqeAJUIpPMQECBAgQqEVA/zYApoInQCk+xQQIECBAoBYB/dsAmAqeAKX4FBMgQIAAgVoE9G8DYCp4ApTiU0yAAAECBGoR0L8NgKngCVCKTzEBAgQIEKhFQP82AKaCJ0ApPsUECBAgQKAWAf3bAJgKngCl+BQTIECAAIFaBPRvA2AqeAKU4lNMgAABAgRqEdC/DYCp4AlQik8xAQIECBCoRUD/NgCmgidAKT7FBAgQIECgFgH92wCYCp4ApfgUEyBAgACBWgT0bwNgKngClOJTTIAAAQIEahHQvw2AqeAJUIpPMQECBAgQqEVA/zYApoInQCk+xQQIECBAoBYB/dsAmAqeAKX4FBMgQIAAgVoE9G8DYCp4ApTiU0yAAAECBGoR0L8NgKngCVCKTzEBAgQIEKhFQP82AKaCJ0ApPsUECBAgQKAWAf3bAJgKngCl+BQTIECAAIFaBPRvA2AqeAKU4lNMgAABAgRqEdC/DYCp4AlQik8xAQIECBCoRUD/NgCmgidAKT7FBAgQIECgFgH92wCYCp4ApfgUEyBAgACBWgT0bwNgKngClOJTTIAAAQIEahHQvw2AqeAJUIpPMQECBAggDehWAAAgAElEQVQQqEVA/zYApoInQCk+xQQIECBAoBYB/bvNBsALLrignH322WXlypVl4sSJ5bzzzitTpkzZaLi+8Y1vlCOPPLL867/+a7n88ss3evvGDQSoaSo3JECAAAECQ0ZA/26jAfDiiy8uRx11VJk/f37Zf//9y7nnnlsuueSSsnTp0jJmzJg+Q3f33XeXgw8+uDz72c8u2223nQFwyPx6eiIECBAgQGBgBAyAbTQAxtA3efLkcv7551dpeeKJJ8q4cePKrFmzysknn9xrgh5//PFyyCGHlGOOOab85Cc/KatWrTIADszvmnslQIAAAQJDRsAA2CYD4KOPPlpGjRpVLr300jJ9+vSugM2YMaMa6q644opeQzdnzpxy8803l8suu6wcffTRGx0AH3nkkRI/3XcBx5C5evXqMnr06CETbE+EAAECBAgQ6FugMQAuX758nf49cuTIEj+dsGy2du3atcP9ha5YsaKMHTu2XH/99WXq1KldL2f27Nnl2muvLYsXL17vJV533XXliCOOKDfddFPZYYcdmhoATzvttDJ37tz17ssAONwT5PkTIECAQCcJNAbAnq85NgxFr++EpSMHwDVr1pS99tqrXHjhheWVr3xltZ5tAeyEuHuNBAgQIECgFFsAO3QXcGz123vvvcvmm2/e9XsQxwzGMmLEiOrEkV133XWjvyOOIdgokRsQIECAAIEhJ6B/t8kAGMmKk0Diki9x6ZdYYqAbP358mTlz5nongfz9738vy5YtWyeQH/nIR0psGfz0pz9ddt9997LllltuNLACtFEiNyBAgAABAkNOQP9uowEwLgMTJ30sWLCgGgTjMjCLFi0qS5YsKTvttFN1iZg4TvDMM8/sNYjN7ALuWShAQ+532hMiQIAAAQIbFdC/22gAjLUdl4BpXAh60qRJZd68edWWwVimTZtWJkyYUBYuXGgA3OivhhsQIECAAIH2FTAAttkAONhRFaDBFvd4BAgQIEAgL6B/GwBTKRKgFJ9iAgQIECBQi4D+bQBMBU+AUnyKCRAgQIBALQL6twEwFTwBSvEpJkCAAAECtQjo3wbAVPAEKMWnmAABAgQI1CKgfxsAU8EToBSfYgIECBAgUIuA/m0ATAVPgFJ8igkQIECAQC0C+rcBMBU8AUrxKSZAgAABArUI6N8GwFTwBCjFp5gAAQIECNQioH8bAFPBE6AUn2ICBAgQIFCLgP5tAEwFT4BSfIoJECBAgEAtAvq3ATAVPAFK8SkmQIAAAQK1COjfBsBU8AQoxaeYAAECBAjUIqB/GwBTwROgFJ9iAgQIECBQi4D+bQBMBU+AUnyKCRAgQIBALQL6twEwFTwBSvEpJkCAAAECtQjo3wbAVPAEKMWnmAABAgQI1CKgfxsAU8EToBSfYgIECBAgUIuA/m0ATAVPgFJ8igkQIECAQC0C+rcBMBU8AUrxKSZAgAABArUI6N8GwFTwBCjFp5gAAQIECNQioH8bAFPBE6AUn2ICBAgQIFCLgP5tAEwFT4BSfIoJECBAgEAtAvq3ATAVPAFK8SkmQIAAAQK1COjfBsBU8AQoxaeYAAECBAjUIqB/GwBTwROgFJ9iAgQIECBQi4D+bQBMBU+AUnyKCRAgQIBALQL6twEwFTwBSvEpJkCAAAECtQjo3wbAVPAEKMWnmAABAgQI1CKgfxsAU8EToBSfYgIECBAgUIuA/m0ATAVPgFJ8igkQIECAQC0C+rcBMBU8AUrxKSZAgAABArUI6N8GwFTwBCjFp5gAAQIECNQioH8bAFPBE6AUn2ICBAgQIFCLgP5tAEwFT4BSfIoJECBAgEAtAvq3ATAVPAFK8SkmQIAAAQK1COjfBsBU8AQoxaeYAAECBAjUIqB/GwBTwROgFJ9iAgQIECBQi4D+bQBMBU+AUnyKCRAgQIBALQL6twEwFTwBSvEpJkCAAAECtQjo3wbAVPAEKMWnmAABAgQI1CKgfxsAU8EToBSfYgIECBAgUIuA/m0ATAVPgFJ8igkQIECAQC0C+rcBMBU8AUrxKSZAgAABArUI6N8GwFTwBCjFp5gAAQIECNQioH8bAFPBE6AUn2ICBAgQIFCLgP5tAEwFT4BSfIoJECBAgEAtAvq3ATAVPAFK8SkmQIAAAQK1COjfBsBU8AQoxaeYAAECBAjUIqB/GwBTwROgFJ9iAgQIECBQi4D+bQBMBU+AUnyKCRAgQIBALQL6twEwFTwBSvEpJkCAAAECtQjo3wbAVPAEKMWnmAABAgQI1CKgfxsAU8EToBSfYgIECBAgUIuA/m0ATAVPgFJ8igkQIECAQC0C+rcBMBU8AUrxKSZAgAABArUI6N8GwFTwBCjFp5gAAQIECNQioH8bAFPBE6AUn2ICBAgQIFCLgP5tAEwFT4BSfIoJECBAgEAtAvq3ATAVPAFK8SkmQIAAAQK1COjfBsBU8AQoxaeYAAECBAjUIqB/GwBTwROgFJ9iAgQIECBQi4D+bQBMBU+AUnyKCRAgQIBALQL6twEwFTwBSvEpJkCAAAECtQjo3202AF5wwQXl7LPPLitXriwTJ04s5513XpkyZUqv4frc5z5XvvKVr5Rbb721+vd99923fOxjH+vz9r3diQDV8nvrQQkQIECAQEpA/26jAfDiiy8uRx11VJk/f37Zf//9y7nnnlsuueSSsnTp0jJmzJj1gvKWt7ylHHTQQeXAAw8sW221VTnrrLPKZZddVm677bYyduzYpoIlQE0xuREBAgQIEBhSAvp3Gw2AMfRNnjy5nH/++VXInnjiiTJu3Lgya9ascvLJJ280eI8//njZdtttq/oYJJtZBKgZJbchQIAAAQJDS0D/bpMB8NFHHy2jRo0ql156aZk+fXpXymbMmFFWrVpVrrjiio0mb82aNdWWwthqePjhh/d6+0ceeaTET2OJAMWQuXr16jJ69OiNPoYbECBAgAABAvULNAbA5cuXr9O/R44cWeKnE5bN1q5du3a4v9AVK1ZUu22vv/76MnXq1K6XM3v27HLttdeWxYsXb/QlHnfcceX73/9+tQs4dgn3tpx22mll7ty56/2TAXCjvG5AgAABAgSGjEBjAOz5hObMmVOi13fCYgAspXz84x8vn/jEJ8o111xT9tprrz7Xuy2AnfAr4TUSIECAQLsL2AJoF3D55Cc/WU4//fRy9dVXl/3222+TMu8Ygk3icmMCBAgQIDAkBPTvNhkAI01xEkhc8iUu/RJLnAQyfvz4MnPmzD5PAomtfmeccUa16/eAAw7Y5FAK0CaTKSBAgAABArUL6N9tNADGZWDipI8FCxZUg2BcBmbRokVlyZIlZaeddqrO7I3jBM8888wqeHHZl49+9KPloosuqi4H01ie+tSnlvhpZhGgZpTchgABAgQIDC0B/buNBsCIVlzCpXEh6EmTJpV58+ZVWwZjmTZtWpkwYUJZuHBh9f/x3/fcc896idyUA0AFaGj9Qns2BAgQIECgGQH9u80GwGZWeitvI0Ct1HRfBAgQIEBgcAT0bwNgKmkClOJTTIAAAQIEahHQvw2AqeAJUIpPMQECBAgQqEVA/zYApoInQCk+xQQIECBAoBYB/dsAmAqeAKX4FBMgQIAAgVoE9G8DYCp4ApTiU0yAAAECBGoR0L8NgKngCVCKTzEBAgQIEKhFQP82AKaCJ0ApPsUECBAgQKAWAf3bAJgKngCl+BQTIECAAIFaBPRvA2AqeAKU4lNMgAABAgRqEdC/DYCp4AlQik8xAQIECBCoRUD/NgCmgidAKT7FBAgQIECgFgH92wCYCp4ApfgUEyBAgACBWgT0bwNgKngClOJTTIAAAQIEahHQvw2AqeAJUIpPMQECBAgQqEVA/zYApoInQCk+xQQIECBAoBYB/dsAmAqeAKX4FBMgQIAAgVoE9G8DYCp4ApTiU0yAAAECBGoR0L8NgKngCVCKTzEBAgQIEKhFQP82AKaCJ0ApPsUECBAgQKAWAf3bAJgKngCl+BQTIECAAIFaBPRvA2AqeAKU4lNMgAABAgRqEdC/DYCp4AlQik8xAQIECBCoRUD/NgCmgidAKT7FBAgQIECgFgH92wCYCp4ApfgUE6hNYO3ateWRRx4pf/3rX8tDDz203k/j7+MJPvWpTy1PecpTqj97/sTfjxw5smy22Wa1vRYPTIDApgvo3wbATU9NtwoBSvEpJtCUwGOPPbbeoLaxwa3nUNfb7R9//PGmHn9jN9p88817HQz7Ghib+fsYLJ/0pCdt7KH9OwEC/RTQvw2A/YzO/y0ToBSf4g4Q+Pvf/17uvvvuctddd5U//vGP1Za23oaxDQ10cR8DuWy11VZ9bt2Lx+3ruQ3G89rQlse+BskddtihPPvZzy4TJkwo8dosBAisL6B/GwBTvxcClOJT3AYCsSv1gQceqAa83/zmN9Wf3f/7vvvua9mrjC1tT3va0za6S7a3waiv3bgxYMX99meJLYit2BLZfShes2ZNadWWyXhNY8eOLbvuums1EMZP9//ecccd7bruz4pX0xYC+rcBMBVkAUrxKR4mAo8++mjXVryeg178fwwwG1pi+IrBY8yYMV0DXPchrdmtXFtuuWXbDywxUId3913YmzJkxgD5hz/8oRrGm1kvvQ2Gsa522WWXEt4WAu0qoH8bAFPZFqAUn+IhIhBDx4MPPtjnVrzly5eXuE1fS5wAsaEtTbFL0kkSg7uyY33FLveeW2QbA3xsmd3QOh0xYkTZeeed+9x6uN1221mng7tKPVqLBfRvA2AqUgKU4lM8iAL/+Mc/yr333tvnQBBZ3tAyatSo9XYhNnYpxtYix5oN4spswUPF8Yv33HNPr7vtY0h8+OGHN/goW2+9da+7lWPr4bhx48oWW2zRgmfpLggMnID+bQBMpUuAUnyKWyywatWqPrfixfC3sWPLnvnMZ/bZ1GP3ra14LV5hQ/TuYsvg73//+z4/LNx///0bfOZxTOX48eP73Hq4zTbbDNFX7ml1koD+bQBM5V2AUnyK+yEQx3j98pe/LL/+9a/X23rz5z//eYP3GFvpnvWsZ/XamOOM0djKZyGwMYHYOhhndnc/6af7saFxfcUNLbH7uOexh7vvvnvZe++9q2NELQQGQ0D/NgCmciZAKT7FGxGIkwFuvvnmcsMNN5Sf//zn1Z+33377Bo/d2mmnnfrcivf0pz+9xLFdFgIDJfDEE0+UlStX9rlrObYs9rXEFubnP//5ZfLkyWXKlCnVn3vttZeTUQZqZXX4/erfBsDUr4AApfgUdxOIxhlb9RqDXvx50003VWeE9lzi4PwXvOAF6w16sXUvzq61EBiqAnFm8m9/+9v1DlW45ZZbyu9+97v1nnaciTxp0qSugTAGw9ha6IPMUF3Dw+d56d8GwFRaBSjF17HFcYxVNLvuW/Z+8YtfVBcW77lsu+226zS/2CoSW/IsBNpNILYcdv+diA9BvR3WMHr06LLffvuts6UwPhQ5RrXdEjGwr0f/NgCmEiZAKb6OKY5LrPRsbL3tCnvyk59c9t1336qxNXaDxbFSGlvHRMUL7SYQH5Ti2MLuW8X/93//t/ztb39bzyk+FHXfdRz/HccaWgj0JaB/GwBTvx0ClOJry+I4QD6aVPeBLw6W77nEmZKxG7dxrFP8Gcc/+f7XtoyFF9Uigfhe6Ntuu22d36/YfdzbGe5xSZruv19xkokTnVq0ItrgbvRvA2AqxgKU4hv2xXFtvWhG3bdQxP/31ox22223dbZQxHFNmtGwj4AXMAQE4kNXHC/b/UPXnXfe2euHrj333HOd38M99tjDNQuHwDqs4yno3wbAVO4EKMU3rIpjd9SyZcvWaTKxpS8uqNtzecYznrHOloc4XimO5bMQIDA4AnHsYBxX2/3DWW/XL4xLI+2zzz7rDIXPec5zHHYxOKup1kfRvw2AqQAKUIpvSBdHs+jePKKZ9HZAenwjQuOYvcYxSPG1aBYCBIaWQHz9Xc/f6dWrV6/3JOPDWnxo6777OD7UWdpLQP82AKYSLUApviFTHE2g59aCaBY9l5EjR1YXq+3eGGJrgUtSDJlV6YkQaFogLr0Uu4q77zqOi6z3diHr+FDX/fc+BsT48GcZvgL6twEwlV4BSvHVVhxb8v77v/+7/PCHPyw/+tGPypIlS9Z7LjHUxfFB3c8sjJM2fMdpbavNAxMYcIG47uatt966zlAYx/XGsNhzed7znlf++Z//ubz4xS8uhx56qMM8BnzttPYB9G8DYCpRApTiG7TiuGzET3/603L11VdXQ9+NN9643rdpxEWUe54x6KLKg7aKPBCBISsQF6+OLYPddx/Hxay7L/GBMY4ljGHwJS95STnooINKXNbJMnQF9G8DYCqdApTiG7DiuFREDHkx7MXQd/3116+3W+ef/umfqjfqeMOON+sddthhwJ6POyZAoL0E/vjHP5brrruueo+JnzvuuGOdFxiHi8T7Sry/xE9c39MlnoZWBvRvA2AqkQKU4mtZcZyhG7txG1v4rrnmmtLz4O74poDGp/PYbfPMZz6zZY/vjggQ6GyBFStWdA2D8T7U8xjiOF5w2rRpXR86Y/exC7zXmxn92wCYSqAApfhSxfFVao0tfPFnz0s8bLPNNl3H58TgF98f6g03Ra6YAIEmBOIDaXyvd+P9KY43XrVq1TqV8QG0sXUw/owPqJbBFdC/DYCpxAlQim+TirufuBFvrEuXLl2nPq7ndfDBB3dt5YuzdePbNiwECBCoUyAuDB/HEDb2UMSu457XD33uc5/b9d4VWwpdN3Tg15j+bQBMpUyAUnwbLO554kZcdLn7mXhx0HVciqFxHN+BBx5YYgi0ECBAYCgLxPAXxyU3thDGJah6vrfFCSXdj1F2Qknr16j+bQBMpUqAUnzrFHc/cSPeGOOs3Z7X44oTNxrH8b3oRS8qsZvXQoAAgeEsELuH47jlZk4oiaEwTiixdyO/xvVvA2AqRQLUf77GiRuNT8G9nbgRF19tfAqOEzd8w0b/vVUSIDA8BOIEkrg+aWOXsRNKBma96d8GwFSyBGjT+BonbjQ+6caZc92X2KIXF1RtDH1O3Ng0X7cmQKC9BBonlDSGQSeUtG796t8GwFSaBGjDfHHiRmzZa7x59TxxI66VFSduNAa+OO7Fro1UJBUTINDGAnFCSRwP3dhzEieU9DxUJk4oabynOqGk7zDo3wbA1FuFAK3L1zhxo/Hm1NeJG43j+Jy4kYqfYgIEOlwg3nMbJ5TE+25vJ5TEMYONS874hpL/Hxj92wCYevsQoFLiK5G+853vVD+xe6Lnp9G44Gn3T6NO3EhFTjEBAgT6FGicUNLY69Lze87jSgmxVfCwww6rfuIrMDt10b8NgKnsd2KA4svS4wzdxtDX8w0mTtTo/o0bTtxIRUwxAQIE+i0QJ5A0jrmOP3ueUBIf0BvDYByOs8UWW/T7sYZbYSf2757raLO1cZSppV8CnRKglStXliuvvLIa+q666qqyZs2aLq84Zi/eOOJN5FWvelV5/vOf7xs3+pUmRQQIEBg4gWj1t99+e/U+/t3vfrf6LuM4prCxjB49urz0pS+t3stf+cpXlqc//ekD92SGwD13Sv/eELUBMBHEdg1QXJQ0jiVpbOW78cYb11EaM2ZM9QYRbxTxhmG3biJESgkQIFCDQOwujg/0MQzGzwMPPLDOs4gL7ceH+nifj/+Oi++309Ku/XtT1pEBcFO0ety2nQLUeDOIoS+29vX2ZtDYyteObwaJGCglQIDAsBbo/qE/hsHYANB92XHHHbs+9L/sZS9riw/97dS/+xs+A2B/5UopwzlA3XcHxNAXx/X13B0Qv+gx9L3iFa9o+90BiRgoJUCAQFsJNA77iWEwthJGr2sscdhPnE3cOHZwuB72M5z7d6vCZgBMSA63AD388MPVmbqNY0DuueeedV59fNVa45c6fsE76YDgRAyUEiBAoG0F/vGPf1THC8YwGL3jjjvuWOe17rLLLl27iuNC/qNGjRoWFsOtfw8EqgEwoTocAhSXaWn84sbwF19E3ljikgDxC9vYtdvJlwRIxEApAQIEOkagcemv6CvxlXXdL/3VvadEX5kwYcKQdRkO/Xug8QyACeGhGKD4tNb9Mi09P62NHz++ayvfcPq0llhNSgkQIEBgAARir1IMgY2NDPfee++w2as0FPv3AKyiDd6lATAhPlQCFMdrfO973+u6TEs7Hq+RWE1KCRAgQGCABeK48ttuu61rGNzQceVxFYmddtppgJ/Rhu9+qPTvOhHaagC84IILytlnn11iIJo4cWI577zzypQpU/r0veSSS8qpp55a7r777rLbbruVs846qzqWodmlrgA1zthqfOra0BlbcZmWbbfdttmX5HYECBAgQCAtEN8F37jMzIauLBG7iuPr6gb7MjN19e80bAvvoG0GwIsvvrgcddRRZf78+WX//fcv5557bokBb+nSpSWuW9dzie9PPOSQQ8qZZ55ZDj/88HLRRRdVA2B8f+2ee+7ZFPFgBqj7NZvil+kPf/jDOs8xfoEaJ3C4TEtTq8+NCBAgQGAQBGKjxQ033NC1dXBD15aNq09svfXWA/6sBrN/D/iL6ecDtM0AGEPf5MmTy/nnn19RRODGjRtXZs2aVU4++eT1eN70pjeVv/71r+Xb3/52178dcMABZdKkSdUQ2cwykAFqXKalsZWv51Xbn/a0p5XGZVo64artzawPtyFAgACBoS9w//33r3PYUvdvl3rSk57UdZmZgfx2qYHs30N/DfzfZ9gWA2B8P22cen7ppZeW6dOnd9nPmDGjxJazK664Yr31ESdDvO997ysnnHBC17/NmTOnXH755eVXv/pVr+svznbqfsbT6tWrS9zP8uXLS3yNTquWGEBjkI377b7EbuqXv/zl1U8Mq1tuuWWrHtL9ECBAgACBQReI/v0///M/5fvf/371c+edd67zHGJDzsyZM8uxxx7b0ucWA2Dcd5y40n2L48iRI0v8dMLSFgPgihUrytixY0vs1p06dWrXeps9e3a59tpry+LFi9dblzE8ffnLXy5HHnlk179deOGFZe7cueX3v/99r+v+tNNOq/7dQoAAAQIECLSfQGwIil7fCYsBcBMGwJ5bAGM384MPPli23377stlmm7U0L41PJ63eutjSJ9kGd8Z5cFYiZ86DIzA4jyLPw985DrP605/+VLbbbrt1TkCxBXBw1m3LHmWwdgG37Ak3cUeOT2gCqQU34dwCxCbugnMTSC24CecWIDZxF5ybQGrBTTi3AHEDd9EWWwDj9cVJIHHJl7j0SyyxdS6Oz4tjB/o6CSQuYvmtb32ri+fAAw8se+21V9MngQzkqhH8gdT9//fNmfPgCAzOo8gz58ERGJxHkeeBdW6bATAuAxMnfSxYsKAaBOMyMIsWLSpLliypLjgZl4iJ4wTjsi+xxPGCL3rRi8rHP/7x6vIp3/jGN8rHPvaxTboMzECuGsEfSF0D4ODocuY82AKD83jenzkPjsDAPkrbDIDBFGfONi4EHZdzmTdvXrVlMJZp06ZV30u4cOHCLtG4TuBHPvKRrgtBf+ITn9ikC0EP5KqJ4w1jWP3Qhz7UMWckDaRnX/fNeXDUOXMeHIHBeRR55jw4AgP7KG01AA4slXsnQIAAAQIECLSHgAGwPdajV0GAAAECBAgQaFrAANg0lRsSIECAAAECBNpDwADYHuvRqyBAgAABAgQINC1gAGyayg0JECBAgAABAu0hYAAcguvxggsu6DqbeeLEidW1DePSNpbWCMTZ1d/85jerSwQ9+clPLnH9x7POOqs897nPbc0DuJdeBeKSS3FW+3vf+97qMk2W1grcd9995aSTTipXXnlliWucPuc5zylf+tKXyn777dfaB+rge3v88cerrwn72te+VlauXFme+cxnlqOPPrq6mkSrvw2qk5h//OMfVz3vxhtvLPfff3+57LLLyvTp07sI4ls74ivaPve5z5VVq1aVgw46qHzmM58pu+22Wycxtfy1GgBbTpq7w7ieYVyzcP78+dUlbKJRxuVqli5dWsaMGZO7c9WVwCte8YpyxBFHlMmTJ5fHHnusnHLKKeXWW28tt99+e3nKU55CaQAEbrjhhvLGN76xjB49uhx66KEGwBYb//nPfy577713Zfvud7+77LjjjuXOO+8su+66a/VjaY1AXCv2nHPOqb5Hfo899ii/+MUvytve9rZyxhlnlOOPP741D9KB9xIfWn7605+Wfffdt7z2ta9dbwCMD+jxwT3cn/WsZ5VTTz213HLLLdV79lZbbdWBYq15yQbA1ji27F5i6IvBJK5pGEt8o8m4cePKrFmzev1Gk5Y9cAff0QMPPFAN19dee2055JBDOlhiYF76Qw89VPbZZ59y4YUXltNPP73ENTptAWytdXzbUTTQn/zkJ629Y/e2jsDhhx9efbHAF77wha6/f93rXlftSYitgpa8QGxJ7b4FMLb+xZbW97///eUDH/hA9QCrV6+u1kNc1zc+zFv6J2AA7J/bgFT15zuNB+SJdNidLlu2rNqVEJ8o99xzzw579QP/cuMbeuIL1z/1qU9VF2Q3ALbe/PnPf355+ctfXn73u99VH2TiW4+OO+648s53vrP1D9bB9xhbAD/72c+Wq666quy+++7lV7/6VXnZy15WbRV8y1ve0sEyrXvpPQfAu+66q9qK/ctf/rJ672gs8U1e8f+f/vSnW/fgHXZPBsAhtMJXrFhRvXHH19RNnTq165nNnj27elNfvHjxEHq27fFUYgvrv/zLv1THlVx33XXt8aKG0KuIr1iM3WOxCzh21RgAB2blNHaDve997ytveMMbKu841jIOJYkB3NIagXi/iENG4lujNt988xLHBEa+49hWS2sEeg6A0Q/jmL/oj894xjO6HiQOKYnbxmFTlv4JGAD75zYgVQbAAWHd4J3G8VJx/EkMfzvvvPPgP4E2fvpOLSgAAASTSURBVMTly5dXJyD84Ac/KHvttVf1Sg2AA7PCt9xyy8o6mmVjiWPSYhD82c9+NjAP2oH3Gh9oPvjBD1YnLMQxgDfddFM54YQTqi2ABu3WBMIA2BrHZu7FANiM0iDdxi7gQYL+fw8zc+bMcsUVV5Q4Ay0OLLa0VuDyyy8vr3nNa6otJY0ltpjEG/yIESNKfJ9q939r7aN31r3tsssu5aUvfWn5/Oc/3/XC4yzJOOYyzg62tEYgjseO4y3f8573dN1hGMfxf3FVAUtewC7gvGGz92AAbFZqkG4XJ4HEJV/i0i+xxC6H8ePHlxhW4o3HkheIg4rjpJo40Piaa65xKYE8aa/3sGbNmnLPPfes829xxuTznve86nIljrdsHfyb3/zmEltcu58EcuKJJ1aHjXTfKti6R+zMe9p+++2roTr2HDSWODs1Lrfz61//ujNRWvyq+zoJJE4AiRNBYvnLX/5SnbjnJJAcvgEw59fy6jieIXYlLFiwoBoE42zJRYsWVZ8u46wnS14gDo6/6KKLqq1/3a/9t/XWW1dn81kGTsAu4IGxjV29cT3LuXPnVpfb+fnPf16dABInLDg5oXXmcc2/q6++unp/jl3AcWLCu971rnLMMcdU1xK19E8grhQQJ+PFEpczil3qcUmjOHksNoCEbVxHtPtlYG6++WaXgekfd1eVATAJOBDlcQmYOMYkLjQaZznNmzevuiagpTUCfV2wNT7Fxxu8ZeAEDIADZ/vtb3+7Ohkhrv8XhzTECSHOAm6td2zVjmvQxd6DP/zhD9XlSY488sjy0Y9+tMRxmJb+CcSemBj4ei6xMSS28jUuBB0faOKEvYMPPri6rFSciW3pv4ABsP92KgkQIECAAAECw1LAADgsV5snTYAAAQIECBDov4ABsP92KgkQIECAAAECw1LAADgsV5snTYAAAQIECBDov4ABsP92KgkQIECAAAECw1LAADgsV5snTYAAAQIECBDov4ABsP92KgkQIECAAAECw1LAADgsV5snTYAAAQIECBDov4ABsP92KgkQIECAAAECw1LAADgsV5snTYAAAQIECBDov4ABsP92KgkQIECAAAECw1LAADgsV5snTYAAAQIECBDov4ABsP92KgkQIECAAAECw1LAADgsV5snTYAAAQIECBDov4ABsP92KgkQIECAAAECw1LAADgsV5snTYAAAQIECBDov4ABsP92KgkQIECAAAECw1LAADgsV5snTYAAAQIECBDov4ABsP92KgkQIECAAAECw1LAADgsV5snTYAAAQIECBDov4ABsP92KgkQIECAAAECw1LAADgsV5snTYAAAQIECBDov4ABsP92KgkQIECAAAECw1LAADgsV5snTYAAAQIECBDov4ABsP92KgkQIECAAAECw1LAADgsV5snTYAAAQIECBDov4ABsP92KgkQIECAAAECw1LAADgsV5snTYAAAQIECBDov4ABsP92KgkQIECAAAECw1LAADgsV5snTYAAAQIECBDov4ABsP92KgkQIECAAAECw1LAADgsV5snTYAAAQIECBDov8D/AfO4PebSCMLhAAAAAElFTkSuQmCC\" width=\"640\">"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/plain": [
       "<matplotlib.animation.FuncAnimation at 0x7fe7e5ca8520>"
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "animate_line(u, t, sampling=1)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "id": "4eb57be9",
   "metadata": {
    "scrolled": false
   },
   "outputs": [
    {
     "data": {
      "application/javascript": [
       "/* Put everything inside the global mpl namespace */\n",
       "/* global mpl */\n",
       "window.mpl = {};\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(\n",
       "            '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",
       "\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 = document.createElement('div');\n",
       "    this.root.setAttribute('style', 'display: inline-block');\n",
       "    this._root_extra_style(this.root);\n",
       "\n",
       "    parent_element.appendChild(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 (fig.ratio !== 1) {\n",
       "            fig.send_message('set_dpi_ratio', { dpi_ratio: fig.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 = document.createElement('div');\n",
       "    titlebar.classList =\n",
       "        'ui-dialog-titlebar ui-widget-header ui-corner-all ui-helper-clearfix';\n",
       "    var titletext = document.createElement('div');\n",
       "    titletext.classList = 'ui-dialog-title';\n",
       "    titletext.setAttribute(\n",
       "        'style',\n",
       "        'width: 100%; text-align: center; padding: 3px;'\n",
       "    );\n",
       "    titlebar.appendChild(titletext);\n",
       "    this.root.appendChild(titlebar);\n",
       "    this.header = titletext;\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._canvas_extra_style = function (_canvas_div) {};\n",
       "\n",
       "mpl.figure.prototype._root_extra_style = function (_canvas_div) {};\n",
       "\n",
       "mpl.figure.prototype._init_canvas = function () {\n",
       "    var fig = this;\n",
       "\n",
       "    var canvas_div = (this.canvas_div = document.createElement('div'));\n",
       "    canvas_div.setAttribute(\n",
       "        'style',\n",
       "        'border: 1px solid #ddd;' +\n",
       "            'box-sizing: content-box;' +\n",
       "            'clear: both;' +\n",
       "            'min-height: 1px;' +\n",
       "            'min-width: 1px;' +\n",
       "            'outline: 0;' +\n",
       "            'overflow: hidden;' +\n",
       "            'position: relative;' +\n",
       "            'resize: both;'\n",
       "    );\n",
       "\n",
       "    function on_keyboard_event_closure(name) {\n",
       "        return function (event) {\n",
       "            return fig.key_event(event, name);\n",
       "        };\n",
       "    }\n",
       "\n",
       "    canvas_div.addEventListener(\n",
       "        'keydown',\n",
       "        on_keyboard_event_closure('key_press')\n",
       "    );\n",
       "    canvas_div.addEventListener(\n",
       "        'keyup',\n",
       "        on_keyboard_event_closure('key_release')\n",
       "    );\n",
       "\n",
       "    this._canvas_extra_style(canvas_div);\n",
       "    this.root.appendChild(canvas_div);\n",
       "\n",
       "    var canvas = (this.canvas = document.createElement('canvas'));\n",
       "    canvas.classList.add('mpl-canvas');\n",
       "    canvas.setAttribute('style', 'box-sizing: content-box;');\n",
       "\n",
       "    this.context = canvas.getContext('2d');\n",
       "\n",
       "    var backingStore =\n",
       "        this.context.backingStorePixelRatio ||\n",
       "        this.context.webkitBackingStorePixelRatio ||\n",
       "        this.context.mozBackingStorePixelRatio ||\n",
       "        this.context.msBackingStorePixelRatio ||\n",
       "        this.context.oBackingStorePixelRatio ||\n",
       "        this.context.backingStorePixelRatio ||\n",
       "        1;\n",
       "\n",
       "    this.ratio = (window.devicePixelRatio || 1) / backingStore;\n",
       "\n",
       "    var rubberband_canvas = (this.rubberband_canvas = document.createElement(\n",
       "        'canvas'\n",
       "    ));\n",
       "    rubberband_canvas.setAttribute(\n",
       "        'style',\n",
       "        'box-sizing: content-box; position: absolute; left: 0; top: 0; z-index: 1;'\n",
       "    );\n",
       "\n",
       "    // Apply a ponyfill if ResizeObserver is not implemented by browser.\n",
       "    if (this.ResizeObserver === undefined) {\n",
       "        if (window.ResizeObserver !== undefined) {\n",
       "            this.ResizeObserver = window.ResizeObserver;\n",
       "        } else {\n",
       "            var obs = _JSXTOOLS_RESIZE_OBSERVER({});\n",
       "            this.ResizeObserver = obs.ResizeObserver;\n",
       "        }\n",
       "    }\n",
       "\n",
       "    this.resizeObserverInstance = new this.ResizeObserver(function (entries) {\n",
       "        var nentries = entries.length;\n",
       "        for (var i = 0; i < nentries; i++) {\n",
       "            var entry = entries[i];\n",
       "            var width, height;\n",
       "            if (entry.contentBoxSize) {\n",
       "                if (entry.contentBoxSize instanceof Array) {\n",
       "                    // Chrome 84 implements new version of spec.\n",
       "                    width = entry.contentBoxSize[0].inlineSize;\n",
       "                    height = entry.contentBoxSize[0].blockSize;\n",
       "                } else {\n",
       "                    // Firefox implements old version of spec.\n",
       "                    width = entry.contentBoxSize.inlineSize;\n",
       "                    height = entry.contentBoxSize.blockSize;\n",
       "                }\n",
       "            } else {\n",
       "                // Chrome <84 implements even older version of spec.\n",
       "                width = entry.contentRect.width;\n",
       "                height = entry.contentRect.height;\n",
       "            }\n",
       "\n",
       "            // Keep the size of the canvas and rubber band canvas in sync with\n",
       "            // the canvas container.\n",
       "            if (entry.devicePixelContentBoxSize) {\n",
       "                // Chrome 84 implements new version of spec.\n",
       "                canvas.setAttribute(\n",
       "                    'width',\n",
       "                    entry.devicePixelContentBoxSize[0].inlineSize\n",
       "                );\n",
       "                canvas.setAttribute(\n",
       "                    'height',\n",
       "                    entry.devicePixelContentBoxSize[0].blockSize\n",
       "                );\n",
       "            } else {\n",
       "                canvas.setAttribute('width', width * fig.ratio);\n",
       "                canvas.setAttribute('height', height * fig.ratio);\n",
       "            }\n",
       "            canvas.setAttribute(\n",
       "                'style',\n",
       "                'width: ' + width + 'px; height: ' + height + 'px;'\n",
       "            );\n",
       "\n",
       "            rubberband_canvas.setAttribute('width', width);\n",
       "            rubberband_canvas.setAttribute('height', height);\n",
       "\n",
       "            // And update the size in Python. We ignore the initial 0/0 size\n",
       "            // that occurs as the element is placed into the DOM, which should\n",
       "            // otherwise not happen due to the minimum size styling.\n",
       "            if (fig.ws.readyState == 1 && width != 0 && height != 0) {\n",
       "                fig.request_resize(width, height);\n",
       "            }\n",
       "        }\n",
       "    });\n",
       "    this.resizeObserverInstance.observe(canvas_div);\n",
       "\n",
       "    function on_mouse_event_closure(name) {\n",
       "        return function (event) {\n",
       "            return fig.mouse_event(event, name);\n",
       "        };\n",
       "    }\n",
       "\n",
       "    rubberband_canvas.addEventListener(\n",
       "        'mousedown',\n",
       "        on_mouse_event_closure('button_press')\n",
       "    );\n",
       "    rubberband_canvas.addEventListener(\n",
       "        'mouseup',\n",
       "        on_mouse_event_closure('button_release')\n",
       "    );\n",
       "    rubberband_canvas.addEventListener(\n",
       "        'dblclick',\n",
       "        on_mouse_event_closure('dblclick')\n",
       "    );\n",
       "    // Throttle sequential mouse events to 1 every 20ms.\n",
       "    rubberband_canvas.addEventListener(\n",
       "        'mousemove',\n",
       "        on_mouse_event_closure('motion_notify')\n",
       "    );\n",
       "\n",
       "    rubberband_canvas.addEventListener(\n",
       "        'mouseenter',\n",
       "        on_mouse_event_closure('figure_enter')\n",
       "    );\n",
       "    rubberband_canvas.addEventListener(\n",
       "        'mouseleave',\n",
       "        on_mouse_event_closure('figure_leave')\n",
       "    );\n",
       "\n",
       "    canvas_div.addEventListener('wheel', function (event) {\n",
       "        if (event.deltaY < 0) {\n",
       "            event.step = 1;\n",
       "        } else {\n",
       "            event.step = -1;\n",
       "        }\n",
       "        on_mouse_event_closure('scroll')(event);\n",
       "    });\n",
       "\n",
       "    canvas_div.appendChild(canvas);\n",
       "    canvas_div.appendChild(rubberband_canvas);\n",
       "\n",
       "    this.rubberband_context = rubberband_canvas.getContext('2d');\n",
       "    this.rubberband_context.strokeStyle = '#000000';\n",
       "\n",
       "    this._resize_canvas = function (width, height, forward) {\n",
       "        if (forward) {\n",
       "            canvas_div.style.width = width + 'px';\n",
       "            canvas_div.style.height = height + 'px';\n",
       "        }\n",
       "    };\n",
       "\n",
       "    // Disable right mouse context menu.\n",
       "    this.rubberband_canvas.addEventListener('contextmenu', function (_e) {\n",
       "        event.preventDefault();\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 toolbar = document.createElement('div');\n",
       "    toolbar.classList = 'mpl-toolbar';\n",
       "    this.root.appendChild(toolbar);\n",
       "\n",
       "    function on_click_closure(name) {\n",
       "        return function (_event) {\n",
       "            return fig.toolbar_button_onclick(name);\n",
       "        };\n",
       "    }\n",
       "\n",
       "    function on_mouseover_closure(tooltip) {\n",
       "        return function (event) {\n",
       "            if (!event.currentTarget.disabled) {\n",
       "                return fig.toolbar_button_onmouseover(tooltip);\n",
       "            }\n",
       "        };\n",
       "    }\n",
       "\n",
       "    fig.buttons = {};\n",
       "    var buttonGroup = document.createElement('div');\n",
       "    buttonGroup.classList = 'mpl-button-group';\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",
       "            /* Instead of a spacer, we start a new button group. */\n",
       "            if (buttonGroup.hasChildNodes()) {\n",
       "                toolbar.appendChild(buttonGroup);\n",
       "            }\n",
       "            buttonGroup = document.createElement('div');\n",
       "            buttonGroup.classList = 'mpl-button-group';\n",
       "            continue;\n",
       "        }\n",
       "\n",
       "        var button = (fig.buttons[name] = document.createElement('button'));\n",
       "        button.classList = 'mpl-widget';\n",
       "        button.setAttribute('role', 'button');\n",
       "        button.setAttribute('aria-disabled', 'false');\n",
       "        button.addEventListener('click', on_click_closure(method_name));\n",
       "        button.addEventListener('mouseover', on_mouseover_closure(tooltip));\n",
       "\n",
       "        var icon_img = document.createElement('img');\n",
       "        icon_img.src = '_images/' + image + '.png';\n",
       "        icon_img.srcset = '_images/' + image + '_large.png 2x';\n",
       "        icon_img.alt = tooltip;\n",
       "        button.appendChild(icon_img);\n",
       "\n",
       "        buttonGroup.appendChild(button);\n",
       "    }\n",
       "\n",
       "    if (buttonGroup.hasChildNodes()) {\n",
       "        toolbar.appendChild(buttonGroup);\n",
       "    }\n",
       "\n",
       "    var fmt_picker = document.createElement('select');\n",
       "    fmt_picker.classList = 'mpl-widget';\n",
       "    toolbar.appendChild(fmt_picker);\n",
       "    this.format_dropdown = fmt_picker;\n",
       "\n",
       "    for (var ind in mpl.extensions) {\n",
       "        var fmt = mpl.extensions[ind];\n",
       "        var option = document.createElement('option');\n",
       "        option.selected = fmt === mpl.default_extension;\n",
       "        option.innerHTML = fmt;\n",
       "        fmt_picker.appendChild(option);\n",
       "    }\n",
       "\n",
       "    var status_bar = document.createElement('span');\n",
       "    status_bar.classList = 'mpl-message';\n",
       "    toolbar.appendChild(status_bar);\n",
       "    this.message = status_bar;\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",
       "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",
       "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], msg['forward']);\n",
       "        fig.send_message('refresh', {});\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_rubberband = function (fig, msg) {\n",
       "    var x0 = msg['x0'] / fig.ratio;\n",
       "    var y0 = (fig.canvas.height - msg['y0']) / fig.ratio;\n",
       "    var x1 = msg['x1'] / fig.ratio;\n",
       "    var y1 = (fig.canvas.height - msg['y1']) / fig.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,\n",
       "        0,\n",
       "        fig.canvas.width / fig.ratio,\n",
       "        fig.canvas.height / fig.ratio\n",
       "    );\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",
       "        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.handle_history_buttons = function (fig, msg) {\n",
       "    for (var key in msg) {\n",
       "        if (!(key in fig.buttons)) {\n",
       "            continue;\n",
       "        }\n",
       "        fig.buttons[key].disabled = !msg[key];\n",
       "        fig.buttons[key].setAttribute('aria-disabled', !msg[key]);\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_navigate_mode = function (fig, msg) {\n",
       "    if (msg['mode'] === 'PAN') {\n",
       "        fig.buttons['Pan'].classList.add('active');\n",
       "        fig.buttons['Zoom'].classList.remove('active');\n",
       "    } else if (msg['mode'] === 'ZOOM') {\n",
       "        fig.buttons['Pan'].classList.remove('active');\n",
       "        fig.buttons['Zoom'].classList.add('active');\n",
       "    } else {\n",
       "        fig.buttons['Pan'].classList.remove('active');\n",
       "        fig.buttons['Zoom'].classList.remove('active');\n",
       "    }\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",
       "            var img = evt.data;\n",
       "            if (img.type !== 'image/png') {\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",
       "                img.type = 'image/png';\n",
       "            }\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",
       "\n",
       "            fig.imageObj.src = (window.URL || window.webkitURL).createObjectURL(\n",
       "                img\n",
       "            );\n",
       "            fig.updated_canvas_event();\n",
       "            fig.waiting = false;\n",
       "            return;\n",
       "        } else if (\n",
       "            typeof evt.data === 'string' &&\n",
       "            evt.data.slice(0, 21) === 'data:image/png;base64'\n",
       "        ) {\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(\n",
       "                \"No handler for the '\" + msg_type + \"' message type: \",\n",
       "                msg\n",
       "            );\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(\n",
       "                    \"Exception inside the 'handler_\" + msg_type + \"' callback:\",\n",
       "                    e,\n",
       "                    e.stack,\n",
       "                    msg\n",
       "                );\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",
       "    }\n",
       "    if (e.target) {\n",
       "        targ = e.target;\n",
       "    } else if (e.srcElement) {\n",
       "        targ = e.srcElement;\n",
       "    }\n",
       "    if (targ.nodeType === 3) {\n",
       "        // defeat Safari bug\n",
       "        targ = targ.parentNode;\n",
       "    }\n",
       "\n",
       "    // pageX,Y are the mouse positions relative to the document\n",
       "    var boundingRect = targ.getBoundingClientRect();\n",
       "    var x = e.pageX - (boundingRect.left + document.body.scrollLeft);\n",
       "    var y = e.pageY - (boundingRect.top + document.body.scrollTop);\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",
       "        }\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",
       "        this.canvas.focus();\n",
       "        this.canvas_div.focus();\n",
       "    }\n",
       "\n",
       "    var x = canvas_pos.x * this.ratio;\n",
       "    var y = canvas_pos.y * this.ratio;\n",
       "\n",
       "    this.send_message(name, {\n",
       "        x: x,\n",
       "        y: y,\n",
       "        button: event.button,\n",
       "        step: event.step,\n",
       "        guiEvent: simpleKeys(event),\n",
       "    });\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",
       "    // Prevent repeat events\n",
       "    if (name === 'key_press') {\n",
       "        if (event.key === this._key) {\n",
       "            return;\n",
       "        } else {\n",
       "            this._key = event.key;\n",
       "        }\n",
       "    }\n",
       "    if (name === 'key_release') {\n",
       "        this._key = null;\n",
       "    }\n",
       "\n",
       "    var value = '';\n",
       "    if (event.ctrlKey && event.key !== 'Control') {\n",
       "        value += 'ctrl+';\n",
       "    }\n",
       "    else if (event.altKey && event.key !== 'Alt') {\n",
       "        value += 'alt+';\n",
       "    }\n",
       "    else if (event.shiftKey && event.key !== 'Shift') {\n",
       "        value += 'shift+';\n",
       "    }\n",
       "\n",
       "    value += 'k' + event.key;\n",
       "\n",
       "    this._key_event_extra(event, name);\n",
       "\n",
       "    this.send_message(name, { key: value, 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",
       "\n",
       "///////////////// REMAINING CONTENT GENERATED BY embed_js.py /////////////////\n",
       "// prettier-ignore\n",
       "var _JSXTOOLS_RESIZE_OBSERVER=function(A){var t,i=new WeakMap,n=new WeakMap,a=new WeakMap,r=new WeakMap,o=new Set;function s(e){if(!(this instanceof s))throw new TypeError(\"Constructor requires 'new' operator\");i.set(this,e)}function h(){throw new TypeError(\"Function is not a constructor\")}function c(e,t,i,n){e=0 in arguments?Number(arguments[0]):0,t=1 in arguments?Number(arguments[1]):0,i=2 in arguments?Number(arguments[2]):0,n=3 in arguments?Number(arguments[3]):0,this.right=(this.x=this.left=e)+(this.width=i),this.bottom=(this.y=this.top=t)+(this.height=n),Object.freeze(this)}function d(){t=requestAnimationFrame(d);var s=new WeakMap,p=new Set;o.forEach((function(t){r.get(t).forEach((function(i){var r=t instanceof window.SVGElement,o=a.get(t),d=r?0:parseFloat(o.paddingTop),f=r?0:parseFloat(o.paddingRight),l=r?0:parseFloat(o.paddingBottom),u=r?0:parseFloat(o.paddingLeft),g=r?0:parseFloat(o.borderTopWidth),m=r?0:parseFloat(o.borderRightWidth),w=r?0:parseFloat(o.borderBottomWidth),b=u+f,F=d+l,v=(r?0:parseFloat(o.borderLeftWidth))+m,W=g+w,y=r?0:t.offsetHeight-W-t.clientHeight,E=r?0:t.offsetWidth-v-t.clientWidth,R=b+v,z=F+W,M=r?t.width:parseFloat(o.width)-R-E,O=r?t.height:parseFloat(o.height)-z-y;if(n.has(t)){var k=n.get(t);if(k[0]===M&&k[1]===O)return}n.set(t,[M,O]);var S=Object.create(h.prototype);S.target=t,S.contentRect=new c(u,d,M,O),s.has(i)||(s.set(i,[]),p.add(i)),s.get(i).push(S)}))})),p.forEach((function(e){i.get(e).call(e,s.get(e),e)}))}return s.prototype.observe=function(i){if(i instanceof window.Element){r.has(i)||(r.set(i,new Set),o.add(i),a.set(i,window.getComputedStyle(i)));var n=r.get(i);n.has(this)||n.add(this),cancelAnimationFrame(t),t=requestAnimationFrame(d)}},s.prototype.unobserve=function(i){if(i instanceof window.Element&&r.has(i)){var n=r.get(i);n.has(this)&&(n.delete(this),n.size||(r.delete(i),o.delete(i))),n.size||r.delete(i),o.size||cancelAnimationFrame(t)}},A.DOMRectReadOnly=c,A.ResizeObserver=s,A.ResizeObserverEntry=h,A}; // eslint-disable-line\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\", \"Left button pans, Right button zooms\\nx/y fixes axis, CTRL fixes aspect\", \"fa fa-arrows icon-move\", \"pan\"], [\"Zoom\", \"Zoom to rectangle\\nx/y fixes axis, CTRL fixes aspect\", \"fa fa-square-o icon-check-empty\", \"zoom\"], [\"\", \"\", \"\", \"\"], [\"Download\", \"Download plot\", \"fa fa-floppy-o icon-save\", \"download\"]];\n",
       "\n",
       "mpl.extensions = [\"eps\", \"jpeg\", \"pgf\", \"pdf\", \"png\", \"ps\", \"raw\", \"svg\", \"tif\"];\n",
       "\n",
       "mpl.default_extension = \"png\";/* global mpl */\n",
       "\n",
       "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.binaryType = comm.kernel.ws.binaryType;\n",
       "    ws.readyState = comm.kernel.ws.readyState;\n",
       "    function updateReadyState(_event) {\n",
       "        if (comm.kernel.ws) {\n",
       "            ws.readyState = comm.kernel.ws.readyState;\n",
       "        } else {\n",
       "            ws.readyState = 3; // Closed state.\n",
       "        }\n",
       "    }\n",
       "    comm.kernel.ws.addEventListener('open', updateReadyState);\n",
       "    comm.kernel.ws.addEventListener('close', updateReadyState);\n",
       "    comm.kernel.ws.addEventListener('error', updateReadyState);\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",
       "        var data = msg['content']['data'];\n",
       "        if (data['blob'] !== undefined) {\n",
       "            data = {\n",
       "                data: new Blob(msg['buffers'], { type: data['blob'] }),\n",
       "            };\n",
       "        }\n",
       "        // Pass the mpl event to the overridden (by mpl) onmessage function.\n",
       "        ws.onmessage(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 = document.getElementById(id);\n",
       "    var ws_proxy = comm_websocket_adapter(comm);\n",
       "\n",
       "    function ondownload(figure, _format) {\n",
       "        window.open(figure.canvas.toDataURL());\n",
       "    }\n",
       "\n",
       "    var fig = new mpl.figure(id, ws_proxy, ondownload, element);\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;\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",
       "    fig.cell_info[0].output_area.element.on(\n",
       "        'cleared',\n",
       "        { fig: fig },\n",
       "        fig._remove_fig_handler\n",
       "    );\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_close = function (fig, msg) {\n",
       "    var width = fig.canvas.width / fig.ratio;\n",
       "    fig.cell_info[0].output_area.element.off(\n",
       "        'cleared',\n",
       "        fig._remove_fig_handler\n",
       "    );\n",
       "    fig.resizeObserverInstance.unobserve(fig.canvas_div);\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.innerHTML =\n",
       "        '<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 / this.ratio;\n",
       "    var dataURL = this.canvas.toDataURL();\n",
       "    this.cell_info[1]['text/html'] =\n",
       "        '<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 () {\n",
       "        fig.push_to_output();\n",
       "    }, 1000);\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._init_toolbar = function () {\n",
       "    var fig = this;\n",
       "\n",
       "    var toolbar = document.createElement('div');\n",
       "    toolbar.classList = 'btn-toolbar';\n",
       "    this.root.appendChild(toolbar);\n",
       "\n",
       "    function on_click_closure(name) {\n",
       "        return function (_event) {\n",
       "            return fig.toolbar_button_onclick(name);\n",
       "        };\n",
       "    }\n",
       "\n",
       "    function on_mouseover_closure(tooltip) {\n",
       "        return function (event) {\n",
       "            if (!event.currentTarget.disabled) {\n",
       "                return fig.toolbar_button_onmouseover(tooltip);\n",
       "            }\n",
       "        };\n",
       "    }\n",
       "\n",
       "    fig.buttons = {};\n",
       "    var buttonGroup = document.createElement('div');\n",
       "    buttonGroup.classList = 'btn-group';\n",
       "    var button;\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",
       "            /* Instead of a spacer, we start a new button group. */\n",
       "            if (buttonGroup.hasChildNodes()) {\n",
       "                toolbar.appendChild(buttonGroup);\n",
       "            }\n",
       "            buttonGroup = document.createElement('div');\n",
       "            buttonGroup.classList = 'btn-group';\n",
       "            continue;\n",
       "        }\n",
       "\n",
       "        button = fig.buttons[name] = document.createElement('button');\n",
       "        button.classList = 'btn btn-default';\n",
       "        button.href = '#';\n",
       "        button.title = name;\n",
       "        button.innerHTML = '<i class=\"fa ' + image + ' fa-lg\"></i>';\n",
       "        button.addEventListener('click', on_click_closure(method_name));\n",
       "        button.addEventListener('mouseover', on_mouseover_closure(tooltip));\n",
       "        buttonGroup.appendChild(button);\n",
       "    }\n",
       "\n",
       "    if (buttonGroup.hasChildNodes()) {\n",
       "        toolbar.appendChild(buttonGroup);\n",
       "    }\n",
       "\n",
       "    // Add the status bar.\n",
       "    var status_bar = document.createElement('span');\n",
       "    status_bar.classList = 'mpl-message pull-right';\n",
       "    toolbar.appendChild(status_bar);\n",
       "    this.message = status_bar;\n",
       "\n",
       "    // Add the close button to the window.\n",
       "    var buttongrp = document.createElement('div');\n",
       "    buttongrp.classList = 'btn-group inline pull-right';\n",
       "    button = document.createElement('button');\n",
       "    button.classList = 'btn btn-mini btn-primary';\n",
       "    button.href = '#';\n",
       "    button.title = 'Stop Interaction';\n",
       "    button.innerHTML = '<i class=\"fa fa-power-off icon-remove icon-large\"></i>';\n",
       "    button.addEventListener('click', function (_evt) {\n",
       "        fig.handle_close(fig, {});\n",
       "    });\n",
       "    button.addEventListener(\n",
       "        'mouseover',\n",
       "        on_mouseover_closure('Stop Interaction')\n",
       "    );\n",
       "    buttongrp.appendChild(button);\n",
       "    var titlebar = this.root.querySelector('.ui-dialog-titlebar');\n",
       "    titlebar.insertBefore(buttongrp, titlebar.firstChild);\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._remove_fig_handler = function (event) {\n",
       "    var fig = event.data.fig;\n",
       "    if (event.target !== this) {\n",
       "        // Ignore bubbled events from children.\n",
       "        return;\n",
       "    }\n",
       "    fig.close_ws(fig, {});\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._root_extra_style = function (el) {\n",
       "    el.style.boxSizing = 'content-box'; // override notebook setting of border-box.\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._canvas_extra_style = function (el) {\n",
       "    // this is important to make the div 'focusable\n",
       "    el.setAttribute('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",
       "    } else {\n",
       "        // location in version 2\n",
       "        IPython.keyboard_manager.register_events(el);\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",
       "\n",
       "    // Check for shift+enter\n",
       "    if (event.shiftKey && event.which === 13) {\n",
       "        this.canvas_div.blur();\n",
       "        // select the cell after this one\n",
       "        var index = IPython.notebook.find_cell_index(this.cell_info[0]);\n",
       "        IPython.notebook.select(index + 1);\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_save = function (fig, _msg) {\n",
       "    fig.ondownload(fig, null);\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(\n",
       "        'matplotlib',\n",
       "        mpl.mpl_figure_comm\n",
       "    );\n",
       "}\n"
      ],
      "text/plain": [
       "<IPython.core.display.Javascript object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/html": [
       "<img src=\"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAyAAAAEsCAYAAAA7Ldc6AAAAAXNSR0IArs4c6QAAIABJREFUeF7tnQewZGWhdb9hyFHJSk5CCRIkSJIHkiQpiCAYAEHJOYNEiyyCoiBIUBCJFmgBKigGUChSgUAJKCA5CEqUJDh/fed/M29mgBnudqZndX+rq6hXT+7u3mft5t67bp/TPWzEiBEjijcJSEACEpCABCQgAQlIQAI9IDBMAekBZR9CAhKQgAQkIAEJSEACEugIKCA+ESQgAQlIQAISkIAEJCCBnhFQQHqG2geSgAQkIAEJSEACEpCABBQQnwMSkIAEJCABCUhAAhKQQM8IKCA9Q+0DSUACEpCABCQgAQlIQAIKiM8BCUhAAhKQgAQkIAEJSKBnBBSQnqH2gSQgAQlIQAISkIAEJCABBcTngAQkIAEJSEACEpCABCTQMwIKSM9Q+0ASkIAEJCABCUhAAhKQgALic0ACEpCABCQgAQlIQAIS6BkBBaRnqH0gCUhAAhKQgAQkIAEJSEAB8TkgAQlIQAISkIAEJCABCfSMgALSM9Q+kAQkIAEJSEACEpCABCSggPgckIAEJCABCUhAAhKQgAR6RkAB6RlqH0gCEpDAxCVwww03lGuuuabsueee5X3ve9/EfbB3uPfnn3++7L///uXyyy8vr7zySllhhRXKN7/5zfLRj350nF3+85//lPPOO69cdtll5fbbby///Oc/ywILLFC22GKLsu+++5app576bfmzzz67nHjiieVvf/tbmWeeecruu+9edtttt7d93eOPP1722muvjkt9nDXWWKOcfPLJZcEFF+w5Hx9QAhKQgAT+PwEFxGeCBCQggQEhUH8h32+//bpfyueff/6eHlX95f7jH/94+dOf/tR1mHXWWctpp51WHn300XLbbbeVRRZZ5F37vPzyy2WGGWYoK664Ytlwww3L7LPPXm688cZy7rnnltVWW6385je/KcOGDRuVP+OMM8qOO+5YNt1007LuuuuW66+/vvzoRz8qxx13XDnggANGfV293yo/L7zwQtlnn33KFFNM0cnHiBEjyh133FFmmWWWnjLywSQgAQlIQAHxOSABCUhgoAhMSgG55JJLyuc+97ly6aWXls9+9rMd12eeeaZ86EMfKuutt1654IIL3pX1G2+8UW699day8sorj/E1X//618vhhx9efvWrX5W11lqr+3evvvpq94pHlZUrr7xy1Nd/8YtfLD/96U874Xn/+9/f/e8nnHBCJyQ333xzWX755bv/7d577y1LLLFE90rNMcccM1D7ezASkIAE+oWAr4D0y1L2lIAEJDAOAkcccUQ58sgj3/YVvXo1ZPPNNy/XXXddeeKJJ8pkk002qscOO+xQzj///O60qqmmmmpIG951111lySWXLKeccsqo06t+/vOflw022KBcddVVZf311x91f/UVkyow9ZWQKiP1Vk8Bq7cqIKPf6qsmDzzwQLn//vuH1McvloAEJCCBCUNAAZkwHL0XCUhAApOUwJ133tmdgnThhRd2pxnVU6DqbZNNNinTTTfdO3ar12nUf8Z3Gz58+KhXFd7ta+spVvWfKgij3+q1Gl/5yldK7feRj3xkfA81xr+vr3yss8463asnW265Zffvjj766HLIIYeUp59+ujtVa+StvooyzTTTdNe/1OtO6ilh0047bdl22227U8FGvx166KHlqKOOKi+++GJ36pc3CUhAAhLoLQEFpLe8fTQJSEACE43AUE/BerdXTcYuON9885WHHnponL2nn3767hSsKhyj30a+YvHLX/6yu15jKLe11167e/Xi4YcfHnVR/a677lpOP/308uabb77trqqQrLnmmp2EPfvss2W22WYr9TSuKhyj36qQ7LLLLt3pWIsuuuhQKvm1EpCABCQwAQgoIBMAonchAQlIgEBgqALy4IMPlvrP+G71lYVVVlllnF9WXyWpp1uN/WpDvYC8SkF9Z6yNN954fA816t/X6zO+9rWvdfe30047jfrft9tuu04w3umVm3nnnbe76HzktSD1/z/++OO76z1Gv51zzjml3k99x62ll176PXfyCyUgAQlIYMIQUEAmDEfvRQISkMAkJzBUAZmQhSfkKyAXX3xxd8pVPX3qrLPOGqOmr4BMyNW8LwlIQAKThoACMmm4+6gSkIAEJjiBoQpIfZva+s/4bvXVjXo607huE+oakHrdR30r3nrtR33VZPLJJx/jYb0GZHxr+e8lIAEJ8AkoIPyNbCgBCUjgPRGoF1/XD+57r+98NSGvAdlss826z+MY+12wtt9++/LjH//4Pb0L1k033dSdrrXUUkuVX//6191F5WPf6rtfVUEZ+12w6ocw1tPE6gcafulLX+pi9a136+eHjP0uWFVu6rtg1X+8SUACEpBA7wkoIL1n7iNKQAISmCgE6sXZ9XqJ93ptw4S8BqSeNlU/uXz0zwGpF4LXV0bqxecXXXTRqGMe+Yv/QgstNOp/u+eee7oPMpxzzjk7kRn5WR5jg6qfAzL33HN3b7l7xRVXjPrXVTrqJ6nXzwGZeeaZu/+9Xv9x4IEHlltuuaUst9xy3f923333lcUXX7wTtfquYd4kIAEJSKD3BBSQ3jP3ESUgAQlMFAL1F+362Rf18zGqDNRP/t5oo43e9W14J2SJt956q6y66qrl7rvvHuOT0B955JFOAEZ/t6mRn9I+8p21XnrppU4KHn/88e7DAeeaa64xqlVRWWmllUb9byPfxap+4OHIT0Kvr3zU07MOPvjgUV9X73eZZZYp9f9W4ag8TjrppFK71k9CH99pZROSj/clAQlIQAL/R0AB8dkgAQlIYIAI1M+3qK+EPPnkk91nYbzX07EmBILnnnuuk4/6LlT1lYp6ClS9LmXkqw8jH2NsAakissACC7xrha233rr88Ic/HOPfn3nmmd3nfdTjq5+MXi9O32OPPbpTrka/PfbYY2WvvfYq11xzTcdj9dVX7z4nZeGFF54Qh+x9SEACEpBAQEABCaAZkYAEJCABCUhAAhKQgAQyAgpIxs2UBCQgAQlIQAISkIAEJBAQUEACaEYkIAEJSEACEpCABCQggYyAApJxMyUBCUhAAhKQgAQkIAEJBAQUkACaEQlIQAISkIAEJCABCUggI6CAZNxMSUACEpCABCQgAQlIQAIBAQUkgGZEAhKQgAQkIAEJSEACEsgIKCAZN1MSkIAEJCABCUhAAhKQQEBAAQmgTczI5FOO+QnAE/OxJsV9Tzfl1JPiYXv6mNNPMdjHOPXwqXrKs9cPNvXwKXr9kD1/vGFlzA/r63mBifyAI8qIifwIk/7uX3vr35O+xERs8Npbr0/Ee5/0d/3yv1+b9CUmcoN/vTHYx/jmG49PZIKDffcKCGxfBQQ2SFBHAQmggSIKCGiMsIoCEoIDxRQQ0BhhFQUkBNdITAGBDa2AwAYJ6iggATRQRAEBjRFWUUBCcKCYAgIaI6yigITgGokpILChFRDYIEEdBSSABoooIKAxwioKSAgOFFNAQGOEVRSQEFwjMQUENrQCAhskqKOABNBAEQUENEZYRQEJwYFiCghojLCKAhKCaySmgMCGVkBggwR1FJAAGiiigIDGCKsoICE4UEwBAY0RVlFAQnCNxBQQ2NAKCGyQoI4CEkADRRQQ0BhhFQUkBAeKKSCgMcIqCkgIrpGYAgIbWgGBDRLUUUACaKCIAgIaI6yigITgQDEFBDRGWEUBCcE1ElNAYEMrILBBgjoKSAANFFFAQGOEVRSQEBwopoCAxgirKCAhuEZiCghsaAUENkhQRwEJoIEiCghojLCKAhKCA8UUENAYYRUFJATXSEwBgQ2tgMAGCeooIAE0UEQBAY0RVlFAQnCgmAICGiOsooCE4BqJKSCwoRUQ2CBBHQUkgAaKKCCgMcIqCkgIDhRTQEBjhFUUkBBcIzEFBDa0AgIbJKijgATQQBEFBDRGWEUBCcGBYgoIaIywigISgmskpoDAhlZAYIMEdRSQABooooCAxgirKCAhOFBMAQGNEVZRQEJwjcQUENjQCghskKCOAhJAA0UUENAYYRUFJAQHiikgoDHCKgpICK6RmAICG1oBgQ0S1FFAAmigiAICGiOsooCE4EAxBQQ0RlhFAQnBNRJTQGBDKyCwQYI6CkgADRRRQEBjhFUUkBAcKKaAgMYIqyggIbhGYgoIbGgFBDZIUEcBCaCBIgoIaIywigISggPFFBDQGGEVBSQE10hMAYENrYDABgnqKCABNFBEAQGNEVZRQEJwoJgCAhojrKKAhOAaiSkgsKEVENggQR0FJIAGiiggoDHCKgpICA4UU0BAY4RVFJAQXCMxBQQ2tAICGySoo4AE0EARBQQ0RlhFAQnBgWIKCGiMsIoCEoJrJKaAwIZWQGCDBHUUkAAaKKKAgMYIqyggIThQTAEBjRFWUUBCcI3EFBDY0AoIbJCgjgISQANFFBDQGGEVBSQEB4opIKAxwioKSAiukZgCAhtaAYENEtRRQAJooIgCAhojrKKAhOBAMQUENEZYRQEJwTUSU0BgQysgsEGCOgpIAA0UUUBAY4RVFJAQHCimgIDGCKsoICG4RmIKCGxoBQQ2SFBHAQmggSIKCGiMsIoCEoIDxRQQ0BhhFQUkBNdITAGBDa2AwAYJ6iggATRQRAEBjRFWUUBCcKCYAgIaI6yigITgGokpILChFRDYIEEdBSSABoooIKAxwioKSAgOFFNAQGOEVRSQEFwjMQUENrQCAhskqKOABNBAEQUENEZYRQEJwYFiCghojLCKAhKCaySmgMCGVkBggwR1FJAAGiiigIDGCKsoICE4UEwBAY0RVlFAQnCNxBQQ2NAKCGyQoI4CEkADRRQQ0BhhFQUkBAeKKSCgMcIqCkgIrpGYAgIbWgGBDRLUUUACaKCIAgIaI6yigITgQDEFBDRGWEUBCcE1ElNAYEMrILBBgjoKSAANFFFAQGOEVRSQEBwopoCAxgirKCAhuEZiCghsaAUENkhQRwEJoIEiCghojLCKAhKCA8UUENAYYRUFJATXSEwBgQ2tgMAGCeooIAE0UEQBAY0RVlFAQnCgmAICGiOsooCE4BqJKSCwoRUQ2CBBHQUkgAaKKCCgMcIqCkgIDhRTQEBjhFUUkBBcIzEFBDa0AgIbJKijgATQQBEFBDRGWEUBCcGBYgoIaIywigISgmskpoDAhlZAYIMEdRSQABooooCAxgirKCAhOFBMAQGNEVZRQEJwjcQUENjQCghskKCOAhJAA0UUENAYYRUFJAQHiikgoDHCKgpICK6RmAICG1oBgQ0S1FFAAmigiAICGiOsooCE4EAxBQQ0RlhFAQnBNRJTQGBDKyCwQYI6CkgADRRRQEBjhFUUkBAcKKaAgMYIqyggIbhGYgoIbGgFBDZIUEcBCaCBIgoIaIywigISggPFFBDQGGEVBSQE10hMAYENrYDABgnqKCABNFBEAQGNEVZRQEJwoJgCAhojrKKAhOAaiSkgsKEVENggQR0FJIAGiiggoDHCKgpICA4UU0BAY4RVFJAQXCMxBQQ2tAICGySoo4AE0EARBQQ0RlhFAQnBgWIKCGiMsIoCEoJrJKaAwIZWQGCDBHUUkAAaKKKAgMYIqyggIThQTAEBjRFWUUBCcI3EFBDY0AoIbJCgjgISQANFFBDQGGEVBSQEB4opIKAxwioKSAiukZgCAhtaAYENEtRRQAJooIgCAhojrKKAhOBAMQUENEZYRQEJwTUSU0BgQysgsEGCOgpIAA0UUUBAY4RVFJAQHCimgIDGCKsoICG4RmIKCGxoBQQ2SFBHAQmggSIKCGiMsIoCEoIDxRQQ0BhhFQUkBNdITAGBDa2AwAYJ6iggATRQRAEBjRFWUUBCcKCYAgIaI6yigITgGokpILChFRDYIEEdBSSABoooIKAxwioKSAgOFFNAQGOEVRSQEFwjMQUENrQCAhskqKOABNBAEQUENEZYRQEJwYFiCghojLCKAhKCaySmgMCGVkBggwR1FJAAGiiigIDGCKsoICE4UEwBAY0RVlFAQnCNxBQQ2NAKCGyQoI4CEkADRRQQ0BhhFQUkBAeKKSCgMcIqCkgIrpGYAgIbWgGBDRLUUUACaKCIAgIaI6yigITgQDEFBDRGWEUBCcE1ElNAYEMrILBBgjoKSAANFFFAQGOEVRSQEBwopoCAxgirKCAhuEZiCghsaAUENkhQRwEJoIEiCghojLCKAhKCA8UUENAYYRUFJATXSEwBgQ2tgMAGCeooIAE0UEQBAY0RVlFAQnCgmAICGiOsooCE4BqJKSCwoRUQ2CBBHQUkgAaKKCCgMcIqCkgIDhRTQEBjhFUUkBBcIzEFBDa0AgIbJKijgATQQBEFBDRGWEUBCcGBYgoIaIywigISgmskpoDAhlZAYIMEdRSQABooooCAxgirKCAhOFBMAQGNEVZRQEJwjcQUENjQCghskKCOAhJAA0UUENAYYRUFJAQHiikgoDHCKgpICK6RmAICG1oBgQ0S1FFAAmigiAICGiOsooCE4EAxBQQ0RlhFAQnBNRJTQGBDKyCwQYI6CkgADRRRQEBjhFUUkBAcKKaAgMYIqyggIbhGYgoIbGgFBDZIUEcBCaCBIgoIaIywigISggPFFBDQGGEVBSQE10hMAYENrYDABgnqKCABNFBEAQGNEVZRQEJwoJgCAhojrKKAhOAaiSkgsKEVENggQR0FJIAGiiggoDHCKgpICA4UU0BAY4RVFJAQXCMxBQQ2tAICGySoo4AE0EARBQQ0RlhFAQnBgWIKCGiMsIoCEoJrJKaAwIZWQGCDBHUUkAAaKKKAgMYIqyggIThQTAEBjRFWUUBCcI3EFBDY0AoIbJCgjgISQANFFBDQGGEVBSQEB4opIKAxwioKSAiukZgCAhtaAYENEtRRQAJooIgCAhojrKKAhOBAMQUENEZYRQEJwTUSU0BgQysgsEGCOgpIAA0UUUBAY4RVFJAQHCimgIDGCKsoICG4RmIKCGxoBQQ2SFBHAQmggSIKCGiMsIoCEoIDxRQQ0BhhFQUkBNdITAGBDa2AwAYJ6iggATRQRAEBjRFWUUBCcKCYAgIaI6yigITgGokpILChFRDYIEEdBSSABoooIKAxwioKSAgOFFNAQGOEVRSQEFwjMQUENrQCAhskqKOABNBAEQUENEZYRQEJwYFiCghojLCKAhKCaySmgMCGVkBggwR1FJAAGiiigIDGCKsoICE4UEwBAY0RVlFAQnCNxBQQ2NAKCGyQoI4CEkADRRQQ0BhhFQUkBAeKKSCgMcIqCkgIrpGYAgIbWgGBDRLUUUACaKCIAgIaI6yigITgQDEFBDRGWEUBCcE1ElNAYEMrILBBgjoKSAANFFFAQGOEVRSQEBwopoCAxgirKCAhuEZiCghsaAUENkhQRwEJoIEiCghojLCKAhKCA8UUENAYYRUFJATXSEwBgQ2tgMAGCeooIAE0UEQBAY0RVlFAQnCgmAICGiOsooCE4BqJKSCwoRUQ2CBBHQUkgAaKKCCgMcIqCkgIDhRTQEBjhFUUkBBcIzEFBDa0AgIbJKijgATQQBEFBDRGWEUBCcGBYgoIaIywigISgmskpoDAhlZAYIMEdRSQABooooCAxgirKCAhOFBMAQGNEVZRQEJwjcQUENjQCghskKCOAhJAA0UUENAYYRUFJAQHiikgoDHCKgpICK6RmAICG1oBgQ0S1FFAAmigiAICGiOsooCE4EAxBQQ0RlhFAQnBNRJTQGBDKyCwQYI6CkgADRRRQEBjhFUUkBAcKKaAgMYIqyggIbhGYgoIbGgFBDZIUEcBCaCBIgoIaIywigISggPFFBDQGGEVBSQE10hMAYENrYDABgnqKCABNFBEAQGNEVZRQEJwoJgCAhojrKKAhOAaiSkgsKEVENggQR0FJIAGiiggoDHCKgpICA4UU0BAY4RVFJAQXCMxBQQ2tAICGySoo4AE0EARBQQ0RlhFAQnBgWIKCGiMsIoCEoJrJKaAwIZWQGCDBHUUkAAaKKKAgMYIqyggIThQTAEBjRFWUUBCcI3EFBDY0AoIbJCgjgISQANFFBDQGGEVBSQEB4opIKAxwioKSAiukZgCAhtaAYENEtRRQAJooIgCAhojrKKAhOBAMQUENEZYRQEJwTUSU0BgQysgsEGCOgpIAA0UUUBAY4RVFJAQHCimgIDGCKsoICG4RmIKCGxoBQQ2SFBHAQmggSIKCGiMsIoCEoIDxRQQ0BhhFQUkBNdITAGBDa2AwAYJ6iggATRQRAEBjRFWUUBCcKCYAgIaI6yigITgGokpILChFRDYIEEdBSSABoooIKAxwioKSAgOFFNAQGOEVRSQEFwjMQUENrQCAhskqKOABNBAEQUENEZYRQEJwYFiCghojLCKAhKCaySmgMCGVkBggwR1FJAAGiiigIDGCKsoICE4UEwBAY0RVlFAQnCNxBQQ2NAKCGyQoI4CEkADRRQQ0BhhFQUkBAeKKSCgMcIqCkgIrpGYAgIbWgGBDRLUUUACaKCIAgIaI6yigITgQDEFBDRGWEUBCcE1ElNAYEMrILBBgjoKSAANFFFAQGOEVRSQEBwopoCAxgirKCAhuEZiCghsaAUENkhQRwEJoIEiCghojLCKAhKCA8UUENAYYRUFJATXSEwBgQ2tgMAGCeooIAE0UEQBAY0RVlFAQnCgmAICGiOsooCE4BqJKSCwoRUQ2CBBHQUkgAaKKCCgMcIqCkgIDhRTQEBjhFUUkBBcIzEFBDa0AgIbJKijgATQQBEFBDRGWEUBCcGBYgoIaIywigISgmskpoDAhlZAYIMEdRSQABooooCAxgirKCAhOFBMAQGNEVZRQEJwjcQUENjQCghskKCOAhJAA0UUENAYYRUFJAQHiikgoDHCKgpICK6RmAICG1oBgQ0S1FFAAmigiAICGiOsooCE4EAxBQQ0RlhFAQnBNRJTQGBDKyCwQYI6CkgADRRRQEBjhFUUkBAcKKaAgMYIqyggIbhGYgpIOPSzzz5bzjnnnHLjjTeWp556qruXOeecs6y88splm222KbPNNlt0zwpIhA0VUkBQcwy5jAIyZGS4gAKCm2TIhRSQISPDBRQQ3CSoQgpIMMctt9xS1l133TLttNOWtdZaq8wxxxzdvTz99NPl2muvLa+88kq5+uqry3LLLTfke1dAhowMF1BAcJMMqZACMiRcyC9WQJCzDKmUAjIkXMgvVkCQs2BKKSDBFCuuuGJZaqmlyumnn16GDRs2xj2MGDGi7LjjjuXOO+/sXh0Z1+31118v9Z/Rb++fZbG33WdQERuZbsqpsd0mVDEFZEKRnDT3o4BMGu4T8lEVkAlJc9LclwIyabhPyEdVQCYkzcG7LwUk2HSaaaYpt99+e1lsscXeMX3vvfeWZZZZprz66qvjvPcjjjiiHHnkkWN8zbDJpi+TDZ8xaNUfEQWkP3YaV8uph0/V/wcxjiNQQPp/XgWk/zdUQPp/QwWk/zecmEeggAR0F1hggU4cttpqq3dMn3feeeWwww4rDz300Djv3VdAAvh9EPEVkD4YSQHp75HG014B6f95FZD+31AB6f8NJ+YRKCAB3VNPPbXss88+ZYcddihrrrnm264BOfPMM8uJJ55Ydt555yHfu9eADBkZLqCA4CYZUiFfARkSLuQXKyDIWYZUSgEZEi7kFysgyFkwpRSQcIqLL764nHzyyeW2224rb731Vncvw4cPL8suu2zZe++9y+abbx7dswISYUOFFBDUHEMuo4AMGRkuoIDgJhlyIQVkyMhwAQUENwmqkALyX87x73//u9S35K23WWedtUwxxRT/1T0qIP8VPkRYAUHMEJdQQGJ0mKACgpkiLqKAxOgwQQUEMwWyiAICm0UBgQ0S1FFAAmigiAICGiOsooCE4EAxBQQ0RlhFAQnBNRJTQGBDKyCwQYI6CkgADRRRQEBjhFUUkBAcKKaAgMYIqyggIbhGYgoIbGgFBDZIUEcBCaCBIgoIaIywigISggPFFBDQGGEVBSQE10hMAYENrYDABgnqKCABNFBEAQGNEVZRQEJwoJgCAhojrKKAhOAaiSkgsKEVENggQR0FJIAGiiggoDHCKgpICA4UU0BAY4RVFJAQXCMxBQQ2tAICGySoo4AE0EARBQQ0RlhFAQnBgWIKCGiMsIoCEoJrJKaAwIZWQGCDBHUUkAAaKKKAgMYIqyggIThQTAEBjRFWUUBCcI3EFBDY0AoIbJCgjgISQANFFBDQGGEVBSQEB4opIKAxwioKSAiukZgCAhtaAYENEtRRQAJooIgCAhojrKKAhOBAMQUENEZYRQEJwTUSU0BgQysgsEGCOgpIAA0UUUBAY4RVFJAQHCimgIDGCKsoICG4RmIKCGxoBQQ2SFBHAQmggSIKCGiMsIoCEoIDxRQQ0BhhFQUkBNdITAGBDa2AwAYJ6iggATRQRAEBjRFWUUBCcKCYAgIaI6yigITgGokpILChFRDYIEEdBSSABoooIKAxwioKSAgOFFNAQGOEVRSQEFwjMQUENrQCAhskqKOABNBAEQUENEZYRQEJwYFiCghojLCKAhKCaySmgMCGVkBggwR1FJAAGiiigIDGCKsoICE4UEwBAY0RVlFAQnCNxBQQ2NAKCGyQoI4CEkADRRQQ0BhhFQUkBAeKKSCgMcIqCkgIrpGYAgIbWgGBDRLUUUACaKCIAgIaI6yigITgQDEFBDRGWEUBCcE1ElNAYEMrILBBgjoKSAANFFFAQGOEVRSQEBwopoCAxgirKCAhuEZiCghsaAUENkhQRwEJoIEiCghojLCKAhKCA8UUENAYYRUFJATXSEwBgQ2tgMAGCeooIAE0UEQBAY0RVlFAQnCgmAICGiOsooCE4BqJKSCwoRUQ2CBBHQUkgAaKKCCgMcIqCkgIDhRTQEBjhFUUkBBcIzEFBDa0AgIbJKijgATQQBEFBDRGWEUBCcGBYgoIaIywigISgmskpoDAhlZAYIMEdRSQABoGQ8OTAAAYv0lEQVQoooCAxgirKCAhOFBMAQGNEVZRQEJwjcQUENjQCghskKCOAhJAA0UUENAYYRUFJAQHiikgoDHCKgpICK6RmAICG1oBgQ0S1FFAAmigiAICGiOsooCE4EAxBQQ0RlhFAQnBNRJTQGBDKyCwQYI6CkgADRRRQEBjhFUUkBAcKKaAgMYIqyggIbhGYgoIbGgFBDZIUEcBCaCBIgoIaIywigISggPFFBDQGGEVBSQE10hMAYENrYDABgnqKCABNFBEAQGNEVZRQEJwoJgCAhojrKKAhOAaiSkgsKEVENggQR0FJIAGiiggoDHCKgpICA4UU0BAY4RVFJAQXCMxBQQ2tAICGySoo4AE0EARBQQ0RlhFAQnBgWIKCGiMsIoCEoJrJKaAwIZWQGCDBHUUkAAaKKKAgMYIqyggIThQTAEBjRFWUUBCcI3EFBDY0AoIbJCgjgISQANFFBDQGGEVBSQEB4opIKAxwioKSAiukZgCAhtaAYENEtRRQAJooIgCAhojrKKAhOBAMQUENEZYRQEJwTUSU0BgQysgsEGCOgpIAA0UUUBAY4RVFJAQHCimgIDGCKsoICG4RmIKCGxoBQQ2SFBHAQmggSIKCGiMsIoCEoIDxRQQ0BhhFQUkBNdITAGBDa2AwAYJ6iggATRQRAEBjRFWUUBCcKCYAgIaI6yigITgGokpILChFRDYIEEdBSSABoooIKAxwioKSAgOFFNAQGOEVRSQEFwjMQUENrQCAhskqKOABNBAEQUENEZYRQEJwYFiCghojLCKAhKCaySmgMCGVkBggwR1FJAAGiiigIDGCKsoICE4UEwBAY0RVlFAQnCNxBQQ2NAKCGyQoI4CEkADRRQQ0BhhFQUkBAeKKSCgMcIqCkgIrpGYAgIbWgGBDRLUUUACaKCIAgIaI6yigITgQDEFBDRGWEUBCcE1ElNAYEMrILBBgjoKSAANFFFAQGOEVRSQEBwopoCAxgirKCAhuEZiCghsaAUENkhQRwEJoIEiCghojLCKAhKCA8UUENAYYRUFJATXSEwBgQ2tgMAGCeooIAE0UEQBAY0RVlFAQnCgmAICGiOsooCE4BqJKSCwoRUQ2CBBHQUkgAaKKCCgMcIqCkgIDhRTQEBjhFUUkBBcIzEFBDa0AgIbJKijgATQQBEFBDRGWEUBCcGBYgoIaIywigISgmskpoDAhlZAYIMEdRSQABooooCAxgirKCAhOFBMAQGNEVZRQEJwjcQUENjQCghskKCOAhJAA0UUENAYYRUFJAQHiikgoDHCKgpICK6RmAICG1oBgQ0S1FFAAmigiAICGiOsooCE4EAxBQQ0RlhFAQnBNRJTQGBDKyCwQYI6CkgADRRRQEBjhFUUkBAcKKaAgMYIqyggIbhGYgoIbGgFBDZIUEcBCaCBIgoIaIywigISggPFFBDQGGEVBSQE10hMAYENrYDABgnqKCABNFBEAQGNEVZRQEJwoJgCAhojrKKAhOAaiSkgsKEVENggQR0FJIAGiiggoDHCKgpICA4UU0BAY4RVFJAQXCMxBQQ2tAICGySoo4AE0EARBQQ0RlhFAQnBgWIKCGiMsIoCEoJrJKaAwIZWQGCDBHUUkAAaKKKAgMYIqyggIThQTAEBjRFWUUBCcI3EFBDY0AoIbJCgjgISQANFFBDQGGEVBSQEB4opIKAxwioKSAiukZgCAhtaAYENEtRRQAJooIgCAhojrKKAhOBAMQUENEZYRQEJwTUSU0BgQysgsEGCOgpIAA0UUUBAY4RVFJAQHCimgIDGCKsoICG4RmIKCGxoBQQ2SFBHAQmggSIKCGiMsIoCEoIDxRQQ0BhhFQUkBNdITAGBDa2AwAYJ6iggATRQRAEBjRFWUUBCcKCYAgIaI6yigITgGokpILChFRDYIEEdBSSABoooIKAxwioKSAgOFFNAQGOEVRSQEFwjMQUENrQCAhskqKOABNBAEQUENEZYRQEJwYFiCghojLCKAhKCaySmgMCGVkBggwR1FJAAGiiigIDGCKsoICE4UEwBAY0RVlFAQnCNxBQQ2NAKCGyQoI4CEkADRRQQ0BhhFQUkBAeKKSCgMcIqCkgIrpGYAgIbWgGBDRLUUUACaKCIAgIaI6yigITgQDEFBDRGWEUBCcE1ElNAYEMrILBBgjoKSAANFFFAQGOEVRSQEBwopoCAxgirKCAhuEZiCghsaAUENkhQRwEJoIEiCghojLCKAhKCA8UUENAYYRUFJATXSEwBgQ2tgMAGCeooIAE0UEQBAY0RVlFAQnCgmAICGiOsooCE4BqJKSCwoRUQ2CBBHQUkgAaKKCCgMcIqCkgIDhRTQEBjhFUUkBBcIzEFBDa0AgIbJKijgATQQBEFBDRGWEUBCcGBYgoIaIywigISgmskpoDAhlZAYIMEdRSQABooooCAxgirKCAhOFBMAQGNEVZRQEJwjcQUENjQCghskKCOAhJAA0UUENAYYRUFJAQHiikgoDHCKgpICK6RmAICG1oBgQ0S1FFAAmigiAICGiOsooCE4EAxBQQ0RlhFAQnBNRJTQGBDKyCwQYI6CkgADRRRQEBjhFUUkBAcKKaAgMYIqyggIbhGYgoIbGgFBDZIUEcBCaCBIgoIaIywigISggPFFBDQGGEVBSQE10hMAYENrYDABgnqKCABNFBEAQGNEVZRQEJwoJgCAhojrKKAhOAaiSkgsKEVENggQR0FJIAGiiggoDHCKgpICA4UU0BAY4RVFJAQXCMxBQQ2tAICGySoo4AE0EARBQQ0RlhFAQnBgWIKCGiMsIoCEoJrJKaAwIZWQGCDBHUUkAAaKKKAgMYIqyggIThQTAEBjRFWUUBCcI3EFBDY0AoIbJCgjgISQANFFBDQGGEVBSQEB4opIKAxwioKSAiukZgCAhtaAYENEtRRQAJooIgCAhojrKKAhOBAMQUENEZYRQEJwTUSU0BgQysgsEGCOgpIAA0UUUBAY4RVFJAQHCimgIDGCKsoICG4RmIKCGxoBQQ2SFBHAQmggSIKCGiMsIoCEoIDxRQQ0BhhFQUkBNdITAGBDa2AwAYJ6iggATRQRAEBjRFWUUBCcKCYAgIaI6yigITgGokpILChFRDYIEEdBSSABoooIKAxwioKSAgOFFNAQGOEVRSQEFwjMQUENrQCAhskqKOABNBAEQUENEZYRQEJwYFiCghojLCKAhKCaySmgMCGVkBggwR1FJAAGiiigIDGCKsoICE4UEwBAY0RVlFAQnCNxBQQ2NAKCGyQoI4CEkADRRQQ0BhhFQUkBAeKKSCgMcIqCkgIrpGYAgIbWgGBDRLUUUACaKCIAgIaI6yigITgQDEFBDRGWEUBCcE1ElNAYEMrILBBgjoKSAANFFFAQGOEVRSQEBwopoCAxgirKCAhuEZiCghsaAUENkhQRwEJoIEiCghojLCKAhKCA8UUENAYYRUFJATXSEwBgQ2tgMAGCeooIAE0UEQBAY0RVlFAQnCgmAICGiOsooCE4BqJKSCwoRUQ2CBBHQUkgAaKKCCgMcIqCkgIDhRTQEBjhFUUkBBcIzEFBDa0AgIbJKijgATQQBEFBDRGWEUBCcGBYgoIaIywigISgmskpoDAhlZAYIMEdRSQABooooCAxgirKCAhOFBMAQGNEVZRQEJwjcQUENjQCghskKCOAhJAA0UUENAYYRUFJAQHiikgoDHCKgpICK6RmAICG1oBgQ0S1FFAAmigiAICGiOsooCE4EAxBQQ0RlhFAQnBNRJTQGBDKyCwQYI6CkgADRRRQEBjhFUUkBAcKKaAgMYIqyggIbhGYgoIbGgFBDZIUEcBCaCBIgoIaIywigISggPFFBDQGGEVBSQE10hMAYENrYDABgnqKCABNFBEAQGNEVZRQEJwoJgCAhojrKKAhOAaiSkgsKEVENggQR0FJIAGiiggoDHCKgpICA4UU0BAY4RVFJAQXCMxBQQ2tAICGySoo4AE0EARBQQ0RlhFAQnBgWIKCGiMsIoCEoJrJKaAwIZWQGCDBHUUkAAaKKKAgMYIqyggIThQTAEBjRFWUUBCcI3EFBDY0AoIbJCgjgISQANFFBDQGGEVBSQEB4opIKAxwioKSAiukZgCAhtaAYENEtRRQAJooIgCAhojrKKAhOBAMQUENEZYRQEJwTUSU0BgQysgsEGCOgpIAA0UUUBAY4RVFJAQHCimgIDGCKsoICG4RmIKCGxoBQQ2SFBHAQmggSIKCGiMsIoCEoIDxRQQ0BhhFQUkBNdITAGBDa2AwAYJ6iggATRQRAEBjRFWUUBCcKCYAgIaI6yigITgGokpILChFRDYIEEdBSSABoooIKAxwioKSAgOFFNAQGOEVRSQEFwjMQUENrQCAhskqKOABNBAEQUENEZYRQEJwYFiCghojLCKAhKCaySmgDQy9Dsd5uuvv16OPfbYctBBB5WppppqIEkM+jF6fP3/tHXD/t5w0Per6wz6MXp8/f3fYAvP0f5f6O1HoIAM4qrv8ZhefPHFMtNMM5UXXnihzDjjjO8x1V9fNujH6PH11/Pxndq6YX9vOOj71XUG/Rg9vv7+b7CF52j/L6SADOKG8TEN+jfdFr4pDfqGg358Pkfjb1+YoM9RzBRxkUHfcNCPr4Xvo/GTGxz0FRDwOBO7mt+UJjbhiX//g77hoB9fCz84B33DQT8+n6MT//v4xH4En6MTm7D3nxBQQBJqA5Lxm1L/DznoGw768fnLnf8N9gOBQf/v0OPrh2fhuDsO+ob9v5CnYA3ihvExDfqFdxXMoB+jxxc//TFBN8RMERUZ9P38Pho9LVAhn6OoOSzzvwR8BcSnggQkIAEJSEACEpCABCTQMwIKSM9Q+0ASkIAEJCABCUhAAhKQgALic0ACEpCABCQgAQlIQAIS6BkBBaRnqH0gCUhAAhKQgAQkIAEJSEAB8TkgAQlIQAISkIAEJCABCfSMgALSM9S8Bzr11FPLN77xjfLUU0+VpZZaqnznO98pK6ywAq9o0Oi6667rju22224rTz75ZLn88svLxhtvHNwTM3LssceWyy67rNx7771lmmmmKSuvvHI5/vjjy6KLLsosPMRW3/ve90r956GHHuqSiy++eDnssMPKeuutN8R76o8vP+6448pBBx1U9thjj/Ktb32rP0qPp+URRxxRjjzyyDG+qj4/63N2UG6PP/54OeCAA8ovfvGL8sorr5SFF164/OAHPyjLLbdc3x/i/PPPXx5++OG3HcfOO+9c6s+OQbi99dZbpT5Pzz///O7n4Ac/+MGyzTbblEMOOaQMGzZsEA6xvPTSS+XQQw/tfgb+/e9/L8sss0z59re/XZZffvm+PL7x/WwfMWJEOfzww8uZZ55Znn/++bLKKqt0P0sWWWSRvjzeQS6tgAzyuuM4tosvvrhstdVW5fTTTy8f+9jHul96Lr300nLfffeV2Wefve+p1F8I/vjHP5Zll122fOYznxk4AfnkJz9Ztthii+6HyJtvvlkOPvjgcvfdd5c///nPZbrppuv7/a644ooyfPjw7odG/YFy7rnndkJ5++23dzIySLdbbrmlbL755mXGGWcsa6yxxkAJyE9+8pPy61//etRck08+eZl11lkHYr7nnnuu+2WubrbTTjuV2Wabrfz1r38tCy20UPdPv9+eeeaZUn9BH3mr31/WXnvt8tvf/rasvvrq/X54Xf9jjjmmnHTSSd33l/p95dZbby1f/vKXy9FHH1123333gTjGz33uc93PhvpLeBWsKlsnn3xy97Nirrnm6rtjHN/P9vqHuPoHurrpAgss0MnXXXfd1R3v1FNP3XfHO8iFFZBBXnccx1alo/7y+t3vfrf7qv/85z9lnnnmKbvttls58MADB4pK/UvWoL0CMvZA9ZeFKo6///3vy2qrrTZQ+408mJlnnrmTkO22225gju/ll18uH/3oR8tpp51WjjrqqLL00ksPlID89Kc/LXfcccfA7DX6gdTvk/WPHNdff/1AHt/YB7XnnnuWK6+8spOsQXl1YMMNNyxzzDFHOfvss0cd7qabbtq9qlx/Ue/326uvvlpmmGGG8rOf/axssMEGow6n/mGuvppcv+f0823sn+31j1VVsvbZZ5+y7777dof2wgsvdBv/8Ic/7P5o541DQAHhbNGzJm+88UaZdtppS/3r5OinJW299dbdS5b1m9Ug3VoQkPvvv797taD+pWeJJZYYpPm6v8LWV+fq87O+AvLhD394YI6vHlMVq/oXyfpX5UETkCqMM800U/eXx5VWWqn7y+S88847EPvV5+G6665bHnvssU7861+T6+lJX/3qVwfi+EY/iPozo/5it/fee3evtg7Krb4C8v3vf79cc8015UMf+lD505/+VNZZZ53uVZEvfOELfX+Y9fSr+spqfRVyzTXXHHU8q666aqmvRv7ud7/r62Mc+2f7gw8+2L36WH9O1O+lI2//8z//0/3/9dQzbxwCCghni541eeKJJ7ofljfccEP3S8HI2/7779/9IL3pppt61qUXDzToAlJfvfrUpz7VyeMf/vCHXiDtyWNUmarPz9dee61MP/305YILLijrr79+Tx67Fw9y0UUXdad61FOw6i/ogyYg9VSJ+gpPve6jXodVrwep10zU00HqX2X7/TbydI76S/lmm23W7Viv4amntVaxHKTbJZdcUj7/+c+XRx55pBORQbnV751VqE444YTulM/6x47632S9HmtQbvX6wCmnnLL7/llfCbjwwgu752e9Xqmect3Pt7F/ttffaeo1H/V3nA984AOjDq2e4lq/tp567o1DQAHhbNGzJgpIz1D35IHq+ef1l70qH3PPPXdPHrMXD1L/6lp/4akvoddX684666xOkAfhFZBHH320u1D5V7/6VVlyySU7nIMmIGM/R6ogzzfffN1flwfhNLr6S13dsP7SM/JWrxuoInLjjTf24j+Rnj1GfaWnHm+9NmuQbvWPAPvtt193ame9BqSeLlhPNavP0UGRyAceeKBsu+22pV68XSWrnvJZX+2pb9Byzz339PWcCkhfz1cUkP7eL2rvKVgRNmRo11137U6Zqz9c6gV3g3xba621upfXzzjjjL4/zHptxCabbNL9QjDyVv/6Wn+gTjbZZOX1118f49/1/QH/7wHU687qjvVUrH6/VZmqF2VXMR55qxf61vPq6ys9g3Kr74S14IILdu+69+lPf3pQDqs7jnrdY72WZ5dddhl1XHW/ev3HIL1bWz24f/3rX+XFF1/sXhmoF6bXVyevuuqqvt7TU7D6ej4FpL/ny9vXi9DrW+7Wt96tt/pSdD03u/5C60XoOddeJevFdvUNA+rF9fU83hbeYvATn/hE9xytFxP2+62emz32W5zWd99ZbLHFurd1HbTreOpe9Reeul9929NBeIehekpSfSVr9IvQ99prr+4U1tFfFen352rdq0p/PdZ63cAg3WaZZZZOGOuryCNvVY7rWyn/5S9/GaRDHXUs9d3b6h+r6mln22+/fV8f47tdhF4vQK8Xotdbla76Bi1ehM6b2ldAeJv0pFE9F7K+xFx/sFQRqW/DW8/zrX/1qeeJ9vut/rJTL8yut/pWmfUl9fp2mfWC30G4CLZe7FrP6a2vfoz+2R/1gt/6Di79fqvnYNd3aalb1V/W67HWt1e8+uqru786D+Jt0E7Bqr8EbLTRRt1pV/W0z/re/PUUl/p2mPUta/v9Vk+1qufX12tb6jnmN998c3cBer2oeRAuYB75h6n6y+qWW25Z6mfVDNqtfuZHvUC7/hysp2DVi5frL+X1lKX6/WYQbvV7Zv2DVf05UX8m1lPO6vVLVZynmGKKvjvE8f1sr7vV5+rob8N75513+ja8wKUVEOAovapU34J35AcR1neIOOWUU7rPBBmEW31VoArH2LcqXYPwF/R3exvM+pe7+kO132/1GoFrr722u3i5SlW9TqK+MjCo8lH3GjQBqW95WU8N/Mc//tEJR33nnXqB7yB8RsbI/77q29JWWa5vTVt/Ua8XpA/Su2DVd4eq13/Ui5XrdQODdhv7Q/rqBfZVtuqHntZrXgbhVv+wWJ+j9d3a6h/g6tsM1/8O6/fVfryN72f7yA8irH8IqNed1e879W3OB/H524/7jd5ZAen3Be0vAQlIQAISkIAEJCCBPiKggPTRWFaVgAQkIAEJSEACEpBAvxNQQPp9QftLQAISkIAEJCABCUigjwgoIH00llUlIAEJSEACEpCABCTQ7wQUkH5f0P4SkIAEJCABCUhAAhLoIwIKSB+NZVUJSEACEpCABCQgAQn0OwEFpN8XtL8EJCABCUhAAhKQgAT6iIAC0kdjWVUCEpCABCQgAQlIQAL9TkAB6fcF7S8BCUhAAhKQgAQkIIE+IqCA9NFYVpWABCQgAQlIQAISkEC/E1BA+n1B+0tAAhKQgAQkIAEJSKCPCCggfTSWVSUgAQlIQAISkIAEJNDvBP4fIqSBIS8WXLsAAAAASUVORK5CYII=\" width=\"800\">"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/plain": [
       "<matplotlib.animation.FuncAnimation at 0x7fe7e5b36fd0>"
      ]
     },
     "execution_count": 8,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "animate_heatmap(u, t, sampling=1)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "4a389c6b",
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "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.9.7"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}