Faza.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "%matplotlib notebook\n",
    "import numpy as np\n",
    "import matplotlib.pyplot as plt"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Posmatraćemo uticaj početnih faza na oblik signala. Signal sa kojim ćemo eksperimentisati je dat jednačinom\n",
    "\n",
    "$$\\tilde{x} \\left( n  \\right) = \\sum_{i=0}^{3} \\frac{1}{2i + 1} \\sin \\left(\\frac{2 \\pi}{64} \\left(2i+1 \\right) n + \\phi_i \\right).$$\n",
    "\n",
    "Ovaj signal se sastoji od prva četiri harmonika razvoja povorke pravougaonih impulsa u Furijeov red.\n",
    "\n",
    "Za početak neka su početne faze svih harmonika jednake nuli."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "phi0 = 0\n",
    "phi1 = 0\n",
    "phi2 = 0\n",
    "phi3 = 0\n",
    "n = np.arange(64)\n",
    "x = np.sin(2*np.pi/64 * n + phi0) + (1/3)*np.sin(2*np.pi/64 * 3 * n + phi1) + \\\n",
    "    (1/5)*np.sin(2*np.pi/64 * 5 * n + phi2) + (1/7)*np.sin(2*np.pi/64 * 7 * n + phi3)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "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,iVBORw0KGgoAAAANSUhEUgAAAyAAAAJYCAYAAACadoJwAAAgAElEQVR4nO3de9RsZ2Ee9qdIMhFQO7ZxsAJCR6Us4bWQUB2DkeLWxBeFVSzqIBO7qrAVC3MzXQZfKlkm9NBKYJsgRW5AcpfNOtTCSQ2iICtCF+ji0sh2HKkNEmrE5YgVRwog2+KgYxsJi+kfe389c+bMfGfP/mbPvPPu32+tvQ5nzzPzvWdrz8f7zOxLAgAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAMAWekqS/Ul+P8l/TDJJcmDJ13hukpuTfLVdbm7X7TXbxbcnuTDJuUm+22KxWCwWi8VS3HJumvnatweS7EtTOh5MU0KWLSDPTvKVJAeTvLFd7k/ycPtY32xXF7ZjtlgsFovFYrGUvVwYSPLEJE9v//eJaXaOA0s8//1JDic5dWrdqe269+0h29W5SSbXX3/95M4777RYLBaLxWKxFLZcf/31OwXk3J7zPSq2bAF5SpJHk7xnzmPvSfK1NrNsdhnfnWRy5513TgAAKM+dd965U0C+u8dcj8otW0DOafOvmfPYa9vHXtgjuwwFBACgYAoIu1m2gFzQ5s+f89hL28cu6JFd5JQce1LThVFAAACKpYCwm2ULyCva/HlzHjuvfeyiHtlF9mfBSU0KCABAmRQQduMbEAAAVkoBYTfOAQEAYKUUEHaz6qtgPZruV8Gazi5DAQEAKJgCwm52KyAnJXlOmsOgpt2Q5j4ez5hat3Nvjxv2kO1KAQEAKJgCwjyvT/KmJG9Os3Pc1f79TUnOajP7Mr+cnJHkUJq7m7+hXQ62687YQ7YrBQQAoGAKCPN8IQuuLpXk4jazL4u/HTkryS1JHmmXD+dIcdlLtgsFBACgYAoItVFAAAAKpoBQGwUEAKBgCgi1UUAAAAqmgFCb0RaQTz9waHLVbfdN9t94z+Sq2+6b3PvgoU0PCargvQWwWgoItRldAbn/ocOTl197x+S0S286Znn5tXdM7n/o8KaHCFvJewtgGAoItRlVATn40OHJ2W+5de4EaWc5+y23mijBkry3AIajgFCbURWQRZ/Ozvu0FujOewtgOAoItRlNAfn0A4c6TZB2FsetQzfeWwDDUkCozWgKyFW33bfUJOnq2+/b9JBhK3hvAQxLAaE2oykg+2+8Z6lJ0v4b79n0kGEreG8BDEsBoTajKSA+pYVheG8BDEsBoTajKSCOU4dheG8BDEsBoTajKSCTyfJX6nFDNTjaoveEq2ABDEcBoTajKiD3d7xXwSc/82U3VIMpx7vJ4Cc/82X3AQEYiAJCbUZVQCYTEylYVtebDCruAMNQQKjN6ArIjnsfPDR5/hW3T0679KbJ86+43aEksMCy74lF7y0A+lFAqE2VBaTruRsv+Y1PTE679KbJS37jE///85xMC0f0fU/Mvrfmva7zqwC6UUCoTVUF5HiHV80eAjI7SXI5UTha3/fEogKy7HsUAAWE+lRTQLoepz49wZmdJLmhGhyt73tiXgHp8x4FQAGhPtUUkD7nbvgGBHa3ym9AnF8F0I8CQm2qKCCrOk7dOSBwNO8tgM1TQKhNFQXEp7Sr5yThxca2bXy7CLBZCgi1qaKArPI49a43K6z1OHUnCS821m3T5z3h/CqA1VFAqE0VBcSVelbDScKLjX3buMIcwOYoINSmigIy1L0KxnZDNYefLWbbNLq+J5wDArA6Cgi1qaKATCarOU59ni6ZbTPvHIZlJ4g3/dsHqjwPYhXbZuc5NW6fyaT/+0aJA+hHAaE21RSQVRynPk9NBWS3w2jOeetHlppkL5o4bushSKvcNue+bX5+m7fPtL7vm7GfXwXQlwJCbaopIJPJ3o9Tn6eWAtLlHIZVLNs4gVzXttnW7TNrL++bsZ5fBbAXCgi1qaqA7Oh7nHrfzDboevjLKpZtO4RmndtmG7fPrFW8b8Z2fhXAXigg1KbKAjKZrK5c1FBAlj2HYRXLtkwoN7Fttmn7zOO9BbBeCgi1UUBGMEla9hKoq1i25TKqm9g227R95vHeAlgvBYTaKCAjmCQtexO4VSzbciO5TWybbdo+83hvAayXAkJtFJA9ZrbhcqvLfsq/6CpONX7Cv6pts+yVskrcPl33ZQUEYL0UEGqjgPTMbNPVfPrcx2L2JOGb/u0DS7/GNljFtul7r5BSbOrqcQoIQDcKCLVRQHpkuly2tbTLra7iRo213khuVTex3Mbt02dfVkAA1ksBoTYKSI/MNk40V3GjxmVfo7TD0xaNZ1U3sdzGG+2tqnzNUkAAVkcBoTYKyJKZsR9q0+U1Sjs8bYgxj/XQvMlEAQFYNwWE2iggS2aWPWl53Scbd/nWYRU3alz0GqUdnrbseFZ1E8sur7Ppb4j67ssKCMB6KSDURgFZMrPsZVvXdbnVTZxIvA2Hp5V4iFEp35L03ZcVEID1UkCojQJSwTcgmzqRuPTD00o8xKikb4h8AwKwHRQQaqOALJkpbZI9mWzuU/7Sy1mJE+ySviEqsaABcCwFhNooICaRS/3bd3u8tMPTSjvESHldPgOAAkJ9FJAemZIut7rJT/l9A7L9h++t6hLEx/u3980AoIBQHwWkZ8aJxOV/wl/aIUalfUO0w53QAcqmgDDrhCSXJflckkfbPy9r1+9mX5odadHy9Zn8gV2yz9jD+BWQPWa6XrZ1KCV9AzKZlHV4Wt/xjOkbkGmrugTxKjMAKCAc611pdoh3J3ll++ckyTuP87wnJ7lozvLL7fNvnMkfaNf/5JznPGkP41dAtnwiVdI5IJPJ5u6WPvRdzpfdNvMypX1D1GXMm84AoIBwtDOTfCPJNTPrr2nXn9njNd+UZgd72cz6A+36E3u85m4UkAomUqVcBWvHOu+Wvs67nK8iU9o3RF3GvMkMAAoIR7syzc5w+sz609v1V/Z4zc8k+dMk3zSz/kD7micl+eYkT+jx2vMoIBVMpDb1KX/fw9NWdS+MTd3lfC+Zki5g0HXMm8wAoIBwtFuTfHHBY19KcsuSr3dump3rf5nz2IH2sa+2f/5lkhuSPGvJnzFLAalkIlXKndC7ZFb1LUBJ53cskynlAgbzlLB9ADiaAsK0u5PcueCxu9rHl/GbaXau75nz2K8meUeSC5NckOTtSb6W5KEkp3V8/VPS7LjTy4VRQAbPDH2ew7R1fsq/yfMgSrvCVZ9Ml/9W69x3VvXvWmUGAAWEo30+yR0LHrsjzRWxunpikoeT3LPEc344zc54oGN+fxZcSUsBKfuT7pLOYdjra/S5EtS8yXNp9/hYdab2fadrBgAFhKOt8huQf5hmx/qlJcdwZ5IHO2Z9A7LGzKbOc1jnv73Payx7L4wXXHn7/PVXzF+/aBn6LuerzIxh3+maAUAB4WirPAfkXyb56zQlYRkfSPLYks+Z5hyQgTJjOM9hHd+ArGrZpm9AxrDvdM0AoIBwtLdmNVfB+s405ePmHmP4VJI/6fG8HQrIAJmxnOewjm2zqqXEc0BWsX22dd/pmgFAAeFoz8vu9wE5q/37SUmek8XfbvxCmp3qHy54/MlpzhGZ9ePt867rPuRjKCADZFZ1x+vSz3MY+hP+VS0lf8I/1n2nawYABYRjXZcjd0K/JEfuhD5dCvZl95PFP5XmBPR5JSNJzk7yQJq7q78hyeuS/E6Sx5N8IcnT+g9fARkis+x5DvtvvGfuidZ9Xmed//a+r9HlXhirWko/x2Gs+07XDAAKCMc6McnlSQ6mORfjYPv36TuW78viAvJftI9du8vP+M4k701zk8LDSR5Nc4Wtq5M8dS+DjwJSxKfY577tI3PXn/PW+etL+RR7L6+x29WZlj3BfNGJ6ttwlaex7jtdMwAoINRHARkgM5bzHFbxGvPuhdHn8KES7nK+isxY9p2uGQAUEOqzdQVkW26oNobzHEqZhJc6ebbv7C2zqpswAmw7BYTabE0B2bYbqo3hPIcSJuGlTp73khnDvrNbZlU3YQSohQJCbbaigGzrDdV2m0gte4z+omP9t+1O6F0zXSbhJf43X1Wm9n1nUWZVN2EEqIkCQm22ooBs+6fh885P6HOVotLOcyjtk/CS/puvKlPrvrMos6qbMALURAGhNsUXkFrPB6jhPg2bnIRvcjybztSw78zLrOomjAC1UUCoTfEFxGSr3GJV0ljGlKlh35mXWdVNGAFqo4BQm+ILSC03VFvV4SabHvMmfo5MffvOvEzf9zpA7RQQalN8Aan1G5DJZPtPtC5pLGPLbPu+My/jGxCA+RQQalN8Aan1cJMd23yidUljGWNmm/edeRnngADMp4BQm+ILyGRS5+Ems7bxROuSxjLmzDbuO4syroIFcCwFhNpsRQGp8XCTGjIljUWmjkyf9zpA7RQQarMVBWQyqe9wkxoyJY1Fpp6MO6EDHE0BoTZbU0B21HS4ybZnShqLTH2Zru91gNopINRm6wrIZFLWJGnMmZLGIjPeDEDtFBBqo4DI9M6UNBaZ8WYAaqeAUBsFRKZ3pqSxyIw3A1A7BYTaKCAyvTMljUVmvBmA2ikg1EYBkemdKWksMuPNANROAaE2CohM70xJY5EZbwagdgoItVFAZHpnShqLzHgzALVTQKiNAiLTO1PSWGTGmwGonQJCbRQQmd6ZksYiM94MQO0UEGqjgMj0zpQ0FpnxZgBqp4BQGwVEpnempLHIjDcDUDsFhNooIDK9MyWNRWa8GYDaKSDURgGR6Z0paSwy480A1E4BoTYKiEzvTEljkRlvBqB2Cgi1UUBkemdKGovMeDMAtVNAqI0CItM7U9JYZMabAaidAkJtFBCZ3pmSxiIz3gxA7RQQaqOAyPTOlDQWmfFmAGqngFAbBUSmd6aksciMNwNQOwWE2iggMr0zJY1FZrwZgNopINRGAZHpnSlpLDLjzQDUTgGhNgqITO9MSWORGW8GoHYKCLVRQGR6Z0oai8x4MwC1U0CojQIi0ztT0lhkxpsBqJ0CQm0UEJnemZLGIjPeDEDtFBBqo4DI9M6UNBaZ8WYAaqeAUBsFRKZ3pqSxyIw3A1A7BYTaKCAyvTMljUVmvBmA2ikg1EYBkemdKWksMuPNANROAWHWCUkuS/K5JI+2f17Wrj+eA2l2pnnLM+bkn5vk5iRfbZeb23V7oYDI9M6UNBaZ8WYAaqeAMOtdaXaIdyd5ZfvnJMk7Ozz3QJv9ySQXzSxPmsk+O8lXkhxM8sZ2uT/Jw+1jfSkgMr0zJY1FZrwZgNopIEw7M8k3klwzs/6adv2Zx3n+gTQ704kdftb7kxxOcurUulPbde/r8PxFFBCZ3pmSxiIz3gxA7RQQpl2ZZmc4fWb96e36K4/z/ANt7qQk35zkCQtyT0lzeNd75jz2niRfazN9KCAyvTMljUVmvBmA2ikgTLs1yRcXPPalJLcc5/kH0uxMX23//MskNyR51kzunPbx18x5jde2j72w04iPpYDI9M6UNBaZ8WYAaqeAMO3uJHcueOyu9vHd/GqSdyS5MMkFSd6e5tuMh5KcNpW7IM1Od/6c13hp+9gFHcZ7Spodd3q5MAqITM9MSWORGW8GoHYKCNM+n+SOBY/dkeaKWMv64TQ72IGpda9o1503J39e+9hFHV57fxZcdUsBkemTKWksMuPNANROAWHaXr8BWeTOJA9O/d03IDNKmwCNNVPSWGTGmwGonQLCtL2eA7LIB5I8NvV354DMKG0CNNZMSWORGW8GoHYKCNPemr1dBWuRTyX5k6m/H+8qWI/GVbBkNpApaSwy480A1E4BYdrzsvt9QM5q/35SkuekOQRqx5OTPHHOa/54mh3supn1N6S558f0HdJ37gNyQ4+x71BAZHpnShqLzHgzALVTQJh1XY7cCf2SHLkT+nSB2JdjTyw/O8kDae6Y/oYkr0vyO0keT/KFJE+b+TlnJDmU5k7ob2iXg+26M/YwfgVEpnempLHIjDcDUDsFhFknJrk8TRl4rP3z8hx9d/N9ObaAfGeS9yb5TJpvMR5Nc9Wsq5M8dcHPOivNeSWPtMuHc+Rblr4UEJnemZLGIjPeDEDtFBBqo4DI9M6UNBaZ8WYAaqeAUBsFRKZ3pqSxyIw3A1A7BYTaKCAyvTMljUVmvBmA2ikg1EYBkemdKWksMuPNANROAaE2CohM70xJY5EZbwagdgoItVFAZHpnShqLzHgzALVTQKiNAiLTO1PSWGTGmwGonQJCbRQQmd6ZksYiM94MQO0UEGqjgMj0zpQ0FpnxZgBqp4BQGwVEpnempLHIjDcDUDsFhNooIDK9MyWNRWa8GYDaKSDURgGR6Z0paSwy480A1E4BoTYKiEzvTEljkRlvBqB2Cgi1UUBkemdKGovMeDMAtVNAqE1RBeTTDxyaXHXbfZP9N94zueq2+yb3Pnhobq60CdBYMyWNRWacma6/MwC2mQJCbYooIPc/dHjy8mvvmJx26U3HLC+/9o7J/Q8dPipf0gRozJmSxiIzrsyyvzMAtpkCQm02XkAOPnR4cvZbbp07kdhZzn7LrUdNKEqYAMmUNRaZ8WT6/M4A2GYKCLXZeAFZ9CnmvE81d2x6AiRT3lhkxpPp8zsDYJspINRmowXk0w8c6jSR2Fl2ju/e9ARIpryxyIwj0/d3BsA2U0CozUYLyFW33bfUZOLq2++bTCblT5LGkilpLDLjyPT9nQGwzRQQarPRArL/xnuWmkzsv/GeyWRS/iRpLJmSxiIzjkzf3xkA20wBoTa+AZHpnSlpLDLjyPgGBBgjBYTaOAdEpnempLHIjCPjHBBgjBQQauMqWDK9MyWNRWY8GVfBAsZGAaE2Gy8g97sPyNZmShqLzHgyfX5nAGwzBYTabLyATCbuhL6tmZLGIjOujDuhA2OigFCbIgrIjnsfPDR5/hW3T0679KbJ86+4feHx2yVMgGTKGovMODNdf2cAbDMFhNoUVUAmk7ImNzIKiEwdGYBtpoBQGwVEpnempLHIyCggQK0UEGqjgMj0zpQ0FhkZBQSolQJCbRQQmd6ZksYiI6OAALVSQKiNAiLTO1PSWGRkFBCgVgoItVFAZHpnShqLjIwCAtRKAaE2CohM70xJY5GRUUCAWikg1EYBkemdKWksMjIKCFArBYTaKCAyvTMljUVGRgEBaqWAUBsFRKZ3pqSxyMgoIECtFBBqo4DI9M6UNBYZGQUEqJUCQm0UEJnemZLGIiOjgAC1UkCojQIi0ztT0lhkZBQQoFYKCLVRQGR6Z0oai4yMAgLUSgGhNgqITO9MSWORkVFAgFopINRGAZHpnSlpLDIyCghQKwWE2iggMr0zJY1FRkYBAWqlgDDrhCSXJflckkfbPy9r1+/mSUlek+TmJH+S5C+T3Jvk15P8zTn5A2l2vHnLM/YwfgVEpnempLHIyCggQK0UEGa9K80O8e4kr2z/nCR553Ge99wk30jy8SSXt8+9NsnXk3w2yTfP5A+0r/uTSS6aWZ60h/ErIDK9MyWNRUZGAQFqpYAw7cw0JeKamfXXtOvP3OW5T01y1pz1P51mB/v5mfUH2vUn9hnoLhQQmd6ZksYiI6OAALVSQJh2ZZqd4fSZ9ae366/s8Zrf3D73t2fWH2jXn9RmntDjtedRQGR6Z0oai4yMAgLUSgFh2q1JvrjgsS8luaXHa56RZgd728z6A+36r7Z//mWSG5I8q8fPmKaAyPTOlDQWGRkFBKiVAsK0u5PcueCxu9rHl/W/pTl863kz6381yTuSXJjkgiRvT/K1JA8lOa3ja5+SZsedXi6MAiLTM1PSWGRkFBCgVgoI0z6f5I4Fj92R5opYy3hVmp3rHR3zP9zmD3TM78+CK2kpIDJ9MiWNRUZGAQFqpYAwbZXfgPxokr9OcmOWO9H8ziQPdsz6BkRmpZmSxiIjo4AAtVJAmLaqc0DOS3M41f+Z5G8sOYYPJHlsyedMcw6ITO9MSWORkVFAgFopIEx7a/Z+FazvT3NC+R8leUqPMXwqzY0M+1JAZHpnShqLjIwCAtRKAWHa87L7fUB27vNxUpLnpDkEatr3JnkkTYn41l1+zpOTPHHO+h9PszNet9Soj6aAyPTOlDQWGRkFBKiVAsKs63LkTuiX5Mid0KdLwb4ce7L4aUn+PM2dz38hx97d/IensmcneSDN3dXfkOR1SX4nyeNJvpDkaXsYvwIi0ztT0lhkZBQQoFYKCLNOTHJ5koNpzsU42P59+kTyfTm2gLwoC65I1S4fm8p+Z5L3JvlMksNJHk1zha2r09xRfS8UEJnemZLGIiOjgAC1UkCojQIi0ztT0lhkZBQQoFYKCLVRQGR6Z0oai4yMAgLUSgGhNgqITO9MSWORkVFAgFopINRGAZHpnSlpLDIyCghQKwWE2iggMr0zJY1FRkYBAWqlgFAbBUSmd6akscjIKCBArRQQaqOAyPTOlDQWGRkFBKiVAkJtFBCZ3pmSxiIjo4AAtVJAqI0CItM7U9JYZGQUEKBWCgi1UUBkemdKGouMjAIC1EoBoTYKiEzvTEljkZFRQIBaKSDURgGR6Z0paSwyMgoIUCsFhNooIDK9MyWNRUZGAQFqpYBQGwVEpnempLHIyCggQK0UEGqjgMj0zpQ0FhkZBQSolQJCbRQQmd6ZksYiI6OAALVSQKiNAiLTO1PSWGRkFBCgVgoItVFAZHpnShqLjIwCAtRKAaE2CohM70xJY5GRUUCAWikg1EYBkemdKWksMjIKCFArBYTaKCAyvTMljUVGRgEBaqWAUBsFRKZ3pqSxyMgoIECtFBBqo4DI9M6UNBYZGQUEqJUCQm0UEJnemZLGIiOjgAC1UkCojQIi0ztT0lhkZBQQoFYKCLVRQGR6Z0oai4yMAgLUSgGhNgqITO9MSWORkVFAgFopINRGAZHpnSlpLDIyCghQKwWE2iggMr0zJY1FRkYBAWqlgFAbBUSmd6akscjIKCBArRQQaqOAyPTOlDQWGRkFBKiVAkJtFBCZ3pmSxiIjo4AAtVJAqI0CItM7U9JYZGQUEKBWCgi1UUBkemdKGouMjAIC1EoB2awTkjwryfe2y7PadfSngMj0zpQ0FhkZBQSolQKyft+c5HVJPprkcJLHZ5bD7WOvTfItGxrjNlNAZHpnShqLjIwCAtRKAVmfb03yjiR/kaZo3J3kt5JcnuTVaQrH5Ul+O8k9beYvkvyT9rl0o4DI9M6UNBYZGQUEqJUCsj4PJ/lskl9I8owO+Wck+aX2OQ8POK7aKCAyvTMljUVGRgEBaqWArM8rkjyhx/OekOSiFY+lZgqITO9MSWORkVFAgFopINRGAZHpnSlpLDIyCghQKwWE2iggMr0zJY1FRkYBAWqlgGzWE5K8OMnPJvnHSd48s/zjzQ1taykgMr0zJY1FRkYBAWqlgGzOWUk+n+ZqV99YsDy+sdFtLwVEpnempLHIyCggQK0UkM35v5L8aZL/Jsnf3PBYpp2Q5LIkn0vyaPvnZel+g8TnJrk5yVfb5eZ23V6zXa2tgHz6gUOTq267b7L/xnsmV9123+TeBw/NzZU2cZFRQGS2N9P19w5AyRSQzfmrJJduehBzvCvNDvHuJK9s/5wkeWeH5z47yVeSHEzyxna5P81lhJ+9h+wyBi8g9z90ePLya++YnHbpTccsL7/2jsn9Dx0+Kl/KxEVGAZHZ3syyv3cASqaAbM6/T/Jzmx7EjDPTHPp1zcz6a9r1Zx7n+e9Pcyf3U6fWndque98esssYtIAcfOjw5Oy33Dp3ErCznP2WW4+aDGx64iLTPVPSWGRkdjJ9fu8AlEwB2Zw3JfmjdD+0aR2uTLMznD6z/vR2/ZW7PPcpaQ7Zes+cx96T5GttZtnssgYtIIs+gZz3ieSO0ic3MmWORUZmJ9Pn9w5AyRSQzTkvyb9O8m/SHOr0w0l+YM6yTrcm+eKCx76U5JZdnntOmh3pNXMee2372At7ZJc1WAH59AOHOk0CdpadY7NLn9zIlDkWGZmX/MYnev/eASiZArI58654Nb1s4ipYdye5c8Fjd7WPL3JBmh3p/DmPvbR97IIe2d2ckmbHnV4uzEAF5Krb7jvm/+zf+YP/7eRDL/jBucsf/NBLJgf/wcsmt33feZMPveAHJ7d933mTg//gZXMXmTIyJY1FRua27ztv8gc/9JKFv2N2+73z8Re9eHLb9503+fiLXrzw58jIDJkpaSwy3TMf+rGfXvn8aR4FZHN+quOyTp9PcseCx+5Ic0WsRV6RZkc6b85j57WPXdQju5v9bfaYZYgCsv/Ge44pIB96wQ9O7j3jORaLxWKxWCxbv9z2feetfP40jwLCNN+A7KLvNyClfsohU/ZYZGQ+/qIX9/4GpORvdWTGkSlpLDLdM9f/yE+ufP40jwLCNOeA7MKx2MC6OfdMZlszJY1FZjWZVVJA1ucfpd8Vr05I8tMrHssib81wV8F6NN2vgjWdXVZxV8EC2AtX35PZxkxJY5FRQMbsT9OcY/GLSZ7eIf+MJL/UPuehAcc17XnZ/T4gZ7V/PynJc9IcAjXthjT38XjG1Lqde3vcsIfsMgYtIPe7Hj+wZn1+75Q2cZEZX6akscgoIGP2LUl+Pc0E+6+TfCrJbyf5lSSva5dfSXPn8bvbzCNJfrV97rpclyN3Qr8kR+6Eft1UZl+77sDMc89IcijN3c3f0C4H23Vn7CG7jOLuhA6wV8v+3ilt4iIzvkxJY5FRQEj+0zTnOdyWpmDMXo73q+1jr2mz63ZiksvTlIHH2j8vb9fv2Jf5BSRpviW5Jc2/7ZEkHxDXxHcAACAASURBVM6Rb072ku1q8AKy494HD02uvv2+yf4b75lcfft9zvkABtf1905pExeZ8WVKGouMAsLRTkgzmX9+u+xL8oQNjqcGaysgAKUqbeIiM75MSWORUUA44uwOmZ8YfBT1UUCA0Stt4iIzvkxJY5FRQDji0TQnpM/zzUnem/XfCb0GCggweqVNXGTGlylpLDIKCEf8XppzPj6ao6+K9V8luT9NQfnlDYxr2ykgwOiVNnGRGV+mpLHIKCAc7afSXPXpz9LcwfvX0lz96t8l+TsbHNc2U0CA0Stt4iIzvkxJY5FRQDjWviT/d5rDrR5P8ptJTt7kgLacAgKMXmkTF5nxZUoai4wCwrFem+Qv0lx+9xtJ/jDJszY6ou2mgACjV9rERWZ8mZLGIqOAcMR3JPn9NKXjxvbvP5bmjumPpLkJIMtTQIDRK23iIjO+TEljkVFAOOJLab75eN3M+qcnuT3N4VgfWPegKqCAAKNX2sRFZnyZksYio4BwxF1JnrPL4z+f5K/WNJaaKCDA6JU2cZEZX6akscgoIBxxUofMcwcfRX0UEGD0Spu4yIwvU9JYZBQQGJoCAoxeaRMXmfFlShqLjAICQ1NAgNErbeIiM75MSWORUUBgaAoIMHqlTVxkxpcpaSwyCggMTQEBRq+0iYvM+DIljUVGAYGhKSDA6JU2cZEZX6akscgoIDA0BQQYvdImLjLjy5Q0FhkFBIamgACjV9rERWZ8mZLGIqOAwNAUEGD0Spu4yIwvU9JYZBQQGJoCAoxeaRMXmfFlShqLjAICQ1NAgNErbeIiM75MSWORUUBgaAoIMHqlTVxkxpcpaSwyCggMTQEBRq+0iYvM+DIljUVGAYGhKSDA6JU2cZEZX6akscgoIDA0BQQYvdImLjLjy5Q0FhkFBIamgACjV9rERWZ8mZLGIqOAwNAUEGD0Spu4yIwvU9JYZBQQGJoCAoxeaRMXmfFlShqLjAICQ1NAgNErbeIiM75MSWORUUBgaAoIMHqlTVxkxpcpaSwyCggMTQEBRq+0iYvM+DIljUVGAYGhKSDA6JU2cZEZX6akscgoIDA0BQQYvdImLjLjy5Q0FhkFBIamgACjV9rERWZ8mZLGIqOAwNAUEGD0Spu4yIwvU9JYZBQQGJoCAoxeaRMXmfFlShqLjAICQ1NAgNErbeIiM75MSWORUUBgaAoIMHqlTVxkxpcpaSwyCggMTQEBRq+0iYvM+DIljUVGAYGhKSDA6JU2cZEZX6akscgoIDA0BQQYvdImLjLjy5Q0FhkFBIamgACjV9rERWZ8mZLGIqOAwNAUEGD0Spu4yIwvU9JYZBQQtsNzk9yc5KvtcnO7rovvSfJPk3wqySNJvpjko0l+aE52X5qdb97yWz3HroAAo1faxEVmfJmSxiKjgFC+Zyf5SpKDSd7YLvcnebh97Hjen+ShJNcleVWSn09yd5qd7LUz2X3t+g8muWhmeWHP8SsgwOiVNnGRGV+mpLHIKCCU7/1JDic5dWrdqe2693V4/t9N8sSZdScnuS/Jnyc5cWr9vjQ73xU9xzqPAgKMXmkTF5nxZUoai4wCQtmekuTRJO+Z89h7knytzfTxjjQ72nSx2ZcjBeTkdtkrBQQYvdImLjLjy5Q0FhkFhLKdk2ZneM2cx17bPtb30Kh/nuTrSZ48tW5f+5qP5Mi5H59N8rM9f0aigAAUN3GRGV+mpLHIKCCU7YI0O8P5cx57afvYBT1e97vSfLNyw8z6Zyb5SJLXtz/z1Un+uP05b+/wuqek2XGnlwujgAAjV9rERWZ8mZLGIqOAsD4nJXlOx+WZ7XNekWZnOG/O653XPnbRkuP4liT3pjmJ/ZnHySbJCUk+nuTxJM86TnZ/FlxFSwEBxqy0iYvM+DIljUVGAWF99mXxJW5nl4+1z1n1NyAnJ/lEkr9K8v1LPO/89me96jg534AAzFHaxEVmfJmSxiKjgLA+Jyd5ccflBe1zVnkOyDcluSXJY0l+ZMmxn9X+rMuXfF7iHBCA4iYuMuPLlDQWGQWEsh3vKliPpttVsE5Mc2+Px5P8RI9xvCzNTnlJj+cqIMDolTZxkRlfpqSxyCgglO+GNPf8eMbUup37gMyeRH5KmnNITppa94Qk/yLNTvUzx/lZ3zZn3clJ7krzzcmpcx4/HgUEGL3SJi4y48uUNBYZBYTynZHkUJo7ob+hXQ62686YyR5Is/Psm1p3VY6cVzJ7d/OLkjxtKvuBNIdp7U/yyiRvTvL59vmX9Ry/AgKMXmkTF5nxZUoai4wCwnY4K00xeKRdPtyum3UgxxaQj2X3E95fNJW9JM1J6l9Kc4+Qh5N8NM0J730pIMDolTZxkRlfpqSxyCggMDQFBBi90iYuMuPLlDQWGQUEhqaAAKNX2sRFZnyZksYio4DA0BQQYPRKm7jIjC9T0lhkFBAYmgICjF5pExeZ8WVKGouMAgJDU0CA0Stt4iIzvkxJY5FRQGBoCggweqVNXGTGlylpLDIKCAxNAQFGr7SJi8z4MiWNRUYBgaEpIMDolTZxkRlfpqSxyCggMDQFBBi90iYuMuPLlDQWGQUEhqaAAKNX2sRFZnyZksYio4DA0BQQYPRKm7jIjC9T0lhkFBAYmgICjF5pExeZ8WVKGouMAgJDU0CA0Stt4iIzvkxJY5FRQGBoCggweqVNXGTGlylpLDIKCAxNAQFGr7SJi8z4MiWNRUYBgaEpIMDolTZxkRlfpqSxyCggMDQFBBi90iYuMuPLlDQWGQUEhqaAAKP26QcOTZ5/xe2T0y69afL8K26f3Pvgobm50iY3MnVlShqLzO6Zrr8zVkkBoTYKCDBK9z90ePLya++YnHbpTccsL7/2jsn9Dx0+Kl/SBEimvkxJY5GZn1n2d8YqKSDURgEBRufgQ4cnZ7/l1rkTiZ3l7LfcetSEooQJkEy9mZLGInNsps/vjFVSQKiNAgKMzqJPMed9qrlj0xMgmbozJY1F5thMn98Zq6SAUBsFBBiVTz9wqNNEYmfZOb570xMgmbozJY1F5uhM398Zq6SAUBsFBBiVq267b6nJxNW33zeZTMqfJMlsd6akscgcnen7O2OVFBBqo4AAo7L/xnuWmkzsv/GeyWRS/iRJZrszJY1F5uhM398Zq6SAUBsFBBgV34DIlJgpaSwyvgGBoSkgwKg4B0SmxExJY5FxDggMTQEBRsdVsGRKy5Q0FpljM66CBaulgACjc7/7gMgUlilpLDLHZvr8zlglBYTaKCDAKLkTukxJmZLGIjM/407osDoKCDBq9z54aHL17fdN9t94z+Tq2+9bePx2CRMgmXozJY1FZvdM198Zq6SAUBsFBKCDkiZAMvVlShqLzPEz66aAUBsFBKCD0iZAMnVlShqLjAICQ1NAADoobQIkU1empLHIKCAwNAUEoIPSJkAydWVKGouMAgJDU0AAOihtAiRTV6akscgoIDA0BQSgg9ImQDJ1ZUoai4wCAkNTQAA6KG0CJFNXpqSxyCggMDQFBKCD0iZAMnVlShqLjAICQ1NAADoobQIkU1empLHIKCAwNAUEoIPSJkAydWVKGouMAgJDU0AAOihtAiRTV6akscgoIDA0BQSgg9ImQDJ1ZUoai4wCAkNTQAA6KG0CJFNXpqSxyCggMDQFBKCD0iZAMnVlShqLjALCdnhukpuTfLVdbm7XdbEvzQ41b/mtOfkTklyW5HNJHm3/vKxd34cCAtBBaRMgmboyJY1FRgGhfM9O8pUkB5O8sV3uT/Jw+9jx7EuzQ30wyUUzywvn5N/V5t+d5JXtn5Mk7+w5fgUEoIPSJkAydWVKGouMAkL53p/kcJJTp9ad2q57X4fn70uzQ13RIXtmkm8kuWZm/TXt+jM7vMYsBQSgg9ImQDJ1ZUoai4wCQtmekuYwqPfMeew9Sb7WZnazL0cKyMntssiVbfb0mfWnt+uvPO6Ij6WAAHRQ2gRIpq5MSWORUUAo2zlpdobXzHnste1j8w6jmravzT2SI+d+fDbJz87J3prkiwte50tJbjnuiI+lgAB0UNoESKauTEljkVFAKNsFaXaG8+c89tL2sQuO8xrPTPKRJK9vX+fVSf64fe7bZ7J3J7lzwevc1T6+m1PS7LjTy4VRQACOq7QJkExdmZLGIqOAsD4nJXlOx+WZ7XNekWZnOG/O653XPnZRj7GckOTjSR5P8qyp9Z9PcseC59yR5opYu9mfBVfcUkAAdlfaBEimrkxJY5FRQFiffVl8OdzZ5WPtc1bxDcgi57fPf9XUOt+AAGxIaRMgmboyJY1FRgFhfU5O8uKOywva56ziHJBFzmqff/nUOueAAGxIaRMgmboyJY1FRgGhbMe7CtajOf5VsBZ5WZod7ZKpdW+Nq2ABbERpEyCZujIljUVGAaF8N6S558czptbt3AfkhpnsKWnOITlpat23zXnNk9McUvVYjr6/yPOy+31Azlpy7IkCAtBJaRMgmboyJY1FRgGhfGckOZTmTuhvaJeD7bozZrIH0uw8+6bWfSDNoVP709zZ/M1pTjafJLlszs+7LkfuhH5JjtwJ/bqe41dAADoobQIkU1empLHIKCBsh7PSlIhH2uXDmf9txIEcW0AuSfKJNOdwfD3Jw0k+muYk9nlOTHNeyME035AcbP9+Ys+xKyAAHZQ2AZKpK1PSWGQUEBiaAgLQQWkTIJm6MiWNRUYBgaEpIAAdlDYBkqkrU9JYZBQQGJoCAtBBaRMgmboyJY1FRgGBoSkgAB2UNgGSqStT0lhkFBAYmgIC0EFpEyCZujIljUVGAYGhKSAAHZQ2AZKpK1PSWGQUEBiaAgLQQWkTIJm6MiWNRUYBgaEpIAAdlDYBkqkrU9JYZBQQGJoCAtBBaRMgmboyJY1FRgGBoSkgAB2UNgGSqStT0lhkFBAYmgIC0EFpEyCZujIljUVGAYGhKSAAHZQ2AZKpK1PSWGQUEBiaAgLQQWkTIJm6MiWNRUYBgaEpIAAdlDYBkqkrU9JYZBQQGJoCAtBBaRMgmboyJY1FRgGBoSkgAB2UNgGSqStT0lhkFBAYmgIC0EFpEyCZujIljUVGAYGhKSAAHZQ2AZKpK1PSWGQUEBiaAgLQQWkTIJm6MiWNRUYBgaEpIAAdlDYBkqkrU9JYZBQQGJoCAtBBaRMgmboyJY1FRgGBoSkgAB2UNgGSqStT0lhkFBAYmgIC0EFpEyCZujIljUVGAYGhKSAAHZQ2AZKpK1PSWGQUEBiaAgLQQWkTIJm6MiWNRUYBgaEpIAAdlDYBkqkrU9JYZBQQGJoCAtBBaRMgmboyJY1FRgGBoSkgAB2UNgGSqStT0lhkFBAYmgIC0EFpEyCZujIljUVGAYGhKSAAx/HpBw5Nnn/F7ZPTLr1p8vwrbp/c++ChubnSJkky25MpaSxjznR9r6+bAkJtFBCABe5/6PDk5dfeMTnt0puOWV5+7R2T+x86fFS+pImUzHZlShrLGDPLvtfXTQGhNgoIwBwHHzo8Ofstt86dkOwsZ7/l1qMmJiVMpGS2M1PSWMaW6fNeXzcFhNooIABzLPo0dN6nozs2PZFaNuPQsnIyJY2lS6amfafPe33dFBBqo4AAzPj0A4c6TUh2lp3J16YnUl0zDi0rL1PSWHbL1Lbv9H2vr5sCQm0UEIAZV91231KTkqtvv28ymZQ/2ZpMHFpWaqaksSzK1Ljv9H2vr5sCQm0UEIAZ+2+8Z6lJyf4b75lMJuVPtiaTcRxato2ZksayKFPjvtP3vb5uCgi1UUAAZtT6DUjth5Ztc6aksczL1Lrv+AYENkMBAZhhsrVdxWqozDpPtC7p3z0vU+u+4xwQ2AwFBGCObT/cZN7kuc/hJjVd7ahrZhMnWpfw794xtn3HVbBg/RQQgDnu39ITbnebPJ/z1o8sNYk8923z89t4taOumU2daL3pf/dkMt59p897fd0UEGqjgAAssG2XHO0yeV7VUmL5WkVmU998bfrfPfZ9x53QYb0UEIDjuPfBQ5Orb79vsv/GeyZX337f1h9Ksqql1MPP+mY2ee7PpreNfafR9b2+bgoItVFAAFZkmybPq1q25QT8WfPOT+hzovWqznOw75SxX5RKAWGe5ya5OclX2+Xmdl0XB9LsUIuWX5nK7tsl91s9x66AAKzIJidSy06eFx2jv+yx/ttytaMdux1q84K2SHRdXnDl/Hyfw/PsO+XsOyVSQJj17CRfSXIwyRvb5f4kD7ePHc85SS6as/xxjt3R9rXrPjgn/8Ke41dAAFZkkxOpPlcpmne4Sc03YSz5PAf7Tjn7TokUEGa9P8nhJKdOrTu1Xfe+nq/5pDTfpHxqZv2+NDvfFT1fdx4FBGBFtulT7EU3VKv1fg+TSdnnOdh3ytl3SqSAMO0pSR5N8p45j70nydfazLL+uzQ72S/MrN+XIwXk5HbZKwUEYEW26Tj+Recq1HoTxtLPc9imbVP7vlMiBYRp56TZGV4z57HXto/1OTTqtiRfT/K0mfX72td8JEfO/fhskp/t8TN2KCAAK7LpidSqbqi27TdhnJdZ9tP5VS1dP+Xf9Paz75RNAWHaBWl2hvPnPPbS9rELlnzNpyd5PMlNcx57ZpKPJHl9+zNfnSPniry9w2ufkmbHnV4ujAICsBKbnkit6oZq23oTxt0yy56fsOgE82VPVO96nsOmt599p2wKSL1OSvKcjssz2+e8Is3OcN6c1zuvfeyiJcdxWfu8H+uYPyHJx9OUlmcdJ7s/C66ipYAA7F0JE6lV3VBt227CeLxMn/MT5p1oPdR5DiVsP/tOuRSQeu3L7pfDnV4+1j5niG9A7k3y50meuMRzzm9/1quOk/MNCMCASppIreqGal1eZ1X3whgyU/p5DuvaNl3+W61z31nVv2uVmRIpIPU6OcmLOy4vaJ+z6nNAnt8+551Ljv2s9nmXL/m8xDkgACsztonUtn3SXfJ5Dtvy7cYQNr1fbAMFhGnHuwrWo1nuKlj/LM3O9YLjBWe8rH3eJUs+L1FAAFZmTBOpLvfU2NSx/os+5d/UeQ5dvnUYctv0+W+1TmN63/SlgDDrhjT3/HjG1Lqd+4DcMJM9Jc05JCfNeZ1vSvKnaQ7BWuTb5qw7OcldSR7L0fci6UoBAViRMU2kSrzaUZdP+dd5nsMyP2sbvvkZypjeN30pIMw6I8mhNHdCf0O7HGzXnTGTPZBm59k353V2vsW4dJef9YEkt6Q5mfyVSd6c5PPt8y7rN3wFBGBVxjKRKvF+D8t+yj/0eQ7Ljqf0c1+GNJb3zV4oIMxzVppi8Ei7fLhdN+tAFheQD6W5ktXf3uXnXJLkE0m+lOY+IQ8n+WiaE977UkAAVmTThxitS4l3vC7tU/5lxzPUtlnVXc6HUuIFDEqkgFAbBQRgRUo4xGgdlr2nRtd7YfTNlPYpf5/xDLVt+v63GlrJFzAokQJCbRQQgBUp6RCjIZX2DUhpn/L3Gc9Yts1kUvYFDEqlgFAbBQRgRcZyiFFp54CU9il/n/GM5duhyaTMCxiUTgGhNgoIwIqYRG5mElnap/wlfQMymSivy2ZKpIBQGwUEYEXGMsGeTPrdU2MsBa2kc0Amk9Xd/2QVSjt8b1soINRGAQFYkbEcYrSjpBOJS/qUv894hp5gu4BB90yJFBBqo4AArMiYvgGZ1uWeGqu63OrQdzlflWXGM/S2mbaq+5/05RuQfhQQaqOAAKzIWA4xWsaqviVZ513OV+V44/nkZ768tm1TCueA9KOAUBsFBGBFxnSIURerutzqpu5yvirzxrOpbVOCki5gsC0UEGqjgACsyFhOJO5qVRPNbSxfxzPmbVPSBQy2hQJCbRQQgBUZy4nEXfQ51GbeOQzbfPjZIrZNWRcw2AYKCLVRQABWZF2TpNIOMZpn2ZONz33bR+auP+et89cvWtZ9An4fts0RXfdlBUQBoS4KCMCKjH2SNG3Zy62ualnXJYj3wrZZ3tjfWwoItVFAAFZk7JOkact+yr+qpcRP+WfZNssb+3tLAaE2CgjAiox9kjRt2fMTVrWUeDjaLNtmeWN/bykg1EYBAViBVd1MrmtmG3S9QtOqlpKu9HQ8ts1yVvG+6foeLZECQm0UEIA92NTVfLZBl8utrmop7RLEx2PbLGcsV49bRAGhNgoIQE+rupncrFoKyGSy++Rv2as4Lboa1LZMImfZNt31fd9s440a51FAqI0CAtDTJu/ovG3mXW61z30stuESxMuybY6v7/tmG2/UOI8CQm0UEIAe+t4Erubj1PuoZYI4BNum0ff8qm2+UeMsBYTaKCAAPSx7KdWdS6DWfJx6H13OhdiGQ2SGMPZts9fzq/q+R0ukgFAbBQSgh2VvJrdzE7iaj1Pva6zlq4uxbptVnF/V9z1aIgWE2iggAD2s8hsQh9o0aj6HYa/Gtm1WcX6Vb0CgXAoIQA+rOgekpuPUYRW8t46lgFAbBQSgJ5/Swur5dvFYCgi1UUAAeupzknDNx6nDKqzy/KpaTuRXQKiNAgKwB67UA6vlCnPHUkCojQICsAJdTxKu+Th1WIWh7rGzzSfyKyDURgEBWJNFN1Sr5Th1WJVl3xO138BTAaE2CgjAwI53CMgnP/PlKo5Th1Xpeu7GJz/z5a0/vKoLBYTaKCAAA+p6Q7WxTKSgK8X9CAWE2iggAANa9lCSbT5OHYaw6D0xpkMXFRBqo4AADMQJ5jCMsb23FBBqo4AADMQldmEYY3tvKSDURgEBGIibDMIwxvbeUkCojQICMJCxfUoL6zK295YCQm0UEICBjO04dViXsb23FBBqo4AADGhMV+qBdRrTe0sBoTYKCMCAut5QrYZ7FcA6jem9pYBQGwUEYGDHu6FaDRMk2ISxvLcUEGqjgACsiZsMwjBqf28pINRGAQEAKJgCQm0UEACAgikg1EYBAQAomALCrPOS/K9J7kzyWJqdY9+Sr3FCksuSfC7Jo+2fl7Xr95LtQgEBACiYAsKsA2mKwL9Jcnf6FZB3tc97d5JXtn9Okrxzj9kuFBAAgIIpIMx6epIntv/7iixfQM5M8o0k18ysv6Zdf2bPbFcKCABAwRQQdtOngFzZPuf0mfWnt+uv7JntSgEBACiYAsJu+hSQW5N8ccFjX0pyS89sVwoIAEDBFBB206eA3J3mBPZ57mof75Od55Q0O+70cmEUEACAYikg9TopyXM6Ls9c8Bp9Csjnk9yx4LE70lzlqk92nv1pxnfMooAAAJRJAanXviyYnM9ZPrbgNXwDAgDASikg9To5yYs7Li9Y8BrOAQEAYKUUEHbTp4C8Nd2vbLVMtisFBACgYAoIuzleATklzTkkJ02te152v7fHWT2zXZ2bZHL99ddP7rzzTovFYrFYLBZLYcv111+/U0DO7THXo0JnJXlTu3wszc7xT9q/v34meyDzC8p1OXJ380ty5O7m1835ectku7gw3c99sVgsFovFYrFsbrkwkOTiLN5JvjCTPdCu3zez/sQklyc5mOSx9s/L2/Wzlsl28e1pduZzc+wJ6kMsO4XnwjX9vLEutrNtXMtiO9vGtSy2s228l+Xc9t/17QGW9t1pfjl896YHUjnbeXi28XrYzsOzjdfDdh6ebQzM5ZfDetjOw7ON18N2Hp5tvB628/BsY2AuvxzWw3Yenm28Hrbz8Gzj9bCdh2cbA3P55bAetvPwbOP1sJ2HZxuvh+08PNsYmOuUJPvbPxmO7Tw823g9bOfh2cbrYTsPzzYGAAAAAAAAAAAAAAAAAAAAAAAAYMNOSHJZks8lebT987J2Pct7SppLDf5+kv+Y5rrnBxZkbfvlfU+Sf5rkU0keSfLFJB9N8kNzsrZvf9+V5F8k+WySw0m+muT/SfKGJE+cydrOq/MDaX5nTJL85zOPnZzk15L8SZKvJbknyc+sdXTbaV+ObNPZ5bdmsvblvTklybVJ/n2a7fdgkg8medpM7mfS7L9fS7M//1qa/RsYkXel+UX87iSvbP+cJHnnJge1xfal2X4PpikhuxUQ235570/yUJLrkrwqyc8nuTvNdnvtTNb27e+8JLcl+Z+TvDrJ65K8N8k30uzX02zn1fimJP9vmsI3r4D8yyRfT3J1mgncB9vcL61xjNtoX5rt9MEkF80sL5zJ2pf7e3aaD90OJnlzkp9Ocmma7f7sqdz/kCP/PX4mzf7810luWudggc06M82E4pqZ9de0689c+4i23xOTPL393ydmcQGx7fv5uzn2E/iTk9yX5M/TbPPE9h3KP0uzT5/R/t12Xp1fTvKlNBOy2QLyI+26N84850NJ/jLJd6xjgFtqX5ptd8Vxcvbl/v6TJP86yZ1pjgJY5DvS7K8fmln/xjT/jV4yyOiA4lyZ5k1/+sz609v1V659RHXZrYDY9qv1jjTb7dT277bvMH4xzfb73vbvtvNqnJbkL5L8ozSHcM4WkPcm+asce5jK32uzDsVabF+OFJCTs/hQH/tyfzuHDv5I+/e/keYbvVk/0+b+3sz6J6XZv9871ACBstya5hj6eb6U5JY1jqVGuxUQ2361/nmaw1Oe3P7d9l2NJyV5apoJ8suTfDnJAzkyibOdV+NDSf4gzSfJ+3NsAbkvyR/Oed7Jbfa6gce3zfal2UaP5Mi5H59N8rMzOftyf7+eZrt+f5J/1f7vb6TZp8+Zyv1m+9i8EvhHSf7dsMMESnF3mq9M57mrfZz+disgtv3qfFeaEx5vmFpn+67G/hx90u4fJjlr6nHbee9+JMnjSf5O+/f9ObaAPJKj9+9pf55jz8vhiGcm+UiS1yc5P805TX+cZhu/fSpnX+5v53ykh5J8IM2HFT+X5M/SfLOxc/ja77fr5vlAmv0cGIHPJ7ljwWN3pLkCCP3tVkBs+9X4liT3Jnk4zURjh+27Gv9ZmiuM/XiaTy//IMmLph63nffm5DQn7f7m1Lr9ObaAPJ7kdxe8xoNpJth0d0KSj6fZrs9q19mX+/tImn12sTAsJAAAA6VJREFUdj/8L9v1v9f+/aNp9td5fjfNyejACPjEZ1i+ARnWyUk+keYTtu+fecz2HcYbkzyW5lunxHbeqyvSfIPx1Kl1++MbkHU4P812flX7d/tyfztXfPypOY99Ic0hbDs534AAjnkdmHNAhvNNabbRYzly4uM023cYT8vRJ+Tazv397TT3QfjVNGVjZ/mNNNv4B9J8A5U4B2QIZ6XZdpe3f7cv93ddmm354jmP/WGa39OJc0CA1lvjqh9D2q2A2Pb9nZjmmOPHk/zEgoztO4x9abbfu9q/2879nZ3FN8jbWQ632d+Nq2Ct2svSbLtL2r/bl/u7JM02euWcx/5DuyTNt03zroJ1clwFC0bledn9uudnHfMMlrFbAbHt+3lCmjt0H2/SZfvuzd9asH7najc7h1rYzv19S5Ifm7P8Xppt/PokP9pmdw4XWnQfkEX/vUi+bc66k9McVvVYjly6277c31PTXEb6D3P0XeN37l+zc8f5v5WmaHxw5vk79wE5f9hhAiXZ+er03Wk+xdi586uv9Pt7fZI3pbkb7CTN/9G9qV2m/0/Mtl/eVWm20cdy7F2NL0pziNAO27e//yPN+TX/U5qi90s5cqLpJ3Pkho+J7bxq+zP/Tui3pLnU9FVpPmneufLQZesc3Bb6QJpttz/NdntzmhPO5207+3J/P5dmW30izf8Hvi1NKflyjpS8pNnmkzS/Y16ZZn/+epIPr3OwwOadmOYY2INpPg062P79xN2exK6+kMWHVFw8lbPtl/ex7H7IyoumsrZvfz+eZkLwQJpt90iaS5f+Yo69E73tvFr7M7+APCnNZWP/Q5rLTt+b5pKy7O6SNJPiL6WZ6D6c5mpML52TtS/vzUVpPnD7WpqLI/zvOXIe07RXp9l/H02zP789zf4NAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAQ7s4ySTJ30/ypiR/kuRrSf5VkudtblgAAECNLk5TQP64Xd6Y5JeT/FmS+5OcuLGRAQAA1bk4TQG5K8lJU+t/tF3/X29gTAAAQKUuTlM0Xj2z/lvb9f/9ugcEAADU6+IcOQdk1iTJ/7jW0QAAAFW7OE3R+KE5j02S7F/nYAAAgLpdHAUEAABYk4ujgAAAAGtycRQQAABgTS6OAgIAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAACzv/wPuYdezfUjL2QAAAABJRU5ErkJggg==\" width=\"640\">"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/plain": [
       "Text(0,0.5,'x(n)')"
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "plt.figure()\n",
    "plt.stem(n, x)\n",
    "plt.xlabel('n')\n",
    "plt.ylabel('x(n)')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Pretpostavimo sada da su početne faze harmonika oblika\n",
    "\n",
    "$$\\phi_i = -5 \\frac{2 \\pi}{64} \\left(2i + 1 \\right), \\; i=0,1,2,3.$$\n",
    "\n",
    "Ovakve početne faze odgovaraju prolasku signala kroz svepropusni filtar sa linearnom fazom i kašnjenjem od 5 odmjeraka. Oblik signala je nepromijenjen."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [],
   "source": [
    "phi0 = -2*np.pi/64 * 5\n",
    "phi1 = -2*np.pi/64 * 3 * 5\n",
    "phi2 = -2*np.pi/64 * 5 * 5\n",
    "phi3 = -2*np.pi/64 * 7 * 5\n",
    "n = np.arange(64)\n",
    "x = np.sin(2*np.pi/64 * n + phi0) + (1/3)*np.sin(2*np.pi/64 * 3 * n + phi1) + \\\n",
    "    (1/5)*np.sin(2*np.pi/64 * 5 * n + phi2) + (1/7)*np.sin(2*np.pi/64 * 7 * n + phi3)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [
    {
     "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,iVBORw0KGgoAAAANSUhEUgAAAyAAAAJYCAYAAACadoJwAAAgAElEQVR4nO3df9AkZ2Ef+O8hyUSYsxPbCVaM0Oo4SqQKCZ0TCCi+s+PYCnU2nINMnNIJR7Ewv8wlYMcnWSHccidhJwQpygUkX9nUchFxLkYcyEQgCa4AX2QnjnRXIHQRP1ZUHCmAHMuL1jYSFnN/dL/Z1uzMuz399sw8/cznU9W12p7vzD7b6rf3+c5MdycAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAATNDTkxxO8mtJ/kOSWZIjK77G85LcluSr7XJbu+6g2T6+PcmlSS5K8t0Wi8VisVgsluKWi9LM1749kORQmtLxUJoSsmoBeU6S30tyNMmb2uWBJI+0jw3N9nVpO2aLxWKxWCwWS9nLpYEkT03yXe1/n55m5ziywvPfl+R4krM7685u1/3qAbJ9XZRkdvPNN8/uvvtui8VisVgsFkthy80337xXQC4aON+jYqsWkKcneSzJexY89p4kX2szq2ZX8d1JZnffffcMAIDy3H333XsF5LsHzPWo3KoF5MVt/rULHntd+9iLBmRXoYAAABRMAWE/qxaQS9r8Sxc89rL2sUsGZJc5Kyef1HRpFBAAgGIpIOxn1QLyyjZ/8YLHLm4fu2xAdpnDWXJSkwICAFAmBYT9+AQEAIBRKSDsxzkgAACMSgFhP2NfBeux9L8KVje7CgUEAKBgCgj72a+AnJHkuWm+BtV1S5r7eDyzs27v3h63HCDblwICAFAwBYRF3pDkzUnekmbnuKf9/ZuTXNBmDmVxOTkvybE0dzd/Y7scbdedd4BsXwoIAEDBFBAW+WKWXF0qyeVt5lCWfzpyQZKPJHm0XT6cE8XlINk+FBAAgIIpINRGAQEAKJgCQm0UEACAgikg1EYBAQAomAJCbRSQCfnMg8dm191x/+zwrffOrrvj/tl9Dx3b9pBg8vxcAaVTQKiNAjIBDzx8fPaKG++anXPlh05aXnHjXbMHHj6+7SHC5Pi5AqZCAaE2Ckjhjj58fHbhW29fOEnaWy586+0mS7ACP1fAlCgg1EYBKdyyd2gXvWML9OPnCpgSBYTaKCAF+8yDx3pNkvYW312HU/NzBUyNAkJtFJCCXXfH/StNlK6/8/5tDxmK5+cKmBoFhNooIAU7fOu9K02UDt9677aHDMXzcwVMjQJCbRSQgnmnFsbn5wqYGgWE2iggBfNddRifnytgahQQaqOAFG7Vq/W4qRqcsOznwVWwgClRQKiNAlK4B3rer+DXP/sVN1WD1qluMvjrn/2K+4AAk6GAUBsFZAJMpqC/vjcZVNqBqVBAqI0CMiH3PXRs9oJr7pydc+WHZi+45k5fJ4EFVv15WPZzBVAKBYTaKCCF6Hvuxg/9o0/OzrnyQ7Mf+kef/E/Pc0ItNIb+PMz/XC16XedWAduigFAbBWTLTvX1qvmvgcxPlFxSFE4Y+vOwrICs+vMJsA4KCLVRQLao73fVu5Oc+YmSm6rBCUN/HhYVkCE/nwDroIBQGwVki4acu+ETEFhuzE9AnFsFlEIBoTYKyJaM9V1154DACX6ugBopINRGAdkS79SOy0nC+9ul7eOTRaA2Cgi1UUC2ZMzvqve9WWGN31V3kvD+dnH7DPl5cG4VUDIFhNooIFviaj0H5yTh/e3y9nF1OaAmCgi1UUC2ZF33K9ilm6r56tn+bJ/+Pw/OAQFKpoBQGwVki8b4rvoifTJTsuj8hSETxFrPg7B99jf0Z0aBA0qhgFAbBWSLxviu+iK1FJD9vkbz4rd9dKUJ9kU/vzg/5a+n2T79DP2Z2eVzq4CyKCDURgHZsoN+V32RGgpIn/MXxlqmOIm0ffo7yM/MLp5bBZRHAaE2Ckghhn5XfWimdH2//jLWMrWv0dg+/Y3xM7NL51YB5VFAqI0CUpCxysXUC8iq5y+MtUxlUmn7rMbPFTB1Cgi1UUAKYqLUWPUSqGMtU7mUqu2zGj9XwNQpINRGASmIiVJj1ZvAjbVM5WZyts9q/FwBU6eAUBsFpCCbmiiVfrnVVd/hX3YFp1WvBDWVd/htnxP67MsKCDB1Cgi1UUAKsu6J0lSu6DPkHhaLThKu9WZyts9q+7ICAkydAkJtFJCCrHOi1OeyrSVdbnWsmzTWejO5Xd4+q+7LCggwdQoItVFACrLOidLUJppj3aRx1dcp7etpy8azre1TglX3ZQUEmDoFhNooIAVZ10Rpql+1GesmjX1ep7Svp61jzFP/at5sNmxfVkCAqVNAqI0CUpB1TZRWPWl5kycb9/nEYaybNC57ndK+nrbqeNa9fbq2/QnRkH1ZAQGmTgGhNgpIQdY1UVr1sq2buNzqWO/ej5Ep7etpY53fMa+GCxgM2ZcVEGDqFBBqo4AUZFc+ARnyicOufD1t6Hh25QIGPgEBdpECQm0UkILsyiS7pHf4SytnQ8ezK58QOQcE2EUKCLVRQAqyC5PI0t7hL+3raUPHsyvldTZzFSxg9ygg1EYBKcg6J0qlXG61tHf4fQKyf6a07TObrb4vKyDA1Ckg1EYBKcgunEjsHf79+YSoH3dCB3aJAsK805JcleTzSR5rf72qXb+fQ2l2pGXL1+fyR/bJPvMA41dACrKpiVLfy7auQ2nv8M9m5Xw97SDj2aVPQLr67MsKCDB1Cgjz3pVmh3h3kle1v86SvPMUz/vmJJctWH6uff6tc/kj7fofX/Ccpx1g/ApIQTY5UdrWZKq0d/hns+3dLX3ddzkfY/uU9glRnzGv+viYGYB1UEDoOj/JN5LcMLf+hnb9+QNe881pdrCXz60/0q4/fcBr7kcBKcguFJDZrKx3+Pds8m7pm7zL+RiZ0j4h6jPmVR4fMwOwDgoIXdem2RnOnVt/brv+2gGv+dkkv5Pkm+bWH2lf84wk35LkKQNeexEFpCC7UkBKeod/3rrvlr6tu5wfJFPKBQxWGfMqj4+ZAVgHBYSu25N8acljX07ykRVf76I0O9f/uuCxI+1jX21//YMktyR59op/xjwFpCC7UkBms7LuhN4nM9anACV++tMnU8IFDJZRQIDaKSB0fTrJ3Useu6d9fBW/mGbn+nMLHvuFJO9IcmmSS5K8PcnXkjyc5Jyer39Wmh23u1waBaQYJUw094xxnkOf1yjhHf5TZcY6D6LE819WzfT5/7Xuc2RWHbMCAkydAkLXF5LcteSxu9JcEauvpyZ5JMm9KzznB9PsjEd65g9nyZW0FJAylDDRHOOd7ql9unGqzJArQS2aPJd4BbAxM5s8R2aVMSsgwNQpIHSN+QnIX02zY/3simO4O8lDPbM+ASnctieRY5znMOQ1tv33PlVm1XthvPDaOxevv2bx+mXLuu+BMmZmW+fI9BmzAgJMnQJC15jngPyLJH+UpiSs4v1JHl/xOV3OASnItieRY5znMNVzHPbLrPrJxVjLlD4BKfkcGQUEmDoFhK63ZZyrYH1nmvJx24AxfCrJbw943h4FpCDbnESOcZ5DDec4jLFtxlpq3T6bPkdGAQGmTgGh6/nZ/z4gF7S/PyPJc7P8042fSbNT/dUlj39zmnNE5v1Y+7yb+g/5JApIQab0Lv+iO17XfI5D33fmx1pq/oRo2d3S17X/KCDA1CkgzLspJ+6EfkVO3Am9WwoOZf+TxT+V5gT0RSUjSS5M8mCau6u/Mcnrk/yTJE8k+WKSZwwfvgJSkm1OIlc9z+Fv/so9J51kveprTOkchz73whhrqf0cmcO33rvwJP117T8KCDB1CgjzTk9ydZKjac7FONr+vnvH8kNZXkD+q/axG/f5M74zyXvT3KTweJLH0lxh6/ok33GQwUcBKcqU3sVetLz4bR9dyzvYpWT2uzrTqieYLztRfReuEnbRzy/eT9a1/yggwNQpINRGASnIlL7HP8YylXMc5i26F8aQrw9N4R4ofTKlnyOjgABTp4BQGwVkA8a6odq6M5s8z2FK5zj0ydR6An7fTMnnyGxq24x1A0aAeQoItVFA1mhqN+Tb1HkOUzvHoW+mxksQ982UfI7MurfNWDdgBFhGAaE2CsiaTPWGfPtNplZZln3Pv8TiNVamzyS8xP/nY2X223dWPb9jzP1nnX/vsW7ACLAfBYTaKCBrMvV3w+fPT/ibv3LPShPIw7feW805Dqtkpvap1zoyi/6/D7nC1Vj7T2mfegGsSgGhNgrIGtR4PkDN9/hYR2YXy9d+mW3uP+v6O411A0aAU1FAqI0CsgY1TtZrLFUym8tsc/8prVQBrEoBoTYKyBrUekO+qX+tTGY3T9Jf199p6M85wKoUEGqjgKxBjZ+AzGZOspY5WGZb+49PQICpU0CojQKyBjV/XclJ1jIHyWxj/1nX38k5IMCmKCDURgFZk9q/ruQka5mDZDa5/6zz7+QqWMAmKCDURgFZk134ulJJY5GR2cZ+OuTnHGBVCgi1UUDWqPavK5U0FhmZbe2n7oQOrJsCQm0UkA2o9etKJY1FRmbb+2nfn3OAVSkg1EYB2ZCSJmRjZUoai4zMVPZTgFUpINRGAdmQ0iZBNU7sZGRKHwvAEAoItVFANqS0SVCNEzsZmdLHAjCEAkJtFJANKW0SVOPETkam9LEADKGAUBsFZENKmwTVOLGTkSl9LABDKCDURgHZkNImQTVO7GRkSh8LwBAKCLVRQDaktElQjRM7GZnSxwIwhAJCbRSQDSltElTjxE5GpvSxAAyhgFAbBWRDSpsE1Tixk5EpfSwAQygg1EYB2ZDSJkE1TuxkZEofC8AQCgi1UUA2pLRJUI0TOxmZ0scCMIQCQm0UkA0pbRJU48RORqb0sQAMoYBQGwVkQ0qbBNU4sZORKX0sAEMoINRGAdmQ0iZBNU7sZGRKHwvAEAoItVFANqS0SVCNEzsZmdLHAjCEAkJtFJANKW0SVOPETkam9LEADKGAUBsFZENKmwTVOLGTkSl9LABDKCDURgHZkNImQTVO7GRkSh8LwBAKCLVRQDaktElQjRM7GZnSxwIwhAJCbRSQDSltElTjxE5GpvSxAAyhgFAbBWRDSpsE1Tixk5EpfSwAQygg1EYB2ZDSJkE1TuxkZEofC8AQCgi1UUA2pLRJUI0TOxmZ0scCMIQCQm0UkA0pbRJU48RORqb0sQAMoYBQGwVkQ0qbBNU4sZORKX0sAEMoINRGAdmQ0iZBNU7sZGRKHwvAEAoItVFANqS0SVCNEzsZmdLHAjCEAkJtFJANKW0SVOPETkam9LEADKGAUBsFZENKmwTVOLGTkSl9LABDKCDMOy3JVUk+n+Sx9ter2vWnciTNzrRoeeaC/POS3Jbkq+1yW7vuIBSQDSltElTjxE5GpvSxAAyhgDDvXWl2iHcneVX76yzJO3s890ib/fEkl80tT5vLPifJ7yU5muRN7fJAkkfax4ZSQDaktElQjRM7GZnSxwIwhAJC1/lJvpHkhrn1N7Trzz/F84+k2ZlO7/FnvS/J8SRnd9ad3a771R7PX0YB2ZDSJkE1TuxkZEofC8AQCghd16bZGc6dW39uu/7aUzz/SJs7I8m3JHnKktzT03y96z0LHntPkq+1mSEUkA0pbRJU48RORqb0sQAMoYDQdXuSLy157MtJPnKK5x9JszN9tf31D5LckuTZc7kXt4+/dsFrvK597EW9RnwyBWRDSpsE1Tixk5EpfSwAQyggdH06yd1LHrunfXw/v5DkHUkuTXJJkren+TTj4STndHKXpNnpXrrgNV7WPnZJj/GelWbH7S6XRgHZiNImQTVO7GRkSh8LwBAKCF1fSHLXksfuSnNFrFX9YJod7Ehn3SvbdRcvyF/cPnZZj9c+nCVX3VJA1q+0SVCNEzsZmdLHAjCEAkLXQT8BWebuJA91fu8TkAqUNgmqcWInI1P6WACGUEDoOug5IMu8P8njnd87B6QCpU2CapzYyciUPhaAIRQQut6Wg10Fa5lPJfntzu9PdRWsx+IqWMUrbRJU48RORqb0sQAMoYDQ9fzsfx+QC9rfn5HkuWm+ArXnm5M8dcFr/liaHeymufW3pLnnR/cO6Xv3AbllwNj3KCAbUtokqMaJnYxM6WMBGEIBYd5NOXEn9Cty4k7o3QJxKCefWH5hkgfT3DH9jUlen+SfJHkiyReTPGPuzzkvybE0d0J/Y7scbdedd4DxKyAbUtokqMaJnYxM6WMBGEIBYd7pSa5OUwYeb3+9Ok++u/mhnFxAvjPJe5N8Ns2nGI+luWrW9Um+Y8mfdUGa80oebZcP58SnLEMpIBtS2iSoxomdjEzpYwEYQgGhNgrIhpQ2CapxYicjU/pYAIZQQKiNArIhpU2CapzYyciUPhaAIRQQaqOAbEhpk6AaJ3YyMqWPBWAIBYTaKCAbUtokqMaJnYxM6WMBGEIBoTYKyIaUNgmqcWInI1P6WACGUECojQKyIaVNgmqc2MnIlD4WgCEUEGqjgGxIaZOgGid2MjKljwVgCAWE2iggG1LaJKjGiZ2MTOljARhCAaE2CsiGlDYJqnFiJyNT+lgAhlBAqI0CsiGlTYJqnNjJyJQ+FoAhFBBqo4BsSGmToBondjIypY8FYAgFhNooIBtS2iSoxomdjEzpYwEYQgGhNgrIhpQ2CapxYicjU/pYAIZQQKiNArIhpU2CapzYyciUPhaAIRQQaqOAHNBnHjw2u+6O+2eHb713dt0d98/ue+jYwlxpk6AaJ3YyMqWPpe/xAqBLAaE2CshADzx8fPaKG++anXPlh05aXnHjXbMHHj7+pHxJk6CxMiWNRUam5P101eMFQJcCQm0UkAGOPnx8duFbb184mdhbLnzr7U+aVJQwCRo7U9JYZGRK3U+HHC8AuhQQaqOADLDsncxF72zu2fYkaB2ZksYiI1PqfjrkeAHQpYBQGwVkRZ958FivycTesvcd721PgtaRKWksMjIl7qdDjxcAXQoItVFAVnTdHfevNKG4/s77Z7NZWROysTIljUVGpsT9dOjxAqBLAaE2CsiKDt9670oTisO33jubzcqakI2VKWksMjIl7qdDjxcAXQoItVFAVuQTkDLHIiNT4n7qExBgDAoItVFAVuQckDLHIiNT4n7qHBBgDAoItVFABnAVrPLGIiNT6n7qKljAQSkg1EYBGeAB9wEpbiwyMqXup0OOFwBdCgi1UUAGcif0ssYiI1PyfupO6MBBKCDURgE5oPseOjZ7wTV3zs658kOzF1xz59LvcJcwCRo7U9JYZGSmsJ/2PV4AdCkg1EYBGUFJE5xNZkoai4xMTfspQJcCQm0UkBGUNnkxsZORKXMsfTMAXQoItVFARlDa5MXETkamzLH0zQB0KSDURgEZQWmTFxM7GZkyx9I3A9ClgFAbBWQEpU1eTOxkZMocS98MQJcCQm0UkBGUNnkxsZORKXMsfTMAXQoItVFARlDa5MXETkamzLH0zQB0KSDURgEZQWmTFxM7GZkyx9I3A9ClgFAbBWQEpU1eTOxkZMocS98MQJcCQm0UkBGUNnkxsZORKXMsfTMAXQoItVFARlDa5MXETkamzLH0zQB0KSDURgEZQWmTFxM7GZkyx9I3A9ClgFAbBWQEpU1eTOxkZMocS98MQJcCQm0UkBGUNnkxsZORKXMsfTMAXQoItVFARlDa5MXETkamzLH0zQB0KSDURgEZQWmTFxM7GZkyx9I3A9ClgFAbBWQEpU1eTOxkZMocS98MQJcCwrzTklyV5PNJHmt/vapdv5+nJXltktuS/HaSP0hyX5K/n+SPL8gfSbPjLVqeeYDxKyAjKG3yYmInI1PmWPpmALoUEOa9K80O8e4kr2p/nSV55yme97wk30jyiSRXt8+9McnXk3wuybfM5Y+0r/vjSS6bW552gPErICMobfJiYicjU+ZY+mYAuhQQus5PUyJumFt/Q7v+/H2e+x1JLliw/ifS7GA/Pbf+SLv+9CED3YcCMoLSJi8mdjIyZY6lbwagSwGh69o0O8O5c+vPbddfO+A1v6V97i/PrT/Srj+jzTxlwGsvooCMoLTJi4mdjEyZY+mbAehSQOi6PcmXljz25SQfGfCa56XZwX5+bv2Rdv1X21//IMktSZ494M/oUkBGUNrkxcRORqbMsfTNAHQpIHR9OsndSx67p318Vf97mq9vPX9u/S8keUeSS5NckuTtSb6W5OEk5/R87bPS7Ljd5dIoIAdW2uTFxE5Gpsyx9M0AdCkgdH0hyV1LHrsrzRWxVvHqNDvXO3rmf7DNH+mZP5wlV9JSQA6mtMmLiZ2MTJlj6ZsB6FJA6BrzE5AfSfJHSW7Naiea353koZ5Zn4CsSWmTFxM7GZkyx9I3A9ClgNA11jkgF6f5OtX/leSPrTiG9yd5fMXndDkHZASlTV5M7GRkyhxL3wxAlwJC19ty8KtgfW+aE8r/VZKnDxjDp9LcyHAoBWQEpU1eTOxkZMocS98MQJcCQtfzs/99QPbu83FGkuem+QpU159P8miaEvEn9vlzvjnJUxes/7E0O+NNK436yRSQEZQ2eTGxk5Epcyx9MwBdCgjzbsqJO6FfkRN3Qu+WgkM5+WTxc5L8bpo7n/9MTr67+Q92shcmeTDN3dXfmOT1Sf5JkieSfDHJMw4wfgVkBKVNXkzsZGTKHEvfDECXAsK805NcneRomnMxjra/755IfignF5Dvy5IrUrXLxzvZ70zy3iSfTXI8yWNprrB1fZo7qh+EAjKC0iYvJnYyMmWOpW8GoEsBoTYKyAhKm7yY2MnIlDmWvhmALgWE2iggIyht8mJiJyNT5lj6ZgC6FBBqo4CMoLTJi4mdjEyZY+mbAehSQKiNAjKC0iYvJnYyMmWOpW8GoEsBoTYKyAhKm7yY2MnIlDmWvhmALgWE2iggIyht8mJiJyNT5lj6ZgC6FBBqo4CMoLTJi4mdjEyZY+mbAehSQKiNAjKC0iYvJnYyMmWOpW8GoEsBoTYKyAhKm7yY2MnIlDmWvhmALgWE2iggIyht8mJiJyNT5lj6ZgC6FBBqo4CMoLTJi4mdjEyZY+mbAehSQKiNAjKC0iYvJnYyMmWOpW8GoEsBoTYKyAhKm7yY2MnIlDmWvhmALgWE2iggIyht8mJiJyNT5lj6ZgC6FBBqo4CMoLTJi4mdjEyZY+mbAehSQKiNAjKC0iYvJnYyMmWOpW8GoEsBoTYKyAhKm7yY2MnIlDmWvhmALgWE2iggIyht8mJiJyNT5lj6ZgC6FBBqo4CMoLTJi4mdjEyZY+mbAehSQKiNAjKC0iYvJnYyMmWOpW8GoEsBoTYKyAhKm7yY2MnIlDmWvhmALgWE2iggIyht8mJiJyNT5lj6ZgC6FBBqo4CMoLTJi4mdjEyZY+mbAehSQKiNAjKC0iYvJnYyMmWOpW8GoEsBoTYKyAhKm7yY2MnIlDmWvhmALgWE2iggIyht8mJiJyNT5lj6ZgC6FBBqo4CMoLTJi4mdjEyZY+mbAehSQKiNAjKC0iYvJnYyMmWOpW8GoEsBoTYKyAhKm7yY2MnIlDmWvhmALgWE2iggIyht8mJiJyNT5lj6ZgC6FBBqo4CMoLTJi4mdjEyZY+mbAehSQKiNAjKC0iYvJnYyMmWOpW8GoEsBoTYKyAhKm7yY2MnIlDmWvhmALgWE2iggIyht8mJiJyNT5lj6ZgC6FJDtOi3Js5P8+XZ5druO4RSQEZQ2eTGxk5Epcyx9MwBdCsjmfUuS1yf5WJLjSZ6YW463j70uybduaYxTpoCMoLTJi4mdjEyZY+mbAehSQDbnTyR5R5LfT1M0Pp3kl5JcneQ1aQrH1Ul+Ocm9beb3k/yD9rn0o4CMoLTJi4mdjEyZY+mbAehSQDbnkSSfS/IzSZ7ZI//MJD/bPueRNY6rNgrICEqbvJjYyciUOZa+GYAuBWRzXpnkKQOe95Qkl408lpopICMobfJiYicjU+ZY+mYAuhQQaqOAjKC0yYuJnYxMmWPpmwHoUkCojQIygtImLyZ2MjJljqVvBqBLAdmupyR5SZKfSvJ3k7xlbvm72xvaZCkgIyht8mJiJyNT5lj6ZgC6FJDtuSDJF9Jc7eobS5Yntja66VJARlDa5MXETkamzLH0zQB0KSDb838n+Z0k/12SP77lsXSdluSqJJ9P8lj761Xpf4PE5yW5LclX2+W2dt1Bs30pIPv4zIPHZtfdcf/s8K33zq674/7ZfQ8dW5grbfJiYicjU+ZY+mT6HneA3aGAbM8fJrly24NY4F1pdoh3J3lV++ssyTt7PPc5SX4vydEkb2qXB9JcRvg5B8iuQgFZ4IGHj89eceNds3Ou/NBJyytuvGv2wMPHn5QvafKyyUxJY5GRmfp+uupxB9gdCsj2/Lskf2vbg5hzfpqvft0wt/6Gdv35p3j++9Lcyf3szrqz23W/eoDsKhSQOUcfPj678K23L5wE7C0XvvX2J00GSpi8bCNT0lhkZKa8nw457gC7QwHZnjcn+Vfp/9WmTbg2zc5w7tz6c9v11+7z3Ken+crWexY89p4kX2szq2ZXpYDMWfYO5KJ3JPdse/KyrUxJY5GRmfJ+OuS4A+wOBWR7Lk7yr5P8mzRfdfrBJN+/YNmk25N8acljX07ykX2e++I0O9JrFzz2uvaxFw3IrkoB6fjMg8d6TQL2lr3vZm978rKtTEljkZGZ6n469LgD7A4FZHsWXfGqu2zjKlifTnL3ksfuaR9f5pI0O9JLFzz2svaxSwZk93NWmh23u1yaDRaQD/7oT8zu+J6LZ5/4vpfMjv6Vly9cPvF9L9la5qNdMqsAACAASURBVDd+4IdmH3zhX+q9/MYP/NDs6F95+eyO77l49sEX/qXZHd9z8dI/q8ZMSWORkZnqfjr0uLPNY6XM+JmSxiLTP/PBH/2JjcyfFJDt+es9l036QpK7ljx2V5orYi3zyjQ70sULHru4feyyAdn9HG6zJy2bKiB3fM/Fs/vOe67FYrFYLBbL5Jc7vufijcyfFBC6fAKyopt/+Me9E1lRpqSxyMhMdT/1yatMaWOR6Z+5+Yd/fCPzJwWELueArMh3sQGezLlnMqWNRWaczJgUkM35Gxl2xavTkvzEyGNZ5m1Z31WwHkv/q2B1s6tSQFyNBtgyV9+TKWksMgrILvudNOdY/O0k39Uj/8wkP9s+5+E1jqvr+dn/PiAXtL8/I8lz03wFquuWNPfxeGZn3d69PW45QHYVCshc5gHX4wc2bMhxZ9vHSplxMyWNRUYB2WXfmuTvp5lg/1GSTyX55SR/J8nr2+XvpLnz+KfbzKNJfqF97qbclBN3Qr8iJ+6EflMnc6hdd2TuueclOZbm7uZvbJej7brzDpBdhQKyIOOOxMCmrXrcKeFYKTNepqSxyCggJP95mvMc7khTMOYvx/vV9rHXttlNOz3J1WnKwOPtr1e36/ccyuICkjSfknwkzd/t0SQfzolPTg6S7UsB2Sdz30PHZtffef/s8K33zq6/837nfABr1/e4U9KxUkYBkVkvBWS7TkszmX9BuxxK8pQtjqcGCsgGDyAAYyntWCmjgMisjwKyPRf2yPy1tY+iPgqIAgJMUGnHShkFRGZ9FJDteSzNCemLfEuS92bzd0KvgQKigAATVNqxUkYBkVkfBWR7/nmacz4+lidfFeu/SfJAmoLyc1sY19QpIAoIMEGlHStlFBCZ9VFAtuuvp7nq039Mcwfvv5fm6lf/Nsmf3eK4pkwBUUCACSrtWCmjgMisjwKyfYeS/D9pvm71RJJfTHLmNgc0cQqIAgJMUGnHShkFRGZ9FJDte12S309z+d1vJPnNJM/e6oimTQFRQIAJKu1YKaOAyKyPArI9fzLJr6UpHbe2v//RNHdMfzTNTQBZnQKigAATVNqxUkYBkVkfBWR7vpzmk4/Xz63/riR3pvk61vs3PagKKCAKCDBBpR0rZRQQmfVRQLbnniTP3efxn07yhxsaS00UEAUEmKDSjpUyCojM+igg23NGj8zz1j6K+iggCggwQaUdK2UUEJn1UUCojQKigAATVNqxUkYBkVkfBYTaKCAKCDBBpR0rZRQQmfVRQKiNAqKAABNU2rFSRgGRWR8FhNooIAoIMEGlHStlFBCZ9VFAqI0CooAAE1TasVJGAZFZHwWE2iggCggwQaUdK2UUEJn1UUCojQKigAATVNqxUkYBkVkfBYTaKCAKCDBBpR0rZRQQmfVRQKiNAqKAABNU2rFSRgGRWR8FhNooIAoIMEGlHStlFBCZ9VFAqI0CooAAE1TasVJGAZFZHwWE2iggCggwQaUdK2UUEJn1UUCojQKigAATVNqxUkYBkVkfBYTaKCAKCDBBpR0rZRQQmfVRQKiNAqKAABNU2rFSRgGRWR8FhNooIAoIMEGlHStlFBCZ9VFAqI0CooAAE1TasVJGAZFZHwWE2iggCggwQaUdK2UUEJn1UUCojQKigAATVNqxUkYBkVkfBYTaKCAKCDBBpR0rZRQQmfVRQKiNAqKAABNU2rFSRgGRWR8FhNooIAoIMEGlHStlFBCZ9VFAqI0CooAAE1TasVJGAZFZHwWE2iggCggwQaUdK2UUEJn1UUCojQKigAATVNqxUkYBkVkfBYTaKCAKCDBBpR0rZRQQmfVRQKiNAqKAABNU2rFSRgGRWR8FhNooIAoIMEGlHStlFBCZ9VFAqI0CooAAE1TasVJGAZFZHwWE2iggCggwQaUdK2UUEJn1UUCojQKigAATVNqxUkYBkVkfBYTaKCAKCDBBpR0rZRQQmfVRQFjkeUluS/LVdrmtXdfHn0vyD5N8KsmjSb6U5GNJfmBB9lCanW/R8ksDx66AKCDABJV2rJRRQGTWRwFh3nOS/F6So0ne1C4PJHmkfexU3pfk4SQ3JXl1kp9O8uk0O9nr5rKH2vUfSHLZ3PKigeNXQBQQYIJKO1bKKCAy66OAMO99SY4nObuz7ux23a/2eP5fSPLUuXVnJrk/ye8mOb2z/lCane+agWNdRAFRQIAJKu1YKaOAyKyPAkLX05M8luQ9Cx57T5KvtZkh3pFmR+sWm0M5UUDObJeDUkAUEGCCSjtWyiggMuujgND14jQ7w2sXPPa69rGhX436lSRfT/LNnXWH2td8NCfO/fhckp8a+GckCogCAkxSacdKGQVEZn0UELouSbMzvHTBYy9rH7tkwOv+mTSfrNwyt/5ZST6a5A3tn/maJL/V/jlv7/G6Z6XZcbvLpVFAFBBgcko7VsooIDLro4DU64wkz+25PKt9zivT7AwXL3i9i9vHLltxHN+a5L40J7E/6xTZJDktySeSPJHk2afIHs6Sq2gpIAoIMC2lHStlFBCZ9VFA6nUoyy9xO798vH3O2J+AnJnkk0n+MMn3rvC8l7Z/1qtPkfMJyIAMQIlKO1bKKCAy66OA1OvMJC/pubywfc6Y54B8U5KPJHk8yQ+vOPYL2j/r6hWflzgHRAEBJqm0Y6WMAiKzPgoIXae6CtZj6XcVrNPT3NvjiSR/bcA4Xp5mp7xiwHMVEAUEmKDSjpUyCojM+iggzLslzT0/ntlZt3cfkPmTyM9Kcw7JGZ11T0nyz9LsVD95ij/r2xasOzPJPWk+OTl7weOnooAoIMAElXaslFFAZNZHAWHeeUmOpbkT+hvb5Wi77ry57JE0O8+hzrrrcuK8kvm7m1+W5Bmd7PvTfE3rcJJXJXlLki+0z79q4PgVEAUEmKDSjpUyCojM+iggLHJBmmLwaLt8uF0370hOLiAfz/4nvH9fJ3tFmpPUv5zmHiGPJPlYmhPeh1JAFBBggko7VsooIDLro4BQGwVEAQEmqLRjpYwCIrM+Cgi1UUAUEGCCSjtWyiggMuujgFAbBUQBASaotGOljAIisz4KCLVRQBQQYIJKO1bKKCAy66OAUBsFRAEBJqi0Y6WMAiKzPgoItVFAFBBggko7VsooIDLro4BQGwVEAQEmqLRjpYwCIrM+Cgi1UUAUEGCCSjtWyiggMuujgFAbBUQBASaotGOljAIisz4KCLVRQBQQYIJKO1bKKCAy66OAUBsFRAEBJqi0Y6WMAiKzPgoItVFAFBBggko7VsooIDLro4BQGwVEAQEmqLRjpYwCIrM+Cgi1UUAUEGCCSjtWyiggMuujgFAbBUQBASaotGOljAIisz4KCLVRQBQQYIJKO1bKKCAy66OAUBsFRAEBJqi0Y6WMAiKzPgoItdnZAvKZB4/NXnDNnbNzrvzQ7AXX3Dm776Fj/f8iAFvU9/hV0jFXRgHZtcyYFBBqs3MF5IGHj89eceNds3Ou/NBJyytuvGv2wMPHV/+LAWzAqsevEo65Mv0yJY1FRgGBddupAnL04eOzC996+8J/vPeWC996uxICFGfI8Wvbx1yZ/pmSxiKjgMC67VQBWfbO4aJ3EgFKMuT4te1jrkz/TEljkVFAYN12poB85sFjvf7x3lucEwKUYujxq7QJmYwCskuZMSkg1GZnCsh1d9y/0j/g1995//C/KMCIhh6/SpuQySggu5QZkwJCbXamgBy+9d6V/gE/fOu9w/+iACMaevwqbUImo4DsUmZMCgi12ZkC4hMQYKp8AlJ/pqSxyCggsG47U0CcAwJMlXNA6s+UNBYZBQTWbWcKyGzmKljAdLkKVt2ZksYio4DAuu1UAXnAfUCAiRpy/Nr2MVemf6akscgoILBuO1VAZjN3Qgemy53Q682UNBYZBQTWbecKyJ77Hjo2u/7O+2eHb713dv2d9zvnA5iMvsevko65MgrIrmXGpIBQm50tIAC1K+2YK6OA7FJmTAoItVFAACpV2jFXRgHZpcyYFBBqo4AAVKq0Y66MArJLmTEpINRGAQGoVGnHXBkFZJcyY1JAqI0CAlCp0o65MgrILmXGpIBQGwUEoFKlHXNlFJBdyoxJAaE2CghApUo75sooILuUGZMCQm0UEIBKlXbMlVFAdikzJgWE2iggAJUq7Zgro4DsUmZMCgi1UUAAKlXaMVdGAdmlzJgUEGqjgABUqrRjrowCskuZMSkg1EYBAahUacdcGQVklzJjUkCojQICUKnSjrkyCsguZcakgFAbBQSgUqUdc2UUkF3KjEkBYZHnJbktyVfb5bZ2XR+H0uxQi5ZfWpA/LclVST6f5LH216va9UMoIACVKu2YK6OA7FJmTAoI856T5PeSHE3ypnZ5IMkj7WOncijNDvWBJJfNLS9akH9Xm393kle1v86SvHPg+BUQgEqVdsyVUUB2KTMmBYR570tyPMnZnXVnt+t+tcfzD6XZoa7pkT0/yTeS3DC3/oZ2/fk9XmOeAgJQqdKOuTIKyC5lxqSA0PX0NF+Des+Cx96T5GttZj+HcqKAnNkuy1zbZs+dW39uu/7aU474ZAoIQKVKO+bKKCC7lBmTAkLXi9PsDK9d8Njr2scWfY2q61CbezQnzv34XJKfWpC9PcmXlrzOl5N85JQjPpkCAlCp0o65MgrILmXGpIDQdUmaneGlCx57WfvYJad4jWcl+WiSN7Sv85okv9U+9+1z2U8nuXvJ69zTPr6fs9LsuN3l0iggAFUq7Zgro4DsUmZMCki9zkjy3J7Ls9rnvDLNznDxgte7uH3ssgFjOS3JJ5I8keTZnfVfSHLXkufcleaKWPs5nCVX3FJAAOpT2jFXRgHZpcyYFJB6Hcryy+HOLx9vnzPGJyDLvLR9/qs763wCMnIGoGalHXNlFJBdyoxJAanXmUle0nN5YfucMc4BWeaC9vlXd9Y5B2TkDEDNSjvmyiggu5QZkwJC16mugvVYTn0VrGVenmZHu6Kz7m1xFaxRMwA1K+2YK6OA7FJmTAoI825Jc8+PZ3bW7d0H5Ja57FlpziE5o7Pu2xa85plpvlL1eJ58f5HnZ//7gFyw4tgTBQSgWqUdc2UUkF3KjEkBYd55SY6luRP6G9vlaLvuvLnskTQ7z6HOuven+erU4TR3Nn9LmpPNZ0muWvDn3ZQTd0K/IifuhH7TwPErIACVKu2YK6OA7FJmTAoIi1yQpkQ82i4fzuJPI47k5AJyRZJPpjmH4+tJHknysTQnsS9yeprzQo6m+YTkaPv70weOXQEBqFRpx1wZBWSXMmNSQKiNAgJQqdKOuTIKyC5lxqSAUBsFBKBSpR1zZRSQXcqMSQGhNgoIQKVKO+bKKCC7lBmTAkJtFBCASpV2zJVRQHYpMyYFhNooIACVKu2YK6OA7FJmTAoItVFAACpV2jFXRgHZpcyYFBBqo4AAVKq0Y66MArJLmTEpINRGAQGoVGnHXBkFZJcyY1JAqI0CAlCp0o65MgrILmXGpIBQGwUEoFKlHXNlFJBdyoxJAaE2CghApUo75sooILuUGZMCQm0UEIBKlXbMlVFAdikzJgWE2iggAJUq7Zgro4DsUmZMCgi1UUAAKlXaMVdGAdmlzJgUEGqjgABUqrRjrowCskuZMSkg1EYBAahUacdcGQVklzJjUkCojQICUKnSjrkyCsguZcakgFAbBQSgUqUdc2UUkF3KjEkBoTYKCEClSjvmyiggu5QZkwJCbRQQgEqVdsyVUUB2KTMmBYTaKCAAlSrtmCujgOxSZkwKCLVRQAAqVdoxV0YB2aXMmBQQaqOAAFSqtGOujAKyS5kxKSDURgEBqFRpx1wZBWSXMmNSQKiNAgJQqdKOuTIKyC5lxqSAUBsFBKBSpR1zZRSQXcqMSQGhNgoIQKVKO+bKKCC7lBmTAkJtFBCASpV2zJVRQHYpMyYFhNooIACVKu2YK6OA7FJmTAoItamygHzmwWOzF1xz5+ycKz80e8E1d87ue+jYgccOMCV9j4OlTdp2NVPSWGQUEFi3qgrIAw8fn73ixrtm51z5oZOWV9x41+yBh4+P9ncBKNGqx8HSJm27milpLDIKCKxbNQXk6MPHZxe+9faF/+juLRe+9XYlBKjWkONgaZO2Xc2UNBYZBQTWrZoCsuwdv0XvAALUaMhxsLRJ26kytX61rKSx1JrZ9L4zJgWE2lRRQD7z4LFe/+juLc4JAWoz9DhY0gRxv0ztXy0raSy1Zba174xJAaE2VRSQ6+64f6V/eK+/8/7R/24A2zT0OFjCBPFUmV34allJY6kps819Z0wKCLWpooAcvvXelf7hPXzrvaP/3QC2aehxcNsTxD6ZXfhqWUljqSmzzX1nTAoItamigPgEBNh1tX4CUvtXy0ocSy2Zbe87Y1JAqE0VBcQ5IMCu2/Zkq7Q3mEr/e63jNWo9QX9oZtv7zpgUEGpTRQGZzVwFC2DqX1VaNIEe8tWyKU7ED/IatZ+g3ydT4r4zJgWE2lRTQB5wHxBgxw05DpYwidxvAv3it310pUnkRT+/OF/6RHzoa+zCCfr7ZUred8akgFCbagrIbOZO6ABTeze8zwR6rKXkifgmP/0v6e99kEzp+86YFBBqU1UB2XPfQ8dm1995/+zwrffOrr/zfud8ADun73Fw25PIvhPosZZSJ+JDXqPW8376Zkrfd8akgFCbKgsIAP1scxK56gR6rKXEifipHl90bsKQk6yneH7MoswU9p0xKSAs8rwktyX5arvc1q7r40iaHWrZ8nc62UP75H5p4NgVEIAdts1J5KoT6GXf0V/1u/4lXilr2eP7faXuhW2R6Lu88NrF+RK/mneqzBT2nTEpIMx7TpLfS3I0yZva5YEkj7SPncqLk1y2YPmtnLyjHWrXfWBB/kUDx6+AAOywbU4ih1ylaNFXy2q4CeOix0s/x8G+o4CwPe9LcjzJ2Z11Z7frfnXgaz4tzScpn5pbfyjNznfNwNddRAEB2GFTehd72U1kt32/h3UVkNLPcbDvKCBsx9OTPJbkPQsee0+Sr7WZVf33aXayn5lbfygnCsiZ7XJQCgjADpvS9/iXna9Qw8nYB902Yy0lbpsxts829p0xKSB0vTjNzvDaBY+9rn1syFej7kjy9STPmFt/qH3NR3Pi3I/PJfmpAX/GHgUEYIdte6I51k1kp3452oO+wz/WUuKnQ8sype87Y1JA6Lokzc7w0gWPvax97JIVX/O7kjyR5EMLHntWko8meUP7Z74mJ84VeXuP1z4rzY7bXS6NAgKws7Y9iRzrJrJTvQnjssdXPTdh2Qnmq56oXuL5Mcsype87Y1JA6nVGkuf2XJ7VPueVaXaGixe83sXtY5etOI6r2uf9aM/8aUk+kaa0PPsU2cNZchUtBQRgN217EjmbjXcT2andhHG/x4ecm7DoJOsazo/ZL1PyvjMmBaReh7L/5XC7y8fb56zjE5D7kvxukqeu8JyXtn/Wq0+R8wkIAE9SwiRyz1g3ke3zOqXdD2P+8Smc47CpTJ//V5vcd/r+vcakgNTrzCQv6bm8sH3O2OeAvKB9zjtXHPsF7fOuXvF5iXNAAHZaSRPNTSj1U5IpnuMwlU831kEBYZtOdRWsx7LaVbD+cZqd64WnCs55efu8K1Z8XqKAAOy0Eiaam9LnvhrbOE9k2Tv82zrHoYRPh4b8v9okBYRtuyXNPT+e2Vm3dx+QW+ayZ6U5h+SMBa/zTUl+J81XsJb5tgXrzkxyT5LH8+R7kfSlgADssF0qIKV9EtDnHf5NnuNQ0qdDY336sy4KCNt2XpJjae6E/sZ2OdquO28ueyTNznNowevsfYpx5T5/1vuTfCTNyeSvSvKWJF9on3fVsOErIAC7bFcKSGnnQqz6Dv+6z3Eo6dOhsc5/WScFhBJckKYYPNouH27XzTuS5QXkg2muZPWn9/lzrkjyySRfTnOfkEeSfCzNCe9DKSAAO2xTx+W+X+tZl9KuBlXaO/wlfTo01l3O12Ub+7ICQm0UEIAdtisnEq96X4113g+jtHf4S/t0aOj/q3Xb5r6sgFAbBQRgh+3KicQlfQJS2jv8JW2bg4xnnba9Lysg1EYBAdhhu3IicUnv8pf2Dn9Jnw7NZmWeA7LtfVkBoTYKCMAOM4ncfxLpE5DNfwIym21/wt9Vwr6sgFAbBQRgh+3S12iG3FfDOSDLx7POf9PHugfKGErYlxUQaqOAAOwwJxI/+d30Td3roqR3+IeOZ93/pk/9AgZjUkCojQICsMN26ROQrj731RjrjuCLXqekd/hns+3dLb3P64x1D5ShStiXFRBqo4AA7LBdOgekr7E+JTnV6/z6Z79SxDv8fcc75t3SS/l0o48S9mUFhNooIAA7bFdOJO5rrDuCr/I6236Hf96675a+7UvaDrHtfVkBoTYKCMAO25UTifsa61yIbU9Y12GXt82292UFhNooIAA7bFdOJO5jyFdtFp3DUMJXdsZm27gTOoxJAQHYYZs6Lpf2NaNFVj3Z+KKf/+jC9S9+2+L1y5ZNn4A/hG1zwjb2ZQWE2iggADvMcfmEVS+3OtayqUsQH4Rts10KCLVRQAB2mOPyCau+yz/WUuK7/PNsm+1SQKiNAgKwwxyXT1j1/ISxlhK/jjbPttkuBYTaTK6A9L35EQD7G+tmcn0zU9D3Ck1jLSVd6elUbJvtUUCozWQKyJSupAJQsrFuJrdqZgr6XG51rKW0SxCfim2zPQoItZlEAZniTYsASjTWzeTm1VJAZrP9C9qqV3FadjWoqb5xZttshwJCbSZRQKZ40yKAEo11M7l5NRWQPYsutzrkPhZTuATxqmybzVJAqM3GCsjQ7xpP+aZFACUZejx1bt6TeVNsOdtmPRQQarP2AnLQ7xqveuk/l+wDWGzo8dS5eU/W51yIXf1asG2zHgoItVlrARnju8ar3vzITYsAFht6PHVu3sl2tXz1YduMTwGhNmstIGN819gnIADjGPMTEF+1aTiHYTnbZjwKCLVZWwEZ67vGzgEBGIfjMkyTAkJt1lZAvNMGUB6fTMP0KCDUZm0FZMzvGjupDWAcQ46nzs2D7VJAqM0kPgGZzZzUBjAWVyeEaVFAqE3x54DMc1IbwDj6Hk+dAwLbpYBQm6KugrVrN7MCKN2y47Jz82BzFBBqs9YC0ve7xr/+2a/4ehVAQU71Na1f/+xXnJsHG6KAUJut3wndP2IAZel7k0FvHsFmKCDUZu0FZM+y7xr7GB+gLKsel52bB+ulgFCbjRWQRZzICFAWx2UojwJCbbZaQFzKEaAsjstQHgWE2my1gLiZFUBZHJehPAoItfEJCAD/ieMylEcBoTbOAQHgP3FchvIoINRmqwVkNnMVLIDSOC5DWRQQarP1AtL3ZoWuJw+wGY7LUBYFhNpsvYDMZqe+WaF/5AA2y3EZyqGAUJsiCsgeN7MCKIvjMmyfAkJtiiogAAA8mQJCbRQQAICCKSDURgEBACiYAsK8i5P8b0nuTvJ4mp3j0IqvcVqSq5J8Pslj7a9XtesPku1DAQEAKJgCwrwjaYrAv0ny6QwrIO9qn/fuJK9qf50leecBs30oIAAABVNAmPddSZ7a/vc1Wb2AnJ/kG0lumFt/Q7v+/IHZvhQQAICCKSDsZ0gBubZ9zrlz689t1187MNuXAgIAUDAFhP0MKSC3J/nSkse+nOQjA7N9KSAAAAVTQNjPkALy6TQnsC9yT/v4kOwiZ6XZcbvLpVFAAACKpYDU64wkz+25PGvJawwpIF9IcteSx+5Kc5WrIdlFDqcZ30mLAgIAUCYFpF6HsmRyvmD5+JLX8AkIAACjUkDqdWaSl/RcXrjkNZwDAgDAqBQQ9jOkgLwt/a9stUq2LwUEAKBgCgj7OVUBOSvNOSRndNY9P/vf2+OCgdm+Lkoyu/nmm2d33323xWKxWCwWi6Ww5eabb94rIBcNmOtRoQuSvLldPp5m5/gH7e/fMJc9ksUF5aacuLv5FTlxd/ObFvx5q2T7uDT9z32xWCwWi8VisWxvuTSQ5PIs30m+OJc90q4/NLf+9CRXJzma5PH216vb9fNWyfbx7Wl25oty8gnq61j2Cs+lG/rzdnWxnW3jWhbb2TauZbGdbeODLBe1f69vD7Cy705zcPjubQ+kcrbz+tnGm2E7r59tvBm28/rZxsBCDg6bYTuvn228Gbbz+tnGm2E7r59tDCzk4LAZtvP62cabYTuvn228Gbbz+tnGwEIODpthO6+fbbwZtvP62cabYTuvn20MLHRWksPtr6yP7bx+tvFm2M7rZxtvhu28frYxAAAAAAAAAAAAAAAAAAAAAAAAAFt2WpKrknw+yWPtr1e161nd09NcavDXkvyHNNc9P7Ika9uv7s8l+YdJPpXk0SRfSvKxJD+wIGv7DvdnkvyzJJ9LcjzJV5P8v0nemOSpc1nbeTzfn+aYMUvyX849dmaSv5fkt5N8Lcm9SX5yo6ObpkM5sU3nl1+ay9qXD+asJDcm+Xdptt9DST6Q5BlzuZ9Ms/9+Lc3+/PfS7N/ADnlXmgPxu5O8qv11luSd2xzUhB1Ks/0eSlNC9isgtv3q3pfk4SQ3JXl1kp9O8uk02+11c1nbd7iLk9yR5H9J8pokr0/y3iTfSLNfd9nO4/imJP9fmsK3qID8iyRfT3J9mgncB9rcz25wjFN0KM12+kCSy+aWF81l7cvDPSfNm25Hk7wlyU8kuTLNdn9OJ/c/5sT/j59Msz//UZIPbXKwwHadn2ZCccPc+hva9edvfETT99Qk39X+9+lZXkBs+2H+Qk5+B/7MJPcn+d002zyxfdflH6fZp89rf287j+fnknw5zYRsvoD8cLvuTXPP+WCSP0jyJzcxwIk6lGbbXXOKnH15uP8syb9OcneabwEs8yfT7K8fnFv/pjT/j35oLaMDinNtmh/6c+fWn9uuv3bjI6rLfgXEth/XO9Jst7Pb39u+6/G302y/P9/+3nYexzlJP2QAAgAABUdJREFUfj/J30jzFc75AvLeJH+Yk7+m8hfbrK9iLXcoJwrImVn+VR/78nB7Xx384fb3fyzNJ3rzfrLN/cW59U9Ls3+/d10DBMpye5rv0C/y5SQf2eBYarRfAbHtx/Urab6e8s3t723fcTwtyXekmSC/IslXkjyYE5M423kcH0zyG2neST6ckwvI/Ul+c8HzzmyzN615fFN2KM02ejQnzv34XJKfmsvZl4f7+2m26/cm+Zftf38jzT794k7uF9vHFpXAf5Xk3653mEApPp3mI9NF7mkfZ7j9CohtP54/k+aEx1s662zfcRzOk0/a/c0kF3Qet50P7oeTPJHkz7a/P5yTC8ijefL+3fW7Ofm8HE54VpKPJnlDkpemOafpt9Js47d3cvbl4fbOR3o4yfvTvFnxt5L8xzSfbOx9fe3X2nWLvD/Nfg7sgC8kuWvJY3eluQIIw+1XQGz7cXxrkvuSPJJmorHH9h3Hf5HmCmM/lubdy99I8n2dx23ngzkzzUm7v9hZdzgnF5AnkvzTJa/xUJoJNv2dluQTabbrs9t19uXhPppmn53fD//rdv0/b3//sTT76yL/NM3J6MAO8I7PevkEZL3OTPLJNO+wfe/cY7bverwpyeNpPnVKbOeDuibNJxjf0Vl3OD4B2YSXptnOr25/b18ebu+Kj399wWNfTPMVtr2cT0AA33ldM+eArM83pdlGj+fEiY9dtu96PCNPPiHXdh7uT6e5D8IvpCkbe8s/SrONvz/NJ1CJc0DW4YI02+7q9vf25eFuSrMtX7Lgsd9Mc5xOnAMCtN4WV/1Yp/0KiG0/3OlpvnP8RJK/tiRj+67HoTTb713t723n4S7M8hvk7S3H2+w/jatgje3labbdFe3v7cvDXZFmG71qwWP/vl2S5tOmRVfBOjOuggU75fnZ/7rnF5z0DFaxXwGx7Yd5Spo7dJ9q0mX7HsyfWrJ+72o3e1+1sJ2H+9YkP7pg+edptvEbkvxIm937utCy+4As+/9F8m0L1p2Z5mtVj+fEpbvty8N9R5rLSP9mnnzX+L371+zdcf5PpSkaH5h7/t59QF663mECJdn76PTdad7F2Lvzq4/0h3tDkjenuRvsLM0/dG9ul+4/Yrb96q5Ls40+npPvanxZmq8I7bF9h/s/05xf8z+nKXo/mxMnmv56TtzwMbGdx3Y4i++E/pE0l5q+Ls07zXtXHrpqk4OboPen2XaH02y3t6Q54XzRtrMvD/e30myrT6b5N/Dn05SSr+REyUuabT5Lc4x5VZr9+etJPrzJwQLbd3qa78AeTfNu0NH296fv9yT29cUs/0rF5Z2cbb+6j2f/r6x8Xydr+w73Y2kmBA+m2XaPprl06d/OyXeit53HdTiLC8jT0lw29t+nuez0fWkuKcv+rkgzKf5ymonuI2muxvSyBVn78sFcluYNt6+luTjC/5ET5zF1vSbN/vtYmv357Wn2bwAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAFi3y5PMkvzlJG9O8ttJvpbkXyZ5/vaGBQAA1OjyNAXkt9rlTUl+Lsl/TPJAktO3NjIAAKA6l6cpIPckOaOz/kfa9f/tFsYEAABU6vI0ReM1c+v/RLv+f9j0gAAAgHpdnhPngMybJfmfNjoaAACgapenKRo/sOCxWZLDmxwMAABQt8ujgAAAABtyeRQQAABgQy6PAgIAAGzI5VFAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAACA1f3/BLDoxIo/ZGkAAAAASUVORK5CYII=\" width=\"640\">"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/plain": [
       "Text(0,0.5,'x(n)')"
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "plt.figure()\n",
    "plt.stem(n, x)\n",
    "plt.xlabel('n')\n",
    "plt.ylabel('x(n)')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "U slučaju da su početne faze harmonika oblika\n",
    "\n",
    "$$\\phi_i = -\\frac{2 i \\pi}{3}$$\n",
    "\n",
    "situacija odgovara prolasku signala kroz sistem sa nelinearnom fazom. Talasni oblik signala je promijenjen iako je njegov amplitudni spektar isti. Sistem je unio linearna izobličenja signala."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [],
   "source": [
    "phi0 = 0\n",
    "phi1 = -2*np.pi*1/3\n",
    "phi2 = -2*np.pi*2/64\n",
    "phi3 = -2*np.pi*3/64\n",
    "n = np.arange(64)\n",
    "x = np.sin(2*np.pi/64 * n + phi0) + (1/3)*np.sin(2*np.pi/64 * 3 * n + phi1) + \\\n",
    "    (1/5)*np.sin(2*np.pi/64 * 5 * n + phi2) + (1/7)*np.sin(2*np.pi/64 * 7 * n + phi3)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {
    "scrolled": false
   },
   "outputs": [
    {
     "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,iVBORw0KGgoAAAANSUhEUgAAAyAAAAJYCAYAAACadoJwAAAgAElEQVR4nO3df7CtV10n6M9A0kyU0Wl/tJ3hRy7DUNBVElP2XFsyzrTttBlqFKdHpOxKRU0LiiBTQo9OME3jYeYC2ui9gzUIdtlUqA7OjAoD14jJPViFMJ22W8NMyUmmD8q9VNs3CmkbL0Ql2PDOH+8+3Jtz9z5n73e/e7/rXe/zVK26OXt/977rrLOTrM9537VWAgAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAALCWJyfZSfKrSf4wSZPkrhVef9fsNfPaU/vrJgAAUIMTacPCw2lDSNcA8r1JbjvUvqS/bgIAADV4UpKnzP75mnQPINf02isAAKB66wSQa5N8WZIn9N4rAACgSusEkE/P/vyzJO9K8sw1+vGVSW5NcnOSr9c0TdM0TdOKazenna99ZWANXQLITyb5mbQfwBcmeVOSzyZ5JMkNS7z++lz9gb4zixe2a5qmaZqmaeW0WwNr6BJA5vnWFd5nJws+0HfffXfzwAMPaJqmaZqmaYW1u++++2DOdvOa80Ymrq8AkiQPpN1Z6zjzroDcmqR54IEHGgAAyvPAAw8cBJCv72HeyIT1GUDeneRzHV/79RFAAACKJYDQlz4DyO8m+YOOrxVAAAAKJoDQl6MCyLVJnpP2lqkDX5r2HJHDvnv2Pm/r2A8BBACgYAII63pFktckeW3aD9KHZ1+/JsmNs5oTuTqc3JTkYpK3JHllkpcn+adJPp/k40m+pmN/BBAAgIIJIKzr41m8tdrts5oTuTqA/NUk70zy0SSPJnksye8nOZPkq9bojwACAFAwAYTaCCAAAAUTQKiNAAIAUDABhNoIIAAABRNAqI0AAgBQMAGE2gggAAAFE0CojQACAFAwAYTaCCAAAAUTQKiNAEInD1681Jw+t9/snN1rTp/bbx56+FKnGgDgaAIItRFAWMmFRx5tXvTW+5sb7rjnqvait97fXHjk0aVqAIDlCCDURgBhaecfebS56XX3zQ0WB+25P3Fvc+PO0TU3ve4+IQQAliSAUBsBhKUtuqrRpb3orfcP/e0AwCgIINRGAGEpD1681Fv4OGjWhADA8QQQaiOAsJTT5/Z7DyBndveH/rYAoHgCCLURQFjKztm93gPIztm9ob8tACieAEJtBBCW4goIAAxDAKE2AghLsQYEAIYhgFAbAYSlbWIXLIcVAsDRBBBqI4CwtAs9ngPyoY9+0mGFALAEAYTaCCCspI+T0D/00U8eG2QcVggALQGE2gggdPLQw5eak6d2mxvuuKc5eWp37q1Ti2qWvZXLYYUAIIBQHwGEuZZZm/FtP/vB5oY77mm+7Wc/uPB9DtesupjdmhAApk4AoTYCCI+zzC1WB7oEkFW387VVLwBTJ4BQGwGELzq/xCLzK9dmdAkgqx5o6LBCAKZOAKE2AghftOraDFdAAGDzBBBqI4DQNE23tRnWgADA5gkg1EYAoWmablcmugSQprELFgCsQgChNgIITdN0W5vRNYAsc6Chc0AAoCWAUBsBhKZptnsFpGlW220LAKZMAKE2AghN02xvDchhyxxoCABTJoBQGwGEL9rGLlhdapY5FBEAaiWAUBsBhC9adW3GpgOI27QAQAChPgIIj7Ppk9CXrVn1UEQAqJUAQm0EEOZaZm3GJgOIrXoBoCWAUBsBhIWOCw+bCiAOKwSAywQQaiOAsNBQAaTLlsAAUCsBhNoIICw0VADpcigiANRKAKE2AggLuQICAMMTQKiNAMJC1oAAwPAEEGojgLDQUAGkaeyCBQAHBBBqI4Cw0JABZNVDEQGgVgIItRFAWGjIANI0TkIHgKYRQKiPAMJCQweQA8sciggAtRJAqI0AwkKlBJBlawCgRgIItRFAWEgAAYDhCSDURgCZoAcvXmpOn9tvds7uNafP7S+8pWlMAWTZ7wkAxkYAoTYCyISsuqh7DAHEQnUAaieAUBsBZCLOd9jWtvQA0uV7AoCxEUCojQAyEV0O9is9gDisEIApEECojQAyAQ9evLTURP2gHayfKDmAdP2eAGBsBBBqI4BMwOlz+ytN1s/s7jdNU3YA6fo9AcDYCCDURgCZgJ2zeytN1nfO7jVNU3YA6fo9AcDYCCDURgCZAFdAXAEBYLwEEGojgEyANSDWgAAwXgIItRFAJsIuWAAwTgIItRFAJuJCheeAdPmeAGBsBBBqI4BMiJPQAWB8BBBqI4BM0EMPX2pOntptbrjjnubkqd2F6yPGEEBW/Z4AYGwEEGojgExUH5P+kgLIsjUAMDYCCLURQCZKAAGAcRBAqI0AMlFTDCAPXrzUnD633+yc3WtOn9t3mxYAoyCAUBsBZKKmFEAsVAdgzAQQaiOATNRUAsh5W/UCMHICCLURQCZqKgHEYYUAjJ0AQm0EkImaQgB58OKlpcLHQbMmBIASCSCs48lJdpL8apI/TPtBumvF9/jaJO9L8ulZe9/ssa4EkImaQgA5fW5/pQByZnd/4XsDwFAEENZxIu2H5+G0IWTVAPKsJH+S5HySV83ahSSfmj3XhQAyUVMIIDtn91YKIDtn9xa+NwAMRQBhHU9K8pTZP1+T1QPIryR5NMnTrnjsabPHfrljnwSQiZpCAHEFBIAaCCD0ZdUA8uQkjyV5x5zn3pHks7OaVQkgEzWFAGINCAA1EEDoy6oB5Hmz+h+a89zLZs99Y4d+CCATNYUA0jR2wQJg/AQQ+rJqAHnhrP4Fc577jtlzLzzmPa5P+8G9st0aAWSSphJALjgHBICRE0Doy6oB5Htm9bfMee6W2XO3HfMeO7O6q5oAMj1TCSBN4yR0AMZNAKEvroAwqCkFkAMPPXypOXlqt7nhjnuak6d2rfkAYBQEEPpiDQiDmmIAWbYGAEoigNCXvnfBeix2wWIFAogAAsA4CCD05agAcm2S56S9ZepK70p75sdTr3js4ByQd3XshwAyUQKIAALAOAggrOsVSV6T5LVpP0gfnn39miQ3zmpOZH44eXaSS2lPQn/lrJ2fPfbsjv0RQCrz4MVLzelz+83O2b3m9Ln9hescBJD5NcuOHwBsiwDCuj6eBTtRJbl9VnMii6+O3Jjk3iSfmbVfz+Xg0oUAUolVd3oSQOyUBcA4CCDURgCpwPkOZ10IIJdruowfAGyLAEJtBJAKdDntWwBxWjoA4yCAUBsBZOQevHhpqcnzQTtY0yCAtDVdxw8AtkUAoTYCyMidPre/0gT6zO5+0zQCyEFN1/EDgG0RQKiNADJyO2f3VppA75zda5pGADmo6Tp+ALAtAgi1EUBGzhWQ9WpcAQGgdAIItRFARs4akPVqrAEBoHQCCLURQCpgF6z1auyCBUDJBBBqI4BU4IJzQNaq6TJ+ALAtAgi1EUAq4SR0J6EDUCcBhNoIIJV56OFLzclTu80Nd9zTnDy1u3DNggAyv2bZ8QOAbRFAqI0AUqFtTdZLCg7brgGAbRFAqI0AUiEBRAABoB4CCLURQCokgAggANRDAKE2AkiFBBABBIB6CCDURgCpkAAigABQDwGE2gggFRJABBAA6iGAUBsBpEICiAACQD0EEGojgFRIABFAAKiHAEJtBJAKCSCbrXnw4qXm9Ln9ZufsXnP63L7DCgHYKAGE2gggFRJANlNz4ZFHmxe99f7mhjvuuaq96K33NxceeXThewJAVwIItRFAKiSA9F9z/pFHm5ted9/c8HHQbnrdfUIIAL0TQKiNAFIhAaT/mkVXPuZdCQGAPgkg1EYAqZAA0m/NgxcvLRU+Dpo1IQD0SQChNgJIhQSQfmtOn9tfKYCc2d1f+N4AsCoBhNoIIBUSQPqt2Tm7t1IA2Tm7t/C9AWBVAgi1EUAqJIC4AgJAPQQQaiOAVEgAsQYEgHoIINRGAKmQAGIXLADqIYBQGwGkQgJI/zUXnAMCwEAEEGojgIzIgxcvNafP7Tc7Z/ea0+f2F97qI4A4CR2Aeggg1EYAGYFVJ74CyGZrHnr4UnPy1G5zwx33NCdP7VrzAcBGCSDURgAp3PkOt/4IIGXUAEAfBBBqI4AUrsviZwGkjBoA6IMAQm0EkIJ13f5VACmjBgD6IIBQGwGkYF0PwBNAyqgBgD4IINRGACnYztm9lQLIztm9pmkEkFJqAKAPAgi1EUAK5grIuGtYz7LbTgPUTgChNgJIwawBGXcN3ThvBeDxBBBqI4AUzi5Y461hdV22nQaonQBCbQSQwl1wDshoa1hdl8ANUDsBhNoIICPgJPRx1rDYvPUdXW85BKidAEJtBJAReejhS83JU7vNDXfc05w8tbtwAiaAlFHD1Y4K0897w/tXCiAHmy4A1E4AoTYCyMiUNMkuqS8l1vB4y6zvWKUdbDsNUDsBhNoIICNT0iS7pL6UVmML2astu77DFRCAxxNAqI0AMjIlTbJL6kspNbaQnW/V9R3WgABcJoBQGwFkZEqYZJfYlxJqbCG72KqHah7X7IIFTIkAQm0EkJEZepJdal9KqLGF7GI7Z/d6Cx9TDXHAdAkg1EYAGZmhJ9ml9mXoGlvIHm3VKyA3v3H+jlhTvo0NmC4BhNoIICNT0kS8pL4MXbPqBHtqC6i7BLRlt50GqJ0AQm0EkJEpaSJeUl+Grln1FqMpbiHb5Ra1ZX4WALUTQKiNADIyJU3ES+rL0DWugBzvQodF+gIIgABCfQSQkSlpIl5SX4ausQZkOatuU3zcz8J5K8AUCCDURgAZmZIm4iX1pYQau2Atb9n1HYt+Fs5bAaZEAKE2AsjIDD3JLrUvJdR0ucVoyrqOs/NWgKkRQKiNAFKIZW8lGXqSXWpfSqnxm/nludIEsBwBhNoIIAPr+574bdaU1JfSamwhezxrbQCWI4BQGwFkQF1uJSlpkl1SX8ZaM2VdxtBuY8AUCSDURgAZ0KbORRBAxlMzZV3G0HkrwBQJINRGABlI11tJSppAl9SXsdZMmSsgAMsRQKiNADKQrhOpkibQJfVlrDU12uSGCtaAAFMkgFAbAWQgXW8lKWkCXVJfxlpTk21tqGAXLGBqBBBqI4AMxBUQNTXZ5oYKzlsBpkYAYR1PTPLqJL+f5LHZn6+ePX6cu9J+8Oa1p67RJwFkINaAqKnJtjdUcN4KMCUCCOv4ubQfnrcnecnszybJW5Z47V2z2u9Nctuh9iVr9EkAGZBdsNTUYMgw7bwVYAoEELp6bpIvJHnzocffPHv8uce8/q60H7xreu6XADKgLreSlDSBLqkvY62pwRhuJwQYMwGErl6f9oPzjEOPP2P2+OuPef1ds7prk3xZkif01C8BZGBOQp92TQ3GsKECwJgJIHR1X5I/WvDcJ5Lce8zr70r7wfv07M8/S/KuJM9cs18CSCGWvZWkpAl0SX0Za00NXAEB2CwBhK4+kuSBBc99ePb8UX4yyc8kuTXJC5O8KclnkzyS5IYl+3B92g/ule3WCCDFKG1yLICYGC9jDBsqAIyZAEJXH0ty/4Ln7k+7I9aqvjXth/GuJet3smAnLQFkszZ5MNuQNSX1Zaw1tSh9QwWAMRNA6GrdKyCLPJDk4SVrXQHZsjGv71impqS+jLWmFqVvqAAwZgIIXa27BmSRdyf5XMfXJtaAbMw2D2Ybqqakvoy1piYlB26AMRNA6OoNWW8XrEV+N8kfrNEvAWRDSr8lRQAZvmbZW/PGpsQNFQDGTAChq6/L0eeA3Dj7+tokz0l7u9SBL03ypDnv+d1pP4xvW6NfAsia5k0ix7AoVwAZrmYKp3iXMM4AtRBAWMfbcvkk9Bfn8knoVwaIE7l6YflNSS6mPTH9lUlenuSfJvl8ko8n+Zo1+iSAHOGo31AfNYl83hvev1IAGWJbUgFkmJout+aN0dDjDFATAYR1XJPkziTn067bOD/7+srTzU/k6gDyV5O8M8lHkzya5LG0u2adSfJVa/ZJAJnjuN9Qf+ijnzx2ErlKG+JgNgFkmJout+aN0dDjDFATAYTaCCCHLPMb6mf++K/1Fj5cAZlOTddb88aotJ8FwJgJINRGADlk2d9Q99msAZlGTdcTw8eotJ8FwJgJINRmsgGkj8XjfTS7YE2nZufs3kqfjYNb88aotJ8FwJgJINRmcgGkz8Xj6zbngEyrxhWQ4WoAxkwAoTaTCiDLrO/ou938xvmhZuiD2QQQa0A2qbSfBcCYCSDUZlIBZKj1HSUezCaA2AVrk4YeZ4CaCCDUZjIBxPqO/mtK6stYai44B2SQGoAxE0CozWQCyKr33x/XjtuKd+zrO5apKakvY6pxEvr2ao46TBRgLAQQajOZALLqDkTHhYsPffSTK00iS5mQ9VlTUl/GWLPsrXljNPQ4TyHkAdMhgFCbyQSQVa+ALLt4vNb1HcvUlNSXmmvGaMgxXGaziRpucwOmQwChNpMJIF12IJpyuFimpqS+1FwzRkOO4VQW+gPTIYBQm8kEkKbpNjEpbTJaUk1Jfam5piTLrqkYagyntNUxMB0CCLWZVADpsgNRaZPRkmpK6kvNNSVYdU3FUGM4pcMegekQQKjNpAJI04xnIjWGmpL6UnPN0LqsqRhqDFfdbGLn7F63QQHYIgGE2kwugBywvkMAGUvN0MZ066IrIECNBBBqM9kA0jTlTTTHVlNSX2quGVLXNRVDjaE1IECNBBBqI4AUNNEcW01Jfam5ZkhdrygMOYZ2wQJqI4BQmyoDSOk79dRSU1Jfaq4ZUtc1FUOOYZfNJgBKJoBQm6oCiAXm260pqS811wxpjFdAmsZJ6EBdBBBqU00AGdNOPbXUlNSXmmuGNLY1IIctu9kEQMkEEGpTTQAZ0049tdSU1Jeaa4Y29n+3xjDGAEcRQKhNFQFk7L+lHWtNSX2puWZoYz/AcwxjDHAUAYTaVBFAxnqf+thrSupLzTUlGPP6qrGMMcAiAgi1qSKAjHGnnhpqSupLrTXL7ui2LWM8wFMAAcZOAKE2VQQQV0CGqSmpL7XVlLyLUwnjs0qNAAKMnQBCbaoIINaADFNTUl9qqumyo9s2DT0+q9YIIMDYCSDUpooA0jTj36lnjDUl9aWmmtJP8h56fFatEUCAsRNAqE01AWTsO/WMsaakvtRS0/Vq3jaVPoZd3gOgZAIItakmgDTNuHfqGWNNSX2ppabreqZtKn0Mu7wHQMkEEGpTVQA5MMadesZYU1JfaqnpuqPbNpU+hl3eA6BkAgi1qTKANE1ZE6Baa0rqSy01roD0XyOAAGMngFAbAURN55qS+lJLjTUg/dcIIMDYCSDURgBR07mmpL7UVGMXrH5rBBBg7AQQaiOAqOlcU1JfaqrpsqNbH5Y9dX3o8Vm1RgABxk4AoTYCiJrONSX1pbaabZ6EXvvucQIIMHYCCLURQNR0rimpL7XWLLujW1ddTl0vaXyWqRFAgLETQKiNAKKmc01JfZl6TVdd1puU9r0LIEDtBBBqI4Co6VxTUl+mXtNF1x23SvveBRCgdgIItRFA1HSuKakvU6/pouuZI6V97wIIUDsBhNoIIGo615TUl6nXdNH11PXSvncBBKidAEJtBBA1nWtK6svUa7pwBQRgHAQQaiOAqOlcU1Jfplyz7Pkd8163SgCxBgRgGAIItRFA1HSuKakvU6zp46wQu2ABlE8AoTYCiJrONSX1ZWo1Xc7vmKfLqetDf++r1gggwNgJINRmdAFk2dtNSpoA1VpTUl+mVtPlysUiTkIHKJsAQm1GE0BqnySNsaakvkyppsvajWWC+7Knrpc+Pl3eA6BkAgi1GUUA6XK7SUkToFprSurLlGpW3b3q5je+f9LBXQABxk4AoTajCCBTWCg7xpqS+jKlmlXP75h6cD/u+a67iAFsiwBCbYoPIFPZKnSMNSX1ZUo1q14BmXpwX/R8H7uIAWyDAEJtig8gUzksbYw1JfVlSjWrhvKpB/d5z/e1ixjANggg1Kb4ALLq7SY7Z/eapilrAlRrTUl9mVrNsrclCu7rjd8yu4gBbJoAQm2KDyCugJRbU1JfplazzPkdgvv857ve1gkwFAGE2hQfQKwBKbempL5MseaoNQzPe8P8na+mGNwPP9/1lxoAQxFAqE3xAaRp7IJVak1JfZlyzbzzOwT3xc93va0TYCgCCLUZRQBZ5naTKWwnWlpNSX1R0886h6H7vIkaV0CAsRNAxuOJSZ6Z5G/M2jNnj/F4owggTeMk9BJrSuqLmm7rRKYQ3A8/bw0IMDYCSNm+LMnLk/xGkkeTfP5Qe3T23MuSfPlAfSzNaALIgXm3m8xT0gSo1pqS+qJm9XUiUwnu8563CxYwJgJImf5ykp9J8qdpg8ZHkvxCkjuTvDRt4LgzyT9Jsjer+dMkPz177ZSNLoA0TVmTmynXlNQXNauvEym9z33VzHu+y9UhgKEIIGX6VJLfS/I/JnnqEvVPTfJjs9d8aoP9GgMBRE3nmpL6okbNqp9TJ6EDYyGAlOl7kjyhw+uekOS2nvsyNgKIms41JfVFjZqun9Nlrw4BDEUAoTYCiJrONSX1RY2aTX9OAYYigFAbAURN55qS+qJGjQAC1EoAGY8nJHl+kh9O8g+TvPZQ+4fDda0oAoiazjUl9UWNGgEEqJUAMg43JvlY2t2uvrCgfX6w3pVFAFHTuaakvqhRI4AAtRJAxuH/TvJvk/x3Sf7jgftypScmeXWS30/y2OzPV2f5AxK/Nsn7knx61t43e2wdAoiazjUl9UWNGgEEqJUAMg5/nuSOoTsxx8+l/fC8PclLZn82Sd6yxGufleRPkpxP8qpZu5B2G+FnrdEnAURN55qS+qJGjQAC1EoAGYd/neRHhu7EIc9Ne+vXmw89/ubZ48895vW/kvYk96dd8djTZo/98hr9EkDUdK4pqS9q1AggQK0EkHF4TZJ/keVvbdqG16f94Dzj0OPPmD3++iNe++S0t2y9Y85z70jy2VlNFwKIms41JfVFjRoBBKiVADIOtyT5l0l+J+2tTt+a5FvmtG26L8kfLXjuE0nuPeK1z0v7ofuhOc+9bPbcN3bslwCipnNNSX1Ro0YAAWolgIzDvB2vrmxD7IL1kSQPLHjuw7PnF3lh2g/dC+Y89x2z5164RB+uT/vBvbLdmi0GkPd+1/c3577pluY3v/n5zfn//jvntt/85ucfW3Pum25p3vsN/3Vz7ptuUTNgTUl9UaNm05/TZf7bpEZN15qS+qJm+Zr3ftf3b2X+JICMw/ct2bbpY0nuX/Dc/Wl3xFrke9J+6G6Z89wts+duW6IPO7Paq9q2Asi5b7qleejZz9E0TdM0TRt9O/dNt2xl/iSA0JUrIE3T3P3t3+s3kRXVlNQXNWo2/Tkt6aqOmvpqSuqLmuVr7v72793K/EkAoStrQBr3YgPjVdK6FjX11ZTUFzXlzVUEkDL9vXTb8eqJSb6/574s8oZsbhesxzKSXbBK/JcaYBmlTW7U1FVTUl/UlDdXEUDK9G/TrrH40SRPWaL+qUl+bPaaRzbYryt9XY4+B+TG2dfXJnlO2tulrvSutGd+PPWKxw7OAXnXGv0SQACWUNrkRk1dNSX1RU15cxUBpExfnuQfpZ2M//skv5vknyT5B0lePmv/IO3J4x+Z1XwmyU/OXrstb8vlk9BfnMsnob/tipoTs8fuOvTaZye5lPYk9FfO2vnZY89eo08CCMASSpvcqKmrpqS+qClvriKAlO0/Srsm4lzagHF4O95Pz577oVnttl2T5M60weFzsz/vnD1+4ETmB5CkvUpyb9rv7TNJfj2Xr5x0JYAALKG0yY2aumpK6oua8uYqAsh4PDHtZP7krJ1I8oQB+1MqAQRgCaVNbtTUVVNSX9SUN1cRQMbhpiVq/u7GezEOAgjAEkqb3Kipq6akvqgpb64igIzDY2kXpM/zZUneme2fhF4qAQRgCaVNbtTUVVNSX9SUN1cRQMbhl9Ku+fiNPH5XrP8qyYW0AeXHB+hXiQQQgCWUNrlRU1dNSX1RU95cRQAZj+9Lu0PUH6c97fun0u5+9a+S/PUB+1WaogLIgxcvNSdP7TY33HFPc/LUbvPQw5e20i+A45Q2uVFTV01JfVEjgLCeE0n+n7S3W30+yc8nuW7IDhWoiABy4ZFHmxe99f7mhjvuuaq96K33NxceeXQr/QNYpLTJjZq6akrqixoBhPW8LMmfpt1+9wtJfivJMwftUXkGDyDnH3m0uel1980NHwftptfdJ4QAgyptcqOmrpqS+tJXzbJ3NZTU5z5r+iSAjMNXJ/nVtKHj7Ozr70p7Yvpn0h4CSGvwALLoyse8KyEAQ5j6REqNALJKzap3NZTQ503U9EkAGYdPpL3y8fJDjz8lyW7a27Heve1OFWrQAPLgxUtLhY+DZk0IsE0mUmoEkNVqutzVMHSfN1XTJwFkHD6c5DlHPP/3k/z5lvpSukEDyOlz+ysFkDO7+1vpJ4CJlJpt1pTUl3VqutzVMHSfN1XTJwFkHK5douZrN96LcRg0gOyc3VspgOyc3dtKPwFMpNRss6akvixTM++2xK53NZT0ffVZ0ycBhNq4AgJwiImUmm3XlNSXo2qOui3xeW94f6f/p5fwfW2ipk8CCLWxBgTgkK6/HCltAqRmPDUl9WVRzTK3Ja7SDu5qGPr72lRNnwQQamMXLIBDut4eWtoESM14akrqy6KaZf9/Lbj3TwChNoMHkAvOAQEK4wqImm3XlNSXeTWr3rGwTKv91sU+CSDUZvAA0jROQgfKYg2Imm3XlNSXeTWrhvLj2hQ2b+iTAEJtigggBx56+FJzZne/2Tm715zZ3bfmA2eDqGIAAB4cSURBVBiMXbDUbLOmpL7Mq1n1tsSj2lS2r+6TAEJtigogAKXocntoaRMgNeOpKakv82pWvQJy8xvn74g1pQM8+ySAUBsBBGABJ6GrEUC6rQF56OFLS93VMO88kdK+9641fRJAqI0AAnCMZW8PLW0CpGY8NSX1ZVFNn7tWTiHc90kAoTYCCEBPSpsAqRlPTUl9WVTT166Vy5wnUsPtjX0SQKiNAALQk9ImQGrGU1NSX46q6WPXyqls8NAnAYTaCCAAPSltAqRmPDUl9WWZmq67Vk5pi+s+CSDURgAB6ElpEyA146kpqS/L1nQxpUM++ySAUBsBBKAnpU2A1IynpqS+LFvTxarnieyc3dv69yWAwOYJIAA9KW0CpGY8NSX1ZdmaLlwB6UYAoTYCCEBPSpsAqRlPTUl9WbamC2tAuhFAqI0AAtCT0iZAasZTU1Jflq3pyi5YqxNAqI0AAtCT0iZAasZTU1Jflj2dvKsu54mUND7L1vRJAKE2AghAT0qbAKkZT00JfenjjI9lOQl9NQIItRFAAHpS2gRIzXhqhu5Ll9PJ+7DseSJDj0+Xmj4JINRGAAHoSWkTIDXjqRm6L13WZWzT0OPTpaZPAgi1EUAAelLaBEjNeGqG7EvXnam2qaSf1bI1fRJAqI0AAtCT0iZAasZTM2Rfup7NsU0l/ayWremTAEJtBBCAHiy7c1BpkyQ1ZdQM2Zeup5NvU0k/q2Vr+iSAUBsBBGANU9jNR03dAcQVkM3U9EkAoTYCCEBHXXYOKm2SpKaMmiH7Yg3IZmr6JIBQGwEEoKOpnOispu4A0jR2wdpETZ8EEGojgAB00PW3xqVNktSUUTN0X7qcTr5NQ49Pl5o+CSDURgAB6KDrffOlTZLUlFFTQl+2eRL6qkoYn1Vr+iSAUBsBBKCDrjsHlTZJUlNGTUl9WfZ08m0Z6w5zfRJAqI0AAtCBKyBq+qwpqS+lGPsOc30SQKiNAALQgTUgavqsKakvJahhh7k+CSDURgAB6MguWGr6qimpLyWo4d+tPgkg1GZrAWTZezgBxqLLzkGlTZLUlFFTUl+GVsvVxT4JINRm4wGk5F01ANY19vvU1ZRRU1JfhlbL+qo+CSDUZqMBpMs9nABjtOzOQaVNktSUUVNSX4ZWyw5zfRJAqM1GA0jpJ6sCbFtpkyQ1ZdSU1JehuQJyNQGE2mwsgHS9hxOgZqVNktSUUVNSX4ZmDcjVBBBqs7EA0vU3GAC1GuuBamoEkG2zC9bjCSDUZmMBpOs9nAC1sVBdjQCymhp2mOuTAEJtXAEB2KAaDlRTI4AMYezBvU8CCLWxBgRgg2q4lUSNADKkse4w1ycBhNrYBQtgQ7r8IsY6kWnWlNSXsap5DAUQarPRANLlHk6AWqx6K+rNb3z/UrebCCn11ZTUl7GqeQwFEGrjJHSADVl1M47jflnzoY9+ctT3xKsZPoAsG17HqKSfZ98EEGqz8QByYNl7OAFqseoVkOPaM3/811a6olzahEzNcAFkCr8MLOnn2TcBhNpsLYAATM2qa0D6aBazj7Nmk39Pl53Yxqikn2ffBBBqI4AAbNCym3H02SxmH1/NJv+eqWwIU9LPs28CCLURQAA2aJnNOPpuyy5mb5ryJm1TrdnU3zOlLfFL+nn2TQChNgIIwIYddf/9894wPyxsqlknUmbNpv6eKR0KXNLPs28CCLURQAC2ZN5mHNaJqNnk37PqTmw7Z/cWvnfpSvp59k0AoTYCCMDAhlon0jTlTdqmWuMKyPpK+nn2TQChNgIIwMCWWSdy3Ba8q7aDiWZpk7ap1mzq77EGZLiaPgkgrOtrk7wvyadn7X2zx5ZxIu2Hb177hY79EUAACnDcOQ0f+ugne13MfnCrTWmTtqnWbPLvsQvWMDV9EkBYx7OS/EmS80leNWsXknxq9txxTqT98L0nyW2H2jd27JMAAlCQow5t7XMxuysgZdVs8u9Z5gqbc0D6r+mTAMI6fiXJo0medsVjT5s99stLvP5E2g/fqR77JIAAjEwfi9mtASmrZtN/j5PQt1/TJwGErp6c5LEk75jz3DuSfHZWc5QTuRxArpu1dQkgAJXocqtNaZO2qdZs6+856grbmJV48GafBBC6el7aD84PzXnuZbPnjruN6sSs7jO5vPbj95L88Br9EkAAKtHlVpuSJuFTrhnrxHhoq17ZGes4CyB09cK0H5wXzHnuO2bPvfCY93h6kvcnecXsfV6a5Ldnr33TEn24Pu0H98p2awQQgGqUPCFTI4D06XzhgbtPAghJcm2S5yzZnj57zfek/eDcMuf9bpk9d1uHvjwxyW8m+XySZx5Tu5MFu2gJIAB1WfZWm5Im4VOuGevEeEil33LYJwGE5OjtcA+3D8xe08cVkEVeMHv9Dx5T5woIAI9T0iR8yjVjnRgPZQybLvRJACFpF38/f8n2DbPX9LEGZJEbZ6+/s8NrrQEBmKgSF+5OtWasE+OhdD3hfazjLIDQ1XG7YD2W43fBWuQ7034oX9zhtQIIwMRYJ1JezVgnxkPZObu3UgAZ4uDNPgkgrONdac/8eOoVjx2cA/KuQ7XXp11Dcu0Vj33FnPe8LsmHk3wujz9fZFkCCMCElL5wd6o1Y50YD8UVEFjes5NcSnsS+itn7fzssWcfqr0r7QftxBWPvTvJvWkXk78kyWuTfGxW9+qOfRJAACak9IW7U60Z68R4KEOuAVn21sU+CSCs68a0IeIzs/brs8cOuytXB5AXJ/lgkk8k+Yskn0ryG2kXsXclgABMxBgW7k61RgBZ3bbD9JCnyQsg1EYAAZiIMdy2MtUaAWR12zx4s8uti30SQKiNAAIwEWNYuDvVGgGkm21tqNDlakufBBBqI4AATIQrIOXWCCDr2eTBm11vXeyTAEJtBBCAibAGpNwaAWQ7uoxh1+DeJwGE2gggABNiF6wyawSQ7egyhl1vXeyTAEJtBBCACdnmwl01y9cIINvhCgiUQQABmBgnoZdX08d7DHE+xdhYAwJlEEAAJmqTC3fVrFazznsMeT7F2HQdZ7tgQb8EEACOVNJEvdaaru8x9PkUY9N1nLvcutgnAYTaCCAAHKmkiXqtNWP9zfyYLHuLWolXmgQQaiOAAHCkkibqtdaMdW3CGPS95mnZWxf7JIBQGwEEgCOVNFGvtabLe5SwO1PputyitszPYtsEEGojgABwpJIm6rXWdHmPEs6nKN2mzr3ZNgGE2gggAByppIl6rTWugPSvyy1qpW5lLIBQGwEEgCOVNFGvtcYakP6tGtBufuP7F14dGXoXMQGE2gggAByppIl6rTVd38MuWIuteovaUW3orYwFEGojgABwpJIm6rXWdH2Poc+nKNmqV0BKDnECCLURQAA4UkkT9Vpr1nkPJ6HPt+otasu0oW5jE0CojQACwJFKmqjXWtPHewxxPkXplr1Fbdk21EJ+AYTaCCAAHKmkiXqtNX39PTzeMreordKG2spYAKE2AggACy27LWlJk/kx1gggm3PULWrPe8P8na9cAYHNEkAAuMqq6wpKmsyPsUYA2bx5t6iNZStjAYTaCCAAPM75DjsrlTSZH2ONADKcMWxlLIBQGwEEgMfpMiEraTI/xhoBZDhj2MpYAKE2AggAX9T1lpSSJvNjrBFAhlX6VsYCCLURQAD4olUPbztYlFvSZH6MNQJIGUrdylgAoTYCCABftHN2b6UAcrAtaUmT+THWCCAcRQChNgIIAF/kCsgwNQIIRxFAqI0AAsAXWQMyTI0AwlEEEGojgADwOHbB2n6NAMJRBBBqI4AA8DhdtiUtaTI/xhoBhKMIINRGAAHgKk5C327Ncc8/ePFSc/LUbnPDHfc0J0/tFrM7E9shgFAbAQSAhZbdlrSkyfwYaxY9X/r5FGyHAEJtBBAA1lbSZH6MNfOePz+CE7rZDgGE2gggAKytpMn8GGvmPd9lMwDqJIBQGwEEgLWVNJkfY83h57tuh0ydBBBqI4AAsLaSJvNjrDn8fNcDIamTAEJtBBAA1lbSZH6MNYef3zm7t1IA2Tm7t/DvZvwEEGojgACwtpIm82OscQWEowgg1EYAAWBtJU3mx1hjDQhHEUCojQACwNpKmsyPscYuWBxFAKE2AggAaytpMj/GmnnPX3AOCDMCCLURQABYW0mT+THWLHreSeg0jQBCfQQQANZW0mR+jDXHPf/Qw5eaM7v7zc7ZvebM7r41HxMjgFAbAQSAtZU0mR9jzTLvwXQJINRGAAFgbSVN5sdYI4BwFAGE2gggAKytpMn8GGsEEI4igFAbAQSAtZU0mR9jjQDCUQQQaiOAALC2kibzY6wRQDiKAEJtBBAA1lbSZH6MNQIIRxFAqI0AAsDaSprMj7FGAOEoAgi1EUAAWFtJk/kx1gggHEUAoTYCCABrefDipebkqd3mhjvuaU6e2l14SF5JE/6SapYdP6ZLAKE2AggAnVx45NHmRW+9v7nhjnuuai966/3NhUcefVx9KRP+UmpWHT+mSwChNgIIACs7/8ijzU2vu2/u5Pmg3fS6+x43iR56wl9STZfxY7oEEGojgACwskW/uZ/3m/wDpYeCbdZ0GT+mSwChNgIIACt58OKlpSbPB+1gTUPpoWBbNV3Hj+kSQKiNAALASk6f219pAn1md79pmrJDwTZruo4f0yWAUBsBBICV7JzdW2kCvXN2r2maskPBNmu6jh/TJYBQGwEEgJW4ArJejSsgrEoAoTYCCAArsQZkvRprQFiVAEJtBBAAVmYXrPVq7ILFKgQQaiOAALCyC84BWaumy/gxXQII67glyT9O8kCSz6X9IJ1Y8T2emOTVSX4/yWOzP189e7wLAQSATpyEvl6Nk9BZlgDCOu5KGxp+J8lH0i2A/NzsdW9P8pLZn02St3TskwACwFoeevhSc2Z3v9k5u9ec2d1fuGZh6Al/qTXLjh/TJYCwjqckedLsn09l9QDy3CRfSPLmQ4+/efb4czv0SQABYCtKmfCXWgOLCCD0pUsAef3sNc849PgzZo+/vkM/BBAAtqK0CX9pNbCIAEJfugSQ+5L80YLnPpHk3g79EEAA2IrSJvyl1cAiAgh96RJAPpJ2Afs8H549f5Tr035wr2y3RgABYAtKm/Bvq+bBi5eak6d2mxvuuKc5eWrXGg9WJoCQJNcmec6S7ekL3qNLAPlYkvsXPHd/2h2xjrIz+zuvagIIAJtWUijYRo1druiLAELShoa5E/k57QML3sMVEAAmpYRQsK2a8875oEcCCElyXZLnL9m+YcF7WAMCwKQMHQq2WeOkc/okgNCXLgHkDbELFgAjNXQo2FbNgxcvLRU+Dpo1IRxHAKEvxwWQ69OuIbn2ise+LkefA3Jjh34IIABsRenBoa+a0+f2VwogZ3b3F743NI0AwnpuTPKaWftA2g/ST8++fsWh2rsyP6C8LZdPQn9xLp+E/raOfRJAANiK0oNDXzU7Z/dWCiA7Z/cWvjc0jQDCem7P4sXqHz9Ue1fmB5BrktyZ5HySz83+vHP2eBcCCABbUXpw6KvGFRD6JoBQGwEEgK0oPTj0VWMNCH0TQKiNAALAVpQeHPqssQsWfRJAqI0AAsBWDB0KtllzwTkg9EgAoTYCCABbMXQo2HaNk9DpiwBCbQQQALaihFAwRM1DD19qzuzuNztn95ozu/vWfLAyAYTaCCAAbEVJoWDbNbAOAYTaCCAAbEVpoUAAYSwEEGojgACwcQ9evNScPLXb3HDHPc3JU7sLb0MqLTgcV7Ps9wXrEECojQACwMasuhC7pHBxVI0F5myTAEJtBBAANuJ8h61oSwgXx9V0+b5gHQIItRFAANiILofxDR0ulqlxyCDbJoBQGwEEgN49ePHSUpP0g3awdmLocHFcTdfvC9YhgFAbAQSA3p0+t7/SRP3M7n7TNOUHkK7fF6xDAKE2AggAvds5u7fSRH3n7F7TNOUHkK7fF6xDAKE2AggAvXMFxBUQ+iOAUBsBBIDeWQNiDQj9EUCojQACwEbYBcsuWPRDAKE2AggAG3Gh0nNAunxfsA4BhNoIIABsjJPQYX0CCLURQADYuIcevtSc2d1vds7uNWd29xeujSghXKxSs+z3BesQQKiNAAJAMUoKF8vWwKYJINRGAAGgGKWFCwGEEggg1EYAAaAYJYWLBy9eak6e2m1uuOOe5uSpXbdXMRgBhNoIIAAUo4QAYoE5pRFAqI0AAkAxhg4g522xS4EEEGojgABQjKEDiEMGKZEAQm0EEACKMWQAefDipaXCx0GzJoRtEUCojQACQDGGDCCnz+2vFEDO7O53/0ZhBQIItRFAACjGkAFk5+zeSgFk5+xe928UViCAUBsBBIBiuAICVxNAqI0AAkAxrAGBqwkg1EYAAaAIyx78ZxcspkYAoTYCCACDWvXgv00GkAvOAaFAAgi1EUAAGEyXg/+chM7UCCDURgABYDBdbnna1jqRhx6+1JzZ3W92zu41Z3b3rflgMAIItRFAABhE10XffYSLZdebQAkEEGojgAAwiK7b3rq9iqkRQKiNAALAILoe/Nc1gHRZbwIlEECojQACwCC2fQXEFruMlQBCbQQQAAaxzTUgDhlkzAQQaiOAADCYbe2C1fVqC5RAAKE2AggAg+ly8F+XANJ1vQmUQAChNgIIAIPaxknoroAwZgIItRFAACjCsgf/WQPC1Agg1EYAAWA0ljlAcFGNXbAYKwGE2gggABRvmdu0jqv50Ec/6RwQRkkAoTYCCABFW+YAwef+xL3NjTvHh4sPffSTTkJndAQQaiOAAFC0ZW+dWuX2qmXXm0AJBBBqI4AAUKxVF49bYE6NBBBqI4AAUKxVt8+1xS41EkCojQACQLFWPUBwmeaQQcZGAKE2AggAxXIFBAQQ6iOAAFAsa0BAAKE+AggARdvELlgwJgIItRFAACjahR7PAXHOB2MkgFAbAQSA4vVxErrwwVgJINRGAAFgNJY5QNAhg9RGAKE2AggAQMEEEGojgAAAFEwAoTYCCABAwQQQaiOAAAAUTAChNgIIAEDBBBBqI4AAABRMAKE2AggAQMEEENZxS5J/nOSBJJ9L+0E6scLrT8xeM6/9Qsc+CSAAAAUTQFjHXUkeS/I7ST6S7gHkPUluO9S+sWOfBBAAgIIJIKzjKUmeNPvnU+keQE712CcBBACgYAIIfVk3gFw3a+sSQAAACiaA0Jd1Ashncnntx+8l+eE1+iGAAAAUTAChL10CyNOTvD/JK5K8IMlLk/z27H3etMTrr0/7wb2y3RoBBACgWAIISXJtkucs2Z6+4D26BJB5npjkN5N8Pskzj6ndyYJdtO6+++7mgQce0DRN0zRN0wprd99998Gc7eY1542M2Iks3g73cPvAgvfoK4Ak7dWQJskPHlM37wrInVn+e9E0TdM0TdOGa7eGybouyfOXbN+w4D36DCA3zt7rzg6v/cq0H+abc3U42US7NZf/BdrG3zfVZpyNcS3NOBvjWppxNsbrtJtn39dXBtbQZwD5ztl7vbiH99q0r0/b168fuiOVM86bZ4y3wzhvnjHeDuO8ecYYjnFcALk+7RqSa6947Cvm1F2X5MNpT1Z/Wo/92xT/cdgO47x5xng7jPPmGePtMM6bZ4xhjhuTvGbWPpD2X5Kfnn39ikO1d+XqgPLuJPemXUz+kiSvTfKxWd2rN9XpnvmPw3YY580zxtthnDfPGG+Hcd48Ywxz3J7Fi4o+fqj2rlwdQF6c5INJPpHkL5J8KslvJPmOjfW4f/7jsB3GefOM8XYY580zxtthnDfPGANzXZ/2Cs71A/ejdsZ584zxdhjnzTPG22GcN88YAwAAAAAAAAAAAAAAAAAAAAAAADCwJ6Y9LPH3kzw2+/PVs8dZ3ZPTbjX4q0n+MO2+53ctqDX2q/vPk/yvSX43yWeS/FHaM3f+9pxa49vdX0vyfyT5vSSPJvl0kv83ySuTPOlQrXHuz7fk8hlU/9mh565L8lNJ/iDJZ5PsJfmBrfZunE5k8Tlfv3Co1md5PdcneWuSf512/B5O8p4kX3Oo7gfSfn4/m/bz/FNpP9/AhPxc2v8Qvz3tKe5vn339liE7NWIn0o7fw2lDyFEBxNiv7leSPJLkbUl+MMnfT/KRtOP2skO1xre7W5KcS/K/JHlpkpcneWeSL6T9XF/JOPfjLyX5/9IGvnkB5NfSHnR7Ju0E7j2zuh/bYh/H6ETacXpPktsOtW88VOuz3N2z0v7S7XyS1yb5/iR3pB33Z11R9z/l8s/jB9J+nv99knu22VlgWM9NO6F486HH3zx7/Llb79H4PSnJU2b/fE0WBxBj381/kat/A39dkv0k/y7tmCfGd1P+t7Sf6WfPvjbO/fnxJJ9IOyE7HEC+ffbYqw695r1J/izJV2+jgyN1Iu3YnTqmzme5u/8gyb9M8kDauwAW+eq0n9f3Hnr8VWl/Rt+2kd4BxXl92n/pn3Ho8WfMHn/91ntUl6MCiLHv18+kHbenzb42vpvxo2nH72/MvjbO/bghyZ8m+Xtpb+E8HEDemeTPc/VtKn9rVutWrMVO5HIAuS6Lb/XxWe7u4NbBb599/R+mvaJ32A/M6v7Woce/JO3n+52b6iBQlvvS3kM/zyeS3LvFvtToqABi7Pv1v6e9PeVLZ18b3358SZKvSjtBflGSTya5mMuTOOPcj/cm+edpf5O8k6sDyH6S35rzuutmtW/bcP/G7ETaMfpMLq/9+L0kP3yozme5u3+Udlz/ZpJ/NvvnL6T9TD/virqfnz03LwT+iyT/arPdBErxkbSXTOf58Ox5ujsqgBj7/vy1tAse33XFY8a3Hzt5/KLd30py4xXPG+f1fXuSzyf567Ovd3J1APlMHv/5vtK/y9Xrcrjs6Unen+QVSV6Qdk3Tb6cd4zddUeez3N3BeqRHkrw77S8rfiTJH6e9snFw+9qvzh6b591pP+fABHwsyf0Lnrs/7Q4gdHdUADH2/fjyJA8l+VTaicYB49uP/zTtDmPfnfa3l/88yTdf8bxxXs91aRft/vwVj+3k6gDy+SS/uOA9Hk47wWZ5T0zym2nH9Zmzx3yWu3t/2s/s4c/hfzl7/JdmX/9G2s/rPL+YdjE6MAF+47NZroBs1nVJPpj2N2x/89BzxnczXpXkc2mvOiXGeV2n0l7B+KorHtuJKyDb8IK04/yDs699lrs72PHx++Y89/G0t7Ad1LkCArjndcOsAdmcv5R2jD6Xywsfr2R8N+Nr8vgFuca5u/8k7TkIP5k2bBy0n007xt+S9gpUYg3IJtyYduzunH3ts9zd29KO5fPnPPdbaf87nVgDAsy8IXb92KSjAoix7+6atPccfz7J311QY3w340Ta8fu52dfGububsviAvIP26Kz2F2MXrL59Z9qxe/Hsa5/l7l6cdoxeMue5fzNrSXu1ad4uWNfFLlgwKV+Xo/c9v/GqV7CKowKIse/mCWlP6D5u0mV81/NXFjx+sNvNwa0Wxrm7L0/yXXPaL6Ud41ck+Tuz2oPbhRadA7Lo50XyFXMeuy7tbVWfy+Wtu32Wu/uqtNtI/1Yef2r8wfk1ByfO/5W0QeM9h15/cA7ICzbbTaAkB5dO3572txgHJ7+6pN/dK5K8Ju1psE3a/9G9Ztau/J+YsV/d6bRj9IFcfarxbWlvETpgfLv7v9Kur/mf0wa9H8vlhaYfyuUDHxPj3LedzD8J/d60W02fTvub5oOdh169zc6N0LvTjt1O2nF7bdoF5/PGzme5ux9JO1YfTPv/wDemDSWfzOWQl7Rj3qT9b8xL0n6e/yLJr2+zs8Dwrkl7D+z5tL8NOj/7+pqjXsSRPp7Ft1TcfkWdsV/dB3L0LSvffEWt8e3uu9NOCC6mHbvPpN269Edz9Un0xrlfO5kfQL4k7bax/ybtttMPpd1SlqO9OO2k+BNpJ7qfSrsb03fMqfVZXs9taX/h9tm0myP8n7m8julKL037+X0s7ef5TWk/3wAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAALBptydpkvw3SV6T5A+SfDbJP0vydcN1CwAAqNHtaQPIb8/aq5L8eJI/TnIhyTWD9QwAAKjO7WkDyIeTXHvF439n9vh/O0CfAACASt2eNmi89NDjf3n2+P+w7Q4BAAD1uj2X14Ac1iT5ia32BgAAqNrtaYPG357zXJNkZ5udAQAA6nZ7BBAAAGBLbo8AAgAAbMntEUAAAIAtuT0CCAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAsLr/H+M3AQO41jGjAAAAAElFTkSuQmCC\" width=\"640\">"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/plain": [
       "Text(0,0.5,'x(n)')"
      ]
     },
     "execution_count": 12,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "plt.figure()\n",
    "plt.stem(n, x)\n",
    "plt.xlabel('n')\n",
    "plt.ylabel('x(n)')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.6.4"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}