12_Designing_Nonlinear_Kalman_Filters.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "[Table of Contents](http://nbviewer.ipython.org/github/rlabbe/Kalman-and-Bayesian-Filters-in-Python/blob/master/table_of_contents.ipynb)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Designing Nonlinear Kalman Filters"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "<style>\n",
       "@import url('http://fonts.googleapis.com/css?family=Source+Code+Pro');\n",
       "@import url('http://fonts.googleapis.com/css?family=Vollkorn');\n",
       "@import url('http://fonts.googleapis.com/css?family=Arimo');\n",
       "\n",
       "    div.cell{\n",
       "        width: 850px;\n",
       "        margin-left: 0% !important;\n",
       "        margin-right: auto;\n",
       "    }\n",
       "    div.text_cell code {\n",
       "        background: transparent;\n",
       "        color: #000000;\n",
       "        font-weight: 600;\n",
       "        font-size: 11pt;\n",
       "        font-style: bold;\n",
       "        font-family:  'Source Code Pro', Consolas, monocco, monospace;\n",
       "   }\n",
       "    h1 {\n",
       "        font-family: 'Open sans',verdana,arial,sans-serif;\n",
       "\t}\n",
       "\t\n",
       "    div.input_area {\n",
       "        background: #F6F6F9;\n",
       "        border: 1px solid #586e75;\n",
       "    }\n",
       "\n",
       "    .text_cell_render h1 {\n",
       "        font-weight: 200;\n",
       "        font-size: 30pt;\n",
       "        line-height: 100%;\n",
       "        color:#c76c0c;\n",
       "        margin-bottom: 0.5em;\n",
       "        margin-top: 1em;\n",
       "        display: block;\n",
       "        white-space: wrap;\n",
       "    } \n",
       "    h2 {\n",
       "        font-family: 'Open sans',verdana,arial,sans-serif;\n",
       "    }\n",
       "    .text_cell_render h2 {\n",
       "        font-weight: 200;\n",
       "        font-size: 20pt;\n",
       "        font-style: italic;\n",
       "        line-height: 100%;\n",
       "        color:#c76c0c;\n",
       "        margin-bottom: 0.5em;\n",
       "        margin-top: 1.5em;\n",
       "        display: block;\n",
       "        white-space: nowrap;\n",
       "    } \n",
       "    h3 {\n",
       "        font-family: 'Open sans',verdana,arial,sans-serif;\n",
       "    }\n",
       "    .text_cell_render h3 {\n",
       "        font-weight: 300;\n",
       "        font-size: 18pt;\n",
       "        line-height: 100%;\n",
       "        color:#d77c0c;\n",
       "        margin-bottom: 0.5em;\n",
       "        margin-top: 2em;\n",
       "        display: block;\n",
       "        white-space: nowrap;\n",
       "    }\n",
       "    h4 {\n",
       "        font-family: 'Open sans',verdana,arial,sans-serif;\n",
       "    }\n",
       "    .text_cell_render h4 {\n",
       "        font-weight: 300;\n",
       "        font-size: 16pt;\n",
       "        color:#d77c0c;\n",
       "        margin-bottom: 0.5em;\n",
       "        margin-top: 0.5em;\n",
       "        display: block;\n",
       "        white-space: nowrap;\n",
       "    }\n",
       "    h5 {\n",
       "        font-family: 'Open sans',verdana,arial,sans-serif;\n",
       "    }\n",
       "    .text_cell_render h5 {\n",
       "        font-weight: 300;\n",
       "        font-style: normal;\n",
       "        color: #1d3b84;\n",
       "        font-size: 16pt;\n",
       "        margin-bottom: 0em;\n",
       "        margin-top: 1.5em;\n",
       "        display: block;\n",
       "        white-space: nowrap;\n",
       "    }\n",
       "    div.text_cell_render{\n",
       "        font-family: 'Arimo',verdana,arial,sans-serif;\n",
       "        line-height: 135%;\n",
       "        font-size: 125%;\n",
       "        width:750px;\n",
       "        margin-left:auto;\n",
       "        margin-right:auto;\n",
       "        text-align:justify;\n",
       "        text-justify:inter-word;\n",
       "    }\n",
       "    div.output_subarea.output_text.output_pyout {\n",
       "        overflow-x: auto;\n",
       "        overflow-y: scroll;\n",
       "        max-height: 50000px;\n",
       "    }\n",
       "    div.output_subarea.output_stream.output_stdout.output_text {\n",
       "        overflow-x: auto;\n",
       "        overflow-y: scroll;\n",
       "        max-height: 50000px;\n",
       "    }\n",
       "    code{\n",
       "      font-size: 70%;\n",
       "    }\n",
       "    .rendered_html code{\n",
       "    background-color: transparent;\n",
       "    }\n",
       "    ul{\n",
       "        margin: 2em;\n",
       "    }\n",
       "    ul li{\n",
       "        padding-left: 0.5em; \n",
       "        margin-bottom: 0.5em; \n",
       "        margin-top: 0.5em; \n",
       "    }\n",
       "    ul li li{\n",
       "        padding-left: 0.2em; \n",
       "        margin-bottom: 0.2em; \n",
       "        margin-top: 0.2em; \n",
       "    }\n",
       "    ol{\n",
       "        margin: 2em;\n",
       "    }\n",
       "    ol li{\n",
       "        padding-left: 0.5em; \n",
       "        margin-bottom: 0.5em; \n",
       "        margin-top: 0.5em; \n",
       "    }\n",
       "    ul li{\n",
       "        padding-left: 0.5em; \n",
       "        margin-bottom: 0.5em; \n",
       "        margin-top: 0.2em; \n",
       "    }\n",
       "    a:link{\n",
       "       font-weight: bold;\n",
       "       color:#447adb;\n",
       "    }\n",
       "    a:visited{\n",
       "       font-weight: bold;\n",
       "       color: #1d3b84;\n",
       "    }\n",
       "    a:hover{\n",
       "       font-weight: bold;\n",
       "       color: #1d3b84;\n",
       "    }\n",
       "    a:focus{\n",
       "       font-weight: bold;\n",
       "       color:#447adb;\n",
       "    }\n",
       "    a:active{\n",
       "       font-weight: bold;\n",
       "       color:#447adb;\n",
       "    }\n",
       "    .rendered_html :link {\n",
       "       text-decoration: underline; \n",
       "    }\n",
       "    .rendered_html :hover {\n",
       "       text-decoration: none; \n",
       "    }\n",
       "    .rendered_html :visited {\n",
       "      text-decoration: none;\n",
       "    }\n",
       "    .rendered_html :focus {\n",
       "      text-decoration: none;\n",
       "    }\n",
       "    .rendered_html :active {\n",
       "      text-decoration: none;\n",
       "    }\n",
       "    .warning{\n",
       "        color: rgb( 240, 20, 20 )\n",
       "    } \n",
       "    hr {\n",
       "      color: #f3f3f3;\n",
       "      background-color: #f3f3f3;\n",
       "      height: 1px;\n",
       "    }\n",
       "    blockquote{\n",
       "      display:block;\n",
       "      background: #fcfcfc;\n",
       "      border-left: 5px solid #c76c0c;\n",
       "      font-family: 'Open sans',verdana,arial,sans-serif;\n",
       "      width:680px;\n",
       "      padding: 10px 10px 10px 10px;\n",
       "      text-align:justify;\n",
       "      text-justify:inter-word;\n",
       "      }\n",
       "      blockquote p {\n",
       "        margin-bottom: 0;\n",
       "        line-height: 125%;\n",
       "        font-size: 100%;\n",
       "      }\n",
       "</style>\n",
       "<script>\n",
       "    MathJax.Hub.Config({\n",
       "                        TeX: {\n",
       "                           extensions: [\"AMSmath.js\"]\n",
       "                           },\n",
       "                tex2jax: {\n",
       "                    inlineMath: [ ['$','$'], [\"\\\\(\",\"\\\\)\"] ],\n",
       "                    displayMath: [ ['$$','$$'], [\"\\\\[\",\"\\\\]\"] ]\n",
       "                },\n",
       "                displayAlign: 'center', // Change this to 'center' to center equations.\n",
       "                \"HTML-CSS\": {\n",
       "                    styles: {'.MathJax_Display': {\"margin\": 4}}\n",
       "                }\n",
       "        });\n",
       "</script>\n"
      ],
      "text/plain": [
       "<IPython.core.display.HTML object>"
      ]
     },
     "execution_count": 1,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "#format the book\n",
    "%matplotlib inline\n",
    "from __future__ import division, print_function\n",
    "import matplotlib.pyplot as plt\n",
    "import book_format\n",
    "book_format.load_style()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Introduction"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "** Author's note: I was initially planning to have a design nonlinear chapter that compares various approaches. This may or may not happen, but for now this chapter has no useful content and I suggest not reading it. **"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We see that the Kalman filter reasonably tracks the ball. However, as already explained, this is a silly example; we can predict trajectories in a vacuum with arbitrary precision; using a Kalman filter in this example is a needless complication."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Kalman Filter with Air Drag"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "I will dispense with the step 1, step 2, type approach and proceed in a more natural style that you would use in a non-toy engineering problem. We have already developed a Kalman filter that does excellently at tracking a ball in a vacuum, but that does not incorporate the effects of air drag into the model. We know that the process model is implemented with $\\textbf{F}$, so we will turn our attention to that immediately.\n",
    "\n",
    "Notionally, the computation that $\\textbf{F}$ computes is\n",
    "\n",
    "$$x' = Fx$$\n",
    "\n",
    "With no air drag, we had\n",
    "\n",
    "$$\n",
    "\\mathbf{F} = \\begin{bmatrix}\n",
    "1 & \\Delta t & 0 & 0 & 0 \\\\\n",
    "0 & 1 & 0 & 0 & 0 \\\\\n",
    "0 & 0 & 1 & \\Delta t & \\frac{1}{2}{\\Delta t}^2 \\\\\n",
    "0 & 0 & 0 & 1 & \\Delta t \\\\\n",
    "0 & 0 & 0 & 0 & 1\n",
    "\\end{bmatrix}\n",
    "$$"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "which corresponds to the equations\n",
    "\n",
    "$$ \n",
    "\\begin{aligned}\n",
    "x &= x + v_x \\Delta t \\\\\n",
    "v_x &= v_x \\\\\n",
    "\\\\\n",
    "y &= y + v_y \\Delta t + \\frac{a_y}{2} {\\Delta t}^2 \\\\\n",
    "v_y &= v_y + a_y \\Delta t \\\\\n",
    "a_y &= a_y\n",
    "\\end{aligned}\n",
    "$$"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "From the section above we know that our new Euler equations must be\n",
    "\n",
    "$$ \n",
    "\\begin{aligned}\n",
    "x &= x + v_x \\Delta t \\\\\n",
    "v_x &= v_x \\\\\n",
    "\\\\\n",
    "y &= y + v_y \\Delta t + \\frac{a_y}{2} {\\Delta t}^2 \\\\\n",
    "v_y &= v_y + a_y \\Delta t \\\\\n",
    "a_y &= a_y\n",
    "\\end{aligned}\n",
    "$$"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Realistic 2D Position Sensors"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The position sensor in the last example are not very realistic. In general there is no 'raw' sensor that provides (x,y) coordinates. We have GPS, but GPS already uses a Kalman filter to create a filtered output; we should not be able to improve the signal by passing it through another Kalman filter unless we incorporate additional sensors to provide additional information. We will tackle that problem later. \n",
    "\n",
    "Consider the following set up. In an open field we put two transmitters at a known location, each transmitting a signal that we can detect. We process the signal and determine  how far we are from that signal, with some noise. First, let's look at a visual depiction of that."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "image/png": [
       "iVBORw0KGgoAAAANSUhEUgAAAt0AAAFpCAYAAABAhJQMAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\n",
       "AAALEgAACxIB0t1+/AAAIABJREFUeJzs3XeUHFeZNvCnqnPP9OSRNEFhNBrlnCUr2ZLlnDBOOGHA\n",
       "NrbxAl4OLCzLsizhY1lYPgNr+ExwNs42DshB0ZKsaGVplOOMwuTUuau+P4TkGVWNJt3q6qp+fuf4\n",
       "HM7tnluvYbj1zr3vvVdSVVUFEREREREZRjY7ACIiIiIiu2PSTURERERkMCbdREREREQGY9JNRERE\n",
       "RGQwJt1ERERERAZzJvuBTU1NyX4kEREREZFQ2dnZPfo+Z7qJiIiIiAzGpJuIiIiIyGBJLy9pr6fT\n",
       "8mQt8YSCJRsO4KXlu9DYGjY7nF5rbGwEAOTk5PS5L1mSkBvwIjfgRV7Ahyy/Bz6PC36vC163Ez63\n",
       "Ez7P2f/s97rgdjogyxJkSYIkne1DVQFVVaGoKuIJBeFoHKFIHMFI7B//OYZQJI5wNI7WUBQNrSHU\n",
       "N4dQ3xJGOBrv87+DXcyfMBh3XT4eA/IyzQ5FqE2bNgEApk6danIkqSl+8gRCG1cjsvNTqLFYt39O\n",
       "cnvgHjYS7uFj4a4YBdmfYWCU6aMvv6+qqiJRcxrRfTsR3bsLsRNHevTzjtx8eKfMgnfSDMh+e40D\n",
       "ZJy+lEmbmnSTPSmKilXbj+L5j7bjVH2b2eEklcspozg/gJKCAEoKstA/NwN5WT7kZ/mRG/AiO8ML\n",
       "WZZMiy8Uif0jAT+biNc2BVFV24Kq2mZU17Va+o+jnlq57ShW7ziGK6cPw+2XjUVOptfskMggaiyG\n",
       "yO5tCG9a06PETHK54Rk1Hp6xk+Eqq4Dk5CszlUiSBGe/AXD2GwD/nEVQWlsQqdyByI7NiB071OXP\n",
       "Jxrq0PbROwguXwLPmInwTrsEzpLBkCTzxmiyNynZ18C3/wuBM932oqoqNu87iWfe34bDpxrNDkcY\n",
       "vZluj8uB8uJclBXlorQwCyUFARTnB1CYk2FqUt1XraEoqs8l4bUtOHq6CQeq6lHTFDQ7NEN53U7c\n",
       "eMkI3DR3FPxel9nh9Alnuj+TaGpEeONqhLesgxLs3gSAJMtwlY+AZ/xUeIaPgeT2GBxlejPq9zXR\n",
       "UIfIzk8R3r4ZidrT3f45Z1EpfNMugWfcVP6RRbr6kscy6SYh9h6rxVPvb8XOwzVmhyJcsLUZpfl+\n",
       "zJkyGsNK8lBenIfSwixLJ9c91dQaxoGqehysbsCBqnrbJuJZfg9uXTAaV8+sgMvpMDucXmHSDcRr\n",
       "zyC0ZikiOzZDTSS69TOuksHwjJ8Cz+iJkDMDBkdI5xj9+6qqKhKnqhDevgmRnZ9CaW3p1s/JgWz4\n",
       "Zs6Hd8osyB6ugtFnmHSTaRpbw3hqyVYs/fSw2aEIk53hwbih/TCurD/GDCnE6WP7IctSWicxeppa\n",
       "w6g8Vovth05jx6EztlrdKCkI4MHrpmBSRZHZofRYOifdserjCK1eisiebd36vuTxwjthGrxTZ8NZ\n",
       "OMDg6EhPMn9fVUVB9MAehDetQXT/nm79jOz1wTtjHnzT57DumwAw6SYTJBIK3lu/H89/tANt4e5v\n",
       "RkpFWf5zSXY/jBvaHwP7ZXWo6UvnJKYnWoIR7Dx8BjsOncHOw/ZIwmePKcVXrpmMwhzrbJpLt99X\n",
       "VVURO3IAodUfIXpoX7d+xjmgBL5pc+AZO4nlIyYz6/c1UV+L0Oa1iGxZDyXU9aqd5HLBO3kWfLMu\n",
       "hSO775vqybr6kseyYIl6bPeRGvz+b5ssnVQNK8nFjFGlmD6yBEMG5KRVqYhRAn4PZo0ZiFljBgI4\n",
       "OxO+ed9JbKiswqf7TiJkwZNT1u46gc37TuLWBWNw09yRli05savooX0ILv97tzZHSg4nN8vReY68\n",
       "AmRefj0yFlyFyO6tCG9cg1jV0U6/r8ZiCK1fhfCmNfCMnwr//CuZfFOPMemmbmtoCeGpJVuxbMsR\n",
       "s0PpMZdTxoTy/pg+sgTTRpagINtvdki2l53pxWWTy3DZ5DLE4gnsOHQGGyqrsH5PFWotVA8eiSXw\n",
       "7IfbsfTTQ3jg2imYMqLY7JDSXqz6OIJL30X00N4uvyv7/PBOnwvftDmQM1geQB1JLtfZEqMJ0xCr\n",
       "OobQ2mWI7O68PElNJBDesh6RHZvhnT4X/ksW8vhI6jaWl1CXFEXFkg0H8PT72xCMWKeUxO10YMao\n",
       "EswZNwiTKgbA5+ndqRTptlxvNFVVcfhkI9buOo4VW4/gdIO1jpWcPaYUD143FXlZPrND0WXn39dE\n",
       "XQ3alr+HyK6tXX5XDmTDP2sBPJNnciNcCkvF39d47Zmzyff2TV1uxJW8XvhnL4Rv5jxILneSIiQz\n",
       "saabDHO6vhWPv74e2w+dMTuUbpEkYPzQ/lgwcQhmjS5Fhq/vg2AqvhTsQlVVVB6rxfItR7B6xzG0\n",
       "hKJmh9QtmT43HrxuCuZPSL0yBTv+viZamhBc+QEiW9ZBVZSLfteRVwD/nIU88s0iUvn3NdHUiNC6\n",
       "FQhvXtvlRUpyIAv+eVfAO2k6JAd/7+yMSTcJd252+y9LtlriFsOyATm4dNIQzBs/GPmCS0dS+aVg\n",
       "J/GEgs17q7Fi2xGs31OFWPziyVUqmDGqBA/fMC2lZr3t9PuqxqIIrlmG0NplXSY9zv7F8M+9HO5R\n",
       "4yHJcpIipL6ywu+rEmxDaMPHCK9fBSUcuuh3HXkFyLjiRniGj0lSdJRs3EhJQllldtvjcmDBxCG4\n",
       "ekYFhhbnmh0O9ZHTIWPG6FLMGF2KtlAUy7Ycxnvr9+NETffO1TXD+j1V2HWkJmVnva1KVVVE9+5E\n",
       "25I3kGhquOh3HbkF8F92FTyjJzLZJkPI/gxkLLgSvhnzEFqzDKH1K6HG9SejEvW1aH7xj3APH4PM\n",
       "K26EI68gydFSKmPSTedZZXa7tDCAq2dU4LJJZULKRyj1ZPjcuG72CFw7azh2HDqD99bvx7rdJ5BQ\n",
       "krow1y2toSh++fInWL3jWMrNeltRoq4GrUteR/RA5UW/J2cG4J+3GN7JM7mcT0kh+/zIWHQtvNPn\n",
       "ILjqA0S2rO+03Cm6bxcaDu6Fb85C+C9ZCMll7ZtuSQyOVAQAqG8O4VevfIJtB7t/XW4yOWQJM0eX\n",
       "4uoZFRg3tB9nFNOEJEkYX94f48v7o64piPc3HsT7Gw+ivuXiS7xmODfr/ehN0zF77ECzw7EcNRpB\n",
       "cPVHCK1dATXR+R/9kscL/yWXwTdjHs/YJlM4snIQuPZW+GYuQHD5e52edqIm4giufB+R7ZuQsfhG\n",
       "uEeM4bsrzbGmm7B5bzX+59V1aGqLmB2KhsflwJXTh+HGOSNNO+bPCjWH6SSeUPDx9qN4bdUeHD3d\n",
       "1PUPmOCamRX40lWT4HYl/1xvK/6+Rip3oG3J60g0dX72v+RwwDttDvxzL+cRbTZixd/XC8WqjqHt\n",
       "o7cRO3Lgot9zV4xC5pWfY8mJxbGmm3olnlDw3Ifb8dqq7l2Hm0yZPjeumzUc184ajqwMzmbRZ5wO\n",
       "GZdOKsP8CUOwsbIKr67ajcpjdWaH1cG76/Zj95EafOeOS1BSmGV2OClLCbai9b3XujwC0D10BDKu\n",
       "+hycBf2SFBlR97lKBiH7nocR3bUVrR+8BaVFfzIgun8PGo78FzIWXQfv1Eu4ByENMelOU6frW/Hf\n",
       "L69NuWQlP8uHG+eMxBXTynt9rjalB1mWMGN0KaaPKsHOw2fwyord2HLglNlhnXf4VCO+8bv38fAN\n",
       "U3HppDKzw0k5kd3b0Pruq1CCrZ1+x5Gdc3ZZftR4LstTSpMkCZ6xk+CqGIXQqg8RWrdCt95bjcXQ\n",
       "+vfXEdmzHYHrbuOsd5ph0p2GPtl1HI+/vgGtKXQmcmG2H7ddOgaXTS7jVdvUI5IkYdzQ/hg3tD8O\n",
       "VNXjr8t2Yv2eKrPDAgCEo3H86pV12H7wNB68fiq8bg65SrD1bNKxc0un35EcTvhmL4B/ziLWbZOl\n",
       "yB4vMi6/Dp5J09H299cRPbRP93uxIwfQ8Ptf/GPWezZnvdME3wBpJJ5Q8Kd3P8U76/abHcp5AZ8b\n",
       "t146BlfPqDCl/pXsZVhJHr5/9zzsOVqDp9/fhl1HaswOCQDw0aeHUXm8Ft/9wlwM6p++e1kiu7eh\n",
       "9b1XobR1PrvtLh+JjCtvYikJWZqzoD+y7voqonu2ofX9t6A0a/crqLEoWv/+GiJ7tnHWO004fvjD\n",
       "H/7QiI5/9rOfYf78+aitrcXVV199vj0S+WyzntfLq3mTpak1jB89sxJrdh43OxQAZzdIfn7+aHz7\n",
       "9kswvrw/HI7U/Su/uroaAFBcXGxyJNRdhTkZWDi5DBWleTh6uhGNreZvEm4ORrF8y2EM6peNUgPr\n",
       "vFPx91UJh9Dy1l8RXPF3qDH9FTbZ60PmdbchY9F1cGRkJjlCMksq/r6KIkkSnIUD4J08E2oohPhJ\n",
       "/fev0liP8Jb1kHwZcBaVspQqxfUljzVkpnvdunV48sknMX486/BSweGTDfjxs6twpjFodihwyBKu\n",
       "nD4Mt106BrkBnmdMxpEkCdNGlmDK8GKs3HYEz3+0A6cb2kyNKRSN4yfPf4y7Fo3HLQtGp8X4GKs6\n",
       "hpbXnkGiofP9I+4RY5B5zS1wBNJ3FYDsS/Z4Ebj2FnhGjUfr23/VPaVHjUXR+u4riB3eh8zrboPs\n",
       "5fvRjoQn3U1NTbjrrrvwl7/8BQZNolMPrN15HL965RNEYgmzQ8G0EcW4/9rJKMoPmB0KpRFZlnDp\n",
       "pDLMGTcI73yyDy8u3YmQiZc/qSrw7IfbceRUI75+8wx4bFrnraoqQutWIrj0HagJ/fFH9vqQcdXn\n",
       "4Bk3JS3+AKH05i4fgZyHvoO2D99GePNa3e9Edm9D/OQJBG6+B66SQUmOkIwmfLR/4IEHcMstt2D+\n",
       "/PlI8hHg1I6iqPjrsp14cdlOs0PBgLwM3H/NFEwfVWJ2KJTGXE4Hbpo7CvPGD8ZflmzFym1HTY3n\n",
       "4x3HUF3Xgu/fPc+0M+iNogRb0fLmi4ju393pdzi7TemoO7PeiYY6NP3lcfgXXgvfzPn8g9RGhCbd\n",
       "Tz75JA4dOoQXXngBAPiLYpJQJIZfv7oOa3edMDUOt9OBz88fhZvnjeYmSUoZ+dl+fOu22bhiWjn+\n",
       "8PZmUy/YOVjdgMd+9z6+e+ccjBpcaFocIsWOHkTza892elax5PEi86qb4RnP2W1KX+dnvd9/E+Et\n",
       "6zWfq4kE2j54C7HD+xG48Q7Ifu5zsANhN1Lu3bsXc+fOxerVqzF8+HAAwIIFCzBu3Dj85je/Of+9\n",
       "9jf57N+fOqdo2EVTMIonP9iPqnpz67fHDsrBjTMGIT/A474odcUTClbvOYP3t1YjHDWvBMspS7h9\n",
       "bhmmlOebFkOfqSo8uzbDu30DpE5eK/H8fgheshgKZ7eJznMd3Q//+hWQOtlkrPgz0DbnCiQKi5Ic\n",
       "GempqKg4/597eiOlsKT7qaeewpe+9CU4HJ/NaCYSCUiSBIfDgba2NrhcLibdBjrTFMYf3t+L+lbz\n",
       "zt/O9rtxy+zBGDMox7QYiHqqKRjFa58cxY6jnV9DngzXTxuIS8cNMDWGXonF4F+3FO5jBzv9Snjk\n",
       "RIQnzgQcXPUiupDc0gT/6vfhrNc/5lSVZYSmz0e0fHSSI6MLpUTS3dTUhKqqzy6kUFUV9913H4YP\n",
       "H47vfe97GD169Pnv9TZY6ty+43X4j6dXojlo3tFoCyeX4SvXTEamz21aDEbYtGkTAGDq1KkmR0JG\n",
       "UlUVq7YdxR/e3owWEy+O+tzckbj3iomQ5d6VXiT79zXRWI/mv/4J8dPVup/LPj8yb7gDnhFjkxIP\n",
       "WQvH18+o8Tjalr6D0LqVnX7HN30uMhbfAIl/vJqmL3mssJru7OxszcP9fj9yc3PPJ9xkjM17q/Gz\n",
       "F1abdkJJXsCHr900DdNGcqMkWZckSZg/cQjGl/fH/761Eet2m3Or5esfV6KxNYxHPzcDzhQ+vx74\n",
       "R/32y091epW7a2AZAjffA0c2V76IuiI5nci84ka4Bg9D61svQAmHNN8JbfgYidrTCNx8D2R/hglR\n",
       "Ul8YOqJLksSNMgZbvuUw/vPZVaYl3JdNGoLffeNqJtxkG7kBH75351x869ZZpq3aLNtyBD95bhXC\n",
       "Jh5t2JXQ5k/Q9OwTnSbc/ksWIvveR5hwE/WQZ+RY5Hz123CVDNb9PHpoHxr/+GvEz5xKcmTUV8LK\n",
       "S7qL5SXivPHxHvz571tNeXZuphdfu2l6WhwDyOXP9NXQEsJv39iADZX6pRNGG16ah3+/dwGyMrq/\n",
       "Idno31c1EUfb+28itHGN7ueSy4XM62+Hd+xkQ55P9sLxtXNqLIbWd19BeNtG3c8ljxeBm+5k6VaS\n",
       "9SWPTe21S9KlqiqeXrLVtIR70rABePyfrkqLhJvSW27Ah+/fPQ8PXjcFLmfyh8t9J+rxnf/3Ieqa\n",
       "zL9NFgCUSBjNf/1Tpwm3nJWDnPv+iQk3kQCSy4XMG+5A5uIbAJ2qAfXc/x83rDYhOuoNJt0Wo6oq\n",
       "/vzeFry6ak/Sn+2QJdx7xQT88IsLkJPpTfrzicwgSRKunTUcv3jwchTnJ/+s3BM1Lfjuk0tRa3Li\n",
       "rQRb0fTsE4geqNT93DVwCHLvfwzOotIkR0ZkX5IkwTdrAbK/8AAkr/57t/Xvr6Ft2Xu8kNACmHRb\n",
       "iKqqePKdT/Hmmr1Jf3Zhth8/u38hPj9/dK9PVSCysvKSPPz6a1di/gT9OksjnaxvxXef/AhnGtqS\n",
       "/mzg7A15jX9+HPGqY7qfeyfNQPY9j0DODCQ5MqL04B42Ejlf/iYcBf10Pw9+/CFa334JasK8+wao\n",
       "a0y6LUJVVfzh7c14+5N9SX/2jFEl+L+PXmmbG/OIesvnceGfb52Ff/rcdHiSfMvqqfo2fPfJj3C6\n",
       "Xn/jolHip6vR+JfHkajTOT9YkpB5xY3IvO42SE6hFxwT0QWcBf2Q8+VvwF0xSvfz8Jb1aHnlKaix\n",
       "WJIjo+5i0m0Bqqri/729Ge+uS+5lQg5ZwleunoR/vWsuAn7eLEkEnF3uvXxqOX718BUoKUjuzO6Z\n",
       "xiD+9U/LUNOYnBnv6JEDaHzqN1BamjWfSU4nsm69D76Z83lKFVGSyF4fsm77MryTZuh+Htm78+yp\n",
       "QqHU2AdCHTHpTnGqquJP723BO0lOuDN9bvzwiwtww5yRfKES6RjUPxu/fGgxJlck9wbJ0w1t+N4f\n",
       "ja/xjlTuQPNzf4AaDms+k70+ZN/1VXhGjjM0BiLSkhwOZF53G/xzF+l+Hjt+GE1P/QaJZnNv2CUt\n",
       "Jt0pTFVVPLVkK95Kcg33wMIs/OrhxZg4zILXURMlUYbPjX+/dwFumjMyqc89Vd+Gf/3jUtQ3ay/P\n",
       "ECG889Ozy9QJ7TnhciAb2V/8GlyDyw15NhF1TZIkZFx2DTKvvEn38/iZU2h66ndINDHxTiVMulPY\n",
       "qyt34/WP9U8KMMr0kcX474cWoyifG6KIukOWJXzp6kn45udnJvVYweq6VvzgL8vRKvjK+vC2TWh5\n",
       "/TmoiqL5zFHQDzlf+ic4+xcLfSYR9Y5vxjwEbr5b91r4REMtmp7+LRINdSZERnqYdKeojzYfwjMf\n",
       "bE/qM2+ZPxr/etc8+L2upD6XyA4um1yGn31lIXKTeJzm0dNN+PGzqxAVdCNteOsGtLz5PKBz9Jir\n",
       "ZDBy7nsUjpw8Ic8iIjG8Yycj6477Ibm1e68SDXVnE+/6WhMiowsx6U5BG/ZU4bdvbEja81xOGd+6\n",
       "dRbuuWICjwMk6oMRgwrwP49cgfLi3KQ9c9eRGvzipTVIJLQz0z0R3rYRLW+9qPuZe+gIZN/zEGR/\n",
       "8s8pJ6KuuctHIPveRyD7MzSfJZoa0fTM7zjjnQKYdKeYPUdr8PMX1yChJOeQe7/HhR/euwDzJw5J\n",
       "yvOI7C4/24+ffmUhJpT3T9oz1+2uwhN/29TryzHC2zeh5c0XdD9zV4xC1u1f1p1FI6LU4SoeiOx7\n",
       "Hr5I4v2/SDQ1mBAZncOkO4UcO92EHz2zCtF4cg63z8n04qdfuQzjk5gcEKUDv9eFf793PuaOG5S0\n",
       "Z76/8SBe+GhHj38usnNL5wn38DHIuvVLkFwsOSOyAmf/4rMz3hnaValEYz2anubmSjMx6U4RtU1B\n",
       "/PtTK4RviurMgLwM/NeDi1BewvpMIiO4nA5867bZuHZmRdKe+dflu7Bmz5lufz+6fzda3nhOt4bb\n",
       "PWIMsm79Ii+9IbIYZ7+is4m3zg2xiYY6ND33BJSgObfbpjsm3SmgLRTFD/683PBzd88pL87FL77K\n",
       "E0qIjCbLEh64bgruvnx80p752rqj2HakvsvvxU4cQfMrT+meUuIePgZZt3wRkoMJN5EVOQsHIPvu\n",
       "h3X3YSRqz6D5hSehRiMmRJbemHSbTFFU/OKltTheo73xzQgTyvvjp19ZiJwknrBAlM4kScKtl47B\n",
       "ozdNh5yEi6ZUFXhh1WEcPtl57Wa85tTZl67OddHuYSOZcBPZgLPfgH9sgNbWeMeqjqL51ad1z+In\n",
       "4zDpNtnT72/F5n0nk/KsqSOK8IN75vNIQCITLJ5Wjm/dNguOJJwQFI0r+M9nVqGxVXubZKK5Ec3P\n",
       "/0H3mmjXkGFna7hZUkJkC87+xci+66uQPNqJtuj+PWj920u93oBNPcek20TLPj2ctMtvpo8sxvfu\n",
       "nAu3S3uAPhElx9zxg/Ht2y9JSuJd0xTE/3lhNeLtjhJUQkE0P/8H3Y1UzgElyLrty9w0SWQzzqLS\n",
       "sycQ6axehbdvQvCjd0yIKj0x6TbJ3mO1+O2byTmLe+boEnz3zrlwOZlwE5lt9tiB+JcvzElK4r3r\n",
       "SA3+8I+jBNVYDM0v/hHxM6c033Pk5iP7Cw9A9rLsjMiO3EOGIfC5uwCdErfg2mUIfrIi+UGlISbd\n",
       "JqhrCuKnz69GLN63yyy6Y8aoEnznjjlwOvg/NVGqmDm6FN+5Izkz3ks2HsTf1+1Dy+vPInb8sOZz\n",
       "OSMT2Xd9FXIgy/BYiMg8ntETkHn1zbqftX3wFiI7tyQ5ovTDTCzJorEEfvr8x6hvCRn+rGkjivGd\n",
       "Oy5hwk2UgmaNGYhv3TY7KYn3p88+j5rNmzTtktuD7DsfhCOvwPAYiMh8vqmXwD//Ct3PWv72ImLV\n",
       "x5McUXphNpZEqqrid29uwL4TXR/n1VeThg3Av3xhDktKiFLYnHGD8Ngts/RWfIUpbz2ByQ2VOFBV\n",
       "j2jss5MKJIcDWbfeB2dRqXEPJ6KU459/BbxTZmva1VgMzS/9GUpLck5TS0dMupNo6aeHsWzLEcOf\n",
       "M7w0D9+9cw43TRJZwLwJg/HAtVMM6Ts/0ohLaz8FACQUBQerG86fVJB54xfgLh9hyHOJKHVJkoTM\n",
       "q2+GZ+Q4zWdKcyOaX/4L1DiPEjQCk+4kOXa6Cb//m3Z5V7SSggB+cM98+Dw8gYDIKq6dNRy3XTpG\n",
       "aJ++RBhXnV4Hl5o439YaiqKqpvnsTNfYyUKfR0TWIckyAjfdCeeAEs1nsRNH0PruKzxK0ABMupMg\n",
       "Gkvgv/66BpFYousv90FewIf/+OICZPPiGyLLuXPROCyeOlRIX7Kq4IrTGxBIaPeOrAlnobL/WCHP\n",
       "ISLrktweZN32Jd1bK8NbNyC8fpUJUdkbk+4kePKdzTh6usnQZ/g9LvzHfQvQP0/7fx4iSn2SJOHh\n",
       "G6ZhxijtzFOPqCoWtexBUaRO81GtOxvLCifjf17fgPpm4zdzE1Fqc+TkIeu2+yA5tOWorR+8hejB\n",
       "vSZEZV9Mug328fajWLLxoKHPcDll/Ns98zBkQI6hzyEiYzkcMr59+yUYPbj3p4mMC1dhfLhK0x6S\n",
       "3VjSbwZisgtNbRH88uW1UBQuHxOlO9egoci4SucoQVVFy2vP6F6mRb3DpNtAJ+ta8Ns3Nhr+nG/d\n",
       "Ohtjy/oZ/hwiMp7b5cC/3TMfAwt7fm52XrQZC1u1M1MKJLzfbwZaXBnn27YfOoOXV+zqU6xEZA++\n",
       "KbPgmz5H066Egmh57RmoCWPLY9MFk26DxBMKfvHXtQhGYoY+5+7Lx2P22IGGPoOIkivT58a/3TMP\n",
       "mT53t3/GqcQx6Pg2qIoCqCr88c/KRz7On4CTPu3s+YtLd2Ln4TNCYiYia8tYfANcQ4Zp2mPHDyO4\n",
       "8n0TIrIfJt0GeWXFLuyvMvY87rnjBuGWBaMNfQYRmaMoP4B/6eatlYqqonVfJX5yph+eDhajON6E\n",
       "8uBJDAjXoTJjIHZnlXX6c//3tXUIGTw5QESpT3I4kXXLvZAD2ZrPgh9/yPpuAZh0G+DwyQa8tNzY\n",
       "Zdvy4lx8/eYZkIy8VYOITDVh2AB85ZqLH+0XjsaxfcdBfNTghgoJbVEFA+ItAIDCaBOmNO1DbrTz\n",
       "yy5O1bfh6fe3CY2biKxJ9mci6+a7oXdjV8ubz/PinD5i0i1YPKHg16+uQ8LADUo5mV58/+558Lid\n",
       "hj2DiFLDNTMrOj1KsKaxDR9vO4ITbQoAwIMEfuHpmEAPDZ7Ed/c9i7HNhzp9xrvr9mPHodPigiYi\n",
       "y3INLkfGgis17UprC1refP5sCRv1CpNuwV5ZsQuHThq309fpkPG9O+egINtv2DOIKHVIkoSHbpjW\n",
       "4UQTRVWx91gt1u+pQiT+2QswAge+HZ2INnT8gzwzEcbDh9/AjdWrIKv6G6Ief309y0yICADgm7MI\n",
       "7rIKTXv00D6E1iw1ISJ7YNItUDLKSh6+YSpGDS409BlElFqcDhnfvXMuCrP9CEfjWLf7hO6ekYnu\n",
       "VtyTdQZLssZiZ2CI5vPFNRvxzQMv65absMyEiM45e2PlXZAztHd/BFcsQaz6uAlRWR+TbkGSUVay\n",
       "cHIZLp9ablj/RJS6cjK9mD22FKu2H9VcbCNDxV2ZZ/Cj3KOAy423syfgibLP4Y2iuUigY21mebC6\n",
       "03ITlpkQ0TlyIAuBm+7StKuKgta3XoCaiJsQlbUx6RbE6LKS0sIAHrxuimH9E1HqiicU/Nufl+GB\n",
       "X76DaKxjeUiuHMOPc4/gjswaQJLxdtZ4RGUnVEnCh/2m49flt6LB1XG26mLlJiwzIaJz3OUj4J+z\n",
       "UNMeP3MKwZUfmBCRtTHpFuBETTNeXrHbsP5dThnfuWMOfB6XYc8gotRUXduCRf/8DH787MdQL1hI\n",
       "m+huxW9Scdr4AAAgAElEQVTyD2KCJwgA2JQ7EqddHS/VOZhZip8Ovxs7A9pjA/XKTU7Vt+HFpTvF\n",
       "/4sQkSX5F1wJZ/9iTXtozVKWmfQQk+4+UlUVf3h7E+IJ43bzPnDtFF7xTpSGPth4EBPv/z1Wbjva\n",
       "oV2GinsCNfhR7lHkOs7OVJ9x52BLtnbjEwC0Of14ouymbpeb/G3tXhw73ST434aIrEhyOBG44Q5I\n",
       "cseU8XyZSZxlJt3FpLuP1u48jq0HjKuBnDNuIK6YxjpuonRyrpzkyu88h5rGYIfP+nuBVye04Qcj\n",
       "E3D8I39OQMbywilQpc6H9J6UmySUs5MJ6oVT60SUlpxFpfDNvVzTHj9zCsFVLDPpLibdfRCKxPDH\n",
       "97YY1v+AvAx87cbpvACHKI1crJzk0rIAlk5uwpzcOLIyvCjODwA4W1ZS787S6U2ru+Um2w+dwcfb\n",
       "j/X534eI7ME/dxHLTPqISXcfvLx8F2qbgl1/sRdkScI/3zobGT63If0TUerptJxElvDDW6bghcEn\n",
       "0M/9WSZeXBBAPL+o07KSznS33OTPf9/CTZVEBKCLMpO3X+KlOd3ApLuXqmqa8eaavYb1f+OcERg5\n",
       "qKDrLxKR5V2snGRAXiY+/MVd+Lr3EBzoOPUtO12Y9PAjcLp6fjttd8pN5la+j5c+3NrzfyEisqVO\n",
       "y0xOVSG8aa0JEVkLk+5eOLt5crNhmydLCwO4c9F4Q/omotRysXKSRVOGYuuTD2K2ow6x44c1P+uf\n",
       "vxgDRw/HnQvH9fr5XZWblD3znzi+Z3+v+ycie+mszCS4/D0oba0mRGQdTLp7Yd3uE9hy4JQhfcuS\n",
       "hK/fPBNul8OQ/okodVysnORH9y3Akp/fiUK/A20fva35WUdBP/hmXwoAuHHOSAwvzet1HBcrNxna\n",
       "Vo3w976M4IaPe90/EdmH5HAi85pbNO1KOIS2pe+YEJF1MOnuoXhCwVNLjFtuZVkJkf11VU7y0X/f\n",
       "jX+7Zz4cDhnBFUugtLZo+si88nOQHGfLShwOGd/4/Ey4nL0f0i9WbuKNBlH3H99E458f5/FgRATX\n",
       "wCHwTpyuaQ9vWY/YiSPJD8gimHT30EebD6G6zpjlE5aVENlfd8pJLp10ttQjfroa4Y2rNX14Rk+A\n",
       "u3xEh7aB/bL7VGZyzsXKTVpeewZn/uUBxGuMWekjIuvIWHgtJK9X09763mvcVNkJJt09EInGDbup\n",
       "TZLAshIim+tOOUn/vLOzzKqq6r68JJcLGYtv1O3/xjkjMagwo89xXqzcJLpnO04/eidCG7R/DBBR\n",
       "+pAzA8i49GpNe/zkCYQ//cSEiFIfk+4e+NvavahvCRnS95XThrGshMimelJOck5k1xbEjh26sCv4\n",
       "5y2GI1v/hlqHQ8Yts4dAFnC2/8XKTZSWJtT+xzdYbkKU5rxTZutvqlz6LpRgmwkRpTYm3d3UEozg\n",
       "tVV7DOk74HPj7sUsKyGyo56Uk5yjJuIILntP05cjvxC+WQsu+rzSfD9mjyzsa9jnsdyEiDojORzI\n",
       "vPrzmnYlHEJozTITIkptTLq76dWVu9EWNuaSiHuumICA32NI30Rknp6Uk7QX3rwOiYY6TXvmVZ9t\n",
       "nryYqyaXIEvgmMJyEyLqjGtQGbzjp2raQxtWIdHcaEJEqYtJdzfUNgXxzifGnFNbUZKHxVPLDemb\n",
       "iMzRm3KSc5RIGMFV72va3cNGwl0+slvP93uc+OKVE3oXfCfal5tEMnM7fMZyE6L05l94DSRnxwkB\n",
       "NR5HcKV2LEtnTLq74a/LdiIaTwjvV5KAr14/FbLc9/pLIkoNvSknaS+8bqXuBRMZl13TozgWTh6K\n",
       "EQPze/Qz3XEwsxQ/qbgbzkkzNZ+x3IQoPTmycuCbPk/THtm6AfHaMyZElJqYdHehtimIpZ9qb4IT\n",
       "YfHUcgw34KVIRObobTnJOUqwFcFPVmjaPWMnw1lU2qNYZFnCV6+fCgF7KjVqFReWz78f2V/8GiB3\n",
       "PHGJ5SZE6cl3yWWQvb4Obaqi6O5PSVdMurvw5upKQ657z/S5cc9iscu/RGSOvpSTtBf8+COokXCH\n",
       "NkmWkXHpVb2Ka1hJHq6cNqxXP9uVt9cdgOv6O9Hv//wBjvx+HT5juQlR+pH9GfBdcpmmPbJnG2JV\n",
       "x0yIKPUw6b6I5rYIlmw4YEjfn583ClkZ3DxJZHV9LSc5J9HUgPDGNZp275RZcOT1/jjR2y8bC48B\n",
       "5/+3hqJ4f+NBeMZMRP/fvADv1Nma77DchCi9+GbMgxzI0rQHeT08ACbdF/XOJ/sQiYmv5c7P8uHa\n",
       "WcOF90tEydXXcpL2gh9/BDXRcVZYcrngn7u4TzHmZflwwyUjuv5iL7y5uhKxeAKO7BwU/PuvWW5C\n",
       "lOYklxv+eVdo2qOH9yN29KAJEaUWJt2dCEViePuTfYb0fftlY+Fxd33sFxGlJlHlJOcoLc2IbN2g\n",
       "affNXKA7a9RTn5s7CgGfu8/9XKiuOYTlW44AOFsGk3XLF1luQpTmvJOm667OBdcsNSGa1MKkuxNL\n",
       "NhxAaygqvN+SggAunzJUeL9ElByiyknaC61bqZ3l9njhm31pX8MFAGT43LhlwWghfV3otVW7oSif\n",
       "/RfBchOi9CY5nPAvuFLTHt2/B/FTVSZElDqYdOuIxRN4a81eQ/q+6/Lx3Z79IqLUIrKc5BwlHEJo\n",
       "81pNu2/qJZqTAPrimpnDUZDtF9bfOdV1rVi763iHNpabEKU3z+iJcORqT2dL99luodnfz372M0yb\n",
       "Ng3Z2dno168frr/+euzatUvkI5Ji+ZYjqGsOCe+3oiQPl4wdKLxfIjKW6HKS9sIbV2tPLHE64Zs5\n",
       "v08xX8jtcuALC8cK7fOcV1fuhnrBtD/LTYjSl+Rw6K7URXZtRaK+1oSIUoPQpHvlypX42te+hk8+\n",
       "+QTLli2D0+nEokWL0NDQIPIxhlJVFe8YVMt9zxUTIBlxaC4RGcaIcpJz1FgUoXUrNe3eSTMgZwZ6\n",
       "1efFXDapDKWF4vs9WN2AymP6L1KWmxClJ++E6dpxTFURXLvMnIBSgNCke8mSJbj33nsxevRojB07\n",
       "Fs8++yxqamqwdq126TRVVR6rxeFTjcL7HV6ahwnl/YX3S0TGMaKcpL3wlvVQgm0d2iRZhm+WmFru\n",
       "CzkcMj4/z5ja7vfW7+/8uSw3IUo7ksulu2IX2boRiZYmEyIyn6HFxc3NzVAUBbm5uUY+RqiLvTj6\n",
       "4vPzR3OWm8gijCwnOUdNJBDSmfHxjJ2kWwspyvyJQ1BoQG33mp3H0dga7vRzlpsQpR/vFO3eFDUR\n",
       "113hSweGJt1f//rXMWnSJMyaNcvIxwjT2BrGmp3Hu/5iDw0szMKMUT27wpmIzGFkOUl70b07kWjS\n",
       "rqr5LlnU574vxumQcdPckcL7jcUVfLT5UJffY7kJUfqQvV54p83RtEc+XQc1Jv6EuFQnqRfufhHk\n",
       "sccew8svv4zVq1djyJAh59ubmj5bUti/35hZ5d5auv0k3tl0Qni/d8wtw/SK3t8oR0TJsW5vDX7w\n",
       "4lY0tHV8GcgScP/lw3HfwmFwyGJWrDKWvgXXqY7jTaxkCNoWXCOk/4uJxhP40cvb0RYWO6ucl+nG\n",
       "v35+POTu/HekKPCv+RCZy96CpCgdP/JloOlzX0R0+Dih8RFR8knhELLefAbSBceiBmdehmj5KJOi\n",
       "6r2Kiorz/zk7O7tHP2vITPc3v/lNvPTSS1i2bFmHhDuVKYqKtZVnhPebm+nG5KF5wvslInHiCQVP\n",
       "LNmLf/rTBk3CnR/w4HcPzMBXLq8QlnDLTQ2ahBsAIiPGC+m/K26nA/NGi99jUt8axZ6qbtZqyjKC\n",
       "c69AwxcfQyIrp+NHoTbkPv87ZH7wOpAQfyswESWP6vUhOniYpt29f6cJ0ZhL+LWIX//61/HKK69g\n",
       "+fLlGD784ledT506VfTje21jZRUU5wHk5IitdXzguimYOYNXvlvZpk2bAKTW7yuJU13bgi/8+DXN\n",
       "ZkngbDnJc9+7qU+bJfW0LnkDoZyOiaYjvxDl19wASe7bXEh3f19Hjoli64m3EIqKne0+0uzEvT35\n",
       "/8rUqUgsuhL1v/oBwps6brrPWPMBcutPIf87P4WzcIDQOCk1cHxND7Gifmj84/90bExEUV7UD66S\n",
       "QeYE1UvtKzZ6SuhM9yOPPIKnnnoKzz//PLKzs3Hq1CmcOnUKbW1tXf+wyYzYQJnl9/D2SaIUZvTp\n",
       "JHrUaAThbTpXvk+Z3eeEuycyfW5cNUM7+9RXm/edxMm6lh79DE83IbI3V8kgOHWS6/DG9Pr/tdAR\n",
       "/oknnkBraysWLlyI4uLi8//88pe/FPkY4eqbQ9i876Twfq+YVg6PW/hiAhH1UTJOJ+lMZOcWqOEL\n",
       "LsNxueCZOF34s7pyzczhEH2okqqevWCsp3i6CZG9+aZeommL7NyiOTbVzoS+URRFQSKRgKIoHf75\n",
       "wQ9+IPIxwq3aflRzSkFfSRJw5XTxs0hE1DfJOp1Ej6qqCOnM7HjGTobsE3+MX1f65WZg2ohi4f2u\n",
       "2HpEc0Nld/F0EyJ78oyZpHt8YHirduXPrpK3lpnCVmw9IrzPaSOK0S83Q3i/RNR7ZpSTtBevOor4\n",
       "qSpNu0/nSK1kuXpGRddf6qGT9a3Yd7yu1z/PchMi+5FcLngmzdC0hzet7fUf6VaT9kn3sdNNOFgt\n",
       "/pp6I15kRNQ7ZpaTtBfZvknT5iodAmeReef4T6oowoA88RMEy/s4mcFyEyL70SsxSTTUIn78SPKD\n",
       "MUHaJ90rtx0R3mdRXiYmVRQJ75eIes7McpL21EQckV1bNe3eyTMNf/bFyLJkyCTBx9uPIZ5Quv5i\n",
       "F1huQmQfjrwCuMu0401kx2YTokm+tE66FUU1pLTkqhnDunc5BBEZyuxykvaiB/dqNgxJTifcoyYk\n",
       "5fkXs2jKULidjq6/2APNwQi27BezQZ3lJkT24RmvPR4ysnsr1IT9V63SOumuPFaLMxcsNfeV2+nA\n",
       "Ih4TSGSqVCknaU9vJsc9YixkrzdpMXQm4Pdg3njxZ+XqnX3eWyw3IbIH98jxkFyuDm1KsA3RA5Um\n",
       "RZQ8aZ10GzHLPXN0CQJ+j/B+iah7UqWcpD0lHEa0coemXW/GxyyXTy0X3ue63ScQisSE9slyEyJr\n",
       "k71euEeM1bSnQ4lJ2ibdiYSCNTuPC+832S9zIvpMKpWTtBet3K6ZgZX9GXCXj0h6LJ0ZOagA/QWf\n",
       "uBSJJbBhj/a0lr5iuQmRtXnGTdG0RffuhHLBHQZ2k7ZJd+WxWjQHI0L7zPJ7MHEYryomSrZULCdp\n",
       "L7LjU02bZ8xESI7UuTxLliXMnzBYeL8bKsUn3QDLTYiszF0+ArK/4x/5ajyOaOV2kyJKjrRNuo14\n",
       "EcwbPwhOk17qROkqFctJ2lNamhE9vE/T7hmXOqUl58yfMER4n5v3nRRyiklnWG5CZD2SwwnPmIma\n",
       "druXmKRthrjegCXPBROHCO+TiDqXquUk7UX27cKFfw04cvPhLBU/q9xXg/pno7w4V2ifbeEYdh0+\n",
       "I7TPC7HchMh69CYeYkcOQAmHTIgmOdIy6a6qaUZVbYvQPovyMjF8YL7QPolIX6qXk7QX3bdL0+YZ\n",
       "PRGSlJrHihoxeWDEJMeFWG5CZC3O0sFwZOd0aFMVBdEDe0yKyHjmv5FMYNQsd6q+RInsJNXLSdpT\n",
       "Y1HEDu3VtOvt3E8V88YPhix4LNtQWZW0a55ZbkJkDZIkwT18jKY9ulc7UWEXaZp0nxDepxEbkIio\n",
       "IyuUk7QXPbhP59SSTDhLxJ+JLUpelg/jhvbr+os9cLqhDUdPNwnt82K6VW6ykeUmRGZzD9dOQEQP\n",
       "7LbtRTlpl3Q3t0VQeaxOaJ+lhQGUFGYJ7ZOIPmOlcpL2ovt2atrcw0dDklMrzgvNHF0qvE8jjg68\n",
       "mC7LTX7IchMis7mGDIPk7ni3iRoOI3bssEkRGSu1R34DbNpbDUXwMueMUeJfUER0lpXKSdpTFUW3\n",
       "njuVS0vOmTaiWHifRh0d2BWWmxClLsnphLt8pKZdb+y0g7RLurcdFD+4zhhVIrxPIrJeOUl78epj\n",
       "UNpaO7RJDifcQ4ebFFH39c/LRNmAnK6/2AP7T9SjLRQV2md3sdyEKHXpTURE9+5M2j6QZEqrpFtV\n",
       "Vew4JPboqiy/ByMGFgjtkyjdWbWcpL3oXm1piWvocM1SaqqaLngyQVFV7D5aI7TPnmC5CVFqcleM\n",
       "Ai7YvJ1oqEOi5rRJERkndd9YBjhV34qapmDXX+yBaSOLIcs8tYRIFKuWk1woeni/pk1vp36qmj5S\n",
       "/Aqe6EmP3mC5CVFqkf0ZcA3Ujul6Jz9ZXVol3TsNuKDBiBcTUbqycjlJe0o4jMRJ7SlJ7mHa2sVU\n",
       "NawkD7mZXqF9bj+UGjNXLDchSi16Y2PsyAETIjFWWiXdogd8l1PGpIoBQvskSkd2KCdpL378EFSl\n",
       "49Xnjtx8OHLyTIqo52RZEl5icuhkg2l13RdiuQlR6nANqdC0xY4e1IyjVmeNN5gARtRzjx5cCJ/H\n",
       "JbRPonRjl3KS9qI6MzSuweUmRNI3kyuKhPanqsCuI+bVdethuQmR+ZzFpZBc7g5tSjiExOlqkyIy\n",
       "Rtok3afqW1HXHBLa57gysRdIEKUbu5STXEhvWdQ1ZJgJkfTNWAPGuB0pUmLSHstNiMwlOZxwDRqq\n",
       "abdbiUnaJN3bD4of6MeX9xfeJ1E6sFs5SXtKOIS4Tj23FZPurAyP8KMDdxiwt0YElpsQmcs1RLsa\n",
       "GD3KpNuSRC9pelwOVJTmC+2TKB3YsZykvdixQ7jwX8yRmw9Hdq5JEfWN6NnuVKrr1sNyEyJz6E1M\n",
       "xI9q98dYWdok3furxF79PmpwAZwWnIUjMpNdy0naix09qGmz4iz3OeOGik26VfVs4p3KWG5ClHzO\n",
       "olLNPQZ2q+tOi6wxFImhqrZFaJ/jylhaQtRddi4nuVD82GFNm2uwdZPusWX9Lry3os8OVNWL7dAA\n",
       "LDchSq6zdd0653Uf146pVmX9N1w3HKpu0Cxj95Xo2R8iu7J7OUl7aiKB+OkqTbtr4JDkByNIwO/B\n",
       "kP5i67oPVqf2THd7LDchSh69zZR6e2SsKi2SbtGzKqznJuqedCgnaS9RdwZqLNahTfb6IOdae7wQ\n",
       "Pcmw/4TYcj+jsdyEKDmcA0o1bfHq4yZEYoy0SLpFz6qUF+eynpvoItKpnKQ9vZeDs3ggJNH1GUk2\n",
       "YmCB0P6q61pTejOlHpabEBnPWTxQ05aoPa2ZzLAqe73xOiF6E+WwEuvcKkeUbOlUTnIhvWVQZ5H2\n",
       "JWI15cXiT15J9c2UnWG5CZFx5IxMOLI7lrOpiqJbtmdFtk+6jdhEyaSbSF+6lZNcSH+mW7tcajVF\n",
       "+QH4Bd++a4XNlJ1huQmRcfRLTOxR1237pNuITZTlxUy6idpL13KS9jrbRGmHmW5ZloTPdls56QZY\n",
       "bkJkFL0Sk/hJe9R12/cN+A9HTjUK7c/rdqK0MEton0RWls7lJO11uokyxx5/pItOukWPzWZhuQmR\n",
       "WHoTFXY5wcT2Sbfo0pKhRTmQZWtviiISJd3LSdqLn9SZ5bbBJspzRJ/YdLK+FYoieBnSJCw3IRLH\n",
       "WaQtL0nUnLLFqpHtk+7qOtZzE4nGchKtRN0ZTZtjQIkJkRhD9Ex3LK6gprFNaJ9mYrkJkRhyZgBy\n",
       "oGNFgaooSDTUmhSROLZ/I1bVNgvtb2iR+F38RFbCchJ9ekm3M98+l2gV5QfgczuF9il6JTIVsNyE\n",
       "qO8u/MMVABJ1TLpTWiyewJmGYNdf7IGB/bKF9kdkJSwn6VyirkbTJueJPd/aTLIsoaQwILRP0SuR\n",
       "qYLlJkR948gv1LQl6k6bEIlYtk66T9a1QhF8dElJgdiXDpEVsJzk4lRFQaJem3Q7C/qbEI1xSgrE\n",
       "biI/USN2JTKVsNyEqPf0Z7q1Y6zV2PoNKXoWJSfTiwyfW2ifRKmO5SRdU1qaNSeXSB4vpAx7zfqL\n",
       "nnSotmF5yYVYbkLUc/oz3Uy6U1qV4FmU4nx7vUCJusJyku5J1GqXPR35hbY5ueQc0TPdovfcpCqW\n",
       "mxD1jCNPL+nW7puxGlsn3aJnukW/cIhSFctJekavtETvpWF1xYJnumuagojGEkL7TFUsNyHqPkdu\n",
       "HiS54/tFaWuFEg6ZFJEYtn5jnm4QexwV67kpHbCcpOf0lj31lketTvQYqKqw1bGB3cFyE6KuSQ6n\n",
       "7sViVp/ttnXS3dAi9i+iEt5ESTbHcpLeURq1V5rbMen2eVzIC/iE9tnQGhbanxWw3ISoa3pjqNLU\n",
       "YEIk4tg66a5vETuYFzHZIJtiOUnfKK3aUjY5YM/jRYsE722pb7b2cnFvsdyE6OL0xlC9sdZKbPsG\n",
       "jcYSaA1FhfaZn+0X2h9RKmA5Sd8prU2aNodNk+78LLEz3fWCVySthuUmRPp0k+5m7VhrJbZNukUP\n",
       "5C6njAyvS2ifRGZjOUnfqYqiO/siZdpzD4jo8pJ0neluj+UmRFpyprakV2m19olHtk26Rddz52f5\n",
       "bHf8F6UvlpOIo4aCUBMdT+CQ3B7IHq9JERkrjzPdhmC5CVFHukl3C5PulFQnePYkN1Psi4bILCwn\n",
       "EUtv5kUO2HfTdS5nug3FchOis/TGUb1SPiuxbdIteqZb9OwOkRlYTiKe3syLXTdRAuLLSxpamXRf\n",
       "iOUmRJ0l3dxImZJEz56IftEQJRPLSYyjtGhnXvSWRe1CeHlJc/odGdgdLDehdCdnaPfFKME2S//O\n",
       "2/YN2yL45JLcgD3rM8n+WE5iLKVN57hAm26iBIDcTLFjYTASQyKhCO3TTlhuQulKcjh0x1K9Mdcq\n",
       "bJt0hyIxof3lCH7RECUDy0mMp4a1M7Wyz77Hi/q9LricYl8doah1Z66SgeUmlK4kr3ZlTY1Yd3XM\n",
       "tkl3WPAg7vfwuECyDpaTJI8ai2jaJLfHhEiSQ5Ik4eOh6PHajlhuQulI0jkFSo1ox1yrsO0bNxQR\n",
       "O/D4mHSTRbCcJLn0XgCSx75JNwD4PE6h/QXDYlcm7YzlJpROZJ0JDDXKpDvlhKJiB3HRLxkiI7Cc\n",
       "JPn0XgCS297laD43Z7rNxHITShd6q4YsL0lBnOmmdMJyEvPoznTbuLwEED8JIXoPTjpguQmlBc50\n",
       "W4PwpNvNmW5KTSwnMZf+TLfbhEiSR/QkBDdS9h7LTcjOdGe6mXSnHtHLlZzpplTEchLz6Sfd9p7p\n",
       "9gqehOBMd9+w3ITsSm9/DJPudv73f/8XZWVl8Pl8mDp1KlavTv7/0VVVZU032RrLSVKHbtKts+Pe\n",
       "TkSv/IlemUxHLDchO+JGyot46aWX8I1vfAPf//73sXXrVsyePRtXXXUVjh8/LvIxXUooqmaZvS8c\n",
       "sgSX09H1F4mSgOUkqUWN6/yB77L3ypjome44L8cRhuUmZCsubameGrPuypjQpPtXv/oV7rvvPnz5\n",
       "y1/GiBEj8Pjjj6OoqAhPPPGEyMd0SRWZcePsUj1RKmA5SQrSGW8kyd4rDKJXUESP2emO5SZkF5Jk\n",
       "r/xL2MgZjUbx6aefYvHixR3aFy9ejLVr14p6TLeIHr9lm/2PTtb0739ZznKSVKQ34Nh8zBD9b6cw\n",
       "6RauW+Umf/kN/+Ch1KY3lirWXRkTtkZYW1uLRCKB/v37d2jv168fTp3SX8ratGmTqMd3EIkl0NjY\n",
       "KKw/j0s2LFayDrN/B05UVWvyu+kVBfjRHRMRSNRh06Y6cwJLc9m1tZAuqDE8tnUrVJPruo38fT16\n",
       "7LjQMXb/gQPY5GsT1h91JH3528h+4y/w7N/Vof1UdTUObN5sUlQdmT2+UmpyHzwI/wVjzZkjRxAy\n",
       "8feloqKi1z9ry2kx/t1OdnT/5RWYPDQPACBLwIOLh+Pxr0xHfsDeJ2VQ6hE+j89B21BqRiYav/AI\n",
       "WhbdBFU++9qPlpahddGNJkdG1AW9mW4Lr84Im+kuKCiAw+HA6dOnO7SfPn0aRUVFuj8zdepUUY/v\n",
       "IByNI+ftQ8L687gchsVKqe/cDEwq/A68XT4KV33nefz6a1dws2SKqFv2BpRQx5KfoZMmQfZnmBJP\n",
       "Mn5fd9Q4kHNc3K1ww4YNw9Spo4T1R52YPh2Rxdeg4Ymfo+gH/4PyfgPMjiilxldKPSElhNb92zu0\n",
       "+cqGINPE35empqZe/6ywmW63240pU6bggw8+6ND+4YcfYvZs7S5qI4muwbbwH1VkM8UFAWz9I08n\n",
       "SSk2m4npDkUR++9n8xL4lHLudBNnCiTcRF3Srd+27oAh9Nynxx57DHfffTemT5+O2bNn4/e//z1O\n",
       "nTqFr371qyIf0yXRA7jKtU9KIXbbzW15aZh0i/6342b15OIYQpahN5bK1q2MFpp033rrrairq8OP\n",
       "f/xjnDx5EuPGjcN7772HgQMHinxMl0QPKImEClVVOVARkZbOuKCq1t1d3x0Jwedqc2wlIj12O11H\n",
       "+DWLDz30EB566CHR3faIQ5YgS5KwY6gUVUUsrsDt4gU5RNSRpHN5Ayx8eUN3hKNibzV0Oa07c0VE\n",
       "BopFNU26Y65F2HKkkyRJ+I1poYi9X6KU+u6//37IsozHHnvM7FCoHUnvmuKIuE2GqSgkOOn2eex9\n",
       "g2eqeOqppyDL8vl/nE4nSktLcdttt2Hfvn1mh0ekoehc+S55rHtily2TbgDweQQn3YJfMkQ9EQqF\n",
       "8PLLL8Pn8+GFF15AIpEwOyT6B8mtnXVR7J50C56E8AmeJKGLe/XVV7Fu3Tp8/PHH+NnPfoYtW7Zg\n",
       "4cKFaG5uNjs0og7UiE7S7WLSnXKEJ92c6SYTvfnmm2hpacHPf/5znDlzBkuWLDE7JPoHya1zCY7O\n",
       "kqidhCKc6bayiRMnYvr06Zg1axbuvvtuPPHEE6iqqsInn3xidmhEHeitGnKmOwX53GIHcdE1jEQ9\n",
       "8fTTT2PUqFF45JFHUFxcjKefftrskOgf9MtLtLMzdhKKCp7pFjxJQj0TCAQAADGb70Ug61FjOjPd\n",
       "OmOuVdg26RZd0x0MczAic1RXV2Pp0qW47bbbIEkSbr31Vrz99ttCr+Gm3tObddF7UdiJ6EkI0eM1\n",
       "Xbpk6+kAACAASURBVFw8Hkc8HkckEsGePXvwve99D/3798eCBQvMDo2oA93yEibdqYc13WQXzz33\n",
       "HBKJBG6//XYAwO23345IJIKXXnrJ5MgI6GSmO2zvmu5gmOUlVjZy5Ei43W74fD6MGTMGlZWVeOed\n",
       "d5CZmWl2aEQdqNxIaQ2iB3HOdJNZnn76aUyYMAHDhw8HAEyfPh1lZWUsMUkRehsp9V4UdqGqqvjy\n",
       "Es50J9Wbb76JTZs2YePGjXjzzTcxevRoXHXVVaisrDQ7NKIOdJNuvX00FmHfpFvwIN7Yau+ZK0pN\n",
       "mzZtwp49e3DttdeisbHx/D/XXXcd1q1bh/3795sdYtqTvD5NmxIOmRBJcrSFY4jFxV6Ow/KS5Bo7\n",
       "diwmT56MKVOm4Prrr8ff/vY3qKqKH/7wh2aHRtSBGgxq2iQPk+6UE/CLXX6ob7HvS5RS17nZ7J/8\n",
       "5CfIy8s7/8/jjz8OAHjmmWfMDI8AyBkBTZvS0mRCJMlR3yx2LMzwuuBw2PZVZAlerxdlZWXYsWOH\n",
       "2aEQnacm4lCCrZp22cJlULYd6fKytLNPfSH6RUPUlWg0ihdffBEzZ87EihUrOvyzfPlyTJw4Ec8+\n",
       "+6zZYaY9OZCtaVNa7HvesegJiLyA2LGaei4YDOLgwYMoLCw0OxSi85RWnYTbnwnJYd2VMetG3gXR\n",
       "A3lDK5NuSq53330X9fX1eOihhzBv3jzN5w8++CAeeughrFixgqcOmEgOZGnalFb7Jt0NopNuwRMk\n",
       "1LUtW7bgzJkzUFUVJ0+exG9/+1s0Njbi0UcfNTs0ovOUVu2Kod54ayX2Tbo5000W98wzzyArKwu3\n",
       "3HKL7ud33HEHHnvsMTzzzDNMuk0kZ+ok3S1NUFUVkiSZEJGxRI+FnOlOnnO/j+3HlMLCQowbNw5L\n",
       "lizB5ZdfblZoRBp6K4Z6462V2Dbpzs0UW2jf0Bq27UuUUtMbb7xx0c+zsrLQ1taWpGioM5LPD8nh\n",
       "hJr47Bg9NRaDGolA8lp3w09nRJeX5Abs999Rqrr33ntx7733mh0GUbfo7Y2x+kw3a7q7KRZX0Bqy\n",
       "99XORNRzkiRBztTZTKmzNGoHDS1iT3LiTDcR6dGd6WbSnZpcTgcCPu35uX1R26Q9uoaIKJ02U9Y1\n",
       "ix0HWdNNRHr09sbImdqx1kpsm3QD4gfzk3XanbRERLqbKZvtOdNdLXgc5Ew3EenRG0P1VhWtxN5J\n",
       "t+DBvKrWnjNXRNQ3cnaupi1RX2NCJMZqC0WFXxTGmW4i0qM3hso5eSZEIo6tk+7+uRlC+6uqbRHa\n",
       "HxHZgyO/n6YtUXfGhEiMJXoMlCUJBdl+oX0SkfWp8TgSjfWadke+tc+St3XSXVwgdhmCM91EpMdR\n",
       "oJN019eaEImxRI+B/XL9cDkdQvskIutLNNQCqtqhTQ5kQ7bwFfCAzZPukgKxu1yra1nTTURaerMv\n",
       "idqzl4/YSbXgmW7RYzQR2UOiTjtpYfVZbsD2SbfYme7mYAQtwYjQPonI+uTMLEiujqclqbGo7jmz\n",
       "Via6vKQ439qboojIGIm605o2Jt0prn9eJhyy2MtsWNdNRBeSJEm/xKTOXpspRZeXlBZyppuItPTG\n",
       "Tkcek+6U5nTIwjdTHj9jr5krIhJDt8TERkm3oqjiZ7oFr0YSkT3oJt06ExtWY+ukGxBfM3iwukFo\n",
       "f0RkD3Y/weRETTMisYTQPkWXABKRPeiNnZzptgDRg/qBKu0RNkREerMw8VNVJkRiDNFjn9vpQH4W\n",
       "jwskoo4SLU1Q2joeXCE5HHDkWvuMbiANkm7Ry5eHTzYikVCE9klE1uccUKJpi588AVWxx3ghOuku\n",
       "LsiELHjPDRFZX7z6uKbN0a8IksNpQjRi2T7pLivS3hTXF9F4AsdreF43EXXkyCuE5PZ0aFMjYSgN\n",
       "dSZFJNbBarFJd9kAsWMzEdlD/OQJTZuzqNSESMSzf9I9IAeyJHY2hSUmRHQhSZZ1Xwyxk9pZG6tR\n",
       "FFX4fpZhJdZfKiYi8fST7oEmRCKe7ZNuj9uJgf3EbqZk0k1EevSSbr2lUqsxYhMlk24i0hPXmajg\n",
       "TLeFVAge3Jl0E5Ee3aRbZ9bGakSPeZIElBXlCO2TiKwv0dIEpaVjCa/kcMDZv9ikiMRKi6S7XHDS\n",
       "fehkA6KCZ32IyPqcxdolUDtsptx3QmxdemlBFnwel9A+icj69CYpHP2KIDmtv4kSSJOkW/QyZiyu\n",
       "CH8JEZH1dbaZMlFfa1JEYuw4pL2SuS8qSllaQkRaeuV4eidDWVVaJN1GbKbcflDsS4iIrK+zzZTx\n",
       "44dNiEaMxtYwjp0Re2JTeTGTbiLSih87pGnTW0G0qrRIuo3YTLnzsH1umiMicVwDh2jaYkcOJD8Q\n",
       "QYwY67iJkogupMbjiJ04oml3DRqa/GAMkhZJNyB+M2Xl8VrWdRORhmvIME1b7OgBqKpqQjR9J7q0\n",
       "xCFL3ERJRBrx6mNQY7EObbI/A46C/iZFJF7aJN1jy7RXNPdFLK6g8pi16zSJSDxn6RBIcsehNdHU\n",
       "aNlLcrYLTrrLi3O5iZKINGKH92vaXIOHacZTK7PPv0kXRCfdAEtMiEhL9njhLBmkaY8dtV6JSUNL\n",
       "CCdqWoT2OW6ofWatiEic2NGDmja9lUMrS5uku39eJvrnZgjtU/QMEBHZg26JiQXrunccEj+xMM6A\n",
       "CRAisjY1HkdMZ8M5k24LEz3YVx6rRVsoKrRPIrI+12DtiyJ65KDl6ro/3X9SaH8OWcLoIYVC+yQi\n",
       "64udOAI1Hu/QJvsz4Si018pYeiXdgpc1E4qKzfvEvpSIyPpcA4dAcjg6tCnN1qrrVhQVGyurhfY5\n",
       "rCSP9dxEpKG3EugaUg5J8HHPZkuzpFv8suaGyirhfRKRtUluD5zF2rru6IE9JkTTO5XHatEcjAjt\n",
       "k6UlRKQndqBS02a30hIgzZLuwpwMDMgTW9e9aW814glrX/FMROK5hg7XtEX37jQhkt5Zv0d7HXNf\n",
       "cRMlEV1IaWlGrOqopt09dIQJ0RgrrZJuABhXJnbQbwvHsPtIjdA+icj6PCPGatpiRw9CCYdMiKbn\n",
       "RK/iOWQJowYXCO2TiKwvemC3ps1R0A+OfPvt/0i7pHvisAHC+2SJCRFdyDGgBHIgu0ObmkggdnCv\n",
       "SRF1X1VNs/CjAkcOKmA9NxFpRHRWAD3DtZMWdpB2SfeU4UVwyGIL89fvOWG5UwmIyFiSJOnOduu9\n",
       "YFKNERMJM0aVCO+TiKxNjUURO7RP0+7WGTvtIO2S7gyfW/hFOafq23D0dJPQPonI+tzDx2jaYvt3\n",
       "Q03Edb6dOtbtFl/PPX0kk24i6ih6aJ/u1e/O0sEmRWSstEu6AWMG/5Vbjwjvk4iszVVWAcnt6dCm\n",
       "hEOIHT9iTkDdcKahDbuP1grts6QggJLCLKF9EpH16W0ud1eMttXV7+3Z89+qC9MNWOZcue0oFIUl\n",
       "JkT0GcnphLt8pKY9lU8xWbntiPA+WVpCRBdSFQXR/dpNlHorhHaRlkn3gLxMDO6f3fUXe6CmKYjd\n",
       "R3mKCRF15B6hfYFEdm2FqqTeUaOqqmL5liPC+2VpCRFdKHbsEJTWjhu2JYcTrnL7HRV4Tlom3YAx\n",
       "My/LtxwW3icRWZveUqnS0oTY0YMmRdS5wycbcbymWWifWX4PRg7iUYFE1FFk+yZNm2vocMgerwnR\n",
       "JEfaJt1GzLys2Xkc0VhCeL9EZF2yPwOuYaM07ZEdm02I5uJWGLA3ZeqIIjgcafuqISIdaiyGyO5t\n",
       "mnbP+CkmRJM8aTsSVpTmIzdT7F9TbeEYNu2tFtonEVmfd5z2RRLZvU2za99MiqJi5TbtrXB9NWNU\n",
       "qfA+icjaovt3Q42EO7RJbo/uMat2krZJtyxLmDNukPB+jZgpIiJrc48YoznFRI2EEd2/y6SItLYf\n",
       "Oo36FrG3Zfo9LkwdUSy0TyKyPr3SEs/IcZBcbhOiSZ60TboB4NJJQ4T3uaGyCg2CX1xEZG2Syw3P\n",
       "qPGa9sj21Ckx+WCj+Brz2WNK4XY5hPdLRNalBNsQPVCpabd7aQmQ5kn3sJI8lBQEhPaZUFRDXl5E\n",
       "ZG0enRKT6IE9UIJtJkTTUUNLCGt3HRfe74KJQ4T3SUTWFtm9TXNBmJwZgGtIhUkRJU9aJ92SJBny\n",
       "Uliy8SASidQ7DoyIzOMqq4Ac6HhBjJpIILJri0kRfeaDjQeREHzPQF7Ah3FD+wvtk4isT7e0ZOxk\n",
       "SA77r4qlddINAPMniL9qtLYpiI3cUElE7UiyDM+YSZr28OZPoKrmXayVSChYYsDq3PwJgyHLkvB+\n",
       "ici64mdOIXZce7yy3kqgHQlLuhsaGvDoo49i1KhR8Pv9GDRoEB5++GHU19eLeoQhivIDGDkoX3i/\n",
       "763bL7xPIrI274Rpmrb46WrETbwWfuPeatQ2BYX3y9ISIrpQeNNqTZuzcACcRelxypGwpLu6uhrV\n",
       "1dX4xS9+gZ07d+K5557DqlWrcMcdd4h6hGEWTBgivM8tB06hSvAlE0Rkbc4BJXANHKJpD21ak/xg\n",
       "/sGICYJB/bJQVpQjvF8isi4lEkZYZ/P4/2/vzuOiuNK9gf+qem+aRXZBFIkgCOICgoILrhGN2dQk\n",
       "OtFMMpkkN3NnsszMnTf3de4kk0w+b2ZJbpKZTDSTRbO7TDb3fUNUwA0QFxSFqIAINPS+VL1/GFFS\n",
       "zSZdVd3wfD8fP0lOFec8w5TdT58+5znazFwwTP/4VsxrSXdqairWrVuHu+66CwkJCZg8eTL+8pe/\n",
       "YPv27TCZTN4aRhQTRw6GQoSvQTcfrvR6n4QQ/6bNzBW0OU4eFxyHLIVLV1twtLLW6/3mjY7vN2+i\n",
       "hJDusZeWCGtzq9TQpGfKFJH0RF3TbTQaodFooNfrxRym14INWlFOqNxWch4Wm+8cfkEIkZ9mxGiw\n",
       "+oB2bbzbBduxQ5LHsr7wjNf7ZBkGU8cM9Xq/hBD/xfM8bEXCb/Q06ZlgtX332PcfEy3pbm5uxu9/\n",
       "/3s88cQTYFnf3685Z7z3S9WYbU5sOkRruwkhNzFKJbRjxgvabcUHwHPSVT1qNtmwtfi81/vNTolF\n",
       "eLBvT7QQQqTlqqmCq/6KoF03TvjNX1/G8F1sm1+2bBleffXVTjvZvXs3Jk+e3PbfJpMJ+fn5UKlU\n",
       "2Lx5M9TqmycMGY3Gtn8/e9Z3ElKO4/HaV2WoN9q6vrkHgnQqLFuYDpXS9z94EEKkwZpaEPTtJ8CP\n",
       "Xn5NeXPhio2XJIaNJd9j23Hhm2BvPXVnEobHBnu9X0KI/9IXbIX6QvuczxUxEKZZ98sU0e1LTLw5\n",
       "SRsc3LPXOmVXNzz33HNYunRpp/fExcW1/bvJZMKcOXPAsizWr1/fLuH2ZSzLIGd4BL4+7N0DIlqs\n",
       "ThyubEBucqRX+yWE+C/OEARnzBCoLl1o1645dUKSpNvmcGN/Rb3X+40M1iJxYFDXNxJC+g3GYoKq\n",
       "WliW1J40UoZo5NVl0h0WFoawsO6V1GttbUV+fj4YhsGmTZu6XMudmelbi+eTUx0orPoadqfbq/2W\n",
       "1brwn4vGQqGg2W5/VFx8vZC/rz2vxL85ggNg/GxF+0ZbK4YNjIQqdvBt99ud53XdnpPQ6AOh8fIq\n",
       "kKVzxiArK9m7nZI+jV5f+z7Tlq9hDWr/YZwNMOCO+xaCUXSZhvqcW1ds9JTXssDW1lbMmjULzc3N\n",
       "+PDDD9Ha2ora2lrU1tbC6fSPzYQGnVqUw3JqG83YX1rt9X4JIf5LdcdwKMIiBO3Wgh2ijutwuvFN\n",
       "wWmv96tRKTAjI8Hr/RJC/BdnMcNWUiho12bm+mXC3VteS7pLSkpw6NAhVFRUICkpCTExMYiJiUFs\n",
       "bCwKC4W/cF81J9v7GyoBYO3ek7KeOkcI8S0My0KfO03Qbq84AVdDnWjjbi85jyaTd/euANdPoAzQ\n",
       "+cdyQkKINKyH94F3Otq1MSoVdFmTZIpIXl5LuvPy8sBxHNxuNziOa/vjdrvbbbL0dXfEhopyQuWF\n",
       "WiMOVVzyer+EEP+lSc8EGyjciGMt2CnKeE6XG+v2nhSlb7EmLAgh/omz22A7tFfQrs3IEZRN7S9o\n",
       "kbEHd41PEqXfj7ceB8fRbDch5DpGoYR+Qp6g3X6iGG5jk9fH21J0DvXN3j/yfcSQcNwRG+r1fgkh\n",
       "/st25CA4m7VdG6NQQDc+T56AfAAl3R5MHDkY0aHe/xRWXd+CXUervN4vIcR/acaOB6trv6OR5zhY\n",
       "D+zy6jhWuxNf7Czzap83LJgyQpR+CSH+iXe5YC3cLWjXjMyAIjhE+oB8BCXdHigULO6bmCJK359u\n",
       "L4XDy9VRCCH+i9VoofWwvtF29CA4s8lr43xTcBpGs91r/d0QHx2MzOExXu+XEOK/7CeKwbUKq3zo\n",
       "cqfLEI3voKS7AzMyEjDA4P2jSa8aLXRKJSGkHV3WJDCq9psQeacTlv3bvdK/0WTDv/dWeKWvH1sw\n",
       "eQQYhhGlb0KI/+GdTlj2bhG0a1JGQRnev88soaS7A2qVAvfkDhel79W7T8Ji848yioQQ8bH6AGgz\n",
       "JgjabUUFcDc39rr/NXtOwupw9bqfH4sODcDEkbdfU5wQ0vfYigvgNjYL2nUeqjX1N5R0dyI/OxEB\n",
       "WpXX+22x2PHVPnFmnQgh/kmfMw2Mqv3rDe92wbJ7c6/6rW8yY6NI367dPymFDv0ihLThbFaP39Cp\n",
       "k1J7dehXX0Gvlp3Qa1WYO16cMlhfF5zGNaP3qwgQQvwTGxjkcVe/7XgRXPVXbrvfj7cdh9PF9SIy\n",
       "zwYYtJg+lg7DIYTcZD2wC5zFLGgPmD5Xhmh8DyXdXZiXMxxqpcLr/docLny4+ZjX+yWE+C9dzlSP\n",
       "9WvNOzfcVn9lVfXYfexib8Py6J7c4VCrvP/aSAjxT1xrC6wHdwvatemZUEYOlD4gH0RJdxdCDFrM\n",
       "zrpDlL73HL+IE+fEO3mOEOJfWK0OuokzBO2O0+VwVves3KjbzWH5d8XeCq2dIL0G+XQYDiHkFpZ9\n",
       "W8E72+9XYxRK6KfmyxSR76GkuxsW5qVCp1aK0vfy74rhcnv/q19CiH/SZeZ6rGNr3v4deL77h2tt\n",
       "OHgWF2qFJbu8YWHeCOhF2O9CCPFP7sYG2EoKBe3acblQhNDBWTdQ0t0NIQYt7puULErf1fUtWF94\n",
       "RpS+CSH+h1GpoJ8yW9DurKmCo+JEt/posTjx6fZSb4cGAIgI1tOR74SQdszbvgXPtZ9AZDRa6D18\n",
       "c9efUdLdTfdOTEZwgEaUvj/fUYbGFmvXNxJC+gXNqHFQRkQL2s1bvwbv6PqAm/XFNbDYxSlL+pMZ\n",
       "I2ktNyGkjaPyFOynhB/y9RPywAYYZIjId1HS3U06jQoPTk0VpW+L3YkPNx8VpW9CiP9hWBZ6D7v9\n",
       "3cbmLg/MOV/XiqLKa6LENTgyCFPHDBWlb0KI/+FdLpg2rRO0s4ZAaMdPkSEi30ZJdw/MzhqGqAHC\n",
       "ygLesPvYRRyrrBWlb0KI/1EnpUKdmCJotx7YDfe1qx5/xulyY+0BcaqVAMDSWaPAsnT6JCHkOmvh\n",
       "brgbGwTtATPuBqvx/qne/o6S7h5QKRX4yYyRovX/1rpDdFIlIQQAwDAMDLPvB6Nov4mbd7tg2vxv\n",
       "j5sqv9xVjitN4ixVSx4chqyUWFH6JoT4H7exCZZ9WwXtqsEJ0KRnyBCR76Oku4emjIpHfHSwKH1f\n",
       "NVrw4SZaZkIIuU4RGg5dzlRBu6PyFByny9q1nf3+GtbuOSlaLI/cORoMQ7PchJDrzFu+EZYIZFkY\n",
       "5syn14oOUNLdQyzL4LH8MaL1v7noHC0zIYS00U+a4bmE4OavwDsdAK4vK3lz3SG4ue6XFOyJ8SNi\n",
       "kTY0UpS+CSH+x3HuNOwVxwXt2sxcKKNiZIjIP1DSfRvGJA5ETuog0fqnZSaEkBsYlRoBs+4VtLuN\n",
       "TbDs3gLg+rKSi3Xi1ORWKxV4fM5YUfomhPgf3umAaeNaQTsbYKCDcLpASfdtenzuWGhEKptFy0wI\n",
       "IbdSp6RDnTBc0G4p3IXK4mOiLitZmDcCUaFU9osQcp1558YONk/OA6vVyRCR/6Ck+zZFhATggTxx\n",
       "SggCtMyEEHITwzAIyPewqZLjcPK95YDLJcq40aEBuH+SsIIKIaR/clZXwXpwj6BdFRcPTXqmDBH5\n",
       "F0q6e+G+ScmICRNvBuiNNQfRbLKJ1j8hxH8owyOhnzKrXVtNvRHKlmsY11whyphP3JVBB+EQQgBc\n",
       "X1bS+s1ngnZGoYRh3oNgWEopu0K/oV5QKRV44i7xyuI0tlrxxppCcCJtjiKE+BddzjQoB17fT9LU\n",
       "akVdkxkAMMpYiUhbo1fHykqOwbhkKhFICLmuo2Ul+qmzPZ6gS4Qo6e6ljOExGD9CvDemI2dr8dU+\n",
       "cWaxCCH+hVEoEHjvYjg4oOpKc1s7Cx7TGo5Awbm9Mo5KyeLnIk4oEEL8S4fLSmKHQDdBWNaUeEZJ\n",
       "txc8Pke8TZUA8PG2EzhVLfx0SQjph8KisIGPhZvj2jUPcLZ6bZnJgskjEE2bJwkh6GJZyb2LaFlJ\n",
       "D9BvyguiQg1YMjNdtP7dHI8/f14Ak9Uh2hiEEP/w6bYTWO8ciHq1sHb3KGMlBjl6t8xkSFQwFuaN\n",
       "6FUfhJC+w7zt2w6WleRDGR4lQ0T+i5JuL5mXMxwjhoSL1v9VowVvrjvo8ehnQkj/cOTMFazdWwGe\n",
       "YbErIgNupv1LOAsec1vKoONu7wO6gmXwzPxsqJS0eZIQAthPHoe1qEDQrhoUD92EPOkD8nOUdHsJ\n",
       "yzL41f3ZUIv4ZnXw5CV8vf+UaP0TQnxXfZMZf1td2PbfjeogFIckC+4L5GzIbykDcxsf0O+flILE\n",
       "QWG9ipMQ0je4m66h9bsvBO2MQgnDPQ/RspLbQL8xL4qNCMLSWeItMwGAjzYfR8npy6KOQQjxLTaH\n",
       "C698vBctFnu79qPBibisESbJCY4GpLdU9miMwZFBWDQ9rVdxEkL6Bt7tQsu6VeBtwrLFAbPuoWUl\n",
       "t4mSbi8Te5kJx/P4y5cHcOlqi2hjEEJ8B8fx+N+1B1FV2yy4xjMstkdmwsqqBdfGN5Z3u4yggmXw\n",
       "7ILxtKyEEAIAMO/YANelakG7JiUd2nG5MkTUN1DS7WVSLDMx25x4+eO9MNPGSkL6vNW7y1FQVtPh\n",
       "dbNSj10RwvJ+LHjMvFoEtdvZ5Ri0rIQQcoP9TDmshbsF7YrgATDc/RAYhpE+qD6Ckm4RSLHM5FJD\n",
       "K/78RQEdnENIH3agrAafbi/t8r6L+mgcDxomaA9yWZDXcAToZH03LSshhNzgNjbD9M3ngnaGZRG4\n",
       "YClYrU6GqPoOSrpFMi9nOEYOjRR1jCNna/HR5mOijkEIkceF2ma8vqaw6xt/cDA0FXWaAYL2OyyX\n",
       "McZ41uPPKBUsnls4gZaVEELAO51oWf0BOItZcE0/bS5Ug+KlD6qPoaRbJCzL4NcPTECQXiPqOF/t\n",
       "P4VtxedEHYMQIq3GFiteXrUHdmf3T5jkGBbbI8bBzigF17KbyjHEckXQ/ujs0RgWG9qrWAkh/o/n\n",
       "ebR+9yVcl4VL2dSJKVQe0Eso6RZRWLAezy8cL/o4//i6CIcrLok+DiFEfGarAy9+tBv1zZYe/2yL\n",
       "KgBbA4UH2zAAZtQXY4Dj5gbs7JRYzMtJ6k2ohJA+wnpgJ+ylJYJ2NjAYgfcspvKAXkK/RZFlDI/B\n",
       "/ZOEtXS9yc3xeO3zAlRcvCrqOIQQcTmcbvzpk30eK5V012ltNIr08YJ2Ne9Cft1BaNwORATr8cz8\n",
       "bNoQRQiB/Uw5zNvXC9oZhRJBD/wUbIBBhqj6Jkq6JbBk1igkDxa3MoDD5cYfV+1FdZ1R1HEIIeLg\n",
       "OB6vrylEaVV9r/vaG5CIizphHd1glxl3Xi3Cbx/IRqDIS98IIb7PdbUWrf/+xOM1w7wHaB23l1HS\n",
       "LQGlgsVvHsiBQSespetNJqsDf/hoNxqMPf9amhAiH57nsWJ9SaelAXvUH8Nge8Q4NKkCBdeyDTbE\n",
       "nTrglXEIIf6Ls1rQ8sX74O3CA3B0E/KgHTVOhqj6Nkq6JRIVasAv78sSfZwGowX/88EutP7o5DpC\n",
       "iO/6clc5Nhz0XGHkdjkUKmyKyoadVbW1Bek1GBhqgPXwPliL9nt1PEKI/+BdLrSu+QjuxgbBNfWw\n",
       "ZATMmCdDVH0fJd0SykmLw7wJ4m9cqrnagj+u2gOrvetDMQgh8tp48Gy3anHfDqMqENsixoEDA5VS\n",
       "gYSYAcAP67hNG9fBfvK4KOMSQnwXz3Fo/fpTOKqEH/QVYREInL+UNk6KhH6rEntszhiMukO41tLb\n",
       "TlVfw4sf7abEmxAftvlwJf75bbGoY9Too1AUkY7E2FBBPe7Wf38Cx4VKUccnhPgOnudh3vI17OXC\n",
       "Mz4YrRZBDz1OB+CIiJJuiSkVLH63KBfRoQGij3XyYgNeXrUXNodL9LEIIT2zrfgc/vF1kSRjTXt0\n",
       "McJyJgnaebcLLV+8D1ctlRwlpD+w7t8B6+F9gnaGZRF0/1Iow8U91K+/o6RbBoF6DZY9PBk6tfAQ\n",
       "C28rrarHy6v2UOJNiA/ZXnIeb391WJKxFkxOwZQxQ2GYuwDqRGENb95ug/HT5XA3XZMkHkKIPGxH\n",
       "D8K8c4PHa4Z5D0KdmCJxRP0PJd0yGRIdgl8/MAFSlMk9cb4eL62kpSaE+IIthyvx5rpD4Hnxxxo3\n",
       "PAZLZo0C8EPN3YWPeCwBxplaYfzkXXBmk/hBEUIkZz9VhtbvVnu8FjDjLmhHi1/ogVDSLavsEYPw\n",
       "8Ix0ScYqq7qKP3y4GxYbJd6EyGXjwbP4u0RLSuIigvCbB3PAsjc/2TMqNYIWPQ5FuHBfibuxAcbP\n",
       "VoCzCcuHEUL8l/PiObSuWwlPn/R146dAlzNNhqj6J0q6ZbYwbwQmjRwsyVgV1Q34v+/vQLOJ3lQJ\n",
       "kRLP81izu1z0TZM36DVKLFsyGXqtSnCN1Qcg+OEnwQYGC665Lteg5bPl4DzU7SWE+B9ndRWMn70H\n",
       "3iVcYqpJG4uAmXfTybQSoqRbZgzD4Jn52UgaFCrJeJWXmvC75dtQ20hfIxMiBY7j8d76I1i17879\n",
       "RgAAHLpJREFU9YQk4ylZBkvzEhATLjwY5wZF8AAEP/yUxyoFzpoLaPl0BSXehPg55/cXYPxsBXiH\n",
       "8NwOdcJwBN67iEoDSox+2z5Ao1bif5ZOQUyYQZLxLl8z4b/e3Ybzl5skGY+Q/srl5vC31QfwXeEZ\n",
       "ycZcNHkohscKZ7F/TBkZjaBFPwejFG7odtZUXT+pzsObNSHE9zkvVcP46XKPp00qY+IQ+MBPwSjE\n",
       "L+ZA2qOk20cEG7T446NTERooTX3MJpMNL7y3A6Xn6yQZj5D+xmp34qWVu7H3RLVkYz4+ZwzGJoR1\n",
       "+37V4KEIeuhnHt98nRcqYaQZb0L8jvP7CzB+/E/wHvZnKKNiEPyTJ8FqtDJERijp9iFRoQa89Gge\n",
       "9BrhOkwxWOxO/OGj3SgolS4pIKQ/aDbZ8N//2oFjldJ9qF0wOQX3TEzu8c+p70hG0IOPek68q8+j\n",
       "5dPltLmSED/hrK6C8ZMOZrgjoxG85D/A6sU/J4R4Rkm3j4mPDsHvl06GSinN/zVOF4fXvijAtwWn\n",
       "wUtRw4yQPu7S1Rb8bvk2VF6SbvnW9LFDsfTOUbf98+rEEQhc+AgYhUJwzVlzAS2fvAvOYu5NiIQQ\n",
       "kTmqzsL4ybseE25FeBSClzwNNkCaZazEM0q6fVDa0Ej85oEcsBLtKOZ54L0NR/D3rw7D6XJLMiYh\n",
       "fVHJ6cv49T+34vI16TYqjxseg/+8L6vXFQg0w9MQuOCnnhPvSxfR/OHbcBubezUGIUQc9vJjaPl0\n",
       "BXinQ3BNGRGNkKVPgzV0vLmaSIOSbh+VkxaH/7gnU9Ixtxafx7L3d1JJQUJ6iOd5fLWvAn9ctRdm\n",
       "CWvhJw8Ow+8W5UKp8M5LuSY5DYEPeF5q4m6oQ/MHb8JVX+uVsQgh3mE9vB8ta1eCdwvLAiojoxG8\n",
       "9GmwgUEyREZ+jJJuHzY7axh+2ouvjG/HyYsNeP4fW6iyCSHd5HC68b9rD+KDTcfASbhEK2FgCP5n\n",
       "6RRo1N6tQKBJSkXQg495TLy5lmYYP3wLzuoqr45JCOk5nudh3rkRpk3rPF5XRsUg+JFf0Ay3D6Gk\n",
       "28fNnzICS2ZKc2rlDVeNFvzX8m3YTxssCelUY4sV//2vHdh59IKk4w6NDsErP5uGQL1GlP7ViSkI\n",
       "evhJMB4qHHA2K4wfvwP76TJRxiaEdI13u2FavxqWfds8XlfFxV+f4dbTGm5fQkm3H3hgaip+MmOk\n",
       "pGPanW689nkBVm05Drebk3RsQvzBqeoGPP/OFpyuuSbpuEOigvHyY1NFS7hvUMcPQ0gHs2S8y4XW\n",
       "1R/CdvSgqDEQQoR4pxOtaz6C7Yjnv3/qpFSqUuKjKOn2Ew9NS8OiaWmSj7tmz0n83/d3osFokXxs\n",
       "QnwRx/FYt+ck/s+K7bjWYpV07MGRQXjlZ9MQbJCmxq5y4CCEPPYMFKHhgms8x6H12y9h3rkBPEcf\n",
       "zAmRgrvVCOOqjr9p0o7JRtADj4JRqSWOjHQHJd1+ZPGMkZLPeANA+YWreObtzSg6dUnysQnxJc0m\n",
       "G15auRsfbTkONydtic346GD86fHpCJEo4b5BMSAMIY/+CsqYOI/XLfu2o2X1B3SIDiEic16qRvN7\n",
       "b8D5/QWP1/UTp8Mw70GPFYiIb6Ck2888NC1N8jXeANBiseOPq/bi/Q1H4KLlJqQfOnGuDr96axOO\n",
       "nJW+esfQ6BD86WfSJ9w3sIZABC99GuqE4R6vO06Xw/jBm3A3NkgcGSH9g63sCIwfvQ2u1ejxumH2\n",
       "fQiYflevS4cScXk96eZ5Hvn5+WBZFuvWed5RS3rngampeCx/tCxjf11wGr9bvg21jdLVISZETm43\n",
       "h0+3ncCyD3aiSYZymomxofjT49MQFCDuGu6usBotghY/Dk3aWI/XXfW1aP7XG3BUnZU4MkL6Lp7j\n",
       "YN6+Hq3rPgbvEpYEZBQKBM5fAl32ZBmiIz3l9aT7b3/7GxQ/fLVBn7jEc9+kFDw7PxsKVvrf8Znv\n",
       "G/HM25uxveQ8nWJJ+rQr11rx3//agS92lUOOR33MsGj86XHxqpT0FKNQIvC+nyAgb7bH65zVgpZP\n",
       "3oX18H56bSCklzibDS1ffgBLwQ6P19kAA4KX/gLaDj4IE9/j1QKvRUVFeOutt1BSUoKoqChvdk08\n",
       "mJ6RgKAADV77vAB2p7QnSVrsTry57hAKyqrxn/dmISxYL+n4hIiJ43isLzyDVVuPS/5364a80UPw\n",
       "zPzxXjv4xlsYloV+yp1QRA5E61efCk7A4zkOpk3r4Kr9Hob8+2lDFyG3wVVfi5Y1H8HdUOfxunLg\n",
       "IAQ9+BgUwQMkjoz0htdezVtbW7F48WK89957iIiI8Fa3pAvjkmPxp59NQ6BOnje24tNX8Is3N9Ks\n",
       "N+kzrlxrxQvvbcd7G47IlnDfmzsczy2Y4HMJ9600KekI+dkzHb7p244eQvO/3qATLAnpAZ7nYTt6\n",
       "EM3/er3DhFuTOhohj/6SEm4/5LVX9Keeegpz5szBnXfe6a0uSTcNHxyO156cgQiZZpvNtuuz3n9c\n",
       "tQfXqLQg8VMcx+PbgtP45VubcPKifBsCH509Go/NGQNWhqVjPaWMikHIz5+HKn6Yx+vX13m/DtvR\n",
       "Q/ShnJAucHYbWr/6FK3ffgne6fR4T8C0OQicv5S+QfJTDN/JK+GyZcvw6quvdtrBrl27UF1djT//\n",
       "+c8oLi6GRqMBz/NQKBRYs2YN5s+f3+5+o/HmztuzZ2nDjTc1mx1YsfUMrjRJWzv4Vjq1AvdkxWHc\n",
       "sHC/SBoIAYB6ow1f7q/C+Tr5NggrWAYPTozHuGHCmtg+z+2GrmQfNGfLO7zFEZ8ES9YUgJIFQgQU\n",
       "jVeh378VitZmj9d5lQrmnJlwDRoqcWTkxxITE9v+PTg4uEc/22nSfe3aNVy71vlpa3FxcXj66aex\n",
       "atUqsOzNiXO32w2WZZGTk4O9e/e2tVPSLS6L3YUPdlTiXG2rrHHERxowf8IQDAqjtd7Ed9mcbmw/\n",
       "fgV7ymrhkrju9q20KgWW5CVgRFyIbDF4g/psOXQl+8C4PS/LcQeGwDJxFtyhtASREAAAz0N9phS6\n",
       "owc6/nsTNADmSXeCCwmTODjiiWhJd3ddvnwZzc03P53xPI+RI0fijTfewD333IP4+Pi2a7cm3T0N\n",
       "lnSPy81h+bfF2Fx0TtY4WIZBfvYwPDwzHQaZ1px7Q3FxMQAgMzNT5kiIt/A8j4KyGry/8ajsp61G\n",
       "hwbg90umYHCUd14P5X5eXXWX0bJ2JdwN9R6vMwolAqbPhTZ7MhjWd9esE2nI/bzKiTO1wrR+dYen\n",
       "SwKAdtQ4BOTfD1YjT41+ItSbPNYr1UtiYmIQExMjaI+Li2uXcBNpKBUsnr53HIYOHIAV60skPznv\n",
       "Bo7nseHgWew7UY2fzh6F6WMTaMkJkV1NvRHLvyvB8XOeNylJadQdUfjdolyfKQnoDcqoGIQ8/hzM\n",
       "m/4N2/EiwXXe7YJp6zewV5xA4D2LoAijWW/Sv/A8D3v5UZg3/RucxezxHkalhmHuQmhH9b8PI32Z\n",
       "V0sGEt/BMAzmjE/EoIgg/L/P9qPV6uj6h0TSYrHjrX8fxpaic3hyXgYSB9FXZER6ZqsDq3eX45uC\n",
       "07J9EL3VvAlJeGzOGJ+uUHK7WI0WgfcuhmpoEkwb1gjKCgKAs6YKTe/+GQHTaNab9B+cqRWmjWth\n",
       "rzjR4T3KqBgELngEyvBICSMjUhAt6eY4OircF6TfEYXXf3EnXvl4Ly7WeT4+Viqna67h+Xe2YuLI\n",
       "ODw8Ix2xEUGyxkP6B4fTjY2HzmLN7pNosdjlDgdKBYun5mXgzizPFT/6Eu2oTChjB6N17Uq46i4L\n",
       "rvMumvUm/QPP83CUH4Np07oOZ7cBQDcuFwEz7wGjUkkYHZGKV9Z09wSt6ZaH1e7E62sKcfDkJblD\n",
       "AXC9UsPMjAQ8NC3N5w/W6c9rDv0Zx/HYebQKn20vxVUfKWUZHKDBC4snInWoeDNYvvi88k4nzNu/\n",
       "g/Xwvg7vYZRKmvXuh3zxefW27sxus/oAGOYuhGbEKAkjI7dD9jXdxPfpNCq8sHgSVu8ux+c7ysDJ\n",
       "XDPXzfHYXHQOO49ewN05SZg/ZYRfb7YkvoPneRyquISPtx5HdX2L3OG0SR4cht8+mIvIAQFyhyI5\n",
       "RqWCIf9+qFPSYfr2C7ibhFWxbsx628qOwDBnPlSxQ2SIlBDv4TkOtpIDsOzcCM7WcSlfTcooGObM\n",
       "B2sIlDA6Igea6e6Hyqrq8dcvD+Bai3z1vH/MoFPjvonJmDs+EQE+lnz3h5mYvoDneRw5cwVf7CrD\n",
       "qerOS51Kbf7kFDw8M12S9du+/rzyDjvMOzZ0OusNANox2QiYcRdYvUGiyIgcfP15vV3OmgswbVoH\n",
       "15XvO7yH1QfAkD8f6tTRYBgqMuAvepPHUtLdT7WY7fjftQdRdFq4zlJOeo0Kc7KH4e7c4RgQqJM7\n",
       "HAB9902hr+A4HgfKa7B2z0mcu9wkdzjtBAdo8NyC8cgYLqzuJBZ/eV4dFyo7nPW+gdXqoJ8+F9qx\n",
       "E2jJSR/lL89rd3GmVph3rIft2OFO76PZbf9FSTe5LRzH49sDp/HR5mM+Uc3hVioli5kZCbh/Ugqi\n",
       "QuWd6eprbwp9hdPlxq6jF7Bu70lcvibfSZIdGTk0Er95MAehQdJ+ePSn57W7s97KgYOuLzkZFC9N\n",
       "YEQy/vS8doZ3u68vJdm1qdOlJDS77f9oTTe5LSzL4N6JyRgxJAJ//qIAdU0d76iWmtPFYeOhSmwp\n",
       "OodJ6YMxf/IIxEf792l9xDusdie2FJ3D1/tP+dQSqRtYhsGi6Wl4IC+V6tJ3gVFrYMi/H5rU0TBt\n",
       "WAtX/RWP97mufI/m99+ENj0T+qn5UISEShwpIZ7xPA9n5SmYd6z3WKHnVtrRWQiYfhfNbvdjNNNN\n",
       "AFyvYfyPr4uwr7Ra7lA6lDY0AnOyEzEhNU7S2sZ9ZSbG312sbcamw5XYeaQKVodL7nA8igjW47mF\n",
       "4zEyIUq2GPz1eeXdbtiK9sO8ezN4u63D+xiFEtrMHOgnzQQbQOu9/Z2/Pq8A4Pz+Asw7NsB5obLT\n",
       "+5TRsTDkz4dq8FCJIiNioplu0msBOjX+a1EuctPi8M9vi2E0y1/P+MfKqq6irOoqBhi0mDXuDszO\n",
       "GoZwHy83SHrH5eZQWF6DjYfOoqzqqtzhdGpWZgIeyx/jcxuB/QWjUEA3fgo0qWNg3v4dbCeKPd7H\n",
       "u12wHtoL29FD0E3Ig25CHh2RTSTluloLy44NnR7fDvywJ2HaHGgzcmhPAgFAM93EA6PJhne/K8b+\n",
       "0hq5Q+kUyzDISonB3PFJSE+IEu2rfH+eifFXV5vN2FJ0DluLzqHJ1PGspy8ID9bjl/dlYWzSQLlD\n",
       "AdB3nldndRVMG9d2+ZU9qzdAP2UWtBkTwChoHsnf+NPz6jY2w7J7E2zHi4AuUiftmOzrS0no25g+\n",
       "hzZSElHsL63GP78p9olT/LoSHqzHlFFDkDc63utrv/3pTcGfma0OHCivwa6jF1B2ob6r9zSf4Iuz\n",
       "233peW3bnLZnS6en+AGAIngAdDnToB2TTaf5+RF/eF7dTddgPbALtqOHwLs7X9qmjB0Mw+z7aNNv\n",
       "H0ZJNxGNv8x63yo+Ohh5o+IxZXS8V5af+MObgr9yutwoPn0Ze45fxOFTl+B0cXKH1C2+Nrt9q774\n",
       "vHI2G6yFu2At3A3e6ej0XjbAAF32FGjH5YLV+kbZUdIxX35eXfVXYNm/A47yo+C5zl+bFOGRCJg2\n",
       "F+rkkVSVpI+jpJuIbn9pNVZ8V+LzX/XfimGAtPhITEofjKzk2Ns+bt6X3xT8kdPlRun5ehwor0FB\n",
       "WQ1M1s6TKF/CMMCdmXfgp7NH+9Ts9q368vPKmVph2bcVtpJC8G53p/cyGi10mbnQjZ9C1SJ8mC8+\n",
       "r86aC7Ds3w7HmfIu72WDQhCQdyc06ePAKBQSREfkRkk3kYTF5sTnO0rxXeEZn6vr3R2JsaHISolF\n",
       "VnIshg4M6fZshC++KfibVosdxacv41DFJRw5c8Vnq490JmlQKJ6cl4mkuDC5Q+lUf3he3Y0NMO/a\n",
       "BHvZkS7vZZRKaMdkQzc+D4rQcAmiIz3hK88rz3FwnjsNy4GdXVYjAa5vktRNmgnduIm0nKmfoaSb\n",
       "SKq6zoh3vy1GaVW93KHctohgfVsCnjY0EmpVxzMUvvKm4E94nsf3V1vaEu2Kiw3g/GGRtgdBeg2W\n",
       "zkrHzMw7/KLudn96Xl1Xvod59+ZuzUgCgPqOZGgzc6FOGkHVJHyE3M8rZzHDduwwbMUH4G5q6PJ+\n",
       "RqOFLmsidDnTaPlSP0VJN5Ecz/PYd6IaH2w66pMHlPSESslieFwYRg6NwsiESAyPC2+XhMv9puAP\n",
       "eJ7H5YZWlFbVo/R8HUrP1/vVUiRPGAaYPW4YlsxKR6BeI3c43dYfn1dX3WVY9u+Avfxol1UlAEAR\n",
       "HAJtRg60Y8bT0hOZyfW8Oi9dhK2oAPayo11ujgSuV8nRjZ8M7biJlGz3c5R0E9lY7U58sbMM3xSc\n",
       "9sslJ56olCyS48IxMiESaUMjYay9AK1a0a+SmK5wHI9LDS2ouNiAE+frUFZV7/cfvm41PC4MT87L\n",
       "QOIg315K4kl/TLpvcDc2wHJgF+zHDncrkWIUCqhTRkGXmQPl4ATaACcDKZ9X3mGHvfwYrMUFcF3u\n",
       "XnEARXDILVVxfHMfB5EWJd1EdpeutuDjbSdQUOY/VU66y2hsRmSQFuNHD8ew2FAMiw1FwsAB0Gv7\n",
       "xzq+Gwl25aXGtj/nrzTD5ofrsrsSNSAAP5kxElNGxfvFUhJP+nPSfYO71Qjbwb2wlhzo9HTLWykG\n",
       "hEGTnglN2lgowyNFjpDcIPbzyrvdcF44C/uJEtgrTnRZ/eYGRXgU9BOnQ5M2huq/k3Yo6SY+4+z3\n",
       "17Byy3EcP1cndyhe09zcDAAICblZ/5thgNjwQAyNHoBBEUGICQ9EbHggYsICfbaqRVdcbg71TWZc\n",
       "amjB5YZWXGpoxcW65j6bYN8qOECDh6alYXbWMCgV/r3Wl5LumzibFbYjB7u9XvcGZexgaEdmQJM6\n",
       "hpafiEyM55XnebiufA97aQnsZUfAmVq7/bO07p90hZJu4lN4nsfxc3X4aPMxnLvcJHc4veYp6e5M\n",
       "iEGL2B+S8KgBBoQG6RAaqENokA4DDFoE6jWyzKLaHS40tlrR2GJFk8mGa0YLGowWXL52PcGubTT1\n",
       "mSVC3aXXqHDfpGTckzscOk3f+OaCkm4hnuPgPH8G1uL9cJw52a113wDAsCxUCUnQpGVAnZgCVh8g\n",
       "cqT9j7eeV57n4W6og6PiBGylJXA3dH+jP6vVQTMmG7rMXKpwQ7rUmzyWvjMhXscwDEYPi8brT9+J\n",
       "A+U1+HjrcVy+ZpI7LMk0m2xoNtlQfuGqx+tKBYsBBi1Cg3QI0mug0yih06igU//wz1v+W6NWggHA\n",
       "sgwYhgEDgON58Pz1Nxmnm4PV7oTN4YLV7oLV7oTV8cM/7S6YbY4fEm0bLHanpL8HX6ZSspiTnYiF\n",
       "U0Yg2KCVOxwiMoZloR6WDPWwZLibG2E7UghbyUFwls5fl3iOg6PyFByVp8CwLJSDE6BOSoVmeBol\n",
       "Zz6Ad7vhrKmC43QZHGfK4W7s/rcZwPVvNHSZudCkjqGyf0QSNNNNROdyc9h5pApr95zElUb/S757\n",
       "OtNNfJdKyWLG2AQsmDICkQP65qwlzXR3D+9ywV5xHPZjh+E4f6bHP68Ij4JmeBrUSalQDhpCSxFu\n",
       "U0+fV85mhfPcadhPl8F59iQ4W882cDNaLTQpo6DNyIEqdnCP4yXEb5eXEEIIIYQQ4o96mnTTR3NC\n",
       "CCGEEEJERkk3IYQQQgghIpN8eQkhhBBCCCH9Dc10E0IIIYQQIjJKugkhhBBCCBGZpEn3ihUrMHXq\n",
       "VISEhIBlWVRXVwvuaWpqwpIlSxASEoKQkBAsXbqUKp4Qn5GXlweWZdv9Wbx4sdxhEdLmnXfewdCh\n",
       "Q6HT6ZCZmYn9+/fLHRIhAi+++KLgtTQmJkbusAgBAOzduxd33303Bg0aBJZlsXLlSsE9L774ImJj\n",
       "Y6HX6zF16lScPHmyy34lTbqtVitmz56Nl156qcN7Fi9ejGPHjmHLli3YvHkzjhw5giVLlkgYJSEd\n",
       "YxgGjz32GGpra9v+LF++XO6wCAEAfPnll3j22WexbNkyHDt2DDk5OcjPz0dNTY3coREikJyc3O61\n",
       "tLS0VO6QCAEAmM1mpKen480334ROpwPDtD9F+rXXXsPrr7+Ov//97ygqKkJkZCRmzpwJk6nzs0hk\n",
       "2UhZXFyMrKwsXLhwAYMH3yxOX1FRgdTUVBQUFGDChAkAgIKCAkyaNAmnTp1CUlKS1KES0s7UqVOR\n",
       "lpaGt99+W+5QCBHIzs7G6NGj230QTEpKwoIFC/Dqq6/KGBkh7b344otYt24dJdrE5wUGBuIf//gH\n",
       "li5dCuD6adAxMTH41a9+hRdeeAEAYLPZEBkZib/+9a944oknOuzLp9Z0FxYWwmAwtCXcAJCTk4OA\n",
       "gAAUFhbKGBkhN33xxReIiIhAWloafvvb33b5yZYQKTgcDhw5cgSzZs1q1z5r1iwcOHBApqgI6dj5\n",
       "8+cRGxuLhIQELFq0CFVVVXKHREiXqqqqUFdX1+61VqvVYvLkyV2+1irFDq4namtrERER0a6NYRhE\n",
       "RkaitrZWpqgIuWnx4sWIj49HTEwMysrK8MILL+DEiRPYsmWL3KGRfq6hoQFutxtRUVHt2un1k/ii\n",
       "8ePHY+XKlUhOTkZdXR1eeeUV5OTkoLy8HKGhoXKHR0iHbryeenqtvXz5cqc/2+uZ7mXLlgk2Q/z4\n",
       "z969e3s7DCGi6ckz/POf/xwzZ85EamoqHnzwQaxevRrbtm3D0aNHZf5fQQgh/mP27NlYsGAB0tLS\n",
       "MH36dGzYsAEcx3ncsEaIv/jx2u8f6/VM93PPPde2zqUjcXFx3eorOjoaV69ebdfG8zzq6+sRHR19\n",
       "2zES0pnePMNjx46FQqFAZWUlxowZI0Z4hHRLeHg4FAoF6urq2rXX1dVh4MCBMkVFSPfo9Xqkpqai\n",
       "srJS7lAI6dSNfLSurg6DBg1qa6+rq+syV+110h0WFoawsLDedgMAmDBhAkwmEwoLC9vWdRcWFsJs\n",
       "NiMnJ8crYxDyY715hktLS+F2uympIbJTq9XIyMjA1q1bMX/+/Lb2bdu2YeHChTJGRkjXbDYbKioq\n",
       "MG3aNLlDIaRTQ4cORXR0NLZu3YqMjAwA15/f/fv3469//WunPyvpmu4bZYHOnDkDACgvL0djYyOG\n",
       "DBmCAQMGICUlBbNnz8aTTz6JFStWgOd5PPnkk5g3bx4SExOlDJUQgfPnz+OTTz7B3LlzERYWhpMn\n",
       "T+LXv/41xo4di9zcXLnDIwTPP/88lixZgqysLOTk5ODdd99FbW0tnnrqKblDI6Sd3/zmN7j77rsR\n",
       "FxeH+vp6vPzyy7BarXjkkUfkDo0QmM1mnD17FgDAcRwuXryIY8eOISwsDHFxcXj22Wfx6quvIjk5\n",
       "GYmJiXjllVcQGBjY9bkdvIT+8Ic/8AzD8AzD8CzLtv1z5cqVbfc0NTXxDz/8MB8UFMQHBQXxS5Ys\n",
       "4Y1Go5RhEuJRTU0NP2XKFD4sLIzXaDT8sGHD+GeffZZvamqSOzRC2rzzzjt8fHw8r9Fo+MzMTH7f\n",
       "vn1yh0SIwEMPPcTHxMTwarWaj42N5RcsWMBXVFTIHRYhPM/z/K5duwT5KsMw/KOPPtp2z4svvsgP\n",
       "HDiQ12q1fF5eHl9eXt5lv7LU6SaEEEIIIaQ/8ak63YQQQgghhPRFlHQTQgghhBAiMkq6CSGEEEII\n",
       "ERkl3YQQQgghhIiMkm5CCCGEEEJERkk3IYQQQgghIqOkmxBCCCGEEJFR0k0IIYQQQojIKOkmhBBC\n",
       "CCFEZP8fqWnOY6OkNLQAAAAASUVORK5CYII=\n"
      ],
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7febf6e4ff98>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "circle1=plt.Circle((-4, 0), 5, color='#004080', \n",
    "                   fill=False, linewidth=20, alpha=.7)\n",
    "circle2=plt.Circle((4, 0), 5, color='#E24A33', \n",
    "                   fill=False, linewidth=5, alpha=.7)\n",
    "\n",
    "fig = plt.gcf()\n",
    "ax = fig.gca()\n",
    "\n",
    "plt.axis('equal')\n",
    "plt.xlim((-10, 10))\n",
    "plt.ylim((-10, 10))\n",
    "\n",
    "plt.plot ([-4, 0], [0, 3], c='#004080')\n",
    "plt.plot ([4, 0], [0, 3], c='#E24A33')\n",
    "plt.text(-4, -.5, \"A\", fontsize=16, horizontalalignment='center')\n",
    "plt.text(4, -.5, \"B\", fontsize=16, horizontalalignment='center')\n",
    "\n",
    "ax.add_artist(circle1)\n",
    "ax.add_artist(circle2)\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Here I have attempted to show transmitter A, drawn in red, at (-4,0) and a second one B, drawn in blue, at (4,0). The red and blue circles show the range from the transmitters to the robot, with the width illustrating the effect of the $1\\sigma$ angular error for each transmitter. Here I have given the blue transmitter more error than the red one. The most probable position for the robot is where the two circles intersect, which I have depicted with the red and blue lines. You will object that we have two intersections, not one, but we will see how we deal with that when we design the measurement function.\n",
    "\n",
    "This is a very common sensor set up. Aircraft still use this system to navigate, where it is called DME (Distance Measuring Equipment). Today GPS is a much more common navigation system, but I have worked on an aircraft where we integrated sensors like this into our filter along with the GPS, INS, altimeters, etc. We will tackle what is called *multi-sensor fusion* later; for now we will just address this simple configuration.\n",
    "\n",
    "The first step is to design our state variables. We will assume that the robot is traveling in a straight direction with constant velocity. This is unlikely to be true for a long period of time, but is acceptable for short periods of time. This does not differ from the previous problem - we will want to track the values for the robot's position and velocity. Hence,\n",
    "\n",
    "$$\\mathbf{x} = \n",
    "\\begin{bmatrix}x\\\\v_x\\\\y\\\\v_y\\end{bmatrix}$$"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The next step is to design the state transition function. This also will be the same as the previous problem, so without further ado,\n",
    "\n",
    "$$\n",
    "\\mathbf{x}' = \\begin{bmatrix}1& \\Delta t& 0& 0\\\\0& 1& 0& 0\\\\0& 0& 1& \\Delta t\\\\ 0& 0& 0& 1\\end{bmatrix}\\mathbf{x}$$\n",
    "\n",
    "The next step is to design the control inputs. We have none, so we set ${\\mathbf{B}}=0$.\n",
    "\n",
    "The next step is to design the measurement function $\\mathbf{z} = \\mathbf{Hx}$. We can model the measurement using the Pythagorean theorem.\n",
    "\n",
    "$$\n",
    "z_a = \\sqrt{(x-x_A)^2 + (y-y_A)^2} + v_a\\\\[1em]\n",
    "z_b = \\sqrt{(x-x_B])^2 + (y-y_B)^2} + v_b\n",
    "$$\n",
    "\n",
    "where $v_a$ and $v_b$ are white noise.\n",
    "\n",
    "We see an immediate problem. The Kalman filter is designed for linear equations, and this is obviously nonlinear. In the next chapters we will look at several ways to handle nonlinear problems in a robust way, but for now we will do something simpler. If we know the approximate position of the robot than we can linearize these equations around that point. I could develop the generalized mathematics for this technique now, but instead let me just present the worked example to give context to that development."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Instead of computing $\\mathbf{H}$ we will compute the partial derivative of $\\mathbf{H}$ with respect to the robot's position $\\mathbf{x}$. You are probably familiar with the concept of partial derivative, but if not, it just means how $\\mathbf{H}$ changes with respect to the robot's position. It is computed as the partial derivative of $\\mathbf{H}$ as follows:\n",
    "\n",
    "$$\\frac{\\partial \\mathbf{h}}{\\partial \\mathbf{x}} = \n",
    "\\begin{bmatrix}\n",
    "\\frac{\\partial h_1}{\\partial x_1} & \\frac{\\partial h_1}{\\partial x_2} &\\dots \\\\\n",
    "\\frac{\\partial h_2}{\\partial x_1} & \\frac{\\partial h_2}{\\partial x_2} &\\dots \\\\\n",
    "\\vdots & \\vdots\n",
    "\\end{bmatrix}\n",
    "$$\n",
    "\n",
    "Let's work the first partial derivative. We want to find\n",
    "\n",
    "$$\\frac{\\partial }{\\partial x} \\sqrt{(x-x_A)^2 + (y-y_A)^2}\n",
    "$$\n",
    "\n",
    "Which we compute as\n",
    "$$\n",
    "\\begin{aligned}\n",
    "\\frac{\\partial h_1}{\\partial x} &= ((x-x_A)^2 + (y-y_A)^2))^\\frac{1}{2} \\\\\n",
    "&= \\frac{1}{2}\\times 2(x-x_a)\\times ((x-x_A)^2 + (y-y_A)^2))^{-\\frac{1}{2}} \\\\\n",
    "&= \\frac{x_r - x_A}{\\sqrt{(x_r-x_A)^2 + (y_r-y_A)^2}} \n",
    "\\end{aligned}\n",
    "$$\n",
    "\n",
    "We continue this computation for the partial derivatives of the two distance equations with respect to $x$, $y$, $dx$ and $dy$, yielding\n",
    "\n",
    "$$\\frac{\\partial\\mathbf{h}}{\\partial\\mathbf{x}}=\n",
    "\\begin{bmatrix}\n",
    "\\frac{x_r - x_A}{\\sqrt{(x_r-x_A)^2 + (y_r-y_A)^2}} & 0 & \n",
    "\\frac{y_r - y_A}{\\sqrt{(x_r-x_A)^2 + (y_r-y_A)^2}} & 0 \\\\\n",
    "\\frac{x_r - x_B}{\\sqrt{(x_r-x_B)^2 + (y_r-y_B)^2}} & 0 &\n",
    "\\frac{y_r - y_B}{\\sqrt{(x_r-x_B)^2 + (y_r-y_B)^2}} & 0 \\\\\n",
    "\\end{bmatrix}\n",
    "$$"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "That is pretty painful, and these are very simple equations. Computing the Jacobian can be extremely difficult or even impossible for more complicated systems. However, there is an easy way to get Python to do the work for you by using the SymPy module [1]. SymPy is a Python library for symbolic mathematics. The full scope of its abilities are beyond this book, but it can perform algebra, integrate and differentiate equations, find solutions to differential equations, and much more. We will use it to compute our Jacobian!\n",
    "\n",
    "First, a simple example. We will import SymPy, initialize its pretty print functionality (which will print equations using LaTeX). We will then declare a symbol for NumPy to use."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/latex": [
       "$$\\phi$$"
      ],
      "text/plain": [
       "\\phi"
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "import sympy\n",
    "from sympy import  init_printing\n",
    "#from sympy.interactive import printing\n",
    "init_printing(use_latex='mathjax')\n",
    "\n",
    "phi, x = sympy.symbols('\\phi, x')\n",
    "phi"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Notice how we use a latex expression for the symbol `phi`. This is not necessary, but if you do it will render as LaTeX when output. Now let's do some math. What is the derivative of $\\sqrt{\\phi}$?"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/latex": [
       "$$\\frac{1}{2 \\sqrt{\\phi}}$$"
      ],
      "text/plain": [
       "   1   \n",
       "───────\n",
       "    ___\n",
       "2⋅╲╱ φ "
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "sympy.diff('sqrt(phi)')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We can factor equations."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/latex": [
       "$$\\left(\\phi - 1\\right) \\left(\\phi^{2} + 1\\right)$$"
      ],
      "text/plain": [
       "        ⎛ 2    ⎞\n",
       "(φ - 1)⋅⎝φ  + 1⎠"
      ]
     },
     "execution_count": 5,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "sympy.factor('phi**3 -phi**2 + phi - 1')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "SymPy has a remarkable list of features, and as much as I enjoy exercising its features we cannot cover them all here. Instead, let's compute our Jacobian."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/latex": [
       "$$\\left[\\begin{matrix}\\frac{x - x_{a}}{\\sqrt{\\left(x - x_{a}\\right)^{2} + \\left(y - y_{a}\\right)^{2}}} & 0 & \\frac{y - y_{a}}{\\sqrt{\\left(x - x_{a}\\right)^{2} + \\left(y - y_{a}\\right)^{2}}} & 0\\\\\\frac{x - x_{b}}{\\sqrt{\\left(x - x_{b}\\right)^{2} + \\left(y - y_{b}\\right)^{2}}} & 0 & \\frac{y - y_{b}}{\\sqrt{\\left(x - x_{b}\\right)^{2} + \\left(y - y_{b}\\right)^{2}}} & 0\\end{matrix}\\right]$$"
      ],
      "text/plain": [
       "⎡           x - xₐ                           y - yₐ              ⎤\n",
       "⎢ ──────────────────────────   0   ──────────────────────────   0⎥\n",
       "⎢    _______________________          _______________________    ⎥\n",
       "⎢   ╱         2           2          ╱         2           2     ⎥\n",
       "⎢ ╲╱  (x - xₐ)  + (y - yₐ)         ╲╱  (x - xₐ)  + (y - yₐ)      ⎥\n",
       "⎢                                                                ⎥\n",
       "⎢          x - x_b                          y - y_b              ⎥\n",
       "⎢────────────────────────────  0  ────────────────────────────  0⎥\n",
       "⎢   _________________________        _________________________   ⎥\n",
       "⎢  ╱          2            2        ╱          2            2    ⎥\n",
       "⎣╲╱  (x - x_b)  + (y - y_b)       ╲╱  (x - x_b)  + (y - y_b)     ⎦"
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "from sympy import symbols, Matrix\n",
    "phi = symbols('\\phi')\n",
    "phi\n",
    "\n",
    "x, y, xa, xb, ya, yb, dx, dy = symbols('x y x_a x_b y_a y_b dx dy')\n",
    "\n",
    "H = Matrix([[sympy.sqrt((x-xa)**2 + (y-ya)**2)], \n",
    "            [sympy.sqrt((x-xb)**2 + (y-yb)**2)]])\n",
    "\n",
    "state = Matrix([x, dx, y, dy])\n",
    "H.jacobian(state)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "In a nutshell, the entry (0,0) contains the difference between the x coordinate of the robot and transmitter A's x coordinate divided by the distance between the robot and A. (2,0) contains the same, except for the y coordinates of the robot and transmitters. The bottom row contains the same computations, except for transmitter B. The 0 entries account for the velocity components of the state variables; naturally the range does not provide us with velocity.\n",
    "\n",
    "The values in this matrix change as the robot's position changes, so this is no longer a constant; we will have to recompute it for every time step of the filter.\n",
    "\n",
    "If you look at this you may realize that this is just a computation of x/dist and y/dist, so we can switch this to a trigonometic form with no loss of generality:\n",
    "\n",
    "$$\\frac{\\partial\\mathbf{h}}{\\partial\\mathbf{x}}=\n",
    "\\begin{bmatrix}\n",
    "-\\cos{\\theta_A} & 0 & -\\sin{\\theta_A} & 0 \\\\\n",
    "-\\cos{\\theta_B} & 0 & -\\sin{\\theta_B} & 0\n",
    "\\end{bmatrix}\n",
    "$$\n",
    "\n",
    "However, this raises a huge problem. We are no longer computing $\\mathbf{H}$, but $\\Delta\\mathbf{H}$, the change of $\\mathbf{H}$. If we passed this into our Kalman filter without altering the rest of the design the output would be nonsense. Recall, for example, that we multiply $\\mathbf{Hx}$ to generate the measurements that would result from the given estimate of $\\mathbf{x}$ But now that $\\mathbf{H}$ is linearized around our position it contains the *change* in the measurement function. \n",
    "\n",
    "We are forced, therefore, to use the *change* in $\\mathbf{x}$ for our state variables. So we have to go back and redesign our state variables. \n",
    "\n",
    ">Please note this is a completely normal occurrence in designing Kalman filters. The textbooks present examples like this as *fait accompli*, as if it is trivially obvious that the state variables needed to be velocities, not positions. Perhaps once you do enough of these problems it would be trivially obvious, but at that point why are you reading a textbook? I find myself reading through a presentation multiple times, trying to figure out why they made a choice, finally to realize that it is because of the consequences of something on the next page. My presentation is longer, but it reflects what actually happens when you design a filter. You make what seem reasonable design choices, and as you move forward you discover properties that require you to recast your earlier steps. As a result, I am going to somewhat abandon my **step 1**, **step 2**, etc.,  approach, since so many real problems are not quite that straightforward."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "If our state variables contain the velocities of the robot and not the position then how do we track where the robot is? We can't. Kalman filters that are linearized in this fashion use what is called a *nominal trajectory* - i.e. you assume a position and track direction, and then apply the changes in velocity and acceleration to compute the changes in that trajectory. How could it be otherwise? Recall the graphic showing the intersection of the two range circles - there are two areas of intersection. Think of what this would look like if the two transmitters were very close to each other - the intersections would be two very long crescent shapes. This Kalman filter, as designed, has no way of knowing your true position from only distance measurements to the transmitters. Perhaps your mind is already leaping to ways of working around this problem. If so, stay engaged, as later sections and chapters will provide you with these techniques. Presenting the full solution all at once leads to more confusion than insight, in my opinion. \n",
    "\n",
    "So let's redesign our *state transition function*. We are assuming constant velocity and no acceleration, giving state equations of\n",
    "$$\n",
    "\\dot{x}' = \\dot{x} \\\\\n",
    "\\ddot{x}' = 0 \\\\\n",
    "\\dot{y}' = \\dot{y} \\\\\n",
    "\\dot{y}' = 0$$\n",
    "\n",
    "This gives us the the *state transition function* of\n",
    "\n",
    "$$\n",
    "\\mathbf{F} = \\begin{bmatrix}0 &1 & 0& 0\\\\0& 0& 0& 0\\\\0& 0& 0& 1\\\\ 0& 0& 0& 0\\end{bmatrix}$$\n",
    "\n",
    "A final complication comes from the measurements that we pass in. $\\mathbf{Hx}$ is now computing the *change* in the measurement from our nominal position, so the measurement that we pass in needs to be not the range to A and B, but the *change* in range from our measured range to our nominal position. \n",
    "\n",
    "There is a lot here to take in, so let's work through the code bit by bit. First we will define a function to compute $\\frac{\\partial\\mathbf{h}}{\\partial\\mathbf{x}}$ for each time step."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "from math import sin, cos, atan2\n",
    "\n",
    "def H_of(pos, pos_A, pos_B):\n",
    "    \"\"\" Given the position of our object at 'pos' in 2D, and two \n",
    "    transmitters A and B at positions 'pos_A' and 'pos_B', return \n",
    "    the partial derivative of H\n",
    "    \"\"\"\n",
    "\n",
    "    theta_a = atan2(pos_a[1]-pos[1], pos_a[0] - pos[0])\n",
    "    theta_b = atan2(pos_b[1]-pos[1], pos_b[0] - pos[0])\n",
    "\n",
    "    return np.array([[0, -cos(theta_a), 0, -sin(theta_a)],\n",
    "                     [0, -cos(theta_b), 0, -sin(theta_b)]])"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Now we need to create our simulated sensor. "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "from numpy.random import randn\n",
    "\n",
    "class DMESensor(object):\n",
    "    def __init__(self, pos_a, pos_b, noise_factor=1.0):\n",
    "        self.A = pos_a\n",
    "        self.B = pos_b\n",
    "        self.noise_factor = noise_factor\n",
    "        \n",
    "    def range_of(self, pos):\n",
    "        \"\"\" returns tuple containing noisy range data to A and B\n",
    "        given a position 'pos'\n",
    "        \"\"\"\n",
    "        \n",
    "        ra = math.sqrt((self.A[0] - pos[0])**2 + (self.A[1] - pos[1])**2)\n",
    "        rb = math.sqrt((self.B[0] - pos[0])**2 + (self.B[1] - pos[1])**2)\n",
    "        \n",
    "        return (ra + randn()*self.noise_factor, \n",
    "                rb + randn()*self.noise_factor)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Finally, we are ready for the Kalman filter code. I will position the transmitters at x=-100 and 100, both with y=-20. This gives me enough space to get good triangulation from both as the robot moves. I will start the robot at (0,0) and move by (1,1) each time step. \n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "image/png": [
       "iVBORw0KGgoAAAANSUhEUgAAAtsAAAFwCAYAAABpQ7cYAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\n",
       "AAALEgAACxIB0t1+/AAAIABJREFUeJzs3XlcVdX+//HXOcwCgiAzgkxi4gw45ixWZpa/bLLyevtm\n",
       "95aV6de8apbV9wZqWpnZdG9Z2TUry2z0OuaQmuGEs4CCM4LM83DO7w+KIhBH5vfz8eCRZ++99vns\n",
       "tnLeLNZey2A2m82IiIiIiMh1Z6zvAkREREREmiqFbRERERGRWqKwLSIiIiJSSxS2RURERERqicK2\n",
       "iIiIiEgtUdgWEREREaklCtsiIiIiIrWkxrC9aNEiunTpgpOTE05OTvTp04fvv/++0jHPP/88Pj4+\n",
       "tGjRgkGDBnHw4MFaLVhEREREpLGoMWy3adOGuXPnsnv3bnbu3MngwYO544472Lt3LwBz5szhlVde\n",
       "4Y033uCXX37B3d2dqKgocnNz66R4EREREZGGzHClK0i6uroye/ZsHn74Yby9vXnyySeZPn06AIWF\n",
       "hbi7uzNv3jweeeSRWilYRERERKSxuOwx22VlZSxbtozCwkL69+/P8ePHSUlJYdiwYRXH2Nra0r9/\n",
       "f7Zu3VorxYqIiIiINCaWlzpg37599O7dm6KiIuzs7Pjss88IDQ2tCNQeHh6Vjnd3d+fMmTO1U62I\n",
       "iIiISCNyybDdvn174uLiyMrK4vPPP+fee+9lw4YNNbYxGAxVtmVlZV19lSIiIiIi9czJyemK21xy\n",
       "GImVlRWBgYF069aN6OhoevXqxaJFi/Dy8gIgJSWl0vEpKSl4enpecSEiIiIiIk3NFc+zXVZWhslk\n",
       "IiAgAE9PT1avXl2xr7CwkC1bttCnT5/rWqSIiIiISGNU4zCSadOmMWLECHx9fcnJyWHp0qVs3LiR\n",
       "VatWAfDUU08RHR1N+/btCQkJ4Z///CeOjo6MGTOmxje9mi54aThiY2MBiIiIqOdK5FrpXjYNuo9N\n",
       "h+5l06F72XRc61DoGsN2SkoKDzzwAOfOncPJyYkuXbqwatUqoqKiAJg6dSoFBQVMmDCBjIwMevXq\n",
       "xerVq7G3t7+mokREREREmoIaw/bixYsveYJZs2Yxa9as61aQiIiIiEhTccVjtkVERERE5PIobIuI\n",
       "iIhIo1dSWkKZqay+y6hCYVtEREREGq2c/CxW/fwpzy8eT1zi9voup4pLLmojIiIiItLQnL1wgh93\n",
       "f0Ps4Y2UlBUDsGH313QL6VvPlVXW4MK22WymuLgYs9lc36XIHxgMBqytratdHVRERESkLpjNZg6f\n",
       "2MOG3V9zOHl3lf1JZ4+QdO4obT3b1UN11WtQYdtkMlFUVIS1tTUWFhb1XY78QVlZGYWFhdjY2NR3\n",
       "KSIiItLMFJcWEXt4Iz/u/oZz6ScrtptNZo7tP8ep+FQG3NkZGytbzmecVti+mOLiYmxtbdV72gBZ\n",
       "WFhga2tLUVFRfZciIiIizUR2Xgab435gy75V5BVkV2w3lZmI332a2LXxpJ/LAeDhvz7Cw2OewM6m\n",
       "Ya330qDCNqCg3YDp3oiIiEhdOJ16nA27v2bn0c2UlZVWbC8tKePwjpPsXB9P9oV8ANw93Jg2bToP\n",
       "3f1Igwva0ADDtoiIiIg0PyaziQPHY/lx9zfEn9pXaV9xUSn7tyaxZ0MCednlv2VvG+DPzGee5cEH\n",
       "H8Ta2ro+Sr4sCtsiIiIiUm+KSgr5+eB6Nu75ltTMM5X2FeYVE7f5OHs3JVKYXwJAWMcOPDvzOUaP\n",
       "Ht0onvFT2K5nSUlJBAYGsnjxYv7yl78A8MEHH/DQQw+RlJSEn59fPVcoIiIicv1l5KSxee/3bN2/\n",
       "mvyi3Er78rIK2bMxkX0/HaekqHyhmp69evLszGcZPnx4oxraqrBdB34Lz9W59dZbMRgMl/xLs3Tp\n",
       "UlJTU5k4cWJtlCgiIiJSJ5LPxfPj7q/ZnbAV059WfMy6kMfu9Qkc/PkEZaUmAKKionjmmWfo379/\n",
       "owrZv1HYrkMvvPACQUFBlbaFhobyxRdfYGlZ861YunQpBw4cUNgWERGRRsdkKmPfsR1s2P01x84c\n",
       "qrI//Vw2sWvjObrrNGZT+Voro0aNYsaMGURERNR1udeVwnYduummm+jRo8dVt6+Nn+YKCgqws7O7\n",
       "7ucVERERKSjKZ/vBtWza8x0XslOq7E85kUHs2niOxZ0FyqcaHvPgGKZNm0aHDh3qutxaYazvApq7\n",
       "pKQkjEYjH3744UWPGThwIN9//33Fsb99/cZsNrNw4UI6deqEnZ0dHh4ePPzww1y4cKHSedq2bcst\n",
       "t9zCunXr6NmzJ3Z2dsydO7fWrk1ERESapwvZKXy56X1mvf8wKza9Xylom81mTieksfKtrXz2yiaO\n",
       "xZ3FxsaGRx99lPj4eD766KMmE7RBPdt1KjMzk7S0tGr31dRrPXPmTKZOncqpU6d47bXXqux/9NFH\n",
       "ef/99xk3bhxPPvkkJ06cYOHChezYsYNffvmlYtVHg8FAQkICd911F4888gjjx4/XA5giIiJyXZjN\n",
       "Zo6fPcKG3SuJS/wZs9lUZX/SwRR2ro3n7PF0ABwcHHj00UeZNGkSXl5e9VF2rWvUYfvJBXfU6vlf\n",
       "n/jVdT3fzTffXOm1wWAgLi7uku2GDh2Kt7c3mZmZjBkzptK+rVu38u6777JkyRLuv//+Su/Vr18/\n",
       "PvroI8aPHw+U/yVPTEzk66+/ZsSIEdfhikRERKS5KysrZU/CNn7c/TXJKfFV9ptMZhL3nmHn2nhS\n",
       "T2cB4OLiwsSJE3n88cdxcXGp65LrVKMO243NwoULueGGGypts7W1vaZzfvbZZzg4ODBs2LBKveah\n",
       "oaG4u7uzYcOGirAN0KZNGwVtERERuWb5hbls3b+aTXu/IzP3QpX9ZaUmjsSeZPf6RNLPly+p7uXl\n",
       "xf/+7//yt7/9DQcHh7ouuV4obNehyMjIKg9IJiUlXdM5jx49Sm5uLh4eHtXuT01NrfQ6MDDwmt5P\n",
       "REREmrfUzLNs3PMN2w+up7iksMr+kuJSDm5LZu/GJLLSy+fPDgwM5B//+Adjx4695o7GxkZhu5Ez\n",
       "mUy4urry6aefVru/VatWlV5r5hERERG5WnGJ21n8wzzKykqr7CvKLyHup+Ps33KC3Kw8AMLCwpg+\n",
       "fTr33HPPJac5bqoa9VVf7zHVDdnFHqAMCgpi7dq19OzZE3t7+zquSkRERJqTNbFfVgna+TlF7Nuc\n",
       "xL6fkijIK+/pjoyM5JlnnuG2226rNINac9S8r74Rsbe3JyMjo8r2e++9F5PJxIsvvlhlX1lZGZmZ\n",
       "mXVRnoiIiDQD2Xm/Z5GcjAK2f3OUj/+5nh2rD1OQV8jgwYNZu3YtP//8M7fffnuzD9rQyHu2m5PI\n",
       "yEg+++wznnrqKXr06IHRaOTee++lX79+TJgwgZdffpm4uDiGDRuGjY0NCQkJfPHFF/zf//0fY8eO\n",
       "re/yRUREpAkoKysl43wuu9bFczj2JKay8tUeb7vtNmbMmEGvXr3qucKGR2G7jlzp6o9/Pv6xxx5j\n",
       "3759fPzxxyxcuBAo79WG8llOunfvzttvv83MmTOxtLTE39+fe+65h8GDB191DSIiIiK/2bt3L1+8\n",
       "s4HDu5LBDAYD3HX3aGY+8yydO3eu7/IaLIPZbDbXxRtlZWVV/NnJyanaYwoLC5vdE6qNTWFhIfv3\n",
       "7wcgIiKinquRaxUbGwvoXjZ2uo9Nh+5l09GU7uXWrVuJjo7mu+++A8BoYeCGHn50HxzMO89+jZ1N\n",
       "i3qusHZdToatiXq2RURERKQSs9nMmjVriI6OZuPGjQC0aNGC9j196DIgAAfn8tnNLCws6rPMRkGj\n",
       "1kVEREQEKJ9S+MsvvyQyMpKbbrqJjRs34uTkxMyZM0lKSuLGOzpWBG0AS6P6bS9F/4dEREREmrmS\n",
       "khKWLVtGTEwMhw4dAsDd3Z3Jkyfz6KOP0rJlS0ymMsxmU0UbAwaMRvVsX4rCtoiIiEgzVVhYyOLF\n",
       "i5k7d27FqtZ+fn5MnTqVhx56qNJieGWmskptLSwUIy+H/i+JiIiINDM5OTm8/fbbzJ8/n5SUFABC\n",
       "Q0OZPn06Y8aMwcrKqkqb0j8tZmOhXu3LorAtIiIi0kxcuHCB119/nYULF1YsltetWzdmzJjBqFGj\n",
       "LvrAY0lpMQeO/1Jpm4VF1UAuVSlsi4iIiDRxZ86cYf78+bzzzjvk5eUB0K9fP2bMmMFNN91U7Voc\n",
       "JlMZ8af2E3tkE3EJ2ygozq+0Xz3bl0dhW0RERKSJOnbsGHPmzOGDDz6guLgYgFtuuYXp06fTr1+/\n",
       "KsebzWaSzh1l19HN7Dq6hZz8zIue29vVv9bqbkoUtkVERESamP379zN79mw++eQTTCYTBoOB0aNH\n",
       "M2PGDLp161bl+LMXTrDzyCZ2HtnMheyUS57fw8WXO/r9tTZKb3IUtkVERESaiB07dhAdHc3KlSsB\n",
       "sLS0ZOzYsfzjH/+gffv2lY69kJ3CriNb2HlkE2cuJF/y3A52TnQL6Ut4aH8CvEKrHXoiVSlsi4iI\n",
       "iDRiZrOZDRs2EB0dzbp16wCwtbXl4YcfZsqUKfj7/z7cIyc/k93xP7HzyGaOnz18yXPbWNvRJagX\n",
       "4aH9adems8ZpXwWFbREREZFGyGQy8d133xEdHc327dsBcHR0ZMKECTz11FN4eHgAUFCUT1zidnYe\n",
       "2cTRk3GY/rAwTXUsLawIC4ggvF0/OgSEY21pU+vX0pQpbDczP/74I4MHD2bZsmXcfffd9V2OiIiI\n",
       "XKHS0lI+//xzYmJi2LdvHwCurq5MmjSJCRMm4OzsTElpMXvit7Lz6GYOHI+ltKykxnMaDEZC23Qm\n",
       "PLQ/nYN6YmdjXxeX0iwobNcBo9F4WcctXryYv/zlL7VcjYiIiDRGRUVFfPTRR8yZM4fExEQAfHx8\n",
       "mDJlCuPHj8fWzpajJ+P4dscm4hJ/pvBPU/VVJ8CrPeGh/ega3JeW9s61fQnNksJ2Hfj4448rvX7n\n",
       "nXfYvn07ixcvrrS9T58+dVmWiIiINAJ5eXm8++67zJs3jzNnzgAQFBTEtGnTeOCBBziTnsQPv/yH\n",
       "3Ud/Iqcg65Ln83b1Jzy0P91Db8S1pUdtl9/sKWzXgTFjxlR6vXr1anbs2FFl+5/l5eVhb69f44iI\n",
       "iDRHGRkZLFq0iNdee40LFy4A0KlTJ6ZNm0avAd3Zn/QLsz95kvTs85c8l2tLD8JD+9G9XT+8W2t+\n",
       "7Lp0eeMbpNaNGzcOOzs7kpOTGTlyJE5OTowYMQKAuLg4/vrXvxIUFISdnR1ubm7cd999nDx5ssp5\n",
       "srKyePrppwkMDMTW1hZfX1/uv//+ip+Eq1NSUsJdd92Fg4NDxVPMIiIiUj9SUlKYNm0a/v7+PPvs\n",
       "s1y4cIFu3bvwwvxp/P2FO9md9Q2vLZ/G2tgvagzaji2cGdB1BJPunsNz495mRJ8HFLTrgXq2GxCT\n",
       "ycSwYcPo2bMn8+bNw9Ky/PasXbuWo0ePMm7cOLy9vUlISODtt99mx44d7N+/Hzs7O6C8J3zAgAEc\n",
       "OHCAv/71r0RERJCWlsYPP/xAYmIi3t7eVd6zqKiI0aNHs3nzZv773//St2/fOr1mERERKZecnMzL\n",
       "L7/Me++9R2FhIQDtOgfQeYAfnoFOpBkOkXai5nPYWbegS3BvwkP7E+zbUVP1NQAK2w1ISUkJt912\n",
       "G/Pmzau0/dFHH2Xy5MmVto0cOZK+ffvy5Zdfcv/99wPw8ssvExcXx+eff86dd95ZceyMGTOqfb/8\n",
       "/Hxuv/12du3axZo1a4iMjLzOVyQiIiKXsjtuJy+8OItvv1pFWVkZAIGdPAkf2g5P/1aXbG9lYU1Y\n",
       "YATh7frToW13rCyta7tkuQI1hu2YmBi+/PJLjh49io2NDb169SImJoawsLCKY8aNG8dHH31UqV2v\n",
       "Xr3YunVr7VT8B7W9cpHZbK7V81fnscceq7Ltt55rgNzcXIqKiggJCcHZ2Zldu3ZVhO3ly5fTsWPH\n",
       "SkH7YrKzs7n55ps5cuQIGzZsoHPnztfvIkREROSicvIzSTx9kNUbfmDJv5exPzYBzGAwGgiN8CV8\n",
       "SAiuXi1rPIeDnRNB3jfQKagnnQJ7YmfToo6qlytVY9jeuHEjjz/+OJGRkZhMJp577jmGDh3KwYMH\n",
       "adWq/Cctg8FAVFQUS5YsqWhnba2fqK6G0Wikbdu2VbZnZGQwbdo0li9fTkZGRqV9WVm/P3WcmJjI\n",
       "qFGjLuu9Jk+eTEFBAbt27aJTp07XVLeIiIhcXEZOKgmnD5J4+gCJpw+yK3YPsWviOXG4fLy10cJI\n",
       "h55+dB8cjFPr6idGcHZwJdinI0E+HQj2CcO9lY+WS28kagzbq1atqvR6yZIlODk5sXXrVm699Vag\n",
       "vPfX2toad3f32qvyIuqj57k2WVtbVzsn9913383WrVuZMmUK3bp1w9HREYB7770Xk+n3VaCu5B/d\n",
       "HXfcwbJly3jppZdYunTpZc8FLiIiIhdnNptJzTxLfMpuUrJO8O2+d0nPPo/ZbCb50Hli1x7l7LF0\n",
       "AKysLejYty1dBwbh4GRX6Txuzt4VwTrIpwMuju4K143UFY3Zzs7OxmQyVfRqQ3nA27JlCx4eHjg7\n",
       "OzNgwABeeukl3NzcrnuxTV11PzxkZGSwbt06XnjhBZ599tmK7YWFhaSnp1c6NigoqGIlqUsZMWIE\n",
       "w4cP54EHHsDe3p733nvv2ooXERFphkxmE+cunCTx9AESfu25zs7//bfQJpOZxLgzxK6JJ+10+W+j\n",
       "bVpY0aVfIJ37B2JnXz4awNvVn6Bfg3WQTwec7F3q5Xrk+ruisD1x4kS6detG7969K7bdfPPN3Hnn\n",
       "nQQEBHD8+HFmzpzJ4MGD2blzp4aT1KC6n06r22ZhUf4U8R97sAFeffXVKuF89OjRvPDCCyxfvpzR\n",
       "o0dfsoZ7772XvLw8xo8fj4ODAwsWLLiSSxAREWl2TGYTp84f+zVYHyDxzCHyC3OqHFdWauLIzpPs\n",
       "XBtPZmoeAC1a2tBtYDCd+wYQ6NeeYJ8OBPmEEejVHnu7msdoS+N12WF78uTJbN26lS1btlQKhffc\n",
       "c0/Fn8PCwggPD8ff35/vvvvuouOHY2Njq93u7++Pra3t5ZbUqFXXi13dtpYtWzJw4EDmzp1LcXEx\n",
       "fn5+bNmyhU2bNuHq6lqpzdNPP80XX3zBfffdx+rVq+nevTuZmZmsWrWKF198kf79+1c5///8z/+Q\n",
       "m5vLpEmTcHBw4KWXXqqx7pyc37+hXOw+SuOje9k06D42HbqXDVNxaSFrD3xCWu7pix5TUlzKwe0n\n",
       "2LU+gdzMAgBaurRgwK09uPmWKHxaB+Hu6IuVpQ0ARelwKP1ondQvVyckJOSa2l9W2J40aRKfffYZ\n",
       "GzZsqPYBvj/y8vLC19eXhISEayqsKTMYDFV6savb9pulS5cyceJE3nnnHUpKShgwYADr169n6NCh\n",
       "ldq0aNGCTZs28fzzz/Pll1/y4Ycf4uHhwYABA2jXrl2l9/qjiRMnkpOTw3PPPYejoyPTpk27jlcr\n",
       "IiLSNBxL3X/RoF1UUMK+LcfZszGRgtxiALx8PbjvgXsYNWI0tjZ21baTps9gvsRThhMnTuTzzz9n\n",
       "w4YNhIaGXvKEqamp+Pr68t577/HAAw9UbP/jrBlOTk7Vti0sLGw2PduNVWFhIfv37wcgIiKinquR\n",
       "a/Vb75nuZeOm+9h06F42bLGHN/LRf1+ttK0gt4j9m0+wZ3MihflFAEREhHP33fcwYMAAevToUR+l\n",
       "ynV0ORm2JjX2bE+YMIGPP/6Yr776CicnJ86dOweAo6Mj9vb25OXlMWvWLEaPHo2npydJSUlMnz4d\n",
       "Dw+Py56CTkRERKQx6BTUE3u7luQVZJOTUcDuDQkc+vkkxUUlAAwcOJAZM2YwdOhQdu7cWc/VSkNR\n",
       "Y9h+6623MBgMDBkypNL2559/nueeew4LCwv279/PkiVLyMzMxMvLi8GDB7N8+XLs7aufJ1JERESk\n",
       "MbKxsiW0dQ9efvllDv9yElNZ+eCAW4bfwsxnZtKnT596rlAaohrD9p9nwPgzW1vbKnNxi4iIiDQ1\n",
       "e/fuJSYmhs8//xyTyYTBACHdfAgfGsJfR/+VPj0UtKV6VzT1n4iIiEhzsm3bNqKjo/n2228BsLKy\n",
       "YtiIgbh1MuDs5gDAht1fM6DrCGys9RCkVKVlA0VERET+wGw2s2bNGgYNGkSfPn349ttvsbOzY+LE\n",
       "iSQmJvLFZyvx9Pl95ey8why27NNv+qV6CtsiIiIilA+fXbFiBT169GDYsGH8+OOPODk58cwzz5Cc\n",
       "nMxrr72Gr68vhUX5BPt2rNR2/c6vKC4pqqfKpSHTMBIRERFp1kpLS/nkk0+YPXs2Bw8eBMDNzY1J\n",
       "k57izvtuJ6swjZ8Of8upTcc4lZZU7YqROQVZ7Du2g/DQfnVdvjRwDS5sm83miy7uIvXrElOyi4iI\n",
       "NCqFhYV88MEHzJkzh6SkJAA8PN24ZfQgQiN9SMvZyxsrd1z2+aoL4SINKmxbW1tTWFiItbU1FhYW\n",
       "9V2O/EFZWRnFxcXY2NjUdykiIiLX5HxaCq8seJl/vf0+6WkZADi7OxA+JJjQ8DZYWBZxJuPYFZ3T\n",
       "3yOEiPYDa6FaaewaVNg2Go3Y2tpSXFxMSUlJfZcjf2AwGLC1tdVvHUREpFHJzsvkVOoxTp1P5HDi\n",
       "flZ8+h3b1sRRlF+eM1r7OBERFUJQZ2+Mxsv7jLOxssXHLQBft8DyL/cAfFoH6DNSqtWgwjaUhzr1\n",
       "noqIiMiVMJvNXMhO4dT5Y5xKPV4esFOPkZ2XQV5WIbt/TGD/T0mUFJcB4BXgQkRUO/xvcK8xJDvY\n",
       "OeHrFoCve1D5f90Cae3sidGgOSbk8jS4sC0iIiJSkzJTGSnpp34N1OXB+nTqcQqK8iodl5WWx671\n",
       "CRz8+QSmsvKF+vzauxMR1Q6fINcq53Vp6f5rb/WvvdbugTjZu6jHWq6JwraIiIg0WMWlRZxNS+ZU\n",
       "6nFOnk/kVOpxzqYlU1JWfNE2F85ms3NdPEd3ncZsMoMBgrp4ETG0He5tnDEYjHi08qkYAuLrFoiP\n",
       "WwD2to51eGXSXChsi4iISIOQX5TLqfO/DwE5nXqclPRTmMymy2qfciKD2DVHObbvHAAGo4EOPfwZ\n",
       "ftcgwrtF4PPrGGuf1m2xttKQVakbCtsiIiJSp8xmM9l5GRWh+rdx1heyU67qXKcTLhC79ignj6QC\n",
       "YGVlye133cqU/51CRJdeWFgo7kj90d8+ERERqVWlZSUcOB7LiZQETqYe4/T5Y+QUZF3TOR3snMhK\n",
       "LmX91zs4vD+hfJuDA4899hiTJk3C09PzepQucs0UtkVERKTWlJQWM/eTyaSkn7rqc7g6eVRMs+fl\n",
       "4kfsln0seHUhcXFxALi4uPDUU0/x+OOP06pVq+tVush1obAtIiIiteboybjLDtpGgxEPF9/K81e7\n",
       "BdDCxoHi4mKWLFnC47MfICGhvCfb29ubKVOmMH78eBwcHGrzMkSumsK2iIiI1Bpnh9bVbreysMa7\n",
       "tX/FFHu+boF4tfbD2rLyg4t5eXkseHsB8+bN49Sp8tAeGBjItGnTGDt2rNbmkAZPYVtERERqjXdr\n",
       "f7xbt+VMWlLFtv5dhjOq//9gYbS4aLvMzEwWLVrEa6+9RlpaGgAdO3Zk+vTp3H333VhaKsJI46Dl\n",
       "j0RERKTWGAwGeocNrbTtwPGdF10o5vz588yYMQN/f39mzpxJWloaPXr0YOXKlezdu5cxY8YoaEuj\n",
       "orAtIiIitSqi/QAsLawqXl/ITiHh1P5Kx5w4cYInn3wSf39/YmJiyM7OZsiQIaxbt47t27czcuRI\n",
       "jEbFFml89LdWREREapW9rSNdgnpV2rZt/xoAjhw5wkMPPURQUBALFy6ksLCQ22+/ne3bt7N27VoG\n",
       "Dx6s5dKlUdPvYURERKTW9e4Yxc6jmyter/nxB1b++ye+WvEVZrMZo9HImDFjmDZtGp06darHSkWu\n",
       "L4VtERERqXXBvh1xbenBvj0HiV1zlORD5wGwtrZm3LhxTJ06laCgoHquUuT6U9gWERGRWmU2m1mz\n",
       "eg3LX9/M/j2HALC0tqDXkM588u7X+Pr61nOFIrVHYVtERERqhclkYsWKFURHR7Nr1y4AbOys6Nwv\n",
       "gC79A7FzsMFkVVjPVYrULoVtERERua5KSkpYunQps2fP5vDhwwB4eHgwefJkrP2yOZby+0wk2/av\n",
       "wc8juL5KFal1mo1ERERErouCggLefPNNQkJCGDduHIcPH8bf359FixZx4NA+Ioa2IyPvXKU2O49u\n",
       "xmQ21VPFIrVPPdsiIiJyTbKzs3n77bd55ZVXSElJAaB9+/ZMnz6dgcP6su3AamKWPk5xaVGVtsWl\n",
       "RZSWlmBtpWXXpWlS2BYREZGrkpaWxuuvv87ChQvJzMwEoHv37kyfPp3gLj5sifueeZ9+WeM5hkWO\n",
       "VtCWJk1hW0RERK7I6dOnmT9/Pu+88w75+fkA9OvXj/99ehL2Xma27FvFpu/PX7S9pYUV4aH96d/l\n",
       "Vtq4B9ZV2SL1QmFbRERELktCQgJz587lww8/pLi4GIDhw4cz/tFxFNqdZ9PhTyg5VnzR9s4OrtzY\n",
       "+RZ6h0Xh2MKprsoWqVcK2yIiIlKjffv2MXv2bJYtW4bJZMJgMHDXXXdx5/23klIaz/rE/9TYPsi7\n",
       "A/27jqBzUE8sjBZ1VLVIw6CwLSIiItXavn07MTExfP311wBYWlrywIP3M3BEJEnZcfyUtOKiba0s\n",
       "rAkP7Uf/rrfi66ahItJ8KWyLiIhIBbPZzPr164mOjmb9+vUA2NnZcd8D99B1YCBJGfvYfXrdRdu3\n",
       "cmhdPlSkYxQOdi3rqmyRBkthW0RERDCZTHzzzTdER0ezY8cOAFq2bMld999BQHgrUnKTiE/dedH2\n",
       "wT5h9O9yK500VESkEoVtERGRZqy0tJRPP/2UmJgYDhw4AIBra1dG3jUMtxssKDBlkZKbVW1bKwtr\n",
       "ItoPoH+X4fi4BdRl2SKNhsK2iIhIM1RUVMQHH3zA3LlzOXbsGABe3p4Mvb0vzsEmsMin4CILO7Zy\n",
       "dKNf51voHTYUew0VEamRwraIiEgzkpuby7vvvsu8efM4e/YsAH5t29Dn5k64BFtgYVl60bbBvh0Z\n",
       "0GUEHQOFVTzcAAAgAElEQVQjNVRE5DIpbIuIiDQD6enpvPHGGyxYsID09HQAAkP86To4AK/QlhiN\n",
       "hmrbWVlaE9l+AP273Ip367Z1WLFI06CwLSIi0oSdPXuWV199lbfeeovc3FwAQjoEENbfG99QVwyG\n",
       "6kO2i6Mb/boMp1fYUOxtHeuyZJEmRWFbRESkCTp+/Dgvv/wy77//PkVFRQC069yWjv198A66eMgO\n",
       "8e3EgK630jEgEqOGiohcM4VtERGRJuTgwYPMnj2bpUuXUlZWBkD77v50HuiPh1+ratuUDxUZ+OtQ\n",
       "Ef+6LFekyVPYFhERaQJiY2OJjo5mxYryVR2NRiM39PCj++AgXDyrnzHEpaU7/ToPp1fYEA0VEakl\n",
       "CtsiIiKNlNlsZtOmTURHR7N69WoALK0suKFHG7oPDqGla4tq27Vr05kBXUcQ1jZcQ0VEapnCtoiI\n",
       "SCNjNpv5/vvviY6OZuvWrQBY21jSsW9bug4Iwt7Jtkoba0sbIm8YRP8uw/Fy9avrkkWaLWNNO2Ni\n",
       "YoiMjMTJyQl3d3dGjhxZsbrUHz3//PP4+PjQokULBg0axMGDB2utYBERkeaqrKyMTz/9lG7dujFi\n",
       "xAi2bt2KbQtret7cnr88F0XfkWFVgnZrJ09G9X+IFx9+j3sG/11BW6SO1dizvXHjRh5//HEiIyMx\n",
       "mUw899xzDB06lIMHD9KqVflDFnPmzOGVV17hww8/pF27drz44otERUVx5MgRHBwc6uQiREREmrLi\n",
       "4mI+/vhjYmbHkBCfAIB9Sxu6DQomrE9brG2qfpx38O9Ovy7DuaFtd4yGGvvWRKQW1Ri2V61aVen1\n",
       "kiVLcHJyYuvWrdx6662YzWZee+01pk+fzqhRowD48MMPcXd3Z+nSpTzyyCO1V7mIiEgTl5+fz7//\n",
       "/W/mzJ3DmdNnAGjp2oLwISG0j2yDpVXl8da21i3o2WEw/ToPx72Vd32ULCJ/ckVjtrOzszGZTBW9\n",
       "2sePHyclJYVhw4ZVHGNra0v//v3ZunWrwraIiMhVyMzMZNGiRcx/ZT4Z6RkAuHg6EjE0hJBuPhgt\n",
       "KvdUe7j40r/zcCJvGISttV19lCwiF3FFYXvixIl069aN3r17A3Du3DkAPDw8Kh3n7u7OmTNnrlOJ\n",
       "IiIizUN6ejr/WfofvvhiOXm5+QC4+zkTGdWOgDBPDH9YUt1gMNIxIIL+XW6lXZvOF12kRkTq12WH\n",
       "7cmTJ7N161a2bNlyWf+gazomNjb2ct9WGjDdx6ZD97Jp0H1svM6dO8f7H7zHt998R0lxCQC+Ia2J\n",
       "iGqHb0jrSp+p1pa2hHh0I9QzHAdbZ3LOl7Dz/M76Kl0uQf8uG7+QkJBran9ZYXvSpEl89tlnbNiw\n",
       "gbZt21Zs9/T0BCAlJQVfX9+K7SkpKRX7REREpHrJycm8+97brF29HlOZCYC2YR5EDG2HV4BLpWNb\n",
       "2XvQ3iuCgNYdsbSwqo9yReQqXDJsT5w4kc8//5wNGzbQrl27SvsCAgLw9PRk9erVhIeHA1BYWMiW\n",
       "LVuYN2/eRc8ZERFxjWVLffrtp3Tdx8ZP97Jp0H1sfPbs2cOsF57jm5XfYjabMRigXXcfwoeG0Nrb\n",
       "qeI4o8FI5+BeDOhyK4HeHTRUpBHRv8umIysr65ra1xi2J0yYwMcff8xXX32Fk5NTxRhtR0dH7O3t\n",
       "MRgMPPXUU0RHR9O+fXtCQkL45z//iaOjI2PGjLmmwkRERJqan376iRf/7wVW/3cNAEYLAx16+NN9\n",
       "cDDObr9Pl+tg50TfTsPo2+lmnB1c66tcEbkOagzbb731FgaDgSFDhlTa/vzzz/Pcc88BMHXqVAoK\n",
       "CpgwYQIZGRn06tWL1atXY29vX3tVi4iINBJms5k1a9bwf//8P7Zs3gKApbUFHXv7021QMA7Ov88e\n",
       "4urgTXuvCP7fsAexstRQEZGmoMawbTKZLusks2bNYtasWdelIBERkabAZDLx1VdfER39Ejt37gLA\n",
       "2taSLv0D6dI/EDsHm4pjfd0Dub3vX8g5X/5wpIK2SNNxRVP/iYiISM1KSkr45JNPmD17NocOHQLA\n",
       "zsGGbgOD6HRjW6xtfw/Srk4ejOj9AN3a9cVoMBJ7XjNXiDQ1CtsiIiLXQWFhIYsXL2bu3LkkJSUB\n",
       "4NjKju6Dg+nQ0x9L699Xe7S3a8nNPe6mb6ebNLOISBOnsC0iInINcnJyePvtt5k/fz4pKSkAOLs7\n",
       "EDEkhHYRvlj8YbVHa0sbBnW/ncHd78DOpkV9lSwidUhhW0RE5CpcuHCB119/nddff53MzEwA3Hyd\n",
       "iBjajsDOXhj/sNqj0WCkd8dh3NLzHlrat6qvkkWkHihsi4iIXIEzZ84wf/583nnnHfLy8gDwCWpN\n",
       "+NBg/Nq7V5kLu0twb27r8wDurXzqo1wRqWcK2yIiIpchMTGRuXPn8sEHH1BcXAxAQAcvug0Jwieo\n",
       "6lzYQT5hjOw7lgCv0LouVUQaEIVtERGRGuzfv5+YmBiWLVuGyWTCYDDQPrwtXQb6497GucrxXq5+\n",
       "3NbnQcICIrTio4gobIuIiFRnx44dREdHs3LlSgAsLS3o2rc9N9zojYuHY5XjnR1cubX3GCLbD8Ro\n",
       "tKiyX0SaJ4VtERGRX5nNZjZs2EB0dDTr1q0DwMbWhsiBYQT1cKWlS9UZROxs7BkWOZp+XYZjbWlT\n",
       "Zb+ING8K2yIi0uyZTCa+/fZboqOj+fnnnwFwcHCg383heHexpYWjbZU2lhZWDOh6K0Mj7sTetmpP\n",
       "t4gIKGyLiEgzVlpaymeffUZMTAz79+8HwMXFhSEj++ASasbarurHpAEDPW4YxC297sOlpVtdlywi\n",
       "jYzCtoiINDtFRUV89NFHzJkzh8TERADcPdwYPnowLQNLMFuUVduuQ9twRvZ9EO/WbeuwWhFpzBS2\n",
       "RUSk2cjISueV1+bx9pvvknb+AgDOrR3oPiSY9pG+WFgWYq6mnb9HCCNvHEuIb6e6LVhEGj2FbRER\n",
       "aZJyC7I5nXqcU6nHOXrsAF9++g1bVu2mMK98jmxXr5ZERIUQ3MUb4x+WVP8jN2dvRvR5gK7BvTWN\n",
       "n4hcFYVtERFp1ExmExeyUjidepzTaeXh+lTqcbJyL5CfU8juDYns++k4JUXlQ0M8/FsRGdWOth08\n",
       "MBirD9COLZy5uec99AmLwsJCH5UicvX0HURERBqNktISzqWf4FTqcU6nHiv/b1oSRcUFlY7LTs9n\n",
       "1/oEDv6cTFmJCYA27dyIiArBJ7h1lV7qFraO+LoF4OsWQBv3IDoGRGJjbVdn1yUiTZfCtoiINGjZ\n",
       "eZmsiV1O/Ml9nMs4hclU/cOLAOnncti5Lp6jO09hMpWPvg7s5En40HZ4+rcCwNXJA9/WAfi4lX/5\n",
       "ugXg7FA1gIuIXA8K2yIi0qD9Z83rHEreVeMx509mErv2KIlxZ8EMBgOEhrdh+F0DCe8eie+vwdqn\n",
       "dVvsbOzrqHIREYVtERFp4M5eSL7ovtOJacSuiefE4fMAWFpZctuom5ky5Wl6dO+NpYVVXZUpIlIt\n",
       "hW0REWnQ/DyCycy9UPHaxsqOovPWrP/qZw7sPQKAvb09f//735k8eTLe3t71VaqISBUK2yIi0qB1\n",
       "De5DXOLPmExmEuPOsGf9cc6dKA/frVq14sknn+SJJ57A1dW1nisVEalKYVtERBq0dr5dOfLLaXas\n",
       "PkRmah4Abu5uTH16Kn/7299wdHSs5wpFRC5OYVtERBqk/Px83nvvPV5++WVOnjwJQEuXFnQfEsxj\n",
       "f3ucu4aMr+cKRUQuTWFbREQalKysLN58801effVVUlNTAQgI8ie4tysh3XywsDBy8EQsZvPDmq5P\n",
       "RBo8hW0REWkQUlNTWbBgAW+88QZZWVkAdO7SkcF39MLomkWJqbji2AvZKZzPPINHK5/6KldE5LIo\n",
       "bIuISL06efIk8+fP591336WgoHwlyG4Rnek+NAg7jzLKDKmUmeq5SBGRq6SwLSIi9SI+Pp45c+bw\n",
       "0UcfUVJSAkC3Xh0J7euBWxsHwARUHSZiMBiJivh/6tUWkUZBYVtEROpUXFwcMTExfPbZZ5hMJoxG\n",
       "I116h3LDjV64+ThdtJ2Tgyu9w4bSO2worRzd6rBiEZGrp7AtIiJ1Ytu2bURHR/Ptt98CYGFpQade\n",
       "AXQdFIizm0O1bQwY6NA2nD6dhtGhbTgWRou6LFlE5JopbIuISK0xm82sXbuW6OhofvzxRwCsrC0J\n",
       "6+1Ht0HBODjbVdvOyd6F3mFR9AobiktL9WKLSOOlsC0iItedyWRi5cqVREdHExsbC4CNnRWdbgyg\n",
       "64BA7BxsqrQxYOAG/2706XQTYQER6sUWkSZBYVtERK6b0tJSli1bRnR0NIcOHQLAzsGargOD6NQ3\n",
       "ABs7qyptWtq3+nUsdhQuLd3rumQRkVqlsC0iItessLCQDz74gJiYaE6cKF/t0cHZju5DgunQ0w8r\n",
       "68ofNwYMtPfvRt9OwwgLiFQvtog0WQrbIiJy1XJyclj05hvMmzePC2npADi72RM+NITQ8DZYWBor\n",
       "Hd+yRSt6hQ2ld8ehuLb0qI+SRUTqlMK2iIhcsfT0dP4Z8yLvvvMv8nLyAWjt40REVAhBnb0xGn+f\n",
       "H9uAgVD/rvTtOIyOAZFYWOijR0SaD33HExGRy5Z8IokZs6ay/JOvKC4qX4jGK8CFiKh2+N/gjsHw\n",
       "e8h2bOFcMRbb1Um92CLSPClsi4jIJW2L3cSsF59l/Q9bKCstXzvdr707EVHt8AlyrXRsqF8X+na8\n",
       "iU6BPdSLLSLNnr4LiojIRX3xw3+IiYlm10+HMJvMYICgLl5EDG2HexvniuMc7ZzoGTaUPh2jaO3k\n",
       "WY8Vi4g0LArbIiJSxS+//MKUaZPYtP4nAAxGA+0j2xA+JAQXT8eK40LbdKFPp5voFBiJpUXVaf1E\n",
       "RJo7hW0REQHKV3vcuHEj0dHRrFmzBgALSyMdevnRfVAILV1bAOBg50SvDkPo3TEKN2ev+ixZRKTB\n",
       "U9gWEWnmzGYz3333HdHR0Wzbtg0AR0dHho8agn1wAfYtbQGwsrTm/qgn6RzUU73YIiKXyXjpQ0RE\n",
       "pCkqKytj2bJldO3aldtuu41t27bh4uLCiy++SHJyMm+98S8cWtpVHF9SWoyXq5+CtojIFVDYFhFp\n",
       "ZoqKivj3v/9N+/btue+++4iLi8Pb25tXXnmF5ORknn32WVq1akUrx9a0a9O5Utsdh9bXU9UiIo2T\n",
       "wraISDORl5fHggULCAoKYvz48SQkJBAYGMi7777LsWPHmDRpEg4ODpXa9OgwqNLrXw5vxGQqq8uy\n",
       "RUQatUuG7U2bNjFy5Eh8fX0xGo18+OGHlfaPGzcOo9FY6atPnz61VrCIiFyZzMxMXnrpJdq2bctT\n",
       "Tz3F6dOn6dixI//5z384cuQI48ePx8bGptq2nYN6YWNlW/E6Oy+DIyfj6qp0EZFG75JhOy8vj86d\n",
       "O7NgwQLs7OwqrQ4GYDAYiIqK4ty5cxVf33//fa0VLCIilyclJYXp06fj5+fHzJkzSUtLo0ePHqxc\n",
       "uZK9e/cyZswYLC1rfk7exsqWrsGVO1B2HNRQEhGRy3XJ2UhuueUWbrnlFqC8F/vPzGYz1tbWuLu7\n",
       "X/fiRETkyp04cYJ58+bxr3/9i8LCQgCGDBnCjBkzGDRoUJVOk0vp0WEQP/9hrHZc4s8UFOVjZ9Pi\n",
       "utYtItIUXfOYbYPBwJYtW/Dw8CA0NJRHHnmE1NTU61GbiIhcgSNHjvDQQw8RFBTEwoULKSwsZOTI\n",
       "kWzbto21a9cyePDgKw7aAEE+YbRydKt4XVJWzJ74n65n6SIiTZbBbDabL/dgR0dHFi1axNixYyu2\n",
       "ffrpp9jb2xMQEMDx48eZOXMmZWVl7Ny5E2tr64rjsrKyKv4cHx9/ncoXEZEjR47wwQcfsG7dOsxm\n",
       "M0ajkaioKMaNG0dwcPAVn6+krJgLuWdIyzlDWu4Z0nJOk1+cU+kYj5b+3NTpwet1CSIiDVZISEjF\n",
       "n52cnK64/TUvanPPPfdU/DksLIzw8HD8/f357rvvGDVq1LWeXkRELmLPnj0sXryYrVu3AmBlZcWI\n",
       "ESMYO3Ysvr6+l3UOk9lEZt75ilCdlnuGrPw0zNTcD5NTmH7N9YuINAfXfQVJLy8vfH19SUhIuOgx\n",
       "ERER1/ttpQ7FxsYCuo9Nge5l42M2m/nvf/9LdHQ0mzdvBsDW1pZRo0bx8ssv4+PjU2Pb9JzzJJ+L\n",
       "J/ncUZJT4jl5PpGS0uIrrqNn2CD9vakF+jfZdOheNh1/HJ1xNa572E5NTeX06dN4eXld71OLiDRb\n",
       "JpOJFStWEB0dza5duwBwdnbmiSeeoH///jg7O1cJ2vmFuSSn/Bqsz8WTnBJPbsHVfWhYWFji6xaI\n",
       "v0cI7f26EhagACEicjkuGbbz8vIqxlibTCaSk5PZs2cPrq6uuLi4MGvWLEaPHo2npydJSUlMnz4d\n",
       "Dw8PDSEREbkOSkpKWLp0KbNnz+bw4cMAeHh4MHnyZP7+97/TsmVLYmNjKTOVknTuaKVgnZp55qrf\n",
       "172VD2092+HnEYK/Rwg+bm21TLuIyFW4ZNj+5ZdfGDx4MFA+88isWbOYNWsW48aN480332T//v0s\n",
       "WbKEzMxMvLy8GDx4MMuXL8fe3r7WixcRaaoKCgp4//33mTt3LidOnADA39+fqVOn8pdxfyG3KJPD\n",
       "p3eSvDOeg8f2kJF3DtM201W9l2MLZ/w92+HvEUJbz3a08QiihY3DpRuKiMglXTJsDxw4EJPp4t/A\n",
       "V61adV0LEhFpzrKzs3nrrbd45ZVXOH/+PADt2rVj7Pj7COsRwOm0Y7z44XgKivOv6vzWlja08Qim\n",
       "rWcIfh7lAbuVY+urmhJQREQu7bqP2RYRkSuXlpbGggULeOONN8jMzATAP8iHHsPa4xFiT4pxDym7\n",
       "91zROQ0GI16ufpWCtadrGyyMFrVxCSIiUg2FbRGRelJmKmPXvh3Mm/cyK5d/R1Fh+awg3kGuRAxt\n",
       "h197tyvqcXZxdMPPMwR/j3a09QzB1z0IGyvb2ipfREQug8K2iEgdMJvNZOSkkZxS/gBj7J7tfL1s\n",
       "NQd+TsJUVj6ntX8HDyKGhuAd6HrJ89lZt8DPs3yMdXGOgdYO3vTrM7CWr0JERK6UwraISC3IL8rl\n",
       "xLmE36feS4knJz+TtDNZ7FwbT/zu05jNgAGCu3oTMTQEN1/nas9lYWGJb+sA/D1D8Pv1IcbWzl4Y\n",
       "DUbg9/l8RUSk4VHYFhG5RiWlJZxJS6rotU5Oied8xulKx5xLSid2TTzHD5wDwGg00D6yDeFDgmnl\n",
       "4VjpWHdn74pea3+PELxbB2BlqWn3REQaI4VtEZErYDabSc08+3uwPneUU2nHKSsrrfbYk0dT2bk2\n",
       "nlPxaQBYWBkJ6+VP98HBOLZqgYOdE/6/Buvf5rRuYatp90REmgqFbRGRGuTkZ1b0ViefO8qJlATy\n",
       "i3JrbGM2mTm2/xyxa49y/kT5zCLWtpZ07RfMbXcPo2No14rhIK0cr+whSBERaVwUtkVEqrHj0AZW\n",
       "/fwpaVnnLruNqcxE/O7TxK6NJ/1cDgCOTg7cN/ZOJj45idDAjpp2T0SkmVHYFhH5k9TMs/xn9euY\n",
       "MV/W8aUlZZyIy+Dn1YdIS0kHwMfHm6efnsrDDz+sFXVFRJoxhW0RkT8pKMqrMWjbWrfA3yME95Zt\n",
       "+Hn9Pj5f/CnnzqUAEBwczLRp03jwwQextrauq5JFRKSBUtgWEfkTH7cAPFx8SUk/VbHN0sKKuwf9\n",
       "nQCvUCzMtix6YxHPvf4P0tPLe7I7d+7MjBkzGD16NBYWGioiIiLlFLZFRP7EwmjBXQMf4Y0vn6vY\n",
       "VlpWQkpKCl9++F/eeustcnPLH5Ls3bs3zzzzDMOHD9eDjiIiUoXCtohINdq16Uy3kL7sjv+JrAt5\n",
       "7F6fwJtTvqGs1ATAsGHDmDFjBv3791fIFhGRi1LYFhG5iPZuvZn3wiIOxZ7AbCofwx15Y1fefPVf\n",
       "RERE1HN1IiLSGChsi4j8SWxsLNHR0axYsQIAg9FAaIQvEUNDcPV0wsOvVT1XKCIijYXCtogI5as9\n",
       "btq0iejoaFavXg2AjY0N48aNo0VwLqVW5WO0zZj5/Md3eequGIwGY32WLCIijYA+KUSkWTObzXz3\n",
       "3XfceOONDBw4kNWrV+Pg4MDTTz/N8ePHefvtt3l49ORKbZLOHuGXQz/WT8EiItKoqGdbRJqlsrIy\n",
       "li9fTkxMDHv37gXAxcWFJ598kieeeALbFtacvZDMlri9nL1wAisLa0rKiivaf73lQ7oE98bW2q6+\n",
       "LkFERBoBhW0RaVaKi4tZsmQJc+bMIT4+HgB3DzfuGfv/6Dm4M1kF53n1yylk5aXXeJ6cgiyOnTlI\n",
       "h7bhdVG2iIg0UgrbItIsZGVnsuCNV1i08G3On0sFoJWbI10HBdI+sg1YnePnI+eu6JzODq1ro1QR\n",
       "EWlCFLZFpEkpKyslNessZy+c4GzaCRJPHOHr5avY8sNuCvPKh4G4eDoSMTSEkG4+GC2u7NEVK0tr\n",
       "vFz9Gdz9drxb+9fGJYiISBOisC0ijZLJbOJCVgrn0k9yNi25PFxfOEFK5mnKykrJzyliz8ZE9m05\n",
       "TnFhKQDufs5ERrUjIMwTg7HmhWgsjJa4t/LGy9UfL1e/ii/Xlu4YjVqOXURELo/Ctog0aGazmczc\n",
       "CxVh+txv/00/SXFpUZXjczIK2L0hgQPbkiktKQPAN6Q1EVHt8A1pXWW1R4PBiJuTJ16ufnhWhGp/\n",
       "3J29sLDQt0gREbk2+iQRkQYlOy+TvQlbOZ2WVBGuC4rzL9ku43wuu9bFczj2JKay8tUeA8I8iYgK\n",
       "wbOtCwAujm6/91S39sPTxQ8PFx+sLW1q9ZpERKT5UtgWkQbBbDaz7cAaVmxeTFFxwWW3Sz2dReya\n",
       "oyTsPQNmMBggrEcAI+8dRkR4D7xc2uDV2h9Plzaapk9EROqcwraI1Lv07FSWrVvE4RN7LrvN2ePp\n",
       "7F5/jMR9pwGwtLTkzrvvYMaMZ+gc1rW2ShUREbkiCtsiUm/Ke7PXsmLz+zX2ZttY2eLlWt47fSr+\n",
       "AsuXfMMv23cC0KJFC/72t78xefJkfH1966p0ERGRy6KwLSL1IiMnlU/Wvcnh5N1V9llb2TK42+34\n",
       "e4bg5eqHk70rK1euJHp6NDt3lodsJycnnnjiCZ588knc3NzqunwREZHLorAtInXKbDaz/cBaVmxe\n",
       "TGE1Dz4G+3bk/qFP4OrkQUlJCcuWLSMmJoZDhw4B4O7uzqRJk3jsscdo2bJlXZcvIiJyRRS2RaTO\n",
       "ZOSksmzdWxxK3lVln7WlDSNv/As3dr6Z4qJi3nrrLebOnUtSUhIAfn5+TJ06lYceegg7Oz3oKCIi\n",
       "jYPCtojUukv2ZvuEMSbqCWyM9syfN5/58+eTkpICQGhoKNOnT2fMmDFYWVnVdekiIiLXRGFbRGpV\n",
       "Rk4ay9a9WUNv9lhu8OnBwlfeYOHChWRkZADQrVs3ZsyYwahRo7Cw0IqNIiLSOClsi0itMJvN/Hxw\n",
       "PSs2vVftojRBPmEMDhvNB//+mOHv3ENeXh4A/fr1Y8aMGdx0001VVnsUERFpbBS2ReS6y8y9wLJ1\n",
       "b3IwaWeVfdaWNoT7R/HjN78w9f5eFBcXA3DLLbcwffp0+vXrV9flioiI1BqFbRG5bi7Vm92izI2E\n",
       "n9J59Yu/YTKZMBgMjB49mhkzZtCtW7d6qFhERKR2KWyLyHVRU2922qlcTsTm8tOPK4Hy1R7Hjh3L\n",
       "P/7xD9q3b1/XpYqIiNQZhW0RuSbFpUVs27+G77ctrdSbbTabOZ2QxoFNZzi6LwkAW1tbHn74YaZM\n",
       "mYK/v389VSwiIlJ3FLZF5KoUlxaxdd9q1u78kuy8jIrtZpOZpIMp7Fwbz9mkdAAcHR2ZMGECTz31\n",
       "FB4eHvVVsoiISJ1T2BaRK1JcWsRP+/7LutgVZOf/HrJNZSYS9pwhdm08F85mA+Dq6sqkSZOYMGEC\n",
       "zs7O9VWyiIhIvVHYFpHLUlxSxJZ9q1i3cwU5+ZkV28tKyzj0y0l2rUsgK618+j4fHx+mTJnC+PHj\n",
       "sbe3r6+SRURE6p3CtojUqKikkJ/2rWJd7ApyCrIqtpcUlbJ/WzK7NySQl1UIQEBgADOmz+DBBx/E\n",
       "xsamvkoWERFpMBS2RaRaRSWFbIn7gXU7vyL3DyG7ML+YfVuOs2fjMQrzyufIDgvrwMyZzzJ69Ggs\n",
       "LfVtRURE5Df6VBSRSoqKC9gc9wPrd62sFLLzcwrZ/WMi+7YkUVJUCkBkj0iee/Y5br31Vq32KCIi\n",
       "Ug3jpQ74/+3de1TVVf7/8SeHO6IoKiKX5BJ4RTPxhqPjNXXya2PNt+wyaWZ2o1RyFLG++Z2p46U0\n",
       "Ky11mgrnO1o5v8nmp84vNctyNO+Kd1FQAQEvIAKCwDmf3x8U305eMOFw4Ph6rMVasM/en/M+a6+9\n",
       "evlpn8/+9ttvGTlyJCEhIZhMJpKTk6/qM3PmTIKDg/Hx8WHAgAEcOnTILsWKiP1cKSth/c5/MPPj\n",
       "p/nnv5dVBe1LeZf55u8pfPzH9ez+6jjlVyoYOGgAX3/9Ndu+38aIESMUtEVERK6j2jvbxcXFdO7c\n",
       "mTFjxvD4449f9R/VOXPmMH/+fJKTk4mOjuaPf/wjQ4YM4ejRo/j6+tqtcBGpOYu1gpyCU6Rt3Mme\n",
       "1H9TXFpY9VpeTiG7vkrl2K5MrFYDgP8YOYJXXv4vunfv7qiSRUREGpRqw/bw4cMZPnw4AGPHjrV5\n",
       "zTAMFixYwPTp0xk1ahQAycnJBAQEsHz5ciZMmFD7FYtIjRSXFnLo5C72p23nYNpOyi1lNq+fzbjI\n",
       "zg3HOJGSDQaYXE2MfvhBXp7xCh07dnRQ1SIiIg1TjfZsp6enk5ubyz333FPV5uXlRb9+/diyZYvC\n",
       "tkg9caEgl5S0bRxI28GJrINYDetVfbJOnGfn+lROHzkLgJu7G2PGjiEpMYmIiIi6LllERMQp1Chs\n",
       "5+TkAFx1IlxAQABnzpypyaVFpAashpWM3BPsT9vOgbTtnLlw6pr9DMPg1OGz7NxwjOy0ytMevbw9\n",
       "eaPjp3gAACAASURBVPrpp5k2NZHWrVvXZdkiIiJOx25PI7nRF6Z27txpr7eVOqR5rF8q91+fJOPC\n",
       "MTLyUykpK7xuX6vV4ETKGXauT+V8VuUXIX19GzF69Ggeemg0TZs2JSsri6ysrLoqX2qB1qTz0Fw6\n",
       "D81lwxcVFVWj8TUK24GBgQDk5uYSEhJS1Z6bm1v1mojYz5XyEjLzU8nIO8aZ/DQqrGU37G+psHJ0\n",
       "VwZ7NqaTl1sZsps3b85jjz3GqFGjdNqjiIhILatR2A4PDycwMJB169bRrVs3AEpLS9m8eTNvvvnm\n",
       "dcfFxsbW5G3FwX78V7rm0THOF+Sw/8R29qdtI+3M4Wvuv/65inIrZw+X8t2aPZzNPQ9AWFgYo0eP\n",
       "ZsSIEfTp08feZYsdaU06D82l89BcOo+CgoLqO93ATT36LzU1FQCr1cqpU6fYu3cvzZs3JzQ0lEmT\n",
       "JmE2m2nXrh1RUVG89tprNG7cmEceeaRGhYlIpcr918fZn7ad/Wnbyb5w+qbGebh70aZFew5vzWTl\n",
       "31Zx/lxlyO7QoQPTp09n9OjR7N27156li4iI3PaqDds7duxg4MCBQOU+7FdffZVXX32VsWPH8uGH\n",
       "HzJ16lRKSkp4/vnnyc/Pp1evXqxbt07/O1qkBsoryjiWkVL5Bcf0HVwqzr+pcU0aNSMmvAdBfney\n",
       "9h8beSXprap/kcfGxjJjxgxGjhyJyVTteVYiIiJSC6oN2/3798dqvfH/pv4xgIvIrSsuucTBk7vY\n",
       "f2Ibh0/vpay89KbGtW5+BzERPYmJ6A5lnrw1/y2WLn2JkpISoHINz5gxg0GDBumkRxERkTpmt6eR\n",
       "iEj1zl3MrtoeknbmMMZN7L82uZiICO5ATEQPYiJ60MIvkNTUVP708hyWLVtGeXk5ACNGjGD69OnE\n",
       "xcXZ+2OIiIjIdShsi9Qhq2HlVE4qB34I2Dl5GTc1ztPdi/Zt7iYmsgcdwrrRyKsxAPv27SN+1iRW\n",
       "rlyJ1WrFZDIxevRoEhMT6dKliz0/ioiIiNwEhW2ROpJ25gh/W/c25wqyb6q/XyN/Ov1w9zoqJAZ3\n",
       "N/eq17Zu3YrZbGb16tUAuLu7M27cOKZOnVrj54GKiIhI7VHYFqkDp3KO8f6qmVypZh92UIswYiK6\n",
       "ExPRk9CASJs91oZhsGHDBsxmM9988w0A3t7eTJgwgZdeeonQ0FB7fgQRERG5BQrbInaWde4k76/6\n",
       "4zWDtsnFxJ3BHYmJ7Emn8O4092t1VR+r1coXX3yB2Wyuem6rn58f8fHxTJw4kZYtW9r9M4iIiMit\n",
       "UdgWsaOz+Vm8t2oml68U2bR3juxJlzvj6BjWDR8v32uOraio4JNPPmHWrFkcOnQIgJYtW5KQkMCz\n",
       "zz6Ln5+f3esXERGRmlHYFrGTvEvnWPSPVym8fNGmfWiP/+Te3o9ed1xpaSkff/wxc+bM4eTJkwCE\n",
       "hobyhz/8gSeffBIfHx97li0iIiK1SGFbxA4uFeez6B//RX7ReZv2X981gt/0uvbpqoWFhSxZsoR5\n",
       "8+aRk5MDQHR0NImJiTz66KN4eHjYvW4RERGpXQrbIrWsuLSQ9z6fedVTR3p2GMSofuOuOljmwoUL\n",
       "vPvuu7zzzjvk51eeFHnXXXeRlJTE/fffj6ura53VLiIiIrVLYVukFpWWlfD+qj9y5sIpm/a7ouJ4\n",
       "eNBzmFz+95j07Oxs5s2bx+LFiykuLgagT58+zJgxg2HDhum0RxERESegsC1SS8oqrrD0n69xOjfV\n",
       "pr1DWDceHzoZk6nyDnVaWhpvvPEGH374IWVlZQAMHTqUGTNm0Ldv3zqvW0REROxHYVukFlRYyvlw\n",
       "zVyOZx20ab8zuCPj7p2Km6s7Bw8eZPbs2axYsQKLxYKLiwsPPPAA06dPp1u3bg6qXEREROxJYVuk\n",
       "hqxWC8u+fItDJ3fZtLdpFcWEkS+zb08KZrOZVatWAeDq6sqYMWOYNm0a7du3d0TJIiIiUkcUtkVq\n",
       "wGpYWfHVe+xN3WLT3tr/Djr4D2DkiPtYv349AJ6enowfP54pU6YQFhbmgGpFRESkrilsi9wiwzD4\n",
       "/NsP2XboK5u2/FPlbFvxPUnb3gXA19eX5557jsmTJxMYGOiockVERMQBFLZFbtH/2/Ypm/auBsBq\n",
       "NTi+N4u9G9PJzcwDwN/fn0mTJhEfH0+zZs0cWaqIiIg4iMK2yC3ILzzH/9v2KZYKC0d2ZLLrq1QK\n",
       "zlc+vi8oKIgpU6bw1FNP4et77aPYRURE5PagsC1yCzJzTrJn03F2bzxOcUEpAG3a3MGMGS/z+OOP\n",
       "4+np6eAKRUREpD5Q2Bb5BS5evMiiRYuYN38e+XmVpz02b92YEQ8O4oM3V+LmpiUlIiIi/0vJQOQm\n",
       "5ObmsmDBAhYtWkRhYSEAre5oSuw90YR3CKR9+F0K2iIiInIVpQORGzh9+jRvvPEGH3zwAaWlldtF\n",
       "Bg0axOix97Pv/JdVR6q7mbSURERE5GomRxcgUh8dPXqUcePGERkZycKFCyktLeW+++7j+++/Z8OG\n",
       "DXTreVdV0AZw/eEodhEREZGf0u04kZ/Ys2cPs2bN4u9//zuGYWAymXjkkUdITEwkJiamqp/VWmEz\n",
       "ztXVva5LFRERkQZAYVsE2Lx5M2azmX/9618AeHh4MHbsWKZOnUpkZORV/SsstmHbzVVLSURERK6m\n",
       "hCC3LcMwWLduHa+//jrfffcdAD4+PjzzzDMkJCQQHBx83bEWq8Xmb5O2kYiIiMg1KGzLbcdqtfL5\n",
       "559jNpvZvXs3AE2bNuWFF17gxRdfpEWLFjccb7FaKLx80aZNX5AUERGRa1FCkNtGeXk5y5cvZ/bs\n",
       "2Rw5cgSAVq1akZCQwDPPPEOTJk1s+ldYyjl3MYecvAxy8jLIzcsg50IGuRezsFh+vmdbS0lERESu\n",
       "poQgTq+kpISPPvqIuXPncurUKQDatGnD1KlTeeKJJ3Bzd+PcxSyOH0sh50JGVbg+e/EM1p9tF7ke\n",
       "dzedGCkiIiJXU9gWp3Xp0iUWL17M/Pnzyc3NBeDOqEgeG/cgXeKiOX8pm3mfvcT5ghyshvWW38fF\n",
       "xUTHsG61VbaIiIg4EYVtcTpZZzKY++YcPvpLMoWXigBo3aYFXQeGExHTmgumQ2zcc6hG79HYpymB\n",
       "/qEE+odyd3QfIoM71EbpIiIi4mQUtqXBKrlymdz8zKqtH0eOHeCLT79k56bDVJRVbv8IimxO7OBo\n",
       "7mjX0uYQmpvl18ifwOahVcG68ieERt5Nqh8sIiIitz2Fban3LpcWVe2j/ume6otFFwC4eK6I3RuP\n",
       "c3j7aawWA4A2HVoROziKoIjmN/UezRq3rArSgf6hBDYPpZV/CD6evnb7XCIiIuL8FLal3igquXRV\n",
       "oM7Jy+BScf41+58/c4ldG46RuicLwwBc4M67gogdHE3LEL9rjmnepFVVmP7xTnUr/xC8PLzt+MlE\n",
       "RETkdqWwLXXKMAwKLxf85E716R9+z6SopOCmrpFzMo+dG1JJP5ADgMnkQvseodw9KIpmAb64uJho\n",
       "4Rdoc5c60D+UgGbBeLp72fPjiYiIiNhQ2Ba7MAyDguK8q+5S5+Rlcrm08Jaul5l6np3rj5GZeh4A\n",
       "N3dXuvfvyG9HD6dj+xgC/e/4IVQH4e7mUdsfSUREROQXU9iWWpeauZ//WfcO+YXnanwtw2pw6tA5\n",
       "9mxMIzOt8vF9vr6NGPfUWKb9YTpBra9/pLqIiIiIoylsS63Kycvgz/93FqVll3/xWDdXd1o1CybQ\n",
       "P5SWfkHs3XqUZR98wtEjRwFo0aIFkydP5rnnnqNp06a1XbqIiIhIrVPYllpzubTopoK2h5snrX7c\n",
       "T/2TPdXNmwRQVlZOcnIy8XNmkJ6eDkBISAh/+MMfGD9+PD4+PnXxUURERERqhcK21AqL1cJH/3qD\n",
       "cxfP2LSHBEQQ3DzM5ukfzZq0xORisulXVFTEggVv8+abb5KdnQ1AVFQUiYmJPPbYY3h4aA+2iIiI\n",
       "NDwK21IrvticzNHT+2zaenYYxCOD4294mExeXh4LFy7k7bffJi8vD4AuXbqQlJTEAw88gKurq13r\n",
       "FhEREbEnhW2psW2HvuKbPf+0aQtr3ZYHBzxz3aCdnZ3NW2+9xfvvv09RUeWR6nFxccyYMYPhw4ff\n",
       "0mmPIiIiIvWNwrbUSHr2ET7Z+L5NW1Pf5oy/NxF3N/er+6en88Ybb/Dhhx9y5coVAIYOHUpSUhJ9\n",
       "+/ZVyBYRERGnorAttyy/8Dx/WT0Hi6Wiqs3d1YPxI6bTpFEzm76HDh1i9uzZLF++HIvFAsD999/P\n",
       "9OnTiY2NrdO6RUREROqKwrbckgpLOR+snsWly7ZHqT8yJJ47Wt1Z9ffOnTsxm818/vnnALi6uvL4\n",
       "448zbdo0OnToUKc1i4iIiNQ1U/VdbmzmzJmYTCabn6CgoNqoTeopwzDYcnw1GWdP2LQPiX2Abm37\n",
       "YRgGmzZtYujQoXTv3p3PP/8cT09PnnvuOY4fP05ycrKCtoiIiNwWauXOdrt27fjmm2+q/tYTJJzb\n",
       "waytnDx/0KatU3h3ftP7EdasWYPZbGbLli0A+Pr68uyzz5KQkEBgYKAjyhURERFxmFoJ266urgQE\n",
       "BNTGpaSeO5C2g92nNtq0BTQNxvtSKN3u7sa+fZWP//P392fixInEx8fj7+/viFJFREREHK5WwnZa\n",
       "WhrBwcF4enrSs2dPzGYz4eHhtXFpqUeyL2SQ/OX8qr8tFVbS951jzZb9pJ14D4DWrVszZcoUJkyY\n",
       "gK+vr6NKFREREakXahy2e/XqRXJyMu3atSM3N5fXXnuNuLg4Dh48qDuaTsRqtfDXL9/iSlkJ5WUV\n",
       "HNp6it1fn6DoYgkAERERTJs2jccffxwvLy8HVysiIiJSP7gYhmHU5gUvX75MeHg4iYmJTJ48uaq9\n",
       "oKCg6vfU1NTafEupA6k5e/hm/ypS/p3Ovk0nKCkqAypD9tixYxkyZAhubnq4jYiIiDiXqKioqt/9\n",
       "/Px+8fhaT0c+Pj507NiR48eP1/alxUFyz+bw9rsL2PPtMcpKK5+pHRrRmhefTaBfv36YTDV+qI2I\n",
       "iIiIU6r1sF1aWsrhw4cZOHDgdfvoEJOGISMjgzfffJPFSxZTdqXyTnZIVAt63NOeF0b/F/1/NdjB\n",
       "FUpN7dy5E9CabOg0j85Dc+k8NJfO46e7M25FjcP2lClTGDlyJKGhoZw9e5Y//elPlJSUMGbMmJpe\n",
       "Whzk2LFjzJkzh2XLllFRUXknO7xTILGDowgM86dzaF98vZo6uEoRERGR+q/GYTsrK4uHH36Y8+fP\n",
       "07JlS3r37s33339PaGhobdQndWjv3r3MmjWLlStXYhgGJpOJHr/uQmTPZrQIqtyj5NfIn47BvR1c\n",
       "qYiIiEjDUOOwvWLFitqoQxxoy5YtvP7666xduxYAd3d3xo4dy4O//y3/3L3Upu9/9Pk9pmIPR5Qp\n",
       "IiIi0uDom223KcMwWLduHf3796dPnz6sXbsWHx8fJk+eTFpaGosXv8/ujA02Y+5oFUVsu187qGIR\n",
       "ERGRhkfParvNWK1WVq1ahdlsZteuXUDlY2xeeOEFXnzxRVq2bAnA1oMbyDyXZjP2/n5PYnLRv89E\n",
       "REREbpbC9m2ivLycFStWMHv2bA4fPgxAQEAACQkJPPvsszRp0qSqb2lZCau3/I/N+Luj+xIR1K5O\n",
       "axYRERFp6BS2nVxpaSkfffQRc+fO5eTJkwDccccdTJ06lXHjxuHt7W3T/1LxRdZs/RuFly9Wtbm7\n",
       "ejCyz+/rsmwRERERp6Cw7aQKCwtZvHgx8+bNIzc3F4C2bduSmJjIo48+iru7e1XfcxezSTmxjf0n\n",
       "tpGefQQD20NFB3a7D/8mAXVav4iIiIgzUNh2MhcuXOCdd97h3XffJT8/H4CuXbuSlJTEqFGjcHV1\n",
       "xTAMTuceZ3/aNlJObCP7wunrXq9Jo2YM7nZ/XZUvIiIi4lQUtp3EmTNnmDdvHkuWLKG4uBiAX/3q\n",
       "V8yYMYOhQ4ditVo4nnWgKmBfLLpQ7TVdXEw8OOBpPD28q+0rIiIiIldT2G7gTpw4wdy5c/n4448p\n",
       "K6s8Un3YsGEkJSXRo1d3Dp/czV/XLeBg+k5KrhTf1DVDAiLoHNGTu6P7EtAsyJ7li4iIiDg1he0G\n",
       "6sCBA8yePZsVK1ZgtVpxcXHhd7/7HS9OisetaTn7T3zLqiXvUm4pq/ZaJhcTdwZ3JCayJzERPfFv\n",
       "0rIOPoGIiIiI81PYbmC2b9+O2Wzmiy++AMDNzY0HR/8nQ+/vS4E1i7/vehvDsFZ7HXc3D9q3uZvO\n",
       "kT3pGB5LI6/G9i5dRERE5LajsN0AGIbB119/jdls5quvvgLA08uTISN+TYc+QZS45LM7c32112nk\n",
       "1ZhOET3oHNmTtqFd8HD3tHfpIiIiIrc1he16zGq1snr1asxmM9u2bQPA28eLu/u3pV3vQHwae1FC\n",
       "/g2v4d8kgM6Rvegc2ZPw1u1wNbnWRekiIiIigsJ2vVRRUcHKlSsxm80cOHAAAO9GnnT5dTgxvwrH\n",
       "y8fjhuODW4bTOaInnSN7EtQiDBcXl7ooW0RERER+RmG7Hrly5QpLP1jC3DlzyMw4A0AjPy/uHnAn\n",
       "HXu3wd3z2tPl4mIiMrgDnSN6EhPZg+ZNWtVl2SIiIiJyHQrb9cDpM2nMnvc6yz/6jIL8IgD8WjSi\n",
       "26Ao2nUPwdXt6q0flV9w7EpMROUXHH29m9R12SIiIiJSDYVtBzAMg6zz6fx7z1f85c8f8d2/dlFa\n",
       "XPmIvuatmxA7JIo7uwRhcjXZjPPxakyn8Fg6R/ak3R1d9QVHERERkXpOYbuO5F06R2pmCscy9rP7\n",
       "4Ha+XbuD/ZtPUn6lAoBWbZrRfUg0YR1a4WL63z3W/o1bEhNZuf86IqiDvuAoIiIi0oAobNtJ4eUC\n",
       "UjP3cywjhdSM/ZwryOZS3mV2bzzOoW2nsJRXPgs7JLoF3YdEE3xni6ovMga1CKPzDwE7uEW4vuAo\n",
       "IiIi0kApbNeSkiuXOZF1kGMZKRzL3M+Z8yerXsvLLWTXhlSO7crEajUAiIgJpNvgaALbNMPFxURE\n",
       "UPuqLzi28At00KcQERERkdqksH2LyiqukH7myA93r/dzOjcV689ObjybcZGdG45xIiUbDHBxgbbd\n",
       "Qug2OIrA0Ba0De1C58hedAyPpbGPn4M+iYiIiIjYi8L2TbJYLZzOTeVYRuXWkPTsI1RYyq/ZN+vE\n",
       "BXauP8bpI2cBMLma6NirDfeNHkrv2H5Eh3amTato3N3c6/IjiIiIiEgdU9i+DqthJfv8qapwffzM\n",
       "Qa6UlVy3v2EYnDp8lp0bjpGdlgeAh6c7w0cNZPLkBHp17Yenu1ddlS8iIiIi9YDC9g8Mw+Dcxewf\n",
       "9lynkJp5gOKSS9WOs1oNTqScYef6VM5nFQDQxK8xz8c/z0uTp9C8eXN7ly4iIiIi9dRtHbbzC8/b\n",
       "PDEkv+j8TY+1VFjJPFjAjvVHyM6s3C4SGBjISy+9xNNPP03jxo3tVbaIiIiINBC3VdguKrlEauaB\n",
       "H8J1CmcvnvlF4xt7+9EmoD2Ht2Xwj/9ZRVZm5fiwsDCmTZvG2LFj8fLSVhERERERqeTUYbu0rMTm\n",
       "cXxZ59J/0XgvDx/uDOlE29DOtGrShs8/Xc0rr7zF2bOVd7I7dOjA9OnTeeihh3B315cdRURERMSW\n",
       "04btr3atYs3Wv133iSHX4u7mQURQe6JDOhMdGkNIQCR5F/J4++23WbjwMQoKKvdkx8bGkpSUxH33\n",
       "3YfJZKrmqiIiIiJyu3LKsH069zhfbP642n4mkythraKJCo0hOrQzYYFtqx7Hl5GRwUsJL7F06VJK\n",
       "SiqfQtK/f3+SkpIYPHiwTnUUERERkWo5Xdg2DIN/bk6+5msuuBDcMpzoH8J1RFAHvDy8bfqkpqYy\n",
       "Z84cli1bRnl55V3xe++9l6SkJOLi4uxev4iIiIg4D6cL24dP7eFY5n6btl4dB9Ohzd1EhXSikXeT\n",
       "a45LSUlh1qxZfPbZZ1itVkwmE6NHjyYxMZEuXbrURekiIiIi4mScKmxbrZar7mq3vaMLjwyOv+6Y\n",
       "rVu3YjabWb16NQDu7u6MGzeOqVOnEhUVZdd6RURERMS5OVXY3nFkE2cunLJpG9lnzFX9DMNgw4YN\n",
       "mM1mvvnmGwC8vb15+umneemllwgJCamLckVERETEyTlN2C6vKGPt1uU2bbFtf01oQETV31arlS++\n",
       "+AKz2czOnTsB8PPzIz4+nokTJ9KyZcs6rVlEREREnJvThO1v962xOQHS1dWNe+MeAaCiooJPPvmE\n",
       "WbNmcejQIQBatmxJQkICzz77LH5+fg6pWUREREScm1OE7eLSQtbt+LtNW9/Ov6GRhx+LFy9m7ty5\n",
       "pKdXHmgTGhrK1KlTGTduHD4+Po4oV0RERERuE04Rttfv+D+UXCmu+ttkdefIv8/w/EPh5OTkABAd\n",
       "HU1iYiKPPvooHh4ejipVRERERG4jDT5s5106y7f71gBQWlzGvm/TOLw1i8JLlXe677rrLpKSkrj/\n",
       "/vtxdXV1ZKkiIiIicptp8GH7X9s+pSCvkD3fnODAv9MpL7MA0KdPH2bMmMGwYcN02qOIiIiIOESD\n",
       "DttpaWl89eku/rpsAxXllSG7Z59Y3pg1n759+zq4OhERERG53TXIsH3w4EFmz57NihUrsFgsuLi4\n",
       "8Jt7h9H3P7oy9ak/YTJpu4iIiIiIOF6DCts7duzAbDazatUqAFxdXRkzZgzTpk2jffv2Dq5ORERE\n",
       "RMRWvQ/bhmGwadMmzGYz69evB8DT05Px48czZcoUwsLCHFugiIiIiMh11NuwbRgGa9aswWw2s3Xr\n",
       "VgAaN27Mc889x6RJkwgMDHRwhSIiIiIiN1bvwrbFYmHlypXMmjWLlJQUAPz9/Zk0aRLx8fE0a9bM\n",
       "wRWKiIiIiNwcU21d6L333iM8PBxvb29iY2PZvHnzLxp/5coVPvjgA9q1a8fDDz9MSkoKQUFBzJ8/\n",
       "n1OnTvHKK68oaIuIiIhIg1IrYfvTTz9l0qRJvPzyy+zdu5e4uDiGDx9ORkZGtWOLi4tZsGABkZGR\n",
       "PPXUUxw/fpyIiAiWLl1KWloakydPxtfXtzbKFBERERGpU7UStufPn88TTzzBk08+Sdu2bXnnnXdo\n",
       "3bo177///nXHXLx4kddff52wsDAmT55MVlYWnTp1Yvny5Rw9epSnnnoKT0/P2ihPRERERMQharxn\n",
       "u6ysjN27dzN16lSb9nvuuYctW7Zcc8z06dNZtGgRhYWFAPTo0YMZM2YwYsQITKZa29kiIiIiIuJQ\n",
       "NQ7b58+fx2Kx0KpVK5v2gIAAcnJyrjlm9uzZAAwaNIikpCQGDBigI9VFRERExOk45GkkFy9etPn7\n",
       "0qVLjihDblFUVBQABQUFDq5Eakpz6Rw0j85Dc+k8NJfyoxrv2WjRogWurq7k5ubatOfm5tK6deua\n",
       "Xl5EREREpMGqcdj28PCgW7durFu3zqZ9/fr1xMXF1fTyIiIiIiINVq1sI0lISOD3v/89PXr0IC4u\n",
       "jsWLF5OTk8MzzzxT1cfPz6823kpEREREpMGolbD94IMPcuHCBV577TWys7OJiYlh7dq1hIaG1sbl\n",
       "RUREREQaJBfDMAxHFyEiIiIi4ozq5KHWNT3KXeqHmTNnYjKZbH6CgoIcXZZU49tvv2XkyJGEhIRg\n",
       "MplITk6+qs/MmTMJDg7Gx8eHAQMGcOjQIQdUKtWpbi7Hjh171RrVd2fqn1mzZtG9e3f8/PwICAhg\n",
       "5MiRHDx48Kp+Wpf1383MpdZlw7Bo0SK6dOmCn58ffn5+xMXFsXbtWps+t7om7R62a3KUu9Q/7dq1\n",
       "Iycnp+pn//79ji5JqlFcXEznzp15++238fb2vuqZ9nPmzGH+/PksXLiQHTt2EBAQwJAhQygqKnJQ\n",
       "xXI91c2li4sLQ4YMsVmjP/+PhTjepk2biI+PZ+vWrWzcuBE3NzcGDx5Mfn5+VR+ty4bhZuZS67Jh\n",
       "CA0NZe7cuezZs4ddu3YxcOBAfvvb37Jv3z6ghmvSsLMePXoYEyZMsGmLiooypk+fbu+3llr26quv\n",
       "Gp06dXJ0GVIDvr6+RnJyctXfVqvVCAwMNMxmc1VbSUmJ0bhxY2PJkiWOKFFu0s/n0jAMY8yYMcaI\n",
       "ESMcVJHcqqKiIsPV1dVYvXq1YRhalw3Zz+fSMLQuGzJ/f39j6dKlNV6Tdr2z/eNR7vfcc49N+42O\n",
       "cpf6LS0tjeDgYCIiInj44YdJT093dElSA+np6eTm5tqsUS8vL/r166c12gC5uLiwefNmWrVqRdu2\n",
       "bZkwYQLnzp1zdFlSjUuXLmG1WmnWrBmgddmQ/XwuQeuyIbJYLHzyySeUlpbSr1+/Gq9Ju4btWznK\n",
       "XeqvXr16kZyczJdffsmf//xncnJyiIuLIy8vz9GlyS36cR1qjTqHYcOG8de//pWNGzcyb948tm/f\n",
       "zsCBAykrK3N0aXIDEydOpGvXrvTu3RvQumzIfj6XoHXZkOzfvx9fX1+8vLyYMGECn332GW3btq3x\n",
       "mnTIce3SMA0bNqzq906dOtG7d2/Cw8NJTk5m8uTJDqxM7OHn+4Gl/nvooYeqfu/YsSPdunWjTZs2\n",
       "rFmzhlGjRjmwMrmehIQEtmzZwubNm29qzWld1l/Xm0uty4ajXbt2pKSkUFBQwMqVKxk9ejRff/31\n",
       "DcfczJq0651tHeXu3Hx8fOjYsSPHjx93dClyiwIDAwGuuUZ/fE0artatWxMSEqI1Wk9NnjyZTz/9\n",
       "lI0bNxIWFlbVrnXZ8FxvLq9F67L+cnd3JyIigq5du2I2m+nVqxeLFi2qyqy3uibtGrZ1lLtzKy0t\n",
       "5fDhw/qHUwMWHh5OYGCgzRotLS1l8+bNWqNO4Ny5c2RlZWmN1kMTJ06sCmfR0dE2r2ldNiw3mstr\n",
       "0bpsOCwWC1artcZr0nXmzJkz7VgnTZo04dVXXyUoKAhvb29ee+01Nm/ezEcffaQj3BuYKVOm4OXl\n",
       "hdVq5dixY8THx5OWlsaSJUs0l/VYcXExhw4dIicnh7/85S/ExMTg5+dHeXk5fn5+WCwWZs+eo+TA\n",
       "TgAAAZtJREFUTdu2bbFYLCQkJJCbm8vSpUvx8PBwdPnyEzeaSzc3N5KSkmjSpAkVFRXs3buX8ePH\n",
       "Y7VaWbhwoeayHnn++edZtmwZK1euJCQkhKKiIoqKinBxccHDwwMXFxetywaiurksLi7WumwgEhMT\n",
       "qzJORkYGCxYsYPny5cydO5fIyMiarUl7PjLlR++9954RFhZmeHp6GrGxscZ3331XF28rtWz06NFG\n",
       "UFCQ4eHhYQQHBxu/+93vjMOHDzu6LKnG119/bbi4uBguLi6GyWSq+v2JJ56o6jNz5kyjdevWhpeX\n",
       "l9G/f3/j4MGDDqxYrudGc1lSUmIMHTrUCAgIMDw8PIw2bdoYTzzxhJGZmenosuVnfj5/P/7893//\n",
       "t00/rcv6r7q51LpsOMaOHWu0adPG8PT0NAICAowhQ4YY69ats+lzq2tSx7WLiIiIiNhJnRzXLiIi\n",
       "IiJyO1LYFhERERGxE4VtERERERE7UdgWEREREbEThW0RERERETtR2BYRERERsROFbRERERERO1HY\n",
       "FhERERGxE4VtERERERE7+f9vZwKceRPTzAAAAABJRU5ErkJggg==\n"
      ],
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7febf5baa630>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "import math\n",
    "from filterpy.kalman import KalmanFilter\n",
    "import numpy as np\n",
    "import book_plots as bp\n",
    "\n",
    "pos_a = (100, -20)\n",
    "pos_b = (-100, -20)\n",
    "\n",
    "f1 = KalmanFilter(dim_x=4, dim_z=2)\n",
    "\n",
    "f1.F = np.array ([[0, 1, 0, 0],\n",
    "                  [0, 0, 0, 0],\n",
    "                  [0, 0, 0, 1],\n",
    "                  [0, 0, 0, 0]], dtype=float)\n",
    "\n",
    "f1.R *= 1.\n",
    "f1.Q *= .1\n",
    "\n",
    "f1.x = np.array([[1, 0, 1, 0]], dtype=float).T\n",
    "f1.P = np.eye(4) * 5.\n",
    "\n",
    "# initialize storage and other variables for the run\n",
    "count = 30\n",
    "xs, ys = [], []\n",
    "pxs, pys = [], []\n",
    "\n",
    "# create the simulated sensor\n",
    "d = DMESensor(pos_a, pos_b, noise_factor=3.)\n",
    "\n",
    "# pos will contain our nominal position since the filter does not\n",
    "# maintain position.\n",
    "pos = [0,0]\n",
    "\n",
    "for i in range(count):\n",
    "    # move (1,1) each step, so just use i\n",
    "    pos = [i, i]\n",
    "    \n",
    "    # compute the difference in range between the nominal track\n",
    "    # and measured ranges\n",
    "    ra,rb = d.range_of(pos)\n",
    "    rx,ry = d.range_of((pos[0] + f1.x[0, 0], pos[1] + f1.x[2, 0]))\n",
    "    z = np.array([[ra-rx], [rb-ry]])\n",
    "\n",
    "    # compute linearized H for this time step\n",
    "    f1.H = H_of (pos, pos_a, pos_b)\n",
    "\n",
    "    # store stuff so we can plot it later\n",
    "    xs.append(f1.x[0, 0]+i)\n",
    "    ys.append(f1.x[2, 0]+i)\n",
    "    pxs.append(pos[0])\n",
    "    pys.append(pos[1])\n",
    "    \n",
    "    # perform the Kalman filter steps\n",
    "    f1.predict()\n",
    "    f1.update(z)\n",
    "\n",
    "\n",
    "bp.plot_filter(xs, ys)\n",
    "bp.plot_track(pxs, pys)\n",
    "plt.legend(loc=2)\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Linearizing the Kalman Filter"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Now that we have seen an example of linearizing the Kalman filter we are in a position to better understand the math. \n",
    "\n",
    "We start by assuming some function $\\mathbf f$"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## References"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "[1] http://sympy.org\n"
   ]
  }
 ],
 "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.4.2"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 0
}