Skip to content

Instantly share code, notes, and snippets.

@flixr
Created September 25, 2016 08:43
Show Gist options
  • Save flixr/3f92702c4086fd3e1a0c7e14816cb2ad to your computer and use it in GitHub Desktop.
Save flixr/3f92702c4086fd3e1a0c7e14816cb2ad to your computer and use it in GitHub Desktop.
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.text.Text at 0x7f3810c728d0>"
]
},
"execution_count": 1,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYEAAAEVCAYAAAAGrllxAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xd8VHW+//FXAGkqKlhWLLCyAjYUUUBddVREXVxB1FWu\nvQCK9eq1IK5kLctPUKzY1oYNFEUsKCDIQGgBRQURKaGsBRXUiBBAYub3x+eMMxknycnknDlT3s/H\nYx45c2bmnO+cJN/P+XYQERERERERERERERERERERERERERERqdY2wDfAe3H7QsBCZ/tw4DFn+zBg\nTJL9tfEIMDiFz2WKesCbwBJgQMJrhcDDCfu6AR87jzXA93HP/+F8Zq3zfD7wKfA+sK8fiZfs1yDo\nBEjOOR3LeA4F2gNfJLx+ALCns/0hcFaS/bURcR7Zak+gO9CUP36PZN9rMtDR2R4MtACuiXt9MDAq\nYd9VwMtYoBWppF7QCZCcMwB4A3gVuC7htT2BO4CjgaeJlRAS9x9LrOQAlUsSzZxjfwGEsUBT4Ly2\nBzAWCy6fAgOrSOMOwHNx7xsO1Hde2wK84hy/U5LnRwOznc/NA05yPncRUAR8BExJcs5kn9sOmICV\nnuYD+1SR3qoUEPvuifvjfYBdJxERX+0PbAJ2xKp6NgLNqZyJXwi87Wy72Z/4/H7gWWe7BbAKuN15\n/gFwqrPd2HkeLWnEewa7OwbL/F8AbnSeVwDnxr03/nkL4Ftid9T7Y1UvrbEg8AOWsSeq6nOtnMcv\nST4DdlefWB1U0+uJ+xoA92FVTiJ/oOog8dIVwHigFLvLXgn0B2bFvSfZnWt1+xOdAFzrbP8AvO5s\nN8VKEDsBdzr7tgUOJtbuEHUqliFf6jxvgmX2UUUJ748+7wIsx+7kAT4HZmJBKgIsADYkSXNVnzsO\nK814qQA4G/ir87wh9rvo6/F5JEcoCIhXtgUuAMqwzB+s6uZKYpmfWxEqB4WGCa/FV2P+5vyMVucc\nAWx2tnfGSiaJ6gFnYo2xYCWX+Pr3xIw8+jxZoKqP/R9tTfK5qOo+57UIMJrKbQIiVVKbgHjlXKyn\nSkvgz85jH6x6ZNe495VjdeCJ4vevBfYGdsEy0F5x75uA3cEXYJl39LVfgDnADc7zHbA7+NOSnGsi\ncL1zjIZYG0Ziz5xkioF2xKp1DsDq+sNUX5Kp7nPVqal0VFV7gNtSlYiCgHjmcqyBNf6O+mfgIayB\nOLp/FtZI+TqVe/bE7/8ceAKrxpiNdTmNvq8Qu+v+AmtD+CzufP8DdMWqZYqxXjKjkqT1GqzkssB5\nfAYMdV6rrofOOqyN4WHncy9hbQHLqb6XUnWfS3bO+HP3wwJc9DEj4fVk6c3m3lIiIiKS7ephPTBm\nANOxorCIiOSJk7G+1WAjHF8LMC0iIlIFv9oENmENcwXOz199Oo+IiGSgBljPhyVYo9gRgaZGRESS\n8qsr2a1Y74tB2JQAHwAHElciaNOmTaSkpMSn04uI5KwS4C9eHcyv6qBtgfXO9k9Y/+/68W8oKSkh\nEonoEYkwePDgwNOQKQ9dC10LXYvqH0AbLzNrv0YMD8PmdynCAsBAko/cFBGRAPkVBEqxKYVFRMQj\nzz/v/TE1YjgDhEKhoJOQMXQtYnQtYnQtoLgY/u//vD9ukHOMRJz6LRERqcaaNdC5M4wYAT17FoCH\nebdKAiIiGWzLFujdG/r1g9OSTYdYRyoJiIhkqEgE+vaFn36CMWOgXj0oKPC2JKD1BEREMtSjj1pb\nwOzZFgD8oJKAiEgGCofhnHNg1izYJ271aa9LAmoTEBHJMKtXQ58+8OKLlQOAHxQEREQySFkZ9OoF\nN90E3br5fz5VB4mIZIhIxEoADRvCyJFQkCSHVsOwiEiOGjoUVqyAadOSBwA/KAiIiGSAd9+Fhx6y\n3kBNmqTvvKoOEhEJ2NKl8Ne/wrhxcOSR1b9XvYNERHLI+vXQsyfcfXfNAcAPKgmIiASkosJ6Au25\npw0Mc0MNwyIiOWLwYPj5Z3jtteDSoCAgIhKA116z9QHmzbMuoUFRdZCISJotWAAnnAATJ8Khh9bu\ns2oYFhHJYj/8YO0ADz1U+wDgB5UERETSpLwcTjoJDjsM7rkntWN4XRJQEBARSZPrroMlS+Cdd6B+\n/dSOkU29gy4ELnK2mwAHA7sB6308p4hIRho50kYFFxenHgD8kK6SwCPAJ8BTcftUEhCRvDB3LvTo\nYXMC7b9/3Y6VjQ3DhwEHUDkAiIjkhTVr4Iwz4Kmn6h4A/JCOIHArUJiG84iIZJQtWywA9OtnU0Nk\nIr8Hi+0ItAWmJXuxsLDw9+1QKEQoFPI5OSIi6RGJwJVXwu67w6BBqR8nHA4TDoc9S1civ9sETgNO\nAK5N8praBEQkZ40YAY8/bovEb7edd8fNpt5BYKWAEp/PISKSUaZNgzvusEXivQwAftA4ARERD61e\nDV27wgsv+LNGcDb2DhIRyQvpXiTeCyoJiIh4wM0i8V7ItjYBEZG8MHQolJTA9OnpWyTeCwoCIiJ1\n9N57wSwS7wVVB4mI1MHSpXD00TB2LBx1lP/nU8OwiEiGiC4Sf9dd6QkAflBJQEQkBaksEu8FNQyL\niGSAwYOhtDTYReK9oCAgIlJLmbJIvBdUHSQiUgt1WSTeC2oYFhEJSKYtEu8FlQRERFyILhLfqZMN\nDAuKFpoXEQnA//4vLF4M48cHu0awegeJiKTZyJHwzju2VnAmLRLvBZUERESq4eUi8V5Qw7CISJpk\n+iLxXlAQEBFJIrpIfN++mbtIvBdUHSQikiASgUsvjY0IrpdBt8tqGBYR8dmtt8KiRTBlSmYFAD/4\n+fUGArOAecCFPp5HRMQzw4fDuHHWFTTTF4n3gl9BIAQcARzpbO/j03lERDzz/PPw4IMwaRLsvHPQ\nqUkPv6qDugMLgXFAM+BGn84jIuKJt9+2BeLDYdhrr6BTkz5+BYFdgL2AU7FSwFtAe5/OJSJSJ0VF\ncMklVgXUPs9yKr+CwDpgMVAOLAU2Azs7+39XWFj4+3YoFCIUCvmUHBGR5D791LqCvvwydO4cdGr+\nKBwOEw6HfTu+X11EewDXYtVCLYFpQFsgvk+ouoiKSKBWrLD1ge+/H/7xj6BT4062dBEdDxwDzMUa\nnwdQOQCIiATq22+he3e47bbsCQB+0GAxEck7P/8Mxx4LvXvD7bcHnZra0VTSIiJ1sGkTnHwydOhg\ni8MUBJkLpkBBQEQkReXl1gi87bbw4ovZORo4W9oEREQySiQC/frZxHBjxmRnAPCDgoCI5IWbb4bP\nP7f5gBo2DDo1maO6IDC7hs9GsGkhREQy2rBhtjJYUZFVBUlMdUGgADibquueRnmfHBERbz37LIwY\nATNmQIsWQacm81TXuNAGKKnD6zVRw7CI+Oqtt6B/f5sPqF27oFPjjSB6B90NXEpssFcEGwVcVwoC\nIuKb6dPhzDPh3XfhsMOCTo13gugd1ANoBWzx6qQiIn765BMLAKNG5VYA8IObTlIfA038ToiIiBdK\nSqBHD3j0UTjhhKBTk/nclAQ+A74BvnOeR9AiMSKSgdassfmAbr/dSgJSMzdB4Bzgz8DPPqdFRCRl\npaU2HcTFF1tjsLjjJgisAsqwNQFERDJOWRn8/e8QCsGgQUGnJru4aWEuxkoCK7CqIK8Gial3kIjU\n2datNhvoDjvYGsG5Ph1EOruI7gT8BLRO8tqquNdTpSAgInVSUWHLQn7/Pbz5JmyzTdAp8l86g0AY\nWxmsqs9NBEJ1OLeCgIikLBKBG2+EWbPg/ffzZzqIdI4TaA0s8epEIiJeGjoUJkywQWH5EgD8UFMQ\nEBHJOE8/DY8/bvMBNW8edGqymxaVEZGs8sYbMGAATJsGbdsGnZr008piIpK3wmFbFP6996BTp6BT\nE4ygVhZrC/wFWICNHq5w8Zn5xAaYrcAmoRMRScn8+RYAXnklfwOAH9wEgauBXkBz4AVsyoiravhM\nY+fncaknTUTELFsGp55q7QDHKVfxlJthFedgXUVLgeFAVxefORhoinUjnQJ0STWBIpLfvvkGTjoJ\n/vUvGxQm3nJTEiigcvWPm+kjNgLDgKeBfYH3sColN9VIIiIA/PSTBYDLLoO+fYNOTW5yEwRGAdOx\nNQXeA8a5+MxSYLmzvQz4Adgd+Dr+TYWFhb9vh0IhQqGQi0OLSD6IzgfUrRsMHBh0aoITDocJh8O+\nHd9NC3N7rNroQOALrHG4Jv2BDsCV2CpkU4ADqFwSUO8gEUlq61bo1cvGAIwcmfvzAdVGEF1EZwJH\n1fK4DYBnsdIDwE3AnIT3KAiIyB9UVMCFF1pV0Btv5Md8QLURRBCYBCzCqngqsFlEn/Tg3AoCIlJJ\nJALXXw/z5sGkSdC0adApyjxBjBOYhWX8u3p1UhGRZIYMgcmTbT4gBYD0cBtNugFtgNlYQ+8mD86t\nkoCI/O4//7EgMGMGtGwZdGoyVxAlgSHAHsD+wFZgINDHqwSIiIwdC4MH23xACgDp5abN/a/ABcAv\nwDPYKmMiIp547TW4/HIYPx723Tfo1OQfN0GgPrFpIOoDv/mXHBHJJ889B1dfDRMnQseOQacmP7mp\nDrof+AjYBZiLTR0hIlInDz8Mw4bB1KnQvn3QqclfbhsXdsJmEV0JrPPo3GoYFslDkYg1AD/zjPUE\nat066BRllyAahk8FrsAmhAPrLnq8VwkQkfwRicAtt8A771g3UDUCB89NNPkYuA74Lm7fFx6cWyUB\nkTxSUQFXXmkDwSZMgJ13DjpF2SmIksAPwDSvTigi+ae8HC6+GFavhg8+gGbNgk6RRFUXBPo7P3/F\npomYj1UFeTVthIjkgS1b4JxzYPNmKwFoJHBmqS4I7I5l+HOdn7ulJUUikjM2boTTT7c7/zffhIYN\ng06RJHJTr1SALS/ZFvgMGO/RudUmIJLDSkuhRw9o29amhGjgdkVzqZbXbQJuBos9CfwPsAW4ELjX\nq5OLSG5au9bWAu7UCZ5+WgEgk7mJJnOBznHP5+BuneGaqCQgkoO++gpOPBHOOAPuvBMKPLtnFQim\nJPBf4E/O9g5U7ioqIvK7khI45hjrCXTXXQoA2cDNr+hToDXwObZo/GbgJ6yxuEMdzq2SgEgOWbTI\nFoUfNAiuuCLo1OSuIFYW84uCgEiO+OgjawS+914477ygU5PbgqgOOhE4BegBrADO9erkIpL9iorg\nlFPg8ccVALKRmyBwN7a+8DXYgvOX+5oiEckaEyZA797w0kvQq1fQqZFUuAkCZcD32Kpia7DF5t3a\nFfgSG2MgIjnk9dfhggtg3DjrDSTZyU0QWA9MAF4FrsQCghvbAE8AG1NLmohkqpEj4aqrbDGYo44K\nOjVSF24aFxpji8wvAg7EFprf4uJzDwDvYmsS98eqlOKpYVgkCz3yCNxzD0yaBPvtF3Rq8k8QDcOb\nsQAANm2EmwBwEbAWmOQ8V29hkRwwZAjcf7+tBaAAkBv8Gsx9MTaOoBtwCDAS6EnCQLPCwsLft0Oh\nEKFQyKfkiEhdRCIwcCC8/bb1BtJiMOkTDocJh8O+Hd/NHfp5wIt1OMdUVB0kkrUqKqz+f+5cLQaT\nCYKoDurn1clEJLuUl8NFF8HChTBligJALnITTYqBRsASrHtoBJtVtK5UEhDJYFu2QJ8+UFYGY8dq\nMZhMEcTykjdjGb+I5InExWAaNQo6ReIXN9VB87GpIy4EmgNf+5oiEQlUaalNBNeyJYwerQCQ69wE\ngWewOYPaYovOP+1rikQkMGvXwvHHQ8eO8MwzWgwmH7gJAi2wQLAVmO7yMyKSZb7+2tYCOOUUeOgh\nqKf/9Lzg5tccAdo723sC5f4lR0SCsGIFHH209QS6+24tBpNP3PyqDwL+A+wHfAFcgbUT1JV6B4lk\ngM8/h+7d4dZbYcCAoFMjNQmid1AbbArp37w6qYhkhuhiMMOGwfnnB50aCYKb6qBu2BKTdwN/9jc5\nIpIu0cVgHntMASCfuS1SNAJOAy7Bpoju5sG5VR0kEpCJE20VsJdf1loA2SaIaSMAOgMnYYvETPHq\n5CKSfmPH2p2/FoMRcBdNPgcWYI3DXgYAlQRE0igSgQcegKFDYfx4OPTQoFMkqQiiYfgYrJtoG2Bn\nYJ1XJxeR9Cgrg379YNEimD0bWrcOOkWSKdxUBx0PzAZuxSaTUxOSSBZZtSq2BOTMmQoAUpmbIHA9\ncCjQC1sg5lpfUyQinpkyBbp2tQXhX3hBM4HKH7mpDvoN2OBs/wJs8i85IuKFSMSWgRw61HoAHX98\n0CmSTOUmCKwE7gOKgKOBEl9TJCJ1UlYGffvC4sVQXAytWgWdIslkbqqDLsYCQTdsNtG+vqZIRFK2\ncqXV/9erZ/X/CgBSkyCniVIXUREPTZ5sA8BuuQWuvVaTwOWqILqIikgGi0Rg+HC4914YNQqOOy7o\nFEk2cRMEtsHWEojaESj1JzkiUhtlZXDZZbBkCcyZo+ofqb3q2gR2B9phDcJtncd+wPsuj10fW4xm\nhnOMA1JPpogkWrkSjjwS6teHGTMUACQ11ZUEugLXYIHgCWdfBTDB5bFPdd7/V+BYbBbSXqklU0Ti\nTZ4M555rawBcc43q/yV1bv50/ga8m+Lx62PjDC4EQlhPoyg1DIvUUiQC991nj1GjIBQKOkWSbkE0\nDNcDxgPRsYYRbCoJN34DngNOB86sbeJEJCa+/r+4GPbeO+gUSS5wEwTuBK4DvkvxHBcBN2PzDu1H\n3IjjwsLC398UCoUI6bZGJKmVK+H006FDB6v/b9Ik6BRJuoTDYcLhsG/Hd1OkmExqi8icjy1MPwRo\nBnyCBYEtzuuqDhJx4f33rf//oEFw9dWq/893XlcHVXeg/s7PnsBX2OLyEefxpItjN8Gqgv6EdTMd\nArwd97qCgEg1IhHr+z98uOr/JSadbQK7Yxl+sfN8t1oeexNwdiqJEsl3Gzda/f+yZar/F3+5iSbH\nJjzfCvwXKx3UhUoCIkmsWGH1/4ccAo8/rvp/qSyd1UFR07FSwYdARywINMaWmxxah3MrCIgkiNb/\n33YbXHWV6v/lj4JYaL4MOAjoAxyMlQIOBM7wKhEi+S4Ssbn/L7gAXn1VDcCSPm66iO5CrEfPr9g6\nw1sIdgZSkZyxcSNceimUlMDcubDXXkGnSPKJm5LAOGz+n/uxOYDeBK4APvMxXSJ5YcUKOOIIaNwY\npk9XAJD0c3s3fzDQHliEZf67AOuw3kOpUpuA5LVJk+D88+Gf/4Qrr1T1j7iTzobhvljj75CE/RHg\nVg/OrSAgeSkSgWHD4IEHYPRoOOaYoFMk2SSd4wT+6/xcQt3u+EXEsXEjXHKJVQMVF6v6R4JXXZvA\nROfni9iI378Aq0h9RlGRvFZSYvX/TZtCUZECgGQGNw3DTwB7AycCOwEjfU2RSA6aONEWgOnfH555\nxhqCRTKBmyDQBrgdmwZiHLCDrykSySGRCNxzD1x8MYwZowZgyTxuxgnUx8YGAGyPrRYmIjX45Rfr\n/79qler/JXO5KQncBswCDsMmk7vD1xSJ5IBJk+Cgg2DHHdX/XzJbbQqmXowNiKcuopJzfvoJbrgB\nPvgAnnwSuncPOkWSa4KYO6gvNkhsLrDCeYhIgjffhAMPtN4/CxcqAEh2cBNN5mNrBMcvL7nZg3Or\nJCA5Ye1am/Bt/nx46ikN/hJ/BVESWAusxjL+6EMk70UituLXQQfZoi+ffqoAINmnut5B0ekiGgKT\nqLy8pBfTRohkrW++gSuusAFgb70FnTsHnSKR1FQXBKLTRXyBFT1UdyN5LxKxwV4DB1oQePVVaNQo\n6FSJpC7IYStqE5CssmoV9O0LP/5ogeDgg4NOkeSjINoEUrEN8AK2NGUx8HefziPiu4oKePhhOOww\n6NbNBn4pAEiucDNiOBXnYg3K52PzDX0CvO3TuUR8s3SpjfqtqICZM6Fdu6BTJOItv0oCY7D5hqLn\nKPfpPCK+KC+3NX+PPBLOOstG/SoASC7yqySw0fm5PRYQBvl0HhHPLVxoE77tuKOt+bvPPkGnSMQ/\nfgUBgL2AscAIYHSyNxQWFv6+HQqFCIVCPiZHpHq//gr//jeMGAFDhlg1kGb8lKCFw2HC4bBvx/fr\nT3w3IAwMAKZW8R71DpKMMW+eZfqtWsFjj8GeewadIpHk0rnGcF08CJyFjTWIOoXKo40VBCRwmzbB\n4MEwciTcfz/06aO7f8ls2RIE3FAQkEAVFdndf8eO1gV0112DTpFIzdK50LxITtqwwUb8jh0LjzwC\np58edIpEguNXF1GRjPT++zbd84YN8NlnCgAiKglIXigttcVeJk+GJ56Ak08OOkUimUElAcl5b71l\nd/+NGtkYAAUAkRiVBCRnrV0L11xj3T9fegmOPTboFIlkHpUEJOdEIjB6tC32sscesGCBAoBIVVQS\nkJwSXexl+XJb87dLl6BTJJLZVBKQnBBd7OWQQ6BDB1vvVwFApGYqCUjWmzoV/vlPKCuDSZMsEIiI\nOyoJSNYqKoLjjrPVvvr3txk/FQBEakclAck6s2fbfD/LlsHtt8P550MD/SWLpEQlAcka8+bB3/4G\nZ58NZ54JS5bYvP8KACKpUxCQjPfxx9Czp03xcOqpVgLo1w8aNgw6ZSLZT0FAMtbChXDGGXb3f/zx\n1u1zwAAb+Ssi3lAQkIyzeLFV+XTrZmv8lpTAtddC48ZBp0wk9ygISMZYuhTOO89G93bsaJn/DTdA\n06ZBp0wkdykISOBWrLAG3qOOgvbtrdrnlltgu+2CTplI7lMQkMCsXm19/Dt3trV9ly2D226DZs2C\nTplI/lAQkLT76iub3+fQQ21Jx6VLobAQdtwx6JSJ5B8FAUmbNWtsaucOHexu/4sv4O67oXnzoFMm\nkr/SFQS6AFPTdC7JMN99B9dfDwccYAO7Fi+Ge+6BXXYJOmUiko4gcBPwH0C9u/PMunVw882w335Q\nXm5r+g4fDrvtFnTKRCQqHUFgOdAbKEjDuSQD/PgjDBoE7drB+vXw6afw0EPQsmXQKRORROkIAmOB\n8jScRwJWWmoNvG3bwvffw0cfwWOPwV57BZ0yEalKoFNvFRYW/r4dCoUIhUKBpUVSt3693ek/+KDN\n7VNcDG3aBJ0qkdwQDocJh8O+HT9dVTStgVHAEXH7IpFIJE2nFz9s2AAjRsB990H37jatc9u2QadK\nJLcVFBSAh3l3OksCyvFzRFmZVfMMGwahEITDsP/+QadKRFKRriCwCjgyTecSH0QiMGsWvPACvPYa\nHHMMvP8+HHRQ0CkTkbrQchxSreXLLeN/8UWbv//8863Bt1WroFMmIl5QEJA/+OEHeOUVy/xXrIBz\nzrHnnTpBgTr6iuSUIP+l1TCcQbZsgXfesYx/6lQ45RS76+/eHbbZJujUiUiU1w3DCgJ5LBKBmTNj\n9fwdOljGf8YZsMMOQadORJLJ5t5BkiGWLYvV8zdubBn/xx/D3nsHnTIRSTcFgTyxbl2snn/lSujT\nB8aMsemcVc8vkr9UHZTDNm+O1fOHw7Zge7Sev4HCv0hWUpuAVKuiIlbP//rrcPDBsXp+rdglkv3U\nJiBJLV0aq+dv2tQy/k8+0eRtIlI9BYEstm4djB5tmf/q1VbP//rr0LGj6vlFxB1VB2WZzZvh7bct\n4582DXr0sLv+E09UPb9IPlCbQB6qqIAZM2L1/B07Wsbfu7fq+UXyjdoE8sQPP9i8/EVFVuWz7baW\n8S9YAHvuGXTqRCRXKAhkgK1bLXOfM8cexcXw7bdw+OFwxBEwdiwccojq+UXEe6oOCsBXX8Uy+zlz\nbLRu69bQtWvssd9+UL9+0CkVkUyjNoEsU1ZmUy9HM/w5c2yytvgM//DDVbcvIu4oCGSwSMTm349m\n9nPmwOLFcOCBltl36WI/99lHVTsikhoFgQxSWgpz51auy99uu8oZfseO0KRJ0CkVkVyhIBCQ8nJY\ntKhyXf6XX9pCK9EMv0sXaNky6JSKSC5TEEiTb7+tnOF/+KF1zYxm+F27WjWPBmiJSDplUxCoBzwK\ndAC2AJcBJXGvBxYEIhHYsMH64v/4Y+zn11/HqnfWr6+c4XfuDDvt5E96wuEwoVDIn4NnGV2LGF2L\nGF2LmGwaLNYLaAgcCXQB7nP2eWrTplhGnpipV7Xvxx9t0fQWLezRvLn93G03m275jjtg333T13ir\nP/AYXYsYXYsYXQv/+BkEjgImONvFwGHVvXnr1uSZdU2Z+2+/xTLz+Ay9eXPYZRdo377yvujPRo18\n/OYiIlnCzyDQDFgf9/w3rIqoIrqjU6dYpl5WZplzYoYd3W7V6o8ZeYsWNm2yuluKiKTGz+zzPmAO\nMMZ5/iUQP7v9cqCNj+cXEclFJcBfgk6EG72BZ53trsD4ANMiIiJpVgA8Bsx0Hm2DTY6IiIiIiASq\nHvA4MAuYSu62C2wDvABMx3pH/R2rx5vh7HuUWJtMX2AeMBvo4exrArzuvHc8sHO6Eu6jXbG2obbk\n97UYiP39zwMuJH+vRT3gGWLfvR35eS26YHkhePP9u2LtsTOA231Oe0p6Y794sC8/LsC0+OkiYLiz\nvRPwX+BN4Bhn32PYuIk/AQuwoNHM2W4IXE/sF3g28EA6Eu2jbYA3gC+wf/a3yM9rEcK+O8C2wL/I\n37+Lk4FXnO1uWIaWb9fiJuz7zHKee/F/8QnwZ2d7PHBIdQmoV9dvkIJajR/IYmOI/YLqAVuBQ7Go\nDfAe9od/ONZmshXrUrscG2Udf50mOO/NZsOwP+o1zvN8vRbdgYXYzc/b2D99J/LzWmwCdsDudncA\nfiX/rsVy7MY4esdf1/+L7bEAsdLZP5EarksQQaCq8QO5ZiOwAfuljAFuo/L3/AX7w28G/FzF/vUJ\n+7LVRcBaYJLzvIDK3ZPz6VrsgmV0ZwKXAy+Tv9diJtAYKx0+ATxE/l2LsUB53PO6fv/E/LXG6xJE\n5rseyxjj01BRxXuz3V7AB8DzwCgqf89mQCl/vB7bJ9kf3ZetLgZOxOo9DwFGYplhVD5di3VYMCwH\nlgKbqfxPmk/X4iYsELTD/i6ex6o8ovLpWkTVNY9IfG/0GFUKIgjMBP7mbHfF6rdy0W7YP/tNwHPO\nvo+BY53y42ZPAAAEGUlEQVTtU7Bi31zgaKARlhnsB3xG5esUfW+2OharCz8Oq6+8ACu+5uO1mIHV\nhQO0BJoCU8jPa7EtsbvWn7AZDPL1fySqrt//F6xabR+sVNGdDLwu+TJ+4EHgG+zuN/roAISxRqCn\niBX9LsN+0R8Cpzv7mgCvAkXAZKxnTS6Yiv3O9yV/r8U9xL7jieTvtdgR6yxQhPVmOYf8vBatiTUM\ne/H9u2C9iOYCd/qachEREREREREREREREREREREREREREcknC52f91N5Rbl4jYBL05McAK4CPgfO\nSuM5471ezWu7ASNcHqcQWIL1q6/JU9hArFwdjyMiGWphzW+hNTaoJV2mAAek8Xx+GQz0q8X7owP0\nRP4gFyduk2A0xSbDmo4tK1rf2R/G5oY5ChsVOh2bMmI7YBCwPza53h7YjJqTsADS0/n8AmxisTCW\nmTXDRlE+gs1C+zFwmvPeIdjoyVnYBG3x+mEzND6NBZ+FzvFuBDo6nws7advLeU8xNvnfPOAG53t9\nAtyd5Pufi43QLMKmSm+ATZw33XkcT2wG1c7Oe6dgc0o9C7QiFhCTfedE0ZGkYWwK4cnYzKQ3YNdw\nLjYiV0QkLa4nNkS9HTbdLVgm1g4YClyHZV49sYw2PuM7gdicKUcQm3F0JTbHFMCL2Lzpp2OZJ1hG\ndwc2H090X2MsOCTOnhi9I24NfI9l1GBD8Ts426dhGX8r5z3bY1U1Zc65GgHfJhy3BbAMmwsHbB2J\nK7EFY+LXy4gGgfnY/C8Ad2FBI/5aJPvO8eJLAlOBPs72e9jMpGDzVfWMe49KApKUSgLilXZYZgpW\nX7027rUI8G/sbn8Kdpe+lcrT5n4L9MdmkrycWAYNlqGDrUrWGMvQohlmKbZuQwdsiuapWGbYAMtY\nq7KS2BS+uxObyLCIWJXRCmxCrp+B75xzbXG+T7x9gEXY9OFgd/7RYyxJcu7dgcVx5ytI8p7E71yd\n+c7PUqzNA6wdoFENnxNREBDPfI5V+YAtGRq/1F8BcB52d3q8895+VF5L4g4sAFyAVXHE/20mZrqL\nsYU2wO7233X2TcVmKj0Ru5tfUU1646fs/QY4yNk+lljGnXjeqqzEqrWaOs9DccdINk36l8RKAkdU\ncUy3567te0UqaVDzW0RceZzYerGrgB/jXotgddRPYXfLv2FB4HtsFaT/h2Xa9wLXYm0Hzas4TwRr\nO+iG3UU3wHrLTMQy3+lYe8NYbFGfqsRnnH2xNoYCrIRyqbMdqeL9iZnuOqyKZiqW6S8DbsZ67yT7\n3ADsWm3Apv39qorjVnU+ERHJYgOIlZTuxBrGa2MwVnXmltoEpEoqCYik33dYw/cGrB7/whSOcT3W\nVjG6hvc9BRycwvFFREREREREREREREREREREREREREQkO/1/2o5l8W8JicwAAAAASUVORK5CYII=\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x7f382c685550>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"%matplotlib inline\n",
"import matplotlib.pyplot as plt\n",
"import numpy as np\n",
"\n",
"# spherical earth assumption with mean radius of 6371.009km\n",
"r = 6371009\n",
"x = np.linspace(0, 10000, 10)\n",
"plt.plot(x, np.sqrt(r**2 + x**2)-r)\n",
"plt.title(\"Altitude error of LTP\")\n",
"plt.xlabel(\"distance from origin [m]\")\n",
"plt.ylabel(\"height over sphere [m]\")"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": []
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 2",
"language": "python",
"name": "python2"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 2
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython2",
"version": "2.7.12"
}
},
"nbformat": 4,
"nbformat_minor": 0
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment