wave.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "2023-02-05 15:26:40.198745: I tensorflow/core/platform/cpu_feature_guard.cc:193] This TensorFlow binary is optimized with oneAPI Deep Neural Network Library (oneDNN) to use the following CPU instructions in performance-critical operations:  SSE4.1 SSE4.2 AVX AVX2 FMA\n",
      "To enable them in other operations, rebuild TensorFlow with the appropriate compiler flags.\n"
     ]
    }
   ],
   "source": [
    "from modules.models import create_dense_model, WavePinn\n",
    "from modules.plots import plot_training_loss_linlog, plot_wave_model, plot_wave_model_slider, plot_wave_at_x\n",
    "from modules.data import simulate_wave\n",
    "import tensorflow as tf\n",
    "import numpy as np\n",
    "import matplotlib.pyplot as plt"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "%matplotlib notebook"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "#https://personal.math.ubc.ca/~feldman/m267/separation.pdf\n",
    "c = 1.0\n",
    "length = 1.0\n",
    "n_samples = 2000\n",
    "\n",
    "def f_u(tx):\n",
    "    t = tx[:, 0:1]\n",
    "    x = tx[:, 1:2]\n",
    "    return tf.sin(5 * np.pi * x) * tf.cos(5 * c * np.pi * t) + \\\n",
    "        2*tf.sin(7 * np.pi * x) * tf.cos(7 * c * np.pi * t)\n",
    "\n",
    "def f_u_init(tx):\n",
    "    x = tx[:, 1:2]\n",
    "    return tf.sin(5 * np.pi * x) + 2*tf.sin(7 * np.pi * x)\n",
    "\n",
    "def f_du_dt(tx):\n",
    "    return tf.zeros_like(tx[:, 0:1])\n",
    "\n",
    "def f_u_bnd(tx):\n",
    "    return tf.zeros_like(tx[:, 1:2])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "2023-02-05 15:26:42.082383: E tensorflow/stream_executor/cuda/cuda_driver.cc:265] failed call to cuInit: CUDA_ERROR_COMPAT_NOT_SUPPORTED_ON_DEVICE: forward compatibility was attempted on non supported HW\n",
      "2023-02-05 15:26:42.082445: I tensorflow/stream_executor/cuda/cuda_diagnostics.cc:169] retrieving CUDA diagnostic information for host: nimooli-fedora\n",
      "2023-02-05 15:26:42.082456: I tensorflow/stream_executor/cuda/cuda_diagnostics.cc:176] hostname: nimooli-fedora\n",
      "2023-02-05 15:26:42.082659: I tensorflow/stream_executor/cuda/cuda_diagnostics.cc:200] libcuda reported version is: 525.85.5\n",
      "2023-02-05 15:26:42.082703: I tensorflow/stream_executor/cuda/cuda_diagnostics.cc:204] kernel reported version is: 525.78.1\n",
      "2023-02-05 15:26:42.082710: E tensorflow/stream_executor/cuda/cuda_diagnostics.cc:313] kernel version 525.78.1 does not match DSO version 525.85.5 -- cannot find working devices in this configuration\n",
      "2023-02-05 15:26:42.083165: I tensorflow/core/platform/cpu_feature_guard.cc:193] This TensorFlow binary is optimized with oneAPI Deep Neural Network Library (oneDNN) to use the following CPU instructions in performance-critical operations:  SSE4.1 SSE4.2 AVX AVX2 FMA\n",
      "To enable them in other operations, rebuild TensorFlow with the appropriate compiler flags.\n"
     ]
    }
   ],
   "source": [
    "(tx_samples, residual), (tx_init, u_init, du_dt_init), (tx_bndry, u_bndry) = simulate_wave(n_samples, f_u_init, f_du_dt, f_u_bnd)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [],
   "source": [
    "inputs = [tx_samples, tx_init, tx_bndry]\n",
    "outputs = [f_u(tx_samples), residual, u_init, du_dt_init, u_bndry]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [],
   "source": [
    "backbone = create_dense_model([128]*3, 'elu', 'he_normal', n_inputs=2, n_outputs=1)\n",
    "pinn = WavePinn(backbone, c)\n",
    "scheduler = tf.keras.optimizers.schedules.ExponentialDecay(1e-3, 500, 0.93)\n",
    "optimizer = tf.keras.optimizers.Adam(scheduler)\n",
    "pinn.compile(optimizer=optimizer)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "metadata": {},
   "outputs": [],
   "source": [
    "history = pinn.fit(inputs, outputs, epochs=100, batch_size=200, verbose=0)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "plot_training_loss_linlog(history.history)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "plot_wave_model(pinn.backbone, 0, 1.0, 1)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "10/10 [==============================] - 0s 3ms/step\n"
     ]
    },
    {
     "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('tabindex', '0');\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",
       "            'z-index: 2;'\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(\n",
       "        'style',\n",
       "        'box-sizing: content-box;' +\n",
       "            'pointer-events: none;' +\n",
       "            'position: relative;' +\n",
       "            'z-index: 0;'\n",
       "    );\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;' +\n",
       "            'left: 0;' +\n",
       "            'pointer-events: none;' +\n",
       "            'position: absolute;' +\n",
       "            'top: 0;' +\n",
       "            '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",
       "            /* This rescales the canvas back to display pixels, so that it\n",
       "             * appears correct on HiDPI screens. */\n",
       "            canvas.style.width = width + 'px';\n",
       "            canvas.style.height = height + 'px';\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",
       "        /* User Agent sniffing is bad, but WebKit is busted:\n",
       "         * https://bugs.webkit.org/show_bug.cgi?id=144526\n",
       "         * https://bugs.webkit.org/show_bug.cgi?id=181818\n",
       "         * The worst that happens here is that they get an extra browser\n",
       "         * selection when dragging, if this check fails to catch them.\n",
       "         */\n",
       "        var UA = navigator.userAgent;\n",
       "        var isWebKit = /AppleWebKit/.test(UA) && !/Chrome/.test(UA);\n",
       "        if(isWebKit) {\n",
       "            return function (event) {\n",
       "                /* This prevents the web browser from automatically changing to\n",
       "                 * the text insertion cursor when the button is pressed. We\n",
       "                 * want to control all of the cursor setting manually through\n",
       "                 * the 'cursor' event from matplotlib */\n",
       "                event.preventDefault()\n",
       "                return fig.mouse_event(event, name);\n",
       "            };\n",
       "        } else {\n",
       "            return function (event) {\n",
       "                return fig.mouse_event(event, name);\n",
       "            };\n",
       "        }\n",
       "    }\n",
       "\n",
       "    canvas_div.addEventListener(\n",
       "        'mousedown',\n",
       "        on_mouse_event_closure('button_press')\n",
       "    );\n",
       "    canvas_div.addEventListener(\n",
       "        'mouseup',\n",
       "        on_mouse_event_closure('button_release')\n",
       "    );\n",
       "    canvas_div.addEventListener(\n",
       "        'dblclick',\n",
       "        on_mouse_event_closure('dblclick')\n",
       "    );\n",
       "    // Throttle sequential mouse events to 1 every 20ms.\n",
       "    canvas_div.addEventListener(\n",
       "        'mousemove',\n",
       "        on_mouse_event_closure('motion_notify')\n",
       "    );\n",
       "\n",
       "    canvas_div.addEventListener(\n",
       "        'mouseenter',\n",
       "        on_mouse_event_closure('figure_enter')\n",
       "    );\n",
       "    canvas_div.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",
       "    canvas_div.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.canvas_div.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",
       "\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",
       "    if (name === 'button_press') {\n",
       "        this.canvas.focus();\n",
       "        this.canvas_div.focus();\n",
       "    }\n",
       "\n",
       "    // from https://stackoverflow.com/q/1114465\n",
       "    var boundingRect = this.canvas.getBoundingClientRect();\n",
       "    var x = (event.clientX - boundingRect.left) * this.ratio;\n",
       "    var y = (event.clientY - boundingRect.top) * 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",
       "    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\", \"home\"], [\"Back\", \"Back to previous view\", \"fa fa-arrow-left\", \"back\"], [\"Forward\", \"Forward to next view\", \"fa fa-arrow-right\", \"forward\"], [\"\", \"\", \"\", \"\"], [\"Pan\", \"Left button pans, Right button zooms\\nx/y fixes axis, CTRL fixes aspect\", \"fa fa-arrows\", \"pan\"], [\"Zoom\", \"Zoom to rectangle\\nx/y fixes axis\", \"fa fa-square-o\", \"zoom\"], [\"\", \"\", \"\", \"\"], [\"Download\", \"Download plot\", \"fa fa-floppy-o\", \"download\"]];\n",
       "\n",
       "mpl.extensions = [\"eps\", \"jpeg\", \"pgf\", \"pdf\", \"png\", \"ps\", \"raw\", \"svg\", \"tif\", \"webp\"];\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,iVBORw0KGgoAAAANSUhEUgAAA+gAAAH0CAYAAACuKActAAAAAXNSR0IArs4c6QAAIABJREFUeF7s3Ql8VPW9//9PEhJIgCyEHUJABFmrgGURQbl1QyubLVz91wpFXPpoUbRq8aqFUsBqW7BaZLEKolB/7uCFXturoIJwLQpVNmVfEkCWhABZyPJ/fL6TGRMIZJLMl3zPzOs8HnnMZObM93zP83O+mjfnzPdElZSUlAgLAggggAACCCCAAAIIIIAAAgjUqkAUAb1W/dk4AggggAACCCCAAAIIIIAAAkaAgM6BgAACCCCAAAIIIIAAAggggIADAgR0B4pAFxBAAAEEEEAAAQQQQAABBBAgoHMMIIAAAggggAACCCCAAAIIIOCAAAHdgSLQBQQQQAABBBBAAAEEEEAAAQQI6BwDCCCAAAIIIIAAAggggAACCDggQEB3oAh0AQEEEEAAAQQQQAABBBBAAAECOscAAggggAACCCCAAAIIIIAAAg4IENAdKAJdQAABBBBAAAEEEEAAAQQQQICAzjGAAAIIIIAAAggggAACCCCAgAMCBHQHikAXEEAAAQQQQAABBBBAAAEEECCgcwwggAACCCCAAAIIIIAAAggg4IAAAd2BItAFBBBAAAEEEEAAAQQQQAABBAjoHAMIIIAAAggggAACCCCAAAIIOCBAQHegCHQBAQQQQAABBBBAAAEEEEAAAQI6xwACCCCAAAIIIIAAAggggAACDggQ0B0oAl1AAAEEEEAAAQQQQAABBBBAgIDOMYAAAggggAACCCCAAAIIIICAAwIEdAeKQBcQQAABBBBAAAEEEEAAAQQQIKBzDCCAAAIIIIAAAggggAACCCDggAAB3YEi0AUEEEAAAQQQQAABBBBAAAEECOgcAwgggAACCCCAAAIIIIAAAgg4IEBAd6AIdAEBBBBAAAEEEEAAAQQQQAABAjrHAAIIIIAAAggggAACCCCAAAIOCBDQHSgCXUAAAQQQQAABBBBAAAEEEECAgM4xgAACCCCAAAIIIIAAAggggIADAgR0B4pAFxBAAAEEEEAAAQQQQAABBBAgoHMMIIAAAggggAACCCCAAAIIIOCAAAHdgSLQBQQQQAABBBBAAAEEEEAAAQQI6BwDCCCAAAIIIIAAAggggAACCDggQEB3oAh0AQEEEEAAAQQQQAABBBBAAAECOscAAggggAACCCCAAAIIIIAAAg4IENAdKAJdQAABBBBAAAEEEEAAAQQQQICAzjGAAAIIIIAAAggggAACCCCAgAMCBHQHikAXEEAAAQQQQAABBBBAAAEEECCgcwwggAACCCCAAAIIIIAAAggg4IAAAd2BItAFBBBAAAEEEEAAAQQQQAABBAjoHAMIIIAAAggggAACCCCAAAIIOCBAQHegCHQBAQQQQAABBBBAAAEEEEAAAQI6xwACCCCAAAIIIIAAAggggAACDggQ0B0oAl1AAAEEEEAAAQQQQAABBBBAgIDOMYAAAggggAACCCCAAAIIIICAAwIEdAeKQBcQQAABBBBAAAEEEEAAAQQQIKBzDCCAAAIIIIAAAggggAACCCDggAAB3YEi0AUEEEAAAQQQQAABBBBAAAEECOgcAwgggAACCCCAAAIIIIAAAgg4IEBAd6AIdAEBBBBAAAEEEEAAAQQQQAABAjrHAAIIIIAAAggggAACCCCAAAIOCBDQHSgCXUAAAQQQQAABBBBAAAEEEECAgM4xgAACCCCAAAIIIIAAAggggIADAgR0B4pAFxBAAAEEEEAAAQQQQAABBBAgoHMMIIAAAggggAACCCCAAAIIIOCAAAHdgSLQBQQQQAABBBBAAAEEEEAAAQQI6BwDCCCAAAIIIIAAAggggAACCDggQEB3oAh0AQEEEEAAAQQQQAABBBBAAAECuqVj4KOPPpKnn35a1q1bJ5mZmfL222/LsGHDzru1lStXygMPPCAbN26Uli1bysMPPyz33HNPuc+8+eab8vjjj8v27dulffv2MnXqVBk+fLilvaBZBBBAAAEEEEAAAQQQQACBCyVAQLckvXz5clm1apX07NlTbrnllkoD+s6dO6Vbt24ybtw4ufvuu81nf/7zn8vixYvN53X59NNPZcCAATJlyhQTyjX0P/HEE/LJJ59Inz59LO0JzSKAAAIIIIAAAggggAACCFwIAQL6BVCOioqqNKA/8sgjsmTJEtm8eXOgR3r2fMOGDSaY6zJq1Cg5fvy4aPj3LzfccIOkpKSYIM+CAAIIIIAAAggggAACCCDgXQEC+gWoXTABfeDAgdKjRw955plnAj3SM+QjR46UU6dOSWxsrLRp00YmTJhgfvzLjBkzZObMmbJ79+4K9yQ/P1/0x78UFxfL0aNHJTU1VbRfLAgggAACCCCAAAIIIBCZAiUlJZKTk2O+XhsdHR2ZCI7tNQH9AhQkmIDesWNHGT16tDz66KOBHq1evVr69+8vGRkZ0qJFC4mLi5P58+fLbbfdFlhn0aJFMmbMmHIhvOwuTZo0SSZPnnwB9pJNIIAAAggggAACCCCAgBcF9u7dK61bt/Zi18OuzwT0C1DSYAO6Bu2JEycGeqTfQ7/yyivNJHPNmzc3AX3BggVy6623BtZ59dVXZezYsZKXl1fhnpx5Bj07O9uciddBmJiYeAH2nk0ggAACCCCAAAIIIICAiwL69dm0tDTJysqSpKQkF7sYcX0ioF+AkgcT0G1d4n7m7ukg1MGnQZ2AfgGKzyYQQAABBBBAAAEEEHBUgGzgXmEI6BegJsEEdJ0kbunSpbJp06ZAj+69915Zv359uUni9Dsiy5YtC6wzePBgSU5ODnqSOAbhBSg4m0AAAQQQQAABBBBAwAMCZAP3ikRAt1STEydOyLZt20zrOvnbn/70Jxk0aJA0atTIXGKul7Lv379fXn75ZbOO/zZreos1vdWaztyus7iXvc2afiddz7Trvc+HDh0q7777rjz22GNVus0ag9BSwWkWAQQQQAABBBBAAAGPCZAN3CsYAd1STVasWGEC+ZnLHXfcYSZ60wnhdu3aJbqef1m5cqWZoX3jxo1mJkU9q64hvezyxhtvmFC+Y8cOad++vQnrI0aMCHovGIRBU7EiAggggAACCCCAAAJhLUA2cK+8BHT3amK1RwxCq7w0jgACCCCAAAIIIICAZwTIBu6VioDuXk2s9ohBaJWXxhFAAAEEEEAAAQQQ8IwA2cC9UhHQ3auJ1R4xCK3y0jgCCCCAAAIIIIAAAp4RIBu4VyoCuns1sdojBqFVXhpHAAEEEEAAAQQQQMAzAmQD90pFQHevJlZ7xCC0ykvjCCCAAAIIIIAAAgh4RoBs4F6pCOju1cRqjxiEVnlpHAEEEEAAAQQQQAABzwiQDdwrFQHdvZpY7RGD0CovjSOAAAIIIIAAAggg4BkBsoF7pSKgu1cTqz1iEFrlpXEEEEAAAQQQQAABBDwjQDZwr1QEdPdqYrVHDEKrvDSOAAIIIIAAAggggIBnBMgG7pWKgO5eTaz2iEFolZfGEUAAAQQQQAABBBDwjADZwL1SEdDdq4nVHjEIrfLSOAIIIIAAAggggAACnhEgG7hXKgK6ezWx2iMGoVVeGkcAAQQQQAABBBBAwDMCZAP3SkVAd68mVnvEILTKS+MIIIAAAggggAACCHhGgGzgXqkI6O7VxGqPGIRWeWkcAQQQQAABBBBAAAHPCJAN3CsVAd29mljtEYPQKi+NI4AAAggggAACCCDgGQGygXulIqC7VxOrPWIQWuWlcQQQQAABBBBAAAEEPCNANnCvVAR092pitUcMQqu8NI4AAggggAACCCCAgGcEyAbulYqA7l5NrPaIQWiVl8YRQAABBBBAAAEEEPCMANnAvVIR0N2ridUeMQit8tI4AggggAACCCCAAAKeESAbuFcqArp7NbHaIwahVV4aRwABBBBAAAEEEEDAMwJkA/dKRUB3ryZWe8QgtMpL4wgggAACCCCAAAIIeEaAbOBeqQjo7tXEao8YhFZ5aRwBBBBAAAEEEEAAAc8IkA3cKxUB3b2aWO0Rg9AqL40jgAACCCCAAAIIIOAZAbKBe6UioLtXE6s9YhBa5aVxBBBAAAEEEEAAAQQ8I0A2cK9UBHT3amK1RwxCq7w0jgACCCCAAAIIIICAZwTIBu6VioDuXk2s9ohBaJWXxhFAAAEEEEAAAQQQ8IwA2cC9UhHQ3auJ1R4xCK3y0jgCCCCAAAIIIIAAAp4RIBu4VyoCuns1sdojBqFVXhpHAAEEEEAAAQQQQMAzAmQD90pFQHevJlZ7xCC0ykvjCCCAAAIIIIAAAgh4RoBs4F6pCOgWazJr1ix5+umnJTMzU7p27SozZ86UAQMGVLjF0aNHy4IFC856r0uXLrJx40bz+vz582XMmDFnrZObmyv16tULak8YhEExsRICCCCAAAIIIIAAAmEvQDZwr8QEdEs1ee211+T2228XDen9+/eXOXPmyAsvvCCbNm2SNm3anLXV7Oxs0aDtXwoLC+XSSy+VX/7ylzJp0qRAQL/vvvtk69at5T7fvHnzoPeCQRg0FSsigAACCCCAAAIIIBDWAmQD98pLQLdUkz59+kjPnj3l+eefD2yhc+fOMmzYMJk+fXqlW33nnXdkxIgRsnPnTklPTw8E9Pvvv1+ysrIq/fy5VmAQVpuODyKAAAIIIIAAAgggEFYCZAP3yklAt1CTgoICSUhIkNdff12GDx8e2IKe/V6/fr2sXLmy0q3efPPNkp+fL++//35gXb3E/c4775RWrVpJUVGRXHbZZTJlyhTp0aPHOdvTNvTHv+ggTEtLEz1jn5iYWGk/WAEBBBBAAAEEEEAAAQTCU4CA7l5dCegWapKRkWFC9KpVq+SKK64IbGHatGnme+ZnXqJ+Zhf0O+saohctWiQjR44MvL1mzRrZtm2bdO/eXXQwPfPMM7Js2TLZsGGDdOjQocI90cvjJ0+efNZ7BHQLhadJBBBAAAEEEEAAAQQ8JEBAd69YBHQLNfEH9NWrV0u/fv0CW5g6daosXLhQtmzZct6t6iXwf/zjH0XbiYuLO+e6xcXF5jL6gQMHyp///OcK1+MMuoUC0yQCCCCAAAIIIIAAAmEgQEB3r4gEdAs1qckl7iUlJdKxY0f54Q9/KDNmzKi0d+PGjZN9+/bJ8uXLK11XV2AQBsXESggggAACCCCAAAIIhL0A2cC9EhPQLdVEJ4nr1auXmcXdv+gt04YOHXreSeJWrFghgwYNki+//FK6det23t5pmO/du7e55P3FF18Mak8YhEExsRICCCCAAAIIIIAAAmEvQDZwr8QEdEs18d9mbfbs2eYy97lz58q8efPMPc11VvaJEyfK/v375eWXXy7XA7012zfffCP6ffMzF/0ued++fc33zXUw6WXtesm8ftddg3owC4MwGCXWQQABBBBAAAEEEEAg/AXIBu7VmIBusSZ69vypp54SnfRNz4brJev6fXFdRo8eLbt27RI9Y+5fdOK2Fi1amMnf9NL1M5cJEybIW2+9JQcOHJCkpCQze7tOAlf2e+6V7Q6DsDIh3kcAAQQQQAABBBBAIDIEyAbu1ZmA7l5NrPaIQWiVl8YRQAABBBBAAAEEEPCMANnAvVIR0N2ridUeMQit8tI4AggggAACCCCAAAKeESAbuFcqArp7NbHaIwahVV4aRwABBBBAAAEEEEDAMwJkA/dKRUB3ryZWe8QgtMpL4wgggAACCCCAAAIIeEaAbOBeqQjo7tXEao8YhFZ5aRwBBBBAAAEEEEAAAc8IkA3cKxUB3b2aWO0Rg9AqL40jgAACCCCAAAIIIOAZAbKBe6UioLtXE6s9YhBa5aVxBBBAAAEEEEAAAQQ8I0A2cK9UBHT3amK1RwxCq7w0jgACCCCAAAIIIICAZwTIBu6VioDuXk2s9ohBaJWXxhFAAAEEEEAAAQQQ8IwA2cC9UhHQ3auJ1R4xCK3y0jgCCCCAAAIIIIAAAp4RIBu4VyoCuns1sdojBqFVXhpHAAEEEEAAAQQQQMAzAmQD90pFQHevJlZ7xCC0ykvjCCCAAAIIIIAAAgh4RoBs4F6pCOju1cRqjxiEVnlpHAEEEEAAAQQQQAABzwiQDdwrFQHdvZpY7RGD0CovjSOAAAIIIIAAAggg4BmB7OxsSU5OFn1MTEz0TL/DuaME9HCubgX7RkCPsIKzuwgggAACCCCAQBgJlJSUSFFxiRQWlz4W6fPicq/53y8uKZHCIv/6xVL29yJ9T9vQ98u0WVzatv9R3yv7mn9d3Yb+aJvmuT6WtqXr+9YT81l/G/7PBj6j72v7pW2Y58USWN+8XuJrw79OScmZ7/s+Y1xMW6XPzWd87etn/Nvxv+9/ryjvlOydOZKA7tAYIaA7VIwL0RUC+oVQZhsIIIAAAggggEDtCviD7OmiEjldXGyC6umiYvPz3XNfuNV1Cs17361rftew6l+/TBv+ABxoq/Q9fd3/moZfbds8lrbhD9W6jlnXhFxffwLr+YN36aN/Xf9n9XMsoRMoziegh04zNC0R0EPj6JlWCOieKRUdRQABBBBAAAHHBDT0alAsKPQFXX0sKH3UcOv/3R+Ev1unRE6Xfsa85g/L/nbKhGffZ0qDbvF3z/3BV7dXLmRrwC70hWHdvi/olph+RdISFSVSJzpKYqKjpE50dOmj73f9iY6KktiY736PiY4260eb9aMkJqr0szG+df2fM6/r50pf0/f8n4uJFvO6tuF/v+z2yj+Xcu0GthEVJdr3suv63vOt719PtxGt65n1/fskZttmndJ2zGdK+6Xr6WfM+mW242vX9/rJnOPSrlVTzqA7NFgI6A4V40J0hYB+IZTZBgIIIIAAAgiESkDPvuZrEC70PeYXFpX53ffc/75/Hd9jkS8wlwnRZj1t7/TZAdt8xh+6C4vM2eSyQVzf02Cslwt7ddFQFhsTXfoTJXX0uQbUmGipExMlsdHRElvHF3A1zOqjeV3fN+uVea30Pd/rvvXMc/9n/Z8pfd8Xjn3t+kP0mW36Xvdtz/88ELhN/8oE8NLtmWBbGpK9Wpfa7DfZoDb1K942Ad29mljtEYPQKi+NI4AAAgggEJYC+h1YDbd5p4skr7BI8k77ngdeK/Ncw6+GY13HH6b96+p7+nn/Ov5wbdYz6/uD+HfB2+VLmvXMZ1xMtO+nTmnwrRNlfteQqa/5n8ea577wGVjXfNb3mr5vHjWglj73v2eCdIyvXf9zf9DWYOx7vbSdM0K2Cdel4VvPtrIgUFaAbODe8UBAd68mVnvEILTKS+MIIIAAAghcMAG93FoDbW5BkeSeLv0p0GBc5rmG6jLv+9/T1zRA6+d8odu3ni98fxfA/e/pmWQXFn9I1YBbt06M1I31BWD/o/91E4x1nTPe8wXmGF9wLrNO4Pcy4dnXVpmgbZ5HSd3Sz/vPBOtlwiwIeFWAbOBe5Qjo7tXEao8YhFZ5aRwBBBBAAIFyAnp59qnTRXIqv0hOFRTKqQJ9LP88t8zrGoj96+Se9q2vAdz/ugniZQJ3bVxurZcg14uNkXqxvpBc9lFf11B7rse6pe/rOhqAfev5gra+Zp6XBmN9zx+S9XV9rpczsyCAQOgEyAahswxVSwT0UEl6pB0GoUcKRTcRQAABBC64gIbpk/lFckIDc36hnMgvNL+fNAFaf9egXSgnNWCbR9/7/uDtf90fwPX3C3XmWc8ia1COj4uReBOefT/6vPxr0YHX/CFa369ngrYvbPs/W+55mdCtl1izIIBAeAiQDdyrIwHdvZpY7RGD0CovjSOAAAIIXGABf6jOyT9tAvWJvELJKX08WRqyc/L8Ydv36AvevsCtj/719HJxW4ueddYgnBAXI/Xj6gSex8fVkYRY3+v+932huk4Fr/kCd4J+XsN0XGnYjo0x30tmQQABBKoqQDaoqpj99Qno9o2d2gKD0Kly0BkEEEAgogX07HJO3mk5nlcox3NPiwZp/V0fj5c++l/TUG2em8fTJojra3q2OtSLfre4ft06Jkg30Me6MeZ3f7g2z/W10hBd9j3/6/oZDdL6GX3Uy7NZEEAAAdcEyAauVUSEgO5eTaz2iEFolZfGEUAAgYgS0Nm1NSxn5/p+jucW+h5LX9PQ7fvdF8D1dX8Q1+c6SVmoFv3ecsN6/kDteyz3u74Xp2Hb916Dev7nvvD9XRgnTIeqJrSDAALuC5AN3KsRAd29mljtEYPQKi+NI4AAAp4T0JnA9Sx0Vu5pyTpVINmnTpvnGqyzzHPfa/7fvwvjp83Z7FAs9eNiJDE+1gTqhvXOfKwjiYHXNFz73teQra9r0NbnnKEORSVoAwEEIk2AbOBexQno7tXEao8YhFZ5aRwBBBCoVYHTRcW+UH2qQI6eLJBjpc/18dipAvO6PtfAbX7XIH7qtBQU1exMtl7GnRQfa340NCfG1zGB2/dcH+v43isN4fq6/u4P2nx/ulYPGzaOAAIRLEA2cK/4BHT3amK1RwxCq7w0jgACCIRMQM9s66XhGrSPnsyXIyc0XOtzX7g2AfxkgRwxQdz3u35Hu7qLzgKelBAryfGxkpygATrOhGh9rq/pe/6Q7Q/j/kDO2evqqvM5BBBAoHYFyAa161/R1gno7tXEao8YhFZ5aRwBBBA4p4AGbp3UTIP2kZP5clgfzU++Cdnm50S+CdomdJ8skMLikiqLRkWJCdIpCXEmXOuj7ydWUur7Qrf/PRO+S9/TWcGj9MMsCCCAAAIRI0A2cK/UBHSLNZk1a5Y8/fTTkpmZKV27dpWZM2fKgAEDKtziihUrZNCgQWe9t3nzZunUqVPg9TfffFMef/xx2b59u7Rv316mTp0qw4cPD3ovGIRBU7EiAgggUKlAcXGJ+W72tyfy5dscDd2+Rw3Yh0t/Dzw/WVCte2Lr97MbNYiTRvXrSqOEWN9jfV/YTq3vC9+N6msY9z1qAI+JJmhXWjxWQAABBBAQsoF7BwEB3VJNXnvtNbn99ttFQ3r//v1lzpw58sILL8imTZukTZs2Z23VH9C3bt0qiYmJgfebNGkiMTEx5vdPP/3UBPwpU6aYUP7222/LE088IZ988on06dMnqD1hEAbFxEoIIBDhArkFRXIoJ8+E7UM5vtDt/12fayA/nFNgAnlVz3Jr4E5tUFdSG/gCdmp933MN177XNID7fteferG+/wewIIAAAgggEGoBskGoRWveHgG95oYVtqCBuWfPnvL8888H3u/cubMMGzZMpk+ffs6AfuzYMUlOTq6wzVGjRpl/5Vq+fHng/RtuuEFSUlJk8eLFQe0JgzAoJlZCAIEwFTiZXygHj+fJweO+wH2o9LH87/nmUvSqLHqpeJMGdaWx/jTUxzjf89LA7X9Nw3d8HIG7KrasiwACCCBgT4BsYM+2ui0T0Ksrd57PFRQUSEJCgrz++uvlLj+/7777ZP369bJy5cpzBvS2bdtKXl6edOnSRR577LFyl73rmfcJEyaYH/8yY8YMc+n87t27g9oTBmFQTKyEAAIeE9D7ceuZ7QPH8+RAdq4J4Pr8YHaeHMzR13yhvCrBu15stDRtWE+aNKwrTRvWNY8awpsm+p5rANdHDd1MkuaxA4buIoAAAggYAbKBewcCAd1CTTIyMqRVq1ayatUqueKKKwJbmDZtmixYsED0MvYzF33to48+kl69ekl+fr4sXLhQZs+eLXrp+8CBA83qcXFxMn/+fLntttsCH1+0aJGMGTPGfKaiRV8v+54OwrS0NMnOzi53Kb0FBppEAAEEQiKgtw7Ts96Z2XmSkZVrwrY+97+mv+sl5xrSg1n0EvNmib7grY/NEjWA1zPB+7vHuube2kyaFowo6yCAAAIIeFWAgO5e5QjoFmriD+irV6+Wfv36BbagE7pp8N6yZUtQW7355pvNH4dLliwJBHQN+Lfeemvg86+++qqMHTvWnHWvaJk0aZJMnjz5rLcI6EGVgJUQQMCygE6yppOoafDWn/3mUQN4rmRoEM/KNeG7JIjsrROjNdPQnVRPmpvg7f+p6/s9yfe7Bm8WBBBAAAEEEOAMuovHAAHdQlWqc4l7Rd3QQP/KK6+IzuSuS3UucecMuoUC0yQCCAQtkF9YJJlZeSZ47z+WK/tKH00gz8417xUUFVfaXmxMlDRPqictkuKlhQZwfZ6oj9/9rpecM3t5pZSsgAACCCCAQECAM+juHQwEdEs10Uni9HJ1ncXdv+j3yocOHVrhJHEVdeNHP/qRHD16VD744APztk4Sl5OTI8uWLQusPnjwYDOpHJPEWSokzSKAwHkFNIDrGe99x07J3qO55nGfBvFjp0wo1xnQKzv7rbfebtawnrRM1p9486MhXMO4vqaPOtt5NLcO42hEAAEEEEAgpAIE9JByhqQxAnpIGM9uxH+bNf0euV7mPnfuXJk3b55s3LhR0tPTZeLEibJ//355+eWXzYd1ojedIE7vl65n4PXM+ZNPPil63/MRI0aYdfSSef0+up5Z16D/7rvvmonkuM2apSLSLAIIiF6CrpOt7T16SvYcPSV7j+Wa5/qjQVwnYKssgOtka62S46VVSoLvMbmetEqJl5YmgMebs+GxMdFoI4AAAggggMAFFiCgX2DwIDZHQA8Cqbqr6Nnzp556SjIzM6Vbt26iM677J3wbPXq07Nq1y0wCp4uupyFeQ3t8fLwJ6hrib7zxxnKbf+ONN0wo37Fjh7Rv396EdX+AD6afDMJglFgHgcgSOFVQaML37iOnZM+RU7L76EnZo2fDS0N4ZZegx8fGSOuUeElrlGAefT+lYTzFd/abydYi65hibxFAAAEEvCFANnCvTgR092pitUcMQqu8NI6AswLH807LrsMnZdeRU6WPJ0vD+Clze7LzLXWio8yZ7jaNEkwIT2sUL2kpvucaxgngzpadjiGAAAIIIHBeAbKBewcIAd29mljtEYPQKi+NI1CrAnqPbw3hOw6flJ3fnpTdR07KziP6eEqOniw4b98S69WR9NT60iY1QdIbJZgwrs81iOv3wetwCXqt1paNI4CfwMmtAAAgAElEQVQAAgggYEOAbGBDtWZtEtBr5ue5TzMIPVcyOoxAOYGCwmLZc/SkbP/2pOzUM+L+QH74ZKVnwnWW83aNE0wQ1xCe3rj0MTVBkhPikEYAAQQQQACBCBMgG7hXcAK6ezWx2iMGoVVeGkcgZAJHTuSbEL7j2xPmjPj2Q75H/a54UfG5bwreuEGctGtcX9qm1pe2gUdfKOf+3yErDw0hgAACCCAQFgJkA/fKSEB3ryZWe8QgtMpL4whUSUBnSNdbkW379oQJ4NtKf7Z/e0KOnTp9zrYS4mLkoib15aLGDUwY9/9oIE+Kj61SH1gZAQQQQAABBCJXgGzgXu0J6O7VxGqPGIRWeWkcgQoF9Iy33pbs64M58k1pCNfnGsTzThdX+Bm9N7jehqx90wZyUeP60l4DeZMG0r5JA2mWWJdZ0TnWEEAAAQQQQKDGAmSDGhOGvAECeshJ3W6QQeh2feidtwX8Z8S3HsiRrQdzfIH84AkTxPMLKw7icTHR5gx4+6b15WIN4E0byMUmlDeQ+LgYb4PQewQQQAABBBBwWoBs4F55COju1cRqjxiEVnlpPIIE9NZkWw4cFw3jGsS3Hjwh3xzMkVMFRRUq1K0TbYJ3B/1p1jDwmJYSzwzpEXTcsKsIIIAAAgi4JEA2cKkavr4Q0N2ridUeMQit8tJ4GArknS4y3w3fnHlcthzIMaF8S2aOHDnHbcv0jLieBb+kmS+IdzQ/DaR1SoLEREeFoRC7hAACCCCAAAJeFSAbuFc5Arp7NbHaIwahVV4a97BASUmJuU3Zpszj5mdzZo5sysg2tzKraNJ0/Y54u9T6cklzXwj3P7ZNTeCMuIePA7qOAAIIIIBAJAmQDdyrNgHdvZpY7RGD0CovjXtEQCdt23n4hGzMOC6b9McE8uNy+ERBhXuQnBArnZsnSqcWDQOPHZo25DviHqk33UQAAQQQQACBigXIBu4dGQR092pitUcMQqu8NO6gQH5hkZmobWNGtny1/7h51LPjuafP/q64XoGuE7Z1aZkknTWMt0iULi0SpWlDZk13sLR0CQEEEEAAAQRqKEA2qCGghY8T0C2gutwkg9Dl6tC3mgpoGNfvh3+5X8N4tvx7X7Z8cyhHTheVnNV0fGyMCeFdTRhPlC4tE+WSZpwVr2kN+DwCCCCAAAIIeEeAbOBerQjo7tXEao8YhFZ5afwCChQUFpsZ1P+9PysQxnU29YrCeFJ8rHRtmSjdWiWZR/1p17gBk7ZdwHqxKQQQQAABBBBwT4Bs4F5NCOju1cRqjxiEVnlp3JKA3l98x+GT8u99WbJhb5Zs2JdtvjeuIf3MJSUh1gTx7qU/+rx1SrxE6axuLAgggAACCCCAAAIBAbKBewcDAd29mljtEYPQKi+Nh0jgUE6erN+TJV9oGN+bJV/uy5ac/MKzWtcz499rTRgPETvNIIAAAggggECECZAN3Cs4Ad29mljtEYPQKi+NV0NA7zOuE7d9URrINZjvz8o9q6V6sdHSrWWSfK91slyaliSXtk6W9NQEzoxXw5yPIIAAAggggAACKkA2cO84IKC7VxOrPWIQWuWl8SAEMrJy5fM9x2Td7mPy+Z4sc6/xM783rlejd2zaUC5LS5bL2iSbMN6xWQPuLx6EL6sggAACCCCAAALBCpANgpW6cOsR0C+ctRNbYhA6UYaI6YR+R1y/K27CuAnkxyQzO++s/W/coK70aJNsAnmPtGTp3jpJGtaLjRgndhQBBBBAAAEEEKgNAbJBbaiff5sEdPdqYrVHDEKrvBHf+PG80yaI/2vXMfnX7qOyfm+W5J0uP5FbTHSUubd4zzbJ0jM9RXq2SWESt4g/cgBAAAEEEEAAgdoQIBvUhjoB3T31WuwRg7AW8cNw0wey82TtziOlgfyYbDlwXErOuOV4ckKs9GqTEgjj+v3xhLg6YajBLiGAAAIIIIAAAt4SIBu4Vy/OoLtXE6s9YhBa5Q3rxktKSmTv0VwTyP9v51FZu/Oo7Dl66qx91onbLk9vJJe3TZHvt02Rixo3kOhobnEW1gcHO4cAAggggAACnhQgG7hXNgK6ezWx2iMGoVXesGpcA7nee3zNjiOydsdRE8oPHC///XHN3V1aJsr32zYyP5enp0jTxHph5cDOIIAAAggggAAC4SpANnCvsgR092pitUcMQqu8nm5cA/muI6dMIP90+xHzeCgnv9w+xcZEmduc9W7XyPxoIGcyN0+Xnc4jgAACCCCAQAQLkA3cKz4B3b2aWO0Rg9Aqr+ca33v0lAnjq7cfljU7zj5DHlcn2kzm1qddqvS5qJH0SEuR+LgYz+0nHUYAAQQQQAABBBA4W4Bs4N5RQUB3ryZWe8QgtMrrfONHTuTL6tJAvmrbkbO+Qx4XE23uO973olTpd1GqufVZvVgCufOFpYMIIIAAAggggEA1BMgG1UCz/BECumVg15pnELpWEbv9OVVQaCZz++Sbw7Jq22HZciCn3Ab1lmd67/Er2vsCud72jEButya0jgACCCCAAAIIuCJANnClEt/1g4DuXk2s9ohBaJW31hsvLi6RjRnH5aNvvjWhfN3uY1JQVP4+5J2aN5T+FzeW/henSu92qdKgLrc8q/XC0QEEEEAAAQQQQKAWBMgGtYBeySYJ6O7VxGqPGIRWeWul8czsXPn468MmlOtZ8mOnTpfrR6vkeBnQQQN5Y+nXPlUaN6hbK/1kowgggAACCCCAAAJuCZAN3KqH9oaA7l5NrPaIQWiV94I0nl9YJP/adUxWfv2trNz6rWw9WP6ydT0jrkFcQ/mADk2kbWqCREVxH/ILUhw2ggACCCCAAAIIeEiAbOBesQjoFmsya9YsefrppyUzM1O6du0qM2fOlAEDBlS4xbfeekuef/55Wb9+veTn55v1J02aJNdff31g/fnz58uYMWPO+nxubq7UqxfcvacZhBYLbrFpnW19xdZDJpTrJG+nCooCW9PsfWnrZBnYsYkJ5fqd8tiYaIu9oWkEEEAAAQQQQACBcBAgG7hXRQK6pZq89tprcvvtt4uG9P79+8ucOXPkhRdekE2bNkmbNm3O2ur9998vLVu2lEGDBklycrK89NJL8oc//EHWrl0rPXr0MOtrQL/vvvtk69at5T7fvHnzoPeCQRg0Va2ueLqoWD7bdVQ+3HJIPthySLZ/e7Jcf/Qy9as6NpGrLmkiAy5uLCn142q1v2wcAQQQQAABBBBAwHsCZAP3akZAt1STPn36SM+ePc1Zcf/SuXNnGTZsmEyfPj2orepZ9FGjRskTTzwRCOga5LOysoL6fEUrMQirTWf9g9/m5Juz5B9uPWS+U56TXxjYps623is9xYTyqy9pIp2bJ0p0NJetWy8KG0AAAQQQQAABBMJYgGzgXnEJ6BZqUlBQIAkJCfL666/L8OHDA1vQs996CfvKlSsr3WpxcbG0bdtWHn74YfnFL34RCOh33nmntGrVSoqKiuSyyy6TKVOmBM6wV9SoXi6vP/5FB2FaWppkZ2dLYmJipf1gBXsCJSUl5rZn/9x0UP65+aBs2JddbmOp9ePMGfL/6NTUfJc8KT7WXmdoGQEEEEAAAQQQQCDiBAjo7pWcgG6hJhkZGSZEr1q1Sq644orAFqZNmyYLFiw46xL1irqg311/8sknZfPmzdK0aVOzypo1a2Tbtm3SvXt30cH0zDPPyLJly2TDhg3SoUOHCvdEv8c+efLks94joFsofBBNFhQWy//tPGoC+T82HZT9WbnlPtW9VZIM6tTUhPLvtUriLHkQpqyCAAIIIIAAAgggUD0BAnr13Gx+ioBuQdcf0FevXi39+vULbGHq1KmycOFC2bJly3m3unjxYtEz5e+++65cc80151xXz7LrZfQDBw6UP//5zxWuxxl0CwWuYpPZuafNpesayHXW9bKXrteLjZYrL24i13T2hfKmicFN9lfFLrA6AggggAACCCCAAAJnCRDQ3TsoCOgWalKTS9x1cjmdqV0vj7/pppsq7d24ceNk3759snz58krX1RUYhEEx1XilQ8fz5P1NB+V/Nh6QNTuOyOmikkCbOsGbBvJrOjcz9yaPj4up8fZoAAEEEEAAAQQQQACBqgqQDaoqZn99ArolY50krlevXmYWd//SpUsXGTp06DknidMz5z/72c9EH3UyucoW/Q5z7969zSXvL774YmWrm/cZhEExVWulXYdPmkCuP1/szZKS7zK5dGjaQK7t0sz86C3RmOCtWsR8CAEEEEAAAQQQQCCEAmSDEGKGqCkCeoggz2zGf5u12bNnm8vc586dK/PmzZONGzdKenq6TJw4Ufbv3y8vv/yy+aiG8p/+9Kfme+UjRowINBcfHy9JSUnmd/0ued++fc33zXUw6WXtesm8ftddg3owC4MwGKXg1/nmYI4s+/KALP8q00z4VnbR+5Ff37W5XN+1mVzUpEHwjbImAggggAACCCCAAAIXQIBscAGQq7gJAnoVwaqyup49f+qppyQzM1O6desmM2bMMN8X12X06NGya9cuWbFihfn96quvrnB29zvuuMPc/1yXCRMmyFtvvSUHDhwwoV3vj66TwJX9nntl/WMQViZ0/vf9M68v/zJTln11QLYdOhH4QJ3oKOl7UaoJ5Nd2aS7Nk/g+ec20+TQCCCCAAAIIIICATQGygU3d6rVNQK+em2c/xSCseuk0lG/MOC7LvsyU5V8dkJ2HTwYaiY2JMrdAG9ytubl8PTkhruob4BMIIIAAAggggAACCNSCANmgFtAr2SQB3b2aWO0RgzB43q8P5sh7GzLkvX9nyo4yoTyuTrRc1bGJ3Ni9ufygczNJrMf9yYNXZU0EEEAAAQQQQAABVwTIBq5U4rt+ENDdq4nVHjEIz8+rZ8c1lC/9d4Z8ffC7y9fr1omWQZc0lRu/18LcDq1B3TpW60TjCCCAAAIIIIAAAgjYFiAb2BauevsE9KqbefoTDMKzy5eZnStLN2TIu+szzKXs/kUvX9cz5T/8Xku5pkszQrmnj3w6jwACCCCAAAIIIHCmANnAvWOCgO5eTaz2iEHo480+ddrMvP7O+v2ydufRwC3RYqKjzL3Jf/i9FnJ9l+aSlMDl61YPSBpHAAEEEEAAAQQQqDUBskGt0Z9zwwR092pitUeRPAjzThfJB1sOyTtf7JcVW7+VgqLigHXvto1kyGUt5cbuLaRRfSZ6s3oQ0jgCCCCAAAIIIICAEwKRnA2cKEAFnSCgu1oZS/2KtEFYXFwin+06Km99vt/Mwp6TXxiQ7dS8oQy9rJXcfGkLaZ2SYEmcZhFAAAEEEEAAAQQQcFMg0rKBm1Uo3ysCuheqFMI+Rsog1Mne3v58n7z1xX7Zdyw3INgqOd6cKR96WUvp1DwxhLI0hQACCCCAAAIIIICAtwQiJRt4qSoEdC9VKwR9DedBqN8rf+/LDHO2fN3uYwEtnXH9pu4tZHjPVqKXskdHR4VAkiYQQAABBBBAAAEEEPC2QDhnA69WhoDu1cpVs9/hNgiLiktk1bbD8v/+tVfe33RQCgp93yvXDD6wYxMZ0bO1XNu5mcTHxVRTjI8hgAACCCCAAAIIIBCeAuGWDcKhSgT0cKhiFfYhXAbh7iMn5Y11++TNdfskIzsvIKDfK7+lZ2tzCXvTxHpVkGFVBBBAAAEEEEAAAQQiSyBcskE4VY2AHk7VDGJfvDwITxUUyvIvD8jr6/bKmh1HA3ubFB8rwy5rKT++PE26tUoKQoFVEEAAAQQQQAABBBBAwMvZIFyrR0AP18qeY7+8OAg3ZmTL4v/bI+9+kRGYhT0qSmRAhyby416t5douzaReLJewR9ihzO4igAACCCCAAAII1FDAi9mghrvs/McJ6M6XKLQd9MogPJFfKEs3ZJhg/u992QGENo0STCi/pVdraZkcH1ocWkMAAQQQQAABBBBAIIIEvJINIqgkQkCPpGqLiMuDsKSkRL7c7ztbvmR9hpwsKDLViY2Jkuu7Npdbe7eRfhelMgt7hB2z7C4CCCCAAAIIIICAHQGXs4GdPXa/VQK6+zUKaQ9dHIT63XIN5K+s3S1f7T8e2N+LGtc3oXxEz1aS2qBuSB1oDAEEEEAAAQQQQACBSBdwMRtEek0I6BF2BLg0CLcdOiGvrNktb36+T3LyCk0l4upEy43dfGfLe7drJFH6ZXMWBBBAAAEEEEAAAQQQCLmAS9kg5Dvn0QYJ6B4tXHW7XduD8HRRsby/8aAJ5p/uOBLYjfTUBPn/+rSRH/dKk5T6cdXdPT6HAAIIIIAAAggggAACQQrUdjYIspsRtRoBPaLKXXvfQT90PE9eXbvHfL/8UE6+UY+OEvlB52byk77pMuDixny3PMKORXYXAQQQQAABBBBAoHYFCOi161/R1gno7tXEao8u9CD8Ys8xmb96lyz7MlNOF5WYfWvcoK7c2jtN/rN3G2nFTOxW603jCCCAAAIIIIAAAgicS+BCZwMqUbkAAb1yo7Ba40IMwvzCIhPI56/eLRv2ZgX8eqWnyB1XtJUbujY33zVnQQABBBBAAAEEEEAAgdoTuBDZoPb2zptbJqB7s27V7rXNQXgoJ09eXbPHXMp++ITvMva4mGi5+dKWMvqKttK9dVK1+80HEUAAAQQQQAABBBBAILQCNrNBaHsaOa0R0COn1mZPbQzCLQeOy18/3invrs+QgqJis51miXXl9r7p5jJ2vaSdBQEEEEAAAQQQQAABBNwSsJEN3NpD7/WGgO69mtWox6EahMXFJbLym29NMP9k2+FAn3q2SZafXdlOru/aXGJjuIy9RsXiwwgggAACCCCAAAIIWBQIVTaw2MWIa5qAHmElr+kgzDtdJG9/sV/++slO0fuY66KzsQ/u3kLGXtlOerZJiTBRdhcBBBBAAAEEEEAAAW8K1DQbeHOv3e41Ad3t+oS8d9UdhMdOFsjLn+6Wlz/dJUdOFph+NahbR/7z+2lm4re0Rgkh7ysNIoAAAggggAACCCCAgD2B6mYDez2iZQJ6hB0DVR2Ee4+eMmfLX/tsr+SeLjJaemu0Mf3byqjvp0nDerERJsjuIoAAAggggAACCCAQHgJVzQbhsddu7wUB3e36hLx3wQ7CjRnZMvejHfLevzOlqNh3//KuLRPlnqvay+BuzaUO3y8PeW1oEAEEEEAAAQQQQACBCykQbDa4kH2K9G0R0CPsCDjfICwpKZHV24/I7JXb5eNvvpv4bUCHxnL3wPbS/+JUiYqKijAxdhcBBBBAAAEEEEAAgfAUIKC7V1cCusWazJo1S55++mnJzMyUrl27ysyZM2XAgAHn3OLKlSvlgQcekI0bN0rLli3l4Ycflnvuuafc+m+++aY8/vjjsn37dmnfvr1MnTpVhg8fHvReVDQINZj/c/Mhee7DbbJhb5ZpSyd+u+l7LeXugRdJt1bcvzxoYFZEAAEEEEAAAQQQQMAjAgR09wpFQLdUk9dee01uv/120ZDev39/mTNnjrzwwguyadMmadOmzVlb3blzp3Tr1k3GjRsnd999t6xatUp+/vOfy+LFi+WWW24x63/66acm4E+ZMsWE8rffflueeOIJ+eSTT6RPnz5B7UnZQVi/QUP57y8zZdaH22TLgRzz+bp1os13y8cNuIiJ34ISZSUEEEAAAQQQQAABBLwpQEB3r24EdEs10cDcs2dPef755wNb6Ny5swwbNkymT59+1lYfeeQRWbJkiWzevDnwnp4937BhgwnmuowaNUp0EC1fvjywzg033CApKSkmyAez+AfhSx9slAXrDsnOwyfNx+rHxchP+qXLnVdeJE0a1g2mKdZBAAEEEEAAAQQQQAABDwsQ0N0rHgHdQk0KCgokISFBXn/99XKXn993332yfv160UvZz1wGDhwoPXr0kGeeeSbwlp4hHzlypJw6dUpiY2PNmfcJEyaYH/8yY8YMc+n87t27g9oT/yBMu///SXTdBEmKjzUzso++oq0kJ8QF1QYrIYAAAggggAACCCCAgPcFCOju1ZCAbqEmGRkZ0qpVK3OZ+hVXXBHYwrRp02TBggWydevWs7basWNHGT16tDz66KOB91avXm0uj9f2WrRoIXFxcTJ//ny57bbbAussWrRIxowZI/n5+RXuib5e9j0dhGlpadLjsXfknmu7yW190s39zFkQQAABBBBAAAEEEEAgsgQI6O7Vm4BuoSb+gK4Bu1+/foEt6IRuCxculC1btlQY0DVoT5w4MfCeBvwrr7zSTDLXvHlzE9A14N96662BdV599VUZO3as5OXlVbgnkyZNksmTJ5/13sHDR6VpaoqFvadJBBBAAAEEEEAAAQQQ8IIAAd29KhHQLdTEpUvcz3UGPTs7WxITEy3sPU0igAACCCCAAAIIIICAFwQI6O5ViYBuqSY6SVyvXr3MLO7+pUuXLjJ06NBzThK3dOlSM8u7f7n33nvNd9bLThKXk5Mjy5YtC6wzePBgSU5OrvIkcQR0S4WnWQQQQAABBBBAAAEEPCJAQHevUAR0SzXx32Zt9uzZ5jL3uXPnyrx588w9ztPT082l7Pv375eXX37Z9MB/mzW9xZreak1Duc7iXvY2a3rJvE4mp5fKa9B/99135bHHHqv2bdY4g26p+DSLAAIIIIAAAggggIAHBAjo7hWJgG6xJnr2/KmnnjLfIdd7nOuM6xqwddEJ4Xbt2iUrVqwI9EBnd9cZ2jXEt2zZUvTWaxrSyy5vvPGGCeU7duyQ9u3bm7A+YsSIoPeCQRg0FSsigAACCCCAAAIIIBDWAmQD98pLQHevJlZ7xCC0ykvjCCCAAAIIIIAAAgh4RoBs4F6pCOju1cRqjxiEVnlpHAEEEEAAAQQQQAABzwiQDdwrFQHdvZpY7RGD0CovjSOAAAIIIIAAAggg4BkBsoF7pSKgu1cTqz1iEFrlpXEEEEAAAQQQQAABBDwjQDZwr1QEdPdqYrVHDEKrvDSOAAIIIIAAAggggIBnBMgG7pWKgO5eTaz2iEFolZfGEUAAAQQQQAABBBDwjADZwL1SEdDdq4nVHjEIrfLSOAIIIIAAAggggAACnhEgG7hXKgK6ezWx2iMGoVVeGkcAAQQQQAABBBBAwDMCZAP3SkVAd68mVnvEILTKS+MIIIAAAggggAACCHhGgGzgXqkI6O7VxGqPGIRWeWkcAQQQQAABBBBAAAHPCJAN3CsVAd29mljtEYPQKi+NI4AAAggggAACCCDgGQGygXulIqC7VxOrPWIQWuWlcQQQQAABBBBAAAEEPCNANnCvVAR092pitUcMQqu8NI4AAggggAACCCCAgGcEyAbulYqA7l5NrPaIQWiVl8YRQAABBBBAAAEEEPCMANnAvVIR0N2ridUeMQit8tI4AggggAACCCCAAAKeESAbuFcqArp7NbHaIwahVV4aRwABBBBAAAEEEEDAMwJkA/dKRUB3ryZWe8QgtMpL4wgggAACCCCAAAIIeEaAbOBeqQjo7tXEao8YhFZ5aRwBBBBAAAEEEEAAAc8IkA3cKxUB3b2aWO0Rg9AqL40jgAACCCCAAAIIIOAZAbKBe6UioLtXE6s9YhBa5aVxBBBAAAEEEEAAAQQ8I0A2cK9UBHT3amK1RwxCq7w0jgACCCCAAAIIIICAZwTIBu6VioDuXk2s9ohBaJWXxhFAAAEEEEAAAQQQ8IwA2cC9UhHQ3auJ1R4xCK3y0jgCCCCAAAIIIIAAAp4RIBu4VyoCuns1sdojBqFVXhpHAAEEEEAAAQQQQMAzAmQD90pFQHevJlZ7xCC0ykvjCCCAAAIIIIAAAgh4RoBs4F6pCOju1cRqjxiEVnlpHAEEEEAAAQQQQAABzwiQDdwrFQHdvZpY7RGD0CovjSOAAAIIIIAAAggg4BkBsoF7pSKgu1cTqz1iEFrlpXEEEEAAAQQQQAABBDwjQDZwr1QEdPdqYrVHDEKrvDSOAAIIIIAAAggggIBnBMgG7pWKgG6hJseOHZPx48fLkiVLTOtDhgyRZ599VpKTkyvc2unTp+Wxxx6TZcuWyY4dOyQpKUmuueYaefLJJ6Vly5aBz1x99dWycuXKcm2MGjVK/va3vwW9FwzCoKlYEQEEEEAAAQQQQACBsBYgG7hXXgK6hZoMHjxY9u3bJ3PnzjWt33XXXdK2bVtZunRphVvLzs6WH/3oRzJu3Di59NJLRQP+/fffL4WFhfKvf/2rXEDv2LGj/Pa3vw28Fh8fbwJ9sAuDMFgp1kMAAQQQQAABBBBAILwFyAbu1ZeAHuKabN68Wbp06SJr1qyRPn36mNb1eb9+/WTLli1yySWXBLXFzz77THr37i27d++WNm3amM/oGfTLLrtMZs6cGVQbFa3EIKw2HR9EAAEEEEAAAQQQQCCsBMgG7pWTgB7imrz44ovywAMPSFZWVrmW9fL2GTNmyJgxY4La4j//+U+57rrrTDuJiYmBgL5x40YpKSmRZs2aiZ6p/81vfiMNGzYMqk1diUEYNBUrIoAAAggggAACCCAQ1gJkA/fKS0APcU2mTZsm8+fPl6+//rpcy3ppuobziRMnVrrFvLw8ufLKK6VTp07yyiuvBNafN2+etGvXTpo3by5fffWVaeviiy+Wf/zjH+dsMz8/X/THv+ggTEtLE72s3h/8K+0QKyCAAAIIIIAAAggggEDYCRDQ3SspAT3ImkyaNEkmT5583rX1svT3339fFixYIFu3bi23bocOHWTs2LHy61//+rxt6IRxP/7xj2XPnj2yYsWK84bodevWyeWXXy762LNnzwrbPVe/CehBFp7VEEAAAQQQQAABBBAIUwECunuFJaAHWZPDhw+L/pxv0YngFi1aVO1L3DWcjxw50szk/sEHH0hqaup5t6eXutetW1cWLlwoOpt7RQtn0IMsMKshgAACCCCAAAIIIBBhAgR09wpOQA9xTfyTxK1du9ZM8qaLPu/bt+95J4nzh/NvvvlGPvzwQ2nSpEmlPdPL3Lt3725uvbNhKUMAACAASURBVDZw4MBK19cVGIRBMbESAggggAACCCCAAAJhL0A2cK/EBHQLNdHJ2zIyMmTOnDmmdb3NWnp6ernbrOn3y6dPny7Dhw83t1O75ZZb5PPPP5f33nvPTADnXxo1aiRxcXGyfft2efXVV+XGG2+Uxo0by6ZNm+TBBx8Uvc2aXlofExMT1J4wCINiYiUEEEAAAQQQQAABBMJegGzgXokJ6BZqcvToURk/frwsWbLEtD5kyBB57rnnRGdy9y9RUVHy0ksvyejRo2XXrl1m8reKFj2brrdX27t3r/zkJz8xk8OdOHHCTPR20003mVncNcQHuzAIg5ViPQQQQAABBBBAAAEEwluAbOBefQno7tXEao8YhFZ5aRwBBBBAAAEEEEAAAc8IkA3cKxUB3b2aWO0Rg9AqL40jgAACCCCAAAIIIOAZAbKBe6UioLtXE6s9YhBa5aVxBBBAAAEEEEAAAQQ8I0A2cK9UBHT3amK1RwxCq7w0jgACCCCAAAIIIICAZwTIBu6VioDuXk2s9ohBaJWXxhFAAAEEEEAAAQQQ8IwA2cC9UhHQ3auJ1R4xCK3y0jgCCCCAAAIIIIAAAp4RIBu4VyoCuns1sdojBqFVXhpHAAEEEEAAAQQQQMAzAmQD90pFQHevJlZ7xCC0ykvjCCCAAAIIIIAAAgh4RoBs4F6pCOju1cRqjxiEVnlpHAEEEEAAAQQQQAABzwiQDdwrFQHdvZpY7RGD0CovjSOAAAIIIIAAAggg4BkBsoF7pSKgu1cTqz1iEFrlpXEEEEAAAQQQQAABBDwjQDZwr1QEdPdqYrVHDEKrvDSOAAIIIIAAAggggIBnBMgG7pWKgO5eTaz2iEFolZfGEUAAAQQQQAABBBDwjADZwL1SEdDdq4nVHjEIrfLSOAIIIIAAAggggAACnhEgG7hXKgK6ezWx2iMGoVVeGkcAAQQQQAABBBBAwDMCZAP3SkVAd68mVnvEILTKS+MIIIAAAggggAACCHhGgGzgXqkI6O7VxGqPGIRWeWkcAQQQQAABBBBAAAHPCJAN3CsVAd29mljtEYPQKi+NI4AAAggggAACCCDgGQGygXulIqC7VxOrPWIQWuWlcQQQQAABBBBAAAEEPCNANnCvVAR092pitUcMQqu8NI4AAggggAACCCCAgGcEyAbulYqA7l5NrPaIQWiVl8YRQAABBBBAAAEEEPCMANnAvVIR0N2ridUeMQit8tI4AggggAACCCCAAAKeESAbuFcqArp7NbHaIwahVV4aRwABBBBAAAEEEEDAMwJkA/dKRUB3ryZWe8QgtMpL4wgggAACCCCAAAIIeEaAbOBeqQjo7tXEao8YhFZ5aRwBBBBAAAEEEEAAAc8IkA3cKxUB3b2aWO0Rg9AqL40jgAACCCCAAAIIIOAZAbKBe6UioLtXE6s9YhBa5aVxBBBAAAEEEEAAAQQ8I0A2cK9UBHT3amK1RwxCq7w0jgACCCCAAAIIIICAZwTIBu6VioDuXk2s9ohBaJWXxhFAAAEEEEAAAQQQ8IwA2cC9UhHQ3auJ1R4xCK3y0jgCCCCAAAIIIIAAAp4RIBu4VyoCuoWaHDt2TMaPHy9LliwxrQ8ZMkSeffZZSU5OPufWRo8eLQsWLCj3fp8+fWTNmjWB1/Lz8+VXv/qVLF68WHJzc+UHP/iBzJo1S1q3bh30XjAIg6ZiRQQQQAABBBBAAAEEwlqAbOBeeQnoFmoyePBg2bdvn8ydO9e0ftddd0nbtm1l6dKl5w3oBw8elJdeeimwTlxcnDRq1Cjw+7333mvamD9/vqSmpsqDDz4oR48elXXr1klMTExQe8IgDIqJlRBAAAEEEEAAAQQQCHsBsoF7JSagh7gmmzdvli5dupgz33oGXBd93q9fP9myZYtccsklFW5Rz6BnZWXJO++8U+H72dnZ0qRJE1m4cKGMGjXKrJORkSFpaWmybNkyuf7664PaEwZhUEyshAACCCCAAAIIIIBA2AuQDdwrMQE9xDV58cUX5YEHHjBhu+yil7fPmDFDxowZc86AruFcz5rruldddZVMnTpVmjZtatb/4IMPzCXtesY8JSUl0Mall14qw4YNk8mTJ1fYrl4Wrz/+RQehhnoN/ImJiSHee5pDAAEEEEAAAQQQQAABrwgQ0N2rFAE9xDWZNm2auQT966+/Ltdyx44dTTifOHFihVt87bXXpEGDBpKeni47d+6Uxx9/XAoLC83l63Xr1pVFixaZz5cN29rQddddJ+3atZM5c+ZU2O6kSZMqDO8E9BAXnuYQQAABBBBAAAEEEPCYAAHdvYIR0IOsybmCbtmPf/bZZ/L++++byd62bt1aruUOHTrI2LFj5de//nVQW8zMzDRh/W9/+5uMGDHinAH92muvlfbt28vs2bMrbJcz6EFxsxICCCCAAAIIIIAAAhEnQEB3r+QE9CBrcvjwYdGf8y06EZye6a7OJe4Vtauh/s4775RHHnmk2pe4n9kugzDIgrMaAggggAACCCCAAAJhLkA2cK/ABPQQ18Q/SdzatWuld+/epnV93rdv3/NOEndmN44cOSKtWrUyM8H/9Kc/Nd8Z10niXnnlFRk5cqRZXc+y6y3WmCQuxEWkOQQQQAABBBBAAAEEIkCAgO5ekQnoFmqit1nTGdb93wvX26zp5eplb7PWqVMnmT59ugwfPlxOnDghegn9LbfcIi1atJBdu3bJo48+Knv27BEN/A0bNjS91Nusvffee+Y77nr7Nb0nugZ5brNmoYg0iQACCCCAAAIIIIBAmAsQ0N0rMAHdQk10pvXx48fLkiVLTOtDhgyR5557zszO7l+ioqLMPc/19mq5ublmJvYvvvjCzP6uIX3QoEEyZcoUM+O6f8nLy5OHHnrIXEavn9FZ3WfNmlVuncp2h0FYmRDvI4AAAggggAACCCAQGQJkA/fqTEB3ryZWe8QgtMpL4wgggAACCCCAAAIIeEaAbOBeqQjo7tXEao8YhFZ5aRwBBBBAAAEEEEAAAc8IkA3cKxUB3b2aWO0Rg9AqL40jgAACCCCAAAIIIOAZAbKBe6UioLtXE6s9YhBa5aVxBBBAAAEEEEAAAQQ8I0A2cK9UBHT3amK1RwxCq7w0jgACCCCAAAIIIICAZwTIBu6VioDuXk2s9ohBaJWXxhFAAAEEEEAAAQQQ8IwA2cC9UhHQ3auJ1R4xCK3y0jgCCCCAAAIIIIAAAp4RIBu4VyoCuns1sdojBqFVXhpHAAEEEEAAAQQQQMAzAmQD90pFQHevJlZ7xCC0ykvjCCCAAAIIIIAAAgh4RoBs4F6pCOju1cRqjxiEVnlpHAEEEEAAAQQQQAABzwiQDdwrFQHdvZpY7RGD0CovjSOAAAIIIIAAAggg4BkBsoF7pSKgu1cTqz1iEFrlpXEEEEAAAQQQQAABBDwjQDZwr1QEdPdqYrVHDEKrvDSOAAIIIIAAAggggIBnBMgG7pWKgO5eTaz2iEFolZfGEUAAAQQQQAABBBDwjADZwL1SEdDdq4nVHjEIrfLSOAIIIIAAAggggAACnhEgG7hXKgK6ezWx2iMGoVVeGkcAAQQQQAABBBBAwDMCZAP3SkVAd68mVnvEILTKS+MIIIAAAggggAACCHhGgGzgXqkI6O7VxGqPGIRWeWkcAQQQQAABBBBAAAHPCJAN3CsVAd29mljtEYPQKi+NI4AAAggggAACCCDgGQGygXulIqC7VxOrPWIQWuWlcQQQQAABBBBAAAEEPCNANnCvVAR092pitUcMQqu8NI4AAggggAACCCCAgGcEyAbulYqA7l5NrPaIQWiVl8YRQAABBBBAAAEEEPCMANnAvVIR0N2ridUeMQit8tI4AggggAACCCCAAAKeESAbuFcqArp7NbHaIwahVV4aRwABBBBAAAEEEEDAMwJkA/dKRUB3ryZWe8QgtMpL4wgggAACCCCAAAIIeEaAbOBeqQjo7tXEao8YhFZ5aRwBBBBAAAEEEEAAAc8IkA3cKxUB3b2aWO0Rg9AqL40jgAACCCCAAAIIIOAZAbKBe6UioLtXE6s9YhBa5aVxBBBAAAEEEEAAAQQ8I0A2cK9UBHT3amK1RwxCq7w0jgACCCCAAAIIIICAZwTIBu6VioBuoSbHjh2T8ePHy5IlS0zrQ4YMkWeffVaSk5PPubWoqKgK33vqqafkoYceMu9dffXVsnLlynLrjRo1Sv72t78FvRcMwqCpWBEBBBBAAAEEEEAAgbAWIBu4V14CuoWaDB48WPbt2ydz5841rd91113Stm1bWbp06Tm3duDAgXLvLV++XMaOHSvbtm2Tiy66KBDQO3bsKL/97W8D68bHx0tSUlLQe8EgDJqKFRFAAAEEEEAAAQQQCGsBsoF75SWgh7gmmzdvli5dusiaNWukT58+pnV93q9fP9myZYtccsklQW1x2LBhkpOTI//7v/8bWF/PoF922WUyc+bMoNqoaCUGYbXp+CACCCCAAAIIIIAAAmElQDZwr5wE9BDX5MUXX5QHHnhAsrKyyrWsl7fPmDFDxowZU+kWDx48KK1bt5YFCxbIbbfdVi6gb9y4UUpKSqRZs2aiZ+p/85vfSMOGDStt078CgzBoKlZEAAEEEEAAAQQQQCCsBcgG7pWXgB7imkybNk3mz58vX3/9dbmW9dJ0DecTJ06sdIv6vfMnn3xSMjIypF69eoH1582bJ+3atZPmzZvLV199Zdq6+OKL5R//+Mc528zPzxf9KRvQ09LSJDs7WxITEyvtCysggAACCCCAAAIIIIBAeAoQ0N2rKwE9yJpMmjRJJk+efN61P/vsM3n//ffNme+tW7eWW7dDhw7mO+W//vWvK91ip06d5NprrzUTy51vWbdunVx++eWijz179qxw1XP1m4BeaRlYAQEEEEAAAQQQQACBsBYgoLtXXgJ6kDU5fPiw6M/5Fp0IbtGiRTW6xP3jjz+WgQMHyvr16+XSSy897/b0Uve6devKwoULRWdzr2jhDHqQBWY1BBBAAAEEEEAAAQQiTICA7l7BCeghrol/kri1a9dK7969Tev6vG/fvkFNEjd69Ghz+fq//vWvSnum63Xv3t3cek1DfTALgzAYJdZBAAEEEEAAAQQQQCD8BcgG7tWYgG6hJjp5m35/fM6cOaZ1vc1aenp6udus6WXs06dPl+HDhwd6oAOkRYsW8sc//lHuueeecj3bvn27vPrqq3LjjTdK48aNZdOmTfLggw+K3mZNL62PiYkJak8YhEExsRICCCCAAAIIIIAAAmEvQDZwr8QEdAs1OXr0qIwfP16WLFliWh8yZIg899xzojO5+5eoqCh56aWXRM+Y+xe9b/r9998vmZmZZ93bfO/evfKTn/zEnF0/ceKE6ERvN910k5nFvVGjRkHvBYMwaCpWRAABBBBAAAEEEEAgrAXIBu6Vl4DuXk2s9ohBaJWXxhFAAAEEEEAAAQQQ8IwA2cC9UhHQ3auJ1R4xCK3y0jgCCCCAAAIIIIAAAp4RIBu4VyoCuns1sdojBqFVXhpHAAEEEEAAAQQQQMAzAmQD90pFQHevJlZ7xCC0ykvjCCCAAAIIIIAAAgh4RoBs4F6pCOju1cRqjxiEVnlpHAEEEEAAAQQQQAABzwiQDdwrFQHdvZpY7RGD0CovjSOAAAIIIIAAAggg4BkBsoF7pSKgu1cTqz1iEFrlpXEEEEAAAQQQQAABBDwjQDZwr1QEdPdqYrVHDEKrvDSOAAIIIIAAAggggIBnBMgG7pWKgO5eTaz2iEFolZfGEUAAAQQQQAABBBDwjADZwL1SEdDdq4nVHjEIrfLSOAIIIIAAAggggAACnhEgG7hXKgK6ezWx2iMGoVVeGkcAAQQQQAABBBBAwDMCZAP3SkVAd68mVnvEILTKS+MIIIAAAggggAACCHhGgGzgXqkI6O7VxGqPGIRWeWkcAQQQQAABBBBAAAHPCJAN3CsVAd29mljtEYPQKi+NI4AAAggggAACCCDgGQGygXulIqC7VxOrPWIQWuWlcQQQQAABBBBAAAEEPCNANnCvVAR092pitUcMQqu8NI4AAggggAACCCCAgGcEyAbulYqA7l5NrPaIQWiVl8YRQAABBBBAAAEEEPCMANnAvVIR0N2ridUeMQit8tI4AggggAACCCCAAAKeESAbuFcqArp7NbHaIwahVV4aRwABBBBAAAEEEEDAMwJkA/dKRUB3ryZWe8QgtMpL4wgggAACCCCAAAIIeEaAbOBeqQjo7tXEao8YhFZ5aRwBBBBAAAEEEEAAAc8IkA3cKxUB3b2aWO0Rg9AqL40jgAACCCCAAAIIIOAZAbKBe6UioLtXE6s9YhBa5aVxBBBAAAEEEEAAAQQ8I0A2cK9UBHT3amK1RwxCq7w0jgACCCCAAAIIIICAZwTIBu6VioDuXk2s9ohBaJWXxhFAAAEEEEAAAQQQ8IwA2cC9UhHQ3auJ1R4xCK3y0jgCCCCAAAIIIIAAAp4RIBu4VyoCuns1sdojBqFVXhpHAAEEEEAAAQQQQMAzAmQD90pFQHevJlZ7xCC0ykvjCCCAAAIIIIAAAgh4RoBs4F6pCOju1cRqjxiEVnlpHAEEEEAAAQQQQAABzwiQDdwrFQHdQk2mTp0q//3f/y3r16+XuLg4ycrKqnQrJSUlMnnyZJk7d64cO3ZM+vTpI3/5y1+ka9eugc/m5+fLr371K1m8eLHk5ubKD37wA5k1a5a0bt260vb9KzAIg6ZiRQQQQAABBBBAAAEEwlqAbOBeeQnoFmrym9/8RpKTk2Xfvn3y17/+NaiA/vvf/1402M+fP186duwov/vd7+Sjjz6SrVu3SsOGDU0v7733Xlm6dKlZJzU1VR588EE5evSorFu3TmJiYoLaEwZhUEyshAACCCCAAAIIIIBA2AuQDdwrMQHdYk00SN9///2VBnQ9e96yZUuz7iOPPGJ6pGfLmzVrJhrc7777bsnOzpYmTZrIwoULZdSoUWadjIwMSUtLk2XLlsn1118f1J4wCINiYiUEEEAAAQQQQAABBMJegGzgXokJ6BZrEmxA37Fjh7Rv314+//xz6dGjR6BHQ4cONWfiFyxYIB988IG5pF3PmKekpATWufTSS2XYsGHm8vhgFgZhMEqsgwACCCCAAAIIIIBA+AuQDdyrMQHdYk2CDeirV6+W/v37y/79+82ZdP9y1113ye7du+V//ud/ZNGiRTJmzBhzZr3sct1110m7du1kzpw5Fe6Jrl/2M3omvk2bNrJ3715JTEy0uPc0jQACCCCAAAIIIIAAAi4LaEDXK3J1zqykpCSXuxoxfSOgB1nqSZMmVXqW+rPPPpPLL7880GJVA7pest6iRYvA58eNG2eC9N///vdzBvRrr73WnH2fPXt2hXsSTL+DJGA1BBBAAAEEEEAAAQQQCEOB7du3y0UXXRSGe+a9XSKgB1mzw4cPi/6cb2nbtq3Uq1evygHd5iXuZ55B138dS09Plz179vCvZEHWntWqJ+D/F1mu1qieH58KXoBjLXgr1qyZAMdazfz4dPACHGvBW7FmzQT8V9fqXaT0q7UstS9AQLdYg2DPoPsniZswYYI8/PDDpkcFBQXStGnTsyaJe+WVV2TkyJFmnczMTHOLNSaJs1hEmq62AN9pqjYdH6yiAMdaFcFYvdoCHGvVpuODVRTgWKsiGKtXW4Bjrdp01j5IQLdAq2endTK3JUuWyNNPPy0ff/yx2crFF18sDRo0MM87deok06dPl+HDh5vfdbZ2/f2ll16SDh06yLRp02TFihVn3WbtvffeM7dZa9Sokbkn+pEjR7jNmoUa0mTNBfgPfs0NaSE4AY614JxYq+YCHGs1N6SF4AQ41oJzYq2aC3Cs1dww1C0Q0EMtKiKjR482M6+fuXz44Ydy9dVXm5ejoqJMGNd1ddGz6DoTu072ppeY9OnTR/7yl79It27dAs3k5eXJQw89ZL6Pnpuba2Z1nzVrlpnYIdiFQRisFOvVVIBjraaCfD5YAY61YKVYr6YCHGs1FeTzwQpwrAUrxXo1FeBYq6lg6D9PQA+9qdMt6nfS9Uz9xIkTpW7duk73lc55W4Bjzdv181LvOda8VC1v95Vjzdv181LvOda8VC1v95Vjzb36EdDdqwk9QgABBBBAAAEEEEAAAQQQiEABAnoEFp1dRgABBBBAAAEEEEAAAQQQcE+AgO5eTegRAggggAACCCCAAAIIIIBABAoQ0COw6OwyAggggAACCCCAAAIIIICAewIEdPdqQo8QQAABBBBAAAEEEEAAAQQiUICAHoZF11uv6f3XMzMzpWvXrjJz5kwZMGDAOfd05cqV8sADD8jGjRulZcuW8vDDD8s999wThjLsUqgFqnKsvfXWW/L888/L+vXrRWcM1WNz0qRJcv3114e6W7QXhgJVOdbK7v6qVavkqquuMres1GOPBYHKBKp6rOl/z37729/KK6+8IgcOHJDWrVvLf/3Xf8nPfvazyjbF+xEuUNVj7dVXX5WnnnpKvvnmG0lKSpIbbrhB/vCHP0hqamqES7L75xL46KOPTCZYt26dyQVvv/22DBs27Lxg5ILaP54I6LVfg5D24LXXXpPbb7/d3B+9f//+5r7qL7zwgmzatEnatGlz1rZ27txp/nAdN26c3H333aJ/zP785z+XxYsXyy233BLSvtFYeAlU9Vi7//77zT8ADRo0SJKTk+Wll14yf1isXbtWevToEV447E1IBap6rPk3np2dLT179pSLL75YDh48SEAPaVXCs7HqHGtDhw41x9fvfvc7c6wdOnRICgsL5YorrghPJPYqJAJVPdY++eQT84+NM2bMkJtvvln2799vTqZ06NDBhC4WBCoSWL58ufnbXv9fqH/XVxbQyQVuHEcEdDfqELJe9OnTxwxCPVPpXzp37mz+tUzvf37m8sgjj8iSJUtk8+bNgbf0P/gbNmyQTz/9NGT9oqHwE6jqsVaRgJ5FHzVqlDzxxBPhB8QehUygusfaf/7nf5o/XmNiYuSdd94hoIesIuHbUFWPtb///e+ix9mOHTukUaNG4QvDnoVcoKrHmv6Dtv5tt3379kBfnn32WXNGfe/evSHvHw2Gn0BUVFSlAZ1c4EbdCehu1CEkvSgoKJCEhAR5/fXXZfjw4YE277vvPvOHqV6ycuYycOBAc/bymWeeCbyl/7o2cuRIOXXqlMTGxoakbzQSXgLVOdbOFCguLpa2bduar1T84he/CC8g9iZkAtU91vQKDb2SSP+hUc9sEtBDVpKwbag6x5pecfb111/L5ZdfLgsXLpT69evLkCFDZMqUKRIfHx+2VuxYzQSqc6ytXr3aXIGmf6MNHjzYXKmhf6vpSZjZs2fXrEN8OiIEggno5AI3DgUCuht1CEkvMjIypFWrVuZSlrKX1k2bNk0WLFggW7duPWs7HTt2lNGjR8ujjz4aeE//J6CXx2t7LVq0CEnfaCS8BKpzrJ0poN+JevLJJ83VG02bNg0vIPYmZALVOdb0+5lXXnmlfPzxx6L/jdO5DgjoIStJ2DZUnWNNvwO8YsUKueaaa8yVQIcPHzZfE/uP//gPefHFF8PWih2rmUB1jjXd4htvvCFjxoyRvLw88zUK/ccgfY2TKTWrR6R8OpiATi5w42ggoLtRh5D0wv8ffA3Y/fr1C7Q5depU8y/7W7ZsqTCg63/sJ06cGHhPA77+cauTSTRv3jwkfaOR8BKozrFWVkDnOLjzzjvl3XffNX/YsiBwLoGqHmtFRUXSt29fGTt2bGCySwI6x1cwAlU91rTN6667zvxDkE4Op5N26aITYv7oRz+SkydPchY9GPgIXKc6x5rOJaT/v5wwYYKZXFX/RnvooYfk+9//vvz1r3+NQEV2uaoCwQZ0ckFVZUO/PgE99Ka11mJ1LpniUpZaK5enN1ydY82/wzoxjv7HX7+KcdNNN3nagc7bF6jqsZaVlSUpKSnme+f+Rb9OUVJSYl57//33zdlNFgTOFKjqsaafv+OOO8xVa9u2bQs0p1cFdenSxVz6rnMgsCAQimNNJwDWM+f6/07/ohPH6V16uOKRYywYgWACOrkgGEn76xDQ7Rtf0C3opCO9evUy3730L/qHgs4ye65J4pYuXWpmefcv9957r/nOOpPEXdDSeW5jVT3WdAf1zLneekgfK7vNh+dA6LA1gaocaxrGy/73TDul/z384IMPzKWg7dq1M98TZkGgIoGqHGv6+blz54reoUK/D9ygQQPTpF4ZNGLECDlx4gRn0DnMzilQ1WNNZ+CuU6eO6D9y+xf9O02/0qgzuutdUlgQOJ9AMAFdJ4kjF9T+cURAr/0ahLQH/tt26IQhepm7/vEwb948c4/z9PR0cym7/of85ZdfNtv1305Bb7Gmt1rT/9jrLO7cZi2kZQnLxqp6rOkx9dOf/tRMSKh/vPoXnUjJf2loWEKxUzUWqOqxduYGucS9xiWImAaqeqxpCNdJuvRrFZMnTzbfQdev7+jtsPT/vSwInEugqsfa/Pnzzd9pf/7znwOXuOs/DkVHR5vblbIgUJGA/jfKf4WPTgr9pz/9yUw2qHed0NsvkwvcPG4I6G7WpUa90rNFetsN/X6S3uNc75mpl6zoohPC7dq1y0xq4190dnf9TpOGeP0XWP3XMw3pLAhUJlCVY+3qq6+u8E4Ceomo/uHBgsD5BKpyrBHQOZZqIlDVY03nd/nlL39pLnVPTU01M2vrnQOYxb0mVYiMz1b1WNPbqukJGD25kpycbL6u8/vf/95MEMyCQEUC+ve+BvIzF//fXuQCN48bArqbdaFXCCCAAAIIIIAAAggggAACESZAQI+wgrO7CCCAAAIIIIAAAggggAACbgoQ0N2sC71CAAEEEEAAAQQQQAABBBCIMAECeoQVnN1FAAEEEEAAAQQQQAABBBBwU4CA7mZd6BUCCCCAAAIIIIAAAggggECECRDQI6zg7C4CCCCAx/IQtQAADztJREFUAAIIIIAAAggggICbAgR0N+tCrxBAAAEEEEAAAQQQQAABBCJMgIAeYQVndxFAAAEEEEAAAQQQQAABBNwUIKC7WRd6hQACCCCAAAIIIIAAAgggEGECBPQIKzi7iwACCCCAAAIIIIAAAggg4KYAAd3NutArBBBAAAEEEEAAAQQQQACBCBMgoEdYwdldBBBAAIHIFfj222+le/fuMn78eHn00UcNxNq1a2XAgAHy3nvvyXXXXRe5OOw5AggggAACDggQ0B0oAl1AAAEEEEDgQgksW7ZMhg0bJqtXr5ZOnTpJjx495KabbpKZM2deqC6wHQQQQAABBBD4/9u7sxivqjsO4D9gRFB2WdIALqEpCkipCgriViqmaQz6YMJDXaJ9KCGoaKgxNpYSlUQpYOK+EAi1idEofdAHQwKairtRhwFNizEgIGBBg4Wwjc25DdOiUmb+MMOZOZ+b+KD533t/5/O7IX455557GAEB3aNBgAABAgQKE5g+fXosX748xo4dGx9++GG888470a1bt8IUDJcAAQIECOQnIKDn1xMVESBAgACBVhXYvXt3jBo1KjZs2BDvvvtujB49ulXv5+IECBAgQIBA8wQE9OY5+RUBAgQIEOgwAg0NDXHeeefFvn374sUXX4wrr7yyw4zNQAgQIECAQHsWENDbc/fUToAAAQIEWiiwd+/eGDduXIwZM6Z6B33+/PlRX18fgwYNauGV/JwAAQIECBA41gIC+rEWdT0CBAgQIJCxwKxZs+L555+v3j3v0aNHXHbZZdGzZ89qF3cHAQIECBAgcHwFBPTj6+/uBAgQIECgzQRWrlwZl19+eaxYsSImTpxY3Xf9+vXVO+hz586NadOmtVktbkSAAAECBAh8X0BA91QQIECAAAECBAgQIECAAIEMBAT0DJqgBAIECBAgQIAAAQIECBAgIKB7BggQIECAAAECBAgQIECAQAYCAnoGTVACAQIECBAgQIAAAQIECBAQ0D0DBAgQIECAAAECBAgQIEAgAwEBPYMmKIEAAQIECBAgQIAAAQIECAjongECBAgQIECAAAECBAgQIJCBgICeQROUQIAAAQIECBAgQIAAAQIEBHTPAAECBAgQIECAAAECBAgQyEBAQM+gCUogQIAAAQIECBAgQIAAAQICumeAAAECBAgQIECAAAECBAhkICCgZ9AEJRAgQIAAAQIECBAgQIAAAQHdM0CAAAECBAgQIECAAAECBDIQENAzaIISCLSWwLfffhv79u2LAwcORJcuXeKEE06ITp06tdbtXJcAAQIECBAgQIAAgaMQENCPAs+pBHIVSKF8y5YtsXHjxti9e3dTmd27d4/BgwfHoEGDqrDuIECAAAECBAgQIEAgHwEBPZ9eqITAMRHYvn17NDQ0RGNjY/Tv3z8GDBgQdXV1sX///ti2bVt8+eWX0blz5xg5cmT069fvmNzTRQgQIECAAAECBAgQOHoBAf3oDV2BQDYCKZzX19dH3759Y/jw4XHiiSd+r7Y9e/bEJ598Ejt27Iizzz5bSM+mewohQIAAAQIECBAoXUBAL/0JMP4OI5CWtb/55pvRu3fvGDVqVDVLfrgjza6vXr06vv7667jgggssd+8wT4GBECBAgAABAgQItGcBAb09d0/tBP5H4PPPP49169ZVgfuHZs6/i5Vm0lOgHzZsWAwZMoQlAQIECBAgQIAAAQLHWUBAP84NcHsCx0Ig7db+9ttvR48ePap3y5t7pHfVv/nmmxg3bpzd3ZuL5ncECBAgQIAAAQIEWklAQG8lWJctU2D27NmxbNmy+OCDD9oUYO/evbFq1aoYMWJEDBw4sNn33rp1a6xZsyYuvPBCy9ybreaHBAgQIECAAAECBFpHQEBvHVdX7YACR/p++PXXXx8PPfRQpKXjp5xySpsKpE+pvfXWWzF69OgWbfqWNpX76KOP4vzzz4/0CTYHAQIECBAgQIAAgZYIPPLII/HAAw/E5s2bq5WcCxcujIsuuqgll/Db/xEQ0D0OBJop8MUXXzT98tlnn42777672g394JECbtqg7XgczZ5B37sronNdRF3Xqkwz6MejW+5JgAABAgQIEOgYAun/ia+99tpIIT2tyHz88cfjqaeeqlZonnrqqR1jkG08CgG9jcHdrmMILF68OG699db46quvDhnQd5e433DDDdVv0jveDz74YDW7PnPmzLjrrrvizjvvjKeffjpOOumkmDNnTtx4441N19q4cWPcdttt8corr1S7sU+cOLE6//TTT/9BwGa/g778jxFrlkX86k8Rw35efS/dO+gd45k0CgIECBAgQIBAWwukVZjnnHNOPProo023Puuss+Kqq66KuXPntnU5HeJ+AnqHaKNBtLVASwL6Cy+8ENddd13MmDEjXn/99bjpppviiiuuiIsvvjiuueaaSH/zmAJ62oF96NChsWvXrhgzZky1NCj9JUBdXV3cc8898d5771XL0bt2/c/s93ePI+7i/q9/RiwYGbF/d8TUv8SeM34Rb7zxRpx88smWt7f1A+R+BAgQIECAAIFMBNLneWs50grONNH03HPPxdVXX910iVtuuaXaj+nVV1+t5bLFnyOgF/8IAKhFoCUBfeXKlfHpp582fZf8zDPPrDZye+2116pbHzhwoFoan5YDTZ06NRYtWhT3339/rF27tmln9fQHYJ8+faoN6CZPnvyDJR/xO+h/Wxix/A8RP/ppNP5mRaxuaIgdO3ZE3759/+8302vxcQ4BAgQIECBAgED7EKg1oG/atCkGDx5cTUBNmDChabD33XdfLFmy5JBXQduHRB5VCuh59EEV7UygJQF927Zt8dJLLzWN8JJLLon0B+HDDz/c9N9OO+20uP322+Pmm2+O6dOnV+/vdOvW7RCVNLOezpk2bdphtdKmb/X19VXoHj58+KHfQ39kfMTWNbHvl/NjbfexkX7bq1evw87It7OWKJcAAQIECBAgQKAGgaMN6OlLQuPHj2+687333htLly6Njz/+uIZqnCKgewYI1CDQkoCe3kFPM98Hj0svvbRawp52uDx4pHfL03L29E8K4O+//34888wz36tswIABR9yILgXv9G55Y2Nj9O/fP9I5XXeujz5LJ0Vjp7pYNWFJHOjaM3r27Cmc19B7pxAgQIAAAQIEOpJArQHdEvfWeQoE9NZxddUOLtCaAf3JJ5+MO+64Iz777LNqhruWIy1337JlS6TN5tIn2IZs+Gv8eN2i2NFvTPxj/LxqZj1tPucgQIAAAQIECBAoW6DWgJ7U0iZx5557brWL+8FjxIgRMWXKFJvE1fhYCeg1wjmtbIHWDOgHN4lL7/SkzeOGDBkS69evj7TZ3KxZs6p/b+6Rdnffv39/dP7twOgytDE2jZ4R238ytbmn+x0BAgQIECBAgEAHFziagH7wM2uPPfZYtcz9iSeeiDTZlFZzplc4HS0XENBbbuYMAtGaAT3xpm+up1n0l19+OXbu3FltwDFp0qSYN29ebbPqv+8VUdcp/j75z7Gn1xk6SIAAAQIECBAgQKASOJqAns5Ps+dpg+PNmzdX11qwYEH1tSJHbQICem1uziLQvgRm947oNThiZkNEp07tq3bVEiBAgAABAgQIEChEQEAvpNGGWbhACug/+3XElP/uHF+4iOETIECAAAECBAgQyE5AQM+uJQoi0AoCKaBfszhi5NWtcHGXJECAAAECBAgQIEDgWAgI6MdC0TUI5C4wp3/E79ZFdOude6XqI0CAAAECBAgQIFCsgIBebOsNnAABAgQIECBAgAABAgRyEhDQc+qGWggQIECAAAECBAgQIECgWAEBvdjWGzgBAgQIECBAgAABAgQI5CQgoOfUDbUQIECAAAECBAgQIECAQLECAnqxrTdwAgQIECBAgAABAgQIEMhJQEDPqRtqIUCAAAECBAgQIECAAIFiBQT0Yltv4AQIECBAgAABAgQIECCQk4CAnlM31EKAAAECBAgQIECAAAECxQoI6MW23sAJECBAgAABAgQIECBAICcBAT2nbqiFAAECBAgQIECAAAECBIoVENCLbb2BEyBAgAABAgQIECBAgEBOAgJ6Tt1QCwECBAgQIECAAAECBAgUKyCgF9t6AydAgAABAgQIECBAgACBnAQE9Jy6oRYCBAgQIECAAAECBAgQKFZAQC+29QZOgAABAgQIECBAgAABAjkJCOg5dUMtBAgQIECAAAECBAgQIFCsgIBebOsNnAABAgQIECBAgAABAgRyEhDQc+qGWggQIECAAAECBAgQIECgWAEBvdjWGzgBAgQIECBAgAABAgQI5CQgoOfUDbUQIECAAAECBAgQIECAQLECAnqxrTdwAgQIECBAgAABAgQIEMhJQEDPqRtqIUCAAAECBAgQIECAAIFiBQT0Yltv4AQIECBAgAABAgQIECCQk4CAnlM31EKAAAECBAgQIECAAAECxQoI6MW23sAJECBAgAABAgQIECBAICcBAT2nbqiFAAECBAgQIECAAAECBIoVENCLbb2BEyBAgAABAgQIECBAgEBOAgJ6Tt1QCwECBAgQIECAAAECBAgUKyCgF9t6AydAgAABAgQIECBAgACBnAQE9Jy6oRYCBAgQIECAAAECBAgQKFZAQC+29QZOgAABAgQIECBAgAABAjkJCOg5dUMtBAgQIECAAAECBAgQIFCsgIBebOsNnAABAgQIECBAgAABAgRyEhDQc+qGWggQIECAAAECBAgQIECgWAEBvdjWGzgBAgQIECBAgAABAgQI5CQgoOfUDbUQIECAAAECBAgQIECAQLECAnqxrTdwAgQIECBAgAABAgQIEMhJQEDPqRtqIUCAAAECBAgQIECAAIFiBQT0Yltv4AQIECBAgAABAgQIECCQk4CAnlM31EKAAAECBAgQIECAAAECxQoI6MW23sAJECBAgAABAgQIECBAICcBAT2nbqiFAAECBAgQIECAAAECBIoVENCLbb2BEyBAgAABAgQIECBAgEBOAgJ6Tt1QCwECBAgQIECAAAECBAgUKyCgF9t6AydAgAABAgQIECBAgACBnAQE9Jy6oRYCBAgQIECAAAECBAgQKFZAQC+29QZOgAABAgQIECBAgAABAjkJCOg5dUMtBAgQIECAAAECBAgQIFCsgIBebOsNnAABAgQIECBAgAABAgRyEhDQc+qGWggQIECAAAECBAgQIECgWIF/A+np/CI1Mt/WAAAAAElFTkSuQmCC\" width=\"1000\">"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "slider = plot_wave_model_slider(backbone, 0, 1, 1)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "1/1 [==============================] - 0s 15ms/step\n"
     ]
    },
    {
     "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('tabindex', '0');\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",
       "            'z-index: 2;'\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(\n",
       "        'style',\n",
       "        'box-sizing: content-box;' +\n",
       "            'pointer-events: none;' +\n",
       "            'position: relative;' +\n",
       "            'z-index: 0;'\n",
       "    );\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;' +\n",
       "            'left: 0;' +\n",
       "            'pointer-events: none;' +\n",
       "            'position: absolute;' +\n",
       "            'top: 0;' +\n",
       "            '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",
       "            /* This rescales the canvas back to display pixels, so that it\n",
       "             * appears correct on HiDPI screens. */\n",
       "            canvas.style.width = width + 'px';\n",
       "            canvas.style.height = height + 'px';\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",
       "        /* User Agent sniffing is bad, but WebKit is busted:\n",
       "         * https://bugs.webkit.org/show_bug.cgi?id=144526\n",
       "         * https://bugs.webkit.org/show_bug.cgi?id=181818\n",
       "         * The worst that happens here is that they get an extra browser\n",
       "         * selection when dragging, if this check fails to catch them.\n",
       "         */\n",
       "        var UA = navigator.userAgent;\n",
       "        var isWebKit = /AppleWebKit/.test(UA) && !/Chrome/.test(UA);\n",
       "        if(isWebKit) {\n",
       "            return function (event) {\n",
       "                /* This prevents the web browser from automatically changing to\n",
       "                 * the text insertion cursor when the button is pressed. We\n",
       "                 * want to control all of the cursor setting manually through\n",
       "                 * the 'cursor' event from matplotlib */\n",
       "                event.preventDefault()\n",
       "                return fig.mouse_event(event, name);\n",
       "            };\n",
       "        } else {\n",
       "            return function (event) {\n",
       "                return fig.mouse_event(event, name);\n",
       "            };\n",
       "        }\n",
       "    }\n",
       "\n",
       "    canvas_div.addEventListener(\n",
       "        'mousedown',\n",
       "        on_mouse_event_closure('button_press')\n",
       "    );\n",
       "    canvas_div.addEventListener(\n",
       "        'mouseup',\n",
       "        on_mouse_event_closure('button_release')\n",
       "    );\n",
       "    canvas_div.addEventListener(\n",
       "        'dblclick',\n",
       "        on_mouse_event_closure('dblclick')\n",
       "    );\n",
       "    // Throttle sequential mouse events to 1 every 20ms.\n",
       "    canvas_div.addEventListener(\n",
       "        'mousemove',\n",
       "        on_mouse_event_closure('motion_notify')\n",
       "    );\n",
       "\n",
       "    canvas_div.addEventListener(\n",
       "        'mouseenter',\n",
       "        on_mouse_event_closure('figure_enter')\n",
       "    );\n",
       "    canvas_div.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",
       "    canvas_div.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.canvas_div.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",
       "\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",
       "    if (name === 'button_press') {\n",
       "        this.canvas.focus();\n",
       "        this.canvas_div.focus();\n",
       "    }\n",
       "\n",
       "    // from https://stackoverflow.com/q/1114465\n",
       "    var boundingRect = this.canvas.getBoundingClientRect();\n",
       "    var x = (event.clientX - boundingRect.left) * this.ratio;\n",
       "    var y = (event.clientY - boundingRect.top) * 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",
       "    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\", \"home\"], [\"Back\", \"Back to previous view\", \"fa fa-arrow-left\", \"back\"], [\"Forward\", \"Forward to next view\", \"fa fa-arrow-right\", \"forward\"], [\"\", \"\", \"\", \"\"], [\"Pan\", \"Left button pans, Right button zooms\\nx/y fixes axis, CTRL fixes aspect\", \"fa fa-arrows\", \"pan\"], [\"Zoom\", \"Zoom to rectangle\\nx/y fixes axis\", \"fa fa-square-o\", \"zoom\"], [\"\", \"\", \"\", \"\"], [\"Download\", \"Download plot\", \"fa fa-floppy-o\", \"download\"]];\n",
       "\n",
       "mpl.extensions = [\"eps\", \"jpeg\", \"pgf\", \"pdf\", \"png\", \"ps\", \"raw\", \"svg\", \"tif\", \"webp\"];\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,iVBORw0KGgoAAAANSUhEUgAAAoAAAAHgCAYAAAA10dzkAAAAAXNSR0IArs4c6QAAIABJREFUeF7s3Qd4FOXaxvE7lYSS0HsvCYTeBRXUgwULiNjwHDQoYjmK2MVjAQtYDypYAKUq2FGwokdBRem99xo6JISSvt81g+Gjs5s3ye7O/Pe6vD6/wzwz8/7mXt7H3Z13Qjwej0e8EEAAAQQQQAABBFwjEEID6JprzUARQAABBBBAAAFbgAaQICCAAAIIIIAAAi4ToAF02QVnuAgggAACCCCAAA0gGUAAAQQQQAABBFwmQAPosgvOcBFAAAEEEEAAARpAMoAAAggggAACCLhMgAbQZRec4SKAAAIIIIAAAjSAZAABBBBAAAEEEHCZAA2gyy44w0UAAQQQQAABBGgAyQACCCCAAAIIIOAyARpAl11whosAAggggAACCNAAkgEEEEAAAQQQQMBlAjSALrvgDBcBBBBAAAEEEKABJAMIIIAAAggggIDLBGgAXXbBGS4CCCCAAAIIIEADSAYQQAABBBBAAAGXCdAAuuyCM1wEEEAAAQQQQIAGkAwggAACCCCAAAIuE6ABdNkFZ7gIIIAAAggggAANIBlAAAEEEEAAAQRcJkAD6LILznARQAABBBBAAAEaQDKAAAIIIIAAAgi4TIAG0GUXnOEigAACCCCAAAI0gGQAAQQQQAABBBBwmQANoMsuOMNFAAEEEEAAAQRoAMkAAggggAACCCDgMgEaQJddcIaLAAIIIIAAAgjQAJIBBBBAAAEEEEDAZQI0gC674AwXAQQQQAABBBCgASQDCCCAAAIIIICAywRoAF12wRkuAggggAACCCBAA0gGEEAAAQQQQAABlwnQALrsgjNcBBBAAAEEEECABpAMIIAAAggggAACLhOgAXTZBWe4CCCAAAIIIIAADSAZQAABBBBAAAEEXCZAA+iyC85wEUAAAQQQQAABGkAygAACCCCAAAIIuEyABtBlF5zhIoAAAggggAACNIBkAAEEEEAAAQQQcJkADaDLLjjDRQABBBBAAAEEaADJAAIIIIAAAggg4DIBGkCXXXCGiwACCCCAAAII0ACSAQQQQAABBBBAwGUCNIAuu+AMFwEEEEAAAQQQoAEkAwgggAACCCCAgMsEaABddsEZLgIIIIAAAgggQANIBhBAAAEEEEAAAZcJ0AC67IIzXAQQQAABBBBAgAaQDCCAAAIIIIAAAi4ToAF02QVnuAgggAACCCCAAA0gGUAAAQQQQAABBFwmQAPosgvOcBFAAAEEEEAAARpAMoAAAggggAACCLhMwFEN4DvvvKNXX31V27dvV8OGDfXGG2/owgsvPO0l/fLLL/Xuu+9q4cKFSk9Pt7cfMGCALr/8cpdFgOEigAACCCCAgNsEHNMAfvLJJ+rZs6esJvD888/X8OHD9f7772v58uWqXr36Kde1X79+qly5si6++GKVLFlSo0eP1muvvaZZs2apefPmbssB40UAAQQQQAABFwk4pgFs27atWrRoYX+ql/tq0KCBrr32Wg0ePNirS2p9CnjTTTfpmWee8Wr7nJwcJSUlqUSJEgoJCfGqho0QQAABBBBAwL8CHo9Hqamp9gdBoaGh/j0ZPx3dEQ1gRkaGihYtqs8++0zdunU7RvnAAw/YX/FOnz79nLxWM1ezZk099thjuu+++067vfVVsfVP7mvbtm1KSEg4577ZAAEEEEAAAQQCT2DLli2qWrVq4J1YIZyRIxpA61O4KlWqaMaMGWrfvv0xtkGDBmns2LFatWrVOSmt3w6+9NJLWrFihcqXL3/a7a3fCA4cOPCUP7MCFBMTc85jsAECCCCAAAII+F/gwIEDqlatmpKTkxUbG+v/E/LDGTiqAfzzzz/Vrl27Y4wvvviixo8fr5UrV56VduLEierdu7e+/vprderU6YzbnvwJYG6AUlJSaAD9EF4OiQACCCCAQF4ErPnbavzcPH87ogE0+QrYunmkV69e9tfHV111lU85IkA+cbExAggggAACASHA/C05ogG00mTdBNKyZUv7LuDcl/X7vK5du57xJhDrk7/bb79d1v+1bhbx9UWAfBVjewQQQAABBPwvwPztoAYwdxmY9957z/4aeMSIERo5cqSWLVumGjVqqH///rJu2hg3bpydPKvpu/XWW/Xmm2/quuuuO5bG6Ohor38PQID8/ybmDBBAAAEEEPBVgPnbQQ2gdfGtT/9eeeUVeyHoRo0aaciQIerQoYOdi8TERG3cuFHTpk2z//+LLrrotHcH33bbbRozZoxXWSJAXjGxEQIIIIAAAgElwPztsAawsNNFgApbnOMhgAACCCBgLsD8TQNolCICZMRHMQIIIIAAAn4RYP6mATQKHgEy4qMYAQQQQAABvwgwf9MAGgWPABnxUYwAAggggIBfBJi/aQCNgkeAjPgoRgABBBBAwC8CzN80gEbBI0BGfBQjgAACCCDgFwHmbxpAo+ARICM+ihFAAAEEEPCLAPM3DaBR8AiQER/FCCCAAAII+EWA+ZsG0Ch4BMiIj2IEEEAAAQROK5CVnaOflu/U6D836qmrGqhJ1ZL5KsX8TQNoFCgCZMRHMQIIIIAAAicI7D+UoY/nbNH4vzYqKSXN/rPrWlTRf29slq9SzN80gEaBIkBGfBQjgAACCCBgC6zemarRMzZo0oJtSsvMsf+3MsUidUvb6vpn2xqqGBuVr1LM3zSARoEiQEZ8FCOAAAIIuFggJ8ej6Wt2a9QfG/T7mj3HJBIqxajX+TV1TdPKiooIKxAh5m8aQKNgESAjPooRQAABBFwokJaZrS/mb7Ubv3W7D9kCoSHSZQkVdfsFtdS6ZimFhIQUqAzzNw2gUcAIkBEfxQgggAACLhLYczBd4/7apA9nbtK+Qxn2yIsXCddNraspsX1NVStdtNA0mL9pAI3CRoCM+ChGAAEEEHCBwNpdB/X+7+v15YJtysg6+vu+KiWj7U/7bmxVVSWiIgpdgfmbBtAodATIiI9iBBBAAAGHCng8Hs3dtF/Dp6/Xzyt2Hhtls2oldeeFtXV5wwoKDwv12+iZv2kAjcJHgIz4KEYAAQQQcJiAdWPH1OU7Nfy3dVqwOdkenfVzvksbVFCfDrXVskbB/77PG1LmbxpAb3Jyxm0IkBEfxQgggAACDhFIz8rWVwu2afhv67X+7xs7IsND1b1FVfW+sJbqlCseUCNl/qYBNAokATLioxgBBBBAIMgFDqZnacKsTfrgjw3aeSDdHk1MVLh6tquhxPa1VK5EkYAcIfM3DaBRMAmQER/FCCCAAAJBKmA9scN6TNvYPzcq5UimPYoKMUXU+4La6tG2un13byC/mL9pAI3ySYCM+ChGAAEEEAgygR0paRr5+3pNmLVZRzKz7bOvXbaY7u5YR12bV1aR8IJZuDm/mZi/aQCNMkWAjPgoRgABBBAIEoEt+w7r3enr9PncrcrIPrqUS8PKMfr3xXV1ecOKCrNWcg6iF/M3DaBRXAmQER/FCCCAAAIBLmCt4ffOtLX6emGSsnM89tm2qVla/76krjrUK1vgT+woKB7mbxpAo2wRICM+ihFAAAEEAlRg1Y5UDf1ljb5dsl2eo32fLqxXVvddXFdta5cJ0LP2/rSYv2kAvU/LabYkQEZ8FCOAAAIIBJjA8qQDduP3/dIdx86sU4MKuv+SumparWSAnW3eT4f5mwYw7+mRRICM+ChGAAEEEAgQgaXbUvTW/9bYizjnvq5sXFH3XVxPCZVjAuQs8+80mL9pAI3SRICM+ChGAAEEEPCzwMmNn/XUjqubVLa/6o2vWMLPZ1dwh2f+pgE0ShcBMuKjGAEEEEDATwLLklL05s///4mf1fh1aVrZ/qq3bnnnNn653MzfNIBGbz0CZMRHMQIIIIBAIQtYN3e88fPqY7/xc1vjRwP4/4EL8Xhy7+8p5BQ64HA0gA64iAwBAQQQcIGAtZzLm/9bo28WJ9l39eZ+1fvAP9zxid/Jl5j5m08Ajd72BMiIj2IEEEAAgQIW2Lz3sN7432p9tWCb/l7GT9bNHf06xSmugvO/6j0TL/M3DaDRW48AGfFRjAACCCBQQALbU45o6C9r9emcLcr6u/O7NKGCHuwU58i7en1lZP6mAfQ1MydsT4CM+ChGAAEEEMhngb0H0/XOtHUaP3OTMrKOPrKtQ1w5PXxpnKPW8TNlY/6mATTKEAEy4qMYAQQQQCCfBFLTMvX+7xv0/u/rdSgj296r9ci2hy+Lc8STO/KJ6dhumL9pAI0yRYCM+ChGAAEEEDAUSMvM1vi/NtnP691/ONPeW6MqMXrksnh1jCsXtM/qNWQ5ZznzNw3gOUNytg0IkBEfxQgggAACeRTIzvHoi/lb9cZPq5WUkmbvpXa5Ynr40nh1blRRoaEhedyzO8qYv2kAjZJOgIz4KEYAAQQQ8FHAWrnt5xW79OqPK7V650G7ulJslPp1qqfuLaoqPCzUxz26c3PmbxpAo+QTICM+ihFAAAEEfBCYt2m/Bn+3QnM37berYqMj9O+L6+jWdjUVFRHmw57YlPmbBtDoXUCAjPgoRgABBBDwQmDd7oN69YdV+mHZDnvrqIhQ3X5+Ld3VsY7dBPLyXYD5mwbQ99QcV0GAjPgoRgABBBA4i8Du1HS9+b/Vmjh7i6zf/Fk/67uxVTV7EeeKsVHYGQgwf9MAGsRHIkBGfBQjgAACCJxG4EhGtr2cy3vT1x1b0qVTg/J6/Ir6qufip3fkZ1iYv2kAjfJEgIz4KEYAAQQQOE4g987e16eu0s4D6fafNK0aq/5XNtB5tctglY8CzN80gEZxIkBGfBQjgAACCPwt8MeaPXrh2+VauSPV/l+qlorWY1fU19WNK7GkSwGkhPmbBtAoVgTIiI9iBBBAwPUCa3elatB3K/XLyl22RUxUuO6/pJ5ubV9DRcK5s7egAsL8TQNolC0CZMRHMQIIIOBagX2HMvTGz6v10azN9g0e4aEh6tmuhvpeUk+likW61qWwBs78TQNolDUCZMRHMQIIIOA6gYysHI37a6Pe/N8apaZl2eO/NKGC+neur9rlirvOw18DZv6mATTKHgEy4qMYAQQQcI2A9QSPn5bv1KDvVmjj3sP2uBMqxeipqxuofZ2yrnEIlIEyf9MAGmWRABnxUYwAAgi4QmDljgN6/pvlmrF2rz3essWL6NHL43R9y2oK45m9fskA8zcNoFHwCJARH8UIIICAowWs3/kN+cn6nd8m5XikyPBQ9b6glu69uK6KFwl39NgDfXDM3zSARhklQEZ8FCOAAAKOFMjMztGHMzfZzd+Bv3/nd2XjiurfuYGqlS7qyDEH26CYv2kAjTJLgIz4KEYAAQQcJ/D7mt16bspyrdl10B5bg0oxevaaBBZyDrArzfxNA2gUSQJkxEcxAggg4BiBzXsP2ws5T12+0x5T6WKReuSyeN3Umt/5BeJFZv6mATTKJQEy4qMYAQQQCHqBwxlZenfaOg3/bb2sJV6smzpubVdD/TrFKTY6IujH59QBMH/TABplmwAZ8VGMAAIIBK2AtazLd0t26MVvlyspJc0ex/l1y+jZaxoqrkKJoB2XW06c+ZsG0CjrBMiIj2IEEEAgKAVW70zVgMnL9Oe6o8u6VCkZraevbqDLG1ZUSEhIUI7JbSfN/E0DaJR5AmTERzECCCAQVAKpaZl64+c1GvPnRvvxbdayLvd0rKO7O9ZRdCTP7Q2mi8n8TQNolFcCZMRHMQIIIBAUAtbXvV8vTNKL363Q7tR0+5ytx7c9fVWCqpdhWZeguIgnnSTzt8MawHfeeUevvvqqtm/froYNG+qNN97QhRdeeMZsTp8+XQ899JCWLVumypUr67HHHtPdd9/tdZYJkNdUbIgAAggEpcCqHal6+uulmr1hn33+NcsU1bNdGuri+PJBOR5O+qgA87eDGsBPPvlEPXv2lNUEnn/++Ro+fLjef/99LV++XNWrVz8l8xs2bFCjRo1055136q677tKMGTN07733auLEierevbtX7xEC5BUTGyGAAAJBJ3AwPUtv/rxao2Yc/bo3KiJU911cV3d2qK0i4XzdG3QXlE8AT7lkIR7rs20HvNq2basWLVro3XffPTaaBg0a6Nprr9XgwYNPGeHjjz+uyZMna8WKFcf+zPr0b9GiRfrrr7+8EqEB9IqJjRBAAIGgEbCmxG+XbLef3bvzwNGve69oWFFPXd1AVUvxdW/QXMhznCjzt0M+AczIyFDRokX12WefqVu3bscu+wMPPKCFCxfK+qr35FeHDh3UvHlzvfnmm8f+aNKkSbrxxht1+PBhRUSce/0mAuSUvwoYBwIIICCt331Qz05ept/X7LE5apQpqoFdGuoivu51XDyYvx3SACYlJalKlSr217jt27c/FtRBgwZp7NixWrVq1SnhjYuLU2Jiop588sljf/bnn3/aXx9b+6tUqdIpNenp6bL+yX1ZAapWrZpSUlIUExPjuDcIA0IAAQTcIJCWma13fl2r96avV0Z2jn13770XHb27NyqCr3udmAEaQIc1gFYD165du2NZffHFFzV+/HitXLnytA1gr1691L9//2N/ZjWQF1xwgX0TScWKFU+pGTBggAYOHHjK/04D6MS/HhgTAgi4QWDaql165utl2rzvsD3cjnHl9FzXhqpRppgbhu/aMdIAOqQBLKyvgPkE0LV/VzBwBBBwmMCOlDQ9980y+2ke1qtiTJSevSZBVzRiMWeHXerTDocG0CENoHV1rZtAWrZsad8FnPtKSEhQ165dz3gTyJQpU+y7hHNf99xzj/2bQW4CccPbnzEigIAbBaw7esf9tVGvT10t605f69m9ie1r6sFL41S8SLgbSVw5ZhpABzWAucvAvPfee/bXwCNGjNDIkSPtNf5q1Khhf9W7bds2jRs3zg577jIw1hIw1lIwVtNn3QXMMjCu/LuAQSOAgAsElm5L0ZOTlmjx1hR7tM2qldSL3RqpYeVYF4yeIR4vQAPooAbQurDWp3+vvPKK/Rs+a42/IUOGyLrb13pZN3xs3LhR06ZNO5YB6+7gBx988NhC0NbSMCwEzV8SCCCAgLMErE/6/jt1tcb8uUE5HqlEVLgev6K+bmlTXaGhPLvXWVfbu9HQADqsAfTusuffVgQo/yzZEwIIIFAQAj8v36lnvl6qpJQ0e/ddmla21/QrXyKqIA7HPoNEgPmbBtAoqgTIiI9iBBBAoMAEdh5I04DJy/T90qM3eVQrHa0Xrm1s3+XLCwHmbxpAo3cBATLioxgBBBDId4GcHI8mzN6sl79fqdS/b/LofWEt9ftHnKIjWdMv38GDdIfM3zSARtElQEZ8FCOAAAL5KrB2V6qe+GKJ5m7ab++3abWSGtytsRIqs1B/vkI7YGfM3zSARjEmQEZ8FCOAAAL5IpCela13p63TO7+us5/kUTQyTI9eHq9b29W0l3nhhcDJAszfNIBG7woCZMRHMQIIIGAsMG/Tfj3xxWKt2XXQ3tcl9cvr+WsbqUrJaON9swPnCjB/0wAapZsAGfFRjAACCORZ4FB6ll79cZXG/rVRHo9Utniknr2moa5uUkkhIXzql2dYlxQyf9MAGkWdABnxUYwAAgjkSWD66t168ssl2pZ8xK7v3qKqnrqqgUoVi8zT/ihynwDzNw2gUeoJkBEfxQgggIBPAsmHM/T8Nyv0xfytdl3VUtEa1K2xOrC0i0+ObCwxf9MAGr0PCJARH8UIIICA1wI/LN2up75apj0H02V9w2s9v/eRy+JVjOf3em3Ihv8vwPxNA2j0fiBARnwUI4AAAucU2J2abi/o/O2S7fa2dcsX18vdm6hljVLnrGUDBM4kwPxNA2j07iBARnwUI4AAAmcU8Hg8mrwoyW7+9h/OtJdzuadjHd13SV1FRbCgM9ExE2D+pgE0ShABMuKjGAEEEDitwK4DafrPV0v10/Kd9p8nVIrRK9c3UaMqsYghkC8CzN80gEZBIkBGfBQjgAACJwhYn/pNWrBNA6csV8qRTEWEhajvJfV090V1FBEWihYC+SbA/E0DaBQmAmTERzECCCBwTMD61O/JSUv084pd9v/WqEqMXruhqepX5DFuxCT/BZi/aQCNUkWAjPgoRgABBGR96vf1wiQ9O3nZsU/9+nWKU58OtfnUj3wUmADzNw2gUbgIkBEfxQgg4HIB6w7f/0xaoql//9bP+tTv9RuaKb5iCZfLMPyCFmD+pgE0yhgBMuKjGAEEXCzw7WJrXb8l9h2+/NbPxUHw09CZv2kAjaJHgIz4KEYAARcK7D+UoWcmL9OURUn26BtUsj71a6qEyvzWz4Vx8NuQmb9pAI3CR4CM+ChGAAGXCfxvxU498eUSWV/9Wuv6/fsia12/eooM5w5fl0XB78Nl/qYBNAohATLioxgBBFwikJqWqRe+WaFP5m6xR2w9zcP61K9ptZIuEWCYgSbA/E0DaJRJAmTERzECCLhAYNb6vXr4s0Xauv+I/Qzf3hfU0sOXxfM0Dxdc+0AeIvM3DaBRPgmQER/FCCDgYIG0zGz996fVGvn7enk8UtVS0fanfm1rl3HwqBlasAgwf9MAGmWVABnxUYwAAg4VWJ50QA9+slCrdqbaI7ypVTU9dXUDlYiKcOiIGVawCTB/0wAaZZYAGfFRjAACDhPIzvHYn/i9PnWVMrM9Kls8Ui9d10SdEio4bKQMJ9gFmL9pAI0yTICM+ChGAAEHCWzZd9j+rd/sDfvsUV2aUEEvXddYZYoXcdAoGYpTBJi/aQCNskyAjPgoRgABBwhYj3L7cv42+1FuB9OzVCwyTM9e01A3tKqqEOuuD14IBKAA8zcNoFEsCZARH8UIIBDkAsmHM/Sfr5bKeqqH9WpZo5SG3NhM1csUDfKRcfpOF2D+pgE0yjgBMuKjGAEEgljgz7V79NCni7TjQJrCQ0PUr1M93d2xjsLDWNQ5iC+ra06d+ZsG0CjsBMiIj2IEEAhCgfSsbL324yqN/H2Dffa1yxbTkJuasahzEF5LN58y8zcNoFH+CZARH8UIIBBkAmt2pqrvxwu1YvsB+8z/2ba6/nNVAxWNDA+ykXC6bhdg/qYBNHoPECAjPooRQCBIBKwbPT6cuUkvfLtC6Vk5Kl0sUi93b2Lf6csLgWAUYP6mATTKLQEy4qMYAQSCQGDvwXQ9/sVi/bxil322HeLK6bUbmqh8iaggOHtOEYHTCzB/0wAavTcIkBEfxQggEOACv63eba/ttzs1XZFhoXqic30ltq+p0FCWdwnwS8fpnUOA+ZsG0OhNQoCM+ChGAIEAFcjIytFrU1dpxG/r7TOsV7643urRXA0qxQToGXNaCPgmwPxNA+hbYk7amgAZ8VGMAAIBKLB+90H1/XiBlm47eqPHv86rrqeuSlBURFgAni2nhEDeBJi/aQDzlpy/qwiQER/FCCAQQALWjR6fz9tqP9HjcEa2ShaN0Cvdm+iyhhUD6Cw5FQTyR4D5mwbQKEkEyIiPYgQQCBCB1LRM/WfSUk1elGSfUbvaZey1/SrGcqNHgFwiTiOfBZi/aQCNIkWAjPgoRgCBABBYuCVZfScu0OZ9hxUWGqKHLo2zn+hh/TsvBJwqwPxNA2iUbQJkxEcxAgj4USAnx6ORv6/Xqz+uUlaOR1VKRuutHs3UskZpP54Vh0agcASYv2kAjZJGgIz4KEYAAT8J7DmYroc/XaTpq3fbZ3Bl44oafF0TxUZH+OmMOCwChSvA/E0DaJQ4AmTERzECCPhB4M91e9Tv44XalZquIuGhevaahurRpppCQvjK1w+Xg0P6SYD5mwbQKHoEyIiPYgQQKESB7ByP3vzfGg39ZY08Hqlu+eIadktz1a/I2n6FeBk4VIAIMH/TABpFkQAZ8VGMAAKFJLDzQJp9o8esDfvsI97YqqoGdGmoopHhhXQGHAaBwBJg/qYBNEokATLioxgBBApBwHqc24OfLNTeQxkqFhmmF7s11rXNqxTCkTkEAoErwPxNA2iUTgJkxEcxAggUoEBWdo6G/Lxa70xbZ3/laz3G7e1bmqt2ueIFeFR2jUBwCDB/0wAaJZUAGfFRjAACBSSwI+XoV76zNx79yvefbavr6at5nFsBcbPbIBRg/qYBNIotATLioxgBBApA4PivfIsXCddL3Rvr6iaVC+BI7BKB4BVg/qYBNEovATLioxgBBPJRwL7L9+fVGvrrWvsr34RKMXrnny1Us2yxfDwKu0LAGQLM3zSARkkmQEZ8FCOAQD4J7EpN0wMTF+qv9XvtPd7Strqe4SvffNJlN04UYP6mATTKNQEy4qMYAQTyQWDm+r26f+IC7U5NV9HIMA2+rrG6NuMu33ygZRcOFmD+pgE0ijcBMuKjGAEEDASsZ/m+99s6vfbjKuV4pPgKJfT2P1vYCzzzQgCBswswf9MAGr1HCJARH8UIIJBHgeTDGfazfP+3cpe9h+4tquqFaxspOjIsj3ukDAF3CTB/0wAaJZ4AGfFRjAACeRBYvDVZ93w4X9uSjygyPFTPd22oG1vxLN88UFLiYgHmbxpAo/gTICM+ihFAwAcBj8ejCbM3a+Dk5crIzlGNMkXtu3wbVo71YS9sigAClgDzNw2g0TuBABnxUYwAAl4KHMnI1n++WqIv52+zKy5NqKDXbmiq2OgIL/fAZgggcLwA8zcNoNE7ggAZ8VGMAAJeCGzYc0j3fDhPK3ekKjREeuyK+rqrQ22FhIR4Uc0mCCBwOgHmbxpAo3cGATLioxgBBM4hMHXZDvtmj9T0LJUtXkRDezRXuzplcEMAAUMB5m+HNID79+9X3759NXnyZDsSXbp00dChQ1WyZMnTRiQzM1NPPfWUvvvuO61fv16xsbHq1KmTXnrpJVWu7P0jkwiQ4TuQcgQQOK2A9VSP16eu0jvT1tl/3qpGKXuJlwoxUYghgEA+CDB/O6QB7Ny5s7Zu3aoRI0bYsejTp49q1qypKVOmnDYmKSkpuv7663XnnXeqadOmshrIfv36KSsrS3PnzvU6WgTIayo2RAABLwX2Hkxzrqy0AAAgAElEQVTXAx8v1B9r99gVt59fS/2vrK+IsFAv98BmCCBwLgHmbwc0gCtWrFBCQoJmzpyptm3b2tfc+vd27dpp5cqVio+PP1cO7D+fM2eO2rRpo02bNql69epe1RAgr5jYCAEEvBRYtMVa4mWeklLSFB0Rppevb6IuTb3/VsLLw7AZAq4XYP52QAM4atQoPfTQQ0pOTj4h0NbXv0OGDFGvXr28CvrPP/+syy67zN5PTEzMaWvS09Nl/ZP7sgJUrVo1WZ8onqnGq4OzEQIIuF7g49mb9czXy+wlXmqXLab3erZUXIUSrncBAIGCEKABdEADOGjQII0ZM0arV68+ISNxcXF289e/f/9zZictLU0XXHCB6tevrw8//PCM2w8YMEADBw485c9pAM9JzAYIIHAGgfSsbA2YvEwTZ2+xt7CWeHn9xqaKiWKJF0KDQEEJ0AAGcAN4pmbr+DBYX9tOnTpVY8eO1apVq07ISb169XTHHXfoiSeeOGt+rBtCbrjhBm3evFnTpk076yd5fAJYUG9F9ouAOwWSko/YX/ku2poia1WXRy6L1z0d6yjUWu+FFwIIFJgADWAAN4B79uyR9c/ZXtaNHhMmTMjzV8BW83fjjTfadwL/8ssvKlPGt+UVCFCBvTfZMQKOF/hr3V7dN2G+9h7KsBd0fqtHc3WMK+f4cTNABAJBgPk7gBtAbwOSexPIrFmz7Js4rJf17+edd95ZbwLJbf7WrFmjX3/9VeXK+f4XLwHy9iqxHQII5ApYj3QbPWOjXvxuhazlXhIqxWh4z5aqVrooSAggUEgCzN8OaACtrFjLwCQlJWn48OF2dKxlYGrUqHHCMjDW7/sGDx6sbt262cu9dO/eXfPnz9c333yjChUqHItc6dKlFRkZ6VUECZBXTGyEAAJ/C1iPdOv/5WJ9tTDJ/l+ubVZZg69roujIMIwQQKAQBZi/HdIA7tu375SFoIcNG3bCQtDWY5NGjx6txMREbdy4UbVq1Tpt1KxPAy+66CKvYkiAvGJiIwQQkLRl32HdNX6elm8/oLDQEP3nygbqdX5NHulGOhDwgwDzt0MaQD9kxz4kAfKXPMdFILgE/ly7R/+eMF/7D2eqTLFIDbulBY90C65LyNk6TID5mwbQKNIEyIiPYgQcL2D93m/UjI0a9Pfv/ZpUjdV7/2qpyiWjHT92BohAIAswf9MAGuWTABnxUYyAowXSMq3f+y3RpAXb7HF2b1FVL3ZrpKgIfu/n6AvP4IJCgPmbBtAoqATIiI9iBBwrsC35iO4eP09LtqXYv/d76qoGSmzP7/0ce8EZWNAJMH/TABqFlgAZ8VGMgCMFZm/Yp3s/mqc9BzNUulik3ub3fo68zgwquAWYv2kAjRJMgIz4KEbAcQIfztxkP9Yt6+/1/Ubc2lJVS7G+n+MuNAMKegHmbxpAoxATICM+ihFwjEBGVo6etZ/nu9ke09VNKunV65uyvp9jrjADcZoA8zcNoFGmCZARH8UIOEJgz8F0+3m+czbut5/n+9jl9XV3x9qs7+eIq8sgnCrA/E0DaJRtAmTERzECQS+wdFuK+oybq6SUNJWICtdbNzfXxfXLB/24GAACThdg/qYBNMo4ATLioxiBoBb4ZnGSHvlskdIyc1S7bDGNvK2V6pQrHtRj4uQRcIsA8zcNoFHWCZARH8UIBKVATo5H//1ptYb9utY+/45x5fRWj+aKjY4IyvFw0gi4UYD5mwbQKPcEyIiPYgSCTuBgepYe/GShflq+0z73Ph1q6/Er6ttr/fFCAIHgEWD+pgE0SisBMuKjGIGgEtiy77DuHDdXK3ekKjI8VC9d11jXtagaVGPgZBFA4KgA8zcNoNF7gQAZ8VGMQNAIzFq/V/d8NF/7DmWoXIkiGtGzpZpXLxU058+JIoDAiQLM3zSARu8JAmTERzECQSHwyZzNeuqrpcrM9qhRlRiNvLWVKsVGB8W5c5IIIHB6AeZvGkCj9wYBMuKjGIGAFsjO8ejFb1do1IwN9nle1aSSXmNx54C+ZpwcAt4KMH/TAHqbldNuR4CM+ChGIGAFUtMydf/EBZq2ard9jg92ilPff9RlceeAvWKcGAK+CTB/0wD6lpiTtiZARnwUIxCQApv3HtYdY+doza6DiooI1es3NLM//eOFAALOEWD+pgE0SjMBMuKjGIGAE5i9YZ/uGj9X+w9nqkJMEb1/a2s1rhobcOfJCSGAgJkA8zcNoFGCCJARH8UIBJTA5/O2qv+Xi+2bPZpUjbVv9qgQExVQ58jJIIBA/ggwf9MAGiWJABnxUYxAQAhYT/Z4deoqvTttnX0+VzauaH/tGx0ZFhDnx0kggED+CzB/0wAapYoAGfFRjIDfBQ5nZOmhTxbph2U77HO5/5K69g0foTzZw+/XhhNAoCAFmL9pAI3yRYCM+ChGwK8COw+kqffYuVqyLUWRYaF6+frG6tacJ3v49aJwcAQKSYD5mwbQKGoEyIiPYgT8JrAsKUV3jJmrHQfSVLpYpP1kj1Y1S/vtfDgwAggUrgDzNw2gUeIIkBEfxQj4ReDn5TvV9+MFOpyRrbrli2vUba1VvUxRv5wLB0UAAf8IMH/TABoljwAZ8VGMQKEKeDweffDHBr343Qp5PNIFdcvq7X+2UGx0RKGeBwdDAAH/CzB/0wAapZAAGfFRjEChCWRl52jAlGX6cOZm+5g92lTXc10bKiIstNDOgQMhgEDgCDB/0wAapZEAGfFRjEChCFiPdfv3hAX6bfVuhYRIT3ZuoN4X1uKxboWiz0EQCEwB5m8aQKNkEiAjPooRKHCBrfsP2zd7rNqZquiIML1xczNd3rBigR+XAyCAQGALMH/TABollAAZ8VGMQIEKLNySbC/zsudgusqVKGLf7MFj3QqUnJ0jEDQCzN80gEZhJUBGfBQjUGACPyzdrn6fLFRaZo7qVyyhUYmtVblkdIEdjx0jgEBwCTB/0wAaJZYAGfFRjEC+C1h3+o78fb0Gf7/SvtP3ovhyGnZLCxUvEp7vx2KHCCAQvALM3zSARuklQEZ8FCOQrwKZ2Tl65utlmjj76J2+Pc+roWevSVA4d/rmqzM7Q8AJAszfNIBGOSZARnwUI5BvAtadvvd+NF+/r9lj3+n71FUJuv38mtzpm2/C7AgBZwkwf9MAGiWaABnxUYxAvghsSz6i20fPOXan75s3N9Nl3OmbL7bsBAGnCjB/0wAaZZsAGfFRjICxwNJtKbp9zBztSuVOX2NMdoCAiwSYv2kAjeJOgIz4KEbASOB/K3bq/olHn+kbX6GERvVqrSrc6WtkSjECbhFg/qYBNMo6ATLioxiBPAuM/XOjBk5ZphyPdGG9o8/0jYnimb55BqUQAZcJMH/TABpFngAZ8VGMgM8C2TkevfjtCo2ascGuvbl1NT1/bSOe6euzJAUIuFuA+ZsG0OgdQICM+ChGwCeBIxnZeuDjBZq6fKdd99gV8bqnYx3u9PVJkY0RQMASYP6mATR6JxAgIz6KEfBaYHdqunqPm6tFW5IVGR6q129oqmuaVva6ng0RQACB4wWYvwugAXzuuefOmrJnnnnGMSkkQI65lAwkgAXW7kpV4ug52rr/iEoWjdDIW1updc3SAXzGnBoCCAS6APN3ATSAzZs3P+G6Z2ZmasOGDQoPD1edOnU0f/78QM+F1+dHgLymYkME8iTw17q9umv8XB1Iy1LNMkU1ulcb1SpbLE/7oggBBBDIFWD+LoAG8HTxsqATExPVrVs39ezZ0zEJJECOuZQMJAAFJi3Yqsc+X6zMbI9a1ihlf/JXulhkAJ4pp4QAAsEmwPxdSA2gFYylS5fq6quv1saNG4MtJ2c8XwLkmEvJQAJIwOPxaNgva/X6T6vts7qqcSW9fmNTRUWEBdBZcioIIBDMAszfhdgA/vHHH7rmmmu0f//+YM7MCedOgBxzKRlIgAhkZufoP5OW6NO5W+0zuqtDbT1+RX2FhoYEyBlyGggg4AQB5u8CaADfeuutE7Jh/df89u3bNX78eHXo0EETJ050QnbsMRAgx1xKBhIAAgfSMvXvj+br9zV7ZPV7A7s2Us/zagTAmXEKCCDgNAHm7wJoAGvVqnVCTkJDQ1WuXDldcskl6t+/v0qUKOGYHBEgx1xKBuJngaTkI/YzfVfuSFXRyDANu6W5Lqlfwc9nxeERQMCpAszfBdAAOjUspxsXAXLT1WasBSWwLCnFbv52HkhXuRJFNDqxtRpViS2ow7FfBBBAgG/wRANo9DagATTioxgB/bpql+77aL4OZWQrrkJxe5mXKiWjkUEAAQQKVID5mwbQKGAEyIiPYpcLTJi1WU9/vVTW833Pr1tG7/6rpWKiIlyuwvARQKAwBJi/aQCNckaAjPgodqlATo5Hr01dpXemrbMFureoqsHXNbYf8cYLAQQQKAwB5m8aQKOcESAjPopdKJCela1HP1usyYuS7NH361RPD/yjnkJCWObFhXFgyAj4TYD5mwbQKHwEyIiPYpcJJB/OUJ/x8zR7wz6Fh4bope5NdH3Lqi5TYLgIIBAIAszfNIBGOSRARnwUu0hg897DShwzW+t3H1KJIuF6r2dLnV+3rIsEGCoCCASSAPM3DaBRHgmQER/FLhFYuCVZvcfO0Z6DGaocG6VRvVqrfsUYl4yeYSKAQCAKMH/TABrlkgAZ8VHsAoGpy3ao78cLlJaZo4RKMRrdq7UqxES5YOQMEQEEAlmA+ZsG0CifBMiIj2KHC4yZsUEDv1kuj0e6KL6cht3SQsWLhDt81AwPAQSCQYD52yEN4P79+9W3b19NnjzZzl2XLl00dOhQlSxZ0qsc3nXXXRoxYoSGDBmifv36eVVjbUSAvKZiQxcJWMu8vPjdCn3wxwZ71Le0ra7nujRUeBjLvLgoBgwVgYAWYP52SAPYuXNnbd261W7irFefPn1Us2ZNTZky5ZwB/OqrrzRgwADt3r1bjz76KA3gOcXYAIEzC6RlZqvfxwv1w7Id9kaPXRGvezrWYZkXQoMAAgElQAPogAZwxYoVSkhI0MyZM9W2bVs7YNa/t2vXTitXrlR8fPwZQ7dt2za75scff9RVV11lN398AhhQ71FOJogE9h5MV+9xc7Vgc7Iiw0L16g1N1LVZlSAaAaeKAAJuEaABdEADOGrUKD300ENKTk4+IbfW17/WV7q9evU6bZ5zcnLUqVMnde3aVQ888ID9ieG5GsD09HRZ/+S+rABVq1ZNKSkpionhrka3/MXBOE8V2LDnkBJHz9amvYcVGx2hET1bqm3tMlAhgAACASlAA+iABnDQoEEaM2aMVq9efULI4uLi7Oavf//+pw3f4MGD9euvv9qf/llPIfCmAbS+Kh44cOAp+6MBDMj3NydVSAJzN+7TnePmav/hTFUrHa3RiW1Ut3zxQjo6h0EAAQR8F6ABDOAG8EzN1vGXec6cOZo6darGjh2rVatWnZCAevXq6Y477tATTzxxSjLmzZtnf+U7f/58Va5c2f5zbxpAPgH0/U1GhbMFvl28XQ9+ulAZWTlqWjVW79/WWuVKFHH2oBkdAggEvQANYAA3gHv27JH1z9leVtM2YcIEn78CfuONN+ya0ND/vysxOzvb/v+tr3Q3btzoVbgJkFdMbORAAY/Ho5G/r9eg71bao7s0oYLevLmZikayzIsDLzdDQsBxAszfAdwAepu23JtAZs2apTZt2thl1r+fd955Z7wJZO/evdq+ffsJh7j88svVs2dP+2vjs904cnwRAfL2KrGdkwSysnP03DfLNe6vTfawEtvX1NNXJygsNMRJw2QsCCDgYAHmbwc0gFY+rWVgkpKSNHz4cDuu1jIwNWrUOGEZmPr168v63V+3bt1OG2lvvgI+uZAAOfhvB4Z2WoHDGVm6f8IC/W/lLoWESP+5soHuuKAWy7yQFwQQCCoB5m+HNID79u07ZSHoYcOGnbAQtHWjx+jRo5WYmEgDGFRvU042UAR2pabpjjFztWRbioqEh+qNm5qpc+NKgXJ6nAcCCCDgtQANoEMaQK+veD5vSIDyGZTdBazAmp2pShw9R9uSj6h0sUiNvLWVWtYoFbDny4khgAACZxNg/qYBNHqHECAjPoqDROCvdXvVZ/xcpaZlqVbZYhqd2Fo1yxYLkrPnNBFAAIFTBZi/aQCN3hcEyIiP4iAQmLRgqx77fLEysz1qVaOURtzayv4EkBcCCCAQzALM3zSARvklQEZ8FAewgLXMy9Bf1uq/Px1dYP2qxpX0+o1NFRURFsBnzakhgAAC3gkwf9MAepeUM2xFgIz4KA5QgczsHP1n0hJ9OnerfYZ3daitx6+or1CWeQnQK8ZpIYCArwLM3zSAvmbmhO0JkBEfxQEokJqWqXs/mq/f1+yR1e8N7NpIPc+rEYBnyikhgAACeRdg/qYBzHt6JBEgIz6KA0wgKfmIbh8zRyt3pKpoZJiG3dJcl9SvEGBnyekggAAC5gLM3zSARikiQEZ8FAeQwLKkFLv523kg3X6W76jbWqtx1dgAOkNOBQEEEMg/AeZvGkCjNBEgIz6KA0Tg11W7dN9H83UoI1txFYprdK82qlIyOkDOjtNAAAEE8l+A+ZsG0ChVBMiIj+IAEJgwa7Oe/nqpsnM8al+njN79V0vFRkcEwJlxCggggEDBCTB/0wAapYsAGfFR7EeBnByPXvlxld6bvs4+i+4tqmrwdY0VGR7qx7Pi0AgggEDhCDB/0wAaJY0AGfFR7CeBtMxsPfLZIn2zeLt9Bg92ilPff9SV9bxsXggggIAbBJi/aQCNck6AjPgo9oPA/kMZunPcXM3dtF/hoSF6uXsTdW9Z1Q9nwiERQAAB/wkwf9MAGqWPABnxUVzIAhv3HFKvMXO0Yc8hlYgK1/B/tVT7umUL+Sw4HAIIIOB/AeZvGkCjFBIgIz6KC1Fg3qb99id/+w5l2Hf4ju7VWnEVShTiGXAoBBBAIHAEmL9pAI3SSICM+CguJIHvlmxXv08WKiMrR42rxOqDxFYqXyKqkI7OYRBAAIHAE2D+pgE0SiUBMuKjuIAFPB6PRv6+XoO+W2kfqVOD8nqrR3MVjQwv4COzewQQQCCwBZi/aQCNEkqAjPgoLkCBrOwcDZiyTB/O3Gwf5bZ2NfTMNQ0VZj3glxcCCCDgcgHmbxpAo7cAATLio7iABA6lZ+n+iQv0y8pdslZ2+c+VDXTHBbVY5qWAvNktAggEnwDzNw2gUWoJkBEfxQUgsPNAmv1M32VJB1QkPFRv3txMVzSqVABHYpcIIIBA8Aowf9MAGqWXABnxUZzPAqt2pKrX6NlKSklTmWKRGnlbK7WoXiqfj8LuEEAAgeAXYP6mATRKMQEy4qM4HwVmrN2ju8fPU2p6lmqXLWYv81KjTLF8PAK7QgABBJwjwPxNA2iUZgJkxEdxPgl8Pm+rnvhisbJyPGpds5RG9GylUsUi82nv7AYBBBBwngDzNw2gUaoJkBEfxYYC1jIvb/y8Rm/+b429p2uaVtar1zdRVESY4Z4pRwABBJwtwPxNA2iUcAJkxEexgYC1qHP/L5foi/lb7b3ce1EdPXJZvEJZ5sVAlVIEEHCLAPM3DaBR1gmQER/FeRRIOZKpez+apxlr99rr+j3ftZFuaVs9j3ujDAEEEHCfAPM3DaBR6gmQER/FeRDYlnzEvtN39c6DKhYZpmH/bKGL48vnYU+UIIAAAu4VYP6mATRKPwEy4qPYR4Gl21LUa8wc7U5NV4WYIhqV2FoNK8f6uBc2RwABBBBg/qYBNHoXECAjPop9EPh15S79e8J8Hc7IVv2KJezmr3LJaB/2wKYIIIAAArkCzN80gEbvBgJkxEexlwIfzdqkp79aqhyPdGG9snr7ny0UExXhZTWbIYAAAgicLMD8TQNo9K4gQEZ8FJ9DICfHo5d/XKnh09fbW97QsqoGXddYEWGh2CGAAAIIGAgwf9MAGsRHIkBGfBSfRSAtM1uPfLZI3yzebm/18KVxuu+SugoJCcENAQQQQMBQgPmbBtAoQgTIiI/iMwjsP5ShPuPnas7G/YoIC9HL3ZvouhZV8UIAAQQQyCcB5m8aQKMoESAjPopPI7B572Eljp6t9XsOqURUuIb/q6Xa1y2LFQIIIIBAPgowf9MAGsWJABnxUXySwILN+9V77FztPZShKiWjNbpXa8VVKIETAggggEA+CzB/0wAaRYoAGfFRfJzAj8t26IGPFygtM0eNqsRo1G2tVT4mCiMEEEAAgQIQYP6mATSKFQEy4qP4b4FRf2zQ898ul8cjXRxfTsNuaaFiRcLxQQABBBAoIAHmbxpAo2gRICM+1xdn53j04rcrNGrGBtvCep7vc10aKpxlXlyfDQAQQKBgBZi/aQCNEkaAjPhcXWwt89Lv44X6YdkO2+GxK+J1T8c6LPPi6lQweAQQKCwB5m8aQKOsESAjPtcW7z2Yrt7j5mrB5mRFhoXqtRubqkvTyq71YOAIIIBAYQswf9MAGmWOABnxubJ4w55D9jIvm/YeVmx0hEb0bKm2tcu40oJBI4AAAv4SYP6mATTKHgEy4nNd8bxN++xlXvYfzlTVUtEa06uN6pYv7joHBowAAgj4W4D5mwbQKIMEyIjPVcXfL9muBz5ZqIysHDWpGqsPbmutciWKuMqAwSKAAAKBIsD8TQNolEUCZMTnmuL3f1+vF79bYS/z0qlBeb3Vo7mKRrLMi2sCwEARQCDgBJi/aQCNQkmAjPgcX2wt8/LCt8s1esZGe6w9z6uhAV0aKiw0xPFjZ4AIIIBAIAswf9MAGuWTABnxObrYWubFerLHj8t22uN88sr6uvPC2izz4uirzuAQQCBYBJi/aQCNskqAjPgcW3zyMi+v39hU17DMi2OvNwNDAIHgE2D+pgE0Si0BMuJzZPHGv5d52fj3Mi8jb22lNrVKO3KsDAoBBBAIVgHmbxpAo+wSICM+xxXP37zfXuZl36EMlnlx3NVlQAgg4CQB5m8aQKM8EyAjPkcVT122Q30/XqC0zBw1rhKrDxJbqXyJKEeNkcEggAACThFg/qYBNMoyATLic0zxuL826tnJy+xlXi6pX15DezRXsSIs8+KYC8xAEEDAcQLM3zSARqEmQEZ8QV+ck+PRyz+s1PDf1ttj6dGmup7v2lDhYaFBPzYGgAACCDhZgPmbBtAo3wTIiC+oi9OzsvXIZ4s1ZVGSPY5HL4/XvRfVYZmXoL6qnDwCCLhFgPmbBtAo6wTIiC9oi1MOZ+rO8XM1e8M+hYeG6JXrm+i6FlWDdjycOAIIIOA2AeZvGkCjzBMgI76gLN66/7ASR8/R2l0HVbxIuIb3bKnz65YNyrFw0ggggIBbBZi/aQCNsk+AjPiCrnhZUop6jZ6jXanpqhgTpdG9WqtBpZigGwcnjAACCLhdgPmbBtDoPUCAjPiCqvi31bt1z4fzdCgjW/EVSmjM7a1VKTY6qMbAySKAAAIIHBVg/nZIA7h//3717dtXkydPti9sly5dNHToUJUsWfKsWV+xYoUef/xxTZ8+XTk5OWrYsKE+/fRTVa9e3av3CAHyiinoN/p83lY98cViZeV41K52Gb3Xs6VioyOCflwMAAEEEHCrAPO3QxrAzp07a+vWrRoxYoSd5T59+qhmzZqaMmXKGbO9bt06tWnTRnfccYd69Oih2NhYWQ1h69atVb58ea/eEwTIK6ag3cjj8WjYL2v1+k+r7TF0bVbZvuGjSHhY0I6JE0cAAQQQ4BNAKwMhHmuWC+KX1bQlJCRo5syZatu2rT0S69/btWunlStXKj4+/rSju/nmmxUREaHx48fnefQ0gHmmC/jCrOwcPf31Uk2cvcU+13suqqNHL4tXaGhIwJ87J4gAAgggcHYB5m8HNICjRo3SQw89pOTk5BOutvX175AhQ9SrV69TUmB93Wt94vfYY4/pjz/+0IIFC1SrVi31799f1157rdfvGwLkNVVQbXg4I0v3TVigX1buUkiINLBLQ93armZQjYGTRQABBBA4swDztwMawEGDBmnMmDFavfro13S5r7i4OLv5s5q6k187duxQpUqVVLRoUb3wwgu6+OKL9cMPP+jJJ5/Ur7/+qo4dO542Nenp6bL+yX1ZAapWrZpSUlIUE8PdoE74y2bPwXTdMWaOFm1NUZHwUL3Vo7kub1jRCUNjDAgggAACfwvQAAZwAzhgwAANHDjwrGGdM2eOpk6dqrFjx2rVqlUnbFuvXj37931PPPHEKftISkpSlSpV7N/+TZgw4difWzePFCtWTBMnTjztcc90TjSAzvg7ZcOeQ7pt1Gxt3ndYpYpG6P3bWqtljVLOGByjQAABBBA44QMc65tAN8/fAfsbwD179sj652wv60YPq4Hz9SvgjIwMu9F79tln9dRTTx07hHVHsPWV8IwZM/gE0GV/USzYvF93jJ2rfYcyVK10tMb2aqPa5Yq7TIHhIoAAAu4Q4BPAAP4E0NsI5t4EMmvWLPuuXutl/ft555131ptA2rdvrzp16pxwE0i3bt0UHR19wqeCZzsPAuTtVQrs7X5avlP3T5yvtMwcNakaqw9ua61yJYoE9klzdggggAACeRZg/nZAA2hdfWsZGOtr3eHDh9thsJaBqVGjxgnLwNSvX1+DBw+W1eRZr0mTJummm27S22+/few3gP369dO0adN0wQUXeBUqAuQVU0Bv9NGsTXr6q6XK8UgXx5fTsFtaqFiR8IA+Z04OAQQQQMBMgPnbIQ3gvn37TlkIetiwYScsBB0SEqLRo0crMTHxWGqsO4itptBaQ9BaLsb6zWHXrl29ThUB8poq4Da0Vj96beoqvf3rOvvcbmpVTS92a6TwsNCAO1dOCAEEEEAgfwWYvx3SAOZvLLzfGwHy3iqQtszIytETXy7Wl/O32afVr1M9PfCPerL+I4EXAggggIDzBZi/aQCNUk6AjPj8Upyalql7P5qv39fsUVhoiAZ3a6wbW1fzy7lwUBz4Bk8AABtkSURBVAQQQAAB/wgwf9MAGiWPABnxFXrxzgNp6jV6jpZvP6CikWF6+58tdHG8d4/9K/ST5YAIIIAAAgUmwPxNA2gULgJkxFeoxWt3peq2UXO0LfmIyhaP1KjE1mpStWShngMHQwABBBAIDAHmbxpAoyQSICO+Qiues3Gfeo+dq5QjmapVtpi9xl/1MkUL7fgcCAEEEEAgsASYv2kAjRJJgIz4CqX4h6Xb1ffjhbJu/GhevaS9xl/pYpGFcmwOggACCCAQmALM3zSARskkQEZ8BV48ZsYGDfxmuTweqVODChrao7miI8MK/LgcAAEEEEAgsAWYv2kAjRJKgIz4Cqw4J8ejl39YqeG/rbeP8c+21TWwS0PW+CswcXaMAAIIBJcA8zcNoFFiCZARX4EUW1/1Pvr5In29MMne/6OXx+vei+qwxl+BaLNTBBBAIDgFmL9pAI2SS4CM+PK9+EBapu75cJ5mrN2r8NAQvdS9ia5vWTXfj8MOEUAAAQSCW4D5mwbQKMEEyIgvX4utNf5uGzVbK3ekqlhkmN79V0t1iCuXr8dgZwgggAACzhBg/qYBNEoyATLiy7fiNTtTlTg6d42/IhrTq7UaVYnNt/2zIwQQQAABZwkwf9MAGiWaABnx5Uvx7A3WGn9zdCAtS7XLHV3jr1pp1vjLF1x2ggACCDhUgPmbBtAo2gTIiM+4+Psl2/XAJ0fX+Gvx9xp/pVjjz9iVHSCAAAJOF2D+pgE0yjgBMuIzKh49Y4Oe+3uNv8sSKuitHs0VFcEaf0aoFCOAAAIuEWD+pgE0ijoBMuLLU7G1xt/g71do5O8b7Pqe59XQgC4NFRYakqf9UYQAAggg4D4B5m8aQKPUEyAjPp+L07Oy9fCni/TN4u127eNX1NfdHWuzxp/PkhQggAAC7hZg/qYBNHoHECAjPp+KUw5nqs/4uZq1YZ8iwkL0yvVN1K05a/z5hMjGCCCAAAK2APM3DaDRW4EAGfF5Xbwt+YgSR83Wml0HVaJIuN7r2VLn1y3rdT0bIoAAAgggcLwA8zcNoNE7ggAZ8XlVvHRbim4fM0e7UtNVIcZa46+NGlSK8aqWjRBAAAEEEDidAPM3DaDRO4MAGfGds/jXVbt030fzdSgjW/EVSmh0r9aqXDL6nHVsgAACCCCAwNkEmL9pAI3eIQTIiO+sxR/P3qz/fLVU2TkenV+3jP1ot5ioiII7IHtGAAEEEHCNAPM3DaBR2AmQEd9piz0ej/7702oN/WWt/efXtaiil65rosjw0Pw/GHtEAAEEEHClAPM3DaBR8AmQEd8pxdYyL499vlhfL0yy/6zvJXX14KVxLPOSv8zsDQEEEHC9APM3DaDRm4AAGfGdUJx8OEN9xs+T9Wzf8NAQDerWWDe2rpZ/B2BPCCCAAAII/C3A/E0DaPRmIEBGfMeKN+89rMQxs7V+9yF7mZd3/tVCF9Yrlz87Zy8IIIAAAgicJMD8TQNo9KYgQEZ8dvH8zft159i52nsoQ5VjozSqV2vVr8gyL+ay7AEBBBBA4EwCzN80gEbvDgJkxKdvF2/XQ58uVHpWjhpWjtGoxNaqEBNltlOqEUAAAQQQOIcA8zcNoNGbhADljc+60/e96ev18g8r7R38o355vdWjuYoVCc/bDqlCAAEEEEDABwHmbxpAH+Jy6qYEyHe+zOwcPTVpqT6Zu8UuTmxfU09fnaCw0BDfd0YFAggggAACeRBg/qYBzENs/r+EAPnGl3IkU//+aL7+WLtHVr9nNX69zq/l207YGgEEEEAAAUMB5m8aQKMIESDv+TbtPWQ/03fd7kMqGhmmt25urk4JFbzfAVsigAACCCCQTwLM3zSARlEiQN7xzd24z17jb9+hDFWMidIHia3UsHKsd8VshQACCCCAQD4LMH/TABpFigCdm++rBdvsp3tkZOeocZVYvX9bK+70PTcbWyCAAAIIFKAA8zcNoFG8CNCZ+XJyPHrj59V66+9n+l7esIKG3NRMRSO509codBQjgAACCBgLMH/TABqFiACdnu9IRrYe+WyRvl2y3d7gro619fjl9RXKnb5GeaMYAQQQQCB/BJi/aQCNkkSATuXbeSBNd46bq8VbUxQRdvSZvje04pm+RkGjGAEEEEAgXwWYv2kAjQJFgE7kW7otRb3HztWOA2kqVTRCw3u2UptapY2MKUYAAQQQQCC/BZi/aQCNMkWA/p/v+yXWY90W6UhmtuqWL65Rt7VW9TJFjXwpRgABBBBAoCAEmL9pAI1yRYAk67FuQ39Zq//+tNq27BBXTsNuaa6YqAgjW4oRQAABBBAoKAHmbxpAo2y5PUBpmUdv9vhm8dGbPW4/v5aevLK+wsNCjVwpRgABBBBAoCAF3D5/W7YhHusjHF55EnBzgI6/2SM8NETPX9tIPdpUz5MjRQgggAACCBSmgJvn71xnGkCDxLk1QAu3JKvPuLnalZpu3+zx7r9a6rzaZQwkKUUAAQQQQKDwBNw6fx8vTANokDc3Bujrhdv0qPVkj6wc1Stf3H6yR40yxQwUKUUAAQQQQKBwBdw4f58sTANokDk3Bch6sserU1fp3WnrbLF/1C+vN25uphLc7GGQIEoRQAABBPwh4Kb5+0y+NIAGyXNLgFLTMvXgJwv184pdttbdHevo0cvjFcaTPQzSQykCCCCAgL8E3DJ/n82XBtAgfW4I0Ka9h+zFndfsOqjI8FC93L2xujWvaqBGKQIIIIAAAv4VcMP8fS5hGsBzCZ3lz50eoBlr9+jej+Yr5UimypcoouE9W6p59VIGYpQigAACCCDgfwGnz9/eCNMAeqN0hm2cGiBrZaAxf27UC9+uUHaOR02rldSIni1VISbKQItSBBBAAAEEAkPAqfO3L7o0gL5onbStEwOUnpWtZ75apk/mbrFHe13zKhp0XWNFRYQZSFGKAAIIIIBA4Ag4cf72VZcG0Fex47Z3WoB2HUjT3R/O0/zNybLu7+jfuYF6X1hLISEhBkqUIoAAAgggEFgCTpu/86JLA5gXtb9rnBSgRVuSddf4edpxIE0xUeEaeksLdYwrZ6BDKQIIIIAAAoEp4KT5O6/CNIB5lZPklAB9OX+rnvhyib24c93yxTXy1laqVZbFnQ2iQSkCCCCAQAALOGX+NiGmATTQC/YAZWXn6KXvV+r9PzbYCp0alNeQm1jc2SASlCKAAAIIBIFAsM/f+UFMA2igGMwB2n8oQ/dNnK8Za/faAvddXFcPXRqnUBZ3NkgEpQgggAACwSAQzPN3fvnSABpIBmuAVmw/oDvHzdXW/UdUNDJMr9/QVJ0bVzKQoBQBBBBAAIHgEQjW+Ts/hWkADTSDMUDfLt6uRz5bpCOZ2apeuqhG3NpS9SvGGChQigACCCCAQHAJBOP8nd/CNIAGosEUIGtB59emrtK709bZI76wXlkN7dFcJYtGGghQigACCCCAQPAJBNP8XVC6jmgA9+/fr759+2ry5Mm2U5cuXTR06FCVLFnyjG4HDx7UE088oa+++kp79+5VzZo17X3cc889XlsHS4BSDmfq/o8X6LfVu+2x9elQW49dHq/wsFCvx8qGCCCAAAIIOEUgWObvgvR2RAPYuXNnbd26VSNGjDja4PTpYzd0U6ZMOaPdnXfeqV9//VXvv/++ve3UqVN177336osvvlDXrl29Mg+GAK3akao+4+dq097DiooI1SvXN1WXppW9Gh8bIYAAAggg4ESBYJi/C9o96BvAFStWKCEhQTNnzlTbtm1tL+vf27Vrp5UrVyo+Pv60ho0aNdJNN92kp59++tift2zZUldeeaWef/55r9wDPUDfLTn6e7/DGdmqWipaw3u2VMPKsV6NjY0QQAABBBBwqkCgz9+F4R70DeCoUaP00EMPKTk5+QQv6+vfIUOGqFevXqd1vPvuuzVv3jz7K+DKlStr2rRp9lfH33//vS644ILT1qSnp8v6J/dlBahatWpKSUlRTEzg3Ehx8u/9zq9bRsN6tFCpYvzerzDeVBwDAQQQQCCwBWgApaBvAAcNGqQxY8Zo9erVJ6QtLi7Obv769+9/2hRmZGTI+hp43LhxCg8PV2hoqP11cM+ePc+Y2gEDBmjgwIGn/HkgNYDJhzN0/8QF+n3NHvs8+b1fYP8lxNkhgAACCBS+AA1gADeAZ2q2jo/JnDlz7N/ujR07VqtWrTohQfXq1dMdd9xh3+hxutdrr72mkSNHyvq/NWrU0G+//WY3i5MmTVKnTp2C8hNAa30/6/d+W/Yd4fd+hf/3CUdEAAEEEAgSARrAAG4A9+zZI+ufs72smzcmTJjg81fAR44cUWxsrN3sXXXVVccO0bt3b/tmkh9++MGrCAdSgL5euE2Pf7FYaZk5qlY6WsP/1UoJlQPna2mvQNkIAQQQQACBQhAIpPm7EIZ72kME/VfAuTeBzJo1S23atLEHaf37eeedd8abQHIv/HfffSfrDuLc11133aUNGzbYnyp68wqEAFnP8x38/Up98PfzfK31/d66uTm/9/PmArINAggggIArBQJh/vY3fNA3gBag1cQlJSVp+PDhtqe1DIz1te7xy8DUr19fgwcPVrdu3extLrroIvsTxmHDhtnbTp8+3V4D8L///a/XawH6O0B7DqbrvgnzNXP9PntM915URw9fFq8wnufr7/cVx0cAAQQQCGABf8/fgUDjiAZw3759pywEbTV2xy8EHRISotGjRysxMdF237Fjh/2bP+vTPqveagKtxvHBBx+Uta03L38GaNGWZN3z4TwlpaSpmPU83xub6opGPM/Xm+vGNggggAAC7hbw5/wdKPKOaAD9hemvAH08e7Oe+XqZMrJzVLtcMY3o2VJ1y5fwFwPHRQABBBBAIKgE/DV/BxISDaDB1SjsAKVnZWvA5GWaOHuLfdaXJlSwP/mLiYowGAWlCCCAAAIIuEugsOfvQNSlATS4KoUZoKTkI/ZXvou2psj6hvqRy+J1T8c6CuX3fgZXkFIEEEAAATcKFOb8Hai+NIAGV6awAvTn2j324s57D2WoZNEIvXlzc3WMK2dw5pQigAACCCDgXoHCmr8DWZgG0ODqFHSAPB6PRvy2Xi//sFI5HimhUoz9PN9qpYsanDWlCCCAAAIIuFugoOfvYNClATS4SgUZoNS0TD362WL9sGyHfYbdW1TVi90aKSoizOCMKUUAAQQQQACBgpy/g0WXBtDgShVUgNbuStVd4+dp3e5DiggL0bPXNNQ/21b3enkagyFRigACCCCAgOMFCmr+DiY4GkCDq1VQAbJ+7zdlUZIqxkTpnX+1UIvqpQzOklIEEEAAAQQQOF6goObvYFKmATS4WgUVoJQjmXr+m+V6/Ir6KleiiMEZUooAAggggAACJwsU1PwdTNI0gAZXiwAZ4FGKAAIIIICAnwSYvyUaQIPwESADPEoRQAABBBDwkwDzNw2gUfQIkBEfxQgggAACCPhFgPmbBtAoeATIiI9iBBBAAAEE/CLA/E0DaBQ8AmTERzECCCCAAAJ+EWD+pgE0Ch4BMuKjGAEEEEAAAb8IMH/TABoFjwAZ8VGMAAIIIICAXwSYv2kAjYJHgIz4KEYAAQQQQMAvAszfNIBGwSNARnwUI4AAAggg4BcB5m8aQKPgESAjPooRQAABBBDwiwDzNw2gUfAIkBEfxQgggAACCPhFgPmbBtAoeATIiI9iBBBAAAEE/CLA/E0DaBQ8AmTERzECCCCAAAJ+EWD+pgE0Ch4BMuKjGAEEEEAAAb8IMH/TABoFLyUlRSVLltSWLVsUExNjtC+KEUAAAQQQQKBwBKwGsFq1akpOTlZsbGzhHDTAjhLi8Xg8AXZOQXM6W7dutQPECwEEEEAAAQSCT8D6AKdq1arBd+L5cMY0gAaIOTk5SkpKUokSJRQSEmKwp1NLc//rhE8X85X1lJ3hXLC+uXvHGefCESico5Dn4He2PvtKTU1V5cqVFRoaWjgDCrCj0AAG2AU5fsK0Ppa2vmbm6+WCu0j8DqTgbI/fM844F45A4RyFPONcOAIFexQawIL1zfPe+Qsmz3Q+FeLsE1eeN8Y5z3Q+FeLsE1eeN8Y5z3Q+FeLsE5fPG9MA+kxWOAUEH+fCESico5BnnAtHoHCOQp5xLhyBgj0KDWDB+uZ57+np6Ro8eLD69++vIkWK5Hk/FJ5dAOfCSQjOOBeOQOEchTzjXDgCBXsUGsCC9WXvCCCAAAIIIIBAwAnQAAbcJeGEEEAAAQQQQACBghWgASxYX/aOAAIIIIAAAggEnAANYMBdEk4IAQQQQAABBBAoWAEawIL1Ze8IIIAAAggggEDACdAA+umSvPPOO3r11Ve1fft2NWzYUG+88YYuvPDCM57N9OnT9dBDD2nZsmX2yuWPPfaY7r77bj+dfXAd1hdr63o8/PDDmjdvntasWaO+ffva14bXuQV8cf7yyy/17rvvauHChbLuqLTeAwMGDNDll19+7gO5fAtfnP/44w89/vjjWrlypQ4fPqwaNWrorrvu0oMPPuhyxXMP3xfn4/c2Y8YMdezYUY0aNbLzzevsAr44T5s2TRdffPEpO1yxYoXq168PtY8CNIA+guXH5p988ol69uwpK/jnn3++hg8frvfff1/Lly9X9erVTznEhg0b7L9M7rzzTvsvb+svmHvvvVcTJ05U9+7d8+OUHLsPX603btyoIUOGqGXLlvb/tf4ipwE8dzx8de7Xr5/9HzLWX+YlS5bU6NGj9dprr2nWrFlq3rz5uQ/o0i18dV6wYIHd/DVp0kTFihWT1RBaf4dY2e7Tp49LFc89bF+dc/doPbmpRYsWqlu3rnbu3EkDeA5qX51zG8BVq1ad8ISscuXKKSws7NwXli1OEKAB9EMg2rZta/8lYX0Ckvtq0KCBrr32Wnvtv5Nf1n/BT548WdZ/5eS+rE//Fi1apL/++ssPIwieQ/pqffzILrroIjVr1owG0IvLbeKcu3vrU8CbbrpJzzzzjBdHdOcm+eF83XXX2c3g+PHj3Ynoxajz6nzzzTerXr16djPy1Vdf0QCew9pX59wGcP/+/fZ/OPIyE6ABNPPzuTojI0NFixbVZ599pm7duh2rf+CBB+y/LKyvek9+dejQwf5U5M033zz2R5MmTdKNN95of60TERHh83m4oSAv1jSAvifD1Nk6Yk5OjmrWrGn/tOG+++7z/SRcUJEfztYngp07d9YLL7yg3r17u0DN9yHm1dn6FNv6Vsf6j3LLlwbw7PZ5cc5tAK2/K9LS0pSQkKCnnnrqtF8L+37l3VdBA1jI1zwpKUlVqlSxv8Zt3779saMPGjRIY8eOlfXR9smvuLg4JSYm6sknnzz2R3/++af99bG1v0qVKhXyKILjcHmxpgH0/dqaOltHtH4P+9JLL9mfcpcvX973k3BBhYlz1apVtXv3bmVlZdm/tXz66addIJa3IebF2fq98AUXXKDff/9d1t/XljEN4Nn98+JszY+//fab/RMd67fD1qfY7733nqzG0PqghJdvAjSAvnkZb50bequBa9eu3bH9vfjii3aYrd/rnK4B7NWrl/1YuNyX1UBaf+FYNy1UrFjR+LycuIO8WNMA+p4EU2frt6zWp1Fff/21OnXq5PsJuKTCxNn6HfHBgwc1c+ZMPfHEExo2bJh69OjhEjnfhumrc3Z2ts477zzdcccdx27MowE8t7mvzmfa4zXXXKOQkBD7Z1K8fBOgAfTNy3jrvHzszVfAeWPPizUNoO/WJs7Wj8Ct/7ixfhJx1VVX+X5wF1WYOB/PZH09af3H5um+bXAR5xmH6qtzcnKySpUqdcJNCNZPGjwej/2/TZ06VZdccgm0Jwn46nwmQOvDkw8//PCE38iD7Z0ADaB3Tvm6lfXDV+sjbOv3Irkv67cMXbt2PeNNIFOmTLHvEs593XPPPfZvBrkJ5OyXxldrGsC8RT0vztYnf7fffrt9N7t1AxSvcwvkxfnkvT7//PP64IMPZN3xzuv0Ar44W83e8X83W3u0/m7/5Zdf9Pnnn6tWrVr2TTe8ThXwxflMftdff7327dtne/PyTYAG0DevfNk699Z367cL1tfAI0aM0MiRI+01/qx1uqyverdt26Zx48bZx8tdBsZavsFaCsZq+qy7gFkG5tyXw1dra4+5a3dZX0vGx8fr0UcfVWRkpP2DY16nF/DV2crurbfeat/YZN2VmvuKjo5WbGwszGcQ8NX57bfftpeWyl0jzVoGxlqC5/7777dvVOCVP3k+eS98BexdsnzNs7Ukl3UDiLVigPUJovXJn/Xb4S+++OKEv0e8Ozpb0QD6KQPWfyG+8sor9m/4rDX+rHW5cn/Eat3wYf3XufXD1tyXdXewtXhr7kLQ1tIwLATt3cXz1dr6PcnJL6sx5xOTs3v74mwtsXO6O95vu+02jRkzxrsL69KtfHEeOnSovc6o9R+R4eHhqlOnzrH1RENDQ10q6N2wfXGmAfTO9HRb+eJszZnWBybWByTWfyxajaD1gcmVV16Z9xNwcSUNoIsvPkNHAAEEEEAAAXcK0AC687ozagQQQAABBBBwsQANoIsvPkNHAAEEEEAAAXcK0AC687ozagQQQAABBBBwsQANoIsvPkNHAAEEEEAAAXcK0AC687ozagQQQAABBBBwsQANoIsvPkNHAAEEEEAAAXcK0AC687ozagQQQAABBBBwsQANoIsvPkNHAAEEEEAAAXcK0AC687ozagQQOIeA9bSSZs2ayXr8FC8EEEDAaQI0gE67oowHAQTyRYAGMF8Y2QkCCASoAA1ggF4YTgsBBPwnYD2Pe+zYsSecgPU8XetB9LwQQAABJwjQADrhKjIGBBDIV4GUlBR17txZjRo10nPPPWfvu1y5cgoLC8vX47AzBBBAwF8CNID+kue4CCAQ0AJ8BRzQl4eTQwABQwEaQENAyhFAwJkCNIDOvK6MCgEEjgrQAJIEBBBA4DQCNIDEAgEEnCxAA+jkq8vYEEAgzwKXXXaZ4uPjNXTo0Dzvg0IEEEAgUAVoAAP1ynBeCCDgV4E+ffpo4cKF+vTTT1W8eHGVLl1aoaGhfj0nDo4AAgj8Xzt3bARACMNA0J1SMymNkH4TH3HbADNeKVDEXwIG4F+S3iFA4CmBc86stWbvPffe8Q3MU/E6hkBewADMVwAAAQIECBAgUBMwAGuJu5cAAQIECBDICxiA+QoAIECAAAECBGoCBmAtcfcSIECAAAECeQEDMF8BAAQIECBAgEBNwACsJe5eAgQIECBAIC9gAOYrAIAAAQIECBCoCRiAtcTdS4AAAQIECOQFDMB8BQAQIECAAAECNQEDsJa4ewkQIECAAIG8gAGYrwAAAgQIECBAoCZgANYSdy8BAgQIECCQFzAA8xUAQIAAAQIECNQEPlNDDTGtJveIAAAAAElFTkSuQmCC\" width=\"640\">"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plot_wave_at_x(backbone, 0.5, 0.5, show=False)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "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.10.8"
  },
  "vscode": {
   "interpreter": {
    "hash": "95b399d55d7ba88e1a788c9d84fe7761de08157d41dd32f183df354a9771fcfe"
   }
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}