{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Interactions\n",
    "\n",
    "This is a collection of interactions, mostly from the book. If you have are reading a print version of the book, or are reading it online via Github or nbviewer you will be unable to run the interactions.\n",
    "\n",
    "So I have created this notebook. Here is how you run an interaction if you do not have IPython installed on your computer.\n",
    "\n",
    "1. Go to try.juptyer.org in your browser. It will launch a temporary notebook server for you.\n",
    "\n",
    "2. Click the **New** button and select `Python 3`. This will create a new notebook that will run Python 3 for you in your browser.\n",
    "\n",
    "3. Copy the entire contents of a cell from this notebook and paste it into a 'code' cell in the notebook on your browser. \n",
    "\n",
    "4. Press CTRL+ENTER to execute the cell.\n",
    "\n",
    "5. Have fun! Change code. Play. Experiment. Hack.\n",
    "\n",
    "Your server and notebook is not permanently saved. Once you close the session your data is lost. Yes, it says it is saving your file if you press save, and you can see it in the directory. But that is just happening in a Docker container that will be deleted as soon as you close the window. Copy and paste any changes you want to keep to an external file.\n",
    "\n",
    "Of course if you have IPython installed you can download this notebook and run it on your own computer. Type\n",
    "\n",
    "    ipython notebook\n",
    "    \n",
    "in a command prompt from the directory where you downloaded this file. Click on the name of this file to open it."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Experimenting with FPF'\n",
    "\n",
    "\n",
    "The Kalman filter uses the equation $P^- = FPF^\\mathsf{T}$ to compute the prior of the covariance matrix during the prediction step, where P is the covariance matrix and F is the system transistion function. For a Newtonian system $x = \\dot{x}\\Delta t + x_0$ F might look like\n",
    "\n",
    "$$F = \\begin{bmatrix}1 & \\Delta t\\\\0 & 1\\end{bmatrix}$$\n",
    "\n",
    "$FPF^\\mathsf{T}$ alters P by taking the correlation between the position ($x$) and velocity ($\\dot{x}$). This interactive plot lets you see the effect of different designs of F has on this value. For example,\n",
    "\n",
    "* what if $x$ is not correlated to $\\dot{x}$? (set F01 to 0)\n",
    "\n",
    "* what if $x = 2\\dot{x}\\Delta t + x_0$? (set F01 to 2)\n",
    "\n",
    "*  what if $x = \\dot{x}\\Delta t + 2*x_0$? (set F00 to 2)\n",
    "\n",
    "* what if $x = \\dot{x}\\Delta t$?  (set F00 to 0)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXAAAAD7CAYAAABzGc+QAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XlYlWX+P/D3h0UUUFxAAVdMTMUFc5lxP65jWspojTrz\ns2zRKdtNrWz66djXmb7NNT8ra7zUaiwrbUqz3EXzKGqpJZLghtsILqmoICD7/fvj5rCDBzjwnOfw\nfl3XfZ3t4TwfQd7c537u535EKQUiIjIfN6MLICKiqmGAExGZFAOciMikGOBERCbFACciMikGOBGR\nSXnU9A5EhPMUiYiqQCklFb1eKz1wpZTTt/nz5xteg6vUaYYaWSfrdPZmDw6hEBGZFAOciMikGOD5\nLBaL0SXYxQx1mqFGgHU6GuusfWLvWEuVdyCianofRESuRkSgnOEgJhEROR4DnIjIpBjgREQmxQAn\nIjIpBjgRkUkxwImITIoBTkRkUgxwIiKTYoATEZkUA5yIyKQY4EREJsUAJyIyKQY4EZFJMcCJiEyK\nAU5EZFIMcCIik2KAExGZVLUCXERai8guEYkTkVgRed5RhRERUcWqdUk1EQkEEKiUOiIivgB+BhCh\nlDpeZBteUo2IqJJq/JJqSqkrSqkj+fdTARwHEFyd9yQiIvs4bAxcRNoB6AnggKPek4iIyueQAM8f\nPvkawAv5PXEiIqphHtV9AxHxBLAWwGdKqfVlbbNgwYKC+xaLBRaLpbq7JSJyKVarFVartVJfU92D\nmALgEwBJSqmXytmGBzGJiCrJnoOY1Q3wgQD2APgFgO2NXlNKbS2yDQOciKiSajzA7SyCAU5EVEk1\nPo2QiIiMwwAnIjIpBjgRkUkxwImITIoBTkRkUgxwIiKTYoATEZkUA5yIyKQY4EREJsUAJyIyKQY4\nEZFJMcCJiEyKAU5EZFIMcCIik2KAExGZFAOciMikGOBERCbFACciMikGOBGRSTHAiYhMigFORGRS\nDHAiIpNigBMRmRQDnIjIpBjgREQmxQAnIjIpBjgRkUkxwImITIoBTkRkUgxwIiKTYoATEZkUA5yI\nyKQY4EREJlXtABeRj0XkVxE56oiCiIjIPo7ogf8bwGgHvA8REVVCtQNcKRUF4KYDaiEiokrwMLoA\nIkdTSuHWrVu4fPlysZaUlIT09HSkp6cjIyMDOTk5yMnJgZubG7y9veHj41Os+fv7o3Xr1gWtQYMG\nRv/TiIqplQBfsGBBwX2LxQKLxVIbuyUXl5eXh3PnziEuLg5xcXGIjY1FXFwcTp48iYyMDIfvr1mz\nZmjTpg26dOmC8PBwhIeHo0ePHggICHD4vqjusVqtsFqtlfoaUUpVe8ci0g7ABqVUtzJeU47YB1Fm\nZiZ++ukn7NmzB3v27MG+fftw+/btMrdt2LAhgoKCEBgYiKCgIAQFBSEgIAC+vr7w9vaGl5cXPD09\n4eHhgdzcXKSlpSE9PR1paWlIS0tDamoqrl69igsXLiAhIQEXL15EdnZ2mfsKDg5G7969MXToUIwY\nMQJhYWEQkZr8VlAdICJQSlX4H4kBTk7tzJkzWLt2LbZu3YoffvihVM86KCgIXbt2RVhYWMFt586d\n4efn59A68vLy8Ouvv+L8+fM4evQojhw5gpiYGMTExCAtLa3Yti1atMCwYcMwYsQIPPjgg+yhU5XU\nSoCLyGoAQwA0A3AVwP9VSv27yOsMcLKbUgrHjh3D2rVrsW7dOsTExBR7PSwsDIMHDy5owcHBBlWq\n5eXl4cyZM9i/fz927NiBnTt34vLlywWvu7u7Y8SIEZgyZQoiIiIc/oeFXFet9cDvUgQDnO4qOTkZ\nq1atwrJlyxAbG1vwfKNGjfDAAw8gIiICFovF6XuzSimcOHECO3bswJYtWxAZGYmcnBwAgJeXF8aM\nGYMnnngC999/P9zceB4dlY8BTk7v559/xtKlS7F69Wqkp6cD0AcLx48fj4kTJ2L48OHw8vIyuMqq\nu379OtatW4fVq1dj9+7dsP0uhIaG4oUXXsCjjz4KX19fg6skZ8QAJ6eklMLWrVuxcOFC/PjjjwXP\nDx06FE8//TQiIiLg6elpYIU149KlS/jss8/wwQcf4MKFCwCAxo0bY/r06XjppZcQFBRkcIXkTBjg\n5FSUUoiMjMT8+fMLgrtx48aYNm0a/vznP6NTp04GV1g7cnJysH79erzzzjvYt28fAMDb2xuzZs3C\nnDlz0KhRI4MrJGdgT4BDKVWjTe+C6ro9e/aoAQMGKAAKgPL391dvv/22Sk1NNbo0Qx08eFBFREQU\nfF8CAgLUe++9pzIzM40ujQyWn50V5it74FSjrl27hjlz5uCTTz4BADRt2hRz5szBs88+y7HfIvbt\n24e5c+di//79AIAOHTpg2bJlGDZsmMGVkVE4hEKGycvLw8cff4y5c+fi5s2b8PLywiuvvILZs2ej\nYcOGRpfnlJRS+Pbbb/Hqq6/i5MmTAIAnnngC//znPzn9sA5igJMhzpw5g0cffbRgfHfkyJH44IMP\nEBoaanBl5pCVlYW3334bb775JrKystCmTRusXLkSQ4cONbo0qkX2BDgnopJDffbZZwgPD8e+ffsQ\nGBiI1atXY9u2bQzvSqhXrx7+8pe/IDo6Gr1798aFCxcwbNgw/PWvf0VeXp7R5ZETYQ+cHCIjIwPP\nPfccPvzwQwDAww8/jGXLlqFJkyYGV2Zu2dnZWLRoERYuXAilFCZOnIiVK1fy+EEdwCEUqhWJiYkY\nP348Dh8+jPr162PJkiV44oknuKCTA23ZsgVTpkxBcnIyunfvjvXr1yMkJMTosqgGMcCpxp04cQKj\nRo1CQkIC2rdvj7Vr1yI8PNzoslzSyZMnMW7cOJw6dQr+/v6IjIzk99qFcQycatTBgwcxcOBAJCQk\noF+/fjh06BADpQbde++9OHDgAEaNGoXr169j+PDhpRb7orqFAU5Vsn37dgwbNgxJSUkYM2YMduzY\ngaZNmxpdlstr3LgxvvvuO4wdOxY3btxgiNdxDHCqtL1792LcuHFIS0vDI488gvXr18Pb29vosuoM\nLy8vrF27FmPHjkVSUhKGDx+OuLg4o8siA3AMnCrl2LFjGDBgAG7duoUZM2Zg6dKlXBbVIJmZmZgw\nYQI2b96M9u3b49ChQ/wU5EJ4EJMcKjExEf3790dCQgIiIiLw9ddfw93d3eiy6rQ7d+5gwIABiI6O\nxujRo7Fx40b+TFwED2KSw6Snp2PMmDFISEhA//798cUXXzAonECDBg3wzTffwN/fH1u3bsUbb7xh\ndElUi9gDN5vsbCAtDUhPr7jZs01GBpCXp5tShfdLPlYKly9dQvLNm/Dy9ESbVq3gLlLutsXex80N\n8PTUrV69sm8res2ebWz3vb0BH5/yW716gIvOTd+1axdGjhyJ3NxcbNmyBaNHjza6JKomDqE4q6ws\n4OpV4Pp1ICmpeCv53I0bQGpqYejmX56LqsDdvXSo+/oCDRsCjRrpW1sr+ris15o00X8QnMhbb72F\n1157De3atUNcXBwPLJscA7y2ZWUB587pduWKbr/+Wvr2xo2q78MWQt7ehT1O2/2KWlnbeXnp9xPR\nPWU3t+L33dxw7fp1THjoISTduoU5c+fisccfL/Z6RV8LWy89O1u3rKzS98u7rew2mZmFnzzKa47+\n4+fjAzRrVro1bVr2882aAX5++ntTA7Kzs9GnTx/ExMTg1Vdfxd///vca2Q/VDgZ4TcjNBRISgPh4\n4NQp3Wz3z5/Xr9+NuzsQEKBbWb/k/v7Fw6BRo8LQreVLjU2YMAHffPMNRo0ahS1btph7xklWVulQ\nT00Fbt8GUlL0bdFW8jnb45QU4Natqv1BcHPTvXd7gz8gAGjeHPDwsOvtDxw4gH79+sHd3R2HDx9G\nt27dKl8jOQUGeFUppXvKRcPZdv/0ad3bK4sI0K4dcM89QFAQEBgItGihm+1+YKD+xTRBEEZFRWHw\n4MHw8fHByZMn0bJlS6NLch5K6SC3DXOVHAor2WzbpKRUfl9ubvr/TnBw8dayZfHH+f+vZs6ciaVL\nl+L+++/H5s2bHf9vp1rBALeXUsCZM4DVWtguXix/+6AgoGNH3UJDC++3b6+HJVyAUgr9+vXDgQMH\nsGDBAsyfP9/oklxDdnb5gV/W81evAteu2ffenp5AUBCymzfHpsOHkZCXh4eefx5BvXoVD3o/P5c9\nmOtKGODlUQo4e7Z4YCcmFt/Gzw/o1Kl0SHfooA9kubj//Oc/mDRpEgIDAxEfH8/lS42UlaWPn1y6\nVHG7edO+92vQoHigt2mjPzW2b69v27Sxe8iGag4DvKSLF4FVq4CVK4H8S1YVaNYMGDIEGDoUsFiA\nLl1MMcxRE5RS6Nq1K44dO4Zly5ZhxowZRpdE9rhzB7h8Gbh0CVePHMFbzz+PYABPjR8P3+TkwqC/\nfbvi93F3B9q21WFeNNhtt3WgA+MMGOCA7m1v3gwsWQJERupZEYA+kGSxFAZ2WFidDeyS9u7di0GD\nBiEwMBAXLlyAZy0fOCXHmDJlCtasWVN6COz27YKgx8WL+uD7mTP6U+mZM/q5in5nAwJKh7rtflAQ\nf48chAF+8iTw4ovA1q36cb16wLhxwLRpwO9+x4+J5XjkkUewatUqzJs3D4sWLTK6HKqizZs3Y+zY\nsejZsycOHz5s/xdmZJQOddv9s2f16+WpXx8ICdGBHhqqO0a2xp57pdTtAP/oI+Dpp/VBIz8/4C9/\nAR5/XE/VonLduHEDwcHByMzMxNmzZ3nVFxPLyMhAQEAAUlNTcf78ebRt27b6b5qXp3vvJYPddlvR\nAde2bXWQd+2qW1gY0LmzHpOnUuwJcNfsgu7cCfz5z3pO9hNPAH/7m55LS3e1bt06ZGZmYuTIkQxv\nk6tfvz7uv/9+fPXVV/j222/x/PPPV/9N3dz09MWWLYFBg0q/npJS2FM/eRKIiwNiY4Hjx4H//le3\nolMb3dz0xIDw8MLWs6eebkt35XoBnpMD/OEPOrxfeQV46y2jKzKV77//HgAwfvx4gyshR3jwwQfx\n1VdfITIy0jEBfjeNGhUGcVE5ObqXHhurmy3YbedYnDoF/Oc/hdu3aFEY5rbb0FBOfyzB9YZQTp/W\nP+iAAD31igdU7KaUQnBwMK5cuYK4uDh06dLF6JKomo4dO4awsDCEhITg7NmzRpdTWmam7p0fOaJb\ndLS+LeuEpxYtgMGD9WyxIUNcfqZY3RwDT0jQ81g9PfWZk44Y96sjTpw4gc6dO6N58+a4cuUKryrv\nArKzs+Hj44Ps7GykpqbCx8fH6JLuTil9ELVooB88qM+OLqpZMz2MYwv07t31FEgXUTfHwFu1AiIi\ngPXrgUcfBbZs4UESO+3btw8AMGTIEIa3i/D09ETHjh0RFxeH48ePo3fv3kaXdHcieiZLSAjw+9/r\n55TSwyy7dwN79ujbxET9e75+vd7Gz694oPfs6fIzzar9+UNERovICRGJF5FXHFFUNQsCli3TBy13\n7wb69tVjbXRX58+fBwB07tzZ2ELIoWw/z1OnThlcSTWIAPfeC8yYAXz2GXDhgh5T//hj3VFr1w5I\nTgY2bgTmzNG/902aAKNHA//4h97eBVUrwEXEHcD7AEYD6AJgiogY/9vfvDmwbZs+9T02FujTB/if\n/6naQkJ1SEJCAgCgdevWBldCjtSsWTMAQHJyssGVOJCIPnHoscf0mdXnzukZLp9+qmeedeigV5rc\ntg2YO1cPpQ4erDt3SUlGV+8w1e2B9wVwWil1XimVDWANAOeYvhAeDvz8s577nZEBvPGG/iEuWGD/\nmhF1DAPcNdnWsUlNTTW4khrWpg0wdSrw4Yf6+NfFi8AXXwCTJukTjKKigKee0meLjhsHfPmlXkPe\nxKob4C0BJBR5nJj/nHPw9dUn9OzYof/63roF/PWvOshnzNBDLLZT6wkX81dgbNWqlcGVkCPZrl2a\na89a9a4kOBiYMgVYs0av6vjJJ8CoUXqK8YYNwOTJembLk0/q8XUTqu4Iv13TSxYsWFBw32KxwGKx\nVHO3lTR8uG579gBvvqkDfcUK3Vq10j/kP/1JH8WuwwfvsrKyAOgTQMh1OM1ickZq2BB45BHdfv1V\n974//1zPbvnoI+Df/9Y99Xnz9FmiBrBarbBarZX7IqVUlRuA3wLYWuTxawBeKbGNcjqxsUq99ppS\nbdsqpY9v69ali1KLFil19qzRFRqiffv2CoCKj483uhRyoKlTpyoA6qOPPjK6FOdz8qRS06cr5elZ\nmANTpyp165bRlan87Kwwg6s7hPITgFARaSci9QBMAvBdNd+z5oWF6dPrz54F9u7Va6Y0awYcOwa8\n/ro+ONKxIzB9ul5+1kWPYJdk63nfuXPH4ErIkXhsowIdOwLLl+sTAJ95Rl+QZdUqfQzthx+Mru6u\nqhXgSqkcAM8C2AbgGIAvlVLHHVFYrXBzAwYMAP71L71Az8aNejjF11cfBPnwQ/2Rq21bPSd12jT9\nUevs2YqX2zSphvmrxaVwto5LScy/WAmPbVSgTRvg/feBo0eBXr30iUSDBhXOMXdSrncmpiPk5ACH\nDxeeNBAVpeeYFtWqVeEJA0OGuMQ6DQ899BDWrl2LTz/9FFOnTjW6HHKAnJwc+Pr6IjMzE7dv3+aV\nleyRlaXnkr/3nv5kHhtryOJa9pyJ6boLCVSHh4c+EWDOHH20OilJB/rixfosz6ZN9Vlgn3+uZ7Pc\ne68+4j1pkv4rvn+/noNqMl3zD97ExcUZXAk5SnR0NDIzM3HPPfcwvO1Vrx7wzjvAyJH6d3/ePKMr\nKpdrn2fqKO7u+rTcnj31BSLy8vRqart3F7YrV/RqarYV1UR0r7zoEpnh4U69TGZYWBgABrgr2bVr\nFwBg6NChBldiMiLASy/pq3idPm10NeVigFeFmxvQrZtuzz6rx8NPnNBB/uOPevGduLiKl8ns1q1w\nUfsuXQBvb+P+PflsAR7LpQdcBgO8ipTSUw0B/enaSXEMvKbYlsm0raZma2UdILSdFmwLdNsVS+69\nV3+cqyU5OTkICAjArVu3cOrUKYSGhtbavsnx7ty5g+bNmyM1NRWXLl1CUFCQ0SWZQ3a2Pl/kzTf1\nrJTdu4Hf/KbWy6ibqxE6Cy+v0gvbF10ms+ii9idPFl6e6ttvC7d3c9NHx8u7Orifn0NL9vDwwNix\nY/H555/j22+/xezZsx36/lS71q5di9TUVPTu3Zvhba/ISOCFF3TnC9DrrBgQ3vZiD9wZZGXpoRZb\noNvC/fTpiqcrNmtW/tXBW7as0mL3X3/9NR5++GEMGDAAe/furcY/iow2ePBgREVFYfny5Zg+fbrR\n5TivpCR9uv3KlcBPP+nnOnQA3n0XGDPGsLLq5gUdXElmpu6xF71obNELyVZ0wk29eoVXBy8Z7u3b\nl7tG+u3btxEQEICsrCwkJiYi2InH/6h8x48fR5cuXeDr64tLly4VzPEnFH4StlqBTZv0TLP8ZSTg\n5we8+qo+gOnlZWSVHEIxPS8vPQ5+772lX1NKz3wp6+rgZ87oxXtOntStLP7++uBMidYwOBgv9O+P\nVbt24YN338Wi//3fmv03Uo1YsmQJAOCPf/wjwxvQy81arYWt6NnVbm563fBHHwXGjzfVBWDYA3dV\nqamFVwcvGfDnz+uTle4iF4AEBsKtZcsyw76g+fu79LUJzebUqVMICwtDXl4efvnll4LZRXVCVpYO\n6/h4PSwZE6MPQv73v8W3a9pUn4BnsQATJ+ohRyfDIRQqW24ucO0acOlS2e3iRdyIi0PT7Gz73s/T\nU6+xXDTUg4L0GH1Zjasd1qiJEydi3bp1ePLJJ7FixQqjy3G8vDx97VtbSNtafLwO77KWzW3SpDCw\nLRY9jdfJOx0McKqyTZs2IeKBB9CjeXPsXr0aPsnJ5Qf+jRuVe3Nvb90DKi/gy3qtcWOXumBtTdm3\nbx8GDhwIb29vxMfHm/MYRm4ucP26Xvb1ypXiYR0frw/uZ2SU/bUieu2i0FC9UFWnTsDAgXqpaCcP\n7JI4Bk5VNmbMGNzXty8OHjyIv2zYgMWLF5e/cUaGXgysZLAnJel240bh/aQkfRWU9HS9HIG9RHQv\nqqxw9/PT6z03aqRvba3k4/r1Tb9eTUUyMjIwc+ZMAMDLL7/sPOGdm1v6/4CtXb1aGNS222vX7n6h\nlcDAwpDu2LHw/j331KlPeOyBU7mio6MLrmJ+6NAh3HfffdV/U6X0+Hx54V7e87duVX/f7u7FQ93H\np3jz9S39XMnWoIGe4ePpWfq25HPu7rX6B+Oll17CO++8gw4dOiA6Orrya5/k5uqZTbY/sEVbWlrZ\nz5f32s2b1fvZNWumz1oODNRDcqGhhSEdGqp/ji6OQyhUbbNmzcLixYvRq1cvHDhwoODyXLUuJ6cw\nFEqGe3IycPt2YUtJKfu+bapYbbI37EveeuR/OM7L002pwvtlPL5x4wbijh6FO4DuXbvC18fn7l+b\nk1M8sMsblqgu26ensobLAgJ0SNvCukULfVFyT8+aqcVEGOBUbampqejSpQsSEhIwb948LFq0yOiS\nqi4zs3jQp6VVvt25o0+1zs7WfxAqujXjNSi9vQubj0/xx/Y+36CBHtZq1kzPUOLxiyphgJND7Nq1\nCyNHjkRubi6+/vprTJw40eiSzCEvz/6wL3mbna17rm5uupVzPz0jAy/NmoUT8fHo3qMH3nnvPbh7\neJT9dSUfu7sXD2AXP0ZgNgxwcpjFixdj1qxZ8PHxwYEDB+rW3GInlZ2djQceeADbt2/HPffcg/37\n96N58+ZGl0UOwgs6kMO8+OKL+NOf/oS0tDSMHz8eV65cMbqkOk0phSeffBLbt29HQEAAtm7dyvCu\ng9gDJ7ulp6dj0KBBOHz4MDp16oRdu3Yh0IkvUOGqcnNzMXPmTCxfvhze3t6wWq3o06eP0WWRg7EH\nTg7l7e2Nbdu2oVu3bjhx4gSGDh3Knngty8jIwB/+8AcsX74c9evXx7p16xjedRh74FRp165dw/Dh\nw3H06FF07twZkZGRaOmEa0m4muTkZERERMBqtcLPzw8bN27EwIEDjS6Lagh74FQjAgICsHPnTnTr\n1g3Hjx9Hnz59cODAAaPLcmnnzp3DkCFDYLVaERQUhKioKIY3McCpagICAvD999/DYrHg8uXLGDx4\nMD755BOjy3JJGzduxH333YeYmBiEhoZi//796Natm9FlkRNggFOV+fv7Y/v27Zg5cyaysrIwbdo0\nzJo1C1lGnPHogjIzMzFr1iw8+OCDuHXrFsaNG4eDBw+iXbt2RpdGToJj4OQQy5cvxzPPPIOcnBz0\n7NkTq1at4lzxavjpp58wffp0HDlyBB4eHli0aBFmz54NN5OtqEdVxzFwqjUzZszA7t27ERISgujo\naPTq1Qt/+9vfkG3vmuIEQB+ofO6559C3b18cOXIEISEh2Lt3L+bOncvwplL4P4Icpn///oiJicGT\nTz6JzMxMvP766+jTpw92795tdGlOTymFNWvWoFOnTnj//ffh5uaG2bNn45dffsFvnPiq6GQsDqFQ\njdixYwdmzJiBc+fOAQDGjh2Lt956C127djW4MueilEJkZCTmz5+PH3/8EYD+Q7h06VJ0797d4OrI\nSBxCIcOMGDECR48excKFC+Hr64tNmzahR48eePzxx5FYmQs5uCilFHbu3IlBgwbhd7/7HX788Uf4\n+/tjxYoViIqKYniTXdgDpxp39epVLFy4EMuWLUNOTg48PT0xadIkvPjii+jVq5fR5dWq7OxsfPfd\nd3j33XcRFRUFAGjatCnmzp2LZ555pvIXYSCXxdUIyanEx8fjjTfewFdffYW8/EtmDRo0CC+++CLG\njx9v3MUiakFiYiJWrFiBFStW4PLlywCAJk2aYPbs2XjuuefQsGFDgyskZ8MAJ6d0/vx5vP/++1ix\nYgVSUlIAAEFBQZg0aRImT56Mvn37QlxgXeqUlBRs3rwZa9aswcaNG5Gbf4GHTp064amnnsJjjz2G\nRnXg0mBUNQxwcmq3b9/GypUrsWTJEsTHxxc8HxISgsmTJ+Phhx9Gjx49TDV9LikpCRs2bMDatWux\nffv2gpOaPDw8MGHCBDz99NMYMmSIS/yBoppVowEuIg8DWACgE4A+SqnD5WzHAKcKKaVw6NAhrF69\nGl9++WXBEAOgz/YcNmwYhg8fjhEjRqB9+/YGVlpaSkoK9u3bh927d2PPnj04ePBgQU9bRDBw4EBM\nmDABkydP5tK7VCk1HeCdAOQBWAbgZQY4OUJubi6ioqKwevVqbNmyBQkJCcVeb9u2LXr37o3w8HCE\nh4ejR48eaNWqVa30aJOTkxEXF4fY2FjExsZi//79iI6OLhjPBwB3d3cMGzYMEyZMQEREBEObqqxW\nhlBEZBcY4FQDlFI4ffo0duzYgZ07d+L777/HzZs3S23XtGlTdO7cGa1bty7VGjduDB8fH/j4+MDb\n27vMA6W5ubm4c+cOUlJScPny5VLt3LlziI2NxcWLF0t9rYeHB/r06YMhQ4Zg8ODB6N+/P/z8/Grk\n+0F1CwOcXEpubi7i4uIQExODI0eOFNwmJSXZ/R5eXl5o0KABlFLIyclBdna23Ytv1a9fH126dEFY\nWBjCwsLQu3dv9OvXD97e3lX9JxGVy54A97jLG0QCKOsz4Dyl1AZ7C1mwYEHBfYvFAovFYu+XEhVw\nd3dH9+7d0b17d0ydOhWA7qVfvHgRZ86cQUJCQrGWmJiI5ORkpKWlIS0tDenp6cjMzERmZmax9xUR\neHt7w8fHB4GBgQgKCkJQUBCCg4MRFBSE1q1bIywsDCEhIS491ZGMZbVaYbVaK/U17IFTnaGUwp07\nd3Dnzh24u7vDw8MDHh4e8PLy4qwQcjrV7oFXZl8Oeh+iGmPraXPIg1xFlSfYisjvRSQBwG8BbBKR\nLY4ri4iI7oYn8hAROSGuRkhE5MIY4EREJsUAJyIyKQY4EZFJMcCJiEyKAU5EZFIMcCIik2KAExGZ\nFAOciMikGOBERCbFACciMikGOBGRSTHAiYhMigFORGRSDHAiIpNigBMRmRQDnIjIpBjgREQmxQAn\nIjIpBjh45aAiAAAFpElEQVQRkUkxwImITIoBTkRkUgxwIiKTYoATEZkUA5yIyKQY4EREJsUAJyIy\nKQY4EZFJMcCJiEyKAU5EZFIMcCIik2KAExGZVJUDXET+ISLHRSRGRNaJiJ8jCyMioopVpwe+HUCY\nUqoHgFMAXnNMSUREZI8qB7hSKlIplZf/8ACAVo4piYiI7OGoMfDHAWx20HsREZEdPCp6UUQiAQSW\n8dI8pdSG/G1eB5CllPqiBuojIqJyVBjgSqmRFb0uItMAjAEwvKLtFixYUHDfYrHAYrHYWx8RUZ1g\ntVphtVor9TWilKrSzkRkNIB/AhiilLpewXaqqvsgIqqrRARKKalwm2oEeDyAegBu5D/1g1JqZhnb\nMcCJiCqpRgO8EkUwwImIKsmeAOeZmEREJsUAJyIyKQY4EZFJMcCJiEyKAU5EZFIMcCIik2KAExGZ\nFAOciMikGOBERCbFACciMikGOBGRSTHAiYhMigFORGRSDHAiIpNigBMRmRQDnIjIpBjgREQmxQAn\nIjIpBjgRkUkxwImITIoBns9qtRpdgl3MUKcZagRYp6OxztrHAM9nlh+qGeo0Q40A63Q01ln7GOBE\nRCbFACciMilRStXsDkRqdgdERC5KKSUVvV7jAU5ERDWDQyhERCbFACciMqlaC3AReVlE8kSkaW3t\nszJE5E0RiRGRIyKyU0RaG11TWUTkHyJyPL/WdSLiZ3RNZRGRh0UkTkRyReQ+o+spSURGi8gJEYkX\nkVeMrqcsIvKxiPwqIkeNrqUiItJaRHbl/7xjReR5o2sqSUTqi8iB/N/vYyLyd6NrqoiIuItItIhs\nqGi7Wgnw/DAcCeC/tbG/KnpbKdVDKRUOYD2A+UYXVI7tAMKUUj0AnALwmsH1lOcogN8D2GN0ISWJ\niDuA9wGMBtAFwBQR6WxsVWX6N3SNzi4bwEtKqTAAvwXwjLN9P5VSGQCG5v9+dwcwVEQGGlxWRV4A\ncAxAhQcpa6sH/v8AzK2lfVWJUup2kYe+AK4bVUtFlFKRSqm8/IcHALQysp7yKKVOKKVOGV1HOfoC\nOK2UOq+UygawBsB4g2sqRSkVBeCm0XXcjVLqilLqSP79VADHAQQbW1VpSqn0/Lv1ALgDuGFgOeUS\nkVYAxgD4EECFs1BqPMBFZDyARKXULzW9r+oSkUUicgHAowDeMroeOzwOYLPRRZhQSwAJRR4n5j9H\n1SQi7QD0hO5cOBURcRORIwB+BbBLKXXM6JrKsRjAHAB5d9vQwxF7E5FIAIFlvPQ69Ef8UUU3d8Q+\nq6KCOucppTYopV4H8LqIvAr9TXysVgvMd7c687d5HUCWUuqLWi2uCHvqdFKcO1sDRMQXwNcAXsjv\niTuV/E+u4fnHjbaJiEUpZTW4rGJE5AEAV5VS0SJiudv2DglwpdTIcorpCiAEQIyIAPrj/s8i0lcp\nddUR+66M8uoswxcwsGd7tzpFZBr0R6zhtVJQOSrx/XQ2FwEUPUjdGroXTlUkIp4A1gL4TCm13uh6\nKqKUShaRTQB6A7AaXE5J/QGME5ExAOoDaCQinyqlHilr4xodQlFKxSqlWiilQpRSIdC/JPcZEd53\nIyKhRR6OBxBtVC0VEZHR0B+vxucfmDEDwz51leMnAKEi0k5E6gGYBOA7g2syLdG9s48AHFNKvWN0\nPWUREX8RaZx/vwH0pAqn+x1XSs1TSrXOz8vJAL4vL7yB2p8H7swfXf8uIkfzx8gsAF42uJ7yLIE+\nyBqZP83oX0YXVBYR+b2IJEDPStgkIluMrslGKZUD4FkA26CP9H+plDpubFWlichqAPsBdBSRBBEx\nZEjPDgMA/B/omR3R+c3ZZs8EAfg+//f7AIANSqmdBtdkjwozk6fSExGZFM/EJCIyKQY4EZFJMcCJ\niEyKAU5EZFIMcCIik2KAExGZFAOciMikGOBERCb1/wEbhM1o5h+khwAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7e73630>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "%matplotlib inline\n",
    "from IPython.html.widgets import interact, interactive, fixed\n",
    "import IPython.html.widgets as widgets\n",
    "import numpy as np\n",
    "import numpy.linalg as linalg\n",
    "import math\n",
    "import matplotlib.pyplot as plt\n",
    "from matplotlib.patches import Ellipse\n",
    "\n",
    "def plot_covariance_ellipse(x, P, edgecolor='k', ls='solid'):\n",
    "    U,s,v = linalg.svd(P)\n",
    "    angle = math.atan2(U[1,0],U[0,0])\n",
    "    width  = math.sqrt(s[0]) * 2\n",
    "    height = math.sqrt(s[1]) * 2\n",
    "\n",
    "    ax = plt.gca()\n",
    "    e = Ellipse(xy=(0, 0), width=width, height=height, angle=angle,\n",
    "                edgecolor=edgecolor, facecolor='none',\n",
    "                lw=2, ls=ls)\n",
    "    ax.add_patch(e)\n",
    "    ax.set_aspect('equal')\n",
    "    \n",
    "    \n",
    "def plot_FPFT(F00, F01, F10, F11, covar):\n",
    "    \n",
    "    dt = 1.\n",
    "    x = np.array((0, 0.))\n",
    "    P = np.array(((1, covar), (covar, 2)))\n",
    "    F = np.array(((F00, F01), (F10, F11)))\n",
    "\n",
    "    plot_covariance_ellipse(x, P)\n",
    "    plot_covariance_ellipse(x, np.dot(F, P).dot(F.T), edgecolor='r')\n",
    "    #plt.axis('equal')\n",
    "    plt.xlim(-4, 4)\n",
    "    plt.ylim(-4, 4)\n",
    "    plt.title(str(F))\n",
    "    plt.xlabel('position')\n",
    "    plt.ylabel('velocity')\n",
    "                 \n",
    "interact(plot_FPFT, \n",
    "         F00=widgets.IntSliderWidget(value=1, min=0, max=2.), \n",
    "         F01=widgets.FloatSliderWidget(value=1, min=0., max=2., description='F01(dt)'),\n",
    "         F10=widgets.FloatSliderWidget(value=0, min=0., max=2.),\n",
    "         F11=widgets.FloatSliderWidget(value=1, min=0., max=2.),\n",
    "         covar=widgets.FloatSliderWidget(value=0, min=0, max=1.));"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.4.3"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 0
}