MotPlotLib Experimenting.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "id": "8c1943d6",
   "metadata": {},
   "source": [
    "<h1>Gaussian Random Data</h1>"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "f1930738",
   "metadata": {},
   "outputs": [],
   "source": [
    "import numpy as np\n",
    "x = np.random.randn(10000)"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "96a787d9",
   "metadata": {},
   "source": [
    "<h1>Using the Artist layer</h1>"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "id": "f3c7d7e7",
   "metadata": {},
   "outputs": [],
   "source": [
    "%matplotlib inline\n",
    "from matplotlib.backends.backend_agg import FigureCanvasAgg as FigureCanvas\n",
    "from matplotlib.figure import Figure\n",
    "\n",
    "fig = Figure()\n",
    "canvas = FigureCanvas(fig)\n",
    "\n",
    "ax = fig.add_subplot(111)\n",
    "ax.hist(x, 100)\n",
    "\n",
    "ax.set_title(r'Normal distribution with $\\mu=0, \\sigma=1$')\n",
    "fig.savefig('matplotlib histogram.png')"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "b177f659",
   "metadata": {},
   "source": [
    "<img src='matplotlib histogram.png' />"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "c324dcf3",
   "metadata": {},
   "source": [
    "<h1>Using the Scripting layer</h1>"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "ec79e6fa",
   "metadata": {},
   "source": [
    "<h2>matplotlib in inline mode</h2>"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "id": "be16a32f",
   "metadata": {},
   "outputs": [
    {
     "ename": "AttributeError",
     "evalue": "'numpy.ndarray' object has no attribute 'mode'",
     "output_type": "error",
     "traceback": [
      "\u001b[1;31m---------------------------------------------------------------------------\u001b[0m",
      "\u001b[1;31mAttributeError\u001b[0m                            Traceback (most recent call last)",
      "\u001b[1;32m~\\AppData\\Local\\Temp\\ipykernel_112380\\2128348373.py\u001b[0m in \u001b[0;36m<module>\u001b[1;34m\u001b[0m\n\u001b[0;32m      1\u001b[0m \u001b[1;32mimport\u001b[0m \u001b[0mmatplotlib\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mpyplot\u001b[0m \u001b[1;32mas\u001b[0m \u001b[0mplt\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m----> 2\u001b[1;33m \u001b[0mmode\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mx\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mmode\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m      3\u001b[0m \u001b[0mplt\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mhist\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mx\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;36m100\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m      4\u001b[0m \u001b[0mplt\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mplot\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mmod\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;36m200\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;34m'o'\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m      5\u001b[0m \u001b[0mplt\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mplot\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;36m0\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;36m300\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;34m'o'\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
      "\u001b[1;31mAttributeError\u001b[0m: 'numpy.ndarray' object has no attribute 'mode'"
     ]
    }
   ],
   "source": [
    "import matplotlib.pyplot as plt\n",
    "plt.hist(x, 100)\n",
    "plt.plot(0, 300, 'o')\n",
    "plt.title(r'Normal distribution with $\\mu=0, \\sigma=1$')\n",
    "plt.savefig('matplotlib histogram.png')\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "da5df47a",
   "metadata": {},
   "source": [
    "<h2>matplotlib in notebook mode</h2>"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "id": "eb42ff44",
   "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_device_pixel_ratio', {\n",
       "                device_pixel_ratio: fig.ratio,\n",
       "            });\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",
       "    fig.rubberband_canvas.style.cursor = msg['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 https://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",
       " * https://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\", \"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",
       "    // 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,iVBORw0KGgoAAAANSUhEUgAAAoAAAAHgCAYAAAA10dzkAAAAAXNSR0IArs4c6QAAIABJREFUeF7tvQm4VdV5//8iDWhkJmCYRIJoFCrBxKASNRQBhyagtECaKBBFi9aIASTBaKDMkhj/bURAU6BQwDqAARHSSACVxyZBTIkoDVEGAVGZnBiC3v/zrl/PfS6Hc7n7nPfsae3Pfh6fRO9ee6/1eb/nrs9Ze7i1KioqKoQNAhCAAAQgAAEIQCAzBGohgJmpNQOFAAQgAAEIQAACjgACSBAgAAEIQAACEIBAxggggBkrOMOFAAQgAAEIQAACCCAZgAAEIAABCEAAAhkjgABmrOAMFwIQgAAEIAABCCCAZAACEIAABCAAAQhkjAACmLGCM1wIQAACEIAABCCAAJIBCEAAAhCAAAQgkDECCGDGCs5wIQABCEAAAhCAAAJIBiAAAQhAAAIQgEDGCCCAGSs4w4UABCAAAQhAAAIIIBmAAAQgAAEIQAACGSOAAGas4AwXAhCAAAQgAAEIIIBkAAIQgAAEIAABCGSMAAKYsYIzXAhAAAIQgAAEIIAAkgEIQAACEIAABCCQMQIIYMYKznAhAAEIQAACEIAAAkgGIAABCEAAAhCAQMYIIIAZKzjDhQAEIAABCEAAAgggGYAABCAAAQhAAAIZI4AAZqzgDBcCEIAABCAAAQgggGQAAhCAAAQgAAEIZIwAApixgjNcCEAAAhCAAAQggACSAQhAAAIQgAAEIJAxAghgxgrOcCEAAQhAAAIQgAACSAYgAAEIQAACEIBAxggggBkrOMOFAAQgAAEIQAACCCAZgAAEIAABCEAAAhkjgABmrOAMFwIQgAAEIAABCCCAZAACEIAABCAAAQhkjAACmLGCM1wIQAACEIAABCCAAJIBCEAAAhCAAAQgkDECCGDGCs5wIQABCEAAAhCAAAJIBiAAAQhAAAIQgEDGCCCAGSs4w4UABCAAAQhAAAIIIBmAAAQgAAEIQAACGSOAAGas4L4Pd86cOTJkyBCpW7eubN68Wdq2bXvckL/+9a/Le++9J3/84x9ThyI3tjfffFPOOuusovs/duxYGTdunFRUVFS2LeWY69atk1/96lcyfPhwadSoUeB+5J8r1593331XPve5zwU+zsl2PFnfShlrWTpV4kEK9be68YXBssRul9Tsww8/lB/96Efyn//5n7Jv3z754he/KD/4wQ9k4MCBJR0vrkYffPCBjB8/Xl555RXZsGGD+13z4x//WLQ+bBBIGgEEMGkVoT8mArlJUw/yne98R+bNm4cA/h+BQgKo8vXnP/9ZunTp4qQ5yPaTn/xERo0aJcWKaP65wpCWk/WtlLEG4RHWPoX6W934wmAZ1rgKHbdXr17yu9/9TqZMmSLnnHOOLFiwQB599FH5j//4D/mHf/iHKLtiOtfWrVvlS1/6knTu3NmNQ8eAAJqQ0jhEAghgiHA5dPQEcgJ41VVXuVWql19+2f0yzm3lXAH8+OOP5bOf/Wxkg7SuYBUSwFI6X6wAVscpDGkptm+ljD/ONj4K4PLly+Xaa6910vetb32rEq9K4auvvirbt2+X2rVrx4k98Llzq+u1atVyq3/NmjVDAAPTY8eoCSCAURPnfKESyEnSqlWrZMCAAXLhhRfKihUrTiqAL7zwgvsl/dvf/lY++eQT9w3+nnvucZNSbsvJyvr162XSpEny3HPPyamnniq7d+92l3f00uof/vAHd/lHxVMnrMGDB8v999/vVtjuvPNOefHFF6Vp06Zy2223yd1331157C1btsjEiRNF+7Fz505p3Lix67ee56//+q8r9ytGAJ955hk3htdee01atmwpt99+u+hltiCXgHXlSds+++yz8s4770iDBg3caoa21T7q/+Zvv/nNb0Tl+mScqrsErJKu3H7961+LTpzf+MY35Gc/+5mbPHObsly9erXoCkvVrarU5v5/dX2rjl8x9ddbByZMmCAqLVp/zYj2tWHDhtXmWiWmU6dO7vLm3//937v9NEdf+cpX5Pzzz3eSk9u++c1vugzoz6vjVWh8ykbrUkr/csdr166dXH755TJ37tzjTvE3f/M37nOxZs2aUD67Q4cOlUWLFsn+/fvlr/7qryrPsXDhQrf6p5+bSy+9tKhzHzt2zNVFrwD86U9/ksOHDx/Xvk2bNk4sw9wQwDDpcuxyEEAAy0GRYySGQG7S1MtJer+UipfKmk5iuuWvAOqk1rNnT7nggguclOll0OnTp8t//dd/iU5AKpG65eRC7ynU+5KuvPJK+eijj6RPnz6VPzv33HPdZeeLL77YtVf5+6d/+icnNip95513nlvlmD17tjz55JNy/fXXu2OvXbtWfvnLX8oll1zipEfvgdJJeOXKle4+Ij2ubkEFUMfbu3dvd7zvf//7bvLWvuzZs8dNejXdA6irpyplKqUqfgcOHHD/rrLSrVs3d6x//dd/laeeekpatGjh+qY/U1E8GafqhEaZ9u/f3zFVGbr33nvdef/7v/9bPvOZz7jjBxHAt95666R9K8Sv2PprLTQTl112mWzcuFF++MMfOkn5t3/7t5N+BlTCVWxnzpzp9ps6daoTtkOHDjnh05+rtOgXhH/8x390P8/v78nG98ADD7jjldq/vXv3uvsw/+Vf/kXuuOOOyrFoVvQLyXe/+13Rc+Rv+nPNV5CtqtxV3V9zqsfQL2BVt5w4K7NbbrklyCkq9xk0aJC7fPy9731PdCVRs6+fb/1Cc9ddd7kvVnqvcG4rxzjyO4gAFlUydo6BAAIYA3ROGR6BqgKoUqfS1aRJEze56OpSvgDq5PPGG2+4Vbp69eq5juVWAVV8VJi0XU5s7rvvvhNWwHI/++lPf+qEK7fpfXV6M7iK0nXXXef+s07yOtmrQKgEFtr0/J9++ql07NhR/vZv/7Zy4g0qgCqgO3bscGPSVSrd9OZ0fXBE5bImAaxfv77cfPPNbgWl0Hayy6wn41SdAOqEXFUuVJK//e1vy/z5893/6hZEAHW/k/WtEL9i66/yq/c/5jZdWVX508vcmpPqthtuuMF9IdGa6KZfOnTF7YknnpAHH3xQbrzxRvdzFWxdQdafF+pvTZeAS+2fftlQ8dc+KJPc9r//+79OKqvWouoYdeWxe/fugT7Q1d0zqrL/hS984biVej2grq7rZ0VXwlW0g265/MyaNUt0dTG35XguXbrUfa7KPY78/iGAQSvGfnERQADjIs95QyFQVQD1ElvuMpJeYtKVm6oCqCt4KjvDhg2Thx566Lj+6EQ6evRodwlVn0iseplXxbLqlvuZPnWsk1lu05UhPa/KQU7E9Gd6Oevo0aPy+9//3u2qUqjn00lWLwf/5S9/qTyGTsp6KVa3IAKoY9KVOF1x1FW6qptKlK4s1iSAPXr0cCt+I0aMcKtyX/7ylytX4vR4QQRQL4fnc6pOAJWDniO3KY/TTjtNdBVHb6LXLQwBLKX+r7/+euWKrPZLV6d0xe7tt9+WM844o9pM58auXzZ01VRX1XSFSv/R+0j1UuU///M/O9nRS6E6/lIEsNT+6WVtXUF8//333blzW+7zk/sc5A9Qv1ho7oNsmoc6deqcsKt+Ztq3b1+Z89wOOQGcPHmyeyI46HbFFVeI3sawadOm45roirJ+OdLPgAp31a0c48jvHwIYtGLsFxcBBDAu8pw3FAL5AqiyoyJ48OBBJ3O6spJ7DYxeUtN7gfT+M30FRdVNZUxXbfT+MF2VyUmeXkKqem+atqnuYQaVFl3h0Xvvqm75q5B6mUoFVIVTJy+Vg1NOOcWtwulrVnSVRbcgApgbk07oeh9f1U0nUb20WJMAKh9tv2TJEtm2bZtbGdUVTJXUz3/+84EEsBCn6gRQ+9yqVavj+qrn0ZWoxYsXu/8ehgCWUv/8V9YEqYn2P3cuXZXSlT8Ve+WsgqX527Vrl7v/Ti956yX86upd0wpgqf3r27eve6pbxb3qpl8CtM/6+dFM5m/luHRazkvAKs96GV1X4pVV1U1X3P/u7/7uuFtCcj8vxzgQwFB+pXPQEAkggCHC5dDRE8gXQO2B3oOn4vfzn/9cHn/88UoBzK2W6QpOdSuAuRWVkz2xahVAvUSt9xLqvYFVt9atW8vZZ59dlADmVrX00mSpK4BV+6CXwPX+RJXHr33ta+4yXZAVwELv9rOsAGqNnn76aXdZsOqm91hq7XJSW8wl4HLUP6gAap/1Uqo+YKSX4vXew5deeklyl1j1/+ttAboKl7vcWcoKYKkCqFnTz0h+BnVl9vTTT3f3qRbaynEJWO/vUxHOfwhEV8/1qeBiHgL5n//5H/fUf/7tGNp3/UKn+dUM5d+PWI5xIIDR/77njDYCCKCNH60TRqCQAGoX9UZwXd3QFT99IjD3Imi9HKsrH3ppLnfpS++/04laJ6T8ewALiY1VAHXFQp8OnTFjRiVNfYpX71PSFcFiVgD1ANZ7AAuVVFcAdSLWlT0VS1211Etseo9l1e1kolzsPYB6WVQfqtFN3w+nK5q6Upa71KqX0fU+Sb1snhPAk/WtkFBZ61+MAKqU65PAmkF9elhX/nTTh2D0NgO990/vVb3ooovcfy907OrGVx33IP3Ty9d6WVpXnB955JHKcqqk6mq1vvC7uvtBy3HpVG9xuOaaa9ztErmHrrQTV199tajQFfMaGN1XeeaPRe9t1BVWfdpfHzLK38oxDgQwYZMB3amRAAJYIyJ2SBOB6gRQn6bV1QwVBZWGnADmngLV166MHDnS3aOkTwHrZFzoKeAwBFDvdXvsscec5Oh9UvoKkGnTprknkvXeqGIFUJ9A1kuMKjd6CU8fKtFLvzrR68MhJ7sErJf69KZ+vX9RpUTvkdQnqnXS1KeW9Z613GrJrbfe6u7T08uWurql+5YigLmngHUFKvcUsI5bZSh3z5hKup5DL8frQxgq8frEqk74+rPcmE7Wt5M9BVxq/YMIVu7zow8D9evXz/2r5k6FRDd9wlZX3vTSv14Wzl1qLXTs6sanK166eljKCqA+FKGvn9FX2agg6aqzPrykK6v69KxmQVcl9VU2YW36BU3vBdWc6vn1s6cyWujhE33YpuoXo/w+6c9U+JTHV7/6VZcjvbdSeS9btqzgpexyjkuFVleXVSq1tvrlTp9y101FN8p3h5ZzXBzLPwIIoH81zfSIqhNAhaJPlOoTglUFUP977j1wepO4rv7pJSRdbar6pGCYl4D1aWOVT52I9X5BlRG98T13X2KxAqhj0mNpe72ErffT6UMh+sqRmt4DeOTIEbfioxOovnNPH0g588wz3atv9DUauclrzJgx7mZ6lUpllv8ewGIuAavwKl89hk7uyl2fjG3evPlxWdaJVc+rDx3oipXe56XnyR9TdX2rTtYs9S9GALXOutqrK826upx7xU3uqVUV7KpPhld37ELjy70HsBQBVOlTQVLpUrnWmuoqssqYfmZUZrRfVZ8OLvcvGc29fuaq/ik4lc78PwWn++kXDf3v2t9Cm65Sa4Y1L/oAlj5koveQ6iuhqnsVTTnHo5f49d7ZQluxfz2nnP3iWBDIJ4AAkgkIQAACGSagq1Iqffrkd9I3fQG3fkHQ2zmqviQ96f2mfxBIIgEEMIlVoU8QgAAEIiKg91TqPZ5V70GN6NRFn0ZXKPXF2bpqygYBCNgIIIA2frSGAAQgkFoCuYcmfvGLX7j71dggAIHsEEAAs1NrRgoBCEAAAhCAAAQcAQSQIEAAAhCAAAQgAIGMEUAAM1ZwhgsBCEAAAhCAAAQQQDIAAQhAAAIQgAAEMkYAAcxYwRkuBCAAAQhAAAIQQAANGdAX4OqfptIXk+oLbNkgAAEIQAACEEg+Af3rQfrXWlq2bBn6X4dJKg0E0FCZt956y/1dTzYIQAACEIAABNJHQP88ZuvWrdPX8TL0GAE0QNS/m9qoUSP391UbNGhgOBJNIQABCEAAAhCIisD777/vFnD0TzTq38HO4oYAGqquAdLgqAgigAaQNIUABCAAAQhESID5m/cAmuJGgEz4aAwBCEAAAhCIhQDzNwJoCh4BMuGjMQQgAAEIQCAWAszfCKApeATIhI/GEIAABCAAgVgIMH8jgKbgESATPhpDAAIQgAAEYiHA/I0AmoJHgEz4aAwBCEAAAhCIhQDzNwJoCh4BMuGjMQQgAAEIQCAWAszfCKApeATIhI/GEIAABCAAgVgIMH8jgKbgESATPhpDAAIQgAAEYiHA/I0AmoJHgEz4aAwBCEAAAhCIhQDzNwJoCh4BMuGjMQQgAAEIQCAWAszfCKApeATIhI/GEIAABCAAgVgIMH8jgKbgESATPhpDAAIQgAAEYiHA/I0AmoJHgEz4aAwBCEAAAhCIhQDzNwJoCh4BMuGjMQQgAAEIQCAWAszfCKApeATIhI/GEEgUgbN+8MwJ/dk65dpE9ZHOQAAC5SHA/I0AmpJEgEz4aAyBRBFAABNVDjoDgVAJMH8jgKaAESATPhpDIFEEEMBElYPOQCBUAszfCKApYATIhI/GEEgUAQQwUeWgMxAIlQDzNwJoChgBMuGjMQQSRQABTFQ56AwEQiXA/I0AmgJGgEz4aAyBRBFAABNVDjoDgVAJMH+nVAAffvhh0X+2bt3qAtKxY0e577775Oqrr3b/PnjwYJk7d+5x4enatau89NJLlf/tyJEjMnLkSFm4cKEcOnRIevToIdOnT5fWrVsHDh0BCoyKHSGQeAIIYOJLRAchUDYCzN8pFcClS5dK7dq15eyzz3ZhUNmbNm2abNiwwcmgCuCePXtk9uzZlWGpU6eONGnSpPLfhw0bJnqcOXPmSNOmTWXEiBGyb98+Wb9+vTt2kI0ABaHEPhBIBwEEMB11opcQKAcB5u+UCmCh4qvcqQTedNNNTgAPHDggS5YsKZiTgwcPSrNmzWTevHkyYMAAt8+uXbukTZs2snz5cundu3egfBGgQJjYCQKpIIAApqJMdBICZSHA/O2BAH7yySfy+OOPy6BBg9wK4Pnnn+8EUOVPV/0aNWokV1xxhUycOFGaN2/ugrNq1Sp3yVdX/Bo3blwZps6dO0vfvn1l3LhxgQJGgAJhYicIpIIAApiKMtFJCJSFAPN3igVw48aNcskll8jhw4elXr16smDBArnmmmtcMB577DH339q2bStvvvmm3HvvvXLs2DF3ebdu3bpu3yFDhojeB1h169Wrl7Rr105mzpxZMGC6f9U2GiBdNdQVxQYNGpQllBwEAhCIhwACGA93zgqBOAgggCkWwKNHj8r27dvdpd4nn3xSHn30UVmzZo1bAczfdu/e7WRw0aJFcv3111crgD179pT27dvLjBkzCuZx7NixBVcHEcA4Pr6cEwLlJYAAlpcnR4NAkgkggCkWwPxgXXnllU7eqlu969Chg9x8880yevToki8BswKY5I8zfYOAjQACaONHawikiQAC6JEA6j19ejlWn+rN3/bu3SutWrWSWbNmyY033ugu2epDIPPnz5f+/fu73XWVUF8Bw0MgafoI01cIlI8AAlg+lhwJAkkngACmVADHjBnj3vmnwvfBBx+4S7tTpkyRFStWuPsC9VJtv379pEWLFu5dgbq/Xi5+7bXXpH79+i6X+hqYZcuWOWHUJ4j1nYAqirwGJukfW/oHgXAIIIDhcOWoEEgiAQQwpQKor3p57rnn3Kpdw4YN5YILLnCXdvUePn2psz7Jq08E6/2BKoHdu3eX8ePHO2HMbfrwyKhRo9z9gFVfBF11n5pCS4BqIsTPIZAeAghgempFTyFgJcD8nVIBtBa+XO0JULlIchwIxE8AAYy/BvQAAlERYP5GAE1ZI0AmfDSGQKIIIICJKgedgUCoBJi/EUBTwAiQCR+NIZAoAghgospBZyAQKgHmbwTQFDACZMJHYwgkigACmKhy0BkIhEqA+RsBNAWMAJnw0RgCiSKAACaqHHQGAqESYP5GAE0BI0AmfDSGQKIIIICJKgedgUCoBJi/EUBTwAiQCR+NIZAoAghgospBZyAQKgHmbwTQFDACZMJHYwgkigACmKhy0BkIhEqA+RsBNAWMAJnw0RgCiSKAACaqHHQGAqESYP5GAE0BI0AmfDSGQKIIIICJKgedgUCoBJi/EUBTwAiQCR+NIZAoAghgospBZyAQKgHmbwTQFDACZMJHYwgkigACmKhy0BkIhEqA+RsBNAWMAJnw0RgCiSKAACaqHHQGAqESYP5GAE0BI0AmfDSGQKIIIICJKgedgUCoBJi/EUBTwAiQCR+NIZAoAghgospBZyAQKgHmbwTQFDACZMJHYwiERiBf5rZOubbGcyGANSJiBwh4Q4D5GwE0hZkAmfDRGAKhEUAAQ0PLgSHgBQHmbwTQFGQCZMJHYwicQKAUcSuEsZTjsAJIICGQHQLM3wigKe0EyISPxhBAAMkABCAQCwHmbwTQFDwCZMJHYwgggGQAAhCIhQDzNwJoCh4BMuGjMQQQQDIAAQjEQoD5GwE0BY8AmfDRGAI1CmAhRKU80VtKGz13kHaUEQIQSB8B5m8E0JRaAmTCR2MIIIBkAAIQiIUA8zcCaAoeATLhozEEEEAyAAEIxEKA+RsBNAWPAJnw0RgCCCAZgAAEYiHA/I0AmoJHgEz4aAwBBJAMQAACsRBg/kYATcEjQCZ8NIYAAkgGIACBWAgwfyOApuARIBM+GkMAASQDEIBALASYvxFAU/AIkAkfjSGAAJIBCEAgFgLM3wigKXgEyISPxhBAAMkABCAQCwHmbwTQFDwCZMJHYwgggGQAAhCIhQDzNwJoCh4BMuGjMQQQQDIAAQjEQoD5GwE0BY8AmfDRGAIIIBmAAARiIcD8jQCagkeATPhoDAEEkAxAAAKxEGD+RgBNwSNAJnw0hgACSAYgAIFYCDB/I4Cm4BEgEz4aQyAyASyEeuuUa4/7z2f94JkTdsvfh5JBAAJ+EGD+RgBNSSZAJnw0hgACSAYgAIFYCDB/I4Cm4BEgEz4aQwABJAMQgEAsBJi/EUBT8AiQCR+NIYAAkgEIQCAWAszfKRXAhx9+WPSfrVu3uuB07NhR7rvvPrn66qvdv1dUVMi4ceNk1qxZsn//funatas89NBDbr/cduTIERk5cqQsXLhQDh06JD169JDp06dL69atA4eRAAVGxY4QCESg0H14+Q2rvS/v009Etq0T+XCPDFz4pvz20y/Kp3JKteflHsBAJWEnCHhJgPk7pQK4dOlSqV27tpx99tkumHPnzpVp06bJhg0bnORNnTpVJk6cKHPmzJFzzjlHJkyYIGvXrpXNmzdL/fr1XZthw4aJHkf3adq0qYwYMUL27dsn69evd8cOshGgIJTYBwLBCZQsgJt+KbJitMj7uypPtquiiYz7y42y8tOvFuwAAhi8LuwJAd8IMH+nVAALBbFJkyZOAr/73e9Ky5YtZfjw4TJ69Gi3q672nXHGGU4Mb731Vjl48KA0a9ZM5s2bJwMGDHD77Nq1S9q0aSPLly+X3r17B8o6AQqEiZ0gEJhASQKo8vefN+ra/3Hn+fT//nXYX4YXlEAEMHBZ2BEC3hFg/vZAAD/55BN5/PHHZdCgQW4F8NRTT5X27dvLyy+/LF26dKkMbZ8+faRRo0ZutXDVqlXukq+u+DVu3Lhyn86dO0vfvn3d5eMgGwEKQol9IBCcQBABrHq0U+RTeaHu9+Tzsk9OqXXieVQC35am8rUj/98Jl4MRwOB1YU8I+EaA+TvFArhx40a55JJL5PDhw1KvXj1ZsGCBXHPNNbJu3Trp1q2b7Ny5060E5rZbbrlFtm3bJitXrnT7DhkyxK0MVt169eol7dq1k5kzZxbMuu5ftY0GSFcNdUWxQYMGvn0+GA8EIidQrABefMomWVRnQo39HHj0R/LSp+cftx8CWCM2doCAtwQQwBQL4NGjR2X79u1y4MABefLJJ+XRRx+VNWvWuH9XAdRLui1atKgM79ChQ2XHjh2yYsWKagWwZ8+ebvVwxowZBUM/duzYgquDCKC3vyMYWMQEihXAb56yTv6lzs9r7OX3jv6T/PLTSxHAGkmxAwSyQQABTLEA5kf0yiuvdPKm9/2FdQmYFcBs/GJglPERKFYAk7ACyF8QiS8vnBkCpRJAAD0SQL2nTy/Hzp492136veuuu+Tuu+922dDVwubNm5/wEMj8+fOlf//+bp/du3e7V8DwEEipHyfaQcBOoFgBTMI9gAigve4cAQJRE0AAUyqAY8aMce/8U+H74IMPZNGiRTJlyhR3eVcv4+rTvpMnT3Yy2KFDB5k0aZKsXr36hNfALFu2zL0GRp8g1ncC7t27l9fARP0p5HyZJVCs7FUHqvcpv5WHP/Og+3HVB0GiegoYAcxshBl4igkggCkVwJtuukmee+45t2rXsGFDueCCC9ylX5U/3XIvgtaHOaq+CLpTp06VcdWHR0aNGuXuB6z6ImiVyqAbAQpKiv0gIFIu4SvEUiXwx5/5d2lZa1/lj3dVNJVxf7mh2vcABqlJtS+drtIYAQxCkn0gkCwCzN8pFcCkxIgAJaUS9CMNBMIUQLf6J5/KV095XZrLAXlHGtX4l0CCMEMAg1BiHwikjwDzNwJoSi0BMuGjccYIhC2AYeBEAMOgyjEhED8B5m8E0JRCAmTCR+OMEUijABYqEe8PzFhwGa6XBJi/EUBTsAmQCR+NM0YAAcxYwRkuBBJMgPkbATTFkwCZ8NE4YwQQwIwVnOFCIMEEmL8RQFM8CZAJH40zRgABzFjBGS4EEkyA+RsBNMWTAJnw0ThjBBDAjBWc4UIgwQSYvxFAUzwJkAkfjTNGAAHMWMEZLgQSTID5GwE0xZMAmfDROGMEEMCMFZzhQiDBBJi/EUBTPAmQCR+NM0YAAcxYwRkuBBJMgPkbATTFkwCZ8NE4YwQQwIwVnOFCIMEEmL8RQFM8CZAJH40zRgABzFjBGS4EEkyA+RsBNMWTAJnw0ThjBBDAjBWc4UIgwQSYvxFAUzwJkAkfjT0m4IvsFSoRfwrO4+AytMwQYP5GAE22wYfNAAAgAElEQVRhJ0AmfDT2mAACeK3H1WVoEEg/AeZvBNCUYgJkwkdjjwkggAigx/FmaB4QYP5GAE0xJkAmfDT2mAACiAB6HG+G5gEB5m8E0BRjAmTCR2OPCSCACKDH8WZoHhBg/kYATTEmQCZ8NPaYAAKIAHocb4bmAQHmbwTQFGMCZMJHY48JIIAIoMfxZmgeEGD+RgBNMSZAJnw09piAzwIYpGz5r4oJ0oZ9IACB6AgwfyOAprQRIBM+GntEIOvCl19KBNCjcDMULwkwfyOApmATIBM+GntEAAE8vpgIoEfhZiheEmD+RgBNwSZAJnw09ogAAogAehRnhpIBAszfCKAp5gTIhI/GHhFAABFAj+LMUDJAgPkbATTFnACZ8NHYIwIIIALoUZwZSgYIMH8jgKaYEyATPhp7RAABRAA9ijNDyQAB5m8E0BRzAmTCR2OPCCCACKBHcWYoGSDA/I0AmmJOgEz4aOwRAQQQAfQozgwlAwSYvxFAU8wJkAkfjT0igAAigB7FmaFkgADzNwJoijkBMuGjsUcEEEAE0KM4M5QMEGD+RgBNMSdAJnw09ogAAogAehRnhpIBAszfCKAp5gTIhI/GHhFAABFAj+LMUDJAgPkbATTFnACZ8NHYIwIIIALoUZwZSgYIMH8jgKaYEyATPhp7RAABRAA9ijNDyQAB5m8E0BRzAmTCR2OPCCCACKBHcWYoGSDA/I0AmmJOgEz4aOwRAQTw5MXcOuVaj6rNUCCQfgLM3wigKcUEyISPxh4RQAARQI/izFAyQID5GwE0xZwAmfDR2CMCCCAC6FGcGUoGCDB/p1QAJ0+eLE899ZS8/vrrctppp8mll14qU6dOlXPPPbcytoMHD5a5c+ceF+OuXbvKSy+9VPnfjhw5IiNHjpSFCxfKoUOHpEePHjJ9+nRp3bp1oPgToECY2CkDBBBABDADMWeIHhFg/k6pAF511VUycOBAueiii+TYsWNyzz33yMaNG2XTpk1y+umnu4iqAO7Zs0dmz55dGdk6depIkyZNKv992LBhsnTpUpkzZ440bdpURowYIfv27ZP169dL7dq1a4w6AaoRETtkhAACiABmJOoM0xMCzN8pFcD8/L377rvSvHlzWbNmjVx++eWVAnjgwAFZsmRJwbgePHhQmjVrJvPmzZMBAwa4fXbt2iVt2rSR5cuXS+/evWuMOQGqERE7ZIQAAogAZiTqDNMTAszfngjgli1bpEOHDm4VsFOnTpUCqPKnq36NGjWSK664QiZOnOhEUbdVq1a5S7664te4cePKSHfu3Fn69u0r48aNOyHmeslY/8ltGiAVRpXJBg0aePKxYBgQKJ4AAogAFp8aWkAgPgIIoAcCWFFRIX369JH9+/fL888/X5mmxx57TOrVqydt27aVN998U+699153uVgv79atW1cWLFggQ4YMOU7otHGvXr2kXbt2MnPmzBOSOXbs2IJiiADG9yHmzMkggAAigMlIIr2AQDACCKAHAnj77bfLM888Iy+88MJJH97YvXu3k8FFixbJ9ddfX60A9uzZU9q3by8zZsxgBTDY54i9ICAIIALIxwACaSKAAKZcAO+44w53j9/atWvdql1Nm14mvvnmm2X06NElXQLOPz4Bqok4P88KAQQQAcxK1hmnHwSYv1MqgHrZV+Vv8eLFsnr1anf/X03b3r17pVWrVjJr1iy58cYb3X17+hDI/PnzpX///q65rhLqK2B4CKQmmvwcAscTQACLTwR/HaR4ZrSAQLkIIIApFcDbbrvNXcJ9+umnj3v3X8OGDd17AT/88EPR+/X69esnLVq0kK1bt8qYMWNk+/bt8tprr0n9+vVdhvQ1MMuWLXOvgdHXw+g7AVUUeQ1MuT5iHCcrBBDA4iuNABbPjBYQKBcBBDClAlirVq2CGdB3/un7//Slzvok74YNG0RfBaMS2L17dxk/frx7aje3HT58WEaNGuVksuqLoKvuc7KwEaByfRQ5TtoJIIDFVxABLJ4ZLSBQLgLM3ykVwHIFwHocAmQlSHtfCCCAxVcSASyeGS0gUC4CzN8IoClLBMiEj8YeEUAAiy8mAlg8M1pAoFwEmL8RQFOWCJAJH409IoAA2ouJENoZcgQIBCXA/I0ABs1Kwf0IkAkfjVNCALmLplAIYDScOQsElADzNwJo+iQQIBM+GqeEAAIYTaEQwGg4cxYIIID/LwO1KvSlemwlEUAAS8JGo5QRQACjKRgCGA1nzgIBBBABNH8KEEAzQg6QMALIXnwFQQDjY8+Zs0eA+ZsVQFPqCZAJH40TSAABjK8oCGB87Dlz9ggwfyOAptQTIBM+GieQAAIYX1EQwPjYc+bsEWD+RgBNqSdAJnw0TiABBDC+oiCA8bHnzNkjwPyNAJpST4BM+GicQAIIYHxFQQDjY8+Zs0eA+RsBNKWeAJnw0TiBBBDA+IqCAMbHnjNnjwDzNwJoSj0BMuGjcQIJIIDxFQUBjI89Z84eAeZvBNCUegJkwkfjBBJAAOMrCgIYH3vOnD0CzN8IoCn1BMiEj8YJJIAAxlcUBDA+9pw5ewSYvxFAU+oJkAkfjRNIAAGMrygIYHzsOXP2CDB/I4Cm1BMgEz4aJ5AAAhhfURDA+Nhz5uwRYP5GAE2pJ0AmfDROIAEEML6iIIDxsefM2SPA/I0AmlJPgEz4aJxAAghgfEVBAONjz5mzR4D5GwE0pZ4AmfDROGICheQuXzoQwIiLUuV0CGB87Dlz9ggwfyOAptQTIBM+GkdMAAGMGHiRp0MAiwTG7hAwEGD+RgAN8REhQCZ8NI6YAAIYMfAiT4cAFgmM3SFgIMD8jQAa4oMAmuDROHICCGDkyIs6IQJYFC52hoCJAAKIABIgEwEap4kAApjsaiGAya4PvfOLAAKIAJoSTYBM+GgcMQEEMGLgRZ4OASwSGLtDwECA+RsBNMSHS8AmeDSOnAACGDnyok6IABaFi50hYCKAACKABMhEgMZpIoAAJrtaCGCy60Pv/CKAACKApkQTIBM+GkdMAAGMGHiRp0MAiwTG7hAwEGD+RgAN8eESsAkejSMngABGjryoEyKAReFiZwiYCCCACCABMhGgcZoIIIDJrhYCmOz60Du/CCCACKAp0QTIhI/GERNAACMGXuTpEMAigbE7BAwEmL8RQEN8uARsgkfjyAkggJEjN58QKTQj5AAQKEgAAUQATR8NAmTCR+OICSCAEQMvw+kQwDJA5BAQKECA+RsBNH0wCJAJH40jJoAARgy8DKdDAMsAkUNAAAEsmIFaFRUVFaSjNAIIYGncaBUNgULCl3/mfMEI0iaa3nMWJYAAkgMIhEOA+ZsVQFOyCJAJH41DJhBE5hDAkItgPDwCaARIcwhUQ4D5GwE0fTgIkAkfjUMmgACGDDiGwyOEMUDnlF4SYP5GAE3BJkAmfDQOmQACGDLgGA6PAMYAnVN6SYD5GwE0BZsAmfDROGQCQQQw5C5w+DITQADLDJTDZZYA83dKBXDy5Mny1FNPyeuvvy6nnXaaXHrppTJ16lQ599xzK8Osz7aMGzdOZs2aJfv375euXbvKQw89JB07dqzc58iRIzJy5EhZuHChHDp0SHr06CHTp0+X1q1bB/pQEKBAmNgpJgIIYEzgQzwtAhgiXA6dKQLM3ykVwKuuukoGDhwoF110kRw7dkzuuece2bhxo2zatElOP/10F2IVwokTJ8qcOXPknHPOkQkTJsjatWtl8+bNUr9+fbfPsGHDZOnSpW6fpk2byogRI2Tfvn2yfv16qV27do0fBgJUIyJ2iJEAAhgj/JBOjQCGBJbDZo4A83dKBTA/qe+++640b95c1qxZI5dffrno6l/Lli1l+PDhMnr0aLe7rvadccYZTgxvvfVWOXjwoDRr1kzmzZsnAwYMcPvs2rVL2rRpI8uXL5fevXvX+IEgQDUiYocYCSCAMcIP6dQIYEhgOWzmCDB/eyKAW7ZskQ4dOrhVwE6dOskbb7wh7du3l5dfflm6dOlSGew+ffpIo0aNZO7cubJq1Sp3yVdX/Bo3bly5T+fOnaVv377u8nH+phKp/+Q2DZAKo8pkgwYNMvcBYsDJJoAAJrs+pfQOASyFGm0gcCIBBNADAdTVPhU7vc/v+eefd1Vet26ddOvWTXbu3OlWAnPbLbfcItu2bZOVK1fKggULZMiQIccJne7Xq1cvadeuncycOfOExIwdO7agGCKA/HpJIgEEMIlVsfUJAbTxozUEqi7gNGzYMNMLOKn/SyC33367PPPMM/LCCy9UPryRE0C9pNuiRYvKxA8dOlR27NghK1asqFYAe/bs6VYPZ8yYwQogvytSTQABTHX5CnYeAfSvpowoHgKsAKZ8BfCOO+6QJUuWuIc7dNUut4V1CTg/pgQong8uZw1GAAEMxilNeyGAaaoWfU0yAebvlAqgXvZV+Vu8eLGsXr3a3f9Xdcs9BHLXXXfJ3Xff7X509OhR96BI/kMg8+fPl/79+7t9du/e7VYReQgkyR9b+haUAAIYlFR69kMA01MreppsAghgSgXwtttuc5dwn3766ePe/afX8/W9gLqp6On7AmfPnu0EcdKkSU4W818Ds2zZMvcamCZNmrh3Au7du5fXwCT7c0vvAhJAAAOCStFuCGCKikVXE00AAUypANaqVatgsFT2Bg8e7H6WexG0PsxR9UXQ+pRwbjt8+LCMGjXKyWTVF0Hrk71BNgIUhBL7xEUAAYyLfHjnRQDDY8uRs0WA+TulApiUmBKgpFSCfhQigAD6lwsE0L+aMqJ4CDB/I4Cm5BEgEz4ah0wAAQwZcAyHRwBjgM4pvSTA/I0AmoJNgEz4aBwyAQQwZMAxHB4BjAE6p/SSAPM3AmgKNgEy4aNxyAQQwJABx3B4BDAG6JzSSwLM3wigKdgEyISPxiETQABDBhzD4RHAGKBzSi8JMH8jgKZgEyATPhqHTAABDBlwDIdHAGOAzim9JMD8jQCagk2ATPhoHDIBBDBkwAk5PFKYkELQjVQRYP5GAE2BJUAmfDQOmQACGDLghBweAUxIIehGqggwfyOApsASIBM+GodMAAEMGXBCDo8AJqQQdCNVBJi/EUBTYAmQCR+NQyaAAIYMOKGHRwgTWhi6lSgCzN8IoCmQBMiEj8YhE0AAQwac0MMjgAktDN1KFAHmbwTQFEgCZMJH45AJIIAhA07o4RHAhBaGbiWKAPM3AmgKJAEy4aNxyAQQwJABJ/TwCGBCC0O3EkWA+RsBNAWSAJnw0ThkAghgyIATengEMKGFoVuJIsD8jQCaAkmATPhoHDIBBDBkwAk9PAKY0MLQrUQRYP5GAE2BJEAmfDQOmQACGDLghB4eAUxoYehWoggwfyOApkASIBM+GodMAAEMGXBCD48AJrQwdCtRBJi/EUBTIAmQCR+NQyaAAIYMOKGHRwATWhi6lSgCzN8IoCmQBMiEj8ZlJIDslRFmyg+FAKa8gHQ/EgLM3wigKWgEyISPxmUkgACWEWbKD4UApryAdD8SAszfCKApaATIhI/GZSSAAJYRZsoPhQCmvIB0PxICzN8IoCloBMiEj8ZlJIAAlhFmyg+FAKa8gHQ/EgLM3wigKWgEyISPxmUkgACWEWbKD4UApryAdD8SAszfCKApaATIhI/GZSSAAJYRZsoPhQCmvIB0PxICzN8IoCloBMiEj8ZlJIAAlhFmyg+FAKa8gHQ/EgLM3wigKWgEyISPxmUkgACWEWbKD4UApryAdD8SAszfCKApaATIhI/GZSSAAJYRZsoPhQCmvIB0PxICzN8IoCloBMiEj8ZlJIAAlhFmyg+FAKa8gHQ/EgLM3wigKWgEyISPxmUkgACWEWbKD4UApryAdD8SAszfCKApaATIhI/GZSSAAJYRZsoPhQCmvIB0PxICzN8IoCloBMiEj8ZlJIAAlhFmyg+FAKa8gHQ/EgLM3wigKWgEyISPxmUkgACWEaaHh0IKPSwqQzIRYP5GAAmQiQCNk0IAAUxKJZLZDwQwmXWhV/ERQAARQFP6CJAJH43LSAABLCNMDw+FAHpYVIZkIsD8jQASIBMBGieFAAKYlEoksx8IYDLrQq/iI4AAIoCm9BEgEz4al5EAAlhGmB4eCgH0sKgMyUSA+RsBJEAmAjROCgEEMCmVSGY/EMBk1oVexUcAAUQATekjQCZ8NC4jAQSwjDA9PBQC6GFRGZKJAPN3SgVw7dq1Mm3aNFm/fr3s3r1bFi9eLH379q0Mw+DBg2Xu3LnHhaNr167y0ksvVf63I0eOyMiRI2XhwoVy6NAh6dGjh0yfPl1at24dOFQEKDAqdgyZAAIYMuCUHx4BTHkB6X7ZCTB/p1QAn332WXnxxRflwgsvlH79+hUUwD179sjs2bMrQ1OnTh1p0qRJ5b8PGzZMli5dKnPmzJGmTZvKiBEjZN++fU4qa9euHShsBCgQJnaKgAACGAHkFJ8CAUxx8eh6KASYv1MqgFXTUKtWrYICeODAAVmyZEnB4Bw8eFCaNWsm8+bNkwEDBrh9du3aJW3atJHly5dL7969AwWOAAXCxE4REEAAI4Cc4lMggCkuHl0PhQDzt8cCqPKnq36NGjWSK664QiZOnCjNmzd3QVq1apW75Ksrfo0bN64MV+fOnd2l5HHjxgUKHAEKhImdIiCAAEYAOcWnQABTXDy6HgoB5m9PBfCxxx6TevXqSdu2beXNN9+Ue++9V44dO+Yu79atW1cWLFggQ4YMEb0PsOrWq1cvadeuncycObNg4HT/qm00QLpqqCuKDRo0CCWkHBQCQQgggEEoZXcfBDC7tWfkhQkggJ4KYH659UERlcFFixbJ9ddfX60A9uzZU9q3by8zZswomJixY8cWXB1EAPkVEzcBBDDuCiT7/AhgsutD76IngABmRAA1Wh06dJCbb75ZRo8eXfIlYFYAo/+QcsZgBBDAYJyyuhcCmNXKM+7qCCCAGRHAvXv3SqtWrWTWrFly4403uku2+hDI/PnzpX///i4fukqor4DhIRB+YaSRAAKYxqpF12cEMDrWnCkdBBDAlArghx9+KFu2bHEp69KlizzwwAPSvXt395oX/Ucv1errYVq0aCFbt26VMWPGyPbt2+W1116T+vXru3b6Gphly5a518BoG30noIoir4FJx4eXXh5PAAEkEScjgACSDwgcTwABTKkArl692glf/jZo0CB5+OGH3ZO8GzZsEH0VjEqg7jt+/Hj3wEZuO3z4sIwaNcrdD1j1RdBV96npA0OAaiLEz6MigABGRTqd50EA01k3eh0eAebvlApgeJEo7sgEqDhe7B0eAQQwPLY+HBkB9KGKjKGcBJi/EUBTngiQCR+NAxIoJHf5EzoCGBBmRndDADNaeIZdLQHmbwTQ9PEgQCZ8NA5IALkLCIrdqiWAABIOCBxPgPkbATR9JgiQCR+NAxJAAAOCYjcEkAxAICAB5m8EMGBUCu9GgEz4aByQAAIYEBS7IYBkAAIBCTB/I4ABo4IAmkDR2EQAATTho7GIcAmYGECAS8D5GahVUVFRQTBKI8A3iNK40ao4AghgcbzY+0QCCCCpgAACiACW8VOAAJYRJoeqlgACSDisBBBAK0Ha+0aA+ZtLwKZMEyATPhoHJIAABgTFbtUSQAAJBwRYAWQFsIyfAgSwjDA5FCuAZCA0AghgaGg5cEoJMH+zAmiKLgEy4aNxQAKsAAYExW6sAJIBCAQkwPyNAAaMSuHdCJAJH40DEkAAA4JiNwSQDEAgIAHmbwQwYFQQQBMoGpsIIIAmfDTmNTBkAAInEEAAEUDTx4IAmfDROCABBDAgKHZjBZAMQCAgAeZvBDBgVFgBNIGisYkAAmjCR2NWAMkABFgBLJABXgRt+GDwDcIAj6aBCSCAgVGxYzUEeAqYaEDgeALM36wAmj4TBMiEj8YBCSCAAUGxG5eAyQAEAhJg/kYAA0aFS8AmUDQ2EUAATfhozCVgMgABLgFzCbi8nwK+QZSXJ0f7fwQQPpJQbgJcAi43UY6XdgLM36wAmjJMgEz4aFwNAQSQaJSbAAJYbqIcL+0EmL8RQFOGCZAJH40RQDIQEQEEMCLQnCY1BJi/EUBTWAmQCR+NEUAyEBMBhDAm8Jw2MQSYvxFAUxgJkAkfjRFAMhATAQQwJvCcNjEEmL8RQFMYCZAJH40RQDIQEwEEMCbwnDYxBJi/EUBTGAmQCR+NEUAyEBMBBDAm8Jw2MQSYvxFAUxgJkAkfjRFAMhATAQQwJvCcNjEEmL8RQFMYCZAJH40RQDIQEwEEMCbwnDYxBJi/EUBTGAmQCR+NEUAykCACSGGCikFXQifA/I0AmkJGgEz4aIwAkoEEEUAAE1QMuhI6AeZvBNAUMgJkwkdjBJAMJIgAApigYtCV0AkwfyOAppARIBM+GiOAZCBBBBDABBWDroROgPkbATSFjACZ8NEYASQDCSKAACaoGHQldALM3wigKWQEyISPxgggGUgQAQQwQcWgK6ETYP5GAE0hI0AmfDRGAMlAgggggAkqBl0JnQDzNwJoChkBMuGjMQJIBhJEAAFMUDHoSugEmL8RQFPICJAJH40RQDKQYAIIYYKLQ9fMBJi/EUBTiAiQCR+NReSsHzwDBwgkkgACmMiy0KkyEWD+RgBNUSJAJnw0RgDJQIIJIIAJLg5dMxNg/kYATSEiQCZ8NEYAyUCCCSCACS4OXTMTYP5OqQCuXbtWpk2bJuvXr5fdu3fL4sWLpW/fvpWBqKiokHHjxsmsWbNk//790rVrV3nooYekY8eOlfscOXJERo4cKQsXLpRDhw5Jjx49ZPr06dK6devAwSJAgVGxYzUEuARMNJJKAAFMamXoVzkIMH+nVACfffZZefHFF+XCCy+Ufv36nSCAU6dOlYkTJ8qcOXPknHPOkQkTJohK4+bNm6V+/fouO8OGDZOlS5e6fZo2bSojRoyQffv2OamsXbt2oHwRoECY2OkkBBBA4pFUAghgUitDv8pBgPk7pQJYtfi1atU6TgB19a9ly5YyfPhwGT16tNtVV/vOOOMMUTG89dZb5eDBg9KsWTOZN2+eDBgwwO2za9cuadOmjSxfvlx69+4dKF8EKBAmdkIAyUAKCSCAKSwaXQ5MgPnbQwF84403pH379vLyyy9Lly5dKsPQp08fadSokcydO1dWrVrlLvnqil/jxo0r9+ncubO7lKyXj4NsBCgIJfY5GQFWAMlHUgkggEmtDP0qBwHmbw8FcN26ddKtWzfZuXOnWwnMbbfccots27ZNVq5cKQsWLJAhQ4a4lcGqW69evaRdu3Yyc+bMgvnS/au20QDpqqGuKDZo0KAcmeQYGSOAAGas4CkfLlKY8gLS/UoCCKDHAqiXdFu0aFFZ7KFDh8qOHTtkxYoV1Qpgz5493erhjBkzCn5Mxo4dW3B1EAHkt0qpBBDAUsnRLg4CCGAc1DlnGAQQQA8FMMxLwKwAhvExzPYxEcBs1z9to0cA01Yx+lsdAQTQQwHMPQRy1113yd133+1qf/ToUWnevPkJD4HMnz9f+vfv7/bR18noK2B4CIRfGOUiUEju8idQBLBctDlOFAQQwCgoc44oCCCAKRXADz/8ULZs2eIyog96PPDAA9K9e3dp0qSJnHnmmU70Jk+eLLNnz5YOHTrIpEmTZPXq1Se8BmbZsmXuNTDaTt8JuHfvXl4DE8UnLyPnQAAzUugMDRMBzFCxPR8qAphSAVSZU+HL3wYNGuSELvciaH2Yo+qLoDt16lTZ5PDhwzJq1Ch3P2DVF0HrQx1BNwIUlFQ292N1L5t193nUCKDP1c3W2Ji/UyqASYkpAUpKJZLZDwQwmXWhV6UTQABLZ0fLZBFg/kYATYkkQCZ83jdGAL0vceYGiABmruTeDpj5GwE0hZsAmfB53xgB9L7EmRsgApi5kns7YOZvBNAUbgJkwud9YwTQ+xJnboAIYOZK7u2Amb8RQFO4CZAJn/eNEUDvS5y5ASKAmSu5twNm/kYATeEmQCZ83jdGAL0vceYGiABmruTeDpj5GwE0hZsAmfB53xgB9L7EmRsgApi5kns7YOZvBNAUbgJkwud9YwTQ+xJnboBB/pINkpi5WKRywMzfCKApuATIhM/7xgig9yXO3AARwMyV3NsBM38jgKZwEyATPu8bI4DelzhzA0QAM1dybwfM/I0AmsJNgEz4vG+MAHpfYgZYgACXgIlFGggwfyOAppwSIBM+7xsjgN6XmAEigGQgpQSYvxFAU3QJkAmf940RQO9LzAARQDKQUgLM3wigKboEyITP+8YIoPclZoAIIBlIKQHmbwTQFF0CZMLnfWME0PsSM0AEkAyklADzNwJoii4BMuHzvjEC6H2JGSACSAZSSoD5GwE0RZcAmfB53xgB9L7EDBABJAMpJcD8jQCaokuATPi8b4wAel9iBogAkoGUEmD+RgBN0SVAJnzeN0YAvS8xA0QAyUBKCTB/I4Cm6BIgEz7vGyOA3peYASKAZCClBJi/EUBTdAmQCZ/3jRFA70vMABFAMpBSAszfCKApugTIhM/7xgig9yVmgAggGUgpAeZvBNAUXQJkwud9YwTQ+xIzQASQDKSUAPM3AmiKLgEy4fO+MQLofYkZIAJIBlJKgPkbATRFlwCZ8HnfGAH0vsQMEAEkAyklwPyNAJqiS4BM+LxrjPB5V1IGVAKBrVOuLaEVTSAQLQHmbwTQlDgCZMLnXWME0LuSMqASCCCAJUCjSeQEmL8RQFPoCJAJn3eNEUDvSsqASiCAAJYAjSaRE2D+RgBNoSNAJnzeNUYAvSspAyqBAAJYAjSaRE6A+RsBNIWOAJnwedcYAfSupAyoBAIIYAnQaBI5AeZvBNAUOgJkwuddYwTQu5IyoBIIIIAlQKNJ5ASYvxFAU+gIkAmfd40RQO9KyoBKIIAAlgCNJpETYP5GAE2hI0AmfN41RgC9KykDKoEAAlgCNJpEToD5GwE0hY4AmfB517tNjjMAABPkSURBVBgB9K6kDKgEAghgCdBoEjkB5m8E0BQ6AmTC511jBNC7kjKgEgjkC2ChzwWSWAJYmpSVAPM3AmgKFAEy4UtV4/xJrNAEhgCmqqR0NkYCCGCM8Dm1I8D8jQCaPgoEyIQvVY0RwFSVi84mnAACmPACZaB7zN8IoCnmBMiEL1WNEcBUlYvOJpwAApjwAmWge8zfCKAp5gTIhC9Vjbm8m6py0dmEE0AAE16gDHSP+dtjARw7dqyMGzfuuBifccYZ8vbbb7v/VlFR4X4+a9Ys2b9/v3Tt2lUeeugh6dixY+DoE6DAqFK/IwKY+hIygAQRQAATVIyMdoX523MBfOKJJ+TXv/51Zbxr164tzZo1c/8+depUmThxosyZM0fOOeccmTBhgqxdu1Y2b94s9evXD/SRIECBMHmxEwLoRRkZREIJIIQJLYzH3WL+9lwAlyxZIq+88soJEdbVv5YtW8rw4cNl9OjR7udHjhwRXSFUMbz11lsDxZ4ABcLkxU4IoBdlZBAJJYAAJrQwHneL+dtzAZw2bZo0bNhQ6tat6y7xTpo0Sb7whS/IG2+8Ie3bt5eXX35ZunTpUhnxPn36SKNGjWTu3LkFY6+SqP/kNg1QmzZt5ODBg9KgQQOPPyoMDQEkAxAIjwACGB5bjlyYAALosQA+++yz8vHHH7vLu3v27HGXeF9//XV59dVX3WXebt26yc6dO91KYG675ZZbZNu2bbJy5cqCiSl0X6HuiAD6/ysGAfS/xowwPgIIYHzss3pmBNBjAcwP9UcffeRW/e6++265+OKLnQDu2rVLWrRoUbnr0KFDZceOHbJixQpWALP6W6GacSOABAIC4RFAAMNjy5FZAawuA7Uq9Ia4jGw9e/aUs88+W0aNGlXSJeB8THyDyEhwRAQBzE6tGWn0BBDA6Jln/YzM3xlaAdR793QFUC/z3nvvve7S71133eVWBHU7evSoNG/enIdAsv5bgRVAEgCByAkggJEjz/wJEUCPBXDkyJHyjW98Q84880x555133D2Aa9askY0bN0rbtm2d6E2ePFlmz54tHTp0cA+IrF69mtfAZP7XQmEArAASDAiERwABDI8tR+YScOYuAQ8cONC91++9995z7/7T+/7Gjx8v559/vmORexH0zJkzj3sRdKdOnQJ/XvgGERhV6ndEAFNfQgaQYAIIYIKL42nXmL89XgGMIrMEKArKyTgHApiMOtALPwkggH7WNcmjYv5GAE35JEAmfKlqjACmqlx0NmUEEMCUFcyD7jJ/I4CmGBMgE75UNUYAU1UuOusBAaTQgyImeAjM3wigKZ4EyIQvVY0RwFSVi856SAAh9LCoMQ6J+RsBNMWPAJnwpaoxApiqctFZDwkggB4WNcYhMX8jgKb4ESATvlQ1RgBTVS46mxECSGFGCh3CMJm/EUBTrAiQCV9iGufLXaFJBQFMTLnoCAQqCSCAhKFUAszfCGCp2XHtCJAJX2IaI4CJKQUdgUBRBBDAonCxcxUCzN8IoOkDQYBM+BLTmNW9xJSCjkCgKAIIYFG42BkBPC4DtSr0T2KwlUQAASwJW+IaIYCJKwkdgkAgAghgIEzsVIAA8zcrgKYPBgEy4UtMYwQwMaWgIxAoigACWBQudmYFkBXAcn0KEMBykYz3OAhgvPw5OwRKJYAAlkqOdszfrACaPgUEyIQvMY0RwMSUgo5AoCgCCGBRuNiZFUBWAMv1KUAAy0Uy3uMggPHy5+wQKJUAAlgqOdoxf7MCaPoUECATvsQ0RgATUwo6AgETAYTQhC9TjZm/EUBT4AmQCV9iGiOAiSkFHYGAiQACaMKXqcbM3wigKfAEyIQvlsbIXizYOSkEIiGAAEaC2YuTMH8jgKYgEyATvlgaI4CxYOekEIiEAAIYCWYvTsL8jQCagkyATPhiaYwAxoKdk0IgEgIIYCSYvTgJ8zcCaAoyATLhi6UxAhgLdk4KgUgIIICRYPbiJMzfCKApyATIhC+WxghgLNg5KQQiIYAARoLZi5MwfyOApiATIBO+WBojgLFg56QQiIQAAhgJZi9OwvyNAJqCTIBM+GJpjADGgp2TQiASAghgJJi9OAnzNwJoCjIBMuGLpTECGAt2TgqBSAgggJFg9uIkzN8IoCnIBMiEr+yNkbuyI+WAEEgVgUICmP97AUlMVUlD6yzzNwJoChcBMuEzN0b4zAg5AAQyRwABzFzJCw6Y+RsBNH0SCJAJn7kxAmhGyAEgkDkCCGDmSo4AVlPyWhUVFRXEoTQCCGBp3MrVCgEsF0mOA4HsEEAAs1Prk42U+ZsVQNMngQCZ8JkbI4BmhBwAApkjgABmruSsALICWP7QI4DlZ1rMERHAYmixLwQgoAQQQHKgBJi/WQE0fRIIkAmfuTECaEbIASCQOQIIYOZKzgogK4DlDz0CWB6mQUQuyOsdytMbjgIBCPhMAAH0ubrBx8b8zQpg8LQU2JMAlYYviPCVdmRaQQACEDg5gSBfJpFE/1PE/I0AmlJOgE7EF+SlqwigKXY0hgAEDASCCGChwyOFBugJbMr8jQCaYkmAEEBTgGgMAQhETqBUAczvKEIYeenKekLmbwTQFCgCVLMAmgDTGAIQgEBCCSCACS1MwG4xfyOAAaNSeDcChACaAkRjCEDAKwLlksIgt9J4BS6GwTB/I4Cm2GUtQEF+KXF/nylSNIYABFJMAAFMT/GyNn8Xqgx/Cs6Q16wFCAE0hIWmEICA9wRKEcAgX5pLOa73sI0DzNr8jQAWIDB9+nSZNm2a7N69Wzp27CgPPvigXHbZZYGilbUAIYCBYsFOEIAABByBIOKGAMYTlqzN3whgHoHHHntMbrjhBlEJ7Natm8ycOVMeffRR2bRpk5x55pk1pjLsAAURrho7GXCHIL+EAh6K3SAAAQhAoIwEgohkGU+XiUOFPX+nAWKmLwF37dpVLrzwQnn44Ycra3XeeedJ3759ZfLkyTXWL+wAIYA1loAdIAABCHhPIIgAFvoSH6Rd0uBFNe+FPX8njSsrgFUIHD16VD772c/K448/Ltddd13lT+6880555ZVXZM2aNSfwOnLkiOg/ue3gwYNupXDHjh3SoEGDste7049XHnfMP47rXZZz5B+3LAflIBCAAAQgEAqBQr/7S/k9Xq45JJRB/t9Bw5r38vusAtimTRs5cOCANGzYMMwhJfbYmV0B3LVrl7Rq1UpefPFFufTSSysLNGnSJJk7d65s3rz5hKKNHTtWxo0bl9hi0jEIQAACEIAABIIT0AWc1q1bB2/g0Z6ZF8B169bJJZdcUlnSiRMnyrx58+T111+vcQXw008/lX379knTpk2lVq1aJ+yf+4YR1gph3Dn0eXyMLe50lXZ+6lYat7hbUbe4K1Da+dNct4qKCvnggw+kZcuWcsopp5QGIOWtMiuApVwCLrbWvt9j4PP4GFuxaU/G/tQtGXUothfUrVhiydjf57olg3C4vcisACpWfQjky1/+snsKOLedf/750qdPn0APgdRUGt8/HD6Pj7HVlO5k/py6JbMuNfWKutVEKJk/97luySRe3l5lWgBzr4GZMWOGuww8a9YseeSRR+TVV1+Vtm3bmkn7/uHweXyMzRz/WA5A3WLBbj4pdTMjjOUAPtctFqARnzTTAqisdfXv/vvvdy+C7tSpk/zsZz+Tyy+/vCxl0CeG9XUyP/zhD6Vu3bplOWaSDuLz+BhbkpIWvC/ULTirJO1J3ZJUjeB98bluwSmkd8/MC2B6S0fPIQABCEAAAhCAQGkEEMDSuNEKAhCAAAQgAAEIpJYAApja0tFxCEAAAhCAAAQgUBoBBLA0brSCAAQgAAEIQAACqSWAAKa2dHQcAhCAAAQgAAEIlEYAASyNm6mVPjml7yD8wx/+IBs2bJAvfelLpuMlofE3v/lN9zeU33nnHWncuLFceeWVMnXqVPeW9TRvW7dulfHjx8uqVavk7bffduP5zne+I/fcc4/UqVMnzUNzfde/fPPMM8+42ul49O9ipnXTJ/qnTZvmnujv2LGjPPjgg3LZZZeldTiV/V67dq0b1/r1693YFi9eLH379k39uHQA+paEp556yv3lpdNOO839WU79vXHuueemfnwPP/yw6D/6O0Q3zeR9990nV199derHlj8AreOYMWPkzjvvdJ87tnQQQABjqJN+SP70pz/Js88+640A6utz9F2KLVq0kJ07d8rIkSMdWf1Te2neVqxYIfq+yG9961ty9tlnyx//+EcZOnSo3HDDDfKTn/wkzUNzff/xj38sjRo1krfeekt+8YtfpFYAc+/0VAns1q2bzJw5Ux599FHZtGmTnHnmmamuk/6e0L9ZfuGFF0q/fv28EsCrrrpKBg4cKBdddJEcO3bMfbHauHGjq9vpp5+e6rotXbpUateu7X5v6KZ/Y15FXr/0qwz6sv3ud7+T/v37S4MGDaR79+4IYIoKiwBGXCz9Zf79739fnnzySfdLwJcVwHyMv/zlL90qha52fuYzn4mYcrin01/i+s3+jTfeCPdEER59zpw5Mnz48NQKoK6oqyBpXXLbeeed5zKoqxO+bPo3x31aAcyvy7vvvivNmzeXNWvWlO19rEmqfZMmTZwE3nTTTUnqVsl9+fDDD93nTr94TZgwwV3NYgWwZJyRN0QAI0S+Z88e96fnlixZIp/73OekXbt2Xgrgvn37ZNiwYW4l8IUXXoiQcDSn+tGPfiS6Mvj73/8+mhNGcJY0C2AUf9c7ghIEOoXvArhlyxbp0KGDWwXUF/P7sn3yySfy+OOPy6BBg9zvfP2Toz5sOh6VWr0C9PWvfx0BTFlREcCIClZRUSHXXHONuzylAqH3hfgmgKNHj5af//zn8vHHH8vFF18sy5Ytk6ZNm0ZEOJrT/PnPf3bfeH/605/KzTffHM1JIzhLmgVw165d0qpVK3eZVO8hy22TJk1yl902b94cAcFoTuGzAOrvSP077Pv375fnn38+GqAhn0VFVm+NOXz4sNSrV08WLFjg5gEftkWLFrlVP/0ifOqppyKAKSwqAmgs2tixY2XcuHEnPYreI6H3wul9SnpDt94XkgYBDDq2r3zlK2787733nujq37Zt2xyThg0bOgnUSStpW7Fj0/6raFxxxRXuH72/LKlbKWPzQQD1M6aTbW7TB1zmzZvnHjDwZfNZAG+//Xb3QJJeNWjdurUXJdPV6e3bt7tbK/S2H/29oZe3074CuGPHDtHf+7/61a+kc+fOrlasAKYvsgigsWYqPfrPybazzjrL3eisNwVXlSG9LKAy+O1vf9utVCRtCzo2/faXv+lDBW3atHHiW3VSTsoYix2byp/e4Kz3mqksnXLKKUkZygn9KHZseoA0CyCXgBMbxcAdu+OOO9ytMfoFWa+M+Lrp2xHat2/vHlJK86a1uu6669z8ldt0PtP5TX836r3fVX+W5rH63HcEMKLq6rfA999/v/JsKhS9e/eWJ554wkmFL994cwPUb4j69OVvfvMb980wzZvey6jyp/dvzp8/38tfbGkWQM2Wfoa0Pnozem7TVRa9pMhDIMn99OllX5U/fbBl9erV7v4/n7cePXq4L8b6eUvz9sEHH7grPVW3IUOGyBe/+EXRW4F8un8zzXWqqe8IYE2EQvp5Gi4BBx36b3/7W9F/vva1r7l3AOrTsfq+K31n2auvvip169YNeqjE7Ze77Ksy++///u/Hyd/nP//5xPW32A7pFxO9bK9PbevTibl7r/TVFXrPUlq23GtgZsyY4VacZ82aJY888ojLX9u2bdMyjIL91Cct9eEI3bp06SIPPPCA+0KiN9+n/RU3t912m7sv7umnnz7u3X96+4i+FzDNm74XT9/5p8KnwqT3zE2ZMsU9QNazZ880D61g37kEnL6SIoAx1cwnAdQbnfXdhvpi648++si9C1Df76UPu+jN+Wne9Ju6frMttOnqRdq3wYMHF7z9II0rt7r6d//997svHroCoU8mXn755WkvkVsZU+HL3/QJzLSvJFV3f/Ds2bNFs5nmTV/18txzz7k8qtBecMEFbnXMR/nTOiGA6UsrApi+mtFjCEAAAhCAAAQgYCKAAJrw0RgCEIAABCAAAQikjwACmL6a0WMIQAACEIAABCBgIoAAmvDRGAIQgAAEIAABCKSPAAKYvprRYwhAAAIQgAAEIGAigACa8NEYAhCAAAQgAAEIpI8AApi+mtFjCEAAAhCAAAQgYCKAAJrw0RgCEIAABCAAAQikjwACmL6a0WMIQAACEIAABCBgIoAAmvDRGAIQgAAEIAABCKSPAAKYvprRYwhAAAIQgAAEIGAigACa8NEYAhCAAAQgAAEIpI8AApi+mtFjCEAAAhCAAAQgYCKAAJrw0RgCEIAABCAAAQikjwACmL6a0WMIQAACEIAABCBgIoAAmvDRGAIQgAAEIAABCKSPAAKYvprRYwhAAAIQgAAEIGAigACa8NEYAhCAAAQgAAEIpI8AApi+mtFjCEAAAhCAAAQgYCKAAJrw0RgCEIAABCAAAQikjwACmL6a0WMIQAACEIAABCBgIoAAmvDRGAIQgAAEIAABCKSPAAKYvprRYwhAAAIQgAAEIGAigACa8NEYAhCAAAQgAAEIpI8AApi+mtFjCEAAAhCAAAQgYCKAAJrw0RgCEIAABCAAAQikjwACmL6a0WMIQAACEIAABCBgIoAAmvDRGAIQgAAEIAABCKSPAAKYvprRYwhAAAIQgAAEIGAigACa8NEYAhCAAAQgAAEIpI8AApi+mtFjCEAAAhCAAAQgYCKAAJrw0RgCEIAABCAAAQikj8D/D2QWrnyOadFlAAAAAElFTkSuQmCC\" width=\"640\">"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "%matplotlib notebook\n",
    "import matplotlib.pyplot as plt\n",
    "plt.hist(x, 100)\n",
    "plt.show()\n",
    "\n",
    "plt.plot(0, 300, 'o')\n",
    "plt.title(r'Normal distribution with $\\mu=0, \\sigma=1$')\n",
    "plt.savefig('matplotlib histogram.png')\n"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "d28d117b",
   "metadata": {},
   "source": [
    "<h1>Canadian Immigration Data</h1>"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 29,
   "id": "985c3e9b",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Requirement already satisfied: xlrd in c:\\users\\mark\\anaconda3\\lib\\site-packages (2.0.1)\n",
      "Requirement already satisfied: openpyxl in c:\\users\\mark\\anaconda3\\lib\\site-packages (3.0.10)\n",
      "Requirement already satisfied: et_xmlfile in c:\\users\\mark\\anaconda3\\lib\\site-packages (from openpyxl) (1.1.0)\n",
      "Collecting js\n",
      "  Downloading js-1.0.tar.gz (2.5 kB)\n",
      "  Preparing metadata (setup.py): started\n",
      "  Preparing metadata (setup.py): finished with status 'done'\n",
      "Collecting fanstatic\n",
      "  Downloading fanstatic-1.2-py2.py3-none-any.whl (31 kB)\n",
      "Requirement already satisfied: setuptools in c:\\users\\mark\\anaconda3\\lib\\site-packages (from js) (63.4.1)\n",
      "Collecting shutilwhich\n",
      "  Downloading shutilwhich-1.1.0.tar.gz (2.3 kB)\n",
      "  Preparing metadata (setup.py): started\n",
      "  Preparing metadata (setup.py): finished with status 'done'\n",
      "Collecting WebOb>=1.2\n",
      "  Downloading WebOb-1.8.7-py2.py3-none-any.whl (114 kB)\n",
      "     -------------------------------------- 115.0/115.0 kB 6.5 MB/s eta 0:00:00\n",
      "Building wheels for collected packages: js, shutilwhich\n",
      "  Building wheel for js (setup.py): started\n",
      "  Building wheel for js (setup.py): finished with status 'done'\n",
      "  Created wheel for js: filename=js-1.0-py3-none-any.whl size=2884 sha256=0a1ebb97d6625025471847caf727549b8ac261da285c742260b0c9c93a8198de\n",
      "  Stored in directory: c:\\users\\mark\\appdata\\local\\pip\\cache\\wheels\\6f\\91\\12\\9fc79cc62b07127faf39b5f3afcc6606e659bb54743a00bebb\n",
      "  Building wheel for shutilwhich (setup.py): started\n",
      "  Building wheel for shutilwhich (setup.py): finished with status 'done'\n",
      "  Created wheel for shutilwhich: filename=shutilwhich-1.1.0-py3-none-any.whl size=2767 sha256=27e07ac4c81c344e2d608db7c9fc6d442eefcc9bb525a9d0cb7dbad6748e2b2d\n",
      "  Stored in directory: c:\\users\\mark\\appdata\\local\\pip\\cache\\wheels\\84\\c7\\f5\\fed66dce1ed897b44e0da776b6a592dfad0a70f7dd61f73a9d\n",
      "Successfully built js shutilwhich\n",
      "Installing collected packages: shutilwhich, WebOb, fanstatic, js\n",
      "Successfully installed WebOb-1.8.7 fanstatic-1.2 js-1.0 shutilwhich-1.1.0\n"
     ]
    }
   ],
   "source": [
    "!pip install xlrd\n",
    "!pip install openpyxl\n",
    "!pip install js\n",
    "import pandas as pd\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 30,
   "id": "9aebd841",
   "metadata": {},
   "outputs": [
    {
     "ename": "ImportError",
     "evalue": "cannot import name 'fetch' from 'js' (C:\\Users\\Mark\\anaconda3\\lib\\site-packages\\js\\__init__.py)",
     "output_type": "error",
     "traceback": [
      "\u001b[1;31m---------------------------------------------------------------------------\u001b[0m",
      "\u001b[1;31mImportError\u001b[0m                               Traceback (most recent call last)",
      "\u001b[1;32m~\\AppData\\Local\\Temp\\ipykernel_112380\\1962386013.py\u001b[0m in \u001b[0;36m<module>\u001b[1;34m\u001b[0m\n\u001b[1;32m----> 1\u001b[1;33m \u001b[1;32mfrom\u001b[0m \u001b[0mjs\u001b[0m \u001b[1;32mimport\u001b[0m \u001b[0mfetch\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m      2\u001b[0m \u001b[1;32mimport\u001b[0m \u001b[0mio\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m      3\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m      4\u001b[0m \u001b[0mURL\u001b[0m \u001b[1;33m=\u001b[0m \u001b[1;34m'https://cf-courses-data.s3.us.cloud-object-storage.appdomain.cloud/IBMDeveloperSkillsNetwork-DV0101EN-SkillsNetwork/Data%20Files/Canada.xlsx'\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m      5\u001b[0m \u001b[0mresp\u001b[0m \u001b[1;33m=\u001b[0m \u001b[1;32mawait\u001b[0m \u001b[0mfetch\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mURL\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
      "\u001b[1;31mImportError\u001b[0m: cannot import name 'fetch' from 'js' (C:\\Users\\Mark\\anaconda3\\lib\\site-packages\\js\\__init__.py)"
     ]
    }
   ],
   "source": [
    "from js import fetch\n",
    "import io\n",
    "\n",
    "URL = 'https://cf-courses-data.s3.us.cloud-object-storage.appdomain.cloud/IBMDeveloperSkillsNetwork-DV0101EN-SkillsNetwork/Data%20Files/Canada.xlsx'\n",
    "resp = await fetch(URL)\n",
    "text = io.BytesIO((await resp.arrayBuffer()).to_py())\n",
    "\n",
    "df_can = pd.read_excel(\n",
    "    text,\n",
    "    sheet_name='Canada by Citizenship',\n",
    "    skiprows=range(20),\n",
    "    skipfooter=2)\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "35fecb98",
   "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.13"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}