Spike_IO-BEE-Sender-Receiver.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "from spynnaker_external_devices_plugin.pyNN.connections\\\n",
    "    .spynnaker_live_spikes_connection import SpynnakerLiveSpikesConnection"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "# plotter in python\n",
    "%matplotlib notebook\n",
    "import matplotlib.pyplot as plt\n",
    "\n",
    "import time as my_time\n",
    "import numpy"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "# Reads the arrays generated on the simulator side to convert from SpiNNaker to BEE\n",
    "from_spinnaker_exc2bee = numpy.load('from_spinnaker_exc2bee.npy')\n",
    "from_spinnaker_inh2bee = numpy.load('from_spinnaker_inh2bee.npy')\n",
    "\n",
    "# Reads the arrays generated on the simulator side to convert from BEE to SpiNNaker\n",
    "from_bee2spinnaker = numpy.load('from_bee2spinnaker.npy')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "n_neurons = from_bee2spinnaker.shape[1]\n",
    "\n",
    "number_of_liquids = from_bee2spinnaker.shape[0]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "# SpynnakerLiveSpikesConnection(self, receive_labels=None, send_labels=None, local_host=None, local_port=19999)\n",
    "\n",
    "# Set up the live connection for sending and receiving spikes\n",
    "live_spikes_connection = SpynnakerLiveSpikesConnection(\n",
    "    receive_labels=[\"LSM_EXC\",\"LSM_INH\"], \n",
    "    send_labels=[\"spike_injector\"],\n",
    "    local_port=12346)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "# Create an initialisation method\n",
    "def init_pop(label, n_neurons, run_time_ms, machine_timestep_ms):\n",
    "    print \"{} has {} neurons\".format(label, n_neurons)\n",
    "    print \"Simulation will run for {}ms at {}ms timesteps\".format(run_time_ms, machine_timestep_ms)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "received_spikes = {'LSM_EXC':[],'LSM_INH':[]}\n",
    "def receive_spikes(label, time, neuron_ids):\n",
    "    received_spikes[label].append([time,neuron_ids])\n",
    "#     my_time.sleep(feedback_delay)\n",
    "#     live_spikes_connection.send_spikes(\"spike_injector\", neuron_ids)\n",
    "#     for neuron_id in neuron_ids:\n",
    "#         print \"Received spike at time\", time, \"from\", label, \"-\", neuron_id\n",
    "    live_spikes_connection.send_spikes(\"spike_injector\", initial_array[range(0,4)+range(700,704)])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "initial_array = numpy.arange(n_neurons*number_of_liquids)\n",
    "# numpy.random.shuffle(initial_array)\n",
    "def started(label, sender):\n",
    "    print \"{}=>Simulation has started!\".format(label)\n",
    "    sender.send_spikes(\"spike_injector\", initial_array[range(0,4)+range(700,704)])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "# Here is the best place to pause SpiNNaker if it is necessary to wait until\n",
    "# something on the external machine is ready.\n",
    "# After the function below returns, SpiNNaker is going to ack the simulation start\n",
    "# and the function supplied for add_start_callback will be called.\n",
    "def database_callback_function(something):\n",
    "    # my_time.sleep(50)\n",
    "    print \"database callback...\"\n",
    "    # print dir(something)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {
    "collapsed": false,
    "scrolled": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "spike_injector has 3000 neurons\n",
      "Simulation will run for 1000.0ms at 2.0ms timesteps\n",
      "LSM_EXC has 2400 neurons\n",
      "Simulation will run for 1000.0ms at 2.0ms timesteps\n",
      "LSM_INH has 600 neurons\n",
      "Simulation will run for 1000.0ms at 2.0ms timesteps\n",
      "database callback...\n",
      "spike_injector=>Simulation has started!\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "spynnaker database connection for None:12346: <Thread(start callback thread for live_event_connection12346:0.0.0.0, initial daemon)>.start(): starting thread\n",
      "start callback thread for live_event_connection12346:0.0.0.0: <Thread(start callback thread for live_event_connection12346:0.0.0.0, started daemon 4570005504)>.__bootstrap(): thread started\n",
      "start callback thread for live_event_connection12346:0.0.0.0: <Thread(start callback thread for live_event_connection12346:0.0.0.0, started daemon 4570005504)>.__bootstrap(): normal return\n"
     ]
    }
   ],
   "source": [
    "live_spikes_connection.add_receive_callback(\"LSM_EXC\", receive_spikes)\n",
    "live_spikes_connection.add_receive_callback(\"LSM_INH\", receive_spikes)\n",
    "\n",
    "live_spikes_connection.add_init_callback(\"LSM_EXC\", init_pop)\n",
    "live_spikes_connection.add_init_callback(\"LSM_INH\", init_pop)\n",
    "live_spikes_connection.add_init_callback(\"spike_injector\",init_pop)\n",
    "\n",
    "live_spikes_connection.add_database_callback(database_callback_function)\n",
    "\n",
    "live_spikes_connection.add_start_callback(\"spike_injector\",started)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "times = []\n",
    "spikes = []\n",
    "\n",
    "for lsm_i in range(number_of_liquids):\n",
    "    times.append([])\n",
    "    spikes.append([])\n",
    "    \n",
    "# If there are spikes, plot using matplotlib\n",
    "if len(received_spikes['LSM_EXC']) != 0 or len(received_spikes['LSM_INH']) != 0:\n",
    "    if len(received_spikes['LSM_EXC']) != 0:\n",
    "        for i,j in received_spikes['LSM_EXC']:\n",
    "            for k in j:\n",
    "                lsm_i, idx = from_spinnaker_exc2bee[int(k)]\n",
    "                (times[lsm_i]).append(i)\n",
    "                (spikes[lsm_i]).append(idx)\n",
    "    if len(received_spikes['LSM_INH']) != 0:\n",
    "        for i,j in received_spikes['LSM_INH']:\n",
    "            for k in j:\n",
    "                lsm_i, idx = from_spinnaker_inh2bee[int(k)]\n",
    "                (times[lsm_i]).append(i)\n",
    "                (spikes[lsm_i]).append(idx)\n",
    "else:\n",
    "    print \"No spikes received\""
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {
    "collapsed": false,
    "scrolled": false
   },
   "outputs": [
    {
     "data": {
      "application/javascript": [
       "/* Put everything inside the global mpl namespace */\n",
       "window.mpl = {};\n",
       "\n",
       "mpl.get_websocket_type = function() {\n",
       "    if (typeof(WebSocket) !== 'undefined') {\n",
       "        return WebSocket;\n",
       "    } else if (typeof(MozWebSocket) !== 'undefined') {\n",
       "        return MozWebSocket;\n",
       "    } else {\n",
       "        alert('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",
       "            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",
       "        this.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 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);\n",
       "        canvas.attr('height', height);\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'];\n",
       "    var y0 = fig.canvas.height - msg['y0'];\n",
       "    var x1 = msg['x1'];\n",
       "    var y1 = fig.canvas.height - msg['y1'];\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;\n",
       "    var y = canvas_pos.y;\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",
       "    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 + '\">');\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 dataURL = this.canvas.toDataURL();\n",
       "    this.cell_info[1]['text/html'] = '<img src=\"' + dataURL + '\">';\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,iVBORw0KGgoAAAANSUhEUgAAAoAAAAHgCAYAAAA10dzkAAAgAElEQVR4Xu3dCZxkZXnv8d+ADgoCilGuLG4RHATxKhGEqIjRTIgroKJRgRjjwiDKuEAwYRUUlEm8MpgEBaPxiguQGAlO9JqIQKJRY4JEwC0KghARZlAMg8D9PNNvyaGo7q7qp+rtU1W/+nwGpqfrOXXO932q+t/v2ZbgQwEFFFBAAQUUUGCqBJZM1da6sQoooIACCiiggAIYAG0CBRRQQAEFFFBgygQMgFM24G6uAgoooIACCihgALQHFFBAAQUUUECBKRMwAE7ZgLu5CiiggAIKKKCAAdAeUEABBRRQQAEFpkzAADhlA+7mKqCAAgoooIACBkB7QAEFFFBAAQUUmDIBA+CUDbibq4ACCiiggAIKGADtAQUUUEABBRRQYMoEDIBTNuBurgIKKKCAAgooYAC0BxRQQAEFFFBAgSkTMABO2YC7uQoooIACCiiggAHQHlBAAQUUUEABBaZMwAA4ZQPu5iqggAIKKKCAAgZAe0ABBRRQQAEFFJgyAQPglA24m6uAAgoooIACChgA7QEFFFBAAQUUUGDKBAyAUzbgbq4CCiiggAIKKGAAtAcUUEABBRRQQIEpEzAATtmAu7kKKKCAAgoooIAB0B5QQAEFFFBAAQWmTMAAOGUD7uYqoIACCiiggAIGQHtAAQUUUEABBRSYMgED4JQNuJurgAIKKKCAAgoYAO0BBRRQQAEFFFBgygQMgFM24G6uAgoooIACCihgALQHFFBAAQUUUECBKRMwAE7ZgLu5CiiggAIKKKCAAdAeUEABBRRQQAEFpkzAADhlA+7mKtASgTuB44ATyvrE348B7gPE99r+6F7/xV7fY4vfxn2syOOAPwP2BG4DPg28Gbipj1qfooACEyJgAJyQgXQzFBgzgd2Ba4Bry3p3Asx9DYALGsltgO2Ar8xT/TDg34H/BE4CHgS8B/gh8NQFvbJFCigwlgIGwLEcNldagYkTMADWGdI/BQ4GHgHcUl7yacAXgf2Bv6mzGr6KAgostoABcLFHwNdXYHwEdgBOBfYCtgBuAP4FeFmZtdsb+EfgRcDzgBcAsUvy74DDgZ82NrV7F2qvAPg7wCeBs0t9lN+/7Dp+MbAt8CPgA8A7gbvK8jcD3gU8H9gaWFtmvQ4DrhoSd/f692uzH/C7xSg+f/+q7H7dAzgNeALwfWAl8A8DrGtnF/pG89R8B/hn4JVdz/uvMna/P8Br+lQFFBhjAQPgGA+eq65AZYFvAzeWcBX/jwAWYeZVwC+BTgC8Gvg88HEgglGEs38FfmuAAHgQcGYJe1EfjwiT/wQsK8cOfhN4Sjn27XTgreV5Ufdc4I+ACDwPBn4T+FQfu0j7Je0OgP3aRNA6D7gQeDrwJ8BfAPsUp9glHv/2xDJL1wzNc61br2MAPwSEYycU3g/4OXA0cErXwi4Afg2IIOpDAQWmQMAAOAWD7CYqMASBCFH/XWbVPjPL8joBMMLNcxrP+T3gI8CzyixTfGuuGcC3ACcCryuzf51FxaxVhJoITpc0lh+BJk4giWPgfgJcBqwBYjmjejTXfxCbDwJ/2FiprwH/uxx/FzNz8Xh8mbGMXbXh1s+jVwCMmdEw26QsII7/ixnTcP3LroXG60SYjsDuQwEFpkDAADgFg+wmKjAkgZhN+59yBmnMxMXXzUcnAMaMYAS1ziNO7PhFCWknl3+cLQCuBl4NvAToDpp/XXY/d4eUJwFfboTTs8rfY5dq7Eb9tz5OLInPwubu01i/zi7lXnzd69+vTey6PrexwI+W3eWxS73pFWfn/jHQ8ZpvCPs5C9gAOJ+i31dgigQMgFM02G6qAkmBR5ZdsrHbN3YXxrFq7wb+vCy3EwBj9i9mAZuPmD38BLBijgAYISZm8GK5zyihsbmMCHMxi9jrEWEtgmccUxfHAMasYIStXy+XN/kw8PYey+wsK44zjBm3eMSyYjmxvNke3QGwX5tY/y80FhqvG7vGH971QrH8d5TQ3M+w9RMA4/jJn7kLuB9On6PA5AsYACd/jN1CBUYhsCsQJ1XEbN2+ZZfrMGYAY3fo54Ary3Jvbaz8x4Anl2DX67Mrjq/rPmZu+3LCRRzzFjOCcVxgr0cEsAi1nUcE0bg0Sr8BsPm8uWwWMwDGOn637D6PYwObjwjdMavrSSCjeLe4TAVaKGAAbOGguEoKjInA5uUM2zj5IsJVJwB+tpwc0tmMlwMxA9fvMYA7llmy2K0a4TJOXIhHzNDFsWtxjNygZ/PGsXZxgkWcnTyMx3wXgu5lEzN/z17EGcDY7rgMTIS/mLHsXAYmrv93kZeBGUZbuAwFxkfAADg+Y+WaKrCYAhG63lvO7I1gFmfkxmxRXDsuTh6I4+yaZwH/P+Ac4LFlV2acBdzcfTvfZWDiOL+YkYqZqbgcTOy6jLuExOxgfC8CZ1zQeCnwmMZlZ+IYxUvL3S3iZJCoi93JcTzdEUCcLTyMR3P9B7GpOQMYJ5xE2ItjMDuPuGD0N8qFoOPs6geWM4Lj5JA4U9qHAgpMiYABcEoG2s1UICnwkHINwLh9WJxtG0ErAlbcTSIu+RKPCIAxy3VAOQnjhSUoxq3G3ti1e/aOcjxhnO0bjziGLS5/EoGucyu4OH4vrisYl5VZXsJcfP8o4KXAo8rsYOzWjBNGYl2iNoLNbwOPLqHxe2XmME4wGdajuf6D2PSaAXxmueRLc91i+XEMYLj08+j4RUjuPOL4wjgLuPlv8b2dgVXlhJr15eLPcca0t4LrR9rnKDAhApMWAA8sB5nHxVQfUH7z7XVf0d3KxVDjIrZxSYnm4/hyXFOclRe7jeKg9csnZLzdDAVGKdAJgN0hZ5Sv6bIVUEABBRYgMGkBMH7wbAVsWu4O0Ou+onFNrNgd9WMgLozaDIBxLFMc2B7HHcWsQvxWHbtQ4pik5sHoC6C2RIGJF+jsAu7ezTnxG+4GKqCAAuMmMGkBsOPfmYnoFQDjxudxva+4PVRcfqEZAGNXUewa6RwnFMc5xYHjcVumuF6XDwUUmF3AGcDRdEd8Ds31iN3FPhRQQIGBBKYtAEbYi9suxW2WjuwKgLHL92YgjnGKi8p2HnFHgTjWaZR3FRho0HyyAgpMjUBnVnW2DY5rFsaxkHNdsmZqsNxQBRToX2CaAmBcHDbOfjukXAcrdu82ZwDjwPb4EN2pXIOsoxhnMq4DXtODNfzirLrO5RT6l/eZCiigwPwCcTjLfLdni3siOws4v6XPUKBbIC7XFHv55rrrz8SqTVMAjLsVxJmLbyqj2R0AFzIDuC1wzcR2hxumgAIKKKDAZAvE5E9cBmnqHtMUAON6Yls2flOO36zjGMHY7RvXMYvj/3odA3hduX5Yr2MAIzSuvfrqq9lii+atPKeuj1qxwUcffTQnn9zvrVNbscoTuxKORbuG1vFoz3g4Fu0Yi3Xr1rH99nGjoA25IPbyTd1j0gJgnNwRoS6Om4l7kcb0buwaiWtdxbW6mtfDenO58Ol+wPXl+mFxnF+cBRz3Mo0weAzwinIx215nAW8IgGvXrjUAtuCts3LlSlatinN4fCy2gGOx2CNwz9d3PNozHo5FO8YiAuCWW0b2MwC2Y0TyaxG3ioqLn3b250fAjb/vU2511HyF7l3Ane8dB7y2hMevznMdQANgfsyGtgQ/WIdGmV6QY5EmHOoCHI+hcqYW5lik+IZWbACESZsBHFpz9LkgA2CfUDWetmbNGpYvjxtG+FhsAcdisUfgnq/veLRnPByLdoyFAdAAmO1EA2BW0HoFFFBAAQUqCxgADYDZljMAZgWtV0ABBRRQoLKAAdAAmG05A2BW0HoFFFBAAQUqCxgADYDZljMAZgWtV0ABBRRQoLKAAdAAmG05A2BW0HoFFFBAAQUqCxgADYDZljMAZgWtV0ABBRRQoLKAAdAAmG05A2BW0HoFFFBAAQUqCxgADYDZljMAZgWtV0ABBRRQoLKAAdAAmG05A2BW0HoFFFBAAQUqCxgADYDZljMAZgWtV0ABBRRQoLKAAdAAmG05A2BW0HoFFFBAAQUqCxgADYDZljMAZgWtV0ABBRRQoLKAAdAAmG05A2BW0HoFFFBAAQUqCxgADYDZljMAZgWtV0ABBRRQoLKAAdAAmG05A2BW0HoFFFBAAQUqCxgADYDZljMAZgWtV0ABBRRQoLKAAdAAmG05A2BW0HoFFFBAAQUqCxgADYDZljMAZgWtV0ABBRRQoLKAAdAAmG05A2BW0HoFFFBAAQUqCxgADYDZljMAZgWtV0ABBRRQoLKAAdAAmG05A2BW0HoFFFBAAQUqCxgADYDZljMAZgWtV0ABBRRQoLKAAdAAmG05A2BW0HoFFFBAAQUqCxgADYDZljMAZgWtV0ABBRRQoLKAAdAAmG05A2BW0HoFFFBAAQUqCxgADYDZljMAZgWtV0ABBRRQoLKAAdAAmG05A2BW0HoFFFBAAQUqCxgADYDZljMAZgWtV0ABBRRQoLKAAdAAmG05A2BW0HoFFFBAAQUqCxgADYDZljMAZgWtV0ABBRRQoLKAAdAAmG05A2BW0HoFFFBAAQUqCxgADYDZljMAZgWtV0ABBRRQoLKAAdAAmG05A2BW0HoFFFBAAQUqCxgADYDZljMAZgWtV0ABBRRQoLKAAdAAmG05A2BW0HoFFFBAAQUqCxgADYDZljMAZgWtV0ABBRRQoLKAAdAAmG05A2BW0HoFFFBAAQUqCxgADYDZljMAZgWtV0ABBRRQoLKAAdAAmG05A2BW0HoFFFBAAQUqCxgADYDZljMAZgWtV0ABBRRQoLKAAdAAmG25ngHwM5+Bo46aWfT++8N5503W39etgyVLYPPNJ3P7JnXcMtvlmE/We3i2XnCcp2Ocm+M/rWN+xx3ruOKKLYMi/rMuGwbGsX7JOK50i9a5ZwDcZRe4/PKZtdxkE7jttsn7e2cMJnX73K5796xjPrnv52a/O87TMc6OeWQ+A2CL8tTYrYoBcEIDrgHQADiJv7j109cGQAPgJE9e3P0eMAA6A5jLnPPuAj7gADj33JkXmZS/N3cZTMo2TdL4jGJbHPPJeg/P1iOO83SMc3P8p3XM3QU8eccAHgisAJ4APAC4L3BnyXi7A38MPBm4P/AD4E+BD3VlwOOBVwMR7r5Wlld26N4rLXoSSC5AW62AAgoooEB1AU8CmbwA+GxgK2BT4ANdAXBf4CHA3wM/AZ4B/C3wSuDTpfveChwGxHO/CxwLHATsCNzao0MNgNXftr6gAgoooIACOQED4OQFwE5H7A18oSsA9uqW84H/Ao4o3/wesAo4vXy9MXAtsBL4qAEw94azWgEFFFBAgTYIGACnOwDG7N0VwJHAR8ou35uBPYEvNxp0DXAZ8BYDYBvetq6DAgoooIACOQED4PQGwDg2MHb7xv9/uxwnuB3wQ2An4MpGa51TrhH0GgNg7g1ntQIKKKCAAm0QMABOZwCME0Bi1+99gOc3ju2LGUFnANvwznQdFFBAAQUUGKGAAXD6AuADgQvKSSAvAm7v6q9exwBeV44RnPUYwBUrVrB06dINi1q+fPmGPz4UUEABBRRQoD0Ca9asIf7EY/369axevTr+6p1A2jNEqTXZqOzWjZNALgQ2B+6IsQYeCnwOiEu6vKL8e/eLxXF+cRbwc4AIg8eU5z7Ws4BT42KxAgoooIACrRFwBnDyZgAPBs4G7ipdFhe6jr/vUy77Epd16VzOpfOcL5XA12nM44DXlvD4Va8D2Jr3qyuigAIKKKDAUAQMgJMXAIfSGAMsxOsADoDlUxVQQAEFFGiDgAHQAJjtQwNgVtB6BRRQQAEFKgsYAA2A2ZYzAGYFrVdAAQUUUKCygAHQAJhtOQNgVtB6BRRQQAEFKgsYAA2A2ZYzAGYFrVdAAQUUUKCygAHQAJhtOQNgVtB6BRRQQAEFKgsYAA2A2ZYzAGYFrVdAAQUUUKCygAHQAJhtOQNgVtB6BRRQQAEFKgsYAA2A2ZYzAGYFrVdAAQUUUKCygAHQAJhtOQNgVtB6BRRQQAEFKgsYAA2A2ZYzAGYFrVdAAQUUUKCygAHQAJhtOQNgVtB6BRRQQAEFKgsYAA2A2ZYzAGYFrVdAAQUUUKCygAHQAJhtOQNgVtB6BRRQQAEFKgsYAA2A2ZYzAGYFrVdAAQUUUKCygAHQAJhtOQNgVtB6BRRQQAEFKgsYAA2A2ZYzAGYFrVdAAQUUUKCygAHQAJhtOQNgVtB6BRRQQAEFKgsYAA2A2ZYzAGYFrVdAAQUUUKCygAHQAJhtOQNgVtB6BRRQQAEFKgsYAA2A2ZYzAGYFrVdAAQUUUKCygAHQAJhtOQNgVtB6BRRQQAEFKgsYAA2A2ZYzAGYFrVdAAQUUUKCygAHQAJhtOQNgVtB6BRRQQAEFKgsYAA2A2ZYzAGYFrVdAAQUUUKCygAHQAJhtOQNgVtB6BRRQQAEFKgsYAA2A2ZYzAGYFrVdAAQUUUKCygAHQAJhtOQNgVtB6BRRQQAEFKgsYAA2A2ZYzAGYFrVdAAQUUUKCygAHQAJhtOQNgVtB6BRRQQAEFKgsYAA2A2ZYzAGYFrVdAAQUUUKCygAHQAJhtOQNgVtB6BRRQQAEFKgsYAA2A2ZYzAGYFrVdAAQUUUKCygAHQAJhtOQNgVtB6BRRQQAEFKgsYAA2A2ZYzAGYFrVdAAQUUUKCygAHQAJhtOQNgVtB6BRRQQAEFKgsYAA2A2ZYzAGYFrVdAAQUUUKCygAHQAJhtOQNgVtB6BRRQQAEFKgsYAA2A2ZYzAGYFrVdAAQUUUKCygAHQAJhtOQNgVtB6BRRQQAEFKgsYAA2A2ZYzAGYFrVdAAQUUUKCygAHQAJhtuXsFwM98Bo46amax++8P5503WX9ftw6WLIHNN5/M7ZvUcctsl2M+We/hfnrBMZ+OMZ/mcT7mmHUceOCWMdDxn3XZMDCO9UvGcaVbtM73CoC77AKXXz6zhptsArfdNnl/7/hP6va5XffuWcd8ct/Ps/W7Yz4dYz6t47xs2TquuMIA2KI8NXarYgCcwIBrADQATuIvboP29bQGg0Gdxv350zrOBkB3AWcT5z0CYOz+PfRQuOEG2HprOOQQOPvsma9jl+mmm8Jdd8EvfgG33DLznL32gvPPn1mN/faDSy5p9/NvvRXWroWNNoLdd4evfGV81r22fT/j389zFrNHYnRvvnlmt//DHta+fm2j32WXQexaG9f3uWPe/s/h7GdZ531z3XUz7+0tN0yEje/PpUE94ufxsceu4yUvcQYwG4Kmuf4eAbC5+3fnneGb34Tufwuszi7ieM44f938zXfct2UU69/P+PfznMXskbaPcRv92rhOg/SQYz7en8v9fJZ192jbx3yQ/u13+z0JZPJmAA8EVgBPAB4A3Be4s5FQdwXeB+wWExvAmcDxXQk2vn41EOHua2V55ai+e2VdA2DZBdzPm27Yb+K2L6+fINDPcxZzO9v+g6GNfm1cp0F6yDE3AA7SL+P42R/vUQPg5AXAZwNbAZsCH+gKgBEIrwLOAk4AdgQuBN4DvLdEu7cChwH7At8FjgUOKs+9tcdU568C4EUXbXGP3b9nnDGzu3e2XcLdu4jH7euYQm/uMmjuihu3bRnF+vYz/v08Z7FdZ9v9PAqzQbe1jX5tXKdBXNv+vh533za+b9o+5oP0b7++554Ld9zhSSCTehbw3sAXugLgwcApwDaNWcHDgTcAO5Rw9z1gFXB6+Xpj4FpgJfDRuQLgXnttcY9du+M+CzDfb3VtnyVY7N9g+xn/fp6z2NvRfP22jXkb/dq4ToP0UNvGuHvdx913vs/VGt+f9t2/d/dUXPnFYwAn8Ri+XgEwgt1OZXavs817AheXLtio7BaOf/tyA2UNcBnwFgNg70vc1PjQGuSHWBvWp58fVP08p03b3bZw0Ea/Nq7TID3UtjE2AA5/d7QBsPNzzAA4TTOAsUt4M+BljSC3LM7JALYHIgD+sITEKxvPOadcJPI1g+wCPvhg+Ku/gmtj/jCmHbeBPfeEv/mbmHqGBz0I7n9/iLNqO2cEj8P3Z6bOZ85qXrp0vNa9pv1DHzpzoewPfejus8K7x7ef5yx2j8RhDHEW8MYb333Wd1v6t41+4/6+jysVxGdSWz+z2jjm4/C53fzs6zaM0N8Z83H9uTToZ3t8rm22mbuApykAjnQGcIsttrjHGb/N36Q7wXG+fxu37y/0ArK9PMZt22Mb5tv+7pnIXtvYz3Pme51Rf7/N/TsOfuPY245574v4T8pnV/f7pntyY9SfKYu9/Nj+Sy9dx5Yz17/xTiATth+41y7gOJnj1K5jAN9YTvqY6xjA64Aj5joGcMWKFSxdupQPfxhuvHE5sPwe4aDNH6aZN2Kmtp8ANe7L7yec9POcxXZoc/+Og58B8J4/XbL9PG5j3sbxn94AGEd0reHBD4aXvnQ9q1evNgBOUPiL3bhx6ZcIgHGGb5yregewvuz+jV27cRbwSeXEjwuA0xpnAcdxfnEW8HOAOCHkGOAVwGOBOc8CjhnA5n2ADzgAYndp816Lc/1bjME4fL9zYevO2VaxjeOy7rONx6jW/5RTZs4E79wbutf49vOcTt2o1nOu5cf6d+793Mb+bLtfHCrRRjfHfObzeTHeU82eWKzXb75vOj+jer3XF2v9Rv2ZF9v/tKc5Azhpu4DjTN+zgbtKqI3ti7/vA1wE7AKcUa4DuBZ4P3BiVwA+DnhtCY9fHeQ6gBMUpN0UBRRQQAEFJlbA6wBO3nUAazfrve4FXHsFfD0FFFBAAQUUGEzAAGgAHKxj7v1sA2BW0HoFFFBAAQUqCxgADYDZljMAZgWtV0ABBRRQoLKAAdAAmG05A2BW0HoFFFBAAQUqCxgADYDZljMAZgWtV0ABBRRQoLKAAdAAmG05A2BW0HoFFFBAAQUqCxgADYDZljMAZgWtV0ABBRRQoLKAAdAAmG05A2BW0HoFFFBAAQUqCxgADYDZljMAZgWtV0ABBRRQoLKAAdAAmG05A2BW0HoFFFBAAQUqCxgADYDZljMAZgWtV0ABBRRQoLKAAdAAmG05A2BW0HoFFFBAAQUqCxgADYDZljMAZgWtV0ABBRRQoLKAAdAAmG05A2BW0HoFFFBAAQUqCxgADYDZljMAZgWtV0ABBRRQoLKAAdAAmG05A2BW0HoFFFBAAQUqCxgADYDZljMAZgWtV0ABBRRQoLKAAdAAmG05A2BW0HoFFFBAAQUqCxgADYDZljMAZgWtV0ABBRRQoLKAAdAAmG05A2BW0HoFFFBAAQUqCxgADYDZljMAZgWtV0ABBRRQoLKAAdAAmG05A2BW0HoFFFBAAQUqCxgADYDZljMAZgWtV0ABBRRQoLKAAdAAmG05A2BW0HoFFFBAAQUqCxgADYDZljMAZgWtV0ABBRRQoLKAAdAAmG05A2BW0HoFFFBAAQUqCxgADYDZljMAZgWtV0ABBRRQoLKAAdAAmG05A2BW0HoFFFBAAQUqCxgADYDZljMAZgWtV0ABBRRQoLKAAdAAmG05A2BW0HoFFFBAAQUqCxgADYDZljMAZgWtV0ABBRRQoLKAAdAAmG05A2BW0HoFFFBAAQUqCxgADYDZljMAZgWtV0ABBRRQoLKAAdAAmG05A2BW0HoFFFBAAQUqCxgADYDZljMAZgWtV0ABBRRQoLKAAdAAmG05A2BW0HoFFFBAAQUqCxgADYDZljMAZgWtV0ABBRRQoLKAAdAAmG05A2BW0HoFFFBAAQUqCxgADYDZljMAZgWtV0ABBRRQoLKAAdAAmG05A2BW0HoFFFBAAQUqCxgADYDZljMAZgWtV0ABBRRQoLKAAdAAmG05A2BW0HoFFFBAAQUqCxgADYDZljMAZgWtV0ABBRRQoLKAAdAAmG05A2BW0HoFFFBAAQUqCxgADYDZljMAZgWtV0ABBRRQoLKAAdAAmG05A2BW0HoFFFBAAQUqCxgApzMAPghYBfwOsAnwLeCPgItK/+0KvA/YDbgZOBM4fpbeNABWftP6cgoooIACCmQFDIDTGQD/Gtge2A+4CTgCOA54OPBL4CrgLOAEYEfgQuA9wHt7NJwBMPsutF4BBRRQQIHKAgbA6QyAlwF/WWb5ouU2A24B9gB2Ak4FtgHuLP14OPAGYAcDYOV3qC+ngAIKKKDACAQMgIsfAL8P3NXH2D66j+f0+5TYnftM4MXAjcBK4FVA7Pp9ZwmB+zYWtidwMbAl8LOuF3EGsF91n6eAAgoooEBLBAyAix8A/6DRC7Fb9lDgw0AEw0cBrwTOmOMYvIW00v2A84HlZZfvT8vu4H8GPlBmBF/WWPAy4PKy2/haA+BCyK1RQAEFFFCgPQIGwMUPgM1u+DzwJ0AEsc7jKcA7gGcNsW2+VI7zi5m/2PX73BI6nw4c4gzgEKVdlAIKKKCAAi0UMAC2KwCuA+IM3TsavbJxOVEjdrUO4/Fg4L+BJwL/3ljg14BzgOt7HAP4RuCwuY4BXLFiBUuXLt2wuOXLl2/440MBBRRQQAEF2iOwZs0a4k881q9fz+rVq+OvcXhX5I+peyxp0Rb/RznT9oONdYpj8+Is3ccPcT2/B8Rs45vLMX3PAT4J/C7wr8CV5Szgk0rouwA4zbOAhzgCLkoBBRRQQIFFFHAGsF0zgDFt9rdAnKUbxwA+sgS/uFzLZ4fYJ48rs3xx1u99gauBPwM6wXOXctxhXAdwLfB+4MRZXt+TQIY4MC5KAQUUUECBGgIGwHYFwBjzOPEjTsDYDrgG+FgJgzX6YSGvYQBciJo1CiiggAIKLKKAAbB9AXAR22FBL20AXBCbRQoooIACCiyegAFw8QNgXOblI6UF4ni/2R5xZ442PgyAbRwV10kBBRRQQIE5BAyAix8AvwnEMXfxiOP+ej3iQtHDvBD0MN8UBsBharosBRRQQAEFKggYABc/AFYY5pG+hAFwpLwuXAEFFFBAgeELGAANgNmuMgBmBa1XQCDXib4AACAASURBVAEFFFCgsoAB0ACYbTkDYFbQegUUUEABBSoLGAANgNmWMwBmBa1XQAEFFFCgsoAB0ACYbTkDYFbQegUUUEABBSoLGAANgNmWMwBmBa1XQAEFFFCgsoABsH0BcG9gd2Dzrl44pnJv9PtyBsB+pXyeAgoooIACLREwALYrAJ4AHAV8A/h5o0fiOoDPbEnPdK+GAbClA+NqKaCAAgooMJuAAbBdAfA6YD/gX8aoZQ2AYzRYrqoCCiiggAIhYABsVwD8CfAQIGb8xuVhAByXkXI9FVBAAQUUKAIGwHYFwPcBFwGfHKMONQCO0WC5qgoooIACCjgDONMDS1rUCh8DXghcAlzbtV4HtWg9m6tiAGzpwLhaCiiggAIKzCbgDGC7AuDZc7Tq77e0jQ2ALR0YV0sBBRRQQAED4Ow90KYZwHHsVAPgOI6a66yAAgooMNUCzgC2awaws0v6ycDDgR8AX235SSEGwKn+CHHjFVBAAQXGUcAA2K4A+DDg08CTgJuABwFfB17Q45jAtvSbAbAtI+F6KKCAAgoo0KeAAbBdAfCcclLK64GfAlsBq4GNgAP7HNPaTzMA1hb39RRQQAEFFEgKGADbFQDjQtCPjeszNsZ1S+BK4H8lx3pU5QbAUcm6XAUUUEABBUYkYABsVwD8MfAY4GeN8Y57An8H2HpEPZBdrAEwK2i9AgoooIAClQUMgO0KgB8HbgcOLbOAMft3OrAJ8JLKvdHvyxkA+5XyeQoooIACCrREwADYrgC4DXAB8PhyDGCcBPJN4LnAj1rSM92rYQBs6cC4WgoooIACCswmYABsVwCMcYoTPnYHtgeuBr4C3NniFjYAtnhwXDUFFFBAAQV6CRgA2xMAl5ZLvvwG8D9j1K4GwDEaLFdVAQUUUECBEDAAticAxnjEbt5HAevHqD0NgGM0WK6qAgoooIACBsCZHmjTreBOAG4DThqj9jQAjtFguaoKKKCAAgoYANsXAL8E7AFcX24D1zz27+ktbVkDYEsHxtVSQAEFFFBgNgF3AbdrBvDYOVr1+Ja2sQGwpQPjaimggAIKKGAAnL0H2rQLeBw71QA4jqPmOiuggAIKTLWAM4DtmgGMS8DM9mjrpWAMgFP9EeLGK6CAAgqMo4ABsF0BMELeXbM00sYtbTADYEsHxtVSQAEFFFDAXcDjsQt4767V3BZ4M3Am8OctbWMDYEsHxtVSQAEFFFDAADgeAbDXWj4a+DDw1Ja2sQGwpQPjaimggAIKKGAAHN8AGLt+bwIiaLXxYQBs46i4TgoooIACCswh4DGA7ToGMGb7mo/NgFcD+wC7trSTDYAtHRhXSwEFFFBAAWcAx2MGsPskkLhEzfeAg4BLW9rGBsCWDoyrpYACCiiggAFwPALgI7pW8xbgpy1vXwNgywfI1VNAAQUUUKBbwF3A7doFPI4dagAcx1FznRVQQAEFplrAANiuABgXgj4SeBXwUGBLYDnw8HIpmDY2qwGwjaPiOimggAIKKDCHgAGwXQHwBOB5wLuBM4AHAr8OfBJ4Uks72QDY0oFxtRRQQAEFFJhNwADYrgD4feDpwNXl2L+tgJgV/AkQf2/jwwDYxlFxnRRQQAEFFHAGcM4eiDNt2/KIoPeQcju4OPkjQt99gR+VXcJtWc/mehgA2zgqrpMCCiiggAIGwLEJgP9Ydv3GLt9OAHwR8IflWMA2NrMBsI2j4jopoIACCihgABybALgH8Hng08B+wP8FIgA+G/jXIXfynsA7gCcDdwCXN243Fxedfh+wG3BzOQHl+Fle3wA45IFxcQoooIACCoxawGMA23UMYIz3Y4FDgR2AH5cZwa8OuREi/P09cBjwKeD2EvYiZD4AuAo4C4iTUnYELgTeA7y3x3oYAIc8OC5OAQUUUECBUQsYANsXAEc95rH8i4AvA2/t8WIHA6cA2wBxZ5J4HA68oYTS7hIDYI0R8zUUUEABBRQYooABsH0BMO4HHLteN+8a55iRG8bj/kDcYeQ04BnlMjNx9vE7gfOAVcBOwL6NF4sZw4vLdQl/1rUSBsBhjIrLUEABBRRQoKKAAbBdAfB1wOnlBJCfN/rgLiCC4TAe25bLzFwPPAf4BvAC4Bxgb+DVwGbAyxovtqwcI7g9cK0BcBjD4DIUUEABBRRYPAEDYLsC4A+AI8pM3Ki6Imbs4sSOdwFHN17ks8C/AZssZAZwxYoVLF26dMPili9fvuGPDwUUUEABBRRoj8CaNWuIP/FYv349q1evjr/GXcfWtWct661Jm64DGMEs7v4x6se3y91FegXAbwGndh0D+MZywkicmNL9cBfwqEfL5SuggAIKKDBkAWcA2zUD+FHgL4EvDnmcuxcXJ3UcVY7z+49y+7nYBRx3IbkCuLKcBXxSOfHjgnLMoGcBj3hgXLwCCiiggAI1BAyA7QqAfwbEWbjn9jjW7pghN8SRwIoy9RszgscBnymvsUu5/EycjLIWeD9w4iyv7wzgkAfGxSmggAIKKDBqAQNguwJg3Amk1yNOAnnmqJthgcs3AC4QzjIFFFBAAQUWS8AA2K4AuFh9kHldA2BGz1oFFFBAAQUWQcAAaADMtp0BMCtovQIKKKCAApUFDIAGwGzLGQCzgtYroIACCihQWcAAaADMtpwBMCtovQIKKKCAApUFDIAGwGzLGQCzgtYroIACCihQWcAAaADMtpwBMCtovQIKKKCAApUFDIAGwGzLGQCzgtYroIACCihQWcAAaADMtpwBMCtovQIKKKCAApUFDIAGwGzLGQCzgtYroIACCihQWcAAaADMtpwBMCtovQIKKKCAApUFDIAGwGzLGQCzgtYroIACCihQWcAAaADMtpwBMCtovQIKKKCAApUFDIAGwGzLGQCzgtYroIACCihQWcAAaADMtpwBMCtovQIKKKCAApUFDIAGwGzLGQCzgtYroIACCihQWcAAaADMtpwBMCtovQIKKKCAApUFDIAGwGzLGQCzgtYroIACCihQWcAAaADMtpwBMCtovQIKKKCAApUFDIAGwGzLGQCzgtYroIACCihQWcAAaADMtpwBMCtovQIKKKCAApUFDIAGwGzLGQCzgtYroIACCihQWcAAaADMtpwBMCtovQIKKKCAApUFDIAGwGzLGQCzgtYroIACCihQWcAAaADMtpwBMCtovQIKKKCAApUFDIAGwGzLGQCzgtYroIACCihQWcAAaADMtpwBMCtovQIKKKCAApUFDIAGwGzLGQCzgtYroIACCihQWcAAaADMtpwBMCtovQIKKKCAApUFDIAGwGzLGQCzgtYroIACCihQWcAAaADMtpwBMCtovQIKKKCAApUFDIAGwGzLGQCzgtYroIACCihQWcAAaADMtpwBMCtovQIKKKCAApUFDIAGwGzLGQCzgtYroIACCihQWcAAaADMtpwBMCtovQIKKKCAApUFDIAGwGzLGQCzgtYroIACCihQWcAAaADMtpwBMCtovQIKKKCAApUFDIAGwGzLGQCzgtYroIACCihQWcAAaADMtpwBMCtovQIKKKCAApUFDIAGwGzLGQCzgtYroIACCihQWcAAaADMtpwBMCtovQIKKKCAApUFDIAGwGzLGQCzgtYroIACCihQWcAAaADMtpwBMCtovQIKKKCAApUFDIAGwGzLGQCzgtYroIACCihQWcAAaADMtpwBMCtovQIKKKCAApUFDIAGwGzLGQCzgtYroIACCihQWcAAaAA8H3gB8CzgC6X/dgXeB+wG3AycCRw/S28aACu/aX05BRRQQAEFsgIGwOkOgAcBvwc8u/yJAPgA4CrgLOAEYEfgQuA9wHt7NJwBMPsutF4BBRRQQIHKAgbA6Q2A2wEXA08FftiYATwYOAXYBriz9OPhwBuAHQyAld+hvpwCCiiggAIjEDAATm8AXAN8AvhgCXqdXcCrgJ2AfRv9tmcJi1sCP+vqQ2cAR/DGdJEKKKCAAgqMUsAAOJ0B8NBy3N/y0lwx09cJgB8ANgNe1mi8ZcDlwPbAtQbAUb4lXbYCCiiggAKjFzAATl8AfHSZzdsDuLpHAFzQDOCKFStYunTphsUtX758wx8fCiiggAIKKNAegTVr1hB/4rF+/XpWr14df429e+vas5b11mRJvZdqxSvFMX5/UQa7s+0PBtYCHwcuBd4NPKxxDOAbgcM8BrAV4+dKKKCAAgookBZwBnD6ZgDvB2zV1TnXAAcCnwN+CVxZzgI+qYS+C4DTPAs4/X5zAQoooIACCrRCwAA4fQGwV+Pd0bgMTHx/F+CMch3AmBl8P3DiLB3rSSCteCu7EgoooIACCvQvYAA0APbfLb2faQDMClqvgAIKKKBAZQEDoAEw23IGwKyg9QoooIACClQWMAAaALMtZwDMClqvgAIKKKBAZQEDoAEw23IGwKyg9QoooIACClQWMAAaALMtZwDMClqvgAIKKKBAZQEDoAEw23IGwKyg9QoooIACClQWMAAaALMtZwDMClqvgAIKKKBAZQEDoAEw23IGwKyg9QoooIACClQWMAAaALMtZwDMClqvgAIKKKBAZQEDoAEw23IGwKyg9QoooIACClQWMAAaALMtZwDMClqvgAIKKKBAZQEDoAEw23IGwKyg9QoooIACClQWMAAaALMtZwDMClqvgAIKKKBAZQEDoAEw23IGwKyg9QoooIACClQWMAAaALMtZwDMClqvgAIKKKBAZQEDoAEw23IGwKyg9QoooIACClQWMAAaALMtZwDMClqvgAIKKKBAZQEDoAEw23IGwKyg9QoooIACClQWMAAaALMtZwDMClqvgAIKKKBAZQEDoAEw23IGwKyg9QoooIACClQWMAAaALMtZwDMClqvgAIKKKBAZQEDoAEw23IGwKyg9QoooIACClQWMAAaALMtZwDMClqvgAIKKKBAZQEDoAEw23IGwKyg9QoooIACClQWMAAaALMtZwDMClqvgAIKKKBAZQEDoAEw23IGwKyg9QoooIACClQWMAAaALMtZwDMClqvgAIKKKBAZQEDoAEw23IGwKyg9QoooIACClQWMAAaALMtZwDMClqvgAIKKKBAZQEDoAEw23IGwKyg9QoooIACClQWMAAaALMtZwDMClqvgAIKKKBAZQEDoAEw23IGwKyg9QoooIACClQWMAAaALMtZwDMClqvgAIKKKBAZQEDoAEw23IGwKyg9QoooIACClQWMAAaALMtZwDMClqvgAIKKKBAZQEDoAEw23IGwKyg9QoooIACClQWMAAaALMtZwDMClqvgAIKKKBAZQEDoAEw23IGwKyg9QoooIACClQWMAAaALMtZwDMClqvgAIKKKBAZQEDoAEw23IGwKyg9QoooIACClQWMAAaALMtZwDMClqvgAIKKKBAZQEDoAEw23IGwKyg9QoooIACClQWMAAaALMtZwDMClqvgAIKKKBAZQEDoAEw23IGwKyg9QoooIACClQWMAAaALMtZwDMClqvgAIKKKBAZQEDoAEw23IGwKyg9QoooIACClQWMABOZwB8J/Ac4BHAz4AvAm8Drmn0367A+4DdgJuBM4Hje/SnAbDym9aXU0ABBRRQICtgAJzOAHgS8CngMmBT4P3A44AnloZ6AHAVcBZwArAjcCHwHuC9XU1nAMy+C61XQAEFFFCgsoABcDoDYHebPQH4OrAVsBY4GDgF2Aa4szz5cOANwA4GwMrvUl9OAQUUUECBIQsYAA2A0VKx+/d1wKNLf60CdgL2bfTbnsDFwJZlt3HnW84ADvlN6eIUUEABBRQYtYAB0AD4LOB8YH/gc6XhPgBsBrys0YDLgMuB7YFrG/9uABz1u9TlK6CAAgooMGQBA+B0B8DnAh8pu3w/3eitgWcAV6xYwdKlSzcsYvny5Rv++FBAAQUUUECB9gisWbOG+BOP9evXs3r16vhr7Nlb1561rLcmS+q9VKte6eXA6cCLgc93rdlBwKldxwC+ETjMYwBbNYaujAIKKKCAAgsScAZwOmcAI8jF2b3PAy7p0TlxFvCV5SzgOGM4Tvy4ADjNs4AX9D6zSAEFFFBAgVYJGACnMwDGmb23A7eVboxZ0LvKSR+dQLgLcEa5DmCcGRyXijmxR/d6DGCr3tKujAIKKKCAAvMLGACnMwDO3xn9P8MA2L+Vz1RAAQUUUKAVAgZAA2C2EQ2AWUHrFVBAAQUUqCxgADQAZlvOAJgVtF4BBRRQQIHKAgZAA2C25QyAWUHrFVBAAQUUqCxgADQAZlvOAJgVtF4BBRRQQIHKAgZAA2C25QyAWUHrFVBAAQUUqCxgADQAZlvOAJgVtF4BBRRQQIHKAgZAA2C25QyAWUHrFVBAAQUUqCxgADQAZlvOAJgVtF4BBRRQQIHKAgZAA2C25QyAWUHrFVBAAQUUqCxgADQAZlvOAJgVtF4BBRRQQIHKAgZAA2C25QyAWUHrFVBAAQUUqCxgADQAZlvOAJgVtF4BBRRQQIHKAgZAA2C25QyAWUHrFVBAAQUUqCxgADQAZlvOAJgVtF4BBRRQQIHKAgZAA2C25QyAWUHrFVBAAQUUqCxgADQAZlvOAJgVtF4BBRRQQIHKAgZAA2C25QyAWUHrFVBAAQUUqCxgADQAZlvOAJgVtF4BBRRQQIHKAgZAA2C25QyAWUHrFVBAAQUUqCxgADQAZlvOAJgVtF4BBRRQQIHKAgZAA2C25QyAWUHrFVBAAQUUqCxgADQAZlvOAJgVtF4BBRRQQIHKAgZAA2C25QyAWUHrFVBAAQUUqCxgADQAZlvOAJgVtF4BBRRQQIHKAgZAA2C25QyAWUHrFVBAAQUUqCxgADQAZlvOAJgVtF4BBRRQQIHKAgZAA2C25QyAWUHrFVBAAQUUqCxgADQAZlvOAJgVtF4BBRRQQIHKAgZAA2C25QyAWUHrFVBAAQUUqCxgADQAZlvOAJgVtF4BBRRQQIHKAgZAA2C25QyAWUHrFVBAAQUUqCxgADQAZlvOAJgVtF4BBRRQQIHKAgZAA2C25QyAWUHrFVBAAQUUqCxgADQAZlvOAJgVtF4BBRRQQIHKAgZAA2C25QyAWUHrFVBAAQUUqCxgADQAZlvOAJgVtF4BBRRQQIHKAgZAA2C25QyAWUHrFVBAAQUUqCxgADQAZlvOAJgVtF4BBRRQQIHKAgZAA2C25QyAWUHrFVBAAQUUqCxgADQAZlvOAJgVtF4BBRRQQIHKAgZAA2C25TYEwI9/fC0nnBB/hXe9a2aRRx018//994fzzpusv69bB0uWwOabT+b2Teq4ZbbLMZ+s9/BsveA4T8c4N8d/Wsf8jjvWccUVWwZF/GddNgyMY/2ScVzpFq3zhgC4bNlarrhiJgDuvPPM2l1++cz/N9kEbrtt8v7eGYNJ3T63694965hP7vu52e+O83SMs2Memc8A2KI8NXarYgCc0IBrADQATuIvbv30tQHQADjJkxd3vwcMgM4A5jLnhgD4tret5YILZmYATzkFvvxlOPXUmQXvtx9ceilcf/3MLtNNN53591tvhVtuga23hr32gvPPH4/nx67fm26aWddtthmvdR93+8Vc/zvumOndCEWz9fEhh8CHPjQ5vb6Y3ov12RGfSzffDBtvDLvvDl/5ynh8Lk3jWA3r50g/7+1x+hnVTy9cdhm4C9hjAHPxD341A/itb80EwHjsssvdu4C7dwmP+9fNGYRx3xbXf6ZfO4crzOcx39h/85uT3fvz+Uza9+cb70F6Z9JsJm17pm2s47PKk0AMgHMFwOOBVzMT8r4GrIiflV0FGwLgttuu5YEPvPskkDgBpHkMYMzy/fCHd+9SG+evYwbwrrsmY1viQ2/UYxFSzV2Jza+33x5uuOHu77f96+bY91rXLbaA73xn9u0N77k8Bv1+27wiFMSsfue9/vCHz2zvuH4ds4AxOxSPNmxLzDxP6udqjc+i7s+65nux+3N9mO/TQd/XNZ7/mMfAMces48ADPQYwOws2ifVvBQ4D9gW+CxwLHATsGHtvGxu8IQBut91arrnm7pNA4kzgF71ofH6wD/qDdL4gsHihZg3bb7980UNVP0EoPoCaYWExfgAMEn4HDf533rmG229fvuGtMooP9Db6nXsuvP71cPXVM58Q8b6KRzu+nnlvDLI+G20Ed97Znm0544xJ+VxdAyzf0B+1Pyu7P5umLfg1g+2yZZ4F7DGAvePr94BVwOnl2xsD1wIrgY92B8Dus4AnfVdYe3cXrGTnnWPY+t+tOYpdOf2Mfz/PWeztaL7+oGP+4x+v5MYbZ8ZisYxr+7V7TAd/bww65qP2brfvIH0eP0ZWjeR9Md8YdBu2bYznW/9hft8A6C7gXvEvpvJuBvYEvtx4Qvzadhnwlu4A2LkOYK/rKXX+LXabdn7bGvf//+IXM7NXnZMBFrptUX/AAXD22TO/CceMVJxIsNCvly5dyRFHrFpwffb1oz5mKcLj0ENn36Z+npNx6LUd8ZrNk5OOPHJmPTsnK833dfOEgPmeG99fs2Yl3/jGTADs5/mDrEss78lPntt4Mfxqr9N829gc81/+ciVvf/uqkY75oGM46PNr+2bfM7P1/cknr+Q+91k1kvfFfKbdht0nd1xyyXA+izufhZnPnFF8bjR93vSmdZxyiruAe8+BTe+/bheH7QA7AVc2GM4pF4t8TXcAvPrqq9liiy3YYw+44oqZ7y5bNvMDt/vf4nvN54zz10uXwvr1d2/voNsybJ8HPehott765EX17Web+nnOoJbzPT/7moOO9fXXH81NN538q/fCfOvX9u9n/eLzoPY23nOdj2bZspnx6PfzZ9AxH2TZi+Ex6PqNbsyPBk7e8DNi0HUa5fOHPd6j81v4z5ym3w47rOPb395wnIYXgp7evHevLR9kBnBb4BrtFFBAAQUUUGAsBWLS50djuebJlfYYwN6AvY4BvA44ousYwPDbBrglOQ6WK6CAAgoooEBdgc3L8f3l2hZ1X3yxX80A2HsE4ji/OAv4OUCEwWOAVwCP7ToLeLHHz9dXQAEFFFBAAQUGFjAAzk52HPBaIH5D+Oos1wEcGNwCBRRQQAEFFFBgsQUMgIs9Ar6+AgoooIACCihQWcAAmAPv524huVew+sAy+/oE4AHAfYFyedoNOLsC7wN2K5fvOROIcWk+HKfh9NE7y2ERjwB+BnwReFvXiVCOx3Cs51tKHJYSF6f/NSDu1xF3KToB+Hyj0LGYT3E03487u78AeBbwhfISjsVorLuXGjdt+JNyqFbkmzi27++AlzsO9x4AA+DCm7Lfu4Us/BWsDIFnA1sBmwIf6AqAEQivAs4qP/ziTi0XAu8B3lv4HKfh9dFJwKfK9TBjPN4PPA54YnkJx2N41vMtaQfghrgTEXAf4PAN1xaBBwM/L78s+d6YT3H4349Q/nvlcys+uyIA+r4YvvNsS4wA+FvA03s8wXHoQjEALrwx+71byMJfwcqmwN7lw7Q5A3gwcEo5E7szKxg/CN8AxA/IeDhOo+ujmJX9egnoEUQcj9FZz7XkuMvy64E3lb6/3bFYlIGIy4lcDDy1XEu2MwPo+6LecMwVAB0HA+BQOnGQawUO5QVdCL0CYNxqIi7YHfds7jziDi7xIRwX99xogLu6SDy4QOz+fR3w6FLqeAxumKn43XJZquj1H5Td8//pWGRIU7Vxt6hPAB8sh6l0AqDvixTrQMURAOMqHreWP5cCbwf+q9ze1Z8XDU5nAAfqrV89eZC7hSzsFazqFugVAGOX8GbAyxpPjuvrx/FQcYn3CID93tVF8cEE4odbHOu0P/C5Uup4DGY4rGc/EDgReF7cernsAnYshqXb33IOLcf9LS9Pjz0SnQDoWPRnOIxnxSEpcV3eq4GHAe8G9gBib8X/8efFPYkNgAtrOWcAF+aWqXIGMKM33NrnAh8puxk/3Vi0Mx3DdR5kafFZ3tkNH8HcsRhEL/fcmAGPvQ4RNCJ4xKMZAB2LnG+meml5X8QvRzFj7gxgQ9MAuPDW6vduIQt/BSubAr0CYBxwfWrXMYBvLBfxnusYwF53dVG7P4E4m+504MVdZ5xGtePRn+EonhUngkQAfGGZkXUsRqHce5lxbNlflHvFd36mxsk4MR4fB2I3ZMxExYxU51hlP6fqjE8nAD6/+PvzwgA4lM7zbiFDYZx3IbEbN078iAAYZ/jGhbnjshfry3T+leUs4DhDNULfBcBpjbOAHad5ift+QtwdJy41Er9NX9KjKs6yczz65kw9MU52OqecCfwQIPo/zn58fDn2ybFI8Q5UfL9yIlSzKO4RH5ewisMjfun7YiDPzJPjF9M48/pGYOsSvOOknHhfRDj388kAmOmve9R6t5ChUc66oPjt+uxyPad4UufaTvsAFwG7AGeU6wDGb9xxaZI4Hqr5cJyGM04xexFnmN5WFtcZizgJpxMIHY/hWM+3lLi22W+US4zcVK7JGNcG/H6j0LGYT3F0349fUjuXgYlXcSxGZ91c8t8CTymTA/G+iJ8RcV3A2GPnOHSNgbuA6zSlr6KAAgoooIACCrRGwADYmqFwRRRQQAEFFFBAgToCBsA6zr6KAgoooIACCijQGgEDYGuGwhVRQAEFFFBAAQXqCBgA6zj7KgoooIACCiigQGsEDICtGQpXRAEFFFBAAQUUqCNgAKzj7KsooIACCiiggAKtETAAtmYoXBEFFFBAAQUUUKCOgAGwjrOvooACCiiggAIKtEbAANiaoXBFFFBgCAJ/VG4b+DtDWJaLUEABBSZWwAA4sUPrhikwcQK3NG4JuAmwcbnvba9b0tXY+FcDLwSeW+PFfA0FFFBgmAIGwGFquiwFFKglEPd7/k3gmbVesMfrXAh8HPjQIq6DL62AAgosSMAAuCA2ixRQYJEFZguAxwLPAp5W1u8fgcuAhwL7AjcDbwWuAt4P7AxcDhwEfLvUbAQcAfwBsE359yOBLzS2eUvgOmBbIG46//0SBPcsj0Y+cAAAA2pJREFUwfRHwOvKLOUq4JHAvwCvBG4oyzkMeBPwkDKTGYHyVYvs6ssroMCUCBgAp2Sg3UwFJkxgrgD4W8DTGwHwCcDzgEuAw4GTgS8Brykh7hxgsxIQo+w44PnAS4DvAC8APgo8vgS9eM4rgENK2IyvIwDGI14nAuW7gZeX0BdB8jbgH4D/AF4PPAb4d+A3gG8BmwJPAi6esHFycxRQoKUCBsCWDoyrpYACcwoMEgAjxP1hWdoWZRbwpcAnyr/tD5wJPLh8HbOELwY+11iDCG//VMJj/PP5JdDFLGInAMYyIlzGY1fg34CYEfxK+beVJRTuVmYEIyhGiPwsEMc3+lBAAQWqCRgAq1H7QgooMESBQQJgzPYdU147Thy5HXgGcFH5t+XA3wFLy67iHwPrgDvL9+Nz8j7Ah4EVZbYuduPGLF48txMAY53OKl//etnN/Cjgh+XfXgu8GdixfB2zhYcCTynP/VMgZiN9KKCAAiMXMACOnNgXUECBEQiMKgBGCIxj+iIUzrY79kXl2L2nNrYrdgEPGgA75XHM4QEl/D227HYeAZmLVEABBe4WMADaDQooMI4CowqAYXFa2XUbJ2RcAdy/HKsXJ33E7uSPld26MWPXefQKgHFSSZz80WsGMGYBY5YwZiF/XgLnBcAOjeMJx3FcXGcFFBgTAQPgmAyUq6mAAvcQ6DcAxpm7MZPX3AW8Hthnll3A8SLxuRhn6MZJItsB/wN8HXgL8F3gv8sxfj9orNH3gHcMsAt4F+DPy1nI8XoREt9ZwqVDrYACCoxcwAA4cmJfQAEFJkggLvocZwnH2bs+FFBAgbEVMACO7dC54goosAgCcYmZOJEkzgr2oYACCoytgAFwbIfOFVdAAQUUUEABBRYmYABcmJtVCiiggAIKKKDA2AoYAMd26FxxBRRQQAEFFFBgYQIGwIW5WaWAAgoooIACCoytgAFwbIfOFVdAAQUUUEABBRYmYABcmJtVCiiggAIKKKDA2AoYAMd26FxxBRRQQAEFFFBgYQIGwIW5WaWAAgoooIACCoytgAFwbIfOFVdAAQUUUEABBRYmYABcmJtVCiiggAIKKKDA2AoYAMd26FxxBRRQQAEFFFBgYQIGwIW5WaWAAgoooIACCoytgAFwbIfOFVdAAQUUUEABBRYm8P8BYRyl/LREg2oAAAAASUVORK5CYII=\">"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "application/javascript": [
       "/* Put everything inside the global mpl namespace */\n",
       "window.mpl = {};\n",
       "\n",
       "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",
       "            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",
       "        this.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 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);\n",
       "        canvas.attr('height', height);\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'];\n",
       "    var y0 = fig.canvas.height - msg['y0'];\n",
       "    var x1 = msg['x1'];\n",
       "    var y1 = fig.canvas.height - msg['y1'];\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;\n",
       "    var y = canvas_pos.y;\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",
       "    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 + '\">');\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 dataURL = this.canvas.toDataURL();\n",
       "    this.cell_info[1]['text/html'] = '<img src=\"' + dataURL + '\">';\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,iVBORw0KGgoAAAANSUhEUgAAAoAAAAHgCAYAAAA10dzkAAAgAElEQVR4Xu3dDbxldV3v8c+AzhgwTFHqDRxNiwdl1NRKuSrmUxOBTzOImYY3MzXGhySsia48mQgo3ExAC4OuXq9Ygq8UXzTp1TS1LLRbMDewxARFIwrnKOQMDHNfvzn/7Sy252Gf89v7P2vt9dmvF8w5c/ZvPbz/v73Pd9Za/7VX4EMBBRRQQAEFFFCgVwIrerW37qwCCiiggAIKKKAABkCbQAEFFFBAAQUU6JmAAbBnA+7uKqCAAgoooIACBkB7QAEFFFBAAQUU6JmAAbBnA+7uKqCAAgoooIACBkB7QAEFFFBAAQUU6JmAAbBnA+7uKqCAAgoooIACBkB7QAEFFFBAAQUU6JmAAbBnA+7uKqCAAgoooIACBkB7QAEFFFBAAQUU6JmAAbBnA+7uKqCAAgoooIACBkB7QAEFFFBAAQUU6JmAAbBnA+7uKqCAAgoooIACBkB7QAEFFFBAAQUU6JmAAbBnA+7uKqCAAgoooIACBkB7QAEFFFBAAQUU6JmAAbBnA+7uKqCAAgoooIACBkB7QAEFFFBAAQUU6JmAAbBnA+7uKqCAAgoooIACBkB7QAEFFFBAAQUU6JmAAbBnA+7uKqCAAgoooIACBkB7QAEFFFBAAQUU6JmAAbBnA+7uKqCAAgoooIACBkB7QAEFFFBAAQUU6JmAAbBnA+7uKqCAAgoooIACBkB7QAEFFFBAAQUU6JmAAbBnA+7uKqCAAgoooIACBkB7QAEFFFBAAQUU6JmAAbBnA+7uKqCAAgoooIACBkB7QAEFFFBAAQUU6JmAAbBnA+7uKtASgXuAM4CzyvbE16cB9wHiZ21/DG//3t7e04vfvotsyJHAa4DHAY8s3ovV7O19c/0KKDABAQPgBFBdpAIKLCrwU8BXgVvKMwcB5r4GwEXt5nrCwcCDgL9ZpPpE4I3ANcADgaMAA+CyyC1SoNsCBsBuj59br8C0CBgA649kBMFTDYD14V2jAm0QMAC2YRTcBgW6IXAocB7wX4EDgVuBvwZeWI7aPQX4BHA88CzgOSVcfLicdvyPxm4On0KdKwD+LPAnwGWlPsq/r5w6fj5wCPA14F3Am4FdZfn7A+cAzy5HubYBfw+8CvjimKiHt39Um+cBP1eM4v33fwK/DjweOB94NPBl4GTgz5ewrYNT6PssocYAuAQsn6rAtAkYAKdtRN0fBSYn8E/Av5dwFX9GAIsw81LgbmAQAG8GPga8H4hgFOHsb4GnLyEAxqnKS0rYi/p4xKnKvwCOKNcOXgc8oVz7diHw+vK8qDsO+C3gn4EfBJ4IfGCEU6Sj6g0HwFFt/gW4ErgaOBp4A/D7wFOLU5wSj797DPAQoBmaF9q2ua4B/CMgHOcLhQbAUUfb5ykwhQIGwCkcVHdJgQkIRIj6t3JU7ap5lj8IgBFujm085xeA9wDPKEcI40cLHQE8pVyn9spy9G+wqF8EItREcPpMY/lxGjMmkMQ1cLcB1wJbgFjOpB7N7V+KzR8Cv9LYqM8DPw48Cfir8vcxOSOOWL6kuI2yD3MFwDgyGmar5lmAAXAUWZ+jwJQKGACndGDdLQUmIBBH074D/G45EhffNx+DABhHBCOoDR4xseM/S0g7u/zlfAHwIuBlwAnAcND8X+X0cxxVbD4eC3yuEU4vLV/HKdU4jfp3I0wsiffC5pGy2L7BKeW5KIe3f1SbOHV9RWOB7y2ny+OUetNrO/DfgYHXYsM56izg5nIMgIup+nMFpljAADjFg+uuKTBmgR8pp2TjtO8PlWvV3gK8s6xnEADj6F8cBWw+4ujhHwObFgiAEWLiCF5cA/fTJTQ2lxFhLo4izvWIsBbBM66pi2sA46hghK0fBW4H3g389hzLHCwrrjOMI27xiGXFcmJ58z2GA+CoNrH9H28sNNYbp8YfPLSiWP7vlNA8yjAaAEdR8jkKKPBdAQOgzaCAAssReFSZVBFH644pp1zHcQQwTod+FLihLPfOxsa9D/jJEuzmeu+K6+uGr5lbWyZcnFsmWcR1gXM9IoBFqB08IojetIQA2HzqQjYGwOV0mzUKKDB2AQPg2EldoAK9EVgNxAzbmHwRp1sHAfDPyuSQAcSLyhG4Ua8BPKwcJYvTqhEu7ygLiiN0f1BuYLzU2bxxrV1MsIjZyeN4LHYj6Lls4sjfMz0COA5+l6GAAlkBA2BW0HoF+iEQExPeVmb2RjCLGbm/BGwoM3HjOrvmLOD/A1wOHF5OZcYs4Obp28VuAxPX+cWM3zgdHLeD+Xb51Io4Ohg/i8AZEyVWAj/WuO1MXKP4WeBDZTJI1MXp5Lie7nVAzBYex6O5/UuxqXkEMCacxCzguAZz8Ijb6MQp/HjE7XvitjRxvWU84ghqBGUfCijQAwEDYA8G2V1UYAwC9y/3AIxPjojZthG0Yrbtm8otX2IVEQDjKNfGMgnjuSUoRhh77dDp2Z3lesKYiBCPuIYtbn8SgW7wUXBx/V7cVzBuK7O+hMD4+Wbg54GHlqODXyoTRmJbojZuG/MzwMNKaLyxHDmMCSbjejS3fyk2cx0BfFq55Utz22L5cQ1guIzyGPjFR+kNHnF9YcwCbv5d3FomQvVcE1wWu+5xlO3wOQoo0BGBaQuA8cYfF6DHm1z8y/+TwG+Uj5waDEn8gohfXnHfstj/eCOMX2pbG2N2ZpmJGDPz4l/EceF68+cdGV43U4GqAoMAOBxyqm6EK1NAAQUUWFxg2gJgHAGIm73GkYn9gHcAjyg3VW0GwJh1F0cW5nrE9UzxiQFx7VEcWYh/WcdplLguqXlB+uK6PkOBfgkMTgEPn+bsl4J7q4ACCnRAYNoC4DB5fKzSF4CDysXq8fM4ArjQL6g4XXRB41qhuNYpLh6Pj2aKe3b5UECBuQU8AjiZzoj3oIUecbrYhwIKKLAkgWkPgHH6Nz5NIK4FGjwiAH6jXBj9lXIPs7hjfjzilO83yynhuLHs4BGfKhBHFSf5yQJLGjifrIACvRAYHFWdb2fjEpa4FnKhW9b0AsqdVECBpQlMcwCMo3wfLLMUY+bg4BGfuRmzBONfzXGtUhzVi3uDxedxxsXt8Ub68HIfskFNzGacAV4+xBt+BwPfWhq7z1ZAAQVGEohLWYY/+WS4MD4T2aOAI3H6JAXuJRC3a4ozfAt96s/Ukk1rAIwPgo/PHo37hsUMxIUe8RmiMWMwPotzqUcADxmaYDK1jeKOKaCAAgooMIUCceDna1O4X4vu0jQGwLjpbNzrKz4G6mOLCsx+iHzcYuKJ5blzXQP49XIPseFrACMwbrv55ps58MDmR3mOsFafMnaBU089lbPPHvWjU8e+ehfYEHAs2tUOjkd7xsOxaMdYzMzMsHZtfFAQa8oZvnZsWMWtmLYAGLN3zyo3hf3MHI6PKbd+iev54lrAmA0cHy8VIXBwj7C4zi+WE7eTiTAYP3txuaHt8Czg3QFw27ZtBsCKTTvfqk4++WQuuCDm7/jY2wKOxd4egXuv3/Foz3g4Fu0YiwiAa9ZE9jMAtmNE8lsRoe4uYHtZ1OA+f3FLlwiEcWr4vHKtX9wHMCaBXAxcMrTqM4BXAHF9wDUL3AfQAJgfs7EtwTfWsVGmF+RYpAnHugDHY6ycqYU5Fim+sRUbAGdvhOxj+QIGwOXbjb1yy5YtrF8fZ/N97G0Bx2Jvj8C91+94tGc8HIt2jIUB0ACY7UQDYFbQegUUUEABBSoLGAANgNmWMwBmBa1XQAEFFFCgsoAB0ACYbTkDYFbQegUUUEABBSoLGAANgNmWMwBmBa1XQAEFFFCgsoAB0ACYbTkDYFbQegUUUEABBSoLGAANgNmWWzQAXnUVbN48u5oNG+DKK7v/9cwMrFgBq1dPzz5N0/hMYl8c8+6/bkfpC8e5H+Pc7IW+jvnOnTNcf733AcyGoD7XLxoA162DrVtniVatgu3lDoVd/3ow6F3fD7d/9N50zKfrNTxf7zvO/Rjn5vj3c8xnyj2gvRF0n0NcZt8NgFMUag2DC4fBfv6SMAxM2z9efZ1/7+u8n69tA6A3gs7EP1g0ADZPAW/cCFdcMbvCLn/dPGXQ5f3o+jjU3H7HvPuv21H6xXHuxzg3e6GvY+4pYK8BzMW/EQJgdgXWK6CAAgoooMB4BZwEYgDMdtSiRwCzK7BeAQUUUEABBcYrYAA0AGY7ygCYFbReAQUUUECBygIGQANgtuUMgFlB6xVQQAEFFKgsYAA0AGZbzgCYFbReAQUUUECBygIGQANgtuUMgFlB6xVQQAEFFKgsYAA0AGZbzgCYFbReAQUUUECBygIGQANgtuUMgFlB6xVQQAEFFKgsYAA0AGZbzgCYFbReAQUUUECBygIGQANgtuUMgFlB6xVQQAEFFKgsYAA0AGZbzgCYFbReAQUUUECBygIGQANgtuUMgFlB6xVQQAEFFKgsYAA0AGZbzgCYFbReAQUUUECBygIGQANgtuUMgFlB6xVQQAEFFKgsYAA0AGZbzgCYFbReAQUUUECBygIGQANgtuUMgFlB6xVQQAEFFKgsYAA0AGZbzgCYFbReAQUUUECBygIGQANgtuUMgFlB6xVQQAEFFKgsYAA0AGZbzgCYFbReAQUUUECBygIGQANgtuUMgFlB6xVQQAEFFKgsYAA0AGZbzgCYFbReAQUUUECBygIGQANgtuUMgFlB6xVQQAEFFKgsYAA0AGZbzgCYFbReAQUUUECBygIGQANgtuUMgFlB6xVQQAEFFKgsYAA0AGZbbsEAeNVVsHnz7Co2bIArr+z+1zMzsGIFrF49Pfs0TeMziX1xzLv/ul1qXzjm/RjzPo/zaafN8IIXrImBjv/NZMNAF+tXdHGjW7TNCwbAdetg69bZrV21CrZvn46vB/7TtE/uy8K96ZhP12t4lH53zPsx5n0d5yOOmOH66w2ALcpTndsUA+CUhNpRfiH2+Tl9/SXhmPcjBDnO/RtnA6CngLOJc+RTwBs3whVXzK6uy183Txl0eT+6Pg41t98x7/7rdqn94pj3Y8z7PM6nnz7DCSd4BDAbgvpc7ySQPo+++66AAgoo0EkBJ4F4BDDbuAbArKD1CiiggAIKVBYwABoAsy1nAMwKWq+AAgoooEBlAQOgATDbcgbArKD1CiiggAIKVBYwABoAsy1nAMwKWq+AAgoooEBlAQOgATDbcgbArKD1CiiggAIKVBYwABoAsy1nAMwKWq+AAgoooEBlAQOgATDbcgbArKD1CiiggAIKVBYwABoAsy1nAMwKWq+AAgoooEBlAQOgATDbcgbArKD1CiiggAIKVBYwABoAsy1nAMwKWq+AAgoooEBlAQOgATDbcgbArKD1CiiggAIKVBYwABoAsy1nAMwKWq+AAgoooEBlAQOgATDbcgbArKD1CiiggAIKVBYwABoAsy1nAMwKWq+AAgoooEBlAQOgATDbcgbArKD1CiiggAIKVBYwABoAsy1nAMwKWq+AAgoooEBlAQOgATDbcgbArKD1CiiggAIKVBYwABoAsy1nAMwKWq+AAgoooEBlAQOgATDbcgbArKD1CiiggAIKVBYwABoAsy1nAMwKWq+AAgoooEBlAQOgATDbcgbArKD1CiiggAIKVBYwABoAsy1nAMwKWq+AAgoooEBlAQOgATDbcgbArKD1CiiggAIKVBYwABoAsy1nAMwKWq+AAgoooEBlAQOgATDbcgbArKD1CiiggAIKVBYwABoAsy1nAMwKWq+AAgoooEBlAQOgATDbcgbArKD1CiiggAIKVBYwABoAsy1nAMwKWq+AAgoooEBlAQPg9AXANwPHAg8Bvg18EvgN4KuN3noU8HbgccA3gUuAM4d6L75/GRAB7/PAJmDrHP1pAKz8onV1CiiggAIKZAUMgNMXAN8EfAC4FtgPeAfwCOAxpVkOAL4IXAqcBRwGXA28FXhbec7rgVcBxwBfAk4HTizPvXOo6QyA2Veh9QoooIACClQWMABOXwAcbqFHA18ADgK2AS8BzgUOBu4pT34N8Grg0PL9jcAFwIXl+32BW4CTgfcaACu/Sl2dAgoooIACYxYwAE5/AIzTv68EHlZ6J4Ldw8vRvUE7HQV8GlgD7FNOC8fffa7Rb1vKUcVTDIBjfhW6OAUUUEABBSoLGACnOwA+A/ggsAH4aOmtdwH7Ay9s9NoR5fq+tSUA3lRC4g2N51wOzAAvNwBWfpW6OgUUUEABBcYsYACc3gB4HPCecsr3Q42+8QjgmF9ELk4BBRRQQIGuCRgApzMAvqhcv/d84GNDTRmTOc4bugbwtWXSx0LXAH4deN181wBu2rSJlStX7l7V+vXrd//nQwEFFFBAAQXaI7Blyxbiv3js2LGDiy66KL6My7/iDF/vHiumbI9j9m7M7n0W8Jk59i1mAcep3ZgFHDOGI/R9BDi/MQs4rvOL5cTtZGJCyGnAi4HDAWcBT1nDuDsKKKCAAv0T8Ajg9B0BjJm9dwHbSztHwN1VJn0MAuE64OJyH8CYGRy3innjUPufAbwCWA1c430A+/fm4B4roIACCkyvgAFw+gJg7W71PoC1xV2fAgoooIACSQEDoAEw2UK7Pylk27Zt2zjwwPjShwIKKKCAAgq0XcAAaADM9qgBMCtovQIKKKCAApUFDIAGwGzLGQCzgtYroIACCihQWcAAaADMtpwBMCtovQIKKKCAApUFDIAGwGzLGQCzgtYroIACCihQWcAAaADMtpwBMCtovQIKKKCAApUFDIAGwGzLGQCzgtYroIACCihQWcAAaADMtpwBMCtovQIKKKCAApUFDIAGwGzLGQCzgtYroIACCihQWcAAaADMtpwBMCtovQIKKKCAApUFDIAGwGzLGQCzgtYroIACCihQWcAAaADMtpwBMCtovQIKKKCAApUFDIAGwGzLGQCzgtYroIACCihQWcAAaADMtpwBMCtovQIKKKCAApUFDIAGwGzLGQCzgtYroIACCihQWcAAaADMtpwBMCtovQIKKKCAApUFDIAGwGzLGQCzgtYroIACCihQWcAAaADMtpwBMCtovQIKKKCAApUFDIAGwGzLGQCzgtYroIACCihQWcAAaADMtpwBMCtovQIKKKCAApUFDIAGwGzLGQCzgtYroIACCihQWcAAaADMtpwBMCtovQIKKKCAApUFDIAGwGzLGQCzgtYroIACCihQWcAAaADMtpwBMCtovQIKKKCAApUFDIAGwGzLGQCzgtYroIACCihQWcAAaADMtpwBMCtovQIKKKCAApUFDIAGwGzLGQCzgtYroIACCihQWcAAaADMtpwBMCtovQIKKKCAApUFDIAGwGzLGQCzgtYroIACCihQWcAAaADMtpwBMCtovQIKKKCAApUFDIAGwGzLGQCzgtYroIACCihQWcAAaADMtpwBMCtovQIKKKCAApUFDIAGwGzLGQCzgtYroIACCihQWcAAaADMtpwBMCtovQIKKKCAApUFDIAGwGzLGQCzgtYroIACCihQWcAAaADMtpwBMCtovQIKKKCAApUFDIAGwGzLGQCzgtYroIACCihQWcAAaADMtpwBMCtovQIKKKCAApUFDIAGwGzLGQCzgtYroIACCihQWcAAaADMtpwBMCtovQIKKKCAApUFDIAGwGzLGQCzgtYroIACCihQWcAAaADMtpwBMCtovQIKKKCAApUFDIAGwGzLGQCzgtYroIACCihQWcAAaADMtpwBMCtovQIKKKCAApUFDIAGwGzLGQCzgtYroIACCihQWcAAaADMtpwBMCtovQIKKKCAApUFDIAGwGzLGQCzgtYroIACCihQWcAAaADMtpwBMCtovQIKKKCAApUFDIAGwGzLGQCzgtYroIACCihQWcAAaADMtpwBMCtovQIKKKCAApUFDIAGwGzLGQCzgtYroIACCihQWcAAaADMtpwBMCtovQIKKKCAApUFDIAGwGzLGQCzgtYroIACCihQWcAAaADMtpwBMCtovQIKKKCAApUFDIAGwGzLGQCzgtYroIACCihQWcAAaADMtpwBMCtovQIKKKCAApUFDIAGwGzLGQCzgtYroIACCihQWcAAaADMtpwBMCtovQIKKKCAApUFDIAGwGzLGQCzgtYroIACCihQWcAAaADMtpwBMCtovQIKKKCAApUFDIAGwGzLGQCzgtYroIACCihQWcAAaADMtpwBMCtovQIKKKCAApUFDIAGwGzLGQCzgtYroIACCihQWcAAaADMtpwBMCtovQIKKKCAApUFDIAGwGzLGQCzgtYroIACCihQWcAAaADMtpwBMCtovQIKKKCAApUFDIAGwGzLGQCzgtYroIACCihQWcAAOH0B8AXAJuDRwAHAfYF7Gn0VX38HuJvZfd8FHAVsbTznTOBlQIS7z5flNX/ebFMDYOUXratTQAEFFFAgK2AAnL4A+EzgIGA/4F3zBMCnA5+Yp3leD7wKOAb4EnA6cCJwGHDnHDUGwOyr0HoFFFBAAQUqCxgApy8ADlroKcDH5wmAzyg/m6vdbgQuAC4sP9wXuAU4GXivAbDyK9TVKaCAAgooMAEBA2A/A+A3SjD8CvDOcqQw2iuO5n2znBL+XKPftgDXAqcYACfwKnSRCiiggAIKVBYwAPYvAD4V+CywE4jTxXFU77eA3wceBNwEPBy4odGLlwMzwMsNgJVfoa5OAQUUUECBCQgYAPsXAIfb6DTgZ4AneQRwAq8wF6mAAgoooEALBQyABsAIgOuBJ5b+nOsawK8Dr1voGsBNmzaxcuXK3YtYv3797v98KKCAAgoooEB7BLZs2UL8F48dO3Zw0UUXxZdrylm+9mxopS2JW6FM02Ofcn1fTAK5GlhdTvfuAH683PolrueL28HEbOD3ARECd3dBuc4vZgEfC0QYjJ+9GDjcWcDT1CbuiwIKKKBAnwU8Ajh9RwBfAlxW7u8XvT24119c+xeTPM4r1/rFfQBjEsjFwCVDL4IzgFeU8HiN9wHs81uE+66AAgooMI0CBsDpC4C1+9T7ANYWd30KKKCAAgokBQyABsBkC+0+qrht27ZtHHhgfOlDAQUUUEABBdouYAA0AGZ71ACYFbReAQUUUECBygIGQANgtuUMgFlB6xVQQAEFFKgsYAA0AGZbzgCYFbReAQUUUECBygIGwL0fAL/cmLG70PA/rHJvjLq6BQPgVVfB5s2zizrnnNk/B99v2ABXXjn7d137emYGVqyA1au7t+1d9G7DNjvm3XytLrV3HOd+jHOzL/o65jt3znD99XELQO8DOGrgGffzfrmxwLXAScC7gQiGDwV+sdyq5cxxr3hMy1swAK5bB1u3zq7pyCNn/xx8v2oVbN8++3dd/Hrg18Vtd5uX13eOeTdfq0vtd8e5H+Pc7It+jnl8wqsBcExZKL2YjwFvAP6qsaQnAL8DPCO99MksYN4AGEf/jj9+T8hbuxZuvXXP93EEbdeu2Y3q4tcDzi5uu9u8vL5zzLv5Wl1qv8/1VrnUZfj85b3G9pZbP1/bBsA2fRJIjMYPlE/uGPTjvsDt5SbOk4lwuaXOGwCbR//iX1oPfCDcdFO33hTmezPq55tFP375O+aO8/Bb4t4KJa63zu+Lvr6fr1w5w44dHgHMRaDxVf8D8DbgDxuLfGn5HN5Hjm81Y13SSAHQ07/LO+W41FNXPn9yzv08RdS/U4GeFuzfmPf1tX3EEV4D2KYjgOuBPwXis3rjGsAfASL4PQ/4s7HGtvEtbMFTwCedNHvaNyZLxL9m42Lbwde3x3FNYL/94M4793wd1wW2+Tk7d86eut5nHzjggNlT2m3e3i4at22bHfPuvD6zvbNjB3zrW7P7G6/xeN+K9yhf57MWXX3fnq8vVq6EO+6YfU+/3/3grrv60esPeACcf/4MJ5zgEcDxxaH8kmLixwvL5/V+FXhfCYP5JU9mCSNPApnrX1nzbdJc/woffm6bntOmbQknt2fhZh+HzziW4VjtuSa4C+8FjvnkX1dte01M65jHWbnPfnaGNWsMgJOJRv1YqgHQwPXd2dxd+CU+rl8w0/qLwf2a/41bGwPgXAJd7AsD4OxI7u1TwHGbl/eUporr/eZ7XNrSPDnSfQCb91nauBGuuGL2dPDgIufmn3E6tc3P+c//nD1FFNsZp4bavr1dNG7bNjvms6fIBve+bPPrM9s7l12257KVeH33Zb+7/J7smC+9T889F578ZI8A7u0AeB2wroS7uO5vrkfcLKWTN4JuaWh1sxRQQAEFFOi1gJ8EsvePAHa9Af0ouK6PoNuvgAIKKNA7AQOgATDb9AbArKD1CiiggAIKVBYwABoAsy1nAMwKWq+AAgoooEBlAQOgATDbcgbArKD1CiiggAIKVBYwABoAsy1nAMwKWq+AAgoooEBlAQOgATDbcgbArKD1CiiggAIKVBYwALYvAD4F+Clg9VAvnFa5N0ZdnQFwVCmfp4ACCiigQEsEDIDtCoBnAZuB/wvc0eiRuA/g01rSM8ObYQBs6cC4WQoooIACCswnYABsVwD8OvA84K871LIGwA4NlpuqgAIKKKBACBgA2xUAbwPuD8QRv648DIBdGSm3UwEFFFBAgSJgAGxXAHw78CngTzrUoQbADg2Wm6qAAgoooIBHAGd7YG9/FnCzE98HPBf4DHDLUIue2NKWNQC2dGDcLAUUUEABBeYT8AhguwLgZQu06i+1tI0NgC0dGDdLAQUUUEABA+D8PdCmI4Bd7FQDYBdHzW1WQAEFFOi1gEcA23UEcHBK+ieBBwNfAa5p+aQQA2Cv30LceQUUUECBLgoYANsVAH8Y+BDwWOB24AeALwDPmeOawLb0mwGwLSPhdiiggAIKKDCigAGwXQHw8jIp5VeB/wAOAi4C9gFeMOKY1n6aAbC2uOtTQAEFFFAgKWAAbFcAjBtBHx73Z2yM6xrgBuC/JMd6UuUGwEnJulwFFFBAAQUmJGAAbFcA/AbwY8C3G+Mdnwn8z8ADJ9QD2cUaALOC1sK5cOcAACAASURBVCuggAIKKFBZwADYrgD4fuAu4KRyFDCO/l0IrAJOqNwbo67OADiqlM9TQAEFFFCgJQIGwHYFwIOBjwCPLNcAxiSQ64DjgK+1pGeGN2PeAHjVVbB58+zTN2yAK6+crq9nZmDFCli9ejr3b1rHban75ThP1+t2lPF3zB3zaft9Ndz355wDRx89w5o1cZyJ+F/z0rOWxo3xb1bb7gMYEz5+ClgL3Az8DXDP+Hd7bEucMwBG+Dv+eNi+fXY997kP3H337Nf77gs7d07H1wPFadon9+V7e9Nxnq7X7Sg97pg75tP2+6rZ9w9+MFx7rQGwLQFwZbnly08A3xlbPJv8guYMgOvWwdate1YeR8p27Zr8xrgGBRRQQAEFFFhYYNUquPVWA2BbAmCMVpzmfSiwo0PNu2gAjEaL06S33Ta7V/vsA/eUY5pd/3owTl3fD7d/4d50nKfrdTtKvzvmjvm0/b5q9v3atXDddQbANgXAs4A4afqmrgfA5vV/5547e/RvcD3gxo1wxRWze9jlr5vXCXV5P7o+DpPefse5+6/VpfaIY+6YT8PvqIX6Pn4vP/nJBsA2BcC/BB4P/Gv5GLjmtX9HtzQUOgu4pQPjZimggAIKKDCfgLOA2zUL+PQFWvXMlraxAbClA+NmKaCAAgooYACcvwfadASwi51qAOziqLnNCiiggAK9FvAIYLuOAMYtYOZ7tPVWMAbAXr+FuPMKKKCAAl0UMAC2KwBGyJvvZin7trTBDIAtHRg3SwEFFFBAAU8Bd+MU8FOGNvMQ4NeBS4B3trSNDYAtHRg3SwEFFFBAAQNgNwLgXFv5MODdwJNa2sYGwJYOjJulgAIKKKCAAbC7ATBO/d4ORNBq48MA2MZRcZsUUEABBRRYQMBrANt1DWAc7Ws+9gdeBjwVeFRLO9kA2NKBcbMUUEABBRTwCGA3jgAOTwKJW9TcCJwIfLalbWwAbOnAuFkKKKCAAgoYALsRAB8ytJnfAv6j5e1rAGz5ALl5CiiggAIKDAt4Crhdp4C72KEGwC6OmtusgAIKKNBrAQNguwJg3Aj6N4GXAg8A1gDrgQeXW8G0sVkNgG0cFbdJAQUUUECBBQQMgO0KgGcBzwLeAlwMfD/wo8CfAI9taScbAFs6MG6WAgoooIAC8wkYANsVAL8MHA3cXK79OwiIo4K3AfF1Gx8GwDaOitukgAIKKKCARwAX7IGYaduWRwS9+5ePg4vJHxH67gt8rZwSbst2NrfDANjGUXGbFFBAAQUUMAB2JgB+opz6jVO+gwB4PPAr5VrANjazAbCNo+I2KaCAAgooYADsTAB8PPAx4EPA84D/DUQAfCbwty3tZANgSwfGzVJAAQUUUGA+Aa8BbNc1gDFOhwMnAYcC3yhHBK9pcQsbAFs8OG6aAgoooIACcwkYANsXALvWqQbAro2Y26uAAgoo0HsBA2D7AmB8HvDjgNVD3XlpS7vVANjSgXGzFFBAAQUU8BTw/D3QplnArwQuLBNA7mhs8i4ggmEbHwbANo6K26SAAgoooMACAh4BbNcRwK8ArwOu7FDXGgA7NFhuqgIKKKCAAiFgAGxXAPxm+fSPLnWnAbBLo+W2KqCAAgooYADc3QNtOgX8XuAPgE92qDsNgB0aLDdVAQUUUEABjwDO9kCbAuDvAi8BrgBuGWrR01rasgbAlg6Mm6WAAgoooMB8Ap4CblcAjE8CmesRk0Ce1tI2NgC2dGDcLAUUUEABBQyA8/dAm44AdrFTvycAXnUVbN48uyvnnDP75+D7DRvgyjLFpctfz8zAihWwejV0eT9ibNz+2R5dzMExH81pFMs2P8dx7sc4N3uwr2O+c+cM11+/JijifzNdDCDZbZ62APgCYBPwaOAA4L7APQ2kRwFvL/cajEknlwBnDiHG9y8DItx9vixv6zzQ9wqAEf6OPx62b5999qpVs38Ovo/QtCuOZ8a5945/PfDo+n64/Qv341x9r9l0vIbnG0df291/f17qa7SfYx6ZzwCYDZFtqo/PDT4I2A9411AAjED4RSBuKn0WcBhwNfBW4G1lJ14PvAo4BvgScDpwYnnunXPs6L0C4Lp1sHW+qNgmJbdFAQUUUECBXgsYAKftCOCgnZ8CfHwoAMYEk3OBgxtHBV8DvLp89nDU3ghcUG5IHd/vWyaknAzELOXhx+4A+P73b+OEEw7EANjrdxN3XgEFFFCgAwJxhPSQQ2b46lc9AtiB4VryJs4VACPYPbwc3Rss8Cjg0+U48D5AnBaOv/tcY41bgGuBU+YLgEccsY1//McDiVPAz3se3H33krfXAgUUUEABBRSoJPCgBxkA+3QEME4J7w+8sNFfRwBx0nYtEAHwphISb2g85/JygejLFwuA8fP73W/PNX+V+tjVKKCAAgoooMASBFaunGHHDo8ALoGsM0+tegTw2GM3cdhhK3fjvOc967nttvW7v95nH7inTEGZpq8HXTBN++S+zI7qfA6O+fS+nh3z/o1tc8z79NpesWILu3bFST3Yf/8d3HHHRfGls4A7E+1G29C5AmBM5jhv6BrA15ZJH4eWxc51DeDXy2cUz3sN4LZt2zjwwLgckN2ngQe3fdm4Ea6I21oD0/J185YB07JP0zQ+k9gXx3y6XsOj9Ihj3o8x7/M4n376DCec4BHA0SJVN54Vp3Hj1i8RAGOG72pgJ7CjnP6NU7sxC/hNZeLHR4DzG7OA4zq/mAV8bJkQEp9A8mLgcGDRWcDdIHIrFVBAAQUU6LeAnwTSrk8CGUc3xkzfy4Byt73dH3UXXz8V+BSwDri43AdwG/AO4I1DKz4DeEUJj9cs5T6A49gBl6GAAgoooIACkxUwAE5fAJxsx3zv0v0ouNrirk8BBRRQQIGkgAHQAJhsod2fFrKteQ1gdoHWK6CAAgoooMBkBQyABsBshxkAs4LWK6CAAgooUFnAAGgAzLacATAraL0CCiiggAKVBQyABsBsyxkAs4LWK6CAAgooUFnAAGgAzLacATAraL0CCiiggAKVBQyABsBsyxkAs4LWK6CAAgooUFnAAGgAzLacATAraL0CCiiggAKVBQyABsBsyxkAs4LWK6CAAgooUFnAAGgAzLacATAraL0CCiiggAKVBQyABsBsyxkAs4LWK6CAAgooUFnAAGgAzLacATAraL0CCiiggAKVBQyABsBsyxkAs4LWK6CAAgooUFnAAGgAzLacATAraL0CCiiggAKVBQyABsBsyxkAs4LWK6CAAgooUFnAAGgAzLacATAraL0CCiiggAKVBQyABsBsyxkAs4LWK6CAAgooUFnAAGgAzLacATAraL0CCiiggAKVBQyABsBsyxkAs4LWK6CAAgooUFnAAGgAzLacATAraL0CCiiggAKVBQyABsBsyxkAs4LWK6CAAgooUFnAAGgAzLacATAraL0CCiiggAKVBQyABsBsyxkAs4LWK6CAAgooUFnAAGgAzLacATAraL0CCiiggAKVBQyABsBsyxkAs4LWK6CAAgooUFnAAGgAzLacATAraL0CCiiggAKVBQyABsBsyxkAs4LWK6CAAgooUFnAAGgAzLacATAraL0CCiiggAKVBQyABsBsyxkAs4LWK6CAAgooUFnAAGgAzLacATAraL0CCiiggAKVBQyABsBsyxkAs4LWK6CAAgooUFnAAGgAzLacATAraL0CCiiggAKVBQyABsBsyxkAs4LWK6CAAgooUFnAAGgAzLacATAraL0CCiiggAKVBQyABsBsyxkAs4LWK6CAAgooUFnAAGgAzLacATAraL0CCiiggAKVBQyABsBsyxkAs4LWK6CAAgooUFnAAGgAzLacATAraL0CCiiggAKVBQyABsBsyxkAs4LWK6CAAgooUFnAAGgAzLacATAraL0CCiiggAKVBQyABsBsyxkAs4LWK6CAAgooUFnAAGgAzLacATAraL0CCiiggAKVBQyABsBsyxkAs4LWK6CAAgooUFnAAGgAzLacATAraL0CCiiggAKVBQyABsBsyxkAs4LWK6CAAgooUFnAAGgAzLacATAraL0CCiiggAKVBQyABsBsyxkAs4LWK6CAAgooUFnAAGgAzLacATAraL0CCiiggAKVBQyABsBsyxkAs4LWK6CAAgooUFnAAGgAzLacATAraL0CCiiggAKVBQyABsBsyxkAs4LWK6CAAgooUFnAAGgAzLacATAraL0CCiiggAKVBQyABsBsyxkAs4LWK6CAAgooUFnAAGgAzLacATAraL0CCiiggAKVBQyABsBsyxkAs4LWK6CAAgooUFnAAGgAzLacATAraL0CCiiggAKVBQyABsBsyxkAs4LWK6CAAgooUFnAAGgAzLacATAraL0CCiiggAKVBQyABsBsyxkAs4LWK6CAAgooUFnAAGgAzLacATAraL0CCiiggAKVBQyABsBsyxkAs4LWK6CAAgooUFnAAGgAzLacATAraL0CCiiggAKVBQyABsBsyxkAs4LWK6CAAgooUFnAAGgAzLacATAraL0CCiiggAKVBQyABsBsyxkAs4LWK6CAAgooUFnAAGgAzLacATAraL0CCiiggAKVBQyABsBsyxkAs4LWK6CAAgooUFnAAGgAzLacATAraL0CCiiggAKVBQyA/QyApwNvAO5kdv93AR8GXlT671HA24HHAd8ELgHOnKc3DYCVX7SuTgEFFFBAgayAAbC/AfDpwNFzNNABwBeBS4GzgMOAq4G3Am+b4/kGwOyr0HoFFFBAAQUqCxgADYDDLfcS4FzgYOCe8sPXAK8GDjUAVn6FujoFFFBAAQUmIGAA7G8APKWcAo7TwJ8Ffhv4F+AC4OHAMY1+Owr4NLAG+PZQH3oEcAIvTBepgAIKKKDAJAUMgP0MgI8AvgXcDPww8Bbg8cCjgd8D9gde2Gi8I4CtwFrgFgPgJF+SLlsBBRRQQIHJCxgA+xkAhztrJbANeBbwc8s5Arhp0yZWrozFwPr163f/50MBBRRQQAEF2iOwZcsW4r947Nixg4suuii+jLN7M+3ZynpbErNg+/4YBMBnlyOC5w1dA/ha4FVeA9j3NnH/FVBAAQWmRcAjgP08Avh84OPAvwMPLKeAnwQ8stwW5oYyC/hNJfR9BDjfWcDT8rJ3PxRQQAEF+i5gAOxnAPxT4AnlWr/bgU+V+wLeWF4Q64CLy30A49TwO4A3zvNicRJI399F3H8FFFBAgc4JGAD7GQDH2agGwHFquiwFFFBAAQUqCBgADYDZNjMAZgWtV0ABBRRQoLKAAdAAmG05A2BW0HoFFFBAAQUqCxgADYDZljMAZgWtV0ABBRRQoLKAAdAAmG05A2BW0HoFFFBAAQUqCxgADYDZljMAZgWtV0ABBRRQoLKAAdAAmG05A2BW0HoFFFBAAQUqCxgADYDZljMAZgWtV0ABBRRQoLKAAdAAmG05A2BW0HoFFFBAAQUqCxgADYDZljMAZgWtV0ABBRRQoLKAAdAAmG05A2BW0HoFFFBAAQUqCxgADYDZljMAZgWtV0ABBRRQoLKAAdAAmG05A2BW0HoFFFBAAQUqCxgADYDZljMAZgWtV0ABBRRQoLKAAdAAmG05A2BW0HoFFFBAAQUqCxgADYDZljMAZgWtV0ABBRRQoLKAAdAAmG05A2BW0HoFFFBAAQUqCxgADYDZljMAZgWtV0ABBRRQoLKAAdAAmG05A2BW0HoFFFBAAQUqCxgADYDZljMAZgWtV0ABBRRQoLKAAdAAmG05A2BW0HoFFFBAAQUqCxgADYDZljMAZgWtV0ABBRRQoLKAAdAAmG05A2BW0HoFFFBAAQUqCxgADYDZljMAZgWtV0ABBRRQoLKAAdAAmG05A2BW0HoFFFBAAQUqCxgADYDZljMAZgWtV0ABBRRQoLKAAdAAmG05A2BW0HoFFFBAAQUqCxgADYDZljMAZgWtV0ABBRRQoLKAAdAAmG05A2BW0HoFFFBAAQUqCxgADYDZljMAZgWtV0ABBRRQoLKAAdAAmG05A2BW0HoFFFBAAQUqCxgADYDZljMAZgWtV0ABBRRQoLKAAdAAmG05A2BW0HoFFFBAAQUqCxgADYDZllswAF51FWzePLuKc86Z/bP5/XHHQdeeM7wP07BPXRyHmtvsmE/Xa3iU3nHMHfNpf2/fuXOG669fE7sZ/5vJhoEu1q/o4ka3aJsXDIDr1sHWrbNbu2oVPOABcPPNe77/wAdmA2GXnjO8D9OwT10ch5rb7JhP12t4lN5xzB3z6X9vj8xnAGxRnurcpowcAAchcPv2Pft45JGzXw8CYBeeE0G2uQ/D33dxn9zmhXvQMd/zmu1Lvzvmjvn097oB0COAucy56Cng44/fE5ge/GD413/d830Ejzg13KXnDO/DNOxTF8eh5jY75tP1Gh6ldxxzx3ya39tXrwZPAXsNYC7+waKTQJrX+J17Luzatec6wPj+2GPvfR1gF54zvA/TsE9dHIea2+yYz75VdOH1Oa6+cMwd82l9b4/XiJNADIATD4DZFVivgAIKKKCAAuMVMAAaALMdtegRwOwKrFdAAQUUUECB8QoYAA2A2Y6aNwAO395llFsvdPE5AThtt7rp4jhMaps3bIArr5x9mcT1qpNaz7TcJmkafBzz/vV7H8fcawANgBMLgM1bwMRF19ddB8N/FysfzADu6nOG92Ea9mlax2o5+9WcCdjVHl3Ofg/PDO/Tvjvms78WHPN7/36avvd2ZwE7CzgXAec9AmgA7Ncb6LQGBsOAYcAAPe1BaPY+tYPbe/Un+BoADYATCoDDs3/HNTOvbctxpuBsA03r7NCNG+GKK/bsY9v6z+25910FxjFr0zHvX7/3ccw9Bewp4Fz8G+E2MNkVWK+AAgoooIAC4xVwEogBMNtRzgLOClqvgAIKKKBAZQEDoAEw23IGwKyg9QoooIACClQWMAAaALMtZwDMClqvgAIKKKBAZQEDoAEw23IGwKyg9QoooIACClQWMAAaALMtZwDMClqvgAIKKKBAZQEDoAEw23IGwKyg9QoooIACClQWMAAaALMtZwDMClqvgAIKKKBAZQEDoAEw23IGwKyg9QoooIACClQWMAAaALMtZwDMClqvgAIKKKBAZQEDoAEw23IGwKyg9QoooIACClQWMAAaALMtZwDMClqvgAIKKKBAZQEDoAEw23IGwKyg9QoooIACClQWMAAaALMtZwDMClqvgAIKKKBAZQEDoAEw23IGwKyg9QoooIACClQWMAAaALMtZwDMClqvgAIKKKBAZQEDoAEw23IGwKyg9QoooIACClQWMAAaALMtZwDMClqvgAIKKKBAZQEDoAEw23IGwKyg9QoooIACClQWMAAaALMtZwDMClqvgAIKKKBAZQEDoAEw23IGwKyg9QoooIACClQWMAAaALMtZwDMClqvgAIKKKBAZQEDoAEw23IGwKyg9QoooIACClQWMAAaALMtZwDMClqvgAIKKKBAZQEDoAFwoZY7E3gZECHv88AmYOtQgQGw8ovW1SmggAIKKJAVMAAaAOfrodcDrwKOAb4EnA6cCBwG3NkoMgBmX4VjrN+yZQvr168f4xJd1HIFHIvlyk2mzvGYjOtylupYLEdt/DUGQAPgfF11I3ABcGF5wr7ALcDJwHsNgON/MY5jiSeffDIXXBDD5mNvCzgWe3sE7r1+x6M94+FYtGMsDIAGwLk6MY7qfRM4Cvhc4wlbgGuBUwyA7XgBD2+Fb6ztGRfHoj1jEVvieLRnPByLdoyFAdAAOFcnPgi4CXg4cEPjCZcDM8DLhwPgzTffzIEHRm70sTcFTj31VM4+++y9uQmuuwg4Fu1qBcejPePhWLRjLCIArl27NjZmTfnd3o4Nq7gVKyquqyurWsoRwEOAr3Zlx9xOBRRQQAEFFLiXQBz0+VofTQyAc4/6XNcAfh143dA1gOF3MPCtPjaP+6yAAgoooECHBVaX6/t3dXgflr3pBsC56eI6v5gFfCwQYfA04MXA4UOzgJcNb6ECCiiggAIKKLC3BAyA88ufAbwCiH8hXDPPfQD31ri5XgUUUEABBRRQYNkCBsBl01mogAIKKKCAAgp0U8AAmBu3UT4tJLcGq19Qjr4+GjgAuC9wT4PlUcDbgceV2/dcAsS4NB+O03j66M3lsoiHAN8GPgn8xtBEKMdjPNaLLSUuS4mb0/8QsLN8StFZwMd8bSxGN/GffxB4DvAM4ONlbb4uJs6+ewXxoQ1vKJdqRb6Ja/s+DLzIcfjeATAALr8pR/20kOWvwcoQeCZwELAf8K6hABiB8IvApUD88otParkaeCvwtsLnOI2vj94EfKDcDzPG4x3AI4DHlFU4HuOzXmxJhwK3AtuA+wCvAeIeSD8I3FH+seRrYzHF8f88QvkvlPeteO+KAOjrYvzO8y0xAuDTgaPneILjMIRiAFx+Y476aSHLX4OVTYGnlDfT5hHAlwDnlpnYg6OC8Yvw1UD8goyH4zS5Poqjsl8oAT2CiOMxOeuFlrwK+FXg10rf3+VY7JWBiNuJfBp4UrmX7OAIoK+LesOxUAB0HAyAY+nEpdwrcCwrdCHMFQDjc9/iht3xmc2DR3yCS7wJx80991nCp7pIvHSBOP37SuBhpdTxWLphpuLnym2pote/Uk7P/z/HIkOaqo1Pi/pj4A/LZSqDAOjrIsW6pOIIgHEXjzvLf58Ffhv4l/Lxrv6+aHB6BHBJvfXdJy/l00KWtwarhgXmCoBxSnh/4IWNJx9RroeKW7xHABz1U10UX5pA/HKLa502AB8tpY7H0gzH9ezvB94IPAs4spwCdizGpTvack4q1/2tL0+PMxKDAOhYjGY4jmfFJSlxX96bgR8G3gI8HoizFb/n74t7ExsAl9dyHgFcnlumyiOAGb3x1h4HvKecZvxQY9Ee6Riv81KWFu/lg9PwEcwdi6Xo5Z4bR8DjrEMEjQge8WgGQMci55upXlleF/GPozhi7hHAhqYBcPmtNeqnhSx/DVY2BeYKgHHB9XlD1wC+ttzEe6FrAOf6VBe1RxOI2XQXAs8fmnEa1Y7HaIaTeFZMBIkA+NxyRNaxmITy3MuMa8t+v3ye7OB3akzGifF4PxCnIeNIVByRGlyr7PtUnfEZBMBnF39/XxgAx9J5flrIWBgXXUicxo2JHxEAY4Zv3Jg7bnuxoxzOv6HMAo4ZqhH6PgKc35gF7DgtSjzyE+LTcWK2dfxr+jNzVMUsO8djZM7UE2Oy0+VlJvD9gej/mP34yHLtk2OR4l1S8f3KRKhmUXxGfNzCKi6PuNvXxZI8M0+Of5jGzOt/Bx5YgndMyonXRYRz358MgJn+uletnxYyNsp5FxT/ur6s3M8pnjS4t9NTgU8B64CLy30A41/ccWuSuB6q+XCcxjNOcfQiZphuL4sbjEVMwhkEQsdjPNaLLSXubfYT5RYjt5d7Msa9Ab/cKHQsFlOc3M/jH6mD28DEWhyLyVk3l/ynwBPKwYF4XcTviLgvYJyxcxyGxsBTwHWa0rUooIACCiiggAKtETAAtmYo3BAFFFBAAQUUUKCOgAGwjrNrUUABBRRQQAEFWiNgAGzNULghCiiggAIKKKBAHQEDYB1n16KAAgoooIACCrRGwADYmqFwQxRQQAEFFFBAgToCBsA6zq5FAQUUUEABBRRojYABsDVD4YYooIACCiiggAJ1BAyAdZxdiwIKKKCAAgoo0BoBA2BrhsINUUCBMQj8VvnYwJ8dw7JchAIKKDC1AgbAqR1ad0yBqRP4VuMjAVcB+5bPvZ3rI+lq7PzLgOcCx9VYmetQQAEFxilgABynpstSQIFaAvF5z08EnlZrhXOs52rg/cAf7cVtcNUKKKDAsgQMgMtis0gBBfaywHwB8HTgGcCTy/Z9ArgWeABwDPBN4PXAF4F3AEcCW4ETgX8qNfsArwN+GTi4/P1vAh9v7PMa4OvAIUB86PyXSxA8qgTTrwGvLEcpLwB+BPhr4BeBW8tyXgX8GnD/ciQzAuVL97Krq1dAgZ4IGAB7MtDupgJTJrBQAHw6cHQjAD4aeBbwGeA1wNnAXwIvLyHucmD/EhCj7Azg2cAJwD8DzwHeCzyyBL14zouB/1bCZnwfATAesZ4IlG8BXlRCXwTJ7cCfA/8A/CrwY8DfAz8B/COwH/BY4NNTNk7ujgIKtFTAANjSgXGzFFBgQYGlBMAIcb9SlnZgOQr488Afl7/bAFwC/GD5Po4SPh/4aGMLIrz9RQmP8dcfLIEujiIOAmAsI8JlPB4F/B0QRwT/pvzdySUUPq4cEYygGCHyz4C4vtGHAgooUE3AAFiN2hUpoMAYBZYSAONo32ll3TFx5C7gp4FPlb9bD3wYWFlOFX8DmAHuKT+P98n7AO8GNpWjdXEaN47ixXMHATC26dLy/Y+W08wPBW4qf/cK4NeBw8r3cbTwJOAJ5bn/A4ijkT4UUECBiQsYACdO7AoUUGACApMKgBEC45q+CIXznY49vly796TGfsUp4KUGwEF5XHO4sYS/w8tp5wmQuUgFFFBgj4AB0G5QQIEuCkwqAIbF+eXUbUzIuB74vnKtXkz6iNPJ7yundeOI3eAxVwCMSSUx+WOuI4BxFDCOEsZRyDtK4PwIcGjjesIujovbrIACHREwAHZkoNxMBRS4l8CoATBm7saRvOYp4B3AU+c5BRwriffFmKEbk0QeBHwH+AJwCvAl4N/KNX5faWzRjcDvLOEU8DrgnWUWcqwvQuKbS7h0qBVQQIGJCxgAJ07sChRQYIoE4qbPMUs4Zu/6UEABBTorYADs7NC54QoosBcE4hYzMZEkZgX7UEABBTorYADs7NC54QoooIACCiigwPIEDIDLc7NKAQUUUEABBRTorIABsLND54YroIACCiiggALLEzAALs/NKgUUUEABBRRQoLMCBsDODp0broACCiiggAIKLE/AALg8N6sUUEABBRRQQIHOChgArvWv5gAAAH1JREFUOzt0brgCCiiggAIKKLA8AQPg8tysUkABBRRQQAEFOitgAOzs0LnhCiiggAIKKKDA8gQMgMtzs0oBBRRQQAEFFOisgAGws0PnhiuggAIKKKCAAssTMAAuz80qBRRQQAEFFFCgswIGwM4OnRuugAIKKKCAAgosT+D/A2wjKXQmn3ugAAAAAElFTkSuQmCC\">"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "for lsm_i in range(number_of_liquids):\n",
    "    if len(times[lsm_i]):\n",
    "        plt.figure()\n",
    "        plt.plot(times[lsm_i],spikes[lsm_i],'.')\n",
    "        plt.ylabel('neuron id')\n",
    "        plt.xlabel('Time/ms')\n",
    "        plt.title('spikes - lsm_i:'+str(lsm_i))\n",
    "        plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 2",
   "language": "python",
   "name": "python2"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 2
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython2",
   "version": "2.7.11"
  },
  "widgets": {
   "state": {},
   "version": "1.1.2"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 0
}