PTF-RRLyrae_Phi31-Check.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import numpy as np\n",
    "import pandas as pd\n",
    "from scipy.optimize import curve_fit\n",
    "import os\n",
    "\n",
    "%matplotlib notebook\n",
    "import matplotlib.pyplot as plt\n",
    "\n",
    "#import warnings\n",
    "#warnings.filterwarnings('ignore')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "54\n"
     ]
    }
   ],
   "source": [
    "with open('KIC_6763132_PTF_R_4505_05_noshift.dat','r') as f:\n",
    "    next(f) # skip first row\n",
    "    df = pd.DataFrame(l.rstrip().split() for l in f)\n",
    "\n",
    "print(len(df))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "t  = df[0].astype('float64').values \n",
    "t0 = t - t[0]\n",
    "mag = df[1].astype('float64').values\n",
    "dmag = df[2].astype('float64').values"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "def func(x,m0,A1,A2,A3,A4,A5,A6,phi1,phi2,phi3,phi4,phi5,phi6):\n",
    "    P = 0.5877887\n",
    "    return m0+ A1*np.sin(2*np.pi*x/P + phi1) \\\n",
    "             + A2*np.sin(4.*np.pi*x/P + phi2) \\\n",
    "             + A3*np.sin(6.*np.pi*x/P + phi3) \\\n",
    "             + A4*np.sin(8.*np.pi*x/P + phi4) \\\n",
    "             + A5*np.sin(10.*np.pi*x/P + phi5) \\\n",
    "             + A6*np.sin(12.*np.pi*x/P + phi6) "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[13.33775594 -0.24653029  0.1187528   0.08803022  0.06160114 -0.03307477\n",
      " -0.02091117  1.57733958 11.89880534  0.51505974  1.80302297 -0.04521184\n",
      "  1.6383884 ]\n"
     ]
    }
   ],
   "source": [
    "popt, pcov = curve_fit(func, t0, mag)\n",
    "print(popt)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "phi3: 0.5150597375878729, phi1: 1.577339584099565\n",
      "phi31: -4.216959014710822\n"
     ]
    }
   ],
   "source": [
    "m = 6\n",
    "print('phi3: {}, phi1: {}'.format(popt[m+3], popt[m+1]))\n",
    "print('phi31: {}'.format(popt[m+3] - 3*popt[m+1]))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {
    "scrolled": false
   },
   "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,iVBORw0KGgoAAAANSUhEUgAAAoAAAAHgCAYAAAA10dzkAAAgAElEQVR4XuydC3hU1dX+38mF3O8EGSAxYDAERCX8FYF6oagF6gWKpNbPFkTb2hbxa2s/qvUCKkpLK5X6KVZstfVTG2zRlibVohUULFWClRoIRIjhMkIuTCDJJGQy83/WGQYnk9sJ+5y5ZN79PDytmbPP7PM768x6z9p7rW1xu91usJEACZAACZAACZAACUQMAQsFYMTca14oCZAACZAACZAACWgEKABpCCRAAiRAAiRAAiQQYQQoACPshvNySYAESIAESIAESIACkDZAAiRAAiRAAiRAAhFGgAIwwm44L5cESIAESIAESIAEKABpAyRAAiRAAiRAAiQQYQQoACPshvNySYAESIAESIAESIACkDZAAiRAAiRAAiRAAhFGgAIwwm44L5cESIAESIAESIAEKABpAyRAAiRAAiRAAiQQYQQoACPshvNySYAESIAESIAESIACkDZAAiRAAiRAAiRAAhFGgAIwwm44L5cESIAESIAESIAEKABpAyRAAiRAAiRAAiQQYQQoACPshvNySYAESIAESIAESIACkDZAAiRAAiRAAiRAAhFGgAIwwm44L5cESIAESIAESIAEKABpAyRAAiRAAiRAAiQQYQQoACPshvNySYAESIAESIAESIACkDZAAiRAAiRAAiRAAhFGgAIwwm44L5cESIAESIAESIAEKABpAyRAAiRAAiRAAiQQYQQoACPshvNySYAESIAESIAESIACkDZAAiRAAiRAAiRAAhFGgAIwwm44L5cESIAESIAESIAEKABpAyRAAiRAAiRAAiQQYQQoACPshvNySYAESIAESIAESIACkDZAAiRAAiRAAiRAAhFGgAIwwm44L5cESIAESIAESIAEKABpAyRAAiRAAiRAAiQQYQQoACPshvNySYAESIAESIAESIACkDZAAiRAAiRAAiRAAhFGgAIwwm44L5cESIAESIAESIAEKABpAyRAAiRAAiRAAiQQYQQoACPshvNySYAESIAESIAESIACkDZAAiRAAiRAAiRAAhFGgAIwwm74QLjc5557Drfccgvef/99/L//9/9OX1JdXR1mzJiB3bt3Y/369diyZQuWLVuG2tpaDB48uNOl/+Uvf8Gvf/1r/Otf/8KxY8eQnJyMCRMmYOHChSguLkZsbKxuVG63GzKmp556Ch9//DGio6ORn5+PBx54ANdff712Hu+Yezrpo48+ih//+MedPn7ttdfw2GOPYceOHejo6EBeXh7uvPNOfOtb3zp93G233YZ//vOfOHjwINra2pCTk6N955IlSzpd84kTJ/DQQw/hww8/1M4nrGR8S5cu7TKk1atX48UXX0RVVRWk31lnnYUpU6bgvvvuw7hx404fv2fPHo3hP/7xD3zyySeIiopCYWEhvv/97+OGG27QzU8OrK6uxsiRI0/3sVgsyMjIwKRJk7TvnTx5MhYsWIDnn3++z/POnz9f433FFVdg06ZN3R6/c+dOnHfeeX2eSw4Q7p9++unpYxMTE1FQUKDZyve+9z3IWM+0ORwOXHDBBdi7dy9WrlyJu+666/Sp/Jn4fsdLL72EG2+8sdPXlpeX43/+5380e4iJicEXv/hF/PznP8eoUaM6HdfTeLuzwb6uq729HY888gh++9vfwmazafdQmNxxxx2duoqdybPo3+Li4tDa2trX1wT9c+/4u/st8R+c2IvYntigt8kz99///d/a83f8+HGsWrUKV199NUpKSjS7lj5sJBBoAhSAgSbO71Mm0J0AFAF01VVX4ciRIygtLcUll1yiiRt/AShiTRy3nGPWrFm46aabNNHU2NioCZlnn30WDz74oCa09Lbbb79dO58InyuvvBJOpxMiMIYPH46vfe1r2mnEcYhI8m/3338//v73v2uiVUSFt61YsQI/+clPIOe+9tprNUEqx8j4Fy1adPo4Ob+IIxGc8fHx+OCDD7B8+XKMGDFCE3qDBg3SjhUxceGFF2pi49xzz8XatWt7FIAiDEXMybEiwvbt2wcZz6FDh7B9+/bT43ziiSfwv//7v/j617+Oiy66SLvuP/zhD5pIE+5ybXqbV+yIcJB7IoJXxLScp76+Hu+99x5SU1M1jt4mgkfEhgiQadOmnf57dnY2zjnnHM0JHzhwAP/3f//XZRhybQkJCbqGJ85ZeIqYknb48GFNmMsLhrC+5557dJ2nu4NE8ImYk3P2JAC9THz7jx49GllZWaf/JLZx8cUXa/dYXiREVAl/ebkR0SFMvE0EoAj0H/7wh52GlJubi2HDhvXrWr75zW/i97//vfZyITbw+uuv46c//SkefvjhTly8z+Lf/vY3pKWlnf4OsTMZd6i3/ghAee7EVsUGvU1eLpubm/H4449rz5TY1Lvvvot58+Zpvztiq2wkEGgCFICBJs7vUybgLwAleiLCS6IR4oDGjx+vfUd3P9o/+9nPtOhYTwLls88+0yJfX/jCF3SN89VXX8WcOXM04SORw/40cQhDhw7VnPY777xzuquILHGKEpGRiE5/m0Qiv/vd7+LNN9/UokDSRDhKE+cv0T8RBD1FALv7vl27dmHs2LFaNE4EsjQ5j4gQ/4jSNddcozm1hoYGSIRHT/MKQH8R9NZbb2H69OmQSOczzzzT6VRvv/22JvzWrVvXbcRRnKqM8T//+Y+eIfR4jDhriRZu2LDh9DESxRHBJGLGNzrYny+S6POll16qCVQRAj0JQP+/d/cdYnveSKyID2kyLhGK8mIioszb5H6JcBYBr9JEoMuzJiL47rvvPn0qiVC/8MILWlQ6MzOzx2dR5bsD3bc/ArC7sckLnIjlJ5988vTHr7zyCgVgoG8kv68TAQpAGkTYEfAVgDLV9aUvfUmbwpVImu90l/+PtghEEVxDhgxBRUWF0tSdF5oIEBEv+/fv7zfH3/zmN7j11lu16KFMW3qbTG+LoJTIl94ole+Xex2LTH9edtllXcZ1JgLQ20dPZE8EoohLiWpZrVZdXHoSgC0tLUhKStKiu2+88UbICEAZiIj0jz766IymME+ePImioiJtGnDx4sXa1OmZCkCJvIro+8Y3voE1a9Z0YiTPhtimTNcbLQBF+N17773a1K88V94m0VpZMiDCVqK50lQFlO9FeQW5RMIl+iovJ2effbb2siTRfd8m4l8i6Zs3b4ZMt48ZM0YTxL7Pm8vl0qLIEsmsqanRXlpE3Muz6Z0J8I5fzifRTZllkIj7l7/8ZW061zeq6TsF3NPSD5kyl+fcv8nfZUqYjQQCQYACMBCU+R2GEvD+qMoPr/wwyxSuiAN/seHvdLyOSSKAMqWp2sTxivCUqWSJGMr0jkQ9xBlJBE6m2HpbHzZ16lQtOiUOVNaVeZtMHaWnp2uOSqbWJCIp13bzzTdr0TfvtK7v+GUssgZQpvskWiZrHiVCJusR/ZteASjTsHJeERAyrSj8ZIpZePfWRBRLdEiuq7vv765vTwJQBJZM14qQ8J/K1RsBFCa+TaYd5Z/e1l0EULgIB4mkyhj720Q4yfXI/Zdp7d4EoERZJeIoLzsiGkXoXHfddae/srKyUhM2Mh0vdufbfvSjH+EXv/gFREiLYJHmXV8pgkjEj0Q3RUx1J0h6uy5ZfiBR5qNHj3Y6TCLb8lxIVFCElTTvsyhCUY4X+xRxKmJKxFZ/mtwPsU0RvWKXskZVljRIJNj3pUe4yLS0vPDJC4lwlMikTLlLRNQbXZffArkf8k9emORFUabU7Xa71s93/LJM46tf/aoWuZVlHnKNYpvyMudtvgLQu/RDlmn4TrvLMy5jFgEr903uqzT5u+90fX+48FgS6C8BCsD+EuPxQSfg+1Ytb94yBdzdj6a/AJSomiyclyjJt7/9beXrkOliEWbiiOSfd+2dOCL5Dvlxl79118TBSMKEjMM/aiOOWkSeOHwRgDL1Ko5WHJU4H38hJIv+xcF4mwjSl19+GSkpKd1+t14BKOMQUSlN1g3KdLeMubcmTk2mukQMS2RLb/MKQHHMP/jBD06vAZQ1kJLs89e//lUT2r5NjwDsLgnkv/7rvzQhoLeJQ5fkF0nKkSaRTREuYoeSbCRRoP40EaQSPZREJBFBPYlfEdAiQCT6KXYm0alf/epXWpKHTIeL0Je2detWyMtEd4khsoxA7NA3GivXL2MWAStiTNa9lpWVaQJI7E1vk+iljEls2b9JFE0iWU8//bT2kUTXZD2mrIUTu5Lpb1mOIVOjsuRB1svqbXI/ZK2vCDyveJQ1j3IOmUr3Pk8iUOX+yO+D70uL2JHYhTCR3w9ZYysvbrJ2r6fm/S2RMYuo9jaZShfxJwLb+7LXXRJId9PunALWe8d5nFkEKADNIsvzmkbAKwAlCvLnP/9ZewP/3e9+1yXaZLYAFAfidVwSHZPEE2+TdYGy4F0iABIN8W/iRCSpwD+TWY4T8SdRCH+HLhHBX/7yl5pDk6QPb5OIi0yDiRMScSFCURyjrJ/zjSx6j9crACXJQqYqJXlFoq3iJEWI+mYC+16XiIjZs2dr0SnJbuxPdmxPGa8S3ZGop2/ms/c79QhAGbOIYd8mkSDfjOO+DNU/C9h7/Jm8SEjkUMSfMBRRJK0nAdjduMQuJDNahJe8gMhLglcAynXKC4Jv8wpA/2la/3OLCBJ7FZvWG4ESASiiTmzPv4kAlIii/8uN73EiAuXFRaKP8sKgt8n9EEEsz5xvk3NJ5FzsUJrYjlQJkJcH3ya2KZzkOKkaIKJXhLa8bEgGvZzHu47S28/7W+KfrCUCV/rJvZDvk0YBqPdO8rhgE6AADPYd4Pf3m4DvGkCJoohAkMieRHV8pxzNngKWKTRZnyaRNski9m1SHkWie9u2beuS5ShOXLJKxYn5T0/KOeTv4lAkiUIyBr1NprklYtRXwol8p4hRyVQV0ejf9ApA335SDkZEp5zXGwnz/VySb0T8ScLGn/70p26nqXu70V4RJGuuZKpbpmjFmYtQ60lI6hGARiWBSARJRLBMPYoAl2QYEZeSeKE3YUiuX0S/RA9lKt1bmkjEnExzixARMST21NvUuURJZepT1rFKRLa/U8Dd3QdvdFzWts2cOVPXM9mfKeCeTijjF7ElNqu3dTclL329mbRiF9JEHEsUUqLSvk2yb2UKV34vJBoqolxKH8l/y/Mo7GUqWDh7y0z1tIbR+1skyyS8pVwoAPXeSR4XbAIUgMG+A/z+fhPwzwL2/jhLJqRMj8oPv7RAJIHImiARa/4C0BsZ6C7CJ9NSX/nKV7TpPN+SLl4QIvJE7PkLQBFZErHoKevV219EikQRJWomGcFGCEA5x+WXX65NGfpHfLziTz4Xcag389d3XP2Jgnn7BVIA+mcBiwgU0SaiWESD3jWFemoZylSkZIb31CTCK2vPvNEobxKIJDb432+xFynj45sE0t15JXoogk6igGJ/epqs75MEC//oondJgm8SSE/nk7WL8pLjH83r7fv1CsC+IoDdXaus+9u4caM2bS5JWBLhlCg6BaAei+Ax4UaAAjDc7hjHe7qosq+48maeyhogKWIsIvBMysCIwBHnLmuq9DRxFDLNJjXhJPPR22QqSaZg5Xz+mbxSJkWmUmW6zTfC5+3rjR76O1CJjknpDnHokmjSU/OWTpFok3+tN+lzJhFA6SPrAIWLRF29TYSqXKtEweTv3kQDPezCWQD6vmCIvXnrPfZ13SLa5IXBt8l/S3+ZSpSpSYk6dbdsQPpI9FimkCX6KP28kULpJ4JYEoa8az8lsugtA9NX0pOsCZQserFJ/6LpPV2TtwyM2L8kVnmbXIcsyfAtA9PdOUQoij3JWlGJruptegWgLA2Rly1ZwuBb39Bbpsi7BrC775UpaSncLNcoa3DNEIDyvMhyif5EXfUy4nEkoIcABaAeSjwmpAj0tBOITK3J1Jxk28n6Ofnv3gpBi9PzLQQtpSJEfEkfvYWgJUonC9slCiNTeDK1K4u7ZZF+dwJMnI6sz+sumcMLWZy8TLWK4xJhKw5IohJyvu985zta5FCa1KWT7xEnIoJQ+snUoqwTlPpr8v99y1PImidZLyjTuVIuQ8Syt3ahLIyXSIdEMiXpQLiIeBDxKtEjcYgiKGTxvHdaTKbSZB2YRFpkIby/0JVx+6+l6smQwi0CKNfR1NSkRQBlqtq7A8yZPCg9Xbskw8g9FZEk2bMSjZJ7Ly8+/uVCRFhKxqtkk/oWghb79C0ELaVmZOpYpurFVr1JICLkReR4s171Xoe3ELQ8a/L9ch4RhP6FoCVaKlP7MuXrTQKRsYjNiZ3qLRck49IrAL1ZwHJuKYotz4S8VMk/32QOWf8oEV6xa1n/KPUT5XipnSkvg5KoYoYAlGljKVslSyfkRU24yJIH3wLfeu8DjyOBMyFAAXgm1NgnqAR6EoAyKO+0lEyxyvSS/Hd32zdJ8ojvVnASNZFpNxFmsni9u1IrPV20OGYpKSFToSKw5HsletBdWQ3v+CRK57t7hf+5xXFLdFEyb+X/i2MQZyuiwDvdKE5fBK8IAsmK9DpHcWgiArxFeL3n7imZQT73rmGSrF+ZlpaIplyXZFeK+JD1VTLtKKLO23ra3sv7eX92OAhHASjXKYV9JRNUdj+ROnxn0nq6dhHVspRAonpSBkZsVKJ/YgMivP2bZNNKJE6mU323gvPdkUKiTiJ+RBzJLiEi2sXupXyM/9Zyeq5FBKpkuosglYik2JjYj/9WcBLhFKEn08WSWCQROVlrKEKrP+KvPwJQjpUyO/IceesAigCVZ9O31p6slf3jH/+ovegIZ7F3eQmSZ8sbaTdDAMr45MXK+3IlSzdYB1CP1fEYowhQABpFkuchARIgARIgARIggTAhQAEYJjeKwyQBEiABEiABEiABowhQABpFkucZUARkOsa7f253FyalSfTucjGgwPTzYoShsOytCcf+1Azs5xB6PFzWbfbW+rtjiNHnM+o6+3ueYN0zPnP9vVM8ngTUCFAAqvFj7wFKQNa8dbeLhPdyZW2QrN1i652At1RLb0cFa91TX6JTyqrIelO9zejz6f1eo4/raf9a3+/pz/pOvePrbY2qnEPKDHlr/Ok9J48jARLomQAFIK2DBLohIIvkJVu2pya17saPH092fRAQhsKytxaszEdJSuitSTkUb3FfPTfa6PPp+U4zjpH6d5IU1FuT+pc9bTV4pmOSvXW9Ww92dw75PvleNhIgAWMIUAAaw5FnIQESIAESIAESIIGwIUABGDa3igMlARIgARIgARIgAWMIUAAaw5FnIQESIAESIAESIIGwIUABqHCrXC6XtnWSrE3pawG4wtewKwmQAAmQAAmQgIEEJNtd1ihLUXK9e3kb+PUhcSoKQIXbIHtd5uTkKJyBXUmABEiABEiABIJFQHY8km0RI7FRACrcddk3VfYBFQPSu+epwtexKwmQAAmQAAmQgAEEZNs/CeDY7fZOe6YbcOqwOQUFoMKtEgNKS0uDCEEKQAWQ7EoCJEACJEACASRA/w1QACoYHA1IAR67kgAJkAAJkECQCNB/UwAqmR4NSAkfO5MACZAACZBAUAjQf1MAKhkeDUgJHzuTAAmQAAmQQFAI0H9TACoZHg1ICR87kwAJkAAJkEBQCNB/UwAqGR4NSAkfO5MACZAACZBAUAjQf1MAKhkeDUgJHzuTAAmQAAmQQFAI0H9TACoZHg1ICR87kwAJkAAJkEBQCNB/UwAqGR4NSAkfO5MACZAACZBAUAjQf1MAKhkeDUgJX/87u91AQwPQ2grExwOZmYDF0v/zsAcJkAAJkEBwCQT595z+O0wE4ObNm7Fy5Ups374dNpsN69evx+zZszXjbW9vx7333ovS0lLs27dP25njyiuvxIoVK7RNnntqsgn0fffdp53r6NGjmDBhAh5//HFcdNFFuh8KGpBuVOoH2mxAeTlQUwO0tQFxcUBuLlBUBFit6ufnGUiABEiABAJDwGZDx/vl2L+pBs0NbUjKjMPIy3MRfVHgfs/pv8NEAJaVlWHLli0oKirC3LlzOwlA2YbthhtuwDe/+U1ccMEFOHbsGP77v/8bTqcTH3zwQY/G/NWvfhX/+c9/8NRTT2lC8YUXXsCqVatQUVGB4cOH63oIaEC6MKkfJOKvrAyw2z1iLyEBcDgA+Xt6OjBzJkWgOmWegQRIgATMJ2Cz4cMVZfhnmR0VdiuaOxKQFO3A2HQbLpmZjgt/HJjfc/rvMBGAvhZpsVg6CcDurPX999/HxRdfjE8//RS5EiXyaw6HAykpKXjttdfw5S9/+fSnF154Ia655ho8/PDDuh4CGpAuTGoHyTRBaSlQWQnXqHzU1lk07ScaMHuwG1H7qoCCAmDWLE4Hq5FmbxIgARIwl4DbjW1LS/H2mkrsas9HeoZFm8yRSR37MTcKY6twxe0FmLTU/N9z+u8BKgA3btyIq6++Gna7HampqV0MWqZ/5e9y3PTp009/PnnyZMTFxeHtt9/W9RDQgHRhUjuovh4oKcHB5nRs/XcyDh6UaX8gNhYYMQKYckETRiTZgeJiICtL7bvYmwRIgARIwDQCziP1ePKKEuz+LB2pw5I7LeGWd/3jh5tQaLXjO/8oRsxZ5v6e038PQAHY2tqKL3zhCxgzZow2rdtTmzJlCgYNGoQXX3wRZ511Fl566SV84xvfwOjRo1FZWdltt7a2Nsg/bxMDysnJgUxDdyc0TXuKIunEhw7hsyfWYX15HhqbojF4sCf/Q/JA6uqAtOQOzCmqxtBF8wCdU/eRhI/XSgIkQAKhQuBf6w/hb7etgz0tD3GJ0V2G1drcgYzj1Zixdh4unqNvKdaZXhsF4AATgJIQMm/ePNTU1GhRvN5E2SeffIKFCxdCEkyio6O19YXnnnsuysvLtXWA3bWlS5di2bJlXT6iADzTR7Dvfq7aerxxWwl2fZaOISO7vjEe3d+EwqF2XL22GFHZ5r4x9j1aHkECJEACJNATgb88V4+t3y9B3JB0tMUmdzksrr0JbUftmLKqGNcuMPf3nAJwAAlAEX/FxcVaJvBbb72FLJ3Tgc3NzRBDsFqtkMSQpqYm/PWvf2UEMER+w6r3u/HKraUY3VGJZmt+53V+bjcSD1ehKqYANzw7C3kjWRImRG4bh0ECJEACXQi8+44bv7mhFONiKlGX3vX3PMtehQpnARa+MgtfuNTc33MKwAEiAL3ib+/evfjHP/6B7Ozsfj96kj08cuRI/OxnP8O3vvUtXf1pQLowKR20cyfwq3tsmIEyWOx2NMRZ4Y5PQNogB1KbbWgalIHXMQN3PGLF+PFKX8XOJEACJEACJhJwOoGbp9sw7KMyjMzw/J6fjErAIJcDmW027D+WAdsFM/D7jVbExJg4EEAL/EjZuEiewbO43bL0MrSbROWqqqq0QUq9vsceewzTpk1DZmamVsJFSsPI1O2GDRu09XzeJp/LOj9pkuwxZ84cLFq0SPvv119/HXLpBQUF2rl/9KMfaQkg7777LmIlw0BHowHpgKR4SHU1cOedgPuwDTl15TirtQZxljbEpcQhNj8XR4cXYX+rFQ8+COTlKX4Zu5MACZAACZhKYMMGYOUPbMitL8eYhBokx7ahqT0Oux25qMkqwo8es+Kaa0wdgnZy+u8wiQDKej4RfP5t/vz5kHV5Ernrrkk08IorrtA+ysvLw4IFC7TjpZWUlODuu+/GwYMHNSEpInL58uXaG4HeRgPSS+rMj/v4Y2DhQuDQIWDEcDcGRzUg6mQrak/E42RyJgZnWyC3eMkSICrqzL+HPUmABEiABAJDQETg6l+60L6rCkNO1sAVEwf3uWOw4K7BuOZac6d+vVdI/x0mAjAwJtn/b6EB9Z9Zf3q4XMCKFcCmTYBUg2lpAaSqjwR1JRlbRKEk/v7mN8C4cf05M48lARIgARIIGgHZCeT1jTjy2j/hsh1BnOzsWTgE0VMmA1deGZDC/vTfFIBK9k8DUsLXZ2eZ/r3/fmilX0Tw7d7tKf0i60hkfUhysuezxx/n9G+fMHkACZAACYQCAdnB6aWXgG3bPAVdhwzxjOroUU+R10mTgK99zXQRSP9NAaj0ONCAlPD12VkSQB56yLPRR3SUG/EtDXAca0VzRzycqZlISrZg717gvvvABJA+afIAEiABEggyAUk5kCobr74KyBTPsGGfV3aQzw4f9qzlmT0bkF26LOZNB9N/UwAqPQ00ICV8fXb2RgBHxttwblM50hprENPRBmd0HBrTcrEnmQkgfULkASRAAiQQKgRkLY+s2dm1C8jM9Ozp6dtkn8+GBqCw0LP4W2c5tzO5PPpvCsAzsZvTfWhASvj67CwviFICJmZjGc7JtONEihXOmATEOB1IOWHDJw3p6LhqJhYttzIBpE+aPIAESIAEgkxAFm4/8wxQU+OZ4vXP3JMffZkizs0FvvlNU3d3ov+mAFR6GmhASvj67ux2o+bpUmz7XSU+Qb6W8evdOLyu1o1zUIVJ3yhA7rfN3zi878HyCBIgARIggV4JMAIYUgYSFnUAQ4qYz2AoAE2+M/JjUVKCg83p2PrvZBw86FkjLOuGc3KAyec3YUSSHSguNnWqwOSr5OlJgARIIDII+K0BdA0dhsYTFpxsk+oObqQ1H0ZUNNcABsoYKAAVSFMAKsDT01WmC9at01J8XZZo1NYCskRElo3IZi9R7g5AFgrOm2fqVIGeofIYEiABEiABHQROZQE3/n0bag7H4sDJIXB2AFntR5GW3I7BX74EQ++8kVnAOlCqHkIBqECQAlABnp6upyKASE/31Hzxb01NgJ0RQD0oeQwJkAAJhAqBvZtteGfpRgyt/ies0UcQGw0cixuCHXGTcWjMlVhwt1XLAzGz0X9zDaCSfdGAlPD13VmmC0pLgcpKIL/rxuGQ7QGlRswsrgHsGyaPIAESIIHgE/AW+N9R7sbEvHqkNH+mDaopeShaErJQscuCoiLzd3ei/6YAVHoaaEBK+PR1lumCsjJPpE+yxpyrOsEAACAASURBVGT+V+aB5e8ZGcCMGaZPFegbKI8iARIgARLoi4BvgX/Z2cm/NTZ6dn4ye393+m8KwL5stdfPaUBK+PR3FrFXXu4pHSBbgkgqsJQJkNdEEYVsJEACJEACYUGgU4H/6K5Dlp2e9uwxv8A//TcFoNIDQwNSwte/zjIdLAVCW1uBeNk4MtPUKvH9GxyPJgESIAES0EOAEUA9lAJzDJNAFDhTACrAY1cSIAESIIGII3B6DeAOYOzYzu/x8p5fUeGZ3FmypGudaCNh0X8zAqhkTzQgJXzsTAIkQAIkEIEEZCe41auBujpgxAggKQloboZW63XwYGDxYs9ucGY2+m8KQCX7ogEp4WNnEiABEiCBCCUgInD9emD37s9X9ojomz3bfPEnyOm/KQCVHj0akBI+diYBEiABEohgAjIdLLl9J04AKSme3D7/7YHNwkP/TQGoZFs0ICV87EwCJEACJEACQSFA/00BqGR4NCAlfOxMAiRAAiQQ6QSCVOGB/psCUOnRowEp4WNnEiABEiCBSCYQxBqv9N8UgEqPHg1ICR87kwAJkAAJRCqB3nZ5kv3fZ840tdA//TcFoNKjRwNSwsfOJEACJEACkUggBPZ5p/+mAFR69GhASvjYmQRIgARIIBIJyGa/JSVAejpcicmorfVs8S5bvWdnA1EtTZ7934uLgawsUwjRf1MAKhkWDUgJHzuTAAmQAAlEIoFDh4B163AgOg/v/StaKwDd3g7ExnoKQ0++uAM5HdXAvHnA8OGmEKL/pgBUMiwakBI+diYBEiABEohEAvX1sD1egtKt6ah1JGu7f8gW77LVu+wOkp3QhFlT7LDeyQigmebBvYAV6FIAKsBjVxIgARIggYgk4Opwo+SWUrR+WInYwnxYoiynObhdbrTvqkLChALM+80sREV//pmRsOi/GQFUsicakBI+diYBEiABEohAAtXVwC/usmHq8TJkWOw4nmyFMyYBMU4HUptsqHdn4L3UGfjhz63IyzMHEP03BaCSZdGAlPCxMwmQAAmQQAQS2LkTeOghYOIwG0YcLUdaYw1iOtrgjI5DY1ouDmQXodxmxX33AePHmwOI/psCUMmyaEBK+NiZBEiABEggAglIBPD++6Gt/UtNcSPB0YAYZyucMfFwJGSi8bgFkij84INgBNBE++AaQAW4FIAK8NiVBEiABEggIgm4XMCKFcCOHcDYsYDFZ5mflAisqACKioAlS4CoKHMQ0X8zAqhkWTQgJXzsTAIkQAIkEKEEdu0CVq/2ZP1K6ZekJKC5GVpJGIkMLl4MFBaaB4f+mwJQybpoQEr42JkESIAESCCCCYgIXL8e2L3bUwJGSsGI6Js921zxJ8jpvykAlR49GpASPnYmARIgARKIcAIyHVxTA5w4AaSkALm55k37+qKm/6YAVHr0aEBK+NiZBEiABEggEgnIQr+GhtNhP1d6JmoOWAIqAum/KQCVHj0akBI+diYBEiABEog0AjYbUF7uCfu1tcHWEIdN1bnY0lKEI1FWbRp4zBhgzhxzp4Hpv8NEAG7evBkrV67E9u3bYbPZsH79esyWRQKn2tKlS/Hyyy/jwIEDGDRoECZOnIjly5dj0qRJvT5aTz75pHZeOee4cePwy1/+Epdeeqnux5EGpBsVDyQBEiABEoh0AiL+ysoAux2wWnGgLgFvvOZA9FEb4oem49OxMzUReOCA+Ykg9N9hIgDLysqwZcsWFBUVYe7cuV0E4IsvvoghQ4Zg1KhRcDgcWLVqFdatW4eqqipkZ2d3+8j94Q9/wNe//nWICJw6dSqefvpprF27FhUVFciVRQg6Gg1IByQeQgIkQAIkQAIy7VtaClRWAvn5cLkteOUVYN8+IGeEG4PtVajLKsDe0bPghsX0UjD032EiAH2fHIvF0kUA+j9Z3hu7ceNGTJ8+vdsHT6KDIiifeuqp058XFhZqkcVHH31U18NKA9KFiQeRAAmQAAlEOgGp7FxSAqSnw5WYrGX+vvJHID0NyMoCBrU3IaHVjo/HFcORmIXGRphaDJr+ewAKwJMnT2L16tV4+OGHtQjgYCko5NfkmMTERC1KOEcWGpxqd955Jz788ENs2rSp20e1ra0N8s/bxIBycnLQ2NiI1NTUSH+8ef0kQAIkQAIk0D2BQ4eAdetwIDoP7/0rGrsrgf37PPX/0tOBEcM6kOeuxsdj5+FE6nA4ncCePTBtOzgKwAEkADds2IAbb7wRLS0tsFqtePXVV3HRRRd1a4iHDx/G8OHDtWnlKVOmnD7mkUcewfPPP49KCVF302St4bJly7p8QgHIXzwSIAESIAES6IVAfT1sj5egdGs6ah3JSEwE9u4FoqOBk21uWKM+w/nD6lF90Q04lplv+nZwFIADSAA2NzdryRx1dXV45pln8NZbb2Hbtm3a2kD/5hWAW7duxeTJk09/LIkjv//977FbYtPdNEYAg/Tz5lcyAJmZnfcOCtKw+LUkQAIkQAL6CLg63Ci5pRStH1YitjBf+w2XQtDO2gYURH8Cq/1jtKdkAhdOQGP62dh0vAijplpN2w6OAnAACUB/Exw9ejQWLlyIu+++u4t1nukUsP+JaED6Hnylo/xKBiAuzlMpVDaKtFqVTs3OJEACJEACgSFQXQ384i4bph4vQ4bFjuPJVjgaHEjcuQ3prZ/hWNxQ/Dt+EnJGJyCuwYaorHRMXjYToy8z53ee/nsAC8D8/HzcfPPNkGnb7pokgUi5GMkC9raxY8fi+uuvZxJIYH4P+v4Wv5IBSEgAHA5A/i6LRmbOpAjsmyKPIAESIIGgE9i5E3joIWDiMBtGHC1Hmv1TDD3yIaKP1eMj11hUOvNha83EyFFAYYEblw2rwlmXFQCzZpky40MBGCYCsKmpSUvokDZhwgQ89thjmDZtGjIzM5GVlaXV/Lvuuuu0tX/19fWaqHvhhRe0uoFS30+aZANLwseiRYu0//aWgVmzZo02DfzrX/9amzr++OOPcfbZZ+t6WGhAujCd2UF+JQNkuuB0k8/EHgrM+3E4s0GzFwmQAAmQQHcEJAJ4//2e+n6pKW5kNFRhXMUraE7IQlPyUNTVWdB4HLhhrqcQdFRLk6deYHGxJ03Y4Eb/HSYC8O2339YEn3+bP38+RMDddNNN2no/Wf8nglCSP+69995OSSB5eXlYsGBBp4igCMWf/exn2trB8847T6sfeNlll+k2MxqQblT9P9CnZACSk7v2bzL3x6H/A2YPEiABEiCBngjInr8rVgA7dgBjxwKpJw5hXMU6HEvLg8sSrRV/PuccYO7cU3sBd3QAohrnzQOGDzccLP13mAhAw++8QSekARkEsrvTnCoZgLw8T5qYfzP5x8HEK+OpSYAESCAiCUjSx+rVQF0dcG5WPS6uLkGjJR2HjydDKqldey2Qk3MKjckv+fTfFIBKDyENSAlf7539iobW1nqW/8kyQNncxezpAROvjKcmARIggYglICJw/Xpg9y43CveX4uxWT1bwJZMtn4u/ACzzof+mAFR6CGlASvh673xqDeCRzZXYdCgfBw9Z0N4OxMYCI4a7cflwcxcIm3hlPDUJkAAJRDQBmQ6uqQEc+2wY/EEZsqLsiBpu9bzhexP9MjKAGTNMS/Sj/6YAVHoIaUBK+PrsvHezDe89UAZXvefHITo5AR1NDrgO2WDJysCUZTNMKxHQ5+B4AAmQAAmQgDqBIJX6ov+mAFQyXhqQEr5eO3sXDO/fasPlqeVIb6xBTEcbnNFxsKflml4k1Lwr45lJgARIgAQ6EQhCsX/6bwpApaeQBqSEr9fO/iUDEhwNiHG2whkTD0dCpunbBJl3ZTwzCZAACZBAsAnQf1MAKtkgDUgJX6+dvUVDpdRfd0nAZm8Ubt6V8cwkQAIkQALBJkD/TQGoZIM0ICV8+iOAqV0PbWwEJFH4wQehVYphIwESIAESIAG9BOi/KQD12kq3x9GAlPD12tm/aKj/RiAVFZ7tgJcsOVU01Lyh8MwkQAIkQAIDjAD9NwWgkknTgJTw9dnZt2joiBFAUhLQ3AwcPOjZTmjxYqCwsM/T8AASIAESIIEQJuAtC3PiBJCSAuTmmv9iT/9NAaj0SNCAlPDp6ny6aOhuoLUViI/3iL7Zsyn+dAHkQSRAAiQQwgS6+42XvYDnzDH3N57+mwJQ6bGgASnh0905GG+HugfHA0mABEiABM6IgO8sj2wB553lkX2BzZ7lof+mADwjo/V2ogEp4WNnEiABEiCBCCXQaZ13oRuJrZ+X+mqJz0TFLoup67zpvykAlR49GpASPnYmARIgARKIUALeWq8j4204t6kcaT7F/hvTcrEnuQj7W62mVXqg/6YAVHr0aEBK+NiZBEiABEggQglIrddf3WPDDJQh6aQdx5OtcMYkIMbpQGqTDc2D0vE3zMQdj1gxfrzxkOi/KQCVrIoGpISPnUmABEiABCKUQPV+N165tRSjOyrRbM0H/Gp9JR6uQlVMAW54dhbyRloMp0T/TQGoZFQ0ICV87EwCJEACJBChBFy19XjjthLs+iwdQ0Ym++s/HN3fhMKhdly9thhR2VmGU6L/pgBUMioakBI+diYBEiABEohUAocO4bMn1mF9eR4am6K1rN+4OKCtDairA9KSOzCnqBpDF80Dhg83nBL9NwWgklHRgJTwsTMJkAAJkECkEpC9PEtKcLA5HVv/nawV+G9vB2JjASkJM/n8JoxIsgPFxUAWI4BmmInF7Xa7zThxJJyTAjAS7jKvkQRIgARIwHACIj1KS4HKSrhG5aO2zgKHA0hIALIHuxG1rwooKABmzeq8PtCggdB/MwKoZEo0ICV87EwCJEACJBDJBGw2oKwMsNsBq9Wj/kQFyt8zMoAZMzx/N6HRf1MAKpkVDUgJHzuTAAmQAAlEOgERe+XlQE2NZwGgLASUzYCLikwTf4Kc/psCUOnRowEp4WNnEiABEiABEgDcbrjqGnDok1acaI9H4ohM5J5tQVSUeXDovykAlayLBqSEj51JgARIgARIALIn8Pr1wO7dQGsrEB8PjBkDzJkDFBaaA4j+mwJQybJoQEr42JkESIAESCDCCYj4W73aU/pFsn+TkoDmZuDAAWilYRYvNkcE0n9TACo9ejQgJXzsTAIkQAIkEMEEXC5gxQpgxw5g7Ngum4GgosKzFHDJEhg+HUz/TQGo9OjRgJTwsTMJkAAJkEAEE6iuBu6/3xPpS03tCqKxEZBygQ8+COTlGQuK/psCUMmiaEBK+NiZBEiABEggggns3Ak89JCn3F90lBsJjgbEOFvhjImHIyETzg4L9uwB7rsPGD/eWFD03xSAShZFA1LCx84kQAIkQAIRTMAbARwZb8O5TeVIa6xBTEcbnNFxaEzLxZ7kIuxvtTICaJKNcCcQBbAUgArw2JUESIAESCCiCcgawF/dY0PMxjKck2nHiRQrnDEJiHE6kHLChk8a0tFx1UwsWm7lGkATLIUCUAEqBaACPHYlARIgARKIbAJuN2qeLsW231XiE+RjcLZFqwMt9aDrat04B1WY9I0C5H7b+O3g6L85Baz08NGAlPCxMwmQAAmQQCQTkAyPkhIcbE7H1n8n4+BBoL0diI31lISZfH4TRiTZgeJiICvLUFL03xSASgZFA1LCx84kQAIkQAKRTODQIWDdOi3F12WJRm2tZytg2RI4OxuIcncAslBw3jxg+HBDSdF/UwAqGRQNSAkfO5MACZAACUQygVMRQKSnA8nJXUk0NQF2RgDNMhGuAVQgSwGoAI9dSYAESIAEIpuA2w2UlgKVlUB+Plxuy+dRwHg3shurEFVYAMziGkAzDCUsBODmzZuxcuVKbN++HTabDevXr8fs2bNP81i6dClefvllHDhwAIMGDcLEiROxfPlyTJo0qUdmfZ1TD2zTBKA8FA0Nn2+KmJnZuUS6nsHxGBIgARIgARIIdQI2G1BWhiOVdmzdb0X1kQRY2hwY4rQh5ewMjLlzBkZfZjX8Kkzz34aP1LwThoUALCsrw5YtW1BUVIS5c+d2EYAvvvgihgwZglGjRsHhcGDVqlVYt24dqqqqkC0LCbppfZ1TD3JTDEgehvJyoKbGkwolKVG5uZ79cKzGPwR6rpPHkAAJkAAJkIBZBPZutuH1FeUYZKtBdmobohPj8FlcLt53FiEmx2rKfsCm+G+zAJl03rAQgL7XbrFYughAfzbeG7tx40ZMnz69T3R6ztndSQw3oFNvQtqaBxF7shJWVsTK32WNxMyZFIF93k0eQAIkQAIkEC4ETu8HXO7GxJENiO34fCcQNyym7QdsuP8OF+A+4xxwAvDkyZNYvXo1Hn74YS0COFg2Geyj6RWAbW1tkH/eJgaUk5ODxsZGpHa3kWFfX+z7ud9aCFgsn38qn1VVefbLMWEtRH+GyWNJgARIgARIwCgCwdoPmAIwDLOAexJrGzZswI033oiWlhZYrVa8+uqruOiii3TZqF4BKGsNly1b1uWchgjAIGZD6YLEg0iABEiABEjAYALB2g+YAnAACcDm5mYtQaSurg7PPPMM3nrrLWzbtk1bG9hX0ysATY0A+tRDQnR01yF3mFcPqS8+/JwESIAESIAEzCAQrP2AKQAHkAD0N8zRo0dj4cKFuPvuu/u0Wb0C0P9EhhqQTwTQlZjctSBmi3n1kPoExANIgARIgARIwAQCwdoP2FD/bQKXQJxywK0B9ELLz8/HzTffDJm27auFhAA8tQbwyOZKbDqUj4OHLKe3xBkx3I3Lh1fhrMu4BrCve8nPSYAESIAEwohAkPYDpgAMkwhgU1OTltAhbcKECXjssccwbdo0ZGZmIisrS6v5d91112lr/+rr6/Hkk0/ihRde0OoGjhs3Tusn2cBz5szBokWLtP/u7Zy5UnZFRzPagCQV/r0HyuCqtyNquBXRyQnoaHLAdcgGS1YGpiwzpx6SjkvlISRAAiRAAiRgPIEg7QdstP82Hoz5ZwyLCODbb7+tCT7/Nn/+fKxZswY33XSTtt5P1v+JIJTkj3vvvbdTEkheXh4WLFhwOiLY2zmfe+45XeSNNCBvKvz+rTZcnlqO9MYaxHS0wRkdB3taLjYdL8KoqVYsWQJERekaHg8iARIgARIggdAi4L/RgZQ6e+WVgO8HbKT/Di3A+kcTFgJQ/+UE9kgjDahTKnyKGwmOBsQ4P6+H1HjcAnlRevBB7TlhIwESIAESIIHwItDdRgdS41YcoDi2AO4HbKT/Dq+b8PloKQAV7pyRBtQpFb6bJGCnE9izB7jvPmD8eIVBsysJkAAJkAAJBJpATxsdHD7s2flKhODEiZ23PTWxBq6R/jvQKI36PgpABZJGGhAjgAo3gl1JgARIgARCl0BfGx3I9qeyA1ZODjBsWOddsDIygBkzDN8Fy0j/Hbrgex8ZBaDCnTPSgLgGUOFGsCsJkAAJkEDoEtCz0cGnnwJnn+0RgrLjVlwcIAmZRUWGiz8BZaT/Dl3wFICm3RujDcg/CzgqKQGuZk8WcFRWOiYvm4nRl1lNux6emARIgARIgAQMJ6B3o4MbbvBE/1pbgfh4IDOz85SwgQMz2n8bOLSAnYoRQAXUhhpQL3UAc0a4cdkw1gFUuFXsSgIkQAIkECwCeiKAEvkrLoYrI0tbEnjiBJCS4gkCmlH5wlD/HSyuit9LAagA0FAD4k4gCneCXUmABEiABEKWQF9rAKXOb0EBdo2chfXrgYMfNaCjpRXRifEYcX4m5nzFgsJCY6/OUP9t7NACdjYKQAXUhhqQ3hD5vHnA8OEKo2ZXEiABEiABEggwgZ6ygOXvGRnYe84MrF0LpO0rx7lxNUiIboOjIw572nLROKoIC+62GioCDfXfAUZp1NdRACqQNNSA+hEiR1aWwqjZlQRIgARIgASCQMBmg2t7OY59WIO2422IS41DxoW5wIQi/OpXQMzGMpyTaceJFCucMQmIcTqQcsKGTxrS0XHVTCxabjVsOthQ/x0ElEZ8JQWgAkVDDcgnRO4alY/aOgukQLqsh80e7EbUPk+IHLNmmbYoVgEFu5IACZAACZBArwR27QLW/8ndZYr34ouBDx8txeiOSjRb87vUAkw8XIWqmALc8Ows5I20GELZUP9tyIgCfxIKQAXmhhuQzYbDz5ahYqsdlcetaHYnIMniQEGqDWOnZmDYQuNrISlcfmC6+m8bZGJWWGAuiN9CAiRAApFHQMTf6tVAXZ2n3F9SEtDcDBw4AKScrMeEvSU4qyAdzvjkLnCiHU04XmPH1MeLMfZSY2bADPffYXhLKQAVbprRBiQPyHOP2gK2BkLh0gPTtbttg0ysCxWYi+K3kAAJkEBkEfDWud2xAxg7tutmH/veOYTzK9chsygPCcldt8JyNHUgzlaNy5+Yh5xLjFkDb7T/Dsc7SgGocNeMNKBOD0ihG4mtn+8F3BKfiYpdFq0e5pIl5qTEK2Awp2tvC4Zly6CZM00pDmrOxfCsJEACJBC5BDrtdJXalcNJWz2S/lqC6Kx05BQmw+IzyyuTQEf3N6FwqB1Xry1GVDYjgEZZEgWgAkkjBWBfD0hjIyB5Ig8+6Nkze0A3nSUDuB5yQFsBL44ESGCAEOhzr/t2Nyx/K0W+sxK2pHxkDbagowNoaQGON7oxLq4Kk+YXIPfbxq2BN9J/h+ttogBUuHNGGlCXB8TtRoLj8yjgidhM7NlrwX33AePHKww6HLqyJmI43CWOkQRIgAR0EdAT4HAesOH2vDLs2WbHthor6poTEO92IC/eBuvYDFx83wxDd8Iy0n/rghCCB1EAKtwUIw3I9wEZZrHBaitHWmMNYjra4IyOw+GYXOyMLcIPf24d+BHAUzURD0Tn4b1/RePgQaC9HYiNBUaMACZf3IGcjmqANREVrJddSYAESCAwBPpaA1hR4dny9yuTbXjjp+Ww1NQgLb4Ng1Li0JiWi/c7ihCTY8XixTCsFqCR/jswFI3/FgpABaZGGpD3Adm/1YZZUWVIbLPjePKpWkjtDjRX2TB4dDpmrZ6JqOEDfD/g+nrYHi9B6dZ01DqSMXiwZ1tI2R5SMsiyE5owa4od1juLwZqICgbMriRAAiQQIAK+WcDyIu/NApYXfPmNX7QIeO01YEe5GxNHNiC2oxXOmHg4EjLhhgVekWjUOngj/XeAEBr+NRSACkiNNqBdFW68eVcpkg9WwpmXj7h4C9raPKInNcWNuRdExn7Arg43Sm4pReuHlYgtzIcl6vMVwW6XG+27qpAwoQDzfjMLUdHG1IRSMAN2JQESIAES0EFAqwO4Hti92/NCLy/2ssXb7Nmemrf33+8Rg6ndJIoYvQ7eaP+t4/JD7hAKQIVbYrgBnYp8/XN3OvbXJp+e9pSaSZdcAuRkNAGnNsweyJEvmQ7/xV02TD1ehgyLTyTU6UBqkw317gy8lzojMqbDFeyTXUmABEgg1AjIbFdNDXDiBJCSAkhlr6gooNM6+KjOa+AlCujssGDPHhi2Dt5w/x1qoHWMhwJQB6SeDjHcgE6tfXPl5uFofTTkP6XJ1r9DhgBR7g5A1NEAX/vm/SGYOMyGEUc7r4WU9SAHsotQbrMa9kOgYALsSgIkQAIkYAAB7zr4kfE2nNvU9Xd/T3IR9rdaDauEYbj/NoBBoE9BAahA3HADOpX9erA5HVv/ndwl+WHKBU0YkWQHigf22rdOGWMpXd8EG49bIqckjoJ9sisJkAAJhAsBiQz+6h5bwPYDNtx/hwton3FSACrcNMMNyO1GzdOl2Pa7SnyCfAzOtnye/FDrxjmowqRvGFsLSeHyTeuqN2PMqMXApl0IT0wCJEACJKCPQDf+Ly4OnnXwJvg/w/23vqsMqaMoABVuh9EG1PcbUAY6rpqBRcut2pqJgdz6yhgzshzAQObIayMBEiCBsCDQywyYrIOffL6xM2BG+++wYOw3SApAhbtmtAEFeg2EwqUHpGtvGWOSOcZGAiRAAiQwQAicWgMvhW5dlmjU1gIOhyc7ODvb+DXwRvvvcLwLFIAKd81oAwp0FpTCpQesa08ZYwEbAL+IBEiABEjAfAI+O0AhObnr9zUZWwXDaP9tPiDjv4ECUIGp0QbUZbscv+3gPjuZifoGi2FZUAqXzq4kQAIkQAIkYByBAO8Bb7T/Ng5E4M5EAajA2mgD8k1+uDjHhmGffZ4K3x4dh52NuUi+tCgi1gAq3BZ2JQESIAESCEcCNhtQVuapd2u1euZ/ZR5Y/p6RAcyY4fm7Ac1o/23AkAJ+CgpABeRmGJCse3vuURtyK8owPNGOtkwrmjoScPyoAyOibZg6Kx3Dbp1p2EOgcPnsSgIkQAIkQALGErDZ4NpejmMf1qDteBviUuOQcWEuoiYWGer3zPDfxoIw/2wUgAqMTTGgU6nwe/5ciY9a8tHutCA2FtB2A5nkRk5bFVBQAMyaBVi4DZrC7WNXEiABEiCBECOgJf/9yY2DHzWgo6UV0YnxGHF+JuZ8xaJtG2dUM8V/GzW4AJ2HAlABtCkGdGohrCs1HbWO5M5ZUFL6xeCFsAqXz64kQAIkQAIkYBgB3/JfEvRISgKam4EDBzx7BBtZ/ssU/20YicCciAJQgbMpBhTgVHiFy2dXEiABEiABEjCEQKcNAArdSGxtQIyzFc6YeLTEZ6JilwVFRYBRGwCY4r8NIRG4k1AAKrA2xYC4HZzCHWFXEiABEiCBcCQQ6Dq4pvjvMANPAahww0wxIN/tcNznYHh6ExJj2tHijMUhezLOsXwSEdvBKdwWdiUBEiABEggzAlIHV/YCnoEyJJ2043iyFc6YBMQ4HUhtsqF5UDr+hpm44xErxo9XvzhT/Lf6sAJ6BgpABdxmGJB3O7jUspdwnqUCzpMudLiA6CggZlAU/uMehxOzbmQpGIX7xq4kQAIkQAKhRaB6vxuv3FqK0R2VaLbmd05ydLuReLgKVTEFuOHZWcgbqZ4AaYb/Di2ifY8mLATg5s2bsXLlSmzfvh02/noSbwAAIABJREFUmw3r16/H7NmzT1/d0qVL8fLLL+PAgQMYNGgQJk6ciOXLl2PSpEk9Enj00Ufxpz/9Cbt370ZCQgKmTJmCn/70pyiQDFudzQwDkjD4L+6y4ZLql5B+uAJtDhdcbiDKAsQlRME+bCz+mfc1/PDnVtkxh40ESIAESIAEwp6Aq7Yeb9xWgl2fpWPIyORORS6kRvTR/U0oHGrH1WuLEZWdpXy9Zvhv5UEF+ARhIQDLysqwZcsWFBUVYe7cuV0E4IsvvoghQ4Zg1KhRcDgcWLVqFdatW4eqqipkyyaC3bQZM2bgxhtvxEUXXQSn04mf/OQn2LlzJyoqKpAkqUc6mhkGtPMjN0puKUXGkUpUuc/B4IQmxEe1o9UVizpHMs7BJ7APLUDxb2dh/Pnqb0E6LpOHkAAJkAAJkIC5BA4dwmdPrMP68jw0NkVrWb9xcUBbG1BXB6Qld2BOUTWGLpoHDB+uPBYz/LfyoAJ8grAQgL5MLBZLFwHoz8x7Yzdu3Ijp06frQlpbW6uJyE2bNuGyyy7T1ccMA6reXo//u74Edc50xA/uuh9ia10TBsfY8V+vFSNvovpbkK4L5UEkQAIkQAIkYCaBXhIgpSTM5PObMCLJDhQXA1nqvs8M/20mHjPOPeAE4MmTJ7F69Wo8/PDDWgRwsLxG6Ghy7OjRo7Uo4HnnnaejB2CGAX269RBK5q3DPlce0rOiu4zDXt+BUVHVKF43D2dPUX8L0nWhPIgESIAESIAEzCTgsxewa1Q+jtZaIFXRpA0f5saQ41WIKjRuEwQz/LeZeMw494ARgBs2bNCmdFtaWmC1WvHqq69q07t6mtvtxvXXX49jx47hnXfe6bFLW1sb5J+3iQHl5OSgsbERqamper6qz2Mq3qnHm98uweGWdDR2JCMxEYiJAZxOoKUFSItuwrBEO6Y/XYyxl6q/BfU5IB5AAiRAAiRAAoEgcGov4COVdmzdb0X1kQRY2hwY4rQh5ewMjLlzBkZfxr2AjboVA0YANjc3awkidXV1eOaZZ/DWW29h27Zt2rRuX+173/se/vrXv+Ldd9/FiBEjejxckk2WLVvW5XMjBaA3E2rYiUpUtOXj+AkLOjqA6GggNcWNsYOqcDjVuEyovtjwcxIgARIgARIIFIG9m214fUU5BtlqkJ3ahujEOHwWl4v3nUWIybEathsII4DAgBGA/sYp07kLFy7E3Xff3avd3nHHHVq0UDKNR44c2euxgYgAesvAxGwsw6gMO45GW9HiTkCixYEhHTbsO5aBjqtmsAxMoH6N+D0kQAIkQAKGEBD/VlMDnDgBpKQAublAlGxxeqqd3g2k3I2JIxsQ2+HZCcSRkAk3LKiogGG7gVAADmABmJ+fj5tvvhkSteuuybSviD8pKfP2229r6//628wyINkP8blHbUjbV47Rg2qQGN2Glo447D2Zi8ZRRVhwt9XQTbH7e908ngRIgARIgAT6Q0D82vr1wO7dQGsrEB8PjBkDzJmD0/7MuxuILN3vblVVYyMguSIPPgjlMmhm+e/+MAn2sWERAWxqatISOqRNmDABjz32GKZNm4bMzExkZWVpNf+uu+46be1ffX09nnzySbzwwgta3cBx48Zp/SQbeM6cOVi0aJH239/97nch5WNee+21TrX/0tLStLqAepqZBqQ9LH9y4+BHDehoaUV0YjxyLsjE7DkWij89N4fHkAAJkAAJhAQB8WerV3vKuUhGr1Raa24GDhyAVu5l8WKPCJTdQB56CJByvLLsyb/JWvg9e4D77oPybiBm+u+QgK5jEGEhACVCJ4LPv82fPx9r1qzBTTfdpK33k/V/Iggl+ePee+/tlASSl5eHBQsWnI4ISjmZ7tpvf/tb7Tg9zWwD6itcrmeMPIYESIAESIAEgkXg9LTuDmDs2C4bfHSa1pXp4fvv94hCRgDNv2NhIQDNx3Bm32C2ADyzUbEXCZAACZAACYQGgf5M68qawBUrgB0iFgvdSGxtQIzTsw6wJT4TFbssXANo4G2lAFSASQGoAI9dSYAESIAEBjyB/k7rBmoNPP13GCaBhNLTQgMKpbvBsZAACZAACYQagf5EALX97W02HH62DBVb7ag8bkWzOwFJFgcK0mwYOzkdw26dCVjVawHSf1MAKj0rATMgqZDe0PB56lRmZueFFEpXwc4kQAIkQAIkYA6B/qwBjLK4gdJSoLISshtIbZ0FDgcgeZnZg92I2lflyRCZNUvZBwbMf5uD1ZCzcgpYAWNADMhmg2t7OY59WIO2422IS41DxoW5iJpYZMhbkMLlsysJkAAJkAAJ9EnANwtY9lrwZgEfPNg5C1ir8VJSAqSnA8nJXc/b1ATYjdkPOCD+u08ywT2AAlCBv+kG1FMoPNWGsVOMC4UrIGBXEiABEiABEuiTQHd1AKX0y+zZn9cB1Db/XbfOU+Svuzowsi2WzCnPmwcMH97nd/Z2gOn+W2l0gelMAajA2VQDcrtR83Qptv2uEp8gH4OzLVrhTCmgWVfrxjmowqRvFCD32+qhcAUE7EoCJEACJEACugj0WdqMEUBdHI06iAJQgaSZAtBVW483bivBrs/SMSQvCR2NJ+Bua4clLhbRaSk4Wt2MwqF2XL22GFHZWQpXwa4kQAIkQAIkEAIEZL37qTWAyM/vWjRQNoTgGkDDbhQFoAJKMwXggX8ewqZF69ASm4r4zz5F3PGjiHa1oyMqFm2pQ9A69Gwkth/H5U/MQ84laqFwBQTsSgIkQAIkQALGEbDZgLIyz1o/yfaVDBDJBJG/Z2QAM2YYsv7dTP9tHAxzz0QBqMDXTAOqeKce7y18BmkN+wGnE61JGXBFxyGqow3xzceAmBg0Zo7E5N98E2MvZQRQ4TayKwmQAAmQQCgRELFXXg7I1iBtbUBcHCBVoouMS34003+HEsrexkIBqHCnzDSgfVUuvD51GUY3fYhDmecBUVGfj9TlwvCG/2Bv8oX40pYHMCrf5zOF62FXEiABEiABEggJAiaXPzPTf4cEPx2DoADUAamnQ8w0oOrt9Xjt2l8ju2k/UhM7cCImA86oOMS42pDiPIbjjhjUJuXh+r98C3kTGQFUuI3sSgIkQAIkEGEEzPTf4YKSAlDhTplpQLs2etYAftaSiuyWapwVXYs4Szva3LE40pGN2sQ8DE30rAEsvJJrABVuI7uSAAmQAAlEGAEz/Xe4oKQAVLhTZhpQzY56vHV7CVoGpeNwYxKc9hOIcrbDFROLmPQUWFObkdRuxxfXFCN3AiOACreRXUmABEiABCKMgJn+O1xQUgAq3CkzDcjV4UbJLaVo/bASsYX5aHFY0N4OxMYCiQlutO+qQsKEAsz7zSxERVsUroJdSYAESIAESCCyCJjpv8OFJAWgwp0y24D2brbhvQfK4Kq3I2q4FVFJCXA1O+A6ZIMlKwNTls3A6MvUN8VWQMCuJEACJEACJBB2BMz23+EAhAJQ4S4FwoBEBG5fW47Wys/T4RPG5KLo1iKKP4V7x64kQAIkQAKRSyAQ/jvU6ZoqAOfMmQOLpev0pPwtPj4e+fn5uOmmm1Aglb3DsAXKgGQ6+OBHDWiub0VSVjxGnJ/Jad8wtBcOmQRIgARIIDQIBMp/h8bVdj8KUwXgggUL8OqrryI9PR0TJ06E2+3Gjh07YLfbcfXVV+Pf//43qqur8eabb2Lq1KmhzKnbsdGAwu6WccAkQAIkQAIkAPpvwFQB+OMf/1iD/MQTTyDqVCFjl8uFO++8EykpKVi+fDluv/12fPzxx3j33XfDziRpQGF3yzhgEiABEiABEqAAhMkCMDs7G1u2bMG5557bydz27NmDKVOmoK6uDjt37sSll16qRQXDrQVUAJpcFT3c2HO8JEACJEACJHCmBALqv890kCb3MzUCmJGRgeeffx7XXXddp8v485//jPnz5+PYsWPYu3cvLr74Yu3/h1sLmAEFYF/EUGTvcnm2gjxxAkhJ8WwF6bsjXiiOmWMiARIgARIIfQIB898hjMJUAbh48WK89NJLuOeee3DRRRdpCSH/+te/8Mgjj2jJH48//jjWrl2L5557jlPAPRmJiL+yMrga7KiLtaLFnYBEiwOD222IykwHZs4ErAOvFMyuXcD69cDu3UBrKxAfD4wZA8yZAxQWhvATxaGRAAmQAAmEPAEKQJOngDs6OrBixQptDeCRI0c0gzjrrLNwxx13YMmSJYiOjkZNTY22PnDEiBEhbzD+AzTdgGTat7QURzZXYtOhfBw89Hkx6BHD3bh8eBXOuqwAmDUL6CbbOuyAnhqwiL/Vq4G6OiAnB0hKApqbgQMHgMGDgcWLKQLD9d5y3CRAAiQQCgRM99+hcJF9jMHUCKDvdwtsaampqWGARd8QTTeg+nrYHi9B6dZ01DqSNfEjkTCJiIk4yk5owqwpdljvLAayBsZ2cDLtu2IFsGMHMHZsZ10reriiAigqApYs4XSwPivlUSRAAiRAAgEP4IQB8oAJwDBg0e8hmi0AXQcOYePt67DjWB5G5EZ3CvKJGDpY04GijGpMXzMPUTnD+z3+UOxQXQ3cf78n0tfdu0JjI1BfDzz4IJCXF4pXwDGRAAmQAAmEOgGz/XeoX7+Mz3QB+Morr6CkpESb6j158mQnJuXl5eHAqMcxmm1ANTvq8dbtJdpav5j0ZMDtxqCTJxDd0Y6O6Fi0tFjgOtaIL64pRu6EgREB3LkTeOghQGqDR0d3Re90Anv2APfdB4wfH9bmw8GTAAmQQGQRCKFqFmb773C4saYKwNWrV+MnP/mJlvH7zDPP4JZbbsEnn3yC999/H9/73ve0OoDh3Mw2oJ0fubHhu6W4IKESrYmZyDi2H0ktRz0CMCoW6GjHPwddjmlrb8b487vuuBKObBkBDMe7xjGTAAmQQB8EQqyahdn+OxzswVQBOGbMGDzwwAP42te+phV+lp0/Ro0ahfvvvx8NDQ1ackg4N7MNSMTQL+6yYXrtS8iv3wanJRbNyUPgBpDSdBTOtnbsSpmEy5/6GnIuHhiZwFwDGM5PBMdOAiRAAt0QOFXNAlLvV6pWJCQADgcgf08PTjULs/13ONiBqQIwMTERu3btwtlnn40hQ4bg73//Oy644AKt9t8ll1yCelnMFcbNbAPSxNCjbiS/9gKmtG+CKyoGzjYn2l2xaEnORlXHKBQObcAXv1OAqC8PnExg3yxgSQ73ZgEfPMgs4DB+XDh0EiCBSCRwqpoFKiuB/PyumX1VVZ41PwGuZmG2/w6HW22qAJRon6wBLCoq0uoA3nbbbfj2t7+NN954AzfeeKMWBQznFggD2vNePT74nxJU29Nwss0NZ0s7WpyxsHekICPDgpuub8LEUXageOBkAotNdFcHUOr/zZ7NEjDh/Mxw7CRAAhFGQAI9JSWeSF9ycteLb2qCFhkMsA8LhP8O9TttqgAUwZeTk6NNA69ZswY/+MEPMHXqVHzwwQf4yle+gmeffTbU+fQ6voAY0KFD+M8D6/Dcpjw02KMRGwvtX2KipySMdUgH5hRVY+iiecDwgZEJ7IXOnUDC+vHg4EmABEgAOHQIWLfOU7ahu8y+jg5A1jvNC6wPC4j/DvH7b6oAdLlckH8xMTEahnXr1uGdd95Bfn4+vvOd7yBWlEwYt0AYkKu2Hm/cVoJdtnQknZUMyYL1CkBBd3R/EwqH2nH12mJEZQ+MTOAwNgkOnQRIgARIwJcAI4Ahaw+mCkC56tbWVnz00Uc4evSoJga9TbaFu/baa0MWjJ6BBUIAVu9345VbSzG6oxLN1q7rJxIPV6EqpgA3PDsLeSMHRiawHvY8hgRIgARIIAwIcA1gyN4kUwXg3/72N3z961/vNtlDBKBsFRfOLRACUOri/eoeG2agDEkn7TiebIUzJgExTgdSm2xoGpSB1zEDdzxiZV28cDYmjp0ESIAEBiqB3rKAMzKAGTMCvqd9IPx3qN9OUwWgTPV+6Utf0sq+yB7AZ9o2b96MlStXYvv27bDZbFi/fj1mSzbAqbZ06VK8/PLLOHDgAAYNGoSJEydqNQYnTZrU41c+9dRTkH/VsvYAwLhx47Rxzpw5U/cwA2FA3rp4I+NtOLepHGmNNYjpaIMzOg6NabnYk1yE/a1W7oyh+67xQBIgARIggYATYB3AgCPv6wtNFYCy7++OHTtwzjnn9DWOXj8vKyvDli1btGziuXPndhGAL774olZmRrKOHQ4HVq1apa03rKqqQnZ2drfn/stf/oLo6GhtPaK0559/XhOZMl4Rg3paIARgp7p4hW4ktjYgxtkKZ0w8WuIzUbHLwr1x9dwsHkMCJEACJBBcAtwJJLj8/b7dVAG4cOFCLev31ltvNeyiZerYPwLof3KvMNu4cSOmT5+u+7szMzM1Eah3vIEQgDJ41sXTfQt5IAmQAAmQAAn0SSBQ/rvPgQTxAFMFYEtLC+bNm6dF4caPH98l63fx4sX9vvS+BKDsNyxb0D388MNaBHDw4MF9foesRZSIoWxZJxHAsWPH9tlHDgikAXnr4sn/SvnEqChPPbz582X6WtdweRAJkAAJkAAJkECA/XeoAjdVAK5duxa33347EhISkJWVBRFv3ib/f9++ff3m0pMA3LBhg1ZcWkSn1WrFq6++qhWf7q3t3LkTkydP1jKVk5OTIVPJs6QaeQ+tra0N8s/bRABKncPGxkbIdLfZ7eOPgeeeA3bvBmRqODPTIwLnzBmgxZFDaLrA7HvL85MACZAACQSOQCADOIG7qv59k6kCcOjQoZAo349//GNEScjKgNaTAGxubtYSROrq6vDMM8/grbfewrZt27S1gT01iRbW1NTAbrfjj3/8I0Swbtq0qccIoCSbLFu2rMvpAiEAT08D17pRkN2A1EGtOH4yHpW1mRicbYEEU0UMDpgWYguGBwxXXggJkAAJkEBAZ/BCFbepAlDW1L3//vvKSSC+8PqaAvYeO3r0aMgaxLvvvls3+yuvvFIb69NPP91tn2BFAL2JIPu32nB5ajnSfTKB7Wm52HS8CKOmWrFkiWdqOOxbCG4cHvZMeQEkQAIkQAKnCTACCJgqAL///e9r6//uuecew8xOrwCU7N6bb74ZErXT2yRhRKZ0n5N5Vh0tUAYkpWB+cZcNU4+XIcPStRZggzsdW1Nn4oc/t2q77YR1C9GioWHNlIMnARIgARLoRCBQ/juUsZsqAGX693e/+x0uuOACnH/++V2SQB577DFdbJqamrSEDmkTJkyA9Js2bRokwihrC6Xm33XXXaet/auvr8eTTz6JF154Qasb6C3pIuJuzpw5WLRokXYeEaVS808E34kTJ7Q6gitWrIAUr77qqqt0jStQBrTzIzc2fLcUFyRUwp7VdTeQ9Loq/Lu1ANc8OQvjzw/z3UBCdNsgXQbBg0iABEiABMKCQKD8dyjDMFUAikjrqUkkT9bp6Wlvv/22Jvj8m2TtrlmzBjfddJO23k/W/4kglOSPe++9t1MSSF5eHhYsWHA6IiilXt58801t3WBaWpomUJcsWaJb/MlYAmVANTvq8dbtJYjKTEdMenIXDk57E1wNdnxxTTFyJ4T5fsAhunG4HjvlMSRAAiRAAuFBIFD+O5RpmCoAQ/nCjRhboAzIdeAQNt6+DjuO5WFEbjR8kqkhM6YHazpQlFGN6WvmISpnuBGXFrxzMAIYPPb8ZhIgARKIEAKB8t+hjJMCUOHuBMyA6uthe7wEpVvTUetIhpQ2jIsDpCJNXR2QndCEWVPssN5ZDGSFeQSQawAVLJJdSYAESIAE9BAImP/WM5ggHUMBqAA+YAZ0ShQd2VyJTYfycfCQBe3tQGwskDPCjcuGVeGsywoAqWHoGx5UuLagdg3BjcODyoNfTgIkQAIkYCiBgPlvQ0dt7MkoABV4BtSATokiWetXF2tFizsBiRYHBrfbEJWVAcyYAVitClcTYl1ZBzDEbgiHQwIkQAIDh0BA/XeIYqMAVLgxATegSBNF3AlEwTrZlQRIgARIoCcCAfffIXgrKAAVbkpQDIiiSOGOsSsJkAAJkAAJBK6KRyizpgBUuDtBEYAK42VXEiABEiABEiABCkCxAQpAhSeBAlABHruSAAmQAAmQQJAI0H9TACqZHg1ICR87kwAJkAAJkEBQCNB/UwAqGR4NSAkfO5MACZAACZBAUAjQf1MAKhkeDUgJHzuTAAmQAAmQQFAI0H9TACoZHg1ICR87kwAJkAAJkEBQCNB/UwAqGR4NSAkfO5MACZAACZBAUAjQf1MAKhkeDUgJHzuTAAmQAAmQQFAI0H9TACoZHg1ICR87kwAJkAAJkEBQCNB/UwAqGR4NSAkfO5MACZAACZBAUAjQf1MAKhkeDUgJHzuTAAmQAAmQQFAI0H9TACoZHg1ICR87kwAJkAAJkEBQCNB/UwAqGR4NSAkfO5MACZAACZBAUAjQf1MAKhkeDUgJHzuTAAmQAAmQQFAI0H9TACoZHg1ICR87kwAJkAAJkEBQCNB/UwAqGR4NSAkfO5MACZAACZBAUAjQf1MAKhkeDUgJHzuTAAmQAAmQQFAI0H9TACoZHg1ICR87kwAJkAAJkEBQCNB/UwAqGR4NSAkfO5MACZAACZBAUAjQf1MAKhkeDUgJHzuTAAmQAAmQQFAI0H9TACoZHg1ICR87kwAJkAAJkEBQCNB/UwAqGR4NSAkfO5MACZAACZBAUAjQf1MAKhkeDUgJHzuTAAmQAAmQQFAI0H9TACoZHg1ICR87kwAJkAAJkEBQCNB/UwAqGR4NSAkfO5MACZAACZBAUAjQf1MAKhkeDUgJHzuTAAmQAAmQQFAI0H9TACoZHg1ICR87kwAJkAAJkEBQCNB/UwAqGR4NSAlfxHd2uYCaGuDECSAlBcjNBaKiIh4LAZAACZCA6QTovykAlYyMBqSEL6I779oFrF8P7N4NtLYC8fHAmDHAnDlAYWFEo+HFkwAJkIDpBOi/w0QAbt68GStXrsT27dths9mwfv16zJ49+7SBLF26FC+//DIOHDiAQYMGYeLEiVi+fDkmTZqky4geffRR3HPPPbjzzjvxy1/+UlcfOYgGpBsVD/QhIOJv9WqgrtaNguwGpA5qxfGT8aiszcTgbAsWL6YIpMGQAAmQgJkE6L/DRACWlZVhy5YtKCoqwty5c7sIwBdffBFDhgzBqFGj4HA4sGrVKqxbtw5VVVXIzs7u1Ybef/99FBcXIzU1FdOmTaMANPOJ47kh074rVgD7t9pweWo50htrENPRBmd0HOxpudh0vAijplqxZAmng2kuJEACJGAWAQrAMBGAvgZgsVi6CEB/A/He2I0bN2L69Ok92k9TU5MmKp988kk8/PDDuPDCCykAzXraeF6NQHU18Iu7bJh6vAwZFjuOJ1vhjElAjNOB1CYbGtzp2Jo6Ez/8uRV5eYRGAiRAAiRgBgEKwAEoAE+ePInVq1drgk4igIMHD+7RdubPn4/MzEwtYnjFFVf0KQDb2tog/7xNDCgnJweNjY1aBJGNBPoisPMjNzZ8txQXJFTCnpUPNyxoaQHa24HYGDeGOarwUWsBrnlyFsafb+nrdPycBEiABEjgDAhQAA4gAbhhwwbceOONaGlpgdVq/f/tnQl4VNXd/78zk2RmspMFmISEAAmBICpBBOF1oUgFWq1IiRSpLFZfV1zqX6RFCoKK9RUt/oso2qJSi4EKfarEvy9SwIJFy6LIEggQsw1Llsk6mSQz9/+cOyQQEmTMmWQmc7/neXh8zNzfufd8zvfO+c5ZsWnTJowYMeKSshBzBsU8QTEEbDKZPDKAYq7h4sWL2+RJA9iBt0+jIQX7yrD1/mzoY6JRqw9HcaGCJls19E2NcAUFw2jUobe5Ej95JwvJw2I1SonFJgESIIHOJUADGEAGsLa2Vl0gUlpaitWrV2Pr1q3YvXu3Ojfw4iQWi1xzzTX49NNPcdVVV6kfswewc1825u4m4Cosxpb712NXcQoUWyV61pxAb8MZhOga4VCCUVgbi9DoYNz07q+Q/qNEYiMBEiABEugEAjSAAWQAL9ZHWloa5syZg/nz57eRjugdnDx5MgwGQ8tnTqcTYn6hXq9Xh3kv/OxS2qOAOuGtDPQsy8pQ8mo23n1XQVzpEcSba1Fn6oF6lxGuegd6K6dgNOtwIutp/Gr5EO4LGOh6YPlIgAR8QoDtdwAbwNTUVMyYMQNi2PbiVF1dje+++67Vn2fPno1BgwZh3rx5uOKKKzwSJAXkESZedCEBRcHpP32M/Ys3oaHeBSsS4FR0MOgBs1nBAFMJdAY9dva8HT//00+Q0o/zACkgEiABEvA2Abbf3cQAitW6YkGHSMOGDcPy5cvVLVvEAo7Y2Fh1Lt9tt92mzv0rKytTV/WuXbtW3TdwyJAhapxYDSx6/R5++OF2deTJEPDFgRSQt19JbeR3dONBWB9bhshwBeXG3nAoRhh1DkQrFWgMCcfp6HScOavDmD9kIeN6zgPUhipYShIgga4kwPa7mxjAbdu2qYbv4iRW8a5atQrTp09X5/uJ+X/CEIrFHwsWLGi1CCQlJQWzZs1qt0dQ5EsD2JWvnrbvVfjvYuz+77cQYWpED2cZDK5GOPXBqAuLR3mPAahwRsFozceN/3cqkkZxHqC21cLSkwAJdAYBGsBuYgA7o/K9kScF5A2K2svDdbYMn/4qG4etUUjqo8DQ1ACD0wGnwYgmQwiKinUYbKnEj9/Kgj6ePYDaUwhLTAIk0NkE2H7TAEppjAKSwqfdYEVBwRubsfvdXJxqiEGq/iTC6s7AWd+IWkcwosIa0W/2jejz1AxAxzmA2hUKS04CJNBZBNh+0wBKaYsCksKn7WCrFade/SuKN+6G9WwwSpp6QlGAhOAzSIhvROIdI9H70V8AFou2ObH0JEACJNAJBNh+0wBKyYoCksKn7WBFQfHv1+LEn7ejsjYIYSFNCDIHozYsHsec/dE7pBwj705H8n9PYi+gtpXC0pPId1QqAAAgAElEQVQACXQCAbbfNIBSsqKApPBpOrjNPEBnI5yGYDSERKjHw505WYPBvW2cB6hplbDwJEACnUWA7TcNoJS2KCApfJoOFiuBtz+8Hg5LCszh5zckb4Zir3FyJbCmFcLCkwAJdCYBtt80gFL6ooCk8Gk6+NDnZdj5aDai+kajyRTehoXBXoOqAhv3AtS0Slh4EiCBziLA9psGUEpbFJAUPk0H559UsOGezUhz5qLWktp6np+iILQkD3lB6fj525N4GoimlcLCk0CAEBCr3MrLgfp6wGQCYmJ8Or+Z7TcNoNSbRQFJ4dN0sMsFvPYbK4K25GBAjA3VERY0BZkR1GRHRLUVx8t7wDl+Ah5+zsLzgDWtFBaeBAKAgNUK7N0LFBQADgdgNALJyUBmps92OmD7TQMo9WZRQFL4NB98+DCw5gUroo7vwRD9YYShFrUIw0HXYFQOGI5Z8y0YPFjzmAiABEigOxMQ5i8nB7DZ3GbPbAbsdkD8PToamDjRJyaQ7TcNoNRrRQFJ4WMwgGM7rPjPm3vg2H8YenstXOYwmIYNxvB7hyPtBu4BSJGQAAl0YwJi2HfzZiA3F0htO9UFeXlAejowqeu3u2L7TQMo9WZRQFL4GGy1ouTtHBzcacOBUguqGs2IDLZjaJwVQ8ZEI+Ee3/wyZsWQAAmQgFcIlJUB2dnunr7wtovdUFPj7hnMygJiu/bYS7bfNIBSGqeApPBpO/iC4+COIxVx8Tp1XrSYH116VsEA5HEjaG0rhKUnge5PoLgYWL8eSEkBDG23u4LTCeTnA1OnAomJXVpett80gFKCo4Ck8Gk6uGUj6FPR6NkvvNWRv2LUhBtBa1oeLDwJBAYB9gD6dT3qFEU0N0wdIUAD2BFqjBEEuBE0dUACJBDwBDgH0K+rmAZQonpoACXgaTyUG0FrXAAsPglohcD3rQLu0QOYMIGrgH2kBRpACfA0gBLwNB7KjaA1LgAWnwS0RID7APplbdMASlQLDaAEPI2HciNojQuAxScBrRHgSSB+V+M0gBJVQgMoAY+haNkI+sRepIUUINTgQJ3TiGMNyajsn8mNoKkREiABEugkAmy/uQpYSloUkBS+ywf72S/Gyz/wD79CmMCNHyoo+roMpspTMIYAva7qjVumx2Jwhu6HZ8gIEiABEiCByxJg+00DeFmRfN8FFJAUvu8P9sM5I51VWlexFac/2YvyfQXQNzoQl2hEbGYy9MN9d05mZ5WV+ZIACZCAPxBg+00DKKVDCkgK36WD/fTsyE4p7bnTQA7tsiG3yoJaxYwwnR3pkVZkjOZpIJ3CnJmSAAlongDbbxpAqZeAApLC136wH+8b5fXS8jQQryNlhiRAAiTgCQG23zSAnujkktdQQFL42g/2453jvV1angbibaLMjwRIwF8JiJ0PCgqA6mogIgJITgb0et89LdtvGkAp9VFAUvjaD/bjsyO9XVqeBuJtosyPBEjAHwmoi902AkeOuM87F+eeDxoETJ4MDB7smydm+00DKKU8CkgKn+Z7AHkaSCfoh1mSAAn4FQFh/lasAEpLgaQkICwMqK0FCguBuDhg7lzfmEC23zSAUi8KBSSFr/1gDc0B5GkgnaAfZkkCJOA3BMSw77JlwL59QEYGoLtgZyvxVX/oEJCZCcyb1/XDwWy/aQClXhQKSArfpYP99OxIb5eWp4F4myjzIwES8CcC+fnAwoXunr7IyLZPVlkJiGnfzz4LpKR07ZOz/aYBlFIcBSSF7/uDNbIP4IWngaQGfQdTfQXs9XqcQAoqrrgRs3+b4LM5Mp1Yu8yaBEhAAwQOHACWLAHS0wGDoW2Bm5qAo0eBZ54Bhg7tWiBsv2kApRRHAUnhu3ywBk4CERCECfzryyVw/O929LDlQwcnmiJiEDe8L37060yk3WC5PCteQQIkQAJ+RoA9gH5WIRc9Ds8ClqgfGkAJeAxtIXBshxW7Fuag4bQNDbEWBEeaYVLsECeE6GOjcd3iiTSB1AsJkEC3I8A5gP5dZTSAEvVDAygBj6EqAZdTQfbszajfn4vgwanqLGmxQq6qCoCiIPJMHiJHpCNrzSToDTwbmLIhARLoXgQuXAXcp8/5VcBFRVwF7OuapAGUqAEaQAl4DFUJFOwrw9b7s6GPiUatPhzH84DTZ4AGhxtQpKEGiWE2jF2ZhRvviCU1EiABEuh2BNrbB1Ds/3f77b7ZAkYAZPvNOYBSLxIFJIXvfLBG5vq1R+vwlmLsfnI97L1S8O1hA8rLFIS5qhCrr0AYalHrNMHcVINjY2Zhwao+XBDiJckxGxIgga4lwJNAupa3J3djD6AnlC5xDQ2gBLzmUI2s9r0UKdED+Nl/Z+NkRTTOFjdgUP1+DFP+g944hWA0oAFGVCAGm+PvRuycyXjkeYtPj0/yQo0zCxIgARLwOQG2392kB3DHjh146aWXsGfPHlitVmzcuBG3i77jc2nRokVYt24dCgsLERISguHDh+O5557DyJEjLykyEbN48eJWn/fq1QunTp3yWJgUkMeo2r/w+/b7i44GJk4ELIG9AlbMAfzz1M0o/eQrJNmPIBN7EY1KKNChAcGqCaxBBE6ahuBE8o0Y9LtfoN9oi8/P0ZSseYaTAAmQgE8JsP3uJgYwJycHO3fuRGZmJqZMmdLGAL7//vvo2bMn+vfvD7vdjldeeQXr169HXl4e4uPj2xWZMIAbNmzAli1bWj43GAyXvL69TCggifdXQyd+XI7S5rdKcPr+Z/Bfzq0Igx1NCIEdRhjhQBCa1F7AU+iNE8bBODzoDhQO/QkGDdb59BzNy5WJn5MACZCAPxNg+91NDOCFItLpdG0M4MUia65YYe7GjRt3SQO4adMm7N+/v8MapYA6jM69/Xt2NiB6+sLD22ZUUwNXuQ1Fo7NQGRSLiAgEbK/Xvz8qRfGUh3BVw39gVA1gsMpDDxcU6BECBxoRhOLgAbBdeQN2jXgMuaWxPj1HU6LmGUoCJEACPifA9jsADWBDQwNWrFiBpUuXqj2AceIMmnaS6AEUw8pRUVEwGo3qcPHzzz+v9iJeKjkcDoh/zUkIKCkpCZWVlYhs75wbn0vcjx+guBhYv149/8elM+DsWcBuB8xmQHTaFhc4cfiTfHxsngqrPhEmEzBoEAKy1+vIhm9hnTUPofYyhCnVcCpBiEC12AcGDoRAB5e6IKQuJBbOvv2x46bfwWrJ9Ok5mn6sLD4aCZAACVyWAA1gABnAjz76CNOmTUNdXR0sFgtE796IESMuKQIxrCyuHThwIE6fPq0axiNHjuDgwYOIjW1/u4325g2KG9AAXvZda3vBuR7Aotpo7Po6HGJPqMZGIDgYam9f3dkamOttKBiVBV1crLo3XmGhb/eN6kApPQop/n/fonDmfOhslTA5bAhRB30b1GFgMRfQgCZEoAYluiQYLPE4MeouHLxiGiqrdD47R9OjgvEiEiABEvBTAjSAAWQAa2tr1QUipaWlWL16NbZu3Yrdu3ercwM9SSJ+wIABeOqpp/DEE0+0G8IeQE9IeniNoqDgjc3Y/W4ujiMVcfE6tZevvh7Yt09BQm0eokem4+yISermyCKJaYOHDgGZmcC8eQiY1bCuM6XY89NnYDh6GKbqUoS7KlCLcDghDs9U1N5AMRcwH/2x3zwaMcP7o+C6O1EdEuuzczQ9rGVeRgIkQAJ+SYAGMIAM4MUKS0tLw5w5czB//nyPxTd+/Hikpqbi9ddf9yiGAvIIU7sXiT2hXvuNFUFbcjAgxobqCAuagsxoqLTjzNdWlLp64ETaBGSMszT7PzWfykr39MFnn1VHjwMjKQqKf78WJSvWw3SqAL1cxbDDDCd0CEUtzKiHHSYUoB8O4EogPg4V0x9Bsa5P4LEIjBplKUiABPycANvvADaAwsjNmDEDYtjWkyR690QP4H333YeFCxd6EsKdxD2i1P5FzYeE9zNZMbBmL6IqCxDkdKC81ogd3yXju5hMlCgWjB3rXifSnJqaEJi9XlYrDvz2ryjK3okBtV+rvX5m1EL0fdYgDKdhQa0uAmalDi69Ad8On4Uvek5G/zGWgOoNlZAUQ0mABEjAYwI0gN3EANbU1KgLOkQaNmwYli9fjrFjxyImJkadryf2/LvtttvUuX9lZWVYuXIl1q5dq+4bOGTIEDVOrAaePHkyHn74YfX/n3zySdx6661ITk7GmTNn1DmA27dvx4EDB9C3b1+PREQBeYSp3YsOHACWLAHS0wGDXoG5rgzhNadQXQ1s3tsbjvBY1NTqcMMNQK9e57MIyB7Ac8Ur/NKKtx/8DwYd+huGOvZAcblwFmIbIwWNCFaHhSNQBT0UnAnvD0fmGIxeMhFpNwT2XokdVxkjSYAESKB9Amy/u4kB3LZtm2r4Lk4zZ87EqlWrMH36dHW+n5j/JwyhWPyxYMGCVotAUlJSMGvWrJYeQbFgRGwwLWLEXoGjRo3CkiVLkJGR4fH7QgF5jKrNhc09gGKRdoLOCovV3QtoaHIgN9+Ib6uTcSIqE1feYmnpAQzUOYDNcMSw+LynFOSu2YV7bMvR6BQLQJyIRRnqYIYR9ahGJErQBxHhCkZOSUTq1Exg0vl5kh2vEUaSAAmQgHYIsP3uJgbQXyVJAXW8ZoTZWbYMOLnLikn6HIQ6bKgKd88DrCu349ReK+ymaOgnTYSrl0VdBSxWCgvDOHeu7w4Q73iJPYs8eBB45I5iXH/0LSSgGGOwU10VLJKYF2jVJeBkUBriw+wYef8wJCXrgaws4BIr1z27K68iARIgAW0RYPtNAyileApICh8OH1Lw2ZObEV6Ui6aUVBhNOohtFktLAb1OQZouD8eD0/Fl7CSYzDoMHgyIEwDFfwM5/WFhGc48txqDXGLJxwE0IQjBuiZ1Y+gQXRMagsJw0tAfaU/cikFx5cDUqUBiYiAjYdlIgARIwKsE2H7TAEoJigKSwqcu57X+IRv/PhKNk2fDW/YBTEoCRo0CEqNqUHHShtM3ZiE0KTZgTwK5mOLX+1zYNm4xrrbtQKxyFmIlSJ0+HIrOAGOQE/GNJSgN6o2Yx2eyB1BSggwnARLQJgG23zSAUsqngKTwAedOA3Elp+BsuaHVSSB6PQCnExCTBTXWw1Wwrwz/nP4G+hVsR0L9CXVb6CpEQQ8njEo9FJ0OJrMeiTcPRsSvpgE/+UnLXomSNcJwEiABEtAEAbbfNIBSQqeApPB5dB4wbDbNzXFzFRZj+5w1sO8/gpDKMvRuLECIrgH1MKMuKAI1ShgSgkuRMKYfol9dDJxb6S5ZGwwnARIgAc0QYPtNAygldgpICp/7aI/Nm4HcXLj6p+Jsqe58L2CcAv2JPPc+MVpb5VpWhlPPvok9aw8hz5EEi7EM8U2nYHLWwunSqdvm6I3BcF5xJUatnQt9Euf/SSqR4SRAAhojwPabBlBK8hSQFD53sNWKkrdzcGiXDblVFtQqZoTp7EiPtCJjTA8kzJkAWDS2z52i4Mxr63DshQ3I16Wg1hkKp1NBqM6OMGMjkkPL1BXTRUH98aM37kTysPbPrvZC7TALEiABEghIAmy/aQClhE0BSeFTgw8fBta8YEXUib0YaCyA2eCA3WnEUUcyKvtnYtZ8S8Cv+m2P4uHPSnD8/hcxQJePkrBUOGCCUedAtFKBxuAw1Bh74CvlWvx05SQMvdJ9VjITCZAACZCAZwTYftMAeqaUS1xFAUnhQ/NegPv2ARmDFYTWlyOoqR5NQSbUmWJw6LAOmZnQ5FFnYu3L6of2Y3zJGsTXF6IxOBQN4p8xEo0hobCGpOCLyAn49f9YAudMZDk5MZoESIAEPCbA9psG0GOxtHchBSSFT13gK45dFps7R0a2zSuQj327HLmWjbJ3lmC8cTtiKvMBxYl6cwxsUX2xvSqT5wBfDiI/JwESIAF24FxSAzpFETPxmTpCgAawI9TOx7Q6D9jQNq+mJuDoUeCZZ4ChQ+Xu1R2jxfD4ihVA6VkFA+PKERlSj6oGE46WxiAuXhfQJ6J0x/riM5MACXQfAmy/2QMopVYKSAofewA9wCdM4MaNwJEjQH09YDK5T0LRwokoHuDhJSRAAiTQIQJsv2kAOySc5iAKSApf6zmAGa33Mhb90ocOQbNzAC8kK4aDCwqA6mogIgKaORFFTl2MJgESIIFLE2D7TQMo9X5QQFL41OCWYc5SoE8fICwMqK0FiorccwPnzg38s3/lKTIHEiABEiCBH0KA7TcN4A/RS5trKSApfC3BHOb0DkfmQgIkQAIk4BkBtt80gJ4p5RJXUUBS+FoFc5jTeyyZEwmQAAmQwPcTYPtNAyj1jlBAUvgYTAIkQAIkQAI+IcD2mwZQSngUkBS+tsFi5Ud5+fnlrjExrVeGePl2zI4ESIAESECbBNh+0wBKKZ8CksLXOthqBfbudS93dTgAo9G93FUcBaK1s4C9iJVZkQAJkAAJtCXA9psGUOq9oICk8J0PFuYvJweuchtKgy2oU8wI1dkR12iFPiYamDiRJtBLqJkNCZAACZAAwPabBlDqPaCApPC5g8Ww7+bNOL0jF9uLU1FUrENjIxAcDPRJVHBjYh563ZAOTJrE4WAv4GYWJEACJEACNIBCAzwKTuJNoAGUgNccWlYG6x+ysXlXNM7aw9W9/8RpF+LUi9JSIN5cg0mjbbA8mgXExnrhhsyCBEiABEhA6wTYftMASr0DFJAUPjXYVViMLfevx76KFPRJNkCnO5+n6BwsKnAis0c+xq2aCn1SovwNmQMJkAAJkIDmCbD9pgGUegkoICl8anDBvjJsvT9bnesXFB3eJsMmW406N/BHq7KQPIw9gPLEmQMJkAAJkADbbxpAqbeAApLCpwYf+EbBRw9uxlXmXNhiU9scCBxdmoev69Px05WTMPTKC7oH5W/NHEiABEiABDRKgO03DaCU9CkgKXxqcH4+8PKTVoypykEPnQ1V4RY0BZkR1GRHZI0VZUoPfBE5Ab/+HwtSUuTvxxxIgARIgARIgO03DaDUW0ABSeFTg8URcMuWASd3WXFj5F5EVxYgyOlAk8EIW1Qytldlov8YC+bNA/R6+fsxBxIgARIgARJg+00DKPUWUEBS+FqCDx8GVqwASs8qGBhXjsiQelQ5jCgpAXpFOzDnQRMGjuKpIN6hzVxIgARIgATYftMASr0FFJAUvlbBwgRu3AgcOQKYKqwYWLsXV0QU4OoMB3on81QQ75FmTiRAAiRAAmy/aQCl3gIKSApfm2AxHFz8Hyt0n+TA7LChR4YF+jAzYLcD4rSQaJ4K4l3izI0ESIAEtEmA7TcNoJTyKSApfG2DFQWujzej4otcVPdKhTlUh/j4c3P/xKaAeXlAOk8F8TJ1ZkcCJEACmiPA9psGUEr0FJAUvjbBR78ow4ll2TheFo0ahLuPg+sDXHcdkJQEoKYGsNmALJ4K4l3yzI0ESIAEtEWA7TcNoJTiKSApfK2CxRzAd54vRvo36+FKToEx1NByHFxkJHDrrUBSgtO9b8zUqUAiTwXxHn3mRAIkQALaIsD2mwZQSvEUkBS+luDmrWCO/bsMdzizUW+KRmOI+1QQMfJbWAgMGABMuaUG+ir2AHqHOnMhARIgAe0SYPvdTQzgjh078NJLL2HPnj2wWq3YuHEjbr/99hblLlq0COvWrUNhYSFCQkIwfPhwPPfccxg5cuT3qru4uBjz5s1DTk4O7HY7Bg4ciLfffluN9yRRQJ5Quvw1olNv4UIgLlbB8NObEVeWi7Ie508FqasDqioVzByTh9jRnAN4eaK8ggRIgARI4PsIsP3uJgZQGLSdO3ciMzMTU6ZMaWMA33//ffTs2RP9+/dXjdwrr7yC9evXIy8vD/FiFUE7qaKiAsOGDcPYsWPxwAMPqPHHjx9HSkoKBojuJg8SBeQBJA8uOXAAWLLEvb4jqs6K1LwcmOvPnwqid9hRf9KKH03pgaR7JwAWiwe58hISIAESIAESaJ8A2+9uYgAvrD6dTtfGAF5cvc0Vu2XLFowbN67d2n/66adVU/n55593+P2ggDqMrlVgSw9gHCDm+4VXW9Hbuhch1gLA4UC9YsQZUzLufDETSdfS/HmHOnMhARIgAe0SYPsdgAawoaEBK1aswNKlS9UewLi4uHYVnpGRgVtuuQVFRUXYvn07EhMT8eCDD+Lee+/1+I2ggDxG9b0XNs8B3LcPyMgASkuBI4cVNJ4uh76hHhV2E3pnxGD5KzoMGeKdezIXEiABEiAB7RJg+x1ABvCjjz7CtGnTUFdXB4vFgk2bNmHEiBGXVLfJZFI/e+KJJzB16lR8+eWXeOyxx/DGG2/g7rvvbjfO4XBA/GtOQkBJSUmorKxEpOi6Yuowgebj4E6ccO/53NgIiCqqr4e6HYwY9e3fH5g7Fxg8uMO3YSAJkAAJkAAJgAYwgAxgbW2tukCktLQUq1evxtatW7F79251bl97SSwWueaaa7Br166Wj+fOnYuvvvoKX3zxRbsxYrHJ4sWL23xGA+idb5ODB4HHHwdyc4GwMLfxEx24gwa5/3voEJCZCcybd25zaO/clrmQAAmQAAlojAANYAAZwIu1m5aWhjlz5mD+/Pntyrpv374YP3483nrrrZbPX3/9dXXoWKwObi+xB7BzvyHEXMBnFijooZQjVF8PfagJwb1ioNPr1BtXVgJlZcCzzwIpKZ37LMydBEiABEggcAnQAAawAUxNTcWMGTMgeu3aS9OnT1e3jblwEcjjjz+u9hpe2Cv4ffKngLz75XB4qxWbFu7FkIgChCgONBmMqIxKhtWSiZoIC5qagKNHgWeeAYYO9e69mRsJkAAJkIB2CLD97iYGsKamRl3QIZLYumX58uXq9i0xMTGIjY1V9/y77bbb1Ll/ZWVlWLlyJdauXavuGzjk3KoBsRp48uTJePjhh9V8xFDv6NGj1SHdrKwsdQ6gWADy5ptv4q677vLoLaCAPMLk2UVWK86+l4P/zbahKd6C4EgzgprsiKyxwm6KRl7qRBS7RP2yB9AzoLyKBEiABEjgUgTYfncTA7ht2zbV8F2cZs6ciVWrVkH05omeOzH/TxhCsfhjwYIFrRaBiP39Zs2a1apHUCwcEUPEx44dQ79+/dQFIVwF7IMvDHHcx+bNcB3OxYb9qThxUqee/asTI7+KgtiKPJyNTcffGyYhc7iOcwB9UEW8JQmQAAkEEgEawG5iAP1VdBSQl2pGdOtlZwPR0SisCMc//gFUVbkXfhiNgFJdg/pTNhy9Ogv3PBXLVcBews5sSIAESECrBNh+0wBKaZ8CksJ3Plgsulm/3r2yw2BQz/4VC7GLitzbwYQYnLgiPB8p/2cq0m5K9NJNmQ0JkAAJkIBWCbD9pgGU0j4FJIXvfPAFPYAID1f/LjaHPnsWsNuBUFcN4oJs0E/LAmJjvXRTZkMCJEACJKBVAmy/aQCltE8BSeE7H3xuDqC6AWBq6rnJf+c+Fp+JBUDioOBJk1p/5qXbMxsSIAESIAFtEWD7TQMopXgKSApf62Bx/EdODmCzuY/9MJvd3X/i7z16ABMmuP/ORAIkQAIkQAKSBNh+0wBKSYgCksLXNliYvb17gYICQBy5J1aAJCe7j/+g+fMybGZHAiRAAtolwPabBlBK/RSQFL52g11OBUXflKO2rB5hsSb0uTIGeoP7JBAmEiABEiABEvAGAbbfNIBSOqKApPC1CT58GNi4ERD/LS93n/crzgGeNQs4t5+3d2/I3EiABEiABDRJgO03DaCU8CkgKXytgoXpW7ECOHECqKkBqquB+nr3SHDv3u7j3376U+/djzmRAAmQAAlolwDbbxpAKfVTQFL4WoLFli/LlgHbt0M96q2uDoiKAoKDgYYG936AiYnAn/7EnkDvEGcuJEACJKBtAmy/aQCl3gAKSApfS3B+vruHT/T+VVQA8fGtd3sRi4FPnQLuvNNtFMXQMBMJkAAJkAAJdJQA228awI5qR42jgKTwtQQfOAA8/TTUE0DCwgBjiILIpnKEuOrRoDfBpo/BqdM6dSvA115zHxjCRAIkQAIkQAIdJcD2mwawo9qhAZQi1zpY9AA+8gjUfaCH9CjB1VXbYanPhw5OVAfFoFDfF/t1mQgdYFF7AIcO9eLNmRUJkAAJkIDmCNAA0gBKiZ4CksLXEizmAM6bBxxYux8zlTWwNBWiXh8KhyEUNUGRKLeHoiq6L6xXTcSTL1vYA+gd7MyFBEiABDRLgO03DaCU+CkgKXytgnP/WYI9015EbGU+zkSlQgkxQd/kgMlegYagUDRFxcD4X9di6p8mcV9A72FnTiRAAiSgSQJsv2kApYRPAUnhOx8szvtdtw4lr23AlhMpqKgLRVAQYDCIw0AUWBQrdHExuPJn/WB59E4gNtZLN2Y2JEACJEACWiTA9psGUEr3FJAUvvPBYu+XN98EDh3CGXMyvj6gV7eDUVzuI4F7R9sx0FKNmCuTgNmz3XvCMJEACZAACZBABwmw/aYB7KB03GEUkBS+88HFxcCaNe4zgCMi4DKaUVkFNDiAECMQFe6CvqgAyMgA7ruPPYBews5sSIAESECrBNh+0wBKaZ8CksLXugfwgw8AsRxYdP1ZLK03AhQ7Q4vPfv5zYNq01p956RGYDQmQAAmQgHYIsP2mAZRSOwUkhe98sJgDuHkz8NVX548C6dFDTAB0nweXlwf06wc89RSQkOClmzIbEiABEiABrRJg+00DKKV9CkgKX+tgqxXIyXEPA4sev8pK93/Fv6QkYNYs4OqrvXhDZkUCJEACJKBVAmy/aQCltE8BSeFrGyxM4N69wHffuc+EE2e+iWM/bryRPX9eRs3sSIAESEDLBNh+0wBK6Z8CksLXfrAYDi4vdw/9mkxATAzn/HUCZmZJAiRAAlomwPabBlBK/xSQFD4GkwAJkAAJkIBPCLD9pgGUEh4FJIWPwSRAAiRAAiTgEwJsv2kApYRHAUnhYzAJkAAJkAAJ+OCTmjQAAAzbSURBVIQA228aQCnhUUBS+BhMAiRAAiRAAj4hwPabBlBKeBSQFD4GkwAJkAAJkIBPCLD9pgGUEh4FJIWPwSRAAiRAAiTgEwJsv2kApYRHAUnhYzAJkAAJkAAJ+IQA228aQCnhUUBS+BhMAiRAAiRAAj4hwPabBlBKeBSQFD4GkwAJkAAJkIBPCLD9pgGUEh4FJIWPwSRAAiRAAiTgEwJsv2kApYRXWVmJ6OhoFBYWIjIyUiovBpMACZAACZAACXQNAWEAk5KSYLPZEBUV1TU39bO76BRFHL7K1BECRUVFqoCYSIAESIAESIAEuh8B0YHTp0+f7vfgXnhiGkAJiC6XCyUlJYiIiIBOp5PIqW1o868TLfYusuxJmu1VZt2z7vmdp63RJF+986Lvq7q6GgkJCdDr9V5tv7tLZjSAflpTWp6fwLJHQUwv0OK0AtY9616L2qfutat7X1oQGkBf0v+ee/MLQZtfCFqud/E6aLn8Wi47675KnYdG86ut3k9f2w8aQF/XwCXur+XGgGXXZkNAE6BdE8C6127da/n73tf2gwbQ1zVwifs7HA688MILmD9/PoxGo58+Zec8FsuuzXoXamLds+75ncfv+85pWZjrxQRoAKkJEiABEiABEiABEtAYARpAjVU4i0sCJEACJEACJEACNIDUAAmQAAmQAAmQAAlojAANoMYqnMUlARIgARIgARIgARpAaoAESIAESIAESIAENEaABtAPK3zlypV46aWXYLVaMWTIELz66qu4/vrr/fBJvftIO3bsUMu9Z88etewbN27E7bff7t2b+GluYsX3hx9+iCNHjsBsNmP06NF48cUXkZ6e7qdP7L3Hev311yH+5efnq5kKzS9cuBATJ0703k26SU5CB7/5zW/w6KOPqu99oKdFixZh8eLFrYrZq1cvnDp1KtCL3lK+4uJizJs3Dzk5ObDb7Rg4cCDefvttDB8+PKAZpKSk4LvvvmtTxgcffBB//OMfA7rs/lI4GkB/qYlzz/HBBx/gl7/8JYQJHDNmDN544w289dZbOHToEJKTk/3sab37OOILcOfOncjMzMSUKVM0ZQAnTJiAadOmYcSIEWhqasJvf/tbHDhwQK33sLAw74L2s9z+8Y9/wGAwIDU1VX2yd955R/0hsG/fPtUMaiV99dVXyMrKUk+AGTt2rGYM4IYNG7Bly5aWahZaiI+P10S1V1RUYNiwYWp9P/DAA+jZsyeOHz8OYY4GDBgQ0AzOnj0Lp9PZUsZvv/0W48ePxz//+U/cdNNNAV12fykcDaC/1MS55xg5cqRqgESPSHMaPHiw2hMmege0ksTZylrqAby4XsWXo2gMtm/fjhtuuEEr1d5SzpiYGNUE3nPPPZooe01Njfreix9+S5cuxdVXX60ZA7hp0ybs379fE/V8cSGffvpp9Ufv559/rsnyX1joxx57DB999BGOHTsG8f3P1PkEaAA7n7HHd2hoaEBoaCjWr1+PyZMnt8SJ4SDxBSnMgFaS1g1gXl4e0tLS1F7AK664QivVrvYICP3PnDlT7QHMyMjQRNlFeYXpfeWVV9TeDy0ZQGH0xTFoYsN78QP4+eefR//+/TVR70Lft9xyC4qKitTv98TERIgh0HvvvVcT5W8upGj7EhIS8MQTT6hTIJi6hgANYNdw9uguJSUl6heA+EUo5oA1J/GFKIbFcnNzPconEC7SsgFUFAU/+9nPIIaHtNIzIIzuddddh/r6eoSHh+P999/HpEmTAkHKly3DunXr8Nxzz0EMAZtMJk0ZQDHto66uTp33dvr0abX3U8yDPXjwIGJjYy/LrrtfIOpbJGF8pk6dii+//BKiJ0xM/bn77ru7e/E8fv7s7GxMnz4dBQUFqhFk6hoCNIBdw9mjuzQbwF27dqmNYXMSjcN7772nfjFqJWnZAD700EP4+OOP8a9//Qt9+vTRRJWLHgDx5W+z2fC3v/1NnfcqekQCvQewsLAQ11xzDT799FNcddVVal1rqQfwYnHX1taqc9+eeuop1RQFegoJCVHrX3znN6e5c+eqPwa++OKLQC9+S/lEL6hgIeYDM3UdARrArmN92TtxCPg8Iq0awEceeQRiTpRYEd2vX7/LaiZQL7j55ptVIyB6QgI5iboW0z3EwofmJIbBhf71er16NvKFnwUyi+ayiYUAYkHQhfOgA7Xcffv2VRc+iB88zUmUW/SEitXBWkhiJbAY8he7IIiRD6auI0AD2HWsPbqTmAMjlv+LyeDNSfSCiBeDi0A8QtgtLxLDvsL8iYUv27ZtU+f/aTmNGzcOSUlJWLNmTUBjqK6ubrMVxuzZszFo0CB1axAtzf8UFS0MrzD+9913n7oVUKAnMewpeoEvnOrx+OOPY/fu3a16BQOZg9gKSPzQExyCgoICuah+VzYaQD+rkuZtYFatWqUOA7/55ptYvXq1OidG/FoM5CRWQorFDyKJrRGWL1+ubo8gJscH+hY4YuK3mPf297//vdXef2JyvNgXMJCTmPQt9vwThk8YIjEnbtmyZfjkk0/U3hGtJS0NAT/55JO49dZb1ff7zJkzas+XGPoXc0ID/ftO6FoM9Yr53mIvRLEFkJgDKBaAiO/9u+66K+Cl73K51JGOX/ziF+o7z9S1BGgAu5a3R3cTvX+///3v1c2QRQ+AWBmoha1ARM+XMHwXJ7FCMtB7gi617cGf//xnzJo1yyPddNeLxFYvn332map3YXivvPJKtfdLi+ZP1KGWDKDY+1JMdygtLVX3/hs1ahSWLFkS8HM/L3xXxdYn8+fPV7c/EWZIzH3UyipgMfdVzP8TCxzFQiCmriVAA9i1vHk3EiABEiABEiABEvA5ARpAn1cBH4AESIAESIAESIAEupYADWDX8ubdSIAESIAESIAESMDnBGgAfV4FfAASIAESIAESIAES6FoCNIBdy5t3IwESIAESIAESIAGfE6AB9HkV8AFIgARIgARIgARIoGsJ0AB2LW/ejQRIgARIgARIgAR8ToAG0OdVwAcgARIgARIgARIgga4lQAPYtbx5NxIgAT8g0LzpeEVFBaKjo/3gifgIJEACJNC1BGgAu5Y370YCJOADAhefrtHQ0IDy8nL06tULlzqFxQePyVuSAAmQQJcRoAHsMtS8EQmQgK8IaOl4NV8x5n1JgAS6FwEawO5VX3xaEiCBH0hAnKX8zjvvtIoSZyzPnj0bzUPA4qzpxx57DGvXrsWvf/1rFBYWYtKkSWrchg0b8Lvf/Q6VlZWYMWMGXn31VRgMBjU/0ZO4YMEC/OUvf4HNZlPP7n7xxRfV83yZSIAESMCfCdAA+nPt8NlIgASkCQjjNnHiRNWcPfvss2p+Bw8exM0339zKAN53332qcRMGrrq6GnfccQeGDx+uzhEUBvDEiROYMmUK3n33Xdx5551qPnfddRfy8/OxbNkyJCQkYOPGjaohPHDgANLS0qSfnRmQAAmQQGcRoAHsLLLMlwRIwG8IXDwEfPEiENEDKHoE8/LyMGDAAPW577//frz33ns4ffo0wsPD1b9NmDABKSkpWLVqFY4fP66avKKiItX8NSdhLK+99lo8//zzflN+PggJkAAJXEyABpCaIAESCHgCnhjAhx56CLW1tS0sRK+fGP4VvYXNaebMmWrv4Icffoj169cjKysLYWFhrfg5HA619/CDDz4IeK4sIAmQQPclQAPYfeuOT04CJOAhAU8MoJgDKObxNadFixZh06ZN2L9/f8vfxHxCcY34uzB4YghYGMTmOYHNF4oew969e3v4dLyMBEiABLqeAA1g1zPnHUmABLqYwI9//GOkp6fjtddeU+/c3hDwDzWAR48eVfPcsWMHrr/++i4uEW9HAiRAAnIEaADl+DGaBEigGxAQCzxET152drY6n++bb77BuHHj2qwC/iE9gKLYYlXwzp078fLLL2PYsGEoLS3F1q1bMXToUHUVMRMJkAAJ+CsBGkB/rRk+FwmQgNcIiN46MX/v66+/ht1ux6W2gfmhBrCxsRFLly5VVwYXFxcjNjYW1113HRYvXqyaQCYSIAES8FcCNID+WjN8LhIgARIgARIgARLoJAI0gJ0EltmSAAmQAAmQAAmQgL8SoAH015rhc5EACZAACZAACZBAJxGgAewksMyWBEiABEiABEiABPyVAA2gv9YMn4sESIAESIAESIAEOokADWAngWW2JEACJEACJEACJOCvBGgA/bVm+FwkQAIkQAIkQAIk0EkEaAA7CSyzJQESIAESIAESIAF/JUAD6K81w+ciARIgARIgARIggU4iQAPYSWCZLQmQAAmQAAmQAAn4KwEaQH+tGT4XCZAACZAACZAACXQSARrATgLLbEmABEiABEiABEjAXwnQAPprzfC5SIAESIAESIAESKCTCPx/OsFgxnvJo7QAAAAASUVORK5CYII=\" width=\"640\">"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.figure(1)\n",
    "plt.plot(t0, mag, 'bo', alpha=0.5)\n",
    "plt.plot(t0, func(t0, *popt), 'ro', alpha = 0.3)\n",
    "plt.title('KIC_6763132_PTF_R_4505_05_noshift')\n",
    "plt.xlabel('time')\n",
    "plt.ylabel('mag')\n",
    "plt.gca().invert_yaxis()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.6.8"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}