WaveCam_SV_noCV2.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Simulate the propagation of sound and its interaction with boundaries in 2D\n",
    "#### The acoustic wave equation is solved using the finite-difference time-domain method to solve for pressure and velocity"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "scrolled": false
   },
   "outputs": [],
   "source": [
    "#### Jupyter notebook magic command to display Matplotlib plots\n",
    "%matplotlib notebook"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "scrolled": false
   },
   "outputs": [],
   "source": [
    "#### Importing the relevant packages\n",
    "import numpy as np\n",
    "import cv2\n",
    "import matplotlib.pyplot as plt\n",
    "import sys\n",
    "import pyaudio"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "scrolled": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Enter 'd' to run the default setup or 'c' for customizing: d\n",
      "\n",
      "Point source of 15000 Hz in a closed domain with air at 295 K using the default camera resolution\n",
      "\n"
     ]
    }
   ],
   "source": [
    "#### Ask for user input\n",
    "plt.ion()\n",
    "dflag = input(\"Enter 'd' to run the default setup or 'c' for customizing: \")\n",
    "print(\"\")\n",
    "if dflag == \"d\":\n",
    "    freq = 15000\n",
    "    nm = 2\n",
    "    c0 = (346.13, 0)\n",
    "    rho = (1.2, 1.0e6)\n",
    "    stype = \"point\"\n",
    "    mflag = \"n\"\n",
    "    sflag = \"n\"\n",
    "    print(\"Point source of 15000 Hz in a closed domain with air at 295 K using the default camera resolution\")\n",
    "    print(\"\")\n",
    "elif dflag == \"c\":\n",
    "    sflag = input(\"Enter 'a' to record an audio source or 'n' to enter a numeric source frequency: \")\n",
    "    print(\"\")\n",
    "    if sflag == \"a\":\n",
    "        # Initialize portaudio\n",
    "        p = pyaudio.PyAudio()\n",
    "        numdev = p.get_device_count() - 1\n",
    "        if (numdev < 1):\n",
    "            sys.exit(\"Error: You do not have the hardware to record audio\")\n",
    "        else:\n",
    "            freq = 0\n",
    "            rt = 80000\n",
    "            rec_sec = 1\n",
    "            chsize = 1024\n",
    "            input(\"Start the audio and then press Enter to record...\")\n",
    "            print(\"\")\n",
    "            astream = p.open(format=pyaudio.paInt16, channels=1, rate=rt, input=True, frames_per_buffer=chsize)\n",
    "            # Record audio\n",
    "            frames = []\n",
    "            for i in range(0, np.int(rt / chsize * rec_sec)):\n",
    "                adata = astream.read(chsize)\n",
    "                frames.append(np.fromstring(adata, dtype=np.int16))\n",
    "            # Close audio stream\n",
    "            astream.stop_stream()\n",
    "            astream.close()\n",
    "            p.terminate()\n",
    "            atemp = np.asfarray(np.hstack(frames))\n",
    "            skipf = int(np.alen(atemp) / 4)\n",
    "            last = int(np.alen(atemp))\n",
    "            asource = ((atemp[skipf:last - skipf] - np.mean(atemp[skipf:last - skipf]))\n",
    "                       / (np.amax(atemp[skipf:last - skipf]) - np.amin(atemp[skipf:last - skipf])))\n",
    "            sfreqspec = np.fft.rfft(asource)\n",
    "            time = np.linspace(0, last, last) * rec_sec / rt\n",
    "            fig1 = plt.figure()\n",
    "            ax1 = fig1.add_subplot(111)\n",
    "            ax1.set_xlabel(\"time (s)\")\n",
    "            ax1.set_ylabel(\"rescaled amplitude (arb units)\")\n",
    "            ax1.set_title(\"Recorded audio\")\n",
    "            ax1.plot(time[skipf:last - skipf], asource)\n",
    "            fig1.tight_layout()\n",
    "            fig1.show()\n",
    "            fig1.canvas.draw()\n",
    "            fig2 = plt.figure()\n",
    "            ax2 = fig2.add_subplot(111)\n",
    "            ax2.set_xlabel('freq (Hz)')\n",
    "            ax2.set_ylabel('power (arb units)')\n",
    "            ax2.set_title(\"Power spectrum of recorded audio\")\n",
    "            ax2.plot(np.abs(sfreqspec)**2)\n",
    "            fig2.tight_layout()\n",
    "            fig2.show()\n",
    "            fig2.canvas.draw()\n",
    "    elif sflag == \"n\":\n",
    "        try:\n",
    "            freq = float(input(\"Enter the sound frequency in Hz (20-20000): \"))\n",
    "            print(\"\")\n",
    "        except ValueError:\n",
    "            sys.exit(\"Error: enter a number between 20 and 20000\")\n",
    "    else:\n",
    "        sys.exit(\"Error: press 'a' to record an audio source or n to enter a numeric source frequency\")\n",
    "    vs = input(\"Enter sound velocity in m/s (If medium is air or water type 'air' or 'water'): \")\n",
    "    print(\"\")\n",
    "    if vs == \"air\":\n",
    "        c0 = (346.13, 0)\n",
    "        rho = (1.2, 1.0e6)\n",
    "    elif vs == \"water\":\n",
    "        c0 = (1481, 0)\n",
    "        rho = (1000, 1.0e6)\n",
    "    else:\n",
    "        try:\n",
    "            c0 = (float(vs), 0)\n",
    "            mdensity = float(input(\"Enter density of medium in kg/m^3: \"))\n",
    "            print(\"\")\n",
    "            rho = (mdensity, 1.0e6)\n",
    "        except ValueError:\n",
    "            sys.exit(\"Error: enter a numeric value, or 'air', or 'water' \")\n",
    "    stype = input(\"Enter 'point' or 'line' for type of source: \")\n",
    "    print(\"\")\n",
    "    if stype != \"line\" and stype != \"point\":\n",
    "        sys.exit(\"Error: type either point or line for type of source\")\n",
    "    mflag = input(\"Do you want to insert a block of medium at a different temperature? (y/n): \")\n",
    "    print(\"\")\n",
    "    if mflag == \"y\":\n",
    "        nm = 3\n",
    "        temparature = float(input(\"Enter the absolute temperature of the block in K (50-500): \"))\n",
    "        print(\"\")\n",
    "        ct = c0[0] * np.sqrt(temparature/293)\n",
    "        c0 = (c0[0], 0, ct)\n",
    "        rho = (rho[0], 1.0e6, rho[0])\n",
    "    elif mflag == \"n\":\n",
    "        nm = 2\n",
    "    else:\n",
    "        sys.exit(\"Error: enter 'y' to insert a block of different temperature or enter 'n' \")\n",
    "else:\n",
    "    sys.exit(\"Error: enter 'd' for default setup or 'c' to customize\")\n",
    "\n",
    "if np.amin(c0) == 0:\n",
    "    wavelmin = 300 / 20000.0\n",
    "else:\n",
    "    wavelmin = np.amin(c0) / 20000.0"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {
    "scrolled": false
   },
   "outputs": [],
   "source": [
    "#### The main class that defines all constants, variables, and functions\n",
    "class fdtdVar:\n",
    "    def __init__(self, rs, cs):\n",
    "        ## Constants\n",
    "        cn = 0.9 / np.sqrt(2.0)  # Courant number\n",
    "        ## Variables\n",
    "        self.r = np.int(rs)  # number of rows\n",
    "        self.c = np.int(cs)  # number of columns\n",
    "        self.freq = freq  # frequency of source\n",
    "        temp = (self.r, self.c + 1)\n",
    "        self.vx = np.zeros(temp)  # velocity along x\n",
    "        self.mvx = np.zeros(temp, dtype=np.int8)\n",
    "        temp = (self.r + 1, self.c)\n",
    "        self.vy = np.zeros(temp)  # velocity along y\n",
    "        self.mvy = np.zeros(temp, dtype=np.int8)\n",
    "        temp = (self.r, self.c)\n",
    "        self.pr = np.zeros(temp)  # pressure\n",
    "        self.gaussamp = np.zeros(temp)\n",
    "        self.mpr = np.zeros(temp, dtype=np.int8)\n",
    "        self.dx = wavelmin/10.0  # grid cell size\n",
    "        self.dt = cn * self.dx / np.amax(c0)  # time step size\n",
    "        self.ca = np.ones(nm)\n",
    "        self.cb = np.ones(nm)\n",
    "        self.da = np.ones(nm)\n",
    "        self.db = np.ones(nm)\n",
    "        for i in range(0, nm, 1):\n",
    "            self.cb[i] = c0[i] ** 2 * rho[i] * self.dt / self.dx\n",
    "            self.db[i] = self.dt / (rho[i] * self.dx)\n",
    "        self.da[1] = 0\n",
    "        temp = (self.r, 2, 2)\n",
    "        self.vxl = np.zeros(temp)\n",
    "        self.vxr = np.zeros(temp)\n",
    "        temp = (self.c, 2, 2)\n",
    "        self.vyb = np.zeros(temp)\n",
    "        self.vyt = np.zeros(temp)\n",
    "        print(\"dx [m] = \", self.dx)\n",
    "        print(\"dt [s] = \", self.dt)\n",
    "        print(\"\")\n",
    "        rtemp = np.arange(0, self.r, 1)\n",
    "        ctemp = np.arange(0, self.c, 1)\n",
    "        rm, cm = np.meshgrid(rtemp, ctemp)\n",
    "        rc = np.int(self.r / 2)\n",
    "        cc = np.int(self.c / 2)\n",
    "        if stype == \"point\":\n",
    "            fwhmc = 2\n",
    "            fwhmr = fwhmc\n",
    "            self.gaussamp = np.exp(-((rm - rc) ** 2 / (2 * fwhmr ** 2) + (cm - cc) ** 2 / (2 * fwhmc ** 2))).T\n",
    "        elif stype == \"line\":\n",
    "            fwhmc = 2\n",
    "            fwhmr = 16\n",
    "            self.gaussamp = np.exp(-((rm - rc) ** 2 / (2 * fwhmr ** 2) + (cm - cc) ** 2 / (2 * fwhmc ** 2))).T\n",
    "    \n",
    "    def source(self, nt):\n",
    "        rm = self.r\n",
    "        cm = self.c\n",
    "        # prs = self.dx * np.sin(2 * pi * self.freq * nt * self.dt) / self.cb[0]\n",
    "        if sflag == \"a\":\n",
    "            if nt < np.alen(asource):\n",
    "                prs = self.dx * asource[nt] / self.cb[0]\n",
    "            else:\n",
    "                prs = 0\n",
    "        else:\n",
    "            prs = self.dx * np.sin(2 * np.pi * self.freq * nt * self.dt) / self.cb[0]\n",
    "\n",
    "        # Update pressure with source\n",
    "        self.pr[1:rm - 1, 1:cm - 1] = (self.pr[1:rm - 1, 1:cm - 1]\n",
    "                                       - self.cb[self.mpr[1:rm - 1, 1:cm - 1]] * prs\n",
    "                                       * self.gaussamp[1:rm - 1, 1:cm - 1] / self.dx)\n",
    "\n",
    "    def fdtd_update(self):\n",
    "        ri = self.r\n",
    "        ci = self.c\n",
    "        self.pr[0:ri, 0:ci] = (self.ca[self.mpr[0:ri, 0:ci]] * self.pr[0:ri, 0:ci]\n",
    "                               - self.cb[self.mpr[0:ri, 0:ci]]\n",
    "                               * ((self.vx[0:ri, 1:ci + 1] - self.vx[0:ri, 0:ci])\n",
    "                                  + (self.vy[1:ri + 1, 0:ci] - self.vy[0:ri, 0:ci])))\n",
    "        self.vx[0:ri, 1:ci] = (self.da[self.mvx[0:ri, 1:ci]] * self.vx[0:ri, 1:ci]\n",
    "                               - self.db[self.mvx[0:ri, 1:ci]] * (self.pr[0:ri, 1:ci] - self.pr[0:ri, 0:ci - 1]))\n",
    "        self.vy[1:ri, 0:ci] = (self.da[self.mvy[1:ri, 0:ci]] * self.vy[1:ri, 0:ci]\n",
    "                               - self.db[self.mvy[1:ri, 0:ci]] * (self.pr[1:ri, 0:ci] - self.pr[0:ri - 1, 0:ci]))\n",
    "\n",
    "    def boundary(self):\n",
    "        ri = self.r\n",
    "        ci = self.c\n",
    "        c1 = (c0[0] * self.dt - self.dx) / (c0[0] * self.dt + self.dx)\n",
    "        c2 = 2 * self.dx / (c0[0] * self.dt + self.dx)\n",
    "        c3 = (c0[0] * self.dt) ** 2 / (2 * self.dx * (c0[0] * self.dt + self.dx))\n",
    "        # Left and right boundaries\n",
    "        self.vx[1:ri - 1, 0] = (-self.vxl[1:ri - 1, 1, 1]\n",
    "                                + c1 * (self.vx[1:ri - 1, 1] + self.vxl[1:ri - 1, 0, 1])\n",
    "                                + c2 * (self.vxl[1:ri - 1, 0, 0] + self.vxl[1:ri - 1, 1, 0])\n",
    "                                + c3 * (self.vxl[2:ri, 0, 0] - 2 * self.vxl[1:ri - 1, 0, 0]\n",
    "                                        + self.vxl[0:ri - 2, 0, 0] + self.vxl[2:ri, 1, 0]\n",
    "                                        - 2 * self.vxl[1:ri - 1, 1, 0] + self.vxl[0:ri - 2, 1, 0]))\n",
    "        self.vx[1:ri - 1, ci] = (-self.vxr[1:ri - 1, 1, 1]\n",
    "                                 + c1 * (self.vx[1:ri - 1, ci - 1] + self.vxr[1:ri - 1, 0, 1])\n",
    "                                 + c2 * (self.vxr[1:ri - 1, 0, 0] + self.vxr[1:ri - 1, 1, 0])\n",
    "                                 + c3 * (self.vxr[2:ri, 0, 0] - 2 * self.vxr[1:ri - 1, 0, 0]\n",
    "                                         + self.vxr[0:ri - 2, 0, 0] + self.vxr[2:ri, 1, 0]\n",
    "                                         - 2 * self.vxr[1:ri - 1, 1, 0] + self.vxr[0:ri - 2, 1, 0]))\n",
    "\n",
    "        # Bottom and top boundaries\n",
    "        self.vy[0, 1:ci - 1] = (-self.vyb[1:ci - 1, 1, 1]\n",
    "                                + c1 * (self.vy[1, 1:ci - 1] + self.vyb[1:ci - 1, 0, 1])\n",
    "                                + c2 * (self.vyb[1:ci - 1, 0, 0] + self.vyb[1:ci - 1, 1, 0])\n",
    "                                + c3 * (self.vyb[2:ci, 0, 0] - 2 * self.vyb[1:ci - 1, 0, 0]\n",
    "                                        + self.vyb[0:ci - 2, 0, 0] + self.vyb[2:ci, 1, 0]\n",
    "                                        - 2 * self.vyb[1:ci - 1, 1, 0] + self.vyb[0:ci - 2, 1, 0]))\n",
    "        self.vy[ri, 1:ci - 1] = (-self.vyt[1:ci - 1, 1, 1]\n",
    "                                 + c1 * (self.vy[ri - 1, 1:ci - 1] + self.vyt[1:ci - 1, 0, 1])\n",
    "                                 + c2 * (self.vyt[1:ci - 1, 0, 0] + self.vyt[1:ci - 1, 1, 0])\n",
    "                                 + c3 * (self.vyt[2:ci, 0, 0] - 2 * self.vyt[1:ci - 1, 0, 0]\n",
    "                                         + self.vyt[0:ci - 2, 0, 0] + self.vyt[2:ci, 1, 0]\n",
    "                                         - 2 * self.vyt[1:ci - 1, 1, 0] + self.vyt[0:ci - 2, 1, 0]))\n",
    "        # Corners\n",
    "        self.vx[0, 0] = self.vxl[1, 1, 1]\n",
    "        self.vx[ri - 1, 0] = self.vxl[ri - 2, 1, 1]\n",
    "        self.vx[0, ci] = self.vxr[1, 1, 1]\n",
    "        self.vx[ri - 1, ci] = self.vxr[ri - 2, 1, 1]\n",
    "        self.vy[0, 0] = self.vyb[1, 1, 1]\n",
    "        self.vy[0, ci - 1] = self.vyb[ci - 2, 1, 1]\n",
    "        self.vy[ri, 0] = self.vyt[1, 1, 1]\n",
    "        self.vy[ri, ci - 1] = self.vyt[ci - 2, 1, 1]\n",
    "\n",
    "        # Store boundary values\n",
    "        for i in range(0, 2, 1):\n",
    "            self.vxl[0:ri, i, 1] = self.vxl[0:ri, i, 0]\n",
    "            self.vxl[0:ri, i, 0] = self.vx[0:ri, i]\n",
    "            self.vxr[0:ri, i, 1] = self.vxr[0:ri, i, 0]\n",
    "            self.vxr[0:ri, i, 0] = self.vx[0:ri, ci - i]\n",
    "            self.vyb[0:ci, i, 1] = self.vyb[0:ci, i, 0]\n",
    "            self.vyb[0:ci, i, 0] = self.vy[i, 0:ci]\n",
    "            self.vyt[0:ci, i, 1] = self.vyt[0:ci, i, 0]\n",
    "            self.vyt[0:ci, i, 0] = self.vy[ri - i, 0:ci]\n",
    "\n",
    "    def update_domain(self):\n",
    "        if mflag == \"y\":\n",
    "            cm = self.c\n",
    "            rm = self.r\n",
    "            c1 = np.int(cm/2) + np.int(cm/8)\n",
    "            c2 = c1 + np.int(cm/8)\n",
    "            self.mvx.fill(0)\n",
    "            self.mvy.fill(0)\n",
    "            self.mpr.fill(0)\n",
    "            self.mvx[40:rm - 40, c1:c2] = 2\n",
    "            self.mvy[40:rm - 40, c1:c2] = 2\n",
    "            self.mpr[40:rm - 40, c1:c2] = 2\n",
    "        elif mflag == \"n\":\n",
    "            self.mvx.fill(0)\n",
    "            self.mvy.fill(0)\n",
    "            self.mpr.fill(0)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {
    "scrolled": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Do you want to read an image (enter i) or capture a frame from webcam (enter w)? i\n",
      "\n",
      "Enter the filename of your image including extension: test1.jpg\n",
      "\n",
      "Resized the image to 640 x 480 pixels\n",
      "\n"
     ]
    }
   ],
   "source": [
    "#### Ask user for option to read an image from file or start webcam to capture a frame\n",
    "\n",
    "vflag = input(\"Do you want to read an image (enter i) or capture a frame from webcam (enter w)? \")\n",
    "print(\"\")\n",
    "if vflag == \"w\":\n",
    "    ## Create VideoCapture object\n",
    "    cap = cv2.VideoCapture(0)\n",
    "\n",
    "    ## Array sizes (resolution dependent)\n",
    "    retval = cap.get(3)\n",
    "    columns = retval\n",
    "    retval = cap.get(4)\n",
    "    rows = retval\n",
    "\n",
    "    ## Set frame resolution\n",
    "    if dflag == \"c\" and columns > 1000:\n",
    "        print(\"Default frame resolution \", np.int(columns), \"x\", np.int(rows))\n",
    "        if input(\"Do you want to reduce resolution to increase speed? (y/n): \") == \"y\":\n",
    "            print(\"\")\n",
    "            print(\"Setting new resolution to 640 x 480\")\n",
    "            print(\"\")\n",
    "            columns = 640\n",
    "            rows = 480\n",
    "\n",
    "    cap.set(3, columns)\n",
    "    cap.set(4, rows)\n",
    "elif vflag == \"i\":\n",
    "    ## Read image\n",
    "    imgname = input(\"Enter the filename of your image including extension: \")\n",
    "    print(\"\")\n",
    "    img = cv2.imread(imgname, 0)\n",
    "    rows, columns = img.shape\n",
    "    ar = 640 / columns\n",
    "    dim = (640, int(rows * ar))\n",
    "    img = cv2.medianBlur(img, 21)\n",
    "    resimg = cv2.resize(img, dim, interpolation = cv2.INTER_AREA)\n",
    "    rows, columns = resimg.shape\n",
    "    print(\"Resized the image to 640 x 480 pixels\")\n",
    "    print(\"\")\n",
    "    threshimg = cv2.adaptiveThreshold(resimg, 255, cv2.ADAPTIVE_THRESH_GAUSSIAN_C, cv2.THRESH_BINARY, 17, 5)\n",
    "else:\n",
    "    sys.exit(\"Error: type either 'i' to read image or 'w' to read from webcam\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {
    "scrolled": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "dx [m] =  0.0015\n",
      "dt [s] =  2.75790643574e-06\n",
      "\n"
     ]
    },
    {
     "data": {
      "application/javascript": [
       "/* Put everything inside the global mpl namespace */\n",
       "window.mpl = {};\n",
       "\n",
       "\n",
       "mpl.get_websocket_type = function() {\n",
       "    if (typeof(WebSocket) !== 'undefined') {\n",
       "        return WebSocket;\n",
       "    } else if (typeof(MozWebSocket) !== 'undefined') {\n",
       "        return MozWebSocket;\n",
       "    } else {\n",
       "        alert('Your browser does not have WebSocket support.' +\n",
       "              'Please try Chrome, Safari or Firefox ≥ 6. ' +\n",
       "              'Firefox 4 and 5 are also supported but you ' +\n",
       "              'have to enable WebSockets in about:config.');\n",
       "    };\n",
       "}\n",
       "\n",
       "mpl.figure = function(figure_id, websocket, ondownload, parent_element) {\n",
       "    this.id = figure_id;\n",
       "\n",
       "    this.ws = websocket;\n",
       "\n",
       "    this.supports_binary = (this.ws.binaryType != undefined);\n",
       "\n",
       "    if (!this.supports_binary) {\n",
       "        var warnings = document.getElementById(\"mpl-warnings\");\n",
       "        if (warnings) {\n",
       "            warnings.style.display = 'block';\n",
       "            warnings.textContent = (\n",
       "                \"This browser does not support binary websocket messages. \" +\n",
       "                    \"Performance may be slow.\");\n",
       "        }\n",
       "    }\n",
       "\n",
       "    this.imageObj = new Image();\n",
       "\n",
       "    this.context = undefined;\n",
       "    this.message = undefined;\n",
       "    this.canvas = undefined;\n",
       "    this.rubberband_canvas = undefined;\n",
       "    this.rubberband_context = undefined;\n",
       "    this.format_dropdown = undefined;\n",
       "\n",
       "    this.image_mode = 'full';\n",
       "\n",
       "    this.root = $('<div/>');\n",
       "    this._root_extra_style(this.root)\n",
       "    this.root.attr('style', 'display: inline-block');\n",
       "\n",
       "    $(parent_element).append(this.root);\n",
       "\n",
       "    this._init_header(this);\n",
       "    this._init_canvas(this);\n",
       "    this._init_toolbar(this);\n",
       "\n",
       "    var fig = this;\n",
       "\n",
       "    this.waiting = false;\n",
       "\n",
       "    this.ws.onopen =  function () {\n",
       "            fig.send_message(\"supports_binary\", {value: fig.supports_binary});\n",
       "            fig.send_message(\"send_image_mode\", {});\n",
       "            if (mpl.ratio != 1) {\n",
       "                fig.send_message(\"set_dpi_ratio\", {'dpi_ratio': mpl.ratio});\n",
       "            }\n",
       "            fig.send_message(\"refresh\", {});\n",
       "        }\n",
       "\n",
       "    this.imageObj.onload = function() {\n",
       "            if (fig.image_mode == 'full') {\n",
       "                // Full images could contain transparency (where diff images\n",
       "                // almost always do), so we need to clear the canvas so that\n",
       "                // there is no ghosting.\n",
       "                fig.context.clearRect(0, 0, fig.canvas.width, fig.canvas.height);\n",
       "            }\n",
       "            fig.context.drawImage(fig.imageObj, 0, 0);\n",
       "        };\n",
       "\n",
       "    this.imageObj.onunload = function() {\n",
       "        fig.ws.close();\n",
       "    }\n",
       "\n",
       "    this.ws.onmessage = this._make_on_message_function(this);\n",
       "\n",
       "    this.ondownload = ondownload;\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._init_header = function() {\n",
       "    var titlebar = $(\n",
       "        '<div class=\"ui-dialog-titlebar ui-widget-header ui-corner-all ' +\n",
       "        'ui-helper-clearfix\"/>');\n",
       "    var titletext = $(\n",
       "        '<div class=\"ui-dialog-title\" style=\"width: 100%; ' +\n",
       "        'text-align: center; padding: 3px;\"/>');\n",
       "    titlebar.append(titletext)\n",
       "    this.root.append(titlebar);\n",
       "    this.header = titletext[0];\n",
       "}\n",
       "\n",
       "\n",
       "\n",
       "mpl.figure.prototype._canvas_extra_style = function(canvas_div) {\n",
       "\n",
       "}\n",
       "\n",
       "\n",
       "mpl.figure.prototype._root_extra_style = function(canvas_div) {\n",
       "\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._init_canvas = function() {\n",
       "    var fig = this;\n",
       "\n",
       "    var canvas_div = $('<div/>');\n",
       "\n",
       "    canvas_div.attr('style', 'position: relative; clear: both; outline: 0');\n",
       "\n",
       "    function canvas_keyboard_event(event) {\n",
       "        return fig.key_event(event, event['data']);\n",
       "    }\n",
       "\n",
       "    canvas_div.keydown('key_press', canvas_keyboard_event);\n",
       "    canvas_div.keyup('key_release', canvas_keyboard_event);\n",
       "    this.canvas_div = canvas_div\n",
       "    this._canvas_extra_style(canvas_div)\n",
       "    this.root.append(canvas_div);\n",
       "\n",
       "    var canvas = $('<canvas/>');\n",
       "    canvas.addClass('mpl-canvas');\n",
       "    canvas.attr('style', \"left: 0; top: 0; z-index: 0; outline: 0\")\n",
       "\n",
       "    this.canvas = canvas[0];\n",
       "    this.context = canvas[0].getContext(\"2d\");\n",
       "\n",
       "    var backingStore = this.context.backingStorePixelRatio ||\n",
       "\tthis.context.webkitBackingStorePixelRatio ||\n",
       "\tthis.context.mozBackingStorePixelRatio ||\n",
       "\tthis.context.msBackingStorePixelRatio ||\n",
       "\tthis.context.oBackingStorePixelRatio ||\n",
       "\tthis.context.backingStorePixelRatio || 1;\n",
       "\n",
       "    mpl.ratio = (window.devicePixelRatio || 1) / backingStore;\n",
       "\n",
       "    var rubberband = $('<canvas/>');\n",
       "    rubberband.attr('style', \"position: absolute; left: 0; top: 0; z-index: 1;\")\n",
       "\n",
       "    var pass_mouse_events = true;\n",
       "\n",
       "    canvas_div.resizable({\n",
       "        start: function(event, ui) {\n",
       "            pass_mouse_events = false;\n",
       "        },\n",
       "        resize: function(event, ui) {\n",
       "            fig.request_resize(ui.size.width, ui.size.height);\n",
       "        },\n",
       "        stop: function(event, ui) {\n",
       "            pass_mouse_events = true;\n",
       "            fig.request_resize(ui.size.width, ui.size.height);\n",
       "        },\n",
       "    });\n",
       "\n",
       "    function mouse_event_fn(event) {\n",
       "        if (pass_mouse_events)\n",
       "            return fig.mouse_event(event, event['data']);\n",
       "    }\n",
       "\n",
       "    rubberband.mousedown('button_press', mouse_event_fn);\n",
       "    rubberband.mouseup('button_release', mouse_event_fn);\n",
       "    // Throttle sequential mouse events to 1 every 20ms.\n",
       "    rubberband.mousemove('motion_notify', mouse_event_fn);\n",
       "\n",
       "    rubberband.mouseenter('figure_enter', mouse_event_fn);\n",
       "    rubberband.mouseleave('figure_leave', mouse_event_fn);\n",
       "\n",
       "    canvas_div.on(\"wheel\", function (event) {\n",
       "        event = event.originalEvent;\n",
       "        event['data'] = 'scroll'\n",
       "        if (event.deltaY < 0) {\n",
       "            event.step = 1;\n",
       "        } else {\n",
       "            event.step = -1;\n",
       "        }\n",
       "        mouse_event_fn(event);\n",
       "    });\n",
       "\n",
       "    canvas_div.append(canvas);\n",
       "    canvas_div.append(rubberband);\n",
       "\n",
       "    this.rubberband = rubberband;\n",
       "    this.rubberband_canvas = rubberband[0];\n",
       "    this.rubberband_context = rubberband[0].getContext(\"2d\");\n",
       "    this.rubberband_context.strokeStyle = \"#000000\";\n",
       "\n",
       "    this._resize_canvas = function(width, height) {\n",
       "        // Keep the size of the canvas, canvas container, and rubber band\n",
       "        // canvas in synch.\n",
       "        canvas_div.css('width', width)\n",
       "        canvas_div.css('height', height)\n",
       "\n",
       "        canvas.attr('width', width * mpl.ratio);\n",
       "        canvas.attr('height', height * mpl.ratio);\n",
       "        canvas.attr('style', 'width: ' + width + 'px; height: ' + height + 'px;');\n",
       "\n",
       "        rubberband.attr('width', width);\n",
       "        rubberband.attr('height', height);\n",
       "    }\n",
       "\n",
       "    // Set the figure to an initial 600x600px, this will subsequently be updated\n",
       "    // upon first draw.\n",
       "    this._resize_canvas(600, 600);\n",
       "\n",
       "    // Disable right mouse context menu.\n",
       "    $(this.rubberband_canvas).bind(\"contextmenu\",function(e){\n",
       "        return false;\n",
       "    });\n",
       "\n",
       "    function set_focus () {\n",
       "        canvas.focus();\n",
       "        canvas_div.focus();\n",
       "    }\n",
       "\n",
       "    window.setTimeout(set_focus, 100);\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._init_toolbar = function() {\n",
       "    var fig = this;\n",
       "\n",
       "    var nav_element = $('<div/>')\n",
       "    nav_element.attr('style', 'width: 100%');\n",
       "    this.root.append(nav_element);\n",
       "\n",
       "    // Define a callback function for later on.\n",
       "    function toolbar_event(event) {\n",
       "        return fig.toolbar_button_onclick(event['data']);\n",
       "    }\n",
       "    function toolbar_mouse_event(event) {\n",
       "        return fig.toolbar_button_onmouseover(event['data']);\n",
       "    }\n",
       "\n",
       "    for(var toolbar_ind in mpl.toolbar_items) {\n",
       "        var name = mpl.toolbar_items[toolbar_ind][0];\n",
       "        var tooltip = mpl.toolbar_items[toolbar_ind][1];\n",
       "        var image = mpl.toolbar_items[toolbar_ind][2];\n",
       "        var method_name = mpl.toolbar_items[toolbar_ind][3];\n",
       "\n",
       "        if (!name) {\n",
       "            // put a spacer in here.\n",
       "            continue;\n",
       "        }\n",
       "        var button = $('<button/>');\n",
       "        button.addClass('ui-button ui-widget ui-state-default ui-corner-all ' +\n",
       "                        'ui-button-icon-only');\n",
       "        button.attr('role', 'button');\n",
       "        button.attr('aria-disabled', 'false');\n",
       "        button.click(method_name, toolbar_event);\n",
       "        button.mouseover(tooltip, toolbar_mouse_event);\n",
       "\n",
       "        var icon_img = $('<span/>');\n",
       "        icon_img.addClass('ui-button-icon-primary ui-icon');\n",
       "        icon_img.addClass(image);\n",
       "        icon_img.addClass('ui-corner-all');\n",
       "\n",
       "        var tooltip_span = $('<span/>');\n",
       "        tooltip_span.addClass('ui-button-text');\n",
       "        tooltip_span.html(tooltip);\n",
       "\n",
       "        button.append(icon_img);\n",
       "        button.append(tooltip_span);\n",
       "\n",
       "        nav_element.append(button);\n",
       "    }\n",
       "\n",
       "    var fmt_picker_span = $('<span/>');\n",
       "\n",
       "    var fmt_picker = $('<select/>');\n",
       "    fmt_picker.addClass('mpl-toolbar-option ui-widget ui-widget-content');\n",
       "    fmt_picker_span.append(fmt_picker);\n",
       "    nav_element.append(fmt_picker_span);\n",
       "    this.format_dropdown = fmt_picker[0];\n",
       "\n",
       "    for (var ind in mpl.extensions) {\n",
       "        var fmt = mpl.extensions[ind];\n",
       "        var option = $(\n",
       "            '<option/>', {selected: fmt === mpl.default_extension}).html(fmt);\n",
       "        fmt_picker.append(option)\n",
       "    }\n",
       "\n",
       "    // Add hover states to the ui-buttons\n",
       "    $( \".ui-button\" ).hover(\n",
       "        function() { $(this).addClass(\"ui-state-hover\");},\n",
       "        function() { $(this).removeClass(\"ui-state-hover\");}\n",
       "    );\n",
       "\n",
       "    var status_bar = $('<span class=\"mpl-message\"/>');\n",
       "    nav_element.append(status_bar);\n",
       "    this.message = status_bar[0];\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.request_resize = function(x_pixels, y_pixels) {\n",
       "    // Request matplotlib to resize the figure. Matplotlib will then trigger a resize in the client,\n",
       "    // which will in turn request a refresh of the image.\n",
       "    this.send_message('resize', {'width': x_pixels, 'height': y_pixels});\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.send_message = function(type, properties) {\n",
       "    properties['type'] = type;\n",
       "    properties['figure_id'] = this.id;\n",
       "    this.ws.send(JSON.stringify(properties));\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.send_draw_message = function() {\n",
       "    if (!this.waiting) {\n",
       "        this.waiting = true;\n",
       "        this.ws.send(JSON.stringify({type: \"draw\", figure_id: this.id}));\n",
       "    }\n",
       "}\n",
       "\n",
       "\n",
       "mpl.figure.prototype.handle_save = function(fig, msg) {\n",
       "    var format_dropdown = fig.format_dropdown;\n",
       "    var format = format_dropdown.options[format_dropdown.selectedIndex].value;\n",
       "    fig.ondownload(fig, format);\n",
       "}\n",
       "\n",
       "\n",
       "mpl.figure.prototype.handle_resize = function(fig, msg) {\n",
       "    var size = msg['size'];\n",
       "    if (size[0] != fig.canvas.width || size[1] != fig.canvas.height) {\n",
       "        fig._resize_canvas(size[0], size[1]);\n",
       "        fig.send_message(\"refresh\", {});\n",
       "    };\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_rubberband = function(fig, msg) {\n",
       "    var x0 = msg['x0'] / mpl.ratio;\n",
       "    var y0 = (fig.canvas.height - msg['y0']) / mpl.ratio;\n",
       "    var x1 = msg['x1'] / mpl.ratio;\n",
       "    var y1 = (fig.canvas.height - msg['y1']) / mpl.ratio;\n",
       "    x0 = Math.floor(x0) + 0.5;\n",
       "    y0 = Math.floor(y0) + 0.5;\n",
       "    x1 = Math.floor(x1) + 0.5;\n",
       "    y1 = Math.floor(y1) + 0.5;\n",
       "    var min_x = Math.min(x0, x1);\n",
       "    var min_y = Math.min(y0, y1);\n",
       "    var width = Math.abs(x1 - x0);\n",
       "    var height = Math.abs(y1 - y0);\n",
       "\n",
       "    fig.rubberband_context.clearRect(\n",
       "        0, 0, fig.canvas.width, fig.canvas.height);\n",
       "\n",
       "    fig.rubberband_context.strokeRect(min_x, min_y, width, height);\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_figure_label = function(fig, msg) {\n",
       "    // Updates the figure title.\n",
       "    fig.header.textContent = msg['label'];\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_cursor = function(fig, msg) {\n",
       "    var cursor = msg['cursor'];\n",
       "    switch(cursor)\n",
       "    {\n",
       "    case 0:\n",
       "        cursor = 'pointer';\n",
       "        break;\n",
       "    case 1:\n",
       "        cursor = 'default';\n",
       "        break;\n",
       "    case 2:\n",
       "        cursor = 'crosshair';\n",
       "        break;\n",
       "    case 3:\n",
       "        cursor = 'move';\n",
       "        break;\n",
       "    }\n",
       "    fig.rubberband_canvas.style.cursor = cursor;\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_message = function(fig, msg) {\n",
       "    fig.message.textContent = msg['message'];\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_draw = function(fig, msg) {\n",
       "    // Request the server to send over a new figure.\n",
       "    fig.send_draw_message();\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.handle_image_mode = function(fig, msg) {\n",
       "    fig.image_mode = msg['mode'];\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.updated_canvas_event = function() {\n",
       "    // Called whenever the canvas gets updated.\n",
       "    this.send_message(\"ack\", {});\n",
       "}\n",
       "\n",
       "// A function to construct a web socket function for onmessage handling.\n",
       "// Called in the figure constructor.\n",
       "mpl.figure.prototype._make_on_message_function = function(fig) {\n",
       "    return function socket_on_message(evt) {\n",
       "        if (evt.data instanceof Blob) {\n",
       "            /* FIXME: We get \"Resource interpreted as Image but\n",
       "             * transferred with MIME type text/plain:\" errors on\n",
       "             * Chrome.  But how to set the MIME type?  It doesn't seem\n",
       "             * to be part of the websocket stream */\n",
       "            evt.data.type = \"image/png\";\n",
       "\n",
       "            /* Free the memory for the previous frames */\n",
       "            if (fig.imageObj.src) {\n",
       "                (window.URL || window.webkitURL).revokeObjectURL(\n",
       "                    fig.imageObj.src);\n",
       "            }\n",
       "\n",
       "            fig.imageObj.src = (window.URL || window.webkitURL).createObjectURL(\n",
       "                evt.data);\n",
       "            fig.updated_canvas_event();\n",
       "            fig.waiting = false;\n",
       "            return;\n",
       "        }\n",
       "        else if (typeof evt.data === 'string' && evt.data.slice(0, 21) == \"data:image/png;base64\") {\n",
       "            fig.imageObj.src = evt.data;\n",
       "            fig.updated_canvas_event();\n",
       "            fig.waiting = false;\n",
       "            return;\n",
       "        }\n",
       "\n",
       "        var msg = JSON.parse(evt.data);\n",
       "        var msg_type = msg['type'];\n",
       "\n",
       "        // Call the  \"handle_{type}\" callback, which takes\n",
       "        // the figure and JSON message as its only arguments.\n",
       "        try {\n",
       "            var callback = fig[\"handle_\" + msg_type];\n",
       "        } catch (e) {\n",
       "            console.log(\"No handler for the '\" + msg_type + \"' message type: \", msg);\n",
       "            return;\n",
       "        }\n",
       "\n",
       "        if (callback) {\n",
       "            try {\n",
       "                // console.log(\"Handling '\" + msg_type + \"' message: \", msg);\n",
       "                callback(fig, msg);\n",
       "            } catch (e) {\n",
       "                console.log(\"Exception inside the 'handler_\" + msg_type + \"' callback:\", e, e.stack, msg);\n",
       "            }\n",
       "        }\n",
       "    };\n",
       "}\n",
       "\n",
       "// from http://stackoverflow.com/questions/1114465/getting-mouse-location-in-canvas\n",
       "mpl.findpos = function(e) {\n",
       "    //this section is from http://www.quirksmode.org/js/events_properties.html\n",
       "    var targ;\n",
       "    if (!e)\n",
       "        e = window.event;\n",
       "    if (e.target)\n",
       "        targ = e.target;\n",
       "    else if (e.srcElement)\n",
       "        targ = e.srcElement;\n",
       "    if (targ.nodeType == 3) // defeat Safari bug\n",
       "        targ = targ.parentNode;\n",
       "\n",
       "    // jQuery normalizes the pageX and pageY\n",
       "    // pageX,Y are the mouse positions relative to the document\n",
       "    // offset() returns the position of the element relative to the document\n",
       "    var x = e.pageX - $(targ).offset().left;\n",
       "    var y = e.pageY - $(targ).offset().top;\n",
       "\n",
       "    return {\"x\": x, \"y\": y};\n",
       "};\n",
       "\n",
       "/*\n",
       " * return a copy of an object with only non-object keys\n",
       " * we need this to avoid circular references\n",
       " * http://stackoverflow.com/a/24161582/3208463\n",
       " */\n",
       "function simpleKeys (original) {\n",
       "  return Object.keys(original).reduce(function (obj, key) {\n",
       "    if (typeof original[key] !== 'object')\n",
       "        obj[key] = original[key]\n",
       "    return obj;\n",
       "  }, {});\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.mouse_event = function(event, name) {\n",
       "    var canvas_pos = mpl.findpos(event)\n",
       "\n",
       "    if (name === 'button_press')\n",
       "    {\n",
       "        this.canvas.focus();\n",
       "        this.canvas_div.focus();\n",
       "    }\n",
       "\n",
       "    var x = canvas_pos.x * mpl.ratio;\n",
       "    var y = canvas_pos.y * mpl.ratio;\n",
       "\n",
       "    this.send_message(name, {x: x, y: y, button: event.button,\n",
       "                             step: event.step,\n",
       "                             guiEvent: simpleKeys(event)});\n",
       "\n",
       "    /* This prevents the web browser from automatically changing to\n",
       "     * the text insertion cursor when the button is pressed.  We want\n",
       "     * to control all of the cursor setting manually through the\n",
       "     * 'cursor' event from matplotlib */\n",
       "    event.preventDefault();\n",
       "    return false;\n",
       "}\n",
       "\n",
       "mpl.figure.prototype._key_event_extra = function(event, name) {\n",
       "    // Handle any extra behaviour associated with a key event\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.key_event = function(event, name) {\n",
       "\n",
       "    // Prevent repeat events\n",
       "    if (name == 'key_press')\n",
       "    {\n",
       "        if (event.which === this._key)\n",
       "            return;\n",
       "        else\n",
       "            this._key = event.which;\n",
       "    }\n",
       "    if (name == 'key_release')\n",
       "        this._key = null;\n",
       "\n",
       "    var value = '';\n",
       "    if (event.ctrlKey && event.which != 17)\n",
       "        value += \"ctrl+\";\n",
       "    if (event.altKey && event.which != 18)\n",
       "        value += \"alt+\";\n",
       "    if (event.shiftKey && event.which != 16)\n",
       "        value += \"shift+\";\n",
       "\n",
       "    value += 'k';\n",
       "    value += event.which.toString();\n",
       "\n",
       "    this._key_event_extra(event, name);\n",
       "\n",
       "    this.send_message(name, {key: value,\n",
       "                             guiEvent: simpleKeys(event)});\n",
       "    return false;\n",
       "}\n",
       "\n",
       "mpl.figure.prototype.toolbar_button_onclick = function(name) {\n",
       "    if (name == 'download') {\n",
       "        this.handle_save(this, null);\n",
       "    } else {\n",
       "        this.send_message(\"toolbar_button\", {name: name});\n",
       "    }\n",
       "};\n",
       "\n",
       "mpl.figure.prototype.toolbar_button_onmouseover = function(tooltip) {\n",
       "    this.message.textContent = tooltip;\n",
       "};\n",
       "mpl.toolbar_items = [[\"Home\", \"Reset original view\", \"fa fa-home icon-home\", \"home\"], [\"Back\", \"Back to  previous view\", \"fa fa-arrow-left icon-arrow-left\", \"back\"], [\"Forward\", \"Forward to next view\", \"fa fa-arrow-right icon-arrow-right\", \"forward\"], [\"\", \"\", \"\", \"\"], [\"Pan\", \"Pan axes with left mouse, zoom with right\", \"fa fa-arrows icon-move\", \"pan\"], [\"Zoom\", \"Zoom to rectangle\", \"fa fa-square-o icon-check-empty\", \"zoom\"], [\"\", \"\", \"\", \"\"], [\"Download\", \"Download plot\", \"fa fa-floppy-o icon-save\", \"download\"]];\n",
       "\n",
       "mpl.extensions = [\"eps\", \"pdf\", \"png\", \"ps\", \"raw\", \"svg\"];\n",
       "\n",
       "mpl.default_extension = \"png\";var comm_websocket_adapter = function(comm) {\n",
       "    // Create a \"websocket\"-like object which calls the given IPython comm\n",
       "    // object with the appropriate methods. Currently this is a non binary\n",
       "    // socket, so there is still some room for performance tuning.\n",
       "    var ws = {};\n",
       "\n",
       "    ws.close = function() {\n",
       "        comm.close()\n",
       "    };\n",
       "    ws.send = function(m) {\n",
       "        //console.log('sending', m);\n",
       "        comm.send(m);\n",
       "    };\n",
       "    // Register the callback with on_msg.\n",
       "    comm.on_msg(function(msg) {\n",
       "        //console.log('receiving', msg['content']['data'], msg);\n",
       "        // Pass the mpl event to the 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,iVBORw0KGgoAAAANSUhEUgAABQAAAAPACAYAAABq3NR5AAAAAXNSR0IArs4c6QAAQABJREFUeAHs3c2PdOl91/92CJjE2B7Psz0PsZ0MOFFAKIAQbBE7kFiggJBYI7JhA2KBnBlLERAlsGOffwKxRyJLpAQJIYfYM7Znxp7nGTsEwkN+97vy+945U1NVXdXVfXdV9+uSzpx6ONd1rut1+ra6P/6eqosLjQABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIELi6wKeu3lXPBP7wQSNBgAABAgQInJfApx6085qx2RIgQIAAAQIECBC4usCPXb2rngQIECBAgAABAgQIECBAgAABAgQInLqAAPDUr5D5ESBAgAABAgQIECBAgAABAgQIEDhCQAB4BJ6uBAgQIECAAAECBAgQIECAAAECBE5dQAB46lfI/AgQIECAAAECBAgQIECAAAECBAgcISAAPAJPVwIECBAgQIAAAQIECBAgQIAAAQKnLiAAPPUrZH4ECBAgQIAAAQIECBAgQIAAAQIEjhAQAB6BpysBAgQIECBAgAABAgQIECBAgACBUxcQAJ76FTI/AgQIECBAgAABAgQIECBAgAABAkcICACPwNOVAAECBAgQIECAAAECBAgQIECAwKkLCABP/QqZHwECBAgQIECAAAECBAgQIECAAIEjBASAR+DpSoAAAQIECBAgQIAAAQIECBAgQODUBQSAp36FzI8AAQIECBAgQIAAAQIECBAgQIDAEQICwCPwdCVAgAABAgQIECBAgAABAgQIECBw6gICwFO/QuZHgAABAgQIECBAgAABAgQIECBA4AgBAeAReLoSIECAAAECBAgQIECAAAECBAgQOHUBAeCpXyHzI0CAAAECBAgQIECAAAECBAgQIHCEgADwCDxdCRAgQIAAAQIECBAgQIAAAQIECJy6gADw1K+Q+REgQIAAAQIECBAgQIAAAQIECBA4QkAAeASergQIECBAgAABAgQIECBAgAABAgROXUAAeOpXyPwIECBAgAABAgQIECBAgAABAgQIHCEgADwCT1cCBAgQIECAAAECBAgQIECAAAECpy4gADz1K2R+BAgQIECAAAECBAgQIECAAAECBI4QEAAegacrAQIECBAgQIAAAQIECBAgQIAAgVMXEACe+hUyPwIECBAgQIAAAQIECBAgQIAAAQJHCAgAj8DTlQABAgQIECBAgAABAgQIECBAgMCpCwgAT/0KmR8BAgQIECBAgAABAgQIECBAgACBIwQEgEfg6UqAAAECBAgQIECAAAECBAgQIEDg1AUEgKd+hcyPAAECBAgQIECAAAECBAgQIECAwBECAsAj8HQlQIAAAQIECBAgQIAAAQIECBAgcOoCAsBTv0LmR4AAAQIECBAgQIAAAQIECBAgQOAIAQHgEXi6EiBAgAABAgQIECBAgAABAgQIEDh1AQHgqV8h8yNAgAABAgQIECBAgAABAgQIECBwhIAA8Ag8XQkQIECAAAECBAgQIECAAAECBAicuoAA8NSvkPkRIECAAAECBAgQIECAAAECBAgQOEJAAHgEnq4ECBAgQIAAAQIECBAgQIAAAQIETl1AAHjqV8j8CBAgQIAAAQIECBAgQIAAAQIECBwhIAA8Ak9XAgQIECBAgAABAgQIECBAgAABAqcuIAA89StkfgQIECBAgAABAgQIECBAgAABAgSOEBAAHoGnKwECBAgQIECAAAECBAgQIECAAIFTFxAAnvoVMj8CBAgQIECAAAECBAgQIECAAAECRwgIAI/A05UAAQIECBAgQIAAAQIECBAgQIDAqQsIAE/9CpkfAQIECBAgQIAAAQIECBAgQIAAgSMEBIBH4OlKgAABAgQIECBAgAABAgQIECBA4NQFBICnfoXMjwABAgQIECBAgAABAgQIECBAgMARAgLAI/B0JUCAAAECBAgQIECAAAECBAgQIHDqAgLAU79C5keAAAECBAgQIECAAAECBAgQIEDgCAEB4BF4uhIgQIAAAQIECBAgQIAAAQIECBA4dQEB4KlfIfMjQIAAAQIECBAgQIAAAQIECBAgcISAAPAIPF0JECBAgAABAgQIECBAgAABAgQInLqAAPDUr5D5ESBAgAABAgQIECBAgAABAgQIEDhCQAB4BJ6uBAgQIECAAAECBAgQIECAAAECBE5dQAB46lfI/AgQIECAAAECBAgQIECAAAECBAgcISAAPAJPVwIECBAgQIAAAQIECBAgQIAAAQKnLiAAPPUrZH4ECBAgQIAAAQIECBAgQIAAAQIEjhAQAB6BpysBAgQIECBAgAABAgQIECBAgACBUxcQAJ76FTI/AgQIECBAgAABAgQIECBAgAABAkcICACPwNOVAAECBAgQIECAAAECBAgQIECAwKkLCABP/QqZHwECBAgQIECAAAECBAgQIECAAIEjBASAR+DpSoAAAQIECBAgQIAAAQIECBAgQODUBQSAp36FzI8AAQIECBAgQIAAAQIECBAgQIDAEQICwCPwdCVAgAABAgQIECBAgAABAgQIECBw6gI/fuoTND8CBAjcV4FvfOMbey395Zdf3us4BxEgQIAAAQIECBAgQIDA/RQQAN7P627VBAjcIYF9g8JZssBwJOwJECBAgAABAgQIECBwPwQEgPfjOlslAQIEHgocEhgKCx+yeUCAAAECBAgQIECAAIGzFfjU2c78RCb+hw/aiUzFNAgQuEMCh4R0t7lsAeFt6jv3MQKfetCO6a8vAQIECBAgQIAAgXMSUAF4TlfLXAkQuPcCE7idSkC4ax4z13t/0QAQIECAAAECBAgQIEDglgX8v99HXgAVgEcC6k6AwEaB9WBtV5i2fuzGAW/pxV3zvqUpOS2BlYAKQD8IBAgQIECAAAEC90ngx+7TYq2VAAEC5yBwyoHeoX53aS2Hrt3xBAgQIECAAAECBAgQOBUBFYBHXgkVgEcC6k6AwMcE7mpgphLwY5fZkxMQUAF4AhfBFAgQIECAAAECBB6ZgArAR0btRAQIENgtcO7hn5Bv9/X1LgECBAgQIECAAAECBG5LQAB4W/LOS4AAgYXAuYd/s5RCQEHgaNgTIECAAAECBAgQIEDgNAR8C/BpXAezIECAwF4Cm8K1Y8PDTWPOZDaNPcevv9fzeW/62xMgQIAAAQIECBAgQIDA7QsIAG//GpgBAQIEjhK4ydBt19i9tx4CHrUQnQkQIECAAAECBAgQIEDgRgTcAnwjrAYlQIDA/gL7hmi7wrj9z3a9R57inK53hUYjQIAAAQIECBAgQIDA+QsIAM//GloBAQJnLHDO4V/s6/Nff37Gl8bUCRAgQIAAAQIECBAgcGcEBIB35lJaCAEC5yawb1h2TlV25zTXc/t5MV8CBAgQIECAAAECBAhcVUAAeFU5/QgQIECAAAECBAgQIECAAAECBAicgYAA8AwukikSIEDgFAX2rWA8xbmbEwECBAgQIECAAAECBO6TgG8Bvk9X21oJEDgJgX2DM7fTnsTlMgkCBAgQIECAAAECBAicvYAKwLO/hBZAgACB0xfYN/Q8/ZWYIQECBAgQIECAAAECBM5PQAB4ftfMjAkQOGOBfYMw1X9nfJFNnQABAgQIECBAgAABAicmIAA8sQtiOgQI3F2BfcO/cxUQWp7rlTNvAgQIECBAgAABAgTuuoAA8K5fYesjQODsBARpZ3fJTJgAAQIECBAgQIAAAQInLSAAPOnLY3IECNwVgX2r/4R/d+WKWwcBAgQIECBAgAABAgROR0AAeDrXwkwIECBAgAABAgQIECBAgAABAgQIXLvAj1/7iAYkQIAAgYcCKv8eUnhAgAABAgQIECBAgAABArckoALwluCdlgABAndJ4LJbly97/y5ZWAsBAgQIECBAgAABAgROTUAAeGpXxHwIECBAgAABAgQIECBAgAABAgQIXKOAW4CvEdNQBAgQWArsc/uvyrilmMcECBAgQIAAAQIECBAgcBMCKgBvQtWYBAjce4F9wr/7hMTjPl1tayVAgAABAgQIECBA4NQEBICndkXMhwCBeyOg+u/eXGoLJUCAAAECBAgQIECAwK0KuAX4VvmdnACBuyag0u2uXVHrIUCAAAECBAgQIECAwPkLqAA8/2toBQQInKGA6r8zvGimTIAAAQIECBAgQIAAgTMVEACe6YUzbQIEzldA+He+187MCRAgQIAAAQIECBAgcI4CAsBzvGrmTIDASQq4/fckL4tJESBAgAABAgQIECBA4N4L+AzAe/8jAIAAgesQ2Cf8U/l3HdLGIECAAAECBAgQIECAAIFDBVQAHirmeAIECKwJ7BP+rXXxlAABAgQIECBAgAABAgQIPDIBAeAjo3YiAgTus8Af/uEf3uflWzsBAgQIECBAgAABAgQI3KKAW4BvEd+pCRA4f4HLqv8K/l555ZXzX6gVECBAgAABAgQIECBAgMDZCqgAPNtLZ+IECNy2wGXhX/MT/t32VXJ+AgQIECBAgAABAgQIEBAA+hkgQIDAFQT2Cf+uMKwuBAgQIECAAAECBAgQIEDg2gXcAnztpAYkQOCuC+wT/vnMv7v+U2B9BAgQIECAAAECBAgQOB8BFYDnc63MlACBMxGY8M/tv2dywUyTAAECBAgQIECAAAECd1xAAHjHL7DlESBwvQKXVf8J/zZ7v/zyy5vf8CoBAgQIECBAgAABAgQI3LiAAPDGiZ2AAIG7InBZ+DfrVPk3EvYECBAgQIAAAQIECBAgcAoCPgPwFK6CORAgcPICl4V/Kv9O/hKaIAECBAgQIECAAAECBO6tgArAe3vpLZwAgX0FLgv/9h3HcQQIECBAgAABAgQIECBA4DYEBIC3oe6cBAicjcA+4Z/qv7O5nCZKgAABAgQIECBAgACBeyngFuB7edktmgCBfQQuC/8Ef3+smJUv+vhjD48IECBAgAABAgQIECBwSgIqAE/papgLAQJnIzDh39lM2EQJECBAgAABAgQIECBA4N4KqAC8t5fewgkQ2CWwq/pvGf5V9bZ8/qlPfWrXsN4jQIAAAQIECBAgQIAAAQKPXEAF4CMnd0ICBE5dYFf4N3Mv9PvlX/7lj4V/vdfrs82x9gQIECBAgAABAgQIECBA4DYFVADepr5zEyBwUgL7Bn//+3//74uvf/3rF+2r+Juqv8v2LXaOOamFmwwBAgQIECBAgAABAgQI3GkBAeCdvrwWR4DAdQrMrb7//J//84+Ffz/2Y39UTF24t3zc87b6LR/PnHpNI0CAAAECBAgQIECAAAECNy0gALxpYeMTIHAWApdV/034984771z8r//1v1ZB34R6hX6bHvfa8r05JpAez5g91ggQIECAAAECBAgQIECAwE0JCABvSta4BAichcBlwV+LKKj77d/+7YuO/VN/6k9d/OhHP9oY7E3gN6Ff+z/xJ/7Ew2OXj+fY9rMNWM81AgQIECBAgAABAgQIECBwXQICwOuSNA4BAndSYKr0/sW/+Ber9fW87f/9v/+3Cu5m0RPaTfg3od6Efr3+4z/+46swcI6Z9zq212o9bvwZb8a3J0CAAAECBAgQIECAAAECVxXwLcBXldOPAIGzF7is+m/Cv//4H//jqvLv05/+9Gr/J//kn7xoK8Brm7Cu4//v//2/F//n//yf1WcE9iUhf/AHf7C6Zbjbhpdb73VcW33aJlhsP9vZI1sAAQIECBAgQIAAAQIECNy6gArAW78EJkCAwKMUuCz0ay6Fb7VXX3314td+7dcufumXfuniJ37iJx6Gf1OtN4FdYd6EeAV+VQf2vH0B34R5BYCFhW2NUYjYvhCxW4t7PNsEi3P8hIyrifkPAQIECBAgQIAAAQIECBA4QEAAeACWQwkQuF8Cv/Ebv7EK8QrwCuTWW6Fcry/fK/ibyr72hX4TDi7DwMLBtvoW+s055nnn6vUJC+fcgsCRsCdAgAABAgQIECBAgACBfQUEgPtKOY4AgbMXOKT678/9uT+3Wm8h3CFtAsFuF64VAFYh2PY//+f/fBgGFvhNlWCh3jIAnEBxwsAJ/ea4eX7IvBxLgAABAgQIECBAgAABAvdXQAB4f6+9lRO4VwKHhH/BfO5zn7sWn4LACQP/x//4H6vPBJzP/5vbhasE7HFh4wSInXweF/gtqwELC4WA13J5DEKAAAECBAgQIECAAIF7ISAAvBeX2SIJ3F+BfYK/dArVav/sn/2zi8985jOrx9f9n5/8yZ+8aJs2XwpS+Pf7v//7q4rAHjeX3isA7JuD//Sf/tMPH8/nBhYAzjbj2RMgQIAAAQIECBAgQIAAgU0CAsBNKl4jQOBOCBwa/v3jf/yPN4Z/N1VxN9WBfTZgFX6Ff8svFKkysNuEawV/PS/0W94aPJWBd+KCWQQBAgQIECBAgAABAgQI3IiAAPBGWA1KgMBtCxwa/hWsPfPMM5+Y9g9/+MPVZ/dVidc39U413icOPOKFxv7CF76wqgKcqsAqAicALBgsCOy45tl+eXvwnLr3NAIECBAgQIAAAQIECBAgsC4gAFwX8ZwAgbMX2Cf8m1t+C83ee++9i1/5lV/52LoL3L75zW9eFADO7bgFgFXt/Zk/82dWFXkFcW291r4qvfZXac1jeYtwnxfYl4Z07sLAKgM7pv2cq/MVSFYFOJsQ8Cr6+hAgQIAAAQIECBAgQOBuC1ztL9W7bWJ1BAicscAh4V/L/K3f+q2Lf/2v//UqbPvsZz+7Wnnh3wcffHDx0UcfrbYq8Kq4mwCwW3YL4drm9fY9/4mf+IlVEFg41/FXbYWBhXtTCTi3BrevMnC2qQRsvwz/lo+vOgf9CBAgQIAAAQIECBAgQOBuCAgA78Z1tAoCBB4IHBL+TUD2T//pP30Ypg3i3IbbvvCvrXBv+lRtV+A2VXe9N69VuTeVgO2rFqxCsPcPbfXpC0kKJKcasDlN+Ne+eXSeKhpnPs3zpj638NA1OJ4AAQIECBAgQIAAAQIEbl9AAHj718AMCBC4BoFDwr+O/fVf//WLJ5544qKqv/ZPPvnkw1lUfVelXeFdYdpU/E2wN2FgIVtbx87jGaQwrr7L24M7V8+r7DukOvBzn/vcRVuB3zvvvLO6LbhzFgY2t+bVeM27uS3DyQktZ172BAgQIECAAAECBAgQIHD/BASA9++aWzGBOyWwT/DXggvoav/m3/ybi9/4jd9YBXF98cazzz578dRTT63eW/7n85///MVjjz22Om5CtoK2wrW2xptKvMK4qvR63rHzeuN1C29hYH16r6CurRCw8x8S0DVOQeB8PmCViY05cyn8m3V27hl79r12k61r8fLLL9/kKYxNgAABAgQIECBAgAABAlcQEABeAU0XAgRuX2Df4K+ZLkOxf//v//0qrCuAe/zxxzeGf7O6qgK79bZAr9BuqusK1BpzQr+OKYhrm6q83mub1+d5YxQkzmf79ZmBbYV7+7TmPZWIjfGjH/1oNb+pCGxebR2zHPNRhYD7rMExBAgQIECAAAECBAgQIPBoBQSAj9bb2QgQuEWBV1999eKLX/ziKhjrdtlu/d3Vqrbrtt0CtWWYtt6ngLAArq3qvHncfsLBqQ4sCOz1+Sy/ji8AbC67zrE8Z9WEza0wsSrAxmv8xpwAsOOnWrHHAsAUNAIECBAgQIAAAQIECNxPgU/dz2Vf36of/LH9R/cVXt+QRiJAYIfAvpV/y3+aBWFf//rXd4x6c2/93u/93sMv8PjhD3+4CuumKrDQrtCv+fV5g1X3FUy2FQru0xqrcxQktq81XhWAhZdTLVhYWAh43UHgvtfDrcH7XE3HPEqBB/8W/A70KMGdiwABAgQIECBA4FYFVADeKr+TEyBwiMC+YdMy/Gv82wr/Onff4ts23yRclV5bz+c24Cr4upW31+b9uVW4MXa1jqsasCCxvgWCVRlm0POxKBScJvcYCXsCBAgQIECAAAECBAjcDwEB4P24zlZJ4OwF9g3/lgv9W3/rb1385b/8l5cv3drjPkOwLxsp9Gvr1uBa4V+hXfv5zMAeT1Vg4eFlrUCvar/5cpAZqwBwWkGhEHA07AkQIECAAAECBAgQIHC/BASA9+t6Wy2BsxM4JPibardf/uVf/sTn6U1lXbfXFoZd1rpdt5CuYwvv9gniLhuzoG5u8e3YvmSk0K7tgw8+WFUAThVg5/70pz+9uhX4mWeeWfW7bPyOmy8GmZCxMLD5d+65HbggsIpBlYCXiXqfAAECBAgQIECAAAECd0Pg8r+C78Y6rYIAgTsuMOFfy1z/Mo1ur/3oo48e3mL7+c9/fhWKbSLpc/Tef//91S25hWeFZgVx3WZbkNbzPqvvOlrjTsjX5wEW2hVUFgJ27vnikM7b86r81te2Po/lZwc2TlWF7dsK/GZrnMyODQH7bL9DQtr1+XpOgAABAgQIECBAgAABAjcvIAC8eWNnIEDgigKHBkuFWf/kn/yTT5zt+9///kUVfd1qWxBWsPbYY4+twrflwe+8887Fu+++u6rGqyqvgKzwrZCuYG4q6AoAqwxs67Wq+o5pjff444+vvsRjvjSkfeFdcy28bO5tl31zcfOYELBAcdY7AWDv51QVYPvrCAEbUyNAgAABAgQIECBAgACB0xUQAJ7utTEzAvdS4JDQb1n194/+0T+6+OIXv/jQrPdef/31i8K/V199dRWeFXp1m+xU2VVR1629VdcVEP7O7/zOKvwrfJvqvwLAtqoI59bZgr8J0aqkKyCcMLBv8+15oV77fVvzmNuMC+5+8IMfrG4Nbq7NrTl9+OGHq0CveXeeXa0QsFuM69t4ba2pORUkLuffY40AAQIECBAgQIAAAQIE7q6AAPDuXlsrI3CnBSb8m/BqGf4Vmr3xxhsXb7755sXbb7+9CvU6bgK5Hk+/qQqsTwFZVXe1gr2q+6bqr+M7ZxV5HVebMaqu6/gCwj67b6oGO98XvvCF1RirDnv+p/CuCsVCxYK/QrzO23mqBuxxW7cl72pVJnZc8y7AnNuBm3dzaz/rmrXsGs97BAgQIECAAAECBAgQIHCeAgLA87xuZk3gTgocUv0XQKHVP/gH/+DipZde+pjHd77znY+Ff4V7hXIFavMlHO0L2grtCvQK8Hq/YKzHtd6ffgWB0wrSCtU2BYIFio054WFjN85sM8Zl+wLACR8L/to6b5V8c972fZ7hrlalYMfVr/5T/Tf7+rbejhEC7pL0HgECBAgQIECAAAECBM5XQAB4vtfOzAncKYFDwr8Jq/7+3//7nwj/qpar+u+tt95aVf4VfBWkFcBVMddn6LUtq+fmNuBure3YquYKxQoJCwAL9No6b9tUChaotRXOTVVejwvS2ten9wsVuyW4885tvvtcvI4tlGwNfSbhVCnOPNq3jsuCu0LAKgfr3zyn6rF+bZ3jsjH2ma9jCBAgQIAAAQIECBAgQOA0BQSAp3ldzIrAvRLYN/wr8KoVVvXts8tWSNbn/fXZeVUATsVbIdqXv/zl1Wf/felLX1p92cayX48L6grnCsom/CsInErA9ePneUFaYWC31xb4tfW413qv7YMPPlgd3jnee++9VRBYGNm8+oy+y1rh4Ve+8pXV5/8VbnZLc+POZ/vVvxCwisFtLa/OVQhYkNj82k842dzaOm62bWN5nQABAgQIECBAgAABAgTOT0AAeH7XzIwJ3BmBfYO/5YILqH7hF35h+dLq8fvvv78KuPrMvILCArwq+Ar1XnzxxYvnnnvu0m/rLWw7pFU511agV5hWhV2Vfu37LMACwYK2Arve73GvzfuFblUdNsZlrVt9CyRbX+HmhIwFjrVcdt0O3NoKHQs482lOjdfW3Orf4/YaAQIECBAgQIAAAQIECNwtgcv/6rxb67UaAgRORODQ8K/QqnCqW3L/9t/+259YxQRjBVuFcIVrBV6PP/74xVe/+tXV8090usYXOt98rmAB3WyFalMdWGhX2NYce60+hYHzbcSXTafj2qY1XkFjLm27AsD6NL9lGDkhYPsZY5znHNexX6/WvI4xjUGAAAECBAgQIECAAAEC+wso9djfauORD/5Y/qN7Eje+60UCBDYJHBL+zT+xAqq/83f+zup21+Utr4VghW3d/jsVdlX/FZR122u39t5mq+KucLLt3XfffVgZOFV3zbXA8plnnlmtrWrCy1pjffjhh6vbgLutt9Y4jdGtwJeNUWiY1wSAVUlWIThfeDIViZnv27ZdU+HfvoKOe9QCD36+9/8Bf9STcz4CBAgQIECAAAEC1yygAvCaQQ1HgMBugW1B0aZey/DvN3/zNy/+7t/9u6tKtfkii/oUpHVcAVZ/zxd+FWgVhC2/6GPT+I/itSoACyObU/OcSsVCuLYqAFtDtzAXZjbny+ZdZWOt8Qo969fW7cBVFbYV5m1rBY69n+NYTr/m27iHZiOCvm3aXidAgAABAgQIECBAgMDtCwgAb/8amAEBAhsECqFqBVH/6T/9p4tf/dVfXYV8VboVVk1A1eNe6/bWgq1e77bfQz/Pb26DbbybaM3r6aefXlXtFdQVBM7twJ27L/Vo39brX/jCF3ZOYz7PrwrArOpXJeQEefUvzNvWsuo83Ypc346tb/sqAJvvVYLAbefzOgECBAgQIECAAAECBAjcnoDbX460f/AHsluAjzTU/X4I7Fv5N/+kCqD6Rt9f+ZVfWQVUPS/kK9ibW1SPlesbel999dVVFV7nrSquEKzxZ18gNqHYvN48Ch13VdntM7fXX399FfxVxVcgWOscVTH2jcVVMV52jioHW0ehYgHg9K+KsC8+2dUKDwseq0SsX+cqGCxcbP2Zt3aNwF0UePDz7Xegu3hhrYkAAQIECBAgQGCjgArAjSxeJEDgOgX2Df/mnP1dXiD267/+66sQqucFUoVu19GqfHvjjTcu3nrrrYs333xzVQlX0DXhYuearXO3TfjX64VkEwAWmhXU9f6hrc8obJz6VonXbbxV4xUIFuoVSjb2rmrGgr7W01aQV/9uK86v8XaZFWR23vrVv75tPR7zQ9fkeAIECBAgQIAAAQIECBA4PYHD/2I9vTWYEQECJyxwSPhX4FXwVFVa/QrbdoVfV1l2VW+Ffn0JRtVzfZnGhF0FYhP0TeVf77X1fPYT2nVsAWCBXaFglXvLb+m9bH6trRCwceZbg+eW3Cr68mh79tlntw6VUecsuMttPg+wUK8xdt1K3Hmbc/36LMAJAdvPemf9WyfgDQIECBAgQIAAAQIECBA4eQEB4MlfIhMkcD8ECrpq3Rb7L//lv1w9LqC67tYtv91a/N57760CsqmSK+iaOcw5J4DreaHYtEK6jp8qwcK2uX228bqFdt/gsqCtzyysYm9Ctyr4CgRnPr3f+Nta5yoEbE0T5DWPxqzffGnIpv4FgDm3pvlCkNba2ppPrXm0Xo0AAQIECBAgQIAAAQIEzlPAX3RHXrcHfxj7DMAjDXW/mwL7Vv7NP6ECpu9+97sX/+pf/atVNV3VeIVQ19Gq9Kvyr/23vvWtVTA2Id7nP//51WcLFoR1O23VfRN+TSi4DMYKxwro5lbZnvd+rcCscQrk2he89Vl++4aBhX5Tlfjuu++uxi2ce+qpp1aVfAWBE8ptcnn77bdXtw4X/LW+5tC3D1cF2K3E21rnrFqw0LBWn+ZceNg25xQCbhP0+jkKPPh59jvQOV44cyZAgAABAgQIELiSwPWX11xpGjoRIHAfBZbhX7flVvlXADch3CEmBXH9PT9/008lXEHYfNlFt+rWCre6Zbd94VqBV4FjAWCB2wSAHdsc2xqn0K+QrHF63uP2BXe91zl7r4CwqsBeLzxr3ILGCdIad1MrsJv5F8hN0Fig15zadt3S2zpmnnn0uHk0p10BYBadqz7r26x/5rVp3l4jQIAAAQIECBAgQIAAgdMWEACe9vUxOwJnKbBP9d8y/Kva7Vd/9VcfVrwdcuvvBHGFb8sAsPEL5SYEK6zrtUKygrjCv6nQa39IxeGEcxP09TmCPZ7nnWsCuM41YWDVgJe1wrgq/foSkFpjTpDYWLsCwKk6nBCw8/Z4gs9t565fxzbvCU6bf48FgNvUvE6AAAECBAgQIECAAIHzEXD7y5HX6sEfx24BPtJQ97sjcGjw97Wvfe3i7/29v7cCKHxqq8qtcO6yVmhYODahVf0KDtuqIKxNmNVxBVq1wr4J/FYvXNN/CukKBvtm4QK3nk8QN2vqyzyqBCzgmzluO31Vf32mX0FgW+FmLi+++OLFE088sVrnZX37rMMMOtfzzz+/6retCrH/KetW4M5bcNqtv52vwLT5t23ru20eXidwygIP/k35HeiUL5C5ESBAgAABAgQIXKuACsBr5TQYgfspsE/wl8zk5f3d/Tf/5t+8+Ot//a8/BOu9wqpCph5v+9u8YwqpCtsK2Hpem9Cw/oWAM05j9V7j9XqVfjfRqqJr67x9q25bQWDVgYWPBYIFegWWBWyFgB2/rRVS1q+ttU41Y+MW6NV/W6tvFX2tN6M5f/Mp0NvU8slpQtU59/Ka7Loum8b0GgECBAgQIECAAAECBAichoAA8DSug1kQOFuBq4R/v/iLv3jxsz/7sw/XXMg01XIFaD2v4my9FUoV/rUVqBWK1aZCrWCsyrVCvsbZFiKuj3udz6vOK2QrhCv8a46Fce2b9wRr7avm29X67MD6zbf75lKI1zoL6zrHttZ7Wcz5mkNmY7OpX2M2dsd1rvq2L/hr0wgQIECAAAECBAgQIEDgPAUEgOd53cyawK0K7Bv6NcllcPTKK69snPcbb7yxqowrsOsW2Qn01g+eW1QLswoMC/kKrPrcvCeffHJnRd36WDf5vCCtLxdp68s3uoW3MLDblqvgKwjstcK1AsPPfvazW6dTpV+B3HxTb/0LBaskrO+2zwTMZL74o/6dM9+2+m0KR6sYbC4TVHaOtjk279o83zppbxAgQIAAAQIECBAgQIDASQkIAE/qcpgMgdMXuGr49/Wvf/0TiysAe+edd1ZhWEFhYV5Vb5sCpqrRCqOmKq2wqq2gq6q2XbfTfuLEj/CFgskq8ZprQVxbayjELNSrtd5tt+b2fqFcgWfrny8gmVBvWwBYv0w6bqor2y8rATtmvRVeThVm12S8l0Hueh/PCRAgQIAAAQIECBAgQOC0BQSAp319zI7AyQgcEvw16fXAaP2W3m5lffvtt1fhX6FWAVlBXmHY+rEz3oxZQNUxE/7tCs/2Aexz7wrHqnyruq7Qa8418ylMmy/GaH9IK9Qs3Kx6r/EK5TrPPG89u9ZQ/+UcJ0Bsvrtac85oPtevALHHhYBZb2tdizEopJ0tk+Y/+2X/V7ZUd257fdnXYwIECBAgQIAAAQIECBC4WQEB4M36Gp3AnRA4Jvz7pV/6pYunn376Yw7dCts31PbZdhNiFUgVVhWUbQoAG6DwqbCsgKpwqz591l19Dm0FXIWQzaFArGBsqvMm8GrMztdWZVzBX1tVd51312fwrc+n+S5vve1cnXfW9MUvfnHnOiYEbL5zC3D9W8O2MC/Pts5VyNqapwpwfX7L5/l2XEFffdp6reebwr9l33ncumqvPAgG2zQCBAgQIECAAAECBAgQuD0BAeDt2TszgZMVODTwayEFQ8u2KfQprOrz/n77t3979Zl4hW6Fa1/60pdWwVqfWbetEq4AqqBvArmOPST4K9Sr4u7NN99chY7NpQBtgrHeL+ja1iao6/xTDdi+EPD5559ffd7eZfNpnYWhb7311urzAGcuhW0FdX1m4LZbegsgCwn7LMHWUTVf4WlVlN1mvC2M7DMEW9dHH320Cg4LAwtYdwV5jVXI2PVp3/FthZ/Tb/bjNYHfPJ/9ttfnfXsCBAgQIECAAAECBAgQuHkBAeDNGzsDgbMSODT8Ww/+Wmyh03orgCrwaiu0KoCrTSVfX/6xLfyasQrfCpQmDJzXL9t37raCs+9///urirgJtwq4ajPunGOCq9bX1nETEk4oVgDXmL3f+M1/09qX82vuBXaN1WcAFsgV5tW/ALGgb1tFXyFcXs19KvPqn+W2ALB1FC523vp03s5X+NlYm1rz6PiOa231ma3xem18NvX3GgECBAgQIECAAAECBAicloAA8LSuh9kQuDWBQ4O/JloQtGy/+Zu/efFrv/Zrq4Br+XqPq3pr60s/puquirICr0KzQrFuc93WJnQqeCqc2rcVlhX6TQA4t9BOtV9zaLy2ArGpdisIrBV8dWzjFIgV/hWetW/r9VphYFtjbAvWVgc++E/n7HbgboNu7CoAC/HGo8Buzj99Zl84WOhXn+bW+Xvetu28EwA2//rMvLcdPwFgc5jQsH7r13vmZE+AAAECBAgQIECAAAECpy2w/1/Rp70OsyNA4ECBqwR+nWIZAhXGvfLKKxf/7t/9u1V49zf+xt9YBXqFesvWbatV/VXxVtA1IVkVc93W+hf/4l9cHr7xceFaoVdB3T7VZx988MHq8/065/e+972HoV39m19zqGquW3LbF6ztqkAs9Jtqu4LMCfx6rXN0vgK9gsYCzcYt1NzWOv9P/dRPrULR+ldJWNhWSNfWe5vafJ5i16EgNc+CuoLL5557blOXVbDYGhu/dTR+fcdhU6c8Oq511q/QsBCwbfxnX/8JaDeN5TUCBAgQIECAAAECBAgQuF0BAeDt+js7gVsRODb8m+Dnr/21v3bxH/7Df1iFSlWvFXwVonU777JNVV2VZVWjdexnP/vZVUi2LbRa9i9kmwCwcy9DyOVx87gwrkrDqv16XEhWiNU8mlufM9i+isMXXnhha7XdjNe+ebe1vsK0Qs3CvvYFgoWTbcsKw875zDPPLIf52OPGqoKvrUCuwK219nhbK+wrPGwuhXltnbf9rpZ54Wf9C/E6V/0aa1PrWuVVy3t9m5+BTX29RoAAAQIECBAgQIAAAQKnJSAAPK3rYTYEblTg2OCvyRX8FCJ9/etfv/iv//W/PqwQK2B69tlnV7e2ri+isKqKsoKu+nZs4VfHt+1qBVWFbAWAhVCFUpcFgH3LcFuVdYVpBXEFWoVdBXIFlW2FkFdpfQFJW1V7E/j1eCoEp1Kuee4KAFtL49Sv+c6tuQWWvbYtnMsy0xw7bvrtWsu2AHBbn67TVFu2juWatvXxOgECBAgQIECAAAECBAicpoAA8DSvi1kRuBaBqwZ+c/Jl0Natvuvtq1/96qqSrCCrUGpbK0iq0q+Ku8KqwreCr33aN7/5zVWlXUFg/TrP+i3GM05BX5/3961vfWvVZ6ri+nbd+nb77LZba2eMQ/YFiG2FfN3+WyXga6+9tqo4LJirQrC19+2/29ZbNWFViIWVHd862zdet0dvWmthZuvJsuCxPoWGBaXr1ZeznuY51YI5zWcHFqw2h02tkHEqLicAnP28rhJwk5zXCBAgQIAAAQIECBAgcFoCAsDTuh5mQ+BogWNDvyawDP56vin86/UCorZ92oRMBUjbwrD1caqu67MDC9MKmvrijCrZCsA2tW73bavqr4CrVoBWAFhIVxB3Ey2DF198cRX2TZA3t+W2hua86/P2eq9gs4CuUK+ArkCvcG5TANgaqg7svSr1ps+sedMaM2ub24A7R8Fh27a2rACcKsD52Zj9tr5eJ0CAAAECBAgQIECAAIHTERAAns61MBMCRwlcR/DXBJbBzk//9E9f/MN/+A+Pmtd0LgQrdFqOP+9t2lfZ9sYbb6wq4+pXyDZB2abjC8+qgGsrCOs8nbPKtyrp2m6yFVBW5VjlXudeVvM1j+b//PPPb51CAel8zmGhXKFnz6vom8/iW3bufHn03gSAhY5tnW+9dXwBYMcXGhbETgi4fuw8Xx7bmpbbHGNPgAABAgQIECBAgAABAqcvIAA8/WtkhgS2ClxX6NcJCnem/bf/9t8u/u2//berW2bntfV9oVPbvhWA9S+EKny6rBXmvf7666swrWq+KtGqeOvLO7bdajzhX8FhwVZ9CtW6Vfamw7/leqo07PwFbM1lbs3dJwBsjfOZfgV5Pa5/Ieam1pits2M7364AsP7rFYDNs2u4q831WoZ/y5+VHnddNQIECBAgQIAAAQIECBA4XQEB4OleGzMj8DGB6wz7GngZ4vT8F3/xFx8GV31WX2FU4dJ66zP2qm4rpKsCrQq1PltvgqL14+d5n49XZVtVZbu++KMKuj5H780331xV/9W/+XQb765+3/3ud1efwVfo1twn+NvnM/8Kz7rV+NVXX10FbnMrbufOoKCtELEvLmkeu9ZaFWC3Ks8aqmKcLwppjPpvar2eT4Fc8yn4K7Dr3NsCwMLQ3i/4W1YNFpRuanN8wWJVkvXZddtwY3S9Cvj6eWlebT1e//nZdD6vESBAgAABAgQIECBAgMBpCHzyr/vTmJdZELj3Atcd+A3oenBTONU3+vbZebXCrYK9gqX1VvgzlXYTABYOdeyuz/Xr2M5T4FSo1bGbxu98fYZeId5UwlXFV6C17cst6lOI1e2ynaeqto7tHIVq285Tv1rBZKHmO++8c/Gd73xndd4J4nq/ALA5N/8CwCeffPLiL/yFv9BbW1vn7NytofEL9Zpbds2r8Ta1ArrsO2fHF+y139Uaq5Cuc7Q193Fe79e1beyOr3U989rV6jMB4AR/y73qv1163iNAgAABAgQIECBAgMBpCAgAT+M6mAWBlcBNhX4Nvh78FRIV/LUvQCvomWCoAGtTcDZB1tya2vMCqEK6qvQaY701diHYu+++uwqbGrewbtP49a3yrT4FU43X2FXA7QoA51bb1lIgVWBYWFcV4K7W8d3uXNVhAV9VgL3Wucer8XJpzoWEraM59aUf24K8ztm5m9d89mFhXs+rnqxCcFOrcjGXxs0h3/rtahMAdkzznhBw29yuGgA2/jL46/m0sZrn9gQIECBAgAABAgQIECBwWgICwNO6HmZzhwVuMtzbxrYezLzyyisPg60Ji+o7IVePC4gK3DZVdhXIVTVWyNTW8wKqArtCwYKxZauqsOq6grVCtoK0wrmOb5z1wHDCq8YowGou3V7ct/duu/23zwhs7PaNWYD2wgsvXOy69beQ73d+53cuvv3tb1/81m/91irc69wFcMsQrnl0bFV4ExBWJdjjvtDjK1/5ysVLL720NQhsDt2a3PobowBxvhBk3apzzS3VhY1VC2ZbENh+0xd71Kf19l6mraHztG37rMTOMeHg8lpmvan1etepn5eO3xYCburrNQIECBAgQIAAAQIECBA4DYHNf/GdxtzMgsBZC9xG4Ddgy+DvlQeh36a2PGbeL+hpm4BoXp99AV7vt28rHGqb1+e42RdKzVaAVN+Cxcbv8abW+703wVNB1qawbPoWkBU+Nn59J0Sb9zftC/Havve9762q8rIoSCvUq3KwALQgsNbYVe4VFrZvPQV6tdbd+fq25G2tuVf1V7+CvMYr4Nu2psYr0GstzWtC0W0BYJbjX0A3od62+XRsW+PXOkd2eW9qHTfbHF+fNo0AAQIECBAgQIAAAQIEzkNg81985zF3syRw6wK3GfJtWvwylNkW/G3rN4HQtmCufh0zAVKBUaHUcluOXRC1rBqbUK/jC7k2tQma5tjOVTA3YdymPgVqVby19sKwju925G2t25H7huG+ZKTbeetXGNdtuX/+z//5VaVhj+eLNKos7LjmNH0K8HreWgooLwsAm9NU8jXX5rytZdM6smhuGVaFuK3N9ej4MR/35rjeem1eb/wJDdePm+cd29i1jtcIECBAgAABAgQIECBA4PwEBIDnd83M+BYETi3oWydYBjOHBn/LvgU9Ew6tn2OeFzgV/hVSTZhXaNXry9ZYszXm8vhtgd4c337O09jbAsMq6grU2teaV2HbrjV0S3K3DBfqTSXeM888c/HlL3/54q/+1b+6XMLqcUFgW5V4BYWd41vf+taqb2P0XiHatnPWpznVrzGWc85vveXU1ngd39i7AsDGzWp5/K4AcI7NuGvf1jm2tbkm8/70mX3vawQIECBAgAABAgQIECBw2gICwNO+PmZ3jQKnHuIdstSXX375kMMvPbYwZ9q2IGveL0AqoKpP4VYBWLfMrreCoWXoV8hUtVxf5tHrm9oyVOo8U2G3LWTqltwq9Kqoa94Fi7uq/zpnn/332muvrT7Hr+c///M/v/pW311VfB3X+22/+7u/u7qlt88BbKt985vfvPja1762erz+n6oJ55uN2zfX5tzcd30ZSOsvACzM6/bhbS3LCQHz6/j6FRr2+nrLsrHHtOvS8dvaHNvxc32Wx/bajLV83WMCBAgQIECAAAECBAgQOB2BT/51eDpzMxMCBBYC1x36zdAFONMKci4Lc+aYArfCoQmgZozlvvfnuMKojt0W/tVvGTB1nvpuCrHmHIVcE5LNsduqBetTtWDBW/uCryrz+oKRy8K/OV/7jn3sscdWIV7jtPWFHe0bb731WnNqHa1vGdCtHzvPJ3TreX2aa/16fb217ra8xm/6rB87z+f4ed7xu9r6z8Rlx+8ay3sECBAgQIAAAQIECBAg8OgFBICP3twZCewtcJOh34RFe09mceCEgAVSy7Bqccjq4YR4BU4dVwi2K9Ar6FqGS9Nvfdx5PuFf/ebYXQFjt/xWgVdw2NwK57Z9u/CcY9O+LwrpG32rzGusPt+vsTcFgPVvThPeTZh3WdVd66nlUZ+2GWP1xv//n46b6zGvd+zScV6ffcdPu+znYP3YZb95bE+AAAECBAgQIECAAAECpy0gADzt62N21yhQmHaKtwHfVMh3CN0yLFoGPtvGmMCp8GmCqk3HLo8rvJoQcNOxvTZhVPvpuyn0mv6FaG0FXs3jsgrDuV24PgWR3br8Uz/1UzPc3vunn376os8SLPRrrMZt23ZL73JezXXmve2ErXnWn0V92ja1uQYdXxvDbcd3zIzd41p9trVDjt02htcJECBAgAABAgQIECBA4HYFBIC36+/s91DgFAK/62CfYGiCp01jzjG9N493Hb8riNo0/qbjd41fKDYBWePtCi83nW9eK6Cbvo23K6Crz6Y5bZr7jL9pf+jxm8bwGgECBAgQIECAAAECBAjcTwEB4P287lb9CAXuSuB3DNmmAGzTePseN303Hb8rKCu0q8/021UlN+fYtO/z+KZvY00V3qZje23TnGYO2/qsv37o8ev9PSdAgAABAgQIECBAgACB+ysgALy/1/5ernzCuGNvBZ5x7griMlzaFFYt17l+7K7jl4Fbx7VNcLYccx5PkNY5Oq7jC9u2tflMwfp1XJ/H17at9Y3FfVNw/f7gD/5g9YUgfSPwobcBv/XWW6u+3crbNxU3btu2tpxXc515bzu+tYzXmLTf1HIaq97Prm3b8R0zY/e4trymf/TKH//3kGP/uJdHBAgQIECAAAECBAgQIHBKAgLAU7oa5vLIBJYB3noYuHzvkU3oEZ9oQqJdwc+2KU3f3l8Ph9b7dOwyiLrs+I5djl+wdVkAOLfjFrJNCLg+j3n+mc98ZhUA9pl8fYFH39zbZ/kdGgC+//77q76tp7F+8id/8qKxt7WZW++3xuZcCLittY7WXhvDpeOy34R/zWXaOM7z9f3y2KX3+nE9Xz92jqmfRoAAAQIECBAgQIAAAQLnIbD9L9DzmL9ZEjha4D4EfruQ1oOcZeCzqd/y+I6doGrTsRNETZ/Ljp8wb3n8rgCw8K0grX7No4q8gr1trW/p7Ys/2vfNvVUBvvHGGxe/+7u/e/HTP/3T27p97PWO/eCDD1Z9W19jff7zn9/6DcCFjM2pubWuCf+a+7bWmuc61Kfz1G9Ta90TAnbsbPXZ1ub4eb8+u9rMZY657Pg5zp4AAQIECBAgQIAAAQIETkNg+1+IpzE/syBA4BEITKBT0LMe9qyffgKmXr/s+EKo2Tp2wqr1MZfPC7pmPh1fcLatdfttWyFg41dptysAbJwvfOELqxCwfgVtb7755ioAfP3117ed5uHrHVMA+MMf/nDVtzEKFBtzW/v93//91ZyaWxbNdea9rU9rbu21LHZVC47pHD/e7be15TVu/F3HNu7yOs+1mbHXn8/r9gQIECBAgAABAgQIECBwOgIqAE/nWpgJgUcuUHgz25x8GQ7Na8t9Ad0ERgVohVVtm0Kqqtza6tO4HV8QtqtNoDe3zVap1+NNFXOFb5/73OcuPvzww4v33ntvdVtujwutZo7r5/r5n//5hyFc4d93vvOdi3fffffiv//3/37xl/7SX7p49tlnL5544omHn+n3ox/9aPX+f/7P/3kVFnZslYNV/X3xi19cVQ4+88wz66d5+Pztt99eza9KwEiUSgQAAEAASURBVNbQrcLNublvaq21NWfVtcm1KsNtrcBzrFpz1uO+qc9cr2VguOnaTd8JAHs+PyvL/RxnT4AAAQIECBAgQIAAAQKnKyAAPN1rY2YEblSgEGfZer6s9Fq+t3xcyNQ2xxcQTZi0PK7HE0gVStUmBFw92fKfCQxn/MKtQq5NAWBD9KUeBWQdX7hV0PbRRx9dPPbYYxvP8PTTT68q+OYW4Kr5fu/3fm/V97/8l/9y8b3vfW8VzjVurQq+jvn2t7+92jeXvvCjkLAA8Lnnntt4nnmxsZtTc5u5zthzzHI/gV5WrWkCveUxy8cFhW3L4+uT/aa2vF6NP9dz07G9th4AbjvO6wQIECBAgAABAgQIECBwugICwNO9NmZG4EYEJlhr8B7P1vN9AsDCpfrUCocmgFq9sOE/VZdNGDXHbzjs4UsFffXpHB1ftd2Ebg8PWjyYL/ZoXoVsHVsV4LYAsK4vvvji6rjmXkXfhHSFf++8884qUJzAsQByQsWOb25PPfXUxfPPP78ap7F2tcZuTtlW3Vj4t+sLQzq2NXd8bp1v5rLpPFP9l9UEf/XZ1uZ6TWCYc/22tfmZaF+bn5f5GdjWz+sECBAgQIAAAQIECBAgcDoC2/9KPJ05mgkBAjcksAxzJtCZkG5bKDQBXeFUYdJU6G2raqs6r22OL9yqqm7b8VXXddttoVvHFaAV6PX6pz/96U9IFKZ122638lb5V7/vfve7q+O2fbtva/i5n/u5i5deemlV7feDH/xgdQtxt+t2zs6/DLyyePzxx1fHtv+FX/iFVfi3K5hrAq+99toqUGxOrf/JJ59czXVbADjhZWuuFRj2DcPtt7XGzrRrkc94bzu+c3TNahMw7goMC1X7mZjjlz8zqxf9hwABAgQIECBAgAABAgROXkAAePKXyAQJ3IzABH4T6MxZCr6mmmxeW+4Li6bSbMLCQqJtbQLDztPYBVUFVtsCwAKs+RzAAqqpwCuY2xQAdt4+U6+gbELDArz333//YWC3a25f+9rXVtWCVf41fufpnLOmqcDr1uG+7KMQb59vDO4zCZvDVPO1pvnMwm3z6dytofNPOLcr/Gucjs201nWZ+a5e2PCf1jXHd45tQe907Rovw9D1n5eOm5+l6WNPgAABAgQIECBAgAABAqclIAA8rethNgQeqcCEObPv5BPSbatuKzCaALAgqW2qxAqU1lvjtPXeHF8otq1NBVvBV4FYY1e11mf2bbutt/EL1+Zz/dp/8MEHq5CxYLAxt7WCvUK9qv+qzKtvIdwyAKx/n/dXANjtv5e1+jdec2icvBqjLw7Z5tqYnbu11ievDHbNvT4TAHbdJvzbdo7CvMa+jgBw+TPTPDQCBAgQIECAAAECBAgQOF0BAeDpXhszI3CjAhPgtC9sautxQVIB3bbgqcq9KuUKmeZ20gKvwqtu093UCtY6vvCprQq9XiucW2+FWAVtUwnXbb1tvf6lL31p/fCHz7/85S+vjum4wre33nprFSDW96/8lb/y8LhND1p73+S769t8N/Xb9Nrrr79+8cYbb6xu/e3LQwrxCv4ae1d42Jz7PMLmm39VjW1ZbGvLsLJr17XZdntxY8zxU2HYNdlWVTnnLCzsZ2L956Tn01555ZV5aE+AAAECBAgQIECAAAECJyjwyXKdE5ykKREgcL0Cy/Cmx8utM02F2KazFrC1za2jU1VWqLStFYLVp6BtAsYCvm2tYLAwq34dX2hYNeCuysE+m6+tsK25TdDYrbhvvvnmtlNd++vzeYKFnM2htTSn5rar9bl/U/HYmlt7BptC0hknw87R8dnWp21b6xp1fNesa55T12VX69ja8mekx+tt02vrx3hOgAABAgQIECBAgAABArcjIAC8HXdnJXAyApuCnblVd9skqxxbBoAdf1kAuOxTkFf1YMHVptbYVbIVgBVsFVp1fBV1u9oTTzyxCtqmX3OqTwFglXnbzrdrzH3fa+zO0WcJds7O3dyr4iv8a2672nxbcGutX2vIYJzX+3a+TCas7bhdAWDHN6eOL9TrHPXpumxrc2zvb/o52dbP6wQIECBAgAABAgQIECBwWgICwNO6HmZD4JEKTNVW+wKheT5VfdsmU2g0wdEcuysAbJzCqQKnztGxhVdVvG1rc4twfQqvCg0L1uq3rRWa9Xl+hW7dwty56tfttd2W2zcF30SrEq+xO8fcwtu5m0PhX3Nqbttaa+r23ObaWltzlX+7bufNrn75d6767LqdN/O2qQDs+OV13DS3jm0+tc6x/Bnp+by+euA/BAgQIECAAAECBAgQIHCyArvv/TrZaZsYAQLHCkyAU6gz1WAT8BQUFWptC5QKpgqsCpGqEiu46vmuVijXcR3f2PO5fttCqM7d5+YV3nVsc+r22kKpvpBjWzjWF4X87M/+7OrYqvHq22frffjhhxff//73V9/M2+fqFcz1xSHHtALJbjH+9re/vQonC+UK5ArvWm/B38/93M/tPEWVf1UoNr/Wlkd9W/s2/yw6b2vrfIWrc6vxtpM112XI2NhTZbitT9epAHDCv/lZmZ+Tbf28ToAAAQIECBAgQIAAAQKnJSAAPK3rYTYEbkWggGdCnvZT1de+sGdTm4q+uU208Kpt22fKVQ1X6FRIVoBVGFjAVCjVZ+RtagVUBVtVunWe+hZkXRZczZd5zJrq39wK2woB+3y+wrOCthdeeGHrGjfNqdf6dt8Cu8ZoXxjXevLqi1Aat1t++4bhXa01NVZr6nHWGbXmXRWDmWXXOVtj5tvCwjl/x2bQHGsFjV3DbW1+BgoAtwV/nVsjQIAAAQIECBAgQIAAgdMXEACe/jUyQwI3KlCIM0HZPC78KZAqqNsWLBUeFTx1zIRFBW3bAsACrULA+hVgFUYVYhXqbQsAW3iBWsdNaFjfQrfnnntup0shYBWKbfVpa4xCu8YrDGyc5jy3G+/6xt2Za+P0DcP1b2v+vZdda6uysK3wrwrAXe39999fzaExC9rqn9O2b1NurKznnBm2vvplu6stA8DmelkAmFU/A7X5uVjud53LewQIECBAgAABAgQIECBwWgICwNO6HmZzBYFvfOMbFy+//PIVet7vLoU506rwKkgqvCv0KVgqMCqY2hYAFlS1zbEdX6C2LVjq9YKxxq/irYCpSrzmUVi3nM/Mq30VevVt61bZ+nSuquTq1xy2tQK4tiryqtJrft0WXOhXiNZrVQS27rbGmtBwKh8nDK1Pc57gcj4br+MKDgvtCjJfeumlnXOauXb+1157bRVIznq6fbn5Pvvss3PYJ/aFj93SnEPzaW2Xfctw4xdWtq9l2Xx32XXtO375s9HjtgkCPzE5LxAgQIAAAQIECBAgQIDASQoIAE/yspjUIQLCv0O0th874U772lT1FTIVGG1qVZ4VnHVMxxcYVQ3X65tagVPbBI2FgQVrhVnbPo+v+RSMdVzBV/vOU5jXOFXabTvfzOH5559fBYZV+jXPwq0J9Jp7W6FWodx6wFXQ11bwN487tnPnUkBaEFkIVxC4K1Sb+TT/1jKhXOMWaGbQWucazPHL/cx9qvM6/7guj1s+7po0/9Y+c98W7NYvj+XxzWdcluN6TIAAAQIECBAgQIAAAQLnISAAPI/rZJYEbkxgqrnaT9DTyfYJAAuRCsEK0wqM2hc2FWQ13qZW0FVg1/GFWIVhBYC9XvXdptZ7BWMFXc2rgKoqwplvIWCB3K42X/pR2FbfuX23czePWe+2MZY+zaM1tC9ULGBsftuC0uWYrbsqxG7/zar51K9KvsZordta8yykzKz55tUcdvVp/AkAe1yfzrfrluFlANjx4zwh4LZru23eXidAgAABAgQIECBAgACB2xXY/Rfz7c7N2QkQuGGBgpxCofYFPT0uZCroKWAqaCoo2xYwVfVWWFSfKtMKpxqnrS/B2NQKzKqU65yNXf8CscZ46qmntt5yXJ9uB55v9C24K9jqM/26lffP/tk/e2klYMHXz/zMzzycVufv3I3TWltD82jtbbUJvwo7C/3amkvraDukvfHGGw9vRc6q1hg5fuUrX9kZIBauzjciFyJmXN/Mtl2fxs8mo+lTYNjtv7tCvFzy6OehYHW2+ky/2XcOjQABAgQIECBAgAABAgROW0AAeNrXx+wIPBKBCXbaFyxNAFgAVHC0rRWGVUlWmDVVgAVHBWq7WsHVhG4FbT0udCpkK9Da1uYzBOvTOZvbBGl9Nl6hYwHXvq15tDXn5t8cGrMQsLXXxqRxW+8EgfueY45rnoWXhXEz94K7ue33surB+jW/5tacpopvV/iX6zjPOuaazbw27TOY9U+gW//ZNvXxGgECBAgQIECAAAECBAicroAA8HSvjZkReCQChTqFPVPpNvt5fVcA2AQLAAvFCpqq5tsnAKwCbYKpxq9P1XeFU7sCwEK4QsACwAKxxqh/fQvXetz73U57SOu8bds+h/CQsTYd21znS0gK/wrxcp7wrznvavWv4rF1dq2qyMu98HJX61zj3PUsNOxa7br9t/Fy7Dz1mZ+H9j2vzX71xH8IECBAgAABAgQIECBA4OQFBIAnf4lMkMCjE5jAp6qvQqCCtraq4wqONrVCuakCnAq1jt/VZ8KrCQ3rNxV4u/p1/sLDwqk+Q6/WGAVqfSlIc24rGGsNp9C6xfiDDz5Yza8AL8/CtAy6lbiqxda0q81nFmZT365F7vsEgPXJJ4+p/qvvttbxzXFC4fotw79t/bxOgAABAgQIECBAgAABAqcrIAA83WtjZgQemcAEf3PCKsUmBCo8KoAqCNr0RRv1rYKtwKgqtaoAqzwr+Crc2lZtVvhVINW5Cp3qV4j33e9+dzXerqq4KudeeumlVVVd4Vrza+u887xj+my9qgG3hZez3uvet56q9lrPD37wg4eVivkV2k3l34svvnjpqascbIzW1jVpTW1PPvnkzgCwYDSL+nUN55y5bqvgKzzNseML/bo2bc2761+/bX0vXYgDCBAgQIAAAQIECBAgQODWBASAt0bvxAROT2ACnsKeqfwqdCqcK9TaFAC2isKiqsoKiwqP2gqequzbFgDWrzCscQvoOk9bfTpP265beauam3m27zbitgkSezxz+NKXvrQKGzvnTbduRZ5Asoq/1tO6sl2Gf7sCzpljIWJbRo2R83xu4GXVf503gwnzug5do8bY1jpPfrWOm5+B+bnY1m/b64XCr7zyyra3vU6AAAECBAgQIECAAAECj0hAAPiIoJ2GwKkLTGVX+2WwVoBUZVjhUOHTHLe+nt4ryCtAKqzq+AKwgqddYVUhX307T30K7SZ8ah67+hZoPf3006uwsH6drzHaZs5VwtUap9Bw1+2v62s65HnnrwKyLyNpHj0vgMuidWRT6Ff1XvvC0l2tAK8KwsLE1pL7fE7hrmC0MevbHOpX6/ytO+dtrbBu3JY/A12Lns+2q//yPcHfUsNjAgQIECBAgAABAgQI3K7AH32i++3O4azP/uCP5j/6qtCzXsX9mvw3vvGNvRf88ssv733suR84P8oFVlNJN2FazwusCtC6dXdbq/qtz+ab6rPCv+nz2GOPbeu2qnLrVte5bbawqcCs0O6rX/3qzirCGbTAa27/ffvtt1dhVqFi65nbWBuz8Kx5tc3zy8K4OcfsC/rGpoCxoG75PMsJMVt3Bt2CWyXiPq0xv/Wtb60cG3eC0Ob+5S9/eecQGXQNmmPjFBrm+DM/8zM7+02fAtz6ZNM2VYMTBO4cxJsEzkjgwf/O+B3ojK6XqRIgQIAAAQIECBwnoALwOD+977jAMiy862FgfwsXXLUv7JkAq/Cp16dCb9clr8KswKoQqX3BYSFUgeCuALBga8K6gqup3utcBYNVzBXY7Wqdu2MKrOpfINi5J8zstXlc4DefP1h41/O2Aq/WOwadrwCx9de3MVpbY89Y7av2672OrW+3Lzd+8+mz+lpflX/7tObcmjPLsDELENsuG6N51G/W3TVsHrsq/5pT/TpX16C1ZtA2oV9r0ggQIECAAAECBAgQIEDgfAUEgOd77cz8CgLLQO/Q7tP3PgSBy/BvArGCqMKvCQk3+RW+FSS1FeRNaFhg1uPG2tYKtzqmCrYCrAK1ximgq1V52Pi7WvPuuPp2/sKwwq35YovW0FyW6+uYArvmVgjYe20TerXe+jVm85s1LZ/3fq1+zXEq56q8e/bZZ1fPd8173muuVeJNRWHjNrfG2SdEbM2tb8K85lIIeVkAmHfXtvO17ixm6/lYzDztCRAgQIAAAQIECBAgQOC8BJR1HHm9HoQDbgE+0vBRdp8Q77rPeZdCwQm8JvSqwq1bcyf8K2Crmq+wa1srVHvttddW4dsEf88888wqnCvM2tX6xttCwMKsgqnmU6BWiPXCCy/s/GKQTeMWhvW5fIVrjde46+Hd+j/j9cBr/f0JCCcka27Lir/W2JwPaX37cc4T4DWHwrsq/77yla+sxt813lXNu86vv/76qqKxcxaCFja2nsLHttq6ya65eI/AOQg8+Jn2O9A5XChzJECAAAECBAgQuBYBFYDXwmiQ+y6wDBbPPQzsb+LZCroKuQqBltVvhWm7grzeq/qsoK2Aqb5TaddYu8KxwsUJ3JpHgV0hXq91a2xjXhZALn8eO1efvVew1rwLFws12xqrsZvjbHPu2U9GMKFf+0KyCccKygrM5vbjPuvvkNZ5m1Nb82tevVbwl2OVkZ1jV2sNVTm2z7o5N8f677pOjdk5c6i1ttbVNZ/19voY9FgjQIAAAQIECBAgQIAAgfMTEACe3zUz4xMXKAy8CyFgAVAh2ASAE+YVUBUaFawVFm1rVcXNbaX1qfpuxtoVAPZeVYaFTm31K9Qq3OrW2ObR1rf/7ts6b8FcW/NqzObUVjA44eZUObbuZQDYPBpjtsK+AraCubYCyV0W2+Y5a6ryr4C089caaz7zr7F3teZe+DeVjR3b3KZ6cFffHLuWBY6zxmUAuAwBd43jPQIECBAgQIAAAQIECBA4bQG3vxx5fR6EBG4BPtLwNrovK/Zu+vznGAZOANa+gKlwakKznhfSFaZV+batFaz1eXaFWxMy1W9uad3Wb14vnKsq7o033nhYpVY41RiFeFX17Tr/jLNrPwFgwVvr6nlh2Ky/vhNEdu4Jxzr/BIC7xr/svW69LbzLKuNa56h6svUV/l1W/ZdxTgWAOU//bJ577rmdU+jadP6ubf06V8FhVYOFna29EFAjcBcFHvx8+x3oLl5YayJAgAABAgQIENgo4JffjSz7v/ggKBAA7s91FkfedDh4LoHg/GgXiFWpVjhWUDRBVWFRAdOugKi+86UW3b7b84Kmqve6tXWf22WnfwFXFYCN0d/thVQFce0LFS+rlDuFH75ltV6Ora0qvKznFuKp/KsK8rJWePe9733vYUA6VX/Z1n/XtZlvGu66TphZ/65rIev0lZFcdhW8f64CAsBzvXLmTYAAAQIECBAgcBUBAeBV1BZ9BIALjDv48L6HgQVTbVXItRUEVjVWkFVA9NRTT10a4hUwVaH25ptvrkLExqs6rfCv0O6yz6jrx6pztvVlHs2h0KzzT8VaFXN9yUjh2akGVq37+9///ipArVKvMLWtVihaYNf8s9mnsrEAMddMujatO8/Cv30qB6us7PwFqgWphX5thYFtjXeqlis0/yFwpMCDn2+/Ax1pqDsBAgQIECBAgMD5CPjl98hr9eCPehWARxqeU/ebDARPsTJwfrwL3NoK8wriCpwKjiasqnJsVyvw+s53vrOqVKtvgVehUyFggdU+IWBjdMtst8sWXBUE9vd7t6oWBBagTUXgPpWFu+Z7ne/lVFjX9u67767mnmMhavMvbMuvADPPwszLWmNVDVkFYNej1jiNUfh32Rj5FUZOkFrwWPiXY9emrSYfWTH4zx0VePDz7XegO3ptLYsAAQIECBAgQOCTAn75/aTJQa8IAA/iunMH35dAsBCrrcBovpW3EKuAcG47vaxqraCq0Kr+BVcTfhVWPf/883uFgP0AFX7V/+23316FaQWKU5E4YWBzaV6N3VYw+KhaJoWVs1Wll1Vbc601zwK7qforsNwnBK1vlXutv/Ebr7FaYyHeCy+80CE7W9eha1CQuvRqPgV/VVbO7b87B/ImgTMXEACe+QU0fQIECBAgQIAAgYMEBIAHcX3yYAHgJ03u4yt3PQgs1JoQsOCo6rvZCo0KoAqxLgvaCgDb3nnnnVWYOCFgtxIX2u1buVcQObfTTjVg86k1ZuFf1Wxzi/AEglW5FXrdRCuEbA6tr33zKvTLq/lm2Lnzah7Nsc9CbM1TcXfZvAr+5ktRCv8as+CvrUrKtl2tuUx42HzzaS4zhwn/MtQI3HUBAeBdv8LWR4AAAQIECBAgsBTwV95S4wqPBYBXQLvjXW4qDLzNW4QLmmabYKuQq0qyWtVjhX8FeZe1KuLm8+cKxgqbCqAKsdpfFmLN+M2j81cRWDVclW2NN7fWFmYVrLVVXTeBYOfp8VS8zXhX2Xe+grjmMVV+Pe613pv5NHbhX0ZTmdicnnzyyb1PW7A4lY+N2/VobY1RwNmtv7ta1yujQsDm1/oL/9qaS/PrWsy2ayzvEbgLAgLAu3AVrYEAAQIECBAgQGBfAQHgvlJbjhMAboHx8kOBuxQIFjpNsFWIVKBU8FV1YIHbE088sQqjCpF2tULA+hbeFeQ1VmFWYVSfY1cQ2ONDWrcEFwQ2XkHZhG/NrXk3p+ZY8FXYNVthYI8nMJwAbNZQ37aCxamCbM09n20M1s/VuIV+raXgL599b/WdtRfcFZhO5WXnzKrxGqvbfmeu02d93/wbo7HmWo3xBIGz7vW+nhO4qwIPfuZ3/w/VXV24dREgQIAAAQIECNxLAb/8HnnZH/xh7UtAjjS8L93vQhDYj3sBUiFUAVvVbhPg9V6B1FTyXXZdq2ibELDArv4FZlUAzhiHhmXzuXvvv//+KuwqWGybgK5zzN/8BX49nvBvHs/7s69PbcK/nrf+Cftm3zH1mYCxfVvrmc8ibH9IW37e4YSOBZYT/lX5t0/FZL7ddt28Z471a36Fie1rs+ZD5uhYAucq8ODn3e9A53rxzJsAAQIECBAgQOBgAb/8Hkz28Q4PwgAB4MdJPLtE4NyDwAm8CpMKqKoqaysQLEwrnCrAa7usFQBOOFVQ19hTLTdh4mW3tm46R+MWTM7cmudUBLbvn+1s9Z8cYH0/Y88/8/V9x7dNiNi+kK+QrnW09UUfE7DNePvsq5KcW5wzyqbQcnm7dAHgZa0xZmuuU2nZOM1r5j5rv2w87xO4KwIPfub9DnRXLqZ1ECBAgAABAgQIXCrgl99LiXYf8CAQEADuJvLuFoFzDQL7kW8rkCpkqzJtgrbCtQKmQrBu5d2nFU69+eabq9t3G6uxJ6Tq9tnGObRybs7bvOaW4MYuZOz5VPC1hh7Pmtpva5MVtC+IW24FfVNJV/A5IWb7q7QqGX/wgx88vO23eXe+gsU+Z3FCwMvGXn7px1yb5tRWwDrhX2NrBO6bgADwvl1x6yVAgAABAgQI3G8BAeCR1/9BYLA9MThybN3vl8A5BYLzYz9VdYVrbVWqFajV5lbefcK76VsY2O27BV6do2CqoGpuee1bc49tjd3W3NsKCScQnH3nmDVO8LcMywriCvxmK6g8tlWlOFshZY8LKDt/oWLnaN/nCO7TChH7XMRsa12Hgr/mPtusbfb7jOsYAndF4MHPvd+B7srFtA4CBAgQIECAAIFLBfzyeynR7gMehAQCwN1E3j1A4KZCwKZw3d8i3I9+AVVbQVqhWrfeFjgVpBU2FToVWu0TAjZOAWKVbwVgjdc4BW8FVo03X6JRIHhdbdbQuVpH65qtc5QRtE0A2L7g7zpb650vMSmQHMPO1W26rTvHKv/2qdYr/CtM7XMWW0vj1DfDxpvgctZ3nWsxFoFzERAAnsuVMk8CBAgQIECAAIHrEBAAHqn44I9rAeCRhrpvFripMPA6g8B+/AvQJjwrdCq8auu9wqZDQsAk5tbXQqzCsILAArhCq0Ksxqsars/WO/fW+rpNt63P/MttQsjsJkTtNuh9qwwn/Bu/wr8C1Owm/Ou1gsRcNQL3VUAAeF+vvHUTIECAAAECBO6ngL/+jrzuAsAjAXXfS+CUw8CCvtmmeq0gqzCwYLCQqeBpqtj2WvCDg7oFtjHee++9VTBWEFjYWHhVGFhFXJWAbRMK7jv2bR2XU8FcW1V/Vei1rkK/rArlsiqwe/zxx1drPOSbkJfVk50rm0LExitAzK7rMduxDpf9XF5n2HzsXPUnsC4gAFwX8ZwAAQIECBAgQOAuC/jk97t8da3tzgjcVJByWYCzD+CESe0L5gqbZitwqhVyFeZ1i2/B1D6t4KsQrM8S7HEhXwFZQVkBY2MVoPWZge++++7D24/3Gfs2jmnO3ebbXJtzcy8I7PXWVPDXGltra27t+4Z/mebReFnXsp/wr+tx3eHfbRg6JwECBAgQIECAAAECBAhcTeB6P8jqanPQiwCBPQSWIeB1BHdzyhlrOf68t+++8K82gV+PC51qBVIFXFW8FVRVxVfAtU8rFOtW34KsbpOt/4zXvsq5ArTG7hyds2MbvzDyFFqVjM2v6shCupl3c8+iUHNZ0VilXltr36c1TufIp8e1rkP9s2jf+Mugdp9xHUOAAAECBAgQIECAAAECd0fALcBHXssHgcR+5UxHnkd3AtsEJsDb9v6hr/cj/corrxzabXV8fWcr8JqArpCu6rTeK/Dq8+iqbtv3c+1mMlURNk7bhGkFfwVptcZu69bXQq/Cr7YJ1CYUm/Gua9/5C/baCuLmcfvCvwn7mmtBXHMspMugkK7Kv26R7vVDWueqorDzNIfWPNWX7XOYyr9Dx95nHpf97B0TKu9zfscQOEbgwb9FvwMdA6gvAQIECBAgQIDAWQmcRonMWZGZLIHTEihkuSyIOWTG/U1cANh2aFv+PT1VZ43R64VUE4QVXBVY9bwQbFk5uOucfe5fYV4BYIHWfOZgYxcyNuYEghOyNXbHN5+2CQd7veeFghOS7Tr3vFeI2fk6T/NvPyFfzyf47PU5pj4ZzBw6ZwFdlYqFf3N785zjsn3jVvXXulr3jF+Y2LiNv1zX8rpcNrb3CRAgQIAAAQIECBAgQODuCfh/v4+8pg/+8FYBeKSh7tcrcF1h4PxoXyUIrG9h3ARyBWYFVoVVEwIWhk2VWp93d2grRCx4a2vcxp9AroBs2gRvEwhOMFbo12s9bz/vd/x6YNZ6GnPWNeeZ9S3Dx1nf+DXWjF3Q15pn3d3evH6umfe2feudtY9n4+dZQLoM/ybYPPQc2869/vqunzXVf+tanp+awIN/F34HOrWLYj4ECBAgQIAAAQI3JqAC8MZoDUzgdgSuqyJw/jYuAGw7pNW3UGrZqk7rtUKrqZibAK3QqkrA9T7L/uuPqwRsm/EKwHpc2Nj4y8Cu12rNqwq9zjNznCBwnrefY9tPkNdce9zW2D2f1+Zcc2z95hytu7m1xm597nlViIWAh7bO99FHH63W0Jp63vwbe4LFzrVc06zn0HM5ngABAgQIECBAgAABAgTujoD/9/vIa/ngD34VgEca6n6zAruqtA458/yoHxIGTp/2VcbNVgVbAdZU0hWWFWIVjlUNWIh11VYo1ucDThjYOSYQLKjr+QR3zWvmeMj5JlRr39Z8W0PBW0Fca5nQr6q8CT8POcfy2EK/1lB42VbrXJ1jPt+wx517gseOmXn2+Cbarp8tFYA3IW7M6xR48O/D70DXCWosAgQIECBAgACBkxa4+l/ZJ70skyNAYAQmiNkV1syxu/bzt/IrB1QE1qeArX2BVa3HE4gVVhXIFQYW2BXQ9dpVPhdv5l7/+ZbhCQELHDtP25ynELBtWb23KxBs3rO1lh53rraCtwnkmnsB4FTlzbyusm9+61940vnmXBOa9ryt95pXbfZXOa8+BAgQIECAAAECBAgQIHC3BPy/30dezweBgQrAIw11f3QCx4aAM9N+7AsCD2lTdde+EG6Cv8K/qtvm/QnO5lbZKtyuo835ll/S0Tk7f+vp8cyh880/7QnSJlybAG4Ct/bNdSr+2l9Hm8/5a76FgM2tuRQ2zrmy6nbi5jTbzPc65nDZGLt+niZ4vmwM7xO4LYEH/1b8DnRb+M5LgAABAgQIECDwyAX88nsk+YOQQAB4pKHutyOwK7zZd0bz479vGDjHT9hW9V1bAeDcElw4OMcValVRV8DWdl1h4L7re5THteaCvrnNd0K/5lDgV9DYPo953mttk2PM/lHNe9fPkADwUV0F57mqgADwqnL6ESBAgAABAgQInKPAxz+l/xxXYM4ECFxJ4DoCmgmc9g0AO362wr2pZJtqv6rZqmqbUKugsHDwRz/60erLL+abfq+04BPuVHXi+++/v1pja52KyKyyqOovm/nykIyyy3A822sECBAgQIAAAQIECBAgQGCTgABwk4rXCNwTgULAY4PACZ72DQGjrc/cslrANbf9TtBVIDhBYBWCc5twX4bxwQcfrIKyu3CJqvabNRVuTiVkax6XCUezmdt/e2/8xv8ueFgDAQIECBAgQIAAAQIECNyMgJKRI10f3LbnFuAjDXU/HYFdt3TuM8vlP4d9A8H6LLfCvrk1uH0h2fJ58yj0qgJutsLCz3zmM6vQcJ953tYxfa7ffDHJfCHJfO7ghKKFe20Ff/O44G+q/drXTiH42/XzcmywfFvXyHnvj8CDf0N+B7o/l9tKCRAgQIAAAQL3XuB6Pq3+3jMCIHA3BCa02RXs7Fppf08vQ8Bdx8578zf49J1bWwu/Cv5qEwAWnk1gVoDW673WsQWH9S0M7LMCJyib89zWvgCzuU0V46wlp9bSupvrVDwuQ78e9177jlta3dZ6nJcAAQIECBAgQIAAAQIEzk/ALcDnd83MmMCNC0wQeJUTTVC1bwVg51gGWxOATZhX9VvVcLP1vLCsPhOiFa71pRndRvvDH/5w9Xl6H3744ep5AdyjbAV8nbO5NIe25tTzCQMnxGytrWXbGjMQ/j3Kq+dcBAgQIECAAAECBAgQuJsCKgDv5nW1KgJHC0wIeNVqwCYwIeDsd01qQsCOqeqtcG9eKzCbIGwq6Hqv1ydMm9erDGxrjKmeW1bXze20jXcdrfMWQM75p0qx582tubSWtlrzbmtuzWvmOPNaznuOnf11zNcYBAgQIECAAAECBAgQIHD/BHz+zZHX/MEf9T4D8EhD3c9D4JggsH8m+4SAS4n5p9V+Hheq9XiCtfbLoK3nsy37zbgTpC33BW49r/V42rw2516O1znmefue1+bY5Rgzfvu5vbnHBZDL573WVpv5zePViyf4n10/ExMgn+C0TYnASuDBvzO/A/lZIECAAAECBAgQuDcCfvk98lI/+INfAHikoe7nJbAr9Nm1kvmncmwQ2DhtE/z1uACuSrsJ/5bvzfvta7OfuS4zgOXjeX+5X/ZdPu6Y6VuI1+PZCvom3FsGfr0/7/V4Gf4tx1ue/9Qe7/pZEACe2tUyn3WBB//u/A60juI5AQIECBAgQIDAnRX443KXO7tECyNA4JQEDg0A52/09m3LMK1ArW0+R28+S2+573Hvt3VsodsEbxO6jU+h3q5tjlvOo7GW81ify/pnF848Zi6znllf5+ixRoAAAQIECBAgQIAAAQIErkvAX5lHSj4IC1QAHmmo+/kK7KoA27Wq+WdzaBg4Y07/ns/j5W24Pe712ZbP5/H0nWPm+XLf49oEcsv9MrArxJvny8e9tv58OcYfjf7H48/zc9nvuv4qAM/lKt7feT74t+h3oPt7+a2cAAECBAgQIHDvBFQA3rtLbsEErk/g2JDnqgFgf7fP3+7zeCrp2i+r8qZCcKoAqwi8bJuqvdmvVxRu6j/jLyv81qv8lmFgV2Hmfn1XxEgECBAgQIAAAQIECBAgQOCTAgLAT5p4hQCBAwSuEgIug6+rhoBNccZZ7pdB4ISBmwLB9WCw4K5gbxnkXfZ4xmg/25xz5tHcetx+fc6rF/yHAAECBAgQIECAAAECBAjcsIAA8IaBDU/gPggUAl41CMznmBBwfCcE7Pk8nuBtwrj28/l/7Web8O6Q/fRdjj2Pl6HfevA387UnQIAAAQIECBAgQIAAAQKPSsDn3xwp/eDzw3wG4JGGut89gV2fDbdttfNP6TrCwE3nmPE3vbd8bdtxE+Qtj930eN/jNvU9t9d2XeerBMLntn7zPW+BB/9W/Q503pfQ7AkQIECAAAECBA4QUAF4AJZDCRDYT+Aq4c9N/y3e+Ju29RVtOmbT3PY9bn18zwkQIECAAAECBAgQIECAwKMWEAA+anHnI0Bgq0Ch2k1VAG476bYg77LXt43ndQIECBAgQIAAAQIECBAgcGoCP35qEzIfAgTuhsCyCnDXraLrq12GgI86DFyfi+cECBAgQIAAAQIECBAgQOAuCKgAvAtX0RoInLjAMgzcZ6pTfScA3EfLMQQIECBAgAABAgQIECBAYLeAAHC3j3cJELgmgUNDwE5bEKgRIECAAAECBAgQIECAAAECxwkIAI/z05sAgQMECgEPDQK7fVgl4AHIDiVAgAABAgQIECBAgAABAmsCAsA1EE8JEDg9geXnAp7e7MyIAAECBAgQIECAAAECBAictoAA8LSvj9kRuJMCV6kEdDvwnfxRsCgCBAgQIECAAAECBAgQeAQCAsBHgOwUBAhsFjj0dmC3Am929CoBAgQIECBAgAABAgQIENglIADcpeM9AgRuXOCQELAqQJ8JeOOXxAkIECBAgAABAgQIECBA4I4JCADv2AW1HALnKHBICNj63A58jlfZnAkQIECAAAECBAgQIEDgtgQEgLcl77wECHxM4NAQ8GOdPSFAgAABAgQIECBAgAABAgS2CggAt9J4gwCBRy1wSAjoVuBHfXWcjwABAgQIECBAgAABAgTOVUAAeK5XzrwJ3FGBQ0LA+UzAO0phWQQIECBAgAABAgQIECBA4FoEBIDXwmgQAgSuU+CQELDzVg2oESBAgAABAgQIECBAgAABApsFBICbXbxKgMAtCxQCHhoE3vKUnZ4AAQIECBAgQIAAAQIECJykgADwJC+LSREgcKhAVYAqAQ9Vu/rxrK9upycBAgQIECBAgAABAgQetYAA8FGLOx8BAgcJHFoFKJg6iNfBBAgQIECAAAECBAgQIHAPBASA9+AiWyKBcxc4NAQ89/WaPwECBAgQIECAAAECBAgQuE4BAeB1ahqLAIEbEzgkBFQFeGOXwcAECBAgQIAAAQIECBAgcIYCAsAzvGimTOC+CggBT+fKH3ItTmfWZkKAAAECBAgQIECAAIH7KSAAvJ/X3aoJECBAgAABAgQIECBAgAABAgTuiYAA8J5caMskcFcEqjzbt/qsW4HdDnxXrrx1ECBAgAABAgQIECBAgMBVBQSAV5XTjwABAgQIECBAgAABAgQIECBAgMAZCAgAz+AimSIBAgQIECBAgAABAgQIECBA4P9r74521dhhKID2SP3/X27hYSQYDRwyE+PEWZWqXg6dOF7mad+EEiBwVkAAeFbOcwQIpAp8eg34vklXgVNHpTgBAgQIECBAgAABAgQIJAsIAJMHoDwBAucFWkLAexXfB3jeev8ky72I1wQIECBAgAABAgQIEBhXQAA47mzsjACBDwRaQ8APlvRXCBAgQIAAAQIECBAgQIBAKQEBYKlxaobAmgItIaCTa2t+RnRNgAABAgQIECBAgACBlQUEgCtPX+8ECgm0hoCCwELD1woBAgQIECBAgAABAgQIvBUQAL7l8SYBAjMJtISA976EgDHTbZ1DzC6sSoAAAQIECBAgQIAAAQKbgABwk/AnAQIECBAgQIAAAQIECBAgQIAAgYICAsCCQ9USgZUFWk+fOQW48qdF7wQIECBAgAABAgQIEFhDQAC4xpx1SWApASHgUuPWLAECBAgQIECAAAECBAj8IiAA/AXI2wQIzCkgBJxzbnZNgAABAgQIECBAgAABAv0FBID9Ta1IgMAgAvcQsDUIHGTrtkGAAAECBAgQIECAAAECBLoJCAC7UVqIAIHZBXwf4OwTtH8CBAgQIECAAAECBAgQOBIQAB6p+BkBAssKCAGXHb3GCRAgQIAAAQIECBAgUFZAAFh2tBojQGATaL0GLATc5PxJgAABAgQIECBAgAABAhUEBIAVpqgHAgR+FWj9PkAh4GtSNq9tvEOAAAECBAgQIECAAIERBQSAI07FnggQIECAAAECBAgQIECAAAECBAh0EhAAdoK0DAEC9QScdGufaet16/YKniBAgAABAgQIECBAgACBVgEBYKuYv0+AwNQCZ64CCwKnHrnNEyBAgAABAgQIECBAYHkBAeDyHwEABNYUcFJtzbnrmgABAgQIECBAgAABAisKCABXnLqeCRAgQIAAAQIECBAgQIAAAQIElhEQAC4zao0SILAXaDkF6BrwXs9rAgQIECBAgAABAgQIEJhFQAA4y6TskwCBEAEhYAirRQkQIECAAAECBAgQIEBgIAEB4EDDsBUCBMYXcBJw/BnZIQECBAgQIECAAAECBAg8CwgAnz28IkBgQYGWU4AL8miZAAECBAgQIECAAAECBCYXEABOPkDbJ0Cgj0BLCOgUYB9zqxAgQIAAAQIECBAgQIDAdwQEgN9xVoUAAQIECBAgQIAAAQIECBAgQIBAioAAMIVdUQIERhRwCnDEqdgTAQIECBAgQIAAAQIECFwVEABeFfQ8AQKlBISApcapGQIECBAgQIAAAQIECBC4CQgAfQwIECCwExAC7kAeXvr+wwcM/0mAAAECBAgQIECAAIFJBASAkwzKNgkQIDC6QEtwOnov9keAAAECBAgQIECAAIFKAgLAStPUCwEC3QRawqz7qTgn47rRW4gAAQIECBAgQIAAAQIEOgv87bye5QgQILCswFEI+Bgk7t9/fG9ZNI0TIECAAAECBAgQIECAQLjAT3iF4gX+3X4Vb1F7BJYW2Id2I2N8I1B85/GN+iP729tcAj+3X3Pt2G4JECBAgAABAgQInBdwBfi8nScJEFhA4FWo9ernmST3cG77HbGPd+FfRD1rEiBAgAABAgQIECBAgEAfAf/3+6KjE4AXAT1OYBKBo/BrCwGP3pukra7b3Dy6LmoxAkECTgAGwVqWAAECBAgQIEBgSAHfATjkWGyKAIEZBAR/M0zJHgkQIECAAAECBAgQIEDAFWCfAQIECBAgQIAAAQIECBAgQIAAAQKFBVwBvjhcV4AvAnqcwGQCTv0dD8z132MXPx1XwBXgcWdjZwQIECBAgAABAv0FnADsb2pFAgQKCwi6Cg9XawQIECBAgAABAgQIECgq4DsAiw5WWwQIEPiGgED0G8pqECBAgAABAgQIECBA4JqAK8DX/P64AnwR0OMEigtUvjIs/Cv+4S3enivAxQesPQIECBAgQIAAgScBJwCfOLwgQIBAX4HMkKxy+Nh3SlYjQIAAAQIECBAgQIBAbQEnAC/O1wnAi4AeJ0CAAAECCQJOACagK0mAAAECBAgQIJAm4B8BSaNXmAABAgQIECBAgAABAgQIECBAgEC8gAAw3lgFAgQIECBAgAABAgQIECBAgAABAmkCAsA0eoUJECBAgAABAgQIECBAgAABAgQIxAsIAOONVSBAgAABAgQIECBAgAABAgQIECCQJiAATKNXmAABAgQIECBAgAABAgQIECBAgEC8gAAw3lgFAgQIECBAgAABAgQIECBAgAABAmkCAsA0eoUJECBAgAABAgQIECBAgAABAgQIxAsIAOONVSBAgAABAgQIECBAgAABAgQIECCQJiAATKNXmAABAgQIECBAgAABAgQIECBAgEC8gAAw3lgFAgQIECBAgAABAgQIECBAgAABAmkCAsA0eoUJECBAgAABAgQIECBAgAABAgQIxAsIAOONVSBAgAABAgQIECBAgAABAgQIECCQJiAATKNXmAABAgQIECBAgAABAgQIECBAgEC8gAAw3lgFAgQIECBAgAABAgQIECBAgAABAmkCAsA0eoUJECBAgAABAgQIECBAgAABAgQIxAsIAOONVSBAgAABAgQIECBAgAABAgQIECCQJiAATKNXmAABAgQIECBAgAABAgQIECBAgEC8gAAw3lgFAgQIECBAgAABAgQIECBAgAABAmkCAsA0eoUJECBAgAABAgQIECBAgAABAgQIxAsIAOONVSBAgAABAgQIECBAgAABAgQIECCQJiAATKNXmAABAgQIECBAgAABAgQIECBAgEC8gAAw3lgFAgQIECBAgAABAgQIECBAgAABAmkCAsA0eoUJECBAgAABAgQIECBAgAABAgQIxAsIAOONVSBAgAABAgQIECBAgAABAgQIECCQJiAATKNXmAABAgQIECBAgAABAgQIECBAgEC8gAAw3lgFAgQIECBAgAABAgQIECBAgAABAmkCAsA0eoUJECBAgAABAgQIECBAgAABAgQIxAsIAOONVSBAgAABAgQIECBAgAABAgQIECCQJiAATKNXmAABAgQIECBAgAABAgQIECBAgEC8gAAw3lgFAgQIECBAgAABAgQIECBAgAABAmkCAsA0eoUJECBAgAABAgQIECBAgAABAgQIxAsIAOONVSBAgAABAgQIECBAgAABAgQIECCQJiAATKNXmAABAgQIECBAgAABAgQIECBAgEC8gAAw3lgFAgQIECBAgAABAgQIECBAgAABAmkCAsA0eoUJECBAgAABAgQIECBAgAABAgQIxAsIAOONVSBAgAABAgQIECBAgAABAgQIECCQJiAATKNXmAABAgQIECBAgAABAgQIECBAgEC8gAAw3lgFAgQIECBAgAABAgQIECBAgAABAmkCAsA0eoUJECBAgAABAgQIECBAgAABAgQIxAsIAOONVSBAgAABAgQIECBAgAABAgQIECCQJiAATKNXmAABAgQIECBAgAABAgQIECBAgEC8gAAw3lgFAgQIECBAgAABAgQIECBAgAABAmkCAsA0eoUJECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECA8X5yuAAAmQElEQVRAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAgAABAgQIECBAoJzAf/M3Z/W2oVjIAAAAAElFTkSuQmCC\" width=\"640\">"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "ename": "KeyboardInterrupt",
     "evalue": "",
     "output_type": "error",
     "traceback": [
      "\u001b[0;31m---------------------------------------------------------------------------\u001b[0m",
      "\u001b[0;31mKeyboardInterrupt\u001b[0m                         Traceback (most recent call last)",
      "\u001b[0;32m<ipython-input-6-27ad5389acfd>\u001b[0m in \u001b[0;36m<module>\u001b[0;34m()\u001b[0m\n\u001b[1;32m     86\u001b[0m     \u001b[0max1\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mclear\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m     87\u001b[0m     \u001b[0max1\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mpcolormesh\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mimgdisp\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mcmap\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;34m\"gray\"\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mvmin\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;34m-\u001b[0m\u001b[0;36m1\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mvmax\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;36m1\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m---> 88\u001b[0;31m     \u001b[0mfig1\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mcanvas\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mdraw\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m     89\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m     90\u001b[0m     \u001b[0mtc\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mtc\u001b[0m \u001b[0;34m+\u001b[0m \u001b[0;36m1\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m/Library/Frameworks/Python.framework/Versions/3.6/lib/python3.6/site-packages/matplotlib/backends/backend_webagg_core.py\u001b[0m in \u001b[0;36mdraw\u001b[0;34m(self)\u001b[0m\n\u001b[1;32m    159\u001b[0m         \u001b[0mbackend_agg\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mRendererAgg\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mlock\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0macquire\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    160\u001b[0m         \u001b[0;32mtry\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 161\u001b[0;31m             \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mfigure\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mdraw\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mrenderer\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m    162\u001b[0m         \u001b[0;32mfinally\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    163\u001b[0m             \u001b[0mbackend_agg\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mRendererAgg\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mlock\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mrelease\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m/Library/Frameworks/Python.framework/Versions/3.6/lib/python3.6/site-packages/matplotlib/artist.py\u001b[0m in \u001b[0;36mdraw_wrapper\u001b[0;34m(artist, renderer, *args, **kwargs)\u001b[0m\n\u001b[1;32m     53\u001b[0m                 \u001b[0mrenderer\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mstart_filter\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m     54\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m---> 55\u001b[0;31m             \u001b[0;32mreturn\u001b[0m \u001b[0mdraw\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0martist\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mrenderer\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m*\u001b[0m\u001b[0margs\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m**\u001b[0m\u001b[0mkwargs\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m     56\u001b[0m         \u001b[0;32mfinally\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m     57\u001b[0m             \u001b[0;32mif\u001b[0m \u001b[0martist\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mget_agg_filter\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m \u001b[0;32mis\u001b[0m \u001b[0;32mnot\u001b[0m \u001b[0;32mNone\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m/Library/Frameworks/Python.framework/Versions/3.6/lib/python3.6/site-packages/matplotlib/figure.py\u001b[0m in \u001b[0;36mdraw\u001b[0;34m(self, renderer)\u001b[0m\n\u001b[1;32m   1297\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m   1298\u001b[0m             mimage._draw_list_compositing_images(\n\u001b[0;32m-> 1299\u001b[0;31m                 renderer, self, artists, self.suppressComposite)\n\u001b[0m\u001b[1;32m   1300\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m   1301\u001b[0m             \u001b[0mrenderer\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mclose_group\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m'figure'\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m/Library/Frameworks/Python.framework/Versions/3.6/lib/python3.6/site-packages/matplotlib/image.py\u001b[0m in \u001b[0;36m_draw_list_compositing_images\u001b[0;34m(renderer, parent, artists, suppress_composite)\u001b[0m\n\u001b[1;32m    136\u001b[0m     \u001b[0;32mif\u001b[0m \u001b[0mnot_composite\u001b[0m \u001b[0;32mor\u001b[0m \u001b[0;32mnot\u001b[0m \u001b[0mhas_images\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    137\u001b[0m         \u001b[0;32mfor\u001b[0m \u001b[0ma\u001b[0m \u001b[0;32min\u001b[0m \u001b[0martists\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 138\u001b[0;31m             \u001b[0ma\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mdraw\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mrenderer\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m    139\u001b[0m     \u001b[0;32melse\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    140\u001b[0m         \u001b[0;31m# Composite any adjacent images together\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m/Library/Frameworks/Python.framework/Versions/3.6/lib/python3.6/site-packages/matplotlib/artist.py\u001b[0m in \u001b[0;36mdraw_wrapper\u001b[0;34m(artist, renderer, *args, **kwargs)\u001b[0m\n\u001b[1;32m     53\u001b[0m                 \u001b[0mrenderer\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mstart_filter\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m     54\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m---> 55\u001b[0;31m             \u001b[0;32mreturn\u001b[0m \u001b[0mdraw\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0martist\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mrenderer\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m*\u001b[0m\u001b[0margs\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m**\u001b[0m\u001b[0mkwargs\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m     56\u001b[0m         \u001b[0;32mfinally\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m     57\u001b[0m             \u001b[0;32mif\u001b[0m \u001b[0martist\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mget_agg_filter\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m \u001b[0;32mis\u001b[0m \u001b[0;32mnot\u001b[0m \u001b[0;32mNone\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m/Library/Frameworks/Python.framework/Versions/3.6/lib/python3.6/site-packages/matplotlib/axes/_base.py\u001b[0m in \u001b[0;36mdraw\u001b[0;34m(self, renderer, inframe)\u001b[0m\n\u001b[1;32m   2435\u001b[0m             \u001b[0mrenderer\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mstop_rasterizing\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m   2436\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m-> 2437\u001b[0;31m         \u001b[0mmimage\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_draw_list_compositing_images\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mrenderer\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0martists\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m   2438\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m   2439\u001b[0m         \u001b[0mrenderer\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mclose_group\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m'axes'\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m/Library/Frameworks/Python.framework/Versions/3.6/lib/python3.6/site-packages/matplotlib/image.py\u001b[0m in \u001b[0;36m_draw_list_compositing_images\u001b[0;34m(renderer, parent, artists, suppress_composite)\u001b[0m\n\u001b[1;32m    136\u001b[0m     \u001b[0;32mif\u001b[0m \u001b[0mnot_composite\u001b[0m \u001b[0;32mor\u001b[0m \u001b[0;32mnot\u001b[0m \u001b[0mhas_images\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    137\u001b[0m         \u001b[0;32mfor\u001b[0m \u001b[0ma\u001b[0m \u001b[0;32min\u001b[0m \u001b[0martists\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 138\u001b[0;31m             \u001b[0ma\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mdraw\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mrenderer\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m    139\u001b[0m     \u001b[0;32melse\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    140\u001b[0m         \u001b[0;31m# Composite any adjacent images together\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m/Library/Frameworks/Python.framework/Versions/3.6/lib/python3.6/site-packages/matplotlib/artist.py\u001b[0m in \u001b[0;36mdraw_wrapper\u001b[0;34m(artist, renderer, *args, **kwargs)\u001b[0m\n\u001b[1;32m     53\u001b[0m                 \u001b[0mrenderer\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mstart_filter\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m     54\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m---> 55\u001b[0;31m             \u001b[0;32mreturn\u001b[0m \u001b[0mdraw\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0martist\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mrenderer\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m*\u001b[0m\u001b[0margs\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m**\u001b[0m\u001b[0mkwargs\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m     56\u001b[0m         \u001b[0;32mfinally\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m     57\u001b[0m             \u001b[0;32mif\u001b[0m \u001b[0martist\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mget_agg_filter\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m \u001b[0;32mis\u001b[0m \u001b[0;32mnot\u001b[0m \u001b[0;32mNone\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m/Library/Frameworks/Python.framework/Versions/3.6/lib/python3.6/site-packages/matplotlib/collections.py\u001b[0m in \u001b[0;36mdraw\u001b[0;34m(self, renderer)\u001b[0m\n\u001b[1;32m   1915\u001b[0m                 \u001b[0mgc\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mtransform\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mfrozen\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_meshWidth\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_meshHeight\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m   1916\u001b[0m                 \u001b[0mcoordinates\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0moffsets\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mtransOffset\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mget_facecolor\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m-> 1917\u001b[0;31m                 self._antialiased, self.get_edgecolors())\n\u001b[0m\u001b[1;32m   1918\u001b[0m         \u001b[0mgc\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mrestore\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m   1919\u001b[0m         \u001b[0mrenderer\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mclose_group\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m__class__\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m__name__\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;31mKeyboardInterrupt\u001b[0m: "
     ]
    }
   ],
   "source": [
    "#### Create FDTD object\n",
    "fs = fdtdVar(rows, columns)\n",
    "tc = 0\n",
    "\n",
    "if dflag == \"c\":\n",
    "    bflag = input(\"Enter 'c' for a closed domain or 'o' for an open domain: \")\n",
    "    print(\"\")\n",
    "    if bflag != \"c\" and bflag != \"o\":\n",
    "        sys.exit(\"Error: enter 'c' for closed domain or 'o' for open domain\")\n",
    "else:\n",
    "    bflag = \"c\"\n",
    "\n",
    "#### Wait to start wave propagation\n",
    "if vflag == \"w\":\n",
    "    print(\"Click on the camera window and press 's' to start wave propagation...\")\n",
    "    print(\"\")\n",
    "    while True:\n",
    "        # Reset domain\n",
    "        fs.update_domain()\n",
    "\n",
    "        # Capture frame\n",
    "        retval, frame = cap.read()\n",
    "\n",
    "        # Convert to grayscale\n",
    "        img = cv2.cvtColor(frame, cv2.COLOR_BGR2GRAY) / 256.0\n",
    "\n",
    "        # Create rigid material\n",
    "        imgtemp = np.pad(img, ((0, 0), (0, 1)), \"constant\", constant_values=1.0)\n",
    "        idx = imgtemp < 0.4\n",
    "        fs.mvx[idx] = 1\n",
    "        imgtemp = np.pad(img, ((0, 1), (0, 0)), \"constant\", constant_values=1.0)\n",
    "        idx = imgtemp < 0.4\n",
    "        fs.mvy[idx] = 1\n",
    "\n",
    "        # Display image\n",
    "        cv2.imshow(\"frame\", img)\n",
    "        if cv2.waitKey(1) & 0xFF == ord(\"s\"):\n",
    "            break\n",
    "elif vflag == \"i\":\n",
    "#     print(\"Click on the image window and press 's' to start wave propagation...\")\n",
    "#     print(\"\")\n",
    "#     while True:\n",
    "    # Reset domain\n",
    "    fs.update_domain()\n",
    "\n",
    "    # Capture frame\n",
    "    img = threshimg / 256.0\n",
    "        \n",
    "    # Create rigid material\n",
    "    imgtemp = np.pad(img, ((0, 0), (0, 1)), \"constant\", constant_values=1.0)\n",
    "    idx = imgtemp < 0.5\n",
    "    fs.mvx[idx] = 1\n",
    "    imgtemp = np.pad(img, ((0, 1), (0, 0)), \"constant\", constant_values=1.0)\n",
    "    idx = imgtemp < 0.5\n",
    "    fs.mvy[idx] = 1\n",
    "\n",
    "#         # Display image\n",
    "#         cv2.imshow(\"frame\", img)\n",
    "#         if cv2.waitKey(1) & 0xFF == ord(\"s\"):\n",
    "#             break\n",
    "\n",
    "## Release the VideoCapture object\n",
    "if vflag == \"w\":\n",
    "    cap.release()\n",
    "    cv2.destroyAllWindows()\n",
    "\n",
    "# print(\"Press 'q' to quit (when the image window is topmost)...\")\n",
    "# print(\"\")\n",
    "\n",
    "fig1 = plt.figure()\n",
    "ax1 = fig1.add_subplot(111)\n",
    "fig1.show()\n",
    "fig1.canvas.draw()\n",
    "while True:\n",
    "    # Update image with FDTD solution\n",
    "    fs.fdtd_update()\n",
    "    fs.source(tc)\n",
    "    if bflag == \"o\":\n",
    "        fs.boundary()\n",
    "    imgdisp = img + fs.pr\n",
    "\n",
    "#     # Display image\n",
    "#     cv2.imshow(\"frame\", imgdisp)\n",
    "#     if cv2.waitKey(1) & 0xFF == ord(\"q\"):\n",
    "#         break\n",
    "    ax1.clear()\n",
    "    ax1.pcolormesh(imgdisp, cmap=\"gray\", vmin=-1, vmax=1)\n",
    "    fig1.canvas.draw()\n",
    "\n",
    "    tc = tc + 1\n",
    "\n",
    "fig3 = plt.figure()\n",
    "ax3 = fig3.add_subplot(111)\n",
    "ax3.imshow(imgdisp, cmap=\"gray\")\n",
    "fig3.tight_layout()\n",
    "fig3.show()\n",
    "fig3.canvas.draw()\n",
    "print(\"If the image window is not responding: Run the cell below (or restart kernel)\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "scrolled": true
   },
   "outputs": [],
   "source": [
    "import os\n",
    "os._exit(00)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "scrolled": false
   },
   "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.2"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}