mats.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "%matplotlib notebook"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "import sys\n",
    "sys.path.append(\"/scratch/tc427/lucy/consensus/code\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "import numpy as np\n",
    "import consensus"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "prog = \"mcfold\"\n",
    "cut = 5\n",
    "seq_cl1 = \"GGAAGAGGAAAGAAGTGCAGAGCCAGTATTCGGGTACATCCAGACGTGCCTCTCTTTTGGATTTCTTGTGCTCGCGAGGCAATCA\""
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "def cutmap(seq, cut, prog=\"mcfold\"):\n",
    "    cpred = consensus.Consensus(prog)\n",
    "    # store the energy of the lowest strcture for each cut\n",
    "    cutmap=np.zeros((cut, cut), dtype=float)\n",
    "    # store the 2D structure in here, in dot-bracket format\n",
    "    cutstruct=np.zeros((cut, cut), dtype=\"U\"+str(len(seq)))\n",
    "    for i in range(cut):\n",
    "        for ji in range(cut):\n",
    "            j = len(seq)-ji\n",
    "            #print(seq[i:j])\n",
    "            p = cpred.run(seq[i:j], outdir=None, rmdir=True)\n",
    "            cutmap[i,ji] = p.ix[0,1]\n",
    "            cutstruct[i,ji] = p.ix[0,0]\n",
    "    return cutmap, cutstruct"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 39,
   "metadata": {
    "collapsed": false,
    "scrolled": true
   },
   "outputs": [],
   "source": [
    "cl1_maps = cutmap(seq_cl1, 30, 'RNAfold')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 33,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "def normalise_len(mat):\n",
    "    # create a matrix of the sequence lengths, to be used as a normalisation factor\n",
    "    l = np.zeros(mat.shape)\n",
    "    for i in range(l.shape[0]):\n",
    "        for j in range(l.shape[1]):\n",
    "            l[i,j] = len(mat[i,j])\n",
    "    return l"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "with sns.axes_style(\"white\"):\n",
    "    fig1, axarr = plt.subplots(len(fprogs))\n",
    "    for ax, folds in zip(axarr, folds_cut30):\n",
    "        ax.imshow(folds, cmap='viridis', interpolation='none')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 34,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "def showmaps(Emat, Smat):\n",
    "    l = normalise_len(Smat)\n",
    "    with sns.axes_style(\"white\"):\n",
    "        fig1, (ax1, ax2) = plt.subplots(ncols=2)\n",
    "        ax1.imshow(Emat, cmap='terrain', interpolation='none')\n",
    "        ax2.imshow(Emat/l, cmap='terrain', interpolation='none')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 37,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "seq_cl2 = \"GGAAGAGGAAAGAAGUGCAGAGCGAUACUCGGGUACAUUCGAUUAAGUCCUCAGAAUCGAAUUUCUUGUGCUCGCGAGGCAAUCA\"\n",
    "cl2_maps = cutmap(seq_cl2, 30, 'RNAfold')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 38,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "seq_cl3 = \"GGAAGAGGAAAGAAGTGCAGATTCCCAAACAATAACGCATCACATGAAGTTCCATCCTCGATTTCTTGTGCTCGCGAGGCAATCA\"\n",
    "cl3_maps = cutmap(seq_cl3, 30, 'RNAfold')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 41,
   "metadata": {
    "collapsed": 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,iVBORw0KGgoAAAANSUhEUgAAAoAAAAHgCAYAAAA10dzkAAAgAElEQVR4nO3dedSddX33+8vKoK22HIX21AQNe1+/33rq6alDT9dRa/r0efr0HNc5esTW4QFBAZEwJyRkghAyJyRkgBABmRJEQAFlNkzKDDJkIkCYA4JhEAijoRX3+cPcLhbJ77eTK8ne+7rzeq31WXonue9epeTb913Sy6IAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA+oGyLL8WQlgZY3wohHBcu18/YMCAD4cQ9m42m58ry/LTZtbZNZvNz4UQ9h4wYMCHO3EjUkq3w6xW65XbQQ/YY489/irGuHrgwIEfKorivSGE60II/5r7nBDC3jHGlpl1dyGEvTt0KjbgdpjVd928HfSIsiz3CSGc9Y6P940xnpn7nGaz+bkYY+tvPnNM6+OD59oW7m8y22XYTdkN+Om+yf3TVd/pyo6+aWh2P717bHLXXXl07XbtFSOy+/HZhyV3zrwDsxv1zS9vdAft+X+1YoytZrP5uW1/JTZuS27HlAsPap1x56gNNuu2Edkde9PQ5IbccFh237j24OT+x9UHZvepy/dLbvef7JPdn160d3I7nzs0v9NOTG/OJdnlbsPF9xyT3Z3Lj0/uqnuPzW7RL0cnN+f2EdndtPS45G742cjsThixV3Injf9WdmfM+k5y5516cHY/+cERyV1z2Yjsbv3FmOSuuuSo7C49/8jkxuz75eSG/Fv3bwc9oizL0SGESX0fNxqNfwkhLG7zOZ+OMbY+Pnhu63/7lx/ZFu7jmf3Z5Gez++vlhyX3rw+N6MpmrJmY3fLX5yS3+sEZtdsTD0zP7u7rxyV3w4+Pzm7mIXttdKO++eVWjLFVluWnt/2VSN6ByrfjjDtHtRavnbnBzvvNtOzmPTspuTG/Oi67A58Ym9wXHz46u//j/qHJDVx2aHY73DMkueKmSflddX56Fy7NLncb7n1jbnZr/3NBcqvenJfdja/MSu7CF6dnt+at+cn96pETsjt3+kHJXXzK4dktPm94crdcMTa7ZTcfn9zjK6dn98Izs5N7aMmU7FbcPiG5Ew/bO7kx+3b/dtAjYoxj333EY4xX5z5HAApAASgAt+R2CEABKAAFIF327n9sU5blPmVZfr/N5whAASgAt/MA3JLbIQAFoAAUgHTZxz72sb+OMT5eluVuRVHsGEK4pizLL+c+RwAKQAEoALfkdghAASgABSA9IITw7zHGFTHGVTHGGe1+vQAUgAJQABZF9dshAAWgABSA1JAAFIACUABWIQAFoAAUgNSYABSA2yoAl996fHaXnzm0X+2S7x2eXX8NwBNvH9H64YvTN1gu8OY9O6k19unxyeUC78Anxra+9MjI5P7hgWHZ7b78sOR2vGdIdpUD76rzW8WFy9L7/lv53XVocletnZndQ789KbkbXz0xu1zgnfTc5OzufX1uci8/n98Zxx+Q3PTvfiu7Sd/+WnJzjvxmdt8/7oDkzp99cHa523DjT0Zld8fPjknu1GP2S27CIV+t5e2gRwhAASgABWAVAlAACkABSI0JQAEoAAVgFQJQAApAAUiNCUABKAAFYBUCUAAKQAFIjQlAASgABWAVAlAACkABSI0JQAEoAAVgFQJQAApAAUiNCUABKAAFYBUCUAAKQAFIjQlAAditAMwd2y3Z6eP2z27CPl/tyibv9z83umFf/f9qecT7bsexNw3d7MBrF3m5wGsXebnAaxd5W/Yuv0zgXbisVZz+VnI7z3gxu+L6+cmd9cLU7KoG3oVtIn7kr47LLhelb//u1Ox+eOLByS2cOiG7U0YNTW7mIQdkl/v9e/bk47Kbf/QRybWL1sn7fz251DeOdf7mkR4hAAWgABSAVQhAASgABSA1JgAFoAAUgFUIQAEoAAUgNSYABaAAFIBVCEABKAAFIDUmAAWgABSAVQhAASgABSA1JgAFoAAUgFUIQAEoAAUgNSYABaAAFIBVCEABKAAFIDUmAAWgABSAVQhAASgABSA1JgAF4PYWgNMP+nZyM4bsl933xgxP7ryZUytt1qjhtTzifbfj4J8f7mXObQKvOP2t1vumv5jcLiMfzq649GfJTXzm+OyqBl67yNvrsdHZnf78lOTWvDU/u8vOODK5BaOHZVc18Cbs89XWiYftnVy78Mw907kzJmX345NnJ5eLzmnDDqvl7aBHCEABKAAFYBUCUAAKQAFIjQlAASgABWAVAlAACkABSI0JQAEoAAVgFQJQAApAAUiNCUABKAAFYBUCUAAKQAFIjQlAASgABWAVAlAACkABSI0JQAEoAAVgFQJQAApAAUiN9R3xT3xyUuvTf3+6bcPtPHNtdgKwMwHYLvJyO3XsiOQWn7ew0s6YNaOWR7zvdvy3qw9sfeHhozfYp+4fml3u7/f33H1QdsWNM9K78qL8zr8/vdPezu59M15MbpeRj2Q3YL/bkyu/dEV2xcW3JHfwk8dktyXv8ssF3uBVw7M77pnjk7v39TnZ3fDjo5Pbknf55QLvxMP2bp0+7oDkThk1LLuqd+PUsSNaN192SXJXLzorudNnTK3l7aBHCEABKAAFYBUCUAAKQAFIjQlAASgABWAVAlAACkABSI0JQAEoAAVgFQJQAApAAUiNCUABKAAFYBUCUAAKQAFIjQlAASgABWAVAlAACkABSI0JQAEoAAVgFQJQAApAAUgPCSFcEGNcFUJYEkJYUpbll3O/XgAKQAEoAIui+u0QgAJQAApAekAI4eFBgwbtsqm/vu+If3Dkna2dT1jb0RUXLM3vjmHb1QRgvQPwnp9fW2kXLTyrJ4541dtx0d1jW/e8PmeDLV47M7vzfjMtuRPXTMzu0pdOSG7NW6dkt+z1ucktXjsru4UvTEtu+q8nZjfsyXHJffXR0dl9/sHhyeVeqP2lR0Zu0cucBz94VHJ7rDg8u/2fGJvcZS/NyO7u68cltyUvc253G3544pDk2t2Gyft/PblZh34nu1X33p3ckhtvSO6SRef0xO2gBwwcOPBDIYSXQwhXxBiXl2U5vt3nCMDemQAUgJ24ExuzJbdDAApAASgA6bIQwt/EGH80aNCgXQYOHPj+GOMNIYT9cp8jAHtnAlAAduhUbGBLbocAFIACUADSY5rN5p4hhItzv0YA9s4EoADs1G1oZ3NuhwAUgAJQANJlzWbz70MIX+z7OITwbzHGC3OfIwB7ZwJQAG77K7FxW3I7BKAAFIACkC4ry/KzIYRHdttttw8URbFTCOHaZrP5jTafIwB7ZAJQAHboVGzsDlS+HQJQAApAAUgPKMvyqBDCAzHGh0IIUzbh1wvAHpkAFICduBGZW1DpdghAASgABSA19McAHPXL1s4zX9nq22nm2uSKC5fmd+ewfrah2QnAegfgvb+4rtIuXtQbAbi5+m7Hmb8c3brmlVkbLBd45/1mWmv2s5OSa/d+utOfn5LcsjfmZnfN2pnJLfzNtOxyvxeGPTUuu689Njq5zz94VHa50Gr3vsWqgTf4waNag1Ycntx77jkou9xdWfDs5OwevGdScu3f5Zf+vX/+7IOzu/zsocnlAm/y/l9vzR26T3KTvv217HKRd9NPL0ru3FPm1fJ20CMEoAAUgAKwCgEoAAWgAKTGBKAAFIACsAoBKAAFoACkxgSgABSAArAKASgABaAApMYEoAAUgAKwCgEoAAWgAKTGBKAAFIACsAoBKAAFoACkxgSgABSAArAKASgABaAApMYEoAAUgAKwCgEoAAWgAKTGBKAA7H8BmN+Jh383udlHHJTdBXNOSG57DcBZt43Y7MBrF3n7Pj4mu4nPHJ9cLvDaRV67v9+rBl67yGv3UuUd7x2S3G7LDs2uauANWnF46z13H5Rccf387MJ9RyTXLvCfeXxmcu2+wcsF3hVnD8vupp+OTm7mIXtld+ox+yfX7q7kIu+qhWckd9qMKbW8HfQIASgABaAArEIACkABKACpMQEoAAWgAKxCAApAASgAqTEBKAAFoACsQgAKQAEoAKkxASgABaAArEIACkABKACpMQEoAAWgAKxCAApAASgAqTEBKAAFoACsQgAKQAEoAKkxASgAt7cAPGn4Ycnd8/NrK297DcAxNw2t9C6/XOC1+/vy4CePSa6O7/Lb8Z4h2RU3Tk+vzV2pGnh/iLxT0rv4luxy0Zp7N+Fej41uvfbyycm1e5dfNvAuHZ3d0pvGJzdv2L7Z5d7ld+NPfpzdled8P7ncN52zx46s5e2gRwhAASgABWAVAlAACkABSI0JQAEoAAVgFQJQAApAAUiNCUABKAAFYBUCUAAKQAFIjQlAASgABWAVAlAACkABSI0JQAEoAAVgFQJQAApAAUiNCUABKAAFYBUCUAAKQAFIjQlAASgABWAVAlAACkABSI31HfGdp9/UKs55dutv4Zr0rv5hfl0Ptt4JwCGrx26z5f6X6VkvTM2uvwVg1YjbktU9AKfcclTrzOenbLApzxyf3eGrj0mu3UuVxzx9XHI9+zLnXOD9YkZ+l/4svcXnZPfJlUcmt8uSQ7IrbhuV3pUXZfe+e4ck99kHj8pu3e++l1zblzlnAm/Zzcdn99h905I7e9KB2eVe5pwLvD9E3ozkTh83KrlJhx1Uy9tBjxCAAlAACsAqBKAAFIACkBoTgAJQAArAKgSgABSAApAaE4ACUAAKwCoEoAAUgAKQGhOAAlAACsAqBKAAFIACkBoTgAJQAArAKgSgABSAApAuKcvyz2OMKxqNxkfXf/xfQwhLY4yrYowLiqL4k034GgJQAArA7SwAt+btEIACUAAKQDqoLMvPxBhXhBDWrT/iO8UYV5dl2SyKogghLIwxfncTvs4fAvCC/VvFHcN6a10Ptt4JwHYhtnrdyV3Z05kJwN4MwK19O25aelxrzVvzN9iy1+dkt3jtzOQWvjA1u+m/npBcz77LLxt4i/PLfSPd5n2quW/wLn5xenbznp2U3LAnx2W35yMjkxu8anh2D705L7l27/LLBd7jK6dn9+Kzc5K7aP7h2VUNvHaRN+vQ7yQ3ep9vCMDtUQjhrGaz+Y8xxscbjcZHG43G52OMN/T9fIxx8Ds/ThGAAlAAbl8BuLVvhwAUgAJQANIFMcYnGo3GR5vN5jdijOf2/XhZls0Y46p2ny8ABaAA3L4CsM/Wuh0CUAAKQAFIF/Qd8RDCXu8+4iGEB9p9vgAUgAJw+w7ALb0dAlAACkABSBf0HfH1/9jm+r4fbzQanw8hXNvu8wWgABSA23cAbuntEIACUAAKQLqg74gXRbFTCOHJRqMRi+IPf84nhDC03ecLQAEoALfvACy28HYIQAEoAAUgXdD3B7nX//vBIYQlIYQHQggLi6LYod3nC0ABKAC32wDcKrdDAApAASgAqSEBKAAF4PYZgFtKAApAASgAqbE/BuAPD20Vtx5r23LbMADfWregO/ttendfPy47Adg/AvDCu8e27nx9zga77OUTsjv9+SnJTXzm+Oy+u/qY5L70yMjsPnX/0OR2W3ZodtlvWK89Pb8f/zK977+V3Q6zX0uuOO3t7L7w8NHJ3fP6nOwufWlGcrk4nPfspNbRT41L7iuPjspu5Rtzk2v3Mudc4L307Nzs3v7dqcldd+HR2VUNvHaRN3n/ryc37GtfrOXtoEcIQAEoAAVgFQJQAApAAUiNCUABKAAFYBUCUAAKQAFIjQlAASgABWAVAlAACkABSI0JQAEoAAVgFQJQAApAAUiNCUABKAAFYBUCUAAKQAFIjQlAASgABWAVAlAACkABSI319HsA+9u2JACfn5rdw0vTu/nSMdldfubQyrtw7qHJbavAE4C9oe92nHHnqK3/Lr8nj83uq4+OSq7du/wGLT8suS16l99PF+d3zprkdjzx1ez+Yszjyf3Z5Gez+y8rj0zulRfnZffQ0inJ3Xvj+OweWTE1ufvvmpjdlrzLr2rgvf27U1tPrjs5uXbf0FYNvHaRd/LIbyU37rv/XsvbQY8QgAJQAArAKgSgABSAApAaE4ACUAAKwCoEoAAUgAKQGhOAAlAACsAqBKAAFIACkBoTgAJQAArAKgSgABSAApAaE4ACUAAKwCoEoAAUgAKQGhOAAlAACsAqBKAAFIACkBoTgAJQAArAKgSgABSAApAa+2MAXrh/q7hzWIeXD6JyxRHb1XIvqG0XgA8tmZLcTT8dnV0u8C4748jsLphzSHICcPsIwBm3jWgt/M20DTZjzcTshj01LrmvPTY6u8Grhid34BNjs/u/HxqR3N+uPDK7D9x7cHLFLePyu+Li9H7waHbvPXldcu+f+nx+9w5Jbt3vvpfdmtWzknt0xdTsnvvVicn96pETsttWL3Neve7k7G5YOyu5R5ZPza5q4E3e/+ut+SO/ldy5JxyU3Akj96rl7aBHCMDemQAUgN2+B5tDAApAASgAqTEB2DsTgAKw2/dgcwhAASgABSA1JgB7ZwJQAHb7HmwOASgABaAApMYEYO9MAArAbt+DzSEABaAAFIDUmADsnQlAAdjte7A5BKAAFIACkBoTgL0zASgAu30PNocAFIACUABSY90NwPz+edXwruyLDx/dlX3pkZHJXfbSjOwEoADstL7bcfRNQzc78NpFXrt3+e2x4vDkhqwem13u9+AnVx6Z3a5LD0kuF4cfuPfgVnHTpPQuuzy/RU8lt8O8N7Lb8d4hyT227uTs3nh1fnLPPz07u7d+uyC5114+ObsteZdfNvBemZXd6c9PSe6FZ2ZnVzXw5o/8VusHM4ckd9XCYcmdNmP/Wt4OeoQAFIACUABWIQAFoAAUgNSYABSAAlAAViEABaAAFIDUmAAUgAJQAFYhAAWgABSA1JgAFIACUABWIQAFoAAUgNSYABSAAlAAViEABaAAFIDUmAAUgAJQAFYhAAWgABSAdElZln8eY1zRaDQ+uv7jaTHGJ0IIS0IIS8qyPGQTvoYAFIACcDsLwK15OwSgABSAApAOKsvyMzHGFSGEdX1HPMZ4fYzxE5v5dXo2AP/bQyO6slyI3frqidnlXqra7ujl9tDSqfkJQAG46b/nt+rt2P+Gwyq9zLlq4O2x4vBs1HzxkaOzy71wfbelh2ZX3D4qvZsnZjf91+ld+fLM7O54bU5yF784I7tHf3tScu2C6JX/XJBcLvDe+u2C1rP/cUpyb/zue9lt0cucM//9nPbc5OyOempccu3Cs2rg/WDmkNZVi45K7o6fHZPcBd8/RABuj0IIZzWbzX+MMT7+jiP+Ygjhkhjj8hDCvKIodmr3dQSgABSA21cAbu3bIQAFoAAUgHRBjPGJRqPx0QEDBnw4hHB5s9ncvSiKPwkhnBdjnNDu8wWgABSA21cA9tlat0MACkABKADpgr4j/u4fDyF8MsZ4T7vPF4ACUABu3wH47h/f3NshAAWgABSAdEHfEW82m2VZlt/s+/GyLD8dQriz3ecLQAEoALfvANzS2yEABaAAFIB0Qd8RbzQaHw0hPLf77rt/pCiK94QQzinLcnS7zxeAAlAAbt8BuKW3QwAKQAEoAOmCd/5B7hDCXiGElTHGVTHGM4qi2KHd5wtAASgAt9sA3Cq3QwAKQAEoAKkhASgABeD2GYBbSgAKQAEoAKkxAbh5u+212dnlAvD1tfOzywXgA3dNyu7+X05M7roLR2QnAAVgFX234/+9ZshWf5nzjvcMya64cXpyucBrF3nFbaPyu/r89M6/P7uDnzwmudwLiE9/fkrr4pdmJDf72UnZ3fP6nOQueWlGdr/5zwXJrXlrfnZ3vj4nufvfnJfdlrzMORd47V5O/oWHj05u1ZvzsqsaeFctOqp15+Jjk8t9c3/VJUfV8nbQIwSgABSAArAKASgABaAApMYEoAAUgAKwCgEoAAWgAKTGBKAAFIACsAoBKAAFoACkxgSgABSAArAKASgABaAApMYEoAAUgAKwCgEoAAWgAKTGBKAAFIACsAoBKAAFoACkxno5APd/Ymxyl740I7tciHVrK26bkF0uxLo1ASgAU/pux/955QHV3uWXC7xfzMjv0p8l1+6buDOfn5LcC/9xSnbLXp+bXLt3+eXCZOzT47PL3ZVHVkzL7rmnZyfX7l1+2cB7bXZ2F744PbnLXj4huy15l1/VwPvCw0dn/3697KUZ2VUNvDsXH9t6eOnU5F5+fm5yd9xyTC1vBz1CAApAASgAqxCAAlAACkBqTAAKQAEoAKsQgAJQAApAakwACkABKACrEIACUAAKQGpMAApAASgAqxCAAlAACkBqTAAKQAEoAKsQgAJQAApAakwACkABKACrEIACUAAKQGpMAApAASgAqxCAAlAACkBqrK4B2O6AvPX2qV1Z7lAvv/X47M6edGDt1q3I21YBuPrB+zu+6668vJZHvO927HHpt6q9zDkbeIvzW7gmuXDfEdlNfOb45Ja/MTe7q14+IbnTnp+S3dhfHZfcXo+Nzu6hN+cl1y4uci+Q//UTs7J747VTklv7nwuyW/b6nOSuXzszu231Muc9lh+eXXH7qORm/HpCdrmbdPEph2d3409GJXffHROTu+yCobW8HfQIASgABaAArEIACkABKACpMQEoAAWgAKxCAApAASgAqTEBKAAFoACsQgAKQAEoAKkxASgABaAArEIACkABKACpMQEoAAWgAKxCAApAASgAqTEBKAAFoACsQgAKQAEoAKmx/hqAT741vztbd3Jyd18/Lrtuh1N/2ZYE4GuvrO347vrlnbU84n23488u2rviu/wygXfOmux2mP1qcjveOyS73F056bnJ2R311LjkvvLoqOw+df/Q5D609JDsVr4xN7llNx+f3eMrpyf30nNzs3v7d6cmt3rdydlVDbwtfpdfLvBuG5Vf5v2S+z42OruJ+349uXnD9s3u3OkHJXf1oqOSO33GAbW8HfQIASgA+9sEYGcIQAEoAAUgNSYABWB/mwDsDAEoAAWgAKTGBKAA7G8TgJ0hAAWgABSA1JgAFID9bQKwMwSgABSAApAaE4ACsL9NAHaGABSAAlAA0gVlWQ4PIdwXY1wRQjirKIodms3m34YQ7gghPBBj/PHAgQPfvwlfRwAKwH41Adj29/xWvR0CUAAKQAFIhzQajX+IMS4fNGjQ+4qiKMqyXFSW5bAQwtJGo/H5oiiKGOPEsiyntftaAlAA9rcJwM7eDgEoAAWgAKRDms1mGWMc3PdxjHFECOHsEMJj7/g1u7/z45T+GoC/ffOUbbM38nvm8ZnJdTuMtpdtSQD+/ve/7/juu+++jh3xbXI7Fh29+YHXJvJ2PPHV7P5izGPJFTdOz+5fHxqRXLsXFC98YWpyE54en91+j49J7uinxmVXNfAeXzm99dKzc5PLBV67yLvhlVnZde1lzrnA++ni/E57O7lPrTwyu5mHHJDctAO/kd1px+6fXO4F0ieN/5YA3J41m82/jDE+EUIYFUK4+R0/9d4Qwrp2ny8ABWB/mwDcNFvtdghAASgABSCdNWjQoEEhhJVlWY4py/KzGznib7b7GgJQAPa3CcD2turtEIACUAAKQDonhPDJEMLTZVkeWhR//Mc2j/b9fFmWA0MID7f7OgJQAPa3CcC8rX47BKAAFIACkM4oy3K3GOOaZrO55zt/PMa4vO/P95RlOT6EMHcTvpYAFID9agIw+/t9698OASgABaAApDNCCFNijK+HEJaEEJau/9fJZVl+PIRwewhhZQjhil133fWD7b6WABSA/W0CsMO3QwAKQAEoAKkfASgA+9sEYGcIQAEoAAUgNfbHI75wTKv4+eye2hFPHptcuwBc9+aCbTIB2PvbkgDshpUrV9byiP/xdsz6Was475EN9t7567L7wHFPJfdXBy/Jrrnn1ckV18/PLtx3RHK5wFv4wtTWxGeOTy4XePs9PqY1eNXw5L726OjsqgZeu8hr9y6/qoHX1Xf55QLv1Lez22Xkw8m95+7vZlc18NpF3k0/HZ3cuacMqeXtoEcIQAHY3yYAO0MACkABKACpMQEoAPvbBGBnCEABKAAFIDUmAAVgf5sA7AwBKAAFoACkxgSgAOxvE4CdIQAFoAAUgNSYABSA/W0CsDMEoAAUgAKQGhOAArC/TQB2hgAUgAJQAFJjAlAA9rcJwM4QgAJQAApAaqzviE+5ZXjrzOenbnSXvnRClzYjuXYvIe3FAJw/6tvWAzt17IjkBOCm67sdM24d0Vr4wrQNlronfZvw9PHJ7ff42Oz+x6oRyZUrjshu9+WHJZcLvHaRlwu8wauGZ//rfur+odmtfnBGcmtfmJfd794+NblHf3tSdovXzkxu3rOTshuyemxy/7xqeHb/67LDktvoS8ffuR/dlVy7F4z/5WHLkituG5ld1cC7+JTDWzdfOjq5B++ZlNyVFw2t5e2gRwhAAbi9TQBuHQJQAApAAUiNCUABuL1NAG4dAlAACkABSI0JQAG4vU0Abh0CUAAKQAFIjQlAAbi9TQBuHQJQAApAAUiNCUABuL1NAG4dAlAACkABSI0JQAG4vU0Abh0CUAAKQAFIjfUd8QV3jEy+U+/WV0/cJrvt1dnZtYut3HoxAOcN2ze76y44N7lfP/GYdWACcNP13Y6ptw5vnfXC1A123m+mZTf91xOSy8Wy/VoAAA2mSURBVMXDkNVjs++Qa/cuvwlPj09uS97lt/uyw7Ir7hya3AeXHJzdmtWzkst9Q3Pq2BGtx+5bltyyW2/M7ubLL0nuyoXfz+78OTOS+97Y4dlN++63kpt+0P/M7vTjDkjuolMOz+4Xl4xKbsXtE7KrGng3Xzq6tereycm9/Pzc5O645Zha3g56hAAUgCYAqxCAAlAACkBqTAAKQBOAVQhAASgABSA1JgAFoAnAKgSgABSAApAaE4AC0ARgFQJQAApAAUiNCUABaAKwCgEoAAWgAKTGBKAANAFYhQAUgAJQAFJjAlAAmgCsQgAKQAEoAKmxviP+w7vGpEPttdld2baKuN4NwB8k9/vf/946MAG46fpux6zbRmw88NZMyO7g1ccklwu8Lzx8dOu/rDwyuW36Mudc4N0xLL8rL0ruPTdOzy53c2YM2S+73N/vr770Yna5b5ZW3Xt3dnddtzi563/0w+wuXjAvuTMnHJPdrEMPSG7agd/Iblu9zDkXeO0i79m35id309Ljank76BECUACaAKxCAApAASgAqTEBKABNAFYhAAWgABSA1JgAFIAmAKsQgAJQAApAakwACkATgFUIQAEoAAUgNSYABaAJwCoEoAAUgAKQLijLcngI4b4Y44oQwllFUewYYzwoxvhMCGHJ+k3ehK8jAAWgbUcBuLVvhwAUgAJQANIhjUbjH2KMywcNGvS+oiiKsiwXlWU5LMZ4ZlmWX96cr9V3xM/85ejWNa/M6ql1O/Z6KQDpvzoZgNvidhxww2Gto54at8G+8uio7D51/9Dk/pelh2RX3Hpsctv0XX4VA6+48qJWce6T6V18S3Zr3pqfXLuoee3ll5Jr917Mh5bcndzd11+TXS7wLlkwL7t84H0nu6qB1y7yttW7/F5+fm72f75XrZ2Z3Gl3jhKA25tms1nGGAf3fRxjHBFjnBNjXB5CuCKEsKwsy0WNRuMv2n0tASgA6a5OBuC2uB0CUAAKQAFIFzSbzb+MMT5RluU/hxCuiDF+oiiKIsY4vSzLRe0+XwAKQLqrW38GcGvdDgEoAAWgAKTDBg0aNCiEsLIsyzEb+bldYowvtvsaAlAA0l3dCMCteTsEoAAUgAKQDgohfDKE8HRZlocWRVGUZblb378viqIYMGDAh0MIz7b7OgJQANJdnQ7ArX07BKAAFIACkA4py3K3GOOaZrO5Z9+PfeQjH/nTEMJze+yxx98VRVHEGCeEEE7dhK8lAAUgXdTJANwWt0MACkABKADpkBDClBjj6+tf17C077UNjUbjX9Z//EAI4Se77rrrB9t9LQEoAOmuTgbgtrgdAlAACkABSA0JQAFId9X9RdACUAAKQAFIDfUd8euWjGutXnfyRtcuerbVuh17ApBOqHsAjrt5WOuk5yZvsOvXzsxu4QtTk5vw9PjsevFlzsWiJ7Pbaeba5N6z4D+yu+O12cnNH/Xt7KoGXrvIa/cy56qB1y7y6vgy51zgrXlrfuvKl09I7rhnjk/usF8cUcvbQY8QgAKQ7hKAAlAACkABSMcJQAFIdwlAASgABaAApOMEoACkuwSgABSAAlAA0nECUADSXQJQAApAASgA6TgBKADpLgEoAAWgABSAdJwAFIB0lwAUgAJQAApAOq7viF+/5LjWk+vmb3Tdjq5emgBka6t7AKbeA5gLvIUvTG1NfOb45Hr2XX6ZwNv5hLXZ7XrEfcn9+bGrs1v0wtTkfjBzSHbe5df9d/nlAu/Kl09ojXt6fHK5v9c/c9V3ank76BECUADSXQJQAApAASgA6TgBKADpLgEoAAWgABSAdJwAFIB0lwAUgAJQAApAOk4ACkC6SwAKQAEoAAUgHScABSDdJQAFoAAUgAKQjhOAApDuEoACUAAKQAFIxwlAAUh3CUABKAAFoACk4wSgAKS76h6AX1h8UOsrj47aYLnAaxd5vfoy51zg7XbEfdk197w6ud33vTm7454en9wV5wzLzsucu/8y51zgjXt6fGvwg8OTK+46NLmdf7RPLW8HPUIACkC6SwAKQAEoAAUgHScABSDdJQAFoAAUgAKQjhOAApDuEoACUAAKQAFIxwlAAUh3CUABKAAFoACk4wSgAKS7BKAAFIACUADScQJQANJdAlAACkABKADpOAEoAOmuugfg312+X+tT9w/dYP3xXX65wCv3vDq7T//96cn974PPze4rj4xK7rarjsnOu/x64F1+mcAb/ODwVvHLQ9M7//7kdj7x8lreDnqEABSAdJcAFIACUAAKQDpOAApAuksACkABKAAFIB0nAAUg3SUABaAAFIACkI4TgAKQ7hKAAlAACkABSMcJQAFIdwlAASgABaAApOMEoACkuwSgABSAAlAA0nHNZvNzMcbWDy4Z37p56Skb3dIlC2wTd9uNc5ObP2FkdhcvPDu5lStXWj/dZZdd1ooxtprN5ue6fQ82R9/t2POsI1vfvG7CBht945TsDrhhYnJfXnxcdv90xYT0Ljw3v7NuTu6/z707u/9n/DXJ/fsRP8rum/ucltxeB5yV3ZE/n5LcZT+alt0li85J7rwF87M7/YRpyc0bNya7E8celdxp08Zk94MF45O79IdTsrvh6lnJ3XXHydndvuyU5C68a25202+dkdwB10/M7p+uHJfewuuS++yMs2p5O+gRIYS9Y4wtM+vuQgh7d/sebA63w6w3VrfbQY8YMGDAh0MIezebzc+VZflpM+vsms3m50IIew8YMODD3b4Hm8PtMOvu6no7AAAAAAAAAAAAAAAAAAAAaKssy6+FEFbGGB8KIRzX7efpE0K4IMa4KoSwJISwpCzLL3frWcqy/PMY44pGo/HR9R//1xDC0hjjqhjjgqIo/qQHnmlajPGJd/z1OqQLzzQ8hHBfjHFFCOGsoih2aDabfxtCuCOE8ECM8ccDBw58fw88144xxoNijM/0/fUKIUzu9HPVXS/eDnej0nO5HZv2TO4G/ccee+zxVzHG1QMHDvxQURTvDSFcF0L4124/V1EURQjh4UGDBu3S7ecoy/Iz6w/AuvUHc6cY4+qyLJtFURQhhIUxxu92+ZmKGOP1McZPdPI53qnRaPxDjHH5oEGD3rf+GReVZTkshLC00Wh8fv0zTizLclovPFeM8cxuxkHd9ertcDc2+7ncjs14JneDfqMsy33Wf2fT9/G+McYzu/lMRVEUAwcO/FAI4eUQwhUxxuVlWY7v1rOEEM5qNpv/GGN8vNFofLTRaHw+xnhD38/HGAe/8+NuPNP653gxhHBJjHF5CGFeURQ7dfKZms1mGWMc3PdxjHFECOHsEMJj7/g1u7/z4249V4xxzvq/TleEEJaVZbmo0Wj8RSefq+568Xa4G5v/XOufxe3YhGdyN+hXyrIcHUKY1Pdxo9H4lxDC4m4+U1EURQjhb2KMPxo0aNAuAwcOfH+M8YYQwn7dfKYY4xONRuOjzWbzGzHGc/t+vCzLZoxxVTefaf1/OsPlzWZz96Io/iSEcF6McUI3nqkoiqLZbP7l+n+kNCqEcPM7fuq9IYR13X6usiz/eX0kfKIoiiLGOL0sy0Xdeq466sXb4W5s/nO5HZv+TO4G/UqMcey7j3iM8epuPtPGNJvNPUMIF3fzGfoOZghhr3cf8hDCA918pnf/eAjhkzHGe7rxTIMGDRoUQlhZluWYsiw/u5Ej/ma3n2sjP7dLjPHFbjxXXdXhdrgb7Z/r3T/udqSfaSM/525QX+/+xzZlWe5TluX3u/lMRVEUzWbz70MIX+z7OITwbzHGC7v5TH0Hc/0/urm+78cbjcbnQwjXdvOZms1mWZblN/t+vCzLT4cQ7uz084QQPhlCeLosy0OL4o//2ObRdzzXwBDCw91+rrIsd+v790Xxx/9822c7/Vx11ou3w93Y/OdyOzb9mdwN+pWPfexjfx1jfLwsy92KotgxhHBNL/wB1/Xf/T2y2267faAoip1CCNc2m81vdPOZ3vEd804hhCcbjUYsij/8mZoQwtBuPtP6/wvDc7vvvvtHiqJ4TwjhnLIsR3fyWcqy3C3GuKbZbO75rmdc3vdnacqyHB9CmNvt5/rIRz7ypyGE5/bYY4+/W/+ME0IIp3byuequF2+Hu7H5z+V2bPozuRv0OyGEf48xroh/eDXBjG4/T5+yLI8Kf/h//38ohDCl28/zrj80PXj9KwAeCCEsLIpih24/Uwhhr/CHV3KsijGe0elnCiFMiTG+vv6vy9L1/zq5LMuPhxBuDyGsDCFcseuuu36wF55r/Z9ZW7r+f4Y/6fRz9Qe9eDvcjc1/Lrdj05/J3QAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAgG3m/we82ald90S6lAAAAABJRU5ErkJggg==\">"
      ],
      "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,iVBORw0KGgoAAAANSUhEUgAAAoAAAAHgCAYAAAA10dzkAAAgAElEQVR4nO3debRdVYHn8eME2iWWCrFaEjTce/Ze3ZZVKtWsVkuqtatta9Vq26HKtpksJwbDkECYCSFhSmQekjAZISooapgiGCYZZZKMQEgYwlAEghKRQQmD9foP81JpYO+TnJfccw/5fNb6LXlJ7uPkPd7mGxO3RQEAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAC8DpRl+aUQwl0xxiUhhMOqfvzw4cM3DyHs0O12P16W5TZm1tt1u92PhxB2GD58+Oa9OCNSSmeHWavWL2cHfWDrrbf+ixjjQyNGjHh3URRvCiFcFUL4dO41IYQdYowDZtbsQgg79OioeBVnh1l71+TZQZ8oy3KnEML0Nd7eOcb4ndxrut3ux2OMAx/80D4Df73NUT1d+enTs3vvjjOTe8/XZmU3bJefN7ItMtt08vX5fXd8H+6w/M4bldxfXvpVW4t9YNo/D8QYB7rd7sc3/Cnx2oZydnx2+o4DX75y91ftC1fslt1nfr5Lcp+6/JvZfe6K3ZL7+jV7ZPe1q0clt8NVu2eX+/n84+xds/vkZd9I7qMVy/18xl8/Orux1+2d3B6/2DO73MfqtT7nay73OfpixQ66fnRye/9iz+yG8sz/e/auyeU+f5+87BsD/2XW15L7wCX/kt1WF+6U3Ltm7pjcsFM/3/jZQZ8oy/LAEMIRg293Op2/DyHMrnjNNjHGgb/e5qiBj2x79nrfhzMrPzsru7/YfW5ym49ZlN279r+vkb0zs+Lcx/O7ekof7rT8bjo0ub+6e7Stxf7y0q8OxBgHyrLcZsOfEslzoPbZ8eUrdx/45oMHv2o7PXBgdp+/b//k/uHe/bLb/oEDk9vnkXHZjXn40OR2e+jg7HI/n3++/4DsPrNkbHLbLd43u9zPZ8ryI7Kb9NiE5MY9Oj673MfqtT7nay73Odp56UHZHf/4xOQmPDo+u6E88/+9/4Dkcp+/zywZO/DRRWOS++Bde2c3Yv6o5Dabu3ty75q5Y+NnB30ixnjwKw/xGOPludcIQAEoAAXgUM4OASgABaAApGGv/G2bsix3KsvyrIrXCEABKAA38gAcytkhAAWgABSANOz973//e2OMS8uyHFYUxVtCCFeUZfm53GsEoAAUgAJwKGeHABSAAlAA0gdCCP8UY1wYY1wcY5xc9eMFoAAUgAKwKOqfHQJQAApAAUgLCUABKAAFYB0CUAAKQAFIiwlAASgABWAdg2fHF67YbWCnpQe9al+4/4DsPnPvfslVBdFn79s/uV0eOqT2qsIkF3hV0fqJe/ZJ7iN3j84u98yTH5+YXS7wqmJ5lwcPTi4XeNs/cGD2c1T1z8aEZYcnN5Rn/r8PHJhd7vP3sXv2yS4XeFst2CO7zebsnlxx6+jkNv1RO88O+oQAFIACUADWIQAFoAAUgLSYABSAAlAA1iEABaAAFIC0mAAUgAJQANYhAAWgABSAtJgAFIACUADWIQAFoAAUgLSYABSAAlAA1iEABaAAFIC0mAAUgAJQANYhAAWgABSAtJgAFIACUADWMXh2/MPsXdc58D4zxCD69JKxyVWFSW5Vd/nVDbyqn1O4c6/scs9cdZdf3Viqirxc4H32vv0HPrVkbHK5z9+nl4wd0jPXDbyqyKu6y69u4G02Z/eB4pYx6V03KblNZ4xt5dlBnxCAAlAACsA6BKAAFIACkBYTgAJQAArAOgSgABSAApAWE4ACUAAKwDoEoAAUgAKQFhOAAlAACsA6BKAAFIACkBYTgAJQAArAOgSgABSAApAWE4ACUAAKwDoEoAAUgAKQFhOAAlAACsA6BKAAFIACkBYTgAJQAArAOgbPjk9e/s11DryhBlHuX9RVlwznoqUqEIYSrbmfz3sX7JFd7pmbusz5k4v3zS73sdh20Zjs6gZe317mnAu8W8YMFNdNTm/2OcltevakVp4d9AkBKAAFoACsQwAKQAEoAGkxASgABaAArEMACkABKABpMQEoAAWgAKxDAApAASgAaTEBKAAFoACsQwAKQAEoAGkxASgABaAArEMACkABKABpMQEoAAWgAKxDAApAASgAabF+DsDhX705uVwc9uuK8+/O7/Lz2rdfnJDc1gv3tLXZxV9p5SE+eHb81599vedBlAvqqn/J51YV67nP47D5o7J765zdkituHZ1dLmqGEnG5u/o+tWRsNtKqIr3qc1g3eLdbvG92uc9fZcTN3T254rZR+d0wMb0rz8zv4p+nl/l3xqbHX9rKs4M+IQAFoAAUgHUIQAEoAAUgLSYABaAAFIB1CEABKAAFIC0mAAWgABSAdQhAASgABSAtJgAFoAAUgHUIQAEoAAUgLSYABaAAFIB1CEABKAAFIC0mAAWgABSAdQhAASgABSB9JITwwxjj4hDC3BDC3LIsP5f78QJQAApAAVgU9c8OASgABaAApA+EEO4dOXLkO9f2x/dzAL7p1JXJFdNX5HfO4/V33t3p3Tguv8wB8T8Wj23d9nvksCFs3AbbUcsOT+6mZ47P7o8vn57c4w8d1/Ndf9WBfXGI1z47Lv1qrcuc3zt/VHK5fxFvNnf3bIiNefjQ7HL/XFW99psPHlx7xz8+sfZy/zzP/t2x2V3828nJ/eDJY7I744kjk6t65qH8facsPyK5SY9NyG7co+OTm1Cxq393bHLPvTwtu6d+fVJyv39mSnYvv3xGcs+9fHpyv1o4oS/ODvrAiBEj3h1CeCqEMCvGuKAsy/FVrxGAAlAACsAhnR0CUAAKQAFIs0II/znGeMHIkSPfOWLEiLfFGK8JIXw19xoBKAAFoAAc0tkhAAWgABSA9Jdut/v5EMJPcz9GAApAASgAX2mdzg4BKAAFoACkWd1u929CCP9r8O0QwhdjjD/KvUYACkABKACHdHYIQAEoAAUgzSrL8mMhhPuGDRv29qIoNgkhXNntdr9c8RoBKAAF4EYegEM6OwSgABSAApDmlWW5TwhhUYxxSQjhqLX48QJQAArAjTwAV50F9c4OASgABaAApH36OQDfeNrK5Irvrsjv3OWZPZ5f7q6+mw7N7/ZRyX16ydjWbf9/Pawvd/RjE5L75bMnZJcLwHvnHd3zXT5z31Ye4oNnR7jkX9Y58Koi7w137Jpd7r69XR46JLvdH05vlwcPzq7qzr3cJi47vPauePq45C5+6tvZnffkMcmd+eujsssF3mHLDs8uF3hTnzgyu8mPT0wuF3jjHh0/sNfDhya3zyPjsrvkqW8nt/yF07J79P5jk/v1oydk98xvT0nu+T9MTW7unMNbeXbQJwSgABSAArAOASgABaAApMUEoAAUgAKwDgEoAAWgAKTFBKAAFIACsA4BKAAFoACkxQSgABSAArAOASgABaAApMUEoAAUgAKwDgEoAAWgAKTFBKAAFIACsA4BKAAFoACkxQSgABSAArAOASgABaAApMUEoAAUgP++BTcdnt3lM/ZZ7ztz8tdbeYgPnh3DL9yp1mXOb/jVrskVN++f3Vvn7JZcVYjtvPSg5Kpe+9n79k/uH+7dL7uh/PNeN/CqIq/qMudc4FXFVN3Aq4q8XODt9fCh2c9v1S8OZjx5THIPPH9Kdg8umpTcjbMOzu72q8YlN/+Gw5O78Pt7tfLsoE8IQAEoAAVgHQJQAApAAUiLCUABKAAFYB0CUAAKQAFIiwlAASgABWAdAlAACkABSIsJQAEoAAVgHQJQAApAAUiLCUABKAAFYB0CUAAKQAFIiwlAASgABWAdAlAACkABSIsJQAEoAAVgHYNnx7tm7rjOgfenyDsgvesm5Zf5OstF2mfv23/gS/cfkFzVaz+1ZGxy2y3eN7vc/YNVa+ouv1zgVcXUhrrLLxd4Oy89KP+5f+DA7HLROue5E7NbtvTY5G67clx2uci7+7aJyV16wehWnh30CQEoAAWgAKxDAApAASgAaTEBKAAFoACsQwAKQAEoAGkxASgABaAArEMACkABKABpMQEoAAWgAKxDAApAASgAaTEBKAAFoACsQwAKQAEoAGkxASgABaAArEMACkABKABpMQEoAAWgAKxDAApAASgAaTEB2D8BuNMDBzay3GHadOgJwP41eHa89cfbr3vg3XzAQHHd5PSuPDO/GyYml4u0T1V8DX5y8b7ZfeTu0bVXFS655QLvB08ek90ZTxyZ3KTHJmSXi7Sqi6A31GXOVZGei/CqMzh3rlz/9PHZ/WbZCcnlAm/+DYcPLLr9iOTuW3B0cj+/aN9Wnh30CQEoAAWgAKxDAApAASgAaTEBKAAFoACsQwAKQAEoAGkxASgABaAArEMACkABKABpMQEoAAWgAKxDAApAASgAaTEBKAAFoACsQwAKQAEoAGlIWZbviDEu7HQ671v19n8LIcyLMS6OMU4tiuKNa/E+BKAAFIAbWQCuz7NDAApAASgA6aGyLD8aY1wYQli56hDfJMb4UFmW3aIoihDCuTHGXdbi/QjAPgnAPR86JLn9HhnXyJoOPQG4/q3vs2PT8xOxlwu86yYPFFecld6ll+R39WnJbbtoTHYfu2ef5KoiLty5V+194f4Dau+E5Uck98MVk7LLxeGU5UdkVzcON+RdflX3Lf7V3aOTq/pnI3fmXPLUt7N7esXJyeUCb9HtRwzcv+Do5B5ePDm5qy/bTwBujEII07vd7t/GGJd2Op33dTqdT8QYrxn8/hjjdmu+nSIABaAA3LgCcH2fHQJQAApAAUgDYowPdjqd93W73S/HGL83+O1lWXZjjIurXi8ABaAA3LgCcND6OjsEoAAUgAKQBgwe4iGE7V95iIcQFlW9XgAKQAG4cQfgUM8OASgABaAApAGDh/iq37a5evDbO53OJ0IIV1a9XgAKQAG4cQfgUM8OASgABaAApAGDh3hRFJuEEB7udDqxKP7053xCCKOrXi8ABaAA3LgDsBji2SEABaAAFIA0YPAPcq/66+1CCHNDCItCCOcWRfHmqtcLQAEoADfaAFwvZ4cAFIACUADSQgJQAArAjTMAh0oACkABKABpscFD/EMfPmJgm7858zUnAHsTgLkDs+ng6rdtqABcMveo7K6duf9634xTd2nlIb46AM89aKD4xQmvXi7wqiLvpzfmd/l5yf2nu/bOLhcIVRH33gV7JDds/qjscl/7VZdXH7bs8OSm/+bo7KY+cWRykx+fmF0/Xub8V3ftnd1WC/ZIbuuFe2a3y0OHJDfjyWOy+91LU5N7/KHjsnvq1ycl9/tnTkvujtsPa+XZQZ8QgAKwjROAzROAAlAACkBaTAAKwDZOADZPAApAASgAaTEBKADbOAHYPAEoAAWgAKTFBKAAbOMEYPMEoAAUgAKQFhOAArCNE4DNE4ACUAAKQFpMAArANk4ANk8ACkABKABpsdX3AH7kyIGP/JezXnsCUAD22TZUAF5/0YHZnTr2K8mduPeOtXbwv3y+lYf46gD8zpH17vL7yU3p/eC+/C68JrlcAFRFwH+cPyq7t87Zrfaq7qDLbZ9HxiV3yhNHZlc38Bq9y69m4G21YI+Bzebunty75n0ruy89cGByVR/n3yw7IbmLz9w7uwtP3yu5WeeMSe70SV9r5dlBnxCAArCNE4DNE4ACUAAKQFpMAArANk4ANk8ACkABKABpMQEoANs4Adg8ASgABaAApMUEoABs4wRg8wSgABSAApAWE4ACsI0TgM0TgAJQAApAWkwACsA2TgA2TwAKQAEoAGkxASgA2zgB2DwBKAAFoACkxVYH4DZHDSn0NpoAvGVMfkMIwP0eGZdc08HVb2sqAE/Zd+fkNtoAnDZ13QOvKvKqvr4zX6O5AKiKgKqIy319F7eOzi4XNVUXUOcuKJ70+MTsXm+XOW82Z/fssp+H3Ofv9lHZy7gnLjs8u2efOjW5y84dk92s76Z3xXn7Jved47/RyrODPiEABWAbJwCbJwAFoAAUgLSYABSAbZwAbJ4AFIACUADSYgJQALZxArB5AlAACkABSIsJQAHYxgnA5glAASgABSAtJgAFYBsnAJsnAAWgABSAtJgAFIBtnABsngAUgAJQANJi/RyAbzr5D8kVZ/wxvw0VgBWH/FACcPeHD7G13MGPjk/uiqePy+65l6cld8vPD8nu9EO+mtxJo3eqtUO++oVWHuKrA/Dk8+rd5Td9RXpTX8wv9/Vb9TWa222j8rthYnrXTcouFzXD5o/KLnc/XVXEfe3Bg5P7wv0HZJcLoo/cPTq73H2LVZGe/RxUfJyLy89Lb/Y52eU+R1X3HubOlV9ednB21/x4v+RuuPSg5L4/bbdWnh30CQEoANs4Adg8ASgABaAApMUEoABs4wRg8wSgABSAApAWE4ACsI0TgM0TgAJQAApAWkwACsA2TgA2TwAKQAEoAGkxASgA2zgB2DwBKAAFoACkxQSgAGzjBGDzBKAAFIACkIaUZfmOGOPCTqfzvlVvHxNjfDCEMDeEMLcsy2+txfsQgAKwdROAQ7M+zw4BKAAFoACkh8qy/GiMcWEIYeXgIR5jvDrG+KF1fD+NBeDIL12VXfYi6HOW59fCAPzZU99O7vGHjsvunjlHbFSbc91hyc06Z0ztnTr2K9nVvey5ny6CXt9nx6bH/XzdA2/6ioE3THkxuU2OfTq77CXRNx2aXzbwjsjvyjPTq4iL3AXUb5mzW3a5c2P2747N7l/v+3Zyd98+Mbv7Fh6d3MNLJmeXO69WLD8xu9ylyi88PzW7lS+fntyLL07L7oG7jklu2dJjs7v6gv2Su/L8sdldd+EByeV+QfrDs74lADdGIYTp3W73b2OMS9c4xFeEEGbGGBeEEE4uimKTqvcjAAVgGycA61vfZ4cAFIACUADSgBjjg51O533Dhw/fPIRwabfb3aooijeGEH4QY5xQ9XoBKADbOAE4dOvr7BCAAlAACkAaMHiIv/LbQwgfjjHeUfV6ASgA2zgBOHTr6+wQgAJQAApAGjB4iHe73bIsyx0Hv70sy21CCLdWvV4ACsA2TgAO3fo6OwSgABSAApAGDB7inU7nfSGEJ7baaqsti6J4QwjhnLIsD6x6vQAUgG2cABy69XV2CEABKAAFIA1Y8w9yhxC2DyHcFWNcHGM8uyiKN1e9XgAKwDZOAA7d+jo7BKAAFIACkBYSgAKwjROAzROAAlAACkBabPAQL//nGQPl5y7r6bqfvzy7bABmA2+oG8Ils0MIwMueOja5qgBcePOE5G64+KDscodPv+7yGfskVxV5l353dHKn7LtzdgLw360OwClnDxQXzX71Lrg9vxmPpHfmC9m96eTfJ1f84oT8soF3Vn6XXprej2/L7g2/2jW54ub9s9t20ZjkfvTkpOye+Nfjk7vnjiOyu3X2ocndftW47ObfcHhyd94yMbv7Fx6T3CNLvp3d8oePT67qtbnA+82yE7K7/qIDk6s6z3KRN+faw5L76Yw9Wnl20CcEoABsOuYEoAAUgAJQAApAekwACsCmY04ACkABKAAFoACkxwSgAGw65gSgABSAAlAACkB6TAAKwKZjTgAKQAEoAAWgAKTHBKAAbDrmBKAAFIACUAAKQHpMAArApmNOAApAASgABaAApMcGD/F37XXtwNsPX9ZXy/7Lo+KwLS47v/5uGZNenwZgLvByB9P1Fx1YeUdVE7uiYkO5608Arh+rA/DH2w8Uv9r11btxXH7XnJhe1dfozGvSy93zd9uo/B2Bua/9W8YMFNdOTu+1PgZrLvf3rYjW/3TX3slN/83R2T3165OSywXerbMPzcZH1WvnXjc+uaoAvHfe0ck9dM/k7B578Ljae/KxE5PL3U347FOnDtx25aHJ5QLvlp8fkv1Y3X3bxOQuvWB0K88O+oQAFIBNx54AFIACUAAKQAFIjwlAAdh07AlAASgABaAAFID0mAAUgE3HngAUgAJQAApAAUiPCUAB2HTsCUABKAAFoAAUgPSYABSATceeABSAAlAACkABSI8JQAHYdOwJQAEoAAWgABSA9JgAFIBNx54AFIACUAAKQAFIj60OwL2vG3j7hMd6uj+bsCy74se3Z1YRgJefn9l5+d06JjMB2JMAPG/f7ARg81YH4HmJqMoF3jUn5r9GL/xFfuffnV5VeF59WnrXTc5v9jnp3TQuv9zft+JMGjZ/VHInLD8iu98/c1pyucCbc+1hAwtuOjy5OdeOzy4XeEvmHZXdUAJvxfITk8sF3pOPnTjw3O9OTe6lF6dlN+/68cnlAm/udeMHFt1+RHIP3HlMcldcMraVZwd9QgAKwKZjTwAKQAEoAAWgAKTHBKAAbDr2BKAAFIACUAAKQHpMAArApmNPAApAASgABaAApMcEoABsOvYEoAAUgAJQAApAekwACsCmY08ACkABKAAFoACkxwSgAGw69gSgABSAAlAACkB6TAC+xm6YuEH2kbtHZ3feiknJPb3i5OxyAVi1qkDsxw0lDnMBeOrYr2RXFYh1Nu4bX2zlIb46AM85eN0Dryryzl+U33dXpHf1lPx+9pP0Zp+b38U/Ty/3fn/2k4Hipzem972H88vcIbj7w4dk9+SLU5Jb+MsJ2S2ec2RyDy6alN2ypccm95tlJ2SXu2/vheenZrfy5WnJVUVcbs+9nN/iO45MrupjlfvFfe4ex1tuPKSVZwd9QgAKwKZjTgAKQAEoAAWgAKTHBKAAbDrmBKAAFIACUAAKQHpMAArApmNOAApAASgABaAApMcEoABsOuYEoAAUgAJQAApAekwACsCmY04ACkABKAAFoABkLZVluW8I4c4Y48IQwvSiKN7c7XY/GEK4JYSwKMb44xEjRrxtLd6PABSArZsArG99nx0CUAAKQAFIj3Q6nW1jjAtGjhz51qIoirIsZ5RlOSaEMK/T6XyiKIoixjixLMtjqt6XABSATcecAOxdAG6Is0MACkABKADpkW63W8YYtxt8O8Y4NoTw3RDCA2v8mK3WfDtFAL7GbhmTXvaS6DED75r3reROWX5Edr97aWpyVQf1UAJvKJcqt3FDCcBLzp6W3P0L59fa7Isv7NkhviHOjk3PPH7dA68q8qavyG/qi+nN+ml+uXPjotn55S6gvuTS/HJ/3wuvye4L9x2Q3OzfHZvdsgePSy4XeFWRN5R4XHrXpOxyFz1XX+Z8WnIvvXR6drnAW7by1OxyH6vjRn0ju7tu+2VyC266LrmLzpshADdm3W73PTHGB0MIB4QQbljju94UQlhZ9XoBKACbDjIB2LsAXNP6OjsEoAAUgAKQHhs5cuTIEMJdZVkeVJblx17jEP9D1fsQgAKw6SATgL0PwPV5dghAASgABSA9FEL4cAjh0bIsRxXF6t+2uX/w+8uyHBFCuLfq/QhAAdh0kAnA3gbg+j47BKAAFIACkB4py3JYjPHxbrf7+TW/Pca4YPDP95RlOT6EcNJavC8BKAA3qm3MAbghzg4BKAAFoACkR0IIR8UYnwshzA0hzFv1n0eWZfmBEMLNIYS7Qgiztthii82q3pcAFIBNB5kA7F0AboizQwAKQAEoAGkhASgAmw4yAdjM/whkqASgABSAApAWGzzEN9/tioE/P2hpj/dAdn15D+Cto7PLBuATR2aXC8AFNx2enQBsPgAfuHNBrV1xyUWtPMRXB+CJMweKH8179c59PL8z/pjcW054Jrv8LxwrzobcfXu5wDv/7vz9g1V3+eU+Fj+4L7tPLRmb3I+enJRd7r69odzll7v3bvEd+cjLvd9lS4+tCLxTsxvKXX6PvnBacgueOym73F1+k3b5Snbzb7ouuTnXXpXcT2dMb+XZQZ8QgAKw6SATgAJQAApAASgA6TEBKACbDjIBKAAFoAAUgAKQHhOAArDpIBOAAlAACkABKADpMQEoAJsOMgEoAAWgABSAApAeE4ACsOkgE4ACUAAKQAEoAOkxASgAmw4yASgABaAAFIACkB4TgAKw6SATgAJQAApAASgA6bHVh/jRN6cvXb3g9mbWlxdBb7gAzB1MGzIALz5z7+QuPH2v191mTtszuRP33jE7AfjvVp8d5xw8UFxz4qt347j8rpqSXtVlzrlz43uP5HfWC+lVXED9xtNWJlecszy/3OXV33k2u4/cPTq5GU8ek93TK05O7rGlx2WXi8d75x2V3f0Lj0nu0fuPze7Xj56Q3DO/PSW7F1ZOS64yAFeemtyC507MLhe0E3f+P9nd8Yurkrv1isuT+9F3zmzl2UGfEIACUAAKwDoEoAAUgAKQFhOAAlAACsA6BKAAFIACkBYTgAJQAArAOgSgABSAApAWE4ACUAAKwDoEoAAUgAKQFhOAAlAACsA6BKAAFIACkBYTgAJQAArAOgSgABSAApAWW32IT1yYvpcrG2JNbQMG4K1jMttwAfjHl09PbkMG4EVn7JVc1b14Te3Yb22fXNWdW0PZFefNSG6jDcApZw0UF//81bt6Sn65wKu8yy8TeGe+kN2bT3g2uTedujK7tx79m+TeMOXF7Dad/Nvkqu493HrhnslVBeALz09NLnff3pOPnTjw7FOnJld1h+DVF+yXXO68uuHigwZuv2pccvNvODy7xXOOTK7qtVeePza5qme+7cpDkztu1Deyy0XeTT+7JLnzzpjayrODPtblLbsAAA6jSURBVCEABaAAFIB1CEABKAAFIC0mAAWgABSAdQhAASgABSAtJgAFoAAUgHUIQAEoAAUgLSYABaAAFIB1CEABKAAFIC0mAAWgABSAdQhAASgABSAtJgAFoAAUgHUIQAEoAAUgLSYABaAAFIB1CEABKAAFIC22+hA/6YKaFzILQAG4/nfKvjtnlwvAqQeOye7S75yR3P0L59feRhuA08fXu8y5buBVRF5VTP3ZhMeSe9vRv8nuzw9amlwu8Dad/NuBzQ59OLk/m7Asu83m7p5c1bny8h9PTy4XeM8+dWo2Hh9/6Ljs6sbSbVcemg+8O47MLhel99xxRHY/PPFbyR399R2zu+MXVyZ3y+zLsrtx1sXJXTvzguRmnHpSK88O+oQAFIACUADWIQAFoAAUgLSYABSAAlAA1iEABaAAFIC0mAAUgAJQANYhAAWgABSAtJgAFIACUADWIQAFoAAUgLSYABSAAlAA1iEABaAAFIA0oCzLfUMId8YYF4YQphdF8ZYY464xxmUhhLmrduRavB8BKAAF4EYUgOv97BCAAlAACkB6o9PpbBtjXDBy5Mi3FkVRlGU5oyzLMTHG75Rl+bl1eV+rD/GTL+iDqBOAAnDDB+Cs6WcmVzfihrJeBuAGOTtOP6Xnd/nlAu/tEx7L7t1jlyT3zoOWZveePRYk945DHspui73uTO7dYxdnlzuTJi47PLsVL01NLhd4Lzw/dWDly9OSe+rXJ2VXN/CqIu+heyZll4vSqrsLZ07bM7kJO/1zdrm7/HKBVxV5V/3w+8lNP+FYAbix6Xa7ZYxxu8G3Y4xjY4wnxhgXhBBmhRDml2U5o9Pp/HnV+xKAAlAAbjwBuEHODgEoAAWgAKT3ut3ue2KMD5Zl+ckQwqwY44eKoihijJPKspxR9XoBKAAF4MYTgGtab2eHABSAAlAA0lsjR44cGUK4qyzLg17j+94ZY1xR9T4EoAAUgBtfAK7Xs0MACkABKADpnRDCh0MIj5ZlOaooiqIsy2GDf10URTF8+PDNQwjLq96PABSAAnDjCsD1fnYIQAEoAAUgvVGW5bAY4+Pdbvfzg9+25ZZb/ocQwhNbb731XxdFUcQYJ4QQTl+L9yUABaAA3EgCcIOcHQJQAApAAUhvhBCOijE+t+q6hnmD1zZ0Op2/X/X2ohDChVtsscVmVe9LAApAAbjxBOAGOTsEoAAUgAKQ9hGAAlAAbjwBuD4JQAEoAAUgLSYABaAAFIB1/P9nxzoG3hAvc84G3r6Ls/uPu/wqufeMmp/diJ1vTG7Ynguzy712y6/fml3xixOS2+eRcdk9sPLU5HKBt/LlaQPLX5iS3PO/z28olznnAm/5w8dnl4vSJ/71+OwuO3dMcsd888vZ1Q28a2deMHDl+d9L7rJzz07ujMlHtfLsoE8IQAEoAAVgHQJQAApAAUiLCUABKAAFYB0CUAAKQAFIiwlAASgABWAdAlAACkABSIsJQAEoAAVgHQJQAApAAUiLCUABKAAFYB0CUAAKQAFIiwlAASgABWAdAlAACkABSIutPsSnThsoLrn09bPXWQDmAq9qAlAAbgj/HoDnDRQ/uenV+/79+eUC8PhnsvuzI5Ynt/m+i7N77zdvT26rnW7IrvzsrOS22vnG7HKv7fzTFdkVP/tJcjsvPSi7bOC9OCW7Jc+fktzLfzw9u6Hc5ZcLvN/95uTscvcaPr3i5Oyu+fF+yR2/xw7Z1Q28qsi75OxpyU09akIrzw76hAAUgAJQANYhAAWgABSAtJgAFIACUADWIQAFoAAUgLSYABSAAlAA1iEABaAAFIC0mAAUgAJQANYhAAWgABSAtJgAFIACUADWIQAFoAAUgLSYABSAAlAA1iEABaAAFIC0mAAUgAJQANYhAAWgABSAtNjqQ/yHXxsobhljt+TiTwC+3gPwt08sz+7FlSvX++bPndvKQ3z12TF9/EBx1ZRXr+oXaRfcnt45y7N7w5QXk3vrpBXZveOQh5Krusw5F4fhs7Oy++B230/uA39/QXbFj+Yl9w/37pddLvAW/+Hk7H757AnJrXhxSnZDucy5buC98PzUgd+/PC25l17M75eXHZzcaQf8S3YP3XN3cn949tnsnvntiuRy59HNN1zfyrODPiEABaAAFIB1CEABKAAFIC0mAAWgABSAdQhAASgABSAtJgAFoAAUgHUIQAEoAAUgLSYABaAAFIB1CEABKAAFIC0mAAWgABSAdQhAASgABSAtJgAFoAAUgHUIQAEoAAUgLbb6EL9gp4HitlFWtdvzGzY/vaEEYFXEXXfhAcldO3P/7GZO2zO5jS0An/r1E9n927/923rfnXfe2cpDfPXZccbxA8Vl5796r3U34NreE5i5I7DqnsC3H74suy1G353c8K/dml33i7OT+8v//qPstvmbM5P78LZnZ1d8d0Vy2y4ak13uLr/5z52Y3TVPH5fcspWnZreh7vJ77uVp2T38wmnJ/e6lqdnNu358cmcd9vXscnd9Vp0BL7/0UnIr//D75ObOuaOVZwd9QgAKQAEoAOsQgAJQAApAWkwACkABKADrEIACUAAKQFpMAApAASgA6xCAAlAACkBaTAAKQAEoAOsQgAJQAApAWkwACkABKADrEIACUAAKQFpMAApAASgA6xCAAlAACkBarNvtfjzGOPCxM/YY+LufjbOqXZbfP85O77Sbj8tu4cLTk7vovKNq78IfHJndD6aNT27KEfs3stMm5nfSuLHJTT1qQnY/POuM5G6/5ebs7rzzzvW+Sy65ZCDGONDtdj/e9HmwLlafHSdNHvi7C8599b53SX7fvTa9M+Zn999PviO5/zn5l9n944SrkvvcQbOy+6fRP07uy7vMyG7Hnc5I7ytnZvd3Z9+W3FeunpDdNfNOS272nFOym/mrk5P75bwp2d18w8nJ/erW07KbP29qcnMWTsvuhvlTkvvVwqnZzb5ocnJnTD4ou+uuuCK5qjNg4cKFyc2bOze5mTNntvLsoE+EEHaIMQ6YWbMLIezQ9HmwLpwdZv2xtp0d9Inhw4dvHkLYodvtfrwsy23MrLfrdrsfDyHsMHz48M2bPg/WhbPDrNm19ewAAAAAAAAAAAAAAAAAAACgUlmWXwoh3BVjXBJCOKzp5xkUQvhhjHFxCGFuCGFuWZafa+pZyrJ8R4xxYafTed+qt/9bCGFejHFxjHFqURRv7INnOibG+OAaH69vNfBM+4YQ7owxLgwhTC+K4s3dbveDIYRbQgiLYow/HjFixNv64LneEmPcNca4bPDjFUI4stfP1Xb9eHY4N2o9l7Nj7Z7JucHrx9Zbb/0XMcaHRowY8e6iKN4UQrgqhPDppp+rKIoihHDvyJEj39n0c5Rl+dFVB8DKVQfmJjHGh8qy7BZFUYQQzo0x7tLwMxUxxqtjjB/q5XOsqdPpbBtjXDBy5Mi3rnrGGWVZjgkhzOt0Op9Y9YwTy7I8ph+eK8b4nSbjoO369exwbqzzczk71uGZnBu8bpRludOqX9kMvr1zjPE7TT5TURTFiBEj3h1CeCqEMCvGuKAsy/FNPUsIYXq32/3bGOPSTqfzvk6n84kY4zWD3x9j3G7Nt5t4plXPsSKEMDPGuCCEcHJRFJv08pm63W4ZY9xu8O0Y49gQwndDCA+s8WO2WvPtpp4rxnjiqo/TrBDC/LIsZ3Q6nT/v5XO1XT+eHc6NdX+uVc/i7FiLZ3Ju8LpSluWBIYQjBt/udDp/H0KY3eQzFUVRhBD+c4zxgpEjR75zxIgRb4sxXhNC+GqTzxRjfLDT6byv2+1+Ocb4vcFvL8uyG2Nc3OQzrfp/Z7i02+1uVRTFG0MIP4gxTmjimYqiKLrd7ntW/ZbSASGEG9b4rjeFEFY2/VxlWX5yVSR8qCiKIsY4qSzLGU09Vxv149nh3Fj353J2rP0zOTd4XYkxHvzKQzzGeHmTz/Raut3u50MIP23yGQYPzBDC9q88yEMIi5p8pld+ewjhwzHGO5p4ppEjR44MIdxVluVBZVl+7DUO8T80/Vyv8X3vjDGuaOK52qoNZ4dzo/q5Xvntzo70M73G9zk3aK9X/rZNWZY7lWV5VpPPVBRF0e12/yaE8L8G3w4hfDHG+KMmn2nwwFz1WzdXD357p9P5RAjhyiafqdvtlmVZ7jj47WVZbhNCuLXXzxNC+HAI4dGyLEcVxerftrl/jecaEUK4t+nnKsty2OBfF8Xq/3/b5b1+rjbrx7PDubHuz+XsWPtncm7wuvL+97//vTHGpWVZDiuK4i0hhCv64Q+4rvrV333Dhg17e1EUm4QQrux2u19u8pnW+BXzJiGEhzudTiyKP/2ZmhDC6CafadV/w/DEVltttWVRFG8IIZxTluWBvXyWsiyHxRgf73a7n3/FMy4Y/LM0ZVmODyGc1PRzbbnllv8hhPDE1ltv/dernnFCCOH0Xj5X2/Xj2eHcWPfncnas/TM5N3jdCSH8U4xxYfzT1QSTm36eQWVZ7hP+9D//XxJCOKrp53nFH5rebtUVAItCCOcWRfHmpp8phLB9+NOVHItjjGf3+plCCEfFGJ9b9XGZt+o/jyzL8gMhhJtDCHeFEGZtscUWm/XDc636M2vzVn0OL+z1c70e9OPZ4dxY9+dydqz9Mzk3AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAABgg/l/ydImzGLG81YAAAAASUVORK5CYII=\">"
      ],
      "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,iVBORw0KGgoAAAANSUhEUgAAAoAAAAHgCAYAAAA10dzkAAAgAElEQVR4nO3de7CeZX3v/7secU/ppkq6qwm48jz3de3Wza6H1vkplb3dv9ZfnY4dbeuJk4ASTgESEiBEjgmnEEggHASEAJEzchAwIIcoKCIiBBIChrOIFChFoWJlz+zt8/vDrJAhua8nPGsl131lvV4zn5GVZK3c0OE775j0tqoAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAACATUBd158LIayIMT4cQji8348fP378u0IIO3S73W3ruv6QmW3cdbvdbUMIO4wfP/5dG+NGNKndDrOi1pbbQQtMnDjxv8QYfzZhwoR3VlX15hDCLSGET6Q+J4SwQ4yxZ2Z5F0LYYSOdirW4HWblLuftoCXqut4phLBwjY93jjGem/qcbre7bYyx9+d/dXDvfR+Zu86N+8rixr399HMK3NfTu3hy4478/tSB9427D9lgu+wnM1u36+49tJ2756uNu/P2mcldc/F+jTt99i4D7Ygp/9yLMfa63e62G/5KrNtIbsfmx07pvX3h7LU27qodk/vHm/Zs3N9/Z4/k/vaG3Rv3PxZ/Jbn/59tfbtwHr9s1uf967S6N6/fzfvMnMxt38d2HJHf7fYc37palhyW3OLV7D03u6p/MbNyldx+S3KIfz2jcuX123l0HN+70Ow9Mbt4Ppzduzh3Tkjvi9imNO6TPblp6WOOW3nNEcldfuF/jrr10SuPOmLdr9ttBS9R1PSOEMHv4406n8zchhO/0+ZwPxRh77/vI3N422124zr1z+srGVd/6ToG7Mb07DmvcGc8fPfBu//eTNtjuemV+6/bgf5zSzv3m5Ma9+Nz85JbdcWTjvnX2/gPt9Nm79GKMvbquP7Thr0TjHRj4drx94exedfPZa+2P79s7uZ0en9G4zz52cHKfeuTAxv3tymnJfeynBzTurx6cktyfrdi/cZ94eHpy97wyv3E//PW85J7936c17mevnprcw79d0LiV/3FKcve/Mr9xP/r1vORue/nExt3UZ7e+NLdx3/rlnOQu/bfjGnfBC8cmd/pzsxs3r88e/+2pjfvtb05P7v7vH9m4B340q3HXXjol++2gJWKMM19/xGOMN6Q+RwAKQAEoAEdyOwSgABSAApDMXv/bNnVd71TX9df7fI4AFIACcIwH4EhuhwAUgAJQAJLZe9/73nfHGJ+o63pcVVVvDSHcVNf1p1OfIwAFoAAUgCO5HQJQAApAAUgLhBD+Oca4PMa4MsY4p9+PF4ACUAAKwKoa/HYIQAEoAAUgBRKAAlAACsBBCEABKAAFIAUbcQAuvqS83XBxenfv07i2BuCK35zcumUPvQ0QgF87ZNfGnTR5h+SO3u3z69zUz32qyCO+OgDPPmmd/x5tdu+eye36xCGN++JjByf3mUcPatzfPTw9uY+vnNa4fgH43xPrF4CpwFvy8onJDRp4D/92QW/Zb05u3GO/XZBc6t/vVNDe88r83h3/flLjUoF360tze9f+6oTGXfbi8cmlAq/fjT7p2VmNm/XMkcmlYvm6c6ckd8Wpkxt34dw9Gzf34B2KvB20hAAUgAJQAA5CAApAASgAKZgAFIACUAAOQgAKQAEoACmYABSAAlAADkIACkABKAApmAAUgAJQAA5CAApAASgAKZgAFIACUAAOQgAKQAEoACmYABSAAlAADkIACkABKAAp2IgD8IZLCtymF4DZo6qkjSAAz5ixS+NO3Gf75DbZADz96+t8n+Yf3LNHcts/PqNxqcD7zKMH9T75yIGNSwXex1dO6334oamN22bF/slNXL5v47ZbOS25VOAtfmlucqnAu/83Jyc3knd13vvK/Mb1e3fhrS+f2LhU4PWLvH7v8hs08PpF3kFPH55c6p2I1503Jblvnja5cZfO37tx82buWOTtoCUEoAAccxOAo0IACkABKAApmAAUgGNuAnBUCEABKAAFIAUTgAJwzE0AjgoBKAAFoACkYAJQAI65CcBRIQAFoAAUgBRMAArAMTcBOCoEoAAUgAKQgglAATjmJgBHhQAUgAJQAFIwASgAx9wE4KgQgAJQAApACiYABeCYmwAcFasDcMFFverKH6y9Ow9KbtDA6xd5qcDrF3mpwJu4fN/eu5dNbly/n3fxr05o3FW/nJNcMvB+PS+51N1IBV6/yEsFXr/Ia+vLnA98+rDG7fXUV5Nb8tLcxt144QHJXXPWfo27fME+jTv5sJ2KvB20hAAUgGNuAnBUCEABKAAFIAUTgAJwzE0AjgoBKAAFoACkYAJQAI65CcBRIQAFoAAUgBRMAArAMTcBOCoEoAAUgAKQgglAATjmJgBHhQAUgAJQAFIwASgAx9wE4KgQgAJQAApACjZ8xN//gdm9D/3l2evcn066u3Hrev9X8bvj0Mb1OyCpbcgAtDewl5v3zBNzk0sFYFPg9VvxAXjcnb3qvBfX3s1nJ5d6Z16/mPqzFfs3biQRt/nSvZJLvddwq2WTk1v4wrGN6xc1qcC77d9PSu6ml+Y2biTv40v9/fT7e5rz7KzkRvI+vlSk7fbkzOQ+9/iMxvX7Rckd/35S4265bHpyV5w6uXFXfW3fxp161JeKvB20hAAUgGNuAnBUCEABKAAFIAUTgAJwzE0AjgoBKAAFoACkYAJQAI65CcBRIQAFoAAUgBRMAArAMTcBOCoEoAAUgAKQgglAATjmJgBHhQAUgAJQAFIwASgAx9wE4KgQgAJQAApAWiSEcGmMcWUIYWkIYWld159O/XgBKADH3ATgOg16OwSgABSAApAWCCE8MjQ0tMX6/vjVR/zMU3rV9d9c9xZfsmltBC+CvujfjkvutpdPbNyGjJrHfrugcbd/a0Zyp07/UuNOmryDrbE5e3yxcZtAAA50O/72hkm9Tz1y0Fr7u4cPTC71Mue/enBKcqmXOQ8tm5zcn96/T+P+8N49k3trYv3C8zsvzW3cj349L7lBA69f5KXu1W0vn9i793uHN27h7N2TO+vQ3Rp3zlFfSS71dc8/dlJy35i7Z+MunrdXcpefOrlx3zx93+S+dfb+jbv6zP2SS0Ve6iXRp8/apRW3gxaYMGHCO0MIvwohXB9jXFbX9RH9PkcACkABKABHcjsEoAAUgAKQzEIIfx5jvHxoaGiLCRMmvCPGuCSEsGvqcwSgABSAAnAkt0MACkABKABpmW63+5kQwpWpHyMABaAAFICv90ZuhwAUgAJQAJJZt9v9yxDCp4Y/DiH8U4zxstTnCEABKAAF4EhuhwAUgAJQAJJZXdcfDSE8Om7cuD+squptIYSbu93uF/p8jgAUgAJwjAfgSG6HABSAAlAA0gJ1XR8QQngoxvhwCOGY9fjxAlAACsAxHoCrbsFAt0MACkABKAAp0OoAPOuUXvXtb657S+ZvWvvuvOQ2u3fPxn37VycktyEjb+AAvGZGcgum7dy43MHVtp24z/ajvhk7fbrIIz58O/568Vd6/+vh6Wvt7/osFYAfeWhqcqkArJfvm1wqALdYuldyIwnAJS+f2LiRvMvvW786IbnUu/z6xeOyO45s3PnH7J7cOUd+uXHnzd49uQuOndS4C0/YI7lL5u3VuMsX7JPclafv27hUiF1z1n69686d0rjFF0xN7qaLpzVuyRUHNu78U3Yv8nbQEgJQAApAATgIASgABaAApGACUAAKQAE4CAEoAAWgAKRgAlAACkABOAgBKAAFoACkYAJQAApAATgIASgABaAApGACUAAKQAE4CAEoAAWgAKRgAlAACkABOAgBKAAFoACkYAJQAApAATgIASgABaAApGDrE4D18v0at82KKZvc/vGxgxvX78Wpm1oAztr58wPv0vknJLf0tiV225LeVYvOL/KID9+OD3/7y72P/vSAtXbI04cPvKlPHZrc7k/ObNxOj89I7rOPHdy4kUTrzk8cklzqbvQLsVTgXfri8cktfOHYxl37yznJPfjjWY1bdNyk5Obvv2PjUi+fP3X6l3pnH3Zw435w/bcG3vevvWrg3Xb1FQPvu1deltytl1/cuCVXNO/8U04q8nbQEgJQAApAATgIASgABaAApGACUAAKQAE4CAEoAAWgAKRgAlAACkABOAgBKAAFoACkYAJQAApAATgIASgABaAApGACUAAKQAE4CAEoAAWgAKRgAlAACkABOAgBKAAFoACkYOsTgOGB/Rr33x+cssktFYC5Am9TDMD7bv+u3f7d3lXfuKDIIz58O/7iul17H3xwylo74OeHJTc1sUk/+2pyqdBKBd5nHzu498lHDmxcKvA+vnJa78MPTW3c5x6fkdxI3uU3aOAtfOHY3hnPH924y/7tuOR+es/sxl180p7Jpe7KaQfvktyZM6c3rm/oXXd14/pF3u3XfLNx/SLve1de1rglV1ySXCoAb7n0wsYtnDe3yNtBSwhAASgABeAgBKAAFIACkIIJQAEoAAXgIASgABSAApCCCUABKAAF4CAEoAAUgAKQgglAASgABeAgBKAAFIACkIIJQAEoAAXgIASgABSAApCCCUABKAAF4CAEoAAUgAKQgglAASgABeAgBKAAFIACkIIJQAEoAAXgIIZvR7h2l3Xehr2e+mpygwZev8hLBV6/yEsF3ocfmpq8G/1+3g31MufTnz86uTnPzmrcwn89NrlHlzXv0vl7J3faQbs07owZu/bZ1MalAq9f5KUCr1/kpQKvX+SlAq9f5N108aLGnXvSnCJvBy0hAAWgABSAgxCAAlAACkAKJgAFoAAUgIMQgAJQAApACiYABaAAFICDEIACUAAKQAomAAWgABSAgxCAAlAACkAKJgAFoAAUgIMQgAJQAApAMqnr+o9ijMs7nc7Wqz7+nyGE+2KMK2OMZ1RV9ab1+BoCUAAKwDEWgKN5OwSgABSAApCNqK7rj8QYl4cQXl11xN8WY/xZXdfdqqqqEMIFMcZJ6/F1BKAAFIBjKABH+3Zsfc3Ova2WTV5ruz05M7k2vstvmxX7Jzdx+b6N227ltORG8i6/QQNvzrOzerOeObJxZzx3dHJPPnR84y5fsE9yXztk18adfdiXkzvtwP0bt6m9y++WSy/s3XTRBY278RvnNe7rJxwnAMeiEMLCbrf71zHGJzqdztadTudjMcYlw98fY9xuzY+bCEABKADHVgCO9u0QgAJQAApAMogxPtnpdLbudrtfiDF+Y/jb67ruxhhX9vt8ASgABeDYCsBho3U7BKAAFIACkAyGj3gIYfvXH/EQwkP9Pl8ACkABOLYDcKS3QwAKQAEoAMlg+Iiv+m2bW4e/vdPpfCyEcHO/zxeAAlAAju0AHOntEIACUAAKQDIYPuJVVb0thPBUp9OJVfX7P+cTQpjS7/MFoAAUgGM7AKsR3g4BKAAFoAAkg+E/yL3qr7cLISwNITwUQrigqqq39Pt8ASgABeCYDcBRuR0CUAAKQAFIgQSgABSAYzMAR0oACkABKAApmAAUgAJQAA5i+Hb86dU79cbdv89a2/7xGcnlepnzoIE3cfm+6/z7HF6/n3dDvcz5qGeOTO6gpw9v3Jx/mZXc04+e0LirvrZvcmcftlvjzjnyK8ktmDa5cUW+zDkReP0iTwCywQhAASgABeAgBKAAFIACkIIJQAEoAAXgIASgABSAApCCCUABKAAF4CAEoAAUgAKQgglAASgABeAgBKAAFIACkIIJQAEoAAXgIASgABSAApCCCUABKAAF4CAEoAAUgAKQgglAASgABeAghm/H5lfu0Nvs3j3X2j88elBybXyX35b37Z3cuv4+h/dnK/ZPbiTv8ksF3oFPH5bcXk99tXGH/+LI5J792YmNu+as/ZI798ivNG7hrN2TO3nKXo3b1N7ld+M3zuvdsGhh4xZfcE7jzppzTJG3g5YQgAJQAArAQQhAASgABSAFE4ACUAAKwEEIQAEoAAUgBROAAlAACsBBCEABKAAFIAUTgAJQAArAQQhAASgABSAFE4ACUAAKwEEIQAEoAAUgBROAAlAACsBBCEABKAAFIAUTgAJQAArAQQhAASgABSAFE4ACUAAKwEGsvh2X79Sr7t5nrX3i4enJ5XqZ86CBt9m9e/aqu/du3FbLJic3kpc5Dxp4ez311d5uT85s3AE/Pyy5F56Z17jrzpuS3Hmzd2/couMmJTdvvz0aV+LLnFOBd8Oihb3F55/TuOsXnt24M4+dXeTtoCUEoAAUgAJwEAJQAApAAUjBBKAAFIACcBACUAAKQAFIwQSgABSAAnAQAlAACkABSMEEoAAUgAJwEAJQAApAAUjBBKAAFIACcBACUAAKQAFIwQSgABSAAnAQAlAACkABSMEE4BsLwGf/92mt3P/9P2c2rt/7ulIRd/rBU5I7Y8bUxvV7b9bPH1lpj6zsLblhcZFHfPXtuGRSr7rzoLX20Z8ekFwb3+VX/Xif9O44tHF/fN/eyV354vGNu+CFY5M76/mjG3fSs7OSO+oXRzTujOdmJ/fTe5q3cPbuyc3Z44uNO2nyDsldfOJxjVt570+Se+juuwbeA3f9cOAt/+EPBt79P7htoF190aIibwctIQAFoAAUgIMQgAJQAApACiYABaAAFICDEIACUAAKQAomAAWgABSAgxCAAlAACkAKJgAFoAAUgIMQgAJQAApACiYABaAAFICDEIACUAAKQAomAAWgABSAgxCAAlAACkAyqev6j2KMyzudztarPj4uxvhkCGFpCGFpXdd7r8fXEIACUACOsQAc1dshAAWgABSAbDx1XX8kxrg8hPDq8BGPMd4aY3z/G/w6AvANBOD/+b9nJvfi8ydvmD03P7knHzq+cadO/1Jy5x19eOO+9fUzkrv2nK817uGl9yT3u9/9zn73u94DDzywUY/4qN+ORdN71W3Hr7UPPjgluWwvc04G3mHpfXde4956757J3fTyiY1Lvcj9sd8u6N3/yvzG/ejX85K79aW5jVvxm5OTu/HCAxp39K5fTO6sQw9q3Jkzpyf3k1tvalzuf1/btI19O2iJEMLCbrf71zHGJ9Y44i+GEK6KMS4LIZxSVdXb+n0dASgABeDYOuKjfjsEoAAUgGPidtAyMcYnO53O1uPHj39XCOG6bre7VVVVbwohXBRjPKrf5wtAASgAx+YRH7XbIQAFoAAcU7eDlhg+4q//9hDCB2KM9/T7fAEoAAXg2Dzio3Y7BKAAFIBj6nbQEsNHvNvt1nVd7zj87XVdfyiEcFe/zxeAAlAAjs0jPmq3QwAKQAE4pm4HLTF8xDudztYhhOe32mqr91RV9QchhPPrup7R7/MFoAAUgGPziI/a7RCAAlAAjqnbQUus+Qe5QwjbhxBWxBhXxhjPqarqLf0+XwAKQAE4No/4qN0OASgABeCYuh1sIgSgABSAjvggBKAAFIBuBwUTgKMbgL98/uQNsg0ZgOcfc0TjUoHXb/0CkN9bsWJFkUd89e1YOLtX3Xz2WkvdjfDAfhlf5pwKvPnpfef85t01Nbmrfjmncct+c3JyqcBb8vKJyS1+aW7j+sXjTRdPa1zqBfKzdv58MvBSL5A/Y8bU3l033dA4XlPq7aAlBKAAFIB5lXrEBaAAFIB5lXo7aAkBKAAFYF6lHnEBKAAFYF6l3g5aQgAKQAGYV6lHXAAKQAGYV6m3g5YQgAJQAOZV6hEXgAJQAOZV6u2gJQSgABSAeZV6xAWgABSAeZV6O2gJASgABWBepR5xASgABWBepd4OWkIACkABmFepR7zf7Xj3ssnJ5XuXXyrwLkjvumubt453Ia65i19sXr8QSwXet391QnKXvXh845a8NDe9Kw5s3OxdPpfcGQdPbdxpB+6fnABcP6XeDlpCAApAAZhXqUdcAApAAZhXqbeDlhCAAlAA5lXqEReAAlAA5lXq7aAlBKAAFIB5lXrEBaAAFIB5lXo7aAkBKAAFYF6lHnEBKAAFYF6l3g5aQgAKQAGYV6lHXAAKQAGYV6m3g5YQgAJQAOZV6hEXgAJQAOZV6u2gJQSgABSAeZV6xAWgABSAeZV6O2gJASgABWBepR7x1bfjtPN71dVL1toW9+2dXL6XOScC79rr0rvix81bx8uw19zZ/3pM4/q9zHnQwLvsxeN7C184tnGLf3VCcrdfM6NxR+/2+eROO3C/xi2YNjm5H914feN4Tam3g5YQgAJQAOZV6hEXgAJQAOZV6u2gJQSgABSAeZV6xAWgABSAeZV6O2gJASgABWBepR5xASgABWBepd4OWkIACkABmFepR1wACkABmFept4OWEIACUADmVeoRF4ACUADmVertoCUEoAAUgHmVesQFoAAUgHmVejtoidVHfMFFverKH6xz26zYv3G5Y00ACsDSlXrEV9+OObf3qgueXWub3bNnctne5Tdo4F3x41510aPNa/gF9PDmPTe7cSN5l9+5Lxyb3ILnj27cZf92XHI/XDyzccft/oXkFhwwuXEnT9krOQG4fkq9HbSEABSAAjCvUo+4ABSAAjCvUm8HLSEABaAAzKvUIy4ABaAAzKvU20FLCEABKADzKvWIC0ABKADzKvV20BICUAAKwLxKPeICUAAKwLxKvR20hAAUgAIwr1KPuAAUgAIwr1JvByNU1/W0EMIDMcblIYSFVVW9pdvtbhNC+FEI4aEY4xUTJkx4x3p8HQEoAAVgRhv7iI/67RCAAlAAZiEAx6BOp/PhGOOyoaGhzaqqquq6XlTX9dQQwn2dTudjVVVVMcZZdV0f1+9rCUABKADz2phHfIPcDgEoAAVgFgJwDOp2u3WMcbvhj2OM00MI54UQHl/jx2y15sdNho/4+67dtbfNiinrXO4gE4ACcFO2MY/4hrgdm39zh97b79lzrRX5MucLH0vv3F83r+EX0MOb9cyRjev3Mudk4D03O7lj/+Woxi3812OS+/HNhzZu7t7bJ5cKvHn77ZHcHd++tnG8RgCOcd1u909ijE+GEA4OIXx/je96cwjh1X6fLwAFoADMK9cRH63bIQAFoADMQwCOYUNDQ0MhhBV1XR9S1/VH13HE/6Pf1xCAAlAA5pXjiI/m7RCAAlAA5iEAx6gQwgdCCL+o63qfqlr92zaPDX9/XdcTQgiP9Ps6AlAACsC8NvYRH+3bIQAFoADMQwCOQXVdj4sxPtvtdj+z5rfHGJcN//meuq6PCCGcvB5fSwAKQAGY0cY84hvidghAASgA8xCAY1AI4ZgY4yshhKUhhPtW/efRdV2/L4RwZwhhRQjh+i233HLzfl9LAApAAZjXxjziG+J2CEABKADzEICMiAAUgAIwr1KPuAAUgAIwr1JvBy0xfMT/23W7Zg+vtkwACsCNqdQjvvo9gJfvVNa7/FKBd86vk3vTaa82LvmOwIse7c38xRGN6/suvwED79h/OSr58/aLx3u/d3jjTpy8Q3In7TupcXP23DW5H1x3deN4Tam3g5YQgAJQAOZV6hEXgAJQAOZV6u2gJQSgABSAeZV6xAWgABSAeZV6O2gJASgABWBepR5xASgABWBepd4OWkIACkABmFepR1wACkABmFept4OWEIACUADmVeoRF4ACUADmVertoCUEoAAUgHmVesQFoAAUgHmVejtoCQEoAAVgXqUecQEoAAVgXqXeDlpCAApAAZhXqUd8dQBePLmslzknAu/Np76a3Nvn/LJxyZdEn/vr3l5PfbVxI3mZ88ynD09uv6cObdysXxyZ3LI7mnfylJ2TSwXe8ZO+lJwAXD+l3g5aQgAKQAGYV6lHXAAKQAGYV6m3g5YQgAJQAOZV6hEXgAJQAOZV6u2gJQSgABSAeZV6xAWgABSAeZV6O2gJASgABWBepR5xASgABWBepd4OWkIACkABmFepR1wACkABmFept4OWEIACUADmVeoRF4ACUADmVertoCUEoAAUgHmVesRXB+AFM4t6l18q8Dab88vk/vDwpxv35tNeTW7nJw5pXN93+SUCb9+nDk0u9fMe9PPDk/vpPbMbt2Dazskdt3vzjv3yjsndfs03G8drSr0dtIQAFIACMK9Sj7gAFIACMK9SbwctIQAFoADMq9QjLgAFoADMq9TbQUsIQAEoAPMq9YgLQAEoAPMq9XbQEgJQAArAvEo94gJQAArAvEq9HbSEABSAAjCvUo+4ABSAAjCvUm8HLSEABaAAzKvUIy4ABaAAzKvU20FLCEABKADzKvWIC0ABKADzKvV20BICUAAKwLxKPeKrA/CcOUW9zDkVeJsf/nRy75z+cOM2m/Nicqm70u9lzoMG3s5PHJL8eSc9+dXkHl12bOPOmLFLcsfstkPjZu38+eS+d9XljeM1pd4OWkIACkABmFepR1wACkABmFept4OWEIACUADmVeoRF4ACUADmVertoCUEoAAUgHmVesQFoAAUgHmVejtoCQEoAAVgXqUecQEoAAVgXqXeDlpCAApAAZhXqUdcAApAAZhXqbeDEarreloI4YEY4/IQwsKqqt4aY9wjxvhMCGHpqh29Hl9HAApAAZjRxj7io307BKAAFIB5CMAxqNPpfDjGuGxoaGizqqqquq4X1XU9NcZ4bl3Xn34jX0sACkABmNfGPOIb4na8/YyvFfUuv1TgvWvayuT+dNJPGrf54T9P7hMPT2/cSN7l95lHD3auj3cAAAt5SURBVErufz08vXHbPzYjuadWzmnc2Yftllwq8I7a6bPJLbniksbxGgE4BnW73TrGuN3wxzHG6THG+THGZSGE60MI99d1vajT6fznfl9LAApAAZjXxjziG+J2CEABKADzEIBjXLfb/ZMY45N1XX88hHB9jPH9VVVVMcbj67pe1O/zBaAAFIB55Trio3U7BKAAFIB5CMAxbGhoaCiEsKKu60PW8X1bxBhf7Pc1BKAAFIB55Tjio3k7BKAAFIB5CMAxKoTwgRDCL+q63qeqqqqu63HDf11VVTV+/Ph3hRCe6/d1BKAAFIB5bewjPtq3QwAKQAGYhwAcg+q6HhdjfLbb7X5m+Nve8573/KcQwvMTJ078i6qqqhjjUSGEM9fjawlAASgAM9qYR3xD3A4BKAAFYB4CcAwKIRwTY3xl1esa7ht+bUOn0/mbVR8/FEK4esstt9y839cSgAJQAOa1MY/4hrgdAlAACsA8BCAjIgAFoADMq9QjLgAFoADMq9TbQUsIQAEoAPMq9YivDsBTLi/qZc6pwHv3pJ8kt/WOtzduy6kPJvfhh6Y2rt/LnJOBt3Jach98cErjPvnIgck988Tcxi2c9ZXkUoF39G6fT+6WSy9sHK8p9XbQEgJQAArAvEo94gJQAArAvEq9HbSEABSAAjCvUo+4ABSAAjCvUm8HLSEABaAAzKvUIy4ABaAAzKvU20FLCEABKADzKvWIC0ABKADzKvV20BICUAAKwLxKPeICUAAKwLxKvR20hAAUgAIwr1KPuAAUgAIwr1JvBy0hAAWgAMyr1CO+OgDnfqeod/mlAu+9O9yWXP0P1zdu/JfvSi48sF/j+r7Lb8DA++CDU3oTl+/buI/+9IDkXnhmXuMWHTcpuVTgzdnji8kJwPVT6u2gJQSgABSAeZV6xAWgABSAeZV6O2gJASgABWBepR5xASgABWBepd4OWkIACkABmFepR1wACkABmFept4OWEIACUADmVeoRF4ACUADmVertoCUEoAAUgHmVesQFoAAUgHmVejtoCQEoAAVgXqUecQEoAAVgXqXeDlpCAApAAZhXqUdcAApAAZhXqbeDlhCAAlAA5lXqEV8dgMfcXdTLnFOBF/7h+uT+28cvadzQ525Obtz9+zRuJC9zTgXexOX79v74vr0b92cr9k/u5RdPadzFJ+2Z3PGTvtC4uXtvn9xNFy9qHK8p9XbQEgJQAArAvEo94gJQAArAvEq9HbSEABSAAjCvUo+4ABSAAjCvUm8HLSEABaAAzKvUIy4ABaAAzKvU20FLCEABKADzKvWIC0ABKADzKvV20BICUAAKwLxKPeICUAAKwLxKvR20hAAUgAIwr1KPuAAUgAIwr1JvBy0hAAWgAMyr1CM+fDvecdjSot7llwq8bT5+SXIf+suzG/dfP3lVcm+9d8/GjeRdflss3Tu56u7mpaJ03P379H77m9Mbd/mCfZJLBd78/XdM7sZvnNc4XlPq7aAlBKAAFIB5lXrEBaAAFIB5lXo7aAkBKAAFYF6lHnEBKAAFYF6l3g5aQgAKQAGYV6lHXAAKQAGYV6m3g5YQgAJQAOZV6hEXgAJQAOZV6u2gJQSgABSAeZV6xAWgABSAeZV6O2gJASgABWBepR5xASgABWBepd4OWqLb7W4bY+x9euH+ve1vOcpuOaq3//eOadzy5Wcmd/edp26Q/fjOBcndfuu8xp0++6Dkzj1pTuMuOftrA2/JjYuTW7Fiha1Y0bv22mt7McZet9vdNvc9eCOGb8e2R1zV+/jp9621v5l/d3L/35wfNO7vZy1J7h9mLm7cP06/JrnPTr60cV/c/YLkdtzprMZ9bq+Lkvt/bzi0cTveclRy/3TTEY375I2HJfc/Fjfv77+T3n1Lz2jcotMOS+7kw6Y3bsGR6V2+8JzG5f73tU0r9XbQEiGEHWKMPTPLuxDCDrnvwRvhdpi1Y6XdDlpi/Pjx7woh7NDtdret6/pDZrZx1+12tw0h7DB+/Ph35b4Hb4TbYZZ3pd4OAAAAAAAAAAAAAAAAAAAA+qrr+nMhhBUxxodDCIfnfp5hIYRLY4wrQwhLQwhL67r+dK5nqev6j2KMyzudztarPv6fIYT7YowrY4xnVFX1phY803ExxifX+Oe1d4ZnmhZCeCDGuDyEsLCqqrd0u91tQgg/CiE8FGO8YsKECe9owXO9Nca4R4zxmeF/XiGEozf2c5WujbfD3RjoudyO9Xsmd4NNx8SJE/9LjPFnEyZMeGdVVW8OIdwSQvhE7ueqqqoKITwyNDS0Re7nqOv6I6sOwKurDubbYow/q+u6W1VVFUK4IMY4KfMzVTHGW2OM79+Yz7GmTqfz4RjjsqGhoc1WPeOiuq6nhhDu63Q6H1v1jLPquj6uDc8VYzw3ZxyUrq23w914w8/ldryBZ3I32GTUdb3Tql/ZDH+8c4zx3JzPVFVVNWHChHeGEH4VQrg+xrisrusjcj1LCGFht9v96xjjE51OZ+tOp/OxGOOS4e+PMW635sc5nmnVc7wYQrgqxrgshHBKVVVv25jP1O126xjjdsMfxxinhxDOCyE8vsaP2WrNj3M9V4xx/qp/TteHEO6v63pRp9P5zxvzuUrXxtvhbrzx51r1LG7HejyTu8Empa7rGSGE2cMfdzqdvwkhfCfnM1VVVYUQ/jzGePnQ0NAWEyZMeEeMcUkIYdeczxRjfLLT6Wzd7Xa/EGP8xvC313XdjTGuzPlMq/7XGa7rdrtbVVX1phDCRTHGo3I8U1VVVbfb/ZNVv6V0cAjh+2t815tDCK/mfq66rj++KhLeX1VVFWM8vq7rRbmeq0RtvB3uxht/Lrdj/Z/J3WCTEmOc+fojHmO8IeczrUu32/1MCOHKnM8wfDBDCNu//pCHEB7K+Uyv//YQwgdijPfkeKahoaGhEMKKuq4Pqev6o+s44v+R+7nW8X1bxBhfzPFcpSrhdrgb/Z/r9d/udjQ/0zq+z92gXK//bZu6rneq6/rrOZ+pqqqq2+3+ZQjhU8MfhxD+KcZ4Wc5nGj6Yq37r5tbhb+90Oh8LIdyc85m63W5d1/WOw99e1/WHQgh3beznCSF8IITwi7qu96mq1b9t89gazzUhhPBI7ueq63rc8F9X1er/fdvnNvZzlayNt8PdeOPP5Xas/zO5G2xS3vve9747xvhEXdfjqqp6awjhpjb8AddVv/p7dNy4cX9YVdXbQgg3d7vdL+R8pjV+xfy2EMJTnU4nVtXv/0xNCGFKzmda9d8wPL/VVlu9p6qqPwghnF/X9YyN+Sx1XY+LMT7b7XY/87pnXDb8Z2nquj4ihHBy7ud6z3ve859CCM9PnDjxL1Y941EhhDM35nOVro23w91448/ldqz/M7kbbHJCCP8cY1wef/9qgjm5n2dYXdcHhN//v/8/HEI4JvfzvO4PTW+36hUAD4UQLqiq6i25nymEsH34/Ss5VsYYz9nYzxRCOCbG+Mqqfy73rfrPo+u6fl8I4c4QwooQwvVbbrnl5m14rlV/Zu2+Vf83vHpjP9emoI23w91448/ldqz/M7kbAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAACwwfz/Yeb1+eSjBVcAAAAASUVORK5CYII=\">"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "showmaps(*cl1_maps)\n",
    "showmaps(*cl2_maps)\n",
    "showmaps(*cl3_maps)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "metadata": {
    "collapsed": false,
    "scrolled": true
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "RNAfold\n",
      "Vfold\n"
     ]
    },
    {
     "ename": "KeyboardInterrupt",
     "evalue": "",
     "output_type": "error",
     "traceback": [
      "\u001b[0;31m---------------------------------------------------------------------------\u001b[0m",
      "\u001b[0;31mKeyboardInterrupt\u001b[0m                         Traceback (most recent call last)",
      "\u001b[0;32m<ipython-input-21-42fd4f24caa4>\u001b[0m in \u001b[0;36m<module>\u001b[0;34m()\u001b[0m\n\u001b[1;32m      3\u001b[0m \u001b[0;32mfor\u001b[0m \u001b[0mfp\u001b[0m \u001b[0;32min\u001b[0m \u001b[0mfprogs\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m      4\u001b[0m     \u001b[0mprint\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mfp\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m----> 5\u001b[0;31m     \u001b[0mfolds_cut30\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mappend\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mcutmap\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mseq_cl1\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;36m30\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mfp\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m",
      "\u001b[0;32m<ipython-input-17-5d352a087bd4>\u001b[0m in \u001b[0;36mcutmap\u001b[0;34m(seq, cut, prog)\u001b[0m\n\u001b[1;32m      9\u001b[0m             \u001b[0mj\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mlen\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mseq\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m-\u001b[0m\u001b[0mji\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m     10\u001b[0m             \u001b[0;31m#print(seq[i:j])\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m---> 11\u001b[0;31m             \u001b[0mp\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mcpred\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mrun\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mseq\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0mi\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0mj\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0moutdir\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;32mNone\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mrmdir\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;32mTrue\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m     12\u001b[0m             \u001b[0mcutmap\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0mi\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0mji\u001b[0m\u001b[0;34m]\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mp\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mix\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0;36m0\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0;36m1\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m     13\u001b[0m             \u001b[0mcutstruct\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0mi\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0mji\u001b[0m\u001b[0;34m]\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mp\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mix\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0;36m0\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0;36m0\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m/scratch/tc427/lucy/consensus/code/consensus.py\u001b[0m in \u001b[0;36mrun\u001b[0;34m(self, seq, outdir, rmdir, nodir)\u001b[0m\n\u001b[1;32m    128\u001b[0m                 \u001b[0mrmdir\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0;32mTrue\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    129\u001b[0m             \u001b[0moutdir\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0msetup_run\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mseq\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0moutdir\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 130\u001b[0;31m             \u001b[0mpreds\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mpredict\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mseq\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m    131\u001b[0m             \u001b[0;32mif\u001b[0m \u001b[0;32mnot\u001b[0m \u001b[0mrmdir\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    132\u001b[0m                 \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mplot_structs\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mpreds\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mseq\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m/scratch/tc427/lucy/consensus/code/consensus.py\u001b[0m in \u001b[0;36mpredict\u001b[0;34m(self, seq)\u001b[0m\n\u001b[1;32m     63\u001b[0m         \u001b[0;32mfor\u001b[0m \u001b[0mbinpath\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mcall\u001b[0m \u001b[0;32min\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mprogs\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m     64\u001b[0m             \u001b[0;32mif\u001b[0m \u001b[0mbinpath\u001b[0m \u001b[0;32mis\u001b[0m \u001b[0;32mnot\u001b[0m \u001b[0;32mNone\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m---> 65\u001b[0;31m                 \u001b[0mpred_data\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mpd\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mmerge\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mpred_data\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mcall\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mbinpath\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mseq\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mmaindir\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mon\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;34m\"struct\"\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mhow\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;34m\"outer\"\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m     66\u001b[0m         \u001b[0;31m#pred_data.fillna(0, inplace=True)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m     67\u001b[0m         \u001b[0;31m# create a global rank, sum of the ranks from each method,\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m/scratch/tc427/lucy/consensus/code/pred.py\u001b[0m in \u001b[0;36mcall_Vfold\u001b[0;34m(binpath, seq, outdir, Npreds)\u001b[0m\n\u001b[1;32m     85\u001b[0m     \u001b[0mos\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mchdir\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mos\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mpath\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mdirname\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mbinpath\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m     86\u001b[0m     \u001b[0;31m# get base pair probabilities from Vfold2D\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m---> 87\u001b[0;31m     \u001b[0msubprocess\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mcall\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0;34m\"./\"\u001b[0m\u001b[0;34m+\u001b[0m\u001b[0mos\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mpath\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mbasename\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mbinpath\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m\"0\"\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m\"37\"\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mos\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mpath\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mjoin\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0moutdir\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m\"seq\"\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mos\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mpath\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mjoin\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0moutdir\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m\"seq.pij\"\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m\"0\"\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mshell\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;32mFalse\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m     88\u001b[0m     \u001b[0;31m##print Vout\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m     89\u001b[0m     \u001b[0mos\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mchdir\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0moutdir\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m/usr/lib/python3.5/subprocess.py\u001b[0m in \u001b[0;36mcall\u001b[0;34m(timeout, *popenargs, **kwargs)\u001b[0m\n\u001b[1;32m    557\u001b[0m     \u001b[0;32mwith\u001b[0m \u001b[0mPopen\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m*\u001b[0m\u001b[0mpopenargs\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m**\u001b[0m\u001b[0mkwargs\u001b[0m\u001b[0;34m)\u001b[0m \u001b[0;32mas\u001b[0m \u001b[0mp\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    558\u001b[0m         \u001b[0;32mtry\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 559\u001b[0;31m             \u001b[0;32mreturn\u001b[0m \u001b[0mp\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mwait\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mtimeout\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mtimeout\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m    560\u001b[0m         \u001b[0;32mexcept\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    561\u001b[0m             \u001b[0mp\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mkill\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m/usr/lib/python3.5/subprocess.py\u001b[0m in \u001b[0;36mwait\u001b[0;34m(self, timeout, endtime)\u001b[0m\n\u001b[1;32m   1656\u001b[0m                         \u001b[0;32mif\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mreturncode\u001b[0m \u001b[0;32mis\u001b[0m \u001b[0;32mnot\u001b[0m \u001b[0;32mNone\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m   1657\u001b[0m                             \u001b[0;32mbreak\u001b[0m  \u001b[0;31m# Another thread waited.\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m-> 1658\u001b[0;31m                         \u001b[0;34m(\u001b[0m\u001b[0mpid\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0msts\u001b[0m\u001b[0;34m)\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_try_wait\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;36m0\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m   1659\u001b[0m                         \u001b[0;31m# Check the pid and loop as waitpid has been known to\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m   1660\u001b[0m                         \u001b[0;31m# return 0 even without WNOHANG in odd situations.\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m/usr/lib/python3.5/subprocess.py\u001b[0m in \u001b[0;36m_try_wait\u001b[0;34m(self, wait_flags)\u001b[0m\n\u001b[1;32m   1606\u001b[0m             \u001b[0;34m\"\"\"All callers to this function MUST hold self._waitpid_lock.\"\"\"\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m   1607\u001b[0m             \u001b[0;32mtry\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m-> 1608\u001b[0;31m                 \u001b[0;34m(\u001b[0m\u001b[0mpid\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0msts\u001b[0m\u001b[0;34m)\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mos\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mwaitpid\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mpid\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mwait_flags\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m   1609\u001b[0m             \u001b[0;32mexcept\u001b[0m \u001b[0mChildProcessError\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m   1610\u001b[0m                 \u001b[0;31m# This happens if SIGCLD is set to be ignored or waiting\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;31mKeyboardInterrupt\u001b[0m: "
     ]
    }
   ],
   "source": [
    "fprogs = [\"RNAfold\", \"Vfold\"]\n",
    "folds_cut30 = list()\n",
    "for fp in fprogs:\n",
    "    print(fp)\n",
    "    folds_cut30.append(cutmap(seq_cl1, 30, fp))"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.5.2"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 1
}