{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/usr/local/lib/python3.4/dist-packages/ipykernel/__main__.py:10: RuntimeWarning: invalid value encountered in true_divide\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAEWCAYAAACdaNcBAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAHoNJREFUeJzt3X+4XFV97/H3h8NvogSMngaChItQjT41mBTw0t6eKN4b\naAVawZJKlCongI0XRaVYNaUofUzVVm9BJaBFg01M7dVGiCLWnFJ/QElK5BoiEPlhAtgYCMghCCT5\n3j/2Hs8+w5wzc+bMr7335/U882T27D17rzUz+cw6a6+9RhGBmZkVy17dLoCZmbWew93MrIAc7mZm\nBeRwNzMrIIe7mVkBOdzNzArI4W49S9JHJW2X9HNJL5U0LKmv2+VqN0mfk/ThbpejFklDks7rdjms\nPod7D5H0gKSn0xCr3K7sdrm6QdJLgfcCsyLiNyLiZxExJSJ2d7tsEyEpJL1snPXnSvpe9rGIuCAi\nPtL+0lmR7d3tAtjzvDEivtPOA0jaOyJ2tfMYLfBS4NGI2NbtgpjlkVvuOVFp4Un6hKQdku6XdEpm\n/cGSPi/pEUkPpV0afZnnfl/S30l6FLhMUp+kT6bdHvdLWpy2MveWdJak9VXHv1jSv4xRtkMl/YOk\nh9OyfT2zblDSZkmPSVot6bDMupB0gaR7JT0u6SolTgZuBg5L/3q5TtLMSvnS5w5J+kharyclfVvS\ntMy+T5T0g3S/P5I0kFk3lL4+P0j3/w1JL5L0ZUm/lHS7pJmZ7V8u6ea0DndLenNm3XVpuW9My3Gb\npKPTdbekm/0oPc4fV71urwA+B7w2Xf94Zp8fTe8PSNoq6RJJ29L39wxJp0q6Jy3TX2T2uZekSyX9\nVNKjklZJOrTmh6oBkt4g6SeSnkj/ilTV+rdL2pS+7zdJOjKzLiS9M31/n0zfr6PT1/2Xadn2zWw/\n5mfFmhARvvXIDXgAOHmMdecCzwGDQB9wIfAwoHT914CrgYOAlwD/AZyfee4u4F0kf60dAFwA3AXM\nAA4BvgNEun4/4DHgFZnj3wG8aYyy3Qh8Jd3PPsDvpY+/DtgOvCbd598Dt2SeF8ANwFSSlvovgPnp\nugFga2bbmZXypctDwE+BY9P6DAEfS9cdDjwKnErSgHlDuvzizHM3A0cDB6evwz3AyWn9vwT8Q7rt\nQcAW4E/TdceldZqVrr8u3ffx6fovAyur6viycd7zc4HvVT12HfDRzOuwC1iSvraD6ev0j8ALgFcC\nTwNHpdtfBNyavq/7pZ+JFU1+HqcBTwJnpsd+T1qW89L1p6ev4yvSun8I+EFV3f8FeGFazmeAfwX+\nW+Z1f1sjnxXfmnj/ul0A3zJvRhLuw8Djmdtguu5cYHNm2wPT/zy/AfSn/3EOyKxfAKzNPPdnVcf6\nLmn4p8snMzo8Pwtckd5/JbAD2K9GmacDe4BDaqz7PPA3meUpJF9QM9PlAH4ns34VcGl6f4D64f6h\nzPp3At9K7/85sLyqLDdlgmQI+GBm3SeBb2aW3whsSO//MfDvVfu6GvjL9P51wLWZdacCP8kstyLc\nnwb60uUXpPs8IbP9euCM9P4m4PVV789zlddtgp/HtwK3ZpYFbGUk3L8JvCOzfi9gJ3Bkpu4nVZXz\nz6te90818lnxbeI3d8v0njMiYmrmdk1m3c8rdyJiZ3p3CnAkScvqkbQb4nGSAHpJ5rlbqo5zWNVj\n1eu/CPyJJAELgVUR8UyN8h4BPBYRO2qsOwx4MFPmYZJW7uG16kQSDFNq7GcsYz33SOCsymuRvh6/\nQxJ0Ff+Vuf90jeXsvk6o2tdbSL5UW1GHRjwaIyeSn07/Ha+8X8uUdROwm6QBMIqkb2rkxP1bahx3\n1GckktTNfk6OBD6dOdZjJF8A2fe30de5kc+KTYBPqBbDFpKW+7QY+0Rp9fSfj5D86V5xxKiNI26V\n9Czwu8CfpLexjn2opKkR8XjVuodJAgAASQcBLwIeGqcurbCFpOU+2KJ9/VtEvKEF+6ql1dOybgHe\nHhHfr3vgiFPqbPIImc9F+kWf/ZxsIfnr7svNFLRKtz4rheWWewFExCPAt4FPSnphelLtaEm/N87T\nVgEXSTpc0lSSroxqXwKuBJ6LiO/VWF859jeBz0g6RNI+kv5HunoF8KeSZkvaD/hr4LaIeKCpijbu\neuCNkv6XkhPH+6cnJmfUfebz3QAcK2lhWrd9JP12ejK0Ef9F0sc83voZ2ROLk/Q54IrKiU1JL5Z0\nepP7uhF4paQ/UnIi+38z+i+WzwEfkPTK9FgHSzqryWN167NSWA733vMNjR7n/rUGn/dWYF+Sk1Q7\ngK8yuhui2jUkXwh3kpwsXUNysiw7jnw58CqSsBzPQpL+0Z8A24B3A0QypPPDwD+TtAKPBs5usD5N\ni4gtJCf7/oLk5OMW4P008XmPiCeB/0lS7odJumCWkpz0a8RlwBfTros311j/XWAj8HNJ2ydavho+\nDawGvi3pSZKTqyc0s6OI2A6cBXyMpIvkGOD7mfVfI3ktVkr6JfBjoN5fA2MdqyuflSKrjLSwklMy\nrPJzEZH90/gAkrB+TUTc27XCmdmEueVeUpIOSMdK7y3pcOAvSYZTZl0I3O5gN8sft9xLStKBwL8B\nLycZtXAjcFFE/DJd/wDJyIczIuKObpXTzJrjcDczKyB3y5iZFVDXxrlPmzYtZs6c2fTzn3rqKQ46\n6KDWFSgHylbnstUXXOeymEyd169fvz0iXlxvu66F+8yZM1m3bl3Tzx8aGmJgYKB1BcqBstW5bPUF\n17ksJlNnSQ/W38rdMmZmheRwNzMrIIe7mVkBOdzNzArI4W5mVkB1w13SF9Kf9/rxGOsl6f+kP491\np6TXtL6YZmY2EY203K8D5o+z/hSS2eKOARaR/IKPmZl1Ud1wj4hbSH5hZSynA1+KxK3AVEnjTTU7\naddcA3fcAVJy6+sb/3699XnYx377wYYNSd3NzOppxUVMhzP6p7e2po89Ur2hpEUkrXv6+/sZGhpq\n6oC/+hUcdtgwn/hEc8/Pq+nTh9m2bYgmX7bcGR4ebvozkleuczl0os4dvUI1IpYBywDmzp0bzV6h\nde+98PDDQ7zvfcnz99oL9uxhzPv11udhH319sHTpEBs2DPCudzX1suWOr1wsB9e5PVoR7g8x+ncV\nZ9Dm3z0cHIShISjThJZHpK9wyRo4ZtakVgyFXA28NR01cyLwRPq7mtZCS5YkrfennnK/u5nV18hQ\nyBXAD4HflLRV0jskXSDpgnSTNcB9wGaS3+V8Z9tKW2KDg0m479gBl1/e7dKYWa+r2y0TEQvqrA/g\nz1pWIhvT9OkwdepI631wsNslMrNe5StUc2TaNJgyxa13M6vP4Z4zAwNJ90zJBheY2QQ53HNmaAh2\n7/aoGTMbn8M9Z5YsGd3vbmZWi8M9ZwYH3e9uZvU53HPI/e5mVo/DPYfc725m9Tjcc8j97mZWj8M9\nh9zvbmb1ONxzyq13MxuPwz2n3Ho3s/E43HPMo2bMbCwO9xyrjJpZscJdM2Y2msM9xypzvO/e7a4Z\nMxvN4Z5jg4OwYIG7Zszs+RzuOecLmsysFod7znlIpJnV4nDPOQ+JNLNaHO4F4Na7mVVzuBeAW+9m\nVs3hXhC+oMnMshzuBeELmswsy+FeEL6gycyyHO4F4QuazCzL4V4gvqDJzCoc7gXiIZFmVuFwL5Ds\nkMgLL3TAm5WZw71gfGLVzMDhXjg+sWpm4HAvJJ9YNTOHewH5xKqZNRTukuZLulvSZkmX1lj/Uklr\nJd0h6U5Jp7a+qNaowcHk3x074JJLulsWM+uOuuEuqQ+4CjgFmAUskDSrarMPAasi4jjgbOAzrS6o\nTYw0+l8zK5dGWu7HA5sj4r6IeBZYCZxetU0AL0zvHww83LoiWjOWLk26ZsBdM2Zl1Ei4Hw5sySxv\nTR/Lugw4R9JWYA3wrpaUzprmaYDNyk0RMf4G0pnA/Ig4L11eCJwQEYsz21yc7uuTkl4LfB54VUTs\nqdrXImARQH9//5yVK1c2XfDh4WGmTJnS9PPzaKJ13r4dtm5N7s+YAdOmtalgbeL3uBxc54mZN2/e\n+oiYW3fDiBj3BrwWuCmz/AHgA1XbbASOyCzfB7xkvP3OmTMnJmPt2rWTen4eNVPnqVMjIPk3b/we\nl4PrPDHAuqiT2xHRULfM7cAxko6StC/JCdPVVdv8DHg9gKRXAPsDv2hg39ZmPrFqVk51wz0idgGL\ngZuATSSjYjZKulzSaelm7wUGJf0IWAGcm37DWJf5xKpZOTU0zj0i1kTEsRFxdERckT62JCJWp/fv\nioiTIuLVETE7Ir7dzkJb4zzm3aycfIVqCbhrxqx8HO4l4K4Zs/JxuJeA53k3Kx+He0l4nnezcnG4\nl4TneTcrF4d7iVTmeV+xwl0zZkXncC8Rd82YlYfDvUQGB+Gzn/UPeZiVgcO9ZHxRk1k5ONxLyBc1\nmRWfw72EfFGTWfE53EvIXTNmxedwLyl3zZgVm8O9pNw1Y1ZsDveS8nwzZsXmcC8xX9RkVlwO9xLz\nfDNmxeVwLznPN2NWTA73kluyJBkxs3u3h0WaFYnDveQGB0dGzXhYpFlxONzNwyLNCsjhbh4WaVZA\nDncD3PduVjQOdwPc925WNA53+zX3vZsVh8Pdfs2zRZoVh8PdRql0yTzxhFvvZnnmcLdRli5NpiOI\n8HwzZnnmcLdR/CPaZsXgcLfncd+7Wf453K0m972b5ZvD3WpaujQJ+Ai33s3yqKFwlzRf0t2SNku6\ndIxt3izpLkkbJf1ja4tpnZa9qMmtd7P8qRvukvqAq4BTgFnAAkmzqrY5BvgAcFJEvBJ4dxvKah3m\nkTNm+dVIy/14YHNE3BcRzwIrgdOrthkEroqIHQARsa21xbRu8MgZs/xSRIy/gXQmMD8izkuXFwIn\nRMTizDZfB+4BTgL6gMsi4ls19rUIWATQ398/Z+XKlU0XfHh4mClTpjT9/DzqVp03bEgmFOvrg9mz\nO3dcv8fl4DpPzLx589ZHxNy6G0bEuDfgTODazPJC4MqqbW4AvgbsAxwFbAGmjrffOXPmxGSsXbt2\nUs/Po27V+ZBDIiBCili2rHPH9XtcDq7zxADrok5uR0RD3TIPAUdklmekj2VtBVZHxHMRcT9JK/6Y\nBvZtOeCRM2b500i43w4cI+koSfsCZwOrq7b5OjAAIGkacCxwXwvLaV3kkTNm+VM33CNiF7AYuAnY\nBKyKiI2SLpd0WrrZTcCjku4C1gLvj4hH21Vo6zyPnDHLl70b2Sgi1gBrqh5bkrkfwMXpzQqoMiXB\nJZeMjJypPGZmvcdXqFrDPOeMWX443G1CPOeMWT443G1Csn3vF17ogDfrVQ53m5DKVatScmGTu2fM\nepPD3SbMQyPNep/D3ZriC5vMepvD3Zri1rtZb3O4W9N8ctWsdzncrWmVk6t9fcnJVV+5atY7HO42\nKYODsGBBEvADA90ujZlVONxt0oaGkpb7ihXumjHrFQ53m7QlSzzu3azXONxt0jxyxqz3ONytJTxy\nxqy3ONytJTwtgVlvcbhby7h7xqx3ONytpTwtgVlvcLhbS7n1btYbHO7Wcj65atZ9Df2GqtlEVH6O\n7/zzk5Or558/+nEzaz+33K0tst0z7n836zyHu7VN5eQquP/drNMc7tY2g4Nw9dXufzfrBoe7tZUv\nbjLrDoe7tZ2HR5p1nsPdOsIXN5l1lsPdOsKtd7POcrhbx/jiJrPO8UVM1jG+uMmsc9xyt46qvrjJ\nLXiz9nC4W8ctXToS8B4eadYeDYW7pPmS7pa0WdKl42z3JkkhaW7rimhFMzgIO3bAIYckyz7BatZ6\ndcNdUh9wFXAKMAtYIGlWje1eAFwE3NbqQloxeXikWfs00nI/HtgcEfdFxLPASuD0Gtt9BFgK/KqF\n5bMCy/a/P/44HHqoW/BmraKIGH8D6UxgfkScly4vBE6IiMWZbV4DfDAi3iRpCHhfRKyrsa9FwCKA\n/v7+OStXrmy64MPDw0yZMqXp5+dREeu8fTs8+ODIcl8fzJ6d3C9ifetxncthMnWeN2/e+oio3/Ud\nEePegDOBazPLC4ErM8t7AUPAzHR5CJhbb79z5syJyVi7du2knp9HRa3zsmURUgQk/y5bljxe1PqO\nx3Uuh8nUGVgXdfI1IhrqlnkIOCKzPCN9rOIFwKuAIUkPACcCq31S1RpVmT3S/e9mrdNIuN8OHCPp\nKEn7AmcDqysrI+KJiJgWETMjYiZwK3Ba1OiWMRtLrf737du7WyazPKsb7hGxC1gM3ARsAlZFxEZJ\nl0s6rd0FtPKoTE8AyVDJrVu7Wx6zPGto+oGIWAOsqXpsyRjbDky+WFZG2ekJIpILnK65xtMTmDXD\nc8vkxbth9tBsmNrtgrTXIPD7x8A998DRn5nNT++Dhz8Bh03vdsk6Y/bjxX+Pq5WuzrOBM9p/GE8/\nYD3nsOmwd6bZce898PAj3SuPWR655Z4Xn4INQxsYGBjodkk6YtM1sO6pDbznPQMA6F64+n3F76Ip\n03tcUcY6M9T+Q7jlbj1pcDC5mKky/4yHSJpNjMPdelpl/hnwBGNmE+Fwt55WucCp8gtOixbBwoXd\nLpVZ73O4W88bHITPfnZk+frr3YI3q8fhbrkwOAjnnDOyfP75Dniz8TjcLTeWLx99gtU/0Wc2Noe7\n5Ur1T/S5BW9Wm8PdcqX6J/p8ktWsNoe75VK2BQ8+yWpWzeFuuVRpwfskq1ltDnfLNZ9kNavN4W65\n55OsZs/ncLfcq3WS1QFvZedwt8LIzkPjUTRWdg53K4zKPDQeRWPmcLeC8Sgas4TD3QqpehSNu2is\nbBzuVli1LnRywFtZONytsGp10Vx/PRx6qLtprPgc7lZ4y5ePDvgdO9wPb8XncLdSWL4cli0bPVTS\nAW9F5nC30qgMlfRYeCsDh7uVylhj4d0Pb0XjcLfSqXWi1f3wVjQOdyst98NbkTncrdTcD29F5XC3\n0hurH94Bb3nmcDdj7Aue9trLIW/51FC4S5ov6W5JmyVdWmP9xZLuknSnpH+VdGTri2rWftUXPEW4\nFW/5VDfcJfUBVwGnALOABZJmVW12BzA3In4L+CrwN60uqFmnVE60HnDAyGMeLml500jL/Xhgc0Tc\nFxHPAiuB07MbRMTaiNiZLt4KzGhtMc06a3AQdu58/nBJn2y1vFBEjL+BdCYwPyLOS5cXAidExOIx\ntr8S+HlEfLTGukXAIoD+/v45K1eubLrgw8PDTJkypenn51HZ6twr9d2+HR58cPRjfX0wYwZMm9ba\nY/VKnTvJdZ6YefPmrY+IuXU3jIhxb8CZwLWZ5YXAlWNsew5Jy32/evudM2dOTMbatWsn9fw8Klud\ne6m+y5ZFHHBARNILP3I755zWHqeX6twprvPEAOuiTr5GREPdMg8BR2SWZ6SPjSLpZOCDwGkR8UwD\n+zXLjUo3TfaiJ/CIGutdjYT77cAxko6StC9wNrA6u4Gk44CrSYJ9W+uLadYbKmPisydbKyNqfMLV\nekndcI+IXcBi4CZgE7AqIjZKulzSaelmHwemAP8kaYOk1WPsziz3sq34bMj7hKv1kobGuUfEmog4\nNiKOjogr0seWRMTq9P7JEdEfEbPT22nj79Es/9xVY73MV6iaTdJ4XTUHHeSuGusOh7tZC4zVVbNz\nZ9JV45C3TnO4m7VQ9uKnbFdNJeR90tU6xeFu1gbLl8OePaOvcIWRk65uyVu7OdzN2qgyT012OmFw\nd421n8PdrM0q0wlX98eDu2usfRzuZh1S6Y+PGLu7Zv16D6G01nC4m3XBWN01kAyh3Hdft+Ztchzu\nZl2S7a6ZOhX23ntk3XPP+eSrTY7D3azLKiH/3HNJaz07hBJG+uUld9lY4xzuZj3kqKOSIZS1Tr5C\n0mUjuTVv9TnczXrQeCdfYaQ17zlsbCwOd7Met3x5EvK1WvOVOWz6+hz0NprD3SwnsvPXVJ+A3bNn\nJOgl98+bw90sd7InYCtz2OxV43+yg77cHO5mOVaZw2b37rFPwoKDvowc7mYFkT0J22jQe9RNcTnc\nzQqo0aDPjqH3VbHF4nA3K7hGgz57Vay7cPLP4W5WItmgH2sMfVa2C8ct+3xxuJuVWGUMfaVVXz3E\nMivbsu/rc+u+1znczQwYPcSyXhfOnj0j97Ote19M1Tsc7mZWU3UXTjbsa42rh9oXUznwu8PhbmYN\nyYb97t21W/e1Qj8b+NnunEro339/5+pQJg53M2tadeu+1sVU2cDPdudUQv+xx0ZC3335reNwN7OW\nqhX4taZJqBX64/Xle9TOxDjczaztstMkVIc+jAT9eH35FWON2qncd/gnHO5m1jXLl8OcOSOhP5G+\nfHh+Nw/UD/+yfAE43M2s59Tq2qk3amesk7m17tf7AijCF4HD3cxypdaonbFa/fXCH2p/AUzki2Cs\n+92elM3hbmaFUi/8G/kCmMgXwVj3s5OyZcP/oINg+/bJ17OehsJd0nxJd0vaLOnSGuv3k/SVdP1t\nkma2uqBmZq1S7wtgot1AE+ke2rkTHnmktfWppW64S+oDrgJOAWYBCyTNqtrsHcCOiHgZ8HfA0lYX\n1MysGxr5S6DW/eykbNnAP/BAmD69/eVupOV+PLA5Iu6LiGeBlcDpVducDnwxvf9V4PWS1Lpimpnl\nS3ZStmz4P/UUTJvW/uM3Eu6HA1syy1vTx2puExG7gCeAF7WigGZmNnFjTO7ZHpIWAYsA+vv7GRoa\nanpfw8PDk3p+HpWtzmWrL7jOZdGJOjcS7g8BR2SWZ6SP1dpmq6S9gYOBR6t3FBHLgGUAc+fOjYGB\ngSaKnBgaGmIyz8+jstW5bPUF17ksOlHnRrplbgeOkXSUpH2Bs4HVVdusBt6W3j8T+G5EROuKaWZm\nE1G35R4RuyQtBm4C+oAvRMRGSZcD6yJiNfB5YLmkzcBjJF8AZmbWJQ31uUfEGmBN1WNLMvd/BZzV\n2qKZmVmzfIWqmVkBqVtd45J+ATw4iV1MAzpwEW9PKVudy1ZfcJ3LYjJ1PjIiXlxvo66F+2RJWhcR\nc7tdjk4qW53LVl9wncuiE3V2t4yZWQE53M3MCijP4b6s2wXogrLVuWz1Bde5LNpe59z2uZuZ2djy\n3HI3M7MxONzNzAoo1+Eu6SxJGyXtkVTYoVT1fgmraCR9QdI2ST/udlk6RdIRktZKuiv9TF/U7TK1\nk6T9Jf2HpB+l9f2rbpepUyT1SbpD0g3tPE6uwx34MfBHwC3dLki7NPhLWEVzHTC/24XosF3AeyNi\nFnAi8GcFf5+fAV4XEa8GZgPzJZ3Y5TJ1ykXApnYfJNfhHhGbIuLubpejzRr5JaxCiYhbSCagK42I\neCQi/jO9/yTJf/7qH8UpjEgMp4v7pLfCj+6QNAP4feDadh8r1+FeEo38EpYVSPoD88cBt3W3JO2V\ndk9sALYBN0dEoeub+hRwCbCn3oaT1fPhLuk7kn5c41bo1quVk6QpwD8D746IX3a7PO0UEbsjYjbJ\nDwAdL+lV3S5TO0n6A2BbRKzvxPE6+jN7zYiIk7tdhi5r5JewrAAk7UMS7F+OiP/b7fJ0SkQ8Lmkt\nyXmWIp9EPwk4TdKpwP7ACyVdHxHntONgPd9yt4Z+CctyTpJIfvRmU0T8bbfL026SXixpanr/AOAN\nwE+6W6r2iogPRMSMiJhJ8v/4u+0Kdsh5uEv6Q0lbgdcCN0q6qdtlarWI2AVUfglrE7AqIjZ2t1Tt\nJWkF8EPgNyVtlfSObpepA04CFgKvk7QhvZ3a7UK10XRgraQ7SRowN0dEW4cGlo2nHzAzK6Bct9zN\nzKw2h7uZWQE53M3MCsjhbmZWQA53M7MCcrhbbkiaKumdmeXDJH21Tcc6Q9KSNux3saS3t3q/ZtU8\nFNJyI51z5YaIaPtl6pJ+AJwWEdtbvN8Dge9HxHGt3K9ZNbfcLU8+BhydXuDzcUkzK3O+SzpX0tcl\n3SzpgbSFfHE6b/atkg5Ntzta0rckrZf075JeXn0QSccCz1SCXdIbJd2W7us7kvprPGcgOz+3pCsl\nnVu9XUTsBB6QdHyrXhSzWhzulieXAj+NiNkR8f4a619FMr//bwNXADvTFvIPgbem2ywD3hURc4D3\nAZ+psZ+TgP/MLH8PODHd10qSWf0mYx3wu5Pch9m4en7iMLMJWJvOhf6kpCeAb6SP/z/gt9IZF/87\n8E/JVC4A7FdjP9OBX2SWZwBfkTQd2Be4f5Ll3AY87y8Gs1ZyuFuRPJO5vyezvIfks74X8Hg6zex4\nngYOziz/PfC3EbFa0gBwWY3n7GL0X8L7j7P//dNjmLWNu2UsT54EXtDsk9P50e+XdBYkMzFKenWN\nTTcBL8ssH8zINMtvG2P3DwKzJO2Xznb4+nGKcizFntrWeoDD3XIjIh4Fvp/+WMvHm9zNW4B3SPoR\nsJHaP1l4C3CcRvpuLiPpylkP/Hr0jKS5kq5Ny7YFWEUS2quAOzLbXS7ptMz+TwJubrL8Zg3xUEiz\nGiR9GvhGRHynxfs9Drg4Iha2cr9m1dxyN6vtr4ED27DfacCH27Bfs1HccjczKyC33M3MCsjhbmZW\nQA53M7MCcribmRWQw93MrID+P/rBXB9URnc3AAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7fa8542c3a58>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "import matplotlib.pyplot as plt\n",
    "import numpy as np\n",
    "\n",
    "# Dle Martin Jerab \"Ziskavani vakua\", slide Prubeh primarni faze \n",
    "data = []\n",
    "for i in range(0,400,1):\n",
    "    data.append([i/100,np.exp(-(i-0)/100/1)])\n",
    "for i in range(-100,0,1):\n",
    "    data.append([i/100,1])    \n",
    "npar=np.array(data)\n",
    "t = np.arange(0, 400, 1)\n",
    "plt.plot(t/100,1/np.exp(t/t),color='magenta')\n",
    "\n",
    "plt.scatter(npar[:,0],npar[0:,1],s=3,color='blue')\n",
    "plt.xlabel('time (a.u.)')\n",
    "plt.ylabel('')\n",
    "plt.title('Energy confinement time - demo')\n",
    "plt.grid(True)\n",
    "plt.savefig(\"test.jpg\")\n",
    "plt.savefig(\"test.svg\")\n",
    "plt.show()\n",
    "\n"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.4.3"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}
