writeup.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Two Dimensional Histogram Filter - 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": {},
   "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": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAD8CAYAAACMwORRAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4xLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvAOZPmwAAIABJREFUeJzt3Xl8FPX9P/DXO+dmswlnlECAlMOTokJE1HprvdEqWjyw1CqKon6lgketrUdttVU8f0U86oG2eKCCxaL1rFqQcGtBroIgV5BACLmzr98fn6WEsMnObGYzs7Pv5+MxD3ezn5l9+2bmvZ/9zHxmhSSUUkr5S5rbASillHKeFnellPIhLe5KKeVDWtyVUsqHtLgrpZQPaXFXSikf0uKulFI+pMVdKaV8SIu7Ukr5UIZbb9y1a1cWFxe79fZKKZWU5s2bt5VkQax2rhX34uJilJaWuvX2SimVlERkrZV2OiyjlFI+pMVdKaV8yLVhmXiQwOefA5MmAatWmb/16QNccw1w3HGAiLvxec2aNcD/+3/A7NlAdTXQtStw6aXARRcBgYDb0XlLZSUwZQrwxhvA9u1AKASccAJw7bVAt25uR+ct4TDw3nvAs88C69YB6enAQQcB118PDBrkdnTe85//AE8+CSxaBNTVAYWFwKhRwLnnAhmJrMAkXVkGDx5MO/7xD7K4mMzNJUVIU+rN41CI7N2bnDnT1iZ9a+1a8qSTyECAzMrakyvA5CoUIn/zG7Kx0e1I3VdbS44dS+bkmH2raa4CATI7mxw2jCwrcztSb5gyhdx/f7MPNc1VWhoZDJKHHEL++99uR+kNS5aQgwebfSsjY+985eWRnTuTTzxBhsP2tguglBZqbFIU9xdeMAlqmpxoS04O+dxz9hLlN8uWmZ0mPb31XAWD5Pnnkw0Nbkfsnpoa8phjYu9bmZlkjx7kd9+5HbG77rvP7DexjsNgUDtaX3yx7wdgS7kaO9ZegfdNcf/Xv6wV9qYF/uOPrSfKTyoqyG7d9v5mE2vHGjfO7ajdc/HF1vetjAzygAPIujq3o3bHq6/aOw6DQfLrr92O2h3r1pH5+fZy9dhj1rdvtbhbPqEqIukiskBE3onyWraITBWRlSIyR0SKnRo2uu02M15sVXW1WScVvfACsHOn2WWsqKoyY/Ll5YmNy4tWrwamT7e+bzU0ABs2ADNmJDYuLyKB8ePtHYc1NcC99yYuJi+bONH8/1tVVQXcdRdQX+9sHHaulrkJwNIWXvsFgHKS/QBMBPBAWwMDzEnTefPsr7doEbB8uRMRJA8S+OMfgV277K2XlgY8/3xCQvK0xx8HGhvtrVNZCTz4YGLi8bLPPgO+/97eOuEw8NZbqddxqKkBnn7anDi1o7HRdDacZKm4i0gRgLMBPNNCk/MAvBB5/DqAU0Tafu3KG2/YPwAB8wn42mttfffksnSp/QMQML2GZ591Ph6ve/nl+HpKCxYA27Y5H4+XvfSS/U4DYK4EeWef7/n+9skn8V21t3Mn8NxzzsZitef+CIAJAMItvN4DwDoAINkAYAeALm0NbuPG+A7A3V+hU0lZWfyXVW3d6mwsyaCiIr71srJMrlPJd99ZH+prqq4u9XJVVma+tcRj0yZnY4lZ3EXkHABbSLY2QBLts2qf3UFERotIqYiUlln4V8/OjtmkRal2HXdmpjvrJqv09PjWI02BTyXx/v+mpaVerjIz459v43SurPTcjwUwTETWAPgbgJNFZEqzNusB9AQAEckA0AHAPl9eSU4mWUKypKAg5n1vcNBBZjKJXcGgWTeV9OkD1NbGt+6BBzobSzLo1Su+9cLh1JvUNHBgfB2tzEygXz/n4/Gy/v3j67mnpwOHHupsLDGLO8nbSRaRLAYwAsCHJC9v1mw6gJ9FHg+PtInji9zeLroovkSRwIgRbX335NKtG3DssfbXy8sDxo1zPh6vGzfOfschI8PsVzk5iYnJq66+Or71srOB005zNhavO+IIoHt3++tlZwM33OBsLHHfW0ZE7hGRYZGnzwLoIiIrAYwD4MjFiLm5wMiR9oYNMjKASy4xRSvVTJhgcmZHTg5w+umJicfLLr3UfschMxO4+ebExONlRUXA8cebYRarAgHgppviH/5KViLArbfaPw779QMOO8zhYKxcDJ+Ixeokpi1byMJCM7051mSAtDQzNXrTJusTAvwkHCYvuMD6ZJOcHHLWLLejds+LL1qbcbl7osn117sdsXtWriQ7dLA+4evAA8mdO92O2h21teRRR5lbV1jJV24uOX++9e3D6UlMbikoMDcL69699ZOkgYAZmvj8c2D//dsvPi8RAV55BTj11NZ7Dunp5rzEiy8CP/5x+8XnNSNHmrkBwWDrJ8Fyc01P/7HH2i82r+nbF/joI6BTp9ZP/OXkmF7oJ5/Ed77MD7KyzI3VDj/c7Fstycw0IwwzZ5rhHMdZ+QRIxGL3xmHl5eTvfkd27WpuurP7Blh5eeZv995Lbttma5O+1dhI/vWv5GGHmd55Xp7pHeTnm5thXXFF6k4Nj+bLL803nkDA5Cg31+QsEDD3npkxw/7Nnfxq40ZywgSTp9252n0sFhWRjz5K7trldpTeUFtLPvUU2bfvnn0qFNqTt7Fjyf/+1/52YbHnLmz7ec+4lJSUMJ5fYmpsBD79FFi/3jzv0cPcmjXVxvas+uors1RXm17XSScBHTq4HZU3lZWZfauiwvTWBw1Kvas9rKqrMz35TZvMea7iYuCYY/S229GQQGkpsGKFyVtBgTkOW+vVt0ZE5pEsidku2Yq7UkqlMqvF3fNj7koppezT4q6UUj6kxV0ppXxIi7tSSvmQFnellPIhLe5KKeVDWtyVUsqHtLgrpZQPaXFXSikf0uKulFI+pMVdKaV8yMpvqAZE5EsRWSQiX4vI3VHajBKRMhFZGFmuSky4SimlrMiw0KYWwMkkK0UkE8BnIvIuydnN2k0lOdb5EJVSStkVs7hH7h9cGXmaGVncuZWkUkopSyyNuYtIuogsBLAFwPsk50RpdqGILBaR10Wkp6NRKqWUssVScSfZSPJwAEUAhojIgGZNZgAoJjkQwD8BvBBtOyIyWkRKRaS0rKysLXErpZRqha2rZUhuB/AxgDOa/f17krWRp08DGNzC+pNJlpAsKSgoiCNcpZRSVli5WqZARDpGHucAOBXAsmZtCps8HQZgqZNBKqWUssfK1TKFAF4QkXSYD4NXSb4jIvfA/FDrdAA3isgwAA0AtgEYlaiAlVJKxaa/oaqUUklEf0NVKaVSmBZ3pZTyIS3uSinlQ1rclVLKh7S4K6WUD2lxV0opH9LirpRSPqTFXSmlfEiLu1JK+ZAWd6WU8iEt7kop5UNWbhzmOfX1wLZt5nHnzkBmprvxeFk4DGzfDtTUAB07AsGg2xF5W2UlsGMHkJsLdOgAiLgdkXfV1gLl5UB6OtCpE5CRlNWkfTQ2mlzV15tcBQKJf8+k6rnPnw9cdhmQlwcUF5slFAIuuQSYN8/t6Lxl0ybgt78FCgqAwkKgXz9TrIYMAd54A2hocDtC76ipAV56CRgwwBx4/fsD++8PdO8O/OlPezoSCiCBTz8Fzj3XHId9+gC9egH5+cA11wDLlsXeRir573+Bm282HauiIqBvX5O3U04BZs0yna+EIenKMnjwYFpVU0MOH04Gg2R6Oml2sT1LWpp57Sc/MW1T3eTJZCBglua5Asi8PLJ3b3LVKrcjdd/ChWRBARkKRc9VMGiWN990O1L37dhBHnecyZXIvrnKyCBzcsjrriMbGtyO1l3hMHnXXeYYzMqKvm+FQuTAgeSWLfa2DXOr9Zg11vPFvaGBPO00s9NES1DTJSeHPPlksr7eXrL85MknTTGKlau0NLJLF3LtWrcjds+SJS0X9Wj71rRpbkfsnl27yB/+kMzOjp2rYJAcOdIUuFQ1fry14zAz03S0tm2zvm2rxd3zwzITJwKffw5UV8duW10NzJ4NPPRQ4uPyomXLgFtuAaqqYrfdPRZ/wQWJj8uLwmHgrLOAXbusta+uBi6/HNi8ObFxedX48cCKFWacPZaqKmDaNGDq1MTH5UUffgg8+aS147C+Hti4EbjqKufjsPIzewER+VJEFonI1yJyd5Q22SIyVURWisgcESl2Irhw2Ix5WknSblVVZp3GRiciSC4PP2x2FqsaG4GlS4HFixMXk1e9/745wUUbv1UTDgNPP524mLxq1y7g+efNuQk769x/f8JC8rT777dXs+rqgJkzne84WOm51wI4meRhAA4HcIaIDG3W5hcAykn2AzARwANOBPfee9Z7Vk3V1gLvvutEBMmjshKYMsX+idLaWvPtKNU8+KDJmR01NcCjj6Zex+Hll+O7amjlSmDhQufj8bJ164DPPotv3aeecjaWmMU9Msyz+zDIjCzN+zvnAXgh8vh1AKeItP0isi++sH8AAsDOnWbdVLJ0aXyXhDY2Ah9/7Hg4njd3bnzrVVWZr9Gp5IMP4utkAcCcOc7G4nVz5wLZ2fbXq6kB/vlPZ2OxNOYuIukishDAFgDvk2z+T9YDwDoAINkAYAeALlG2M1pESkWktKysLOb7bt9uJbroUu3ytcrK+K/JtvMV0i/sDDE0lZ5uOg+ppKIivvUaGlIvV5WV8V/e6HSuLBV3ko0kDwdQBGCIiAxo1iRaWdlnNJPkZJIlJEsKCgpivm/nzvEXrC77fLT4W36+vfHjpkIhZ2NJBjk58a1XX29ynUo6dYpvvczM1MtVfj6QFudlKh06OBuLrTBIbgfwMYAzmr20HkBPABCRDAAdALS573zyyWamoF2hkJkkkEoOPTS+4p6ZCZx5pvPxeN2JJ8bXcejc2UxuSiVnnx1fB4AETjjB+Xi87JhjzAlSu4JBMzHMSVaulikQkY6RxzkATgXQfB7adAA/izweDuDDyPWYbXLccfH1wDt2BE46qa3vnlwCAXM5VVaWvfXS04GbbkpMTF52yy32b8UQDJr1Uu2WBMOHx/f/PHAgcOCBzsfjZfvtB5x+uv18hcPAlVc6G4uVnnshgI9EZDGAuTBj7u+IyD0iMizS5lkAXURkJYBxAG5zIjgR4Fe/sncQBoPAHXek3gEIADfeaO/+HllZwNFHm+n2qeZHPwJ69jQfblalpwM//3niYvKq7Gxg7Fh7x2FuLvDrXycuJi+7/XZ7w36BADBiRPzDXy2yMtMpEYvVGarhMHnppdZmewWD5E9/mtoz41591dps3sxMsmdPcutWtyN2z5o1ZOfOZraulX3ro4/cjtg99fXkCSdY27eCQTNDM5U99JC1mhUImFsQVFVZ3zb8MkNVxNzU6eqrzSdctMuMsrLMa1deGf81uX5x0UVmZmBubvRxUhHz2g9/CJSWpt6J56Z69zY3nOvTp+Ux5VDIDPPNmmXG6VNVRobJwbnnml5ptMtuAwGz3HEH8IAjM12S17hxZv5IIBD9G096uvn7cceZGfjxnuBvlZVPgEQsdm4cttu335K33WbuiZKWZpbOnckJE1L7HinRVFaSzzxDHnDAnput5eSQF1xAfvFFan+7aa6xkXz/fXMPo933TsnIIA8/nHzlFb0ZXXPffGNuDpafb24glpZGdutG3ncfuXmz29F5S3k5+fDD5ltyWprJV24uOWoUuWhRfNuExZ67sO3nPeNSUlLC0tLSuNfffS1pvJcdpRLS5MvO+HIqa2zUXFkVDptvg6n8bdkqp45DEZlHsiRWu6S9vb4WdetEtFjZobmyTo9D69r7ONR/GqWU8iEt7kop5UNa3JVSyoe0uCullA9pcVdKKR/S4q6UUj6kxV0ppXxIi7tSSvmQFnellPIhLe5KKeVDWtyVUsqHrPwSU08R+UhElorI1yKyz+/2iMiJIrJDRBZGlrsSE65SSikrrNw4rAHAL0nOF5E8APNE5H2S/2nW7l8kz3E+RKWUUnbF7LmT3EhyfuTxTgBLAfRIdGBKKaXiZ2vMXUSKARwBYE6Ul48WkUUi8q6IHOpAbEoppeJk+X7uIhIC8AaA/yNZ0ezl+QB6k6wUkbMAvAVgn59dFpHRAEYDQK9eveIOWimlVOss9dxFJBOmsL9Mclrz10lWkKyMPJ4JIFNEukZpN5lkCcmSgoKCNoaulFKqJVaulhEAzwJYSvLhFtp0i7SDiAyJbPd7JwNVSillnZVhmWMBjASwREQWRv52B4BeAEByEoDhAMaISAOAagAj6NaPsyqllIpd3El+BqDVn78l+QSAJ5wKSimlVNvoDFWllPIhLe5KKeVDWtyVUsqHtLgrpZQPaXFXSikf0uKulFI+pMVdKaV8SIu7Ukr5kBZ3pZTyIS3uSinlQ1rclVLKhyzfz90rdu0C3n4bWL/ePO/RAzjvPCAUcjcuLwqHgQ8/BJYsAaqrgU6dgLPOAnr3djsyb1q2DPjgA2DHDiA3J4ySIWk45hhAWr2zUmoqLzfH4ebNQHo6UFwMnHsukJ3tdmTe09AAvPsusHw5UFcHFBSYXO2/f4LfmKQry+DBg2nHunXkNdeQwSCZl0dmZJglFCJzcsirrya//dbWJn2rtpb805/Ibt1MfrKyyLQ0k6dAgDzpJPLTT92O0jv+/nfyyCP35KerbOV36M6i4PcsLiYnTSIbG92O0huWLSMvucTkKTeXTE83x2Fenll++Uty61a3o/SGigry178mO3UyucnMNPkKBk3+hg0jFyywv10ApbRQY5OiuC9caBKUkWEijrZkZJAdO5Lz59tPlp9UVJBHHWUKVUu5AswO9sQTbkfrvl/9yuSiaW5G4TkS4M/wl//l6qyzzIdmKvvwQ1PQ09Ja3q+yssju3clVq9yO1l2bNpH9+5PZ2S3nSsTsW2+8YW/bvinua9eawt5aoWq6dOxI/ve/9pLlFw0N5PHHt75DNV1ycsi//c3tqN3z0EP7FnaA/BjHkwA/wgl75eqii8hw2O2o3bFwoSnsVvartDSyR4/U7cHv2kUedFDrndHmx+HHH1vfvtXibuWXmHqKyEcislREvhaRm6K0ERF5TERWishiERnk1LDR7bcDFc1/sbUVFRXArbc69e7JZcYMYP58oLbWWvvqauDaa4H6+sTG5UXbtwO/+hVQVbX33/NQgaMwGwAwFLORB7PzVVcDf/878OWX7R2pN4wZY853WREOA2VlwB//mNiYvOqZZ4BvvzVj7VZUVwNXXWVKvZOsXC3TAOCXJA8GMBTA9SJySLM2Z8L8IHZ/mB/A/rMTwZWXA9OmAY2N1tcJh4Hp04HvU/BH/h54AKistLdOYyPw1luJicfLnn8e6CzlGIhFey1X4RnUwZwVrEMWrsIz/3vtgOpFeO3ORWbHTCErVgALFthbp64OmDTJ/DeVkOZDrXmnIZaNG4E5cxwPxt5wCoC3AZzW7G9PAbikyfNvABS2th0rwzKTJkX/2mzla06qjSd/+605SWM3VwB57LFuR9/++vQh/4AJJMAqBLgd+f9bmian6d+rYBJce/MEt8NvVxMmmJOBdvervDzy7bfdjr59ffGFuYjBbq7S0sjLL7f2HnBqWKYpESkGcASA5p8xPQCsa/J8feRvbbJypf1PQMB8zVmxoq3vnlzWro3/MrQ1axwNJSls2ADcjt/jdtwPQtABFf9bmmr6d0Lwm8z7sW7M712K2h3LlsU3dFdXl3r71po18V06Gw6bSyWdZLm4i0gIwBsA/o9k81HwaP87+4wgichoESkVkdKysrKY72lnOMbJdZNROOzOusmKBIg0/AG340R8jI3ohmoEoratRjY2oBAn4BM8nH07wik29y/eY4lMvX0rHI5/7NzpmmVpLxWRTJjC/jLJaVGarAfQs8nzIgAbmjciOZlkCcmSgoKCmO/bs2d8vdGsLKBXL/vrJbPCwvhPjCZ8MoUHde685/FcDMGB+AYbURi17UZ0x0FYhlIc+b9JKKmkuBhIi+PzLDsb6NbN8XA8rbAwvlwBpt45ycrVMgLgWQBLST7cQrPpAK6IXDUzFMAOkhvbGtzFF8f3FSctDfjpT9v67smlf//4do7cXHPFTKq58sq9Ow51yEIhou+y3bAJtZGTrMceC3Ts2B4ReseoUUAg+peaVjU0AOec43g4nnb88WbGrl2hEDB6tLOxWPmMORbASAAni8jCyHKWiFwrIrvLwkwAqwGsBPA0gOucCK5HD+CEE+yvd+yxqddzB8wloLm59tYhgcsuS0w8XjZmzN7PT8P7qEMWAGAXgqhFFnYhCACoRyZOw/sIhYAJE9o7UveVlNg/njIyzH6VarcFycgAxo61P+IQDAKnn+5wMFbOuiZisTqJacECe1fMBIPkvHnWzjr7TVUV2bev9ckTwSD5u9+5HbV7xozZs29NxUUMA6xEkI/gBuZiJx/FWFYih2GAr8rFPPLI1L0NwaxZsWc9N13y88nVq92O2h1btpAFBWYGqtXj8OWXrW8ffpmhSpIzZ1or8Dk55IwZ1pPkR+vXk0VFsS9dCwbJa69N3RmXpJnRe/75ZIecWlYiyO3I54/xj73ydDre5Q7kc5fkctvmOrdDdtUzz8Qu8CLmEsgvvnA7Wnd99RXZubO5l0ys4/D3v7e3bV8Vd5KcO5c8+mhzLXfTwpWZaXa4o44iZ8+2lyS/2rqVHDnS5Kr5h2IoZO79MXlyahf23RobyYfuq+Iz2WPYJ3fTXrnKzTX5u3XUJtZdPcZ8NUpx779P/vCHJi9N7zGTnW32t1NPJZcudTtKb1i71twcLBDYew6KiNm3+vYlp02zv12rxV1M2/ZXUlLC0tJS2+utWAH85S/AqlUmVX36AD//OXDggQkIMslt3w68+CIwd66ZOr7ffsDw4cDJJ8d/Rt+vGhrM7QVmzDCzm/PzzcmxESPsn8dIBQsXmn3r22+BzEzgoIPMSWqnr/jwg82bTc1asgSoqQG6dzfnI446Kr4LRkRkHsmSmO2SrbgrpVQqs1rctf+mlFI+pMVdKaV8SIu7Ukr5kBZ3pZTyIS3uSinlQ1rclVLKh7S4K6WUD2lxV0opH9LirpRSPqTFXSmlfEiLu1JK+ZAWd6WU8iErP7P3nIhsEZGvWnj9RBHZ0eRXmu5yPkyllFJ2ZFho8zyAJwC82Eqbf5FMsV9LVEop74rZcyf5KYBt7RCLUkophzg15n60iCwSkXdF5NCWGonIaBEpFZHSsrIyh95aKaVUc04U9/kAepM8DMDjAN5qqSHJySRLSJYUFBQ48NZKKaWiaXNxJ1lBsjLyeCaATBHp2ubIlFJKxa3NxV1EuomYXwIUkSGRbX7f1u0qpZSKX8yrZUTkrwBOBNBVRNYD+A2ATAAgOQnAcABjRKQBQDWAEXTrh1mVUkoBsFDcSV4S4/UnYC6VVEop5RE6Q1UppXxIi7tSSvmQFnellPIhLe5KKeVDWtyVUsqHtLgrpZQPaXFXSikf0uKulFI+pMVdKaV8SIu7Ukr5kBZ3pZTyoaQq7rW1wCuvAEcfDRQVmWXoUGDKFKCmxu3ovIUEPv8cuOAC4Ac/ALp3Bw4+GLjnHmDzZrej8561a4Hx44EDDgAKC4E+fYCRI4EFC9yOzHsqK4GnngIGDQJ69AB69QJOPBF46y2gocHt6LwlHAZmzQLOOAPo3dvka+BAYOJEYPv2BL85SVeWwYMH045nniHz8sxiSteeJRQyf580ydYmfWvBArJvXzI3lxTZO1eBAJmdTY4cSVZXux2p+3bsIM8+2+QlK2vvXKWnk8EgOXAguWqV25G6Lxwm773X5CQ3d9/jMC+P7NyZfPNNtyP1hg8/JAsLo9esYNDsc+PGkQ0N9rYLoJQWamxSFPf77jPJaJ6gaAn77W/tJcpvZs+OfuA1X3JyyCFDUrvAl5eT/fubD7vWcpWWRnbsSC5d6nbE7gmHyauvtnYc5uSQzz7rdsTumj7des067zyysdH6tn1T3GfMsJakpsl66y3rifKTsjKyQwfrucrJIS+7zO2o3XPSSbEL++5FhOzWjayqcjtqdzz5pP3jcPZst6N2xzffWOtgNc3VnXda377V4h5zzF1EnhORLSLyVQuvi4g8JiIrRWSxiAxyctjozjuBqirr7auqzDqpaPJkoK7OevvqauCNN4BNmxIXk1ctWQLMnm3O41hBmrHm115LbFxe1NgI3H23veOwutqc30lFDz5o7xxgVRXwyCMmZ06yckL1eQBntPL6mQD6R5bRAP7c9rCMxYuBFSvsr7d6deqdCGtsjH8Heeop5+PxuokT7X0QAqa4P/BAYuLxslmz7O9XJPDBB8CGDYmJyat27jQXfTQ22l936lRnY4lZ3El+CmBbK03OA/Bi5BvDbAAdRaTQieDeecf+AQiYT83p052IIHl89VV8hb2mxuyMqebtt+M7AFetArZscT4eL3vtNVO07MrIMB8MqeSTT4DMTPvrVVY6fxw6cSlkDwDrmjxfH/lbm23eHN+lVeFw6g01lJcD6enxrbtjh7OxJIPKyvjWy8oCtrXW1fGheI+lurrUy1V5uak/8di61dlYnCjuEuVvUX8gW0RGi0ipiJSWlZXF3HBubvxB5eXFv24yCgTiXzc727k4kkU8vSvAHLg5Oc7G4nXBYHzrpaenXq4CASAtzqrqdK6cKO7rAfRs8rwIQNSRNpKTSZaQLCkoKIi54cMOA0Ih+wGFQmbdVNK/f3wTuUSAww93Ph6vO/DA+NYTMZOcUslRR8VXeDIygAEDnI/HywYMAOrr7a+XmQkceaSzsThR3KcDuCJy1cxQADtIbnRguzj//Pg/BS+80IkIkkeXLsBZZ5niY0cwCPzyl4mJycvGj7f/7S4zE/jFL8zQTCq58kpzgtSuTp2A445zPh4vO/hg4KCD7K+Xng6MHetsLFYuhfwrgH8DOFBE1ovIL0TkWhG5NtJkJoDVAFYCeBrAdU4Fl50NXHONvSGH7Gzg6qvbNkyRrMaPt9/DKihIvQMQMLdlsNtxSE8HbrghMfF4WdeuwNln2xvKCgbN/mi3s+EHt91mb0g5LQ0oKQH69XM4ECsXwydisTqJaedO8qCDyMzM2JMBMjPJAw4gKyqsTwjwmxtusD7ZJDeXnDfP7Yjd8957ZiKX1YkmDzzgdsTu2bTJTOJKS4udq0CA/NFO6khbAAAOc0lEQVSPyLo6t6N2R2OjmXVqZd8SMbOfV6+2vn04NYnJbaEQ8NlnwCGHtP5pmJtrxlE/+yz1TqY29cgjZuggGGy51xQIAPn55jK1QY5OOUsup50GvPqq2XdaGmrZfVLwzjtNTzRV7b+/mfRVVNT6cRgKmRv7vftu/Cetk11amtmvzjyz9VwFg+Zb0eefm5v7Oc7KJ0AiFrs3DqutJV9+2dzEKSfHTLPv0ME8HjCAnDKFrKmxtUlf+/xz8vzzzfT63bnKyyO7dCHvuYfcvNntCL1j7Vpy/HiTn/z8PfkKBMjLLzc3YlNGZaW5QV+fPuabX9NcDR1qbv1h90ZYfhUOk//4B3nyyXsfh6EQ2b07OXEiuX27/e3CYs9dGM+ZEgeUlJSwtLQ0rnVXrQLWrzePe/RIwFiVj2zdCixfbiY4depkbjeakeF2VN5UW2tmRVdUmB7XwQcDHTq4HZU3kcCyZeYa+IwMczvbXr3cjsq7NmwwM+fr6sx5rgED4j8fISLzSJbEbJeMxV0ppVKV1eLu+TF3pZRS9mlxV0opH9LirpRSPqTFXSmlfEiLu1JK+ZAWd6WU8iEt7kop5UNa3JVSyoe0uCullA9pcVdKKR/S4q6UUj5kqbiLyBki8o2IrBSR26K8PkpEykRkYWS5yvlQlVJKWRXz/oAikg7gSQCnwfxe6lwRmU7yP82aTiXp8A9FKaWUioeVnvsQACtJriZZB+BvAM5LbFhKKaXawkpx7wFgXZPn6yN/a+5CEVksIq+LSE9HolNKKRUXK8U92i3lm98EfgaAYpIDAfwTwAtRNyQyWkRKRaS0rKzMXqRKKaUss1Lc1wNo2hMvArChaQOS35OsjTx9GsDgaBsiOZlkCcmSgoKCeOJVSillgZXiPhdAfxH5gYhkARgBYHrTBiJS2OTpMABLnQtRKaWUXTGvliHZICJjAcwCkA7gOZJfi8g9MD/UOh3AjSIyDEADgG0ARiUwZqWUUjHob6gqpVQS0d9QVUqpFKbFXSmlfEiLu1JK+ZAWd6WU8iEt7kop5UNa3JVSyoe0uCullA9pcVdKKR/S4q6UUj6kxV0ppXxIi7tSSvlQzBuHedGGDcB335nHPXoA3bu7G4+XVVQAq1cDNTVAp05A//5Amn6kR9XQACxfDuzYAeTmAv36AcGg21F519q1wObNQHo60LMnsN9+bkfkXdu2AWvWAPX1QNeuQJ8+gET7pQwnkXRlGTx4MO1oaCDffJMcOpQMBMj8fLNkZ5NHHUVOm0bW19vapK8tWEBefvmeXHXoQObmkoWF5J/+RG7b5naE3rFpE3n33WSXLmQoZHKVn08Gg+S115LLlrkdoXfU1JAvvUQOGEDm5OzJVXY2ecop5HvvkeGw21F6QzhM/utf5LBhJj+7j8OcHLJvX/Kpp8jKSvvbhbkbb8wamxTFvbzcFPBQyEQcbQmFyJISLVrhMHn77WYHSk+PnqtgkOzYkfzyS7ejdd9775l9JxCInquMDJPLRx5xO1L3rV9P9uvX+nGYm0ueeSZZVeV2tO6qrydHjjT5EGk5V4WF5PLl9rbtm+JeVUUedpj55Gtph9q9ZGWZHsWuXfaS5Se33mp2mli52v2BuHix2xG755NPzAedlVwFg+Rjj7kdsXu2biWLisyHXaxc5eSYXnxDg9tRuyMcJkeMsLZvpaWZb4zr1lnfvtXi7vnR1/vvB775Bqitjd22rg5YuRK4557Ex+VFpaXA448Du3ZZa19ZCVxwgdnNUk19PfCTnwBVVdbaV1UBEyaYceZUdOONZny9oSF22+pq4N//Bp5+OvFxedFbbwEzZljbt8JhYPt24IornI/DUnEXkTNE5BsRWSkit0V5PVtEpkZenyMixU4EV18PPPGEORloVU0N8Oc/m0Kfav74R3u5AoCNG4HZsxMTj5dNn272LzvCYeDJJxMTj5eVlwPTptnLV1UV8OCDqdlx+MMfrHewAKCx0XwYrlnjbBwxi7uIpAN4EsCZAA4BcImIHNKs2S8AlJPsB2AigAecCG76dPM/blc4DLz5phMRJI/ycpOvcNjeelVVwEMPJSYmL3vwQWDnTnvr1NUBTz1l/0Mh2f3lL/FdYbVliylaqWT5cmDJEvvrJaLjYOWfbAiAlSRXk6wD8DcA5zVrcx6AFyKPXwdwikjbL/SZP9/+AQiY4Yb589v67sll+XIgO9v+eqQZzkk1//lPfOs1NppvO6nkiy+sD1811dgILFrkfDxetngxkJlpf726Ouc/CK0U9x4A1jV5vj7yt6htSDYA2AGgS/MNichoESkVkdKysrKYbxxPYd+toiL+dZNRPAffbtXVzsWRLKycw4kmLa1tuU5GdoYYmmpsTL1cVVXZ//a8W7x5bomV4h6tB958JM1KG5CcTLKEZElBQUHMN95vPzNBwi6R1JtQ0alT/DtVhw7OxpIMQqH41qurAzp2dDYWr7NwqEaVlZV6uerUKf5Jgp07OxuLlTDWA+jZ5HkRgA0ttRGRDAAdAGxra3BnnRXfUEMwCJxzTlvfPbkMGBDf18HsbODii52Px+vOPju+jkPv3sD++zsfj5cNHw7k5dlfr6EBOP105+PxsuOPj+9ijlAIuOQSZ2OxUtznAugvIj8QkSwAIwBMb9ZmOoCfRR4PB/Bh5HrMNhk0CCgutr9eURFw5JFtfffkkpFhLlcLBOytJwKMGZOYmLxs3Dj7HYfcXHM5ZMKnjXvM2WebXrhdxx9vjsVU0qEDcNFF9jsO4bALxT0yhj4WwCwASwG8SvJrEblHRIZFmj0LoIuIrAQwDsA+l0vG69577d3fIxg066Sia6+113sPBMw3nB7Nz6CkgCOOAA4/3F7Ryslx/gBMBunpwB132D8O77orcTF52a232us4BIPAddeZzoOjrMx0SsRi5/YDt99ubbZXMEjecovlzfrSxx9by1UgQA4cmNqzebduJXv3NjObW8uVCJmXRy5a5HbE7gmHyUsvtX4cPv642xG7a+pUM1PXSq5+/GN798WCX24/sNvEiSYR0abW5+aa1x56yNYmfWvOHLJ79+j3AMnKMoX9Jz+J76ZFfvP99+SJJ5oDsfnUehGzb/XvTy5d6nak7mtsNLe3CASiF/lQyHwITpnidqTe8Pe/k507m5xE61wFAuTo0fZveGi1uItp2/5KSkpYavMC6127gFdeMbNWN2wwaSosBMaOBS67LP4rIPwoHAY++MDMWl282Jzkyc83J0+vv96cGFR7fPMN8OijZtr4rl1myGroUOCWW4Cjj069cfbWlJebiU2TJwNlZWbYpndv4OabgQsvjO8iCL9qaDD71EMPAStWmAlwnTsDo0YBV18d38l5EZlHsiRmu2Qq7kopleqsFnfP3zhMKaWUfVrclVLKh7S4K6WUD7k25i4iZQDaenfsrgC2OhCO07wYlxdjAjQuO7wYE6Bx2eFETL1JxrwphGvF3QkiUmrlxEJ782JcXowJ0Ljs8GJMgMZlR3vGpMMySinlQ1rclVLKh5K9uE92O4AWeDEuL8YEaFx2eDEmQOOyo91iSuoxd6WUUtEle89dKaVUFElR3EXkDBH5RkRWisg+txMWkWwRmRp5fY6IFHsgplEiUiYiCyPLVe0Q03MiskVEvmrhdRGRxyIxLxaRQYmOyWJcJ4rIjia5SvjNYkWkp4h8JCJLReRrEbkpSpt2z5fFuNzIV0BEvhSRRZG47o7Spl2PQ4sxtftx2OS900VkgYi8E+W1xOfKyt3F3FwApANYBaAPgCwAiwAc0qzNdQAmRR6PADDVAzGNAvBEO+fqeACDAHzVwutnAXgX5mcRhwKY45G4TgTwTjvnqhDAoMjjPADLo/wbtnu+LMblRr4EQCjyOBPAHABDm7Vp7+PQSkztfhw2ee9xAF6J9m/VHrlKhp77EAArSa4mWQfgbwDOa9bmPAAvRB6/DuAUkYTex89KTO2O5Kdo/ecNzwPwIo3ZADqKSKEH4mp3JDeSnB95vBPmh2ia/2xJu+fLYlztLpKDysjTzMjS/IRdux6HFmNyhYgUATgbwDMtNEl4rpKhuPcAsK7J8/XYd2f/XxuaX47aAaCLyzEBwIWRr/Ovi0jPKK+3N6txu+HoyNfrd0Xk0PZ848hX4iNgen5NuZqvVuICXMhXZJhhIYAtAN4n2WK+2uk4tBIT4M5x+AiACQBa+tn6hOcqGYp7tE+z5p/OVto4ycr7zQBQTHIggH9iz6e0m9o7T1bNh5lSfRiAxwG81V5vLCIhAG8A+D+SFc1fjrJKu+QrRlyu5ItkI8nDARQBGCIiA5o1afd8WYip3Y9DETkHwBaS81prFuVvjuYqGYr7egBNP22LAGxoqY2IZADogMQOA8SMieT3JGsjT58GMDiB8VhlJZftjmTF7q/XJGcCyBSRrol+XxHJhCmgL5OcFqWJK/mKFZdb+Wry/tsBfAzgjGYvtfdxGDMml47DYwEME5E1MEO2J4vIlGZtEp6rZCjucwH0F5EfiEgWzMmH6c3aTAfws8jj4QA+ZORMhVsxNRubHQYzduq26QCuiFwFMhTADpIb3Q5KRLrtHm8UkSEw++X3CX5Pgflh96UkH26hWbvny0pcLuWrQEQ6Rh7nADgVwLJmzdr1OLQSkxvHIcnbSRaRLIapDR+SvLxZs4TnKsPJjSUCyQYRGQtgFsxVKs+R/FpE7oH5LcHpMAfDSyKyEubTb4QHYrpRRIYBaIjENCqRMQGAiPwV5kqKriKyHsBvYE4ygeQkADNhrgBZCaAKwM8THZPFuIYDGCMiDQCqAYxI8IczYHpXIwEsiYzZAsAdAHo1icuNfFmJy418FQJ4QUTSYT5MXiX5jpvHocWY2v04bEl750pnqCqllA8lw7CMUkopm7S4K6WUD2lxV0opH9LirpRSPqTFXSmlfEiLu1JK+ZAWd6WU8iEt7kop5UP/H31hrXWveiVWAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7f0ad1683bd0>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "# This code defines a 5x5 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": 3,
   "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": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAD8CAYAAACMwORRAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4xLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvAOZPmwAAIABJREFUeJzt3Xl8FPX9P/DXO5tjs9mEM0o4Uw5FpagQEbXetd5oFS0eWGoVL9RvqeDR1lZt/X498fxV8agHWrGKFiwWrWfVgoRbC3IJQrmCBELInX39/vhsSgyb7MxmNjM7+34+HvN47GY/M/vmzcx7P/uZ+cwKSSillPKXDLcDUEop5Twt7kop5UNa3JVSyoe0uCullA9pcVdKKR/S4q6UUj6kxV0ppXxIi7tSSvmQFnellPKhTLfeuHv37iwuLnbr7ZVSKiUtWLBgO8nCeO1cK+7FxcUoLS116+2VUiolich6K+10WEYppXxIi7tSSvmQa8MyiSCBTz8FnngCWLPG/K1/f+Cqq4BjjwVE3I3Pa9atA/7f/wPmzgWqq4Hu3YGLLwYuuAAIBt2OzlsqK4Fp04DXXwd27gTCYeD444GrrwZ69HA7Om+JRIB33gGeeQbYsAEIBIDBg4HrrgOGDXM7Ou/597+Bxx8HliwB6uqAoiJg3Djg7LOBzGRWYJKuLMOHD6cdf/87WVxM5uWRIqQp9eZxOEz260fOnm1rk761fj154olkMEhmZ+/NFWByFQ6Tv/0t2djodqTuq60lJ0wgc3PNvtU8V8EgmZNDjhpFlpW5Hak3TJtG7r+/2Yea5yojgwyFyIMPJv/1L7ej9IZly8jhw82+lZn53Xzl55Ndu5KPPUZGIva2C6CUFmpsShT35583CWqenFhLbi757LP2EuU3K1aYnSYQaDtXoRB57rlkQ4PbEbunpoY8+uj4+1ZWFtmrF/mf/7gdsbt+/3uz38Q7DkMh7Wh99tm+H4Ct5WrCBHsF3jfF/Z//tFbYmxf4Dz+0nig/qagge/T47jebeDvWxIluR+2eCy+0vm9lZpIHHEDW1bkdtTtefdXecRgKkV9+6XbU7tiwgSwosJerRx6xvn2rxd3yCVURCYjIIhF5K8ZrOSIyXURWi8g8ESl2atjollvMeLFV1dVmnXT0/PPA7t1ml7GiqsqMyZeXJzcuL1q7Fpg50/q+1dAAbNoEzJqV3Li8iAQmTbJ3HNbUAHfdlbyYvGzKFPPvt6qqCrj9dqC+3tk47FwtcyOA5a289nMA5SQHApgC4J72BgaYk6YLFthfb8kSYOVKJyJIHSRw333Anj321svIAJ57LikhedqjjwKNjfbWqawE7r03OfF42SefAN9+a2+dSAR488306zjU1ABPPWVOnNrR2Gg6G06yVNxFpDeAMwE83UqTcwA8H338GoCTRdp/7crrr9s/AAHzCfiXv7T33VPL8uX2D0DA9Bqeecb5eLzupZcS6yktWgTs2OF8PF724ov2Ow2AuRLkrX2+5/vbRx8ldtXe7t3As886G4vVnvtDACYDiLTyei8AGwCAZAOAXQC6tTe4zZsTOwCbvkKnk7KyxC+r2r7d2VhSQUVFYutlZ5tcp5P//Mf6UF9zdXXpl6uyMvOtJRFbtjgbS9ziLiJnAdhGsq0BklifVfvsDiIyXkRKRaS0zML/ek5O3CatSrfruLOy3Fk3VQUCia1HmgKfThL992ZkpF+usrISn2/jdK6s9NyPATBKRNYBeAXASSIyrUWbjQD6AICIZALoBGCfL68kp5IsIVlSWBj3vjcYPNhMJrErFDLrppP+/YHa2sTWPfBAZ2NJBX37JrZeJJJ+k5qGDk2so5WVBQwc6Hw8XjZoUGI990AAOOQQZ2OJW9xJ3kqyN8liAGMAvE/y0hbNZgL4afTx6GibBL7IfdcFFySWKBIYM6a9755aevQAjjnG/nr5+cDEic7H43UTJ9rvOGRmmv0qNzc5MXnVlVcmtl5ODnDKKc7G4nWHHw707Gl/vZwc4PrrnY0l4XvLiMidIjIq+vQZAN1EZDWAiQAcuRgxLw8YO9besEFmJnDRRaZopZvJk03O7MjNBU49NTnxeNnFF9vvOGRlAb/4RXLi8bLevYHjjjPDLFYFg8CNNyY+/JWqRICbb7Z/HA4cCBx6qMPBWLkYPhmL1UlM27aRRUVmenO8yQAZGWZq9JYt1icE+EkkQp53nvXJJrm55Jw5bkftnhdesDbjsmmiyXXXuR2xe1avJjt1sj7h68ADyd273Y7aHbW15JFHmltXWMlXXh65cKH17cPpSUxuKSw0Nwvr2bPtk6TBoBma+PRTYP/9Oy4+LxEBXn4Z+OEP2+45BALmvMQLLwA/+lHHxec1Y8eauQGhUNsnwfLyTE//kUc6LjavGTAA+OADoEuXtk/85eaaXuhHHyV2vswPsrPNjdUOO8zsW63JyjIjDLNnm+Ecx1n5BEjGYvfGYeXl5B/+QHbvbm6603QDrPx887e77iJ37LC1Sd9qbCT//Gfy0ENN7zw/3/QOCgrMzbAuuyx9p4bH8vnn5htPMGhylJdnchYMmnvPzJpl/+ZOfrV5Mzl5sslTU66ajsXevcmHHyb37HE7Sm+orSWffJIcMGDvPhUO783bhAnk11/b3y4s9tyF7T/vmZCSkhIm8ktMjY3Axx8DGzea5716mVuzptvYnlVffGGW6mrT6zrxRKBTJ7ej8qayMrNvVVSY3vqwYel3tYdVdXWmJ79liznPVVwMHH203nY7FhIoLQVWrTJ5Kyw0x2Fbvfq2iMgCkiVx26VacVdKqXRmtbh7fsxdKaWUfVrclVLKh7S4K6WUD2lxV0opH9LirpRSPqTFXSmlfEiLu1JK+ZAWd6WU8iEt7kop5UNa3JVSyoe0uCullA9Z+Q3VoIh8LiJLRORLEbkjRptxIlImIoujyxXJCVcppZQVmRba1AI4iWSliGQB+ERE3iY5t0W76SQnOB+iUkopu+IW9+j9gyujT7Oiizu3klRKKWWJpTF3EQmIyGIA2wC8S3JejGbni8hSEXlNRPo4GqVSSilbLBV3ko0kDwPQG8AIERnSosksAMUkhwL4B4DnY21HRMaLSKmIlJaVlbUnbqWUUm2wdbUMyZ0APgRwWou/f0uyNvr0KQDDW1l/KskSkiWFhYUJhKuUUsoKK1fLFIpI5+jjXAA/BLCiRZuiZk9HAVjuZJBKKaXssXK1TBGA50UkAPNh8CrJt0TkTpgfap0J4AYRGQWgAcAOAOOSFbBSSqn49DdUlVIqhehvqCqlVBrT4q6UUj6kxV0ppXxIi7tSSvmQFnellPIhLe5KKeVDWtyVUsqHtLgrpZQPaXFXSikf0uKulFI+pMVdKaV8yMqNwzynvh7YscM87toVyMpyNx4vi0SAnTuBmhqgc2cgFHI7Im+rrAR27QLy8oBOnQARtyPyrtpaoLwcCASALl2AzJSsJh2jsdHkqr7e5CoYTP57plTPfeFC4JJLgPx8oLjYLOEwcNFFwIIFbkfnLVu2AL/7HVBYCBQVAQMHmmI1YgTw+utAQ4PbEXpHTQ3w4ovAkCHmwBs0CNh/f6BnT+D++/d2JBRAAh9/DJx9tjkO+/cH+vYFCgqAq64CVqyIv4108vXXwC9+YTpWvXsDAwaYvJ18MjBnjul8JQ1JV5bhw4fTqpoacvRoMhQiAwHS7GJ7l4wM89qPf2zaprupU8lg0CwtcwWQ+flkv37kmjVuR+q+xYvJwkIyHI6dq1DILG+84Xak7tu1izz2WJMrkX1zlZlJ5uaS115LNjS4Ha27IhHy9tvNMZidHXvfCofJoUPJbdvsbRvmVutxa6zni3tDA3nKKWaniZWg5ktuLnnSSWR9vb1k+cnjj5tiFC9XGRlkt27k+vVuR+yeZctaL+qx9q0ZM9yO2D179pDf/z6ZkxM/V6EQOXasKXDpatIka8dhVpbpaO3YYX3bVou754dlpkwBPv0UqK6O37a6Gpg7F3jggeTH5UUrVgA33QRUVcVv2zQWf955yY/LiyIR4IwzgD17rLWvrgYuvRTYujW5cXnVpEnAqlVmnD2eqipgxgxg+vTkx+VF778PPP64teOwvh7YvBm44grn47DyM3tBEflcRJaIyJcickeMNjkiMl1EVovIPBEpdiK4SMSMeVpJUpOqKrNOY6MTEaSWBx80O4tVjY3A8uXA0qXJi8mr3n3XnOCijd+qiUSAp55KXkxetWcP8Nxz5tyEnXXuvjtpIXna3Xfbq1l1dcDs2c53HKz03GsBnETyUACHAThNREa2aPNzAOUkBwKYAuAeJ4J75x3rPavmamuBt992IoLUUVkJTJtm/0Rpba35dpRu7r3X5MyOmhrg4YfTr+Pw0kuJXTW0ejWweLHz8XjZhg3AJ58ktu6TTzobS9ziHh3maToMsqJLy/7OOQCejz5+DcDJIu2/iOyzz+wfgACwe7dZN50sX57YJaGNjcCHHzoejufNn5/YelVV5mt0OnnvvcQ6WQAwb56zsXjd/PlATo799WpqgH/8w9lYLI25i0hARBYD2AbgXZIt/8t6AdgAACQbAOwC0C3GdsaLSKmIlJaVlcV93507rUQXW7pdvlZZmfg12Xa+QvqFnSGG5gIB03lIJxUVia3X0JB+uaqsTPzyRqdzZam4k2wkeRiA3gBGiMiQFk1ilZV9RjNJTiVZQrKksLAw7vt27Zp4weq2z0eLvxUU2Bs/bi4cdjaWVJCbm9h69fUm1+mkS5fE1svKSr9cFRQAGQleptKpk7Ox2AqD5E4AHwI4rcVLGwH0AQARyQTQCUC7+84nnWRmCtoVDptJAunkkEMSK+5ZWcDppzsfj9edcEJiHYeuXc3kpnRy5pmJdQBI4PjjnY/Hy44+2pwgtSsUMhPDnGTlaplCEekcfZwL4IcAWs5Dmwngp9HHowG8H70es12OPTaxHnjnzsCJJ7b33VNLMGgup8rOtrdeIADceGNyYvKym26yfyuGUMisl263JBg9OrF/89ChwIEHOh+Pl+23H3DqqfbzFYkAl1/ubCxWeu5FAD4QkaUA5sOMub8lIneKyKhom2cAdBOR1QAmArjFieBEgF/9yt5BGAoBt92WfgcgANxwg737e2RnA0cdZabbp5sf/ADo08d8uFkVCAA/+1nyYvKqnBxgwgR7x2FeHvCb3yQvJi+79VZ7w37BIDBmTOLDX62yMtMpGYvVGaqRCHnxxdZme4VC5E9+kt4z41591dps3qwssk8fcvt2tyN2z7p1ZNeuZraulX3rgw/cjtg99fXk8cdb27dCITNDM5098IC1mhUMmlsQVFVZ3zb8MkNVxNzU6corzSdcrMuMsrPNa5dfnvg1uX5xwQVmZmBeXuxxUhHz2ve/D5SWpt+J5+b69TM3nOvfv/Ux5XDYDPPNmWPG6dNVZqbJwdlnm15prMtug0Gz3HYbcI8jM11S18SJZv5IMBj7G08gYP5+7LFmBn6iJ/jbZOUTIBmLnRuHNfnmG/KWW8w9UTIyzNK1Kzl5cnrfIyWWykry6afJAw7Ye7O13FzyvPPIzz5L7283LTU2ku++a+5h1HTvlMxM8rDDyJdf1pvRtfTVV+bmYAUF5gZiGRlkjx7k739Pbt3qdnTeUl5OPvig+ZackWHylZdHjhtHLlmS2DZhsecubP95z4SUlJSwtLQ04fWbriVN9LKjdEKafNkZX05njY2aK6siEfNtMJ2/LVvl1HEoIgtIlsRrl7K319eibp2IFis7NFfW6XFoXUcfh/pfo5RSPqTFXSmlfEiLu1JK+ZAWd6WU8iEt7kop5UNa3JVSyoe0uCullA9pcVdKKR/S4q6UUj6kxV0ppXxIi7tSSvmQlV9i6iMiH4jIchH5UkT2+d0eETlBRHaJyOLocntywlVKKWWFlRuHNQD4JcmFIpIPYIGIvEvy3y3a/ZPkWc6HqJRSyq64PXeSm0kujD7eDWA5gF7JDkwppVTibI25i0gxgMMBzIvx8lEiskRE3haRQxyITSmlVIIs389dRMIAXgfwPyQrWry8EEA/kpUicgaANwHs87PLIjIewHgA6Nu3b8JBK6WUapulnruIZMEU9pdIzmj5OskKkpXRx7MBZIlI9xjtppIsIVlSWFjYztCVUkq1xsrVMgLgGQDLST7YSpse0XYQkRHR7X7rZKBKKaWsszIscwyAsQCWicji6N9uA9AXAEg+AWA0gGtEpAFANYAxdOvHWZVSSsUv7iQ/AdDmz9+SfAzAY04FpZRSqn10hqpSSvmQFnellPIhLe5KKeVDWtyVUsqHtLgrpZQPaXFXSikf0uKulFI+pMVdKaV8SIu7Ukr5kBZ3pZTyIS3uSinlQ5bv5+4Ve/YAf/0rsHGjed6rF3DOOUA47G5cXhSJAO+/DyxbBlRXA126AGecAfTr53Zk3rRiBfDee8CuXUBeHlBSAhx9NCBt3lkpPZWXm+Nw61YgIBEU98/A2WcDOTluR+Y9DQ3A228DK1cCdXVAYSFw9tnA/vsn+Y1JurIMHz6cdmzYQF51FRkKkfn5ZGamWcJhMjeXvPJK8ptvbG3St2pryfvvJ3v0MPnJziYzMkyegkHyxBPJjz92O0rv+NvfyCOO2JufjAwyJ4fMyyOLi8knniAbG92O0htWrCAvusjkKS+PLMzYzo3oyb7hb5mfT/7yl+T27W5H6Q0VFeRvfkN26WJqVlYWGQiYGhYMkqNGkYsW2d8ugFJaqLEpUdwXLzYJysw0EcdaMjPJzp3JhQvtJ8tPKirII480haq1XAFmB3vsMbejdd+vfmVyES9XZ5xhPjTT2fvvm4KekbE3N+PwLAnwp/gTAdOR6NmTXLPG7WjdtWULOWiQ6SS0tl+JmH3r9dftbds3xX39elPY2zr4mi+dO5Nff20vWX7R0EAed1zbO1TzJTeXfOUVt6N2zwMPxC/szXN1wQVkJOJ21O5YvNgU9pZ5+RDHkQA/wPH//VtGBtmrV/r24PfsIQcPbrsz2nLf+vBD69u3Wtyt/BJTHxH5QESWi8iXInJjjDYiIo+IyGoRWSoiw5waNrr1VqCi5S+2tqGiArj5ZqfePbXMmgUsXAjU1lprX10NXH01UF+f3Li8aOdO4Fe/AqqqrLWvrgb+9jfg88+TG5dXXXONOd/VXD4qcCTmAgBGYi7yYQ7USAQoKwPuu6+jo/SGp58GvvnGjLVbUV0NXHGFKfVOsnK1TAOAX5I8CMBIANeJyMEt2pwO84PYg2B+APuPTgRXXg7MmAE0NlpfJxIBZs4Evk3DH/m75x6gstLeOo2NwJtvJiceL3vuOSDD5rViNTXA/fcnJRxPW7UK+HphOYZiyXeWK/A06mDOoNYhG1fg6f++NrhuCT55fAnqtpa7HH3HIs2HmtVOQ5PNm4F58xwPxt5wCoC/Ajilxd+eBHBRs+dfAShqaztWhmWeeML61+aWX3PSbTz5m2/MSRq7uQLIY45xO/qO179/YrnKySErK92OvmNNnkzeK5NJgFUIcicK/rs0T07zv1fB7Iwrz5vsdvgd6rPPzEUMdverjAzy0kutvQcsDsvYuhRSRIoBHA6g5WdMLwAbmj3fGP3bZtufNs2sXm3/ExAwX3NWrWrPO6ee9evNZWg1NfbXXbfO8XA8b9OmxNbLygK2bAEGDHA2Hi9bsQKYxf/FDnTGb3AXOiH2OGnzv1chF78O3I39jr0ZgzoqUA9Yty6xS2cjEXOppJMsfzEVkTCA1wH8D8mW/7ux/jn7jCCJyHgRKRWR0rKysrjvaWc4xsl1U1Ek4s66qao945vplq/GRoDIwP/hVpyAD7EZPVCNYMy21cjBJhTheHyE+wK3IpJm8yQjkcT3LadrlqXMi0gWTGF/ieSMGE02AujT7HlvAPv0jUhOJVlCsqSwsDDu+/bpk9ikiOxsoG9f++ulsqKixE+MJn0yhQd17ZrYek2TUNJJcfHe8xPzMQIH4itsRlHMtpvRE4OxAqU4Ajk5QI8eHRenFxQV2T+X06RPn/ht7LBytYwAeAbAcpIPttJsJoDLolfNjASwi2S7hmQA4MILE/uKk5EB/OQn7X331DJoUGI7R16euWIm3Vx+eWIdh2OOATp3dj4eLxs3Dgg266jXIRtFrYy49sAW1EZPsjY0AGed1QEBeshxxwGBgP31wmFg/HhnY7HyGXMMgLEAThKRxdHlDBG5WkSaysJsAGsBrAbwFIBrnQiuVy/g+OPtr3fMMenXcwfMJaB5efbWIYFLLklOPF52zTX21wmHgcmTnY/F60pKvns8nYJ3UYdsAMAehFCLbOxBCABQjyycgneRmWn2q3S7LUhmJjBhgv2OQygEnHqqw8FYOeuajMXqJKZFi+xdMRMKkQsWWDvr7DdVVeSAAdYnT4RC5B/+4HbU7rnmGuv7Vna2uUVBut6GYM6cvbOep+MCRgBWIsSHcD3zsJsPYwIrkcsIwFdwIQsKyLVr3Y7aHdu2kYWFZgaq1ePwpZesbx9+maFKkrNnWzsIc3PJWbOsJ8mPNm4ke/c297GIt0NdfXX6zrgkzYzec8+Nv28Fg+RBB5E7drgdsbuefposCNayEiHuRAF/hL9/J0+n4m3uRAErkcd/fVzndriu+uILsmtXcy+ZeMfh//6vvW37qriT5Pz55FFHmQOteeHKyjJF/cgjyblz7SXJr7ZvJ8eONblqWbjCYXPvj6lT07uwN2lsJO+5h+zWzdzcqXmu8vJM/q69Nv2ubW/Ne29V8c9dr2Fx7pbv3GMmJ8fsb6OP3cIdY64xXyPT3Pr15uZgweB356CImH1rwAByxgz727Va3MW07XglJSUsLS21vd6qVcCf/gSsWWNS1b8/8LOfAQcemIQgU9zOncALLwDz55up4/vtB4weDZx0UuJn9P2qocHcXmDWLDO7uaDAnBwbM8b+eYx0sHix2be++cZc+z94sDlJ7fQVH36wdaupWcuWmXkoPXua8xFHHpnYBSMisoBkSdx2qVbclVIqnVkt7tp/U0opH9LirpRSPqTFXSmlfEiLu1JK+ZAWd6WU8iEt7kop5UNa3JVSyoe0uCullA9pcVdKKR/S4q6UUj6kxV0ppXxIi7tSSvmQlZ/Ze1ZEtonIF628foKI7Gr2K023Ox+mUkopOzIttHkOwGMAXmijzT9JptmvJSqllHfF7bmT/BjAjg6IRSmllEOcGnM/SkSWiMjbInJIa41EZLyIlIpIaVlZmUNvrZRSqiUnivtCAP1IHgrgUQBvttaQ5FSSJSRLCgsLHXhrpZRSsbS7uJOsIFkZfTwbQJaIdG93ZEoppRLW7uIuIj1EzC8BisiI6Da/be92lVJKJS7u1TIi8mcAJwDoLiIbAfwWQBYAkHwCwGgA14hIA4BqAGPo1g+zKqWUAmChuJO8KM7rj8FcKqmUUsojdIaqUkr5kBZ3pZTyIS3uSinlQ1rclVLKh7S4K6WUD2lxV0opH9LirpRSPqTFXSmlfEiLu1JK+ZAWd6WU8iEt7kop5UMpVdxra4GXXwaOOgro3dssI0cC06YBNTVuR+ctJPDpp8B55wHf+x7Qsydw0EHAnXcCW7e6HZ33rF8PTJoEHHAAUFQE9O8PjB0LLFrkdmTeU1kJPPkkMGwY0KsX0LcvcMIJwJtvAg0NbkfnLZEIMGcOcNppQL9+Jl9DhwJTpgA7dyb5zUm6sgwfPpx2PP00mZ9vFlO69i7hsPn7E0/Y2qRvLVpEDhhA5uWRIt/NVTBI5uSQY8eS1dVuR+q+XbvIM880ecnO/m6uAgEyFCKHDiXXrHE7UvdFIuRdd5mc5OXtexzm55Ndu5JvvOF2pN7w/vtkUVHsmhUKmX1u4kSyocHedgGU0kKNTYni/vvfm2S0TFCshP3ud/YS5Tdz58Y+8FouubnkiBHpXeDLy8lBg8yHXVu5ysggO3cmly93O2L3RCLklVdaOw5zc8lnnnE7YnfNnGm9Zp1zDtnYaH3bvinus2ZZS1LzZL35pvVE+UlZGdmpk/Vc5eaSl1zidtTuOfHE+IW9aREhe/Qgq6rcjtodjz9u/zicO9ftqN3x1VfWOljNc/XrX1vfvtXiHnfMXUSeFZFtIvJFK6+LiDwiIqtFZKmIDHNy2OjXvwaqqqy3r6oy66SjqVOBujrr7aurgddfB7ZsSV5MXrVsGTB3rjmPYwVpxpr/8pfkxuVFjY3AHXfYOw6rq835nXR07732zgFWVQEPPWRy5iQrJ1SfA3BaG6+fDmBQdBkP4I/tD8tYuhRYtcr+emvXpt+JsMbGxHeQJ590Ph6vmzLF3gchYIr7PfckJx4vmzPH/n5FAu+9B2zalJyYvGr3bnPRR2Oj/XWnT3c2lrjFneTHAHa00eQcAC9EvzHMBdBZRIqcCO6tt+wfgID51Jw504kIUscXXyRW2GtqzM6Ybv7618QOwDVrgG3bnI/Hy/7yF1O07MrMNB8M6eSjj4CsLPvrVVY6fxw6cSlkLwAbmj3fGP1bu23dmtilVZFI+g01lJcDgUBi6+7a5WwsqaCyMrH1srOBHW11dXwo0WOpri79clVebupPIrZvdzYWJ4q7xPhbzB/IFpHxIlIqIqVlZWVxN5yXl3hQ+fmJr5uKgsHE183JcS6OVJFI7wowB25urrOxeF0olNh6gUD65SoYBDISrKpO58qJ4r4RQJ9mz3sDiDnSRnIqyRKSJYWFhXE3fOihQDhsP6Bw2KybTgYNSmwilwhw2GHOx+N1Bx6Y2HoiZpJTOjnyyMQKT2YmMGSI8/F42ZAhQH29/fWysoAjjnA2FieK+0wAl0WvmhkJYBfJzQ5sF+eem/in4PnnOxFB6ujWDTjjDFN87AiFgF/+MjkxedmkSfa/3WVlAT//uRmaSSeXX25OkNrVpQtw7LHOx+NlBx0EDB5sf71AAJgwwdlYrFwK+WcA/wJwoIhsFJGfi8jVInJ1tMlsAGsBrAbwFIBrnQouJwe46ip7Qw45OcCVV7ZvmCJVTZpkv4dVWJh+ByBgbstgt+MQCADXX5+ceLyse3fgzDPtDWWFQmZ/tNvZ8INbbrE3pJyRAZSUAAMHOhyIlYvhk7FYncS0ezc5eDCZlRV/MkBWFnnAAWRFhfUJAX5z/fXWJ5vk5ZELFrgdsXsZ0FfSAAAOj0lEQVTeecdM5LI60eSee9yO2D1btphJXBkZ8XMVDJI/+AFZV+d21O5obDSzTq3sWyJm9vPatda3D6cmMbktHAY++QQ4+OC2Pw3z8sw46iefpN/J1OYeesgMHYRCrfeagkGgoMBcpjbM0SlnqeWUU4BXXzX7TmtDLU0nBX/9a9MTTVf7728mffXu3fZxGA6bG/u9/XbiJ61TXUaG2a9OP73tXIVC5lvRp5+am/s5zsonQDIWuzcOq60lX3rJ3MQpN9dMs+/UyTweMoScNo2sqbG1SV/79FPy3HPN9PqmXOXnk926kXfeSW7d6naE3rF+PTlpkslPQcHefAWD5KWXmhuxKaOy0tygr39/882vea5GjjS3/rB7Iyy/ikTIv/+dPOmk7x6H4TDZsyc5ZQq5c6f97cJiz12YyJkSB5SUlLC0tDShddesATZuNI979UrCWJWPbN8OrFxpJjh16WJuN5qZ6XZU3lRba2ZFV1SYHtdBBwGdOrkdlTeRwIoV5hr4zExzO9u+fd2Oyrs2bTIz5+vqzHmuIUMSPx8hIgtIlsRtl4rFXSml0pXV4u75MXellFL2aXFXSikf0uKulFI+pMVdKaV8SIu7Ukr5kBZ3pZTyIS3uSinlQ1rclVLKh7S4K6WUD2lxV0opH9LirpRSPmSpuIvIaSLylYisFpFbYrw+TkTKRGRxdLnC+VCVUkpZFff+gCISAPA4gFNgfi91vojMJPnvFk2nk3T4h6KUUkolwkrPfQSA1STXkqwD8AqAc5IbllJKqfawUtx7AdjQ7PnG6N9aOl9ElorIayLSx5HolFJKJcRKcY91S/mWN4GfBaCY5FAA/wDwfMwNiYwXkVIRKS0rK7MXqVJKKcusFPeNAJr3xHsD2NS8AclvSdZGnz4FYHisDZGcSrKEZElhYWEi8SqllLLASnGfD2CQiHxPRLIBjAEws3kDESlq9nQUgOXOhaiUUsquuFfLkGwQkQkA5gAIAHiW5JcicifMD7XOBHCDiIwC0ABgB4BxSYxZKaVUHPobqkoplUL0N1SVUiqNaXFXSikf0uKulFI+pMVdKaV8SIu7Ukr5kBZ3pZTyIS3uSinlQ1rclVLKh7S4K6WUD2lxV0opH9LirpRSPhT3xmFetGkT8J//mMe9egE9e7obj5dVVABr1wI1NUCXLsCgQUCGfqTH1NAArFwJ7NoF5OUBAwcCoZDbUXnX+vXA1q1AIAD06QPst5/bEXnXjh3AunVAfT3QvTvQvz8gsX4pw0kkXVmGDx9OOxoayDfeIEeOJINBsqDALDk55JFHkjNmkPX1tjbpa4sWkZdeujdXnTqReXlkURF5//3kjh1uR+gdW7aQd9xBdutGhsMmVwUFZChEXn01uWKF2xF6R00N+eKL5JAhZG7u3lzl5JAnn0y+8w4ZibgdpTdEIuQ//0mOGmXy03Qc5uaSAwaQTz5JVlba3y7M3Xjj1tiUKO7l5aaAh8Mm4lhLOEyWlGjRikTIW281O1AgEDtXoRDZuTP5+eduR+u+d94x+04wGDtXmZkmlw895Hak7tu4kRw4sO3jMC+PPP10sqrK7WjdVV9Pjh1r8iHSeq6KisiVK+1t2zfFvaqKPPRQ88nX2g7VtGRnmx7Fnj32kuUnN99sdpp4uWr6QFy61O2I3fPRR+aDzkquQiHykUfcjtg927eTvXubD7t4ucrNNb34hga3o3ZHJEKOGWNt38rIMN8YN2ywvn2rxd3zo6933w189RVQWxu/bV0dsHo1cOedyY/Li0pLgUcfBfbssda+shI47zyzm6Wb+nrgxz8Gqqqsta+qAiZPNuPM6eiGG8z4ekND/LbV1cC//gU89VTy4/KiN98EZs2ytm9FIsDOncBllzkfh6XiLiKnichXIrJaRG6J8XqOiEyPvj5PRIqdCK6+HnjsMXMy0KqaGuCPfzSFPt3cd5+9XAHA5s3A3LnJicfLZs40+5cdkQjw+OPJicfLysuBGTPs5auqCrj33vTsOPzf/1nvYAFAY6P5MFy3ztk44hZ3EQkAeBzA6QAOBnCRiBzcotnPAZSTHAhgCoB7nAhu5kzzD7crEgHeeMOJCFJHebnJVyRib72qKuCBB5ITk5fdey+we7e9derqgCeftP+hkOr+9KfErrDats0UrXSyciWwbJn99ZLRcbDyXzYCwGqSa0nWAXgFwDkt2pwD4Pno49cAnCzS/gt9Fi60fwACZrhh4cL2vntqWbkSyMmxvx5phnPSzb//ndh6jY3m2046+ewz68NXzTU2AkuWOB+Ply1dCmRl2V+vrs75D0Irxb0XgA3Nnm+M/i1mG5INAHYB6NZyQyIyXkRKRaS0rKws7hsnUtibVFQkvm4qSuTga1Jd7VwcqcLKOZxYMjLal+tUZGeIobnGxvTLVVWV/W/PTRLNc2usFPdYPfCWI2lW2oDkVJIlJEsKCwvjvvF++5kJEnaJpN+Eii5dEt+pOnVyNpZUEA4ntl5dHdC5s7OxeJ2FQzWm7Oz0y1WXLolPEuza1dlYrISxEUCfZs97A9jUWhsRyQTQCcCO9gZ3xhmJDTWEQsBZZ7X33VPLkCGJfR3MyQEuvND5eLzuzDMT6zj06wfsv7/z8XjZ6NFAfr799RoagFNPdT4eLzvuuMQu5giHgYsucjYWK8V9PoBBIvI9EckGMAbAzBZtZgL4afTxaADvR6/HbJdhw4DiYvvr9e4NHHFEe989tWRmmsvVgkF764kA11yTnJi8bOJE+x2HvDxzOWTSp417zJlnml64XccdZ47FdNKpE3DBBfY7DpGIC8U9OoY+AcAcAMsBvErySxG5U0RGRZs9A6CbiKwGMBHAPpdLJuquu+zd3yMUMuuko6uvttd7DwbNN5xeLc+gpIHDDwcOO8xe0crNdf4ATAWBAHDbbfaPw9tvT15MXnbzzfY6DqEQcO21pvPgKCsznZKx2Ln9wK23WpvtFQqRN91kebO+9OGH1nIVDJJDh6b3bN7t28l+/czM5rZyJULm55NLlrgdsXsiEfLii60fh48+6nbE7po+3czUtZKrH/3I3n2x4JfbDzSZMsUkItbU+rw889oDD9japG/Nm0f27Bn7HiDZ2aaw//jHid20yG++/ZY84QRzILacWi9i9q1Bg8jly92O1H2Njeb2FsFg7CIfDpsPwWnT3I7UG/72N7JrV5OTWJ2rYJAcP97+DQ+tFncxbTteSUkJS21eYL1nD/Dyy2bW6qZNJk1FRcCECcAllyR+BYQfRSLAe++ZWatLl5qTPAUF5uTpddeZE4Nqr6++Ah5+2Ewb37PHDFmNHAncdBNw1FHpN87elvJyM7Fp6lSgrMwM2/TrB/ziF8D55yd2EYRfNTSYfeqBB4BVq8wEuK5dgXHjgCuvTOzkvIgsIFkSt10qFXellEp3Vou7528cppRSyj4t7kop5UNa3JVSyodcG3MXkTIA7b07dncA2x0Ix2lejMuLMQEalx1ejAnQuOxwIqZ+JOPeFMK14u4EESm1cmKho3kxLi/GBGhcdngxJkDjsqMjY9JhGaWU8iEt7kop5UOpXtynuh1AK7wYlxdjAjQuO7wYE6Bx2dFhMaX0mLtSSqnYUr3nrpRSKoaUKO4icpqIfCUiq0Vkn9sJi0iOiEyPvj5PRIo9ENM4ESkTkcXR5YoOiOlZEdkmIl+08rqIyCPRmJeKyLBkx2QxrhNEZFezXCX9ZrEi0kdEPhCR5SLypYjcGKNNh+fLYlxu5CsoIp+LyJJoXHfEaNOhx6HFmDr8OGz23gERWSQib8V4Lfm5snJ3MTcXAAEAawD0B5ANYAmAg1u0uRbAE9HHYwBM90BM4wA81sG5Og7AMABftPL6GQDehvlZxJEA5nkkrhMAvNXBuSoCMCz6OB/Ayhj/hx2eL4txuZEvARCOPs4CMA/AyBZtOvo4tBJThx+Hzd57IoCXY/1fdUSuUqHnPgLAapJrSdYBeAXAOS3anAPg+ejj1wCcLJLU+/hZianDkfwYbf+84TkAXqAxF0BnESnyQFwdjuRmkgujj3fD/BBNy58t6fB8WYyrw0VzUBl9mhVdWp6w69Dj0GJMrhCR3gDOBPB0K02SnqtUKO69AGxo9nwj9t3Z/9uG5pejdgHo5nJMAHB+9Ov8ayLSJ8brHc1q3G44Kvr1+m0ROaQj3zj6lfhwmJ5fc67mq424ABfyFR1mWAxgG4B3Sbaarw46Dq3EBLhzHD4EYDKA1n62Pum5SoXiHuvTrOWns5U2TrLyfrMAFJMcCuAf2Psp7aaOzpNVC2GmVB8K4FEAb3bUG4tIGMDrAP6HZEXLl2Os0iH5ihOXK/ki2UjyMAC9AYwQkSEtmnR4vizE1OHHoYicBWAbyQVtNYvxN0dzlQrFfSOA5p+2vQFsaq2NiGQC6ITkDgPEjYnktyRro0+fAjA8ifFYZSWXHY5kRdPXa5KzAWSJSPdkv6+IZMEU0JdIzojRxJV8xYvLrXw1e/+dAD4EcFqLlzr6OIwbk0vH4TEARonIOpgh25NEZFqLNknPVSoU9/kABonI90QkG+bkw8wWbWYC+Gn08WgA7zN6psKtmFqMzY6CGTt120wAl0WvAhkJYBfJzW4HJSI9msYbRWQEzH75bZLfU2B+2H05yQdbadbh+bISl0v5KhSRztHHuQB+CGBFi2YdehxaicmN45DkrSR7kyyGqQ3vk7y0RbOk5yrTyY0lA8kGEZkAYA7MVSrPkvxSRO6E+S3BmTAHw4sishrm02+MB2K6QURGAWiIxjQumTEBgIj8GeZKiu4ishHAb2FOMoHkEwBmw1wBshpAFYCfJTsmi3GNBnCNiDQAqAYwJskfzoDpXY0FsCw6ZgsAtwHo2ywuN/JlJS438lUE4HkRCcB8mLxK8i03j0OLMXX4cdiajs6VzlBVSikfSoVhGaWUUjZpcVdKKR/S4q6UUj6kxV0ppXxIi7tSSvmQFnellPIhLe5KKeVDWtyVUsqH/j/PtaazA2j0JQAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7f0ad1683750>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "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 (remember - you can see other files in the current directory by clicking on the `jupyter` logo in the top left of this notebook).\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": 4,
   "metadata": {},
   "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": 5,
   "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": 11,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Tests pass! Your sense function is working as expected\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": 12,
   "metadata": {},
   "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": 24,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAW4AAAD8CAYAAABXe05zAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4xLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvAOZPmwAAIABJREFUeJzt3Xl4lOW9xvHvj4SQBXApFEFF1CpuuBG1Kq5o1WLdqFZrqZWjnGq19nJBabXH5RyLvWrVo9VTBXGjWk8rYl1Qq4iigCSuuHGEYgURUEFMAll/548nlMWEzMw7M++84f5c11wJycy89zCZe5553s3cHRERSY4ucQcQEZH0qLhFRBJGxS0ikjAqbhGRhFFxi4gkjIpbRCRhVNwiIgmj4hYRSRgVt4hIwhTn4k579erlAwYMyMVdi4h0StXV1Z+5e+9UrpuT4h4wYABVVVW5uGsRkU7JzD5K9bqaKhERSRgVt4hIwqi4RUQSRsUtIpIwKm4RkYRJqbjNbHMz+4uZvW9m75nZgbkOJiIibUt1c8BbgCnu/n0zKwHKc5hJREQ2osMRt5n1BA4FxgO4e4O7r8hFmLo6ePxx6AxnU5s8Gf72t7hTRDdpEvTqBddcE3eS6ObPh8ceiztFdO5w003wyitxJ4lu6VK4+WZYvTruJNEtWwZLluRnWalMlewALAMmmNnrZjbOzCo2vJKZjTKzKjOrWrZsWUZhXnkFfvUr+PLLjG5eUC6+GC65JO4U0d13H3z+OYwfH3eS6O66KzwvSbdyZfjbGjs27iTRPfEEXHYZzJkTd5LoZs+GWbPysyzr6GTBZlYJzAQOdvdZZnYLsNLdr2rvNpWVlZ7pnpO1tVDxtbeF5Fm0CMygX7+4k0Qzdy5cdx2MHAlHHBF3mmgaG8OgoFevuJNEN2dO+Nvacsu4k0TT3Bz+xnbdNe4k8TOzanevTOm6KRT3VsBMdx/Q+u9DgCvcfVh7t4lS3CIim6J0irvDqRJ3/xT42MwGtv5oKPBuhHwiIhJBqluVXAhMbN2iZD5wdu4iiYjIxqRU3O7+BpDSEF5ERHJLe06KiCSMiltEJGFU3CIiCaPiFhFJGBW3iEjCqLhFRBJGxS0ikjAqbhGRhFFxi4gkjIpbRCRhVNwiIgmj4hYRSRgVt4hIwqi4RUQSRsUtIpIwKm4RkYRRcYuIJIyKW0QkYVTcIiIJo+IWEUkYFbeISMKouEVEEkbFLSKSMCpuEZGEKU7lSma2APgKaAaa3L0yl6GSbtEieOQRMINTToF+/eJOJCKdSToj7iPcfe9clnZTEyxcCO65WkLuzZ0Lu+8Oo0fDZZeF7+fNiztVZmpr4ZhjoKQE9tgD/vnPuBOJCBTYVMmQIbD99qH0kuqXv4SvvoLVq8Nl5Uq48sq4U2Vm7FiYNg0aG+G99+Dss+NOFM1pp8HWW8OHH8adJJolS8Ib6YUXxp0kulmz4Ljj4OOP404S3c03ww035GdZqRa3A8+YWbWZjWrrCmY2ysyqzKxq2bJlGYV5/fUw6p4+PaObF4QlS6ClZe2/W1rCz5JowQKorw/ft7Qkf8T99NOwdCm8/XbcSaKZOxfeeQcmTYo7SXR//ztMmQJvvhl3kujuuANuvz0/y0q1uA92932B44CfmdmhG17B3e9090p3r+zdu3dGYSZNgjPPhPHjM7p5QTj5ZCgvX/vv8vLwsyQ699yQv6QEKiqSP8KbNg3uuw9OPDHuJNEMGRLehKZNiztJdJdfDnPmwLBhcSeJbvZseOut/CzLPM0JZTO7Gqhx99+1d53KykqvqqqKGC2ZWlrgiivCu68ZXHAB/Nd/he+T6LXX4PnnYdCgMN8tIrlhZtWprkPssLjNrALo4u5ftX7/LHCtu09p7zabcnGLiGQineJOZXPAPsAkC0PGYuBPGyttERHJrQ6L293nA3vlIYuIiKSgoDYHFBGRjqm4RUQSRsUtIpIwKm4RkYRRcYuIJIyKW0QkYVTcIiIJo+IWEUkYFbeISMKouEVEEkbFLSKSMCpuEZGEUXGLiCSMiltEJGFU3CIiCaPiFhFJGBW3iEjCqLhFRBJGxS0ikjAqbhGRhEnlLO8iUiA+/xz+8hdYtAhKSmDwYPjOd6CoKO5kkk8qbpEE+OwzuOACmDwZunSBurrwtbw8XH79azj/fDCLO6nkQ8EUtzu89hosXQo77QTf+lbciWSNxkYoLlYpxGXxYthvv/DaaGxc+/OWFqipCZfLL4e334Y77tDztClIeY7bzIrM7HUzezzbIaZNg+23h8MPhzPOgD33hAMPhI8+yvaSJFXNzXDjjbDVVtCtG/TsCb/4BXz1VdzJMjNnDlx2GYwcCQ8/vH4BFjJ3OO44WLJk45lra+H++2H8+Pxli6qxER59FG66CaZODY9VUuTuKV2Ai4E/AY93dN3Bgwd7qqqq3MvL3cPTtvZSVOT+zW+6f/FFynclWXTGGV9/Xrp1cx80yH3VqrjTpefGG93LysLfFLh37x4ex5dfxp2sYy+/7F5R8fXXR3uXbbd1b2mJO3XH5s9379fPvUcP95KS8Jzsu28ynpNcAao8xT5OacRtZtsAw4Bx2X7jGDMmzNdtqLk5jO7GZX2J+bF4MXz6adwpMvPGG2EudcPnpb4e5s+HBx+MJ1cm5s2DX/0KVq0Kf1MQphbmzoVrr403Wyp+//u2Xx/tWb4cXnwxd3my5bTTwuvjq6+goSE8J++8Ez4VJdW774bpqnxIdarkZmA00JLNhTc3w/PPt//7VavgvvuyucT8ePfdMEe/447wwQdxp0nfww/D6tVt/662Fu6+O795ovjzn9cW9rrq6+Hee/OfJ13V1elNITQ1wVtv5S5PNixaFKauWjZok/p6mDgxnkxRzZgBlZVwwAFh2ifXOixuMzseWOru1R1cb5SZVZlZ1bJly1JaeHPz15+8DdXXp3RXBWXp0rXfp/hfUVDq6jb+vKxalb8sUdXWhjJrS3tvToWkrTedjWlpaf/xFoq6uvY3X0zi6x3CJ+w1K4U/+ST3y0tlxH0wcIKZLQAeAo40swc2vJK73+nule5e2bt375QWXlICAwe2//viYhg6NKW7KiiHHRZGc/ffD0OGxJ0mfUcfDT16tP270lIYNiy/eaI49tiwudyGunQJ2z8Xuv7907t+SUn6t8m3HXcMK7s3ZAZHHJH/PNlw0klwww1w/fVw+um5X555Gp/DzOxw4FJ3P35j16usrPSqqqqU7vPhh+Hss9uexysvD/OtO+2UckTJguZm2H33MJ+94ZYMm20Wpn/69IknW7rc4fjj4YUX1v6NFRVBRQXMng077xxrvDa9/nqYp37xxbDF1eefp3f7QYPgkEPgoIPgmGOgV6/c5IxiyhQYPjyMsJubwxtOaSnMnAm77hp3uniYWbW7V6Z05VTXYrYW/OFkeasSd/err3YvLQ1rlyGsRa+ocH/iiQxWzUpWLFnifthh4XnZbLPwfOy8s/tbb8WdLH2Nje633OI+cKD7Vlu5jxjh/uGHcadaX22t+913h4wVFeH/PdUtSdq6mIUtNUpL3YcPd58xo/C2NnnnHfdzznE/5BD30aPdFy6MO1G8SGOrkrRG3KlKZ8S9xsKFYcXEJ5+E0d4ZZ7T/cV3yZ/78sAVG375h+3rt3JFd//d/8N//DRMmhP/bmprsL6NLFygrC8/h6NHwwx+GTxxSWNIZcRdMcYtsSlasgPPOC5tdNjXlb4eg7t3DG8Rtt8GIEXojLiTpFLeODiiSZ08+GVbQTZoUttDJ516cNTVh2+nzz4ejjgpbQ0jyqLhF8mTFirDFwamnwhdfxLvpW21tWPm5885hXwntbp4sKm6RPHjxxTDKfvTR9PaEzKWmpjACP//8sGnkypVxJ5JUqbhFcuxvfwsHiop7lN2e2lp46SXYf/9k7jC2KVJxi+TQxInwgx8Uzii7PWuOQ1NZmZ89/yQaFbdIjjz6KJx7bnIOEdDYGEr7oIPS3+lH8kvFLZIDzz0HZ56ZnNJeo6kpbGlyyCGa8y5kKm6RLFu8GE4+ufCnR9rT0BCmTUaOjDuJtEfFLZJF7vDjHyfjyIMbU18PTz0VVqxK4VFxi2TRxInh2MxJOTXaxtTVwVlnha1hpLCouEWyZPHisE10bW3cSbKnrg5GjYo7hWxIxS2SBZ1limRDa6ZMHnss7iSyLhW3SBZMndp5pkg2tGbU3dHZqiR/VNwiWfDb33auKZIN1dXB00/HnULWUHGLRLRwYThTTWf21VfhzUkKg4pbJKLbb980jq43cyb84x9xpxBQcYtE0tgYirsQDx6VbS0t4QQMEj8Vt0gEkydvOivtGhrgrrs2jTepQqfiFongscfC/O+mwgzeeCPuFKLilk6tpgb++EfYdddw8umyMvjmN+GnP4X3349+/6+8Ev0+MmXkf6jf2AjV1XlfrGxAxS2dkjvceCP06QOXXhpKuqYm7CCzbBmMHw/77AOHHw6ffZbZMurr4Z//zGrslG3J53zMtmxBfvdHX7UqnHRB4qXilk7p4ovh178O2x/X1Hz9901NocRnzAgFvmRJ+suYMwdKS6NnzcQJPMbWfMIJ5H+Xxlmz8r5I2YCKWzqdu+6CO+9M7bCqDQ2htI86Kv2VjNXV0NycWcaofsI9633Np48/jn6c8YYG+POf4bDDwjTW8cfDs89uGptVZkOHxW1mpWb2qpm9aWbvmNk1+QgmkomWFrjyyvSOhd3YCAsWhOJIx/Tp8RxzuwcrOYCZAHybmfQgv2c8KC+Ht97K/PYrVsDgwXDOOeEkyu+/D088EY5hfvLJ4dOQbFxxCtepB4509xoz6wpMN7On3H1mjrOJpG3KlMxGgzU1Yc/AY45J/TYLFqS/nHRtznL6s/5E+lCeo4FulNJAAyWcwzieY+h61/kn/VnBFjnJZAaffpr57X/4Q5g7N4y611VbG948r7kGrrsuWsZOz91TvgDlwGvAARu73uDBgz0Tzc3uixdndNOC0tzsPmGC+z33hO+TqqnJ/fvfd+/a1X3wYPfPP487UcdOPNE9fOBO/1JS4r5iRerLqqzMfFmpXsYy2h28jlJfQc9/Xda90ro/r6PUHXwso3OWqUcP94ceyuz5WbDAvbR04/ffs6f76tWZ3X+cXnrJferUzG8PVHmKXZzSHLeZFZnZG8BS4Fl3z8nqiZEjoV8/uO++XNx7/tx9dzgu83nnJfuxPPpoOKRnY2P4aHzDDXEn6liUrTxKStJbSZmPQ7iO4TeM4XocYzNW/uuyrnV/7hhjuJ4x/CZnmVpaMn/s06ZBcQef893h3Xczu/+4TJ8O3/kODBuW/pRbJlIqbndvdve9gW2A/c1sjw2vY2ajzKzKzKqWLVuWUZjGxvAxLOlzXOvmj2vlVTasu7LOPRl7COZz5VY+luV0YSxjOJwXWMxWrKLtzVhW0Y1P6MthTGMsY0hxTJZ5rgiPvaPbmiVvJWVTU8i95vucS3VovuYC/Adw6cauk+lUSWOj+/z5Gd20oDQ1ud9+u/sdd4Tvk6qhwX3YMPcuXdx339196dK4E3Xs+OOjTZUsX576svIxVbLeFAVf+jy2b/OX89jee/BlfnJEmCqZN6/jqZIePdxXrcrs/uP0zDPuTz6Z+e1JY6qkw5WTZtYbaHT3FWZWBhwF5ORDc3ExbL99Lu45v4qKwjRJ0nXtCo8/HkbaXRKy4eioUeHjeCa7oR94IGy+eerXLytLfxlRNFBCXxa3+but+JR6uuUlR5cumW+/vsMOMGRIeI7aOulEWVl47cS1fXwURx+dv2Wl8nLsC0w1s7eA2YQ57sdzG0sKSVJKG+C734VuGfRXjx4wenR6t9l22/SXE8XRPEsDJQDUUk49JdRSDkAjXTmaPEyuEsbFffpkfvuHHgoDtIqK9X/evTsccgj8539Gy7cp6PAl6e5vufs+7r6nu+/h7tfmI5hIJoqK4Oqrw7bGqSouDivF09kUEELJ5HPUPYL76clKailnHP/GN/ic8YykljJ6spIR3J+XHHV1sNdemd/+G9+AN9+EP/wB9tsPBgyAI46AP/0prAzv2jVrUTst8xysBaisrPSqqqqs369IKtzDx+377+94B5ni4lAkr70Wyjsds2eHPS5X5mH/l640sJwtaKKY03iYZ1j7LnMMU/gzP6CIZrZgOU3ktvm22y4/27Bvasys2t0rU7lugj4Ei6TGDO64A664IsyVtjX6LioKP9933zD6S7e0AQYNyt+ek0U0cw9nsTNz1yttgKc5lp2Zy738mGJyv0nDAQfkfBHSARW3dEpmcNVVsHgxXH99mFPt1i0U9mabwZlnwssvhwMmZTpfW1oK/ftnN3d7VlPGBdzOUtoOu5Q+XMDtrCa3czdlZWGKSOKVyi7vIom1+eZw0UXhkgvf/jbMn5+b+y5EXbtCZUof5iWXNOIWieB73wtbpGwqWlpg773jTiEqbpEITj557R5znV3XruGwFEncxrqzUXGLRNCtG/z7v4fjnHR2RUXw85/HnUJAxS0S2QUXJGsnpUwNHgw77hh3CgEVt0hk/fvDwQfHnSK3evSAyy+PO4WsoeIWyYLRo7++C3dnUloaDicghUHFLZIFRx8dphI6OtZ0ElVUhB2aioriTiJrqLhFssAMJk7M7ABXhaykJBxHZPjwuJPIulTcIlmyzTZwyy2da8qkrCyc0UkKi4pbJItGjuw8Uybl5TB+PPTuHXcS2ZCKWySLOsuUSUkJDB2qKZJCpeIWybJttoH//d/8nyEnW4qLwyaOST7RdWen4hbJgeOOC3PDSSvvoqJwtMTp09M7jZvkl4pbJEdOPx1uuy055V1cHEr7lVeinZpMck/FLZJDI0fCvfcWfnmXlITpkerq/B1jXDKn4hbJsVNPhcmTwwkcCvFgVBUV4dyPs2fDVlvFnUZSoeIWyYOjj4YPP4Rjj03vRMa5VFwcSvvGG+Gll2DLLeNOJKlScYvkSa9eYeR9//3xj74rKmD//eHdd8NhaTeVY4p3Fipu2WS4hzO4xO2UU9aOvsvK8nsMkPLytaPs6dM1n51UKm7p9GbOhGHDwgi3uBgGDAhbezQ2xpdpzej71VfhrLNCgedqV3kz6N4dtt0Wxo6FRYviH2WvWAG//S3ssEOYotlzT5gwARoa4suUKO6+0QuwLTAVeA94B7ioo9sMHjzYN1UffeQ+YoR7z57um23mftZZ7h9/HHeqTdcDD7iXlbmbuYcxd7iUl7sfeqh7fX3cCYMvv3S/7Tb37bZz797dvWvX9fNmcqmocC8tdf/ud92nTnVvaYn7UQYLF7pvvXV4XjbMu//+7nV1cSeMB1DlHXTrmouF67fPzPoCfd39NTPrAVQDJ7n7u+3dprKy0quqqrL01pIcn3wSRg4rVkBzc/hZcXHYkeHtt7XGPt8WLw5nbFm1qu3fl5XBL38JV16Z31wb4w4zZsCLL8K0aTBrFixfnt597LRTOLHDwQeHY2j365ebrJk66KCwBUtT09d/V1oK55wDt96a/1xxM7Nqd69M5bodTpW4+2J3f631+68II++to0Vs3+rVubrn3Bs7FlauXFvaEP44V66E3/0uvlxRvPdeODXXc8/FnSR9//M/oQjbs2oV3Hzz+s9X3MxCsV1xBTz1VNgDM50pjYqKsN34hAmhAAuttN9/H954o+3ShvD6v/tuqKvLb65saGnJ399SWnPcZjYA2AeYlYswN90URkEzZuTi3nNv8uS2500bGmDSpPznyYYTToA//CHMESftxTR1ascDgbo6WLIkP3ky8cEHG3/z2VCXLrBgQc7iRDZrVscrY4uLQ8Enzd57wy675GdZKR980sy6A38FfuHuK9v4/ShgFED/DFdVb7fd2pUVSbSxI8IV4o4XqejVK6zM6tYteYcq7dq14+u0tBT248okWyGfqSaVx9PSktpzV2j23BPq6/OzrJRG3GbWlVDaE939kbau4+53unulu1f2zvAAvqecAvPmwcCBGd08dj/5Sdu7NpeVhV2fk+ipp+CPf4TXX0/em8/3v9/xlhrbbFPYx5veb7/0SqypCfbYI3d5ojryyI635unWDXbdNT95sumBB8JRIfOhw+I2MwPGA++5++9zHym5LroIvvWt9feMKy8Pb0Q/+1l8uaLYfHMYMSJsQpc0P/rRxkd45eVw1VWFvfPJhRemN+oeOBB22y13eaLq2xdOPjmshGxLeXmY3y/kT0GFIJUR98HACOBIM3uj9aLzPbehoiLM4d14Y1ijP2RImLefMaNwdnPelPToAc88Az17rv9JyCw8H+eeCz/+cXz5UrHzzuHvKJUTM5SXw7XX5j5TVOPHh702KyrWvmkWFYX8P/oRXHJJvPmSoMPNATOxqW4OKIVp2TK4667wUXb16nBqsUsugW9/O+5kqfnyy5D1H/9ofw51zUj1qqvymy1T7vDCC3DnnWGzzYED4fzzYa+94k4Wn3Q2B1RxiyRATU0o5nvuCVuO1NeH0eqaY2j/5jdw2mlxp5QoVNwinVRdHTz2WNjSp6QEKivDaLyQ5+klNekUt1YBiCRIeXk4s45s2nSQKRGRhFFxi4gkjIpbRCRhVNwiIgmj4hYRSRgVt4hIwqi4RUQSRsUtIpIwKm4RkYRRcYuIJIyKW0QkYVTcIiIJo+IWEUkYFbeISMKouEVEEkbFLSKSMCpuEZGEUXGLiCSMiltEJGFU3CIiCaPilo1avhyqq+GTT+JOIiJrdFjcZna3mS01szn5CCSFYdUqOOss6NcPjjwSdtwxfFWBi8QvlRH3PcCxOc4hBebUU+Hhh2H1ali5Mnx96SU46CBoaIg7ncimrcPidvcXgS/ykIWmJpg7Nx9Lyr0ZM2DWrLhTZOaDD+D550NZr6upCT7/HB55JJ5cUU2dCrfeCvX1cSeJprkZ7rknuX9f66qpgeeeC39bSeceLvlQUHPcv/41DBwIU6bEnSSaRYvg0ENhyBBYsiTuNOl7+WXo0s5fRk0N/P3v+c2TDbW1cMwxcNllMG5c3GmiefBBOO88OOQQaGyMO000558PRx0Fd9wRd5LozjgDhg/Pz7KKs3VHZjYKGAXQv3//jO7jiCPghRdgl12ylSoeW2wR3oC6dIHNNos7Tfp69my/uIuKYMst85snG0pLwzz9vHkwaFDcaaLZfffwda+9oDhrr+B4DB0aPgkNHhx3kuiGD8/fG6l5CmN7MxsAPO7ue6Ryp5WVlV5VVRUtmcSmthb69AlfN1RWBq++Cnuk9JdQWFpawjRJWVncSaJraQGzcJHOwcyq3b0ylesW1FSJFIaKChg/HsrL1x95V1TAz3+ezNKG8Fg6Q2lDeCwq7U1XKpsDPgjMAAaa2UIz+7fcx5K4/eAHYa779NNht93guOPCSsmxY+NOJiIdzpC5+xn5CCKFZ++9YeLEuFOIyIY0VSIikjAqbhGRhFFxi4gkjIpbRCRhVNwiIgmj4hYRSRgVt4hIwqi4RUQSRsUtIpIwKm4RkYRRcYuIJIyKW0QkYVTcIiIJo+IWEUkYFbeISMKouEVEEkbFLSKSMCpuEZGEUXGLiCSMiltEJGFU3CIiCaPiFhFJGBW3iEjCqLhFRBImpeI2s2PN7AMz+9DMrsh1KBERaV+HxW1mRcAfgOOA3YAzzGy3XAcTySZ3ePxxGDcOPvss7jQi0aQy4t4f+NDd57t7A/AQcGJuYyXf8uWwYkXcKaJraoLJk5P/WC6/HE4/HS66CPbcE1aujDtRNI2N4c2oM2hujjtBdrS0hEs+pFLcWwMfr/Pvha0/y7qlS+GVV3Jxz/l3+OEwdGjcKaJ74AE46SS4IuETZBMmQG0t1NWFrzNmxJ0omuHD4dZb404R3auvhtfK8uVxJ4luzBi49NL8LKs4hetYGz/72nu9mY0CRgH0798/ozC9e0NDQ0Y3LTjXXQddOsGq3yOPhO99D0aMiDtJNLvsAjNnhk8QTU2w445xJ4rmootg4MC4U0Q3aBBccglsvnncSaIbMSJ/nx7MO/i8ZWYHAle7+zGt/x4D4O6/ae82lZWVXlVVlc2cIpEsWQI//Sl8/DFcdRWcqMk+KTBmVu3ulalcN5UR92xgJzPbHlgEnA78MEI+kbzr0wcmTYo7hUh2dFjc7t5kZhcATwNFwN3u/k7Ok4mISJtSGXHj7k8CT+Y4i4iIpKATrD4TEdm0qLhFRBJGxS0ikjAqbhGRhFFxi4gkTIc74GR0p2bLgI8yvHkvoLMcBkiPpfB0lscBneexdJbHAdEey3bu3juVK+akuKMws6pU9x4qdHoshaezPA7oPI+lszwOyN9j0VSJiEjCqLhFRBKmEIv7zrgDZJEeS+HpLI8DOs9j6SyPA/L0WApujltERDauEEfcIiKyEQVV3J3lpMRmdreZLTWzOXFnicLMtjWzqWb2npm9Y2YXxZ0pU2ZWamavmtmbrY/lmrgzRWFmRWb2upk9HneWKMxsgZm9bWZvmFmiD+JvZpub2V/M7P3W18yBOVtWoUyVtJ6UeC5wNOH0aLOBM9z93ViDZcDMDgVqgPvcfY+482TKzPoCfd39NTPrAVQDJyX0OTGgwt1rzKwrMB24yN1nxhwtI2Z2MVAJ9HT34+POkykzWwBUunvit+M2s3uBl9x9nJmVAOXunpOztRbSiLvTnJTY3V8Evog7R1TuvtjdX2v9/ivgPXJ0vtFc86Cm9Z9dWy+FMWpJk5ltAwwDxsWdRQIz6wkcCowHcPeGXJU2FFZx5+2kxJI+MxsA7APMijdJ5lqnF94AlgLPuntSH8vNwGggT+cUzykHnjGz6tbz1ibVDsAyYELrFNY4M6vI1cIKqbhTOimx5J+ZdQf+CvzC3VfGnSdT7t7s7nsD2wD7m1niprHM7HhgqbtXx50lSw52932B44CftU4zJlExsC9wh7vvA9QCOVtPV0jFvRDYdp1/bwN8ElMWadU6H/xXYKK7PxJ3nmxo/Qj7AnBszFEycTBwQuvc8EPAkWb2QLyRMufun7R+XQpMIkyZJtFCYOE6n+L+QijynCik4v7XSYlbJ/ZPBx6LOdMmrXWF3njgPXf/fdx5ojCz3ma2eev3ZcBRwPvxpkqfu49x923cfQDhNfK8u/8o5lgZMbOK1pXetE4rfAepYLQzAAAAo0lEQVRI5JZY7v4p8LGZDWz90VAgZyvxUzrnZD50ppMSm9mDwOFALzNbCPyHu4+PN1VGDgZGAG+3zg0D/LL1HKRJ0xe4t3XrpS7Aw+6e6E3pOoE+wKQwPqAY+JO7T4k3UiQXAhNbB57zgbNztaCC2RxQRERSU0hTJSIikgIVt4hIwqi4RUQSRsUtIpIwKm4RkYRRcYuIJIyKW0QkYVTcIiIJ8/9wNTMgqGJCAwAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7f0acf397190>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "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. Though, if your sense function is implemented correctly\n",
    "# and this output is not converging as expected.. it may have to do\n",
    "# with the `move` function bug; your next task!\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": 28,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAD8CAYAAACMwORRAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4xLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvAOZPmwAAIABJREFUeJzt3Xl4VdW9//H3l5A5RIogIoM4gThcsUTEquiVap2u2p9Di16sVau2trWtV6v41OlpLb21V2u1erk4WxWuVi8OrbZOoEUkDA6AWkqt4FARkCkhIeH7+2Md2xgCZ8jeOefsfF7Pcx6SnJ11vpuT89l7r7322ubuiIhIsvTIdwEiIhI9hbuISAIp3EVEEkjhLiKSQAp3EZEEUriLiCSQwl1EJIEU7iIiCaRwFxFJoJ75euG+ffv60KFD8/XyIiJFae7cuR+7e790y+Ut3IcOHUp9fX2+Xl7kM1pbYfFimDsXFi2CDRugpgb22QdGjYLhw6GHjnOlAJjZ3zJZLm/hLlIIVq6EW2+FX/4SNm4MP1u//p/P19SAO9TWwg9+AOeeC71756dWkWxoX0S6JXeYOhV23RV+8hP4+OMQ6m2DHcL3GzbABx/AlVeG5R9/PD81i2RD4S7dTksLnHkmnHMOrF37zz32dBobYfVq+MpX4Nvfhs2b461TpDPULSPdijuccUbY+25oyK2Nhga4886wkbj1VjCLtkaRKGjPXbqVX/+6c8H+qYYGuPdeePDBaOoSiZrCXbqNd96BSy/tfLB/qqEBzj8f/v73aNoTiZLCXbqNq6+GpqZo22xqgp//PNo2RaKgcJduYe3aMDqmtTXadpubYfLk6DcaIp2VNtzNrMLMXjGzV81soZld08Ey5WY21cyWmNlsMxsaR7Eiufr976G0NJ62zeCFF+JpWyRXmey5NwFHuPt+wEjgaDMb026Zc4DV7r47cAPws2jL7FhS7u2dlCF17oW7LrNmbTmGPSoNDaCLraXQpA13Dz79WJSmHu1j9UTg7tTXDwHjzOIbILZsGQwZAuXlcNddcb1K1/jWt6BnTxgzJoyjLlavvBKu3KypgT/+Md/VbGn27Ph2Blpa4E9/iqftzkrKDpBkL6M+dzMrMbMFwEfAH9x9drtFBgLLANy9BVgDbB9loW3dd1+4YnDTJrjiirheJX5r14b+WndYuLC4D+1/9rOwPo2NcNVV+a5mS3HttX9q7dp428/VFVfAbbfluwppa9o0mD8//tfJKNzdvdXdRwKDgNFmtk+7RTraS99in8HMzjOzejOrX7FiRfbVphxwAJSVQXU1fOELOTeTd9XVMHhw+BdgxIj81tMZY8dCVVVYl7Fj813NluLqb/9UWVm87efq1FPh2GPzXYW0td9+YSK6uGV1haq7f2JmzwNHA2+0eWo5MBhYbmY9ge2AVR38/mRgMkBdXV3OB4xf/GLYy122DI4/PtdW8q+kJPTVPv001NXBzjvnu6Lcffe7sNdeYdRIIYbJnnvCvHnxtG0WZo8sRPvvn+8KpL2uCHbIbLRMPzPrnfq6Evgi8Ga7xaYDX0t9fQrwrHu8vX11dfDlL8e/Rxa37beH8eNhjz3yXUnnmMGRR4aNbSFOjXvooeHIIg41NeGciUghyWTPfQBwt5mVEDYG09z9cTO7Fqh39+nA7cC9ZraEsMf+1dgqFsnBEUfEd3Jx06bC7IqS7i1tuLv7a8AWB3fufmWbrzcCp0Zbmkh0hg0LXSdz5kTf9mGHwU47Rd+uSGcU4AG0SDyuuSb6rpnKyjDPu0ihUbhLt3HMMXDUUdGNbKmogNNPL+4RW5JcCnfpVqZMCSexO3vSt2dPGDAAbrwxmrpEoqZwl25l++3hpZdghx1yH2lVVgYDB4Z2amqirU8kKgp36XZ22QUWLAgnQj+9gCxTVVVhHP+8eWHPXaRQKdylW+rfP1w8NnkyDB0a9sC31lXTo0fYCAwbBvffD488An36dGm5IlnTPVSl2zILJ0THjw9dLNOnw4wZ8Je/hHnay8rC1YRjx8JJJ8Ho0fmuWCRzCnfp9szgkEPCQyQp1C0jIpJACncRkQRSuIuIJJDCXUQkgRTuIiIJpHAXEUkghbuISAIp3EVEEkjhLiKSQAp3EZEEUriLiCSQwl1EJIEU7iIiCaRwFxFJIIW7iEgCpQ13MxtsZs+Z2WIzW2hmF3WwzOFmtsbMFqQeV8ZTroiIZCKTm3W0ABe7+zwz6wXMNbM/uPuidsvNdPfjoy9RRESylXbP3d0/cPd5qa/XAYuBgXEXJiIiucuqz93MhgL7A7M7ePogM3vVzH5nZntv5ffPM7N6M6tfsWJF1sWKiEhmMg53M6sBHga+5+5r2z09D9jZ3fcDfgU82lEb7j7Z3evcva5fv3651iwiImlkFO5mVkoI9t+4+2/bP+/ua919ferrJ4FSM+sbaaUiIpKxTEbLGHA7sNjd/2sry+yYWg4zG51qd2WUhYqISOYyGS1zMDABeN3MFqR+NhEYAuDutwGnAN80sxagEfiqu3sM9YqISAbShru7vwhYmmVuBm6OqigREekcXaEqIpJACncRkQRSuIuIJJDCXUQkgRTuIiIJpHAXEUkghbuISAIp3EVEEkjhLiKSQAp3EZEEUriLiCSQwl1EJIEU7iIiCaRwFxFJIIV7F3OHWbPg+9+HUaOgthYqKqC6GoYPh7PPhkcfhZaWfFfafWzeDL/7HZx/PuyzD9TUQHk59OoF++0HF14Izz4b3juRYmH5uqdGXV2d19fX5+W182X6dPjBD+DDD6GxMYRKR3r1gtJSuOyysHxJSdfW2V24w+TJcNVVsGEDrF/f8XJmYePbuzf89KdwxhnhZyL5YGZz3b0u7XIK9/h98gl8/evw9NPQ0JD571VXwy67wMMPw7Bh8dXXHS1bBqeeCm+8EYI9U9XVcOCBcP/90L9/fPWJbE2m4a5umZitWAF1dfDkk9kFO4TQWbgQDjgA5s6Np77u6K23YORIqK/PLtghLD9zZvj9v/0tnvpEoqBwj9HGjTB2LLz7LjQ359aGO6xdC0ccAUuXRltfd/T3v8Mhh8Dq1dDamlsbmzaFjfbBB8OaNdHWJxIVhXuMJk4Mwb5pU+fb2rABvvKVrffTS3rucNZZIZA72xvZ2gorV8K3vhVJaSKRU7jH5LXX4Lbbsu+K2ZrWVli8GKZMiaa97uj//i90qUSxsYVwZPboo/D889G0JxIlhXtMJk2CpqZo29ywAX78Yw3Jy9U112Tfx55OQwNce220bYpEQeEegzVr4JFH4ulCWb067H1KdhYtCidS4zBrVhh9I1JI0oa7mQ02s+fMbLGZLTSzizpYxszsJjNbYmavmdnn4ym3OMyeDWVl8bTd2Ah//GM8bSfZc8/F13ZpKcyYEV/7IrnIZM+9BbjY3UcAY4ALzWyvdsscA+yRepwH3BpplR144QW4554QdoWmvj66vvb2WlvDuhca93C08uCDuY9CidPMmfH9raxfH/beRQpJ2nB39w/cfV7q63XAYmBgu8VOBO7x4GWgt5kNiLzalCeegGOPhW9+M1yIUmgWL453+oB33omv7Vz9/OcwYQKccw5cckm+q9nS22/H17Y7vPlmfO13xvPPh79HKRyLFsG6dfG/TlZ97mY2FNgfmN3uqYFA217H5Wy5AcDMzjOzejOrX7FiRXaVtvHaa2HceEMDzJ+fczOxiWo0xtYU4rwzc+aEk5UNDfDKK/muZktx/5/leh1D3B56KMybI4Vj8WL461/jf52emS5oZjXAw8D33H1t+6c7+JUtxnS4+2RgMoTpB7Ko8zPOOgvuvhuWL4df/CLXVuLTp0+87dfUxNt+Lq64IvQ7NzeH+VcKTW1tvO1/7nPxtp+rm2/OdwXS3sknd83rZBTuZlZKCPbfuPtvO1hkOTC4zfeDgPc7X17HBgwo3MNgCNMNVFdHP+yubfuFZuTIcPVnoTrooNAvHscIptLS0L5IIclktIwBtwOL3f2/trLYdODM1KiZMcAad/8gwjqLSl1dfLMGVlWFy+clO2PGhA1uHCorC3ODK91bJnvuBwMTgNfNbEHqZxOBIQDufhvwJHAssARoAL4efanFY999Ybvttj6FbGe0tsKXvxx9u0n3pS/F1+9uFuYQEikkacPd3V+k4z71tss4cGFURRU7szBiZOLEaIdEmsFRR8GOO0bXZndRUwP//u9w553Rhnx5eRi1Fdd1DSK50hWqMTn33LD3HqWKCvjZz6Jtszu56qrwfxiligr4j/+Itk2RKCjcY1JdHS7oqayMpr2qKvjRj2DEiGja644GDoQbb4yu772qCu66C7bfPpr2RKKkcI/R2LFhsqqqqs61U1UF48bBpZdGU1d3dvbZcNpp0bwnF1wAJ50UTV0iUVO4x+ySS8KsgbnuwVdXw3HHhVvt6V6qnWcWpk2eMCH3gK+qgu98B66/PtraRKKkcO8CF18cJq4aMiTzLoHKynCj7MmTYerUMJZaotGjR5hrf+rUcPFRpiFfXQ077ACPPRamdNZNsqWQKdy7yIEHhvlNbrsN9tknjLKorYWeqfFKZiHMq6rCjZevvjpconz66QqRuBx/fJinZ9IkGDw4bFBra0P4QzhSqq0NJ0133TVcDb10abjloUihM8/TnR/q6uq8vr4+L69dCD74INz0euHCMB6+vBz22ANGjYLddlOgdzX3cMPruXPD1c8bN4YN7YgR4T0ZNEjviRQGM5vr7mkvm1O4i4gUkUzDXd0yIiIJpHAXEUkghbuISAIp3EVEEkjhLiKSQAp3EZEEUriLiCSQwl1EJIEU7iIiCaRwFxFJIIW7iEgCKdxFRBJI4S4ikkAKdxGRBEob7mZ2h5l9ZGZvbOX5w81sjZktSD2ujL5MERHJRs8MlrkLuBm4ZxvLzHT34yOpSEREOi3tnru7zwBWdUEtIiISkaj63A8ys1fN7HdmtndEbYqISI4y6ZZJZx6ws7uvN7NjgUeBPTpa0MzOA84DGDJkSAQvLSIiHen0nru7r3X39amvnwRKzazvVpad7O517l7Xr1+/zr60iIhsRafD3cx2NAv3hTez0ak2V3a2XRERyV3abhkzewA4HOhrZsuBq4BSAHe/DTgF+KaZtQCNwFfd3WOrWERE0kob7u4+Ps3zNxOGSoqISIHQFaoiIgmkcBcRSSCFu4hIAincRUQSSOEuIpJACncRkQRSuIuIJJDCXUQkgRTuIiIJpHAXEUkghbuISAIp3EVEEkjhLiKSQAp3EZEEiuI2e12qoQEWLIDly6G1FXr3hv32g512yndl2XGHv/wFFi4M61RaCrvvDnvvHb6WrtfaCm++CW+9BRs3QmUljBgBw4ZBD+0GSZEpinBvaIAHHoDrr4clS6CqKoQjhA9dU1P4IJ55Jnz3u7Drrvmtd1vmzAnr8dhjYAY9e4Z1MQv/btwIo0fDpZfCccdBSUm+K042d3jmGfjP/4QXXoCysvA39el70toKLS1w5JHhPTnkkPBzkUJn+bppUl1dndfX16dd7uGH4Zxzwods/fptL1taGsJwwgS44Qaoro6o2Ai89x6ccUYI940bYfPmbS9fUwN9+8K0aXDAAV1TY3ezaBGceiq8+276vy2z8Pc0fDhMnQq77dY1NYq0Z2Zz3b0u3XIFe7DZ1ASnnAJf+xqsWZP+wwewaVMIzvvuC10cCxfGX2cmnngC9twTXnopHIWkC3YI6/vOO3DYYXDNNf88UpFo/PrXUFcHixdn9rflHpabPx/+5V/g3nvjr1GkMwoy3Jub4Zhj4MknYcOG7H+/sRH+/nf4whfgjTeiry8bjz4a9g7Xrw+H99lqbAxdBpdeGn1t3dX118Mll4T/22w3mps3hw30+efD7bfHU59IFAoy3L//fXj55fDhy5U7rFsH48aFf/Nh0aLQFdOZ9YAQJr/+dTgikc55+mm46qrwf7otxrYPrxob4TvfgVmzIixOJEIFF+4vvgh33tn5QIQQ8GvXhpOsXa2lBU47LZr1gBBG3/oWfPBBNO11R2vWwOmnpw/2PqxkGYP5HKu2uVxjY7TvsUiUCi7cv/GNaD8sGzeGE2CLFkXXZiYeeCD0mUfZV75xI1x+eXTtdTeTJmXWzXcC0xnI+5zA9LTLrloVjqpECk1BhfucObBsWfTtbtoEN94YfbvbkmmQZGPTpjB6Zu3aaNvtDpqb4dZbwwYynbO46zP/bktDA/ziF5mdJBfpSmnD3czuMLOPzKzDU5MW3GRmS8zsNTP7fK7FTJkSzyFuS0vor+6qD+Cf/wx//Ws8bZeUwCOPxNN2kj3zTGZHUb1Yy4G8DMAYXqYX6bek69fD7NmdrVAkWplcxHQXcDNwz1aePwbYI/U4ELg19W/WZsyIL4B79AgXQA0bFk/7bb3ySnwXH61fDzNnhiGikrmXX97ySKo3qxnCu5/52TieoZlyKmimmTLOZQrPMO4zy7zLED7hc//4ftOm8J4fdFBs5YtkLe2eu7vPgG2eWToRuMeDl4HeZjYg20I2bw6X48elRw+YNy++9tuaMyezsdO5KsQRGsuWhaGno0aFS/gLzcyZ4UK4ti5jEq8ykpcZwwzGMoOxXMU11BKGV9Wyjqu45h/PvcwYXmUklzHpM+1s3BgGAogUkij63AcCbXvKl6d+lpWNG+O9UKelBVavjq/9tj78MN72P/kk3vZzcdlloWti/vwwqqfQrFy55c8u56dcznU4xnas/cejrbY/d4zLuY7L+ekWbX30UVyVd86vfhWutZDCcc01oZciblGEe0czbXQY02Z2npnVm1n9ihUr2j0X/1WYXTX5U9xzjxTiJFY77BCmfygtDV8Xmo7eE6cHk7icw3meD9iRRio6/N1GynmfARzGC0zicryDj00hvicAQ4cW36R6Sbf77rDjjvG/ThQThy0HBrf5fhDwfkcLuvtkYDKEuWXaPldRAeXl6ccg56q0tOv+yHfbLfS5t+8GiMqArDu94nfddbDdduEIbOLEfFezpSFD4NVXO35uDqMZzlssYCS7suWZ8A/YiZEsYB21W21/l12iqjRa//Zv+a5A2jvjjK55nSj2N6YDZ6ZGzYwB1rh71pfamMFee0VQzVY0N4f+4K5wwAHxTVpmBmPHxtN2Z1RWwtVXhyGgtVvPwLw57LAw4+PWNFPGADr+s92RD2mifKu/W10NBx/c2QpFopXJUMgHgFnAcDNbbmbnmNkFZnZBapEngaXAEuB/gJx7XL/0pW1/ADujV6+u23MfMyZMfBaHmpoQVJKdQw7Z9t/WkfyBZsICG6iiiTI2UAXAJko5kj9s9Xc3b1a4S+HJZLTMeHcf4O6l7j7I3W9399vc/bbU8+7uF7r7bu6+r7unn8d3K77xjXj6LisrwzwgXaV/fzj00Hja7tkzTKom2Rk9OkyhvDUTuJda1rKBKqZwDtuzkts5mw1UUstaJrD1aSCHDw+zfooUkoI6DbTzzqHLoWcMtxA577zo29yWK66Ivmumqgouuiie/5+kMwtTN3T0npTSzHE8wVpq+X/8lu9xExuo4SJ+xcn8lrXUchxP0JNNW/xudXV4r0UKTUGFO8DkyeHEalSqq+HHPw57013p8MOj72bq2xd++MPo2utuzjknjFRof3RYQit38TWG8TZP86XPPPcURzOMt7mbM+nJZ+ds7tkTPv95OPnkuCsXyV5B3olpypSwh9rZkTNlZTByZLjoJx9D1VatCofsK1d2fphnVRU891zoXpDcvf12COQo5v2prQ33Cxg8OP2yIlEp6jsxnXsuXHxxCLRclZeHe6k+9VT+xiD36RMuVujdu3M1VFbCb36jYI/CsGHw+OOd+9uCcIL+mWcU7FK4CjLcAa69Nsy2V1WV/Twt1dXhJh2zZ4dgzacRI8J0BHvumX2glJeHDcTjj8NJJ8VTX3d0+OHw/PPheoHKyux+t6oqXBj00kvhNn0ihapgwx3gggvg9dfDnCWVleFCpG3p1Qv69YM77gj3LS2U8da77RYuoJk4MWx4amq2vXxFRXicemqYb+eII7qmzu7kgAPCRHJnnx3+r9NteKurw9/ghReGuXP23bdr6hTJVUH2uXfkrbfgllvCbdKWLg0ns8zCjHx9+oQP6/nnh2GCcc3IGIWGhnDzkLvuCoHf2BjW5dPZMPfaC7785TAstKtPAndXq1eHHYJp08JN1Vtawt9QS0s4etp333AHpzPPDDsQIvmUaZ970YR7W5s2hZOUra1h77yYP3CrVoXALysLo2EKdY6S7sIdPv44XIRWWRl2HOKeK0gkG5mGe1GOmC4t7ZqJd7pCnz7hIYXBLHTtiRQ77SeKiCSQwl1EJIEU7iIiCaRwFxFJIIW7iEgCKdxFRBJI4S4ikkAKdxGRBFK4i4gkkMJdRCSBFO4iIgmkcBcRSSCFu4hIAincRUQSKKNwN7OjzewtM1tiZpd18PxZZrbCzBakHudGX6qIiGQq7XzuZlYC3AIcCSwH5pjZdHdf1G7Rqe7+7RhqFBGRLGWy5z4aWOLuS929GXgQODHeskREpDMyCfeBwLI23y9P/ay9k83sNTN7yMwGR1KdiIjkJJNw7+gOku1vvPoYMNTd/wX4I3B3hw2ZnWdm9WZWv2LFiuwqFRGRjGUS7suBtnvig4D32y7g7ivdvSn17f8AozpqyN0nu3udu9f1040qRURik0m4zwH2MLNdzKwM+Cowve0CZjagzbcnAIujK1FERLKVdrSMu7eY2beBp4AS4A53X2hm1wL17j4d+K6ZnQC0AKuAs2KsWURE0jD39t3nXaOurs7r6+vz8toiIsXKzOa6e1265XSFqohIAincRUQSSOEuIpJACncRkQRSuIuIJJDCXUQkgRTuIiIJpHAXEUkghbuISAIp3EVEEkjhLiKSQAp3EZEEUriLiCSQwl1EJIHSzudeSF59FWbOhBkzYOlSaG2FPn3gkEPgwANh3DgoL893lel98gk8/TTMmgVz5sCGDVBWBnvvHdZl3DjYeed8V5mee1iHP/0pvCfvvRd+tuOOcOihcNBB4d+SknxXKtINuXteHqNGjfJMbN7sfv/97iNGuFdVuVdUuIcI+eejpMS9tjY8LrnEfeXKjJrucn/+s/vpp4d16NXL3WzLdamuDs8fdpj7Cy/ku+KONTe7/+pX7oMGudfUuJeVbbkepaVhHXfYwX3SJPeGhnxXLZIMhJskpc3Ygr5Zx7JlMH48LFgQ9m4zUV4OlZVw991wwgkRFBqB1la44Qa48kpobg7fZ6KqCk47DW66CXr1irfGTL3+Opx6Kixfnvl7UlUVjrCmTQt78yKSu6K/Wce8ebDvvjB7duYhAtDUFLo9xo8PYZpvzc1hI3P11dDYmHmwAzQ0wIMPwn77wYcfxlZixp54AsaMgbffzu49aWgIG4Nx4+Cee+KrT0T+qSDDffFi+Nd/hTVroKUltzYaGuAXv4Drrou2tmy4h73c557LLgzb2rgxHMF84Qtho5Uvzz4b1qWhIaxXLhob4YIL4KGHoq1NRLZUcOG+aROcdBKsW9f5thoa4Mc/Dict82HyZHjmmRBqndHSAu+/H4IxH1atglNO6fx6QGjjrLPCnryIxKfgwv2668IHP6pTAY2NYY9z06Zo2svUe+/BxRfnvsfeXlMTPPYY/P730bSXjW9+M2woo9LUBBMmRNeeiGypoMK9sRGuvz7aIIGw5zl9erRtpnPDDdFvUBoaYOLEaNtM5913w/9dU1N0bba0wCuvhKGtIhKPggr3adPiaXfdOpg0KZ62O9LUFLpkmpujb/vNN+GNN6Jvd2tuuQU2b46+3aamcE5EROKRUbib2dFm9paZLTGzyzp4vtzMpqaen21mQ3MpZupUWL8+l99M79VXo+siSWfOHDCLp+3WVnjyyXja7sjDD8ezkWptDaNvRCQeacPdzEqAW4BjgL2A8Wa2V7vFzgFWu/vuwA3Az3IpZt68XH4rM5WVYbx8V5g7N9pujLaam+GFF+Jpu72mptAtE5cNG+Cjj+JrP1cNDTB/fhipJFKsMtlzHw0scfel7t4MPAic2G6ZE4G7U18/BIwzy27fddMmWLEim9/ITmsrvPVWfO23NX9+fOEOYahoV1i6FCoq4mu/oqLr3pNMbdgAe+0VpoHYf/+uPxEftdWroxnlJNFZtiyeo+H2Mgn3gcCyNt8vT/2sw2XcvQVYA2zfviEzO8/M6s2sfkW7JG9uhh4xngHYvLnr9sSiPiHcXpwbjvavE1f30qcKbe948WL4+OPwHr7zTngUs//+79DdKYXjiivgf/83/tfJZOKwjj7e7QcqZrIM7j4ZmAxh+oG2z5WVxXPi7lM9esS7F9pWTU287XfV5GgVFfG+J5++RiHZe2/YaacwHHfYMBg6NN8Vdc6FFxbHZHrdyc9/DttvsesbvUz2lZcDg9t8Pwh4f2vLmFlPYDtgVTaFlJbCDjtk8xvZKSmBESPia7+tkSPjDa292p/xiMkuu8R7lNDY2HXvSaYqK8P8OfPmhRPjpaX5rqhzevUKO05SOPr3h55dMB9vJuE+B9jDzHYxszLgq0D7UePTga+lvj4FeNZzmJFs1KhsfyNzDQ1hjpauUFcX3weqrAwOOyyettsrL493z7VXL+jbN772c1VeDnvuWfzBLt1b2nBP9aF/G3gKWAxMc/eFZnatmX067+LtwPZmtgT4AbDFcMlMjB8fX5dGXV2YnbAr1NXF11ddUgLHHx9P2x057bR4DutLSuDE9qflRSQyBTXlb1MT9OsXzbwybdXUwL33hjlrusoPfwg33hj9WfEDDghXd3aV996D3XeP/sRnZWVYj332ibZdkaQryil/y8vhssugujradvv379q9XYDvfS/6rpmqKvjpT6NtM52BA+Hkk6M9h1BaGoYaKthF4lNQ4Q5w6aXhFnNRdWtUVoYpZrviBEZbAwbAL38Z3YaqoiLMzDhuXDTtZeOWW6Ld4FZUhJupiEh8Ci7ce/aERx+F2trOt1VVBT/5SRi9kg9f/zoce2zn+/pLS2HIELj55mjqytZ228Ejj0RzzqKyEu6/P2z8RCQ+BRfuAHvsEW6E/bnP5T5ioaoKfvQj+P73o60tG2YhyDoT8JWVsOuu8NJL+b0TdEkzAAAEsUlEQVTV3qGHho1udXXuF5tVVYU99q7uIhPpjgoy3CHcYm/RIjj88Oy6BCoqwvC6Rx4J/ff51rNnmO3y+uvDemTaPWQWgv3ss8OY60IYMnjkkWHenH32ye49qaqC3XaDF18Mc+uLSPwKNtwBdtwRnnoK7rsvjIGvrOw4VMrKQjdOnz5hlMpf/gJHHdX19W6NWbjhxcKFIayrqkK9JSVbLterVzixfOyx8PzzoSumq4ZwZmL48LCxuemmENjV1eF9aa+iIoxSGjQonARetCjM1SIiXaOghkKm8/bbobvmpZdCgLe2hkA/+GAYPTp0HXT1idNcrF8f7kn6yivhsWFD2EDtuy8cdFA4WimGPmn3MNPmrFnhfXnvvfCz/v1h7Fg48MDwvsQ9P41Id5LpUMiiCncRke6u4MPdzFYAf+tkM32BjyMopxAkZV20HoUnKeuSlPWAzq3Lzu7eL91CeQv3KJhZfSZbsGKQlHXRehSepKxLUtYDumZdCvqEqoiI5EbhLiKSQMUe7pPzXUCEkrIuWo/Ck5R1Scp6QBesS1H3uYuISMeKfc9dREQ6ULThbmZHm9lbZrbEzApgooHcmNkdZvaRmb2R71o6w8wGm9lzZrbYzBaa2UX5rikXZlZhZq+Y2aup9bgm3zV1hpmVmNl8M3s837V0hpm9Y2avm9kCMyvaC2TMrLeZPWRmb6Y+KwfF9lrF2C1jZiXA28CRhPu3zgHGu/uivBaWAzMbC6wH7nH3op3h3MwGAAPcfZ6Z9QLmAicV23tiZgZUu/t6MysFXgQucveX81xaTszsB0AdUOvuRTtlm5m9A9S5e1GPczezu4GZ7j4lddvSKnf/JI7XKtY999HAEndf6u7NwINAUd60zd1nkOXNxAuRu3/g7vNSX68j3JJxYH6ryp4H61PflqYexbcHBJjZIOA4YEq+axEws1pgLOG2pLh7c1zBDsUb7gOBZW2+X04RBklSmdlQYH9gdn4ryU2qK2MB8BHwB3cvyvUAbgQuBTbnu5AIOPC0mc01s/PyXUyOdgVWAHemusqmmFnE9537p2IN946moirKvaukMbMa4GHge+6+Nt/15MLdW919JDAIGG1mRdddZmbHAx+5+9x81xKRg93988AxwIWp7sxi0xP4PHCru+8PbABiO19YrOG+HBjc5vtBwPt5qkVSUn3UDwO/cfff5ruezkodMj8PHJ3nUnJxMHBCqq/6QeAIM7svvyXlzt3fT/37EfAIoWu22CwHlrc5EnyIEPaxKNZwnwPsYWa7pE5KfBWYnueaurXUicjbgcXu/l/5ridXZtbPzHqnvq4Evgi8md+qsuful7v7IHcfSvh8POvu/57nsnJiZtWpk/SkujGOAopudJm7fwgsM7PhqR+NA2IbcFAEs59vyd1bzOzbwFNACXCHuy/Mc1k5MbMHgMOBvma2HLjK3W/Pb1U5ORiYALye6q8GmOjuT+axplwMAO5OjcjqAUxz96IeRpgA/YFHwv4DPYH73f33+S0pZ98BfpPaKV0KfD2uFyrKoZAiIrJtxdotIyIi26BwFxFJIIW7iEgCKdxFRBJI4S4ikkAKdxGRBFK4i4gkkMJdRCSB/j+LHjrKFys1bwAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7f0acf505b90>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "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": 29,
   "metadata": {},
   "outputs": [
    {
     "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-29-019c317be5a5>\u001b[0m in \u001b[0;36m<module>\u001b[0;34m()\u001b[0m\n\u001b[0;32m----> 1\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[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     46\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     47\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---> 48\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     49\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": [
    "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": 30,
   "metadata": {},
   "outputs": [
    {
     "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-30-7f72fc6f19d8>\u001b[0m in \u001b[0;36m<module>\u001b[0;34m()\u001b[0m\n\u001b[1;32m     29\u001b[0m         \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[1;32m     30\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m---> 31\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;32m<ipython-input-30-7f72fc6f19d8>\u001b[0m in \u001b[0;36mtest_robot_works_in_rectangle_world\u001b[0;34m()\u001b[0m\n\u001b[1;32m     27\u001b[0m     \u001b[0;32mfor\u001b[0m \u001b[0mi\u001b[0m \u001b[0;32min\u001b[0m \u001b[0mrange\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;36m1000\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m     28\u001b[0m         \u001b[0msimulation\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0msim\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mSimulation\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mgrid\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mblur\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mp_hit\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m---> 29\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     30\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m     31\u001b[0m \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[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     46\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     47\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---> 48\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     49\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": [
    "# 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": 31,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "> /home/workspace/localizer.py(48)move()\n",
      "-> new_G[int(new_i)][int(new_j)] = cell\n",
      "(Pdb) c\n",
      "> /home/workspace/localizer.py(47)move()\n",
      "-> pdb.set_trace()\n",
      "(Pdb) c\n",
      "> /home/workspace/localizer.py(48)move()\n",
      "-> new_G[int(new_i)][int(new_j)] = cell\n",
      "(Pdb) c\n",
      "> /home/workspace/localizer.py(47)move()\n",
      "-> pdb.set_trace()\n",
      "(Pdb) c\n",
      "> /home/workspace/localizer.py(48)move()\n",
      "-> new_G[int(new_i)][int(new_j)] = cell\n",
      "(Pdb) c\n",
      "> /home/workspace/localizer.py(47)move()\n",
      "-> pdb.set_trace()\n",
      "(Pdb) c\n",
      "> /home/workspace/localizer.py(48)move()\n",
      "-> new_G[int(new_i)][int(new_j)] = cell\n",
      "(Pdb) c\n",
      "> /home/workspace/localizer.py(47)move()\n",
      "-> pdb.set_trace()\n",
      "(Pdb) c\n",
      "> /home/workspace/localizer.py(48)move()\n",
      "-> new_G[int(new_i)][int(new_j)] = cell\n",
      "(Pdb) c\n",
      "> /home/workspace/localizer.py(47)move()\n",
      "-> pdb.set_trace()\n",
      "(Pdb) c\n",
      "> /home/workspace/localizer.py(48)move()\n",
      "-> new_G[int(new_i)][int(new_j)] = cell\n",
      "(Pdb) c\n",
      "> /home/workspace/localizer.py(47)move()\n",
      "-> pdb.set_trace()\n",
      "(Pdb) c\n",
      "> /home/workspace/localizer.py(48)move()\n",
      "-> new_G[int(new_i)][int(new_j)] = cell\n",
      "(Pdb) c\n",
      "> /home/workspace/localizer.py(47)move()\n",
      "-> pdb.set_trace()\n",
      "(Pdb) c\n",
      "> /home/workspace/localizer.py(48)move()\n",
      "-> new_G[int(new_i)][int(new_j)] = cell\n",
      "(Pdb) c\n",
      "> /home/workspace/localizer.py(47)move()\n",
      "-> pdb.set_trace()\n",
      "(Pdb) c\n",
      "> /home/workspace/localizer.py(48)move()\n",
      "-> new_G[int(new_i)][int(new_j)] = cell\n",
      "(Pdb) c\n",
      "> /home/workspace/localizer.py(47)move()\n",
      "-> pdb.set_trace()\n",
      "(Pdb) c\n",
      "> /home/workspace/localizer.py(48)move()\n",
      "-> new_G[int(new_i)][int(new_j)] = cell\n",
      "(Pdb) c\n",
      "> /home/workspace/localizer.py(47)move()\n",
      "-> pdb.set_trace()\n",
      "(Pdb) c\n",
      "> /home/workspace/localizer.py(48)move()\n",
      "-> new_G[int(new_i)][int(new_j)] = cell\n",
      "(Pdb) c\n",
      "> /home/workspace/localizer.py(47)move()\n",
      "-> pdb.set_trace()\n",
      "(Pdb) c\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-31-9d9590a250fe>\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;32m<ipython-input-30-7f72fc6f19d8>\u001b[0m in \u001b[0;36mtest_robot_works_in_rectangle_world\u001b[0;34m()\u001b[0m\n\u001b[1;32m     27\u001b[0m     \u001b[0;32mfor\u001b[0m \u001b[0mi\u001b[0m \u001b[0;32min\u001b[0m \u001b[0mrange\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;36m1000\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m     28\u001b[0m         \u001b[0msimulation\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0msim\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mSimulation\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mgrid\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mblur\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mp_hit\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m---> 29\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     30\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m     31\u001b[0m \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[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     45\u001b[0m             \u001b[0mnew_i\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0;34m(\u001b[0m\u001b[0mi\u001b[0m \u001b[0;34m+\u001b[0m \u001b[0mdy\u001b[0m \u001b[0;34m)\u001b[0m \u001b[0;34m%\u001b[0m \u001b[0mwidth\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m     46\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[0;32m---> 47\u001b[0;31m             \u001b[0mpdb\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mset_trace\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m     48\u001b[0m             \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[1;32m     49\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": [
    "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": null,
   "metadata": {},
   "outputs": [],
   "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.14"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}