diff --git a/.gitignore b/.gitignore index 14e3cc80..70c52100 100644 --- a/.gitignore +++ b/.gitignore @@ -38,3 +38,4 @@ nosetests.xml *.sublime-project *.sublime-workspace docs/_build/ +FDEMMagDipoleWholeSpace-checkpoint.ipynb diff --git a/notebooks/FDEMMagDipoleWholeSpace.ipynb b/notebooks/FDEMMagDipoleWholeSpace.ipynb new file mode 100644 index 00000000..3b546a88 --- /dev/null +++ b/notebooks/FDEMMagDipoleWholeSpace.ipynb @@ -0,0 +1,184 @@ +{ + "metadata": { + "name": "", + "signature": "sha256:3527637bdb920eb3c32f99dc1242bd648d175b1ce5e156380d61f57c52af9874" + }, + "nbformat": 3, + "nbformat_minor": 0, + "worksheets": [ + { + "cells": [ + { + "cell_type": "code", + "collapsed": false, + "input": [ + "from SimPEG import *\n", + "import simpegEM as EM" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "output_type": "stream", + "stream": "stdout", + "text": [ + "Efficiency Warning: Interpolation will be slow, use setup.py!\n", + "\n", + " python setup.py build_ext --inplace\n", + " \n" + ] + } + ], + "prompt_number": 1 + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "cs = 5\n", + "ncore = 20\n", + "pad = 5\n", + "padfactor = 1.3\n", + "\n", + "h = [(cs,pad,-padfactor),(cs,ncore),(cs,pad,padfactor)]\n", + "mesh = Mesh.TensorMesh([h, h, h],'CCC')" + ], + "language": "python", + "metadata": {}, + "outputs": [], + "prompt_number": 18 + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "x = mesh.gridCC[:,0]\n", + "y = mesh.gridCC[:,1]\n", + "z = mesh.gridCC[:,2]\n", + "sig0 = 1e-2 \n", + "f = 100 " + ], + "language": "python", + "metadata": {}, + "outputs": [], + "prompt_number": 19 + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "Bx,By,Bz = EM.Analytics.FEM.AnalyticMagDipoleWholeSpace(x,y,z,sig0,f)" + ], + "language": "python", + "metadata": {}, + "outputs": [], + "prompt_number": 21 + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "mesh.plotSlice(Bx.real,grid = True)\n", + "mesh.plotSlice(Bx.imag,grid = True)\n", + "mesh.plotSlice(By.real,grid = True)\n", + "mesh.plotSlice(By.imag,grid = True)\n", + "mesh.plotSlice(Bz.real,grid = True)\n", + "mesh.plotSlice(Bz.imag,grid = True)" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "metadata": {}, + "output_type": "pyout", + "prompt_number": 28, + "text": [ + "(,\n", + " )" + ] + }, + { + "metadata": {}, + "output_type": "display_data", + "png": "iVBORw0KGgoAAAANSUhEUgAAAYkAAAEZCAYAAABiu9n+AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAFOhJREFUeJzt3X/wZXV93/Hna/mpCxEQs6xIIbbAiOmoRRcZbdnJjy3p\nTETSCeJMEkyp42gkmWmSiT+mKYkzNGQaJk2j/lFRaaKQnbaYTSsJCy2JGiliU0MloDiuZXdgQcUq\nRnBh3/3jHsqX3fu59/u993vvPfu9z8fMd773ns857/O55977fX/PeX/OOakqJEkaZtOiOyBJ6i+T\nhCSpySQhSWoySUiSmkwSkqQmk4QkqckkIa1Bkrck+dSK599JctbieiTNlklCOkSS1yf5yyTfSvKN\nJJ9O8uph81bViVW1Z53X/84kdyd5IslHDmk7K8nBLjk98/Pe9Vy/tNLRi+6A1CdJfgD4L8DbgJ3A\nccA/BJ6cYzf2Ae8D/jHwvMY8P1CeCas5cE9Ceq5zgKqqP6qBJ6pqd1XdM2zm7r/6l3aPn5fkd5Ls\n6fZCPpXk+K7ttd3eyWNJ/leSi1odqKqbq+qPgW+M6KffXc2FHzTpue4Hnk7y0SQXJzl5Dcv+G+BV\nwIXAKcCvAgeTnM5g7+Q3q+pk4FeA/5Tk1DHxMqLta0keTPLhJC9cQx+lNTFJSCtU1XeA1wMF/Hvg\nkSR/nOQHRy2XZBPw88AvVdVDVXWwqu6squ8DPwN8sqr+tFvHbcDdwD8Z150h0x4FXg38HeB84ETg\nY6t+gdIamSSkQ1TVfVX181V1BvDDwIuB3x2z2KnA8cBXhrSdCfx0d6jpsSSPAa8DThsT87A9iar6\nblX9zy4JPQK8E9iRZPOYWNJETBLSCFV1P3ADg2QxyteBJ4C/N6Tt/wB/UFUnr/g5sap+e9zq19BV\nv8uaCT9Y0gpJzk3yL7o6AknOAN4MfHbUclV1EPgwcF2SrUmOSnJhkmOBPwR+MsmObvrxSbY/s44h\nfTiqK3gfDRyV5LgkR3Vt27o+bupqEb8H/PfuMJm07kwS0nN9B7gA+B9JHmeQHP4a+OWuvXjuf/gr\nH/8KcA/wOQYjk/41sKmq9gKXAO8BHmGwZ/HLtL9//xL4W+DXGNQzvgc8cy7ES4FbgG936/oegyQm\nzUQcai1JanFPQpLUZJKQJDWZJCRJTSYJSVLThrvAXxIr8ZI0gao67ATODZckAPjiGvPE+6+GX7h6\nsnW9PGtf37ziThJjtcuMm2+a9iM19iLXvayx1zrftMvMIsZ6x53079nLh18qzMNNkqQmk4QkqWmh\nSaK7zPH+JPesmHZKkt1JvpTk1iQnrWh7d5IvJ7kvyY5168hrtq9bKElaqHX+e7boPYmPABcfMu1d\nwO6qOge4vXtOkvOANwHndct8oLs88/S2bV+XMJK0cOv892yhSaKqPgU8dsjkNzC46ibd7zd2jy8B\nbqyqA909hR8Ats2jn5K0rBa9JzHMlqra3z3eD2zpHr8Y2Ltivr3A0KtoSpLWRx+TxP/X3eh91Dgw\nz4mQpBnq43kS+5OcVlUPJ9nK4NLKAPuAM1bM95Ju2uHef/Wzj1+z3ZqDJB3qrjvgc3eMna2PSWIX\ncAVwbff7EyumfzzJdQwOM50N3DU0wqQnxknSsti2/bn/QH/gN4bOttAkkeRG4CLg1CQPAr8O/Baw\nM8mVwB7gMoCqujfJTuBe4CngHeXNMCRpphaaJKqqdUetH2vMfw1wzex6JElaqdeFa0nSYpkkJElN\nJglJUpNJQpLUZJKQJDWZJCRJTdlopxp4+1JJmoy3L52FPt6+cJoY3r70yFz3ssZe63zTLjOLGPOM\nO26dQ3i4SZLUZJKQJDWZJCRJTSYJSVKTSUKS1GSSkCQ1mSQkSU0mCUlSk0lCktRkkpAkNZkkJElN\nJglJUpNXgZUkAV4Fdjb6fBVIrwI739iLXPeyxl7rfNMuM4sY84w7bp1DeLhJktRkkpAkNZkkJElN\nJglJUpNJQpLUZJKQJDWZJCRJTSYJSVKTSUKS1GSSkCQ1ee0mSRLgtZtmo8/XbvHaTfONvch1L2vs\ntc437TKziDHPuOPWOYSHmyRJTSYJSVKTSUKS1GSSkCQ1mSQkSU0mCUlSk+dJSJIAz5OYjT6Pk/Y8\nicmW/cKI9ldM0T7NsquJvVHfD8+TmA/Pk5AkrVVv9ySS7AG+DTwNHKiqbUlOAf4IOBPYA1xWVd9a\nWCclaYPr855EAdur6lVVta2b9i5gd1WdA9zePZckzUifkwTAoQfJ3gDc0D2+AXjjfLsjSculz0mi\ngNuS3J3krd20LVW1v3u8H9iymK5J0nLo7RDYJFur6qEkLwJ2A1cBu6rq5BXzfLOqTjlkuX6+IEnq\nuSNqCGxVPdT9fjTJzcA2YH+S06rq4SRbgUeGLuwQ2MljOATWIbDrueysY691vmmXmUWMecYdt84h\nenm4Kcnzk5zYPd4M7ADuAXYBV3SzXQF8YjE9lKTl0Nc9iS3AzUlg0MePVdWtSe4Gdia5km4I7OK6\nKEkbXy+TRFV9FXjlkOnfBH5s/j2SpOXU28L1pCxcS9JkjqjC9VQsXE8eYxkK17MqHo9rX3Tsvr4f\nFq7nF3fcOofoZeFaktQPJglJUpNJQpLUZOFakgRYuJ6NPheuLFwPb7NwfTgL15Pp8/d/knUO4eEm\nSVKTSUKS1GSSkCQ1WbiWJAEWrmejz4UrC9fD2yYtDj8FnB/4/IjlR7VPs+zRY/pm4drC9XqscwgP\nN0mSmkwSkqQmk4QkqcnCtSQJsHA9G30uXFm4nmzZSQvXB4DXBu5stI9qO2ZM7NUUrjfq+2Hhej4s\nXEuS1sokIUlqMklIkppMEpKkJkc3SZIARzfNRp9HNyxqdNPRwLmB+0fEGdU+zbKriT3p6KbHgdcF\nPtNofwL40cDtjfZRbcePiX0C049uOlLfD0c3zYejmyRJa2WSkCQ1mSQkSU0WriVJgIXr2ehz4WrW\nhetZFjO/crC97N/dNHn7apZtvf6vAxcF/rzR/hBweeCmRvuotq1jYp/K+MtbLHKbzbIobuF6Pixc\nS5LWyiQhSWoySUiSmixcS5IAC9ez0efC1aTFvFGFxGdMU5g+uqYvlH7twNCmTcc9ycHTTmDTw48P\nbT942gkc+43/2wz9/Re+oL3s/ZtHF5dvAj4YeHujfVTb5YyMvenc787udT15HJx5THObcuYxkxe9\nn8psz+ZeOZ+F6+nXOYSHmyRJTSYJSVKTSUKS1GThWpIEWLiejT4Xro7kwnWriApw5jHNIuzzT/we\nj29+ESd899Gh7Y9vfhFn1t80Q38tL+P1devQtqc5ms/mR7iw/tvQ9s/mL4Cru59h2m0X1j8aGfso\nnuLT2dHs26ezY+zram2Tv/3O88YWxce9Hxau15mFa0nSkcAkIUlqMklIkposXEuSgA1SuE5yMfC7\nwFHAh6rq2sNmsnA9eYz1KFyPK0zD+ML1vifay55+fPPs4pNO+RaP5Ex+sL42tP2RnNksDgN8Nj/C\nR+uyoW17snN0WfqLNf5y3o22q1+ekbHPqst4S3Y2+/aW7Bz7ulrb5FvfPInvv/AFzW36/Re+YOz7\nMdV7beF6fnHHrXOII+pwU5KjgN8HLgbOA96c5GWL7ZUkbVxjk0SSX0xy8jw6swrbgAeqak9VHWBw\ntZxLFtwnSdqwVrMnsQX4XJKdSS5OMnyfZD5OBx5c8XxvN02SNAOrKlwn2QTsAN4CvBrYCVxfVV+Z\nae8O78c/BS6uqrd2z38GuKCqrloxj4VrSZrAxIXrqjqY5GFgP/A0cDLwH5PcVlW/ur7dHGkfcMaK\n52cw2Jt4LgvXk8ewcG3h+hAWrmccY55xx61ziNXUJH4pyeeB3wY+A/xwVb0dOB/4qfXs4yrcDZyd\n5KwkxwJvAnbNuQ+StDRWsydxCvBTVc/9N6Tbu/jJ2XRruKp6Ksk7gT9jMAT2+qoRF6yRJE1lbJKo\nqn81ou3e9e3OeFV1C3DLvNcrScvIM64lScAGOeN6VSxcTx7DS4V7qfBDeKnwGceYZ9xx6xziiDrj\nWpI0XyYJSVKTSUKS1GThWpIEWLiejT4Xro7kwvW4M3gbhdRNxz05tgjbOrMYBmcXN5e9fzNcFPjz\nxuu6Cfhg4O2N9lFtlzMy9qZzvzu71/XkcYPic6s4PaowDaPfLwvX/Yo7bp1DeLhJktRkkpAkNZkk\nJElNJglJUpOjmyRJgKObZqPPoxsmHd20mmXGjYKaZtTKuZl8NM249tUs23r9X2f06KaHgMsDN7VG\nMI1o2zom9qmMv1fFIrfZNO/1uGVX+3l0dNP06xzCw02SpCaThCSpySQhSWqycC1JAixcz0afC1ez\nLly35jua6QvXsyyKj3p9Lw98odH+OPC6wGca7U8APxq4vdE+qu34MbFPAF4xom+vWMXrOlLfDwvX\n82HhWpK0ViYJSVKTSUKS1GThWpIEWLiejT4XrhZVuJ62fdGxW8Xhp4DzA59vtB8AXhu4s9E+qu2Y\nMbGPZvrC9ZH6fli4ng8L15KktTJJSJKaTBKSpCYL15IkwML1bPS5cGXhenhbq/gLo4vD4wrXjGmf\nZtnVFK7Hva6+vh8WrucXd9w6h/BwkySpySQhSWoySUiSmixcS5IAC9ez0efClYXr4W2TFq6nbV90\n7L6+Hxau5xd33DqH8HCTJKnJJCFJajJJSJKaLFxLkgAL17PR58KVhevhbRauD2fhejJ9/v5Pss4h\nPNwkSWoySUiSmkwSkqQmC9eSJOAIKVwnuRr458Cj3aT3VNUtXdu7gX8GPA38YlXdOjSIhevJYyxD\n4XrSe1xDvwvXG/X9sHA9H43Cde+SBFDAdVV13cqJSc4D3gScB5wO3JbknKo6uIA+StJS6GtNYlhK\nuwS4saoOVNUe4AFg21x7JUlLpq9J4qokX0hyfZKTumkvBvaumGcvgz0KSdKMLKRwnWQ3cNqQpvcC\nd/JsPeJ9wNaqujLJvwPurKqPdTE+BHyyqv7zIbEtXEvSBHpTuK6qH1/NfF0i+JPu6T7gjBXNL+mm\nHc7C9eQxLFxbuF7PZWcde63zTbvMLGLMM+64dQ7Ru8NNSbaueHopcE/3eBdweZJjk/wQcDZw17z7\nJ0nLpI+jm65N8koGo5y+CrwNoKruTbITuBd4CnhHbbSTPCSpZ3qXJKrq50a0XQNcM8fuSNJS693h\nJklSf5gkJElNXrtJkgT0aAjszDkEdvIYDoE9Mte9rLHXOt+0y8wixjzjjlvnEB5ukiQ1mSQkSU0m\nCUlSk0lCktRkkpAkNZkkJElNJglJUpMn00mSAE+mm40+n0zjyXTzjb3IdS9r7LXON+0ys4gxz7jj\n1jmEh5skSU0mCUlSk0lCktRkkpAkNZkkJElNJglJUpNJQpLUZJKQJDWZJCRJTV6WQ5IEeFmO2ejz\naflelmO+sRe57mWNvdb5pl1mFjHmGXfcOofwcJMkqckkIUlqMklIkppMEpKkJpOEJKnJJCFJajJJ\nSJKaTBKSpCaThCSpySQhSWoySUiSmkwSkqQmk4QkqclLhUuSAC8VPht9vlSwlwqfb+xFrntZY691\nvmmXmUWMecYdt84hPNwkSWoySUiSmkwSkqSmhSSJJD+d5ItJnk7yDw5pe3eSLye5L8mOFdPPT3JP\n1/Zv599rSVo+i9qTuAe4FPiLlROTnAe8CTgPuBj4QJJnqikfBK6sqrOBs5NcPMf+StJSWkiSqKr7\nqupLQ5ouAW6sqgNVtQd4ALggyVbgxKq6q5vvPwBvnE9vJWl59a0m8WJg74rne4HTh0zf102XJM3Q\nzM6TSLIbOG1I03uq6k9mtV5J0vqZWZKoqh+fYLF9wBkrnr+EwR7Evu7xyun7mlHef/Wzj1+zHbZt\nn6ArkrSB3XUHfO6OsbP14Yzrlaf57QI+nuQ6BoeTzgbuqqpK8u0kFwB3AT8L/F4z4i9cPbveStJG\nsG37c/+B/sBvDJ1tUUNgL03yIPBa4L8muQWgqu4FdgL3ArcA76hnLy71DuBDwJeBB6rqT+ffc0la\nLgvZk6iqm4GbG23XANcMmf554O/PuGuSpBX6NrpJktQjJgkYFHAkaSNY579nJglYVYVfko4I6/z3\nzCQhSWoySUiSmrx9qSQJGH770g2XJCRJ68fDTZKkJpOEJKlp6ZKEd8WbTpKrk+xN8lfdz0+saBu6\n/QRJLu62y5eT/Nqi+9NXSfYk+evus3VXN+2UJLuTfCnJrUlOWnQ/FynJh5PsT3LPimnNbTTt93Lp\nkgTeFW9aBVxXVa/qfm6B5vZbxs/XYZIcBfw+g+1yHvDmJC9bbK96q4Dt3WdrWzftXcDuqjoHuL17\nvsw+wuCztNLQbbQe38ul+xJ7V7x1cdgICIZvv21D5ltG2xhclHJPVR0AbmKwvTTcoZ+vNwA3dI9v\nYMm/f1X1KeCxQya3ttHU38ulSxIjeFe81bsqyReSXL9it7a1/TTYDg+ueO62aSvgtiR3J3lrN21L\nVe3vHu8Htiyma73W2kZTfy/7cD+Jdedd8aYzYvu9l8Ght9/snr8P+B3gykYox1cPuB1W73VV9VCS\nFwG7k9y3srG7t4zbc4RVbKM1bb8NmSQWele8DWC12y/Jh4Bnku6w7beht9MaHLptzuC5/92pU1UP\ndb8fTXIzg0Mj+5OcVlUPd4d/H1loJ/uptY2m/l4u++GmQ++Kd3mSY5P8EM/eFe9h4NtJLugK2T8L\nfGIBfe2F7gP4jEsZDASAxvabd/966m4GAx7OSnIsg0LirgX3qXeSPD/Jid3jzcAOBp+vXcAV3WxX\nsMTfvxFa22jq7+WG3JMYJcmlDG59eiqDu+L9VVX9RFXdm+SZu+I9xeF3xfso8Dzgk0t+V7xrk7yS\nwS7rV4G3weCugiO231KrqqeSvBP4M+Ao4Pqq+psFd6uPtgA3d4MKjwY+VlW3Jrkb2JnkSmAPcNni\nurh4SW4ELgJO7e7w+evAbzFkG63H99LLckiSmpb9cJMkaQSThCSpySQhSWoySUiSmkwSkqQmk4Qk\nqckkIUlqMklIkppMEtKMJXlNd9Xc45JsTvK/u+v8S73nGdfSHCR5H3A8g0u7PFhV1y64S9KqmCSk\nOUhyDIML/X0PuNDrWulI4eEmaT5OBTYDJzDYm5COCO5JSHOQZBfwceClwNaqumrBXZJWZekuFS7N\nW5KfA56sqpu6m9D/ZZLtVXXHgrsmjeWehCSpyZqEJKnJJCFJajJJSJKaTBKSpCaThCSpySQhSWoy\nSUiSmkwSkqSm/wf4j2LS/kQNqQAAAABJRU5ErkJggg==\n", + "text": [ + "" + ] + }, + { + "metadata": {}, + "output_type": "display_data", + "png": "iVBORw0KGgoAAAANSUhEUgAAAYkAAAEZCAYAAABiu9n+AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAGhBJREFUeJzt3X+sJeV52PHvY+7CYgyYH+GXwWFR2SooEbjQBWQqTBsw\nQcpigmJjCRso3jhFQKSSyDZ2ExpHNEQ1Sl3kRjEYE8eGoLrYmxYCi0tdm9ghuMGhJQQTs4HFsLhA\nwHFYvBee/nHm1oe7552Ze+b84p7vR1rtOfPOvPOed945z53zzI/ITCRJGuQN026AJGl2GSQkSUUG\nCUlSkUFCklRkkJAkFRkkJElFBglpBSLiwoj4Wt/7H0TEkdNrkTReBglpmYg4JSL+NCL+LiKejYiv\nR8QJg+bNzL0zc+uI139pRNwfETsi4sZlZUdGxKtVcFr699FRrl/qtzDtBkizJCL2Af4r8EHgVmAP\n4J8BL0+wGU8CHwfeCexZmGef9EpYTYBHEtJrrQcyM/8oe3Zk5pbMfHDQzNVf9UdVr/eMiE9ExNbq\nKORrEbG2KjupOjp5PiIeiIhTSw3IzNsy88vAszXtdN/VRDjQpNf6a+CViPhsRJwZEfutYNl/D7wN\nOBnYH/g14NWIeAu9o5PfzMz9gF8FvhgRBzbUFzVlfxsRT0TEZyLigBW0UVoRg4TUJzN/AJwCJPBp\n4JmI+HJEHFS3XES8AbgI+JXMfCozX83Mb2bmj4Dzgdsz80+qddwN3A+c1dScAdO+D5wAvBU4Htgb\n+HzrDyitkEFCWiYzH87MizLzCOCngcOA321Y7EBgLfA3A8p+EvjF6qem5yPieeDtwCENde5yJJGZ\nP8zM/1UFoWeAS4EzImKvhrqkoRgkpBqZ+dfATfSCRZ3/C+wA/tGAsseBz2Xmfn3/9s7M32la/Qqa\n6r6ssXBgSX0i4h9HxL+u8ghExBHAe4Fv1C2Xma8CnwGujYhDI2K3iDg5InYH/hD4+Yg4o5q+NiLe\nsbSOAW3YrUp4LwC7RcQeEbFbVbahauMbqlzEJ4F7qp/JpJEzSEiv9QPgRODPIuLv6QWHvwSuqMqT\n1/6F3//6V4EHgT+nd2bSvwPekJnbgLOBK4Fn6B1ZXEF5//s3wD8AH6KXz3gJWLoW4ijgDuDFal0v\n0Qti0liEp1pLkko8kpAkFRkkJElFBglJUpFBQpJUtOpu8BcRZuIlaQiZucsFnKsuSABctcL57wFO\n67Cula5vUvUOU0fbZZrm61L+eq17muue17pXOl/XZcZRx6jrHfb7rLQ+f26SJBUZJCRJRVMNEtVt\njrdHxIN90/aPiC0R8UhE3BURb+4r+0hEfCciHo6IM0bVjiNHVZEkTdmRI65v2kcSNwJnLpv2YWBL\nZq4HvlK9JyKOAd4DHFMt86nq9sydrRtFJZI0A0b9fTbVIJGZXwOeXzZ5I727blL9/67q9dnAzZm5\ns3qm8KPAhkm0U5Lm1bSPJAY5ODO3V6+3AwdXrw8DtvXNtw0YeBdNSdJozGKQ+P+qB73XXffgNRGS\nNEazeJ3E9og4JDOfjohD6d1aGeBJ4Ii++Q6vpu3inr7XR2LOQZKWewzY2mK+WQwSm4ELgGuq/7/U\nN/0LEXEtvZ+ZjgbuG1TBsBfGSdK8WMdr/4D+amG+qQaJiLgZOBU4MCKeAH4d+G3g1oi4mF6gezdA\nZj4UEbcCDwGLwCXpwzAkaaymGiQys/RErZ8tzH81cPX4WiRJ6jfTiWtJ0nQZJCRJRQYJSVKRQUKS\nVGSQkCQVGSQkSUWx2i418PGlkjScuXl86W9NcF0fo3dp+Kh9aAT1DlNH22Wa5utSPu6668bHx8ZY\nPu26Z3V7dKl7pfN1XWYcdZTqneT3GPTGyCD+3CRJKjJISJKKDBKSpCKDhCSpyCAhSSoySEiSigwS\nkqQig4QkqcggIUkqMkhIkooMEpKkIoOEJKnIu8BKkoA5ugvsJya4rivGtL4rgE92rOPyIepou0zT\nfF3Kx1133fZq2p5dyqdd96xujy51r3S+rssMqmNc+/8kv8eW1jmIPzdJkooMEpKkIoOEJKnIICFJ\nKjJISJKKDBKSpCKDhCSpyCAhSSoySEiSigwSkqQi790kSQLm6N5NXe/HshKXA58aQ72XjKDeS4BP\nr3CZTS2XaZqvS/m4667r16Z+71I+7bpndXt0qXul8y1fZhT72Lj2/0l+j0Hvu2wQf26SJBUZJCRJ\nRQYJSVKRQUKSVGSQkCQVGSQkSUVeJyFJArxOYiwuZ+XnZ7exCbixYx0XDVHHRcAftpjv/Ib5upRP\nu+66Pmvq07ryLsu2qXu1bo+243GYsT6KfWxc+7/XSUiSZt7MHklExFbgReAVYGdmboiI/YE/An4S\n2Aq8OzP/bmqNlKRVbpaPJBJ4R2a+LTM3VNM+DGzJzPXAV6r3kqQxmeUgAbA8ibIRuKl6fRPwrsk2\nR5LmyywHiQTujoj7I2JTNe3gzNxevd4OHDydpknSfJjZU2Aj4tDMfCoifgLYAlwGbM7M/frmeS4z\n91+23Gx+IEmaca+rU2Az86nq/+9HxG3ABmB7RBySmU9HxKHAM4OWHcete0suoftpdIO0PRW1zvnA\nLStc5jzgiy3mO7dhvi7l5wKba5bd2LHupmXr+uy8DuVdlm1T92rdHm3H4zBjfRT72Lj2/0l+j0Hv\nu2yQmfy5KSLeGBF7V6/3As4AHqQ3Vi+oZrsA+NJ0WihJ82FWjyQOBm6LCOi18fOZeVdE3A/cGhEX\nU50CO70mStLqN5NBIjMfA44bMP054Gcn3yJJmk8zm7gelolrSRrO6ypx3cU47qVSMop7LA1yEStP\nxC3XNgndrylJuWQjcHtN+Vk15WuA0+mdsjZIXVnX8tNr2gX17V6k97mbkril8i7LLjS0ra5sqbxL\nnw277M6GtrVpd9vxOMxYH8U+Nq79f5LfY9D7LhtkJhPXkqTZYJCQJBUZJCRJRSauJUmAieux2ET3\nqzYHOZ+VJ+KWO5f6pOAgTQnOJacD99SUn1ZTvgY4Bfh6ofwU4Js1dZ/UofykmvUurbtUvpP6z0VD\neZdl2/RZ0+fq0mdNyw7bZ236pO14HGasj2IfG9f+b+JakjTzDBKSpCKDhCSpyMS1JAkwcT0Wm+h+\n1eYg59HuStM6G2mX9Ot3OvUJ0CV1idA1wPHAtwrlC8CxwLcL5ccCD9Ws+5ia8gVgPfBIoXx9h7p3\nUt/uReo/d5c+WdPQtrqypfJSnyyOoO4ufTJsUrzfKQw31kexj41r/zdxLUmaeQYJSVKRQUKSVGTi\nWpIEzFHiehy37i25iO5XbQ4yzNXSy51F/dWsg5xGfSJxyUmUk5V7Up88XgOsAx4rlK8Dnq5Z9yHA\ns4WRu7AA++6AF9YOLq8rWyov1b1zsbfuUtsWgcOBbYXyurIF6uteswAHLJbbdsBi8+cqlS8uNtfd\ntD1K23InzScSlMYR9JL5bcfjMGN9FPvYuPb/SX6PQe+7bBB/bpIkFRkkJElFBglJUpFBQpJU5NlN\nkiTAs5vGYpxnN630NgPLnU67M0P61Z211O9Yymes7EP9mTp7rm0+A+lH+5bXvfsLkAcVCtdCPA75\n1sHF8TjkUeW647s15Tsgvgd5WKF8EeKZctvqylhoqHttfdtq271UXugTdjS3u2l7lLblSzvqx8Ih\nlMcRNJ/9tKTtWVD9TmI0+5hnN0mS5pZBQpJUZJCQJBWZuJYkAXOUuB7Hg8lLzqf7PekH2cjKbzOw\n3GmU7+Nfcjz1icQl6xk+Mb3PPi0SvE1J2OMKhXtB3Av59sKy90K+s6buOyE3FgpfrspLy++A+Crk\nqYW6a8pY21D3HhCby22LzS0+V6FP+CHEA+U+jQdabI/CtnzxxeaTFJpu+dF2PA4z1kexj41r/5/k\n9xj0vssG8ecmSVKRQUKSVGSQkCQVmbiWJAEmrsdi1hPXba5W7Xcs5WcD9FtHTWJ6L4hnIQ8oLLxv\ni6uHTyivO+6HPKtQeBDEZyEvLCz7Wcgra+q+GvK6QuECxC9D/l6h/DmIKyGvLtRdU8b+DXUvQlxa\nbltc2uJzXVgofAbi9nKfxu0ttkcpsf1C/ViIZ5ufg9F2PA4z1k1c/5iJa0nSihkkJElFBglJUpGJ\na0kSYOJ6LMaZuP56xzpOAR5a4TLHANtazHc48A97DS7bc9+G214fAPEg5M8MLo4HIf9Fed3xlZok\n7PoWyeM7a+p+J/z9jsEH2PftcSL/PL7Bf8+TB5bvyUucHA/wjcKly3VlL7Fnbd0bXv4z3rT21WLb\n3rT21cbPVUyaP9Ii2d+0PQrbkmfrx0J8rzyOAN74w/bjcZixPop9zMS1JGluGSQkSUUGCUlSkYlr\nSRKwShLXEXEm8LvAbsD1mXnN8nlMXPcMm7iuu3XzkkMoP/d4zb71z5lm/xa3pi5dUU11BfClhcLj\nID4AeX1h2Q9Afq+m7sPg9/N9A8t+afMfwNkBXy78HXL2vwWuqv4NUlP25d+orfv3N76fX4rPldsW\nn2v8XKU+4QGI68p9Gte12B6lW7c/1/zM8abnZ7cdjyauu1kVieuI2A24DjiT3jZ+b0T81HRbJUmr\nV2OQiIjLI2K/STSmhQ3Ao5m5NTN3ArcAZ0+5TZK0arU5kjgY+POIuDUizoyIXX6zmqC3AE/0vd9W\nTZMkjUGrxHVEvAE4A7gQOAG4FbghM/9mrK3btR3nAmdm5qbq/fnAiZl5Wd88Jq4laQhDJ64z89WI\neBrYDrwC7Af854i4OzN/bbTNrPUkcETf+yMYcEGmieseE9cDljVxvSsT10MzcQ1ExK9ExLeA3wHu\nBX46M/8VveeI/8LomtjK/cDREXFkROwOvIfxbCNJEu2OJPYHfiEz/7Z/YnV08fPjadZgmbkYEZcC\nd9I7BfaGzPyrSbZBkuZJY5DIzN+oKVvpEV5nmXkHcMek1ytJ88grriVJwCq54roNE9c93ip8wLLe\nKnxX3ip8aCauJUlzzSAhSSoySEiSikxcS5IAE9djMc7E9T0d6zgN+PYKlzkWeKzFfOuAF9YOLttn\nL4hnIQ8oLLwvxHchjxpcHN+FPKG87ri/5grgg1okYa+sqftqyOsKhQsQvwz5e4Xy51okzUvJ4/0b\n6l6EuLTctri0xee6sFD4THXVdKFP4/YW26OwLXmhfizEs+VxBLDvjvbjcZixPop9zMS1JGluGSQk\nSUUGCUlSkUFCklTk2U2SJMCzm8Zi1s9u+tYKlzkeeKTFfOsp3+d/z7W9s1KKZz/tA/EM5EGDy+OZ\nmrNlqM5+Kj2/YC+IeyHfXlj2Xsh31tR9J+TGQuHLVXlp+R0QX4U8tVB3TRlrG+reA2JzuW2xucXn\nKvQJP2zxfI+m7VHYli++WD8W9t1R/7yIQ2g/HocZ657d9GOe3SRJWjGDhCSpyCAhSSoycS1JAuYo\ncX3jBNd1EfDFMdR7LrClYx2nA99c4TIn0e72BsdSTijuQy/hOGxie98d8KN9y+ve/YVyopS1EI9D\nvnVwcTzeIglbKt/R8JyMxRYJ+VK7FxrqXtviViZNn6vQJ+xobnfT9ihty5d21I+FpsT0etqPx2HG\n+ij2sXHt/5P8HoPed9kg/twkSSoySEiSigwSkqQiE9eSJMDE9ViMM3F9e8c6zmLlV5SeRrsEYF2C\ne096CcdSQnINvfv/l54TsI7mq3CfLYzchYXmpHjT8wtKde9crE/CLgKHA9sK5XVlC9TXvWYBDlgs\nt+2AxebPVSpfXGyuu2l7lLblTurHQlNium1C+iSGG+uj2MdMXEuS5pZBQpJUZJCQJBWZuJYkAXOU\nuP70BNe1CbhlDPWeR/dbEG9k5VeUng58vcV8p1BOKK6hdxvm0q2bF+glJEsJy2OBh2rWfUxN+QLN\nidJh695JfbsXqf/cXfpkTUPb6sqWykt9sjiCurv0SV1i+iTaj8dhxvoo9rFx7f+T/B6D3nfZIP7c\nJEkqMkhIkooMEpKkIhPXkiTAxPVYbGI8z6I9n+5Xcg5z1fZZtEsAnk79Fa6n1ZSvoZdoLCUk65Li\n0EtmDlvelAita9dO6j8XDeVdlm3TZ02fq0ufDZtcbuqzNn3SdjwOM9ZHsY+Na/83cS1JmnkGCUlS\nkUFCklRk4lqSBJi4HotNjOeWvhfR/UrO81h5Yu5c2l2FupH6ROFZNeVr6CUaSwnJurKu5U0Jzrp2\nL9L73HX9U1feZdmFhrbVlS2Vd+mzYZfd2dC2Nu1uOx6HGeuj2MfGtf+buJYkzTyDhCSpyCAhSSoy\ncS1JAl4nieuIuAr4APD9atKVmXlHVfYR4F8CrwCXZ+Zdg+r41ATaueQSxpe46nol5/msPDHXNtl9\nbsN8XcqbkudNScqmupuWreuz8zqUd1m2Td2rdXu0HY/DjPVR7GPj2v8n+T0Gve+yQWYuSAAJXJuZ\n1/ZPjIhjgPfQu339W4C7I2J9Zr46hTZK0lyY1ZzELoc8wNnAzZm5MzO3Ao8CGybaKkmaM7MaJC6L\niG9HxA0R8eZq2mHAtr55ttE7opAkjclUEtcRsQU4ZEDRR+ndcHIpH/Fx4NDMvDgi/iPwzcz8fFXH\n9cDtmflfltVt4lqShjAzievMPL3NfFUg+OPq7ZPAEX3Fh1fTdvHJTq1bmcsZz5WRo7iS+6Ih6mib\nMD+/Yb4u5dOuu67Pmvq0rrzLsm3qXq3bo+14HGasj2IfG9f+P8nvMeh9lw0ycz83RcShfW/PAR6s\nXm8GzouI3SNiHXA0cN+k2ydJ82QWz266JiKOo3eW02PABwEy86GIuBV4iN5tdC7J1XaRhyTNmJkL\nEpn5/pqyq4GrJ9gcSZprM/dzkyRpdhgkJElF3rtJkgTM0Cmw4zbpU2DHcY+VS0ZQ7yWs/PS8TS2X\naZqvS/m4667r16Z+71I+7bpndXt0qXul8y1fZhT72Lj2f0+BlSTNPIOEJKnIICFJKjJISJKKDBKS\npCKDhCSpyCAhSSryYjpJEjBHF9N9YoLrumJM67uC7hfTXD5EHW2XaZqvS/m4667bXk3bs0v5tOue\n1e3Rpe6Vztd1mUF1jGv/n+T32NI6B/HnJklSkUFCklRkkJAkFRkkJElFBglJUpFBQpJUZJCQJBUZ\nJCRJRQYJSVKRt+WQJAFzdFuO35rguj4GXDOGej80gnqHqaPtMk3zdSkfd9114+NjYyyfdt2zuj26\n1L3S+bouM446SvVO8nsMemNkEH9ukiQVGSQkSUUGCUlSkUFCklRkkJAkFRkkJElFBglJUpFBQpJU\nZJCQJBUZJCRJRQYJSVKRQUKSVGSQkCQVeatwSRIwR7cKv2rC6xrH+kZR7zB1tF2mab4u5a/Xuqe5\n7nmte6XzdV1mHHVMst6mdQ7iz02SpCKDhCSpyCAhSSqaSpCIiF+MiP8TEa9ExD9ZVvaRiPhORDwc\nEWf0TT8+Ih6syv7D5FstSfNnWkcSDwLnAP+zf2JEHAO8BzgGOBP4VEQsZdv/E3BxZh4NHB0RZ06w\nvZI0l6YSJDLz4cx8ZEDR2cDNmbkzM7cCjwInRsShwN6ZeV813x8A75pMayVpfs1aTuIwYFvf+23A\nWwZMf7KaLkkao7FdJxERW4BDBhRdmZl/PK71SpJGZ2xBIjNPH2KxJ4Ej+t4fTu8I4snqdf/0J0uV\n3NP3+khg3RANkaTV7DFga4v5ZuGK6/7LwDcDX4iIa+n9nHQ0cF9mZkS8GBEnAvcB7wM+WarwtHG2\nVpJWgXW89g/orxbmm9YpsOdExBPAScB/i4g7ADLzIeBW4CHgDuCS/PHNpS4Brge+AzyamX8y+ZZL\n0nyZypFEZt4G3FYouxq4esD0bwE/M+amSZL6zNrZTZKkGWKQoJfAkaTVYNTfZwYJ2mX4Jen1YOuI\n6zNISJKKDBKSpCIfXypJAgY/vnTVBQlJ0uj4c5MkqcggIUkqmrsg4VPxuomIqyJiW0T8RfXv5/rK\nBvafICLOrPrlOxHxoWm3Z1ZFxNaI+MtqbN1XTds/IrZExCMRcVdEvHna7ZymiPhMRGyPiAf7phX7\nqOt+OXdBAp+K11UC12bm26p/d0Cx/+ZxfO0iInYDrqPXL8cA742In5puq2ZWAu+oxtaGatqHgS2Z\nuR74SvV+nt1Ibyz1G9hHo9gv524n9ql4I7HLGRAM7r8NA+abRxvo3ZRya2buBG6h118abPn42gjc\nVL2+iTnf/zLza8DzyyaX+qjzfjl3QaKGT8Vr77KI+HZE3NB3WFvqP/X64Ym+9/ZNWQJ3R8T9EbGp\nmnZwZm6vXm8HDp5O02ZaqY8675ez8DyJkfOpeN3U9N9H6f309pvV+48DnwAuLlTl+dU99kN7b8/M\npyLiJ4AtEfFwf2H1bBn7s0aLPlpR/63KIDHNp+KtBm37LyKuB5aC7qD+W9X9tALL++YIXvvXnSqZ\n+VT1//cj4jZ6P41sj4hDMvPp6uffZ6bayNlU6qPO++W8/9y0/Kl450XE7hGxjh8/Fe9p4MWIOLFK\nZL8P+NIU2joTqgG45Bx6JwJAof8m3b4ZdT+9Ex6OjIjd6SUSN0+5TTMnIt4YEXtXr/cCzqA3vjYD\nF1SzXcAc7381Sn3Ueb9clUcSdSLiHHqPPj2Q3lPx/iIzfy4zH4qIpafiLbLrU/E+C+wJ3D7nT8W7\nJiKOo3fI+hjwQeg9VbCm/+ZaZi5GxKXAncBuwA2Z+VdTbtYsOhi4rTqpcAH4fGbeFRH3A7dGxMX0\nbnL67uk1cfoi4mbgVODA6gmfvw78NgP6aBT7pbflkCQVzfvPTZKkGgYJSVKRQUKSVGSQkCQVGSQk\nSUUGCUlSkUFCklRkkJAkFRkkpDGLiH9a3TV3j4jYKyL+d3Wff2nmecW1NAER8XFgLb1buzyRmddM\nuUlSKwYJaQIiYg29G/29BJzsfa30euHPTdJkHAjsBbyJ3tGE9LrgkYQ0ARGxGfgCcBRwaGZeNuUm\nSa3M3a3CpUmLiPcDL2fmLdVD6P80It6Rmf9jyk2TGnkkIUkqMichSSoySEiSigwSkqQig4Qkqcgg\nIUkqMkhIkooMEpKkIoOEJKno/wGw2I6lWJlX/wAAAABJRU5ErkJggg==\n", + "text": [ + "" + ] + }, + { + "metadata": {}, + "output_type": "display_data", + "png": "iVBORw0KGgoAAAANSUhEUgAAAYkAAAEZCAYAAABiu9n+AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAFEtJREFUeJzt3X+spuVd5/H3h+kOpSymIHagLTtjIzRONGm33cFGt8wf\nyo4mlKJxahOVKjaNTdHEan/Q4I6SsNIsxBXj/rGlLbu21IlKHbQoQyOlVSviT9yRUkxnZAgMXcW0\nu6X8GL77x3NTDjP3dZ7nnOf3Oe9XMpnnua/r+l4Xd87hO/f9vX+kqpAkqc8p816AJGlxmSQkSU0m\nCUlSk0lCktRkkpAkNZkkJElNJglpDZK8LclnV3z/apId81uRNF0mCekESb4nyZ8m+dck/5zkc0le\n39e3qs6oqsMTnv9dSe5N8vUkHzmhbUeSZ7vk9NyfD0xyfmmlF817AdIiSfJNwO8D7wD2A6cC/xF4\ncobLeBi4BvhPwGmNPt9U3gmrGfBIQnqhC4Cqqt+qga9X1cGquq+vc/ev+ld1n09Lcn2Sw91RyGeT\nvLhr+67u6OTxJH+T5KLWAqrq1qr6PeCfV1mnv7uaCX/QpBf6AnA8yUeT7Ely5hrG/lfgtcAbgLOA\nXwCeTfIKBkcnv1xVZwI/D/xOkrOHxMsqbUeSPJTkw0m+eQ1rlNbEJCGtUFVfBb4HKOB/AI8l+b0k\nL1ttXJJTgJ8AfraqHqmqZ6vq81X1FPCjwKeq6g+7Oe4E7gV+YNhyerZ9GXg98O+A1wFnAB8b+T9Q\nWiOThHSCqrq/qn6iqs4DvgN4OfCrQ4adDbwY+Meetu3AD3enmh5P8jjw3cA5Q2KedCRRVf+vqv6q\nS0KPAe8CLk5y+pBY0rqYJKRVVNUXgJsZJIvV/B/g68C39bT9E/C/qurMFX/OqKoPDpt+DUv1d1lT\n4Q+WtEKSVyf5ua6OQJLzgLcCf7bauKp6FvgwcEOSc5NsSfKGJFuB3wQuSXJxt/3FSXY/N0fPGrZ0\nBe8XAVuSnJpkS9e2q1vjKV0t4teAP+5Ok0kTZ5KQXuirwIXAnyf5vwySw98B7+7aixf+C3/l558H\n7gP+gsGVSf8FOKWqjgKXAlcBjzE4sng37d+/q4GvAe9lUM94AnjuXohXAbcDX+nmeoJBEpOmIl5q\nLUlq8UhCktRkkpAkNZkkJElNJglJUtOGe8BfEivxkrQOVXXSDZwbLkkAXFVXr6n/3fs+wxv3NZ+3\ntqprc82a55tV3PXEGHXMsH7jtC9r7FHmvrquao69Jteuu32cscsce639xh0zjRiTjvuZfXdz0b43\nrmvOPp5ukiQ1mSQkSU1zTRLdY46PJblvxbazkhxM8kCSO5K8dEXb+5N8Mcn9SS6e1Dq2794+qVCS\nNFeT/v/ZvI8kPgLsOWHb+4CDVXUB8OnuO0l2Am8BdnZjfqN7PPPYtu/eMYkwkjR3OzZSkqiqzwKP\nn7D5TQyeukn395u7z5cCt1TV0907hR8Eds1inZK0Wc37SKLPtqo61n0+BmzrPr8cOLqi31Gg9yma\nkqTJWMQk8Q3di95Xu+/BeyIkaYoW8T6JY0nOqapHk5zL4NHKAA8D563o98pu20nu3veZb3zevnu7\nNQdJOsHhu45w5K4jQ/stYpI4AFwOXNf9/ckV2z+e5AYGp5nOB+7pC7DeG+MkabPYsXv7C4rcd//S\n53r7zTVJJLkFuAg4O8lDwC8CvwLsT3IFcBjYC1BVh5LsBw4BzwDvLF+GIUlTNdckUVWtN2p9b6P/\ntUD/veOSpIlb6MK1JGm+TBKSpCaThCSpySQhSWoySUiSmkwSkqSmbLRbDXx9qSStj68vnQJfXzr5\n9mWNPcrci/oa0GWNvdZ+446ZRoxZxh02Zx9PN0mSmkwSkqQmk4QkqckkIUlqMklIkppMEpKkJpOE\nJKnJJCFJajJJSJKaTBKSpCaThCSpySQhSWryKbCSJMCnwE6FT4GdfPuyxh5l7kV9muqyxl5rv3HH\nTCPGLOMOm7OPp5skSU0mCUlSk0lCktRkkpAkNZkkJElNJglJUpNJQpLUZJKQJDWZJCRJTSYJSVKT\nz26SJAE+u2kqfHbT5NuXNfYocy/qM5CWNfZa+407ZhoxZhl32Jx9PN0kSWoySUiSmkwSkqQmk4Qk\nqckkIUlqMklIkpq8T0KSBHifxFR4n8Tk25c19ihzL+r9Bssae639xh0zjRizjDtszj6ebpIkNS3s\nkUSSw8BXgOPA01W1K8lZwG8B24HDwN6q+te5LVKSNrhFPpIoYHdVvbaqdnXb3gccrKoLgE933yVJ\nU7LISQLgxCLKm4Cbu883A2+e7XIkaXNZ5CRRwJ1J7k3y9m7btqo61n0+Bmybz9IkaXNY2Etgk5xb\nVY8k+RbgIHAlcKCqzlzR51+q6qwTxi3mf5AkLbilugS2qh7p/v5ykluBXcCxJOdU1aNJzgUe6xvr\nJbDrj+ElsF4Cu0yx19pv3DHTiDHLuMPm7LOQp5uSvCTJGd3n04GLgfuAA8DlXbfLgU/OZ4WStDks\n6pHENuDWJDBY48eq6o4k9wL7k1xBdwns/JYoSRvfQiaJqvoS8Jqe7f8CfO/sVyRJm9PCFq7Xy8K1\nJK3PUhWux2Hhev0xLFxbuF6m2GvtN+6YacSYZdxhc/ZZyMK1JGkxmCQkSU0mCUlSk4VrSRJg4Xoq\nLFxPvn1ZY48y96IWgJc19lr7jTtmGjFmGXfYnH083SRJajJJSJKaTBKSpCYL15IkwML1VFi4nnz7\nssYeZe5FLQAva+y19ht3zDRizDLusDn7eLpJktRkkpAkNZkkJElNFq4lSYCF66mwcD359mWNPcrc\ni1oAXtbYa+037phpxJhl3GFz9vF0kySpySQhSWoySUiSmkwSkqQmr26SJAFe3TQVXt00+fZljT3K\n3It6ldCyxl5rv3HHTCPGLOMOm7OPp5skSU0mCUlSk0lCktRk4VqSBFi4ngoL15Nvn3bs99a+5tjr\nsm/d7Vs4vu7C9RaOsy/Xsa/e25x7tfZRxlq4Hn/MNGLMMu6wOft4ukmS1GSSkCQ1mSQkSU0WriVJ\ngIXrqbBwPfn2eReu313X9LY9xancmPdwZX2wt30rT3J9rm6Ovz5Xc029u7ftVJ7iPbmRD9aVve1P\nspWrc31z/NW53sK1heux5+zj6SZJUpNJQpLUZJKQJDVZuJYkARaup8LC9eTb5124bhWmj/Ey9udt\n7K2P9rY/wUu4LXu5pPb3tt+WveyvS3rbvsZpvC37+Wjt7W1/jG2rFrbfkxstXFu4HnvOPp5ukiQ1\nmSQkSU0mCUlSk4VrSRKwQQrXSfYAvwpsAT5UVded2MfC9fpjbIbC9bCxrTumH+DVQwrTh4B93Z8+\n+5otO+sS9ua2ZmH7C1ww9I7rRS0uW7henLjD5uyzVKebkmwBfh3YA+wE3prk2+e7KknauIYmiSQ/\nk+TMWSxmBLuAB6vqcFU9DXwCuHTOa5KkDWuUI4ltwF8k2Z9kT5KTzlnN0CuAh1Z8P9ptkyRNwUiF\n6ySnABcDbwNeD+wHbqqqf5zq6k5exw8Be6rq7d33HwUurHr+DiML15K0PusuXFfVs0keBY4Bx4Ez\ngd9OcmdV/cJkl7mqh4HzVnw/j8HRxAvMsuCzyIWraRbzFrmYOW7sVnH41TywanH5UG4bUram2XpJ\n7Vy1KH4BXxj6GPKNeiHBJC6kmNSYacSYZdxhc/YZpSbxs0n+Evgg8CfAd1TVTwOvA35wkoscwb3A\n+Ul2JNkKvAU4MOM1SNKmMcqRxFnAD1bVkZUbu6OL/n9STUlVPZPkXcAfMbgE9qaq+odZrkGSNpOh\nSaKq/vMqbYcmu5zhqup24PZZzytJm5F3XEuSgA1yx/UoLFyvP8ZmKFwPe6R263HcL+PYqo/zfglP\nrFrY3pvbmoXp0/jaqo8h38Zjq75f+8a8Z+gj0C1cjz9mGjFmGXfYnH2W6o5rSdJsmSQkSU0mCUlS\nk4VrSRJg4XoqLFxPvn3ehevWHden8tSq75l+kq1DH+fdumP6KU5dtTC9lSeH3nFt4drC9bhz9vF0\nkySpySQhSWoySUiSmkwSkqQmr26SJAFe3TQVXt00+fZ5X9203vbjbBk6d+uKleNs4brsG3qFUqt9\nlLFe3TT+mGnEmGXcYXP28XSTJKnJJCFJajJJSJKaLFxLkgAL11Nh4Xry7csae5S5F7UAvKyx19pv\n3DHTiDHLuMPm7OPpJklSk0lCktRkkpAkNVm4liQBFq6nwsL15NuXNfYocy9qAXhZY6+137hjphFj\nlnGHzdnH002SpCaThCSpySQhSWqycC1JAixcT4WF68m3L2vsUeZe1ALwssZea79xx0wjxizjDpuz\nj6ebJElNJglJUpNJQpLUZOFakgRYuJ4KC9eTb1/W2KPMvagF4GWNvdZ+446ZRoxZxh02Zx9PN0mS\nmkwSkqQmk4QkqcnCtSQJsHA9FRauJ9++rLFHmXtRC8DLGnut/cYdM40Ys4w7bM4+nm6SJDWZJCRJ\nTSYJSVKThWtJErAkhesk+4CfAr7cbbqqqm7v2t4P/CRwHPiZqrqjL4aF6/XHsHBt4XqZYq+137hj\nphFjlnGHzdln4ZIEUMANVXXDyo1JdgJvAXYCrwDuTHJBVT07hzVK0qawqDWJkw55gEuBW6rq6ao6\nDDwI7JrpqiRpk1nUJHFlkr9NclOSl3bbXg4cXdHnKIMjCknSlMylcJ3kIHBOT9MHgM/zfD3iGuDc\nqroiyY3A56vqY12MDwGfqqrfPSG2hWtJWoeFKVxX1feN0q9LBLd1Xx8GzlvR/Mpu20ksXK8/hoVr\nC9fLFHut/cYdM40Ys4w7bM4+C3e6Kcm5K75eBtzXfT4A/EiSrUm+FTgfuGfW65OkzWQRr266Lslr\nGFzl9CXgHQBVdSjJfuAQ8AzwztpoN3lI0oJZuCRRVT++Stu1wLUzXI4kbWoLd7pJkrQ4TBKSpCaf\n3SRJAhboEthp8xLY9cfwElgvgV2m2GvtN+6YacSYZdxhc/bxdJMkqckkIUlqMklIkppMEpKkJpOE\nJKnJJCFJajJJSJKavJlOkgR4M91UeDPd5NuXNfYocy/qTWnLGnut/cYdM40Ys4w7bM4+nm6SJDWZ\nJCRJTSYJSVKTSUKS1GSSkCQ1mSQkSU0mCUlSk0lCktRkkpAkNflYDkkS4GM5psLHcky+fVljjzL3\noj7eYlljr7XfuGOmEWOWcYfN2cfTTZKkJpOEJKnJJCFJajJJSJKaTBKSpCaThCSpySQhSWoySUiS\nmkwSkqQmk4QkqckkIUlqMklIkppMEpKkJh8VLkkCfFT4VPio8Mm3L2vsUeZe1EduL2vstfYbd8w0\nYswy7rA5+3i6SZLUZJKQJDWZJCRJTXNJEkl+OMn/TnI8yb8/oe39Sb6Y5P4kF6/Y/rok93Vt/232\nq5akzWdeRxL3AZcBd6/cmGQn8BZgJ7AH+I0kz1Xb/ztwRVWdD5yfZM8M1ytJm9JckkRV3V9VD/Q0\nXQrcUlVPV9Vh4EHgwiTnAmdU1T1dv/8JvHk2q5WkzWvRahIvB46u+H4UeEXP9oe77ZKkKZrafRJJ\nDgLn9DRdVVW3TWteSdLkTC1JVNX3rWPYw8B5K76/ksERxMPd55XbH24F+cy+50sd23dvZ8fu7etY\niiRtXEfuOsyRu44M7bcId1yvvA38APDxJDcwOJ10PnBPVVWSryS5ELgH+DHg11oBL9r3xmmuV5KW\n3vbdO9i+e8c3vn/ul+7u7TevS2AvS/IQ8F3AHyS5HaCqDgH7gUPA7cA76/mHS70T+BDwReDBqvrD\n2a9ckjaXuRxJVNWtwK2NtmuBa3u2/yXwnVNemiRphUW7ukmStEBMEsDhEYo3krQMjtx1eKLxTBIw\nUoVfkpbBpP9/ZpKQJDWZJCRJTb6+VJIE9L++dMMlCUnS5Hi6SZLUZJKQJDVtuiThW/HGk2RfkqNJ\n/rr78/0r2nr3nyDJnm6/fDHJe+e9nkWV5HCSv+t+tu7ptp2V5GCSB5LckeSl817nPCX5cJJjSe5b\nsa25j8b9vdx0SQLfijeuAm6oqtd2f26H5v7bjD9fJ0myBfh1BvtlJ/DWJN8+31UtrAJ2dz9bu7pt\n7wMOVtUFwKe775vZRxj8LK3Uu48m8Xu56X6JfSveRJx0BQT9+29XT7/NaBeDh1IerqqngU8w2F/q\nd+LP15uAm7vPN7PJf/+q6rPA4ydsbu2jsX8vN12SWIVvxRvdlUn+NslNKw5rW/tPg/3w0Irv7pu2\nAu5Mcm+St3fbtlXVse7zMWDbfJa20Fr7aOzfy0V4n8TE+Va88ayy/z7A4NTbL3ffrwGuB65ohPL6\n6gH3w+i+u6oeSfItwMEk969s7N4t4/5cxQj7aE37b0MmiXm+FW8jGHX/JfkQ8FzS7dt/G3o/rcGJ\n++Y8XvivO3Wq6pHu7y8nuZXBqZFjSc6pqke707+PzXWRi6m1j8b+vdzsp5tOfCvejyTZmuRbef6t\neI8CX0lyYVfI/jHgk3NY60LofgCfcxmDCwGgsf9mvb4FdS+DCx52JNnKoJB4YM5rWjhJXpLkjO7z\n6cDFDH6+DgCXd90uZxP//q2itY/G/r3ckEcSq0lyGYNXn57N4K14f11V319Vh5I891a8Zzj5rXgf\nBU4DPrXJ34p3XZLXMDhk/RLwDhi8VXCV/bepVdUzSd4F/BGwBbipqv5hzstaRNuAW7uLCl8EfKyq\n7khyL7A/yRXAYWDv/JY4f0luAS4Czu7e8PmLwK/Qs48m8XvpYzkkSU2b/XSTJGkVJglJUpNJQpLU\nZJKQJDWZJCRJTSYJSVKTSUKS1GSSkCQ1mSSkKUvyH7qn5p6a5PQkf989519aeN5xLc1AkmuAFzN4\ntMtDVXXdnJckjcQkIc1Akn/D4EF/TwBv8LlWWhaebpJm42zgdODfMjiakJaCRxLSDCQ5AHwceBVw\nblVdOeclSSPZdI8Kl2YtyY8DT1bVJ7qX0P9pkt1VddeclyYN5ZGEJKnJmoQkqckkIUlqMklIkppM\nEpKkJpOEJKnJJCFJajJJSJKaTBKSpKb/D9nKmwOrstnaAAAAAElFTkSuQmCC\n", + "text": [ + "" + ] + }, + { + "metadata": {}, + "output_type": "display_data", + "png": "iVBORw0KGgoAAAANSUhEUgAAAYkAAAEZCAYAAABiu9n+AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3X+QXfV53/H3x6tqjbCCkBWBjNU4/HDHDJ1i48rgqIOG\nJAI6YwtHEzk4IibdWh4xUrLTTYSFR/E6mhERUzFqxKxntpYtV0Y4mqr8ai0jbBfXKlAFHBxccAzE\npKCRJWN+WLbkFVqe/nHv2pflPOfcvbt396L9vGYEu+c5v/Zwl0f3Pt/v81VEYGZmVuQtU30DZmbW\nuZwkzMws5SRhZmYpJwkzM0s5SZiZWcpJwszMUk4SZmMg6XpJ3274/qikd03dHZm1l5OE2SiSFkt6\nUNLLkn4iab+k9xftGxGzI+LZCb7+GkmPSPqFpC+Oir1L0mv15DTy59MTeX2zRjOm+gbMOomkXwP+\nO/BJYDfQDfwbYGgSb+MgsBG4Ejgt2efXwjNhbRL4nYTZ670biIj4m6j5RUTcHxGPF+1c/1v9ufWv\nT5O0RdKz9Xch35b01nrs0vq7k5ckPSbp8uwGIuLOiLgb+EnJffp31yaFX2hmr/cPwLCkHZKuknTm\nGI79j8B7gcuAucCfA69JOofau5O/jIgzgT8D9kiaV3E+lcT+SdJzkr4g6e1juEezMXGSMGsQEUeB\nxUAA/xk4IuluSfPLjpP0FuCPgT+NiEMR8VpEPBwRJ4CVwFcj4mv1a3wdeAT4t1W3U7Dtx8D7gX8O\nXALMBm5v+gc0GyMnCbNRIuL7EfHHEbEQuAh4B7C14rB5wFuBZwpivwH8fv2jppckvQT8FnB2xTnf\n8E4iIn4eEd+pJ6EjwBpgqaTTK85l1hInCbMSEfEPwJeoJYsyLwC/AM4viP0/YGdEnNnwZ3ZE3FJ1\n+THcqn+XrS38wjJrIOlfSPoP9ToCkhYC1wIPlR0XEa8BXwBulbRAUpekyyTNBL4MfEjS0vr2t0pa\nMnKNgnvoqhe8ZwBdkrolddVji+r3+JZ6LeKvgf9Z/5jMbMI5SZi93lHgA8D/kfQzasnh74G+ejx4\n/d/wG7/+M+Bx4G+pjUy6GXhLRDwPLANuAo5Qe2fRR/77twE4BtxIrZ5xHBiZC3EusBf4af1ax6kl\nMbO2kIdam5lZxu8kzMws5SRhZmYpJwkzM0s5SZiZWeqUa/AnyZV4M7MWRMQbJnCeckkCYEvcMKb9\n7+s/wJX9i1q6Vp8G2BqrWjq2TK8Gx33eVs7R7DG9Gix9zn0aaDnepwFuibXpseu0LY0P08V6beXm\n6C2Mr9dWNkZfYQxgg7ak8WFm0K/N9MeN6bU3ahMb4qbC+EZt4qbYkB67Wf3cGP2F8ZN0sUUb6IuN\nhfGy2Ei8N25Or71N61ibzO0ri43Eb4gtaXxAfWm8LDYSXxVVk91hUL1N7TfeY9pxjuy8Zc+lzIH+\n+1jUf+WYjxtQ8eveHzeZmVnKScLMzFJTmiTqbY4PS3q8YdtcSfdL+oGkfZLmNMTWS3pK0vclLZ2o\n+zhvSWF3BDOzN51zlpw3oeeb6ncSXwSuGrXtU8D9EfFu4Bv175F0IfBR4ML6MQP19szjdr6ThJmd\nIs5ZUtRjsnVTmiQi4tvAS6M2f5ha103q/76m/vUy4I6IeLW+pvDTQGvVZjMza8pUv5MoclZEHK5/\nfRg4q/71O4DnG/Z7HvBbADOzNurEJPFL9YXey+Y9eE6EmVkbdeI8icOSzo6IH0laQK21MsBBYGHD\nfu+sb3uD+/oP/PLr85ac45qDmdkoBx94moMPFC2k+HqdmCTuAT4ObK7/+66G7bsk3UrtY6YLgANF\nJ2h1YpyZ2XRxzpLzX1fkfuSz+wr3m9IkIekO4HJgnqTngL8A/grYLakHeBZYARART0jaDTwBnARu\nCC+GYWbWVlOaJCIiW1Hrd5L9NwGb2ndHZmbWqKML12ZmNrWcJMzMLOUkYWZmKScJMzNLOUmYmVnK\nScLMzFI61aYaePlSM7PWTJvlS9uxnGimV4PcFj0Tft412j7u87ZyjmaPqdqvLH6SrtJlUquWUC1b\nOnWYrtLlTddpW7q0KZQvb9rpy5dmy5MCbNX6dAnSYboqlxgtW6KzbAnPYWawXWvoidsK42WxZuJj\n3W+8x7TjHNl527EsaplBFf9e+OMmMzNLOUmYmVnKScLMzFJOEmZmlnKSMDOzlJOEmZmlnCTMzCzl\nJGFmZiknCTMzSzlJmJlZyknCzMxSThJmZpZyF1gzMwPcBbYt2tkF9nNx/bjOsVo7xnyO1drBYFxX\nud8q7Uz3G6ar8tpl8WaOHU+H2ayDLECfBtIusSfpYoO2pF1iT9DdchfYIWaWdnkdZgZbtT7t9LpV\n69MurlDeybWZTq3Xx+fSc+/Q6pbjO7Sa62IwPXanVpXGG/cru4ex3tdYzuEusGZmNm05SZiZWcpJ\nwszMUk4SZmaWcpIwM7OUk4SZmaWcJMzMLOUkYWZmKScJMzNLOUmYmVnKvZvMzAxw76a26PTeTc30\nYWq0SjvZHh+r3K9Hu0r3G0+8mWNb7Ru1WjtKXx9lvZ1OMJP12pr2dhqim35tpj9uLIz3azM3Rn96\nbFlvphPMZEB9aX+mst5MUOvLk/UpGqartEfSTq3iY7E9Pfcu9bQcH8+xo/drpsdTo1b6PY3m3k1m\nZjatOUmYmVnKScLMzFJOEmZmlnKSMDOzlJOEmZmlPE/CzMwAz5Noi3bOkxjrHIfRmp3z0KhHu9gR\nKyr3u167S/cbT3w8xw4zo3IORtX62dnr5wTdrNM2bom1hfFjzGp5nsQxZrFN61gbtxTGh5jJoHrT\nsfNl8yCgNp4/m29wki5263pWxI7CeFlsvPHxnrtxv2bmUzRqZW7FaDu1yvMkzMxs+urYdxKSngV+\nCgwDr0bEIklzgb8BfgN4FlgRES9P2U2amZ3iOvmdRABLIuK9EbGovu1TwP0R8W7gG/XvzcysTTo5\nSQCMLqJ8GPhS/esvAddM7u2YmU0vnZwkAvi6pEckfaK+7ayIOFz/+jBw1tTcmpnZ9NCxQ2AlLYiI\nQ5J+HbgfWAvcExFnNuzzYkTMHXVcZ/5AZmYd7k01BDYiDtX//WNJdwKLgMOSzo6IH0laABwpOjZr\n9dwOfRqYtkNgvxzL0/hK7UnjXQxzre7ijij+tLAsVhUfpqv02iu1p7INeTZE9hiz6NNA+vo6yuyW\nh8AeZXZpK/BjzGKHVqfDXMuGuEJtuOfy+HJhbJgu7tK1XBN3FMbLYlXxLobZo5XptctiI/HpOgQ2\ney20y4D6Crd35MdNkmZJml3/+nRgKfA4cA/w8fpuHwfumpo7NDObHjr1ncRZwJ2SoHaPt0fEPkmP\nALsl9VAfAjt1t2hmdurryCQRET8ELi7Y/iLwO5N/R2Zm01PHFq5b5cK1mVlr3lSF6/E4VQrXE7HG\n9VQUrqsK010Ms0L3sjs+VBgvi403vkL3Vhbcs2d2lLfRq8G0t9NRZrNBW9gYxQXADdpCX2xMjy3r\nzfQzZleuFV1VAP5Q7E7j92pFGi+LVcW7GK4sik9l4bqT17h24drMzDqek4SZmaWcJMzMLOXCtZmZ\nAS5ct0WfBtqyyFGvBju+cN1qYbqLkyzXXvbE1YXx5drL3bE0vfYy7Ws5vkz70uuOXDsrbB9lNqu1\nI/3v8gJvb7lw/RPmlc6oPsrsypnLV8eewhjAXi1nadydxvdpWRovi43Es2vPYLiyKF41m7vTC9ft\nWBxoUL0uXJuZWedzkjAzs5SThJmZpVy4NjMzYBoVrm+JtZN2rXXa5sL1KFWF626GKovLe2NJeu2r\n9QD7YnEaX6r9aXyp9leeOytsv8yZ9GhX+kyPMJ/12srN0VsYX6+t9MbN6bFlM6pfYg57tTwtEO/V\ncpbE3sIYwAO6msWxL43v19I0vl9LK8+dFba7GXLhugWD6mVt3DLh5y2zTesKt/vjJjMzSzlJmJlZ\nyknCzMxSLlybmRngwnVbrNO2tszwnogW5K2sk71KOzuicF1VmP5WLEqvvViPsj8uKYyXxUbi34zL\nCmMvMK/053qOhaVrYPdpIJ1F+xwLSwu8LzCPh3QFl8U3C+MP6Qouif2FMYBHtTiNdzHMAV3OovhW\nYfyALq8sendC4Xqs61VPxPrU7WrpPaA+F67NzKzzOUmYmVnKScLMzFJOEmZmlvLoJjMzA6bR6Kas\nLUI7rNfWtoymmoh2H6209mi2lUePdqXrLszkRNtHN2UjlGZxjPfpSb4T7ymMv09P8r04Lz33RXom\njR9hPlfooXT00xNcyBptT0elrdH2dDTNE1xYOnrpCPN5RhdxXnyvMF4WG4m/J75TGJvFsdLRT49q\ncVtHN2VrZEBtnYxm2m200mJjIlpqtKt9xjatS1u4tMtWrS/c7o+bzMws5SRhZmYpJwkzM0u5cG1m\nZsA0KlxnC9G3wwZtaUuhfL22jrvdR58Gxlz8brbYXVbg7maIldqTFrZP4xjLtTddt2G59lau+ZAV\nrufwcmnx+SI9w6E4Iz33Ar3Ci3FaYexl5nCuDvGPsaAw/iAfLP25V2pPWqR9kA9ySOeyIP6x+No/\nncPxM+Zy2isvFsaPnzGXM4YOFcYAXulekBa25/ByZeG6aj2JbJ2LWRxnj1amP3dVYbrZgnQrRehB\n9Y67pcaA+tpSYN6q9fTFxgk/b5kt2lC43R83mZlZyknCzMxSThJmZpZy4drMzAAXrttig7a05XoT\nURBvpfjdbLG7rMA9kxP0aFfLhe2V2pMWtaFW2M5mPc/nSGXhOitMA8zVceKpJPhz0MUQjxWHH/5X\nF3OZHuOhuLgwfpke4+J4qDD22Dcvhd8WfCP5O87pwKWCh5P4pUqL2lArbGeF6/kcqVyrIitMA+zV\n8rQw3c0Qu9STFqerCtPNFqRbKUJPRNG5XQXmLdrgwrWZmXU+JwkzM0s5SZiZWcqFazMzA6ZR4bo/\nbpy8a2lzW67Xr83jLoi3UvxutvX5Om1Li+JdnKRXg2kBvJsTrNaOtPBd1a68R7vSNuTzeKG0nfcV\neiidMQ1wrg6lhWnuAG2G9D/3FaArIe4rDutK4Pbk7zBfA3YKrkvi11Ne2P5tpbO1AQ7p3LQwPY8X\nKtt5tzoreiZDDKo3LT5XFZwH1NdUK+5WWmtPRNF5izZwY/SP6xxFNqu/LeetumYRf9xkZmYpJwkz\nM0s5SZiZWcqFazMzA06RwrWkq4CtQBfw+YjYPHqfDXHTpN3PRm1qy/Um4rwbtWnMRfVmC/Fl+3Vx\nsnImelm8quBeNpO8m6HKdaazmd5Qm+2dzZi+dN9jpYXpz14J/dT+FOlv+Ocb3P4Z+EOlhe2LP/Yw\nj+myfMa2LqtcKzpbX7ubIQbUlxaQq2YmVxWAy2YPVxV+my3gtlLo3ax+boriWcbN2qSN4z7HZJ63\n6ppF3lQfN0nqAm4DrgIuBK6VVLzivZmZjVtlkpD0J5LOnIybacIi4OmIeDYiXgW+Aiyb4nsyMztl\nNfNO4izgbyXtlnSVpDd8ZjWJzgGea/j++fo2MzNrg6YK15LeAiylNqXn/cBuYHtEPNPWu3vjfSwH\nroqIT9S/Xwl8IOJXs79cuDYza03LheuIeE3Sj4DDwDBwJvBfJX09Iv58Ym+z1EFgYcP3C6m9m3id\nySz4dHLhapM2tlTMG2+h8CRdla2Oq4qZVYXSrMg6RDfbtSYt0m7XmsoCb1oc3nVpaXGZP/wsVaXr\nLPKZ+8pnaz+8tLoNeVVBPivmD9FNnwbSwQB9GqgcSDCeQQplAyUmYiBF2TETMTikUweutHLNIs3U\nJP5U0qPALcD/Bi6KiNXAJcDvTeRNNuER4AJJ75I0E/gocM8k34OZ2bTRzDuJucDvRcQ/NW6sv7so\nbqDTJhFxUtIa4D5qQ2C3R8STk3kPZmbTSWWSiIjPlMSemNjbqRYRe4G9k31dM7PpyDOuzcwMOEVm\nXDdjMlvstqul72b1T0gb41baJzfbmjkrHp+kq7Q99Am62aHVaXvpHVpd2Zo6a2v9AvMq12uuaqmd\ntuPeQXk776uoKGwrLUzzzYo25NeWr6+ti6lsgZ61T3+BeazQvWn79RW6t7J1e9b2fYiZpW3jezVY\nug57nwaabl3fSlv8iWjH366lAiZzyYORaxZ5U824NjOzyeUkYWZmKScJMzNLuXBtZmbANCpcj7fg\nOxZVM4vHc96xFp1HK5uZnBlQX1pwbjSo3rTwfIKZ7FJPWnweops9WpnOfN6jlVwde9Jr79XytDB9\nhPk8o4s4L75XGH9GF3HaKy+m5z5+xlx4OPl7xs8pXWf64iuq23mnbci/+1hpYZrTQRdAPFUc1gXw\nYpyWHAxzdZzvxXmFsSPMr1wXfE9cnZ57ufams72H6KZHu9LCd1nRG2rrnWdF70ZVBfAiVTPJmzER\nxe8iVa3222GDiv9f4Y+bzMws5SRhZmYpJwkzM0s5SZiZWcqjm8zMDPDopraYiFFIRVoZmTRasyOV\nGpWNWmpU1jqjavTSMU5jr5anI5j2ajlLIu/h+ICu5pLYXxh7mTmVo5vOGDqUnvuV7gXp6Kc5v/Yy\nh3Ru2tbjgzxYOWorGwX0QR7kXB1KW2vM4WXm6ng6gmmujnMoziiMASzQK+noppeZw2I9yv64pDC+\nWI+yN5ak575aD6Sjn44zi5Xak/7cK7Wn5ZYfjZodBdWolRFRo03ECKki7Ro1Vcajm8zMbMycJMzM\nLOUkYWZmKReuzcwMmEaF63YUkjPNrr8wVtu0bsxF59GaLUI3qlrLYcQu9aQF2hPM5F6tSNd8GKKb\nfVrG0ri7ML5Py1gc+9Jr79fStHB9jFk8qffxnvhOYfxJvS8tagOlRe/5HCldq+JCnmC71tATtxXG\nt2sNt0VPemxZa4z5HOEiPZMWn8tiI/HvxHsKY8eYVVm43heL03Mv1X7ujqWFsSG6K9eqyIraUF3Y\nHtFsgbtRK8Xu0Xo12NR6F2PVyvoY47Vexf+/8cdNZmaWcpIwM7OUk4SZmaVcuDYzM2AaFa7bUUjO\nbNO6cc+MLjKgvrQA2qztWsN1MTimY3ZqFStiR+V+u3U918QdhbFhutpeuF4U30qv/agWp4XtsthI\nPCtMz+OF0p9rIc8xoL709TCgvnSG70KeKy3wzuOFyjUfssIzUFqYHqaLy3WAb8WiwvjlOtDWwvUd\ncU167mt1FztiRRofcb12MxjXVe7XaJV2pgMJmrVG28c9a7tInwbaUhAvs07bCrf74yYzM0s5SZiZ\nWcpJwszMUi5cm5kZ4MJ1W0zEzOgircyWHq3Z2dONdqmnIwrXVa3CqwrbWXy/llaeO2thfiYvs0s9\n6TOdzxG2an0643+r1qezaOdzhB7tSmcXz+Ellmtv2pJ7ufZWtvOuKj5n8aXaX3nuTihcNzMzu1Er\ns7RHm4hZ20XaNZO7jAvXZmY2Zk4SZmaWcpIwM7OUC9dmZgZMo8J1O2ZAZ1pZR7oZp3Lh+iRdlWtc\nZ0VtoLToXRXfp2XpdUeunbVAn81Rdmh1+t/l7bzAFm1I11jfog3pusXz+AmrtSMtpM7maOVa0VlR\nG2qF7ay4DLBM+9J4WWwknl17mBkuXLdgItbfHqs+DRRu98dNZmaWcpIwM7OUk4SZmaVcuDYzM8CF\n67aYiJbeRbZrTccXrrMC7zBd3KVrWy5sl8XGG79XK9L7BtijlekzextHGVRvOlBhNkdbLlzP5ii9\nGkyLoLP5WemM7B7tqlwrOiseA5XF5VaPHaaLa3VXWpy+VndV3nenF67H2268SLtakJdx4drMzMbM\nScLMzFJOEmZmlnLh2szMgDdJ4VpSP/DvgR/XN90UUevtLGk98O+AYeBPIor7QZ8qheuxrk892k6t\nmpLCNdQKwK0WtstiVfEuhkuvXVaYhtrPnw0YmMWx0jWsZ3OUzernxugvjG9WP/1xY3psnwbSYuUs\njpXOyF6tHaWF27LCdhfDlcXlqlnRWXyYrsqZ4lNZuB7rutijTcQ62UU6qXDdcUkCCODWiLi1caOk\nC4GPAhcC5wBfl/TuiHhtCu7RzGxa6NSaxBve8gDLgDsi4tWIeBZ4Glg0qXdlZjbNdGqSWCvpu5K2\nS5pT3/YO4PmGfZ6n9o7CzMzaZEoK15LuB84uCH0aeJhf1SM2AgsiokfSNuDhiLi9fo7PA1+NiP82\n6twuXJuZtaBjCtcR8bvN7FdPBPfWvz0ILGwIv7O+7Q3a0bo7M6jeaVu4LttvPPHxHDuD4dJ1qMsK\n01CbpZ69fro5wTatS9dQn8WxlgvXszjGOm1L1zWeyVDpjOxeDZbOHi4rbHdxkuu1Oy0Ql8XGGx/v\nuRv3O9UK1+1oQV6mV8X/v+m4j5skLWj49iPA4/Wv7wH+QNJMSb8JXAAcmOz7MzObTjpxdNNmSRdT\nG+X0Q+CTABHxhKTdwBPASeCGONUmeZiZdZiOSxIR8UclsU3Apkm8HTOzaa3jPm4yM7PO4SRhZmYp\n924yMzOgg4bAttupMgR2IhYdGusw2maHzZYNMx1vvJljs5+ri2F2aHX67MqGuELtv2fWm2kmJ9iq\n9fTGzYXxboZaHgLbzRDrtZWboze9dllvpz4NlA6ZLBsi28Uwq7QzHQ66Sjsr+0K1Gh/PsaP3G+tw\n1lXa2dGLDnkIrJmZdTwnCTMzSzlJmJlZyknCzMxSThJmZpZykjAzs5SThJmZpTyZzszMAE+ma4tO\nn0w31nM0OwFvp1a1PKGt6t6aOTZ75l2cZFC96WugbLIcwID60slyXZxkizbQFxsL492cYJM2clNs\nKIxv0kY2xE2FsZkMsUFb2Bh9hfEZDJdOtluvrelEOyifbDeDYdZoezopbI22V65V0Wp8tXaUToIr\nm+Q3er+xToyruu9mz+HJdGZmNm05SZiZWcpJwszMUk4SZmaWcpIwM7OUk4SZmaWcJMzMLOUkYWZm\nKScJMzNLuS2HmZkBbsvRFu1syzHe87ZyjmaPqdqvLN5M64xW16HuYphtWsfauKUwvk3r0rYbAFu1\nPm27MYPh0jWsuxhuuS1HF8P0a3O6BnYXJ0vbdmzQlrRlB9TadtwSa9Nrj3f97PG0/Chra1EVH+t+\n4z2mHefIzuu2HGZm1vGcJMzMLOUkYWZmKScJMzNLOUmYmVnKScLMzFJOEmZmlnKSMDOzlJOEmZml\nnCTMzCzlJGFmZiknCTMzSzlJmJlZyq3CzcwMmEatwrM20u0woL62tCavapfdrnM0e0xZu26oPZdW\n4wPqS1t9A6WtwLsYZqvWp+3Ay1qBA2zRhjdtq/AsNhLPWol3Mcw6bUtbiZfFRuJZm3Ggsg151bHN\ntMwua1c+kce04xzZecueSzv0aaBwuz9uMjOzlJOEmZmlnCTMzCw1JUlC0u9L+r+ShiW9b1RsvaSn\nJH1f0tKG7ZdIerwe+0+Tf9dmZtPPVL2TeBz4CPC/GjdKuhD4KHAhcBUwIGmk2v45oCciLgAukHTV\nJN6vmdm0NCVJIiK+HxE/KAgtA+6IiFcj4lngaeADkhYAsyPiQH2//wJcMzl3a2Y2fXVaTeIdwPMN\n3z8PnFOw/WB9u5mZtVHb5klIuh84uyB0U0Tc267rmpnZxGlbkoiI323hsIPAwobv30ntHcTB+teN\n2w9mJznQf98vvz5nyXmcs+T8Fm7FzOzU9fQDB3nmgfR/o7/UCTOuG6eB3wPsknQrtY+TLgAORERI\n+qmkDwAHgOuAv85OuKj/ynber5nZm975S87h/CW/+tR+32cfKdxvqobAfkTSc8ClwP+QtBcgIp4A\ndgNPAHuBG+JXzaVuAD4PPAU8HRFfm/w7NzObXqbknURE3AncmcQ2AZsKtj8K/Ms235qZmTXotNFN\nZmbWQZwkgIMPPD3Vt2BmNiGebqIYPRZOEsDBB56Z6lswM5sQzYxYGgsnCTMzSzlJmJlZysuXmpkZ\nULx86SmXJMzMbOL44yYzM0s5SZiZWWraJQmvijc+kvolPS/p7+p/rm6IFT4/A0lX1Z/LU5JunOr7\n6VSSnpX09/XX1oH6trmS7pf0A0n7JM2Z6vucSpK+IOmwpMcbtqXPaLy/l9MuSeBV8cYrgFsj4r31\nP3shfX7T8fX1BpK6gNuoPZcLgWslvWdq76pjBbCk/tpaVN/2KeD+iHg38I3699PZF6m9lhoVPqOJ\n+L2cdr/EXhVvQrxhBATFz29RwX7T0SJqTSmfjYhXga9Qe15WbPTr68PAl+pff4lp/vsXEd8GXhq1\nOXtG4/69nHZJooRXxWveWknflbS94W1t9vys9hyea/jezyYXwNclPSLpE/VtZ0XE4frXh4GzpubW\nOlr2jMb9e9kJ60lMOK+KNz4lz+/T1D56+8v69xuBLUBPciqPr67xc2jeb0XEIUm/Dtwv6fuNwfra\nMn6eJZp4RmN6fqdkkpjKVfFOBc0+P0mfB0aSbtHzO6Wf0xiMfjYLef3f7qwuIg7V//1jSXdS+2jk\nsKSzI+JH9Y9/j0zpTXam7BmN+/dyun/cNHpVvD+QNFPSb/KrVfF+BPxU0gfqhezrgLum4F47Qv0F\nOOIj1AYCQPL8Jvv+OtQj1AY8vEvSTGqFxHum+J46jqRZkmbXvz4dWErt9XUP8PH6bh9nGv/+lcie\n0bh/L0/JdxJlJH2E2tKn86itivd3EXF1RDwhaWRVvJO8cVW8HcBpwFen+ap4myVdTO0t6w+BT0Jt\nVcGS5zetRcRJSWuA+4AuYHtEPDnFt9WJzgLurA8qnAHcHhH7JD0C7JbUAzwLrJi6W5x6ku4ALgfm\n1Vf4/Avgryh4RhPxe+m2HGZmlpruHzeZmVkJJwkzM0s5SZiZWcpJwszMUk4SZmaWcpIwM7OUk4SZ\nmaWcJMzMLOUkYdZmkv51vWtut6TTJX2v3uffrON5xrXZJJC0EXgrtdYuz0XE5im+JbOmOEmYTQJJ\n/4xao7/jwGXua2VvFv64yWxyzANOB95G7d2E2ZuC30mYTQJJ9wC7gHOBBRGxdopvyawp065VuNlk\nk/RHwFBEfKW+CP2DkpZExANTfGtmlfxOwszMUq5JmJlZyknCzMxSThJmZpZykjAzs5SThJmZpZwk\nzMws5SQHnaGhAAAAD0lEQVRhZmYpJwkzM0v9f+5iSBL3nlObAAAAAElFTkSuQmCC\n", + "text": [ + "" + ] + }, + { + "metadata": {}, + "output_type": "display_data", + "png": "iVBORw0KGgoAAAANSUhEUgAAAYkAAAEZCAYAAABiu9n+AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAE71JREFUeJzt3W+sZVd53/Hvz0MNxHWEHadjA+5MUOwoo1SCQsdBSfG8\nSNxJJTBOFROkJE7qIlSEk6gkARs5nWQkE6Ni0RClL4Ihbgsmo7Qm4wYnHqMYQ1LqOE0TtxMDjpip\nx7LHNHEEbYH6z9MXZxuuZ/a655x7zj5/5n4/0mjO2WutZy1v3etn9n72n1QVkiT1OWvZC5AkrS6T\nhCSpySQhSWoySUiSmkwSkqQmk4QkqckkIU0hyU8m+fSG719Jsnt5K5KGZZKQTpHk+5P8UZK/SfJX\nST6T5DV9favq3Ko6Nuf5357kgSRfS/LhU9p2J3m2S07P/Xn3POeXNnrBshcgrZIk3wr8J+CtwCHg\nhcA/BL6+wGU8ChwE/hHw4kafby3vhNUCeCQhPd+lQFXVb9XI16rqSFU92Ne5+1f9K7rPL07yviTH\nuqOQTyd5Udf2vd3RyZNJ/luSy1sLqKo7qup3gL/aZJ3+7moh/EGTnu9zwDNJfjPJ/iTnTTH2XwGv\nAl4LnA/8PPBskpcxOjr55ao6D/g54D8kuWBMvGzSdjzJI0k+lOTbplijNBWThLRBVX0F+H6ggN8A\nnkjyO0n+zmbjkpwF/BTwM1X1WFU9W1Wfrar/B/wY8Imq+r1ujnuAB4B/PG45Pdu+BLwG+LvAq4Fz\ngY9M/B8oTckkIZ2iqh6qqp+qqouB7wFeCrx/zLALgBcBf9nTtgv4ke5U05NJngS+D7hwTMzTjiSq\n6v9U1X/tktATwNuBK5KcMyaWtCUmCWkTVfU54DZGyWIz/wv4GvCdPW3/E/h3VXXehj/nVtV7x00/\nxVL9XdYg/MGSNkjyXUn+RVdHIMnFwJuB/7zZuKp6FvgQcEuSi5LsSPLaJGcD/x54fZIruu0vSrLv\nuTl61rCjK3i/ANiR5IVJdnRte7s1ntXVIn4V+IPuNJk0dyYJ6fm+AlwG/Jck/5tRcvhz4B1de/H8\nf+Fv/PxzwIPAHzO6Muk9wFlVdQK4ErgBeILRkcU7aP/+3Qj8X+CdjOoZXwWeuxfiFcBdwJe7ub7K\nKIlJg4iXWkuSWjySkCQ1mSQkSU0mCUlSk0lCktR0xj3gL4mVeEnagqo67QbOMy5JANxQN07V/74D\nn+J1B5rPW9vUTTk49XyLiruVGJOOGddvlvZ1jT3J3DfWDc2xB3PTlttnGbvOsaftN+uYIWLMO+6n\nDtzH5Qdet6U5+3i6SZLUZJKQJDUtNUl0jzk+meTBDdvOT3IkyeeT3J3kJRvark/yhSQPJbliXuvY\ntW/XvEJJ0lLN+/9nyz6S+DCw/5Rt7wKOVNWlwCe77yTZA7wJ2NON+fXu8cwz27Vv9zzCSNLS7T6T\nkkRVfRp48pTNb2D01E26v9/Yfb4SuL2qnureKfwwsHcR65Sk7WrZRxJ9dlbVye7zSWBn9/mlwIkN\n/U4AvU/RlCTNxyomiW/oXvS+2X0P3hMhSQNaxfskTia5sKoeT3IRo0crAzwKXLyh38u7bae578Cn\nvvF5175d1hwk6RTH7j3O8XuPj+23ikniMHANcHP398c3bP9oklsYnWa6BLi/L8BWb4yTpO1i975d\nzyty3/dLn+ntt9QkkeR24HLggiSPAL8I/ApwKMm1wDHgaoCqOprkEHAUeBp4W/kyDEka1FKTRFW1\n3qj1A43+NwH9945LkuZupQvXkqTlMklIkppMEpKkJpOEJKnJJCFJajJJSJKacqbdauDrSyVpa3x9\n6QB8fen829c19iRzr+prQNc19rT9Zh0zRIxFxh03Zx9PN0mSmkwSkqQmk4QkqckkIUlqMklIkppM\nEpKkJpOEJKnJJCFJajJJSJKaTBKSpCaThCSpySQhSWryKbCSJMCnwA7Cp8DOv31dY08y96o+TXVd\nY0/bb9YxQ8RYZNxxc/bxdJMkqckkIUlqMklIkppMEpKkJpOEJKnJJCFJajJJSJKaTBKSpCaThCSp\nySQhSWry2U2SJMBnNw3CZzfNv31dY08y96o+A2ldY0/bb9YxQ8RYZNxxc/bxdJMkqckkIUlqMklI\nkppMEpKkJpOEJKnJJCFJavI+CUkS4H0Sg/A+ifm3r2vsSeZe1fsN1jX2tP1mHTNEjEXGHTdnH083\nSZKaVvZIIskx4MvAM8BTVbU3yfnAbwG7gGPA1VX1N0tbpCSd4Vb5SKKAfVX1qqra2217F3Ckqi4F\nPtl9lyQNZJWTBMCpRZQ3ALd1n28D3rjY5UjS9rLKSaKAe5I8kOQt3badVXWy+3wS2LmcpUnS9rCy\nl8AmuaiqHkvy7cAR4DrgcFWdt6HPX1fV+aeMW83/IElacWt1CWxVPdb9/aUkdwB7gZNJLqyqx5Nc\nBDzRN9ZLYLcew0tgvQR2nWJP22/WMUPEWGTccXP2WcnTTUm+Jcm53edzgCuAB4HDwDVdt2uAjy9n\nhZK0PazqkcRO4I4kMFrjR6rq7iQPAIeSXEt3CezylihJZ76VTBJV9UXglT3b/xr4gcWvSJK2p5Ut\nXG+VhWtJ2pq1KlzPwsL11mNYuLZwvU6xp+0365ghYiwy7rg5+6xk4VqStBpMEpKkJpOEJKnJwrUk\nCbBwPQgL1/NvX9fYk8y9qgXgdY09bb9ZxwwRY5Fxx83Zx9NNkqQmk4QkqckkIUlqsnAtSQIsXA/C\nwvX829c19iRzr2oBeF1jT9tv1jFDxFhk3HFz9vF0kySpySQhSWoySUiSmixcS5IAC9eDsHA9//Z1\njT3J3KtaAF7X2NP2m3XMEDEWGXfcnH083SRJajJJSJKaTBKSpCaThCSpyaubJEmAVzcNwqub5t++\nrrEnmXtVrxJa19jT9pt1zBAxFhl33Jx9PN0kSWoySUiSmkwSkqQmC9eSJMDC9SAsXM+/fV1jTzL3\nqhaA1zX2tP1mHTNEjEXGHTdnH083SZKaTBKSpCaThCSpycK1JAmwcD0IC9fzb1/X2JPMvaoF4HWN\nPW2/WccMEWORccfN2cfTTZKkJpOEJKnJJCFJarJwLUkCLFwPwsL1/NuXHfuddaC37eu8kPfnen62\n3tPb/gw7+EB+gevqvb3tH8gv8N66rjn2+ryf99TPNuc+kJs5UO/sbT+Qm1e2uGzhenXijpuzj6eb\nJElNJglJUpNJQpLUZOFakgScIYXrJPuB9wM7gA9W1c2n9rFwvfUYFq4P8o462Nv2eb6LO3M1r69D\nve135ihwoPvT50CzZU+9nqtzJ4fq9b3tn+NSbsz7OFjv6G2/Me9b2eKyhevViTtuzj5rdbopyQ7g\n14D9wB7gzUm+e7mrkqQz19gkkeSnk5y3iMVMYC/wcFUdq6qngI8BVy55TZJ0xprkSGIn8MdJDiXZ\nn+S0c1YL9DLgkQ3fT3TbJEkDmKhwneQs4ArgJ4HXAIeAW6vqLwdd3enr+CfA/qp6S/f9x4DLqr55\nh5KFa0nami0Xrqvq2SSPAyeBZ4DzgN9Ock9V/fx8l7mpR4GLN3y/mNHRxPNYuN56DAvXFq7nOXbo\n2NP2m3XMEDEWGXfcnH0mqUn8TJI/Ad4L/CHwPVX1z4FXAz88z0VO4AHgkiS7k5wNvAk4vOA1SNK2\nMcmRxPnAD1fV8Y0bu6OL/n/2DKSqnk7yduD3GV0Ce2tV/cUi1yBJ28nYJFFV/3KTtqPzXc54VXUX\ncNei55Wk7cg7riVJwBlyx/UkLFxvPYaFax8VPs+xQ8eett+sY4aIsci44+bss1Z3XEuSFsskIUlq\nMklIkposXEuSAAvXg7BwPf/2dY09ydyrWgBe19jT9pt1zBAxFhl33Jx9PN0kSWoySUiSmkwSkqQm\nk4QkqcmrmyRJgFc3DcKrm+bfvq6xJ5l7Va8SWtfY0/abdcwQMRYZd9ycfTzdJElqMklIkppMEpKk\nJgvXkiTAwvUgLFzPv31dY08y96oWgNc19rT9Zh0zRIxFxh03Zx9PN0mSmkwSkqQmk4QkqcnCtSQJ\nsHA9CAvX829f19iTzL2qBeB1jT1tv1nHDBFjkXHHzdnH002SpCaThCSpySQhSWqycC1JAixcD8LC\n9fzb1zX2JHOvagF4XWNP22/WMUPEWGTccXP28XSTJKnJJCFJajJJSJKaLFxLkgAL14OwcD3/9nWN\nPcncq1oAXtfY0/abdcwQMRYZd9ycfTzdJElqMklIkppMEpKkJgvXkiTAwvUgLFzPv31dY08y96oW\ngNc19rT9Zh0zRIxFxh03Zx9PN0mSmkwSkqQmk4QkqcnCtSQJWJPCdZIDwD8DvtRtuqGq7urargf+\nKfAM8NNVdXdfDAvXW49h4drC9TrFnrbfrGOGiLHIuOPm7LNySQIo4JaqumXjxiR7gDcBe4CXAfck\nubSqnl3CGiVpW1jVmsRphzzAlcDtVfVUVR0DHgb2LnRVkrTNrGqSuC7JnyW5NclLum0vBU5s6HOC\n0RGFJGkgSylcJzkCXNjT9G7gs3yzHnEQuKiqrk3yAeCzVfWRLsYHgU9U1X88JbaFa0nagpUpXFfV\nD07Sr0sEd3ZfHwUu3tD88m7baSxcbz2GhWsL1+sUe9p+s44ZIsYi446bs8/KnW5KctGGr1cBD3af\nDwM/muTsJN8BXALcv+j1SdJ2sopXN92c5JWMrnL6IvBWgKo6muQQcBR4GnhbnWk3eUjSilm5JFFV\nP7FJ201A/zGRJGnuVu50kyRpdZgkJElNPrtJkgSs0CWwQ/MS2K3H8BJYL4Fdp9jT9pt1zBAxFhl3\n3Jx9PN0kSWoySUiSmkwSkqQmk4QkqckkIUlqMklIkppMEpKkJm+mkyQB3kw3CG+mm3/7usaeZO5V\nvSltXWNP22/WMUPEWGTccXP28XSTJKnJJCFJajJJSJKaTBKSpCaThCSpySQhSWoySUiSmkwSkqQm\nk4QkqcnHckiSAB/LMQgfyzH/9nWNPcncq/p4i3WNPW2/WccMEWORccfN2cfTTZKkJpOEJKnJJCFJ\najJJSJKaTBKSpCaThCSpySQhSWoySUiSmkwSkqQmk4QkqckkIUlqMklIkppMEpKkJh8VLkkCfFT4\nIHxU+Pzb1zX2JHOv6iO31zX2tP1mHTNEjEXGHTdnH083SZKaTBKSpCaThCSpaSlJIsmPJPkfSZ5J\n8vdPabs+yReSPJTkig3bX53kwa7tXy9+1ZK0/SzrSOJB4Crgvo0bk+wB3gTsAfYDv57kuWr7vwGu\nrapLgEuS7F/geiVpW1pKkqiqh6rq8z1NVwK3V9VTVXUMeBi4LMlFwLlVdX/X798Cb1zMaiVp+1q1\nmsRLgRMbvp8AXtaz/dFuuyRpQIPdJ5HkCHBhT9MNVXXnUPNKkuZnsCRRVT+4hWGPAhdv+P5yRkcQ\nj3afN25/tBXkvgOf+sbnXft2sWvf7i0sRZLOXMfuPc7xe4+P7bcKd1xvvA38MPDRJLcwOp10CXB/\nVVWSLye5DLgf+HHgV1sBX3fg8iHXK0lrb/e+Xezet+sb3+/7pc/09lvWJbBXJXkE+F7gd5PcBVBV\nR4FDwFHgLuBt9c2HS70N+CDwBeDhqvq9xa9ckraXpRxJVNUdwB2NtpuA0x4iUlV/Avy9gZcmSdpg\n1a5ukiStEJMEcPzeY8tegiTNxbEJitHTMEnARBV+SVoH8/7/mUlCktRkkpAkNfn6UkkS0P/60jMu\nSUiS5sfTTZKkJpOEJKlp2yUJ34o3myQHkpxI8qfdnx/a0Na7/wRJ9nf75QtJ3rns9ayqJMeS/Hn3\ns3V/t+38JEeSfD7J3Ulesux1LlOSDyU5meTBDdua+2jW38ttlyTwrXizKuCWqnpV9+cuaO6/7fjz\ndZokO4BfY7Rf9gBvTvLdy13VyipgX/eztbfb9i7gSFVdCnyy+76dfZjRz9JGvftoHr+X2+6X2Lfi\nzcVpV0DQv//29vTbjvYyeijlsap6CvgYo/2lfqf+fL0BuK37fBvb/Pevqj4NPHnK5tY+mvn3ctsl\niU34VrzJXZfkz5LcuuGwtrX/NNoPj2z47r5pK+CeJA8keUu3bWdVnew+nwR2LmdpK621j2b+vVyF\n90nMnW/Fm80m++/djE69/XL3/SDwPuDaRiivrx5xP0zu+6rqsSTfDhxJ8tDGxu7dMu7PTUywj6ba\nf2dkkljmW/HOBJPuvyQfBJ5Lun3774zeT1M4dd9czPP/dadOVT3W/f2lJHcwOjVyMsmFVfV4d/r3\niaUucjW19tHMv5fb/XTTqW/F+9EkZyf5Dr75VrzHgS8nuawrZP848PElrHUldD+Az7mK0YUA0Nh/\ni17finqA0QUPu5OczaiQeHjJa1o5Sb4lybnd53OAKxj9fB0Grum6XcM2/v3bRGsfzfx7eUYeSWwm\nyVWMXn16AaO34v1pVf1QVR1N8txb8Z7m9Lfi/SbwYuAT2/yteDcneSWjQ9YvAm+F0VsFN9l/21pV\nPZ3k7cDvAzuAW6vqL5a8rFW0E7iju6jwBcBHquruJA8Ah5JcCxwDrl7eEpcvye3A5cAF3Rs+fxH4\nFXr20Tx+L30shySpabufbpIkbcIkIUlqMklIkppMEpKkJpOEJKnJJCFJajJJSJKaTBKSpCaThDSw\nJP+ge2ruC5Ock+S/d8/5l1aed1xLC5DkIPAiRo92eaSqbl7ykqSJmCSkBUjytxg96O+rwGt9rpXW\nhaebpMW4ADgH+NuMjiakteCRhLQASQ4DHwVeAVxUVdcteUnSRLbdo8KlRUvyE8DXq+pj3Uvo/yjJ\nvqq6d8lLk8bySEKS1GRNQpLUZJKQJDWZJCRJTSYJSVKTSUKS1GSSkCQ1mSQkSU0mCUlS0/8HcCKP\nA2mbunsAAAAASUVORK5CYII=\n", + "text": [ + "" + ] + }, + { + "metadata": {}, + "output_type": "display_data", + "png": "iVBORw0KGgoAAAANSUhEUgAAAYkAAAEZCAYAAABiu9n+AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAFb1JREFUeJzt3X+QZWV95/H3h2ZBYUkBITuAsjOxFrZCmSpc3UEqWZk/\nEsStikisaHQlJDsRaylIrJ0EBYukDVVjxlooNlikinU0LAKGWhcluxIZrCC60RCyMcElREk5LLAw\nGAOluyhC890/7hnpmbnPvd19+/7o7verqul7z3PO93nm1L18+5zvc85JVSFJUj+HTXsAkqTZZZKQ\nJDWZJCRJTSYJSVKTSUKS1GSSkCQ1mSSkZUjyy0m+uOj9d5Nsmd6IpPEySUgHSfLTSf40yTNJvp3k\nS0le12/dqjqmqvaucv+XJLk/yfeTfPygti1JXuyS0/6fD6xm/9Jih097ANIsSfIjwH8D3gPcBhwJ\n/CvguQkO43HgKuCNwMsb6/xIeSWsJsAjCelApwFVVX9YPd+vqj1V9UC/lbu/6l/VvX55kquT7O2O\nQr6Y5GVd2+u7o5Onk3w1ydmtAVTV7VX1GeDbA8bpd1cT4QdNOtDfAgtJ/iDJuUmOW8a2/wF4DXAW\ncDzwm8CLSV5B7+jkd6rqOOA3gE8lOWFIvAxoeyTJo0k+luRHlzFGaVlMEtIiVfVd4KeBAv4T8FSS\nzyT5J4O2S3IY8CvAr1fVE1X1YlV9pap+ALwL+GxV/XHXx93A/cC/HjacPsu+BbwO+KfAa4FjgJuX\n/A+UlskkIR2kqh6qql+pqlOAVwMnA9cO2ewE4GXA3/Vp2wz8Qneq6ekkTwM/BZw4JOYhRxJV9f+q\n6n92Segp4BLgnCRHD4klrYhJQhqgqv4WuJFeshjk74HvA/+sT9v/Bm6qquMW/RxTVR8e1v0yhup3\nWWPhB0taJMk/T/LvuzoCSU4B3gF8edB2VfUi8DHgmiQnJZlLclaSI4BPAD+X5Jxu+cuSbNvfR58x\nzHUF78OBuSRHJpnr2rZ2Yzysq0X8HvAn3WkyadWZJKQDfRc4E/izJP+XXnL4a2BH114c+Bf+4te/\nATwA/Dm9mUkfAg6rqseA84ArgKfoHVnsoP39uxJ4FngfvXrG94D910K8CrgT+E7X1/foJTFpLOJU\na0lSi0cSkqQmk4QkqckkIUlqMklIkprW3Q3+kliJl6QVqKpDLuBcd0kC4Mq6Ylnrf2H+Xs6ef8OK\n+roqO5fd36TiriTGUrcZtt4o7Ws19lL6vqKubG67M1etuH2Ubddy7OWuN+o244ix2nHvnf8Cb5hv\n3j9yYJ/9eLpJktRkkpAkNU01SXS3Od6X5IFFy45PsifJ15PcleTYRW2XJ/lGkoeSnLNa49i8bfNq\nhZKkqVrt/59N+0ji48C5By17P7Cnqk4DPt+9J8npwNuB07ttru9uzzyyLSYJSevE5m1bVjXeVJNE\nVX0RePqgxW+md9dNut9v6V6fB9xaVc93zxR+GNg6iXFK0kY17SOJfjZV1b7u9T5gU/f6ZOCxRes9\nBvS9i6YkaXXMYpL4oe5B74Oue/CaCEkao1m8TmJfkhOr6skkJ9G7tTLA48Api9Z7ZbfsEF+Yv/eH\nrzdv22zNQZIO8sg9e3nknkeGrjeLSeIO4EJgV/f704uW35LkGnqnmU4F7usXYKUXxknSRrF525YD\nitxf+uC9fdebapJIcitwNnBCkkeB3wJ+F7gtyXZgL/A2gKp6MMltwIPAC8DF5cMwJGmsppokqqr1\nRK2faay/E9g5vhFJkhab6cK1JGm6TBKSpCaThCSpySQhSWoySUiSmkwSkqSmrLdLDXx8qSStjI8v\nHQMfX7r67Ws19lL6ntXHgK7V2Mtdb9RtxhFjknGH9dmPp5skSU0mCUlSk0lCktRkkpAkNZkkJElN\nJglJUpNJQpLUZJKQJDWZJCRJTSYJSVKTSUKS1GSSkCQ1eRdYSRLgXWDHwrvArn77Wo29lL5n9W6q\nazX2ctcbdZtxxJhk3GF99uPpJklSk0lCktRkkpAkNZkkJElNJglJUpNJQpLUZJKQJDWZJCRJTSYJ\nSVKTSUKS1OS9myRJgPduGgvv3bT67Ws19lL6ntV7IK3V2Mtdb9RtxhFjknGH9dmPp5skSU0mCUlS\nk0lCktRkkpAkNZkkJElNJglJUpPXSUiSAK+TGAuvk1j99rUaeyl9z+r1Bms19nLXG3WbccSYZNxh\nffbj6SZJUtPMHkkk2Qt8B1gAnq+qrUmOB/4Q2AzsBd5WVc9MbZCStM7N8pFEAduq6jVVtbVb9n5g\nT1WdBny+ey9JGpNZThIABxdR3gzc2L2+EXjLZIcjSRvLLCeJAu5Ocn+Sd3fLNlXVvu71PmDTdIYm\nSRvDzE6BTXJSVT2R5MeAPcClwB1Vddyidf6hqo4/aLvZ/AdJ0oxbU1Ngq+qJ7ve3ktwObAX2JTmx\nqp5MchLwVL9tnQK78hhOgXUK7FqKvdz1Rt1mHDEmGXdYn/3M5OmmJEclOaZ7fTRwDvAAcAdwYbfa\nhcCnpzNCSdoYZvVIYhNwexLojfHmqroryf3AbUm2002Bnd4QJWn9m8kkUVXfBM7os/wfgJ+Z/Igk\naWOa2cL1Slm4lqSVWVOF61FYuF55DAvXFq7XUuzlrjfqNuOIMcm4w/rsZyYL15Kk2WCSkCQ1mSQk\nSU0WriVJgIXrsbBwvfz2ORaYzy7m63192we1jdq+lG3XauH6fTXf3HZX5lfcPsq2cyxYuJ6huMP6\n7MfTTZKkJpOEJKnJJCFJarJwLUkCLFyPxUYuXE+zuDyN2AvMTa1wvcDc2IrLw9qnHdvC9WRYuJYk\nLZtJQpLUZJKQJDVZuJYkARaux2I9Fq4HFXb3m9Xi8hwvcGWu5qra0bd9UNv+9nEVtgcVrocVpmF4\nAXhH9S88AlydK5vth7Mw04XrQe2L17NwPXqf/Xi6SZLUZJKQJDWZJCRJTSYJSVKTs5skSYCzm8bC\n2U2r3z6fXUNnIK20/cpczYfqvc1tL8+1zfYXmBsY+wccueLZTc9xxMAZSAsczrW5nPfWh/q2D2rb\n377S2U+D2va3O7tpdTm7SZK0JpgkJElNJglJUpOFa0kSYOF6LCxcL799NW6dsdLi8+W5lg/Xpc1t\nL8t1XF0X921bYI7Lcl1z+2c5amhBvlWEfZajuC6XcWl9uNn39dnBxXV13/brs6O5LcB1uWxoYXtQ\nUXxWb/mxeD0L16P32Y+nmyRJTSYJSVKTSUKS1GThWpIEWLgei1kvXC+lCL3YsIL0UtYbVpiG0a+K\nbhWP51hgR65vFp935Ho+UtubsS/Jbm6oC5rtF+WmZvszHDuwsH1ZrmsWl5/hWG7KRVxQNzT7HtR+\nUy5ie32kue3uXNIses+xMLBovpSi9yhXc69W4Xop6x28jYXrA/vsx9NNkqQmk4QkqckkIUlqsnAt\nSQIsXI+FhetDTbNwfQTP8d7cwLV1Ud/29+YGdtc7m7G35xY+U+f0bTuJ/8PWfI376tV927/NCbwp\n93Bnbevb/qbcw7a6s2/b3/OjfC1beXXd17f9ZJ7grpzHOfWZvu135TzeWbv7tgHcku1cVNf2bTuS\nH1i4XiEL15KkDc0kIUlqMklIkposXEuSgHVSuE5yLnAtMAd8tKp2HbyOheuXYqzXwnXriuojeY5L\nsrt5VfUl2c0n6q3N2O/Kp/hyndG37fV3fZW8Eepz/bf94Bthnt5PP/OL/nuIm38b/k3g5v5/45zx\nzq/w1ZzFGfXlvu1fzVm8tT7R6Bk+lXc1r8g+kueG3obcwnV/Fq5nTJI54CPAucDpwDuS/MR0RyVJ\n69fQJJHk15IcN4nBLMFW4OGq2ltVzwOfBM6b8pgkad1aypHEJuDPk9yW5Nwkh5yzmqBXAI8uev9Y\nt0ySNAZLKlwnOQw4B/hl4HXAbcDuqvq7sY7u0HG8FTi3qt7dvX8XcGbVS1dWWbiWpJVZceG6ql5M\n8iSwD1gAjgP+S5K7q+o3V3eYAz0OnLLo/Sn0jiYOYOH6pRgWrg9k4fpQFq5XzsI1kOTXk/wF8GHg\nfwCvrqp/B7wW+PnVHOQS3A+cmmRLkiOAtwN3THgMkrRhLOVI4njg56vqkcULu6OLnxvPsPqrqheS\nXAJ8jt4U2N1V9TeTHIMkbSRDk0RV/faAtgdXdzjDVdWdQP+7pEmSVpVXXEuSgHVyxfVSWLh+KcZ6\nLVx7q/ADeatwC9er0Wc/a+qKa0nSZJkkJElNJglJUpOFa0kSYOF6LCxcH2qahes5FtiR65tXZO/I\n9c2rsaF3RfYNdUGz/aLc1Gx/hmO5LNc1x3ZZrmsWh5/hWG7KRVxQNzT7HtR+Uy5qXlENsDuXNK+o\nnmPBwvUKWbiWJG1oJglJUpNJQpLUZJKQJDU5u0mSBDi7aSxmfXbTcmMsdUbUsFlQo8x+WsrMqGGz\nn1rtg2ZGQW8GUmtm1AJzA2cvPctRA//d89nVnIHzLEcNnGG0wNzQZz60toWlzVBqtY8ye+lwFgbO\nPBo2K2k5s5uWOxtoNWYQObtJkrShmSQkSU0mCUlSk4VrSRJg4XosLFyvfvt8dq34lh7D2kcper/A\n3MDYP+DIgfv8quxsFiOf44iBBeAFDh9aXB7nrTOGbbvS4rOF68nGHdZnP55ukiQ1mSQkSU0mCUlS\nk4VrSRJg4XosLFyvfvu0r+Zead8LzK24cL3A3EhF3F2Zn9mroi1cz07cYX324+kmSVKTSUKS1GSS\nkCQ1WbiWJAEWrsdiPRaul7LNsAL3tArX44w9rDANg/ffoMI1DC5WDitsz3JxedTYSyngrqTQa+H6\n0D778XSTJKnJJCFJajJJSJKaLFxLkgAL12OxkQvXKy3gzrEw04Xrlf67hrWPUrge1r4zV81k4XqO\nhaHjHmWfLHe9UbcZR4xJxh3WZz+ebpIkNZkkJElNJglJUpOFa0kSYOF6LCxcr377Wo29lL7HWbje\niLGXu96o24wjxiTjDuuzH083SZKaTBKSpCaThCSpycK1JAlYI4XrJPPArwLf6hZdUVV3dm2XA/8W\nWAB+raru6hfDwvXKY1i4tnC9lmIvd71RtxlHjEnGHdZnPzOXJIACrqmqaxYvTHI68HbgdOAVwN1J\nTquqF6cwRknaEGa1JnHIIQ9wHnBrVT1fVXuBh4GtEx2VJG0ws5okLk3yV0l2Jzm2W3Yy8NiidR6j\nd0QhSRqTqRSuk+wBTuzT9AHgK7xUj7gKOKmqtie5DvhKVd3cxfgo8Nmq+q8HxbZwLUkrMDOF66r6\n2aWs1yWCP+rePg6csqj5ld2yQ1i4XnkMC9cWrtdS7OWuN+o244gxybjD+uxn5k43JTlp0dvzgQe6\n13cAv5jkiCQ/DpwK3Dfp8UnSRjKLs5t2JTmD3iynbwLvAaiqB5PcBjwIvABcXOvtIg9JmjEzlySq\n6pcGtO0Edk5wOJK0oc3c6SZJ0uwwSUiSmrx3kyQJmKEpsOPmFNiVx3AKrFNg11Ls5a436jbjiDHJ\nuMP67MfTTZKkJpOEJKnJJCFJajJJSJKaTBKSpCaThCSpySQhSWryYjpJEuDFdGPhxXSr375WYy+l\n71m9KG2txl7ueqNuM44Yk4w7rM9+PN0kSWoySUiSmkwSkqQmk4QkqckkIUlqMklIkppMEpKkJpOE\nJKnJJCFJavK2HJIkwNtyjIW35Vj99rUaeyl9z+rtLdZq7OWuN+o244gxybjD+uzH002SpCaThCSp\nySQhSWoySUiSmkwSkqQmk4QkqckkIUlqMklIkppMEpKkJpOEJKnJJCFJajJJSJKaTBKSpCZvFS5J\nArxV+Fh4q/DVb1+rsZfS96zecnutxl7ueqNuM44Yk4w7rM9+PN0kSWoySUiSmkwSkqSmqSSJJL+Q\n5H8lWUjyLw5quzzJN5I8lOScRctfm+SBru0/Tn7UkrTxTOtI4gHgfODexQuTnA68HTgdOBe4Psn+\navvvA9ur6lTg1CTnTnC8krQhTSVJVNVDVfX1Pk3nAbdW1fNVtRd4GDgzyUnAMVV1X7fefwbeMpnR\nStLGNWs1iZOBxxa9fwx4RZ/lj3fLJUljNLbrJJLsAU7s03RFVf3RuPqVJK2esSWJqvrZFWz2OHDK\novevpHcE8Xj3evHyx1tBvjD/Uqlj87bNbNm2eQVDkaT165F79vLIPY8MXW8WrrhefBn4HcAtSa6h\ndzrpVOC+qqok30lyJnAfcAHwe62AZ8+/YZzjlaQ1b/O2LWzetuWH77/0wXv7rjetKbDnJ3kUeD3w\n35PcCVBVDwK3AQ8CdwIX10s3l7oY+CjwDeDhqvrjyY9ckjaWqRxJVNXtwO2Ntp3Azj7L/wL4yTEP\nTZK0yKzNbpIkzRCTBLB3CcUbSVoLHrln76rGM0nAkir8krQWrPb/z0wSkqQmk4QkqcnHl0qSgP6P\nL113SUKStHo83SRJajJJSJKaNlyS8Kl4o0kyn+SxJH/Z/bxpUVvf/SdIcm63X76R5H3THs+sSrI3\nyV93n637umXHJ9mT5OtJ7kpy7LTHOU1JPpZkX5IHFi1r7qNRv5cbLkngU/FGVcA1VfWa7udOaO6/\njfj5OkSSOeAj9PbL6cA7kvzEdEc1swrY1n22tnbL3g/sqarTgM937zeyj9P7LC3Wdx+txvdyw32J\nfSreqjhkBgT999/WPuttRFvp3ZRyb1U9D3yS3v5Sfwd/vt4M3Ni9vpEN/v2rqi8CTx+0uLWPRv5e\nbrgkMYBPxVu6S5P8VZLdiw5rW/tPvf3w6KL37pu2Au5Ocn+Sd3fLNlXVvu71PmDTdIY201r7aOTv\n5Sw8T2LV+VS80QzYfx+gd+rtd7r3VwFXA9sboZxf3eN+WLqfqqonkvwYsCfJQ4sbu2fLuD8HWMI+\nWtb+W5dJYppPxVsPlrr/knwU2J90++2/db2fluHgfXMKB/51p05VPdH9/laS2+mdGtmX5MSqerI7\n/fvUVAc5m1r7aOTv5UY/3XTwU/F+MckRSX6cl56K9yTwnSRndoXsC4BPT2GsM6H7AO53Pr2JANDY\nf5Me34y6n96Ehy1JjqBXSLxjymOaOUmOSnJM9/po4Bx6n687gAu71S5kA3//Bmjto5G/l+vySGKQ\nJOfTe/TpCfSeiveXVfWmqnowyf6n4r3AoU/F+wPg5cBnN/hT8XYlOYPeIes3gfdA76mCA/bfhlZV\nLyS5BPgcMAfsrqq/mfKwZtEm4PZuUuHhwM1VdVeS+4HbkmwH9gJvm94Qpy/JrcDZwAndEz5/C/hd\n+uyj1fheelsOSVLTRj/dJEkawCQhSWoySUiSmkwSkqQmk4QkqckkIUlqMklIkppMEpKkJpOENGZJ\n/mV319wjkxyd5Gvdff6lmecV19IEJLkKeBm9W7s8WlW7pjwkaUlMEtIEJPlH9G709z3gLO9rpbXC\n003SZJwAHA38Y3pHE9Ka4JGENAFJ7gBuAV4FnFRVl055SNKSbLhbhUuTluSXgOeq6pPdQ+j/NMm2\nqrpnykOThvJIQpLUZE1CktRkkpAkNZkkJElNJglJUpNJQpLUZJKQJDWZJCRJTSYJSVLT/we2csUD\ncixwRgAAAABJRU5ErkJggg==\n", + "text": [ + "" + ] + } + ], + "prompt_number": 28 + }, + { + "cell_type": "code", + "collapsed": false, + "input": [ + "Bx" + ], + "language": "python", + "metadata": {}, + "outputs": [ + { + "metadata": {}, + "output_type": "pyout", + "prompt_number": 10, + "text": [ + "array([ 8.24267786e-15 -2.83223098e-15j,\n", + " -6.93872399e-16 -1.71019341e-15j, -6.93872399e-16 -1.71019341e-15j])" + ] + } + ], + "prompt_number": 10 + }, + { + "cell_type": "code", + "collapsed": false, + "input": [], + "language": "python", + "metadata": {}, + "outputs": [] + } + ], + "metadata": {} + } + ] +} \ No newline at end of file diff --git a/simpegEM/Analytics/FEM.py b/simpegEM/Analytics/FEM.py index 78c63ac4..a0fdf42e 100644 --- a/simpegEM/Analytics/FEM.py +++ b/simpegEM/Analytics/FEM.py @@ -36,11 +36,17 @@ def hzAnalyticDipoleF(r, freq, sigma, secondary=True): return hz -def AnalyticDipoleH(x,y,z,sig,f,xs=0.,ys=0.,zs=0.,m=1.,orientation='X'): +def AnalyticMagDipoleWholeSpace(x,y,z,sig,f,xs=0.,ys=0.,zs=0.,m=1.,orientation='X'): """ Analytical solution for a dipole in a whole-space. Equation 2.57 of Ward and Hohmann + + TODOs: + - set it up to instead take a mesh & survey + - add E-fields + - handle multiple frequencies + - add divide by zero safety """ dx = x-xs @@ -69,17 +75,8 @@ def AnalyticDipoleH(x,y,z,sig,f,xs=0.,ys=0.,zs=0.,m=1.,orientation='X'): Hy = front*( (dy*dz/r**2.) * mid ) Hz = front*( (dz/r)**2. * mid + (kr**2. - 1j*kr - 1.) ) - return Hx, Hy, Hz - - -if __name__ == '__main__': - - x = np.arange(-100.,102.,2.) - y = 50. - z = 0. - - sig = 1. - f = 1. - - Hx, Hy, Hz = AnalyticDipoleH(x,y,z,sig,f) + Bx = mu_0*Hx + By = mu_0*Hy + Bz = mu_0*Hz + return Bx, By, Bz