diff --git a/Applied-Plotting-Charting-And-Data-Representation-in-Python/week3/week3_Assignment.ipynb b/Applied-Plotting-Charting-And-Data-Representation-in-Python/week3/week3_Assignment.ipynb
index 34205a0..8442664 100644
--- a/Applied-Plotting-Charting-And-Data-Representation-in-Python/week3/week3_Assignment.ipynb
+++ b/Applied-Plotting-Charting-And-Data-Representation-in-Python/week3/week3_Assignment.ipynb
@@ -49,16 +49,827 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 57,
-   "metadata": {
-    "collapsed": false
-   },
+   "execution_count": 151,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "# Use the following data for this assignment:\n",
+    "\n",
+    "import pandas as pd\n",
+    "import numpy as np\n",
+    "import matplotlib.pyplot as plt\n",
+    "import matplotlib\n",
+    "%matplotlib inline\n",
+    "%matplotlib notebook\n",
+    "\n",
+    "np.random.seed(12345)\n",
+    "\n",
+    "df = pd.DataFrame([np.random.normal(32000,200000,3650), \n",
+    "                   np.random.normal(43000,100000,3650), \n",
+    "                   np.random.normal(43500,140000,3650), \n",
+    "                   np.random.normal(48000,70000,3650)], \n",
+    "                  index=[1992,1993,1994,1995])\n",
+    "\n",
+    "# confidence_interval_1 = stats.norm.interval(0.95, loc=32000, scale=2000)\n",
+    "# confidence_interval_2 = stats.norm.interval(0.95, loc=43000, scale=1000)\n",
+    "# confidence_interval_3 = stats.norm.interval(0.95, loc=43500, scale=1400)\n",
+    "# confidence_interval_4 = stats.norm.interval(0.95, loc=48000, scale=700)\n",
+    "\n",
+    "#Get year average (year_avg), year standard devition (year_std), number of year (year_count) and standard error (yerr):\n",
+    "year_avg = df.mean(axis = 1)\n",
+    "year_std = df.std(axis = 1)\n",
+    "year_count = df.count(axis = 1)\n",
+    "yerr = year_std*1.96/np.sqrt(year_count)\n",
+    "\n",
+    "#Get the lower and upper bound error value in the plot coordinates:\n",
+    "lower_bound = year_avg - yerr\n",
+    "upper_bound = year_avg + yerr"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "**Easiest option implementation:** Implement the bar coloring as described above - a color scale with only three colors, (e.g. blue, white, and red). Assume the user provides the y axis value of interest as a parameter or variable:"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 253,
+   "metadata": {},
    "outputs": [
     {
      "data": {
-      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYcAAAEICAYAAAC0+DhzAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAFsxJREFUeJzt3X20XXV95/H3hwCCDxAeUpYSamhFOsG1REyBmdrKyAwE\nnBRmLWVwHEkdNEzFWa0zjkXrkhbFaqdKq1VLWrJ4sEqpjoItGiPCdOqIcKkKRQUiVQPyEAmCimiB\n7/xxftc55HeTe5Pc3HNv7vu11lln7+/+7b1/e+fkfs5+OOekqpAkadhuo+6AJGn2MRwkSR3DQZLU\nMRwkSR3DQZLUMRwkSR3DQfNGklcm+eyo+zEVSZYkqSS7j7ovmp8MB82oJN9K8uMkPxx6PGsm1l1V\nf1lVJ8zEupIcl+SuGVrXc5NcmWRjkk1J1iY5fLM2b0hyb5KHk6xJ8pShaW9PckuSx5L83lbWs6YF\n1nN24uZoljAcNAorqurpQ4/vbsvMGdhtsto2LnMuv0NfCFwFHA4cBNwAXDk+McmJwDnA8cCzgV8A\nfn9o/vXAm4C/3dIKkrwI+MXp7rhmL8NBs0aSY5P83yTfT/LVJMcNTbsuyflJvgA8AvzCFmr7Jrko\nyT1J7k7yjiQL2jJ+I8nfDy2zkpyd5A7gjlb7kyQb2jvsm5L86lb6+5Qkf5TkO0nuS/JnSfZO8jTg\n08CztnZ0lOSlSb7c1rVha+/at6aqbqiqi6pqU1X9M3ABcHiSA1qTlcBFVXVrVT0IvB34jaH5L6mq\nTwM/2MJ27g68H/iv29M/zU2Gg2aFJAczeOf6DmB/4I3Ax5MsGmr2KmAV8Azg21uoXQw8BjwHeAFw\nAvCaraz6VOAYYGkbvxE4svXhI8BfJ9lrC/O+C3hua/8c4GDgbVX1I+Ak4LuTHB39CDiDwTv/lwK/\nmeTUrfR1qn4NuLeqHmjjRwBfHZr+VeCgofCYzBuAv6uqm6ehb5ojDAeNwifb0cH3k3yy1f4TcHVV\nXV1VT1TVOmAMOHlovovbu9/H2jvkJ9UY/EE/GfjtqvpRVd3P4F306Vvpyx+0d9w/BqiqD1fVA20d\n7wGewuB0zZMkCYNQekOb/wfAOydZ15NU1XVVdUvb3puBjwIvnur8E0myGPgA8N+Gyk8HHhoaHx9+\nxhSWdwhwFvC2HemX5p65fJ5Vc9epVfW5zWrPBl6eZMVQbQ/g2qHxDRMsa7j27DbPPYO/3cDgDdBE\n8024zCRvBM4EngUUsA9w4ATzLQKeCtw0tK4AC7ayridJcgyDo4/nAXsyCKK/nur8EyxvEfBZ4INV\n9dGhST9ksB3jxocnPI20mT8GzquqhyZtqV2KRw6aLTYAl1XVwqHH06rqXUNtJvoK4eHaBuAnwIFD\ny9inqo7Yynp/Nn+7vvAm4DRgv6payOBddiaY73vAj4Ejhta1b1U9fSt93dxHGFxIPqSq9gX+bAvr\nmlSS/RgEw1VVdf5mk28Fnj80/nzgvqHTTltzPPA/251O97baF5P8x+3pp+YOw0GzxYeBFUlOTLIg\nyV7tdtDFU11AVd3D4A/ke5Lsk2S3JL+YZKqnap7B4HrFRmD3JG/jye+4h9f1BPDnwAVJfg4G103a\nnUEA9wEHJNl3kvVtqqpHkxwNbPEPbpLfS3LdFqbtA6wFvlBV50zQ5FLgzCRLkywE3srg2sz4/Hu0\n6yq7te3ea/wiPoNrKs9ncF3lyFZbAXxiK9ulXYDhoFmhqjYApwBvYfDHeQPwP9j21+gZDE7RfA14\nEPgY8MwpzrsW+AxwO4OL24+y9VNSv8PgNtDrkzwMfI52faKqvsHgGsKd7drKRJ/leB1wXpIfMDin\nf8VW1nUI8IUtTPv3wC8Dr86TPz/y860vnwH+kMEpuu+0bTt3aP4/Z3AU9Argd9vwq9q891fVveOP\n1v5749dotOuKP/YjzX5JvgIcP8VTQdIOMxwkSR1PK0mSOoaDJKljOEiSOnP2Q3AHHnhgLVmyZNTd\nkKQ546abbvpeVS2avOUcDoclS5YwNjY26m5I0pyR5NuTtxrwtJIkqWM4SJI6hoMkqWM4SJI6UwqH\nDH7395YkX0ky1mr7J1mX5I72vF+rJ8n7kqxPcnOSo4aWs7K1vyPJyqH6C9vy17d5t+ubKSVJ02Nb\njhz+dVUdWVXL2vg5wDVVdRhwTRuHwS9gHdYeq4APwSBMGHzZ1zHA0cC544HS2rx2aL7l271FkqQd\ntiOnlU4BLmnDlzD4ucXx+qU1cD2wMMkzgROBde1Xsx4E1gHL27R9qur6GnzR06VDy5IkjcBUw6GA\nz7YfXF/Vage1788HuBc4qA0fzJO/5viuVtta/a4J6p0kq5KMJRnbuHHjFLsuSdpWU/0Q3Iuq6u72\noybrknxjeGJVVZKd/vWuVbUaWA2wbNkyv05WknaSKR05VNXd7fl+Br8AdTRwXzslRHu+vzW/m8EP\nk4xb3Gpbqy+eoC5Js8pxxx3HcccdN+puzIhJwyHJ05I8Y3wYOAH4Rwa/fTt+x9FK4Mo2fBVwRrtr\n6VjgoXb6aS1wQpL92oXoE4C1bdrDSY5tdymdMbQsSdIITOW00kHAJ9rdpbsDH6mqzyS5EbgiyZkM\nfnbwtNb+auBkBj+f+AjwaoCq2pTk7cCNrd15VbWpDb+OwW/a7g18uj0kSSMyaThU1Z0MfmB88/oD\nwPET1As4ewvLWgOsmaA+BjxvCv2VJM0APyEtSeoYDpKkjuEgSeoYDpKkzpz9JThJ89Ntt501snU/\n8sjtI+/D4YdfOCPr8chBktQxHCRJHcNBktQxHCRJHcNBktQxHCRJHcNBktTxcw6SNEWXXbZi1F2Y\nMR45SJI6hoMkqWM4SJI6hoMkqWM4aE6ZTz/wLo2S4SBJ6hgOkqSO4SBJ6hgOkqSO4SDNI17Q11QZ\nDpKkjuEgSeoYDpKkjuEgSerM2a/sfuCBB7j44otH3Q3NsCOPPBLAf/vttCvsv0ceOWDUXRipL37x\n4hlZz5wNB43O7Q/cPrJ1//ixH4+8D8894LkjW7c0U1JVo+7Ddlm2bFmNjY2Nuhvz0lmfOmtk6/7U\nWz4FwIp3ju5HVy5cceHI1r2jxm9jve6660bajx1x222je/3NBocfvv2vvyQ3VdWyqbT1moMkqWM4\nSJI6hoMkqTPlcEiyIMmXk/xNGz80yZeSrE/yV0n2bPWntPH1bfqSoWW8udVvS3LiUH15q61Pcs70\nbZ4kaXtsy5HDbwFfHxp/N3BBVT0HeBA4s9XPBB5s9QtaO5IsBU4HjgCWAx9sgbMA+ABwErAUeEVr\nK0kakSmFQ5LFwEuBv2jjAV4CfKw1uQQ4tQ2f0sZp049v7U8BLq+qn1TVPwHrgaPbY31V3VlVPwUu\nb20lSSMy1SOHPwbeBDzRxg8Avl9Vj7Xxu4CD2/DBwAaANv2h1v5n9c3m2VK9k2RVkrEkYxs3bpxi\n1yVJ22rScEjy74D7q+qmGejPVlXV6qpaVlXLFi1aNOruSNIuayqfkP4V4NeTnAzsBewD/AmwMMnu\n7ehgMXB3a383cAhwV5LdgX2BB4bq44bn2VJd2uWcNcLPcN1+++j7cOHc/QzhvDLpkUNVvbmqFlfV\nEgYXlD9fVa8ErgVe1pqtBK5sw1e1cdr0z9fgY9hXAae3u5kOBQ4DbgBuBA5rdz/t2dZx1bRsnXY5\nK965YqSfjpbmix35bqXfAS5P8g7gy8BFrX4RcFmS9cAmBn/sqapbk1wBfA14DDi7qh4HSPJ6YC2w\nAFhTVbfuQL8kSTtom8Khqq4DrmvDdzK402jzNo8CL9/C/OcD509Qvxq4elv6MlftCt9tI2nX5yek\nJUkdw0GS1DEcJEkdw0GS1DEcJEkdw0GS1DEcJEkdw0GS1DEcJEkdw0GS1DEcJEkdw0GS1NmRb2WV\nNMesWHHdqLugOcIjB0lSx3CQJHUMB0lSx3CQJHUMB0lSx3CQJHUMB0lSx3CQJHUMB0lSZ35+Qvqs\ns0a37ttvH30fLrxwdOuWNCd45CBJ6hgOkqSO4SBJ6hgOkqSO4SBJ6hgOkqSO4SBJ6hgOkqSO4SBJ\n6kwaDkn2SnJDkq8muTXJ77f6oUm+lGR9kr9KsmerP6WNr2/Tlwwt682tfluSE4fqy1ttfZJzpn8z\nJUnbYipHDj8BXlJVzweOBJYnORZ4N3BBVT0HeBA4s7U/E3iw1S9o7UiyFDgdOAJYDnwwyYIkC4AP\nACcBS4FXtLaSpBGZNBxq4IdtdI/2KOAlwMda/RLg1DZ8ShunTT8+SVr98qr6SVX9E7AeOLo91lfV\nnVX1U+Dy1laSNCJTuubQ3uF/BbgfWAd8E/h+VT3WmtwFHNyGDwY2ALTpDwEHDNc3m2dLdUnSiEwp\nHKrq8ao6EljM4J3+L+3UXm1BklVJxpKMbdy4cRRdkKR5YZvuVqqq7wPXAv8SWJhk/Cu/FwN3t+G7\ngUMA2vR9gQeG65vNs6X6ROtfXVXLqmrZokWLtqXrkqRtMJW7lRYlWdiG9wb+LfB1BiHxstZsJXBl\nG76qjdOmf76qqtVPb3czHQocBtwA3Agc1u5+2pPBReurpmPjJEnbZyo/9vNM4JJ2V9FuwBVV9TdJ\nvgZcnuQdwJeBi1r7i4DLkqwHNjH4Y09V3ZrkCuBrwGPA2VX1OECS1wNrgQXAmqq6ddq2UJK0zSYN\nh6q6GXjBBPU7GVx/2Lz+KPDyLSzrfOD8CepXA1dPob+SpBngJ6QlSR3DQZLUMRwkSR3DQZLUmcrd\nSppG161YMeouSNKkPHKQJHUMB0lSx3CQJHUMB0lSx3CQJHUMB0lSx3CQJHUMB0lSx3CQJHUMB0lS\nx3CQJHUMB0lSx3CQJHUMB0lSx3CQJHUMB0lSx3CQJHUMB0lSx3CQJHUMB0lSx3CQJHUMB0lSx3CQ\nJHUMB0lSx3CQJHUMB0lSx3CQJHUMB0lSZ9JwSHJIkmuTfC3JrUl+q9X3T7IuyR3teb9WT5L3JVmf\n5OYkRw0ta2Vrf0eSlUP1Fya5pc3zviTZGRsrSZqaqRw5PAb896paChwLnJ1kKXAOcE1VHQZc08YB\nTgIOa49VwIdgECbAucAxwNHAueOB0tq8dmi+5Tu+aZKk7TVpOFTVPVX1D234B8DXgYOBU4BLWrNL\ngFPb8CnApTVwPbAwyTOBE4F1VbWpqh4E1gHL27R9qur6qirg0qFlSZJGYJuuOSRZArwA+BJwUFXd\n0ybdCxzUhg8GNgzNdlerba1+1wT1ida/KslYkrGNGzduS9clSdtgyuGQ5OnAx4HfrqqHh6e1d/w1\nzX3rVNXqqlpWVcsWLVq0s1cnSfPWlMIhyR4MguEvq+p/tfJ97ZQQ7fn+Vr8bOGRo9sWttrX64gnq\nkqQRmcrdSgEuAr5eVe8dmnQVMH7H0UrgyqH6Ge2upWOBh9rpp7XACUn2axeiTwDWtmkPJzm2reuM\noWVJkkZg9ym0+RXgVcAtSb7Sam8B3gVckeRM4NvAaW3a1cDJwHrgEeDVAFW1KcnbgRtbu/OqalMb\nfh1wMbA38On2kCSNyKThUFV/D2zpcwfHT9C+gLO3sKw1wJoJ6mPA8ybriyRpZvgJaUlSx3CQJHUM\nB0lSx3CQJHUMB0lSx3CQJHUMB0lSx3CQJHUMB0lSx3CQJHUMB0lSx3CQJHUMB0lSx3CQJHUMB0lS\nx3CQJHUMB0lSx3CQJHUMB0lSx3CQJHUMB0lSx3CQJHUMB0lSx3CQJHUMB0lSx3CQJHUMB0lSx3CQ\nJHUMB0lSx3CQJHUMB0lSx3CQJHUMB0lSZ9JwSLImyf1J/nGotn+SdUnuaM/7tXqSvC/J+iQ3Jzlq\naJ6Vrf0dSVYO1V+Y5JY2z/uSZLo3UpK0baZy5HAxsHyz2jnANVV1GHBNGwc4CTisPVYBH4JBmADn\nAscARwPnjgdKa/Paofk2X5ckaYZNGg5V9XfAps3KpwCXtOFLgFOH6pfWwPXAwiTPBE4E1lXVpqp6\nEFgHLG/T9qmq66uqgEuHliVJGpHtveZwUFXd04bvBQ5qwwcDG4ba3dVqW6vfNUF9QklWJRlLMrZx\n48bt7LokaTI7fEG6veOvaejLVNa1uqqWVdWyRYsWzcQqJWle2t5wuK+dEqI939/qdwOHDLVb3Gpb\nqy+eoC5JGqHtDYergPE7jlYCVw7Vz2h3LR0LPNROP60FTkiyX7sQfQKwtk17OMmx7S6lM4aWJUka\nkd0na5Dko8BxwIFJ7mJw19G7gCuSnAl8GzitNb8aOBlYDzwCvBqgqjYleTtwY2t3XlWNX+R+HYM7\novYGPt0ekqQRmjQcquoVW5h0/ARtCzh7C8tZA6yZoD4GPG+yfkiSZo6fkJYkdQwHSVLHcJAkdQwH\nSVLHcJAkdQwHSVLHcJAkdQwHSVLHcJAkdQwHSVLHcJAkdQwHSVLHcJAkdQwHSVLHcJAkdQwHSVLH\ncJAkdQwHSVLHcJAkdQwHSVLHcJAkdQwHSVLHcJAkdQwHSVLHcJAkdQwHSVLHcJAkdQwHSVLHcJAk\ndQwHSVLHcJAkdQwHSVLHcJAkdWZNOCRZnuS2JOuTnDPq/kjSfDYrwiHJAuADwEnAUuAVSZaOtleS\nNH/NinAAjgbWV9WdVfVT4HLglBH3SZLmrVTVqPtAkpcBy6vqNW38VcAxVfX6zdqtAla10cOB22a0\no9PnQOB7o+7EHOb+2zHuvx0zl/ffs6tq0VQa7r6zezKdqmo1sHrU/dhRScaqatmo+zFXuf92jPtv\nx8yX/TdbTivdDRwyNL641SRJIzBbwuFG4LAkhybZEzgduGrEfZKkeWtWnFaqqseSvB5YCywA1lTV\nrSPu1s4050+NjZj7b8e4/3bMvNh/s+KCtCRpdpktp5UkSbOI4SBJ6syKaw5zXZK3Ab8KPFhVp7Xa\nc4E3A08Fvgu8tap+lGQP4C0MPgn+BPBHVXVTkr2AdzO4U+tx4P9U1ftnfmtm3nTsvzbP+xncg74A\n+Arwrqp6Yqa3Z6ZN1/4bWt57gcXjy9rVTePrbzWD19+jbdGvr6pNM7ox06mqfOzgAzgK+CXgiqHa\npcBRbfjXgd9sw6cB57bh/YEPMziC2wtY1up7AH8B/KtRb9tc2X9t/GntOcAfAieMetvm0v5rtZcA\n5w8va1d/TOPrbzWwdNTbM10PTytNg6r6B+DhzcrPBr7chr/E4D8dwKEMbt2lBu8qfgD8i6p6tKrG\nWv2fgW8AB+3krs8K07H/2viPWpsFDAJ2Xpiu/ZfkqcArgYt2cpdnlenaf7saw2Hn+Sbw4jb8b/j/\nf+jvAF6cZEGSZzF4YT0pBJI8A/g14IYZ6utstF37L8mfAuuAR4DPzVx3Z53t2X//hcE74UfR9v7/\nPTfJR5K8JklmrrvTz3DYec4DXp7kw8DTgMda/UrgPuAy4I3AzQzOXQI/+4ba84HLq2o+f0p8u/Zf\nDb6P60QGRw6/PJMdnmW2af+1c+yLq+raUXR2Ftqe199bq+o/AK8BXgCcPKM9nmZekN5JqupbwNkA\nSX4eeFGrPw68d7xdkjXAd4Zm/V1gQ1V9ZMY6OwvtwP6jqn6a5H8zeOf3pRnq8qyyHfvvKGBpkk8x\nOC23f5LVVbWKeWh7Xn9VdX97fiTJZ4AjgL+d0Y5PI48cdpIk+7fn3Ri8k/hYG98ryd5t+Bjg8aq6\ns42/Dng68J6RdHoW2db9l+SpSQ5s9QUM/jN/axR9nw22df9V1ceqanlVrWjtvzNfgwG26/W3IMnC\nVt+dwd1P3xxJ56eJRw7TIMk7gRcCC5NcDVwIPDXJy1uTa4FPteH9gT9N8gRwP/C2toyfA/4zgz9o\nH26nK6+oqk/O1HaMynTsPwZ3e723fTfXbsAY8PEZ2oSRmqb9N29N0/7bs9V3Z/D6uwH4xAxtwk7h\n12dIkjqeVpIkdQwHSVLHcJAkdQwHSVLHcJAkdQwHSVLHcJAkdf4fplY266m4yjcAAAAASUVORK5C\nYII=\n",
+      "application/javascript": [
+       "/* Put everything inside the global mpl namespace */\n",
+       "window.mpl = {};\n",
+       "\n",
+       "\n",
+       "mpl.get_websocket_type = function() {\n",
+       "    if (typeof(WebSocket) !== 'undefined') {\n",
+       "        return WebSocket;\n",
+       "    } else if (typeof(MozWebSocket) !== 'undefined') {\n",
+       "        return MozWebSocket;\n",
+       "    } else {\n",
+       "        alert('Your browser does not have WebSocket support. ' +\n",
+       "              'Please try Chrome, Safari or Firefox ≥ 6. ' +\n",
+       "              'Firefox 4 and 5 are also supported but you ' +\n",
+       "              'have to enable WebSockets in about:config.');\n",
+       "    };\n",
+       "}\n",
+       "\n",
+       "mpl.figure = function(figure_id, websocket, ondownload, parent_element) {\n",
+       "    this.id = figure_id;\n",
+       "\n",
+       "    this.ws = websocket;\n",
+       "\n",
+       "    this.supports_binary = (this.ws.binaryType != undefined);\n",
+       "\n",
+       "    if (!this.supports_binary) {\n",
+       "        var warnings = document.getElementById(\"mpl-warnings\");\n",
+       "        if (warnings) {\n",
+       "            warnings.style.display = 'block';\n",
+       "            warnings.textContent = (\n",
+       "                \"This browser does not support binary websocket messages. \" +\n",
+       "                    \"Performance may be slow.\");\n",
+       "        }\n",
+       "    }\n",
+       "\n",
+       "    this.imageObj = new Image();\n",
+       "\n",
+       "    this.context = undefined;\n",
+       "    this.message = undefined;\n",
+       "    this.canvas = undefined;\n",
+       "    this.rubberband_canvas = undefined;\n",
+       "    this.rubberband_context = undefined;\n",
+       "    this.format_dropdown = undefined;\n",
+       "\n",
+       "    this.image_mode = 'full';\n",
+       "\n",
+       "    this.root = $('<div/>');\n",
+       "    this._root_extra_style(this.root)\n",
+       "    this.root.attr('style', 'display: inline-block');\n",
+       "\n",
+       "    $(parent_element).append(this.root);\n",
+       "\n",
+       "    this._init_header(this);\n",
+       "    this._init_canvas(this);\n",
+       "    this._init_toolbar(this);\n",
+       "\n",
+       "    var fig = this;\n",
+       "\n",
+       "    this.waiting = false;\n",
+       "\n",
+       "    this.ws.onopen =  function () {\n",
+       "            fig.send_message(\"supports_binary\", {value: fig.supports_binary});\n",
+       "            fig.send_message(\"send_image_mode\", {});\n",
+       "            if (mpl.ratio != 1) {\n",
+       "                fig.send_message(\"set_dpi_ratio\", {'dpi_ratio': mpl.ratio});\n",
+       "            }\n",
+       "            fig.send_message(\"refresh\", {});\n",
+       "        }\n",
+       "\n",
+       "    this.imageObj.onload = function() {\n",
+       "            if (fig.image_mode == 'full') {\n",
+       "                // Full images could contain transparency (where diff images\n",
+       "                // almost always do), so we need to clear the canvas so that\n",
+       "                // there is no ghosting.\n",
+       "                fig.context.clearRect(0, 0, fig.canvas.width, fig.canvas.height);\n",
+       "            }\n",
+       "            fig.context.drawImage(fig.imageObj, 0, 0);\n",
+       "        };\n",
+       "\n",
+       "    this.imageObj.onunload = function() {\n",
+       "        fig.ws.close();\n",
+       "    }\n",
+       "\n",
+       "    this.ws.onmessage = this._make_on_message_function(this);\n",
+       "\n",
+       "    this.ondownload = ondownload;\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype._init_header = function() {\n",
+       "    var titlebar = $(\n",
+       "        '<div class=\"ui-dialog-titlebar ui-widget-header ui-corner-all ' +\n",
+       "        'ui-helper-clearfix\"/>');\n",
+       "    var titletext = $(\n",
+       "        '<div class=\"ui-dialog-title\" style=\"width: 100%; ' +\n",
+       "        'text-align: center; padding: 3px;\"/>');\n",
+       "    titlebar.append(titletext)\n",
+       "    this.root.append(titlebar);\n",
+       "    this.header = titletext[0];\n",
+       "}\n",
+       "\n",
+       "\n",
+       "\n",
+       "mpl.figure.prototype._canvas_extra_style = function(canvas_div) {\n",
+       "\n",
+       "}\n",
+       "\n",
+       "\n",
+       "mpl.figure.prototype._root_extra_style = function(canvas_div) {\n",
+       "\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype._init_canvas = function() {\n",
+       "    var fig = this;\n",
+       "\n",
+       "    var canvas_div = $('<div/>');\n",
+       "\n",
+       "    canvas_div.attr('style', 'position: relative; clear: both; outline: 0');\n",
+       "\n",
+       "    function canvas_keyboard_event(event) {\n",
+       "        return fig.key_event(event, event['data']);\n",
+       "    }\n",
+       "\n",
+       "    canvas_div.keydown('key_press', canvas_keyboard_event);\n",
+       "    canvas_div.keyup('key_release', canvas_keyboard_event);\n",
+       "    this.canvas_div = canvas_div\n",
+       "    this._canvas_extra_style(canvas_div)\n",
+       "    this.root.append(canvas_div);\n",
+       "\n",
+       "    var canvas = $('<canvas/>');\n",
+       "    canvas.addClass('mpl-canvas');\n",
+       "    canvas.attr('style', \"left: 0; top: 0; z-index: 0; outline: 0\")\n",
+       "\n",
+       "    this.canvas = canvas[0];\n",
+       "    this.context = canvas[0].getContext(\"2d\");\n",
+       "\n",
+       "    var backingStore = this.context.backingStorePixelRatio ||\n",
+       "\tthis.context.webkitBackingStorePixelRatio ||\n",
+       "\tthis.context.mozBackingStorePixelRatio ||\n",
+       "\tthis.context.msBackingStorePixelRatio ||\n",
+       "\tthis.context.oBackingStorePixelRatio ||\n",
+       "\tthis.context.backingStorePixelRatio || 1;\n",
+       "\n",
+       "    mpl.ratio = (window.devicePixelRatio || 1) / backingStore;\n",
+       "\n",
+       "    var rubberband = $('<canvas/>');\n",
+       "    rubberband.attr('style', \"position: absolute; left: 0; top: 0; z-index: 1;\")\n",
+       "\n",
+       "    var pass_mouse_events = true;\n",
+       "\n",
+       "    canvas_div.resizable({\n",
+       "        start: function(event, ui) {\n",
+       "            pass_mouse_events = false;\n",
+       "        },\n",
+       "        resize: function(event, ui) {\n",
+       "            fig.request_resize(ui.size.width, ui.size.height);\n",
+       "        },\n",
+       "        stop: function(event, ui) {\n",
+       "            pass_mouse_events = true;\n",
+       "            fig.request_resize(ui.size.width, ui.size.height);\n",
+       "        },\n",
+       "    });\n",
+       "\n",
+       "    function mouse_event_fn(event) {\n",
+       "        if (pass_mouse_events)\n",
+       "            return fig.mouse_event(event, event['data']);\n",
+       "    }\n",
+       "\n",
+       "    rubberband.mousedown('button_press', mouse_event_fn);\n",
+       "    rubberband.mouseup('button_release', mouse_event_fn);\n",
+       "    // Throttle sequential mouse events to 1 every 20ms.\n",
+       "    rubberband.mousemove('motion_notify', mouse_event_fn);\n",
+       "\n",
+       "    rubberband.mouseenter('figure_enter', mouse_event_fn);\n",
+       "    rubberband.mouseleave('figure_leave', mouse_event_fn);\n",
+       "\n",
+       "    canvas_div.on(\"wheel\", function (event) {\n",
+       "        event = event.originalEvent;\n",
+       "        event['data'] = 'scroll'\n",
+       "        if (event.deltaY < 0) {\n",
+       "            event.step = 1;\n",
+       "        } else {\n",
+       "            event.step = -1;\n",
+       "        }\n",
+       "        mouse_event_fn(event);\n",
+       "    });\n",
+       "\n",
+       "    canvas_div.append(canvas);\n",
+       "    canvas_div.append(rubberband);\n",
+       "\n",
+       "    this.rubberband = rubberband;\n",
+       "    this.rubberband_canvas = rubberband[0];\n",
+       "    this.rubberband_context = rubberband[0].getContext(\"2d\");\n",
+       "    this.rubberband_context.strokeStyle = \"#000000\";\n",
+       "\n",
+       "    this._resize_canvas = function(width, height) {\n",
+       "        // Keep the size of the canvas, canvas container, and rubber band\n",
+       "        // canvas in synch.\n",
+       "        canvas_div.css('width', width)\n",
+       "        canvas_div.css('height', height)\n",
+       "\n",
+       "        canvas.attr('width', width * mpl.ratio);\n",
+       "        canvas.attr('height', height * mpl.ratio);\n",
+       "        canvas.attr('style', 'width: ' + width + 'px; height: ' + height + 'px;');\n",
+       "\n",
+       "        rubberband.attr('width', width);\n",
+       "        rubberband.attr('height', height);\n",
+       "    }\n",
+       "\n",
+       "    // Set the figure to an initial 600x600px, this will subsequently be updated\n",
+       "    // upon first draw.\n",
+       "    this._resize_canvas(600, 600);\n",
+       "\n",
+       "    // Disable right mouse context menu.\n",
+       "    $(this.rubberband_canvas).bind(\"contextmenu\",function(e){\n",
+       "        return false;\n",
+       "    });\n",
+       "\n",
+       "    function set_focus () {\n",
+       "        canvas.focus();\n",
+       "        canvas_div.focus();\n",
+       "    }\n",
+       "\n",
+       "    window.setTimeout(set_focus, 100);\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype._init_toolbar = function() {\n",
+       "    var fig = this;\n",
+       "\n",
+       "    var nav_element = $('<div/>');\n",
+       "    nav_element.attr('style', 'width: 100%');\n",
+       "    this.root.append(nav_element);\n",
+       "\n",
+       "    // Define a callback function for later on.\n",
+       "    function toolbar_event(event) {\n",
+       "        return fig.toolbar_button_onclick(event['data']);\n",
+       "    }\n",
+       "    function toolbar_mouse_event(event) {\n",
+       "        return fig.toolbar_button_onmouseover(event['data']);\n",
+       "    }\n",
+       "\n",
+       "    for(var toolbar_ind in mpl.toolbar_items) {\n",
+       "        var name = mpl.toolbar_items[toolbar_ind][0];\n",
+       "        var tooltip = mpl.toolbar_items[toolbar_ind][1];\n",
+       "        var image = mpl.toolbar_items[toolbar_ind][2];\n",
+       "        var method_name = mpl.toolbar_items[toolbar_ind][3];\n",
+       "\n",
+       "        if (!name) {\n",
+       "            // put a spacer in here.\n",
+       "            continue;\n",
+       "        }\n",
+       "        var button = $('<button/>');\n",
+       "        button.addClass('ui-button ui-widget ui-state-default ui-corner-all ' +\n",
+       "                        'ui-button-icon-only');\n",
+       "        button.attr('role', 'button');\n",
+       "        button.attr('aria-disabled', 'false');\n",
+       "        button.click(method_name, toolbar_event);\n",
+       "        button.mouseover(tooltip, toolbar_mouse_event);\n",
+       "\n",
+       "        var icon_img = $('<span/>');\n",
+       "        icon_img.addClass('ui-button-icon-primary ui-icon');\n",
+       "        icon_img.addClass(image);\n",
+       "        icon_img.addClass('ui-corner-all');\n",
+       "\n",
+       "        var tooltip_span = $('<span/>');\n",
+       "        tooltip_span.addClass('ui-button-text');\n",
+       "        tooltip_span.html(tooltip);\n",
+       "\n",
+       "        button.append(icon_img);\n",
+       "        button.append(tooltip_span);\n",
+       "\n",
+       "        nav_element.append(button);\n",
+       "    }\n",
+       "\n",
+       "    var fmt_picker_span = $('<span/>');\n",
+       "\n",
+       "    var fmt_picker = $('<select/>');\n",
+       "    fmt_picker.addClass('mpl-toolbar-option ui-widget ui-widget-content');\n",
+       "    fmt_picker_span.append(fmt_picker);\n",
+       "    nav_element.append(fmt_picker_span);\n",
+       "    this.format_dropdown = fmt_picker[0];\n",
+       "\n",
+       "    for (var ind in mpl.extensions) {\n",
+       "        var fmt = mpl.extensions[ind];\n",
+       "        var option = $(\n",
+       "            '<option/>', {selected: fmt === mpl.default_extension}).html(fmt);\n",
+       "        fmt_picker.append(option);\n",
+       "    }\n",
+       "\n",
+       "    // Add hover states to the ui-buttons\n",
+       "    $( \".ui-button\" ).hover(\n",
+       "        function() { $(this).addClass(\"ui-state-hover\");},\n",
+       "        function() { $(this).removeClass(\"ui-state-hover\");}\n",
+       "    );\n",
+       "\n",
+       "    var status_bar = $('<span class=\"mpl-message\"/>');\n",
+       "    nav_element.append(status_bar);\n",
+       "    this.message = status_bar[0];\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype.request_resize = function(x_pixels, y_pixels) {\n",
+       "    // Request matplotlib to resize the figure. Matplotlib will then trigger a resize in the client,\n",
+       "    // which will in turn request a refresh of the image.\n",
+       "    this.send_message('resize', {'width': x_pixels, 'height': y_pixels});\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype.send_message = function(type, properties) {\n",
+       "    properties['type'] = type;\n",
+       "    properties['figure_id'] = this.id;\n",
+       "    this.ws.send(JSON.stringify(properties));\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype.send_draw_message = function() {\n",
+       "    if (!this.waiting) {\n",
+       "        this.waiting = true;\n",
+       "        this.ws.send(JSON.stringify({type: \"draw\", figure_id: this.id}));\n",
+       "    }\n",
+       "}\n",
+       "\n",
+       "\n",
+       "mpl.figure.prototype.handle_save = function(fig, msg) {\n",
+       "    var format_dropdown = fig.format_dropdown;\n",
+       "    var format = format_dropdown.options[format_dropdown.selectedIndex].value;\n",
+       "    fig.ondownload(fig, format);\n",
+       "}\n",
+       "\n",
+       "\n",
+       "mpl.figure.prototype.handle_resize = function(fig, msg) {\n",
+       "    var size = msg['size'];\n",
+       "    if (size[0] != fig.canvas.width || size[1] != fig.canvas.height) {\n",
+       "        fig._resize_canvas(size[0], size[1]);\n",
+       "        fig.send_message(\"refresh\", {});\n",
+       "    };\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype.handle_rubberband = function(fig, msg) {\n",
+       "    var x0 = msg['x0'] / mpl.ratio;\n",
+       "    var y0 = (fig.canvas.height - msg['y0']) / mpl.ratio;\n",
+       "    var x1 = msg['x1'] / mpl.ratio;\n",
+       "    var y1 = (fig.canvas.height - msg['y1']) / mpl.ratio;\n",
+       "    x0 = Math.floor(x0) + 0.5;\n",
+       "    y0 = Math.floor(y0) + 0.5;\n",
+       "    x1 = Math.floor(x1) + 0.5;\n",
+       "    y1 = Math.floor(y1) + 0.5;\n",
+       "    var min_x = Math.min(x0, x1);\n",
+       "    var min_y = Math.min(y0, y1);\n",
+       "    var width = Math.abs(x1 - x0);\n",
+       "    var height = Math.abs(y1 - y0);\n",
+       "\n",
+       "    fig.rubberband_context.clearRect(\n",
+       "        0, 0, fig.canvas.width / mpl.ratio, fig.canvas.height / mpl.ratio);\n",
+       "\n",
+       "    fig.rubberband_context.strokeRect(min_x, min_y, width, height);\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype.handle_figure_label = function(fig, msg) {\n",
+       "    // Updates the figure title.\n",
+       "    fig.header.textContent = msg['label'];\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype.handle_cursor = function(fig, msg) {\n",
+       "    var cursor = msg['cursor'];\n",
+       "    switch(cursor)\n",
+       "    {\n",
+       "    case 0:\n",
+       "        cursor = 'pointer';\n",
+       "        break;\n",
+       "    case 1:\n",
+       "        cursor = 'default';\n",
+       "        break;\n",
+       "    case 2:\n",
+       "        cursor = 'crosshair';\n",
+       "        break;\n",
+       "    case 3:\n",
+       "        cursor = 'move';\n",
+       "        break;\n",
+       "    }\n",
+       "    fig.rubberband_canvas.style.cursor = cursor;\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype.handle_message = function(fig, msg) {\n",
+       "    fig.message.textContent = msg['message'];\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype.handle_draw = function(fig, msg) {\n",
+       "    // Request the server to send over a new figure.\n",
+       "    fig.send_draw_message();\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype.handle_image_mode = function(fig, msg) {\n",
+       "    fig.image_mode = msg['mode'];\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype.updated_canvas_event = function() {\n",
+       "    // Called whenever the canvas gets updated.\n",
+       "    this.send_message(\"ack\", {});\n",
+       "}\n",
+       "\n",
+       "// A function to construct a web socket function for onmessage handling.\n",
+       "// Called in the figure constructor.\n",
+       "mpl.figure.prototype._make_on_message_function = function(fig) {\n",
+       "    return function socket_on_message(evt) {\n",
+       "        if (evt.data instanceof Blob) {\n",
+       "            /* FIXME: We get \"Resource interpreted as Image but\n",
+       "             * transferred with MIME type text/plain:\" errors on\n",
+       "             * Chrome.  But how to set the MIME type?  It doesn't seem\n",
+       "             * to be part of the websocket stream */\n",
+       "            evt.data.type = \"image/png\";\n",
+       "\n",
+       "            /* Free the memory for the previous frames */\n",
+       "            if (fig.imageObj.src) {\n",
+       "                (window.URL || window.webkitURL).revokeObjectURL(\n",
+       "                    fig.imageObj.src);\n",
+       "            }\n",
+       "\n",
+       "            fig.imageObj.src = (window.URL || window.webkitURL).createObjectURL(\n",
+       "                evt.data);\n",
+       "            fig.updated_canvas_event();\n",
+       "            fig.waiting = false;\n",
+       "            return;\n",
+       "        }\n",
+       "        else if (typeof evt.data === 'string' && evt.data.slice(0, 21) == \"data:image/png;base64\") {\n",
+       "            fig.imageObj.src = evt.data;\n",
+       "            fig.updated_canvas_event();\n",
+       "            fig.waiting = false;\n",
+       "            return;\n",
+       "        }\n",
+       "\n",
+       "        var msg = JSON.parse(evt.data);\n",
+       "        var msg_type = msg['type'];\n",
+       "\n",
+       "        // Call the  \"handle_{type}\" callback, which takes\n",
+       "        // the figure and JSON message as its only arguments.\n",
+       "        try {\n",
+       "            var callback = fig[\"handle_\" + msg_type];\n",
+       "        } catch (e) {\n",
+       "            console.log(\"No handler for the '\" + msg_type + \"' message type: \", msg);\n",
+       "            return;\n",
+       "        }\n",
+       "\n",
+       "        if (callback) {\n",
+       "            try {\n",
+       "                // console.log(\"Handling '\" + msg_type + \"' message: \", msg);\n",
+       "                callback(fig, msg);\n",
+       "            } catch (e) {\n",
+       "                console.log(\"Exception inside the 'handler_\" + msg_type + \"' callback:\", e, e.stack, msg);\n",
+       "            }\n",
+       "        }\n",
+       "    };\n",
+       "}\n",
+       "\n",
+       "// from http://stackoverflow.com/questions/1114465/getting-mouse-location-in-canvas\n",
+       "mpl.findpos = function(e) {\n",
+       "    //this section is from http://www.quirksmode.org/js/events_properties.html\n",
+       "    var targ;\n",
+       "    if (!e)\n",
+       "        e = window.event;\n",
+       "    if (e.target)\n",
+       "        targ = e.target;\n",
+       "    else if (e.srcElement)\n",
+       "        targ = e.srcElement;\n",
+       "    if (targ.nodeType == 3) // defeat Safari bug\n",
+       "        targ = targ.parentNode;\n",
+       "\n",
+       "    // jQuery normalizes the pageX and pageY\n",
+       "    // pageX,Y are the mouse positions relative to the document\n",
+       "    // offset() returns the position of the element relative to the document\n",
+       "    var x = e.pageX - $(targ).offset().left;\n",
+       "    var y = e.pageY - $(targ).offset().top;\n",
+       "\n",
+       "    return {\"x\": x, \"y\": y};\n",
+       "};\n",
+       "\n",
+       "/*\n",
+       " * return a copy of an object with only non-object keys\n",
+       " * we need this to avoid circular references\n",
+       " * http://stackoverflow.com/a/24161582/3208463\n",
+       " */\n",
+       "function simpleKeys (original) {\n",
+       "  return Object.keys(original).reduce(function (obj, key) {\n",
+       "    if (typeof original[key] !== 'object')\n",
+       "        obj[key] = original[key]\n",
+       "    return obj;\n",
+       "  }, {});\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype.mouse_event = function(event, name) {\n",
+       "    var canvas_pos = mpl.findpos(event)\n",
+       "\n",
+       "    if (name === 'button_press')\n",
+       "    {\n",
+       "        this.canvas.focus();\n",
+       "        this.canvas_div.focus();\n",
+       "    }\n",
+       "\n",
+       "    var x = canvas_pos.x * mpl.ratio;\n",
+       "    var y = canvas_pos.y * mpl.ratio;\n",
+       "\n",
+       "    this.send_message(name, {x: x, y: y, button: event.button,\n",
+       "                             step: event.step,\n",
+       "                             guiEvent: simpleKeys(event)});\n",
+       "\n",
+       "    /* This prevents the web browser from automatically changing to\n",
+       "     * the text insertion cursor when the button is pressed.  We want\n",
+       "     * to control all of the cursor setting manually through the\n",
+       "     * 'cursor' event from matplotlib */\n",
+       "    event.preventDefault();\n",
+       "    return false;\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype._key_event_extra = function(event, name) {\n",
+       "    // Handle any extra behaviour associated with a key event\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype.key_event = function(event, name) {\n",
+       "\n",
+       "    // Prevent repeat events\n",
+       "    if (name == 'key_press')\n",
+       "    {\n",
+       "        if (event.which === this._key)\n",
+       "            return;\n",
+       "        else\n",
+       "            this._key = event.which;\n",
+       "    }\n",
+       "    if (name == 'key_release')\n",
+       "        this._key = null;\n",
+       "\n",
+       "    var value = '';\n",
+       "    if (event.ctrlKey && event.which != 17)\n",
+       "        value += \"ctrl+\";\n",
+       "    if (event.altKey && event.which != 18)\n",
+       "        value += \"alt+\";\n",
+       "    if (event.shiftKey && event.which != 16)\n",
+       "        value += \"shift+\";\n",
+       "\n",
+       "    value += 'k';\n",
+       "    value += event.which.toString();\n",
+       "\n",
+       "    this._key_event_extra(event, name);\n",
+       "\n",
+       "    this.send_message(name, {key: value,\n",
+       "                             guiEvent: simpleKeys(event)});\n",
+       "    return false;\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype.toolbar_button_onclick = function(name) {\n",
+       "    if (name == 'download') {\n",
+       "        this.handle_save(this, null);\n",
+       "    } else {\n",
+       "        this.send_message(\"toolbar_button\", {name: name});\n",
+       "    }\n",
+       "};\n",
+       "\n",
+       "mpl.figure.prototype.toolbar_button_onmouseover = function(tooltip) {\n",
+       "    this.message.textContent = tooltip;\n",
+       "};\n",
+       "mpl.toolbar_items = [[\"Home\", \"Reset original view\", \"fa fa-home icon-home\", \"home\"], [\"Back\", \"Back to previous view\", \"fa fa-arrow-left icon-arrow-left\", \"back\"], [\"Forward\", \"Forward to next view\", \"fa fa-arrow-right icon-arrow-right\", \"forward\"], [\"\", \"\", \"\", \"\"], [\"Pan\", \"Pan axes with left mouse, zoom with right\", \"fa fa-arrows icon-move\", \"pan\"], [\"Zoom\", \"Zoom to rectangle\", \"fa fa-square-o icon-check-empty\", \"zoom\"], [\"\", \"\", \"\", \"\"], [\"Download\", \"Download plot\", \"fa fa-floppy-o icon-save\", \"download\"]];\n",
+       "\n",
+       "mpl.extensions = [\"eps\", \"jpeg\", \"pdf\", \"png\", \"ps\", \"raw\", \"svg\", \"tif\"];\n",
+       "\n",
+       "mpl.default_extension = \"png\";var comm_websocket_adapter = function(comm) {\n",
+       "    // Create a \"websocket\"-like object which calls the given IPython comm\n",
+       "    // object with the appropriate methods. Currently this is a non binary\n",
+       "    // socket, so there is still some room for performance tuning.\n",
+       "    var ws = {};\n",
+       "\n",
+       "    ws.close = function() {\n",
+       "        comm.close()\n",
+       "    };\n",
+       "    ws.send = function(m) {\n",
+       "        //console.log('sending', m);\n",
+       "        comm.send(m);\n",
+       "    };\n",
+       "    // Register the callback with on_msg.\n",
+       "    comm.on_msg(function(msg) {\n",
+       "        //console.log('receiving', msg['content']['data'], msg);\n",
+       "        // Pass the mpl event to the overridden (by mpl) onmessage function.\n",
+       "        ws.onmessage(msg['content']['data'])\n",
+       "    });\n",
+       "    return ws;\n",
+       "}\n",
+       "\n",
+       "mpl.mpl_figure_comm = function(comm, msg) {\n",
+       "    // This is the function which gets called when the mpl process\n",
+       "    // starts-up an IPython Comm through the \"matplotlib\" channel.\n",
+       "\n",
+       "    var id = msg.content.data.id;\n",
+       "    // Get hold of the div created by the display call when the Comm\n",
+       "    // socket was opened in Python.\n",
+       "    var element = $(\"#\" + id);\n",
+       "    var ws_proxy = comm_websocket_adapter(comm)\n",
+       "\n",
+       "    function ondownload(figure, format) {\n",
+       "        window.open(figure.imageObj.src);\n",
+       "    }\n",
+       "\n",
+       "    var fig = new mpl.figure(id, ws_proxy,\n",
+       "                           ondownload,\n",
+       "                           element.get(0));\n",
+       "\n",
+       "    // Call onopen now - mpl needs it, as it is assuming we've passed it a real\n",
+       "    // web socket which is closed, not our websocket->open comm proxy.\n",
+       "    ws_proxy.onopen();\n",
+       "\n",
+       "    fig.parent_element = element.get(0);\n",
+       "    fig.cell_info = mpl.find_output_cell(\"<div id='\" + id + \"'></div>\");\n",
+       "    if (!fig.cell_info) {\n",
+       "        console.error(\"Failed to find cell for figure\", id, fig);\n",
+       "        return;\n",
+       "    }\n",
+       "\n",
+       "    var output_index = fig.cell_info[2]\n",
+       "    var cell = fig.cell_info[0];\n",
+       "\n",
+       "};\n",
+       "\n",
+       "mpl.figure.prototype.handle_close = function(fig, msg) {\n",
+       "    var width = fig.canvas.width/mpl.ratio\n",
+       "    fig.root.unbind('remove')\n",
+       "\n",
+       "    // Update the output cell to use the data from the current canvas.\n",
+       "    fig.push_to_output();\n",
+       "    var dataURL = fig.canvas.toDataURL();\n",
+       "    // Re-enable the keyboard manager in IPython - without this line, in FF,\n",
+       "    // the notebook keyboard shortcuts fail.\n",
+       "    IPython.keyboard_manager.enable()\n",
+       "    $(fig.parent_element).html('<img src=\"' + dataURL + '\" width=\"' + width + '\">');\n",
+       "    fig.close_ws(fig, msg);\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype.close_ws = function(fig, msg){\n",
+       "    fig.send_message('closing', msg);\n",
+       "    // fig.ws.close()\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype.push_to_output = function(remove_interactive) {\n",
+       "    // Turn the data on the canvas into data in the output cell.\n",
+       "    var width = this.canvas.width/mpl.ratio\n",
+       "    var dataURL = this.canvas.toDataURL();\n",
+       "    this.cell_info[1]['text/html'] = '<img src=\"' + dataURL + '\" width=\"' + width + '\">';\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype.updated_canvas_event = function() {\n",
+       "    // Tell IPython that the notebook contents must change.\n",
+       "    IPython.notebook.set_dirty(true);\n",
+       "    this.send_message(\"ack\", {});\n",
+       "    var fig = this;\n",
+       "    // Wait a second, then push the new image to the DOM so\n",
+       "    // that it is saved nicely (might be nice to debounce this).\n",
+       "    setTimeout(function () { fig.push_to_output() }, 1000);\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype._init_toolbar = function() {\n",
+       "    var fig = this;\n",
+       "\n",
+       "    var nav_element = $('<div/>');\n",
+       "    nav_element.attr('style', 'width: 100%');\n",
+       "    this.root.append(nav_element);\n",
+       "\n",
+       "    // Define a callback function for later on.\n",
+       "    function toolbar_event(event) {\n",
+       "        return fig.toolbar_button_onclick(event['data']);\n",
+       "    }\n",
+       "    function toolbar_mouse_event(event) {\n",
+       "        return fig.toolbar_button_onmouseover(event['data']);\n",
+       "    }\n",
+       "\n",
+       "    for(var toolbar_ind in mpl.toolbar_items){\n",
+       "        var name = mpl.toolbar_items[toolbar_ind][0];\n",
+       "        var tooltip = mpl.toolbar_items[toolbar_ind][1];\n",
+       "        var image = mpl.toolbar_items[toolbar_ind][2];\n",
+       "        var method_name = mpl.toolbar_items[toolbar_ind][3];\n",
+       "\n",
+       "        if (!name) { continue; };\n",
+       "\n",
+       "        var button = $('<button class=\"btn btn-default\" href=\"#\" title=\"' + name + '\"><i class=\"fa ' + image + ' fa-lg\"></i></button>');\n",
+       "        button.click(method_name, toolbar_event);\n",
+       "        button.mouseover(tooltip, toolbar_mouse_event);\n",
+       "        nav_element.append(button);\n",
+       "    }\n",
+       "\n",
+       "    // Add the status bar.\n",
+       "    var status_bar = $('<span class=\"mpl-message\" style=\"text-align:right; float: right;\"/>');\n",
+       "    nav_element.append(status_bar);\n",
+       "    this.message = status_bar[0];\n",
+       "\n",
+       "    // Add the close button to the window.\n",
+       "    var buttongrp = $('<div class=\"btn-group inline pull-right\"></div>');\n",
+       "    var button = $('<button class=\"btn btn-mini btn-primary\" href=\"#\" title=\"Stop Interaction\"><i class=\"fa fa-power-off icon-remove icon-large\"></i></button>');\n",
+       "    button.click(function (evt) { fig.handle_close(fig, {}); } );\n",
+       "    button.mouseover('Stop Interaction', toolbar_mouse_event);\n",
+       "    buttongrp.append(button);\n",
+       "    var titlebar = this.root.find($('.ui-dialog-titlebar'));\n",
+       "    titlebar.prepend(buttongrp);\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype._root_extra_style = function(el){\n",
+       "    var fig = this\n",
+       "    el.on(\"remove\", function(){\n",
+       "\tfig.close_ws(fig, {});\n",
+       "    });\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype._canvas_extra_style = function(el){\n",
+       "    // this is important to make the div 'focusable\n",
+       "    el.attr('tabindex', 0)\n",
+       "    // reach out to IPython and tell the keyboard manager to turn it's self\n",
+       "    // off when our div gets focus\n",
+       "\n",
+       "    // location in version 3\n",
+       "    if (IPython.notebook.keyboard_manager) {\n",
+       "        IPython.notebook.keyboard_manager.register_events(el);\n",
+       "    }\n",
+       "    else {\n",
+       "        // location in version 2\n",
+       "        IPython.keyboard_manager.register_events(el);\n",
+       "    }\n",
+       "\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype._key_event_extra = function(event, name) {\n",
+       "    var manager = IPython.notebook.keyboard_manager;\n",
+       "    if (!manager)\n",
+       "        manager = IPython.keyboard_manager;\n",
+       "\n",
+       "    // Check for shift+enter\n",
+       "    if (event.shiftKey && event.which == 13) {\n",
+       "        this.canvas_div.blur();\n",
+       "        // select the cell after this one\n",
+       "        var index = IPython.notebook.find_cell_index(this.cell_info[0]);\n",
+       "        IPython.notebook.select(index + 1);\n",
+       "    }\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype.handle_save = function(fig, msg) {\n",
+       "    fig.ondownload(fig, null);\n",
+       "}\n",
+       "\n",
+       "\n",
+       "mpl.find_output_cell = function(html_output) {\n",
+       "    // Return the cell and output element which can be found *uniquely* in the notebook.\n",
+       "    // Note - this is a bit hacky, but it is done because the \"notebook_saving.Notebook\"\n",
+       "    // IPython event is triggered only after the cells have been serialised, which for\n",
+       "    // our purposes (turning an active figure into a static one), is too late.\n",
+       "    var cells = IPython.notebook.get_cells();\n",
+       "    var ncells = cells.length;\n",
+       "    for (var i=0; i<ncells; i++) {\n",
+       "        var cell = cells[i];\n",
+       "        if (cell.cell_type === 'code'){\n",
+       "            for (var j=0; j<cell.output_area.outputs.length; j++) {\n",
+       "                var data = cell.output_area.outputs[j];\n",
+       "                if (data.data) {\n",
+       "                    // IPython >= 3 moved mimebundle to data attribute of output\n",
+       "                    data = data.data;\n",
+       "                }\n",
+       "                if (data['text/html'] == html_output) {\n",
+       "                    return [cell, data, j];\n",
+       "                }\n",
+       "            }\n",
+       "        }\n",
+       "    }\n",
+       "}\n",
+       "\n",
+       "// Register the function which deals with the matplotlib target/channel.\n",
+       "// The kernel may be null if the page has been refreshed.\n",
+       "if (IPython.notebook.kernel != null) {\n",
+       "    IPython.notebook.kernel.comm_manager.register_target('matplotlib', mpl.mpl_figure_comm);\n",
+       "}\n"
+      ],
       "text/plain": [
-       "<matplotlib.figure.Figure at 0x7f4b461ecf60>"
+       "<IPython.core.display.Javascript object>"
       ]
      },
      "metadata": {},
@@ -66,8 +877,847 @@
     },
     {
      "data": {
+      "text/html": [
+       "<img src=\"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAdsAAAE8CAYAAACIOxsOAAAgAElEQVR4Xu2dDdRlV1nf/28GZoJhBkgAp5nEDsQg5gMrmHaMaBtQRBJYQPADDIQEFaHlI120CIWCClUEFWkRUQoSaOsXiqC2okKU2E79QNoMsSQoE+iEhJDBzGBgpuFN13+yT3K4c+9973vv89yzzz2/s9a7VibnnH32/u2P33n22efcNbFBAAIQgAAEIJBKYC01dRKHAAQgAAEIQEDIlkYAAQhAAAIQSCaAbJMBkzwEIAABCEAA2dIGIAABCEAAAskEkG0yYJKHAAQgAAEIIFvaAAQgAAEIQCCZALJNBkzyEIAABCAAAWRLG4AABCAAAQgkE0C2yYBJHgIQgAAEIIBsaQMQgAAEIACBZALINhkwyUMAAhCAAASQLW0AAhCAAAQgkEwA2SYDJnkIQAACEIAAsqUNQAACEIAABJIJINtkwCQPAQhAAAIQQLa0AQhAAAIQgEAyAWSbDJjkIQABCEAAAsiWNgABCEAAAhBIJoBskwGTPAQgAAEIQADZ0gZWkcB+Sa+W9MurWDjKdIxAU7fPhgcE+kAA2fahlurP41WS/mnJ5pck3Sjpw5J+VtL/2kT2Lch/Vv42cdpxhy4q292SPinpIZKcVl+250p6pqRHSvqsJJejvZ0g6V9J+iFJp0r6iKTnj9TR+ZJ+StI/knRQ0hskvamVyMslXSzp6yT9naTflOT/94UlQ9pItt8m6UcknSfpgRPq8iJJP17K8n8lvVLSr7bKcVppwxdIMjvvu0KS27g3i/4dI+X+bUlPXjILLtcDAsi2B5XUgyxatn9RBuatks6Q9AOSvrsMzO+fsQzIdkZQEw57kSTz31m4j8r2hZL+TZHE9ZIaOVuct0myXD4m6Y2S3iXprBJBWsi/Uq75u5L+k6S/lPTVkn5R0v+QdNliWd/02RvJ9rskfbOkT5c8jt44+YZkr6QXS/p9Sd8i6a2SHivpvxe5/lm5cbS0t0t6i6T/Kel5Ldla1hZ6s1nEvglhg8BXEEC2NIgIApat/yzL9vY2SY8vUcX/k/SkMth7ED8s6X2SXlKionFRggfIL0p6s6RHSzpJ0kclWSqOyiZtjkZ9jiMS//nfjuYcbTfbt0v6SUlnlwHZUbgHU293jiT8oyW/j5P0nWWfj/dAbeE4inyUpKsl3U/SUUkPlfRz5fqHJL1H0r8u5XESXyXp9ZK+R9K9S97+uaRPlfQtE/fPA2Vwv13S60aizEnlN0vXxahsLZHfkfTvyolO3xHda0rZzcj18bBWwj7WAvonEy7mGypL6uQp9WF5u+7/QZkxsKD+S+t4149ZuY4tSf/bQvvj1jEvkPQKSSdK+o8lWl0vNw7T2vCkWQqX65skuU6b7T9L2iLpewuDj0vaVYTrY3ys+T2o3JxM4jwtP+wbKAFkO9CKDy72JNl+Q5HjnhIRWCyWp6MnD2K/IOm/lgH+PmVKz9OYTy35u0XS6WWg/pCkI0VwT5R0Zms6b7Q4Hqx3lGnEP5FkiX2/pH9YpOlIzlGLZekB3fJ/e5lS/Y0SqXj/Py4i9hSppfyHku4v6cuS/rzcRFhQnkr1DcBTyhS4o0uX8b2SfMPhvPyHcs6/KJm9sgzanrr0jcdLy/Svpe30LVtz+PmSt++Q9O8lnVvSnlaFkyTgG5R3S/qZ1sl/W24SniXJka9nJB7R2v9vy/Sqb3R8EzG6WYqemvbNxaTNZXT0+DlJTTksuv9dTnB9uf6djqNkT0v7ON9s+SbNjyj+QJLZuT6dT0+X+wZmo2e2k2RrBl9b2laTb7cB34R9TWHgRyCnlOl0H+OpabcX33x8sFzbNxoulyNatw9HwZ8P7l8ktwIEkO0KVGIFRZgkW0chluv3jTwLa7L8tPJ8sBmoZ5lG9rMzT9P5eZsH3nGbB29HmZeUnY5WLJWfKIL3oOpo1ANjs3mA92DqSHzcAO3o09f1Mf+nnO+IzRGqn+P9uiRHQo6+LC5Pvfomo9l8E+HB2NKy9J3GgyU56vXWpO+B3NOblq2nOtvi89TvT5cyTKv2SbJ1JO3niU8oPHyjYJl/oETsvpaFfHmZKvZNyX8rNzx+xvuZkYs6ivfx5vnaTbRDR4e+mfmxlmz9qMHRqzffiDni/vrCyc9KLd12ffqZeiO8aZeeJNsLy02S29EflcjdNwTbyp9vmNxmnK9/WerZka+j22eUyNz16/R9Y+UpeLcv16dvDkZnRzaBh0NXkQCyXcVaXX6ZJsnWA5fv+D0t92tlas4Dkp9xPUDSvcqfReNtnGx9jKMri9kDvv9twTmy8bPDcZtl68ilvbDH0ecNRYyOSh0h3tE62el6CtdTqJMGaEfXXgBzTZl6dTTq6zhK86IwP7f0gG2pOWp2JN5s7mvOt0XyjWUQ99Rwe/N+C8WDumVrLo7Im81RlcXQSGpSTU+S7X3Lgh5HzJaBb0h8M+TN07fefrDk38c6QnMEbtn4OfDNrQtaRn5+25zbZjmaL8vJPHxT5RswtwuXz9fy5vqyrH+p/NusPEVsafmGyo8OLPTR+rTY5o1sfSm3q5cVzn6267q7tETZ3u8bpHeWfFv2nnr2I4VJN4+OiF0Wz4h4DQMbBO4mgGxpDBEENppG9vM+RzKO/Dx16EHTkaUHMw+iTTscJ1sPho4SPU17XZG3pxo95Tjp1R4PeI4APe3abG3Z/rUkTxd7IG1vHlAt5Emy9bNGTzm7DJaRp0dvLeVwhOcbCAvAz369SMz5Ht2cN6/m9cIiTxmPbhaap5XHLQCaxHk0jY2eJXrK1lHpTWUxkKdLm0U/Tsv14eernsb3DYRvKHy8p7e9+cbEkbxvfhyJT1uJ7Dr2TYqfBVvuLptnBDz12ohy3Opx3wz4ebvL7Px5Or5dn7/Vem46rQ1vtLLcsx6+kTAL168fBTiibm9+RuubEpfXbbh5LDLuum7XfmxhPmwQuJsAsqUxRBDYaIGUBzw/6/Tg7YhyX7moRepooWmHXinradxvbWXKU46ecvVg7c0DnqcYPdU5TbZeDOXo15unnj0l6AVRfk7slbX+f36GPG5rpjH9TO9vWgd4cZQjPcvWA78XeFkgloUl7Glfbz9cFoL5/HZ02yT1cEkWvp8DXzshD5mybS7pSMzls2DMedzmaM+C9MxCw9LPfV2Pjjz9etC0zfXmSNDRvzfXtcvsVb2zytazIubY1GfzWMASXySybefbaXrG4vdabW20XG6vfl5sbr4xG90sbc9wWMa+uWSDwN0EkC2NIYLALK/+eDCzbD1F2ky9+h1OL4Bq2qGnUL3Pz0Ud+Xgg96pfPydrplP9DqijJUeN02TrRUle/Wvp+lg/R/WzUkeejij/tKzubcTr6NtTo5axIzcLxtOnjoD/XpKnfP36h6dWPcXpyM9RrafFHWX72acHYm8+zkL+q3Iz4We9lrFX2zpP3ry4x89EfQ0/i/UAbvn7hsPHzyNbD/b+88pfC9/PZr15Gtab82DBOwr3lK5nGHzjYr7NM0az8tSt68v/7SlgT4v65sCbI0yn6z9Hg83mum0i33absshd55b1J8oiJ0/VeqZhVtk6wrX0nR/Xp5/tOo1pC6Q88+CbHd+cebrbz2gtQufBkbin6M3IC69cX35mb0n6ubXL4s314ellt0WX123PC8j8WpQ317vZOk1fx48PnK7TGceizYX/HhgBZDuwCk8qbvujFo5A2h+1aAZ6X9rPBT3AezGJp4I9jdu84uL9fmbpxTD+sIUHS69Gteh8nKMoT9F5ULTgpn0hqnn1xytLfZ7/7YHV+Ww2C91RtcXrKUJHNU7XC4K8eSrQi538ao+fkTavNVlUfu5ocXnzIOzBvHku3aRvsftmwqtqPQB7QHZZffPgzWn4OaWjPn90wa/4WCjNRxPmka3z+Koxddz0c4vEz7m9kts3HY4YvUisPQ3s6U/n2c9VvVDLQmk/f5y08GfSB0B8bb/W46jUNymWtW8InM6ssnWR/BjBNyLm5g9JeMp+2qs/bkOOfEe3ZmraN1Te7493ePN/+8bH9dRsvjFyffi1Jt8Qub34xqHZ/MqUb+LcRizoZjVy+yZkTBb4X0MkgGyHWOuUGQIQgAAElkoA2S4VNxeDAAQgAIEhEkC2Q6x1ygwBCEAAAkslgGyXipuLQQACEIDAEAkg2yHWOmWGAAQgAIGlEkC2S8XNxSAAAQhAYIgEBiXbbdu23fmlLzU/RTnE6qbMEIAABPpLYG1trbfO6m3G52kuW7ZsufOOO6Z9wnWeVDkHAhCAAASWQQDZLoNywDWQbQBEkoAABCDQEQFk2xH4zV4W2W6WGMdDAAIQqIcAsq2nLqbmBNn2pKLIJgQgAIExBJBtT5oFsu1JRZFNCEAAAsi2v20A2fa37sg5BCAAASLbnrQBZNuTiiKbEIAABIhsp7YB/5SZX2RtXmZ9afnZMP++41sl3af83Jl/t9Q/l+YtY9/YTCJb+i8EIACB/hIgsr2n7izbiyTta1Wn3+X1b0H6tyuvLr8R6h+uvrz8aHj0voktCdn2t5ORcwhAAALIdrpszys/9nxOOcw/lG0p+8fBM/YhW/okBCAAgRUkgGy/Ura3lYjVUezLJT22RLEXtur+dkmnSbogYd/BSW2MyHYFex9FggAEBkMA2d5T1adL+rSkbZLeKGm7pN9KEOo0Sbdle4Uk/x3b1tbWTl9fXx9Mw6SgEIDAcAmce+65xwp/zTXXrAwEZDu+Kl3T75P0PWOmkW+QdNKEaeRF9zGNvDJdi4JAAALzEkC285LLOS/yhwgsz3tJaqaRXyPpLEkXlwVSl7YWSJ0h6TJJJyTsQ7Y5bYVUIQCBHhFAtnVVVqRsvcL4PZK2lL9rJb1Q0mcknV9e/Tmx9erPzQVFxr6xlHlmW1fjIzcQgEAeAWSbx3aelCNlO8/1l3oOsl0qbi4GAQh0SADZdgh/zKWRbV31QW4gAAEIhBBAtiEYwxJBtmEoSQgCEIBAPQSQbT114Zwg27rqg9xAAAIQCCGAbEMwhiWCbMNQkhAEIACBeggg23rqgsi2rrogNxCAQM8JHDlyREePHq2iFHv2+DdepL1791aRn61bt2rbNn/vaP6Nj1rMz26pZ7Iaeam4uRgEBkXAoj311N06ePCmQZV71sKefPJO3Xjj/oWEi2xnpd3xcci24wrg8hBYYQKHDx/Wjh07JB0oX6rturB3RbZSDZHtYUm7dOjQIW3f7q/4zrch2/m4Lf0sZLt05FwQAoMhcI9sD1Ui27u+jSzV8G1ky3YHsh1Kb0C2Q6lpygmB5RNAttOYI1tWIy+/T3JFCEBgBQkgW2Q7jQCyXcFOT5EgAIHlE0C2yBbZFgJMIy9/AOKKEBgKAWSLbJEtsh3KeEc5IdAZAWSLbJEtsu1sAOLCEBgKAWSLbJEtsh3KeEc5IdAZAWSLbJEtsu1sAOLCEBgKgfpkWxN5Xv1hNXJN7ZG8QAACvSWAbIlsiWyJbHs7gHWV8VX8xZSuWA7lusgW2SJbZDuU8S6snMg2DOVgEkK2yBbZItvBDHhRBUW2USSHkw6yRbbIFtkOZ8QLKimyDQI5oGSQLbJFtsh2QENeTFGRbQzHIaWCbJEtskW2QxrzQsqKbEMwTk1k1RgjW2SLbJFt/si5YldYNRHUWD2rxhjZdi7bV0l6tST/kO+1kv5U0leVXH1G0nMl3VD+vUfSWyXdR9J+SZdI+uyC+6Z2M96zrXEUIk+dE1g1EXQOdEwGVo0xsu1Uto+U9FpJZ0m6UNI+SfeTdFvJ1YslfaukiyXZe9dLerakqyW9QtJDJV2+wL4Nuxiy3RARBwyRwKqJoMY6XDXGyLYz2W6TdJWkZ0j6kKSLimzbGXqlpG+Q9DRJ50l6h6RzygEPLNHtfRfYt2EXG5RsJd15zjkN3w3ZcMCACVx/vW98pTPPPHPAFHKLvmqM19fXde21nr10cHVCLrzepb5+bGb3rLPO0gknzM9m375945z1OkmfkvTmIs22bH9PkqPez0l6nKQbS3TrKNYRcLPdLuk0SReUCHez+w5uVCXIdiNC7B8kgVUTQY2VuGqMke20VpYm228u08ePdTA1RrbOlO3+Mkm7JD0f2S5hNODH45cAeUUusWpTnDVWy6oxZhq5k2nkH5H0QklHy9Udnd4s6TJJH2jlaGd5Trt9wlSxF06dtMC+DbvYoCJbZLthe+CAQmDVRFBjxa4aY2TbiWxHL+qVxZ5GvqksdrqlHPACSU+XdH6JdP2c6NLWAqkziqAdBc+zb8Muhmw3RMQBQySwaiKosQ5XjTGyrUq2nlK+UtK9i3Q/Kckrkv+25NLS9as/J7Ze/XFE7G3efVO7GbKtcRQacJ6OHDmio0ebGaHuQOzZ49fwpL1793aXiZErb926Vdu2eeHlamzIdjXqcbZS8Hu2yHa2lsJRSyBg0e7evVs33eQZILZRAjt37tT+/ftXRrjIdkhtHNki2yG198rL2kzDHThwQNu3ex1Dd1ttka3Z7Nq1S4cOHeqcTVStINsokn1IB9ki2z6004HksZFtDUKpTQQ1sYlqjrUxXrRcPLOdRhDZIttFexjnhxGoSSi1iaAmNlEVXhvjRcuFbJHtNAIZsm1/DNrfp8z44PO0NCeWl1d/Fh1Ocs+vSSi1iaAmNlGtoDbGi5YL2SLbZcp29GPQH0v44PO0j0hP7S/IdtHhJPf8moRSmwhqYhPVCmpjvGi5kC2yXZZsx30M2j9fFP3B52kfkUa2i44YHZ5fk1BqE0FNbKKaSG2MFy0XskW2y5LtuI9Bf92cH3VO+Rg0ke2iw0nu+TUJpTYR1MQmqhXUxnjRciFbZLsM2U76GHTXsr1Ckv+ObWtra6f7Y+FsdRKoSSi1iaAmNlGtpzbGi5YL2SLbZch20seg3yTpmSO/G7joB5/HTSM3aQ5mGnnVBipXXE1CqY1vFJtavtDl+q7tXeZFv9CFbJHtMmQ7eo3mY9D+ccd5Puqc8jHoVZpGrk0Gi0YFyHY6wQjZ8oWu6YwX/UIXskW2XcrWr/7M+1Hnec+bWF5kG6HEvDQihBKVu9puZiLY1PSFrtoi24gvdCFbZNuFbKPGvNB0kG0ozvDEIoQSlalVlm0NX+hyPdXEOKLtIVtki2wLAWQbpaKcdCIGvKic1SSCqCn2mvgi26iW2pd0+Fxjxhekqq19ZFtt1RzLWG0yqIlWBJuINCKZ1HRDE8GGyJbIlsiWyDZyjExLK2LAS8tcxwlHsIlIIxIDso2kWXtaRLZEtrW30Qn5q2mgikJYmwyiyhWRTgSbiDQiytKkUVMbjmBDZEtkS2RLZBs5RqalFTHgpWWu44Qj2ESkEYkB2UbSrD0tIlsi29rbKJFtT2soNtsRooxII7JUyDaSZu1pIVtkW3sbRbY9raHYbEeIMiKNyFIh20iataeFbJFt7W0U2fa0hmKzHSHKiDQiS4VsI2nWnhayRba1t1Fk29Mais12hCgj0ogsFbKNpFl7WsgW2dbeRpFtT2soNtsRooxII7JUyDaSZu1pIVtkW3sbRbY9raHYbEeIMiKNyFIh20iataeFbJFt7W0U2fa0hmKzHSHKiDQiS4VsI2nWnhayRba1t1Fk29Mais12hCgj0ogsFbKNpFl7WsgW2dbeRpFtT2soNtsRooxII7JUyDaSZu1pIVtkW3sbRbY9raHYbEeIMiKNyFIh20iataeFbJFt7W0U2fa0hmKzHSHKiDQiS4VsI2nWnhayRba1t1Fk29Mais12hCgj0ogsFbKNpFl7WsgW2dbeRpFtT2soNtsRooxII7JUyDaSZu1pIVtkW3sbRbY9raHYbEeIMiKN2FLVk1oEmyYN6ZCk7fUUroqcIFtkW0VD3HwmaooKNp/78WdEDHhReaktnQg2EWnUxiUqPxFskO202kC2yDaqty45HWS7ZOAdXy5SBocOHdL27URe7SqN5EtkO66zIFtk2/EgOu/lke285Pp5XqQMkO3xbSCSL7JFtuMIINt+jr1Ctj2tuDmzHSkDZIts52yGC5xGZItsF2g+XZ6KbLukv/xrI9tc5pF8iWyJbIlst2y584477sjttUtKHdkuCXQll4mUAZEtke3ymzWRLZHtJlvdkSNHdPTo0U2eFX/4nj17jiW6d+/e+MTnTHHr1q3atm3bnGdLEUKZ++KVnxjBJiKNyjHNnb0INqxGnoYf2SLbTXRPi3b3qafqpoMHN3HWcA7defLJ2n/jjXMLN2LAW1XaEWwi0oDvZALIFtlOI4BsNzF6NJ3pQAWvrN8V10q1xLW+b90laZEpSmSw8UAO30102E0cGtH2kC2yRbaFwJYFn9ne3SErkO25pUzXbGJAyTz0rkkiZJvFOFIGiwg7q3xdpxvJlwVS42qTaWQi2030cmQ7JfJCtptoSZs/NFIGyPZ4/pF8kS2yHUcA2W5i3EO2yHYTzSX00EgZIFtkG9o4Z0qMyBbZztRQ7joI2SLbTTSX0EORbSjO4xKL5EtkS2RLZMsz27QRi2e2aWi/8kZvge8aRwglt5TdpR7BhgVS0+qPyJbIdhP9m8iWyHYTzSX00EgZMI3MNHJo45wpMWSLbGdqKEwjb4SJyHYjQovtR7aL8dvo7Ei+TCMzjcw0MtPIG405c+9HtnOjm+nESBkQ2RLZztToQg8isiWy3USDYhqZaeRNNJfQQ5FtKM7jEovkS2RLZLuMyPa9kh4iaV3SFyS9QNJHJT1M0jslnSLpVknPknR9yVDGvrE9k49a5A1YRLZ5bJ1ypAyIbIlsc1srsl2GbO8n6bZyoSdLeqWkR0n6oKS3S3q3pEskXS7pMeW4jH3Idsm9CdnmAke2/eFLZItslyHb9jWeKelFkp4g6boS1X5Z0pYS3Z4pydPY0ftumdQtiWzzBixkm8eWyDaXbTRfZItslyXbt0l6XBHp4yWdKOlKSWe3MnBtiXAt2+h9H0G2+YPT6BWQbS5zItv+8EW2yHZZsm2u48j26WUqOVqo0yTdlu0Vkvx3bFtbWzt9fd2Pk+fbWCA1mRuyna9NzXoWsp2V1HzHRfJFtsh22bL19b4oabekjzONPN8gMOksfvUnlmftqUXKgAVSx9d2JF9ki2yzZXtSEeqnyoWeJOktkk6T9CFJnl5uFkg9R9IF5birEvaNHTt5ZpunFCLbPLZOOVIGyBbZ5rZWZJst2wdJer8kS9cLoQ5KeokkT+s+vLz68wBJny+v/jja9ZaxD9kuuTch21zgyLY/fIlskW22bHN7Q0DqRLYBECckgWzz2BLZ5rKN5otskS2y5XONaaMWsk1DeyxhItv+8EW2yBbZItu0EQvZpqFFtrlow/kiW2SLbJFt2rCFbNPQhsuABVLH11XkzAGyRbbIFtmmGQHZpqFFtrlow/kiW2SLbJFt2rCFbNPQhsuAyJbINre1Iltki2zT+hiyTUOLbHPRhvMlskW2yBbZpg1byDYNbbgMiGyJbHNbK7JFtsg2rY8h2zS0yDYXbThfIltki2yRbdqwhWzT0IbLgMiWyDa3tSJbZLtCsl1+Z5l+RWSbWyORr6YgW2Sb21qRLbJFtml9DNmmoSWyzUUbzpdpZGSLbJFt2rCFbNPQhsuAyJbINre1Iltki2zT+hiyTUOLbHPRhvMlskW2yBbZpg1byDYNbbgMiGyJbHNbK7JFtsg2rY8h2zS0yDYXbThfIltki2yRbdqwhWzT0IbLgMiWyDa3tSJbZIts0/oYsk1Di2xz0YbzJbJFtsgW2aYNW8g2DW24DIhsiWxzWyuyRbbINq2PIds0tMg2F204XyJbZItskW3asIVs09CGy4DIlsg2t7UiW2SLbNP6GLJNQ4tsc9GG8yWyRbbIFtmmDVvINg1tuAyIbIlsc1srskW2yDatjyHbNLTINhdtOF8iW2SLbJFt2rCFbNPQhsuAyJbINre1Iltki2zT+hiyTUOLbHPRhvMlskW2yBbZpg1byDYNbbgMiGyJbHNbK7JFtsg2rY8h2zS0yDYXbThfIltki2yRbdqwhWzT0IbLgMiWyDa3tSJbZIts0/oYsk1Di2xz0YbzJbJFtsgW2aYNW8g2DW24DIhsiWxzWyuyRbbINq2PIds0tMg2F204XyJbZItskW3asIVs09CGy4DIlsg2t7UiW2SLbNP6GLJNQ4tsc9GG8yWyRbbIFtmmDVvINg1tuAyIbIlsc1srskW2yDatjyHbNLTINhdtOF8iW2SLbJFt2rCFbNPQhsuAyJbINre1Itts2Z4i6V2SzpB0RNJ1kp4r6VZJeyS9VdJ9JO2XdImkz5YMZewb25a2INu0PoZs09Ai21y04XyJbJFttmxPlvQISVeVC71e0v0l/ZCk6yU9W9LVkl4h6aGSLpe0lrBvYtdEtnmjFrLNY+uUDx8+rB07dmiRqDQijdxSdpd6BJsmDWSLbLNlO5r+xZKeJ+llkt4h6ZxywANLdHtfSecl7EO2HYxZyDYXeqQMFhF2bim7Sz2SL7JFtsuU7QmSPiDpfZIOlCj2wlYGbpd0mqQLEvYdnNRliWzzBjNkm8eWyDaXbTRfZItslynbN0vaJempkp6SINRpkm7L9gpJ/ju2ra2tnb6+vj53z7377lfS9rlTWc0TkW1uvUZGXkS2x9dVJF9ki2yXJds3lGe3TywLpcZNFd8g6aQJ08iL7mMaOXfcH5s6ss2FHikDZItsc1srsl2GbF8r6XxJnjL2VLE3Tyl7gdSlrQVSXrF8WdI+ZLv8niRkmwsd2faHL5Etss2W7dmS9pVXfr5YLvYJSU8rAvarPye2Xv25uRxjOUfvG9szeWabN2Ah2zy2ThnZ9ocvskW22bLN7Q0BqSPbAIgTkkC2eWyRbS7baL7IFtkiWz5qkTZqIds0tMcSJrLtD19ki2yRLbJNG7GQbRpaZJuLNpwvskW2yBbZpg1byDYNbbgMWI18fF1FzhwgW2SLbJFtmhGQbRpaZJuLNpwvskW2yBbZpg1byDYNbbgMiGyJbHNbK5j758gAAAzeSURBVLJFtsg2rY8h2zS0yDYXbThfIltki2yRbdqwhWzT0IbLgMiWyDa3tSJbZIts0/oYsk1Di2xz0YbzJbJFtsgW2aYNW8g2DW24DIhsiWxzWyuyRbbINq2PIds0tMg2F204XyJbZItskW3asIVs09CGy4DIlsg2t7UiW2SLbNP6GLJNQ4tsc9GG8yWyRbbIFtmmDVvINg1tuAyIbIlsc1srskW2yDatjyHbNLTINhdtOF8iW2SLbJFt2rCFbNPQhsuAyJbINre1Iltki2zT+hiyTUOLbHPRhvMlskW2yBbZpg1byDYNbbgMiGyJbHNbK7JFtsg2rY8h2zS0yDYXbThfIltki2yRbdqwhWzT0IbLgMiWyDa3tSJbZIts0/oYsk1Di2xz0YbzJbJFtsgW2aYNW8g2DW24DIhsiWxzWyuyRbbINq2PIds0tMg2F204XyJbZItskW3asIVs09CGy4DIlsg2t7UiW2SLbNP6GLJNQ4tsc9GG8yWyRbbIFtmmDVvINg1tuAyIbIlsc1srskW2yDatjyHbNLTINhdtOF8iW2SLbJFt2rCFbNPQhsuAyJbINre1Iltki2zT+hiyTUOLbHPRhvMlskW2yBbZpg1byDYNbbgMiGyJbHNbK7JFtsg2rY8h2zS0yDYXbThfIltki2yRbdqwhWzT0IbLgMiWyDa3tSJbZIts0/oYsk1Di2xz0YbzJbJFtsgW2aYNW8g2DW24DIhsiWxzWyuyRbbINq2PIds0tMg2F204XyJbZItskW3asIVs09CGy4DIlsg2t7UiW2SLbNP6GLJNQ4tsc9GG8yWyRbbIFtmmDVvINg1tuAyIbIlsc1srss2W7RskXSxpt6RzJe0rF3yYpHdKOkXSrZKeJen6xH0T29EWZJvWx5BtGlpkm4s2nC+RLbLNlu2jJd0g6cOSLmrJ9oOS3i7p3ZIukXS5pMeUzGTsQ7ZLGJxGL4Fsc6EfPnxYO3bs0CJRaUQauaXsLvUINk0ayBbZZsu2SX9/S7YPlnRdiWq/LGlLiW7PlLSWsO+Wad2VyDZvMEO2eWydcqQMFhF2bim7Sz2SL7JFtl3I9lGSrpR0duvi15YI17KN3vcRZNvNgIVsc7lHygDZHl9XkXyRLbIdomyvkOS/Y9va2trp6+vrc4+Kd3dISdvnTmU1T0S2ufUaKQNki2xzWyuy7UK2TCMvv1V3ckVkm4sd2faHL5Etsu1Ctr7mVZLe1log9RxJF5TMZOyb2Ct5Zps3YCHbPLZOGdn2hy+yRbbZsn2TpKdK2inpc2UhlJ/VPry8+vMASZ8vr/58vGQmYx+yzR2XxqaObHOhI9v+8EW2yDZbtrm9ISB1ItsAiBOSQLZ5bIlsc9lG80W2yBbZ8lGLtFEL2aahPZYwkW1/+CJbZItskW3aiIVs09Ai21y04XyRLbJFtsg2bdhCtmlow2XAqz/H11XkzAGyRbbIFtmmGQHZpqFFtrlow/kiW2SLbJFt2rCFbNPQhsuAyJbINre1Iltki2zT+hiyTUOLbHPRhvMlskW2yBbZpg1byDYNbbgMiGyJbHNbK7JFtsg2rY8h2zS0yDYXbThfIltki2yRbdqwhWzT0IbLgMiWyDa3tSJbZIts0/oYsk1Di2xz0YbzJbJFtsgW2aYNW8g2DW24DIhsiWxzWyuyRbbINq2PIds0tMg2F204XyJbZItskW3asIVs09CGy4DIlsg2t7UiW2SLbNP6GLJNQ4tsc9GG8yWyRbbIFtmmDVvINg1tuAyIbIlsc1srskW2yDatjyHbNLTINhdtOF8iW2SLbJFt2rCFbNPQhsuAyJbINre1Iltki2zT+hiyTUOLbHPRhvMlskW2yBbZpg1byDYNbbgMiGyJbHNbK7JFtsg2rY8h2zS0yDYXbThfIltki2yRbdqwhWzT0IbLgMiWyDa3tSJbZIts0/oYsk1Di2xz0YbzJbJFtsgW2aYNW8g2DW24DIhsiWxzWyuyRbbINq2PIds0tMg2F204XyJbZItskW3asIVs09CGy4DIlsg2t7UiW2SLbNP6GLJNQ4tsc9GG8yWyRbbIFtmmDVvINg1tuAyIbIlsc1srskW2yDatjyHbNLTINhdtOF8iW2SLbJFt2rCFbNPQhsuAyJbINre1Iltki2zT+hiyTUOLbHPRhvMlskW2yBbZpg1byDYNbbgMiGyJbHNbK7JFtsg2rY8h2zS0yDYXbThfIltki2yRbdqwhWzT0IbLgMiWyDa3tSJbZIts0/oYsk1Di2xz0YbzJbJFtsgW2aYNW8g2DW24DIhsiWxzWyuyRbbINq2PIds0tMg2F204XyJbZLtqsn2YpHdKOkXSrZKeJen6af1yC7JNG7aQbRracBkQ2RLZ5rZWZLtqsv2gpLdLerekSyRdLukxyHb53chXRLa53A8fPqwdO3ZoEVFGpJFbyu5Sj2DTpEFki2xXSbYPlnRdiWq/LGlLiW7PlHTLpC5LZJs3mCHbPLbHbmaQbSrgSL7IFtmukmwfJelKSWe3CnVtiXA/ki3bA5K2p3b9/iVu2e6SQiKvAwcOaPt2CLdbgWWwa9cu+CZ1jUi+EiPEmIn6YyPEIjMzTnNtbW0tqQmkJ9vXjM8q2ysk+a/ZTpfkSHhVthMkra9KYSosB3zzKwXGuYxXje8dkk7MRZaTel9lO9c0cg7CTlP9lKSv6TQHq31x+ObXL4xzGcM3l+/MqfdVti7gVZLe1log9RxJF8xc8tU4kI6UW4/wzeXr1GGcyxi+uXxnTr3Psn14efXnAZI+X179+fjMJV+NA+lIufUI31y+yBa++QQquUKfZVsJwk6z4efRP9tpDlb74vDNr18Y5zKGby7fmVNHtjOj4kAIQAACEIDAfASQ7XzcOAsCEIAABCAwMwFkOzMqDoQABCAAAQjMRwDZzsct66w3SLpY0m5J50raVy50oaQfl3RvSQclPVvSJzfY529Gv0vSGZKOlC9uPbd8aSsr/7WnG8nXZX2vpIeUd52/IOkFkj5aO4TE/EXzbbL6KkmvHukTicWoNulovvslfan8udAvlfT71Za+5xlDtnVV4KMl3SDpw5IuKrL1amv/wML5RZhPl3SppMdLmrbvZEmPKK9IuZSvl3R/ST9YV5GXmptIvs74/STdVkrwZEmvlOQPrgx1i+Zrjo+U9FpJZ0nyTWdzAzpExtF8LdtmnBkiz6WWGdkuFffMF2t3gvMk/XLr05SWqH/l6EElqpq073MjV3PE/DxJ3z5zLlb3wAy+z5T0IknftLrYZi5ZFN9t5WbxGZI+hBju5h/FF9nO3KQXPxDZLs4wI4V2J3D09DeSvkvSn0t6vqQ3lwjK/3/SvvY3ov3Jtg9Iep+kN2VkuGdpRvL1h1Ue58+2ltmGj/WMRUZ2o/i+rnz0wu0dMdxTU1F8nY5nZtx2r5b0ckl/l9EgSPMuyGz1ERgdWPzTgT9avgn6uyWC+jZJ15SfFZy0rymZByv/TsBT+ZbyMSTRfJ2mI1tP8T+hvua09BxF8L1vmT5+rKQ7ke1X1GEEX48d/lb8pyV5BuGN5fdV/HOlbAkEkG0C1IAkp93Ff3V5rusFUH8/cq1x+7yows9un1gWSgVkr/dJRPJtw/iipNMGvght3M1Mm9Gs7deLzV4o6Wg52VxvlnRZmaXpfSNcoAAZ7dcLMj3z5QV/bAkEkG0C1IAkRzvTTkk3SfJ08C+V32p/cbnOtH1eWOKFVV5YcntAvlYliQi+J5XfU/YnHb09SdJbimwdiQ15i+A7yo9p5HuIRPB1+71Xaxr5NWUR2lOG3HAzy45sM+luPm0/T/VUrwXqBU5eCOXf7PVzwW+RtLXc1fsTbF6y723SPp/nlZvXSXLE5e0Tkp62+WytzBmRfL1A7f2SPGj5Zxv9StZLJE38PeWVoTi5IJF8ke3xnCP5PlTSeyRtKX/+PXDPJHxmAO20kyIi206wc1EIQAACEBgSAWQ7pNqmrBCAAAQg0AkBZNsJdi4KAQhAAAJDIoBsh1TblBUCEIAABDohgGw7wc5FIQABCEBgSASQ7ZBqm7JCAAIQgEAnBJBtJ9i5KAQgAAEIDIkAsh1SbVNWCEAAAhDohACy7QQ7F4UABCAAgSERQLZDqm3KCgEIQAACnRBAtp1g56IQgAAEIDAkAsh2SLVNWSEAAQhAoBMCyLYT7FwUAhCAAASGRADZDqm2KSsEIAABCHRCANl2gp2LQgACEIDAkAgg2yHVNmWFAAQgAIFOCCDbTrBzUQhAAAIQGBIBZDuk2qasEIAABCDQCQFk2wl2LgoBCEAAAkMigGyHVNuUFQIQgAAEOiGAbDvBzkUhAAEIQGBIBJDtkGqbskIAAhCAQCcEkG0n2LkoBCAAAQgMiQCyHVJtU1YIQAACEOiEwP8HUIFulfa6JBIAAAAASUVORK5CYII=\" width=\"431.818172458775\">"
+      ],
+      "text/plain": [
+       "<IPython.core.display.HTML object>"
+      ]
+     },
+     "metadata": {},
+     "output_type": "display_data"
+    }
+   ],
+   "source": [
+    "#Initialize figure and axes:\n",
+    "fig, ax = plt.subplots()\n",
+    "\n",
+    "#Set thereshold (manual input):\n",
+    "threshold = 43000\n",
+    "\n",
+    "#Set the bar colour:\n",
+    "colour_list=[]\n",
+    "for l,u in zip(lower_bound,upper_bound):\n",
+    "    if threshold < l:\n",
+    "        colour_list.append(\"blue\")\n",
+    "    elif threshold > u:\n",
+    "        colour_list.append(\"red\")\n",
+    "    else:\n",
+    "        colour_list.append(\"white\")\n",
+    "\n",
+    "def plot(colour_list):\n",
+    "    #Set plotting:   \n",
+    "    plot = ax.bar(df.index, year_avg, align='center',yerr= yerr,color = colour_list, edgecolor =\"black\")\n",
+    "\n",
+    "    #set thereshold line as a 2nd y axes and keep the thereshold value:\n",
+    "    ax2 = ax.twinx() \n",
+    "    ax2.set_ylim(ax.get_ylim())                           \n",
+    "    lhor = ax2.axhline(y = threshold, zorder = 1, color = 'black') \n",
+    "    ax2.tick_params(axis = 'both', which = 'both', length = 0) \n",
+    "    ytick = plt.gca().get_yticks()\n",
+    "    ytick = []\n",
+    "    ytick = np.append(ytick,threshold)\n",
+    "    plt.gca().set_yticks(ytick)\n",
+    "\n",
+    "    # Add xticks, titles and remove left and top edges:\n",
+    "    plt.xticks(df.index, ('1992', '1993', '1994', '1995'))\n",
+    "    plt.title('Data between 1992 and 1995')\n",
+    "    ax_list = fig.axes # list of axis\n",
+    "    for i in ax_list:\n",
+    "        [i.spines[loc].set_visible(False) for loc in ['top', 'right']]\n",
+    "    \n",
+    "    #return plot, lhor\n",
+    "\n",
+    "plot(colour_list)\n",
+    "\n",
+    "#show figure:\n",
+    "plt.show()"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "**Harder option implementation:** Implement the bar coloring as described in the paper, where the color of the bar is actually based on the amount of data covered (e.g. a gradient ranging from dark blue for the distribution being certainly below this y-axis, to white if the value is certainly contained, to dark red if the value is certainly not contained as the distribution is above the axis)."
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 244,
+   "metadata": {},
+   "outputs": [
+    {
+     "data": {
+      "application/javascript": [
+       "/* Put everything inside the global mpl namespace */\n",
+       "window.mpl = {};\n",
+       "\n",
+       "\n",
+       "mpl.get_websocket_type = function() {\n",
+       "    if (typeof(WebSocket) !== 'undefined') {\n",
+       "        return WebSocket;\n",
+       "    } else if (typeof(MozWebSocket) !== 'undefined') {\n",
+       "        return MozWebSocket;\n",
+       "    } else {\n",
+       "        alert('Your browser does not have WebSocket support. ' +\n",
+       "              'Please try Chrome, Safari or Firefox ≥ 6. ' +\n",
+       "              'Firefox 4 and 5 are also supported but you ' +\n",
+       "              'have to enable WebSockets in about:config.');\n",
+       "    };\n",
+       "}\n",
+       "\n",
+       "mpl.figure = function(figure_id, websocket, ondownload, parent_element) {\n",
+       "    this.id = figure_id;\n",
+       "\n",
+       "    this.ws = websocket;\n",
+       "\n",
+       "    this.supports_binary = (this.ws.binaryType != undefined);\n",
+       "\n",
+       "    if (!this.supports_binary) {\n",
+       "        var warnings = document.getElementById(\"mpl-warnings\");\n",
+       "        if (warnings) {\n",
+       "            warnings.style.display = 'block';\n",
+       "            warnings.textContent = (\n",
+       "                \"This browser does not support binary websocket messages. \" +\n",
+       "                    \"Performance may be slow.\");\n",
+       "        }\n",
+       "    }\n",
+       "\n",
+       "    this.imageObj = new Image();\n",
+       "\n",
+       "    this.context = undefined;\n",
+       "    this.message = undefined;\n",
+       "    this.canvas = undefined;\n",
+       "    this.rubberband_canvas = undefined;\n",
+       "    this.rubberband_context = undefined;\n",
+       "    this.format_dropdown = undefined;\n",
+       "\n",
+       "    this.image_mode = 'full';\n",
+       "\n",
+       "    this.root = $('<div/>');\n",
+       "    this._root_extra_style(this.root)\n",
+       "    this.root.attr('style', 'display: inline-block');\n",
+       "\n",
+       "    $(parent_element).append(this.root);\n",
+       "\n",
+       "    this._init_header(this);\n",
+       "    this._init_canvas(this);\n",
+       "    this._init_toolbar(this);\n",
+       "\n",
+       "    var fig = this;\n",
+       "\n",
+       "    this.waiting = false;\n",
+       "\n",
+       "    this.ws.onopen =  function () {\n",
+       "            fig.send_message(\"supports_binary\", {value: fig.supports_binary});\n",
+       "            fig.send_message(\"send_image_mode\", {});\n",
+       "            if (mpl.ratio != 1) {\n",
+       "                fig.send_message(\"set_dpi_ratio\", {'dpi_ratio': mpl.ratio});\n",
+       "            }\n",
+       "            fig.send_message(\"refresh\", {});\n",
+       "        }\n",
+       "\n",
+       "    this.imageObj.onload = function() {\n",
+       "            if (fig.image_mode == 'full') {\n",
+       "                // Full images could contain transparency (where diff images\n",
+       "                // almost always do), so we need to clear the canvas so that\n",
+       "                // there is no ghosting.\n",
+       "                fig.context.clearRect(0, 0, fig.canvas.width, fig.canvas.height);\n",
+       "            }\n",
+       "            fig.context.drawImage(fig.imageObj, 0, 0);\n",
+       "        };\n",
+       "\n",
+       "    this.imageObj.onunload = function() {\n",
+       "        fig.ws.close();\n",
+       "    }\n",
+       "\n",
+       "    this.ws.onmessage = this._make_on_message_function(this);\n",
+       "\n",
+       "    this.ondownload = ondownload;\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype._init_header = function() {\n",
+       "    var titlebar = $(\n",
+       "        '<div class=\"ui-dialog-titlebar ui-widget-header ui-corner-all ' +\n",
+       "        'ui-helper-clearfix\"/>');\n",
+       "    var titletext = $(\n",
+       "        '<div class=\"ui-dialog-title\" style=\"width: 100%; ' +\n",
+       "        'text-align: center; padding: 3px;\"/>');\n",
+       "    titlebar.append(titletext)\n",
+       "    this.root.append(titlebar);\n",
+       "    this.header = titletext[0];\n",
+       "}\n",
+       "\n",
+       "\n",
+       "\n",
+       "mpl.figure.prototype._canvas_extra_style = function(canvas_div) {\n",
+       "\n",
+       "}\n",
+       "\n",
+       "\n",
+       "mpl.figure.prototype._root_extra_style = function(canvas_div) {\n",
+       "\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype._init_canvas = function() {\n",
+       "    var fig = this;\n",
+       "\n",
+       "    var canvas_div = $('<div/>');\n",
+       "\n",
+       "    canvas_div.attr('style', 'position: relative; clear: both; outline: 0');\n",
+       "\n",
+       "    function canvas_keyboard_event(event) {\n",
+       "        return fig.key_event(event, event['data']);\n",
+       "    }\n",
+       "\n",
+       "    canvas_div.keydown('key_press', canvas_keyboard_event);\n",
+       "    canvas_div.keyup('key_release', canvas_keyboard_event);\n",
+       "    this.canvas_div = canvas_div\n",
+       "    this._canvas_extra_style(canvas_div)\n",
+       "    this.root.append(canvas_div);\n",
+       "\n",
+       "    var canvas = $('<canvas/>');\n",
+       "    canvas.addClass('mpl-canvas');\n",
+       "    canvas.attr('style', \"left: 0; top: 0; z-index: 0; outline: 0\")\n",
+       "\n",
+       "    this.canvas = canvas[0];\n",
+       "    this.context = canvas[0].getContext(\"2d\");\n",
+       "\n",
+       "    var backingStore = this.context.backingStorePixelRatio ||\n",
+       "\tthis.context.webkitBackingStorePixelRatio ||\n",
+       "\tthis.context.mozBackingStorePixelRatio ||\n",
+       "\tthis.context.msBackingStorePixelRatio ||\n",
+       "\tthis.context.oBackingStorePixelRatio ||\n",
+       "\tthis.context.backingStorePixelRatio || 1;\n",
+       "\n",
+       "    mpl.ratio = (window.devicePixelRatio || 1) / backingStore;\n",
+       "\n",
+       "    var rubberband = $('<canvas/>');\n",
+       "    rubberband.attr('style', \"position: absolute; left: 0; top: 0; z-index: 1;\")\n",
+       "\n",
+       "    var pass_mouse_events = true;\n",
+       "\n",
+       "    canvas_div.resizable({\n",
+       "        start: function(event, ui) {\n",
+       "            pass_mouse_events = false;\n",
+       "        },\n",
+       "        resize: function(event, ui) {\n",
+       "            fig.request_resize(ui.size.width, ui.size.height);\n",
+       "        },\n",
+       "        stop: function(event, ui) {\n",
+       "            pass_mouse_events = true;\n",
+       "            fig.request_resize(ui.size.width, ui.size.height);\n",
+       "        },\n",
+       "    });\n",
+       "\n",
+       "    function mouse_event_fn(event) {\n",
+       "        if (pass_mouse_events)\n",
+       "            return fig.mouse_event(event, event['data']);\n",
+       "    }\n",
+       "\n",
+       "    rubberband.mousedown('button_press', mouse_event_fn);\n",
+       "    rubberband.mouseup('button_release', mouse_event_fn);\n",
+       "    // Throttle sequential mouse events to 1 every 20ms.\n",
+       "    rubberband.mousemove('motion_notify', mouse_event_fn);\n",
+       "\n",
+       "    rubberband.mouseenter('figure_enter', mouse_event_fn);\n",
+       "    rubberband.mouseleave('figure_leave', mouse_event_fn);\n",
+       "\n",
+       "    canvas_div.on(\"wheel\", function (event) {\n",
+       "        event = event.originalEvent;\n",
+       "        event['data'] = 'scroll'\n",
+       "        if (event.deltaY < 0) {\n",
+       "            event.step = 1;\n",
+       "        } else {\n",
+       "            event.step = -1;\n",
+       "        }\n",
+       "        mouse_event_fn(event);\n",
+       "    });\n",
+       "\n",
+       "    canvas_div.append(canvas);\n",
+       "    canvas_div.append(rubberband);\n",
+       "\n",
+       "    this.rubberband = rubberband;\n",
+       "    this.rubberband_canvas = rubberband[0];\n",
+       "    this.rubberband_context = rubberband[0].getContext(\"2d\");\n",
+       "    this.rubberband_context.strokeStyle = \"#000000\";\n",
+       "\n",
+       "    this._resize_canvas = function(width, height) {\n",
+       "        // Keep the size of the canvas, canvas container, and rubber band\n",
+       "        // canvas in synch.\n",
+       "        canvas_div.css('width', width)\n",
+       "        canvas_div.css('height', height)\n",
+       "\n",
+       "        canvas.attr('width', width * mpl.ratio);\n",
+       "        canvas.attr('height', height * mpl.ratio);\n",
+       "        canvas.attr('style', 'width: ' + width + 'px; height: ' + height + 'px;');\n",
+       "\n",
+       "        rubberband.attr('width', width);\n",
+       "        rubberband.attr('height', height);\n",
+       "    }\n",
+       "\n",
+       "    // Set the figure to an initial 600x600px, this will subsequently be updated\n",
+       "    // upon first draw.\n",
+       "    this._resize_canvas(600, 600);\n",
+       "\n",
+       "    // Disable right mouse context menu.\n",
+       "    $(this.rubberband_canvas).bind(\"contextmenu\",function(e){\n",
+       "        return false;\n",
+       "    });\n",
+       "\n",
+       "    function set_focus () {\n",
+       "        canvas.focus();\n",
+       "        canvas_div.focus();\n",
+       "    }\n",
+       "\n",
+       "    window.setTimeout(set_focus, 100);\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype._init_toolbar = function() {\n",
+       "    var fig = this;\n",
+       "\n",
+       "    var nav_element = $('<div/>');\n",
+       "    nav_element.attr('style', 'width: 100%');\n",
+       "    this.root.append(nav_element);\n",
+       "\n",
+       "    // Define a callback function for later on.\n",
+       "    function toolbar_event(event) {\n",
+       "        return fig.toolbar_button_onclick(event['data']);\n",
+       "    }\n",
+       "    function toolbar_mouse_event(event) {\n",
+       "        return fig.toolbar_button_onmouseover(event['data']);\n",
+       "    }\n",
+       "\n",
+       "    for(var toolbar_ind in mpl.toolbar_items) {\n",
+       "        var name = mpl.toolbar_items[toolbar_ind][0];\n",
+       "        var tooltip = mpl.toolbar_items[toolbar_ind][1];\n",
+       "        var image = mpl.toolbar_items[toolbar_ind][2];\n",
+       "        var method_name = mpl.toolbar_items[toolbar_ind][3];\n",
+       "\n",
+       "        if (!name) {\n",
+       "            // put a spacer in here.\n",
+       "            continue;\n",
+       "        }\n",
+       "        var button = $('<button/>');\n",
+       "        button.addClass('ui-button ui-widget ui-state-default ui-corner-all ' +\n",
+       "                        'ui-button-icon-only');\n",
+       "        button.attr('role', 'button');\n",
+       "        button.attr('aria-disabled', 'false');\n",
+       "        button.click(method_name, toolbar_event);\n",
+       "        button.mouseover(tooltip, toolbar_mouse_event);\n",
+       "\n",
+       "        var icon_img = $('<span/>');\n",
+       "        icon_img.addClass('ui-button-icon-primary ui-icon');\n",
+       "        icon_img.addClass(image);\n",
+       "        icon_img.addClass('ui-corner-all');\n",
+       "\n",
+       "        var tooltip_span = $('<span/>');\n",
+       "        tooltip_span.addClass('ui-button-text');\n",
+       "        tooltip_span.html(tooltip);\n",
+       "\n",
+       "        button.append(icon_img);\n",
+       "        button.append(tooltip_span);\n",
+       "\n",
+       "        nav_element.append(button);\n",
+       "    }\n",
+       "\n",
+       "    var fmt_picker_span = $('<span/>');\n",
+       "\n",
+       "    var fmt_picker = $('<select/>');\n",
+       "    fmt_picker.addClass('mpl-toolbar-option ui-widget ui-widget-content');\n",
+       "    fmt_picker_span.append(fmt_picker);\n",
+       "    nav_element.append(fmt_picker_span);\n",
+       "    this.format_dropdown = fmt_picker[0];\n",
+       "\n",
+       "    for (var ind in mpl.extensions) {\n",
+       "        var fmt = mpl.extensions[ind];\n",
+       "        var option = $(\n",
+       "            '<option/>', {selected: fmt === mpl.default_extension}).html(fmt);\n",
+       "        fmt_picker.append(option);\n",
+       "    }\n",
+       "\n",
+       "    // Add hover states to the ui-buttons\n",
+       "    $( \".ui-button\" ).hover(\n",
+       "        function() { $(this).addClass(\"ui-state-hover\");},\n",
+       "        function() { $(this).removeClass(\"ui-state-hover\");}\n",
+       "    );\n",
+       "\n",
+       "    var status_bar = $('<span class=\"mpl-message\"/>');\n",
+       "    nav_element.append(status_bar);\n",
+       "    this.message = status_bar[0];\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype.request_resize = function(x_pixels, y_pixels) {\n",
+       "    // Request matplotlib to resize the figure. Matplotlib will then trigger a resize in the client,\n",
+       "    // which will in turn request a refresh of the image.\n",
+       "    this.send_message('resize', {'width': x_pixels, 'height': y_pixels});\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype.send_message = function(type, properties) {\n",
+       "    properties['type'] = type;\n",
+       "    properties['figure_id'] = this.id;\n",
+       "    this.ws.send(JSON.stringify(properties));\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype.send_draw_message = function() {\n",
+       "    if (!this.waiting) {\n",
+       "        this.waiting = true;\n",
+       "        this.ws.send(JSON.stringify({type: \"draw\", figure_id: this.id}));\n",
+       "    }\n",
+       "}\n",
+       "\n",
+       "\n",
+       "mpl.figure.prototype.handle_save = function(fig, msg) {\n",
+       "    var format_dropdown = fig.format_dropdown;\n",
+       "    var format = format_dropdown.options[format_dropdown.selectedIndex].value;\n",
+       "    fig.ondownload(fig, format);\n",
+       "}\n",
+       "\n",
+       "\n",
+       "mpl.figure.prototype.handle_resize = function(fig, msg) {\n",
+       "    var size = msg['size'];\n",
+       "    if (size[0] != fig.canvas.width || size[1] != fig.canvas.height) {\n",
+       "        fig._resize_canvas(size[0], size[1]);\n",
+       "        fig.send_message(\"refresh\", {});\n",
+       "    };\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype.handle_rubberband = function(fig, msg) {\n",
+       "    var x0 = msg['x0'] / mpl.ratio;\n",
+       "    var y0 = (fig.canvas.height - msg['y0']) / mpl.ratio;\n",
+       "    var x1 = msg['x1'] / mpl.ratio;\n",
+       "    var y1 = (fig.canvas.height - msg['y1']) / mpl.ratio;\n",
+       "    x0 = Math.floor(x0) + 0.5;\n",
+       "    y0 = Math.floor(y0) + 0.5;\n",
+       "    x1 = Math.floor(x1) + 0.5;\n",
+       "    y1 = Math.floor(y1) + 0.5;\n",
+       "    var min_x = Math.min(x0, x1);\n",
+       "    var min_y = Math.min(y0, y1);\n",
+       "    var width = Math.abs(x1 - x0);\n",
+       "    var height = Math.abs(y1 - y0);\n",
+       "\n",
+       "    fig.rubberband_context.clearRect(\n",
+       "        0, 0, fig.canvas.width / mpl.ratio, fig.canvas.height / mpl.ratio);\n",
+       "\n",
+       "    fig.rubberband_context.strokeRect(min_x, min_y, width, height);\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype.handle_figure_label = function(fig, msg) {\n",
+       "    // Updates the figure title.\n",
+       "    fig.header.textContent = msg['label'];\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype.handle_cursor = function(fig, msg) {\n",
+       "    var cursor = msg['cursor'];\n",
+       "    switch(cursor)\n",
+       "    {\n",
+       "    case 0:\n",
+       "        cursor = 'pointer';\n",
+       "        break;\n",
+       "    case 1:\n",
+       "        cursor = 'default';\n",
+       "        break;\n",
+       "    case 2:\n",
+       "        cursor = 'crosshair';\n",
+       "        break;\n",
+       "    case 3:\n",
+       "        cursor = 'move';\n",
+       "        break;\n",
+       "    }\n",
+       "    fig.rubberband_canvas.style.cursor = cursor;\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype.handle_message = function(fig, msg) {\n",
+       "    fig.message.textContent = msg['message'];\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype.handle_draw = function(fig, msg) {\n",
+       "    // Request the server to send over a new figure.\n",
+       "    fig.send_draw_message();\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype.handle_image_mode = function(fig, msg) {\n",
+       "    fig.image_mode = msg['mode'];\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype.updated_canvas_event = function() {\n",
+       "    // Called whenever the canvas gets updated.\n",
+       "    this.send_message(\"ack\", {});\n",
+       "}\n",
+       "\n",
+       "// A function to construct a web socket function for onmessage handling.\n",
+       "// Called in the figure constructor.\n",
+       "mpl.figure.prototype._make_on_message_function = function(fig) {\n",
+       "    return function socket_on_message(evt) {\n",
+       "        if (evt.data instanceof Blob) {\n",
+       "            /* FIXME: We get \"Resource interpreted as Image but\n",
+       "             * transferred with MIME type text/plain:\" errors on\n",
+       "             * Chrome.  But how to set the MIME type?  It doesn't seem\n",
+       "             * to be part of the websocket stream */\n",
+       "            evt.data.type = \"image/png\";\n",
+       "\n",
+       "            /* Free the memory for the previous frames */\n",
+       "            if (fig.imageObj.src) {\n",
+       "                (window.URL || window.webkitURL).revokeObjectURL(\n",
+       "                    fig.imageObj.src);\n",
+       "            }\n",
+       "\n",
+       "            fig.imageObj.src = (window.URL || window.webkitURL).createObjectURL(\n",
+       "                evt.data);\n",
+       "            fig.updated_canvas_event();\n",
+       "            fig.waiting = false;\n",
+       "            return;\n",
+       "        }\n",
+       "        else if (typeof evt.data === 'string' && evt.data.slice(0, 21) == \"data:image/png;base64\") {\n",
+       "            fig.imageObj.src = evt.data;\n",
+       "            fig.updated_canvas_event();\n",
+       "            fig.waiting = false;\n",
+       "            return;\n",
+       "        }\n",
+       "\n",
+       "        var msg = JSON.parse(evt.data);\n",
+       "        var msg_type = msg['type'];\n",
+       "\n",
+       "        // Call the  \"handle_{type}\" callback, which takes\n",
+       "        // the figure and JSON message as its only arguments.\n",
+       "        try {\n",
+       "            var callback = fig[\"handle_\" + msg_type];\n",
+       "        } catch (e) {\n",
+       "            console.log(\"No handler for the '\" + msg_type + \"' message type: \", msg);\n",
+       "            return;\n",
+       "        }\n",
+       "\n",
+       "        if (callback) {\n",
+       "            try {\n",
+       "                // console.log(\"Handling '\" + msg_type + \"' message: \", msg);\n",
+       "                callback(fig, msg);\n",
+       "            } catch (e) {\n",
+       "                console.log(\"Exception inside the 'handler_\" + msg_type + \"' callback:\", e, e.stack, msg);\n",
+       "            }\n",
+       "        }\n",
+       "    };\n",
+       "}\n",
+       "\n",
+       "// from http://stackoverflow.com/questions/1114465/getting-mouse-location-in-canvas\n",
+       "mpl.findpos = function(e) {\n",
+       "    //this section is from http://www.quirksmode.org/js/events_properties.html\n",
+       "    var targ;\n",
+       "    if (!e)\n",
+       "        e = window.event;\n",
+       "    if (e.target)\n",
+       "        targ = e.target;\n",
+       "    else if (e.srcElement)\n",
+       "        targ = e.srcElement;\n",
+       "    if (targ.nodeType == 3) // defeat Safari bug\n",
+       "        targ = targ.parentNode;\n",
+       "\n",
+       "    // jQuery normalizes the pageX and pageY\n",
+       "    // pageX,Y are the mouse positions relative to the document\n",
+       "    // offset() returns the position of the element relative to the document\n",
+       "    var x = e.pageX - $(targ).offset().left;\n",
+       "    var y = e.pageY - $(targ).offset().top;\n",
+       "\n",
+       "    return {\"x\": x, \"y\": y};\n",
+       "};\n",
+       "\n",
+       "/*\n",
+       " * return a copy of an object with only non-object keys\n",
+       " * we need this to avoid circular references\n",
+       " * http://stackoverflow.com/a/24161582/3208463\n",
+       " */\n",
+       "function simpleKeys (original) {\n",
+       "  return Object.keys(original).reduce(function (obj, key) {\n",
+       "    if (typeof original[key] !== 'object')\n",
+       "        obj[key] = original[key]\n",
+       "    return obj;\n",
+       "  }, {});\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype.mouse_event = function(event, name) {\n",
+       "    var canvas_pos = mpl.findpos(event)\n",
+       "\n",
+       "    if (name === 'button_press')\n",
+       "    {\n",
+       "        this.canvas.focus();\n",
+       "        this.canvas_div.focus();\n",
+       "    }\n",
+       "\n",
+       "    var x = canvas_pos.x * mpl.ratio;\n",
+       "    var y = canvas_pos.y * mpl.ratio;\n",
+       "\n",
+       "    this.send_message(name, {x: x, y: y, button: event.button,\n",
+       "                             step: event.step,\n",
+       "                             guiEvent: simpleKeys(event)});\n",
+       "\n",
+       "    /* This prevents the web browser from automatically changing to\n",
+       "     * the text insertion cursor when the button is pressed.  We want\n",
+       "     * to control all of the cursor setting manually through the\n",
+       "     * 'cursor' event from matplotlib */\n",
+       "    event.preventDefault();\n",
+       "    return false;\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype._key_event_extra = function(event, name) {\n",
+       "    // Handle any extra behaviour associated with a key event\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype.key_event = function(event, name) {\n",
+       "\n",
+       "    // Prevent repeat events\n",
+       "    if (name == 'key_press')\n",
+       "    {\n",
+       "        if (event.which === this._key)\n",
+       "            return;\n",
+       "        else\n",
+       "            this._key = event.which;\n",
+       "    }\n",
+       "    if (name == 'key_release')\n",
+       "        this._key = null;\n",
+       "\n",
+       "    var value = '';\n",
+       "    if (event.ctrlKey && event.which != 17)\n",
+       "        value += \"ctrl+\";\n",
+       "    if (event.altKey && event.which != 18)\n",
+       "        value += \"alt+\";\n",
+       "    if (event.shiftKey && event.which != 16)\n",
+       "        value += \"shift+\";\n",
+       "\n",
+       "    value += 'k';\n",
+       "    value += event.which.toString();\n",
+       "\n",
+       "    this._key_event_extra(event, name);\n",
+       "\n",
+       "    this.send_message(name, {key: value,\n",
+       "                             guiEvent: simpleKeys(event)});\n",
+       "    return false;\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype.toolbar_button_onclick = function(name) {\n",
+       "    if (name == 'download') {\n",
+       "        this.handle_save(this, null);\n",
+       "    } else {\n",
+       "        this.send_message(\"toolbar_button\", {name: name});\n",
+       "    }\n",
+       "};\n",
+       "\n",
+       "mpl.figure.prototype.toolbar_button_onmouseover = function(tooltip) {\n",
+       "    this.message.textContent = tooltip;\n",
+       "};\n",
+       "mpl.toolbar_items = [[\"Home\", \"Reset original view\", \"fa fa-home icon-home\", \"home\"], [\"Back\", \"Back to previous view\", \"fa fa-arrow-left icon-arrow-left\", \"back\"], [\"Forward\", \"Forward to next view\", \"fa fa-arrow-right icon-arrow-right\", \"forward\"], [\"\", \"\", \"\", \"\"], [\"Pan\", \"Pan axes with left mouse, zoom with right\", \"fa fa-arrows icon-move\", \"pan\"], [\"Zoom\", \"Zoom to rectangle\", \"fa fa-square-o icon-check-empty\", \"zoom\"], [\"\", \"\", \"\", \"\"], [\"Download\", \"Download plot\", \"fa fa-floppy-o icon-save\", \"download\"]];\n",
+       "\n",
+       "mpl.extensions = [\"eps\", \"jpeg\", \"pdf\", \"png\", \"ps\", \"raw\", \"svg\", \"tif\"];\n",
+       "\n",
+       "mpl.default_extension = \"png\";var comm_websocket_adapter = function(comm) {\n",
+       "    // Create a \"websocket\"-like object which calls the given IPython comm\n",
+       "    // object with the appropriate methods. Currently this is a non binary\n",
+       "    // socket, so there is still some room for performance tuning.\n",
+       "    var ws = {};\n",
+       "\n",
+       "    ws.close = function() {\n",
+       "        comm.close()\n",
+       "    };\n",
+       "    ws.send = function(m) {\n",
+       "        //console.log('sending', m);\n",
+       "        comm.send(m);\n",
+       "    };\n",
+       "    // Register the callback with on_msg.\n",
+       "    comm.on_msg(function(msg) {\n",
+       "        //console.log('receiving', msg['content']['data'], msg);\n",
+       "        // Pass the mpl event to the overridden (by mpl) onmessage function.\n",
+       "        ws.onmessage(msg['content']['data'])\n",
+       "    });\n",
+       "    return ws;\n",
+       "}\n",
+       "\n",
+       "mpl.mpl_figure_comm = function(comm, msg) {\n",
+       "    // This is the function which gets called when the mpl process\n",
+       "    // starts-up an IPython Comm through the \"matplotlib\" channel.\n",
+       "\n",
+       "    var id = msg.content.data.id;\n",
+       "    // Get hold of the div created by the display call when the Comm\n",
+       "    // socket was opened in Python.\n",
+       "    var element = $(\"#\" + id);\n",
+       "    var ws_proxy = comm_websocket_adapter(comm)\n",
+       "\n",
+       "    function ondownload(figure, format) {\n",
+       "        window.open(figure.imageObj.src);\n",
+       "    }\n",
+       "\n",
+       "    var fig = new mpl.figure(id, ws_proxy,\n",
+       "                           ondownload,\n",
+       "                           element.get(0));\n",
+       "\n",
+       "    // Call onopen now - mpl needs it, as it is assuming we've passed it a real\n",
+       "    // web socket which is closed, not our websocket->open comm proxy.\n",
+       "    ws_proxy.onopen();\n",
+       "\n",
+       "    fig.parent_element = element.get(0);\n",
+       "    fig.cell_info = mpl.find_output_cell(\"<div id='\" + id + \"'></div>\");\n",
+       "    if (!fig.cell_info) {\n",
+       "        console.error(\"Failed to find cell for figure\", id, fig);\n",
+       "        return;\n",
+       "    }\n",
+       "\n",
+       "    var output_index = fig.cell_info[2]\n",
+       "    var cell = fig.cell_info[0];\n",
+       "\n",
+       "};\n",
+       "\n",
+       "mpl.figure.prototype.handle_close = function(fig, msg) {\n",
+       "    var width = fig.canvas.width/mpl.ratio\n",
+       "    fig.root.unbind('remove')\n",
+       "\n",
+       "    // Update the output cell to use the data from the current canvas.\n",
+       "    fig.push_to_output();\n",
+       "    var dataURL = fig.canvas.toDataURL();\n",
+       "    // Re-enable the keyboard manager in IPython - without this line, in FF,\n",
+       "    // the notebook keyboard shortcuts fail.\n",
+       "    IPython.keyboard_manager.enable()\n",
+       "    $(fig.parent_element).html('<img src=\"' + dataURL + '\" width=\"' + width + '\">');\n",
+       "    fig.close_ws(fig, msg);\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype.close_ws = function(fig, msg){\n",
+       "    fig.send_message('closing', msg);\n",
+       "    // fig.ws.close()\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype.push_to_output = function(remove_interactive) {\n",
+       "    // Turn the data on the canvas into data in the output cell.\n",
+       "    var width = this.canvas.width/mpl.ratio\n",
+       "    var dataURL = this.canvas.toDataURL();\n",
+       "    this.cell_info[1]['text/html'] = '<img src=\"' + dataURL + '\" width=\"' + width + '\">';\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype.updated_canvas_event = function() {\n",
+       "    // Tell IPython that the notebook contents must change.\n",
+       "    IPython.notebook.set_dirty(true);\n",
+       "    this.send_message(\"ack\", {});\n",
+       "    var fig = this;\n",
+       "    // Wait a second, then push the new image to the DOM so\n",
+       "    // that it is saved nicely (might be nice to debounce this).\n",
+       "    setTimeout(function () { fig.push_to_output() }, 1000);\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype._init_toolbar = function() {\n",
+       "    var fig = this;\n",
+       "\n",
+       "    var nav_element = $('<div/>');\n",
+       "    nav_element.attr('style', 'width: 100%');\n",
+       "    this.root.append(nav_element);\n",
+       "\n",
+       "    // Define a callback function for later on.\n",
+       "    function toolbar_event(event) {\n",
+       "        return fig.toolbar_button_onclick(event['data']);\n",
+       "    }\n",
+       "    function toolbar_mouse_event(event) {\n",
+       "        return fig.toolbar_button_onmouseover(event['data']);\n",
+       "    }\n",
+       "\n",
+       "    for(var toolbar_ind in mpl.toolbar_items){\n",
+       "        var name = mpl.toolbar_items[toolbar_ind][0];\n",
+       "        var tooltip = mpl.toolbar_items[toolbar_ind][1];\n",
+       "        var image = mpl.toolbar_items[toolbar_ind][2];\n",
+       "        var method_name = mpl.toolbar_items[toolbar_ind][3];\n",
+       "\n",
+       "        if (!name) { continue; };\n",
+       "\n",
+       "        var button = $('<button class=\"btn btn-default\" href=\"#\" title=\"' + name + '\"><i class=\"fa ' + image + ' fa-lg\"></i></button>');\n",
+       "        button.click(method_name, toolbar_event);\n",
+       "        button.mouseover(tooltip, toolbar_mouse_event);\n",
+       "        nav_element.append(button);\n",
+       "    }\n",
+       "\n",
+       "    // Add the status bar.\n",
+       "    var status_bar = $('<span class=\"mpl-message\" style=\"text-align:right; float: right;\"/>');\n",
+       "    nav_element.append(status_bar);\n",
+       "    this.message = status_bar[0];\n",
+       "\n",
+       "    // Add the close button to the window.\n",
+       "    var buttongrp = $('<div class=\"btn-group inline pull-right\"></div>');\n",
+       "    var button = $('<button class=\"btn btn-mini btn-primary\" href=\"#\" title=\"Stop Interaction\"><i class=\"fa fa-power-off icon-remove icon-large\"></i></button>');\n",
+       "    button.click(function (evt) { fig.handle_close(fig, {}); } );\n",
+       "    button.mouseover('Stop Interaction', toolbar_mouse_event);\n",
+       "    buttongrp.append(button);\n",
+       "    var titlebar = this.root.find($('.ui-dialog-titlebar'));\n",
+       "    titlebar.prepend(buttongrp);\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype._root_extra_style = function(el){\n",
+       "    var fig = this\n",
+       "    el.on(\"remove\", function(){\n",
+       "\tfig.close_ws(fig, {});\n",
+       "    });\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype._canvas_extra_style = function(el){\n",
+       "    // this is important to make the div 'focusable\n",
+       "    el.attr('tabindex', 0)\n",
+       "    // reach out to IPython and tell the keyboard manager to turn it's self\n",
+       "    // off when our div gets focus\n",
+       "\n",
+       "    // location in version 3\n",
+       "    if (IPython.notebook.keyboard_manager) {\n",
+       "        IPython.notebook.keyboard_manager.register_events(el);\n",
+       "    }\n",
+       "    else {\n",
+       "        // location in version 2\n",
+       "        IPython.keyboard_manager.register_events(el);\n",
+       "    }\n",
+       "\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype._key_event_extra = function(event, name) {\n",
+       "    var manager = IPython.notebook.keyboard_manager;\n",
+       "    if (!manager)\n",
+       "        manager = IPython.keyboard_manager;\n",
+       "\n",
+       "    // Check for shift+enter\n",
+       "    if (event.shiftKey && event.which == 13) {\n",
+       "        this.canvas_div.blur();\n",
+       "        // select the cell after this one\n",
+       "        var index = IPython.notebook.find_cell_index(this.cell_info[0]);\n",
+       "        IPython.notebook.select(index + 1);\n",
+       "    }\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype.handle_save = function(fig, msg) {\n",
+       "    fig.ondownload(fig, null);\n",
+       "}\n",
+       "\n",
+       "\n",
+       "mpl.find_output_cell = function(html_output) {\n",
+       "    // Return the cell and output element which can be found *uniquely* in the notebook.\n",
+       "    // Note - this is a bit hacky, but it is done because the \"notebook_saving.Notebook\"\n",
+       "    // IPython event is triggered only after the cells have been serialised, which for\n",
+       "    // our purposes (turning an active figure into a static one), is too late.\n",
+       "    var cells = IPython.notebook.get_cells();\n",
+       "    var ncells = cells.length;\n",
+       "    for (var i=0; i<ncells; i++) {\n",
+       "        var cell = cells[i];\n",
+       "        if (cell.cell_type === 'code'){\n",
+       "            for (var j=0; j<cell.output_area.outputs.length; j++) {\n",
+       "                var data = cell.output_area.outputs[j];\n",
+       "                if (data.data) {\n",
+       "                    // IPython >= 3 moved mimebundle to data attribute of output\n",
+       "                    data = data.data;\n",
+       "                }\n",
+       "                if (data['text/html'] == html_output) {\n",
+       "                    return [cell, data, j];\n",
+       "                }\n",
+       "            }\n",
+       "        }\n",
+       "    }\n",
+       "}\n",
+       "\n",
+       "// Register the function which deals with the matplotlib target/channel.\n",
+       "// The kernel may be null if the page has been refreshed.\n",
+       "if (IPython.notebook.kernel != null) {\n",
+       "    IPython.notebook.kernel.comm_manager.register_target('matplotlib', mpl.mpl_figure_comm);\n",
+       "}\n"
+      ],
       "text/plain": [
-       "<matplotlib.figure.Figure at 0x7f4b46390668>"
+       "<IPython.core.display.Javascript object>"
       ]
      },
      "metadata": {},
@@ -75,57 +1725,904 @@
     },
     {
      "data": {
+      "text/html": [
+       "<img src=\"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAdsAAAE8CAYAAACIOxsOAAAgAElEQVR4Xu2dD7glZ13fvycXd7MJd4EEQswmIZAGMclSC0bXiNaARSSBJxCqBgNLgkLQEkylovwRVNKKQKVpKQEpGEBjW9EI2ioqxJLWtbaUNstWEqyb4IYEyCJ3Mcuul7t9vnvfSQ5n55x73jnv75w5Zz7zPPeBzbzzzszn/c185veb95zTEwsEIAABCEAAAqEEeqG90zkEIAABCEAAAkK2BAEEIAABCEAgmACyDQZM9xCAAAQgAAFkSwxAAAIQgAAEggkg22DAdA8BCEAAAhBAtsQABCAAAQhAIJgAsg0GTPcQgAAEIAABZEsMQAACEIAABIIJINtgwHQPAQhAAAIQQLbEAAQgAAEIQCCYALINBkz3EIAABCAAAWRLDEAAAhCAAASCCSDbYMB0DwEIQAACEEC2xAAEIAABCEAgmACyDQZM9xCAAAQgAAFkSwxAAAIQgAAEggkg22DAdA8BCEAAAhBAtsQABCAAAQhAIJgAsg0GTPcQgAAEIAABZEsMLCKBvZLeIOlXF/HkOKejBKqxfRE8IDAPBJDtPIxS+4/xFkn/MB3mVyXdLenjkn5Z0v/OOHwL8nvSX8ZmxzSdVLZnSforSY+V5L7mZXmppBdIepKkz0vyefQvx0n6Z5JeIuk0SZ+Q9GMDY3ShpF+S9C2S9kt6i6Tr+zp5taTLJH2TpL+R9FuS/N++MmVIG8n2uyX9tKQLJD1yyFheIukX0rn8taTXSfr3fedxeorhiySZndddK8kx7sWif+/Aef+OpEunzILdzQEBZDsHgzQHh2jZ/o90Y94k6WxJPyLpH6cb84fHPAdkOyaoIc1eIcn8T03cB2V7jaTXJEncIamSs8X5ZUmWy6ckvU3S+yWdmzJIC/k30j5/T9KvSfqfkh4t6V2S/lTSlZMdevbWG8n2+yV9h6TPpmMcfHDyA8kuST8h6Q8kfaekd0p6mqT/luT639ODo6W9LOkdkv5M0sv6ZGtZW+jVYhH7IYQFAl9HANkSECUIWLb+syz7l3dLekbKKv5O0rPTzd438QOSPiTplSkrqssSfIM8KOntkp4i6URJn5RkqTgrG7Y4G/U2zkj85387m3O2XS3fK+kXJZ2XbsjOwn0z9XJkoOOfS8f7dEnfl9a5vW/UFo6zyCdLulXSwyQdlvQ4Sf8q7X9F0gcl/VQ6H3dxgqQ3S/oBSd+Qju3HJd2V+rdMfH3uSzf3+yW9aSDLHHb+ZumxGJStJfK7kv552tD9O6N7Yzp3M/J4PL6vY7e1gL59yM78QGVJnTRiPCxvj/03poqBBXVTX3uPj1l5jC1J/9tC+5O+Ni+X9FpJx0v6dylbXUsPDqNieFiVwuf1rZI8ptXy65KWJP1gYvBpSduScN3Gbc3vUenhZBjnUcfDuo4SQLYdHfjCpz1Mtn8/yXFHyggsFsvT2ZNvYjdI+s/pBr8llfRcxnxuOr4vSDoj3ag/JulQEtyzJJ3TV84bPB3frLemMuJ/kWSJ/bCkxyRpOpNz1mJZ+oZu+b8nlVR/M2UqXv9tScQukVrKfyTp4ZK+JunP00OEBeVSqh8AnpNK4M4ufY43S/IDh4/l36Rt/kk62Pelm7ZLl37weFUq/1ra7t+yNYd/m47tH0n615K2p75HDeEwCfgB5QOS/mXfxv8vPSS8UJIzX1cknti3/mdTedUPOn6IGFwsRZem/XAxbPE5Onv8oqTqPCy6/5M28Hh5/N2Ps2SXpd3OD1t+SPMrij+UZHYeTx+ny+V+gNnone0w2ZrB30uxVR23Y8APYWcmBn4FcnIqp7uNS9OOFz98fDTt2w8aPi9ntI4PZ8FfKnx90d0CEEC2CzCILTiFYbJ1FmK5/tDAu7DqkJ+X3g9WN+pxysh+d+Yynd+3+cZbt/jm7SzzirTS2Yql8i+S4H1TdTbqG2O1+Abvm6kz8bobtLNP79dt/iJt74zNGarf4/1HSc6EnH1ZXC69+iGjWvwQ4ZuxpWXpu49TJDnr9VL17xu5y5uWrUud/eJz6fet6RxGDfsw2TqT9vvEZyYeflCwzD+SMnbvy0K+KpWK/VDy++mBx+94PzewU2fxbm+e12XEobNDP8z8fJ9s/arB2asXP4g54/7mxMnvSi3d/vH0O/VKeKN2PUy2F6eHJMfRH6fM3Q8Em9OfH5gcMz6uf5rG2Zmvs9vnp8zc4+v+/WDlErzjy+Pph4PB6kgGHpouIgFku4ijOv1zGiZb37j8xO+y3H9IpTnfkPyO6xGSHpL+LBovdbJ1G2dXFrNv+P63BefMxu8O6xbL1plL/8QeZ593JjE6K3WGuNq3sft1Cdcl1GE3aGfXngBzWyq9Ohv1fpyleVKY31v6hm2pOWt2Jl4tvtZ83BbJP0g3cZeG+xevt1B8U7dszcUZebU4q7IYKkkNG+lhsn1omtDjjNky8AOJH4a8uHzr5UfT8butMzRn4JaN3wPf27dDy8jvb6tt+1kOHpflZB5+qPIDmOPC5+d9efF4Wda/kv5tVi4RW1p+oPKrAwt9cDwttqaZrXfluPqZxNnvdj12O1OW7fV+QLoxHbdl79KzXykMe3h0RuxzcUXEcxhYIPAAAWRLMJQgsFEZ2e/7nMk483Pp0DdNZ5a+mfkmWsVhnWx9M3SW6DLt7UneLjW65Djsoz2+4TkDdNm1Wvpl+38luVzsG2n/4huqhTxMtn7X6JKzz8Eycnn0vnQezvD8AGEB+N2vJ4n5uAcXH5tn83pikUvGg4uF5rJy3QSgYZwH+9joXaJLts5K70mTgVwurSb9uC+Ph9+vuozvBwg/ULi9y9te/GDiTN4PP87ER81E9hj7IcXvgi13n5srAi69VqKsmz3uhwG/b/c5+/hcju8fz9/ue286KoY3mlnuqocfJMzC4+tXAc6o+xe/o/VDic/XMVy9Fqnbr+Pary3MhwUCDxBAtgRDCQIbTZDyDc/vOn3zdka5O+3UInW2UMWhZ8q6jPtdfQflkqNLrr5Ze/ENzyVGlzpHydaToZz9enHp2SVBT4jye2LPrPV/8zvkuqUqY/qd3l/2NfDkKGd6lq1v/J7gZYFYFpawy75erk4Twbx9f3ZbdfUESRa+3wPvGXIMkbKtdulMzOdnwZhz3eJsz4J0ZaFi6fe+Hkdnnv540KjF4+ZM0Nm/F4+1z9mzeseVrasi5liNZ/VawBKfJLPtP2736YrFf+qLtcHzcrz6fbG5+cFscLG0XeGwjP1wyQKBBwggW4KhBIFxPvrjm5ll6xJpVXr1Zzg9AaqKQ5dQvc7vRZ35+EbuWb9+T1aVU/0ZUGdLzhpHydaTkjz719J1W79H9btSZ57OKP9rmt1bidfZt0ujlrEzNwvG5VNnwH8rySVff/zDpVWXOJ35Oat1WdxZtt99+kbsxe0s5P+VHib8rtcy9mxbH5MXT+7xO1Hvw+9ifQO3/P3A4fZNZOubvf8889fC97tZLy7DevExWPDOwl3SdYXBDy7mW71jNCuXbj1e/v8uAbss6ocDL84w3a//nA1Wi8e2ynz7Y8oi95hb1p9Jk5xcqnWlYVzZOsO19H08Hk+/23UfoyZIufLghx0/nLnc7Xe0FqGPwZm4S/Rm5IlXHi+/s7ck/d7a5+LF4+HysmPR5+vY8wQyfyzKi8fdbN2n9+PXB+7X/dSx6OfC/+8YAWTbsQEPOt3+L7VwBtL/pRbVjd679ntB3+A9mcSlYJdxq4+4eL3fWXoyjL/YwjdLz0a16NzOWZRLdL4pWnCjviGq+uiPZ5Z6O//bN1YfZ7VY6M6qLV6XCJ3VuF9PCPLiUqAnO/mjPX5HWn2syaLye0eLy4tvwr6ZV++lq/4tdj9MeFatb8C+Iftc/fDgxX34PaWzPn/pgj/iY6FUX5rQRLY+xtfXjHF1nVskfs/tmdx+6HDG6Eli/WVglz99zH6v6olaFkr/+8dhE3+GfQGI9+2P9Tgr9UOKZe0HAvczrmx9Sn6N4AcRc/MXSbhkP+qjP44hZ76DS1Wa9gOV1/vLO7z4//vBx+NULX4w8nj4Y01+IHK8+MGhWvyRKT/EOUYs6Go2cv9DSM0h8J+6SADZdnHUOWcIQAACEJgqAWQ7VdzsDAIQgAAEukgA2XZx1DlnCEAAAhCYKgFkO1Xc7AwCEIAABLpIANl2cdQ5ZwhAAAIQmCoBZDtV3OwMAhCAAAS6SKBTst28efORr361+inKLg435wwBCEBgfgn0er25ddbcHniTcFlaWjqyujrqK1yb9Mo2EIAABCAwDQLIdhqUC+wD2RaASBcQgAAEZkQA2c4IfO5ukW0uMdpDAAIQaA8BZNuesRh5JMh2TgaKw4QABCBQQwDZzklYINs5GSgOEwIQgACynd8YQLbzO3YcOQQgAAEy2zmJAWQ7JwPFYUIAAhAgsx0ZA/4pM3+Qtfow66vSz4b59x3fKWlL+rkz/26pfy7NS8S62oNEtly/EIAABOaXAJntg2Nn2V4iaXffcPqzvP4tSP925a3pN0L9w9VXpR8NL71uaCQh2/m9yDhyCEAAAsh2tGwvSD/2fH5q5h/KtpT94+AR65At1yQEIACBBSSAbL9etl9OGauz2FdLelrKYi/uG/v7JZ0u6aKAdfuHxRiZ7QJefZwSBCDQGQLI9sGhPkPSZyVtlvQ2ScuSfjtAqKMk3S/bayX57+jS6/XOWFtb60xgcqIQgEB3CWzfvv3oyd92220LAwHZ1g+lR/pDkn6gpox8p6QTh5SRJ11HGXlhLi1OBAIQaEoA2TYlF7NdyR8isDwfIqkqI79R0rmSLksTpHb2TZA6W9KVko4LWIdsY2KFXiEAgTkigGzbNVglZesZxh+UtJT+9ki6RtLnJF2YPvpzfN9Hf+5NKCLW1VLmnW27go+jgQAE4ggg2zi2TXouKdsm+5/qNsh2qrjZGQQgMEMCyHaG8Gt2jWzbNR4cDQQgAIEiBJBtEYzFOkG2xVDSEQQgAIH2EEC27RkLHwmybdd4cDQQgAAEihBAtkUwFusE2RZDSUcQgAAE2kMA2bZnLMhs2zUWHA0EIDDnBA4dOqTDhw+34ix27PBvvEi7du1qxfFs2rRJmzf7+46aL3ypRXN2U92S2chTxc3OINApAhbtaaedpf377+nUeY97sieddKruvnvvRMJFtuPSnnE7ZDvjAWD3EFhgAgcOHNDWrVsl7UvfVDvrk13PbKU2ZLYHJG3TysqKlpf9Lb7NFmTbjNvUt0K2U0fODiHQGQIPynalJbJd/25kqQ3fjWzZbkW2XbkakG1XRprzhMD0CSDbUcyRLbORp39NskcIQGABCSBbZDuKALJdwIueU4IABKZPANkiW2SbCFBGnv4NiD1CoCsEkC2yRbbItiv3O84TAjMjgGyRLbJFtjO7AbFjCHSFALJFtsgW2Xblfsd5QmBmBJAtskW2yHZmNyB2DIGuEGifbNtEno/+MBu5TfHIsUAAAnNLANmS2ZLZktnO7Q1sVge+iL+YMiuWXdkvskW2yBbZduV+V+w8kW0xlJ3pCNkiW2SLbDtzwyt1osi2FMnu9INskS2yRbbdueMVOlNkWwhkh7pBtsgW2SLbDt3yypwqsi3DsUu9IFtki2yRbZfueUXOFdkWwTiyk0VjjGyRLbJFtvF3zgXbw6KJoI3Ds2iMkS2yRbbIto332lYf06KJoI2wF40xskW2yBbZtvFe2+pjWjQRtBH2ojFGtsgW2SLbNt5rW31MiyaCNsJeNMbIFtkiW2Tbxnttq49p0UTQRtiLxhjZIltki2zbeK9t9TEtmgjaCHvRGCNbZItskW0b77WtPqZFE0EbYS8aY2SLbJEtsm3jvbbVx7RoImgj7EVjjGyRLbJFtm2819Ye06FDh3T48OGZH++OHTuOHsOuXbtmfizVAWzatEmbN29uzfFMeiDIdlKC87Q9v2fL79nOU7wu+LFatI95zFm69957FvxMm53eox99qu68c+/CCBfZNouD+dxqKrJ9vaQ3SNouabekX5N0kaRvlLQs6St97Pw0/U5JWyTtlXSFpM+n9U3XjRwaZDufkbuQR12V4W66aZ+2bPG1MbvlmmvWM9vrr29HZnvw4AFdfvk2raysaHl5tmxKjQqyLUVyHvoJl+2TJF0n6VxJFyfZPjX9770DsrX37pD0Ikm3SnqtpMdJukpS03UbDgKy3RARDaZFoJLtzTev6IQTZiuUl7zED8fSu95127ROf+R+7r//gC69dCuybcVo1B8E72xHDU6obP1u5RZJz5f0MUmXJMlWB3RkQLYXSHqvpPNTg0em7Pahkpqu2zAyOyVbSUfOP7/iuyEbGkyZwNramvbs2aMzzzxXxx133JT3/vW727fPD77Stm3nzPQ4qp2bzV137dG5586eTSkgd9yxzvicc9rBeNLzquJ3PbmabfxOei7lt1+TNHn87t69u85Zb5J0l6S3J2luJNvLUhbrDLha7pd0eio7O8PNXbd/I2YRsh2smzetfzfdbtQ5I9uNImKG65HtcPjIdoaBOeauke0oUGGy/Y5UPn6ak6kuyXawbv6pgNr4qJr6yMtiaWnpyOrq6piXDs2mTYAy8nDilJGnHY35+6OMPIpZWBn5pyVdI6n6CIOzU7+jvVLSR9IRjVNGvlPSiUPKyOOs2zBgSma2dXVzz/QqXRsfVVNHthsOeXsbIFtk297o3PjIkO1MZDu4U88s3qiM7Bq/32Hs7JsgdXYSdNN1GwZISdnW1c2/KaA27qncw2rqI+vmZLYbxsNMGyBbZDvTAJxw58i2dbL9LUnf5qkXku5Ok6a+Lx3lhemjP8f3ffTHGbGXputGRlAp2Q6rm89attdK8t/RpdfrneH3KiztJIBskW07I3O8o0K2rZDteIM1g1alZDusbn69pBcMTLEep/5dVyoeZ7uRCBcps120zyh64JBtvGzb8g1dPtO2fUvXpN/QhWyR7SgCpWQ7uI+qbr4noDY+qqaObGfwxFZql8g2VrYW7Vlnnql7Pl99UU6pkVuMfk495RTtveuuxt/QhWyR7Sxl66/Malr/brrd0PMls233TRHZxsq24rvvhhu0vMVzF2e77HjNa44ewK7r/MU/s10OHDyobVdfPdGXhiBbZDsL2c72yhmyd2TbymF54KCQ7XRku3Ljja2Q7faf/MmjJ3zbW98688C0bLfu3Ilsw0Yi7KM/YUdcuuOoMnLp4yzSH7ItgjGsE2SLbMOCa4OOkW00eWSLbKNjLKh/JkgFgW1ptyW+1KJ6mCGzPXaQkW104CNbZBsdY0H9I9sgsC3tFtnGDgyyjeUrIVtkGx1jQf0j2yCwLe0W2cYODLKN5Yts13+7rzML72zbPdRtemfbNlLINnZEkG0sX2SLbKMjLKx/MtswtK3sGNnGDguyjeWLbJFtdISF9Y9sw9C2smNkGzssyDaWL7JFttERFtY/sg1D28qOkW3ssCDbWL7IFtlGR1hY/8g2DG0rO0a2scOCbGP5IltkGx1hYf0j2zC0rewY2cYOC7KN5YtskW10hIX1j2zD0LayY2QbOyzINpYvskW20REW1j+yDUPbyo6RbeywINtYvsgW2UZHWFj/yDYMbSs7Rraxw4JsY/kiW2QbHWFh/SPbMLSt7BjZxg4Lso3li2yRbXSEhfWPbMPQtrJjZBs7LMg2li+yRbbRERbWP7INQ9vKjpFt7LAg21i+yBbZRkdYWP/INgxtKztGtrHDgmxj+SJbZBsdYWH9I9swtK3sGNnGDguyjeWLbJFtdISF9Y9sw9C2suNFlG2bQCPb6NHg92z5ib3oGAvqH9kGgW1pt8g2dmCQbSxfMlsy2+gIC+sf2YahbWXHyDZ2WJBtLF9ki2yjIyysf2QbhraVHSPb2GFBtrF8kS2yjY6wsP6RbRjaVnaMbGOHBdnG8kW2yDY6wsL6R7ZhaFvZMbKNHRZkG8sX2SLb6AgL6x/ZhqFtZcfINnZYkG0sX2SLbLMj7NChQzp8+HD2dqU32LFjx9Eud+3aVbrrxv1t2rRJmzdvbrz9gQMHtHXrVt1884pOOGG5cT+LuCGyjR1VZBvLF9ki26wIs2jPOu003bN/f9Z2XWl86kknae/ddzcWLrIdHinINvYqQraxfJEtss2KsEoG+yTNOu9az2ultuS1/sj6NkkrKytaXm5GB9ki26wLsmBjZFsQZm1XfKkFX2qREWOVDFZaINvt6bhvyzj+yKbrlxKyjWJMZhtFdr1fZBvLl8yWzDYrwpDtcFzINiuUshsj22xkWRsg2yxcDRqT2ZLZZoQNskW2GeFStCmyLYrzmM6QbSxfMlsy26wIQ7bINitgCjZGtgVh1nSFbGP5IltkmxVhyBbZZgVMwcbItiBMZBsLs7Z3ysiUkTPCDtki24xwKdoU2RbFSRk5Fmfd48zRKZSTfFrBnfZ6vbl11tweeJNYWVpaOrK6utpk06PbIFtk2zh4JtwQ2U4IcIPNKSPH8qWMTBk5K8KQLbLNCpiCjZFtQZiUkWNhUkauJUBmmxF2yBbZZoRL0abItihOysixOCkj1xAoLdubJT1W0pqkr0h6uaRPSnq8pBslnSzpPkkvlHRHOp6IdbWhRBk57grjc7ZxbN0zso3lSxk5li9l5PJl5IdJ+nIatkslvU7SkyV9VNJ7JH1A0hWSrpL01NQuYh2yjb52BvpHtrHAkW0sX2QbyxfZlpdt/4i9QNIrJD1T0u0pq/2apKWU3Z7jyWUB674wLGzIbOMuKGQbx5bMNpate0e20Yz56E/pMrJH7N2Snp5E+gxJx0t6n6Tz+oZzT8pwvf/S6z6BbKMvnGP7R7axzMlsY/ki21i+ZLbxme3lqZRcWqijJN0v22sl+e/o0uv1zlhb8+vkZgsTpIZzQ7bNYmrcrZDtuKSatUO2zbiNvxWZbURm28//oKSzJH2aMvL4YTlOS371ZxxKi9MG2caOJbKN5UtmWzazPTEJ9a40bM+W9A5Jp0v6WCovVxOkXizpotTuloB1tZHDO9u4C4rMNo6te0a2sXyRbSxfZFtWto+S9GFJlq4nQu2X9EpJLus+IX305xGSvpQ++uNs10vEOmQbfe0M9I9sY4Ej21i+yDaWL7ItK9vo0Zq4fzLbiREO7QDZxrEls41l696RbTRj3tlGv7ONHsGs/pFtFq6sxsg2C1d2YzLbbGRZGyDbLFwNGiNbZJsRNsxGHg4L2WYEUoOmyLYBtIxNkG0GrEZNkS2yzQgcZItsM8KlaFNkWxTnMZ0h21i+vLPlnW1WhCFbZJsVMAUbI9uCMGu6QraxfJEtss2KMGSLbLMCpmBjZFsQJrKNhVnbO2VkysgZYYdskW1GuBRtimyL4qSMHIuz7nFG0latrKxoeXm58d57vd7cOmtuD7zJaDEbuQm18bZhgtR4nJq2QrZNyY23HWXk8Tg1b0Vmi2wzoofMlsw2I1yKNkW2RXGS2cbiJLOtIYBsM4IO2SLbjHAp2hTZFsWJbGNxIltku3RkdXW1cZi1SbaNTyJoQ8rIQWBTt8g2li9l5Fi+zEZmNnJWhCFbMtusgCnYGNkWhFmXdx08qK07d040gae6P0grkppPAoo901n1zjtbysgZsYdskW1GuBRtimyL4qSMHIuTMjJlZMrIUdcYZeQosuv9IttYvpSRY/lSRqaMnBVhZLZktlkBU7Axsi0IkzJyLMza3ikjU0bOCDtki2wzwqVoU2RbFCdl5FiclJEpI1NGjrrGKCNHkaWMHEt2vXfKyNGUyWzJbDNijMyWzDYjXIo2JbMtipPMNhYnmS2ZLZlt1DVGZhtFlsw2liyZ7TT4MkGKCVJZcUZmS2abFTAFG5PZFoRZl3fxOdtYwKKMTBk5I8SQLbLNCJeiTZFtUZyUkWNxUkamjEwZOeoao4wcRZYycixZysjT4EsZmTJyVpyR2ZLZZgVMwcZktgVhUkaOhVnbO2VkysgZYYdskW1GuBRtimyL4qSMHIuTMjJlZMrIUdcYZeQospSRY8lSRp4GX8rIlJGz4ozMlsw2K2AKNiazLQiTMnIsTMrItQQoI2eEHbJFthnhUrQpsi2KkzJyLE7KyJSRKSNHXWOUkaPIUkaOJUsZeRp8KSNTRs6KMzJbMtusgCnYmMy2IEzKyLEwKSNTRl5aIrONusrIbKPIktnGkiWznQZfMlsy26w4I7Mls80KmIKNyWwLwiSzjYVJZktmS2Ybd42R2caxdc/INpYvP7EXy5fMlsw2K8LIbMlsswKmYGNkWxAmmW0sTDJbMlsy27hrjMw2ji2ZbSxb905mG82Yr2vkc7YZMUZmS2abES5Fm5LZFsV5TGfINpYvZWTKyFkRhmyRbVbAFGyMbAvCpIwcC5MyMmVkyshx1xhl5Di2lJFj2VJGjudLZls2sz1Z0vslnS3pkKTbJb1U0n2Sdkh6p6QtkvZKukLS59MQR6yrjR5kG3dRIds4tsg2li2yjeeLbMvK9iRJT5R0Sxq6N0t6uKSXSLpD0osk3SrptZIeJ+kqre+/9LqhkYNs4y4qZBvHFtnGskW28XyRbVnZDo7YZZJeJulnJL1X0vmpwSNTdvtQSRcErEO207h2BvaBbGOh8842li8TpGL5Its42R4n6SOSPiRpX8piL+4bzvslnS7pooB1+4eFDZlt3AWFbOPYktnGsiWzjeeLbONk+3ZJ2yQ9V9JzAoQ6StL9sr1Wkv+OLr1e74y1tbXGkcVs5OHokG3jsBprQzLbsTA1bkRm2xjdmBvyOduIz9m+Jb27fVaaKFVXKr5T0olDysiTrqOMPGb4l2yGbEvSPLYvZBvLF9nG8iWzLZ/ZXifpQkkuGbtU7MUlZU+C2tk3Qcozlq8MWodso6+bmv6RbSx0ZBvLF9nG8kW2ZWV7nqTd6SM/B9PQfUbS85KA/dGf4/s++nNvamM5l15XGzm8s427oJBtHFv3jGxj+SLbWL7Itqxso0dr4v6R7cQIh3aAbOPYIttYtu4d2UYz5p1txCi+Kv8AAA8ASURBVDvb6FFr3D+ybYxuww2R7YaIJmpAZjsRvg03RrYbIpqwAbJFthkhxGzk4bCQbUYgNWiKbBtAy9gE2WbAatQU2SLbjMBBtsg2I1yKNkW2RXEe0xmyjeXLO1ve2WZFGLJFtlkBU7Axsi0Is6YrZBvLF9ki26wIQ7bINitgCjZGtgVhIttYmLW9U0amjJwRdsgW2WaES9GmyLYoTsrIsTjrHmckbdXKyoqWl5cb773X682ts+b2wJuMFrORm1AbbxsmSI3HqWkrZNuU3HjbUUYej1PzVmS2yDYjeshsyWwzwqVoU2RbFCeZbSxOMtsaAsg2I+iQLbLNCJeiTZFtUZzINhYnskW2S0dWV1cbhxmyRbaNg2fCDZHthAA32JwycixfZiMzGzkrwpAtss0KmIKNkW1BmHV518GD2rpz50QTeKr7g7QiqfkkoNgznVXvvLOljJwRe8gW2WaES9GmyLYoTsrIsTgpI1NGpowcdY0xGzmK7Hq/yDaWL2XkWL6UkSkjZ0UYmS2ZbVbAFGyMbAvCpIwcC7O2d8rIlJEzwg7ZItuMcCnaFNkWxUkZORYnZWTKyJSRo64xyshRZCkjx5Jd750ycjRlMlsy24wYI7Mls80Il6JNyWyL4iSzjcVJZktmS2YbdY2R2UaRJbONJUtmOw2+TJBiglRWnJHZktlmBUzBxmS2BWHW5V18zjYWsCgjU0bOCDFki2wzwqVoU2RbFCdl5FiclJEpI1NGjrrGKCNHkaWMHEuWMvI0+FJGpoycFWdktmS2WQFTsDGZbUGYlJFjYdb2ThmZMnJG2CFbZJsRLkWbItuiOCkjx+KkjEwZmTJy1DVGGTmKLGXkWLKUkafBlzIyZeSsOCOzJbPNCpiCjclsC8KkjBwLkzJyLQHKyBlhh2yRbUa4FG2KbIvipIwci5MyMmVkyshR1xhl5CiylJFjyVJGngZfysiUkbPijMyWzDYrYAo2JrMtCJMycixMysiUkZeWyGyjrjIy2yiyZLaxZMlsp8GXzJbMNivOyGzJbLMCpmBjMtuCMMlsY2GS2ZLZktnGXWNktnFs3TOyjeXLT+zF8iWzJbPNijAyWzLbrIAp2BjZFoRJZhsLk8yWzJbMNu4aI7ONY0tmG8vWvZPZRjPm6xr5nG1GjJHZktlmhEvRpmS2RXEe0xmyjeVLGZkyclaEIVtkmxUwBRsj24IwKSPHwqSMTBmZMnLcNUYZOY4tZeRYtpSR4/mS2ZLZZkUZmS2ZbVbAFGxMZlsQJpltLEwy2/DM9i2SLpN0lqTtknanPT5e0o2STpZ0n6QXSrojcN3QQCKzjbvGyGzj2JLZxrIls43nS2ZbNrN9iqQ7JX1c0iV9sv2opPdI+oCkKyRdJempaXgj1iHbaVw7A/tAtrHQyWxj+TJBKpYvsi0r22q09vbJ9hRJt6es9muSllJ2e47W91163RdGhQyZbdwFhWzj2JLZxrIls43ni2zjZftkSe+TdF7fcO5JGa5lW3rdJ5DtNC6cY/eBbGO5k9nG8iWzjeWLbBdfttdK8t/RpdfrnbG2ttY4qpggNRwdsm0cVmNtiGzHwtS4EbJtjG7MDflSi4gvtaCMPGb4LVIzZBs7msg2li+yjeVLZhuf2XoEb5H07r4JUi+WdFEa2oh1Q6OGd7ZxFxSyjWPrnpFtLF9kG8sX2ZaV7fWSnivpVElfTBOh/K72CemjP4+Q9KX00Z9Pp6GNWIdso6+bmv6RbSx0ZBvLF9nG8kW2ZWUbPVoT909mOzHCoR0g2zi2ZLaxbN07so1mzDvbiHe20aPWuH9k2xjdhhsi2w0RTdSAzHYifBtujGw3RDRhA2SLbDNCiNnIw2Eh24xAatAU2TaAlrEJss2A1agpskW2GYGDbJFtRrgUbYpsi+I8pjNkG8uXd7a8s82KMGSLbLMCpmBjZFsQZk1XyDaWL7JFtlkRhmyRbVbAFGyMbAvCRLaxMGt7p4xMGTkj7JAtss0Il6JNkW1RnJSRY3HWPc5I2qqVlRUtLy833nuv15tbZ83tgTcZLWYjN6E23jZMkBqPU9NWyLYpufG2o4w8HqfmrchskW1G9JDZktlmhEvRpsi2KE4y21icZLY1BJBtRtAhW2SbES5FmyLbojiRbSxOZItsl46srq42DjNki2wbB8+EGyLbCQFusDll5Fi+zEZmNnJWhCFbZJsVMAUbI9uCMOvyroMHtXXnzokm8FT3B2lFUvNJQLFnOqveeWdLGTkj9pAtss0Il6JNkW1RnJSRY3FSRqaMTBk56hpjNnIU2fV+kW0sX8rIsXwpI1NGzoowMlsy26yAKdgY2RaESRk5FmZt75SRKSNnhB2yRbYZ4VK0KbItipMycixOysiUkSkjR11jlJGjyFJGjiW73jtl5GjKZLZkthkxRmZLZpsRLkWbktkWxUlmG4uTzJbMlsw26hojs40iS2YbS5bMdhp8mSDFBKmsOCOzJbPNCpiCjclsC8Ksy7v4nG0sYFFGpoycEWLIFtlmhEvRpsi2KE7KyLE4KSNTRqaMHHWNUUaOIksZOZYsZeRp8KWMTBk5K87IbMlsswKmYGMy24IwKSPHwqztnTIyZeSMsEO2yDYjXIo2RbZFcVJGjsVJGZkyMmXkqGuMMnIUWcrIsWQpI0+DL2VkyshZcUZmS2abFTAFG5PZFoRJGTkWJmXkWgKUkTPCDtki24xwKdoU2RbFSRk5FidlZMrIlJGjrjHKyFFkKSPHkqWMPA2+lJEpI2fFGZktmW1WwBRsTGZbECZl5FiYlJEpIy8tkdlGXWVktlFkyWxjyZLZToMvmS2ZbVackdmS2WYFTMHGZLYFYZLZxsIksyWzJbONu8bIbOPYumdkG8uXn9iL5UtmS2abFWFktmS2WQFTsDGyLQiTzDYWJpktmS2Zbdw1RmYbx5bMNpateyezjWbM1zXyOduMGCOzJbPNCJeiTclsi+I8pjNkG8uXMjJl5KwIQ7bINitgCjZGtgVhUkaOhUkZmTIyZeS4a4wychxbysixbCkjx/Mls53vzPbxkm6UdLKk+yS9UNIdo8IG2cZdVMg2ji2yjWWLbOP5Itv5lu1HJb1H0gckXSHpKklPRbbTuHCO3QeyjeVOGTmWL+9sY/ki2/mV7SmSbk9Z7dckLaXs9hxJXxgWNmS2cRcUso1jS2Yby5bMNp4vsp1f2T5Z0vskndcXJntShvuJaNnuk7Q8jfico31YttskraysaHm5GZ1qAtpNN+3Tli3N+pgjZFmHevDgAV1++bYifPfdcIOWt2zJ2v+iN3Zmu+3qq4vwlbhD1Ne+Jotf99nr9eb2EzTzeuDjyvZaSf6rljMkORNelOU4SWuLcjItPA/4xg8KjGMZLxrfVUnHxyKL6X1eZduojByDcKa93iXpzJkewWLvHL7x4wvjWMbwjeU7du/zKluf4C2S3t03QerFki4a+8wXoyEXUuw4wjeWr3uHcSxj+MbyHbv3eZbtE9JHfx4h6Uvpoz+fHvvMF6MhF1LsOMI3li+yhW88gZbsYZ5l2xKEMz0Mv4/+5ZkewWLvHL7x4wvjWMbwjeU7du/IdmxUNIQABCAAAQg0I4Bsm3FjKwhAAAIQgMDYBJDt2KhoCAEIQAACEGhGANk24xa11VskXSbpLEnbJe1OO7pY0i9I+gZJ+yW9SNJfbbDO3xn9fklnSzqUvnHrpembtqKOv+39luTrc71Z0mPTZ52/Iunlkj7ZdgiBx1eab3Wor5f0hoFrIvA0Wtt1ab57JX01/fmkXyXpD1p79nN+YMi2XQP4FEl3Svq4pEuSbD3b2j+wcGES5uWSdkp6hqRR606S9MT0ESmf5ZslPVzSj7brlKd6NCX5+sAfJunL6QwulfQ6Sf7Cla4upfma45MkXSfpXEl+6KweQLvIuDRfy7a6z3SR51TPGdlOFffYO+u/CC6Q9Kt9X01pifpXjh6Vsqph6744sDdnzC+T9L1jH8XiNozg+wJJr5D0rYuLbewzK8V3c3pYfL6kjyGGB/iX4otsxw7pyRsi28kZRvTQfxE4e/pLSd8v6c8l/Zikt6cMyv992Lr+74j2V7Z9RNKHJF0fccBz1mdJvv5ilaf7a1tTteFTc8Yi4nBL8X1T+tILxztieHCkSvF1P67MOHZvlfRqSX8TERD0Ob8/RLDoYzd4Y/FPB/5c+k7Q30sZ1HdLui39rOCwdRUn36z8OwHP5buUjyIpzdd9OrN1if+Zix6cY5xfCb4PTeXjp0k6gmy/jnoJvr53+LviPyvJFYS3pd9X8c+VsgQQILMNgFqgy1FP8Y9O73U9AepvB/ZVt86TKvzu9llpolSBw5v7Lkry7YdxUNLpHZ+EVvcw089o3Pj1ZLNrJB1OG5vrvZKuTFWauQ/CCU4gIn49IdOVL0/4YwkggGwDoBbocvBiOlXSPZJcDv4Vrf845E+k/Yxa54klnljliSX3FziuRemiBN8T0+8p+ysdvTxb0juSbJ2JdXkpwXeQH2XkB4mU4Ov4fUhfGfmNaRLac7ocuJHnjmwj6eb37fepLvVaoJ7g5IlQ/s1evxf8Tkmb0lO9v4LNU/a9DFvn7Txz83ZJzri8fEbS8/IPa2G2KMnXE9Q+LMk3Lf9soz+S9UpJQ39PeWEoDj+RknyR7bGcS/J9nKQPSlpKf/49cFcSPteBOJ3JKSLbmWBnpxCAAAQg0CUCyLZLo825QgACEIDATAgg25lgZ6cQgAAEINAlAsi2S6PNuUIAAhCAwEwIINuZYGenEIAABCDQJQLItkujzblCAAIQgMBMCCDbmWBnpxCAAAQg0CUCyLZLo825QgACEIDATAgg25lgZ6cQgAAEINAlAsi2S6PNuUIAAhCAwEwIINuZYGenEIAABCDQJQLItkujzblCAAIQgMBMCCDbmWBnpxCAAAQg0CUCyLZLo825QgACEIDATAgg25lgZ6cQgAAEINAlAsi2S6PNuUIAAhCAwEwIINuZYGenEIAABCDQJQLItkujzblCAAIQgMBMCCDbmWBnpxCAAAQg0CUCyLZLo825QgACEIDATAgg25lgZ6cQgAAEINAlAsi2S6PNuUIAAhCAwEwIINuZYGenEIAABCDQJQLItkujzblCAAIQgMBMCPx/VpxwlWuz/doAAAAASUVORK5CYII=\" width=\"431.818172458775\">"
+      ],
       "text/plain": [
-       "<matplotlib.text.Text at 0x7f4b45de3198>"
+       "<IPython.core.display.HTML object>"
       ]
      },
-     "execution_count": 57,
      "metadata": {},
-     "output_type": "execute_result"
+     "output_type": "display_data"
     }
    ],
    "source": [
-    "# Use the following data for this assignment:\n",
+    "#Initialize figure and axes:\n",
+    "fig, ax = plt.subplots()\n",
     "\n",
-    "import pandas as pd\n",
-    "import numpy as np\n",
-    "import matplotlib.pyplot as plt\n",
+    "#Set thereshold (manual input):\n",
+    "threshold = 41000\n",
     "\n",
-    "np.random.seed(12345)\n",
+    "#Set the bar colour:\n",
+    "colormap = matplotlib.colors.LinearSegmentedColormap.from_list(\"colormap\",[\"blue\", \"white\", \"red\"])\n",
+    "def colour_list(threshold):\n",
+    "    colour_list=[]\n",
+    "    for l,u in zip(lower_bound,upper_bound):\n",
+    "        cnorm = matplotlib.colors.Normalize(vmin=l,vmax=u) # normalize based on l and u\n",
+    "        cpick = matplotlib.cm.ScalarMappable(norm=cnorm,cmap=colormap) # scale it based on norm and cmap\n",
+    "        cpick.set_array([]) # set an array\n",
+    "        colour_list.append(cpick.to_rgba(threshold)) #transform it to colour and append it to norm\n",
+    "    return colour_list\n",
     "\n",
-    "df = pd.DataFrame([np.random.normal(32000,200000,3650), \n",
-    "                   np.random.normal(43000,100000,3650), \n",
-    "                   np.random.normal(43500,140000,3650), \n",
-    "                   np.random.normal(48000,70000,3650)], \n",
-    "                  index=[1992,1993,1994,1995])\n",
+    "#plot\n",
+    "plot(colour_list(threshold))\n",
+    "    \n",
+    "#show figure:\n",
+    "plt.show()"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "**Even Harder option implementation:** Add interactivity to the above, which allows the user to click on the y axis to set the value of interest. The bar colors should change with respect to what value the user has selected."
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 243,
+   "metadata": {},
+   "outputs": [
+    {
+     "data": {
+      "application/javascript": [
+       "/* Put everything inside the global mpl namespace */\n",
+       "window.mpl = {};\n",
+       "\n",
+       "\n",
+       "mpl.get_websocket_type = function() {\n",
+       "    if (typeof(WebSocket) !== 'undefined') {\n",
+       "        return WebSocket;\n",
+       "    } else if (typeof(MozWebSocket) !== 'undefined') {\n",
+       "        return MozWebSocket;\n",
+       "    } else {\n",
+       "        alert('Your browser does not have WebSocket support. ' +\n",
+       "              'Please try Chrome, Safari or Firefox ≥ 6. ' +\n",
+       "              'Firefox 4 and 5 are also supported but you ' +\n",
+       "              'have to enable WebSockets in about:config.');\n",
+       "    };\n",
+       "}\n",
+       "\n",
+       "mpl.figure = function(figure_id, websocket, ondownload, parent_element) {\n",
+       "    this.id = figure_id;\n",
+       "\n",
+       "    this.ws = websocket;\n",
+       "\n",
+       "    this.supports_binary = (this.ws.binaryType != undefined);\n",
+       "\n",
+       "    if (!this.supports_binary) {\n",
+       "        var warnings = document.getElementById(\"mpl-warnings\");\n",
+       "        if (warnings) {\n",
+       "            warnings.style.display = 'block';\n",
+       "            warnings.textContent = (\n",
+       "                \"This browser does not support binary websocket messages. \" +\n",
+       "                    \"Performance may be slow.\");\n",
+       "        }\n",
+       "    }\n",
+       "\n",
+       "    this.imageObj = new Image();\n",
+       "\n",
+       "    this.context = undefined;\n",
+       "    this.message = undefined;\n",
+       "    this.canvas = undefined;\n",
+       "    this.rubberband_canvas = undefined;\n",
+       "    this.rubberband_context = undefined;\n",
+       "    this.format_dropdown = undefined;\n",
+       "\n",
+       "    this.image_mode = 'full';\n",
+       "\n",
+       "    this.root = $('<div/>');\n",
+       "    this._root_extra_style(this.root)\n",
+       "    this.root.attr('style', 'display: inline-block');\n",
+       "\n",
+       "    $(parent_element).append(this.root);\n",
+       "\n",
+       "    this._init_header(this);\n",
+       "    this._init_canvas(this);\n",
+       "    this._init_toolbar(this);\n",
+       "\n",
+       "    var fig = this;\n",
+       "\n",
+       "    this.waiting = false;\n",
+       "\n",
+       "    this.ws.onopen =  function () {\n",
+       "            fig.send_message(\"supports_binary\", {value: fig.supports_binary});\n",
+       "            fig.send_message(\"send_image_mode\", {});\n",
+       "            if (mpl.ratio != 1) {\n",
+       "                fig.send_message(\"set_dpi_ratio\", {'dpi_ratio': mpl.ratio});\n",
+       "            }\n",
+       "            fig.send_message(\"refresh\", {});\n",
+       "        }\n",
+       "\n",
+       "    this.imageObj.onload = function() {\n",
+       "            if (fig.image_mode == 'full') {\n",
+       "                // Full images could contain transparency (where diff images\n",
+       "                // almost always do), so we need to clear the canvas so that\n",
+       "                // there is no ghosting.\n",
+       "                fig.context.clearRect(0, 0, fig.canvas.width, fig.canvas.height);\n",
+       "            }\n",
+       "            fig.context.drawImage(fig.imageObj, 0, 0);\n",
+       "        };\n",
+       "\n",
+       "    this.imageObj.onunload = function() {\n",
+       "        fig.ws.close();\n",
+       "    }\n",
+       "\n",
+       "    this.ws.onmessage = this._make_on_message_function(this);\n",
+       "\n",
+       "    this.ondownload = ondownload;\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype._init_header = function() {\n",
+       "    var titlebar = $(\n",
+       "        '<div class=\"ui-dialog-titlebar ui-widget-header ui-corner-all ' +\n",
+       "        'ui-helper-clearfix\"/>');\n",
+       "    var titletext = $(\n",
+       "        '<div class=\"ui-dialog-title\" style=\"width: 100%; ' +\n",
+       "        'text-align: center; padding: 3px;\"/>');\n",
+       "    titlebar.append(titletext)\n",
+       "    this.root.append(titlebar);\n",
+       "    this.header = titletext[0];\n",
+       "}\n",
+       "\n",
+       "\n",
+       "\n",
+       "mpl.figure.prototype._canvas_extra_style = function(canvas_div) {\n",
+       "\n",
+       "}\n",
+       "\n",
+       "\n",
+       "mpl.figure.prototype._root_extra_style = function(canvas_div) {\n",
+       "\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype._init_canvas = function() {\n",
+       "    var fig = this;\n",
+       "\n",
+       "    var canvas_div = $('<div/>');\n",
+       "\n",
+       "    canvas_div.attr('style', 'position: relative; clear: both; outline: 0');\n",
+       "\n",
+       "    function canvas_keyboard_event(event) {\n",
+       "        return fig.key_event(event, event['data']);\n",
+       "    }\n",
+       "\n",
+       "    canvas_div.keydown('key_press', canvas_keyboard_event);\n",
+       "    canvas_div.keyup('key_release', canvas_keyboard_event);\n",
+       "    this.canvas_div = canvas_div\n",
+       "    this._canvas_extra_style(canvas_div)\n",
+       "    this.root.append(canvas_div);\n",
+       "\n",
+       "    var canvas = $('<canvas/>');\n",
+       "    canvas.addClass('mpl-canvas');\n",
+       "    canvas.attr('style', \"left: 0; top: 0; z-index: 0; outline: 0\")\n",
+       "\n",
+       "    this.canvas = canvas[0];\n",
+       "    this.context = canvas[0].getContext(\"2d\");\n",
+       "\n",
+       "    var backingStore = this.context.backingStorePixelRatio ||\n",
+       "\tthis.context.webkitBackingStorePixelRatio ||\n",
+       "\tthis.context.mozBackingStorePixelRatio ||\n",
+       "\tthis.context.msBackingStorePixelRatio ||\n",
+       "\tthis.context.oBackingStorePixelRatio ||\n",
+       "\tthis.context.backingStorePixelRatio || 1;\n",
+       "\n",
+       "    mpl.ratio = (window.devicePixelRatio || 1) / backingStore;\n",
+       "\n",
+       "    var rubberband = $('<canvas/>');\n",
+       "    rubberband.attr('style', \"position: absolute; left: 0; top: 0; z-index: 1;\")\n",
+       "\n",
+       "    var pass_mouse_events = true;\n",
+       "\n",
+       "    canvas_div.resizable({\n",
+       "        start: function(event, ui) {\n",
+       "            pass_mouse_events = false;\n",
+       "        },\n",
+       "        resize: function(event, ui) {\n",
+       "            fig.request_resize(ui.size.width, ui.size.height);\n",
+       "        },\n",
+       "        stop: function(event, ui) {\n",
+       "            pass_mouse_events = true;\n",
+       "            fig.request_resize(ui.size.width, ui.size.height);\n",
+       "        },\n",
+       "    });\n",
+       "\n",
+       "    function mouse_event_fn(event) {\n",
+       "        if (pass_mouse_events)\n",
+       "            return fig.mouse_event(event, event['data']);\n",
+       "    }\n",
+       "\n",
+       "    rubberband.mousedown('button_press', mouse_event_fn);\n",
+       "    rubberband.mouseup('button_release', mouse_event_fn);\n",
+       "    // Throttle sequential mouse events to 1 every 20ms.\n",
+       "    rubberband.mousemove('motion_notify', mouse_event_fn);\n",
+       "\n",
+       "    rubberband.mouseenter('figure_enter', mouse_event_fn);\n",
+       "    rubberband.mouseleave('figure_leave', mouse_event_fn);\n",
+       "\n",
+       "    canvas_div.on(\"wheel\", function (event) {\n",
+       "        event = event.originalEvent;\n",
+       "        event['data'] = 'scroll'\n",
+       "        if (event.deltaY < 0) {\n",
+       "            event.step = 1;\n",
+       "        } else {\n",
+       "            event.step = -1;\n",
+       "        }\n",
+       "        mouse_event_fn(event);\n",
+       "    });\n",
+       "\n",
+       "    canvas_div.append(canvas);\n",
+       "    canvas_div.append(rubberband);\n",
+       "\n",
+       "    this.rubberband = rubberband;\n",
+       "    this.rubberband_canvas = rubberband[0];\n",
+       "    this.rubberband_context = rubberband[0].getContext(\"2d\");\n",
+       "    this.rubberband_context.strokeStyle = \"#000000\";\n",
+       "\n",
+       "    this._resize_canvas = function(width, height) {\n",
+       "        // Keep the size of the canvas, canvas container, and rubber band\n",
+       "        // canvas in synch.\n",
+       "        canvas_div.css('width', width)\n",
+       "        canvas_div.css('height', height)\n",
+       "\n",
+       "        canvas.attr('width', width * mpl.ratio);\n",
+       "        canvas.attr('height', height * mpl.ratio);\n",
+       "        canvas.attr('style', 'width: ' + width + 'px; height: ' + height + 'px;');\n",
+       "\n",
+       "        rubberband.attr('width', width);\n",
+       "        rubberband.attr('height', height);\n",
+       "    }\n",
+       "\n",
+       "    // Set the figure to an initial 600x600px, this will subsequently be updated\n",
+       "    // upon first draw.\n",
+       "    this._resize_canvas(600, 600);\n",
+       "\n",
+       "    // Disable right mouse context menu.\n",
+       "    $(this.rubberband_canvas).bind(\"contextmenu\",function(e){\n",
+       "        return false;\n",
+       "    });\n",
+       "\n",
+       "    function set_focus () {\n",
+       "        canvas.focus();\n",
+       "        canvas_div.focus();\n",
+       "    }\n",
+       "\n",
+       "    window.setTimeout(set_focus, 100);\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype._init_toolbar = function() {\n",
+       "    var fig = this;\n",
+       "\n",
+       "    var nav_element = $('<div/>');\n",
+       "    nav_element.attr('style', 'width: 100%');\n",
+       "    this.root.append(nav_element);\n",
+       "\n",
+       "    // Define a callback function for later on.\n",
+       "    function toolbar_event(event) {\n",
+       "        return fig.toolbar_button_onclick(event['data']);\n",
+       "    }\n",
+       "    function toolbar_mouse_event(event) {\n",
+       "        return fig.toolbar_button_onmouseover(event['data']);\n",
+       "    }\n",
+       "\n",
+       "    for(var toolbar_ind in mpl.toolbar_items) {\n",
+       "        var name = mpl.toolbar_items[toolbar_ind][0];\n",
+       "        var tooltip = mpl.toolbar_items[toolbar_ind][1];\n",
+       "        var image = mpl.toolbar_items[toolbar_ind][2];\n",
+       "        var method_name = mpl.toolbar_items[toolbar_ind][3];\n",
+       "\n",
+       "        if (!name) {\n",
+       "            // put a spacer in here.\n",
+       "            continue;\n",
+       "        }\n",
+       "        var button = $('<button/>');\n",
+       "        button.addClass('ui-button ui-widget ui-state-default ui-corner-all ' +\n",
+       "                        'ui-button-icon-only');\n",
+       "        button.attr('role', 'button');\n",
+       "        button.attr('aria-disabled', 'false');\n",
+       "        button.click(method_name, toolbar_event);\n",
+       "        button.mouseover(tooltip, toolbar_mouse_event);\n",
+       "\n",
+       "        var icon_img = $('<span/>');\n",
+       "        icon_img.addClass('ui-button-icon-primary ui-icon');\n",
+       "        icon_img.addClass(image);\n",
+       "        icon_img.addClass('ui-corner-all');\n",
+       "\n",
+       "        var tooltip_span = $('<span/>');\n",
+       "        tooltip_span.addClass('ui-button-text');\n",
+       "        tooltip_span.html(tooltip);\n",
+       "\n",
+       "        button.append(icon_img);\n",
+       "        button.append(tooltip_span);\n",
+       "\n",
+       "        nav_element.append(button);\n",
+       "    }\n",
+       "\n",
+       "    var fmt_picker_span = $('<span/>');\n",
+       "\n",
+       "    var fmt_picker = $('<select/>');\n",
+       "    fmt_picker.addClass('mpl-toolbar-option ui-widget ui-widget-content');\n",
+       "    fmt_picker_span.append(fmt_picker);\n",
+       "    nav_element.append(fmt_picker_span);\n",
+       "    this.format_dropdown = fmt_picker[0];\n",
+       "\n",
+       "    for (var ind in mpl.extensions) {\n",
+       "        var fmt = mpl.extensions[ind];\n",
+       "        var option = $(\n",
+       "            '<option/>', {selected: fmt === mpl.default_extension}).html(fmt);\n",
+       "        fmt_picker.append(option);\n",
+       "    }\n",
+       "\n",
+       "    // Add hover states to the ui-buttons\n",
+       "    $( \".ui-button\" ).hover(\n",
+       "        function() { $(this).addClass(\"ui-state-hover\");},\n",
+       "        function() { $(this).removeClass(\"ui-state-hover\");}\n",
+       "    );\n",
+       "\n",
+       "    var status_bar = $('<span class=\"mpl-message\"/>');\n",
+       "    nav_element.append(status_bar);\n",
+       "    this.message = status_bar[0];\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype.request_resize = function(x_pixels, y_pixels) {\n",
+       "    // Request matplotlib to resize the figure. Matplotlib will then trigger a resize in the client,\n",
+       "    // which will in turn request a refresh of the image.\n",
+       "    this.send_message('resize', {'width': x_pixels, 'height': y_pixels});\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype.send_message = function(type, properties) {\n",
+       "    properties['type'] = type;\n",
+       "    properties['figure_id'] = this.id;\n",
+       "    this.ws.send(JSON.stringify(properties));\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype.send_draw_message = function() {\n",
+       "    if (!this.waiting) {\n",
+       "        this.waiting = true;\n",
+       "        this.ws.send(JSON.stringify({type: \"draw\", figure_id: this.id}));\n",
+       "    }\n",
+       "}\n",
+       "\n",
+       "\n",
+       "mpl.figure.prototype.handle_save = function(fig, msg) {\n",
+       "    var format_dropdown = fig.format_dropdown;\n",
+       "    var format = format_dropdown.options[format_dropdown.selectedIndex].value;\n",
+       "    fig.ondownload(fig, format);\n",
+       "}\n",
+       "\n",
+       "\n",
+       "mpl.figure.prototype.handle_resize = function(fig, msg) {\n",
+       "    var size = msg['size'];\n",
+       "    if (size[0] != fig.canvas.width || size[1] != fig.canvas.height) {\n",
+       "        fig._resize_canvas(size[0], size[1]);\n",
+       "        fig.send_message(\"refresh\", {});\n",
+       "    };\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype.handle_rubberband = function(fig, msg) {\n",
+       "    var x0 = msg['x0'] / mpl.ratio;\n",
+       "    var y0 = (fig.canvas.height - msg['y0']) / mpl.ratio;\n",
+       "    var x1 = msg['x1'] / mpl.ratio;\n",
+       "    var y1 = (fig.canvas.height - msg['y1']) / mpl.ratio;\n",
+       "    x0 = Math.floor(x0) + 0.5;\n",
+       "    y0 = Math.floor(y0) + 0.5;\n",
+       "    x1 = Math.floor(x1) + 0.5;\n",
+       "    y1 = Math.floor(y1) + 0.5;\n",
+       "    var min_x = Math.min(x0, x1);\n",
+       "    var min_y = Math.min(y0, y1);\n",
+       "    var width = Math.abs(x1 - x0);\n",
+       "    var height = Math.abs(y1 - y0);\n",
+       "\n",
+       "    fig.rubberband_context.clearRect(\n",
+       "        0, 0, fig.canvas.width / mpl.ratio, fig.canvas.height / mpl.ratio);\n",
+       "\n",
+       "    fig.rubberband_context.strokeRect(min_x, min_y, width, height);\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype.handle_figure_label = function(fig, msg) {\n",
+       "    // Updates the figure title.\n",
+       "    fig.header.textContent = msg['label'];\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype.handle_cursor = function(fig, msg) {\n",
+       "    var cursor = msg['cursor'];\n",
+       "    switch(cursor)\n",
+       "    {\n",
+       "    case 0:\n",
+       "        cursor = 'pointer';\n",
+       "        break;\n",
+       "    case 1:\n",
+       "        cursor = 'default';\n",
+       "        break;\n",
+       "    case 2:\n",
+       "        cursor = 'crosshair';\n",
+       "        break;\n",
+       "    case 3:\n",
+       "        cursor = 'move';\n",
+       "        break;\n",
+       "    }\n",
+       "    fig.rubberband_canvas.style.cursor = cursor;\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype.handle_message = function(fig, msg) {\n",
+       "    fig.message.textContent = msg['message'];\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype.handle_draw = function(fig, msg) {\n",
+       "    // Request the server to send over a new figure.\n",
+       "    fig.send_draw_message();\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype.handle_image_mode = function(fig, msg) {\n",
+       "    fig.image_mode = msg['mode'];\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype.updated_canvas_event = function() {\n",
+       "    // Called whenever the canvas gets updated.\n",
+       "    this.send_message(\"ack\", {});\n",
+       "}\n",
+       "\n",
+       "// A function to construct a web socket function for onmessage handling.\n",
+       "// Called in the figure constructor.\n",
+       "mpl.figure.prototype._make_on_message_function = function(fig) {\n",
+       "    return function socket_on_message(evt) {\n",
+       "        if (evt.data instanceof Blob) {\n",
+       "            /* FIXME: We get \"Resource interpreted as Image but\n",
+       "             * transferred with MIME type text/plain:\" errors on\n",
+       "             * Chrome.  But how to set the MIME type?  It doesn't seem\n",
+       "             * to be part of the websocket stream */\n",
+       "            evt.data.type = \"image/png\";\n",
+       "\n",
+       "            /* Free the memory for the previous frames */\n",
+       "            if (fig.imageObj.src) {\n",
+       "                (window.URL || window.webkitURL).revokeObjectURL(\n",
+       "                    fig.imageObj.src);\n",
+       "            }\n",
+       "\n",
+       "            fig.imageObj.src = (window.URL || window.webkitURL).createObjectURL(\n",
+       "                evt.data);\n",
+       "            fig.updated_canvas_event();\n",
+       "            fig.waiting = false;\n",
+       "            return;\n",
+       "        }\n",
+       "        else if (typeof evt.data === 'string' && evt.data.slice(0, 21) == \"data:image/png;base64\") {\n",
+       "            fig.imageObj.src = evt.data;\n",
+       "            fig.updated_canvas_event();\n",
+       "            fig.waiting = false;\n",
+       "            return;\n",
+       "        }\n",
+       "\n",
+       "        var msg = JSON.parse(evt.data);\n",
+       "        var msg_type = msg['type'];\n",
+       "\n",
+       "        // Call the  \"handle_{type}\" callback, which takes\n",
+       "        // the figure and JSON message as its only arguments.\n",
+       "        try {\n",
+       "            var callback = fig[\"handle_\" + msg_type];\n",
+       "        } catch (e) {\n",
+       "            console.log(\"No handler for the '\" + msg_type + \"' message type: \", msg);\n",
+       "            return;\n",
+       "        }\n",
+       "\n",
+       "        if (callback) {\n",
+       "            try {\n",
+       "                // console.log(\"Handling '\" + msg_type + \"' message: \", msg);\n",
+       "                callback(fig, msg);\n",
+       "            } catch (e) {\n",
+       "                console.log(\"Exception inside the 'handler_\" + msg_type + \"' callback:\", e, e.stack, msg);\n",
+       "            }\n",
+       "        }\n",
+       "    };\n",
+       "}\n",
+       "\n",
+       "// from http://stackoverflow.com/questions/1114465/getting-mouse-location-in-canvas\n",
+       "mpl.findpos = function(e) {\n",
+       "    //this section is from http://www.quirksmode.org/js/events_properties.html\n",
+       "    var targ;\n",
+       "    if (!e)\n",
+       "        e = window.event;\n",
+       "    if (e.target)\n",
+       "        targ = e.target;\n",
+       "    else if (e.srcElement)\n",
+       "        targ = e.srcElement;\n",
+       "    if (targ.nodeType == 3) // defeat Safari bug\n",
+       "        targ = targ.parentNode;\n",
+       "\n",
+       "    // jQuery normalizes the pageX and pageY\n",
+       "    // pageX,Y are the mouse positions relative to the document\n",
+       "    // offset() returns the position of the element relative to the document\n",
+       "    var x = e.pageX - $(targ).offset().left;\n",
+       "    var y = e.pageY - $(targ).offset().top;\n",
+       "\n",
+       "    return {\"x\": x, \"y\": y};\n",
+       "};\n",
+       "\n",
+       "/*\n",
+       " * return a copy of an object with only non-object keys\n",
+       " * we need this to avoid circular references\n",
+       " * http://stackoverflow.com/a/24161582/3208463\n",
+       " */\n",
+       "function simpleKeys (original) {\n",
+       "  return Object.keys(original).reduce(function (obj, key) {\n",
+       "    if (typeof original[key] !== 'object')\n",
+       "        obj[key] = original[key]\n",
+       "    return obj;\n",
+       "  }, {});\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype.mouse_event = function(event, name) {\n",
+       "    var canvas_pos = mpl.findpos(event)\n",
+       "\n",
+       "    if (name === 'button_press')\n",
+       "    {\n",
+       "        this.canvas.focus();\n",
+       "        this.canvas_div.focus();\n",
+       "    }\n",
+       "\n",
+       "    var x = canvas_pos.x * mpl.ratio;\n",
+       "    var y = canvas_pos.y * mpl.ratio;\n",
+       "\n",
+       "    this.send_message(name, {x: x, y: y, button: event.button,\n",
+       "                             step: event.step,\n",
+       "                             guiEvent: simpleKeys(event)});\n",
+       "\n",
+       "    /* This prevents the web browser from automatically changing to\n",
+       "     * the text insertion cursor when the button is pressed.  We want\n",
+       "     * to control all of the cursor setting manually through the\n",
+       "     * 'cursor' event from matplotlib */\n",
+       "    event.preventDefault();\n",
+       "    return false;\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype._key_event_extra = function(event, name) {\n",
+       "    // Handle any extra behaviour associated with a key event\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype.key_event = function(event, name) {\n",
+       "\n",
+       "    // Prevent repeat events\n",
+       "    if (name == 'key_press')\n",
+       "    {\n",
+       "        if (event.which === this._key)\n",
+       "            return;\n",
+       "        else\n",
+       "            this._key = event.which;\n",
+       "    }\n",
+       "    if (name == 'key_release')\n",
+       "        this._key = null;\n",
+       "\n",
+       "    var value = '';\n",
+       "    if (event.ctrlKey && event.which != 17)\n",
+       "        value += \"ctrl+\";\n",
+       "    if (event.altKey && event.which != 18)\n",
+       "        value += \"alt+\";\n",
+       "    if (event.shiftKey && event.which != 16)\n",
+       "        value += \"shift+\";\n",
+       "\n",
+       "    value += 'k';\n",
+       "    value += event.which.toString();\n",
+       "\n",
+       "    this._key_event_extra(event, name);\n",
+       "\n",
+       "    this.send_message(name, {key: value,\n",
+       "                             guiEvent: simpleKeys(event)});\n",
+       "    return false;\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype.toolbar_button_onclick = function(name) {\n",
+       "    if (name == 'download') {\n",
+       "        this.handle_save(this, null);\n",
+       "    } else {\n",
+       "        this.send_message(\"toolbar_button\", {name: name});\n",
+       "    }\n",
+       "};\n",
+       "\n",
+       "mpl.figure.prototype.toolbar_button_onmouseover = function(tooltip) {\n",
+       "    this.message.textContent = tooltip;\n",
+       "};\n",
+       "mpl.toolbar_items = [[\"Home\", \"Reset original view\", \"fa fa-home icon-home\", \"home\"], [\"Back\", \"Back to previous view\", \"fa fa-arrow-left icon-arrow-left\", \"back\"], [\"Forward\", \"Forward to next view\", \"fa fa-arrow-right icon-arrow-right\", \"forward\"], [\"\", \"\", \"\", \"\"], [\"Pan\", \"Pan axes with left mouse, zoom with right\", \"fa fa-arrows icon-move\", \"pan\"], [\"Zoom\", \"Zoom to rectangle\", \"fa fa-square-o icon-check-empty\", \"zoom\"], [\"\", \"\", \"\", \"\"], [\"Download\", \"Download plot\", \"fa fa-floppy-o icon-save\", \"download\"]];\n",
+       "\n",
+       "mpl.extensions = [\"eps\", \"jpeg\", \"pdf\", \"png\", \"ps\", \"raw\", \"svg\", \"tif\"];\n",
+       "\n",
+       "mpl.default_extension = \"png\";var comm_websocket_adapter = function(comm) {\n",
+       "    // Create a \"websocket\"-like object which calls the given IPython comm\n",
+       "    // object with the appropriate methods. Currently this is a non binary\n",
+       "    // socket, so there is still some room for performance tuning.\n",
+       "    var ws = {};\n",
+       "\n",
+       "    ws.close = function() {\n",
+       "        comm.close()\n",
+       "    };\n",
+       "    ws.send = function(m) {\n",
+       "        //console.log('sending', m);\n",
+       "        comm.send(m);\n",
+       "    };\n",
+       "    // Register the callback with on_msg.\n",
+       "    comm.on_msg(function(msg) {\n",
+       "        //console.log('receiving', msg['content']['data'], msg);\n",
+       "        // Pass the mpl event to the overridden (by mpl) onmessage function.\n",
+       "        ws.onmessage(msg['content']['data'])\n",
+       "    });\n",
+       "    return ws;\n",
+       "}\n",
+       "\n",
+       "mpl.mpl_figure_comm = function(comm, msg) {\n",
+       "    // This is the function which gets called when the mpl process\n",
+       "    // starts-up an IPython Comm through the \"matplotlib\" channel.\n",
+       "\n",
+       "    var id = msg.content.data.id;\n",
+       "    // Get hold of the div created by the display call when the Comm\n",
+       "    // socket was opened in Python.\n",
+       "    var element = $(\"#\" + id);\n",
+       "    var ws_proxy = comm_websocket_adapter(comm)\n",
+       "\n",
+       "    function ondownload(figure, format) {\n",
+       "        window.open(figure.imageObj.src);\n",
+       "    }\n",
+       "\n",
+       "    var fig = new mpl.figure(id, ws_proxy,\n",
+       "                           ondownload,\n",
+       "                           element.get(0));\n",
+       "\n",
+       "    // Call onopen now - mpl needs it, as it is assuming we've passed it a real\n",
+       "    // web socket which is closed, not our websocket->open comm proxy.\n",
+       "    ws_proxy.onopen();\n",
+       "\n",
+       "    fig.parent_element = element.get(0);\n",
+       "    fig.cell_info = mpl.find_output_cell(\"<div id='\" + id + \"'></div>\");\n",
+       "    if (!fig.cell_info) {\n",
+       "        console.error(\"Failed to find cell for figure\", id, fig);\n",
+       "        return;\n",
+       "    }\n",
+       "\n",
+       "    var output_index = fig.cell_info[2]\n",
+       "    var cell = fig.cell_info[0];\n",
+       "\n",
+       "};\n",
+       "\n",
+       "mpl.figure.prototype.handle_close = function(fig, msg) {\n",
+       "    var width = fig.canvas.width/mpl.ratio\n",
+       "    fig.root.unbind('remove')\n",
+       "\n",
+       "    // Update the output cell to use the data from the current canvas.\n",
+       "    fig.push_to_output();\n",
+       "    var dataURL = fig.canvas.toDataURL();\n",
+       "    // Re-enable the keyboard manager in IPython - without this line, in FF,\n",
+       "    // the notebook keyboard shortcuts fail.\n",
+       "    IPython.keyboard_manager.enable()\n",
+       "    $(fig.parent_element).html('<img src=\"' + dataURL + '\" width=\"' + width + '\">');\n",
+       "    fig.close_ws(fig, msg);\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype.close_ws = function(fig, msg){\n",
+       "    fig.send_message('closing', msg);\n",
+       "    // fig.ws.close()\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype.push_to_output = function(remove_interactive) {\n",
+       "    // Turn the data on the canvas into data in the output cell.\n",
+       "    var width = this.canvas.width/mpl.ratio\n",
+       "    var dataURL = this.canvas.toDataURL();\n",
+       "    this.cell_info[1]['text/html'] = '<img src=\"' + dataURL + '\" width=\"' + width + '\">';\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype.updated_canvas_event = function() {\n",
+       "    // Tell IPython that the notebook contents must change.\n",
+       "    IPython.notebook.set_dirty(true);\n",
+       "    this.send_message(\"ack\", {});\n",
+       "    var fig = this;\n",
+       "    // Wait a second, then push the new image to the DOM so\n",
+       "    // that it is saved nicely (might be nice to debounce this).\n",
+       "    setTimeout(function () { fig.push_to_output() }, 1000);\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype._init_toolbar = function() {\n",
+       "    var fig = this;\n",
+       "\n",
+       "    var nav_element = $('<div/>');\n",
+       "    nav_element.attr('style', 'width: 100%');\n",
+       "    this.root.append(nav_element);\n",
+       "\n",
+       "    // Define a callback function for later on.\n",
+       "    function toolbar_event(event) {\n",
+       "        return fig.toolbar_button_onclick(event['data']);\n",
+       "    }\n",
+       "    function toolbar_mouse_event(event) {\n",
+       "        return fig.toolbar_button_onmouseover(event['data']);\n",
+       "    }\n",
+       "\n",
+       "    for(var toolbar_ind in mpl.toolbar_items){\n",
+       "        var name = mpl.toolbar_items[toolbar_ind][0];\n",
+       "        var tooltip = mpl.toolbar_items[toolbar_ind][1];\n",
+       "        var image = mpl.toolbar_items[toolbar_ind][2];\n",
+       "        var method_name = mpl.toolbar_items[toolbar_ind][3];\n",
+       "\n",
+       "        if (!name) { continue; };\n",
+       "\n",
+       "        var button = $('<button class=\"btn btn-default\" href=\"#\" title=\"' + name + '\"><i class=\"fa ' + image + ' fa-lg\"></i></button>');\n",
+       "        button.click(method_name, toolbar_event);\n",
+       "        button.mouseover(tooltip, toolbar_mouse_event);\n",
+       "        nav_element.append(button);\n",
+       "    }\n",
+       "\n",
+       "    // Add the status bar.\n",
+       "    var status_bar = $('<span class=\"mpl-message\" style=\"text-align:right; float: right;\"/>');\n",
+       "    nav_element.append(status_bar);\n",
+       "    this.message = status_bar[0];\n",
+       "\n",
+       "    // Add the close button to the window.\n",
+       "    var buttongrp = $('<div class=\"btn-group inline pull-right\"></div>');\n",
+       "    var button = $('<button class=\"btn btn-mini btn-primary\" href=\"#\" title=\"Stop Interaction\"><i class=\"fa fa-power-off icon-remove icon-large\"></i></button>');\n",
+       "    button.click(function (evt) { fig.handle_close(fig, {}); } );\n",
+       "    button.mouseover('Stop Interaction', toolbar_mouse_event);\n",
+       "    buttongrp.append(button);\n",
+       "    var titlebar = this.root.find($('.ui-dialog-titlebar'));\n",
+       "    titlebar.prepend(buttongrp);\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype._root_extra_style = function(el){\n",
+       "    var fig = this\n",
+       "    el.on(\"remove\", function(){\n",
+       "\tfig.close_ws(fig, {});\n",
+       "    });\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype._canvas_extra_style = function(el){\n",
+       "    // this is important to make the div 'focusable\n",
+       "    el.attr('tabindex', 0)\n",
+       "    // reach out to IPython and tell the keyboard manager to turn it's self\n",
+       "    // off when our div gets focus\n",
+       "\n",
+       "    // location in version 3\n",
+       "    if (IPython.notebook.keyboard_manager) {\n",
+       "        IPython.notebook.keyboard_manager.register_events(el);\n",
+       "    }\n",
+       "    else {\n",
+       "        // location in version 2\n",
+       "        IPython.keyboard_manager.register_events(el);\n",
+       "    }\n",
+       "\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype._key_event_extra = function(event, name) {\n",
+       "    var manager = IPython.notebook.keyboard_manager;\n",
+       "    if (!manager)\n",
+       "        manager = IPython.keyboard_manager;\n",
+       "\n",
+       "    // Check for shift+enter\n",
+       "    if (event.shiftKey && event.which == 13) {\n",
+       "        this.canvas_div.blur();\n",
+       "        // select the cell after this one\n",
+       "        var index = IPython.notebook.find_cell_index(this.cell_info[0]);\n",
+       "        IPython.notebook.select(index + 1);\n",
+       "    }\n",
+       "}\n",
+       "\n",
+       "mpl.figure.prototype.handle_save = function(fig, msg) {\n",
+       "    fig.ondownload(fig, null);\n",
+       "}\n",
+       "\n",
+       "\n",
+       "mpl.find_output_cell = function(html_output) {\n",
+       "    // Return the cell and output element which can be found *uniquely* in the notebook.\n",
+       "    // Note - this is a bit hacky, but it is done because the \"notebook_saving.Notebook\"\n",
+       "    // IPython event is triggered only after the cells have been serialised, which for\n",
+       "    // our purposes (turning an active figure into a static one), is too late.\n",
+       "    var cells = IPython.notebook.get_cells();\n",
+       "    var ncells = cells.length;\n",
+       "    for (var i=0; i<ncells; i++) {\n",
+       "        var cell = cells[i];\n",
+       "        if (cell.cell_type === 'code'){\n",
+       "            for (var j=0; j<cell.output_area.outputs.length; j++) {\n",
+       "                var data = cell.output_area.outputs[j];\n",
+       "                if (data.data) {\n",
+       "                    // IPython >= 3 moved mimebundle to data attribute of output\n",
+       "                    data = data.data;\n",
+       "                }\n",
+       "                if (data['text/html'] == html_output) {\n",
+       "                    return [cell, data, j];\n",
+       "                }\n",
+       "            }\n",
+       "        }\n",
+       "    }\n",
+       "}\n",
+       "\n",
+       "// Register the function which deals with the matplotlib target/channel.\n",
+       "// The kernel may be null if the page has been refreshed.\n",
+       "if (IPython.notebook.kernel != null) {\n",
+       "    IPython.notebook.kernel.comm_manager.register_target('matplotlib', mpl.mpl_figure_comm);\n",
+       "}\n"
+      ],
+      "text/plain": [
+       "<IPython.core.display.Javascript object>"
+      ]
+     },
+     "metadata": {},
+     "output_type": "display_data"
+    },
+    {
+     "data": {
+      "text/html": [
+       "<img src=\"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAdsAAAE8CAYAAACIOxsOAAAgAElEQVR4Xu2dD7hkZ13fv5OLu9mEu5AEYsxm40KaiEkWFYxdI9UGlCIJPECQsjQQdlERa4Jp+SMUBBVaIVQxlFIqTQxQSf8gNGgVVIiStmtVpGbZlgR1k3RD/pBF7mI2u97c9fnufU8yzJ6ZO++Z9zdzZuZznuc+Tzbznvec83l/53zO73feM9MRCwQgAAEIQAACoQQ6ob3TOQQgAAEIQAACQrYEAQQgAAEIQCCYALINBkz3EIAABCAAAWRLDEAAAhCAAASCCSDbYMB0DwEIQAACEEC2xAAEIAABCEAgmACyDQZM9xCAAAQgAAFkSwxAAAIQgAAEggkg22DAdA8BCEAAAhBAtsQABCAAAQhAIJgAsg0GTPcQgAAEIAABZEsMQAACEIAABIIJINtgwHQPAQhAAAIQQLbEAAQgAAEIQCCYALINBkz3EIAABCAAAWRLDEAAAhCAAASCCSDbYMB0DwEIQAACEEC2xAAEIAABCEAgmACyDQZM9xCAAAQgAAFkSwzMIoG9kt4q6ddm8eA4pqMEqrF9OTwgMA0EkO00jFL79/EmST+QdvNBSXdJ+qykX5b0fzJ234L8h+kvY7Vjmo4q2y2S/krSEyS5r2lZXinppZKeIuleST6O7uU4Sa+V9OOSTpf0OUk/2TNGF0p6p6TvlLRf0rskXdPVyRslXSrp2yT9taTfkOT/9/UxQ1pLtt8v6WckXSDpcX3G8hJJv5CO5f9LerOk/9R1HGekGL5Iktn5s6skOca9WPTX9Rz3f5P0vDGzYHNTQADZTsEgTcEuWrZ/ki7M6ySdJelHJf1IujB/YshjQLZDgurT7NWSzP+0xL1XtldK+hdJErdJquRscX5NkuXyBUnvlvQhSeemDNJCviFt87ck/UdJfyrpmyX9e0n/S9KO0XY9e+21ZPvDkr5X0p1pH3tvnHxDskvST0v6pKTvk/R+Sc+Q9D+TXP93unG0tBclvU/SH0l6VZdsLWsLvVosYt+EsEDgGwggWwKiBAHL1n+WZffyAUnPSlnF30p6brrY+yJ+QNKNkl6TsqK6LMEXyIOS3ivpaZJOlPR5SZaKs7J+i7NRr+OMxH/+t7M5Z9vV8oOSflHSeemC7CzcF1MvR3o6/rm0v8+U9I/SZ27vC7WF4yzyqZJulvQYSYclPVHSr6TtL0n6qKTXpeNxFydIulrSiyR9U9q3fyrpjtS/ZeLzc1+6uD8g6R09WWa/4zdLj0WvbC2R35T0L9OK7t8Z3dvSsZuRx+Ocro7d1gL6+3025hsqS+rkAeNheXvsvyVVDCyoj3S19/iYlcfYkvS/LbQ/6GpzhaQ3STpe0n9I2epKunEYFMP9qhQ+ru+W5DGtll+XtCDpHycGX5S0KQnXbdzW/B6fbk76cR60P3w2pwSQ7ZwOfOHD7ifb70hy3JYyAovF8nT25IvYv5P02+kCvyGV9FzGfEHav/skbU4X6s9IOpQE9xxJZ3eV83oPxxfrjamM+IeSLLF/IulbkzSdyTlrsSx9Qbf8r00l1f+aMhV//j1JxC6RWsq/J+mxkh6S9MfpJsKCcinVNwDPTyVwZ5c+xo9L8g2H9+XfpHV+Ku3sB9NF26VL33i8PpV/LW33b9maw79N+/ZDkt4jaWvqe9AQ9pOAb1A+LOmXulb+y3ST8DJJznxdkXhy1+c/m8qrvtHxTUTvYim6NO2bi36Lj9HZ41ckVcdh0f15WsHj5fF3P86SXZZ2O99s+SbNjyh+V5LZeTy9ny6X+wZmrWe2/WRrBn8vxVa1344B34SdmRj4EcgpqZzuNi5NO1588/HptG3faPi4nNE6PpwFf7Xw+UV3M0AA2c7AILbgEPrJ1lmI5frinmdh1S6/MD0frC7Uw5SR/ezMZTo/b/OFt27xxdtZ5mXpQ2crlsq/SoL3RdXZqC+M1eILvC+mzsTrLtDOPr1dt/l/aX1nbM5Q/Rzvv0hyJuTsy+Jy6dU3GdXimwhfjC0tS999nCrJWa+Xqn9fyF3etGxd6uwWn0u//zodw6Bh7ydbZ9J+nvjsxMM3Cpb5p1LG7m1ZyDtTqdg3Jb+Tbnj8jPfLPRt1Fu/25vn2jDh0duibmZ/vkq0fNTh79eIbMWfc3544+Vmppds9nn6mXglv0Kb7yfbidJPkOPr9lLn7hmB9+vMNk2PG+/XP0jg783V2+5KUmXt83b9vrFyCd3x5PH1z0FsdycBD01kkgGxncVTHf0z9ZOsLl+/4XZb7z6k05wuSn3GdJOlR6c+i8VInW7dxdmUx+4Lvf1twzmz87LBusWyduXRP7HH2eXsSo7NSZ4jLXSu7X5dwXULtd4F2du0JMLek0quzUW/HWZonhfm5pS/YlpqzZmfi1eJzzfttkXxXuoi7NNy9+HMLxRd1y9ZcnJFXi7Mqi6GSVL+R7ifbR6cJPc6YLQPfkPhmyIvLt15+LO2/2zpDcwZu2fg58D1dG7SM/Py2WrebZe9+WU7m4Zsq34A5Lnx83pYXj5dl/avp32blErGl5RsqPzqw0HvH02Jrmtl6U46rNyTOfrbrsbs8Zdn+3DdI16f9tuxdevYjhX43j86IfSyuiHgOAwsEHiaAbAmGEgTWKiP7eZ8zGWd+Lh36ounM0hczX0SrOKyTrS+GzhJdpr01ydulRpcc+73a4wueM0CXXaulW7b/V5LLxb6Qdi++oFrI/WTrZ40uOfsYLCOXR+9Px+EMzzcQFoCf/XqSmPe7d/G+eTavJxa5ZNy7WGguK9dNAOrHubePtZ4lumTrrPTuNBnI5dJq0o/78nj4+arL+L6B8A2F27u87cU3Js7kffPjTHzQTGSPsW9S/CzYcvexuSLg0mslyrrZ474Z8PN2H7P3z+X47vH8WNdz00ExvNbMclc9fCNhFh5fPwpwRt29+Bmtb0p8vI7h6rFI3XYd135sYT4sEHiYALIlGEoQWGuClC94ftbpi7czyt1poxaps4UqDj1T1mXcf9C1Uy45uuTqi7UXX/BcYnSpc5BsPRnK2a8Xl55dEvSEKD8n9sxa/z8/Q65bqjKmn+n9RVcDT45ypmfZ+sLvCV4WiGVhCbvs6+Un0kQwr9+d3VZdPUmShe/nwHv67EOkbKtNOhPz8Vkw5ly3ONuzIF1ZqFj6ua/H0ZmnXw8atHjcnAk6+/fisfYxe1bvsLJ1VcQcq/GsHgtY4qNktt377T5dsfjvXbHWe1yOVz8vNjffmPUulrYrHJaxby5ZIPAwAWRLMJQgMMyrP76YWbYukValV7/D6QlQVRy6hOrP/FzUmY8v5J716+dkVTnV74A6W3LWOEi2npTk2b+Wrtv6OaqflTrzdEb5P9Ls3kq8zr5dGrWMnblZMC6fOgP+G0ku+fr1D5dWXeJ05ues1mVxZ9l+9ukLsRe3s5D/LN1M+FmvZezZtt4nL57c42ei3oafxfoCbvn7hsPtm8jWF3v/eeavhe9ns15chvXifbDgnYW7pOsKg29czLd6xmhWLt16vPzfLgG7LOqbAy/OMN2v/5wNVovHtsp8u2PKIveYW9ZfSpOcXKp1pWFY2TrDtfS9Px5PP9t1H4MmSLny4Jsd35y53O1ntBah98GZuEv0ZuSJVx4vP7O3JP3c2sfixePh8rJj0cfr2PMEMr8W5cXjbrbu09vx4wP3637qWHRz4b/njACynbMBDzrc7i+1cAbS/aUW1YXem/ZzQV/gPZnEpWCXcatXXPy5n1l6Moy/2MIXS89GtejczlmUS3S+KFpwg74hqnr1xzNLvZ7/7Qur97NaLHRn1RavS4TOatyvJwR5cSnQk538ao+fkVavNVlUfu5ocXnxRdgX8+q5dNW/xe6bCc+q9QXYF2Qfq28evLgPP6d01ucvXfArPhZK9aUJTWTrfXxLzRhX57lF4ufcnsntmw5njJ4k1l0GdvnT++znqp6oZaF0P3/sN/Gn3xeAeNt+rcdZqW9SLGvfELifYWXrQ/JjBN+ImJu/SMIl+0Gv/jiGnPn2LlVp2jdU/txf3uHF/+0bH49TtfjGyOPh15p8Q+R48Y1DtfiVKd/EOUYs6Go2cvdNSM0u8L/mkQCyncdR55ghAAEIQGCsBJDtWHGzMQhAAAIQmEcCyHYeR51jhgAEIACBsRJAtmPFzcYgAAEIQGAeCSDbeRx1jhkCEIAABMZKANmOFTcbgwAEIACBeSQwV7Jdv379kQcfrH6Kch6Hm2OGAAQgML0EOp3O1Dprane8SbgsLCwcWV4e9BWuTXplHQhAAAIQGAcBZDsOygW2gWwLQKQLCEAAAhMigGwnBD53s8g2lxjtIQABCLSHALJtz1gM3BNkOyUDxW5CAAIQqCGAbKckLJDtlAwUuwkBCEAA2U5vDCDb6R079hwCEIAAme2UxACynZKBYjchAAEIkNkOjAH/lJlfZK1eZn19+tkw/77j+yVtSD935t8t9c+leYn4rHYnkS3nLwQgAIHpJUBm+8jYWbaXSNrdNZx+l9e/Benfrrw5/Uaof7h6Z/rR8NKf9Y0kZDu9Jxl7DgEIQADZDpbtBenHns9PzfxD2Zayfxw84jNkyzkJAQhAYAYJINtvlO3XUsbqLPaNkp6RstiLu8b+AUlnSLoo4LP9/WKMzHYGzz4OCQIQmBsCyPaRod4s6U5J6yW9W9KipI8FCHWQpLtle5Uk/x1dOp3O5pWVlbkJTA4UAhCYXwJbt249evC33HLLzEBAtvVD6ZG+UdKLasrIt0s6sU8ZedTPKCPPzKnFgUAAAk0JINum5GLWK/lDBJbnoyRVZeS3STpX0qVpgtTlXROkzpK0Q9JxAZ8h25hYoVcIQGCKCCDbdg1WSdl6hvFHJS2kvz2SrpT0ZUkXpld/ju969eeehCLis1rKPLNtV/CxNxCAQBwBZBvHtknPJWXbZPtjXQfZjhU3G4MABCZIANlOEH7NppFtu8aDvYEABCBQhACyLYKxWCfIthhKOoIABCDQHgLItj1j4T1Btu0aD/YGAhCAQBECyLYIxmKdINtiKOkIAhCAQHsIINv2jAWZbbvGgr2BAASmnMChQ4d0+PDhVhzFtm3+jRdp165drdifdevWaf16f99R84UvtWjObqxrMht5rLjZGATmioBFe/rpW7R//91zddzDHuzJJ5+mu+7aO5Jwke2wtCfcDtlOeADYPARmmMCBAwe0ceNGSfvSN9VO+mBXM1upDZntAUmbtLS0pMVFf4tvswXZNuM29rWQ7diRs0EIzA2BR2S71BLZrn43stSG70a2bDci23k5G5DtvIw0xwmB8RNAtoOYI1tmI4//nGSLEIDADBJAtsh2EAFkO4MnPYcEAQiMnwCyRbbINhGgjDz+CxBbhMC8EEC2yBbZItt5ud5xnBCYGAFki2yRLbKd2AWIDUNgXgggW2SLbJHtvFzvOE4ITIwAskW2yBbZTuwCxIYhMC8E2ifbNpHn1R9mI7cpHtkXCEBgagkgWzJbMlsy26m9gE1qx2fxF1MmxXJetotskS2yRbbzcr0rdpzIthjKuekI2SJbZIts5+aCV+pAkW0pkvPTD7JFtsgW2c7PFa/QkSLbQiDnqBtki2yRLbKdo0temUNFtmU4zlMvyBbZIltkO0/XvCLHimyLYBzYyawxRrbIFtki2/gr54xtYdZE0MbhmTXGyBbZIltk28Zrbav3adZE0EbYs8YY2SJbZIts23itbfU+zZoI2gh71hgjW2SLbJFtG6+1rd6nWRNBG2HPGmNki2yRLbJt47W21fs0ayJoI+xZY4xskS2yRbZtvNa2ep9mTQRthD1rjJEtskW2yLaN19pW79OsiaCNsGeNMbJFtsgW2bbxWtvqfZo1EbQR9qwxRrbIFtki2zZea2v36dChQzp8+PDE93fbtm1H92HXrl0T35dqB9atW6f169e3Zn9G3RFkOyrBaVqf37Pl92ynKV5nfF8t2jPP3KJ77717xo+02eGdeuppuuOOvTMjXGTbLA6mc62xyPYtkt4qaauk3ZJ8x/x+SRsk7ZV0maR7E78dkq6StCDpLyVdLml/+qzfeuek/r5F0t9K+iNJPyXpwWHGBNkOQ4k2YyFQleGuu26fNmxYHMs2+23kta9dzWyvvrodme3Bgwe0Y8cmLS0taXFxsmxKDQyyLUVyGvoJl+1TJL1d0rmSLpb0BUm3SXq5pJslvUnSEyXtlPTtkn5P0ndKuk/SGyR9q6SfkGQn9ltvi6STJP2ZpOMkfUTSn6ftrjkIyHZNRDQYF4FKtjfcsKQTTpisUK64wjfH0nvec8u4Dn/gdh544IBe/OKNyLYVo1G/EzyzHTQ4obL1s5WbJL1E0mckXZKy2esknZ/26nEpu320pB9JEraUvVjUfyDJF50LJPVbr/cA/3kS948OE5ZzJVtJR84/v2I/DB7ajJPAysqK9uzZo82bz9Vxx/nGcXLLXXf55lY6/fSzJ7cTXVs2mzvv3KNzz508m1JAbrttlfHZZ7eD8ajHVcXvanI12fgd9VjKr78iafT43b17d52z3iHpDknvTUK1bL8tZbGVUH1ID0g6Q9IpKdt1+crlZa//2vT/LxqwXlVmdl8uTf9JyopvHIZXhGxz6uaDaupNPxt03Mh2mKiYUBtk2x88sp1QUGZsFtkOghUm2+9NZdxnOJkaUraW5oskvUaSd+xjkn5R0mMk/dAQsn2UpI8mwV8xbIiUlm1O3XxQbbzpZwOPe2Fh4cjy8vKwbGg3ZgKUkfsDp4w85mBssDnKyIOghZWRf0bSlZKqVxicud4j6RpJL+0pI98u6cSavfweSb+Rst66MnL3ep5Q9euSDkryJCsLfqilpGxz6+aDauNNP0O2Qw17OxshW2Tbzsgcbq+Q7URk27tRl4VdRt6TJjp5lnE1QeqsJEivc5okv/ZwfMpsPynp3an+7+cbdev52cD1aQazRf7QcJGx2qqkbHPr5oNq400/666pH8OBzDYnNMbfFtki2/FHXbktIttWydav/lyYXtWxUKtXf5z1evntNAN5naQbJP1sKin7s37r+fnvb6bXiirR/mHKrNcMpFKybVI3byrUYR9g++D9HpX/Vu8sOp3Nfq7C0k4CyBbZtjMyh9srZNsK2Q43WBNoVUq2Termg2rjTT8biHCWMttZe0fRA4ds42Xblm/o8pG27Vu6Rv2GLmSLbAcRKCXb3m0MUzd3/XtQbbzJZ8h2AndspTaJbGNla9FuOfNM3X1v9SU6pUZuNvo57dRTtfeOOxp/QxeyRbaTlO1adfNBNfWmn/U9XjLbdl8UkW2sbCu++669Vosb/JrgZJdtr3vd0R3Y9c53TnZHXFU5eFCbdu4c6UtDkC2ynYRsJ37y1O0Asm3lsDy8U8h2PLJd+shHWiHbrVf6jQ3plmv8lsZkF8t24/btyDZsGMJe/Qnb49IdR5WRS+9nkf6QbRGMYZ0gW2QbFlxrdIxso8kjW2QbHWNB/TNBKghsS7st8aUW1c0Mme2xg4xsowMf2SLb6BgL6h/ZBoFtabfINnZgkG0sXwnZItvoGAvqH9kGgW1pt8g2dmCQbSxfZFv2G6SiR2vk/nlmOzLC0A7a9Mw29EAbdI5sG0DLWAXZZsBq1JTMlsy2UeBMfiUy28mPwTj3ANnG0ka2sXzJbMlsoyMsrH9kG4a2lR0j29hhQbaxfJEtso2OsLD+kW0Y2lZ2jGxjhwXZxvJFtsg2OsLC+ke2YWhb2TGyjR0WZBvLF9ki2+gIC+sf2YahbWXHyDZ2WJBtLF9ki2yjIyysf2QbhraVHSPb2GFBtrF8kS2yjY6wsP6RbRjaVnaMbGOHBdnG8kW2yDY6wsL6R7ZhaFvZMbKNHRZkG8sX2SLb6AgL6x/ZhqFtZcfINnZYkG0sX2SLbKMjLKx/ZBuGtpUdI9vYYUG2sXyRLbKNjrCw/pFtGNpWdoxsY4cF2cbyRbbINjrCwvpHtmFoW9kxso0dFmQbyxfZItvoCAvrH9mGoW1lx8g2dliQbSxfZItsoyMsrH9kG4a2lR3PomzbBBrZRo8Gv/rDr/5Ex1hQ/8g2CGxLu0W2sQODbGP5ktmS2UZHWFj/yDYMbSs7Rraxw4JsY/kiW2QbHWFh/SPbMLSt7BjZxg4Lso3li2yRbXSEhfWPbMPQtrJjZBs7LMg2li+yRbbRERbWP7INQ9vKjpFt7LAg21i+yBbZRkdYWP/INgxtKztGtrHDgmxj+SJbZJsdYYcOHdLhw4ez1yu9wrZt2452uWvXrtJdN+5v3bp1Wr9+feP1Dxw4oI0bN+qGG5Z0wgmLjfuZxRWRbeyoIttYvsgW2WZFmEW75fTTdff+/VnrzUvj004+WXvvuquxcJFt/0hBtrFnEbKN5YtskW1WhFUy2Cdp0nnXal4rtSWv9SvrmyQtLS1pcbEZHWSLbLNOyIKNkW1BmLVd8aUWfKlFRoxVMlhqgWy3pv2+JWP/I5uunkrINooxmW0U2dV+kW0sXzJbMtusCEO2/XEh26xQym6MbLORZa2AbLNwNWhMZktmmxE2yBbZZoRL0abItijOYzpDtrF8yWzJbLMiDNki26yAKdgY2RaEWdMVso3li2yRbVaEIVtkmxUwBRsj24IwkW0szNreKSNTRs4IO2SLbDPCpWhTZFsUJ2XkWJx1tzNHp1CO8raCO+10OlPrrKnd8SaxsrCwcGR5ebnJqkfXQbbItnHwjLgish0R4BqrU0aO5UsZmTJyVoQhW2SbFTAFGyPbgjApI8fCpIxcS4DMNiPskC2yzQiXok2RbVGclJFjcVJGriFQWrYfl/QESSuSvi7pCkmfl3SOpOslnSLpfkkvk3Rb2p+Iz2pDiTJy3BnGe7ZxbN0zso3lSxk5li9l5PJl5MdI+loatudJerOkp0r6tKRrJX1Y0mWSdkp6emoX8RmyjT53evpHtrHAkW0sX2QbyxfZlpdt94i9VNKrJT1b0q0pq31I0kLKbs/25LKAz+7rFzZktnEnFLKNY0tmG8vWvSPbaMa8+lO6jOwR+4CkZyaRPkvS8ZI+KOm8ruHckzJcb7/0Z59DttEnzrH9I9tY5mS2sXyRbSxfMtv4zHZ7KiWXFuogSXfL9ipJ/ju6dDqdzSsrfpzcbGGCVH9uyLZZTA27FrIdllSzdsi2Gbfh1yKzjchsu/kflLRF0hcpIw8flsO05Fd/hqE0O22QbexYIttYvmS2ZTPbE5NQ70jD9lxJ75N0hqTPpPJyNUHqFZIuSu1uCvisNnJ4Zht3QpHZxrF1z8g2li+yjeWLbMvK9vGSPiHJ0vVEqP2SXiPJZd0npVd/TpL01fTqj7NdLxGfIdvoc6enf2QbCxzZxvJFtrF8kW1Z2UaP1sj9k9mOjLBvB8g2ji2ZbSxb945soxnzzDb6mW30CGb1j2yzcGU1RrZZuLIbk9lmI8taAdlm4WrQGNki24ywYTZyf1jINiOQGjRFtg2gZayCbDNgNWqKbJFtRuAgW2SbES5FmyLbojiP6QzZxvLlmS3PbLMiDNki26yAKdgY2RaEWdMVso3li2yRbVaEIVtkmxUwBRsj24IwkW0szNreKSNTRs4IO2SLbDPCpWhTZFsUJ2XkWJx1tzOSNmppaUmLi4uNt97pdKbWWVO7401Gi9nITagNtw4TpIbj1LQVsm1Kbrj1KCMPx6l5KzJbZJsRPWS2ZLYZ4VK0KbItipPMNhYnmW0NAWSbEXTIFtlmhEvRpsi2KE5kG4sT2SLbhSPLy8uNw6xNsm18EEErUkYOApu6RbaxfCkjx/JlNjKzkbMiDNmS2WYFTMHGyLYgzLq86+BBbdy+faQJPNX1QVqS1HwSUOyRTqp3ntlSRs6IPWSLbDPCpWhTZFsUJ2XkWJyUkSkjU0aOOscoI0eRXe0X2cbypYwcy5cyMmXkrAgjsyWzzQqYgo2RbUGYlJFjYdb2ThmZMnJG2CFbZJsRLkWbItuiOCkjx+KkjEwZmTJy1DlGGTmKLGXkWLKrvVNGjqZMZktmmxFjZLZkthnhUrQpmW1RnGS2sTjJbMlsyWyjzjEy2yiyZLaxZMlsx8GXCVJMkMqKMzJbMtusgCnYmMy2IMy6vIv3bGMBizIyZeSMEEO2yDYjXIo2RbZFcVJGjsVJGZkyMmXkqHOMMnIUWcrIsWQpI4+DL2VkyshZcUZmS2abFTAFG5PZFoRJGTkWZm3vlJEpI2eEHbJFthnhUrQpsi2KkzJyLE7KyJSRKSNHnWOUkaPIUkaOJUsZeRx8KSNTRs6KMzJbMtusgCnYmMy2IEzKyLEwKSPXEqCMnBF2yBbZZoRL0abItihOysixOCkjU0amjBx1jlFGjiJLGTmWLGXkcfCljEwZOSvOyGzJbLMCpmBjMtuCMCkjx8KkjEwZeWGBzDbqLCOzjSJLZhtLlsx2HHzJbMlss+KMzJbMNitgCjYmsy0Ik8w2FiaZLZktmW3cOUZmG8fWPSPbWL78xF4sXzJbMtusCCOzJbPNCpiCjZFtQZhktrEwyWzJbMls484xMts4tmS2sWzdO5ltNGO+rpH3bDNijMyWzDYjXIo2JbMtivOYzpBtLF/KyJSRsyIM2SLbrIAp2BjZFoRJGTkWJmVkysiUkePOMcrIcWwpI8eypYwcz5fMtmxme4qkD0k6S9IhSbdKeqWk+yVtk/R+SRsk7ZV0maR70xBHfFYbPcg27qRCtnFskW0sW2QbzxfZlpXtyZKeLOmmNHRXS3qspB+XdJukl0u6WdKbJD1R0k6tbr/0Z30jB9nGnVTINo4tso1li2zj+SLbsrLtHbFLJb1K0hskXSfp/NTgcSm7fbSkCwI+Q7bjOHd6toFsY6HzzDaWLxOkYvki2zjZHifpU5JulLQvZbEXdw3nA5LOkHRRwGf7+4UNmW3cCYVs49iS2cayJbON54ts42T7XkmbJL1A0vMDhDpI0t2yvUqS/44unU5n88rKSuPIYjZyf3TItnFYDbUime1QmBo3IrNtjG7IFXnPNuV+ViQAAA+0SURBVOI923elZ7fPSROl6krFt0s6sU8ZedTPKCMPGf4lmyHbkjSP7QvZxvJFtrF8yWzLZ7Zvl3ShJJeMXSr24pKyJ0Fd3jVByjOWdwR9hmyjz5ua/pFtLHRkG8sX2cbyRbZlZXuepN3plZ+Daei+JOmFScB+9ef4rld/7kltLOfSn9VGDs9s404oZBvH1j0j21i+yDaWL7ItK9vo0Rq5f2Q7MsK+HSDbOLbINpate0e20Yx5ZhvxzDZ61Br3j2wbo1tzRWS7JqKRGpDZjoRvzZWR7ZqIRmyAbJFtRggxG7k/LGSbEUgNmiLbBtAyVkG2GbAaNUW2yDYjcJAtss0Il6JNkW1RnMd0hmxj+fLMlme2WRGGbJFtVsAUbIxsC8Ks6QrZxvJFtsg2K8KQLbLNCpiCjZFtQZjINhZmbe+UkSkjZ4QdskW2GeFStCmyLYqTMnIszrrbGUkbtbS0pMXFxcZb73Q6U+usqd3xJqPFbOQm1IZbhwlSw3Fq2grZNiU33HqUkYfj1LwVmS2yzYgeMlsy24xwKdoU2RbFSWYbi5PMtoYAss0IOmSLbDPCpWhTZFsUJ7KNxYlske3CkeXl5cZhhmyRbePgGXFFZDsiwDVWp4wcy5fZyMxGzoowZItsswKmYGNkWxBmXd518KA2bt8+0gSe6vogLUlqPgko9kgn1TvPbCkjZ8QeskW2GeFStCmyLYqTMnIsTsrIlJEpI0edY8xGjiK72i+yjeVLGTmWL2VkyshZEUZmS2abFTAFGyPbgjApI8fCrO2dMjJl5IywQ7bINiNcijZFtkVxUkaOxUkZmTIyZeSoc4wychRZysixZFd7p4wcTZnMlsw2I8bIbMlsM8KlaFMy26I4yWxjcZLZktmS2UadY2S2UWTJbGPJktmOgy8TpJgglRVnZLZktlkBU7AxmW1BmHV5F+/ZxgIWZWTKyBkhhmyRbUa4FG2KbIvipIwci5MyMmVkyshR5xhl5CiylJFjyVJGHgdfysiUkbPijMyWzDYrYAo2JrMtCJMycizM2t4pI1NGzgg7ZItsM8KlaFNkWxQnZeRYnJSRKSNTRo46xygjR5GljBxLljLyOPhSRqaMnBVnZLZktlkBU7AxmW1BmJSRY2FSRq4lQBk5I+yQLbLNCJeiTZFtUZyUkWNxUkamjEwZOeoco4wcRZYycixZysjj4EsZmTJyVpyR2ZLZZgVMwcZktgVhUkaOhUkZmTLywgKZbdRZRmYbRZbMNpYsme04+JLZktlmxRmZLZltVsAUbExmWxAmmW0sTDJbMlsy27hzjMw2jq17RraxfPmJvVi+ZLZktlkRRmZLZpsVMAUbI9uCMMlsY2GS2ZLZktnGnWNktnFsyWxj2bp3MttoxnxdI+/ZZsQYmS2ZbUa4FG1KZlsU5zGdIdtYvpSRKSNnRRiyRbZZAVOwMbItCJMycixMysiUkSkjx51jlJHj2FJGjmVLGTmeL5ktmW1WlJHZktlmBUzBxmS2BWGS2cbCJLMNz2zfJelSSVskbZW0O23xHEnXSzpF0v2SXibptsDP+gYSmW3cOUZmG8eWzDaWLZltPF8y27KZ7dMk3S7ps5Iu6ZLtpyVdK+nDki6TtFPS09PwRnyGbMdx7vRsA9nGQiezjeXLBKlYvsi2rGyr0drbJdtTJd2astqHJC2k7PZsrW679Gf3DQoZMtu4EwrZxrEls41lS2YbzxfZxsv2qZI+KOm8ruHckzJcy7b0Z59DtuM4cY7dBrKN5U5mG8uXzDaWL7KdfdleJcl/R5dOp7N5ZWWlcVQxQao/OmTbOKyGWhHZDoWpcSNk2xjdkCvypRYRX2pBGXnI8JulZsg2djSRbSxfZBvLl8w2PrP1CN4k6QNdE6ReIemiNLQRn/WNGp7Zxp1QyDaOrXtGtrF8kW0sX2RbVrbXSHqBpNMkfSVNhPKz2ielV39OkvTV9OrPF9PQRnyGbKPPm5r+kW0sdGQbyxfZxvJFtmVlGz1aI/dPZjsywr4dINs4tmS2sWzdO7KNZswz24hnttGj1rh/ZNsY3ZorIts1EY3UgMx2JHxrroxs10Q0YgNki2wzQojZyP1hIduMQGrQFNk2gJaxCrLNgNWoKbJFthmBg2yRbUa4FG2KbIviPKYzZBvLl2e2PLPNijBki2yzAqZgY2RbEGZNV8g2li+yRbZZEYZskW1WwBRsjGwLwkS2sTBre6eMTBk5I+yQLbLNCJeiTZFtUZyUkWNx1t3OSNqopaUlLS4uNt56p9OZWmdN7Y43GS1mIzehNtw6TJAajlPTVsi2Kbnh1qOMPByn5q3IbJFtRvSQ2ZLZZoRL0abItihOMttYnGS2NQSQbUbQIVtkmxEuRZsi26I4kW0sTmSLbBeOLC8vNw4zZItsGwfPiCsi2xEBrrE6ZeRYvsxGZjZyVoQhW2SbFTAFGyPbgjDr8q6DB7Vx+/aRJvBU1wdpSVLzSUCxRzqp3nlmSxk5I/aQLbLNCJeiTZFtUZyUkWNxUkamjEwZOeocYzZyFNnVfpFtLF/KyLF8KSNTRs6KMDJbMtusgCnYGNkWhEkZORZmbe+UkSkjZ4QdskW2GeFStCmyLYqTMnIsTsrIlJEpI0edY5SRo8hSRo4lu9o7ZeRoymS2ZLYZMUZmS2abES5Fm5LZFsVJZhuLk8yWzJbMNuocI7ONIktmG0uWzHYcfJkgxQSprDgjsyWzzQqYgo3JbAvCrMu7eM82FrAoI1NGzggxZItsM8KlaFNkWxQnZeRYnJSRKSNTRo46xygjR5GljBxLljLyOPhSRqaMnBVnZLZktlkBU7AxmW1BmJSRY2HW9k4ZmTJyRtghW2SbES5FmyLbojgpI8fipIxMGZkyctQ5Rhk5iixl5FiylJHHwZcyMmXkrDgjsyWzzQqYgo3JbAvCpIwcC5Myci0BysgZYYdskW1GuBRtimyL4qSMHIuTMjJlZMrIUecYZeQospSRY8lSRh4HX8rIlJGz4ozMlsw2K2AKNiazLQiTMnIsTMrIlJEXFshso84yMtsosmS2sWTJbMfBl8yWzDYrzshsyWyzAqZgYzLbgjDJbGNhktmS2ZLZxp1jZLZxbN0zso3ly0/sxfIlsyWzzYowMlsy26yAKdgY2RaESWYbC5PMlsyWzDbuHCOzjWNLZhvL1r2T2UYz5usaec82I8bIbMlsM8KlaFMy26I4j+kM2cbypYxMGTkrwpAtss0KmIKNkW1BmJSRY2FSRqaMTBk57hyjjBzHljJyLFvKyPF8yWynO7M9R9L1kk6RdL+kl0m6bVDYINu4kwrZxrFFtrFskW08X2Q73bL9tKRrJX1Y0mWSdkp6OrIdx4lz7DaQbSx3ysixfHlmG8sX2U6vbE+VdGvKah+StJCy27Ml3dcvbMhs404oZBvHlsw2li2ZbTxfZDu9sn2qpA9KOq8rTPakDPdz0bLdJ2lxHPE5RduwbDdJWlpa0uJiMzrVBLTrrtunDRua9TFFyLJ29eDBA9qxY1MRvvuuvVaLGzZkbX/WGzuz3bRzZxG+EleI+trXaPHrPjudztS+QTOtOz6sbK+S5L9q2SzJmfCsLMdJWpmVg2nhccA3flBgHMt41vguSzo+FllM79Mq20Zl5BiEE+31DklnTnQPZnvj8I0fXxjHMoZvLN+he59W2foAb5L0ga4JUq+QdNHQRz4bDTmRYscRvrF83TuMYxnDN5bv0L1Ps2yflF79OUnSV9OrP18c+shnoyEnUuw4wjeWL7KFbzyBlmxhmmXbEoQT3Q0/j/7lie7BbG8cvvHjC+NYxvCN5Tt078h2aFQ0hAAEIAABCDQjgGybcWMtCEAAAhCAwNAEkO3QqGgIAQhAAAIQaEYA2TbjFrXWuyRdKmmLpK2SdqcNXSzpFyR9k6T9kl4u6a/W+MzfGf0hSWdJOpS+ceuV6Zu2ova/7f2W5Otj/bikJ6R3nb8u6QpJn287hMD9K8232tW3SHprzzkReBit7bo0372SHkx/PujXS/pka49+yncM2bZrAJ8m6XZJn5V0SZKtZ1v7BxYuTMLcLulySc+SNOizkyU9Ob0i5aO8WtJjJf1Yuw55rHtTkq93/DGSvpaO4HmS3izJX7gyr0tpvub4FElvl3SuJN90Vjeg88i4NF/LtrrOzCPPsR4zsh0r7qE31n0SXCDp17q+mtIS9a8cPT5lVf0++0rP1pwxv0rSDw69F7PbMILvSyW9WtJ3zy62oY+sFN/16WbxJZI+gxge5l+KL7IdOqRHb4hsR2cY0UP3SeDs6S8k/bCkP5b0k5LemzIo//9+n3V/R7S/su1Tkm6UdE3EDk9ZnyX5+otVnumvbU3Vhi9MGYuI3S3F9x3pSy8c74jhkZEqxdf9uDLj2L1Z0hsl/XVEQNDn9P4QwayPXe+FxT8d+HPpO0F/K2VQ3y/plvSzgv0+qzj5YuXfCXgB36V8FElpvu7Tma1L/M+e9eAc4vhK8H10Kh8/Q9IRZPsN1Evw9bXD3xV/pyRXEN6dfl/FP1fKEkCAzDYAaoEuB93Ff3N6rusJUH/Ts626zzypws9un5MmShXYvanvoiTfbhgHJZ0x55PQ6m5muhkNG7+ebHalpMNpZXO9R9KOVKWZ+iAc4QAi4tcTMl358oQ/lgACyDYAaoEue0+m0yTdLcnl4F/V6o9D/nTazqDPPLHEE6s8seSBAvs1K12U4Hti+j1lf6Wjl+dKel+SrTOxeV5K8O3lRxn5ESIl+Dp+H9VVRn5bmoT2/HkO3MhjR7aRdPP79vNUl3otUE9w8kQo/2avnwt+n6R16a7eX8HmKfte+n3m9Txz81ZJzri8fEnSC/N3a2bWKMnXE9Q+IckXLf9so1/Jeo2kvr+nPDMU+x9ISb7I9ljOJfk+UdJHJS2kP/8euCsJX56DOJ3IISLbiWBnoxCAAAQgME8EkO08jTbHCgEIQAACEyGAbCeCnY1CAAIQgMA8EUC28zTaHCsEIAABCEyEALKdCHY2CgEIQAAC80QA2c7TaHOsEIAABCAwEQLIdiLY2SgEIAABCMwTAWQ7T6PNsUIAAhCAwEQIINuJYGejEIAABCAwTwSQ7TyNNscKAQhAAAITIYBsJ4KdjUIAAhCAwDwRQLbzNNocKwQgAAEITIQAsp0IdjYKAQhAAALzRADZztNoc6wQgAAEIDARAsh2ItjZKAQgAAEIzBMBZDtPo82xQgACEIDARAgg24lgZ6MQgAAEIDBPBJDtPI02xwoBCEAAAhMhgGwngp2NQgACEIDAPBFAtvM02hwrBCAAAQhMhACynQh2NgoBCEAAAvNEANnO02hzrBCAAAQgMBECyHYi2NkoBCAAAQjMEwFkO0+jzbFCAAIQgMBECPwduZZxlYIxI5MAAAAASUVORK5CYII=\" width=\"431.818172458775\">"
+      ],
+      "text/plain": [
+       "<IPython.core.display.HTML object>"
+      ]
+     },
+     "metadata": {},
+     "output_type": "display_data"
+    }
+   ],
+   "source": [
+    "#Comment 22/08/2002: could net use the \"plot\" function because cannot control/update the artist with the event handler.\n",
     "\n",
-    "# confidence_interval_1 = stats.norm.interval(0.95, loc=32000, scale=2000)\n",
-    "# confidence_interval_2 = stats.norm.interval(0.95, loc=43000, scale=1000)\n",
-    "# confidence_interval_3 = stats.norm.interval(0.95, loc=43500, scale=1400)\n",
-    "# confidence_interval_4 = stats.norm.interval(0.95, loc=48000, scale=700)\n",
+    "#Initialize figure and axes:\n",
+    "fig, ax = plt.subplots()\n",
     "\n",
-    "# print(confidence_interval_1)\n",
-    "# df\n",
+    "#Initializate threshold\n",
+    "threshold = 41000\n",
     "\n",
-    "from scipy import stats\n",
-    "year_avg = df.mean(axis = 1)\n",
-    "year_std = df.std(axis = 1)\n",
-    "yerr = year_std / np.sqrt(df.shape[1]) * stats.t.ppf(0.95, df.shape[1]-1)\n",
+    "# #bar plot           \n",
+    "bar_plot = ax.bar(df.index, year_avg, align='center',yerr=yerr,color = colour_list(threshold), edgecolor =\"black\")\n",
     "\n",
-    "plt.figure()\n",
-    "plt.show()\n",
+    "#set line and thereshold in the graph:\n",
+    "ax2 = ax.twinx() # \n",
+    "ax2.set_ylim(ax.get_ylim()) # set the y axis of ax = to y axis of ax2\n",
+    "ax2.tick_params(axis= 'both', which= 'both',length=0) # remove the ticket length\n",
+    "lhor = ax2.axhline(y = threshold, zorder=1, color = 'black') # set the horizontal line\n",
+    "ytick = plt.gca().get_yticks()\n",
+    "ytick = []\n",
+    "ytick = np.append(ytick,threshold)\n",
+    "plt.gca().set_yticks(ytick)\n",
     "\n",
-    "bars = plt.bar(range(df.shape[0]), year_avg, yerr= yerr, alpha = 0.6, color = 'rgby')\n",
-    "# bars = plt.bar(range(df.shape[0]), year_avg, color = 'rgby')\n",
-    "# plt.errorbar([0,1,2,3], list(year_avg), yerr=stats.t.ppf(0.95, )*list())\n",
+    "# Add xticks, titles and remove left and top edges:\n",
+    "plt.xticks(df.index, ('1992', '1993', '1994', '1995'))\n",
+    "plt.title('Data between 1992 and 1995')\n",
+    "ax_list = fig.axes # list of axis\n",
+    "for i in ax_list:\n",
+    "    [i.spines[loc].set_visible(False) for loc in ['top', 'right']]\n",
     "\n",
-    "threshold=42000\n",
-    "plt.axhline(y = threshold, color = 'grey', alpha = 1)\n",
+    "#event handler\n",
+    "def onclick(event):\n",
+    "    threshold = event.ydata#get threshold based on the event\n",
+    "    lhor.set_ydata(threshold)#update horizontal line\n",
+    "    ytick = ax2.get_yticks()#set ax2 for new thereshold\n",
+    "    ytick = []#set ax2 for new thereshold\n",
+    "    ytick = np.append(ytick,threshold)#set ax2 for new thereshold\n",
+    "    ax2.set_yticks(ytick) #update ax2\n",
+    "    for ind,bar in enumerate(bar_plot):\n",
+    "        bar.set_facecolor(colour_list(threshold)[ind])#update colours for each bar\n",
     "\n",
-    "plt.xticks(range(df.shape[0]), ['1992', '1993', '1994', '1995'], alpha = 0.8)\n",
-    "plt.title('Ferreira et al, 2014')\n",
+    "cid = fig.canvas.mpl_connect('button_press_event', onclick)\n",
     "\n",
-    "# plt.show()"
+    "#show figure:\n",
+    "plt.show()"
+   ]
+  },
+  {
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "**Hardest option implementation:** Allow the user to interactively set a range of y values they are interested in, and recolor based on this (e.g. a y-axis band, see the paper for more details)."
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "#not implemented"
    ]
   }
  ],
@@ -145,7 +2642,7 @@
    "name": "python",
    "nbconvert_exporter": "python",
    "pygments_lexer": "ipython3",
-   "version": "3.5.2"
+   "version": "3.8.3"
   }
  },
  "nbformat": 4,