the-laplace-eqn.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": [
    "'''@Author: Subhadeep Bej'''"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "import matplotlib.pyplot as plt\n",
    "import numpy as np"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 127,
   "metadata": {},
   "outputs": [],
   "source": [
    "#%matplotlib notebook"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 130,
   "metadata": {},
   "outputs": [],
   "source": [
    "xvalues = np.linspace(0, 20, 20) \n",
    "yvalues = xvalues"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 131,
   "metadata": {},
   "outputs": [],
   "source": [
    "X, Y = np.meshgrid(xvalues, yvalues)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 132,
   "metadata": {},
   "outputs": [],
   "source": [
    "Phi = np.zeros_like(X)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 133,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Boundary conditions\n",
    "\n",
    "Phi[0, :] = 20                 #top boundary\n",
    "Phi[-1, :] = 40                #bottom boundary\n",
    "Phi[:, 0] = 10                 #bottom boundary\n",
    "Phi[:, -1] = 30                #bottom boundary"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 134,
   "metadata": {},
   "outputs": [],
   "source": [
    "for k in range(500):\n",
    "    for i in range(1, len(xvalues)-1):\n",
    "        for j in range(1, len(yvalues)-1):\n",
    "            Phi[i][j]  =  0.25 * (Phi[i+1][j] + Phi[i-1][j] + Phi[i][j+1] + Phi[i][j-1])\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 135,
   "metadata": {},
   "outputs": [],
   "source": [
    "colormap = plt.cm.jet\n",
    "colorinterpolation = 50"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 142,
   "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, 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\", \"jpeg\", \"pdf\", \"png\", \"ps\", \"raw\", \"svg\", \"tif\"];\n",
       "\n",
       "mpl.default_extension = \"png\";var comm_websocket_adapter = function(comm) {\n",
       "    // Create a \"websocket\"-like object which calls the given IPython comm\n",
       "    // object with the appropriate methods. Currently this is a non binary\n",
       "    // socket, so there is still some room for performance tuning.\n",
       "    var ws = {};\n",
       "\n",
       "    ws.close = function() {\n",
       "        comm.close()\n",
       "    };\n",
       "    ws.send = function(m) {\n",
       "        //console.log('sending', m);\n",
       "        comm.send(m);\n",
       "    };\n",
       "    // Register the callback with on_msg.\n",
       "    comm.on_msg(function(msg) {\n",
       "        //console.log('receiving', msg['content']['data'], msg);\n",
       "        // Pass the mpl event to the 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,iVBORw0KGgoAAAANSUhEUgAAAoAAAAHgCAYAAAA10dzkAAAgAElEQVR4XuydB9gVxfX/D10Q3xdRMah0FEFsKLZfbFFiV1CjhkRFUWKNGguKChgLKsaSGDWGWBM1ioL6V1ETFSyxETTYBQGJgp0mLSD/58x957J3392dM23v7t6zz8Mj3p05c+Y7s7sfzrQma9asWQN8sQKsACvACrACrAArwArUjAJNGABrpq25oqwAK8AKsAKsACvACggFGAC5I7ACrAArwAqwAqwAK1BjCjAA1liDc3VZAVaAFWAFWAFWgBVgAOQ+wAqwAqwAK8AKsAKsQI0pwABYYw3O1WUFWAFWgBVgBVgBVoABkPsAK8AKsAKsACvACrACNaYAA2CNNThXlxVgBVgBVoAVYAVYAQZA7gOsACvACrACrAArwArUmAIMgDXW4FxdVoAVYAVYAVaAFWAFGAC5D7ACrAArwAqwAqwAK1BjCjAA1liDc3VZAVaAFWAFWAFWgBVgAOQ+wAqwAqwAK8AKsAKsQI0pwABYYw3O1WUFWAFWgBVgBVgBVoABkPsAK8AKsAKsACvACrACNaYAA2CNNThXlxVgBVgBVoAVYAVYAQZA7gOsACvACrACrAArwArUmAIMgDXW4FxdVoAVYAVYAVaAFWAFGAC5D7ACrAArwAqwAqwAK1BjCjAA1liDc3VZAVaAFWAFWAFWgBVgAOQ+wAqwAqwAK8AKsAKsQI0pwABYYw3O1WUFWAFWgBVgBVgBVoABkPsAK8AKsAKsACvACrACNaYAA2CNNThXlxVgBVgBVoAVYAVYAQZA7gOsACvACrACrAArwArUmAIMgDXW4FxdVoAVYAVYAVaAFWAFGAC5D7ACrAArwAqwAqwAK1BjCjAA1liDc3VZAVaAFWAFWAFWgBVgAOQ+wAqwAqwAK8AKsAKsQI0pwABYYw3O1WUFWAFWgBVgBVgBVoABkPsAK8AKsAKsACvACrACNaYAA2CNNThXlxVgBVgBVoAVYAVYAQZA7gOsACvACrACrAArwArUmAIMgDXW4FxdVoAVYAVYAVaAFWAFGAC5D7ACrAArwAqwAqwAK1BjCjAA1liDc3VZAVaAFWAFWAFWgBVgAOQ+wAqwAqwAK8AKsAKsQI0pwABYYw3O1WUFWAFWgBVgBVgBVoABkPsAK8AKsAKsACvACrACNaYAA2CNNThXlxVgBVgBVoAVYAVYAQZA7gOsACvACrACrAArwArUmAIMgDXW4FxdVoAVYAVYAVaAFWAFGAC5D7ACrAArwAqwAqwAK1BjCjAA1liDc3VZAVaAFWAFWAFWgBVgAOQ+wAqwAqwAK8AKsAKsQI0pwABYYw3O1WUFWAFWgBVgBVgBVoABkPsAK8AKsAKsACvACrACNaYAA2CNNThXlxVgBVgBVoAVYAVYAQZA7gOsACvACrACrAArwArUmAIMgDXW4FxdVoAVYAVYAVaAFWAFGAC5D7ACrAArwAqwAqwAK1BjCjAA1liDc3VZAVaAFWAFWAFWgBVgAOQ+wAqwAqwAK8AKsAKsQI0pUAgAHDNmDDzyyCPwwQcfQOvWrWG33XaDa665Bnr16pXYnA8//DBceumlMHPmTOjRowdceeWVMGjQoBrrAlxdVoAVYAVYAVaAFag1BQoBgPvvvz8cc8wx0L9/f1i1ahVcfPHFMH36dHjvvfdg3XXXjWzTf/3rX7D77rvD5ZdfLqBvwoQJMHLkSHjppZdg5513rrV+wPVlBVgBVoAVYAVYgRpSoBAAGG6vr776Cjp06ACTJ0+GPfbYI7I5jz76aFi0aBE89dRT5fsIkuuvvz7cf//9NdQFuKqsACvACrACrAArUGsKFBIAZ8yYAZtvvrmIAvbt2zeyTTt37gznnHOO+COvG264AW688UaYM2dOrfUDri8rwAqwAqwAK8AK1JAChQPANWvWwGGHHQbfffcdvPjii7FN2bJlS7jrrrtg8ODB5TT33XcfnHDCCbBixYrIfPh78N4PP/wA3377LWywwQbQpEmTGuo2XFVWgBVgBVgBXQXw+7R48WLYZJNNoGnTprrZyemXL18OK1euJKePS4jfyXXWWcfaDhvIpgKFA8DTTz8dnnjiCTGXb7PNNksEwLvvvht+/vOfl9P87W9/g6FDhwI+PFHX6NGj4bLLLstmS7JXrAArwAqwArlQYO7cuYnfJ5tK4PerQ+vWsNjGSEPeH/3oRzBr1iyGQAdaZtFEoQDwzDPPhIkTJ8KUKVOgW7duiXqbDAGHI4ALFy4EtAMwC6B9ewBcdLw5lP67C8A6W3wHnes+hc4wB/rBW9AZPoUuMBu2humw/szlAG8CwMyGP28CfD4b4N2A119mscdo+PSVRlpOygqwAqxA1hTYqMGhDgCAfxf/rSv92KI9AKwPAPK/+GPU3/E3vNYHWNQCoNMeAAsWLID6+nov1cW57Wj7IgCwid1hGGQMAOB3rq6uodJePGaj1VKgEACIYXWEP1zJ+8ILL4j5f6oLF4FgKP7JJ58sJz3ggAOgXbt25EUg8kEDeA+gSW+AHzfAX28A2BL/rIGO3WdDV5gN3WAW7ABToSvMgm4wG7b79n1o8jEAvAYA+N+PAf73KsD0hdGef6GqkMF9HzYN3MhVlrxDea7EZmdZAQ0FEM58XBsDAP7BC//boR6gxQYNsIc/yr/jf/FC4Av/vQECF7UEqO/nF6rkdwnHqmwBcBQDoI8ulRmbhQDA0047DXD+3qOPPlqx9x/+Kwj3BcTruOOOg0033RRwz0C8XnnlFbFCGPf+wzmDmPeSSy7R2gamEQAK6Av86QXQuvd30LUOIbAEfgiBCIPbwVuw/oxlJfh7tQECXwP47JP4vmEKbKb5THtp2uWZ+sn5ohVgyOWeYaOALxCz8ckmrwRACX9oqwyAUSAYhMLQ3xkAbVqC87pWoBAAGLcA484774QhQ4YIzfbaay/o2rWrWPghr/Hjxwvo++STT8obQR9++OFkjXUBcH1YAD1gJmwH0yqjgOFI4DcxH+aY6GBUahcQ5sIGWUxFwiz54qpObIcVYAWyoYCM8EV5YxX9YwDMRgOzF5EKFAIAq9W2kQC4YfIwcDtYUI4CiqHhGfMqo4AIg3h9G12r/8XAoUz9pQYkYh7XYOXans+2zZOvPnVg26xAkRVIgjtKvWUEsDz0i5nkMK+c35c0/CvTAwBHACmKc5q0FGAAtFC6EQDiSwDn/+EiEPxvAwwGh4HDUUCEwPVfaxgKlvMBpU9RsBcDhphFBYemkOgDFMOyFwXGilIPi8eCs7IC3hWwhTqVg0H78u8Vc//igA8NJ9xjAFQpz/fTVIAB0ELtCgCEngAbtlgLgLhkDOcDIgxuVLkYBKOAOBS8PUwrLRCRUUA5HzDoUxj44iKACWDoGw6D7qYJQGmWZdFNnGWttfo6E44N5UYB32AXFoJSXsXCD53oXxgGU44A3gIApRnwZtcyADiNF4GYiZeTXAyAFg21FgAnA8CuJQCU4IfwhxHAhtXArTdcUF4MkhgFRAgMwlwY+OJAzxAMdeAQ0+oOMcfJWy2YqVa5Ft0slaysSyoy11QhFLhKSxBTX6yif2EA3ABg0eL0VgEzAKbVu/JbDgOgRdslAqAcDkYA3KC0Grhd3YLyljDhKGDXb+dBk+AQcBDoVFFAXSjEOisihrpgKGV0BYhoLytQkhU/LLpqprKynplqjlhnTKGpmrWz9TkqP3nuX0TEr7xFTMM9BsBq9g4uO6wAA6BFn2gEgE1alCJ+wYUg8u8Nw8AIfrgljIwC4hCwHAouzwUMA5oqCqgzVzBpEQkBCqVc1PmGQXldwmHQblaBIqt+WXT5mspa7fazhZmiNpYrXSh2EP7wqtj2BX9QLf6IgUEGwKL2ynzWiwHQot2UACiHg8VQ8BqQw8AIgbgXoIwCIhBuD2+BiAIGh4DjooAmQJgU9VOsLKZEC23A0Ef0MNys1f6YU7tZXvyk1ofTsQKmClAATce2jj0JfhXwJ6EuCv5igE/4F0jPAKjTYpzWtwIMgBYKRwKgHPqV8//kf0PDwO3gu+goIG4OLYEsbi6g7hxBk/mBDqEQJTaJGAabxlf0MK+gGNVtGR4tHmbOWhUFdKBMx0Fdu0HgC5cTuekzEfgYAHVajdOmrQADoIXilQDYHaDJZqUtAOQWMMG/B1YDY+QvKQqILjWRkKcDg5hRFR2MSiM1SBoCVgEhxUaE1rZgKE2mBYiyvKLBVtHqY/FYc1ZHCuhCmE2xJmUlQR/6IsAPr2DETxX9C6cP/T9HAG1amfO6VoAB0ELRWACUQ7/4AokYBpaLQcJRQHlEHLrU7ltchB8AwSC4USKDcaAXBXkmEcIkkAxrqjG3UGZ1BYZoL204LCokJj0qDJAWL5IcZjUBLlfVNC1bBXwV0Bd0NrzpM96LgsMg7MWkYQB01QvYjgsFGAAtVGwEgHhUuNwKprz4I7AdjADC0p6AcVFAPC8Y7wkIlP8Nw6CruYG6q4cpIOcpUhhsJpdwWE1ArEVQVD1uDJIqhfzcN4UqP96stWrrlzH0SRck/AXhLgx/SUAYypcmAI4HgDYWDbQUAI7kfQAtFMx+VgZAizbSAsDAauDgYpCoKKAEvwoIVMGg7lYx1JXDKqBTQaEqf1h/lb2E9nINhrKoakUQo6rKgGTxwAayFkVHW0Byo6YbKy7qQgE+9LY8vBvluhzmlfcQ+IK/xUX/oiCRAdBN52ArXhRgALSQNRIAg1vBIPQF5wEGVgPjMHBwS5jgimB0CbeJiQTAEAhimsT5gpjAZF6gbnRQ6kgBuBShEN3yBYZoO0twGOzKRQEci8eTs2ZYARewh9WjAp829AXhjwB2Zanj5ghigpQ3gk47AnjrrbcC/pk9e7aQY6uttoKRI0fCAQccIP5//vz5cP7558Ozzz4Lixcvhl69esGIESPgyCMxzqi+xowZI9KfddZZcOONN6ozcAqlAgyASoniE5ABMDgPMGIYOBwFlCXi73hRYbARCGJmVWTQNI10kgJzKiik2IhqBpVdQtv6hMMsAyLDIqFzcBIrBVxBXtAJHeCzhj5ZcFL0rwHsRFJCuiIPAT/++OPQrFkz6Nmzp5Dj7rvvhrFjx8K0adMEDA4YMAAWLlwIN998M2y44YZw3333wahRo+DNN9+E7bffPrGvvfHGG3DUUUdBXV0d7L333gyAVk/m2swMgBZCxgJgeCUwAiBG/0LDwFFRwA3gm8ZzAEMgiC43mieoGiKOAj38zWV0kAqFVHAzAUOqbUW7+wZDWXxWI4hJ8nB00eKlUbCsPiBPSqQLezKf1vBusD3k0G64jcJRvSDsBfPUOABGde327dsLCBw6dCi0bdtWRAiPPfbYctINNtgArr32WnE/7lqyZAn069cPbrnlFrjiiitgu+22YwB09B5hALQQci0AYrB9ZwBcBIJDwHL1b3ghiITAwJ6AcjEIRvtwAQgCYDTglaKBeFlFBF2CYJytoKYqiNMBNpWtqLbUsZ8RKAy6kUdADMvIwGjxkslQVp+wF6ymKfihjVj4C8/ro4BfEPTw70mLP4hQmMcI4Ny5c0XkTV6tWrUC/JN0rV69Gh566CE4/vjjRQSwT58+sP/++0Pz5s3hnnvugXbt2sGDDz4IJ510Erz99tvQo0ePWHNoA0HyhhtugL322osB0OE7gQHQQkwSAMp5gBHDwLgYJLwlDAKhFwjEerqK9vnaNkYX1kyAkAKtmn0irWhh2K0iwGGS1AyOmh3RMHlaYJfknhfokwW6hD8i6JXrGooK5hEAw+2Gw7ajR4+ObM7p06fDrrvuCsuXLxcRPxzmPfDAA0VaHP49+uij4emnnxYg2KZNGxg/frwYGo67HnjgARH1w2HiddZZhwHQ8BmPy8YAaCFoJQB2wynJazeDDg77BiOCMgrYcDScHAYObgxd+sdmdCQQ7xnPDZR19Q2CFMiiwFtaQEjx16CfVAsMpatFB0RKk9QyRGYB7HxBH9pNHOq1gb8oYKRG/0ov77VXaPg4TQCcDABtKQ9JTJolALAnAOhEAFeuXAmffvopLFiwAB5++GEYN24cTJ48WUQAzzzzTHj99dfhqquuEnMAJ06cKKJ6L774Imy99daNvMByd9xxR3jmmWdg2223Ffc5AmjRoBFZGQAt9IwEQDkMvGVg3l8UADYsBonaEkbO7wtDYAn+5OpgvQUimLe8SMQ1CKI906hgUt5g2+jCINVuUvublEnsT9WGQ4ZEYkNxMmcK2ET5pBMk6FPBX9xcP8yXBH/h+wZQmEcAxMhdcAhYp0Psu+++Ynj3ggsuEItD3nnnHbEgRF54H3+/7bbbGplFQBw0aJBYWCIvHFpu0qQJNG3aFFasWFFxT8cvTltSgAHQoickAmAY+sL/33A0nBwGDm4JEwQ91xBoBILUPQNV0EUBKkpkEMuh2Aq3LdV2laAQi80KGAYl4EiixUuiRrO6gD0j6IuDOGlMF/4wX0JET5iNWvwR83utAeA+++wDnTp1gnPPPRe22WYbeO+996B3bzwrtXTtt99+0KVLF7j99tsbPSm4VcycOXMqfj/hhBNgyy23hOHDh0Pfvn1r9OlyV20GQAstrQCwPDew8mQQuSVMFAQGf9ONBIq8DSeKyCprRwR1QDBNGDQFQpWP1L5hAqNU24F0WYRD6R5DokGDFiSLS9gLSqIV7ZMZXcz3C7eLQaQvDgyLDIC4Rx/u+YfAh/CG8/euvvpqmDRpkhi6xWHgjh07wnXXXQe4+hcjfLgv4P/7f/+vPE8QgRGjfmeccUbk08FDwG5fGgyAFnq6AsDwYhB0KbzSF39TLQ5RrQ6WVfUCgioIS4q+UQBKN3pHsRnV9rrlJPUfUx80+2SWwTBcFQZFzcbNWHJfsCeraQR9mDkJ/Eovz+hLBxiTtoOJKiMiMlhkAMStXP75z3/CvHnzoL6+XkT8MFInF3l8/PHHcOGFF8JLL70EuLULDv2ed955FdvCdO3aFYYMGRK7yIQB0O0LgQHQQk9tAMSygotDGqKArXt/J1YDB7eEcQ2BaC94xJz4/0BEUBkNjAM83bl/KsCiQpPKTrhdqXbD+XTLUfUnUz9UdmPu5wkOGRYNG9lhNt+AF+WqMfTZgJ8qbxgYNYZ5y3WsMQB02A3ZVEoKMABaCK0EwDDwKRaDyC1hSrBWucgjCHCmkUBZ1YqzhnWHhV2BINpRwRUVllR2otqYajsqr0l5Sf3MxheL/ptnOIyrNkcY4ztENeAuzhuv0CcLNZnvh3mj8qmif+F8UecJA0CRI4AWryLOWiUFGAAthE8EQPlCSNoOBtNEnA8cnt+HycLDu7oQGATIqoBgHDhKZyhQRQUliq1gu1PtJvUV3TIp/c6FX5RyItIUEQ51pMgqSGYJ4nT0tAI+LEg1xBt2xmTIVwV/RMgTrmQAAN8HgPV0GimUdjEA4HINm1XAFsVz1hQUYAC0EDkWANFm+ESQ4IbQEvwCJ4PARmtAbgkThDUZCXQFgWEQTBoWFtUwOUtYBXUqsKHAlMqGyoekdqfaThsIZXku/LPo95i11gHRUr5CZ7eGvaA6rsBPBZGUaCF1WJgBsND9u0iVYwC0aM1oAESDEUfCIezhRRgGDp/zW4K26H3/qJHAMPhlAgTRCRXMuIRBLI9iL9wnVD5S+5BJ2VTbrnyklqdIx4DoSMgcmHEKfCpQS9KDAnFR+an54s79le/1JHBtyJvmEDBHAHPw8FTZRQZAiwZoDIBorENjAMSfK6J9Df/f6PfSljBxcOYCApXg17DRtJRFe8WwzBgHJLqLRoLtQwEoHRCi2IvqHzplqPqXqQ8qu6p2oOb3lI4B0ZOwns06h70kaKLUJQneZH6dlb5x/iRF/1QAGMjLAEhpVE6TlgIMgBZKkwAQ7UfNAwzCX2BPQLkljGsIRHtyHmEZ7gKwV7EwRBcCpUGdfQLT2BZGF9RsYEy3rLh+Z+ODTl925a9OmZppGRI1BXOY3Cvo2UJfFHBF1d0U/NBWzBw+UUzSPcV9BkCHnZRNWSvAAGghoREAypdXEAolDG4AgFvCCPirKx35FgRB1ebPlPODwyDYaA5gDBQKP0IrhvG3RnME8UcdEIxLLyuvAhUdYFLZCvcFHdvhvLplJfVDGz90+7dLv3XLtkjPsKgnXmqApwtmSdWgRPtkflP4Ux0FxwCo19E4daYVYAC0aJ5KAGxTGvoNDgGHYS8AeuWooPwtFAUsw18DCIbnBeoOB5dAsgSX8gpGBBMjgIqIoBYIYmLd4eG8w2BSnU37X5pQ6MN/03o7zFdkaKwq4IXbSHchRzi/DvhFAVrQHnW+X1IejgA6fArZVDUVYAC0UD8RAAUVtahc9BEGwCD8yUUiW64RHuFQsCsILMFfMKK4FgSpw8JhG+L/bSKCKqhQAY4qUqXKH2x3la2oPqJjPyq/SZmqvmrrk8p+3H0fdTH1JcP54oAzU7DmQj9b4JM+pAV+OtCogr8oW1WaA/h1HUBdE/MGXbQGYMNFvA2MuYLZz8kAaNFGWgAYhD35QogZBsYtYeQlQVAOCZtGAiMBLhARjIsGNobHtSBZBlQqCGIG3eHhuDyycAp86IIRxWa43+iWEc5vUial79r6RSkjKY2vetn6xfndKeAK+NAjXehTwZvKpsp33cUfDIDu+hVb8q4AA6CFxGQAxDLCsBfcDiYKDkWeymigKQQmQRx1WNgUBDEfeZ4gJk4CBhXMUGFDZSfYJ6g2XQOhSguLfiuy6mhgW1ZcflNtffnDdmkKqKCJZqUylQ/ws4G/pNNAgp6rADHkQ5qLQDgCaNIRaysPA6BFe8cDIBoN7AUYBEAJexIAg/AXvBcBgALCDOcEqgAuasNpKY3O/EBRTkRE0CkIUgCGChc6IES1GdWndMqpFjC58NHieWqU1UZvl37Uui0fwCc11QU/HV9M5vslQaPJAhEGwFp/ejJdfwZAi+YhAaB8AcgIYCPIC+wRKO8FXxrEKGAJ8JI3i04TAuNAUCsaiEZUEKCCFlV+2f4qO+F+QrXrCwYp2lj07XJWXV1clKljw6YddMopelodsLLVQhf6ZHk6PrqGP/TBZP4fA6Btb+H8HhXIBQBOmTIFxo4dC1OnToV58+bBhAkTYODAgbGyDBkyBO6+++5G9/v06QPvvvuu+H306NFw2WWXVaTZeOONYf78+WS5tQFQAl7UqSBR96QnRAiMA8DSO2jtlzx56xezBSJhuJSuO4kGUmBHBSo6oKCylUUYpGhE7tmKhLr6uCrX1I5O25uWkYd8OgDluj6m0BcFXirf0oC/ENiVXVJAIg8BqxqP76epQC4A8KmnnoKXX34Z+vXrB0cccYQSAPHw6mXLlpV1XLVqFWy77bZw5plnCvCTADh+/Hj4xz/+UU7XrFkz2GgjDNXRrrUAeAsAdA1sA4P5G4aA5YtCmpUngiT9HgRETBcCQAFbmkPBpeLcQqDwI7RFTOxvtgtFZJOoPuYUOFHZCDY/xV4wvY7tqG6mW15SV7X1hfYYrE3l0nfdsl2mT1s3l75XE/LC9bCBPrRlUhdT+JPv47i2cBD9Q9MMgC47O9uyVSAXABisZJMmTZQAGBZl4sSJcPjhh8OsWbOgS5cu4jaCIP7+1ltvGWsYDYBoruE4OPyr3ApGgFzoCLik1cDSq3Ia2oKQEoBFDwVnEQKFRCb7Aqo+0hQYUdmwAUHMq2PfNwxK+7Y+mTwtlLYwsct5sqeALfTJGmUZ/uJgkTBHkAEwe122lj2qCQA85JBDYMWKFfDMM8+U2xoBEIeV6+vroVWrVrDzzjvDVVddBd27dyf3B2sAlC+SMBji7zFRQLyluzVMCQrXbt9iEglEG0lbxYTLkCJGRgh1ooFoKAkgVEBDhQ+VnWrDoEoHcq8NJdSpt2kZacGtS//YFk0BV8DnC/zQrgomVXUggJ1wn5AuTQBc2BWgrimtGaNSLfoBoH427wNormD2cxYeAHHOYKdOneC+++6Do446qtwiOKy8dOlS2GKLLeCLL76AK664Aj744AMxR3CDDaLfCAiR+EdeCIBoGyA4BIx3QxHAIOhJuAvPAwz/Lu/LvPhfi6HgPEAg+hgbDVQBEAVkKDBIsWMLg5hft5w0IcqFbzbvPko72djnvHYKqIDJxLoK0qJsUv0wPRYuDkqpW8QE390NthgATToH5/GlQOEBcMyYMfC73/0OPv/8c2jZsmWsjt9//z306NEDLrjgAvjNb34TmS5q4UgpYRQA4u8J8wAl8MUNA8v74ZcIcUFICfgqh4KTIDB8L+7kEEyXdHqIFM5mbmAmQNAE0kzBxQVwmZZNebO48I9SDiWNz3pSyq/VNFTY0tXHJ/ihL67hLwLqRJXjygnpxgCo20E4vU8FCg2Aa9asERG+gw8+GG644QaljgMGDICePXvCrbfeGplWLwIYAYD4U9Rwb9Sm0NIDx0PBWYRA4ZPO3oFSG5uhYbShAxO6EKRjO9jbdMuJ69Wm5SufkoYErvyklkdJ57vOFB/ynMYX5MVpUhT4YwDMc6+vad8LDYAvvPAC7L333jB9+nTo27dvYkMj3GEEcNiwYTBy5EhSp0ieA5gAgHgraTWwvB/8L/49ZkGIACjFqmBMEzeHLzgnMAyIOpHAcF4pom00EO14HRrWgUET8LEBE5PyonqvjQ+kp8HRsDa1LJt0aWhh45/vvGmDnqyPCfDJvLo++4j8Bd/B4TYizP/DLBwB9N252b6OArkAwCVLlsCMGTNEvbbffnu4/vrrBdi1b98eOnfuDBdddBF89tlncM8991TU/dhjj4WPP/4YXn311UaanHfeeYCLQzD/l19+KeYATp48WcCiXCmsErIxADZAn5gDKP8esRIYb8nIHv43PAwcvB/8e/AlaDAUjKZsIbAEeWv3CgzbLN2POC+Y+JvIHxMNxHvGIIiZKTClAwcUe+FOpGM/nNekvLhObOOH6sHw6bdu2Tbp09TIxs+ovLrA5Lp8ac8G+tCGST2yAH8xvgC/Kq8AACAASURBVDMA+upobNdEgVwAoIzkhSt4/PHHw1133QW48fPs2bMB08kL9wLs2LEj3HTTTXDyySc30uaYY44B3GD666+/Fnv/7bLLLnD55ZcDbhZNvbQAUL4QgvsBhuEuODwcvOdoQYisVxoQmHsQxAroAIAJnOnY9w1VNr5QH5hgOhO9TMrJap6g3iaQk9V6oV+20GcKfqqyqToT5/OVm0AjPQNgljtu7fmWCwDMarMYASBWJmkeYBz4hecCCjvRewMG4UtG45KiduGIXdw2MY1PEKmMBGK5LheIyHY3mh9IATgqhOjAEdVmsFPr2I97GEzKVT1YLvxSlcFQqKNQNtO6AL5gzaigFlbDNvJnCpBZBcDtAeqamXeZRasB6qfxNjDmCmY/JwOgRRvFAyAajdkKRgKgBL1GizwaHAoOEUsfDRaEVMJg/NCtTwgM+hCUmzpU7B0EsQAKQOkCEcVmVP/TLSdsw7Rc1bNg65fKftJ9X3Wy8amW87qGPtTSB/jJNqLYNgFIzTypRgAZAGv5CSXVnQGQJFN0IicAKEFQvgDDQ8TB+5ZDwSUQcw+BUXbxN9NoYBwwit9N5wdiZgrAUEGDYivYbah2w11Nt5yormpaturZcOGbqgzVfV91U5Vbi/d9QJ9v8KPaV9UtDiAZAGvxSShMnRkALZpSDYBofONSCeEj4aLATm4HI32iRAExrcZQsDQdBMGkEz6op4ZQITAO7qKigbkAQSpYZgUG0Q+f0JQFKJRa+6ynxXsj81lVMOSqApSoXFRZOv5RylDZMz1fOCIfRwBddR6240IBBkALFbUBUP5rNGkhiAC6BqeCEb8oGAxtC4O54o6Ji4IpHxBIBUHt4d+IVcSiLJuIIBXeqCChCz9Uu1F9VLesuH5u4wPl2XHlJ6UskzS+62/iUxp5VNDjywcKkMWVreszpSyKTRMAjMnDAOirY7FdEwUYAE1Ua8hjDIBRkBd1NBymM4wCCjgK7Q3oCwJVdqXE4SHhqHxxv0kbsZHCLIEgFSyDfc8GRFxClo0f1GfJpb/UMm3SpaGJjX+ugMmHD2iTAmKqsimgFrRBKZNq09HwL7rHAKhqaL6fpgIMgBZqewPAICBSIbBhGBiTh6OAYahKWs2bNG8vaTg4HuaS9wxMAj7dYWFhK2sgmDYMmpSX9AykBT95g0KL90bsEDwVSGzKTisvBcAovphoQi2bYtsk+pcAvQyAlEbnNGkpwABoobRXAAxCoIOhYB8QGA99lRtBhxeeYD5qNNAEEFUgiPcTN5TWgSgdQDKBHB37UX3ZpMwsAKH0wbX/Fs87Z01QgApdVBEpcBZli+oHxb7KluYCEHQ3VQDcB6CuOVXwxukWrQKo/ydvA2OuYPZzMgBatFFqAIg+OhgKNoXApC1iqg2BiYCYEA0kQaAvENSxG+yfWYNB9M3WJ5Pnj6HQRDX3eVSAZFIiBczi7FL9oZZhGv1D/7IwB5AB0KQH1lQeBkCL5tYCQEEcLUqllRdvNBQehrvw/SAAxsEgYUGIrGrl+b5ro3WmW8SkBYGJsBezSETkKRoIugAvXxBVDSA0BWqLZ7+ms1JBS1ckKphF2aX6RC1DZU9lpwYB8NZbbwX8g6dy4bXVVlvByJEj4YADDqhosTVr1sCBBx4IkyZNggkTJsDAgQNjewoeA3vhhRfCxIkT4ZtvvoGuXbvCr3/9azj11FN1exenj1CAAdCiWyQDIBqOOBNYBwAxbdKK4CAMBl84oW1hBAQ1LAiJgqhKIKTtExicDxgFlnG/UYeDjWHPYLWw9FU5LKwDGrogZApkuuVE9XfTsinPjgv/KOXEpfFZNxu/8pZXBUQ29VHBVJJtHb+o5VBsGgz/YjWKPAT8+OOPQ7NmzaBnz56ixe6++24YO3YsTJs2TcCgvG644QZ49tln4amnnlICIB7j+vzzz8O4ceME/D3zzDNw2mmnwcMPPwyHHXaYTa/jvBiTWoM4zpeRAt4AMAoSg7AXNSewIs/aJo1aEKIDgdRFIToQWCpfvTjENQQKe4poIKbJJQii465gyxc0ufLP6GmNyeSrri59TNMWBX5c+kOFsnCZun5Sy6HYVdlKsFFkAIzqFu3btxcQOHToUHH77bffhoMPPhjeeOMN6NixoxIA+/btC0cffTRceumlZfM77LCDiCBefvnlLntiTdpiALRo9mgARIMNmz+bRgCDMId/9xAFDMMVdWVw1MrccDSQusefLQSaAiIFBEkQiIZ0AEIXgHRsh/uxbllxz4GND6pny5WPqnJM7/usu6lPLvJRIMdFOUk2VBCVlNfEf0p5VLsqWwyAsHr1anjooYfg+OOPFxHAPn36wNKlS2HHHXeEMWPGiOhdkyZNlAB4yimnwNSpU8UQ8CabbAIvvPACHHrooSJ6+OMf/9h3Ly28fQZAiya2BsAouAuf9xtMg39PWgwSBMeIbWEE+ASGgpMgUGc+YKnYyq+lKQSiragVwqawF7eVjGx2VUSQQTDwgKQBRFmHwqT3RRr66LyvqECjY9MmrQqcVLZN60Mtl2KfYqtgADh37lyoq6srt06rVq0A/0Rd06dPh1133RWWL18Obdu2hfvuu09E6/D61a9+JcAQh3PxogDgypUrAYeB77nnHmjevDk0bdpU5D/22GNVvYXvExRgACSIFJfEGAArQK3BeuJJHwEP4oZ/I8GxNBQsh4HL0ON4PqC0W00INAVEkc/VsDAa04EAE9jRsR/uuCblZQF4XPtt8cxzVk0FKMBEMUmBsyg71PKp9in2sgKARwLUNaw7pEgcTrPofwD14xvnHDVqFIwePTrSJALbp59+CgsWLBDz9BDWJk+eDDNmzIBzzz1XRAMRDKkAeN1118Gf//xnwP926dIFpkyZAhdddJGIHO67774m1eI8AQUYAC26QyUAtg4M/cYMAYteT1wJHIRE/LscBsa/e4wChkEq7ri4OOBKAwJNYc82GiiajwojOqBGtRnuqzplFAUGZT1MNbN43jkrUQEKJBFNWZ0iQvWDCn/os8qmwlaqcwAdAaBOBDDcrAhpPXr0gNatW8Pvf/97EcGTF0YD8f933313MbQbvpYtWwb19fUC9g466KDy7ZNOOgn++9//ilXEfNkpwABooZ9TAIwDO+qWMMH8ESuC8XYwEkgdCi7B1toFG0mLQqSUlBXC0UPElQtDpD2XQ8JJ8CjLy1VEEJ22AUHM7wOmbH3SfS591EHXh1pOrwIjE210wCxsn+qPThkUmwUEwIULF1YMAes05T777AOdOnWCq6++Gr7++uuKrFtvvTXcdNNNcMghh0C3bt0amZXf1yeffLJiKxkcSp41a5ZYEcyXnQIMgBb60QEQC2mICoYjgHgrvMgjOMxbAXMBZy2jgAKEiPMB8wiBKtBTRQNF/moNDdtAmS14+QApW58snlEvcGvjT5HyUoDIpL46UBZlX8cvnbIodmsYAEeMGCFADYFv8eLF8MADDwjww0jdgAEDGrVU1BzALbfcUiwSGTRokEi/1157CXC8+eabxRAwDifjHoDXX3897wVo8myF8jAAWoiYCQBE/6PmBRKigGEIbLwSmLZJdBxsmUYCw8AZbKIsRgPRPy9Dw9UEQZuyVc9UNYFQ+uYDdFX1LsJ9CgTZ1FMHyNKEPyyLUvcaBkDc6uWf//wnzJs3TwzdbrPNNjB8+PBI+BPvzIhVwPjbnXfeCUOGDBGtO3/+fDHnD6N93377rYDAYcOGwTnnnCPy82WnAAOghX5WAIjlhod3Vfv7YR7qljAV9iu3eqQOBYfBTrVtC2WLmChYjIvGRW0ajfl1ITAOUINNr4oIUqKBXkHQBshcAJdPYHLhn8Vz3Cirz7q69NO3LQrwuPLBFvykH1SfdctzZDePcwBthoBddQ+240cBBkALXdUAiMZjTgOJAkD8LWk1cBAAg2lJ+cwgMGl/wCgYM4XAOEjThUAV7ClBL+FIOWGbMCycWRBEx1zBlk9IcuWjxbMdm9VnvX34q7JJBRuVHdP7uiAWV45OPUzKpNgn2GUANO0onM+HAgyAFqoaAaCgg9BK4CDYqaKAUauB42AwCJmijMYnhMjqm84HTAMCS1DnZoGIChAp972AoAmc2cKIK9Cy9YPyDLrylVKWyzRpaJPkLwVcXNZXxxYBmJTmTOqnWy61DIJdBkBli3KCFBVgALQQuyoAGITFOPCLXUQSHwUUUGO4P2AUBEaBFGVOYBKAuYRACuipooWZAUF0xAY0XMKVjR86z6JLn3XK5bR2ChAgiVwAFcyCBk3Kp5ZDsJ0qAA4DqGtJVrNRwkUrAepvB+AhYHMNs56TAdCihTILgGEwjFkQgslczQdMCwJL4OYuGsggGHoAXIMVA6HFG6YAWQlQpF1LKpCFDZv6Qi2PYJ8BULu1OYNHBRgALcTVA0AsKLQVDP6ks8+fakuYIPjFRQExTWAoOMsQmARnuhCoAj1VtE91X3Yj53MEpWETMLOFL5Myk54nW390n1XX/uuWX6vpCSBkJA0VxKKMm/pELZNonwHQqOU5kycFGAAthI0HwADslReBKABQgFmDM6p5gMG0QehL+nv4BeVpPiC6ELVKNwxQ1OHgNCFQBYmU+4UEQayUD5hKGwhtYNriPVETWYkAZKwFFcTiCjDxT6dMon0GQOMewBk9KMAAaCFqrgAQ60kcChagQ5wPWIKixkOyFAgsudSYLOK3hVm7L2Gw2VxHAymgl8uIIFbMBXQVCQaDHclHvSzeL5nOSgQe6zroQJhL+ENbOmUT9WAAtO4RbMChAgyAFmIaAyCWmbQSGO+nHAXEIk3nA9pAYBQAJgGYqz0DZbMngRwF8ihpRH1cbx8jK2ADLVmFQVegavFsl7Pa6Oui/CzZIEKOE5d14CupQFOfdcsnlsMA6KR3sBFHCjAAWgjpBACxfF/zAMMgGSxL1tvRUHBRITAJRikQGexeVAgU/z7QhTMbUNEtK+6ZsfFB9Ry68lFVju59n3XW9cVVeiLMuCqukR1d+IpzxLQeuuVrlMMA6K3XsGEDBbwC4Ny5c2H27NmwdOlS2GijjWCrrbaCVq1aGbiZzSw0AETfQ5tBiy98w16AugAYhjibfQHRVsKCEAE/GkPBVAikbhadBF+6kUC0FXeCiCvIcx0NTB0EsUBXoOUTjFz5mMZrxacOpv5rAItpEUb5dMErqRDTOpr4oFFWqgA4HKDO4nO7aAVA/TW8DYxRX85JJucAOGfOHLjtttvg/vvvBwTANWvW7j3XsmVL2H333cVZfkcccQQ0bdo0JzJFu+kcAAWQBcryMQwcBkgFBCZtEB0HTpQ5gVmEQBUIUgCPkkaUQxwSlr0h1YhgXkBQipMnIHT5xouDSw0gcemOsS0T6MoK/EW9TxN8YwA07iWc0YMCTgHwrLPOEgc5//SnP4VDDz0UdtppJ9h0002hdevW4iDnd955B1588UUBh82bNxdp+/fv76Fa6ZjUB0D0K2ErGFsAxPxx0Bj8PfzS0ogCRkFSFPSozg2Ogy3fcwKx3KRIoAoCKfepaVIBQSzEJgrlEq5s/NB5pF36rFMup9VXoIbgD8VhANTvIpzDnwJOAfD888+HCy64QAz3qq4nn3xSDA0feeSRqqSZvW8FgFgr04UgjQAuIJEJAArwNDsrWJacNwhUgaAqkqe6n6RLVIf2HhG0BUHM7xKs0oJB135n9m2UM8dcg5+svmn009QfzfIYAHPWTwvurlMALLhWjarnDACDQBdmZ1fDwFhGUhQwBIHBFcF4SzUULNJA5TYtlKHgqHwlOaIJwdUWMbIxsxQNFFr4HhrOGgi68MfkxeMSZk3Kr7U8poBF1UkTxCrM2vimWS4DILVBOV0aCngDwGXLlon5f23atBH1wLmBEyZMgN69e8N+++2XRt28l5EZABTwFqiuhyhgGAKjI36N9+nLAgSWIDP6+Di8lwYExoFuVCdNBQJdgZcPkEozOhhsAB918f4WynABNmBFrZYmgDUya+OjQdkMgNSG5XRpKOANAHEe4OGHHw6nnHIKLFiwALbcckto0aIFfP3113D99dfDqaeemkb9vJbhBQCpMIfpdI6Gk0poRAExi87egHGQQ4HAuMieq0igbwjUArxQpDSuk6YGggyD8e8JhkL6O9QGpuillFIawJdT+DP0gQFQt6E5vU8FvAHghhtuCJMnTxZbv4wbNw7+8Ic/wLRp0+Dhhx+GkSNHwvvvv++zXqnYTh0AsVZxZ/xSh44tAFBATmBbmHjgM4sEVhsCsT620UAqCFLnEAp7aQwLyyfGRfTNFzS58M3Fm8FX/Vz4lqaNNIFP1ssF+KEtF74b+JIqAF4DULeOeYdYtBygfjhvA2OuYPZzegNAHPr94IMPoHPnznDUUUcJEBw1apTYGqZXr15iAYjLa8qUKTB27FiYOnUqzJs3Tww3Dxw4MLaIF154Afbee+9G9xFMMVpJuVIBwDD0xQEgpqMMA4ftRf1L1nJVcDwYVg7DUo+LK7kY//XXnRdY8s/fkDAVAnXS5RIE0WmfsJQVIAy/LHzWmfJi8pnGBTiZ+mcAXJFFuaiDoS8MgKaNz/l8KOANALfZZhs46aSTYNCgQdC3b1+YNGkS7LrrrgLQDjroIJg/f77T+jz11FPw8ssvQ79+/cQeg1QA/PDDD6Gurq7sC65gbtasGcm3xgCI2Rq2eREWgn+Xm0GHfo/aEFoH5kyGgTUhULUgJB741JHAuIhbtaOBLiKBVMDTiQaagCDm0d5HMPgEuAItn2DkykfSk2+RyKcGFm41yuoCklz5YwhbscW7qJuhTwyArjoF23GhgDcAHD9+PAwePBhWr14NP/nJT+DZZ58V/o4ZMwYwWofA5utq0qQJGQC/++47aNeunZErVQHAMLzFAWA4XXjo1/FQMBUCKfMBZWNUGwLRDxcgSAU8arqyPppDw5kBQXTENwjlBQh13jymmrkAHh0/XaU1hKzE4l1oYeEXA6CrzsF2XCjgDQDROYzy4XDstttuWz714/XXXxcRN+owq0kldQCwa9eusHz5cujTpw9ccsklkcPC0ocVK1YA/pEXAmCnTp0A4BZcLtHws4cIYBLMhV9GccPAYRtR/9/IFn1vQCoAltI1Hn6tZiQwzifZzioIjKt7uO9SAY+ajkFQ4+1QRCDUqH6ukloAlnf4wwIs/GMAzFVPLLyzXgEQ1ZsxYwbMnDkT9thjD3EiCG4Ng4Dm86IAIA79YiRyhx12EFB37733iiPscG4g+hp1jR49Gi677LKIW1UEwPALiboYxBIABfSEFoRkHQJVoGYzL1BlO9hpqIBHTVc1EMSCXYKVaYTL5GXi0m+T8jlPYwUswEopp4vInyzEwk8GQGVLcYIUFfAGgN98841Y/PH8888L4Pv444+he/fuMHToUDHk+rvf/c5bNSkAGFX4IYccInx97LHHIn1LNQKIHlCjeVWKAtYaBGJ9044G6oBlVUHQNQyivTSB0If/3t5wBTNsAVQkJTICf+grAyCpxThRSgp4A8DjjjsOvvzyS7EFDG7+/PbbbwsAfOaZZ+Ccc86Bd99911sVTQHwyiuvhL/+9a/kLWqczwFERUzn9FEWg6B9y7mAaMLnopAkyHI1J5ACVVmLBlJ8Dj5QulvHBPNaLRjxBVJpw6Cvenh76+XUcJ7gL/x+NpA8VQC8HaBOzkwy8XUZQP0w3gbGQLrcZPEGgD/60Y/g6aefFvP/1ltvvTIAzpo1C7beemtYsmSJN5FMARDPJf7222/hueeeI/nmBACxJMpK4DC8JS3iSBoGdgCBYQAUcFKw4eAScMVvFZMEqhUgRtj0WXuol2DTBQxag6AviKoGDEpBeeiY9G5MTOQb+oKFu4z+MQDatz1byJQC3gAQoe/f//43bL755hUA+MYbb8D+++8POETs8kKgxPmGeG2//fbitBHc5699+/ZiL8KLLroIPvvsM7jnnntEmhtvvBFwAQjuT7hy5UoR+bv66qvFRtV4ggnlMgdAtB5YLBIHgJisGsPAUS86xd6AVAiMPkIuGrR8LwyRbZwEYC4gsASTjbfECfcxShqKz3F91zQq6AQEiwiDQaHdvs4or598pEkT9sKKZAz+0D2OAOaj29aKl94AEPf6wz35Lr/8cgGA//nPf6BLly5wzDHHwA8//AC4TYzLK25j5+OPPx7uuusuGDJkCMyePVss8sDr2muvhdtvv11AIS5OQRBESDzwwAPJbmUKAMPQVoAoIFbJFQSqIEwFXy5AUFWGCdxRbQY7dSFBECtYzchg3Fuj1sCwmsAXbAPX8Bf1j2Lyl2JtQgZAA9E4izcFvAHge++9B3vttZdYZYtDqoceeqiY94dDrLhhc48ePbxVKi3DegCIXmluBo1ZdFb1UucPot3wEHLUbxrbwpThxfFQMEMgrTczCIZ0yiIMBl0sAhhmBfbCj0hG4Q/dZACkvc84VToKeANAdB/3Abz11lvF6R8Y9cOI4Omnnw4dO3ZMp3aeS/ECgOF/aepE8pJWA0cBnocFIViM6/mARYNAoRFhSFgnnW7a4KNhGhFEG5keHg5WMutAGPWuyhokZhX4gtr5gL/wO9niu8IAaCEeZ3WugBcAXLVqFeCK2hNPPLFho2TnfmfCoDMAFF/SFmvrZLOtS86jgCWQ8TsnUAVLKjhTDQcnAWsFeBEhUOVv+GFQ+R/18NhAYK5AUFY+j0CoeutFAWMeoE1VL8p9X+DnEP7QFAMgpTE5TVoKeAFAdL5t27bwzjvviIUWRb2qAoAoZjByp7MaOJyX8v+KYWA0kcaqYNmH8jInMI8QKEDT4Hi54PPtLCKIRtOKgBURBov60g3Wyyf0BctxCNEMgLXQMfNTR28AOHDgQMA/uPiiqJcaALHmwaPhYuYAihAKMQKogjbXw8BR/wIOrQhOEwJ1AbDkfjxJJK8ATl656yoSKMArQ9HAzIFgmjCIZTEQZvuVnRb4Rb37LJVJFQDHA9S1MXd40VKA+iN5H0BzBbOf0xsA/ulPfwI8Ou0Xv/iFWAiy7rrrVqiBi0LyftkBYAgOgwAYfvHozAPUzasCyioBYAmK3AwF20AgBc5cgaAvCKTUIe5ZzFREUDqZVmSQYTBbr+g0wU/W3GH0D00yAGarS9W6N94AsGnTprHa4kbNq1evzr32qQEgKkXdDzAS2EJSqxZ/RK0QrsJQcC1CoC6s6UCjru1gr7EFQbTldHgYDaYJglIMjg6m/94uAPhJ0YoMgLjgE//gdmt44dZqI0eOhAMOOED8P267dt9994n9gRcvXgzfffedOBY26RozZgw88sgj8MEHH4jt2nbbbTe45pproFevXun3wwKW6A0AC6hVoypVDQDRE+o8wDA8hvPKWqmgkACAaMr1fECGQPWTpAuBhQPBasEgQ6G6c5qmqAb0SV8dR/2CEhQZAB9//HFo1qwZ9OzZU1T57rvvhrFjx8K0adMEDOLhC8uXLxf3cM9dCgDioRG4d3D//v0BF5defPHFMH36dMBt5sKjiqZdrZbzeQNAPHHj6KOPhlatWlXoi6duPPDAA4BnBef9cgqAKIareYBoK2k1cBQEqgAwbFOA5ZpGTWgDgElgUqThYKxn3FzGsKA6cKeTVpZjkke0k+ViEdHdfUXTqhEZDDacr3rl/YVJ8b+a4Bf1jqP4rJGmyAAYJQOexIUQOHTo0PJteWgDBQDDNr/66ivo0KEDTJ48GfbYYw8N5TlplALeABD/JTBv3jzRWMELj4DD32pzCBiVSGEhSCSshZqfshF00gpjac4xBJocy5bmwpAkSA0q7GpeILU8G6AzhcDMgyA6WG0YlA3DUJj8Ba42+KUAf1hEHgFw7ty5UFdXV24/DOqEAzvhxsXv+0MPPQR4EhdGAPv06eMEAPG4VzxeFqOAffv2ZaqzVMAbAOIcwC+++AI22qhyBcPbb78tzujFE0HyfulHADUAUBfiVLAWXkiC9lVRP8O5gGjaJhIYByRJUFV0CKwFEMQ6eosKZgkGGQrXvvqzAH4pwV9eATD8nR41apRY4Bl1IZjtuuuuYqgXt4LDOX/h41VNI4Br1qyBww47TAwdv/jii3nHh0z47xwAt99+e8BFHgh6OO7fvHnzckXxXwWzZs0CHNd/8MEHMyGAjROZAsAooNMdBqZAYdT8GIsoYBzY6EJg0pBqUoQraYsYFXRRImcuI4Eqf6L6MsXHcD6TPNKGi6FhaaumYDDcCPn/93Hj7pgV2At75nHOX7ioVCOA/wSoq9x8Q+tzt+h7gPp9AHQigDjF69NPP4UFCxbAww8/DOPGjRPDtS4igHiK2BNPPAEvvfQSbLbZZlp14cTRCjgHwMsuu0yUhP8999xzxb8C5NWyZUuxMfQRRxwB+Pe8X/YAiAoE9glM2goGk9psBxOVnwJ8lCigsJ3d+YAqcLKBQJXtMhjFbGkTfAao8wKpZa4tO3k/w7jnsCZAECuflWFinRdiVgExq5AXp22K8Icu5BEAFy5cWDEErNNN9913X+jRowfgtnDyMokAnnnmmTBx4kSYMmUKdOvWTccFTpuggHMAlGXhCiBcBLLOOusUtgFqAgCjIJEYBcSsNkPBSaCjuyhEBU1Fh0BV/ZMeUhsQFOU6WDAi/fMaFZSF5BEIqW9ZW3DMG+Al6ZIy/NUiAO6zzz7iONi77rrLCABx2Bfhb8KECYDgiPP/+HKngDcARBcxDDx+/HiYOXMmnH/++YArgnAPoI033hg23XRTd7WokqVoAERngqd/BP+O9yoXxVSmDa0ExuS2J3vo5o8CvhSigMmwFx3BKioEoha+ooGmIGgLgbkEQXS6yDBYpfdmJoqtAvwVHQBHjBgh9vxD4MN9/nC3j6uvvhomTZoEAwYMgPnz54s/b775Jpx88skimrfeeutB586dBRvghcA4aNAgOOOMM8T/n3baaWIe4aOPPlqx9199fb3YF5AvOwW8AeB//vMfazS58wAAIABJREFUwPAvNhRuDPnhhx9C9+7d4dJLL4U5c+YAbhOT96sQAEgFPtUiE9mYHoaCk+FQ77QQFQDZRgJV9kv3o32Oeh6oIGgCaGnlCdfLZUQQbacSFWQYzPvreq3/VYK/ogMgbvXyz3/+U+z+gd/9bbbZBoYPHy7gDy9cOCKniAU705133lk+MhaniOHxsXKRCa4niLqCeYrTMdOviTcARJLHI+CuvfZaQfm4KAQB8JVXXoHBgweXdwtPv8ruSswlAGL1VXMJqVCYgaHgJKCqxsIQ2bsocEUFQSoEUuAzqvdTfG0EcRpnF8c9cbkFwWCFOELo7oXq21IVwU9WrdbmAPpuUrZvp4A3AMR/AeBwL04ADQIgRv/wGBe5I7id+9XNTQNA9DE4DOx5CDgK3lTDwFTgczwULIClrvHwbvwK4PjFDK6Hg9G3pGggBZpoafIbDdSB3aQn1TUIYlmpRQXDFWMgrO5LOa70DMAfusYAmM3uUateeQNAnOeHY/+4LUwQAJ955hmxKzguLc/7ZQaAWGviZtAlCml8qSJ4lOFalQ0qFMa9WC2GggUYxkSYdLeHQVumkcA8Q2CShokwZhHZowCv6pkvFAxiZRkIVU3u/35G4C91APw3QN165vKmAavm3nFOFwp4A8Bhw4YBHtuC+/3hBE+cE4ingwwcOFAc4YLnAub9cgOAqELCVjBREKiCNxMApAJfSlHAvEAgFbQocORjSJjqX/hZpPgb9/za5JU2fYAg2q5aZFBWjIEw3dd+huCPATDdpufS1Ap4A0CEI9wB/N133xUrgjbZZBOxAgh3CX/yyScLcZCzFwAUX6kWlS2nGsJVnehBgcgoAKRCYcpRwGQ4jB9SNY0EqhaGUCGLAkZUCMQyszg3sAxwFpHEYOcvLAhiJRkG1V8o0xQZAz9ZjTSiavK7tJAjgKa9p2byeQNAqeBzzz0n5gL+8MMP0K9fP7EyuChXZgCQAmomQ8lUKIwCTPwtYhgYf6buDZgMem7nA6ogrggQqKpj3HNJAVcfedMAQfHvLdu98Vy+0BgK7dTMKPgxANo1K+f2o4B3APTjdjasVg0ABVyFNFBFASkASAU+6jBwDARGAaCAk4gFISYQaHJmsFTT5tg4KmBRgMpXJJDqY/gJo/ic9FTa5i+3j8NNpaP8zRQMSgcZCpNf+BmHvqDzHAHMxrebvSgp4BUAX3/9dbF795dffikigMHr+uuvz30b5AoAUW0KBEbBHfU3D9vCqMDM9aIQFSClGQks+eJ+lbBKU58gxyDo6LVXy1CYI+ALtzYDoKP+z2acKOANAK+66iq45JJLxJYvuCI4uKEj/h2HhvN+pQaAUfCmigBinrQXg8RCJu2cYAE8mlHAJGArQiQwixCogmTKc+0KBIUvtRgVjBK5qFCYY+BjAKS8DThNtRTwBoAIfddcc015h+9qVdBnud4AEJ3WXQgSBXyuAJBiWwptuSDEBAKTYMIUAlWAUqRIoCnQqTRSPXu2+cP2fYOgeCyzNF9QJXD4fh4AsUCwF9U8aUYAF8wCqKvT7SRr0y9aBNCuG8DChQuhzsaQuQuc07MC3gCwY8eO4qy/Ih/e7A4AsZVDZwarABCzqKKA1Ll6KjtRABj3G/6ek6FgdNV0ZXCpmuovKgVyKGnke4A6JKyzQnit7fiFNUnvIR3/o+zY5q8GCOYeBuMaVN2l3X6SCg58YbEYAN12H7Zmp4A3AMQj4D7//PNC7PcXJ3FhAVAH+KiQKYA1e0PBRYZAVd3i+rUpkJnmswXQ2Hp4HhoOlpvryKDdN4RzayjAAKghFif1roA3AMRFHwcddBB89NFH0KdPH2jRonJvu0ceecR75XwXYA6A6JnmkXClkFPlZRq5U9mxBcAoX2MAEH/2vSoYy1BFzooaCUwbAktam0USfYGg8ClFGMTyGAh9v33zaZ8BMJ/tVlSvvQHg6aefDn/5y19g7733brQIBMW88847c69p5gAwCtyoEbowTNpCoMe5gEmQYTofUAVKKqhxNRysA1AqqA0/YHkaFtbRQedFkjYIMgzqtE7x0zIAFr+N81RDbwCI5/8+8MADIgpY1IsOgKhAaI6fjwhglgDQcxTQBwSqACmLEEiJbgafP1Udo55VVb2Tnm+bvD4jgtWICjIMFvVLQK8XAyBdK07pXwFvANilSxd4+umnYcstt/RfiyqV4BUAxddCcSQcpnE1DBxliwqUUeniAFCUQ58LKD7UDreGUQGTCpAoQKOKBlJsJAFuNKTR9wvE/Kp6Mgj6f6nwMLF/jbNWAgNg1lqktv3xBoA4xDtp0iQx1NumTZtCqlyTABgHe1FDzY4gMA4AkyBJBVmm28NQwUwFgVQ7eumyDYE6dVG9MFTtq8ofd78aQ8TSFwZC01bLTz4GwPy0VS146g0At99+e5g5cyasWbMGunbt2mgRCJ4PnPcrdQCMAirK3L005gHGgaHGXEA0obsghCGw8inSnReY12igrHURQZCBMO9fhnj/0wTA2QvXgbq6JsZiLlq0BrrWL+d9AI0VzH5GbwB42WWXJdZ+1KhR2VdH4aFbAMTCFHsBmgKgDpxRgDLOnscooAC9FIeCKWBEgY9qRAJLUOw/GpgE35SHm6IfxY6tH6oyqhkVDPrGEUJVS2X/PgNg9tuoljz0BoC1IGImABCFTnseYBwAJv1uuTl0OeLDEEh+tHQhkAK9UYXbgpxt/qBPLm2F65oVEJR+MRCSH4XMJGQAzExTsCO4zGANjtFW6cKig2cEx7mBJ4qMHTsWpk6dCvPmzYMJEybAwIEDY73GPQZvvfVWeOutt2DFihWw1VZbwejRo2G//fYr58H/D0cp8fi6+fPnk9XwDoDoiclCkCgQsxkGjgO7KJsOooBYnMuhYFWESAVKqgUTFOhwGQlU1acRuOQgGqhbJ9VDSmkTlY2k+1mDQfGqyPMxdTaNkaO8DIA5aqwacNUpAPbu3RsuvfRSOPLII6Fly5ax8n388cdw/fXXA64UvvDCC5UyP/XUU/Dyyy9Dv3794IgjjlAC4Nlnnw2bbLKJ2IOwXbt2YiHKddddB6+99hrg3ES8EADHjx8P//jHP8rlN2vWDDbaKGoMNNrFeADE9OFtX8L/j2kUm0G7BMA4iIuMzEXUVwf2dCAwYkVwEgAKUIiJAiZBhAoIGAIbt7kKfOMeXJXWqgfeNn/QvktbsfVNeZNplX7yPgMhVan00jEApqc1l6RWwCkAPvfcczB8+HCYMWMG/PSnP4Udd9xRgNg666wD3333Hbz33nvw0ksvif+eccYZMGLECO1DpjFiqIoARlUbo4BHH300jBw5sgyAEydOFFFC06sqAIjOujzJwwYA46BSBwDRhkMITN4IOvmEiqxAYBLIhvuqDuCo6hf1HJhCoE4dfIGkjVam74QsRgbDdWEwNG1d+3wMgPYasgV3CjgFQOnWK6+8An//+98Bh25nz54Ny5Ytgw033FBE33AY9pe//KWIzJlcJgCIx9LhSuQLLrhAgCdeGAHEYeX6+npo1aoV7LzzznDVVVdB9+7dY93C4WT8Iy8EwE6dOgHALThoGcrnKQIYBYACokLFUyN2sat0I2Sg2owDw1jf42chpDkUXIKW+MUTFBiiABllOFgHoChlBltTFwQp9fYJcbr1U71XXNuLKi8PICj9ZiBU9Rg399e0B1i0CKBdN/C6slYGJngVsJt2K7IVLwDoUzATAETQu/rqq+H999+HDh1Kw644rLx06VLYYost4IsvvoArrrgCPvjgA3j33Xdhgw2iqShq3mCprhkEwDgIs5kLGBfZ0wFDR9vCCEAyGApWgZUKjigwRAGMvEEg6kapuy8QpGiq+97xYTPvMIj+MxDq9qTk9Ah+8koTAN9euDGsV9fUuDKLF/0A29Z/4RVWjZ3jjE4UKDwA3n///XDSSSfBo48+Cvvuu2+saN9//z306NFDRAl/85vfRKbzHwHEYg22gsFsNtu3+BgGjgNQ/F1jRTAmj4sCJkGg6uPu68xg2XFU5Zdk+Ib0EFNsqcA2Ekw0F4dUGwJN6kgRmKovxVZSmjxFBYP1YCA0a/kg+DEAmmnIufwqUGgAxGHoE044AR566CHSmcQDBgyAnj17ihXElMv9HMAIAMSfXK0EjoOyLAwDC4jNx1AwFYQoYFFtCCxBVTp7BroEOIq2lGc4mMaHzTgf8gqDDIbqXhUFfgyAat04RfoKFBYAMfJ34oknAv43acsYKTlG9zACOGzYsPJCEVVzVA0A0TGXC0Gi7Akgi1BAZxg4Djhjy9MHQAEVBnsDUmBEBUaUIVEKVOQRAqkQHAtAkLwgR/XsyfsUfam2fNosMgiG61arEcMk8GMA1H0KOX0aCuQCAJcsWSJWFuOFC0lwCxnc4qV9+/bQuXNnuOiii+Czzz6De+65R6RB6DvuuOPgpptugsMPP7ysY+vWrcWiD7zOO+88OOSQQ0T+L7/8UswBnDx5MkyfPl1sT0O59AAQLYYXhoS3gYlKExEBtAU2m3mAcVCnC4aacwGxWNcLQlQQqAJAKgSpIIUKgCp/w31WVW7j9HqRQGr98wiBulpT3heqNEWICsbVschQSAE/BkBV7+f71VAgFwD4wgsvCOALX8cffzzcddddMGTIELHaGNPhtddeewmYi0uPvx9zzDFilfLXX38t9v7bZZdd4PLLL4c+ffqQ26HwAIhKpB0FFGVGRwJN5gJSPuR5mg9IqU+wA6cBgQyC5FcGOWGRYTAsQl7hUAf8GADJXZ8TpqiAVwDE7VcwcocRNvx78Npjjz1SrKafogoFgChRFhaDJAAg3jKBQBUEqe8nR8ZcDQWXmsDtwhBdYJRPCiX6GX6qKDokPYmqdqA+xa7shMvzZVdVr1qCwSgtsgaIJuDHAKjq5Xy/Ggp4A8BXX30VBg8eDHPmzIHwaXO4lcvq1aurUV+nZdoDILpjcBqIL1ijAiCWrzvk62hz6CQINN0WhgJJKiCiwA8VIPIMgaglRYu4B5GqEeVBdmkrWJ4vu6Q6ZfTUEYrvvtKkCYg28If1T3MbmCkLt4C2dc2MZV+yaDXsUf8RbwNjrGD2M3oDwO22207ssYfn7Xbs2LHRmb9yLl72JYr3MDUARBfCK4FtIZA6DxDLydAwcBIACpAz3BtQBYEqAKSCDxUefECgqo5RPZ1S76h8NhBo4mfSe4Sque67yJddqh+1Hhmk6hSVThcabcFP+sAAaNNqnNe1At4AcN1114W3335bbKtS1MsPAKJaEecGpwGAtlApG1o3OuhwQQhDoPpp0wUXUwikQnFa8KZbb7WSa1P4tK3yg0FQpVB27hcZAMeMGQOPPPKIOFABF1zutttucM0110CvXr1EA+A8/W7dukU2xoMPPgg/+9nPYhsKD3HAY2Zxbj9OJ8OjXTEPLuLky1wBbwD4k5/8RGyqvP/++5t7l/GcuQZA1JYaBYyKAMblN/k9DgDRVoYWhAi4JOyZR4l+UYEhK5FAat2jHlmKHkWAwJJGbra3MX31MQyaKpdOviIDIH7rcXFl//79YdWqVXDxxReLXTXee+89wIAQTvv66quvKoS+/fbb4dprr4X58+dD27ZtIxth5syZsNNOO8HQoUPh5z//udjJA4EQy5Ene6XTesUrxRsATpgwAS655BI4//zzYeutt4YWLVpUqLfNNtvkXs2aAUABYhHNpRvpS4LDnEQBqSBEgR4qLFAhUBdAqOUHW54CwHEPNkUTBkF3r0WGQXdaurJUZAAMa4Swh4CGUbu4RZ+4rVu/fv3gL3/5S6zECJXID/fee6+rZmA7DQp4A8CmTRufQYiLP3BBCC8CCfa/FPYCjAMvnYUZGYwCYrVM9gakgJIKjiggRAEeVTmyp/iCQIoW4bclpe6+INDE37SgMq4cahv7/CoxDPpUl267lgAQdwDZfPPNRRSwb9++jUSaOnUq7LjjjvDyyy+L4eKoC4d7MeKHo4kvvfQSTJs2TQwj496/lAMe6C1Tmym9ASCu/k26qJstZ7lZqh4BRHGoK3ejYE8HALEsF1FA3TJFuW5PCJF9SvVhVt9Xb5xcVAgsgZi6/r5AUNU2uu8N1/aiyk+jDEq9GQYpKvlJk0cAnDt3LtTV1ZUFadWqFeCfpAsDPYcddhh899138OKLL0YmPe2008TevThEHHfh0DAuIm3Tpo04rAH3A540aRKMGDECnn/+edhzzz39NFSNWPUGgLWgX6oAiIK6XgiCNqnzAOMAMM6Gye8GcwGxGJO9AUsAkzxfS32fBkAMgdFvA4ouqveIqo1U+cP3XdvLMgiKZ4C3ldHtIlbp0wTAxxf2h3Xrmhv7+/2iVXBI/RuN8o8aNQpGjx6daPf000+HJ554QkTtNttss0Zply1bJsDu0ksvhXPPPTfW1ueffw6bbrqpmPt33333ldMdeuihYl4hnvrFl7kCTgHwsccegwMOOECM1+Pfky5swLxfhQVAbBhqZNEE9OKigHHlCvgsdhSQAqTyecnScLD0yTQa6AICdbSjvnNqDQQZBqk9wy5dHgFQNwJ45plnwsSJE8VJW3GrfnE+Hy7qwCNc8SSuuGvlypUC9BA6cU2BvHBFMMIlDh/zZa6AUwDEeX8YssWJn1FzAKWbPAcw2GBRcwDxvuOtYOJATXdINmdzAcWHLWFvQAo8qGCAAj9U0FGVlQYEUjSJeuVQdIh7VVH1Ub3qqPqp7KwF23RW9br2m1q/uHQcGbRVMDp/HgFw4cKFFUPAccrgsC/CHy4AxaFdnP8Xd+FxrRtuuCGMHz9eKTTOD+zRo0fFIpBBgwaJrWaCUUGlIU7QSAGnAFhr+rqJAKJqjheCoEnb+Xqxq3JjWll3RbBJFFDUiyOBJhBoAnWmUFJtEDT1O+n95cNmLHxVeSuZsF8Mg+6+bEUGQJzTh0D26KOPlvf+Q+VwEQfCmrxwcQgeEvHkk09GbhO35ZZbAu4piJCHFwLl0UcfDX/84x/LcwDPPvtsAZk//vGP3TVODVpiALRodH0AxMIiIntZBEB01cUwMNrRhcOU5wJS4IgCABTwoUa6KOWVmoh2drDs5lS7punX5qPNj4x6/KgaqR5d3bqq7FH6CcUGNY0P/6llx8Ipzxm0krDIAIgje1HXnXfeCUOGDCnfwgUcOASMC0XjdgsJ57njjjsEFP73v/8VcIknjOEiE77sFPAKgN9//73YA+jTTz8FHMsPXr/+9a/tPM9A7kwAoA6ouQIx38PAcXWSbZ7RKGAJENTgQwUcKgAUEQJRS6pOSa8Cqoa6rxNfdqP8SLMsqg4cFaQqVZmuyABopgjnqqYC3gAQ9+s58MADYenSpYAg2L59e/j666/Fcm6cI/jJJ59Us95Oys4dAGKtfW0HE2c76fekex6igALQPM8HZAisfLQoQJz0MLqAwFKb+JnL58sug6CTV3TmjDAAZq5JatohbwCIkzxxnP/WW2+Fdu3aiXOBcXXwL3/5SzjrrLPg8MMPz73wqQMgKmazFYwJpFVrGJijgOL50AGMrEYCqVCcVxDUaSMXL720y6P6zFFBtVIMgGqNOEV6CngDQIS+1157TYzX49//9a9/Qe/evcVvxx9/vDgwOu9XMgBi7aLm+1nMAdQBQExruxAkDsLiV+3rz/dLeTFIkaOApebK5pzAokOgLqy7ePcxCLpQMV0baQLgvQv3hTZ1lUew6tR26aL/wbH1/wDqKmAd25w2Gwp4A0Dc2wf36MEoIELg73//e9hvv/0E+OHZfzg0nPerZgEwDi6TIoxJ9xgClY+CzsfeNwTawE5WhoRt6qBqLJ22Utmi3E+7PIpPMg1HBSvVYgDU6T2c1rcC3gDwpz/9qVj5M3jwYDjllFPEGX648ANX/+DxMBgJzPvlFwBRHeJegLqROp3FILrbwZiAXkEBkBr1os5z0/nQFxkCUVeqZpR3jI6uFHtl+PE07zDOB1/10KlzrG+8elhIwwDoojexDVcKeAPAN998ExYvXiz27fnqq6/EsC/u3N2zZ0/AJd7bbrutqzpUzU5NAKAuXJoAYFKeuPJlqxuuCBZwlrMFISWgpC9myDIEUuE46eHOAwTqtpmLl5lOH3FRno6NWo8IMgDq9BZO61sBbwDo2/Es2K9pAMQG0N0OxhQOM7oimPpxpwx76sCMzge+6BCIbaCjneq9oaOtylb4vk/bUb6kXZ62HjUYFWQA1O0lnN6nAgyAFupWBQDRX+pK4DhI0xkCTorApbUYJMNRQIbAtQ+QKXBQAFn1mLqEQGqbqnyKu2+qU17K0/Gz1iKCDIA6vYPT+lbAGwB+8803MHLkSHj++efhyy+/hB9++KGiLt9++63vunm37w4A0VXicXAuANAkEqc7FzBpXp8ugHoEQPGxT2EouAQV7jaJ1oWUNCKBuj6FH1CKPqqH2iUI+gQ1n7azAp6qtgrerxUQZADU6RWc1rcC3gDwgAMOgJkzZ8LQoUNh4403hvAxMTgnMO8XA2CKw8AeIdAWAKngQwUcHYjRAQmGQLM3jo7GuiX4tM0gqNsa/tOnCYA3L/wZtLbYBmbZov/BGfUP8TYw/rtF1UrwBoDrrbeeWPRRhMUeca2TKQCMAyTdeXpx0TndCKBJlDEpjwUAYtbWGyYvnmAIjH8HmUCKSZ6KiBAhWqp6a+qAtMqWbX2qbT+qfN91UtU5FlALPDeQAdC0V3A+Hwp4A8D+/fvDH/7wB9hll118+J0JmwyADc2gC5mmcGgBgSoAFFE8y6HgokYCqfVyDRnUiGnSy8AlBNroQH1hVQPKqlEmRY8iDgszAFJantOkpYA3AHzjjTfgwgsvFPMA+/btK46BC151dXVp1dFbOf4BEF233AsQTegCGkcBI/sM5UNJSVMCCfV8QEynAzDUskscrXdiiA386PgVFp6qk+oh19FRZctGC4rtNOy7hnVqvUzSFQkEGQBNegDn8aWANwD8+OOP4ec//7nYADp4rVmzRswHXL16ta86pWaXAbBBapPVwJjVZDFIiV7ir4R9ATGTKhKYVhSQIVDvMXUBgq4hMA1Qs4FnPYXXpq5GmRRfiwCCDICUlnaX5tlnn4U2bdrA//3f/wmjt912G4wbNw769OkjTifDY2pr+fIGgDvttBM0b94czjrrrMhFIHvuuWfudTcDQKw29TzguLQOtoIxAbBE8IppTpPVwEm+VRkAqR99ykdUB2qo8EIpN9hSeYkE6gCz6sVC1VJlJ3hfV3cd29Q+p2tTld53nVTlR93POwQyAJq0unmebbbZBsaMGQMHHXQQvPvuu7DDDjuIE8leeOEFMTJ5xx13mBsvQE5vAIjUjdE/PAe4qFfmADAOjnSHgJMALAuLQaoMgdQPIyVdFiCwJKfekDClbnHPvU3eLIOgbb0o78k0ygj7UY0yVVrkFQQZAFUt6/Z+27Zt4Z133oGuXbvCZZddBm+//TY88sgjMHXqVDj44INh3rx5bgvMmTVvALjHHnuI+X/77rtvziShu5sbAMQq6UJgGvMAk0BTdS8pGinquya2IVXDwAIyHCwIoUZuGALpz5xMqaNZknWOBtK1ZxCkaxWXMk0AvHLhSbBOXUtjp5cvWgkX14/L9TYw7du3F7uR4JDv7rvvDr/85S/hV7/6FcyePVv8tnTpUmN9ipDRGwA+9NBDMHr0aDj//PNh6623brQIBEOzeb8YAEMtqAuZKshLGj4uha3iL89zAelwRzu7lwo0OsBi8sHWjQRSdYhqKBP/gnaomlHeMzq6UuzZ6JIV+3F+2LYbtX7UdHmKBjIAUlvVTbpDDjlEHELx4x//GEaNGgWffPIJbLbZZoBzA0877TTAtQq1fHkDwKZNmzbSFRd/8CIQlMXTHMAkKNKFsyT4cjkMbAOBFgCIxaoigaooIPUjT/lg6sCMDqxQyg4/qHmCwFIb0FZUq170OrqqbMn7JvpTbadZRtinNOqlo0NeIJABUKdV7dNipO+UU06BuXPnivUIw4YNE0bx76tWrYI//vGP9oXk2II3AJwzZ06iLF26dMmxbCXX3UYA0aLlcXAuATAJzLICgEn1xXsZiQLSQZEOMjqwovuxNgFAah2jHnpd/6Jt0LVLevHo6KrzAnNRx6TyfNuPK7ta5Ub2gRxsIM0AqPPUcFrfCngDQN+OZ8F+OgCINfW4F2AS6OnOAxTQFdMyPlYDWwIgZs9SFLAEUXSQ0YEV3Q91LUMgtoOOttR3kW4bUO0G06VRRtivapQZC6QZh0AGQJNebZ6nZcuW8Nlnn8FGG1V+mL799lv40Y9+BCtXrjQ3XoCcXgEQzwK+8cYb4f333xd7//Xu3VuEXnv06FEA6dKKAMYAIP7cpHJz7bKoURE6k736XAJgEmja3LNYDEIBQAFlKS4IqWUILNWdNmcyFgA0AFr1EvIBgS7qqPLbVkOV/aj71SgzjxDIAGjSu8zz4FS0+fPnQ4cOlaNrn3/+OXTv3h2WL19ubrwAOb0B4NNPPw2HHnoobLfddmITRpz798orr4hl2I8//jgMGDDAqXxTpkyBsWPHiuXduLR7woQJMHDgwMQyJk+eDL/5zW/E/kCbbLIJXHDBBWK+APWqagRQFwAxvW50Lg/zALFeOZkLqPPxp0YCdSFF90OddiRQRyOGwOQ3lW5bU997SemqUWYkkGY0EsgA6KKXqW3ccsstItGZZ54p9gHE7WDkhYdQ4LcfF4Agj9Ty5Q0At99+e9hvv/3g6quvrtAXj4d75pln4N///rdT3Z966il4+eWXoV+/fnDEEUcoAXDWrFliI8iTTz5ZLAvHvLgq6P777xf5KVfhARBFcBkFNB0GTvJDBYACfOO3hMHbqmFgASWKKCAVXKgfSCoAYrkMgdFPq46GquddV2OVPXmf2h+o9sLpfNuPhC/LKK5pXYP5srooJE0AvGDhudCqrpWxnCsWrYBr63+Xy21gOnXqJOqNw78dO3aE4KJUHBbGfQEVgbOqAAAgAElEQVQvv/xy2G233Yz1KUJGbwC4zjrrwPTp02HzzTev0Omjjz4C3ALGZ+gVh5tVEcDhw4fDY489Joan5YXRP/wXwb/+9S9S2zIAxshkMtysgjyPW8JQINAVAFJBsZTOz3xAHR9kC+cxEqiroeqhzysEmrS3SgvK/WrAZ9YhkAGQ0nPcpcG9//A7v/7667szWiBL3gAQCfz666+Hn/3sZxVyPfjgg3DeeefBp59+6k1GCgDiRtUYpbzpppvKfiA0HnXUUWJzyBYtGs+vW7FiBeAfeSEAlv6lgeHm1hH1idruBZPF/R61Cjghvc4cQDRjAmYuI4C2kGcyJC1bJUNRQJ0Pcq1DoI5WSS8UHR1VL6a8gmA1gKwaZWYZAhkAVU8X309TAW8A+Nvf/hZuuOEGwCFfDLMilOGO3Ndccw2ce+65cMkll3irJwUAt9hiCxgyZAiMGDGi7AfOUcT5ijhBFMPG4Qs3tsbjZBpfNQiAKEKaZwOr4FHpT/IwMGZXDQVTooBUYNH5MFLhxQRMdPwoSax3ZJx8VnTLCT9jtvlL7UKPqKpeTiZaq2xS+w7FTjIM2y200S3fRdvplplVCGQAtGlJWl6cy0+9rr32WmrSQqbzBoC46ANXAP/ud78TQIUXLrTAk0HwMGaENF8XFQBPOOEEuOiii8pu4DxA3DEcF5HgEvHwlfsIIFYojYUgSeWoQM5mnmAGVgTrfMSpH0YdcNEFE6oPwWeBIXCtGrp6U955Jm1CsVsBRVWYp5dGveJ0yMqcQAZA3Z6qnx6HfSkXcgIuHq3lyxsABkVdvHix+N/11lsvFa0pAGgyBBx2PrNzANFR3c2aTcDLJAJoA4CqvJYAiObTjALqwSI9eqULJSYf5jxDYEl3up6ql5au3ip78r5Ju1Btp1lGtcFTtHdGVgUzAOr2UE7vU4FUANBnBaJsUwAQF4HgdjTvvfde2cSpp54Kb731Vn4WgaDnruYBFgEAk8BXtnJKcwGpcKfzkdeBFl0o0fFDSmkKgVRt4t4bJr5G2dLRU/UO09VbZS9NQHOlJ7VOtu2vU04FfGYAAhkATVvPPh/uB4hssPHGcXPw7cvImwWnAIiLKqhDu663gVmyZAnMmDFD6I9+4AKUvffeG9q3bw+dO3cWQ724JPyee+4RaeQ2MLgFDG4Fgyt/cRWw2TYwOC+wc0Tbu1oEgqZjbLkCQCwiDgJNF1+YLDpJ8kMqbOoP5lcAICZRRQHFB8zRtjC6H0MdaNGFEhMQqBYE6uoWD5PuIoFYhq7m1A+GSdtQbacJmxVAVo1h6CpDIAOgbq+0S49T0XAfwOuuu05sZyPeG+3aiYWouD6Byit2XmQ3t1MADC6QwG1ecDPGPn36wK677ioUePXVV8Wmy7jfHjaKy+uFF14QwBe+jj/+eLjrrrvEgg88GBrTyQs3gzznnHPKG0FjVNBsI+g4AEwAtyKsBMbqZW0YWOWTIwikAKAOpOh84BkC1z7lOrolvW90NFW9txgCVQq5bz96idUdDk4TAE9aeCW0rFtHR5qKtCsXLYdx9Rfnch9AWRFcbHrbbbfBqFGjygdS4Fx/XKSKI364F2AtX04BMCjkSSedJFbShgXGhpg7dy7ccccdudd97RzAHAGgAKAY6bMyDIzumfgiq5XCXEABd4QoIBUCdUBGF1Z0gUTHl7WSm60OpuqTDG9uVrXq6qp6genqrrIn75u0D9V2mmVUoyzR36oYBWQA1O2Jduk33XRT+OMf/9joVDDc8u2MM84Qo4K1fHkDwPr6enjzzTcbbQSNx6/suOOO5XBsnsVnAGxoPY4CJnZj6gebmq4ETfShSxMQ0fElWHnTIWHT8lxDhI6ulHeXifYUu7Z6ZaUM1+1HqVc1IbDIAIijeo888gh88MEH0Lp1a7H9G2771qtXr0bNgkOzBx54IEyaNEl5aAOmxdHF22+/Hb777jvYeeedBdRttdVWyubGAyn+85//AG77Frw+/PBDcUztsmXLlDaKnMAbAOI2KtghcKuV4HXnnXeKsfcvvvgi97pmGgBRXd2VwJinCPMAk+oue12KcwFLwEaLUlHTMQRGvz509Et6ATEIrlXHlaaqF35a5ZShswqRwCID4P777w/HHHMM9O/fH1atWgUXX3yxOA0MF1quu+66Fc2PewQ/++yzgEe4qk7tQoi88sorxVQuBLkrrrhCbN+CEKfaWQR9wa3dsLzgdfbZZwPu+/v666+rumWh73sDQDwDGDdOxqHgXXbZRYiIcwBx6HfkyJECAvN+MQA2tKCPCGASjMqOY7MYBG0wBCofQdOPcrUigTqwrao8QyBDoKqP6N4vMgCGtfjqq6+gQ4cOgHPtcds1eeFxqwcffDC88cYbYppYEgBi9A/3D0Zgwzn6eOF+vLiSF8EQF3EmXc8//zwcdNBB0KNHj/KBFDgH8JNPPoEnn3wS9txzT90mLFR6bwCIKuGxb3jUmjxvt3fv3nDWWWeJ49aKcGUCAFHINFYCJwGZcs5dQmvbzPWzOR9YupSxbWF0AUYXUkyGJNOGQF0NonqXqc9hW7r6Ut5rJm1AseuqznFl+bYvy02rnHJ5KUYC8wiAOGe/rq6u3C1atWoF+Ed14a4cm2++uYgC9u3bVyTHY1ZxChiODh522GFiFW4SACKoIbzhriG4u4e8MC+u5r377rtVbsB///tfuPnmmwWHIFDiwlSc/7fZZpsp8xY9gVcALLp4uQVAbBiXC0HQXjWigCkAIFbN1bYwOmCj8xHUhRQTANHxJ/jcm0YCdbRKA1h0NVa9+0zaQGXThWaqMkz7gcpuY/CmTZnQtRv5jwUGwEgZ5Srg8E1cyImje0kXghZCGs7Ze/HFF8tJMWK3evVqGDdunPhNBYDyeFZcrIGRQHkNGzYM5syZA08//XSkG7ifL87x4ytZAQZAix7CABgQrxoAiMXbQqCjYWDx8SWsCtb5gOqlpS8KQV91AUTHl/AjxRAY/ZLRbQPqq8qmrShl+LYvfUitnIIC4M8W3gwt6lpTmjQyzf8WLYOH6s8Qu3boRgBPP/10eOKJJ+Cll14qR9oee+wxOPfcc2HatGnQtm1bLQDE42RxuFheuHcv+oWLSKKupk2bioghTkEbPHgw4KJUvhorwABo0SvSBUB01NFm0GjKZJNm0zl3SWWpIE4FeKr7quFpocUaZS/IehRQAKjGyuBagsCSNm4iSroaqzoWQ2CyQq7aTdUOaW0Nk+YQsCsAxA2UgwCo0vLMM8+EiRMnioUa3bp1KyfHeXy///3vAeFMXhgNxP/H83uDe/TK+6ZDwHiwA643wGlo//vf/+Dwww+HoUOHRu4VrKpPke8zAFq0LgNgQDwVaJkApzSvgjxTMJX2HQGgAI0CRwFtQaqakUBb34OvCdcQaALj1NeWT4DyaVvWL40yRN9IKQpYZADEYV+EP5zThzCH8/+CFx7F9vXXX1f8tvXWW4t1AoccckgFLMpEchEIHthwwQUXiJ9XrlwpFpdQFoHgNi8Igbj7CA5Fd+3aFU488UTAAyJ4DiAAAyD1TRqRruYAEDUwha1qASD6rIJTTONgMQgVAHVhROcjqAsnJhEoHX+Cj40NAOpqFvdYm/oeZU9Xa9WrxqQtVDZd6ZaGntUsIy0ILDIA4glf9913Hzz66KMVe//h8CvuCxh1Rc0B3HLLLcUikUGDBoksCHr4/whxCJVXXXWVAEzKNjDBMmfOnCls4HGw8+bNgwEDBoiVwLV8eQdApHU8dxdX8jRv3rxQWmceAJPgxxTIahQAUUrKMLAPCNSFFl0wMQEPXZ/kg88QmPwKNGkLykvVtL2qbVuW79P/chkpRAGLDIBx5+oidOFRrFQARDvBPHIj6D/96U8VG0HLlcWUPirTLFmyBP72t7/BiBEjYMGCBWJBSi1f3gAQl3tjOFgu0/7oo4+ge/fu8Otf/1qs5uF9AGPm80GHhP6oOQcwLwCIftoM81LycxQw8T1nAh6mH2WGQIZAk4+uaX+jlpXGMHCRAZCqczXS4V6EOCfw4YcfhmbNmomt6HBOoNyjuBo+ZaFMbwCI+/3hhos33ngj4A7heBwLAiCuBMJl5LgSKO9XriOAKH6aW8EklUcBOFtAdACA1YwCisiixkIG3Sgg2q81CNTVNOl9ZaK36v1n0h4qmzp9SGUrfN+nbZdtldiOnqOADIC6vco8Pa4SxtND8A+OQuLRdAh9CH/hk0nMS8l3Tm8A2KVLF/j73/8uCBuPa8HdvxEAcXPIfv36AcJT3i8/AIiqxEUB46KGBptBmwKZCsSyuh0M1tcBBFZrGNjkA2gCJSbQYfrhz0Ik0ETXuPeWid6qd6BJe6hsuqxz0SDQdxQwTQDcd+G90KKuDaU7RKb536Kl8I/6Y0F3FbBxgQ4z4vw+PAVko402guOOO04s/Ig6k9hhkbk05Q0A27RpA++8846AviAAIgjisTDYqfJ+1SQAYqNVYx6gqlzKfQcAiMXkBQJNgMQUOBgC177NTHRPeheatonq/WraZtWyK8v15XfZvscoIAOgqve4uX/ooYeKaB8eOYdDvnxFK+ANAPGMvSOPPFLMA0QAxCFg3BMIj2DBKGDcBo55aigGwIjWMo0AUgBOFX1U3S8AAOpGb0xgxBQ4TD/MtpFAXU3i3jGm/kfZM9GdIZD+9nfZVuFSfUYBGQDpbcwp/SvgDQDxCBec+/eLX/xCjMHjETDvvvsu4AaNOCFzhx128F87zyUwAGoCICY3XX3sAhDRBkMg6algCCTJlJjINQRiYabtkuSoL5jyZdcV8Mf+Q4AjgEKaPA8B2z+9tWHBGwCifHgI9HXXXQdTp06FH374Qcz9Gz58OODmj0W4GABjWtE0CqiK4LmAwJQBUHysCJtD637UdD+uJjBiChu6vslexJFA2lvRtF0YAmn6imfREwRyBJDeBpzSvwJeAdC/+9UtIRcAqIp6uV4JbFqebEoVBNreV/kn/Uj5dBBdADRLr3desE3EKe8QaKJvMmDpa59kjyGwpI5pP1N9ORgAOQKo6iNFuO8VADHqh/P9vvzySxEBDF64ECTvlx0AYu2rvBcgumA6JOtjIQj6Ywt4qvxVAEDxofIQBdT9AKYZBdT1LfwusI0GugIDV3ZKemQfAm3bLe6d7lLHcBk+bDMAMgDmnU8o/nsDwFdffRUGDx4Mc+bMAdzJO3jhTt9F2IGbATCmi6mGWU2hUxangjzVfZV/spwqRAF1P8C6Hz8TCDGNNun6xhBIeWU3TmPaPkml2bZdlG0fNnWfFx2FfUBgmkPA/Rc+Ds3r1tWpckXaVYu+hzfqD8nlNjDGla6xjN4AcLvttoMtttgCLrvsMujYsSOEj4nB8wHzfjEAegDAnEUB0d1qbgtj8gFkCNR/87iEFxP9VR4zBC5QSaR9nwGQAVC70+QsgzcAxJ22cc+/nj175kwSursMgAlamS4EYQAUouoCh356/eFIU8jQ9a0WIoGlNtZvg6S3k2n7cCQwXgHXEMgRQPr3lVP6V8AbAP7kJz+BCy64QGwFU9QrUwCIIjdpES21KYwlDaXaDrNWexgYlaIMBROGgX1EAXUh0ASyTADEFDJM/At2Zts5gbp6pgVFJm3AEBitgG0fi7LKAMhDwEXlF4EMa8IT9Cxqi5s9y2vmzJlwySWXwPnnny+2fWnRohJOttlmG4uSspG1pgFQFalTwVUSAKpsu7hfMAA0ARxT+GAI1I/QJgMlRwJdvdFdQyADIAOgq76ZRTtOAbBp06Zirl8cU8p7vAhEdgWHq4BNI4CYzzQaZ7oSWFWmC8BTRSipACh8rVzEFPcgu54LqAt1Jh8/Ewg0BUDd+kTpzJFA2mfEpo3iSjDpXypvXdt0bs/xfoA8BKzqEXw/TQWcAiCu+KVeXbp0oSbNbLrcRABVsFNEAKRApEoX2fNyAoCmgMUQaP6KcQkcJu2Q5HkeINClflIL1zZdRgEZAM2fNc7pXgGnABh0b8qUKbDbbrtB8+bNK7xetWoV4DFxvA8gypJSBFAFOj4A0KZMKryponyq+yofPQGgADXivoAmUGfyATSBDxvAMPEx+CLhSCDtY2DTRmlFAm37QthP5/YcRgEZAGn9llOlo4A3AGzWrBnMmzcPOnToUFGTb775RvzG+wAyACZ2cQq8qdKo7ksHVPMVMV0Vo4C6EGjyATQBQPTLBjBM/KwFCCy1t7t5gTZtxBDo9li4NAFwi4VToFldW2OSWL1oCXxUvwfvA2isYPYzegNAnA/4xRdfwEYbVYaXPvroI9hxxx0Bh0/zfvEQsKIFVWDleyEINZKo8tMDABYlCmgDgbYAWArgfmP9GnHhhy6kU5zOMgS60kzqkGV7eR0CZgCkPGW1ncY5AB5++OFC0UcffVRsAdOqVauywhj1w5XCvXr1gkmTJuVeeX8AiNJURk7XihU3bJywDUzpSxl/VWMIWECVoguoIniq+y4B0AME6gwD6wKG6QfVFDpMo0ymfgZ7TpYgULedVC9B0/aIsmvaRnE+umi7oG2X9lzaEm3qaBiYI4CqHs/301TAOQCecMIJwv+7774bjjrqKGjdunW5Pi1btoSuXbvCySefDBtuSPl6pymFflk1D4AqwFJF1mwBUFU+5b5sdpWvHgBQfFgyNhewBDBmw4+mgOHiY80QSHt/mbZRrUMgA2AdrYNxqlwp4BwAZe3xCLjzzjsP8ESQol4MgACg4njTyCMV3lTlq+57AEA06WNLGJPokglcpQ2AJvWKeqcwBNLetFmGQJP+mgacMgAyANKernyl8gaA+ZLBzNuaAEAViKkAywYAVWVT7lPTUCKAGYgCmsCSyUeVIdDsnRDOZaJ9PNCYRWaj7DEE6rUvAyADoF6PyUdqBkCLdqoOAKLDMfMA446Cwyw2IFatI+Eo8KYCUIoNjgJGPgUMgRYvh0DWrEIguugSBN3Wc4ET8V35xADIAOikQ2bMCAOgRYMUBgBFZCtBiCwDIAXwKJCogmQpD3E7GExOHQYWUT2NuYBZjwLagoWLj7aL4WATneOjd26ApuSTu0igbVuF6+ui7Vzq7tKWCwhMcxHIxgvfhqZ16xl/4X5YtBi+qN+Wt4ExVjD7GXMDgLfccguMHTtW7C241VZbwY033gi77757pMJ77bUXTJ48udG9Aw88EJ544gnx+5AhQ8RCleC18847w6uvvkpuNQZAyzmAKvCkwJ3LNI6HgXUg0DcAmn4IbWDDJrrkAiQYAsmvskxGAl30AdN+H6UcAyC9P3HKfCjgFADbt28PuM8frvA98cQT4aabboL11jP/F4iU8O9//zsce+yxgBD4f//3f/CnP/0Jxo0bB++99x507ty5kdLffvstrFy5svw7bj697bbbijwIfhIAcZ/CO++8s5wOVyljHahX5gAQHfcxDKyKoPk8E5gCeCr/KDZko1MgsIaigKUPqHnEiSGw8m3iCmqkVZu2Cb/nbNoqaMtlHV3ZcmGHAZD6ZeR0eVHAKQC2bdtW7PPXvXt3wJNA5s+f32gjaBNhMDLXr18/uPXWW8vZe/fuDQMHDoQxY8YoTWK0cOTIkSJ6KFclIwguWLAAJk6cqMwfl4ABsEEZFYDZzD+kwJuqfIoNHQDEtEQIzNIwcAnmzIYiTUHDFipM/Q0+s1mLBNq0Q9S7yLRtomzZttdaMDXrZ2GfXLS/K70ZAI0/lZwxowo4BcABAwaI0z922GEHMbx69NFHV+wDGNTgjjvuIEmCkbw2bdrAQw89BIMGDSrnOeuss+Ctt96KHOoNG956661h1113hdtvv718CwEQ4Q+jfu3atYM999wTrrzyykZH1yU5yQDoAAAFTCm6givAo9ihRAA1ABCTFgECbSDDFipcQABDIOl1KxLZthdDYLzWPAeQ3g85pX8FnAIgwt8NN9wAM2fOhEceeQT222+/ipNAgtWZMGECqXaff/45bLrppvDyyy/DbrvtVs5z1VVXCcj88MMPE+28/vrrgBHE1157DXbaaadyWhxWxohlly5dYNasWXDppZfCqlWrYOrUqbE+r1ixAvCPvBAAO3XqBACXAUDjoehSuoSTOxLvxZ0EorCZtSFgdFcFVVkCQIq/sgNkIAqIrpgAklme6gwFm9Yx/GJgCCS9cgsLgSZ9PqyYbRSQAZDeBzmlfwWcAmDQ3W7dusGbb74JG2yg+vonV1IC4CuvvCKiePLCaN29994LH3zwQaKBX/3qV4B5p0+fnpgOh4cRBh944AGQx9mFM4wePRpwg+vGV8EBECtczZXAqvIp96lpPAAgmixCFLAEYgyBroBUvkdcgMlaW+btE36vZSkS6EIjJzYsj4RjAPQPNVwCXQFvAEh3ITmlzRDw0qVLoWPHjvDb3/4WcMhYdW2++eZw0kknwfDhwyOT1mwEUAVPqqFV1b8BVBFAVfmytVR+UO2o/JXlESOAugAoAMPztjA2EMMQWOoALoCCIVD1Vl5734XetjbyFAFcZ+5saFJnvn/fmkWLYHmnrrwNDL2L5i6lVwDErViuu+46eP/996FJkyaACzfOP//82O1b4tTDIVycV4irgOXVp08fOOywwxIXgdx1111wyimnwGeffaaMROJKYRxqxnmCxx13HKkha2YOoAqcVOClAqqsASDWV+WzZwhMAwBtIIYhkCGQ9JIMJLKFL5v+6gq0GQB1W53TZ1kBbwD417/+FU444QQxnIpbt6xZs0YMxeLcPwSzwYMHk3WR28Dcdttt5cUcf/7zn+Hdd98Vw7YIbAhv4RXBuE8g/o7DusFryZIlgMO5RxxxhIgQzp49G0aMGAGffvqpgFXq1jU0AMSS4+YBJs0PzNAcQAbA+L7KUUDycxxMaDu86AImeE4grels28oVfLkAQBc2bCAwzSHgtCOA+P3Fuf84Lat169Zizv4111wDvXr1Knc0HEU777zz4P7774dly5bBPvvsIwI7m222Gakz4pQuDNLgWoOzzz6blIcTxSvgDQAx2jds2DA455xzKkq//vrrAeENQUvnwk5y7bXXiq1c+vbtKzrAHnvsIUzgxs9du3YVYCkv3I8QO94zzzwDuDo5eGHHwy1kpk2bJraCQQjce++94fLLL29Y1EHzjAEwoJPvKKDKvgpSpasUO5iWGgHEtJ4gkKOAyc9hliDQBVi4BKW1ttzMCSwSBNr2GwbA6Ody//33h2OOOQb69+8vFlRefPHFYu497tcrt1879dRT4fHHHxffalwfcO655wLu24uLL3HruKQLd+3AwM1XX30lRhIZAGmckpTKGwC2atVKROh69uxZUf6MGTMEwC1fvtze+ypbKBQACpBJEFQFTqr7KqBSDQOr7FMBUCedymcplycAFFCRwlxAG3ip5lCwjd/Bnu4qEujKH5d2SraKBYHWAGe4F2YZqi0WghQ5Ahj+eiCodejQQWzVhsGahQsXin2BcfEmbhGHFy7yxJ00nnzySbFrSNyF07hwKtjTTz8NBx10kIA/BkB7APIGgAh+SOkYsg1eeIoHzgv8+OOP7b2vsgUGwEADqABNBVMqAKSAm8oH6S41ncpnBkChgA1guIgs2QIB1oEhkPYyddFeLgDXts1t8nMEkNZXMNiDCysxCohBn+eee04M+WLEb/311y8bwVO6cEQueocNgB9++AH23XdfMecfF3PiaB8DIK0NVKm8ASCe2oGNhEfC4VwAXATy0ksvidDv/2/vPMCuKK4+ftSgAkEESyIqKCgW1GAXe48tCmoM6icoKlZiSaJRRETFhj0qikZFEAuxoGLJhwWwJ3assYJi7IJg+0C/5z8vc9l33907Z8ru3b33zPPwJL47MzvzP7O7v3vOFBwRFwdDU0OLeF0AsGAAyIFEbp4mKuAn8QLytYrlDAEVPh903RyBQJ4Ji2AvX3v7lG80AJwxYwYtFVlNjOge/lVLmPMPYPvqq69o6tSpKuu4cePUuoDoXrr4+84770zYNg7OoaSEuYWPPvqo8v6BIwQAec8pJ1dmAIibY8HHRRddVJnvp1cBY2DUQxIAFACsKFAgAPTxsLh+GH28gGhvEaCiifm/CPZqctUy3oBQ9TSNC/9wcAhb+YxRrY+PLl5lSxICphe/Jmrnvg0MfTObqOfSLZ6HoUOHqrl41dIxxxxDEydOVE4fvcAjDQAxR79bt26ERZ7xhLmBCPk+//zz1KlTJ3VZADDYK4oyBcBwzSxmTQKAJQVANLtEYWD1sSz4XMAQcBECLHw+7Ho0CwSa37dFsZWPvb3KOkJgnnMAQwGgrQdw0KBB6pjVKVOmKM+eTi4h4EsvvZROPPFEWnTRRSv1zJ8/X/035g5iBw9J7goIALprRwKAAoDNhk+DewGLAIE+H/WoLQUCzS/GIkCgj729yjYQAGLxRjQEnDYyEPYF/CHy99hjj6n5f9GkF4Fgi7j99ttPXcKuHvAQpi0Cwf68yBNNWCxy0EEHqXBydIsZ84iVHHEFBAA9xoQAoABgXgCYpxewCeS+dn4yfMKMRYCKLEDQR8+4IULU5WOjaHuKYC9XPVzLqedDALDF++Hoo49W8/wmTJjQDMzat2+v9gVEwjYw9913n1oL0LFjR7UnICAvug3Mmmuuqfb07dOnT+I7SELAzq/mFgUFAD20FAAsIACiSZzwLicP6spoIQiqtjkf2AUAfUDO6+PoOc+sCFAhEMh/Mfray2es+Yxxn7ICgC3HBxZoJKUbbriBDj74YHUJ279hdxCAYnQjaIRzdUI90TLxOgUA+c+mKacAoEmhKtcFAAMCIKoKsRUMFwBt8mUEgbYA6AKBPh9Xv7J+iw18ocLn4570yIcKCftoGm9XiLpCeAKLYCtXLZzLiQfQ48spRYuiQGYAeOaZZyr3bps2bZr1FdQ/YsQIOv3004uigXM7BAAFAFsMngznAboAoC8IuX4gm+4rEJj0cvHRVCAw+XXtqqlzOQFA5++mFCyOApkBII51weRN7Leqo7UAACAASURBVAQeTYj3429YyVP2JABoAYDIavKkNZgHEJLk4QX0gUDXD6QeGQKBYYElK6D0tRPa5esJ9B9rbvNWXe/rEgbOdRXww7OI2npsAzN3NtEO7dUJHpxFIGX/njdi+zMDQCzT/uSTT9TRL9GEpeA4BgbHxJQ9CQA2CABy4DU6mMUL2OzR9oULX7BYCKNugBB/T0k4OP3N7WsrVxirxY8cAcCyf8Gl/cEBEEe8YBKn/tUQnRgKr9+cOXPoyCOPpCuvvLL06jcUAMJa1RZOcBZVlNUDmCEAiheQ9xrwBQuBQJ7OvrDu6wn0AUBXCHS9pwAgb0xJruIqEBwAR48eTdgPCEfAYRNHLAHXafHFF1e7ePfq1au4ili0TAAwJpYJAvMCQBOsRpttarPOa2q7owcwTwB0/UD6lFsIXn7zAX3BImoe1w9+Vp7AEPqGBFxfCPSFdR/7uJZ1KScAaPGxlKyFVCA4AOpeTp48WZ0B3KpVq0J2PESjBABrAIBcuOOCHTefDQCijQUMA/uChstHsjl4CQRWe+/46isQ6L5/pYv2AoAhvqJSRy0VyAwA0amffvqJ3n77bfr000/V/4+mrbfeupb9DnLv0gEgel0NZEyLMEywZLpugijT/bXVTPfhQqJNPpN28RFlAYAomtdiEB8IdPlIxmXx9S6hPl8PU0hQahoW9Xl+sK+tfOzkM9ZcyjqVcVgJLItAgnx6pZJACmQGgE8//TQdcMAB9MEHH6iQcDRhXqCsAoYiv6pixuarp5tnrFJuEYPHtcgAiE5yIFAAUA0H2/OBQ4CPy4cyOnZ9oULX5QMXzdtTrIUhPoDeErb9++ZrLx87uY61PMvZegEFAAORi1QTRIHMALBnz57UvXt3GjZsGK2wwgpqYUg0RecGBulJDSopnQfQ1wNnAi/TddP9awGAuKep3XpscdofHYd16AUMASi+UCEQyH/ZucJQSGh3hUCftruUdSpj6QUUAOSPXcmZvQKZAWDbtm3ppZdeotVWWy37XtToDgKAMeFNIMUBqLw9gA0KgL4g5/KxDAkVAoH8l56vrZrGivv8TVcA9BmjLn12KlNkABw/i6iNxz6A384m+r3sA8h/0sqXMzMA3H777emkk06iXXbZpXyqMFvMA8BqYd6cQ8AmADPBlwnwTNdN9xcPIHPkNc9WxlCwL1SEBkAf2IgbrR7nBPoAIPTJGwKdYI7sQ+aFDgELADq9TxupUFAAfPnllyvavfPOO3Taaaepg5/XXXfdFquB11tvvdLrLADYYB5AdJcDsVFZChwG9oUel49s/KH3BQtfuGjZHnsISHuRhQLBEDr72joEsLtCoGv/XcrZlhEALP1nvKE7EBQAcfoH5vrFF31ohfW1xloEIh7AyhPGgSeTFxKVmTyN+oah8xUUANXHeSk3cLH94IUGphAAKBDI/4b529s9FOxjJ5d251FGAJA/9iRn8RQICoBY8ctNXbp04WYtbL6G8gByYMqURwCQNZbz3BImjGfIDT61GAKBrGFBLkCTVLNvPb72cvEEurbZtpx1fpkDyBu8kquQCgQFwEL2MMNGCQDGxBUAbDnaLEPAqMAFAH28gL4QaPvRTIYSP8+SrtMFLtJeESH6pesuUjg4RL98INDVRi7tzrqMeAAz/MBK1ZkrkBkA3nPPPYmNR/h3ySWXVKuDV1111cw7mOUNBAAj6prgD1mL6gFE2zjt5/YhPugsIbBsAOgLkOIJtHtLuUBN/A5h6nCHdhcIdGlzHmVsIDDXbWBkEYjdg9WAuTMDwLT5gNF5gFtuuSXdfffd1KFDh1JKLwAoAMgauJYAWEYvoEAgayQEOzXEBWxCQ6CPFxBtKSoE2mpbWAC8ehZRa49tYL6bTXSkbAPDe7LLmSszAHz44Ydp8ODBNHz4cNpkk02UOs8++6xaGTxkyBDCRtBHHHEEbbrppvT3v/+9lOoJADYgALp4ARsEAAUCea+xegoH+0CgAKAHnFUZapXvkgAg74Fs4FyZAeA666xDo0aNos0337yZvE888QQNHDiQXn31VZo0aRINGDCApk+fXkoTCAAKALIGbkkAMATA2XpP0vTzAYtonS6Qkd4mv8Uu0XoFApvUcLGPyxizLWOTXzyArLegZCqgApkBYOvWrelf//oXAQSj6ZVXXlEewe+++06dE7zWWmvRt99+W0BpzE0SABQANI+SBTkEAtlS6YwCgXzJbIAlqVb/8sWeD2jbP5v8AoD8cSo5i6VAZgCI+X3t2rWjm266iZZbrmlzt88++4z69etHc+fOpSlTpigP4NFHH01vvfVWsVRhtqauANC0/55pkYTpOjSth0Ug3H7Ex1ADAWAIT6JAIPMlFMlmAy1FgsA8vIC22tjkFwC0H6tSohgKZAaAb775Ju2111703nvv0corr6w2iEaot2vXrjRhwgTq3r27WgDyzTff0EEHHVQMNSxbIQDYoB5AFwh0AEDcphYrgkMBnM1HtNqjJ55A/ovJV3Of8j52KhoE2uggAMgfn5KzWApkBoDoJk4Eeeihh5SHD/9/zTXXpJ122omwQrgekgCgAKDVOHaAQAFAIh+wiNvHBTSqw2n9zQu0gZ+4Nq62crGLbTuzzM+FwFy3gZFFIFav50bMnCkA1rugAoA1AkDclhNy5uSJDlKb/JxwdvwByBEAlRfP8Xg43WzbD2YWYcWFbXGfYyYQaP8m9rF9USHQtk82+QsJgMNnES3psdL4+9lEg2UbGPunpzwlggLg5Zdfrlb4YqNn/P9q6Y9//GN5VEppaW0AsMrZwou0qq5pNWgp0xzABgJAdLVWXkAFkeTv4QpRR1NbBAJtXpq+uvuUd7WVrSfQto02+a3yMo+Ey9UDKABo87g0ZN6gAIiTPf7973/TMsssU/WUD8wHfPfdd0sveHYAuHwVbQQAlThcbx03n02dyJuTB9AHAIviBQwFkgKB9q9MG4gJ7cF1gUBbALQdW7Z6cPOLB9B+bEqJ2isQFABr3518W1AqADQBSz16AG2hzgYWSwKBvmFg2w9s2hPI/ZBynmAXsKhWrwt05NHPouwV6Go7VzvZ2sO2fTb5uXkFADlPruQpmgKZA+CPP/6oVgJ369aNfvGLX2Ta/6uuuopGjBhBH3/8MfXo0YMuvfRS2mqrrRLveeONN9IhhxzS4hr2J0QIm5MEACMqceDJBKGozgSi+pac+wkAKrV8IZD7ETQ9M6HqCe0JRH220FGtryH72YgQ6GILG82zyCsAaHr65XoRFcgMALG586BBg2j06NGq31gJjC1gMPevU6dO9Ne//jWoHrfddpvaTgYQuMUWW9A111xD1113Hb322mvUuXPnFvcCAB533HGE7Wqi6de//jW7XQKAgQEwNPwJAAYBwCbg8p8LGLIegUD2a6qS0deGruVdPIFZQqBNP7h5BQDtx6OUqL0CmQEg4ArHvsELt8suu9DLL7+sAPCee+6hoUOH0gsvvBC09zhTeIMNNqCRI0dW6sUpI71796Zzzz03EQCPP/54+vpr9w9bwwAgx9tmylN0758tLCI/p0/xkeewElhX4boYRMGS54rgkPDG/ahyXhAucFGtXhfwSKsvZD/FE2geDTZ6c/Oy88kiELOBJEfhFMgMALt06ULwym222WbqRJCXXnpJAeDbb7+tQA3wFCohzNymTRsaP3489enTp1ItIPTFF1+kyZMnJwLgYYcdRiuuuCLNnz+fevbsSWeddRatv/767GYJAC6QygR/XFgK7QHktCtqbdv8JQLAEBDI/RhyHqCwdYVbHYy2h4TAkOBcBAj0sZstrNvawaZtofMW0gN44iyiJTy2gflhNtHFsg0M531W1jyZASCAbNq0aQr6ogAIENx6661p1qxZwTSbOXOmAjl4HDfffPNKveecc44KQcfDvMjw9NNPKxhdd911FYxedtlldP/99ytQXX311RPb9sMPPxD+6YRyOOWEaBgRtQwzN+VLW7VbZTUvZbAKuNG2gIH0tkBnm98FANEuRy+gjwcwBACGhBmbDzDnRWELF6Y6beHDVF/I/oYAQZ/2+JW1g3VbO3Dbxs1nM+Y5EJjrNjACgKbHsuGvZwaA22yzDe27775qHiAAECFgbBNz7LHHKvB68MEHg4mvAfDJJ5+kXr16VeodPnw4jRkzht544w3jvX766SflmQScpu1heMYZZ9CwYYC9eCo4AJpAxeR5M4GR6TrkMrVBgZHRTE0ZOPezyadvy61X5+f0KalLjgCIqgQC08eIQCDz+VmQzQaC4jX7lc0OAm3axc3LzscIAwsA2o1RyZ2tApkBIGAMc/8OPPBAwoKLI444gl599VV66qmnVEh2ww03DNYzlxBw0s0PP/xw+vDDD+mBBx5IbFuhPYDVNoE2gYoJvExgZLouANh8PAkAVvTgfly5LwuBQK5STfl89Hcta2ujWnsBuf1sdA/glClT1C4czz33nNqJ46677lJz8KPp9ddfp5NPPlkxAJwu2K3j9ttvT1yoqcthHQHm9k+fPp2WXXZZ5VjCvH7ubh12T0Rj5c4MACHjK6+8QhdeeKEaENrDBuMj7Bo6YREIoBKrgHVae+21aa+99kpcBBK/P84q3mSTTVTbrr/+elbz/OYAuoSAHTeBFgBk2ZPtWdS1mXRNu2sNAVB99Au0IMQXQpIktgUM0+CwBRBTfSH7XNZwsK2NbGzABjaLle2cOhsdAOE4wTQsRNL22WefFgD4zjvvqG/soYceSvvvvz+1b9+eAIQbb7wxLb988rSnm2++WeXHNxnTu7CbyMEHH0x/+MMf6JJLLuE8apKnigKZAmCeyuttYK6++moVBh41ahRde+21yuuIBSn9+vVT8wT1imCEcrFABfP9AHII+yJcjAGMQcpJ2QBgzvP/0NFqHkCTd890HfVzQMnkhYwahHNP5Ofm03Xb5uf2LT6YPAAQVRUhDBwSYkLX1VSfXZiR87zbQAinPg5UcOppBAi01Z6rbch8jQ6A0bGK077iHsC+fftSq1at1HeWmzBlDJD48MMPV4r86U9/omeffZamTp3KrUbypSgQHAC5q3uXWspjdVJKZ+D9u+CCC5T7eZ111lG/EDCnD2nbbbelVVZZRYWjkU444QS688476b///a/6JYLVv5jjF51DaBo1dQGAJvAyQZHpOheSTO2oJwBU0P2zaXilXvcFQAVIAbyAocGN+zHmCtdIENj0mH3BlSYxn4/+rmVtbGQDgdz2hMxXrwA4Y8YMin6vl1hiCcK/aikOgIgA4jt70kkn0eOPP662gcOagFNOOaVFmDha76233kpHHnkk/fOf/1SOGRwhu/vuu1P//v2D7yXs9fCUtHBwAFx00UUJxk9LCLXiOrZeKXsSAGR62UJ6ADnAqQeWTV6Usc3PhdukgS4AGBxCkiq0AQzu+8gGRDh1ciGEU1fZINDWPjbac3Xl5GPlqdNFIPFxh3184SyxAUA4WlZYYQW1XdvZZ59N2223nVoIeuqpp9Kjjz5KWDSalv72t78RvH5gh3nz5tFRRx3VbKoX57mQPMkKBAfA6J57MNhuu+2mTuRA+DWaqhm8LMYSAGRCkwBgyyHtAYCorF69gOgb52Nr846whQxO3TYgwqkvZJ/rGQJtdOdqysnHylM0AOw3i2hxj0jbj7OJbmpPITyAeqcOzP0bN25c5ZHYc889qW3btnTLLbckPiaPPfYYIXQMaMQ8f+wggv19sWBzyJAhnEdL8lRRIDgAxu8V3QOw3ixRCgA0wZcp9Gryipmuc71kpnbowcO5n0telLGpW9/DpG/aoK8jAMwG2txP6EmTvJFAsFYQyAGlJPvY2IYLgdy2hMpXuBBwIADEnr22U7biIWDs1AHQg/fwtNNOqwwBLApFSBhz75PSVlttpebqY3WxTmPHjqWBAwfSnDlzCBFHSe4KCAC6a6cWj2BeQ9iNoB0WgdRqCxguMHEgSQDQaiSG8AAqcAs0F1Ag0Mp8VTNzgcR0xzJBoA0Aot+1gECOXUwQmOs+gAUCQNgMq3i7devWbBEITu5q3bp1M69gdFxjZ48dd9yRzj///Mqf4S0cMGCAAsDFFlvM9BjI9SoKCAB6DI+6B0AT4Jmuc71/yNdoAKj67L4QBMVDQGCRATALqGyqU1YI27z2OOATr8+tDN8uAoDpFqx8l3IGQAAZQrRIWFR58cUXq7l+HTt2VPv8YVUwtm+58sorK3MAjz/+eEKYd8stt1Tl4rt1YK4h6sGuHjoEjDmAAEPs/CHJT4FcAFCfAuLX1OKVLj0AmqDLBHim61wANLUjanrOPXV+m7woY5uf27+0oVsAAFRAVGAvoECg/3uvXj2BISGQA6ysPIZ5gPXsAQTIAfjiCSt29e4b2M8PW7HhwIU11lhDnayFvXp1iu/WgUUf+kSvjz76iJZbbjn63e9+p/629NJL+z8cDV5DcADce++9m0l677330vbbb6/i/9GELVjKnvIFwAw2gTaBlwmITNe5gGRqhwBg4qMSwgMYGgCzAzaZE+jzviwLBNp4Z0MCIHfcmiBQQsA+o1TK5q1AcAA85JBDWH244YYbWPmKnEkAkGGdkPP/bL10HEB1hUtdjtO/jDyAqFYgkDEGq2SxAQ6bO3HhhFunCTw49QgEVleJo7EpjwAgZyRKnqIoEBwAi9KxPNpReAA0wYnJ82YCKNP1RvAAcvuYNCA9Q8CNBoBcL43tsy8QyFfMBEBpNdmWs7EJB7Y59w+Rp1AAuMcsolYe28D832yi+9qTyypg/oiSnLVUQADQQ30zAKaFbdP+7rACGO1PWwXsA4AmuDNdtwEjE4hqG3Hu6ZLXp4xNPwUCPZ62hUU5H2qXG9lAB7d+Dpxw6woFwGXwBNrYgqsxZ9yY8hivF2kOoACgzaPVkHkFAD3MvhAAr0IwLqGmHADQdQsYE3SZYMt0nQtGpnZEVeXc0wfmbOqPtssE2tXGWIG8gAowAi4ICQUsSfKZPsSuj7UNeHDvwQUUbn2h+u4Dgq5tsClnYwuOxpx7m/IYrwsAcoex5CuAAgKAHkYQADSIxwGjsgMgp48NDIACgQuNz4EUm9eRCUY4dQkEtlTJpGu16xIC5ow6yVMUBQQAPSyRHwA6rAA2gUk18DJ5wkzXud4/5OMCIOeert5ClLOtX9/LpLNpfAXwAOIWoRaDZOEFLCMENrWZvy+dycz6ukBgkxImyIrqybUDR1vOfU15jNereAFz3QZGQsDcx7Jh8wkAepi+tABogi4TDJmucwHQ1A5XoOO0L253lzLcflYbYw0CgLYffZvH0vRBtqkrnpcLH7b34MCKTZ0hNHD1Brrem1vOxgYmXTn3NOUxXhcAtBm6kreGCggAeohfcwCsxfw/DihxvWICgN6ngejhG9IDWDYvYJZwmZUnEPWaYMX21WQCE059eUKgTXu5EMjRlHPfqmFeqr4nZbUwsHgAOaNQ8uSlgACgh9J1CYAmwDNdt/GKcQGQc09Xb6EuZ3sPXY4Lu2njLJAHENULBIbfLDpqNi6E2L5SONDCrZMDN6a66h0CORrVBQD2mkX0C49tYObNJnpKtoExPS9lvi4A6GG9sACYtgVM4Pl/JugygZDpOhcATe1wBTpO+5Js7lpOAND6CeJ8gK0rXVAgy7obyRtYZgjkADVnnLhCYGE8gAKArq+RhiknAOhh6oYDQA4kcYGIC4Cce7rCoq/3jwu7pjEWyAsY2gOogCfwtjBaCs4H2CRbtevZ1x9+gQj6w4EXri4hNMgLBLlt5XphTTpy7icAKB5A7rNW1nwCgB6WcwNA202gU/K7zP8zQZcJtkzXbYDI1BZXOOO0MW5zlzK6Di7wVhtngQAQtwgNgVkBYJM3LeuQbdb1CwSGhHDueOBAoAkAOeNPAFAA0AMPSlFUANDDTAKACeJxgIgLf6jeFs5s87vcI9ptTn9NY6zAAJilF5DzETZJZ7rOhQpTPdVBp/5BMA9PoI2tQkAg535pearCYVFWAUsI2OexboiyAoAeZi4kAFYDEhN4VYMnDlhxYcjUjjy9fwKArCegzJ7AfEBTIDBtIHFAS5e1y1td81p5AWUOIOuVIpkKoIAAoIcRsgfAgAtATNBlAjzTdegoAOg2mgJ6ANGA0GHgrL2AeQBafvcID4IckOEOPBvASqqzSJ7APLyAocPAuW4DIx5A7mPRsPkEAD1MHw4ALVcAF3H+X2j4c/HMcSA1yd6u5Wygt9o4KwEACgTyXxQcMOHXtjBnUUCwTBBo0swExKUGwO6ziBbz2AZm/myit2QbGJdntSxlBAA9LFU4AGzk8K8LMLqWiY4ZLvgKALKeNNMHmVUJI1Me98kCBE1Aw+h6JYuPBgKB1RcxpYWBc/UACgDaPA4NmVcA0MPspQHAooR/Te2I2sLWK2ebX9/LtZwuX0AARNOyCAPn4QVU98h4dbA2XR73yQIC0f5QIOijQdYQaNM2k87V9DLdJ+RCEAFAjw+uFA2ugACgh6TVATBt/l7S35NCwCWa/8eFIAHA5NEWOAScJQDWGwTmBZwmQHF5DYWCQF8NXEDQBF0ukF5NY5NWzqHelB8r4gF0GdFSJm8FBAA9FK8JAIae/2fygJmuQ7/QAMi5Z9xuLmVQh2s5m35zxlhgCMzKA1j5MGe0QXRcKi4ocCSulieP+2QBgWX2BnI15+dzXxVsuoetF1AA0PeJlPJ5KCAA6KGyPQDabAJtuQF0GoSZvG4mADJdDw1/rlBmameanV3LFRwAs/YC5uUJ9PVO2TzeJgiwqas6bBZ3pbCPBkXwBObtBUwFw5S9ACUEHOopknpCKCAA6KFiGAAMtAK4XgDQBchcyriCZnS8cOGXM8YCewDrCQDzhMC87lVkb2CeEMi9Fz9fNl5A8QB6rCbmvP8kT00UEAD0kD07ALT0/lXzRlXzAJrAyXSdC0AmL2TUBqZ7xu1lm9/nXvF7c/vPGWMZAKBAIEf45Dxc4HC/Q1PJLEDQNN+N22ZXDWw9gdz78POlQ6BJG2vQkzmA3OEk+QqogACgh1EKA4Au3j8TOJmu28APFwBN90yylUuZEN6/atDtMqZKCoAKYnKaD5iXdy5qPi50uJi8+X3qKyxcSwg0QbXrimAbOCzEHMCOs4gW9fDc/TSb6EvZB9D32S5yeQFAD+tkA4A5rf41gZPpOhcAiwh/IQCQ23/u+MoIAPPwAuYNgfUKgiZw4Q6leD6T14tTrysI24Ag9x6cfFnMBbQBQDVGE+YB5joHUACQM7QbOo8AoIf5cwXAtNW/Rfb+ceHPBchMgFrNrj5ldb0CgC0UztMTWK8Q2NSvxvUGcuBODzxTXlcItAY9izCwAKDHB1eKBldAANBDUn8AtNj/LxQAmuDHdJ0LPlnCnwswRu1s6iNnTHB14NSFPBl6APPyAtbCEyggyB1gTfkayRtoAmmXUHASHNqsBBYAtBuvkjtbBQQAPfS1A0CPDaBDwR8HnExwxAUfLgCa7pdkH5cyuh6fsll5/+oIACvemRznBdYCAvO8pwlkXF5hviBo8r6ltSnvkLCrFzDNvgKALqNNyhRVAQFAD8uEB0DL1b+24V8T/Jiuh4Y/DpDG7WNqYzV7+pSN1svVwWZsZewBzNML2EgQKCD4tc0or+TlgiAXNKvlc4VALuyJB9BpCEihAihQGgC86qqraMSIEfTxxx9Tjx496NJLL6WtttoqUcJrr72WbrrpJpo2bZq6vuGGG9I555xDm2yySSX/wQcfTKNHj25WftNNN6Wnn36abZaaAmDe8AdVOODD9fy5wJ9rGW1RAUD22A6VMe95gXkCWVwjLqz4ahvaI1h0byBX1/S5e+G2hbGZHyiLQHxHupTPWoFSAOBtt91GBx10EAECt9hiC7rmmmvouuuuo9dee406d+7cQqMDDzxQ5dt8881pySWXpAsuuIDuvPNOevXVV2nFFVdU+QGAn3zyCd1www2V8osvvjh17NiRrXkuABgq/GuCH9N1gb+mccHRgT2CIhlz8ADiblkfEZfU9VpAoICg/SCsBQiG9AQWzQtYcwCkWUSLeGwD8/NsIpJtYOyfpPKUKAUAwjO3wQYb0MiRIyvKrrXWWtS7d28699xzjWrPnz+fOnToQFdccQX169evAoBff/013X333cbyaRnSATAplBv/G3MBiA0ApnnfTHBnus6Fniy9f6Y2mqzoW17Xz9XC1J749ToGQAVjOc8JjMrL9SDZmsyUP4/7hvYGok8+IOjS5yJDoE8YWADQ9ITI9VorUHgA/PHHH6lNmzY0fvx46tOnT0Wv4447jl588UWaPHmyUcNvvvmGll9+eVXHHnvsUQFAwB+8fksvvTRts802NHz4cJWPm8ICYAI02sAfGu0CgBww4kBPlvCHvnHamWY4n7LxOjlacAdQNF9OAIhb1sIL2KgQqPqdsk2IyzBJKxMaBH0g0KXPeUBgqLmAbCis9T6AOXsAp0yZoqZpPffcc2qq1l133aWcNDqdccYZdOutt9KMGTPUdxdTs/DNhYMnLcHBg+jdG2+8Qa1bt1ZRvfPPP5/WWGONkI9Pw9ZVeACcOXOmCts+8cQTyvg6YU4f5vC9+eabRuMdc8wx9NBDD6k5gQgJIyGs/Mtf/pK6dOlC7733Hg0ZMoTmzZunBu8SSyyRWOcPP/xA+KcTAHDllVcmoqvwWY2UcfEAZgSAJvgxXecAT5HhzxceoyOBo4VxNKZkaAAArDUEuoCJqzmTyjUaCNr2NxQEhg4Fc4AvMU+DAeADDzygvtOI1u2zzz4tAHDcuHHKwdK1a1f67rvv6JJLLlFOmbfffpuWWy75I7LLLrtQ3759aeONN1bf58GDB9Mrr7yipn+1bds25OPZkHWVBgCffPJJ6tWrV8VI+OUwZswY9cugWsL8v/POO48ee+wxWm+99VKz4hcLYBC/UPbee+/EfPgFM2zYsIRrJgDkbAETy1ME7x8HeGzgzwXGTIBqemx9iOmZMgAAIABJREFUy9chAKJLtfICFgECBQRND03L6z4eQRsQrCUEpvWRA4BJY6qRQ8CLLLJICwCMjyodQZs0aRLtsMMOrEH52WefKYhE5G/rrbdmlZFM6QoUHgB9QsAXXnghnX322YQBttFGGxnHweqrr06HHXYYnXzyyYl53T2ADvP/QgBgNfjhgFFoAOTcM6q8bf641XzLx+vj6GEcZSkZcvQA1hoAiwKBAoJ2g7VeINAlFByHQDYUxryAuW4EnXMIODqaTACI7/rll1+uvs/wAC67LO9ljbz4TsMLuM4669gNYMndQoHCAyBajDkCmC+AVcA6rb322rTXXnulLgLBXAQMLoR+N9tsM6Ppv/jiCxVqHjVqVGWhiKkQfw6gCQCZ4d+QW7+YnjcO7Nh4/0z3SxLbpUxIgIzWxdHDNGCqXc8ZAAUCmxvDxkvlY+aksnncO+QcwTxAMIQnMGQouJEBEHP2llpq4WpiTJFKmyalx3caAN53330qpPvtt9/SCiusoBZhIrzLST///LP65n/11Vc0depUThHJY1CgFACot4G5+uqrVRgYkIa9/rCtC8K2WNkLeNMrghH2xZw+zDnAdjA6Yc4f/s2ZM4cQzsU8BQzC999/n0499VSaPn06vf7669SuXTvWwMkMAMvg/bOBP6hpC3O2+eMW8y0fr08AkPVM2Gaq5ergWsFYmkZlAsF6hsCkvnE8fpx5gPl6AD8nIo9tYAjbwLR8kQ4dOlR9P6ulNACcO3euWiDy+eefq2/4I488Qs888wxr8SXm8k+cOJEef/xxWmmllWxfNZI/QYFSACDaDe8fwA6DB65fTCDVcwC23XZbWmWVVejGG29UXcT//+CDD1p0Vw9cTEDF6qQXXniBsBUMIHC77bajs846a8GiDt5YyRUAQ3n/OGBkgp2iw58LcFYzuUkP3nCpnqsGHkA0qJZzAbUgRYNAtCsPGGtkEOTqm6Un0DYUbPIC1isAhvQAxsc8wrkDBgygU045per7cdCgQcpbiJXGq666aog3rtRBRKUBwCJaKwwAZrD4w2fuHwd2bACQA5yhvXcu9xQArPkjJiDY3ARcSHI1XKiwsKs3kNs/Xwh0CQVzvIAsr2BN5wCG8QDOmjWrWQiYM95McwB1Hautthr9z//8T6pHEWFfwB+2lMFCTgCjpHAKCAB6aOkGgPF9Bj0A0HbfPw4YmQDQBv5cPHGcNppsFqKO6D1Mmpjaw7leIw9gUbyARfYG1rtHsJYgyIVA2IADgmn1hYBAFvDF9nyMrwQuYwiYC4CYWoVFGkjrr78+XXzxxSqyhtO1lllmGbXn35577qkibphzj6je2LFj1dZrON4VCauBsd/vscceq/776KOPVlO5JkyY0Gzvv/bt26t9ASX5KSAA6KEfDwAtF4Akzf9LApCiw58LhLmUidrPt3zSWKhz+CsaAAoIpr+QbGDJ5bUWAgSL4A20hUCbULBtGLiRABAeOgBfPPXv358wf/+AAw5Q8/0w/w9AiMUfp512WrNFIJi+hWNa9RxDeBKTEo5wRT5JfgoIAHrolwyAceCzAEAu/KHNSQCYBkAmMOJAjo3nz3S/uOa2+ZNsFqKOeL0cXTzGT6VoDb1/ug1FmAuYJGURw8JoZ9YwZhpWWd/fFwZdQJDbJ5MnsLq37+tEadP6G++HAKBpZMr1MikgAOhhLX8ANIR/bRZ+CPx5WDKhaF7wp2D+57Btd6ytqBCogKuGZwmb5OSCi6ke1+tZ3r8eQTDdQ/hVCxMIAC5L3BCw6/iVcrVTQADQQ3t7ALSY/1cE+MvS6wfdQ3jtQtRRK89fgeCv6J5A3b4ig2C9ewbzhEEbqK3mEQwRDraBwBYewpouAnmNiHhbmiV/Br8horUFAD0YoehFBQA9LBQUAKPhXxv4SwMpExiZPFxZwp+pbRybhKgj6T4mXThts8lTEO9ftMlF9gSWBQRrDYM2AGUzXJv61dJTZlOHTXiY248sIdAGAON2r+0cQAFAm3HZiHkFAD2s7geAKeHfrOGPAzhc+HOBMJcycRuFqKMI8FdAD2BZPIECgvwXFxei+DU25SwaCOYFgT7zAPNdBSwAaDumGy2/AKCHxc0AWG0BSOSayfsXasWvwF+6tTnaeIyVxKIF9P6VzRNYJhCstVcwq/uXBQR9wsHVvIBGIIyEgQUAQ79EpT4fBQQAPdSzA8Aq8/+qAWBe8Mf1+kEvWw+cbf4km4SooyheP90OAUCPpy+5aNHnCEZbnZVnjitqFvd3hcHQYeE0b6ArBNqEgavNAxQA5I5OyZeHAgKAHiq7AyAz/Fs0+HOBMJcycZuEqKNo8If2FBwA0cQyzAdMe4QFBnkvt3oFwRAQGAXaKARW8/oJAPLGneSqvQICgB42CAKAad6/ssNfKGgLVU/czrUI+UbbUAL4080tMwSWLURcaW/sRAmP15RV0dAwWGuPYEgIdA0DRxeCiAfQajhK5owVEAD0EDgoAMaBJMRGz9Ugp1FDvrB3reGvJN6/6KNRDxAoMGj3sisCDHJDw9Xa6guBHC9gVY+gzAG0G3iSOzcFBAA9pM4MAIsCfy7eN5cyeYV8iwJ/JQTAsoeD6yFMjD6EhjLu6y/kfV28gr4gaAuBca11m13CwLXzAE4mol9yTZyQbw4RbSP7AHooWPSiAoAeFmoJgNG5fdH/H10AYlj9W1b4CwF+sEWoepLsWgTPX0nhT8tZT57A+BAp05zBWsFgLUGwlhBo6wWM6iQA6PGRlaKZKiAA6CFvMACMgkkcAG2PeAsR9rWFMNv8aZqHqqfI8FdyAGwEEFRwVeCj51K9mTnPGwwFg7YeQQ4IprXN1hOo6xEAXMrjSylFi6qAAKCHZbwAMGnxB9f7lwZKtptIx/tuC2C2+WsBfrhnUTx/uv8lWgBiejzq2SOo+y4waBoFYULToUEwBAQ28+QtOAElLQzcPO/XFdHEA2geP5KjNgoIAHrobg+AluHfUN4/zoIPW5izzZ+kc4g6TPYT+DMpFOR6I4CgAKF5qITwCoYEQVsIVN7fmCfV1wsoAGgeN5KjNgoIAHrozgNAw/y/tPCvwJ+HZRYULRr8oVl15P2LG6iRIDDa97J5CENAGufh9L1PKBD0hcCQXsB8t4GRRSCccdrIeQQAPazvDICm8G9e8GfrgbPNL16/lgrUMQDqzjYqCJbVQ+gLapxXqOs9igaBSXMBU0O/EU+i9gIKAHJGi+TJSwEBQA+l0wEwaQVwQvg3yfsn8OdhkQLO94v2pgHgL9rdRgfBMgKhK6jZPLQu97ABwbRFIj6eQK4X0DQPMF8A/AcRtbExTSzvt0S0r2wD46Fg0YsKAHpYKBgARufoJQFgyEUftl482/zi9UsfUQ0GgOIRTB4KZQoZu8Aa95XqUjcXBENB4ML5f5FFHQsWg6Cf+j5J4Je0FYwAIHd0SL48FBAA9FDZDIAJ8/+qhX+jsCXQZ7ZMEef4JbW6QcEvLoV4BNOHdBmg0AXYzA+x+wpiEwxW2y6G6w1MBruvKt0SAORYWPIUVQEBQA/LeAGghpe496/M4BfCW2iyR1mgT/dD4C/RogKDpoFejr0Is4BClzpdYbDavfSegdW8gElbwiTmX3AcnHgAzeNecuSngACgh9Z8AKwy/08DoIanKES57OtnC2G2+eN6+Zbn6F826EOfBPw4llV5BAbZUqmMRfYWusBbtd671BcSBgGBtl7AxPwCgHaDXHLnooAAoIfMzQFwlUhNGvh0CHjBf+vwb9riDw1TtuBnC2G2+WsBfbingJ/H6CxvUQFCe9sVFQpdAC4kEFaDQU6I2MYLmOwpbJo7WJtVwLIIxP5JaqwSAoAe9k4GwCorgOMAmOT9i0NPtU2cbUDOJm+tgK+s0Kf1Eq+fx9NUvahAoZu0RQTDUFBoW48LDOIeaRCo6+PMAxQAdBu/UipbBQQAPfStDoAx7x/uUw0AAWjVzgTW7bQBOZu8UR1cy7loWUYvX7yfAn4ulvcqI0DoJV+hwsi2IJfWc5t6uDDIAUC0Jw6BcW9gbQDwKkyw8Bgo3xHR0bINjIeCRS8qAOhhISsA5Hr/uOcBJ7XbB9x8ytpoWA/AF+2vwJ+N9TPNK1DoL29RPIY2MJfUa5vyaTAYhTpbL6AAoP9YlBqyV0AA0EPjYAAY9f6Z9gSMt9cV3FzL2epVb8Cn+y/gZzsSapJfoDCc7LWEQxug8wHCOAwCAvW9owtCFgJe05Yw4gEMN86kpvwUEAD00NoaAONbv0RX/uIa4M8EZqbr1frjU5ajU73Cnnj8ONYvZR4BRH+z1QoMfaCQUzY6x88UChYA9B9HUkP+CggAemhuBsCU1b/RxR9x75/NSSCmtgvwmRQyXxdPn1mjOs4hgOhn3DzhkAN1tt5BQGAIL6DMAfQbR1I6GwUEAD10TQfAKtu/pK38jXr/XMDNpYxN38W7Z6OW5G0gBQQS3YydBxy6QGHzffwWhnjjXkCbMLAAoNsYkVLZKiAA6KFvSwCssv9fUvg3yfvHATlOHpd+NQLkSTjXZWRImUAKCCzyhMwSDm2gEHm5XkBOGDjfk0BkFTBvtDVuLgFAD9uzADBp9W+1uX9xuMsC9hoN9GBjCeV6jHQpWksFBBqbq58VHCaBoQZAtECHgqMrgpO8gC1WANPCzaDzBcBhRLSkx9D9noiGyjYwHgoWvagAoIeFqgNgwvy/pLl/8b+FAr5GhDwBPY/RLEXrSYFGhsZQgBgFOe0FxBjRoeD4dQ2JOk/SVjACgPX0lJW/LwKAHjZkA2A0/AvAix75puf+Rb2CpjY1ItyJB880KuS6KBBMgXoGSBtAbPIALvDgRRaEpHkBk8LAlfJffkcCgMGGqFQUQIG6AsCrrrqKRowYQR9//DH16NGDLr30Utpqq61SZbrjjjtoyJAh9M4771C3bt1o+PDh1KdPH7asCwFQn7kIr19kAUg0/Gvj/WsEwBOgY48zySgKlE2BMgNkHBA1BNp4AZt7BxsnBPzNN9+ob+pdd91Fn376Ka2//vp02WWX0cYbb5w6hH/44Qc688wzaezYsfTf//6XVlppJRo8eDANGDCgbMO+dO2tGwC87bbb6KCDDiJA4BZbbEHXXHMNXXfddfTaa69R586dWxjmqaeeUnB41llnKejDgD399NPp8ccfp0033ZRlSBYA2nr/ygB/Am+s8SGZRAFRwE+BWoMkYNDHC9isbAN4AP/whz/QtGnTaOTIkdSpUycFdZdccon6Dq+44oqJg2GvvfaiTz75hM4++2xabbXVFDjOmzePNt98c7/BI6WNCtQNAALaNthgAzXwdFprrbWod+/edO6557YQAgMVAPfAAw9Uru2yyy7UoUMHuuWWW4zCIUM6AEbm/+kNnlFAh3+T/hY/C5jVglgmATMX1aSMKCAKNIACLjCpvYE2XsBoGLiZJ7DOAfC7776jdu3a0YQJE2j33XevjKiePXvSHnvsoQAvnh588EHq27cvvfvuu9SxY8cGGIXF6mJdAOCPP/5Ibdq0ofHjxzcL4R533HH04osv0uTJk1uoDq/gCSecoP7phF8qCBt/8MEHiVaCqxr/dJo1a9YC7+JNRNRlwZ8R60UYuFXTf2NMa7jT3r2kv+Ea/r7Mz8UaIdIaUUAUEAUaVIEll/ma2i81i9pT0z+k9vQ1dVjw3x3pywV/03ma9g1cesH1hWVm0dJffk+zvyHqvB7R119/Te3bt89E1YWOiVMCrAI+l2bMmEFLLbVUpa1LLLEE4V88IfyLfJMmTaIddtihcrlXr14q/2OPPdaizNFHH01vvfUWbbTRRjRmzBhq27Yt7bnnnioy17p160z0kUoXKlAXADhz5kzlXn7iiSeauY3POeccGj16NL355pstbL744ovTjTfeSAcccEDl2rhx4+iQQw5pBnnRgmeccQYNG4al9ZJEAVFAFBAFRAE3BTDvvGvXrm6FDaW+//57WnXVVdV8Ot/0y1/+kubMmdOsmqFDhxK+hUkJYVt8W/Et/dWvfqWiaf369aPVV1898TuMqBvAcMcdd1RTsD7//HMCFG6//fZ0/fXX+zZfyhsUqCsAfPLJJwm/NnTCog78qnjjjTcSARBwuP/++1eu3XzzzXTooYcSHqCkFPcA4ldcly5daPr06Zn9mstjBOMX48orr9zil14e9w51j3roA7Soh37UQx/qxRb10o96GVM6avTVV1/R0ksvHer116IefMMQGfNNP//8My2yyCLNqknzACITwBaLN6ZMmUKLLbaYmpbVvXt3ev7559U8wHjaeeedaerUqQpWtUf0zjvvpH333Zfmzp0rXkBfAzYCAOYVAo5rqV3teKijLvKMbRa8+nroRz30QX+s8SIs85gSWwR/RL0qrAd71EMf6uX55gxGwBtstsIKKxDm28OLOHHixBZF+/fvryJ3b7/9duXa66+/TmuvvbYKDcNzKCk7BerCAwh5sAhkww03VKuAdcIgwgqjtEUgmLNw//33V/Lvuuuu6leZ7SKQMn+s6+WlJB+I7F4StjWLLWwVyzZ/PdijHvpQL+9am9EKTyfC0RdccAENHDiwRdFRo0bR8ccfr1b+ItyMhEUke++9t4JGmQdoo7Z93roBQL0NzNVXX63CwBhY1157Lb366qsqTIt5CJgnqGEQ4eKtt95a7f0HSMSgO+2005y2gREAtB94oUvIByK0ou71iS3ctcuiZD3Yox760AgA+NBDDxHCxmussYby6v3lL39RC0CwvVqrVq3olFNOoY8++ohuugkLJ0lBHnbr2GyzzdT8eswBPOyww2ibbbZR329J2SpQNwAImeD9wy8NbAS9zjrrqP2HAHlI2267La2yyipq4YdO//jHPxT0YQm63ggavzy4CXMCAZQY1Emrorj11DpfPfSjHvqAcVAP/aiHPtSLLeqlHzKmav2V4N3/9ttvV9/DDz/8UG3rss8++ygni57fd/DBB9P777/fbEUw5ugPGjRIhYKXWWYZ2m+//dSWMeL942nuk6uuANBHCCkrCogCooAoIAqIAqJAoyggANgolpZ+igKigCggCogCooAosEABAUAZCqKAKCAKiAKigCggCjSYAgKADWZw6a4oIAqIAqKAKCAKiAICgDIGRAFRQBQQBUQBUUAUaDAFBAANBsfK4hEjRqiVxT169FBnBW+11Vappe644w4aMmSI2hFdryzu06dPrsPKps1Yao8l+dOmTVNtxF6KOEJvk002qbQZK7dwako0Yd/Fp59+urD9wmpvHOsXTziwfMkll8yt3Ta2wEr1pHOrd9ttt8omqkWwBXb5xzPx3HPPqefirrvuot69e1fVFP068cQT1bZMnTp1opNOOomOPPLI3Oxg22acRjBy5Eh1ljhWoOLZx/FXv/3tbyttTjoaEsdfhTiCiyuMbb9w7NZ2223Xonpsvrvmmmtyb+uVz7bNSWMeDcA+rxhPSEWwBXaEwLjBqlasYMWxaOeff77aEqVaKsI3w8ugUri0CggAVjGd3lsQH/EtttiCrrnmGrruuuvUkTadO3duUfKpp55ScIiDrAF9+DDifEPsgQRgyiPZtvnAAw9UfcPLCmCEbXTwEsOLFfsmIuEF/Mknn9ANN9xQ6QLOe8Qy/7ySbb8AgMcdd1yL8yd//etf59Vksm3zl19+2ez4pi+++IJ+85vfqDEHGxTFFg888IDasgHHPGGbBxMAvvfee2pbpsMPP5yOOOIIVRbnfWLDdZTPI9m2GZvTAlQBS9gcHmP/wgsvpGeeeYbWX3/9CnRgK6lJkyZVuoDjr5Zbbrk8uqTuYdsvDYA4Hz16ehHajLbnkWzbjH1W8cNNp3nz5qnnAluH6DNp8b+1tgXOte3bty9tvPHGhDYOHjyYXnnlFfW9aNu2baK0Rfhm5GFzuUcxFRAArGIXQBs+cvAE6IRNK+HtSDtdBBuW4gWnE14KHTp0YJ8u4jtMbNscv9/8+fNVe6+44gq1ebaGDpx7fPfdd/s2z7m8bb8AgPiIo921SrZtjrcT3mb8gICXTX9AAIK1tkW0nTgn1ASAJ598Mt1zzz0EL5NO8P699NJLhA9g3onT5qQ2wQuIY61gEyRAB54JeAmLkDj90gCY9Vm0XD04bY7XBc2xXyt+WGCT/yLaAm367LPPaPnll1defb0fbbwvGE+1/mZwbSX56k8BAcAUm9bqfGGfIebS5vj9cDweXlrjx4+nPfbYowKAeOnC6wdvCHZpx+aeyJdHcukXABA7ysOLCajt2bOn8sxq703W7XZpc7xN6667buVUG30NAFhLW8TbyPmA4+MH3S+77LJKcUAjNnz99ttv1QkBeSZOm+Pt+emnn9RG8ghdH3vssRXoQCgcm9xiI3gAP6ZPdO3aNc/uVO7F6ZcGQPTl+++/V2FUbIafFBbOoxOcNsfb8bvf/U6F5f/5z39WLgHGi2QLNAwnYeAsW3gB4QFPSogknXDCCeqfTjjAAD/+PvjggzxMIPdoYAUEAFOMP3PmTAUPCFchPKoTXvCYD4cQSjwBkAAeBxxwQOXSuHHj1Fw0vLCyTi5tjrfpmGOOIRzngzmBeq4cQpk4pxG/tvGrG3McEeLAHLA8TkBx6RfmJ+IFDIjCL2zAB859htcpjwPGXdoctcWzzz6rgAIhx+h8zFrbwgUAu3fvrkLYp556aqU4jmLE1APohAPj80wu0AG4OO+885QXU//wgacfAIv+YYoETi/A/C9Mn8CJBnknTr/w3sIcPMz1xTtpzJgxhOMzAYZpXqos+8Fpc/T+8IavvPLKhPcqfkDoVDRb4Dg0HDEKT+vUqVNTJaz1NyNL20rdxVdAANAAgPhQ4WxhneD5wksTL/okAAQc7r///pVLN998Mx166KHq13bWSUOHTZujbcL8P3zk8DFYb731UpuLlzBg8NZbb1WhmKyTb7/QPnhwEM7HR+7yyy/PuskKbPADwtUWmCuHsvAeVEt528IVAPEjCEdE6YQfVltuuaUKb+c5LxP3t4UOzFWENxnnhe+4446p5pg7d65a+AUvIRa85J1s+6XbB48ayiJMn3eybTOm3lx00UXq+QI8paVa2wI/pCdOnKjmf6+00kpVAbCW34y87S33K5YCAoAp9nAJ4dXane/SZt19THCHBwMT2jfaaCPjKIUXDR9FzO/KOvn0K9o2LELAGZXROZpZtd2nzfAqwSt25plnqoUsppSnLVwAsMwhYHhcAa+YErH77rubTEE77bQTrbbaas3mDRsLBcpgC1P6tvhRO3bs2GZzNAM1yViNTZvhVYO3FVNTECY1pVrZAotTME0DntZVV121ajNr/c0waSjX61sBAcAq9kUIDqESrALWCXNm4NpPWwSCOXQINeq06667qnlz8CLkkWzbjDYhvAX4Q+h3s802MzYTq1Ph3Ro1alRloYixkGcGl35Fb4mPB0KpCAlff/31nq3hFXdtM6YRYJHERx99ZAwl1sIW0d5zPuD4kXDvvfeq1ZA6HXXUUWrxRJEXgeCZHTBggHp2TdvcoF8IqcIDOHDgwMpCEd5ICZOLY4ukO+27776EFeiPPPJImIZY1GLTZj1/sdqcOn3rWtgC7xjAH+a3oq2cqSZYBFLrb4aFuSRrnSkgAFjFoHobD8yRQRgYwIN98zDHByFQrJIFCGkYRMgO3g78ogYkImSECda12AaG22aEfTGnD3NqMCdLJ8z5w785c+ao1Y7YrgNeqffff1/N5Zo+fbryGLRr1y6XR8LWFsOGDVMwi5cw5gAi7IvQPUKP0Tl1WTbets26LdhKCOMKIfZoKoot0A7Mr0TC4o6LL75YLSLAtkDwaCDUC3jF/pJIehsYhLXhhQX0AXDz3AbGts1oG55vzB2NTnPA/m5Y9IH05z//mRA+RZ8//fRT9SMKKz4BKHp1apbjC3Xb9guLC7AABCua4aWG5w/TPrAXXR7TOVzarDU86KCD6D//+U/i/qNFsAW2NsJ7FO/96N5/GC8YN0hF/GZkPUal/uIqIABosA28f4AkzFXCSi6EHvRkaWzci5cpPDY6YS8qQN+7775b2Qg6rxerboNNm9H+pNVmQ4cOVeCH/bfg/XjhhRfU9iOAQHzssaIWk7HzTDb9wqo67GeITXnxAgaooD/R+Zx5tN2mzWjPW2+9pT4eWOGIEFY0FcUWaZsJ9+/fXz0LWPCBHwrIpxPACDbRG0HDK5jnRtC2bU7blFv3Ef3Cnm8I833++edq7z/84MBzgShBXsm2X3iX4YcsAB1QAhAEsGOz8bySbZvRLuwFiHcPgBw/IuKpCLaANzMpYQ9JvY9nUb8Zedle7lMsBQQAi2UPaY0oIAqIAqKAKCAKiAKZKyAAmLnEcgNRQBQQBUQBUUAUEAWKpYAAYLHsIa0RBUQBUUAUEAVEAVEgcwUEADOXWG4gCogCooAoIAqIAqJAsRQQACyWPaQ1ooAoIAqIAqKAKCAKZK6AAGDmEssNRAFRQBQQBUQBUUAUKJYCAoDFsoe0RhQQBUQBUUAUEAVEgcwVEADMXGK5gSggCogCooAoIAqIAsVSQACwWPaQ1ogCooAoIAqIAqKAKJC5AgKAmUssNyiiAqYzSHGSBQ5yxwkoPXv2TO0CdvbHdRyxFTK98cYb6vQAnJe75pprqv/lJH3KwldffaXOoM4yoX04HQYH39cycW3FaePf//53whF+OImlXhLO+t18883pxBNPrJcuST9EAVEggAICgAFElCpqowDOO/7LX/5CgJ1f/OIXqhE4G7VDhw7qWK6pU6dWGob/jyP83nzzTerevbs6Ig75llhiicTGx6EiDayyAkAcEo8jxq6//np1JvMyyyzTop1J925EAJw/fz599tlntOyyy1bGgcuI/OGHH6hr167qDGacx5xlyhOeX375ZXV8I85kXmqppbLsltQtCogCJVJAALBExpKmNlcAMAfv2FNPPaWAD+mBBx6ggQMHKiD48ssvqU2bNurvOKNcTquKAAALb0lEQVQVwIgzUDmp1gC40UYb0e67707Dhg1Lba4AIMeS/Dzjxo1TemNcZZ24APh///d/1KpVK+/mbLjhhnTYYYfRUUcd5V2XVCAKiAL1oYAAYH3YsWF7seKKK9KgQYPor3/9q9Lg5JNPprlz59Kjjz6qDo7fcccd1d932GEHdZj82LFj1X/HQ8DPPvssHXHEEfT666/TOuusQ4MHD6a9995bhYARSkU4OJr69+9PN954IwHC1ltvPVpyySXpuuuuo8UXX5yOPPJIOuOMM1Jt8tNPP9HZZ59No0aNUqC61lpr0XnnnUe77LJLpW3RwkOHDm1RHwBi9OjRze4BDw/AFd6eSZMmKS1ee+01FaLGgfRrrLFGJf+9996r6nz11VepU6dOhP6gz9qTGm88vGzwtsIjudhii9Ghhx5Kn3zyCc2aNasSAoYHDXngQZs9ezYBYi+55BLaeOONVXXaO/nggw8qeyHM3atXL5X/ueeeUyFKADrAF6FYDe/ID72mTZum7o0ysG23bt1UvWmwbtIg3sc999xTaTRixIjKJQ1qm2yyibon+njCCScorU455ZRKO88880waMGBApRz6gf4glLzooovSlltuqcqvssoqSvc42GO84hrGGULQV111FT399NM0cuRIOuSQQ+jJJ59Umv3rX/9Sns4+ffrQueeeS23btlX3RH5oPWPGDGrfvr3yYP7jH/+otAf3e/jhh2nKlCkN+66QjosCokBzBQQAZUSUWoEDDzxQhUofeugh1Q98qE866SR65JFHVIh3+PDh9OOPPyqI+9vf/qbAJQ6AAEZ8eLfffnsaMmSICpUdd9xx9O677yoAXHfddWnChAm0zz77KO8QwmitW7dWH1oAIPLgY3/AAQcobySgAe3ZaaedErXFhxoQcM0119D666+voAp/A4ytvvrqKjwNcAUQ/vnPf1YhYPyLJoDXrrvuqmAV8IG03HLLqbA3AHDTTTel888/X/0NQAqAe+KJJ1Q+tG2//fajyy+/XIHCO++8o7ymaDdgMyldcMEFdM455yjgWXvttemiiy6i22+/XWmm5wBCM0AHQLhLly6EMvfccw+9/fbb1LFjxwoAwlt74YUXKsBDOwDxCMUDghHCB9wAJAGwSHfccYcCdtgBtjr99NMV9GFeJOAqDQCraZDUR4wXeIkRftcJmtx5553Ur18/9UMDGmIM/fa3v1VTCn7/+98rYIMNoOPKK69M3377rYJuaHv88ccrqAbAAnIRjsV4RB2AZIA5EvSZOXOmGocAQeiLsQFdvvjiCzWHD15swDF+NBx77LH0m9/8RpX/97//rTzgY8aMUfng+cY4+OMf/1jpBzzj0BXjJm3aQ6lfBNJ4UUAUsFZAANBaMilQJAWuvfZa5ZHBYoTvvvtOfUjhfYFHBYCDDza8Httss436QGOOVxwA4YmDNwfeE+11AgggXKYXgVSbAwi4is43BIQCjAA0SQnAc8wxx9Cpp55auYwy8JRdeeWV6m8AiN69e1f1JFYLAcP7Ba8n0v3336/AAfrAUwlwATyizzrBMwpwBoQkJXgJAXgayubNm6dgBaFFACDADAAFryhAGAnhS8AMIAhApzWMtg0aoR1R2wBYAXXw/CUlANDyyy9Pr7zyigLgah7ANA3i9WL8oP0YK9H5fwBAtBs/BgCbSJh2gPtrbxrsjx8DAN++ffsqoAf8wpsMcEXSP0Kg1c4776xgO76ARvcDC4qgtU6AT/zgwA8GnR5//HE1pqE77Asv4Ycffkjt2rVL1AzgCWDEPQDnkkQBUUAUEACUMVBqBeBdgtcMITIsBgFowJMGL1rnzp3V3+BNgefqgw8+qPQ1GgIGQL700kvKa6gT/hsQxgHAHj16VMAN5ffaay+1aAMgEE/w+gAWABX4gOsUb4MvAH766afK+4eEPmywwQaq/9AEYUOEoRFO1QkQ8/333yug0BCsr8FrBA/q5MmTFTzqBI/Szz//rAAwDTCQB2AFLTQARtsGDxa8WbivTvBCIkT9/PPPqz8BDuGZRUgU3l60HfknTpxIu+22WyoAVtMgbpePP/5YhcIxFUCHrJEHoAbgxL10gt0AnhrW8XdA1Z/+9CfldQPcA9YA29EEzyDK4IdFNQAE3G2xxRaVohhfGOfRuYDQHfUhxL/SSiup/OgDvMb4B92jdvzPf/6jFj8hP6YcSBIFRAFRQABQxkDpFUDYDR9VwB7AAPOhkDCfC2FfzJWCJ0qH23AtCoDwUAFgXAEwvg0MPHcAJnjD4kkDYBym0AZ4tDBPC8kXAKPbwCBUinAiQtvQAd4kzAnDHMd4godUe7psAFADs4ZMXRZaAIYB4EleVGiEvsMbphPC44BKvfUNQs6wMTyUgDQAIADsrrvuUl5SzoKduAbxfsNDB2CD1xEeOp2SQC3J86o9negLxiLg9eabb26hL6AcPwCqAWB86yEAG6YTREO6umIAPeadwiMLfTHnECFz2BDzBfVWQM8884wKE+vV0qV/6KUDooAo4K2AAKC3hFJBrRVAiAweP+0BxLwyJMxrQ0gMQAiPDPLpFAVAHQJGCA1whIT8CEXqjzE8jPCywAMV3ZIlCQaqASDqTgsBIwx8xRVXsAEQoKIhV/crCbLi8IN+IIwJKOMmgBfgBhCGBOAALMKzqEPACL8DsqMhYISJUQ5zGV0AEPPfsOghGpqFhwxh2pAAiD4BKrFSFu31AUBMS0CoHGCatu0KxiY8dvB06pS2nyHmuWJ86x8HJpvhRxDAD3MTNeTD1gBrTHOQJAqIAqIAFBAAlHFQegUAHQi7Yc4ZIO5Xv/qV6hM8MPDGfPPNNzR9+nTlRUoCQCw8AKjAy3LaaaepDzfmYCHspgEQ8wpRHvdC2BGgiIUZLgCIOV4IcwI89Qrdiy++uLIIBG3keAABEYA7LMZAWwBgACUsAqnmAcQikD322EOtZMUiBniL4AGFBxKLFZISFpTgH0ACHim0F6t3o4tAAE7jx49XeeCZ0otAEMJFGNgFAOHtw3w7zFmEZrCjXg0bGgARwoUHM7p61sUDqBeBAPSxOAQhWrQbi0kwRQH/jQU1+JEBjx1+UMAriDGWtPk4bAPvHeb5HX744SqEj/mF//u//6s83Pfdd5+ao4jwPHTGnECE1VEO4WMk9AMhfxvoL/2LQTogCogCVRUQAJQBUnoFtOcEXi18GHUCDALasF0IYC6a4tvAYH4ZPH4oj5Aj5pxh1W80HIdVmPAmYvsTeBP1NjA2IWC0IboNDOap4X7RbWC4APjWW2+p7VsQfsUCj+g2MNUAEPUDAgEn6B/mlkE7eL8AGEkJHj948QDAAEZseQJvaHQbGMwhhIfwlltuUdCdtg1MtG2cEDAWjSD8CciBxxOLewDeoQEQ29LAownPHIBMg1N8sYYpBIxy8NjBCwgYgxaAQSxIwepneAURioVnD6vG8QMkug1M0ukzCOcC2JEf8/8wprFaGQuJ4BHFDxcAH2yAObHIqz3h+Bt+FMHmer/M0j/00gFRQBTwVkAA0FtCqUAUEAXqRQFAE+ZLRldIl71vWHiCbYzq6Xi7sttE2i8KFEEBAcAiWEHaIAqIAoVQACFg7F2IPf/qJWGqAVYuRzcCr5e+ST9EAVHAXQEBQHftpKQoIAqIAqKAKCAKiAKlVEAAsJRmk0aLAqKAKCAKiAKigCjgroAAoLt2UlIUEAVEAVFAFBAFRIFSKiAAWEqzSaNFAVFAFBAFRAFRQBRwV0AA0F07KSkKiAKigCggCogCokApFRAALKXZpNGigCggCogCooAoIAq4KyAA6K6dlBQFRAFRQBQQBUQBUaCUCggAltJs0mhRQBQQBUQBUUAUEAXcFRAAdNdOSooCooAoIAqIAqKAKFBKBQQAS2k2abQoIAqIAqKAKCAKiALuCvw/tXdViFNrDysAAAAASUVORK5CYII=\" width=\"640\">"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.contourf(X, Y, Phi, colorinterpolation, cmap = colormap)\n",
    "cbar = plt.colorbar()\n",
    "cbar.set_label('Volts ')\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 161,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "Text(33.847222222222214, 0.5, 'Height of the domain (metres)')"
      ]
     },
     "execution_count": 161,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "ticklabels = [i for i in np.arange(0, 2.1, 0.25)]\n",
    "plt.gca().xaxis.set_ticklabels(ticklabels)\n",
    "plt.gca().yaxis.set_ticklabels(ticklabels)\n",
    "plt.xlabel('Width of the domain (metres)')    \n",
    "plt.ylabel('Height of the domain (metres)')    \n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 162,
   "metadata": {},
   "outputs": [],
   "source": [
    "plt.savefig('potential_distribution.png') #Saving the image"
   ]
  },
  {
   "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.2"
  },
  "toc": {
   "base_numbering": 1,
   "nav_menu": {},
   "number_sections": true,
   "sideBar": true,
   "skip_h1_title": false,
   "title_cell": "Table of Contents",
   "title_sidebar": "Contents",
   "toc_cell": false,
   "toc_position": {},
   "toc_section_display": true,
   "toc_window_display": false
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}