Untitled1.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "from scipy.optimize import least_squares"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "With wiki data cost equals =  0.21900131665664999 | With our approximate 0.11386136722898951\n"
     ]
    },
    {
     "data": {
      "application/javascript": [
       "/* Put everything inside the global mpl namespace */\n",
       "window.mpl = {};\n",
       "\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('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",
       "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 = $('<div/>');\n",
       "    this._root_extra_style(this.root)\n",
       "    this.root.attr('style', 'display: inline-block');\n",
       "\n",
       "    $(parent_element).append(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 (mpl.ratio != 1) {\n",
       "                fig.send_message(\"set_dpi_ratio\", {'dpi_ratio': mpl.ratio});\n",
       "            }\n",
       "            fig.send_message(\"refresh\", {});\n",
       "        }\n",
       "\n",
       "    this.imageObj.onload = function() {\n",
       "            if (fig.image_mode == 'full') {\n",
       "                // Full images could contain transparency (where diff images\n",
       "                // almost always do), so we need to clear the canvas so that\n",
       "                // there is no ghosting.\n",
       "                fig.context.clearRect(0, 0, fig.canvas.width, fig.canvas.height);\n",
       "            }\n",
       "            fig.context.drawImage(fig.imageObj, 0, 0);\n",
       "        };\n",
       "\n",
       "    this.imageObj.onunload = function() {\n",
       "        fig.ws.close();\n",
       "    }\n",
       "\n",
       "    this.ws.onmessage = this._make_on_message_function(this);\n",
       "\n",
       "    this.ondownload = ondownload;\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._init_header = function() {\n",
       "    var titlebar = $(\n",
       "        '<div class=\"ui-dialog-titlebar ui-widget-header ui-corner-all ' +\n",
       "        'ui-helper-clearfix\"/>');\n",
       "    var titletext = $(\n",
       "        '<div class=\"ui-dialog-title\" style=\"width: 100%; ' +\n",
       "        'text-align: center; padding: 3px;\"/>');\n",
       "    titlebar.append(titletext)\n",
       "    this.root.append(titlebar);\n",
       "    this.header = titletext[0];\n",
       "}\n",
       "\n",
       "\n",
       "\n",
       "mpl.figure.prototype._canvas_extra_style = function(canvas_div) {\n",
       "\n",
       "}\n",
       "\n",
       "\n",
       "mpl.figure.prototype._root_extra_style = function(canvas_div) {\n",
       "\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._init_canvas = function() {\n",
       "    var fig = this;\n",
       "\n",
       "    var canvas_div = $('<div/>');\n",
       "\n",
       "    canvas_div.attr('style', 'position: relative; clear: both; outline: 0');\n",
       "\n",
       "    function canvas_keyboard_event(event) {\n",
       "        return fig.key_event(event, event['data']);\n",
       "    }\n",
       "\n",
       "    canvas_div.keydown('key_press', canvas_keyboard_event);\n",
       "    canvas_div.keyup('key_release', canvas_keyboard_event);\n",
       "    this.canvas_div = canvas_div\n",
       "    this._canvas_extra_style(canvas_div)\n",
       "    this.root.append(canvas_div);\n",
       "\n",
       "    var canvas = $('<canvas/>');\n",
       "    canvas.addClass('mpl-canvas');\n",
       "    canvas.attr('style', \"left: 0; top: 0; z-index: 0; outline: 0\")\n",
       "\n",
       "    this.canvas = canvas[0];\n",
       "    this.context = canvas[0].getContext(\"2d\");\n",
       "\n",
       "    var backingStore = this.context.backingStorePixelRatio ||\n",
       "\tthis.context.webkitBackingStorePixelRatio ||\n",
       "\tthis.context.mozBackingStorePixelRatio ||\n",
       "\tthis.context.msBackingStorePixelRatio ||\n",
       "\tthis.context.oBackingStorePixelRatio ||\n",
       "\tthis.context.backingStorePixelRatio || 1;\n",
       "\n",
       "    mpl.ratio = (window.devicePixelRatio || 1) / backingStore;\n",
       "\n",
       "    var rubberband = $('<canvas/>');\n",
       "    rubberband.attr('style', \"position: absolute; left: 0; top: 0; z-index: 1;\")\n",
       "\n",
       "    var pass_mouse_events = true;\n",
       "\n",
       "    canvas_div.resizable({\n",
       "        start: function(event, ui) {\n",
       "            pass_mouse_events = false;\n",
       "        },\n",
       "        resize: function(event, ui) {\n",
       "            fig.request_resize(ui.size.width, ui.size.height);\n",
       "        },\n",
       "        stop: function(event, ui) {\n",
       "            pass_mouse_events = true;\n",
       "            fig.request_resize(ui.size.width, ui.size.height);\n",
       "        },\n",
       "    });\n",
       "\n",
       "    function mouse_event_fn(event) {\n",
       "        if (pass_mouse_events)\n",
       "            return fig.mouse_event(event, event['data']);\n",
       "    }\n",
       "\n",
       "    rubberband.mousedown('button_press', mouse_event_fn);\n",
       "    rubberband.mouseup('button_release', mouse_event_fn);\n",
       "    // Throttle sequential mouse events to 1 every 20ms.\n",
       "    rubberband.mousemove('motion_notify', mouse_event_fn);\n",
       "\n",
       "    rubberband.mouseenter('figure_enter', mouse_event_fn);\n",
       "    rubberband.mouseleave('figure_leave', mouse_event_fn);\n",
       "\n",
       "    canvas_div.on(\"wheel\", function (event) {\n",
       "        event = event.originalEvent;\n",
       "        event['data'] = 'scroll'\n",
       "        if (event.deltaY < 0) {\n",
       "            event.step = 1;\n",
       "        } else {\n",
       "            event.step = -1;\n",
       "        }\n",
       "        mouse_event_fn(event);\n",
       "    });\n",
       "\n",
       "    canvas_div.append(canvas);\n",
       "    canvas_div.append(rubberband);\n",
       "\n",
       "    this.rubberband = rubberband;\n",
       "    this.rubberband_canvas = rubberband[0];\n",
       "    this.rubberband_context = rubberband[0].getContext(\"2d\");\n",
       "    this.rubberband_context.strokeStyle = \"#000000\";\n",
       "\n",
       "    this._resize_canvas = function(width, height) {\n",
       "        // Keep the size of the canvas, canvas container, and rubber band\n",
       "        // canvas in synch.\n",
       "        canvas_div.css('width', width)\n",
       "        canvas_div.css('height', height)\n",
       "\n",
       "        canvas.attr('width', width * mpl.ratio);\n",
       "        canvas.attr('height', height * mpl.ratio);\n",
       "        canvas.attr('style', 'width: ' + width + 'px; height: ' + height + 'px;');\n",
       "\n",
       "        rubberband.attr('width', width);\n",
       "        rubberband.attr('height', height);\n",
       "    }\n",
       "\n",
       "    // Set the figure to an initial 600x600px, this will subsequently be updated\n",
       "    // upon first draw.\n",
       "    this._resize_canvas(600, 600);\n",
       "\n",
       "    // Disable right mouse context menu.\n",
       "    $(this.rubberband_canvas).bind(\"contextmenu\",function(e){\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 nav_element = $('<div/>')\n",
       "    nav_element.attr('style', 'width: 100%');\n",
       "    this.root.append(nav_element);\n",
       "\n",
       "    // Define a callback function for later on.\n",
       "    function toolbar_event(event) {\n",
       "        return fig.toolbar_button_onclick(event['data']);\n",
       "    }\n",
       "    function toolbar_mouse_event(event) {\n",
       "        return fig.toolbar_button_onmouseover(event['data']);\n",
       "    }\n",
       "\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",
       "            // put a spacer in here.\n",
       "            continue;\n",
       "        }\n",
       "        var button = $('<button/>');\n",
       "        button.addClass('ui-button ui-widget ui-state-default ui-corner-all ' +\n",
       "                        'ui-button-icon-only');\n",
       "        button.attr('role', 'button');\n",
       "        button.attr('aria-disabled', 'false');\n",
       "        button.click(method_name, toolbar_event);\n",
       "        button.mouseover(tooltip, toolbar_mouse_event);\n",
       "\n",
       "        var icon_img = $('<span/>');\n",
       "        icon_img.addClass('ui-button-icon-primary ui-icon');\n",
       "        icon_img.addClass(image);\n",
       "        icon_img.addClass('ui-corner-all');\n",
       "\n",
       "        var tooltip_span = $('<span/>');\n",
       "        tooltip_span.addClass('ui-button-text');\n",
       "        tooltip_span.html(tooltip);\n",
       "\n",
       "        button.append(icon_img);\n",
       "        button.append(tooltip_span);\n",
       "\n",
       "        nav_element.append(button);\n",
       "    }\n",
       "\n",
       "    var fmt_picker_span = $('<span/>');\n",
       "\n",
       "    var fmt_picker = $('<select/>');\n",
       "    fmt_picker.addClass('mpl-toolbar-option ui-widget ui-widget-content');\n",
       "    fmt_picker_span.append(fmt_picker);\n",
       "    nav_element.append(fmt_picker_span);\n",
       "    this.format_dropdown = fmt_picker[0];\n",
       "\n",
       "    for (var ind in mpl.extensions) {\n",
       "        var fmt = mpl.extensions[ind];\n",
       "        var option = $(\n",
       "            '<option/>', {selected: fmt === mpl.default_extension}).html(fmt);\n",
       "        fmt_picker.append(option)\n",
       "    }\n",
       "\n",
       "    // Add hover states to the ui-buttons\n",
       "    $( \".ui-button\" ).hover(\n",
       "        function() { $(this).addClass(\"ui-state-hover\");},\n",
       "        function() { $(this).removeClass(\"ui-state-hover\");}\n",
       "    );\n",
       "\n",
       "    var status_bar = $('<span class=\"mpl-message\"/>');\n",
       "    nav_element.append(status_bar);\n",
       "    this.message = status_bar[0];\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",
       "\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",
       "\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]);\n",
       "        fig.send_message(\"refresh\", {});\n",
       "    };\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_rubberband = function(fig, msg) {\n",
       "    var x0 = msg['x0'] / mpl.ratio;\n",
       "    var y0 = (fig.canvas.height - msg['y0']) / mpl.ratio;\n",
       "    var x1 = msg['x1'] / mpl.ratio;\n",
       "    var y1 = (fig.canvas.height - msg['y1']) / mpl.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, 0, fig.canvas.width, fig.canvas.height);\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",
       "    {\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.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",
       "            fig.imageObj.src = (window.URL || window.webkitURL).createObjectURL(\n",
       "                evt.data);\n",
       "            fig.updated_canvas_event();\n",
       "            fig.waiting = false;\n",
       "            return;\n",
       "        }\n",
       "        else if (typeof evt.data === 'string' && evt.data.slice(0, 21) == \"data:image/png;base64\") {\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(\"No handler for the '\" + msg_type + \"' message type: \", msg);\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(\"Exception inside the 'handler_\" + msg_type + \"' callback:\", e, e.stack, msg);\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",
       "    if (e.target)\n",
       "        targ = e.target;\n",
       "    else if (e.srcElement)\n",
       "        targ = e.srcElement;\n",
       "    if (targ.nodeType == 3) // defeat Safari bug\n",
       "        targ = targ.parentNode;\n",
       "\n",
       "    // jQuery normalizes the pageX and pageY\n",
       "    // pageX,Y are the mouse positions relative to the document\n",
       "    // offset() returns the position of the element relative to the document\n",
       "    var x = e.pageX - $(targ).offset().left;\n",
       "    var y = e.pageY - $(targ).offset().top;\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",
       "    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",
       "    {\n",
       "        this.canvas.focus();\n",
       "        this.canvas_div.focus();\n",
       "    }\n",
       "\n",
       "    var x = canvas_pos.x * mpl.ratio;\n",
       "    var y = canvas_pos.y * mpl.ratio;\n",
       "\n",
       "    this.send_message(name, {x: x, y: y, button: event.button,\n",
       "                             step: event.step,\n",
       "                             guiEvent: simpleKeys(event)});\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",
       "\n",
       "    // Prevent repeat events\n",
       "    if (name == 'key_press')\n",
       "    {\n",
       "        if (event.which === this._key)\n",
       "            return;\n",
       "        else\n",
       "            this._key = event.which;\n",
       "    }\n",
       "    if (name == 'key_release')\n",
       "        this._key = null;\n",
       "\n",
       "    var value = '';\n",
       "    if (event.ctrlKey && event.which != 17)\n",
       "        value += \"ctrl+\";\n",
       "    if (event.altKey && event.which != 18)\n",
       "        value += \"alt+\";\n",
       "    if (event.shiftKey && event.which != 16)\n",
       "        value += \"shift+\";\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,\n",
       "                             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",
       "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\", \"Pan axes with left mouse, zoom with right\", \"fa fa-arrows icon-move\", \"pan\"], [\"Zoom\", \"Zoom to rectangle\", \"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\";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 overriden (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 = $(\"#\" + id);\n",
       "    var ws_proxy = comm_websocket_adapter(comm)\n",
       "\n",
       "    function ondownload(figure, format) {\n",
       "        window.open(figure.imageObj.src);\n",
       "    }\n",
       "\n",
       "    var fig = new mpl.figure(id, ws_proxy,\n",
       "                           ondownload,\n",
       "                           element.get(0));\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.get(0);\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",
       "\n",
       "    var output_index = fig.cell_info[2]\n",
       "    var cell = fig.cell_info[0];\n",
       "\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_close = function(fig, msg) {\n",
       "    var width = fig.canvas.width/mpl.ratio\n",
       "    fig.root.unbind('remove')\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).html('<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/mpl.ratio\n",
       "    var dataURL = this.canvas.toDataURL();\n",
       "    this.cell_info[1]['text/html'] = '<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 () { fig.push_to_output() }, 1000);\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._init_toolbar = function() {\n",
       "    var fig = this;\n",
       "\n",
       "    var nav_element = $('<div/>')\n",
       "    nav_element.attr('style', 'width: 100%');\n",
       "    this.root.append(nav_element);\n",
       "\n",
       "    // Define a callback function for later on.\n",
       "    function toolbar_event(event) {\n",
       "        return fig.toolbar_button_onclick(event['data']);\n",
       "    }\n",
       "    function toolbar_mouse_event(event) {\n",
       "        return fig.toolbar_button_onmouseover(event['data']);\n",
       "    }\n",
       "\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) { continue; };\n",
       "\n",
       "        var button = $('<button class=\"btn btn-default\" href=\"#\" title=\"' + name + '\"><i class=\"fa ' + image + ' fa-lg\"></i></button>');\n",
       "        button.click(method_name, toolbar_event);\n",
       "        button.mouseover(tooltip, toolbar_mouse_event);\n",
       "        nav_element.append(button);\n",
       "    }\n",
       "\n",
       "    // Add the status bar.\n",
       "    var status_bar = $('<span class=\"mpl-message\" style=\"text-align:right; float: right;\"/>');\n",
       "    nav_element.append(status_bar);\n",
       "    this.message = status_bar[0];\n",
       "\n",
       "    // Add the close button to the window.\n",
       "    var buttongrp = $('<div class=\"btn-group inline pull-right\"></div>');\n",
       "    var button = $('<button class=\"btn btn-mini btn-primary\" href=\"#\" title=\"Stop Interaction\"><i class=\"fa fa-power-off icon-remove icon-large\"></i></button>');\n",
       "    button.click(function (evt) { fig.handle_close(fig, {}); } );\n",
       "    button.mouseover('Stop Interaction', toolbar_mouse_event);\n",
       "    buttongrp.append(button);\n",
       "    var titlebar = this.root.find($('.ui-dialog-titlebar'));\n",
       "    titlebar.prepend(buttongrp);\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._root_extra_style = function(el){\n",
       "    var fig = this\n",
       "    el.on(\"remove\", function(){\n",
       "\tfig.close_ws(fig, {});\n",
       "    });\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._canvas_extra_style = function(el){\n",
       "    // this is important to make the div 'focusable\n",
       "    el.attr('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",
       "    }\n",
       "    else {\n",
       "        // location in version 2\n",
       "        IPython.keyboard_manager.register_events(el);\n",
       "    }\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",
       "    // Check for shift+enter\n",
       "    if (event.shiftKey && event.which == 13) {\n",
       "        this.canvas_div.blur();\n",
       "        event.shiftKey = false;\n",
       "        // Send a \"J\" for go to next cell\n",
       "        event.which = 74;\n",
       "        event.keyCode = 74;\n",
       "        manager.command_mode();\n",
       "        manager.handle_keydown(event);\n",
       "    }\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_save = function(fig, msg) {\n",
       "    fig.ondownload(fig, null);\n",
       "}\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('matplotlib', mpl.mpl_figure_comm);\n",
       "}\n"
      ],
      "text/plain": [
       "<IPython.core.display.Javascript object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/html": [
       "<img src=\"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAp8AAAH3CAYAAAD9pzSiAAAAAXNSR0IArs4c6QAAAARnQU1BAACxjwv8YQUAAAU0SURBVHhe7cEBDQAAAMKg909tDjcgAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAeFQDnMcAAdG9ztIAAAAASUVORK5CYII=\" width=\"639.0476480687299\">"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "application/javascript": [
       "/* Put everything inside the global mpl namespace */\n",
       "window.mpl = {};\n",
       "\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('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",
       "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 = $('<div/>');\n",
       "    this._root_extra_style(this.root)\n",
       "    this.root.attr('style', 'display: inline-block');\n",
       "\n",
       "    $(parent_element).append(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 (mpl.ratio != 1) {\n",
       "                fig.send_message(\"set_dpi_ratio\", {'dpi_ratio': mpl.ratio});\n",
       "            }\n",
       "            fig.send_message(\"refresh\", {});\n",
       "        }\n",
       "\n",
       "    this.imageObj.onload = function() {\n",
       "            if (fig.image_mode == 'full') {\n",
       "                // Full images could contain transparency (where diff images\n",
       "                // almost always do), so we need to clear the canvas so that\n",
       "                // there is no ghosting.\n",
       "                fig.context.clearRect(0, 0, fig.canvas.width, fig.canvas.height);\n",
       "            }\n",
       "            fig.context.drawImage(fig.imageObj, 0, 0);\n",
       "        };\n",
       "\n",
       "    this.imageObj.onunload = function() {\n",
       "        fig.ws.close();\n",
       "    }\n",
       "\n",
       "    this.ws.onmessage = this._make_on_message_function(this);\n",
       "\n",
       "    this.ondownload = ondownload;\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._init_header = function() {\n",
       "    var titlebar = $(\n",
       "        '<div class=\"ui-dialog-titlebar ui-widget-header ui-corner-all ' +\n",
       "        'ui-helper-clearfix\"/>');\n",
       "    var titletext = $(\n",
       "        '<div class=\"ui-dialog-title\" style=\"width: 100%; ' +\n",
       "        'text-align: center; padding: 3px;\"/>');\n",
       "    titlebar.append(titletext)\n",
       "    this.root.append(titlebar);\n",
       "    this.header = titletext[0];\n",
       "}\n",
       "\n",
       "\n",
       "\n",
       "mpl.figure.prototype._canvas_extra_style = function(canvas_div) {\n",
       "\n",
       "}\n",
       "\n",
       "\n",
       "mpl.figure.prototype._root_extra_style = function(canvas_div) {\n",
       "\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._init_canvas = function() {\n",
       "    var fig = this;\n",
       "\n",
       "    var canvas_div = $('<div/>');\n",
       "\n",
       "    canvas_div.attr('style', 'position: relative; clear: both; outline: 0');\n",
       "\n",
       "    function canvas_keyboard_event(event) {\n",
       "        return fig.key_event(event, event['data']);\n",
       "    }\n",
       "\n",
       "    canvas_div.keydown('key_press', canvas_keyboard_event);\n",
       "    canvas_div.keyup('key_release', canvas_keyboard_event);\n",
       "    this.canvas_div = canvas_div\n",
       "    this._canvas_extra_style(canvas_div)\n",
       "    this.root.append(canvas_div);\n",
       "\n",
       "    var canvas = $('<canvas/>');\n",
       "    canvas.addClass('mpl-canvas');\n",
       "    canvas.attr('style', \"left: 0; top: 0; z-index: 0; outline: 0\")\n",
       "\n",
       "    this.canvas = canvas[0];\n",
       "    this.context = canvas[0].getContext(\"2d\");\n",
       "\n",
       "    var backingStore = this.context.backingStorePixelRatio ||\n",
       "\tthis.context.webkitBackingStorePixelRatio ||\n",
       "\tthis.context.mozBackingStorePixelRatio ||\n",
       "\tthis.context.msBackingStorePixelRatio ||\n",
       "\tthis.context.oBackingStorePixelRatio ||\n",
       "\tthis.context.backingStorePixelRatio || 1;\n",
       "\n",
       "    mpl.ratio = (window.devicePixelRatio || 1) / backingStore;\n",
       "\n",
       "    var rubberband = $('<canvas/>');\n",
       "    rubberband.attr('style', \"position: absolute; left: 0; top: 0; z-index: 1;\")\n",
       "\n",
       "    var pass_mouse_events = true;\n",
       "\n",
       "    canvas_div.resizable({\n",
       "        start: function(event, ui) {\n",
       "            pass_mouse_events = false;\n",
       "        },\n",
       "        resize: function(event, ui) {\n",
       "            fig.request_resize(ui.size.width, ui.size.height);\n",
       "        },\n",
       "        stop: function(event, ui) {\n",
       "            pass_mouse_events = true;\n",
       "            fig.request_resize(ui.size.width, ui.size.height);\n",
       "        },\n",
       "    });\n",
       "\n",
       "    function mouse_event_fn(event) {\n",
       "        if (pass_mouse_events)\n",
       "            return fig.mouse_event(event, event['data']);\n",
       "    }\n",
       "\n",
       "    rubberband.mousedown('button_press', mouse_event_fn);\n",
       "    rubberband.mouseup('button_release', mouse_event_fn);\n",
       "    // Throttle sequential mouse events to 1 every 20ms.\n",
       "    rubberband.mousemove('motion_notify', mouse_event_fn);\n",
       "\n",
       "    rubberband.mouseenter('figure_enter', mouse_event_fn);\n",
       "    rubberband.mouseleave('figure_leave', mouse_event_fn);\n",
       "\n",
       "    canvas_div.on(\"wheel\", function (event) {\n",
       "        event = event.originalEvent;\n",
       "        event['data'] = 'scroll'\n",
       "        if (event.deltaY < 0) {\n",
       "            event.step = 1;\n",
       "        } else {\n",
       "            event.step = -1;\n",
       "        }\n",
       "        mouse_event_fn(event);\n",
       "    });\n",
       "\n",
       "    canvas_div.append(canvas);\n",
       "    canvas_div.append(rubberband);\n",
       "\n",
       "    this.rubberband = rubberband;\n",
       "    this.rubberband_canvas = rubberband[0];\n",
       "    this.rubberband_context = rubberband[0].getContext(\"2d\");\n",
       "    this.rubberband_context.strokeStyle = \"#000000\";\n",
       "\n",
       "    this._resize_canvas = function(width, height) {\n",
       "        // Keep the size of the canvas, canvas container, and rubber band\n",
       "        // canvas in synch.\n",
       "        canvas_div.css('width', width)\n",
       "        canvas_div.css('height', height)\n",
       "\n",
       "        canvas.attr('width', width * mpl.ratio);\n",
       "        canvas.attr('height', height * mpl.ratio);\n",
       "        canvas.attr('style', 'width: ' + width + 'px; height: ' + height + 'px;');\n",
       "\n",
       "        rubberband.attr('width', width);\n",
       "        rubberband.attr('height', height);\n",
       "    }\n",
       "\n",
       "    // Set the figure to an initial 600x600px, this will subsequently be updated\n",
       "    // upon first draw.\n",
       "    this._resize_canvas(600, 600);\n",
       "\n",
       "    // Disable right mouse context menu.\n",
       "    $(this.rubberband_canvas).bind(\"contextmenu\",function(e){\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 nav_element = $('<div/>')\n",
       "    nav_element.attr('style', 'width: 100%');\n",
       "    this.root.append(nav_element);\n",
       "\n",
       "    // Define a callback function for later on.\n",
       "    function toolbar_event(event) {\n",
       "        return fig.toolbar_button_onclick(event['data']);\n",
       "    }\n",
       "    function toolbar_mouse_event(event) {\n",
       "        return fig.toolbar_button_onmouseover(event['data']);\n",
       "    }\n",
       "\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",
       "            // put a spacer in here.\n",
       "            continue;\n",
       "        }\n",
       "        var button = $('<button/>');\n",
       "        button.addClass('ui-button ui-widget ui-state-default ui-corner-all ' +\n",
       "                        'ui-button-icon-only');\n",
       "        button.attr('role', 'button');\n",
       "        button.attr('aria-disabled', 'false');\n",
       "        button.click(method_name, toolbar_event);\n",
       "        button.mouseover(tooltip, toolbar_mouse_event);\n",
       "\n",
       "        var icon_img = $('<span/>');\n",
       "        icon_img.addClass('ui-button-icon-primary ui-icon');\n",
       "        icon_img.addClass(image);\n",
       "        icon_img.addClass('ui-corner-all');\n",
       "\n",
       "        var tooltip_span = $('<span/>');\n",
       "        tooltip_span.addClass('ui-button-text');\n",
       "        tooltip_span.html(tooltip);\n",
       "\n",
       "        button.append(icon_img);\n",
       "        button.append(tooltip_span);\n",
       "\n",
       "        nav_element.append(button);\n",
       "    }\n",
       "\n",
       "    var fmt_picker_span = $('<span/>');\n",
       "\n",
       "    var fmt_picker = $('<select/>');\n",
       "    fmt_picker.addClass('mpl-toolbar-option ui-widget ui-widget-content');\n",
       "    fmt_picker_span.append(fmt_picker);\n",
       "    nav_element.append(fmt_picker_span);\n",
       "    this.format_dropdown = fmt_picker[0];\n",
       "\n",
       "    for (var ind in mpl.extensions) {\n",
       "        var fmt = mpl.extensions[ind];\n",
       "        var option = $(\n",
       "            '<option/>', {selected: fmt === mpl.default_extension}).html(fmt);\n",
       "        fmt_picker.append(option)\n",
       "    }\n",
       "\n",
       "    // Add hover states to the ui-buttons\n",
       "    $( \".ui-button\" ).hover(\n",
       "        function() { $(this).addClass(\"ui-state-hover\");},\n",
       "        function() { $(this).removeClass(\"ui-state-hover\");}\n",
       "    );\n",
       "\n",
       "    var status_bar = $('<span class=\"mpl-message\"/>');\n",
       "    nav_element.append(status_bar);\n",
       "    this.message = status_bar[0];\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",
       "\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",
       "\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]);\n",
       "        fig.send_message(\"refresh\", {});\n",
       "    };\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_rubberband = function(fig, msg) {\n",
       "    var x0 = msg['x0'] / mpl.ratio;\n",
       "    var y0 = (fig.canvas.height - msg['y0']) / mpl.ratio;\n",
       "    var x1 = msg['x1'] / mpl.ratio;\n",
       "    var y1 = (fig.canvas.height - msg['y1']) / mpl.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, 0, fig.canvas.width, fig.canvas.height);\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",
       "    {\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.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",
       "            fig.imageObj.src = (window.URL || window.webkitURL).createObjectURL(\n",
       "                evt.data);\n",
       "            fig.updated_canvas_event();\n",
       "            fig.waiting = false;\n",
       "            return;\n",
       "        }\n",
       "        else if (typeof evt.data === 'string' && evt.data.slice(0, 21) == \"data:image/png;base64\") {\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(\"No handler for the '\" + msg_type + \"' message type: \", msg);\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(\"Exception inside the 'handler_\" + msg_type + \"' callback:\", e, e.stack, msg);\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",
       "    if (e.target)\n",
       "        targ = e.target;\n",
       "    else if (e.srcElement)\n",
       "        targ = e.srcElement;\n",
       "    if (targ.nodeType == 3) // defeat Safari bug\n",
       "        targ = targ.parentNode;\n",
       "\n",
       "    // jQuery normalizes the pageX and pageY\n",
       "    // pageX,Y are the mouse positions relative to the document\n",
       "    // offset() returns the position of the element relative to the document\n",
       "    var x = e.pageX - $(targ).offset().left;\n",
       "    var y = e.pageY - $(targ).offset().top;\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",
       "    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",
       "    {\n",
       "        this.canvas.focus();\n",
       "        this.canvas_div.focus();\n",
       "    }\n",
       "\n",
       "    var x = canvas_pos.x * mpl.ratio;\n",
       "    var y = canvas_pos.y * mpl.ratio;\n",
       "\n",
       "    this.send_message(name, {x: x, y: y, button: event.button,\n",
       "                             step: event.step,\n",
       "                             guiEvent: simpleKeys(event)});\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",
       "\n",
       "    // Prevent repeat events\n",
       "    if (name == 'key_press')\n",
       "    {\n",
       "        if (event.which === this._key)\n",
       "            return;\n",
       "        else\n",
       "            this._key = event.which;\n",
       "    }\n",
       "    if (name == 'key_release')\n",
       "        this._key = null;\n",
       "\n",
       "    var value = '';\n",
       "    if (event.ctrlKey && event.which != 17)\n",
       "        value += \"ctrl+\";\n",
       "    if (event.altKey && event.which != 18)\n",
       "        value += \"alt+\";\n",
       "    if (event.shiftKey && event.which != 16)\n",
       "        value += \"shift+\";\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,\n",
       "                             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",
       "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\", \"Pan axes with left mouse, zoom with right\", \"fa fa-arrows icon-move\", \"pan\"], [\"Zoom\", \"Zoom to rectangle\", \"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\";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 overriden (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 = $(\"#\" + id);\n",
       "    var ws_proxy = comm_websocket_adapter(comm)\n",
       "\n",
       "    function ondownload(figure, format) {\n",
       "        window.open(figure.imageObj.src);\n",
       "    }\n",
       "\n",
       "    var fig = new mpl.figure(id, ws_proxy,\n",
       "                           ondownload,\n",
       "                           element.get(0));\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.get(0);\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",
       "\n",
       "    var output_index = fig.cell_info[2]\n",
       "    var cell = fig.cell_info[0];\n",
       "\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_close = function(fig, msg) {\n",
       "    var width = fig.canvas.width/mpl.ratio\n",
       "    fig.root.unbind('remove')\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).html('<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/mpl.ratio\n",
       "    var dataURL = this.canvas.toDataURL();\n",
       "    this.cell_info[1]['text/html'] = '<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 () { fig.push_to_output() }, 1000);\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._init_toolbar = function() {\n",
       "    var fig = this;\n",
       "\n",
       "    var nav_element = $('<div/>')\n",
       "    nav_element.attr('style', 'width: 100%');\n",
       "    this.root.append(nav_element);\n",
       "\n",
       "    // Define a callback function for later on.\n",
       "    function toolbar_event(event) {\n",
       "        return fig.toolbar_button_onclick(event['data']);\n",
       "    }\n",
       "    function toolbar_mouse_event(event) {\n",
       "        return fig.toolbar_button_onmouseover(event['data']);\n",
       "    }\n",
       "\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) { continue; };\n",
       "\n",
       "        var button = $('<button class=\"btn btn-default\" href=\"#\" title=\"' + name + '\"><i class=\"fa ' + image + ' fa-lg\"></i></button>');\n",
       "        button.click(method_name, toolbar_event);\n",
       "        button.mouseover(tooltip, toolbar_mouse_event);\n",
       "        nav_element.append(button);\n",
       "    }\n",
       "\n",
       "    // Add the status bar.\n",
       "    var status_bar = $('<span class=\"mpl-message\" style=\"text-align:right; float: right;\"/>');\n",
       "    nav_element.append(status_bar);\n",
       "    this.message = status_bar[0];\n",
       "\n",
       "    // Add the close button to the window.\n",
       "    var buttongrp = $('<div class=\"btn-group inline pull-right\"></div>');\n",
       "    var button = $('<button class=\"btn btn-mini btn-primary\" href=\"#\" title=\"Stop Interaction\"><i class=\"fa fa-power-off icon-remove icon-large\"></i></button>');\n",
       "    button.click(function (evt) { fig.handle_close(fig, {}); } );\n",
       "    button.mouseover('Stop Interaction', toolbar_mouse_event);\n",
       "    buttongrp.append(button);\n",
       "    var titlebar = this.root.find($('.ui-dialog-titlebar'));\n",
       "    titlebar.prepend(buttongrp);\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._root_extra_style = function(el){\n",
       "    var fig = this\n",
       "    el.on(\"remove\", function(){\n",
       "\tfig.close_ws(fig, {});\n",
       "    });\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._canvas_extra_style = function(el){\n",
       "    // this is important to make the div 'focusable\n",
       "    el.attr('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",
       "    }\n",
       "    else {\n",
       "        // location in version 2\n",
       "        IPython.keyboard_manager.register_events(el);\n",
       "    }\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",
       "    // Check for shift+enter\n",
       "    if (event.shiftKey && event.which == 13) {\n",
       "        this.canvas_div.blur();\n",
       "        event.shiftKey = false;\n",
       "        // Send a \"J\" for go to next cell\n",
       "        event.which = 74;\n",
       "        event.keyCode = 74;\n",
       "        manager.command_mode();\n",
       "        manager.handle_keydown(event);\n",
       "    }\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_save = function(fig, msg) {\n",
       "    fig.ondownload(fig, null);\n",
       "}\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('matplotlib', mpl.mpl_figure_comm);\n",
       "}\n"
      ],
      "text/plain": [
       "<IPython.core.display.Javascript object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/html": [
       "<img src=\"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAt4AAALeCAYAAACdlBZfAAAAAXNSR0IArs4c6QAAAARnQU1BAACxjwv8YQUAAH6oSURBVHhe7d0HnFxV+cbxd7MtDTa9dyCQIFUEBAEJgqGKoqK0AAKCgAICAgoJRSmhg4LwV4g0QREpwQgSwIJIh0BCDQmb3nuy/f8+t5DJZsvs7uzulN+Xz8OcuXcyO3vv7O57z5x7rgEAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAABA4KueGs+JwT3kgmEe7fMJwb32pddxb9hsFO9VAHXqEN0CSE9ben7uedWz0lPmmeW5z7OvB21LhZQKqmrPF7WgltM9Wq/Cq7aungs8L3mWeSo8iz1TPed5SjzJ2tbzV4/+/SrP656fepoiLmqvCe7VL/6elRO0oA7ne+LHjNWCHNDNowOCrwX3ACAJFN5A+hrtmea50jPfM95zhucPnp09L3r28uSKf3o6eXTQ0d7yPFeFzaRs53nHc61nhedXntM8Knp1/2rP055kqED/h2d/z288P/Oo8B7naU0bPPV9DRXkWp/N9N47NWwGVHjrZ5LCGwCADKfe0U88KmYO0oJaVPid7Nk9uJfdOnvSpZMg7v1Voavbr3gS1dXj3cXzsWe9p659KQM8+mQjGd/w6Gv8ILi3UXF0m6ym9ng/7FFP/2BPol08Wv/H6DaVPd7tPdSkwKOCuy4NbT+GmgCoEz3eQHpSb+gIj3pIn9GCWvRH/feeV4J7oSLPZZ4PPRqSssCjx6ioSxQXUuqpU6/tXM9az989QzxytkfFogr//3l28iSKn+NAjwqPeR4Vlv/yfMmTSAcR+jqveTTEQs+pnvwfeWpTgaXn3dGj3tyFnjUeDbmpq5iJX8cBnks8sz16fhXGdQ3F6eNRj/lyj4ZoPOXZyqPhOy94knW9R0N/fhnca5j2pb6GerXr2pei7ZfMc4mKX9HBVyLt89b0Z4/28fHBvY3UC65PZJ4N7iXnux5t+zmeco++///zaP8koyn7UdtJ72e95/TeWOr5k0fDdRLF76+TPBq2E7//v+4RrYvHeOuxn4bN4BMHrVPqeg/90POBR/tnhudbnkSJBxfHeqZ7tJ3f9uhTDdHPqn7WtXym53ue2k7xvOlZHeV9z689AACgESpg9cdYf5STpTG/+je6PdNzs0dFjYrRnp5YXKyqOP2P58ceFc8ac6xlKiRUpGjc8S88+iOuP/b5nlj8HPpDrwJBj1XhsMSjx2uYTGxrjwqrWz16Xed6/u3Rv7/UkyguvDUsQ2OfVTCpoFavd0OFt4oSRc99sUfFlYZwaDhATD3Cet4qz50eDdtR8abto7HSyRTe8ddTr662jdqJvdh19XhriIyWxQc1LaUeWBV9pZ6BWtBMTe3x1vd8v0cFXazQs8gz0ZP4uMZo36qQv8ijA5PfevRefc+T2HNfV493U/fjLR49x/Mevdc1zEcHTSrat/HE4veXnvsjj95Her9qWJdoXVx49/Wc49GyJz3HRdGBqMTPpfektpd+pvTe1Gus9OhnIhZ/j/rZ03p9XRX+OnDWAfG3o7a2gX4eVMTrORJfe7zt9bOv7aH34XUe/RwDAIBGqHBUcZCsgz36w3t3cG+jozxaflNwLxT/kf6vJ7GYVvGk5SqyEz9eV/Gh5foasfg5NBxGPdqxL3hUED0e3AupOFMSqRdShZC+R/XUx+LCe4qn9idyDRXeOvlUwwJiR3q0XEVILP4+VHwlUk+0lje18Nb3raJTXztWV+EdHwTUpu+7V60k8ymktvFnHvXGqgdVRWBzNKfwVmGp9h4eOcKj+3pNiY9rjA6kalNPs/59Ym9uXYV3U/ajDgC1TO+nxPe6PpXRJwePBfdC8ftLn7IkHrDFtC5xVpOGtl/8XDo4Svz5UBGv5fokKxY/j3ru+2tB5FCPluuAWNs3tr2n9nPo+9BBC4A0x1ATID1paIX+ECdLBZCoNy/Rox4NPdG44NpUpKtIjqn3W9R7qI+0Y/HyxF66mHoqNRQk9q5HQ1b08XxcbKtwUETFZg+PeuA17ELfp048rO12TzykIhnq+VQvYEwnnkriaz7cs86j15zoxui2qfR9q+jazfNNLaiHvkd9ClCbhhyohzYxjfWKa9iQTqzUQZNmVdF2fM7T2xOLi031JKeavpaGJsUnWepWn3povzeF9oPob5BOFtVBR1wwN3beQlP2ox4rKlIT3+s6WNJ21MFk4oGf6P1f14FSc6hQT/z5eMujn+u6fpbUW60hO7H4506z4CRuXxXYOmBNfA693kGevYN7ANIWhTeQnvTHeYuwmZThHhXL8bjTROoVHeqpPSZYH2sniosN9aYmiperYK4tcdhBTB+FazhAPLZcX1cfy+ux8RhbFZnxQUJdvYvqSW8Kje1NpGEEkvia1bOo77n2WGi9luYWWhqHrkJUM8/U9/u0vn2pHn/1ICs6MTEZGoag7aXx8SrANMxFw01UEKt4FfWIigq2VNPBkIabHO1R7+xhnkmeplJxrTHhGkqhba99oE9apK73Q6Km7Ef9XIjGTdemZXqfJvYyS1Pfew2p/b4UvTfr+llK9udRtC7xOdTbr+fVEC71smvmo+94Env5AaQBCm8gPakoUE+gCubWktgDmKi+5bUL92RpjmcNddHJapqJ5RCPis14+Etdv4cSe9yTkcxrbu7rb4gOJFR0q9j9vhbUQQVyXftSQxrU66rULrrqs59HBzY6eBH1oKrXVgWmnkfFmIpyncja1F7oZKnQ1tdRz7D23YOeptB20EHHSI/GyevTGL0f4iEqjf1dSvV+1KcDiZr63mtIU36WWvLzqE+19MmRtqWGnehTmEc8L3vqGtYDoJ1QeAPpKR57qnGvyVBPt8Zlxz18iUZ5VNjVLjBSoa5hIpotQr2R8cfmx3j0+vSxvz56/5tHRaKK1rak16CiL/HkPdEwjcZ6WRuimWPUW3u5J3GceSzel6mYZzuezi+x6FJxpaFG2u4aQqHxxeoZby36BEWFvWaS0b5UT3NTaPy9ikHN/X2D5wmP3g91fVpTl6bsx/g5E0/2jennQu9TnbjYHK3x89QS+nnSttTYd32/OkFTBbg+nQCQJii8gfR0l0dFg2ZDqO8CHSrK4/Gw+oMrenwijT1Wz2LiyY6ppKnSNE91TCeBaXy3xm9rlgqJe+wSf99obLJ6v9uSppxTwafXnEgzsrSExq/r5D9NZ5d44mdM+1LDF1QM1zePd7K9uBo/392j2UASqQdZ0xFqCkpN0acx4K1J7zMdaKi3v6nqej+IruqZjKbsR804InruxK+n8fHqZdeBQ/w+bap47Lb2R3uLhxkleiO6TZzRCEA7o/AG0pP+qKuHWLNmqIjViVc/8ahYVZGnIQbqaY3HcKqAUHGtYuQvHg030MlmGjusMaLJzhHdVBq/rLHEmipNV/HTCXLqedMUgDG9dvXEq2DS69MUgnr9GhvdljRPtIbwaIhLPA2dxsKqR17TILakB/MBj567rsvIaxyzxkLrEwDNrqErVKo3UgdOKgg1p7S2n16DHtsQjeXV0BWNj5/s0XtC34cubqP3hU4q1dh6TdVXV+97XXT1Uw35qB3NKV0fTQeor5c4o0uytA30HtG2v9BzVrSsrm1Xl6bsRz1O01hqOI561TUdn+aU15h4naBY+0C1KTTcRwdUmjlIvcyajWWMpz3od4TGzOtnUBdW0sWYtE10EmrizC0AAKABGhusIkgf7avIVe+cTtjSONval4vX7Az6w6ux1HqcxhDf44lPcozF074lTnkn8RRotXtt65rSLX4O9eBqxggVlSqmNBNDPNVcTAWgekf1uvUYDVVQsVXX69DX0DJ9zdrqen31fS+i5YnTv0k/j4pkFV2aaSS+8IoKtmQu2R5/vbqmzIunbqzv9egESxWa6o3WyXHqKdcwDRWxKryTHe6i51HhrXG92s96LhWSmu9Z1Buu16DCq6Ge9Hi/1hcVqtLQ95wo2ceJhqnowkw60NA20H7SUBH9+7qm7Et870lT9qO2gQpujXnX0BJdxEkHJvVdQCfx/ZWo9muTL3u0P1Xgan08M0tDz1X7Ij/1fY9S19eU2s+hS9nrPaCfeX2P+tRDB2M7eAAAQIaLi6y6CsxMoxMF9b3Unp4OmYX9CCDtMdQEQC5JvDBQLB5b3JRLnqN9sR8BZKRkT+gBgETq8dYwlv09iR95pzuNd9VH8ZoJRDOEaEyupmDT/X08iRfhQfpiPwIAgJyRqUNNdDXH1z0aF63x0Zo5Rieh6uqSyBzsRwAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAKBN5UW3aBv5nm08yzw1WgAAAJAFVFP28HzkqdICbI7Cu21t55kRNgEAALLOKM/7YRO1UXi3rT6ehdOnT7devXqFS1pRWVmZTZ061caMGWPFxcXRUqQr9lfmYZ9lFvZXZmF/ZZZ58+bZzjvvrGZfzyI1sDkK7+S94NkvbG5CQ0b6eZJ5k/X2LHLWu7eareejdWb3zqu0Nz4ptV23GmwnDiiwbTpHK5GW9EdmypQpNnbsWP7IZAj2WWZhf2UW9ldmmTNnjg0ePFhNdTIuVgOb6xDdonHf8vSvlbc9z3nS5siuvNrs5PdrbOQrZjfMyrMZa/sEt7qv5VoPAACAtkfhnTydELkgIerl3slzhydtnP5hjd073+yKx961+ec+abMueNrmnfdkcF/LtR4AAABtj8K7+c7wzPM8EdxLAx+uM7tnQZ5d/tf37NKnZlj3dRXB8h5rK4L7Ex5/L1ivYSgAAABoW4zxbp4tPXM9N3gmaEE9ungSR1b39MwoLS1tlTHeV5TmB8NK1NMdF92JlnUptAE3Hm7nD6uxSwcz00+6KS8v//xEoqKiomgp0hn7LDPU1NRYVVWVrV+/3l555RXbfffdraCgIFqLdFVZWcn+ShN5eXnBPsjPzw/addHJlSNGjFCTMd4NoPBunjM9t3iGelSA10dF+fiwudGkSZOspKQkupc6d3bcMRjTreEl9Rk68RAb3WWRnb7hnWgJAGS3Lbfc0rp27WodOvAhL9Bc1dXVtmbNGlu1alW0ZFMrV660cePGqUnh3QAK7+aZ5tEE8TrhsiHt0uOtMd0aXlIbPd7pjd7TzMM+S2/q6dYsUioWNIVrp06dbN26ddali341IxOsXbuW/ZUm9LOzZMmS4CC2T58+m/V80+OdHArvptvH80/PQZ5ntaAJWnU6QY3d1uwlOpFSY7pru+LwUTb+yC/Yh7sbUwumIabOyjzss/SmoQofffRRUHTrd6567NRbpx5wer/TH/sr/SxevDgovrfZZpvNhv8wnWByeCc3nU6q/Njzj+BeGlExfVK/Gi+utw+KbPVwi251f8I3tg/WU3QDyAUVFeEnf50780sPSIX4Zyn+2ULTUXg3jbqpj/L81pOW8/LdOTLPTuxvQc+2hpUMmXiI9bvxiOD+QT3C9QCQCzTUROo7GQxA08Q/S/HPFpqOwrtp9NGJPk++PriXhop8j/5+u7xgOInGcm/fZZEN7xquW1+dF6wHAABA26MMy1IaTqITKDV7yb3bVAbL/rnSsyJoAgCQUsOGDbMJExqaYRfpQvtJ+wttj8I7B+zatcYO6RG2r5wd3gIAkqeT18d/anbGh+FtW12I7MQTT7SvfvWr0T20lIZKKI888ki0ZKOzzjorWJeL25tCvO1QeOeISzXjuPvHcrOXV4ZtAEDDyqvNTn6/Jpgx6tpPqm3Ke+uCW93Xcq1H69qwYUPUSo2hQ4fab3+rU7U20lR5999/f7CupTQbSzqcfJjq7YbUoPDOEXuWmH2te9im1xsAknP6hzV27/xwmtb55z5hn54/2ead90RwX8u1vj3ddttttt1221nHjh2DKd5++ctfBtMoyvjx44NCsvaJcL/4xS+C+Zbj5R9//LEdddRR1r17d+vWrZsdeOCBNm2aLlcReuGFF4Ke4GeeeSaYs14zW2y99db2hz/8IXrERppi84wzzgieSxeK++EPf7hZAdjQaxb1OOtCLFdeeaWNHj3aevQIP7LVlUdPPfXU4Hn1/Po6F198cZN7ak8++WR78cUXg6kmYw899FAwN/V+++0XLQm98cYbdsghh1jfvn2D+cR33HFHe+CBB6K1obi3+C9/+Yttv/32wdSir732WrDulltusUGDBgXb7Otf/3qwzbQtZ82aFayPt218P1Z72I5e35e//OXg+1b2339/e/XVV6O1oXi7XXHFFTZw4MBgX4q2v/ZD4nbTforde++9dvnll9vs2bOD16Lokxa0DgrvHBL3ej+9zH+ZrA7bAJBLqrzWXFDu8bqjsfxnhdk9C/Ls8r++F1wbofu6sBdTFyjT/QmPvxesf8kfV9e/ryv6+qmiwmzixIl29dVX24wZM+zWW28NenJVcIuKMF2wTcVdTMW2enZPOOGEoMBauHChfeUrXwnmOVcx+r///S8oilWA6poTic4555wg77zzjh1++OFBAfvhhx9Ga0O//vWvbYsttrCXX37Z7rnnHvvzn/9sP/vZz6K1jb/m2GOPPRYUo/r3//3vf4Nl559/vj3++ONB8arnV2H5m9/8JljXFEOGDAmK4LvuuitaYsFrOO2004JtkkjziH/729+25557zt5+++1gm2rbPf/889EjQtpWN9xwg919993B96VtqEL8vPPOC7bZW2+9ZUcffbRdcMEF0b9oGh106CBDxfY///nPYL5sXT9g6dKl0SNC8Xb7+9//HuxLufDCC+3RRx8NCmxtNx1AaD/F9Lp++tOfBgcI8+fPD6IDBiAbaDrCGv8B9d99rc+Pcmv++te/Brexfd7w37rP19R8c1q0AGmjrv2F9MY+S29r166tmT59enArVVVVNe8vWhH8Dkw2hc9W1Sz7wmM1NSMe2SxLd3gsWF/Xv6sv85v4VvFCr8aL4OjeRvqeOnXqVPPUU09FS0KTJk2q6dq1a011dXVwf5999gmeIzZ16tQaLy5rPvnkk+C+F7w1u+22W9CO6d+OGDGixgvk4L4XmTpcqHnooYeC+1JWVlbTsWPHmjvvvDNaUlMzdOjQmi996UvRvdBtt91WU1RUVLNmzZqkX7O+3+HDh9eUl5fXLF++PNhvq1evDp7njjvuCB4T22OPPYKvmyx9H35AUPPEE0/U9OrVK/g+3njjjZri4uKaxYsX17u9Ex188ME1p5xySnQv3IZ63nibxvbaa6+a73znO9G9kBe4wWM//fTT4H68beP7MX1Pet76VFZWBtvMD6KiJRu3m9bFtN213bzQjpaEvvjFL26y3fS1ktmOtX+mEvlBXvC9eFJ/hcAsQo93jol7vR9bYjZtTdgGANSt76oNn/d016ae7z6r22cc7XvvvRf0gn73u98NLuEdR0MKdIn8BQsWBI/TkAH1dnqhFNxXb7EX4/GlvYMeVPXGJj6HeqzVa5o4FEO8WItaZl7MBUMz1GOeaO+9945aIfWmewFtXpQm/ZpFXys/Pz+6Z8G/1/NouEWi2vcTnzeOhn/UpuEjGuqiXmk/eAiG2ugKp7UtW7Ys6LHWkBcN09DzPfvss8GwjET6t/E2jXmBWuf2aI5333032G5bbbVVcCVPDRvRPq39OmpvNw0j0nZL1etAy1F45xiN895ji7D9S8Z6A8gxvQpqbO6e1Tbf67XGct4gs8Vbdvz8KsC1afmSLTraT/1xdf37utK7KPrHLaQT+ERjf1U4x9HYbBXMGjoi3/nOd4LHqsDUCYQaupE4flfrNKwk8TmUDz74IBhjnaiwcNPtoGEZ8euI1R6qEUt8bGOvWTQcoi71PX8s8XnjPP3009HajVScnnLKKXbzzTfbgw8+GBT/ddG20rANjUPXUBw9n4apqJhNVNfr1Wtt7PXGl8KvCTrjN0p8fu23r33ta8FYbQ1l0XARvY6ePXsm9TqksdeBtkPhnWP0s3dpdB7KI4vN3g87QQAgJ+T778B+Xvz2K248pw8wKyvoYL/ef+voX2/q9jFbB+t/6I+r69/XFX39VFAvrnps1ROsEx1rp6CgIHiceq+/9a1v2aRJk4KebxW/GrMc22233YKe2QEDBmz2HOrRbqr//Oc/USuk+zrZUL3Byb7muqinV73s8XjvmIrQRHU979B6ZipR4a2TIDW2ed99942WbkpjuX/0ox/ZN7/5zeDEyuHDhwcHJclQL3ld2yNRvI3nzp0b3IrGWCf2/mvMuD5ZUPGvk1v1vNpWS5YsiR5RP33/2m61v+5LL70UtUJ6TFVVVXQPrYnCOwdpTu9du4YDsX71WbgMALApXYjspH41Nv7I7e2Kw0d93vOtW92f8I3tg/V6XGvSMIzEHlxFQ0EuueQSu+yyy4JZQlQMaijHww8/vNkJfDohUAXk9ddfHwypUDEe09zV8o1vfMP+9a9/Bc+rIu3SSy8N7jeVisSLLrooeD1//etfg5Mm1ZusnlgN00j2Ndemf6tCWf/2iSeeCE7q/PnPfx58veb25mrmj8WLF282O0ginSSpHnqdWKmeefWA1z7ptD46YVEHOzfeeGPQo6+TG++7775obUiFsYp5bScdAL3yyit23HHHBQcoMc1wovu33357cNDy73//277//e9bp06dokfUT9v99NNPD55f+0PbXCdbvv/++9EjQjqwUbGvglwFvd5zaB0U3jlIv6N+EXUAPLjQ7JP1YRsAsKk7R+bZif3Ni+8v2IAbj7BhNxwa3Oq+lmt9a3v99ddtl1122SSaOk7FsYZK/O53v7OddtopGLerIq/29HrqJVWvrmYjURGeSNPkqRdZtyrKt912WzvmmGNs5syZQS94U6l3WOOid9999+Brqbf9mmuuidZa0q+5LpoNRbOpqDDV8y9fvjwohBOL1KaKx23XR8WyhoNoLLlmERk5cmTQ+50MPU4znSjqLddsMtdee220NqSeax14rF69Ovj0QTOmaEx54qcNGlKi4TCanUa93Zp9RQc3yX4ioe2v16L9oe22YsUKO/PMM6O1Ie0nFfNHHHFEMOQnPiBD6rX+bwwk0gA2P1hetMlYttaieTqnTJkS/LLQR32JqmvMdnrN7N21Zqf4H4+7t41WoN00tL+Qnthn6U1jY3XymYYaaB5lDbPQ9HA6OS0eW5ssXanygYVmCyu8WC00O7Zv2COO1pPM/tJBhXrwNc1gJlDxrDm4P/3006QONtJN7Z+pRHPmzAmmOXQ6IlisBjZHj3eO6uCHXD8fErYnLTD7jAtcAUC9VGRPGG52x8jwlqK77WlYinqgNcxEs3xofnANodFFdYBMQeGdw77jx6TbdjKrqDG7lrHeAIA0pl5vXejmS1/6UjD0Qxe10Rjqww47LHoEkP4ovHOYzq6/JBrr/bv5ZvPKwjYAAOlm1KhRwXj0lStXBmOiNSOJxiZnEo3N19SBmTjMBKlB4Z3jjuljNqKjWVmN2cTSaCEAAABSjsI7xxX4O+DiaKz3b+eZLdp0Ln4AAACkCIU37IR+ZoOLzdZXm91IrzeALBHP71z7qoAAmif+WeJKmM1H4Q0r8nfBz6Je71/PM1taEbYBIJPFlzjXFGgAWi7+WYp/ttB0HLK0rbSZx7u2DVVmI/5nNr/c7NKhZlcMj1agzTAndOZhn6U39c7NmzcvmAu6V69ewZX+1q5dG1zNr6nzeKPtaR5v9ld60M+Sim5d1VLzquviSrV7vZnHOzkU3m0rbQtvuanU7LxPzEryzWZ/2W8LohVoExRxmYd9lv4Si2+1169fHxTgfFSe/thf6ae+olsovJPDO7ltpXXhva7KbNjL/tNSYXblMLNfMNtRm6KIyzzss8ygAq6qqiroPf3Pf/5je+21lxUVFUVrka7Ky8vtpZdeYn+lARXaGl6Sn59f70EQhXdyKLzbVloX3qIL6Vw006xHgdmsPc22oNe7zVDEZR72WWZhf2UW9ldmofBODoOmmqan5w7PPI8uNzPLc7ona/xoQFh0L6s0u1PfJQAAAFKCwjt5XT3/8mzt+b5n2+h2uidrqIf7nEFh+/rScPgJAAAAWo7CO3kXeDp7DvO86FFv9389//RklbMHmm2Zb7aowuzu+dFCAAAAtAhjvJP3rudNzwrPtz2rPJM9l3rWeurSxaNiPaahKjNKS0vbZIy3TkyZOnWqjRkzpsknpoz/LN+unVtgAwprbPqu5daRQ7RW15L9hfbBPsss7K/Mwv7KLJo9aMSIEWoyxrsBFN7JW+/R9nrYc5tngOd2j3q9j/bUZYJnfNjcaNKkSVZSUhLdS0+r8orstC0OtA15BfbD9W/bweXq4AcAANjcypUrbdy4cWpSeDeAwjt5OplyqWeoJ76243c8j3j6ehZpQS0Z2+Mtl8zOtxvnFdjgohp7b5fy4AqXaD307mQe9llmYX9lFvZXZqHHOzkU3slTl6/y1eBeaJRHJ1fu7nlVCxqR9tMJJlpYbjb8ZbP11Wb/t63ZD/pHK9AqmDor87DPMgv7K7OwvzIL0wkmhz7M5GlGEx3K5Qf3QiOj26wch9G3yOy0qNi+erZZpRfgAAAAaB4K7+Rd79GQkls9mkpQPd8TPQ96svbI7oIhZkV5Zp9sMPtjXYNpAAAAkBQK7+S97TnEo2Elat/redJzmidrDSw2Oznq9f7lZ2ZVNWEbAAAATUPh3TTPeb7k6egZ5vmpp76pBLPGRUPMCvLM3l9n9iijtgAAAJqFwhuNGuqHGSdokI27arZZNb3eAAAATUbhjaRcPCR8s0xba/bEknAZAAAAkkfhjaRs3dnsmKjX+8rZZjX0egMAADQJhTeSdsmQcOL3N9aYTVkWLgMAAEByKLyRtFFdzL4dXfeHXm8AAICmofBGk/xCF8x3/11lNnVF2AYAAEDjKLzRJDt2NftGz7B9ZVZerxMAAKB1UHijyS7VDObuxZVm/6LXGwAAICkU3miyL25hdnCPsK2x3gAAAGgchTea5dJorPezy83+typsAwAAoH4U3miWL5eYHdAtbOtqlgAAAGgYhTeaLZ7h5KmlZm+uDtsAAACoG4U3mm2/bmZfKQnb9HoDAAA0jMIbzZaXt3Gs91+WmL27JmwDAABgcxTeaJEDu5vtvkXY/uVn4S0AAAA2R+GNFkns9X54kdkH68I2AAAANkXhjRY7tKfZLl3Narx9NWO9AQAA6kThjRZTr3c8w8n9C81mrg/bAAAA2IjCGylxZC+z7TubVXn7GsZ6AwAAbIbCGynRIc/s51Gv970LzD7bELYBAAAQovBGyny3j9nITmYVNWbX0esNAACwCQpvpEx+ntklUa/3/803m18WtgEAAEDh3RQTPJq4o3aGeRA5po/Z8I5mZb5lJpZGCwEAAEDh3URzPP1rhfIyQaG/oy4eErbvnGe2uDxsAwAA5DoK76bRpB0LakXLkOCEfmaDis3WV5vdqEMVAAAAUHg3kZeUQQ+3ysm/efbyoJZif1f9bHDYvn2u2bKKsA0AAJDL8qJbNO5gzxae9z0lntM8x3jGep711KWLp3PYDPT0zCgtLbXevXuHS1pReXm5TZ061caMGWNFRUXR0raxvsps1JtFtqAizy4ZVGmXDeaDgca05/5C87DPMgv7K7OwvzLLvHnzbMSIEWr28SxWA5uj8G6Zf3o0Y/VBwb3N6YTM8WFzo0mTJllJiWr37PZ40VZ2T6cvWOeaCrt71TN+FFIZrQEAANlk5cqVNm7cODUpvBtA4d0yEz3f8IwM7m0uZ3u8ZW2V2bZvFNmSyjybMLjSLhpEr3dD6N3JPOyzzML+yizsr8xCj3dyKLxb5rno9oDotjGqthe5Nim8y8rKbMqUKTZ27FgrLi6Olrata2abXfypH3EUmM3a06yr36Ju6bC/0DTss8zC/sos7K/MMmfOHBs8ODjBi8K7AZxcmbwbPPt7hnt29twR3b/Zg3r8aKBZdy+2l1aG0wsCAADkKgrv5A3w3OfRyZXPeDS85GueJz2ox5ZedP9kUNi+vjQ86RIAACAXUXgn7/selZD6vEsfo2h4yVQPGvHjgWZb5JstrDC7e360EAAAIMdQeKPVdS80O9uLb7nuM7Oy6rANAACQSyi80SbOHWTWxd9tc8vN7qHXGwAA5CAKb7SJXkVmZ0S93td8ZlZBrzcAAMgxFN5oMz8dZNbR33Gzy8zuXxgtBAAAyBEU3mgz/YrNTusftn/1mVklvd4AACCHUHijTV0w2Kwoz+zj9WYPM70+AADIIRTeaFODOpqd1C9s/3K2WXVN2AYAAMh2FN5ocxcNMSvIM5uxzuxRer0BAECOoPBGmxvWyez4vmH7Knq9AQBAjqDwRru4ZEj45ntnrdlTS8NlAAAA2YzCG+1i685m39eF992Vs81q6PUGAABZjsIb7eaSoWZ5fvvaarO/LwuXAQAAZCsKb7Sb0V3Mjuodtun1BgAA2Y7CG+3qF0PD25dWmT2/ImwDAABkIwpvtKudupod0TNsq9cbAAAgW1F4o91dGvV6v7DC7N/0egMAgCxF4Y12t9uWZmN7hG3N6w0AAJCNKLyRFuJe778vN3tlVdgGAADIJhTeSAt7lZjt3y1s0+sNAACyEYU30kbc6/3kUrO3VodtAACAbEHhjbTx1W5me28Ztun1BgAA2YbCG2kjL8/s0mFh+9ElZu+tDdsAAADZgMIbaeWg7mZf2iJs/5JebwAAkEUovJvnAE+VZ1ZwDykT9HpHY70fXmT24bqwDQAAkOkovJtugGeS55ngHlLusJ5mO3Uxq/b21Z+FywAAADIdhXfT5Hse8tzi+Z8WIPXU6/2LqNf7vgVmn64P2wAAAJnMSxw0wbWeHTyHesZ7TvREpwPWqYunc9gM9PTMKC0ttd69e4dLWlF5eblNnTrVxowZY0VFRdHSzFBdY/bFtwttxvoO9oM+VfbrrSqjNdkrk/dXrmKfZRb2V2Zhf2WWefPm2YgRI9Ts41msBjZH4Z08Fdu/9ezi0RtqgqexwluPUYG+iUmTJllJSUl0D/V5sXCg3dR5NyuoqbY7Vj9rvWs2RGsAAEA6WblypY0bN05NCu8GUHgnZ5DnDc93PS9ogUum8KbHuwWqasx2eqvQPt7Qwc7oV2U3Dc/uXm96dzIP+yyzsL8yC/srs9DjnRwK7+Qc6XnMo5lMYhofr+2nZad5fu9pjKrtRa5NCu+ysjKbMmWKjR071oqLi6OlmeXe+WYnfWBW7Ft61p5m/TLz20hKNuyvXMM+yyzsr8zC/sosc+bMscGDB6tJ4d0ATq5MznMeje3eOSF3euZFbRXlaAXH9jUb1tF/AdeYXV8aLQQAAMhAFN7JWe15t1YWeSqi9nIPWkGhv0MvGhK27/DDnMXlYRsAACDTUHgj7Z3Yz2xQsdm6arOb5kQLAQAAMgyFd/Pp5MqGTqxEihT7u/TCYNiY2e1zzZbpcwYAAIAMQ+GNjHBKf7O+hWarq8xupdcbAABkoGwvvDXjiM6urU3T+iXOUII01ynf7IJorPctc81WZf/1dAAAQJbJ9sK7vukSCzyUbhnm9AFmvQrNVvie+7UX3wAAAJkkWwvvfaPIlz3xfWV/z7me2R5kkC75ZufpUkbuxjlma/nMAgAAZJBsLbx1dcnnPTUezbGt+3E0J/eZns0u5Y70d+ZAs24FZksqzO7ULOoAAAAZIlsL7+GerTwaarJ7dD+Ol25W4nnIgwyzpRfdP9EedBM/M1tPrzcAAMgQ2Vp4axjJLI++v9ei+3Hme6o9yFA/GWS2Rb7Zwgqz/9PeBAAAyADZWnjHDvTsHTYDp3ne9DzgUa83MlD3QrOzol7vaz8zK+MwCgAAZIBsL7yv82jqQBnpud3zsmcnz0QPMtS5g8w6+7t3brnZvQuihQAAAGks2wtvjfN+L2zaEZ5nPWd4TvUc4kGG6l3kO3JA2L7mM7MKer0BAECay/bCO9EYz9/DpmkW6LgnHBnq/MFmxXlmszaYPbAwWggAAJCmsr3wnubR1IFf8Wj+7n94xEs2Wxw2kan6FZudGvV6/+ozsypNHgkAAJCmsr3wvsRzoudFz32e6R453PNq2EQmu9APoQrzzD5ab/bwomghAABAGsr2wlsFd2+PhpVoRpPY3R71hCPDDe5odlK/sP3L2WbV9HoDAIA0lQtjvHWJlfWe7T2jPcWeTzzMhZElLhpilu+309eZ/YUBRAAAIE1le+Gt7++XnuWedzzvRu2rPLlw0JEThncyOz7q9b5qtlkNvd4AACANZXvxqQJb0wde5Nk1ys88p3uu8CBLXDwkfDO/vdbsqaXhMgAAgHSS7YX3CR6N7b7V87bnLc9tHhXe4zzIEiM7mx3dJ2xfSa83AABIQ9leeOvEShXbtWmZ1iGL/HxoePvqarNnNKAIAAAgjWR74f2x56iwuYlvenSCJbLI9l18Z/cK21fOotcbAACkl2wvvK/z/MrzJ8+PPWd7/ui52nONB1nmF1Gv939Wmb2wImwDAACkg2wvvCd5jvD09Yz3TPDoWoe6gI4uqIMss/MWvnM1a7vTWG8AAIB0ka2FdzePerhLPJM9+3pUjikqxLfxbOlpim97XvEs82zwaBiLZk0p8iCNXBr1ej+/wuw/K8M2AABAe8vWwltTCH7NU1fZpQEIYzya2aQpVHBr6Mo+nu08F3g0Y8r1HqSRL/kh1de7h23N6w0AAJAOsrXw1gmVd4bNOt3l+W7YTNpUz58973lmeR7zPODZ34M0E4/1nuKHS6+uCtsAAADtKVsL7608ukplfVQ86zEtMcpzsOf54B7Syle6mX1VA44cvd4AACAd5EW32Wa1ZxePxmHXZWfPvz1dg3tNs8ZT6NHY7t96zvRUeerSxdM5bAY0xnxGaWmp9e7d+tOIl5eX29SpU23MmDFWVJR7Q9GfX5lnB08Pv+9Xdiy3Hbuk9/yCub6/MhH7LLOwvzIL+yuzzJs3z0aMGKGmLme3WA1sLlsL79c9ukLlvcG9zWl8t8Zn6xLyTbW1p6Pni55rPSq+NWNKXTSLymbrJk2aZCUlOu8TrUll9sVdvmLvF/S0vSrm2oXrXgtXAACAlFq5cqWNGxdcFJzCuwHZWnj/zPMjj8Zfz9SCBMM9//Xc5FHh3BK6JP3vPaqi12pBLfR4t7NnlufZEe8X+Ru9xt7YqcJGdU7fXm/2V+Zhn2UW9ldmYX9lFnq8k5OthXex51+eL3ju90z3yPaeYz0a4/0VT5mnJXRod49HBXUyFylXtb3ItUnhXVZWZlOmTLGxY8dacbE2Se7R1St3f8PstdW+4/1Xwf2joxVpiP2VedhnmYX9lVnYX5llzpw5NnjwYDUpvBuQrSdXqqBWb7eGgejy8DdGOTJapnVNLbov9Rzk0VATzQOuAl495o97kim60Q7y/NAyntf7oUVmH60L2wAAAG0tWwtv0dCPcz3qWtaVK/tFbS3TCZJN1cmjcePveDSGXMNZNFzlGA/SmK5kuWMXs2pvX/1ZuAwAAKCtZXPhnUgfeSwKm812iWdbj8Zs66qXO3rU473egzSmXu94Xu/7FprNYo8BAIB2kCuFN3LcUb3NRvkhU2WN2TX0egMAgHZA4Y2c0CHP7OdRr/c9C8zmbAjbAAAAbYXCGznj6N5mW3cyK68xm1gaLQQAAGgjFN7IGQX+br9kSNi+a77ZgpZOJgkAANAEFN7IKcf1NRtabLah2uyGOdFCAACANkDhjZxS6O/4i6Je7zvmmi0pD9sAAACtjcIbOeek/mYDi8zWVpvdRK83AABoIxTeyDnF/q6/MOr1vm2u2fKKsA0AANCaKLyRk07tb9a30Gx1ldmtXnwDAAC0Ngpv5KRO+WbnDw7bt8wxW1UZtgEAAFoLhTdy1ukDzHoWmC33ovs39HoDAIBWRuGNnNXVi+5zo15vTS24tipsAwAAtAYKb+S0swaadfMCfEmF2W/nRQsBAABaAYU3clqJF90/9uJbdBn59fR6AwCAVkLhjZz3k0FmXfPNFpSb/W5+tBAAACDFKLyR83oUhkNO5NpSs/LqsA0AAJBKFN6AO2+QWWf/aZhTZjZpQbQQAAAghSi8Ade7KJxeUK7+zKyCXm8AAJBiFN5ARBfUKc4z+3SD2YOLooUAAAApQuENRPoXm53SP2z/arZZVU3YBgAASAUKbyDBz4aYFeaZfbje7BF6vQEAQApReAMJBnc0O7Ff2P7lbLNqer0BAECKUHgn5yTP857FntWe1z3HepCFLhpilu+3760z++uScBkAAEBLUXgn5wDPE55DPLt4/ui5z3O0B1lmRCez4/qG7atmm9XQ6w0AAFKAwjs5x3lu8rzq+dgz0TPF810PstAlQ83y/PbNNWaTlwaLAAAAWoTCu/m6eRiIkKVGdjY7uk/YvpJebwAAkALq1EPTjfPc7dnD86YW1KOLx0u4z/X0zCgtLbXevXuHS1pReXm5TZ061caMGWNFRUXRUiTrvXV59sW3w+325KhyO7Bb61bf7K/Mwz7LLOyvzML+yizz5s2zESNGqKluK50ThzpQeDfdEZ6HPT/y3KMFDZjgGR82N5o0aZKVlJRE95DOrun8JXu5cICNqlxqv1r7b35gAACow8qVK23cOPVLUng3hDqiab7nuddzenTbGHq8M9xba/Nsz3fCbff30eW2X0nr9XqzvzIP+yyzsL8yC/srs9DjnRwK7+Sd6rnNc6JHs5o0h6rtRa5NCu+ysjKbMmWKjR071oqLi6OlaKrDp5k9tdTsgG5m/9g5WtgK2F+Zh32WWdhfmYX9lVnmzJljgwcPVpPCuwGcXJmccz13eH7iecGjS6worV89o939Ymh4+9wKs/+uDNsAAABNReGdHBXcuqbKnZ75CdH0gshye2xpdmD3sK0ZTgAAAJqDwjs5wzwallM7Wo4ccGnU6/23ZWavrQrbAAAATUHhDSRhn25m+0UT0ehqlgAAAE1F4Q0k6dLo843Hl5q9syZsAwAAJIvCG0jSmG5mX94ybP+SXm8AANBEFN5AkvLyNo71/tNisxlrwzYAAEAyKLyBJhjbw+yLXc10GZ1ffRYuAwAASAaFN9AE6vWO5/V+cKHZx+vCNgAAQGMovIEmOqKX2Q5dzKq9fTW93gAAIEkU3kATdUjo9f7DQrNZ68M2AABAQyi8gWY4qrfZdp3NKmvMri2NFgIAADSAwhtohvw8s58PCdu/n282tyxsAwAA1IfCG2im7/Ux26qjWXmN2UTGegMAgEZQeAPNVOA/PZdEY71/O99sYXnYBgAAqAuFN9ACx/U1G1JstqHa7AbGegMAgAZQeAMtUOQ/QRdFY71/M9dsCb3eAACgHhTeQAud1M9sQJHZ2mqzm+dECwEAAGqh8AZaqGO+2YVRr/dtc81WVIRtAACARBTeQAqc2t+sT6HZqqqw+AYAAKiNwhtIgc75ZucPDtsabrK6MmwDAADEKLyBFDl9gFmPArNlXnT/Zl60EAAAIELhDaTIFl50nzsobGtqwbVVYRsAAEAovIEUOtsL75J8s8UVZnfR6w0AABJQeAMpVFJg9uOo13tiqdkGer0BAECEwhtIsXO88O6abza/3Oz3C6KFAAAg51F4J29fz+Oe2Z4azwQPsJkehWZnDgjb13xmVl4dtgEAQG6j8E5eV890z4Ue+jHRoPMGm3Xyn67SMrM/8G4BAACOwjt5T3su9jzs8XIKqF+fIrMfRr3eV39mVkmvNwAAOS8vukXTzPLc62lsuEkXT+ewGejpmVFaWmq9e/cOl7Si8vJymzp1qo0ZM8aKirwSRJuaV2426o0iK6vJs//busKO691w9c3+yjzss8zC/sos7K/MMm/ePBsxYoSafTyL1cDmKLybJ9nCW+vHh82NJk2aZCUlJdE9ZLPfdtzR/lY83AZUrbHb1jxn+dFyAACyycqVK23cuHFqUng3gMK7eejxRlI+KzPb/s0iq6jJsz9sU2Hf7VV/rzf7K/OwzzIL+yuzsL8yCz3eyaHwbp5kC+/aVG0vcm1SeJeVldmUKVNs7NixVlxcHC1FWzv1A7P/m2/2BT8Me3s3sw71/NSxvzIP+yyzsL8yC/srs8yZM8cGDx6sJoV3Azi5EmhlFw+xYIjJu2vNHl8SLgMAALmHwjt5mk5w5yj6zKtf1B7tAeo1opPZsX3D9pWzzWo0CzwAAMg5FN7J283zZpT+nh9GbU0zCDRIvd4aYfLmGn/DLAuXAQCA3ELhnbwXPKqdameYB2jQdl3MvhsN679yFr3eAADkIgpvoI38fGh4+7/VZv9YHrYBAEDuoPAG2sgOXc2+2StsXzU7vAUAALmDwhtoQ7+Ier3/udKzImwDAIDcQOENtKFdtzA7tEfY1gwnAAAgd1B4A20s7vXWOO+XV4ZtAACQ/Si8gTa2Z4nZ17qHbXq9AQDIHRTeQDu4NOr11pzer68O2wAAILtReAPtYN9unpKwzQwnAADkBgpvoJ3Evd5/XWI2bU3YBgAA2YvCG2gnB3Q323PLsP1Ler0BAMh6FN5AO8nL29jr/fBis5/MLLA7O+5oV5Tm20frwuUAACB7UHgD7eiAbmY9C2qC9j1z82zG2j52w6w8G/mK2cnv11h5dbAKAABkAQpvoB2d8VGNLaswu+Kxd23BuU/YrAuetnnnPRncv3e+2ekfhkU5AADIfBTeQDv5cJ3ZPQvy7PK/vmeXPjXDuq/zCtz1WFsR3J/w+HvBeoadAACQHSi8gXbywEKz4spqO2vqx9GSTWm51v/8U7Pnl5u9tdrssw1mayrNaugIBwAg41B4A+1kUYVZ/9UbPu/prk093718/Z8Wm41522yX182Gvmy2xb+9YP+nWb//mI16xWzvN8wOn2Y2bobZuV7DXzHL7PY5Zg96YT9lqdn/VlnQa77Uv0wVBTsAAO2GwhtoJ30KzeZv0dGWdfFGHbR8sa/vlm9W4klU4QX0Qi+k3/eC+iUvrJ/yAvsPXmjf7AX3eC+8z/YC/FgvxA/2gnxPL8x1smYvL9QLXjTr7oX7CC/gd3vN7CAv6L/3ntkZH5r9fKbZDaVm98w3e3yJ2T9XmL27xmxemdmGqugLAwCAZsuLbtE2ensWOevdW83WVVZWZlOmTLGxY8dacXFxtBTpQr3QKoh1IqXGdNd2xeGjbPyRX7APdzfbprNZZbXZSi+AdTLmssqNt8tr309ox7eVKejp7uSH6d0LzHr4cUKP6Da4X9eyhPtb+kGDpk7MRvyMZRb2V2Zhf2WWOXPm2ODBg9Xs41msBjZH4d22KLyxCU0ZqNlLdCKlxnRreIl6um8fs7VN+Mb2dmJ/s99v17IfU40HX6uCvVZBHhToCctqF/C6XZuC6QzVWd+tkeJct3UV9YVp/JmcDpzunVdpb3xSartuNdhOHFAQHCAhPbG/Mgv7K/NQeCeHwrttUXhjE5qnW1MGavaS4soq67tqgy3csqOVFeTbSf1q7M6ReVbUjsWnXl/tHvS6etRrL1vht6mYgryrV+0NFefBstr3PV1asZe99j7r5/tsQRrtM2yK/ZVZ2F+Zi8I7ORTebYvCG3VS786keZX2+iel9sWtBtu4DO/dqa4xW6UiPKEgr7Ngr3Vft2UpGBZT6L/ZGirOdRsU8bUeo575/EZ+K8afUmgaSH1KoZNj9SnFr/ff2sYfmZpPKZA67K/Mwv7KXBTeyeHd27YovFEv9ldofa1hMXUV57WHyeh2VYpOAA2GxUQFee3iXGPlry1tfFz+33YwG9YxWoh28+kGs0OmNb6/puzI/koHn64PTwhP9rwXpBcK7+RQeCfvEM+vPKM8fjxut3uu9zQFhTfqxf5qGZ18qiEuDRXn9Q2TacrJp4X+hRae+0Sd00CqZ67fjUdYRQGfhacL9ldmaWx/DfD9ddFWHWzC8Ggh0gaFd3L4bZOc3TyPe57x7OyZ4LnKc5YHQBpQ7dSryGxkZ7M9S8wO7ml2bF+zsweZjR9mdss2Zn/ww+andjR7aVez9/fwo+C9zcr3NVv1FbNZe5q98UWzf+xk9qfRZr8daXa1/3G/wP+OnNLf7Fu9zAb482scfl1Fgejk2D6rN0T3kA7YX5mlsf3Vb82GYCpVIFPR452cBz1be3YP7oUmer7nGeJJtr+MHm/Ui/2V/sZ/anbtJ9U277wngiKgtrhH7vwRHewXQ6OFaDdXzja7YWbj++un7K+0oP11YxL7ix7v9ESPd3Lo8U7O3p4pYfNzuj/Iw69rIEcc19cPkAo6BCd61UXTQGr9uH5mHfNJe+dE3w/J7C89rpM/nrRvTkpyf+mTLCBT0eOdnHLPOZ7fBPdC23ve9ezl+a8W1KGLJ/EUkJ6eGaWlpW3S411eXm5Tp061MWPGWFFRUbQU6Yr9lRlO+zjf7luUX+/c68f3qbK7tuZSn+mC/ZVZ2F+Za968eTZixAg16fFuAIV3choqvL/seVkL6qCx4OPD5kaTJk2ykpKS6B6ATFLhvzbv7LiTPVc8dLO51w8om22nb3jbCpMefYbWxv7KLOyvzLVy5UobN26cmhTeDaDwTs5szyTPZcG90BjPcx6NNJulBXWgxxtJY39llo/X59n9C2vsrVnzbOdhA+y4vnm2dScKgnTF/sos7K/MQ493cii8k6OTK7fy7BHcC13n+b6HkyuREuyvzMM+yyzsr8zC/sosnFyZHE6uTM5Nnl0913i285zgOdtzrYdDcAAAADSKwjs5r3qO9BzsedtzpUfDTnQRHQAAAKBRFN7Jm+zZyaPPuzSFoObxBgAAAJJC4Q0AAAC0AQpvAAAAoA1QeAMAAABtgMIbAAAAaAMU3gAAAEAboPAGAAAA2gCFNwAAANAGKLwBAACANkDhDQAAALSBvOgWbaOPZ+H06dOtV69e4ZJWVFZWZlOnTrUxY8ZYcbEuuIl0xv7KPOyzzML+yizsr8wyb94823nnndXs61mkBjZH4d22tvPMCJsAAABZZ1fPm2ETtVF4t618zzaeZZ4aLWhlPT0q9Ed5lmoB0hr7K/OwzzIL+yuzsL8yiz7Kn+7p51moBUCu6e1Rga9bpD/2V+Zhn2UW9ldmYX9lFvZXEji5EgAAAGgDFN4AAABAG9CYY2QvjeGv8LwY3SK9sb8yD/sss7C/Mgv7K7OwvwAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAASIlDPG95yjyzPOd7kL729Tzume3RVb8meJC+TvI871nsWe153XOsB+np255XPMs8Gzwfe67yFHmQ/g7wVHn0twzpSX+z9LerdoZ5gKy3m0fzZ17nGeU50aM/Nmd5kJ50oHS152jPfA+Fd3q733Ou50uerT0XeKo92n9IP2M8Kr6396gQ+KZnkedWD9LbAM8cz988FN7pS3+zSj39aoVrxSAnPOhR706iiR79UGhye6Q3/XGh8M48T3seDZvIADd5poVNpCkVbboQiw5s9TuRwjt9sX+agEvGZ5+9PVPC5ud0f5BnaHAPQKp18ywJm0hz+iTwYI+GCyF9/cqz1nN9cA/pTj3c6uCLP6HYy4M6UHhnn/6eBWHzc/F9rQOQWuM8GuJ1Z3AP6WqNR+e9TPe84NFwIaSnQz06b0I/WxorjPT2P88JHu2373vUCfEvz4EeIOuVe34UNj+nsY365bVncA/pjKEmmeUIz3qPTrhEetN4/C94VMypM+JyD9KPPp3VGPyvBvdCDGXIPP/0PBM2geymmTGuCJuf08lFnGGcGSi8M8f3PDpxWScwI7Ood67S0yW4h3RypEd/r7R/4ujk5XjZyR6kP51b9mHYRCKGmmSf/3i+HjY/N9ajcVcqygG03Kmeez0qunWLzKITzfX3jykF089znh08OydEw7jmRe3HPEh/u3o05hvIepriTNMJXuPZzqOeHX0UznSC6aurJ/4Doz8u+iOj9mgP0o/GBqvn7YeexKmzenuQfi71HOTRUJNtPBo7rKEmFHCZg6Em6e0Gz/6e4R797brDo08pDvcAOUEnOLzt0YlE6uXWdExIXxrLqI9Ra4c/NOlJ+4X9lTk0O8YHnnWeVZ53PD/zdPIgM1B4p7eHPPpUXTWHxufrUwsNcQUAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAIDtN8NTUkWEeAAAAoEEdolskZ46nf62UegAAAIAGUXg3TZVnQa1oGQAAANAgCu+m6edRD7d6vv/m2csDAAAANCovukXjDvZs4XnfU+I5zXOMZ6znWU9dung6h82ADnQGeOZ6ND4cAAAgG6im7OH5yMNogHpQeLfMPz0bPAcF9zanEzLHh00AAICsN8qjTkrUgcK7ZSZ6vuEZGdzbXO0e716e6a+++qr166dRK8n5eenP7c11b9rxPY+37/f6frQUUlFRYf/6179sn332scLCwmgpWoJtmnps09Rjm6Ye2zT1cmmbLl261HbeeWc1+3oWqYHNUXi3zHPR7QHRbWN6exaVlpbaoEGDwiVJuL70ento8UP29e5ft18N/1W0FFJWVmZTpkyxsWPHWnFxcbQULcE2TT22aeqxTVOPbZp6ubRNFy9ebH369FFT/1usBjbHyZXJu8Gzv2e4R4d0d0T3b/a0qhGdRgS3MzfMDG4BAACQeSi8k6eTIu/zaNzSMx4NL/ma50lPqxreUbW+2ewNs62ypjJoAwAAILNQeCdPg6s1PkSfFeljFA0vmeppdXHhXV5TbvPK5gVtAAAAZBbGeLetZo3xloPeOciWVi61G0bcYF/t9tVoKRiTmHps09Rjm6Ye2zRUU1NjVVVVwUl8ardEeXm5vfTSS7bXXntZUVFRtBQtkS3bNC8vLzg5ND8/P2jXhTHeyaHwblvNLrx/+OEP7bU1r9mZA860k/udHC0Ff3xTj22aemzT1GObhkX3vHnzbNWqVdGSltHzrV+/3jp16lRvcYWmybZtuuWWW9qAAQPq/F4ovJPDT1bbanbhfW3ptfbI4kfs4O4H21XDr4qWgj++qcc2TT22aerl+jZNLLp79eplnTt3bnFhV11dbWvXrrUuXbpYhw6MRE2FbNmmer+tW7fOlixZUm/xTeGdHArvttXswltFt4rvbTttaw+OejBaCgqa1GObph7bNPVyfZtWVlbaRx99FBTdvXvrT0vLqUhUIa/CisI7NbJtm6q4VvG9zTbbWEFBQbQ0ROGdHH6yMsRWHbcKbmdtmGVVNVyJFQBymcZ0i3q6gbYSv9/i9x+ajsI7Q4zoGM7lXVZTZvPL5wdtAEBu0kf/wlhstKX4/Ra//9B0FN4Zonthd+tW0C1ocyEdAACAzEPhnUGGFQ8Lbv+46I82eelkW1+9PrgPAACA9EfhnSGmrZ1m76+ZEbTfWP6aXTb7Mjv87cOC5QAAZJoTTzwxGLpw9tlnR0s20vJ77703updaL7zwQvD8s2bNipa0nblz59pPfvIT23bbbYMpBnv27Gm77rqrXXrppTZnzpzoUchmFN4ZQD3b5354jvWft4Xdd/E4e/nYC4LbHnMKg+X0fAMAmkt/Q/Qp6t0L7rZnVz9rG6o3RGtan4rPO++8095///1oSfZ65513bKeddgoK/yuvvNJeffVV+89//mNXXXWVLVy40CZOnBg9EtmMwjsDTF0+1ZbXrLArbzvMRs/sHyzT7YTbDwmWaz0AAE2lT0316ak+Rb13zu/tmsXX2BHTDm+zT1N3331322OPPez888+PlmxuzZo1QS/xwIEDg1k1dtllF/vLX/4SrTU7/vjj7aijjoruWbBOPdoqaGOXXXZZ8LXUy73//vsHy4YPHx48btiwcBinTJo0yUaPHh1cZVJf7+KLLw6mboypl/6rX/2q3XPPPTZixAjr2rVr8Hwff/xx9Ii66WTEE044wfr372+vvfaaffe737UvfOELtt1229khhxxid911l918883Ro83eeOONYHnfvn2DOcB33HFHe+CBB6K1ofi1JKrdm6/ZR+Jtp2k39fW/973vBetk2rRpduCBB1q3bt2CbavXc99990Vr0RoovDPA3PK5VlRZYKM+7RctCan4LqzMD9YDAHKbpppdUrEk6cwpm2PnfPgT6z6nMPgU9b/Hnh/c6r6Wa31d/66+NHeq25tuusmefvpp+8c//hEt2UgF6+GHH25vvfWWPfzww/buu+/aGWecERSPzz77bPCYAw44wJ5//vlgzmyZOnVqMLf5c889F9wXtfW4wYMH26OPPhose+WVV2z+/PlBz7NMnjzZTjrpJDvuuOOCr6PXpd54DQNJ9Oabb9rjjz8eFPgvvvhiMK/1ySc3fEVp9Xa//fbbdtFFFwWXXq+LCuaY5v7+9re/Hbxu/btx48YFhbu+z6a49dZb7c9//rP94Q9/COZ9f/LJJ+3LX/5ytNbsmGOOCYa7qOdd37OK/+7du0dr0RqYh6htNesCOvoIUL0R91904ibF9/QR8+34qyfZFUOvsEN7HhotzS1cmCT12KapxzZNvVzfprqK4OzZs23o0KGfz62s4vfr074etJuivr8tTfX3Hf5uvQp7Rfcapx5b9cyql/bYY48NCj8VtbrQjIpQ9SqrN/qggw4KhmIkFoQqdBctWmRPPfWU6W/qkCFDgp7kL37xizZq1Cg77bTTgt7q5cuXW1VVVfBvVdyrd1dfT73Un3766Sa93fvss09wAZi4MJdbbrnFfvaznwXPo2Exes1PPPFEMB473u7qIVZhrH3SsWPHYFksvoCO3qvf//737fXXXw/GdMc09OSTTz4J2tqX7733XtCui3rA1XN99913B/cTt1+s9vem3m5tVx3UJBb2sZKSErvtttuCoj4Zdb3vYlxAJzn0eGeAMd3HWPe8bjb+7MnBL0TR7YSznrYeed2D9QAANFVBRYc6P00tqGzb8uCaa64JemR///vfR0tC6o3WUA/1VGtYR5z7778/eLxona6kqN5h9WBruQpzDatQT656pVXMf+UrXwkeXx8VvSq+E+23337BQV7iUBIV9olFpzrS1DOvAwHZfvvtP3+dumKl1DfvtXqg1Zuv16tLy8eWLVtm55xzTjDsRQcNei718KvobQr14Ku3feTIkcEnBeql1/cT0xCfH/zgB0GxfsUVVwRDXNC6KLwzQKcOneymkTfbskEVQS/EHg9cF9wuH1RpN468KVgPAMht3Qu6B73OyeaCQRdYZWG1zRi+IHqGkDp2Kguqg/V1/bv6oq/fXCqezz333GBYh8Z0x9RjrDHOKk4TM336dHvmmWeiR4XDTVR4K7vttlvQkztmzJjgvnp7NbxCPdaNqd0rHBfMictrDxWJ18VDXdSzHr/OuJDV2GnR606knvqtt97aevToES0JqTf773//u/3yl78MDhz0XF//+tetvLw8eoQXcH4wUbugT1wvO++8c9D7rRM31Rv/4x//OBgjv3LlymC9tvcHH3xgRx99dDDee88997RLLrkkWIfWQeGdIXbosoM9udNTNq7PuOAXYmFeoT2241+D5QAA5OflB0M9ks03en2jzk9Tx0efpmp9Xf+uvujrt4SGhoh6v2MqolWIr169OihQE6PhDjEV3v/617/sb3/7m33ta1/7fJnGe6v4VjumEydFQ1ASqadaPeSJ/v3vfwdDmbbaaqtoSeP0uhJfp+jkSOXqq69O6nLrGsv9ox/9yL75zW8G/04ngqpATqRhHZqeMFFdPdbqLT/yyCODMesajjNjxoxgu8R0kujpp59uf/rTn4Je71//+tfRGrQGCu8Mop7tE/udGLQraipsQfmmvRQAACSr9qepez4wMbhd0U6fpqpAVOF34403Rkss6LVWIf2d73zHHnvsMZs5c2YwTlrF4W9/+9voURYMldiwYYM98sgjnxfZ+rd6rMY4qx1TEaveYp1MqbHjGr8tKvw1vvu6664LhqvopMQJEyYEPfHJ9JY3RL3imjFFQ2E0xlsnimpoi8Z3a5y6esnz8zceuKiH/KGHHgpOrFRPtHrA46EsMY1917+/4YYbgtsHH3zQ7rjjjmhtSD3dmg1FX0s937/73e+C711DT3RAc+aZZwZFuNZpfL0OXDSUBq2HwjvDbFmwpQ0qDk/MnLEuvKAOAADNEX+aqpP0Txr0A7uo90X2+A5PtNunqRpvrPHaMRWsOplRhfd5550XFKSHHnpoMDY6sRdaM3NoWIWGgey1117Bsn79+gUXqtF4bE0lGNPYbxXXKkp1sqKGXohOXtQYc53Uqd5vnZiokzQ153Yq6PVpvLWmAPzFL34RnAiqKQU104nGn6tnPqaLB6lA1hAZnUCsQlm934l0MHHttdcGM5HssMMOQaGe+GmBaIy5DmT0PHqMDl50QKHvr6CgIDjo0DZXsa2hLJq+UAcFaD3MatK2mjWrSW0XzbzInl3xrB3d+2i7cPCF0dLcxGwRqcc2TT22aerl+jZtaHaJ5opn4FCxpqIPLZdt25RZTVqOn6wMNLrL6OB2+rpNT9IAAABA+qLwzkCjOofjrz5c96FV1my8ohYAAADSF4V3BtquczgtUVlNmX26/tOgDQAAgPRG4Z2BtsjfwoYUDwnaDDcBAADIDBTeGSoebkLhDQAAkBkovDPU6M7hCZZMKQgAuSe+WmJ9lyIHWkNdV/JE01B4Z6i48P5o/UdWUd34VbAAANkjvmy5pncD2kr8fqt92Xwkj0OWtpWSebxlbdVa2+/t/azG/3tguwc+P+Ey1zA/cuqxTVOPbZp6ub5N1fM4b968YI7oXr16BXMqt7QXUnNOr1271rp06cI83imSLdtU7zcV3UuWLAnmJB8wYMBm7zfm8U4OhXfbSlnhLUe9d5TNKptlPx/yc/tWr29FS3MLBU3qsU1Tj22aemzTTYvvVNDzrV+/Prg8OkMJUiPbtml9RbdQeCeHn6y2ldLC+9JZl9rTy562o3odZZcMuSRamlv445t6bNPUY5umHts0pMKuqqrKKioqgnZLlJeX20svvRRccr2oqChaipbIlm2qQlvDS/Lz8+s9gKDwTg6Fd9tKaeH94KIH7YY5NwQznNy/3f3R0tzCH9/UY5umHts09dimqcc2Tb1c2qYU3slhEFcGi6cU1AmW5dXlQRsAAADpicI7g23baVvfgR2Cy8Z/vP7jaCkAAADSEYV3Buuc39mGdxwetJnPGwAAIL1ReGe4eBpBrmAJAACQ3ii8MxxXsAQAAMgMFN4ZLi68Nca7rLosaAMAACD9UHhnuG06b2P5/l+V/6fZTQAAAJCeKLwzXKcOnT4/wZJx3gAAAOmLwjsLjO7COG8AAIB0R+HdPAd4qjyzgnvtLL6QzvS19HgDAACkKwrvphvgmeR5JriXBuITLD/d8Kmtr14ftAEAAJBeKLybJt/zkOcWz/+0IB1s0ynhBMt1nGAJAACQjvKiWyTnWs8OnkM94z0neoZ56tPF0zlsBnp6ZsycOdMGDFDHeeqM+2ScfbThIzu337n2nZ7fiZZmv/Lycps6daqNGTPGioqKoqVoCbZp6rFNU49tmnps09TLpW26ePFiGzx4sJp9dFcNbI7CO3kqtn/r2cWjN9QET2OFtx6jAn0TkyZNspKSkuheajzZ/Ul7s+ubtuPaHe3IZUdGSwEAAFrfypUrbdy4cWpSeDeAwjs5gzxveL7reUELXDKFd5v1eD+27DGbOH+iDS8ebg9s/UC0NPvRQ5N6bNPUY5umHts09dimqZdL25Qe7+RQeCdHXciPeTSTSUzj47X9tOw0z+89jentWVRaWmqDBqmWTx3NaHL8B8f7i+pgL+70onXOT6z3s1dZWZlNmTLFxo4da8XFxdFStATbNPXYpqnHNk09tmnq5dI2VeHdp49qbgrvhnByZXKe82hs984JudMzL2qrKG9XW3fa2gryCqza//tg/QfRUgAAAKQLCu/krPa8WyuLPBVRe7mnXRV1KApmNxEupAMAAJB+KLyzSHwhHQpvAACA9EPh3Xw6ubKhEyvb3NYdtw5uX1r5kk1eOpmL6QAAAKQRCu8sMW3tNLtzzh1Be82GVXbZ7Mvs8LcPC5YDAACg/VF4ZwH1bJ/74TnWa25Hu+/icfa/4y4MbnvMKQyW0/MNAADQ/ii8s8DU5VNtec0Ku+K2Q230zP7BMt1OuP2QYLnWAwAAoH1ReGeBueVzraiywEZ92i9aElLxXViZH6wHAABA+6LwzgIDiwZaeUGlzRi+IFoSmj5ivlUUVAXrAQAA0L4ovLPAmO5jrHteNxt/9uSg2Bbdjj/raeuR1z1YDwAAgPZF4Z0FOnXoZDeNvNmWDaqw46+eZHs8cF1wu3DgWrtx5E3BegAAALQvCu8ssUOXHezJnZ6yK4ZeYdtuGV5IZ5ctdw2WAwAAoP1ReGcR9Wwf2vNQO7HvicH9d9a+Y9U11UEbAAAA7YvCOwvtssUuwe2qqlX28fqPgzYAAADaF4V3Fupe0N226rhV0H5tzWvBLQAAANoXhXeW+uIWXwxu31j9RnALAACA9kXhnaW+2DUqvNe8wThvAACANEDhnaV27bprcLuyaqV9suGToA0AAID2Q+GdpXoU9rARHUcE7ddXvx7cAgAAoP1QeGexxOEmAAAAaF8U3lls1y3C4Savr3ndampqgjYAAADaB4V3Fot7vFdUrrCZG2YGbQAAALQPCu8s1rOwpw0rHha01esNAACA9kPhneXi+bw5wRIAAKB9UXhnucQTLBnnDQAA0H4ovLNcfILlssplNmvDrKANAACAtkfhneV6F/a2ocVDg/Zra14LbgEAAND2KLxzAOO8AQAA2h+Fdw6ILx/POG8AAID2Q+GdA+ITLJdWLrVZZYzzBgAAaA8U3jmgT1EfG1w8OGi/sZrLxwMAALQHCu8cEfd6cyEdAACA9kHhnSM+P8HSC2/GeQMAALQ9Cu8cEfd4L6lYYp+VfRa0AQAA0HYovHNE36K+Nqh4UNBmuAkAAEDbo/DOIZ+P82Y+bwAAgDZH4Z1D4sKb+bwBAADaHoV3Dtl1i/BCOosqFtn1c663yUsn2/rq9cEyAAAAtC4K7xyiEysLqvOD9qPz/2SXzb7MDn/7MJu2dlqwDAAAAK2HwjtHqGf73A/PscHzu9l9F4+zl4+9ILjtMacwWE7PNwAAQOui8M4RU5dPteU1K+zK2w630TP7B8t0O+H2Q4LlWg8AAIDWQ+GdI+aWz7WiygIb9Wm/aElIxXdhZX6wHgAAAK2Hwjt53/a84lnm2eD52HOVp8iT9gYWDbTygkqbMXxBtCQ0fcR8qyioCtYDAACg9VB4J08F93WefTzbeS7wnOa53pP2xnQfY93zutn4sycHxbbodsJZT1uPvO7BegAAALQeCu/kaRD0nz3veWZ5HvM84Nnfk/Y6dehkN4282ZYNqrDjr55kezxwXXC7dFC53TjypmA9AAAAWg+Fd/ON8hzseT64lwF26LKDPbnTU3b50MutuCgstE8b+MNgOQAAAFpXXnSL5K3xFHo0tvu3njM9VZ66dPF0DpuBnp4ZM2fOtAEDBoRL2snEeRPtseWP2V5d97Lrh2bEaJk6lZeX29SpU23MmDFWVJQRw+3THts09dimqcc2TT22aerl0jZdvHixDR48WM0+uqsGNkfh3XRbezp6dP31az0qvsd76jLBs9m6SZMmWUlJSXSvfXzc8WN7sPeDll+TbxfMvcCKavglCwAAmmflypU2btw4NSm8G0Dh3TIneH7vURW9VgtqSdse7/Lqcjvkg0NsXfU6u3rw1bbflvtFazILPTSpxzZNPbZp6rFNU49tmnq5tE3p8U4OhXfL6NDuHo8K6uVa0IjenkWlpaU2aNCgcEk7unDmhfbciufsiJ5H2Pih9XXap7eysjKbMmWKjR071oqLi6OlaAm2aeqxTVOPbZp6bNPUy6VtqsK7Tx/V3BTeDeHkyuRd6jnIo6Em23iO9WioyeOeZIrutLNvyb7B7b9W/suqauobpg4AAIBUoPBOnqYBuc3zjud1z888N3mO8WSkvUv29jdAB1teudzeXftutBQAAACtgcI7eZd4tvVozPaWnh096vFe78lI3Qu6205ddwraL658MbgFAABA66DwznHxcJN/rvxncAsAAIDWQeGd4/YrCWcz+XTDp/bZhs+CNgAAAFKPwjvHDe041IYWDw3a9HoDAAC0HgpvfN7rTeENAADQeii8Yft1Cwvvt9a8ZSsqVwRtAAAApBaFN2yHLjtYt4JuVuX/vbTypWgpAAAAUonCG5afl2/7bLlP0GZaQQAAgNZB4Y3Avt3CaQVfWvWSlVeXB20AAACkDoU3AntusacV5RXZuup19tqa16KlAAAASBUKbwQ653e2L23xpaD9zxXMbgIAAJBq2V54V3n6hM1N9PRoHRIkTitYU1MTtAEAAJAa2V5450W3tRV4KsMmYvuUhCdYLqxYaB+s/yBoAwAAIDWytfDWmYLh2YJmX/bE95X9Ped6ZnuQoE9RH9uu03ZB+8Y5N9rkpZNtffX64D4AAABaJlsL7xc8z3s0XuKx6H6c5zxnesZ7kGDa2mn22brweOSdFW/ZZbMvs8PfPixYDgAAgJbJ1sJ7uGcrj4aa7B7djzPQU+J5yIOIerbP/fAc6zevq9138Th7+dgLgtsecwqD5fR8AwAAtEy2Ft7qtp3l0fenufF0P858T7UHCaYun2rLa1bYFbcdaqNn9g+W6XbC7YcEy7UeAAAAzZethXfsQM/eYTNwmudNzwMe9XojMrd8rhVVFtioT/tFS0Iqvgsr84P1AAAAaL5sL7yv82jqQBnpud3zsmcnz0QPIgOLBlp5QaXNGL4gWhKaPmK+VRRUBesBAADQfNleeGuc93th047wPOs5w3Oq5xAPImO6j7Hued1s/NmTg2JbdHvZWZOtR173YD0AAACaL9sL70SqHP8eNk3jJuKecLhOHTrZTSNvtmWDKuz4qyfZHg9MDG5L+y+367e5IVgPAACA5sv2wlvz4GnqwK94NH/3Pzwy2LM4bCK2Q5cd7MmdnrIrhl5hxww81t8cHayyQ5WtrVobPQIAAADNle2F9yWeEz0veu7zTPfI4Z5XwyYSqWf70J6H2k8G/uTzS8g/uezJ4BYAAADNl+2Ftwru3h4NK9GMJrG7PeoJRwMO76njE7MXVrxgqypXBW0AAAA0Ty6M8a7y6Oov23tGe4o9n3g2nb4Dm9mrZC/rXtDdymvK7Znlz0RLAQAA0BzZXnjr+/ulZ7nnHc+7UfsqTy4cdLRIYV6hHdIjnPzlyaUMNwEAAGiJbC8+VWBr+sCLPLtG+ZnndM8VHjTi8B7hcJN3171rn67/NGgDAACg6bK98D7Bo7Hdt3re9rzluc2jwnucB43YpvM2tm2nbYM2J1kCAAA0X7YX3jqxUsV2bVqmdUjCYT0PC26fXva0VdVoyDwAAACaKtsL7489R4XNTXzToxMskYSDux9s+f7f4orF9r9V/4uWAgAAoCmyvfC+zvMrz588P/ac7fmj52rPNR4koXthd9unZJ+gzXATAACA5sn2wnuS5whPX894zwTPAI/OGNQFdZAk5vQGAABomWwtvLt51MNd4pns2deji+goKsS38WzpQZL2Ltn78zm9n13+bLQUAAAAycrWwltTCH7NszK4t6kVnjEezWyCJGlO74N7HBy0GW4CAADQdNlaeOuEyjvDZp3u8nw3bCJZ8Zze09ZOs083MKc3AABAU2Rr4b2VR1eprM97Hj0GTTCy88jP5/R+aulTwS0AAACSk62Fd4GnKGzWqbunMGyiKeI5vR9d8qjdOe9Om7x0sq2vXh8sAwAAQP2ytfD+0POVsFmnPT16DJpoWPEw61CVZ6urVtu9c35vl82+zA5/+7Bg+AkAAADql62F9yOeyz0jgnubGu7RtIIPB/eQNPVsXzbzUhu6oIfdd/E4e/nYC4LbHnMK7dwPz6HnGwAAoAHZWnjf7Fno0ThvnUh5TpS7PRrfXerRY9AEU5dPteU1K+zK2w630TP7B8t0O+H2Q4LlWg8AAIC6ZWvhXebZ3/Nbjy4Pf2OUI6NlWqfHJOskz/OexZ7Vntc9x3pyytzyuVZUWWCjPu0XLQmp+C6szA/WAwAAoG7ZWnjLWs+5nt4eXblS1aLaWrbG0xQHeJ7wHOLZxaPLzuvKl0d7csbAooFWXlBpM4YviJaEpo+YbxUFVcF6AAAA1C2bC+9E6qleFDab5TjPTZ5XPR97JnqmeHJqLvAx3cdY97xuNv7syUGxLbq97KzJ1iOve7AeAAAAdcuLbtF0L3k0lccPg3t16+LpHDYDumT9jJkzZ9qAAQPCJRnmvXXv2QWzLrAVNSusoDLfKguqrKA6324f8WvbscuO0aPaTnl5uU2dOtXGjBljRUUNzSCJZLFNU49tmnps09Rjm6ZeLm3TxYsX2+DBg9Xso7tqYHMU3s0zzqMTNffwvKkF9dDsKePD5kaTJk2ykpKS6F7mqcirsBmdZtjCwoX2v67/s+oO1Xbk0iNtx3VtX3gDAID2t3LlShs3TuURhXdDKLyb7giPpiL8keceLWhA1vV413b9vOvtL8v/YiOKR9h9W91neXlt+5aihyb12KapxzZNPbZp6rFNUy+Xtik93smh8G6a73nu9Zwe3TaVTu5cVFpaaoMGDQqXZLjSslL71nvfsmr/79atbrW9S/aO1rSNsrIymzJlio0dO9aKi4ujpWgJtmnqsU1Tj22aemzT1MulbarCu08f1dwU3g3JlZMrU+FUj4rtE6NbuMHFg+2Abpr0xWzSwknBLQAAADZH4Z0cTUF4h+cnnhc8mpownp4w553Q94Tg9vU1r9u7a3XNIgAAANRG4Z0cFdz5njs9mkcvjqYXzHmju4y2L23xpaD9h4V/CG4BAACwKQrv5AzzaDx87Wg5XNzrPXXFVPtsw2dBGwAAABtReCMlvrzFl22bTttYjf93/6L7o6UAAACIUXgjJTSNYNzr/eTSJ21pxdKgDQAAgBCFN1LmwO4HWr+iflZeU24PL9ZU5wAAAIhReCNlCvMK7dg+xwbtBxY+YL+e+2ubvHSyra9eHywDAADIZRTeSKltOm5j+dUdbEPNBrtv7iS7bPZldvjbh9m0tdOiRwAAAOQmCm+kjHq2L/7kIhsyv4fdd/E4e/nYC4LbHnMK7dwPz6HnGwAA5DQKb6TM1OVTbXnNCrvytsNs9Mz+wTLdTrj9kGC51gMAAOQqCm+kzNzyuVZUWWCjPtVFPTdS8V1YmR+sBwAAyFUU3kiZgUUDrbyg0mYMXxAtCU0fMd8qCqqC9QAAALmKwhspM6b7GOue183Gnz05KLZFt5ee+ZR1yysJ1gMAAOQqCm+kTKcOneymkTfbskEVdvzVk2zPByYGt5/1X2ajt9g+WA8AAJCrKLyRUjt02cGe3Okpu2LoFXby4FPs6F5HW3V+jb206iWbvnZ69CgAAIDcQ+GNlFPP9qE9D7XT+p9mFwy+wHbtumuw/Po511tNTU3QBgAAyDUU3mhVeXl5dt6g8yzP/3t77dv2jxX/iNYAAADkFgpvtLpRnUfZYT0OC9q3zr3VyqrLgjYAAEAuofBGmzhzwJnBEJR55fPsoUUPRUsBAAByB4U32kTvot42ru+4oP37Bb+3pRVLgzYAAECuoPBGmzmu73HWt7Cvra1ea3fMuyNaCgAAkBsovNFmNNTk7IFnB+2/Lv2r/W7+7+yu+XfZ5KWTbX31+mA5AABAtqLwRpsa232sjSgeYXlVZr+Z/xu7p/R3dtnsy+zwtw+zaWunRY8CAADIPhTeaFMbajbYkrLFNmR+D7vv4nH232PPD257zCm0cz88h55vAACQtSi80aamLp9qq2y1XXX74TZ6Zv9gmW4n3H6ILa9ZEawHAADIRhTeaFNzy+daUWWBjfq0X7QkpOK7sDI/WA8AAJCNKLzRpgYWDbTygkqbMXxBtCQ0fcR8qyioCtYDAABkIwpvtKkx3cdY97xuNv7syUGxLbq99KynfHlJsB4AACAbUXijTWlKwZtG3mzLBlXY8VdPsj0fmBjcftZvmY3osrV1zOsYPRIAACC7UHijze3QZQd7cqen7IqhV9jJg0+xY3ofY9X5Nfb6mtftqWVPRY8CAADILhTeaBfq+T6056F2Wv/T7KeDf2rf7PnNYPnE0ok2t4wTLAEAQPah8EZaOG/QeTa4eHBwOfnLZl1mVTVV0RoAAIDsQOGNtNA5v7NdOexKy/f/3lr7lv1h4R+iNQAAANmBwhtpQ2O/T+53ctC+Y94ddvf8u+2u+XfZ5KWTuaIlAADIeBTeSCs/6P8DG1483Gqqqu3O+XfaPaW/s8tmX2aHv32YTVs7LXoUAABA5qHwRlqprKm0ZeVLbcj8HnbfxePsv8eeH9z2mFNo5354Dj3fAAAgY1F4I61MXT7VVtassqtuPzy4jLzodsLth9jymhXBegAAgExE4Y20Mrd8rhVVFtioT/tFS0Iqvgsr84P1AAAAmYjCG2llYNFAKy+otBnDF0RLQrqsfEVBVbAeAAAgE1F4I62M6T7Guud1s/FnTw6KbdHtpWc+ZYXVBbbXlnsFywAAADINhTfSiq5oedPIm23ZoAo7/upJtucDE4Pbz/ovs4oOlXblZ1dycR0AAJCRKLyTt6/ncc9sT41nggetQPN5P7nTU3bF0Cvs5MGnBLfjR4Sb+8WVL9otc28J2gAAAJmEwjt5XT3TPRd6Nh2AjJRTz/ehPQ+10/qfFtwe1vMwO2vAWcG6BxY9YI8sfiRoAwAAZAoK7+Q97bnY87CnTAvQtk7se6J9o+c3gvbE0onB1IJ/W/E3e3HLF4Nb5vgGAADpjMIbGSMvL88uHnKx7b7F7lbt/130yc/syrlX2kud/x3ccnVLAACQzvKiWzTNLM+9nsbGeXfxdA6bgZ6eGTNnzrQBAwaES9BkiysW21Hvf8sGzetmV/76sGCOb818Mv6sp4OTMv+y3V+sY4eO0aPRVOXl5TZ16lQbM2aMFRUVRUvREmzT1GObph7bNPVyaZsuXrzYBg8erGYf3VUDm6Pwbp5kC2+tHx82N5o0aZKVlJRE99BU73R+x/7a8692/0UnbnKhHRXfmgHlyKVH2o7rdoyWAgCA1rZy5UobN26cmhTeDaDwbh56vNvR7xb9zu6bP8n+e+z50ZKNNP3gCf1PtB/0+UG0BE1Fr1fqsU1Tj22aemzT1KPHG7VReDdPsoV3bb09i0pLS23QoEHhEjTZ5KWT7bLZl9Xb463pBzUTCpqnrKzMpkyZYmPHjrXi4uJoKVqCbZp6bNPUY5umXi5tUxXeffqo5qbwbggnVyZP0wnuHEWHrar41B7tQRtq6OqW+dUdbFjHYcEyAACAdELhnbzdPG9G6e/5YdTWNINoQ/Vd3XLOgOVW1aHazv74bHtv7XvRowEAANIDhXfyXvBoaE7t0L3aDuKrW1468FLba+1XgtsHt/+jDSgaYCurVtoZH51h/13532BYyl3z7wpumecbAAC0JwpvZCz1fB/c7WDbb9V+we1Wnbay3438nQ0rHmZrq9fajz86OxgLfk/p74Jb5vkGAADticIbWaVPUR+7fevbrbC6wIbM62H3XTwumP1Etz3mFNq5H55DzzcAAGgXFN7IOm+secMqOlTaVbcfHlxcR3Q74fZDbHnNiuBS8wAAAG2NwhtZZ275XCuqLNhkqkFR8V1YmR+sBwAAaGsU3sg6A4sGWnlBpc0YviBaEtKUgxUFVZbv/wEAALQ1Cm9knYbm+e5QlWf/N///7ImlTwTLAQAA2gqFN7JOffN8rxhcYQM7D7Jy/+/y2ZfblbOvtBWVK5hyEAAAtAkKb2SleJ5vXT7+5MGnBLdP7fS0/XHUH+2InkcEj/nr0r/a2Le+zpSDAACgTVB4I2up5/vQnofaaf1PC251v2OHjjZ+6Hi7eNDF1qE6zwbOK2HKQQAA0CYovJGTOuV3suoONUw5CAAA2gyFN3JSQ1MOFlTmW2lZabQEAAAgNSi8kZMamnKwsqDKnlr6lE1fOz0YcsLJlwAAIBUovJGT6pty8LKzJgdTDs6vmG8nfHCCfe3NAzj5EgAApASFN3JSfVMOrhxUZVdudZXt2HlHy6sy6ze3KydfAgCAlKDwRs6qa8rBJ3Z60sb2GGtH9T7KqvM5+RIAAKQOhTdyWl1TDsq88nkNnHzZwV5Z/Uq0xBgHDgAAkkLhDdSh4ZMvq+2pZU/ZaR+eZn9e/Odg3DfjwAEAQGMovIE61Hfy5fiznrbimqLg/utrXrdrZ11j3ebkMw4cAAA0isIbqEN9J1+uGFRpv93uLvvDtn+wbTpuE4wDv/I2xoEDAIDGUXgD9ajv5Est377L9kGveEMX4flg/QfBfcaAAwAAofAGGlDfyZfS2EV4Hlz0oJ3ywSl28FtjGQMOAAAovIHmaugiPAXV+Vbj/7296i3rOaeYMeAAAIDCG2iu+saB6yI8d2/3f3ZCnxManAv8L4v/EiwThqMAAJD9KLyBFqhvHPiOXXe0TvmdGpwL/Ma5N9qJH5xoE0sn2qFvH8JwFAAAshyFN9BC9Y0Db2wucFFx/ciCh617aeEmw1G61xqOQo84AACZj8IbaCX1jQGfcNbT1iOvu/1u5O/sK1t+pc7hKJdHw1HUG/7Cihe4SA8AAFmAwhtoJfWNAV8+qNJuHHmT7dx152BawoaGozy+9HG74KPzrSSJi/TQKw4AQHqj8AZaUUNzgUtjw1E65nUMe8TruUjP1bOvtulrp9ubq99stFecwhwAgPZF4Q20svrGgEtjw1GO6XNMgz3ik5dPtuM/ON5O++DUzS5dnzhOXAV4MsNV9Ni/rfibvbjli8EtxTkAAKlD4Q20o8aGowzrOKzBHvG+hX0tz/+r7rD5pevjceJnfHiG/ej9M7wQL2hwuEpcnF8590p7qfO/g1t6zQEASB0Kb6CdNTQcpbEe8Ue3f9RO7HuiFVbm19srPm3dNFuXt96uuO2wTQrzeLjKTXNusjdWv2HnfPiToBhXUf7ycRdsVpw3pdec4hwAgM1ReANpoL7hKI31iGv98I7DraKgqt5e8QFFAxoszB9d8qid+tGptqJmpV1+26F1FueTFkzapDBvrNe8oeI8mcKc4h0AkI0ovIE019gJmo31ip/c7+QGC/Pu+d2D4SoFFR3qLc7vXnB3g4X5rXNutZdXvWw/+fDHwdjy+orzZArztirek3kOAABSKS+6Rdvo7VlUWlpqgwYNCpegRcrKymzKlCk2duxYKy4ujpbmHhWlKm5VBKt3W4W2im71im/daeugcFUBrEJZBXNcmKvnXEX8s8uetcs/u9zuv+jETYpvPU497KLC/H/HXRi0E+3xwHWfXxBI6nuO3bvubu+unmb95nWxy28PC3itG++vY9mgcntixyetQ16Hel+rev11APLx+o8//1514qnGwOvAQ58MxAcjidujrsc0tj6mYnzq8qk2t3xuMAONDnLiTyOksfWix/x98d/tn9P/afuO3te+3vvrzXqOVLyOdHiOVOFnP/XYpqmXS9t08eLF1qdPHzX1v8VqYHMU3m2LwjvF+EOxUUNFT0OFuQpN/duGivPTBv3Qrim9pt6iWr3mq6pWWV5FTaPFeUPFfVFekZXXlNf7mH223MdeX/Wa9Ztbu3ifbEsGldnd2/2fFXcothOnq6e9KDjBtHbx/qcd/mzfmfbtBot7bbdUFO+bbPcK3+6FVS16jlS8jvZ8DklF8a7HcDCzUaqeg226UaqeI1u2aTIovJND4d22KLxTjMI7eY39ck0snJrTa/7c8uds/Ozx9RbNmoFl2fol9vKxF0RrNkoszJPpWW+oeI/V95iS/BJbWbWy3vV7bLGHbdVxK3ts0V+s37yudkWtAn/hwHX2k0Hn2K1zbrE+czttfgAwsMxuHnmLFeQV2NkfnmU95xRvcgAQ9/DfN/p+O3768cE2resAQQcA0tAnAMkcRKTLc+i9luoDAA5mUv8cbFO2ae3nSBaFd3IovNsWhXeKUXinlorz+npoEn9BN6fX/EeDz7SrPruq3oL3jP5n2IbqDXbPwnvqfUyvgl62csPyFhfvLR02E2tsvTT3OdT738H/21Czod7HdO3Q1dZUr6l3vQ6w9Byl5aX1PqZPYR9bVLGo3vWjOo2yvLw8m75uer2PGVw0uMGvsWvXXYOTfP+x9FnrP2+LTQ5mLvODlfkDVtuBPQ+yZ5c+U+/6I/p8I3jOJxc9sdkBkR6zYMAaO6z34fbU4ifrXf+tvkcFz/GXhY/W+5hv9vlWnQdd8frv9jva/3Dm2cML/ljvY77d9zv254V/qnf99/sfE7yOh+Y/WO9jvtPvu/anBY/UuX7hgLWbPEffeV3qfIxe6yMLHq53/TH9jw2e48H5D9T7mKP7fS/4Xutbf1z/44PnuH/+ffU+5nv9vm9/XPAQz5Ehz1H7gDlZFN7JofBuWxTeKUbhnXoNbdPW7jWXlhTvlwy+xCpqKmzinIkNFoFvrHmj3vU7d9nZVletts/WzKq3wM8vLLTq8spg2sXaUn0AkE3PIfVt91hj64Xn4DliPEfrPYdO5NcsW8mi8E4Os5oASJqK7Pquwinq+VYvSV0zsOix+vhSPSn6pV7X1IiNPWZsj7HBx6Dj65nB5ZCeh9g3en2jwcdcN+K6Btffvs3tNq7vuOCgob6ZYA7tcWjw0XF96y8fermNHzreKgur633MQd0OanD9SX1PCuZob+gxOmG1ofWH9TgseK0NPWb7zts3uH7fkn2DNPSYrTtu3ejrHN15dHAwlvgHXnRwpZlzehb0bHD9Np22CdLQY3oU9GhwvYYQKQ09pltBtwbXDyseFqShx2g4U0PrhxYPDdKS5xhSPCRIaz/HlvlbNrh+cPHgIDxHdj2HlquDBalH4Z28Qzxveco8szznewDU0lBx3lBhHmvt4r17QfdGn0M9+Q0V5z8a+KMG1x/Q/QA7sPuBDT7mwiEXNrj+B/1/YKf0P6XBx/xqxK8aXH/RkIvs4iEXN/iYW7a+peGvMfxXQRp6zJ0j72xw/Y1b32jf6/29Bg9m9ivZr8H1x/c5PkhDj/lqyVcbXK8DqsYOqsaUjGlwvabnbGyKzgO6HdDg+h/0+0GQljzHKf1OCdLaz/G1bl9rcP2p/U4NwnNk13NouT7VROox1CQ5u3n+67nJc49nD8+dHhXft3uSxVCTFGOoSeplwjZN5iz8xh7TkmEzOghobD3Psen6xs4BeGSHP9V7gqbWJzMUiefgOXiO1D5H7d+rDWGoSXIovJPzoGdrz+7BvdBEz/c8Qzw1WpAECu8Uo/BOPbbpRi0t3kWPYZq2UDocAPAcPAfP0bTnSBaFd3IovJMz26MzES4L7oUO8PzDM9yjoSd16eLpHDYDPT0zZs6caQMGDAiXoEXKy8tt6tSpNmbMGCsqKoqWoiXYpqnHNt1IM9c8v+p5m1c+L5jpZP8t97eOHTpGaxtfL3rMP5b9w/7z4X9s75F729d6fK1Zz5GK15FNz8E23ShVz5Et2zQZKrwHDx6sJoV3Ayi8k1PuOcfzm+BeaHvPu569PBqGUpcJnvFhc6NJkyZZSUlJdA8AACCzrVy50saNG6cmhXcDKLyT01Dh/WXPy1pQB3q8Wxk9ianHNk09tmnqsU1Tj22aerm0TenxTg6Fd3LqGmoyxvOcp6GhJrUxxjvFGI+cemzT1GObph7bNPXYpqmXS9uUMd7JYTrB5PzH8/Ww+bmxnjkeFeUAAABAgyi8k6NpBHf1XOPZznOC52zPtZ5kZzQBAABADqPwTs6rniM9B3ve9lzp0bCTpszhDQAAgBxG4Z28yZ6dPBqkNdSjebwBAACApFB4AwAAAG2AwhsAAABoA0wn2LY0xc7Ct956i3m8U0RTNcVzpDL9VWqwTVOPbZp6bNPUY5umXi5t0yVLltjo0aPV7OtZpAY2R+HdtnbxvBE2AQAAss4oz/thE7VReLctHQUu8OiQcIkWoMWCq4F69IO+VAvQYmzT1GObph7bNPXYpqmXS9tUNWUPz0eeKi0A2puuXKl5v3WL1GCbph7bNPXYpqnHNk09tmnqsU2xCU6uBAAAANoAhTcAAADQBvKjW7QNjX+q8LwY3aLl2KapxzZNPbZp6rFNU49tmnpsUwAAAAAAAAAAAAAAAAAAAAAAAAAAAACt6xDPW54yzyzP+R4kb1/P457ZHl2IYIKntpGev3vWeXR1sLs8XT2o20me5z2LPas9r3uO9STq73nEsyqK2rr6Kur2bc8rnmWeDZ6PPVd5ijwx3qfNd4BHV8PT79BEbNOm0e9P/R6tnWGe2B6elzx6H8/3XOsp8KB+ukLlHZ55nvhv/emeRGxXoA3s5tEUQtd5dMnYEz36oTvLg+TowOVqz9Ee/bKqXXjrj2ypZ4pnF88Yj37p/dmDut3vOdfzJc/Wngs81R5tY9Ec/yrGlT2jvOFRYcn8/3XT+07F9/YeFTHf9Czy3OoR3qfNN8Azx/M3T2LhzTZtOv3+1DbrVyvx9MKDPTrQnuT5gucbHh1MXu9B3fQ+nO551rOfRz//X/ao0yjGdgXayIMeFSuJJnr0i0/ze6Jp9Ee1duF9mkc9DN2De6FDPerFUVGJ5DzteTRs2kEebT8VkTG1texrwT0k4ybPtLDJ+7SZVBBqDmQdHOpnP7HwZps2Xe1tWNuvPOrgSLzOx5me9Z4tgnuo7XKPtmlxcK9ubFcE6LlqfXt71BuTSPcHeYYG99BS2sb/8ywP7oWe8agH9yvBPSSjm2dJ2Ay2qQ4O3wvuhdRWryPbNDn6hOtgj4b0CO/T5lHBstZTV88g27R51MOtn+/4U4S9PDFtU/XcalhPTH+zOnq+GNxDbUd5/uXRe1TF9QeeGz1dPDG2KwIU3q1P42QXhM3Pxfe1Di1X1zbW8B59jMc2Ts44j4ZF3Rncq3ubipaxTRu2xqNeWH30/IJHQ3qE92nTqfda5x7o/ale7NrYpk2nA5UTPNq23/foYFtF44EeqWubxvfZpnXbyvMdz5aewz36dEbDzn7vibFdEaDwbl91/SFBarGNG3eERwX3Dz1vakEj2KYN29mjHiydz3Gk5zJPY9imm9Ongvd4jvPoJOCmYpvWTT3cOlH6HY8K7uM9//GoWKxPvC3ZpnVTLaWDvVM8r3me8PzU811PH0992K45iMK79eljJ32slyieGaL20S+ap65tXOjp4WEbN+x7Hv0RPsOjIidW1zYVvXfZpg3TbCbvenQS1YWen3v0kTPv06bRJzC9Pf/wVEbRQYyG6Kl9sodtmhrqBY9nNalrm8b32aZ10zb70KNPW2L6HSDxkFK2KwIU3q1PPQlfD5ufG+vR2DpNj4eW0zbWNE0aoxzTx6Z6f2sd6naq516PemZ1m0jbTWfha5xyTG0t+3dwD8nQCdR6H2pKQd6nTfOcZwePPkGIo09mNF2b2o952KapsatHY75F200nUCfWB/qbpdm4NMsRNqdPDkZ4Ek+c1DSXEp/IynYF2oima9NR8DWe7TwaW6ezmJlOMHmaqin+w6s/uvrjq/Zoj8RTimlWDi3/qmemhynF6qdxx+o11PAS9brEUQ+j6I+D/hi86lFho+gjVKYTrN+lHs0Go9k0tvFobLJ6slQgCu/Tlqs9IwfbtOlu8OzvGe7RNtPc0zoZVWOTJZ72Tp+AaSYjDUXT/OhMe1e/nTw6r+PXnm09eh+qB/wBT4ztCrQhncTytkc/mOrlbmgsHTanX2IaA1c7iX+A9ctOsxnoIhoaa3e3R3+UUTdtu8a2qU74+ZNHF9hR1K79USk20uwbms1A70H9gdUY2p95OnlivE9bpq6p8NimTfOQR5+46u+R5pnXJwua/zyR5u2PL/Sig0cu9NI4XeBJHRXaZnqP6gAncVYTYbsCAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAQDa71/NC2Ew7umrpiWETAAAAyGwlnh5hs92KcF3KWpddr62fJ/HS9gCAVtAhugUAtK6VnmVhM6U6RrctscCzPmwCAAAAmS3u5VaPs4Z2JCbuhS7wqP2pR4Xwe54fehKp1/pqzy2exZ4PPPJ9z389y6M87/mSJ6avXfvrftUjaicONenv+aNnhUev40XP7p7YMI/+zfc8j3nWeko9P/MAAAAA7SouvLt6Hva85NEQD0XLRI95x3OQZ7jnux4Vv6d6Yiq8V3l+5dnWs71HTvYc4dnas4PnD56lnp4e0TAXFcfXe+KvW+SRxMI7z/M/j17HVzx6LhXh6rHv45G48P7Mc7xHX/OCaNkYDwAAANBuEsd11zXGW4V2tScupGOXed4NmwEV3v8Mmw3K96z2HBvcC9U3xjux8D7Ao/squGPFnvmeK4J7Gwvvi4J7G33iUW88AKAOjPEGgPSwmyfubV6TkEs823gS6TG1fcHziEfFr3rE1UPdxTPU0xQq/PVvpwX3QmUefc3aBwWvR7exOZ6+YRMAUBuFNwCkh/j38b6enROigrp2wasx1Yk6e/7h0YmWGpayp0f/VkNN4uEkTaHe7Np0UFB7eUV0G9N6/q4AQD34BQkAba/co6EgieLeY/VQf1xHGjLKo57mn3umeqZ7Kj29PInq+rq16YTObp7Rwb2QhproRE2tAwA0E4U3ALQ9DQdRsazebBXH6rFWcf17z52eEzw6YXEnz0keDTdpiMZub/Cc5dnKo5MiH/LUniJQX3dvj4p7fd1CT20q3F/x3O/RY/Ua7/Fo2MqvPQCAZqLwBoC2d5dHU//9y6MpAS/0yGmeGzwqtNVr/ZxHJz2qYG6IhpQc49H0gPp3ev5rPIs8idQjvoVHj9HXVWFdm4aLHOn50DPZ86pnoOdAT+3nAwAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAIA0YPb/oLVs0mJ01GEAAAAASUVORK5CYII=\" width=\"699.0476507935139\">"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "from cosmology import f,j,data"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "this = lambda args: f(data[:,0], *args) - data[:,1]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "import numpy as np"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 37,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "39.3 ms ± 53.3 µs per loop (mean ± std. dev. of 7 runs, 10 loops each)\n"
     ]
    }
   ],
   "source": [
    "%timeit least_squares(fun=this,  x0 = np.array([50/70, 0.5]), jac = jac, xtol= 1e-8, method='lm',gtol=0.01)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 36,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "36.4 ms ± 456 µs per loop (mean ± std. dev. of 7 runs, 10 loops each)\n"
     ]
    }
   ],
   "source": [
    "%timeit least_squares(fun=this,  x0 = np.array([50/70, 0.5]), xtol= 1e-8, method='lm',gtol=0.01)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "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.3"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}