tri_inh_exc_stdp_jpnb2.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Tripartite with excitatory STDP"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "On the excitatory synapses"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 69,
   "metadata": {},
   "outputs": [],
   "source": [
    "from brian2 import *\n",
    "%matplotlib notebook"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 70,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Constants\n",
    "tau = 10*ms\n",
    "I = 1.1/ms\n",
    "gamma = 1/ms\n",
    "\n",
    "# Non-dimensionalized constants\n",
    "tp = tau/ms\n",
    "curr = I*ms\n",
    "g = gamma*ms\n",
    "out = curr/g  #should be > the threshold of firing\n",
    "\n",
    "# equation (1)\n",
    "eqs = '''\n",
    "dv/dt = I - gamma*v : 1\n",
    "'''\n",
    "\n",
    "# equation (2) \n",
    "eqs2 = '''\n",
    "dv/dt = ((I/gamma) - v)/tau : 1\n",
    "'''"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 71,
   "metadata": {},
   "outputs": [],
   "source": [
    "start_scope()\n",
    "\n",
    "n_i = 90 # 3 groups of n_i/3\n",
    "n_e = 30\n",
    "\n",
    "seed(21)\n",
    "\n",
    "Gi = NeuronGroup(n_i, eqs2, threshold='v>1', reset='v = 0', method='exact')\n",
    "Ge = NeuronGroup(n_e, eqs2, threshold='v>1', reset='v = 0', method='exact')\n",
    "\n",
    "group1init = randint(0, 9, size=n_i)/10\n",
    "group2init = randint(0, 9, size=n_e)/10\n",
    "\n",
    "Gi.v = group1init\n",
    "Ge.v = group2init"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 72,
   "metadata": {},
   "outputs": [],
   "source": [
    "matrix = zeros((n_i, n_i))\n",
    "\n",
    "def matcon(i,j):\n",
    "        matrix[i, j] = 1\n",
    "\n",
    "# each group of 30 to everything else\n",
    "for i in range(n_i//3):\n",
    "        for j in range(n_i//3, n_i):\n",
    "                matcon(i,j)\n",
    "for i in range(n_i//3, 2*n_i//3):\n",
    "        for j in range(2*n_i//3, n_i):\n",
    "                matcon(i,j)\n",
    "        for j in range(0, n_i//3):\n",
    "                matcon(i,j)\n",
    "for i in range(2*n_i//3, n_i):\n",
    "        for j in range(0, 2*n_i//3):\n",
    "                matcon(i,j)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 73,
   "metadata": {},
   "outputs": [],
   "source": [
    "# imshow(matrix, origin='lower');"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 74,
   "metadata": {},
   "outputs": [],
   "source": [
    "inh = 0.1\n",
    "exc = 0.01\n",
    "\n",
    "p = 0.8\n",
    "exc_matrix = choice(2, (n_e, n_i), p=[1-p, p])\n",
    "inh_matrix = exc_matrix.T\n",
    "\n",
    "\n",
    "\n",
    "si, ti = matrix.nonzero()\n",
    "se, te = exc_matrix.nonzero()\n",
    "si2, ti2 = inh_matrix.nonzero()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 75,
   "metadata": {},
   "outputs": [],
   "source": [
    "# imshow(exc_matrix, origin='lower');"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 76,
   "metadata": {},
   "outputs": [],
   "source": [
    "# imshow(inh_matrix, origin='lower');"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 77,
   "metadata": {},
   "outputs": [],
   "source": [
    "tau_stdp = 100*ms\n",
    "# gmax = .01\n",
    "# dApre = .01\n",
    "# dApost = -dApre * taupre / taupost * 1.05\n",
    "# dApost *= gmax\n",
    "# dApre *= gmax\n",
    "\n",
    "gmax = 1.0\n",
    "\n",
    "\n",
    "I1 = Synapses(Gi, Gi,\n",
    "              '''w : 1\n",
    "              dApre/dt = -Apre/tau_stdp : 1 (event-driven)\n",
    "              dApost/dt = -Apost/tau_stdp : 1 (event-driven)\n",
    "              ''', \n",
    "              on_pre='''Apre += 0.4\n",
    "              w = clip(w + Apost, 0, gmax)\n",
    "              v = out * (1-exp(-(tp * log(1/(1 - g*(v_post - inh*w)/curr)))/tp))''',\n",
    "              on_post='''Apost -= 0.4\n",
    "              w = clip(w + Apre, 0, gmax)\n",
    "              ''')\n",
    "\n",
    "# Excitatory synapses\n",
    "E1 = Synapses(Ge, Gi,\n",
    "              '''w : 1\n",
    "              dApre/dt = -Apre/tau_stdp : 1 (event-driven)\n",
    "              dApost/dt = -Apost/tau_stdp : 1 (event-driven)\n",
    "              ''',\n",
    "              on_pre='''Apre += 0.4\n",
    "              w = clip(w + Apost, 0, gmax)\n",
    "              v = out * (1-exp(-(tp * log(1/(1 - g*(v_post + exc*w)/curr)))/tp))''',\n",
    "              on_post='''Apost -= 0.4\n",
    "              w = clip(w + Apre, 0, gmax)\n",
    "              ''')\n",
    "\n",
    "# Inhibitory recurrent synapses\n",
    "I2 = Synapses(Gi, Ge,\n",
    "              '''w : 1\n",
    "              dApre/dt = -Apre/tau_stdp : 1 (event-driven)\n",
    "              dApost/dt = -Apost/tau_stdp : 1 (event-driven)\n",
    "              ''', \n",
    "              on_pre='''Apre += 0.4\n",
    "              w = clip(w + Apost, 0, gmax)\n",
    "              v = out * (1-exp(-(tp * log(1/(1 - g*(v_post - inh*w)/curr)))/tp))''',\n",
    "              on_post='''Apost -= 0.4\n",
    "              w = clip(w + Apre, 0, gmax)\n",
    "              ''')\n",
    "\n",
    "\n",
    "I1.connect(i=si, j=ti)\n",
    "\n",
    "E1.connect(i=se, j=te)\n",
    "\n",
    "I2.connect(i=si2, j=ti2)\n",
    "\n",
    "M1 = StateMonitor(Gi, 'v', record=True)\n",
    "Sp1 = SpikeMonitor(Gi)\n",
    "WmonI1 = StateMonitor(I1, 'w', record=True)\n",
    "# WmonE1 = StateMonitor(E1, 'w', record=True)\n",
    "# WmonI2 = StateMonitor(I2, 'w', record=True)\n",
    "I1.w = 'rand() * gmax'\n",
    "E1.w = 'rand() * gmax'\n",
    "I2.w = 'rand() * gmax'"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 78,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "WARNING    Came across an abstract code block that may not be well-defined: the outcome may depend on the order of execution. You can ignore this warning if you are sure that the order of operations does not matter. 8 lines of abstract code, first line is: '_Apost := Apost * exp(_lio_1 * (- (t - lastupdate))) (constant)'\n",
      " [brian2.codegen.generators.base]\n",
      "WARNING    Came across an abstract code block that may not be well-defined: the outcome may depend on the order of execution. You can ignore this warning if you are sure that the order of operations does not matter. 8 lines of abstract code, first line is: '_Apost := Apost * exp(_lio_1 * (- (t - lastupdate))) (constant)'\n",
      " [brian2.codegen.generators.base]\n",
      "WARNING    Came across an abstract code block that may not be well-defined: the outcome may depend on the order of execution. You can ignore this warning if you are sure that the order of operations does not matter. 8 lines of abstract code, first line is: '_Apost := Apost * exp(_lio_1 * (- (t - lastupdate))) (constant)'\n",
      " [brian2.codegen.generators.base]\n"
     ]
    }
   ],
   "source": [
    "run(800*ms)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 79,
   "metadata": {},
   "outputs": [],
   "source": [
    "# plot(Sp1.t/ms, Sp1.i, '.b', markersize=5)\n",
    "# xlabel('Time (ms)')\n",
    "# ylabel('Neuron index')\n",
    "# text(30, n_i+8, 'e = {}    ||    i = {}    ||    n = {}    ||    p = {}'.format(exc, inh, n_i, p))\n",
    "# show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 80,
   "metadata": {},
   "outputs": [],
   "source": [
    "# hist(I2.w / gmax, 20)\n",
    "# xlabel('Outer Inh Weights / gmax');"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 81,
   "metadata": {},
   "outputs": [],
   "source": [
    "# hist(I1.w / gmax, 20)\n",
    "# xlabel('Inner Inh Weights / gmax');"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 82,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(5400, 8000)"
      ]
     },
     "execution_count": 82,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "shape(WmonI1.w) #synapses * time"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 83,
   "metadata": {},
   "outputs": [],
   "source": [
    "# hist(WmonI1.w[:,3000], 20)\n",
    "# xlabel('Weight / gmax')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 84,
   "metadata": {},
   "outputs": [],
   "source": [
    "from matplotlib.animation import FuncAnimation"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 90,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "application/javascript": [
       "/* Put everything inside the global mpl namespace */\n",
       "/* global mpl */\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(\n",
       "            '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",
       "\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 = document.createElement('div');\n",
       "    this.root.setAttribute('style', 'display: inline-block');\n",
       "    this._root_extra_style(this.root);\n",
       "\n",
       "    parent_element.appendChild(this.root);\n",
       "\n",
       "    this._init_header(this);\n",
       "    this._init_canvas(this);\n",
       "    this._init_toolbar(this);\n",
       "\n",
       "    var fig = this;\n",
       "\n",
       "    this.waiting = false;\n",
       "\n",
       "    this.ws.onopen = function () {\n",
       "        fig.send_message('supports_binary', { value: fig.supports_binary });\n",
       "        fig.send_message('send_image_mode', {});\n",
       "        if (fig.ratio !== 1) {\n",
       "            fig.send_message('set_device_pixel_ratio', {\n",
       "                device_pixel_ratio: fig.ratio,\n",
       "            });\n",
       "        }\n",
       "        fig.send_message('refresh', {});\n",
       "    };\n",
       "\n",
       "    this.imageObj.onload = function () {\n",
       "        if (fig.image_mode === 'full') {\n",
       "            // Full images could contain transparency (where diff images\n",
       "            // almost always do), so we need to clear the canvas so that\n",
       "            // there is no ghosting.\n",
       "            fig.context.clearRect(0, 0, fig.canvas.width, fig.canvas.height);\n",
       "        }\n",
       "        fig.context.drawImage(fig.imageObj, 0, 0);\n",
       "    };\n",
       "\n",
       "    this.imageObj.onunload = function () {\n",
       "        fig.ws.close();\n",
       "    };\n",
       "\n",
       "    this.ws.onmessage = this._make_on_message_function(this);\n",
       "\n",
       "    this.ondownload = ondownload;\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._init_header = function () {\n",
       "    var titlebar = document.createElement('div');\n",
       "    titlebar.classList =\n",
       "        'ui-dialog-titlebar ui-widget-header ui-corner-all ui-helper-clearfix';\n",
       "    var titletext = document.createElement('div');\n",
       "    titletext.classList = 'ui-dialog-title';\n",
       "    titletext.setAttribute(\n",
       "        'style',\n",
       "        'width: 100%; text-align: center; padding: 3px;'\n",
       "    );\n",
       "    titlebar.appendChild(titletext);\n",
       "    this.root.appendChild(titlebar);\n",
       "    this.header = titletext;\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._canvas_extra_style = function (_canvas_div) {};\n",
       "\n",
       "mpl.figure.prototype._root_extra_style = function (_canvas_div) {};\n",
       "\n",
       "mpl.figure.prototype._init_canvas = function () {\n",
       "    var fig = this;\n",
       "\n",
       "    var canvas_div = (this.canvas_div = document.createElement('div'));\n",
       "    canvas_div.setAttribute(\n",
       "        'style',\n",
       "        'border: 1px solid #ddd;' +\n",
       "            'box-sizing: content-box;' +\n",
       "            'clear: both;' +\n",
       "            'min-height: 1px;' +\n",
       "            'min-width: 1px;' +\n",
       "            'outline: 0;' +\n",
       "            'overflow: hidden;' +\n",
       "            'position: relative;' +\n",
       "            'resize: both;'\n",
       "    );\n",
       "\n",
       "    function on_keyboard_event_closure(name) {\n",
       "        return function (event) {\n",
       "            return fig.key_event(event, name);\n",
       "        };\n",
       "    }\n",
       "\n",
       "    canvas_div.addEventListener(\n",
       "        'keydown',\n",
       "        on_keyboard_event_closure('key_press')\n",
       "    );\n",
       "    canvas_div.addEventListener(\n",
       "        'keyup',\n",
       "        on_keyboard_event_closure('key_release')\n",
       "    );\n",
       "\n",
       "    this._canvas_extra_style(canvas_div);\n",
       "    this.root.appendChild(canvas_div);\n",
       "\n",
       "    var canvas = (this.canvas = document.createElement('canvas'));\n",
       "    canvas.classList.add('mpl-canvas');\n",
       "    canvas.setAttribute('style', 'box-sizing: content-box;');\n",
       "\n",
       "    this.context = canvas.getContext('2d');\n",
       "\n",
       "    var backingStore =\n",
       "        this.context.backingStorePixelRatio ||\n",
       "        this.context.webkitBackingStorePixelRatio ||\n",
       "        this.context.mozBackingStorePixelRatio ||\n",
       "        this.context.msBackingStorePixelRatio ||\n",
       "        this.context.oBackingStorePixelRatio ||\n",
       "        this.context.backingStorePixelRatio ||\n",
       "        1;\n",
       "\n",
       "    this.ratio = (window.devicePixelRatio || 1) / backingStore;\n",
       "\n",
       "    var rubberband_canvas = (this.rubberband_canvas = document.createElement(\n",
       "        'canvas'\n",
       "    ));\n",
       "    rubberband_canvas.setAttribute(\n",
       "        'style',\n",
       "        'box-sizing: content-box; position: absolute; left: 0; top: 0; z-index: 1;'\n",
       "    );\n",
       "\n",
       "    // Apply a ponyfill if ResizeObserver is not implemented by browser.\n",
       "    if (this.ResizeObserver === undefined) {\n",
       "        if (window.ResizeObserver !== undefined) {\n",
       "            this.ResizeObserver = window.ResizeObserver;\n",
       "        } else {\n",
       "            var obs = _JSXTOOLS_RESIZE_OBSERVER({});\n",
       "            this.ResizeObserver = obs.ResizeObserver;\n",
       "        }\n",
       "    }\n",
       "\n",
       "    this.resizeObserverInstance = new this.ResizeObserver(function (entries) {\n",
       "        var nentries = entries.length;\n",
       "        for (var i = 0; i < nentries; i++) {\n",
       "            var entry = entries[i];\n",
       "            var width, height;\n",
       "            if (entry.contentBoxSize) {\n",
       "                if (entry.contentBoxSize instanceof Array) {\n",
       "                    // Chrome 84 implements new version of spec.\n",
       "                    width = entry.contentBoxSize[0].inlineSize;\n",
       "                    height = entry.contentBoxSize[0].blockSize;\n",
       "                } else {\n",
       "                    // Firefox implements old version of spec.\n",
       "                    width = entry.contentBoxSize.inlineSize;\n",
       "                    height = entry.contentBoxSize.blockSize;\n",
       "                }\n",
       "            } else {\n",
       "                // Chrome <84 implements even older version of spec.\n",
       "                width = entry.contentRect.width;\n",
       "                height = entry.contentRect.height;\n",
       "            }\n",
       "\n",
       "            // Keep the size of the canvas and rubber band canvas in sync with\n",
       "            // the canvas container.\n",
       "            if (entry.devicePixelContentBoxSize) {\n",
       "                // Chrome 84 implements new version of spec.\n",
       "                canvas.setAttribute(\n",
       "                    'width',\n",
       "                    entry.devicePixelContentBoxSize[0].inlineSize\n",
       "                );\n",
       "                canvas.setAttribute(\n",
       "                    'height',\n",
       "                    entry.devicePixelContentBoxSize[0].blockSize\n",
       "                );\n",
       "            } else {\n",
       "                canvas.setAttribute('width', width * fig.ratio);\n",
       "                canvas.setAttribute('height', height * fig.ratio);\n",
       "            }\n",
       "            canvas.setAttribute(\n",
       "                'style',\n",
       "                'width: ' + width + 'px; height: ' + height + 'px;'\n",
       "            );\n",
       "\n",
       "            rubberband_canvas.setAttribute('width', width);\n",
       "            rubberband_canvas.setAttribute('height', height);\n",
       "\n",
       "            // And update the size in Python. We ignore the initial 0/0 size\n",
       "            // that occurs as the element is placed into the DOM, which should\n",
       "            // otherwise not happen due to the minimum size styling.\n",
       "            if (fig.ws.readyState == 1 && width != 0 && height != 0) {\n",
       "                fig.request_resize(width, height);\n",
       "            }\n",
       "        }\n",
       "    });\n",
       "    this.resizeObserverInstance.observe(canvas_div);\n",
       "\n",
       "    function on_mouse_event_closure(name) {\n",
       "        return function (event) {\n",
       "            return fig.mouse_event(event, name);\n",
       "        };\n",
       "    }\n",
       "\n",
       "    rubberband_canvas.addEventListener(\n",
       "        'mousedown',\n",
       "        on_mouse_event_closure('button_press')\n",
       "    );\n",
       "    rubberband_canvas.addEventListener(\n",
       "        'mouseup',\n",
       "        on_mouse_event_closure('button_release')\n",
       "    );\n",
       "    rubberband_canvas.addEventListener(\n",
       "        'dblclick',\n",
       "        on_mouse_event_closure('dblclick')\n",
       "    );\n",
       "    // Throttle sequential mouse events to 1 every 20ms.\n",
       "    rubberband_canvas.addEventListener(\n",
       "        'mousemove',\n",
       "        on_mouse_event_closure('motion_notify')\n",
       "    );\n",
       "\n",
       "    rubberband_canvas.addEventListener(\n",
       "        'mouseenter',\n",
       "        on_mouse_event_closure('figure_enter')\n",
       "    );\n",
       "    rubberband_canvas.addEventListener(\n",
       "        'mouseleave',\n",
       "        on_mouse_event_closure('figure_leave')\n",
       "    );\n",
       "\n",
       "    canvas_div.addEventListener('wheel', function (event) {\n",
       "        if (event.deltaY < 0) {\n",
       "            event.step = 1;\n",
       "        } else {\n",
       "            event.step = -1;\n",
       "        }\n",
       "        on_mouse_event_closure('scroll')(event);\n",
       "    });\n",
       "\n",
       "    canvas_div.appendChild(canvas);\n",
       "    canvas_div.appendChild(rubberband_canvas);\n",
       "\n",
       "    this.rubberband_context = rubberband_canvas.getContext('2d');\n",
       "    this.rubberband_context.strokeStyle = '#000000';\n",
       "\n",
       "    this._resize_canvas = function (width, height, forward) {\n",
       "        if (forward) {\n",
       "            canvas_div.style.width = width + 'px';\n",
       "            canvas_div.style.height = height + 'px';\n",
       "        }\n",
       "    };\n",
       "\n",
       "    // Disable right mouse context menu.\n",
       "    this.rubberband_canvas.addEventListener('contextmenu', function (_e) {\n",
       "        event.preventDefault();\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 toolbar = document.createElement('div');\n",
       "    toolbar.classList = 'mpl-toolbar';\n",
       "    this.root.appendChild(toolbar);\n",
       "\n",
       "    function on_click_closure(name) {\n",
       "        return function (_event) {\n",
       "            return fig.toolbar_button_onclick(name);\n",
       "        };\n",
       "    }\n",
       "\n",
       "    function on_mouseover_closure(tooltip) {\n",
       "        return function (event) {\n",
       "            if (!event.currentTarget.disabled) {\n",
       "                return fig.toolbar_button_onmouseover(tooltip);\n",
       "            }\n",
       "        };\n",
       "    }\n",
       "\n",
       "    fig.buttons = {};\n",
       "    var buttonGroup = document.createElement('div');\n",
       "    buttonGroup.classList = 'mpl-button-group';\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",
       "            /* Instead of a spacer, we start a new button group. */\n",
       "            if (buttonGroup.hasChildNodes()) {\n",
       "                toolbar.appendChild(buttonGroup);\n",
       "            }\n",
       "            buttonGroup = document.createElement('div');\n",
       "            buttonGroup.classList = 'mpl-button-group';\n",
       "            continue;\n",
       "        }\n",
       "\n",
       "        var button = (fig.buttons[name] = document.createElement('button'));\n",
       "        button.classList = 'mpl-widget';\n",
       "        button.setAttribute('role', 'button');\n",
       "        button.setAttribute('aria-disabled', 'false');\n",
       "        button.addEventListener('click', on_click_closure(method_name));\n",
       "        button.addEventListener('mouseover', on_mouseover_closure(tooltip));\n",
       "\n",
       "        var icon_img = document.createElement('img');\n",
       "        icon_img.src = '_images/' + image + '.png';\n",
       "        icon_img.srcset = '_images/' + image + '_large.png 2x';\n",
       "        icon_img.alt = tooltip;\n",
       "        button.appendChild(icon_img);\n",
       "\n",
       "        buttonGroup.appendChild(button);\n",
       "    }\n",
       "\n",
       "    if (buttonGroup.hasChildNodes()) {\n",
       "        toolbar.appendChild(buttonGroup);\n",
       "    }\n",
       "\n",
       "    var fmt_picker = document.createElement('select');\n",
       "    fmt_picker.classList = 'mpl-widget';\n",
       "    toolbar.appendChild(fmt_picker);\n",
       "    this.format_dropdown = fmt_picker;\n",
       "\n",
       "    for (var ind in mpl.extensions) {\n",
       "        var fmt = mpl.extensions[ind];\n",
       "        var option = document.createElement('option');\n",
       "        option.selected = fmt === mpl.default_extension;\n",
       "        option.innerHTML = fmt;\n",
       "        fmt_picker.appendChild(option);\n",
       "    }\n",
       "\n",
       "    var status_bar = document.createElement('span');\n",
       "    status_bar.classList = 'mpl-message';\n",
       "    toolbar.appendChild(status_bar);\n",
       "    this.message = status_bar;\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",
       "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",
       "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], msg['forward']);\n",
       "        fig.send_message('refresh', {});\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_rubberband = function (fig, msg) {\n",
       "    var x0 = msg['x0'] / fig.ratio;\n",
       "    var y0 = (fig.canvas.height - msg['y0']) / fig.ratio;\n",
       "    var x1 = msg['x1'] / fig.ratio;\n",
       "    var y1 = (fig.canvas.height - msg['y1']) / fig.ratio;\n",
       "    x0 = Math.floor(x0) + 0.5;\n",
       "    y0 = Math.floor(y0) + 0.5;\n",
       "    x1 = Math.floor(x1) + 0.5;\n",
       "    y1 = Math.floor(y1) + 0.5;\n",
       "    var min_x = Math.min(x0, x1);\n",
       "    var min_y = Math.min(y0, y1);\n",
       "    var width = Math.abs(x1 - x0);\n",
       "    var height = Math.abs(y1 - y0);\n",
       "\n",
       "    fig.rubberband_context.clearRect(\n",
       "        0,\n",
       "        0,\n",
       "        fig.canvas.width / fig.ratio,\n",
       "        fig.canvas.height / fig.ratio\n",
       "    );\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",
       "    fig.rubberband_canvas.style.cursor = msg['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.handle_history_buttons = function (fig, msg) {\n",
       "    for (var key in msg) {\n",
       "        if (!(key in fig.buttons)) {\n",
       "            continue;\n",
       "        }\n",
       "        fig.buttons[key].disabled = !msg[key];\n",
       "        fig.buttons[key].setAttribute('aria-disabled', !msg[key]);\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_navigate_mode = function (fig, msg) {\n",
       "    if (msg['mode'] === 'PAN') {\n",
       "        fig.buttons['Pan'].classList.add('active');\n",
       "        fig.buttons['Zoom'].classList.remove('active');\n",
       "    } else if (msg['mode'] === 'ZOOM') {\n",
       "        fig.buttons['Pan'].classList.remove('active');\n",
       "        fig.buttons['Zoom'].classList.add('active');\n",
       "    } else {\n",
       "        fig.buttons['Pan'].classList.remove('active');\n",
       "        fig.buttons['Zoom'].classList.remove('active');\n",
       "    }\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",
       "            var img = evt.data;\n",
       "            if (img.type !== 'image/png') {\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",
       "                img.type = 'image/png';\n",
       "            }\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",
       "\n",
       "            fig.imageObj.src = (window.URL || window.webkitURL).createObjectURL(\n",
       "                img\n",
       "            );\n",
       "            fig.updated_canvas_event();\n",
       "            fig.waiting = false;\n",
       "            return;\n",
       "        } else if (\n",
       "            typeof evt.data === 'string' &&\n",
       "            evt.data.slice(0, 21) === 'data:image/png;base64'\n",
       "        ) {\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(\n",
       "                \"No handler for the '\" + msg_type + \"' message type: \",\n",
       "                msg\n",
       "            );\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(\n",
       "                    \"Exception inside the 'handler_\" + msg_type + \"' callback:\",\n",
       "                    e,\n",
       "                    e.stack,\n",
       "                    msg\n",
       "                );\n",
       "            }\n",
       "        }\n",
       "    };\n",
       "};\n",
       "\n",
       "// from https://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",
       "    }\n",
       "    if (e.target) {\n",
       "        targ = e.target;\n",
       "    } else if (e.srcElement) {\n",
       "        targ = e.srcElement;\n",
       "    }\n",
       "    if (targ.nodeType === 3) {\n",
       "        // defeat Safari bug\n",
       "        targ = targ.parentNode;\n",
       "    }\n",
       "\n",
       "    // pageX,Y are the mouse positions relative to the document\n",
       "    var boundingRect = targ.getBoundingClientRect();\n",
       "    var x = e.pageX - (boundingRect.left + document.body.scrollLeft);\n",
       "    var y = e.pageY - (boundingRect.top + document.body.scrollTop);\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",
       " * https://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",
       "        }\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",
       "        this.canvas.focus();\n",
       "        this.canvas_div.focus();\n",
       "    }\n",
       "\n",
       "    var x = canvas_pos.x * this.ratio;\n",
       "    var y = canvas_pos.y * this.ratio;\n",
       "\n",
       "    this.send_message(name, {\n",
       "        x: x,\n",
       "        y: y,\n",
       "        button: event.button,\n",
       "        step: event.step,\n",
       "        guiEvent: simpleKeys(event),\n",
       "    });\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",
       "    // Prevent repeat events\n",
       "    if (name === 'key_press') {\n",
       "        if (event.key === this._key) {\n",
       "            return;\n",
       "        } else {\n",
       "            this._key = event.key;\n",
       "        }\n",
       "    }\n",
       "    if (name === 'key_release') {\n",
       "        this._key = null;\n",
       "    }\n",
       "\n",
       "    var value = '';\n",
       "    if (event.ctrlKey && event.key !== 'Control') {\n",
       "        value += 'ctrl+';\n",
       "    }\n",
       "    else if (event.altKey && event.key !== 'Alt') {\n",
       "        value += 'alt+';\n",
       "    }\n",
       "    else if (event.shiftKey && event.key !== 'Shift') {\n",
       "        value += 'shift+';\n",
       "    }\n",
       "\n",
       "    value += 'k' + event.key;\n",
       "\n",
       "    this._key_event_extra(event, name);\n",
       "\n",
       "    this.send_message(name, { key: value, 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",
       "\n",
       "///////////////// REMAINING CONTENT GENERATED BY embed_js.py /////////////////\n",
       "// prettier-ignore\n",
       "var _JSXTOOLS_RESIZE_OBSERVER=function(A){var t,i=new WeakMap,n=new WeakMap,a=new WeakMap,r=new WeakMap,o=new Set;function s(e){if(!(this instanceof s))throw new TypeError(\"Constructor requires 'new' operator\");i.set(this,e)}function h(){throw new TypeError(\"Function is not a constructor\")}function c(e,t,i,n){e=0 in arguments?Number(arguments[0]):0,t=1 in arguments?Number(arguments[1]):0,i=2 in arguments?Number(arguments[2]):0,n=3 in arguments?Number(arguments[3]):0,this.right=(this.x=this.left=e)+(this.width=i),this.bottom=(this.y=this.top=t)+(this.height=n),Object.freeze(this)}function d(){t=requestAnimationFrame(d);var s=new WeakMap,p=new Set;o.forEach((function(t){r.get(t).forEach((function(i){var r=t instanceof window.SVGElement,o=a.get(t),d=r?0:parseFloat(o.paddingTop),f=r?0:parseFloat(o.paddingRight),l=r?0:parseFloat(o.paddingBottom),u=r?0:parseFloat(o.paddingLeft),g=r?0:parseFloat(o.borderTopWidth),m=r?0:parseFloat(o.borderRightWidth),w=r?0:parseFloat(o.borderBottomWidth),b=u+f,F=d+l,v=(r?0:parseFloat(o.borderLeftWidth))+m,W=g+w,y=r?0:t.offsetHeight-W-t.clientHeight,E=r?0:t.offsetWidth-v-t.clientWidth,R=b+v,z=F+W,M=r?t.width:parseFloat(o.width)-R-E,O=r?t.height:parseFloat(o.height)-z-y;if(n.has(t)){var k=n.get(t);if(k[0]===M&&k[1]===O)return}n.set(t,[M,O]);var S=Object.create(h.prototype);S.target=t,S.contentRect=new c(u,d,M,O),s.has(i)||(s.set(i,[]),p.add(i)),s.get(i).push(S)}))})),p.forEach((function(e){i.get(e).call(e,s.get(e),e)}))}return s.prototype.observe=function(i){if(i instanceof window.Element){r.has(i)||(r.set(i,new Set),o.add(i),a.set(i,window.getComputedStyle(i)));var n=r.get(i);n.has(this)||n.add(this),cancelAnimationFrame(t),t=requestAnimationFrame(d)}},s.prototype.unobserve=function(i){if(i instanceof window.Element&&r.has(i)){var n=r.get(i);n.has(this)&&(n.delete(this),n.size||(r.delete(i),o.delete(i))),n.size||r.delete(i),o.size||cancelAnimationFrame(t)}},A.DOMRectReadOnly=c,A.ResizeObserver=s,A.ResizeObserverEntry=h,A}; // eslint-disable-line\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\", \"Left button pans, Right button zooms\\nx/y fixes axis, CTRL fixes aspect\", \"fa fa-arrows icon-move\", \"pan\"], [\"Zoom\", \"Zoom to rectangle\\nx/y fixes axis\", \"fa fa-square-o icon-check-empty\", \"zoom\"], [\"\", \"\", \"\", \"\"], [\"Download\", \"Download plot\", \"fa fa-floppy-o icon-save\", \"download\"]];\n",
       "\n",
       "mpl.extensions = [\"eps\", \"jpeg\", \"pgf\", \"pdf\", \"png\", \"ps\", \"raw\", \"svg\", \"tif\"];\n",
       "\n",
       "mpl.default_extension = \"png\";/* global mpl */\n",
       "\n",
       "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.binaryType = comm.kernel.ws.binaryType;\n",
       "    ws.readyState = comm.kernel.ws.readyState;\n",
       "    function updateReadyState(_event) {\n",
       "        if (comm.kernel.ws) {\n",
       "            ws.readyState = comm.kernel.ws.readyState;\n",
       "        } else {\n",
       "            ws.readyState = 3; // Closed state.\n",
       "        }\n",
       "    }\n",
       "    comm.kernel.ws.addEventListener('open', updateReadyState);\n",
       "    comm.kernel.ws.addEventListener('close', updateReadyState);\n",
       "    comm.kernel.ws.addEventListener('error', updateReadyState);\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",
       "        var data = msg['content']['data'];\n",
       "        if (data['blob'] !== undefined) {\n",
       "            data = {\n",
       "                data: new Blob(msg['buffers'], { type: data['blob'] }),\n",
       "            };\n",
       "        }\n",
       "        // Pass the mpl event to the overridden (by mpl) onmessage function.\n",
       "        ws.onmessage(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 = document.getElementById(id);\n",
       "    var ws_proxy = comm_websocket_adapter(comm);\n",
       "\n",
       "    function ondownload(figure, _format) {\n",
       "        window.open(figure.canvas.toDataURL());\n",
       "    }\n",
       "\n",
       "    var fig = new mpl.figure(id, ws_proxy, ondownload, element);\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;\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",
       "    fig.cell_info[0].output_area.element.on(\n",
       "        'cleared',\n",
       "        { fig: fig },\n",
       "        fig._remove_fig_handler\n",
       "    );\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_close = function (fig, msg) {\n",
       "    var width = fig.canvas.width / fig.ratio;\n",
       "    fig.cell_info[0].output_area.element.off(\n",
       "        'cleared',\n",
       "        fig._remove_fig_handler\n",
       "    );\n",
       "    fig.resizeObserverInstance.unobserve(fig.canvas_div);\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.innerHTML =\n",
       "        '<img src=\"' + dataURL + '\" width=\"' + width + '\">';\n",
       "    fig.close_ws(fig, msg);\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.close_ws = function (fig, msg) {\n",
       "    fig.send_message('closing', msg);\n",
       "    // fig.ws.close()\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.push_to_output = function (_remove_interactive) {\n",
       "    // Turn the data on the canvas into data in the output cell.\n",
       "    var width = this.canvas.width / this.ratio;\n",
       "    var dataURL = this.canvas.toDataURL();\n",
       "    this.cell_info[1]['text/html'] =\n",
       "        '<img src=\"' + dataURL + '\" width=\"' + width + '\">';\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.updated_canvas_event = function () {\n",
       "    // Tell IPython that the notebook contents must change.\n",
       "    IPython.notebook.set_dirty(true);\n",
       "    this.send_message('ack', {});\n",
       "    var fig = this;\n",
       "    // Wait a second, then push the new image to the DOM so\n",
       "    // that it is saved nicely (might be nice to debounce this).\n",
       "    setTimeout(function () {\n",
       "        fig.push_to_output();\n",
       "    }, 1000);\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._init_toolbar = function () {\n",
       "    var fig = this;\n",
       "\n",
       "    var toolbar = document.createElement('div');\n",
       "    toolbar.classList = 'btn-toolbar';\n",
       "    this.root.appendChild(toolbar);\n",
       "\n",
       "    function on_click_closure(name) {\n",
       "        return function (_event) {\n",
       "            return fig.toolbar_button_onclick(name);\n",
       "        };\n",
       "    }\n",
       "\n",
       "    function on_mouseover_closure(tooltip) {\n",
       "        return function (event) {\n",
       "            if (!event.currentTarget.disabled) {\n",
       "                return fig.toolbar_button_onmouseover(tooltip);\n",
       "            }\n",
       "        };\n",
       "    }\n",
       "\n",
       "    fig.buttons = {};\n",
       "    var buttonGroup = document.createElement('div');\n",
       "    buttonGroup.classList = 'btn-group';\n",
       "    var button;\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",
       "            /* Instead of a spacer, we start a new button group. */\n",
       "            if (buttonGroup.hasChildNodes()) {\n",
       "                toolbar.appendChild(buttonGroup);\n",
       "            }\n",
       "            buttonGroup = document.createElement('div');\n",
       "            buttonGroup.classList = 'btn-group';\n",
       "            continue;\n",
       "        }\n",
       "\n",
       "        button = fig.buttons[name] = document.createElement('button');\n",
       "        button.classList = 'btn btn-default';\n",
       "        button.href = '#';\n",
       "        button.title = name;\n",
       "        button.innerHTML = '<i class=\"fa ' + image + ' fa-lg\"></i>';\n",
       "        button.addEventListener('click', on_click_closure(method_name));\n",
       "        button.addEventListener('mouseover', on_mouseover_closure(tooltip));\n",
       "        buttonGroup.appendChild(button);\n",
       "    }\n",
       "\n",
       "    if (buttonGroup.hasChildNodes()) {\n",
       "        toolbar.appendChild(buttonGroup);\n",
       "    }\n",
       "\n",
       "    // Add the status bar.\n",
       "    var status_bar = document.createElement('span');\n",
       "    status_bar.classList = 'mpl-message pull-right';\n",
       "    toolbar.appendChild(status_bar);\n",
       "    this.message = status_bar;\n",
       "\n",
       "    // Add the close button to the window.\n",
       "    var buttongrp = document.createElement('div');\n",
       "    buttongrp.classList = 'btn-group inline pull-right';\n",
       "    button = document.createElement('button');\n",
       "    button.classList = 'btn btn-mini btn-primary';\n",
       "    button.href = '#';\n",
       "    button.title = 'Stop Interaction';\n",
       "    button.innerHTML = '<i class=\"fa fa-power-off icon-remove icon-large\"></i>';\n",
       "    button.addEventListener('click', function (_evt) {\n",
       "        fig.handle_close(fig, {});\n",
       "    });\n",
       "    button.addEventListener(\n",
       "        'mouseover',\n",
       "        on_mouseover_closure('Stop Interaction')\n",
       "    );\n",
       "    buttongrp.appendChild(button);\n",
       "    var titlebar = this.root.querySelector('.ui-dialog-titlebar');\n",
       "    titlebar.insertBefore(buttongrp, titlebar.firstChild);\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._remove_fig_handler = function (event) {\n",
       "    var fig = event.data.fig;\n",
       "    if (event.target !== this) {\n",
       "        // Ignore bubbled events from children.\n",
       "        return;\n",
       "    }\n",
       "    fig.close_ws(fig, {});\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._root_extra_style = function (el) {\n",
       "    el.style.boxSizing = 'content-box'; // override notebook setting of border-box.\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._canvas_extra_style = function (el) {\n",
       "    // this is important to make the div 'focusable\n",
       "    el.setAttribute('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",
       "    } else {\n",
       "        // location in version 2\n",
       "        IPython.keyboard_manager.register_events(el);\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure.prototype._key_event_extra = function (event, _name) {\n",
       "    // Check for shift+enter\n",
       "    if (event.shiftKey && event.which === 13) {\n",
       "        this.canvas_div.blur();\n",
       "        // select the cell after this one\n",
       "        var index = IPython.notebook.find_cell_index(this.cell_info[0]);\n",
       "        IPython.notebook.select(index + 1);\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_save = function (fig, _msg) {\n",
       "    fig.ondownload(fig, null);\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(\n",
       "        'matplotlib',\n",
       "        mpl.mpl_figure_comm\n",
       "    );\n",
       "}\n"
      ],
      "text/plain": [
       "<IPython.core.display.Javascript object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/html": [
       "<img src=\"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAyAAAAJYCAYAAACadoJwAAAAAXNSR0IArs4c6QAAIABJREFUeF7t3Qu0lFd1OPBjQdIYta1miVbQAPEtIGhjAvHRaqS2aqzGN6sGvQohUlIVVGjXH1wCCqkWkQDxWrSNtb7SamKL2kStgagJFYINPhCiYECiVduYtAHMf31fFpTXZWbYZ8bzXX6zVlbbO+fu2fPbO+nZ93zfzP3uvffee5MHAQIECBAgQIAAAQIEeiBwPwNID5S9BAECBAgQIECAAAECtYABRCMQIECAAAECBAgQINAzAQNIz6i9EAECBAgQIECAAAECBhA9QIAAAQIECBAgQIBAzwQMID2j9kIECBAgQIAAAQIECBhA9AABAgQIECBAgAABAj0TMID0jNoLESBAgAABAgQIECBgANEDBAgQIECAAAECBAj0TMAA0jNqL0SAAAECBAgQIECAgAFEDxAgQIAAAQIECBAg0DMBA0jPqL0QAQIECBAgQIAAAQIGED1AgAABAgQIECBAgEDPBAwgPaP2QgQIECBAgAABAgQIGED0AAECBAgQIECAAAECPRMwgPSM2gsRIECAAAECBAgQIGAA0QMECBAgQIAAAQIECPRMwADSM2ovRIAAAQIECBAgQICAAUQPECBAgAABAgQIECDQMwEDSM+ovRABAgQIECBAgAABAgYQPUCAAAECBAgQIECAQM8EDCA9o/ZCBAgQIECAAAECBAgYQPQAAQIECBAgQIAAAQI9EzCA9IzaCxEgQIAAAQIECBAgYADRAwQIECBAgAABAgQI9EzAANIzai9EgAABAgQIECBAgIABRA8QIECAAAECBAgQINAzAQNIz6i9EAECBAgQIECAAAECBhA9QIAAAQIECBAgQIBAzwQMID2j9kIECBAgQIAAAQIECBhA9AABAgQIECBAgAABAj0TMID0jNoLESBAgAABAgQIECBgANEDBAgQIECAAAECBAj0TMAA0jNqL0SAAAECBAgQIECAgAFEDxAgQIAAAQIECBAg0DMBA0jPqL0QAQIECBAgQIAAAQIGED1AgAABAgQIECBAgEDPBAwgPaP2QgQIECBAgAABAgQIGED0AAECBAgQIECAAAECPRMwgPSM2gsRIECAAAECBAgQIGAA0QMECBAgQIAAAQIECPRMwADSM2ovRIAAAQIECBAgQICAAUQPECBAgAABAgQIECDQMwEDSM+ovRABAgQIECBAgAABAgYQPUCAAAECBAgQIECAQM8EDCA9o/ZCBAgQIECAAAECBAgYQPQAAQIECBAgQIAAAQI9EzCA9IzaCxEgQIAAAQIECBAgYADRAwQIECBAgAABAgQI9EzAANIzai9EgAABAgQIECBAgIABRA8QIECAAAECBAgQINAzAQNIz6i9EAECBAgQIECAAAECBhA9QIAAAQIECBAgQIBAzwQMID2j9kIECBAgQIAAAQIECBhA9AABAgQIECBAgAABAj0TMID0jNoLESBAgAABAgQIECBgANEDBAgQIECAAAECBAj0TMAA0jNqL0SAAAECBAgQIECAgAFEDxAgQIAAAQIECBAg0DMBA0jPqL0QAQIECBAgQIAAAQIGED1AgAABAgQIECBAgEDPBAwgPaP2QgQIECBAgAABAgQIGED0AAECBAgQIECAAAECPRMwgPSM2gsRIECAAAECBAgQIGAA0QMECBAgQIAAAQIECPRMwADSM2ovRIAAAQIECBAgQICAAUQPECBAgAABAgQIECDQMwEDSM+ovRABAgQIECBAgAABAgYQPUCAAAECBAgQIECAQM8EDCA9o/ZCBAgQIECAAAECBAgYQPQAAQIECBAgQIAAAQI9EzCA9Iy69Qv95Cc/SZ///OfTWWedlU4//fTWv2AFAQIECBAgQIBATwXuvvvudNttt6UpU6akM888s6evPVhezABSUCU/+tGPpqlTpxaUkVQIECBAgAABAgSOJ3DVVVel17zmNXBOQsAAchJo3fqV9evXp/PPPz9VDf2EJzyhWy8jLgECBAgQIECAwEkKbN26tf6D8Q033JAmT558klFO7V8zgBRU/3//939PT33qU9PGjRvTxIkTC8pMKgQIECBAgAABApWA/Vq8DwwgccNsETR0NkqBCBAgQIAAAQJdEbBfi7MaQOKG2SJo6GyUAhEgQIAAAQIEuiJgvxZnNYDEDbNF0NDZKAUiQIAAAQIECHRFwH4tzmoAiRtmi6Chs1EKRIAAAQIECBDoioD9WpzVABI3zBZBQ2ejFIgAAQIECBAg0BUB+7U4qwEkbpgtgobORikQAQIECBAgQKArAvZrcVYDSNwwWwQNnY1SIAIECBAgQIBAVwTs1+KsBpC4YbYIGjobpUAECBAgQIAAga4I2K/FWQ0gccNsETR0NkqBCBAgQIAAAQJdEbBfi7MaQOKG2SJo6GyUAhEgQIAAAQIEuiJgvxZnNYDEDbNF0NDZKAUiQIAAAQIECHRFwH4tzmoAiRtmi6Chs1EKRIAAAQIECBDoioD9WpzVABI3zBZBQ2ejFIgAAQIECBAg0BUB+7U4a2MHkJ07d6Z3vvOd6brrrku7d+9Ow4cPT8961rPS/Pnz02Mf+9hDMgcOHEjLli1L/f39qfqdkSNHpr6+vjRnzpw0ZMiQIwS7tbbdMmnodqWsI0CAAAECBAj8egTs1+LujRxAfvrTn6YnP/nJ6Z577kmXXHJJGjVqVNq2bVtatWpVut/97pe2bNmSRowYUevMnDmz/vm0adPSpEmT0oYNG9LatWvrn69cufIIwW6tbbdMGrpdKesIECBAgAABAr8eAfu1uHsjB5ArrrgiXXrppemzn/1seuELX3hI4dOf/nS66KKL0vve97502WWX1YPI+PHj06xZs9Ly5csPrZs9e3ZasWJF2rx5cxo7dmz9826t7aREGroTLWsJECBAgAABAr0XsF+LmzdyAHn3u9+d3vGOd6SbbropPe1pTzukcOONN9anHGvWrElvfOMb68uxFi9enLZv316fkhx87NixI40ePTrNmzcvLVq0qP5xt9Z2UiIN3YmWtQQIECBAgACB3gvYr8XNGzmAVIPHOeeck84999x0+eWXp7POOqu+BOvNb35zuvPOO+vB5MEPfnCaMmVKfcqxZ8+eY6Sqe0YmTJiQ1q1bVz/XrbWdlEhDd6JlLQECBAgQIECg9wL2a3HzRg4g1du+8sor6xOM6n6Qg49nPvOZ6eqrr04PfehD6x9Vl1cNGzYsbdy48RipiRMnpn379tWXXnVz7UAlqm6cr/45/LF169Y0derUOt8qPw8CBAgQIECAAIGyBAwg8Xo0dgD5x3/8x7R69er65OLss8+uB4nq066qm9M///nPpzPOOCONGTOm/nSs6sbzox/VpVp79+6tT06qR7fWDlSiBQsWpIULFx73aQNIvLFFIECgHIGz3v65cpLpQia3vfuPuxBVSAIEShUwgMQr08gBpDrleMUrXpGqBjh4E3lF8YUvfKEeSN7znvekuXPnOgGJ94cIBAgQCAsYQMKEAhAgUJCAASRejEYOINX3fdxxxx3p1ltvPUaguvejuhTr2muv7dp9HZ3cL9JJiTR0J1rWEiDQFAEDSFMqJU8CBNoRsF9rR+nEaxo5gDzucY+r39V3vvOdI97dvffemx70oAel888/v765vLpHZMmSJW19Cla31nZSIg3diZa1BAg0RcAA0pRKyZMAgXYE7NfaURqEA8iFF16YrrnmmrR+/fp03nnnHXqHn/rUp9LLXvay+iN13/Wud9WfgFV90tVA3wOyadOmNG7cuPr3u7W2kxJp6E60rCVAoCkCBpCmVEqeBAi0I2C/1o7SIBxAqpvKn/3sZ6fTTjut/kbz6gby6ib06pOxHvKQh6RvfvOb6eEPf3j9zmfMmFF/L0j1TeiTJ0+uh5bqm9CnT59e38R++KNba9stk4ZuV8o6AgSaJGAAaVK15EqAQCsB+7VWQq2fb+QlWNXbuuWWW9I73/nOdPPNN6fbb7+9HjwuuOCC+uTj0Y9+9KF3vn///rR06dLU39+fdu3alUaMGJH6+vrqm9SHDh16hFC31rYuw30rNHS7UtYRINAkAQNIk6olVwIEWgnYr7USav18YweQ1m+teSs0dPNqJmMCBFoLGEBaG1lBgEBzBOzX4rUygMQNs0XQ0NkoBSJAoCABA0hBxZAKAQJhAfu1MGEygMQNs0XQ0NkoBSJAoCABA0hBxZAKAQJhAfu1MKEBJE6YL4KGzmcpEgEC5QgYQMqphUwIEIgL2K/FDZ2AxA2zRdDQ2SgFIkCgIAEDSEHFkAoBAmEB+7UwoROQOGG+CBo6n6VIBAiUI2AAKacWMiFAIC5gvxY3dAISN8wWQUNnoxSIAIGCBAwgBRVDKgQIhAXs18KETkDihPkiaOh8liIRIFCOgAGknFrIhACBuID9WtzQCUjcMFsEDZ2NUiACBAoSMIAUVAypECAQFrBfCxM6AYkT5ougofNZikSAQDkCBpByaiETAgTiAvZrcUMnIHHDbBE0dDZKgQgQKEjAAFJQMaRCgEBYwH4tTOgEJE6YL4KGzmcpEgEC5QgYQMqphUwIEIgL2K/FDZ2AxA2zRdDQ2SgFIkCgIAEDSEHFkAoBAmEB+7UwoROQOGG+CBo6n6VIBAiUI2AAKacWMiFAIC5gvxY3dAISN8wWQUNnoxSIAIGCBAwgBRVDKgQIhAXs18KETkDihPkiaOh8liIRIFCOgAGknFrIhACBuID9WtzQCUjcMFsEDZ2NUiACBAoSMIAUVAypECAQFrBfCxM6AYkT5ougofNZikSAQDkCBpByaiETAgTiAvZrcUMnIHHDbBE0dDZKgQgQKEjAAFJQMaRCgEBYwH4tTOgEJE6YL4KGzmcpEgEC5QgYQMqphUwIEIgL2K/FDZ2AxA2zRdDQ2SgFIkCgIAEDSEHFkAoBAmEB+7UwoROQOGG+CBo6n6VIBAiUI2AAKacWMiFAIC5gvxY3dAISN8wWQUNnoxSIAIGCBAwgBRVDKgQIhAXs18KETkDihPkiaOh8liIRIFCOgAGknFrIhACBuID9WtzQCUjcMFsEDZ2NUiACBAoSMIAUVAypECAQFrBfCxM6AYkT5ougofNZikSAQDkCBpByaiETAgTiAvZrcUMnIHHDbBE0dDZKgQgQKEjAAFJQMaRCgEBYwH4tTOgEJE6YL4KGzmcpEgEC5QgYQMqphUwIEIgL2K/FDZ2AxA2zRdDQ2SgFIkCgIAEDSEHFkAoBAmEB+7UwoROQOGG+CBo6n6VIBAiUI2AAKacWMiFAIC5gvxY3dAISN8wWQUNnoxSIAIGCBAwgBRVDKgQIhAXs18KETkDihPkiaOh8liIRIFCOgAGknFrIhACBuID9WtzQCUjcMFsEDZ2NUiACBAoSMIAUVAypECAQFrBfCxM6AYkT5ougofNZikSAQDkCBpByaiETAgTiAvZrcUMnIHHDbBE0dDZKgQgQKEjAAFJQMaRCgEBYwH4tTOgEJE6YL4KGzmcpEgEC5QgYQMqphUwIEIgL2K/FDZ2AxA2zRdDQ2SgFIkCgIAEDSEHFkAoBAmEB+7UwoROQOGG+CBo6n6VIBAiUI2AAKacWMiFAIC5gvxY3dAISN8wWQUNnoxSIAIGCBAwgBRVDKgQIhAXs18KETkDihPkiaOh8liIRIFCOgAGknFrIhACBuID9WtzQCUjcMFsEDZ2NUiACBAoSMIAUVAypECAQFrBfCxM6AYkT5ougofNZikSAQDkCBpByaiETAgTiAvZrccNGnoBcfPHF6SMf+ciA7/5d73pXmj9/fv38gQMH0rJly1J/f3/auXNnGjlyZOrr60tz5sxJQ4YMOSJGt9a2WyYN3a6UdQQINEnAANKkasmVAIFWAvZrrYRaP9/IAeTGG29M3//+9495d8uXL08333xz2rhxY5o4cWL9/MyZM9OqVavStGnT0qRJk9KGDRvS2rVr65+vXLnyiBjdWtu6DPet0NDtSllHgECTBAwgTaqWXAkQaCVgv9ZKqPXzjRxAjve27rrrrvTwhz88nXXWWemWW26pl2zZsiWNHz8+zZo1K1XDycHH7Nmz04oVK9LmzZvT2LFju7q2dQn+b4WG7kTLWgIEmiJgAGlKpeRJgEA7AvZr7SideM2gGUA++tGPpqlTp6bLL788veUtb6nfdXUZ1uLFi9P27dvTqFGjDkns2LEjjR49Os2bNy8tWrSoq2s7KZGG7kTLWgIEmiJgAGlKpeRJgEA7AvZr7SidIgPI8573vPSlL30p7dq1Kw0fPrx+11OmTKlPOfbs2XOMQrVmwoQJad26dV1d20mJNHQnWtYSINAUAQNIUyolTwIE2hGwX2tH6RQYQH70ox+lRz3qUen5z39+uvbaaw+94+ryqmHDhtX3hBz9qO4R2bdvX32ZVvXo1tqB+Hfv3p2qfw5/bN26tT7FOfwelniJRSBAgMCvV8AA8uv19+oECOQVMIDEPQfFJVjvfve70zve8Y70yU9+Ml100UWHVMaMGVOfhlQ3nh/9qG5I37t3b9q2bVv9VLfWDlSiBQsWpIULFx73aQNIvLFFIECgHAEDSDm1kAkBAnEBA0jccFAMIE984hPry6yqE4XTTjvtkEq3TjU6iTtQiZyAxJtXBAIEmiFgAGlGnWRJgEB7AgaQ9pxOtKrxA8hNN92UzjnnnON+rK57QOINIgIBAgSiAgaQqKDfJ0CgJAEDSLwajR9A3vSmN9Xf5/H1r3+9HkQOf1SfcrVkyZK2PgWrW2s7KZGG7kTLWgIEmiJgAGlKpeRJgEA7AvZr7SideE2jB5B77rkn/e7v/m562MMelm699dZj3mn1CVjVJ10N9D0gmzZtSuPGjat/r1trOymRhu5Ey1oCBJoiYABpSqXkSYBAOwL2a+0oDeIB5Oqrr04vfelLU3UT+tve9rbjvtMZM2akNWvW1N+EPnny5LR+/fr6m9CnT5+eVq9efcTvdGttu2XS0O1KWUeAQJMEDCBNqpZcCRBoJWC/1kqo9fONPgG58MIL64/d3blzZ30ScrzH/v3709KlS1N/f3/9HSEjRoxIfX19ae7cuWno0KFH/Eq31rYuw30rNHS7UtYRINAkAQNIk6olVwIEWgnYr7USav18oweQ1m+vWSs0dLPqJVsCBNoTMIC052QVAQLNELBfi9fJABI3zBZBQ2ejFIgAgYIEDCAFFUMqBAiEBezXwoTJABI3zBZBQ2ejFIgAgYIEDCAFFUMqBAiEBezXwoQGkDhhvggaOp+lSAQIlCNgACmnFjIhQCAuYL8WN3QCEjfMFkFDZ6MUiACBggQMIAUVQyoECIQF7NfChE5A4oT5ImjofJYiESBQjoABpJxayIQAgbiA/Vrc0AlI3DBbBA2djVIgAgQKEjCAFFQMqRAgEBawXwsTOgGJE+aLoKHzWYpEgEA5AgaQcmohEwIE4gL2a3FDJyBxw2wRNHQ2SoEIEChIwABSUDGkQoBAWMB+LUzoBCROmC+Chs5nKRIBAuUIGEDKqYVMCBCIC9ivxQ2dgMQNs0XQ0NkoBSJAoCABA0hBxZAKAQJhAfu1MKETkDhhvggaOp+lSAQIlCNgACmnFjIhQCAuYL8WN3QCEjfMFkFDZ6MUiACBggQMIAUVQyoECIQF7NfChE5A4oT5ImjofJYiESBQjoABpJxayIQAgbiA/Vrc0AlI3DBbBA2djVIgAgQKEjCAFFQMqRAgEBawXwsTOgGJE+aLoKHzWYpEgEA5AgaQcmohEwIE4gL2a3FDJyBxw2wRNHQ2SoEIEChIwABSUDGkQoBAWMB+LUzoBCROmC+Chs5nKRIBAuUIGEDKqYVMCBCIC9ivxQ2dgMQNs0XQ0NkoBSJAoCABA0hBxZAKAQJhAfu1MKETkDhhvggaOp+lSAQIlCNgACmnFjIhQCAuYL8WN3QCEjfMFkFDZ6MUiACBggQMIAUVQyoECIQF7NfChE5A4oT5ImjofJYiESBQjoABpJxayIQAgbiA/Vrc0AlI3DBbBA2djVIgAgQKEjCAFFQMqRAgEBawXwsTOgGJE+aLoKHzWYpEgEA5AgaQcmohEwIE4gL2a3FDJyBxw2wRNHQ2SoEIEChIwABSUDGkQoBAWMB+LUzoBCROmC+Chs5nKRIBAuUIGEDKqYVMCBCIC9ivxQ2dgMQNs0XQ0NkoBSJAoCABA0hBxZAKAQJhAfu1MKETkDhhvggaOp+lSAQIlCNgACmnFjIhQCAuYL8WN3QCEjfMFkFDZ6MUiACBggQMIAUVQyoECIQF7NfChE5A4oT5ImjofJYiESBQjoABpJxayIQAgbiA/Vrc0AlI3DBbBA2djVIgAgQKEjCAFFQMqRAgEBawXwsTOgGJE+aLoKHzWYpEgEA5AgaQcmohEwIE4gL2a3FDJyBxw2wRNHQ2SoEIEChIwABSUDGkQoBAWMB+LUzoBCROmC+Chs5nKRIBAuUIGEDKqYVMCBCIC9ivxQ2dgMQNs0XQ0NkoBSJAoCABA0hBxZAKAQJhAfu1MKETkDhhvggaOp+lSAQIlCNgACmnFjIhQCAuYL8WN3QCEjfMFkFDZ6MUiACBggQMIAUVQyoECIQF7NfChE5A4oT5ImjofJYiESBQjoABpJxayIQAgbiA/Vrc0AlI3DBbBA2djVIgAgQKEjCAFFQMqRAgEBawXwsTOgGJE+aLoKHzWYpEgEA5AgaQcmohEwIE4gL2a3FDJyBxw2wRNHQ2SoEIEChIwABSUDGkQoBAWMB+LUzoBCROmC+Chs5nKRIBAuUIGEDKqYVMCBCIC9ivxQ0bfQKye/fu9M53vjN97nOfSz/+8Y/TQx/60HTOOeekNWvWpOHDh9c6Bw4cSMuWLUv9/f1p586daeTIkamvry/NmTMnDRky5AjBbq1tt0waul0p6wgQaJKAAaRJ1ZIrAQKtBOzXWgm1fr6xA8j3vve99MxnPjOdfvrp6eKLL04jRoxId9xxR7rxxhvrgeMxj3lM/e5nzpyZVq1alaZNm5YmTZqUNmzYkNauXVv/fOXKlUcIdWtt6zLct0JDtytlHQECTRIwgDSpWnIlQKCVgP1aK6HWzzdyALn33nvT05/+9Pp04ytf+Up64AMfeNx3umXLljR+/Pg0a9astHz58kNrZs+enVasWJE2b96cxo4dW/+8W2tbl+D/VmjoTrSsJUCgKQIGkKZUSp4ECLQjYL/WjtKJ1zRyALn++uvTc57znHTNNdekF7zgBel//ud/0m/8xm+kYcOGHfFu58+fnxYvXpy2b9+eRo0adei5HTt2pNGjR6d58+alRYsW1T/v1tpOSqShO9GylgCBpggYQJpSKXkSINCOgP1aO0qDcACZO3dufZnVl7/85XqIqC6rut/97lefirz3ve9N5513Xv2up0yZUp9y7Nmz5xiF6h6RCRMmpHXr1nV1bScl0tCdaFlLgEBTBAwgTamUPAkQaEfAfq0dpUE4gLz4xS9On/nMZ9KZZ56ZnvGMZ6RXvepV6fbbb69vSL/rrrvSN77xjfrSquqf6lRk48aNxyhMnDgx7du3r770qnp0a+1A/NUN9NU/hz+2bt2apk6dWudb5edBgACBwSBgABkMVfQeCBA4KGAAifdCIy/Beu5zn5uuu+66+jKsf/3Xfz2k8NWvfrW+Mf1lL3tZ+sQnPpHGjBlTfxpWdUJy9KO6IX3v3r1p27Zt9VPdWjtQiRYsWJAWLlx43KcNIPHGFoEAgXIEDCDl1EImBAjEBQwgccNGDiAvfOEL07XXXps+/OEPp9e+9rVHKJx11lnp7rvvrj+Wt1unGp3EHahETkDizSsCAQLNEDCANKNOsiRAoD0BA0h7Tida1cgBZMaMGfV3ffzLv/xL+sM//MMj3t+5555bf5ztPffc4x6QeH+IQIAAgbCAASRMKAABAgUJGEDixWjkAPKhD32o/jLBD37wg/X/PPxRfR9I9di1a1d9g/qSJUva+hSsbq3tpEQauhMtawkQaIqAAaQplZInAQLtCNivtaN04jWNHEB+8pOfpEc/+tH1JVbr168/9I3m1WVZ1eVZr3/96+tvPq8+Aav6pKuBvgdk06ZNady4cbVQt9Z2UiIN3YmWtQQINEXAANKUSsmTAIF2BOzX2lEahANI9ZaqLxa87LLL6k/BevnLX55+9KMfpfe///3pjDPOqD9FauTIkfU7P3i5VvVN6JMnT64Hluqb0KdPn55Wr159hE631rZbJg3drpR1BAg0ScAA0qRqyZUAgVYC9muthFo/38gTkINv66qrrqq/9+PWW29ND3jAA9IFF1xQX3JVfcngwcf+/fvT0qVL6xOR6rKs6hKt6rKt6rtEhg4deoRQt9a2LsN9KzR0u1LWESDQJAEDSJOqJVcCBFoJ2K+1Emr9fKMHkNZvr1krNHSz6iVbAgTaEzCAtOdkFQECzRCwX4vXyQASN8wWQUNnoxSIAIGCBAwgBRVDKgQIhAXs18KEyQASN8wWQUNnoxSIAIGCBAwgBRVDKgQIhAXs18KEBpA4Yb4IGjqfpUgECJQjYAAppxYyIUAgLmC/Fjd0AhI3zBZBQ2ejFIgAgYIEDCAFFUMqBAiEBezXwoROQOKE+SJo6HyWIhEgUI6AAaScWsiEAIG4gP1a3NAJSNwwWwQNnY1SIAIEChIwgBRUDKkQIBAWsF8LEzoBiRPmi6Ch81mKRIBAOQIGkHJqIRMCBOIC9mtxQycgccNsETR0NkqBCBAoSMAAUlAxpEKAQFjAfi1M6AQkTpgvgobOZykSAQLlCBhAyqmFTAgQiAvYr8UNnYDEDbNF0NDZKAUiQKAgAQNIQcWQCgECYQH7tTChE5A4Yb4IGjqfpUgECJQjYAAppxYyIUAgLmC/Fjd0AhI3zBZBQ2ejFIgAgYIEDCAFFUMqBAiEBezXwoROQOKE+SJo6HyWIhEgUI6AAaScWsiEAIG4gP1a3NAJSNwwWwQNnY1SIAIEChIwgBRUDKkQIBAWsF/NnqRZAAAgAElEQVQLEzoBiRPmi6Ch81mKRIBAOQIGkHJqIRMCBOIC9mtxQycgccNsETR0NkqBCBAoSMAAUlAxpEKAQFjAfi1M6AQkTpgvgobOZykSAQLlCBhAyqmFTAgQiAvYr8UNnYDEDbNF0NDZKAUiQKAgAQNIQcWQCgECYQH7tTChE5A4Yb4IGjqfpUgECJQjYAAppxYyIUAgLmC/Fjd0AhI3zBZBQ2ejFIgAgYIEDCAFFUMqBAiEBezXwoROQOKE+SJo6HyWIhEgUI6AAaScWsiEAIG4gP1a3NAJSNwwWwQNnY1SIAIEChIwgBRUDKkQIBAWsF8LEzoBiRPmi6Ch81mKRIBAOQIGkHJqIRMCBOIC9mtxQycgccNsETR0NkqBCBAoSMAAUlAxpEKAQFjAfi1M6AQkTpgvgobOZykSAQLlCBhAyqmFTAgQiAvYr8UNnYDEDbNF0NDZKAUiQKAgAQNIQcWQCgECYQH7tTChE5A4Yb4IGjqfpUgECJQjYAAppxYyIUAgLmC/Fjd0AhI3zBZBQ2ejFIgAgYIEDCAFFUMqBAiEBezXwoROQOKE+SJo6HyWIhEgUI6AAaScWsiEAIG4gP1a3NAJSNwwWwQNnY1SIAIEChIwgBRUDKkQIBAWsF8LEzoBiRPmi6Ch81mKRIBAOQIGkHJqIRMCBOIC9mtxQycgccNsETR0NkqBCBAoSMAAUlAxpEKAQFjAfi1M6AQkTpgvgobOZykSAQLlCBhAyqmFTAgQiAvYr8UNnYDEDbNF0NDZKAUiQKAgAQNIQcWQCgECYQH7tTChE5A4Yb4IGjqfpUgECJQjYAAppxYyIUAgLmC/Fjd0AhI3zBZBQ2ejFIgAgYIEDCAFFUMqBAiEBezXwoROQOKE+SJo6HyWIhEgUI6AAaScWsiEAIG4gP1a3NAJSNwwWwQNnY1SIAIEChIwgBRUDKkQIBAWsF8LEzoBiRPmi6Ch81mKRIBAOQIGkHJqIRMCBOIC9mtxQycgccNsETR0NkqBCBAoSMAAUlAxpEKAQFjAfi1M6AQkTpgvgobOZykSAQLlCBhAyqmFTAgQiAvYr8UNG3kCctttt6VRo0Yd992//vWvT/39/YeeO3DgQFq2bFn9s507d6aRI0emvr6+NGfOnDRkyJAjYnRrbbtl0tDtSllHgECTBAwgTaqWXAkQaCVgv9ZKqPXzjR5ALrzwwnTRRRcd8S7PPvvsdO655x762cyZM9OqVavStGnT0qRJk9KGDRvS2rVrU/XzlStXHvG73Vrbugz3rdDQ7UpZR4BAkwQMIE2qllwJEGglYL/WSqj1840eQObPn5/e9a53Dfgut2zZksaPH59mzZqVli9ffmjd7Nmz04oVK9LmzZvT2LFj6593a23rEvzfCg3diZa1BAg0RcAA0pRKyZMAgXYE7NfaUTrxmsYPINUQUj1OP/30Y95p9dzixYvT9u3bj7hka8eOHWn06NFp3rx5adGiRfXvdWttJyXS0J1oWUuAQFMEDCBNqZQ8CRBoR8B+rR2lQTyAPPCBD0x33nln/Q6rS68uu+yydOmllx56x1OmTKlPOfbs2XOMwvDhw9OECRPSunXr6ue6tbaTEmnoTrSsJUCgKQIGkKZUSp4ECLQjYL/WjtIgHEB++MMfpte97nXpxS9+cXr0ox+dbr/99vom85tvvjm99a1vrW86rx7V5VXDhg1LGzduPEZh4sSJad++ffWlV91cOxD/7t27U/XP4Y+tW7emqVOn1vlW+XkQIEBgMAgYQAZDFb0HAgQOChhA4r3QyEuwjve2q0+w+oM/+IN0ww03pO9+97tpzJgx9T/VSUd14/nRj+qG9L1796Zt27bVT3Vr7UAlWrBgQVq4cOFxnzaAxBtbBAIEyhEwgJRTC5kQIBAXMIDEDQfNAFJRXHPNNelFL3pRWrNmTXrjG9/oBCTeHyIQIEAgLGAACRMKQIBAQQIGkHgxBtUAcsstt9SfelXdWF7dYN6t+zo6idtJiTR0J1rWEiDQFAEDSFMqJU8CBNoRsF9rR+nEawbVAHL11Venl770pfX9INUXElZDyJIlS9r6FKxure2kRBq6Ey1rCRBoioABpCmVkicBAu0I2K+1ozQIB5D//M//TA95yEOOeGd33313mjx5cvrWt76Vvv/979ffeF59Alb1SVcDfQ/Ipk2b0rhx4+o43VrbSYk0dCda1hIg0BQBA0hTKiVPAgTaEbBfa0dpEA4gL3nJS9Jdd91Vf+P5iBEj6k/B+shHPlKfdFQnHm9/+9sPvesZM2bU94RU34ReDSjr16+vvwl9+vTpafXq1UfodGttu2XS0O1KWUeAQJMEDCBNqpZcCRBoJWC/1kqo9fONvATrQx/6UD1wfOc730nVaUj1fSDVx9ZW33Be3YR++GP//v1p6dKl9WVZu3btqgeWvr6+NHfu3DR06NCerG1dhvtWaOh2pawjQKBJAgaQJlVLrgQItBKwX2sl1Pr5Rg4grd9WM1do6GbWTdYECJxYwACiQwgQGEwC9mvxahpA4obZImjobJQCESBQkIABpKBiSIUAgbCA/VqYMBlA4obZImjobJQCESBQkIABpKBiSIUAgbCA/VqY0AASJ8wXQUPnsxSJAIFyBAwg5dRCJgQIxAXs1+KGTkDihtkiaOhslAIRIFCQgAGkoGJIhQCBsID9WpjQCUicMF8EDZ3PUiQCBMoRMICUUwuZECAQF7Bfixs6AYkbZougobNRCkSAQEECBpCCiiEVAgTCAvZrYUInIHHCfBE0dD5LkQgQKEfAAFJOLWRCgEBcwH4tbugEJG6YLYKGzkYpEAECBQkYQAoqhlQIEAgL2K+FCZ2AxAnzRdDQ+SxFIkCgHAEDSDm1kAkBAnEB+7W4oROQuGG2CBo6G6VABAgUJGAAKagYUiFAICxgvxYmdAISJ8wXQUPnsxSJAIFyBAwg5dRCJgQIxAXs1+KGTkDihtkiaOhslAIRIFCQgAGkoGJIhQCBsID9WpjQCUicMF8EDZ3PUiQCBMoRMICUUwuZECAQF7Bfixs6AYkbZougobNRCkSAQEECBpCCiiEVAgTCAvZrYUInIHHCfBE0dD5LkQgQKEfAAFJOLWRCgEBcwH4tbugEJG6YLYKGzkYpEAECBQkYQAoqhlQIEAgL2K+FCZ2AxAnzRdDQ+SxFIkCgHAEDSDm1kAkBAnEB+7W4oROQuGG2CBo6G6VABAgUJGAAKagYUiFAICxgvxYmdAISJ8wXQUPnsxSJAIFyBAwg5dRCJgQIxAXs1+KGTkDihtkiaOhslAIRIFCQgAGkoGJIhQCBsID9WpjQCUicMF8EDZ3PUiQCBMoRMICUUwuZECAQF7Bfixs6AYkbZougobNRCkSAQEECBpCCiiEVAgTCAvZrYUInIHHCfBE0dD5LkQgQKEfAAFJOLWRCgEBcwH4tbugEJG6YLYKGzkYpEAECBQkYQAoqhlQIEAgL2K+FCZ2AxAnzRdDQ+SxFIkCgHAEDSDm1kAkBAnEB+7W4oROQuGG2CBo6G6VABAgUJGAAKagYUiFAICxgvxYmdAISJ8wXQUPnsxSJAIFyBAwg5dRCJgQIxAXs1+KGTkDihtkiaOhslAIRIFCQgAGkoGJIhQCBsID9WpjQCUicMF8EDZ3PUiQCBMoRMICUUwuZECAQF7Bfixs6AYkbZougobNRCkSAQEECBpCCiiEVAgTCAvZrYUInIHHCfBE0dD5LkQgQKEfAAFJOLWRCgEBcwH4tbugEJG6YLYKGzkYpEAECBQkYQAoqhlQIEAgL2K+FCZ2AxAnzRdDQ+SxFIkCgHAEDSDm1kAkBAnEB+7W4oROQuGG2CBo6G6VABAgUJGAAKagYUiFAICxgvxYmdAISJ8wXQUPnsxSJAIFyBAwg5dRCJgQIxAXs1+KGTkDihtkiaOhslAIRIFCQgAGkoGJIhQCBsID9WpjQCUicMF8EDZ3PUiQCBMoRMICUUwuZECAQF7Bfixs6AYkbZougobNRCkSAQEECBpCCiiEVAgTCAvZrYUInIHHCfBE0dD5LkQgQKEfAAFJOLWRCgEBcwH4tbugEJG6YLYKGzkYpEAECBQkYQAoqhlQIEAgL2K+FCZ2AxAnzRdDQ+SxFIkCgHAEDSDm1kAkBAnEB+7W44aA5Abn++uvTc57znFrke9/7Xjr77LMP6Rw4cCAtW7Ys9ff3p507d6aRI0emvr6+NGfOnDRkyJAjFLu1tp1Saeh2lKwhQKBpAgaQplVMvgQInEjAfi3eH4NiALnnnnvS+PHj6+Hil7/85TEDyMyZM9OqVavStGnT0qRJk9KGDRvS2rVrU/XzlStXHqHYrbXtlEpDt6NkDQECTRMwgDStYvIlQMAA0t0eGBQDyJIlS9Jf//Vfp1e/+tX1/zz8BGTLli31cDJr1qy0fPnyQ5qzZ89OK1asSJs3b05jx46tf96tte2W0ADSrpR1BAg0ScAA0qRqyZUAgVYC9muthFo/3/gB5Ac/+EF64hOfmD7wgQ+k6n9fuHDhEQPI/Pnz0+LFi9P27dvTqFGjDons2LEjjR49Os2bNy8tWrSo/nm31rYuw30rNHS7UtYRINAkAQNIk6olVwIEWgnYr7USav184weQCy+8MO3du7e+rKoaPo4eQKZMmVKfcuzZs+cYjeHDh6cJEyakdevW1c91a23rMhhA2jWyjgCB5gkYQJpXMxkTIDCwgAEk3h2NHkCuvfbaVA0g3/jGN9JTn/rUtGDBgmMGkOryqmHDhqWNGzceozVx4sS0b9+++tKr6tGttccr0+7du1P1z+GPrVu3pqlTp9a5Vrl5ECBAYDAIGEAGQxW9BwIEDgoYQOK90NgB5O67705PetKT0gUXXJDWrFlTSxxvABkzZkyqTjqqE5KjH9UN6dXpybZt2+qnurX2eGU6mOvxnjOAxBtbBAIEyhEwgJRTC5kQIBAXMIDEDRs7gPzFX/xFuuKKK9J3v/vddOaZZw44gHTrVKOTuMcrkxOQePOKQIBAMwQMIM2okywJEGhPwADSntOJVjVyALn99tvrG8gvu+yy+vs8Dj7e//73159sdd1116WzzjqrXtOt+zo6idtumTR0u1LWESDQJAEDSJOqJVcCBFoJ2K+1Emr9fCMHkE2bNtU3j5/occYZZ6Q777yz/pSr6mN62/kUrG6tbV2G+1Zo6HalrCNAoEkCBpAmVUuuBAi0ErBfayXU+vlGDiC/+MUv0he/+MVj3t0nPvGJ9MlPfrI+BRkxYkR68YtfXH8CVjWsDPQ9INUwM27cuDpWt9a2LoMBpF0j6wgQaJ6AAaR5NZMxAQIDCxhA4t3RyAFkoLd9vJvQq7UzZsyob1Svvgl98uTJaf369fU3oU+fPj2tXr36iHDdWttOqTR0O0rWECDQNAEDSNMqJl8CBE4kYL8W749TYgDZv39/Wrp0aerv70+7du2qT0eqe0fmzp2bhg4deoRit9a2UyoN3Y6SNQQINE3AANK0ismXAAEDSHd7YFANIN2l6n50A0j3jb0CAQK9FzCA9N7cKxIg0D0B+7W4rQEkbpgtgobORikQAQIFCRhACiqGVAgQCAvYr4UJkwEkbpgtgobORikQAQIFCRhACiqGVAgQCAvYr4UJDSBxwnwRNHQ+S5EIEChHwABSTi1kQoBAXMB+LW7oBCRumC2Chs5GKRABAgUJGEAKKoZUCBAIC9ivhQmdgMQJ80XQ0PksRSJAoBwBA0g5tZAJAQJxAfu1uKETkLhhtggaOhulQAQIFCRgACmoGFIhQCAsYL8WJnQCEifMF0FD57MUiQCBcgQMIOXUQiYECMQF7Nfihk5A4obZImjobJQCESBQkIABpKBiSIUAgbCA/VqY0AlInDBfBA2dz1IkAgTKETCAlFMLmRAgEBewX4sbOgGJG2aLoKGzUQpEgEBBAgaQgoohFQIEwgL2a2FCJyBxwnwRNHQ+S5EIEChHwABSTi1kQoBAXMB+LW7oBCRumC2Chs5GKRABAgUJGEAKKoZUCBAIC9ivhQmdgMQJ80XQ0PksRSJAoBwBA0g5tZAJAQJxAfu1uKETkLhhtggaOhulQAQIFCRgACmoGFIhQCAsYL8WJnQCEifMF0FD57MUiQCBcgQMIOXUQiYECMQF7Nfihk5A4obZImjobJQCESBQkIABpKBiSIUAgbCA/VqY0AlInDBfBA2dz1IkAgTKETCAlFMLmRAgEBewX4sbOgGJG2aLoKGzUQpEgEBBAgaQgoohFQIEwgL2a2FCJyBxwnwRNHQ+S5EIEChHwABSTi1kQoBAXMB+LW7oBCRumC2Chs5GKRABAgUJGEAKKoZUCBAIC9ivhQmdgMQJ80XQ0PksRSJAoBwBA0g5tZAJAQJxAfu1uKETkLhhtggaOhulQAQIFCRgACmoGFIhQCAsYL8WJnQCEifMF0FD57MUiQCBcgQMIOXUQiYECMQF7Nfihk5A4obZImjobJQCESBQkIABpKBiSIUAgbCA/VqY0AlInDBfBA2dz1IkAgTKETCAlFMLmRAgEBewX4sbOgGJG2aLoKGzUQpEgEBBAgaQgoohFQIEwgL2a2FCJyBxwnwRNHQ+S5EIEChHwABSTi1kQoBAXMB+LW7oBCRumC2Chs5GKRABAgUJGEAKKoZUCBAIC9ivhQmdgMQJ80XQ0PksRSJAoBwBA0g5tZAJAQJxAfu1uKETkLhhtggaOhulQAQIFCRgACmoGFIhQCAsYL8WJnQCEifMF0FD57MUiQCBcgQMIOXUQiYECMQF7Nfihk5A4obZImjobJQCESBQkIABpKBiSIUAgbCA/VqY0AlInDBfBA2dz1IkAgTKETCAlFMLmRAgEBewX4sbOgGJG2aLoKGzUQpEgEBBAgaQgoohFQIEwgL2a2FCJyBxwnwRNHQ+S5EIEChHwABSTi1kQoBAXMB+LW7oBCRumC2Chs5GKRABAgUJGEAKKoZUCBAIC9ivhQmdgMQJ80XQ0PksRSJAoBwBA0g5tZAJAQJxAfu1uKETkLhhtggaOhulQAQIFCRgACmoGFIhQCAsYL8WJnQCEifMF0FD57MUiQCBcgQMIOXUQiYECMQF7Nfihk5A4obZImjobJQCESBQkIABpKBiSIUAgbCA/VqY0AlInDBfBA2dz1IkAgTKETCAlFMLmRAgEBewX4sbNvIEZOvWrWnhwoVp48aNaffu3ek3fuM30ujRo9PFF1+cLrnkknTaaacdkjlw4EBatmxZ6u/vTzt37kwjR45MfX19ac6cOWnIkCFHCHZrbbtl0tDtSllHgECTBAwgTaqWXAkQaCVgv9ZKqPXzjRxAvvCFL6TLL788Pf3pT08jRoxI1eCwfv369LGPfSz98R//cbrmmmsOvfOZM2emVatWpWnTpqVJkyalDRs2pLVr16bq5ytXrjxCqFtrW5fhvhUaul0p6wgQaJKAAaRJ1ZIrAQKtBOzXWgm1fr6RA8hAb+tNb3pTPVR8+9vfTo973OPSli1b0vjx49OsWbPS8uXLD/3a7Nmz04oVK9LmzZvT2LFj6593a23rEvzfCg3diZa1BAg0RcAA0pRKyZMAgXYE7NfaUTrxmkE1gFSnItWlVV/72tfq05H58+enxYsXp+3bt6dRo0YdktixY0d9yda8efPSokWL6p93a20nJdLQnWhZS4BAUwQMIE2plDwJEGhHwH6tHaVBPIDcddddqfrnl7/8ZfrGN76RLr300nT/+98/bdu2LZ1++ulpypQp9SnHnj17jlEYPnx4mjBhQlq3bl39XLfWdlIiDd2JlrUECDRFwADSlErJkwCBdgTs19pRGsQDyIIFC+qb0Q8+qlOPK6+8Mo0bN67+UXV51bBhw+qb1Y9+TJw4Me3bt6++9Kqbawfir26er/45/FHdXD916tQ63yo/DwIECAwGAQPIYKii90CAwEEBA0i8Fxp9CVZ1aVX1z09/+tN0/fXXp1tuuSUtWbIkPfvZz65lxowZk6qTjurG86Mf1Q3pe/furU9Lurl2oBIdPTwdvs4AEm9sEQgQKEfAAFJOLWRCgEBcwAASN2z0AHL023/f+96X3va2t9WXXT3hCU9wAhLvDxEIECAQFjCAhAkFIECgIAEDSLwYg2oA+fGPf5we/vCHH7q5vFv3dXQSt5MSaehOtKwlQKApAgaQplRKngQItCNgv9aO0onXDKoB5Lbbbqs/7ar6MsIrrriiHkSqS7La+RSsbq3tpEQauhMtawkQaIqAAaQplZInAQLtCNivtaM0CAeQ6t6Nhz3sYce8s7lz59bfev7hD384vfa1r60vxao+6Wqg7wHZtGnToRvWu7W2kxJp6E60rCVAoCkCBpCmVEqeBAi0I2C/1o7SIBxA/uRP/qS+8by62XzkyJHp5z//efr85z+frrvuunT++eenL33pS2no0KH1O58xY0Zas2ZN/U3okydPrr8xvfom9OnTp6fVq1cfodOtte2WSUO3K2UdAQJNEjCANKlaciVAoJWA/VorodbPN/ISrI9//OP1KUf1qVd33HFHOu2009LjH//49IpXvKI+7aj+74OP/fv3p6VLl6b+/v60a9euNGLEiNTX15eq05KDQ0q317Yuw30rNHS7UtYRINAkAQNIk6olVwIEWgnYr7USav18IweQ1m+rmSs0dDPrJmsCBE4sYADRIQQIDCYB+7V4NQ0gccNsETR0NkqBCBAoSMAAUlAxpEKAQFjAfi1MmAwgccNsETR0NkqBCBAoSMAAUlAxpEKAQFjAfi1MaACJE+aLoKHzWYpEgEA5AgaQcmohEwIE4gL2a3FDJyBxw2wRNHQ2SoEIEChIwABSUDGkQoBAWMB+LUzoBCROmC+Chs5nKRIBAuUIGEDKqYVMCBCIC9ivxQ2dgMQNs0XQ0NkoBSJAoCABA0hBxZAKAQJhAfu1MKETkDhhvggaOp+lSAQIlCNgACmnFjIhQCAuYL8WN3QCEjfMFkFDZ6MUiACBggQMIAUVQyoECIQF7NfChE5A4oT5ImjofJYiESBQjoABpJxayIQAgbiA/Vrc0AlI3DBbBA2djVIgAgQKEjCAFFQMqRAgEBawXwsTOgGJE+aLoKHzWYpEgEA5AgaQcmohEwIE4gL2a3FDJyBxw2wRNHQ2SoEIEChIwABSUDGkQoBAWMB+LUzoBCROmC+Chs5nKRIBAuUIGEDKqYVMCBCIC9ivxQ2dgMQNs0XQ0NkoBSJAoCABA0hBxZAKAQJhAfu1MKETkDhhvggaOp+lSAQIlCNgACmnFjIhQCAuYL8WN3QCEjfMFkFDZ6MUiACBggQMIAUVQyoECIQF7NfChE5A4oT5ImjofJYiESBQjoABpJxayIQAgbiA/Vrc0AlI3DBbBA2djVIgAgQKEjCAFFQMqRAgEBawXwsTOgGJE+aLoKHzWYpEgEA5AgaQcmohEwIE4gL2a3FDJyBxw2wRNHQ2SoEIEChIwABSUDGkQoBAWMB+LUzoBCROmC+Chs5nKRIBAuUIGEDKqYVMCBCIC9ivxQ2dgMQNs0XQ0NkoBSJAoCABA0hBxZAKAQJhAfu1MKETkDhhvggaOp+lSAQIlCNgACmnFjIhQCAuYL8WN3QCEjfMFkFDZ6MUiACBggQMIAUVQyoECIQF7NfChE5A4oT5ImjofJYiESBQjoABpJxayIQAgbiA/Vrc0AlI3DBbBA2djVIgAgQKEjCAFFQMqRAgEBawXwsTOgGJE+aLoKHzWYpEgEA5AgaQcmohEwIE4gL2a3FDJyBxw2wRNHQ2SoEIEChIwABSUDGkQoBAWMB+LUzoBCROmC+Chs5nKRIBAuUIGEDKqYVMCBCIC9ivxQ2dgMQNs0XQ0NkoBSJAoCABA0hBxZAKAQJhAfu1MKETkDhhvggaOp+lSAQIlCNgACmnFjIhQCAuYL8WN3QCEjfMFkFDZ6MUiACBggQMIAUVQyoECIQF7NfChE5A4oT5ImjofJYiESBQjoABpJxayIQAgbiA/Vrc0AlI3DBbBA2djVIgAgQKEjCAFFQMqRAgEBawXwsTOgGJE+aLoKHzWYpEgEA5AgaQcmohEwIE4gL2a3FDJyBxw2wRNHQ2SoEIEChIwABSUDGkQoBAWMB+LUzoBCROmC+Chs5nKRIBAuUIGEDKqYVMCBCIC9ivxQ2dgMQNs0XQ0NkoBSJAoCABA0hBxZAKAQJhAfu1MKETkDhhvggaOp+lSAQIlCNgACmnFjIhQCAuYL8WN3QCEjfMFkFDZ6MUiACBggQMIAUVQyoECIQF7NfChE5A4oT5ImjofJYiESBQjoABpJxayIQAgbiA/VrcsJEnIDfffHO66qqr0vXXX5927NiRzjjjjPSkJz0pveMd70jPfe5zj1A5cOBAWrZsWerv7087d+5MI0eOTH19fWnOnDlpyJAhPVnbbpk0dLtS1hEg0CQBA0iTqiVXAgRaCdivtRJq/XwjB5CLLroofeUrX0kvfelL08SJE9Odd96Z1q5dm771rW+lK664Il1yySWH3vnMmTPTqlWr0rRp09KkSZPShg0b6rXVz1euXHmEULfWti7DfSs0dLtS1hEg0CQBA0iTqiVXAgRaCdivtRJq/XwjB5D169enpz3taem000479A7vvvvu9JSnPCXdcccdae/evWno0KFpy5Ytafz48WnWrFlp+fLlh9bOnj07rVixIm3evDmNHTu2/nm31rYuwf+t0NCdaFlLgEBTBAwgTamUPAkQaEfAfq0dpROvaeQAMtBbestb3pLe+973ph/+8If1pVbz589PixcvTtu3b0+jRo069GvVZVujR49O8+bNS4sWLap/3q21nZRIQ3eiZS0BAk0RMIA0pVLyJECgHQH7tXaUTqEB5FWvelX61Kc+lX7+85/X94VMmTKlPuXYs2fPMQrDhw9PEyZMSOvWrauf69baTkqkoTvRspYAgaYIGECaUil5EiDQjoD9WjtKp8gAsl7QbIIAAB+oSURBVHXr1voSrBe84AXp05/+dP2uq8urhg0bljZu3HiMQnXvyL59++pLr7q5diD+3bt3p+qfwx/Ve5g6dWqdb5WfBwECBAaDgAFkMFTReyBA4KCAASTeC4PiEqxf/OIX6bzzzqs39NWJx6Me9ahaZsyYMak66ahuPD/6Ud2QXt0rsm3btq6uHahECxYsSAsXLjzu0waQeGOLQIBAOQIGkHJqIRMCBOICBpC4YeMHkOrm8+ryqZtuuqm+nOpZz3rWIRUnIPEGEYEAAQJRAQNIVNDvEyBQkoABJF6NRg8g99xzT3rRi15Ufx/I1VdfXV9+dfijW/d1dBK3kxJp6E60rCVAoCkCBpCmVEqeBAi0I2C/1o7Sidc0dgDZv39/qr4P5Jprrkkf/ehH0ytf+cpj3mn1KVdLlixp61OwurW2kxJp6E60rCVAoCkCBpCmVEqeBAi0I2C/1o7SIBxAfvWrX6VXv/rV6eMf/3i68sor0xve8IbjvsvqfpDqk64G+h6QTZs2pXHjxtW/2621nZRIQ3eiZS0BAk0RMIA0pVLyJECgHQH7tXaUBuEA8uY3vzm9733vq+/36OvrO+YdXnDBBfXN59VjxowZac2aNfU3oU+ePDlVX2JYfRP69OnT0+rVq4/43W6tbbdMGrpdKesIEGiSgAGkSdWSKwECrQTs11oJtX6+kZdgPfvZz05f+cpXBnx3X/rSl1K1pnpUl2otXbo09ff3p127dqURI0bUQ8vcuXPrb0s//NGtta3LcN8KDd2ulHUECDRJwADSpGrJlQCBVgL2a62EWj/fyAGk9dtq5goN3cy6yZoAgRMLGEB0CAECg0nAfi1eTQNI3DBbBA2djVIgAgQKEjCAFFQMqRAgEBawXwsTJgNI3DBbBA2djVIgAgQKEjCAFFQMqRAgEBawXwsTGkDihPkiaOh8liIRIFCOgAGknFrIhACBuID9WtzQCUjcMFsEDZ2NUiACBAoSMIAUVAypECAQFrBfCxM6AYkT5ougofNZikSAQDkCBpByaiETAgTiAvZrcUMnIHHDbBE0dDZKgQgQKEjAAFJQMaRCgEBYwH4tTOgEJE6YL4KGzmcpEgEC5QgYQMqphUwIEIgL2K/FDZ2AxA2zRdDQ2SgFIkCgIAEDSEHFkAoBAmEB+7UwoROQOGG+CBo6n6VIBAiUI2AAKacWMiFAIC5gvxY3dAISN8wWQUNnoxSIAIGCBAwgBRVDKgQIhAXs18KETkDihPkiaOh8liIRIFCOgAGknFrIhACBuID9WtzQCUjcMFsEDZ2NUiACBAoSMIAUVAypECAQFrBfCxM6AYkT5ougofNZikSAQDkCBpByaiETAgTiAvZrcUMnIHHDbBE0dDZKgQgQKEjAAFJQMaRCgEBYwH4tTOgEJE6YL4KGzmcpEgEC5QgYQMqphUwIEIgL2K/FDZ2AxA2zRdDQ2SgFIkCgIAEDSEHFkAoBAmEB+7UwoROQOGG+CBo6n6VIBAiUI2AAKacWMiFAIC5gvxY3dAISN8wWQUNnoxSIAIGCBAwgBRVDKgQIhAXs18KETkDihPkiaOh8liIRIFCOgAGknFrIhACBuID9WtzQCUjcMFsEDZ2NUiACBAoSMIAUVAypECAQFrBfCxM6AYkT5ougofNZikSAQDkCBpByaiETAgTiAvZrcUMnIHHDbBE0dDZKgQgQKEjAAFJQMaRCgEBYwH4tTOgEJE6YL4KGzmcpEgEC5QgYQMqphUwIEIgL2K/FDZ2AxA2zRdDQ2SgFIkCgIAEDSEHFkAoBAmEB+7UwoROQOGG+CBo6n6VIBAiUI2AAKacWMiFAIC5gvxY3dAISN8wWQUNnoxSIAIGCBAwgBRVDKgQIhAXs18KETkDihPkiaOh8liIRIFCOgAGknFrIhACBuID9WtzQCUjcMFsEDZ2NUiACBAoSMIAUVAypECAQFrBfCxM6AYkT5ougofNZikSAQDkCBpByaiETAgTiAvZrcUMnIHHDbBE0dDZKgQgQKEjAAFJQMaRCgEBYwH4tTOgEJE6YL4KGzmcpEgEC5QgYQMqphUwIEIgL2K/FDZ2AxA2zRdDQ2SgFIkCgIAEDSEHFkAoBAmEB+7UwoROQOGG+CBo6n6VIBAiUI2AAKacWMiFAIC5gvxY3dAISN8wWQUNnoxSIAIGCBAwgBRVDKgQIhAXs18KETkDihPkiaOh8liIRIFCOgAGknFrIhACBuID9WtzQCUjcMFsEDZ2NUiACBAoSMIAUVAypECAQFrBfCxM6AYkT5ougofNZikSAQDkCBpByaiETAgTiAvZrcUMnIHHDbBE0dDZKgQgQKEjAAFJQMaRCgEBYwH4tTOgEJE6YL4KGzmcpEgEC5QgYQMqphUwIEIgL2K/FDZ2AxA2zRdDQ2SgFIkCgIAEDSEHFkAoBAmEB+7UwoROQOGG+CBo6n6VIBAiUI2AAKacWMiFAIC5gvxY3bOwJyJ133pkuv/zytHHjxnTzzTenPXv2pNe+9rXpwx/+8DEqBw4cSMuWLUv9/f1p586daeTIkamvry/NmTMnDRky5Ij13VrbTqk0dDtK1hAg0DQBA0jTKiZfAgROJGC/Fu+Pxg4gt912Wxo1alR6xCMekZ761Kema6+9dsABZObMmWnVqlVp2rRpadKkSWnDhg1p7dq1qfr5ypUrj1Ds1tp2SqWh21GyhgCBpgkYQJpWMfkSIGAA6W4PNHYA+d///d/0k5/8JD3ykY9M+/fvT/e///2PO4Bs2bIljR8/Ps2aNSstX778kObs2bPTihUr0ubNm9PYsWPrn3drbbslNIC0K2UdAQJNEjCANKlaciVAoJWA/VorodbPN3YAOfytnWgAmT9/flq8eHHavn17fWJy8LFjx440evToNG/evLRo0aL6x91a27oM963Q0O1KWUeAQJMEDCBNqpZcCRBoJWC/1kqo9fODfgCZMmVKfcpR3SNy9GP48OFpwoQJad26dfVT3VrbugwGkHaNrCNAoHkCBpDm1UzGBAgMLGAAiXfHoB9Aqsurhg0bVt+sfvRj4sSJad++ffWlV9WjW2uPV6bdu3en6p/DH1u3bk1Tp06tc61y8yBAgMBgEDCADIYqeg8ECBwUMIDEe2HQDyBjxoxJ1UlHdeP50Y/qhvS9e/embdu21U91a+3xyrRgwYK0cOHC41bQABJvbBEIEChHwABSTi1kQoBAXMAAEjcc9ANIt041OonrBCTeqCIQINBcAQNIc2sncwIEjhUwgMS7YtAPIN26r6OTuO2WSUO3K2UdAQJNEjCANKlaciVAoJWA/VorodbPD/oBpPqUqyVLlrT1KVjdWtu6DPet0NDtSllHgECTBAwgTaqWXAkQaCVgv9ZKqPXzg34AqT4Bq/qkq4G+B2TTpk1p3LhxtVS31rYugwGkXSPrCBBonoABpHk1kzEBAgMLGEDi3dHoAeQDH/hA+vnPf55+9atfpf/3//5fPWi85CUvqVVe9KIXHRosZsyYkdasWVN/E/rkyZPT+vXr629Cnz59elq9evURit1a206pNHQ7StYQINA0AQNI0yomXwIETiRgvxbvj0YPIGeddVb6wQ9+cFyFasC4+OKL6+eqLypcunRp6u/vT7t27UojRoxIfX19ae7cuWno0KFH/H631rZTKg3djpI1BAg0TcAA0rSKyZcAAQNId3ug0QNId2l6H90A0ntzr0iAQPcFDCDdN/YKBAj0TsB+LW5tAIkbZougobNRCkSAQEECBpCCiiEVAgTCAvZrYcJkAIkbZougobNRCkSAQEECBpCCiiEVAgTCAvZrYUIDSJwwXwQNnc9SJAIEyhEwgJRTC5kQIBAXsF+LGzoBiRtmi6Chs1EKRIBAQQIGkIKKIRUCBMIC9mthQicgccJ8ETR0PkuRCBAoR8AAUk4tZEKAQFzAfi1u6AQkbpgtgobORikQAQIFCRhACiqGVAgQCAvYr4UJnYDECfNF0ND5LEUiQKAcAQNIObWQCQECcQH7tbihE5C4YbYIGjobpUAECBQkYAApqBhSIUAgLGC/FiZ0AhInzBdBQ+ezFIkAgXIEDCDl1EImBAjEBezX4oZOQOKG2SJo6GyUAhEgUJCAAaSgYkiFAIGwgP1amNAJSJwwXwQNnc9SJAIEyhEwgJRTC5kQIBAXsF+LGzoBiRtmi6Chs1EKRIBAQQIGkIKKIRUCBMIC9mthQicgccJ8ETR0PkuRCBAoR8AAUk4tZEKAQFzAfi1u6AQkbpgtgobORikQAQIFCRhACiqGVAgQCAvYr4UJnYDECfNF0ND5LEUiQKAcAQNIObWQCQECcQH7tbihE5C4YbYIGjobpUAECBQkYAApqBhSIUAgLGC/FiZ0AhInzBdBQ+ezFIkAgXIEDCDl1EImBAjEBezX4oZOQOKG2SJo6GyUAhEgUJCAAaSgYkiFAIGwgP1amNAJSJwwXwQNnc9SJAIEyhEwgJRTC5kQIBAXsF+LGzoBiRtmi6Chs1EKRIBAQQIGkIKKIRUCBMIC9mthQicgccJ8ETR0PkuRCBAoR8AAUk4tZEKAQFzAfi1u6AQkbpgtgobORikQAQIFCRhACiqGVAgQCAvYr4UJnYDECfNF0ND5LEUiQKAcAQNIObWQCQECcQH7tbihE5C4YbYIGjobpUAECBQkYAApqBhSIUAgLGC/FiZ0AhInzBdBQ+ezFIkAgXIEDCDl1EImBAjEBezX4oZOQOKG2SJo6GyUAhEgUJCAAaSgYkiFAIGwgP1amNAJSJwwXwQNnc9SJAIEyhEwgJRTC5kQIBAXsF+LGzoBiRtmi6Chs1EKRIBAQQIGkIKKIRUCBMIC9mthQicgccJ8ETR0PkuRCBAoR8AAUk4tZEKAQFzAfi1u6AQkbpgtgobORikQAQIFCRhACiqGVAgQCAvYr4UJnYDECfNF0ND5LEUiQKAcAQNIObWQCQECcQH7tbihE5C4YbYIGjobpUAECBQkYAApqBhSIUAgLGC/FiZ0AhInzBdBQ+ezFIkAgXIEDCDl1EImBAjEBezX4oZOQOKG2SJo6GyUAhEgUJCAAaSgYkiFAIGwgP1amNAJSJwwXwQNnc9SJAIEyhEwgJRTC5kQIBAXsF+LGzoBiRtmi6Chs1EKRIBAQQIGkIKKIRUCBMIC9mthQicgccJ8ETR0PkuRCBAoR8AAUk4tZEKAQFzAfi1u6AQkbpgtgobORikQAQIFCRhACiqGVAgQCAvYr4UJnYDECfNF0ND5LEUiQKAcAQNIObWQCQECcQH7tbihE5C4YbYIGjobpUAECBQkYAApqBhSIUAgLGC/FiZ0AhInzBdBQ+ezFIkAgXIEDCDl1EImBAjEBezX4oZOQAKGBw4cSMuWLUv9/f1p586daeTIkamvry/NmTMnDRkypOPIGrpjMr9AgEADBAwgDSiSFAkQaFvAfq1tqgEXGkAChjNnzkyrVq1K06ZNS5MmTUobNmxIa9euTdXPV65c2XFkDd0xmV8gQKABAgaQBhRJigQItC1gv9Y2lQEkTnVkhC1btqTx48enWbNmpeXLlx96cvbs2WnFihVp8+bNaezYsR29rIbuiMtiAgQaImAAaUihpEmAQFsC9mttMZ1wkROQkzScP39+Wrx4cdq+fXsaNWrUoSg7duxIo0ePTvPmzUuLFi3qKLqG7ojLYgIEGiJgAGlIoaRJgEBbAvZrbTEZQOJMx0aYMmVKfcqxZ8+eY54cPnx4mjBhQlq3bl1HL62hO+KymACBhggYQBpSKGkSINCWgP1aW0wGkDjTsRGqy6uGDRuWNm7ceMyTEydOTPv27UvVZVoDPXbv3p2qfw5/bNq0Kb3+9a9PV111VXrCE57QjbTFJECAQM8F/vj9X+35a/byBT/3Z8/o5ct5LQIEfs0CW7duTVOnTk033HBDmjx58q85m2a+vEuwTrJuY8aMSdVJR3Xj+dGP6ob0vXv3pm3btg0YfcGCBWnhwoUn+ep+jQABAgQIECBA4NcpUP3B+DWvec2vM4XGvrYB5CRL140TkJ/97Gepmqqry7dOP/30k8zMr51I4OBfLZwyNbdP1LC5tasyV79m108Nm18/NYzX8O6770633XZbqi7HP/PMM+MBT8EIBpCTLHo37gE5yVT8WgcCrtvsAKvQpWpYaGHaTEv92oQqeJkaFlycNlNTwzahLOuagAHkJGmrT7lasmRJ1k/BOslU/FoHAv6j2wFWoUvVsNDCtJmW+rUJVfAyNSy4OG2mpoZtQlnWNQEDyEnSVp+AVV0qNdD3gFQ3lI8bN+4ko/u1bgn4j263ZHsXVw17Z92NV1K/bqj2NqYa9ta7G6+mht1QFbMTAQNIJ1pHrZ0xY0Zas2ZN/U3o1acgrF+/vv4m9OnTp6fVq1cHIvvVbgn4j263ZHsXVw17Z92NV1K/bqj2NqYa9ta7G6+mht1QFbMTAQNIJ1pHrd2/f39aunRp6u/vT7t27UojRoxIfX19ae7cuWno0KGByH61WwLVRx9XQ2M1JD7iEY/o1suI20UBNewibg9Cq18PkLv8EmrYZeAehFfDHiB7iRMKGEA0CAECBAgQIECAAAECPRMwgPSM2gsRIECAAAECBAgQIGAA0QMECBAgQIAAAQIECPRMwADSM2ovRIAAAQIECBAgQICAAUQPECBAgAABAgQIECDQMwEDSM+ovRABAgQIECBAgAABAgYQPTDoBb71rW/VH418ww031O/1/PPPrz8++clPfnLL937zzTenq666Kl1//fVpx44d6YwzzkhPetKT0jve8Y703Oc+t+XvW9CewIEDB9KyZcvqj7TeuXNnGjlyZP2R1nPmzElDhgxpGSRS45bBLWhL4GRreNddd6W//du/TZ/97GfTli1b0k9/+tN01llnpRe84AVp3rx56bd/+7fben2LYgInW7/jveqznvWs9G//9m/pNa95Tf3fT4/eCOSo4T/8wz+kD3zgA+mWW25J9957bxozZky65JJL6o+u9yCQU8AAklNTrOIEvve976Xf+73fSw95yEPqb62vHu9///vTz3/+8/SNb3wjPeYxjzlhzhdddFH6yle+kl760pemiRMnpjvvvLP+sslqw3vFFVfU/2H2iAvMnDkzrVq1qv5Sz0mTJqUNGzbUztXPV65cecIXiNY4nr0IlcDJ1rD6d2ncuHHpGc94RpoyZUp62MMeljZu3FgPo9UgUv3vD37wgyF3WeBk63d0WtUwWcX65S9/aQDpcs2ODh+t4Zvf/Oa0fPny9PKXvzxVQ2Q1gFT/ff3N3/zNtHjx4h6/Gy832AUMIIO9wqf4+6sGiHXr1qWtW7fWf1WvHtVf2J/whCek5z//+emTn/zkCYWqb7d/2tOelk477bRD6+6+++70lKc8Jd1xxx1p7969vnQy2GPVX73Hjx9fD4jV//M7+Jg9e3ZasWJF2rx5cxo7duyArxKtcTB9v55SfXJxsjX8yU9+km6//fZ6CDn88Td/8zfp9a9/ffqrv/qrVG2MPLonEKnf4VlVf9h53OMel/78z/+8PiV2AtK9mh0dOVrDz33uc/Wp49///d+nV73qVb1L3CudsgIGkFO29IP/jVenFQ996EPTK1/5yvSRj3zkiDf82te+Nn384x9P1ebngQ98YMcYb3nLW9J73/ve9MMf/vDQYNNxEL9QC8yfP7/+69r27dvTqFGjDqlUl7yNHj26vgxn0aJFx9XqZo2Vp32BSA0HepX/+q//Sr/1W7+VXve616UPfehD7SdjZccCuep36aWXpi9+8Yv1CXH1RxsDSMelOOlfiNbwmc98Zqouh6wuO65OPqr/tj7oQQ866Xz8IoFWAgaQVkKeb6zAjTfeWF/OU13aM2PGjCPeR/Wz6ri6WnPuued2/B6rvxB96lOfqi/lqu4L8Th5geqym+qUY8+ePccEGT58eJowYUJ9inW8RzdrfPLv6NT7zUgNB9L6zne+kx7/+Ment7/97WnJkiWnHmoP33GO+lWXyp1zzjnp2muvrU+X73e/+xlAGlLDatiohv3qkuIzzzyzvkz5Zz/7Wfqd3/md+hSy+gPR/e9//x6+Gy91KggYQE6FKp+i7/HTn/50qi7PqW5ufeELX3iEQvWzCy+8sB4iqvs7OnlUl3NVl2BVx9XVa3jEBKrLq4YNG1Zf63/0o7rvZt++ffUlPsd7dKvGsXd06v12pIYDaf3pn/5pfQPzN7/5zfryLo/uCUTr96tf/ar+Q84jHvGI9JnPfKZO1ADSvXodL3Kkhps2bar/0FMNH1Ut//Iv/zKNGDEifexjH0tXX321QbK3pTxlXs0AcsqUutlvtNqEfv/732/rTTzgAQ9Ij3rUo9Lf/d3fpWoT8/nPfz4973nPO+J3v/CFL9Q3vFZrpk6d2lbcatEvfvGLdN5556Xdu3fXf7WvXscjJlB9ykp10lHdeH70ozrBqu6z2bZt23FfpBs1jr2bU/O3IzU8ntiVV15Zf+pOde9HdQ+IR3cFovVbvXp1fd/HrbfeeugySgNId2t2dPRIDatPiKw+BKJ6VJ9edvB/r/7v5zznOfWnQP7Hf/xHeuITn9jbN+XVBrWAAWRQl3fwvLnbbrvtiPsDTvTOqk/v+PKXv1yfTuQ8AaluPq+Glptuuqm+JKh6HY+4QOQvd7lrHH83p2aESA2PFvunf/qn+t/bP/qjP6r/+jp06NBTE7WH7zpSv+rDOKobz9/0pjeld77znYeyNoD0sIAp1R/UcbInydXpc/VhK9WnzlX33h3++PCHP1x/OqFPfextPU+FVzOAnApVHgTvsdr8Vx+H286j+sjd6lrknPcH3HPPPelFL3pR/ZegalNUXX7lkUcgcv15zhrneTenZpRIDQ8Xq04mq3/PqpOvf/7nf64//tOj+wKR+lWfXld9clL1R5/TTz/9ULLVR5xXl7lefvnl9aU9vs+lu3WM1LC6/666fO7pT396+trXvnZEotUf26p7eqoPAqk+EMSDQC4BA0guSXGKE2j1CUnVFy5VX3rW6lOw9u/fX/9F9pprrkkf/ehH60/V8sgnUP0/teom4258Cla7Nc73bk7NSJEaHhSr/sBQbXSqv+Red911Lf+9PDWlu/OuI/V78YtffOi+j4Gyq75k9K1vfWt3khe1FojUsPr9wz+m/nDS6vt43vCGN6QPfvCD9ZfDehDIJWAAySUpTpEC1Q3m1T0g3/72t+ub6qrHwe8Bqf5idPhN5NV9HdU9HtW1tAc/8aO6Ie/Vr351/ZG91XXp1X+IPfIKVPfSVDdADvQ9INUNktV3RBy8D6j6tJbqr3UHH53UOG/moh0UiNbw61//enruc59bX2ZZDSLVp+949E4gUr/qFPJHP/rRMcm+7GUvq+8l+LM/+7P639/HPvaxvXtDp+ArRWpYcb3tbW9LS5curf/QdvCEv/pm9eo0srpEq7oPr7pEy4NALgEDSC5JcYoUqD7Ks7ocq/o+kOr/EVaP6iMGq5OP6pvQq2uXDz4uvvji+vtCqmtgD/6HtroJ9n3ve199v8fx/vpzwQUX1DdQe8QEqo9JXrNmTX2t8eTJk1P1BZDVN6FXNyJXN7hWj4P3AVXf4VJdl3zw0UmNY1n67RMJnGwNf/CDH9QD6H//93+nd7/73cf8+1T9+1X9e+bRXYGTrd9AWbkHpLv1Ol70SA2rj92t7gOp/hBXfQnsIx/5yPSJT3wiffWrX/VR2L0v5SnxigaQU6LMp/abvOWWW9LcuXPrTW31OP/889N73vOeY755+XgDyLOf/ewT3nvypS99KVVrPGIC1WVu1V/fquP+Xbt21adV1cBX1e3gTcgDDSDVK7db41iWfvtEAidbw+regd///d8fMPTBD5Wg312Bk62fAaS7dekkerSG1fBRfYN9df9VdTXA2WefXX+4QPX9IB4EcgsYQHKLikeAAAECBAgQIECAwIACBhDNQYAAAQIECBAgQIBAzwQMID2j9kIECBAgQIAAAQIECBhA9AABAgQIECBAgAABAj0TMID0jNoLESBAgAABAgQIECBgANEDBAgQIECAAAECBAj0TMAA0jNqL0SAAAECBAgQIECAgAFEDxAgQIAAAQIECBAg0DMBA0jPqL0QAQIECBAgQIAAAQIGED1AgAABAgQIECBAgEDPBAwgPaP2QgQIECBAgAABAgQIGED0AAECBAgQIECAAAECPRMwgPSM2gsRIECAAAECBAgQIGAA0QMECBAgQIAAAQIECPRMwADSM2ovRIAAAQIECBAgQICAAUQPECBAgAABAgQIECDQMwEDSM+ovRABAgQIECBAgAABAgYQPUCAAAECBAgQIECAQM8EDCA9o/ZCBAgQIECAAAECBAgYQPQAAQIECBAgQIAAAQI9EzCA9IzaCxEgQIAAAQIECBAgYADRAwQIECBAgAABAgQI9EzAANIzai9EgAABAgQIECBAgIABRA8QIECAAAECBAgQINAzAQNIz6i9EAECBAgQIECAAAECBhA9QIAAAQIECBAgQIBAzwQMID2j9kIECBAgQIAAAQIECBhA9AABAgQIECBAgAABAj0TMID0jNoLESBAgAABAgQIECBgANEDBAgQIECAAAECBAj0TMAA0jNqL0SAAAECBAgQIECAgAFEDxAgQIAAAQIECBAg0DMBA0jPqL0QAQIECBAgQIAAAQIGED1AgAABAgQIECBAgEDPBAwgPaP2QgQIECBAgAABAgQIGED0AAECBAgQIECAAAECPRMwgPSM2gsRIECAAAECBAgQIGAA0QMECBAgQIAAAQIECPRMwADSM2ovRIAAAQIECBAgQICAAUQPECBAgAABAgQIECDQM4H/D/ivmH24FCDyAAAAAElFTkSuQmCC\" width=\"640\">"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "data = WmonI1.w #shape(A x 8000 = running time)\n",
    "number_of_frames = WmonI1.w.shape[1]\n",
    "\n",
    "def update_hist(num, data):\n",
    "    plt.cla()\n",
    "    plt.hist(data[num])\n",
    "\n",
    "fig = plt.figure()\n",
    "hist = plt.hist(data[0])\n",
    "\n",
    "animation = FuncAnimation(fig, update_hist, number_of_frames, fargs=(data, ))\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 86,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(2181, 8000)"
      ]
     },
     "execution_count": 86,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "shape(WmonE1.w)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 87,
   "metadata": {},
   "outputs": [],
   "source": [
    "# data = WmonE1.w #shape(A x 8000 = running time)\n",
    "# number_of_frames = 8000\n",
    "\n",
    "# def update_hist(num, data):\n",
    "#     plt.cla()\n",
    "#     plt.hist(data[num])\n",
    "\n",
    "# fig = plt.figure()\n",
    "# hist = plt.hist(data[0])\n",
    "\n",
    "# animation = FuncAnimation(fig, update_hist, number_of_frames, fargs=(data, ))\n",
    "# plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 88,
   "metadata": {},
   "outputs": [],
   "source": [
    "# data = WmonI2.w #shape(A x 8000 = running time)\n",
    "# number_of_frames = 8000\n",
    "\n",
    "# def update_hist(num, data):\n",
    "#     plt.cla()\n",
    "#     plt.hist(data[num])\n",
    "\n",
    "# fig = plt.figure()\n",
    "# hist = plt.hist(data[0])\n",
    "\n",
    "# animation = FuncAnimation(fig, update_hist, number_of_frames, fargs=(data, ))\n",
    "# plt.show()"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "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.10.5"
  },
  "vscode": {
   "interpreter": {
    "hash": "01f4427c4855e802c456af4062c13258c3be347db629675187ef9df5c74fa0bc"
   }
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}