transtivity.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "127.4(8.32)\n"
     ]
    }
   ],
   "source": [
    "import math\n",
    "def avg(e1,s1,e2,s2):\n",
    "    e=(e1*s2*s2+e2*s1*s1)/(s1*s1+s2*s2)\n",
    "    s=s1*s2/math.sqrt(s1*s1+s2*s2)\n",
    "    print(\"%.1f(%.2f)\"%(e,s))\n",
    "\n",
    "avg(125.3,10,132,15)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "  R&    0&  194&  275&  244&  319\\\\\n",
      "  I& -194&    0&   80&   41&  127\\\\\n",
      "  G& -275&  -80&    0&    3&   60\\\\\n",
      " RT& -244&  -41&   -3&    0&   99\\\\\n",
      "ZTS& -319& -127&  -60&  -99&    0\\\\\n",
      "{'R': {'I': (194, 4), 'RT': (244, 8), 'G': (275, 5), 'ZTS': (319, 16), 'R': (0, 0)}, 'I': {'RT': (41, 8), 'G': (80, 6), 'ZTS': (127, 8), 'R': (-194, 4), 'I': (0, 0)}, 'RT': {'G': (-3, 9), 'ZTS': (99, 18), 'R': (-244, 8), 'I': (-41, 8), 'RT': (0, 0)}, 'G': {'ZTS': (60, 5), 'R': (-275, 5), 'I': (-80, 6), 'G': (0, 0), 'RT': (3, 9)}, 'ZTS': {'R': (-319, 16), 'I': (-127, 8), 'G': (-60, 5), 'RT': (-99, 18), 'ZTS': (0, 0)}}\n"
     ]
    }
   ],
   "source": [
    "#AIs=[\"R\",\"I\",\"RT\",\"G\",\"ZTS\"]\n",
    "AIs=[\"R\",\"I\",\"G\",\"RT\",\"ZTS\"]\n",
    "def antisym(d):\n",
    "    for i in AIs:\n",
    "        for j in AIs:\n",
    "            if i==j:\n",
    "                d[j][i]=(0,0)\n",
    "            if i not in d[j]:\n",
    "                d[j][i]=(-1*d[i][j][0],d[i][j][1])\n",
    "    return d\n",
    "\n",
    "def totable(d):\n",
    "    for i in AIs:\n",
    "        print(\"%3s& \"%(i)+\"& \".join([\"%4d\"%(d[i][j][0],) for j in AIs])+\"\\\\\\\\\")\n",
    "\n",
    "d={\"R\":{\"I\":(194,4),\"RT\":(244,8),\"G\":(275,5),\"ZTS\":(319,16)},\n",
    "   \"I\":{\"RT\":(41,8),\"G\":(80,6),\"ZTS\":(127,8)},\n",
    "   \"RT\":{\"G\":(-3,9),\"ZTS\":(99,18)},\n",
    "   \"G\":{\"ZTS\":(60,5)},\n",
    "   \"ZTS\":{}}\n",
    "d=antisym(d)\n",
    "totable(d)\n",
    "print(d)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 36,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[('R,I,G', 1.0, 162.20894646124296, 0.006164888076866542), ('R,I,ZTS', 2.0, 198.82243025498838, 0.010059227208092237), ('R,G,ZTS', 16.0, 173.95017730614614, 0.09198036039848678), ('I,G,ZTS', 13.0, 84.7907192881195, 0.1533186663486826)]\n",
      "0.09673059595964928\n"
     ]
    }
   ],
   "source": [
    "import itertools,random,numpy\n",
    "d={'R': {'I': (194, 4), 'RT': (244, 8), 'G': (275, 5), 'ZTS': (319, 16), 'R': (0, 0)}, \n",
    "   'I': {'RT': (41, 8), 'G': (80, 6), 'ZTS': (127, 8), 'R': (-194, 4), 'I': (0, 0)}, \n",
    "   'RT': {'G': (-3, 9), 'ZTS': (99, 18), 'R': (-244, 8), 'I': (-41, 8), 'RT': (0, 0)}, \n",
    "   'G': {'ZTS': (60, 5), 'R': (-275, 5), 'I': (-80, 6), 'G': (0, 0), 'RT': (3, 9)}, \n",
    "   'ZTS': {'R': (-319, 16), 'I': (-127, 8), 'G': (-60, 5), 'RT': (-99, 18), 'ZTS': (0, 0)}}\n",
    "AIs=[\"R\",\"I\",\"G\",\"ZTS\"]\n",
    "def get_transitivity(gamma=0.0):\n",
    "    l=[]\n",
    "    for i,j,k in itertools.combinations(AIs,3):\n",
    "        x=random.gauss(d[i][j][0],d[i][j][1]*gamma)\n",
    "        y=random.gauss(d[j][k][0],d[j][k][1]*gamma)\n",
    "        z=random.gauss(d[i][k][0],d[i][k][1]*gamma)\n",
    "        n=abs(x+y-z)\n",
    "        #m=abs(x)+abs(y)+abs(z)\n",
    "        m=(abs(x)*abs(y)*abs(z))**(1/3)\n",
    "        name=\",\".join([i,j,k])\n",
    "        l.append((name,n,m,n/m))\n",
    "        #print(\"%8s %2d %3d %.3f\"%(name,n,m,n/m))\n",
    "    print(l)\n",
    "    #return sum([i[3]*i[1] for i in l])/sum([i[1] for i in l])\n",
    "    return sum([i[3]*i[1]*i[2] for i in l])/sum([i[1]*i[2] for i in l])\n",
    "\n",
    "print(get_transitivity())\n",
    "\n",
    "#l=[get_transitivity() for i in range(10000)]\n",
    "#print(numpy.mean(l))\n",
    "#print(numpy.sqrt(numpy.var(l)))\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 40,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0.07167746756153773"
      ]
     },
     "execution_count": 40,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "\"\"\"trans_data=[('R,I,G', 1.0, 549.0, 0.0018214936247723133), \n",
    "            ('R,I,ZTS', 2.0, 640.0, 0.003125), \n",
    "            ('R,G,ZTS', 16.0, 654.0, 0.024464831804281346), \n",
    "            ('I,G,ZTS', 13.0, 267.0, 0.04868913857677903),\n",
    "            ('origin',0,0,0)]\"\"\"\n",
    "\"\"\"trans_data=[('R,I,G', 1.0, 549.0, 0.0018214936247723133), \n",
    "            ('R,I,RT', 9.0, 479.0, 0.018789144050104383), \n",
    "            ('R,I,ZTS', 2.0, 640.0, 0.003125), \n",
    "            ('R,G,RT', 34.0, 522.0, 0.06513409961685823), \n",
    "            ('R,G,ZTS', 16.0, 654.0, 0.024464831804281346), \n",
    "            ('R,RT,ZTS', 24.0, 662.0, 0.03625377643504532), \n",
    "            ('I,G,RT', 42.0, 124.0, 0.3387096774193548), \n",
    "            ('I,G,ZTS', 13.0, 267.0, 0.04868913857677903), \n",
    "            ('I,RT,ZTS', 13.0, 267.0, 0.04868913857677903), \n",
    "            ('G,RT,ZTS', 42.0, 162.0, 0.25925925925925924),\n",
    "            ('origin',0,0,0)]\"\"\"\n",
    "trans_data=[('R,I,G', 1.0, 162.20894646124296, 0.006164888076866542),\n",
    "            ('R,I,ZTS', 2.0, 198.82243025498838, 0.010059227208092237),\n",
    "            ('R,G,ZTS', 16.0, 173.95017730614614, 0.09198036039848678),\n",
    "            ('I,G,ZTS', 13.0, 84.7907192881195, 0.1533186663486826),\n",
    "            ('origin',0,0,0)]\n",
    "trans_data.sort(key=lambda x:x[2])\n",
    "\n",
    "def int_area(trans_data):\n",
    "    s=0\n",
    "    for i in range(len(trans_data)-1):\n",
    "        yi=trans_data[i][1]\n",
    "        yj=trans_data[i+1][1]\n",
    "        xi=trans_data[i][2]\n",
    "        xj=trans_data[i+1][2]\n",
    "        s+=(yi+yj)*(xj-xi)\n",
    "    return s\n",
    "        \n",
    "#int_area([(None,0,0),(None,1,1),(None,1,2)])\n",
    "int_area(trans_data)/max([i[2] for i in trans_data])**2"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 39,
   "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 overriden (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,iVBORw0KGgoAAAANSUhEUgAAAlgAAAHCCAYAAAAzc7dkAAAgAElEQVR4XuydB5hU1fmHX9pSlipFAQELdlEiolixa1QUFaKxgsFu8jeWKBqjSVSiRqMmxoKIPSqI2I29KypWrIjSQZCysAts5f98e+/AuNlyZ+fO3PY7z7OPuHPuOd/3fmdnfnPKd5qgIgIiIAIiIAIiIAIi4CuBJr62psZEQAREQAREQAREQASQwNIgEAEREAEREAEREAGfCUhg+QxUzYmACIiACIiACIiABJbGgAiIgAiIgAiIgAj4TEACy2egak4EREAEREAEREAEJLA0BkRABERABERABETAZwISWD4DVXMiIAIiIAIiIAIiIIGlMSACIiACIiACIiACPhOQwPIZqJoTAREQAREQAREQAQksjQEREAEREAEREAER8JmABJbPQNWcCIiACIiACIiACEhgaQyIgAiIgAiIgAiIgM8EJLB8BqrmREAEREAEREAEREACS2NABERABERABERABHwmIIHlM1A1JwIiIAIiIAIiIAISWBoDIiACIiACIiACIuAzAQksn4GqOREQAREQAREQARGQwNIYEAEREAEREAEREAGfCUhg+QxUzYmACIiACIiACIiABJbGgAiIgAiIgAiIgAj4TEACy2egak4EREAEREAEREAEJLA0BkRABERABERABETAZwISWD4DVXMiIAIiIAIiIAIiIIGlMSACIiACIiACIiACPhOQwPIZqJoTAREQAREQAREQAQksjQEREAEREAEREAER8JmABJbPQNWcCIiACIiACIiACEhgaQyIgAiIgAiIgAiIgM8EJLB8BqrmREAEREAEREAEREACS2NABERABERABERABHwmIIHlM1A1JwIiIAIiIAIiIAISWBoDIiACIiACIiACIuAzAQksn4GqOREQAREQAREQARGQwNIYEAEREAEREAEREAGfCUhg+QxUzYmACIiACIiACIiABJbGgAiIgAiIgAiIgAj4TEACy2egak4EREAEREAEREAEJLA0BkRABERABERABETAZwISWD4DVXMiIAIiIAIiIAIiIIGlMSACIiACIiACIiACPhOQwPIZqJoTAREQAREQAREQAQksjQEREAEREAEREAER8JmABJbPQNWcCIiACIiACIiACEhgaQyIgAiIgAiIgAiIgM8EJLB8BqrmREAEREAEREAEREACS2NABERABERABERABHwmIIHlM1A1JwIiIAIiIAIiIAISWBoDIiACIiACIiACIuAzAQksn4GqOREQAREQAREQARGQwNIYEAEREAEREAEREAGfCUhg+QxUzYmACIiACIiACIiABJbGgAiIgAiIgAiIgAj4TEACyztQY9UDWOn9EdUUAREQAREQARFwCbQD5gNrk0BEAst7lHsCc71XV00REAEREAEREIEaBDYG5iWBigSW9yi3B4rmzJlD+/b2z+xKSVkJPW6wCTGYf8F8CgsKs2tQT4uACIiACIhASAmsWLGCXr16mXUdgBUhNdNXsySwvOOsFlhFRUW+Cay2Y9pW9148ulgCy3scVFMEREAERCBiBExgdehg2koCK2Khy4u5Elh5waxOREAEREAE4kZAAituEfXXHwksf3mqNREQAREQgYQQkMBKSKAb6aYEViPB6TEREAEREIFkE5DASnb8G/JeAqshQnpdBERABERABGohIIGlYVEfAQksjQ8REAEREAERaAQBCaxGQEvQIxJYCQq2XBUBERABEfCPgASWfyz9bGk0cDSwNbAaeAe4GPgmrZOWwN+BXwOtgZeBs2skBu0N3Ars57bzEHAhUObRWAksj6BUTQREQAREQATSCUhghXM8PA88DHwANAeuBvoB2wIlrsm3AUOAEcAS4AZgA2AAUAk0Az4BFgMXAJ2Be4FJwG89ui2B5RGUqomACIiACIiABFb0xkBXYBEwGHjDTVpmwukk4BHXHUuRPgc4FPgv8EvgacDSyNo9SFaOA+4BunnMKiuBFb2xIotFQAREQARCQEAzWCEIggcT+gLT3Vmsae6Sny0J2ozVsrTnPwUmA1cAfwGOBHZMe70TsNR9/lUP/UpgeYCkKiIgAiIgAiJQk4AEVvjHhF3t8wRg4mgv19zjgfGA7cNKLy8APwBnAHcCmwAH1ahT6i4r/qcW16299DbtFvC5uion/INEFoqACIiACISLgARWuOJRmzW2Sf0wYM+0Dex1CawXgRnAma7A6gMcXKNR2+B+srvHq2Z/V7qzXz/7vQRW+AeJLBQBERABEQgXAQmscMWjpjX/BIYCe7szU6nX7VRgLpYINYMV7vEg60RABERABPwkUFkOr18HA0ZAh55+towElq84fWvMlgVNXB0F7OPuv0pv3K7ntk3uJwKPui90d2e4am5y3xhY4NY51j1JqE3uvoVKDYmACIiACESSwPLZMPFUmPsB9N4NRjwLTZv65ooElm8ofW3o34AtA9om9fTcV0VuPivrzNI0HO7up7KN65YTy1Ix1EzT8CNwkbsh3k4Q2iZ4pWnwNVxqTAREQAREIFIEvn4GJp8Na5ZDyw5w5L9g2yN8dUECy1ecvjW2to6WRrppFuzlVsD1rhBLTzRqqRpSxRKNmlirmWjUNrp7KTpF6IWS6oiACIiACESDQEUZvPgnmGJzFECPnWD4eOhkZ8L8LRJY/vKMW2sSWHGLqPwRAREQgaQSWLkQ/vNrmP+RQ2DQOXDAldC8ICdEJLBygjU2jUpgxSaUckQEREAEEk6gbBXctT+smA9Db4Otbcty7ooEVu7YxqFlCaw4RFE+iIAIiEBSCVSUQtMW6zevL5kBzQqgo11yktsigZVbvlFvXQIr6hGU/SIgAiKQVAImpiaMcDav721nvfJbJLDyyztqvUlgRS1islcEREAERAA+nwhPnQdlK6GwG/zuI2hpl5Pkr0hg5Y91FHuSwIpi1GSzCIiACCSVQPlqeP4SmGpZiYDeu8Mxd/meRNQLXgksL5SSW0cCK7mxl+ciIAIiEC0CP013lgR/nAY0gb0ugH1GQ7PmgfghgRUI9sh0KoEVmVDJUBEQARFIMIHSYripH6xeCoVd4eg7YXNLARlckcAKjn0UepbAikKUZKMIiIAIiAB8cBd8MdlZEmy3UeBEJLACD0GoDZDACnV4ZJwIiIAIJJjAoq+gsgy67+hAWLsW1lZB02ahgCKBFYowhNYICazQhkaGiYAIiEBCCZiQ+uRBeOZCaNsNznwTWnUIHQwJrNCFJFQGSWCFKhwyRgREQAQSTsD2Wj1zPnz2iANis32dJcHCLqEDI4EVupCEyiAJrFCFQ8aIgAiIQIIJLJzmnBJcMh2aNIV9L4M9z1+fpT1kaCSwQhaQkJkjgRWygMgcERABEUgcAVsStLxWz10MlaXQrgcMGwd9dg81CgmsUIcncOMksAIPgQwQAREQgYQTMIH10LEw/b/Q90A46g4o7Bx6KBJYoQ9RoAZKYAWKX52LgAiIgAhUEyhZAtMmwsDTQrskWDNSElgau/URkMDS+BABERABEcgvAZuxspxWlpF9yM357dvH3iSwfIQZw6YksGIYVLkkAiIgAqElsHo5PPlb+OpJx8STHg88I3tjWUlgNZZcMp6TwEpGnOWlCIiACARPYN5UmDASls+Cpi3gwL/AoLOgSZPgbWuEBRJYjYCWoEcksBIUbLkqAiIgAoEQsCXB926DF/8EVeXQsQ8MHw89BwRijl+dSmD5RTKe7UhgxTOu8koEREAEwkPgiXPh4/sde7Y5Ao74J7TuGB77GmmJBFYjwSXkMQmshARaboqACIhAYAS+exkePh4OugoGjorskmBNfhJYgY2oSHQsgRWJMMlIERABEYgQgaoqWDoDumyx3ujiRc69gjEqElgxCmYOXJHAygFUNSkCIiACiSVg+awmnwmz34Mz3oANNo0tCgms2IbWF8cksHzBqEZEQAREQASY9Q5M/A2snA/NWjrX3WwzJLZgJLBiG1pfHJPA8gWjGhEBERCBBBOwJcG3boRXr4G1ldB5Cxh+D2y0fayhSGDFOrxZOyeBlTVCNSACIiACCSZQvBgmnQbfv+pA2OE4OOwGaNk29lAksGIf4qwclMDKCp8eFgEREIGEE3jlKnjjemjeGg77O/Q/ITanBBuKrARWQ4SS/boEVrLjL+9FQAREIDsCFaXO1Td7/h66bZNdWxF7WgIrYgHLs7kSWHkGru5EQAREINIEVi6Ed/8F+18JzZpH2pVsjZfAypZgbp7fG7gIsHsCugNHAZPTulpbR7d/AK53X5sJ9KlR71rgkgxMlsDKAJaqioAIiECiCcx4BSadDiWLYfDFsO+licYhgRXO8P8S2AP4CHisFoG1UQ2zrf44oC/wfZrAst+NTatbDNiP1yKB5ZWU6omACIhAUglUVsBrY+DNG4C10G0755Rg1y2TSqTabwms8IffZqtqzmDVtNpmt9oB+6e9YDNYN7k/jfVSAqux5PScCIiACCSBQNE8eGwUzH7H8XbASDhkDLRonQTv6/VRAiv8Q6AhgbUhMBc4BXiohsBqCRQAc4AJ7vJhWT0uW337SRUTbXOLiopo3960VnalpKyEtmOco7nFo4spLCjMrkE9LQIiIAIiEByBH96ER0+G1UuhoB0MuQn6DQvOnpD1LIEVsoDUYk5DAsv2Xdm+qh7AmrTnf+8uMS4DdgHGAE8Ao+px+UrgipqvS2CFf5DIQhEQARHIO4FFX8PYfaFzX2dJsPPmeTchzB1KYIU5Oo5tDQmsr4EXgd824MoxwESgC7CkjrqawQr/eJCFIiACIhAcgbJVUNBmff9z3oeNdoAWrYKzKaQ9S2CFNDBpZtUnsPYC3gD6A5824EpPdylxEDDFo9vag+URlKqJgAiIQOwJfP2sk9Pq2Aegz26xdzdbByWwsiWY++frE1j3AHaZ084ezDgceMpN3TDbQ32rIoHlEZSqiYAIiEBsCVSUwUtXwnu3Oi5ufTgc92Bs3fXLMQksv0j6247tBLeUC1Y+Bs4H7CKnpUBKHJn4WQBcANxeo3v7amEzVfZMETAQ+AfwIXBkBqZKYGUAS1VFQAREIHYEls2EiafCvKmOa4POhgP+DM3t/JRKfQQksMI5PvZxxVFN6+4FRri/PN1NwWCJSE1EpZedgH/b9wz3VOAs4GHgOmBVBi5LYGUAS1VFQAREIFYEvnwSnjgXSougVUcYehtsfWisXMylMxJYuaQb/bYlsKIfQ3kgAiIgApkTmPUOjLcc1sDGA2HY3dCxd+btJPgJCawEB9+D6xJYHiCpigiIgAjEjsDatTBhBHTqA/tdDs1axM7FXDskgZVrwtFuXwIr2vGT9SIgAiLgncBXT8Ome0GrDs4zVVXQtKn351XzZwQksDQg6iMggaXxIQIiIAJxJ1C+Gp4fDVPHw7ZDnaShTZrE3euc+yeBlXPEke5AAivS4ZPxIiACItAAgZ+mO0uBP04DmsBe58O+l0HTZkKXJQEJrCwBxvxxCayYB1juiYAIJJjAp4/A07+H8hJo0wWOvhP67p9gIP66LoHlL8+4tSaBFbeIyh8REAERsOtunrsIPn7AYbHJXnDMXdBuI7HxkYAElo8wY9iUBFYMgyqXREAEEk6gZAncvgesXAiDL4bBf9CSYA6GhARWDqDGqEkJrBgFU66IgAiIwDoCs96FyjLYbLCg5IiABFaOwMakWQmsmARSboiACCSYQGkxPHshbDoY+v86wSDy67oEVn55R603CayoRUz2ioAIiEA6gYXTYOJI+OlbaNkezvscWncUozwQkMDKA+QIdyGBFeHgyXQREIEEE7BM7FPvgecvgYo10K47HDMONtkjwVDy67oEVn55R603CayoRUz2ioAIiMCaFfD0eTDtMYdF3wPhqNuhsIvY5JGABFYeYUewKwmsCAZNJouACCSYgKVguGMvWPIdNGkG+/8Jdv+drrwJYEhIYAUAPUJdSmBFKFgyVQREQASqCbx0JXw2AYbdDb13FZSACEhgBQQ+It1KYEUkUDJTBEQgwQTWFIGdFOzQ04FQWQGlK6DNBgmGErzrEljBxyDMFkhghTk6sk0EREAE5k2FCSOhTWc49b/QvEBMQkJAAiskgQipGRJYIQ2MzBIBEUg4ATslOOV2eOFyqCqHjr3h5Cdhg00TDiY87ktghScWYbREAiuMUZFNIiACySawaik8cS5884zDYZshcMS/lN8qZKNCAitkAQmZORJYIQuIzBEBEUg4gTkfOIlDi+ZAswI46GrY5TRo0iThYMLnvgRW+GISJosksMIUDdkiAiKQbAK2LDjuQJj7AXTaFIaPhx6/SDaTEHsvgRXi4ITANAmsEARBJoiACIjAOgJLv4c3boBDxkAre4tWCSsBCaywRiYcdklghSMOskIERCCpBGa9C/M/gt3OSSqByPotgRXZ0OXFcAmsvGBWJyIgAiJQg0BVFbz9D3jlalhbBac8BZvuJUwRIiCBFaFgBWCqBFYA0NWlCIhAwgkUL4bHT4cZrzggdjgWDrsRWrZNOJhouS+BFa145dtaCax8E1d/IiACySbww5vw2CgoXgjNW8Oh18MvTtQpwQiOCgmsCAYtjyZLYOURtroSARFIOIG3boKX/+wsCXbZCn51L3TbJuFQouu+BFZ0Y5cPyyWw8kFZfYiACIiAEfj4AXjiHOh/Ihx6HRQUikuECUhgRTh4eTBdAisPkNWFCIhAggmUlawXUpbnava70Gf3BAOJj+sSWPGJZS48kcDKBVW1KQIiIAKVFfD63+DTR+CM16HNBmISMwISWOEM6N7ARcAAoDtwFDA5zdR7gFNqmD4FGJT2u5bA34FfA62Bl4GzgbkZuCyBlQEsVRUBERABTwRWzHc2ss9626l+2A0wcJSnR1UpOgQksMIZq18CewAfAY/VIbA2BEammV8GLE37/9uAIcAIYAlwA2BfkUy0VXp0WwLLIyhVEwEREAFPBKa/5KRgWLUECtrCkJuh3zBPj6pStAhIYIU/XmvrEFgdgaF1mN8BWAycBDzi1ukBzAEOBf7r0W0JLI+gVE0EREAE6iVQWQ6vXAVv3+RU26gfDL8XOm8ucDElIIEV/sDWJbBMXNms1XLgdeAyYJHrzn7ukqDNWC1Lc/FTd6nxijrctmVF+0mVdrakWFRURPv22d95VVJWQtsxTqK84tHFFOqETPhHnywUARHwh4BlZH/jOqetgafBQVdBi1b+tK1WQklAAiuUYfmZUbUJrGNNowCzgE2BvwLN3eW/UuB4YHwNsWSNvgD8AJxRh9tXAv8jviSwwj9IZKEIiEDICaxeDvccDntfCNvVtfgQch9kXkYEJLAywhVI5doEVk1DbCO8ia3jgEn1CKwXgRnAmXV4ohmsQEKsTkVABGJHoKIMvpjkXHPTpInjnt0v2LRp7FyVQ7UTkMAK/8jwIrDMi+nAXcC1QGOXCGvS0B6s8I8PWSgCIhA2AstmwsRTYd5UOPwfsPOpYbNQ9uSBgARWHiBn2YUXgdUZmAecDtwHpDa5nwg86vZvs1yWokGb3LMMiB4XAREQgToJfPUUTD4HSougVQcYehtsfZiAJZCABFY4g247wfu6pn0MnA+86qZhsFQMtlfK0jcsADYBrgF6A3Zp1Ur3OUvTcLibpsGesZxYJsSUpiGcMZdVIiACUSZQUQovXA7v3+F4sfFAGHY3dLS3ZpUkEpDACmfU93EFVU3r7gXOck8C/gKwVA0mskx8Xe6mYUg9Y8dTrnf3Y6UnGrVUDV6Llgi9klI9ERCB5BJY+j1MGAkLPnEY7P5b2P8KaNYiuUzkORJYGgT1EZDA0vgQAREQgYYIzHwL7h0CrTrCUbfDlgc39IReTwABCawEBDkLFyWwsoCnR0VABGJMwC5mTp0ONDc/fRg22Qs69Iyx03ItEwISWJnQSl5dCazkxVwei4AINETgp+/gyXNhyC3QdcuGauv1hBKQwEpo4D26LYHlEZSqiYAIJITAZ4/CU+dBeQlsOhhOeTIhjsvNTAlIYGVKLFn1JbCSFW95KwIiUBeBslXw3B/g4/udGrYcePRYaG8ZcFRE4H8JSGBpVNRHQAJL40MEREAEFn8DE0bAoi+BJjD4DzD4YmjaTGxEoE4CElgaHBJYGgMiIAIiUBeBeR/BPYdB+Soo7AbHjIXNLJOOigjUT0ACSyNEAktjQAREQATqIlBZ7gisFq2dJcG23cRKBDwRkMDyhCmxlbREmNjQy3ERSDCBxd9Cp02geYEDYdVS59obLQkmeFBk7roEVubMkvSEBFaSoi1fRSDpBCy31Uf3OZvZB46Cg69OOhH5nwUBCaws4CXgUQmsBARZLoqACAClK530C9MmOjj6Hgi/fhiaNRceEWgUAQmsRmFLzEMSWIkJtRwVgQQTWPCZc0pw6Qxo0gz2vxx2/z9o2jTBUOR6tgQksLIlGO/nJbDiHV95JwLJJmBLgh+Og+cvhcpSaN8Tht0NvQclm4u894WABJYvGGPbiARWbEMrx0RABFixAP41EMpWwpaHwNDboM0GAiMCvhCQwPIFY2wbkcCKbWjlmAiIQDWBaZNgxXzY7ZyfX94sPCKQJQEJrCwBxvxxCayYB1juiUCiCNiS4JQ7oNvWShaaqMAH46wEVjDco9KrBFZUIiU7RUAE6iewehk8cS58/bSTkf2cKVoO1JjJKQEJrJzijXzjEliRD6EcEAERYO6HMGEkFM2GZgVw0FWwy+laEtTQyCkBCayc4o184xJYkQ+hHBCBBBOoqoL3boWXroSqCic7+/B7oMcvEgxFrueLgARWvkhHsx8JrGjGTVaLgAiUr4FHT4bp/3VYbHcUDLnZufJGRQTyQEACKw+QI9yFBFaEgyfTRSDRBGxD+6TT4csn4JAxsPOpWhJM9IDIv/MSWPlnHqUeJbCiFC3ZKgJJJ2BLghWroaDQIVFaDMtmwkbbJ52M/A+AgARWANAj1KUEVoSCJVNFINEEihfD42dA85Zw3EOarUr0YAiH8xJY4YhDWK2QwAprZGSXCIjAegIz34KJv4HihdC8FZz2Kmy4rQiJQKAEJLACxR/6ziWwQh8iGSgCCSZQVQlv3gCvjYG1VdBlK+eUoMRVggdFeFyXwApPLMJoiQRWGKMim0RABGDljzDpNPjhdYdG/xPg0OvX778SIxEImIAEVsABCHn3ElghD5DME4FEErATguMOhLkfQIs2cNiN0P/XiUQhp8NLQAIrvLEJg2USWGGIgmwQARH4XwKWnf2Z8+HosdB1KxESgdARkMAKXUhCZZAEVqjCIWNEIMEEViyABZ/AVr9cD8HSMjRtmmAocj3MBCSwwhmdvYGLgAFAd+AoYLJragvgKuBQYDOgCHgJuASYn+bOTKBPDfeudet59VoCyysp1RMBEcgdge9ecpKGlpXAqJeV1yp3pNWyjwQksHyE6WNT9hVtD+Aj4LEaAsvueZgIjAU+BToBNwHNgZ1rCKxxbr3Ur4sB+/FaJLC8klI9ERAB/wlUVsCrV8Fb/3Da3rAf/Ope6Ly5/32pRRHwmYAEVmZAfwDWZvZIdW0TQLc04jl7xPpLn8GqrZmBwPvujNVst4LNYFm/9tPYIoHVWHJ6TgREIDsCRXOd3FZz3nPaGTgKDroaWrTKrl09LQJ5IiCBlRnowZlVX1fbxM6sRj7rRWAdALwAdARWpAmslkABMAeYAFwPlGVghwRWBrBUVQREwCcC3zwPk8+E1cugZXs44hbnsmYVEYgQAQms8AerIYFlX+feAr4GTkxz5/fuEuMyYBdgDPAEMKoel02Q2U+qtAPmFhUV0b69aa3sSklZCW3HtK1upHh0MYWp+8Kya1ZPi4AIxI3Aa9fCa9dA9/4wfDxsYNtNVUQgWgQksPyL10hgvH/NrWupPoFlG95tZqo3sE/a7FVtZhzj7t3qAiypw84rgStqviaBlYOoqkkREIGfE7DcVk2aOL+zDO0fjIMBpzh3C6qIQAQJSGD5F7SngBuBV90mWwN3ASdk2UVdAsvE1aPuScL96hFNqe572mwUMAiYUodNmsHKMlh6XAREoBEEvnoaptwOJ0zUHqtG4NMj4SQggeVfXOx033PAqUAl8B/gjhqn+BrTW20CKyWutgD2BRZ7aPhwwESgpW5IbYRv6DHtwWqIkF4XARFoPIGKUnjxT464snLAlbCn7W5QEYHoE5DAyj6GNmv1iftjG8gfAmye24TWx41s3jYq9XWftTbOd2fGlrq5rix1w06AiaYf0/qw182G3dyZKptNszxZdsrQzjl/CByZgU0SWBnAUlUREIEMCCz9HiaMdJKHWtntXNj/Cmhu53JURCD6BCSwso+hCZYd3J+tgR6AnSt+A5gGPNOILmw/VWqpMf3xewHbJ2XpImorNpv1miu+/g2YPbbsZycYHwauA1ZlYI8EVgawVFUERMAjgS8ehyd/B6UroHUnGHo7bHWIx4dVTQSiQUACy/842d6r7V3BZf+N8ny3BJb/40MtikCyCbx3GzxvF08AvQbBsHHQYeNkM5H3sSQggZV9WP8PuBmw20an2/mX7JsMTQsSWKEJhQwRgZgQWD4b7hgMA0bAvpdCM9tSqiIC8SMggZV9TPcHXgaeBWzT+UrgC3d5sLFLhNlb5U8LElj+cFQrIpBsAgs+g+62k8Itq5ZCmw2SzUTex56ABJb/ITZRYkuD9rMtcJ7/XeStRQmsvKFWRyIQQwJlq+D5i+Gj++D4R2HLg2PopFwSgdoJSGBlNjKedvNa2cm8JBQJrCREWT6KQC4ILP4GJoyARV86B6sP/Avs8btc9KQ2RSCUBCSwMguL7a+yeW5b+ksVS6fwXS3NWKqGxlwMnZlFua0tgZVbvmpdBOJJ4JOH4JkLoHwVFHaDY8bCZnY4WkUEkkNAAiuzWJvAshtH7U4/KyaiyoH+NUTXPcBJwFT3fsBvM+smNLUlsEITChkiAhEgUFYCz1wIn1o6QGDTwXD0WGi3YQSMl4ki4FmaD+0AACAASURBVC8BCazMeJrAuh84xX3MzhZbVvQDgFfc31lGd0v4ORTYC9gd2DOzbkJTWwIrNKGQISIQAQJfPwMPHw9NmsI+o2GvC6BpswgYLhNFwH8CEliZMTWBZYLqb+49g5cBltDFEnimLknexp25agPYj91FODazbkJTWwIrNKGQISIQEQIv/wU23w82ier3yohwlpmhJyCBlVmIbPbqbuA+oBvQHDjHza5us1WWB+tCYIR7ijCz1sNXWwIrfDGRRSIQHgKlK8EE1d4XQVt7S1QRARFIEZDAatxYsMx4tu9quSuqLgKuAj4DdnRns8Y0rulQPSWBFapwyBgRCBEBy201cSQs+Q76HggnTgyRcTJFBIInIIHlXwz2ACzJiy0hjovBCUIjI4Hl3/hQSyIQDwJr18KH4+D5S6GyFNr3hGF3Q+9B8fBPXoiATwQksHwCGdNmJLBiGli5JQKNIrCmCJ76P7DLmq1seQgMvU1Z2RsFUw/FnYAEVtwjnJ1/EljZ8dPTIhAfAj99Bw8eA8tmQtPmcMCfYbdzoIllq1ERARGoSUACS2OiPgISWBofIiACDgHb0G6XNFeWw/DxsPHOIiMCIlAPAQksf4bHTm6i0TJ/mgtNKxJYoQmFDBGBAAisWQEt262fpfppOhR2gdadAjBGXYpAtAhIYPkTr0rA8l9FNWN7XRQksPwZH2pFBKJHYO6HzinBXc+C3c6Onv2yWAQCJiCB5U8ALAHp1hJY9cMsKSuh7Zi21ZWKRxdTWFDoD321IgIi4B8BOyX47q3w0hVQVQFdtoQz34bmBf71oZZEIAEEJLD8CbIElgeOElgeIKmKCARJYNVSmHwWfPu8Y8W2Q+GIW6CV3QCmIgIikAkBCaxMaNVdVwLLA0cJLA+QVEUEgiIwewpMPBVWzIVmLeGQa2Dn3+iUYFDxUL+RJyCB5U8IJbA8cJTA8gBJVUQgCALFi+GmflCxGjbYHIbfA913CMIS9SkCsSEggeVPKCWwPHCUwPIASVVEICgCb98CCz6FITc5JwdVREAEsiIggZUVvnUPS2B54CiB5QGSqohAvgjMfBtad4QNt3N6tM3tVpQ4NF8RUD8xJyCB5U+AJbA8cJTA8gBJVUQg1wSqKuHNG+G1a6BzXzjtVWjpnO5VEQER8I+ABJY/LK8A/gks9ae50LSiPFihCYUMEQEfCBQvgkmnwfevOY3teDwc9ndQyhQf4KoJEfg5AQksjYj6CEhgaXyIQFwImKh67DQoWQQt2sBhN0D/4+PinfwQgdARkMAKXUhCZZAEVqjCIWNEoBEEbEnw9Wvh9etsoxV029Y5Jdh1q0Y0pkdEQAS8EpDA8koqmfUksJIZd3kdJwImsO47Ema+CTudDIdcCwVt4uShfBGBUBKQwAplWEJjlARWaEIhQ0QgQwJ2KjB1InDlQpj5FvQblmEjqi4CItBYAhJYjSX3v8+NBMb711woWpLACkUYZIQIZECgsgJevRpKVzob2FVEQAQCISCB5R/2p4AbgVfdJlsDdwEnNKKLvYGLgAFAd+AoYHJaO00AO7l4OtAJmAKcA3yRVsd+fwtwhPu7J4HfAsszsEcCKwNYqioCgRMomgeP/QZmv+uYcvpr0OMXgZslA0QgiQQksPyLut2G+hxwKlAJ/Ae4AxjbiC5+CewBfAQ8VovAuhi4DBgBfAv8ETBRZrtWV7r9mS0buyLMfnUnMBMYkoE9ElgZwFJVEQiUwLf/hcfPhNVLoWV7GHIzbH90oCapcxFIMgEJrOyjb7NWn7g/ZcBDlgvZFVofZ9+8Hfv52QyWtT0fuAm41m2/JfAjYMLLRN02wJfAIHd2y6rZv+1r7dbANx7tksDyCErVRCAwApXl8PKf4R1LxWdz3v1h+HjYYLPATFLHIiACIIGV/Sg4ErBbUe3HxEsP4D3gDWAa8EyWXdQUWPauOQPYCUgXcE+4y3+nuOLOhF/HGn3b8uDv69krZkLNflLFLiSbW1RURPv2prWyK8rknh0/PS0CtRJ4cDhMf8F5adez4MA/Q/P0P2NxEwERCIKABFb21G3WyGaSdnRnjZoD27uCy/5rgiabUlNg7Q68DfR0Z7JSbdsSYB/gYOBSd/lwyxod23KibcQfU4dBV7p7u372sgRWNuHTsyKQYwJfPwOTz4Yjb4VtDs9xZ2peBETAKwEJLK+k6q5ne6VM8ExyZ7Aq3JmrzwH7+cBdvmtsT3UJLJspW5DWqO316gUc4gosm8mqmUlwOjAO+FsdxmgGq7FR0nMikC8CFWWwZPr6S5qt39XLnYubVURABEJDQALL/1BYBj+buern/uwFPAtc3siu8rlEWNNE7cFqZND0mAjkhMDSH2DiSFg2C858CzrYRLaKCIhAGAlIYGUWlafdtAtFmT3GVDflQoaPVVeva5P7PwC7+8JKAbColk3uuwLvu3Xs37Y3TJvcGxMFPSMCQRP4YjI8+VsoXQGtO8GxD8ImNoGuIgIiEEYCEliZRaXK3Vtlm9dTpS/wXS3N2Gk/E0dWWgDlGXTVFrB2rdhG9vPd/FpLgdmukBoNWHJTW/azPVf71JKmwZYRz3DbsT1as5SmIYMoqKoIhIFA+Rp44TL4wNLq2UaAQTBsHHSwLCwqIiACYSUggZVZZExgWdJPO7FnxUSUCaf+7r6rVGv3ACdB9czViW6uqkx6MrGUSlia/ty97ub1VKJRE0/piUbThd8GtSQaPVeJRjMJg+qKQMAElsyACafAQtvOCex5Pux7KTSz72wqIiACYSYggZVZdExg3Q/YBnIr9hXSZpQOAF5xf2cJR22maShg+6/s1N+emXUTmtragxWaUMiQRBJ49iJ4/05o0xmOvhP62luNigiIQBQISGBlFiUTWCao7BSezddbNvVL3L1QdnWNFUvyaTNXttndfuyqnMZkc8/MstzUlsDKDVe1KgLeCJSvhucuhn0ugfa24q8iAiIQFQISWJlFymav7gbuA7oBlvPK7gC0/FE2W2X7oS50l/HsJGHUiwRW1CMo+6NFYPG3zl6rQ/4GTZtGy3ZZKwIi8DMCEliNGxC2AcL2XVlmdBNVdjHzVcBnbsJRm82qK5ln43oM5ikJrGC4q9ckEvjkP/DM+VC+Cg66Gna3LZMqIiACUSUggeVf5Oy8tGVRtyVES+aZOkHoXw/5b0kCK//M1WPSCJSVgO21+uRBx/NNB8PRY6HdhkkjIX9FIFYEJLBiFU7fnZHA8h2pGhSBNAKLvoIJI2Dx19CkKewzGva6AJo2EyYREIGIE5DAingAc2y+BFaOAav5BBOwxKGPnwkVq6FddzjmLtgkqgeOExxHuS4CdRCQwNLQqI+ABJbGhwjkisCCT+GuA2CTvZwUDIVdctWT2hUBEQiAgARWANAj1KUEVoSCJVMjQKDmpczzP4GNdtCJwQiETiaKQKYEJLAyJVZ7fdvg/iFQ6k9zoWlFAis0oZAhkSawdi1MHQ8vXgEnPwE9d4q0OzJeBESgYQISWA0z8lJjhZu24XsvlSNURwIrQsGSqSElsGYFPPU7+OJxx8CdfwOH3xhSY2WWCIiAXwQksPwhudLNfyWBVQ/PkrIS2o6xe6yheHQxhQWF/tBXKyIQVgK2BGinBJf9AE2bwwFXwqBztCQY1njJLhHwkYAElj8wJbA8cJTA8gBJVeJBwJYE3x8LL1wGlWXQoTcMuxt6DYyHf/JCBESgQQISWA0iqrXCyTV+ewdwObAo7fd2nU7Ui5YIox5B2R8Mga+egkdOdPre+nA48l/QulMwtqhXERCBQAhIYDUO+6s1HrN7CG2T+2r395bFfb/GNR2qpySwQhUOGRMZAlVV8OhJTgqGXc+AJk0iY7oMFQER8IeABJY/HLVE6IGjlgg9QFKVaBKwJcGP74ftj4HU3kL7nYRVNOMpq0XABwISWD5ABCSwPHCUwPIASVWiR2DVUph8Nnz7HPQ/AYb+O3o+yGIREAHfCUhg+YNUAssDRwksD5BUJVoEZk+BiafCirnQrCUcco2ThkEzV9GKo6wVgRwQkMDyB+po4DZguT/NhaYV7cEKTShkSKgI2B6rd26Bl/8Caythg81h+D3QfYdQmSljREAEgiMggRUc+yj0LIEVhSjJxvwSKFkCj58B373o9NtvOBz+D2jZLr92qDcREIFQE5DACnV4AjdOAivwEMiA0BFYsQBu3xPKiuGX18FOJ2tJMHRBkkEiEDwBCazgYxBmCySwwhwd2ZY/AjVPBM58y8lrteF2+bNBPYmACESKgARW9uG6GLgWsM0XXwHl2TcZmhYksEITChkSGIHiRTDpdBgwArYbGpgZ6lgERCBaBCSwso/XHsDbwCRgG1dgTQM+d38+AH7MvptAWpDACgS7Og0Nge9fh0mnQfGP0HYj+L9PoUWr0JgnQ0RABMJLQALL/9i0Aba3ra/uj2V5f9a9Ssf/3nLbogRWbvmq9bASqKqE16+D121yei103cY5Jdht67BaLLtEQARCRkACK/uA/B9wM7AVMB2oqqXJqcCA7LvKewsSWHlHrg4DJ7ByITw2Cma+6Zhim9gPuRYK7LuTigiIgAh4IyCB5Y1TfbX2B152Z6m2cLO6fwHYMqH9PAO0iOjeLAms7MeHWogSAcvKfusuULIYCtrC4TfBDsOj5IFsFQERCAkBCSz/A2GixJYI7Wdb4Dz/u8hbixJYeUOtjkJD4LmLYebbzpJgl76hMUuGiIAIRIuABFZm8XoaOAEoyuyxyNaWwIps6GS4ZwJF86BJU2jf3XmkohQsLYM2s3tGqIoiIAL/S0ACK7NRYfurLB2DLf2lin3F/a6WZpo4u2NzVmYCfWpp3W6aPQd4DRhc4/VHgOMysEgCKwNYqhpBAt++4GRl77o1nPIUNGseQSdksgiIQBgJSGBlFhUTWEcBT7iPmYiyvFf9a4iue4CTANvcfiLwbWbdeKrdFWiWVtOWJO3ujn1dcWUCy/r9U1qd1RnOvklgeQqFKkWOQGW5c4+g3SdopXt/OGEitLU/KxUREAERyJ6ABFZmDE1g3Q+c4j62MTAbOAB4xf1dB2ApYBkJLUXD7sCemXXTqNo3AYcDttHeZs5MYH2S5R4wCaxGhUIPhZrA8tkw8VSYaynqgF3OgIP+Cs1bhtpsGScCIhAtAhJYmcXLBJYJqr8BdwGXAZcA1wFXuE1ZslGbubIz3fZje7bGZtZNxrULgPnAjcA17tMmsOweD5tls0SnzwF/dk85eu1AAssrKdWLBoGvn4HJZ8Oa5dCqAxx5K2wzJBq2y0oREIFIEZDAyixcNnt1N3Af0A2wDRu23+lKd7bK8mBdCIxwTxFm1nrja/8KeAjo7Qota+k04AdgoWvLGHev2IH1dGNf4dO/xrcD5hYVFdG+vWmt7EpJWQltx7StbqR4dDGFBYXZNainRSATApUVcMdesOhL6DkAho2HTrVtY8ykUdUVAREQgdoJSGA1bmRYXivbd7XcTS56EXAV8BmwozubZYImX+W/QBlQ31dxS3T6oZvw9KM6DDOhmJqJW1dFAitfYVQ/OSew6Gv49CHY94/Q3CZ+VURABEQgNwQksPzjancSHuwuIY7L8QnCdKvtK/j3wNFpm+9r88qWCkvdzfd2mrC2ohks/8aDWgoDgS+fgBXzYdBZYbBGNoiACCSIgARW9INts05nAL0sg0897tgpQ7uA2lI3vOHRbe3B8ghK1UJGoHwNvPBH+GAsNGkGo16CnjuFzEiZIwIiEGcCEliZRdf2NDUmt5Wd8HPPg2fWYQO1m7r7rP7jbrZPVd/c3Vxvl0z/5GaUvwGwNA0DgUqPVkhgeQSlaiEisGQGTBgBC23FHtjjPNjvj9DMVvZVREAERCA/BCSwMuNcM3Gn16ctKegsr5UzqHcQYPuv7KLp9FxbNpv1gLu53XaVz3HvRLRThJZCwmuRwPJKSvXCQeDzifDUeVC2Etp0hqPuhC0si4qKCIiACOSXgARWfnlHrTcJrKhFLMn22h2CU253CPTZA465C9r3SDIR+S4CIhAgAQmsAOFHoGsJrAgESSa6BN4fC89eBHtfCIMv0bU3GhgiIAKBEpDAChR/6DuXwAp9iBJu4Opl0LqTA8EuaLZ9V90tU4qKCIiACARLQAIrWP5h710CK+wRSqp9ZSXw7B9g1ttwxutOVnYVERABEQgRAQmsEAUjhKZIYIUwKIk3adFXzinBxV9Dk6Yw/F7Y9ojEYxEAERCBcBGQwApXPMJmjQRW2CKSZHtsCfCTB+GZC6FiNbTdyNnIvqndqa4iAiIgAuEiIIEVrniEzRoJrLBFJKn2lBbDM+fDZ+4lBJvv56RgaNs1qUTktwiIQMgJSGCFPEABmyeBFXAA1L1L4Onz4cNxTlb2/S6DPX4PTS3ProoIiIAIhJOABFY44xIWqySwwhKJpNtR8hM8dCwcdBX02S3pNOS/CIhABAhIYEUgSAGaKIEVIPxEd71mBUx7DHYeuR6D7cFqYneWq4iACIhA+AlIYIU/RkFaKIEVJP2k9j3/E+eU4LIf4Kg7YMfjkkpCfouACESYgARWhIOXB9MlsPIAWV24BGyGyrKxv3AZVJZBh94w7G7oZfeTq4iACIhAtAhIYEUrXvm2VgIr38ST2t/q5fDkufDVUw6BrQ6Dobeuz9KeVC7yWwREILIEJLAiG7q8GC6BlRfMCe9k3lSYMBKWz4KmLeCgv8KuZ2q/VcKHhdwXgagTkMCKegRza78EVm75qnUjMP0lePAY6NgHho+HngPERQREQAQiT0ACK/IhzKkDElg5xZvgxmueCPx8IvQ9AFp3TDAUuS4CIhAnAhJYcYqm/75IYPnPVC3OngJP/x5+/RB02kQ8REAERCCWBCSwYhlW35ySwPINpRqiqgreuQVe/gusrYTtjnaWBFVEQAREIIYEJLBiGFQfXZLA8hFmopsqWQKPnwHfvehg2P4YOPwmaGVDTEUEREAE4kdAAit+MfXTIwksP2kmta1Z78DE38DK+dC8FfzyWtjpFJ0STOp4kN8ikBACElgJCXQj3ZTAaiQ4PeYS+O5leHC4syTYZUsYfg9suJ3wiIAIiEDsCUhgxT7EWTkogZUVPj1MRSmMOxC6bQuH/h1athUUERABEUgEAQmsRIS50U5KYDUaXYIfnDsVuu8IzZo7EEpXQst2CQYi10VABJJIQAIriVH37rMElndWqllVCa9fB69fC3tfBPtdJiYiIAIikFgCEliJDb0nxyWwPGFSJVYuhMdGwcw3HRi2iX3IzdrIrqEhAiKQWAISWIkNvSfHJbA8YUp4pRmvwKTToWQxtCiEITfBDr9KOBS5LwIikHQCElhJHwH1+y+BpfFRN4HKCnhtDLx5A7AWNtzeOSXYZQtREwEREIHEE5DASvwQqBeABJbGR90ElsyA2/aAitWw86lw8DXQorWIiYAIiIAIABJYGgb1EZDA0vion8CnjzinBS0zu4oIiIAIiMA6AhJYGgwSWBoD3ghUlsMrf4WtDoXeg7w9o1oiIAIikFACEljRDfyVwBU1zP8R2Mj9XRP39dOBTsAU4Bzgiwxc1gxWBrBiXXX5bJh4Ksz9ANpvDOd+AAVtYu2ynBMBERCBbAhIYGVDL9hnTWANAw5IM6MSWOz+/8WAJSIaAXwL/BHYG9gKWOnRdAksj6BiXe3rZ2Dy2bBmObTsAEf+E7Y9MtYuyzkREAERyJaABFa2BIN73gTWUKB/LSbY7NV84CbgWvf1loDNcJnwusOj2RJYHkHFslpFGbx0Bbz3b8e9HjvB8PHQaZNYuiunREAERMBPAhJYftLMb1smsC4CiuwyEncJ8FLge2AzYIalewQ+TjPrCWA5cIpHUyWwPIKKXbU1RXDfUJj/kePaoHPggCuheUHsXJVDIiACIpALAhJYuaCanzZ/CdgmGFv+29BdAtwa2M5dBnwb6OnOZKUsuhPoAxxch4k2y2U/qWIXyM0tKiqifXvTWtmVkrIS2o5xLvstHl1MYUFhdg3q6dwRWLsWHj0JfngTht4GWx+au77UsgiIgAjEkIAEVnyCamrFZq2uA94DTGD1ABakuTgW6AUcUofbtW2cRwIrPoOkXk/K10BVBbR0RDCrlzsXNXe0IaMiAiIgAiKQCQEJrExohb/ui8B3wPWNXCLUDFb4Y5wbCy1p6IQR0LkvDLtbdwjmhrJaFQERSBABCaz4BNvEkc1g2TLgX92lwX+4M1rmpW2eWaRN7vEJuG+efD4RnjoPylZC6w3gjDc0a+UbXDUkAiKQVAISWNGN/N+Bp4DZQDd3D9ZgoB8wyxVSo4GRwHTANsDvozQN0Q2475aXr4bnL4Gp9zhN994djrkLOtjWPRUREAEREIFsCEhgZUMv2GcfdvNadXFzX9m+q8uBL12zUolGz6iRaHRaBmbrFGEGsCJVdfG3zpLgIss72wT2ugD2Ge1ce6MiAiIgAiKQNQEJrKwRxroBCaw4hreqEm7dBZZ8B4Vd4eg7YfP94uipfBIBERCBwAhIYAWGPhIdS2BFIkyNMHLmW/DG9XDUHdAudbtSI9rRIyIgAiIgArUSkMDSwKiPgARWXMbHoq9g6Q8/z2dlua6a2EqyigiIgAiIgN8EJLD8Jhqv9iSwoh5PE1GfPAjPXAhNmsIZr0OXLaLulewXAREQgdATkMAKfYgCNVACK1D8WXZeWgzPnA+fPeI0ZPusjroT2nbNsmE9LgIiIAIi0BABCayGCCX7dQmsqMZ/4TSYcIqzkd1mrva9DPY8H5o2japHslsEREAEIkVAAitS4cq7sRJYeUfuQ4cfjofnLobKUmjXA4aNgz67+9CwmhABERABEfBKQALLK6lk1pPAimLcX/wTvH0zbHEQDL0dCjtH0QvZLAIiIAKRJiCBFenw5dx4CaycI/apg6qq9ct/leVg19/scKyWBH3Cq2ZEQAREIFMCEliZEktWfQmssMfbTgm+Pxa+nAwnTYbmduWkigiIgAiIQNAEJLCCjkC4+5fACnN8Vi+HJ38LXz3pWHnEP2Gnk8NssWwTAREQgcQQkMBKTKgb5agEVqOw5eGheVNhwkhYPguatoAD/wKDzlLi0DygVxciIAIi4IWABJYXSsmtI4EVttjbkuB7t4FtZK8qh469Yfg90HNA2CyVPSIgAiKQaAISWIkOf4POS2A1iCjPFV68At6+yel0myFwxL+gdcc8G6HuREAEREAEGiIggdUQoWS/LoEVtvgvmQHjDoR9RsPAUVoSDFt8ZI8IiIAIuAQksDQU6iMggRX0+LD0C3Pfh96D1ltiV+C0bBu0ZepfBERABESgHgISWBoeElhhHQMlS2DymTD9RTh5Mmy2T1gtlV0iIAIiIAI1CEhgaUhIYIVxDMx6Byb+BlbOh+atYMgtsOOxYbRUNomACIiACNRCQAJLw0ICK0xjwJYE37oRXr0G1lZC5y2cU4IbbR8mK2WLCIiACIhAAwQksDREJLDCMgaKF8Gk0+H7Vx2LdjgODrtB+63CEh/ZIQIiIAIZEJDAygBWAqtqk3s+g/7Jf5w9V81bw2F/h/4n6JRgPvmrLxEQARHwkYAElo8wY9iUBFY+g2pJRF+9GrY/Brptk8+e1ZcIiIAIiIDPBCSwfAYas+YksHIZ0JULnYzsh/wN2myQy57UtgiIgAiIQJ4JSGDlGXjEupPAylXAZrzi7LcqWezMWA27O1c9qV0REAEREIEACEhgBQA9Ql1KYPkdrMoKeG0MvHkDsBa6beecEuy6pd89qT0REAEREIEACUhgBQg/Al1LYPkZpKJ58NgomP2O0+qAkXDIGGjR2s9e1JYIiIAIiEAICEhghSAIITZBAsuv4MydCg8Og9VLoaAdDLkJ+g3zq3W1IwIiIAIiEDICElghC0jIzJHA8isgq5fB7XtD647OkmDnzf1qWe2IgAiIgAiEkIAEVgiDEiKTJLCyCUbJT9Cm8/pcVktmQPue0KJVNq3qWREQAREQgQgQkMCKQJACNFECq7Hwv34WJp8FB/4ZBoxobCt6TgREQAREIKIEJLAiGjhgNHA0sDWwGrCd0xcD36S59BowuIaLjwDHeXRbAssjqHXVKsrgpSvhvVudX/XZA055Gpo2zbQl1RcBERABEYgwAQms6AbveeBh4AOgOXA10A/YFihx3TKB9S3wpzQ3TYwVeXRbAssjqOpqy2bChJEw/yPnqUFnwwF/huYFmbSiuiIgAiIgAjEgIIEVgyC6LnQFFrkzVm+kCaxPgPMa6aYElldwXz4JT5wLpUXQqgMMvQ22Pszr06onAiIgAiIQMwISWPEJaF9gujuLNS1NYG0HNAF+BJ4D/gysrMPtloD9pEo7YG5RURHt25vWyq6UlJXQdkzb6kaKRxdTWFCYXYNhedo2r/9rZ1hbBRsPdLKyd+wdFutkhwiIgAiIQAAEJLACgJ6DLk1APQF0AvZKa/804AdgIbA9MAb4DjiwDhuuBK6o+ZoEloeIvX4dlK6E/f8EzVp4eEBVREAEREAE4kxAAise0bUd1bYetafNONXj0gDgQ8shDrgbhX5WWzNYXsfDtEmw0Q7QxSYOVURABERABETg5wQksKI/Iv4JDAX2dmer6vPIZrpKgZMAO03YUNEerJqEylfD86Nh6njYsB+MelFX3TQ0ivS6CIiACCSQgARWdINuYsnE1VHAPu7+q4a8sWXCz2tshK/vGQmsdDo/TYcJI+BH2+LWBPY6H/a5FJrZIU4VERABERABEVhPQAIruqPh38DxwJE1cl9ZCgZLxWB3sZwAPAv85KZvuMF9bSBQ6cF1CawUpM8ehafOg/ISaNMFjr4T+u7vAaGqiIAIiIAIJJGABFZ0o762DtNHAvcAvYAH3M3tdnRvDvCMe4pwqUe3JbDK18CzF8LH9zvINtkLjrkL2m3kEaGqiYAIiIAIJJGABFYSo+7dZwmsynIYfyjM/QAGXwyD/wBNm3knqJoiIAIiteTMIAAAHilJREFUIAKJJCCBlciwe3Y6mQJr7Vqwn9T1NsvnwNLvYbOatw555qiKIiACIiACCSMggZWwgGfobvIEVmkxPHMBtO0KB12VIS5VFwEREAEREAGHgASWRkJ9BJIlsBZOg4kj4advoUkz+O1U2GBTjRAREAEREAERyJiABFbGyBL1QDIEli0HTr0Hnr8EKtZAu+5wzDjYZI9EBVvOioAIiIAI+EdAAss/lnFsKf4Ca80KePo8mPaYE7++B8BRd0BhlzjGUz6JgAiIQFYEqqrWUrS6nI5tWrB8VTkdWregaVNLy6hSk4AElsZEfQTiLbCqqmDsPrDgU2dJ0O4R3P136ze3a2yIgAiIgAisI2Di6tLHP+fhD+awU++OfDR7OccN7MU1R/WTyKplnEhg6Y8nuQLLPP98Irx4BQy7G3rvqtEgAiIgAiJQB4FlJWX84q8vVr9qc1apZIwfX34gnQoLxK0GAQksDYlkCaw1RbBsFnTfYb3fZSVQUKiRIAIiIAIiUA+B8opKdh3zCktLytbV6ti6Bbee8At237wLTZpoqTAdnwSW/pySI7DmTYUJI52N7Ge+BW27KfoiIAIiIAIeCJRVVHHWA1N5+etFtdbeoltbTt6tD0fttDFtW+p+VoMkgeVhYCW4Sjz2YNkpwSm3wwuXQ1U5dOwNv34YNtwuwaGV69kS0GbfbAnq+agQWFNeyZkPTOW1bxavWxpM7cHq27WQectXs7q8qtodE1dH79SzWmz17dYuKi7mxE4JrJxgjU2j0RdYq5bCE+fCN3YNI7DNEDjiX9C6Y2yCJEfyT0CbffPPXD0GQ6C4tIJR937Ae98vpVWLptx24gD6b9zxZ6cIi8sqmDR1Lve9N4vvF5esM3T3zTtz8m6bcMA23WjerGkwDgTYqwRWgPAj0HW0BdacD5zEoUVzoFkBHHQ17HIaaJ9ABIZeuE1M3+ybbukpu/Vh0y6F9OjYet1PpzYttDcl3OGUdXUQKFpVzinj3+eTOcurZ6bGnbIzu27WuU5ea9eu5e3vlnDfuzN56asfqXJ3wffo0IoTBvXh2IG96NK2ZWJ4S2AlJtSNcjTaAmvSGfDZw9BpUxh+D/To3ygIekgEahKwD5Jjbnun+ph6Q8W+9ffo4Aiu7h1aucLL+W/36t+3ok2B9qw0xFGv55fAT8WlnDTufb5asKI619V9p+7Cjr28z/zbsuGD782qTumQ2hRf0Kwph/bbiJN226Q6zUPcN8VLYOV3zEatt2gLrNKV8OoY2OcSaGWuqIiAPwR+dly9iXM3uJXjd+ld/WGyoGg185avwT6kvBSb5XLEliO4UmKspzsT1q1dy0QusXhhpzr+E1hYtIYT7nqPGYtL6NK2gPt/syvbdG/ce6jt33r28wXc9+6s6pmwVNm+Z3tOHrQJR/TvQasWzfx3IgQtSmCFIAghNiFaAmvWu/D5o3DYjVoGDPGgioNpXvdglVZUYh9W85fbz2rnp/r/V1eLMPu97XFpqFii7I3at6J7aunRnQlLzYiZELPM2nGfEWiIk17PnsCcpas4/q73mLN0dfWM64OjdmWzrm2zbxj4bO7yaqH15KfzsVOJVmzc/mrnXpy4ax96d27jSz9haUQCKyyRCKcd0RBYlpH97X/AK1fD2koYcgsMOCWcRGVVbAj4dYpwxZry9eLLFWILitZUn8wyEbZg+RoqUptZ6qFXvRRpAsxddrQZsdQMWHebFevQmtYF8ZwpiM2gCtiR7xYVc+JdU1i4Yg19Orfhgd/sSq8N/Bc9Nsv76IdzeOC9Wcxdtrraa9sau+9W3Thptz4M3qJrLDLDS2AFPKBD3n34BVbxYnj8dJjxioOy36/g8BuhZbKPB4d8XMm8DAiYkLOlRhNcNuPlLD86wmu+OwuWyVJkau9Xz9RSZEcTYq0wQaalyAwCE7OqX85fwUnjprCkpAzLafXAqF3ZsH2rnHpZWbWW175ZxL3vzuKNbxev68vE3UmD+jB8QC86tGmRUxty2bgEVi7pRr/tcAusH96Ex0ZB8UJo3hoOvR5+caKWB6M/7uRBhgRsn0v1UqQruBZUL0U6giy1NFlSVtlgq82aNmHDdi3XnYC0mS+bBUttxrdZMC1FNogxchU+nr2MU+5+nxVrKtiuR/vqPVcb5Pnqmx9+Kqme0bKZrZVrnGVzm5Ud2r9n9azWdj06RI6rBFbkQpZXg8MrsN4fC8/9AdZWQdetnVOC3bbJKxx1JgJRIWCnHu3DM7X3yzbgV4uwtD1hJtC8LEW2btFs3UZ8E1zVy49pS5P277huWo5KvDOx890ZS6rzXJkAH9CnE3ePGFh9ajCosqqsgic+mc+978zk64Ur15lhtlny0l9u352C5tHIqSWBFdQoika/4RVYluNq/CGww3Fw6HW6SzAa40lWhpiALdekliKrlx/XzYKtX5r8qXj9HXT1uWKzH3YaMrUPrGZ6im7tWmGzZSrBEnj1m0Wcef9USiuqsKSgY0/emcKQXHNjXwo+nLWselP8c58vWCf+LY/Wr3fpxfG79q4eX2EuElhhjk7wtoVLYBUv+vn9gYu/ga5bBU9JFohAQgisW4pMm/lKnwWzf6/yuBRppyJTImxdeoq0VBU2i6JTkbkbWCZafvfwx5RXrmW/rbvx7xN2Cu3M46IVa/jP+3N4cMosFq10Up+YQD9o2w2rlw9326xzKMeKBFbuxm8cWg6HwKqsgNfGwHu3waiXYMNt48BWPohA7AhUL0WurnD3f61PSVG9NLncORn54wpvS5FtCpqtT8yaStTq7glL5QnTUmTjhtCkj+Zy4YRPqzOtH9avO/84tn8klt3KK6t44YsfqzPFT/lh6Trnw3rRtARW48ZnUp4KXmCtmO9sZJ/1tsN83z/C4IuSwl9+ikDsCNhS5OKVpetFmHs60tkf5ixN2kk2L6VzYYGzB+x/krQ6KSq6tmuppcgaIG0W6I+Tp1Unxx02YGOuPWaHSDL6ZuHKaqH1+Mfz1s2a2nU+x+zkbIoPw0XTElhe/oqTWydYgTX9JScFw6olUNAWhtwM/YYlNxryXAQSQsCWIk1s2Ub8/0lP4YowL0uRze1UpLsUue5+yHVJWh0R1r5181AuL+Ui1GPf+J6rn/2qumm7N/OKIdtFPt+U5ZF7bOpc7n93Ft//tP6i6T36duakQcFeNC2BlYtRHJ82gxFYleXwylXw9k0OyY36wfB7ofPm8SErT0RABBpNwJYii1Zbgtb1WfGrT0ZWp6ZwNuVbskybLWuo2FJkzauJbFN+dXoK9/7IqC9FGq+bX57OTS9Nr8Zx1j6b84eDt4qVsLR8cW/P+Kl6U/zLIbloWgKrob++ZL8ejMCaeg889X8O+YGnwUFXQYvcJrxLdpjlvQjEj4CJq0Ur119TlLqaKHU60vaEZbIUmdqIv+5kZFp6ijAvRZq4uubZrxj75g/VQb7o4K04Z9++8Qt4mkdzl63ioSmza71o+uTdN+EXvfJz0bQEVqyHWdbOBSOwqirh0ZOh33DYbmjWTqgBERABEaiNgC1FpvZ+rcuOXyM9xeryhhO02lLkRrb0mLqmqMadkSbO2rfK/1Kkzepc/sQ0Hpwyu9r9Px2+LafuuWliBkPqomnLFP9pzYumd9uEI3bM7UXTEliJGWqNcjQ/AquiDKbcDrucrpmqRoVJD4mACOSCQGopsubVRKns+LZPzOtSZGFqKTLtaqL09BQm0PxciqyorOIPEz9j0sfzqu/5+9vR/Th2YO9cYIpEm/VdNH38Lr2rk6vaLQXLV5VX/7upD3naJLAiMTSyMvJsmxUGugNfAOcBb3psMfcCa9lMmHgqzJsKA0fBYTd4NE3VREAERCB4AiZkFheXrtv7lT4jlvq3XW7spXRpawla7Wqi9Oz4rddlzu/atmW9H/ypC8hbFzTl7Ac+5pVvFlWfELzxVztyZP+eXkyIfZ3URdO2Kd6Ec3rZvGshMxaXcNzAXlxzVL+sRZYEVryH07HA/YCJLMtzcAYwCrBEUs6ccf0ltwLLLmiefA6UFkGrDnDkv2GbwxuySa+LgAiIQKQIrC6zU5Hu3ZDrNuI7aSlSs2NeliJbNEudimxNj5QIq16ObFW9RHnXmz8w6aN51cuRdjWS5cq/7cSdOGR7+36tkk7A9ui9+vUixr39A3ZdUM3y8eUH0inL+xglsOI95qYAH9mhkTQ37YzuZGC0B9dzJ7B2+j8Kp453TNh4IAy7Gzomd/raQyxURQREIKYEbCnSlqZ+dkF3+sXdlqB1ZamnU5G5EAoxxV7tlrE/7Ja3+GrBCuzMqa0M9u/VkcfO2j3rU5YSWPEdOQXAKmA48HiamzcD/YHBHlyvFljzF8+nfXv7Z3alpKyEDW/YsLqRH9e2pdC+X+16JuxzKTQL7nLR7LzS0yIgAiKQewKppciFJryKSqtnxBZWb8gvZeEKJz1F0eqKdYbYvqsdN+7AQ6MGZS0Ucu9dcD0sKylj92tfqTbAxFUqs8c7F+/nywxWj649rOkOwIrgvMxfz0m5YdSiOg/YA3gnDe+llmMOqO0Sv5aA/aRKO2AulwDKkpC/EaqeREAEREAEok9gDfC3ajcksKIfzZ95kBJYuwPvpr1yGXASsHUt/l4JXPE/v5fAitnQkDsiIAIiIAI5JyCBlXPEQXXQmCXCWmewcrJEeMGPFBYUBsVG/YqACIhA7AikThH6nW4gdqDy5JDtwdISYZ5gB9CNbXKf6p4iTHX/JfBE4JvcRxdLYAUwINSlCIiACIhAfghok3t+OAfVSypNw5nuMuHpwGnAdsAsD0bl7hShBJYH/KoiAiIgAiIQVQISWFGNnHe7LQfWH9xEo9OA3wNveHxcAssjKFUTAREQAREQgXQCElgaD/URkMDS+BABERABERCBRhCQwGoEtAQ9IoGVoGDLVREQAREQAf8ISGD5xzKOLUlgxTGq8kkEREAERCDnBCSwco440h1IYEU6fDJeBERABEQgKAISWEGRj0a/EljRiJOsFAEREAERCBkBCayQBSRk5khghSwgMkcEREAERCAaBCSwohGnoKyUwAqKvPoVAREQARGINAEJrEiHL+fGS2DlHLE6EAEREAERiCMBCaw4RtU/nySw/GOplkRABERABBJEQAIrQcFuhKsSWI2ApkdEQAREQAREQAJLY6A+AtUCa86cObRvb//MrpSUldDjhh7Vjcy/YL4ue84Op54WAREQAREIMQETWL169TILOwArQmyqb6Y18a2l+DfUE5gbfzfloQiIgAiIgAjkjMDGwLyctR6ihiWwvAfDWNmU00rvjzRYs50r2mzA+dlugx3nsIJ8yiFcn5tWrHwGmqPmFKccgfW5WcWpYaDGaD6wtuGq0a8hgRVsDKuXHWM2ZSqfgh1TmfSuWGVCK7i6ilNw7DPpWXHKhFYC6kpgBRtk/UEGy99r73GMk/keR7/kk9dRHWw9xSlY/l57j2OcvPqedT0JrKwRZtVAHAevfMpqSOT1YcUqr7gb3Zni1Gh0eX1Qccor7vB3JoEVbIxaAqOBMUBpsKb41rt88g1lzhtSrHKO2JcOFCdfMOa8EcUp54ij1YEEVrTiJWtFQAREQAREQAQiQEACKwJBkokiIAIiIAIiIALRIiCBFa14yVoREAEREAEREIEIEJDAikCQZKIIiIAIiIAIiEC0CEhgBRuvs4GLgO7AF8B5wJvBmuSpd9uYfzSwNbAaeAe4GPgm7enXgME1WnsEOM5TD8FUuhK4okbXPwIbub+zvxd7/XSgEzAFOMeNXTAWN9zrTKBPLdX+7doehTjt7f6dDHD/Vo4CJqf55CUuFq9bgCPc554EfgssbxhhTmrU51ML4CrgUGAzN1feS8AlbpLGlEG1xfZat15OjG6g0YbidA9wSo027G9oUNrvbKP434FfA62BlwF7nwzyFo2G/KoraeYfgOtd38IUKy/v317i0Bu4FdjP/Rx4CLgQKAti8IWxTwms4KJyLHC/++bxNnAGMArYFpgdnFmeen4eeBj4AGgOXA30c20vcVuwD+5vgT+ltWhizBKrhrWYwBoGHJBmYCWw2P1/E5GXASNc3/4I2JvvViHOxN8VaJbmz/bAi8C+gMUoCnH6JbAH8BHwGFBTYHmJy3OA3Zhg4tjKnYB96A0JaDDW55Pd1TYRGAt86or5m9y/tZ3T7DX7x7n1Ur8uBuwniNJQnExgbQiMTDPOPoyXpv3/bW5M7G9sCXADsAFg4tr+FoMoDfmV+gKWss3qW1z6At+7vwxTrLy8fzcUB3tP+cR9b7wA6AzcC0xyv7gEEafQ9SmBFVxI7JubfWCclWbCV+43c/uGEaViH+KL3BmrN1zD7YPb/gBtVi4qxQTWUKB/LQbb34pd8WAfdDZLYMW+5dkMl33A3xERJ83+w4Et3OsqohYnmy1IF1he4rIN8KU7U2J/d1Zs1uRddxY2feY1iDDW9Kk2GwYC77uzkakvYPahbfG0n7CV2nwygdXR/RurzV4TlvZl5iTAZrut2PVkc9zZvP+GwEkvsbLZVbsSZv80e8Mcq5rv317iYCLyacBub7b3RSu2OmEx7paUy5wbGo8SWA0Rys3rBcAqYDjweFoXN7sf7jWX1nJjhX+t2je16e4s1jS3Wfvg3g6wMWYixGYQ/hzimR4z2wSWLdnaLJvlJbMP40vdb6G2VDMD2An4OA3dE+4yU82lD//o+teSjTt7M7wRuCaicar5AeclLqe6PtuHe3qx5cHfA+P9Q9yolrx8aNus6guuQFnh9mIf2ibyLa4mQia4S1JhWKKpS2DZFxizz9i/7s4I25czK7bUZEuCNmO1LI2kzeKZaKm5fN8o2Fk+1FCsbIbOljPt/cCWzFIlzLGq+f7tJQ5/AY4Edkzz0ZbhbTbSnn81S86xeFwCK5gw2rcyu03clj1s/1Kq2Ie5/WHaklNUio0hExn2x7VXmtGnAT8ACwFblrJkqt8BB4bYMftW1sZd/rM3SlsCtH1mJhQtJraU27PGPhhbarI9TgeH2K+Uab9y3/Rt70TqW2fU4lTzA253D3GxvytbctqyRoxsCdvElY3NIEtDH9qtgLeAr4ET0ww1cWiz4CZGdnH9sL9F22oQdKnNJ9sWYcuXs4BNgb+6y562/GdfaI5342GiMb2YsLT3EttGEXRpKFa278r2ytl7/JoIxKq2928vcbD3vU2Ag2oExOJof2v/CTpQYehfAiuYKKQEln042DJFqtj+Hpsetw/1qBTb5HgYsGcDG1HtTfRDdy+FfShEoRS6s1bXAe+5H+QWuwVpxts+GZsmPyQCDtkSi80e1LfvKOxxqktg1ReXur642Kyr7ZX5W8Cxq+9D2za828yUieJ9Glh6Ocbdu9XF3b8UpFsNCRGzzQ73mNiypSXbu1PXB7vtGbTZ4zODdMjtuyG/TASbvXaAor4SlljV9v7tJQ51fbG095eT3T26IQhXsCZIYAXDPy5LhP9091PYRm/7hllfsbFm327S91cEQz+zXu3N0mbe7DRQlJcIbZbNNtza6U+b5airhD1OSVoiNHH1qHuS0JZdbNN3fcVmV215yvaXpfaaZTba/avdkBBJ9WQi9y53X6OXpSn/LGxcS/X5ZTP4tgfV9nDasmbYY1XX+7eXOGiJ0MP4kcDyAClHVewNcKp7ijDVhW3EtQ+/sG9yt3Fjf5y22di+VdubZEPFlgk/r7ERvqFngn7dlipMVNm3NVvOsGW1fwA2o2XFhLLtH4nCJnfbX2ZLLDbbVlEP2LDHqa5N7vXFJbXJfVd3o7i5b/+2WUmbLQ7jJveUuLLDCHbiM3WStb6/CTu88FSNjfBB/Q15EVh28sy2StjJzvuA1OZqWwY1YWnFZrlMNFrKirBvcrcN3vb3k37Ssy7+QcaqofdvL3H4//bOJeS6KY7DT18ilyglAxNRMnCNyH1iIiUTRQYkl8jEVwbU5xIydCmRYiBmIvRJLiEpJlLIJZeSS6EoSe761dpfu9N7zr6c/Z29z3ueXafe3rPP2r/1rHX2+p21/uu/qyD37MytZvSzBJydhAa5l1bXYI11+4EqTUOmvbNMmJtM4mES75Np8ykfyaGUaeQEOdYHpwSHJxXD0cDlwIvATyV9Q7Zb573shhpru3UT0+TfyQCVXVq5SSQGKxsOkoIibRIjFfObbeYxlVl6isGccpqG1HlHmWFMXETiQ6pjXdrpoLLlPbqzwWBnCaJNQG3aqk27ZJNFlhGrOJ6Y5rTpWGkaFtUpRj7pKLKhIgNxNolUR+qcZZgzykxVgonzvcv3KiYzy/D5Xo5xLKpTdMfkp14ZkBO/k40WWfqMAf61CE56gNQ5cTz5TL6TMWJjpmlo6n+RfnCpV1IWPDIDf2pt1XT/jvymdqjSNKRvZmNQNibEYGYzQtPy6Bh9c5RrarBGwb7nokmgl6DI/ErL7rsErVZpDsZVtvjq8xLrxXjkS5ZZkifLr7ncnLLDaXfZRVjPeTO1Oia3V5Y7E8OSGYPMcOwqW/yjtUpomUG6nmi02jk5tfpUehKIml//MYIJ7K6OdWmnmNitdiXl13IG4jbtkgFgNtHojSMmGl1UpxiReUvuVf6ymK8MlJmBy0xrzGL6b2ZXs0N5jGNRnZKOJoPvyWUnZExW2jTfr9wfqiMB/VmOzw+4eqLR+jmrrltT/4ue/EBOuozcy2dz/U2trZru36lPm3aIOU4fnE00mlAQj3JjEoQEJCABCUhAAhKQwIAEnMEaEKZFSUACEpCABCQggRDQYNkPJCABCUhAAhKQwMAENFgDA7U4CUhAAhKQgAQkoMGyD0hAAhKQgAQkIIGBCWiwBgZqcRKQgAQkIAEJSECDZR+QgAQkIAEJSEACAxPQYA0M1OIkIAEJSEACEpCABss+IAEJSEACEpCABAYmoMEaGKjFSUACEpCABCQgAQ2WfUACm0sgjwDJo43yXLjqyCOO8hzJ+pHn+OVhu/OOrcpZluq8Mvvoq2tZF617Q+eybeLnJSCBDgQ0WB1geaoEthmBrQbxPFcspivPxMyDeb8oz4r7eSIGq4++sQzWMlo1WNvsy2Z1No+ABmvz2twaS6AiMDuI7wt8CdwLPNQBUxcz8EYxcDFxi46tyuyrr6/BGlNrF6YdmspTJSCBVRHQYK2KtNeRwPQIzA7ipwFvAwcCf3aQ28UMLGNa+uobw2Atq7UL0w5N5akSkMCqCGiwVkXa60hgegRmB/FjgE+BO4HHgG+Bf1vI7mIGljFYffWNYbCW1dqFaYsm8hQJSGDVBDRYqybu9SQwHQJbDeLXA/cB+wH/AccCnxXJ5wPHl/fPBs4py4ldzMAyBisyFumra7oYOA+4aQb3VLRWLL8aQOd0epRKJCCBPQQ0WHYGCWwugVmzcTOQ16NAjNAPwAfFaC2itMi03ArkVR37A38Bf9f+dwHwVgsj1EXfbcDXJd6r7QzWGFr76NzcHmvNJbBGBDRYa9RYSpXAwATqxugs4HXgBOCTOddJuoadwMfAs2Up8X1gkcE6FMirOp4quxOfqf0vS5G/NxisNvrqmp4G7gYyQ5Rr5r0se05Fa8Xyrh46B+4GFicBCewNAhqsvUHVMiWwHgTqZiMDfXJdnbJA+udlyTCzTzFhJwJ/NJiW2eL6LhG20VfX9CFwCfA4sAt4tQhZxRJhG60VyxjUZXWuR29TpQQ2jIAGa8Ma3OpKoEagbjauAR4GbgF2A98Bv9TOPQR4Ezip7DJ8FzhuhaalSV92PlaaDgC+KUuESZr60Zw6N3WGvmawSWvF8syBdDbVw/clIIERCGiwRoDuJSUwEQJ1g7WjmKsYkqNKkPsTwBVFa5bobihZ3k8vweOXrtBgNemra0qKhDuAI8rsUBWkH7mrmMFq0lqxfGAgnRPpTsqQgATqBDRY9gcJbC6BRWbjIuA5YB/gH+C6EkuVJKRXA4cD9/QwWG1pNxmhWX11Tfn7MODlErB/LvDbhLRWLH9csc627D1PAhIYgIAGawCIFiGBNSUwz8Rkue12II96ObXU7UHgFeAFIDMvrwHPj2RattJX15S/E3MVrTEzqedlE9JasUyqhlXqXNNuqmwJrCcBDdZ6tpuqJTAEgXkG61rgSuCqOTsK3wMuBL4fybRspW9W0zw+TTNjfbguKrOJ5Sp19qmbn5GABHoS0GD1BOfHJLANCBwJJCHn/S3rkmcBvgO8NJPbqms5bS7Xtsx5muZdo225bTRW56xLmV3q5LkSkMCSBDRYSwL04xKQgAQkIAEJSGCWgAbLPiEBCUhAAhKQgAQGJqDBGhioxUlAAhKQgAQkIAENln1AAhKQgAQkIAEJDEzgf3zQxu5YovBRAAAAAElFTkSuQmCC\" width=\"600\">"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "%matplotlib notebook\n",
    "trans_data=[('R,I,G', 1.0, 162.20894646124296, 0.006164888076866542), ('R,I,ZTS', 2.0, 198.82243025498838, 0.010059227208092237), ('R,G,ZTS', 16.0, 173.95017730614614, 0.09198036039848678), ('I,G,ZTS', 13.0, 84.7907192881195, 0.1533186663486826)]\n",
    "trans_data.sort(key=lambda x:x[2])\n",
    "import matplotlib.pyplot as plt\n",
    "fig=plt.figure()\n",
    "fig.set_size_inches(6,4.5)\n",
    "ax1=fig.subplots(1)\n",
    "xs=[i[2] for i in trans_data]\n",
    "ys=[i[1] for i in trans_data]\n",
    "ax1.plot(xs,ys,\"X-\",ms=3)\n",
    "ax1.plot([0,max(xs)],[0,max(xs)],\"--\")\n",
    "ax1.axhline(y=0, xmin=0, xmax=1,c=\"green\")\n",
    "ax1.axvline(x=0, ymin=0, ymax=1,c=\"green\")\n",
    "ax1.set_xlabel(\"$|\\\\xi_{ij}|+|\\\\xi_{jk}|+|\\\\xi_{ik}|$\")\n",
    "ax1.set_ylabel(\"$|\\\\xi_{ij}+\\\\xi_{jk}-\\\\xi_{ik}|$\")\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.6.9"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}