Supervised classification.ipynb

       "    );\n",
       "\n",
       "    fig.rubberband_context.strokeRect(min_x, min_y, width, height);\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_figure_label = function (fig, msg) {\n",
       "    // Updates the figure title.\n",
       "    fig.header.textContent = msg['label'];\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_cursor = function (fig, msg) {\n",
       "    var cursor = msg['cursor'];\n",
       "    switch (cursor) {\n",
       "        case 0:\n",
       "            cursor = 'pointer';\n",
       "            break;\n",
       "        case 1:\n",
       "            cursor = 'default';\n",
       "            break;\n",
       "        case 2:\n",
       "            cursor = 'crosshair';\n",
       "            break;\n",
       "        case 3:\n",
       "            cursor = 'move';\n",
       "            break;\n",
       "    }\n",
       "    fig.rubberband_canvas.style.cursor = cursor;\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_message = function (fig, msg) {\n",
       "    fig.message.textContent = msg['message'];\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_draw = function (fig, _msg) {\n",
       "    // Request the server to send over a new figure.\n",
       "    fig.send_draw_message();\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_image_mode = function (fig, msg) {\n",
       "    fig.image_mode = msg['mode'];\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_history_buttons = function (fig, msg) {\n",
       "    for (var key in msg) {\n",
       "        if (!(key in fig.buttons)) {\n",
       "            continue;\n",
       "        }\n",
       "        fig.buttons[key].disabled = !msg[key];\n",
       "        fig.buttons[key].setAttribute('aria-disabled', !msg[key]);\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_navigate_mode = function (fig, msg) {\n",
       "    if (msg['mode'] === 'PAN') {\n",
       "        fig.buttons['Pan'].classList.add('active');\n",
       "        fig.buttons['Zoom'].classList.remove('active');\n",
       "    } else if (msg['mode'] === 'ZOOM') {\n",
       "        fig.buttons['Pan'].classList.remove('active');\n",
       "        fig.buttons['Zoom'].classList.add('active');\n",
       "    } else {\n",
       "        fig.buttons['Pan'].classList.remove('active');\n",
       "        fig.buttons['Zoom'].classList.remove('active');\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.updated_canvas_event = function () {\n",
       "    // Called whenever the canvas gets updated.\n",
       "    this.send_message('ack', {});\n",
       "};\n",
       "\n",
       "// A function to construct a web socket function for onmessage handling.\n",
       "// Called in the figure constructor.\n",
       "mpl.figure.prototype._make_on_message_function = function (fig) {\n",
       "    return function socket_on_message(evt) {\n",
       "        if (evt.data instanceof Blob) {\n",
       "            /* FIXME: We get \"Resource interpreted as Image but\n",
       "             * transferred with MIME type text/plain:\" errors on\n",
       "             * Chrome.  But how to set the MIME type?  It doesn't seem\n",
       "             * to be part of the websocket stream */\n",
       "            evt.data.type = 'image/png';\n",
       "\n",
       "            /* Free the memory for the previous frames */\n",
       "            if (fig.imageObj.src) {\n",
       "                (window.URL || window.webkitURL).revokeObjectURL(\n",
       "                    fig.imageObj.src\n",
       "                );\n",
       "            }\n",
       "\n",
       "            fig.imageObj.src = (window.URL || window.webkitURL).createObjectURL(\n",
       "                evt.data\n",
       "            );\n",
       "            fig.updated_canvas_event();\n",
       "            fig.waiting = false;\n",
       "            return;\n",
       "        } else if (\n",
       "            typeof evt.data === 'string' &&\n",
       "            evt.data.slice(0, 21) === 'data:image/png;base64'\n",
       "        ) {\n",
       "            fig.imageObj.src = evt.data;\n",
       "            fig.updated_canvas_event();\n",
       "            fig.waiting = false;\n",
       "            return;\n",
       "        }\n",
       "\n",
       "        var msg = JSON.parse(evt.data);\n",
       "        var msg_type = msg['type'];\n",
       "\n",
       "        // Call the  \"handle_{type}\" callback, which takes\n",
       "        // the figure and JSON message as its only arguments.\n",
       "        try {\n",
       "            var callback = fig['handle_' + msg_type];\n",
       "        } catch (e) {\n",
       "            console.log(\n",
       "                \"No handler for the '\" + msg_type + \"' message type: \",\n",
       "                msg\n",
       "            );\n",
       "            return;\n",
       "        }\n",
       "\n",
       "        if (callback) {\n",
       "            try {\n",
       "                // console.log(\"Handling '\" + msg_type + \"' message: \", msg);\n",
       "                callback(fig, msg);\n",
       "            } catch (e) {\n",
       "                console.log(\n",
       "                    \"Exception inside the 'handler_\" + msg_type + \"' callback:\",\n",
       "                    e,\n",
       "                    e.stack,\n",
       "                    msg\n",
       "                );\n",
       "            }\n",
       "        }\n",
       "    };\n",
       "};\n",
       "\n",
       "// from http://stackoverflow.com/questions/1114465/getting-mouse-location-in-canvas\n",
       "mpl.findpos = function (e) {\n",
       "    //this section is from http://www.quirksmode.org/js/events_properties.html\n",
       "    var targ;\n",
       "    if (!e) {\n",
       "        e = window.event;\n",
       "    }\n",
       "    if (e.target) {\n",
       "        targ = e.target;\n",
       "    } else if (e.srcElement) {\n",
       "        targ = e.srcElement;\n",
       "    }\n",
       "    if (targ.nodeType === 3) {\n",
       "        // defeat Safari bug\n",
       "        targ = targ.parentNode;\n",
       "    }\n",
       "\n",
       "    // pageX,Y are the mouse positions relative to the document\n",
       "    var boundingRect = targ.getBoundingClientRect();\n",
       "    var x = e.pageX - (boundingRect.left + document.body.scrollLeft);\n",
       "    var y = e.pageY - (boundingRect.top + document.body.scrollTop);\n",
       "\n",
       "    return { x: x, y: y };\n",
       "};\n",
       "\n",
       "/*\n",
       " * return a copy of an object with only non-object keys\n",
       " * we need this to avoid circular references\n",
       " * http://stackoverflow.com/a/24161582/3208463\n",
       " */\n",
       "function simpleKeys(original) {\n",
       "    return Object.keys(original).reduce(function (obj, key) {\n",
       "        if (typeof original[key] !== 'object') {\n",
       "            obj[key] = original[key];\n",
       "        }\n",
       "        return obj;\n",
       "    }, {});\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.mouse_event = function (event, name) {\n",
       "    var canvas_pos = mpl.findpos(event);\n",
       "\n",
       "    if (name === 'button_press') {\n",
       "        this.canvas.focus();\n",
       "        this.canvas_div.focus();\n",
       "    }\n",
       "\n",
       "    var x = canvas_pos.x * this.ratio;\n",
       "    var y = canvas_pos.y * this.ratio;\n",
       "\n",
       "    this.send_message(name, {\n",
       "        x: x,\n",
       "        y: y,\n",
       "        button: event.button,\n",
       "        step: event.step,\n",
       "        guiEvent: simpleKeys(event),\n",
       "    });\n",
       "\n",
       "    /* This prevents the web browser from automatically changing to\n",
       "     * the text insertion cursor when the button is pressed.  We want\n",
       "     * to control all of the cursor setting manually through the\n",
       "     * 'cursor' event from matplotlib */\n",
       "    event.preventDefault();\n",
       "    return false;\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._key_event_extra = function (_event, _name) {\n",
       "    // Handle any extra behaviour associated with a key event\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.key_event = function (event, name) {\n",
       "    // Prevent repeat events\n",
       "    if (name === 'key_press') {\n",
       "        if (event.which === this._key) {\n",
       "            return;\n",
       "        } else {\n",
       "            this._key = event.which;\n",
       "        }\n",
       "    }\n",
       "    if (name === 'key_release') {\n",
       "        this._key = null;\n",
       "    }\n",
       "\n",
       "    var value = '';\n",
       "    if (event.ctrlKey && event.which !== 17) {\n",
       "        value += 'ctrl+';\n",
       "    }\n",
       "    if (event.altKey && event.which !== 18) {\n",
       "        value += 'alt+';\n",
       "    }\n",
       "    if (event.shiftKey && event.which !== 16) {\n",
       "        value += 'shift+';\n",
       "    }\n",
       "\n",
       "    value += 'k';\n",
       "    value += event.which.toString();\n",
       "\n",
       "    this._key_event_extra(event, name);\n",
       "\n",
       "    this.send_message(name, { key: value, guiEvent: simpleKeys(event) });\n",
       "    return false;\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.toolbar_button_onclick = function (name) {\n",
       "    if (name === 'download') {\n",
       "        this.handle_save(this, null);\n",
       "    } else {\n",
       "        this.send_message('toolbar_button', { name: name });\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.toolbar_button_onmouseover = function (tooltip) {\n",
       "    this.message.textContent = tooltip;\n",
       "};\n",
       "\n",
       "///////////////// REMAINING CONTENT GENERATED BY embed_js.py /////////////////\n",
       "// prettier-ignore\n",
       "var _JSXTOOLS_RESIZE_OBSERVER=function(A){var t,i=new WeakMap,n=new WeakMap,a=new WeakMap,r=new WeakMap,o=new Set;function s(e){if(!(this instanceof s))throw new TypeError(\"Constructor requires 'new' operator\");i.set(this,e)}function h(){throw new TypeError(\"Function is not a constructor\")}function c(e,t,i,n){e=0 in arguments?Number(arguments[0]):0,t=1 in arguments?Number(arguments[1]):0,i=2 in arguments?Number(arguments[2]):0,n=3 in arguments?Number(arguments[3]):0,this.right=(this.x=this.left=e)+(this.width=i),this.bottom=(this.y=this.top=t)+(this.height=n),Object.freeze(this)}function d(){t=requestAnimationFrame(d);var s=new WeakMap,p=new Set;o.forEach((function(t){r.get(t).forEach((function(i){var r=t instanceof window.SVGElement,o=a.get(t),d=r?0:parseFloat(o.paddingTop),f=r?0:parseFloat(o.paddingRight),l=r?0:parseFloat(o.paddingBottom),u=r?0:parseFloat(o.paddingLeft),g=r?0:parseFloat(o.borderTopWidth),m=r?0:parseFloat(o.borderRightWidth),w=r?0:parseFloat(o.borderBottomWidth),b=u+f,F=d+l,v=(r?0:parseFloat(o.borderLeftWidth))+m,W=g+w,y=r?0:t.offsetHeight-W-t.clientHeight,E=r?0:t.offsetWidth-v-t.clientWidth,R=b+v,z=F+W,M=r?t.width:parseFloat(o.width)-R-E,O=r?t.height:parseFloat(o.height)-z-y;if(n.has(t)){var k=n.get(t);if(k[0]===M&&k[1]===O)return}n.set(t,[M,O]);var S=Object.create(h.prototype);S.target=t,S.contentRect=new c(u,d,M,O),s.has(i)||(s.set(i,[]),p.add(i)),s.get(i).push(S)}))})),p.forEach((function(e){i.get(e).call(e,s.get(e),e)}))}return s.prototype.observe=function(i){if(i instanceof window.Element){r.has(i)||(r.set(i,new Set),o.add(i),a.set(i,window.getComputedStyle(i)));var n=r.get(i);n.has(this)||n.add(this),cancelAnimationFrame(t),t=requestAnimationFrame(d)}},s.prototype.unobserve=function(i){if(i instanceof window.Element&&r.has(i)){var n=r.get(i);n.has(this)&&(n.delete(this),n.size||(r.delete(i),o.delete(i))),n.size||r.delete(i),o.size||cancelAnimationFrame(t)}},A.DOMRectReadOnly=c,A.ResizeObserver=s,A.ResizeObserverEntry=h,A}; // eslint-disable-line\n",
       "mpl.toolbar_items = [[\"Home\", \"Reset original view\", \"fa fa-home icon-home\", \"home\"], [\"Back\", \"Back to previous view\", \"fa fa-arrow-left icon-arrow-left\", \"back\"], [\"Forward\", \"Forward to next view\", \"fa fa-arrow-right icon-arrow-right\", \"forward\"], [\"\", \"\", \"\", \"\"], [\"Pan\", \"Left button pans, Right button zooms\\nx/y fixes axis, CTRL fixes aspect\", \"fa fa-arrows icon-move\", \"pan\"], [\"Zoom\", \"Zoom to rectangle\\nx/y fixes axis, CTRL fixes aspect\", \"fa fa-square-o icon-check-empty\", \"zoom\"], [\"\", \"\", \"\", \"\"], [\"Download\", \"Download plot\", \"fa fa-floppy-o icon-save\", \"download\"]];\n",
       "\n",
       "mpl.extensions = [\"eps\", \"jpeg\", \"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.close = function () {\n",
       "        comm.close();\n",
       "    };\n",
       "    ws.send = function (m) {\n",
       "        //console.log('sending', m);\n",
       "        comm.send(m);\n",
       "    };\n",
       "    // Register the callback with on_msg.\n",
       "    comm.on_msg(function (msg) {\n",
       "        //console.log('receiving', msg['content']['data'], msg);\n",
       "        // Pass the mpl event to the overridden (by mpl) onmessage function.\n",
       "        ws.onmessage(msg['content']['data']);\n",
       "    });\n",
       "    return ws;\n",
       "};\n",
       "\n",
       "mpl.mpl_figure_comm = function (comm, msg) {\n",
       "    // This is the function which gets called when the mpl process\n",
       "    // starts-up an IPython Comm through the \"matplotlib\" channel.\n",
       "\n",
       "    var id = msg.content.data.id;\n",
       "    // Get hold of the div created by the display call when the Comm\n",
       "    // socket was opened in Python.\n",
       "    var element = document.getElementById(id);\n",
       "    var ws_proxy = comm_websocket_adapter(comm);\n",
       "\n",
       "    function ondownload(figure, _format) {\n",
       "        window.open(figure.canvas.toDataURL());\n",
       "    }\n",
       "\n",
       "    var fig = new mpl.figure(id, ws_proxy, ondownload, element);\n",
       "\n",
       "    // Call onopen now - mpl needs it, as it is assuming we've passed it a real\n",
       "    // web socket which is closed, not our websocket->open comm proxy.\n",
       "    ws_proxy.onopen();\n",
       "\n",
       "    fig.parent_element = element;\n",
       "    fig.cell_info = mpl.find_output_cell(\"<div id='\" + id + \"'></div>\");\n",
       "    if (!fig.cell_info) {\n",
       "        console.error('Failed to find cell for figure', id, fig);\n",
       "        return;\n",
       "    }\n",
       "    fig.cell_info[0].output_area.element.on(\n",
       "        'cleared',\n",
       "        { fig: fig },\n",
       "        fig._remove_fig_handler\n",
       "    );\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_close = function (fig, msg) {\n",
       "    var width = fig.canvas.width / fig.ratio;\n",
       "    fig.cell_info[0].output_area.element.off(\n",
       "        'cleared',\n",
       "        fig._remove_fig_handler\n",
       "    );\n",
       "    fig.resizeObserverInstance.unobserve(fig.canvas_div);\n",
       "\n",
       "    // Update the output cell to use the data from the current canvas.\n",
       "    fig.push_to_output();\n",
       "    var dataURL = fig.canvas.toDataURL();\n",
       "    // Re-enable the keyboard manager in IPython - without this line, in FF,\n",
       "    // the notebook keyboard shortcuts fail.\n",
       "    IPython.keyboard_manager.enable();\n",
       "    fig.parent_element.innerHTML =\n",
       "        '<img src=\"' + dataURL + '\" width=\"' + width + '\">';\n",
       "    fig.close_ws(fig, msg);\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.close_ws = function (fig, msg) {\n",
       "    fig.send_message('closing', msg);\n",
       "    // fig.ws.close()\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.push_to_output = function (_remove_interactive) {\n",
       "    // Turn the data on the canvas into data in the output cell.\n",
       "    var width = this.canvas.width / this.ratio;\n",
       "    var dataURL = this.canvas.toDataURL();\n",
       "    this.cell_info[1]['text/html'] =\n",
       "        '<img src=\"' + dataURL + '\" width=\"' + width + '\">';\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.updated_canvas_event = function () {\n",
       "    // Tell IPython that the notebook contents must change.\n",
       "    IPython.notebook.set_dirty(true);\n",
       "    this.send_message('ack', {});\n",
       "    var fig = this;\n",
       "    // Wait a second, then push the new image to the DOM so\n",
       "    // that it is saved nicely (might be nice to debounce this).\n",
       "    setTimeout(function () {\n",
       "        fig.push_to_output();\n",
       "    }, 1000);\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._init_toolbar = function () {\n",
       "    var fig = this;\n",
       "\n",
       "    var toolbar = document.createElement('div');\n",
       "    toolbar.classList = 'btn-toolbar';\n",
       "    this.root.appendChild(toolbar);\n",
       "\n",
       "    function on_click_closure(name) {\n",
       "        return function (_event) {\n",
       "            return fig.toolbar_button_onclick(name);\n",
       "        };\n",
       "    }\n",
       "\n",
       "    function on_mouseover_closure(tooltip) {\n",
       "        return function (event) {\n",
       "            if (!event.currentTarget.disabled) {\n",
       "                return fig.toolbar_button_onmouseover(tooltip);\n",
       "            }\n",
       "        };\n",
       "    }\n",
       "\n",
       "    fig.buttons = {};\n",
       "    var buttonGroup = document.createElement('div');\n",
       "    buttonGroup.classList = 'btn-group';\n",
       "    var button;\n",
       "    for (var toolbar_ind in mpl.toolbar_items) {\n",
       "        var name = mpl.toolbar_items[toolbar_ind][0];\n",
       "        var tooltip = mpl.toolbar_items[toolbar_ind][1];\n",
       "        var image = mpl.toolbar_items[toolbar_ind][2];\n",
       "        var method_name = mpl.toolbar_items[toolbar_ind][3];\n",
       "\n",
       "        if (!name) {\n",
       "            /* Instead of a spacer, we start a new button group. */\n",
       "            if (buttonGroup.hasChildNodes()) {\n",
       "                toolbar.appendChild(buttonGroup);\n",
       "            }\n",
       "            buttonGroup = document.createElement('div');\n",
       "            buttonGroup.classList = 'btn-group';\n",
       "            continue;\n",
       "        }\n",
       "\n",
       "        button = fig.buttons[name] = document.createElement('button');\n",
       "        button.classList = 'btn btn-default';\n",
       "        button.href = '#';\n",
       "        button.title = name;\n",
       "        button.innerHTML = '<i class=\"fa ' + image + ' fa-lg\"></i>';\n",
       "        button.addEventListener('click', on_click_closure(method_name));\n",
       "        button.addEventListener('mouseover', on_mouseover_closure(tooltip));\n",
       "        buttonGroup.appendChild(button);\n",
       "    }\n",
       "\n",
       "    if (buttonGroup.hasChildNodes()) {\n",
       "        toolbar.appendChild(buttonGroup);\n",
       "    }\n",
       "\n",
       "    // Add the status bar.\n",
       "    var status_bar = document.createElement('span');\n",
       "    status_bar.classList = 'mpl-message pull-right';\n",
       "    toolbar.appendChild(status_bar);\n",
       "    this.message = status_bar;\n",
       "\n",
       "    // Add the close button to the window.\n",
       "    var buttongrp = document.createElement('div');\n",
       "    buttongrp.classList = 'btn-group inline pull-right';\n",
       "    button = document.createElement('button');\n",
       "    button.classList = 'btn btn-mini btn-primary';\n",
       "    button.href = '#';\n",
       "    button.title = 'Stop Interaction';\n",
       "    button.innerHTML = '<i class=\"fa fa-power-off icon-remove icon-large\"></i>';\n",
       "    button.addEventListener('click', function (_evt) {\n",
       "        fig.handle_close(fig, {});\n",
       "    });\n",
       "    button.addEventListener(\n",
       "        'mouseover',\n",
       "        on_mouseover_closure('Stop Interaction')\n",
       "    );\n",
       "    buttongrp.appendChild(button);\n",
       "    var titlebar = this.root.querySelector('.ui-dialog-titlebar');\n",
       "    titlebar.insertBefore(buttongrp, titlebar.firstChild);\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._remove_fig_handler = function (event) {\n",
       "    var fig = event.data.fig;\n",
       "    if (event.target !== this) {\n",
       "        // Ignore bubbled events from children.\n",
       "        return;\n",
       "    }\n",
       "    fig.close_ws(fig, {});\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._root_extra_style = function (el) {\n",
       "    el.style.boxSizing = 'content-box'; // override notebook setting of border-box.\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._canvas_extra_style = function (el) {\n",
       "    // this is important to make the div 'focusable\n",
       "    el.setAttribute('tabindex', 0);\n",
       "    // reach out to IPython and tell the keyboard manager to turn it's self\n",
       "    // off when our div gets focus\n",
       "\n",
       "    // location in version 3\n",
       "    if (IPython.notebook.keyboard_manager) {\n",
       "        IPython.notebook.keyboard_manager.register_events(el);\n",
       "    } else {\n",
       "        // location in version 2\n",
       "        IPython.keyboard_manager.register_events(el);\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._key_event_extra = function (event, _name) {\n",
       "    var manager = IPython.notebook.keyboard_manager;\n",
       "    if (!manager) {\n",
       "        manager = IPython.keyboard_manager;\n",
       "    }\n",
       "\n",
       "    // Check for shift+enter\n",
       "    if (event.shiftKey && event.which === 13) {\n",
       "        this.canvas_div.blur();\n",
       "        // select the cell after this one\n",
       "        var index = IPython.notebook.find_cell_index(this.cell_info[0]);\n",
       "        IPython.notebook.select(index + 1);\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_save = function (fig, _msg) {\n",
       "    fig.ondownload(fig, null);\n",
       "};\n",
       "\n",
       "mpl.find_output_cell = function (html_output) {\n",
       "    // Return the cell and output element which can be found *uniquely* in the notebook.\n",
       "    // Note - this is a bit hacky, but it is done because the \"notebook_saving.Notebook\"\n",
       "    // IPython event is triggered only after the cells have been serialised, which for\n",
       "    // our purposes (turning an active figure into a static one), is too late.\n",
       "    var cells = IPython.notebook.get_cells();\n",
       "    var ncells = cells.length;\n",
       "    for (var i = 0; i < ncells; i++) {\n",
       "        var cell = cells[i];\n",
       "        if (cell.cell_type === 'code') {\n",
       "            for (var j = 0; j < cell.output_area.outputs.length; j++) {\n",
       "                var data = cell.output_area.outputs[j];\n",
       "                if (data.data) {\n",
       "                    // IPython >= 3 moved mimebundle to data attribute of output\n",
       "                    data = data.data;\n",
       "                }\n",
       "                if (data['text/html'] === html_output) {\n",
       "                    return [cell, data, j];\n",
       "                }\n",
       "            }\n",
       "        }\n",
       "    }\n",
       "};\n",
       "\n",
       "// Register the function which deals with the matplotlib target/channel.\n",
       "// The kernel may be null if the page has been refreshed.\n",
       "if (IPython.notebook.kernel !== null) {\n",
       "    IPython.notebook.kernel.comm_manager.register_target(\n",
       "        'matplotlib',\n",
       "        mpl.mpl_figure_comm\n",
       "    );\n",
       "}\n"
      ],
      "text/plain": [
       "<IPython.core.display.Javascript object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/html": [
       "<img src=\"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAA+gAAAPoCAYAAABNo9TkAAAgAElEQVR4nOzda5DedXn/8d2FZDchITQRCCRycpjKYSAyjAKlxNAgKoeRokNRxIwcH4QCzii1WAJoISnWwWItOj2AlQK1gB1tBYFplDJTESMQKBASoiyJY6tpEoYSJOT6P7BuuyTcfLtJc/3yv16vme8M/Pbe3e+Tz4M3u9zbFwAAAEC6vuwLAAAAAAIdAAAAOkGgAwAAQAcIdAAAAOgAgQ4AAAAdINABAACgAwQ6AAAAdIBABwAAgA4Q6AAAANABAh0AAAA6QKADAABABwh0AAAA6ACBDgAAAB0g0AEAAKADBDoAAAB0gEAHAACADhDoAAAA0AECHQAAADpAoAMAAEAHCHQAAADoAIEOAAAAHSDQAQAAoAMEOgAAAHSAQAcAAIAOEOgAAADQAQIdAAAAOkCgAwAAQAcIdAAAAOgAgQ4AAAAdINABAACgAwQ6AAAAdIBABwAAgA4Q6AAAANABAh0AAAA6QKADAABABwh0AAAA6ACBDgAAAB0g0AEAAKADBDoAAAB0gEAHAACADhDoAAAA0AECHQAAADpAoAMAAEAHCHQAAADoAIEOAAAAHSDQAQAAoAMEOgAAAHSAQAcAAIAOEOgAAADQAQIdAAAAOkCgAwAAQAcIdAAAAOgAgQ4AAAAdINABAACgAwQ6AAAAdIBABwAAgA4Q6AAAANABAh0AAAA6QKADAABABwh0AAAA6ACBDgAAAB0g0AEAAKADBDoAAAB0gEAHAACADhDoAAAA0AECHQAAADpAoAMAAEAHCHQAAADoAIEOAAAAHSDQAQAAoAMEOgAAAHSAQAcAAIAOEOgAAADQAQIdAAAAOkCgAwAAQAcIdAAAAOgAgQ4AAAAdINABAACgAwQ6AAAAdIBABwAAgA4Q6AAAANABAh0AAAA6QKADAABABwh0AAAA6ACBDgAAAB0g0AEAAKADBDoAAAB0gEAHAACADhDoAAAA0AECHQAAADpAoAMAAEAHCHQAAADoAIEOAAAAHSDQAQAAoAMEOgAAAHSAQAcAAIAOEOgAAADQAQIdAAAAOkCgAwAAQAcIdAAAAOgAgQ4AAAAdINABAACgAwQ6AAAAdIBABwAAgA4Q6AAAANABAh0AAAA6QKADAABABwh0AAAA6ACBDgAAAB0g0AEAAKADBDoAAAB0gEAHAACADhDoAAAA0AECHQAAADpAoAMAAEAHCHQAAADoAIEOAAAAHSDQAQAAoAMEejGvvvpqDA8Px9q1a2PdunWO4ziO4ziO4xQ9a9eujeHh4Xj11VezM4X/ItCLGR4ejr6+PsdxHMdxHMdxnOjr64vh4eHsTOG/CPRi1q5dOzLC7P9i5ziO4ziO4zhO3vnVD+/Wrl2bnSn8F4FezLp166Kvry/WrVuXfRUAACCRNugegV6MEQIAABHaoIsEejFGCAAARGiDLhLoxRghAAAQoQ26SKAXY4QAAECENugigV6MEQIAABHaoIsEejFGCAAARGiDLhLoxRghAAAQoQ26SKAXY4QAAECENugigV6MEQIAABHaoIsEejFGCAAARGiDLhLoxRghAAAQoQ26SKAXY4QAAECENugigV6MEQIAABHaoIsEejFGCAAARGiDLhLoxRghAAAQoQ26SKAXY4QAAECENugigV6MEQIAABHaoIsEejFGCAAARGiDLhLoxRghAAAQoQ26SKAXY4QAAECENugigV6MEQIAABHaoIsEejFGCAAARGiDLhLoxRghAAAQoQ26SKAXY4QAAECENugigV6MEQIAABHaoIsEejFGCAAARGiDLhLoxRghAAAQoQ26SKAXY4QAAECENugigV6MEQIAABHaoIsEejFGCAAARGiDLhLoxRghAAAQoQ26SKAXY4QAAECENugigV6MEQIAABHaoIsEejFGCAAARGiDLhLoxRghAAAV7XvZNztzukIbdI9AL8YIAQCoKDvKBTotBHoxRggAQEXZUS7QaSHQizFCAAAqyo5ygU4LgV6MEQIAUFF2lAt0Wgj0YowQAICKsqNcoNNCoBdjhAAAVJQd5QKdFgK9GCMEAKCi7CgX6LQQ6MUYIQAAFWVHuUCnhUAvxggBAKgoO8oFOi0EejFGCABARdlRLtBpIdCLMUIAACrKjnKBTguBXowRAgBQUXaUC3RaCPRijBAAgIqyo1yg00KgF2OEAABUlB3lAp0WAr0YIwQAoKLsKBfotBDoxRghAAAVZUe5QKeFQC/GCAEAqCg7ygU6LQR6MUYIAEBF2VEu0Gkh0IsxQgAAKsqOcoFOC4FejBECAFBRdpQLdFoI9GKMEACAirKjXKDTQqAXY4QAAFSUHeUCnRYCvRgjBACgouwoF+i0EOjFGCEAABVlR7lAp4VAL8YIAQCoKDvKBTotBHoxRggAQEXZUS7QaSHQizFCAAAqyo5ygU4LgV6MEQIAUFF2lAt0Wgj0YowQAICKsqNcoNNCoBdjhAAAVJQd5QKdFgK9GCMEAKCi7CgX6LQQ6MUYIQAAFWVHuUCnhUAvxggBAKgoO8oFOi0EejFGCABARdlRLtBpIdCLMUIAACrKjnKBTguBXowRAgBQUXaUC3RaCPRijBAAgIqyo1yg00KgF2OEAABUlB3lAp0WAr0YIwQAoKLsKBfotBDoxRghAAAVZUe5QKeFQC/GCAEAqCg7ygU6LQR6MUYIAEBF2VEu0Gkh0IsxQgAAKsqOcoFOC4FejBECAFBRdpQLdFoI9GKMEACAirKjXKDTQqAXY4QAAFSUHeUCnRYCvRgjBACgouwoF+i0EOjFGCEAABVlR7lAp4VAL8YIAQCoKDvKBTotBHoxRggAQEXZUS7QaSHQizFCAAAqyo5ygU4LgV6MEQIAUFF2lAt0Wgj0YowQAICKsqNcoNNCoBdjhAAAVJQd5QKdFgK9GCMEAKCi7CgX6LQQ6MUYIQAAFWVHuUCnhUAvxggBAKgoO8oFOi0EejFGCABARdlRLtBpIdCLMUIAACrKjnKBTguBXowRAgBQUXaUC3RaCPRijBAAgIqyo1yg00KgF2OEAABUlB3lAp0WAr0YIwQAoKLsKBfotBDoxRghAAAVZUe5QKeFQC/GCAEAqCg7ygU6LQR6MUYIAEBF2VEu0Gkh0IsxQgAAKsqOcoFOC4FejBECAFBRdpQLdFoI9GKMEACAirKjXKDTQqAXY4QAAFSUHeUCnRYCvRgjBACgouwoF+i0EOjFGCEAABVlR7lAp4VAL8YIAQCoKDvKBTotBHoxRggAQEXZUS7QaSHQizFCAAAqyo5ygU4LgV6MEQIAUFF2lAt0Wgj0YowQAICKsqNcoNNCoBdjhAAAVJQd5QKdFgK9GCMEAKCi7CgX6LQQ6MUYIQAAFWVHuUCnhUAvxggBAKgoO8oFOi0EejFGCABARdlRLtBpIdCLMUIAACrKjnKBTguBXowRAgBQUXaUC3RaCPRijBAAgIqyo1yg00KgF2OEAABUlB3lAp0WAr0YIwQAoKLsKBfotBDoxRghAAAVZUe5QKeFQC/GCAEAqCg7ygU6LQR6MUYIAEBF2VEu0Gkh0IsxQgAAKsqOcoFOC4FejBECAFBRdpQLdFoI9GKMEACAirKjXKDTQqAXY4QAAFSUHeUCnRYCfTt55ZVX4vLLL4/99tsvhoaGYv/994+rrroqXn311ZHXbNq0KRYsWBB77bVXDA0NxezZs+Pxxx8f9XU2bNgQ8+fPj2nTpsXEiRPjlFNOieHh4eZ7GCEAABVlR7lAp4VA304+85nPxLRp0+Kb3/xmrFy5Mr72ta/FpEmT4vrrrx95zcKFC2Py5Mlxxx13xNKlS+OMM86IvfbaK9avXz/ymgsvvDBmzJgR9957byxZsiTmzJkThx9+eGzcuLHpHkYIAEBF2VEu0Gkh0LeTk046KT760Y+Oevbbv/3bcdZZZ0XEL396Pn369Fi4cOHIxzds2BBTpkyJG2+8MSIi1q5dG+PGjYvbbrtt5DWrVq2KgYGBuPvuu5vuYYQAAFSUHeUCnRYCfTu59tprY999942nn346IiIeeeSR2GOPPeJv/uZvIiJixYoV0dfXF0uWLBn1eaeeemqcffbZERFx//33R19fX6xZs2bUaw477LC44oortvh9N2zYEOvWrRs5w8PDRggAQDnZUS7QaSHQt5NNmzbF7/3e70V/f3/svPPO0d/fH9dcc83Ixx988MHo6+uLVatWjfq88847L971rndFRMQtt9wS48eP3+xrn3DCCXH++edv8fsuWLAg+vr6NjtGCABAJdlRLtBpIdC3k1tvvTVmzpwZt956azz22GPxla98JaZOnRo33XRTRPx3oK9evXrU55177rlx4oknRsTrB/rcuXPjggsu2OL39RN0AAAQ6Fsi0LtHoG8nM2fOjC984Qujnn3605+OX//1X4+I/7tfcX8tIwQAoKLsKBfotBDo28nUqVPji1/84qhn11xzTRx44IER8d9vErdo0aKRj7/88stbfJO422+/feQ1q1ev9iZxAADwBrKjXKDTQqBvJx/5yEdixowZI39m7c4774w3velN8YlPfGLkNQsXLowpU6bEnXfeGUuXLo0zzzxzi39mbebMmXHffffFkiVL4vjjj/dn1gAA4A1kR7lAp4VA307Wr18fF198ceyzzz4xNDQUBxxwQFx++eXx8ssvj7xm06ZNsWDBgpg+fXoMDg7GcccdF0uXLh31dV566aWYP39+TJ06NSZMmBAnn3xyPPfcc833MEIAACrKjnKBTguBXowRAgBQUXaUC3RaCPRijBAAgIqyo1yg00KgF2OEAABUlB3lAp0WAr0YIwQAoKLsKBfotBDoxRghAAAVZUe5QKeFQC/GCAEAqCg7ygU6LQR6MUYIAEBF2VEu0Gkh0IsxQgAAKsqOcoFOC4FejBECAFBRdpQLdFoI9GKMEACAirKjXKDTQqAXY4QAAFSUHeUCnRYCvRgjBACgouwoF+i0EOjFGCEAABVlR7lAp4VAL8YIAQCoKDvKBTotBHoxRggAQEXZUS7QaSHQizFCAAAqyo5ygU4LgV6MEQIAUFF2lAt0Wgj0YowQAICKsqNcoNNCoBdjhAAAVJQd5QKdFgK9GCMEAKCi7CgX6LQQ6MUYIQAAFWVHuUCnhUAvxggBAKgoO8oFOi0EejFGCABARdlRLtBpIdCLMUIAACrKjnKBTguBXowRAgBQUXaUC3RaCPRijBAAgIqyo1yg00KgF2OEAABUlB3lAp0WAr0YIwQAoKLsKBfotBDoxRghAAAVZUe5QKeFQC/GCAEAqCg7ygU6LQR6MUYIAEBF2VEu0Gkh0IsxQgAAKsqOcoFOC4FejBECAFBRdpQLdFoI9GKMEACAirKjXKDTQqAXY4QAAFSUHeUCnRYCvRgjBACgouwoF+i0EOjFGCEAABVlR7lAp4VAL8YIAQCoKDvKBTotBHoxRggAQEXZUS7QaSHQizFCAAAqyo5ygU4LgV6MEQIAUFF2lAt0Wgj0YowQAICKsqNcoNNCoBdjhAAAVJQd5QKdFgK9GCMEAKCi7CgX6LQQ6MUYIQAAFWVHuUCnhUAvxggBAKgoO8oFOi0EejFGCABARdlRLtBpIdCLMUIAACrKjnKBTguBXowRAgBQUXaUC3RaCPRijBAAgIqyo1yg00KgF2OEAABUlB3lAp0WAr0YIwQAoKLsKBfotBDoxRghAAAVZUe5QKeFQC/GCAEAqCg7ygU6LQR6MUYIAEBF2VEu0Gkh0IsxQgAAKsqOcoFOC4FejBECAFBRdpQLdFoI9GKMEACAirKjXKDTQqAXY4QAAFSUHeUCnRYCvRgjBACgouwoF+i0EOjFGCEAABVlR7lAp4VAL8YIAQCoKDvKBTotBHoxRggAQEXZUS7QaSHQizFCAAAqyo5ygU4LgV6MEQIAUFF2lAt0Wgj0YowQAICKsqNcoNNCoBdjhAAAVJQd5QKdFgK9GCMEAKCi7CgX6LQQ6MUYIQAAFWVHuUCnhUAvxggBAKgoO8oFOi0EejFGCABARdlRLtBpIdCLMUIAACrKjnKBTguBXowRAgBQUXaUC3RaCPRijBAAgIqyo1yg00KgF2OEAABUlB3lAp0WAr0YIwQAoKLsKBfotBDoxRghAAAVZUe5QKeFQC/GCAEAqCg7ygU6LQR6MUYIAEBF2VEu0Gkh0IsxQgAAKsqOcoFOC4FejBECAFBRdpQLdFoI9GKMEACAirKjXKDTQqAXY4QAAFSUHeUCnRYCvRgjBACgouwoF+i0EOjFGCEAABVlR7lAp4VAL8YIAQCoKDvKBTotBHoxRggAQEXZUS7QaSHQizFCAAAqyo5ygU4LgV6MEQIAUFF2lAt0Wgj0YowQAICKsqNcoNNCoBdjhAAAVJQd5QKdFgK9GCMEAKCi7CgX6LQQ6MUYIQAAFWVHuUCnhUAvxggBAKgoO8oFOi0EejFGCABARdlRLtBpIdCLMUIAACrKjnKBTguBXowRAgBQUXaUC3RaCPRijBAAgIqyo1yg00KgF2OEAABUlB3lAp0WAr0YIwQAoKLsKBfotBDoxRghAAAVZUe5QKeFQC/GCAEAqCg7ygU6LQR6MUYIAEBF2VEu0Gkh0IsxQgAAKsqOcoFOC4FejBECAFBRdpQLdFoI9GKMEACAirKjXKDTQqAXY4QAAFSUHeUCnRYCvRgjBACgouwoF+i0EOjFGCEAABVlR7lAp4VAL8YIAQCoKDvKBTotBHoxRggAQEXZUS7QaSHQizFCAAAqyo5ygU4LgV6MEQIAUFF2lAt0Wgj0YowQAICKsqNcoNNCoBdjhAAAVJQd5QKdFgK9GCMEAKCi7CgX6LQQ6MUYIQAAFWVHuUCnhUAvxggBAKgoO8oFOi0EejFGCABARdlRLtBpIdCLMUIAACrKjnKBTguBXowRAgBQUXaUC3RaCPRijBAAgIqyo1yg00KgF2OEAABUlB3lAp0WAr0YIwQAoKLsKBfotBDoxRghAAAVZUe5QKeFQC/GCAEAqCg7ygU6LQR6MUYIAEBF2VEu0Gkh0IsxQgAAKsqOcoFOC4FejBECAFBRdpQLdFoI9GKMEACAirKjXKDTQqAXY4QAAFSUHeUCnRYCvRgjBACgouwoF+i0EOjFGCEAABVlR7lAp4VAL8YIAQCoKDvKBTotBHoxRggAQEXZUS7QaSHQizFCAAAqyo5ygU4LgV6MEQIAUFF2lAt0Wgj0YowQAICKsqNcoNNCoBdjhAAAVJQd5QKdFgK9GCMEAKCi7CgX6LQQ6MUYIQAAFWVHuUCnhUAvxggBAKgoO8oFOi0EejFGCABARdlRLtBpIdCLMUIAACrKjnKBTguBXowRAgBQUXaUC3RaCPRijBAAgIqyo1yg00KgF2OEAABUlB3lAp0WAr0YIwQAoKLsKBfotBDo29Hzzz8fH/rQh2Lq1KkxYcKEOPzww+Phhx8e+fimTZtiwYIFsddee8XQ0FDMnj07Hn/88VFfY8OGDTF//vyYNm1aTJw4MU455ZQYHh5uvoMRAgBQUXaUC3RaCPTtZM2aNbHvvvvGvHnz4nvf+16sXLky7rvvvli+fPnIaxYuXBiTJ0+OO+64I5YuXRpnnHFG7LXXXrF+/fqR11x44YUxY8aMuPfee2PJkiUxZ86cOPzww2Pjxo1N9zBCAAAqyo5ygU4Lgb6dXHbZZXHssce+7sc3bdoU06dPj4ULF44827BhQ0yZMiVuvPHGiIhYu3ZtjBs3Lm677baR16xatSoGBgbi7rvvbrqHEQIAUFF2lAt0Wgj07eSggw6KSy65JN7//vfH7rvvHrNmzYovf/nLIx9fsWJF9PX1xZIlS0Z93qmnnhpnn312RETcf//90dfXF2vWrBn1msMOOyyuuOKKpnsYIQAAFWVHuUCnhUDfTgYHB2NwcDA++clPxpIlS+LGG2+MoaGhuPnmmyMi4sEHH4y+vr5YtWrVqM8777zz4l3veldERNxyyy0xfvz4zb72CSecEOeff/4Wv++GDRti3bp1I2d4eNgIAQAoJzvKBTotBPp2Mm7cuDj66KNHPbvoooviqKOOioj/DvTVq1ePes25554bJ554YkS8fqDPnTs3Lrjggi1+3wULFkRfX99mxwgBAKgkO8oFOi0E+nayzz77xDnnnDPq2Re/+MXYe++9I+L/7lfc/QQdAAAE+pYI9O4R6NvJmWeeudmbxF1yySUjP1X/1ZvELVq0aOTjL7/88hbfJO72228fec3q1au9SRwAALyB7CgX6LQQ6NvJQw89FDvvvHP84R/+YTzzzDNxyy23xMSJE+OrX/3qyGsWLlwYU6ZMiTvvvDOWLl0aZ5555hb/zNrMmTPjvvvuiyVLlsTxxx/vz6wBAMAbyI5ygU4Lgb4dfeMb34hDDz00BgcH461vfeuod3GP+OVP0RcsWBDTp0+PwcHBOO6442Lp0qWjXvPSSy/F/PnzY+rUqTFhwoQ4+eST47nnnmu+gxECAFBRdpQLdFoI9GKMEACAirKjXKDTQqAXY4QAAFSUHeUCnRYCvRgjBACgouwoF+i0EOjFGCEAABVlR7lAp4VAL8YIAQCoKDvKBTotBHoxRggAQEXZUS7QaSHQizFCAAAqyo5ygU4LgV6MEQIAUFF2lAt0Wgj0Hr71rW/FAw88MPLvX/jCF+Lwww+PM888M9asWZN4s7EzQgAAKsqOcoFOC4Hew6GHHhr/8A//EBERjz32WAwODsYnP/nJeMc73hHz5s1Lvt3YGCEAABVlR7lAp4VA72GXXXaJlStXRkTEggUL4vTTT4+IiB/84Aex5557Jt5s7IwQAICKsqNcoNNCoPfwa7/2a/HEE09ERMRv/MZvxJe+9KWIiFi5cmVMmDAh82pjZoQAAFSUHeUCnRYCvYdTTjklTjzxxLj66qtj3Lhx8fzzz0dExD333BMHHnhg8u3GxggBAKgoO8oFOi0Eeg8//vGP46STTorDDjss/vzP/3zk+SWXXBIXXXRR4s3GzggBAKgoO8oFOi0EejFGCABARdlRLtBpIdDfwPLly+Pyyy+P3/md34mf/vSnEfHLP7/2+OOPJ99sbIwQAICKsqNcoNNCoPewePHimDBhQsydOzfGjx8fK1asiIiIRYsWjbyj+47GCAEAqCg7ygU6LQR6D0cddVT88R//cURETJo0aSTQH3roodh7770zrzZmRggAQEXZUS7QaSHQe9hll13i2WefjYjRgb5y5coYHBzMvNqYGSEAABVlR7lAp4VA72HGjBnx4IMPRsToQL/zzjvjgAMOyLzamBkhAAAVZUe5QKeFQO/h4x//eBx77LHxk5/8JCZPnhzPPPNM/PM//3MccMABceWVV2Zfb0yMEACAirKjXKDTQqD38Itf/CI++MEPxsDAQPT398e4ceNiYGAgzjrrrNi4cWP29cbECAEAqCg7ygU6LQR6gxUrVsTXvva1uP3222PZsmXZ19kqRggAQEXZUS7QaSHQizFCAAAqyo5ygU4Lgd7D6aefHtdee+1mz//oj/4o3v/+9yfcaOsZIQAAFWVHuUCnhUDv4U1velM89thjmz1/7LHHYo899ki40dYzQgAAKsqOcoFOC4Hew9DQUDz11FObPX/yySdjaGgo4UZbzwgBAKgoO8oFOi0Eeg9HHnlkXHXVVZs9X7BgQRxxxBEJN9p6RggAQEXZUS7QaSHQe/j7v//72HnnnePss8+Om266KW666ab48Ic/HDvvvHPcdddd2dcbEyMEAKCi7CgX6LQQ6G/gm9/8ZhxzzDExceLEmDZtWsyZMycWL16cfa0xM0IAACrKjnKBTguBXowRAgBQUXaUC3RaCPQGL7/8cgwPD8ePf/zjUWdHZIQAAFSUHeUCnRYCvYdly5bFscceGwMDA6NOf39/DAwMZF9vTIwQAICKsqNcoNNCoPdwzDHHxHHHHRf/+I//GD/84Q/jkUceGXV2REYIAEBF2VEu0Gkh0HuYOHFiPPnkk9nX2KaMEACAirKjXKDTQqD3cOSRR8YDDzyQfY1tyggBAKgoO8oFOi0Eeg/3339/HH300fFP//RP8bOf/SzWrVs36uyIjBAAgIqyo1yg00Kg99Df3z/yhnDeJA4AAHZc2VEu0Gkh0HtYvHhxz7MjMkIAACrKjnKBTguBXowRAgBQUXaUC3RaCPQ38N3vfjc+9KEPxdFHHx3PP/98RER85Stf2WHfPM4IAQCoKDvKBTotBHoPf/d3fxcTJkyIc889NwYHB2PFihUREfGnf/qn8Z73vCf5dmNjhAAAVJQd5QKdFgK9h1mzZsXNN98cERGTJk0aCfQf/vCHseeee2ZebcyMEACAirKjXKDTQqD3MGHChFi5cmVEjA70FStWxODgYOLNxs4IAQCoKDvKBTotBHoPBxxwQNx7770RMTrQb7755jjooIMyrzZmRggAQEXZUS7QaSHQe1i0aFEcfPDB8S//8i8xefLkeOCBB+KrX/1q7L777nHDDTdkX29MjBAAgIqyo1yg00Kgv4Hf//3fjwkTJkR/f3/09/fH0NBQfOpTn8q+1pgZIQAAFWVHuUCnhUB/HRs3bozFixfHz3/+83jxxRfj+9//fnzve9+LF154IftqW8UIAQCoKDvKBTotBHoPg4OD8eyzz2ZfY5syQgAAKsqOcoFOC4Hew5FHHhn33Xdf9jW2KSMEAKCi7CgX6LQQ6D3cc889MWvWrPjGN74Rq1evjnXr1o06OyIjBACgouwoF+i0EOg9/OqN4fr7+2NgYGDk/Orfd0RGCABARdlRLtBpIdB7WLx4cc+zIzJCAAAqyo5ygU4LgV6MEQIAUFF2lAt0Wgj0Hr7zne/0PDsiIwQAoKLsKBfotBDoPfzP/wd9S/8v+o7ICAEAqCg7ygU6LQR6D2vXrh11/v3f/z2+/e1vxzve8Y4d9s+vGSEAABVlR7lAp4VAH4PvfOc7ccQRR2RfY0yMEACAirKjXKDTQqCPwfmixNUAACAASURBVL/+67/GLrvskn2NMTFCAAAqyo5ygU4Lgd7Do48+Ouo88sgj8a1vfStmz54dxxxzTPb1xsQIAQCoKDvKBTotBHoPv3pTuNe+UdzRRx8dTz75ZPb1xsQIAQCoKDvKBTotBHoPP/rRj0ad5557Ll566aXsa20VIwQAoKLsKBfotBDoxRghAAAVZUe5QKeFQO/hoosuis9//vObPb/hhhvi4osvTrjR1jNCAAAqyo5ygU4Lgd7D3nvvHQ8//PBmz3/wgx/EjBkzEm609YwQAICKsqNcoNNCoPcwODgYzzzzzGbPn3nmmRgcHEy40dYzQgAAKsqOcoFOC4HewyGHHBI33HDDZs//5E/+JA466KCEG209IwQAoKLsKBfotBDoPfzFX/xFTJgwIa644opYvHhxLF68OP7gD/4gJk6cGF/+8pezrzcmRggAQEXZUS7QaSHQ38AXv/jFmDFjxsjfQN9///3j5ptvzr7WmBkhAAAVZUe5QKeFQG/0b//2b/HCCy9kX2OrGSEAABVlR7lAp4VA7+HZZ5+NZcuWbfZ82bJlsXLlyu1/oW3ACAEAqCg7ygU6LQR6D8cdd1zcdNNNmz3/67/+65g9e/b2v9A2YIQAAFSUHeUCnRYCvYfJkye/7p9ZmzJlSsKNtp4RAgBQUXaUC3RaCPQedt1111iyZMlmzx9++OGYNGlSwo22nhECAFBRdpQLdFoI9B5OOumk+MAHPhAbN24cebZx48Y4/fTT493vfnfizcbOCAEAqCg7ygU6LQR6D0888URMmzYt3vKWt8S8efNi3rx58Za3vCV23333WLp0afb1xsQIAQCoKDvKBTotBPobWLVqVXzyk5+M9773vXH66afHVVddFT//+c+zrzVmRggAQEXZUS7QaSHQizFCAAAqyo5ygU4Lgf4G/uM//iM++9nPxjnnnBPnnntufO5zn4u1a9dmX2vMjBAAgIqyo1yg00Kg9/D9738/pk6dGjNmzIjTTjst3ve+98XMmTNj2rRp8YMf/CD7emNihAAAVJQd5QKdFgK9h2OPPTbmzZsXr7zyysizV155JT7ykY/Eb/7mbybebOyMEACAirKjXKDTQqD3MDQ0FE8++eRmz5944omYMGFCwo22nhECAFBRdpQLdFoI9B722GOPuOeeezZ7fvfdd8cee+yRcKOtZ4QAAFSUHeUCnRYCvYeLLrooZs6cGbfddls899xzMTw8HLfeemvMnDkzLr744uzrjYkRAgBQUXaUC3RaCPQeXn755fjd3/3dGD9+fAwMDMTAwEAMDg7GJZdcEhs2bMi+3pgYIQAAFWVHuUCnhUBv8OKLL8Zjjz0Wjz76aLz44ovZ19kqRggAQEXZUS7QaSHQizFCAAAqyo5ygU4LgV6MEQIAUFF2lAt0Wgj0YowQAICKsqNcoNNCoBdjhAAAVJQd5QKdFgL9Nd72trfFmjVrIiLiqquu2uHfFO61jBAAgIqyo1yg00Kgv8bQ0FAMDw9HRMTAwED89Kc/Tb7RtmWEAABUlB3lAp0WAv01jjrqqJg7d25ceeWV0d/fHx//+Mfjqquu2uLZERkhAAAVZUe5QKeFQH+Np556Ks4444w48sgjY2BgIA499NCYNWvWZudtb3tb9lXHxAgBAKgoO8oFOi0Eeg/9/f1+xR0AAP4/kB3lAp0WAr0YIwQAoKLsKBfotBDob2D58uUxf/78+K3f+q2YO3duXHTRRbF8+fLsa42ZEQIAUFF2lAt0Wgj0Hu6+++4YP358vP3tb49LL700Lrnkknj7298eg4OD8e1vfzv7emNihAAAVJQd5QKdFgK9h1mzZsVll1222fPLLrvMm8QBAMAOJDvKBTotBHoPg4ODsWzZss2eP/300zE4OJhwo61nhAAAVJQd5QKdFgK9h5kzZ8bf/u3fbvb89ttvjze/+c0JN9p6RggAQEXZUS7QaSHQe7jqqqtit912i4ULF8Z3v/vdeOCBB+Laa6+N3XbbLT796U9nX29MjBAAgIqyo1yg00Kg97Bp06b43Oc+FzNmzIj+/v7o7++PGTNmxPXXXx+bNm3Kvt6YGCEAABVlR7lAp4VAb7R+/fpYv3599jW2mhECAFBRdpQLdFoI9GKMEACAirKjXKDTQqAXY4QAAFSUHeUCnRYCvRgjBACgouwoF+i0EOjFGCEAABVlR7lAp4VAfx2/+MUv4p3vfGc8/fTT2VfZpowQAICKsqNcoNNCoPfwpje9KZYtW5Z9jW3KCAEAqCg7ygU6LQR6Dx/72Mfisssuy77GNmWEAABUlB3lAp0WAr2H+fPnx6677hpHHHFEnH/++XHppZeOOjsiIwQAoKLsKBfotBDoPbzzne983TNnzpzs642JEQIAUFF2lAt0Wgj0YowQAICKsqNcoNNCoDd45pln4u67747//M//jIiITZs2Jd9o7IwQAICKsqNcoNNCoPfws5/9LI4//vjo7++PgYGBWLFiRUREfPSjH42PfexjybcbGyMEAKCi7CgX6LQQ6D18+MMfjhNPPDGGh4dj0qRJI4F+zz33xMEHH5x8u7ExQgAAKsqOcoFOC4Hew5577hmPPPJIRMSoQH/22Wdjl112ybzamBkhAAAVZUe5QKeFQO9h0qRJsWzZspF//lWgP/TQQzF16tTMq42ZEQIAUFF2lAt0Wgj0Ht773vfGpz71qYj4ZaA/++yz8eqrr8YHPvCBOP3005NvNzZGCABARdlRLtBpIdB7eOKJJ2L33XePd7/73TF+/Ph4//vfHwcddFDsueeesXz58uzrjYkRAgBQUXaUC3RaCPQ38JOf/CSuuOKKOOmkk+I973lPXH755bF69ersa42ZEQIAUFF2lAt0Wgj0YowQAICKsqNcoNNCoL+BNWvWxHXXXRcf/ehH45xzzonPfvaz8fOf/zz7WmNmhAAAVJQd5QKdFgK9h8WLF8eUKVPizW9+c5x22mlx2mmnxT777BO77rprLF68OPt6Y2KEAABUlB3lAp0WAr2HQw45JM4777zYuHHjyLONGzfG+eefH4ccckjizcbOCAEAqCg7ygU6LQR6D0NDQ/HUU09t9vypp56KoaGhhBttPSMEAKCi7CgX6LQQ6D0cc8wxcdddd232/K677oqjjjoq4UZbzwgBAKgoO8oFOi0E+ms8+uijI+e2226LffbZJ6677rp44IEH4oEHHojrrrsu9ttvv7jtttuyrzomRggAQEXZUS7QaSHQX6O/vz8GBgaiv7+/5xkYGMi+6pgYIQAAFWVHuUCnhUB/jR/96EfNZ0dkhAAAVJQd5QKdFgK9GCMEAKCi7CgX6LQQ6G/g+eefj9tvvz1uuOGG+PznPz/q7IiMEACAirKjXKDTQqD38Jd/+Zcxfvz4mDRpUuy7776x3377jZz9998/+3pjYoQAAFSUHeUCnRYCvYeZM2fGZz7zmXj11Vezr7LNGCEAABVlR7lAp4VA72Hq1KmxfPny7GtsU0YIAEBF2VEu0Gkh0Hv4+Mc/Htdee232NbYpIwQAoKLsKBfotBDoPWzcuDHe/e53x+zZs2P+/Plx6aWXjjo7IiMEAKCi7CgX6LQQ6D1cffXV0d/fH29961tj9uzZ8c53vnPkzJkzJ/t6Y2KEAABUlB3lAp0WAr2H3XbbLf7qr/4q+xrblBECAFBRdpQLdFoI9B723HPPWLZsWfY1tikjBACgouwoF+i0EOg9XHPNNXHRRRdlX2ObMkIAACrKjnKBTguB3sP73ve+2HXXXWP//fePk08+OU477bRRZ0dkhAAAVJQd5QKdFgK9h3nz5vU8OyIjBACgouwoF+i0EOjFGCEAABVlR7lAp4VAL8YIAQCoKDvKBTotBHoP++23X+y///6ve3ZERggAQEXZUS7QaSHQe7j++utHneuuuy4++MEPxtSpU+Paa6/Nvt6YGCEAABVlR7lAp4VAH4MvfOEL3iQOAAB2INlRLtBpIdDHYMWKFTF58uQxf/4111wTfX19cfHFF48827RpUyxYsCD22muvGBoaitmzZ8fjjz8+6vM2bNgQ8+fPj2nTpsXEiRPjlFNOieHh4f/V9zZCAAAqyo5ygU4LgT4GixYtin333XdMn/vQQw/FfvvtF4cddtioQF+4cGFMnjw57rjjjli6dGmcccYZsddee8X69etHXnPhhRfGjBkz4t57740lS5bEnDlz4vDDD4+NGzc2f38jBACgouwoF+i0EOg9zJo1K972treNnFmzZsX06dNjp512ii996Uv/66/3wgsvxIEHHhj33ntvzJ49eyTQN23aFNOnT4+FCxeOvHbDhg0xZcqUuPHGGyMiYu3atTFu3Li47bbbRl6zatWqGBgYiLvvvrv5DkYIAEBF2VEu0Gkh0Hu48sorR52rr746/uzP/iyefPLJMX29s88+Oy655JKIiFGBvmLFiujr64slS5aMev2pp54aZ599dkRE3H///dHX1xdr1qwZ9ZrDDjssrrjiitf9nhs2bIh169aNnOHhYSMEAKCc7CgX6LQQ6NvJrbfeGoceemi89NJLETE60B988MHo6+uLVatWjfqc8847L971rndFRMQtt9wS48eP3+zrnnDCCXH++ee/7vddsGBB9PX1bXaMEACASrKjXKDTQqBvB88991zsscce8cgjj4w821Kgr169etTnnXvuuXHiiSdGxOsH+ty5c+OCCy543e/tJ+gAACDQt0Sgd49A34L+/v4YGBjoeXbaaafmr3fXXXdFX19f7LTTTiOnr68v+vv7Y6eddorly5f/n/2K+2sZIQAAFWVHuUCnhUDfgq9//euvez7xiU/EhAkTYmhoqPnrrV+/PpYuXTrqHHnkkXHWWWfF0qVLR94kbtGiRSOf8/LLL2/xTeJuv/32kdesXr3am8QBAECD7CgX6LQQ6I2efPLJeN/73hc77bRTnH322fHjH/94q77e//wV94hf/pm1KVOmxJ133hlLly6NM888c4t/Zm3mzJlx3333xZIlS+L444/3Z9YAAKBBdpQLdFoI9DewatWqOPfcc2PcuHFx8sknx9KlS7fJ131toG/atCkWLFgQ06dPj8HBwTjuuOM2+14vvfRSzJ8/P6ZOnRoTJkyIk08+OZ577rn/1fc1QgAAKsqOcoFOC4H+OtauXTvy6+xHH310fPe7382+0jZhhAAAVJQd5QKdFgJ9CxYtWhRTp06Ngw8+OL7+9a9nX2ebMkIAACrKjnKBTguBvgX9/f0xceLEOPXUU+O000573bMjMkIAACrKjnKBTguBvgUf+chHYt68eW94dkRGCABARdlRLtBpIdCLMUIAACrKjnKBTguBXowRAgBQUXaUC3RaCPRijBAAgIqyo1yg00KgF2OEAABUlB3lAp0WAr0YIwQAoKLsKBfotBDoxRghAAAVZUe5QKeFQC/GCAEAqCg7ygU6LQR6MUYIAEBF2VEu0Gkh0IsxQgAAKsqOcoFOC4FejBECAFBRdpQLdFoI9GKMEACAirKjXKDTQqAXY4QAAFSUHeUCnRYCvRgjBACgouwoF+i0EOjFGCEAABVlR7lAp4VAL8YIAQCoKDvKBTotBHoxRggAQEXZUS7QaSHQizFCAAAqyo5ygU4LgV6MEQIAUFF2lAt0Wgj0YowQAICKsqNcoNNCoBdjhAAAVJQd5QKdFgK9GCMEAKCi7CgX6LQQ6MUYIQAAFWVHuUCnhUAvxggBAKgoO8oFOi0EejFGCABARdlRLtBpIdCLMUIAACrKjnKBTguBXowRAgBQUXaUC3RaCPRijBAAgIqyo1yg00KgF2OEAABUlB3lAp0WAr0YIwQAoKLsKBfotBDoxRghAAAVZUe5QKeFQC/GCAEAqCg7ygU6LQR6MUYIAEBF2VEu0Gkh0IsxQgAAKsqOcoFOC4FejBECAFBRdpQLdFoI9GKMEACAirKjXKDTQqAXY4QAAFSUHeUCnRYCvRgjBACgouwoF+i0EOjFGCEAABVlR7lAp4VAL8YIAQCoKDvKBTotBHoxRggAQEXZUS7QaSHQizFCAAAqyo5ygU4LgV6MEQIAUFF2lAt0Wgj0YowQAICKsqNcoNNCoBdjhAAAVJQd5QKdFgK9GCMEAKCi7CgX6LQQ6MUYIQAAFWVHuUCnhUAvxggBAKgoO8oFOi0EejFGCABARdlRLtBpIdCLMUIAACrKjnKBTguBXowRAgBQUXaUC3RaCPRijBAAgIqyo1yg00KgF2OEAABUlB3lAp0WAr0YIwQAoKLsKBfotBDoxRghAAAVZUe5QKeFQC/GCAEAqCg7ygU6LQR6MUYIAEBF2VEu0Gkh0IsxQgAAKsqOcoFOC4FejBECAFBRdpQLdFoI9GKMEACAirKjXKDTQqAXY4QAAFSUHeUCnRYCvRgjBACgouwoF+i0EOjFGCEAABVlR7lAp4VAL8YIAQCoKDvKBTotBHoxRggAQEXZUS7QaSHQizFCAAAqyo5ygU4LgV6MEQIAUFF2lAt0Wgj0YowQAICKsqNcoNNCoBdjhAAAVJQd5QKdFgK9GCMEAKCi7CgX6LQQ6MUYIQAAFWVHuUCnhUAvxggBAKgoO8oFOi0EejFGCABARdlRLtBpIdCLMUIAACrKjnKBTguBXowRAgBQUXaUC3RaCPRijBAAgIqyo1yg00KgF2OEAABUlB3lAp0WAr0YIwQAoKLsKBfotBDoxRghAAAVZUe5QKeFQC/GCAEAqCg7ygU6LQR6MUYIAEBF2VEu0Gkh0IsxQgAAKsqOcoFOC4FejBECAFBRdpQLdFoI9GKMEACAirKjXKDTQqAXY4QAAFSUHeUCnRYCvRgjBACgouwoF+i0EOjFGCEAABVlR7lAp4VAL8YIAQCoKDvKBTotBHoxRggAQEXZUS7QaSHQizFCAAAqyo5ygU4LgV6MEQIAUFF2lAt0Wgj0YowQAICKsqNcoNNCoBdjhAAAVJQd5QKdFgK9GCMEAKCi7CgX6LQQ6MUYIQAAFWVHuUCnhUAvxggBAKgoO8oFOi0EejFGCABARdlRLtBpIdCLMUIAACrKjnKBTguBXowRAgBQUXaUC3RaCPRijBAAgIqyo1yg00KgF2OEAABUlB3lAp0WAr0YIwQAoKLsKBfotBDoxRghAAAVZUe5QKeFQC/GCAEAqCg7ygU6LQR6MUYIAEBF2VEu0Gkh0IsxQgAAKsqOcoFOC4FejBECAFBRdpQLdFoI9GKMEACAirKjXKDTQqAXY4QAAFSUHeUCnRYCvRgjBACgouwoF+i0EOjFGCEAABVlR7lAp4VAL8YIAQCoKDvKBTotBHoxRggAQEXZUS7QaSHQizFCAAAqyo5ygU4LgV6MEQIAUFF2lAt0Wgj0YowQAICKsqNcoNNCoBdjhAAAVJQd5QKdFgK9GCMEAKCi7CgX6LQQ6MUYIQAAFWVHuUCnhUAvxggBAKgoO8oFOi0EejFGCABARdlRLtBpIdCLMUIAACrKjnKBTguBXowRAgBQUXaUC3RaCPRijBAAgIqyo1yg00KgF2OEAABUlB3lAp0WAr0YIwQAoKLsKBfotBDoxRghAAAVZUe5QKeFQC/GCAEAqCg7ygU6LQR6MUYIAEBF2VEu0Gkh0IsxQgAAKsqOcoFOC4FejBECAFBRdpQLdFoI9GKMEACAirKjXKDTQqAXY4T8v/buPUiq8k78cAvMDBeBBUQBUcQUKm7UhFLiVhKpXfGSEtEY16xZlU0WFUTFyxrZuBV0Sy0viVEj5mJcdSteEkWtrfVSUSN4S3RXUEFEvKASZRNFua2iK3x/f/CbKWbomX57AnMOvs9TdarkOH36Pd3ve6Y/TNMDAJCjoqNcoJNCoGfGIgQAIEdFR7lAJ4VAz4xFCABAjoqOcoFOCoGeGYsQAIAcFR3lAp0UAj0zFiEAADkqOsoFOikEemYsQgAAclR0lAt0Ugj0zFiEAADkqOgoF+ikEOiZsQgBAMhR0VEu0Ekh0DNjEQIAkKOio1ygk0KgZ8YiBAAgR0VHuUAnhUDPjEUIAECOio5ygU4KgZ4ZixAAgBwVHeUCnRQCPTMWIQAAOSo6ygU6KQR6F7n00ktj//33j+233z4GDx4cRx11VCxevLjV12zYsCFmzpwZQ4cOjZ49e8a4ceNi4cKFrb5m3bp1cfrpp8egQYOid+/eceSRR8ayZcuSx2ERAgCQo6KjXKCTQqB3kcMOOyxuuummWLhwYTz33HNxxBFHxK677hpr165t+ZrLLrss+vbtG7Nnz44FCxbEN7/5zRg6dGisXr265WumTJkSO++8czz00EMxb968+Ou//uvYb7/94tNPP00ah0UIAECOio5ygU4KgV6QP/3pT1GpVGLu3LkRsfGn50OGDInLLrus5WvWrVsX/fv3j5/+9KcREbFy5cpoaGiIO+64o+Vr3n777ejWrVs8+OCDSfdrEQIAkKOio1ygk0KgF+SVV16JSqUSCxYsiIiI1157LSqVSsybN6/V102cODFOOumkiIh45JFHolKpxPvvv9/qa/bdd9/4/ve/n3S/FiEAADkqOsoFOikEegE2bNgQRx55ZHzlK19p2ffkk09GpVKJt99+u9XXnnzyyXHooYdGRMStt94ajY2Nmx3vkEMOiVNOOaXqfa1bty5WrVrVsi1btswiBAAgO0VHuUAnhUAvwGmnnRYjRoxo9eFuzYH+zjvvtPrayZMnx2GHHRYR7Qf6+PHj49RTT616XzNnzoxKpbLZZhECAJCToqNcoJNCoHex008/PYYPHx6vv/56q/1b6y3ufoIOAAACvRqBXj4CvYts2LAhpk2bFsOGDYslS5ZU/f9DhgyJyy+/vGXfxx9/XPVD4n71q1+1fM0777zjQ+IAAKCGoqNcoJNCoHeRqVOnRv/+/WPOnDmxfPnylu3DDz9s+ZrLLrss+vfvH3fffXcsWLAgjj/++Kq/Zm348OHx8MMPx7x58+Jv/uZv/Jo1AACooegoF+ikEOhdpNq/A69UKnHTTTe1fM2GDRti5syZMWTIkGhqaoqDDjqo5VPem3300Udx+umnx8CBA6NXr14xYcKEeOutt5LHYRECAJCjoqNcoJNCoGfGIgQAIEdFR7lAJ4VAz4xFCABAjoqOcoFOCoGeGYsQAIAcFR3lAp0UAj0zFiEAADkqOsoFOikEemYsQgAAclR0lAt0Ugj0zFiEAADkqOgoF+ikEOiZsQgBAMhR0VEu0Ekh0DNjEQIAkKOio1ygk0KgZ8YiBAAgR0VHuUAnhUDPjEUIAECOio5ygU4KgZ4ZixAAgBwVHeUCnRQCPTMWIQAAOSo6ygU6KQR6ZixCAAByVHSUC3RSCPTMWIQAAOSo6CgX6KQQ6JmxCAEAyFHRUS7QSSHQM2MRAgCQo6KjXKCTQqBnxiIEACBHRUe5QCeFQM+MRQgAQI6KjnKBTgqBnhmLEACAHBUd5QKdFAI9MxYhAAA5KjrKBTopBHpmLEIAAHJUdJQLdFII9MxYhAAA5KjoKBfopBDombEIAQDIUdFRLtBJIdAzYxECAJCjoqNcoJNCoGfGIgQAIEdFR7lAJ4VAz4xFCABAjoqOcoFOCoGeGYsQAIAcFR3lAp0UAj0zFiEAADkqOsoFOikEemYsQgAAclR0lAt0Ugj0zFiEAADkqOgoF+ikEOiZsQgBAMhR0VEu0Ekh0DNjEQIAkKOio1ygk0KgZ8YiBAAgR0VHuUAnhUDPjEUIAECOio5ygU4KgZ4ZixAAgBwVHeUCnRQCPTMWIQAAOSo6ygU6KQR6ZixCAAByVHSUC3RSCPTMWIQAAOSo6CgX6KQQ6JmxCAEAyFHRUS7QSSHQM2MRAgCQo6KjXKCTQqBnxiIEACBHRUe5QCeFQM+MRQgAQI6KjnKBTgqBnhmLEACAHBUd5QKdFAI9MxYhAAA5KjrKBTopBHpmLEIAAHJUdJQLdFII9MxYhAAA5KjoKBfopBDombEIAQDIUdFRLtBJIdAzYxECAJCjoqNcoJNCoGfGIgQAIEdFR7lAJ4VAz4xFCABAjoqOcoFOCoGeGYsQAIAcFR3lAp0UAj0zFiEAADkqOsoFOikEemYsQgAAclR0lAt0Ugj0zFiEAADkqOgoF+ikEOiZsQgBAMhR0VEu0Ekh0DNjEQIAkKOio1ygk0KgZ8YiBAAgR0VHuUAnhUDPjEUIAECOio5ygU4KgZ4ZixAAgBwVHeUCnRQCPTMWIQAAOSo6ygU6KQR6ZixCAAByVHSUC3RSCPTMWIQAAOSo6CgX6KQQ6JmxCAEAyFHRUS7QSSHQM2MRAgCQo6KjXKCTQqBnxiIEACBHRUe5QCeFQM+MRQgAQI6KjnKBTgqBnhmLEACAHBUd5QKdFAI9MxYhAAA5KjrKBTopBHpmLEIAAHJUdJQLdFII9MxYhAAA5KjoKBfopBDombEIAQDIUdFRLtBJIdAzYxECAJCjoqNcoJNCoGfGIgQAIEdFR7lAJ4VAz4xFCABAjoqOcoFOCoGeGYsQAIAcFR3lAp0UAj0zFiEAADkqOsoFOikEemYsQgAAclR0lAt0Ugj0zFiEAADkqOgoF+ikEOiZsQgBAMhR0VEu0Ekh0DNjEQIAkKOio1ygk0KgZ8YiBAAgR0VHuUAnhUDPjEUIAECOio5ygU4KgZ4ZixAAgBwVHeUCnRQCPTMWIQAAOSo6ygU6KQR6ZixCAAByVHSUC3RSCPTMWIQAAOSo6CgX6KQQ6JmxCAEAyFHRUS7QSSHQM2MRAgCQo6KjXKCTQqBnxiIEACBHRUe5QCeFQM+MRQgAQI6KjnKBTgqBnhmLEACAHBUd5QKdFAI9MxYhAAA5KjrKBTopBHpmLEIAAHJUdJQLdFII9MxYhAAA5KjoKBfopBDombEIAQDIUdFRLtBJIdAzYxECAJCjoqNcoJNCoGfGIgQAIEdFR7lAJ4VAz4xFCABAjoqOcoFOCoGeGYsQAIAcFR3lAp0UAj0zFiEAADkqOsoFqtQpfAAAHy9JREFUOikEemYsQgAAclR0lAt0Ugj0zFiEAADkqOgoF+ikEOiZsQgBAMhR0VEu0Ekh0DNjEQIAkKOio1ygk0KgZ8YiBAAgR0VHuUAnhUDPjEUIAECOio5ygU4KgZ4ZixAAgBwVHeUCnRQCPTMWIQAAOSo6ygU6KQR6ZixCAAByVHSUC3RSCPTMWIQAAOSo6CgX6KQQ6JmxCAEAyFHRUS7QSSHQM2MRAgCQo6KjXKCTQqBnxiIEACBHRUe5QCeFQM+MRQgAQI6KjnKBTgqBnhmLEACAHBUd5QKdFAI9MxYhAAA5KjrKBTopBHpmLEIAAHJUdJQLdFII9MxYhAAA5KjoKBfopBDombEIAQDIUdFRLtBJIdAzYxECAJCjoqNcoJNCoGfGIgQAIEdFR7lAJ4VAz4xFCABAjoqOcoFOCoGeGYsQAIAcFR3lAp0UAj0zFiEAADkqOsoFOikEemYsQgAAclR0lAt0Ugj0zFiEAADkqOgoF+ikEOiZsQgBAMhR0VEu0Ekh0DNjEQIAkKOio1ygk0KgZ8YiBAAgR0VHuUAnhUDPjEUIAECOio5ygU4KgZ4ZixAAgBwVHeUCnRQCPTMWIQAAOSo6ygU6KQR6ZixCAAByVHSUC3RSCPTMWIQAAOSo6CgX6KQQ6JmxCAEAyFHRUS7QSSHQM2MRAgCQo6KjXKCTQqBnxiIEACBHRUe5QCeFQM+MRQgAQI6KjnKBTgqBnhmLEACAHBUd5QKdFAI9MxYhAAA5KjrKBTopBHpmLEIAAHJUdJQLdFII9MxYhAAAdIWiI7jMW1log/IR6JmxCAEA6ApFR3CZt7LQBuUj0DNjEQIA0BWKjuAyb2WhDcpHoGfGIgQAoCsUHcFl3spCG5SPQM+MRQgAQFcoOoLLvJWFNigfgZ6Zsi3Coi+OZbxQAgB8FhT92q7MW1mUrQ0Q6Nkp2yIs+uJYxgslAMBnQdGv7cq8lUXZ2gCBnp2yLcKiL45lvFACAHwWFP3arsxbWZStDRDo26RZs2bFbrvtFk1NTTFmzJh47LHHkm9btkVY9MWxjBdKAIDPgqJf25V5K4uytQECfZtzxx13RENDQ9xwww2xaNGimD59evTp0yfefPPNpNuXbREWfXEs44USAOCzoOjXdmXeyqJsbYBA3+aMHTs2pkyZ0mrfXnvtFTNmzEi6fdkWYdEXxzJvUK+i56xt21vXRT8eZX5syqLo58Vms22drSzK1gYI9G3Kxx9/HN27d4+777671f4zzzwzDjrooKq3WbduXaxataple+utt6JSqcSyZcta7S9q2+WsX9u2ga3oeWLO2Gw2m81m+yxtRb+eat6WLVsWlUolVq5c2RU5QwKBvg15++23o1KpxJNPPtlq/yWXXBJ77LFH1dvMnDkzKpWKzWaz2Ww2m81ms1Xdli1b1hU5QwKBvg1pDvSnnnqq1f6LL7449txzz6q3afsT9A8++CBee+21WLlyZWn+xq4sP823lX8zZ2z1buaMrd7NnLHVu5kztnq3Ms2ZlStXxrJly2L9+vVdkTMkEOjbkM68xb3MVq3yb16ojzlDvcwZ6mXOUC9zhnqZM3REoG9jxo4dG1OnTm21b/To0ckfElcmLk7Uy5yhXuYM9TJnqJc5Q73MGToi0Lcxzb9m7cYbb4xFixbFWWedFX369Ik33nij6KHVzcWJepkz1MucoV7mDPUyZ6iXOUNHBPo2aNasWTFixIhobGyMMWPGxNy5c4seUqesW7cuZs6cGevWrSt6KGwjzBnqZc5QL3OGepkz1MucoSMCHQAAAEpAoAMAAEAJCHQAAAAoAYEOAAAAJSDQAQAAoAQEOlvVrFmzYrfddoumpqYYM2ZMPPbYYx1+/Zw5c2LMmDHR1NQUI0eOjJ/85CddNFLKop45M3v27Bg/fnzssMMO0bdv3zjwwAPjwQcf7MLRUgb1XmeaPfHEE9G9e/fYb7/9tvIIKZt658y6devie9/7Xuy6667R2NgYu+++e9x4441dNFrKoN4588tf/jL23Xff6NWrVwwZMiT+4R/+Id57770uGi1Fmjt3bkyYMCGGDh0alUol7rnnnpq38fqXTQl0tprm39l+ww03xKJFi2L69OnRp0+fePPNN6t+/euvvx69e/eO6dOnx6JFi+KGG26IhoaGuOuuu7p45BSl3jkzffr0uPzyy+OZZ56JJUuWxD//8z9HQ0NDzJs3r4tHTlHqnTPNVq5cGbvvvnsceuihAj0znZkzEydOjC996Uvx0EMPxdKlS+Ppp5+OJ598sgtHTZHqnTOPP/54dOvWLa655pp4/fXX4/HHH4+//Mu/jKOPPrqLR04R7r///rjgggti9uzZSYHu9S9tCXS2mrFjx8aUKVNa7dtrr71ixowZVb/+u9/9buy1116t9p166qlx4IEHbrUxUi71zplq9t5777jooou29NAoqc7OmW9+85vxL//yLzFz5kyBnpl658wDDzwQ/fv3jxUrVnTF8CiheufMlVdeGbvvvnurfddee20MHz58q42RckoJdK9/aUugs1V8/PHH0b1797j77rtb7T/zzDPjoIMOqnqbr371q3HmmWe22nf33XdHjx494pNPPtlqY6UcOjNn2lq/fn3ssssu8eMf/3hrDJGS6eyc+bd/+7fYf//94//+7/8EemY6M2emTp0aBx98cJx//vkxbNiwGDVqVJx77rnx4YcfdsWQKVhn5syTTz4ZjY2Ncd9998WGDRvif/7nf+Kggw6KU089tSuGTImkBLrXv7Ql0Nkq3n777ahUKpu9BfCSSy6JPfbYo+ptRo0aFZdcckmrfU8++WRUKpV45513ttpYKYfOzJm2rrjiihg4cGD88Y9/3BpDpGQ6M2eWLFkSO+64Y7z88ssREQI9M52ZM4cddlg0NTXFEUccEU8//XTcd999MWLEiPj2t7/dFUOmYJ393nTnnXfG9ttvHz169IhKpRITJ04UWxlKCXSvf2lLoLNVNH9De+qpp1rtv/jii2PPPfeseptRo0bFpZde2mrfE088EZVKJZYvX77Vxko5dGbObOq2226L3r17x0MPPbS1hkjJ1DtnPv3009h///1bffiOQM9LZ64zhxxySPTs2TNWrlzZsm/27Nmx3Xbb+Sl6BjozZ1588cUYOnRoXHHFFfH888/Hgw8+GPvss0985zvf6YohUyKpge71L5sS6GwV3uJOvf6ct7jfcccd0atXr/jP//zPrTlESqbeOfPBBx9EpVKJ7t27t2zbbbddy75HHnmkq4ZOQTpznTnppJPic5/7XKt9ixYtikqlEkuWLNlqY6UcOjNnTjjhhDj22GNb7Xv88cf9RDRD3uJOZwh0tpqxY8fG1KlTW+0bPXp0hx8SN3r06Fb7pkyZ4kMyMlLvnInY+JPznj17Jv0aEz576pkz69evjwULFrTapk6dGnvuuWcsWLAg1q5d21XDpkD1Xmd+9rOfRa9evWLNmjUt++69997o1q2bn6Bnot45c8wxx8Rxxx3Xat9TTz0VlUol3n777a02Tson9UPivP5lUwKdrab515LceOONsWjRojjrrLOiT58+8cYbb0RExIwZM+LEE09s+frmXzNx9tlnx6JFi+LGG2/0ayYyU++cue2226JHjx4xa9asWL58ecu26VtR+Wyrd8605S3u+al3zqxZsyaGDx8exx57bLz44osxd+7cGDVqVEyePLmoU6CL1TtnbrrppujRo0dcf/318dprr8UTTzwR+++/f4wdO7aoU6ALrVmzJubPnx/z58+PSqUSV111VcyfP7/l1/J5/UstAp2tatasWTFixIhobGyMMWPGxNy5c1v+36RJk2LcuHGtvn7OnDnxxS9+MRobG2O33XZr9W9FyUM9c2bcuHFRqVQ22yZNmtT1A6cw9V5nNiXQ81TvnHnppZdi/Pjx0atXrxg+fHicc845fnqemXrnzLXXXht777139OrVK4YOHRp///d/H3/4wx+6eNQU4dFHH+3wtYnXv9Qi0AEAAKAEBDoAAACUgEAHAACAEhDoAAAAUAICHQAAAEpAoAMAAEAJCHQAAAAoAYEOAAAAJSDQAQAAoAQEOgAAAJSAQIdMjRs3LqZPn77Vj9H2a2r9eWvqyvvKTVfNp5x4PLaeIq9DW9Nn5Ty2tM48Lu+9914MHjw4li5d2iW368ijjz4aI0aMqPl1XfX8f+Mb34gf/vCHyfuB+gh0KLlJkyZFpVKJSqUSPXr0iJEjR8a5554ba9eu/bOO21VBtWLFili9enW7t6n1/7fk2LbmsbeEsr24vvTSS6NSqSSNKadA76pxluXxqGcebCtqXYe66n639HHKMmeqKXJsnbnvc889N77zne/UfV/t3W758uVx+umnx8iRI6OxsTGGDx8eEyZMiIcffrjmMbdWoHf2OXn++edj4MCBsWrVqqT9QH0EOpTcpEmT4vDDD4/ly5fHW2+9Fbfeemv06tUrpkyZUvXrP/7446TjFhVUtW4j0GsfN/U5/nM888wzsdtuu8W+++4r0NvIKdDrnQdl0tE6KeqxFehdO7bmOVDvfX/44YfxF3/xF/HUU0/VdX/t3W7p0qUxbNiw2HvvvePOO++Ml19+ORYuXBg//OEPY88996x53LIFekTEmDFj4vrrr0/eD6QT6FBykyZNiqOOOqrVvsmTJ8eQIUMiYuM32GnTpsXZZ58dgwYNioMOOigiItatWxdnnHFGDB48OJqamuLLX/5yPPPMMy3HaL7dtGnTon///jFw4MC44IILYsOGDS1f88ADD8SXv/zllv9/xBFHxKuvvlrXMep5i/um7xZo3i666KIYOHBgrFu3rtVjcMwxx8SJJ57Y7mPW9jhLly6NcePGxRlnnBHnnXdeDBgwIHbaaaeYOXNmq9vWOuf2jl1NrWOljHnTx7ntczxixIj40Y9+1Oo4++23X8s5bdiwIS6//PIYOXJk9OzZM/bdd9+48847q97/ptasWROjRo2Khx56KPkFXK25kDKWtveVModrPZ+rV6+Ob33rW9G7d+8YMmRIXHXVVZvdz5133hmf//zno2fPnjFw4MA4+OCD232HSnvPU62xVrN+/fq47LLL4nOf+1w0NjbGLrvsEhdffHHVxyNlLnV0HvWcY7N650HKGDs651qPR+o1re06Wbt2bZx44onRp0+fGDJkSPzgBz9Iui5tibm1qS01d2pdKzoad2euCSmPRa3jVhvztddeG/3794/169dHRMT8+fOjUqnEP/3TP7Xc7pRTTom/+7u/i4jaz3/zWKvNgbbPywMPPBD9+vWLW265peo5z549O3bYYYeqj0VH17n2bve1r30tdt5556pr7oMPPqg6hk1VC/SUed3RmuxoHqWs5QsvvDC++tWvbjbW9vYD6QQ6lFy1QD/jjDNi0KBBEbHxBcP2228f5513XixevDheeumliIg488wzY9iwYXH//ffHiy++GJMmTYoBAwbEihUrWt1u+vTpsXjx4vjlL38ZvXv3jp///Oct93PXXXfF7NmzY8mSJTF//vw48sgjY5999ml5QZVyjHoCfeXKlfFXf/VXcfLJJ8fy5ctj+fLlsXbt2ujfv3/8+te/brnNu+++G42NjfHb3/626mNW7TiffvppjBs3Lvr16xcXXnhhLFmyJG655ZbYbrvt4je/+U3yObd37GpqHStlzB09x7UC/Xvf+17stdde8eCDD8Zrr70WN910UzQ1NcWcOXOqjrfZSSedFGedddZmz09Has2FlLG0va+UOVzr+Zw8eXKMGDEiHn744ViwYEF8/etfj759+7bczzvvvBM9evSIq666KpYuXRovvPBCzJo1K9asWVP1PNt7nmqNtZrvfve7MWDAgLj55pvj1VdfjccffzxuuOGGqo9HrbnU0XnUe47N6p0HKfO9o3Ou9XikXtParpOpU6fG8OHD4ze/+U288MILMWHChJa5Wu2xbv7znzu32tpSc6fWtaKjcXfmmpDyWNQ6brUxr1y5Mrp16xb//d//HRERV199deywww5xwAEHtBx3jz32iJ/85CdJz39Hc2DT5/f222+Pvn37xr333tvuOU+fPj0OP/zwqo9FR9e5ardbsWJFbLfddnHppZe2e3+1VAv0lHnd0ZrsaB6lrOX7778/mpqaNvvL8/b2A+kEOpRc20B/+umnY9CgQXHcccdFxMYXDF/4whda3Wbt2rXR0NAQt956a8u+Tz75JIYNGxZXXHFFy+1Gjx7d6qfd559/fowePbrdsfzpT3+KSqUSCxYsSD5GPYFe7c8RG1+IfO1rX2v589VXXx277757q/ttq9pxxo0bF1/5ylda7TvggAPi/PPPTz7n9o6dotqxao25eX/b5zii40Bfu3Zt9OzZc7O3Wv7jP/5jHH/88e2O8fbbb4/Pf/7z8dFHH3U4pmpjbG8upI5l0/tKncMdPZ+rV6+OhoaGVj/JW7lyZfTu3bvlfp599tmoVCrxxhtv1DzHauNMHWtbq1evjqamplYBWut+NtV2LnV0Hp05x87Og47G2NE513o8UudD23WyZs2aaGxsjDvuuKNl34oVK6JXr141r0N/7tyqZkvMnWrHSRl3Z68JtR6LzqzvZmPGjIkf/OAHERFx9NFHxyWXXBKNjY2xevXqWL58eVQqlXjppZeSH6f2rpXN9z1r1qzo379/u3+52+yoo46q+u/Ia33Pq3a7p59+OiqVStx9990d3mez//3f/41dd901zj333JZ9bQM9dV63Ve17eMq6rva96/nnn696XWlvP5BOoEPJTZo0Kbp37x59+vSJpqam6NatW3z961+PP/7xjxGx8Rvs5MmTW92mvW+QRx99dHz7299uuV3zfze79957o0ePHi1/i/7qq6/G8ccfHyNHjoy+fftGnz59olKpxH333Zd8jC0R6PPmzYvu3bvHH/7wh4jYGKH/+q//2uHj1l6gn3baaa32TZw4sdU51Drn9o5dTcqxao25eX/b5zii40B/5plnolKpRJ8+fVptDQ0NMXbs2Kr3/9Zbb8WOO+4Yzz33XN3n2tFc+P3vf580lk3vK3UOd/R8Pvfcc1GpVOLNN99s9TVf/OIXW+7n008/jYMPPjj69u0bxx57bPz85z+P999/v+a5bvqYpIy1reYX7a+//nrS/dSaSx2dR73n2Nl5UGuMHZ1zrccjdT60XSftzYEvfOELNa9Df+7cqmZLzJ1qx0kZd2euCbWOGRHJx6025nPOOScmTJgQGzZsiEGDBsXChQtjzJgxcf/998dtt90WO+20U12PU3vXynHjxsXw4cOjoaEhnn766XbPtdmhhx662Tk3H6ej73nVbtd8/bvnnntq3m/Exncj/O3f/m2HgZ46r1O+h1ebRynfu5YsWRKVSiUWLVrU6rbt7QfSCXQouUmTJsX48ePjlVdeiTfeeCM++eSTVv+/2jfY9r55b/q3+ylxPXr06Dj00EPj4YcfjkWLFsXChQtbvdDoqkCP2PiTlksvvTSeffbZ6NatW7z11lsdPGrtB3rbfUcddVRMmjSp5c+1zrmjMbaVcqxa4+to/8iRI+Oqq65qtW/vvfeOmTNntrwonDNnTrzyyiuttvYeu3vuuScqlUp07969ZatUKrHddttF9+7d230rf/MYawV6rbFsep6pc7ij57P537S2Pd+2L2I3bNgQTzzxRHz/+9+PffbZJwYPHpwczqljbeuFF16oK9BT5lJH51HPOXZ2HtQaY0fnXOvx6Ox8aJ4DnQn0LTG32toSc6e98dUad2euCbWOGRGdWt/N/uM//iP69+8f8+fPj8GDB8eGDRvi7LPPjvPPPz9OOeWUlneKpT5OHT0uEyZMiGHDhsWUKVM6fPdVRMS3vvWtqu8qqPU9r9rt6nmL+5IlS+KYY46Jm266qcNAT53XKd/Dqz1eKdeb5uf93XffbXXb9vYD6QQ6lFy1f4O+qWrfYNeuXRuNjY2bvR1w5513jiuvvLLldm3fzj5jxoyWfe+9915UKpV47LHHWv7/448/vtk3946OUW18tf58yCGHxOmnn77ZeV5//fUxatSomDZtWhx66KHtPh4dHafWC82Uc+5ojJtKPVatMbc37oiIsWPHxnnnndfy51WrVkWvXr1i5syZLW8Z/vd///cOx7mp1atXx4IFC1pt+++/f5xwwgntvi1/0zG2NxdSx7LpeabO4Y6ez+a3Id91110t/3/VqlXRp0+fdiPq008/jZ133rnD3+Xb9nlKGWtbH330UfTq1SvpLe6dmUsdnUetc+zMPEgZY0fnXOvx6Ox8WLNmTTQ0NMSvfvWrln3vv//+Zm9FrzfQOzO3IrbM3Kl2nJRxd+aaUOuYEZF83Gpjbv536JMmTYpjjz02IjYG75e+9KXYY489YtasWRGR/jjV+ouLl19+OYYOHRrTpk3rcKxXXnll7LffflWP09H3vPZud/jhhyd9SNzEiRPj5ZdfrhnoKfM6ZU1We05Srze/+MUvYvjw4ZudT3v7gXQCHUquM4EesfHDaoYNGxYPPPBAqw/UaX5r67j//2E3Z599dixevDhuu+226NOnT/z0pz+NiI2fqDxo0KA44YQT4pVXXolHHnkkDjjggM0CvaNjVBtfrT+ffPLJccABB8TSpUvj3XffbflQmlWrVkXv3r03+3d37al2nFovNFPOuaMxbir1WLXGXO0xajZjxowYMmRIPPbYY7FgwYI4+uijY/vtt2/5kLgLLrggBg0a1PKhW/PmzYvrrrsubr755pqPX7PUdwvUmgspY2l7XylzuNY7IiZPnhwjR46M3/72t7Fw4cL4xje+EX379m358LPf//73cckll8R//dd/xZtvvhm//vWvo7GxMe6///52z7Xa81RrrNVceOGFMWDAgLjlllvi1Vdfjd/97nfxi1/8YrPHI2UudXQenTnHtmrNg9T53tE513o8OjMfIiKmTJkSu+66a8uHuU2cODHpQ+L+3LlVzZaaO/VcKzYdd2euCSmPRcpx2xvzmDFjonv37nHddddFxMbQbGhoiEqlEi+++GLL7VMep5R3FixevDiGDBnS4Xx+4YUXokePHps9B7Wuc+3d7vXXX48hQ4bE3nvvHXfddVcsWbIkFi1aFNdcc03stddeEbHxLyaaP8G+VqBH1J7XKWuy2nOSupYnTZpU9V0e7e0H0gl0KLnOBvpHH30UZ5xxRuywww7t/kqi0047LaZMmRL9+vWLAQMGxIwZM1q99e+hhx6K0aNHR1NTU+y7774xZ86czQK91jHqDfSXX345DjzwwOjVq9dmv8LsxBNPrPor16qpdpyUF5q1zrnWGDeVcqxaY672GDVbtWpVHHfccdGvX7/YZZdd4uabb97s16xdc801seeee0ZDQ0MMHjw4DjvssJg7d27Nx69ZPYHe0VxIGUvb+0qZw7Wez2q/Cmvs2LExY8aMiIhYtGhRHHbYYS2/ummPPfaIH//4xx2ea7XnqdZYq1m/fn1cfPHFMWLEiGhoaIhdd9211dtgNz2/WnOpo/PozDm2lTIPUuZ7R+dc6/HozHyI2PjTxhNOOCF69+4dO+20U1xxxRV1X5ci6p9b1WypuVPPtWLTcXfmmpDyWKQct70xn3vuuVGpVGLhwoUtX7vffvu1vOW9WcrjlBLoERvXxI477hjnnHNOu+d94IEHtvrL5ubj1PqeV+12ERt/08K0adNixIgR0djYGDvvvHNMnDgxHn300YjY+Beuw4cPjxEjRsSgQYOiX79+cdFFF0VE9UBPmde11mR7z0mt23300UfRr1+/+N3vftdqTO3tB+oj0IFtxvjx4+OMM84oehhsw5p/bd+mP5mFLcHc+my57777YvTo0a3eIZXyF1XVblevlJ+gF+m6666LQw45JHk/UB+BDpTeihUr4vbbb49u3brF4sWLix4O25B58+bFbbfdFq+++mo8++yzcdRRR0X//v19gBF/NnPrs+/qq6+u+UF3KberV9kD/Wc/+1nV78Xt7QfqI9CB0hsxYkT069evww9NgmrmzZsXY8aMiT59+sSAAQNi/Pjx8cILLxQ9LD4DzK38pAb6lla2QAe2LoEOAAAltXTp0vjRj35U9DCALiLQAQAAoAQEOgAAAJSAQAcAAIASEOgAAABQAgIdAAAASkCgAwAAQAkIdAAAACgBgQ4AAAAlINABAACgBAQ6AAAAlIBABwAAgBIQ6AAAAFACAh0AAABKQKADAABACQh0AAAAKAGBDgAAACUg0AEAAKAEBDoAAACUgEAHAACAEhDoAAAAUAICHQAAAEpAoAMAAEAJCHQAAAAoAYEOAAAAJSDQAQAAoAQEOgAAAJSAQAcAAIASEOgAAABQAgIdAAAASkCgAwAAQAkIdAAAACgBgQ4AAAAlINABAACgBAQ6AAAAlIBABwAAgBIQ6AAAAFACAh0AAABKQKADAABACQh0AAAAKAGBDgAAACUg0AEAAKAEBDoAAACUgEAHAACAEhDoAAAAUAICHQAAAEpAoAMAAEAJCHQAAAAoAYEOAAAAJSDQAQAAoAQEOgAAAJSAQAcAAIASEOgAAABQAgIdAAAASkCgAwAAQAkIdAAAACgBgQ4AAAAlINABAACgBAQ6AAAAlIBABwAAgBIQ6AAAAFACAh0AAABKQKADAABACQh0AAAAKAGBDgAAACXw/wBLJhMnk5qYpAAAAABJRU5ErkJggg==\" width=\"1000\">"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "fig, ax = plt.subplots(figsize=(10, 10))\n",
    "plt.hist(probs_4[:, 4], bins=20)\n",
    "ax.set(xlabel=r\"Probability that a true 4 belongs to class 4 according to the network (p(C$_4$|data))\",\n",
    "       ylabel=\"Number of occurences\",\n",
    "       title=\"\")\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "bff1f8db",
   "metadata": {},
   "source": [
    "### Contact us at the EuXFEL Data Analysis group at any time if you need help analysing your data!\n",
    "\n",
    "#### Data Analysis group: da@xfel.eu"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "0879d703",
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.6.13"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}