Lab9.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Решить однородное уравнение теплопроводности с граничными условиями Дирихле \n",
    "\\begin{equation} \\label{c9eq1}\n",
    "\t\\begin{cases}\n",
    "\t\t\\displaystyle \\frac{\\partial u}{\\partial t} = k \\frac{\\partial^2 u}{\\partial x^2} \\\\\n",
    "\t\t\\displaystyle u(x,0) = e^{-(x-5)^4} + 0.01x \\\\\n",
    "\t\tu(0,t) = u(0,0) \\\\  \n",
    "\t\tu(a,t) = u(a,0)\n",
    "\t\\end{cases}\n",
    "\\end{equation}\n",
    "на отрезке $x \\in [0; a]$ при $t \\in [0; T]$ Выбрать $a = 20$ и $T = 20$, а коэффициент теплопроводности $k = 2$. Шаг по пространству $h = 0.01$, шаг по времени $\\tau = 0.05$. Расчет проводить с помощью комплексной схемы Розенброка. После каждого временного слоя выводить решение на текущем слое для получения анимационной картинки. Чтобы избежать постоянного изменения масштаба графика рекомендуется после команды $\\texttt{plot}$ вставить команду $\\texttt{axis([0 a 0 1])}$, не забыв далее поставить команду $\\texttt{pause(1e-6)}$.\n",
    "Тот же расчет необходимо повторить с граничными условиями Неймана\n",
    "\\begin{equation} \\label{c9eq2}\n",
    "\tu_x(0,t) = u_x(a,t) = 0.\n",
    "\\end{equation}\n",
    "\n",
    "Счет будет идти быстрее, если использовать аппарат разреженных матриц MATLAB (смотри документацию к функции $\\texttt{spdiags}$).\n",
    "\n",
    "Граничные условия надо включить в оператор пространственного дифференцирования $\\Lambda_x$, видоизменив в нем первую и последнюю строки. Поскольку в данной задаче коэффициент теплопроводности постоянен, то оператор $\\Lambda_x$ можно вычислить один раз до начала расчета.\n",
    "\n",
    "Для аппроксимации граничных условий Неймана со 2 порядком следует использовать метод фиктивных точек.\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 25,
   "metadata": {
    "code_folding": []
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Wall time: 7.77 s\n"
     ]
    }
   ],
   "source": [
    "%%time\n",
    "import numpy as np\n",
    "import scipy.sparse.linalg\n",
    "from numpy import exp\n",
    "from scipy import sparse\n",
    "\n",
    "def matprint(mat, fmt=\"g\"):\n",
    "    col_maxes = [max([len((\"{:\"+fmt+\"}\").format(x)) for x in col]) for col in mat.T]\n",
    "    for x in mat:\n",
    "        for i, y in enumerate(x):\n",
    "            print((\"{:\"+str(col_maxes[i])+fmt+\"}\").format(y), end=\"  \")\n",
    "        print(\"\")\n",
    "\n",
    "def mult(a, b):\n",
    "    res = np.matrix(np.zeros(a.shape))\n",
    "    for i in range(a.shape[0]):\n",
    "        res[i] = a[i]*b[i]\n",
    "    return res\n",
    "\n",
    "def ux0(x):\n",
    "    return [0.01*t + exp(-(t-5)**4) for t in x]\n",
    "\n",
    "#global a, T, k, h, tau, x, t, sigma\n",
    "\n",
    "a = 20\n",
    "sigma = (1+1j)/2\n",
    "T = 20\n",
    "k = 2\n",
    "h = 0.01\n",
    "tau = 0.05\n",
    "\n",
    "x = np.arange(0, a+h, h)\n",
    "t = np.arange(0, T+tau, tau)\n",
    "u = np.zeros((len(t), len(x)))\n",
    "\n",
    "for i,ui in enumerate(ux0(x)):\n",
    "    u[0][i] = ui \n",
    "    \n",
    "for i in range(len(t)):\n",
    "    u[i][0] = u[0][0]\n",
    "\n",
    "for i in range(len(t)):\n",
    "    u[i][-1] = u[0][-1]\n",
    "    \n",
    "diff_op = np.zeros((len(x),len(x)),dtype= np.float64)\n",
    "for n in range(1, len(x)-1):\n",
    "    diff_op[n][n-1] = 1\n",
    "    diff_op[n][n] = -2\n",
    "    diff_op[n][n+1] = 1\n",
    "diff_op = k/h**2 * diff_op\n",
    "A = np.eye(len(x),dtype= np.float64) - sigma*tau*diff_op;\n",
    "A = sparse.csc_matrix(A)\n",
    "\n",
    "for k in range(len(t)-1):\n",
    "    B = np.matrix(diff_op,dtype= np.float64)*np.matrix(u[k],dtype= np.float64).T\n",
    "    w = scipy.sparse.linalg.spsolve(A, B)\n",
    "    u[k+1] = np.matrix(u[k]) + tau*np.real(w)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 26,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "Traceback (most recent call last):\n",
      "  File \"D:\\Programs\\Anaconda3\\lib\\site-packages\\matplotlib\\cbook\\__init__.py\", line 388, in process\n",
      "    proxy(*args, **kwargs)\n",
      "  File \"D:\\Programs\\Anaconda3\\lib\\site-packages\\matplotlib\\cbook\\__init__.py\", line 228, in __call__\n",
      "    return mtd(*args, **kwargs)\n",
      "  File \"D:\\Programs\\Anaconda3\\lib\\site-packages\\matplotlib\\animation.py\", line 1499, in _stop\n",
      "    self.event_source.remove_callback(self._loop_delay)\n",
      "AttributeError: 'NoneType' object has no attribute 'remove_callback'\n"
     ]
    },
    {
     "data": {
      "application/javascript": [
       "/* Put everything inside the global mpl namespace */\n",
       "window.mpl = {};\n",
       "\n",
       "\n",
       "mpl.get_websocket_type = function() {\n",
       "    if (typeof(WebSocket) !== 'undefined') {\n",
       "        return WebSocket;\n",
       "    } else if (typeof(MozWebSocket) !== 'undefined') {\n",
       "        return MozWebSocket;\n",
       "    } else {\n",
       "        alert('Your browser does not have WebSocket support.' +\n",
       "              'Please try Chrome, Safari or Firefox ≥ 6. ' +\n",
       "              'Firefox 4 and 5 are also supported but you ' +\n",
       "              'have to enable WebSockets in about:config.');\n",
       "    };\n",
       "}\n",
       "\n",
       "mpl.figure = function(figure_id, websocket, ondownload, parent_element) {\n",
       "    this.id = figure_id;\n",
       "\n",
       "    this.ws = websocket;\n",
       "\n",
       "    this.supports_binary = (this.ws.binaryType != undefined);\n",
       "\n",
       "    if (!this.supports_binary) {\n",
       "        var warnings = document.getElementById(\"mpl-warnings\");\n",
       "        if (warnings) {\n",
       "            warnings.style.display = 'block';\n",
       "            warnings.textContent = (\n",
       "                \"This browser does not support binary websocket messages. \" +\n",
       "                    \"Performance may be slow.\");\n",
       "        }\n",
       "    }\n",
       "\n",
       "    this.imageObj = new Image();\n",
       "\n",
       "    this.context = undefined;\n",
       "    this.message = undefined;\n",
       "    this.canvas = undefined;\n",
       "    this.rubberband_canvas = undefined;\n",
       "    this.rubberband_context = undefined;\n",
       "    this.format_dropdown = undefined;\n",
       "\n",
       "    this.image_mode = 'full';\n",
       "\n",
       "    this.root = $('<div/>');\n",
       "    this._root_extra_style(this.root)\n",
       "    this.root.attr('style', 'display: inline-block');\n",
       "\n",
       "    $(parent_element).append(this.root);\n",
       "\n",
       "    this._init_header(this);\n",
       "    this._init_canvas(this);\n",
       "    this._init_toolbar(this);\n",
       "\n",
       "    var fig = this;\n",
       "\n",
       "    this.waiting = false;\n",
       "\n",
       "    this.ws.onopen =  function () {\n",
       "            fig.send_message(\"supports_binary\", {value: fig.supports_binary});\n",
       "            fig.send_message(\"send_image_mode\", {});\n",
       "            if (mpl.ratio != 1) {\n",
       "                fig.send_message(\"set_dpi_ratio\", {'dpi_ratio': mpl.ratio});\n",
       "            }\n",
       "            fig.send_message(\"refresh\", {});\n",
       "        }\n",
       "\n",
       "    this.imageObj.onload = function() {\n",
       "            if (fig.image_mode == 'full') {\n",
       "                // Full images could contain transparency (where diff images\n",
       "                // almost always do), so we need to clear the canvas so that\n",
       "                // there is no ghosting.\n",
       "                fig.context.clearRect(0, 0, fig.canvas.width, fig.canvas.height);\n",
       "            }\n",
       "            fig.context.drawImage(fig.imageObj, 0, 0);\n",
       "        };\n",
       "\n",
       "    this.imageObj.onunload = function() {\n",
       "        fig.ws.close();\n",
       "    }\n",
       "\n",
       "    this.ws.onmessage = this._make_on_message_function(this);\n",
       "\n",
       "    this.ondownload = ondownload;\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._init_header = function() {\n",
       "    var titlebar = $(\n",
       "        '<div class=\"ui-dialog-titlebar ui-widget-header ui-corner-all ' +\n",
       "        'ui-helper-clearfix\"/>');\n",
       "    var titletext = $(\n",
       "        '<div class=\"ui-dialog-title\" style=\"width: 100%; ' +\n",
       "        'text-align: center; padding: 3px;\"/>');\n",
       "    titlebar.append(titletext)\n",
       "    this.root.append(titlebar);\n",
       "    this.header = titletext[0];\n",
       "}\n",
       "\n",
       "\n",
       "\n",
       "mpl.figure.prototype._canvas_extra_style = function(canvas_div) {\n",
       "\n",
       "}\n",
       "\n",
       "\n",
       "mpl.figure.prototype._root_extra_style = function(canvas_div) {\n",
       "\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._init_canvas = function() {\n",
       "    var fig = this;\n",
       "\n",
       "    var canvas_div = $('<div/>');\n",
       "\n",
       "    canvas_div.attr('style', 'position: relative; clear: both; outline: 0');\n",
       "\n",
       "    function canvas_keyboard_event(event) {\n",
       "        return fig.key_event(event, event['data']);\n",
       "    }\n",
       "\n",
       "    canvas_div.keydown('key_press', canvas_keyboard_event);\n",
       "    canvas_div.keyup('key_release', canvas_keyboard_event);\n",
       "    this.canvas_div = canvas_div\n",
       "    this._canvas_extra_style(canvas_div)\n",
       "    this.root.append(canvas_div);\n",
       "\n",
       "    var canvas = $('<canvas/>');\n",
       "    canvas.addClass('mpl-canvas');\n",
       "    canvas.attr('style', \"left: 0; top: 0; z-index: 0; outline: 0\")\n",
       "\n",
       "    this.canvas = canvas[0];\n",
       "    this.context = canvas[0].getContext(\"2d\");\n",
       "\n",
       "    var backingStore = this.context.backingStorePixelRatio ||\n",
       "\tthis.context.webkitBackingStorePixelRatio ||\n",
       "\tthis.context.mozBackingStorePixelRatio ||\n",
       "\tthis.context.msBackingStorePixelRatio ||\n",
       "\tthis.context.oBackingStorePixelRatio ||\n",
       "\tthis.context.backingStorePixelRatio || 1;\n",
       "\n",
       "    mpl.ratio = (window.devicePixelRatio || 1) / backingStore;\n",
       "\n",
       "    var rubberband = $('<canvas/>');\n",
       "    rubberband.attr('style', \"position: absolute; left: 0; top: 0; z-index: 1;\")\n",
       "\n",
       "    var pass_mouse_events = true;\n",
       "\n",
       "    canvas_div.resizable({\n",
       "        start: function(event, ui) {\n",
       "            pass_mouse_events = false;\n",
       "        },\n",
       "        resize: function(event, ui) {\n",
       "            fig.request_resize(ui.size.width, ui.size.height);\n",
       "        },\n",
       "        stop: function(event, ui) {\n",
       "            pass_mouse_events = true;\n",
       "            fig.request_resize(ui.size.width, ui.size.height);\n",
       "        },\n",
       "    });\n",
       "\n",
       "    function mouse_event_fn(event) {\n",
       "        if (pass_mouse_events)\n",
       "            return fig.mouse_event(event, event['data']);\n",
       "    }\n",
       "\n",
       "    rubberband.mousedown('button_press', mouse_event_fn);\n",
       "    rubberband.mouseup('button_release', mouse_event_fn);\n",
       "    // Throttle sequential mouse events to 1 every 20ms.\n",
       "    rubberband.mousemove('motion_notify', mouse_event_fn);\n",
       "\n",
       "    rubberband.mouseenter('figure_enter', mouse_event_fn);\n",
       "    rubberband.mouseleave('figure_leave', mouse_event_fn);\n",
       "\n",
       "    canvas_div.on(\"wheel\", function (event) {\n",
       "        event = event.originalEvent;\n",
       "        event['data'] = 'scroll'\n",
       "        if (event.deltaY < 0) {\n",
       "            event.step = 1;\n",
       "        } else {\n",
       "            event.step = -1;\n",
       "        }\n",
       "        mouse_event_fn(event);\n",
       "    });\n",
       "\n",
       "    canvas_div.append(canvas);\n",
       "    canvas_div.append(rubberband);\n",
       "\n",
       "    this.rubberband = rubberband;\n",
       "    this.rubberband_canvas = rubberband[0];\n",
       "    this.rubberband_context = rubberband[0].getContext(\"2d\");\n",
       "    this.rubberband_context.strokeStyle = \"#000000\";\n",
       "\n",
       "    this._resize_canvas = function(width, height) {\n",
       "        // Keep the size of the canvas, canvas container, and rubber band\n",
       "        // canvas in synch.\n",
       "        canvas_div.css('width', width)\n",
       "        canvas_div.css('height', height)\n",
       "\n",
       "        canvas.attr('width', width * mpl.ratio);\n",
       "        canvas.attr('height', height * mpl.ratio);\n",
       "        canvas.attr('style', 'width: ' + width + 'px; height: ' + height + 'px;');\n",
       "\n",
       "        rubberband.attr('width', width);\n",
       "        rubberband.attr('height', height);\n",
       "    }\n",
       "\n",
       "    // Set the figure to an initial 600x600px, this will subsequently be updated\n",
       "    // upon first draw.\n",
       "    this._resize_canvas(600, 600);\n",
       "\n",
       "    // Disable right mouse context menu.\n",
       "    $(this.rubberband_canvas).bind(\"contextmenu\",function(e){\n",
       "        return false;\n",
       "    });\n",
       "\n",
       "    function set_focus () {\n",
       "        canvas.focus();\n",
       "        canvas_div.focus();\n",
       "    }\n",
       "\n",
       "    window.setTimeout(set_focus, 100);\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._init_toolbar = function() {\n",
       "    var fig = this;\n",
       "\n",
       "    var nav_element = $('<div/>')\n",
       "    nav_element.attr('style', 'width: 100%');\n",
       "    this.root.append(nav_element);\n",
       "\n",
       "    // Define a callback function for later on.\n",
       "    function toolbar_event(event) {\n",
       "        return fig.toolbar_button_onclick(event['data']);\n",
       "    }\n",
       "    function toolbar_mouse_event(event) {\n",
       "        return fig.toolbar_button_onmouseover(event['data']);\n",
       "    }\n",
       "\n",
       "    for(var toolbar_ind in mpl.toolbar_items) {\n",
       "        var name = mpl.toolbar_items[toolbar_ind][0];\n",
       "        var tooltip = mpl.toolbar_items[toolbar_ind][1];\n",
       "        var image = mpl.toolbar_items[toolbar_ind][2];\n",
       "        var method_name = mpl.toolbar_items[toolbar_ind][3];\n",
       "\n",
       "        if (!name) {\n",
       "            // put a spacer in here.\n",
       "            continue;\n",
       "        }\n",
       "        var button = $('<button/>');\n",
       "        button.addClass('ui-button ui-widget ui-state-default ui-corner-all ' +\n",
       "                        'ui-button-icon-only');\n",
       "        button.attr('role', 'button');\n",
       "        button.attr('aria-disabled', 'false');\n",
       "        button.click(method_name, toolbar_event);\n",
       "        button.mouseover(tooltip, toolbar_mouse_event);\n",
       "\n",
       "        var icon_img = $('<span/>');\n",
       "        icon_img.addClass('ui-button-icon-primary ui-icon');\n",
       "        icon_img.addClass(image);\n",
       "        icon_img.addClass('ui-corner-all');\n",
       "\n",
       "        var tooltip_span = $('<span/>');\n",
       "        tooltip_span.addClass('ui-button-text');\n",
       "        tooltip_span.html(tooltip);\n",
       "\n",
       "        button.append(icon_img);\n",
       "        button.append(tooltip_span);\n",
       "\n",
       "        nav_element.append(button);\n",
       "    }\n",
       "\n",
       "    var fmt_picker_span = $('<span/>');\n",
       "\n",
       "    var fmt_picker = $('<select/>');\n",
       "    fmt_picker.addClass('mpl-toolbar-option ui-widget ui-widget-content');\n",
       "    fmt_picker_span.append(fmt_picker);\n",
       "    nav_element.append(fmt_picker_span);\n",
       "    this.format_dropdown = fmt_picker[0];\n",
       "\n",
       "    for (var ind in mpl.extensions) {\n",
       "        var fmt = mpl.extensions[ind];\n",
       "        var option = $(\n",
       "            '<option/>', {selected: fmt === mpl.default_extension}).html(fmt);\n",
       "        fmt_picker.append(option)\n",
       "    }\n",
       "\n",
       "    // Add hover states to the ui-buttons\n",
       "    $( \".ui-button\" ).hover(\n",
       "        function() { $(this).addClass(\"ui-state-hover\");},\n",
       "        function() { $(this).removeClass(\"ui-state-hover\");}\n",
       "    );\n",
       "\n",
       "    var status_bar = $('<span class=\"mpl-message\"/>');\n",
       "    nav_element.append(status_bar);\n",
       "    this.message = status_bar[0];\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.request_resize = function(x_pixels, y_pixels) {\n",
       "    // Request matplotlib to resize the figure. Matplotlib will then trigger a resize in the client,\n",
       "    // which will in turn request a refresh of the image.\n",
       "    this.send_message('resize', {'width': x_pixels, 'height': y_pixels});\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.send_message = function(type, properties) {\n",
       "    properties['type'] = type;\n",
       "    properties['figure_id'] = this.id;\n",
       "    this.ws.send(JSON.stringify(properties));\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.send_draw_message = function() {\n",
       "    if (!this.waiting) {\n",
       "        this.waiting = true;\n",
       "        this.ws.send(JSON.stringify({type: \"draw\", figure_id: this.id}));\n",
       "    }\n",
       "}\n",
       "\n",
       "\n",
       "mpl.figure.prototype.handle_save = function(fig, msg) {\n",
       "    var format_dropdown = fig.format_dropdown;\n",
       "    var format = format_dropdown.options[format_dropdown.selectedIndex].value;\n",
       "    fig.ondownload(fig, format);\n",
       "}\n",
       "\n",
       "\n",
       "mpl.figure.prototype.handle_resize = function(fig, msg) {\n",
       "    var size = msg['size'];\n",
       "    if (size[0] != fig.canvas.width || size[1] != fig.canvas.height) {\n",
       "        fig._resize_canvas(size[0], size[1]);\n",
       "        fig.send_message(\"refresh\", {});\n",
       "    };\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_rubberband = function(fig, msg) {\n",
       "    var x0 = msg['x0'] / mpl.ratio;\n",
       "    var y0 = (fig.canvas.height - msg['y0']) / mpl.ratio;\n",
       "    var x1 = msg['x1'] / mpl.ratio;\n",
       "    var y1 = (fig.canvas.height - msg['y1']) / mpl.ratio;\n",
       "    x0 = Math.floor(x0) + 0.5;\n",
       "    y0 = Math.floor(y0) + 0.5;\n",
       "    x1 = Math.floor(x1) + 0.5;\n",
       "    y1 = Math.floor(y1) + 0.5;\n",
       "    var min_x = Math.min(x0, x1);\n",
       "    var min_y = Math.min(y0, y1);\n",
       "    var width = Math.abs(x1 - x0);\n",
       "    var height = Math.abs(y1 - y0);\n",
       "\n",
       "    fig.rubberband_context.clearRect(\n",
       "        0, 0, fig.canvas.width, fig.canvas.height);\n",
       "\n",
       "    fig.rubberband_context.strokeRect(min_x, min_y, width, height);\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_figure_label = function(fig, msg) {\n",
       "    // Updates the figure title.\n",
       "    fig.header.textContent = msg['label'];\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_cursor = function(fig, msg) {\n",
       "    var cursor = msg['cursor'];\n",
       "    switch(cursor)\n",
       "    {\n",
       "    case 0:\n",
       "        cursor = 'pointer';\n",
       "        break;\n",
       "    case 1:\n",
       "        cursor = 'default';\n",
       "        break;\n",
       "    case 2:\n",
       "        cursor = 'crosshair';\n",
       "        break;\n",
       "    case 3:\n",
       "        cursor = 'move';\n",
       "        break;\n",
       "    }\n",
       "    fig.rubberband_canvas.style.cursor = cursor;\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_message = function(fig, msg) {\n",
       "    fig.message.textContent = msg['message'];\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_draw = function(fig, msg) {\n",
       "    // Request the server to send over a new figure.\n",
       "    fig.send_draw_message();\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_image_mode = function(fig, msg) {\n",
       "    fig.image_mode = msg['mode'];\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.updated_canvas_event = function() {\n",
       "    // Called whenever the canvas gets updated.\n",
       "    this.send_message(\"ack\", {});\n",
       "}\n",
       "\n",
       "// A function to construct a web socket function for onmessage handling.\n",
       "// Called in the figure constructor.\n",
       "mpl.figure.prototype._make_on_message_function = function(fig) {\n",
       "    return function socket_on_message(evt) {\n",
       "        if (evt.data instanceof Blob) {\n",
       "            /* FIXME: We get \"Resource interpreted as Image but\n",
       "             * transferred with MIME type text/plain:\" errors on\n",
       "             * Chrome.  But how to set the MIME type?  It doesn't seem\n",
       "             * to be part of the websocket stream */\n",
       "            evt.data.type = \"image/png\";\n",
       "\n",
       "            /* Free the memory for the previous frames */\n",
       "            if (fig.imageObj.src) {\n",
       "                (window.URL || window.webkitURL).revokeObjectURL(\n",
       "                    fig.imageObj.src);\n",
       "            }\n",
       "\n",
       "            fig.imageObj.src = (window.URL || window.webkitURL).createObjectURL(\n",
       "                evt.data);\n",
       "            fig.updated_canvas_event();\n",
       "            fig.waiting = false;\n",
       "            return;\n",
       "        }\n",
       "        else if (typeof evt.data === 'string' && evt.data.slice(0, 21) == \"data:image/png;base64\") {\n",
       "            fig.imageObj.src = evt.data;\n",
       "            fig.updated_canvas_event();\n",
       "            fig.waiting = false;\n",
       "            return;\n",
       "        }\n",
       "\n",
       "        var msg = JSON.parse(evt.data);\n",
       "        var msg_type = msg['type'];\n",
       "\n",
       "        // Call the  \"handle_{type}\" callback, which takes\n",
       "        // the figure and JSON message as its only arguments.\n",
       "        try {\n",
       "            var callback = fig[\"handle_\" + msg_type];\n",
       "        } catch (e) {\n",
       "            console.log(\"No handler for the '\" + msg_type + \"' message type: \", msg);\n",
       "            return;\n",
       "        }\n",
       "\n",
       "        if (callback) {\n",
       "            try {\n",
       "                // console.log(\"Handling '\" + msg_type + \"' message: \", msg);\n",
       "                callback(fig, msg);\n",
       "            } catch (e) {\n",
       "                console.log(\"Exception inside the 'handler_\" + msg_type + \"' callback:\", e, e.stack, msg);\n",
       "            }\n",
       "        }\n",
       "    };\n",
       "}\n",
       "\n",
       "// from http://stackoverflow.com/questions/1114465/getting-mouse-location-in-canvas\n",
       "mpl.findpos = function(e) {\n",
       "    //this section is from http://www.quirksmode.org/js/events_properties.html\n",
       "    var targ;\n",
       "    if (!e)\n",
       "        e = window.event;\n",
       "    if (e.target)\n",
       "        targ = e.target;\n",
       "    else if (e.srcElement)\n",
       "        targ = e.srcElement;\n",
       "    if (targ.nodeType == 3) // defeat Safari bug\n",
       "        targ = targ.parentNode;\n",
       "\n",
       "    // jQuery normalizes the pageX and pageY\n",
       "    // pageX,Y are the mouse positions relative to the document\n",
       "    // offset() returns the position of the element relative to the document\n",
       "    var x = e.pageX - $(targ).offset().left;\n",
       "    var y = e.pageY - $(targ).offset().top;\n",
       "\n",
       "    return {\"x\": x, \"y\": y};\n",
       "};\n",
       "\n",
       "/*\n",
       " * return a copy of an object with only non-object keys\n",
       " * we need this to avoid circular references\n",
       " * http://stackoverflow.com/a/24161582/3208463\n",
       " */\n",
       "function simpleKeys (original) {\n",
       "  return Object.keys(original).reduce(function (obj, key) {\n",
       "    if (typeof original[key] !== 'object')\n",
       "        obj[key] = original[key]\n",
       "    return obj;\n",
       "  }, {});\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.mouse_event = function(event, name) {\n",
       "    var canvas_pos = mpl.findpos(event)\n",
       "\n",
       "    if (name === 'button_press')\n",
       "    {\n",
       "        this.canvas.focus();\n",
       "        this.canvas_div.focus();\n",
       "    }\n",
       "\n",
       "    var x = canvas_pos.x * mpl.ratio;\n",
       "    var y = canvas_pos.y * mpl.ratio;\n",
       "\n",
       "    this.send_message(name, {x: x, y: y, button: event.button,\n",
       "                             step: event.step,\n",
       "                             guiEvent: simpleKeys(event)});\n",
       "\n",
       "    /* This prevents the web browser from automatically changing to\n",
       "     * the text insertion cursor when the button is pressed.  We want\n",
       "     * to control all of the cursor setting manually through the\n",
       "     * 'cursor' event from matplotlib */\n",
       "    event.preventDefault();\n",
       "    return false;\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._key_event_extra = function(event, name) {\n",
       "    // Handle any extra behaviour associated with a key event\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.key_event = function(event, name) {\n",
       "\n",
       "    // Prevent repeat events\n",
       "    if (name == 'key_press')\n",
       "    {\n",
       "        if (event.which === this._key)\n",
       "            return;\n",
       "        else\n",
       "            this._key = event.which;\n",
       "    }\n",
       "    if (name == 'key_release')\n",
       "        this._key = null;\n",
       "\n",
       "    var value = '';\n",
       "    if (event.ctrlKey && event.which != 17)\n",
       "        value += \"ctrl+\";\n",
       "    if (event.altKey && event.which != 18)\n",
       "        value += \"alt+\";\n",
       "    if (event.shiftKey && event.which != 16)\n",
       "        value += \"shift+\";\n",
       "\n",
       "    value += 'k';\n",
       "    value += event.which.toString();\n",
       "\n",
       "    this._key_event_extra(event, name);\n",
       "\n",
       "    this.send_message(name, {key: value,\n",
       "                             guiEvent: simpleKeys(event)});\n",
       "    return false;\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.toolbar_button_onclick = function(name) {\n",
       "    if (name == 'download') {\n",
       "        this.handle_save(this, null);\n",
       "    } else {\n",
       "        this.send_message(\"toolbar_button\", {name: name});\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.toolbar_button_onmouseover = function(tooltip) {\n",
       "    this.message.textContent = tooltip;\n",
       "};\n",
       "mpl.toolbar_items = [[\"Home\", \"Reset original view\", \"fa fa-home icon-home\", \"home\"], [\"Back\", \"Back to  previous view\", \"fa fa-arrow-left icon-arrow-left\", \"back\"], [\"Forward\", \"Forward to next view\", \"fa fa-arrow-right icon-arrow-right\", \"forward\"], [\"\", \"\", \"\", \"\"], [\"Pan\", \"Pan axes with left mouse, zoom with right\", \"fa fa-arrows icon-move\", \"pan\"], [\"Zoom\", \"Zoom to rectangle\", \"fa fa-square-o icon-check-empty\", \"zoom\"], [\"\", \"\", \"\", \"\"], [\"Download\", \"Download plot\", \"fa fa-floppy-o icon-save\", \"download\"]];\n",
       "\n",
       "mpl.extensions = [\"eps\", \"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 overridden (by mpl) onmessage function.\n",
       "        ws.onmessage(msg['content']['data'])\n",
       "    });\n",
       "    return ws;\n",
       "}\n",
       "\n",
       "mpl.mpl_figure_comm = function(comm, msg) {\n",
       "    // This is the function which gets called when the mpl process\n",
       "    // starts-up an IPython Comm through the \"matplotlib\" channel.\n",
       "\n",
       "    var id = msg.content.data.id;\n",
       "    // Get hold of the div created by the display call when the Comm\n",
       "    // socket was opened in Python.\n",
       "    var element = $(\"#\" + id);\n",
       "    var ws_proxy = comm_websocket_adapter(comm)\n",
       "\n",
       "    function ondownload(figure, format) {\n",
       "        window.open(figure.imageObj.src);\n",
       "    }\n",
       "\n",
       "    var fig = new mpl.figure(id, ws_proxy,\n",
       "                           ondownload,\n",
       "                           element.get(0));\n",
       "\n",
       "    // Call onopen now - mpl needs it, as it is assuming we've passed it a real\n",
       "    // web socket which is closed, not our websocket->open comm proxy.\n",
       "    ws_proxy.onopen();\n",
       "\n",
       "    fig.parent_element = element.get(0);\n",
       "    fig.cell_info = mpl.find_output_cell(\"<div id='\" + id + \"'></div>\");\n",
       "    if (!fig.cell_info) {\n",
       "        console.error(\"Failed to find cell for figure\", id, fig);\n",
       "        return;\n",
       "    }\n",
       "\n",
       "    var output_index = fig.cell_info[2]\n",
       "    var cell = fig.cell_info[0];\n",
       "\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_close = function(fig, msg) {\n",
       "    var width = fig.canvas.width/mpl.ratio\n",
       "    fig.root.unbind('remove')\n",
       "\n",
       "    // Update the output cell to use the data from the current canvas.\n",
       "    fig.push_to_output();\n",
       "    var dataURL = fig.canvas.toDataURL();\n",
       "    // Re-enable the keyboard manager in IPython - without this line, in FF,\n",
       "    // the notebook keyboard shortcuts fail.\n",
       "    IPython.keyboard_manager.enable()\n",
       "    $(fig.parent_element).html('<img src=\"' + dataURL + '\" width=\"' + width + '\">');\n",
       "    fig.close_ws(fig, msg);\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.close_ws = function(fig, msg){\n",
       "    fig.send_message('closing', msg);\n",
       "    // fig.ws.close()\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.push_to_output = function(remove_interactive) {\n",
       "    // Turn the data on the canvas into data in the output cell.\n",
       "    var width = this.canvas.width/mpl.ratio\n",
       "    var dataURL = this.canvas.toDataURL();\n",
       "    this.cell_info[1]['text/html'] = '<img src=\"' + dataURL + '\" width=\"' + width + '\">';\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.updated_canvas_event = function() {\n",
       "    // Tell IPython that the notebook contents must change.\n",
       "    IPython.notebook.set_dirty(true);\n",
       "    this.send_message(\"ack\", {});\n",
       "    var fig = this;\n",
       "    // Wait a second, then push the new image to the DOM so\n",
       "    // that it is saved nicely (might be nice to debounce this).\n",
       "    setTimeout(function () { fig.push_to_output() }, 1000);\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._init_toolbar = function() {\n",
       "    var fig = this;\n",
       "\n",
       "    var nav_element = $('<div/>')\n",
       "    nav_element.attr('style', 'width: 100%');\n",
       "    this.root.append(nav_element);\n",
       "\n",
       "    // Define a callback function for later on.\n",
       "    function toolbar_event(event) {\n",
       "        return fig.toolbar_button_onclick(event['data']);\n",
       "    }\n",
       "    function toolbar_mouse_event(event) {\n",
       "        return fig.toolbar_button_onmouseover(event['data']);\n",
       "    }\n",
       "\n",
       "    for(var toolbar_ind in mpl.toolbar_items){\n",
       "        var name = mpl.toolbar_items[toolbar_ind][0];\n",
       "        var tooltip = mpl.toolbar_items[toolbar_ind][1];\n",
       "        var image = mpl.toolbar_items[toolbar_ind][2];\n",
       "        var method_name = mpl.toolbar_items[toolbar_ind][3];\n",
       "\n",
       "        if (!name) { continue; };\n",
       "\n",
       "        var button = $('<button class=\"btn btn-default\" href=\"#\" title=\"' + name + '\"><i class=\"fa ' + image + ' fa-lg\"></i></button>');\n",
       "        button.click(method_name, toolbar_event);\n",
       "        button.mouseover(tooltip, toolbar_mouse_event);\n",
       "        nav_element.append(button);\n",
       "    }\n",
       "\n",
       "    // Add the status bar.\n",
       "    var status_bar = $('<span class=\"mpl-message\" style=\"text-align:right; float: right;\"/>');\n",
       "    nav_element.append(status_bar);\n",
       "    this.message = status_bar[0];\n",
       "\n",
       "    // Add the close button to the window.\n",
       "    var buttongrp = $('<div class=\"btn-group inline pull-right\"></div>');\n",
       "    var button = $('<button class=\"btn btn-mini btn-primary\" href=\"#\" title=\"Stop Interaction\"><i class=\"fa fa-power-off icon-remove icon-large\"></i></button>');\n",
       "    button.click(function (evt) { fig.handle_close(fig, {}); } );\n",
       "    button.mouseover('Stop Interaction', toolbar_mouse_event);\n",
       "    buttongrp.append(button);\n",
       "    var titlebar = this.root.find($('.ui-dialog-titlebar'));\n",
       "    titlebar.prepend(buttongrp);\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._root_extra_style = function(el){\n",
       "    var fig = this\n",
       "    el.on(\"remove\", function(){\n",
       "\tfig.close_ws(fig, {});\n",
       "    });\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._canvas_extra_style = function(el){\n",
       "    // this is important to make the div 'focusable\n",
       "    el.attr('tabindex', 0)\n",
       "    // reach out to IPython and tell the keyboard manager to turn it's self\n",
       "    // off when our div gets focus\n",
       "\n",
       "    // location in version 3\n",
       "    if (IPython.notebook.keyboard_manager) {\n",
       "        IPython.notebook.keyboard_manager.register_events(el);\n",
       "    }\n",
       "    else {\n",
       "        // location in version 2\n",
       "        IPython.keyboard_manager.register_events(el);\n",
       "    }\n",
       "\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._key_event_extra = function(event, name) {\n",
       "    var manager = IPython.notebook.keyboard_manager;\n",
       "    if (!manager)\n",
       "        manager = IPython.keyboard_manager;\n",
       "\n",
       "    // Check for shift+enter\n",
       "    if (event.shiftKey && event.which == 13) {\n",
       "        this.canvas_div.blur();\n",
       "        event.shiftKey = false;\n",
       "        // Send a \"J\" for go to next cell\n",
       "        event.which = 74;\n",
       "        event.keyCode = 74;\n",
       "        manager.command_mode();\n",
       "        manager.handle_keydown(event);\n",
       "    }\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_save = function(fig, msg) {\n",
       "    fig.ondownload(fig, null);\n",
       "}\n",
       "\n",
       "\n",
       "mpl.find_output_cell = function(html_output) {\n",
       "    // Return the cell and output element which can be found *uniquely* in the notebook.\n",
       "    // Note - this is a bit hacky, but it is done because the \"notebook_saving.Notebook\"\n",
       "    // IPython event is triggered only after the cells have been serialised, which for\n",
       "    // our purposes (turning an active figure into a static one), is too late.\n",
       "    var cells = IPython.notebook.get_cells();\n",
       "    var ncells = cells.length;\n",
       "    for (var i=0; i<ncells; i++) {\n",
       "        var cell = cells[i];\n",
       "        if (cell.cell_type === 'code'){\n",
       "            for (var j=0; j<cell.output_area.outputs.length; j++) {\n",
       "                var data = cell.output_area.outputs[j];\n",
       "                if (data.data) {\n",
       "                    // IPython >= 3 moved mimebundle to data attribute of output\n",
       "                    data = data.data;\n",
       "                }\n",
       "                if (data['text/html'] == html_output) {\n",
       "                    return [cell, data, j];\n",
       "                }\n",
       "            }\n",
       "        }\n",
       "    }\n",
       "}\n",
       "\n",
       "// Register the function which deals with the matplotlib target/channel.\n",
       "// The kernel may be null if the page has been refreshed.\n",
       "if (IPython.notebook.kernel != null) {\n",
       "    IPython.notebook.kernel.comm_manager.register_target('matplotlib', mpl.mpl_figure_comm);\n",
       "}\n"
      ],
      "text/plain": [
       "<IPython.core.display.Javascript object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/html": [
       "<img src=\"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAoAAAAHgCAYAAAA10dzkAAAgAElEQVR4Xu2dCZSU5ZW/b4vNpuwiILKpuLCMgom4RAb/Dhk1JwNqJkbOuHCMmpjJKGYmxkQNxKCJnBFMJi7oGJCIcVyIOkIiM4kQAWeMDkbZDCpGFhe2Vmxptv6f+0F1mu6qrq/7dlN16z6fp09j9/dW3ff3/u69T31bl1VXV1cLGwqgAAqgAAqgAAqgQBgFygDAMGvNRFEABVAABVAABVAgUQAAxAgogAIogAIogAIoEEwBADDYgjNdFEABFEABFEABFAAA8QAKoAAKoAAKoAAKBFMAAAy24EwXBVAABVAABVAABQBAPIACKIACKIACKIACwRQAAIMtONNFARRAARRAARRAAQAQD6AACqAACqAACqBAMAUAwGALznRRAAVQAAVQAAVQAADEAyiAAiiAAiiAAigQTAEAMNiCM10UQAEUQAEUQAEUAADxAAqgAAqgAAqgAAoEUwAADLbgTBcFUAAFUAAFUAAFAEA8gAIogAIogAIogALBFAAAgy0400UBFEABFEABFEABABAPoAAKoAAKoAAKoEAwBQDAYAvOdFEABVAABVAABVAAAMQDKIACKIACKIACKBBMAQAw2IIzXRRAARRAARRAARQAAPEACqAACqAACqAACgRTAAAMtuBMFwVQAAVQAAVQAAUAQDyAAiiAAiiAAiiAAsEUAACDLTjTRQEUQAEUQAEUQAEAEA+gAAqgAAqgAAqgQDAFAMBgC850UQAFUAAFUAAFUAAAxAMogAIogAIogAIoEEwBADDYgjNdFEABFEABFEABFAAA8QAKoAAKoAAKoAAKBFMAAAy24EwXBVAABVAABVAABQBAPIACKIACKIACKIACwRQAAIMtONNFARRAARRAARRAAQAQD6AACqAACqAACqBAMAUAwGALznRRAAVQAAVQAAVQAADEAyiAAiiAAiiAAigQTAEAMNiCM10UQAEUQAEUQAEUAADxAAqgAAqgAAqgAAoEUwAADLbgTBcFUAAFUAAFUAAFAEA8gAIogAIogAIogALBFAAAgy0400UBFEABFEABFEABABAPoAAKoAAKoAAKoEAwBQDAYAvOdFEABVAABVAABVAAAMQDKIACKIACKIACKBBMAQAw2IIzXRRAARRAARRAARQAAPEACqAACqAACqAACgRTAAAMtuBMFwVQAAVQAAVQAAUAQDyAAiiAAiiAAiiAAsEUAACDLTjTRQEUQAEUQAEUQAEAEA+gAAqgAAqgAAqgQDAFAMBgC850UQAFUAAFUAAFUAAAxAMogAIogAIogAIoEEwBADDYgjNdFEABFEABFEABFAAA8QAKoAAKoAAKoAAKBFMAAAy24EwXBVAABVAABVAABQBAPIACKIACKIACKIACwRQAAIMtONNFARRAARRAARRAAQAQD6AACqAACqAACqBAMAUAQMOC79mzR9avXy8dOnSQsrIywysxFAVQAAVQAAVQ4EApUF1dLR9//LEcccQRctBBBx2oty2q9wEADcuxdu1a6dOnj+EVGIoCKIACKIACKFAoBd5991058sgjC/X2BX3fkgDA22+/XZ588klZuXKltGvXTk4//XT58Y9/LMcdd1yD4j7xxBNy8803y5tvvilHH320TJ48Wc4///zUC1JRUSGdO3eWt99+W7p27Zp6nPcdd+7cKc8995x8/vOfl/Lycu/TSR0/82a9U5vF8Y74HJ87tm/q0Ddv3iwDBgyQrVu3SqdOnVKPK6UdSwIAzznnHPnKV74in/3sZ2XXrl3yve99T1577TVZvny5HHLIIVnXa8mSJXLmmWfKrbfemkDfnDlz5JZbbpEXXnhBRowYkWqNP/roo8Q4GzdulG7duqUaUwo7aYOYO3eunHfeeeEAkHnHAn7Wm/UuhZrd0Byi1vNNmzbJYYcdJnogp2PHjqW+zFnnVxIAWHdmH374oRx++OGyYMECGTlyZNaJX3TRRaIAN2/evJrfK0h26dJFHnnkkVRmAAABwFRGcb5T1AbBvMlv56mbKvyoPgcARUoSAFevXi0DBw5MjgIOGTIkaxL07dtXJkyYkHxltqlTp8q0adPknXfeSZU4ACANIpVRnO8UtUEwb/LbeeqmCj+qzwHAEgRAvbNnzJgxsmXLFvn973+fMwFat24tM2bMkHHjxtXsM3v2bBk/frxUVVVlHac/r/07BUC9CWTDhg3hTgHPnz9fRo8eHe4UMPOOdUqQ9Wa9U1GU450UACP6XAGwV69enAJ27N16oX/jG9+QZ599NrmWr6E7exQAZ86cKRdffHHNazz88MNyxRVXyPbt27NKMnHiRJk0aVK93yk4tm/fvpRkZC4ogAIogAIoULIKVFZWJgeAuAawRJb4m9/8pvzqV7+ShQsXJnf3NLQ15RQwRwD3Khr1EyPz5ohviZTKBqeBz/F5BJ9zBLBETgHraV+FP72T9/nnn0+u/8u36U0g+hBIvcsvs5177rnJY124CaRh9aJeM8K8uSYsX10phd/jc3xeCj7ONweuASwRALzmmmtET8M+9dRT+z37Tx/Ros8F1O3SSy+V3r17iz4zULfFixcndwjrs//0mkEde9NNN/EYmHxZs+8III/HiHVtFOvNeqcoDa53AXxjgS8AWCIAmOvPsP385z+Xyy+/PClKo0aNkv79+yc3fmS2xx9/PIG+t956q+ZB0BdccEHqIsZdwLEKBg2C9U5dHBzviM/xuWP7pg4dACwRAEy94s28IwBIoWxmSxXlywEE+LwojdnMQeHzWD4HAAFAUwkBAGMVDBoE620qGE4G43N87sSqpjABQADQZCAAkEJpMpCTwQABPndiVVOY+DyWzwFAANBUMADAWAWDBsF6mwqGk8H4HJ87saopTAAQADQZCACkUJoM5GQwQIDPnVjVFCY+j+VzABAANBUMADBWwaBBsN6mguFkMD7H506sagoTAAQATQYCACmUJgM5GQwQ4HMnVjWFic9j+RwABABNBQMAjFUwaBCst6lgOBmMz/G5E6uawgQAAUCTgQBACqXJQE4GAwT43IlVTWHi81g+BwABQFPBAABjFQwaBOttKhhOBuNzfO7EqqYwAUAA0GQgAJBCaTKQk8EAAT53YlVTmPg8ls8BQADQVDAAwFgFgwbBepsKhpPB+ByfO7GqKUwAEAA0GQgApFCaDORkMECAz51Y1RQmPo/lcwAQADQVDAAwVsGgQbDepoLhZDA+x+dOrGoKEwAEAE0GAgAplCYDORkMEOBzJ1Y1hYnPY/kcAAQATQUDAIxVMGgQrLepYDgZjM/xuROrmsIEAAFAk4EAQAqlyUBOBgME+NyJVU1h4vNYPgcAAUBTwQAAYxUMGgTrbSoYTgbjc3zuxKqmMAFAANBkIACQQmkykJPBAAE+d2JVU5j4PJbPAUAA0FQwAMBYBYMGwXqbCoaTwfgcnzuxqilMABAANBkIAKRQmgzkZDBAgM+dWNUUJj6P5XMAEAA0FQwAMFbBoEGw3qaC4WQwPsfnTqxqChMABABNBgIAKZQmAzkZDBDgcydWNYWJz2P5HAAEAE0FAwCMVTBoEKy3qWA4GYzP8bkTq5rCBAABQJOBAEAKpclATgYDBPjciVVNYeLzWD4HAAFAU8EAAGMVDBoE620qGE4G43N87sSqpjABQADQZCAAkEJpMpCTwQABPndiVVOY+DyWzwFAANBUMADAWAWDBsF6mwqGk8H4HJ87saopTAAQADQZCACkUJoM5GQwQIDPnVjVFCY+j+VzABAANBUMADBWwaBBsN6mguFkMD7H506sagoTAAQATQYCACmUJgM5GQwQ4HMnVjWFic9j+RwABABNBQMAjFUwaBCst6lgOBmMz/G5E6uawgQAAUCTgQBACqXJQE4GAwT43IlVTWHi81g+BwABQFPBAABjFQwaBOttKhhOBuNzfO7EqqYwAUAA0GQgAJBCaTKQk8EAAT53YlVTmPg8ls8BQADQVDAAwFgFgwbBepsKhpPB+ByfO7GqKUwAEAA0GQgApFCaDORkMECAz51Y1RQmPo/lcwAQADQVDAAwVsGgQbDepoLhZDA+x+dOrGoKEwAEAE0GAgAplCYDORkMEOBzJ1Y1hYnPY/kcAAQATQUDAIxVMGgQrLepYDgZjM/xuROrmsIEAAFAk4EAQAqlyUBOBgME+NyJVU1h4vNYPgcAAUBTwQAAYxUMGgTrbSoYTgbjc3zuxKqmMAFAANBkIACQQmkykJPBAAE+d2JVU5j4PJbPAUAA0FQwAMBYBYMGwXqbCoaTwfgcnzuxqilMABAANBkIAKRQmgzkZDBAgM+dWNUUJj6P5XMAEAA0FQwAMFbBoEGw3qaC4WQwPsfnTqxqChMABABNBgIAKZQmAzkZDBDgcydWNYWJz2P5HAAEAE0FAwCMVTBoEKy3qWA4GYzP8bkTq5rCBAABQJOBAEAKpclATgYDBPjciVVNYeLzWD4HAAFAU8EAAGMVDBoE620qGE4G43N87sSqpjABQADQZCAAkEJpMpCTwQABPndiVVOY+DyWzwHAEgLAhQsXypQpU+Tll1+WDRs2yJw5c2Ts2LE5C8Lzzz8vZ511Vr3fr1ixQo4//vhUhQQAjFUwaBCsd6rC4HwnfI7PnVs4VfgAYAkB4Lx582TRokUyfPhwufDCC1MD4KpVq6Rjx441hunevbu0atUqlYEAQAplKqM43wkgwOfOLZwqfHwey+cAYAkBYO0MLysrSw2AW7Zskc6dO6cqEHV3AgBjFQwaBOvdpELhbBA+x+fOLNukcAFAAFD69+8v27dvl0GDBslNN92U9bRwxl1VVVWiX5lNAbBPnz7JKedu3bo1yYQeB2mDmD9/vowePVrKy8s9TqFJMTNv1rtJxnE2CJ/jc2eWbVK4CoC9evWSioqK/c4CNunFnA4qq66urnYae86w0xwB1FO/et3gySefnEDdrFmz5N577xW9NnDkyJFZX3vixIkyadKker+bPXu2tG/fvtRkZD4ogAIogAIoUJIKVFZWyrhx4wDAUlvdNACYbc5f/OIXRcc+/fTTWSXhCOBeWThCwBGCUqsZ2eaDz/E5Pi9dBTgCGPgUcDZbT548WX7xi1+I3gmcZuMaQK6VSeMT7/twTRg+9+7hNPHj81g+5xpAAHC/uvClL31JNm/eLL/97W/T1AsBAGMVDBoE652qMDjfCZ/jc+cWThU+AFhCALht2zZZvXp1svDDhg2TO++8M7mho2vXrtK3b1+58cYbZd26dfLQQw8l+0ybNi25AWTw4MGyY8eO5Mjfj370I3niiSfkggsuSGUgAJBCmcoozncCCPC5cwunCh+fx/I5AFhCAJjrwc6XXXaZzJgxQy6//HJZs2ZNcpOHbnfccYdMnz49gcJ27dolIKiQeN5556UqFroTABirYNAgWO/UxcHxjvgcnzu2b+rQAcASAsDUq96MOwKAFMpmtFPRvhRAgM+L1pzNGBg+j+VzABAANJUPADBWwaBBsN6mguFkMD7H506sagoTAAQATQYCACmUJgM5GQwQ4HMnVjWFic9j+RwABABNBQMAjFUwaBCst6lgOBmMz/G5E6uawgQAAUCTgQBACqXJQE4GAwT43IlVTWHi81g+BwABQFPBAABjFQwaBOttKhhOBuNzfO7EqqYwAUAA0GQgAJBCaTKQk8EAAT53YlVTmPg8ls8BQADQVDAAwFgFgwbBepsKhpPB+ByfO7GqKUwAEAA0GQgApFCaDORkMECAz51Y1RQmPo/lcwAQADQVDAAwVsGgQbDepoLhZDA+x+dOrGoKEwAEAE0GAgAplCYDORkMEOBzJ1Y1hYnPY/kcAAQATQUDAIxVMGgQrLepYDgZjM/xuROrmsIEAAFAk4EAQAqlyUBOBgME+NyJVU1h4vNYPgcAAUBTwQAAYxUMGgTrbSoYTgbjc3zuxKqmMAFAANBkIACQQmkykJPBAAE+d2JVU5j4PJbPAUAA0FQwAMBYBYMGwXqbCoaTwfgcnzuxqilMABAANBkIAKRQmgzkZDBAgM+dWNUUJj6P5XMAEAA0FQwAMFbBoEGw3qaC4WQwPsfnTqxqChMABABNBgIAKZQmAzkZDBDgcydWNYWJz2P5HAAEAE0FAwCMVTBoEKy3qWA4GYzP8bkTq5rCBAABQJOBAEAKpclATgYDBPjciVVNYeLzWD4HAAFAU8EAAGMVDBoE620qGE4G43N87sSqpjABQADQZCAAkEJpMpCTwQABPndiVVOY+DyWzwFAANBUMADAWAWDBsF6mwqGk8H4HJ87saopTAAQADQZCACkUJoM5GQwQIDPnVjVFCY+j+VzABAANBUMADBWwaBBsN6mguFkMD7H506sagoTAAQATQYCACmUJgM5GQwQ4HMnVjWFic9j+RwABABNBQMAjFUwaBCst6lgOBmMz/G5E6uawgQAAUCTgQBACqXJQE4GAwT43IlVTWHi81g+BwABQFPBAABjFQwaBOttKhhOBuNzfO7EqqYwAUAA0GQgAJBCaTKQk8EAAT53YlVTmPg8ls8BQADQVDAAwFgFgwbBepsKhpPB+ByfO7GqKUwAEAA0GQgApFCaDORkMECAz51Y1RQmPo/lcwAQADQVDAAwVsGgQbDepoLhZDA+x+dOrGoKEwAEAE0GAgAplCYDORkMEOBzJ1Y1hYnPY/kcAAQATQUDAIxVMGgQrLepYDgZjM/xuROrmsIEAAFAk4EAQAqlyUBOBgME+NyJVU1h4vNYPgcAAUBTwQAAYxUMGgTrbSoYTgbjc3zuxKqmMAFAANBkIACQQmkykJPBAAE+d2JVU5j4PJbPAUAA0FQwAMBYBYMGwXqbCoaTwfgcnzuxqilMABAANBkIAKRQmgzkZDBAgM+dWNUUJj6P5XMAEAA0FQwAMFbBoEGw3qaC4WQwPsfnTqxqChMABABNBgIAKZQmAzkZDBDgcydWNYWJz2P5HAAEAE0FAwCMVTBoEKy3qWA4GYzP8bkTq5rCBAABQJOBAEAKpclATgYDBPjciVVNYeLzWD4HAAFAU8EAAGMVDBoE620qGE4G43N87sSqpjABQADQZCAAkEJpMpCTwQABPndiVVOY+DyWzwFAANBUMADAWAWDBsF6mwqGk8H4HJ87saopTAAQADQZCACkUJoM5GQwQIDPnVjVFCY+j+VzABAANBUMADBWwaBBsN6mguFkMD7H506sagoTACwhAFy4cKFMmTJFXn75ZdmwYYPMmTNHxo4d26BBFixYINdff70sW7ZMjjjiCPn2t78tX/va11KbCgCkUKY2i+MdAQJ87ti+qUPH57F8DgCWEADOmzdPFi1aJMOHD5cLL7wwLwC+/fbbMmTIELnyyivl6quvTsZec8018sgjjyTj02wAYKyCQYNgvdPUBe/74HN87t3DaeIHAEsIAGsveFlZWV4AvOGGG+Tpp5+WFStW1AzVo3+vvvqqLFmyJI1/BACkUKYyivOdAAJ87tzCqcLH57F8DgAGBsCRI0fKsGHD5K677qopDnra+Mtf/rJUVlZKeXl53qIBAMYqGDQI1jtvUSiBHfA5Pi8BG+edAgAYGACPPfZYufzyy+W73/1ujVEWL14sZ5xxhqxfv1569epVz0BVVVWiX5lNAbBPnz7JNYfdunXLa7hS2UEbxPz582X06NGpQJl5+1aA9cbnvh2cLnp8HsvnCoDa5ysqKqRjx47pTFJie5VVV1dXl9icJM0pYAXA8ePHy4033lgzfb0O8HOf+1wCdD179qwny8SJE2XSpEn1fj579mxp3759qcnIfFAABVAABVCgJBXQM33jxo0DAEttddMAYFNOAXMEcK9T+KQc65My6816l1qPyDYffB7L5xwBDHwKWG8CeeaZZ2T58uU1teDrX/+6LF26lJtA8lR7rhHiGqEoQDB37lw57zzWm/UuXQWi1nOuASwhANy2bZusXr06yVK9uePOO++Us846S7p27Sp9+/ZNTvWuW7dOHnrooWSfzGNg9BEw+igYvfNX7wLmMTD5C13UgsG8AaH82eF/D3yOz/27OP8MAMASAsDnn38+Ab6622WXXSYzZsxIbvhYs2aN6H6ZTR8EPWHChJoHQetRQR4EnT9xaBA0iPwu8b8HPsfn/l2cfwZRfQ4AlhAA5rd58+/BY2BoEM3vquJ7xagNgnmT38WXjc0fUVSfA4AAoCmbAEAahMlATgZHbRDMm/x2kqKmMKP6HAAEAE2JAwDSIEwGcjI4aoNg3uS3kxQ1hRnV5wAgAGhKHACQBmEykJPBURsE8ya/naSoKcyoPgcAAUBT4gCANAiTgZwMjtogmDf57SRFTWFG9TkACACaEgcApEGYDORkcNQGwbzJbycpagozqs8BQADQlDgAIA3CZCAng6M2COZNfjtJUVOYUX0OAAKApsQBAGkQJgM5GRy1QTBv8ttJiprCjOpzABAANCUOAEiDMBnIyeCoDYJ5k99OUtQUZlSfA4AAoClxAEAahMlATgZHbRDMm/x2kqKmMKP6HAAEAE2JAwDSIEwGcjI4aoNg3uS3kxQ1hRnV5wAgAGhKHACQBmEykJPBURsE8ya/naSoKcyoPgcAAUBT4gCANAiTgZwMjtogmDf57SRFTWFG9TkACACaEgcApEGYDORkcNQGwbzJbycpagozqs8BQADQlDgAIA3CZCAng6M2COZNfjtJUVOYUX0OAAKApsQBAGkQJgM5GRy1QTBv8ttJiprCjOpzABAANCUOAEiDMBnIyeCoDYJ5k99OUtQUZlSfA4AAoClxAEAahMlATgZHbRDMm/x2kqKmMKP6HAAEAE2JAwDSIEwGcjI4aoNg3uS3kxQ1hRnV5wAgAGhKHACQBmEykJPBURsE8ya/naSoKcyoPgcAAUBT4gCANAiTgZwMjtogmDf57SRFTWFG9TkACACaEgcApEGYDORkcNQGwbzJbycpagozqs8BQADQlDgAIA3CZCAng6M2COZNfjtJUVOYUX0OAAKApsQBAGkQJgM5GRy1QTBv8ttJiprCjOpzABAANCUOAEiDMBnIyeCoDYJ5k99OUtQUZlSfA4AAoClxAEAahMlATgZHbRDMm/x2kqKmMKP6HAAEAE2JAwDSIEwGcjI4aoNg3uS3kxQ1hRnV5wAgAGhKHACQBmEykJPBURsE8ya/naSoKcyoPgcAAUBT4gCANAiTgZwMjtogmDf57SRFTWFG9TkACACaEgcApEGYDORkcNQGwbzJbycpagozqs8BQADQlDgAIA3CZCAng6M2COZNfjtJUVOYUX0OAAKApsQBAGkQJgM5GRy1QTBv8ttJiprCjOpzABAANCUOAEiDMBnIyeCoDYJ5k99OUtQUZlSfA4AAoClxAEAahMlATgZHbRDMm/x2kqKmMKP6HAAEAE2JAwDSIEwGcjI4aoNg3uS3kxQ1hRnV5wAgAGhKHACQBmEykJPBURsE8ya/naSoKcyoPgcAAUBT4gCANAiTgZwMjtogmDf57SRFTWFG9TkACACaEgcApEGYDORkcNQGwbzJbycpagozqs8BQADQlDgAIA3CZCAng6M2COZNfjtJUVOYUX0OAAKApsQBAGkQJgM5GRy1QTBv8ttJiprCjOpzABAANCUOAEiDMBnIyeCoDYJ5k99OUtQUZlSfA4AAoClxAEAahMlATgZHbRDMm/x2kqKmMKP6HAAEAE2JAwDSIEwGcjI4aoNg3uS3kxQ1hRnV5wAgAGhKHACQBmEykJPBURsE8ya/naSoKcyoPgcAAUBT4gCANAiTgZwMjtogmDf57SRFTWFG9TkACACaEgcApEGYDORkcNQGwbzJbycpagozqs8BQADQlDgAIA3CZCAng6M2COZNfjtJUVOYUX0OAAKApsQBAGkQJgM5GRy1QTBv8ttJiprCjOpzABAANCUOAEiDMBnIyeCoDYJ5k99OUtQUZlSfA4AAoClxAEAahMlATgZHbRDMm/x2kqKmMKP6HAAEAE2JAwDSIEwGcjI4aoNg3uS3kxQ1hRnV5wBgiQHg3XffLVOmTJENGzbI4MGDZdq0aXLmmWdmTY4ZM2bI+PHj6/3u008/lbZt26ZKKACQBpHKKM53itogmDf57Tx1U4Uf1ecAYAkB4KOPPiqXXHKJKASeccYZct9998kDDzwgy5cvl759+9ZLBAXAa6+9VlatWrXf73r27JkqaXQnAJAGkdosjneM2iCYN/ntOG1Thx7V5wBgCQHgiBEjZPjw4XLPPffUGP+EE06QsWPHyu23354VAK+77jrZunVr6kSpuyMASINosnkcDYzaIJg3+e0oTZscalSfA4AlAoA7duyQ9u3by2OPPSbnn39+TSLoEb6lS5fKggULsgLgV7/6Vendu7fs3r1bTjrpJLn11ltl2LBhOROpqqpK9CuzKQD26dMnOeXcrVu3Jiegt4FaMObPny+jR4+W8vJyb+E3OV7mzXo32TyOBuJzfO7Irk0OVQGwV69eUlFRIR07dmzy63geWFZdXV3teQIa+/r16xOQW7RokZx++uk107nttttk5syZ9U7z6g4vvviirF69WoYOHZqcyr3rrrtk7ty58uqrr8rAgQOzSjJx4kSZNGlSvd/Nnj07AVA2FEABFEABFECB4legsrJSxo0bBwAW/1I1HGEGABcvXiynnXZazc6TJ0+WWbNmycqVK/NOcc+ePckp5JEjR8pPfvKTrPtzBHCvLBwh4AhB3oQqgR3wOT4vARvnnUJUn3MEMPAp4GxZceWVV8ratWtl3rx5eZNGd+AaQK4RSmUU5ztFvUaIeZPfzlM3VfhRfc41gCUCgOpyvQnk5JNPTu4CzmyDBg2SMWPGZL0JpG5m6JnwU045JTkl/OCDD6ZKHACQBpHKKM53itogmDf57Tx1U4Uf1ecAYAkBYOYxMPfee29yGnj69Oly//33y7Jly6Rfv35y6aWXJtcJZu4I1mv5Tj311OR6PwU5Pe2rp4v1OkIFwY01I+AAACAASURBVDQbAEiDSOMT7/tEbRDMm/z2nrtp4o/qcwCwhABQja5H/+64447krtwhQ4bI1KlTk2v6dBs1apT0799f9Pl/uk2YMEGefPJJee+996RTp07J3b96k0ftawjzJQ8ASIPI55FS+H3UBsG8ye9SyN98c4jqcwCwxAAwn9Gb+/cAIA2iuT1VjK8XtUEwb/K7GPOxuWOK6nMAEAA05RIASIMwGcjJ4KgNgnmT305S1BRmVJ8DgACgKXEAQBqEyUBOBkdtEMyb/HaSoqYwo/ocAAQATYkDANIgTAZyMjhqg2De5LeTFDWFGdXnACAAaEocAJAGYTKQk8FRGwTzJr+dpKgpzKg+BwABQFPiAIA0CJOBnAyO2iCYN/ntJEVNYUb1OQAIAJoSBwCkQZgM5GRw1AbBvMlvJylqCjOqzwFAANCUOAAgDcJkICeDozYI5k1+O0lRU5hRfQ4AAoCmxAEAaRAmAzkZHLVBMG/y20mKmsKM6nMAEAA0JQ4ASIMwGcjJ4KgNgnmT305S1BRmVJ8DgACgKXEAQBqEyUBOBkdtEMyb/HaSoqYwo/ocAAQATYkDANIgTAZyMjhqg2De5LeTFDWFGdXnACAAaEocAJAGYTKQk8FRGwTzJr+dpKgpzKg+BwABQFPiAIA0CJOBnAyO2iCYN/ntJEVNYUb1OQAIAJoSBwCkQZgM5GRw1AbBvMlvJylqCjOqzwFAANCUOAAgDcJkICeDozYI5k1+O0lRU5hRfQ4AAoCmxAEAaRAmAzkZHLVBMG/y20mKmsKM6nMAEAA0JQ4ASIMwGcjJ4KgNgnmT305S1BRmVJ8DgACgKXEAQBqEyUBOBkdtEMyb/HaSoqYwo/ocAAQATYkDANIgTAZyMjhqg2De5LeTFDWFGdXnACAAaEocAJAGYTKQk8FRGwTzJr+dpKgpzKg+BwABQFPiAIA0CJOBnAyO2iCYN/ntJEVNYUb1OQAIAJoSBwCkQZgM5GRw1AbBvMlvJylqCjOqzwFAANCUOAAgDcJkICeDozYI5k1+O0lRU5hRfQ4AAoCmxAEAaRAmAzkZHLVBMG/y20mKmsKM6nMAEAA0JQ4ASIMwGcjJ4KgNgnmT305S1BRmVJ8DgACgKXEAQBqEyUBOBkdtEMyb/HaSoqYwo/ocAAQATYkDANIgTAZyMjhqg2De5LeTFDWFGdXnACAAaEocAJAGYTKQk8FRGwTzJr+dpKgpzKg+BwABQFPiAIA0CJOBnAyO2iCYN/ntJEVNYZa6z3fv2S3bd22v9/X+h+/L2UPOloqKCunYsaNJQ6+Dy6qrq6u9Bl/ouAFAGkShPXgg3r/UG0QuDZk3+X0g8qvQ79GSPle82LF7h1TtrkoArGpXVfJv/Z78f45/6+8+3fVpVnDL+fOd2fffuWdndom3i8iPBAAstAG9vj8ASIMohHe1qGqBbFfeLnn7XXt2yac7P5UObTq0SDh1G4S+v76nFtbk++693/P9TD+JV+t/1dWyp3pPzb/1Z8n/Vye/3e/fdffV/XQrKyuT5L+yMjmo7KB6/05+tm+fbP/ONkb3b1XWSlod1Cr5vmf3Hln0+0Vy1qizpF2bdvv9Tvc5+KCD6/1Mx+nreN5aEgiKWRfP89bcUtDSnNTvyb937/135meZ/6/7s0+qPpH/feV/5fjBx8su2ZUb1PKAW13I0xql71VMm+Zs24PbJl+td7aW9d9bDwAW0wJ5igUABAAb61ctiK1btU4gRyHEuunrVO6slG07tiVfn+z4ZO/3nXu/1/5Ztk/OmaKd7RRJ7Z9VVlVKdVl10kwyEGaNvVTH67rWBslcsJgNIHNBZQZKk9feB6g5v5e1SryVb79c+8gekRUrVsiJQ06U1uWt975WLTDO9rqZfWqDeQLq+yA9873uz+rCvO7XmJ/V/iCR+XeuDxOZn+cas3PXTnnppZdk+MnDpVWrVjUfRnR/BSz9gLO7eu/35N91fpbv/7OO2fd6mbF1YS0t0OkHJw+b1r42rdpIm4PbJBCW+bd+T/7/4DbS7uB2NZCWgbXM9wZ/V559XGaMvrbmV2bjGkCuATTlDAAYGwAVtg5pfUgCXfq9MZse2VIgq9heIVu3b5WKqn3fs/1/1da9+2yvkI+qPvoL2O38JGlEWuB0vG5aSPXIoO5/oDctrvpVflB58j35d6u//DtzdKz20Tdt+LmO1uU6uqc/z3ZkMNvRxWzNPtuRxgwcaIPXZpz5vn3Hdjmo1UE1zT/zu3zaZuA+CizrfOvONbPOmfUttqNB2dawlbSSPbL36LTCStr1zueHA/17zTXNPZ2DfmlOJt/r/Cz5/4Nay0ebP5Ijex2Z1I4EzlrthbHaYJYP3PLtq+/fHB96m0tLABAANHkJACxNAFRAyHYaL9cpIt1fAWxj5cbka1Plpr3fP937vd6/KzfJ5k83J0fTFN4UBK2bgtShrQ9NQDT5Xr7ve63/b1/evt4n60zRrvtJO/P/2gRaVbeSJS8skdH/b3RyKrQ25GUArxROfdZdg4ZOCdY9KlQbGjNHivL9LHM0qfZ+uX5We5/kveuAakP/n4m1wTG1Xm/X7l2ybv066d6jewJC+d6rdjwZCFctM//O9bNkn32n/TPfG/szBYq0p/gz++W6ZEDf+6OKj6RL5y5y0EH7XnffB5TkyGyt0/6ZDzh1f6Z5kG3fzH6ZI7y59smAWm1YywVv2YAu87N8oJW59F/rnOdT35a6CQACgBb/CADoBwC1KTcGUrSpKaS9t+29mq/1H62XJX9cIu17tJcPKj9Ifv7BJx8kwKcgpKdMG7tpoe7UppN0bttZOrXd933f/yc/q/O7jm06JoBXF/IUJFvq2rOoDYJ5587vXB+SGuv/Yto/33rXrSHNdRlHoTXIN+9Cx9dS7w8AAoAmbwGAxQWAWqBrX+ORbXG1cSnYrf1obfK17uN1e79/tE7e++QvsKdgp6/XmE2Puh3W/rDkq1v7bnv/3a7Wv/Xn7fb+vGu7rgn0Kci1FLg1JvaG9o3aIJh3ceV3c/k51+uw3rHWGwAEAE01BQA8cAVDrx/KXIiea9EU7hTcMnBX8/3jvbCX+dLTm2mvkVNg63loz+Sre/vuUvlBpZw29DTp3bG39Di0hxx+yOHJzxX49HVLcaMxHjifF4N/WG/Wuxh82NIxAIAAoMljAOCBLZQKeO9/8r6s2bom69c7Fe/U3KGXb2EV3Hp36C1Hdjwy+dJ/9+rQS3oc0qMG+HQfPbWb2WiMB3a9861hS/+e9Wa9W9pjxfD6UX0OAAKApvwDAJu/Qejp2dWbV8tbW96qB3kKeHoxfM4He+577ISCXAbsjuywF/Bqfx3R4YjkDrfGblELJfNufp831nsHcn/Wm/U+kH4r1HsBgACgyXsAYPZCWfsOs7oCZ07TKuTVfG35y7/znZrVmyb0aF3/zv2zfino6Z1wLbHRGGmMLeGrYntNfI7Pi82TLREPAAgAmnwVHQDPOfccadO6/pE0hbwN2zbInzb9KSvk6cOJG9oU8AZ0GSADOg+oB3ktCXj5zEBjpDHm80gp/B6f4/NS8HG+OQCAAGA+jzT4+ygAWPeRD5kGcfbnz5Z3Pn5HVm5cmXyt2rSq5t8KeXq0Tp8hVnfTZ9b17dRXjul6TL2vo7ocJfq8umLcaIw0xmL0ZXPHhM/xeXN7qhhfDwAEAE2+LDUAzPZsL/2ZPsi4NuSt+HCFvPLnV+SDnR8kj13J9oR/vWNXT9Nmgzw9steUa/BMi9UMg2mMNMZmsFHRvwQ+x+dFb9JmCBAABABNNiolAFTQ08ekLPtwmSz7YFnyPXNUT2/MUKDLdjRPH1R8/GHHy3GHHSfHdzs++bd+Hd316Ba7Fs+0aIbBNEYao8E+bobic3zuxqyGQAFAANBgH3H5l0Ay1+e9/sHrNaCnsLf8w+XJ35nV07N1/7C4/qxf53574a7b8TKwy0DZ/KfNctkXLpMjOx9Z9A8yNi1yrcE0Rhpjc3mpmF8Hn+PzYvZnc8UGAAKAJi8V8xHAzDPzMkfzar5/uKzmIch1YU9P5x7b7VgZ3H1w8nVC9xMS6BvYdWDyR8IzGw2CBmFKHCeD8Tk+d2JVU5hRfQ4AAoCmxCkWANzy6RZ57YPX5LX3X9t7CnffadxNn27KOj89navX5g0+fC/o6deQw4fIwG4DU522jVowmDdAYCoYTgbjc3zuxKqmMAFAANBkoAMNgPq3afXZea++96r88f0/yh8/+GPy/c8Vf846Dz3Cp9fiZSBPgU9B77hux5luwqBB0CBMieNkMD7H506sagozqs8BQADQlDgtCYCbKjfthbx9X6++/2pyZG/7ru1ZY9bHqgw9fGgCeAnwHT44OX3bEo9UiVowmDdAYCoYTgbjc3zuxKqmMAFAANBkoOYAwJ27d8obm95IQE8hLwN86z5elzU2BToFvb/q8VfJ14k9TpShPYZK57adTXNpzGAaBA2iMX7xui8+x+devduYuKP6HAAsMQC8++67ZcqUKbJhwwYZPHiwTJs2Tc4888ycufDEE0/IzTffLG+++aYcffTRMnnyZDn//PNT505jAfDDTz6sB3p6VC/bc/Q0CH1eXgbyMsCnD0pudVCr1DG2xI5RCwbzBghaIp+K7TXxOT4vNk+2RDwAYAkB4KOPPiqXXHKJKASeccYZct9998kDDzwgy5cvl759+9bzz5IlSxI4vPXWWxPomzNnjtxyyy3ywgsvyIgRI1L5LRcAKtDpM/Rqn77Vf7+37b2sr3to60P3HtE7fO9RPf3So3od23RMFceB3okGQYM40J4rxPvhc3xeCN8d6PeM6nMAsIQAUKFt+PDhcs8999TkzwknnCBjx46V22+/vV5OXXTRRclz/ObNm1fzu3POOUe6dOkijzzySKoczADgY394TNZUrak5hat/KWPnnp31XiNzU0YG9k7seWICe/oXM/TPpnnZohYM5g0QeMlRS5z4HJ9b/ONlLABYIgC4Y8cOad++vTz22GP7ncK99tprZenSpbJgwYJ6ntSjghMmTEi+MtvUqVOT08bvvPNOVg9XVVWJfmU2BcA+ffqIfEdE2u4/RI/eDe2+91o9vWZPv/TmDD3a533TBjF//nwZPXq0lJeXe59O6viZN+ud2iyOd8Tn+NyxfVOHrgDYq1cvqaiokI4di/NsW+rJNHHHsmp9YrDzbf369dK7d29ZtGiRnH766TWzue2222TmzJmyatWqejNs3bq1zJgxQ8aNG1fzu9mzZ8v48eP3g7zaAydOnCiTJk2q91q9vt9LBnQZIP3b9pf+7fZ+dS/vHuYvZDi3D+GjAAqgAAoEU6CysjLp/wCg84XPAODixYvltNNOq5mN3tQxa9YsWblyZVYAVDi8+OKLa3738MMPyxVXXCHbt2d/1EquI4B600m3bt2cq5g+fI4QcIQgvVv87onP8blf96aPPKrPOQLIKeBGnQKum1KNvQs4fUoW955cI8Q1QsXt0OaJDp/j8+ZxUnG/SlSfcw1giQCgppfeBHLyyScndwFntkGDBsmYMWNy3gTy8ccfy9y5c2v2P/fcc6Vz586Nvglk48aN4Y4Aqm7nnUeDKO7S3jzRRW0QzJv8bp4MKu5XiepzALCEADDzGJh77703OQ08ffp0uf/++2XZsmXSr18/ufTSS5PrBDN3BOvp4pEjRybP/lNIfOqpp+Smm25qlsfAFHe626OLWjCYN0Bgz57ifwV8js+L36X2CAHAEgJAtYMe/bvjjjuSB0EPGTJE9K5ehTzdRo0aJf37909u/Mhsjz/+eAJ9b731Vs2DoC+44ILUzuIUMIUytVkc7wgQ4HPH9k0dOj6P5XMAsMQAMHWmN9OOAGCsgkGDYL2bqXQU9cvgc3xe1AZtpuAAQADQZCUAkEJpMpCTwQABPndiVVOY+DyWzwFAANBUMADAWAWDBsF6mwqGk8H4HJ87saopTAAQADQZCACkUJoM5GQwQIDPnVjVFCY+j+VzABAANBUMADBWwaBBsN6mguFkMD7H506sagoTAAQATQYCACmUJgM5GQwQ4HMnVjWFic9j+RwABABNBQMAjFUwaBCst6lgOBmMz/G5E6uawgQAAUCTgQBACqXJQE4GAwT43IlVTWHi81g+BwABQFPBAABjFQwaBOttKhhOBuNzfO7EqqYwAUAA0GQgAJBCaTKQk8EAAT53YlVTmPg8ls8BQADQVDAAwFgFgwbBepsKhpPB+ByfO7GqKUwAEAA0GQgApFCaDORkMECAz51Y1RQmPo/lcwAQADQVDAAwVsGgQbDepoLhZDA+x+dOrGoKEwAEAE0GAgAplCYDORkMEOBzJ1Y1hYnPY/kcAAQATQUDAIxVMGgQrLepYDgZjM/xuROrmsIEAAFAk4EAQAqlyUBOBgME+NyJVU1h4vNYPgcAAUBTwQAAYxUMGgTrbSoYTgbjc3zuxKqmMAFAANBkIACQQmkykJPBAAE+d2JVU5j4PJbPAUAA0FQwAMBYBYMGwXqbCoaTwfgcnzuxqilMABAANBkIAKRQmgzkZDBAgM+dWNUUJj6P5XMAEAA0FQwAMFbBoEGw3qaC4WQwPsfnTqxqChMABABNBgIAKZQmAzkZDBDgcydWNYWJz2P5HAAEAE0FAwCMVTBoEKy3qWA4GYzP8bkTq5rCBAABQJOBAEAKpclATgYDBPjciVVNYeLzWD4HAAFAU8EAAGMVDBoE620qGE4G43N87sSqpjABQADQZCAAkEJpMpCTwQABPndiVVOY+DyWzwFAANBUMADAWAWDBsF6mwqGk8H4HJ87saopTAAQADQZCACkUJoM5GQwQIDPnVjVFCY+j+VzABAANBUMADBWwaBBsN6mguFkMD7H506sagoTAAQATQYCACmUJgM5GQwQ4HMnVjWFic9j+RwABABNBQMAjFUwaBCst6lgOBmMz/G5E6uawgQAAUCTgQBACqXJQE4GAwT43IlVTWHi81g+BwABQFPBAABjFQwaBOttKhhOBuNzfO7EqqYwAUAA0GQgAJBCaTKQk8EAAT53YlVTmPg8ls8BQADQVDAAwFgFgwbBepsKhpPB+ByfO7GqKUwAEAA0GQgApFCaDORkMECAz51Y1RQmPo/lcwAQADQVDAAwVsGgQbDepoLhZDA+x+dOrGoKEwAEAE0GAgAplCYDORkMEOBzJ1Y1hYnPY/kcAAQATQUDAIxVMGgQrLepYDgZjM/xuROrmsIEAAFAk4EAQAqlyUBOBgME+NyJVU1h4vNYPgcAAUBTwQAAYxUMGgTrbSoYTgbjc3zuxKqmMAFAANBkIACQQmkykJPBAAE+d2JVU5j4PJbPAUAA0FQwAMBYBYMGwXqbCoaTwfgcnzuxqilMABAANBkIAKRQmgzkZDBAgM+dWNUUJj6P5XMAEAA0FQwAMFbBoEGw3qaC4WQwPsfnTqxqChMABABNBgIAKZQmAzkZDBDgcydWNYWJz2P5HAAEAE0FAwCMVTBoEKy3qWA4GYzP8bkTq5rCBAABQJOBAEAKpclATgYDBPjciVVNYeLzWD4HAAFAU8EAAGMVDBoE620qGE4G43N87sSqpjABQADQZCAAkEJpMpCTwQABPndiVVOY+DyWzwHAEgHALVu2yD/90z/J008/nRSAv/u7v5Of/vSn0rlz55wFYdSoUbJgwYL9fn/RRRfJL3/5y9RFBACMVTBoEKx36uLgeEd8js8d2zd16ABgiQDgueeeK2vXrpXp06cni3/VVVdJ//795ZlnnmkQAI899lj5wQ9+ULNPu3btpFOnTqkNBABSKFObxfGOAAE+d2zf1KHj81g+BwBLAABXrFghgwYNkhdffFFGjBiRJLv++7TTTpOVK1fKcccdl7UA6BHAk046SaZNm5a6QNTdEQCMVTBoEKx3k4uFo4H4HJ87smuTQwUASwAAH3zwQbn++utl69at+xlBT/9OnTpVxo8fnxMAly1bJtXV1dKjRw/Ro4jf//73pUOHDqkNBQBSKFObxfGOAAE+d2zf1KHj81g+BwBLAABvu+02mTFjhrzxxhv7Jbqe3lX4u/HGG7MWgPvvv18GDBggPXv2lNdffz3Z75hjjpH58+fnLBhVVVWiX5lNAbBPnz6yYcMG6datW+pC431HLZSq0+jRo6W8vNz7dFLHz7xZ79RmcbwjPsfnju2bOnQFwF69eklFRYV07Ngx9bhS2rGsWg+BFeE2ceJEmTRpUoORvfTSS/Lcc8/JzJkzZdWqVfvtO3DgQLniiivkO9/5TqrZvfzyy/KZz3xG9Pvw4cOzjskV0+zZs6V9+/ap3oedUAAFUAAFUAAFCqtAZWWljBs3DgAs7DJkf/eNGzeKfjW06Y0eCl9NOQVc93WVg9u0aSOzZs0SvRs428YRwL2qcISAIwTFWDOaOyZ8js+b21PF+HpRfc4RwBI4BZy5CeR//ud/5JRTTknyS/996qmnNngTSN1E1NPAQ4cOTR4NM3LkyFR5yjWAsa4Z4Roh1jtVYXC+Ez7H584tnCp8rgEsAQDUldYbONavXy/33XdfsvD6GJh+/frVPAZm3bp1cvbZZ8tDDz2UQOKbb74pDz/8sJx33nly2GGHyfLly+Vb3/qW6GNg9LRyq1atUhkIAKRQpjKK850AAnzu3MKpwsfnsXwOAJYIAG7evLneg6D/7d/+reZB0GvWrElu+Pjd734n+viXd999V/7hH/4huflj27ZtyY0cX/jCF5K7gLt27ZqqWOhOAGCsgkGDYL1TFwfHO+JzfO7YvqlDBwBLBABTr3gz7wgAUiib2VJF+XIAAT4vSmM2c1D4PJbPAUAA0FRCAMBYBYMGwXqbCoaTwfgcnzuxqilMABAANBkIAKRQmgzkZDBAgM+dWNUUJj6P5XMAEAA0FQwAMFbBoEGw3qaC4WQwPsfnTqxqChMABABNBgIAKZQmAzkZDBDgcydWNYWJz2P5HAAEAE0FAwCMVTBoEKy3qWA4GYzP8bkTq5rCBAABQJOBAEAKpclATgYDBPjciVVNYeLzWD4HAAFAU8EAAGMVDBoE620qGE4G43N87sSqpjABQADQZCAAkEJpMpCTwQABPndiVVOY+DyWzwFAANBUMADAWAWDBsF6mwqGk8H4HJ87saopTAAQADQZCACkUJoM5GQwQIDPnVjVFCY+j+VzABAANBUMADBWwaBBsN6mguFkMD7H506sagoTAAQATQYCACmUJgM5GQwQ4HMnVjWFic9j+RwABABNBQMAjFUwaBCst6lgOBmMz/G5E6uawgQAAUCTgQBACqXJQE4GAwT43IlVTWHi81g+BwABQFPBAABjFQwaBOttKhhOBuNzfO7EqqYwAUAA0GQgAJBCaTKQk8EAAT53YlVTmPg8ls8BQADQVDAAwFgFgwbBepsKhpPB+ByfO7GqKUwAEAA0GQgApFCaDORkMECAz51Y1RQmPo/lcwAQADQVDAAwVsGgQbDepoLhZDA+x+dOrGoKEwAEAE0GAgAplCYDORkMEOBzJ1Y1hYnPY/kcAAQATQUDAIxVMGgQrLepYDgZjM/xuROrmsIEAAFAk4EAQAqlyUBOBgME+NyJVU1h4vNYPgcAAUBTwQAAYxUMGgTrbSoYTgbjc3zuxKqmMAFAANBkIACQQmkykJPBAAE+d2JVU5j4PJbPAUAA0FQwAMBYBYMGwXqbCoaTwfgcnzuxqilMABAANBkIAKRQmgzkZDBAgM+dWNUUJj6P5XMAEAA0FQwAMFbBoEGw3qaC4WQwPsfnTqxqChMABABNBgIAKZQmAzkZDBDgcydWNYWJz2P5HAAEAE0FAwCMVTBoEKy3qWA4GYzP8bkTq5rCBAABQJOBAEAKpclATgYDBPjciVVNYeLzWD4HAAFAU8EAAGMVDBoE620qGE4G43N87sSqpjABQADQZCAAkEJpMpCTwQABPndiVVOY+DyWzwFAANBUMADAWAWDBsF6mwqGk8H4HJ87saopTAAQADQZCACkUJoM5GQwQIDPnVjVFCY+j+VzABAANBUMADBWwaBBsN6mguFkMD7H506sagoTAAQATQYCACmUJgM5GQwQ4HMnVjWFic9j+RwABABNBQMAjFUwaBCst6lgOBmMz/G5E6uawgQAAUCTgQBACqXJQE4GAwT43IlVTWHi81g+BwABQFPBAABjFQwaBOttKhhOBuNzfO7EqqYwAUAA0GQgAJBCaTKQk8EAAT53YlVTmPg8ls8BQADQVDAAwFgFgwbBepsKhpPB+ByfO7GqKUwAEAA0GQgApFCaDORkMECAz51Y1RQmPo/lcwAQADQVDAAwVsGgQbDepoLhZDA+x+dOrGoKEwAEAE0GAgAplCYDORkMEOBzJ1Y1hYnPY/kcAAQATQUDAIxVMGgQrLepYDgZjM/xuROrmsIEAAFAk4EAQAqlyUBOBgME+NyJVU1h4vNYPgcAAUBTwQAAYxUMGgTrbSoYTgbjc3zuxKqmMAFAANBkIACQQmkykJPBAAE+d2JVU5j4PJbPAUAA0FQwAMBYBYMGwXqbCoaTwfgcnzuxqilMABAANBkIAKRQmgzkZDBAgM+dWNUUJj6P5XMAsEQAcPLkyfLss8/K0qVLpXXr1rJ169a8haC6ulomTZok06dPly1btsiIESPkZz/7mQwePDjv2MwOAGCsgkGDYL1TFwfHO+JzfO7YvqlDBwBLBAC///3vS+fOnWXt2rXy7//+76kA8Mc//rEoOM6YMUOOPfZY+eEPfygLFy6UVatWSYcOHVKZCACkUKYyivOdAAJ87tzCqcLH57F8DgCWCABmslth7rrrrssLgHr074gjjkj2veGGG5LhVVVV0qNHD1EwvPrqq1MVDAAwVsGgQbDeqQqD853wOT53buFU4QOAQQHwrbfekqOPPlpeeeUVGTZsWI1ZxowZkxxJnDlzZlYDKSTqV2arqKiQvn37yhtvvCFdu3ZNZbpS2EkbxO9+ylBL4wAAD5pJREFU9zs566yzpLy8vBSmlGoOzJv1TmUU5zvhc3zu3MKpwt+8eXNy9k8vGevUqVOqMaW2U1m1Hg4rkS3tEcDFixfLGWecIevWrUuOBGa2q666St555x35zW9+k1WRiRMnJtcNsqEACqAACqAACvhX4M0335SjjjrK/0SaMIOiBcA0sPXSSy/JZz7zmZppNxYA169fL7169aoZf+WVV8q7774rv/71r1MdAdRPDv369ZM///nPoT5B6KnvPn36JFp17NixCbbzOYR5s94+ndu4qPE5Pm+cY3zunTmDpzeB6pm/iFvRAuDGjRtFvxra+vfvL23btm00ADb1FHDdWDLXAKqRooGQHjJn3jFKBj4nvyM4HZ/j8wg+rz3HogXApixE2iOAmZtAJkyYIN/+9reTt9qxY4ccfvjhTboJBBBqymr5G0ODoEH4c23jI8bn+LzxrvE3IqrPSw4A9RSsXtD59NNPy5QpU+T3v/99MsdjjjlGDj300OTfxx9/vNx+++1y/vnnJ/+vd/vq///85z+XgQMHym233SbPP/98kx4DAwD6S/6mRBy1YDBvgKAp+eJtDD7H5948a423JI4AXn755Vnv3NU7VUeNGpVoVFZWlsCe7qtb5kHQ9913334Pgh4yZEhqTfWOYIXIG2+8Udq0aZN6nPcdmTfr7d3DaeLH5/g8jU+874PPY/m85I4Aek9A4kcBFEABFEABFECBA6lASRwBPJCC8V4ogAIogAIogAIo4F0BAND7ChI/CqAACqAACqAACjRSAQCwkYKxOwqgAAqgAAqgAAp4VwAA9L6CxI8CKIACKIACKIACjVQAAMwj2N133508WmbDhg0yePBgmTZtmpx55pk5Rz3xxBNy8803i/55Gf17w5MnT6559Ewj16Ygu+tdzU8++aSsXLlS2rVrJ6effnryyJzjjjsuZzz6/MXx48fX+/2nn36634O6CzKhlG+a7S/P9OjRQ957772cr7BgwQK5/vrrZdmyZcmfFNRnSn7ta19L+Y7FsZs+TF3//GHd7ZprrpGf/exn9X7uda0XLlyY5PHLL7+c5PKcOXNk7NixNfPLPBVg+vTp+z0VQHO+oa2x9eFAr3pD89a/+XvTTTfJ3LlzRR+Orw93/5u/+Rv50Y9+tN+fyKwbc1NypZjmrbFke3LEiBEj5MUXX2ww1GKv7/l8rk/DyLbdcccd8i//8i9Zf+dhvdP0Lb3b+Z//+Z/lkUceEe1NZ599tmj+HnnkkTnXvKl14UD7vanvBwA2oNyjjz4ql1xySWIS/dvB+siYBx54QJYvXy59+/atN3LJkiUJHN56660J9GmTueWWW+SFF14QLS4etnPOOUe+8pWvyGc/+1nZtWuXfO9735PXXnstmfMhhxySdQoKBddee23yDMXaW8+ePT1MOYlRi9zjjz8u//Vf/1UTc6tWraR79+5Z5/D222+LPjJI/3zg1VdfLYsWLRKFJi0uF154oZt5f/jhh7J79+6aeF9//XUZPXq01H6EUu3JeF3refPmJWs0fPjwZH3qAqB+yNEPazo//QPxP/zhD0WbqXq6Q4cOWdezsfWhEKZoaN76/NIvfelLiYdPPPHEBHyvu+66JO//8Ic/5Ay3sblSbPPOAOD777+fPBoss7Vu3Vq6du2aM1wP9T2fz+t+oNX9r7jiClm9enXOv4frYb3T9K2vf/3r8swzzyQ53q1bN/nWt76VPD9YPxRqrc+2NaUuFMLvTX1PALAB5RTatGHcc889NXudcMIJyZED/cRRd7voootEHyaqSZXZ1JhdunRJwMDjpoCgfyFFj3aNHDkyJwBq49C/jex10yL3q1/9SpYuXZpqCjfccEPy4PEVK1bU7K9H/1599VXRRuF103X8z//8T/nTn/6UPDuz7pb2r+0U8/x1XrUBMPOXgXTuuq666dECPQKsDUABP9vW2PpQaE3qzjtbPPr31U855ZTkqHC2D7k6prG5Uozz1iOAWq8059Nu3up7mvXWXvbxxx/Lf//3f+eUwdt660Tq9i39sKMf5mfNmiW6jrqtX78++Zv2egT8b//2b+vNv6l1Ia2fimE/ADDHKuifhmvfvr089thj+53C1SNdCgkKRHU3LZj65+X0K7NNnTo1OW2c7TRbMRggXwz6yVD/UooeBcz1kGyFgq9+9avSu3fv5GjSSSedlBwFHTZsWL6XL5rfa5HTU4R6Gkwf6q3NXf86zFFHHZU1RoVhnd9dd91V83uFii9/+ctSWVkp5eXlRTO3tIGo5/VUtp7W/u53v5t1WCmsdd3G2JS/Dd6U+pB2HVpqvzRAoEfAP//5zydwlOvvmzc2V1pqPmlfN9u8FQAV/vSoX+fOneWv//qvkyPA+mE31+atvudbbz0Cqqc/Z86cKePGjWsQABtTG9OuS0vuV7dv/fa3v01O+eoRPz0gk9n0yLdC8KRJk+qF05S60JJzaonXBgBzqKqfDhRo9LSRXgeX2RQKNGHqnu7U32sx0QZZO5lmz56dXB+nRxS8bfoJaMyYMcmpocyf18s2B71uRhNu6NChyRFQhSL9VKVHwxQePWx61FbBTU//aWHUU4B6HaRe36enC+puup82kdqgtHjx4uRSAfVOr169PEx7vxj/4z/+I/Gu/mlFBcFSXeu6jTGzbuvWrdtv3ldddVXywe03v/lNPSmaUh8KbYh8QLB9+3b53Oc+l/zZzF/84hc5w21srhTjvPX0vf6Z0H79+olezqHXbeupbz0dmOuvOnmr7/nWW6/70+s91ctt27YtmfXO1rdy9WH9sDNgwIDk8q66W1PqQqG93tj3BwDzAKCa4LTTTqvZSz8l6mFkhYO6mxYIhcOLL7645lcPP/xwco2FFldv2ze+8Q159tlnk2sYG7pQtu689uzZk5w616NkP/nJT7xNO4n3k08+SW7i0Rs79IhYNgBUsNc/A5jZ9MOCNlC9ycDT9Y+Z+PU0iHpYr5NJu3lc61wAWBfc9dq4d999V37961/nBMDG1Ie0mrbUfg0Bgd4Q8vd///cJ/OvfRM919C9bbPlypaXmk/Z184GQvo7mrMLgL3/5S7nggguyvrS3+p5v3gr6er3vT3/607RSpqqNjXqxFtg5W9/KBYA6f63z9957b71Icn2gb6gutMB0WvQlAcAc8jblFI+3UwQNOeub3/xmcopEL4TXT0iN3TRJ1q5du9/1kI19jULvr8XhmGOO2e8a0ExMpXYKWI906eluvQNcj/o2ZvO21pwC/svdz7rOCn966YKe8tJTZdmOeOfzQ0O5km9sS/8+Hwhl3l/PVuilLJnrQOvG5a2+NzRvPaOjNUwvZ9LToI3dinW9c/UtTgFnX2EAsAHn63VgJ598cnIXcGYbNGhQ0iBz3QSiF9Tq6c/Mdu655ybXmHi5CUQPn2sS6fVseiSgKadw9TX0QnI9Jfzggw82trYUxf56yl4/GeppQL2Tu+6mTUKPlOnd0ZlN7zLTgurxJhC9rktPg+gRr4MPPjj1Gnhc61w3gei1u3rEVzf9AKjXg+W7CaQx9SG1qC20YzYgyMCf3vSjd37nuuu9oZDy5UoLTSf1y6YBwE2bNiWX/OhjgC699NKsr603D3iq7w3NWy9f0Tv+G7rbO5fAxbje+fpW5iYQvbRBP+zopkd99cxWvptAGlsXUhuzCHYEABtYhMxjHvTwsJ4G1uJw//33J9eF6ekCLRRaNDIwqIeM9VOVniZWSHzqqaeS52x5egyMPspED5dr7LWf/ac3R+hzAXWrO2+9gPbUU09NYFGvAdTTvnqaXE+JKgh62PT5UF/84heTOx8/+OCD5BpAvdFHb37RtdZTvXqN2EMPPZRMJ/MYGL1DVI+AKfTpXcDeHgOjc9HTuHqUVy9d0GuCam+lstbbtm1LrlPVTW/eufPOO+Wss85KHvuha66gp3msjwVRH+u1vvoBqPZjYPQicn280z/+4z8mr5OvPhSD7xuat17nqY/EeeWVV5I7v/Wu58ymuugpT93qzjtfrhT7vHVu+oFH567X6q5Zsya5lldPf+td/ZnH/nis7/l8rmujNVrn/a//+q9Zn1vqcb3T9C39gK4+1+v01QPqYwX/2o+B0dPiWgc0z3VLUxeKwe9NjQEAzKOcHv3Ti2X104LeBat39WYehzJq1CjRB+mqoTKbPktOoS9zB5HCYK5rSpq6aC05LteDQrUx6qdG3erOWz8h6alDfcaUgqI2WC2wta+dbMmYm+O19dmHerp748aNyVEQBVq9k1mP+Oqmc9dGoVCQ2RQQde6ZB0HrUUFvD4LWuTz33HPJYxAUdvTmltpbqay1rpsCX93tsssuS/I388BXPQqqNz3p0X99EHbtO98119UH6u3M1lB9aA5fWl+joXnrPHJd3lH7OZB1550vV6wxN8f4huatj/XSOz//7//+L7nbWWFIvaH5ro8FyWwe63s+n+vc9ECGPvJIe5rW67qbx/VO07f0Onx92LUe4Kj9IOjaa66vU7vXpakLzeHXQr0GAFgo5XlfFEABFEABFEABFCiQAgBggYTnbVEABVAABVAABVCgUAoAgIVSnvdFARRAARRAARRAgQIpAAAWSHjeFgVQAAVQAAVQAAUKpQAAWCjleV8UQAEUQAEUQAEUKJACAGCBhOdtUQAFUAAFUAAFUKBQCgCAhVKe90UBFEABFEABFECBAikAABZIeN4WBVAABVAABVAABQqlAABYKOV5XxRAARRAARRAARQokAIAYIGE521RAAVQAAVQAAVQoFAKAICFUp73RQEUQAEUQAEUQIECKQAAFkh43hYFUAAFUAAFUAAFCqUAAFgo5XlfFEABFEABFEABFCiQAgBggYTnbVEABVAABVAABVCgUAoAgIVSnvdFARRAARRAARRAgQIpAAAWSHjeFgVQAAVQAAVQAAUKpQAAWCjleV8UQAEUQAEUQAEUKJACAGCBhOdtUQAFUAAFUAAFUKBQCgCAhVKe90UBFEABFEABFECBAikAABZIeN4WBVAABVAABVAABQqlAABYKOV5XxRAARRAARRAARQokAIAYIGE521RAAVQAAVQAAVQoFAKAICFUp73RQEUQAEUQAEUQIECKQAAFkh43hYFUAAFUAAFUAAFCqUAAFgo5XlfFEABFEABFEABFCiQAgBggYTnbVEABVAABVAABVCgUAoAgIVSnvdFARRAARRAARRAgQIpAAAWSHjeFgVQAAVQAAVQAAUKpQAAWCjleV8UQAEUQAEUQAEUKJACAGCBhOdtUQAFUAAFUAAFUKBQCgCAhVKe90UBFEABFEABFECBAikAABZIeN4WBVAABVAABVAABQqlwP8HDycDeUuSsRoAAAAASUVORK5CYII=\" width=\"640\">"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "from matplotlib import pyplot as plt\n",
    "from matplotlib import animation\n",
    "%matplotlib notebook\n",
    "fig, ax = plt.subplots()\n",
    "ax.set_xlim(0,a)\n",
    "maxv=-1e10\n",
    "minv=1e10\n",
    "minv = [min(minv,min(t)) for t in u]\n",
    "minv=min(minv)\n",
    "maxv = [max(maxv,max(t)) for t in u]\n",
    "maxv=max(maxv)\n",
    "ax.set_xlim(0,a)\n",
    "ax.set_ylim(minv-1,maxv+1)\n",
    "ax.grid(True)\n",
    "line, = plt.plot(x, u[0], color='g')\n",
    "def animate(n):\n",
    "    #line, = plt.plot(x, u[n], color='g')\n",
    "    line.set_data(x, u[n])\n",
    "    return line,\n",
    "anim = animation.FuncAnimation(fig, animate, frames=len(t), interval=100,repeat = True)\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 27,
   "metadata": {
    "code_folding": []
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Wall time: 7.73 s\n"
     ]
    }
   ],
   "source": [
    "%%time\n",
    "\n",
    "def ux0(x):\n",
    "    return [0.01*t + exp(-(t-5)**4) for t in x]\n",
    "\n",
    "a = 20\n",
    "sigma = (1+1j)/2\n",
    "T = 20\n",
    "k = 2\n",
    "h = 0.01\n",
    "tau = 0.05\n",
    "\n",
    "x = np.arange(0, a+h, h)\n",
    "t = np.arange(0, T+tau, tau)\n",
    "u = np.zeros((len(t), len(x)))\n",
    "\n",
    "for i,ui in enumerate(ux0(x)):\n",
    "    u[0][i] = ui \n",
    "    \n",
    "diff_op = np.zeros((len(x),len(x)),dtype= np.float64)\n",
    "diff_op[0][0] = -2\n",
    "diff_op[0][1] = 2\n",
    "diff_op[-1][-2] = 2\n",
    "diff_op[-1][-1] = -2\n",
    "for n in range(1, len(x)-1):\n",
    "    diff_op[n][n-1] = 1\n",
    "    diff_op[n][n] = -2\n",
    "    diff_op[n][n+1] = 1\n",
    "diff_op = k/h**2 * diff_op\n",
    "A = np.eye(len(x),dtype= np.float64) - sigma*tau*diff_op;\n",
    "A = sparse.csc_matrix(A)\n",
    "\n",
    "for k in range(len(t)-1):\n",
    "    B = np.matrix(diff_op,dtype= np.float64)*np.matrix(u[k],dtype= np.float64).T\n",
    "    w = scipy.sparse.linalg.spsolve(A, B)\n",
    "    u[k+1] = np.matrix(u[k]) + tau*np.real(w)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 28,
   "metadata": {
    "scrolled": false
   },
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "Traceback (most recent call last):\n",
      "  File \"D:\\Programs\\Anaconda3\\lib\\site-packages\\matplotlib\\cbook\\__init__.py\", line 388, in process\n",
      "    proxy(*args, **kwargs)\n",
      "  File \"D:\\Programs\\Anaconda3\\lib\\site-packages\\matplotlib\\cbook\\__init__.py\", line 228, in __call__\n",
      "    return mtd(*args, **kwargs)\n",
      "  File \"D:\\Programs\\Anaconda3\\lib\\site-packages\\matplotlib\\animation.py\", line 1499, in _stop\n",
      "    self.event_source.remove_callback(self._loop_delay)\n",
      "AttributeError: 'NoneType' object has no attribute 'remove_callback'\n"
     ]
    },
    {
     "data": {
      "application/javascript": [
       "/* Put everything inside the global mpl namespace */\n",
       "window.mpl = {};\n",
       "\n",
       "\n",
       "mpl.get_websocket_type = function() {\n",
       "    if (typeof(WebSocket) !== 'undefined') {\n",
       "        return WebSocket;\n",
       "    } else if (typeof(MozWebSocket) !== 'undefined') {\n",
       "        return MozWebSocket;\n",
       "    } else {\n",
       "        alert('Your browser does not have WebSocket support.' +\n",
       "              'Please try Chrome, Safari or Firefox ≥ 6. ' +\n",
       "              'Firefox 4 and 5 are also supported but you ' +\n",
       "              'have to enable WebSockets in about:config.');\n",
       "    };\n",
       "}\n",
       "\n",
       "mpl.figure = function(figure_id, websocket, ondownload, parent_element) {\n",
       "    this.id = figure_id;\n",
       "\n",
       "    this.ws = websocket;\n",
       "\n",
       "    this.supports_binary = (this.ws.binaryType != undefined);\n",
       "\n",
       "    if (!this.supports_binary) {\n",
       "        var warnings = document.getElementById(\"mpl-warnings\");\n",
       "        if (warnings) {\n",
       "            warnings.style.display = 'block';\n",
       "            warnings.textContent = (\n",
       "                \"This browser does not support binary websocket messages. \" +\n",
       "                    \"Performance may be slow.\");\n",
       "        }\n",
       "    }\n",
       "\n",
       "    this.imageObj = new Image();\n",
       "\n",
       "    this.context = undefined;\n",
       "    this.message = undefined;\n",
       "    this.canvas = undefined;\n",
       "    this.rubberband_canvas = undefined;\n",
       "    this.rubberband_context = undefined;\n",
       "    this.format_dropdown = undefined;\n",
       "\n",
       "    this.image_mode = 'full';\n",
       "\n",
       "    this.root = $('<div/>');\n",
       "    this._root_extra_style(this.root)\n",
       "    this.root.attr('style', 'display: inline-block');\n",
       "\n",
       "    $(parent_element).append(this.root);\n",
       "\n",
       "    this._init_header(this);\n",
       "    this._init_canvas(this);\n",
       "    this._init_toolbar(this);\n",
       "\n",
       "    var fig = this;\n",
       "\n",
       "    this.waiting = false;\n",
       "\n",
       "    this.ws.onopen =  function () {\n",
       "            fig.send_message(\"supports_binary\", {value: fig.supports_binary});\n",
       "            fig.send_message(\"send_image_mode\", {});\n",
       "            if (mpl.ratio != 1) {\n",
       "                fig.send_message(\"set_dpi_ratio\", {'dpi_ratio': mpl.ratio});\n",
       "            }\n",
       "            fig.send_message(\"refresh\", {});\n",
       "        }\n",
       "\n",
       "    this.imageObj.onload = function() {\n",
       "            if (fig.image_mode == 'full') {\n",
       "                // Full images could contain transparency (where diff images\n",
       "                // almost always do), so we need to clear the canvas so that\n",
       "                // there is no ghosting.\n",
       "                fig.context.clearRect(0, 0, fig.canvas.width, fig.canvas.height);\n",
       "            }\n",
       "            fig.context.drawImage(fig.imageObj, 0, 0);\n",
       "        };\n",
       "\n",
       "    this.imageObj.onunload = function() {\n",
       "        fig.ws.close();\n",
       "    }\n",
       "\n",
       "    this.ws.onmessage = this._make_on_message_function(this);\n",
       "\n",
       "    this.ondownload = ondownload;\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._init_header = function() {\n",
       "    var titlebar = $(\n",
       "        '<div class=\"ui-dialog-titlebar ui-widget-header ui-corner-all ' +\n",
       "        'ui-helper-clearfix\"/>');\n",
       "    var titletext = $(\n",
       "        '<div class=\"ui-dialog-title\" style=\"width: 100%; ' +\n",
       "        'text-align: center; padding: 3px;\"/>');\n",
       "    titlebar.append(titletext)\n",
       "    this.root.append(titlebar);\n",
       "    this.header = titletext[0];\n",
       "}\n",
       "\n",
       "\n",
       "\n",
       "mpl.figure.prototype._canvas_extra_style = function(canvas_div) {\n",
       "\n",
       "}\n",
       "\n",
       "\n",
       "mpl.figure.prototype._root_extra_style = function(canvas_div) {\n",
       "\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._init_canvas = function() {\n",
       "    var fig = this;\n",
       "\n",
       "    var canvas_div = $('<div/>');\n",
       "\n",
       "    canvas_div.attr('style', 'position: relative; clear: both; outline: 0');\n",
       "\n",
       "    function canvas_keyboard_event(event) {\n",
       "        return fig.key_event(event, event['data']);\n",
       "    }\n",
       "\n",
       "    canvas_div.keydown('key_press', canvas_keyboard_event);\n",
       "    canvas_div.keyup('key_release', canvas_keyboard_event);\n",
       "    this.canvas_div = canvas_div\n",
       "    this._canvas_extra_style(canvas_div)\n",
       "    this.root.append(canvas_div);\n",
       "\n",
       "    var canvas = $('<canvas/>');\n",
       "    canvas.addClass('mpl-canvas');\n",
       "    canvas.attr('style', \"left: 0; top: 0; z-index: 0; outline: 0\")\n",
       "\n",
       "    this.canvas = canvas[0];\n",
       "    this.context = canvas[0].getContext(\"2d\");\n",
       "\n",
       "    var backingStore = this.context.backingStorePixelRatio ||\n",
       "\tthis.context.webkitBackingStorePixelRatio ||\n",
       "\tthis.context.mozBackingStorePixelRatio ||\n",
       "\tthis.context.msBackingStorePixelRatio ||\n",
       "\tthis.context.oBackingStorePixelRatio ||\n",
       "\tthis.context.backingStorePixelRatio || 1;\n",
       "\n",
       "    mpl.ratio = (window.devicePixelRatio || 1) / backingStore;\n",
       "\n",
       "    var rubberband = $('<canvas/>');\n",
       "    rubberband.attr('style', \"position: absolute; left: 0; top: 0; z-index: 1;\")\n",
       "\n",
       "    var pass_mouse_events = true;\n",
       "\n",
       "    canvas_div.resizable({\n",
       "        start: function(event, ui) {\n",
       "            pass_mouse_events = false;\n",
       "        },\n",
       "        resize: function(event, ui) {\n",
       "            fig.request_resize(ui.size.width, ui.size.height);\n",
       "        },\n",
       "        stop: function(event, ui) {\n",
       "            pass_mouse_events = true;\n",
       "            fig.request_resize(ui.size.width, ui.size.height);\n",
       "        },\n",
       "    });\n",
       "\n",
       "    function mouse_event_fn(event) {\n",
       "        if (pass_mouse_events)\n",
       "            return fig.mouse_event(event, event['data']);\n",
       "    }\n",
       "\n",
       "    rubberband.mousedown('button_press', mouse_event_fn);\n",
       "    rubberband.mouseup('button_release', mouse_event_fn);\n",
       "    // Throttle sequential mouse events to 1 every 20ms.\n",
       "    rubberband.mousemove('motion_notify', mouse_event_fn);\n",
       "\n",
       "    rubberband.mouseenter('figure_enter', mouse_event_fn);\n",
       "    rubberband.mouseleave('figure_leave', mouse_event_fn);\n",
       "\n",
       "    canvas_div.on(\"wheel\", function (event) {\n",
       "        event = event.originalEvent;\n",
       "        event['data'] = 'scroll'\n",
       "        if (event.deltaY < 0) {\n",
       "            event.step = 1;\n",
       "        } else {\n",
       "            event.step = -1;\n",
       "        }\n",
       "        mouse_event_fn(event);\n",
       "    });\n",
       "\n",
       "    canvas_div.append(canvas);\n",
       "    canvas_div.append(rubberband);\n",
       "\n",
       "    this.rubberband = rubberband;\n",
       "    this.rubberband_canvas = rubberband[0];\n",
       "    this.rubberband_context = rubberband[0].getContext(\"2d\");\n",
       "    this.rubberband_context.strokeStyle = \"#000000\";\n",
       "\n",
       "    this._resize_canvas = function(width, height) {\n",
       "        // Keep the size of the canvas, canvas container, and rubber band\n",
       "        // canvas in synch.\n",
       "        canvas_div.css('width', width)\n",
       "        canvas_div.css('height', height)\n",
       "\n",
       "        canvas.attr('width', width * mpl.ratio);\n",
       "        canvas.attr('height', height * mpl.ratio);\n",
       "        canvas.attr('style', 'width: ' + width + 'px; height: ' + height + 'px;');\n",
       "\n",
       "        rubberband.attr('width', width);\n",
       "        rubberband.attr('height', height);\n",
       "    }\n",
       "\n",
       "    // Set the figure to an initial 600x600px, this will subsequently be updated\n",
       "    // upon first draw.\n",
       "    this._resize_canvas(600, 600);\n",
       "\n",
       "    // Disable right mouse context menu.\n",
       "    $(this.rubberband_canvas).bind(\"contextmenu\",function(e){\n",
       "        return false;\n",
       "    });\n",
       "\n",
       "    function set_focus () {\n",
       "        canvas.focus();\n",
       "        canvas_div.focus();\n",
       "    }\n",
       "\n",
       "    window.setTimeout(set_focus, 100);\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._init_toolbar = function() {\n",
       "    var fig = this;\n",
       "\n",
       "    var nav_element = $('<div/>')\n",
       "    nav_element.attr('style', 'width: 100%');\n",
       "    this.root.append(nav_element);\n",
       "\n",
       "    // Define a callback function for later on.\n",
       "    function toolbar_event(event) {\n",
       "        return fig.toolbar_button_onclick(event['data']);\n",
       "    }\n",
       "    function toolbar_mouse_event(event) {\n",
       "        return fig.toolbar_button_onmouseover(event['data']);\n",
       "    }\n",
       "\n",
       "    for(var toolbar_ind in mpl.toolbar_items) {\n",
       "        var name = mpl.toolbar_items[toolbar_ind][0];\n",
       "        var tooltip = mpl.toolbar_items[toolbar_ind][1];\n",
       "        var image = mpl.toolbar_items[toolbar_ind][2];\n",
       "        var method_name = mpl.toolbar_items[toolbar_ind][3];\n",
       "\n",
       "        if (!name) {\n",
       "            // put a spacer in here.\n",
       "            continue;\n",
       "        }\n",
       "        var button = $('<button/>');\n",
       "        button.addClass('ui-button ui-widget ui-state-default ui-corner-all ' +\n",
       "                        'ui-button-icon-only');\n",
       "        button.attr('role', 'button');\n",
       "        button.attr('aria-disabled', 'false');\n",
       "        button.click(method_name, toolbar_event);\n",
       "        button.mouseover(tooltip, toolbar_mouse_event);\n",
       "\n",
       "        var icon_img = $('<span/>');\n",
       "        icon_img.addClass('ui-button-icon-primary ui-icon');\n",
       "        icon_img.addClass(image);\n",
       "        icon_img.addClass('ui-corner-all');\n",
       "\n",
       "        var tooltip_span = $('<span/>');\n",
       "        tooltip_span.addClass('ui-button-text');\n",
       "        tooltip_span.html(tooltip);\n",
       "\n",
       "        button.append(icon_img);\n",
       "        button.append(tooltip_span);\n",
       "\n",
       "        nav_element.append(button);\n",
       "    }\n",
       "\n",
       "    var fmt_picker_span = $('<span/>');\n",
       "\n",
       "    var fmt_picker = $('<select/>');\n",
       "    fmt_picker.addClass('mpl-toolbar-option ui-widget ui-widget-content');\n",
       "    fmt_picker_span.append(fmt_picker);\n",
       "    nav_element.append(fmt_picker_span);\n",
       "    this.format_dropdown = fmt_picker[0];\n",
       "\n",
       "    for (var ind in mpl.extensions) {\n",
       "        var fmt = mpl.extensions[ind];\n",
       "        var option = $(\n",
       "            '<option/>', {selected: fmt === mpl.default_extension}).html(fmt);\n",
       "        fmt_picker.append(option)\n",
       "    }\n",
       "\n",
       "    // Add hover states to the ui-buttons\n",
       "    $( \".ui-button\" ).hover(\n",
       "        function() { $(this).addClass(\"ui-state-hover\");},\n",
       "        function() { $(this).removeClass(\"ui-state-hover\");}\n",
       "    );\n",
       "\n",
       "    var status_bar = $('<span class=\"mpl-message\"/>');\n",
       "    nav_element.append(status_bar);\n",
       "    this.message = status_bar[0];\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.request_resize = function(x_pixels, y_pixels) {\n",
       "    // Request matplotlib to resize the figure. Matplotlib will then trigger a resize in the client,\n",
       "    // which will in turn request a refresh of the image.\n",
       "    this.send_message('resize', {'width': x_pixels, 'height': y_pixels});\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.send_message = function(type, properties) {\n",
       "    properties['type'] = type;\n",
       "    properties['figure_id'] = this.id;\n",
       "    this.ws.send(JSON.stringify(properties));\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.send_draw_message = function() {\n",
       "    if (!this.waiting) {\n",
       "        this.waiting = true;\n",
       "        this.ws.send(JSON.stringify({type: \"draw\", figure_id: this.id}));\n",
       "    }\n",
       "}\n",
       "\n",
       "\n",
       "mpl.figure.prototype.handle_save = function(fig, msg) {\n",
       "    var format_dropdown = fig.format_dropdown;\n",
       "    var format = format_dropdown.options[format_dropdown.selectedIndex].value;\n",
       "    fig.ondownload(fig, format);\n",
       "}\n",
       "\n",
       "\n",
       "mpl.figure.prototype.handle_resize = function(fig, msg) {\n",
       "    var size = msg['size'];\n",
       "    if (size[0] != fig.canvas.width || size[1] != fig.canvas.height) {\n",
       "        fig._resize_canvas(size[0], size[1]);\n",
       "        fig.send_message(\"refresh\", {});\n",
       "    };\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_rubberband = function(fig, msg) {\n",
       "    var x0 = msg['x0'] / mpl.ratio;\n",
       "    var y0 = (fig.canvas.height - msg['y0']) / mpl.ratio;\n",
       "    var x1 = msg['x1'] / mpl.ratio;\n",
       "    var y1 = (fig.canvas.height - msg['y1']) / mpl.ratio;\n",
       "    x0 = Math.floor(x0) + 0.5;\n",
       "    y0 = Math.floor(y0) + 0.5;\n",
       "    x1 = Math.floor(x1) + 0.5;\n",
       "    y1 = Math.floor(y1) + 0.5;\n",
       "    var min_x = Math.min(x0, x1);\n",
       "    var min_y = Math.min(y0, y1);\n",
       "    var width = Math.abs(x1 - x0);\n",
       "    var height = Math.abs(y1 - y0);\n",
       "\n",
       "    fig.rubberband_context.clearRect(\n",
       "        0, 0, fig.canvas.width, fig.canvas.height);\n",
       "\n",
       "    fig.rubberband_context.strokeRect(min_x, min_y, width, height);\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_figure_label = function(fig, msg) {\n",
       "    // Updates the figure title.\n",
       "    fig.header.textContent = msg['label'];\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_cursor = function(fig, msg) {\n",
       "    var cursor = msg['cursor'];\n",
       "    switch(cursor)\n",
       "    {\n",
       "    case 0:\n",
       "        cursor = 'pointer';\n",
       "        break;\n",
       "    case 1:\n",
       "        cursor = 'default';\n",
       "        break;\n",
       "    case 2:\n",
       "        cursor = 'crosshair';\n",
       "        break;\n",
       "    case 3:\n",
       "        cursor = 'move';\n",
       "        break;\n",
       "    }\n",
       "    fig.rubberband_canvas.style.cursor = cursor;\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_message = function(fig, msg) {\n",
       "    fig.message.textContent = msg['message'];\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_draw = function(fig, msg) {\n",
       "    // Request the server to send over a new figure.\n",
       "    fig.send_draw_message();\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_image_mode = function(fig, msg) {\n",
       "    fig.image_mode = msg['mode'];\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.updated_canvas_event = function() {\n",
       "    // Called whenever the canvas gets updated.\n",
       "    this.send_message(\"ack\", {});\n",
       "}\n",
       "\n",
       "// A function to construct a web socket function for onmessage handling.\n",
       "// Called in the figure constructor.\n",
       "mpl.figure.prototype._make_on_message_function = function(fig) {\n",
       "    return function socket_on_message(evt) {\n",
       "        if (evt.data instanceof Blob) {\n",
       "            /* FIXME: We get \"Resource interpreted as Image but\n",
       "             * transferred with MIME type text/plain:\" errors on\n",
       "             * Chrome.  But how to set the MIME type?  It doesn't seem\n",
       "             * to be part of the websocket stream */\n",
       "            evt.data.type = \"image/png\";\n",
       "\n",
       "            /* Free the memory for the previous frames */\n",
       "            if (fig.imageObj.src) {\n",
       "                (window.URL || window.webkitURL).revokeObjectURL(\n",
       "                    fig.imageObj.src);\n",
       "            }\n",
       "\n",
       "            fig.imageObj.src = (window.URL || window.webkitURL).createObjectURL(\n",
       "                evt.data);\n",
       "            fig.updated_canvas_event();\n",
       "            fig.waiting = false;\n",
       "            return;\n",
       "        }\n",
       "        else if (typeof evt.data === 'string' && evt.data.slice(0, 21) == \"data:image/png;base64\") {\n",
       "            fig.imageObj.src = evt.data;\n",
       "            fig.updated_canvas_event();\n",
       "            fig.waiting = false;\n",
       "            return;\n",
       "        }\n",
       "\n",
       "        var msg = JSON.parse(evt.data);\n",
       "        var msg_type = msg['type'];\n",
       "\n",
       "        // Call the  \"handle_{type}\" callback, which takes\n",
       "        // the figure and JSON message as its only arguments.\n",
       "        try {\n",
       "            var callback = fig[\"handle_\" + msg_type];\n",
       "        } catch (e) {\n",
       "            console.log(\"No handler for the '\" + msg_type + \"' message type: \", msg);\n",
       "            return;\n",
       "        }\n",
       "\n",
       "        if (callback) {\n",
       "            try {\n",
       "                // console.log(\"Handling '\" + msg_type + \"' message: \", msg);\n",
       "                callback(fig, msg);\n",
       "            } catch (e) {\n",
       "                console.log(\"Exception inside the 'handler_\" + msg_type + \"' callback:\", e, e.stack, msg);\n",
       "            }\n",
       "        }\n",
       "    };\n",
       "}\n",
       "\n",
       "// from http://stackoverflow.com/questions/1114465/getting-mouse-location-in-canvas\n",
       "mpl.findpos = function(e) {\n",
       "    //this section is from http://www.quirksmode.org/js/events_properties.html\n",
       "    var targ;\n",
       "    if (!e)\n",
       "        e = window.event;\n",
       "    if (e.target)\n",
       "        targ = e.target;\n",
       "    else if (e.srcElement)\n",
       "        targ = e.srcElement;\n",
       "    if (targ.nodeType == 3) // defeat Safari bug\n",
       "        targ = targ.parentNode;\n",
       "\n",
       "    // jQuery normalizes the pageX and pageY\n",
       "    // pageX,Y are the mouse positions relative to the document\n",
       "    // offset() returns the position of the element relative to the document\n",
       "    var x = e.pageX - $(targ).offset().left;\n",
       "    var y = e.pageY - $(targ).offset().top;\n",
       "\n",
       "    return {\"x\": x, \"y\": y};\n",
       "};\n",
       "\n",
       "/*\n",
       " * return a copy of an object with only non-object keys\n",
       " * we need this to avoid circular references\n",
       " * http://stackoverflow.com/a/24161582/3208463\n",
       " */\n",
       "function simpleKeys (original) {\n",
       "  return Object.keys(original).reduce(function (obj, key) {\n",
       "    if (typeof original[key] !== 'object')\n",
       "        obj[key] = original[key]\n",
       "    return obj;\n",
       "  }, {});\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.mouse_event = function(event, name) {\n",
       "    var canvas_pos = mpl.findpos(event)\n",
       "\n",
       "    if (name === 'button_press')\n",
       "    {\n",
       "        this.canvas.focus();\n",
       "        this.canvas_div.focus();\n",
       "    }\n",
       "\n",
       "    var x = canvas_pos.x * mpl.ratio;\n",
       "    var y = canvas_pos.y * mpl.ratio;\n",
       "\n",
       "    this.send_message(name, {x: x, y: y, button: event.button,\n",
       "                             step: event.step,\n",
       "                             guiEvent: simpleKeys(event)});\n",
       "\n",
       "    /* This prevents the web browser from automatically changing to\n",
       "     * the text insertion cursor when the button is pressed.  We want\n",
       "     * to control all of the cursor setting manually through the\n",
       "     * 'cursor' event from matplotlib */\n",
       "    event.preventDefault();\n",
       "    return false;\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._key_event_extra = function(event, name) {\n",
       "    // Handle any extra behaviour associated with a key event\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.key_event = function(event, name) {\n",
       "\n",
       "    // Prevent repeat events\n",
       "    if (name == 'key_press')\n",
       "    {\n",
       "        if (event.which === this._key)\n",
       "            return;\n",
       "        else\n",
       "            this._key = event.which;\n",
       "    }\n",
       "    if (name == 'key_release')\n",
       "        this._key = null;\n",
       "\n",
       "    var value = '';\n",
       "    if (event.ctrlKey && event.which != 17)\n",
       "        value += \"ctrl+\";\n",
       "    if (event.altKey && event.which != 18)\n",
       "        value += \"alt+\";\n",
       "    if (event.shiftKey && event.which != 16)\n",
       "        value += \"shift+\";\n",
       "\n",
       "    value += 'k';\n",
       "    value += event.which.toString();\n",
       "\n",
       "    this._key_event_extra(event, name);\n",
       "\n",
       "    this.send_message(name, {key: value,\n",
       "                             guiEvent: simpleKeys(event)});\n",
       "    return false;\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.toolbar_button_onclick = function(name) {\n",
       "    if (name == 'download') {\n",
       "        this.handle_save(this, null);\n",
       "    } else {\n",
       "        this.send_message(\"toolbar_button\", {name: name});\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.toolbar_button_onmouseover = function(tooltip) {\n",
       "    this.message.textContent = tooltip;\n",
       "};\n",
       "mpl.toolbar_items = [[\"Home\", \"Reset original view\", \"fa fa-home icon-home\", \"home\"], [\"Back\", \"Back to  previous view\", \"fa fa-arrow-left icon-arrow-left\", \"back\"], [\"Forward\", \"Forward to next view\", \"fa fa-arrow-right icon-arrow-right\", \"forward\"], [\"\", \"\", \"\", \"\"], [\"Pan\", \"Pan axes with left mouse, zoom with right\", \"fa fa-arrows icon-move\", \"pan\"], [\"Zoom\", \"Zoom to rectangle\", \"fa fa-square-o icon-check-empty\", \"zoom\"], [\"\", \"\", \"\", \"\"], [\"Download\", \"Download plot\", \"fa fa-floppy-o icon-save\", \"download\"]];\n",
       "\n",
       "mpl.extensions = [\"eps\", \"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 overridden (by mpl) onmessage function.\n",
       "        ws.onmessage(msg['content']['data'])\n",
       "    });\n",
       "    return ws;\n",
       "}\n",
       "\n",
       "mpl.mpl_figure_comm = function(comm, msg) {\n",
       "    // This is the function which gets called when the mpl process\n",
       "    // starts-up an IPython Comm through the \"matplotlib\" channel.\n",
       "\n",
       "    var id = msg.content.data.id;\n",
       "    // Get hold of the div created by the display call when the Comm\n",
       "    // socket was opened in Python.\n",
       "    var element = $(\"#\" + id);\n",
       "    var ws_proxy = comm_websocket_adapter(comm)\n",
       "\n",
       "    function ondownload(figure, format) {\n",
       "        window.open(figure.imageObj.src);\n",
       "    }\n",
       "\n",
       "    var fig = new mpl.figure(id, ws_proxy,\n",
       "                           ondownload,\n",
       "                           element.get(0));\n",
       "\n",
       "    // Call onopen now - mpl needs it, as it is assuming we've passed it a real\n",
       "    // web socket which is closed, not our websocket->open comm proxy.\n",
       "    ws_proxy.onopen();\n",
       "\n",
       "    fig.parent_element = element.get(0);\n",
       "    fig.cell_info = mpl.find_output_cell(\"<div id='\" + id + \"'></div>\");\n",
       "    if (!fig.cell_info) {\n",
       "        console.error(\"Failed to find cell for figure\", id, fig);\n",
       "        return;\n",
       "    }\n",
       "\n",
       "    var output_index = fig.cell_info[2]\n",
       "    var cell = fig.cell_info[0];\n",
       "\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.handle_close = function(fig, msg) {\n",
       "    var width = fig.canvas.width/mpl.ratio\n",
       "    fig.root.unbind('remove')\n",
       "\n",
       "    // Update the output cell to use the data from the current canvas.\n",
       "    fig.push_to_output();\n",
       "    var dataURL = fig.canvas.toDataURL();\n",
       "    // Re-enable the keyboard manager in IPython - without this line, in FF,\n",
       "    // the notebook keyboard shortcuts fail.\n",
       "    IPython.keyboard_manager.enable()\n",
       "    $(fig.parent_element).html('<img src=\"' + dataURL + '\" width=\"' + width + '\">');\n",
       "    fig.close_ws(fig, msg);\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.close_ws = function(fig, msg){\n",
       "    fig.send_message('closing', msg);\n",
       "    // fig.ws.close()\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.push_to_output = function(remove_interactive) {\n",
       "    // Turn the data on the canvas into data in the output cell.\n",
       "    var width = this.canvas.width/mpl.ratio\n",
       "    var dataURL = this.canvas.toDataURL();\n",
       "    this.cell_info[1]['text/html'] = '<img src=\"' + dataURL + '\" width=\"' + width + '\">';\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.updated_canvas_event = function() {\n",
       "    // Tell IPython that the notebook contents must change.\n",
       "    IPython.notebook.set_dirty(true);\n",
       "    this.send_message(\"ack\", {});\n",
       "    var fig = this;\n",
       "    // Wait a second, then push the new image to the DOM so\n",
       "    // that it is saved nicely (might be nice to debounce this).\n",
       "    setTimeout(function () { fig.push_to_output() }, 1000);\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._init_toolbar = function() {\n",
       "    var fig = this;\n",
       "\n",
       "    var nav_element = $('<div/>')\n",
       "    nav_element.attr('style', 'width: 100%');\n",
       "    this.root.append(nav_element);\n",
       "\n",
       "    // Define a callback function for later on.\n",
       "    function toolbar_event(event) {\n",
       "        return fig.toolbar_button_onclick(event['data']);\n",
       "    }\n",
       "    function toolbar_mouse_event(event) {\n",
       "        return fig.toolbar_button_onmouseover(event['data']);\n",
       "    }\n",
       "\n",
       "    for(var toolbar_ind in mpl.toolbar_items){\n",
       "        var name = mpl.toolbar_items[toolbar_ind][0];\n",
       "        var tooltip = mpl.toolbar_items[toolbar_ind][1];\n",
       "        var image = mpl.toolbar_items[toolbar_ind][2];\n",
       "        var method_name = mpl.toolbar_items[toolbar_ind][3];\n",
       "\n",
       "        if (!name) { continue; };\n",
       "\n",
       "        var button = $('<button class=\"btn btn-default\" href=\"#\" title=\"' + name + '\"><i class=\"fa ' + image + ' fa-lg\"></i></button>');\n",
       "        button.click(method_name, toolbar_event);\n",
       "        button.mouseover(tooltip, toolbar_mouse_event);\n",
       "        nav_element.append(button);\n",
       "    }\n",
       "\n",
       "    // Add the status bar.\n",
       "    var status_bar = $('<span class=\"mpl-message\" style=\"text-align:right; float: right;\"/>');\n",
       "    nav_element.append(status_bar);\n",
       "    this.message = status_bar[0];\n",
       "\n",
       "    // Add the close button to the window.\n",
       "    var buttongrp = $('<div class=\"btn-group inline pull-right\"></div>');\n",
       "    var button = $('<button class=\"btn btn-mini btn-primary\" href=\"#\" title=\"Stop Interaction\"><i class=\"fa fa-power-off icon-remove icon-large\"></i></button>');\n",
       "    button.click(function (evt) { fig.handle_close(fig, {}); } );\n",
       "    button.mouseover('Stop Interaction', toolbar_mouse_event);\n",
       "    buttongrp.append(button);\n",
       "    var titlebar = this.root.find($('.ui-dialog-titlebar'));\n",
       "    titlebar.prepend(buttongrp);\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._root_extra_style = function(el){\n",
       "    var fig = this\n",
       "    el.on(\"remove\", function(){\n",
       "\tfig.close_ws(fig, {});\n",
       "    });\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._canvas_extra_style = function(el){\n",
       "    // this is important to make the div 'focusable\n",
       "    el.attr('tabindex', 0)\n",
       "    // reach out to IPython and tell the keyboard manager to turn it's self\n",
       "    // off when our div gets focus\n",
       "\n",
       "    // location in version 3\n",
       "    if (IPython.notebook.keyboard_manager) {\n",
       "        IPython.notebook.keyboard_manager.register_events(el);\n",
       "    }\n",
       "    else {\n",
       "        // location in version 2\n",
       "        IPython.keyboard_manager.register_events(el);\n",
       "    }\n",
       "\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._key_event_extra = function(event, name) {\n",
       "    var manager = IPython.notebook.keyboard_manager;\n",
       "    if (!manager)\n",
       "        manager = IPython.keyboard_manager;\n",
       "\n",
       "    // Check for shift+enter\n",
       "    if (event.shiftKey && event.which == 13) {\n",
       "        this.canvas_div.blur();\n",
       "        event.shiftKey = false;\n",
       "        // Send a \"J\" for go to next cell\n",
       "        event.which = 74;\n",
       "        event.keyCode = 74;\n",
       "        manager.command_mode();\n",
       "        manager.handle_keydown(event);\n",
       "    }\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_save = function(fig, msg) {\n",
       "    fig.ondownload(fig, null);\n",
       "}\n",
       "\n",
       "\n",
       "mpl.find_output_cell = function(html_output) {\n",
       "    // Return the cell and output element which can be found *uniquely* in the notebook.\n",
       "    // Note - this is a bit hacky, but it is done because the \"notebook_saving.Notebook\"\n",
       "    // IPython event is triggered only after the cells have been serialised, which for\n",
       "    // our purposes (turning an active figure into a static one), is too late.\n",
       "    var cells = IPython.notebook.get_cells();\n",
       "    var ncells = cells.length;\n",
       "    for (var i=0; i<ncells; i++) {\n",
       "        var cell = cells[i];\n",
       "        if (cell.cell_type === 'code'){\n",
       "            for (var j=0; j<cell.output_area.outputs.length; j++) {\n",
       "                var data = cell.output_area.outputs[j];\n",
       "                if (data.data) {\n",
       "                    // IPython >= 3 moved mimebundle to data attribute of output\n",
       "                    data = data.data;\n",
       "                }\n",
       "                if (data['text/html'] == html_output) {\n",
       "                    return [cell, data, j];\n",
       "                }\n",
       "            }\n",
       "        }\n",
       "    }\n",
       "}\n",
       "\n",
       "// Register the function which deals with the matplotlib target/channel.\n",
       "// The kernel may be null if the page has been refreshed.\n",
       "if (IPython.notebook.kernel != null) {\n",
       "    IPython.notebook.kernel.comm_manager.register_target('matplotlib', mpl.mpl_figure_comm);\n",
       "}\n"
      ],
      "text/plain": [
       "<IPython.core.display.Javascript object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/html": [
       "<img src=\"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAoAAAAHgCAYAAAA10dzkAAAgAElEQVR4Xu2dC5RV5XXH9wjDSxkQYhSVh28FXBGM4iNQXAmp2mVATWNk1QfLqIlpqpg2xkQN1OCLVcGk8YHGgESM9UHUCom0UYigXVaLK+FlQbG8bFVghCDDa7r2wTOZuXPv3DN3D97Zd/8OizUw93z37v3//vvbvzmvqaqvr68XNhRAARRAARRAARRAgTAKVAGAYeaaRFEABVAABVAABVAgUQAAxAgogAIogAIogAIoEEwBADDYhJMuCqAACqAACqAACgCAeAAFUAAFUAAFUAAFgikAAAabcNJFARRAARRAARRAAQAQD6AACqAACqAACqBAMAUAwGATTroogAIogAIogAIoAADiARRAARRAARRAARQIpgAAGGzCSRcFUAAFUAAFUAAFAEA8gAIogAIogAIogALBFAAAg0046aIACqAACqAACqAAAIgHUAAFUAAFUAAFUCCYAgBgsAknXRRAARRAARRAARQAAPEACqAACqAACqAACgRTAAAMNuGkiwIogAIogAIogAIAIB5AARRAARRAARRAgWAKAIDBJpx0UQAFUAAFUAAFUAAAxAMogAIogAIogAIoEEwBADDYhJMuCqAACqAACqAACgCAeAAFUAAFUAAFUAAFgikAAAabcNJFARRAARRAARRAAQAQD6AACqAACqAACqBAMAUAwGATTroogAIogAIogAIoAADiARRAARRAARRAARQIpgAAGGzCSRcFUAAFUAAFUAAFAEA8gAIogAIogAIogALBFAAAg0046aIACqAACqAACqAAAIgHUAAFUAAFUAAFUCCYAgBgsAknXRRAARRAARRAARQAAPEACqAACqAACqAACgRTAAAMNuGkiwIogAIogAIogAIAIB5AARRAARRAARRAgWAKAIDBJpx0UQAFUAAFUAAFUAAAxAMogAIogAIogAIoEEwBADDYhJMuCqAACqAACqAACgCAeAAFUAAFUAAFUAAFgikAAAabcNJFARRAARRAARRAAQAQD6AACqAACqAACqBAMAUAwGATTroogAIogAIogAIoAADiARRAARRAARRAARQIpgAAGGzCSRcFUAAFUAAFUAAFAEA8gAIogAIogAIogALBFAAAg0046aIACqAACqAACqAAAIgHUAAFUAAFUAAFUCCYAgBgsAknXRRAARRAARRAARQAAPEACqAACqAACqAACgRTAAAMNuGkiwIogAIogAIogAIAIB5AARRAARRAARRAgWAKAIDBJpx0UQAFUAAFUAAFUAAAxAMogAIogAIogAIoEEwBADDYhJMuCqAACqAACqAACgCAeAAFUAAFUAAFUAAFgikAAAabcNJFARRAARRAARRAAQAQD6AACqAACqAACqBAMAUAwGATTroogAIogAIogAIoAADiARRAARRAARRAARQIpgAAaJjwPXv2yPr166V79+5SVVVleCeGogAKoAAKoAAKfFoK1NfXy5YtW+TQQw+V/fbb79P62Hb1OQCgYTrWrl0rffv2NbwDQ1EABVAABVAABcqlwJo1a+Twww8v18eX9XMrAgBvv/12efrpp2X58uXStWtXOeOMM+TOO++U4447rkVxn3rqKbn55ptl1apVctRRR8mkSZPk/PPPzzwhtbW10rNnT3nnnXekV69emcd533Hnzp3ywgsvyJe//GWprq72nk7m+Mmb+c5sFsc74nN87ti+mUPfuHGjHHHEEbJ582bp0aNH5nGVtGNFAODZZ58tX//61+WUU06RXbt2yQ9/+EP5wx/+IEuXLpX9998/73y98sorMnz4cLn11lsT6Js9e7bccsst8vLLL8uwYcMyzfFHH32UGOeDDz6Q3r17ZxpTCTtpg5gzZ46ce+654QCQvGMBP/PNfFfCmt1SDlHX8w8//FA+85nPiB7IqampqfRpzptfRQBgbmbvv/++fPazn5X58+fLiBEj8iZ+0UUXiQLc3LlzG15XkDzwwAPlsccey2QGABAAzGQU5ztFbRDkTX07L91M4Uf1OQAoUpEAuHLlSjnmmGOSo4CDBw/OWwT9+vWT8ePHJ3/TbcqUKTJ16lR59913MxUOAEiDyGQU5ztFbRDkTX07L91M4Uf1OQBYgQCod/aMHj1aNm3aJL///e8LFkCnTp1k+vTpMnbs2IZ9Zs2aJePGjZO6urq84/T7jV9TANSbQDZs2BDuFPC8efNk1KhR4U4Bk3esU4LMN/OdiaIc76QAGNHnCoB9+vThFLBj7zYL/dvf/rY8//zzybV8Ld3ZowA4Y8YMufjiixve49FHH5UrrrhCtm/fnleSCRMmyMSJE5u9puDYrVu3SpKRXFAABVAABVCgYhXYtm1bcgCIawArZIq/853vyK9//WtZsGBBcndPS1spp4A5ArhX0ag/MZI3R3wrZKlsMQ18js8j+JwjgBVyClhP+yr86Z28L730UnL9X7FNbwLRh0DqXX7pds455ySPdeEmkJbVi3rNCHlzTVixdaUSXsfn+LwSfFwsB64BrBAAvOaaa0RPwz7zzDNNnv2nj2jR5wLqdumll8phhx0m+sxA3RYtWpTcIazP/tNrBnXsTTfdxGNgilXNJ0cAeTxGrGujmG/mO8PS4HoXwDcW+AKAFQKAhX4N2y9+8Qu5/PLLk0Vp5MiRMmDAgOTGj3R78sknE+h7++23Gx4EfcEFF2RexLgLONaCQYNgvjMvDo53xOf43LF9M4cOAFYIAGae8TbeEQBkoWxjS7XLtwMI8Hm7NGYbB4XPY/kcAAQATUsIABhrwaBBMN+mBcPJYHyOz51Y1RQmAAgAmgwEALJQmgzkZDBAgM+dWNUUJj6P5XMAEAA0LRgAYKwFgwbBfJsWDCeD8Tk+d2JVU5gAIABoMhAAyEJpMpCTwQABPndiVVOY+DyWzwFAANC0YACAsRYMGgTzbVownAzG5/jciVVNYQKAAKDJQAAgC6XJQE4GAwT43IlVTWHi81g+BwABQNOCAQDGWjBoEMy3acFwMhif43MnVjWFCQACgCYDAYAslCYDORkMEOBzJ1Y1hYnPY/kcAAQATQsGABhrwaBBMN+mBcPJYHyOz51Y1RQmAAgAmgwEALJQmgzkZDBAgM+dWNUUJj6P5XMAEAA0LRgAYKwFgwbBfJsWDCeD8Tk+d2JVU5gAIABoMhAAyEJpMpCTwQABPndiVVOY+DyWzwFAANC0YACAsRYMGgTzbVownAzG5/jciVVNYQKAAKDJQAAgC6XJQE4GAwT43IlVTWHi81g+BwABQNOCAQDGWjBoEMy3acFwMhif43MnVjWFCQACgCYDAYAslCYDORkMEOBzJ1Y1hYnPY/kcAAQATQsGABhrwaBBMN+mBcPJYHyOz51Y1RQmAAgAmgwEALJQmgzkZDBAgM+dWNUUJj6P5XMAEAA0LRgAYKwFgwbBfJsWDCeD8Tk+d2JVU5gAIABoMhAAyEJpMpCTwQABPndiVVOY+DyWzwFAANC0YACAsRYMGgTzbVownAzG5/jciVVNYQKAAKDJQAAgC6XJQE4GAwT43IlVTWHi81g+BwABQNOCAQDGWjBoEMy3acFwMhif43MnVjWFCQACgCYDAYAslCYDORkMEOBzJ1Y1hYnPY/kcAAQATQsGABhrwaBBMN+mBcPJYHyOz51Y1RQmAAgAmgwEALJQmgzkZDBAgM+dWNUUJj6P5XMAEAA0LRgAYKwFgwbBfJsWDCeD8Tk+d2JVU5gAIABoMhAAyEJpMpCTwQABPndiVVOY+DyWzwFAANC0YACAsRYMGgTzbVownAzG5/jciVVNYQKAAKDJQAAgC6XJQE4GAwT43IlVTWHi81g+BwABQNOCAQDGWjBoEMy3acFwMhif43MnVjWFCQACgCYDAYAslCYDORkMEOBzJ1Y1hYnPY/kcAAQATQsGABhrwaBBMN+mBcPJYHyOz51Y1RQmAAgAmgwEALJQmgzkZDBAgM+dWNUUJj6P5XMAEAA0LRgAYKwFgwbBfJsWDCeD8Tk+d2JVU5gAIABoMhAAyEJpMpCTwQABPndiVVOY+DyWzwFAANC0YACAsRYMGgTzbVownAzG5/jciVVNYQKAAKDJQAAgC6XJQE4GAwT43IlVTWHi81g+BwABQNOCAQDGWjBoEMy3acFwMhif43MnVjWFCQACgCYDAYAslCYDORkMEOBzJ1Y1hYnPY/kcAAQATQsGABhrwaBBMN+mBcPJYHyOz51Y1RQmAAgAmgwEALJQmgzkZDBAgM+dWNUUJj6P5XMAEAA0LRgAYKwFgwbBfJsWDCeD8Tk+d2JVU5gAIABoMhAAyEJpMpCTwQABPndiVVOY+DyWzwFAANC0YACAsRYMGgTzbVownAzG5/jciVVNYQKAAKDJQAAgC6XJQE4GAwT43IlVTWHi81g+BwABQNOCAQDGWjBoEMy3acFwMhif43MnVjWFCQACgCYDAYAslCYDORkMEOBzJ1Y1hYnPY/kcAAQATQsGABhrwaBBMN+mBcPJYHyOz51Y1RQmAAgAmgwEALJQmgzkZDBAgM+dWNUUJj6P5XMAsIIAcMGCBTJ58mR5/fXXZcOGDTJ79mwZM2ZMwQXhpZdekrPOOqvZ68uWLZPjjz8+00ICAMZaMGgQzHemhcH5Tvgcnzu3cKbwAcAKAsC5c+fKwoULZejQoXLhhRdmBsAVK1ZITU1Ng2EOOugg6dChQyYDAYAslJmM4nwngACfO7dwpvDxeSyfA4AVBICNK7yqqiozAG7atEl69uyZaYHI3QkAjLVg0CCY75IWCmeD8Dk+d2bZksIFAAFAGTBggGzfvl0GDhwoN910U97Twqm76urqRP+mmwJg3759k1POvXv3LsmEHgdpg5g3b56MGjVKqqurPaZQUszkzXyXZBxng/A5Pndm2ZLCVQDs06eP1NbWNjkLWNKbOR1UVV9fX+809oJhZzkCqKd+9brBk08+OYG6mTNnyv333y96beCIESPyvveECRNk4sSJzV6bNWuWdOvWrdJkJB8UQAEUQAEUqEgFtm3bJmPHjgUAK212swBgvpzPO+880bHPPvtsXkk4ArhXFo4QcISg0taMfPngc3yOzytXAY4ABj4FnM/WkyZNkl/+8peidwJn2bgGkGtlsvjE+z5cE4bPvXs4S/z4PJbPuQYQAGyyLnz1q1+VjRs3yu9+97ss64UAgLEWDBoE851pYXC+Ez7H584tnCl8ALCCAHDr1q2ycuXKZOKHDBkid999d3JDR69evaRfv35y4403yrp16+SRRx5J9pk6dWpyA8igQYNkx44dyZG/O+64Q5566im54IILMhkIAGShzGQU5zsBBPjcuYUzhY/PY/kcAKwgACz0YOfLLrtMpk+fLpdffrmsXr06uclDt7vuukumTZuWQGHXrl0TEFRIPPfcczMtFroTABhrwaBBMN+ZFwfHO+JzfO7YvplDBwArCAAzz3ob7ggAslC2oZ3a7VsBBPi83ZqzDQPD57F8DgACgKblAwCMtWDQIJhv04LhZDA+x+dOrGoKEwAEAE0GAgBZKE0GcjIYIMDnTqxqChOfx/I5AAgAmhYMADDWgkGDYL5NC4aTwfgcnzuxqilMABAANBkIAGShNBnIyWCAAJ87saopTHwey+cAIABoWjAAwFgLBg2C+TYtGE4G43N87sSqpjABQADQZCAAkIXSZCAngwECfO7EqqYw8XksnwOAAKBpwQAAYy0YNAjm27RgOBmMz/G5E6uawgQAAUCTgQBAFkqTgZwMBgjwuROrmsLE57F8DgACgKYFAwCMtWDQIJhv04LhZDA+x+dOrGoKEwAEAE0GAgBZKE0GcjIYIMDnTqxqChOfx/I5AAgAmhYMADDWgkGDYL5NC4aTwfgcnzuxqilMABAANBkIAGShNBnIyWCAAJ87saopTHwey+cAIABoWjAAwFgLBg2C+TYtGE4G43N87sSqpjABQADQZCAAkIXSZCAngwECfO7EqqYw8XksnwOAAKBpwQAAYy0YNAjm27RgOBmMz/G5E6uawgQAAUCTgQBAFkqTgZwMBgjwuROrmsLE57F8DgACgKYFAwCMtWDQIJhv04LhZDA+x+dOrGoKEwAEAE0GAgBZKE0GcjIYIMDnTqxqChOfx/I5AAgAmhYMADDWgkGDYL5NC4aTwfgcnzuxqilMABAANBkIAGShNBnIyWCAAJ87saopTHwey+cAIABoWjAAwFgLBg2C+TYtGE4G43N87sSqpjABQADQZCAAkIXSZCAngwECfO7EqqYw8XksnwOAAKBpwQAAYy0YNAjm27RgOBmMz/G5E6uawgQAAUCTgQBAFkqTgZwMBgjwuROrmsLE57F8DgACgKYFAwCMtWDQIJhv04LhZDA+x+dOrGoKEwAEAE0GAgBZKE0GcjIYIMDnTqxqChOfx/I5AAgAmhYMADDWgkGDYL5NC4aTwfgcnzuxqilMABAANBkIAGShNBnIyWCAAJ87saopTHwey+cAIABoWjAAwFgLBg2C+TYtGE4G43N87sSqpjABQADQZCAAkIXSZCAngwECfO7EqqYw8XksnwOAAKBpwQAAYy0YNAjm27RgOBmMz/G5E6uawgQAAUCTgQBAFkqTgZwMBgjwuROrmsLE57F8DgACgKYFAwCMtWDQIJhv04LhZDA+x+dOrGoKEwAEAE0GAgBZKE0GcjIYIMDnTqxqChOfx/I5AAgAmhYMADDWgkGDYL5NC4aTwfgcnzuxqilMABAANBkIAGShNBnIyWCAAJ87saopTHwey+cAIABoWjAAwFgLBg2C+TYtGE4G43N87sSqpjABQADQZCAAkIXSZCAngwECfO7EqqYw8XksnwOAAKBpwQAAYy0YNAjm27RgOBmMz/G5E6uawgQAAUCTgQBAFkqTgZwMBgjwuROrmsLE57F8DgACgKYFAwCMtWDQIJhv04LhZDA+x+dOrGoKEwAEAE0GAgBZKE0GcjIYIMDnTqxqChOfx/I5AAgAmhYMADDWgkGDYL5NC4aTwfgcnzuxqilMABAANBkIAGShNBnIyWCAAJ87saopTHwey+cAIABoWjAAwFgLBg2C+TYtGE4G43N87sSqpjABQADQZCAAkIXSZCAngwECfO7EqqYw8XksnwOAAKBpwQAAYy0YNAjm27RgOBmMz/G5E6uawgQAAUCTgQBAFkqTgZwMBgjwuROrmsLE57F8DgACgKYFAwCMtWDQIJhv04LhZDA+x+dOrGoKEwCsIABcsGCBTJ48WV5//XXZsGGDzJ49W8aMGdOiQebPny/XX3+9LFmyRA499FD53ve+J9/85jczmwoAZKHMbBbHOwIE+NyxfTOHjs9j+RwArCAAnDt3rixcuFCGDh0qF154YVEAfOedd2Tw4MFy5ZVXytVXX52Mveaaa+Sxxx5LxmfZAMBYCwYNgvnOsi543wef43PvHs4SPwBYQQDYeMKrqqqKAuANN9wgzz77rCxbtqxhqB79e/PNN+WVV17J4h8BAFkoMxnF+U4AAT53buFM4ePzWD4HAAMD4IgRI2TIkCFyzz33NCwOetr4a1/7mmzbtk2qq6uLLhoAYKwFgwbBfBddFCpgB3yOzyvAxkVTAAADA+Cxxx4rl19+ufzgBz9oMMqiRYvkzDPPlPXr10ufPn2aGaiurk70b7opAPbt2ze55rB3795FDVcpO2iDmDdvnowaNSoTKJO3bwWYb3zu28HZosfnsXyuAKh9vra2VmpqarKZpML2qqqvr6+vsJwkyylgBcBx48bJjTfe2JC+Xgf4hS98IQG6Qw45pJksEyZMkIkTJzb7/qxZs6Rbt26VJiP5oAAKoAAKoEBFKqBn+saOHQsAVtrsZgHAUk4BcwRwr1P4STnWT8rMN/NdaT0iXz74PJbPOQIY+BSw3gTy3HPPydKlSxvWgm9961uyePFibgIpstpzjRDXCEUBgjlz5si55zLfzHflKhB1PecawAoCwK1bt8rKlSuTKtWbO+6++24566yzpFevXtKvX7/kVO+6devkkUceSfZJHwOjj4DRR8Honb96FzCPgSm+0EVdMMgbECpeHf73wOf43L+Li2cAAFYQAL700ksJ8OVul112mUyfPj254WP16tWi+6WbPgh6/PjxDQ+C1qOCPAi6eOHQIGgQxV3ifw98js/9u7h4BlF9DgBWEAAWt3nb78FjYGgQbe+q9veOURsEeVPf7a8a2z6iqD4HAAFAUzUBgDQIk4GcDI7aIMib+nZSoqYwo/ocAAQATYUDANIgTAZyMjhqgyBv6ttJiZrCjOpzABAANBUOAEiDMBnIyeCoDYK8qW8nJWoKM6rPAUAA0FQ4ACANwmQgJ4OjNgjypr6dlKgpzKg+BwABQFPhAIA0CJOBnAyO2iDIm/p2UqKmMKP6HAAEAE2FAwDSIEwGcjI4aoMgb+rbSYmawozqcwAQADQVDgBIgzAZyMngqA2CvKlvJyVqCjOqzwFAANBUOAAgDcJkICeDozYI8qa+nZSoKcyoPgcAAUBT4QCANAiTgZwMjtogyJv6dlKipjCj+hwABABNhQMA0iBMBnIyOGqDIG/q20mJmsKM6nMAEAA0FQ4ASIMwGcjJ4KgNgrypbyclagozqs8BQADQVDgAIA3CZCAng6M2CPKmvp2UqCnMqD4HAAFAU+EAgDQIk4GcDI7aIMib+nZSoqYwo/ocAAQATYUDANIgTAZyMjhqgyBv6ttJiZrCjOpzABAANBUOAEiDMBnIyeCoDYK8qW8nJWoKM6rPAUAA0FQ4ACANwmQgJ4OjNgjypr6dlKgpzKg+BwABQFPhAIA0CJOBnAyO2iDIm/p2UqKmMKP6HAAEAE2FAwDSIEwGcjI4aoMgb+rbSYmawozqcwAQADQVDgBIgzAZyMngqA2CvKlvJyVqCjOqzwFAANBUOAAgDcJkICeDozYI8qa+nZSoKcyoPgcAAUBT4QCANAiTgZwMjtogyJv6dlKipjCj+hwABABNhQMA0iBMBnIyOGqDIG/q20mJmsKM6nMAEAA0FQ4ASIMwGcjJ4KgNgrypbyclagozqs8BQADQVDgAIA3CZCAng6M2CPKmvp2UqCnMqD4HAAFAU+EAgDQIk4GcDI7aIMib+nZSoqYwo/ocAAQATYUDANIgTAZyMjhqgyBv6ttJiZrCjOpzABAANBUOAEiDMBnIyeCoDYK8qW8nJWoKM6rPAUAA0FQ4ACANwmQgJ4OjNgjypr6dlKgpzKg+BwABQFPhAIA0CJOBnAyO2iDIm/p2UqKmMKP6HAAEAE2FAwDSIEwGcjI4aoMgb+rbSYmawozqcwAQADQVDgBIgzAZyMngqA2CvKlvJyVqCjOqzwFAANBUOAAgDcJkICeDozYI8qa+nZSoKcyoPgcAAUBT4QCANAiTgZwMjtogyJv6dlKipjCj+hwABABNhQMA0iBMBnIyOGqDIG/q20mJmsKM6nMAEAA0FQ4ASIMwGcjJ4KgNgrypbyclagozqs8BQADQVDgAIA3CZCAng6M2CPKmvp2UqCnMqD4HAAFAU+EAgDQIk4GcDI7aIMib+nZSoqYwo/ocAAQATYUDANIgTAZyMjhqgyBv6ttJiZrCjOpzABAANBUOAEiDMBnIyeCoDYK8qW8nJWoKM6rPAUAA0FQ4ACANwmQgJ4OjNgjypr6dlKgpzKg+BwABQFPhAIA0CJOBnAyO2iDIm/p2UqKmMKP6HAAEAE2FAwDSIEwGcjI4aoMgb+rbSYmawozqcwAQADQVDgBIgzAZyMngqA2CvKlvJyVqCjOqzwFAANBUOAAgDcJkICeDozYI8qa+nZSoKcyoPgcAAUBT4QCANAiTgZwMjtogyJv6dlKipjCj+hwABABNhQMA0iBMBnIyOGqDIG/q20mJmsKM6nMAEAA0FQ4ASIMwGcjJ4KgNgrypbyclagozqs8BQADQVDgAIA3CZCAng6M2CPKmvp2UqCnMqD4HAAFAU+EAgDQIk4GcDI7aIMib+nZSoqYwo/ocAAQATYUDANIgTAZyMjhqgyBv6ttJiZrCjOpzALDCAPDee++VyZMny4YNG2TQoEEydepUGT58eN7imD59uowbN67Zax9//LF06dIlU0EBgDSITEZxvlPUBkHe1Lfz0s0UflSfA4AVBICPP/64XHLJJaIQeOaZZ8oDDzwgDz30kCxdulT69evXrBAUAK+99lpZsWJFk9cOOeSQTEWjOwGANIjMZnG8Y9QGQd7Ut+OyzRx6VJ8DgBUEgMOGDZOhQ4fKfffd12D8E044QcaMGSO33357XgC87rrrZPPmzZkLJXdHAJAGUbJ5HA2M2iDIm/p2VKYlhxrV5wBghQDgjh07pFu3bvLEE0/I+eef31AIeoRv8eLFMn/+/LwA+I1vfEMOO+ww2b17t5x00kly6623ypAhQwoWUl1dnejfdFMA7Nu3b3LKuXfv3iUXoLeBumDMmzdPRo0aJdXV1d7CLzle8ma+SzaPo4H4HJ87smvJoSoA9unTR2pra6Wmpqbk9/E8sKq+vr7ecwIa+/r16xOQW7hwoZxxxhkN6dx2220yY8aMZqd5dYdXX31VVq5cKSeeeGJyKveee+6ROXPmyJtvvinHHHNMXkkmTJggEydObPbarFmzEgBlQwEUQAEUQAEUaP8KbNu2TcaOHQsAtv+pajnCFAAXLVokp59+esPOkyZNkpkzZ8ry5cuLprhnz57kFPKIESPkJz/5Sd79OQK4VxaOEHCEoGhBVcAO+ByfV4CNi6YQ1eccAQx8CjhfVVx55ZWydu1amTt3btGi0R24BpBrhDIZxflOUa8RIm/q23npZgo/qs+5BrBCAFBdrjeBnHzyycldwOk2cOBAGT16dN6bQHIrQ8+En3rqqckp4YcffjhT4QCANIhMRnG+U9QGQd7Ut/PSzRR+VJ8DgBUEgOljYO6///7kNPC0adPkwQcflCVLlkj//v3l0ksvTa4TTO8I1mv5TjvttAXQ11wAACAASURBVOR6PwU5Pe2rp4v1OkIFwSwbAEiDyOIT7/tEbRDkTX17r90s8Uf1OQBYQQCoRtejf3fddVdyV+7gwYNlypQpyTV9uo0cOVIGDBgg+vw/3caPHy9PP/20vPfee9KjR4/k7l+9yaPxNYTFigcApEEU80glvB61QZA39V0J9Vssh6g+BwArDACLGb2tXwcAaRBt7an2+H5RGwR5U9/tsR7bOqaoPgcAAUBTLQGANAiTgZwMjtogyJv6dlKipjCj+hwABABNhQMA0iBMBnIyOGqDIG/q20mJmsKM6nMAEAA0FQ4ASIMwGcjJ4KgNgrypbyclagozqs8BQADQVDgAIA3CZCAng6M2CPKmvp2UqCnMqD4HAAFAU+EAgDQIk4GcDI7aIMib+nZSoqYwo/ocAAQATYUDANIgTAZyMjhqgyBv6ttJiZrCjOpzABAANBUOAEiDMBnIyeCoDYK8qW8nJWoKM6rPAUAA0FQ4ACANwmQgJ4OjNgjypr6dlKgpzKg+BwABQFPhAIA0CJOBnAyO2iDIm/p2UqKmMKP6HAAEAE2FAwDSIEwGcjI4aoMgb+rbSYmawozqcwAQADQVDgBIgzAZyMngqA2CvKlvJyVqCjOqzwFAANBUOAAgDcJkICeDozYI8qa+nZSoKcyoPgcAAUBT4QCANAiTgZwMjtogyJv6dlKipjCj+hwABABNhQMA0iBMBnIyOGqDIG/q20mJmsKM6nMAEAA0FQ4ASIMwGcjJ4KgNgrypbyclagozqs8BQADQVDgAIA3CZCAng6M2CPKmvp2UqCnMqD4HAAFAU+EAgDQIk4GcDI7aIMib+nZSoqYwo/ocAAQATYUDANIgTAZyMjhqgyBv6ttJiZrCjOpzABAANBUOAEiDMBnIyeCoDYK8qW8nJWoKM6rPAUAA0FQ4ACANwmQgJ4OjNgjypr6dlKgpzKg+BwABQFPhAIA0CJOBnAyO2iDIm/p2UqKmMKP6HAAEAE2FAwDSIEwGcjI4aoMgb+rbSYmawozqcwAQADQVDgBIgzAZyMngqA2CvKlvJyVqCjOqzwFAANBUOAAgDcJkICeDozYI8qa+nZSoKcyoPgcAAUBT4QCANAiTgZwMjtogyJv6dlKipjCj+hwABABNhQMA2htEfX297K7fLXW76mTH7h1St/uTr3n+39JrOlb/6nvt2rNLdu/Z3eTfyffqdyffT/9d8Ht7dsue+j3NvKGxvv/++3LQQQdJVVVVUe/sV7WfdNyvo3TYr8Per1WffNX/V/35+81e+2TfTh06SeeOnaVzh87Nvnbp2KXga+mYrtVdk7FZYm0pmagNgrzt9V20SNrRDsx3rPkGAAFA0/LjCQAVXhR4FKIUthrDVOPvFXtNIWvbjm2y7K1lcviAw2Xnnp1/fq8W4C0X7hp/fr3Um+aBwYUVqJIq6VbdTfbvtH/yNfl3daN/f/L9Jt/Leb1rh66y5I0l8qURX5Le+/eW7p27S/dO3UUh1AqX7XnuAIJYQMB8x5pvABAANPWfFADvnHendO3eVRRk9MiRwlbyNef/+b6X7pMemVKgUlDbufuTr+n/G32/4D6NxiTAlQN6XkCr8RGv9ChY8rVDZyn0f32ter/qhiNuTY62fXJELd+RuEL76dE7hafG267du+TNN9+Uz33uc9KxQ8cWvZM7r42PTDY+ApkerWxyVPKTo5QpNKeAvn3X9iYAn85vvq+fxlzrkUsFwRQIazrXNP33J68l38/Zr2eXntKjSw/p0blH8lXfq71tAEEsIGC+Y803AAgAmnpOCoDyfRHpYnqrT32wAk6+U4stnXZMIUxBa/2a9XLc0cclR5Qaj2kJ0LLAnIJAez2q5KVB6A8g+kOCAuO2nduSv3/a8ae9X3fu/dr4e42/n+y368+v6/4fbf9I3tv0ntRX18tHdR8l79HWmx6BVBBMwPATKNSvjf9f6N+6n0KoerotNy/z3ZY563uRdywQijrfACAAaFo7UwA87+HzpMsBXZIGlBw5qqra+1X2fm387/S13O/pkajqDnuPYKVHslrz/9x9c2Er91oyPRpW6hZ1wSDvvY1Rj1YnYFj3kWyp2yJbdmxp8m/9XvLaji3J6x/taLqfvla7vVY2b98sH+/6uFQbNhmntaZHGnMhsgEaP4HKfBCZjlMIbfzDB/MNCLWJOdv5m0T1OQAIAJpK09M1gKZEcwZHXTDIu+2BQC91qK2rbQDCFv9dYD892tkWm/5QpqeqFSIVCms61cj22u1ydN+j9x6JTL+voNl57z7p99L/6/e8XxuJz9ve523hz331HlHnGwAEAE01BQCyUJoM5GRwe24QerpbT3UrOOoRRT2yWAwic/fTI5J6jWZbbXoEPx8cJrDYqSZ5TU9Z6xHHAzodkPzVm3TSfyf/b/SaHt3/NLf2PN/7UgfyjrWeA4AAoGk9AQBjLRg0iMqcb4VIvQ5SQTA5PV1Xm3z98E8fysLXFsqA4wbI1l1b9762vTY5pd3w70b76+nufXHzjV7OkQUWu3bsKvron8Zf9Whk7vf0/8n3c/ZPr6HE55Xp80LNrlLnO70ZM735Mjd/BcBDDz5UamtrpaamxsQCXgdX1atKbCUpAACyUJZkHGeDKrVBFJuG1uatjWbrjkagmAOUCTx+Aox6s83WnVuT/ZN/79j77+T/O/f+X+8C/zQ3Bc0UDOt31kuPA3okN3jl+6tHJZt8f78/75f7Wu7d+Xr9sV7znPVr+izNxmMa3/CTXreZ3rVf6v937dolC+YvkOEjhkvHjh2TJzlkfZpDS094yPdkiELfSz/z0/y6c9dOWf3uajm87+GiDz5IPlv/fJK/169Fa2e7iNwhAGBRodghrwIAIAAYoTRaC0KVokm581YAVDhMgbABEHOAMf2+3lDz8c6PZfvu7cnX9P8N39+1vdn32ur6yUqZc/IIpAAAKBwBNPgdAAQADfZxM7TcIFQuoSLkrXd06zWUjWHxo48/khcXvCinnHaK7KnakxyJbPxXb9zJ/V6zffb8eR99jmX6m3lK+dr4N/akv81HT7U3PnmVnnpPv1fq/3fs2CGdO+397TnNnuTQwveKPemh4NMhqvS4ZVVyNDR9YkTJX2XvEyda+1c1W7VylRx/7PFS3bG61eNb+3mf9v7p0zhy15GNH26UIw8/kiOA5VpgvX8uAAgAevdwlvgjgFA+Hcib+s5SH973iepzbgLhJhBT7QKANAiTgZwMjtogyJv6dlKipjCj+hwABABNhQMA0iBMBnIyOGqDIG/q20mJmsKM6nMAEAA0FQ4ASIMwGcjJ4KgNgrypbyclagozqs8BQADQVDgAIA3CZCAng6M2CPKmvp2UqCnMqD4HAAFAU+EAgDQIk4GcDI7aIMib+nZSoqYwo/ocAAQATYUDANIgTAZyMjhqgyBv6ttJiZrCjOpzABAANBUOAEiDMBnIyeCoDYK8qW8nJWoKM6rPAUAA0FQ4ACANwmQgJ4OjNgjypr6dlKgpzKg+BwABQFPhAIA0CJOBnAyO2iDIm/p2UqKmMKP6HAAEAE2FAwDSIEwGcjI4aoMgb+rbSYmawozqcwAQADQVDgBIgzAZyMngqA2CvKlvJyVqCjOqzwFAANBUOAAgDcJkICeDozYI8qa+nZSoKcyoPgcAAUBT4QCANAiTgZwMjtogyJv6dlKipjCj+hwArDAAvPfee2Xy5MmyYcMGGTRokEydOlWGDx9esDieeuopufnmm2XVqlVy1FFHyaRJk+T888/PXEwAIA0is1kc7xi1QZA39e24bDOHHtXnAGAFAeDjjz8ul1xyiSgEnnnmmfLAAw/IQw89JEuXLpV+/fo1K4ZXXnklgcNbb701gb7Zs2fLLbfcIi+//LIMGzYsU/EAgDSITEZxvlPUBkHe1Lfz0s0UflSfA4AVBIAKbUOHDpX77ruvwfQnnHCCjBkzRm6//fZmhXDRRReJAtzcuXMbXjv77LPlwAMPlMceeyxT4QCANIhMRnG+U9QGQd7Ut/PSzRR+VJ8DgBUCgDt27JBu3brJE0880eQU7rXXXiuLFy+W+fPnNysEPSo4fvz45G+6TZkyJTlt/O677+YtnLq6OtG/6aYA2Ldv3+SUc+/evTMVWyXspAvGvHnzZNSoUVJdXV0JKWXKgbyZ70xGcb4TPsfnzi2cKXwFwD59+khtba3U1NRkGlNpO1XV19fXe09q/fr1cthhh8nChQvljDPOaEjntttukxkzZsiKFSuapdipUyeZPn26jB07tuG1WbNmybhx45pAXuOBEyZMkIkTJzZ7Lx2nAMqGAiiAAiiAAijQ/hXYtm1b0v8BwPY/Vy1GmALgokWL5PTTT2/YV2/qmDlzpixfvjwvACocXnzxxQ2vPfroo3LFFVfI9u3b834eRwD3ysIRAo4QOF8yMoWPz/F5JqM43ymqzzkCyCngVp0Czq1zrgHkGiHna3+m8KNeI0Te1HemAnG+U1Sfcw1ghQCg1p/eBHLyyScndwGn28CBA2X06NEFbwLZsmWLzJkzp2H/c845R3r27MlNIEUWtKgLBnkDBM57fabw8Tk+z2QU5zsBgBUEgOljYO6///7kNPC0adPkwQcflCVLlkj//v3l0ksvTa4TTO8I1tPFI0aMSJ79p5D4zDPPyE033cRjYDIUNQ2CBpHBJu53wef43L2JMyQQ1ecAYAUBoPpcj/7dddddyV25gwcPFr2rVyFPt5EjR8qAAQOSGz/S7cknn0yg7+233254EPQFF1yQoWT27sIpYBpEZrM43jFqgyBv6ttx2WYOParPAcAKA8DMjm+jHQFAGkQbWaldv03UBkHe1He7Lsw2Ci6qzwFAANBUQgAgDcJkICeDozYI8qa+nZSoKcyoPgcAAUBT4QCANAiTgZwMjtogyJv6dlKipjCj+hwABABNhQMA0iBMBnIyOGqDIG/q20mJmsKM6nMAEAA0FQ4ASIMwGcjJ4KgNgrypbyclagozqs8BQADQVDgAIA3CZCAng6M2CPKmvp2UqCnMqD4HAAFAU+EAgDQIk4GcDI7aIMib+nZSoqYwo/ocAAQATYUDANIgTAZyMjhqgyBv6ttJiZrCjOpzABAANBUOAEiDMBnIyeCoDYK8qW8nJWoKM6rPAUAA0FQ4ACANwmQgJ4OjNgjypr6dlKgpzKg+BwABQFPhAIA0CJOBnAyO2iDIm/p2UqKmMKP6HAAEAE2FAwDSIEwGcjI4aoMgb+rbSYmawozqcwAQADQVDgBIgzAZyMngqA2CvKlvJyVqCjOqzwFAANBUOAAgDcJkICeDozYI8qa+nZSoKcyoPgcAAUBT4QCANAiTgZwMjtogyJv6dlKipjCj+hwABABNhQMA0iBMBnIyOGqDIG/q20mJmsKM6nMAEAA0FQ4ASIMwGcjJ4KgNgrypbyclagozqs8BQADQVDgAIA3CZCAng6M2CPKmvp2UqCnMqD4HAAFAU+EAgDQIk4GcDI7aIMib+nZSoqYwo/ocAAQATYUDANIgTAZyMjhqgyBv6ttJiZrCjOpzABAANBUOAEiDMBnIyeCoDYK8qW8nJWoKM6rPAUAA0FQ4ACANwmQgJ4OjNgjypr6dlKgpzKg+BwABQFPhAIA0CJOBnAyO2iDIm/p2UqKmMKP6HAAEAE2FAwDSIEwGcjI4aoMgb+rbSYmawozqcwAQADQVDgBIgzAZyMngqA2CvKlvJyVqCjOqzwFAANBUOAAgDcJkICeDozYI8qa+nZSoKcyoPgcAAUBT4QCANAiTgZwMjtogyJv6dlKipjCj+hwABABNhQMA0iBMBnIyOGqDIG/q20mJmsKM6nMAEAA0FQ4ASIMwGcjJ4KgNgrypbyclagozqs8BQADQVDgAIA3CZCAng6M2CPKmvp2UqCnMqD4HAAFAU+EAgDQIk4GcDI7aIMib+nZSoqYwo/ocAAQATYUDANIgTAZyMjhqgyBv6ttJiZrCjOpzABAANBUOAEiDMBnIyeCoDYK8qW8nJWoKM6rPAUAA0FQ4ACANwmQgJ4OjNgjypr6dlKgpzKg+BwABQFPhAIA0CJOBnAyO2iDIm/p2UqKmMKP6HAAEAE2FAwDSIEwGcjI4aoMgb+rbSYmawozqcwAQADQVDgBIgzAZyMngqA2CvKlvJyVqCjOqzwFAANBUOAAgDcJkICeDozYI8qa+nZSoKcyoPgcAAUBT4QCANAiTgZwMjtogyJv6dlKipjCj+hwABABNhQMA0iBMBnIyOGqDIG/q20mJmsKM6nMAEAA0FQ4ASIMwGcjJ4KgNgrypbyclagozqs8BQADQVDgAIA3CZCAng6M2CPKmvp2UqCnMqD4HAAFAU+EAgDQIk4GcDI7aIMib+nZSoqYwo/ocAAQATYUDANIgTAZyMjhqgyBv6ttJiZrCjOpzABAANBUOAEiDMBnIyeCoDYK8qW8nJWoKM6rPAUAA0FQ4ACANwmQgJ4OjNgjypr6dlKgpzKg+BwABQFPhAIA0CJOBnAyO2iDIm/p2UqKmMKP6HAAEAE2FAwDSIEwGcjI4aoMgb+rbSYmawozqcwAQADQVDgBIgzAZyMngqA2CvKlvJyVqCjOqzwFAANBUOAAgDcJkICeDozYI8qa+nZSoKcyoPgcAKwQAN23aJH/3d38nzz77bFIIX/nKV+SnP/2p9OzZs2BhjBw5UubPn9/k9Ysuukh+9atfZS4mAJAGkdksjneM2iDIm/p2XLaZQ4/qcwCwQgDwnHPOkbVr18q0adMS01911VUyYMAAee6551oEwGOPPVb+8R//sWGfrl27So8ePTIXDgBIg8hsFsc7Rm0Q5E19Oy7bzKFH9TkAWAEAuGzZMhk4cKC8+uqrMmzYsMT0+u/TTz9dli9fLscdd1zeQtAjgCeddJJMnTo1c6Hk7ggA0iBKNo+jgVEbBHlT347KtORQo/ocAKwAAHz44Yfl+uuvl82bNzcpAD39O2XKFBk3blxBAFyyZInU19fLwQcfLHoU8Uc/+pF07949cyEBgDSIzGZxvGPUBkHe1Lfjss0celSfA4AVAIC33XabTJ8+Xd56660mhtfTuwp/N954Y95CePDBB+WII46QQw45RP74xz8m+x199NEyb968goVTV1cn+jfdFAD79u0rGzZskN69e2cuOO876oKhOo0aNUqqq6u9p5M5fvJmvjObxfGO+ByfO7Zv5tAVAPv06SO1tbVSU1OTeVwl7VhVr4fA2uE2YcIEmThxYouRvfbaa/LCCy/IjBkzZMWKFU32PeaYY+SKK66Q73//+5mye/311+Xzn/+86NehQ4fmHVMoplmzZkm3bt0yfQ47oQAKoAAKoAAKlFeBbdu2ydixYwHA8k5D/k//4IMPRP+2tOmNHgpfpZwCzn1f5eDOnTvLzJkzRe8GzrdxBHCvKhwh4AhBe1wz2jomfI7P29pT7fH9ovqcI4AVcAo4vQnkP/7jP+TUU09N6kv/fdppp7V4E0huIepp4BNPPDF5NMyIESMy1SnXAHKNUCajON8p6jVC5E19Oy/dTOFH9TnXAFYAAKrD9QaO9evXywMPPJAYXh8D079//4bHwKxbt06++MUvyiOPPJJA4qpVq+TRRx+Vc889Vz7zmc/I0qVL5bvf/a7oY2D0tHKHDh0yFQ4ASIPIZBTnO0VtEORNfTsv3UzhR/U5AFghALhx48ZmD4L+53/+54YHQa9evTq54ePFF18UffzLmjVr5G/+5m+Smz+2bt2a3MjxV3/1V8ldwL169cpUNLoTAEiDyGwWxztGbRDkTX07LtvMoUf1OQBYIQCY2eltvCMASINoY0u1y7eL2iDIm/pulwXZxkFF9TkACACaSgkApEGYDORkcNQGQd7Ut5MSNYUZ1ecAIABoKhwAkAZhMpCTwVEbBHlT305K1BRmVJ8DgACgqXAAQBqEyUBOBkdtEORNfTspUVOYUX0OAAKApsIBAGkQJgM5GRy1QZA39e2kRE1hRvU5AAgAmgoHAKRBmAzkZHDUBkHe1LeTEjWFGdXnACAAaCocAJAGYTKQk8FRGwR5U99OStQUZlSfA4AAoKlwAEAahMlATgZHbRDkTX07KVFTmFF9DgACgKbCAQBpECYDORkctUGQN/XtpERNYUb1OQAIAJoKBwCkQZgM5GRw1AZB3tS3kxI1hRnV5wAgAGgqHACQBmEykJPBURsEeVPfTkrUFGZUnwOAAKCpcABAGoTJQE4GR20Q5E19OylRU5hRfQ4AAoCmwgEAaRAmAzkZHLVBkDf17aRETWFG9TkACACaCgcApEGYDORkcNQGQd7Ut5MSNYUZ1ecAIABoKhwAkAZhMpCTwVEbBHlT305K1BRmVJ8DgACgqXAAQBqEyUBOBkdtEORNfTspUVOYUX0OAAKApsIBAGkQJgM5GRy1QZA39e2kRE1hRvU5AAgAmgoHAKRBmAzkZHDUBkHe1LeTEjWFGdXnACAAaCocAJAGYTKQk8FRGwR5U99OStQUZlSfA4AAoKlwAEAahMlATgZHbRDkTX07KVFTmFF9DgACgKbCAQBpECYDORkctUGQN/XtpERNYUb1OQAIAJoKBwCkQZgM5GRw1AZB3tS3kxI1hRnV5wAgAGgqHACQBmEykJPBURsEeVPfTkrUFGZUnwOAAKCpcABAGoTJQE4GR20Q5E19OylRU5hRfQ4AAoCmwgEAaRAmAzkZHLVBkDf17aRETWFG9TkACACaCgcApEGYDORkcNQGQd7Ut5MSNYUZ1ecAIABoKhwAkAZhMpCTwVEbBHlT305K1BRmVJ8DgACgqXAAQBqEyUBOBkdtEORNfTspUVOYUX0OAAKApsIBAGkQJgM5GRy1QZA39e2kRE1hRvU5AAgAmgoHAKRBmAzkZHDUBkHe1LeTEjWFGdXnACAAaCocAJAGYTKQk8FRGwR5U99OStQUZlSfA4AAoKlwAEAahMlATgZHbRDkTX07KVFTmFF9DgACgKbCAQBpECYDORkctUGQN/XtpERNYUb1OQAIAJoKBwCkQZgM5GRw1AZB3tS3kxI1hRnV5wAgAGgqHACQBmEykJPBURsEeVPfTkrUFGZUnwOAAKCpcABAGoTJQE4GR20Q5E19OylRU5hRfQ4AAoCmwgEAaRAmAzkZHLVBkDf17aRETWFG9TkACACaCgcApEGYDORkcNQGQd7Ut5MSNYUZ1ecAIABoKhwAkAZhMpCTwVEbBHlT305K1BRmVJ8DgACgqXAAQBqEyUBOBkdtEORNfTspUVOYUX0OAAKApsIBAGkQJgM5GRy1QZA39e2kRE1hRvU5AAgAmgoHAKRBmAzkZHDUBkHe1LeTEjWFGdXnACAAaCocAJAGYTKQk8FRGwR5U99OStQUZlSfA4AAoKlwAEAahMlATgZHbRDkTX07KVFTmFF9DgACgKbCAQBpECYDORkctUGQN/XtpERNYUb1OQAIAJoKBwCkQZgM5GRw1AZB3tS3kxI1hRnV5wAgAGgqHACQBmEykJPBURsEeVPfTkrUFGZUnwOAAKCpcABAGoTJQE4GR20Q5E19OylRU5hRfQ4AAoCmwgEAaRAmAzkZHLVBkDf17aRETWFG9TkACACaCgcApEGYDORkcNQGQd7Ut5MSNYUZ1ecAIABoKhwAkAZhMpCTwVEbBHlT305K1BRmVJ8DgBUCgJMmTZLnn39eFi9eLJ06dZLNmzcXLYj6+nqZOHGiTJs2TTZt2iTDhg2Tn/3sZzJo0KCiY9MdAEAaRGazON4xaoMgb+rbcdlmDj2qzwHACgHAH/3oR9KzZ09Zu3at/PznP88EgHfeeacoOE6fPl2OPfZY+fGPfywLFiyQFStWSPfu3TMVDwBIg8hkFOc7RW0Q5E19Oy/dTOFH9TkAWCEAmLpcYe66664rCoB69O/QQw9N9r3hhhuS4XV1dXLwwQeLguHVV1+dqXAAQBpEJqM43ylqgyBv6tt56WYKP6rPAcCgAPj222/LUUcdJW+88YYMGTKkoUhGjx6dHEmcMWNG3sJRSNS/6VZbWyv9+vWTt956S3r16pWp2CphJ10wXnzxRTnrrLOkurq6ElLKlAN5M9+ZjOJ8J3yOz51bOFP4GzduTM7+6SVjPXr0yDSm0naqqtfDYRWyZT0CuGjRIjnzzDNl3bp1yZHAdLvqqqvk3Xffld/+9rd5FZkwYUJy3SAbCqAACqAACqCAfwVWrVolRx55pP9ESsig3QJgFth67bXX5POf/3xD2q0FwPXr10ufPn0axl955ZWyZs0a+c1vfpPpCKD+5NC/f3/5n//5n1A/Qeip7759+yZa1dTUlGA7n0PIm/n26dzWRY3P8XnrHONz7/QMnt4Eqmf+Im7tFgA/+OAD0b8tbQMGDJAuXbq0GgBLPQWcG0t6DaAaKRoI6SFz8o6xZOBz6juC0/E5Po/g88Y5tlsALGUish4BTG8CGT9+vHzve99LPmrHjh3y2c9+tqSbQAChUmbL3xgaBA3Cn2tbHzE+x+etd42/EVF9XnEAqKdg9YLOZ599ViZPniy///3vkxyPPvpoOeCAA5J/H3/88XL77bfL+eefn/xf7/bV///iF7+QY445Rm677TZ56aWXSnoMDADor/hLiTjqgkHeAEEp9eJtDD7H5948a423Io4AXn755Xnv3NU7VUeOHJloVFVVlcCe7qtb+iDoBx54oMmDoAcPHpxZU70jWCHyxhtvlM6dO2ce531H8ma+vXs4S/z4HJ9n8Yn3ffB5LJ9X3BFA7wVI/CiAAiiAAiiAAijwaSpQEUcAP03B+CwUQAEUQAEUQAEU8K4AAOh9BokfBVAABVAABVAABVqpAADYSsHYHQVQAAVQAAVQAAW8KwAAep9B4kcBFEABFEABFECBVioAABYR7N57700eLbNhwwYZNGiQTJ06VYYPH15w1FNPPSU333yz6K+X0d83PGnSpIZHz7Rybsqyu97V/PTTT8vy5cul/cGdrgAADO9JREFUa9eucsYZZySPzDnuuOMKxqPPXxw3blyz1z/++OMmD+ouS0IZPzTfb545+OCD5b333iv4DvPnz5frr79elixZkvxKQX2m5De/+c2Mn9g+dtOHqeuvP8zdrrnmGvnZz37W7Pte53rBggVJHb/++utJLc+ePVvGjBnTkF/6VIBp06Y1eSqA1nxLW2vXh0971lvKW3/n70033SRz5swRfTi+Ptz9S1/6ktxxxx1NfkVmbsyl1Ep7yltjyffkiGHDhsmrr77aYqjtfX0v5nN9Gka+7a677pJ/+Id/yPuah/nO0rf0bue///u/l8cee0y0N33xi18Urd/DDz+84JyXui582n4v9fMAwBaUe/zxx+WSSy5JTKK/O1gfGfPQQw/J0qVLpV+/fs1GvvLKKwkc3nrrrQn0aZO55ZZb5OWXXxZdXDxsZ599tnz961+XU045RXbt2iU//OEP5Q9/+EOS8/777583BYWCa6+9NnmGYuPtkEMO8ZByEqMuck8++aT827/9W0PMHTp0kIMOOihvDu+8847oI4P01wdeffXVsnDhQlFo0sXlwgsvdJP3+++/L7t3726I949//KOMGjVKGj9CqXEyXud67ty5yRwNHTo0mZ9cANQfcvSHNc1Pf0H8j3/8Y9Fmqp7u3r173vls7fpQDlO0lLc+v/SrX/1q4uHPfe5zCfhed911Sd3/53/+Z8FwW1sr7S3vFAD/93//N3k0WLp16tRJevXqVTBcD+t7MZ/n/kCr+19xxRWycuXKgr8P18N8Z+lb3/rWt+S5555Larx3797y3e9+N3l+sP5QqGt9vq2UdaEcfi/1MwHAFpRTaNOGcd999zXsdcIJJyRHDvQnjtztoosuEn2YqBZVuqkxDzzwwAQMPG4KCPobUvRo14gRIwoCoDYO/d3IXjdd5H7961/L4sWLM6Vwww03JA8eX7ZsWcP+evTvzTffFG0UXjedx3/913+V//7v/06enZm7Zf1tO+05f82rMQCmvxlIc9d51U2PFugRYG0ACvj5ttauD+XWJDfvfPHo71c/9dRTk6PC+X7I1TGtrZX2mLceAdT1Sms+6+Ztfc8y39rLtmzZIv/+7/9eUAZv862J5PYt/WFHf5ifOXOm6Dzqtn79+uR32usR8L/8y79sln+p60JWP7WH/QDAArOgvxquW7du8sQTTzQ5hatHuhQSFIhyN10w9dfL6d90mzJlSnLaON9ptvZggGIx6E+G+ptS9ChgoYdkKxR84xvfkMMOOyw5mnTSSSclR0GHDBlS7O3bzeu6yOkpQj0Npg/11uauvx3myCOPzBujwrDmd8899zS8rlDxta99TbZt2ybV1dXtJresgajn9VS2ntb+wQ9+kHdYJcx1bmMs5XeDl7I+ZJ2HfbVfFiDQI+Bf/vKXEzgq9PvNW1sr+yqfrO+bL28FQIU/PerXs2dP+Yu/+IvkCLD+sFto87a+F5tvPQKqpz9nzJghY8eObREAW7M2Zp2Xfblfbt/63e9+l5zy1SN+ekAm3fTIt0LwxIkTm4VTyrqwL3PaF+8NABZQVX86UKDR00Z6HVy6KRRoweSe7tTXdTHRBtm4mGbNmpVcH6dHFLxt+hPQ6NGjk1ND6a/Xy5eDXjejBXfiiScmR0AVivSnKj0apvDoYdOjtgpuevpPF0Y9BajXQer1fXq6IHfT/bSJNAalRYsWJZcKqHf69OnjIe0mMf7Lv/xL4l391YoKgpU617mNMZ23devWNcn7qquuSn5w++1vf9tMilLWh3IbohgQbN++Xb7whS8kvzbzl7/8ZcFwW1sr7TFvPX2vvya0f//+opdz6HXbeupbTwcW+q1O3tb3YvOt1/3p9Z7q5S5dulTMfOfrW4X6sP6wc8QRRySXd+VupawL5fZ6az8fACwCgGqC008/vWEv/SlRDyMrHORuukAoHF588cUNLz366KPJNRa6uHrbvv3tb8vzzz+fXMPY0oWyuXnt2bMnOXWuR8l+8pOfeEs7ifdPf/pTchOP3tihR8TyAaCCvf4awHTTHxa0gepNBp6uf0zj19Mg6mG9Tibr5nGuCwFgLrjrtXFr1qyR3/zmNwUBsDXrQ1ZN99V+LQGB3hDy13/91wn86+9EL3T0L19sxWplX+WT9X2LgZC+j9aswuCvfvUrueCCC/K+tbf1vVjeCvp6ve9Pf/rTrFJmWhtb9Wb7YOd8fasQAGr+us7ff//9zSIp9AN9S+vCPkhnn74lAFhA3lJO8Xg7RdCSs77zne8kp0j0Qnj9Cam1mxbJ2rVrm1wP2dr3KPf+ujgcffTRTa4BTWOqtFPAeqRLT3frHeB61Lc1m7e55hTwn+9+1nlW+NNLF/SUl54qy3fEu5gfWqqVYmP39evFQCj9fD1boZeypNeB5sblbX1vKW89o6NrmF7OpKdBW7u11/ku1Lc4BZx/hgHAFpyv14GdfPLJyV3A6TZw4MCkQRa6CUQvqNXTn+l2zjnnJNeYeLkJRA+faxHp9Wx6JKCUU7j6HnohuZ4Sfvjhh1u7trSL/fWUvf5kqKcB9U7u3E2bhB4p07uj003vMtMF1eNNIHpdl54G0SNeHTt2zDwHHue60E0geu2uHvHVTX8A1OvBit0E0pr1IbOo+2jHfECQwp/e9KN3fhe6672lkIrVyj5KJ/PbZgHADz/8MLnkRx8DdOmll+Z9b715wNP63lLeevmK3vHf0t3ehQRuj/NdrG+lN4HopQ36w45uetRXz2wVuwmktetCZmO2gx0BwBYmIX3Mgx4e1tPAujg8+OCDyXVherpAFwpdNFIY1EPG+lOVniZWSHzmmWeS52x5egyMPspED5dr7I2f/ac3R+hzAXXLzVsvoD3ttNMSWNRrAPW0r54m11OiCoIeNn0+1HnnnZfc+fh///d/yTWAeqOP3vyic62nevUasUceeSRJJ30MjN4hqkfAFPr0LmBvj4HRXPQ0rh7l1UsX9JqgxlulzPXWrVuT61R105t37r77bjnrrLOSx37onCvoaR3rY0HUx3qtr/4A1PgxMHoRuT7e6W//9m+T9ym2PrQH37eUt17nqY/EeeONN5I7v/Wu53RTXfSUp265eRerlfaet+amP/Bo7nqt7urVq5NrefX0t97Vnz72x+P6XsznOje6Rmve//RP/5T3uaUe5ztL39If0NXnep2+ekB9rODf+DEwelpc1wGtc92yrAvtwe+lxgAAFlFOj/7pxbL604LeBat39aaPQxk5cqTog3TVUOmmz5JT6EvvIFIYLHRNSamTti/HFXpQqDZG/alRt9y89SckPXWoz5hSUNQGqwts42sn92XMbfHe+uxDPd39wQcfJEdBFGj1TmY94qub5q6NQqEg3RQQNff0QdB6VNDbg6A1lxdeeCF5DILCjt7c0nirlLnWeVPgy90uu+yypH7TB77qUVC96UmP/uuDsBvf+a61rj5Qb6dbS+tDW/jS+h4t5a15FLq8o/FzIHPzLlYr1pjbYnxLeetjvfTOz//6r/9K7nZWGFJvaL3rY0HSzeP6XsznmpseyNBHHmlP0/U6d/M431n6ll6Hrw+71gMcjR8E3XjO9X0a97os60Jb+LVc7wEAlkt5PhcFUAAFUAAFUAAFyqQAAFgm4flYFEABFEABFEABFCiXAgBguZTnc1EABVAABVAABVCgTAoAgGUSno9FARRAARRAARRAgXIpAACWS3k+FwVQAAVQAAVQAAXKpAAAWCbh+VgUQAEUQAEUQAEUKJcCAGC5lOdzUQAFUAAFUAAFUKBMCgCAZRKej0UBFEABFEABFECBcikAAJZLeT4XBVAABVAABVAABcqkAABYJuH5WBRAARRAARRAARQolwIAYLmU53NRAAVQAAVQAAVQoEwKAIBlEp6PRQEUQAEUQAEUQIFyKQAAlkt5PhcFUAAFUAAFUAAFyqQAAFgm4flYFEABFEABFEABFCiXAgBguZTnc1EABVAABVAABVCgTAoAgGUSno9FARRAARRAARRAgXIpAACWS3k+FwVQAAVQAAVQAAXKpAAAWCbh+VgUQAEUQAEUQAEUKJcCAGC5lOdzUQAFUAAFUAAFUKBMCgCAZRKej0UBFEABFEABFECBcikAAJZLeT4XBVAABVAABVAABcqkAABYJuH5WBRAARRAARRAARQolwIAYLmU53NRAAVQAAVQAAVQoEwKAIBlEp6PRQEUQAEUQAEUQIFyKQAAlkt5PhcFUAAFUAAFUAAFyqQAAFgm4flYFEABFEABFEABFCiXAgBguZTnc1EABVAABVAABVCgTAoAgGUSno9FARRAARRAARRAgXIpAACWS3k+FwVQAAVQAAVQAAXKpAAAWCbh+VgUQAEUQAEUQAEUKJcCAGC5lOdzUQAFUAAFUAAFUKBMCgCAZRKej0UBFEABFEABFECBcinw/24UHWqO4XChAAAAAElFTkSuQmCC\" width=\"640\">"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "%matplotlib notebook\n",
    "fig, ax = plt.subplots()\n",
    "ax.set_xlim(0,a)\n",
    "maxv=-1e10\n",
    "minv=1e10\n",
    "minv = [min(minv,min(t)) for t in u]\n",
    "minv=min(minv)\n",
    "maxv = [max(maxv,max(t)) for t in u]\n",
    "maxv=max(maxv)\n",
    "ax.set_xlim(0,a)\n",
    "ax.set_ylim(minv-1,maxv+1)\n",
    "ax.grid(True)\n",
    "line, = plt.plot(x, u[0], color='g')\n",
    "def animate(n):\n",
    "    #line, = plt.plot(x, u[n], color='g')\n",
    "    line.set_data(x, u[n])\n",
    "    return line,\n",
    "anim = animation.FuncAnimation(fig, animate, frames=len(t), interval=100,repeat = True)\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.6.4"
  },
  "varInspector": {
   "cols": {
    "lenName": 16,
    "lenType": 16,
    "lenVar": 40
   },
   "kernels_config": {
    "python": {
     "delete_cmd_postfix": "",
     "delete_cmd_prefix": "del ",
     "library": "var_list.py",
     "varRefreshCmd": "print(var_dic_list())"
    },
    "r": {
     "delete_cmd_postfix": ") ",
     "delete_cmd_prefix": "rm(",
     "library": "var_list.r",
     "varRefreshCmd": "cat(var_dic_list()) "
    }
   },
   "types_to_exclude": [
    "module",
    "function",
    "builtin_function_or_method",
    "instance",
    "_Feature"
   ],
   "window_display": false
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}