{ "metadata": { "name": "" }, "nbformat": 3, "nbformat_minor": 0, "worksheets": [ { "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "This is a IPython notebook tutorial for defining a ricker wavelet function, and plotting it. The iPython notebook is designed to be used seamlessly with Matplotlib. When you start your own notebook, type: ipython notebook --pylab inline" ] }, { "cell_type": "code", "collapsed": false, "input": [ "%pylab inline" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stdout", "text": [ "Populating the interactive namespace from numpy and matplotlib\n" ] } ], "prompt_number": 26 }, { "cell_type": "code", "collapsed": false, "input": [ "import numpy as np\n", "import matplotlib.pyplot as plt\n", "\n", "def ricker(f, length=0.512, dt=0.001):\n", " t = np.linspace(-length/2, (length-dt)/2, length/dt)\n", " y = (1.-2.*(np.pi**2)*(f**2)*(t**2))*np.exp(-(np.pi**2)*(f**2)*(t**2))\n", " return t, y" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 27 }, { "cell_type": "markdown", "metadata": {}, "source": [ "Here we define the function \"ricker\" and state that it needs 3 input paramaters; frequency, f, the length of time over which we want it to be defined, length, and the sample rate of the signal, dt. Calling the function returns two arrays, the time axis (t), and the value of the function, y.\n", "\n", "To create a 5 Hz ricker wavelet, assign the value of 5 to the variable f, and pass it into the function like so," ] }, { "cell_type": "code", "collapsed": false, "input": [ "f = 5\n", "t, y = ricker (f)" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 28 }, { "cell_type": "markdown", "metadata": {}, "source": [ "To plot the result:" ] }, { "cell_type": "code", "collapsed": false, "input": [ "plt.plot(t, y)\n", "savefig( '5Hz_ricker_basic.png', facecolor='white')" ], "language": "python", "metadata": {}, "outputs": [ { "metadata": {}, "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAX8AAADICAYAAAAA/OEeAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3X1YVdWeB/DvIV/wpVSC0BJE0BRQXgRELr5gijIyvjMZ\nGnqVZtKbr13L7Nb1ZmOTYzdyqNFsxicdcG5jZnazNCkBX64gimKAJaipZCVyVRRQgTV/rEAUkMPZ\nZ5999j7fz/OcR/Dss9dvC/7OOmuvtX4mIYQAERE5FCetAyAiIttj8icickBM/kREDojJn4jIATH5\nExE5ICZ/IiIHpDj5z5kzB+7u7hg4cGCzxyxfvhze3t4ICQnByZMnlTZJREQKKU7+s2fPxq5du5p9\nPjs7G/v27UNOTg6WLl2KpUuXKm2SiIgUUpz8hw0bhm7dujX7fFZWFuLi4uDi4oL4+HgUFhYqbZKI\niBRSfcw/Ozsbfn5+9d+7ubmhuLhY7WaJiOg+2qjdgBAC9+4gYTKZGh3Xp08fvikQEbWSj48PioqK\nWv061Xv+4eHhKCgoqP/+0qVL8Pb2bnRccXFx/RuFER8rVqzQPAZem/mPlSsFBg0SOH1a4NVXV+Ct\ntwR69xb4+WftY+PPj9fX8GFpp9kmyX/btm24fPkytmzZAl9fX7WbJFLkm2+A9euBzz8HevcGnJyA\n3/8emDFDPgS3QiQDUDzsEx8fj4yMDJSWlsLDwwOvvfYabt++DQB49tlnMXjwYAwdOhShoaFwcXFB\nSkqK4qCJ1FJdDSxYAKxbB/TocfdzK1YAwcHA9u3AlCnaxEdkLSYh7KMfYzKZYCehqCI9PR1RUVFa\nh6EKI13bxo3Apk1AejpQd2uq4fV9/TUwdy5w8iTwwAOahWlVRvr5NcXo12dp7mTyJ/qVEICfn+z1\nN5crhACGDAFeegmYPNmm4RE1ydLcye0diH6Vlga0bQuMGNH8MSaTHP9/+23bxUWkBiZ/ol/9538C\n8+ffGe5pzpQpQFER8N13tomLSA0c9iECUFYmZ/acPw889FDLxy9dCjg7A//6r+rHRnQ/HPYhUuCT\nT4AxY8xL/ACQkAD8z/8AtbXqxkWkFiZ/IgD/+7/A9OnmHx8YCHTqBOTkqBcTkZqY/Mnh/f3vwOHD\nQExM6143aRLw6afqxESkNiZ/cnhffQUMHw506NC61zH5k54x+ZPD27kTiI1t/etCQ4ErV4BTp6wf\nE5HamPzJodXUAF9+CYwb1/rXOjnJm8R79lg/LiK1MfmTQ8vJAdzdgV69LHt9dLRcHEakN0z+5NAs\n7fXXGTUK2LtXbghHpCeKk39mZiZ8fX3Rt29fJCcnN3q+srISs2bNQnBwMEaMGIEdO3YobZLIavbu\nBZ54wvLXd+8OeHgAR45YLyYiW1Cc/BctWoT3338faWlpeO+991BaWnrX85s2bUKnTp2Qm5uLzZs3\n4/nnn+dKXrILVVUyaUdGKjvP6NEc+iH9UZT8r169CgAYPnw4evXqhTFjxiArK+uuY7p06YLy8nLc\nvn0bZWVl6NixY5NlHIlsLSsL8PcHHnxQ2Xmio3nTl/RHUfI/fPgw+vfvX/+9n58fDh06dNcx8fHx\nqKmpgaurK4YOHYrU1FQlTRJZTXp681s3t8bQofLG8a1bys9FZCuqF3B/99130aZNG1y8eBEnTpxA\nbGwsfvjhBzg5NX7f+dOf/lT/dVRUlKELMJD2MjLkBm1KPfgg0LcvcPSo3OufSE3p6elIT09XfB5F\nu3pevXoVUVFRyM3NBQAsWLAAMTExiG2wYubJJ59EYmIixo4dC0DW9N20adNdnxgA7upJtnXzJuDq\nCpSUmL+Z2/3Mnw94eVnnzYSoNTTZ1bNLly4A5Iyfs2fPYs+ePQgPD7/rmFGjRuGvf/0ramtrcfr0\naZSVlTVK/ES2duQI8Pjj1kn8gBz6OXDAOucisgXFwz7vvPMOnn32Wdy+fRsLFy6Eq6sr3n//fQCy\ngPtTTz2FgoIChIaGws3NDWvXrlUcNJFShw4BERHWO19kJLBwoSzzyPkMpAcs5kIO6ckngQkTgKef\ntt45PT1lgfe+fa13TqKWsJgLUSscOmT9m7ORkcD+/dY9J5FamPzJ4ZSUABUVgI+Pdc8bESHXDhDp\nAZM/OZysLNnrt/bYfFgYK3uRfjD5k8NRY8gHkKUdCwrkNFIie8fkTw7n0CHgnhnJVtGxo5w+mpdn\n/XMTWRuTPzmUmhq5Ejc0VJ3zh4bKesBE9o7JnxzK99/L4i3duqlz/tBQjvuTPjD5k0M5ehQICVHv\n/GFh7PmTPjD5k0M5cgQYNEi98w8cCBQXAzduqNcGkTUw+ZNDOXpU3eTfrh0wYADw616HRHaLyZ8c\nRm2tTMpqJn+A4/6kD0z+5DCKi4GuXeVWzmrijB/SA9ULuAOy4ldYWBh8fX1ZoIU0o/aQT52QEA77\nkP1TvKtncHAw1q5di169emHs2LHYv38/XBt0rYQQCAgIQFJSEkaPHo3S0tK7nq8PhLt6ksqWLZNV\nt155Rd12bt0CunQBysqADh3UbYtIk109zSngnpOTg4CAAIwePRoAmkz8RLag9kyfOu3ayZW+336r\nfltEllK9gPvu3bthMpkwbNgwjB8/Hrt371bSJJFFhFB/jn9DQUHA8eO2aYvIEqoXcK+qqsKxY8eQ\nlpaGiooKREdH49tvv0WHJj4Ps4A7qeXcOaB9e7m61xYCA4Fjx2zTFjkWaxVwV5T8w8LC8MILL9R/\nn5+fj5iYmLuOiYiIwM2bN9G9e3cAQGhoKDIzM+sLujfUMPkTWVNenkzIthIUBGzfbrv2yHHc2zF+\n7bXXLDqP6gXchwwZgoyMDFRUVKCsrAy5ubmIjIxU0ixRq504IVff2kpgoBz2qa21XZtEraF6AfeH\nH34Ys2fPri/gvnLlSnTu3Flx4EStkZcHxMbarr2HH5Yzfs6eBby9bdcukblYwJ0cgr8/kJoqh2Ns\n5R//EZgzB5gyxXZtkuNhAXeiZty8CZw+Dfj62rZdzvghe8bkT4ZXWCiHXtq3t227deP+RPaIyZ8M\n78QJICDA9u0GBXG6J9kvJn8yvLw82870qePjA1y+DFy5Yvu2iVrC5E+Gp1XP38lJ7u3PoR+yR0z+\nZHha9fwB3vQl+8XkT4ZWWipLKnp6atN+YKB88yGyN0z+ZGh1K3tNJm3aDwhg8if7xORPhmbrbR3u\nNWAAkJ8P1NRoFwNRU5j8ydDy8rS52VvnoYfkTqJFRdrFQNQUJn8yNK17/gCHfsg+MfmTYdXWyiEX\nJn+ixmxSwB2QVb/atGmDTz75RGmTRGY5ffrO7ppaYvIne6Q4+S9atAjvv/8+0tLS8N5776G0tLTR\nMTU1NVi2bBliYmK4cyfZjFaLu+7F5E/2SPUC7gCQnJyMuLg4uLm5KWmOqFW0XNzVkI8PcOkS8Ot/\nFyK7oHoB95KSEuzYsQPz5s0DIPeeJrIFe+n5P/CArCfw7bdaR0J0h+oF3BcvXow333yzvuDA/YZ9\nWMCdrCkvD7CwvKnV1Q39sIIpKWWtAu6KKnldvXoVUVFRyM3NBQAsWLAAMTExiG1QL8/b27s+4ZeW\nlqJjx4744IMPMGHChLsDYSUvsqKKCnmz99o1oG1braMBkpOBggJg3TqtIyGj0aSSlzkF3E+fPo0z\nZ87gzJkziIuLw7p16xolfiJry88H+vWzj8QP8KYv2R/VC7gTacEeFnc1NHCgjKm2Vm71TKQ1FnAn\nQ1qyBHj0UeCFF7SO5A4PDyAzE+jdW+tIyEhYwJ2oAXuZ5tkQh37InjD5k+EIof2Gbk1h8id7wuRP\nhvPzz3JsvUcPrSO5G5M/2RMmfzKcusVd9raekMmf7AmTPxmOPY73A8DjjwPnz8s1CERaY/Inw7GX\nbR3u1batXHuQn691JERM/mRA9trzBzj0Q/aDyZ8MpboaOHlSbqRmj5j8yV4w+ZOhFBXJxV2dO2sd\nSdOY/MleMPmTodjzkA8gk//x43ItApGWmPzJUPLygMBAraNonru7vPFbUqJ1JOTomPzJUOxxZe+9\nAgM59EPaU72Ae2pqKgIDAxEYGIjp06fj+++/V9okUbP0kPw57k/2QPUC7t7e3sjMzMTx48cxduxY\nvP7660qbJGrS1atAaSng7a11JPfH5E/2QPUC7hEREfVFX2JjY5GRkaGkSaJmnTghp3ja+375TP5k\nD1Qv4N7Qhg0bMH78eCVNEjVLD0M+AODrCxQXAzdvah0JOTLVC7jXSUtLQ0pKCg4ePNjsMSzgTkrY\n67YO92rfHvDxAQoLgaAgraMhvdFNAXcAyMvLw5QpU7Br1y706dOn6UBYyYsUiowE3ngDGDFC60ha\nNn06EBMDzJypdSSkd3ZbwP3cuXOYOnUqUlNTm038RErV1tpf3d774bg/aU31Au4rV65EWVkZ5s6d\nCwBo27YtsrOzlTZLdJcffgC6dAFcXLSOxDwBAcA772gdBTkyFnAnQ9ixA9iwAdi5U+tIzHPhAhAa\nCvz0k9aRkN6xgDs5NL3M9Knz2GPArVuy5CSRFpj8yRD0lvxNJhnviRNaR0KOismfDEFvyR/gTV/S\nFpM/6V5FBXDunKyRqydM/qQlJn/Svfx8WRu3bVutI2kdJn/SEpM/6Z4eh3wAuQ/RyZOy9CSRrTH5\nk+7pNfl36gT07Alwl3PSApM/6V5uLhAcrHUUluHQD2mFyZ90rbYWOHZMv8mfVb1IK0z+pGvFxXJL\nB71s63Av9vxJK0z+pGtHj+q31w/I5H/8uNZRkCNi8iddy80FBg3SOgrL9eoly0+WlWkdCTka1Qu4\nA8Dy5cvh7e2NkJAQnDx5UmmTRPX03vN3cpLbUHObB7I11Qu4Z2dnY9++fcjJycHSpUuxdOlSpU0S\nAQCE0H/PH+C4P2lD9QLuWVlZiIuLg4uLC+Lj41FYWKikSaJ6Fy7InnOPHlpHogyTP2lBUTGX5gq4\nNyzjmJ2djYSEhPrv3dzcUFxcDB8fHyVNt8q1a8C778qvTaaWH9Y8ru6YNm0AZ+fGjw4d5J+dOwPd\nuulviwIt1fX66/6d9SogAPjwQ62j0BchgBs35L2SigqgshKoqpKPuq9v3gRqauSjtvbOnw2/bvh3\nDbfEb+7rlr6/97m4OPvdc0r1Au5CiEaFBkzN/G9Vq4C7EMD16/LPlh51x1v7uOpq+ctY9wt67+Pa\nNeDKFaBjxztTF93dAS8veVPQy0v+Evn6yjcL0v94f52BA+X+RDU1wAMPaB2N9oQALl6U90FOn5ab\n9tU9Sktlwv/732VHqVs32XG6tzPVoQPQrp3sdDk5yX/Xhn/e+3cmk/y6oYZp6t6Udb/vG36txtYd\nuingnpycjOrqaixZsgQA4OPjg+Li4saBsJIXamuB8nLg8mX5C37xoixP+MMPwNmzQGGhnNfu4yN7\nvMOHAyNHAt7e+u/9WmLiRCAhQfau9K5vX+DTT+V+P47mxg3g4EEgMxM4cEAOgdXdCO/TB/D0lA8P\nD+CRR2THqFs3oH17rSO3D5bmTkU9/4YF3D09PbFnzx6sWLHirmPCw8Px/PPPY+bMmdi9ezd8fX2V\nNGloTk6yDm2XLjKhN6WqCigoAHJygG++AV59Ve4RExcH/NM/GWMYxFxHjgBJSVpHYR2DBwPZ2Y6T\n/MvKZOnNbduAjAy50nnECODFF+WnOXd3rSM0PtULuA8ePBhDhw5FaGgoXFxckJKSojhoR+bsLBP8\noEHAv/yL/Ih87BiwdSswbZp8fsEC4Omn5ZuCUZWUyGG03r21jsQ6Bg8GDh8GZs/WOhJ15eQA770H\nbN8OPPEEEB8PpKbKDg/ZFgu4G4gQ8tNAcjKwfz+weLF8dO6sdWTWt3078F//pZ+C7S05eBBYuFAm\nRyM6cAB45RU5hj9vHpCYCLi5aR2VMbCAO8FkAkaNkmPHf/ubHB7q21cmydparaOzruxs2Vs2iuBg\n+fOqqtI6Eus6cwaYMAGYMQOYOVPes3rpJSZ+e8Dkb1B9+wJbtsie8YYN8sbwd99pHZX1ZGUZK/l3\n6AD07y+H8Iygpgb4938HwsKAiAj5uzd7tpx9Q/aByd/gBg2SnwKmTAGGDgU2bdI6IuVqauTwiJGS\nPyAT5eHDWkehXEkJMHo08MUX8hPa8uWcmWOPmPwdwAMPAIsWAXv3Am++KXtgN29qHZXlvvtOTvl7\n+GGtI7Guuhk/epaeDoSEyOHHr79uftYaaY/J34EMGCB7zOXlQHS0fneSNNqQTx29J/9Nm+SMs9RU\neXOXC9bsG5O/g+nUCfi//wPCw4Hf/EaumtSb7GwZv9H4+gI//ihXeuvNqlXAypVyzv6oUVpHQ+Zg\n8ndATk7AmjVyncDIkfp7AzDaTJ86bdrIWT966v0LAaxYIScXHDggb1qTPjD5O7DnnwfmzweiovTz\nBlBZKbe5CArSOhJ1DB0q12joxYoVcs3F3r1A9+5aR0OtweTv4JYsAZ57DoiJ0cc9gKwsuQtmhw5a\nR6KOYcOAffu0jsI8ycnARx/JG7uPPKJ1NNRaTP6E3/8e+Id/ACZNsv9FRvv3y96xUUVEyOmet25p\nHcn9bd0KrF4N7NrFBVt6xeRPAOQ9gEcfldNA7XmXjX37ZO/YqLp2lTtZHj2qdSTNO3wY+N3v5AJC\no+yt5IiY/AmAvAn84YdAURHw9ttaR9O06mrg0CFj9/wB+x76+eUXYOpUuWo8MFDraEgJRcm/vLwc\nEydOhKenJyZNmoTr1683Oub8+fMYOXIk/P39ERUVhS1btihpklTk7Cy32F2zRm4QZ2/y8oCePY23\nuOtew4bZ503f27fltuGzZgGTJ2sdDSmlKPmvW7cOnp6eOHXqFHr27In169c3OqZt27ZISkpCfn4+\nPv74Y7zyyisoLy9X0iypyNNTLtKZMcP+ZgDt22f8Xj9wZ8aPvW3G9+KLcofYBgX3SMcUJf/s7Gwk\nJiaiffv2mDNnTqPi7QDQvXt3BP06L8/V1RX+/v7IMeq+tQYxapTcDmLGDLmPjr0w+nh/nUcflWP/\nhYVaR3LHzp3AJ58AKSlcuWsUipJ/wwLu/fv3R3YLq1OKioqQn5+PwUZcoWMwL74oa6D+279pHYlU\nUyPnko8cqXUktmFP4/4XLwLPPCMTf7duWkdD1tLiBqvR0dH46aefGv39qlWrWlVAoLy8HNOmTUNS\nUhI6NVNiSq0C7tR6Tk5yr5aQELkPkNbbKRw7Jkv7PfaYtnHYysiRwOefA3PnahtHbS3w298C//zP\njvGpSw+sVcAdQoEpU6aIo0ePCiGEyMnJEVOnTm3yuFu3bono6GiRlJTU7LkUhkIq2bZNCB8fIa5d\n0zaON98UYsECbWOwpZISIVxchKiu1jaOd94RYsgQIW7f1jYOap6luVPRsE94eDg2btyIyspKbNy4\nEUOGDGnqzQWJiYkYMGAAFi9erKQ50sCUKcDw4cCyZdrGsWeP3CPeUTz6qHwcOaJdDMXFwOuvA5s3\nswiLESlK/vPmzcO5c+fQr18/lJSUYO6vn1F//PFHxMbGAgAOHDiAlJQUfPPNNwgODkZwcDB27dql\nPHKymbffBj77TO7VroXKSrmtg6ONAkZHA199pU3bQsihnpdeklXhyHhYwJ3M8vnncgZQXp7cFtqW\n0tLkBmIHDti2Xa3t2gW88QaQmWn7tjdskLWfDx5kr9/eWZo7mfzJbAkJcoHVO+/Ytt0lS+Qskz/+\n0bbtaq2iQm6Y9uOPwEMP2a7d8+dl+c+9e2UBILJvluZObu9AZlu7VhaCsWUPXAhgxw5g4kTbtWkv\nOnYEIiOB3btt16YQcobRggVM/EbH5E9mc3EB3n0XmDNHjsPbQn6+nG4YEGCb9uzN5Mlyv3xbSU2V\nPf+XXrJdm6QNDvtQq02bBnh5yS191bZqFfDzz8B//If6bdmjixcBPz/5b9Cunbpt/fyzfJPduRMI\nDVW3LbIeDvuQzSQnyx1ADx9Wv63PPnPMIZ86PXrI2r5796rf1vz5cktvJn7HwORPrfbII3L6Z2Ki\nukVHLlyQW0wPH65eG3pgi6GfTz6RM7lWrFC3HbIfTP5kkenTgV691N37JzUViIsD2rZVrw09mDJF\nJv/bt9U5f1mZ7PX/938btzwmNcbkTxYxmYD16+UN4BMnrH9+IeTK0oQE659bb3x85EKrL79U5/yL\nFsk3WUfYLpvuYPIniz32mFyENGeOrLJlTbm5ckZRZKR1z6tXv/2tvM9ibZ99Bvztb/azeyvZDpM/\nKfLMM3IBUlKSdc9b1+s3max7Xr168klZXe3SJeuds6wMmDcP2LjR9qu2SXuc6kmKnT4NDB4stwJ4\n/HHl56uokBXFsrMBb2/l5zOKhAS58nbJEuudr1s3x51GaxSc6kma8fYGXn1VDk1YY/gnJQX4zW+Y\n+O/1u9/JeyzWqK62YweHexydxcnfnOLtdWpqahAcHIzx48db2hzZuQULgAcfVF7ftbZW7h20aJFV\nwjKUiAhZ0ObTT5Wd58IF4NlnZbEeDvc4LouTvznF2+usXbsWfn5+MHEA17CcnOQ4/caNcmzaUlu3\nynsITzxhvdiM5MUX5apnS4u7V1cD8fHyzZU30x2bxcnfnOLtAHDhwgV88cUXeOaZZzimb3Du7rI3\nmZAAlJS0/vXV1fKTw2uv8UZvcyZOlAXUt2617PV//KOcy691cR7SnsXJ39zi7UuWLMGaNWvg5MTb\nC44gOhp47jlg0qTWb/6WnCynj44Zo05sRmAyAW++CSxfDty40brXpqYCW7bIeyr870j3LdOgtHj7\n559/jkceeQTBwcFmFRxmAXdjWL4cKCgAZswAPvrIvBW6Z87I4YwDB9jrb8moUfKG+Kuvym02zJGR\nIWcJffON3J6D9EvzAu7mFG9fvny56Nmzp/Dy8hLdu3cXHTt2FAkJCU2eT0EoZIeqqoSIjRUiLk6I\nW7fuf2xlpRChoUK8/bZtYjOCS5eE8PAQYvv2lo/9+mshXF3ln2Q8luZOiz/8mVO8/Y033sD58+dx\n5swZ/OUvf8ETTzyBzZs3K3irIr1o3x74+GOgqkoOBTXxARKAfH7qVLl9weLFto1Rz1xdgW3bZJ3d\nr79u/riUFLkF99atvIlOd7M4+ZtTvP1enO3jWJyd5bTE4cPl4qQPPrgzTi0EsH8/MGQI0KWLvFHM\nX4/WCQuTbwDx8cDKlcDVq3eey8uTb6qvvy6HejiCSvfiCl+yiexsOZNn/365eresDOjcGXj5ZWDW\nLCZ+Jc6dk1NAv/wS6NlTvgmYTHJR2JIl8k2YjIsF3EkXrl4Fzp4FunaVbwJM+tZz7Rrwww/yTbVX\nL87ocRRM/kREDoh7+xARkdmY/ImIHBCTPxGRA2LyJyJyQEz+REQOiMmfiMgBMfkTETkgJn8iIgfE\n5E9E5ICY/ImIHJDqBdxv3LiBWbNm4fHHH4efnx8OHTpkcbB6ZpXiC3bKyNcG8Pr0zujXZynVC7iv\nWLECnp6eyMvLQ15eHnx9fS0OVs+M/Ato5GsDeH16Z/Trs5TqBdzT0tLw8ssvw9nZGW3atEGXLl0s\nDpaIiKxD1QLuFy5cQFVVFebNm4fw8HCsXr0aVVVVlkdLRETWcb8aj6NHjxYDBgxo9NixY4fw8PAQ\nlZWVQgghbty4ITw9PRu9/tSpU8JkMonPPvtMVFRUiISEBLFp06Ym2/Lx8REA+OCDDz74aMXDx8fH\nohq+bXAfe/bsafa5TZs2obCwEMHBwSgsLERYWFijY/r06YN+/fph/PjxAID4+Hhs3rwZM2fObHRs\nUVHR/UIhIiIrUrWAOwD07dsXWVlZqK2txc6dOzF69GiLgyUiIutQvYD7W2+9hUWLFmHQoEFwdnbG\nU089pTxqIiJSxG7KOBIRke1otsLXnEViVVVVCA8PR1BQEIYMGYKkpCQNIrWMOdd3/vx5jBw5Ev7+\n/oiKisKWLVs0iLT1zF3gN2fOHLi7u2PgwIE2jtAymZmZ8PX1Rd++fZGcnNzkMcuXL4e3tzdCQkJw\n8uRJG0eoTEvXd/LkSURERMDZ2Rl//vOfNYhQmZauLzU1FYGBgQgMDMT06dPx/fffaxClZVq6th07\ndiAwMBBBQUGIjY3F4cOHWz6pRbeJrWD16tVi/vz5oqqqSjz33HNizZo1TR5348YNIYQQVVVVwt/f\nX5w6dcqWYVrMnOu7ePGiyM3NFUIIcenSJdG7d29x7do1W4faaub+7DIzM8XRo0fFgAEDbByhZYKC\ngkRGRoY4e/as6Nevn7h06dJdz2dlZYnIyEhx+fJlsWXLFhEbG6tRpJZp6fp++eUXcfjwYfGHP/xB\nvPXWWxpFabmWru/gwYPiypUrQgghPvzwQ/H0009rEaZFWrq269ev13+dnp4uhg0b1uI5Nev5m7tI\nrGPHjgCA69evo7q6Gu3bt7dlmBYz5/q6d++OoKAgAICrqyv8/f2Rk5Nj61Bbzdyf3bBhw9CtWzcb\nR2eZq1evAgCGDx+OXr16YcyYMY2uKysrC3FxcXBxcUF8fDwKCwu1CNUi5lyfm5sbQkND0bZtWy1C\nVMSc64uIiKhfZBobG4uMjAybx2kJc66tU6dOdx3v7Ozc4nk1S/7mLBIDgNraWgQGBsLd3R3z58+H\nh4eHLcO0mLnXV6eoqAj5+fkYPHiwLcJTpLXXpgcNrwlAk/tQZWdnw8/Pr/57Nzc3FBcX2yxGJcy5\nPj1r7fVt2LChfgq6vTP32rZv3w4vLy/MmTMHH3zwQYvnve88f6Wio6Px008/Nfr7VatWQZh5n9nJ\nyQnHjx/H2bNnMW7cOERGRiI4ONjaoVrEGtcHyDH0adOmISkp6a53cC1Z69qMRAjR6NpNJpNG0ZCl\n0tLSkJKSgoMHD2odilVNnjwZkydPxkcffYRJkyYhNzf3vsermvyVLhJryMvLC+PGjUNWVpbdJH9r\nXN/t27cxdepUJCQkYOLEiWqF2mrW/NnpQVhYGF544YX67/Pz8xETE3PXMeHh4SgoKMDYsWMBAJcu\nXYK3t7dN47SUOdenZ+ZeX15eHubOnYtdu3aha9eutgzRYq392U2bNg0LFy5EZWUlOnTo0Oxxmg37\nmLNIrLRxDpx7AAABmUlEQVS0FFeuXAEAXL58GV999ZVdJcj7Mef6hBBITEzEgAEDsHjxYg2itIy5\nC/z0pG4sODMzE2fPnsWePXsQHh5+1zHh4eHYtm0bLl++jC1btuhqh1pzrq+OHj/ZmXN9586dw9Sp\nU5Gamoo+ffpoEaZFzLm24uLi+p/bF198gZCQkPsmfgDazfa5du2amDBhgvDw8BATJ04U5eXlQggh\nSkpKxLhx44QQQhw/flwEBweLgIAAMWbMmGb3BbJH5lzfvn37hMlkEoGBgSIoKEgEBQWJL7/8Usuw\nzWLOtQkhxFNPPSV69Ogh2rVrJ3r27Ck2btyoVchmSU9PF/379xc+Pj5i7dq1Qggh1q9fL9avX19/\nzLJly4SXl5cYNGiQKCgo0CpUi7R0fRcvXhQ9e/YUDz30kOjatavw8PCo/9nqQUvXl5iYKFxcXOr/\nr4WFhWkZbqu0dG2rV68W/v7+IigoSMyePVucOHGixXNykRcRkQNiGUciIgfE5E9E5ICY/ImIHBCT\nPxGRA2LyJyJyQEz+REQOiMmfiMgB/T+dpHdm/Gp2vAAAAABJRU5ErkJggg==\n", "text": [ "" ] } ], "prompt_number": 47 }, { "cell_type": "markdown", "metadata": {}, "source": [ "But with a few more commands, we can improve the cosmetics," ] }, { "cell_type": "code", "collapsed": false, "input": [ "plt.figure(figsize=(6,4))\n", "plt.plot( t, y, lw=2, color='black', alpha=0.5)\n", "plt.fill_between(t, y, 0, y > 0.0, interpolate=False, hold=True, color='blue', alpha = 0.5)\n", "plt.fill_between(t, y, 0, y < 0.0, interpolate=False, hold=True, color='red', alpha = 0.5)\n", "\n", "# Axes configuration and settings (optional)\n", "plt.title('%d Hz Ricker wavelet' %f, fontsize = 16 )\n", "plt.xlabel( 'two-way time (s)', fontsize = 14)\n", "plt.ylabel('amplitude', fontsize = 14)\n", "plt.ylim((-1.1,1.1))\n", "plt.xlim((min(t),max(t)))\n", "plt.grid()\n", "savefig( '5Hz_ricker_pretty.png', facecolor='white')" ], "language": "python", "metadata": {}, "outputs": [ { "metadata": {}, "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAYoAAAEfCAYAAABf1YHgAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3Xd4FOXaP/DvbHrvvRJaCiW0hJ4EpUgRlBop0quAIJ4D\nL6LoeS2o/MTXoxQL4AE8iB4VUcFQQpCYBqGlQEjvjfSElN3n98eSPcQUNsnO7szm/lzXXjC7szP3\n3jzsvfM8M89wjDEGQgghpB0STQdACCFE2KhQEEII6RAVCkIIIR2iQkEIIaRDVCgIIYR0iAoFIYSQ\nDlGhIF0SHh4OiUTS6mFtbf3E92ZkZEAikeDLL79s8/WxY8ciJCSk2zHu3r27RWxWVlYYPnw4jhw5\n0u66ndGcg4sXL3Y7Vm115MgRSCQSZGVldep9FRUV2L17N+Lj43mKjHSGrqYDIOL2ySefYMSIEYpl\nXV3lmxTHce0+395rXXH16lXo6OggISEBR48exfLly9G/f3+MGjVKsc6qVaswdepUle2TdE9ZWRne\neustuLu7Y8iQIZoOp8ejQkG6xcfHBwEBASrdpqqvAQ0MDIREIkFAQAAWLFgAe3t7HDt2rEWhcHFx\ngYuLi0r3210NDQ3Q19fXdBgaRdcDCwN1PZFuUcd/ZE9Pzza7uTrbVQQARkZGGDduHP7zn/+0eL6t\nrifGGPbv34+AgABYWFjA3t4ezzzzDO7evdvu9tPS0tC3b1+MGzcO5eXlAID09HTMmzcPnp6esLGx\nwcSJE3Hr1q0W71u6dCnc3NyQnJyMqVOnws7ODn/729/a3MfevXthYmKCpqYmxXOzZ8+GRCLBhQsX\nFM99/vnn0NPTQ3V1NQAgNjYWc+bMgZubGywtLbFw4UL88ssvLbazYcMGODo6QiqVtthnfX09rKys\nsGXLFsVzNTU12Lp1K3x9fWFubo6AgAB8++237ebmcadOncKYMWNgZWWF3r17Y9OmTairqwMg75r0\n8vICID/Sa/63/vrrr5XaNlE9KhSkWxYuXAh9fX306tULK1euRHFxsdLvlUqlaGpqavX4q59++glR\nUVGKR1hYGOzs7ODj49PpeKVSKeLi4hRfRI/7a3fXihUrsGHDBhgbG+Ozzz7D559/Dj8/PxQUFLS5\n7fj4eIwePRp+fn44f/48LC0tce/ePQwaNAi3b9/G+vXrsXfvXpSWlmLkyJGtCk5FRQUmT54MZ2dn\nfPPNN1i4cGGb+5kwYQLq6urw559/ApAXtPDwcBgbG7cYL7l48SKGDx8OU1NTAEBWVhYGDx6Mzz77\nDN9++y369OmDJUuWYMOGDYr3LFmyBEVFRfj9999b7PPMmTOoqKjAkiVLAMiPdkaOHIljx45hypQp\nOHLkCBwdHREaGtrmGNDj3nnnHbzwwguwsbHBwYMHMXfuXBw7dkzR9efs7Kwo5P/zP/+j+HenrkEN\nYoR0QXx8PHv11VfZmTNn2G+//caWLVvGLCwsmKenJ6uqqurwvenp6YzjuA4fISEhbb5XKpWyGTNm\nMHt7e5aent7hft544w3GcRyrr69njY2N7NatW+yll15iHMexffv2tblus8jISMZxHNu+fXu72790\n6RLjOI5duHCBnT9/npmZmbGVK1cymUymWGf69OnMwsKC1dTUKJ6rr69n5ubmbP78+YrnXnzxRcZx\nHDtw4ECHn6k5B9bW1uzNN99kjMn/LSQSCXvllVfYqFGjFOs5OjqyHTt2tLkNmUzGGhsb2Z49e5iO\njg6rrq5WvNavXz8WGhraYv2ZM2cyPz8/xfKHH37IOI5jV65cabHe5MmTmbOzsyIHhw8fZhzHsczM\nTMYYY9nZ2czQ0JCtWLGixfu++eYbxnEc++233xhj/20jX3755RPzQfhHhYKozLfffss4jmNvv/12\nh+s1fwm8/vrr7Nq1ay0ecXFxzN/fv91CsW3bNmZoaMgiIyOfGE/zl/9fH5999lm76zZ75513GMdx\nLC4urt3tNxeKNWvWMAMDA7Zz584Wr8tkMmZsbMzWrVvHGhsbWzzmzZvHPDw8FOs2F4rmL9Qnee65\n51hQUBBjjLG9e/eyIUOGsD///JPp6emx6upqlpCQwDiOY7///rviPdXV1ezQoUNszpw5zNraWpEP\niUTCbt68qVjvH//4BzM2NlYU/JKSEqavr8/ee+89xTrPPPMM8/X1ZU1NTS0+16effso4jmMZGRmM\nsdaF4uTJk4zjOBYZGdnifaWlpczAwIC98cYbjDEqFEJDXU9EZWbNmgVbW1tcv35dqfU9PDwwdOjQ\nFo9hw4bBxMSkzfW//PJL/L//9/9w+PDhFgPRTxIdHY2YmBgcOnQIfn5+2LNnD/Ly8jp8T3Z2NoyM\njDBs2LAnbv/777+HsbExXnzxxRbPl5SUoK6uDgcOHIC+vn6Lx6lTp5Cbm9tifScnJ7i7uyv1mYKD\ngxEVFYWHDx/i0qVLCAkJwYgRI2BoaIiIiAhcunQJenp6GDt2rOI9n332GTZu3IiRI0fi1KlTiI2N\nxb59+8AYU4ynAMCiRYvw8OFDfPfddwCAkydPoqmpCYsWLVKsk5mZiaSkJOjp6bX4XC+99BI4jkNm\nZmabcTc/P2bMmBbvs7W1RWNjY6dPoyXqQWc9EZVh8iNUlZ7a2uzy5ctYv3493njjDSxYsKBT7x02\nbBgkEgmGDx+OCRMmYOjQoVi7di1Onz7d7nvc3NxQV1eHa9euPbFYfP755/jggw8QHByMS5cuoV+/\nfgAAGxsbGBkZ4fnnn28xCNyezgzOh4SEoKGhAREREbhy5QrWrl0LHR0djBs3DhcvXkRaWhoCAwNh\nZGQEAJDJZNi/fz8WL16MV155RbGdkydPttq2p6cnxowZg2PHjmHp0qU4duwYQkJCWpwV5ubmBplM\nhuPHj7cZX3MO/srV1RUA8M0336Bv376tXrexsVE6B0R9qFAQlfnhhx9QWlqq1K/wzrh//z6ef/55\nzJs3D6+//nq3tuXl5YV//OMf2Lx5M8LCwjBx4sQ215swYQIA+dHCkz6Pubk5zp07h2eeeQbBwcG4\nePEivL29IZFIMH78eISFheGrr76Cnp5eh9vpTIEdMGAA7Ozs8MEHH6CmpgZBQUGKuI8dO4acnBys\nX79esX5DQwMePHgACwsLxXMymQw///xzm9tfsmQJ1q5di/DwcERFReHw4cMtXg8JCcGOHTugo6OD\nwYMHKx33mDFjYGBggKioKMyfP7/d9czMzAAAhYWFSm+b8IcKBemSRYsWoU+fPvD394dEIsF3332H\nH3/8EX379sXmzZu7vX322Gm306dPh4mJCVatWoWoqKgW640cObLT216zZg0++OAD7Nq1q91CERgY\niGXLlmHPnj2Ijo7GihUrYGZmhoiICEyfPl3xxdzM1NQUZ8+exbRp0xASEoKLFy/Cx8cHn3zyCYYN\nG4ZRo0bh2WefxZgxY1BeXo7r169DJpPh3XffbfMzPwnHcQgODsapU6cQEBCgOLMpJCQEr776KoD/\nFjsAMDQ0xNy5c3H48GE4OjrCxsYGH3zwQbvXacydOxcbN27EokWLYGxsjDlz5rR4fevWrfj3v/+N\nZ599FrNmzVLkMTk5GX/88Qd+/PHHNrfr7u6ON998E2+88QaysrIwdepUeHp6Ijs7G+fPn8fKlSsR\nHBwMGxsbODo64tSpUxgwYACcnJzg5eWl1JX/hAcaHSEhovXuu++yQYMGMQsLC2ZsbMz8/PzY1q1b\nWVlZ2RPf+6SByrFjx7YYzG4ecP3roLREIulwP7t372YSiYRJpdJWrx06dIhJJBL2888/M8bkg9lt\nbe+f//wnGz58ODM1NWV2dnZs2rRp7N69e4wx+WC2RCJhFy5cUKxfU1PDJkyYwBwdHVlCQgJjjLGc\nnBz24osvsr59+zJDQ0Pm6urKZs6cqTjDhzHGli5dytzc3Dr8PH+1f/9+JpFIWpzZJJPJmLW1NTMy\nMmL19fUt1s/NzWXbt29n7u7uzNfXl+3YsYNdvHiRSSQSdvny5Vbbnzt3LpNIJGzhwoVt7r+uro7t\n2rWLDR48mJmYmDB7e3s2fvx49vHHHyvWOXz4MJNIJK0G6X/55RcWEhLC7O3tmampKfPx8WEbN25k\nubm5inXOnTvH/P39mYmJCZNIJOzo0aOdyg9RHY4xuvSREEJI++isJ0IIIR2iQkEIIaRDVCgIIYR0\niAoFIYSQDmnl6bHBwcG4fPmypsMghBBRCQoKQnh4eKvntfKI4vLly4qrhIXweOONNzQegxAelAfK\nAeVB2Dlo7we2VhYKocnIyNB0CIJAeaAcNKM8iCsHVCgIIYR0iAqFGixdulTTIQgC5YFy0IzyIK4c\nCOrK7OXLl+OXX36Bvb09bt++3eY6O3bswMmTJ2FlZYXjx4/D29u71Tocx0FAH4sQQkShve9OQR1R\nLFu2DGfPnm339ZiYGFy5cgVxcXHYtm0btm3bpsbouq6tswh6IsoD5aAZ5UFcORBUoRg3bhysrKza\nfT06Ohpz5syBtbU1QkNDkZSUpMboCCGkZxJUoXiSmJgY+Pr6Kpbt7OyQmpqqwYiUExwcrOkQBKEn\n54ExhpKSEvTu3Rt1dXWaDkfjenJbaCamHIjqgrvmc30fx8fd1AhRpdzcXHz11Wlcv16IujrAwUEH\nM2cOxfTpE9u9HwQhQiKqQhEYGIjExERMnjwZAFBcXAwvL6821126dCk8PT0BAJaWlvD391dU8Oa+\nQXUt79u3T6P7F8py83NCiUcdy0lJydi06QOkpUlha+uH+voC3L9fitjYVKSk5OLll1/En3/+KZh4\n1bV848YNvPzyy4KJRxPLzc9pMp7w8HAcOXIEABTfl20R1FlPgPwilBkzZrR51lNMTAy2bt2Kn376\nCefOncOJEydw5syZVusJ7ayn8PBwxT9ST9bT8lBQUIBXX/0CN282wdt7KLy9pyIr6w/Y2vogOvrf\nqK8vw4oV/bFq1fxO3S9bG/S0ttAWIeagve9OQRWK0NBQXL58GSUlJXBwcMCbb76JxsZGAPLbVwLA\n9u3bcfLkSVhbW+PYsWPw8fFptR2hFQrS8zQ0NODNNw/gzJkH6NdvCHx9n23RTVpX9wCXL38OQ8M6\nvPXWJIwfP1qD0RIiJ4pCoSpUKIimnTt3Cbt3X4aRkQPGjFkJHR29VuuUlNzDlSsnMHSoPvbt2wAL\nCwsNRErIf4niOgpt9XifZE/WU/JQXl6O48evorYWGDx4WosikZERrvi7rW0/9Orli5s3G/DTT+c1\nEKnm9JS20BEx5YAKBSEq9vvvV5CQ0AQvrwGwsHDvcF1v78kAdPDzz3dQUlKingAJ6STqeiJEhaqr\nq7F27Ue4c0eGSZM2wNjY9onvSUz8GffvX8PWrYOxaNFzaoiSkLZR1xMhahAREYW7d6Xw8PBWqkgA\nQK9eY8GYBL/9dhtVVVU8R0hI51GhUAMx9UXySdvz0NTUhN9+u47GRsDLq+2zmB4fo2hmZGQFd3dv\npKTIEBNznecohUHb24IyxJQDKhSEqEhy8l3cvl0La2tHmJu7duq9np7D0dAA/PrrNchkMp4iJKRr\naIyCEBX5+OOvcfhwGoYPnwpX14BOvZcxhkuX/glDw1L8858voF+/fjxFSUj7aIyCEB5VVFQgMjIN\nEokuHBwGdvr9HMfBy8sfRUVAZOQtHiIkpOuoUKiBmPoi+aTNebh9OxEZGYCbWz/o6Rm1u15bYxTN\nHB3lBSYi4i4aGhpUHKGwaHNbUJaYckCFghAVuHIlEfX1gIuL75NXboehoSXs7d2RmtqIpKRkFUZH\nSPdQoVADoU38pSnamoeKigpcu5YNXV092Nh0PLbg6Rnc4evu7gNQUQH8+WeCCiMUHm1tC50hphxQ\noSCkmxISkpCVBbi49IWOTvfuL2FnJ5/kMiYmTTEhJiGaRoVCDcTUF8knbc1DREQCGhqU63bqaIwC\nAAwMzGBr64K0tEakpaWpKELh0da20BliygEVCkK6obq6GjduZEMi0X1it5OyXF37o7wcuHaNximI\nMFChUAMx9UXySRvzkJJyH9nZgJNTL6W6nZ40RgEAtrb9AQCRkfe09uI7bWwLnSWmHFChIKQbrl1L\nQU0N4OTUV2XbNDGxh6mpJVJTa5Cbm6uy7RLSVVQo1EBMfZF80rY8yGQyxMWlAgBsbPoo9Z4njVEA\n8ovv3Ny8UVQE3LlztzshCpa2tYWuEFMOqFAQ0kU5OTlIT38Ic3MbGBlZq3TbdnZ9IZMBsbGpKt0u\nIV1BhUINxNQXySdty0NiYgqKiwEXF+WOJgDlxigAwMLCHRyni4SEfNTW1nYxQuHStrbQFWLKARUK\nQrooJuY+GJP/+lc1HR09ODi4Iy8PWn2aLBEHKhRqIKa+SD5pUx5qamqQlJQPiUQXFhYeSr9PmTGK\nZo6OXqiqAm7f1r5CoU1toavElAMqFIR0QWZmJvLyAFtbN+jo6PGyD2vr3gCAuLhUmjafaBQVCjUQ\nU18kn7QpD0lJGaiokF8/0RnKjlEAgKmpIwwMjJGaWoEHDx50MkJh06a20FViygEVCkK64Pr1dACA\nlZUnb/vgOA5OTl4oKADu36ezn4jmUKFQAzH1RfJJW/JQU1ODe/eKoaurBzMz5069tzNjFABgb98L\nDx8Ct29ndOp9QqctbaE7xJQDKhSEdFJGRgby8wE7OzdIJLq87svSUj5QfutWFo1TEI2hQqEGYuqL\n5JO25CExMQPV1YCDg2en39uZMQoAMDKygaGhCTIzq7VqnEJb2kJ3iCkHVCgI6aQbNzIAAFZWnRvI\n7gqO4+Dg4IHCQiAjI5P3/RHSFioUaiCmvkg+aUMeqquruzw+AXR+jAIAbG09UF8PJCRoT6HQhrbQ\nXWLKARUKQjohIyMDhYWAnZ07JBIdtezTysoDjAE3b2pPoSDiQoVCDcTUF8knbchDcnIGKiu7Nj4B\ndH6MApBPO66vb4j09HJUVFR0ab9Cow1tobvElAMqFIR0wp07OeA4wNLSXW375DgJ7O3dUVgovyKc\nEHWjQqEGYuqL5JPY89DQ0ICUlEIAki6NTwBdG6MA5F1dNTVAcnJWl94vNGJvC6ogphxQoSBESXl5\necjPZ7C2duBtfqf2WFp6gOOAW7foiIKoHxUKNRBTXySfxJ6HzMxclJcDdnYuXd5GV8YoAMDMzBkS\niR7u3y9GTU1Nl/cvFGJvC6ogphxQoSBESbdv54AxwNLSVe37lkh0YGPjgsJC+Z31CFEnKhRqIKa+\nSD6JPQ+JifJCYW7e9ULR1TEKALC3d0V5OZCeLv5CIfa2oApiygEVCkKUUFlZiczMKhgZGcLIyEYj\nMVhaugEAbt/O1sj+Sc9FhUINxNQXyScx5yE7OwfFxYCtrQs4juvydro6RgHIj2QYA+7ezYNMJuvy\ndoRAzG1BVcSUAyoUhCghKSkHjY2AtbX6xyea6eubwNTUCnl5DSgqKtJYHKTnoUKhBmLqi+STmPOQ\nkCAfF+jO+ATQvTEKALCzc0VJCZCVJe7uJzG3BVURUw6oUBDyBFKpFMnJ+eA4wNy866fGqoKtrRsa\nG4HkZPEPaBPxEFShiIiIgI+PD/r27YtPPvmk1evh4eGwsLDAkCFDMGTIEPzv//6vBqLsPDH1RfJJ\nrHkoKipCfn4jLCysoadn3K1tdWeMAvjvEU1iorgLhVjbgiqJKQf83p6rkzZv3oyDBw/Cw8MDkydP\nRmhoKGxtbVusExQUhNOnT2soQtITpafnoKIC6NNHc+MTzUxMHCCR6CI9vRS1tbUwNu5e4SJEGYI5\nomieFXP8+PHw8PDApEmTEB0d3Wo9Md4OUkx9kXwSax7u3MkFoJoL7bo7RiG/8M4ZRUXivvBOrG1B\nlcSUA8EUitjYWHh7eyuWfX19ERUV1WIdjuMQGRkJf39/bN26FampqeoOk/RAzRfamZlpdnyimb29\nG8rLgYwM8RYKIi6C6np6kqFDhyI7Oxt6eno4evQoNm/ejDNnzrS57tKlS+Hp6QkAsLS0hL+/v6JP\nsLmSq2u5+TlN7Z+Wu75cV1eH6Og41NfrwNTUEcB/jwqaxxs6s+zpGdyt9wNAdXUeqqoycOeOFyZN\nEla+OrPcTCjx9MTl8PBwHDlyBAAU35dt4ZhA+nIqKioQHByM+Ph4AMDGjRsxZcoUTJs2rc31GWNw\ndHREVlYWDAwMWrzGcZwou6iI8KSk3MfixcdgZeWGwMAVmg4HAFBfX4XfftuLwEAD7N//d0gkgukY\nICLX3nenYFqYhYUFAPmZTxkZGQgLC0NgYGCLdQoLCxUf4ueff8agQYNaFQkh+usvqJ5KjHm4dy8X\nDx8C1taq6Xbq7hgFABgYmMHY2AL5+fUoLi7uflAaIMa2oGpiyoGgup727duHNWvWoLGxEZs2bYKt\nrS0OHjwIAFizZg2+++477N+/H7q6uhg0aBD27t2r4YiJtmu+o52FhebPeHqcvb0bCgsrkJOTAwcH\nB02HQ7ScYLqeVIm6nogqMMawfPkHiI+vxZQpL8PQ0FLTISlkZ0fh2rWz2L59CBYsmKnpcIiWEHzX\nEyFCU1ZWhpycWpiamsLAwELT4bRgbu4CjgMSE3M1HQrpAahQqIGY+iL5JLY8ZGXl4MED+R3tujNj\n7ONUMUYBAGZmTuA4HaSnF6O+vl4l21QnsbUFPogpB1QoCGlHQkIOZDLN3NHuSSQSXVhYOKCoiCEv\nL0/T4RAtR4VCDZrPX+7pxJYHVdzR7q+6O9fT4+ztXVFaKr9XhtiIrS3wQUw5oEJBSBsaGxuRklIA\nHR0OZmbOmg6nTVZWLmhqkt8rgxA+UaFQAzH1RfJJTHkoKChAUZEMlpZ20NVV3bU6qhqjAORHOhwH\n3L2bK7qz/MTUFvgiphxQoSCkDc0zxtraCm98opmRkTV0dY2Qm1utmFSTED5QoVADMfVF8klMebhz\nR96do+qBbFWOUXAcBxsbFxQXA7m54jpNVkxtgS9iygEVCkLaIO/OEc6Mse2xs3NFeTmQmUnjFIQ/\nVCjUQEx9kXwSSx6qq6uRkVEOfX19mJjYqXTbqhyjAABLSxcwJr4L78TSFvgkphxQoSDkL3JyclBS\nAtjauoDjhP1fxNxcXihSUvIglUo1HQ7RUsL+X6AlxNQXySex5OH+/VzU1AA2NqofyFblGAUA6OkZ\nw8TEBkVFTSgsLFTptvkklrbAJzHlgAoFIX+RkNA8Y6ywxyea2dq6oKREfAPaRDw6VShqamrw448/\norq6GoC8L7exsZGXwLSJmPoi+SSGPMhkMsVAtiqvyG6m6jEKQD6gXVUFpKWJZ0BbDG2Bb2LKgVKF\noqKiAqGhoXBwcMDs2bNRVFQEAHjllVewbds2XgMkRJ2Ki4uRn98AMzNL6OubajocpdBMsoRvShWK\nY8eOoaCgAFFRUTA2NlY8P3fuXJw7d4634LSFmPoi+SSGPOTk5KK0lL8L7VQ9RgHg0b28dZGZWYK6\nujqVb58PYmgLfBNTDpQqFF9++SV27NiBAQMGtHje09MT2dnZvARGiCbcvZuDhgbV3fpUHSQSHVhY\nOKKkBDSTLOGFUoUiPz8fJiYmrZ6PioqCra2tyoPSNmLqi+STGPKQmCgfyOZjfALgZ4wCkM8kW1Ii\nnplkxdAW+CamHChVKEaPHo2jR4+2mHisoqICX3zxBcaPH89bcISoU319PdLSisFxOjAzc9J0OJ1i\nY+OKhw+BlBRxFAoiLkrdMzszMxOBgYGwtrbG/fv3MXLkSMTExMDGxgYxMTFwcRHWYTrdM5t0RXp6\nOlatOgrGnDFu3GpNh9MpdXVl+P33jzFhgjE+/PBVld2Rj/Qs3bpntoeHB+7evYuFCxdi1qxZsLOz\nw86dO5GYmCi4IkFIV2Vm5qCsTH66qdgYGlpCV9cEBQW1KCsr03Q4RMsofR2FhYUFdu7ciW+//Rbf\nf/89du3aBQsLYd1wXqjE1BfJJ6HnITFRfutTKyv+CgVfYxQcx8HaWjwzyQq9LaiDmHKg294LERER\nSm+EximI2DHGeL3QTh3s7V2RmHgPWVk5GDhwoKbDIVqk3TEKiUS5gw2O4wQ3GRmNUZDOKi8vx4oV\n+5CTY4wpU8TZx//gQSquXPkXFi92xdatKzUdDhGh9r472z2iaL76GgCSk5OxfPlyLFu2DE899RQA\n4MKFC/j6669x6NAhHsIlRL1ycnJQXAzY2LiIskgA8iu0AQ4pKfloamqCrm67/70J6ZR2DxtsbW0V\nj127dmH9+vXYsWMHAgICEBAQgB07dmDdunXYuXOnOuMVJTH1RfJJyHlIT89BVRX/A9l8jVEAgK6u\nIYyNbVFSIhX8TLJCbgvqIqYcKNW/dOPGjTb7PP38/HDr1i2VB0WIuiUkyAeALSzEOT7RzMbGBUVF\n8iMkQlRFqUIxYMAA7N69u8WvlIKCArz11ls0aKYEMc3pwieh5kEqlSIlJf/RQDa/p3vzMdfT4+zt\nXVFZCWRkCLtQCLUtqJOYcqBUJ+bBgwcxbdo0ODk5wc/PD4wxJCYmwt3dHb/88gvfMRLCq4KCAhQU\nNMHMzBa6uoaaDqdbmu+hQTPJElVS6ojC19cX6enpuHHjBkJDQ/HCCy8gPj4eGRkZ8PPz4ztG0RNT\nXySfhJqHrKxslJQA9vZuvO+LzzEKADA1dQDH6SEr6wFqa2t53Vd3CLUtqJOYctCp0yIGDRqEQYMG\n8RULIRqRnJyN+np+bn2qbhwngYWFE4qLs5Cbm4u+fftqOiSiBZQqFP/5z386fP35559XSTDaSkx9\nkXwSah6Skppvfcr/EQXfYxSAfJwiMzML2dk5gi0UQm0L6iSmHChVKObMmdPh6zKZTCXBEKJulZWV\nyMqqgI6OAYyN7TQdjkrY2Lji7l0gJSUXEyZoOhqiDZQao5DJZC0elZWV+OabbzBq1CjEx8fzHaPo\niakvkk9CzENOTg6KiuRfruq40I7vMQrgv7dGTUrKEewMBUJsC+omphx06dJNU1NTzJ8/H42Njdiw\nYQP++OMPVcdFiFqkp2ejvBzw8+O/20ldDAzMoadnhoKCKpSWltLNxUi3KT17bFtcXFyQkpKiqli0\nlpj6IvnDHZ6OAAAgAElEQVQkxDzcuZMNxvidMfZx6hij4DhOceGdUGeSFWJbUDcx5UCpI4rr16+3\nWJbJZIiPj8f+/fuxYMECXgIjhG9NTU1ISckHwIl2xtj22Nu74vbtZGRk5GDw4MGaDoeInFJHFMOH\nD2/xCAgIwOuvvw4/Pz9s3bqV7xhFT0x9kXwSWh7y8/NRUCCFqamd2i60U8cYBSA/QmJMuBfeCa0t\naIKYcqDUEUVaWlqLZR0dHbi6qmfwjxC+ZGVlo7gYcHXVrqMJADAzc4Z8JtkCNDY2Qk9PT9MhERFT\n6ogiKysLLi4u8PT0hKenJ9zc3MBxHJqamjp1g6OeSkx9kXwSWh7u3s15dKGd+gay1TFGAQA6Ovow\nM7NHUZEMBQUFatlnZwitLWiCmHKgVKEIDg5u8z685eXlovqwhDRjjCEhIRsAYGmpPWc8Pc7eXj6g\nnZ0t7AkCifB166ynjIwMeHh4qCoWrSWmvkg+CSkPFRUVyMqqgp6eEYyMbNS2X3WNUQDya0MaGoDE\nxGy17VNZQmoLmiKmHHQ4RjFjxgzF3xcvXgx9fX0A8tPvamtrcePGjSdetd0ZERERWLNmDZqamrBp\n0yZs3Lix1To7duzAyZMnYWVlhePHj8Pb21tl+yc9R3Z2NoqK5Dcq0taxNgsLdwBAQkIWGGNa+zkJ\n/zosFDY2//2lZWVlBUPD/54Z4unpiXXr1mHixIkqC2bz5s04ePAgPDw8MHnyZISGhra4WCgmJgZX\nrlxBXFwczp07h23btuHMmTMq2z9fqHtOTkh5SEnJQkUF0KuXu1r3q64xCgAwMrKBgYEJsrKq8eDB\ngxb/nzVNSG1BU8SUgw4LxZEjRwDIi8Krr74KExMT3gKpqKgAAIwfPx4AMGnSJERHR2PatGmKdaKj\nozFnzhxYW1sjNDQUr732Gm/xEO1282YmOA6wstLerlP5hXceqKhIRGZmpqAKBREXpU6P3b17N89h\nALGxsS26kXx9fREVFdWiUMTExGDx4sWKZTs7O6SmpqJ3796ttndwwwY0H2g/fsjd4u8dvc5xT36/\nRNL6PRwHPFrmJBKA45CUlQW/Xr2gq6cHHV1d6OjqQldfX/53PT353/X0oKerCwN9fRgaGMj/1NeH\ngb4+jI2MYGhgIPqug/DoaAQHBmo6DNTW1SErMQ0PH1qivt4ZJSXq23dubjhcXILVtj/GPID6eGRG\nRmIojz/0OksobaGrpFIpaurqUPfwIeobGlDf0ICH9fXyvzc2ounRQ9rQgKaGBkibmiBtbERTUxOk\nTU1gMhkS0tPh4+YGxhiYTCZ/API/GWvz0aFHrys7v1ertTr4IdFuoRg4cCAiIiJgZWXV4e1OOY5T\n232z20pWe1+e+7//HpaPusoMdXXhaGoKT0tLAEDGozO4FMvl5d1fZqzd5dicHGTfuAFPc/OWr1tY\nAIwho6ICkMngaWb232XG4GlsDMhkyKiqgg6AQWZmMNPTQ1ZDA8yMjTHB1RU2BgZILC+HoY4Ogp2c\nAADh+fkAILhlAMDlyxqP59SdO6itqMfQYS/Dz08X9+6FAwD69QsGAF6XdXQAQH37c3T0QMx30biw\nNx3PPZrAU9P5D8/Px43SUgRfviyYeP66LGMMg6ytUVpfj/O5uaisq4Obri6q6utxs6ICDxmDp6kp\noKODjNpaQCKR///mOGRUVwMcJ///zXHIqKxsvfyITkaG/PuA41p/v1hZtVzu7vdVG9vLKC/HjUen\nT1s++n5qC8faKT+7d+9WdDd1dETBcRzeeOONdl9XVkVFBYKDgxWz0W7cuBFTpkxpcUTxySefoKmp\nCVu2bAEA9O7dG6mpqW3GlJub26KosDaqrTpfl8lk8l8TUimkUmnrvzc0oPHhQ9TX1OBhdbXiz4fV\n1aitrER9dTVQWwvU1MgftbVA8/TujMHcyAguNjZwtrWFu6UlXM3NoSPp1kltWuvsrVuIevAAwbt3\ni6qfuCtkMhneX78eD6OisGXKFFgYivtWr3ypbmhAelkZcsrKkFdcjILycjQ+9v8LBgaAiYn8YWwM\nzsQEJubmMDIzg6GJCQxMTWFgYgJDExPoGxtDz8hI3nOgqwsdHZ1Wf3KPeh86+2iLsj0Nyqzn4uLS\n5hFJu0cUjxcHdXQ9WVhYAJCf+eTu7o6wsLBWBSgwMBBbt27FkiVLcO7cOfj4+LS7PWdnZ17jVbeG\nhgZUV1ejqqoKlZWVeFBaipKsLJRmZ6MkOxuVRUWoLClB0qNJ4PR1dNDLwQF9nZ3hY2cHk0dnrBEg\ns7wc0NPrEad2SyQSuPv54V50NDIrKjCICgUA+Y+4vKoqJBQUIDU3F4UVFfIuYx0deRdM//6wcHaG\nrbs7bDw8YO3gACsrK5iZmcHMzAwmJiaQ9KAfYl2aZpwv+/btw5o1a9DY2IhNmzbB1tYWBw8eBACs\nWbMGAQEBGDt2LIYPHw5ra2scO3ZMwxErJzw8vNu/XPX19WFtbQ1ra+tWrzHGUFpaitzcXOSmpSE9\nPh7F9+7hbm4u7ubm4leOQ29HRwz08ICvvT10NdTAwzMyEOzpqZF9N3vY1ISCsjLoDByokak7VNEW\nOst90CDc09FBZmkpBjk4qHXf7dFUW3hQV4f43FzcSU9HWU0NIJEADg7Q7d8fHgMHwnPQIDi7usLJ\nyQnGxsa8xqKJttBVHY5RKEOVYxRBQUFISkpq8dyaNWtaLL/33nt47733VLI/bcFxHGxtbWFrayuf\nKfS551BRUYHU+/eRFBOD1OhopKSmIiUqCucMDDC8Xz+McHODaQ88ysiqqAAD4Ozj02PmP/Lw8gLs\n7JBVWAj4+mo6HLVjjCG9vBzR9+/jXl6efBDXxQVmo0bBLygI/QcMgJubG3R1BfW7WVA6HKNQagMq\nGqNQJY7jBHtnL02oqalBwp07uB4WhoLYWCArC7oSCQL69cNYLy8Y95AvTAAIS0nB1cREjHvnHTz1\nzDOaDkctpFIp3lu7Fo2xsXh16tQe1Q2ZXlaGi4mJyC4qAoyNoevjgwETJ8J/1Ci4u7v3qO4jZbT3\n3dluoRAzKhRtY4whKysLf4aFIfnsWSA1FQa6uhjn54dR7u49YvD7iytXkMNxWHTgAPr06aPpcNTm\n6EcfIf3rrzF31Cj42dtrOhzeldbW4tc7d5CamwuYmcF46FCMfO45DAsI4PV6MLFr77uz098M0dHR\niI6OVklQPYVQ5nThOA4eHh5YsHIlVv/f/6HP6tWod3LC+fh4HLh8WT7Iy6PwjAxet/8k9U1NyCsr\nA+fkBDc3zUwEqKm24DVkCKCjgzR1XjTSAb7aQpNMhov37+Oz8+eRWlICw9GjMeEf/8DmDz7A+JAQ\nQRUJoXwvKEOpQlFVVYVVq1bB0dERo0ePxujRo+Ho6IhVq1ah8rFzgol4ODs7Y9G6dVj80UewnjkT\nxTIZDoeH49zdu2hqPi1Qy2RWVEDGGFx8fWFgYKDpcNTKq29fwN4e6Y9f06JlCqqrcSgiAhG3b0Pq\n5YUh27Zh44cfYnxISI/791Y1pbqeFi5ciMjISLzzzjuYMGECOI7DhQsXsHPnTowcORInTpxQR6xK\no66nzmlqasKVS5dw5fBhyBIT4WhtjTnDh8OW57M+1O235GRE37uHoD17EKLCOcrE4PHrKTZPngwr\nIyNNh6QyjDFEZWfj/I0bkBoZwebppzFzxQq4u6t3Hi9t0N53p1LD/GfPnsXevXsRGhqqeC40NBSN\njY2Ki9+IeOnq6iJk4kT09fHB9wcPoiAsDIcuXsTswED0t7PTdHgqk5afD9jZwatfP02HonYSiQSe\n/v5IjopCenm51hSKRqkUP92+jTtpaUCvXhi+ZAkmzZihmOmaqIZSXU+urq6wenT59+MsLS217sI2\nPoilL9LV1RVrX3sNAzZsQIO1Nf79xx+4kp6usqMzTY5RVNbXo7iiAvru7hq99akm24LXkCGAgQHS\nCgs1FkMzVbSF8ocP8eXVq7iTmQn9kSMx/913MX32bNEUCbF8LwBKHlG89tpreOutt+Do6IiAgAAA\n8gn63nrrLbz++uu8BkjUy8DAALMXLYKDmxsu/vOfuBAfjwc1NZjh5weJiCclTCsrk8+nM3gwdOQT\nLvU4Xr17A87OSM/PF/39KYpqavCvq1dR1dQE6xkzsOCll2DfA87m0hSlxigGDhyIzMxMVFdXw/zR\nxFGVlZUwNTVtMQ2COicI7AiNUajG3eRkfPf++2i8dg0+Hh6Y7e+vsau6u+s/N2/iVl4epnz0EUaO\nHq3pcDSCMYaPtm1D5YULWPPUU3AyM9N0SF2SXVGB41ev4qGeHjzmz8eCVatgpCVdaZrWrTGK2bNn\nK70Toj36e3tjyZtv4viePUj64w+caGpC6LBh0BPZL3LGmHx8wskJvfv21XQ4GsNxHPoEBOD6hQu4\nX1oqykKR/uABTkRGotHUFN5LlmDOkiV0RbUaCOZ+FNpMTHO6/JWbmxuW7dqFf+3Zg7SLF/ENgBeG\nD+/SkYWm5vfJr65G9cOHMO/du8UdEzVB022h76BBuG5piXu5uRinwXm3utIWMsvKcOLqVTRaWMB/\n9Wo8O2+eqK+s1nRb6IxOZ7m+vh61tbUtHkS7OTg4YOn//A9Mn3oKaXl5+Pb6dUhFdK3FveJigOPQ\nb+TIHn/U6+XlBR13d+SUlKC2sVHT4Sgt61F3U6OFBYasXYuZ8+eLukiIjVKZLioqwpYtWzBgwAAY\nGRnB1NRU8TAT4eGruonlV0NHbG1tseTvf4dxUBDuZWfj+5s3IevkOJCmZo69l5sLWFmhv7+/Rvb/\nOE23BQMDA3gMGQLGGO4/eKCxODrTFgqrq3H86lU0mJpi8KpVmDFnjlYUfE23hc5QquvpvffeQ1hY\nGJYsWYIBAwZQJe+h7O3tsXj7dhytr0fin3/idyMjTHns9rVCVFVfj7zSUugNGwZPDU9xLhR9AwKQ\n9s03SCkoEMy04+2pePgQxyMjUa+vD9+lS+lIQkOUKhQnTpzAN998g5CQEL7j0Upi6ot8EicnJyzY\nvh3/2rULUTdvwsLICKOUvAGQJsYo7j2a28hr2DBBTCsuhLbQz9sb51xccD83F7JBgzRy2rMybeFh\nUxOOR0WhUiqFx4IFeH7RIq0qEkJoC8pSKuvDhg1Deno637EQkfD09MSsV18F+vTB7/HxSCwq0nRI\n7bqXnw8YG6Pfo+t/CGBjYwNrb2/UNTQgq6JC0+G0SSqT4d9xcSiqrITd9OlYsGoVnd2kQUpdR5Ga\nmorVq1djxowZmDRpEkxNTVu8LrQ5Veg6CvX4Izwc5997D7rFxVgWFASXDm7OrgkNUik++PVXNPbp\ng62ff664BogAYb/8gqu7diGgXz9MFWD34ZmEBMTduwezp5/GytdfV9wqmfCrW9dRSCQS1NbWYuvW\nrW1uWCqVdj9CIjpjgoJQVlKCax9/jJNRUVgdHCyou+bdKy1Fo1QKtyFDqEj8hd+QIbjq5ISkrCw8\n07+/oAaH43NzEXf3LnSHDUPoK69QkRAApbqe1qxZA6lUii+++AJRUVGIiYlRPOjeFE8mpjldOoPj\nOEx97jm4z56NyoYGfBsX1+Fps+qe6ykhNxcwMIDf+PFq3W9HhNIWnJycYOnjg6raWmRr4FYB7bWF\n3IoKnImLA7y8MG3jRq2eS04obUEZSh1RXL9+HT/++CPGjh3LdzxEZHR0dDBv5UocKihA1m+/4dek\nJMzw89N0WKhvakJKTg7g5QXfAQM0HY7gcBwH36AgRF64gMSCArgL4Fd7TUMDTv75J6RWVhixdCmG\nDB2q6ZDII0odUYwcORI3btzgOxatJZYzG7rK1NQUCzZtgu7Qobh29y7icnLaXE+dZzzdKy1Fk1QK\n9+HDBdXtJKS24DtkCODsjMTMTLWP6f21LUhlMpyKi0MlY3CbNQtTZsxQazyaIKS28CRKHVGsX78e\nW7duRUFBQZuD2UOp8vd4zs7OmLF5M37YtQu/XrsGezMzjf5KvZOTAxgaCqrbSWhcXFxg4euLirAw\nZFZUwNPSUmOxhN29i4yiIphOmYJ5y5b12Bl+hUqpI4rp06fj3r17eOeddxAcHIzhw4crHiNGjOA7\nRtETU19kdwz298fIlSshMzPDt9HRqKqvb/G6usYoahoacD83F5wAu52E1BY4jsPgp58GdHVxo52j\nQL483hZu5+cjKikJEn9/zFu/vsfM9iCktvAkSh1RpKWl8R0H0RITp09HQUYGMr7+Gt/Gx2NpQAB0\n1HyR1M2CAkhlMvQbO7bHfOl0lf+IEYhwd0dCRgae8faGgZqvVSiorsbp2FjA0xPPrF8vuFPtiZxS\nraJ56oPKykpcuHABVVVVbb5O2iamvsju0tHRwdyVK3EwKwvZYWE4e/cupvn4AFDPGAVjDNfT0gAr\nKwydMIH3/XWW0NqCtbU1PEaNQmZaGhKKijBUTWcZBXt6oq6xESejotBoYgL/hQsxvIf1TgitLXRE\nqUKRk5ODZcuW4cqVK2hoaICOjg6kUil0dXVhYGCAJUuW8B0nERETExPM37oVh/PyEJuYCGcLCwxR\n0xdQdmUlSsrLYTp+PPr24HtPdMaQp59G5unTiM/IUFuhkDGG7+LjUVZXB+cFCzB99mxBXctBWlKq\nT+Do0aPQ1dVFYmIiTExMEB8fj5MnT2LgwIE4f/483zGKnpj6IlXFxcUF07ZuBezt8cv168irqlLL\nGMX1rCxAVxf+U6YIckBUiG3B188PBr6+yC4qQkF1tVr2+dHVq0jNyYHx2LGYv3p1j5yeQ4htoT1K\nFYpjx45hy5Yt8rnsdXTQ1NSEuXPnYt26dXjllVf4jpGI1JARIzB85Uo06eriZGws6ni+/0F1QwPu\nZGQAXl4YOnIkr/vSJvr6+vCfMgXQ0cGfapjTLbGoCLfT08H5+GDuyy/TldcioFShKCwsRO/evQEA\ngwcPRnJyMgD5abGJiYn8RaclxNQXqWrPzJkDt1mzUFFZicLy8k7fw6IzYrKz0SSVov9TT8Ha2pq3\n/XSHUNvCqJAQcL1743Z6OioePuRtP4XV1fgxNhaevXtj4oYN6NWrF2/7EjqhtoW2KFUo/Pz8cOvW\nLQDyD7d371788MMPeP3110X1YYn66ejoYN6GDTAdNQoZ2dkIS0nhZT8NUili794FXF0x5plneNmH\nNrO0tITflCmQMYborCxe9lHb2Ih/x8SgQVcXg5Yswahx43jZD1E9pQrFzp07FVdurly5Er6+vggN\nDcWDBw+wfft2XgPUBmLqi+SDmZkZ5m3fjkwLC/yZmIjbBQUq30dMdjbqGhrgFhQENzc3lW9fVYTc\nFkZPmQK4u+NaSorKb5MqYwyn4uNRVlUF55kzYe7o2OMHr4XcFv5KqUIxZcoUPP/88wAANzc3HD16\nFGVlZYiMjEQAzfNPlODu4YHABQsAW1ucjotT6aDpw6Ym/JGYCDg7I4TOnukyZ2dn9J44EfWNjbii\nwmunGGM4m5yM9JwcmI4fjwXr1/fIwWsxU+p+FGJD96MQJsYYfvrmG9z49FOYMoYV48bBysio29v9\n/d49RCYmoteyZViyaRMVim7Iz8/HwZdegk5GBl56+mmV/PtcycjAhevXoTtoEF58911BH/H1dO19\nd2rPfQWJ4HEch+nz5qHX/PmobmzEv6KiUNPQ0K1tFlZXIyopCVyfPpg4dy4ViW5ycnLCoOeeg5Tj\n8GtCQrd/cN3Iz8eF+HhwXl54fvt2KhIiRYVCDcTUF8mn8PBw6OrqYsGaNXCaNQsPqqpwLCamy6fN\nSmUynL55EzJ9fYx44QVR3LtADG1h0qxZMBw2DCnZ2bhdWNjl7SQUFuJ0TAzg7Ixn/va3FvNuiSEP\nfBNTDqhQELUzMDDAws2bYT11KvIfPMDRLh5ZXExNRW5RESzGj8eEqVN5iLRnMjU1xeTVqwEHB/xy\n7RpKams7vY1bBQX4LioKMjs7jN+8GQGjR/MQKVEXGqMgGlNRUYGvP/wQpb/+ChtTUywICICdiYlS\n743Py8NP0dGQ9O+PpR9+SJPJqRhjDN8dPYqE/ftha2SEZaNHw0SJ29wyxnA1MxMXbtwAc3BA0Nat\nCJ40iboERaK9704qFESjampq8K+PP0bB6dPQl0oxY9gwDHBw6PCLJT4vD6djYsCcnDB11y76tcqT\nhoYGfLFnD4p+/BEOFhZYGBgIcwODdteva2zEmcREJKSlAZ6eeHrLFoyl66xEhQazNUhMfZF8aisP\nJiYmWL5tGwauW4cGCwt8/+ef+Ob6deT/ZYZiQH6fiZ8SEvBTdDSYkxOCt2wRXZEQU1vQ19fHki1b\nYDttGgorK3EwPBx3CgtbfZE0yWSIy83FpxcvIiEjAwYjRiB0z54Oi4SY8sAXMeWATmYmGqevr4/n\nlyyBh7c3wr74Avdu3MC97GzYW1rC1c4Oerq6eFBZibRH95nQGTgQUzdtwjC6hod3pqamWPb3v+N7\nOzuk/fADvouMxHlTU3g6OsLQwAAV1dVIy89HfWMjYGcHj3nz8OySJbCxsdF06ESFqOuJCEp1dTWu\nXryI67/+ivq0NKC8HJBKAWNjcC4u6BcSgomzZ8PW1lbTofYoMpkM8deuIeLkSVQkJgIlJfJ/F0ND\nwMEBTkOHYuysWfD186PxCBGjMQoiKk1NTcjOzkZpaSkaGhpgbm4OT0/PVvdrJ+olk8mQl5eHgoIC\nNDQ0wMzMDC4uLoKdhJF0jqALRVVVFRYtWoT4+HgMHToUx44da/MLwdPTE+bm5tDR0YGenh5iYmLa\n3J7QCkV4eDhNngjKA0A5aEZ5EGYOBD2YvX//fri7uyMlJQWurq44cOBAm+txHIfw8HDEx8e3WyQI\nIYSoliCOKObMmYPXXnsN/v7+uH79Ot59912cOnWq1Xq9evVCXFzcEwfKhHZEQQghYiDoI4rY2Fh4\ne3sDALy9vTvsUpowYQJmzZqF06dPqzNEQgjpsdR2euzEiRNR0MZ9CN5++22lf/1fvXoVTk5OSEpK\nwowZMxAQEABHR8c21126dCk8PT0ByG/K4u/vr+gPbD5/WV3L+/bt0+j+hbLc/JxQ4tHE8l9zoel4\nNLV848YNvPzyy4KJRxPLzc9puj0eOXIEABTfl21iAvD888+z69evM8YYi4uLY7Nnz37ie7Zs2cIO\nHTrU5msC+VgKly5d0nQIgkB5oBw0ozwIMwftfXcKYozi/fffR3Z2Nt5//31s27YNvXr1wrZt21qs\nU1tbC6lUCjMzMxQXFyM4OBhnz55tc9piGqMghJDOE/QYxbp165CVlYX+/fsjNzcXa9euBQDk5eVh\n2rRpAICCggKMGzcO/v7+WLBgAV555RWa254QQtRAEEcUqia0I4pwAZ4vrQmUB8pBM8qDMHMg6CMK\nQgghwkVHFIQQQgDQEQUhhJAuokKhBo+fN92TUR4oB80oD+LKARUKQgghHaIxCkIIIQBojIIQQkgX\nUaFQAzH1RfKJ8kA5aEZ5EFcOqFAQQgjpEI1REEIIAUBjFIQQQrqICoUaiKkvkk+UB8pBM8qDuHJA\nhYIQQkiHaIyCEEIIABqjIIQQ0kVUKNRATH2RfKI8UA6aUR7ElQMqFIQQQjpEYxSEEEIA0BgFIYSQ\nLqJCoQZi6ovkE+WBctCM8iCuHFChIIQQ0iEaoyCEEAKAxigIIYR0ERUKNRBTXySfKA+Ug2aUB3Hl\ngAoFIYSQDtEYBSGEEAA0RkEIIaSLqFCogZj6IvlEeaAcNKM8iCsHVCgIIYR0iMYoCCGEAKAxCkII\nIV1EhUINxNQXySfKA+WgGeVBXDmgQkEIIaRDNEZBCCEEAI1REEII6SIqFGogpr5IPlEeKAfNKA/i\nygEVCkIIIR2iMQpCCCEAaIyCEEJIF1GhUAMx9UXyifJAOWhGeRBXDqhQEEII6ZAgxihOnTqF3bt3\nIzk5GbGxsRg6dGib60VERGDNmjVoamrCpk2bsHHjxjbXozEKQgjpPEGPUQwcOBA//PADxo8f3+F6\nmzdvxsGDB3H+/Hl8+umnKCkpUVOEhBDScwmiUHh7e6Nfv34drlNRUQEAGD9+PDw8PDBp0iRER0er\nI7xuE1NfJJ8oD5SDZpQHceVAEIVCGbGxsfD29lYs+/r6IioqSoMREUJIz6Crrh1NnDgRBQUFrZ5/\n5513MGPGDHWFoRHBwcGaDkEQKA+Ug2aUB3HlQG2FIiwsrFvvHzFiBF599VXFckJCAqZMmdLu+kuX\nLoWnpycAwNLSEv7+/op/mOZDPlqmZVqm5Z68HB4ejiNHjgCA4vuyTUxAgoODWVxcXLuv+/v7s8uX\nL7P09HTWv39/Vlxc3OZ6AvtY7NKlS5oOQRAoD5SDZpQHYeagve9OQYxR/PDDD3Bzc0NUVBSmTZuG\nZ555BgCQl5eHadOmKdbbt28f1qxZg6effhrr16+Hra2tpkImhJAeQxDXUagaXUdBCCGdJ+jrKAgh\nhAgXFQo1aB486ukoD5SDZpQHceWACoUa3LhxQ9MhCALlgXLQjPIgrhxQoVCD8vJyTYcgCJQHykEz\nyoO4ckCFghBCSIeoUKhBRkaGpkMQBMoD5aAZ5UFcOdDK02ODg4Nx+fJlTYdBCCGiEhQU1OYgu1YW\nCkIIIapDXU+EEEI6RIWCEEJIh6hQ8KCqqgozZ86Eu7s7Zs2aherq6lbrZGdnIyQkBH5+fggODsaJ\nEyc0ECm/lMkDACxfvhwODg4YOHCgmiPkT0REBHx8fNC3b1988sknba6zY8cOeHl5YdiwYUhOTlZz\nhOrxpDwkJydj1KhRMDQ0xN69ezUQIf+elIPjx49j8ODBGDx4MF544QXcu3dPA1E+gZomJexR9uzZ\nw1566SX28OFDtmHDBvbBBx+0Wic/P5/Fx8czxhgrLi5mvXr1YpWVleoOlVfK5IExxiIiItj169fZ\ngAED1Bwhf5pnOs7IyGhzpuPo6Gg2ZswYVlpayk6cOMGmTZumoUj59aQ8FBUVsdjYWLZz50724Ycf\nahPuyZcAAAuzSURBVChKfj0pB5GRkay8vJwxxtiRI0fYokWLNBFmh+iIggcxMTFYsWIFDAwMsHz5\n8jZv2ero6Ah/f38AgK2tLfz8/BAXF6fuUHmlTB4AYNy4cbCyslJzdPxR5ra90dHRmDNnDqytrREa\nGoqkpCRNhMorZfJgZ2eH4cOHQ09PTxMh8k6ZHIwaNQoWFhYAgGnTpgnyjE0qFDx4/Lat3t7eiImJ\n6XD9+/fvIyEhAQEBAeoIT206mwdtocxte2NiYuDr66tYtrOzQ2pqqtpiVAe6fXHnc3Do0CFB3vFT\nbXe40zbt3dr17bff7tQU51VVVZg/fz4++ugjmJiYqDJEtVBVHnoaxlir/HAcp6FoiBCcP38ex44d\nQ2RkpKZDaYUKRRd1dGvXo0ePIikpCUOGDEFSUhJGjBjR5nqNjY2YPXs2Fi9ejJkzZ/IVKq9UkQdt\no8xtewMDA5GYmIjJkycDAIqLi+Hl5aXWOPnW2dsXayNlc3Dr1i2sXbsWZ8+ehaWlpTpDVAp1PfEg\nMDAQX331Ferq6vDVV19h5MiRrdZhjGHFihUYMGAAXn75ZQ1EyT9l8qCNmvubIyIikJGRgbCwMAQG\nBrZYJzAwEN9//z1KS0tx4sQJ+Pj4aCJUXimTh2baevSpTA6ysrIwe/ZsHD9+HH369NFEmE+myZF0\nbVVZWcmeffZZ5ubmxmbOnMmqqqoYY4zl5uayqVOnMsYYu3LlCuM4jg0ePJj5+/szf39/9ttvv2ky\nbJVTJg+MMbZgwQLm5OTE9PX1maurK/vqq680FbLKhIeHM29vb9a7d2/28ccfM8YYO3DgADtw4IBi\nnb///e/M09OTDR06lCUmJmoqVF49KQ/5+fnM1dWVmZubM0tLS+bm5qZoJ9riSTlYsWIFs7a2VnwP\njBgxQpPhtomm8CCEENIh6noihBDSISoUhBBCOkSFghBCSIeoUBBCCOkQFQpCCCEdokJBCCGkQ1Qo\nCBGwpUuXanTun2vXrsHV1RUymeyJ6+bm5sLKygplZWVqiIyoExUKIhjBwcHYuHGjpsPQiPDwcEgk\nEjx48KDF85988gmOHz+uoaiAPXv2YO3atZBInvxV4eLigqeeegqfffaZGiIj6kSFghAB+ev1r2Zm\nZjA3N9dILPfv38eZM2ewevVqpd+zYcMGfPLJJ6ivr+cxMqJuVCiIICxduhQRERH49NNPIZFIoKOj\ng8zMTIwcORJ79uxRrLdo0SJIJBIUFhYCAGpra2FgYKCYcbOxsRGvvfYa+vTpA0tLS0yYMAE3b97s\ncN8LFizAunXrFMuvvfYaJBJJi/sGuLm5Ke5CGBsbi0mTJsHOzg6enp7YsWMH7t69q1h3+fLlrbqL\nZDIZ3N3dsW/fvlb7z8jIwIQJEwDIpxuXSCRYvny5Ii+Pb6v5qOvNN99Er1694O3tja+//hoymQx/\n+9vf0KtXL0ycOBEREREt9lFaWoply5bBw8MD/fv3x549exQ5bM/hw4cxduxY2NvbK54rLi7Gpk2b\n4OvrCxMTE/j4+ODkyZOK10NCQiCTyfDzzz93uG0iLlQoiCD83//9H0aNGoXly5ejoKAA+fn5cHNz\nQ0hICMLDwxXrXb58GXZ2dornIiMjoaenp7iXx0cffYQDBw5g7dq1OH36NPT19TFhwgRIpdJ29/3X\nfYSHh7fYx/3795Gbm4vg4GAAQHV1NV588UX88ccfOH78OLKyslp8ma9evRpnz55tMf16WFgYCgsL\nsXjx4lb7d3d3x/fffw8ASExMREFBAT7++GMA8qnH/zr9+IkTJ5CTk4NTp05h+vTpWL16NUJDQyGT\nyfDdd9+hV69eeP755xVHJ01NTQgKCkJlZSU+//xzHD58GLGxsZg1a1YH/yLAnTt3Wk1S9+GHHyI6\nOhqHDh1CQkIC9u3b12q20z59+uDOnTsdbpuIjGanmiLkv4KDg9nGjRtbPPfbb78xU1NTJpVKWUpK\nCjM3N2e7du1ia9asYYwxtnPnTjZx4kTF+oGBgeztt99WLNfV1TEjIyMWFhbW7n6TkpIYx3GsoKCA\n1dTUMAMDA/b++++zyZMnM8YY+/zzz1nfvn3bfX9lZSUzMzNjd+7cUTw3YMAA9t577ymW582bx+bO\nndvuNi5dusQ4jmOlpaUtnn/xxRfZ9OnTFctBQUHM1dVVsfzw4UOmr6/PBg8erHiusLCQcRzH4uLi\nGGOMHTt2jNna2raYbK+kpIQZGRmxtLS0dmPq169fq9uTBgcHs1dffbXd9zDG2OLFi9m8efM6XIeI\nCx1REEEbO3Ys6uvrERMTg/DwcIwbNw5PPfWU4td+eHi44pd+fn6+oluomaGhIYKCgvCf//wHALB2\n7VqYmZkpHoD87nuOjo64dOkSIiMj0adPH8ybNw9Xr15FU1NTi30AQFlZGd59911MnjwZNjY2cHZ2\nRm1tLc6dO6dYZ9WqVTh8+DAA4MGDBzh9+jRWrFjR7XxwHIeJEycqlg0MDNCnT58Wn9ne3h52dnaK\nW+vevHkT5eXlcHJyUnxuT09PNDQ04Keffmp3Xzk5OXBycmrx3MaNG3Hw4EH4+/vjzTffRG5ubqv3\nOTk5ITMzs7sflQgIFQoiaKamphg2bBguXbqEy5cvIyQkBCNHjkRWVhZSU1MRFxfX4ku8LYwxRffN\nP/7xD9y8eVPxaBYUFNRiHx4eHrC1tUVsbCwiIiJa7GPfvn344osvsGjRIpw7dw43btyAjY0Nqqur\nFessWrQI/7+9+wlJJIrjAP61hCCIqKB/lIaUXlRMwsiKXgRBEFQTGBGRhw4djEi6RAgdgiA7BFFY\nEF4UomMX6WIqCEIEFZ1CzUhBKCLoHxHiHpbcHbVp2XW3dvf3AQ8z88b3ex78zbw/vIuLC/j9fjid\nTpSXl6c2KcrFb/K9vLy8rOdeu9sSiQQkEgmv3cfHxzg7O8PY2Nib9UgkkoxEwHEcotEoTCYTAoEA\nGhoaYLPZeGVisRjq6up+oYXks6FEQT6NoqIiPD4+ZpxnjMHtdsPr9YIxhoKCAjQ3N2NhYYE3PlFV\nVQWdTsd7sn96eoLP5wPHcQC+DhbLZLLU5/s69vf3eW8PjDFsbm4iGo3yEoXdbofJZMLo6Ciamppw\nc3ODq6srXsylpaXgOA5bW1uw2+2Cf8ivbX+NN9c0Gg3Oz8/x/PzMa7tMJkNJScmb9ykUCoTD4ayx\njo+Pw+VywWw2Y21tjXc9FAr9kxsx/c8oUZBPo6OjA36/H6enp7i+vk4NxjLG4PF4cHd3B61Wmzrn\ncDjQ0tICsfjbjr4cx2FlZQXLy8vwer3o7+9HYWEhOjs7BetmjCEYDOLg4ICXKBwOB+rr61FdXZ0q\n29bWhu3tbQQCAezs7MBkMmU80QNfu5+cTidOTk5Ss5jeotFoUFxcDKfTiXg8joeHh6zlkln22k4/\nTjc8PAy1Wo2RkRFsbGwgGAzC5/NhZmYGwWDwzfuUSmXG9ampKezt7SEUCsHj8cDlcqGnp4dXJhQK\nQalUCsZE/i6UKMinYTAYoNPp0N3djYqKClxeXgIAWltbIRKJ0N7enupCYowhkUhkdDtNT09jYmIC\nNpsNfX19eHl5gdvtfnfBmEKhQGVlJeRyOcrKygTrmJ+fh1wuR29vL6xWKywWC6RSacZ3MsZQU1MD\nxti7XTH5+flYWlrC7u4uamtrUwsP02c9ZZsFlX6cTiwWw+PxQK/Xw2q1QqVSwWg04vb2VvCNwmg0\nwu/386bRJpNJTE5OQq1Ww2w2o6urC7Ozs6nrbrcbIpHoQ1eTk9yjHe4I+U3u7+8hlUqxvr6OoaGh\njw7npxgMBqhUKlgslh8qPzg4CK1Wi7m5ud8cGfmTKFEQkmPJZBLxeByLi4vw+Xw4Ojr66JB+2uHh\nIQYGBhCJRN59K4vFYlCpVAiHwxlrK8jfjRIFITkWiUQgk8nQ2NiI1dVV6PX6jw6JkF9CiYIQQogg\nGswmhBAiiBIFIYQQQZQoCCGECKJEQQghRBAlCkIIIYIoURBCCBH0BavI3/TYY7zNAAAAAElFTkSu\nQmCC\n", "text": [ "" ] } ], "prompt_number": 54 }, { "cell_type": "code", "collapsed": false, "input": [], "language": "python", "metadata": {}, "outputs": [] } ], "metadata": {} } ] }