Korean Flag in Python

Korean Flag.ipynb

{
 "metadata": {
  "name": "",
  "signature": "sha256:511d0c108a8543da73b1b6f822817a270cec46679b334b547986444a85818f92"
 },
 "nbformat": 3,
 "nbformat_minor": 0,
 "worksheets": [
  {
   "cells": [
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "from utils3pt import *\n",
      "from SimPEG import *"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 8
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "C = Vector(1.5,1,0)\n",
      "V1 = Vector(0,2,0)\n",
      "V2 = Vector(3,0,1)\n",
      "P = Plane(C,V1,V2)\n",
      "Cd = (V1 - V0)\n",
      "Cd.length = 0.25\n",
      "C1 = V0 + Cd\n",
      "C2 = V0 - Cd"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 37
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "M = Mesh.TensorMesh([[(0.01,300)],[(0.01,200)],[0]])"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 38
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "Cd.length = 1./2 + 1./4 + 1./24\n",
      "ltb1 = C + Cd\n",
      "ltP = Plane(Cd,ltb1)"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 54
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "def getBar(D, dash, i=None):\n",
      "    if i is None:\n",
      "        D = Vector(D)\n",
      "        return getBar(D,dash=dash[0],i=0) | getBar(D,dash=dash[1],i=1) | getBar(D,dash=dash[2],i=2)\n",
      "    D.length = 1./2 + 1./4 + 1./24 + i*(3./24)\n",
      "    PN = Plane(D, C + D)\n",
      "    P  = Plane(C + Vector(0,0,1),C,C+D)\n",
      "    if dash:\n",
      "        return (np.abs(P.dist(Vector(M.gridCC))) < 0.25) & (np.abs(PN.dist(Vector(M.gridCC))) < 1./24) & (np.abs(P.dist(Vector(M.gridCC))) > 0.5/24)\n",
      "    else:\n",
      "        return (np.abs(P.dist(Vector(M.gridCC))) < 0.25) & (np.abs(PN.dist(Vector(M.gridCC))) < 1./24)\n",
      "\n",
      "black = (getBar(V1 - V0, [False,False,False]) | \n",
      "            getBar(V0 - V1, [True,True,True]) |\n",
      "            getBar(C - Vector(), [True,False,True]) |\n",
      "            getBar(Vector()-C, [False,True,False]) )\n",
      "\n",
      "blue = ((Vector(M.gridCC) - C).length < 0.5)\n",
      "\n",
      "red = ((P.dist(Vector(M.gridCC)) > 0) & ((Vector(M.gridCC) - C).length < 0.5) | ((Vector(M.gridCC) - C1).length < 0.25)) & ((Vector(M.gridCC) - C2).length > 0.25)"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 76
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "kf = np.zeros(M.nC)\n",
      "kf[blue] = 1\n",
      "kf[red] = 2\n",
      "kf[black] = 3\n",
      "M.plotImage(kf)"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "metadata": {},
       "output_type": "pyout",
       "prompt_number": 79,
       "text": [
        "<matplotlib.collections.QuadMesh at 0x10d9e3210>"
       ]
      },
      {
       "metadata": {},
       "output_type": "display_data",
       "png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAEKCAYAAADpfBXhAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAGyBJREFUeJzt3X9s1fW9x/HXwTKWAuOH0yJtHWpxKwhtEdctV+K4bmN0\nocPAFkhAruDWMYQwbzKvNy6cGkbc/QPjhiO6ZBjjFdhcDHW2JGZX0AGVq6vLXHUUR6+nFWoKloGK\n1PK9f2hrT3+cX98f5/P9fJ+PptFyvuecz7ff9/fV9/l8P6eNOY7jCABglTH5HgAAwHuEOwBYiHAH\nAAsR7gBgIcIdACxEuAOAhQh3ALAQ4Y5IefLJJzV//nxNnDhR06dPV01NjQ4dOiRJOnbsmL773e/q\niiuu0OTJk1VRUaEHH3xQly5dyvOogewR7oiM7du368c//rHuu+8+vfPOO0okEtqwYYMaGhr05ptv\nqrq6Wl/4whf02muvqaenR7/73e/0yiuv6Ny5c/keOpC1GO9QRRScPXtWJSUleuyxx7Rs2bJht69a\ntUpnz57VM888k4fRAd6jc0ckHDlyRBcuXNBtt9024u1//OMftXz58oBHBfiHcEcknD59Wp///Oc1\nZszIJX/69GldddVVAY8K8A/hjki4/PLL1d3dPerF0csvv1xvv/12wKMC/EO4IxK++tWvaty4cXr6\n6adHvP3rX/+6fv/73wc8KsA/hDsiYdKkSbr//vu1YcMG7du3T++//756e3vV1NSke+65R/X19Tp8\n+LB+8pOfqKurS5J0/PhxrV69WmfPns3z6IHsEe6IjLvvvlvbt2/X1q1bdeWVV+rqq6/Wr371K912\n22269tprdeTIEbW3t2v27NmaPHmyli9frptuukkTJ07M99CBrKVcCplIJHT77bfrnXfeUSwW0w9+\n8ANt2rRp2HabNm1SU1OTCgsL9dhjj6mqqsrXQQMAUitIdePYsWP14IMPqrKyUufPn9eNN96ob3zj\nGyovLx/YprGxUcePH1dbW5teeuklrV+/Xs3Nzb4PHAAwupTTMtOmTVNlZaUkacKECSovLx+2oqCh\noUFr1qyRJFVXV6unp2dgzhIAkB8Zz7m3t7erpaVF1dXVSf/e2dmp0tLSga9LSkrU0dHh3QgBAFnL\nKNzPnz+v5cuX66GHHtKECROG3T502j4Wi3kzOgBATlLOuUtSb2+vli1bplWrVmnp0qXDbi8uLlYi\nkRj4uqOjQ8XFxcO2Kysr05tvvulyuAAQLRUVFXr11Vezvl/KcHccR+vWrdOsWbO0efPmEbepra3V\njh07tGLFCjU3N2vy5MkqKioatt3HwR7PeoDh8bykhfkehE9s3jeJ/Qs7e/fPcbbkPBOSMtwPHTqk\nJ554QnPnzh1Y3rht2za99dZbkqS6ujrV1NSosbFRZWVlGj9+vHbt2pXTQAAA3kkZ7jfffHNGf6hg\nx44dng0IAOAe71D1zIx8D8BHM/I9AJ/NyPcAfDYj3wPw2Yx8D8BIhLtnrsn3AHxk875J7F/Y2b5/\nuSHcAcBChDsAWIhwBwALEe4AYCHCHQAsRLgDgIUIdwCwEOEOABYi3AHAQoQ7AFiIcAcACxHuAGAh\nwh0ALES4A4CFCHcAsBDhHiJxq/8GLWxCreZfyj+zB/P0nzScPDANNWkWOncAsBDhHhJDuyK6JJhk\npHqkRvMr5jiOE8gTxWISBzsn6U4STiLkSya1R33mznG2KBaLKZeYpnO3QPyTDyBI1JzZCHeDEdqw\nATWcH4S7obI5IfghgHzItkYRLMLdYpxQ8AJ1FE6Eu4G86IhYDw8vpKsjundzsVrGMF4Gu9vHR3Rl\nW0OEvD9YLYOMcWIhHWrEDoQ7AFiIaRmD+D0l42Zb2M9t7VB73mNaxgJBBzvghhfhTr36i3CPKE4s\nDEVN2IVwN4AXL4eDfEmN8PNrWs+LWoQ3CPc8MymU44M+YCevjy/XfsxFuEdMpj8gOOHsM/SYcu3G\nbqyWySMTVidk8hic7OGW7vj52X3TwbvDapmQCdscOSddOAU5xUa4m4dwBwALEe4RwNwqUhmtww+6\nPqhHbzHnHjBT5jF5uRwdYbi2Q42NjDn3kAg6lIGgZfsqgKW3/iHcLUbXDinYzj2b56G2/MW0TEBM\nmUoh2KMr6JDPtX7670f9MS1jPKZjEEW51Ozg+9Ddu0O4W2i0rouuPdq86Nz9qguC3HuEu89MmErJ\n9v6cZPYyoR6z3YZ6zE3acF+7dq2Kioo0Z86cEW8/cOCAJk2apKqqKlVVVWnr1q2eDzKsTO6UgKCk\nquFM65vzIHsF6Ta44447tHHjRt1+++2jbnPLLbeooaHB04EBAHKXtnNfsGCBpkyZknKbgBbcIA0v\npmRgPy7wR4PrOfdYLKbDhw+roqJCNTU1am1t9WJcoWfCPDnBjtGYMvdOjfrHdbjPmzdPiURCf/nL\nX7Rx40YtXbo0xdbPD/o84fapjeTFqhSKGDYi4DN1Qv05GY/Hc36UtHPu6UycOHHg/xcvXqwf/ehH\nOnPmjKZOnTrC1gvdPh2G4IRJz3k77ur+senu7m+qbBqRVK8Oba+f4F3zyacUj29RfX19To+S0TtU\n29vbtWTJEv31r38ddltXV5euvPJKxWIxHT16VN/73vfU3t4+/Iki8A7VbIvcr6mUqAa72xDPhk2B\nT92ay807VNN27itXrtTBgwfV3d2t0tJS1dfXq7e3V5JUV1enp556Sjt37lRBQYEKCwu1Z8+e7PfA\nAlykCl6QYZ7Jc9sU+CPxqj69eBXAK4b00ob77t27U96+YcMGbdiwwbMBITNR7n7yGeqp9I8rbCHf\nXwep6iFdrRC25uEXh3nAlJektga7qWGeibAG/WhfZ3v/XLa3tY5zwS8OyyOmY/zjvB0PdbBL4d6H\nIOqScPcP4Q4AFnK9FBLBctvp5LJ90MLa6aYyeJ9MnqrJ5oLmSNtnMn+PYNC5u2DCPDnBHj5h38fB\nNWRCXXMxd2R07jlirt07YQ+7XIRtZY3ftTpSQGf7KoCQT0bnHhK2djdRDPbBwrD/mSyDzOV+8BdL\nIXNgwlRK2FcUhCHUgmZiF2/C1KNXjxFGLIUMWLYXnXK5r80I9pHxfXEf7pxfnyLccxTENMdozxHm\nLoYAS820dfFhuLYUH/SBTxHuLmVT0FF/eWpSaJnOpO+VF9d2/Kp902rcJIQ7AFiIcPfA4G7Fi87F\nRiZ1omER5e9Zuu6daZj0CHcPeVVstk3JRDmk3DLlexf0Ci8vxxJVLIUMQFTn2k0JJhuYtEzSbS16\n8YMiKlgKabCoFibB7q0ofz+jeg65RecekEwK1JauPcpB5DdTOvgguncTajnf3HTuhHuAUhWrLS9f\nCXb/2RLwufx71Pj6N1ThncFX+gEMx7nhHTr3PEoV9mGbipHo2oNkY/duSh2bhGkZyxDsyJQJIR/G\neg0LVssAAJIQ7gBgIcLdMGF8icuUTP6Y8L0PY81GAeFumEyL35iT5G9b8j2CyAtTwBtTtxHAUkgD\nhWn1gDMllu8hIATCUMu2oXM3WPyTj5H+3Qh07cYwuXs3pl4jhs49BAZ38iadKHTtGMrUWo0i1rkj\ne5907IS7eUxY9w7vsM4dAJCEcAcACxHuAGAhwh3Z+dsWOVNizLcbyoRVMzAD4Q4AFiLckTnWtYcC\n3Tskwh0ArES4A4CFCHcAsBDhjsww3x4qzLuD3y0DuBR7N/u3hrOUFH6jcwcAC9G5I2v9nWqUu89c\nuvVU94/y9xL+INyBLLgN9XSPS8jDK4Q7kAG/Qj3V8xD0cINwR3oRXSkTVKCnen43Ad+/Yobf8R5N\nXFBFzvIdfn4yZd9i7zrGjAXhQrgDQ5gYpiaOCWZLG+5r165VUVGR5syZM+o2mzZt0syZM1VRUaGW\nlhZPBwgEyeQQNXlsME/acL/jjju0f//+UW9vbGzU8ePH1dbWpkcffVTr16/3dIAwm02BE4Z9CcMY\nYYa04b5gwQJNmTJl1NsbGhq0Zs0aSVJ1dbV6enrU1dXl3QiBAIQpNMM0VuSP6zn3zs5OlZaWDnxd\nUlKijo4Otw8LAHDBk6WQjpPcScRioy3fen7Q/8+QdI0XT488c7tkL9/ohGGWE5LaJUnxeO616bpz\nLy4uViKRGPi6o6NDxcXFo2y9cNAnwW4TAjJYLJG02TXqz8l4PJ7zo7gO99raWj3++OOSpObmZk2e\nPFlFRUVuHxYIRNgDMuzjh3/STsusXLlSBw8eVHd3t0pLS1VfX6/e3l5JUl1dnWpqatTY2KiysjKN\nHz9eu3bt8n3QMFOYpmcIRdgu5gydMPfriWIxSfEgngpey/LXD5ge8DYGe6rvOb9+ILwcZ4tisdiw\n65qZ4B2q8JzJ4Wny2AAvEe7whYkhauKYAL8Q7vBN/4oOE0LVhDH4yfb9Q/YIdwCwEOGO9GbXu36I\nfHWWprxyyIfY9DgXUyOMP9aBwAT5p+SiGuhAP8IdgfMr5Al04FOEO/JmpDDONvAJ9E+F6U1k8B/h\nDqMQ1oA3uKCKzHhwURXB4UIqCHcAsBDhDgAWItwBwEKEOwBYiHBH5rioGgpcTIVEuAOAlQh3ZIfu\n3VjOlBhdOwYQ7gBgIcIdACxEuIdI3JS/QTu7nukZA5k0JWNMrUYYfyA7BIaeKMacOFn+4Wz4x5T5\n9pFq05h6DSH+QLal4p98jPTvRqB7N4apwZ7q3+EvfiukgTgZkA1nSkwxw2umv6ap7eDQuQOAhQh3\nw2Ta2RjTATE1k3cmT8nkuh3cI9wNk03xG3OiEPD5Y8D3PpQ1GwGEOwBYiHDPIy9WwxjTCRnQQUaK\nIe818KJWRzsP4A7r3PMg00DPtuCNOEFY++4/A0Jd8i7Yc33MKHCzzp2lkAGicIHUWDLpHTr3gGRS\nrF68CcSIk4Lu3T+WdO1uzococdO5E+4+82pqhYBHlILdi/vYgF8/AABIQrj7LKodhyldpjUi/P2M\n7DnkEtMyAfFjZYFXz+c7pmjcMSTYvZhiDG0N5wnTMobwqhitC3dDwimUDPneBR3sXo4lqgh3Dwx+\nE4YXb0yykiEhFSoR/p6la3B441N6hLtL2XTZkZ6akSIdVlkz6HuVbc0FWfvG1bhBCPccBdE5+Hmi\n5I1BoWUkQ36tQL+gG5JcxAd94FOEew4yLSI6jVEYFF5G4fvi+gcE59enWC2TIxOmWELdwfdjJY3x\noZ5J7fg5jWhFneeI1TIR4EW4G1n4hgeb70Kw/+nqJl/TMUiNcAcACxXkewBhlU0nnOrCKN2NkrvX\nqEzThKBjH2xwnfpRs15cn+JcSkbn7oIJ84nWnQAhC72chHwfMwn6oK8nGV/XeUDnHjIjFXL/125f\nSRijP/xs6uJDEuiZdujMs5uPzt2loLv3SDFszXdOQrwPQdSltavADJB2KeT+/fu1efNm9fX16c47\n79Q999yTdPuBAwf0ne98R9dee60kadmyZbrvvvuGP5FlSyGHMmGKxqvHMJ7pHb2FYe5mOWQ220ai\nfrPg21LIvr4+3XXXXdq/f79aW1u1e/duvf7668O2u+WWW9TS0qKWlpYRgx3ecnsS5rKtUUzthk0d\nVxqZTNN5He7wX8o596NHj6qsrEwzZsyQJK1YsUL79u1TeXl50nYBvQ/KaNmunhn830wfg5NniKFB\nGmRHH8IQdyuTaztuLrZGojEJUMrOvbOzU6WlpQNfl5SUqLOzM2mbWCymw4cPq6KiQjU1NWptbfVn\npCHgxcqVoAvfqpOkv3Me+unH41rCz1UpJtR3lKXs3D+eJ09t3rx5SiQSKiwsVFNTk5YuXapjx46N\nsvXzg/5/hqRrMh5oFGX6Uni0kyibVxJWnzAWhbGXvJojz6STt7q+PHdCUrskKR7PfVYkZedeXFys\nRCIx8HUikVBJSUnSNhMnTlRhYaEkafHixert7dWZM2dGecSFgz4JdgAY7hr152Q8Hs/5UVKG+/z5\n89XW1qb29nZdvHhRe/fuVW1tbdI2XV1dA3PuR48eleM4mjp1as4DCjuvXrpyAQth4eYVZqb3z2Xb\nqEu7FLKpqWlgKeS6det077336pFHHpEk1dXV6eGHH9bOnTtVUFCgwsJCbd++XV/5yleGP5HlSyFH\nEvQ8eTZznG7HgHAzod68GofN3CyF5Ff++siPwk23HQGPdPwK9lyeg9pMjV/5aygvC9f6i56witdT\ni9R+9gh3nwXdmXjRuXMi2SuIrj1bXr0KQDLC3UJuLtKmegyEW1AX+nNBU+E9wj0AJnZLgN9yqdnB\n9yHw3eGCaoCCDnk/HwPh4kXX7vb5hgY30uOCKgAgCeFuMTp3SMF37Zk+D7XlL8I9QMy9w3ajhXaq\neqam/cGcex54PX+Z6/PRwdstbHVGjQ3HnDtS4lUAUsm22/ZzHPAO4Z4HXqw5D/JE4KQLpyCnPHh1\naB6mZfLMtCVqbh4bZvPr2BPs/mFaBhnL9ATnZLPP0GPKD3W7Ee55ZtL0THzQB+zk9fGlazcX0zKG\n4GUwTGDCNCE+xbQMACAJ4R5RdEwYipqwC+FuiKCnZQA3vHjDE/XqL+bcDRNEyHNSIZVc6ofrPf5g\nzh0Z48RCOtSIHejcDeR1987JilxlUkd07f6hc7eMl+HOyQQ30tURwW4uwt1inEzwAnUUToS7oVh5\nANPRtZuNcAcACxHuBqMjhw1squFTp/5d8+ZdJUk6ePDftGLFDUm3f/nLxTp0aK3ef/8/1dl5t372\ns39VLJaPkRLuVuCHAPIhajV33XVTVFg4Vi0tJzV27BjNnz9df/rTWwO3l5R8Ts89t1qvv96tefMe\n1fr1z6qu7kb97Ge35mW8hHsIpFqpELUTDGZJV4M21ee//MvVeumlTjmOdNNNxTp9+n11dPxz4Pb1\n6+erp+eC7ryzQW+80a2Ghr/rpz99Xhs3flmf/WxB4ONlnXuIsHYdJhtak7bU6Lvv3iPHcTRuXIHG\njInpgw96NXbsZRo37jKdP39RjiNdfvl/6cCBNWprO6Pvf/+Zgftee+0UtbVt1M03/0ZHjnRk/dxu\n1rkH/+MErthywsA+tjYfc+fuVCwWU3PzOv3wh8/q1VdPac+eZXryyde0b98bA9tNmzZBL774VtJ9\nT506L0m66qqJgY5ZItxDxaYTBnazqVYTiX9qzpwrNXbsZXrmmb9rwoTPqLJymmpr96i7+/18D29U\nhDsAjOK119br6qsnqaBgjMaOvUxnz/6HxoyJady4Av3jH5skSeXlD6uz85xOnjw/rEMvKhovSTp5\n8lzgYyfcAWAU3/rWf+szn7lMv/lNrZqajuu3v/2btmy5RR9+2KcHHviTJOnkyY+nXg4dSmj16rlD\n7l+m9967qJaWU4GPndUyADCKjo5/qr29R3PnFunpp9/QiRM9mjOnSH/4wzGdONGjEyd6dOnSxxc7\nd+78X02aNE6//vUSzZp1hZYsuV73379Qv/zlUV248FHgY6dzB4AUqqqm6cMP+3Ts2Gl97nPjNHv2\nFXrhhf8btl1n5zl985tPaPv2b+rll7+vnp4LeuSRV3Tfff+Th1GzFBIAjMWv/AUAJCHcAcBChDsA\nWIhwBwALEe4AYCHCHQAsRLgDgIUIdwCwEOEOABYi3AHAQmnDff/+/frSl76kmTNn6uc///mI22za\ntEkzZ85URUWFWlpaPB8kACA7KcO9r69Pd911l/bv36/W1lbt3r1br7/+etI2jY2NOn78uNra2vTo\no49q/fr1vg7YXCfyPQAf2bxvEvsXdrbvX25ShvvRo0dVVlamGTNmaOzYsVqxYoX27duXtE1DQ4PW\nrFkjSaqurlZPT4+6urr8G7Gx2vM9AB+153sAPmvP9wB81p7vAfisPd8DMFLKcO/s7FRpaenA1yUl\nJers7Ey7TUdH9n8IFgDgnZS/z/3jX9Ob3tBfRzna/RxnS4bDCp943FE8buf+2bxvEvsXdrbvX86c\nFI4cOeIsWrRo4Ott27Y5DzzwQNI2dXV1zu7duwe+/uIXv+icOnVq2GNdd911jiQ++eSTTz6z+Kyo\nqEgV06NK2bnPnz9fbW1tam9v1/Tp07V3717t3r07aZva2lrt2LFDK1asUHNzsyZPnqyioqJhj3X8\n+PFUTwUA8FDKcC8oKNCOHTu0aNEi9fX1ad26dSovL9cjjzwiSaqrq1NNTY0aGxtVVlam8ePHa9eu\nXYEMHAAwusD+zB4AIDiev0PV5jc9pdu3AwcOaNKkSaqqqlJVVZW2bt2ah1HmZu3atSoqKtKcOXNG\n3Sasx01Kv39hPnaSlEgktHDhQs2ePVs33HCDfvGLX4y4XViPYSb7F9ZjeOHCBVVXV6uyslKzZs3S\nvffeO+J2WR+7nGbqR/HRRx851113nXPixAnn4sWLTkVFhdPa2pq0zbPPPussXrzYcRzHaW5udqqr\nq70cgm8y2bfnn3/eWbJkSZ5G6M4LL7zg/PnPf3ZuuOGGEW8P63Hrl27/wnzsHMdxTp486bS0tDiO\n4zjnzp1zrr/+emvOPcfJbP/CfAzfe+89x3Ecp7e316murnZefPHFpNtzOXaedu42v+kpk32TlNNf\nKTfBggULNGXKlFFvD+tx65du/6TwHjtJmjZtmiorKyVJEyZMUHl5ud5+++2kbcJ8DDPZPym8x7Cw\nsFCSdPHiRfX19Wnq1KlJt+dy7DwNd5vf9JTJvsViMR0+fFgVFRWqqalRa2tr0MP0TViPW6ZsOnbt\n7e1qaWlRdXV10r/bcgxH278wH8NLly6psrJSRUVFWrhwoWbNmpV0ey7HLuVqmWx5/aYnk2Qyxnnz\n5imRSKiwsFBNTU1aunSpjh07FsDoghHG45YpW47d+fPntXz5cj300EOaMGHCsNvDfgxT7V+Yj+GY\nMWP06quv6uzZs1q0aJEOHDigr33ta0nbZHvsPO3ci4uLlUgkBr5OJBIqKSlJuU1HR4eKi4u9HIYv\nMtm3iRMnDry8Wrx4sXp7e3XmzJlAx+mXsB63TNlw7Hp7e7Vs2TKtWrVKS5cuHXZ72I9huv2z4RhO\nmjRJ3/72t/Xyyy8n/Xsux87TcB/8pqeLFy9q7969qq2tTdqmtrZWjz/+uCSlfNOTaTLZt66uroGf\nrkePHpXjOMPmzsIqrMctU2E/do7jaN26dZo1a5Y2b9484jZhPoaZ7F9Yj2F3d7d6enokSR988IGe\ne+45VVVVJW2Ty7HzdFrG5jc9ZbJvTz31lHbu3KmCggIVFhZqz549eR515lauXKmDBw+qu7tbpaWl\nqq+vV29vr6RwH7d+6fYvzMdOkg4dOqQnnnhCc+fOHQiGbdu26a233pIU/mOYyf6F9RiePHlSa9as\n0aVLl3Tp0iWtXr1at956q+vc5E1MAGAh/sweAFiIcAcACxHuAGAhwh0ALES4A4CFCHcAsBDhDgAW\nItwBwEL/Dw0oRuIv2wV8AAAAAElFTkSuQmCC\n",
       "text": [
        "<matplotlib.figure.Figure at 0x112097290>"
       ]
      }
     ],
     "prompt_number": 79
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [],
     "language": "python",
     "metadata": {},
     "outputs": []
    }
   ],
   "metadata": {}
  }
 ]
}