2016-09-20 Exp preliminar.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Experimento preliminar: últimas pruebas"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "import numpy as np\n",
    "import seaborn, math, os, random\n",
    "from scipy.integrate import odeint\n",
    "from matplotlib import pyplot as plt\n",
    "from copy import deepcopy"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "%run modelo.py\n",
    "%run medidas_biodiversidad.py"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "files = []\n",
    "with open('iteraciones.txt', 'r') as file:\n",
    "    rows = file.readlines()\n",
    "    for row in rows:\n",
    "        files.append([row[:20], int(row[-2])])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "all_condiciones_iniciales = [np.loadtxt(\"./comunidades/\" + file[0][:20] + \"-c.txt\") for file in files]\n",
    "all_tasas_reproduccion = [np.loadtxt(\"./comunidades/\" + file[0][:20] + \"-v.txt\") for file in files]\n",
    "all_matriz_interacciones = [np.loadtxt(\"./comunidades/\" + file[0][:20] + '-m.txt') for file in files]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "def modificar_matriz(m):\n",
    "    for i in range(10):\n",
    "        m[i][i] *= 15\n",
    "    return m\n",
    "\n",
    "def modificar_tasas(tasa):\n",
    "    tasa[tasa > 0] = tasa[tasa > 0] * 20\n",
    "    return tasa\n",
    "\n",
    "def modificar_condiciones0(c0):\n",
    "    c0 = c0 * 10\n",
    "    return c0"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "matriz_interacciones = deepcopy(all_matriz_interacciones)\n",
    "matriz_interacciones = [modificar_matriz(m) for m in matriz_interacciones]\n",
    "\n",
    "tasas_reproduccion = deepcopy(all_tasas_reproduccion)\n",
    "tasas_reproduccion = [modificar_tasas(tasa) for tasa in tasas_reproduccion]\n",
    "\n",
    "condiciones_iniciales = deepcopy(all_condiciones_iniciales)\n",
    "condiciones_iniciales = [modificar_condiciones0(c) for c in condiciones_iniciales]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {
    "collapsed": false,
    "scrolled": true
   },
   "outputs": [],
   "source": [
    "random.seed(17)\n",
    "comunidades = random.sample(range(201), 100)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "[133, 106, 77, 93]"
      ]
     },
     "execution_count": 8,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "comunidades = comunidades[:4]\n",
    "comunidades"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Experimento: comunidades originales\n",
    "Parámetros"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "paisajes = np.load(\"./paisajes_h2.npy\")\n",
    "\n",
    "Dispersion = {\"b\": 0.3,\n",
    "              \"m\": 1,\n",
    "              \"i\": 1}\n",
    "\n",
    "Mortalidad = {\"b\": 0,\n",
    "              \"m\": 0.3,\n",
    "              \"i\": 0.6}\n",
    "t_total = 100"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "for c in comunidades:\n",
    "    res = [None] * len(paisajes)\n",
    "    for p in range(len(paisajes)):\n",
    "        res[p] = matriz_agroecologica(paisajes[p].tolist(),\n",
    "                                      matriz_interacciones[c], \n",
    "                                      tasas_reproduccion[c], \n",
    "                                      condiciones_iniciales[c],\n",
    "                                      t_total, \n",
    "                                      Dispersion,\n",
    "                                      Mortalidad)\n",
    "        #res[p] = res[p][-3:]\n",
    "        \n",
    "    res = np.array(res)\n",
    "    np.save(\"./datos/exp_prelim1_c{:03d}\".format(c), res)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Análisis de datos"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "poblacion = np.load('./datos/exp1_h0_c077.npy')\n",
    "poblacion = poblacion[0]\n",
    "paisaje = paisajes[0].tolist()\n",
    "\n",
    "biomasa_min = 5\n",
    "\n",
    "# Biomasa paisaje\n",
    "biomasa = [None] * len(poblacion)\n",
    "riqueza = [None] * len(poblacion)\n",
    "\n",
    "for i in range(len(poblacion)):\n",
    "    biomasa[i] = riqueza_paisaje(poblacion, i, biomasa_min)[0] \n",
    "    riqueza[i] = riqueza_paisaje(poblacion, i, biomasa_min)[1]\n",
    "\n",
    "riqueza_milpa = [None] * len(poblacion)\n",
    "biomasa_milpa = [None] * len(poblacion)\n",
    "\n",
    "for i in range(len(poblacion)):\n",
    "    riqueza_milpa[i] = riqueza_agricola(poblacion, paisaje, i, biomasa_min)[1]   \n",
    "    biomasa_milpa[i] = riqueza_agricola(poblacion, paisaje, i, biomasa_min)[0]   "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "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\", \"pdf\", \"png\", \"ps\", \"raw\", \"svg\"];\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,iVBORw0KGgoAAAANSUhEUgAAAoAAAAHgCAYAAAA10dzkAAAgAElEQVR4nOy9e3Bc9ZXv+zMDwxy4mVMXCAJLoFbv31qVgRQnVIXUDOcmtxKok3PriphAIIBDeISER5Jxpo1fgME8jAEb846Jgw0Y8zAYaMtvW5ZsS5YsWX6/JEuW9Va3JdsSM5OZOecO6/7RrZYwtlv27tbuvdf3U/WpQt0t9RZf91rLWq2fjQEAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAOQq1trJzFxLRFXMvMAYc4619lZmriSiciIqdxzn28YYk/x4OzOXMXOZtXZ68mucS0QLk59TY62909NvCgAAAAAAnBjHca4loh3GmLOMMYaIPmHm3zBz70UXXfQNY4xh5l8y88rk/eXM/P3jvw4RTbDWzjXGmPz8/AuJqNVxnItH8FsBAAAAAADDZFReXt75Ax8Q0RwimsjMzcwcTt42jojeSf53ORH94Pgvkvwp4Y+GfLzQWjs2+5cPAAAAAADOGMdxLBG1hcPhy8Ph8HXMfJSI9hJRUygUusSY1AC4hIjWEtHqcDj8XWOMYeYGa+0VA1+LiGZaayd59b0AAAAAAIA0FBUVXUVEjeFw+Lq8vLzzmbk5FAp9yxhjrLVjiWh58r/HDLwf0HGcHzJzpzHmrOMHQGvtLCKa6Mk3AwAAAAAATo219mpmbgiHw9cYY0w4HL6GiLYOuf9vmfmfT/S5RBS/7LLLRjNzKTNfP3A7My9i5tvSPfeXX34pAAAAAPAXmZg/gIeMHj36PGaut9ZeOeS2i4gobq39pjHGWGt/zMy1JvGTvgpmzjfGGCL6DhF1GGNGJd8nON8YY0Kh0CVE1FZQUHBBuuf/y7//b+nt/UJ6emDQ7e39QkQEeSsReesSeeuyt/cLDIB+x1p7HzP3MHNZ8jd8y5h5ChHdTETVyY9Lh6x9b2TmuuRjK621/5D8UucQ0Xxm3kRE1dbaMcN5/p8/skzWVjXL4cNfwIDb05NoED093l8LRN4QeUN3eWdtMAE6KI5EpTgSlTcXb5dYrN/zP9QwuwUDDUKPyFuXyFuXGACBa1ZVH5IbJyyR4khUJr26QQ629Hr+Bxtmr2CgQegReesSeesSAyBwjYhIzc4OuWvaSimOROWuaSulZmeH53+4YXYKBhqEHpG3LpG3LjEAAtcMFIyDLb0y6bWNUhyJyo0TlsjHa/ZLPI6VcJBEg9Al8tYl8tYlBkDgmqEFIxbrlzcX75CB9wXOeHuzdHQe8/wPOsxcwUCD0CPy1iXy1iUGQOCaExWMNVXNcsuUZVIcicpDz5XKngNxz/+ww8wUDDQIPSJvXSJvXWIABK45WcHY0xCTB58rleJIVG6dskxKq3FUjN9Fg9Al8tYl8tYlBkDgmlMVjPbOYzJj/ubUSnjupztwVIyPRYPQJfLWJfLWJQZA4Jp0BSMe75ePV++XMQ8njoqZ/NpGaW494vkffnhmBQMNQo/IW5fIW5cYAIFrhlswqne0y51PDBwVs0pqd+GoGL+JBqFL5K1L5K1LDIDANadTMJpaemXiqxuSR8WUyOK19Z6/CODpFQw0CD0ib10ib11iAASuOd2CEYv1y5xPtqfeF/j8OzXS2YWjYvwgGoQukbcukbcuMQAC15xpwVhdeVBumbJUiiNR+e3z62QfjorJedEgdIm8dYm8dYkBELjGTcHYXR+TB2asleJIVH7+yDJZt/mQ5y8KeOqCgQahR+StS+StSwyAwDVuC0Z7xzGZPq9aiiNRuSESlbc+2ykx/BNyOSkahC6Rty6Rty4xAALXZKJgxOP98uHKfamjYh55owJHxeSgaBC6RN66RN66xAAIXJPJglG9ffComLufXCVbdnV6/iKBXy0YaBB6RN66RN66xAAIXJPpgtF0qFcefnm9FEei8tOJJfJpKY6KyRXRIHSJvHWJvHWJARC4JhsFozvWJ298vC11VMzMBbU4KiYHRIPQJfLWJfLWJQZA4JpsFoyVFQflZ5MTR8X87oV1sq8RR8V4XTDQIPSIvHWJvHWJARC4JtsFY1d9TO5/dvComLIaHBXjZcFAg9Aj8tYl8tYlBkDgmpEoGO0dx+TptwaPipn3OY6K8apgoEHoEXnrEnnrEgMgcM1IFYx4vF8+WLE3dVTMo29UyKE2HBUz0gUDDUKPyFuXyFuXGACBa0a6YGza1ia/eHyFFEeics9Tq6RuN46KGcmCgQahR+StS+StSwyAwDVeFIzG5h4Z/9LgUTGfr2vw/MWkQTQIXSJvXSJvXWIABK7xqmB0x/rk9Y8Gj4p58b0t0tnV5/mLKsiiQegSeesSeesSAyBwjdcFY8XGJrk5eVTM72eWSX3TYc9fWEEVDUKXyFuXyFuXGACBa3KhYOzc3y2/nr5GiiNRue3R5VJe2+L5iyuIokHoEnnrEnnrEgMgcE2uFIy2jqPy5NyqxFEx46Py9pJdOComCwUjV/KGyBsib+gub6/nB+BzcqlgxOP9snD54FExj/2xUlpwVExGC0Yu5Q2RN0Te8Mzz9np+AD4nFwtG5dY2GTs1cVTMvU+vlq17ujy/piCIBqFL5K1L5K1LDIDANblaMA4090hkdrkUR6Jy08QSWVJ2wPNr8rtoELpE3rpE3rrEAAhck8sFo6u7T179cGvqqJjZC7dIVzeOinFTMHI5b4i8IfKGw8/b6/kB+Bw/FIxlGxrl5kmJo2LGzcJRMW4Khh/yhsgbIm+YPm+v5wfgc/xSMLbvGzwq5vbHlsuGLa2eX5PfRIPQJfLWJfLWJQbAgGCtnczMtURUxcwLjDHnWGtvZeZKIiononLHcb5tjDGFhYVFzFxGRBuZucJae3Xya5xLRAuTn1Njrb1zOM/tp4LR2n5Upv1pU+qomHdLdkscR8WcVsHwU94QeUPkDU+edxbHEjASOI5zLRHtMMacZYwxRPQJM/+GmXsvuuiibxhjDDP/kplXJu9fTkS3J2//PjPvTN4+0Vo71xhj8vPzLySiVsdxLk73/H4rGPF4vyxYukd+Mj5xVMzjb26S1vajnl+XH0SD0CXy1iXy1iUGwGAwKi8v7/yBD4hoDhFNZOZmZg4nbxtHRO8YY84mor8YY84Z8vh2a21B8qeEPxpy+0Jr7dh0T+7XglFR1yp3JI+K+dUza2T7XhwVM5yC4de8IfKGyBt+Ne9MDiLAYxzHsUTUFg6HLw+Hw9cx81Ei2ktETaFQ6JLCwsJLiejw0M9h5lpr7d8zc4O19oqB24loprV2Urrn9HPBaDjYI394MXlUzKQSKVmPo2LSFQw/5w2RN0TecDDvbMwhwAOKioquIqLGcDh8XV5e3vnM3BwKhb5ljDHW2rFEtLywsPBSZu4Z+nnMXMfM3zt+ALTWziKiiemeV0Sktzfxh8mPdsf65OUP6lJHxbz8QZ10x/o8v65ctLc30SD8nDdE3hB5w0TemZ9EwIhjrb2amRvC4fA1xhgTDoevIaKtQ+7/W2b+Z2PMXxHRvxpj/nrgPiKKhUKhS5i5lJmvH7idmRcx823pnlsCwvqt7XLz5MRRMRNe3Si9fX/x+pIAAACArJHZSQSMOKNHjz6PmeuttVcOue0iIopba79pjDHW2h8zc60xxhDRZ8x818DtRLQ5efs4IppvjDGhUOgSImorKCi4IN3ziwTnb4w79nXJfc8kjoq5Y+oKqdja6vk15ZL4CYEukbcukbcu8RPAAGCtvY+Ze5JHu5QzcxkzTyGim4moOvlx6cAxMMlf+FhLRBuJqHxgTWyMOYeI5jPzJiKqttaOGc7ziwTrPSOt7Ufl8TcTR8X8ZPwSWbB0D46KGfKekaDlDZE3RN4a7enBAAhcEsSCEY/3y7tLd8sN4xPvC5z2JxwVM1Awgpg3RN4QeWsTAyBwTZALxoYtrXL7Y8ulOBKV+55ZI9v3dXt+TV4XjCDnDZG3ZpG3LjEAAtcEvWDUNx2WcbPKpDgSlZsnLZVlGxo9vyYvC0bQ84bIW6vIW5cYAIFrNBSMru4+eWnh4FExr3ywVbq6+zy/Li8Khoa8IfLWKPLWJQZA4BpNBWNJ+QG5aWKJFEei8k+zy+VAc4/n1zTSBUNT3tpF3rpE3rrEAAhco61gbNvTJfc+vVqKI1EZO3WFVG5t8/yaRrJgaMtbs8hbl8hblxgAgWs0FoyWtqMydU5l6qiYhct0HBWDBqFL5K1L5K1LDIDANVoLRizeL28v2ZU6KubJuVXS1hHso2LQIHSJvHWJvHWJARC4RnvBWF/bIrc9mjgq5tfT18jO/cE9KgYNQpfIW5fIW5cYAIFrUDASR8X8fmbyqJjJS2XFxibPrylbBQN56xF56xJ56xIDIHANCkbCzq4+efG9LamjYl77aKt0x4J1VAwahC6Rty6Rty4xAALXoGB81WhZQ+qomPEvrZfGAB0VgwahS+StS+StSwyAwDUoGF+3bnen3PPUqsRRMY+vkE3bgnFUDBqELpG3LpG3LjEAAtegYJzYQ21H5NE3KqQ4EpUxDy+R91fs9f1RMWgQukTeukTeusQACFyDgnFyY/F+mRfdJTck3xf49FvVvj4qBg1Cl8hbl8hblxgAgWtQMNJbVnModVTM/c+ulV31Mc+v6UwLBvLWI/LWJfLWJQZA4BoUjOG5rzEuv3thnRRHovKzyUtlZcVBz6/pTAoG8tYj8tYl8tYlBkDgGhSM4dvZdUxmLqhNHRXzxqJtvjoqBg1Cl8hbl8hblxgAgWtQME7fT0vr5afJo2Iefnm9NB3q9fyahlswkLcekbcukbcuMQAC16BgnJlbdnXK3U8mjoq584mVUrW93fNrGk7BQN56RN66RN66xAAIXIOCceYeajsijww5KubDlfty+qgYNAhdIm9dIm9dYgAErkHBcGcs3i9vfbYzdVTMM/Oqpb3jmOfXdbKCgbz1iLx1ibx1iQEQuAYFIzOu23xIfv7IMimOROWBGWtldw4eFYMGoUvkrUvkrUsMgMA1KBiZc9+BuPz2+cRRMbdMWSqrK3PrqBg0CF0ib10ib11iAASuQcHIrJ1dx+T5d2pSR8XM+WS7xGK58b5ANAhdIm9dIm9dYgAErkHByI6L19bLjRMSR8VMeGWDNLV4f1QMGoQukbcukbcuMQAC16BgZM/aXR1y17TBo2I27+jwvGAgbz0ib10ib11iAASuQcHIrs2tR2Ty6xtTR8UsWr3fs6Ni0CB0ibx1ibx1iQEQuAYFI/vGYv0y99MdqfcFPjt/s7R3jvxRMWgQukTeukTeusQACFyDgjFyllY3y61TBo6KKZU9DSN7VAwahC6Rty6Rty4xAALXoGCMrHsOxOWh50qTR8UskzVVzSNaMJC3HpG3LpG3LjEAAtegYIy8HZ3HZMbbm1Mr4TcXj8xRMWgQukTeukTeusQACFyDguGN8Xi/fLxmv9w4YYkUR6Iy6dUNcjDLR8WgQegSeesSeesSAyBwDQqGt9bs7JC7pq2U4khU7pq2Ump2Zu+oGDQIXSJvXSJvXWIADAjW2snMXEtEVcy8wHGcy4ionJnLmLmMiKqZud4YY4hoPRFtH7jPWjs9+TXOJaKFzFxJRDXW2juH89woGN57sKVXJr26QYojUblxwhL5eE12jopBg9Al8tYl8tYlBsAA4DjOtUS0wxhzljHGENFiIrpn6GOI6Clr7YPJ/y5n5u8f/3WIaIK1dq4xxuTn519IRK2O41yc7vlRMHLDWKxf3lw8eFTMc2/XSEeGj4pBg9Al8tYl8tYlBsBgMCovL+/8gQ+IaA4RTRj4OBQKhYhoqzFmVPL+ciL6wfFfhIjKHcf50ZCPF1prx6Z7chSM3HLNpma5JXlUzEPPl8qeA/GMFgzkrUfkrUvkrUsMgAHDcRxLRG3hcPjygduIaM7QQS45AC4horVEtDocDn/XGGOYucFae8WQx8201k5K95woGLnnnoaYPJg8KubnjyyT0urMHBWDBqFL5K1L5K1LDIABoqio6CoiagyHw9cN3HbRRRd9g4haTXI9bIwx1toxjuN82xhjHMf5ITN3GmPOOn4AtNbOIqKJ6Z5XRKS3N/GHCeaOHV3H5Nn5iaNibohEZe6nOyQe73f1NXt7Ew0CeesQeesSeeuytxcDYCCw1l7NzA3hcPiaobcz821E9PapPpeI4pdddtloZi5l5uuHfO4iZr4t3XMLyFm+/PJLWbGpOXVUzONzq6T/X/7D68sCAACQA5z51AFygtGjR5/HzPXW2iuPv89a+5q19qEhN53FzBXMnG+MMUT0HSLqMMaMIqJxRDTfGGNCodAlRNRWUFBwQbrnF8HfGHPdzTva5c4nkkfFPLlKtuzqOOO/MSJvPSJvXSJvXeIngAHAWnsfM/ckj3sZOPplijHGMPPnRHTT0Mc7jnMjM9clH1tprf2H5F3nENF8Zt5ERNXW2jHDeX4RvGfEDza19MqEVwaOiimRxWvrz+g9I8hbj8hbl8hblz09GACBS1Aw/GMs1i9zPtmeOirm+XdqpLNr+EfFoEHoEnnrEnnrEgMgcA0Khv9cXXlQbpmyVIojUfntC+tk3zCPikGD0CXy1iXy1iUGQOAaFAx/urs+Jg/MWJs6Kmbd5kPDKhjIW4/IW5fIW5cYAIFrUDD8a3vHMXlmXnXqqJi3PtspsVP8E3JoELpE3rpE3rrEAAhcg4Lhb+Pxfvlw5T4Z83DiqJhH3qiQQ21HTlowkLcekbcukbcuMQAC16BgBMOq7YNHxdz95CrZsqvzhAUDeesReesSeesSAyBwDQpGcGw61CsPv7xeiiNR+enEEvmstOFrBQN56xF56xJ56xIDIHANCkaw7I71yRuLtqWOipm5oFY6u/pSBQN56xF56xJ56xIDIHANCkYwXVnRJD+bnDgq5vcz18n+xsNoEMpE3rpE3rrEAAhcg4IRXHfu75b7n00cFXPbo8tlfW0LGoQiMRDoEnnrEgMgcA0KRrBt6zgqT781eFTMR2vrJX745EfFwOCIgUCXyFuXGACBa1Awgm883i/vr9ibOirmsT9WSstJjoqBwREDgS6Rty4xAALXoGDosWp7m/xyWuKomHufWi11u79+VAwMjhgIdIm8dYkBELgGBUOPPT1fyNH+f5PxyaNibppYIp+va/D8umD28sbrW4/IW5cYAIFrUDD0ONAgYvE+ef2jwaNiXnxvS+qoGBgcMRDoEnnrEgMgcA0Khh6PbxArNjbJzamjYsqkvumw59cIs5c3DLbIW5cYAIFrUDD0eKIGsXN/t/x6+pqvHBXj9XXC7OUNgyvy1iUGQOAaFAw9nqxBtHUclSfnViWOihkflbeX7JJ4HEfF+F0MBLpE3rrEAAhcg4Khx1M1iHi8XxYu2yM/GZ84KmbqnEppaTvq+TXD7OQNgyfy1iUGQOAaFAw9DqdBVG5tk7FTVySOinl6tWzb0+X5dcPs5Q2DI/LWJQZA4BoUDD0Ot0EcaO6Rf5pdnjoqZkn5Ac+vHWYvbxgMkbcuMQAC16Bg6PF0GkRXd5+8+uHW1FExLy2sk65uHBXjJzEQ6BJ56xIDIHANCoYez6RBLNvQKDdPShwVM+5FHBXjJzEQ6BJ56xIDIHANCoYez7RBbN/XLfc9kzgq5o7HlsuGLa2efy8we3lDf4q8dYkBELgGBUOPbhpEa/tReeJPm1JHxbxbshtHxeS4GAh0ibx1iQEQuAYFQ49uG0Q83i8Llg4eFfPEm5uktR1HxeSqGAh0ibx1iQEQuAYFQ4+ZahAVda1yx2PLpTgSlV89s0a278VRMbkoBgJdIm9dYgAErkHB0GMmG0TDwR75w4vJo2ImlcjS9Y2ef38we3nD3Bd56xIDIHANCoYeM90gurr75OX361JHxbz8AY6KySUxEOgSeesSAyBwDQqGHrPVIErWH5CbJpVIcSQqf5hdLg0Hezz/XiEGAm0ib11iAASuQcHQYzYbxLa9XfKrp1cnjoqZukIq6nBUjNdiINAl8tYlBkDgGhQMPWa7QbS0HZXH30wcFfOT8UtkwdI9OComwHnD3BJ56xIDIHANCoYeR6JBxOP98m7JbrlhfOJ9gdP+hKNigpw3zB2Rty4xAALXoGDocSQbxIbaFrk9eVTMr6evkR37uj3//rWJgUCXyFuXGACBa1Aw9DjSDaK+6bCMm1UmxZGo3DxpqSzbgKNigpw3RN5wZPP2en4AGcBaO5mZa4moipkXOI5zGRGVM3MZM5cRUTUz1xtjTGFhYVHyto3MXGGtvTr5Nc4looXMXElENdbaO4fz3CgYevSiQXR198nshVtSR8W8+uFWHBUT4Lwh8oYjl3c25xIwAjiOcy0R7TDGnGWMMUS0mIjuGfoYInrKWvtg8r+XE9HtxhjDzN9n5p3J2ydaa+caY0x+fv6FRNTqOM7F6Z4fBUOPXjaIJWUH5KaJiaNiIrPL5UAzjooJct4QecPs553hcQR4wKi8vLzzBz4gojlENGHg41AoFCKircaYUcaYs4noL8aYc4Y8vt1aW0BE5Y7j/GjI7QuttWPTPTkKhh69bhBb93TJvcmjYsZOXSGVW9s8/38SZL3OGyJvmN28MzaFAO9xHMcSUVs4HL584DYimjMwyBUWFl5KRIeHfg4z11pr/56ZG6y1Vwz5vJnW2knpnhMFQ4+50CBa2o7I1DmVUhyJypiHl8jC5XtxVEyA84bIG2Yv70zOH8BDioqKriKixnA4fN3AbRdddNE3iKjVJNfDJxkA65j5e8cPgNbaWUQ0Md3zioj09ib+MMFg29ubaBBe5x0/3C9vL9mVOirmqT9XSXvnUc///wTNXMkbIm+YnbwzOIIAr7DWXs3MDeFw+JqhtzPzbUT09pCb/oqZ/8UY89cDNxBRLBQKXcLMpcx8/ZDPXcTMt6V7bgHAI7Y3xOWOqSukOBKVB59fJ22xL7y+JAAA8A2ZmD+Ah4wePfo8Zq631l55/H3W2testQ8NvY2IPmPmu5L3/5iINidvH0dE840xJhQKXUJEbQUFBReke34R/I1Ri7n4E4L6g4flH2cmjor52eSlsrKiyfNrCoq5mDdE3jBzeWdmCgGeYa29j5l7kke7DBz9MsUYY5j5cyK66bjHFxDRWiLaSETloVDoW8m7ziGi+cy8iYiqrbVjhvP8InjPiBZ7enLzPUKdXX0y673Bo2Je/2ibdMdwVExQ84bIG2Ym7wyPI0AbKBh6zPUG8fm6Bvlp8qiYh19eL42Hej2/Jj+b63lD5A3d5e31/AB8DgqGHv3QIOp2d8o9T62S4khUfvH4Cqna1u75NflVP+QNkTc887y9nh+Az0HB0KNfGsShtiPy6BsVqaNiPliBo2KCnDdE3vDM8vZ6fgA+BwVDj35qELF4v8z7fKfckHxf4NNvVUt7xzHPr8tP+ilviLzh6eft9fwAfA4Khh792CDKag7Jzx9ZJsWRqNz/7FrZVR/z/Jr8oh/zhsgbDj9vr+cH4HNQMPTo1waxrzEuv3th3ZCjYg56fk1+0K95Q+QNh5e31/MD8DkoGHr0c4Po7DomMxfUpo6KeeNjHBUT5Lwh8obp8/Z6fgA+BwVDj0FoEJ+W1suNEwaPimnCUTGBzhsib3jyvL2eH4DPQcHQY1AaxJZdnXL3k4mjYu58YqVUb8dRMUHOGyJveOK8vZ4fgM9BwdBjkBpEc+sRmfL64FExH67ch6NiApw3RN7w63l7PT8An4OCocegNYhYvF/+/NngUTHT51VLeyeOiglq3hB5w6/m7fX8AHwOCoYeg9og1m0ePCrmgRlrZXcDjooJct4QeUMMgCADoGDoMcgNYu+BuPz2+cRRMbdMWSarN+GomCDnDZG3djEAAtegYOgx6A2is+uYPP9OTeqomDmfbJdYTO/7AoOeN0TemsUACFyDgqFHLQ1i8drBo2ImvrpBmlp0HhWjJW+IvDWKARC4BgVDj5oaRO2uDrlrWuKomF9OWymbd3R4fk3IGyJvmMm8vZ4fgM9BwdCjtgbR3HpEJr+2UYojUblxwhJZtHq/qqNitOWtXeStSwyAwDUoGHrU2CBisX6Z++mO1PsCZ8zfrOaoGI15axZ56xIDIHANCoYeNTeItVXNcuuUxFExDz5XKnsUHBWjOW+NIm9dYgAErkHB0KP2BrHnQFweeq5UiiNRuXXKMllT1ez5NSFviLzhmebt9fwAfA4Khh7RIL6Qjs5jMmP+5tRK+M3FOwJ7VAzy1iXy1iUGQOAaFAw9okEkjMf75ePV++XGCUukOBKVSa9tlObWI55fF/KGyBueTt5ezw/A56Bg6BEN4qvW7OyQX05bKcWRqNw1baXU7AzWUTHIW5fIW5cYAIFrUDD0iAbxdQ+29MqkVzekjor5ZE1wjopB3rpE3rrEAAhcg4KhRzSIExuL9cubi7en3hf43Ds10hGAo2KQty6Rty4xAALXoGDoEQ3i1K7Z1Cy3JI+Keej5Utl7IO75NSFviLzhyfL2en4APgcFQ49oEOnd3RCTB2Ykjor5+SPLpLTav0fFIG9dIm9dYgAErkHB0CMaxPBs7zwm0+cljoq5IRKVuZ/ukJgP3xeIvHWJvHWJARC4BgVDj2gQwzce75dFq/bJmIcTR8VMft1/R8Ugb10ib11iAASuQcHQIxrE6Vu9o13ufCJxVMzdT66SLbs6Pb8m5A2RN8QACFyDgqFHNIgzs6mlVya8MnBUTIl8ujtZGpcAACAASURBVLbe82tC3hB56xYDIHANCoYe0SDO3O5Yn/zx48GjYl54t0Y6u3L7qBjkrUvkrUsMgMA1KBh6RINw76rKg3LLlKVSHInKb19YJ/sac/eoGOStS+StSwyAwDUoGHpEg8iMu+tjcv+MtamjYso2H/L8mpA3RN66xAAIXIOCoUc0iMzZ3nFMnnmrOnVUzFuf78y5o2KQty6Rty4xAAYEa+1kZq4loipmXmCMObugoOACIlpKRDVEVG2tvcIYY4ionIi2M3MZM5dZa6cnv8a5RLSQmSuJqMZae+dwnhsFQ49oEJk1Hu+XD1cOHhXzyBsVcqgtd46KQd66RN66xAAYABzHuZaIdhhjzjLGGCJabK29l4jeZuYpxhhjrf0xEU1N3l/OzN8//usQ0QRr7VxjjMnPz7+QiFodx7k43fOjYOgRDSI7Vm1vl188vkKKI1G556lVUrc7N46KQd66RN66xAAYDEbl5eWdP/ABEc0hoolEdKygoOCC4x+c/AngD050u+M4Pxry8UJr7dh0T46CoUc0iOzZeKhXHn55vRRHovLTiSXyWWmD59eEvHWJvHWJATBgOI5jiagtHA5fw8xfWGsfSv7Eb0UoFPqWMakBcAkRrSWi1eFw+LvGGMPMDQNr4uTjZlprJ6V7ThQMPaJBZNfuWJ+8vmhb6qiYWQtqpbOrD3lD5A2zknf2phEwohQVFV1FRI3hcPg6x3EuJqL/dBznvxtjjLV2LDPXJf97jOM43zbGGMdxfsjMncaYs44fAK21s4hoYrrnFRHp7U38YYLBtrc30SCQd3ZdWdkkP5ucOCrm9zPXSX3TYeQNkTfMeN5ZGkfASGKtvZqZG8Lh8DXJm0Yx8xcD9xcUFPwXZv6PE30uEcUvu+yy0cxcyszXD9zOzIuY+bZ0zy0AgIzTHv9CHnx+nRRHonL7YytkW33c60sCAAQM18MH8JbRo0efx8z11torh95ORJ9Ya/+HMcYw8/XJnwCexcwVzJyffMx3iKjDGDOKiMYR0XxjjAmFQpcQUduJ3kN4PCL4G6MW8ROCkbW986g8/eeq1FExby/ZJfHD/cgbIm+YkbwzP5GAEcVaex8z9zBzWfL9fmXMPKWoqCiPiFYzcwURbWDm/2aMMY7j3MjMdcnHVlpr/yH5pc4hovnMvCl5bMyY4Ty/CN4zosWeHrxHaKSNx/vl/RV7U0fFPPbHSmkZoaNikLcukbcue3owAAKXoGDoEQ3COyu3tsnY5FEx9z61Wrbu6ULeEHlDV3l7PT8An4OCoUc0CG9tbO6RyEuJo2Jumlgi0bLsHhWDvHWJvHWJARC4BgVDj2gQ3tvV3SevfbQ1dVTMi+9tydpRMchbl8hblxgAgWtQMPSIBpE7Lt/YJDcnj4r5x1llUt90GHlD5A1PK2+v5wfgc1Aw9IgGkVvu2Nctv56+JnlUzHLZUNuCvCHyhsPO2+v5AfgcFAw9okHknq3tR2Xa3ORRMeOTR8XE+5E3RN4wbd5ezw/A56Bg6BENIjeNx/vlvWV75CfjE0fFTJ1TKS1tR5E3RN7wlHl7PT8An4OCoUc0iNy2oq5Vxk5NHBXzq6dXyzaXR8Ugb10ib11iAASuQcHQIxpE7nuguUf+aXZ54qiYSSWypPwA8obIG54wb6/nB+BzUDD0iAbhD7u6++SVDwaPinlpYZ10dZ/+UTHIW5fIW5cYAIFrUDD0iAbhL5dtaJSbJyWOihn3Ypk0HOxB3hB5w1TeXs8PwOegYOgRDcJ/bt/bJfc9kzgq5o7HlsvGLa3IGyJviAEQuAcFQ49oEP60tf2oPPHmJimOROUn45fIu0t3D+uoGOStS+StSwyAwDUoGHpEg/Cv8Xi/LFi6R24Yn3hf4BNvbpLW9lMfFYO8dYm8dYkBELgGBUOPaBD+d+OWVrnjseWJo2KeWSPb9578qBjkrUvkrUsMgMA1KBh6RIMIhg0He+QPLyaOirl50lJZtqEReUPkrUwMgMA1KBh6RIMIjl3dffLS+3Wpo2Je/uDrR8Ugb10ib11iAASuQcHQIxpE8CwpPyA3TSqR4khU/jC7XA40Dx4Vg7x1ibx1iQEQuAYFQ49oEMF0294u+dXTqxNHxUxdIRV1rchbochblxgAgWtQMPSIBhFcW9qOytQ5lamjYhYs3SOHD/cjb0Xi9a1LDIDANSgYekSDCLbxeL+8U7I7dVTMtLmb5F//7X8hbyXi9a1LDIDANSgYekSD0OGG2ha5PXlUzAPPlcr2fSc/KgYGR7y+dYkBELgGBUOPaBB6rG86LONmlUlxJCo3TSyRT0vrh/Wvh0D/ite3LjEAAtegYOgRDUKX3bE+eWvJ7tRRMc/O3yztHcc8vy6YHfH61iUGQOAaFAw9okHociDv1ZsOyq1TlklxJCq/nn7qfz0E+le8vnWJARC4BgVDj2gQuhya974Dcfn9zORKeFKJfL6uwfPrg9nL2+trgSOTt9fzA/A5KBh6RIPQ5fF5d3X3ySsfbE2thGe8vVnaO7ESDop4fesSAyBwDQqGHtEgdHmyvFdWHJRbkivh+59dKzv2dXt+rTB7ecNgigEQuAYFQ49oELo8Vd57D8Tldy+sS62Eo2VYCftdvL51iQEQuAYFQ49oELpMl3dnV5+8/H5daiX83Ds1WAn7WLy+dYkBELgGBUOPaBC6HG7eKzY2yS1TliZWwjPWys79WAn7Uby+dYkBELgGBUOPaBC6PJ289zTE5LfPJ1bCN09aKkvKD3h+/TB7eUP/iwEQuAYFQ49oELo83bw7u/rkpYWDK+EX3q2VDqyEfSNe37rEAAhcg4KhRzQIXZ5p3ss3NsnPJidWwg/MKJVd9THPvxeYvbyhP8UAGBCstZOZuZaIqph5gTHm7IKCgguIaCkR1RBRtbX2CmOMKSwsLGLmMiLayMwV1tqrk1/jXCJayMyVRFRjrb1zOM+NgqFHNAhdusl7d0NMHnquNLESnrxUStZjJZzr4vWtSwyAAcBxnGuJaIcx5ixjjCGixdbae4nobWaeYowx1tofE9HU5P3Lieh2Y4xh5u8z887k7ROttXONMSY/P/9CImp1HOfidM+PgqFHNAhdus27s+uYvPjeltRKeNaCWunswko4V8XrW5cYAIPBqLy8vPMHPiCiOUQ0kYiOFRQUXHDcY88mor8YY84Z8vh2a20BEZU7jvOjIbcvtNaOTffkKBh6RIPQZabyXrahMbUSfvC5UtmNlXBOite3LjEABgzHcSwRtYXD4WuY+Qtr7UNEVM7MK0Kh0LcKCwsvJaLDQz+HmWuttX/PzA0Da2JjjCGimdbaSemeEwVDj2gQusxk3rvrY/LgkJXwsg2Nnn9/MHt5w9wXA2CAKCoquoqIGsPh8HWO41xMRP/pOM5/N8YYa+1YZq4rLCy8lJl7hn4eM9cx8/eOHwCttbOIaGK65xUR6e1N/GGCwba3N9EgkLcOM513Z9cxmfVebWol/OJ7W6Sr+5jn3yfMTt4wt+3txQAYCKy1VzNzQzgcviZ50yhm/mLg/oKCgv/CzP9hjPkrIvpXY8xfD9xHRLFQKHQJM5cy8/UDtzPzIma+Ld1zCwAAnAbrt7anDo7+xxfLpbPnn72+JABUksExBHjB6NGjz2PmemvtlUNvJ6JPrLX/wxhjmPl6Zq5L3v4ZM99lTOqXQzYnbx9HRPONMSYUCl1CRG0neA/h1xDB3xi1iJ8Q6DKbee+q75YHZqyV4khUfjZ5qazY2OT596tdvL51iZ8ABgBr7X3M3JM82qWcmcuYeUpRUVEeEa1m5goi2sDM/y35+AIiWktEG4moPBQKfSv5pc4hovnMvCl5bMyY4Ty/CN4zosWeHrxHSJPZzruj85i88O7gSnj2wi3S2dXn+fetVby+ddnTgwEQuAQFQ49oELocqbyXlB+QmyclVsK/fWGd7GnAbwkHOW+YG2IABK5BwdAjGoQuRzLvnfu75f7kSviWKYmVsNffvzbx+tYlBkDgGhQMPaJB6HKk827vPCbPv1OTWgm//H4dVsIBzht6n7fX8wPwOSgYekSD0KVXeX++rkFumlQixZGo/H7mOtl7IO75/wsN4vWtSwyAwDUoGHpEg9Cll3lv39ctv3l2YCW8TFZWHPT8/0fQxetblxgAgWtQMPSIBqFLr/Nu7zgmM97enFoJv/LBVunqxko4qHnDkc/b6/kB+BwUDD2iQegyV/L+rLRBbpo4sBIuk31YCQc6bzhyeXs9PwCfg4KhRzQIXeZS3tv3dsmvp6+R4khUbp2yTFZXYiUc5LzhyOTt9fwAfA4Khh7RIHSZa3m3dRyV6fMGV8KvfoiVcJDzhtnP2+v5AfgcFAw9okHoMhfzjsf75dO19fLT5Ep43Kwy2deIlXBQ84bZzdvr+QH4HBQMPaJB6DKX8962p0vueyaxEv75I8tkTVWz59fkd3M5b5idvL2eH4DPQcHQIxqELnM977aOo/LMvOrUSvj1j7ZhJRzgvGHm8/Z6fgA+BwVDj2gQuvRD3vF4vyweuhJ+sUzqmw57fl1+1A95w8zm7fX8AHwOCoYe0SB06ae8t+7pkl8NWQmXVmMlHOS8YWby9np+AD4HBUOPaBC69Fvere1H5em3BlfCb3yMlXCQ84bu8/Z6fgA+BwVDj2gQuvRj3vF4v3y8Zr/cOCGxEv7D7HKshAOcN3SXt9fzA/A5KBh6RIPQpZ/zrtvdKfc+vVqKI1G57dHlUlp9yPNrynX9nDc8s7y9nh+Az0HB0CMahC79nndr+1F56s9VqZXwnE+2S3cMK+Gg5g1PP2+v5wfgc1Aw9IgGocsg5B2P98uiVfvkxglLpDgSlcjscmk42OP5deWiQcgbnl7eXs8PwOegYOgRDUKXQcp7y65OufepwZVw2WashIOcNxxe3l7PD8DnoGDoEQ1Cl0HLu6XtqEybm1gJ3xCJypuLsRIOct4wfd5ezw/A56Bg6BENQpdBzDse75cPVw6uhMe/tB4r4QDnDU+dt9fzA/A5KBh6RIPQZZDzrt3VIXc/uUqKI1G5/bHlUl7T4vk1eW2Q84Ynztvr+QH4HBQMPaJB6DLoebe0HZEn/rQptRL+0+IdEov1e35dyBuOVN5ezw/A56Bg6BENQpca8o7H++WDFXtlzMOJlfDDL6+XA806V8Ia8oZfzdvr+QH4HBQMPaJB6FJT3jU7B1fCd0xdIetr9a2ENeUNMQCCDICCoUc0CF1qy7ul7Yg8/ubgSvjPn+1UtRLWlrd2MQAC16Bg6BENQpca847H++X95YMr4QmvbJDGQ72eXxfyhtnI2+v5AfgcFAw9okHoUnPem3d0yF3TEivhsVNXyIYtrZ5fE/KGmc7b6/kB+BwUDD2iQehSe97NrUfksT9WJlbC46Py1ufBXglrz1ubGACBa1Aw9IgGoUvknVgJv7dsT2olPOnVDdIU0JUw8tYlBkDgGhQMPaJB6BJ5D1q9vV3umrYysRJ+fIVU1AVvJYy8dYkBELgGBUOPaBC6RN5ftbn1iDz6RkVqJTwvukti8eCshJG3LjEAAtegYOgRDUKXyPvrxuP9smDpHvnJ+ORK+LWN0tQSjJUw8tYlBkDgGhQMPaJB6BJ5n9yqbe1y5xOJlfAvHl8hlVvbPL8m5A1PN2+v5weQAay1k5m5loiqmHmBMeYcIvqSmcuIqJyZy6y1txhjDBGtJ6LtzFyWvH168mucS0QLmbmSiGqstXcO57lRMPSIBqFL5H1qm1uPyCPJlfBPxi+Rt5f4eyWMvHWJATAAOI5zLRHtMMacZYwxRPSJtfZeIvrPEz0+ORB+/wS3T7DWzjXGmPz8/AuJqNVxnIvTPT8Khh7RIHSJvNMbi/fLuyW7Uyvhya9vlIM+XQkjb11iAAwGo/Ly8s4f+ICI5hDRxFMNgET0gxPd7jjOj4Z8vNBaOzbdk6Ng6BENQpfIe/hu2taWWgnf+cRK2bTNfyth5K1LDIABw3EcS0Rt4XD4ciL6kojmEdFGIlo48NO85AC4hIjWEtHqcDj8XWOMYeYGa+0VA1+LiGZaayele04UDD2iQegSeZ+eB1t6ZfLrG1Mr4XdLdvtqJYy8dYkBMEAUFRVdRUSN4XD4OmOMIaIH8vPzLzTGGGaexswfG2OMtXaM4zjfNsYYx3F+yMydxpizjh8ArbWziGhiuucVEentTfxhgsG2tzfRIJC3DpH36Rs/3C/vlOxKrYQfeaNCDrUe8fy6kDc8Ud5ZGUbAyGKtvZqZG8Lh8DUnuf8KIjpwovuIKH7ZZZeNZuZSZr5+4HZmXsTMt6V7bgEAAPAV9hzslV8m/y3hu59cJfuaj3h9SQB8jTOfOkBOMHr06POYud5ae+XAbUT0d0S0zBhzdvLjccy8yCR+0lfBzPnJ279DRB3GmFFENI6I5htjTCgUuoSI2goKCi5I9/wi+BujFvETAl0ib3cebO2Vya8lVsJjHl4iC5btkfjhfs+vC3nDgbyzMpSAkcNaex8z9ww98oWZp1hrnyCircn3/C0tLCy81BhjHMe5kZnrko+ttNb+Q/JLnUNE85l5ExFVW2vHDOf5RfCeES329OA9QppE3u6NxftlXnSX3DA+KsWRqDz6RoU0tx7x/LqQN+zpwQAIXIKCoUc0CF0i78xZUdcqYx9fIcWRqNw1baVUb2/3/JqQt24xAALXoGDoEQ1Cl8g7szYd6pVJr25IrYTfW7ZH4jn0W8LIW5cYAIFrUDD0iAahS+SdeWOxfnnr852plfBjf6yUQ225sRJG3rrEAAhcg4KhRzQIXSLv7LlhS6uMnTqwEl4lm3d0eH5NyFuXGACBa1Aw9IgGoUvknV0bD/XKhFcGV8ILl+/1dCWMvHWJARC4BgVDj2gQukTe2TcW65c/f7ZTbogkVsJT53i3EkbeusQACFyDgqFHNAhdIu+Rc31ti9yRXAnf/eQqqdk58ith5K1LDIDANSgYekSD0CXyHlkbm3vk4ZfXp1bCH6wY2ZUw8tYlBkDgGhQMPaJB6BJ5j7yxWL/M/XRHaiX8xJubpGWEVsLIW5cYAIFrUDD0iAahS+TtneU1LXL7Y8ulOBKVe55aJbW7sr8SRt66xAAIXIOCoUc0CF0ib29tONgj419KrIRvnLBEPly5L6srYeStSwyAwDUoGHpEg9Al8vbe7lifvLl4e2olPG1ulbS0HUXeMCN5ez0/AJ+DgqFHNAhdIu/csWzzIbnt0cRK+N6nVsuWXZ3IG7rO2+v5AfgcFAw9okHoEnnnlg0HeyQyuzy1El60KrMrYeStSwyAwDUoGHpEg9Al8s49u2N9MueT7VKcXAk/9ecqaW3PzEoYeesSAyBwDQqGHtEgdIm8c9d1Q1bCv3p6tdTtdr8SRt66xAAIXIOCoUc0CF0i79y2vumw/CG1Ei6Rj1fvd7USRt66xAAIXIOCoUc0CF0i79y3q7tP3vh4W2ol/PRb1We8EkbeusQACFyDgqFHNAhdIm//WFrdLD9/ZFliJfzMGtm6pwt5w7R5ez0/AJ+DgqFHNAhdIm9/Wd90WMa9WCbFkaj8dGKJLF5bf1orYeStSwyAwDUoGHpEg9Al8vafXd198vqiwZXwM/Oqpa1jeCth5K1LDIDANSgYekSD0CXy9q9rNg2uhO97Zo1sG8ZKGHnrEgMgcA0Khh7RIHSJvP3t/sbDMm5WYiV808QS+TTNShh56xIDIHANCoYe0SB0ibz9b1d3n7z64dbUSnj6vM3S3nEMeUMMgMA9KBh6RIPQJfIOjqsrD8qtUxIr4V9PXyPb9n59JYy8dYkBELgGBUOPaBC6RN7Bct+BuPx+5uBK+LPSBuStWAyAwDUoGHpEg9Al8g6eXd198soHgyvhGfMHV8LIW5cYAIFrUDD0iAahS+QdXFdVHpRbkivh3zy7Vrbv60beysQACFyDgqFHNAhdIu9gu/dAXH73wrrESnhSiUTLGpC3IjEAAtegYOgRA4EukXfw7ezqk5ffr0uthF/5aJt0dJ34t4RhsMQACFyDBqFHDAS6RN56XFnRlFoJ3z9jrezc3+35NcHsigEQuAYNQo8YCHSJvHW590Bc/jC7XIojUbl50lJZUnbA82uC2RMDIHANGoQeMRDoEnnrsqfnC/lf//v/k5eGrIRfeLdGOjqxEg6iGACBa9Ag9IiBQJfIW5dD816+sUl+NnmpFEei8gBWwoEUAyBwDRqEHjEQ6BJ56/L4vHc3xOSh50tTK+GScqyEgyQGwIBgrZ3MzLVEVMXMC4wx5xDRl8xcRkTlzFxmrb3FGGMKCwuLkrdvZOYKa+3Vya9xLhEtZOZKIqqx1t45nOdGg9AjBgJdIm9dnijvzq5j8uJ7W1Ir4ZkLarESDogYAAOA4zjXEtEOY8xZxhhDRJ9Ya+8lov880eOJaDkR3W6MMcz8fWbembx9orV2rjHG5OfnX0hErY7jXJzu+dEg9IiBQJfIW5enynvZhsbUSvjB50plV33M8+uF7vPO2CACPGNUXl7e+QMfENEcIpp4kgHwbCL6izHmnCGPb7fWFhBRueM4Pxpy+0Jr7dh0T44GoUcMBLpE3rpMl/fu+pg89FxyJTx5qSzb0Oj5NUN3ebsdPkAO4TiOJaK2cDh8ORF9SUTziGgjES10HOfiwsLCS4no8NDPYeZaa+3fM3ODtfaKgduJaKa1dlK650SD0CMGAl0ib10OJ+/OrmMya0FtaiU8670t0omDo30pBsAAUVRUdBURNYbD4euMMYaIHsjPz7/QGGOYeRozLyosLLyUmXuGfh4z1zHz944fAK21s4hoYrrnFRHp7U38YYLBtrc30SCQtw6Rty5PJ++l6w+kVsIPPV8qexpinl8/PP28Mz2HAA+w1l7NzA3hcPiak9x/BREdMMb8FRH9qzHmrwfuI6JYKBS6hJlLmfn6gduZeREz35buuQUAAIA62uNfyO9mlklxJCq3TFkmFTs6vL4kcJpkYv4AHjJ69OjzmLneWnvlwG1E9HdEtMwYc3by43HMvCj5358x813GGGOt/TERbR54DBHNN8aYUCh0CRG1FRQUXJDu+UXwEwIt4idCukTeujyTvDu7jsnMISvh2Qu3SFd3n+ffCxxe3lkYScBIYq29j5l7hh75wsxTrLVPENFWIionoqWFhYWXJh9fQERrk+8NLA+FQt9KfqlziGg+M28iompr7ZjhPL8I3iOkxZ4evCdMk8hbl27yXlJ+QG6elFgJ//b5dbKnAb8lnOv29GAABC5Bg9AjBgJdIm9dus175/5uuX/G2uRKeKms2Njk+fcET5231/MD8DloEHrEQKBL5K3LTOTd3nlMnn+nJrUSfun9Ouns6vP8e4Mnztvr+QH4HDQIPWIg0CXy1mUm8/58XYPcNKlEiiNR+d0L62TPgbjn3x/8et5ezw/A56BB6BEDgS6Rty4znfeOfd1y/7NrU78lvLLioOffI/xq3l7PD8DnoEHoEQOBLpG3LrORd3vHMZnx9ubUSvjlD+qkqxsr4VwQAyBwDRqEHjEQ6BJ56zKbeX9W2iA3TUyshH8/c53sw0rYczEAAtegQegRA4Eukbcus5339r1d8uvpa6Q4EpVbpyyT1ZVYCXudt9fzA/A5aBB6xECgS+Sty5HIu73jmDw7f3Al/OqHW7ES9jBvr+cH4HPQIPSIgUCXyFuXI5V3PN4vn5bWp1bC/zirTPY1YiXsRd5ezw/A56BB6BEDgS6Rty5HOu9tQ1bCP39kmazZ1Oz5/wNNYgAErkGD0CMGAl0ib116kXdbx1F5Zl51aiX8+kfbsBIewby9nh+Az0GD0CMGAl0ib116lXc83i+L19bLT5Mr4XEvlsn+xsOe//8IuhgAgWvQIPSIgUCXyFuXXue9dU+X/OqZwZXw2iqshLOdt9fzA/A5aBB69LpBQOQNg513a/tRefqtwZXwG4uwEs5m3l7PD8DneF0w4MgWDOStR+Sty1zJOx7vl4/X7JcbJyRWwn+YXS71TVgJZyNvr+cH4HNyoWDAkSsYyFuPyFuXuZZ33e5O+dXTq6U4EpXbHl0updWHPL+mIIkBELgmlwoGzH7BQN56RN66zMW8W9uPylN/rkqthP/48XbpjmElnKm8vZ4fgM/JtYIBs1swkLcekbcuczXveLxfFq0eXAn/0+xyaTjY4/l1+V0MgMA1uVgwYPYKBvLWI/LWZa7nvWVXp9z71OBKuGwzVsJu8/Z6fgA+J5cLBsx8wUDeekTeuvRD3q3tR2Xa3MGV8JxPsBJ2k7fX8wPwObleMGBmCwby1iPy1qVf8o7H++WjVfvkxglLpDgSlchL67ESPsO8vZ4fgM/xQ8GAmSsYyFuPyFuXfsu7dleH3PPUKimOROX2x5ZLeU2L59fkJzEAAtf4qWBA9wUDeesReevSj3m3tB2VJ/60SYojUbkhEpU/Ld4hsVi/59flBzEAAtf4rWBAdwUDeesReevSr3nH4/3y4cp9MubhxEr44ZfXy4FmrISHk7fX8wPwOX4sGPDMCwby1iPy1qXf867Z2SF3P5lYCd/x2HJZX4uVcLq8vZ4fgM/xc8GAp18wkLcekbcug5B3S9sRefzNwZXw3E+xEj5V3l7PD8Dn+L1gwNMrGMhbj8hbl0HJOx7vl/eX702thCe8skEasRI+Yd5ezw/A5wShYMDhFwzkrUfkrcug5b15R4fcNS25Ep66QjZgJfy1vL2eH4DPCVLBgOkLBvLWI/LWZRDzPtR2RKbOqUyshMdH5a3PdmIlPCRvr+cH4HOCVjDgqQsG8tYj8tZlUPOOx/tl4bI9qZXwxFc3SOOhXs+vy2sxAALXBLFgwJMXDOStR+Sty6DnXb29Xe6atlKKI1EZO3WFbNzS6vk1eZ231/MD8DlBLhjw6wUDeesReetSQ97NrUfksT8OroTnfb5TYnGdK2EMgMA1QS8Y8KsFA3nrEXnrUkve8Xi/LFg6uBKe9OoGaVK4EsYACFyjoWDAwYKBdPsnQwAADa1JREFUvPWIvHWpLe+q7e3yy+RK+BePr5CKOl0rYQyAAcFaO5mZa4moipkXGGPOHriPiG4moi+HfFxORNuZuYyZy6y105Nf41wiWsjMlURUY629czjPralgaFdbg9Au8talxrybW4/II29UpFbCby/ZpWYljAEwADiOcy0R7TDGnGWMMUS0mIjuSd53MRFtJqKOgccTUTkzf//4r0NEE6y1c40xJj8//0IianUc5+J0z6+tYGhWY4PQLPLWpda8Y/F+eXfpbvnJ+MRKePJrG6WpJfgrYQyAwWBUXl7e+QMfENEcIppgjDHMHCWiHzBz85D7y4noB8d/ESIqdxznR0M+XmitHZvuyTUWDK1qbRBaRd661J73pm1tcucTiZXwnU+slMqtbZ5fU7bzdj9+gJzBcRxLRG3hcPhya+29zDzbGGOY+dDAY5ID4BIiWktEq8Ph8HeTj2mw1l4x5HEzrbWT0j2n5oKhTe0NQpvIW5fI+ws52NIrU15PrIR/Mn6JvFOyO7ArYQyAAaKoqOgqImoMh8PXFRUVFTLzllAo9DfGfHUAtNaOcRzn28YY4zjOD5m50xhz1vEDoLV2FhFNTPe8IiK9vYk/TDDY9vYmGgTy1iHy1iXyThg/3C/vlOxKrYQfeaNCmlt7Pb+ubOSdhVEEjDTW2quZuSEcDl9jjDFENI6IdhNRFRFVE9G/E1HVRRdd9I3jP5eI4pdddtloZi5l5usHbmfmRcx8W7rnFgAAACBg7DnYK79M/lvCdz+5SvY293p9SRkns5MIGHFGjx59HjPXW2uvPNljhvwE8CxmrmDmfGOMIaLvJH9BZFRyaJxvjDGhUOgSImorKCi4IN3zi+BvjFrETwh0ibx1iby/7sHWXpn82sbUSnjB0t0SP9zv+XVlKu9MzCDAQ6y19zFzDzOXJX/Dt4yZpwx9zNBfAnEc50Zmrks+ttJa+w/Ju84hovnMvImIqq21Y4bz/CK63zOiyZ4evEdIk8hbl8j7xMbi/TI/uktuGB8dshI+4vl1ZSLvDI4iQCMoGHpEg9Al8tYl8j61FXWtMvbxFVIcicovp62U6u3tnl+T27y9nh+Az0HB0CMahC6Rty6Rd3qbDvXKpFc3SHEkKmMeXiILlu6RuE9/SxgDIHANCoYe0SB0ibx1ibyHZyzeL/M+35laCT/2x0pfroQxAALXoGDoEQ1Cl8hbl8j79NywpVXGTk2shO+atkqqd/hrJYwBELgGBUOPaBC6RN66RN6nb+OhXpnwyuBKeOHyvb5ZCWMABK5BwdAjGoQukbcukfeZGYv1y1uf7ZQbIomV8NQ5lXKoLfdXwhgAgWtQMPSIBqFL5K1L5O3ODbUtckdyJXz3k6ukZmeH59eULm+v5wfgc1Aw9IgGoUvkrUvk7d7G5h55+OX1qZXw+ytydyWMARC4BgVDj2gQukTeukTemTEW65e5n+5IrYSfeHOTtOTgShgDIHANCoYe0SB0ibx1ibwza3lti9zx2PKcXQljAASuQcHQIxqELpG3LpF35m042CPjXxpcCX+QQythDIDANSgYekSD0CXy1iXyzo7dsT55c/HgSnjanzZJS9tRz68LAyBwDQqGHtEgdIm8dYm8s2tZzSG57dHESviep1bJll2dnuft9fwAfA4Khh7RIHSJvHWJvLNvw8EeiSRXwjdOWCIfrdrn2UoYAyBwDQqGHtEgdIm8dYm8R8buWJ/M+WS7FCdXwk/OrZLW9pFfCWMABK5BwdAjGoQukbcukffIum7z4Er43qdXj/hKGAMgcA0Khh7RIHSJvHWJvEfe+qbD8k+zy5Mr4RJZtHr/iK2EMQAC16Bg6BENQpfIW5fI2xu7Y33yx48HV8JP/XlkVsIYAIFrUDD0iAahS+StS+TtraXVzfLzR5ZJcSQqv3p6tdTtzu5KGAMgcA0Khh7RIHSJvHWJvL23vumw/OHFwZXwx2uytxLGAAhcg4KhRzQIXSJvXSLv3LCru09eX7QttRJ+5q1qaevI/EoYAyBwDQqGHtEgdIm8dYm8c8s1VYMr4fueWSNb93RlPG+v5wfgc1Aw9IgGoUvkrUvknXvubzws42aVSXEkKj+dWCKL19ZnbCWMARC4BgVDj2gQukTeukTeuWlXd5+89tHW1Ep4+rzMrIQxAALXoGDoEQ1Cl8hbl8g7t1296aDcOiWxEv719DWyba+7lTAGQOAaFAw9okHoEnnrEnnnvvsOxOX3MxMr4ZsmlsinpfWu8vZ6fgA+BwVDj2gQukTeukTe/rCru09e/XBwJTxj/mZp7zh2Rnl7PT8An4OCoUc0CF0ib10ib3+5qvKg3JJcCf/m2bWyfV/3aeft9fwAfA4Khh7RIHSJvHWJvP3n3gNx+f3MdYmV8KQS+Xxdw2nl7fX8AHwOCoYe0SB0ibx1ibz9aVd3n7z8QV1qJfzc2zXS3pl+JYwBELgGBUOPaBC6RN66RN7+dmVFU2olfP+za2VHmpUwBkDgGhQMPaJB6BJ56xJ5+989B+Ly2xcGV8LRspOvhDEAAtegYOgRDUKXyFuXyDsYdnb1yUsLB1fCz79TIx0nWAljAASuQcHQIxqELpG3LpF3sFyxsUlumbJUiiNReWDGWtm5/6srYQyAwDUoGHpEg9Al8tYl8g6euxti8tDzpVIcicrNk5ZKSfmBr+Tt9fwAMoC1djIz1xJRFTMvMMacPXAfEd1MRF8OfFxYWFjEzGVEtJGZK6y1Vye/xrlEtJCZK4moxlp753CeGwVDj2gQukTeukTewbSz65i8uHBLaiX8wru10tF5DANgEHAc51oi2mGMOcsYY4hoMRHdk7zvYiLaTEQdA48nouVEdLsxxjDz95l5Z/L2idbaucYYk5+ffyERtTqOc3G650fB0CMahC6Rty6Rd7BdtqFRfjZ5YCVcKrvrYxgAA8CovLy88wc+IKI5RDTBGGOYOUpEP2Dm5uTdZxPRX4wx5wx5fLu1toCIyh3H+dGQ2xdaa8eme3IUDD2iQegSeesSeQff3fUxeei5xEr4Z5OXYgAMEo7jWCJqC4fDl1tr72Xm2cYYw8yHjDGmsLDwUiI6PPRzmLnWWvv3zNxgrb1i4HYimmmtnZTuOVEw9IgGoUvkrUvkrcPOrmMy673ESjjTMwjwiKKioquIqDEcDl9XVFRUyMxbQqHQ3xiTdgCsY+bvHT8AWmtnEdHEdM8rItLbmygeMNj29iYaBPLWIfLWJfLW5ca6FgyAQcBaezUzN4TD4WuMMYaIxhHRbiKqIqJqIvp3Iqqy1v4tM/+LMeavBz6XiGKhUOgSZi5l5usHbmfmRcx8mwffDgAAAAAAOBWjR48+j5nrrbVXnuwxAz8BNMYYIvqMme8yxhhr7Y+JaHPy9nFENN8YY0Kh0CVE1FZQUHBBtq8fAAAAAACcJtba+5i5J3m0SzkzlzHzlKGPGfJLICb5Cx9riWgjEZWHQqFvJe86h4jmM/MmIqq21o4Z0W8EAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAcJow8xQiqkn+lvBLx9/vOM7/k/x3hzcQUUk4HP6vXlwnyAzp8iaiL4f+5rm19hYvrhNkhnA4/F+Z+SNm7j7R/Xh9B4t0eeP1HSystZOZuZaIqph5gTHm7KH34/UNTko4HL6GiHaYxAHSo4hojeM4Nw7cb609l5k7Q6FQyBhjiGgqEb3i0eUCl6TL2xhjiOg/vbk6kA2YeQUz/4aZu46/D6/v4HGqvI3B6ztIOI5zbbKen2WMMUS0mIjuGbgfr29wSqy1TxDRU0M+vo+Z/zzk4/+biDYOfJz8N4mbRvo6QWZIl7cxiZ8QjPyVgWxx0UUXfSP5T0ieaADE6ztgnCpvY/D6Dhij8vLyzh/4gIjmENGEgY/x+ganhIj+ZK19aMjH/y8zrxzy8e3M/PHAx3l5eecT0b+N9HWCzJAu7+RtXxLRvOSB4guLioryRv5KQSY52UCA13cwSTcA4vUdPJLDXVs4HL584Da8vsEpOcFAUMzMK4Z8/JU/QN/85jf/DyL6y0hfJ8gM6fJO3vZAfn7+hcYYw8zThuYP/MlwB0C8voNBmgEQr++AUVRUdBURNYbD4euG3o7XNzglRPQYET0z8LG19kEienPg43A4/H8x86Yhj/87Zq4f6esEmSFd3sdjrb2CiA6MzNWBbHGygQCv72ByqgFwKHh9+x9r7dXM3BAOh685/j68vsEpIaLvENHuUCj0N8aYs4log+M4/3PIQ84holbHcawxxlhrn2Xm57y5WuCWdHkT0d8R0TKT/E0yIhrHzIu8ul6QGUKhUOgkvxWK13cAOVneeH0Hi9GjR5/HzPXW2itP8hC8vsGpIaJxRFST/DXyJ5O3vRQOh79rjDGO4/yQiKqJaCMzfzz0TafAf6TLO/mLIluJqJyIlhYWFl7q7RWDM+Xyyy//P5M5VhHRvyeP/3iFmWfj9R08hpM3Xt/BIflLfD1Dj/Vh5il4fQMAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAADg/28PDkgAAAAABP1/3Y5ABQAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAACAiwA+9L+IbYlIKgAAAABJRU5ErkJggg==\">"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/plain": [
       "[<matplotlib.lines.Line2D at 0x7fc136ffab38>]"
      ]
     },
     "execution_count": 8,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "%matplotlib notebook\n",
    "fig = plt.figure()\n",
    "plt.plot(range(len(riqueza_milpa)), biomasa)\n",
    "#plt.xlim(0, 4)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {
    "collapsed": false,
    "scrolled": false
   },
   "outputs": [
    {
     "data": {
      "application/javascript": [
       "/* Put everything inside the global mpl namespace */\n",
       "window.mpl = {};\n",
       "\n",
       "mpl.get_websocket_type = function() {\n",
       "    if (typeof(WebSocket) !== 'undefined') {\n",
       "        return WebSocket;\n",
       "    } else if (typeof(MozWebSocket) !== 'undefined') {\n",
       "        return MozWebSocket;\n",
       "    } else {\n",
       "        alert('Your browser does not have WebSocket support.' +\n",
       "              'Please try Chrome, Safari or Firefox ≥ 6. ' +\n",
       "              'Firefox 4 and 5 are also supported but you ' +\n",
       "              'have to enable WebSockets in about:config.');\n",
       "    };\n",
       "}\n",
       "\n",
       "mpl.figure = function(figure_id, websocket, ondownload, parent_element) {\n",
       "    this.id = figure_id;\n",
       "\n",
       "    this.ws = websocket;\n",
       "\n",
       "    this.supports_binary = (this.ws.binaryType != undefined);\n",
       "\n",
       "    if (!this.supports_binary) {\n",
       "        var warnings = document.getElementById(\"mpl-warnings\");\n",
       "        if (warnings) {\n",
       "            warnings.style.display = 'block';\n",
       "            warnings.textContent = (\n",
       "                \"This browser does not support binary websocket messages. \" +\n",
       "                    \"Performance may be slow.\");\n",
       "        }\n",
       "    }\n",
       "\n",
       "    this.imageObj = new Image();\n",
       "\n",
       "    this.context = undefined;\n",
       "    this.message = undefined;\n",
       "    this.canvas = undefined;\n",
       "    this.rubberband_canvas = undefined;\n",
       "    this.rubberband_context = undefined;\n",
       "    this.format_dropdown = undefined;\n",
       "\n",
       "    this.image_mode = 'full';\n",
       "\n",
       "    this.root = $('<div/>');\n",
       "    this._root_extra_style(this.root)\n",
       "    this.root.attr('style', 'display: inline-block');\n",
       "\n",
       "    $(parent_element).append(this.root);\n",
       "\n",
       "    this._init_header(this);\n",
       "    this._init_canvas(this);\n",
       "    this._init_toolbar(this);\n",
       "\n",
       "    var fig = this;\n",
       "\n",
       "    this.waiting = false;\n",
       "\n",
       "    this.ws.onopen =  function () {\n",
       "            fig.send_message(\"supports_binary\", {value: fig.supports_binary});\n",
       "            fig.send_message(\"send_image_mode\", {});\n",
       "            fig.send_message(\"refresh\", {});\n",
       "        }\n",
       "\n",
       "    this.imageObj.onload = function() {\n",
       "            if (fig.image_mode == 'full') {\n",
       "                // Full images could contain transparency (where diff images\n",
       "                // almost always do), so we need to clear the canvas so that\n",
       "                // there is no ghosting.\n",
       "                fig.context.clearRect(0, 0, fig.canvas.width, fig.canvas.height);\n",
       "            }\n",
       "            fig.context.drawImage(fig.imageObj, 0, 0);\n",
       "        };\n",
       "\n",
       "    this.imageObj.onunload = function() {\n",
       "        this.ws.close();\n",
       "    }\n",
       "\n",
       "    this.ws.onmessage = this._make_on_message_function(this);\n",
       "\n",
       "    this.ondownload = ondownload;\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._init_header = function() {\n",
       "    var titlebar = $(\n",
       "        '<div class=\"ui-dialog-titlebar ui-widget-header ui-corner-all ' +\n",
       "        'ui-helper-clearfix\"/>');\n",
       "    var titletext = $(\n",
       "        '<div class=\"ui-dialog-title\" style=\"width: 100%; ' +\n",
       "        'text-align: center; padding: 3px;\"/>');\n",
       "    titlebar.append(titletext)\n",
       "    this.root.append(titlebar);\n",
       "    this.header = titletext[0];\n",
       "}\n",
       "\n",
       "\n",
       "\n",
       "mpl.figure.prototype._canvas_extra_style = function(canvas_div) {\n",
       "\n",
       "}\n",
       "\n",
       "\n",
       "mpl.figure.prototype._root_extra_style = function(canvas_div) {\n",
       "\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._init_canvas = function() {\n",
       "    var fig = this;\n",
       "\n",
       "    var canvas_div = $('<div/>');\n",
       "\n",
       "    canvas_div.attr('style', 'position: relative; clear: both; outline: 0');\n",
       "\n",
       "    function canvas_keyboard_event(event) {\n",
       "        return fig.key_event(event, event['data']);\n",
       "    }\n",
       "\n",
       "    canvas_div.keydown('key_press', canvas_keyboard_event);\n",
       "    canvas_div.keyup('key_release', canvas_keyboard_event);\n",
       "    this.canvas_div = canvas_div\n",
       "    this._canvas_extra_style(canvas_div)\n",
       "    this.root.append(canvas_div);\n",
       "\n",
       "    var canvas = $('<canvas/>');\n",
       "    canvas.addClass('mpl-canvas');\n",
       "    canvas.attr('style', \"left: 0; top: 0; z-index: 0; outline: 0\")\n",
       "\n",
       "    this.canvas = canvas[0];\n",
       "    this.context = canvas[0].getContext(\"2d\");\n",
       "\n",
       "    var rubberband = $('<canvas/>');\n",
       "    rubberband.attr('style', \"position: absolute; left: 0; top: 0; z-index: 1;\")\n",
       "\n",
       "    var pass_mouse_events = true;\n",
       "\n",
       "    canvas_div.resizable({\n",
       "        start: function(event, ui) {\n",
       "            pass_mouse_events = false;\n",
       "        },\n",
       "        resize: function(event, ui) {\n",
       "            fig.request_resize(ui.size.width, ui.size.height);\n",
       "        },\n",
       "        stop: function(event, ui) {\n",
       "            pass_mouse_events = true;\n",
       "            fig.request_resize(ui.size.width, ui.size.height);\n",
       "        },\n",
       "    });\n",
       "\n",
       "    function mouse_event_fn(event) {\n",
       "        if (pass_mouse_events)\n",
       "            return fig.mouse_event(event, event['data']);\n",
       "    }\n",
       "\n",
       "    rubberband.mousedown('button_press', mouse_event_fn);\n",
       "    rubberband.mouseup('button_release', mouse_event_fn);\n",
       "    // Throttle sequential mouse events to 1 every 20ms.\n",
       "    rubberband.mousemove('motion_notify', mouse_event_fn);\n",
       "\n",
       "    rubberband.mouseenter('figure_enter', mouse_event_fn);\n",
       "    rubberband.mouseleave('figure_leave', mouse_event_fn);\n",
       "\n",
       "    canvas_div.on(\"wheel\", function (event) {\n",
       "        event = event.originalEvent;\n",
       "        event['data'] = 'scroll'\n",
       "        if (event.deltaY < 0) {\n",
       "            event.step = 1;\n",
       "        } else {\n",
       "            event.step = -1;\n",
       "        }\n",
       "        mouse_event_fn(event);\n",
       "    });\n",
       "\n",
       "    canvas_div.append(canvas);\n",
       "    canvas_div.append(rubberband);\n",
       "\n",
       "    this.rubberband = rubberband;\n",
       "    this.rubberband_canvas = rubberband[0];\n",
       "    this.rubberband_context = rubberband[0].getContext(\"2d\");\n",
       "    this.rubberband_context.strokeStyle = \"#000000\";\n",
       "\n",
       "    this._resize_canvas = function(width, height) {\n",
       "        // Keep the size of the canvas, canvas container, and rubber band\n",
       "        // canvas in synch.\n",
       "        canvas_div.css('width', width)\n",
       "        canvas_div.css('height', height)\n",
       "\n",
       "        canvas.attr('width', width);\n",
       "        canvas.attr('height', height);\n",
       "\n",
       "        rubberband.attr('width', width);\n",
       "        rubberband.attr('height', height);\n",
       "    }\n",
       "\n",
       "    // Set the figure to an initial 600x600px, this will subsequently be updated\n",
       "    // upon first draw.\n",
       "    this._resize_canvas(600, 600);\n",
       "\n",
       "    // Disable right mouse context menu.\n",
       "    $(this.rubberband_canvas).bind(\"contextmenu\",function(e){\n",
       "        return false;\n",
       "    });\n",
       "\n",
       "    function set_focus () {\n",
       "        canvas.focus();\n",
       "        canvas_div.focus();\n",
       "    }\n",
       "\n",
       "    window.setTimeout(set_focus, 100);\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._init_toolbar = function() {\n",
       "    var fig = this;\n",
       "\n",
       "    var nav_element = $('<div/>')\n",
       "    nav_element.attr('style', 'width: 100%');\n",
       "    this.root.append(nav_element);\n",
       "\n",
       "    // Define a callback function for later on.\n",
       "    function toolbar_event(event) {\n",
       "        return fig.toolbar_button_onclick(event['data']);\n",
       "    }\n",
       "    function toolbar_mouse_event(event) {\n",
       "        return fig.toolbar_button_onmouseover(event['data']);\n",
       "    }\n",
       "\n",
       "    for(var toolbar_ind in mpl.toolbar_items) {\n",
       "        var name = mpl.toolbar_items[toolbar_ind][0];\n",
       "        var tooltip = mpl.toolbar_items[toolbar_ind][1];\n",
       "        var image = mpl.toolbar_items[toolbar_ind][2];\n",
       "        var method_name = mpl.toolbar_items[toolbar_ind][3];\n",
       "\n",
       "        if (!name) {\n",
       "            // put a spacer in here.\n",
       "            continue;\n",
       "        }\n",
       "        var button = $('<button/>');\n",
       "        button.addClass('ui-button ui-widget ui-state-default ui-corner-all ' +\n",
       "                        'ui-button-icon-only');\n",
       "        button.attr('role', 'button');\n",
       "        button.attr('aria-disabled', 'false');\n",
       "        button.click(method_name, toolbar_event);\n",
       "        button.mouseover(tooltip, toolbar_mouse_event);\n",
       "\n",
       "        var icon_img = $('<span/>');\n",
       "        icon_img.addClass('ui-button-icon-primary ui-icon');\n",
       "        icon_img.addClass(image);\n",
       "        icon_img.addClass('ui-corner-all');\n",
       "\n",
       "        var tooltip_span = $('<span/>');\n",
       "        tooltip_span.addClass('ui-button-text');\n",
       "        tooltip_span.html(tooltip);\n",
       "\n",
       "        button.append(icon_img);\n",
       "        button.append(tooltip_span);\n",
       "\n",
       "        nav_element.append(button);\n",
       "    }\n",
       "\n",
       "    var fmt_picker_span = $('<span/>');\n",
       "\n",
       "    var fmt_picker = $('<select/>');\n",
       "    fmt_picker.addClass('mpl-toolbar-option ui-widget ui-widget-content');\n",
       "    fmt_picker_span.append(fmt_picker);\n",
       "    nav_element.append(fmt_picker_span);\n",
       "    this.format_dropdown = fmt_picker[0];\n",
       "\n",
       "    for (var ind in mpl.extensions) {\n",
       "        var fmt = mpl.extensions[ind];\n",
       "        var option = $(\n",
       "            '<option/>', {selected: fmt === mpl.default_extension}).html(fmt);\n",
       "        fmt_picker.append(option)\n",
       "    }\n",
       "\n",
       "    // Add hover states to the ui-buttons\n",
       "    $( \".ui-button\" ).hover(\n",
       "        function() { $(this).addClass(\"ui-state-hover\");},\n",
       "        function() { $(this).removeClass(\"ui-state-hover\");}\n",
       "    );\n",
       "\n",
       "    var status_bar = $('<span class=\"mpl-message\"/>');\n",
       "    nav_element.append(status_bar);\n",
       "    this.message = status_bar[0];\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.request_resize = function(x_pixels, y_pixels) {\n",
       "    // Request matplotlib to resize the figure. Matplotlib will then trigger a resize in the client,\n",
       "    // which will in turn request a refresh of the image.\n",
       "    this.send_message('resize', {'width': x_pixels, 'height': y_pixels});\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.send_message = function(type, properties) {\n",
       "    properties['type'] = type;\n",
       "    properties['figure_id'] = this.id;\n",
       "    this.ws.send(JSON.stringify(properties));\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.send_draw_message = function() {\n",
       "    if (!this.waiting) {\n",
       "        this.waiting = true;\n",
       "        this.ws.send(JSON.stringify({type: \"draw\", figure_id: this.id}));\n",
       "    }\n",
       "}\n",
       "\n",
       "\n",
       "mpl.figure.prototype.handle_save = function(fig, msg) {\n",
       "    var format_dropdown = fig.format_dropdown;\n",
       "    var format = format_dropdown.options[format_dropdown.selectedIndex].value;\n",
       "    fig.ondownload(fig, format);\n",
       "}\n",
       "\n",
       "\n",
       "mpl.figure.prototype.handle_resize = function(fig, msg) {\n",
       "    var size = msg['size'];\n",
       "    if (size[0] != fig.canvas.width || size[1] != fig.canvas.height) {\n",
       "        fig._resize_canvas(size[0], size[1]);\n",
       "        fig.send_message(\"refresh\", {});\n",
       "    };\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_rubberband = function(fig, msg) {\n",
       "    var x0 = msg['x0'];\n",
       "    var y0 = fig.canvas.height - msg['y0'];\n",
       "    var x1 = msg['x1'];\n",
       "    var y1 = fig.canvas.height - msg['y1'];\n",
       "    x0 = Math.floor(x0) + 0.5;\n",
       "    y0 = Math.floor(y0) + 0.5;\n",
       "    x1 = Math.floor(x1) + 0.5;\n",
       "    y1 = Math.floor(y1) + 0.5;\n",
       "    var min_x = Math.min(x0, x1);\n",
       "    var min_y = Math.min(y0, y1);\n",
       "    var width = Math.abs(x1 - x0);\n",
       "    var height = Math.abs(y1 - y0);\n",
       "\n",
       "    fig.rubberband_context.clearRect(\n",
       "        0, 0, fig.canvas.width, fig.canvas.height);\n",
       "\n",
       "    fig.rubberband_context.strokeRect(min_x, min_y, width, height);\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_figure_label = function(fig, msg) {\n",
       "    // Updates the figure title.\n",
       "    fig.header.textContent = msg['label'];\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_cursor = function(fig, msg) {\n",
       "    var cursor = msg['cursor'];\n",
       "    switch(cursor)\n",
       "    {\n",
       "    case 0:\n",
       "        cursor = 'pointer';\n",
       "        break;\n",
       "    case 1:\n",
       "        cursor = 'default';\n",
       "        break;\n",
       "    case 2:\n",
       "        cursor = 'crosshair';\n",
       "        break;\n",
       "    case 3:\n",
       "        cursor = 'move';\n",
       "        break;\n",
       "    }\n",
       "    fig.rubberband_canvas.style.cursor = cursor;\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_message = function(fig, msg) {\n",
       "    fig.message.textContent = msg['message'];\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_draw = function(fig, msg) {\n",
       "    // Request the server to send over a new figure.\n",
       "    fig.send_draw_message();\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_image_mode = function(fig, msg) {\n",
       "    fig.image_mode = msg['mode'];\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.updated_canvas_event = function() {\n",
       "    // Called whenever the canvas gets updated.\n",
       "    this.send_message(\"ack\", {});\n",
       "}\n",
       "\n",
       "// A function to construct a web socket function for onmessage handling.\n",
       "// Called in the figure constructor.\n",
       "mpl.figure.prototype._make_on_message_function = function(fig) {\n",
       "    return function socket_on_message(evt) {\n",
       "        if (evt.data instanceof Blob) {\n",
       "            /* FIXME: We get \"Resource interpreted as Image but\n",
       "             * transferred with MIME type text/plain:\" errors on\n",
       "             * Chrome.  But how to set the MIME type?  It doesn't seem\n",
       "             * to be part of the websocket stream */\n",
       "            evt.data.type = \"image/png\";\n",
       "\n",
       "            /* Free the memory for the previous frames */\n",
       "            if (fig.imageObj.src) {\n",
       "                (window.URL || window.webkitURL).revokeObjectURL(\n",
       "                    fig.imageObj.src);\n",
       "            }\n",
       "\n",
       "            fig.imageObj.src = (window.URL || window.webkitURL).createObjectURL(\n",
       "                evt.data);\n",
       "            fig.updated_canvas_event();\n",
       "            fig.waiting = false;\n",
       "            return;\n",
       "        }\n",
       "        else if (typeof evt.data === 'string' && evt.data.slice(0, 21) == \"data:image/png;base64\") {\n",
       "            fig.imageObj.src = evt.data;\n",
       "            fig.updated_canvas_event();\n",
       "            fig.waiting = false;\n",
       "            return;\n",
       "        }\n",
       "\n",
       "        var msg = JSON.parse(evt.data);\n",
       "        var msg_type = msg['type'];\n",
       "\n",
       "        // Call the  \"handle_{type}\" callback, which takes\n",
       "        // the figure and JSON message as its only arguments.\n",
       "        try {\n",
       "            var callback = fig[\"handle_\" + msg_type];\n",
       "        } catch (e) {\n",
       "            console.log(\"No handler for the '\" + msg_type + \"' message type: \", msg);\n",
       "            return;\n",
       "        }\n",
       "\n",
       "        if (callback) {\n",
       "            try {\n",
       "                // console.log(\"Handling '\" + msg_type + \"' message: \", msg);\n",
       "                callback(fig, msg);\n",
       "            } catch (e) {\n",
       "                console.log(\"Exception inside the 'handler_\" + msg_type + \"' callback:\", e, e.stack, msg);\n",
       "            }\n",
       "        }\n",
       "    };\n",
       "}\n",
       "\n",
       "// from http://stackoverflow.com/questions/1114465/getting-mouse-location-in-canvas\n",
       "mpl.findpos = function(e) {\n",
       "    //this section is from http://www.quirksmode.org/js/events_properties.html\n",
       "    var targ;\n",
       "    if (!e)\n",
       "        e = window.event;\n",
       "    if (e.target)\n",
       "        targ = e.target;\n",
       "    else if (e.srcElement)\n",
       "        targ = e.srcElement;\n",
       "    if (targ.nodeType == 3) // defeat Safari bug\n",
       "        targ = targ.parentNode;\n",
       "\n",
       "    // jQuery normalizes the pageX and pageY\n",
       "    // pageX,Y are the mouse positions relative to the document\n",
       "    // offset() returns the position of the element relative to the document\n",
       "    var x = e.pageX - $(targ).offset().left;\n",
       "    var y = e.pageY - $(targ).offset().top;\n",
       "\n",
       "    return {\"x\": x, \"y\": y};\n",
       "};\n",
       "\n",
       "/*\n",
       " * return a copy of an object with only non-object keys\n",
       " * we need this to avoid circular references\n",
       " * http://stackoverflow.com/a/24161582/3208463\n",
       " */\n",
       "function simpleKeys (original) {\n",
       "  return Object.keys(original).reduce(function (obj, key) {\n",
       "    if (typeof original[key] !== 'object')\n",
       "        obj[key] = original[key]\n",
       "    return obj;\n",
       "  }, {});\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.mouse_event = function(event, name) {\n",
       "    var canvas_pos = mpl.findpos(event)\n",
       "\n",
       "    if (name === 'button_press')\n",
       "    {\n",
       "        this.canvas.focus();\n",
       "        this.canvas_div.focus();\n",
       "    }\n",
       "\n",
       "    var x = canvas_pos.x;\n",
       "    var y = canvas_pos.y;\n",
       "\n",
       "    this.send_message(name, {x: x, y: y, button: event.button,\n",
       "                             step: event.step,\n",
       "                             guiEvent: simpleKeys(event)});\n",
       "\n",
       "    /* This prevents the web browser from automatically changing to\n",
       "     * the text insertion cursor when the button is pressed.  We want\n",
       "     * to control all of the cursor setting manually through the\n",
       "     * 'cursor' event from matplotlib */\n",
       "    event.preventDefault();\n",
       "    return false;\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._key_event_extra = function(event, name) {\n",
       "    // Handle any extra behaviour associated with a key event\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.key_event = function(event, name) {\n",
       "\n",
       "    // Prevent repeat events\n",
       "    if (name == 'key_press')\n",
       "    {\n",
       "        if (event.which === this._key)\n",
       "            return;\n",
       "        else\n",
       "            this._key = event.which;\n",
       "    }\n",
       "    if (name == 'key_release')\n",
       "        this._key = null;\n",
       "\n",
       "    var value = '';\n",
       "    if (event.ctrlKey && event.which != 17)\n",
       "        value += \"ctrl+\";\n",
       "    if (event.altKey && event.which != 18)\n",
       "        value += \"alt+\";\n",
       "    if (event.shiftKey && event.which != 16)\n",
       "        value += \"shift+\";\n",
       "\n",
       "    value += 'k';\n",
       "    value += event.which.toString();\n",
       "\n",
       "    this._key_event_extra(event, name);\n",
       "\n",
       "    this.send_message(name, {key: value,\n",
       "                             guiEvent: simpleKeys(event)});\n",
       "    return false;\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.toolbar_button_onclick = function(name) {\n",
       "    if (name == 'download') {\n",
       "        this.handle_save(this, null);\n",
       "    } else {\n",
       "        this.send_message(\"toolbar_button\", {name: name});\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.toolbar_button_onmouseover = function(tooltip) {\n",
       "    this.message.textContent = tooltip;\n",
       "};\n",
       "mpl.toolbar_items = [[\"Home\", \"Reset original view\", \"fa fa-home icon-home\", \"home\"], [\"Back\", \"Back to  previous view\", \"fa fa-arrow-left icon-arrow-left\", \"back\"], [\"Forward\", \"Forward to next view\", \"fa fa-arrow-right icon-arrow-right\", \"forward\"], [\"\", \"\", \"\", \"\"], [\"Pan\", \"Pan axes with left mouse, zoom with right\", \"fa fa-arrows icon-move\", \"pan\"], [\"Zoom\", \"Zoom to rectangle\", \"fa fa-square-o icon-check-empty\", \"zoom\"], [\"\", \"\", \"\", \"\"], [\"Download\", \"Download plot\", \"fa fa-floppy-o icon-save\", \"download\"]];\n",
       "\n",
       "mpl.extensions = [\"eps\", \"pdf\", \"png\", \"ps\", \"raw\", \"svg\"];\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,iVBORw0KGgoAAAANSUhEUgAAAjAAAAIwCAYAAACY8VFvAAAgAElEQVR4nO3babCuV1Xg8QcbxDJYKUtunZJUJee+77N3G2gRNNgO0RhplAgKjaDGGY0IOERRFBQIIMaAkBBlRkAwCFxoRIJBERMTzABKG8UhQFSChEFwwAHwg3n6Sy51Oga56F5nPfvs36/q/+Hk3qqz31u91ruqWqYJAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAJq5TSnl50spV5dSriilvPLQoUN32PsXttvtGaWUa0opl5dSXrvZbI7PeiwAwLTZbE4tpVx09Od5nl88z/Mj9vx8+1rrjbu7u7vTNE2llMeWUi5MeCoAwK369FLKb8/z/MCj/2Ge59NKKVcc/Xm73c6llOtzngcAsEcp5cm11veWUi64xX8/s9Z65OjPOzs7x5VSPrr/LwQAuHW3K6W8bO//F9ItD5hDhw7doZTykZznAQCr8fknfuUS0bH87nme73z48OG7Hv25lPKAUsrFR3/ebDan1lqv3PPnJ9dar2v7LwAAdCf5gHngzf83Lv9tmqap1vq0Wut5e/7K7UopN2y32/nmv3/uLf4cABjRXU86bYnoWH//PM/nllKuqbW+qZTyqs1mc3wp5YLNZnPKNE3Tdrs9/ej/zLrWemRnZ+e4uH8NAKAL2QcMAMCn7At2T18iyv5cAMAB5oABALrjgAEAunP33a9eIsr+XADAAeaAAQC684WH77lElP25AIADzAEDAHTnlM29loiyPxcAcIA5YACA7txj+zVLRNmfCwA4wBwwAEB3HDAAQHe+ePu1S0TZnwsAOMAcMABAd75kPmOJKPtzAQAHmAMGAOjOl5avWyLK/lwAwAH25fU+S0TZnwsAOMAcMABAdxwwAEB3HDAAQHdOrV+/RJT9uQCAA8wBAwB05yv++zcsEWV/LgDgADvt8+63RJT9uQCAA8wBAwB056s+7/5LRNmfCwA4wBwwAEB3HDAAQHccMABAd7765AcsEWV/LgDgALvnnb9xiSj7cwEAB5gDBgDozr3u8sAlouzPBQAcYA4YAKA7X3OXBy0RZX8uAOAAc8AAAN1xwAAA3bn3//jmJaLszwUAHGAOGACgO2d8/rcsEWV/LgDgAHPAAADduc9dz1wiyv5cAMAB5oABALpz37t+6xJR9ucCAA4wBwzD+tcPf2iRJMUWtcMdMAwre6glaYSidvjXf8G3LRFFvReayR5qSRqhqB3ugGFY2UMtSSMUtcPvd7dvXyKKei80kz3UkjRCUTvcAcOwsodakkYoaoff/27fsUQU9V5oJnuoJWmEonb4A+7+XUtEUe+FZrKHWpJGKGqHO2AYVvZQS9IIRe1wBwzDyh5qSRqhqB3ugGFY2UMtSSMUtcO/8Qu/e4ko6r3QTPZQS9IIRe1wBwzDyh5qSRqhqB3+oC968BJR1HuhmeyhlqQRitrh3/RF37NEFPVeaCZ7qCVphKJ2uAOGYWUPtSSNUNQO/5ZTvneJKOq90Ez2UEvSCEXtcAcMw8oeakkaoagd7oBhWNlDLUkjFLXDHTAMK3uoJWmEonb4mfc4a4ko6r3QTPZQS9IIRe3wb/vi71siOtbfP8/zo2qtbymlXFVrfck0Tbfd++ellJtqrZeWUi6rtV46z/ODmv8jMKbsoZakEYra4ZkHzHa7/bJSyrXTNH3aNE1TKeVVpZQH7/07pZR/C/jY4ICRpP0oaod/+xc/ZInoGH/9bXZ2do47+kMp5dmllEfu/QullJuafmA4KnuoJWmEonZ48gHzcdvtdi6lvHuz2Zy497+XUm4qpbyglHJFKeWiw4cP77T79Awte6glaYSidvh3/s/vXyL6VN5w+PDhu5ZS3rnZbO55yz8rpTz0hBNO+JxpmqZa6+NrrUdafXYGlz3UkjRCUTs8+4CZ5/nutda3bzabexzD371zKeUd/7VPDDfLHmpJGqGoHZ55wNzpTnf6zFrrdfM83+XW/ryUcnIp5XXTzf/LpFLK2bXWVzT8+Iwse6glaYSidvh3fclDl4iO5XfP83xWrfWDe/9n0rXWR9daz99sNqfc/HfOKaW8tZRyWSnl4pNOOulzo/4tGEz2UEvSCEXt8MwDBlJlD7UkjVDUDn/wlz5siSjqvdBM9lBL0ghF7XAHDMPKHmpJGqGoHf49X/rwJaKo90Iz2UMtSSMUtcMdMAwre6glaYSidvhZX/YDS0RR74VmsodakkYoaoc7YBhW9lBL0ghF7XAHDMPKHmpJGqGoHf59X/6DS0RR74VmsodakkYoaoc7YBhW9lBL0ghF7fDvP/WHloii3gvNZA+1JI1Q1A53wDCs7KGWpBGK2uEP/YofXiKKei80kz3UkjRCUTv8YV9x9hJR1HuhmeyhlqQRitrhDhiGlT3UkjRCUTvcAcOwsodakkYoaoc7YBhW9lBL0ghF7fAf+Mqzl4ii3gvNZA+1JI1Q1A53wDCs7KGWpBGK2uE/eNqPLBFFvReayR5qSRqhqB3+w1/1o0tEUe+FZrKHWpJGKGqHO2AYVvZQS9IIRe3ws09/xBJR1HuhmeyhlqQRitrhDhiGlT3UkjRCUTvcAcOwsodakkYoaoc7YBhW9lBL0ghF7fAfPf3Hloii3gvNZA+1JI1Q1A5/xFf/+BJR1HuhmeyhlqQRitrhDhiGlT3UkjRCUTv8x+/5yCWiqPdCM9lDLUkjFLXDHTAMK3uoJWmEonb4I//XI5eIot4LzWQPtSSNUNQOd8AwrOyhlqQRitrhDhiGlT3UkjRCUTv8J+/1E0tEUe+FZrKHWpJGKGqHO2AYVvZQS9IIRe3wR93rJ5eIot4LzWQPtSSNUNQOd8AwrOyhlqQRitrhj/6aRy0RRb0Xmskeakkaoagd7oBhWNlDLUkjFLXDf/prH71EFPVeaCZ7qCVphKJ2+GPu/egloqj3QjPZQ60++tPnv1wrL/v/jeg/LmqHO2AYVvZQq4+yv5zlgOm9qB3ugGFY2UOtPsr+cpYDpveidrgDhmFlD7X6KPvLWQ6Y3ova4Y8746eWiKLeC81kD7X6KPvLWQ6Y3ova4Q4YhpU91Oqj7C9nOWB6L2qHn3PGTy8RRb0XmskeavVR9pezHDC9F7XDn/B1j1kiinovNJM91Oqj7C9nOWB6L2qHO2AYVvZQq4+yv5zlgOm9qB3+xPs8Zoko6r3QTPZQq4+yv5zlgOm9qB3ugGFY2UOtPsr+cpYDpveidrgDhmFlD7X6KPvLWQ6Y3ova4Q4YhpU91Oqj7C9nOWB6L2qH/8x9H7tEFPVeaCZ7qNVH2V/OcsD0XtQO/9mvf9wSUdR7oZnsoVYfZX85ywHTe1E73AHDsLKHWn2U/eUsB0zvRe3wc7/hcUtEUe+FZrKHWn2U/eUsB0zvRe1wBwzDyh5q9VH2l7McML0XtcPPu985S0RR74VmsodafZT95SwHTO9F7XAHDMPKHmr1UfaXsxwwvRe1wx0wDCt7qNVH2V/OcsD0XtQOf/L9Hr9EFPVeaCZ7qNVH2V/OcsD0XtQOd8AwrOyhVh9lfznLAdN7UTv8Kfd/whJR1HuhmeyhVh9lfznLAdN7UTvcAcOwsodafZT95SwHTO9F7fCn/u8nLBEd6++f5/lRtda3lFKuqrW+ZJqm2+798+12e0Yp5ZpSyuWllNduNpvjm/8jMKbsoVYfZX85ywHTe1E7PPOA2W63X1ZKuXaapk+bpmkqpbyqlPLgo38+z/Pta6037u7u7t78548tpVwY8e/AgLKHWn2U/eUsB0zvRe3wpz3giUtEx/jrb7Ozs3Pc0R9KKc8upTzy6M/zPJ9WSrni6M/b7XYupVzf8OMzsuyhVh9lfznLAdN7UTs8+YD5uJuPk3dvNpsTj/63UsqZtdYjR3/e2dk5rpTy0VafncFlD7X6KPvLWQ6Y3ova4Ws4YA4fPnzXUso7N5vNPff+91seMIcOHbpDKeUjrT47g8seavVR9pezHDC9F7XDz//Gn1kiOtbfP8/z3Wutb99sNve45Z9tNptTa61XHv25lHJyrfW6Vp+dwWUPtfoo+8tZDpjei9rhmQfMne50p8+stV43z/NdPsFfuV0p5YbtdjtP0zTN83xurfW8Zh+esWUPtfoo+8tZDpjei9rhT3/gzywRHcvvnuf5rFrrB2utl5ZSLqu1XlprfXSt9fzNZnPKNE3Tdrs9vZRydSnlilrrkb3/R7/wX5I91Oqj7C9nOWB6L2qHZx4wkCp7qNVH2V/OcsD0XtQOv/BBT1oiinovNJM91Oqj7C9nOWB6L2qH/8I3PWmJKOq90Ez2UK+x7C8iSW3K3iUOGAiUPdRrLHvpSmpT9i5xwECg7KFeY9lLV1KbsneJAwYCZQ/1GsteupLalL1L9uOAecY3/+wSUdR7oZnsoV5j2UtXUpuyd4kDBgJlD/Uay166ktqUvUv244B55jefu0QU9V5oJnuo11j20pXUpuxdsh8HzLO+5dwloqj3QjPZQ73GspeupDZl7xIHDATKHuo1lr10JbUpe5fsxwHznDN/boko6r3QTPZQr7HspSupTdm7xAEDgbKHeo1lL11JbcreJQ4YCJQ91Gsse+lKalP2LnHAQKDsoV5j2UtXUpuyd8l+HDDP/dbzloii3gvNZA/1GsteupLalL1L9uOAed63nbdEFPVeaCZ7qNdY9tKV1KbsXeKAgUDZQ73GspeupDZl75L9OGCe/+1PXiKKei80kz3Uayx76UpqU/YuccBAoOyhXmPZS1dSm7J3yX4cMC/4jicvEUW9F5rJHuo1lr10JbUpe5c4YCBQ9lCvseylK6lN2bvEAQOBsod6jWUvXUltyt4l+3HAvPA7n7JEFPVeaCZ7qNdY9tKV1KbsXeKAgUDZQ73GspeupDZl75L9OGBe9F1PWSKKei80kz3Uayx76UpqU/YuccBAoOyhXmPZS1dSm7J3yX4cMC/+7p9fIop6LzSTPdRrLHvpSmpT9i5xwECg7KFeY9lLV1KbsnfJfhwwL3nwU5eIot4LzWQP9RrLXrqS2pS9SxwwECh7qNdY9tKV1KbsXeKAgUDZQ73GspeupDZl75L9OGB+5XueukQU9V5oJnuo11j20pXUpuxd4oCBQNlDvcayl66kNmXvkv04YC763qctEUW9F5rJHuo1lr10JbUpe5c4YCBQ9lCvseylK6lN2btkPw6Yl571tCWiqPdCM9lDvcayl66kNmXvkv04YF72fecvEUW9F5rJHuo1lr109/b5J36l1FXZM+OAccAwiOyhXmPZS9cBo57LnhkHjAOGQWQP9RrLXroOGPVc9sw4YBwwDCJ7qNdY9tJ1wKjnsmdmtAPm5Q85f4ko6r3QTPZQr7HspeuAUc9lz4wDxgHDILKHeo1lL10HjHoue2ZGO2Be8f0XLBFFvReayR7qNZa9dB0w6rnsmRntgDny0AuWiKLeC81kD/Uay166Dhj1XPbMOGAcMAwie6jXWPbSdcCo57JnZrQD5pUPe/oSUdR7oZnsoV5j2UvXAaOey54ZB4wDhkFkD/Uay166Dhj1XPbMOGAcMAwie6jXWPbSdcCo57JnxgHjgGEQ2UO9xrKXrgNGPZc9M6MdMK96+NOXiKLeC81kD/Uay166Dhj1XPbMjHbAvPoHLlwiinovNJM91Gsse+k6YNRz2TPjgHHAMIjsoV5j2UvXAaOey56Z0Q6YX/vBC5eIot4LzWQP9RrLXroOGPVc9sw4YBwwDCJ7qNdY9tJ1wKjnsmdmtAPmNT/0C0tEUe+FZrKHeo1lL10HjHoue2YcMA4YBpE91Gsse+k6YNRz2TPjgHHAMIjsoV5j2UvXAaOey56Z0Q6YX//hX1giinovNJM91Gsse+k6YNRz2TPjgHHAMIjsoV5j2UvXAaOey56Z0Q6Y1579i0tEUe+FZrKHeo1lL10HjHoue2YcMA4YBpE91Gsse+k6YNRz2TMz2gFz8Y/84hJR1HuhmeyhXmPZS9cBo57LnhkHjAOGQWQP9RrLXroOGPVc9syMdsC87kefsUQU9V5oJnuo11j20nXAqOeyZ8YB44BhENlDvcayl64DRj2XPTMOGAcMg8ge6jWWvXQdMOq57JkZ7YD5jUc8Y4noWH//ZrM5vtb68lrr+27tz0spN9VaLy2lXFZrvXSe5we1+/QMLXuo11j20nXAqOeyZ8YBs78HTK31klrrQ2qt7721Py+l/Fu7Twt7ZA/1Gsteug4Y9Vz2zIx2wFzyY89cIjrW33/HO97xsw4fPnzSf3DA3NTu08Ie2UO9xrKXrgNGPZc9Mw6Y/T1gpmmaPtkBU0p5QSnlilLKRYcPH95p88kZXvZQr7HspeuAUc9lz8xoB8zrf/yZS0Sfyhs+yQHz0BNOOOFzpmmaaq2Pr7UeafG5wQHjgJGalj0zox0wv/XIZy0RfSpv+I8OmL3meb5zKeUd//lPC3tkD/Uay166Dhj1XPbMOGD2/4DZ3d3dvbX/FVIp5eRSyuumabrtzT+fXWt9RaOPzuiyh3qNZS9dB4x6LntmRjtg3vATz1oiOpbffeKJJ352KeWyUspVpZSP3fw/l76w1nr+ZrM5ZZqmaZ7nc0opb73571180kknfW7UvwWDyR7qNZa9dB0w6rnsmXHA7N8BA6myh3qNZS9dB4x6LntmHDAOGAaRPdRrLHvpOmDUc9kz44BxwDCI7KFeY9lL1wGjnsuemdEOmN/+yWcvEUW9F5rJHuo1lr10HTDqueyZGe2AeeOjnr1EFPVeaCZ7qNdY9tJ1wKjnsmfGAeOAYRDZQ73GspeuA0Y9lz0zox0wv/NTz1kiinovNJM91Gsse+k6YNRz2TPjgHHAMIjsoV5j2UvXAaOey56Z0Q6YS3/6OUtEUe+FZrKHeo1lL10HjHoue2YcMA4YBpE91Gsse+k6YNRz2TPjgHHAMIjsoV5j2UvXAaOey56Z0Q6Yyx7z3CWiqPdCM9lDvcayl64DRj2XPTMOGAcMg8ge6jWWvXQdMOq57JkZ7YD53cc+d4ko6r3QTPZQr7HspeuAUc9lz4wDxgHDILKHeo1lL10HjHoue2ZGO2Auf9zzloii3gvNZA/1Gsteug4Y9Vz2zDhgHDAMInuo11j20nXAqOeyZ2a0A+aKc563RBT1Xmgme6jXWPbSdcCo57JnxgHjgGEQ2UO9xrKXrgNGPZc9Mw4YBwyDyB7qNZa9dB0w6rnsmRnugHn885eIot4LzWQP9RrLXroOGPVc9sw4YBwwDCJ7qNdY9tJ1wKjnsmdmtAPmTU/4pSWiqPdCM9lDvcayl64DRj2XPTMOGAcMg8ge6jWWvXQdMOq57JkZ7YD5vSf+0hJR1HuhmeyhXmPZS9cBo57LnpnRDpgrn/SCJaKo90Iz2UO9xrKXrgNGPZc9Mw4YBwyDyB7qNZa9dB0w6rnsmXHAOGAYRPZQr7HspeuAUc9lz4wDxgHDILKHeo1lL10HjHoue2ZGO2Cu+tkXLBFFvReayR7qNZa9dB0w6rnsmXHAOGAYRPZQr7HspeuAUc9lz8xoB8zV575wiSjqvdBM9lCvseylK6lN2btkPw6Ya37uhUtEUe+FZrKHeo1lL11JbcreJQ4YCJQ91Gsse+lKalP2LtmPA+bN571oiSjqvdBM9lCvseylK6lN2bvEAQOBsod6jWUvXUltyt4lDhgIlD3Uayx76UpqU/YuccBAoOyhXmPZS1dSm7J3yX4cMG958ouWiKLeC81kD/Uay166ktqUvUv244D5/af88hJR1HuhmeyhXmPZS1dSm7J3iQMGAmUP9RrLXrqS2pS9S/blgHnqLy8RRb0Xmske6jWWvXQltSl7lzhgIFD2UK+x7KUrqU3Zu2Q/Dpg/eNqLl4ii3gvNZA/1GsteupLalL1LHDAQKHuo11j20pXUpuxd4oCBQNlDvcayl66kNmXvkv04YN56/ouXiKLeC81kD/Uay166ktqUvUscMBAoe6jXWPbSldSm7F2yHwfM/73gJUtEUe+FZrKHeo1lL11JbcreJQ4YCJQ91Gsse+lKalP2LtmPA+YPn/6SJaKo90Iz2UO9xrKXrqQ2Ze8SBwwEyh7qNZa9dCW1KXuX7MsBc+GvLBFFvReayR7qNZa9dCW1KXuXOGAgUPZQr7HspSupTdm7xAEDgbKHeo1lL11JbcreJftxwFz7ixctEUW9F5rJHuo1lr10JbUpe5c4YCBQ9lCvseylK6lN2btkPw6YP3rGRUtEUe+FZrKHeo1lL11JbcreJQ4YCJQ91Gsse+lKalP2LtmPA+aPn/nSJaKo90Iz2UO9xrKXrqQ2Ze+S/Thg3vasly4RRb0Xmske6jWWvXQltSl7lzhgIFD2UK+x7KUrqU3Zu8QBA4Gyh3qNZS9dSW3K3iUOGAiUPdRrLHvpSmpT9i7ZlwPm2b+6RBT1Xmgme6jXWPbSldSm7F3igIFA2UO9xrKXrqQ2Ze+S/Thg/uS5L1siinovNJM91Oqj7C8i9fVlrX9f1A7/0+e9bInoWH//ZrM5vtb68lrr+27tz7fb7RmllGtKKZeXUl672WyOb/fpGVr2UKuPsr+c5YDpvagdnn3A1FovqbU+pNb63lv+2TzPt6+13ri7u7s7TdNUSnlsKeXChh+fkWUPtfoo+8tZDpjei9rhUf/v6Vh//x3veMfPOnz48Emf4IA5rZRyxdGft9vtXEq5vtVnZ3DZQ60+yv5ylgOm96J2ePYBM03T9IkOmFLKmbXWI0d/3tnZOa6U8tEWnxscMDqmsr+c5YDpvagd3tMBc+jQoTuUUj7S4nODA0bHVPaXsxwwvRe1w9d8wGw2m1NrrVce/bmUcnKt9boWnxscMDqmsr+c5YDpvagd/me/9PIlok/lDbu7u7uf4H+FdLtSyg3b7Xaepmma5/ncWut5TT44ZA+1+ij7y1kOmN6L2uF//sJXLBEdy+8+8cQTP7uUclkp5apSysdqrZeWUi6stZ6/2WxOmaZp2m63p5dSri6lXFFrPbKzs3Nc1L8Fg8keavVR9pezHDC9F7XDMw8YSJU91Oqj7C9nOWB6L2qHX/eiI0tEUe+FZrKHWn2U/eUsB0zvRe1wBwzDyh5q9VH2l7McML0XtcOv++UjS0RR74VmsodafZT95SwHTO9F7XAHDMPKHmr1UfaXsxwwvRe1wx0wDCt7qNVH2V/OcsD0XtQOf/uLX7lEFPVeaCZ7qNVH2V/OcsD0XtQOd8AwrOyhVh9lfznLAdN7UTv8HS951RJR1HuhmeyhVh9lfznLAdN7UTvcAcOwsodafZT95SwHTO9F7fB3XPSqJaKo90Iz2UOtPsr+cpYDpveidrgDhmFlD7X6KPvLWQ6Y3ova4e986f9ZIop6LzSTPdTqo+wvZzlgei9qhztgGFb2UKuPsr+c5YDpvagd7oBhWNlDrT7K/nKWA6b3onb49b/66iWiqPdCM9lDrT7K/nKWA6b3ona4A4ZhZQ+1+ij7y1kOmN6L2uHXv+zVS0RR74VmsodafZT95SwHTO9F7XAHDMPKHmr1UfaXsxwwvRe1w//iFb+2RBT1Xmgme6jVR9lfznLA9F7UDv+LI69ZIop6LzSTPdTqo+wvZzlgei9qhztgGFb2UKuPsr+c5YDpvagd7oBhWNlDrT7K/nKWA6b3ona4A4ZhZQ+1+ij7y1kOmN6L2uF/+crXLBFFvReayR5q9VH2l7McML0XtcMdMAwre6jVR9lfznLA9F7UDv+rV/36ElHUe6GZ7KGWpBGK2uF/9erXLhFFvReayR5qSRqhqB3ugGFY2UMtSSMUtcPf9WsXLxFFvReayR5qSRqhqB3ugGFY2UMtSSMUtcMdMAwre6glaYSidrgDhmFlD7UkjVDUDn/Xay5eIop6LzSTPdSSNEJRO/yGX3/dElHUe6GZ7KGWpBGK2uEOGIaVPdSSNEJRO/zdF//GElHUe6GZ7KGWpBGK2uEOGIaVPdSSNEJRO/zdr7tkiSjqvdBM9lBL0ghF7XAHDMPKHmpJGqGoHf7Xv/H6JaKo90Iz2UMtSSMUtcMdMAwre6glaYSidrgDhmFlD7UkjVDUDv/rS16/RBT1Xmgme6glaYSidrgDhmFlD7UkjVDUDn/Pb/7mElHUe6GZ7KGWpBGK2uEOGIaVPdSSNEJRO/w9v/VbS0RR74VmsodakkYoaoff+IY3LBFFvReayR5qSRqhqB3ugGFY2UMtSSMUtcMdMAwre6glaYSidrgDhmFlD7UkjVDUDr/xjb+9RBT1Xmgme6glaYSidrgDhmFlD7UkjVDUDn/v77xxiSjqvdBM9lBL0ghF7fD3XvrGJaKo90Iz2UMtSSMUtcMdMAwre6glaYSidvj7LvudJaKo90Iz2UMtSSMUtcMdMAwre6glaYSidrgDhmFlD7UkjVDUDnfAMKzsoZakEYra4e+7/NIloqj3QjPZQy1JIxS1w99/xWVLRFHvhWayh1qSRihqhztgGFb2UEvSCEXt8Pe/6XeXiKLeC81kD7UkjVDUDnfAMKzsoZakEYra4R/4vcuXiI7ld9daH11KeXOt9cpSygW3/PNSyk211ktLKZfVWi+d5/lB7f8FGFb2UEvSCEXt8KwDZrPZ3KOUcu00TZ8+TdNtSilv2G6399/7d0op/xb1ucEBI0n7UNQOzzpg5nk+p5TyxD0/n1Vrff7ev1NKuSniM8M0TQ4YSdqPonb4B666Yonok/3eUspz53l++J6f71Nrff0t/s5NpZQXlFKuKKVcdPjw4Z2IfwMGlT3UkjRCUTt8RQfMfWutl9zi7zz0hBNO+JxpmqZa6+NrrUfa/wswrOyhlqQRitrhf3P1m5aIPtnvLaU8ppTypKM/z/P8sFLKcz7R35/n+c6llHe0+tzggJGkfShqh9+cSv8AAAUESURBVCceMHcrpbxtd3f3M6Zpum0p5fLtdnvvPX9+cinlddM03fbmn8+utb4i6t+BAWUPtSSNUNQO/5trfm+J6Fh+dynl7FLKm0spV9Van3Dzf7tgs9mcMk0f/z/0fWsp5bJSysUnnXTS50b9OzCg7KGWpBGK2uGZBwykyh5qSRqhqB3+wbdcuUQU9V5oJnuoJWmEona4A4ZhZQ+1JI1Q1A53wDCs7KGWpBGK2uEf/P2rloii3gvNZA+1JI1Q1A53wDCs7KGWpBGK2uEf+oOrl4ii3gvNZA+1JI1Q1A53wDCs7KGWpBGK2uEfeus1S0RR74VmsodakkYoaof/7R++eYko6r3QTPZQS9IIRe1wBwzDyh5qSRqhqB3ugGFY2UMtSSMUtcMdMAwre6glaYSidvjfXvuWJaKo90Iz2UMtSSMUtcMdMAwre6glaYSidvjf/tHvLxFFvReayR5qSRqhqB3+d3/8B0tEUe+FZrKHWpJGKGqHO2AYVvZQS9IIRe3wv/uTty4RRb0Xmskeakkaoagd7oBhWNlDLUkjFLXDHTAMK3uoJWmEona4AwYA6M7f/+kfLhFlfy4A4AD7+z+7doko+3MBAAeYAwYA6M4//PkfLRFlfy4A4ABzwAAA3fmHt//xElH25wIADjAHDADQHQcMANCdD7/jbUtE2Z8LADjAHDAAQHc+/M4/WSLK/lwAwAHmgAEAuvPh6/9siSj7cwEAB5gDBgDozj/+5Z8vEWV/LgDgAHPAAADdccAAAN35x7+6boko+3MBAAeYAwYA6M4/vevtS0TZnwsAOMAcMABAd/7phncuEWV/LgDgAPvnv75+iSj7cwEAB5gDBgDojgMGAOiOAwYA6M4/v+cvloiyPxcAcIA5YACA7vzLjX+5RJT9uQCAA+xf3vuuJaLszwUAHGAOGACgO//y/huWiLI/FwBwgDlgAIDuOGAAgO44YACA7nzkA+9eIsr+XADAAfaRv3nPElH25wIADjAHDADQnY9+8MYlouzPBQAcYA4YAKA7H/3Q+5aIsj8XAHCAOWAAgO44YACA7nzs796/RJT9uQCAA8wBAwB052N//4ElomP53bXWR5dS3lxrvbKUcsEt/3y73Z5RSrmmlHJ5KeW1m83m+Pb/AgBAd7IOmM1mc49SyrXTNH36NE23KaW8Ybvd3v/on8/zfPta6427u7u70zRNpZTHllIuDPuHAAD68a//8MElok/2e+d5PqeU8sQ9P59Va33+np9PK6VccfTn7XY7l1Kub/8vAAB0J+uAKaU8d57nh+/5+T611tfv+fnMWuuRoz/v7OwcV0r5aPt/AQCgO//64Q8tEX2y33srB8x9a62X7Pn5/ztgDh06dIdSykfa/wsAAByjUspjSilPOvrzPM8PK6U85+jPm83m1FrrlXv+/sm11uv2+50AAB9XSrlbKeVtu7u7nzFN021LKZdvt9t77/krtyul3LDdbudpmqZ5ns+ttZ6X81oAgJuVUs4upby5lHJVrfUJN/+3CzabzSnTNE3b7fb0UsrVpZQraq1HdnZ2jst9MQAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAD/zv8DL3B3vpjy3DUAAAAASUVORK5CYII=\">"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "biomasa_celda = medidas_por_celda(poblacion, -1, biomasa_min)[0]\n",
    "riqueza_celda = medidas_por_celda(poblacion, -1, biomasa_min)[1]\n",
    "\n",
    "fig = plt.figure(figsize=(7,7))\n",
    "seaborn.heatmap(riqueza_celda, yticklabels=False,\n",
    "               xticklabels=False, vmin = 0.0, annot=True, fmt=\".1f\", square=True) \n",
    "\n",
    "#seaborn.heatmap(poblacion[-2, :, :, 1], yticklabels=False,\n",
    "#               xticklabels=False, vmin = 0.0, annot=True, fmt=\".1f\", square=True) \n",
    "plt.show() "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([[   0.50606321,    0.57567066,    0.75465859,    0.97389628,\n",
       "           1.14783573,    1.21434742,    1.14783573,    0.97389628,\n",
       "           0.75465859,    0.57567066],\n",
       "       [   0.86320076,    1.04976219,    1.53448155,    2.11756442,\n",
       "           2.57851362,    2.75016445,    2.57851362,    2.11756442,\n",
       "           1.53448155,    1.04976219],\n",
       "       [   1.56108911,    2.30455896,    4.01088187,    6.16620143,\n",
       "           7.67377767,    8.30377385,    7.67377767,    6.16620143,\n",
       "           4.01088187,    2.30455896],\n",
       "       [   2.41306739,    4.10400462,   10.1984485 ,   16.95695822,\n",
       "          22.12408321,   23.11952831,   22.12408321,   16.95695822,\n",
       "          10.1984485 ,    4.10400462],\n",
       "       [   3.01578709,    5.56444052,   15.8098524 ,  174.50555664,\n",
       "         198.17070759,  200.48874982,  198.17070759,  174.50555664,\n",
       "          15.8098524 ,    5.56444052],\n",
       "       [   3.01578709,    5.56444052,   15.8098524 ,  174.50555664,\n",
       "         198.17070759,  200.48874982,  198.17070759,  174.50555664,\n",
       "          15.8098524 ,    5.56444052],\n",
       "       [   2.41306739,    4.10400462,   10.1984485 ,   16.95695822,\n",
       "          22.12408321,   23.11952831,   22.12408321,   16.95695822,\n",
       "          10.1984485 ,    4.10400462],\n",
       "       [   1.56108911,    2.30455896,    4.01088187,    6.16620143,\n",
       "           7.67377767,    8.30377385,    7.67377767,    6.16620143,\n",
       "           4.01088187,    2.30455896],\n",
       "       [   0.86320076,    1.04976219,    1.53448155,    2.11756442,\n",
       "           2.57851362,    2.75016445,    2.57851362,    2.11756442,\n",
       "           1.53448155,    1.04976219],\n",
       "       [   0.50606321,    0.57567066,    0.75465859,    0.97389628,\n",
       "           1.14783573,    1.21434742,    1.14783573,    0.97389628,\n",
       "           0.75465859,    0.57567066]])"
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "biomasa_celda"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array([[  15.49003116,   14.6359444 ,   11.6140351 ,    5.5259693 ,\n",
       "           2.79645506,    2.79645506,    5.5259693 ,   11.6140351 ,\n",
       "          14.6359444 ,   15.49003116],\n",
       "       [  19.34824593,  152.4156092 ,   15.4097418 ,    6.46276533,\n",
       "           3.56008371,    3.56008371,    6.46276533,   15.4097418 ,\n",
       "         152.4156092 ,   19.34824593],\n",
       "       [  18.68193668,  151.27762938,   15.8567421 ,    8.29613437,\n",
       "           6.59802872,    6.59802872,    8.29613437,   15.8567421 ,\n",
       "         151.27762938,   18.68193668],\n",
       "       [  12.26233819,   12.09899826,   11.60349213,   11.98267135,\n",
       "          15.21784467,   15.21784467,   11.98267135,   11.60349213,\n",
       "          12.09899826,   12.26233819],\n",
       "       [   5.96749494,    6.03673085,    7.09049507,   12.28132359,\n",
       "         149.25151631,  149.25151631,   12.28132359,    7.09049507,\n",
       "           6.03673085,    5.96749494],\n",
       "       [   5.94068144,    5.93896323,    6.46728098,   10.70900551,\n",
       "          14.92791147,   14.92791147,   10.70900551,    6.46728098,\n",
       "           5.93896323,    5.94068144],\n",
       "       [  12.20876506,   11.87603857,   10.62387976,    6.93666901,\n",
       "           6.24725518,    6.24725518,    6.93666901,   10.62387976,\n",
       "          11.87603857,   12.20876506],\n",
       "       [  18.63380888,  150.74166231,   14.92594207,    6.14245045,\n",
       "           3.3562364 ,    3.3562364 ,    6.14245045,   14.92594207,\n",
       "         150.74166231,   18.63380888],\n",
       "       [  19.32423843,  152.11135275,   15.08306349,    5.85039891,\n",
       "           2.74547877,    2.74547877,    5.85039891,   15.08306349,\n",
       "         152.11135275,   19.32423843],\n",
       "       [  15.4830273 ,   14.60186188,   11.52463267,    5.36555352,\n",
       "           2.58840879,    2.58840879,    5.36555352,   11.52463267,\n",
       "          14.60186188,   15.4830273 ]])"
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "biomasa_celda"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "389.62040827576391"
      ]
     },
     "execution_count": 20,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "sum(condiciones_iniciales[77])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.5.2"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 0
}