{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/Users/mitchellleibowitz/anaconda/lib/python3.5/site-packages/matplotlib/figure.py:1744: UserWarning: This figure includes Axes that are not compatible with tight_layout, so its results might be incorrect.\n",
      "  warnings.warn(\"This figure includes Axes that are not \"\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAagAAAEYCAYAAAAJeGK1AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XecHWXZ//HPlWwaKZsQehCkiCIkbER8RIhZRPhJsQAW\nVBALTYwU4Sf1UTQIqIA0pT4YRFCagiBFETaQB0EgrHQISACRkrabnpDd6/njnsM5e3La7pk5M5t8\n36/XvKbPXHty7nPlvueeGXN3REREsmZA2gGIiIiUogQlIiKZpAQlIiKZpAQlIiKZpAQlIiKZpAQl\nIiKZlHiCMrNNzKzDzI4use6bZtZdZngw6dhERCS7mpI8uJkNB/4AjCyzyQ6AA2cDK4rW/TvB0ERE\nJOMSS1BmtjkhOU0kJKFSJgDz3f3UpOIQEZH+KZEmPjM7FngCGA/8rcKm44Enk4hBRET6t6SuQR0D\nvAxMAn4LWPEGZjYOWJeQyERERHpIKkEdDkx094crbDMhGg82sz+a2VtmttDM7jKznRKKS0RE+olE\nEpS7/9WrP4U2l6COAIYAVwF/AT4BPGBmeyQRm4iI9A+J9uKrYgAwGzjF3X+fW2hmk4B7gV+b2Zbu\nvjKl+EREJEWpJSh3Pws4q8TyB8zsWuBgYDLw11qOZ2aPAx8EllXYrL0PoYqISNBSYd0w4Bl3nxjX\nydKsQVUyk5CgtujFPh8EBjc3Nw8ut0FLS8vkegOLS3t7yJUtLZX+vbND8SZL8Savv8WcxXhzMZXS\n2dkJ4Xc4NqklKDObCIxw9wdKrB4WjZf34pDLmpubB3d0dNQfXAO0trYC0NbWlmoctVK8yVK8yetv\nMfe3eEePHk1nZ2elFqxeS7MGdSuwsZlt6O7zi9ZNisaPNjgmERHJiDQfFntDdP4zCxea2ReAvYHp\n7v5MGoGJiEj60qxBTQX2Ag4zsx2AGcAHCMnpdeCbjQ6oqwvuvRcmTIANN2z02UVEpFCjalCr3RPl\n7p3AzsD5wEbAdwnP7bsC+LC7z25QbO866STYc08YPx6WxdqSKiIivZV4DcrdrwauLrNuIXB8NKRu\nxowwnjMH2tth553TjUdEZG2mFxYWGDo0P60alIhIuqz6E4n6BzNrmzx58uR6umTusw/ccUeYvv32\nMC8iItW1trYyffr06e7eGtcxVYMqsM46+emlS9OLQ0RElKB6GDYsP60EJSKSLiWoAqpBiYhkhxJU\ngcIEpU4SIiLpUoIqoBqUiEh2KEEV0DUoEZHsUIIqoCY+EZHsUIIqoCY+EZHsUIIqoCY+EZHsUIIq\noCY+EZHsUIIqoCY+EZHsUIIqoCY+EZHsUIIqoBqUiEh2KEEV0DUoEZHsUIIqoBqUiEh2KEEV0DUo\nEZHsUIIqoCY+EZHsUIIqoCY+EZHsUIIqMGgQDIg+kXfeCYOIiKRDCaqAmZr5RESyIvEEZWabmFmH\nmR1dZv3XzGymmS02s9fM7FwzG550XOUoQYmIZEOiCSpKNH8ARpZZfzIwDTDgQqAdOA6428yakoyt\nHPXkExHJhsQSlJltDtwPfKTM+s2AHwH/C3zY3U9x908DU4GdgcOTiq0SdZQQEcmGRBKUmR0LPAGM\nB/5WZrMjgIHAme7eVbD8TGARcGgSsVWjJj4RkWxIqgZ1DPAyMAn4LaEJr9ikaDy9cKG7rwD+Duxg\nZiWbBpOkGpSISDYklaAOBya6+8MVttkKeMvdS6WB2dF4m7gDq0bXoEREsiGRBOXuf3V3r7LZWKCj\nzLrOaNwcX1S1UQ1KRCQb0rwPahCwosy63PKhDYrlXboGJSKSDal05Y4sAwaXWTckGi/pzQHb29tp\nbW0tu76tra3qMdTEJyJSWqXf1/b29tjPl2YNagHlm/ByyzvLrE+MmvhERLIhzRrUC8DHzWxI1HOv\n0BZANzCrNwdsaWmpqZZUiZr4RERKq/T72trayvTp08uu74s0a1AzovNPKlxoZkOAjwJPuXuvmvji\noBqUiEg2pJmgriPUkk43s8JrUacSHo10eRpB6RqUiEg2pNbE5+7Pm9k5wPeBx83sNmB7YG/gAeDK\nNOJSE5+ISDY0qgZV8p4odz8ZmEKoSR0NfBA4F9jX3VN5G5Oa+EREsiHxGpS7Xw1cXWH9JcAlScdR\nKzXxiYhkg15YWEQ1KBGRbFCCKqJrUCIi2aAEVURNfCIi2aAEVURNfCIi2aAEVURNfCIi2aAEVUQ1\nKBGRbFCCKqJrUCIi2aAEVUQ1KBGRbFCCKjK04BWJy5dDd3d6sYiIrM2UoIoMGLB6khIRkcZTgipB\nzXwiIulTgipBXc1FRNKnBFWCevKJiKRPCaoENfGJiKRPCaoEJSgRkfQpQZWga1AiIulTgipB16BE\nRNKnBFWCmvhERNKnBFWCmvhERNKnBFWCmvhERNKnBFWCmvhERNKnBFWCmvhERNKXiQRlZlPNrLvM\ncF2j41ENSkQkfU1pBxCZACwHzgKsaN1TjQ5G16BERNKXpQT1jLtPTTsQUA1KRCQLUm/iM7ORwObA\nE2nHkqNrUCIi6Us9QRFqT5ChBKUmPhGR9GUlQTmwgZn9xczmR8ONZrZNGgGpiU9EJH1ZSVAGHA90\nApcDDwH7Aw+b2YQK+yZCTXwiIunLQieJLmA2cIi7P5BbaGZfBq4FrgI+3MiA1MQnIpK+1BOUu08p\ns/x3ZnYEMMnM3ufus6odq729ndbW1rLr29raaopJTXwiIqur9Pva3t4e+/my0MRXycxovEUjT6oE\nJSKSvlRrUGY2EJgIDHD3f5TYJNfYtryW47W0tNRcS6pE16BERFZX6fe1tbWV6dOnx3q+tJv4BgIP\nAgvNbH1396L1uwCrgPjrjhXoGpSISPpSbeJz95XAn4AxwEmF68zsBGB74Fp3X9jIuIoT1GppU0RE\nEpd2DQpC9/KPAWeY2W7AP4EdgVbCc/iOb3RAgwaF4Z13oLs7jAcPbnQUIiJrt9Q7Sbj7K4SEdBWw\nHfBd4L3Az4Fd3H1BGnGpmU9EJF1ZqEHh7m8Ah6UdR6F11oGFUcPi0qUwenS68YiIrG1Sr0Fllbqa\ni4ikSwmqjMImPnU1FxFpPCWoMlSDEhFJlxJUGUpQIiLpUoIqQ0+TEBFJlxJUGepmLiKSLiWoMtTE\nJyKSLiWoMtTEJyKSLiWoMtTEJyKSLiWoMtTEJyKSLiWoMpSgRETSpQRVhq5BiYikSwmqjOHD89OL\nFqUXh4jI2koJqozCp5d3dKQXh4jI2koJqowxY/LTC1J5I5WIyNpNCaqMwgQ1f356cYiIrK2UoMpQ\nDUpEJF1KUGUoQYmIpEsJqozCThKdndDVlV4sIiJrIyWoMgYOhObm/HxnZ3qxiIisjZSgKlAzn4hI\nepSgKlBPPhGR9GQiQZnZQDM7zsyeNrOlZvaSmZ1mZk1pxqUalIhIejKRoIBfAecCc4DzgX8DPwau\nSzMoJSgRkfSkWkMBMLOPAYcBN7j7gQXLpwEHm9ne7n5HGrGtu25+WglKRKSxslCD+g7gwI+Klp8c\njQ9tbDh5qkGJiKQnCwlqEjDX3Z8tXOjubwAvAJNTiQolKBGRNKWaoMxsMLAp8FKZTWYDo81sbMOC\nKqBefCIi6Um7BpW7ylPuhRa522Oby6xPlGpQIiLpSbuTxKBovKLM+tzyoQ2IZTVKUCKrc4fu7vD4\nr0pDbpvu7srTvRly5+7NuHi6cCi1PLcs97dWGoq3KZzPTZc6TrXtCpcVfu6Vtq+0XbnjFU9X26/4\nGIVeKtcOVoe0E1TuZeqDy6wfEo2X1HKw9vZ2Wltby65va2urNS5ACUoao7sbli+HZct6DsuX9xxW\nrMiPV64M48LplSvzwzvv9JwuNaxaFYbC6eKhq2v1sZ5LuTZrrbCuPfazpZ2gOoFuyjfhNRds13Dq\nZi6lrFoVns3Y0RHGnZ2wcGHPYdGiMCxenB+WLMkPS5fmx8uXp/0XiWRTqgnK3d8xs1eALcpssgUw\nx91reul6S0tLr2tJlagGtebr7oZ58+Ctt+Dtt8MwZ04Y5s4Nw/z5YZg3L3wPFi1KO+r0mYUHKlcb\nBgwIQ2564MCe++amc9sVzhdPF25TPF9ueeGywvlal+f+1nJDqfWFy3LTfdmueFnhuNz2hesq7VNq\nvpZzQdtqx8i56KJWZs2aXnplH6VdgwKYARxkZlu7+4u5hWa2MbANcGtagTU3h38Y9/C/4lWroCkL\nn5jUZOVK+Pe/4dVXw/Daa/D662H4z3/gjTdCYlq1Ku1IYehQGDas5zBkSBgPHRqmhw7NT+eGwYPD\nkJseNKjn9KBBPadzQ1NTz+ncfFNTSAqDBuWTSG594Xy5HylZe918M8yaFe8xs/Bz+xvgYOBMM/uS\n+7uX4c4m3MB7eVqBDRgQklRHVH/r6ID11ksrGill0SJ44YVQMGbNChdq//UvePnlkIgqXdStR3Nz\neGfY6NFhOjeMHAmjRoXxyJEwYkR+GD6857DOOmEYOjR810Skp9QTlLv/zcyuB74I/N3M7gN2iYYb\n3f3ONOMbMyafoBYsUIJKy9Kl8PTT8MQT8OST8Mwz8OyzoYZUr9GjYcMNw7DBBrD++mFYbz0YOzaM\n1103DGPGhAQ0cGD95xWRylJPUJGDgKeArwPHAK8C/w38PMWYgPCD9PLLYVrXoRpj2TJ47DF45JEw\nfuwxeP753teGzGDcONhsszBsumkYxo2DTTYJw0YbhRqMiGRPJhKUu3cBP4mGTFFPvuTNmQP33x+G\nBx+E9vbarws1NcHWW8P73gfbbBOmt9wyDJttFq6/iEj/lIkElWXqyRe/JUugrQ3uvhvuuw+eeqr6\nPgMGhAQ0YQKMHw/bbQfbbgtbbRUu6IvImkcJqgo9jy8er70Gt9wCf/pTqCmtXFl5+223hf/6L9hx\nxzDssEPoUCAiaw8lqCpUg+q7V16B3/0ObropXEcqp6kpJKPJk2HXXeGjH+35uYvI2kkJqgolqN7p\n7AxJ6dprYcaM8ttttx186lOw556wyy6h27WISCElqCqUoKpzh4cegssvh+uvD73wijU1wSc+Afvt\nB/vuG3rTiYhUogRVhXrxlbdyJdxwA5x3Hjz++OrrBw6EPfaAL38ZPvOZcL+RiEitlKCqUA1qdYsX\nwyWXwPnnh0cGFZswAQ49FL70pXDjq4hIXyhBVaFefHmLF8OvfgU//3l4iGqhYcPgK1+Bww+HnXbS\ns9pEpH5KUFWoBhXeF3TZZfDjH4ebagtttBFMmQJHHKHHQIlIvJSgqlibE5Q73HEHnHACPPdcz3Wb\nbw6nngpf+1p4eraISNyUoKoYNSr/yo3Fi0NtYm14csHs2XDUUXBn0aN6N9sMTjsNDjlEjxESkWTp\nIf9VDBjQsxbVUdOrE/uvVatCr7zttuuZnEaOhJ/+NDy09bDDlJxEJHmqQdVgzJh8B4kFC8KrGNZE\nzz0HBx8Mjz6aX2YWEtLUqeqRJyKNpQRVgzW9J597uMn2uON63mS7/fZwxRXh0UMiIo2mJr4arMkd\nJebNg/33hyOPzCenIUPgJz+BmTOVnEQkPapB1WBNTVD//Cd87nOhQ0TOdtuFZ+mNH59aWCIigGpQ\nNVkTE9T118POO/dMTlOmhLfYKjmJSBYoQdVgTUpQ3d1w8slw4IH5Jr2RI+GPf4SLLgpPhBARyQI1\n8dVgTXlg7MqV8I1vwHXX5Ze9731w663hBYEiIlmiBFWDNaEGtXgxfP7z4TXrOXvtFZKVnjIuIlmk\nJr4a9Pdu5nPnhncxFSanI46A225TchKR7FKCqkF/rkHNmwe77x46P+T88IfhdRkDB6YXl4hINakm\nKDN7zcy6SwxdZrZnmrEV6q8Jav58+OQn4YknwrxZeF3G6afrdRgikn2pXYMyszHAOOAh4K4Sm7zY\n2IjK648JasEC2HNPaG8P82YwbVp4+riISH+QZieJCdH4One/OMU4qirsxTd/fng0UJZrIEuWhA4Q\njz2WX/Y//6PkJCL9S5pNfBMAB55IMYaajBwJw4eH6aVLs12LWrUq3OP08MP5ZVdcEbqXi4j0J2kn\nKOgHCcoM3vve/Hzh0xeyxD08DeL22/PLLroIDj00vZhERPoq7QQ1DzjMzJ4ys6Vm9pKZnW5mmXvb\nUH9IUGefHV7NnnPiiSFhiYj0R6kkKDMzYDtgLHAscB9wJfAO8APgdjPLVBf4rCeoG2+EU07Jz3/l\nK3DmmenFIyJSr1g7SZjZbGCzKptdDJwBvADMB/Zz90XR/oOBm4B9gKOibTOhMEG9/HJqYZT05JPw\n9a/n53fbDa66KrwNWESkv4q7F9/NQLX3zf7D3d8GPlS8wt1XmtnRwL7Al+llgmpvb6e1tbXs+ra2\ntt4croes1qAWLID99gudNwC23hpuvjm800lEJE6Vfl/bc/e0xCjWBOXux8dwjNlmtgDYIoaQYrNF\nQTRZSVBdXfDVr8JLL4X54cPhllt63rclItJfpXIflJmtD7wfmO3u/y6xyVCgs7fHbWlpqauWVElx\nDSoL90L96Edw5535+auvDi8cFBFJQqXf19bWVqZPnx7r+dK6SrEvcD9wQvEKM/swMAx4pHhdmtZd\nF0aMCNOLF6f/0Nj77oMzzsjPn3wyHHBAevGIiMQtrQR1O7AM+IaZbZNbaGYjgQsIN/D+MqXYSiq+\nFyrNjhLz5sHBB4daHIQnlU+dml48IiJJSCVBufsc4HvACOBRM7vCzC4GngQ+Cpzt7venEVslWego\n4Q7f+ha8/nqYHzsWrrlGTyYXkTVPas/ic/fLom7pJwJfAoyQoE509+vTiquSLHSUuPTS8AbcnGnT\nYJNN0olFRCRJqb5R193vBu6uumFGpF2Deu45+N738vPf/S7su2/j4xARaQTdytkLaSaorq7wwNfl\ny8P8hAnws581NgYRkUZSguqFNDtJ/OIX8NBDYXrQoHDdaejQxsYgItJISlC9UOpeqEZ47jk47bT8\n/H//d6hBiYisyZSgemHMGBg1KkwvXQpz5yZ/zlzT3ooVYf5DH4KTTkr+vCIiaVOC6oU03gt1wQU9\nm/amTQtjEZE1nRJULzUyQc2eHZrzcn7wAxg/PtlziohkhRJULzWqo4Q7fOc7+aeUjx8fXkAoIrK2\nUILqpUbVoG66Ce64I0ybhTflqmlPRNYmSlC91IinSXR0wNFH5+ePPBJ23jmZc4mIZJUSVC81ogZ1\nyinw5ptheqON4KyzkjmPiEiWKUH1UtL3Qj3ySHjeXs6FF0Jzc7znEBHpD5Sgemn06HzCWLYM3n47\nvmN3dcG3v51PenvtBZ//fHzHFxHpT5Sg+iCpZr7LLoPHHgvTQ4bARRel/9ZeEZG0KEH1QWGCev75\neI751lvh2lPOKafAVlvFc2wRkf5ICaoPdtwxP93WFs8xv/996OwM01ttFeZFRNZmSlB98MlP5qfv\nuaf+jhL33w+/+U1+/uKL9aRyERElqD7YaScYOTJMv/YazJrV92OtXBnuc8o54AD41Kfqi09EZE2g\nBNUHTU2w2275+Xvu6fuxzjkHnn02TI8YAeefX19sIiJrCiWoPipu5uuLf/0Lpk7Nz0+dCptuWl9c\nIiJrCiWoPtpjj/z0vffCqlW92z/3MNjcK9wnToQpU+KLT0Skv1OC6qP3vx/GjQvTnZ35+5dqdeON\ncNddYdosPD2iqSneGEVE+jMlqD4y63sz35tvhtpTzlFHwUc+El9sIiJrgkQTlJlNMbNuMxtVZv0Y\nM7vYzF42syVm9qiZfTHJmOLUlwTlDocemn9d/LhxcMYZ8ccmItLfJZagzOzjwE+BkncJmdk6wD3A\nEcCDwEVAM/B7MzsqqbjitPvu+ekHH4QlS6rvc+WV8Oc/5+enTQvP9xMRkZ4SSVBmdiBwB1DpdtNj\ngRbgu+7+VXc/KZp/Gvipma2XRGxx2nhj2H77ML1yJcyYUXn7F1+E447Lzx9zTM9amIiI5MWaoMxs\nrJn9EbgOeAt4qcLm3462uSy3wN2XAD8BhgNfiTO2pBQmmNwbcEtZuhQOOihfy9p2W73nSUSkkrhr\nUNsDnwauAiYCr5fayMy2BMYBD7iv9qCg+6Lx5JhjS0Rhd/Nf/rL0s/mWLYPPfhYefjjMNzXBNdfA\nsGENCVFEpF+KO0G9COzg7oe6+8IK2+We071aDcvd3wKWA9vEHFsi9twz3wOvqwu+8AV49dX8+hUr\nYP/9e3aiOOecng+cFRGR1cWaoNz9dXd/uoZNx0bjjjLrFxI6TGReUxPcfDNssEGYnzsX9tsPXnkF\nbrsNPve5/P1OEHrsHXNMOrGKiPQnVW8NNbPZwGZVNrvY3Y/uxXkHReMVZdavAPpNA9imm4Yktdtu\n4YkSM2f2fGdUzg9+AKee2vDwRET6pVqeXXAzsH6Vbf7Ry/Mui8aDy6wfAtTQabun9vZ2Wltby65v\ni+vlTSXsuitceGG46baUE0+E009P7PQiIomr9Pva3t4e+/mqJih3Pz72s8KCaFyuGW8U8GYC503U\nkUeGruTnnRfe57TjjuH61N57h3um9Pp2EZHapfX0txei8RbFK8xsI8L9U71+mXpLS0uitaRqzODc\nc0Mz3siRMGhQ9X1ERPqLSr+vra2tTJ8+PdbzpZKg3P01M3sV2LXE6tyblh5sYEixWnfdtCMQEen/\n0nxY7DXAe8zs3ZdMmNlI4FRgKfDbtAITEZH0pfmCh58BXwQuMLNWwj1RBxCa/aa4+7wUYxMRkZQ1\nogZV8mGx7r6I0MR3VTQ+itB54kB3v6QBcYmISIYlWoNy992qrJ8DHJZkDCIi0j/phYUiIpJJSlAi\nIpJJSlAiIpJJSlAiIpJJSlAiIpJJSlAiIpJJSlAiIpJJSlAiIpJJSlAiIpJJSlAiIpJJSlAiIpJJ\nSlAiIpJJSlAiIpJJSlAiIpJJSlAiIpJJSlAiIpJJSlAiIpJJSlAiIpJJSlAiIpJJSlAiIpJJSlAi\nIpJJSlAiIpJJiSYoM5tiZt1mNqrM+gei9aWGw5OMTUREss3cPZkDm30cuBMYCoxx94UltukA/gP8\nvsQhbnf3mb04X0dzc3NzR0dHX0NuqNbWVgDa2tpSjaNWijdZijd5/S3m/hbv6NGj6ezs7HT30XEd\nsymuAxUyswOBKwnJqdw2mwOjgCvd/cdJxCEiIv1XrE18ZjbWzP4IXAe8BbxUYfMJ0fiJOGMQEZE1\nQ9zXoLYHPg1cBUwEXq+w7QTAUYISEZES4k5QLwI7uPuhpa45FcklqElm9piZLTaz18zsF+U6VYiI\nyNoj1gTl7q+7+9M1bj4hOv+PgMeAy4G3gWOAB8xsRJyxiYhI/1I1QZnZ7ApdwXPDhb05qZkZ0AHM\nBD7o7oe7+/fcfUfgMkJT4el9+HtERGQNUbWbuZmdC6xf5Th/cfffltj3PuDjlOlmXuZ8wwkdLDrd\nfVwt+0T7rQAGNzc3l92mpaWl1sMlrr29HchWTJUo3mQp3uT1t5izGG8uplI6OzsBVrr7kLjOl9h9\nUNC3BBXtNxPYARjm7itr3Odx4IPAsgqblf90RUSkmkrZchjwjLtPjOtkidwHVY2ZNROSyVx3n1Vi\nk2FAN/BOrceM80MREZH0pfUsvg8B/wv8vHiFmW0EbAnM9CSrdyIikmlpJagZwJvA3mY2KbfQzAYB\nvyTU7H6ZUmwiIpIBqTTxufs7ZnYY8AfgHjO7AZgH7AF8APidu/8mjdhERCQbGlGDKtlM5+5/BiYB\nfwH2AQ4HVgJT3P2gBsQlIiIZlmgvPhERkb7SCwtFRCSTlKBERCSTlKBERCSTlKBERCSTMp+gzGyg\nmR1nZk+b2VIze8nMTjOzXneRN7MBZvaQmXUnEWt0jrriNbMhZvYDM3vezJaZ2Ytmdl709I0sxjve\nzG41s/lRvP+MbiFIlJltYmYdZnZ0L/YZY2YXm9nLZrbEzB41sy8mGWfR+fsS845mdouZzTWzFdH3\n4SwzWyfJWKNz9zreov0TL29F5+vL59vQ8lZ07r7E29DyZmYbmtmlZvZq9P17w8yuMbMtaty/vjLn\n7pkeCE837wbagDOB6dH8DX041vHRvl1ZjJdwX9p9QBdwL3B2NO4GHgSaMhbvBGAx4faAa4BzCe8E\n6wbOSvAzHg48FH1OR9e4zzqE17q8A1wbfbazoliPasD3uC8x7wYsJzxf8rfAz6JjdEfjwVmKt8Qx\nEi9vdX6+DS9vdcbb0PIGbAi8GsV4F/BT4JZofg6wVZX96y5ziX5pYviAPhb9Mb8vWj4t+pD27sWx\ntgaWRPslUmDqjRc4odSXDbgo2v+gjMV7W7Tdp4u+lM9FX8rNE/iMN4++9N29LNynRNsfWbBsOPAk\nsAhYL4nvRJ0xPwusAHYsWn5ZdJxjsxRv0TESL28xfL4NLW8xxNvQ8gZcGp3vmKLlX41iv6XK/nWX\nucS+NDF9QNdGf+C2Rcs3jpb/oRfHagNeiP4xk0pQdcULvAy8BAwsWr45cBVwQMbifYPwwN/i5T+M\n9v9CzPEeC3QS/gf5114W7teA/xDd+1ew/MCosPWplpBUzMC2UVzXl1g3Llp3d1biLXGcxMtbDN+J\nhpa3GOJtdHl7E3ijzLpZwLIq+9dd5rJ+DWoS4R/k2cKF7v4G4cs/uZaDmNmR0bEOo/LrOOrV53jN\nbFtCwbjV3buK9n/F3b/p7jdnJd7IPGBUifb6TaPxnFiizDuG8KMyidDkZbXsZGZbEn7UH/CohBS4\nLxrX9F3qgz7FDCwEvg/8usS6FdE4ibdO9zXedzWwvEHfvxNplDeo7/NtWHkzswHATwhvPC9lBTA4\nen5qqf1jKXOZTVBmNpjwwb9UZpPZwGgzG1vlOO8htJ1e6e7TYw2y53nqjXd7wmOhnjGzvc1sRnRR\n8XUzOyfui+Ixfb6XAgOB68xsKzMbYWbfBA4hNGHE/XkfDkx094d7ud9W0Xi1v9Xd3yJc59mmztjK\n6VPM7v66u5/j7neVWL1/NH6q7uhW19fPGGhceSvQ13gbWt4K1PP5Nqy8uXu3u1/k7pcWrzOzDxCe\nmfqiu5d7JVIsZS6zCQpYNxp3lFnfGY2r9ba5jNDeeUIcQVVQb7ybEP439RngdmABcAmhWv894E4z\nGxhPqEAt7LaUAAAE4ElEQVQMn6+7Xwx8B/gkocq/ELiScKF5zxL/c6qLu/+1j8fMJdlyf+tCqn+P\n+qSOmEsysw2BHxN+XK+I67g5McTbqPIG1BVvo8sbUN/n2+jyVoqZGXAx4bO7vMKmsZS5LCeoXNVx\nRZn1ueVDyx3AzL4G/D/CA2gXxRhbKfXGOzwa7wMc6u6fdvcTgJ2AG4FdgaPiCDQSx+f7UeDkaNtp\nwAXAM4QCNDWWKONRy99a9u/MCjMbBfwZWB+4wN0fTTmkHhpc3urV6PJWt4yUt8uBTwCPROcvJ5Yy\nl+UElWu7Hlxmfe6990tKrTSzDYDzCBf6b4k5tlLqipdw0RDgcXe/Krcw+l/R/yf8jyXOe3bq/XxH\nEn4sRxCaLL7p7scRXgn9B+A70bWILKjlby3375IJZrYeoe1+IqE3V0NqKLVKobzVq9HlrS5pl7fo\nfslfA98idG3/nLuvqrBLLGUuywmqk/AlKlcNbC7YrpRfEf6+KTHHVU698eaWzyxe4e6vEqrKWxWv\nq0O98X4GGEP4n/y77czRlzb3mX+9/jBjsSAal/tbR1H+70ydmW1FuGemBbiV0FurITe/9kKjy1u9\nGl3e6pVaeTOzYcCfCNe6ngd2c/c3q+wWS5lL5YWFtfDwUsNXgHJ3LG8BzHH3cm2c+xPa6d8IzaY9\nDx/d3T7b3bfMSLyzonG5/3E0AUvrCLGHGOJ9D+Hzfa7Esd82s7nAZrEEW78XovFqf6uZbURoani+\noRHVyMxagLuB9QjNOodlMDlBg8tbDBpa3mKQSnkzs9GEm3Q/QuiIsZe7z61h11jKXGYTVGQGcJCZ\nbe3uL+YWmtnGhB4gt1bY9/Qyy78NbEC4dyDu/zXXE+8/CPdGTDYzK7zgGfWaGUG4byIr8b5FaAZZ\nrSdO9KUeC/wz3nD7xt1fM7NXCdcViu0WjR9sYEg1MbOtCS/0HAuc6+7fTzmkSk4vszzJ8laPNMpb\nPRpe3sxsCKFZcSdC8/Jn3X1xLfvGVuaq3SiV5gDsTvTYHQpu9gKuJtyYtlcfjvk4yd2oW1e8Bdud\nUrCsiVC97gL2yUq8hAv1iwkFZ4uC5QMI9+50ASck+N04hF7cYEu4iNxNuICfWzaS0FV7MTA2qVj7\nEjPhx+iR6HM8N+nY4viMyxwjsfIWw3eioeWtzu9Dw8sb4ZpiN/AAMKQP+9dd5jJdg3L3v5nZ9YSL\nlX83s/uAXaLhRne/M7etmZ0edvFyN5YlLoZ4TwB2BqaaWSvhf0S7AzsQHkf056zE6+5zzGwKoatz\nu5ndRGi3/wThmWFtVO7lk5gyn+3PCH/nBdFn+xJwAKEJYoq7z2t0nIVKxLwfsCPhfpGlZvbDEru9\n6e6XNSjEHrJQ3nojC+WtN9Iub9HtDEcRmhWfB04q0XQL4TFRKxMrc0n+DyGmLD4QOJXQZryU0AZ7\nCjCoaLtuYFUNx3u8lu3SipdwIfQXwCvR/s8Q7suwjMY7GbgTmE/oufM0oSvsoCTiLTjvIZR5TEyF\nWNcnFPA3CffqPELMj4eJK+boO9BVZZiZlXgrHCPR8hbDd6Kh5S2GeBtS3oDP1vD96wJGVYm3rjJn\n0UFEREQyJcvdzEVEZC2mBCUiIpmkBCUiIpmkBCUiIpmkBCUiIpmkBCUiIpmkBCUiIpmkBCUiIpmk\nBCUiIpn0f4dx/D4VGgDYAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x1094157f0>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "%matplotlib inline\n",
    "import matplotlib.pyplot as plt\n",
    "import potentials as p\n",
    "from potentials import *\n",
    "r, v_lj = lennard_jones()\n",
    "plt.figure()\n",
    "plt.plot(r,v_lj)\n",
    "x1,x2,y1,y2 = plt.axis()\n",
    "plt.axis((x1,x2,-15,y2))\n",
    "plt.show()\n",
    "\n",
    "\n",
    "\n",
    "#plt.plot([1,2,3,4], [1,4,9,16], 'ro')\n",
    "#plt.axis([0, 6, 0, 20])\n",
    "#plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.5.3"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 1
}
