Skip to content

Instantly share code, notes, and snippets.

@mattions
Created July 26, 2013 15:35
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save mattions/6089841 to your computer and use it in GitHub Desktop.
Save mattions/6089841 to your computer and use it in GitHub Desktop.
Quick NGS file
{
"metadata": {
"name": "reading stuff in"
},
"nbformat": 3,
"nbformat_minor": 0,
"worksheets": [
{
"cells": [
{
"cell_type": "code",
"collapsed": false,
"input": [
"%load_ext rmagic "
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
"The rmagic extension is already loaded. To reload it, use:\n",
" %reload_ext rmagic\n"
]
}
],
"prompt_number": 12
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"%%R\n",
"X=c(1,4,5,7); sd(X); mean(X)\n",
"Y = c(2,4,3,9)\n",
"print(summary(lm(Y~X)))"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "display_data",
"text": [
"\n",
"Call:\n",
"lm(formula = Y ~ X)\n",
"\n",
"Residuals:\n",
" 1 2 3 4 \n",
" 0.88 -0.24 -2.28 1.64 \n",
"\n",
"Coefficients:\n",
" Estimate Std. Error t value Pr(>|t|)\n",
"(Intercept) 0.0800 2.3000 0.035 0.975\n",
"X 1.0400 0.4822 2.157 0.164\n",
"\n",
"Residual standard error: 2.088 on 2 degrees of freedom\n",
"Multiple R-squared: 0.6993,\tAdjusted R-squared: 0.549 \n",
"F-statistic: 4.651 on 1 and 2 DF, p-value: 0.1638 \n",
"\n"
]
}
],
"prompt_number": 13
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"cd /home/participant/NGS/velvet/part1/SRS004748/run_25"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
"/home/participant/NGS/velvet/part1/SRS004748/run_25\n"
]
}
],
"prompt_number": 14
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"%%R\n",
"library(plotrix)\n",
"data = read.table(\"stats.txt\", header=TRUE)\n",
"weighted.hist(data$short1_cov, data$lgth, breaks=0:50)"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "display_data",
"png": "iVBORw0KGgoAAAANSUhEUgAAAeAAAAHgCAMAAABKCk6nAAAC9FBMVEUAAAABAQECAgIDAwMEBAQF\nBQUGBgYHBwcICAgJCQkKCgoLCwsMDAwNDQ0ODg4PDw8QEBARERESEhITExMUFBQVFRUWFhYXFxcY\nGBgZGRkaGhobGxscHBwdHR0eHh4fHx8gICAhISEiIiIjIyMkJCQlJSUmJiYnJycoKCgpKSkqKior\nKyssLCwtLS0uLi4vLy8wMDAxMTEyMjIzMzM0NDQ1NTU2NjY3Nzc4ODg5OTk6Ojo7Ozs8PDw9PT0+\nPj4/Pz9AQEBBQUFCQkJDQ0NERERFRUVGRkZHR0dISEhJSUlKSkpLS0tMTExNTU1OTk5PT09QUFBR\nUVFSUlJTU1NUVFRVVVVWVlZXV1dYWFhZWVlaWlpbW1tcXFxdXV1eXl5fX19gYGBhYWFiYmJjY2Nk\nZGRlZWVmZmZnZ2doaGhpaWlqampra2tsbGxtbW1vb29wcHBxcXFycnJzc3N0dHR1dXV2dnZ3d3d4\neHh5eXl6enp8fHx9fX1+fn5/f3+AgICBgYGCgoKDg4OEhISFhYWGhoaHh4eIiIiJiYmKioqLi4uM\njIyNjY2Ojo6Pj4+QkJCRkZGSkpKTk5OUlJSVlZWWlpaXl5eYmJiZmZmampqbm5ucnJydnZ2enp6f\nn5+goKChoaGioqKjo6OkpKSlpaWmpqaoqKipqamqqqqrq6usrKytra2urq6wsLCxsbGysrKzs7O0\ntLS1tbW2tra3t7e4uLi5ubm6urq7u7u8vLy9vb2+vr6/v7/AwMDBwcHCwsLDw8PExMTFxcXGxsbH\nx8fIyMjJycnKysrLy8vMzMzNzc3Ozs7Pz8/Q0NDR0dHS0tLT09PU1NTV1dXW1tbX19fY2NjZ2dna\n2trb29vc3Nzd3d3e3t7f39/g4ODh4eHi4uLj4+Pk5OTl5eXm5ubn5+fo6Ojp6enq6urr6+vs7Ozt\n7e3u7u7v7+/w8PDx8fHy8vLz8/P09PT19fX29vb39/f4+Pj5+fn6+vr7+/v8/Pz9/f3+/v7///8o\n6bIMAAAUPklEQVR4nO2dfXxU1ZnHh4SEJogJTW1WqGaLwFABeduurwvR4FrsxhfYkF3LLhi17hZ8\nK76sXaWlrBF1WS1Lm5oVbG3Z+IKFbgWCJGYJkrAqKxpchKwSAjJTIikvCRNy/tnJuSeROffM3DOZ\nTMjz7O/7R57PPM+59z4z35mbc3Pv3PgEYI3vXDcAkgsEMweCmQPBzIFg5kAwcyCYORDMHAhmDgQz\nB4KZA8HMgWDmQDBzIJg5EMwcCGYOBDMHgpkDwcyBYOZAMHMgmDkQzBwIZg4EMweCmQPBzIFg5kAw\ncyCYORDMHAhmDgQzB4KZA8HMgWDmQDBzIJg5EMwcCGYOBDMHgpkDwcyBYOZAMHMgmDkQzBwIZg4E\nMweCmQPBzIFg5kAwcyCYORDMHAhmDgQzB4KZA8HMgWDmQDBzIJg5EMwcCGYOBDMHgpkDwcyBYOZA\nMHMgmDkQzBwIZg4EMweCmQPBzIFg5kAwcyCYORDMHAhmDgQzB4KZA8HMgWDmQDBzLAUfaUluGyBZ\neAu+4bBouiolraC5H7oBfY63YF+jKJ53vP2B2f3QDehzrASP2iNEIKsfuonOwmmSP8V+JE4sBG8L\nFb4hRNXYfugmOvlOuPe9c9oFQbwFT79oSPYkUZ1d1g/dRAeCe4nNLLp9b62oq0l6KzGB4F5idZjU\neexMsvvwAoJ7ibfgE8vGpPtSRy9p64duogPBvcRb8PyZVcHTwZrCkn7oJjoQ3Eu8BWcdlKE1J9mt\nxASCe4m34MnlMlRMSXYrMYHgXuItuD5vXFHJ3PEjdvZDN9GB4F5iMYsOVZaXllWGkt9LLCC4l1ie\nTVqT3C68geBeYik4N7ldeAPBvaQ3gj+pkLz4UTIaigIE9xIbwW/5/Sl+v7/n8btlku/+MHltuYDg\nXpLALvrf/7VPO4kNBPcSCGZOArNoCKaAneCO5k53EoIpYCO4KT8zu/6a/XoagilgI/j6RafyQosL\n9DQEU8BGcEaLyBOBTD0NwRSwEXzZ+rDgDRP0NARTwEZwVc6cYbMv2KinIZgCVrPo4Oql5YdcWQim\ngOVx8OeG7yZBMAW8BX947ZzgXwxJvfagXoBgClhc+L5wcc5DbScX3awXIJgC3oIzDrcOOilE4Hy9\nAMEU8BZ8ceV233YhNn1DL0AwBbwFr0o7f9WFdyzIek0vQDAFLGbRe5tFw/InPnDlIZgCCdyjA4Ip\nAMHMgWDmQDBzIJg5EMwcCGYOBDMHgpkDwcyBYOZAMHMgmDkQzBwIZg4EMweCmQPBzIFg5kAwcyCY\nORDMHAhmDgQzB4KZA8HMgWDmQDBzIJg5EMwcCGYOBDMHgpkDwcyBYOZAMHMgmDkQzBwIZg4EM8dS\n8BHcTpgo3oJvOCyarkpJK2jWCxBMAW/BvkZRPO94+wOz9QIEU8BK8Kg9QgSy9AIEU8BC8LZQ4RtC\nVI3VCxBMAYs7vl80JHuSqM4u0wsQTAGbWXT73lpRV+NKQzAFrA6TOo+dMWQhmALegk8sG5PuSx29\npE0vQDAFvAXPn1kVPB2sKSzRCxBMAW/BWc4/TGrN0QsQTAFvwZPLZaiYohcgmALeguvzxhWVzB0/\nYqdegGAKWMyiQ5XlpWWVIVcegilgeTZpjSEHwRSwFJxryEEwBSCYOTaC3/L7U/x+f8/j12dKJj6Y\nvLZcQHAv6c0nuP2o5IWfJKOhKEBwL8EumjmYRTPHTnBHc6c7CcEUsBHclJ+ZXX/Nfj0NwRSwEXz9\nolN5ocUFehqCKWAjOKNF5IlApp6GYArYCL5sfVjwhgl6GoIpYCO4KmfOsNkXbNTTEEwBq1l0cPXS\n8kOu7LkQXPBN+Ue06W/046ZpQ+y7SZOdT3DFyn7cNG28Be/J3/XpFalp+Qf0AgRTwFvw5fe3zV54\nqu2+b+sFCKaAt+DzDovR/yNEYJhegGAKeAv+1jOdd/1UiJem6gUIpoC34INTx92Ukj8jd4degGAK\nWMyiO+tWL1+53vXFBggmAbF7dEBwvEAwcyCYORDMHAhmDgQzB4KZowTfU90R96IQTAEl+LHLcv9+\ni/sLhDGBYAr07KL3PfNnuXdtPB3HohBMgR7BLRULho+/cuQ6+0UhmAJK8PIZQ//8uf1CbP0j+0Uh\nmAJK8O2vtsp4/DX7Rc+h4Cf8zhccV/RjB0RRgtvKq8UvV7THteg5FHy/8wlu/ut+7IAoSvDd094T\ntd903QorJhBMASX4y/vCP/a57hgcEwimgBJ8ydvhH7WXxLUoBFNACV6bc8/T9+X8Iq5FIZgC3cfB\ne5be/fiu+BaFYArQPNkAwdYowZVX+v1n30jHBgimgBJ88cPvN4SJa1EIpoASfOHJ+BeFYAoowU+V\nkjofDMHWKMFXDz1/LH4Hc0QJbnCIa1EIpkDPYZLxVlgxgWAKKMFRboUVk34R3HF7keRrzkMIjhcl\nOMqtsGLSL4JP5O+TnO88hOB4UYKj3AorJv0j+EYnZjsBguNFCY5yK6yYQDAFlOAot8KKCQRToHsW\nbb4VVkwgmAID/WwSBCeIEny5Q1yLQjAFlODt27fXVsx4Oa5FIZgCZ++iAxPjWhSCKXC24Hdc9zqL\nCQRT4KzfwdMG3xvXohBMgS9+B2/f3hDf6QYIpgAOk5ijBI/M6sY87JzdLxqCE0QJXnXdtuba655t\naXGLvOGwaLoqJa2gWS9AMAW6P8EHwz+avmYc0SiK5x1vf2C2XoBgCijBF1WHf1SPMI5oFKP2hA+S\nXTtvCKaAElw2/JHnHxn+lHHEtlDhG0JUjdULEEyB7ll03feLF20xHiZNv2hI9iRRnV2mFyCYAjYX\n3bXvrRV1Na40BFPA6qK7zmNnDFkIpoD3RXcnlo1J96WOXuK65TsEU8D7orv5M6uCp4M1ha4beEAw\nBbwvuss6KENrjl6AYAp4X3Q3uVyGiil6AYIp4H3RXX3euKKSueNH7NQLEEwBJXji7uhDQpXlpWWV\n7lvRQjAFlOClC9z/F+ls1hhyEEwBJXhGVvolsb4fnGvIQTAFlOD3HaIOixDcflTywk+S2pkDBCeI\nFDy0RYhf/SHqmLf8/pSzP94bnG90XvFg8tuD4ESRgn1hwVmNsYZhF00VCGaOpWDMoqniCH6zvn7o\n6/X19dFGGc8lQjAFpOCcbsxjopxLhGAK2FwXHeVcIgRTwEZwlHOJEEwBG8FRziVCMAVsBEc5lwjB\nFLD6bpL5XCIEU8D2y2fb3WebIJgCtoJzDrhSEEwBb8FDU7vwpaTqBQimgLfgD68s2hcIDN8V0AsQ\nTAGLXXTHCv9vsYumitXv4I/z5w2DYJrYTbLOlBUHXUkIpgDu0cEcCGYOBDMHgpkDwcyBYOZAMHMg\nmDkQzBwIZg4EMweCmQPBzIFg5kAwcyCYORDMHAhmDgQzB4KZA8HMgWDmQDBzIJg5EMwcCGYOBDMH\ngpkDwcyBYOZAMHMgmDkQzBwIZg4EMweCmQPBzIFg5kAwcyCYOZaCj7S4cxBMAW/BNxwWTVelpBU0\n6wUIpoC3YF+jKJ53vP2B2XoBgilgJXjUHiECWXoBgilgIXhbqPANIarG6gUIpoC34OkXDcmeJKqz\ny/QCBFPAZhbdvrdW1NW40hBMAavDpM5jZwxZCKaAt+ATy8ak+1JHL3H9Xx0IpoC34Pkzq4KngzWF\nJXoBgingLTjroAytrv89C8EU8BY8uVyGiil6AYIp4C24Pm9cUcnc8SN26gUIpoDFLDpUWV5aVhly\n5SGYApZnk9YYchBMAUvBuYYcBFMAgpljI/gtvz/F7/f3PH63VHLbj5LXVg8QnCC9+QQf2Cz5wZPJ\naEgDghMEu2jmYBbNHDvBHc2d7iQEU8BGcFN+Znb9Nfv1NARTwEbw9YtO5YUWF+hpCKaAjeCMFpEn\nApl6GoIpYCP4svVhwRsm6GkIpoCN4KqcOcNmX7BRT0MwBaxm0cHVS8sPubIQTAEbwUe7jpE6Anoa\ngingLXj3pYNGrRei0TUSgingLfjqf2yvHlkPwUTxFpx5TIh10zogmCbegie9LETnLY9CME28BW86\n74rPRGDKZAgmicUsunltqxBtax/S8xBMgQF7j44Tr1R08UsITowBK3jTLWVdPAPBiTFgBW9cJkMA\nghMDgpkDwcyBYOZAMHMgmDkQzBwIZg4EMweCmQPBzIFg5kAwcyCYORDMHAhmDgQzB4KZA8HMgWDm\nQDBzIJg5EMwcCGYOBDMHgpkDwcwZeIKPOPeyLYXgPmHgCX78toe7mG4h+MBU587Vv0pKIzwYeIIf\nq5bhnywEv3eV/IppxZ8kpREe0Bb8V042PymN8ACCmQPBzIFg5kAwcyCYORDMHAhmDgQzx0pw57Ez\nhiwEU8Bb8IllY9J9qaOXtOkFCKaAt+D5M6uCp4M1hSV6AYIp4C0466AMrTl6AYIp4C14crkMFVP0\nAgRTwFtwfd64opK540fs1AsQTAGLWXSosry0rDLkyve14IXTJF+F4L5kAB0mXedsZDoE9yUD6DAJ\ngpPBADpMguBkMIAOkyA4GQyEw6SJzuxqKAQngYFwmKT8DIfg2Pz+5/LfWJTVx7PQQDhMgmA7Xr1d\nXiT8/Nx4FrI8XbjGkEtY8AffmCnpY8G/dy6Wrvg4IvuOyrYk2PQ55NVnZfgsGYJzDbmEBb/9sBMT\nFjzB2XWt/GhfF6XF8usOi0YXSZY4Y25dIcf85W/ko443nS/I7InV3++ct8SWxJ5l30FHcPloOavK\n6yvB5zsmRtwojf7xL2TydwuOSqY7Y249KsNTv5bJLZfL78c8PNb57stdzqN7XpHWN9XJN0qDKo5y\n3gtvOu3tct42RSuctX/mhA+cMZuPx3jWJ3Y69HofkjTBb/n9KX6/v+fxy86sd8zSiFFfyZQMzpAh\nPV2GjMEqGxmcYnpkdpATUr/kFNPUw4hiilNMc4oZxqJa+xBVTE+VpKjVpTgP1TadB6mDIkKKGtO9\n9shi6pAu0lNiLTl4iDMoYtNpkZtOd8aoZIrzUD2xzDSnmKpCRsRL8vW+FizMn2AN9Tmc/78y/Pwl\nGRr+LqJ4i/PufdrZT2q7aHUcrE42qK+Pdv/vQjXmXucTXOF8gg9F/g5e0ChDufMJ3nN3RFF9gp95\nXYYdD0UUCzpkeLxKhk0/luHUtyLG3PeuDK88J8Ph4oji7ftl+LcXZfjouxHFOUEZ/nmdDHUPRhRn\nOrPXx7fKsNn50LTdEDHm/ndk6O9dtAYEUxZsmkVrQDBhwR3NnZ5jIJis4Kb8zOz6a/Z7jIJgsoKv\nX3QqL7S4wGMUBJMVnNEi8kQg02OU6qmkSYYX1sqwd1FEsahVhn/5Dxl2PhpRnOUI/uE2GbY8KcOp\nmyPGPLBbhtd+JkNgXkTxjgMyrPm1DB8vjNz0MRme/a0M//UPEcUbnU3/6D9leLNUhvabIsZ8/30Z\n1q1yNv2diOKdn8rwonMziX3fiygWO2/r5zbI8M4jEcVvO++tpTUybH3C2XRhxJjF/y3D684fHoK3\niTiwEXzZ+rDgDRM8Rp2KCKEOU1aF006xs81UbHde7DPtpmKbMxfoOB1j7aHQ/49N22EjuCpnzrDZ\nF2yMZ7VgoGA1iw6uXlp+KNmdgKSQwEV3gAIJXHQHKJDARXeAAglcdAcokMBFd4ACCVx0ByiQwEV3\nEdRPyV5wUk92PnZhxjUf6NnmWcMud18q85p/+N/o10PMagj/2Dwp8+rdhi3KonldahFjT1+sr2Go\ne6Wq+LO8jBkN0Z6NsaHuorkhtTFzQz2dmBpSWWNDVlieLvQilPd8U8EyPVt58YdHSmZpyc6pTx+6\n33UdZEP21sCt90Yufacv/Jyaz6v4/AeXujboFM3rUouYe+pZX8eVqa61quLetMpDC6+N8mzMDami\nuSG1MXNDPZ2YGlJZc0NW9JHgynFCbB2jZ/fXdX7+6HwtWTeuU7Tt0oc+tUCIPZHzuKcWZoYdrr1C\niPZBR13jZdG8LrWIuaee9a0ocr+eqtg87O3WxXOiPBtzQ6pobkhtzNxQTyemhlTW3JAVfSS4vEiI\nYLr7pPFa31catdTqm+8aM+dTfeCP7wy/SD5tHz0y7LD1MyGqv244Hd1VNK9LLWLuqXt9H/v3uV/P\n7uIq36CcoHuT8tlEa0gWzQ2pjUV5kVQnxoa6s1EasqCPBJeGj5JP+465CycW6zus5YNWffS9K/Vx\nO4bXBef6jkQmR8pfO53rRv7GsMmuonldapFoPcnimWs3BEwfGFlsuLD25IOFhqrzbKI01FU0N6Q2\nZm5IFc0NqWzUhrzpI8Flc8NvzsH6HzQ/+qTrLav9CWzlDCFOpgb0Nfx0ZO6y1I7InBQcvHWa8bsa\nXcUo63IWMffUXfyOMAl2isvvEKJtyOdRno2xIVU0N6Q2Zm6ou2hsSGXNDVnRR4I3jxeiZrSefTI8\na/r0S5q19eHXoG2wfnFwa5MQtRO1ZJfDtqmPais4q2hel1rE3JMqFg/LGe7L2W4uPhGeEJxyrVY9\nG3NDqmhuSG3M3JAqmhtSWXNDVvTVLHrEy3+4aYmerf/qjsDf6tcftOW+FLx/hj70w6E1n0xfpSW7\nHK6d1BjG8JJK+8Z1qUXMPali8MCBXSkH9L+vq+LurM3Be1yTVvVszA2porkhtTFzQ6pobkhlzQ1Z\n0UeCRf2kLy9wn4544dKsua4d6I6p58064Bq68oKRS/T5R5fDh3xduFai9t/GdXUvYuzpi/UZ9ojd\nxVf9wwqbojybKA2pp2p+cmpj5heppxPjpMDJmhuyoa8EgwEKBDMHgpkDwcyBYOZAMHMgmDkQzBwI\nZg4EMweCmQPBzIFg5kAwcyCYORDMHAhmDgQzB4KZA8HMgWDmQDBzIJg5EMwcCGYOBDMHgpkDwcyB\nYOZAMHMgmDkQzBwIZg4EMweCmQPBzIFg5kAwcyCYORDMHAhmDgQzB4KZA8HMgWDmQDBzIJg5EMyc\n/wMFeAHUJaP2CgAAAABJRU5ErkJggg==\n"
}
],
"prompt_number": 15
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"pwd()"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "pyout",
"prompt_number": 16,
"text": [
"u'/home/participant/NGS/velvet/part1/SRS004748/run_25'"
]
}
],
"prompt_number": 16
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"cd /home/participant/NGS/velvet/part2/SRX008042"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "stream",
"stream": "stdout",
"text": [
"/home/participant/NGS/velvet/part2/SRX008042\n"
]
}
],
"prompt_number": 18
},
{
"cell_type": "code",
"collapsed": false,
"input": [
"%%R\n",
"library(plotrix)\n",
"data = read.table(\"run_21/stats.txt\", header=TRUE)\n",
"weighted.hist(data$short1_cov, data$lgth, breaks=0:50)\n",
"title(\"Normal\")\n",
"data2 = read.table(\"run_21trim/stats.txt\", header=TRUE)\n",
"weighted.hist(data2$short1_cov, data2$lgth, breaks=0:50)\n",
"title(\"Trimmed\")"
],
"language": "python",
"metadata": {},
"outputs": [
{
"output_type": "display_data",
"png": "iVBORw0KGgoAAAANSUhEUgAAAeAAAAHgCAMAAABKCk6nAAAC61BMVEUAAAABAQECAgIDAwMEBAQF\nBQUGBgYHBwcICAgJCQkKCgoLCwsMDAwNDQ0ODg4PDw8QEBARERESEhITExMUFBQVFRUWFhYXFxcY\nGBgZGRkaGhobGxscHBwdHR0eHh4fHx8gICAhISEiIiIjIyMkJCQlJSUmJiYnJycoKCgpKSkqKior\nKyssLCwtLS0uLi4vLy8wMDAxMTEyMjIzMzM1NTU2NjY3Nzc4ODg5OTk6Ojo7Ozs8PDw9PT0+Pj4/\nPz9AQEBBQUFCQkJDQ0NERERFRUVGRkZHR0dISEhJSUlKSkpLS0tMTExNTU1OTk5PT09QUFBRUVFS\nUlJTU1NUVFRVVVVWVlZXV1dYWFhZWVlaWlpbW1tcXFxdXV1eXl5fX19gYGBhYWFiYmJjY2NkZGRl\nZWVmZmZnZ2doaGhpaWlqampra2tsbGxtbW1vb29wcHBxcXFycnJzc3N0dHR1dXV2dnZ3d3d4eHh5\neXl6enp8fHx9fX1+fn5/f3+AgICBgYGCgoKDg4OEhISFhYWGhoaHh4eIiIiJiYmKioqLi4uMjIyN\njY2Ojo6Pj4+QkJCRkZGSkpKUlJSVlZWWlpaXl5eYmJiZmZmampqbm5ucnJydnZ2enp6fn5+goKCh\noaGioqKjo6OkpKSlpaWmpqanp6eoqKipqamqqqqrq6usrKytra2urq6vr6+wsLCxsbGysrKzs7O0\ntLS1tbW2tra4uLi5ubm6urq7u7u8vLy9vb2+vr6/v7/AwMDBwcHCwsLDw8PExMTFxcXGxsbHx8fI\nyMjJycnKysrLy8vMzMzNzc3Ozs7Pz8/Q0NDR0dHS0tLT09PU1NTV1dXW1tbX19fY2NjZ2dnb29vc\n3Nzd3d3e3t7f39/g4ODh4eHi4uLj4+Pk5OTl5eXm5ubn5+fo6Ojp6enq6urr6+vs7Ozu7u7v7+/w\n8PDx8fHy8vLz8/P09PT19fX29vb39/f4+Pj5+fn6+vr7+/v8/Pz9/f3+/v7///8YLzCuAAATLklE\nQVR4nO2dfXwUxRnHlwQCEimJUSmmQgVtlFRepFW0BV9oKdj4RiG+SwGx9QURVIoFW0pBsfWtVJTW\n1tYWQ6VYX6qgQKmJEGq1IgEMVEhBuRBEXhIImT97NzMXsrtz3Ca7s3vz5Pf9fGS8eWZ3nrvv7d7u\n3W4eiwHSWFEnAPQCwcSBYOJAMHEgmDgQTBwIJg4EEweCiQPBxIFg4kAwcSCYOBBMHAgmDgQTB4KJ\nA8HEgWDiQDBxIJg4EEwcCCYOdcGDLcEE8fCQle11yQ+sIl1JhQl1wWOLi/Osk4qLfywe1hf397ok\nBJvCBGtqWxaDYFMQgndZBeWDlyR20VusPr8s7PVYxaDcYdX2B2zblQUFpTsY+/SavP7PQLAhJAXn\nnm5JwVZOUfy/vrnW1fYHn/fsWDLCOmN/Y3/rSwOzINgQkoKtubX1UnAVG2pd1/SW1df+4HFrfCx2\nifXCMuuc/U0TINgQkoK7HBVH0VusQsZusJ6J/09v+4OJ4oj7J7Os+DHZKgg2hObPYJYU3Dvh9Fkp\nuMWDm63Jy+NsfsCaydg/INgQvAuemzhdfmrqv5ZY/Q40TYJgQ/AuOJZvfe9yq8++w32sXoMsCDYE\n74LZxhF5p9xUw9iOku5F8yAYGAAEEweCiQPBxIFg4kAwcSCYOBBMHAgmDgQTB4KJA8HEgWDiQDBx\nIJg4EEwcCCYOBBMHgokDwcSBYOJAMHEgmDgQTBwIJg4EEweCiQPBxIFg4kAwcSCYOBBMHAgmDgQT\nB4KJA8HEgWDiQDBxIJg4EEwcCCYOBBMHgokDwcSBYOJAMHEgmDgQTBwIJg4EE8ej4N11etMAukgv\neMQnrObCrE6X7QwhGxA46QVb21jpjfsb7rkmhGxA4HgS3GcTY7HuIWQDAseD4LePlPydsVVfCSEb\nEDjpBQ89vXPeALY6b2EI2ST482DB0pDmI46Xo+iGLeVs3RrtqUhmrubN63PCmpA2nk6Tmj47qjuP\nY0BwoKQXfGDOWTlW9pkP1oeQTQIIDpT0gm8Zvqr2cO2akvEhZJMAggMlveDu/+PNvgLdqUggOFDS\nCx64iDdlgzRnkgSCAyW94MreZ48ZP7b4tPUhZJMAggPFw1H0kRWL5i1ccUR/LgIIDhSPvyb9Tm8W\nLYHgQPEouIfeLFoCwYECwcTxIvgfRUVZRUXHKtq/PJxz/tNaMpKCn//6JM5zWiZpP/jYgl/4VaCZ\nJJGCfz6tOsGGUVomaT9krmCxiz4Awf7wcRQNwSbgTXDjziZ3JwSbgBfBNRd3zav8xlZnNwSbgBfB\n37rzUO8j0y5zdkOwCXgRfEId681iXZ3dEGwCXgT3fzku+JWvOruDFrxmIacEgoPEi+BVBaO7XXPK\n687uoAVf+kJZgn4QHCSejqJrn529aJerN3DB4sKvoRAcJD7uTYJgE/BxbxIEm4CPe5Mg2AR83JsE\nwSbg494kCDYBH/cmQbAJ+Lg3CYJNwMff6IBgE4Bg4kAwcSCYOBBMHAgmDgQTB4KJA8HEgWDiQDBx\nIJg4EEwcCCYOBBMHgokDwcSBYOJAMHEgmDgQTBwIJg4EEweCiQPBxIFg4uAGcOLgBnDi4AZw4uAG\ncOLgBnDi4AZw4uAGcOL4KE4JwSbgozglBJuAj+KUEGwCPopTQrAJ+ChOCcEm4KM4JQSbgI/ilBBs\nAhlUVgeCdZBBhbEgWAcQTJy2FKdcKopTnjs12FQgWAfYgokDwcTBUTRxMqg4JQTrIIOKU0KwDjKo\nOCUE6yCDilNCsA4yqDglBOsgk4pTQrAGUJySOF5vPqtwXZIFwUbgVXDBDlcXBJtAesG52QmsrGxn\nAIJNIL3gjUPGVMdi+e/HnAEINgEPu+jGR4te1bqLLhU/P54KwRrw9Bn80cU3dtMo+GLR5EOwBrwd\nZB1dWFrr6oRgE8iEv9EBwRqBYOJAMHEgmDgQTBwIJg4EEweCiQPBxMl0wXu/Kr6oviCgudodmS44\nNso2BrQWCCYOBBMHgokDwcSBYOJAMHEgmDgQTBwIJg4EEweCiQPBxJGC71rd2OpFIdgEpOCZ/Xv8\n8E33H5Q9LhBsAs276OpffLPHra8fbsWiEGwCzYLrysblFw8pXKYeprU4JQRrRAp+eFjut5/YytjK\nL7pG6C9OCcEakYK//+I+3u5f6h6hvTglBGtECq5ftJr94dEG5QjtxSkhWCNS8G2D32PlX3OVRuIj\ntBenhGCNSMEnVcf/qXZtpAn0F6eEYI1IwX3fif9T3lc9RndxSgjWiBS8uOCuR+4ueC7FIM3FKSFY\nI8nz4E2zb5v1vnKE/uKUEKyRTChOCcEakYJXDCkqallYpQX6i1NCsEak4F73f1AVRzVCf3FKCNaI\nFNzzYMoR+otTQrBGpOD581L/Hqy9OCUEa0QKvij3C19J8RnM0VpWB4I1IgVXCVIO01oYC4I10nya\npCyN1IxN8N71nLmPBZQDBGtECk5RGkngLE75xiTO8OkB5QDBGpGCU5RGaga7aFORglOURmoGgk1F\nCk5RGqkZHEWbihScojRSEr3FKSFYI8mjaHVpJIHu4pQQrJFMKE4JwRqRgs8XqMfoLk4JwRqRgisq\nKsrLhi1Rj9FdnBKCNdJyFx07Vz1Gd3FKCNZIS8HvdksxSHNxSgjWSIvP4MEdJ6ccVq/4NRGCTeDY\nZ3BFRZXy54YPR96y7eKcLqX6SttBsEbSnyZdOH7GyffWbr/5WmcAgk1ACi7snsQ1osvuOusAY7vz\nnAEINgEpeMGlb+8sv/Txujr3fcCnbmxaHG8qip0BCDaB5BacuDa25kuqEQ/0WsfYx5NPfdYZgGAT\nkIJPT7yqq09TjWhavpWxzQ9VugIQbAJS8ML86c9Mz5/fqkUh2ASSR9Hrppbe+ebxrspyA8Em4PGi\nOxUQbAJeLrpLAQSbgMeL7lRAsAl4vOhOBQSbgMeL7lRAsAl4u+hOCQSbgJeL7lIAwSYgBZ+7ofWL\nQrAJSMGzx7n+xkpaINgEpOBh3XP6Hu/+YBUQbAJS8AeCVi0KwSbABefWMfb8561dFIJNgAu24oK7\nb2vtohBsAhBMHAgmjhD8VmVl7kuVle6rNo4HBJsAF1yQpFWLQrAJoPIZcSCYOBBMHAgmDgQTB4KJ\nA8HE8SgYxSlNJb1gFKc0mvSCUZzSaDwJRnFKc/EgGMUpTSa9YBSnNBovR9EoTmkwnk6TUJzSXNIL\nRnFKo0FxSuKkF4zilEaTXjCKUxpNesEoTmk0Ho6iUZzSZDz+moSyOqbiUTAKY5lKWwRvWsi59ecB\n5QDBGvEi2FmccnMZZ8qcgHKAYI1gF00cCCYOjqKJ400wilMaixfBKE5pMChOSRwvglGc0mC8CEZx\nSoPxIhjFKQ3G01E0ilOaixfBexLnSI2oXWgk6QVv6Nehz8uMbXONhGATSC/4oh83rC6shGBDSS+4\n62eMLRvcCMFmkl7wgCWMNV01Q4PgmhHDORCskfSC3zjxgk9ZbNDA4AW/c79oIVgjHo6idy7ex1j9\n4vuc/RBsAlH+jQ4IDgEIJg4EEweCiQPBxIFg4kAwcSCYOBBMHAgmjimCTxK/S5xR7XPOdocpgvNE\nM/k9n3O2OyCYOBBMHAgmDgQTB4KJA8HEgWDiQDBxIJg4EEwcCCYOBBMHgokTZe1CCA6BKGsXQnAI\nRFm7EIJDIMrahRAcAlHWLoTgEIiydiEEh0CUtQshOASirF0IwSEQZe1CCA6BKGsXQnAIRFm7EIJD\nIMrahRAcAlHWLoTgEIiydiEEh0CUVVcgOASirJsEwSEAwcRpS+3CNyZxhk/3OTcEh0BbtuDPqzlP\nPOZzbggOAeyiiYOjaOJEWbsQgkMgytqFEBwCUdYuhOAQiLJ2IQSHQJS1CyE4BKKsXQjBIRBl7UII\nDgHcfEYcCCZOesFVSZwBCDaB9IJHWl0LOc4ABJuAh130xDvU/RBsAh4Er5iv7odgE8BBFnEgmDgQ\nTBwIJg4EEweCiWOY4JLR4pLdtT6nbj8YJnjgK/yK3Sef9Dl1+8E0wWIXXQbBXoFg4kAwcSCYOBBM\nHAgmDgQTB4KJA8HEgWDiRCL4+Qt5vfaBEKyfSAQ/8jfe/AaC9QPBxIFg4kAwcSCYOBBMHAgmDgQT\nJ5KqKxAcHpFUXYHg8Iik6goEh0ckVVcgODwiqboCweERSdUVCA6PSKquQHB44DSJOGaeJs0/m18x\nMNzvNSXtADNPk6aILXjndW3NoP1g5mlSUvCo9ZytbU2kHWDmaZIU/J9z7ucMaGsi7YBwT5M2lHFu\nDkjwe9faVkCdgyuWc/7bmoXCPU26ata8BOdAcFv46xX81Xvg2tYsFHhZnZj4WFzfoApeVceb70Jw\nW3jxcd58OrY1CwVeGGv89fyPaPTvM5jza9F7ozit6aFH8BfFXEP2pn8iDaJqW/XB9EMzjswQfIv4\nhLjpj7xZ+WVhVjrsqUewDPYfyuf6cjFvLhgijsDuE5/73zuPvwlOG8HffyWXL+S8tocT29OST+Rn\nXbl4L2yS7wkR3C5WV7bR2wvnlaNybtU3SqxO7BQf1iXYWZxyidhgzpptG3VyV052F950zOFN52zR\n20k0WXKMCOZ0FA9l08G2ghy5SLYtmCWCnUTwBGUwOXVWNqejaJJrl1M7grYmSwazPQSz7E2nzgk6\ndeRN52zZyBXkiGByBTIh25gs+7NOvkD24BlBC2bqLdiB3A7lFvy02IKrfmALys9geZrkuDfpUvHu\nnSm24NfFFnxglG3MZNtXlbvsn8HjtvFm0XO82XSbLXj1Ht784iXerL3PFryskTezVvHmjZ/x5tB3\nbGPu/jdv/vIEbz4ptQW/v5U3v/k9bzZPsgVH1/Lml8t4s+5eW3C4OHqdtZI3y8VGUz/CNmbKu7wJ\nexftAIJNFqw6inYAwQYLVhandADBxgpOUZzSAQQbKzhFcUoHEGys4BTFKR3InMbX8Oa3i3mz5U5b\ncMw+3jz2Gm/Wz7AFRwrBP3mbN28+xJtDV9rG3LOBN0uf4k3sRltwwg7e/O5PvPnoDvvUn/Hm8Vd5\n868f2YKjxNQ//Sdv3prHm4YrbGOmfsCbZQvE1DfYghO38+b3z/Om+nZbsFS8rZ94hTfvTrcFLxfv\nrdlreLNyrpi6xDZm2n9485L44qH2etYKfBSndHDI1hxpVPXK5rAINtWrgg3ixT7aoArWi2OBxsPH\nWfuRI+1jam/4KE4JTMBHcUpgAj4uugMm4OOiO2ACPi66Aybg46I7YAI+LroDJuDjojtgAj4uurNR\nOShvnOsymKaZPU/4xofO3p0ju52/ybWCpUX5N+139I1MFKVePqDrRRsUM/Kgel1yEWVOx9ZXlete\nqQw+1fuEYa6C2Mlno0woGVQnJCdTJ9SciSoh2atMyBM+/kZHS470fqbmsjnO3hW9Nu4eP9LR2XTe\nI7umuC6Tq8pbGbt6sn3piVb8Oe08sWzvA/1cE4qgel1yEXVOzetrHJLtWqsMbum0Ytcdl6R4NuqE\nZFCdkJxMnVBzJqqEZK86IU8EJHjF2YytPMvZu3Vd094Ztzg6153dxOrfdw6dP46xTfbjuPl3dI07\nXHwBYw0d9rjG86B6XXIRdU7N63t0jPv1lMGd3d7ZN210imejTkgG1QnJydQJNWeiSkj2qhPyRECC\nF41hrDbH/aPxYuvkbY6uZ6+89azR250DfzYx/iJZjn10Ydzhvk8ZW32G4ufoRFC9LrmIOqfk+j4q\nqna/nsngAqtDQa17Sv5sUiXEg+qE5GQpXiSZiTKhZG+KhDwQkOB58bPkw9Zn7sCBac4d1sMdFmy+\nfYhz3Nr8dbVjrd32zkL+sdO0rPBviikTQfW65CKpcuLBo5e8ElNtMDxY1bP84L0liqh4NikSSgTV\nCcnJ1AnJoDoh2ZsyofQEJHjh2Pibs6PzC83NHyfeso6vwJ4cxtjB7JhzDb8u7DEnu9HexwXXXj24\nUjVlIphiXWIRdU7J4A1MJVgEH57AWH1n51XWyWejTEgG1QnJydQJJYPKhGSvOiFPBCR4eTFja850\n9j4UP2ra3sVh7eX4a1Dfsc4xdF8NY+XnOjoTDuvPm+FYQYugel1yEXVOMljarSDfKqhQB+fGDwgO\nuVYrn406IRlUJyQnUyckg+qEZK86IU8EdRR92pLPr3jQ2Vt56trYzc7rD+p7/LF2yjDn0I25az4e\nusDRmXC4eMC2OIqXlNtXrksuos5JBmt37Hg/a4fz+3UZ3NB9ee1droNW+WzUCcmgOiE5mTohGVQn\nJHvVCXkiIMGscsBJ49w/R/y2X/exrh3o2vNOHLnDNfTJUwofdB5/JBzeZyVwrUTuv5XrSi6izOnY\n+hR7xGTwxaJuJTUpnk2KhORTVT85OZn6RWrORHlQIHrVCXkhKMEgQ4Fg4kAwcSCYOBBMHAgmDgQT\nB4KJA8HEgWDiQDBxIJg4EEwcCCYOBBMHgokDwcSBYOJAMHEgmDgQTBwIJg4EEweCiQPBxIFg4kAw\ncSCYOBBMHAgmDgQTB4KJA8HEgWDiQDBxIJg4EEwcCCYOBBMHgokDwcSBYOJAMHEgmDgQTBwIJg4E\nEweCifN/UplSUs85VSoAAAAASUVORK5CYII=\n"
},
{
"output_type": "display_data",
"png": "iVBORw0KGgoAAAANSUhEUgAAAeAAAAHgCAMAAABKCk6nAAAC+lBMVEUAAAABAQECAgIDAwMEBAQF\nBQUGBgYHBwcICAgJCQkKCgoLCwsMDAwNDQ0ODg4PDw8QEBARERESEhITExMUFBQVFRUWFhYXFxcY\nGBgZGRkaGhobGxscHBwdHR0eHh4fHx8gICAhISEiIiIjIyMkJCQlJSUmJiYnJycoKCgpKSkqKior\nKyssLCwtLS0uLi4vLy8wMDAxMTEyMjIzMzM0NDQ1NTU2NjY3Nzc4ODg5OTk6Ojo7Ozs8PDw9PT0+\nPj4/Pz9AQEBBQUFCQkJDQ0NERERFRUVGRkZHR0dISEhJSUlKSkpLS0tMTExNTU1OTk5PT09QUFBR\nUVFSUlJTU1NUVFRVVVVWVlZXV1dYWFhZWVlaWlpbW1tcXFxdXV1eXl5fX19gYGBhYWFiYmJjY2Nk\nZGRlZWVmZmZnZ2doaGhpaWlqampra2tsbGxtbW1vb29wcHBxcXFycnJzc3N0dHR1dXV2dnZ3d3d4\neHh5eXl6enp7e3t8fHx9fX1+fn5/f3+AgICBgYGCgoKDg4OEhISFhYWGhoaHh4eIiIiJiYmKioqL\ni4uMjIyNjY2Ojo6Pj4+QkJCRkZGSkpKTk5OUlJSVlZWWlpaXl5eYmJiZmZmampqbm5ucnJydnZ2e\nnp6fn5+goKChoaGioqKjo6OkpKSlpaWmpqanp6eoqKipqamqqqqrq6usrKytra2urq6wsLCxsbGy\nsrKzs7O0tLS1tbW2tra3t7e4uLi5ubm6urq7u7u8vLy9vb2+vr6/v7/AwMDBwcHCwsLDw8PExMTF\nxcXGxsbHx8fIyMjJycnKysrLy8vMzMzNzc3Ozs7Pz8/Q0NDR0dHS0tLT09PU1NTV1dXW1tbX19fY\n2NjZ2dna2trb29vc3Nzd3d3e3t7f39/g4ODh4eHi4uLj4+Pk5OTl5eXm5ubn5+fo6Ojp6enq6urr\n6+vs7Ozt7e3u7u7v7+/w8PDx8fHy8vLz8/P09PT19fX29vb39/f4+Pj5+fn6+vr7+/v8/Pz9/f3+\n/v7///8vXNKMAAAUxUlEQVR4nO2dC3wV1Z3Hh1weQoiBRkshIFXAUK3y2hZQS4LYRVTwwQJ10UIR\nnyhQEKzuqi0qFFQQKUqhoC3dUqhiwV15CSw1EYKurEiwGKi8NVleIgQIOZ/P3nvOuZGZOZc59zG5\nc/739/185Mj5nznzn3yZycydufO3GCCNle4EgL9AMHEgmDgQTBwIJg4EEweCiQPBxIFg4kAwcSCY\nOBBMHAgmDgQTB4KJA8HEgWDiQDBxIJg4EEwcCCYOBBMHgomTGYLPWNYz4Wao9YTO6JNWyOd86pCM\nEdz0IATTJSzYegCC6RIW3DK0jQs+Pq5DdpfXatgBK6+425Id1mUv5l8yo6RLdmE5Y6t6NG12/QcQ\nbB5hwXOtWyKCa/pYrW+6wJoRFpzdxgoLthoWhP9rl23dznY3DhV1t9rWQLBxhAXvGGCtCQteb+Uf\nZmusZmcPWNbkyqqw4DLWy7qz5l2rHVvX91fscH2rAoKNIyK4LNTlTuuJl60Hw3/Ptz4/YF1wlrEd\nVn7kV/Pc8P+0ZWzLhN5NLesABBtHRDB7yMq1npjJBbexdoZ/BzMmvA61FvD/+Vuo+di/XATBBsIF\nf5FjWU+ss1ofYWut3LMKweOsiWx3PQg2EC6YPRcWXNPbanNzY36S5RL8vNX0lu9Y1j4INg4h+ETr\n8GXSsTHtmsjLJOYQfHxoTruXrrEWQDAwBwgmDgQTB4KJA8HEgWDiQDBxIJg4EEwcCCYOBBMHgokD\nwcSBYOJAMHEgmDgQTBwIJg4EEweCiQPBxIFg4kAwcSCYOBBMHAgmDgQTB4KJA8HEgWDiQDBxIJg4\nEEwcCCYOBBMHgokDwcSBYOJAMHEgmDgQTBwIJg4EEweCiaMp+MvD/qYB/MJb8PaiLbt7hBoU7amD\nbEDK8RbcfWzVHaNOVo25uQ6yASnHW3DTg6z93xmryKmDbEDK8RZ84ws1977C2MKudZANSDnegvd1\n7Tggq6iwxcY6yAakHI2z6JpNC6bOWlblfy7AB7Quk2qOnvU7D+AT3oK/frZDQyvU/mnswkbiLXjY\nDesqT1du6D+iDrIBKcdbcO4+3hzL8zsV4AfegjvP483iLnoTTu8mKE4mK5AyvAWXtu04aMTgK1tt\n1ptw2D9489uFyWQFUobGWfSZ1fOmzFl9RnNCCA4WmneTXtOeEIKDhabgFtoTQnCwgGDi6Aj+74KC\nrIKCgtq/r72P85MlqsEQHCwS2YOPlnNmTleNhOBgkcQh+s+/UY2E4GCRxFk0BJuAnuDq/TXuTgg2\nAR3Be4uaNCu9bqezG4JNQEfwjx8+2fbM+D7Obgg2AR3BjQ+ztqyiibMbgk1AR/DVy8KCl3/f2Q3B\nJqAjeF3ewJw7Ll7h7IZgE9A6i65cMGneAVcvBJtAEl8+g2ATgGDieAsui+IMQLAJeAvuZzXJ5zgD\nEGwCGofokaPU/RBsAhqCV09T90OwCeAkizgQTBwIJg4EEweCiQPBxIFg4kAwcSCYOBBMHAgmDgQT\nB4KJA8HEgWDiQDBxIJg4EEwcCCYOBBMHgokDwcSBYOJAMHEgmDgQTBwIJg4EEyeJ6qMQbALegvse\nZHuvyWrQZ78zAMEm4C3Y2sWG3HX81M/vcAYg2AS0BF+2nbGKXGcAgk1AQ/B7Z/q/w9i6y50BCDYB\nb8G92jRq1omtbzbHGYBgE9A5iz61o5ht2uDqhmATSKL6KASbQBLVRyHYBJKoPgrBJpBE9VEINoEk\nqo9CsAkkUX0Ugk0gieqjEGwCqJtEHFQ+Iw4EEyeR6qPvDOL0eEw1GIKDRSJ78Fey+ugM1UgIDhY4\nRBMHZ9HEQfVR4qD6KHFQfZQ4qD5KHFQfJQ6qjxIH1UeJo3kdfATfTTIUb8Hbeg+svKVRqPc+ZwCC\nTUDjwfdR4/MmVJ14+FZnAIJNwFtw44PH6p1grOJCZwCCTcBb8CWrS6wSxlZ+zxmAYBPwFjy7wYWz\nW94zPPdNZwCCTUDjLHrHflY2dfInrn4INgG8o4M4EEwcCCYOBBMHgokDwcSBYOJAMHEgmDgQTBwI\nJg4EEweCiQPBxIFg4kAwcSCYOBBMHAgmDgQTB4KJA8HEgWDiQDBxIJg4EEwcCCYOqo8SRwp+ZH11\nrBGoPmo0UvCTV7d4cI37rf18BKqPmkztIbr8hR+1uHfFacUIVB81mVrBhxcPb35lz/yl7hGoPmoy\nUvDUwux/nrmTsbXfcY1A9VGjkYJ/9sYx3h53vaiBofqo0UjBVfPWsz9MPxVjEKqPmosUfH+3j1jx\nD1z1JyOg+qjRSMHfKg//Ue46UY6A6qNGIwW3ez/8R3E71QhUHzUaKXhR3iPPj8n7vWoEqo8aTfQ6\nePuk+5/aohyB6qNGg+qjxJGCV/csKDi3ep0T1E0yFSn4kokfl4WJOQyVz0xFCm554vzDINhUpOBp\nU2LeD3ZXH10uq48+qhoMwcFCCr42+8LLz/c72LYHnzrEmf+yaiQEBwspuEwQcxgO0aZSe5mkrD9Z\nC86iTUUKjlF/MgqqjxqLFByj/qQA1UcNRgqOUX9SgOqjBiMFx6g/KUD1UYORgmPUnxSg+qjBRM+i\n1fUnBag+ajCoPkocKbi7IOawuL+b9JtJmzn/l3yKIBmk4JKSkuLFhUtUI7YXbdndI9SgaI8zcF7B\nw/pPjDByVMoyBQlx7iG64irViO5jq+4YdbJqzM3OwHkF3y0O0dseTDpDkBTnCv4wRzWi6UHW/u9h\n+64gBJvAOb+Du9UfrRpx4ws1977C2MKuzgAEm8A3v4NLSsqUtxv2de04IKuosMVGZwCCTUDjMqlm\n04Kps5a5vtgAwUYgBefnRtFfFIJNQAqeff17+4uvf+nwYcX1biwg2ASie3Dk+yl7W8e1KASbgBTc\nZn34j/Wt4loUgk1ACp7T/LG5jzWfFteiEGwC0bPoTeOGPLzmfE9luYFgE9B86E4FBJuA3kN3SiDY\nBHQeuosBBJuAzkN3MYBgE9B56C4GEGwCOg/dxQCCTUDnobsYQLAJSMFXbY1/UQg2ASl40nD37UAv\nINgEpODC3Ibtzvf9YBUQbAJS8MeCuBaFYBPggrMPM/bHr+JdFIJNgAu2woJzd8W7KASbAAQTB4KJ\nIwS/W1qa/VZpaWlci0KwCXDBeVHiWhSCTcCv0nYQHBAgmDgQTBwIJo5f1UchOCB4C06s+igEBwRv\nwYlVH4XggKAlOIHqoxAcEDQEJ1R9FIIDgrfgxKqPQnBA0DmLTqT6KAQHBK3LpASqj0JwQPAWnFj1\nUQgOCN6CE6s+CsEBwVtwYtVHITggeAtOrPooBAcEb8GJVR+F4IDgV/VRCA4ImneT4q6bBMEBQVNw\n3JXPIDggQDBxdAQ7q48uu4Fz1QTVYAgOFtiDiQPBxMFZNHH0BMdffRSCA4KO4ESqj0JwQNARnEj1\nUQgOCDqCE6k+CsEBQUdwItVHITgg6AhOpPooBAcEv6qPQnBA0P3yWQnqJpmJruA8V/FRCDYCb8HZ\noQhWVsgZgGAT8Ba8reeg8oqK5lsqnAEINgGNQ3T19IK3cYg2Fa3fwZ8V3ZUDwWaid5J1ds6QSlcn\nBJsA3tFBHAgmDgQTB4KJA8HEgWDiQDBxIJg4EEwcCCYOBBMHgokDwcSBYOIEXfDpboLeSc6TsQRd\n8Nc3ibYoyXkyFggmjimCO9zHefBIkvNlHKYIblbOGf5RkvNlHIEVvOcBvsuOiAoWzWgIjpPACl4x\nnu+ym+2CR8xZxXE/AgjUBFfws7ypsAvuPHpKhKHKVQMFpgkWh+jFsxKdNuOAYOJAMHF8FrzuIvHa\nw3+Je3YITg0+C17+gOiN/4MoCE4NPlcfheB043P1UQhONz5XH4XgdONz9VEITjc+Vx9NQPCC6/h5\n95UQnBJ8rj6agOAn1/PmOQhOCT5XH/VJ8MwfT+Qs1Z82U/G5+qhPgsdO3BzhnTv1p81UfK4+6pdg\ncYjeD8Ge+Fx9FILTjc/VRyE43fhcfRSC043P1UchON34XH0UgtONz3WTIDjd+Fz5DILTDQQTJ5Hq\no2/J6qOPqgZDcLDAHkwcCCaO0WfRn+aIb4fju6Wx8bn6aDyCJwpbrfQFf/QT/dkzFZ+rj8Yj+Hpx\nT7IXBKcSn6uPQnC68bn6KASnG5+rj0JwuvG5+igEpxufq49CcLrREXwoco1UnVBpOwhON96Ct15R\n77JljO1yjYRgE/AWfO2/n1qfXwrBhuItuMlRxpZ2qw6y4G7i1TtrFB+3ZTzegjstYazmtseDLDiH\nv3pnSrd/aKwl0/AWvLJpjy9YRZfOARYsg8N3aawl09A4i96/6BhjVYsmOPsh2AQC9I4OCPYDCCZO\nEASvXMzpDME+EATBHedwmkCwDwRBsAw2h2AfqCPBncVReJlyIgj2kToSfLU4CndWvE8Ngn2lbg/R\nt/kqeNBC8ZHl8bg2gziUBHd4gL+Zpd8f49oM4lAS3H4Xb+b+QXsTMgAIJg4EEweCiQPBxIFg4kAw\ncSCYOOkU/Mn3xLsgUiz4Fz8cxIHnCOkU/P5EWzBVgoe9eijCpvvj2hqqUBT8e95sh+AIEEycuhXc\nZwG/LbxE3O+B4DqgbgW3nM5vC982g7+w/TUI9p86FiwO0QP68urehRDsP2kRfMtfefM7XwW/d7G4\nCPvX+LaKGnQF/9f9ttkzFQgmjpbglJXVgeA6x1twKsvqQHCd4y04lWV1ILjO8RacyrI66RB8sXgF\n5o9cR6DMwFtwKsvqpHEPvv2QcuvI4y04lWV1ILjO0TiLTmFZnTQK7v9heYSdqusBymTMZVLrW/lT\nAJ1cb2QkTsZcJrUSh+jn/6rcSLpkzGUSBMfC9Msku+BxA0Xt8GJV7hTJmMskKbj/i/xW9Ouul0JR\nJWMuk6KC3+LN25eKe4lPqjaBFBlzmWQXPF/uwe3EAfs51ZYEjrM7+ZVe+bF4FtK8XZiiukmBE/xP\n/IC9+fLJ/GWXU4+otigwLL+GPwhz913xLKQpOEWVzwInOLrqFfzdD33785/gyJfFqyA+U21fGnnj\nJd58MTiehSCYr7qaN4V/4ofAib34cXtsgbhN8V1xGP/pL/lePmm+sL9S7PsfiKNmmWjKK/gz94fs\nv89Oi07l13biwTfBzuqjS8SGd5hkG3VRE07oAt7Ub8ibRiHR20A0WXKMCDasL/4qm3q2CRrKRUK2\nYJYINhDBxspgylcd4sgmK0u2tt5oYw82tA+1BxvY/hZtGjTiyCYkm4a8qS/zujTVgpl6D3Ygdwa5\nB/9W7MFlygff5acNjuei5bsqZe3CFWIP/vom25jRYg9eLPbgA/Y3vsvXKM2zP1Vpv9nwgtiDN9r3\n4D5iD35qHW9WPsObkzfaxoz5H978ZSZvDg6xBX8m6sL97nXefHqfLTiwkjcvLuXNpkdtwRvE3v7U\nWt6sEjtNVV/bmLEf8qauD9EOINhkwaqzaAcQbLBgZfVRBxBsrOAY1UcdQLCxgmNUH3UAwcYKjlF9\n1IHMacRe3sxfxJsdD9uCg8SnbDP+kzebH7cF+wnBv3yPN2t+zZuTt9rG/Hwrb958lTcVd9mC9+zh\nzWv/wZvPRtlXfZQ3L73Nmw9+YQveJFb9q7/x5t0pvDk1wDZm3Me8WTpbrHqoLThyN29eF29RLH/I\nFhwi/lnPXM6bDx+zBW8W/7YmbeDN2sli1f1tY8b/L2/eEh88VMb1ZRwdwTGqjzo4aWvOVKt6ZXNa\nBGuqVMFT4od99pQqWCXOBapPn2f2M2cyY9V66AiOUX0UmIDWWbS6+igwAS3B6ofugAl4C4750B0w\ngSQeugMmkMRDd8AEknjoDphAEg/dARNI4qE7G6Vdmg0/4eysebJl4+s+cfbu75fTfbtrgjcLmt/t\nLJfSryz8x6pOTa7dqlgjD6rnkosoc/pmvrJs96Qy+GrbxoVlsbZGmVA0qE5IrkydUG0mqoRkrzIh\nLTRvF3pxpu3cvX2edfauvmTblyP6OTpruj5/YGyRc2hZs7UVt4+2Lz3SCm/T/qaLj/zbFa4ViqB6\nLrmIOqfa+ap7hlyzyuCOBqsPjOodY2vUCcmgOiG5MnVCtZmoEpK96oS0SJHg1R0ZW9vB2btzU82R\nx4c5Ojd1rGFVW5xDpw1nbLv9PG7aqCZhh4t6MHaqnuu7nyKonksuos6pdr7pg9w/Txncn/P+sfED\nY2yNOiEZVCckV6ZOqDYTVUKyV52QFikSPG8QY5UN3TeNF1kX7XJ0Lbj13g4DdzsHPjMy/EOyHMfo\n/LDDY18wtv5Sxe3oSFA9l1xEnVN0vs8Kyt0/z2hwtlUvr9K9Sr41sRLiQXVCcmUxfkgyE2VC0d4Y\nCWmQIsFTwlfJp62j7sDX450HrKn1Zn/6UE/nuI3NN1UOtr60d+bzXzs1S/NV3xiLBNVzyUVi5cSD\nZ3svr1DtMDxY1rL4xKP9FVGxNTESigTVCcmVqROSQXVCsjdmQt6kSPCcweF/nPWdH2h++nnkn6zj\nI7BZhYydCFU4Z3glv8WzoWp7HxdceXu3UtUqI8EYc4lF1DlFg0OZSrAITr2HsapGzqfgo1ujTEgG\n1QnJlakTigaVCcledUJapEjwqisZ29De2fvr8FnT7gsc1paFfwZV9Z0PCR/by1jxVY7OiMOqro87\nJjgnqJ5LLqLOSQaH5OQ1t/JK1MHJ4ROCk65p5daoE5JBdUJyZeqEZFCdkOxVJ6RFqs6iWy35asDT\nzt7Sb2+s+Knz+YOqFgsrxxY6h27L3vB5r9mOzojDRZ12hVH8SLl95VxyEXVOMli5Z8+WrD3Oz9dl\ncGvuqspHXCetcmvUCcmgOiG5MnVCMqhOSPaqE9IiRYJZaadvDXffjph/Re5g1wF0Y9em/fa4hs66\nOP9p5/lHxOEEK4JrEnn8Vs4VXUSZ0zfzKY6I0eAbBTn998bYmhgJyU1Vb5xcmfqHVJuJ8qRA9KoT\n0iFVgkFAgWDiQDBxIJg4EEwcCCYOBBMHgokDwcSBYOJAMHEgmDgQTBwIJg4EEweCiQPBxIFg4kAw\ncSCYOBBMHAgmDgQTB4KJA8HEgWDiQDBxIJg4EEwcCCYOBBMHgokDwcSBYOJAMHEgmDgQTBwIJg4E\nEweCiQPBxIFg4kAwcSCYOBBMHAgmDgQTB4KJ8/+y4MwDC5QnbQAAAABJRU5ErkJggg==\n"
}
],
"prompt_number": 22
},
{
"cell_type": "code",
"collapsed": false,
"input": [],
"language": "python",
"metadata": {},
"outputs": []
}
],
"metadata": {}
}
]
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment