Video5.4-MakingMovies.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import numpy as np\n",
    "# import matplotlib as mpl\n",
    "import matplotlib.pyplot as plt"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "%matplotlib notebook\n",
    "#magic calls weren't allowed in Thonny (but they are in j.notebook)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "from matplotlib.animation import FuncAnimation, FFMpegWriter"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "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 / mpl.ratio, fig.canvas.height / mpl.ratio);\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,iVBORw0KGgoAAAANSUhEUgAAAoAAAAHgCAYAAAA10dzkAAAAAXNSR0IArs4c6QAAQABJREFUeAHsnQd8lFXW/w+k90JIQgm9S+8dFKQormVVXFbAfUHX11XX9u4uq/uu4qvsa1v9rwvuK6KuomKvgIBK0dCbgKG3AAkppBLS+d9zZ+6YSSaTSTIzT/vdz2eYmee5z73nfO/w5Dz33nNOi8uiEAoIgAAIgAAIgAAIgIBlCLS0jKZQFARAAARAAARAAARAQBKAAYgfAgiAAAiAAAiAAAhYjAAMQIsNONQFARAAARAAARAAARiA+A2AAAiAAAiAAAiAgMUIwAC02IBDXRAAARAAARAAARCAAYjfAAiAAAiAAAiAAAhYjAAMQIsNONQFARAAARAAARAAARiA+A2AAAiAAAiAAAiAgMUIwAC02IBDXRAAARAAARAAARCAAYjfAAiAAAiAAAiAAAhYjAAMQIsNONQFARAAARAAARAAARiA+A2AAAiAAAiAAAiAgMUIwAC02IBDXRAAARAAARAAARCAAYjfAAiAAAiAAAiAAAhYjAAMQIsNONQFARAAARAAARAAARiA+A2AAAiAAAiAAAiAgMUIwAC02IBDXRAAARAAARAAARCAAYjfAAiAAAiAAAiAAAhYjAAMQIsNONQFARAAARAAARAAARiA+A2AAAiAAAiAAAiAgMUIwAC02IBDXRAAARAAARAAARCAAYjfAAiAAAiAAAiAAAhYjAAMQIsNONQFARAAARAAARAAARiA+A2AAAiAAAiAAAiAgMUIwAC02IBDXRAAARAAARAAARCAAYjfAAiAAAiAAAiAAAhYjAAMQIsNONQFARAAARAAARAAARiA+A2AAAiAAAiAAAiAgMUIwAC02IBDXRAAARAAARAAARCAAYjfAAiAAAiAAAiAAAhYjAAMQIsNONQFARAAARAAARAAARiA+A2AAAiAAAiAAAiAgMUIwAC02IBDXRAAARAAARAAARCAAYjfAAiAAAiAAAiAAAhYjAAMQIsNONQFARAAARAAARAAARiA+A2AAAiAAAiAAAiAgMUIwAC02IBDXRDQC4E33niDWrRoQSdPnvSaSOvXr5dt8rteSlZWFt1xxx2UkJBA4eHhNGrUKPrmm288Fu+jjz6iMWPGUHx8PMXGxtLw4cPprbfe8vh6VAQBEAABVwRgALqigmMgAAI+J3DttdfS5s2bqU2bNj7vS6sOysrKaNKkSdLge+mll+izzz6jpKQkmjZtGm3YsKFBsZYtW0Y333yzZLR8+XJ67733qGvXrjRnzhz6+9//3uD1qAACIAAC9RFocVmU+k7iOAiAAAgYiQDP/F155ZX03Xff0cSJEzUXffHixfS73/2OUlNT5cwfC1RZWUkDBgygyMhI2rp1q1sZx44dS2fOnKHjx49Ty5a253W+Zffp04eCg4Np7969bq/HSRAAARCojwBmAOsjg+MgAAI+JeBqCZiNtr59+9L27dtp3Lhxcsm0S5cu9Le//Y2qq6ud5Dl48KCcSeNlVV5evfvuu6moqMipjvqybt06ORMXHR0t2+Ql1ZrLsKWlpTRo0CDq1q0bFRQUqMsoMzOTkpOTpTFZVVXlOO7ph08++YR69uzpMP74usDAQLr99ttp27ZtdPbsWbdNBQUFSUNRGX9cmZfNWY/Q0FC31+IkCIAACLgjAAPQHR2cAwEQ8DsBNrp+/etfSyPp888/p+nTp9OCBQvo7bffdshy/vx5mjBhAu3fv594lo33xBUXF9O9997rqKM+8HVTpkyRRtObb75J77//vtxPN3XqVIcRyMYUH+f9ev/xH/8hL2WDk+XgGbd3332XAgIC5HHes8hGGO/ra6iwfP37969TTR07cOBAnXM1D9x3332UlpZGTz31FGVnZ1NOTg4999xztHPnTnrkkUdqVsVnEAABEGgUgcBG1UZlEAABEPAxgdzcXFq5cqV0duCuJk+eTLy0+84778i9b3yM97+xQbR79265nMrH2FBkQ+/06dP8VZaSkhL6/e9/TzNmzCCejVPlmmuuocGDB9Of//xnxzJs9+7daenSpTRz5kzi/XoXLlyQ/a5evdppnyIbf2wMKoNQtenqnXVh543aRR3j8+7KTTfdRB9//DHNnTuXHnvsMVk1LCyM2JC95ZZb3F2KcyAAAiDglgAMQLd4cBIEQMDfBHjJlT1daxaeMduzZ4/jEO/xu+KKKxzGnzoxa9YsWrt2rfoq996xIccGFO+9q1nYEeOZZ56hixcvUkREhDx16623SqPvv/7rv4iXfNlAvPrqq2teRh07dqzTllOFWl/YYKyvuDvH17DxycvFbOyxbLx8zLOiPPtYXl5Ov/nNb+prGsdBAARAwC0BGIBu8eAkCICAvwm0atWqTpchISF06dIlx3GeOevcubPju/rAxmPNwkvFXNiTtr7CBqIyALkOLwEvWbJEOlncf//99V3m0XHWxdUsH/fJRc0EumqMl55ZlvHjxxN7A6vCM6K8T5GXh9korCm7qoN3EAABEGiIAPYANkQI50EABHRHgA0r3itYu9Q+xs4hXP7xj39IxxJ2Lqn94rAsqvBs4OzZs6lHjx7ES63z589Xp5r03q9fP9q3b1+da9Uxdnipr7DxmpGRUWc2lOsPGzZMzlx6M4ZifXLgOAiAgDkJwAA057hCKxAwNQEO9cIOFLXDoPA+wZqFvX05ePJPP/1EQ4cOdfnicCqqsCcx7yHkfXevvfaaXG5tTry9G2+8kdhbuWa4F16KZseUESNGUNu2bVXXdd7j4uKkp++WLVvqnOP4iewZbOYYinWUxgEQAAGvEsASsFdxojEQAAF/EHjggQfksigHk/6f//kfGVyZAyWzsVWzcKw9nv3jPYC87MpLwYmJidKBhI1HdiTh5V4u7ADChtnrr78u9xfyHkP2Kv7jH/8oM3GofYmnTp2SwZi5TTYS3RVewv3nP/8p9/BxKBvum72WDx06RByapmbhgNEcHFrtVeRl73vuuYdeeOEF6fzCzinsePLpp59Kh5h58+a5XUKu2TY+gwAIgEBtApgBrE0E30EABHRPgPf6sbHEAZH/8z//UzpKcCiXl19+uY7s7ETBTiMcJua3v/2t9Cpmz+Bdu3bJ2IB8AS/J8n4/NurYwUIVDrnCDihsfOXn58vDvDePHUQ8iQvIRhzHG+QZS96zd91118ll3VWrVskwNqoffnfV5rPPPkuvvvqqDAXDerAcHD+Q9VSGa8028BkEQAAEPCWATCCekkI9EAABEAABEAABEDAJAcwAmmQgoQYIgAAIgAAIgAAIeEoABqCnpFAPBEAABEAABEAABExCAAagSQYSaoAACIAACIAACICApwRgAHpKCvVAAARAAARAAARAwCQEYACaZCChBgiAAAiAAAiAAAh4SgAGoKekUA8EQAAEQAAEQAAETEIABqBJBhJqgAAIgAAIgAAIgICnBJAJxFNSLupVV1fTuXPnKCoqilq0aOGiBg6BAAiAAAiAAAjojQAHdC8qKpLpGDmtohULDMBmjDobfykpKc1oAZeCAAiAAAiAAAhoRSA9PZ3at2+vVfea9gsDsBn4eeaPC/+AoqOjm9ESLgUBEAABEAABEPAXgcLCQjmBo/6O+6tfPfUDA7AZo6GWfdn4gwHYDJC4FARAAARAAAQ0IKD+jmvQteZdWnPhW3PsEAAEQAAEQAAEQAAEtCMAA1A79ugZBEAABEAABEAABDQhAANQE+zoFARAAARAAARAAAS0IwADUDv26BkEQAAEQAAEQAAENCEAA1AT7OgUBEAABEAABEAABLQjAANQO/boGQRAAARAAARAAAQ0IQADUBPs6BQEQAAEQAAEQAAEtCMAA1A79ugZBEAABEAABEAABDQhAANQE+zoFARAAARAAARAAAS0I2AIA3Djxo103XXXyaTNHLX7008/bZDYhg0baMiQIRQaGkpdunShV155pc41ixcvps6dO8s6XHfTpk116uAACIAACIAACIAACJiNgCEMwIsXL9KAAQPo5Zdf9oj/iRMn6JprrqFx48bR7t276c9//jPdf//99NFHHzmuX7FiBT3wwAP06KOPyjpcd/r06XT69GlHHXwAARAAARAAARAAATMSaHFZFCMpxjOAn3zyCd1www31iv3HP/6RPv/8c0pLS3PUufvuu2nv3r20efNmeWzEiBE0ePBgWrJkiaNO7969ZbuLFi1yHHP3gZNJx8TEUEFBAXIBuwOFcyAAAiAAAiCgIwL4+00UqKPx8JoobORNmTLFqb2pU6fSa6+9RhUVFcQ2786dO+lPf/qTUx2+JjU11elYzS9lZWXEL1X4B2TGcrGskraeyKWc4nLKuyheJRXyPS4imCb2bE1DOsZRUIAhJo/NODzQCQRMSYDvy3vPFNCW47lUXllNldWXqVq8+L1jq3Ca3jeZYsODTak7lAIBLQiY0gDMzMykpKQkJ578vbKyknJycqQBWFVV5bIOX1tf4ZnBJ554or7Thj9+Lv8SvZl6kt7ZdpqKSitd6vPKhmMUFRpI43u0pkm9Eumafm0oNCjAZV0cBAEQAIGGCOQUl9Gnu8/S+zvS6fD54nqr//dn+2lCj0S6fmBbmtw7icKCcd+pFxZOgIAHBExpALLevFRcs6iVbj5e83PtOrWvq3l+wYIF9NBDDzkO8QxgSkqK47tRP+w/W0D/t/E4fbUvg6rE0zaXlPgw6to6kuLFEzc/dceFB9HxnIu0/lCWnBH86scM4tc/vj1Kz9zcn4Z1ijeq+pAbBEBAAwJZhaX0xJc/0df7M+UsH4sQEtiSruyZSLzaENiyBQWIFxeeFTyYWUTr0s7LV2RIID10dQ/6zZhOde718gL8AwIg0CABUxqAycnJVHsmLysriwIDA6lVq1bSAAwICHBZp/bMYU2CISEhxC+zFDaE/yUMv2dWHyS73UejurSieWM701Vidq+l/eZbU182EPeeyafvDmbJJ/YTwii89V+bae6oTvRfU3tShLgxo4AACICAOwKbjmTTgyv2yG0mXG9ASizdOrQ9zejflmLCglxeekgYgJ/tOSte5+isWK1YKIzHjaKdZ28eQK2jzHNfdqk8DoKADwiYciPXqFGjaO3atU641qxZQ0OHDqWgoCAKDg6WIWJq1+Hvo0ePdrrOrF8ulVfR79/bQ39bZTP+rumXTF/eN5bevWskTe6T5NL4Yxb8RD64Qxw9PKUnrXlwAs0cmiIMaqI3xNLx1Bc3UurRHLMig14gAALNJFBZVU3PfX2I5izbJo2/XslR8r7z2e/G0K9HdKzX+ONue4q6f5jWizb94Up68vor5Gzh+kPZNP2ljXJlopmi4XIQsBwBQxiAxcXFtGfPHvniEeIwL/xdhWzhpdk5c+Y4Bo89fk+dOiWXa9kTeNmyZdIB5JFHHnHU4aXcpUuXynNc58EHH5Tt8bVmL2fySuiXS1Lp873n5DLLkzf0pX/OGkx928U0SnV+Uv9fsfz71rzh1C42jM7kXaLbX9sqn9Ib1RAqgwAImJ7AebHkO2vpVnr5u6PyofFXwzvQp8Lwa+x9h1cmZosVh8/vHUs9k6KkIXnH69vpqa9+kk4jpgcJBUHASwQMEQZm/fr1dOWVV9ZRee7cufTGG2/QHXfcQSdPniSupwoHgmaj7sCBAzKANIeGqW3ccSDoZ555hjIyMqhv377097//ncaPH6+aaPDdiG7kO05eoLve2kkXhHdvK7HPZvGvB9MIsezb3FIsPIf/+9P99LHYzM0rx//vV4Pkck5z28X1IAACxieQXVQmHzpPXyihCOG8seiX/ekXA9o2W7HSiiq5isErEFx4T+B/z+iDfYGSBv5xR8CIf7/d6dOUc4YwAJuimD+uMdoP6Mj5IrpJzPyxh2/fdtH0r9lD5cydt1hxyIY/fvQjfbDzjFwqflkYgdOFlzAKCICAdQkUlVbQbf+3hQ6cE05zwrns3/8xgjonRHgVyEfinvPwB3tlm+wccv+k7l5tH42Zj4DR/n77YgQMsQTsC8Wt1iY/gfMyCRt/wzrF0Qe/He1V44958tLM38ST/U2D20lv4vve3U1rDtQfVsdqYwB9QcBqBMoqq+i3YsWBjT9ecXjLB8YfM/3lkPb01+v6SLwvrD0sw1lZjTX0BYHGEoAB2FhiBqzPDh/z/71Des51EgFVeebPVzG02EmEvfI4VhcHcP3dO7voOxE6BgUEQMBaBDhiwEMr9lLqsVy57PvGb4ZTJy/P/NUk+psxnemBybaZv79+fkDGFqx5Hp9BAAScCcAAdOZhum+8LMvhFvam54t4fkH0urgJx4sncV8WNgKfv2WA2APYhiqqLtPvxUwgh21AAQEQsAYBDjH1xBcHZGzRoIAW8qGzX/vGOZk1hdTvxdLvHaM7yUt5SRgPn02hiGusQgAGoMlH+n9FjL/VYhk2WKRu+z8x8+ftvTf14QsU/f195kAZ36tQLDuzEcghIFBAAATMT2DF9nT69+ZTwhmD6IVbB9LY7gl+UZoD+bMTyE2DbNtQHhIPvxxwGgUEQKAuARiAdZmY5sgqkdmDAz1z4Wwdwzv7N1sH5wv+x22DKEoEh95xKo9e+uaIadhCERAAAdcETueWyCDNfJaDw1/nBW9f1z25Psp7kRf9sh/1bhMtsxb914c/OrI/ub4CR0HAmgRgAJp03DnMy2MiLAuXuyd0pRvEE7EWpYPYc/jUTf1k1xz/K/UYAkVrMQ7oEwT8QUDu+3t/D5WIfcf8wPnb8V390W2dPkICA+il2wZSsEgtt+FwNr215VSdOjgAAlYnAAPQpL8A3gSdK4zAHkmR9ODV2oZE4HhfKmPIAyL7SK5I/o4CAiBgPgKvbjouZ/s5Vy/vA1a5fLXQtIcIEr1gei/Z9VNfpdHRrCItxECfIKBbAjAAdTs0TRds9f4M+kJk+eCb73PiJsxPw1qXv/6iD3VtHUFZIhzNI2JzNjunoIAACJiHQFpGIb2w5rBU6L9FSJaU+HDNleMc5ePE/sOyymp6QOwHLBfvKCAAAjYCMABN9kvIq7H0+9vxXah/+1hdaBgeHEgvi3RzvCTzncjf+YnIGIICAiBgDgIc74+jDZQLR6/JvZPoFhGXTw+F9wPyQzBHQNh/tpBeXGczUPUgG2QAAa0JwADUegS83P/jIvRCTnE5dU+MpN/bY2J5uYsmN8ebslWcrr8J72TOEIACAiBgfAJ/X3uEDmYWyWDPfxMOGOyNq5eSFB1Ki2607UNesuGYMAQL9CIa5AABTQnAANQUv3c7/1qEe/lszzmZi/dZnSz91tZw3tjOMhQNZyb5x7dHa5/GdxAAAYMROJZdTEvF3j8uTwlDKyEyRHcacEpK3osswhNKD2WOU4gCAlYnAAPQJL8ATor+uHD84HKX8LwbmKKPpd/aeHk/Isfp4rLs+xNiY3Zx7Sr4DgIgYCACi1YelFl/JvVKpGl9k3Ur+Z+EQ0hoUEvaduICrd6PFJW6HSgI5jcCMAD9htq3Hb3+w0nKKCiV+X3VMqtve2x661eKPxT8x4JTxS388ifE6Go6SlwJApoSSD2aQ+vSzkuHswXX9NZUloY6bxsbJh+Oud7Tq9KIH5pRQMDKBGAAmmD02fFj8XrbcurDU3qIp1ztvX4bwvoXMQvI2Uk2ihhd69KQK7ghXjgPAnojwDH//keEV+Fy+4gO1E3sO9Z7uXtCF0qKDqH0C5eIH5pRQMDKBGAAmmD0OcBykUi3xk4WNwzUJuBzYzFyUvj54zrLy54Us4B4Gm8sQdQHAW0JfLzrDP0kQr9EhQYKh7Me2grjYe8cjeCP02yxAf8p7ptZRUgT5yE6VDMhARiABh/U9AslIufmSakFBz3lsAdGKb+7shslCw+900KHV+0p64wiO+QEASsTKCmvpGe/PiQR3HdVN4qPCDYMDn5IHtA+horLKh1xCw0jPAQFAS8SgAHoRZhaNPXcmkNUUXWZxnZLoPE9WmshQpP7jBDZAhZcY3sa/z9hABaUICxMk2HiQhDwI4F/bTgug7qnxIfR3NGd/Nhz87vih2QOVM1lxY50OnAOYWGaTxUtGJEADEAjjppd5n1nCmTYF/7KHm5GLNf1b0u9kqOoSDyNv/bDCSOqAJlBwFIEMoWz2b82HpM6/2lab11kGmrsAAzpGE/X2cPC/G3VwcZejvogYAoCMAANOowcx2qR8GTjcsPAttS3XYwhNeGn8fsn2XIVvy4MwIJLmAU05EBCaMsQYIez0opqGtIxjq7pp9+wLw0NyB+m9pTey5uO5NDu03kNVcd5EDAdARiABh3S70X4hdRjudKT9uEpPQ2qhU3saVckU0+RuJ0dWdgIRAEBENAnAXaaeG97uhSOIw7oKeNHY4lxruIbB9mc5tghBAUErEYABqBBR1zdsGaJ8At6SLreHIw8C3jfpG6yCQ4OXYgUcc3BiWtBwGcEXhP/P8srq2lQh1ga1aWVz/rxV8P/ObGrMGJJhqL66Vyhv7pFPyCgCwIwAHUxDI0Tgpcrthy/QIHCcLprfJfGXazT2tP7tpFxxArFLOAbP5zUqZQQCwSsS4CdtN7efEoC+N3Eboae/VOj2LV1JF0r0sRx+ac9lqo6h3cQMDsBGIAGHOFXREJzLteLcAYc3d4MJYBnAUU4CS48y1CEWUAzDCt0MBGBNzefpIvlVdJpa1LvRNNoxuGouKzcl4HUlKYZVSjiCQEYgJ5Q0lEdzp275qfzUiKOam+mMkN4BHdtHSEdQd5MPWkm1aALCBiawEXhpb/Mvj/3HmEwGXnvX+2B4AD6k3sniZSUREvW2x6ua9fBdxAwIwEYgAYb1f8T4Rf4RsU3rO7CccJMxTYLaPMIXipmATlQKwoIgID2BN7ddpryxRJwp1bhjiVT7aXyngT32lcfPt1zVqSJK/Few2gJBHRMAAagjgentmgcf+uT3Wfl4f+caK7ZP6Urx+bqItLE8R+b98QfHRQQAAFtCZRVVhEHaufCThP8oGa2MjAllsZ1TyDOb7zEvsXGbDpCHxCoTQAGYG0iOv7OSzCc9WN4p3gRgytex5I2XTT+4zJ/nM245T1HfENGAQEQ0I7ARzvPyqwfbWJCRdiU9toJ4uOe77vKtvrw4Y4zxA/bKCBgdgIwAA0ywuyBt3yLzQPvbpPO/qmh4NhcMWFBYinmEn17MEsdxjsIgICfCfADmHI6u1M8mAUHmvdPxvDO8TSsUxyVV1XT8q22e62fcaM7EPArAfP+b/YrRt939ra4IbEHHgdMvrKneTzwXJELCw6g24anyFNvpCIwtCtGOAYC/iDwnXgAOy32xPEDmfo/6Y9+terjP8Z0ll3znkde+kYBATMTgAFogNHlwKuv/3BSSsqzf2bywKsP/+yRHYm3Gv1wNJcOZRbVVw3HQQAEfEjg3/ZVh5nDUig8ONCHPemj6av7JBEvdecUl8uwMPqQClKAgG8IwAD0DVevtrr6QKa4IZVRYlQIcagUK5T2ceE0VaSI4/JG6kn5jn9AAAT8R+B4djFtPJwtM2XcPqKj/zrWsKfAgJb0a5FdicubqVgG1nAo0LUfCBjGAFy8eDF17tyZQkNDaciQIbRp06Z68UycOFHOkvFMWc3Xtdde67jmjjvucDrH9UaOHOk4r6cPb9ufwm8b3oGCxA3KKuWO0Z2kqp/sPiO8gsutojb0BAFdEHh7y2kpB2856SDCv1il8H02WNxn96Tn017xQgEBsxIwhDWxYsUKeuCBB+jRRx+l3bt307hx42j69Ol0+rTtBlV7cD7++GPKyMhwvPbv308BAQF0yy23OFWdNm2aow7XX7lypdN5PXw5fL6Itp24IEMv/Mq+L04PcvlDBt6U3UcEaS2tqHYkoPdHv+gDBKxOoKS8kj7YmS4xzB5ljdk/NeYJkSF0bX9bejiORIACAmYlYAgD8IUXXqB58+bR/PnzqXfv3vTiiy9SSkoKLVmyxOW4xMfHU3JysuO1du1aCg8Pr2MAhoSEOOpwfb5Ob0V5/k4WqZfaxJgj7ZunjHlW9jdjOsnqb4kcpJXCOw8FBEDA9wQ+23NOpGOspI5i5m9C99a+71BnPcy1rz58uTeDcsX2GxQQMCMB3RuA5eXltHPnTpoyZYoTf/6emprqdKy+L6+99hrddtttFBER4VRl/fr1lJiYSD169KA777yTsrLchxwpKyujwsJCp5dTg17+wumXPt5lC/x8u3CKsGLhwNCtIoLpbP4lWmtPgWdFDtAZBPxF4LJINfRv8cDFhff+tTRh4OeGWHJg6AHtY2RImPe222ZCG7oG50HAaAR0bwDm5ORQVVUVJSUlObHl75mZmU7HXH3Ztm0b8RIwzx7WLLyEvHz5cvr222/p+eefp+3bt9NVV11FbOTVVxYtWkQxMTGOF89C+rLIp3BhBHL6pTFdE3zZlW7bDg0KoFn2TdmvwxlEt+MEwcxDYMepPErLKKQQEfPvlqHmDfzc0IjNGdVJVuFVGKw+NEQL541IQPcGoILKy4E1Cz+l1j5W87z6zLN/ffv2peHDh6tD8n3mzJnETiF87rrrrqNVq1bR4cOH6auvvnKqV/PLggULqKCgwPFKT/fdkyHrp5w/fm3Rp3DFnvXnDCG8F/KY8ExEAQEQ8B0BNft3w8B2FBse7LuOdN4y7wOMF6sP50RWkHVp7leHdK4KxAMBlwR0bwAmJCRIB47as328XFt7VrC2hiUlJfTee+/Vmf2rXY+/t2nThjp27EhHjhxxdVoe4z2D0dHRTq96KzfzxG7hffaTeArnyPs3D7HuUzhjTBZxuSb2sO1Den+H74zuZg4ZLgcBwxPIKiql1fszpB5Wc/6oPXi8+qAc7/69+WTt0/gOAoYnoHsDMDg4WIZ9YUeOmoW/jx49uuahOp/ff/99uaR7++231zlX+0Bubi7xjB4bgnooavbvOhH3L048hVq93CoC0XLhvKQVcAax+s8B+vuIwIpt6TLf+OAOsdS3XYyPejFOs7PE6gMvPqUey6XTuSXGERySgoAHBHRvALIODz30EC1dupSWLVtGaWlp9OCDD8oQMHfffbdUcc6cOcTLs7ULL//ecMMN1KpVK6dTxcXF9Mgjj9DmzZvp5MmTxM4gvAzMs4033nijU10tvuRdLKcvf7Q9hd8+0haUVAs59NTnVb0SicMzcEBsTk+FAgIg4F0CvO3kg51nZKO87QKFqF1sGI3tZtt//eEuGxtwAQGzEDCEAcj79Tj0y8KFC2ngwIG0ceNGGbOPl2y5cDxAjuNXs/B+vu+//16Gj6l5nD9zTMB9+/bR9ddfLz2A586dK9/ZIIyKiqpd3e/fPxQ3YU7/dkXbaGJvNBSSAbB/ObidRIFlYPwiQMD7BHiPLef9jRC5uKf3s2Xh8X4vxmtRbcH5SNyXq6svG08BSAwC9RAIrOe47g7fc889xC9XhWfwahcO7cJPtK5KWFgYff31165O6eIY7/tLiAwmDv3iiaOLLoT2gxC3DE2hf208Tt8dyqaswlJKjA71Q6/oAgSsQYAfPLmw84MV8v56OqqckjIqNFCGotp8PJfG2GcEPb0e9UBArwQMMQOoV3i+kouDkKb+aRLdZJ/x8lU/Rmu3W2IkDekYR1XiKfwje3xEo+kAeUFAjwQ45uhX+2yrKPyghfIzAXYG+YWIR8rlAzih/QwGnwxPAAagToeQZwFDAgN0Kp12Ys20/3HiG3F9M7zaSYeeQcCYBFbtz6SS8ioZc3SoeMhCcSZwq/2+w5wKLlU4n8Q3EDAoARiABh04q4p9jVyeCqDjORdp+8k8q2KA3iDgVQIf2vP+8n43bDupi7a/yArSIymSysTe7C9/PFe3Ao6AgAEJwAA04KBZWeTIkECaYU/UvgIpmqz8U4DuXiKQLhw/thy/IMOd3DjY2jFH60PKRvEtQ2xL4x/sgDdwfZxw3FgEYAAaa7wgrSAw0x4TcKXYs1RUiuUY/ChAoDkElPMHp5vksCcorgncMKidzEi0RwTpP5pV5LoSjoKAgQjAADTQYEFUG4HBHeKoS+sIulRR5YiXCDYgAAKNJ8BhTT6yx7ezct5fT8i1jgqhK3smyqqYBfSEGOronQAMQL2PEOSrQ6Dmcsxne87WOY8DIAACnhHYciKXzuRdoiixtWJKH8T+a4iaMpI/3n2WKpGRqCFcOK9zAjAAdT5AEM81gesG2FL2bRXBazMKLrmuhKMgAAJuCajl3xni/1OYCACN4p4AZyRqJVJzZheV0YbD2e4r4ywI6JwADECdDxDEc02gfVw4DesUJ0LBEH251zkLjOsrcBQEQKAmAY79t2pfpjx0s93BoeZ5fK5LICigJV0/0JaR6LM98AauSwhHjEQABqCRRguyOhFQN+JPsQzsxAVfQMATAmt/Oi/30XZqFU6DOyDlpCfMuM4vBtqCQjO/kvJKTy9DPRDQHQEYgLobEgjkKYFr+rWhwJYt6MC5QnjleQoN9UDATuCLvbYZLM5ygdh/nv8sBoiYgB2F0cxOaOvSsjy/EDVBQGcEYADqbEAgjucE4sVenPE9WssLPsdyjOfgUNPyBPJLymnjEdsetuvsac4sD8VDAGwsX9ffNguI+46H0FBNlwRgAOpyWCCUpwSuty/HfCZmM5AazlNqqGd1AqtFSrOKqsvUKzmKuidFWR1Ho/VXy8AbDmdRQQlikTYaIC7QBQEYgLoYBgjRVAKTeydRmEjWfiq3hPaeKWhqM7gOBCxF4At7OjPM/jVt2HsIo5mNZzaiVx+AE1rTKOIqrQnAANR6BNB/swhEiPhlV/dJkm0gJmCzUOJiixDIKiqlzcdypba8/w+laQSU8fwFohA0DSCu0pwADEDNhwACNJfADYNsf8T4RlwlMhuggAAI1E9g5Y8ZxP9NBqbEUkp8eP0VccYtAbUPMPVYDrFRjQICRiMAA9BoIwZ56xAY1701xYUHUU5xmWNmo04lHAABEJAEvhAGIBc1gyW/4J9GE+ggPIHZiGZjmo1qFBAwGgEYgEYbMchbhwAHZ+WQMFywDFwHDw6AgIPAmbwS2nkqT4R9IZrR3/Z/xnESHxpNQC2hf24PqdPoBnABCGhIAAaghvDRtfcIqKDQ7N1YKuJzoYAACNQl8KV9pmpE53hKig6tWwFHGkWAjWgRipR2nc6n9AsljboWlUFAawIwALUeAfTvFQJDO8ZRsviDViTSW/1wNMcrbaIREDAbgZ+DP9vSmZlNP3/rkyjuOSO7tJLdKs9qf8uA/kCgqQRgADaVHK7TFYGW4jF8Wt9kKdNKe35TXQkIYUBAYwLHsotl1hzOnqP+r2gskim6dywDIxi9KcbTSkrAALTSaJtc1+l2A3DtT5lUXlltcm2hHgg0joCa/RvbPYE4iw6KdwhM72tLSXkws4hO5Fz0TqNoBQT8QAAGoB8gowv/EBjaKZ4SIkOosLSSNh+3xTnzT8/oBQT0T2CVfWZ8hj2Nmf4lNoaEMSICwaiutmXgVfvhDWyMUYOUTAAGIH4HpiEQIJa2pl5hCwq9ah9uxKYZWCjSbALHxfLvofNFxMu/V4vsOSjeJcCzgFy+Fk5oKCBgFAIwAI0yUpDTIwIqHMyan85TZRWWgT2ChkqmJ7D6gM0w4ZkqnrFC8S6BKeLBk0PrcDrKs/mXvNs4WgMBHxGAAegjsGhWGwIc3oKDQl+4WE7bTlzQRgj0CgI6I6BmptRMlc7EM7w4vPVkmNiCwoVDUaGAgBEIwAA0wihBRo8JBIqg0FP62L2BsR/HY26oaF4CPCPFM1M8Q6XyZptXW+00U05oq3Hf0W4Q0HOjCMAAbBQuVDYCgen9bAbg1wfOIzewEQYMMvqUgJr9G9YxnlpHhfi0Lys3PvUK231nh8i0gtzAVv4lGEd3GIDGGStI6iGB0V0TKCo0kLKLymTaKw8vQzUQMCUBtf8Psf98O7xtY8NogMgNfFnkBl4jHj5RQEDvBGAA6n2EIF+jCQQHtnQsda2EN3Cj+eEC8xDgh6DtJ217Yafa42SaRzv9afLzMjD2AepvdCBRbQIwAGsTwXdTEFCb3b8W3o/V1eKRHAUELEhgrfCG5xmp/u1jqJ2YoULxLYFp9mVgjkOaJxzRUEBAzwQMYwAuXryYOnfuTKGhoTRkyBDatGlTvVzfeOMNseG5RZ1XaWmp0zWNadPpQnzRPYFxIttBRHAAZRSU0p4z+bqXFwKCgC8IYPnXF1Trb7NTQgT1So6Se4/XpWEZuH5SOKMHAoYwAFesWEEPPPAAPfroo7R7924aN24cTZ8+nU6fPl0vw+joaMrIyHB6sfGoSlPaVNfiXf8EQoMCaJI94C3CMuh/vCCh9wkUXKqg1KM5smE1M+X9XtBibQJq9QH3ndpk8F1vBAxhAL7wwgs0b948mj9/PvXu3ZtefPFFSklJoSVLltTLk2cAk5OTnV41KzelzZrX47P+CSivPNsyGJaB9T9ikNCbBL4RM1CVYvtDj6RI6tI60ptNoy03BFQUgk1HcqiotMJNTZwCAW0J6N4ALC8vp507d9KUKVOcSPH31NRUp2M1vxQXF1PHjh2pffv2NGPGDDlzqM43tc2ysjIqLCx0eqk28a4/AhN6tqZgEReQE7QfE6mwUEDASgTUDBRm//w76t0ThcEtloLLRSaibw9m+bdz9AYCjSCgewMwJyeHqqqqKCnJOX8lf8/MdO1p1atXL+J9gJ9//jm9++67ct/gmDFj6MiRIxJNU9rkCxctWkQxMTGOF89CouiXQGRIII20J2nn1HAoIGAVAiXllbThcLZUd5o9T61VdNdaT159Uh7XvPqAAgJ6JaB7A1CB4/9UNctl4dpW+5g6P3LkSLr99ttpwIABcr/g+++/Tz169KB//OMfqop8r329uzb5ggULFlBBQYHjlZ6e7tQevuiPgMp8sA43Yv0NDiTyGYGNh3OorLKaOsSHU+82UT7rBw27JqDuOxsOZVO5GAcUENAjAd0bgAkJCRQQEFBnti8rK6vOrGB9gFu2bEnDhg1zzAA2tc2QkBBi55Kar/r6xHF9ELja7giyOz1fBobWh1SQAgR8S0B5oE4Wv//aD7q+7RmtM4GB7WOJ8wMXlVUiJzl+EroloHsDMDg4WIZ9Wbt2rRNE/j569GinY/V94Zm9PXv2UJs2bWQVb7RZX184ri8CyTGhMgYax0LjTfEoIGB2AlXC8UPtPZvcJ9Hs6upSv5YtW9Dk3jb2yhjXpaAQytIEdG8A8ug89NBDtHTpUlq2bBmlpaXRgw8+KEPA3H333XLw5syZI5dn1Ug+8cQT9PXXX9Px48el4ccexGwAqvqetKnawrvxCfAsCBfciI0/ltCgYQK7T+fRBRGEOFqkQxzWKb7hC1DDJwTUfQdRCHyCF416gUCgF9rweRMzZ86k3NxcWrhwoYzr17dvX1q5cqX08uXOOR4gL/Oqkp+fT3fddZdcNmanjUGDBtHGjRtp+PDhqgo11KajIj4YngDvx3lh7WHisAy8OT482BA/e8NzhwLaEFiXZvM8ndgzkYKEFzyKNgTGdEug0KCWdDb/Eh3MLBJ7MaO1EQS9gkA9BFqI5VEESKsHTkOHOSQMG5jsGML7AlH0SYB/4uOe+Y7O5F2if80eQio+oD6lhVQg0DwCk1/YQEeziun//WoQ/WJA2+Y1hqubRWD+mzvkysPDV/eg+yZ1b1ZbuNi7BPD3mwiPh979TaE1HRLgTfDKKw9hGXQ4QBDJawROipiXbPwFij1oE3q09lq7aKhpBK6278HE9pOm8cNVviUAA9C3fNG6Tggob2DeHM+b5FFAwIwElKExvHM8xYQFmVFFQ+l0VS/2wibae6aAzhc656I3lCIQ1pQEYACaclihVG0Cw8QfRN4Uz5vjd4lN8iggYEYCygBUDghm1NFIOrWOCqGBKbFS5G/sezONJD9kNTcBGIDmHl9oZyfAm+Gv6mULy4BlYPwszEggv6Sctp+0PdzAANTPCKuxUMa5fiSDJFYnAAPQ6r8AC+l/dZ9kqS3CMlho0C2k6nqRdYK3N/RIiqQOrcItpLm+VVX7j78/aotCoG9pIZ2VCMAAtNJoW1zX8T0SKFjMBJ4QG+WPZV+0OA2obzYCaoZJzTiZTT+j6tM9URjkIiUfp4TjUFQoIKAXAjAA9TISkMPnBKJCg2hEF1tg3O+EMwgKCJiFABsXnHeWy2QR9xJFPwQ4CoEyypGTXD/jAkkQBga/AYsRuFIEx+WiUmVZTH2oa1IC205ckHlnEyKDZR5ak6ppWLVUSr5vEIXAsGNoRsExA2jGUYVO9RJQjiDbT16gwtKKeuvhBAgYiYBa/uXfN+ehRdEXAU7Jp6IQ7EnP15dwkMayBGAAWnboral4p4QI6iJelWKz/PfYj2PNH4HJtOZMN98cPC+1mmTPe20yFQ2vDkchGGcPzL3+ELafGH5ATaIADECTDCTU8JzAlfZwMFgG9pwZauqXADs0pV+4JB2cxor8syj6JHAVtp/oc2AsLBUMQAsPvlVVV8vA/CRejawgVv0ZmEZvNaPEDk4RIYGm0ctsikzo2VpmBTlwrhBZQcw2uAbVBwagQQcOYjedAO/HiRR/KHOKy2nf2YKmN4QrQUAHBNRM9kT7DJMORIIILggkRIZQ//a2rCDKaHdRDYdAwG8EYAD6DTU60guB4MCWpJbK1B9PvcgGOUCgMQSKhCMTOzRxUTPbjbkedf1LAMvA/uWN3twTgAHong/OmpSA+mMJA9CkA2wRtX4Q2SUqqi5TJ5H5o7NwbkLRN4Ere7WWArIDGsduRAEBLQnAANSSPvrWjMBE+42Yl4CzCks1kwMdg0BzCHx30Bb8WTk2NactXOt7An3bxhAvBV8sr3LM3Pq+V/QAAq4JwAB0zQVHTU4gMSpU7MeJkVpyDlUUEDAaAQ7/8p09pIgKcG40HawmL8donCicQbggG5HVRl9/+sIA1N+YQCI/EVB/NLEM7Cfg6MarBNibNKuojMKCAmh4Z1uKQ692gMZ8QsCx/QTxAH3CF416TgAGoOesUNNkBNSN+Huxjwr7cUw2uBZQR3mSjhGx/0KFEYhiDAJjuydQoJgJPC7iN57KvWgMoSGlKQnAADTlsEIpTwj0a2fbj1NcVon9OJ4AQx1dEVAz18qxQFfCQZh6CUSHBtHQTnHyPJaB68WEE34gAAPQD5DRhT4J1NyPo/6Y6lNSSAUCzgQuXCyn3facsmorg3MNfNMzAbX68C32H+t5mEwvGwxA0w8xFHRHQN2I1XKau7o4BwJ6IbDpSDYJHxDqlRxFbWPD9CIW5PCQgDLatxzPpZLySg+vQjUQ8C4BGIDe5YnWDEaA908FiP04tnyqJQaTHuJalYCasUb4F2P+ArolRlL7uDC59zj1aK4xlYDUhicAA9DwQwgFmkMgJiyIBnewpWfacBjhYJrDEtf6h0CVyF+tfqtqJsk/PaMXbxFo0aIFqbH7Ft7A3sKKdhpJAAZgI4GhuvkITOhhi8uFeIDmG1szarRH7P3LL6mg6NBAx8OLGfU0u07KeWejePDkmI4oIOBvAjAA/U0c/emOwMSeiVKm1GMIB6O7wYFAdQhssM8YjRMPLoEBuIXXAWSQAyO7tKJgMX5n8i7JLSgGERtimogA7h4mGkyo0jQCfdpEi/RMwWIzdhXtOHWhaY3gKhDwEwG1/Ktmrv3ULbrxMoHw4EBHAG81pl7uAs2BgFsCMADd4sFJKxDgcDDj7cvAGxCWwQpDblgdc4vL6EeRv5oLDEDDDqNDcDWGMAAdSPDBjwRgAPoRNrrSLwG1DIwbsX7HCJIRcdYaFf4lKToUSAxOQOUF3irCwZRWVBlcG4hvNAIwAI02YpDXJwTGiXAwYiKQDmYWUUbBJZ/0gUZBoLkE1Az1hJ42x6XmtofrtSXA4WDaxoRSWWU1cUxAFBDwJwEYgP6kjb50SyAuIpgGpNjCwbBXHgoI6I1AtQj/slEEgOailg71JiPkaRwBDgejjHmsPjSOHWo3n4BhDMDFixdT586dKTQ0lIYMGUKbNm2qV/tXX32Vxo0bR3FxcfI1efJk2rZtm1P9O+64g/g/X83XyJEjnergi7UIqD+qCAdjrXE3irY/ZRRSTnE5hQcH0NCO8UYRG3I2QEDdd2AANgAKp71OwBAG4IoVK+iBBx6gRx99lHbv3i2Nu+nTp9Pp06ddAlm/fj396le/ou+++442b95MHTp0oClTptDZs2ed6k+bNo0yMjIcr5UrVzqdxxdrEVD7AL8/kkMVVdXWUh7a6p6AMhBGd02g4EBD3Lp1z1QPAo62ZyM6nn2R0i8gG5EexsQqMhjiLvLCCy/QvHnzaP78+dS7d2968cUXKSUlhZYsWeJynJYvX0733HMPDRw4kHr16kU8I1hdXU3ffPONU/2QkBBKTk52vOLj8VTtBMhiX/q1i6G48CAqKqskDraLAgJ6IqAMQLVkqCfZIEvTCUSHBtGQDnGyATXGTW8NV4KA5wR0bwCWl5fTzp075QxeTbV4Ri81NbXmoXo/l5SUUEVFBdU28HimMDExkXr06EF33nknZWVl1dsGnygrK6PCwkKnl9sLcNJQBDgnsAoHsx7pmQw1dmYXtrC0gnadypNqTugOBxCzjbcy6mEAmm1k9a2P7g3AnJwcqqqqoqSkJCeS/D0zM9PpWH1f/vSnP1G7du2I9wKqwkvIPFP47bff0vPPP0/bt2+nq666Shp5qk7t90WLFlFMTIzjxbOQKOYigP045hpPs2iTejSXKoUTSJeECOrQKtwsakEPOwF130kVYX7KhUcwCgj4g4DuDUAFgZ01ahbOnVj7WM3z6vMzzzxD7777Ln388cfSgUQdnzlzJl177bXUt29fuu6662jVqlV0+PBh+uqrr1SVOu8LFiyggoICxys9Pb1OHRwwNgE1A7j/bCFlFZUaWxlIbxoCamZI/T5NoxgUkQRs2YhC6KLIRrTTPtMLNCDgawK6NwATEhIoICCgzmwfL9fWnhWsDeu5556jp59+mtasWUP9+/evfdrpe5s2bahjx4505MgRp+M1v/CewejoaKdXzfP4bHwCCZEhxHsBubAzCAoIaE2AH3ZVaCK1VKi1TOjfuwRs2YgSZKPrD7vfiuTdntGalQno3gAMDg6WYV/Wrl3rNE78ffTo0U7Han559tln6cknn6TVq1fT0KFDa55y+Tk3N5d4Ro8NQRRrExjfw3YjVn90rU0D2mtN4Fh2MZ3NvyQ9f0d2bqW1OOjfRwTUMrAK9u2jbtAsCDgI6N4AZEkfeughWrp0KS1btozS0tLowQcflCFg7r77bqnInDlziJdnVeFl38cee0zW79Spk5w95P2CxcXFsgq/P/LIIzJEzMmTJ4mdQXgZmGcbb7zxRtUM3i1KYJx9k/0mMQPIwXdRQEBLAiou5YjO8RQmYgCimJMA33d4pxNnIzpfiO0n5hxlfWllCAOQ9+tx6JeFCxfK0C4bN24kjtnHS7ZcOB4gx/NThYNGs/fwzTffLGf0eFaPX7wkzIWXlPft20fXX3+99ACeO3eufOeYgVFRUaoZvFuUwGARkiFC/KHNvVhOHHwXBQS0JKD2/6kZIi1lQd++IxAvshH1b49sRL4jjJZrEwisfUCv3zmuH79cFZ7Bq1l4Vs9dCQsLo6+//tpdFZyzMAEOsjtKBNtdl3Zept7qa98TaGEkUF0jAqUVVbT1xAXZOwxAjQbBj92O755Ae0UMUl59uGUookz4Eb0luzLEDKAlRwZKa0oA+wA1xY/O7QS2CeOPw4K0iQmlbomR4GJyAsrL+3sRDgbbT0w+2DpQDwagDgYBIuiPwHj7PkAOyXBRZAZBAQEtCChHpHFiZsiTsFdayIg+vUdgYEosRYYE0gWx/eTAOWw/8R5ZtOSKAAxAV1RwzPIEOnHA3fhwkRP4Mm05nmt5HgCgDQFeCuSiZoa0kQK9+otAUEBLGt3V5um98Ui2v7pFPxYlAAPQogMPtRsmwLMuXNQsTMNXoAYIeI8Ae4IeOl8kPUPHiD2pKNYgMK6HLdUf7jvWGG8ttYQBqCV99K1rAmrWZSMCQut6nMwqnDIA+gsnpDjhIYpiDQLsCMKFt58UY/uJNQZdIy1hAGoEHt3qnwAvxQS0bEEnci5S+oUS/QsMCU1FAMu/phpOj5Xp2CqCOop8z5z7ecsxbD/xGBwqNpoADMBGI8MFViEQFRpEgzvY43JhP45Vhl0XerIHKHuCclGByXUhGITwCwHH9hPcd/zC26qdwAC06shDb48IKG9gtRzn0UWoBALNJMAeoOwJyh6hg+wPIc1sEpcbiIC676hZYAOJDlENRAAGoIEGC6L6n4DaB5h6NFd4BFf7XwD0aEkCygN0lNiGwJ6hKNYiwOMeiO0n1hp0DbTFnUUD6OjSOAQ4C0hseBAVic3Ye0SEfhQQ8AcBNeOsHAL80Sf60A8B2/aTOCmQehjQj3SQxCwEYACaZSShh08IsBPI2G42r7xNhxGXyyeQ0agTAfb83HU6Tx5TM9BOFfDFEgQc+wBx37HEeGuhJAxALaijT0MRUH+ENyAcjKHGzajCbhWBxzkAOQciZ49QFGsSUPcd3n5Sie0n1vwR+FhrGIA+BozmjU9APYn/eCaf8kvKja8QNNA1Acfybw8Ef9b1QPlYOGw/8TFgNE8wAPEjAIEGCLSJCaPuiZF0+TJRKuJyNUALp5tLQHl+IvxLc0ka+/qa20/UQ4GxNYL0eiMAA1BvIwJ5dElgrD06/ybE5dLl+JhFKA44flwEHuc//ionrFl0gx6NJ6DCwSAbUePZ4YqGCcAAbJgRaoAAOW7Eh3PETKCYCkQBAR8QULN/HICcPUFRrE1APXjy9pOCkgprw4D2XicAA9DrSNGgGQmM6BIv4rG1oLP5l+hkLtLCmXGM9aDT90dtnuZju7XWgziQQWMCbWPDqGvrCBKJYcT2E1tmGI1FQvcmIgAD0ESDCVV8RyA8OJCGdLTF5cIysO84W7nlKk7/Zvc0Hw8HECv/FJx0V3tBsQzshAVfvEAABqAXIKIJaxBw3IjFMjAKCHibAC/zFZZWUnRoIPVvb8tB7e0+0J7xCKiHAX7wxPYT442fniWGAajn0YFsuiKg9gFukXHakBZOV4NjAmHU7N8YEXicnUBQQIAJjOjM6QBb0Jm8S3QK20/wo/AiARiAXoSJpsxN4Iq20RQn0sJxpgakhTP3WGuhnXIAURv/tZABfeqPQERIIA3uYN9+chSrD/obIeNKBAPQuGMHyf1MoKWYleHZGS5IC+dn+Cbvzin9W3c4gJh8uButnsoKgvtOo9HhAjcEYAC6gYNTIFCbgFoGxobs2mTwvTkEtogA45XCCaRjq3BKESngUECgJgGVj3wz/06QFq4mGnxuBgEYgM2Ah0utR0AtzyEul/XG3pcaf29f2lNpB33ZF9o2HoGaaeH2CmchFBDwBgEYgN6giDYsQ4DjcnUTaeEQl8syQ+4XRTfaM8woT3O/dIpODEOAnYLGdLVtP9mIKASGGTe9CwoDUO8jBPl0R0Atx2AZWHdDY0iBOLj48Wxb+rdRXVsZUgcI7XsCanZYzRb7vkf0YHYCMADNPsLQz+sEEJfL60gt3eD39tm/gSmxIgYg0r9Z+sfgRnm1/YQjEBRcQlo4N6hwykMCMAA9BIVqIKAIIC6XIoF3bxBQM8lqZtkbbaIN8xFoHxdOXRIiiDPGsDMICgg0lwAMwOYSxPWWI+AUl8s+e2M5CFDYKwT4j/kPdgcQNbPslYbRiCkJ/LwMbMsZbUoloZTfCMAA9BtqdGQmAo64XPbcrWbSDbr4j8CBcwWUX1JBUSLY7wCkf/MfeIP2NNYeI1IFDTeoGhBbJwRgAOpkICCGsQio5TrE5TLWuOlNWvWHnJ0/AgNwO9bb+OhNnpFd4ilQeARzSrjTSAunt+ExnDyGuuMsXryYOnfuTKGhoTRkyBDatGmTW+AfffQR9enTh0JCQuT7J5984lSfE4j0VbsAAEAASURBVGs//vjj1LZtWwoLC6OJEyfSgQMHnOrgCwi4IoC4XK6o4FhjCWxS4V96IPtHY9lZsX6UcBL6OS0cloGt+Bvwps6GMQBXrFhBDzzwAD366KO0e/duGjduHE2fPp1Onz7tksfmzZtp5syZNHv2bNq7d698v/XWW2nr1q2O+s888wy98MIL9PLLL9P27dspOTmZrr76aioqKnLUwQcQcEUAcblcUcGxxhAoKa+knafy5CXj7CkGG3M96lqTgPIG3oR4gNb8AXhRa8MYgGyozZs3j+bPn0+9e/emF198kVJSUmjJkiUucfB5NuYWLFhAvXr1ku+TJk2S1/EFPPvHddigvOmmm6hv37705ptvUklJCb3zzjsu28RBEKhJQN2IEZerJhV89pTA1uMXqKLqskj9FiZTwHl6HepZm4C676Qey5EewdamAe2bQ8AQBmB5eTnt3LmTpkyZ4qQrf09NTXU6pr7wDGDt+lOnTnXUP3HiBGVmZjrV4aXiCRMmOOqottR7WVkZFRYWOr3UObxbj4DaB8hxuQpLEZfLer+A5mms9v+N7daaWrRo0bzGcLVlCPRvFyPiRQaKe04lcUpKFBBoKgFDGIA5OeJJp6qKkpKSnPTk72zEuSp83F19dZ27OrXbXbRoEcXExDhePAOJYl0CKfHh1Blxuaz7A2im5mr/3/juthRfzWwOl1uEADsLjbanhVMPERZRHWp6mYAhDEClc+2nZF7GrX1M1eX32udc1fekjmqTl5MLCgocr/T0dHUK7xYl4IjLhXAwFv0FNE3tzIJSOpJVTMKh0/HHvGkt4SorEhjXw/bQ8D3uO1Ycfq/pbAgDMCEhgQICAurM9mVlZdWZ5VNk2KFDzfKpYzXr83ku7uqo69Q7LxFHR0c7vdQ5vFuTgFoGVrM51qQArRtLQP1e+ovYfzHhSP/WWH5Wrz9ObBvgsut0HhWXVVodB/RvIgFDGIDBwcEy7MvatWud1OTvo0ePdjqmvowaNYpq11+zZo2jPoeTYSOwZh3ea7hhwwZHHdUW3kGgPgIcv409gk+KmFzpF0rqq4bjIOBEQDkOqRlkp5P4AgINEOjQKlw6DlWKTDJbkBauAVo4XR8BQxiALPxDDz1ES5cupWXLllFaWho9+OCDMgTM3XffLXWbM2eO9PRViv7+978nNvj+93//lw4ePCjf161bJ0PJcB1e+uWwMk8//TRxfMD9+/fTHXfcQeHh4TRr1izVDN5BwC0Bjss1KCVW1sF+HLeocNJOoFr80VZLd2oGGXBAoLEE1G9HzSY39nrUB4FAoyDgmH65ubm0cOFCysjIkGFbVq5cSR07dpQqcDzAli1/tmd5ZvC9996jxx57jP7yl79Q165diWMJjhgxwqHyH/7wB7p06RLdc889lJeXJ8+x0RgVFeWogw8g0BABDsuwQ8Rz+/5oNs0a0aGh6jhvcQJpmYWUe7GcIoIDaFCHOIvTgPpNJcCzx8u3nqZN9lzSTW0H11mXQAvhGHHZuuo3T3MOCcNewewYwnsDUaxJgIP5/nJJKsWEBdGuv1wtl4StSQJae0LglQ3H6G+rDtKkXon02h3DPLkEdUCgDoGCSxU0aOEaEhPK9MOfrqJ2sWF16uBA/QTw95vo5ymz+jnhDAiAgBsCA9rHUJSIy8U35H1nC9zUxCkQIMfyL/b/4dfQHAL8wDnAvv3ke3tKwea0h2utRwAGoPXGHBp7mYAtLlcr2SpuxF6Ga7LmSiuqaNvJC1Krsd2R/9dkw+t3dcbZf0PYf+x39KboEAagKYYRSmhNQP0x34i4XFoPha7733biApVXVlPbmFDq2jpC17JCOP0TULPIP4h9gOxchAICjSEAA7AxtFAXBOohoLI57EZcrnoI4TATUB6b7DhUOwg9CIFAYwkMFEvAkSGBlFdSQQfOFTb2ctS3OAEYgBb/AUB97xDo2CqCUuLDqKLqMm09nuudRtGK6QiopTq1dGc6BaGQXwkEibRwI7vYtp9sxD5Av7I3Q2cwAM0witBBFwTUH3X1R14XQkEI3RDIKiqlg5lFYuaPaEw35P/VzcAYXJDx9rRwanbZ4OpAfD8SgAHoR9joytwE1DIwbsTmHuemasf7tLhc0Taa4iOCm9oMrgMBJwIqIDSHoyopR1o4Jzj44pYADEC3eHASBDwnMKprAomscHQs+yKdzb/k+YWoaQkCamZ4rD2PqyWUhpI+J9A5IULGALRtP7F5mPu8U3RgCgIwAE0xjFBCDwQ4LhdvyuaCcDB6GBH9yMDx9pUBqJbs9CMdJDEyAXYmUr8p7AM08kj6X3YYgP5njh5NTEDtA0Q4GBMPchNU471/2UVlFBYUQEM6Iv1bExDiEjcE1H1HPWS4qYpTIOAgAAPQgQIfQKD5BGrG5apCXK7mAzVJC2pf6Mgu8RQSGGASraCGXgiM7tpKbj85mlVM57D9RC/Dons5YADqfoggoJEIcGqmKBGXK1/G5UJaOCONnS9lVTMzaqbGl32hbesRiA0Ppv7t1fYTm7OR9ShA48YSgAHYWGKoDwJuCHBcrlHiaZyL+qPvpjpOWYAAp3/bKjKAcFF7tSygNlT0MwEVhQD7AP0M3sDdwQA08OBBdH0SGNfDluN14+FsfQoIqfxKQKV/ayPTv0X6tW90Zh0C6r7zvQg3hO0n1hn35mgKA7A59HAtCLggoJ7Ed4m0cBfLEJfLBSJLHVL7/3h/KNK/WWro/aqsSguH7Sd+xW7ozmAAGnr4ILweCTilhTuBtHB6HCN/yqS2AmD/nz+pW68vbD+x3pg3V2MYgM0liOtBwAUB9cd+42FsyHaBxzKHsgqR/s0yg60DRdXqA7af6GAwDCACDEADDBJENB4BdSNWy3/G0wASe4OAmv3r1y4G6d+8ARRtuCWgHjyx/cQtJpy0E4ABiJ8CCPiAANLC+QCqAZtUDwAqPqQBVYDIBiLQsVU4pcSHkUwLh+0nBho5bUSFAagNd/RqcgJIC2fyAfZAvWoRCJw9MrmomRkPLkMVEGgyAXYyUr81bD9pMkbLXAgD0DJDDUX9TQA3Yn8T11d/aZmFlFNcTuHBATS4A9K/6Wt0zCsNtp+Yd2y9rRkMQG8TRXsgYCeggv4iLpc1fxJq/9+oLq0oOBC3Wmv+CvyvNbaf+J+5UXvEXcmoIwe5dU9ggEjNFBUaSAWXKmjvmXzdywsBvUsA+/+8yxOteUag5vaTTQhG7xk0i9aCAWjRgYfavicQKNLCje2WIDtCWAbf89ZTD5fKq2j7iTwpksrQoCf5IIu5CTi2nxxBNiJzj3TztIMB2Dx+uBoE3BIYj7RwbvmY9eSW47lUXlVN7WLDqEtChFnVhF46JaDuO98fyaFK8TtEAQFXBGAAuqKCYyDgJQLqRrwnPZ8KSiq81Cqa0TuBDfalNx5/pH/T+2iZT74B7WMoWmw/KSytxPYT8w2v1zSCAeg1lGgIBOoS4Bmgrq0jSEQEoR+OIStIXULmPLLRvvQ2oYdtC4A5tYRWeiXA20/UMvAGZCPS6zBpLhcMQM2HAAKYnYCaBcQ+QLOPtE2/9AsldDz7IgW0bEGj7XtAraE5tNQTARWFAPcdPY2KvmSBAaiv8YA0JiSgDEAOC3L5spgKRDE1ATX7N7hDrFiGCzK1rlBOvwTUfedHEYEgv6Rcv4JCMs0IwADUDD06tgqBkZ1tceDO5l+iY2JmCMXcBDYcsnleju/e2tyKQjtdE2gTE0bdEyPl9hOVkUbXAkM4vxOAAeh35OjQagTCRCaI4Z3ipdpYjjH36FcIj8vUY7lSSTUDY26NoZ2eCUxAFAI9D4/msuneAMzLy6PZs2dTTEyMfPHn/Pz6g+peuHCB7rvvPurZsyeFh4dThw4d6P7776eCggIn2OyZV/v1yiuvONXBFxDwFgHHfhzE5fIWUl22s/t0PhWXVVJ8RDD1axejSxkhlHUIqIcQ9krH9hPrjLunmureAJw1axbt2bOHVq9eLV/8mY3A+sq5c+eIX8899xzt27eP3njjDXndvHnz6lzy+uuvU0ZGhuM1d+7cOnVwAAS8QUDdiDk+XGlFlTeaRBs6JLDhcJaUigOAtxROICggoCWB4Z3jKUSkITxfWEaHzxdrKQr61iGBQB3K5BApLS1NGm9btmyhESNGyOOvvvoqjRo1ig4dOiRn+RyV7R/69u1LH330keNw165d6amnnqLbb7+dKisrKTDwZ5VjY2MpOTnZURcfQMBXBHomRVFSdIi8Ee84mUdjuyM8iK9Ya9nuRnvIDWXwaykL+gaB0KAAGiFyUfPWE371TI4CFBBwEND1DODmzZvlsq8y/ljqkSNHymOpqakOJRr6wMu/0dHRTsYfX3PvvfdSQkICDRs2jHj5t7rafcT0srIyKiwsdHo11DfOgwAT4O0GKi6X8hIFGXMRyCkuo31nbVtNxsPAN9fgGlgbxz5AbD8x8Cj6RnRdG4CZmZmUmJhYR3M+xuc8Kbm5ufTkk0/Sb3/7W6fqfOyDDz6gdevW0W233UYPP/wwPf300051an9ZtGiRYy8i70lMSUmpXQXfQaBeAmpWCI4g9SIy9AlOu8Wld5toSowONbQuEN48BFQw8q0nLhDnqEYBAUVAEwPw8ccfr+OAUdshY8eOHVJGPl678GZWV8dr1+PZumuvvZb69OlDf/3rX51OP/bYY3IpeeDAgdL4W7hwIT377LNOdWp/WbBggXQm4RlFfqWnp9eugu8gUC+BcWJfGP+cD2YWUWZBab31cMKYBJRhrxx+jKkFpDYbga6tI6ltTCiVV1bT1hM2D3Wz6Qh9mkbg5w1xTbu+SVfx0ivPurkrnTp1oh9//JHOnz9fp1p2djYlJSXVOV7zQFFREU2bNo0iIyPpk08+oaAg9wFZeWmZDUbur762Q0JCiF8oINAUAnHCM3RA+1jivMDsLDBzWIemNINrdEigWuT6U0v7ExD/T4cjZF2ReLKEVx/e254u7jvZNLFn3VU169KxtuaaGIC8745fDRV29uCZtm3bttHw4cNl9a1bt8pjo0ePrvdyNuSmTp0qjbXPP/+cQkMbXo7ZvXu3rMeOISgg4CsCE3u2lgbgehEsGAagryj7v92fMgopp7icwkXMxyGd4vwvAHoEATcEeB8gG4BqltpNVZyyEAFNloA95du7d285i3fnnXcSewLziz/PmDHD4QF89uxZ6tWrlzQSuV2e+ZsyZQpdvHiRXnvtNTmrx/sF+VVVZdv/8MUXXxB7E+/fv5+OHTtGS5cupUcffZTuuusuzPB5Ojio1yQC6umb94tx0GAUcxBQs3+jhMdlSGCAOZSCFqYhwDmpOTc1ZyLiXNUoIMAEdG0AsoDLly+nfv36SaOODbv+/fvTW2+9xadkqaiokCFhSkpsP+qdO3cSzxJyDMBu3bpRmzZtHC+1Z4+XgxcvXiz3AHJ7L730EvEewOeff141i3cQ8AmB/iI4MAcJLhLBgjloMIo5CPCMLhfl6GMOraCFWQjEhAUR56bmsl4sA6OAABPQZAm4Mejj4+Pp7bffrvcS3itYM8L5xIkTnb67upD3BvILBQT8TYCDA48TIUI+23OO1h/KIg7UimJsAoWlFbTzVJ5U4krsrzL2YJpYel592C5ikG4Q953ZIzuaWFOo5ikB3c8AeqoI6oGAUQjwPkAuatbIKHJDTtcEfhDL+VXCCaRLQgR1aBXuuhKOgoDGBNR9h3NVl1UiHIzGw6GL7mEA6mIYIISVCIwXXqIcDoYdB7IKEQ7G6GP/nZhR4aL2dxpdH8hvTgJ9OD5lVAiViFiA20/YZqzNqSm08pQADEBPSaEeCHiJQKvIEOK9gFywH8dLUDVqhrefqJlcNcOikSjoFgTcEuBwMCoriHpocXsBTpqeAAxA0w8xFNQjgQn2vWIb7M4DepQRMjVMQM7iFpVRmMi5iv2cDfNCDW0JXNnLFgOQ9x+jgAAMQPwGQEADAupJfJPIz1mJcDAajIB3ulSzf6O7tqJQYQSigICeCYxBOBg9D4/fZYMB6Hfk6BAEiAamxFJseBAVllbKwNBgYkwCaiZlon1mxZhaQGqrEOBwMEM62AKVY/uJVUa9fj1hANbPBmdAwGcEOCjrOHvKMDWL5LPO0LBPCBRcqqBd9liOE0WmBRQQMAKBCSoKwUEsAxthvHwpIwxAX9JF2yDghoAyGtaLvMAoxiPA2Vw4/EvX1hGUEh9uPAUgsSUJKGclDgdTWoFwMJb8EdiVhgFo5dGH7poSUFkj9p8V4WCKEA5G08FoQufKkxLBn5sAD5doRkCFg7kkjL/tJy9oJgc61p4ADEDtxwASWJRAaxGTq589HMzGwzkWpWBMtavFzN8Ge0otxP8z5hhaVWoOB6NmAbH9xKq/ApveMACtPf7QXmMC6kasZpM0Fgfde0iAw79ki/Av4cEBNKyzbVO9h5eiGghoTkA9tCgnJs0FggCaEIABqAl2dAoCNgIqLtdGMZtUgXAwhvlZqD+co7smUEggwr8YZuAgqCSAcDD4ITABGID4HYCAhgQGtI+lVhHBVCTCwewQidpRjEHgO3sA7yt7wfvXGCMGKWsScAoHg6DQNdFY6jMMQEsNN5TVGwEOB6PCMmAZWG+j41qe/JJy2n3aZqyrpTTXNXEUBPRLQK0+fItwMPodJB9LBgPQx4DRPAg0RGBSryRZ5Zu08w1VxXkdEGDnD+EDQj2SIqldbJgOJIIIINB4ApN629LC/SDCwZSUVza+AVxheAIwAA0/hFDA6ATG9UigQDETeCz7Ip3KvWh0dUwv/zdptriNV9kNd9MrDAVNSaB7YiS1jwuj8spq+uForil1hFLuCcAAdM8HZ0HA5wSiQ4NoWKd42Q+WY3yOu1kdcN5m5QAy2T6D0qwGcTEIaESAw8FMsqcw/PYgVh80GgZNu4UBqCl+dA4CNgJqOQYGoL5/ETtO5cn8zXEij/Mge05VfUsM6UCgfgKTeqvtJ1l0+bLY14BiKQIwAC013FBWrwTUhuwtx3OpuAz7cfQ6TspA5+wf7MCDAgJGJjCiS7yMZZklYlpyRiIUaxGAAWit8Ya2OiXQJSGCOrUKF7EALxPnmEXRJ4F1dkedq7D8q88BglSNIsAxLMd1T5DXfINl4EaxM0NlGIBmGEXoYHgCvB9HORVgP44+h/NEzkU6Lhx12GFH5XHWp6SQCgQ8J6CiEKjZbc+vRE2jE4ABaPQRhPymIXCVfUM2BxnmXLMo+iKgwvQM7xxP7LiDAgJmIDDRHsz8xzMFlFVYagaVoIOHBGAAeggK1UDA1wTYsIgQuWU5x+z+cwW+7g7tN5KAmiFRG+cbeTmqg4AuCSRGhdKAlFgpm/qN61JQCOV1AjAAvY4UDYJA0wgEB7Z0LC2qWHNNawlXeZtAYWkFbTtxQTarQmd4uw+0BwJaEVC/6W+QFUSrIdCkXxiAmmBHpyDgmoDyBkZaONd8tDq6UWT/qBTL8l1bC2cd4bCDAgJmIqC2n7ADWmlFlZlUgy5uCMAAdAMHp0DA3wQ4vAgX3o9zHvtx/I2/3v7UjCyWf+tFhBMGJnBF22hKjg6lS8L42yxCUaFYgwAMQGuMM7Q0CIHWUSE00L4fR4UcMYjophWzSsz8qRlZtVRmWmWhmCUJyCgE9tBG39pTHVoShMWUhgFosQGHuvonMOUKW3T+tT8hPZMeRmvX6TzKL6mgmLAgGtIxTg8iQQYQ8DoB9XDD3u7ICuJ1vLpsEAagLocFQlmZwJQ+NgMwVSRoR1YQ7X8Javl3Ys/WFBiAW6b2IwIJfEFgdNcECg1qSecKSikto8gXXaBNnRHA3UxnAwJxQKBr60jizCDlVdW0QcQERNGWgFqKVxvltZUGvYOAbwiEiRBU47q3lo1/fSDTN52gVV0RgAGoq+GAMCBAxPtxrrbPAq75CTdiLX8Tx7KL6WhWMQUFtCDloa2lPOgbBHxJQK0+rMH2E19i1k3bujcA8/LyaPbs2RQTEyNf/Dk/P98twIkTJ8o/ovyHVL1uu+02p2ua0q5TA/gCAj4koAxADsxaIWYCUbQhsOaAbR/myC6tkP1DmyFAr34kwF7uItOhWAIupPQLJX7sGV1pQUD3BuCsWbNoz549tHr1avniz2wENlTuvPNOysjIcLz+9a9/OV3S1HadGsEXEPARgUEd4ighMpiKSitp63FbAGIfdYVm3RBQS2FTr0h2UwunQMAcBOIjgokzEnHBLKA5xtSdFoHuTmp9Li0tTRp9W7ZsoREjRkhxXn31VRo1ahQdOnSIevbsWa+I4eHhlJzs+qbdnHbr7RAnQMCLBALEYzgnaV+xI53WimXgsd0TvNg6mvKEQKbYDL8n3bbaoGZkPbkOdUDAyASm9EmmLeKhc43YBzhvbGcjqwLZGyCg6xnAzZs3y2VfZfyxLiNHjpTHUlNT3aq2fPlySkhIoCuuuIIeeeQRKir62aupqe2WlZVRYWGh08utEDgJAs0gUDMcDMIyNANkEy9dK8JhcBnUIZaSRJBcFBCwAgH1sLP95AXKLS6zgsqW1VHXBmBmZiYlJtoyI9QcIT7G5+orv/71r+ndd9+l9evX01/+8hf66KOP6KabbnJUb2q7ixYtcuxF5D2JKSkpjjbxAQS8TWBMtwQKCwqQYRkOnCv0dvNorwECPAPCBcu/DYDCaVMRSIkPJ84MIuKfE3IDm2po6yijiQH4+OOPO5wzlJNG7fcdO3ZIYfl47cKzIa6Oq3q8/2/y5MnUt29fYuePDz/8kNatW0e7du1SVVxe31C7CxYsoIKCAscrPT3d0R4+gIC3CYQK429CD1tYBmWMeLsPtOeaQMGlCtp8zJYSS3lGuq6JoyBgPgK8DMxFOUGZT0NoxAQ02QN47733SsPM3RB06tSJfvzxRzp/vm42hOzsbEpKsgXLddeGOjd48GAKCgqiI0eOEH/mvYFNaTckJIT4hQIC/iLAyzGrxUwUb8h+aEr9e179JY9V+vlOeF9XiimQ7okiJqOIy4gCAlYiwNtP/r7uMG06kk0l5ZUUHqyJqWAl5Jroqsmo8t48fjVU2NmDZ9y2bdtGw4cPl9W3bt0qj40ePbqhyx3nDxw4QBUVFdSmTRt5zFvtOjrABxDwEQEOPswOIQczi2RYBl6eQfE9AeX9q/Zh+r5H9AAC+iHQKzmKUuLDxD3nEm08nEPT+rp2qNSPxJCkKQQ0WQL2VNDevXvTtGnTiJd02ROYX/x5xowZDg/gs2fPUq9evaSRyO0eO3aMFi5cSLyEfPLkSVq5ciXdcsstNGjQIBozZozs2pN2PZUR9UDAlwTiRFiGYZ1s+WcRlsGXpH9uu7SiijYctmVgwf6/n7ngk3UI8BarqY5l4Pr321uHiDk11bUByMjZm7dfv340ZcoU+erfvz+99dZbjtHgmT0OCVNSYgtaGRwcTN988w1NnTpVGon333+/vI73AAYEBDiua6hdR0V8AAGNCaj9OF/vx43YH0Px/ZEcsexVRW1iQqlfuxh/dIk+QEB3BKbYY1+yIwiC0etueLwikCZLwI2RPD4+nt5+++16L+G9gjVDZLBn7oYNG+qtr0401K6qh3cQ0JoAL78s/PIn2n7qAmUVllIiQpL4dEgcy79i/6U7ZzOfCoHGQUBjAkM6xhEHhr5wsZy2n7hAo0VUAhRzEdD9DKC5cEMbEGg8gbaxYTIWnXB+lw4hjW8BV3hKoFKk3Vtnj/+H5V9PqaGeGQnw3uPJvW1h2NgRDcV8BGAAmm9MoZEJCVzT1+bAtHJfhgm1049KO07lUV5JBcWEBdEwe0os/UgHSUDAvwSm2+87q8X2kyoODIhiKgIwAE01nFDGrASm97N54W0TSzHZRYjO76txVgb25N5JFBSA26OvOKNdYxDgYPTRoYGUJe45O0RmEBRzEcAdzlzjCW1MSqB9XDgNaB8jo/NjOcY3g8wzHCv32Za6Zgywzbj6pie0CgLGIBAc2JKUM4h6ODKG5JDSEwIwAD2hhDogoAMC1/SzGSWrsAzsk9Hg2dUckfuUl3/HdMWGd59ARqOGI3Btf/v2EywDG27sGhIYBmBDhHAeBHRCQBmAW47nSkNFJ2KZRowvfzwndZkqsiDwzAcKCIAAyYchXgbmrSdYBjbXLwJ3OXONJ7QxMQHOAsJx6XgvNnJ0eneg2fuXN7pzubZ/W+82jtZAwMAE+GFIecR/hdUHA49kXdFhANZlgiMgoFsCahYQ+3G8O0RbxfJvroh3FhseRKO7tvJu42gNBAxOwLEMLPbIwhvY4INZQ3wYgDVg4CMI6J3ANXZv4M1iGZgDtKJ4h8CXP9rC60wT2Q/g/esdpmjFPATYG5j3xvIe2e3wBjbNwMIANM1QQhErEOjYKoKuaBstn8LXIDirV4bctvxrMwDVTIdXGkYjIGASAvxQxHtjuXxlf1gyiWqWVgMGoKWHH8obkYBjGRi5gb0yfDybysGfOe3VqC5Y/vUKVDRiOgLqvrMK3sCmGVsYgKYZSihiFQLqRpx6NIfysAzc7GFXMxqcczkQwZ+bzRMNmJNAzWVgDpmEYnwCMACNP4bQwGIEOidEUJ820VQp3IH5aRyl6QQq2PvXvpQ+wx5nsemt4UoQMC8BXgbmPbJcvtpnC5lkXm2toRkMQGuMM7Q0GYHrB9pClXy6+6zJNPOvOqnHcilfLP8mRAbTcOT+9S989GY4AtfYg0IjN7Dhhs6lwDAAXWLBQRDQN4FfCAOwRQuibcIj70xeib6F1bF0X9mDP2P5V8eDBNF0Q4BDJMWJUEk5xeW0WTw8oRibAAxAY48fpLcogTYxYTSys81h4fO9WI5pys+gtKLq5+DP/RD8uSkMcY21CPAysPKU/3j3GWspb0JtYQCacFChkjUI3DDo52Xgy5dFehCURhH49mAWFZZWUpuYUBqB5d9GsUNl6xK4cVA7qTwvA5eUV1oXhAk0hwFogkGECtYkMK1vG5mz9vD5YkrLKLImhGZo/fEu2/7JG8QftJYtxXo6CgiAQIMEBneIo46twoXxV4WUlA3S0ncFGID6Hh9IBwL1EuDI/JN6Jcrzn+6BM0i9oFycyBUZDdYfypJnbrLPaLiohkMgAAK1CLQQm49vGGibBfwETmi16BjrKwxAY40XpAUBJwLX22/En+85hxydTmTcf+HUbxxGp1+7GOqeFOW+Ms6CAAg4EeBZcy6bjmRTVlGp0zl8MQ4BGIDGGStICgJ1CFzZqzVFhwZSZmEpbRUZLVA8I/DxLtsGdrWfybOrUAsEQIAJcCzSQR1iSTxDET98ohiTAAxAY44bpAYBSSAkMMDhlYdlYM9+FMeyi2nvmQIKEPv+OJwOCgiAQOMJqK0TWAZuPDu9XAEDUC8jATlAoIkE1H6cVfsyiUOboLgn8Ind+WNCj9YiAHSI+8o4CwIg4JLAtf3bUqB4iDpwrpAOn4cTmktIOj8IA1DnAwTxQKAhAsM6xVNbEcqkqKySOLQJSv0EqsWalZqxwPJv/ZxwBgQaIhAfEUwTe9qc0NT/qYauwXl9EYABqK/xgDQg0GgCHMLkF3ZnEBXapNGNWOQCzpxyNv8SRYUE0tV9kiyiNdQEAd8QuGmwzRnkM+ENzA9XKMYiAAPQWOMFaUHAJYFf2m/E34nQJvDKc4lIHlTOH9f0a0OhQQH1V8QZEACBBglcJcJQRQkntHMFpbTlBJzQGgSmswowAHU2IBAHBJpCgEOZDBZeeVXiKfzDnUjR5Ioh74/kfZJc1MyFq3o4BgIg4BkBfoi6VjxMccHqg2fM9FQLBqCeRgOygEAzCNw2vIO8esX2dEJquLog1/x0Xu6TbBcbRrxvEgUEQKD5BH45pL1s5CsRW7OwtKL5DaIFvxGAAeg31OgIBHxLYEb/NhQp9radyi2hzYgJWAf2u1tPy2O8XI7Ub3Xw4AAINInA0I5x1C0xki6JGfbPEBOwSQy1uggGoFbk0S8IeJlAeHCgI64dzwKi/EzguIj9x0Yxp/ydaZ8p/fksPoEACDSVAKeG+5X9/9Q74iELqw9NJen/62AA+p85egQBnxG4bViKbHvV/kzKLyn3WT9Ga/jdbbbZvytF2ApeAkYBARDwHgGeVQ8ObElpGYUyyLr3WkZLviSgewMwLy+PZs+eTTExMfLFn/Pz8+tlcvLkSeInElevDz74wHGdq/OvvPKK4zw+gIARCXBu2z5toqm8stoR786IenhTZnb++MDuGDNrhG2fpDfbR1sgYHUCseHBDmcQtdXC6kyMoL/uDcBZs2bRnj17aPXq1fLFn9kIrK+kpKRQRkaG0+uJJ56giIgImj59utNlr7/+ulO9uXPnOp3HFxAwGgF+sLltuG0W8L1tcAbh8VstZ0MrZLBsFbjWaOMKeUFA7wTUMvDne89REZxB9D5cUr5APUuZlpYmjb4tW7bQiBEjpKivvvoqjRo1ig4dOkQ9e/asI35AQAAlJyc7Hf/kk09o5syZFBkZ6XQ8Nja2Tl2nCvgCAgYkcL0ICv3UV2l0SKRn2pOeL5K2xxlQC++JzPuSuLCXNOf/RQEBEPA+gWGdbM4gR7OK6VPhDDJ7ZEfvd4IWvUpA1zOAmzdvlsu+yvhjzUeOHCmPpaamegRi586dcgZx3rx5derfe++9lJCQQMOGDSNe/q2urq5Tp+aBsrIyKiwsdHrVPI/PIKAHAjFhQY7lGJ4FtHI5Ioxgzv7Bht9M+/5IK/OA7iDgKwK8+qBmAeEM4ivK3m1X1wZgZmYmJSbacg3WVJuP8TlPymuvvUa9e/em0aNHO1V/8sknifcErlu3jm677TZ6+OGH6emnn3aqU/vLokWLHHsReU8iLzejgIAeCaiYgF/8eI6KRY5gq5bl9tm/SSJjQVJ0qFUxQG8Q8AuBms4gP54p8Euf6KTpBDQxAB9//HGXTho1HTN27NghteJjtQu7mbs6XrvepUuX6J133iFXs3+PPfaYXEoeOHCgNP4WLlxIzz77bO0mnL4vWLCACgoKHK/0dGvPrjjBwRddEeDlmC6tI6ikvMqyziDs/KFSv8H5Q1c/TwhjUgI1nUHU1guTqmoKtTQxAHnplff3uXv17dtX7s87f/58HdDZ2dmUlNRwIvcPP/yQSkpKaM6cOXXaqH2Al5Z5eddVf6puSEgIRUdHO73UObyDgJ4I8APSHPsenNe/P2HJRO1fyswEldQ+LozGd2+tp+GBLCBgWgJqGZidQZAZRN/DrIkTCO+741dDhZ09eMZt27ZtNHz4cFl969at8ljtJV1XbfHy7y9+8Qtq3brhm//u3bspNDSU2DEEBQTMQODmoSn0/JrDdDznIq0/nEVX9Wr4ockMeisd3tl6Sn7kP0jI/KGo4B0EfEugpjPIBzvO0LyxnX3bIVpvMgFNZgA9lZb37k2bNo3uvPNOYk9gfvHnGTNmODyAz549S7169ZJGYs12jx49Shs3bqT58+fXPCw/f/HFF8TexPv376djx47R0qVL6dFHH6W77rqLeJYPBQTMQIDTwqmQMK+JWUArld2n82jX6XwKCmhBtwy15Sq1kv7QFQS0IsCrD78Z00l2v0zcdyqr3DtXaiUn+iXStQHIA7R8+XLq168fTZkyRb769+9Pb731lmPsKioqZEgYXuqtWZYtW0bt2rWT19Q8zp+DgoJo8eLFcg8gt/fSSy8R7wF8/vnna1fFdxAwNIG5ozvJ9Gc/HM2VUfoNrUwjhH9103FZm0PiJEbB+aMR6FAVBJpN4JeD21N8RDCdzb9Eqw945rDZ7E7RQKMJtBAOFZcbfRUukAR4zyB7A/MyNe8NRAEBPRL43fJd9NW+DLpVzIQ9c/MAPYroVZlO5V6kK59b///buxP4KKo8geN/uQOGIGc4wuUg9w1yyfEZAQ/QcYfFRe5hRhFEBFxExXHAERBcGBxRWC9kVUbcURDcRQFFQBFBEJAbBAFRTiMkgImBzPu/0DHQkHRMdbqq+vc+n5BOdfWr977VdP371TvkvPlk+2BEB6kdH+to/mSGAAI5C/xt6S555sPd0rhKnCy4r11IAzdzztW5Pbh+e6AF0LnTTU4IRKfAoAt9cHRy1uPJKb5H0NvdGvx1ql2O4M/3Z5sKulWgX5tqdn3gTWY6mHXfJLq1mFFdLtffAo7qs0PlEXBAoHm1a6RJQim7PvDrazIGRjiQrSuz+OF0qrz1Rcb0TPd0qOnKMlIoBKJBoOzVRUXnBdQU6JIRDfX2Uh0JAL10tigrAr9SINAKqAGgzo/n15RRv/PSoHJJaVOzjF+rSb0Q8ITAH2/I+BK2bPsR2Xss2RNljqZCEgBG09mmrlErcEuDeKkYV8zcAk4VnZ/Lj0kD2zmrv7FVu6fDta7rc+RHc+qEQHYCvyl/tegqPDrSINpmIsjOxS3PEQC65UxQDgTCKFC4YAHREcGaXjIjZM9rJzmfpbc3fCsnzC3gyqVi5FYT8JIQQCDyAn9qn9EK+M/134p20SC5R4AA0D3ngpIgEFYBnRA51swNuOtIsize4q+pGTSgfWnVPuunE88WMgEvCQEEIi/QumZp2yUjJe28+L0PcuS1c1cCPiVz58XeCHhWIC6msAT6Ak5ftkvO+agVcKnpY7TPrHhSslgh+Y+WCZ49RxQcAb8J6MTQd19oBXzl032S9NPPfquiZ+tDAOjZU0fBEci9wB/b17BB0u6jyfLeZn/0BdTWP51zTJNOPVHCtHKSEEDAPQLdGlaUmuVKyI9nfqYvoHtOi/tXAnGRFUVBwPMCJYsVlsD0KM8s2+2LZZoWmUB2x+EkiTWtf4GWBs+fKCqAgI8EtEvGg11q2xppVw36Arrj5NIC6I7zQCkQyDeBge1qSKnihWWvuWX6rpkc2ssp1fQrmroko/Xv3o7XmnoV8XJ1KDsCvhXQmQjqVSwpySlpMmvF176tp5cqRgDopbNFWRFwQOBqc4t0sJkmRdPfP9otP3t4sfZ5ZtLnAz+cEZ10NrAAvQNEZIEAAg4LFChwlYy+KaMVUKdrOnzyJ4ePQHa5FSAAzK0Y+yPgA4EBbatJGbNY+/4TZ2T+hkOerNHZ1HPyd7PWqKbhN/5Giheh758nTySFjhoBXZ6xhVmZSEcEP2u+fJIiK0AAGFl/jo5ARAQ0WBrS6ZdWQL2V6rU0e/U+OZaUIgmlY6RXy6peKz7lRSDqBHREcKAVcN4603pvvoCSIidAABg5e46MQEQF+rSqJuVii8q3iWdl7ufeWiP4pBlNOOvjjH5EIztfZxedjygmB0cAgZAEWpklGjtcV07SzOh9nY6KFDkBAsDI2XNkBCIqEFOkoDxwYy1bhqlmGhVtTfNKmrXyazn1U5rUrhArv2tS2SvFppwIIGAERnfN6As4f+MhM4L/FCYREiAAjBA8h0XADQK6OkjDynFmctY0mfT/291QpBzLcNAM+phtJpTV9J+mU3lB07mchAAC3hFoWCVObm0Yb9cIfmz+Fl8uTemFs0EA6IWzRBkRCJOABk9P3tFATNcceefLQ7Jm74kwHcmZbNPNqvJjF2yRn34+L61qlJbOdcs7kzG5IIBAvgqM7VbPDNwqKF/sT5Q3TX9AUv4LEADmvzlHRMBVAo0TSklv0xKo6c8muHLztDA6b+HKXcdsn7+Jv29oAlda/1z1ZqIwCIQoULlUjDx44VbwpMXb5WgS08KESOfYbgSAjlGSEQLeFdCReaXNtDC6RNwrn2TcXnVbbU4kp8j4RVttsbTv4rXlrnZbESkPAgjkQmBg2+qZXVCeWLQtF69kVycECACdUCQPBDwuoCtoPHJLHVuL6WaJuO9+POu6Gj35f9sl0Yz+rRMfm7mcnesKSYEQQCBkAe2CMsm05Gs33vc2fy/Ldx4N+bXsmHcBAsC8G5IDAr4Q6NGsirSsfo2c/fmcjFu41XTQTndNvT42F4b5po+iXiie6tFICpu1RUkIIOB9gQZmENogszylJh0QciY1zfuV8kgN+BT1yImimAiEW0CXavqrGRCi38qXbDsic9ceCPchQ8r/tFk7dKy5MGga2LaGNDF9FkkIIOAfgZFdrhPtE3jI3Hn4m5mSipQ/AgSA+ePMURDwhECd+JLy0IX1Oscv3CZffXsy4uXWW796YcjoNH5dxMtDARBAwFmBEmZ98r/eUd9m+pLpg/zRjiPOHoDcLitAAHhZFjYiEL0C93SoaaZXqSCp587L0LnrRVfdiFR6bc1++YdpidTBvk/1aCh6oSAhgID/BH5bp4L0bV3Vzg34wD82yt5jyf6rpMtqRADoshNCcRCItIBOrTK1Z2O7xu7BH87Kg/+7MSITtX729QkZb/oianropjrSvla5SNNwfAQQCKPA493r237ISabbxz2vrTcT1Efuy2cYq+marAkAXXMqKAgC7hGIK15YZvZpbufbW7b9qLywam++Fk5X+xj6xnq7XujvmlSSezvWzNfjczAEEMh/gSKFCshzfZpJfMlissdMSfXgW5si8uUz/2semSMSAEbGnaMi4HoBHZ037raMfjlPf7BTVu0+li9l1kEfd//PF3bKF12mbrIZ9cuEz/lCz0EQiLhA+dhiMrNvMyliRvrrYLTnlu+JeJn8WgACQL+eWeqFgAMCd12fIL9vWlnOnU+XP835IuzzdKWZfoej3tpoFohPkrJXF5UX+jeXYoULOlATskAAAa8INK16jV2iUss7bdkuBoWE6cQRAIYJlmwR8IOAtrxNMoMvdFBIStp5uce0zL2/5fuwVE1b/gabfj8fbD1iv/3/d7/mUjEuJizHIlMEEHC3wJ0tE6Rf62p2UMjKXcfdXViPlu4qM9mre2Z79RjiqVOnJC4uTk6ePCklS5b0WOkpLgKhC+j6wCPnbbSz9es8gf/Vs5H8W9MqoWeQw55HT/0kg+asky2HTklR0w/omV5N5eYG8Tm8iqcRQMDPAqnmS+eSbYelW8OKjncD4fotwpwKfv7fQ90QcEhAV97QoExvx/5z/bfmNu0mSU45J31bVc3zB/OuI0nyh9nr7Fx/uh7xSwNaSDNzC4iEAALRLaCDQro3qhTdCGGsvetvAU+YMEHatm0rxYsXl1KlQlsBQBs1x40bJ5UqVZKYmBjp1KmTbN2aMZ1EwDIxMVH69etnW/C0FU8f//jjj4Gn+Y0AApcIaMvfFDMgo3+bjNsyf16wRfq/sla+/pXzden/08VffS89Zq62wV/NsiVk/tC2BH+XuPMnAgggEA4B1weAqamp0rNnTxkyZEjI9Z8yZYpMmzZNZsyYIevWrZP4+Hjp0qWLJCUlZebRu3dv2bhxo7z//vv2Rx9rEEhCAIErC+hyceNvry+jzWohOkpv1e7jcvP0lfLU4h2iffhCSRr4fbj9iNw24xMZ8sYGM9dXmp376+0hbaVamRKhZME+CCCAAAJ5FPBMH8BXX31VRowYkWMrnV5ctOVP9x0zZozlSUlJkQoVKsjkyZNl8ODBsn37dqlXr56sWbNGWrVqZffRx23atJEdO3ZI7dq1Q2KlD0FITOzkU4Fvjp+WJ97bZkboHbU11Lm7/r15FWle/RrbihcXUziz5ufNKOLDpp/fZrO03MyP98imC0vMFS9SUP7Qrrrc/9tajPbN1OIBAgiEW4Drtw/7AO7bt08OHz4sXbt2zXz/FC1aVDp27CirV6+2AeBnn31mb/0Ggj/dsXXr1nab7nOlAFADSf0JJH0DkRCIVoHq5pbtKwNbyjIzV9f497aKrhoy48KcXbp023XlY6VSqWJyMPGsHDATO2uH7kCKMX0JB7StLrrsnPb7IyGAAAII5K+A7waBaPCnSVv8sib9e//+/XaT7lO+fPmsT9vHui3w+qAnzYZJkybJ+PHjL/cU2xCIWoHO9SrIDbXKysKN38nn+36QDQcSZZ9pHdxpBnfoTyAVMrePE0oXN1PKlJfBHa+18/wFnuM3AggggED+CkQkANQBGjkFUtp3r0WLFr9a49KVA/TWcNZtWR8HDnLpPoHtgd+PPPKIjBo1KvCnaAtgQkJC5t88QCBaBXR0sM7bpT+ajienyPr9iXIiOdWuKVzd9O2rGFdMCpl+gyQEEEAAgcgLRCQAHDZsmPTq1Svb2levXj3b56/0pA740KQteRUrVszc7ejRo5mtgrrPkSNHMp8LPDh27FjmPoFtWX/rrWT9ISGAQPYCuorHTfWZxy97JZ5FAAEEIicQkQCwbNmyoj/hSDVq1LCjfpcuXSpNmza1h9CRxCtWrLCDQHSDDvbQyZvXrl0r119/vd3n888/t9t0yhkSAggggAACCCDgZwHX3485cOCAna5Ff587d84+1ilbkpOTM89LnTp1ZP78+fZvvbWrI4AnTpxot23ZskUGDhxo5xHUqV801a1bV26++Wa5++677UhgHQGsj7t3737FASD2hfyDAAIIIIAAAgj4QCAiLYC5cXv88cdlzpw5mS8JtOotX77cTvCsT+zcudO23gV2euihh+Ts2bMydOhQ0QmfdbTvkiVLJDY2NrCLvPHGGzJ8+PDM0cK33367nTcwcwceIIAAAggggAACPhXwzDyAbvRnHiE3nhXKhAACCCCAQPYCXL9FXH8LOPtTyLMIIIAAAggggAACuRUgAMytGPsjgAACCCCAAAIeFyAA9PgJpPgIIIAAAggggEBuBQgAcyvG/ggggAACCCCAgMcFCAA9fgIpPgIIIIAAAgggkFsBAsDcirE/AggggAACCCDgcQECQI+fQIqPAAIIIIAAAgjkVoAAMLdi7I8AAggggAACCHhcgADQ4yeQ4iOAAAIIIIAAArkVcP1ScLmtUH7un56ebg+nM4qTEEAAAQQQQMAbAoHrduA67o1SO1tKAsA8eCYlJdlXJyQk5CEXXooAAggggAACkRDQ63hcXFwkDh3xY7IWcB5Owfnz5+W7776T2NhYueqqq/KQU/BL9duJBpYHDx6UkiVLBu/AFscEsHaMMseMsM6RyLEdsHaMMseMsM6RyLEdnLLWlj8N/ipVqiQFCkRnbzhaAPPwttQ3TZUqVfKQQ84v1eCPADBnJyf2wNoJxdDywDo0Jyf2wtoJxdDywDo0Jyf2csI6Wlv+Av7RGfYGas9vBBBAAAEEEEAgCgUIAKPwpFNlBBBAAAEEEIhugYLjTIpuAvfWvmDBgtKpUycpVIg79eE+S1iHW/iX/LH+xSLcj7AOt/Av+WP9i0W4H2HtjDCDQJxxJBcEEEAAAQQQQMAzAtwC9sypoqAIIIAAAggggIAzAgSAzjiSCwIIIIAAAggg4BkBAkDPnCoKigACCCCAAAIIOCNAAOiMI7kggAACCCCAAAKeESAAdOGpev7556VGjRpSrFgxad68uaxatcqFpfR2kSZNmiQtW7a0q7iUL19e7rjjDtm5c6e3K+WR0qu9rpwzYsQIj5TYW8U8dOiQ9O3bV8qUKSPFixeXJk2ayPr1671VCQ+UNi0tTR577DH7WR0TEyM1a9aUJ554QnSFKFLeBVauXCm33XabXalDPy8WLFhwUaa6kodOYqIreai/zpixdevWi/bhj+wFCACz98n3Z+fNm2cvjGPHjpUvv/xS2rdvL7fccoscOHAg38vi5wOuWLFC7rvvPlmzZo0sXbpU9MO8a9eucvr0aT9XO+J1W7dunbzwwgvSqFGjiJfFjwVITEyUdu3aSeHChWXx4sWybds2mTp1qpQqVcqP1Y1onSZPniyzZs2SGTNmyPbt22XKlCny9NNPy7PPPhvRcvnl4PpZ3LhxY+t7uTqp97Rp0+zz+rkSHx8vXbp0scu7XW5/tgULMA1MsElEt7Rq1UqaNWsmM2fOzCxH3bp1bQuVtpyQwiNw7Ngx0ZZADQw7dOgQnoNEea7Jycn2va0t3E8++aRtmZo+fXqUqzhb/Ycfflg+/fRT7ho4y3rZ3Lp37y4VKlSQl19+OfP5Hj162FbX1157LXMbD/IuoC2A8+fPt9dBzU1b/7TlT+8ijBkzxh4gJSXFng8NzAcPHpz3g0ZBDrQAuugkp6am2ls12hKVNenfq1evzrqJxw4LnDx50uZYunRph3Mmu4CAtrh269ZNOnfuHNjEb4cFFi5cKC1atJCePXvaLzRNmzaVF1980eGjkJ0K3HDDDfLhhx/Krl27LMimTZvkk08+kVtvvRWgMAvs27dPDh8+bO/aBA5VtGhR6dixI9fKAEgIv1liIgSk/Nrl+PHjcu7cOfstJusx9VumvtlJ4RHQb5OjRo2yH+gNGjQIz0GiPNc333xTNmzYIHqrhhQ+gb1799q7B/p+fvTRR2Xt2rUyfPhw0Ytj//79w3fgKMxZW570i2OdOnVEV6bQz+4JEybIXXfdFYUa+VvlwPVQr41Zk/69f//+rJt4nI0AAWA2OJF6Spu7syYNUC7dlvV5HudNYNiwYbJ582b77T1vOfHqywkcPHhQHnjgAVmyZIkd2HS5fdjmjIAOQNAWwIkTJ9oMtQVQO8ZrlxICQGeMA7lof+3XX39d5s6dK/Xr15eNGzfaW5J6a3LAgAGB3fgdRoFLr4tcK3OHTQCYO6+w7l22bFn7TTLw7SZwsKNHjwa1Cgae43feBO6//37R22Y64qxKlSp5y4xXX1ZAR6Dqe1hHtAeStpaouXag17472oJCyrtAxYoVpV69ehdlpH2I33777Yu28UfeBUaPHi3a57JXr142s4YNG9rWJ+2rTQCYd9/sctABH5r0Wqnv+UDiWhmQCO03fQBDc8qXvYoUKWIvkjoqNWvSv9u2bZt1E4/zKKDfFLXl75133pGPPvrITuWQxyx5+RUEbrzxRvnqq69sC4m2kuiPtlL16dPHPib4uwLcr9isI4Avnc5I+6hVq1btV+TGS7ITOHPmjBQocPElVN/LTAOTnZozz+k0aRoEZr1Wah96HcTHtTJ0Y1oAQ7fKlz21706/fv3sBbJNmzZ2ygydAubee+/Nl+NHy0F0QILeunn33XftXICBVte4uDg7p1S0OORHPWNjY+XSvpUlSpSw89Rduj0/yuPnY4wcOdJeAPUW8J133mn7AOq0O/pDclZA56jTPn9Vq1a1t4B12i6dlmTQoEHOHihKc9NZA/bs2ZNZex34oV8edaCemusIYH2f16pVy/7oY533snfv3pmv4UEOAqYlhOQygeeeey7dfGNPNy2C6WZKmHTzrcZlJfR+ccx/i/TL/cyePdv7lfNADcxovXTTL9ADJfVeERctWpRuAut0M/Aj3QxQSDfBn/cq4YESnzp1yr6HTTCSbibtTzcTQaeb+VvTTZcGD5Te/UVcvnz5ZT+jze11W3jT0pr+l7/8Jd20BNr3upm+K93caXB/xVxUQuYBzCFA5mkEEEAAAQQQQMBvAhd3YPBb7agPAggggAACCCCAQJAAAWAQCRsQQAABBBBAAAF/CxAA+vv8UjsEEEAAAQQQQCBIgAAwiIQNCCCAAAIIIICAvwUIAP19fqkdAggggAACCCAQJEAAGETCBgQQQAABBBBAwN8CBID+Pr/UDgEEEEAAAQQQCBIgAAwiYQMCCCCAAAIIIOBvAQJAf59faocAAggggAACCAQJEAAGkbABAQQQQAABBBDwtwABoL/PL7VDAAEEEEAAAQSCBAgAg0jYgAACCCCAAAII+FuAANDf55faIYAAAggggAACQQIEgEEkbEAAAQQQQAABBPwtQADo7/NL7RBAAAEEEEAAgSABAsAgEjYggAACCCCAAAL+FiAA9Pf5pXYIIIAAAggggECQAAFgEAkbEEAAAQQQQAABfwsQAPr7/FI7BBBAAAEEEEAgSIAAMIiEDQgggAACCCCAgL8FCAD9fX6pHQIIIIAAAgggECRAABhEwgYEEEAAAQQQQMDfAgSA/j6/1A4BBBBAAAEEEAgSIAAMImEDAggggAACCCDgbwECQH+fX2qHAAIIIIAAAggECRAABpGwAQEEEEAAAQQQ8LcAAaC/zy+1QwB6sUfIAAAAU0lEQVQBBBBAAAEEggQIAINI2IAAAggggAACCPhbgADQ3+eX2iGAAAIIIIAAAkECBIBBJGxAAAEEEEAAAQT8LUAA6O/zS+0QQAABBBBAAIEggX8BIPUoR/QJSFIAAAAASUVORK5CYII=\" width=\"640\">"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "#basic animation\n",
    "nums = np.arange(0, 10, 0.1)\n",
    "#print(nums)\n",
    "sin, = plt.plot(nums, np.sin(nums))\n",
    "frames = np.arange(0, 2*np.pi, 0.1) #overrides the default frames behavior\n",
    "interval = 60 #200 ms between frames by default\n",
    "def update(i): #by default, i counts up in integers from 0\n",
    "    sin.set_ydata(np.sin(nums+i))\n",
    "    plt.title('index: {:^2.1f}'.format(i))\n",
    "#assignment of FuncAnimation to a variable is required for some reason in Thonny:\n",
    "progress_callback = lambda i, n: print(f'Saving frame {i} of {n}')\n",
    "animation = FuncAnimation(plt.gcf(), #get current figure or gcf\n",
    "                          update, #the function to update\n",
    "                          frames=frames, #overrides the default frames behavior of counting up i in integers from 0\n",
    "                          interval=interval,\n",
    "                          #repeat = False, #option to stop animation after one iteration\n",
    "                          repeat_delay = 2e3, #optional delay in ms; 1e3 = 1000 ms = 1 s\n",
    "                          #progress_callback=progress_callback\n",
    "                          )\n",
    "\n",
    "# ani = animation.FuncAnimation(fig, animate, range(1, len(y)), interval=dt*1000, blit=True, init_func=init)\n",
    "# from matplotlib.animation import FFMpegWriter\n",
    "writer = FFMpegWriter(fps=15, metadata=dict(artist='Stephen Shadle'), bitrate=1800)\n",
    "#ani.save(\"/home/pi/Documents/Thonny/dynamic_image.mp4\", writer=writer) #explicit directory\n",
    "animation.save(\"sine_waveA.mp4\", writer=writer) #current directory"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [],
   "source": [
    "# animation.save('sine_wave1.mp4', writer = 'ImageMagick')\n",
    "#animation.save('sine_wave.gif', writer = 'ImageMagick') #didn't work; only a few frames saved to .gif\n",
    "#animation.save('sine_wave1.gif', writer = 'ffmpeg') #didn't work (not installed?)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [],
   "source": [
    "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.7.4"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}