writeup.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Project 1 - Your First Feature (and your first bug).\n",
    "Writing code is important. But a big part of being on a self driving car team is working with a **large** existing codebase. On high stakes engineering projects like a self driving car, you will probably have to earn the trust of your managers and coworkers before they'll let you make substantial changes to the code base. \n",
    "\n",
    "A typical assignment for someone new to a team is to make progress on a backlog of bugs. So with that in mind, that's what you will be doing for your first project in the Nanodegree.\n",
    "\n",
    "You'll go through this project in a few parts:\n",
    "\n",
    "1. **Explore the Code** - don't worry about bugs at this point. The goal is to get a feel for how this code base is organized and what everything does.\n",
    "2. **Implement a Feature** - write code that gets the robot moving correctly.\n",
    "3. **Fix a Bug** - Implementing motion will reveal a bug which hadn't shown up before. Here you'll identify what the bug is and take steps to reproduce it. Then you'll identify the cause and fix it.\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Part 1: Exploring the code\n",
    "In this section you will just run some existing code to get a feel for what this localizer does.\n",
    "\n",
    "You can navigate through this notebook using the arrow keys on your keyboard. You can run the code in a cell by pressing **`Ctrl + Enter`**\n",
    "\n",
    "Navigate through the cells below. In each cell you should\n",
    "\n",
    "1. Read through the code. It's okay to not understand everything at this point. \n",
    "2. Make a guess about what will happen when you run the code. \n",
    "3. Run the code and compare what you see with what you expected. \n",
    "4. When you get to a **TODO** read the instructions carefully and complete the activity.\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "# This code \"imports\" code from some of the other files we've written\n",
    "# in this directory. Specifically simulate.py and helpers.py\n",
    "import simulate as sim\n",
    "import helpers\n",
    "import localizer\n",
    "\n",
    "# Don't worry too much about this code for now...\n",
    "from __future__ import division, print_function\n",
    "%load_ext autoreload\n",
    "%autoreload 2"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAD8CAYAAACMwORRAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4wLCBo\ndHRwOi8vbWF0cGxvdGxpYi5vcmcvpW3flQAAIABJREFUeJzt3Xl8FPX9P/DXO+dmswlnlECAlMOT\nokJE1HprvdEqWjyw1CqKon6lgketrUdttVU8f0U86oG2eKCCxaL1rFqQcGtBroIgV5BACLmzr98f\nn6WEsMnObGYzs7Pv5+MxD3ezn5l9+2bmvZ/9zHxmhSSUUkr5S5rbASillHKeFnellPIhLe5KKeVD\nWtyVUsqHtLgrpZQPaXFXSikf0uKulFI+pMVdKaV8SIu7Ukr5UIZbb9y1a1cWFxe79fZKKZWU5s2b\nt5VkQax2rhX34uJilJaWuvX2SimVlERkrZV2OiyjlFI+pMVdKaV8yLVhmXiQwOefA5MmAatWmb/1\n6QNccw1w3HGAiLvxec2aNcD/+3/A7NlAdTXQtStw6aXARRcBgYDb0XlLZSUwZQrwxhvA9u1AKASc\ncAJw7bVAt25uR+ct4TDw3nvAs88C69YB6enAQQcB118PDBrkdnTe85//AE8+CSxaBNTVAYWFwKhR\nwLnnAhmJrMAkXVkGDx5MO/7xD7K4mMzNJUVIU+rN41CI7N2bnDnT1iZ9a+1a8qSTyECAzMrakyvA\n5CoUIn/zG7Kx0e1I3VdbS44dS+bkmH2raa4CATI7mxw2jCwrcztSb5gyhdx/f7MPNc1VWhoZDJKH\nHEL++99uR+kNS5aQgwebfSsjY+985eWRnTuTTzxBhsP2tguglBZqbFIU9xdeMAlqmpxoS04O+dxz\n9hLlN8uWmZ0mPb31XAWD5Pnnkw0Nbkfsnpoa8phjYu9bmZlkjx7kd9+5HbG77rvP7DexjsNgUDta\nX3yx7wdgS7kaO9ZegfdNcf/Xv6wV9qYF/uOPrSfKTyoqyG7d9v5mE2vHGjfO7ajdc/HF1vetjAzy\ngAPIujq3o3bHq6/aOw6DQfLrr92O2h3r1pH5+fZy9dhj1rdvtbhbPqEqIukiskBE3onyWraITBWR\nlSIyR0SKnRo2uu02M15sVXW1WScVvfACsHOn2WWsqKoyY/Ll5YmNy4tWrwamT7e+bzU0ABs2ADNm\nJDYuLyKB8ePtHYc1NcC99yYuJi+bONH8/1tVVQXcdRdQX+9sHHaulrkJwNIWXvsFgHKS/QBMBPBA\nWwMDzEnTefPsr7doEbB8uRMRJA8S+OMfgV277K2XlgY8/3xCQvK0xx8HGhvtrVNZCTz4YGLi8bLP\nPgO+/97eOuEw8NZbqddxqKkBnn7anDi1o7HRdDacZKm4i0gRgLMBPNNCk/MAvBB5/DqAU0Tafu3K\nG2/YPwAB8wn42mttfffksnSp/QMQML2GZ591Ph6ve/nl+HpKCxYA27Y5H4+XvfSS/U4DYK4EeWef\n7/n+9skn8V21t3Mn8NxzzsZitef+CIAJAMItvN4DwDoAINkAYAeALm0NbuPG+A7A3V+hU0lZWfyX\nVW3d6mwsyaCiIr71srJMrlPJd99ZH+prqq4u9XJVVma+tcRj0yZnY4lZ3EXkHABbSLY2QBLts2qf\n3UFERotIqYiUlln4V8/OjtmkRal2HXdmpjvrJqv09PjWI02BTyXx/v+mpaVerjIz459v43SurPTc\njwUwTETWAPgbgJNFZEqzNusB9AQAEckA0AHAPl9eSU4mWUKypKAg5n1vcNBBZjKJXcGgWTeV9OkD\n1NbGt+6BBzobSzLo1Su+9cLh1JvUNHBgfB2tzEygXz/n4/Gy/v3j67mnpwOHHupsLDGLO8nbSRaR\nLAYwAsCHJC9v1mw6gJ9FHg+PtInji9zeLroovkSRwIgRbX335NKtG3DssfbXy8sDxo1zPh6vGzfO\nfschI8PsVzk5iYnJq66+Or71srOB005zNhavO+IIoHt3++tlZwM33OBsLHHfW0ZE7hGRYZGnzwLo\nIiIrAYwD4MjFiLm5wMiR9oYNMjKASy4xRSvVTJhgcmZHTg5w+umJicfLLr3UfschMxO4+ebExONl\nRUXA8cebYRarAgHgppviH/5KViLArbfaPw779QMOO8zhYKxcDJ+Ixeokpi1byMJCM7051mSAtDQz\nNXrTJusTAvwkHCYvuMD6ZJOcHHLWLLejds+LL1qbcbl7osn117sdsXtWriQ7dLA+4evAA8mdO92O\n2h21teRRR5lbV1jJV24uOX++9e3D6UlMbikoMDcL69699ZOkgYAZmvj8c2D//dsvPi8RAV55BTj1\n1NZ7Dunp5rzEiy8CP/5x+8XnNSNHmrkBwWDrJ8Fyc01P/7HH2i82r+nbF/joI6BTp9ZP/OXkmF7o\nJ5/Ed77MD7KyzI3VDj/c7Fstycw0IwwzZ5rhHMdZ+QRIxGL3xmHl5eTvfkd27WpuurP7Blh5eeZv\n995Lbttma5O+1dhI/vWv5GGHmd55Xp7pHeTnm5thXXFF6k4Nj+bLL803nkDA5Cg31+QsEDD3npkx\nw/7Nnfxq40ZywgSTp9252n0sFhWRjz5K7trldpTeUFtLPvUU2bfvnn0qFNqTt7Fjyf/+1/52YbHn\nLmz7ec+4lJSUMJ5fYmpsBD79FFi/3jzv0cPcmjXVxvas+uors1RXm17XSScBHTq4HZU3lZWZfaui\nwvTWBw1Kvas9rKqrMz35TZvMea7iYuCYY/S229GQQGkpsGKFyVtBgTkOW+vVt0ZE5pEsidku2Yq7\nUkqlMqvF3fNj7koppezT4q6UUj6kxV0ppXxIi7tSSvmQFnellPIhLe5KKeVDWtyVUsqHtLgrpZQP\naXFXSikf0uKulFI+pMVdKaV8yMpvqAZE5EsRWSQiX4vI3VHajBKRMhFZGFmuSky4SimlrMiw0KYW\nwMkkK0UkE8BnIvIuydnN2k0lOdb5EJVSStkVs7hH7h9cGXmaGVncuZWkUkopSyyNuYtIuogsBLAF\nwPsk50RpdqGILBaR10Wkp6NRKqWUssVScSfZSPJwAEUAhojIgGZNZgAoJjkQwD8BvBBtOyIyWkRK\nRaS0rKysLXErpZRqha2rZUhuB/AxgDOa/f17krWRp08DGNzC+pNJlpAsKSgoiCNcpZRSVli5WqZA\nRDpGHucAOBXAsmZtCps8HQZgqZNBKqWUssfK1TKFAF4QkXSYD4NXSb4jIvfA/FDrdAA3isgwAA0A\ntgEYlaiAlVJKxaa/oaqUUklEf0NVKaVSmBZ3pZTyIS3uSinlQ1rclVLKh7S4K6WUD2lxV0opH9Li\nrpRSPqTFXSmlfEiLu1JK+ZAWd6WU8iEt7kop5UNWbhzmOfX1wLZt5nHnzkBmprvxeFk4DGzfDtTU\nAB07AsGg2xF5W2UlsGMHkJsLdOgAiLgdkXfV1gLl5UB6OtCpE5CRlNWkfTQ2mlzV15tcBQKJf8+k\n6rnPnw9cdhmQlwcUF5slFAIuuQSYN8/t6Lxl0ybgt78FCgqAwkKgXz9TrIYMAd54A2hocDtC76ip\nAV56CRgwwBx4/fsD++8PdO8O/OlPezoSCiCBTz8Fzj3XHId9+gC9egH5+cA11wDLlsXeRir573+B\nm282HauiIqBvX5O3U04BZs0yna+EIenKMnjwYFpVU0MOH04Gg2R6Oml2sT1LWpp57Sc/MW1T3eTJ\nZCBglua5Asi8PLJ3b3LVKrcjdd/ChWRBARkKRc9VMGiWN990O1L37dhBHnecyZXIvrnKyCBzcsjr\nriMbGtyO1l3hMHnXXeYYzMqKvm+FQuTAgeSWLfa2DXOr9Zg11vPFvaGBPO00s9NES1DTJSeHPPlk\nsr7eXrL85MknTTGKlau0NLJLF3LtWrcjds+SJS0X9Wj71rRpbkfsnl27yB/+kMzOjp2rYJAcOdIU\nuFQ1fry14zAz03S0tm2zvm2rxd3zwzITJwKffw5UV8duW10NzJ4NPPRQ4uPyomXLgFtuAaqqYrfd\nPRZ/wQWJj8uLwmHgrLOAXbusta+uBi6/HNi8ObFxedX48cCKFWacPZaqKmDaNGDq1MTH5UUffgg8\n+aS147C+Hti4EbjqKufjsPIzewER+VJEFonI1yJyd5Q22SIyVURWisgcESl2Irhw2Ix5WknSblVV\nZp3GRiciSC4PP2x2FqsaG4GlS4HFixMXk1e9/745wUUbv1UTDgNPP524mLxq1y7g+efNuQk769x/\nf8JC8rT777dXs+rqgJkzne84WOm51wI4meRhAA4HcIaIDG3W5hcAykn2AzARwANOBPfee9Z7Vk3V\n1gLvvutEBMmjshKYMsX+idLaWvPtKNU8+KDJmR01NcCjj6Zex+Hll+O7amjlSmDhQufj8bJ164DP\nPotv3aeecjaWmMU9Msyz+zDIjCzN+zvnAXgh8vh1AKeItP0isi++sH8AAsDOnWbdVLJ0aXyXhDY2\nAh9/7Hg4njd3bnzrVVWZr9Gp5IMP4utkAcCcOc7G4nVz5wLZ2fbXq6kB/vlPZ2OxNOYuIukishDA\nFgDvk2z+T9YDwDoAINkAYAeALlG2M1pESkWktKysLOb7bt9uJbroUu3ytcrK+K/JtvMV0i/sDDE0\nlZ5uOg+ppKIivvUaGlIvV5WV8V/e6HSuLBV3ko0kDwdQBGCIiAxo1iRaWdlnNJPkZJIlJEsKCgpi\nvm/nzvEXrC77fLT4W36+vfHjpkIhZ2NJBjk58a1XX29ynUo6dYpvvczM1MtVfj6QFudlKh06OBuL\nrTBIbgfwMYAzmr20HkBPABCRDAAdALS573zyyWamoF2hkJkkkEoOPTS+4p6ZCZx5pvPxeN2JJ8bX\ncejc2UxuSiVnnx1fB4AETjjB+Xi87JhjzAlSu4JBMzHMSVaulikQkY6RxzkATgXQfB7adAA/izwe\nDuDDyPWYbXLccfH1wDt2BE46qa3vnlwCAXM5VVaWvfXS04GbbkpMTF52yy32b8UQDJr1Uu2WBMOH\nx/f/PHAgcOCBzsfjZfvtB5x+uv18hcPAlVc6G4uVnnshgI9EZDGAuTBj7u+IyD0iMizS5lkAXURk\nJYBxAG5zIjgR4Fe/sncQBoPAHXek3gEIADfeaO/+HllZwNFHm+n2qeZHPwJ69jQfblalpwM//3ni\nYvKq7Gxg7Fh7x2FuLvDrXycuJi+7/XZ7w36BADBiRPzDXy2yMtMpEYvVGarhMHnppdZmewWD5E9/\nmtoz41591dps3sxMsmdPcutWtyN2z5o1ZOfOZraulX3ro4/cjtg99fXkCSdY27eCQTNDM5U99JC1\nmhUImFsQVFVZ3zb8MkNVxNzU6eqrzSdctMuMsrLMa1deGf81uX5x0UVmZmBubvRxUhHz2g9/CJSW\npt6J56Z69zY3nOvTp+Ux5VDIDPPNmmXG6VNVRobJwbnnml5ptMtuAwGz3HEH8IAjM12S17hxZv5I\nIBD9G096uvn7cceZGfjxnuBvlZVPgEQsdm4cttu335K33WbuiZKWZpbOnckJE1L7HinRVFaSzzxD\nHnDAnput5eSQF1xAfvFFan+7aa6xkXz/fXMPo933TsnIIA8/nHzlFb0ZXXPffGNuDpafb24glpZG\ndutG3ncfuXmz29F5S3k5+fDD5ltyWprJV24uOWoUuWhRfNuExZ67sO3nPeNSUlLC0tLSuNfffS1p\nvJcdpRLS5MvO+HIqa2zUXFkVDptvg6n8bdkqp45DEZlHsiRWu6S9vb4WdetEtFjZobmyTo9D69r7\nONR/GqWU8iEt7kop5UNa3JVSyoe0uCullA9pcVdKKR/S4q6UUj6kxV0ppXxIi7tSSvmQFnellPIh\nLe5KKeVDWtyVUsqHrPwSU08R+UhElorI1yKyz+/2iMiJIrJDRBZGlrsSE65SSikrrNw4rAHAL0nO\nF5E8APNE5H2S/2nW7l8kz3E+RKWUUnbF7LmT3EhyfuTxTgBLAfRIdGBKKaXiZ2vMXUSKARwBYE6U\nl48WkUUi8q6IHOpAbEoppeJk+X7uIhIC8AaA/yNZ0ezl+QB6k6wUkbMAvAVgn59dFpHRAEYDQK9e\nveIOWimlVOss9dxFJBOmsL9Mclrz10lWkKyMPJ4JIFNEukZpN5lkCcmSgoKCNoaulFKqJVaulhEA\nzwJYSvLhFtp0i7SDiAyJbPd7JwNVSillnZVhmWMBjASwREQWRv52B4BeAEByEoDhAMaISAOAagAj\n6NaPsyqllIpd3El+BqDVn78l+QSAJ5wKSimlVNvoDFWllPIhLe5KKeVDWtyVUsqHtLgrpZQPaXFX\nSikf0uKulFI+pMVdKaV8SIu7Ukr5kBZ3pZTyIS3uSinlQ1rclVLKhyzfz90rdu0C3n4bWL/ePO/R\nAzjvPCAUcjcuLwqHgQ8/BJYsAaqrgU6dgLPOAnr3djsyb1q2DPjgA2DHDiA3J4ySIWk45hhAWr2z\nUmoqLzfH4ebNQHo6UFwMnHsukJ3tdmTe09AAvPsusHw5UFcHFBSYXO2/f4LfmKQry+DBg2nHunXk\nNdeQwSCZl0dmZJglFCJzcsirrya//dbWJn2rtpb805/Ibt1MfrKyyLQ0k6dAgDzpJPLTT92O0jv+\n/nfyyCP35KerbOV36M6i4PcsLiYnTSIbG92O0huWLSMvucTkKTeXTE83x2Fenll++Uty61a3o/SG\nigry178mO3UyucnMNPkKBk3+hg0jFyywv10ApbRQY5OiuC9caBKUkWEijrZkZJAdO5Lz59tPlp9U\nVJBHHWUKVUu5AswO9sQTbkfrvl/9yuSiaW5G4TkS4M/wl//l6qyzzIdmKvvwQ1PQ09Ja3q+yssju\n3clVq9yO1l2bNpH9+5PZ2S3nSsTsW2+8YW/bvinua9eawt5aoWq6dOxI/ve/9pLlFw0N5PHHt75D\nNV1ycsi//c3tqN3z0EP7FnaA/BjHkwA/wgl75eqii8hw2O2o3bFwoSnsVvartDSyR4/U7cHv2kUe\ndFDrndHmx+HHH1vfvtXibuWXmHqKyEcislREvhaRm6K0ERF5TERWishiERnk1LDR7bcDFc1/sbUV\nFRXArbc69e7JZcYMYP58oLbWWvvqauDaa4H6+sTG5UXbtwO/+hVQVbX33/NQgaMwGwAwFLORB7Pz\nVVcDf/878OWX7R2pN4wZY853WREOA2VlwB//mNiYvOqZZ4BvvzVj7VZUVwNXXWVKvZOsXC3TAOCX\nJA8GMBTA9SJySLM2Z8L8IHZ/mB/A/rMTwZWXA9OmAY2N1tcJh4Hp04HvU/BH/h54AKistLdOYyPw\n1luJicfLnn8e6CzlGIhFey1X4RnUwZwVrEMWrsIz/3vtgOpFeO3ORWbHTCErVgALFthbp64OmDTJ\n/DeVkOZDrXmnIZaNG4E5cxwPxt5wCoC3AZzW7G9PAbikyfNvABS2th0rwzKTJkX/2mzla06qjSd/\n+605SWM3VwB57LFuR9/++vQh/4AJJMAqBLgd+f9bmian6d+rYBJce/MEt8NvVxMmmJOBdvervDzy\n7bfdjr59ffGFuYjBbq7S0sjLL7f2HnBqWKYpESkGcASA5p8xPQCsa/J8feRvbbJypf1PQMB8zVmx\noq3vnlzWro3/MrQ1axwNJSls2ADcjt/jdtwPQtABFf9bmmr6d0Lwm8z7sW7M712K2h3LlsU3dFdX\nl3r71po18V06Gw6bSyWdZLm4i0gIwBsA/o9k81HwaP87+4wgichoESkVkdKysrKY72lnOMbJdZNR\nOOzOusmKBIg0/AG340R8jI3ohmoEoratRjY2oBAn4BM8nH07wik29y/eY4lMvX0rHI5/7NzpmmVp\nLxWRTJjC/jLJaVGarAfQs8nzIgAbmjciOZlkCcmSgoKCmO/bs2d8vdGsLKBXL/vrJbPCwvhPjCZ8\nMoUHde685/FcDMGB+AYbURi17UZ0x0FYhlIc+b9JKKmkuBhIi+PzLDsb6NbN8XA8rbAwvlwBpt45\nycrVMgLgWQBLST7cQrPpAK6IXDUzFMAOkhvbGtzFF8f3FSctDfjpT9v67smlf//4do7cXHPFTKq5\n8sq9Ow51yEIhou+y3bAJtZGTrMceC3Ts2B4ReseoUUAg+peaVjU0AOec43g4nnb88WbGrl2hEDB6\ntLOxWPmMORbASAAni8jCyHKWiFwrIrvLwkwAqwGsBPA0gOucCK5HD+CEE+yvd+yxqddzB8wloLm5\n9tYhgcsuS0w8XjZmzN7PT8P7qEMWAGAXgqhFFnYhCACoRyZOw/sIhYAJE9o7UveVlNg/njIyzH6V\narcFycgAxo61P+IQDAKnn+5wMFbOuiZisTqJacECe1fMBIPkvHnWzjr7TVUV2bev9ckTwSD5u9+5\nHbV7xozZs29NxUUMA6xEkI/gBuZiJx/FWFYih2GAr8rFPPLI1L0NwaxZsWc9N13y88nVq92O2h1b\ntpAFBWYGqtXj8OWXrW8ffpmhSpIzZ1or8Dk55IwZ1pPkR+vXk0VFsS9dCwbJa69N3RmXpJnRe/75\nZIecWlYiyO3I54/xj73ydDre5Q7kc5fkctvmOrdDdtUzz8Qu8CLmEsgvvnA7Wnd99RXZubO5l0ys\n4/D3v7e3bV8Vd5KcO5c8+mhzLXfTwpWZaXa4o44iZ8+2lyS/2rqVHDnS5Kr5h2IoZO79MXlyahf2\n3RobyYfuq+Iz2WPYJ3fTXrnKzTX5u3XUJtZdPcZ8NUpx779P/vCHJi9N7zGTnW32t1NPJZcudTtK\nb1i71twcLBDYew6KiNm3+vYlp02zv12rxV1M2/ZXUlLC0tJS2+utWAH85S/AqlUmVX36AD//OXDg\ngQkIMslt3w68+CIwd66ZOr7ffsDw4cDJJ8d/Rt+vGhrM7QVmzDCzm/PzzcmxESPsn8dIBQsXmn3r\n22+BzEzgoIPMSWqnr/jwg82bTc1asgSoqQG6dzfnI446Kr4LRkRkHsmSmO2SrbgrpVQqs1rctf+m\nlFI+pMVdKaV8SIu7Ukr5kBZ3pZTyIS3uSinlQ1rclVLKh7S4K6WUD2lxV0opH9LirpRSPqTFXSml\nfEiLu1JK+ZAWd6WU8iErP7P3nIhsEZGvWnj9RBHZ0eRXmu5yPkyllFJ2ZFho8zyAJwC82Eqbf5FM\nsV9LVEop74rZcyf5KYBt7RCLUkophzg15n60iCwSkXdF5NCWGonIaBEpFZHSsrIyh95aKaVUc04U\n9/kAepM8DMDjAN5qqSHJySRLSJYUFBQ48NZKKaWiaXNxJ1lBsjLyeCaATBHp2ubIlFJKxa3NxV1E\nuomYXwIUkSGRbX7f1u0qpZSKX8yrZUTkrwBOBNBVRNYD+A2ATAAgOQnAcABjRKQBQDWAEXTrh1mV\nUkoBsFDcSV4S4/UnYC6VVEop5RE6Q1UppXxIi7tSSvmQFnellPIhLe5KKeVDWtyVUsqHtLgrpZQP\naXFXSikf0uKulFI+pMVdKaV8SIu7Ukr5kBZ3pZTyoaQq7rW1wCuvAEcfDRQVmWXoUGDKFKCmxu3o\nvIUEPv8cuOAC4Ac/ALp3Bw4+GLjnHmDzZrej8561a4Hx44EDDgAKC4E+fYCRI4EFC9yOzHsqK4Gn\nngIGDQJ69AB69QJOPBF46y2gocHt6LwlHAZmzQLOOAPo3dvka+BAYOJEYPv2BL85SVeWwYMH045n\nniHz8sxiSteeJRQyf580ydYmfWvBArJvXzI3lxTZO1eBAJmdTY4cSVZXux2p+3bsIM8+2+QlK2vv\nXKWnk8EgOXAguWqV25G6Lxwm773X5CQ3d9/jMC+P7NyZfPNNtyP1hg8/JAsLo9esYNDsc+PGkQ0N\n9rYLoJQWamxSFPf77jPJaJ6gaAn77W/tJcpvZs+OfuA1X3JyyCFDUrvAl5eT/fubD7vWcpWWRnbs\nSC5d6nbE7gmHyauvtnYc5uSQzz7rdsTumj7des067zyysdH6tn1T3GfMsJakpsl66y3rifKTsjKy\nQwfrucrJIS+7zO2o3XPSSbEL++5FhOzWjayqcjtqdzz5pP3jcPZst6N2xzffWOtgNc3VnXda377V\n4h5zzF1EnhORLSLyVQuvi4g8JiIrRWSxiAxyctjozjuBqirr7auqzDqpaPJkoK7OevvqauCNN4BN\nmxIXk1ctWQLMnm3O41hBmrHm115LbFxe1NgI3H23veOwutqc30lFDz5o7xxgVRXwyCMmZ06yckL1\neQBntPL6mQD6R5bRAP7c9rCMxYuBFSvsr7d6deqdCGtsjH8Heeop5+PxuokT7X0QAqa4P/BAYuLx\nslmz7O9XJPDBB8CGDYmJyat27jQXfTQ22l936lRnY4lZ3El+CmBbK03OA/Bi5BvDbAAdRaTQieDe\necf+AQiYT83p052IIHl89VV8hb2mxuyMqebtt+M7AFetArZscT4eL3vtNVO07MrIMB8MqeSTT4DM\nTPvrVVY6fxw6cSlkDwDrmjxfH/lbm23eHN+lVeFw6g01lJcD6enxrbtjh7OxJIPKyvjWy8oCtrXW\n1fGheI+lurrUy1V5uak/8di61dlYnCjuEuVvUX8gW0RGi0ipiJSWlZXF3HBubvxB5eXFv24yCgTi\nXzc727k4kkU8vSvAHLg5Oc7G4nXBYHzrpaenXq4CASAtzqrqdK6cKO7rAfRs8rwIQNSRNpKTSZaQ\nLCkoKIi54cMOA0Ih+wGFQmbdVNK/f3wTuUSAww93Ph6vO/DA+NYTMZOcUslRR8VXeDIygAEDnI/H\nywYMAOrr7a+XmQkceaSzsThR3KcDuCJy1cxQADtIbnRguzj//Pg/BS+80IkIkkeXLsBZZ5niY0cw\nCPzyl4mJycvGj7f/7S4zE/jFL8zQTCq58kpzgtSuTp2A445zPh4vO/hg4KCD7K+Xng6MHetsLFYu\nhfwrgH8DOFBE1ovIL0TkWhG5NtJkJoDVAFYCeBrAdU4Fl50NXHONvSGH7Gzg6qvbNkyRrMaPt9/D\nKihIvQMQMLdlsNtxSE8HbrghMfF4WdeuwNln2xvKCgbN/mi3s+EHt91mb0g5LQ0oKQH69XM4ECsX\nwydisTqJaedO8qCDyMzM2JMBMjPJAw4gKyqsTwjwmxtusD7ZJDeXnDfP7Yjd8957ZiKX1YkmDzzg\ndsTu2bTJTOJKS4udq0CA/NFO6khbAAAOc0lEQVSPyLo6t6N2R2OjmXVqZd8SMbOfV6+2vn04NYnJ\nbaEQ8NlnwCGHtP5pmJtrxlE/+yz1TqY29cgjZuggGGy51xQIAPn55jK1QY5OOUsup50GvPqq2Xda\nGmrZfVLwzjtNTzRV7b+/mfRVVNT6cRgKmRv7vftu/Cetk11amtmvzjyz9VwFg+Zb0eefm5v7Oc7K\nJ0AiFrs3DqutJV9+2dzEKSfHTLPv0ME8HjCAnDKFrKmxtUlf+/xz8vzzzfT63bnKyyO7dCHvuYfc\nvNntCL1j7Vpy/HiTn/z8PfkKBMjLLzc3YlNGZaW5QV+fPuabX9NcDR1qbv1h90ZYfhUOk//4B3ny\nyXsfh6EQ2b07OXEiuX27/e3CYs9dGM+ZEgeUlJSwtLQ0rnVXrQLWrzePe/RIwFiVj2zdCixfbiY4\ndepkbjeakeF2VN5UW2tmRVdUmB7XwQcDHTq4HZU3kcCyZeYa+IwMczvbXr3cjsq7NmwwM+fr6sx5\nrgED4j8fISLzSJbEbJeMxV0ppVKV1eLu+TF3pZRS9mlxV0opH9LirpRSPqTFXSmlfEiLu1JK+ZAW\nd6WU8iEt7kop5UNa3JVSyoe0uCullA9pcVdKKR/S4q6UUj5kqbiLyBki8o2IrBSR26K8PkpEykRk\nYWS5yvlQlVJKWRXz/oAikg7gSQCnwfxe6lwRmU7yP82aTiXp8A9FKaWUioeVnvsQACtJriZZB+Bv\nAM5LbFhKKaXawkpx7wFgXZPn6yN/a+5CEVksIq+LSE9HolNKKRUXK8U92i3lm98EfgaAYpIDAfwT\nwAtRNyQyWkRKRaS0rKzMXqRKKaUss1Lc1wNo2hMvArChaQOS35OsjTx9GsDgaBsiOZlkCcmSgoKC\neOJVSillgZXiPhdAfxH5gYhkARgBYHrTBiJS2OTpMABLnQtRKaWUXTGvliHZICJjAcwCkA7gOZJf\ni8g9MD/UOh3AjSIyDEADgG0ARiUwZqWUUjHob6gqpVQS0d9QVUqpFKbFXSmlfEiLu1JK+ZAWd6WU\n8iEt7kop5UNa3JVSyoe0uCullA9pcVdKKR/S4q6UUj6kxV0ppXxIi7tSSvlQzBuHedGGDcB335nH\nPXoA3bu7G4+XVVQAq1cDNTVAp05A//5Amn6kR9XQACxfDuzYAeTmAv36AcGg21F519q1wObNQHo6\n0LMnsN9+bkfkXdu2AWvWAPX1QNeuQJ8+gET7pQwnkXRlGTx4MO1oaCDffJMcOpQMBMj8fLNkZ5NH\nHUVOm0bW19vapK8tWEBefvmeXHXoQObmkoWF5J/+RG7b5naE3rFpE3n33WSXLmQoZHKVn08Gg+S1\n15LLlrkdoXfU1JAvvUQOGEDm5OzJVXY2ecop5HvvkeGw21F6QzhM/utf5LBhJj+7j8OcHLJvX/Kp\np8jKSvvbhbkbb8wamxTFvbzcFPBQyEQcbQmFyJISLVrhMHn77WYHSk+PnqtgkOzYkfzyS7ejdd97\n75l9JxCInquMDJPLRx5xO1L3rV9P9uvX+nGYm0ueeSZZVeV2tO6qrydHjjT5EGk5V4WF5PLl9rbt\nm+JeVUUedpj55Gtph9q9ZGWZHsWuXfaS5Se33mp2mli52v2BuHix2xG755NPzAedlVwFg+Rjj7kd\nsXu2biWLisyHXaxc5eSYXnxDg9tRuyMcJkeMsLZvpaWZb4zr1lnfvtXi7vnR1/vvB775Bqitjd22\nrg5YuRK4557Ex+VFpaXA448Du3ZZa19ZCVxwgdnNUk19PfCTnwBVVdbaV1UBEyaYceZUdOONZny9\noSF22+pq4N//Bp5+OvFxedFbbwEzZljbt8JhYPt24IornI/DUnEXkTNE5BsRWSkit0V5PVtEpkZe\nnyMixU4EV18PPPGEORloVU0N8Oc/m0Kfav74R3u5AoCNG4HZsxMTj5dNn272LzvCYeDJJxMTj5eV\nlwPTptnLV1UV8OCDqdlx+MMfrHewAKCx0XwYrlnjbBwxi7uIpAN4EsCZAA4BcImIHNKs2S8AlJPs\nB2AigAecCG76dPM/blc4DLz5phMRJI/ycpOvcNjeelVVwEMPJSYmL3vwQWDnTnvr1NUBTz1l/0Mh\n2f3lL/FdYbVliylaqWT5cmDJEvvrJaLjYOWfbAiAlSRXk6wD8DcA5zVrcx6AFyKPXwdwikjbL/SZ\nP9/+AQiY4Yb589v67sll+XIgO9v+eqQZzkk1//lPfOs1NppvO6nkiy+sD1811dgILFrkfDxetngx\nkJlpf726Ouc/CK0U9x4A1jV5vj7yt6htSDYA2AGgS/MNichoESkVkdKysrKYbxxPYd+toiL+dZNR\nPAffbtXVzsWRLKycw4kmLa1tuU5GdoYYmmpsTL1cVVXZ//a8W7x5bomV4h6tB958JM1KG5CcTLKE\nZElBQUHMN95vPzNBwi6R1JtQ0alT/DtVhw7OxpIMQqH41qurAzp2dDYWr7NwqEaVlZV6uerUKf5J\ngp07OxuLlTDWA+jZ5HkRgA0ttRGRDAAdAGxra3BnnRXfUEMwCJxzTlvfPbkMGBDf18HsbODii52P\nx+vOPju+jkPv3sD++zsfj5cNHw7k5dlfr6EBOP105+PxsuOPj+9ijlAIuOQSZ2OxUtznAugvIj8Q\nkSwAIwBMb9ZmOoCfRR4PB/Bh5HrMNhk0CCgutr9eURFw5JFtfffkkpFhLlcLBOytJwKMGZOYmLxs\n3Dj7HYfcXHM5ZMKnjXvM2WebXrhdxx9vjsVU0qEDcNFF9jsO4bALxT0yhj4WwCwASwG8SvJrEblH\nRIZFmj0LoIuIrAQwDsA+l0vG69577d3fIxg066Sia6+113sPBMw3nB7Nz6CkgCOOAA4/3F7Ryslx\n/gBMBunpwB132D8O77orcTF52a232us4BIPAddeZzoOjrMx0SsRi5/YDt99ubbZXMEjecovlzfrS\nxx9by1UgQA4cmNqzebduJXv3NjObW8uVCJmXRy5a5HbE7gmHyUsvtX4cPv642xG7a+pUM1PXSq5+\n/GN798WCX24/sNvEiSYR0abW5+aa1x56yNYmfWvOHLJ79+j3AMnKMoX9Jz+J76ZFfvP99+SJJ5oD\nsfnUehGzb/XvTy5d6nak7mtsNLe3CASiF/lQyHwITpnidqTe8Pe/k507m5xE61wFAuTo0fZveGi1\nuItp2/5KSkpYavMC6127gFdeMbNWN2wwaSosBMaOBS67LP4rIPwoHAY++MDMWl282Jzkyc83J0+v\nv96cGFR7fPMN8OijZtr4rl1myGroUOCWW4Cjj069cfbWlJebiU2TJwNlZWbYpndv4OabgQsvjO8i\nCL9qaDD71EMPAStWmAlwnTsDo0YBV18d38l5EZlHsiRmu2Qq7kopleqsFnfP3zhMKaWUfVrclVLK\nh7S4K6WUD7k25i4iZQDaenfsrgC2OhCO07wYlxdjAjQuO7wYE6Bx2eFETL1JxrwphGvF3QkiUmrl\nxEJ782JcXowJ0Ljs8GJMgMZlR3vGpMMySinlQ1rclVLKh5K9uE92O4AWeDEuL8YEaFx2eDEmQOOy\no91iSuoxd6WUUtEle89dKaVUFElR3EXkDBH5RkRWisg+txMWkWwRmRp5fY6IFHsgplEiUiYiCyPL\nVe0Q03MiskVEvmrhdRGRxyIxLxaRQYmOyWJcJ4rIjia5SvjNYkWkp4h8JCJLReRrEbkpSpt2z5fF\nuNzIV0BEvhSRRZG47o7Spl2PQ4sxtftx2OS900VkgYi8E+W1xOfKyt3F3FwApANYBaAPgCwAiwAc\n0qzNdQAmRR6PADDVAzGNAvBEO+fqeACDAHzVwutnAXgX5mcRhwKY45G4TgTwTjvnqhDAoMjjPADL\no/wbtnu+LMblRr4EQCjyOBPAHABDm7Vp7+PQSkztfhw2ee9xAF6J9m/VHrlKhp77EAArSa4mWQfg\nbwDOa9bmPAAvRB6/DuAUkYTex89KTO2O5Kdo/ecNzwPwIo3ZADqKSKEH4mp3JDeSnB95vBPmh2ia\n/2xJu+fLYlztLpKDysjTzMjS/IRdux6HFmNyhYgUATgbwDMtNEl4rpKhuPcAsK7J8/XYd2f/Xxua\nX47aAaCLyzEBwIWRr/Ovi0jPKK+3N6txu+HoyNfrd0Xk0PZ848hX4iNgen5NuZqvVuICXMhXZJhh\nIYAtAN4n2WK+2uk4tBIT4M5x+AiACQBa+tn6hOcqGYp7tE+z5p/OVto4ycr7zQBQTHIggH9iz6e0\nm9o7T1bNh5lSfRiAxwG81V5vLCIhAG8A+D+SFc1fjrJKu+QrRlyu5ItkI8nDARQBGCIiA5o1afd8\nWYip3Y9DETkHwBaS81prFuVvjuYqGYr7egBNP22LAGxoqY2IZADogMQOA8SMieT3JGsjT58GMDiB\n8VhlJZftjmTF7q/XJGcCyBSRrol+XxHJhCmgL5OcFqWJK/mKFZdb+Wry/tsBfAzgjGYvtfdxGDMm\nl47DYwEME5E1MEO2J4vIlGZtEp6rZCjucwH0F5EfiEgWzMmH6c3aTAfws8jj4QA+ZORMhVsxNRub\nHQYzduq26QCuiFwFMhTADpIb3Q5KRLrtHm8UkSEw++X3CX5Pgflh96UkH26hWbvny0pcLuWrQEQ6\nRh7nADgVwLJmzdr1OLQSkxvHIcnbSRaRLIapDR+SvLxZs4TnKsPJjSUCyQYRGQtgFsxVKs+R/FpE\n7oH5LcHpMAfDSyKyEubTb4QHYrpRRIYBaIjENCqRMQGAiPwV5kqKriKyHsBvYE4ygeQkADNhrgBZ\nCaAKwM8THZPFuIYDGCMiDQCqAYxI8IczYHpXIwEsiYzZAsAdAHo1icuNfFmJy418FQJ4QUTSYT5M\nXiX5jpvHocWY2v04bEl750pnqCqllA8lw7CMUkopm7S4K6WUD2lxV0opH9LirpRSPqTFXSmlfEiL\nu1JK+ZAWd6WU8iEt7kop5UP/H31hrXWveiVWAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7fce55bf4f50>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# This code defines a 4x4 robot world as well as some other parameters\n",
    "# which we will discuss later. It then creates a simulation and shows \n",
    "# the initial beliefs. \n",
    "R = 'r'\n",
    "G = 'g'\n",
    "grid = [\n",
    "    [R,G,G,G,R],\n",
    "    [G,G,R,G,R],\n",
    "    [G,R,G,G,G],\n",
    "    [R,R,G,R,G],\n",
    "    [R,G,R,G,R],\n",
    "]\n",
    "blur = 0.05\n",
    "p_hit = 200.0\n",
    "simulation = sim.Simulation(grid, blur, p_hit)\n",
    "simulation.show_beliefs()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Run the code below multiple times by repeatedly pressing Ctrl + Enter.\n",
    "\n",
    "After each run observe how the state has changed. "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "NOTE! The robot doesn't have a working sense function at this point.\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAD8CAYAAACMwORRAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4wLCBo\ndHRwOi8vbWF0cGxvdGxpYi5vcmcvpW3flQAAIABJREFUeJzt3Xl8FfW5+PHPkz0nCXuACCguaFWs\nC1FxqSKuaCvVeq2+6lI3qq1X+9O2Wlvr0t629nJtr7Xqyxattnrd20YLVaxYtIoYlLghFRUFAYns\nIRCSnOf3xzepIR6SmckcZs6c5/16nZcnySwPjzPP+Z7vfOc7oqoYY4xJloKoAzDGGBM+K+7GGJNA\nVtyNMSaBrLgbY0wCWXE3xpgEsuJujDEJZMXdGGMSyIq7McYkkBV3Y4xJoKKodjxkyBAdPXp0VLs3\nxpicNG/evE9Utbq35SIr7qNHj6a+vj6q3RtjTE4SkQ+8LGfdMsYYk0BW3I0xJoEi65YJiyq8+CLc\ncQcsWuR+Hj0avvENOPJIEIk6wvh4+WW47TZYuBDSadhxR7joIjj6aCiwj/l/mz8ffvMbeOstaGuD\nkSPhggvg+OOhsDDq6OLjjTdcnl5/HbZsgR12gPPOg5NOgqKcryzhWbgQbr3VHVctLTBsGJx7Lkye\nDMXFWdyxqkbyGjdunPbV00+r7rKLakWFqoiqK+3ufWWl6qhRqnV1fd5NznvuOdU99nB5Kij4NE/g\n8lRTo/rgg1FHGb25c1XHjlVNpVQLC7fOU1WV6tChqr//fdRRRm/+fNX991ctL8+cp8GDVe+4I+oo\no7dggerBB7s8FRV9Nk8DBqjefLNqOu1vu0C9eqixOVvc77/fJa1rwjK9ysvz+0D7y1+85SmVUp06\nNepoo/PUUy4HXvJ0/fVRRxud555zjQQvebriiqijjc7cua6Ad210bitPU6b4K/CJLu5z5ng7EbsW\n+JkzA+8uZzU0+MtTKuU+DPLNwoXeClbXPN13X9RRb3+LF7uC5SdP+diwWr7ctcr95OkXv/C+fa/F\n3XNPq4gUisirIvJEhr+VisiDIrJIRF4SkdGh9RtlcM010NzsfflNm+Cqq7IXT1xdd537t3vV3Azf\n+Y475PLJT34Cmzd7X765Gb73vfzL09Sp/vP0wx9Ce3v2YoqjW27xf979+MeuPz5Mfi6jXQ4s2Mbf\nLgDWqOpuwC+Bm/oa2LZ88AG88IL/9RYsgDffDD+euFq5Ev72N/8FaNkyd+E1X6xbBw8/7L8ArVsH\ns2ZlJ6Y4am6Gu++G1lZ/67W0wPTp2Ykpjlpb3aCFIIX6scfCjcVTcReRkcBJwO+2schk4J6O948A\nR4tkZ5zKn/4UrMXU2goPPhh+PHFVVxdsZMemTXDffeHHE1czZgQbsdDUBPfeG348cfXMM8GOpw0b\n3IdCvvjnP4PVpw0bYNq0cGPx2nL/FfA9IL2Nv48AlgCoahuwDhjc5+gy+PjjYJ+KbW3w0UfhxxNX\njY3+vkJ3Sqdh6dLw44mrxkb/rdFO+XY8Be1eWbYs3FjirLExeHfdihXhxtJrcReRLwIrVXVeT4tl\n+N1n/okiMkVE6kWkvrGx0UeYnyopCbQaAGVlwdfNNcXFwceu51uegn7HLC0NN5Y460ue+nLO5po4\n5cnL6X8YcLKILAYeACaKyB+7LbMUGAUgIkVAf2B19w2p6p2qWquqtdXVvc57k9Huu0Nlpf/1UinY\na69Au8xJY8YEK9JlZTB2bPjxxNWYMcFuuCkpyb88BWmRFhXBPvuEH09cjRnjegn8KiiAvfcON5Ze\ni7uqfl9VR6rqaOAM4BlVPavbYnXAuR3vT+tYJitjCU49Ndh66TR87WvhxhJnJ54Y/G7K884LN5Y4\nO+oo98HvV0EBXHxx+PHE1UEHQZD2WHExXHpp+PHE1d57wy67+F+vvBwuvzzcWALfdC4iN4rIyR0/\nTgMGi8gi4Arg6jCCy6S83BUfP19hiorgK1+BAQOyFVX8FBfDt77lr+ugoACOPRaGD89eXHFTUABX\nXOGOK69EYPx4N81FvhBxwz/9fhDutRfsuWd2Yoqrq66Cigp/64waBbW14cYhWWpg96q2tlaDTvm7\nejXsuy8sX977RR4RGDIEGhqgpibQ7nLW+vVwwAFu+KiXr4oDB8Krr8JOO2U/tjhpbnYt03/9y9vF\n1X793HDR3XfPfmxx0tIChx8Or73m5pLpTUWFm/cpn7plwJ1rEyfC3LneBn+kUvDss3Dggd62LyLz\nVLXXj4KcnC5q0CB4/nk3oVNPLa6yMhg61C2bb4UdXBF67jnXwuwpT6WlMHgw/OMf+VfY4dOTa/fd\ne26ZlpS4b39//3v+FXZwx8nMma5Y95Sn4mJ37D35ZP4VdnA9BdOnu2LdUwu+qMhdP6yr817Y/cjJ\n4g6uCDU0uLswhw6FqiqXqMpK937QILj6anfjUj6eiJ1qalxr/Cc/ce87c9SZpwED4Mor3QyI+Xgi\ndhoyxLXGb7rJfUXuzFFFhStUVVVw2WXueAr763MuGTDAjeX+n/+BnXd2+ek89/r1c//95jfdTJGH\nHRZ1tNGprHT3Bvz6167+dOap83hKpdyMrA0NblbWbMjJbpnu0mnXQl2yxL0fMcJN92vTjm4tnXYn\n5gcfuPfDh7sLilmddjQHqcKcOfDee+4r9rBhLk/5NPTRC1X3gbhokeumGTrU5cnP9Yt8oOoaWAsX\num6awYNdt43ffvlOXrtlElHcjTEmXyS6z90YY0zPrLgbY0wCWXE3xpgEsuJujDEJZMXdGGMSyIq7\nMcYkkBV3Y4xJICvuxhiTQFbcjTEmgay4G2NMAllxN8aYBPLyDNUyEZkrIg0i8qaI3JBhma+LSKOI\nzO94XZidcI0xxnjhZd7EFmCiqjaJSDHwvIjMUNU53ZZ7UFXz6IFaxhgTX70W945noTZ1/Fjc8Ypm\nKkljjDGeeOpzF5FCEZkPrARmqupLGRb7ioi8JiKPiMioUKM0xhjji6firqrtqrofMBI4SETGdlvk\ncWC0qn4eeBq4J9N2RGSKiNSLSH1jY2Nf4jbGGNMDX6NlVHUt8CxwQrffr1LVzkfB/hYYt43171TV\nWlWtra6uDhCuMcYYL7yMlqkWkQEd78uBY4C3uy3T9fHTJwMLwgzSGGOMP15Gy9QA94hIIe7D4CFV\nfUJEbgTqVbUOuExETgbagNXA17MVsDHGmN7ZM1SNMSaH2DNUjTEmj1lxN8aYBLLibowxCWTF3Rhj\nEsiKuzHGJJAVd2OMSSAr7sYYk0BW3I0xJoGsuBtjTAJZcTfGmASy4m6MMQnkZeKwnNDaCmvWgCoM\nGgTFxVFHFE9tbbB6NaTTLk8lJVFHFE9tbe54amtzeSotjTqieGpvd3lqbYWBA6GsLOqI4qm9Hdau\nhZYWl6fy8uzvM+db7g0NcM45UFUFO+0Eo0dDZSWcfjrMnRt1dPHx5ptwwQWf5mmXXVyeTjkF/vlP\n96Fo4F//gksugf79YdQo2HVXl7OTToJnn7U8dXrvPbj8cpenkSM/zdOxx8LTT7vGg4EPP4TvftcV\n9BEjYLfdoF8/mDABZszIcp5UNZLXuHHjtC9aWlTPPFM1lVItLFR1p92nr4IC1YoK1ZNOUm1u7tOu\nclprq+oFF6iWl2fOk4jL08SJqhs2RB1tdNrbVS+7zOWpqChzniorVQ85RHXt2qijjU46rXrNNapl\nZarFxZ/NE7g8HXCA6iefRB1tdNJp1f/6L5enkpJt52mvvVSXL/e3bdxU673W2Jws7u3tqiee6E7E\nTEnr+iorU/3CF1S3bAm8u5yVTquefrr7APSSp3HjVDdvjjrq7S+dVj3/fG95Ki11J2RTU9RRR+Pb\n3/aWp5IS1V12yd8Pwh/+0FueiotVR45UbWz0vm2vxT0nu2Vuu819Rd60qfdlN2+G+nr42c+yHlbs\n3Hsv/PWv0Nzc+7KbN8Nbb8EPfpD9uOLmscfgwQe95amlxXVJ/L//l/244ubJJ+G3v/WWpy1bYOlS\n+MY3sh9X3Dz/PNx8s7c8tbbCxx/DueeGH0evD+sQkTJgNlCKuwD7iKpe122ZUuBe3LNTVwFfVdXF\nPW036MM6VF1f6Ecf+Vtv4EBYuRKKEnMJuWeqsMce8M47/tarrITGxvy6MHbAAfDqq/7WKS+HFStc\n/2m+OPJImD3b3zplZbBkCQwZkp2Y4uiLX4Tp0/1dnykthUWL3PWL3oT5sI4WYKKq7gvsB5wgIuO7\nLXMBsEZVdwN+CdzkYbuBPPMMrFvnf722NqirCz+euJo71/8HYKeHHgo3ljh78014++3el+uuoAD+\n8Ifw44mr998PNkBBBKZNCz+euFqxwl1Q9lPYO912W7ix9FrcO7p5mjp+LO54dQ99MnBPx/tHgKNF\nREKLsos5c7x93eluwwb3dSlfzJnjhl/51dQEs2aFH09czZnjCrVfGzfC3/8efjxxNXdusOHFmzbB\nzJnhxxNX8+YFGzbb0uI+FMLk6bAWkUIRmQ+sBGaq6kvdFhkBLAFQ1TZgHTA4w3amiEi9iNQ3NjYG\nCnjduuDDh1avDrZeLmpqcv15QaxdG24scdbU5L7VBRHkG2SuamoKft6tXx9uLHHW1BR8uGxTU+/L\n+OGpuKtqu6ruB4wEDhKRsd0WydRK/8w/UVXvVNVaVa2trq72Hy2u77ywMNCqedXv169f8BuUBg0K\nN5Y469cv+HWYgQPDjSXO+vUL9g0H3Fj4fNGXPIV9/cZXGKq6FngWOKHbn5YCowBEpAjoD2SlnTxh\nQrCLfZWVcMwxoYcTWxMmuP5Ov6qqYNKk0MOJrSOOCNZ9VVnpLpzli8MOcyNg/KqogJNPDj+euDr4\nYNfF4ld5OXzpS+HG0mtxF5FqERnQ8b4cOAbofgmqDugczHMa8Iz2NgwnoPHjoabG/3qVlXDcceHH\nE1f77ONGy/hVWAiTJ4cfT1ztuivU9jruILOvfjXcWOJshx3g6KP9Nxja27MzzC+uBg1y54/f1rsq\nTJkSbixeQqgBZonIa8DLuD73J0TkRhHp/EyeBgwWkUXAFcDV4Yb5KRG49lpIpbyvk0rBVVcF/7qU\nq370I/95uuKK/JuXx+/xVF4OF1+8feYHiZNrrvH3by4vd1OD5NNwUXC1xs9F1bIy+PKXIWBP9bZ5\nudMpG6++3KGaTqued563O8BSKdVTTnF3teajyy7zlqfyctXjjlNta4s64mh4vaOwvFz1sMPy845n\nVdWbbvJ+x/MBB6hu2hR1xNG49Vbvdzzvuae/O55J8vQDqq5Yf+c77iAqLf1s0kpK3N+mTMnfgqXq\nPgivu87loqzss3kqLna/P+us/C1Yqi5Pv/iFy0WmaS2KitzvTz01fwtWp1tvdXnKVLw68zRpUv5O\n0dDp7rtdLjLlqbDQ/f6oo1TXrfO33cQX905Ll6r+4AeqQ4a4ycIKClQHDnSF//33Q9lFIixfrnrD\nDarDhn2ap/79Xcv+nXeiji4+Vq5U/elPVWtqXI5EVKuqVC++WHXBgqiji49Vq1SnTnXzonTmqbLS\nzdHz+utRRxcfa9eq3nKL6k47uYLeOVHfWWepvvJKsG16Le69Tj+QLUGnH+hJ5zjcfOtb96uz/WB5\n6pnlyRtVd+4FHaKcL8LKk9fpBxI104qdhN6IBBsmmW8sT96IWGH3YnvnycqhMcYkkBV3Y4xJICvu\nxhiTQFbcjTEmgay4G2NMAllxN8aYBLLibowxCWTF3RhjEsiKuzHGJJAVd2OMSSAr7sYYk0BensQ0\nSkRmicgCEXlTRC7PsMwEEVknIvM7Xj/KTrjGGGO88DJxWBtwpaq+IiJVwDwRmamqb3Vb7jlVzaOn\nShpjTHz12nJX1eWq+krH+w3AAmBEtgMzxhgTnK8+dxEZDewPvJThz4eISIOIzBCRvUOIzRhjTECe\n53MXkUrgUeDbqrq+259fAXZS1SYRORH4MzAmwzamAFMAdtxxx8BBG2OM6ZmnlruIFOMK+32q+lj3\nv6vqelVt6ng/HSgWkSEZlrtTVWtVtbY69Ed9G2OM6eRltIwA04AFqnrzNpYZ3rEcInJQx3ZXhRmo\nMcYY77x0yxwGnA28LiLzO353DbAjgKreAZwGXCIibcAm4AyN6uGsxhhjei/uqvo80OOTJFX1VuDW\nsIIyxhjTN3aHqjHGJJAVd2OMSSAr7sYYk0BW3I0xJoGsuBtjTAJZcTfGmASy4m6MMQlkxd0YYxLI\nirsxxiSQFXdjjEkgK+7GGJNAnudzj7PmZqirgyVLQNvT7DCygMmToaoq6sjiZfNmePxxWLwY0mkY\nPhwmT4YBA6KOLF5aWuCvf4X33oO2Nhg2DE4+GQYPjjqyeGlthenTYdEi97662uXJZvPeWlsbPPkk\nLFzojq0hQ+CLX4SamizvWFUjeY0bN0776qOPVC+5RLWiQrWqSnVo4Se6lB10VMUqLS9XPf981cWL\n+7ybnPfxx6qXX65aWenyVFSkWljofi4rUz3rLNV33ok6yuitWqX63e+q9uvn8lRc7PJUUeHydPrp\nqm+9FXWU0Vu3TvUHP1AdMCBznk45RbWhIeooo9fUpHrDDaqDB2+dp1RKtbRU9aSTVOvr/W8XqFcP\nNTZni/sbb7ikFRW5fwWofp27VEHP5W4Fl8j+/VXnzu3TrnLav/6lOmyYO7A689T9VVjoDr7nnos6\n2ugsXqw6cqRqScm281RQ4ArYzJlRRxudZctUd9nFFadt5UnEFbC6uqijjc7Klap77uk+7HrKU3m5\n6gMP+Nt2oov70qWqgwZ9NlnPcoQq6CyO3Or3/frlZ8v0449dYRfZ9gHW9VVZ6T40882aNaqjRrni\n7SVPFRXBWly5bsMG1d1227pB1dMrlVJ9/vmoo97+Nm1SHTu25wZV11d5ub8Gg9fi7uVJTKNEZJaI\nLBCRN0Xk8gzLiIjcIiKLROQ1ETkgCz1I/3bttbC+21Ncq1jPwcwBYDxzqOLTBZqa4LvfzWZE8fTT\nn8Lq1e4Q8qKpCS67LLsxxdHNN8PKle46hBcbN8Ill2Q3pji64w746CPXh+xFczNcdFF2Y4qje+5x\n12taW70tv2kTXHCB9/PUK9FetigiNUCNqr4iIlXAPODLqvpWl2VOBP4TOBE4GPhfVT24p+3W1tZq\nfX2974DXr4fPDVtD9eYPt/r90fyd67mefmxgPVVcz/X8naP//ffSEvjr6ztSvftA3/vMRZs2wdCh\nrmD7UVYGb78NO+2UnbjiprXVXTBds8bfeuXlMG8e7LlnduKKm3QaRo6E5cv9rZdKwezZMG5cduKK\nG1XYdVd4/31/61VWwhNPwJFH9r6siMxT1dreluu15a6qy1X1lY73G4AFwIhui00G7u341jAHGNDx\noRC6Rx6B77T9nAb2Yw7jmc0RzOYIruMG+rEBgH5s4Dpu+Pff5jCeuVv24/2Lf56NkGLp8cdBenx+\nVmbpNNx9d/jxxNXMmd5bol21tsKdd4YfT1w9/7z/hgK4EVq33RZ+PHH16qvuW6BfGzfCrSE/y87X\nUEgRGQ3sD7zU7U8jgCVdfl7a8Tufn/O9e/99+K+2n9HIAK7lx/Rnfcbluv6+mXK+z09ZO+YqDgo7\noJhavNi13v3assW13PPF4sXBintbmxvali8WLw7WbZBO51+eCgv9r6cK77wTbiyeb2ISkUrgUeDb\nqtq9omZqI37mUBCRKSJSLyL1jY2N/iLt0N4OSgE/5/tM4FmWM5xNlGVcdhOlLKOGI/kHP+f7tKXz\n556tdDp4H57Xvuck6Eue2tvDjSXO7Hjypi//1rCPJ0/VTkSKcYX9PlV9LMMiS4FRXX4eCSzrvpCq\n3qmqtapaWx3wTocRI1x/J8DLHMQeLGQ5mXuAlrMDn+Nt6jmQ4uL86UcGd4NEZ578KCzMvzwVF/tf\nr6Ag//IUpEUKsOOO4cYSZzU1wT8ER43qfRk/vIyWEWAasEBVb97GYnXAOR2jZsYD61Q19C4ZgK98\nZevkbaGEmm30/gxnBS2UAu7APPPMbEQUT5MnB+tuKCmBc84JP564mjQpWIupvBzOPz/8eOLqqKOC\nXcOprIQLLww/nrgaP94NSvCrqgqmTAk3Fi8t98OAs4GJIjK/43WiiFwsIhd3LDMdeA9YBPwW+Ga4\nYX5q+HA47rhPD7RjmckWSgDYSIoWSthICoBWijmWmQAceKC7ip0vBgxwH4R+W1t77AH77JOdmOIo\nlYJzz/Xfeq+pgYN7HA+WLCUlcPHFUFrqb71+/WDixOzEFEeFhfDtb/sv8MXFbkqCUHkZDJ+NV19u\nYnr9dXcjCag+yH9oGrSJlP6K/9QKNuj/cqk2Ua5p0Ac4XVMp1TlzAu8uZ73zjrvz1MuNFJ03U8ya\nFXXU29+HH7o7mf3kafr0qKPe/lasyHzzYE83MT30UNRRb3+rVvm7eTCVUr3rLu/bJ8l3qKqqPvWU\nav/yFm0ipWvpp8fxt60SdjwzdC39tIkKfezBLX3aVy6bPdvdeeqlYP3hD1FHG525c92dzL2dkOXl\nqnfcEXW00WlocHPK9HY3b3m56tSpUUcbnbffdtOjFBb2Xtivv97fthNf3FVVX32hWR8dfonuWLpi\nq1t9i4rcwXX8fit02SmXqDY393lfuez111WPOsrNc9F17pTCQpen/fbLzxZ7dwsXqh5/vMtT17lT\nCgrcSbj33qp/+1vUUUbvvffcpFelpVvPndI5986YMap/+UvUUUZv6VLVU091OSov/zRPIi5PO++s\n+uCD/rfrtbj3eodqtgS9QzWTd9+Fu+5yU4+qwujRcN55+XP3oFeLF7sblBYudBdbR492/c351Mfu\nxdKl7nhasMCN+99xRzj7bDggq5Nq5J4VK9zx9PrrbirbESPgrLPc9a0gF1+TqrHR5em119y9JzU1\nbnDHoYcGy5PXO1QTUdyNMSZfhDb9gDHGmNxjxd0YYxLIirsxxiSQFXdjjEkgK+7GGJNAVtyNMSaB\nrLgbY0wCWXE3xpgEsuJujDEJZMXdGGMSyIq7McYkkBV3Y4xJIC+P2btLRFaKyBvb+PsEEVnX5SlN\nPwo/TGOMMX4UeVjm98CtwL09LPOcqob9kChjjDEB9dpyV9XZwOrtEIsxxpiQhNXnfoiINIjIDBHZ\ne1sLicgUEakXkfrGxsaQdm2MMaa7MIr7K8BOqrov8Gvgz9taUFXvVNVaVa2trq4OYdfGGGMy6XNx\nV9X1qtrU8X46UCwiQ/ocmTHGmMD6XNxFZLiIexKgiBzUsc1Vfd2uMcaY4HodLSMi/wdMAIaIyFLg\nOqAYQFXvAE4DLhGRNmATcIZG9WBWY4wxgIfirqpn9vL3W3FDJY0xxsSE3aFqjDEJZMXdGGMSyIq7\nMcYkkBV3Y4xJICvuxhiTQFbcjTEmgay4G2NMAllxN8aYBLLibowxCWTF3RhjEsiKuzHGJFDOF/ct\nW+CBB+Cww2DUKBg5Eg4+GO69FzZvjjq6+GhthUcfhSOOcHkaMQJqa2HaNGhujjq6+Ghvh7o6mDgR\ndtzR5Wn//eH222HDhqiji490GmbMgOOOg512cnnad1+45RZYty7q6OIjnYann4YTT4TRo2GHHWDs\nWJg6FVZn+/l2qhrJa9y4cdpXv/+9ar9+qlVVqrD1q7LSvX796z7vJuc9+KDqwIHbzlNFhepNN6mm\n01FHGq3HH1cdMiRznioqVMvLVa+/3vI0c6bqsGHu2MmUp7Iy1e99T7W9PepIozV7tuqIEZnzlEq5\nPF16qWprq7/tAvXqocbmbHH/7/92CeqetExJvOaaPu0qp91+uytKXvJ06aVRRxude+/1nqevfz1/\nC/yjj3rP02mn5W+BnzHDe32aNEm1rc37thNd3J980lviuibw4YcD7y5nPfectxOxa56mTYs66u2v\nvt7f8VRRofrLX0Yd9fb35pv+z7sbbog66u3v3XfdMeInT9/5jvfth1bcgbuAlcAb2/i7ALcAi4DX\ngAO87Lgvxf3AA70nrvM1Zkz+tbYmTvSfpxEj8i9PkyerivjL0+DB/lpbSXD22aqFhf7yVFWlunlz\n1JFvX5deqlpc7C9PqZTqhg3etu+1uHu5oPp74IQe/j4JGNPxmgLcHqTv36sFC+CNN/yvt2wZzJ0b\nfjxx9cEH8MIL/tdbtw6eeSb8eOLq44/hySfdKebHli3wxBPZiSmO1q6Fhx92F5z9UHUX8vNFczPc\nfbcbwOCHCNx3X7ix9FrcVXU20NN13cnAvR0fKnOAASJSE1aA3T3xhP8DDFzS//zn8OOJq+nToSDA\nWKimpvw6GZ96Cop6fR7ZZ23YAA89FH48cTVrFpSU+F+vqQnuvz/8eOLqn/+EwkL/623cGEFx92AE\nsKTLz0s7fpcVn3ziWk1+qbpWWr5Ysyb4UNAVK8KNJc7WrPHfyuqUb8dTkEYVQGNjuLHE2Zo1/r8F\ndlq1KtxYwijukuF3Gf95IjJFROpFpL4x4P/xigr3FSaIyspg6+WisrJgLVLIvzwFaWmBOxbzRVlZ\nsG+CAKlUuLHEWVlZ8PpUXh5uLGEU96XAqC4/jwSWZVpQVe9U1VpVra2urg60s7FjgxWfigp3M0q+\nGDs22NfoVMrd3JQvxo4NdjKWlcGBB4YfT1yNHRus5V5SAgcdFH48cbX33sF6FoqKwj/vwijudcA5\n4owH1qnq8hC2m9GXvhS8pfXVr4YbS5wdc0ywD8F0Gs45J/x44uqQQyBgO4OLLgo3ljj7/Odhl138\nr1dQAN/8ZvjxxNWuu8J++/lfr7gYLr883Fh6Le4i8n/Ai8AeIrJURC4QkYtF5OKORaYD7+GGQv4W\nyOr/yuJiuPRS13LyqqQEzj03v74eFhTAFVf4+6pXXAynnQYDBmQvrrgRgauu8ndsFBW52+6HDcte\nXHF09dX+uqIKC92H5047ZS+mOPKbJxH3zWjPPcONQzRo738f1dbWan19faB1m5th/HhYuLD3r0BF\nRW6OkFdegf79A+0uZ7W0wBe+AK+95t73pLAQampg/nwYPHj7xBcXbW1w7LEwZ07vF6ELCmDIEGho\ngOHDt098cZFOu2/Os2bBpk09LyviGgnz57vzL5+oul6CJ57oPU/g6tLLL8OYMd62LyLzVLXXTpyc\nnDgslYJ//AP22afnT8iKCpcrRRxSAAAMRElEQVSwF17Iv8IOUFrqJi2qre25i6aiwrWuXnwx/wo7\nuAbAX/8Khx/ee5522MHlKd8KO7gPtscecx+EPeUplYKhQ915l2+FHT4dsz55cs/1KZVy59vs2d4L\nuy9e7nTKxiuMicO2bFF94AHV/fd3t9n37+9eqZTqXnup3nNP/t0dl0lrq5sT5MAD3WRFXfO0226q\nv/udanNz1FFGr71dta5O9dBDt85TRYXq6NGqt93m/S7CJEunVadPV50w4bN5GjVK9ZZbVNetizrK\n6KXTqk8/rXrssVvnqbJStaZGdepU1dWr/W8Xj3eo5mS3TCbvvQdLl7qvjiNGZOmTMAEWL4YPP3R5\nGj4c9tgj+NCtJFuyxOWqrc31re+5p+Upk48+gvffd92jQ4e60SKWp89asQLefdd1jw4e7Hodgg4t\n9dotk5jibowx+SDRfe7GGGN6ZsXdGGMSyIq7McYkkBV3Y4xJICvuxhiTQFbcjTEmgay4G2NMAllx\nN8aYBLLibowxCWTF3RhjEsiKuzHGJJCn4i4iJ4jIQhFZJCJXZ/j710WkUUTmd7wuDD9UY4wxXvX6\nCGURKQR+AxyLe17qyyJSp6pvdVv0QVW9NAsxGmOM8clLy/0gYJGqvqeqW4AHgMnZDcsYY0xfeCnu\nI4AlXX5e2vG77r4iIq+JyCMiMiqU6IwxxgTipbhnmnq/+yTwjwOjVfXzwNPAPRk3JDJFROpFpL6x\nsdFfpMYYYzzzUtyXAl1b4iOBZV0XUNVVqtr5CObfAuMybUhV71TVWlWtra6uDhKvMcYYD7wU95eB\nMSKys4iUAGcAdV0XEJGaLj+eDCwIL0RjjDF+9TpaRlXbRORS4EmgELhLVd8UkRtxD2qtAy4TkZOB\nNmA18PUsxmyMMaYX9gxVY4zJIfYMVWOMyWNW3I0xJoGsuBtjTAJZcTfGmASy4m6MMQlkxd0YYxLI\nirsxxiSQFXdjjEkgK+7GGJNAVtyNMSaBrLgbY0wC9TpxWK5Yvhw++ghUYYcdYESmx4kYPv4YliyB\ndBpqamCUPVYlo8ZG+PBDaGuDYcNg9OioI4qnVavggw+gtRWqq2HnnUEyPQEiz61ZA4sXQ0sLDBkC\nu+66HfKkqpG8xo0bp33V1qZaV6d62GGqpaWq/fqp9u+vWlamWlur+vDDqq2tfd5NzmtvV50xQ3XC\nhM/mad99Ve+/X7WlJeooo5dOqz79tOqxx342T3vtpXrPPaqbNkUdZfTSadV//EP1pJO2zlN5ueqY\nMaq/+53qxo1RRxkPL7ygesopW+cplVLdeWfV229X3bDB/zZxs/H2WmNztrivW6d6+OGqlZXuX5Hp\nVVmput9+qp980qdd5bSmJtVjjuk9T5/7nOry5VFHG51Nm1RPPlm1oqLnPO2yi+qHH0YdbXS2bFE9\n4wyXJ5Ft52nkSNV334062ui0taleeGHPeaqoUB02THXBAn/b9lrcc7LPvaUFJk6El1+GpqZtL9fU\nBG+9BYcf3vNySdXaCpMmwfPP956nRYvgkENg7drtF19ctLfDl78MM2fCxo3bXq6pyXVBjB8Pn3yy\n/eKLC1U480z4y19cnrY1W3hTEyxbBgcf7P6bb1Th/PPh/vt7ztPGjbByJRx6qOuyCVtOFvdf/MIV\n7ZaW3pfdsgXefx+uvTb7ccXNbbfBvHmweXPvy7a1uRPxyiuzH1fc3HOP+wDctKn3ZdvbXX/8JZdk\nP664efhh+NvfvOUpnXb9zOefn/244mbGDHj0UWhu7n1ZVVi3Dr72tfDj8PSwDhE5Afhf3JOYfqeq\nP+/291LgXtyzU1cBX1XVxT1tM+jDOtraYPhwdyHHj8pK9ylZXu57lzlJFUaO9N9yKi+HFSugX7/s\nxBU3qrD77u6bix9lZe7C9JAh2Ykrjg44AF591d86ZWXwzjvuWMwXRxwBzz3nb53ycmhogDFjel82\ntId1iEgh8BtgErAXcKaI7NVtsQuANaq6G/BL4KbeQwxm+nTXGg/ikUfCjSXOZs2C9ev9r1dQAH/8\nY/jxxNXLL7uRVn6JwLRp4ccTV2+9BQsX+l9PFW6/Pfx44mrxYndM+dXWBrfcEm4sXrplDgIWqep7\nqroFeACY3G2ZycA9He8fAY4Wyc5An4aGYP3nTU2QT0/1a2gI9iG4cSO89FL48cRVQ8O2+0R7smkT\nvPBC+PHEVUMDFBb6X6+lJb/y9PrrUFLif73WVnjxxXBj8VLcRwBLuvy8tON3GZdR1TZgHTC4+4ZE\nZIqI1ItIfWNjY6CAN2wIdjJCsJZsrmpudq2BIDZsCDeWOGtudv3oQfR08TVpmptdP3oQ+ZanoPXJ\ny7UMP7wU90wt8O7he1kGVb1TVWtVtba6utpLfJ8xZAgUF/tfTwSGDg20y5w0cCCUlgZbN+D/mpw0\ncGCw4wnyq7994MBgLXeAwZ9p5iXXwIHBb04aODDcWLwU96VA1/sYRwLdL9P9exkRKQL6A6vDCLC7\nSZOgKMB9takUnHxy+PHE1fHHB2tBVFbCqaeGH09cHXOM+0rsV1UV/Md/hB9PXE2YEKybr7LSDZ/M\nF4ceGuwbc0VF+HnyUtxfBsaIyM4iUgKcAdR1W6YOOLfj/WnAM+plGE4A++wDe+zhf72hQ13i88Wu\nu8KBB/pfr7ISjj02/Hjiaocd4Oij/be2Cgvzq7EwaJC7F6AgwODp008PP564qqyEs87y3wBNp+Hs\ns8ONpdf/VR196JcCTwILgIdU9U0RuVFEOg/vacBgEVkEXAFcHW6YW/vxj11L3KtUCm68Mf/mvLj+\nev95uvbaYCdwLrv2Wn9DZFMpuPrq4N05uer733dDG71KpeCyy/ytkwRXXunvomoqBRdckIXhx15u\nY83Gq6/TD1x/vZujYVu3ine+UinVyy7r065y2tSp3vN0/vlu3pB8dMcdbm4UL3k6/fT8zdMf/+g9\nTyee6G7Dz0d/+pP3PB15pJvWwSuSPreMqupvfuPmZ8g0H0hFhUvuz36Wvydip7vvVq2qyjy/TCrl\nJsa69lrL0wMPuMmdMuWpvNzl6cor3URs+ayuTnXgQHdMbStP3/ymTdr31FOqgwdnzlNZmZtM7Lzz\n/E/a57W4e7pDNRuC3qHaXXMzPPCAuwGgc8rf4cPhW99yfV9VVSEEmwCbN7vbx3/1KzeVraobFXPJ\nJXDuudC/f9QRxsOWLfDYY3Dzze6GlPZ2NyrmwgvdV+dBg6KOMB5aW6GuzuVp0SL38+DBcN55cNFF\n+TXiqidtbe7Gy6lT3U1gra1uVMw558A3vuFqlV9e71DN+eJujDH5JLTpB4wxxuQeK+7GGJNAVtyN\nMSaBIutzF5FG4IMsbHoIkIePUvDN8uSN5ckby5M3YeRpJ1Xt9ZJ1ZMU9W0Sk3svFhnxnefLG8uSN\n5cmb7Zkn65YxxpgEsuJujDEJlMTifmfUAeQIy5M3lidvLE/ebLc8Ja7P3RhjTDJb7sYYk/cSVdxF\n5AQRWSgii0Qkq9MO5yoRuUtEVorIG1HHEmciMkpEZonIAhF5U0QujzqmOBKRMhGZKyINHXm6IeqY\n4kxECkXkVRF5Itv7SkxxF5FC4DfAJGAv4EwR2SvaqGLp98AJUQeRA9qAK1V1T2A88C07njJqASaq\n6r7AfsAJIjI+4pji7HLcczGyLjHFHTgIWKSq76nqFuABYHLEMcWOqs4mS49ATBJVXa6qr3S834A7\nIbs/GD7vdcxC29TxY3HHyy7kZSAiI4GTgN9tj/0lqbiPAJZ0+XkpdjKaEIjIaGB/4KVoI4mnjq6G\n+cBKYKaqWp4y+xXwPSC9PXaWpOKe6SF61oIwfSIilcCjwLdVdX3U8cSRqrar6n7ASOAgERkbdUxx\nIyJfBFaq6rzttc8kFfelwKguP48ElkUUi0kAESnGFfb7VPWxqOOJO1VdCzyLXdPJ5DDgZBFZjOsy\nnigif8zmDpNU3F8GxojIziJSApwB1EUck8lRIiK4B78vUNWbo44nrkSkWkQGdLwvB44B3o42qvhR\n1e+r6khVHY2rTc+o6lnZ3GdiiruqtgGXAk/iLn49pKpvRhtV/IjI/wEvAnuIyFIRuSDqmGLqMOBs\nXAtrfsfrxKiDiqEaYJaIvIZrYM1U1awP8zO9sztUjTEmgRLTcjfGGPMpK+7GGJNAVtyNMSaBrLgb\nY0wCWXE3xpgEsuJujDEJZMXdGGMSyIq7McYk0P8H3jzEnjUD6u0AAAAASUVORK5CYII=\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7fce55936150>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "simulation.run(1)\n",
    "simulation.show_beliefs()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "What do you think this call to `run` is doing? Look at the code in **`simulate.py`** to find out.\n",
    "\n",
    "Spend a few minutes looking at the `run` method and the methods it calls to get a sense for what's going on."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### What am I looking at?\n",
    "\n",
    "The red star shows the robot's true position. The blue circles indicate the strength of the robot's belief that it is at any particular location.\n",
    "\n",
    "Ideally we want the biggest blue circle to be at the same position as the red star."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "# We will provide you with the function below to help you look\n",
    "# at the raw numbers.\n",
    "\n",
    "def show_rounded_beliefs(beliefs):\n",
    "    for row in beliefs:\n",
    "        for belief in row:\n",
    "            print(\"{:0.3f}\".format(belief), end=\"  \")\n",
    "        print()\n",
    "        \n",
    "# The {:0.3f} notation is an example of \"string \n",
    "# formatting\" in Python. You can learn more about string \n",
    "# formatting at https://pyformat.info/"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "0.040  0.040  0.040  0.040  0.040  \n",
      "0.040  0.040  0.040  0.040  0.040  \n",
      "0.040  0.040  0.040  0.040  0.040  \n",
      "0.040  0.040  0.040  0.040  0.040  \n",
      "0.040  0.040  0.040  0.040  0.040  \n"
     ]
    }
   ],
   "source": [
    "show_rounded_beliefs(simulation.beliefs)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "_____"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Part 2: Implement a 2D sense function.\n",
    "As you can see, the robot's beliefs aren't changing. No matter how many times we call the simulation's sense method, nothing happens. The beliefs remain uniform."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Instructions\n",
    "1. Open `localizer.py` and complete the `sense` function.\n",
    "3. Run the code in the cell below to import the localizer module (or reload it) and then test your sense function.\n",
    "4. If the test passes, you've successfully implemented your first feature! Keep going with the project. If your tests don't pass (they likely won't the first few times you test), keep making modifications to the `sense` function until they do!"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Dimensionality error! Incorrect height\n"
     ]
    }
   ],
   "source": [
    "reload(localizer)\n",
    "def test_sense():\n",
    "    R = 'r'\n",
    "    _ = 'g'\n",
    "\n",
    "    simple_grid = [\n",
    "        [_,_,_],\n",
    "        [_,R,_],\n",
    "        [_,_,_]\n",
    "    ]\n",
    "\n",
    "    p = 1.0 / 9\n",
    "    initial_beliefs = [\n",
    "        [p,p,p],\n",
    "        [p,p,p],\n",
    "        [p,p,p]\n",
    "    ]\n",
    "\n",
    "    observation = R\n",
    "\n",
    "    expected_beliefs_after = [\n",
    "        [1/11, 1/11, 1/11],\n",
    "        [1/11, 3/11, 1/11],\n",
    "        [1/11, 1/11, 1/11]\n",
    "    ]\n",
    "\n",
    "    p_hit  = 3.0\n",
    "    p_miss = 1.0\n",
    "    beliefs_after_sensing = localizer.sense(\n",
    "        observation, simple_grid, initial_beliefs, p_hit, p_miss)\n",
    "\n",
    "    if helpers.close_enough(beliefs_after_sensing, expected_beliefs_after):\n",
    "        print(\"Tests pass! Your sense function is working as expected\")\n",
    "        return\n",
    "    \n",
    "    elif not isinstance(beliefs_after_sensing, list):\n",
    "        print(\"Your sense function doesn't return a list!\")\n",
    "        return\n",
    "        \n",
    "    elif len(beliefs_after_sensing) != len(expected_beliefs_after):\n",
    "        print(\"Dimensionality error! Incorrect height\")\n",
    "        return\n",
    "              \n",
    "    elif len(beliefs_after_sensing[0] ) != len(expected_beliefs_after[0]):\n",
    "        print(\"Dimensionality Error! Incorrect width\")\n",
    "        return\n",
    "    \n",
    "    elif beliefs_after_sensing == initial_beliefs:\n",
    "        print(\"Your code returns the initial beliefs.\")\n",
    "        return\n",
    "        \n",
    "    total_probability = 0.0\n",
    "    for row in beliefs_after_sensing:\n",
    "        for p in row:\n",
    "            total_probability += p\n",
    "    if abs(total_probability-1.0) > 0.001:\n",
    "        \n",
    "        print(\"Your beliefs appear to not be normalized\")\n",
    "        return\n",
    "    \n",
    "    print(\"Something isn't quite right with your sense function\")\n",
    "\n",
    "test_sense()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Integration Testing\n",
    "Before we call this \"complete\" we should perform an **integration test**. We've verified that the sense function works on it's own, but does the localizer work overall?\n",
    "\n",
    "Let's perform an integration test. First you you should execute the code in the cell below to prepare the simulation environment."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "from simulate import Simulation\n",
    "import simulate as sim\n",
    "import helpers\n",
    "reload(localizer)\n",
    "reload(sim)\n",
    "reload(helpers)\n",
    "\n",
    "R = 'r'\n",
    "G = 'g'\n",
    "grid = [\n",
    "    [R,G,G,G,R,R,R],\n",
    "    [G,G,R,G,R,G,R],\n",
    "    [G,R,G,G,G,G,R],\n",
    "    [R,R,G,R,G,G,G],\n",
    "    [R,G,R,G,R,R,R],\n",
    "    [G,R,R,R,G,R,G],\n",
    "    [R,R,R,G,R,G,G],\n",
    "]\n",
    "\n",
    "# Use small value for blur. This parameter is used to represent\n",
    "# the uncertainty in MOTION, not in sensing. We want this test\n",
    "# to focus on sensing functionality\n",
    "blur = 0.1\n",
    "p_hit = 100.0\n",
    "simulation = sim.Simulation(grid, blur, p_hit)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "NOTE! The robot doesn't have a working sense function at this point.\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAW4AAAD8CAYAAABXe05zAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4wLCBo\ndHRwOi8vbWF0cGxvdGxpYi5vcmcvpW3flQAAIABJREFUeJzt3Xd4HdWZP/DvkXR1myQ3CReMbVoA\nB+MigSGYsnYAx0AoG0wnwCaGLLC0YGDDbzeQbEIJbMguJUAoCd0YCBg2lADBoQTL2JhiYwyxjXGT\nLRdJVr/f3x+vFLlId85caXRmzPt5nnmscjT3PfPOeWfuzBlfQxJKKaWiI891AEoppfzRwq2UUhGj\nhVsppSJGC7dSSkWMFm6llIoYLdxKKRUxWriVUipitHArpVTEaOFWSqmIKQhipaWlpRwxYkQQq1ZK\nqZ3S3Llz15Ess2kbSOEeMWIEKisrg1i1UkrtlIwxy2zb6qUSpZSKGC3cSikVMVq4lVIqYgK5xp2L\nr74CnngCWLkSKC0FTjkF2HNP11H519oKvPgi8Pbb8v2BBwLHHw/EYm7jysWyZcCTTwKrVwODBgFT\npwLDh7uOyr/mZuD554E5c+T7Qw4Bjj0WyM93G1cuPv8cmDEDWLcOGDIEOPVUYNddXUflX2Mj8Oyz\nwPvvSx4OPxw4+mggL4KnkosWAU8/DVRXA7vtBpx2GjBwYMAvSrLHl/LyctqqqyOnTiUTCTIeJwEy\nFpPvv/1tcv1661U59+KLZGkpWVws/QDk6379yBkzXEdnb/Nm8vjjJQeFhdKPwkL5/rjj5PdRMWOG\nbP/tc1JaSr7wguvo7K1fTx51lOQgFpN+xOPy/dSpMo6i4sEHyZISsqioIydFReTAgeRrr7mOzt6a\nNeSECWQySRYUSD8SCVnOPZdsbPS3PgCVtKyxdo2AvgCeArAIwEIAh2Rrb1u4m5ul44lERwK3XgoL\nyW98g6yp8bcBXHjxRUlgZ/0A5HdPPOE6Sm/19eSYMR0H0e2XeFx+X1/vOlJvTz7pnZMXX3Qdpbea\nGnKffToOotsviQT5rW+RTU2uI/X2u9+RqVTXOUmlyDfecB2lt+pqcvjwjoNoZ/vW0UeTra326wyi\ncD8E4AdtXxcC6JutvW3hfvRRMp3uOontG+Cmm+w770JLC7nLLtn7AZB9+vg/Cve2u+7KPrDaB9dd\nd7mONLvGRjmr88rJLrtI/sLsppuyH4AAGUePPuo60uxqa733LYDcfXcyk3EdbXbXXNP1yc3WOZk1\ny36dPVq4AZQA+DsAY7tS28I9erR3EgF5C+XnyNXbXnhh27fiXS3FxeTjj7uONrvdd7fLye67u440\nu8ces89JmC+ZZDKy/9vkZPRo19Fm99vfep+otRe82bNdR9u1piY5CbPJyWGH2a/XT+G2uRWwB4Aq\nAA8YY+YZY+4zxqR74vr6xx/btdu4EdiwoSdeMRhz5gC1td7tamo6blqGUXMzsHSpXdulS4GmpiCj\n6Z533pHt7aW2tuOmZRhVV9vv+7bjyZU33wTq6rzbNTeHOydffSUx2pg3L5gYbAp3AYBxAO4iORZA\nHYBrtm9kjJlmjKk0xlRWVVX1cJgqbIxxHUHPkDeVSvnjev+3KdwrAKwg+be275+CFPJtkLyHZAXJ\nirIyq8ft8c1v2gXZpw/Qr59dWxcOPBAoKvJuV1wsU9HCKhYDbP+LmeHDgcLCQMPplkMOke3tpbgY\nqKgIPp5c9e9vv++PHBlsLN11+OFA2uK9eiwW7pzsuitQYDmReuzYYGLwLNwkVwP40hizT9uPJgH4\npCdefPp070Qmk8AVV4R7fucxxwCJhHc7Y4CTTw4+nu646ioglcreJpWSdmF20kl2Z0WJBDB5cvDx\n5MoY4PLLZRxkk04DV1/dOzHl6owz7N7hlJUBEyYEH0+uYjHgwguBeDx7u6IiqXGBsLkQDmAMgEoA\nCwA8C6Bftva2NyebmmQaU7bpgHvvvfNMBwz7jUnSbjrg6NHRmA74xBPeOQnzjcl2NTUyLXZnmA54\n3307z3TAYcM65m93tm8ddZTj6YB+F78P4Hzve50/gDNpUrQewJk1ixwwYNvZDEVFZN++0ZjD3W7T\nJvLYYzt/AOfYY+X3UfHkk7L9t8/JgAH+pmq5tn69jIftH8CJx2X8ROkBnAcekHx09gDOn//sOjp7\nq1fLAXP7B3DicfKcc8iGBn/r81O4DQO4O1NRUUG//63rihXyyPtXX3U88r733j0eWuBaWoAXXpBZ\nDaRc/z7hhGg+8r506Y6PvEfxv1lvbgb++EegfZdsf+Td9jplmHz2GfDUU0BVlVxrPfVUYOhQ11H5\n19goj4nPny+PvE+YIJeswnxJtCsLFwIzZ8rsn/ZH3gcN8r8eY8xcklZX90NTuJVS6uvMT+GO4PFN\nKaW+3rRwK6VUxGjhVkqpiNHCrZRSEaOFWymlIkYLt1JKRYwWbqWUihgt3EopFTFauJVSKmK0cCul\nVMRo4VZKqYjRwq2UUhGjhVsppSJGC7dSSkWMFm6llIoYLdxKKRUxWriVUipitHArpVTEaOFWSqmI\n0cKtlFIRE4rCXVcH3HMPMHasfHL1yJHAbbfJpyZHzdy5wJlnAsOGyfK973V84nuUbN4M3HEHcMAB\nkpNRo4D//V/5eZSQwLvvAqec0pGTM86QPEXNhg0yLkaOlJyMHSvjpq7OdWT+kMCbbwInnCCfij58\nOHDeecCCBa4j82/dOuDGG4H99pOcVFQADz4I1NcH/MIke3wpLy+nrQ8/JEtLyXSalJTKkkqRxcXk\nm29ar8qp1lbyX/5F4s7L6+hHXp70bepUsrnZdZR25swh+/bdMSfptPx8zhzXEdppbiZPPVXi3j4n\nqRR5/vmStyiYPVvGQyq1Y05KS2UcRUFDA3nssRK3MR39yM8nk0nyssvITMZ1lHZefln6kUxum5Oi\nInLQIHLxYn/rA1BJyxpr1whYCuBDAPNtVm5buNeuJfv337bT2y/ptP8N4MLVV+84qLY/EF14oeso\nvX35JVlSkj0nJSXk8uWuI/V24YXeOZk+3XWU3j77TIpBtpz07y/jKezOOGPHQrd9Tn7xC9dReluw\nIPu+ZQw5cCC5aZP9OoMq3KW2K7Ut3DfcQCYS2XfI/Hw5MwqzTZu8+wFImzVrXEeb3RVXkIWF2ftR\nWCjtwmzNGjIe985JMulvcLlw/vkyDrz2rRtucB1pdn//u904KSoi6+tdR5vdySdv+46hq4PQb35j\nv04/hdvpNe477gAaGrK3aW0FHn0UaG7unZhyMWMGkJ9v1/bhh4ONpTtI4N57gaam7O2amqQdQ3zd\n/uGHgTyLvTsvT/IXVs3Nsv+3tmZv19Ag4ynMHngAyGS82xkDPPdc8PHkqrYWeOEF7/1/yxbg9tuD\nicG2cBPAy8aYucaYaT314lVV9m2rq3vqVXve0qV2N4gaGoAlSwIPJ2eNjbKz2diyxfug69Lnn9vd\nIKqrk/yFlZ8b9H7Gkwuffup9UgDIfrVsWfDx5GrNGiAWs2u7alUwMRRYtjuU5EpjzC4AXjHGLCL5\n5tYN2gr6NAAYNmyY1UoLC+0Gf0sLkEpZRupAURFQUCBxZmMMUFLSOzHlorDQ/iyaBOLxYOPpjuJi\n2d5e/SkokPyFVSrlvV+1KywMNpbust33CwqAdDrYWLojlbK/ApBIBBOD1Rk3yZVt/64F8AyAgzpp\ncw/JCpIVZWVlVi8+ZYoMLi9jxshADKvjjrM7AqdSwEknBR9PrvLygCOPtGt75JF2lyJcOfFEu4N9\nLCb5C6uiItn/vRgj4ynMTjnF7iBJAsceG3w8uRo0SKYweikokP0wCJ5DzxiTNsYUt38N4GgAH/XE\ni0+fDiST2duk08C11/bEqwXnm9+Uec7ZrnPn5cmc1YN2OOSFy7XXep/tpNPANdf0Tjy5Gj9etne2\ng0t+PrD//pK/MPv3f/fOSTIJXHVV78STq0mTgH79sp+sxWLAoYfaFUZXjLHLSSwGXH55QEF43b0E\nsAeAD9qWjwH8xOtv/Mzj/tnPup5Wk0qRP/xhNOZ1Ll9O7rILWVDQ+cyY/v2jMa2RJK+6KntOfvxj\n1xHa+ewz2e6dzcgoKJB8RWFaYyZDTpu247z6rXPy85+7jtLOggUynXTrefXtSyxGDh1Krl7tOkpv\nmQx52mld5ySZ9DejhPQ3q8Sqkd/FT+EmyZkzyZEjZQfs00c2xvDh5H33RaNot1u1ijzvPElaSYn0\nJZGQuavLlrmOzp9HHiH32qsjJ6mUfP/II64j82f5cvLMMyUPffpIXpJJ8txzyZUrXUdnL5OR8TB8\nuIyP9pyMHCnjJ0qWLCG/971tc5JKkRdcQFZVuY7OXmsr+T//Qw4ZIlMY+/SRfWvMGHLWLP/r81O4\nDQOY01VRUcHKykrff/fZZ3LHtn9/eYTU5vp3GNXUAJ98IsfeffcF+vZ1HVFuSJkJsG4dUFoK7LNP\ndHOycaP0BZBHxsN8zyQbEli4UGZZDRwI7L2364hyV10NLF4sl7O++c1w35DMhgQ++kj+O4jBg4E9\n9shtPcaYuSQrrNqGqXArpdTXlZ/CHeJ5AUoppTqjhVsppSJGC7dSSkWMFm6llIoYLdxKKRUxWriV\nUipitHArpVTEaOFWSqmI0cKtlFIRo4VbKaUiRgu3UkpFjBZupZSKGC3cSikVMVq4lVIqYrRwK6VU\nxGjhVkqpiNHCrZRSEaOFW9nJZFxHoJRqo4VbeVu/HthtN/mQQKWUc6Eq3O0fHrp2retIuqe5Gfji\nC1mamlxH0z3r1gGr73kOWLkSeO451+HkrKmpIyfNza6j6Z61a2WcRP042tgIfP458Pe/Ay0trqPp\nnjVrJCcbNvTO64WicL/xBnDkkfIJyeXlcnJXURG9OrFpE3DttUBZGXDAAbKUlgJXXikFMEpefhk4\n5BBg112BJdc9CAB4/98exEsvuY3Lr/XrgR//WPLQnpOyMsnTxo2uo/Pn+edlXOy2m4yTwYNl3Lzx\nhuvI/Fm9GrjkEmDAAGD0aGD//eUT63/6U6C21nV09khgxgzZp4YPl5wMGgQcdRTwzjuBvzh7fCkv\nL6etu+8mUylSNsO2SypFXned9aqcWruW3H13Mh7fsR+FheSQIeSKFa6jtHPLLR05KcYm1qOQBFiP\nOAcmN/GWW1xHaGfFCtnuhYU75iQeJ0eMkLxFwf/7f9nHyd13u47Qzuefk6WlZCy2Yz8SCXLffcmN\nG11H6S2TIf/t37rOSTJJPvKIv3UCqKRljbUuxgDyAcwDMMurrW3h/uCDrjvevqTT5J/+5G8DuPBP\n/9T5zti+FBSQFRWuo/RQXc15D87n+MR8HgBZLset3IRiEuAmFPNy3Mrxifl8/4H55Py2pbradeSd\nKi+X7d5VTmIxyVvYvfSS9zhJpSQVYZbJkPvsQ+bldd2PeJw86STXkXqbMUNqU7acJJPkZ5/ZrzOo\nwn0FgEd7snCfdRaZn5+98wB52GH2nXdhyRI5W/DqR+gH1/TpJMAtSHAjSv6xbN2J9p/VFpR0dHr6\ndNeR72DePO9i136Wt2SJ62izO+ww737k58t4CrPZs72LXXtOVq50HW12o0d79yMWIy++2H6dfgq3\n1TVuY8xQAMcCuK8nL9M88wzQ2urd7p13gLq6nnzlnvXMM3az5Rob5ZpYWGX+65f4ifkFCIM+2PyP\nZWvtP0u3bAaNAX7xC+CXv3QUcddmzpTt7SWTkfyF1ZYtdtdLW1vD3Q8AeOwx6Y+XvLxw399auxZY\ntMi7XXMz8PjjwcRge3Py1wCmA+iyPBljphljKo0xlVVVVVYrra+3e/GCgnDftNi40W72SGur3CwL\nq8bmPNyUdy2OxBtYhUGoR6LTdvWIYxUGo/Glv8hdvrxQ3OPexvr1dicFTU1yUzmsampk/7dhO55c\nWbdOzkW9NDUBmzd7t3Nl82YgFrNra3OgyoXniDPGHAdgLcm52dqRvIdkBcmKsrIyqxfv188uyEwG\n6NvXrq0LQ4cCqZR3u3gcGDEi8HBylkgAhYXAHByEffApVmFwp+1WYQjGJBYhPuHAXo7Q3vDhsr29\npFLAkCHBx5Orvn3tn32yHU+u7L67XcFLJsOdk7Iy+2m+paXBxGBzqnQogO8aY5YCeBzARGPMwz3x\n4j/4gffgyssDTj7ZbhC6csopdmd3xgBnnhl8PLkyBjjnHDnDa0IhBmNVp+0GYTVOOSsOY3o5QB/O\nPBNW8bW2AlOnBh9PruJx2f+93tTE4zKewuy88+zePbS2AiecEHw8uerTB5g40btdMgn8678GFITt\nxXC5do4j0YM3J1esIIuLvW/oLVhgf4HflQsuyH4zLJEgTzvNdZTeliyRG0jH4bl/3JisRYoNKGQt\npIObUMKvfvu861A9nXaa3NnPtm9dcIHrKL19+KH3jdbi4mhMN50ypfMps1vnJIT3unfw7rvZ9y2A\n7NePXLfOfp0IYlYJAyjcJPn222RJyY6zMgoLJYnPh78+kCSbmsjvfrfzu+ZFReTEiWR9veso7fz5\nz+TM/FOYaSvav8YlTKOG/5t3MeuQZAYgp051Haan+npy0iTZ/tvnJJ0mjz9e8hYFs2ZJzNvPSU8m\nZfy8/bbrCO3U1JAHH9x1Tk4/nWxpcR2lnSeekBq1/TTgZFKK9vvv+1tfYIXbdvFTuElyzRryZz8j\nhw6VDTFokBx1ly3z13HXMhmZcz5xogym4mJywgTyueeiszOSJBsb2ZpMsb6whKf3/xNTKXmQ5brr\nyHUP/590Lp2ORNVraZGD/4QJEnZJieTn//5P8hUly5bJuBg0SMbJ0KEybtascR2ZP83N5MyZ5Pjx\nMkZKSsjJk+WEIWo5+fxz8tJLybIyGRLDh5M33eTvTLudn8JtpH3PqqioYGVlZY+vV/WS+np5Tv8/\n/1OeRd7emjXA9dcDt94qF/KUUt1mjJlLssKqrRZupZRyz0/hDt8EXKWUUllp4VZKqYjRwq2UUhGj\nhVsppSJGC7dSSkWMFm6llIoYLdxKKRUxWriVUipitHArpVTEaOFWSqmI0cKtlFIRo4VbKaUiRgu3\nUkpFjBZupZSKGC3cSikVMVq4lVIqYrRwK6VUxGjhVkqpiNHCrZRSEVPgOoB2mzYBzz4LrF0L9OsH\nfPe7wC67uI7KPxJ45x3gb3+T78eNAw4/HDDGbVy5qK4G/vhHYN06YMAA4MQTgf79XUflHwm8+Sbw\n/vvy/fjxwCGHRDMna9cCzz0HbNgg4+PEE4E+fVxH5V8mA7z2GrBgAZCXBxx6KHDgga6jys2qVcDz\nz0sNGzxYclJUFPCL2n4cvJ+lvLzc+iPpm5rIiy4iEwmyqIgsKJCPuU8kyNNPJ2tr/X/MvSt//Su5\n554Sf2GhLEVF5G67ka++6jo6e/X15Pnnd56T886T30fFq6/K9i8q6shJOk3usYfkKypqa8kzzpAc\npNOSk6Ii+f6ii2QcRcXzz5ODB3fkJB4nUyly333JykrX0dnbuJE86aSO+NtzkkySV11FtrT4Wx+A\nSlrWWO8GQALAewA+APAxgOu9/sa2cLe2ksceK52W86Jtl0SCLC8nGxr8bQAXZs/uuh+A/O5Pf3Id\npbfmZvKII2Tn66wfyaT8vrnZdaTeXnrJOyezZ7uO0ltDA1lRIeOhq35MmSLjKeyefrrrfQuQg1IU\nindtLTlypBTtrnJy6qlkJmO/zp4u3AZAUdvXMQB/A3Bwtr+xLdxPPy2J6iqJ7Rvg9tvtO+9CJiNn\nddn6AZADBoS/4D34oHdO0mnygQdcR5pdc7Nsb6+cDB3qb3C58JvfZD8Atedk5kzXkWZXX08WF3vn\nZN99XUfq7YYbuj6Qbp2TV16xX6efwu15c7JtnbVt38baFvbEZZqbbgLq6rK32bIF+NWvZFOE1euv\nyzVHL01NwKxZwcfTHTY5qauTdmE2a5Zsby8bN8q11rAigVtukXGQTRRyMmOG3ThevhyYMyf4eHLV\n2grcfjvQ0JC9XV0dcPPNwcRgNavEGJNvjJkPYC2AV0j+rSdevP1mkZc1a+wKoytvveVd7ACgpgb4\ny1+CjydXzc3AokV2bT/9VNqH1V/+ItvbS10d8PbbwceTq40bZf+3MW9esLF016uvArW13u1aW2VM\nhdWKFUB9vV3bd98NJgarwk2yleQYAEMBHGSM2X/7NsaYacaYSmNMZVVVldWLZzJ2QRpj39aF1lb7\ndwQtLcHG0h2k/UyLsOfEdjuTkr+wam21z0mY8wH4y0mY+5LJ2OckqCsFvuZxk9wI4A0Akzv53T0k\nK0hWlJWVWa1vr73sXjeVCvc0tNGjgeJi73ZFRUBFRfDx5KqwUKYz2Rg0CIjHg42nOyoq7KZkFRcD\nBxwQfDy56t8fSKft2tqOJ1cOPljGspfCwnDnZMgQ+8K9337BxOBZuI0xZcaYvm1fJwF8G4DlG+rs\npk/33injceCSS2SuZ1gddxyQn2/XdurUYGPpriuuAJLJ7G2SSeDKK3snnlydcord2U5BAXD88cHH\nk6u8PODii70Pkum0jKcwO/tsuzPp4mJg4sTg48lVPA6cfz4Qi2VvF2hOvO5eAjgAwDwACwB8BOA/\nvP7GdlZJfT05apTM5ezsrmxBAbnrrmR1tf2dWVceeyz7NKdUirz7btdRequpIffai4zFOu9HLCZz\n1WtqXEfq7e67s8/GSCbJRx91HaW36mqZ/VJQ0Hk/CgvJ/fePxvz6m2/2zsnzz7uO0tvq1WRZGZmX\n13k/4nHyoIP8zSJDT04HzGXx8wDOhg3kkUdKwvLzJaK8PJlKM3o0uWKFfcdde+ghibuoaNspQckk\neccdrqOzt3YtOX68DLD2HTMvT74fP55cs8Z1hPbuvFO2/9ZTHIuKpC8PPeQ6OnsrVsh4SKc7cpKf\n3zGvPgonN6RMvbz5ZplKt3UBLyqSJexTGrf2xRcydTGdJo3pOLFJJsnJk/2f3Pgp3IYBXD2vqKhg\nZWWlr7/56CPgoYdkKtAuu8jbqgMPjN5jyVu2AI8/DrzxhrwtPPRQ4Kyz7K6Bh828ecDvfw+sXCnX\n9c4+Wx7hj5qaGuDhh2X2iDHAEUcAp59ud701TEiZJvfwwzLTZNgw4PvfB/bfYapA+G3cKPvWe+/J\nZcZJk+QyYiLhOjJ/SJkB89hjwPr1wIgRwHnnAfvs439dxpi5JK3ugoWmcCul1NeZn8Id4lt+Siml\nOqOFWymlIkYLt1JKRYwWbqWUihgt3EopFTFauJVSKmK0cCulVMRo4VZKqYjRwq2UUhGjhVsppSJG\nC7dSSkWMFm6llIoYLdxKKRUxWriVUipitHArpVTEaOFWSqmI0cKtlFIRo4VbKaUiRgu3UkpFjBZu\npZSKmFAU7pYW4OmngaOPBkaNkk/h/sMfgIYG15H59/e/A5dfDowdC4wZA1x0EbB4seuo/GtqAp54\nQj59e//9gYkT5dPrm5pcR+bf4sWShzFjJC+XXy55ipqGBhkXRxwh4+Too2XctLS4jsy/jz8Gpk0D\nRo8Gxo0Drr4aWL7cdVT+bdkC3H8/cNhhkpPvfAeYNQtobQ34hUn2+FJeXk5bX3xBDh9OFhWR8mH3\nshQVkaWl5Lx51qtyKpMhr72WTCTIwsKOfsRiZDJJXnwx2drqOko7CxeSgwaRxcXb5qS4WH6+cKHr\nCO20tpKXXCLbPxbr6EdhoeTpmmskb1Ewf76Mh+3HSXGxjJ8vvnAdoZ2WFvLccyUn+fkd/YjHJSe/\n+IXrCO29+y7Zt2/nOdl7b3LFCn/rA1BJyxrr3QDYDcDrABYC+BjApV5/Y1u4N2wgBw8m8/K27fjW\nS58+5PLl/jaAC7/8JZlKdd2PVEoKe9itWUMOGEAa03k/jCH795d2YXfNNd45+eUvXUfpbflyGQdd\n9SMvT8bRhg2uI/X2ox955+TOO11H6W3x4h0L9tZLfr4cUOvq7NfZ04V7MIBxbV8XA1gMYGS2v7Et\n3LfcIkferjrffsZ68cX2nXehtjb7zti+JBJkdbXraLP7yU/k7CdbP+JxaRdm1dWyvb1ykkpJ/sLs\n4ou3fcfQ2ZJMyngKs6++8t63ADmLbWpyHW12Z5657TuGzpZ0mrz3Xvt1+incnte4Sa4i+X7b1zVt\nZ9679sRlmttvB+rrs7dpbgYeeCDc1/GefhrIs7hbkJcHPPJI8PF0x513Ao2N2ds0Nkq7MHv0Ufuc\nzJwZfDy5ammR/b+5OXu7+noZT2F2//2AMd7tWluBF14IPp5cbdki+4zXdey6OuDWW4OJwdfNSWPM\nCABjAfytJ1585Uq7dq2tQHV1T7xiMD77DKit9W63ZQuwcGHw8eSqoQHYtMmu7ebN4b55/Mknsr29\n1NYCS5YEH0+uqqvtb3TZjidXPvrIbp9paAh3TlatAvLz7doGdcPVunAbY4oAzARwGcnNnfx+mjGm\n0hhTWVVVZbXOWMzutVtbgXjcNtLel0zaJdIYIJUKPp5c2eYDkDeDftr3tnTa7uwuPx9IJIKPJ1eJ\nhP27zTDnA5Cc2IhCTmwPpoWFwcRgVbiNMTFI0X6E5NOdtSF5D8kKkhVlZWVWLz5xol2Q++wD9Olj\n19aFY46xO7Ck08BxxwUfT67y84Hx4+3ajh9vf9bhwpQpdoUiHgcmTw4+nlyVlAD77mvX1nY8uXLi\niUBxsXc7Y2SqY1gNGQIMHOjdLi9PpgcGwbNwG2MMgN8BWEjytp588enTvc9A02mZ4xlm48YBe+zh\nfYZXWgocfnjvxJSra67xLnhRyMkRR8j2zsYYYPfdJX9hZpOTVAq46qreiSdXU6bIu9Ns8vNlvv03\nvtE7MeXCGNnWXrUrkQCuvDKgILzuXgKYAIAAFgCY37ZMyfY3fuZxX3pp1zMyUinyn/85GvOfP/1U\n7oZ3NrXRGJnbOX++6yi9ZTLkeedlz8m550Zj/vP8+bLdO5vamJcn+Vq0yHWU3lpbZRxky8mll7qO\n0s4778hsi66m0JWWkkuXuo7SW0sLecwxXc+KS6XI//gPf+tET04HzGXxU7gzGfKuu2QeanExWVIi\n//brR/7Xf0WjaLf7/HNyyhSZ8lRSInNv43Fy0iTyk09cR2cvkyFvu00G0dY5KS2Vn0ehaLf75BPZ\n/vG45KOkRL7+znfIJUtcR2d1RRk1AAARA0lEQVSvtVXGQ79+2+Zk0CAZP1HKyfz55GGHyXTNkhJZ\nEgnyxBOj8cxGu+Zm8rrrOvarkhKZ2z10KPngg/7X56dwG2nfsyoqKlhZWenrbzIZ4L33gDVrgH79\ngG99Cygo6PHQesWqVcC8eXLsPeAAYLfdXEeUm9ZW4N13gXXr5LLDwQeH+7p2Nl9+CXz4oXw9diww\neLDbeHLV0gK8847MNhk4EDjoILtpj2G0dKk8+p6XB5SXA7vs4jqi3DQ1AW+/LTOtBg8GKirsboxv\nzxgzl2SFVduwFG6llPo681O4I3qsVkqpry8t3EopFTFauJVSKmK0cCulVMRo4VZKqYjRwq2UUhGj\nhVsppSJGC7dSSkWMFm6llIoYLdxKKRUxWriVUipitHArpVTEaOFWSqmI0cKtlFIRo4VbKaUiRgu3\nUkpFjBZupZSKGC3cSikVMVq4lVIqYkJVuJubgfXrgcZG15F0Dwls2iRLAB/p2auamiQnTU2uI+me\nnSknjY2Sk+Zm15F0Dwls3Cgfsqs58ScUhXv+fODUU4GiImDXXeXfKVOAv/7VdWT+1NcDv/kNMGwY\nUFYmy5AhwM03A7W1rqPzZ84c4MQTt83JCScA773nOjJ/amuBW26RPLTnZNgw4PbbJV9R8tZbMi62\nzsnUqTJ+omTjRuDnP5dPdR84EBgwANhjD+Cuu6J3gvDaa8CkSR05KS4GzjoL+OSTgF+YZI8v5eXl\ntPXkk2QqReblkXLc7VhSKfLXv7ZelVObNpGjRknM2/cjmST33ptct851lHbuv1/6Ycy2/TBGfn7/\n/a4jtLNunWz3ZLLzfWvUKMlbFPz6153vW3l58vMnn3QdoZ2VK8lhw8hEovOcjB9P1tW5jtLOz37W\neU7y8+XnL7zgb30AKmlZY70bAPcDWAvgI9uV2hbuRYs67/j2yXzzTX8bwIXjjyfj8a77UVhIHnGE\n6yi9zZ1rl5PKSteRejviCDIW67of8bjkLezefNMuJwsXuo7U29ixZEFB1/1IJMizz3YdpbcXX/TO\nSTpNLltmv04/hdvmUsmDACYHcbZ/223eb422bJG3VWH25ZfAK69kvzbf1CSXGRYt6r24cnHjjUBD\nQ/Y2DQ1y+SfMFi6U7Z3tmmNjo+Ttyy97L65c/PznMg6yaWoC/vu/eyeeXM2ZAyxeDLS0dN2moQGY\nMQNYt6734srF9dd756S5GbjjjmBe37Nwk3wTQHUQL/7YY9mT2O7118N9PXLmTLubK83NwOOPBx9P\nrkjg2WeBTCZ7u0wGeOaZcN9QeuIJuxtFJPDUU8HHk6v6etn/vbS0yHgKsz/8wW4c5+cDf/xj8PHk\nav16YN4873ZNTcBDDwUTQ4/dnDTGTDPGVBpjKquqqqz+pq7Obt35+TIbIKzWrbObCdPSAqxeHXw8\nuWpo8C7a7TKZcB9M1661Oylonw0QVps3y/5vw3Y8ubJmjd3+FfacbNgAFBbatd28OZgYeqxwk7yH\nZAXJirKyMqu/KSmxW3cmA/Tt243gAjZwIJBIeLcrLJQ7z2GVSAAFBXZtCwqAZDLYeLpjyBC7wZVM\nSv7Cqk8f+4Op7XhyZehQu4NQIiEzTsJqwAD72S/9+gUTg9PpgOecA8Ri2dsYA0yebFcYXTnlFLvL\nBnl5wJlnBh9ProyR6WVegys/X6ZvGtM7ceXi9NNle3vJZCR/YZVIyP7vta1jMRlPYfb97wPxuHe7\nlhaZihpW/foBhxzi3S6RAH7wg2BicFq4L7vM+6womQSuu6534snVoEHAySdnP7jE48DEiTJfNcyu\nvto7J/E4MH1678STqz33lO2drVAkk8BJJ0n+wuy667zf3RQWyngKswMOAMrLs+9fqRRw7rnhfocN\nAD/9qXdOYjHgRz8KKACvaScAHgOwCkAzgBUA/sXrb/zM437pJZk2s/20rfa5kA89ZD+dxqUtW8gJ\nE6QvnU0LGjuW3LzZdZR2nn5atv3207YKCuTnM2e6jtDO5s3kuHFd52TCBMlbFPz+97Lt8/O37Ucs\nJn156SXXEdpZv57cb7/Op9Kl0+TRR5ONja6jtHPnnZ0/g1JYSBYV+Z/GjJ6cx53L4qdwk+Tnn5MX\nX0wWF/MfD6yccw65YIG/jrvW1EQ+/LA82JGXJ8u++8oDKw0NrqPzZ9Eictq0jqKXTpM//KH8PEoa\nGsgHHpA8tOdk1CjyD3+QfEXJhx/KuGh/oKi4WMbNkiWuI/Nnyxby7rvJPffsyMm4ceQTT5AtLa6j\n82fuXPK00zqe4ejTh7ziCn/zt9v5KdyGAczpqqioYGVlZU5/m8nYXZsMu/bNGubrwLY0J+GjOQmf\n7ubEGDOXZIVNW8s5BL1nZ9gZgZ1jR2ynOQkfzUn49GZOdpL0K6XU14cWbqWUihgt3EopFTFauJVS\nKmK0cCulVMRo4VZKqYjRwq2UUhGjhVsppSJGC7dSSkWMFm6llIoYLdxKKRUxWriVUipitHArpVTE\naOFWSqmI0cKtlFIRo4VbKaUiRgu3UkpFjBZupZSKGC3cSikVMaH5zMnGRuDVV4E1a4D+/YGjjgLS\naddR5ebTT4H33pMPQh03Dth/f9cR5aahAXjlFaCqCigtlZwkk66jys1HHwHz5snXBx0E7LOP23hy\nVVcnOamuBgYOlJwUFrqOKjfz5wMLFshnNR5yCLDnnq4jys3mzZKTTZuAIUOASZOAWCzgF7X9OHg/\nS3l5ufVH0re2kjfcQJaUyJJOk8XFZCpFXnop2djo/2PuXVmwgDzoIDKZJIuKZEkmydGjyffecx2d\nvZYW8tprJf6tc1JUJD9vaXEdob05c8gxY3bMyYEHSr6iorGRvOwyGRfFxZKT9jFzww0yjqJi9mxy\n5EjpSzotOUkkyEMPJT/91HV09urryQsukP2pPSfFxWTfvuSvfkVmMv7WB6CSljXWrhEwGcCnAJYA\nuMarvW3hzmTIs8+WBMr56bZLMklOmkQ2N/vbAC7Mmyc7YGf9AKSPb73lOkpvra3kCSd0nZNUSn4f\nhULx1ltd9wOQfL3/vusovTU3yzhIJrvOydln+y8ULrzyStc5MUYORAsXuo7SW0MDefDB2XNy0UX+\n1tmjhRtAPoDPAewBoBDABwBGZvsb28L94otylOpqYAHy+3vv9bcBelsmQ+61V/Z+AOTgweEveI8/\nbpeTxx93HWl2ra2yvb1ysuee4S94995rl5MXX3QdaXZNTWS/ftn7YQw5dqzrSL3dckvXRXvrnPz1\nr/br9FO4bW5OHgRgCckvSDYBeBzACT1xmebmm+WaXTZ1dcCNN8qmCKt33gFWrfJuV1sLvPxy8PF0\nx0032eckzF5+Wba3l9WrgbffDj6eXJH2Obnppt6JKVfPPgu0tGRvQwKLFgEfftg7MeUikwFuuw2o\nr8/ebssWqXFBsCncuwL4cqvvV7T9rNveesuu3fLlcuE/rF5/3TuJAFBTIzcxwqq1VW4Y2fjgA+9B\n6NIrr8j29lJfD7zxRuDh5GzTJmDZMru2YT4AAcALL9jlJJORMRVWX30FbNzo3Y4Mbt+yKdymk5/t\ncP5rjJlmjKk0xlRWVVVZvXgmY9UM+flAU5NdWxcaG+370tAQbCzd0dICmM6y3Qljwl24Gxvt2mUy\n9m1daG6WWRc2bPdBV2zHcCYT7vHuJydBjRGbl18BYLetvh8KYOX2jUjeQ7KCZEVZWZnViw8bZtUM\nsRgwYIBdWxdGjgSKirzbpdPA6NHBx5OreNx+Ow8YACQSwcbTHQccYDedtKhI8hdW/ftLXmzYjidX\nxo61m06aSAD77Rd8PLkaPNj+IBnUFEebwj0HwN7GmN2NMYUATgPwXE+8+JVXAqlU9jaFhcCFF8pZ\nd1iddJLdmWomA5x+evDxdMell3oX5ERC2oXZ6afbDS5jJH9hlZ8PXHCB91ztVErGU5idd57dvap4\nHJg8Ofh4cpVMAmeeCRR4PAVTVARMnx5QEDZ3MAFMAbAYMrvkJ17tbWeV1NSQu+9OFhR0fYd5wABy\n9Wr7O7Ou3HVX9qlnqRR5442uo/RWXU0OGULm53fej/x8+X11tetIvd14o3dO7rzTdZTeVq8mS0vJ\nvLzO+1FQIOOopsZ1pN6uuSb7DJlUinz4YddRelu2TOZrG9N5P2Ixcr/9ZNqgLfT0PG6/i58HcFau\nJEeNkjm1W2+E4mJy2DBy8WL7jrt2663yIMHW04QSCVl++tPwTztrt2wZuffeO85LLyqSny9d6jpC\nO5kMef31HTlo70cyKd//6leuI7S3eDE5fPi2OTFGvh81SsZRFGQy5I9/LNs/Ht+2YCcS5G9/6zpC\nex9+KFNOi4u3zUk6LQ94VVX5W5+fwm0YwDy7iooKVlZW+jjrB2bPBu65B/jyS2CXXYDzzweOOcb+\nJkBYrFsH3Hcf8Oc/S78OO0ze6g4a5Doyf0jgtdeAe++VqY6DBwM/+IE8zmt7AzMsVq+WfWv2bPl+\n4kTghz+Ux/ijJJMBXnoJuP9+YO1aYLfdgGnTZB+LWk5WrADuvht4910Z48ccI2O+Xz/XkfnT0iKz\nZR56SP4bguHD5dLuwQf7z4kxZi7JCqu2YSjcSin1deencEfsfFYppZQWbqWUihgt3EopFTFauJVS\nKmICuTlpjKkCYPk/LHSqFMC6HgrHpZ2lH8DO0xftR/jsLH3pbj+Gk7R67DyQwt1dxphK27urYbaz\n9APYefqi/QifnaUvvdkPvVSilFIRo4VbKaUiJqyF+x7XAfSQnaUfwM7TF+1H+Owsfem1foTyGrdS\nSqmuhfWMWymlVBdCVbiNMZONMZ8aY5YYY65xHU+ujDH3G2PWGmM+ch1LdxhjdjPGvG6MWWiM+dgY\nE/L/hbtrxpiEMeY9Y8wHbX253nVM3WGMyTfGzDPGzHIdS3cYY5YaYz40xsw3xkT2PzgyxvQ1xjxl\njFnUNl4OCfT1wnKpxBiTD/k/v4+CfOrOHACnk/zEaWA5MMYcDqAWwO9J7u86nlwZYwYDGEzyfWNM\nMYC5AE6MaE4MgDTJWmNMDMBfAVxK8l3HoeXEGHMFgAoAJSSPcx1ProwxSwFUkIz0PG5jzEMAZpO8\nr+0DZ1IkLT6ZMjdhOuMO7NPkexvJNwFUu46ju0iuIvl+29c1ABaihz4oure1/ZfH7Z/7HmtbwnHW\n4pMxZiiAYwHc5zoWBRhjSgAcDuB3AECyKciiDYSrcAf2afKq+4wxIwCMBfA3t5Hkru3ywnwAawG8\nQjKqffk1gOkAQv7xwFYI4GVjzFxjzDTXweRoDwBVAB5ou3x1nzHG4hNPcxemwm31afKq9xljigDM\nBHAZyc2u48kVyVaSYyAfeH2QMSZyl7GMMccBWEtyrutYesihJMcB+A6Ai9ouM0ZNAYBxAO4iORZA\nHYBA79GFqXBbfZq86l1t14NnAniE5NOu4+kJbW9j3wAQ4o+k7dKhAL7bdm34cQATjTEPuw0pdyRX\ntv27FsAzkEumUbMCwIqt3sE9BSnkgQlT4Q7s0+RVbtpu6P0OwEKSt7mOpzuMMWXGmL5tXycBfBvA\nIrdR+UfyWpJDSY6AjJHXSJ7lOKycGGPSbTe90XZp4WgAkZuJRXI1gC+NMfu0/WgSgEBv4Ht8wHzv\nIdlijLkYwEsA8gHcT/Jjx2HlxBjzGIAjAZQaY1YA+E+Sv3MbVU4OBXA2gA/brg0DwL+TfNFhTLka\nDOChttlLeQCeJBnpqXQ7gYEAnpHzAxQAeJTkn9yGlLNLADzSdtL5BYDzgnyx0EwHVEopZSdMl0qU\nUkpZ0MKtlFIRo4VbKaUiRgu3UkpFjBZupZSKGC3cSikVMVq4lVIqYrRwK6VUxPx/TMRQQSh6BNsA\nAAAASUVORK5CYII=\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7fce55956750>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# Use control+Enter to run this cell many times and observe how \n",
    "# the robot's belief that it is in each cell (represented by the\n",
    "# size of the corresponding circle) changes as the robot moves.\n",
    "# The true position of the robot is given by the red star.\n",
    "\n",
    "# Run this cell about 15-25 times and observe the results\n",
    "simulation.run(1)\n",
    "simulation.show_beliefs()\n",
    "\n",
    "# If everything is working correctly you should see the beliefs\n",
    "# converge to a single large circle at the same position as the \n",
    "# red star.\n",
    "#\n",
    "# When you are satisfied that everything is working, continue\n",
    "# to the next section"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Part 3: Identify and Reproduce a Bug\n",
    "Software has bugs. That's okay.\n",
    "\n",
    "A user of your robot called tech support with a complaint\n",
    "\n",
    "> \"So I was using your robot in a square room and everything was fine. Then I tried loading in a map for a rectangular room and it drove around for a couple seconds and then suddenly stopped working. Fix it!\"\n",
    "\n",
    "Now we have to debug. We are going to use a systematic approach.\n",
    "\n",
    "1. Reproduce the bug\n",
    "2. Read (and understand) the error message (when one exists)\n",
    "3. Write a test that triggers the bug.\n",
    "4. Generate a hypothesis for the cause of the bug.\n",
    "5. Try a solution. If it fixes the bug, great! If not, go back to step 4."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Step 1: Reproduce the bug\n",
    "The user said that **rectangular environments** seem to be causing the bug. \n",
    "\n",
    "The code below is the same as the code you were working with when you were doing integration testing of your new feature. See if you can modify it to reproduce the bug."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "NOTE! The robot doesn't have a working sense function at this point.\n"
     ]
    },
    {
     "ename": "IndexError",
     "evalue": "list index out of range",
     "output_type": "error",
     "traceback": [
      "\u001b[0;31m\u001b[0m",
      "\u001b[0;31mIndexError\u001b[0mTraceback (most recent call last)",
      "\u001b[0;32m<ipython-input-13-792590274d38>\u001b[0m in \u001b[0;36m<module>\u001b[0;34m()\u001b[0m\n\u001b[1;32m     22\u001b[0m \u001b[0;31m# remember, the user said that the robot would sometimes drive around for a bit...\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m     23\u001b[0m \u001b[0;31m# It may take several calls to \"simulation.run\" to actually trigger the bug.\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m---> 24\u001b[0;31m \u001b[0msimulation\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mrun\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;36m1\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m     25\u001b[0m \u001b[0msimulation\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mshow_beliefs\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m/home/workspace/simulate.pyc\u001b[0m in \u001b[0;36mrun\u001b[0;34m(self, num_steps)\u001b[0m\n\u001b[1;32m    103\u001b[0m                         \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0msense\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    104\u001b[0m                         \u001b[0mdy\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mdx\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mrandom_move\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 105\u001b[0;31m                         \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mmove\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mdy\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0mdx\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m",
      "\u001b[0;32m/home/workspace/simulate.pyc\u001b[0m in \u001b[0;36mmove\u001b[0;34m(self, dy, dx)\u001b[0m\n\u001b[1;32m     55\u001b[0m                 \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mtrue_pose\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0;34m(\u001b[0m\u001b[0mnew_y\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mnew_x\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m     56\u001b[0m                 \u001b[0mbeliefs\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mdeepcopy\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mbeliefs\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m---> 57\u001b[0;31m                 \u001b[0mnew_beliefs\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mlocalizer\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mmove\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mdy\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mdx\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mbeliefs\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mblur\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m     58\u001b[0m                 \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mbeliefs\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mnew_beliefs\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m     59\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m/home/workspace/localizer.pyc\u001b[0m in \u001b[0;36mmove\u001b[0;34m(dy, dx, beliefs, blurring)\u001b[0m\n\u001b[1;32m     33\u001b[0m             \u001b[0mnew_j\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0;34m(\u001b[0m\u001b[0mj\u001b[0m \u001b[0;34m+\u001b[0m \u001b[0mdx\u001b[0m \u001b[0;34m)\u001b[0m \u001b[0;34m%\u001b[0m \u001b[0mheight\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m     34\u001b[0m             \u001b[0;31m# pdb.set_trace()\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m---> 35\u001b[0;31m             \u001b[0mnew_G\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0mint\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mnew_i\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0mint\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mnew_j\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m]\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mcell\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m     36\u001b[0m     \u001b[0;32mreturn\u001b[0m \u001b[0mblur\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mnew_G\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mblurring\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;31mIndexError\u001b[0m: list index out of range"
     ]
    }
   ],
   "source": [
    "from simulate import Simulation\n",
    "import simulate as sim\n",
    "import helpers\n",
    "reload(localizer)\n",
    "reload(sim)\n",
    "reload(helpers)\n",
    "\n",
    "R = 'r'\n",
    "G = 'g'\n",
    "\n",
    "grid = [ \n",
    "    [R,G,G,G,R,R,R],\n",
    "    [G,G,R,G,R,G,R],\n",
    "    [G,R,G,G,G,G,R],\n",
    "    [R,R,G,R,G,G,G],\n",
    "]\n",
    "\n",
    "blur = 0.001\n",
    "p_hit = 100.0\n",
    "simulation = sim.Simulation(grid, blur, p_hit)\n",
    "\n",
    "# remember, the user said that the robot would sometimes drive around for a bit...\n",
    "# It may take several calls to \"simulation.run\" to actually trigger the bug.\n",
    "simulation.run(1)\n",
    "simulation.show_beliefs()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "simulation.run(1)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Step 2: Read and Understand the error message\n",
    "If you triggered the bug, you should see an error message directly above this cell. The end of that message should say:\n",
    "\n",
    "```\n",
    "IndexError: list index out of range\n",
    "```\n",
    "\n",
    "And just above that you should see something like\n",
    "\n",
    "```\n",
    "path/to/your/directory/localizer.pyc in move(dy, dx, beliefs, blurring)\n",
    "     38             new_i = (i + dy ) % width\n",
    "     39             new_j = (j + dx ) % height\n",
    "---> 40             new_G[int(new_i)][int(new_j)] = cell\n",
    "     41     return blur(new_G, blurring)\n",
    "```\n",
    "\n",
    "This tells us that line 40 (in the move function) is causing an `IndexError` because \"list index out of range\".\n",
    "\n",
    "If you aren't sure what this means, use Google! \n",
    "\n",
    "Copy and paste `IndexError: list index out of range` into Google! When I do that, I see something like this:\n",
    "\n",
    "![Search Results](http://i.imgur.com/gleBmBy.png)\n",
    "\n",
    "Browse through the top links (often these will come from stack overflow) and read what people have said about this error until you are satisfied you understand how it's caused."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Step 3: Write a test that reproduces the bug\n",
    "This will help you know when you've fixed it and help you make sure you never reintroduce it in the future. You might have to try many potential solutions, so it will be nice to have a single function to call to confirm whether or not the bug is fixed"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "# According to the user, sometimes the robot actually does run \"for a while\" \n",
    "# - How can you change the code so the robot runs \"for a while\"?\n",
    "# - How many times do you need to call simulation.run() to consistently\n",
    "#   reproduce the bug?\n",
    "# Modify the code below so that when the function is called \n",
    "# it consistently reproduces the bug.\n",
    "def test_robot_works_in_rectangle_world():\n",
    "    from simulate import Simulation\n",
    "    import simulate as sim\n",
    "    import helpers\n",
    "    reload(localizer)\n",
    "    reload(sim)\n",
    "    reload(helpers)\n",
    "\n",
    "    R = 'r'\n",
    "    G = 'g'\n",
    "\n",
    "    grid = [ \n",
    "        [R,G,G,G,R,R,R],\n",
    "        [G,G,R,G,R,G,R],\n",
    "        [G,R,G,G,G,G,R],\n",
    "        [R,R,G,R,G,G,G],\n",
    "    ]\n",
    "\n",
    "    blur = 0.001\n",
    "    p_hit = 100.0\n",
    "    for i in range(1000):\n",
    "        simulation = sim.Simulation(grid, blur, p_hit)\n",
    "        simulation.run(1)\n",
    "    \n",
    "test_robot_works_in_rectangle_world()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Step 4: Generate a Hypothesis\n",
    "In order to have a guess about what's causing the problem, it will be helpful to use some Python debuggin tools\n",
    "\n",
    "The `pdb` module (`p`ython `d`e`b`ugger) will be helpful here!\n",
    "\n",
    "#### Setting up the debugger \n",
    "1. Open `localizer.py` and uncomment the line to the top that says `import pdb`\n",
    "2. Just before the line of code that is causing the bug `new_G[int(new_i)][int(new_j)] = cell`, add a new line of code that says `pdb.set_trace()`\n",
    "3. Run your test by calling your test function (run the cell below this one)\n",
    "4. You should see a text entry box pop up! For now, type `c` into the box and hit enter to **c**ontinue program execution. Keep typing `c` and enter until the bug is triggered again"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "test_robot_works_in_rectangle_world()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### Using the debugger\n",
    "The debugger works by pausing program execution wherever you write `pdb.set_trace()` in your code. You also have access to any variables which are accessible from that point in your code. \n",
    "\n",
    "Try running your test again. This time, when the text entry box shows up, type `new_i` and hit enter. You will see the value of the `new_i` variable show up in the debugger window. Play around with the debugger: find the values of `new_j`, `height`, and `width`. Do they seem reasonable / correct?\n",
    "\n",
    "When you are done playing around, type `c` to continue program execution. Was the bug triggered? Keep playing until you have a guess about what is causing the bug."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Step 5: Write a Fix\n",
    "You have a hypothesis about what's wrong. Now try to fix it. When you're done you should call your test function again. You may want to remove (or comment out) the line you added to `localizer.py` that says `pdb.set_trace()` so your test can run without you having to type `c` into the debugger box."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [
    {
     "ename": "NameError",
     "evalue": "name 'test_robot_works_in_rectangle_world' is not defined",
     "output_type": "error",
     "traceback": [
      "\u001b[0;31m\u001b[0m",
      "\u001b[0;31mNameError\u001b[0mTraceback (most recent call last)",
      "\u001b[0;32m<ipython-input-14-9c7d38501e24>\u001b[0m in \u001b[0;36m<module>\u001b[0;34m()\u001b[0m\n\u001b[0;32m----> 1\u001b[0;31m \u001b[0mtest_robot_works_in_rectangle_world\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m",
      "\u001b[0;31mNameError\u001b[0m: name 'test_robot_works_in_rectangle_world' is not defined"
     ]
    }
   ],
   "source": [
    "test_robot_works_in_rectangle_world()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Congratulations!\n",
    "You've implemented your first feature and successfully debugged a problem the robot was having with rectangular environments. Well done."
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 2",
   "language": "python",
   "name": "python2"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 2
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython2",
   "version": "2.7.13"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}