multi_slope.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "e6a9e2c4",
   "metadata": {},
   "outputs": [],
   "source": [
    "import numpy as np\n",
    "import matplotlib\n",
    "import matplotlib.pyplot as plt\n",
    "%matplotlib notebook"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 52,
   "id": "972d686b",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "application/javascript": [
       "/* Put everything inside the global mpl namespace */\n",
       "/* global mpl */\n",
       "window.mpl = {};\n",
       "\n",
       "mpl.get_websocket_type = function () {\n",
       "    if (typeof WebSocket !== 'undefined') {\n",
       "        return WebSocket;\n",
       "    } else if (typeof MozWebSocket !== 'undefined') {\n",
       "        return MozWebSocket;\n",
       "    } else {\n",
       "        alert(\n",
       "            'Your browser does not have WebSocket support. ' +\n",
       "                'Please try Chrome, Safari or Firefox ≥ 6. ' +\n",
       "                'Firefox 4 and 5 are also supported but you ' +\n",
       "                'have to enable WebSockets in about:config.'\n",
       "        );\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure = function (figure_id, websocket, ondownload, parent_element) {\n",
       "    this.id = figure_id;\n",
       "\n",
       "    this.ws = websocket;\n",
       "\n",
       "    this.supports_binary = this.ws.binaryType !== undefined;\n",
       "\n",
       "    if (!this.supports_binary) {\n",
       "        var warnings = document.getElementById('mpl-warnings');\n",
       "        if (warnings) {\n",
       "            warnings.style.display = 'block';\n",
       "            warnings.textContent =\n",
       "                'This browser does not support binary websocket messages. ' +\n",
       "                'Performance may be slow.';\n",
       "        }\n",
       "    }\n",
       "\n",
       "    this.imageObj = new Image();\n",
       "\n",
       "    this.context = undefined;\n",
       "    this.message = undefined;\n",
       "    this.canvas = undefined;\n",
       "    this.rubberband_canvas = undefined;\n",
       "    this.rubberband_context = undefined;\n",
       "    this.format_dropdown = undefined;\n",
       "\n",
       "    this.image_mode = 'full';\n",
       "\n",
       "    this.root = document.createElement('div');\n",
       "    this.root.setAttribute('style', 'display: inline-block');\n",
       "    this._root_extra_style(this.root);\n",
       "\n",
       "    parent_element.appendChild(this.root);\n",
       "\n",
       "    this._init_header(this);\n",
       "    this._init_canvas(this);\n",
       "    this._init_toolbar(this);\n",
       "\n",
       "    var fig = this;\n",
       "\n",
       "    this.waiting = false;\n",
       "\n",
       "    this.ws.onopen = function () {\n",
       "        fig.send_message('supports_binary', { value: fig.supports_binary });\n",
       "        fig.send_message('send_image_mode', {});\n",
       "        if (fig.ratio !== 1) {\n",
       "            fig.send_message('set_device_pixel_ratio', {\n",
       "                device_pixel_ratio: fig.ratio,\n",
       "            });\n",
       "        }\n",
       "        fig.send_message('refresh', {});\n",
       "    };\n",
       "\n",
       "    this.imageObj.onload = function () {\n",
       "        if (fig.image_mode === 'full') {\n",
       "            // Full images could contain transparency (where diff images\n",
       "            // almost always do), so we need to clear the canvas so that\n",
       "            // there is no ghosting.\n",
       "            fig.context.clearRect(0, 0, fig.canvas.width, fig.canvas.height);\n",
       "        }\n",
       "        fig.context.drawImage(fig.imageObj, 0, 0);\n",
       "    };\n",
       "\n",
       "    this.imageObj.onunload = function () {\n",
       "        fig.ws.close();\n",
       "    };\n",
       "\n",
       "    this.ws.onmessage = this._make_on_message_function(this);\n",
       "\n",
       "    this.ondownload = ondownload;\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._init_header = function () {\n",
       "    var titlebar = document.createElement('div');\n",
       "    titlebar.classList =\n",
       "        'ui-dialog-titlebar ui-widget-header ui-corner-all ui-helper-clearfix';\n",
       "    var titletext = document.createElement('div');\n",
       "    titletext.classList = 'ui-dialog-title';\n",
       "    titletext.setAttribute(\n",
       "        'style',\n",
       "        'width: 100%; text-align: center; padding: 3px;'\n",
       "    );\n",
       "    titlebar.appendChild(titletext);\n",
       "    this.root.appendChild(titlebar);\n",
       "    this.header = titletext;\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._canvas_extra_style = function (_canvas_div) {};\n",
       "\n",
       "mpl.figure.prototype._root_extra_style = function (_canvas_div) {};\n",
       "\n",
       "mpl.figure.prototype._init_canvas = function () {\n",
       "    var fig = this;\n",
       "\n",
       "    var canvas_div = (this.canvas_div = document.createElement('div'));\n",
       "    canvas_div.setAttribute(\n",
       "        'style',\n",
       "        'border: 1px solid #ddd;' +\n",
       "            'box-sizing: content-box;' +\n",
       "            'clear: both;' +\n",
       "            'min-height: 1px;' +\n",
       "            'min-width: 1px;' +\n",
       "            'outline: 0;' +\n",
       "            'overflow: hidden;' +\n",
       "            'position: relative;' +\n",
       "            'resize: both;'\n",
       "    );\n",
       "\n",
       "    function on_keyboard_event_closure(name) {\n",
       "        return function (event) {\n",
       "            return fig.key_event(event, name);\n",
       "        };\n",
       "    }\n",
       "\n",
       "    canvas_div.addEventListener(\n",
       "        'keydown',\n",
       "        on_keyboard_event_closure('key_press')\n",
       "    );\n",
       "    canvas_div.addEventListener(\n",
       "        'keyup',\n",
       "        on_keyboard_event_closure('key_release')\n",
       "    );\n",
       "\n",
       "    this._canvas_extra_style(canvas_div);\n",
       "    this.root.appendChild(canvas_div);\n",
       "\n",
       "    var canvas = (this.canvas = document.createElement('canvas'));\n",
       "    canvas.classList.add('mpl-canvas');\n",
       "    canvas.setAttribute('style', 'box-sizing: content-box;');\n",
       "\n",
       "    this.context = canvas.getContext('2d');\n",
       "\n",
       "    var backingStore =\n",
       "        this.context.backingStorePixelRatio ||\n",
       "        this.context.webkitBackingStorePixelRatio ||\n",
       "        this.context.mozBackingStorePixelRatio ||\n",
       "        this.context.msBackingStorePixelRatio ||\n",
       "        this.context.oBackingStorePixelRatio ||\n",
       "        this.context.backingStorePixelRatio ||\n",
       "        1;\n",
       "\n",
       "    this.ratio = (window.devicePixelRatio || 1) / backingStore;\n",
       "\n",
       "    var rubberband_canvas = (this.rubberband_canvas = document.createElement(\n",
       "        'canvas'\n",
       "    ));\n",
       "    rubberband_canvas.setAttribute(\n",
       "        'style',\n",
       "        'box-sizing: content-box; position: absolute; left: 0; top: 0; z-index: 1;'\n",
       "    );\n",
       "\n",
       "    // Apply a ponyfill if ResizeObserver is not implemented by browser.\n",
       "    if (this.ResizeObserver === undefined) {\n",
       "        if (window.ResizeObserver !== undefined) {\n",
       "            this.ResizeObserver = window.ResizeObserver;\n",
       "        } else {\n",
       "            var obs = _JSXTOOLS_RESIZE_OBSERVER({});\n",
       "            this.ResizeObserver = obs.ResizeObserver;\n",
       "        }\n",
       "    }\n",
       "\n",
       "    this.resizeObserverInstance = new this.ResizeObserver(function (entries) {\n",
       "        var nentries = entries.length;\n",
       "        for (var i = 0; i < nentries; i++) {\n",
       "            var entry = entries[i];\n",
       "            var width, height;\n",
       "            if (entry.contentBoxSize) {\n",
       "                if (entry.contentBoxSize instanceof Array) {\n",
       "                    // Chrome 84 implements new version of spec.\n",
       "                    width = entry.contentBoxSize[0].inlineSize;\n",
       "                    height = entry.contentBoxSize[0].blockSize;\n",
       "                } else {\n",
       "                    // Firefox implements old version of spec.\n",
       "                    width = entry.contentBoxSize.inlineSize;\n",
       "                    height = entry.contentBoxSize.blockSize;\n",
       "                }\n",
       "            } else {\n",
       "                // Chrome <84 implements even older version of spec.\n",
       "                width = entry.contentRect.width;\n",
       "                height = entry.contentRect.height;\n",
       "            }\n",
       "\n",
       "            // Keep the size of the canvas and rubber band canvas in sync with\n",
       "            // the canvas container.\n",
       "            if (entry.devicePixelContentBoxSize) {\n",
       "                // Chrome 84 implements new version of spec.\n",
       "                canvas.setAttribute(\n",
       "                    'width',\n",
       "                    entry.devicePixelContentBoxSize[0].inlineSize\n",
       "                );\n",
       "                canvas.setAttribute(\n",
       "                    'height',\n",
       "                    entry.devicePixelContentBoxSize[0].blockSize\n",
       "                );\n",
       "            } else {\n",
       "                canvas.setAttribute('width', width * fig.ratio);\n",
       "                canvas.setAttribute('height', height * fig.ratio);\n",
       "            }\n",
       "            canvas.setAttribute(\n",
       "                'style',\n",
       "                'width: ' + width + 'px; height: ' + height + 'px;'\n",
       "            );\n",
       "\n",
       "            rubberband_canvas.setAttribute('width', width);\n",
       "            rubberband_canvas.setAttribute('height', height);\n",
       "\n",
       "            // And update the size in Python. We ignore the initial 0/0 size\n",
       "            // that occurs as the element is placed into the DOM, which should\n",
       "            // otherwise not happen due to the minimum size styling.\n",
       "            if (fig.ws.readyState == 1 && width != 0 && height != 0) {\n",
       "                fig.request_resize(width, height);\n",
       "            }\n",
       "        }\n",
       "    });\n",
       "    this.resizeObserverInstance.observe(canvas_div);\n",
       "\n",
       "    function on_mouse_event_closure(name) {\n",
       "        return function (event) {\n",
       "            return fig.mouse_event(event, name);\n",
       "        };\n",
       "    }\n",
       "\n",
       "    rubberband_canvas.addEventListener(\n",
       "        'mousedown',\n",
       "        on_mouse_event_closure('button_press')\n",
       "    );\n",
       "    rubberband_canvas.addEventListener(\n",
       "        'mouseup',\n",
       "        on_mouse_event_closure('button_release')\n",
       "    );\n",
       "    rubberband_canvas.addEventListener(\n",
       "        'dblclick',\n",
       "        on_mouse_event_closure('dblclick')\n",
       "    );\n",
       "    // Throttle sequential mouse events to 1 every 20ms.\n",
       "    rubberband_canvas.addEventListener(\n",
       "        'mousemove',\n",
       "        on_mouse_event_closure('motion_notify')\n",
       "    );\n",
       "\n",
       "    rubberband_canvas.addEventListener(\n",
       "        'mouseenter',\n",
       "        on_mouse_event_closure('figure_enter')\n",
       "    );\n",
       "    rubberband_canvas.addEventListener(\n",
       "        'mouseleave',\n",
       "        on_mouse_event_closure('figure_leave')\n",
       "    );\n",
       "\n",
       "    canvas_div.addEventListener('wheel', function (event) {\n",
       "        if (event.deltaY < 0) {\n",
       "            event.step = 1;\n",
       "        } else {\n",
       "            event.step = -1;\n",
       "        }\n",
       "        on_mouse_event_closure('scroll')(event);\n",
       "    });\n",
       "\n",
       "    canvas_div.appendChild(canvas);\n",
       "    canvas_div.appendChild(rubberband_canvas);\n",
       "\n",
       "    this.rubberband_context = rubberband_canvas.getContext('2d');\n",
       "    this.rubberband_context.strokeStyle = '#000000';\n",
       "\n",
       "    this._resize_canvas = function (width, height, forward) {\n",
       "        if (forward) {\n",
       "            canvas_div.style.width = width + 'px';\n",
       "            canvas_div.style.height = height + 'px';\n",
       "        }\n",
       "    };\n",
       "\n",
       "    // Disable right mouse context menu.\n",
       "    this.rubberband_canvas.addEventListener('contextmenu', function (_e) {\n",
       "        event.preventDefault();\n",
       "        return false;\n",
       "    });\n",
       "\n",
       "    function set_focus() {\n",
       "        canvas.focus();\n",
       "        canvas_div.focus();\n",
       "    }\n",
       "\n",
       "    window.setTimeout(set_focus, 100);\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._init_toolbar = function () {\n",
       "    var fig = this;\n",
       "\n",
       "    var toolbar = document.createElement('div');\n",
       "    toolbar.classList = 'mpl-toolbar';\n",
       "    this.root.appendChild(toolbar);\n",
       "\n",
       "    function on_click_closure(name) {\n",
       "        return function (_event) {\n",
       "            return fig.toolbar_button_onclick(name);\n",
       "        };\n",
       "    }\n",
       "\n",
       "    function on_mouseover_closure(tooltip) {\n",
       "        return function (event) {\n",
       "            if (!event.currentTarget.disabled) {\n",
       "                return fig.toolbar_button_onmouseover(tooltip);\n",
       "            }\n",
       "        };\n",
       "    }\n",
       "\n",
       "    fig.buttons = {};\n",
       "    var buttonGroup = document.createElement('div');\n",
       "    buttonGroup.classList = 'mpl-button-group';\n",
       "    for (var toolbar_ind in mpl.toolbar_items) {\n",
       "        var name = mpl.toolbar_items[toolbar_ind][0];\n",
       "        var tooltip = mpl.toolbar_items[toolbar_ind][1];\n",
       "        var image = mpl.toolbar_items[toolbar_ind][2];\n",
       "        var method_name = mpl.toolbar_items[toolbar_ind][3];\n",
       "\n",
       "        if (!name) {\n",
       "            /* Instead of a spacer, we start a new button group. */\n",
       "            if (buttonGroup.hasChildNodes()) {\n",
       "                toolbar.appendChild(buttonGroup);\n",
       "            }\n",
       "            buttonGroup = document.createElement('div');\n",
       "            buttonGroup.classList = 'mpl-button-group';\n",
       "            continue;\n",
       "        }\n",
       "\n",
       "        var button = (fig.buttons[name] = document.createElement('button'));\n",
       "        button.classList = 'mpl-widget';\n",
       "        button.setAttribute('role', 'button');\n",
       "        button.setAttribute('aria-disabled', 'false');\n",
       "        button.addEventListener('click', on_click_closure(method_name));\n",
       "        button.addEventListener('mouseover', on_mouseover_closure(tooltip));\n",
       "\n",
       "        var icon_img = document.createElement('img');\n",
       "        icon_img.src = '_images/' + image + '.png';\n",
       "        icon_img.srcset = '_images/' + image + '_large.png 2x';\n",
       "        icon_img.alt = tooltip;\n",
       "        button.appendChild(icon_img);\n",
       "\n",
       "        buttonGroup.appendChild(button);\n",
       "    }\n",
       "\n",
       "    if (buttonGroup.hasChildNodes()) {\n",
       "        toolbar.appendChild(buttonGroup);\n",
       "    }\n",
       "\n",
       "    var fmt_picker = document.createElement('select');\n",
       "    fmt_picker.classList = 'mpl-widget';\n",
       "    toolbar.appendChild(fmt_picker);\n",
       "    this.format_dropdown = fmt_picker;\n",
       "\n",
       "    for (var ind in mpl.extensions) {\n",
       "        var fmt = mpl.extensions[ind];\n",
       "        var option = document.createElement('option');\n",
       "        option.selected = fmt === mpl.default_extension;\n",
       "        option.innerHTML = fmt;\n",
       "        fmt_picker.appendChild(option);\n",
       "    }\n",
       "\n",
       "    var status_bar = document.createElement('span');\n",
       "    status_bar.classList = 'mpl-message';\n",
       "    toolbar.appendChild(status_bar);\n",
       "    this.message = status_bar;\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.request_resize = function (x_pixels, y_pixels) {\n",
       "    // Request matplotlib to resize the figure. Matplotlib will then trigger a resize in the client,\n",
       "    // which will in turn request a refresh of the image.\n",
       "    this.send_message('resize', { width: x_pixels, height: y_pixels });\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.send_message = function (type, properties) {\n",
       "    properties['type'] = type;\n",
       "    properties['figure_id'] = this.id;\n",
       "    this.ws.send(JSON.stringify(properties));\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.send_draw_message = function () {\n",
       "    if (!this.waiting) {\n",
       "        this.waiting = true;\n",
       "        this.ws.send(JSON.stringify({ type: 'draw', figure_id: this.id }));\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_save = function (fig, _msg) {\n",
       "    var format_dropdown = fig.format_dropdown;\n",
       "    var format = format_dropdown.options[format_dropdown.selectedIndex].value;\n",
       "    fig.ondownload(fig, format);\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_resize = function (fig, msg) {\n",
       "    var size = msg['size'];\n",
       "    if (size[0] !== fig.canvas.width || size[1] !== fig.canvas.height) {\n",
       "        fig._resize_canvas(size[0], size[1], msg['forward']);\n",
       "        fig.send_message('refresh', {});\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_rubberband = function (fig, msg) {\n",
       "    var x0 = msg['x0'] / fig.ratio;\n",
       "    var y0 = (fig.canvas.height - msg['y0']) / fig.ratio;\n",
       "    var x1 = msg['x1'] / fig.ratio;\n",
       "    var y1 = (fig.canvas.height - msg['y1']) / fig.ratio;\n",
       "    x0 = Math.floor(x0) + 0.5;\n",
       "    y0 = Math.floor(y0) + 0.5;\n",
       "    x1 = Math.floor(x1) + 0.5;\n",
       "    y1 = Math.floor(y1) + 0.5;\n",
       "    var min_x = Math.min(x0, x1);\n",
       "    var min_y = Math.min(y0, y1);\n",
       "    var width = Math.abs(x1 - x0);\n",
       "    var height = Math.abs(y1 - y0);\n",
       "\n",
       "    fig.rubberband_context.clearRect(\n",
       "        0,\n",
       "        0,\n",
       "        fig.canvas.width / fig.ratio,\n",
       "        fig.canvas.height / fig.ratio\n",
       "    );\n",
       "\n",
       "    fig.rubberband_context.strokeRect(min_x, min_y, width, height);\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_figure_label = function (fig, msg) {\n",
       "    // Updates the figure title.\n",
       "    fig.header.textContent = msg['label'];\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_cursor = function (fig, msg) {\n",
       "    fig.rubberband_canvas.style.cursor = msg['cursor'];\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_message = function (fig, msg) {\n",
       "    fig.message.textContent = msg['message'];\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_draw = function (fig, _msg) {\n",
       "    // Request the server to send over a new figure.\n",
       "    fig.send_draw_message();\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_image_mode = function (fig, msg) {\n",
       "    fig.image_mode = msg['mode'];\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_history_buttons = function (fig, msg) {\n",
       "    for (var key in msg) {\n",
       "        if (!(key in fig.buttons)) {\n",
       "            continue;\n",
       "        }\n",
       "        fig.buttons[key].disabled = !msg[key];\n",
       "        fig.buttons[key].setAttribute('aria-disabled', !msg[key]);\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_navigate_mode = function (fig, msg) {\n",
       "    if (msg['mode'] === 'PAN') {\n",
       "        fig.buttons['Pan'].classList.add('active');\n",
       "        fig.buttons['Zoom'].classList.remove('active');\n",
       "    } else if (msg['mode'] === 'ZOOM') {\n",
       "        fig.buttons['Pan'].classList.remove('active');\n",
       "        fig.buttons['Zoom'].classList.add('active');\n",
       "    } else {\n",
       "        fig.buttons['Pan'].classList.remove('active');\n",
       "        fig.buttons['Zoom'].classList.remove('active');\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.updated_canvas_event = function () {\n",
       "    // Called whenever the canvas gets updated.\n",
       "    this.send_message('ack', {});\n",
       "};\n",
       "\n",
       "// A function to construct a web socket function for onmessage handling.\n",
       "// Called in the figure constructor.\n",
       "mpl.figure.prototype._make_on_message_function = function (fig) {\n",
       "    return function socket_on_message(evt) {\n",
       "        if (evt.data instanceof Blob) {\n",
       "            var img = evt.data;\n",
       "            if (img.type !== 'image/png') {\n",
       "                /* FIXME: We get \"Resource interpreted as Image but\n",
       "                 * transferred with MIME type text/plain:\" errors on\n",
       "                 * Chrome.  But how to set the MIME type?  It doesn't seem\n",
       "                 * to be part of the websocket stream */\n",
       "                img.type = 'image/png';\n",
       "            }\n",
       "\n",
       "            /* Free the memory for the previous frames */\n",
       "            if (fig.imageObj.src) {\n",
       "                (window.URL || window.webkitURL).revokeObjectURL(\n",
       "                    fig.imageObj.src\n",
       "                );\n",
       "            }\n",
       "\n",
       "            fig.imageObj.src = (window.URL || window.webkitURL).createObjectURL(\n",
       "                img\n",
       "            );\n",
       "            fig.updated_canvas_event();\n",
       "            fig.waiting = false;\n",
       "            return;\n",
       "        } else if (\n",
       "            typeof evt.data === 'string' &&\n",
       "            evt.data.slice(0, 21) === 'data:image/png;base64'\n",
       "        ) {\n",
       "            fig.imageObj.src = evt.data;\n",
       "            fig.updated_canvas_event();\n",
       "            fig.waiting = false;\n",
       "            return;\n",
       "        }\n",
       "\n",
       "        var msg = JSON.parse(evt.data);\n",
       "        var msg_type = msg['type'];\n",
       "\n",
       "        // Call the  \"handle_{type}\" callback, which takes\n",
       "        // the figure and JSON message as its only arguments.\n",
       "        try {\n",
       "            var callback = fig['handle_' + msg_type];\n",
       "        } catch (e) {\n",
       "            console.log(\n",
       "                \"No handler for the '\" + msg_type + \"' message type: \",\n",
       "                msg\n",
       "            );\n",
       "            return;\n",
       "        }\n",
       "\n",
       "        if (callback) {\n",
       "            try {\n",
       "                // console.log(\"Handling '\" + msg_type + \"' message: \", msg);\n",
       "                callback(fig, msg);\n",
       "            } catch (e) {\n",
       "                console.log(\n",
       "                    \"Exception inside the 'handler_\" + msg_type + \"' callback:\",\n",
       "                    e,\n",
       "                    e.stack,\n",
       "                    msg\n",
       "                );\n",
       "            }\n",
       "        }\n",
       "    };\n",
       "};\n",
       "\n",
       "// from https://stackoverflow.com/questions/1114465/getting-mouse-location-in-canvas\n",
       "mpl.findpos = function (e) {\n",
       "    //this section is from http://www.quirksmode.org/js/events_properties.html\n",
       "    var targ;\n",
       "    if (!e) {\n",
       "        e = window.event;\n",
       "    }\n",
       "    if (e.target) {\n",
       "        targ = e.target;\n",
       "    } else if (e.srcElement) {\n",
       "        targ = e.srcElement;\n",
       "    }\n",
       "    if (targ.nodeType === 3) {\n",
       "        // defeat Safari bug\n",
       "        targ = targ.parentNode;\n",
       "    }\n",
       "\n",
       "    // pageX,Y are the mouse positions relative to the document\n",
       "    var boundingRect = targ.getBoundingClientRect();\n",
       "    var x = e.pageX - (boundingRect.left + document.body.scrollLeft);\n",
       "    var y = e.pageY - (boundingRect.top + document.body.scrollTop);\n",
       "\n",
       "    return { x: x, y: y };\n",
       "};\n",
       "\n",
       "/*\n",
       " * return a copy of an object with only non-object keys\n",
       " * we need this to avoid circular references\n",
       " * https://stackoverflow.com/a/24161582/3208463\n",
       " */\n",
       "function simpleKeys(original) {\n",
       "    return Object.keys(original).reduce(function (obj, key) {\n",
       "        if (typeof original[key] !== 'object') {\n",
       "            obj[key] = original[key];\n",
       "        }\n",
       "        return obj;\n",
       "    }, {});\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.mouse_event = function (event, name) {\n",
       "    var canvas_pos = mpl.findpos(event);\n",
       "\n",
       "    if (name === 'button_press') {\n",
       "        this.canvas.focus();\n",
       "        this.canvas_div.focus();\n",
       "    }\n",
       "\n",
       "    var x = canvas_pos.x * this.ratio;\n",
       "    var y = canvas_pos.y * this.ratio;\n",
       "\n",
       "    this.send_message(name, {\n",
       "        x: x,\n",
       "        y: y,\n",
       "        button: event.button,\n",
       "        step: event.step,\n",
       "        guiEvent: simpleKeys(event),\n",
       "    });\n",
       "\n",
       "    /* This prevents the web browser from automatically changing to\n",
       "     * the text insertion cursor when the button is pressed.  We want\n",
       "     * to control all of the cursor setting manually through the\n",
       "     * 'cursor' event from matplotlib */\n",
       "    event.preventDefault();\n",
       "    return false;\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._key_event_extra = function (_event, _name) {\n",
       "    // Handle any extra behaviour associated with a key event\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.key_event = function (event, name) {\n",
       "    // Prevent repeat events\n",
       "    if (name === 'key_press') {\n",
       "        if (event.key === this._key) {\n",
       "            return;\n",
       "        } else {\n",
       "            this._key = event.key;\n",
       "        }\n",
       "    }\n",
       "    if (name === 'key_release') {\n",
       "        this._key = null;\n",
       "    }\n",
       "\n",
       "    var value = '';\n",
       "    if (event.ctrlKey && event.key !== 'Control') {\n",
       "        value += 'ctrl+';\n",
       "    }\n",
       "    else if (event.altKey && event.key !== 'Alt') {\n",
       "        value += 'alt+';\n",
       "    }\n",
       "    else if (event.shiftKey && event.key !== 'Shift') {\n",
       "        value += 'shift+';\n",
       "    }\n",
       "\n",
       "    value += 'k' + event.key;\n",
       "\n",
       "    this._key_event_extra(event, name);\n",
       "\n",
       "    this.send_message(name, { key: value, guiEvent: simpleKeys(event) });\n",
       "    return false;\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.toolbar_button_onclick = function (name) {\n",
       "    if (name === 'download') {\n",
       "        this.handle_save(this, null);\n",
       "    } else {\n",
       "        this.send_message('toolbar_button', { name: name });\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.toolbar_button_onmouseover = function (tooltip) {\n",
       "    this.message.textContent = tooltip;\n",
       "};\n",
       "\n",
       "///////////////// REMAINING CONTENT GENERATED BY embed_js.py /////////////////\n",
       "// prettier-ignore\n",
       "var _JSXTOOLS_RESIZE_OBSERVER=function(A){var t,i=new WeakMap,n=new WeakMap,a=new WeakMap,r=new WeakMap,o=new Set;function s(e){if(!(this instanceof s))throw new TypeError(\"Constructor requires 'new' operator\");i.set(this,e)}function h(){throw new TypeError(\"Function is not a constructor\")}function c(e,t,i,n){e=0 in arguments?Number(arguments[0]):0,t=1 in arguments?Number(arguments[1]):0,i=2 in arguments?Number(arguments[2]):0,n=3 in arguments?Number(arguments[3]):0,this.right=(this.x=this.left=e)+(this.width=i),this.bottom=(this.y=this.top=t)+(this.height=n),Object.freeze(this)}function d(){t=requestAnimationFrame(d);var s=new WeakMap,p=new Set;o.forEach((function(t){r.get(t).forEach((function(i){var r=t instanceof window.SVGElement,o=a.get(t),d=r?0:parseFloat(o.paddingTop),f=r?0:parseFloat(o.paddingRight),l=r?0:parseFloat(o.paddingBottom),u=r?0:parseFloat(o.paddingLeft),g=r?0:parseFloat(o.borderTopWidth),m=r?0:parseFloat(o.borderRightWidth),w=r?0:parseFloat(o.borderBottomWidth),b=u+f,F=d+l,v=(r?0:parseFloat(o.borderLeftWidth))+m,W=g+w,y=r?0:t.offsetHeight-W-t.clientHeight,E=r?0:t.offsetWidth-v-t.clientWidth,R=b+v,z=F+W,M=r?t.width:parseFloat(o.width)-R-E,O=r?t.height:parseFloat(o.height)-z-y;if(n.has(t)){var k=n.get(t);if(k[0]===M&&k[1]===O)return}n.set(t,[M,O]);var S=Object.create(h.prototype);S.target=t,S.contentRect=new c(u,d,M,O),s.has(i)||(s.set(i,[]),p.add(i)),s.get(i).push(S)}))})),p.forEach((function(e){i.get(e).call(e,s.get(e),e)}))}return s.prototype.observe=function(i){if(i instanceof window.Element){r.has(i)||(r.set(i,new Set),o.add(i),a.set(i,window.getComputedStyle(i)));var n=r.get(i);n.has(this)||n.add(this),cancelAnimationFrame(t),t=requestAnimationFrame(d)}},s.prototype.unobserve=function(i){if(i instanceof window.Element&&r.has(i)){var n=r.get(i);n.has(this)&&(n.delete(this),n.size||(r.delete(i),o.delete(i))),n.size||r.delete(i),o.size||cancelAnimationFrame(t)}},A.DOMRectReadOnly=c,A.ResizeObserver=s,A.ResizeObserverEntry=h,A}; // eslint-disable-line\n",
       "mpl.toolbar_items = [[\"Home\", \"Reset original view\", \"fa fa-home icon-home\", \"home\"], [\"Back\", \"Back to previous view\", \"fa fa-arrow-left icon-arrow-left\", \"back\"], [\"Forward\", \"Forward to next view\", \"fa fa-arrow-right icon-arrow-right\", \"forward\"], [\"\", \"\", \"\", \"\"], [\"Pan\", \"Left button pans, Right button zooms\\nx/y fixes axis, CTRL fixes aspect\", \"fa fa-arrows icon-move\", \"pan\"], [\"Zoom\", \"Zoom to rectangle\\nx/y fixes axis\", \"fa fa-square-o icon-check-empty\", \"zoom\"], [\"\", \"\", \"\", \"\"], [\"Download\", \"Download plot\", \"fa fa-floppy-o icon-save\", \"download\"]];\n",
       "\n",
       "mpl.extensions = [\"eps\", \"jpeg\", \"pgf\", \"pdf\", \"png\", \"ps\", \"raw\", \"svg\", \"tif\"];\n",
       "\n",
       "mpl.default_extension = \"png\";/* global mpl */\n",
       "\n",
       "var comm_websocket_adapter = function (comm) {\n",
       "    // Create a \"websocket\"-like object which calls the given IPython comm\n",
       "    // object with the appropriate methods. Currently this is a non binary\n",
       "    // socket, so there is still some room for performance tuning.\n",
       "    var ws = {};\n",
       "\n",
       "    ws.binaryType = comm.kernel.ws.binaryType;\n",
       "    ws.readyState = comm.kernel.ws.readyState;\n",
       "    function updateReadyState(_event) {\n",
       "        if (comm.kernel.ws) {\n",
       "            ws.readyState = comm.kernel.ws.readyState;\n",
       "        } else {\n",
       "            ws.readyState = 3; // Closed state.\n",
       "        }\n",
       "    }\n",
       "    comm.kernel.ws.addEventListener('open', updateReadyState);\n",
       "    comm.kernel.ws.addEventListener('close', updateReadyState);\n",
       "    comm.kernel.ws.addEventListener('error', updateReadyState);\n",
       "\n",
       "    ws.close = function () {\n",
       "        comm.close();\n",
       "    };\n",
       "    ws.send = function (m) {\n",
       "        //console.log('sending', m);\n",
       "        comm.send(m);\n",
       "    };\n",
       "    // Register the callback with on_msg.\n",
       "    comm.on_msg(function (msg) {\n",
       "        //console.log('receiving', msg['content']['data'], msg);\n",
       "        var data = msg['content']['data'];\n",
       "        if (data['blob'] !== undefined) {\n",
       "            data = {\n",
       "                data: new Blob(msg['buffers'], { type: data['blob'] }),\n",
       "            };\n",
       "        }\n",
       "        // Pass the mpl event to the overridden (by mpl) onmessage function.\n",
       "        ws.onmessage(data);\n",
       "    });\n",
       "    return ws;\n",
       "};\n",
       "\n",
       "mpl.mpl_figure_comm = function (comm, msg) {\n",
       "    // This is the function which gets called when the mpl process\n",
       "    // starts-up an IPython Comm through the \"matplotlib\" channel.\n",
       "\n",
       "    var id = msg.content.data.id;\n",
       "    // Get hold of the div created by the display call when the Comm\n",
       "    // socket was opened in Python.\n",
       "    var element = document.getElementById(id);\n",
       "    var ws_proxy = comm_websocket_adapter(comm);\n",
       "\n",
       "    function ondownload(figure, _format) {\n",
       "        window.open(figure.canvas.toDataURL());\n",
       "    }\n",
       "\n",
       "    var fig = new mpl.figure(id, ws_proxy, ondownload, element);\n",
       "\n",
       "    // Call onopen now - mpl needs it, as it is assuming we've passed it a real\n",
       "    // web socket which is closed, not our websocket->open comm proxy.\n",
       "    ws_proxy.onopen();\n",
       "\n",
       "    fig.parent_element = element;\n",
       "    fig.cell_info = mpl.find_output_cell(\"<div id='\" + id + \"'></div>\");\n",
       "    if (!fig.cell_info) {\n",
       "        console.error('Failed to find cell for figure', id, fig);\n",
       "        return;\n",
       "    }\n",
       "    fig.cell_info[0].output_area.element.on(\n",
       "        'cleared',\n",
       "        { fig: fig },\n",
       "        fig._remove_fig_handler\n",
       "    );\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_close = function (fig, msg) {\n",
       "    var width = fig.canvas.width / fig.ratio;\n",
       "    fig.cell_info[0].output_area.element.off(\n",
       "        'cleared',\n",
       "        fig._remove_fig_handler\n",
       "    );\n",
       "    fig.resizeObserverInstance.unobserve(fig.canvas_div);\n",
       "\n",
       "    // Update the output cell to use the data from the current canvas.\n",
       "    fig.push_to_output();\n",
       "    var dataURL = fig.canvas.toDataURL();\n",
       "    // Re-enable the keyboard manager in IPython - without this line, in FF,\n",
       "    // the notebook keyboard shortcuts fail.\n",
       "    IPython.keyboard_manager.enable();\n",
       "    fig.parent_element.innerHTML =\n",
       "        '<img src=\"' + dataURL + '\" width=\"' + width + '\">';\n",
       "    fig.close_ws(fig, msg);\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.close_ws = function (fig, msg) {\n",
       "    fig.send_message('closing', msg);\n",
       "    // fig.ws.close()\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.push_to_output = function (_remove_interactive) {\n",
       "    // Turn the data on the canvas into data in the output cell.\n",
       "    var width = this.canvas.width / this.ratio;\n",
       "    var dataURL = this.canvas.toDataURL();\n",
       "    this.cell_info[1]['text/html'] =\n",
       "        '<img src=\"' + dataURL + '\" width=\"' + width + '\">';\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.updated_canvas_event = function () {\n",
       "    // Tell IPython that the notebook contents must change.\n",
       "    IPython.notebook.set_dirty(true);\n",
       "    this.send_message('ack', {});\n",
       "    var fig = this;\n",
       "    // Wait a second, then push the new image to the DOM so\n",
       "    // that it is saved nicely (might be nice to debounce this).\n",
       "    setTimeout(function () {\n",
       "        fig.push_to_output();\n",
       "    }, 1000);\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._init_toolbar = function () {\n",
       "    var fig = this;\n",
       "\n",
       "    var toolbar = document.createElement('div');\n",
       "    toolbar.classList = 'btn-toolbar';\n",
       "    this.root.appendChild(toolbar);\n",
       "\n",
       "    function on_click_closure(name) {\n",
       "        return function (_event) {\n",
       "            return fig.toolbar_button_onclick(name);\n",
       "        };\n",
       "    }\n",
       "\n",
       "    function on_mouseover_closure(tooltip) {\n",
       "        return function (event) {\n",
       "            if (!event.currentTarget.disabled) {\n",
       "                return fig.toolbar_button_onmouseover(tooltip);\n",
       "            }\n",
       "        };\n",
       "    }\n",
       "\n",
       "    fig.buttons = {};\n",
       "    var buttonGroup = document.createElement('div');\n",
       "    buttonGroup.classList = 'btn-group';\n",
       "    var button;\n",
       "    for (var toolbar_ind in mpl.toolbar_items) {\n",
       "        var name = mpl.toolbar_items[toolbar_ind][0];\n",
       "        var tooltip = mpl.toolbar_items[toolbar_ind][1];\n",
       "        var image = mpl.toolbar_items[toolbar_ind][2];\n",
       "        var method_name = mpl.toolbar_items[toolbar_ind][3];\n",
       "\n",
       "        if (!name) {\n",
       "            /* Instead of a spacer, we start a new button group. */\n",
       "            if (buttonGroup.hasChildNodes()) {\n",
       "                toolbar.appendChild(buttonGroup);\n",
       "            }\n",
       "            buttonGroup = document.createElement('div');\n",
       "            buttonGroup.classList = 'btn-group';\n",
       "            continue;\n",
       "        }\n",
       "\n",
       "        button = fig.buttons[name] = document.createElement('button');\n",
       "        button.classList = 'btn btn-default';\n",
       "        button.href = '#';\n",
       "        button.title = name;\n",
       "        button.innerHTML = '<i class=\"fa ' + image + ' fa-lg\"></i>';\n",
       "        button.addEventListener('click', on_click_closure(method_name));\n",
       "        button.addEventListener('mouseover', on_mouseover_closure(tooltip));\n",
       "        buttonGroup.appendChild(button);\n",
       "    }\n",
       "\n",
       "    if (buttonGroup.hasChildNodes()) {\n",
       "        toolbar.appendChild(buttonGroup);\n",
       "    }\n",
       "\n",
       "    // Add the status bar.\n",
       "    var status_bar = document.createElement('span');\n",
       "    status_bar.classList = 'mpl-message pull-right';\n",
       "    toolbar.appendChild(status_bar);\n",
       "    this.message = status_bar;\n",
       "\n",
       "    // Add the close button to the window.\n",
       "    var buttongrp = document.createElement('div');\n",
       "    buttongrp.classList = 'btn-group inline pull-right';\n",
       "    button = document.createElement('button');\n",
       "    button.classList = 'btn btn-mini btn-primary';\n",
       "    button.href = '#';\n",
       "    button.title = 'Stop Interaction';\n",
       "    button.innerHTML = '<i class=\"fa fa-power-off icon-remove icon-large\"></i>';\n",
       "    button.addEventListener('click', function (_evt) {\n",
       "        fig.handle_close(fig, {});\n",
       "    });\n",
       "    button.addEventListener(\n",
       "        'mouseover',\n",
       "        on_mouseover_closure('Stop Interaction')\n",
       "    );\n",
       "    buttongrp.appendChild(button);\n",
       "    var titlebar = this.root.querySelector('.ui-dialog-titlebar');\n",
       "    titlebar.insertBefore(buttongrp, titlebar.firstChild);\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._remove_fig_handler = function (event) {\n",
       "    var fig = event.data.fig;\n",
       "    if (event.target !== this) {\n",
       "        // Ignore bubbled events from children.\n",
       "        return;\n",
       "    }\n",
       "    fig.close_ws(fig, {});\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._root_extra_style = function (el) {\n",
       "    el.style.boxSizing = 'content-box'; // override notebook setting of border-box.\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._canvas_extra_style = function (el) {\n",
       "    // this is important to make the div 'focusable\n",
       "    el.setAttribute('tabindex', 0);\n",
       "    // reach out to IPython and tell the keyboard manager to turn it's self\n",
       "    // off when our div gets focus\n",
       "\n",
       "    // location in version 3\n",
       "    if (IPython.notebook.keyboard_manager) {\n",
       "        IPython.notebook.keyboard_manager.register_events(el);\n",
       "    } else {\n",
       "        // location in version 2\n",
       "        IPython.keyboard_manager.register_events(el);\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._key_event_extra = function (event, _name) {\n",
       "    // Check for shift+enter\n",
       "    if (event.shiftKey && event.which === 13) {\n",
       "        this.canvas_div.blur();\n",
       "        // select the cell after this one\n",
       "        var index = IPython.notebook.find_cell_index(this.cell_info[0]);\n",
       "        IPython.notebook.select(index + 1);\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_save = function (fig, _msg) {\n",
       "    fig.ondownload(fig, null);\n",
       "};\n",
       "\n",
       "mpl.find_output_cell = function (html_output) {\n",
       "    // Return the cell and output element which can be found *uniquely* in the notebook.\n",
       "    // Note - this is a bit hacky, but it is done because the \"notebook_saving.Notebook\"\n",
       "    // IPython event is triggered only after the cells have been serialised, which for\n",
       "    // our purposes (turning an active figure into a static one), is too late.\n",
       "    var cells = IPython.notebook.get_cells();\n",
       "    var ncells = cells.length;\n",
       "    for (var i = 0; i < ncells; i++) {\n",
       "        var cell = cells[i];\n",
       "        if (cell.cell_type === 'code') {\n",
       "            for (var j = 0; j < cell.output_area.outputs.length; j++) {\n",
       "                var data = cell.output_area.outputs[j];\n",
       "                if (data.data) {\n",
       "                    // IPython >= 3 moved mimebundle to data attribute of output\n",
       "                    data = data.data;\n",
       "                }\n",
       "                if (data['text/html'] === html_output) {\n",
       "                    return [cell, data, j];\n",
       "                }\n",
       "            }\n",
       "        }\n",
       "    }\n",
       "};\n",
       "\n",
       "// Register the function which deals with the matplotlib target/channel.\n",
       "// The kernel may be null if the page has been refreshed.\n",
       "if (IPython.notebook.kernel !== null) {\n",
       "    IPython.notebook.kernel.comm_manager.register_target(\n",
       "        'matplotlib',\n",
       "        mpl.mpl_figure_comm\n",
       "    );\n",
       "}\n"
      ],
      "text/plain": [
       "<IPython.core.display.Javascript object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/html": [
       "<img src=\"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAoAAAAHgCAYAAAA10dzkAAAAAXNSR0IArs4c6QAAAARzQklUCAgICHwIZIgAACAASURBVHhe7J0HfFVF2sYfUklIhZAAIQlIDWChKIbewUqxFxSR4qerWLbIt7t+WIBdWQVXpVsAyyKKXVlEDSARVDSKhA6hJCGVdNL53jkhuQkEcpO5Lfc+8/vlx7Kcd845//vczOPMvO80OyMNbCRAAiRAAiRAAiRAAi5DoBkNoMt81nxREiABEiABEiABEjAI0ABSCCRAAiRAAiRAAiTgYgRoAF3sA+frkgAJkAAJkAAJkAANIDVAAiRAAiRAAiRAAi5GgAbQxT5wvi4JkAAJkAAJkAAJ0ABSAyRAAiRAAiRAAiTgYgRoAF3sA+frkgAJkAAJkAAJkAANIDVAAiRAAiRAAiRAAi5GgAbQxT5wvi4JkAAJkAAJkAAJ0ABSAyRAAiRAAiRAAiTgYgRoAF3sA+frkgAJkAAJkAAJkAANIDVAAiRAAiRAAiRAAi5GgAbQxT5wvi4JkAAJkAAJkAAJ0ABSAyRAAiRAAiRAAiTgYgRoAF3sA+frkgAJkAAJkAAJkAANIDVAAiRAAiRAAiRAAi5GgAbQxT5wvi4JkAAJkAAJkAAJ0ABSAyRAAiRAAiRAAiTgYgRoAF3sA+frkgAJkAAJkAAJkAANIDVAAiRAAiRAAiRAAi5GgAbQxT5wvi4JkAAJkAAJkAAJ0ABSAyRAAiRAAiRAAiTgYgRoAF3sA+frkgAJkAAJkAAJkAANIDVAAiRAAiRAAiRAAi5GgAbQxT5wvi4JkAAJkAAJkAAJ0ABSAyRAAiRAAiRAAiTgYgRoAF3sA+frkgAJkAAJkAAJkAANIDVAAiRAAiRAAiRAAi5GgAbQxT5wvi4JkAAJkAAJkAAJ0ABSAyRAAiRAAiRAAiTgYgRoAF3sA+frkgAJkAAJkAAJkAANIDVAAiRAAiRAAiRAAi5GgAbQxT5wvi4JkAAJkAAJkAAJ0ABSAyRAAiRAAiRAAiTgYgRoAF3sA+frkgAJkAAJkAAJkAANIDVAAiRAAiRAAiRAAi5GgAbQxT5wvi4JkAAJkAAJkAAJ0ABSAyRAAiRAAiRAAiTgYgRoAF3sA+frkgAJkAAJkAAJkAANIDVAAiRAAiRAAiRAAi5GgAbQxT5wvi4JkAAJkAAJkAAJ0ABSAyRAAiRAAiRAAiTgYgRoAF3sA+frkgAJkAAJkAAJkAANIDVAAiRAAiRAAiRAAi5GgAbQxT5wvi4JkAAJkAAJkAAJ0ABSAyRAAiRAAiRAAiTgYgRoAF3sA+frkgAJkAAJkAAJkAANIDVAAiRAAiRAAiRAAi5GgAbQxT5wvi4JkAAJkAAJkAAJ0ABSAyRAAiRAAiRAAiTgYgRoAF3sA+frkgAJkAAJkAAJkAANIDVAAiRAAiRAAiRAAi5GgAbQxT5wvi4JkAAJkAAJkAAJ0ABSAyRAAiRAAiRAAiTgYgRoAF3sA+frkgAJkAAJkAAJkAANIDVAAiRAAiRAAiRAAi5GgAZQ4wOvqKhAcnIy/P390axZM42eGEoCJEACJEACJGArAmfOnEFeXh7atWsHNzc3W93Woe5DA6jxcZw4cQIREREaPTCUBEiABEiABEjAXgSOHz+O9u3b2+v2dr0vDaAG/pycHAQFBUEJKCAgQKMnhpIACZAACZAACdiKQG5urjGBk52djcDAQFvd1qHuQwOo8XEoASnhKCNIA6gBkqEkQAIkQAIkYEMCHL8BGkANwVFAGvAYSgIkQAIkQAJ2IsDxmwZQS3oUkBY+BpMACZAACZCAXQhw/KYB1BIeBaSFj8EkQAIkQAIkYBcCHL9pALWERwFp4WMwCZAACZAACdiFAMdvGkAt4VFAWvgYTAIkQAIkQAJ2IcDxmwZQS3gUkBY+BpMACZAACZCAXQhw/KYB1BIeBaSFj8EkQAIkQAIkYBcCHL9pALWERwFp4WMwCZAACZAACdiFAMdvGkAt4VFAWvgYTAIkQAIkQAJ2IcDxmwZQS3gUkBY+BpMACZAACZCAXQhw/KYB1BIeBaSFj8EkQAIkQAIkYBcCHL9pALWERwFp4WMwCZAACZAACdiFAMdvGkAt4VFAWvgYTAIkQAIkQAJ1EsjIL8amhFRs2H0SM4d0QkynVhYlxfGbBlBLUBSQFj4GkwAJkAAJkEA1geNZhfivGL6Nu1Px09EsVJyp/Kd7YqLwzPheFiXF8ZsGUEtQFJAWPgaTAAmQAAm4MIEzZ85gf2o+Nvx+0jB+CSm5tWhc1j4QY3u2wbhebdCptZ9FSXH8pgHUEhQFpIWPwSRAAiRAAi5GoEKm9X45ni2zfJWmLzGzsJqAWzPgqo4tDdM3Rn7Cg3ysRofjNw2glrgoIC18DCYBEiABEnABAqXlFdh+OLN6eTctr7j6rb083DCkS4hh+EZFh6FlCy+bEOH4TQOoJTQKSAsfg0mABEiABJyUQGFJGbbsTxfTl4qv96Qit6is+k39vT0wvHuoMdM3tFtr+Mnfbd04ftMAammOAtLCx2ASIAESIAEnIpBdWCJmL82Y6dtyIB1FpRXVbxfi54XRPdqI6QszMnq9Pdzt+uYcv2kAtQRIAWnhYzAJkAAJkEATJ3AypwgbEyr3820/nIXyqtRdea+Ilj4Yq0yfJHH0iQyGu9rk5yCN4zcNoJYUKSAtfAwmARIgARJoggQOp+cbS7vK9MVLQkfN1r2Nv7Gfb5z8RLf1R7NmjmP6aj4nx28aQK2vHgWkhY/BJEACJEACTYCAKteyOzm3ulzLgbT86qdW/k7N7qmlXbWnL6pViybwRgDHbxpALaFSQFr4GEwCJEACJOCgBNRS7o+JWdWZu0nZp6uf1EOWctU+PqNcS48whAY0d9C3uPBjcfymAdQSLQWkhY/BJEACJEACDkSgqLQccYcyjJm+TZLMkVVQUv10Pp7uGCYZu8r0qQzeQB9PB3ryhj8Kx28awIarpkYEBaSFj8EkQAIkQAJ2JpBXVIpv96lyLScRuzcNBSXl1U8U5OuJkd3V0m4YhnRtjeZiAp2lcfymAdTSMgWkhY/BJEACJEACdiCQkV+MTQmp2CCmL+5gJkqkUHNVayPLuVX7+dSpHB7ubnZ4QuvfkuM3DaCWyiggLXwMJgESIAESsBGB41mF1fv5fjqahRrVWnBJ6xbG0q76uSw8EG4OVK7FWng4ftMAammLAtLCx2ASIAESIAErEVCZu/tTVbmWk8aevoSU3Fp3ulSM3jipz6dm+zqH+lvpKRy3W47fNIBa6qSAtPAxmARIgARIwIIEKmRa7xepy7dRTJ8yfomZhdW9q0k9taRrZO7KT3iQjwXv3PS64vhNA6ilWgpICx+DSYAESIAENAmUyv697Yczq5d30/KKq3v08nDD4M4hhukbGR2KVn7emndznnCO305mABcvXowFCxYgJSUFPXv2xKJFizB48OA6FTtlyhSsWrXqvH/r0aMHdu/ebZbKKSCzMPEiEiABEiABCxIoLCnDlv0qczdVzt5NRW5RWXXv/t4eRpkWZfqGStkWP/k72/kEOH47kQFcu3YtJk+eDGUCBw4ciGXLlmHlypVISEhAZGTkeZ9+Tk4OTp82FbYsKyvD5Zdfjocffhhz5swx6/tCAZmFiReRAAmQAAloEsguLBGzl2bM9G05kI6iUlPmboifF0ZLQWa1tDtACjR7ezhPuRZNbBcM5/jtRAawf//+6NOnD5YsWVL9gUdHR2PChAmYP39+vRr66KOPMGnSJBw5cgRRUVH1Xq8uoIDMwsSLSIAESIAEGkHgZE4RNiZU7ufbfjgL6nSOqhbR0gdje0gShyRyqKPY3F0gc7cRCGkALwKtmWQKmRRlSbo27KukpAS+vr5Yt24dJk6cWH3nWbNmIT4+Hps3b673aW644QYUFxdj48aNF7xW/bv6qWrKAEZEREDNJgYEBNR7D15AAiRAAiRAAhcjcDhdZe6mGqYvXhI6arbubfyNWb5x8hPd1h/N1EG8bI0iwAkcJ5kBTE5ORnh4OLZt24YBAwZUi2HevHnGPr99+/ZdVCBqz6Aycu+88w5uvfXWC16rloaffvrp8/6dBrBR3z8GkQAJkIDLE1BzMLuTc41SLcr0HUjLr8WkT2SQUa5ljMz2dQhp4fK8LAWABtDJDGBcXBxiYmKq9TF37lysWbMGe/fuvahm1BLxCy+8AGUkvby8OANoqW8Y+yEBEiABEjiPgFrK/TExqzpzNynbtB/dQ5ZyY2Qfn1GuRfb1hcrJHGyWJ0AD6CQGUGcJWP3XV9euXXH99ddj4cKFDVIZBdQgXLyYBEiABFyWQFFpOeIOZRgzfZskmSOroKSahY+csTtMMnaV6VMZvIE+ni7LyVYvzvHbSQygEoxKAunbt6+RBVzVVEmX8ePHXzQJJDY2FsOHD8euXbvQq1evBmmPAmoQLl5MAiRAAi5FIK+oFN/uU+VaTiJ2bxoKSsqr3z/I1xMju4cZJ3EM6doazcUEstmOAMdvJzKAVWVgli5daiwDL1++HCtWrDBq+qms3tmzZyMpKQmrV6+upTBVOubAgQPYvn17g5VHATUYGQNIgARIwKkJZOQXY1NCZRLHtoOZKJFCzVWtjSznjhHDp5I4rpRTOTzd3ZyahSO/HMdvJzKASmhq9u/55583CkGr2Ty1pDtkyBBDg6rwc2JiItSMX1VTyRtt27bFSy+9hOnTpzdYqxRQg5ExgARIgAScjsDxrMLq/Xw/Hc1CjWotuKR1C2NpV/1cJufvurFci0N8/hy/ncwA2lpVFJCtifN+JEACJGB/Amrv+P5UVa6lMnNXZfHWbJeK0VNLuyp7t3Oov/0fuCk/Qd5JwEMSYXyCLPoWHL9pALUERQFp4WMwCZAACTQZAhUyrfeL1OXbeNb0JWYWVj+7mtS7SpZ0jcxd+QkP8mky7+VwD5qfDhz9DjiyFUiUn4z9wPWSoNlvqkUfleM3DaCWoCggLXwMJgESIAGHJlAq+/e2H86sXt5NyzMdBODl4YbBnUMM0zcyOhSt/Lwd+l0c9uEKs8ToieFTZk+ZvvQ95zyquOuBjwCjn7HoK3D8pgHUEhQFpIWPwSRAAiTgcAROS6bu5v2Vmbtf70lFblFZ9TP6eXtghJRpUaZvqJRtUX9nayCB03K6ydE4k+FL/V06OOdAstCeQMfBQAf5iZLDHXxbNvAm9V/O8ZsGsH6VXOQKCkgLH4NJgARIwCEIZBeWiNlLM0zflgPpKCo1Ze6G+HlhtBRkVku7A6RAs7cHy7U06EMrzhPD970Yvi2VM3wnfxO/Z+Jr9BXSzWT4OgwCWoQ06BaNuZjjNw1gY3RTHUMBaeFjMAmQAAnYjcDJnCJsTKhM4th+OAvqdI6q1j7YxyjVMlaSOPpEBsOdmbvmf04lBcAxKatWtaSb/IsYPlP9Q6OjVp0rZ/eU2VN/+oeZ37+FruT4TQOoJSUKSAsfg0mABEjApgQOp6vM3coaffGS0FGzdW/jb8zyqezdHm0D0KyZ7D1jq59AqRxjd3xH5T4+NcOXtBOoKK0dF9yh0uh1lLJsyvQFtKu/XytfwfGbBlBLYhSQFj4GkwAJkIBVCahyLapEizJ86gi2A2n5te7XJzKoukZfh5AWVn0Wp+m8TBJhTvxoytJV/7vcdKyd8Z6BEWcN39l9fEHydwdrHL9pALUkSQFp4WMwCZAACVicgFrK/TExqzpzNylbZqjONg9Zyo2RfXxGuRbZ1xcqJ3Ow1UOgTMxd8s9nDZ/s4zv+A1BWVDvIv21tw6dm/Bx8BpXjNw2g1nefAtLCx2ASIAESsAiBotJyxB3KwH9/T8VXkrmbVWCakfKRM3aHylm7Y3uFYUS3MATKGbxsFyFQLlnPKfFi+MTsqX18aj9fqanmoRHZIrRG0obM8rXq5PCG79w35vhNA6j1e4AC0sLHYBIgARJoNIG8olJ8u6+yXEvs3jQUSPmWqhYkJm9k9zBjP9/gLq3h48XM3QuCrhBuKjO3qvCyytgtkczdms23lSlhQ+3jC+na5AwfDeD5CmgmeyTOKcDT6O+jywXSALrcR84XJgESsCOBjPxibEqoTOLYdjATJVKouaq1keXcMWL41PKuOpXD093Njk/qwLeuEGaq9l5Vlq6qyVecU/uBm8uxa1UZuqoeX+toyCHGDvxSDX80jt+cAWy4ampEUEBa+BhMAiRAAvUSOJ5VWL2f76ejWahRrQWXSOKGKtWiTN9lcv6uG8u1nM9TzfGkyekahuGTZd2j24DTp2pf5x1QWXDZyNSVn7BeYvice9aU4zcNYL2/fC52AQWkhY/BJEACJHAeAbUotT9VlWuprNGnsnhrtkvF6KmlXWX6Oof6sVzLuQSV4cs4YCq8rMqzFGbUvsrLD4iMMe3ja3MZ4O5ap5pw/KYB1Pr1SwFp4WMwCZAACRgEKmRaL/5EtiRxVJq+xExT0oGa1FNLukbmrvyEB/mQWk0CyvBlHTYt6aqZvvzU2ow8hFnk1WcNn+zha3eFGD7XTobh+E0DqPWLhALSwsdgEiABFyZQKvv3th/OrF7eTcuT+nJnm5eHGwZ3DjFM38joULTy83ZhUnW8+qlEU9KGmuHLTap9kbvwirjqbOFlWdIN7wt4eJFhDQIcv2kAtb4QFJAWPgaTAAm4GIHTkqm7eX9l5u7XUq4lt0hKjpxtft4eGN491DiCbWi31lB/ZztLIOeEyfCpbN2cY7XRuMlsXvsrTUu66n97ssbhxfTD8ZsGUOv3CwWkhY/BJEACLkAgu7BEzF6aYfq2HEhHUakpczfEzwujpSCzWtodIAWavT2cO/HA7I8776Sp8LIyfKeOnGP4xBy362MyfBH9AS9fs7vnhTJpmpuLwMBA5OTkICBAkmBcsLEMjMaHTgFpwGMoCZCA0xI4mVOErxLk+DUxfdsPZ0GdzlHV2gf7GEu74yR7t09kMNyZuSt79tIqz9KtKs2SKUkcNVszKcHSrrcpSzdC9vN5SyIHW6MJcPzmDGCjxaMCKSAtfAwmARJwIgKH01XmbmWNvvjj2bXerHsbf2OWT2Xv9mgbwMzdwqzaSRvpe89RgmS+tJXMXKMsiyRtqIzd5q45S2WtrwjHbxpALW1RQFr4GEwCJNCECahyLapES1W5FlW6pWbrExlkzPSpnw5Sr8+lm6q7pwouV522oQoxn9tU7b2qOnyqJp9PsEsjs/bLc/ymAdTSGAWkhY/BJEACTYyAWsr9MTGrOnM3Kft09Rt4yFJujOzjM8q1yL6+UDmZw2VbkdQuPCZHqlWdp5siR63hnEO3WnevYfgGyfm6ctwam80IcPymAdQSGwWkhY/BJEACTYBAUWk54g5lSI2+VHwlmbtZBSXVT+3j6Y6hXVvLaRxhGNEtDIFyBq9LtmKZ/Ty2vXJZV/0kx4vfM51NbDBp1cWUtKGOWfMLdUlUjvLSHL9pALW0SAFp4WMwCZCAgxLIKypF7L50I4kjdm8aCqR8S1UL9PHEqGh1EkcYBndpDR8vF8zcLZFC1cd3mPbxJf8s1axNJW0MVsEdTYWXleELaOugn7ZrPhbHbxpALeVTQFr4GEwCJOBABDLyi7EpoTKJY9vBTJRIoeaq1kaWc8ecPX5Nncrh6S5Zqa7USouAEz+aDF/ST0C5aSbUQBEYaZrhU+fpBrZ3JUJN7l05ftMAaomWAtLCx2ASIAE7EzieVVi9n++no1moUa0Fl0jihsrcVeVaLpPzd91cqVxLmZi7pJ1nDd+WSvNXJiawZvNvV9vwBXew86fJ2zeEAMdvGsCG6OW8aykgLXwMJgESsDEBlbmrsnWrMndVFm/NdqkYPbW0qxI5Oof6uU65lvLSyn17iWL2VKauWt4tNZ1HbDDyCzMlbahs3ZaXyAgq5VrYmiQBjt80gFrCpYC08DGYBEjABgQqZFov/kR2pen7/SQSM03GRk3qXdmhZWXmrhi/9sEucppEhexpTBHDV1WWRSVwlNQuYwPfEDF8sndPLed2kFp8IZLEQcNnA8Xa5hYcv2kAtZRGAWnhYzAJkICVCJTK/r3thzOrl3fT8oqr7+Tl4YbBnUMM0zcyOhSt/Lyt9BQO1G2F7GdM3WUyfKomX3Ht2U+j7l7UwMrCy2qGLzSahs+BPkJLPwrHbxpALU1RQFr4GEwCJGBBAqclU3fz/nTD9H0t5Vpyi0xZqX7eHhjePdRY3h3WLRTq707dlOFL32MyfOqYtaLap5PAO1AMnxRcNmb45EcVYnZzseQWpxbBxV+O4zcNoJb8KSAtfAwmARLQJJBdWCJmL80wfVsOpKOo1JS526qFl7GsqxI5BkiBZm8PJy7XInsbkbHfVHhZGb7CzNp0veTsXGX4qk7baCNHrbk5MRNNbTl7OMdvGkAtjVNAWvgYTAIk0AgCqblF2Kj288m5u9/LMq86naOqtQ/2qT5+rW9UMNydNXNXGb7MQ6bCy8rw5afWpukp+xkjrzadp9v2CsDdyWc+G6EnVw3h+E0DqKV9CkgLH4NJgATMJHA4XWXuVtboiz9eeymzW5i/nMShztwNQ4+2Ac6ZuasM36lEUx0+ZfjykmvT85Cj5yKuqkzYUMu67foAHl5mEuZlrkaA4zcNoJbmKSAtfAwmARK4AAFVrkWVaKkq16JKt9RsfSKDqmf6Oki9Pqds2cdrGD4pzZIjf6/Z3MXctb/StKQb3g/wdOHzh51SBNZ7KY7fNIBa6qKAtPAxmARIoAYBtZT7U2KWcfzaRpntS8o+Xf2vHrKUGyP7+NR+vjE9whAmJ3M4XctNMRVeVufpqhm/ms1Nlm/D+5oMX0R/MXw+ToeBL2QbAhy/ncwALl68GAsWLEBKSgp69uyJRYsWYfBgWQq4QCsuLsYzzzyDt956CydPnkT79u3x17/+FVOnTjVLgRSQWZh4EQmQwAUIFJWWI+5QhtTnS8UmydzNLDAdL+bj6Y6hXVvL8m4YRnQLQ6Cvp3NxzE+rPcOXebD2+zWTBI12vU1ZusrweUsiBxsJWIAAx28nMoBr167F5MmToUzgwIEDsWzZMqxcuRIJCQmIjJQzGuto48ePR2pqKp577jl07twZaWlpKCsrw4ABkilmRqOAzIDES0iABGoRyCsqRey+dGOmL3ZvGgqkfEtVC/TxNGrzqRp9Q7q0ho+XE2WpFkhWrprZUz+qAHPGvnMMn5RgUZm5VYWXVQJH8wCqhwSsQoDjtxMZwP79+6NPnz5YsmRJtViio6MxYcIEzJ8//zwBbdiwAbfffjsOHz6Mli1bNkpgFFCjsDGIBFyOQEZ+MTYlVCZxbDuYiRIp1FzV2shyrirXokzfVR1bwtPdSWrRnT4lZm+byfCl7T7/cw+71DTDp0q0+AS5nDb4wvYhwPHbSQxgSUkJfH19sW7dOkycOLFaTbNmzUJ8fDw2b958nsIefPBB7N+/H/369cOaNWvQokUL3HjjjXj22Wfh41P3vhK1ZKx+qpoSUEREBHJychAQwP9Stc/XmHclAcckcDyrEBuV6ZPj1346moUa1VpwiSRuqP184yR79zI5f9fNGcq1FOUAR7837eM7KSdvwFSixviUWsvpGlWFl9Uxa76N+49vx/zE+VRNiQANoJMYwOTkZISHh2Pbtm21lm/nzZuHVatWYd++c5YaRKXjxo1DbGwsRo0ahaeeegoZGRlQpnDEiBF4/fXX69TxnDlz8PTTT5/3bzSATelrz2clAesQUJm7Klu3KnNXZfHWbL3CAzBOTJ+a6esc6tf0y7UU5wHqDN2qJV11tu4Z08ym8e4hXU1JG1Fi+PxaWwc+eyWBBhKgAXQyAxgXF4eYmJhqGcydO9eY3du7d+950hgzZgy2bt1qJH8EBsqRQNLWr1+Pm2++GQUFBXXOAnIGsIHfMF5OAk5OoEKm9eJPZFeaPpnpS8wsrH5jNal3ZYeWhuFTS7ztg6UwcVNuJfJux8Xwqf17yvQl/SyGz7R/0Xi1lpeYCi+rGT7/Nk35jfnsTkyABtBJDGBjloDvvfdeY8bw4EFT5tmePXvQo0cPY2m4S5cu9UqfAqoXES8gAacjUCr797bLCRzK9H0lS7ypuaZtIV6yf29wlxDD9KlkjlZ+3k33/UuLgBM/mAzfiZ+AitLa7xMkCXZVhZfVEWuB4U33ffnkLkWA47eTGEClWpUE0rdvXyMLuKopM6cyfetKAlm+fDkeffRRI/PXz6+ytMDHH3+MSZMmIT8//4L7AGt+Qyggl/p9wZd1YQKnJVN38/504wg2Va4lt6ismoaftweGd1eZu2EY1i0U6u9NspWJkU3aaTJ8x8X8lZvMrfFOAWLwqs7SVX8GRzXJV+VDkwDHbycygFVlYJYuXWosAyuDt2LFCuzevRtRUVGYPXs2kpKSsHr1akP5yuSpLOGrr77a2Nen9gBOmzYNQ4cONeLMaRSQOZR4DQk0TQLZhSX4ek+aMdO35UA6ikpN+9tatfDCaCnIrI5gGyAFmr09mmC5lnKZzVPLuIlb5EeOVju2AygzFZ82PjU/WcKtStpQfwZ3lFFD1rbZSKCJE+D47UQGUGlRzf49//zzRiHoXr16YeHChRgyRM6FlDZlyhQkJiYaiR9VTe0NfPjhh42l4FatWuHWW281agJeKAv4XL1TQE38NwAfnwTOIZCaW2TM8qlzd7+XZV51OkdVax/sU338Wt+oYLg3tczdcpm1TPm10vCpfXwqgaO0oDYB35Aahk9+d7bqTMPHb4lTEuD47WQG0NYqpYBsTZz3IwHLEzicrjJ3K2v0xR/PrnWDbmH+xtKumunr0TagaWXuVkiChirFUpWle0xKtBTXzkyGT7As6UqyRtU+vtbdafgsLzH26IAEOH7TAGrJkgLSwsdgErALAVWuRZVoqSrXokq31Gx9IoOqZ/o6SL2+JtMqZIk6LcFk+I7Ksq6qzVezeUvFgw4DTfv4QntCihA2mVfkg5KApQhw/KYB1NISBaSFj8EkYDMCain3p8Qs4/i1jTLbl5Rt2uvmIUu5MbKPTxVmHiP7+sLkzIsBIAAAIABJREFUZI4m0cTIIl1KXKn9e0fO7uM7nVX70b38AXXCRtU+vjZy8oZbE9yv2CQ+ED5kUyLA8ZsGUEuvFJAWPgaTgFUJFJeVy7FrGVKfL9XI3M0sKKm+X3NPNwzrKpm7vcIwolsYAn09rfosFulcGb5MKVtlmD1Vi0+MX0F67a49ZcZSnaFbdZ5u28sB9yaalWwRaOyEBOomwPGbBlDru0EBaeFjMAlYnEBeUSli96UbM32xe9NQIOVbqlqgj6dRm0/V6BvSpTV8vBx8JkwZvlNHTGVZlOHLS6nNzENmKyP6mwxfeB8xfE3AzFr8k2eHJNAwAhy/aQAbpphzrqaAtPAxmAQsQiAjvxib1Jm7Yvq2HcxEiRRqrmphAd7V+/mu6tgSnlKo2aFb9jGT4VOZurknaj+uuxfQ/irTkm77foBHEy427dAfBh/OmQlw/KYB1NI3BaSFj8Ek0GgCJ04VVmfuqr19Naq14BJJ3FD7+VT27uXtgyTHwYHr1uUmnzV8Z0uzZB+tzcRNZvOUyVNFl1W2boSYP0+fRnNjIAmQQCUBjt80gFrfBQpICx+DScBsAipzV2XrVmXuqizemq1XeADG9miDcVKupXOon+OWa8lLPZulezZpI+tQbQbNZFlaLeNWnbahlne9mlAmstmfKC8kAfsS4PhNA6ilQApICx+DSeCiBCpkWi/+RHal6fv9JBIzC6uvV5N6V3ZoaSzvjpGZvvbBvo5JsyDDVJZFJW5k7D/H8MmStErUMAyfFF5WCRzekrnLRgIkYFUCHL9pALUERgFp4WMwCZxHoFT27+04rMq1pOAr2deXmms6i9ZL9u8N6hJiLO2Oig5DKz8H3PtWKGVYjm4z7eNTdflqNXGubXqZCi9HxgA+QVQCCZCAjQlw/KYB1JIcBaSFj8EkYBA4LZm6m/enG0ewfS2Zuzmn5Yzas83P2wPDu6vM3TAM6xYK9XeHaqrQ8tE40z6+k7/L45mOjzOeNbSHaUk3Soow+7Z0qFfgw5CAKxLg+E0DqKV7CkgLH4NdmEBOYalRm08t7245kI6iUlPmbqsWXhgtBZnV8WsDpECzt4cDlWspzhPDJ0eqGXX45EedrXvG9OzGRxrSzZSlqxI3Wsj5umwkQAIORYDjNw2gliApIC18DHYxAqm5RcYsnzp3d/vhTJTVSN0ND/IxEjjUnr6+UcFwd5TM3ZIC4Nh20z6+5F/E8JlqCxofYctONQyf7OXzD3OxT5avSwJNjwDHbxpALdVSQFr4GOwCBA6nq8zdypm++OPZtd64W5i/sbSrSrb0bBfgGJm7pXJE3PEfTIYvaSdQYVqSNl4gKMpUeFmduBHQzgU+Sb4iCTgXAY7fNIBaiqaAtPAx2AkJqHItqkRLVbkWVbqlZusdGYRxRuZuG3SUen12b2WSZHLiJ5PhOyHmr9x0ZJzxfAHtTTN8yvAFRdr9sfkAJEACegQ4ftMAaimIAtLCx2AnIVAuS7mqGLM6fm2jzPYlZcss2tnmIUu5MbKPTxm+MbKvLyxAji6zZysTc5f8sylL9/gOoKyo9hP5tzUlbajyLMEd5DelAxeTtidP3psEmigBjt80gFrSpYC08DG4CRMoLiuXY9cypD5fqpHMkVlgmjVr7umGoV1bG/v5RnYPQ6CvHc+mLS+TRI14MXyq8LIkbaj9fKWmeoLGR9CidQ3DJ7X4WsmePhq+JqxOPjoJ1E+A4zcNYP0qucgVFJAWPgY3MQL5xWX4Vsq0qOVd9WeBlG+paoE+nhgZrcq1tMGQLq3h42WnzN0KeaaTv5lm+FTGbolk7tZsPlKGRWXnqsLLaoavtWTt0vA1MTXycUlAjwDHbxpALQVRQFr4GNwECGTkF2OTFGRWpm/bwUyUSKHmqhYW4C3LupXHr13VsSU8pVCzzVuFPE/a7hqGT4owq9p8NVvzQCBKGT51nq78qLp8bnZ4VpvD4Q1JgAQuRIDjNw2g1reDAtLCx2AHJXDiVGF15q7a21ejWouRuKFm+VT27uXtg8RH2XhvnCSZIG2P6TxdderG6VO1SXoHiOEbYFrWDZOTN9zsNCPpoJ8xH4sEXJ0Ax28aQK3vAAWkhY/BDkJAZe4eSMvHBjlvV830qSzemq1XeADGykyfKszcJdTPtuValOHLOCCGT/bwHVHFl78DCuV83ZrNU7KJo+RINeM8XflpI2frujvYiSEO8lnzMUiABCoJcPymAdT6LlBAWvgYbEcCyvSpunxVmbtHMqTg8dmmJvWu7NDSmOkbIzN97YN9bfekyvBlHTaVZVGGL/9k7ft7+ACR/c8aPtnH1663GD47JprYjg7vRAIkYCECHL9pALWkRAFp4WOwnQgczyrEjDU7sSfFNNPnJfv3BnUJMZZ2R0WHoZWft+2f7tRR4L3Jlcer1Wzu8iwRV5lm+ML7Ah52eD7bE+EdSYAErESA4zcNoJa0KCAtfAy2AwE183fnih34Xo5i8/JwM2rzqSSOYd1C4edtx2XTEinN8tpoIPV32a8nz9FeDJ+RtCHJG+p/e9q5fqAdPivekgRIwHoEOH7TAGqpiwLSwsdgOxBY99Nx/On936Bq9X328GB0lj19dm9q2Xf9dGDXusqafNO/4Wkbdv9Q+AAk4NwEOH7TAGopnALSwsdgGxPIlJIuI1/cjOzCUjx5TXc8MFQKHjtC274U2PAX+W0kmbr3fiqzfgMd4an4DCRAAk5MgOM3DaCWvCkgLXwMtjGBx9bG48NfkhDdNgCf/GGgfer2nfvOiVLGZfWNQIWc2DF2PhDzoI2p8HYkQAKuSIDjNw2glu4pIC18DLYhga0H0jH5tR+MAy8+fHAgrogIsuHdL3Cr3BRgmWTxFqQBvW4GblrJEzns/6nwCUjAJQhw/KYB1BI6BaSFj8E2InBajmwbu2gLjkn275QBHTDnxp42uvNFblMmZwe/eR1w4gc5mUOeZ9pXgJfU82MjARIgARsQ4PhNA6glMwpICx+DbUTgnxv2YknsIbQNbI6vHh9q32zfqnf+/AngR5nxU8e0zYgFWl5iIxq8DQmQAAmwELTSQDMpCyEpeGyNIUAD2BhqjLElAVXr7/qXv0O5nOe24p5+GC1lX+zefnkb+Fjt9ZP16DvfA7qOsfsj8QFIgARciwDHbxpALcVTQFr4GGxlAsr03bQkzjjxY5yc6rF0shRQtndLjpd6f2L4youBYf8rP5L9y0YCJEACNibA8ZsGUEtyFJAWPgZbmcCquET83ye74S8FntXSbxtZArZrK8gElg8Dco7JrN844PZ3peizm10fiTcnARJwTQIcv2kAtZRPAWnhY7AVCaTknMboF7cgv7gMz47vickxHax4NzO6rigH3poEHI4FgjtW7vvzcYBMZDMenZeQAAk4HwGO3zSAWqqmgLTwMdiKBGas/gkbE1LROzIIHzwwQCbaZL+dPdumOcB3C+VIN1/J+N0EhDlAJrI9efDeJEACdiXA8dvJDODixYuxYMECpKSkoGfPnli0aBEGD5bzROtosbGxGD58+Hn/smfPHnTv3t0sYVJAZmHiRTYmsOH3k3jgrZ3wENP32SOD0L1NgI2f4JzbJXwCvDe58v+86TXgUqn5x0YCJEACdiTA8duJDODatWsxefJkKBM4cOBALFu2DCtXrkRCQgIiIyPPk1mVAdy3bx8CAkwDZOvWreHuLkdSmdEoIDMg8RKbEsgrKjWWfk/mFuHBYZ3w53Hm/ceM1R4yfT+wQv5DqyRfTvn4g5z2Mddqt2LHJEACJGAuAY7fTmQA+/fvjz59+mDJkiXVn390dDQmTJiA+fPliKlzWpUBPHXqFIKCGrcXiQIy96vG62xF4P8+/h2rvj+KqFa++O+jQ9Dc07z/mLHK8xXlAitHAhliAjvITPzkjwB3D6vcip2SAAmQQEMIcPx2EgNYUlICX19frFu3DhMnTqzWwKxZsxAfH4/Nmzdf0AB26NABRUVF6NGjB/72t7/VuSxcFVxcXAz1U9WUgCIiIpCTk1NrFrEhIuS1JGApAj8fO2WUfVGVPd+e1h8DO4dYquuG96MeYu3dwN7PAP92wMwtgF/rhvfDCBIgARKwAgEaQCcxgMnJyQgPD8e2bdswYMCAaqnMmzcPq1atglrmPbep/2/Lli3o27evYerWrFmDpUuXQs0MDhki55PW0ebMmYOnn376vH+hAbTCt5NdNohAaXkFbpCCz3tP5mFSn3C8eOsVDYq3+MVbXwS+lu+Kuxdw35dA+34WvwU7JAESIIHGEqABdDIDGBcXh5iYmGo9zJ071zB2e/fuNUsjN9xwA5o1a4ZPPpFN63U0zgCahZEX2YHA4tiDeH7DPgT7euLrJ4ahZQsxXvZqh76Rki83AWcqgOsXAf3us9eT8L4kQAIkUCcBGkAnMYCNWQKuSxHKML711ltQmcDmNArIHEq8xtoEjmYWYMzCLSguq8ALt1yOm/q2t/YtL9z/qaOVxZ5PZwG9JfP3xpflt4ydS9DYjwbvTAIk4KAEOH47iQFU+lJJIGo5V2UBVzW1r2/8+PF1JoHUpcmbb74ZWVlZ+OYbmcEwo1FAZkDiJVYloI7yvuf1H7D1QIbs+WuFt+7vb8xi26WVngZeHwuk/Aq06y1Lvxuk7p+dTx+xCwjelARIwNEJcPx2IgNYVQZG7eNTy8DLly/HihUrsHv3bkRFRWH27NlISkrC6tWrDV2qGoEqAUTVC1QziGrm7x//+Ac++OADTJokJxaY0SggMyDxEqsS+PCXE3hs7a/w9nAzsn47hLSw6v0u2LlK+vjoQeDXdwDfVnLShyReBUXY51l4VxIgARKohwDHbycygOqzVrN/zz//vFEIulevXli4cGF1QseUKVOQmJhoJHmopq5TJlGZQh8fH8MIKpN47bXXmv3FoYDMRsULrUDgVEEJRr64GVny55/GdsNDwztb4S5mdvnjSuDzJ+Q3ipztq8q9XDLUzEBeRgIkQAK2J8Dx28kMoK0lRAHZmjjvV5PAH9f9ivd3nkC3MH98+vAgeMksoF3asR3Am9cBFaXA6GeAgbPs8hi8KQmQAAmYS4DjNw2guVqp8zoKSAsfgzUIxB3MwJ0rdxj5Fe/LWb99o4I1etMIzUsFlknZpPyTQI/xwC2rmPShgZOhJEACtiHA8ZsGUEtpFJAWPgY3kkBRaTmueWkrjmQUYPLVUXh2Qq9G9qQZVi4zfqtuBI7FAa3lyLlpmwBvf81OGU4CJEAC1ifA8ZsGUEtlFJAWPgY3ksALG/fh5W8OItTfG5ueGIqA5p6N7Ekz7MsngR1y9KK3nKU9XTLnQ7podshwEiABErANAY7fNIBaSqOAtPAxuBEE9qfm4VqZ/SurOIOld/fBuF5tG9GLBUJ+ew9YP72yo9sl87e77AFkIwESIIEmQoDjNw2gllQpIC18DG4ggQoxfbcs+x47j57CqOgwrLinr31q/p3cBawcDZRJ3b8hfwJG/K2Bb8LLSYAESMC+BDh+0wBqKZAC0sLH4AYSeGv7Ufzto9/RwssdXz0+FO2CfBrYgwUuL5QTPtRJH9ly4kenkcBd6wA3dwt0zC5IgARIwHYEOH7TAGqpjQLSwsfgBhBIzS3CqBc2I6+4DP93Qw/cN7BjA6ItdGmFnO37zq3Awa+kyHOUFHuOlaLPLS3UObshARIgAdsR4PhNA6ilNgpICx+DG0Dgwbd34otdJ3F5+0Csf3Ag3N3scNzbt/OAzf8EPOR4t/vFBLa9rAFvwEtJgARIwHEIcPymAdRSIwWkhY/BZhLYlJCKaat/MkzfJ38YiJ7tAs2MtOBl+74E3r29ssOJy4HLb7Ng5+yKBEiABGxLgOM3DaCW4iggLXwMNoNAgSz5jpbj3pJzijBzyCWYfW20GVEWviTzUOW+v+Jc4KqZwLXPW/gG7I4ESIAEbEuA4zcNoJbiKCAtfAw2g8Aznybg9W1HENHSB/99dAh8vTzMiLLgJcX5kvE7CkjfA0TGAPd8IkvAXha8AbsiARIgAdsT4PhNA6ilOgpICx+D6yHw24lsTHh1G6T6C1ZNvQpDu7a2LbMzcuP37wN2fwj4hQEztwD+bWz7DLwbCZAACViBAMdvGkAtWVFAWvgYfBECZeUVuPGVbUhIycX4K9rhpdt7255X3MvARqnx5yazjlM+lxnAq23/DLwjCZAACViBAMdvGkAtWVFAWvgYfBECK7Ycxtwv9iDQxxNfy3FvIX7etuV1RGb7Vk8AzpTLnr9/yd6/s6d+2PYpeDcSIAESsAoBjt80gFrCooC08DH4AgSOZxVizMItOF1ajudvugy3XhlhW1Y5J4BlQ4HCDMn2vQOYIOf9NrND2RnbvjXvRgIk4EIEOH7TAGrJnQLSwsfgOgickX13U974EZv3p6N/x5b4z4yrbXvcW2kR8MY1QPLPQJtLK+v9edrhxBGqgwRIgASsSIDjNw2glrwoIC18DK6DwCe/JuORd3+Bl7sbvnx0MDq19rMtp08eAX5eBfgEV570EdzBtvfn3UiABEjABgQ4ftMAasmMAtLCx+BzCOQUlmLki7HIyC/BY6O6YtaoLrZltPNN4NNZck9Z7r37A6CznPXLRgIkQAJOSIDjNw2glqwpIC18DD6HwJMf/Ib//HgcnUP98Pkjg+Dt4W47Rid2ytLvOKC8BBjxd2DIH213b96JBEiABGxMgOM3DaCW5CggLXwMrkFgx+FM3LZ8u/H/rHsgBld2aGk7PvnpctKHJH3kJgHdrwduXSOlX9xsd3/eiQRIgARsTIDjNw2gluQoIC18DD5LoLisHNe8tBWH0wtwx1WRmD9Jki9s1crLgDVS7iVxK9CqMzD9W6B5gK3uzvuQAAmQgF0IcPymAdQSHgWkhY/BZwks2rQfizYdMGr9ff34UAT6etqOjSr0rAo+e0myybSvgdDutrs370QCJEACdiLA8ZsGUEt6FJAWPgYLgYNp+bhWZv9K5OSPl+/ojRsub2c7Lr+vrzzqTbVbVwM9xtvu3rwTCZAACdiRAMdvGkAt+VFAWvhcPrhCDvm9Xfb9/ZCYheHdWuP1KVfaruZfagKwchRQWgAMlMzf0c+4/OdBACRAAq5DgOM3DaCW2ikgLXwuH/yfH47hyfW74OPpjo2PDUFES1/bMCnKkaSP4UDWIaCjJH/cLTOB7nLeLxsJkAAJuAgBjt80gFpSp4C08Ll0cHpeMUa+EIvcojL87bpoTBt8iW14VFQAa+8C9n0BBMoRczNigRYhtrk370ICJEACDkKA4zcNoJYUKSAtfC4drE77UKd+9AoPwEcPDoSHnPxhk7Z5AfDtczLj5w1M3QCE97HJbXkTEiABEnAkAhy/aQC19EgBaeFz2eDYfWnGeb9ucuDGxw8NwqXtA23D4oCc6/v2LXKvM8D4V4Hed9vmvrwLCZAACTgYAY7fNIBakqSAtPC5ZHBhSRnGLNyCE6dO4/5BHfH363vYhkPWYdn3NwxQ+//6TQWuX2ib+/IuJEACJOCABDh+0wBqyZIC0sLnksHzvtiD5VsOIzzIx0j8aOFtg+SLkkLgtdFA6u+y5NsPuE/2/3nIEjAbCZAACbgoAY7fNIBa0qeAtPC5XPDvSTmy8roN5VL+5fUp/TCie5j1GZyR5d71M4Bd70myR2tJ+tgsyR/h1r8v70ACJEACDkyA4zcNoJY8KSAtfC4VrEzfxMXb8NuJHFx3WVu8eqeNki92LAO+/LN8092Bez8BOgxyKe58WRIgARKoiwDHbxpArW8GBaSFz6WCX//uCJ75LAH+zT3w9RNDEerf3PrvfzQOWHUDUCHn/Y6dD8Q8aP178g4kQAIk0AQIcPymAdSSKQWkhc9lgpOyT2P0i5tRWFKOeRMvxZ39I63/7rkpwLIhQEEa0Otm4KaV8m2XtGM2EiABEiABcPx2MgO4ePFiLFiwACkpKejZsycWLVqEwYMH1yv1bdu2YejQoejVqxfi4+Prvb7qAgrIbFQue+EZ2YM3bdVP+HpvGvpFBeO9mTFwU/VfrNnKSoA3rwNO/ACE9gSmSfkXrxbWvCP7JgESIIEmRYDjtxMZwLVr12Ly5MlQJnDgwIFYtmwZVq5ciYSEBERGXnjGJScnB3369EHnzp2RmppKA9ikvsKO/7Bf7ErBg2//DE/3ZvjikcHoEuZv/Yf+/AngR5nx85b6gjO+BVp1sv49eQcSIAESaEIEaACdyAD279/fMHJLliyplmB0dDQmTJiA+fNl/9MF2u23344uXbrA3d0dH330EQ1gE/oCO/qj5pwuNZZ+0+TYt0dGdMbjY7pZ/5Hj3wE++p/K+9wpmb9dx1r/nrwDCZAACTQxAjSATmIAS0pK4Ovri3Xr1mHixInVMpw1a5Zh6DZvltIXdbQ33njDmDH8/vvv8dxzz9VrAIuLi6F+qpoSUEREBNQsYkBAQBOTPx/X2gT++uEuvL3jGC4JaYEvZg1Gc0/JxLVmS5btC6+L4SsrAobNlp8nrXk39k0CJEACTZYADaCTGMDk5GSEh4dD7eUbMGBAtSDnzZuHVatWYd++feeJ9MCBAxg0aBC2bt2Krl27Ys6cOfUaQHXN008/fV5fNIBN9neA1R78p8Qs3Lz0e6P/d6dfjZhOrax2L6PjwixJ+hgK5ByTWb9xwO3vQjYbWvee7J0ESIAEmigBGkAnM4BxcXGIiYmpluPcuXOxZs0a7N27t5ZEy8vLcfXVV+P+++/HAw88YPybOQaQM4BN9Jtu48cuKavAdf/eigNp+bilb3ssuOVy6z5BRTnw1k3AYdnvF9xR9v3FAj5B1r0neycBEiCBJkyABtBJDGBDl4Czs7MRHBxs7PurahUVFVAZm+r/27hxI0aMGFGvtCmgehG55AWvfHMA/9q4H61aeGHT40MRLH9atW2SWenvXgQ8fSXjdxMQJpm/bCRAAiRAAhckwPHbSQyg+oRVEkjfvn2NPX1VrUePHhg/fvx5SSDK7Kns4JpNxX3zzTd4//330bFjR7RoUX/ZDAqIv13OJXAkowBjF22BmgVcdNsVmNDbyseuJcjpHu9NrnyMm14DLpWaf2wkQAIkQAIXJcDx24kMYFUZmKVLlxrLwMuXL8eKFSuwe/duREVFYfbs2UhKSsLq1avrFIU5S8DnBlJA/A1Tk4CaQb5zxQ58fzgTg7uEYPXUq6T2shVr/qXvB1YMB0rygasfAsbN4wdCAiRAAiRgBgGO305kANXnrWbxnn/+eaMQtCrqvHDhQgwZIqchSJsyZQoSExMRGxtLA2jGl4OXNJzA+ztP4I/rfpVsXzdsfHQoIlvJkqy1WnGemD/ZppAhJjBKzve95yPA3dNad2O/JEACJOBUBGgAncwA2lqdFJCtiTvu/TLzizFSav5lF5biyWu644GhViy+LDONxrLvnk8B/3bATClz5BfquHD4ZCRAAiTgYAQ4ftMAakmSAtLC51TBj6+Nx/pfktC9jT8+fXiQnPxhxRIs3y0ENs2RGT9JLrnvS6B9P6diyZchARIgAWsT4PhNA6ilMQpIC5/TBG89kI7Jr/0g+/2ADx8ciCsirFiC5dA3lSVfzlQA1y8C+t3nNBz5IiRAAiRgKwIcv2kAtbRGAWnhc4rg0yXlRtbvsaxCTBnQAXNutGIJllNHgeXDgNNS9Ln33cCNr8g32IpJJk7xCfElSIAESOB8Ahy/aQC1vhcUkBY+pwj+54a9WBJ7CG0Dm+Mrqfnn5+1hnfcqPV15zFvKr0C73rL0u0Hq/jW3zr3YKwmQAAk4OQGO3zSAWhKngLTwNfngPSm5uOHl71BWcQbLJ/fFmJ5trPNOKunjYynzEv824CtHys2QpI+gCOvci72SAAmQgAsQ4PhNA6glcwpIC1+TDi4X03fTkjjEH8/GODF+S8UAWq39KAWeP39cvq2SWDL5Q+CSYVa7FTsmARIgAVcgwPGbBlBL5xSQFr4mHbz6+0Q89fFuY8lXHffWRpaArdKO/wC8cS1QUQqMfgYYOMsqt2GnJEACJOBKBDh+0wBq6Z0C0sLXZINP5hRhlNT8yy8uw7Pje2JyTAfrvEteqiR9DAXyUoAe44FbVjHpwzqk2SsJkICLEeD4TQOoJXkKSAtfkw2eueYn/Hd3KnpHBuGDBwbAzc0KmbjlMuO36kbgWBwQ0g2Y/jXg7d9kmfHBSYAESMCRCHD8pgHU0iMFpIWvSQb/d/dJzFyzEx5i+j57ZJAUfg6wznt8+SSwYwngJaZvxrdiArtY5z7slQRIgARckADHbxpALdlTQFr4mlxwXlEpRr+4BSdzi/DgsE7487ju1nmH39YB66dV9n2bZP5GX2+d+7BXEiABEnBRAhy/aQC1pE8BaeFrcsFzPtmNN+MSEdXKF/99dAiae7pb/h1O/g6sHAWUSd2/wX8ERv7d8vdgjyRAAiTg4gQ4ftMAan0FKCAtfE0q+JdjpzBJyr6oknxv3d8fg7qEWP75T5+qPOnjVCLQaSRwl8wEulnBZFr+ydkjCZAACTQpAhy/aQC1BEsBaeFrMsGl5RVGwee9J/MwqXc4XrztCss/e4Wc7fvOrcDBr6TIc5Ts+4uVos8tLX8f9kgCJEACJACO3zSAWl8DCkgLX5MJVke9qSPfgn09jZp/rfy8Lf/s384DNv8T8JB6gvdvBNpebvl7sEcSIAESIAGDAMdvGkCtrwIFpIWvSQQfzSzA2EVbUFRagX/dcjlu7tve8s+970vg3dsr+524DLj87P+2/J3YIwmQAAmQAA2goYFmZ6RRDY0jQAPYOG5NJUp9Ne55/QdsPZCBAZ1a4e1p/dGsmYVr/mUekn1/w4HiHOCqGcC1C5oKHj4nCZAACTRZAhy/aQC1xEsBaeFz+OCPfknCo2vj4eXhZmT9dgxpYdlnLs4HXhsNpCUAEVcD934qS8Belr0HeyMBEiABEjiPAMdvGkCtrwUFpIXPoYNPFZRgpBz3liV//mlsNzw0vLNXuaL6AAAgAElEQVRln1dNvL8/Fdi9HvALA2ZuAfzbWPYe7I0ESIAESKBOAhy/aQC1vhoUkBY+hw7+07pfsW7nCXQN88NnDw82ZgEt2uJeATb+Vcq8eABTPgciZQaQjQRIgARIwCYEOH7TAGoJjQLSwuewwXGHMnDnih2y308m6eSs375RwZZ91iMy27d6AnCmHLhG9vz1l71/bCRAAiRAAjYjwPGbBlBLbBSQFj6HDC4qLcc1L23FkYwC3H11JJ6bcKllnzPnBLBsKFCYAVwm2b4Tl8q30MKJJZZ9YvZGAiRAAk5HgOM3DaCWqCkgLXwOGfzCxn14+ZuDCPX3xqYnhiKguaflnrOsGHjjGiBpJ9BGjOVUqffn5Wu5/tkTCZAACZCAWQQ4ftMAmiWUC11EAWnhc7jg/al5uO7fW1FafgZL7uqDay5ta9ln/HQWsPNNoHmQJH1sBoI7WLZ/9kYCJEACJGAWAY7fNIBmCYUGUAtTkwiuqDiDW5Z9j51HT2FUdBhW3NPXsjX/dq4CPn1EWMhy793vA51HNQkufEgSIAEScEYCNIA0gFq6poC08DlU8Ns7juKvH/6OFl7u+EqOe2sX5GO55zshS75vjAPKS4ARfweG/NFyfbMnEiABEiCBBhPg+E0D2GDR1AyggLTwOUxwWm6RUfMvr6gMT13fA1MHdbTcs+Wny0kfkvSRmwR0uw647S0p/WLhkjKWe1r2RAIkQAIuQYDjNw2gltApIC18DhP80Ns/4/NdKbisfSA+fHAg3N0slJVbXgaskXIviVuBVlJIevo3sv8v0GHemw9CAiRAAq5KgOM3DaCW9ikgLXwOEfz1nlTcv+onw/R98oeB6NnOggZt49+AuJcBTzlCTpm/0O4O8c58CBIgARJwdQIcv2kAtb4DFJAWPrsHFxSXYczCLUjKPo2ZQy7B7GujLfdMuz8E1k2p7O8WSQDpKTOBbCRAAiRAAg5BgOM3DaCWECkgLXx2D372swS89t0RtA/2wcbHhsDXS45ls0RL2wOsGAmUFgADpfTL6Gcs0Sv7IAESIAESsBABjt80gFpSooC08Nk1eNeJHIx/9TtI9Re8ed+VGNYt1DLPU5QjSR/DgaxDQEdJ/rh7PeBuIWNpmSdkLyRAAiTg8gQ4ftMAan0JKCAtfHYLLiuvEPO3DbuTc3Hj5e3w7zt6W+ZZKiqAtXcB+74AAtpXFntuEWKZvtkLCZAACZCAxQhw/KYB1BITBaSFz27BK7Ycxtwv9iDQxxObpOZfazn2zSJt8wLg2+dkxk/6m7oBCO9jkW7ZCQmQAAmQgGUJcPx2MgO4ePFiLFiwACkpKejZsycWLVqEwYMH16ma7777Dn/5y1+wd+9eFBYWIioqCjNnzsRjjz1mtsooILNROcyFx7MKjcSP06Xl+OdNl+K2KyMt82wHvgLevkX6kjXlG18B+ky2TL/shQRIgARIwOIEOH47kQFcu3YtJk+eDGUCBw4ciGXLlmHlypVISEhAZOT5g/wvv/ximL/LLrsMLVq0gDKEygAuXLgQM2bMMEtsFJBZmBzmojNnzuC+N39E7L509O/YEv+ZcbVljnvLOiL7/oYBRdlA3/uAGxY5zDvzQUiABEiABM4nwPHbiQxg//790adPHyxZsqT6k46OjsaECRMwf/58s/Q/adIkwwyuWbPGrOspILMwOcxFn/6ajIff/QVe7m748tHB6NTaT//ZSgqB18YAqbtkybcfcJ/s//Ow0JKy/tOxBxIgARIggToIcPx2EgNYUlICX19frFu3DhMnTqz+qGfNmoX4+Hhs3iyb8etpakbwmmuuwXPPPYdp06bVeXVxcTHUT1VTAoqIiEBOTg4CAgLquwX/3Y4EcgpL5bi3WGTkl+CxUV0xa1QX/aeRGUWsl9niXe9JskdrYIboLDBcv1/2QAIkQAIkYFUCNIBOYgCTk5MRHh6Obdu2YcCAAdWimTdvHlatWoV9+/ZdUEjt27dHeno6ysrKMGfOHPz973+/4LXq359++unz/p0G0KrfU4t0Pnv9b3j3h+My69cCX8waDG8Pd/1+dywDvvyzfIukr3s/AToM0u+TPZAACZAACVidAA2gkxnAuLg4xMTEVAtn7ty5xnKu2ut3oXbkyBHk5+dj+/btePLJJ/HKK6/gjjvu4Ayg1b9+trvBD0eycOuy740bvjczBlfJ/j/tdjQOWHUDUCHn/Y6dB8Q8pN0lOyABEiABErANARpAJzGAllgCVpJTy7/KMF5sxrCmNCkg23xRde5SXFaOa1/aikPpBbjjqgjMn3SZTneVsbkpwLIhQEEa0Osm4KbX5JvUTL9f9kACJEACJGATAhy/ncQAKrWoJJC+ffsaWcBVrUePHhg/frzZSSDPPvssXnvtNSQmJpolQArILEx2vWjRpv1YtOkAQvy88bXU/Av09dR7nrISmfm7Hji+AwjtAUzbBHi10OuT0SRAAiRAAjYlwPHbiQxgVRmYpUuXGsvAy5cvx4oVK7B7926jxt/s2bORlJSE1atXGyJ79dVXjfIw3bt3N/6uysA8+uijePjhh42ZQHMaBWQOJftdczAt35j9K5GTP16W0z5ukFM/tNvnfwR+XAF4B0rSx7dAq07aXbIDEiABEiAB2xLg+O1EBlBJR83+Pf/880Yh6F69ehk1/YYMkaU6aVOmTDFm9mJjY42/v/zyy0atQLUH0MPDA506dcL06dONWoBubm5mKZECMguTXS6qkEN+b1+xHWr/37BurfHGlCv1a/7FvwN89D+V73OnZP52HWuXd+NNSYAESIAE9Ahw/HYyA6gnh4ZHU0ANZ2ariLU/HsNfPtgFH093bHxsCCJa+urdOjkeeF0MX1kRMGy2/Dyp1x+jSYAESIAE7EaA4zcNoJb4KCAtfFYLTs8rxqgXNyPndCn+em00pg+5RO9ehVmS9DEUyDkGdBETeMd/INPEen0ymgRIgARIwG4EOH7TAGqJjwLSwme14EfktI9P5NSPnu0C8PFDA+EhJ380ulWUA29Jpu9h2e8X3LFy359PcKO7YyAJkAAJkID9CXD8pgHUUiEFpIXPKsGx+9Iw5Y0f4SZVWT5+aBAubS/JGjrt62eArS8AnrKEfP9XQJteOr0xlgRIgARIwAEIcPymAdSSIQWkhc/iwYUlZRizcAtOnDqN+wd1xN+vlzItOm3Pp8Dauyt7ULX+Lr1ZpzfGkgAJkAAJOAgBjt80gFpSpIC08Fk8eP4Xe7Bsy2GEB/kYiR8tvD0af4/0/cCKEUBJHnC1nPIxTk77YCMBEiABEnAKAhy/aQC1hEwBaeGzaPDu5Bzc+Mo2lEv5l9en9MOI7mGN779YTJ8yfxliAqPkfN97PgLcNQtIN/5pGEkCJEACJGBhAhy/aQC1JEUBaeGzWLAyfRMXb8NvJ3Jw3aVt8epdfRrf95kzcmDwZEAt//pL4eiZmwG/0Mb3x0gSIAESIAGHI8DxmwZQS5QUkBY+iwW//t0RPPNZAvybexjHvYUGNG98398tBDbNkTIvMuN335dAxJWN74uRJEACJEACDkmA4zcNoJYwKSAtfBYJTso+jdFS86+wpBxzJ/bCXf2jGt/vISnx8tYk4EwFcL0YwX5TG98XI0mABEiABByWAMdvGkAtcVJAWvi0g8/Icu301T9h05409IsKxnszY6Q+s9R/aUzLliLPqtjzaSn63Fsyf298Rb4djeyrMfdnDAmQAAmQgM0IcPymAdQSGwWkhU87+MtdKfift3+Gp3szfPHIYHQJ829cn6WnK495S/kVaNdbln43SN0/jWXkxj0Fo0iABEiABGxEgOM3DaCW1CggLXxawblFpRj1wmakybFvj4zojMfHdGtcfyrp42Mp8xL/NuDbSk76kKSPoIjG9cUoEiABEiCBJkGA4zcNoJZQKSAtfFrBf/toF97afgyXhLTAF7MGo7mne+P6+1EKPH/+uHwT5Li4yR8ClwxrXD+MIgESIAESaDIEOH7TAGqJlQLSwtfo4J1Hs3Dz0u+hJu/emd4fAzqFNK6v4z8Ab1wLVJQCo54GBj3auH4YRQIkQAIk0KQIcPymAdQSLAWkha9RwSVlFbj+5a3Yn5qPW/q2x4JbLm9UP8hPk6SPIUBeChB9I3DraiZ9NI4ko0iABEigyRHg+E0DqCVaCkgLX6OCX/32IBb8dx9atvAyav4Fy58NbuUy47d6PHB0GxAiewenfw14NzKBpME3ZwAJkAAJkIC9CXD8pgHU0iAFpIWvwcFHMgowdtEWqFnARbddgQm9wxvchxGwYTawfTHgJaZvhtT+C+nSuH4YRQIkQAIk0CQJcPymAdQSLgWkha9Bwarm310rdyDuUCYGdwnB6qlXSZm+RtTp+20dsH5a5b1vk8zf6Osb9By8mARIgARIoOkT4PhNA6ilYgpIC1+Dgt/feQJ/XPerZPu6YeOjQxHZyrdB8cbFJ38HVo4CyqTu3+AngJFPNbwPRpAACZAACTR5Ahy/aQC1REwBaeEzOzgzvxij5Li3U4Wl+Mu47vifYZ3Mjq2+8PQpYPkw4FQi0GkEcNf7ct5vI0vHNPzujCABEiABEnAgAhy/aQC15EgBaeEzO/jxtfFY/0sSurfxx6cPD5KTP6RmX0NahZzt++5twIGNUuQ5srLYs2/LhvTAa0mABEiABJyIAMdvGkAtOVNAWvjMCv7uQAbufm2HcSzvhw8OxBURQWbF1bro2/nA5n8AHnK82/1iAts2snRMw+/MCBIgARIgAQckwPGbBlBLlhSQFr56g4tKy42s36OZhZgyoAPm3Niz3pjzLtgn5/qq2T/VJi4DLr+94X0wggRIgARIwKkIcPymAdQSNAWkha/e4Oc37MXi2ENoE9AcXz0+BP7NPeuNqXVB5iHZ9zccKM4BrpoBXLugYfG8mgRIgARIwCkJcPymAdQSNgWkhe+iwXtP5uL6f3+HsoozWD65L8b0bNOwmxXnA6+NBtISgIj+wL2fyRJwI4pGN+yuvJoESIAESKAJEOD4TQOoJVMKSAvfBYPLxfTdtCQO8cezMbZnGJZN7tewG6lDgt+fCuxeD/iFVSZ9BLRtWB+8mgRIgARIwGkJcPymAdQSNwWkhe+Cwau/T8RTH++Gn7cHNslxb20CJXmjIe37V4H//q+UefGonPmLimlINK8lARIgARJwcgIcv2kAtSROAWnhqzP4ZE6RUfMvv7gMz4zviXtiOjTsJke2Vp7ze6YcuEb2/PWXvX9sJEACJEACJFCDAMdvGkCtLwQFpIWvzuAH1uzEht0n0TsyCO8/MADubg047i0nCVg2BCjMAC6TbN+JS0XhDYi3/OuwRxIgARIgAQckwPGbBlBLlhSQFr7zgjeK8ZshBtBDTN9njwySws8B5t+grBh44xogaSfQ5lJgqtT782rEcXHm35FXkgAJkAAJNFECHL9pALWkSwFp4asVnFdUitEvbsHJ3CLjqDd15FuD2qezgJ1vAs2lUPRMSfoI7tCgcF5MAiRAAiTgOgQ4ftMAaqmdAtLCVyt4zie78WZcIqJa+eK/jw5Bc88GnNP782rgk4elP1nuVWf8dhlluQdjTyRAAiRAAk5HgOM3DaCWqCkgLXzVwarcy8TF26Cqt7x1f38M6hJifsdqyfd1WfotlyXgEX8DhvzJ/FheSQIkQAIk4JIEOH7TAGoJnwLSwmcEl5ZX4IaXv8Pek3mY1DscL952hfmdFkiyx7KhQO4JoNt1wG1vSekXN/PjeSUJkAAJkIBLEuD47WQGcPHixViwYAFSUlLQs2dPLFq0CIMHD65T3OvXr8eSJUsQHx+P4uJi4/o5c+Zg7NixZn8ZKCCzUV3wwqWbD+EfX+5FsK+nUfOvlZ+3eZ2Wl8l04UTgyBagVWdg+jey/y/QvFheRQIkQAIk4NIEOH47kQFcu3YtJk+eDGUCBw4ciGXLlmHlypVISEhAZGTkeUJ/9NFH0a5dOwwfPhxBQUF444038K9//Qs7duxA7969zfpiUEBmYbrgRccyCzFm0WYUlVbgX7dcjpv7tje/w41/B+L+DXi2qDR/oQ1MGjH/TrySBEiABEjAyQhw/HYiA9i/f3/06dPHmNWratHR0ZgwYQLmz59vlnTVLOBtt92Gp556yqzrKSCzMNV50RnZ8HfP6z9g64EMDOjUCm9P6y8l+8ys2bf7Q2DdlMp+b3kT6CkzgWwkQAIkQAIkYCYBjt9OYgBLSkrg6+uLdevWYeJEkxmYNWuWscS7ebOUBamnVVRUoEOHDvjzn/+MP/zhD3VerZaK1U9VUwKKiIhATk4OAgIaULOuvodxgX//6JckPLo2Hl4ebkbWb8cQmckzp6XtAVaMlM2DBcCAR4Axz5oTxWtIgARIgARIoNb4HRgY6NLjdzOZiZHcy6bdkpOTER4ejm3btmHAgAHVLzNv3jysWrUK+/btq/cF1d7Bf/zjH9izZw9CQ0PrvF7tEXz66afP+zcawHrx1rrgVEGJcdxbpvz5xzFd8YcRXczroChHzN8IIPMg0FFO/LhbZgLd5bxfNhIgARIgARJoAAHOADrJDGCVAYyLi0NMTEy1BObOnYs1a9Zg7969F5XFu+++i2nTpuHjjz/GqFEXriHHGcAGfLsucumf1v2KdTtPoGuYHz57eLAxC1hvkxlarL0b2Pc5ECB7BVWx5xYNKBdT7w14AQmQAAmQgKsQoAF0EgOoswSskkfuu+8+Y/n4uuuklEgDGgXUAFhnL407lIE7V+ww/vbB/8Sgb1RL8zrZsgD45jmZ8ZMs4akbgPA+5sXxKhIgARIgARI4hwDHbycxgOpzVUkgffv2NbKAq1qPHj0wfvz4CyaBqJm/qVOnQv2pkkUa2iighhErKi3HNS9txZGMAtx9dSSemyBn9prTDmwC3r5ZrpTdCje+AvSZbE4UryEBEiABEiCBOglw/HYiA1hVBmbp0qXGMvDy5cuxYsUK7N69G1FRUZg9ezaSkpKwerUcGyZNmb577rkHL730EiZNmlQtEB8fH6iNoeY0CsgcSqZrXty4D//+5iBC/b2x6YmhCGjuWX8HWUeA5cOAomyg7xTghpfqj+EVJEACJEACJHARAhy/ncgAqs9Zzf49//zzRiHoXr16YeHChRgyRJIFpE2ZMgWJiYmIjY01/j5s2LA6s4PvvfdevPnmm2Z9cSggszAZFx1IzcO1/94qJ3+cwZK7+uCaS9vWH1xSCLw2BkjdJUu+fYH7vgQ8zCwUXX/vvIIESIAESMBFCXD8djIDaGsdU0DmEa+oOINbl32Pn46ewqjoUKy4p1/9Nf9UcvqHDwC//QfwlWSPmXLiR2C4eTfkVSRAAiRAAiRwEQIcv2kAtb4gFJB5+N7ZcQz/++EutPByx1dy3Fu7IJ/6A3csB778kyjUHbjnYyn7UveRfvV3xCtIgARIgARIoDYBjt80gFrfCQqofnxpuUUYKTX/8orK8NT1PTB1UMf6g45+D6y6HqiQ837HzgNiHqo/hleQAAmQAAmQgJkEOH7TAJoplbovo4Dqx/fQOz/j899ScFn7QHz44EC4u9Vz3FtuiiR9DAXyU4FeNwE3vSYqNfOIuPofh1eQAAmQAAmQADh+0wBqfQ0ooIvj+2ZvKqa++ZNh+j5+aCB6hdeTXV1WUjnzd1zqBIb2AKZJ+RcvM4+I0/okGUwCJEACJOBKBDh+0wBq6Z0CujC+guIyjFm4BUnZpzFjyCX432uj62f9+R+BH1cA3mIUZ3wLtOpUfwyvIAESIAESIIEGEuD4TQPYQMnUvpwCujC+Zz9LwGvfHUH7YB9sfGwIfL3qObM3/l3gI8n6Ve2OtUC3cVqfDYNJgARIgARI4EIEOH7TAGp9OyiguvHtOpGD8a9+B6n+gjfvuxLDuoVenHPKr5X1/sqKgKFPAsNna30uDCYBEiABEiCBixHg+E0DqPUNoYDOx1dWXiHmbxt2J+fixsvb4d939L4448KsyqSP7GNAFzGBavbPzU3rc2EwCZAACZAACdAAXlwDzc5Io0waR4AG8HxuK7cexnOf75Fj3jzw9RPD0FqOfbtgqyivPOP30DdAsJSHUfv+fIIb92EwigRIgARIgATMJMDxmzOAZkql7ssooNpcjmcVGokfp0vL8c+bLsVtV0ZenO/XzwBbX5Dj3aQwtMr4bdNL6/NgMAmQAAmQAAmYQ4DjNw2gOTq54DUUkAmNmki+780fEbsvHVd1bIm1M66++HFvez4D1t5V2cGklcBlt2h9FgwmARIgARIgAXMJcPymATRXK3VeRwGZsHz6azIefvcXeLm74YtZg9E51O/CbDMOyL6/4UBJHnD1g8C4+VqfA4NJgARIgARIoCEEOH7TADZEL+ddSwFVIskpLDWOe8vIL8ajo7rIT9cLcy0W07diJJCxD4gaWHnOr7un1ufAYBIgARIgARJoCAGO3zSADdELDeAFaM1evwvv/nAMnVq3MGb/vD3c675S5Ru9dw+w5xPAvy0wcwvgV0+JGK1PiMEkQAIkQAIkcD4BGkAaQK3vBQUE/HAkC7cu+97g+N7MGGP/3wXbd4uATf8nZV5kxu++L4GIK7X4M5gESIAESIAEGkOA4zcNYGN0Ux3j6gIqLivHdf/+DgfT8nHHVRGYP+myC/M8JCVe3poEnKkArnsRuPJ+LfYMJgESIAESIIHGEnD18VtxYx3AxqpH4lxdQC9tOoCFm/YjxM8bXz8+FIG+F9jLp4o8L5Niz6el6PMVdwPjXxHlNdMgz1ASIAESIAESaDwBVx+/aQAbrx0j0pUFdCg9H9cs2ooSOfnjZTnt4wY59aPOVirHu70+FkiJB9peAUzdAHhK3T82EiABEiABErATAVcev6uQcwZQQ3yuKiBV8+/25duxQ/b/DevWGm9MubLumn8q6ePjPwDxb8kJH7I3cOZmIKie4tAanwdDSYAESIAESMAcAq46ftdkQwNojlIucI2rCui9H4/jzx/8Bh9Pd2x8bAgiWvrWTein14HPHpPlXjnb9+71QCep/cdGAiRAAiRAAnYm4KrjNw2ghYTnigJStf5GvrAZOadL8ddrozF9yCV10zz+I/DGNUBFKTDqaWDQoxaizm5IgARIgARIQI+AK47f5xLjDKCGhlxRQLP+8ws+jk9Gz3YB+PihgfCQkz/Oa/lpkvQxBMhLAaJvBG5dzaQPDZ0xlARIgARIwLIEXHH8pgG0oIZcTUCx+9Iw5Y0f4SYJvB+J+busfdD5NMtlxm/1eODoNiBETgSZ/g3g7W9B6uyKBEiABEiABPQIuNr4XRctzgBqaMiVBFRYUoYxC7fgxKnTmDqwI566oUfd5Db8L7D9VcBLTJ8yf60vciycBnuGkgAJkAAJkEBjCbjS+H0hRjSAjVWPxLmSgOZ/sQfLthxGeJCPkfjRwtvjfHK73gc+OFvg+TbJ/I2+QYMuQ0mABEiABEjAOgRcafymAbSChlxFQLuTc3DjK9tQXnEGr93bDyOjw86nefJ3YOUooOw0MPgJYORTViDOLkmABEiABEhAn4CrjN8XI8UZQA0duYKAlOmbuHgbfjuRg+subYtX7+pzPrHTp4DlUuLl1BEp9TICuEtmAt3cNcgylARIgARIgASsR8AVxu/66NEA1kfoIv/uCgJ6Y9sRPP1pAvybexjHvYUGNK9NpELO9n33NuDAxsoizzOk2LOvFH1mIwESIAESIAEHJeAK43d96GkA6yPkwgYwOfs0Rr+4GQUl5Zg7sRfu6h91Po1v5wOb/wF4iDGc+l+gnRz3xkYCJEACJEACDkyABlCqs8mxXnJeF1tjCDizgJQspq/eiU17UtEvKhjvzYyBm6r/UrPtk3N91eyfahOWAFfc2RiMjCEBEiABEiABmxJw5vHbXJA0gOaSquM6ZxbQht9T8MBbP8PTvRk+f2QwuoadU8sv81Dlvr/iHODK6cB1/9IgyVASIAESIAESsB0BZx6/zaVIA2guKRcygLlFpRglx72l5RXj4RGd8cSYbrXfvqSgMuM3LQGI6A/c+5ksAXtpkGQoCZAACZAACdiOAA0gl4C11OasAvr7R79jzfaj6BjSAl/OGozmnjUyetWOAVXr7/cPAD8pB6OSPgLaanFkMAmQAAmQAAnYkoCzjt8NYcgZwIbQOudaZxTQzqOncPPSOCif9870/hjQKaT2W/9/e2cCVlW19vE/KGiY880J0azUELupfGlOaLey8mZITnVvTjnbLdNSH2349DO1K59jOZOlVA5klmampYJj3bT0elH81CQV1K6mZHAFEb71bjqEwEEOe3P2wH89j4/lOWvt9f7ed5/3v9e096q3fGxWb/vwVQdBy8hfw7Y6CLIqCZAACZAACXifgBPzt6cUHSUAFyxYgMjISJw9exYhISGYM2cOOnbsWCgT+c5LL72E/fv349ixY3jhhRe073tSnBZAGZlZ6PbWLhw9fwW9Qusjste9N+I4uTPnPb/Z14HHZgBthnmCi98lARIgARIgAUsQcFr+LglUxwjA1atXo2/fvhAR2L59eyxevBhRUVE4fPgwGjRQ59PlK4mJiZg9ezZCQ0O1vzt16lTmBeD87ccRufkoalTy1878q67+zi0pScDiMCDtAvBHtfM3YrFaQJBvV3BJIpB1SIAESIAESMDLBCgAHbQGsE2bNmjVqhUWLlTHkfxWgoOD0b17d0yfrs6qK6J07twZLVq0KNMC8OSFVDwyZwdkFHB2n3sR0bL+78Qy04F3uwJJ+4Da9wCD1KHP/gFevl15ORIgARIgARIwhgAFoEMEYEZGBgICAhATE4OIiIjc6Bg1ahQOHDiAuDi1UcEAAZieng754yoSQEFBQUhJSUGVKlWMiUoTWpEz//4a9Q32nLiIjo3/gBXPtlaDe3lG9za8COx/F6hYTW36iAVqNDKhl7wkCZAACZAACRhDgALQIQIwOTkZgYGB2L17N9q1a5cbHdOmTcPy5YDxInwAABw4SURBVMtx9OhRQwTgpEmTMHny5AJt2V0Art1/Bi/FHESF8r7YMjoMDWtW+t3G71YA659X/68Eobzjt7E6/oWFBEiABEiABGxMgALQYQJwz549aNv2912pU6dORXR0NBISEgwRgE4cAfw5NQMPzozFpbRrGP/o3RjR+c7fWSXtB5Y9BlxXo55/ehUIG2vj251dJwESIAESIIEcAhSADhGA3poCzn/jOCGAxqw5gI+/S8LddSpjw/Md1Js/fHPMTFWbPRZ3UnfJGaCpWv/X5wN19Mtvn/EXhARIgARIgARsTMAJ+VsvfsfsApZNILKjV3YBu0qzZs0QHh7OTSBuomTXsQt45p1vtM28H49oh5YNqud883om8L5aS3lyh1rvp0YEh25X6/+q6o011icBEiABEiABSxCgAHTICKBEk+sYmEWLFmnTwEuWLMHSpUsRHx+Phg0bYsKECUhKSsKKFWpN229FNohIGTx4MJo2bYqxY8fC398fIhyLU+wcQFevXdd2/f54MQ392zbE5PDmv5v85evA7rmAn1oLOGQrUCu4ODj4HRIgARIgARKwBQE752+jADtmBFCAyOjfjBkztIOgmzdvrp3vFxamzq5TZcCAAZCz/2JjY3PZ3bDT9bd/FbEo3ytOsXMARW5OwPztJ1CnSkV8OSYMlSv65Zgcvw6IGZDz373eA0J+31VdHCb8DgmQAAmQAAlYnYCd87dRbB0lAI2CUtx27BpACed+wePzdiEzKxuL+4bikZA6OSb/pDbLLP0TcC0VaPcC0GVKcVHweyRAAiRAAiRgGwJ2zd9GAqYA1EHTjgGUpURfD/Wu3+9PXVbCr7YSgP+VQ+BqSo74u3gcaKRGTZ9RI4Hl1Pt+WUiABEiABEjAYQTsmL+NdgEFoA6idgyg6L2JeO3TeNxaoTy+Uq97q1O1IpCVpRZRPgMc3QhUUW8AGaYOzq70Bx1kWJUESIAESIAErEvAjvnbaJoUgDqI2i2AzqVcxcOz4nAlPRP/Ex6Cfm1vz7F+RySw7Q014qfe/fvsF0BgqA4qrEoCJEACJEAC1iZgt/xdGjQpAHVQtVsADY/ejy/iz6FFUDWsVce+lPNV578c/0od+dJTUcgGnngLaNVPBxFWJQESIAESIAHrE7Bb/i4NohSAOqjaKYC2KOE3VAnA8kr0yYHPwXXVu4svJeYc9nz1MhA6AOimjn5hIQESIAESIAGHE7BT/i4tV1AA6iBrlwD6VU35ytTvWTUFLK96k1e+ISNNveatC3DuUM6U78BNQPkKOmiwKgmQAAmQAAnYg4Bd8ndp0qQA1EHXLgE0aX083tuTiAY1ArBldBgqllevdFs3HPjnKiBAbfYYpt74UTVQBwlWJQESIAESIAH7ELBL/i5NohSAOujaIYAOnL6MiAW7ka2W+EUPao2OjW8DvlkCbBqr3gNTDuj3qTr2paMOCqxKAiRAAiRAAvYiYIf8XdpEKQB1ELZ6AF27noUn3t6NI2d/QUTLQMzu0wL4cS+w/HF19It632+XqerA57/pIMCqJEACJEACJGA/AlbP394gSgGog7LVA2hR3Am8uSkB1QL8sFWd+Vcz+5La9KEOef71vHrF25NAz2VqFFDtBGYhARIgARIggTJEwOr52xuuoADUQdnKAXTqYhq6zInD1WtZiOz5R/RqUVuN/HUDTn8N1GoGDFbHv/hX0mE9q5IACZAACZCAPQlYOX97iygFoA7SVg2gbLXgr9+yf2DnsQtoe0dNfDikDXw2jQP+odb+VagKDN0O1LxTh+WsSgIkQAIkQAL2JWDV/O1NohSAOmhbNYA+PZCEUasOwF/t9t38YhganVkPfKJ2/Up5ejXQ9FEdVrMqCZAACZAACdibgFXztzepUgDqoG3FALqcloEHZ8bhYmoGXu7SBH8LVuf9vaPO+8u8CnQaDzwwUYfFrEoCJEACJEAC9idgxfztbaoUgDqIWzGAxn10EGv2nUGT2rfis0Eh8F/2AHD5FNBYiUAZ/fNVZwCykAAJkAAJkEAZJmDF/O1td1AA6iButQDae+Iinl6qNnmosnZYa4TuGgac2ApUv12t+4sFbqmuw1pWJQESIAESIAFnELBa/jaDKgWgDupWCqCr166j69yd+OFCKv7apgGmVvkE2Pm/6vVut+Ts+K3TXIelrEoCJEACJEACziFgpfxtFlUKQB3krRRAs778P8zbegy1KlfA9sdTUWldvxzLnowC/thLh5WsSgIkQAIkQALOImCl/G0WWQpAHeStEkDHzl9B13k7ce16NpY/UQ2dYvsAGVeA+0cCj07XYSGrkgAJkAAJkIDzCFglf5tJlgJQB30rBFBWVjb6LNmLbxMv4c9Nb8XbqWPhc+Eo0LB9znt+y/npsJBVSYAESIAESMB5BKyQv82mSgGowwNWCKAPvzmFiesOIcDfF/ubROOW4xuBynXVpo849bd6+wcLCZAACZAACZDADQSskL/NdgkFoA4PmB1AP/1yFQ/OisOVq5lY0/wfaH18jjrmRY34DfwcCGqtwzJWJQESIAESIAHnEjA7f1uBLAWgDi+YHUDPffgdNv7zLJ6pdRJTrrwGn+ws4M+zgPsG6bCKVUmABEiABEjA2QTMzt9WoEsBqMMLZgbQtoTzePa9fQjyvYDtlf8b5dMvAS2eAcLfBnx8dFjFqiRAAiRAAiTgbAJm5m+rkKUA1OEJswIoNT0TXWbvwIXLKYit8SbqpiUAdVsAz34B+Klz/1hIgARIgARIgATcEjArf1vJJRSAOrxhVgC98dlhRO36AW9XegePX9+m3vBRAximNn1Ua6DDGlYlARIgARIggbJBwKz8bSW6FIA6vGFGAB06k4Lw+bvwlO9WTPN7R033qnf7PrMWuPNPOixhVRIgARIgARIoOwTMyN9Wo0sBqMMj3g6gzOtZ6L5gN/yS9yOmwhSURybw0CSgw2gdVrAqCZAACZAACZQtAt7O31akSwGowyveDqConT9g0cavsbHiK6iNn4HgJ4DeK7jpQ4cPWZUESIAESKDsEfB2/rYiYQpAHV7xZgCduZSGx2Ztw1KfN3C/7xHgD02AIWr9X4XKOixgVRIgARIgARIoewS8mb+tSpcCUIdnvBVA2dnZGLR8H9odn4nB5Tch278yfET83aZEIAsJkAAJkAAJkIBHBLyVvz3qlJe/TAGoA7i3AuizfyZjy6r5mOevzviT0ud9Nf3bTUfPWZUESIAESIAEyi4Bb+VvKxOmANThHW8EUEraNQyduQLvZk5AgE+62vAxRm38+G8dvWZVEiABEiABEijbBLyRv61OmAJQh4e8EUCTY/ag/6EBuN33PLIadYZv34/V+37L6eg1q5IACZAACZBA2SbgjfxtdcKOEoALFixAZGQkzp49i5CQEMyZMwcdO3Z064O4uDiMGTMG8fHxqFevHsaNG4fhw4cX22elHUDfnryAX5b1xIPlvkd6pUBUeG4XEKAOfWYhARIgARIgARIoMYHSzt8l7pgXKzpGAK5evRp9+/aFiMD27dtj8eLFiIqKwuHDh9GgQcE3ZJw8eRLNmzfHkCFDMGzYMOzevRsjR47EypUr0aNHj2K5oDQDKD3zOlZFPof+6StxzccffkO+BOqp172xkAAJkAAJkAAJ6CJQmvlbV8e8WNkxArBNmzZo1aoVFi5cmIsvODgY3bt3x/Tp0wsgHT9+PNavX48jR9SRKr8VGf07ePAg9u7dWywXlGYArV+zDE8czjngOa3r2who3bdYfeKXSIAESIAESIAEiiZQmvnbLuwdIQAzMjIQEBCAmJgYRERE5LIfNWoUDhw4AJnqzV/CwsLQsmVLzJ07N/ejdevWoXfv3khLS4Ofn1+BOunp6ZA/riIBFBQUhJSUFFSpUsUwn/947BCqv98FVXzScLLR02jUf5FhbbMhEiABEiABEijrBCgA1Tsk1Blz2XYPhOTkZAQGBmrTuO3atcs1Z9q0aVi+fDmOHj1awMQmTZpgwIABmDhxYu5ne/bs0aaPpb26desWqDNp0iRMnjy5wL8bLQC/j3wcLVN34liFZrhrbCx8ylewu4vYfxIgARIgARKwDAEKQIcJQBFwbdu2zQ2wqVOnIjo6GgkJCYUKwIEDB2LChAm5n4mA7NChg7aJpE6dOqaNAKZcuohj0S+gbvcpCGxwh2VuGHaEBEiABEiABJxAgALQIQLQW1PA+YOeAeSEnwHaQAIkQAIkUNYIMH87RABK4MomkNDQUG0XsKs0a9YM4eHhbjeBbNiwQdsl7CojRozQ1gxaYRNIWbsZaS8JkAAJkAAJeIsABaCDBKDrGJhFixZp08BLlizB0qVLtTP+GjZsqE31JiUlYcWKFVp8uY6BkSNg5CgYEX2yC9gqx8B46ybgdUiABEiABEigrBGgAHSQAJTgldG/GTNmaGv45Iy/2bNnQ3b7SpENH4mJiYiNjc2Nc9kdPHr06NyDoOVoGCsdBF3WbkjaSwIkQAIkQALeIEAB6DAB6I2gyXsNBpC3ifN6JEACJEACJKCfAPM3BaCuKGIA6cLHyiRAAiRAAiRgCgHmbwpAXYHHANKFj5VJgARIgARIwBQCzN8UgLoCjwGkCx8rkwAJkAAJkIApBJi/KQB1BR4DSBc+ViYBEiABEiABUwgwf1MA6go8BpAufKxMAiRAAiRAAqYQYP6mANQVeAwgXfhYmQRIgARIgARMIcD8TQGoK/AYQLrwsTIJkAAJkAAJmEKA+ZsCUFfgMYB04WNlEiABEiABEjCFAPM3BaCuwGMA6cLHyiRAAiRAAiRgCgHmbwpAXYGXkpKCatWq4fTp06hSpYqutliZBEiABEiABEjAOwREAAYFBeHy5cuoWrWqdy5qsav4ZKtisT7ZpjtnzpzRAoiFBEiABEiABEjAfgRkAKd+/fr267gBPaYA1AExKysLycnJqFy5Mnx8fHS0VLCq6+nEqaOLtM/QcDGlMfrQFOyGXdTp/hNQTreR9pX8dpCxrytXrqBevXrw9fUteUM2rkkBaFHnOX19Au2zaOB50C360ANYFvyq0/3nEoAyvSfLdZy4TMfpPnS6fWb/LFAAmu0BN9d3euDTPosGngfdog89gGXBrzrdfxSAFgw6D7tUFmLUQySGfp0C0FCcxjXm9MCnfcbFilkt0YdmkTfmuk73HwWgMXFiZitlIUbN5EsBaCb9Iq6dnp6O6dOnY8KECahQoYJFe1nybtG+krOzSk360CqeKFk/nO4/oeJ0G2lfyWKftXIIUAAyEkiABEiABEiABEigjBGgACxjDqe5JEACJEACJEACJEAByBggARIgARIgARIggTJGgAKwjDmc5pIACZAACZAACZAABSBjgARIgARIgARIgATKGAEKQJMcvmDBAkRGRuLs2bMICQnBnDlz0LFjR7e9iYuLw5gxYxAfH6+dXD5u3DgMHz7cpN4XfVnZvfzxxx8jISEBt9xyC9q1a4e///3vaNq0qduKsbGxeOCBBwp8fuTIEdx9992WsnPSpEmYPHnyDX2qXbs2zp075wj/iRG33347fvzxxwL2jBw5EvPnzy/w71b3344dO7T7bf/+/do9t27dOnTv3j3XDnkrgPh0yZIluHTpEtq0aaPZKfdmUWXt2rV47bXXcOLECdx5552YOnUqIiIivB6vRdl37do1vPrqq/j888/xww8/aO89feihh/Dmm29qvyXuynvvvYeBAwcW+Pg///kPKlasaCkbpTMDBgzA8uXLb+iX+PHrr7+2vQ/FAHdvm5oxYwbGjh1bqI1W8mFx8oLd70Ov3xQ6L0gBqBNgSaqvXr0affv2hYjA9u3bY/HixYiKisLhw4fRoEGDAk2ePHkSzZs3x5AhQzBs2DDs3r0bkohXrlyJHj16lKQLpVrn0UcfxVNPPYX77rsPmZmZeOWVV3Do0CHNvkqVKhV6bZeAOHr06A0n9t92220oV65cqfbX08ZFAH700Uf46quvcqtKH6WvhRW7+U9s+Pe//43r16/nmvOvf/0LDz/8MLZv347OnTsXMNPq/tu0aZN237Rq1Uq7Z/ILQHlAEfEmCbNJkyZ44403IKJK4lFe9VhY2bt3r/bQNmXKFE30SZuvv/46du3apQlIb5ai7JO3YPTs2VP7/bj33ns1gfviiy9q9+a+ffvcdlNYjBo1SmOQt9SpU8ebpuVe62Y+FAF4/vx5vPvuu7l1/P39UaNGDbf9tYsPxYD8D5jCY9CgQTh+/DjuuOMOtwLQKj4sTl6w+31oyo2h46IUgDrglbSqJAdJRAsXLsxtIjg4WBuRkKek/GX8+PFYv349ZDTMVWT07+DBg5AfMKsXERO1atWCjGKGhYUV2l2XgJDkVK1aNUubJALwk08+wYEDB4rVT7v7T4wUwfDZZ5/h2LFjhY5E2Ml/MpKSVwDKqIOMhImN4ispcr6ajOpKQpKHrsJKnz59tHfNSiJ2FUly1atX1x7OzCr57SusH99++y1at26tjfIW9tApdUQACpPLly+bZYrb6xZmowhA6avcm8Utdvah5At5l+3WrVvdmmtlH+bPC067D4sbg2Z+jwLQy/QzMjIQEBCAmJiYG6aK5ClNBIWIpPxFRFPLli0xd+7c3I8kgfXu3RtpaWnw8/PzshWeXU6eUBs3bqyNAspIZmHFJSBk6vHq1ato1qyZNm1V2LSwZ1c3/tsiAGU6UabS5JBuEfTTpk1z+xRud/9JzIpAkiUIEydOtL3/8osHmRaV6dvvvvtOu89cJTw8XHsYyT+t6PpchNPo0aO1P64ye/ZsbTlHYdPnxkdi4S0WRwDK6HWXLl00weTuHbkiHgYPHozAwEBtNLhFixbaaGdeRt6yKf913AlAEX8y6id+69SpkzaqKw+f7opdfSgjnfXr19di8y9/+UuRAtCqPsyfF5x2H5p1b3hyXQpAT2gZ8N3k5GTtB1Wmo2RtnKuIgJCbOf90i3wuU1LydJs3+e7Zs0ebPpb26tata0DPSqcJeaqTRCojezt37nR7EbFbptxCQ0O10Zfo6GgsWrQIIgzdjRqWTo9v3qqM+IjwFr/ID7FMF8p6R1mfWbNmzQIN2Nl/YsyaNWu0JHPq1Cm3a8bs5L/84sF1LyUlJd1g39ChQzUht3nz5kKDQoSGiKS8CfjDDz/U1s1JDJtVbiYA5QGrQ4cO2tra999/3203Ze2cJOl77rlHG+mUB1BZRygzD/JAZ2YpzEZZWnPrrbeiYcOGkGUXsjZTprll3ae7tynZ1Yey7k/WcMrvf1HrMa3qw8LygtPuQzPvj+JemwKwuKQM+p5LAEqwt23bNrdVeVIV0SNCIn8RASFJRV4L5yoiIOVHXBa0m7UmpzhInnvuOWzcuFFbFyVPrJ6Ubt26adONMv1t5ZKamqqNIMnGHBklc5L/xJZHHnlEG1XZsGGDR26wqv/cCcD8D1OyZu706dP44osv3ApAeWh7+umncz//4IMPtHVZIrLMKkUJQNkQ0qtXL03My8OVu9G/wvqelZWlLV2RB7J58+aZZZ523ZuJXPmO/DaKGFy1ahWefPJJx/hQDBHxLmty33rrLY/8YBUfFpYX3A1q2PU+9MgxJn2ZAtDL4MvSFPDzzz+vrceRkb1GjRp5TFpEsYxQ5F376HEjXqogP8Z33XXXDes6XZe28xSwjIDJAnPZ1S0juZ4Uq/qvrE4Bi/iTZSMy1bZt27ZCR6tv5l9JxmfOnLlh3ePN6pTG58URgHJdGamUKVDX2s78fbHjFLDMpMhviiwZkk09nhazfeguL3AK2FNP6v8+BaB+hh63IGvGZKpTdgG7iqx5kwTrbhOIjL7ILlpXGTFihPYDYMVNIDK8Lze5rFOUUYaSThfJzsWff/5ZS1ZWLjLdJyOAMmUou0DzF0k+dvJf3v7LekfZpS4jYeXLl/fIDVb1n7tNILKWT0ZxpciDmqwdu9kmEFmEL9OirvLYY49p68+stgnEJf5kE4/s5Ha3Y70oB8t9LRtHZEp42bJlHsWC0V8ujgC8ePGittxGjvbp169foV2QTSB28aHLAFkOJLvyi9rB7Y63mT68WV5wbQJxyn1odMyXRnsUgKVB9SZtuo6BkTVuMg0sP1BLly7V1pDJlIVM9cp6pBUrVmgtuY4Rkd2I8vQmok92AVv1GBg5okbWQn366ac3nP0nmybkXEAp+W2UhfOyAUTOXZPkKyN/ssZFzllzN31jguu0S7788suQ6U0ZPfjpp5+0NYCyeUc2uTjBfy6uMl0kI7cyxSm+yFvs5r9ff/1VW88mRTYxzJo1S9tgJEeEiB9F6MnDlxwhIg8ssiZXHl7yHgMjIkIEheshTaasZCRGRjrl4U3iXTYumXEMTFH2yQYeOfpGNrnITm7Z3ewqYr9M70vJb5+ci3j//fdrPGQNoEz7yjIVWX4iQtDbpSgbxQ55WBE7ZU10YmKitmZaprplBsF1lI9dfejaqS1+EPtmzpxZ6DmwVvZhcfKC3e9Db98Teq9HAaiXYAnry+ifLOSVdSqyM1Z2D7o2O8gTnvyASQJyFREY8mTkOghaRpWsehC0uwNLJbmKbVLy2ygsRAiL8BWRKEJQREbXrl1LSLj0qskZhzKtfeHCBW0kRZKk7I6UUdzCbJN/s5P/XOS2bNmirf8TESTrUPMWu/nP3UHV/fv31zZyuA6gldHOvAdB5921LucfykOKfN9V5DxIEX2u6SsRg2Y8sBRlnwgjd0sw8p7rmN8++b2RqX85f04e3kQ4S1t51y6X3l1WsOWibJQjteRYlO+//17b2SwiSQS+3JdBQUG5jdnVh66Yk99IOZpH8ob4JH+xsg+Lkxfsfh96834w4loUgEZQZBskQAIkQAIkQAIkYCMCFIA2cha7SgIkQAIkQAIkQAJGEKAANIIi2yABEiABEiABEiABGxGgALSRs9hVEiABEiABEiABEjCCAAWgERTZBgmQAAmQAAmQAAnYiAAFoI2cxa6SAAmQAAmQAAmQgBEEKACNoMg2SIAESIAESIAESMBGBCgAbeQsdpUESIAESIAESIAEjCBAAWgERbZBAiRAAiRAAiRAAjYiQAFoI2exqyRAAiRAAiRAAiRgBAEKQCMosg0SIAESIAESIAESsBEBCkAbOYtdJQESIAESIAESIAEjCFAAGkGRbZAACZAACZAACZCAjQhQANrIWewqCZAACZAACZAACRhBgALQCIpsgwRIgARIgARIgARsRIAC0EbOYldJgARIgARIgARIwAgCFIBGUGQbJEACJEACJEACJGAjAhSANnIWu0oCJEACJEACJEACRhCgADSCItsgARIgARIgARIgARsRoAC0kbPYVRIgARIgARIgARIwggAFoBEU2QYJkAAJkAAJkAAJ2IgABaCNnMWukgAJkAAJkAAJkIARBCgAjaDINkiABEiABEiABEjARgQoAG3kLHaVBEiABEiABEiABIwgQAFoBEW2QQIkQAIkQAIkQAI2IkABaCNnsaskQAIkQAIkQAIkYAQBCkAjKLINEiABEiABEiABErARAQpAGzmLXSUBEiABEiABEiABIwhQABpBkW2QAAmQAAmQAAmQgI0IUADayFnsKgmQAAmQAAmQAAkYQYAC0AiKbIMESIAESIAESIAEbESAAtBGzmJXSYAESIAESIAESMAIAv8PMpMw5fWM7qIAAAAASUVORK5CYII=\" width=\"640\">"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/plain": [
       "[<matplotlib.lines.Line2D at 0x7f70469f8220>]"
      ]
     },
     "execution_count": 52,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "exposure_tot = 10 #ms\n",
    "exp_kp1 = 2 #ms\n",
    "\n",
    "Vhigh = 127\n",
    "Vlow = 64\n",
    "\n",
    "Vtfl2 = 100\n",
    "\n",
    "\n",
    "def model(x, rate):\n",
    "    kp1_ms = exposure_tot - exp_kp1\n",
    "    before_kneepoint = np.clip(Vhigh - x * rate, Vtfl2, Vhigh)\n",
    "    at_kp = np.clip(Vhigh - kp1_ms * rate, Vtfl2, Vhigh)\n",
    "    after_kneepoint = np.clip(at_kp - (x - kp1_ms) * rate, Vlow, Vhigh)\n",
    "    return np.where(x < kp1_ms, before_kneepoint, after_kneepoint)\n",
    "\n",
    "def at_end(rate, tot=exposure_tot, kp1=exp_kp1):\n",
    "    kp1_ms = tot - kp1\n",
    "    at_kp = np.clip(Vhigh - kp1_ms * rate, Vtfl2, Vhigh)\n",
    "    at_end = np.clip(at_kp - kp1 * rate, Vlow, Vhigh)\n",
    "    return at_end\n",
    "\n",
    "def norm(v):\n",
    "    return 1.0 - (v - Vlow) / (Vhigh - Vlow)\n",
    "    \n",
    "plt.figure()\n",
    "total_exposure = np.linspace(0, 20, 1000)\n",
    "kp1 = total_exposure * 0.2\n",
    "plt.plot(total_exposure, norm(at_end(5, total_exposure, kp1)))\n",
    "plt.plot(total_exposure, norm(at_end(4, total_exposure, kp1)))"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "language": "/gnu/store/kg33ii73xhnl64ld34dv2afv66ayn7g8-python-3.9.9/bin/python3",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.9.9"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}