{ "metadata": { "name": "", "signature": "sha256:86da8b82c45841d020bbb210437247025be78a7d463316f2db1b3104fce74c8b" }, "nbformat": 3, "nbformat_minor": 0, "worksheets": [ { "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "An incomplete example showing attempts to fit the RADEX parameter space with low-order polynomials" ] }, { "cell_type": "code", "collapsed": false, "input": [ "%matplotlib inline\n", "import pylab as pl" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 1 }, { "cell_type": "code", "collapsed": false, "input": [ "import os\n", "if not os.path.exists('oh2co-h2.dat'):\n", " import urllib\n", " urllib.urlretrieve('http://home.strw.leidenuniv.nl/~moldata/datafiles/oh2co-h2.dat')" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 2 }, { "cell_type": "code", "collapsed": false, "input": [ "import pyradex\n", "import numpy as np\n", "tau = np.empty([20,4])" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 3 }, { "cell_type": "code", "collapsed": false, "input": [ "abundance = 10**-8.5\n", "R = pyradex.Radex(species='oh2co-h2', abundance=abundance, column=1e14, temperature=20)\n", "R.run_radex()" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "stream", "stream": "stderr", "text": [ "/Users/adam/repos/astroquery/astroquery/lamda/__init__.py:21: UserWarning: Experimental: LAMDA has not yet been refactored to have its API match the rest of astroquery.\n", " warnings.warn(\"Experimental: LAMDA has not yet been refactored to have its API match the rest of astroquery.\")\n", "WARNING" ] }, { "output_type": "stream", "stream": "stderr", "text": [ ": Assuming the density is n(H_2). [pyradex.core]\n", "WARNING:astropy:Assuming the density is n(H_2).\n" ] }, { "output_type": "stream", "stream": "stderr", "text": [ "/Users/adam/repos/pyradex/pyradex/core.py:483: UserWarning: Using a default ortho-to-para ratio (which will only affect species for which independent ortho & para collision rates are given)\n", " warnings.warn(\"Using a default ortho-to-para ratio (which \"\n", "WARNING: Assuming the density is n(H_2). [pyradex.core]\n", "WARNING:astropy:Assuming the density is n(H_2).\n" ] }, { "output_type": "stream", "stream": "stderr", "text": [ "/Users/adam/repos/pyradex/pyradex/core.py:800: RuntimeWarning: invalid value encountered in divide\n", " frac_level_diff = level_diff/self.level_population\n" ] }, { "metadata": {}, "output_type": "pyout", "prompt_number": 4, "text": [ "88" ] } ], "prompt_number": 4 }, { "cell_type": "code", "collapsed": false, "input": [ "fortho = 1e-3\n", "densities = np.logspace(0,8,20)\n", "for jj,dd in enumerate(densities):\n", " R.density = {'oH2':dd*fortho,'pH2':dd*(1-fortho)}\n", " R.abundance = abundance # reset column to the appropriate value\n", " R.run_radex(reuse_last=False, reload_molfile=True)\n", " tau[jj,:] = R.tau[:4]" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 5 }, { "cell_type": "code", "collapsed": false, "input": [ "pl.loglog(densities,tau,'o')" ], "language": "python", "metadata": {}, "outputs": [ { "metadata": {}, "output_type": "pyout", "prompt_number": 6, "text": [ "[,\n", " ,\n", " ,\n", " ]" ] }, { "metadata": {}, "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAX8AAAEHCAYAAABGNUbLAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3X+Q1NWZ7/H3kxl+5Jq4It69JGSyaCdZzca60VVIblYZ\nEegxUxHjza6LJptrBEw2zJhUEk34NTMiKqVWZJq9iThA4i8wqY0RtwMMJY6wmi0xa66aiAWzmiAI\nGn9EJQEZPfeP7oGeme6Z/nZ/u7+nuz+vKqvoM92nH2fg6e+c73meY845RESktrwn6gBERKT8lPxF\nRGqQkr+ISA1S8hcRqUFK/iIiNUjJX0SkBin5i4jUICV/EZEaVNLkb2aNZrbdzH5gZlNL+V4iIpK/\nUl/5vwu8CYwBXijxe4mISJ4CJ38zW2NmB8zsqUHjTWa208x2mdk16eHtzrnPAt8FOkKIV0REQlDI\nlf9aoClzwMzqgJXp8Y8Ds83sNHescdDrpK7+RUTEA/VBX+Cc225mkwYNTwZ2O+eeBzCz9cAsMzsV\niAMnAImiIhURkdAETv45TAT2ZDx+AZjinLsRuG+4F5qZ2oqKiBTAOWeFvjasG75FJXDnXKD/2tra\nQn1Ntq8NHgv7PSs5zpGeX8txluNnrjjDjbNS/60XK6zkvxdoyHjcQIDdPe3t7fT09OT9Zo2NjXk/\nN5/XZPtaIe8RxhyVEOdIz6/lOMvxMy/0fYp9fbXGWWn/1nt6emhvbw/8fkME/bRKf+JMAp7KeFwP\n9KbHRwO/Bk7Lcy5XCdra2qIOIS+KM1yKM1yVEGclxOicc+ncWVAOd84VtNVzHfAo8DEz22Nmlzvn\n+oD5wGbgt8C9zrlnivxc8koYVwfloDjDpTjDVQlxVkKMYTAXwtpRUQGYuba2NhobG2vmmy4iUqie\nnh56enro6OjAFXHD14vkH3UMIiKVxsyKSv5q7CYiUoOU/EVEapAXyT/oVk8RkVoV1lZPrfmLiFQg\nrfmLiEhgSv4iIjVIyV9EpAYp+YuI1CAvkr92+4iI5Ee7fUREapj3u33M7Dgz22FmzaV+LxERyU9Y\nJ3kN52rg3jK8j4hI1Utu3Urnz39e9DyFtHReY2YHzOypQeNNZrbTzHaZ2TXpsRmkWjy/XHSkIiI1\nLrl1K3NWd9F98cVFzxV4zd/MzgHeAu5wzp2eHqsDngWmkzrVawcwG7gMOA74OPBn4PODF/i15i8i\nkp8zL5vNE3OvTD0477yi1vwDL/s457ab2aRBw5OB3c655wHMbD0wyzm3KP34y8DLyvIiIoV77pX9\noc0V1pr/RGBPxuMXgCn9D5xzPx7uxZnblnSoi4jIUD09PRx6dhf86EehzBdW8i/qij6Uw4hFRKpY\nY2Mjp50+iSfe/B20tMGPh72mHlFYWz33Ag0ZjxtIXf3nRUVeIiIjW9ryPcb917Mw738XPVdBRV7p\nNf8HMm741pO64Xs+sA94DJidzyHuuuErIpK/5JYkiXUJNq/dXN4zfM1sHTAVGA+8BCxxzq01swuA\nW4E6YLVz7oY851PyFxEJqNgK30J2+8zOMb4R2FhIEO3t7brRKyKSh56enlCWydXbR0SkTNpvWc7K\n5Ab6Ro+i/u0jzG++kPZvXVPQXGW/8hcRkeDab1nOsocepG/JsqNjy26+PvW1Aj8AiqErfxGRMjhp\n2md4JSPx9xu/dCF/ePCRwPN539UzH9rqKSLVrm/0qOzjo7KP56J+/iIiFURX/iIiNWh+84XUp9f4\n+9XftIz5n70wknh05S8iUibttyxn5S820DdqFPVHjjD/s9Ht9lHyFxGpQFr2ERGRwJT8RURG0H7L\nck6a9hlOaGrkpGmfof2W5VGHVDQvirzU3kFEfOVbcZbaO4iIlEHYWzTD4vWav5mdamY/MLOfmNkV\npXwvEZFSCKs4yzclTf7OuZ3Oua8B/wjES/leIiKlUP/2kezjR7KPV4rAyd/M1pjZATN7atB4k5nt\nNLNdZnZNxvjngCSwvvhwRUTKy7firLAUcpjLOcBbwB0ZJ3nVkTrJazqpIx13MOgkLzO73zk3K8t8\nWvMXEa+FWZwVlkiKvLIc4/hpoM0515R+/N30U38JXAyMBZ5xzt2aZS7X1tZ29LF2/YiIDDV4l09H\nR4cXyf8LQNw5Nzf9+IvAFOdcSx5z6cpfRCQgXw5zUfYWEe8kt26l8+c/57AZY5yj9aKLaJ42Leqw\nvBBW8t8LNGQ8bgBeyPfFKvISkbAlt25lzuou9s+dd3TsydtX0QUV/QEQVpFXWFs9Hwc+amaTzGw0\ncAmwId8X9yd/EZGwLF59+4DED7B/7jyWrOmKKKJwvOfgQfp++cvi5wn6AjNbBzwKfMzM9pjZ5c65\nPmA+sBn4LXBv5k6fkegkLxEJ23Ov7M8+/ofs45VgWzLJqnnzqO/uLnquwMs+zrnZOcY3AhsLCSKM\nI8lERDLZ24eyfyHXeAXo7uzkrn37AOgoci519RSRqjTpfXWQGJQiOzuY9L7KTXt/3HcgtLm86Oop\nIhK2pS3fY84Nrey/+Uqofy/0/ZkJB/tYuqAz6tAKtuvFN0Kby4vkr90+IhK25hnNdAGJdQkOvXuI\nse+ZQMvsFppnNEcdWsH+OOFTnPfKYaayr+i51NJZRKRCxOOLeLT705xKgsfZ7G9LZxERCU9r60z+\nR+xRHmdT0XMp+YuIl6rx6MRiNTefy4oVceLxxUXPpWUfEfHO0aMTv73g6Fj9zdez8LzzI++m6YtI\nunqGSclfRAbz9ejEYm1LJunu7KT+8GH6xoxhZmsr5zYXdgPal8ZuIiKhqcajE7clk/xszjxu3X9s\np843nnwaulYV/AFQDC/W/NXeQUQyVePRiXcv7hiQ+AFu3b+Pu5dcG2ienp6eULoiaNlHRLyTdc3/\npmUsnDa9Ytf8Pz/ur7jv9d8PGb943If52au/Czyf18s+ZjYLaAaOB1Y757aU8v1EpDr0J/iVSxd6\ndXRiMQ5a9nR7kGiWsspy5W9mJwA3O+fmZPmarvxFpOpNPfMSJjzxK+6l9+jYPxDjwBln8fB/rg88\nX7FX/oW0dF5jZgfM7KlB401mttPMdpnZ4I/nRcDKQoMUEal0Vy/9Ok9MOJOzidPIVM4mzq8nnMHV\nS/85kngCX/mb2TnAW8AdGWf41gHPAtNJneq1A5gN7ARuBLqdcw/mmE9X/iJVREcn5pZMbiOR2MKh\nQ3WMHfsOLS0zaG4+t6C5Itnnn+UA908Dbc65pvTj76afehD4MqkPg187527LMpeSv0iVyHZ04oTb\nV9F1xRx9AITMlxu+E4E9GY9fAKY451qAxEgvzty2pO6eIpUrdXTilQPG9s+dx5KuVRWd/MMszipU\nWGf39gsr+Rd16a6TvESqQ7Uenbj5qqtY1nvsRu3C9J/L+QEw+MK4o6O4s7zCKvLaCzRkPG4gdfWf\nFxV5iVSHaj06MTPxAyzr7WVLYsRFjZIIq8grrOT/OPBRM5tkZqOBS4AN+b64/zAXEalstXR04ut7\no/ltprGxMZTkH3jZx8zWAVOB8Wa2B1jinFtrZvOBzUAdqYKuZ4qOTkQqSi0dnbj7xTfLHEm4Aid/\n59zsHOMbgY2FBKFjHEWqQ7UenXjJK+8ZUpz1+oTJkcQT1o1f9fYRERlG5tGJx3GIg4xlJy18Jv4f\nbNq0NLK4fNnqKSJVQAVaQ7W2zqS3dzOP9x47OjEWW0BLS1OEURVPyV9EgOwFWk/evoouqOkPgP4K\n3ERicUZlblPBlbm+8GLZp62tTWv+IhE787LZPDGoQAvgzK5V/OqueyKIqHg+FGeFrX/Nv6OjQ8c4\nikjxxjWdx+vfbRs6fuO1vLppawQRFSdrcVYsRnzFior/AIAIunqKSHWqtgIt34qzfKPkLyJA9RVo\n+Vac5Rvd8BURoPoKtKq1OCssSv4iVSCMLZrVVqDlW3GWb7xI/qrwFSlcmFs0m2c0V2yyH+z4iafw\ni99cxtmDi7M+9B9Rh1YUVfiKCFCdWzTDkExu46qrNtPbu+zoWCy2gBUrKn+PPnhe4WtmJwMLgb9w\nzv19Kd9LpFZVYw/9MFRrcVZYSpr8nXPPAXPM7KelfB+RShXGWn21bdEMU3PzuUr2ORTS0nkN0Ay8\n1H+Gb3q8CbiVVEvnLufc8tCiFKlCYa3VT3pfHa8lOqAlo0Crgrdo9qvG6lyfBF7zN7NzgLeAOzIO\ncK8DngWmkzrVawcwu7+nv5n9NNeyj9b8pVK137KclckN9I0eRf3bR5jffCHt37om79eHtVaf3JJM\nbdH8b/UDtmh2Leis2Ju325JJfjZnHrfu33d07BsTPsjFXav0AZBW9jV/59x2M5s0aHgysNs593w6\nqPXALDM7AFwPfNLMrtFvA1It2m9ZzrKHHqRvybGbictuvj71tTw/AMJaq6+2LZoAdy/u4LaMxA9w\n6/59XLnkWiX/kIS15j8R2JPx+AVginPuVeCrI70480gybfmUSrAyuWFA4gfo+/YCVi5dmHfyD3Ot\nvpq2aAK89Fz26tyXn6vdm9hhbfHsF1byL2rdJozzKEXyldySpPOeTg67w4yxMbRe2ho4cfaNHpV9\nfFT28Wyqda0+DActe2o6SP7f32oz+MK4o6Mj95PzEFby3ws0ZDxuIHX1nxcVeUm5JLckuepfrqL3\njGNVn73/kvpzkA+A+rePZB8/kn08m2prpxCmw5PO4pLXbEh17qFJZ0UYlR8iLfJKr/k/kHHDt57U\nDd/zgX3AY2Tc8B1hLt3wlbKJXx6nu247vPhhGDUWjhyCD/ye+LvnsmnNppEnSDu65v/tBUfH6m9a\nxsJp0wPd9E1uSWas1Y+t+LX6sCST2/jmnJX8xf43jlbn/nHC+/l+V4u2bqYVe8O3kN0+64CpwHjg\nJWCJc26tmV3Asa2eq51zN+Q5n5K/lM0nPjeZ34w+buBSS6KDv3n7IE8/8FigudpvWc7KX2ygb9Qo\n6o8cYf5ng+32keElk9tIJLZkFGjNUOLPUPbkHzad5CXldNK0z/DKoBu1AOOXLuQPDz4SQUQiwegk\nL5ECnP5PX+Lpr1wxZPwTa1bz1B13RhBR9VFxVnl43dtHxDcfPGEcT2cZnzjuxLLHUo2yHp2Y/rM+\nAPyiPWVSMZJbksQvj9P4fxqJXx4nuSUZeI7Wiy4idvfdA8Zid91Fy6xZYYVZ03R0YuXw4spfWz1l\nJGFt0ezvmZO47z4OAWOBlksvDdxMTbLT0Ymlp37+UlPil8fpntQ9dPx38UBbNKW0mk46hU2vPDdk\n/ILxp7DxD71ZXiGF0pq/1ITD7jDsee+Q/fmH3lXbYp/o6MTKoeQvFeGNV/tg1BT4zsD9+W8EqKiV\n0qvWoxOrkZK/VIb3N8Dg9sctbVjXqsBTJZPb6Ozs5vDhesaM6aO1dWZBxUNhzVNNWltn0tu7mcd7\njy3FxWILaGlpijAqyUbJXyrC8Sf9Zdbx94//74HmyXaua2/vQoBAiTuseaqNjk6sHEr+UhHG5NgU\nMDbgPJ2d3QMSNkBv7zISicWBElRY8/gmjAItHZ1YGbxI/trqKSNpvegieu++m97LLjs6FrvrLlou\nvTTQPIcPZ/8rf+hQXSTz+EQFWpUhrK2e3iR/qW7F9tAPa3/+mDF9WcfHjn0nknl8kqtAa3EioeTv\nkf4LZV/6+WdlZscB/xc4DPQ45/I/mFSqRpgFWsUWY6VuSC4csGRTyA3JsObxqQ+OCrRqS6mv/C8G\nfuKcS6bP9VXyr0Gd93QOSPwAvWf0kliXCJT8w9hdE9YNyebmc3lmxy95cGWMMX3vcLi+jvO/OC/Q\nPFkPKX/yaYjokPJdL76RdXz3i2+WORIph8DJ38zWAM3AS/2HuaTHmzjWz78rfVj7ROD/pZ9Sub8P\nS1HCKNAKc3dNGDcktyWTvHbX7Wx85b+Oji2863a2nf2JvBN3mIeUh/EbhAq0akshV/5rgQRwR/+A\nmdUBK4HppI503GFmG0gd5dgAPImayNWsMAq0fNtdE8b6eFiHlId1o1YFWrUlcPJ3zm1PH+OYaTKw\n2zn3PEB6iWcW0AmsNLNmYEOuOTNv+GrXTxUKoUDLt9019YcPZx2vO5T/bzNhHVIe1o1aFWj5Laxd\nPv3CWvOfCOzJePwCMMU59yfgKyO9WLt9qlsYBVq+7a7pGzMm6/g7Y/OvPAjrkPIwPohABVq+G3xh\n7MtuH7XllJzCKNAKa3dNWGa2trKwt3fAFfeCWIymlpa857h66df55pyVnL3/IwMPKV/6z4FiefGN\nP2Uff/PPgeYBFWjVkrCS/15Sa/v9Gkhd/edFRV7VLYwCLd+uSvuXUxYnEtQdOsQ7Y8fS1NISaJml\nuflc6OLoIeXjx75DewGHlO/kZC7h1SG/QRxwJweaRypDpP3802v+D/Tv9jGzeuBZ4HxgH/AYMNs5\n90wec6mffw1Ibt1K4v77jxVozZqlA1RC0tjYzq8ePptTB92o/dupO+jpaY86PCmRsvfzN7N1wFRg\nvJntAZY459aa2XxgM6mtnqvzSfz9dOXvr2Irc/uFUaAl2Y0Z08dbNPM4A38uY8dql0410kleUnLZ\nKnNjT8RY8fUVgT8AfGp/7FNVbRiy1UDEYgtYsUI3a6uZTvKSkum8p5Pek/bBY399tDir9wO/L6gy\n15f2x9XYvMy3+yFSGXTlLzl94nOT+c3o46BlYHHW37x9kKcfeCzveeLxRXR3X5dlfDGbNi0NI9S8\nLYrHua576FnAi+Nxlm7SWcBSOYq98vei6ra9vT3U4gUJx/6DowYmfoCWNvb/KVgRkk8FWmHtiReJ\nSk9PTyi1Ud4kf93s9c8HPnRK9vGJ2cdz8alAK4ziLJEoNTY2Vk/yFz998IRxWccnjjsx0DytrTOJ\nxRYOGEsVaM0oOLZCzWxtZWEsNmBsQSzGjADFWSLVQDd8JaewTs/y6YZkGMVZItVAN3xlWCrOEvFT\nsTd8lfyrWFgFWiLiH+3zl6zCOjpRRKqTFzd8tdUzfAMKtJ74n/DYX9N70j4S6xJRhyYiRQhrq6eW\nfapUWAVaIuInr4u8zOxkM+sys5+W8n1kqLAKtJLJbcTji2hsbCceX0QyuS3EKEUkKiVd83fOPQfM\nUfIvvw986BReyTYeoEDLp548IhKuvK78zWyNmR0ws6cGjTeZ2U4z22Vm15QmRClEGAVauQ9N31JU\nbIXalkyyKB6nvbGRRfE425LJSOIQqQb5XvmvBRLAHf0DZlYHrASmkzrJa4eZbQDOAs4EbnLO7Qs3\nXMlXGAVaPvXkqcZunCJRyiv5O+e2p0/vyjQZ2O2cex7AzNYDs5xzNwJ3psdOBK4HPmlm1zjnlocU\nt4ygvxArcd99xwq0Lr00UIGWTz15ujs7ByR+gGW9vSxOJJT8RQpQzJr/RGBPxuMXgCmZT3DOvQp8\ndaSJMrct6USv8BR7epZPh6arG6fUurBO8OpXTPIPbX9mGHtWq4kvlbk+9eRRN06pdYMvjDs6Ooqa\nr5jkvxdoyHjcQOrqPzCd4XtMckuSOTe0sv+9o46envXkDa10EU1lbnPzuV7s7JnZ2srC3t4BSz8L\nYjGa1I1TakzZz/BNr/k/4Jw7Pf24HngWOB/YBzwGzA5ycHt6HhV5ZTjzwr/jiVGjhhRnnXHkCP+5\n4d+jC8wD25JJtmR045yhbpxSw8rS2M3M1gFTgfHAS8AS59xaM7sAuBWoA1Y7524IHICS/wAnTvs0\nry0Z+m0cd+0CXt36aAQRiYiPytLYzTk3O8f4RmBjoW/eT8s+x7jROdawc42LSE0p+7JPqejKf6Az\nL5vNE3OvHDretYpf3XVPBBGJiI+87u0jwS29Yi4Tbl81YGzCqtu49itzAs+lvjwikov6+Xumedo0\nuhhUnDVnbuD9+urLIyLD8WLZp62tTWv+IYvHF9HdfV2W8cVs2rS07PFsSybp7uyk/vBh+saMYWZr\nq3bqiBSgf82/o6Oj8k/yqpYiL1+Ks0B9eUSqVf+FcrFFXlrzD0l/cVb3S7/j4ddfp/ul3zHnhlaS\nW6LpPFkJfXm2JHSqmEhUlPxDsjhxA/vHfRi+80P4xq3wnR+yf9yHWZwIXPoQitbWmcRiCweMpfry\nzCh7LOrLI+IfL5Z9qsHzb70DSwatsbe08fy1CyKJR315RGQ4Sv4h8bE4S315RCQXL5J/NVT4njx+\nAk9kGz9pQtlj8U3/Td3FGX15mtSXR6QgqvD1THLrVuas7mL/3HlHxyasuo2uAvboi4iMpCyN3Yph\nZrOAZuB4Us3ftgz6elUkf0h9ACTuv/9YcdasWUr8IlIS3if/o29kdgJws3NuzqDxqkn+IiLlUrbe\nPma2xswOmNlTg8abzGynme0ys2uGmWIRqQPfRUQkYkEOczkHeAu4I+NAlzpSB7pMJ3Wy1w5gNnAW\ncCZwE/AicCPQ7Zx7MMu8Xlz5+1Sdm0xuo7Ozm8OH6xkzpo/W1pmR7dpRWwYRP5Wlnz+Ac257+jSv\nTJOB3c6559PBrAdmOeduBO5Mj7WSOu3reDP7iHPutkKDLRWfjk70qSGb2jKIVK9Aa/5ZjnL8AhB3\nzs1NP/4iMMU5l/cG7v7Gbv2i2PLp09GJPjVkWxSPc11395DxxfE4SzdtKmssIrVu8BbPqBu7hbJe\nE3VjN5+qc31qyKa2DCL+GHxhHHVjt71AQ8bjBuCFoJO0t7eHUrRQKJ+qc31qyKa2DCL+6enpCeWC\nudjk/zjwUTObZGajgUuADUEn6a/wjcrJ47NX4UZRnetTQ7aZra0sjMUGjC2IxZihtgwikWlsbAwl\n+ee97GNm64CpwHgz2wMscc6tNbP5wGagjlQR1zNFR1VmS6+Yy5zbVw2pzr12ztyyx+JTQza1ZRCp\nXl60d/DhJC9V54pIJQjrJC8vkn/UMYiIVJqy7fP3lU/FWSIilaKik79PxVkiIpXEi2WfQtf8fSrO\nEhEpB635AydO+zSvLRl6Ru64axfw6tZHiw0tMJ968oD68ohUs5pe8/epOMunnjygvjwiMrxii7wi\n5VNxVmdn94DED9Dbu4xEYkuOV5RWd2fngMQPsKy3ly2JRCTxiIhfKjr5L71iLhNuXzVgbMKq27j2\nK3NyvKJ0fOrJA+rLIyLDq+hln+Zp0+gCEvfdd6w4K6Izc33qyQPqyyMiw6vo5A+pDwAfKnFbW2fS\n27twwNJPqidPUyTxzGxtZWFv74ClnwWxGE3qyyMieJL8+xu7RdneoVg+9eQB9eURqVaD+/oXyout\nnuPP+1/Mb76Q9m8NdwSwiIj0K9sB7oUws1PN7Adm9hMzuyLX815ZsoxlDz1I+y3LSxmOiIikleXK\n38zeA6x3zv1Dlq85HnoIgPFLF/KHBx8peTwiIpWuLFf+ZrbGzA6Y2VODxpvMbKeZ7TKzrGs2ZvY5\nIAmsH+l9+kaNyieckkgmtxGPL6KxsZ14fBHJ5LbIYhERKbV8b/iuBRLAHf0DZlYHrASmkzrOcYeZ\nbQDOAs4EbnLO7XPOPQA8YGb3Az8bNpgjR4L/H4TAx+pctWUQkVLKK/k757ab2aRBw5OB3c655wHM\nbD0wyzl3I3BnemwqcDGpLfgP5XyDH/0Ie2Qbnzj5FHp6esq+6yd3de7isid/tWUQkWzC2uXTr5it\nnhOBPRmPXwCmZD7BOfcw8PBIE43fs4v5X70yst0+PlXn5mrLsDiRUPIXqWGDt8N3dHQUNV8xyT+0\nO8Xzz5lB499OGfmJJeJTda7aMojIcML6DaCYrZ57gYaMxw2krv4D6y/yikpr60xisYUDxlLVuTPK\nHovaMojIcBobG2lvby96nmKu/B8HPpq+F7APuASYXXREEfCpOldtGUSkHPLa529m64CpwHjgJWCJ\nc26tmV0A3ArUAaudc0NPVhl57oJP8qpW25JJtmS0ZZihtgwikqaTvEREapjX7R1ERMRPXnT1LIZv\n5+aKiFQCL5J/oS2dVZkrIrWmqlo6FxpDPL6I7u7rsowvZtOmpcWGFkjWytxYjPiKFfoAEJHQ1fSa\nfyVU5urAdBHxUUUnf1XmiogUpqKTvypzRUQK48UN30KpMldEpDAVfcPXN6rMFZFyKfaGrxfJX+0d\nRETyo/YOIiI1zPutnmZ2nJntMLOc6x86M1dEpLzKccP3auDe4Z7Q3X1dpJW5oOpcEakteSV/M1sD\nNAMvOedOzxhv4lhL5y7n3PJBr5sB/JbUGb7DiurMXNC5uSJSe/Jd9lkLNGUOmFkdsDI9/nFgtpmd\nZmZfMrPvm9kHSZ0B8CngUmCumQ27PhVFZS6oOldEak9eV/7Oue3pE7syTQZ2O+eeBzCz9cAs59yN\nwJ3p5yxKf+3LwMsj3dmNojIXVJ0rIrWnmDX/icCejMcvAFlPYXfO/Xj4qdoZN+5BTjzxI/T09JR9\ny6eqc0XEd2F18+xXTPIPbX9mPP4OLS3LIrvZq+pcEfHd4Fqojo6OouYrJvnvBRoyHjeQuvoP7FOf\nquO4494tIpTi9N/UXZxRnduk6lwR8VDZ+/mn1/wf6N/tY2b1wLPA+cA+4DFgtnPumUABqMhLRCSw\nshR5mdk64FHgY2a2x8wud871AfOBzaS2c94bNPGLiEg08t3tMzvH+EZgY7FBFHqMI6g4S0Rqi45x\nREcnikjt8r63TympOEtEpDAVnfxVnCUiUhgvkn97e3tBa1gqzhKRWtPT00N7e3vR81Tdmv+CWIwm\nrfmLSJWripO8iolBRyeKSC2q+eQvIlKLanq3j4iIFEbJX0SkBnmR/M+NxVhxww1RhyEi4r3q2u2D\nKnNFRILwes3fzBrNbLuZ/cDMpg73XFXmioiUT6mXfd4F3gTGkEevf58rc8M8QaeUFGe4FGe4KiHO\nSogxDPm2dF5jZgfM7KlB401mttPMdpnZNVleut0591ngu8CIx874XJlbKX8hFGe4FGe4KiHOSogx\nDPle+a8FmjIHzKwOWJke/zgw28xOM7Mvmdn3zeyDGRv4Xyd19Z/TgliMGXkem1jID2e412T7Whh/\nAao1zpGeX8txluNnXuj7FPv6ao2zlv6tZ8or+TvntgOvDRqeDOx2zj3vnDsCrAdmOefudM590zm3\nz8w+b2bqbw9mAAADzElEQVQ/BO4Aci7oL47HA7VkqJRvdLXGqeRf+HsGfU0lJ9VC36fY11fCv6GR\nXlOO5F/MMY5fAOLOubnpx18EpjjnAp16bmYq7xURKUAxu32KOcA9lKRdTPAiIlKYYnb77AUaMh43\nkMeOHhERiV4xyf9x4KNmNsnMRgOXABvCCUtEREop362e64BHgY+Z2R4zu9w51wfMBzYDvwXudc49\nU7pQRUQkLJG3dxARkfLzorGbiIiUl3fJ38yOM7Mfm9kqM7s06nhyMbOTzazLzH4adSzDMbNZ6e/l\nejObEXU8uZjZqekeUD8xsyuijieX9N/PHWbmbQfCID21omQpy8ys08z+Kep4cjGzv0t/L283s0ei\njicXM/uQmf3MzFbn6LgwgHfJH7gY+Ilzbh5wYdTB5OKce845NyfqOEbinLs//b38Kqmb8l5yzu10\nzn0N+EcgHnU8w7gauDfqIEYQqKdWhC4CJgJv43Gczrl/T//d/DfgRxGHM5zTgX91zl0BnDHSk8uS\n/AP2BpoI7En/+Z1yxFdgnJEpMM5FpNpxlE3QOM3sc0CSVLW4dzGmf3P6LfByueIrJE4C9tSKMM6P\nAY84574NfM3jOPtdCtxTvigDx/koMM/MHgQ2jTi5c67k/wHnkPokeipjrA7YDUwCRgG/Bk4Dvgg0\np5+zrhzxFRJnxtd/Ws4YC/h+GrAcON/nOAe97n4fYwSuA75Paofbz0lvmPAtzoyvjy7338+A38/L\ngL9PP+deX+NMf+3DwKpyxljA9/MbwDnp54z4cy+mwjdvzrnt6fYQmY72BgIws/XALKATWJleUy1r\n3UCQOM3sAHA98Ekzu8Y5t9zHOIHpwPnA8Wb2EefcbT7GaWZ/SWrJbyzwkI8xOucWpR9/GXjZpf+V\n+RanmZ1KaunsBIbpqVUKAf9urgASZnYO0FO+KAPH+QzwFWBNGUMEAsf5C2BJ+l7pcyPNXZbkn0Pm\n8g6k1vymOOf+ROob7Ytccb5Kah3dF7nibKHMCWAEueJ8GHg4mpCGyBpj/wPn3I/LHlF2ub6XNwL3\nRRNSVrni/DPg032znD9351x7FAHlkOv7+STwhXwnifKGb6UUGCjOcFVCnJUQIyjOsNVUnFEm/0rp\nDaQ4w1UJcVZCjKA4w1ZTcUaZ/CulN5DiDFclxFkJMYLiDFttxVmmO9brgH3AYVJrVZenxy8AniV1\n5/p75b6TrjgVZyXEqDgVZyniVG8fEZEa5GOFr4iIlJiSv4hIDVLyFxGpQUr+IiI1SMlfRKQGKfmL\niNQgJX8RkRqk5C8iUoP+P/O787v1lBsuAAAAAElFTkSuQmCC\n", "text": [ "" ] } ], "prompt_number": 6 }, { "cell_type": "code", "collapsed": false, "input": [ "fits_7 = [np.polyfit(np.log10(densities),np.log10(t),7) for t in tau.T]\n", "fits_11 = [np.polyfit(np.log10(densities),np.log10(t),11) for t in tau.T]\n", "fits_7" ], "language": "python", "metadata": {}, "outputs": [ { "metadata": {}, "output_type": "pyout", "prompt_number": 7, "text": [ "[array([ 4.86015567e-04, -1.45822371e-02, 1.69408821e-01,\n", " -9.47638215e-01, 2.57322263e+00, -3.05652747e+00,\n", " 2.30492326e+00, -3.97239753e+00]),\n", " array([ 2.39146864e-04, -7.13724406e-03, 8.25362832e-02,\n", " -4.60721266e-01, 1.25940091e+00, -1.55001397e+00,\n", " 1.66597043e+00, -2.90593168e+00]),\n", " array([ 5.87632284e-05, -2.77338874e-03, 4.41794386e-02,\n", " -3.16452339e-01, 1.05884924e+00, -1.50480111e+00,\n", " 1.74907731e+00, -4.85548112e+00]),\n", " array([ 1.47711558e-04, -4.58305035e-03, 5.49716696e-02,\n", " -3.17766208e-01, 9.00971669e-01, -1.16643424e+00,\n", " 1.52717654e+00, -2.93472028e+00])]" ] } ], "prompt_number": 7 }, { "cell_type": "code", "collapsed": false, "input": [ "pl.loglog(densities,tau[:,0],'o')\n", "pl.loglog(densities,10**np.polyval(fits_7[0],np.log10(densities)))\n", "pl.loglog(densities,10**np.polyval(fits_11[0],np.log10(densities)))" ], "language": "python", "metadata": {}, "outputs": [ { "metadata": {}, "output_type": "pyout", "prompt_number": 8, "text": [ "[]" ] }, { "metadata": {}, "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAX8AAAEHCAYAAABGNUbLAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XmcjfX7x/HXZTBjydZmLZKK0pfsfSu7GSQSaZAQLd8I\nv3bRjKJosc1UlIjKjLKODmMGsxAylmQPobGE7GLGLJ/fH2doaIZZzpz7LNfz8TgPc+65z32/5xjX\nuX3uzyLGGJRSSnmXQlYHUEop5Xxa/JVSygtp8VdKKS+kxV8ppbyQFn+llPJCWvyVUsoLafFXSikv\npMVfKaW8UIEWfxG5R0Q+F5HvReTZgjyXUkqpnBNnjPAVkUJAuDHmyQI/mVJKqevK9ZW/iEwVkSMi\nsvmq7QEiskNEdonIG5m2dwBsQHj+4yqllHKEXF/5i8jDwDlghjGmdsY2H2An0Ao4CCQAgcaY7Zle\nt8AY09FRwZVSSuVd4dy+wBizQkSqXrW5IbDbGLMPQETCgY4icgvQGfADYvKVVCmllMPkuvhnoxKQ\nmOn5AaCRMSYOiLvWC0VEpxVVSqk8MMZIXl/rqN4++SrgQUFBxMTEYIzJ0SMoKCjH++bkNVl97+pt\njj6nO+e83v7enNMZf+ea07E53e3fekxMDEFBQfku2o668j8IVMn0vAr2q/8cCQ4OztXJmjVrlqv9\nr/earL6Xl3M44hjukPN6+3tzTmf8nef1PPl9vafmdLd/682aNaNZs2aMGDEi1+e8Qm4/rYwxAFWB\nzZmeFwb2ZGwvCvwC1MzhsUxQUJCJiYkxriwoKMjqCDmiOR1LczqWO+R09YwxMTEmKCjI2Mt37uv3\npUdeunqGAauAu0QkUUT6GGNSgQHAEmAbMMtk6ulzPcHBwQ759C1Irp7vEs3pWJrTsdwhp6tnbNas\nWa5bS7LilEFe1wwgYoKCgi7/V0YppVT2YmNjiY2NZcSIEZh83PB1ieJvdQallHI3IpKv4u8SE7sF\nBwcTGxtrdQyllHJ5sbGxntPsY3UGpZRyN3rlr5RSXuSDDyZQvfoj+T6OXvkrpZSbsNniGTRoCXv2\njAI84MpfKaXU9U2cGMWePaN4gHX5PpZLFH9t9lFKqeu7mORD18K9aeHTKN/H0mYfpZRyB0ePsrLO\n/fhwgu63P8e+NZ9qs49SSnmysz/O5UTNqmy4+zyBlZ9i35rQfB9Ti79SSrmqlBT2Pv8kf/foysxX\n/ekX9SefBvXD3394vg/tEsVf2/yVUupK53dsYd99ldm7YiG/Lf2eAW/No3iR4pQokU7jxj75Pr62\n+SullIv5LXQEN735Lou61ePR0CWUKVb2X/vkd5CXFn+llHIRF0+f4NduzSmzYSt/fD6aFk+8mu2+\nHjHCVymlvN1vS2dx8O6KnE46Salff7tm4XcELf5KKWWhtLRUlg7qQLlOgfzx8jO0iNnPLeXvKPDz\nOmoZR6WUUrm097e1HOoaQKWzqSStiKFp3aZOO7dLFP9LK3npYi5KKU9ms8UzcWIUSck+3FEkkpFr\nE/B9rDl3T7FRyNcvR8e4tJhLfukNX6WUcoJLk7LtP/Icwbe3oPfeRF6v2Iru49+kffvcz9KpN3yV\nUsoNTJwYxRmpQUzpGtQ7JDxw4Xdm7l5ESEi0JXm0+CulVAG7kHKBC6UWsPpYf2LPdKfdyd84aioD\nkJSU/wFbeeESbf5KKeWpthzdwsejH2OO7Q/eSBnPtNQBV3zfzy/NklwFeuUvIh1F5AsRCReR1gV5\nLqWUciXGGCatm0TI4CZ8/vVf7Bk6hvjbD1+xT/XqQxk40JrS6JQbviJSBvjYGNMvi+/pDV+llEc5\nceEE/RY8S8uwn3lug1Bk8RK47z5stnhCQqJJSvLBzy+NgQNb5+lmL1gwvYOITAXaA0eNMbUzbQ8A\nxgM+wBRjzJhM3/sY+NYY80sWx9Pir5TyGPH74+n9Qw/CokvT4FhRCtlsUKGCw89jRW+faUDAVSF8\ngNCM7bWAQBGpKXZjgMVZFX6llPIUqempBMUE0e+brqyddxONilajUHx8gRR+R8j1DV9jzAoRqXrV\n5obAbmPMPgARCQc6Aq2AlkApEbnTGDM5X2mVUsoF7T+1nx5ze3D7aWFrWDmKNP8vTJgAPtb05MkJ\nR/X2qQQkZnp+AGhkjBkIhFzvxcHBwZe/1pG+Sil3MnvbbP5n+x8f3RhIr0/mIK+8AoMHg+S5RSZL\njhrZe4mjin++G+216Cul3Mn5lPMMjhzM8r3LWXnLm9z16gcweTJ07lwg57tUIy2d3iGj2WfhpRu+\nItIYCDbGBGQ8fwtIz3zT9xrH0hu+Sim3sunPTQTOCaR+xfp8kVgHv9Efw7x50KiR0zK4yvQO64Aa\nIlJVRIoC3YCInL5Yl3FUSrkDYwwhP4fQ6ptWDP3vm8xYeTN+n30BK1c6rfDHxsZe0VSeV3np6hkG\nNAVuBI4C7xhjpolIW/7p6vmVMeaDHB5Pr/yVUi7NZovn48/ms/kOGym+Zxn3wCj6zvkRjh+HuXOh\nXDmnZ8rvlX9eevsEZrN9MbA4LyF0SmellKuy2eJ5cfhMEh9eCjs6cfPSwTxQpBEHm9xLpSVLwNfX\nqXl0SmellHKCBl17sa5qFCx7n7s3PoiN9nxHD1a1SSdyyUjLcrlKm3++aJu/UsoVzdw8k013zYW5\n3+K/sSJxNGUkwwjiXZKSrZkX07I2f0fTK3+llKsxxjBqxSimbJjCfUseokdkGo34mb5MJY5mAPj7\nDycy8j3LMjq9zb8gaJu/UspVpKSl8PyPz7P58CY2mRfwWzWaaWVqcd+pLVygOHBpNs6A6xypYGib\nv1JKOdippFN0+b4LtQ6lMG7eBXwKF4HJk7HtP+Gw2TgdxemzejqaFn+llCvYf2o/Xaa1ZfRPfrSI\nP4CMGgXPPguFXOLW6L9os49SSuXTukPrCB3uz7JFhlKt2sOWSLjlFqtjZUmbfZRSygGi4qZyccCL\nNDt7IyW/+gZatrQ6Uo54RFdPpZRyutRU4od0pn67ftRt1YuSO353m8LvCC7R7KOUUs6UtvZnDvV4\nDF85x/nYaCo38J6if4lLXPnrIC+llFOcOUPKSy9ypk1Tvmt+I3dvOuB2hV8HeSml1HXYbPFMnBhF\ncpIPLU9t4f8OxhJdsxDR/VsyrscMivoUtTpinnlEbx+llHI0my2eQYOWcHbPy3zFs1Tz2UH7rslU\nbvEk3zz9FeLglbbcjUs0+yillKNNnBhF+T1tWU89NpYpQ91Bp4nb9jl/zank9YUftPgrpTyRMXTY\nk8AcnuC5qk/xTr9oUub+AL8+TVKS6y6q7kza7KOU8ixnz0LfvrQ59iuN6vVm/8PhMD0GjtUCwM8v\nzeKArsElrvy1t49SyiG2bYMGDUgvW5a33mvE4YemwdSfLhd++4RsrS0OmT/a20cppTILC4OXXyZl\n9Pt0LxXFX+f/4oUyrzLt8zUuNSGbo+jEbkop73bxIrz6KixaxJnvptFh+zDKlyzPjE4z8C3s3CUW\nnUmnd1BKea8DB6BpU9i/n4PLF/Dg+hepW74uYU+EeXThdwQt/kop97RsGTRoAJ06sWXyezSZ05be\ndXozzn8chURL2/UUaG8fEakGvA2UNsZ0LchzKaW8RHo6jB4NoaEwcybxd/jQ9dvWjG0zlh7397A6\nndso0I9HY8xeY0y/gjyHUsqLnDwJHTuCzQYJCcy+9Thdvu/Cd52/08KfS7ku/iIyVUSOiMjmq7YH\niMgOEdklIm84LqJSSgEbN0L9+nDnnRAbS8iBuQyKHMSSnktodUcrq9O5nbxc+U8Drli5WER8gNCM\n7bWAQBGpmf94SikFTJ0KbdrA++9jxo7lzbjhhCaEsrLPSupWqGt1OreU6zZ/Y8wKEal61eaGwG5j\nzD4AEQkHOorIEeB9oI6IvGGMGZO/uEopb2GzxTNpnI0XtkZT++wf7P14LE26PE6/+c+w68Qufur7\nEzcVv8nqmG7LUTd8KwGJmZ4fABoZY04AL1zvxZlHq+lavkopmy2eEQPmErpvFXupwb3EcfP4EZS5\nEEql8reyrNcyihcpbnVMp3LU2r2XOKr453uUlhZ9pdQl4WN+YMa+aGbTheG8ByWPcK5pDJV3GdYO\nmkfhQt43LdmlGmnpAu4ZzT4LjTG1M543BoKNMQEZz98C0nPSzKMjfJVSV9i6laP1HmRU8rtMZBCU\n2wU9A2BTLx4hnbjYEVYndAmuMsJ3HVBDRKqKSFGgGxCR0xfrxG5KKQBWr4YWLZhWo5m98Ff6Gfo8\nAivfhLggivmlW53Qco6a2C3X/3cSkTCgKXCjiCQC7xhjponIAGAJ4AN8ZYzZnu90SinvERkJTz8N\n06dzb3oJbvq0LX/VXgcLpsJvHTJm5Ay4/nFUjujEbkop682cCUOGwLx5nK1Xm+d/fJ7Ve9ZSeVUb\nfE7d4nEzcjqCruGrlHJvEyfCRx/BsmVsvtnQ9csGPHzbw2wbsplirxezOp3HconZj7TNXykvZAwM\nH26fo2fFCqalJNBiRgvefvhtvnzsS4oV0cKfFV3MRSnlvtLS4KWXICGBvxfM5qV1I1h7cC2zn5xN\nrZtrWZ3OLbhKb5980St/pbxIcjI89RTs2sXOHybRKKID6SadhP4JWvhzQK/8lVLu5+xZ6NQJypYl\n7M1HeTnmNUa3HE3fun0RyfNFrFfSZRyVUu7h2DFo25bUuv9hQPtCxByI54euP3D/rfdbncwteURv\nn+DgYJ3eQSkPYbPFM3FiFMnJhfH1TeXll9vQ/r7boU0bjndoRcvqK7kntSYJ/RMo5VvK6rhux9Lp\nHRxJr/yV8hw2WzyDBi1hz55Rl7cFVHmOuRfms6Nfe/zL2QhuFsyL9V/UZp588ogrf6WUZ5g4MeqK\nwt+EVXyduIA3W5TkxwrxLO6ymHoV61mYUF2ixV8p5TDJyf+UlAAWM0N60sv/RtZV92HXc6sp41fG\nwnQqM+3qqZRyGF/fVMDwPz7la59AHgs0RJ76Hw/sfkILv4NoV0+llMuJmhkB/3uZG+UMgR2LsWvZ\nHKr7RjBhQoDOy+Ng2tVTKeUS9s75ilLPDWRmzXQmVW9M2QNNKFmksE7IVkC0+CuVS1l2RcxDcXLU\ncdzd5gMb2D2wO42X72JFcB/avTSekkVLWh3L42lvH6VyIauuiHv2vA2Qq8LtqOO4s01/bmLy96/T\n7+Pl/KfynZTaupsnK1ezOpbKKWOMpQ97BKWco02bt419OskrH/7+wyw5jjvaeHijeTz8cfNSYGnz\nd5kSJvmTD41JT7c6ltfJqJ15rr0uceWvI3yVs6QmFSKAxfRiBlXZxyEqcpgK+O7cDl9/DRUr/vMo\nWxayGYiUuUtjZklJPgWY3lobDm/g3bh32bZnDQtWVObufeUpFBMOdepYHc2r6AhfpXJj0yaYMYPj\noZPZffFeZtCLTfyH8vxJRQ7x32qz6PZQdTh0CA4ftv954YL9Q6BChSs/FCpW5I0JkUSsH8YhKnKG\n0pdP4+8/nMjI9yz8QR1v/aH1jIgbwfrD6xlbuhtd359PoVatYexYKFHC6nheS9v8lcrO4cP25QFn\nzIBTp+Dpp9k8cRL9PtrOnj0vXd6tevWh+E/4AK5uqz9//p8PgkuPw4dh82ZeS9nOc0Ue4taUZLZR\ni1l0Y+1tvzNwYDcn/5AFJ+FgAiPiRvDLn7/wVpPXmfN7A4qMCYXPPoMnnrA6nsonvfJXnuX8eViw\nwF7w16yBxx+HXr3gkUegkH1Mo80WT0hINElJPvlaG9Zmi+fTCUuo+ecftDqxhZZn91L03lrw5JPQ\ntStUquTon67ARfy4nJFfT2dPxVWcL3mUPjX6MrbFAPz69LcvwPLtt1ClitUxFdrVUylIT4f4eHvB\nnzcPGjWyF/xOnaB48St2Ncaw68QuVieu5vC5w/j6+FLUpyi+hX2v+LqoT1F8fXyv+7VfYT+K+hS1\nHzwlBZYuhe+/t38A3Xef/YOgSxcoX96CN+ba0tLT+O34b6w7tI51h9YRvXU5O07txJysAWsHwsY+\nPH9zT8YnLcfv/wbD0KHg47n3NNyNSxd/ESkBfAYkA7HGmJlZ7KPFX+XNzp3wzTf2R+nS8Mwz0L27\nvY0+w9nksyQcSmB14mpWH1jNmgNrKF6kOE2qNOH20rdzMe0iyanJ9j/T/vkz87Z/fT81+fLXSalJ\nlCxakiqlqlCldBWqlKpC5VKVud2vPPf/+ifVohK4IToeqVvX/kHwxBNw881Of6uMMew+sftyoV93\neB0bD2/klhK3UL9ifepXrM+scZtYt/BTSC6FHxf4mFdpxyImNmzKuJ+/dnpmdW2uXvyfBk4YY2wi\nEm6MeSqLfbT4qxyx2eIJHR9J6/2/0vboRqpJEn59e9uv8v/znyuu6lcfsD92n9hNnfJ1aFypMU2q\nNKFJ5SZUKuW45hhjDMcvHCfxdCKJZxI5cObA5a8TzySSeDqR4ycO0GmfH4HbfHho2zkS7y7P3tYN\n+fvRNvx+8CI/ztwJ50pTrAgMetk/3+MEjDHsP73/n0J/aB3rD6+ntG/py4W+fsX6PFDhAcoVK3f5\ndc2aBRMXF8y9bCGcp9hMbV7kc+o0HU9sbHA+3ynlaE4v/iIyFWgPHDXG1M60PQAYD/gAU4wxY0Tk\nTWCRMeZXEfnOGNMji+Np8VfXZbPFEzRgHh/s20JhUvmQ19l9VzS9g6tApfP/uqpvUtn+qFO+Dr6F\nfS3Nnm7SOfr3URJPJ3LoyG6KREZRefFP3LHud1beUoTwe4qzoPZZTt2QAumF8CvsS7Gifpebli41\nNV3dJJX5eVGfohQtVJTEM4msO7QOv8J+VxT6ehXqcXOJbP7HceoUrFjBnIEjqbY/hWrsZQjjmM4z\ngHhkDyZPYEXxfxg4B8y4VPxFxAfYCbQCDgIJQCBQDziZceUfZowJzOJ4WvzVdT3/4HO8vnoZc3xb\n8lYLH9JvXwXldlMmqRzP+nejSeUmNK7c2KFX9QWtY8vXKb68Lt2YRWuiOUgl1ktdzt51hG6fDOLC\nfXeTVLo4F9MuZtk8dfW25NRkKt5QkfoV61PhhgrZn/j0aVixAmJiIDYWfvsNGjdmZ4WqDF+WzPxD\nU0jBfh+jevWhOimbi3J6V09jzAoRqXrV5obAbmPMvoxQ4UBHYCIQKiLtgYjsjpl5elId7KX+Zc4c\nRq/7lpfua0FYu7mw9iVY+AX8WZf/PPQ+H38SbHXCPDmdVpwIAgknEB9SuYcdPGA28Ni5zyg9Zhyl\nf/kFbrwRHnjgysett+byRKdh5cp/iv3Onfab4s2awYQJ0KAB+PpyN/CMLZ4zIe9l6gmlhd9VOGpw\n1yV5avPPKP4LM135dwH8jTH9M573BBoZYwbm4Fh65a+ylp5O+vBh/D11Mq3anWftX91h+Ug4989V\nrTs3Sfj7DyMqamQW2zN+pvR02L0bNmywP9avt/9ZvDjUq3flB0KlSv+MRj5z5spiv337P8W+WTNo\n2BB8rW0KU/nnKoO88lW9dXoH9S+nT3OscwD7/9jMyLfup0eFvhx/ay97MhX+6tWHMnBggIUh8+fl\nl9uwZ8/bV0wOd8XPVKgQ3HWX/fFURl8JY2Dfvn8+ED7/3P6hAFC3Lpw8aS/2DRvaC/3YsVrsPYyj\n/gfgqOJ/EMg88qMKcMBBx1ZeZu+aSIp07sLyO324Yc5U5tXuiohQvXg8ISHDPaZJ4lL2XP1MIlCt\nmv1xaZStMSydMZf48bM4bSqwp1ELXhzSzq3fG1XwHNXsUxj7Dd+WwCFgLRBojNmeg2Nps48C4OSF\nk8z+qA+dPoxg/YDONB/5neU9ddxBVtNLV6/+NhMm5L/bqHJd+W32yfUaviISBqwC7hKRRBHpY4xJ\nBQYAS4BtwKycFP5LdA1f75aSlkLozyF83rkKXSdE4xOxkIDRs7Xw59DEiVFXFH6APXtGERISbVEi\nVZB0DV/lESJ3R/J2xGDGzjpJg5RbKL4w0i3nxLHSpcFZV2vaNFgHZ3kwV7nhmy96w9f7bDu2jVei\nXuHirp3EhBtuaNIO+fxz8POzOprb8fVNzXK7n1+ak5MoZ9D5/JVbOn7+OMGxwczaOovPinXliZFz\nkKFDYeDAbBdOUdeWdZu/Ds7ydC49t0+OAmjx9xpTNkzhrWVv8VStbozeWoESY0MgLAyaN7c6mttz\n1DTVyn1os49yecYYAr98hoW/26izsQsdxq8g1fe8fb79qlWtjucR2rd/RIu9l9BmH+UW0tLTeHRS\nZ2K2J3DTlMXMTerP79zBqGpVGB3SQQuWUnnk9K6eSuVUcmoyT815irW7fqX85Ch+SnqMeTxOIGFs\n2fuRdkVUykIuUfy1n7/nOZt8lnYz2wHQfE0Hlqc8xke8xmjeAuwXK0lJuiqUUrml/fyVyzr29zHa\nfteW+hXr82ndYRy+6z+MOz+UsbxyxX7uPCmbUlbTZh/lUvaf2s9D0x6iXY12fN5gBD6t2/D3E51Z\nUP3EFfvZJzBrbVFKpZRe+SuH2Xp0KwHfBfDag6/x8p097F04O3WCd9/VrohKOZhH9PMPCgrSrp5u\nbnXiah6f9Thj/cfSvUo7aNECWreG0aN18JZSDnSpq+eIESPcv/hbnUHlT+TuSHrN68X0TtNpe+t/\n7UW/SRMYN04Lv1IFxCOu/K3OoPJu5uaZDFkyhPnd5tOkbG0ICID774dPP9XCr1QB8ogRvso9hfwc\nwoerPmR5r+XcW7IatG8P99wDoaFa+JVycVr8Va4ZYwiKDSJ8Szgr+qygql956NgRqlSByZPtyw8q\npVyaFn+VK2npaby06CUSDiWwsu9KbilSBjp3hjJlYOpU8NGBW0q5A5co/jqxm3tITk2m57yenLhw\ngphnYihVqBh06wZFisC330Jhl/h1Usqj6cRuyqnOJp/l8VmPU9qvNN91/g4/CkPPnnD2LMydC766\n5KJSzqQjfFWBO37+OC1mtKB62ep83+V7/KQI9O0LJ07AnDla+JVyQ1r81TUlpSbxWPhjPHLbI0x6\ndBI+CLzwAiQmwvz5uuyiUm5Km31UttJNOs1CWrNr15/ctbkLfkVTmVx0G1VPHYPISChZ0uqISnkt\nl+7nLyLVgLeB0saYrgV5LuV43b7oRcLO30j6Yhd/pvryCa9wynctS76Zgr8WfqXcWoE2+xhj9hpj\n+hXkOVTB+HL9l9j2LiLp642Q6sv7DKUZsTRP3sK4r1ZZHU8plU85Kv4iMlVEjojI5qu2B4jIDhHZ\nJSJvFExE5WxRe6IYHjOc+zb1gPM3EcQI2mOjNdGcoqwuwqKUB8jplf80ICDzBhHxAUIzttcCAkWk\npog8LSLjRKSiY6MqZ9h8ZDM95/Zk9pOzKZtemqGM4km+pxVLOcGNAPj5pVmcUimVXzlq8zfGrBCR\nqldtbgjsNsbsAxCRcKCjMWY08E3GtnLA+0AdEXnDGDMmq+NnXpJMB3tZ59DZQzwa9igTAibw0G0P\nMb7Sd/gW+Z7/pmzhGLcAlxZhCbjOkZRSjuaowV2X5Li3T0bxX2iMqZ3xvAvgb4zpn/G8J9DIGDMw\nVwG0t49LOHfxHE2/bkrnezrz9iNv26dj/uwzlg4bzcdhv+giLEq5GCt7+zisYuv0DtZKS0+j+5zu\n1Lm1DkMfHmqfjjkkBGJjaXXbbbR65gmrIyqlMjjqfwD5Kf4HgSqZnlcBDuQvjrLCkCVDOJ9ynkmP\nTkK+/BI+/BDi4uC226yOppQqIPlp9ikM7ARaAoeAtUCgMWZ7rgJos4+lJqyZwBcbvuCnvj9RZuZc\nCA6GmBioXt3qaEqpa3DK3D4iEgasAu4SkUQR6WOMSQUGAEuAbcCs3Bb+S4KDgx16I0PlzIIdC/hw\n1Ycs6r6IMj8shHfegaVLtfAr5cJiY2Ov6CSTVzq9g5dad2gdbb9ry6Lui2iwYg/83//BsmVQs6bV\n0ZRSOeDS0zvklN7wda79p/bTMbwjUzpMocHq/TBkCERHa+FXyg3ofP4qT04nnea/U/9L/wf6M+hI\nVXj+efskbXXqWB1NKZUL+b3y1+LvRVLSUmg3sx01b6rJxHR/+5z8ixZBvXpWR1NK5ZI2+6gcMcbw\nwo8vUKxwMcYRAH16w8KFWviVcjPa7KNy5f0V7zNn+xxW3jaCYk/3gQUL4MEHrY6llMojj7jyVwUr\nbHMYk9dPZsPd4yjWq6996UUt/Ep5NZco/trsUzBstnhGTJvCL3fPoc+KlpTc1BfmzYVHdG4epdyV\nNvuoa7LZ4vnf2+H80XoujcLeIuLgKF6v0JKuX76oE7Mp5QGcMsJXuZ9PPo3gj4diqWfrTcTBUfTm\na6YfDiMkJNrqaEopF6DF3wOlpaex6c4FNNxYC9v2afTnSxbTDkBX4VJKAS5S/HVuH8d6Lfo1Hjp0\njoWr4unLVCLoePl7ugqXUu5N5/ZRWZqyYQoJXwQxIfxvni3xKDMPfXv5e9WrD2XChABt81fKA2hX\nT3VZ3L441o59hc+iilA4Korux5I4HjI80ypcWviVUnZ65e8hfj/5O6HP1eH9FUXxi46B2rWtjqSU\nKkA6t4/iTPIZPu9xFy+uSqFU3BqoUcPqSEqpAqbNPl4uLS2ViO716b0miVJrftWlF5VSOeISxV9H\n+OaRMax4sjFNEg5TLmE7VKxsdSKlVAHTEb7eLi2NnV2ac2FjArev2krZindYnUgp5UQ6wtcbpaRw\ntFMbjm1eQ/HYlVr4lVK55hLNPioXLlzg/OOP8ssfP+ETMY+7qup8/Eqp3CvwK38R6SgiX4hIuIi0\nLujzebSzZ0lt60/sX+vYPWUMLWu1tzqRUspNOa3NX0TKAB8bY/pdtV3b/HPixAlMu7ZEFTvMgsFt\n+fSxSYjkublPKeXmnNbmLyJTReSIiGy+anuAiOwQkV0i8sY1DjEMCM1rUK925Ag0b058pTQ+7H0n\nEx4N1cKvlMqX3DT7TAMCMm8QER/sBT0AqAUEikhNEXlaRMaJSEWxGwMsNsb84rDk3uKPP+Dhh9nw\n4B30a3qKH56cTRGfIlanUkq5uRzf8DXGrBCRqldtbgjsNsbsAxCRcKCjMWY08E3GtpeBlkApEbnT\nGDPZAbkAHggxAAAOxklEQVQ9ms0Wz8SJUdx08jSfbP6GrV1bEHjHSuK7x1OuWDmr4ymlPEB+e/tU\nAhIzPT8ANMq8gzFmIjDxWgfJPD2ptw/2stniGTRoCcX2PMVUAhhW7E2+rjCSoEojuOeme6yOp5Sy\niKMGd12Sqxu+GVf+C40xtTOePwEEGGP6ZzzvCTQyxgzMxTH1hm8m/v7DOBn1GBE8xmCfD5j13FjY\n0B//MseIjHzP6nhKKRdh9SCvg0CVTM+rYL/6zxVdzOUf9xw5wI88Sj8mM+vJuZD4IPw8UFfgUkoB\njlvMJb/Ffx1QQ0SqikhRoBsQke9U3mrxYkZs+4GnCMMWsByKnoNFoYDoClxKKYfKcbOPiIQBTYEb\ngaPAO8aYaSLSFhgP+ABfGWM+yFUAbfax++EHGDCAla8F0WHtN5wqexFmLIWksroCl1LqX5w2pbMx\nJjCb7YuBxXkNADqrJ19/DUOHwpIlLDkxhzLpx6gb35H0RhN0BS6l1BV0Vk9PERoKH34I0dG89+f3\nhG8NJ/aZWG4ucbPVyZRSLswjFnPx2iv/Dz6Ar76C+HhGHwjnu83fEdtbC79SKnt65e/OjLE38yxc\nCNHRjN0XxufrPieudxwVb6hodTqllBvwiCt/r5KeDi+/DGvWQGwsob+HE7o2VAu/UsqpXKL4e02z\nT2oqPPss/P47LFvG5F3hfLTqI+J6x1GldJXrv14p5fW02cfdJCdD9+5w7hzMm8fUneEExQYR+0ws\n1ctVtzqdUsrNWD3CV+XE+fPQqZP964gIvt09l+Exw1n69FIt/EopS7hE8ffo6R3OnIGAALj5Zpg1\ni1m75vN69OtEPx3N3TfdbXU6pZSbcdT0DtrsU5COH7cX/gYNIDSUuTvn89Kil4jqGUXtW2tbnU4p\n5ca02cdVHT4MTZtCy5bw6acs3GXjRduLLO6xWAu/UspyWvwdyGaLx99/GN0aD+bgnfeyo14TGD2a\nxbsj6bewH7buNuqUr2N1TKWU0q6ejnJpEZZCe54hmtaMIZhFPx2hT9hYJhwYTURgBPUr1rc6plLK\nzWlXTxfj7z+MQ1HdiCSA4bzHNPpC1ViK9mjPsueW8NBtD1kdUSnlQXSEr4uo+tdRptOaQUzge7rB\nbSuga1dqbumqhV8p5XK0zd8RVq5kzObv6M+X9sJfeTV0ewLmzKR8ko7cVUq5Hr3yz6/oaOjRgx1B\n77F92hooVAg69oV506kuMQwcGGB1QqWU+heXKP5ue8M3IgL69YO5c6nTuD61yvXgYGIod2/tTPnq\nq3QRFqWUw+kNX6uFh8PgwfDjj2y5zY/AOYHUurkWk9pPomyxslanU0p5OB3kZYWpU+GVVzBRUYSm\nr6H59Oa80uQVwp8I18KvlHILLtHs41ZCQuDjjzlum03vrW/z57k/WdV3FTVurGF1MqWUyjG98s+N\n0aNhwgTipr/L/cu7cN/N9/FT35+08Cul3E6BtvmLyD3AIOBGYIkx5qss9nH9Nn9jYPhw0ufO4d1h\nD/PVkcXM6DSD5tWaW51MKeWl8tvm75QbviJSCAg3xjyZxfdcu/gbA0OGkLRsCW17+VDu9rv54tEv\nuLH4jVYnU0p5Mafc8BWRqSJyREQ2X7U9QER2iMguEXkjm9d2AGxAeF5DWiYtDdO/P38uj+Derkfp\n0XIws7vO1sKvlHJ7ObryF5GHgXPADGNM7YxtPsBOoBVwEEgAAoH6wAPAR8aYQ5mOscAY0zGLY7vm\nlX9KCkk9A/ltSxwvPleRqT2+18VXlFIuwylz+xhjVohI1as2NwR2G2P2ZQQJBzoaY0YD32Rsawp0\nBvyAmLyGdAabLZ6JE6NITi5MySJJfHLKxoG/d7F0zAssbzsG38K+VkdUSimHyU9Xz0pAYqbnB4BG\nmXcwxsQBcdc7UOYlyawY6XtpOuY9e0ZRrNAp5pWty5YbD7Dz7ZF80CHL1iyllHIqR43svSTHN3wz\nrvwXZmr2eQIIMMb0z3jeE2hkjBmYqwAu0Ozj7z+MqKiR3FB2Az/6NmWfXxn6Hkmg1SOfEhn5nqXZ\nlFIqK1aO8D0IZJ6ysgr2q/9cs3oB97/TUrnjv/2JLdKQrWl16b1vH2kXypOU5GNZJqWUyoqjFnDP\nT/FfB9QQkaoiUhToBkTkO5ETXUy7SMjPIdxRcjxrEr7l61NB/O9YHAZ70ffzS7M4oVJKFYyc9vYJ\nA5piH6x1FHjHGDNNRNoC4wEf4CtjzAe5DmBBs48xhvk75jPC9hofRiTRYHcqz/i0Z+GBf8agVa8+\nlAkTdFZOpZRrclZvn8Bsti8GFuf15Jc4c0rntQfX8mrUq9y66zA/zUqmxCOtYVEIz8dt4GLIcJKS\nfPDzS9PpmJVSLkmndM6lfaf2MXTZUOL2xTL7aHMaT41Cxo+HHj0K/NxKKeVoHjGlc0He8D2VdIo3\not+g3hf1qFO4Cvtj69Jk2U5k9Wot/Eopt+OoG74ee+WfkpbCpHWTGLliJB3u6sDoou246bnB8OST\n8P77ULSow8+plFLO4pQ2f3dijGHBzgW8Hv061cpWI7p7JPd/GQGTXrIvwtK2rdURlVLKci5R/B11\nwzfhYAKvRL3CyaSThLQNwd/vXgjsAUWKwIYNUKGCYwIrpZRF9IZvJn+e+5NXol4hdl8s7zZ7l951\neuPzow3697evs/v66+CjA7aUUp7D65t9/jr/Fy2mt6B9jfbsHLCTkumFYfAQiIiAefPgwQetjqiU\nUi7HrXv7nEk+Q5NPHyLpl7IkvF+CQS1e43St2nD4MPzyixZ+pZTH8frePhdSLtBw4oP8sb4QZ8IS\n6M10PuR1xt/cgAenvkH7R5sWQFqllHINXtnsk5KWQrfZ3fhrbwrJYWuZyrM0IIHmxLD12H34hw7X\n4q+UUtfgEs0+uZFu0ukb0Zc0k8aD61sSTzOKc55G/MxW7gPQ2TiVUuo6XKL457TN3xjDoMWD2H9q\nP3MqDuGLjVOZwxM8RTjnKXF5P52NUynlqbyyzf+dmHf4cedCfkrpRbGRo1n70qt0n3GKPXtGXd5H\nZ+NUSnmD/Lb5u03xH7d6HF+t+YyEXxtTLGEjzJ8Pd96JzRZPSEh0ptk4W2vhV0p5PK8o/tM2TuOz\niOGsXHQrvpVuh+nT4YYbnJRQKaVcj0fM6nkt87bPY/a0V1k1ORXfDo/D7Nla+JVSKp9coqtndnP7\nLP19KXHv9GJ+TBGKTPsSOnSwJqBSSrkIj5/bZ+3vK9nSszVPHbmF4rYlcM89FqRTSinX5JHNPtu3\nxpHWqgXt/WpTfMOvWviVUsrBXK74H4hZwA0Pt+SGNo9y69I1ULq01ZGUUsrjuFTxPzl5AsU6dGbL\nm325b9JcKORS8ZRSymMU+A1fESkBxALBxhhbVvu0bT2U8cW24rcmkqhPB9D/mQkFHUsppbyaMy6t\nXwdmXWuHV5etYt+mJQx/41H69RrvhEi5V1ALzDua5nQszelY7pDTHTI6Qo6Kv4hMFZEjIrL5qu0B\nIrJDRHaJyBtZvK41sA04dq3jJ9TaR7v/dOdIdE1E8nzzukC5yy+E5nQszelY7pDTHTI6Qk6v/KcB\nAZk3iIgPEJqxvRYQKCI1ReRpERknIhWBpkBjoDvQX7Kp7G/VakD6j1+SnJSzVqi8/OVc6zVZfc8R\nvwCemvN6+3tzTmf8nef1PPl9vafm9KZ/65nlqPgbY1YAJ6/a3BDYbYzZZ4xJAcKBjsaYb4wxQ4wx\nh4wxw4wxQ4CZwBfZzuMw91swPjmejdNd3mhPzanFP+/nzO1r3Lmo5vU8+X29O/wbut5rnFH8czzI\nS0SqAguNMbUznncB/I0x/TOe9wQaGWMG5iqAiLWjzJRSyk1ZtZKXQ4p2fsIrpZTKm/z09jkIVMn0\nvApwIH9xlFJKOUN+iv86oIaIVBWRokA3IMIxsZRSShWknHb1DANWAXeJSKKI9DHGpAIDgCXYu3PO\nMsZsL7ioSimlHMXyWT2VUko5n06eo5RSXsjlir+IlBCR6SLyhYh0tzpPdkSkmohMEZEfrM5yLSLS\nMeO9DM8Yce2SROQeEflcRL4XkWetzpOdjN/PBBFpb3WW7IhIMxFZkfF+NrU6T3bEbpSITBSRXlbn\nyY6IPJTxXn4pIj9ZnSc7IlJZROaKyFdZzbhwNZcr/kBn4HtjzHPAY1aHyY4xZq8xpp/VOa7HGLMg\n4718AftNeZdkjNlhjHkReArwtzrPNVx3rioXkA6cBXxx7R54nYBKwEVcOKcxZmXG7+aPwNcWx7mW\n2sAcY8yzQN3r7eyU4p/LuYEqAYkZX+dsyK81OS2Tx5zDsE/H4TS5zSkiHQAb9tHiLpcxp3NVWZ0T\nWGGMaQe8CYxw4Zx3AT8ZY14FXnThnJd0xz5bgdPkMucq4DkRWQZEXvfgxpgCfwAPY/8k2pxpmw+w\nG6gKFAF+AWoCPYH2GfuEOSNfXnJm+v4PzsyYh/dTgDFAS1fOedXrFrhiRmAkMA57D7f5ZHSYcLWc\nmb5f1Nm/n7l8P3sAXTP2meWqOTO+dxv2KWqcljEP7+dg4OGMfa779+6UBdyNMSsypofI7PLcQAAi\nEg50BCYCoRltqk4dN5CbnCJyBHgfqCMibxhjxrhiTqAV0BIoJSJ3GmMmu2JOEbkFe5OfHxDjihmN\nMcMynj8DHDMZ/8pcLaeI3IO96awMEOKsjJDr380JQIiIPIx9zQ+nyWXO7UBfYKoTIwK5zrkIeCfj\nXune6x3bKcU/G5mbd8De5tfIGHMe+xvtKrLLeQJ7O7qryC7nQJxcAK4ju5xxQJw1kf4ly4yXnhhj\npjs9Udayey9HA/OsiZSl7HJeAFzpvlm2f+/GmGArAmUju/fzV6BLTg9i5Q1fdxlgoDkdyx1yukNG\n0JyO5lU5rSz+7jI3kOZ0LHfI6Q4ZQXM6mlfltLL4u8vcQJrTsdwhpztkBM3paN6V00l3rMOAQ0Ay\n9raqPhnb2wI7sd+5fsvZd9I1p+Z0h4yaU3MWRE6d20cppbyQK47wVUopVcC0+CullBfS4q+UUl5I\ni79SSnkhLf5KKeWFtPgrpZQX0uKvlFJeSIu/Ukp5of8HOP4gmf5WnhIAAAAASUVORK5CYII=\n", "text": [ "" ] } ], "prompt_number": 8 }, { "cell_type": "code", "collapsed": false, "input": [ "pl.semilogx(densities,abs(tau[:,0]-10**np.polyval(fits_7[0],np.log10(densities)))/tau[:,0])\n", "pl.semilogx(densities,abs(tau[:,0]-10**np.polyval(fits_11[0],np.log10(densities)))/tau[:,0])" ], "language": "python", "metadata": {}, "outputs": [ { "metadata": {}, "output_type": "pyout", "prompt_number": 9, "text": [ "[]" ] }, { "metadata": {}, "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAXkAAAEFCAYAAAAG45eHAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xd4VNXWB+DfIvTeewClCdKkSw1FQBABRRFF7CJXBL2f\nilgusSIWOiIoKKJSRCkqLUCCQAiEIjX0kgqEnkBIm/X9sTNhkkw5Z+ZMS9b7PHlMZk5ZhsyaPWs3\nYmYIIYTInwp5OwAhhBDuI0leCCHyMUnyQgiRj0mSF0KIfEySvBBC5GOS5IUQIh9zmOSJqB8RHSWi\nE0Q03srzbxLRvqyvg0SUQUTl3ROuEEIIPcjeOHkiCgBwDEBvAHEAIgEMZ+YoG8c/BOB1Zu7thliF\nEELo5Kgl3x7ASWY+y8zpAJYAGGTn+CcBLDYqOCGEEK5xlORrAYix+Dk267E8iKgkgL4AfjcmNCGE\nEK5ylOT1rHkwEMA2Zr7mQjxCCCEMVNjB83EAAi1+DoRqzVvzBOyUaohIFskRQggnMDM5e66jlvxu\nAA2JqB4RFQUwDMDq3AcRUTkA3QCschCorq+JEycaery153M/ZvQ99Z5j6zl/iFPrY74Wp9H39Oc4\nvfEa8pc4vfVad5XdljwzZxDRGADrAQQAmM/MUUQ0Kuv5uVmHDgawnplTXI7IQlBQkKHHW3te7z2M\nON/eObae84c4XY3R2WtInM7d05njC3Kc/vJaz0PvO5OzX+pWvm/ixIneDkETidM4/hAjs8RpNH+J\nMyt3Op17ZcZrLoa/i7qJxGkcf4gRkDiN5i9xusruZChDb0TEnrqXEELkF0QEdmPHqxBCCD8mSV4I\nIfIxSfJCCJGPSZIXQoh8TJK8EELkY5LkhRAiH5MkL4QQ+ZgkeSGEyMckyQshRD4mSV4IIfIxSfJC\nCJGPSZIXQoh8TJK8EELkY5LkhRAiH5MkL4QQ+ZgkeSGEyMckyQvhR1JSgHff9XYUwp/IzlBC+JGD\nB4HWrYHbt4GAAG9HIzxBdoYSogCJjwcyMoCEBG9HIvyFJHkh/Eh8vPpvdLR34xD+Q5K8EH4kLk79\nV5K80MphkieifkR0lIhOENF4G8cEEdE+IjpERGGGRymEAKBa8qVKAefOeTsS4S8K23uSiAIAzALQ\nG0AcgEgiWs3MURbHlAcwG0BfZo4losruDFiIgiwuDmjfXlryQjtHLfn2AE4y81lmTgewBMCgXMc8\nCeB3Zo4FAGa+ZHyYQghAteTvv19a8kI7R0m+FoAYi59jsx6z1BBARSIKJaLdRPS0kQEKIe4wJ3lp\nyQut7JZrAGgZ2F4EQGsAvQCUBLCDiCKY+YSrwQkh7sjIABITpVwj9HGU5OMABFr8HAjVmrcUA+AS\nM6cASCGifwC0BJAnyQcHB2d/HxQUhKCgIP0RC1FAXbwIVKoEVKkCZGYC168D5cp5OyphtLCwMISF\nhRl2PbszXomoMIBjUK30eAC7AAzP1fF6D1TnbF8AxQDsBDCMmY/kupbMeBXCBbt3A6NGAXv2APfe\nCyxZAjRv7u2ohLu5dcYrM2cAGANgPYAjAJYycxQRjSKiUVnHHAWwDsABqAT/Xe4EL4RwXVwcULOm\n+r5uXel8Fdo4KteAmdcCWJvrsbm5fv4KwFfGhiaEsBQfD9TKGvZQp47U5YU2MuNVCD8RH3+nJV+n\njrTkhTaS5IXwE3Fxd1rydetKS15oI0leCD+RuyUvSV5oIUleCD8hHa/CGZLkhfATlh2vNWuqcfPp\n6d6NSfg+SfJC+IHbt4HkZDUZCgAKFwaqVwdic09NFCIXSfJC+IH4eKBGDYAspsRI56vQQpK8EH7A\nslRjJp2vQgtJ8kL4AcuRNWbS+Sq0kCQvhB+wHFljJi15oYUkeSH8gLVyjbTkhRaS5IXwA9bKNdKS\nF1pIkhfCD1guaWBmTvKygrewR5K8EH7AWku+TBmgWDHg8mXvxCT8gyR5IXwcs/WOV0BKNsIxSfJC\n+LgbN4CAANVyz006X4UjkuSF8HHWSjVm0pIXjkiSF8LH2SrVALJ5iHBMkrwQPs7aGHkzWb9GOCJJ\nXggfJ+Ua4QpJ8kL4OHvlGul4FY5IkhfCx9kr11SrBly7BqSkeDYm4T8kyQvh4+yVawoVAmrXBmJi\nPBuT8B8OkzwR9SOio0R0gojGW3k+iIiuE9G+rK/33ROqEAWTtSUNLEnnq7CnsL0niSgAwCwAvQHE\nAYgkotXMHJXr0C3M/LCbYhSiwDKZgPPn1VZ/tsgwSmGPo5Z8ewAnmfksM6cDWAJgkJXjyMpjQggX\nJSYC5cqpNWpskZa8sMdRkq8FwLLaF5v1mCUG0ImI9hPRGiJqamSAQhRk9jpdzWQYpbDHUZLXsojp\nXgCBzNwSwEwAK12OSggBwH6nq5kMoxT22K3JQ9XhAy1+DoRqzWdj5iSL79cS0TdEVJGZr+S+WHBw\ncPb3QUFBCAoKciJkIQoOe2PkzaQln7+EhYUhLCzMsOsR29lxgIgKAzgGoBeAeAC7AAy37HglomoA\nLjIzE1F7AMuYuZ6Va7G9ewkh8goOVksNf/ih7WNSUoAKFYBbt9SQSpG/EBGY2el+T7steWbOIKIx\nANYDCAAwn5mjiGhU1vNzAQwFMJqIMgDcAvCEs8EIIXKKjwfatLF/TIkSqnP2wgWgRg3PxCX8h6Ny\nDZh5LYC1uR6ba/H9bACzjQ9NCBEXBwwc6Pg4c8lGkrzITT7cCeHDtIyuAaTzVdgmSV4IH6ZldA0g\nna/CNknyQviotDTg6lWgShXHx8qsV2GLJHkhfNT582qVyYAAx8fKrFdhiyR5IXyUljHyZtKSF7ZI\nkhfCR2ntdAWkJS9skyQvhI/S2ukKAJUqAampQFKS42NFwSJJXggfpadcQyQjbIR1kuSF8FF6yjWA\nJHlhnSR5IXyUnnINIJ2vwjpJ8kL4KD3lGkA6X4V1kuSF8FFSrhFGkCQvhA9KTgbS09XqklrJ+jXC\nGknyQvggcz2edKwiLi15YY0keSF8kN5SDaCOT0gAMjLcE5PwT5LkhfBBejtdAaBoUaBqVfUGIYSZ\nJHkhfJAzLXlAhlGKvCTJC+GD9I6RN5NhlCI3SfJC+CBnyjWAdL6KvCTJC+GDpFwjjCJJXggfJOUa\nYRRJ8kL4GGbnk7y05EVukuSF8DFXrgAlSwIlSug/19ySZzY+LuGfHCZ5IupHREeJ6AQRjbdzXDsi\nyiCiR4wNUYiCxdlWPKCWQShUCLh2zdiYhP+ym+SJKADALAD9ADQFMJyImtg4bjKAdQB0TMQWQuTm\n7MgaM1nDRlhy1JJvD+AkM59l5nQASwAMsnLcawCWA0g0OD4hChxnR9aYyTBKYclRkq8FIMbi59is\nx7IRUS2oxD8n6yGpBgrhAldb8tL5KiwVdvC8loQ9DcA7zMxERLBTrgkODs7+PigoCEFBQRouL0TB\nEh8PNGvm/PkyjNK/hYWFISwszLDrEdvphieijgCCmblf1s8TAJiYebLFMadxJ7FXBnALwEvMvDrX\ntdjevYQQyqBBwHPPAYMHO3f+4sXAihXAsmXGxiW8g4jAzE73dTpqye8G0JCI6gGIBzAMwHDLA5j5\nbotgfgDwZ+4EL4TQzoiOV2nJCzO7SZ6ZM4hoDID1AAIAzGfmKCIalfX8XA/EKESBIh2vwkh2yzWG\n3kjKNUI4lJGhJkLdugUUdvQ524bMTHWNGzeAYsWMjU94nqvlGpnxKoQPuXABqFTJ+QQPAAEBqtwT\nG2tcXMJ/SZIXwoe4Wqoxk2GUwkySvBA+xJUlDSxJ56swkyQvhA9xdWSNmXS+CjNJ8kL4EKPKNbJ+\njTCTJC+ED5GWvDCaJHkhfIhRNXnpeBVmkuSF8CFGjq6JiZHNQ4QkeSF8ilHlmlKl1NfFi65fS/g3\nSfJC+IiUFDXTtVIlY64nwygFIEleCJ+RkADUqAGQQXurSeerACTJC+EzjCrVmEnnqwAkyQvhM4zq\ndDWTco0AJMkL4TOMGj5pJi15AUiSF8JnGF2ukZa8ACTJC+EzjC7XSMerACTJC+EzjG7JV6kCJCcD\nN28ad03hfyTJC+EjjK7JFyoEBAaqma+i4JIkL4QPYDY+yQPS+SokyQvhE65fV9v2lSlj7HWl81VI\nkhfCBxjd6WomLXkhSV4IH+COUg0gLXnhZ0l+3z5ZOlXkT0aPrDGTYZTCYZInon5EdJSIThDReCvP\nDyKi/US0j4j2EFFPdwSalAS0awccP+6OqwvhXVKuEe5S2N6TRBQAYBaA3gDiAEQS0WpmjrI4bCMz\nr8o6vjmAFQAaGB3o1q1AZiawezfQuLHRVxcir+Rk9fcWEQE0aAAMHeq+e8XHAw0bGn/dwEB17cxM\n1bErCh5HLfn2AE4y81lmTgewBMAgywOY2XKqRWkAl4wNUQkNBSpXBiIj3XF1UdCZTMDRo8CPPwKv\nvAK0agVUqwZMmAAcPgwEB7v3/u4q1xQrBlSsqJYxFgWT3ZY8gFoALKdSxALokPsgIhoMYBKAGgD6\nGBadhc2bgbFjgXXr3HF1UdBcuQLs2qVa6RERwM6dQPnyQMeO6uu551SiL1YMyMhQs0cvXFCJ3x3c\nVa4B7nS+1q7tnusL3+YoyWvq5mTmlQBWElFXAIsAWC2oBFs0h4KCghAUFKQpyKtXVS1+9Ghg8mT1\noivsKHIhsmRkAIcO3UnoERGq5dy2rUroo0erFnz16tbPL1wY6NoVCAsDhg1zT4zuaskDdzpfO3Vy\nz/WFscLCwhAWFmbY9YjtDFchoo4Agpm5X9bPEwCYmHmynXNOAWjPzJdzPc727mXPqlXA7NnAhg3A\nPfcAy5YBLVo4dSlRgDADjz0GrF+vatPmVnqHDsC99+prKEydChw7Bnz7rfFxmkxA8eKqD6BoUeOv\n/+ab6pPI+DzDJoQ/ICIws9P7hTn6M98NoCER1QMQD2AYgOG5AqgP4DQzMxG1BoDcCd5VmzcDPbPG\n7LRrp+rykuSFI2fPAtu3q1ZshQquXatnT/ckeABITFSlInckeECVa44edc+1he+z2/HKzBkAxgBY\nD+AIgKXMHEVEo4hoVNZhjwI4SET7AEwH8ITRQYaGAj16qO/NSV4IRyIjVavd1QQPAM2bA5cvq7KK\n0dxZqgFkGGVB5/ADKzOvBbA212NzLb7/AsAXxoemJCaqP9A2bdTPbdsCCxe6624F20MPqRbliy8C\n3bsbt6G0t+zapRoFRihUCAgKUg2OESOMuaaZu2a7msms14LN52e8hoUBXbrcqZ+2agVERQG3b3s1\nrHwnIQEID1dvomPGAI0aAZMm+ffQu8hI45I8oD5Nbt5s3PXM3DmyBpCWfEHn80k+NPROPR4ASpZU\nCejAAe/FlB9t2AD06gW8/jpw8CDw88/A6dNA06bAoEHA6tVqlIq/yMwE9u5Vb1pG6dFD/T0azd3l\nmgoVVOfu9evuu4fwXX6R5M31eDOpyxtvwwagb1/1PZGqZX/3ndpwYtAgNXS1Th01OejkSe/GqkVU\nlBoSWbGicdds0gRISQHOnDHumoD7W/JEsoZNQebTST4hQU1Aadky5+Nt26rp5sIYJhMQEgI88EDe\n50qXBp5/Xo1S2bgRSEtT46179AB++UUlPV9kdKkGUMnSHa15d9fkASnZFGQ+neRDQ1UHYO41N6Ql\nb6z9+1WLt25d+8c1bQp8/TUQGwu8+iqwaJGaRTlmjFoh1JdERgLt2xt/XXckeXeXawDpfC3IfD7J\n5y7VAECzZuojc3Ky52PKjzZsAProWIyiaFG1WNe6dSq5V6kCDB6sRkDNmeMbtV93tOQB1T8UGmrs\nktfuLtcA0pIvyPwyyRctqhL93r2ejyk/0pvkLdWpA0ycqDppJ01S1+rd29j49EpNVYuK3Xef8deu\nX1+VbU6cMOZ6qanAtWvqjdKdpCVfcPlsko+OBm7cUNPPrWnXTuryRrh5U40n17iMkE0BAeqN4rff\n1OzKq1cNCc8p+/erEVglSxp/baPr8ufPq0XPCrn5lSgdrwWXzyb50FCVeGz98Utd3hhbtqiO7NKl\njble4cJqZM727cZczxnuKtWYmUs2RvBEqQaQck1B5tNJ3lqpxkySvDFcKdXY0qULsG2bsdfUw8iZ\nrtaYW/JG1OU90ekKqDeSixeB9HT330v4Fp9M8sw5FyWzpnFjNbzyyhXPxZUfuSPJd+3q3STvrpE1\nZnXrAqVKAUeOuH4tTwyfBNQnrOrV1cgoUbD4ZJI/fVrNrmzUyPYxAQFA69bAnj2eiyu/iYlRawMZ\n3UHZoYMadeONpSdu3FC1Z1t9OUYxqmTjqXINIJ2vBZVPJnlzqcbRAllSsnGNeSSM0Z1+pUurMfXe\n6Bjfs0ctQ12kiHvvY9Q6Np4q1wBSly+ofDLJb95svx5vJkneNe4o1Zh16aI2X/c0d5dqzHr0UJ3W\nJpNr15GWvHA3n0vyzHkXJbNFljdwXmamWqbAnUneG3V5d4+sMatZU41td3WhPE/V5AEZRllQ+VyS\nP3ZMTXa66y7Hx959N3DrlhprLPTZuxeoUcN9rcguXdTSxa62dPVy98gaS0aUbKRcI9zN55K8uVSj\nZcMKItWal5KNfu4s1QBqgk+VKmrmqadcvKg6Xhs08Mz9XO18TUpSn6jKlTMuJnukXFMw+VyS11qq\nMZO6vHPcneQBz5dsIiPVm767Z4+aBQWpfgdn19k3l2o8tQOXuVxj5Lo7wvf5VJI3mdROUFo6Xc2k\nLq/fjRuqXNOtm3vv4+nOV0/V482qVAECA51fQ8mT9XgAKFMGKFZM7VUrCg6fSvKHD6uProGB2s8x\nt+SldaJdWJgay+6OtV0sebolv2uXZ0bWWHKlZOPMyJpMUyZCzzhfI5LO14LHp5K8o1mu1tSqpcZE\nS4eSdpa7QLlTw4ZqlUVPJBVmz7fkAdc6X53pdA2PCUevn3oh9oZzU1el87Xg8akk72i9GlukLq+P\nJ+rxgKo1e6o1f+6cerP31Jhzs+7dgR071I5ZejlTrgk5HYIiAUXwy4Ff9N8Q0vlaEPlMks/MBP75\nx7klb6Uur92ZM6om37y5Z+7nqSTvyaGTlipUUJ9YnGlkOFOuCTkdgve7vo+F+xeCnahRSku+4NGU\n5ImoHxEdJaITRDTeyvNPEdF+IjpARNuJqIXeQP79Vy2gVKOG3jOlJa/Hhg1qL1dPjUDxVOerp2a6\nWuNsyUZvueZqylUcungIb3V+C2mZadiToH/hJmnJFzwOX+pEFABgFoB+AJoCGE5ETXIddhpAN2Zu\nAeBjAPP0BuJsqQZQLfk9ezw78cZfO3o9Vaoxa9UKOHvW/ZuIeKMeb+bsJiJ6W/KhZ0PRObAzihcu\njpEtR2Lhvwt131Na8gWPlvZcewAnmfksM6cDWAJgkOUBzLyDmc07e+4EUFtvIK4k+cqV1UbURm3J\npsVjjwGzZnnufkbIyFAtzgce8Nw9ixRRI3nCw913j8xM9Sbftq377mFP167qTUbPqpvMQEKCvk+u\nIadC8MDd6h9vRIsRWHJ4CdIy9XUGyOiagkdLkq8FIMbi59isx2x5AcAaPUGkp6u6rStb0HmyZHPx\nIhASAnz4oetrl3hSZKT6uF69umfv6+66/NGj6v+pYkX33cOesmXV0sY7dmg/5/JlNYS1RAnt54Sc\nDkGf+upj2N0V7kaTyk2w9sRaXbFWr672lE1J0XWa8GOFNRyjuTBBRD0APA+gs7Xng4ODs78PCgpC\nUFZW37MHqFdPtcidZU7yI0Y4fw2tli8HHnpIlT2efFLdV8+L1Vs8Xaox69JFvSG6izdLNWbmko3W\nT6N6SzWnr57GzfSbaFa1WfZjz7R8Bgv3L8SgewbZOTOnQoWA2rXVXgL29msQ3hMWFoawsDDjLsjM\ndr8AdASwzuLnCQDGWzmuBYCTABrYuA7b8tlnzOPG2Xxak9BQ5vvvd+0aWnXtyrx6NbPJxPz448xj\nx3rmvq7q1Ik5JMTz971xg7lkSeaUFPdc/z//YZ4yxT3X1mrDBuYuXbQfv2YNc58+2o//NvJbHvHH\niByPXUu5xuUmleNLNy9pvxAz9+jhnb8D4Zys3OkwV9v60lKu2Q2gIRHVI6KiAIYBWG15ABHVAfAH\ngBHMfFLvG43W9ePtad0a2L/f+XVEtIqNVTNz+/RR48C//RZYsQJYt86993XVtWuqtNSli+fvXaYM\n0KSJ+4a5emOma26dO6vdsG7e1Ha83jHyIadD0OfunB/DyhUvh/4N+2Pp4aU6IpXO14LGYZJn5gwA\nYwCsB3AEwFJmjiKiUUQ0Kuuw/wGoAGAOEe0jol1aA0hNBSIi1KQSV5Qtq/543b3q4bJlwODBag0Q\nQI2TXrgQeOEFtZWer9q8WSWi4sW9c3937fuamqr2WjV6C0O9SpZUMWzfru14PeWaTFMmNp/ZjN53\n987z3MiWI7Fwv75RNjKMsmDRNFqamdcyc2NmbsDMk7Iem8vMc7O+f5GZKzHzfVlfmttVu3apTbnL\nl3fuf8CSJzpflywBnngi52M9eqi+gBde8N2hld6qx5u5q/P1wAG1tLC71+HRQs86NnrGyO+O341a\nZWuhRpm8Q3F6390bMddjcPTSUc1xygibgsXrM16NKNWYuTvJnzqlPuZai/fjj9ULd57uGQLuxwys\nX++Z9Wps6dzZPZuI+EKpxkzPpCg95RprpRqzwoUK46nmT+Gn/T9pjFLKNQWN15O83vXj7XH38gZL\nlwJDhwKFrYxJKloU+OUX4P331ZA+X3LqlFpbpWlT78VQvTpQqZLx5TRfGFlj1rGjKh1dv+74WD3l\nmg2nNuCB+rYnN4xsORKLDiyCibW9g3q6XJOaqibEhYerkWkzZgDjxwNPPw306qX6ax55xHPxFDRa\nhlC6TUqKSspGdQa2agVERalJKe6oPS9ZAsyebfv5e+5RLfqnnlJjposWNT4GZ6xff6ej2JvMJRsj\n182JjATGjTPueq4oXlx9qti6VQ2xtUdruSYpNQn7zu9Dt7q2F/9vXq05qpSsgrCzYeh5l+MWU2Cg\nGkBgMrm+vMXFi+pTQXy8+oqLu/O9+evGDTXpq2bNnF9Nm6r/Vq+ukv2pU0D9+q7FI/LyapIPD1cv\n+DJljLleiRKqvr9/v5plaaTDh9XU/M5WZwDcMWoUsGYNMHEiMGmSsTE4a8OGvP0I3tC1K7BpEzB6\ntDHXS0pSLcRmzRwe6jHm8fL2knxGBnDpktoi0ZEt57agXc12KFnEfqeDuQNWS5IvUULt23Dhgr4Z\nt0lJak7Lrl13vpKT1V7Llsm7U6ecP1eu7PjNZPhwYNEiwGIqjTCIV5O8kaUaM3Nd3ugkv3QpMGyY\n4z9WIuD779Wnir59XZvFa4T0dLVJyPffezcOwPhJUXv2AC1bqqUTfEXPnsCYMfaPuXBB7SplreyX\nW8ipO7Nc7RnebDiCw4KRnJaM0kVLOzzeXJe3leTT04FDh+4k85071Rtqy5bq08qjjwKTJ6sEb8Qn\nxJEj1VIhEyd6/xNnfuPVmrwr69XY4o66PLP1UTW2VK0KzJ8PPPOM+xfmciQiQo0+qVLFu3EAakne\nlBTj6sG+VI83a9cOOHkSuHLF9jF6RtZsOL0he70ae6qVroaudbvij6g/NF3XcoQNM3D6NLB4MfDG\nG+rTaoUKqma+c6f6f/rpJ/W3vH07MHWqannXr29cQm7dWn3C8OROYgWF15J8crIqq3TqZOx13THC\nZt8+Vb9s00b7OQ8+CAwapEoT3hxW6aldoLQwehMRby4vbEuRIupvessW28do7XSNvRGLxJuJuK+G\ntkkAI1uM1DzKpm5dYO5coH9/1QDo3l11itaoAXz6qVo87dAhYMEC4JVXVBJ25ycmItUo+kn7ICGh\nkdeS/LZt6g/H6PHNzZqpj5VJScZdc/Fi1YrX22qZPFm9UH7+2bhY9PL2+PjcjEzy3tooxBFH4+W1\ntuRDToWg9929UYi0vUwHNh6Ifef3IeZ6jMNjH3tMjQYaNUrNNYiJAX7/HXj7bVViNKqfTI+nnlIx\nyOJpxvJakndHPR5QrY0WLVTr2wgmk6rHO9NxWaKEGlb53/+qHZk87fJlNdro/vs9f29bjJr5mpio\nlmpo0MD1axnN0Xh5rWPktZZqzIoXLo7Hmj6Gnw84blXcf79qsQ8apH8LQnepVUuVW1evdnys0M6r\nSd7oerxZ27bGlWx27FBLJjg7gqNlS2DCBDUj1t3r6uS2aRPQrdudJRh8QatW6g3P1b4Kcz3eUztc\n6XHffWqI4sWL1p/XUq4xsQmbTm+yOz7empEtR+KnAz85tTWgLxg5Ui0TIozjlZfI9etq0kjHju65\nvpF1eT0drra8/roqS332mTExaeVrpRpAfdJq317f2uvW+GqpBlCjZrp1U6OarNFSrtl/fj8qlKiA\nOuXq6Lr3/bXvR4YpA5Hx/rkf5pAh6m/j/HlvR5J/eCXJ//OPSvDuamEaleQzMoDfflNDJ11RqJBq\nncyerUa7eAKzbyZ5wJh9X31xZI0leyUbLeWaDac22FzKwB4i0tUB62tKlVILAP76q7cjyT+8kuTd\nWaoB1GYIiYn2h7FpsWWL2mChYUPXY6pZUy1LPGKEsZ3Cthw7pjqKGzd2/730crXzldk3R9ZYsrfv\nq5ZyTcjpEN2lGrOnWz6NpYeX6t4a0FeMHCmjbIzklSRv5KJk1gQEqJE7ro6XN6JUY2nIEDVywRPT\n8H1lKQNrOnZUHeN69kS1dO6cKono2VnJ01q0UB3fcXE5H09JAW7dsr9VYUp6CnbG7URQvSCn7l2v\nfD3cW+Ve/H38b6fO97bu3VWfzf793o4kf/B4kr98WU28cPdHbVdLNmlpwB9/AI8/blxMADBtmipV\nLF9u7HVz89VSDaCG591zj5qx6gxPl2rWnVyHdze9q+ucQoVUssrdmo+PV2PR7b35bo3eilbVW6Fs\nsbJORKuYO2D9UaFCaiKWtOaN4fEkv2WLmlHn7qnorib5kBC1Ol4dff1eDpUurYZVvvqqGoHhDqmp\n6o2kVy/3XN8IrpRsPJnk0zPTMXbtWEzfOR3R1/VN1bVWstFSqtlwSt/QSWuGNh2K0DOhuHTrkkvX\n8Zann1beJuKlAAAfk0lEQVR1eU+PSMuPPJ7k3V2qMXN1eQOjSzWW2rcHxo5VM/yMXl8dUAu/NWli\nvyTgba50vnpyDfkF+xagTrk6GN12NKZFTNN1bs+eeTtftYysCTkd4nKSL1usLAY0GoAlh5a4dB1v\nadxYzcoNCfF2JP7P40neXZOgcrvrLlXzTUjQf25KCvDnn2rteHd55x1VEpoyxfhr+9JSBrZ06eLc\nJiKZmcDevepN3N1upd/CR/98hM97f45xHcbhx39/xLXb1zSf36SJ+ls6e/bOY45a8ueTzyP6ejTa\n1XL9o4o/j7IBVCNIxsy7zqNJ/sIF9Ufuif04iZyfFLVmjTq3enXj4zILCFBLq06ZosbPG9mi9+V6\nvJl5E5EjR/Sdd+yYWgDOE59SpkdMR+fAzmhbsy0CywViYOOB+Hb3t5rPJ8pbsnHUkt94eiN61OuB\nwoVcXyC29929EZcUh6jEKJev5Q3DhgFr16qZzcJ5Hk3yoaFqWntAgGfu166dcyWbJUvUKnvuVq+e\nehNaswZ4+GHXh3wCapblqVPGL7XsDs7U5T01dPJKyhVMiZiCT3p+kv3Ym/e/iRk7ZyA1I1XzdXKP\nl3c0Rt6IUo1ZQKEA3VsD6hVyKgQp6e5ZbKZiRaB3b/cPUsjvPJ7kPVGPN3OmJZ+UpFrCQ4a4J6bc\natVSv5fGjdUql64O+9y4UQ3T9KU11m1xJsl7aqbrpK2T8GiTR9GoUqPsx5pXa46W1VtqWhvGzNyS\nN68yYK9cw8wIOeX8+HhrRrYciZ8P/oxMU6Zh1wTUsgv/t/7/8MiyRzDg1wFITks29PpmMmbedR5P\n8p6ox5uZR9joWcZj9Wo1Jd2TnZZFigBffw189ZVa+vXbb51fntgfSjVmznS+emJkTfT1aMzfNx//\n6/6/PM+91ektfLXjK837qTZooMo2J0+qn+2Vaw4nHkbxwsVRv4Jxe+A1q9oMVUtVRehZO8ti6pSa\nkYoRf4zArvhdODPuDOpXqI8+i/ro6q/Q6sEH1Z7Jp08bfukCw6NJ/soVY/f3dKRmTbV0gmXHlyPu\nHFXjyKOPqk0ZvvlGDSG7eVPf+b68lIE1jRrp20QkNVVtw+juPp3gsGC80vYV1CyTNxv3qNcDJYuU\n1DzRyFyX37xZ/fvYK9eEnFKlGjJ4BpuRHbA3Um+g/6/9kZqZipCnQ1C5ZGXMGzgP7Wu1R4+FPZB4\nM9GQ+5gVLapej4sWGXrZAkVTkieifkR0lIhOENF4K8/fQ0Q7iOg2Ef2fret07+75VQP11OWvXFHr\n6jz8sHtjsqdhQ7W+TeHCqq5+9Kj2cw8fVssb++Lyu9aYNxHZvl3b8QcOqP+3UqXcF9ORxCP46/hf\neLvz21afJyK81ektfBn+peZrmks216+rT22lbezO58pSBvYMbz4cq4+tRlKqa+tpxCfFo9sP3dCk\nchMsG7oMxQsXB6B+J1P7TsWAhgPQ/cfuiE+KNyLsbOaSjZ8urOl1DlMuEQUAmAWgH4CmAIYTUZNc\nh10G8BqAr+xdy5OlGjM9dfk//lCtYG9smGCpZEnghx/U6pVduwLLlmk7z7yUgT/RU5f3RKnm3U3v\nYnzn8ShfvLzNY4Y2HYqYGzHYGbtT0zXNm4jExtpuxadmpGJb9DZNG3HrVbVUVXSr203z1oDWHL10\nFJ0XdMawe4dh5oMzEVAo5+gJIsInPT/ByJYj0e2Hbjh77ayLUd/Rpo36RB4ebtglbcrMBObNy18b\nl2hpV7cHcJKZzzJzOoAlAAZZHsDMicy8G0C6vQt5stPVTM/MV2+WanIjAl58UZVfJkxQ692kOVhv\nyp9KNWbWknx6Zjomhk7EiqgVOR5398ia8Jhw7E3Yi1fbv2r3uMKFCuO/Hf+ruTVft6769LFxo+0k\nHx4TjiZVmqBiCfd0BrmyzEF4TDiCfgxCcPdgTOg6wW456Z0u72Bch3Ho9kM3HL983NlwczBvDeiJ\nMfOTJqndsfr2zUdDN5nZ7heAoQC+s/h5BICZNo6dCOD/bDzHJhN73KVLzGXLMmdm2j8uIYG5XDnm\nW7c8E5ceV68yP/wwc4cOzNHR1o+5dYu5dGl1rD9JS8sZd8z1GO40vxP3XNiTK39RmffG780+tmlT\n5r17bVzIRSaTibsu6MoL9i7QdHxyajJX/qIyH790XNPxzz/P3LYt89NPW39+wsYJ/N6m97SGq9vt\n9NtccXJFPnftnK7zVkat5CpfVOG1J9bqOm/+3vlc8+uafOD8AV3n2RITw1yhgntfn+HhzFWrqtfY\n2LHMzZszx8W5735aqTRtP0/b+9LSkjesEvbhh8EIDlZfYbZ2VDBYpUpA5crAcQeNiuXLgYEDVU3b\n15QvD6xcCTzyiPpksn593mO2bVMrH5a3XWXwSUWKqP+n8HC1EFjbeW0xsNFAhDwdgm/6f4MhS4fg\n4s2LSEpSHejO7tDlyJoTa3A55TJGthyp6fhSRUvhlTavYMoObVOWe/ZUfUM2O10NHB9vTbHCxfB4\n08d1Df+cu3suRv89GmueWoN+Dfrput/z9z2PKX2m4IFFDyAyzvXNHWrXVmWbP/90+VJWXb+u9pid\nOxcIDFQLCT75pFpny1HuMFpYWFh2ngwODnb9go7eBQB0BLDO4ucJAMbbONZuS95bHn+c+aef7B/T\nuTPzX395Jh5XhIUx16zJPHEic0bGncfffJM5ONhrYbnk/f+lc6f33+NaX9fisDNhOZ57d+O73O2H\nbrxxcxrff7977p+RmcHNv2nOK6NW6jrvfNJ5rvB5Bb6QfMHhsXFxzADz9Ol5n7t08xKX+awMp2ak\n6rq/XjtidnCjmY3Y5OAjtclk4g82f8ANZjTgk5dPunTP1UdXc5UvqvA/Z/9x6TrM6jU8YIDLl8nD\nZGJ+8knmV17J+9z8+czVqzPv2mX8fbWCB1ryuwE0JKJ6RFQUwDAAtrba9cHVyx3X5aOj1YbXD7iv\nIWWY7t1VizAsTI2pv5S1yKA/rFdjTUJSAv4s/wCOXN+JPS/vQfd63XM8/3HPj1G2WFn8L/wNt3W6\n/nLwF5QpVgYPN9Y3rKpa6Wp4rOljmL1rtsNja9ZUQ0atteQ3ndmEbnW7oWhAUV3316tDrQ5gZuyK\n22XzmAxTBl5c/SLWnlyL7c9vR/2Kro3ZH9h4IH599Fc8uuxRbDi1waVrPfKIGol14YJLlwGgGreX\nb13GoYuH8M68EIRdXYgaj3+OcWvHYfRfo7NX73z+edW6HzBAvcb8EbGGcUlE9CCAaQACAMxn5klE\nNAoAmHkuEVUHEAmgLAATgCQATZk52eIarOVe7rBli+q8tNU7/9VXak2U777zbFyuyMgA3n9fLcc6\nfTrwwgtqSYPCri954jGhZ0IxYsUIPNPsZUx/9H1cuRRgdUvI67evI/CjDni89pv4/j8vGhpDakYq\nGs9qjEVDFqFr3a66zz9++Ti6LOiCs6+fRckiJe0eGxEBNG2qNoa39NLql9CsajOM6+j+3WQ++ecT\nJCQlYPaAvG9MN9Nu4vHlj4OZseyxZShd1MZYTydsj96OIUuH4LuB32HQPYMcn2DDs88CLVsCb7xh\n/flMUyYSbyUiISkBCckJ2f+NT4rP8fP55PMoWaQkKhWrgehDNdGnUw00DayBGqVr4MSVE9gavRUb\nn96IaqWrAVDl0EcfBaZOVWUcTyIiMLPzDWhXPgbo+YIXyzU3bjCXKqU6+axp04Z540bPxmSUVauY\ny5dnfuQRb0eiXaYpkz8K+4irf1WdQ06FMDNz69bM27bZPqdWi6NccVIV3h693dBYpu6YygN+ca0G\nMGTJEJ61c5ZT55pMJq4ztQ4fvnjYpRi0Onv1LFeaXIlvp9/O8fjF5Ivc/rv2/NzK5zgtw8YLxUW7\n43ZztS+r8a8HfnX6Gps2MbdqdefnjMwM3npuK7+94W1uOrspF/6oMFf+ojK3mNOC+y7qy8+ufJYn\nbJzAMyJm8PLDy3l79HY+feU030q7xWlpzO3b5y2hmUwm/jDsQ248szHHXo/NfvzgQebatZmnTXM6\nfKfAxXJNgUjyzMxNmjDv25f38ePHmatVy1nf9jfR0cxnz3o7Cm0uJl/kPov6cLcfunHcjTtDF8aO\nZf78cxvnXFRvZH8e/Ztrfl0zxwvPFddvX+eqX1Z1eQRIeHQ43zXtLk7PTNd97rFLx7jW17Uc1smN\nFPRjEP9+5Pfsn09dOcUNZzTk9za95/Y4Dl44yDW/rsnf7fnOqfMzM5lr3X2Dv1rzG49cMTI7ob+3\n6T2OiInQ1a8xYQJz//5sc9Tf5G2Tuf70+nz26p0X19mzzI0bM7/zju3zjCZJXqORI5nnzcv7+Mcf\nM48Z4/l4CqKt57Zy7Sm1ecLGCXkS4m+/MT/0kPXz/v6buVcv9f1n/3zG7ea145T0FJfj+WDzB/z0\nHzbGNOrUeX5nXnpoqe7zZu2cxc+ufNaQGLRasHcBP7z4YWZWresaX9Xgb3Z947H7H790nOtMrcPT\ndmhvEp+9epZn7pzJfRb14aITS3Pd9/vwrJ2zciRgPTZvVgMYLjjoM58eMZ3rTq3LJy6fyH4sMVF9\nAnj+eeZ0/e/rukmS12jmTOaXX877+L332i8TCNdlmjJ58rbJXO3Lavz38b+tHhMfr8ZBW5vPEBys\nWl3M6qP04789zs+seMalVuf5pPNccXJFPnP1jNPXsLQyaiW3nddWd0yDFg/iXw78YkgMWt24fYPL\nTSrHP+//mat8UYVXRK3w6P2ZVdJuMKMBf/rPp1afzzRl8s7Ynfzepve4xZwWXPmLyjxyxUhefng5\n7z5wg6tXdz7BXrqkyi7r1mk7fu7uuVx7Sm2OSozKfiwpiblvX+aBA5lv3nQuDq0kyWu0Ywfzfffl\nfOzgQeY6dRxPlBLOu3zrMj/060Pc8fuODifi1K/PfOhQ3scHDGD+4487PyenJnPLOS11tQRze/Xv\nV/n1ta87fX5umaZMbjyzMYeeCdV8TlpGGpebVE7TEEyjPfX7U1zh8wq87Zz3WjjxN+L53tn38oSN\nE9hkMvHNtJu86ugqfmHVC1z9q+rcZFYTfnvD27zt3DbOyMxZT23fnnmtvvlZzKxKLIMHM//3v/rO\nW/jvQq7xVY0cpb3UVOannlLDr69c0R+LVpLkNUpJYS5RQv3X7L33mN96y3sx5XcRMRFcd2pd/u+6\n/2qqlT7zDPOcOTkfM5mYq1RRMx4tnbl6hqt9WY03ntLfY37y8kmuNLkSJ95M1H2uPfN2z+P+v/TX\nfPy2c9u41betHB/oBglJCU6XOoyUeDORW89tza3ntuYyn5XhHj/24CnhU3KUR6yZNYt5+HD99/v2\nW9XYu33b8bG5LTm4hKt9WY13x+3Ofiwzk/mNN1RFINaYrqI8JMnrcN99zFu23+aTl09yZqaJ69dn\n3rPH21HlPyaTiaftmMZVvqjCfxz5w/EJWb77TrWMLJ09y1yjhvVOrs2nN3O1L6vx6SundcU3fPlw\n/ijsI13naJGSnsLVv6rOhy5Y+ThixcTQifzWBmllXEu5xquOruIrt7Q3hy9dUsuQXLum/T6HDzNX\nrsx89KgTQWZZEbWCq3xRhcOjw7MfM5mYJ09mrluXOSrK9rnOcjXJ+9GoatfV7RSJYZufhSn8ElLT\nMpDRoxPW3+yMm+fUPp4livjgmgYuMrEJt9Jv4WbaTSSnJef4upme87GbaTfBYBSiQiAQClEh9T3d\n+V7Lc+tPrUf09WhEvBiBuyvcrTnWrl2BTz7J+Zh55Ulra2L1uKsH3uv6HgYvHYzw58NRqqjjNYj3\nJuxF6NlQzBs4T3NcWhUvXByvtX8NX+34Cj8M+sHh8RtObcCHQR8aHoe/KVe8nO6JaJUqqaUili9X\nc0QcuX1bbek5aZLahc1Zg+8ZjGIBxTBoySAsf3w5utXtBiK1qFnVqmpXtlWrfGv7TU2ToQy5kRcn\nQ6VmpOLDLR9iVvh8tEqYhi2zn8DotxNwqcR21O60HeEx4TiceBjNqzZH58DO6FynMzoFdkL10m7c\nydtFJjbh2KVj2Bm3ExGxEThz7UyehJ2cloxb6bdQokgJlC5aOvurVJFSOX42P1aqaCkUokIwsQkm\nNoGZs783sQmMXD/beL522doY33k8ihW2MrPJDmb1Qtm7V60fAqgXT9myauKX9XMYL6x+AUlpSVg2\ndJnDDTf6/twXDzd62OFKk866mnIV9WfUx8HRB1GrrI19/qAmeNWeWhsX37yYLxsXnrBypZqctGWL\n42Nff13tyrVsmfUGg16bz2zGE8ufwC+P/JJjD4C//gKee05tctJP33I/Nrk6GSrfJ/nIuEg8u+pZ\nNKrUCK/WnYOxz1fHoUNq+df169UMREDN9ouMj8T26O0Ijw3HjpgdqFCigkr6gSrp31v1XhQiD+96\nkuVKyhXsituFHTE7EBEXgV1xu1CheAV0rN0RHWt3RKNKjVCmaBmUKpozgZcsUtJrMTtjyBBg2LA7\nSz736KFmK9tbQjk1IxXdf+yOgY0G4r1u79k8bvOZzXjpz5cQ9WqUW5cQeH3d6ygaUBRfPPCFzWNW\nHl2JbyK/wYan/XSuvA9IS1P75e7aBdx1l+3j1qwBRo8G/v0XqFDBuPtvi96GR5Y+ggWDFuChRg9l\nPx4erv6O33hD/R3Xq+fafSTJ22Buvc/fNx/T+k7DE82eQEYGoUIF9W7+zjtqpyFbTGxCVGIUtseo\nlv72mO24dOsSOtbumJ3061eoj0olK6FUkVKGbtmWYcrAoYuHEBEbkf0VlxSHdjXbZSf1DrU6ZE+5\nzk++/lrt5zl7NmAyqRflmTOO99yNT4pH++/a49uHvs3xgjNjZnT4vgPe6PgGhjcf7qbolXPXzqH1\nvNY4M+4MyhYra/WY//z9H9xV/i681fktt8aS3732mvr098EH1p8/f15tF7l0qdq72Wi74nZh4OKB\nmDNgDh5p8kj241FRwBdfAH//reJ76CG1ym3HjkBAgJ0LWiFJ3gpz671xpcaYM2BOjmTYqROQnKze\nYd99V991LyRfQHhMuPqKDUf09WhcvnUZmZyJSiUqoXLJyqhUshIqlcj6Kpn1WNb3lo+VL14+u4Wd\nkJSQXXaJiI3AnoQ9CCwbmJ3QO9buiHur3JtnN578aOdO4OWXgf371Qtl4MA7m2A7EhEbgYcXP4x/\nnvsH91S+J8dzy48sx2dbP8Pul3d75JPNk78/idY1WuPNTm9afb7hzIZY/thytKze0u2x5GeRkWot\nmePH85ZhTCa1EXiHDsBHH7kvhn0J+9D/1/6Y0mdKngaEyaQ+afz1l/qKjVUxDRyoFhQsV87x9SXJ\nW7DWes/dwh47Fpg5UyWO+q4tsJctJT0Fl1Mu49KtS7h86zIup1zO8d9LKXkfT0pNQvni5VE0oChu\nZ9zOkdDb12pvd/u5/Cw9XbXaY2JUB9a6dcDixdrPX7BvASZvn4ydL+7M/h1mmDJw7zf3YuaDM9Gn\nvme2ztqXsA8PL3kYp8aeylMaOnvtLDp83wEJ/5fgV6U0X8SsSq7z56sGnKUpU1TH7D//uH/hvsMX\nD6PPz33wac9P8WyrZ20eFx2tWvd//QVs3aq2JzW38hs2tH6OLFCWZVfsLm46uykPWTKEzyedt3nc\nL7+oiRTelp6ZzheSL/C5a+c8um6JP+jRg3nNGuZXX2WeMkX/+WP+HsP9f+mfPYFm7u653HNhT4//\nnnv/1JsX/rswz+Pzds/jJ39/0qOx5GeffcY8alTOx/buVcMlT+sbXeuSY5eOceCUQJ4TOcfxwcyc\nnMy8ejXzSy+pYcKNGqlJWqGhORdTREEfJ387/Ta/E/IOV/2yKi8+uNjhCzk9XY2xFb7rgw+Y331X\nvRlv3ar//LSMNA76MYgnbJzAN9Nucs2va/KuWM/v+rD+5Hpu9k2zPH+TQ5cN5R/2/eDxePKr6Gjm\nihXvTHRMTlaLiP3q/GKXTjt15RTXm1aPp+6Yqus8k0nN2QkOVttEVqjAPGwY86JFBTzJa229C/+y\nfr3az7ZkSfWCdUbizUSuN60e91nUh4cuG2psgBqZTCZuOaclrzm+JvuxjMwMrji5omEraQqlVy/m\nZcvU9y++qBYk9Jboa9HccEZDfnvD23wrzblNaePjmb//nnnQoAKa5PW23oV/uX6duVAh5hYtXLvO\nvoR9XO3Lanw00YUpji5atH8R9/ixR/bP5oaJMNbChWoV099+U2sg3bjh3XgSkhL4sWWPcd2pdfmP\nI3+4lKNcTfJ+1/Fqb+SMyD/atFFD377/3rXrmNjk1c7N9Mx01J9RHyuGrUCbmm3w6T+fIvFWIqb1\nm+a1mPKj5GS12XfRoqpTs317b0ekbDq9CWPXjUWtMrUw48EZeUZ9aeFqx6vfdO1fvHkRb4e8jYcW\nP4QPun2A3x//XRJ8PjZ8uBp14Cpvj14pElAEr3d8HV+GfwkACDkd4rERPgVJ6dJqP9Z33/WdBA8A\nve7uhX9H/Yv+Dfuj6w9d8eaGN3Ej9YZHY/D5lvy/5//F9J3TsfLoSjzW9DF83ONjSe7CrySlJuGu\n6Xch7NkwdPy+Iy68eUHTOjsif7mQfAETNk3AupPr8HnvzzGixQhNjZB8OU4+w5SBVUdXYfrO6Thz\n7QxebfcqXmr9EiqVrOTmKIVwjwkbJ+DP43+icsnKCHs2zNvhCC+KiI3Aa2tfQ5FCRTCr/yy0rtHa\n7vH5KslfSbmC7/d+j9mRsxFYNhDjOozD4HsGo0hAEY/EKIS7JCQloN70evhft//ZXV9HFAwmNmHB\nvgV4f/P7GHzPYHza81Objdh8keSPJB7BjJ0zsPTwUgxsNBDjOoxDm5ptPBKXEJ4yf+98PFD/AdQp\nV8fboQgfcTXlKiaGTcSSQ0sQHBSMUW1G5Vm+xO1Jnoj6AZgGIADA98w82coxMwA8COAWgGeZeZ+V\nY3IkeRObsObEGkzfOR2HLh7CK21ewai2o3x6eV8hhHCHgxcO4rW1r+Ha7WuY1X8WutTpkv2cW5c1\ngErsJwHUA1AEwL8AmuQ6pj+ANVnfdwAQYeNazMx8/fZ1nrZjGtefXp/bzG3DP/37E99Od2IvLjcJ\nDQ31dgiaSJzG8YcYmSVOo/lanCaTiZccXMK1p9Tmp35/iuNuxDGz6+PkHXXttgdwkpnPMnM6gCUA\nBuU65mEAC7Oy+E4A5YnI6vCXcWvHod60egiPDcdPQ35C5EuReLrl07o3l3CnsLAwb4egicRpHH+I\nEZA4jeZrcRIRhjUbhqOvHkXdcnXRYk4LfLHd9p4EWjlK8rUAxFj8HJv1mKNjalu7WMkiJbH/lf1Y\nOnQpOgV2crgGu95/BEfHW3ve1X9oZ863d46t5/whTiNeNBKna/dw9RreeA05c438/FqPDI/Ep70+\nRcSLEdgavVX3PXNzlOS19srmztZWz5vUexICywVqvKT8w7t6H1fP94ek5OgciVPfNSTJO39Pvec4\neq03qNgAfw7/U/c9c7Pb8UpEHQEEM3O/rJ8nADCxRecrEX0LIIyZl2T9fBRAd2a+kOta3tngVQgh\n/By70PHqaCn93QAaElE9APEAhgHIvXfaagBjACzJelO4ljvBuxqkEEII59hN8sycQURjAKyHGmkz\nn5mjiGhU1vNzmXkNEfUnopMAbgJ4zu1RCyGE0MRjk6GEEEJ4nt+sQimEEEI/SfJCCJGPeS3JE1Ep\nIlpIRPOI6ElvxeEIEd1FRN8T0W/ejsUWIhqU9XtcQkQPeDseW4joHiKaQ0TLiOgFb8djT9bfZyQR\nDfB2LLYQURARbc36nXb3djy2kPIpEc0gopHejscWIuqS9bv8joi2ezsea4ioNhH9QUTziWi8lnO8\n2ZJ/BMAyZn4ZatasT2LmM8z8orfjsIeZV2X9Hl+BGgHlk5j5KDOPBvAEgL7ejseBtwEs9XYQDpgA\nJAEoBjUJ0VcNhpo0mQYfjpOZt2X9ff4F4Ecvh2NLcwC/M/MLAO7TcoKhSZ6IFhDRBSI6mOvxfkR0\nlIhOWLz7WM6UzTQyDoPj9AonY3wfwCzPRak/TiIaCOBvqCUyfDLOrE9DRwAkejJGvXEC2MrM/QG8\nA+BDH46zEYDtzPwmgNE+HKfZkwB+9dEYwwG8TESbAKzTdANXFr7J/QWgK9S7y0GLx6wucgZgBIAB\nWccsNjIOI+O0eP43X40RasbxZAC9PBmjs7/LrGNW+WqcAD4BMBVq6PBKZI1C87U4LZ4v6uN/n08B\neCzrmKW+GmfWc3UAzPPVGAG8DqBr1jGa/s0dTYbShZm3Zk2cspS9yBkAEJF5kbMZAGZl1TxXGxmH\nI3riJKILAD4D0IqIxrOVpZa9HSOA3gB6AShLRA2Yea4nYtQbJxFVhSrTFQcQ6qkYAX1xMvP7WT8/\nAyCRs15RvhYnEd0DVfYqD2Cmp2IEdP99Tgcwk4i6AgjzXJS644wC8DyABR4MUW+MawD8L6sf84yW\n6xua5G2wtoBZB2a+BfUL9RW24rwCVev2BbZifA0efpE7YCvOLQC2eCckq6zGaf6BmRd6PCLrbP0+\nPwewwjshWWUrzhQAvtSvZfPfnZmDvRGQFbZ+lwcADNVzIU90vPrLbCt/iNMfYgQkTqNJnMbyhzgN\ni9ETST4OgOXSk4HwzR52f4jTH2IEJE6jSZzG8oc4DYvRE0k+e5EzIioKNcTPozV4jfwhTn+IEZA4\njSZxGssf4jQuRoN7iRdDrVaZClVPei7r8QcBHIPqLZ7gyZ5rf43TH2KUOCVOidP3Y5QFyoQQIh+T\ntWuEECIfkyQvhBD5mCR5IYTIxyTJCyFEPiZJXggh8jFJ8kIIkY9JkhdCiHxMkrwQQuRjkuSFECIf\n+39Vs+ydqSMVCgAAAABJRU5ErkJggg==\n", "text": [ "" ] } ], "prompt_number": 9 }, { "cell_type": "code", "collapsed": false, "input": [ "def buildgrid(densities=np.logspace(0,8,20), abundance=10**-8.5, fortho=1e-3, **kwargs):\n", "\n", " tau = np.zeros_like(densities)\n", "\n", " R = pyradex.Radex(species='co', abundance=abundance, collider_densities={'oH2':10,'pH2':10},\n", " temperature=20)\n", " for jj,dd in enumerate(densities):\n", " R.density = {'oH2':dd*fortho,'pH2':dd*(1-fortho)}\n", " R.abundance = abundance # reset column to the appropriate value\n", " R.run_radex(**kwargs)\n", " tau[jj] = R.tau[0]\n", "\n", " return tau\n" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 22 }, { "cell_type": "code", "collapsed": false, "input": [ "tau1 = buildgrid()" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 23 }, { "cell_type": "code", "collapsed": false, "input": [ "pl.loglog(np.logspace(0,8,20), tau1)" ], "language": "python", "metadata": {}, "outputs": [ { "metadata": {}, "output_type": "pyout", "prompt_number": 24, "text": [ "[]" ] }, { "metadata": {}, "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAX8AAAEHCAYAAABGNUbLAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAHfRJREFUeJzt3XmUVeWV9/HvFoMhKE5rsWyxDLHVBgNGEwWzBCmDNhgT\nsR2DL8YGwcQEHN+ARlvL9nXAeUDjwBCkW0DbkWhEY1uIGBWcmA0YiQwJJgEjKAGq2O8fz61wLauo\nO5x7zzn3/j5rseCcqnvupig2m32esx9zd0REpLrsFHcAIiJSfkr+IiJVSMlfRKQKKfmLiFQhJX8R\nkSqk5C8iUoWU/EVEqpCSv4hIFdq5lBc3s0HAiUAnYIK7v1DK9xMRkdxYOZ7wNbM9gFvcfXjJ30xE\nRNqUd9vHzCaa2VozW9Ds/EAzW2pmy8xsTLOXXQmMKyZQERGJTiE9/0nAwOwTZtaOkNwHAocAg82s\nuwVjgV+7+ztFRysiIpHIu+fv7rPNrGuz072A5e6+AsDMpgGDgOOA/kAnMzvQ3e8vKloREYlEVDd8\nuwArs45XAb3dfRRw945eaGYaKyoiUgB3t0JfG9VSz6ISuLvn9ePqq6+O9DUtfaz5uajfM81xtvX5\n1RxnOf7MFWe0cab173qxokr+q4GarOMaQvWfk7q6Ourr63N+s9ra2pw/N5fXtPSxQt4jimukIc62\nPr+a4yzHn3mh71Ps6ys1zrT9Xa+vr6euri7v9/uCfP+1yvyL0xVYkHW8M/B+5nx74B2ge47X8jS4\n+uqr4w4hJ4ozWoozWmmIMw0xurtncmdBOdzdC1rqORV4FTjYzFaa2VB3bwBGAjOBxcB0d19S5L9L\niRJFdVAOijNaijNaaYgzDTFGoSwPee0wADOPOwYRkbQxMzwBN3yLkm/PX0SkWkXV81flLyKSQhVR\n+YuISHmVdKqniIhEq7ERbrut+Oso+YuIpMSSJTB0KHzlK8VfKxFtH93wFRFpXWMj3HQTHHMMfPvb\n9fTtW1f0NXXDV0QkwZqq/Y4dYcIE6No1nNcNXxGRCpRd7Z9zDrzwwvbEHwX1/EVEEia72p87N9qk\n30SVv4hIQpS62s+myl9EJAHKUe1nU+UvIhKjclb72VT5i4jEpNzVfrZEVP5a5y8i1aSYal+D3URE\nUqi1dfv50jp/EZEUaKr2+/Ytb2+/Ner5i4iUWHa1P29evEm/SUkrfzP7mpmNN7NHS/k+IiJJFNdK\nnlyUtPJ39w+A4Ur+IlJt4lzJk4tCNnCfaGZrzWxBs/MDzWypmS0zszHRhSgikh5JrvazFVL5TwLu\nBh5qOmFm7YBxwHHAamCumT3t7ksiiVJEJAWSXu1ny7vyd/fZwPpmp3sBy919hbtvBaYBg8xsLzO7\nDzhM/xsQkUqVlmo/W1Q9/y7AyqzjVUBvd18H/LitF2c/sFBbW0ttbW1EYYmIlFa5qv36+vpIH4Yt\n6CEvM+sKzHD3npnjU4GB7j4iczyEkPxH5XAtPeQlIqnT2Ai33go33wz/+Z/wox/BTmV8cqrYh7yi\nqvxXAzVZxzWE6j8ndXV1qvhFJDXi7O1H9T+AqCr/nYH3gP7AGuANYHAuN3xV+YtIWsRd7Wcre+Vv\nZlOBfsDeZrYSuMrdJ5nZSGAm0A6YoJU+IlJJ0rSSJxca7CYisgMNDdur/Wuvjbfaz5aUnn9R1PMX\nkSRasgT+/d9h112TM5Mn1p5/lFT5i0jSJLXaz1YRlb+ISFIksdovhYT9WyYiEo+GBhg7NjylO3Qo\n/OY3lZv4QZW/iAiLF4eEv+uulbGSJxeJqPy1h6+IxKGp2u/XLz3VvvbwFREpQna1X8xeunHRHr4i\nInlIY7VfCur5i0jVaKr2d9utenr7rVHlLyIVL7vaHzYsHfP2S02Vv4hUNFX7LVPlLyIVqaEBbrxx\ne29f1f7nqfIXkYqzaFFI+J06qdpvTSIqf63zF5EoNDTADTeEav/ccyuz2tc6fxGRLE3V/u67w/jx\n8NWvxh1RaWmdv4hUtebV/vPPV37ij4J6/iKSWosWhQmce+wBb76ppJ+Pklb+ZtbRzCab2QNmdlYp\n30tEqkdDA1x/PdTWwogRqvYLUerK/xTgEXd/xsymAQ+X+P1EpMItXBiq/T33DPP2lfQLk3flb2YT\nzWytmS1odn6gmS01s2VmNiZzuguwMvPrxiJjFZEq1lTtH3ts2FlL1X5xCmn7TAIGZp8ws3bAuMz5\nQ4DBZtYdWAXUFPFeIiIsXAhHHQX19aG3P2IEWMHrXAQKSMjuPhtY3+x0L2C5u69w963ANGAQ8Dhw\nqpndCzxdbLAiUl0aGuC667ZX+zNnwv77xx1VZYiq55/d3oFQ8fd298+AYW29OPuBhdraWmprayMK\nS0TSasGCsG5/r71CtV/tSb++vj7Sh2ELesjLzLoCM9y9Z+b4VGCgu4/IHA8hJP9ROVxLD3mJyD9s\n3RomcN55Z1i/f+65avG0pNiHvKKq/FezvbdP5tercn1xXV2dKn4RYf78sJKnc2dV+62J6n8AUVX+\nOwPvAf2BNcAbwGB3X5LDtVT5i1S5rVtDlX/33WES57BhqvbbUvbK38ymAv2Avc1sJXCVu08ys5HA\nTKAdMCGXxC8i8u67odrfZx946y2oqWnzJRIBDXYTkVhs2RKq/XHj4Kabwj8AqvZzl5Sef1HU8xep\nLu+8E5J9ly7w9tuw335xR5Qesfb8o6TKX6R6bNkS1u3/4hdw883wwx+q2i9URVT+IlL53n47VPs1\nNeHXXbrEHVF108gFESmpLVvgqqtgwAC49FKYMUOJPwlU+YtIybz5ZnhK96tfDX3+ffeNOyJpkojK\nX3v4ilSWzZvhyivhhBNg9Gh4+mkl/qhoD18RSaS5c0O1f8ABcP/98E//FHdElUl7+IpIIvz973D5\n5fC978HPfw5PPaXEn2Tq+YtI0V57LYxk6N49PLG7zz5xRyRtUfIXkYJt2hRW8kyZAnfdBaefrnX7\naaG2j4gUZM4cOOww+PDDMHv/jDOU+NNElb+I5OWzz+CKK2DatDCX59RT445ICqHKX0Ry9vLL8I1v\nwEcfhWpfiT+9ElH5a7CbSLJ9+mlYyfPYY3DvvTBoUNwRVS8NdhORsnjppbCVYp8+cMcdYU9diZ8G\nu4lISWzYAGPGhKdz77svrN+XyqGev4h8wW9+Az17hge3Fi5U4q9EJa38zexrwBXA7u5+einfS0SK\n97e/wc9+Bs89Bw88AAMHxh2RlEpJK393/8Ddh5fyPUQkGs8+G6p9CCt5lPgrW07J38wmmtlaM1vQ\n7PxAM1tqZsvMbExpQhSRUlq3Ds45B0aOhEmTQsW/++5xRyWllmvlPwn4XB1gZu2AcZnzhwCDzay7\nmZ1tZrebmQa4iiTcE09Ajx6wxx4wfz707x93RFIuOfX83X22mXVtdroXsNzdVwCY2TRgkLvfCEzJ\nnNsLuB44zMzGuPvYlq6fPZta6/1FSu+jj2DUqLDByiOPhGWckmxRre9vkvM6/0zyn+HuPTPHpwED\n3H1E5ngI0NvdR+UVgNb5i5SNexjLcPHFodVTVwcdOsQdlRQiznX+ytgiKbJmDZx/Prz/fthH98gj\n445I4lTMap/VQE3WcQ2wqpALaRtHkdJxh4kTw0yeww4L++oq8adX2bdxbKHtszPwHtAfWAO8AQx2\n9yV5BaC2j0jJ/OEPMGIE/OUvYSXPN74Rd0QSlbJs42hmU4FXgYPNbKWZDXX3BmAkMBNYDEzPN/GL\nSGls2xYGsB1xBBx7LLz+uhK/fJ4Gu4lUmOXLYfhw2Lw5tHu6d487IimFitjAXT1/keI1NsLtt8NR\nR8HJJ8MrryjxV6Ky9/xLRZW/SPEWLw5jl3fZBcaPhwMPjDsiKbWKqPxFpDBbt8K110K/fmHd/v/+\nrxK/5Ebz/EVS6s03Ydgw2G8/eOstqKlp+zUiTVT5i6TMpk1w2WXw3e+G8cu/+pUSv+QvEZW/9vAV\nyc3s2aG3f/jhYexy585xRyTlpj18RarIhg2h2n/ySbjnnrCaR6qbbviKVLiZMz+/paISv0QhEW0f\nEfmidevgkktg1ix48EE4/vi4I5JKospfJIEeeyxssrL77qG3r8QvUVPlL5Igf/oT/PSnsGgRPPoo\nHH103BFJpVLlL5IA7jB5Mhx6KHTrFnbYUuKXUlLlLxKzDz+E886DtWvDzd3DD487IqkGiaj8NdhN\nqpE73H8/fOtbcMwx8MYbSvzSNg12E0mxDz4IY5c3bgxjl7/+9bgjkrTROn+RFNm2LTykdeSRMGAA\nzJmjxC/xUM9fpEyaNlnZsiXM2u/WLe6IpJqVvPI3s0Fm9oCZTTMzrVaWqtPYCHfcETZZOemkMJ9H\niV/iVraev5ntAdzi7sObnVfPXyrWe++Fscs77RR6+wcdFHdEUinK1vM3s4lmttbMFjQ7P9DMlprZ\nMjMbs4NLXAmMKzRQkTRpbIRbbglr9X/wgzCiQYlfkiTnyt/M+gIbgYfcvWfmXDvgPeA4YDUwFxgM\nHAF8E7gZ+CNwI/C8u7/YwnVV+UtFWbw4VPsdOsCECXDAAXFHJJWo2Mo/5xu+7j7bzLo2O90LWO7u\nKzLBTAMGufuNwJTMuQuA/kAnMzvQ3e9vfu3sNaua6y9p1dAAN98Mt94atlb80Y9Cu0ckClHN8W+S\nV88/k/xnZFX+pwED3H1E5ngI0NvdR+VxTVX+knoLF8LQobDHHmECZ9eucUcklS7udf7K2lLVmjZQ\nP/bYMKLh+eeV+CUdil3nvxrI3j20BliV70W0jaOk0SuvwMiRsM8+YTP1/fePOyKpBrFs49hC22dn\nwg3f/sAa4A1gsLsvyeOaavtIqqxeDWPGQH093HQTDB4MVvB/vkUKU86lnlOBV4GDzWylmQ119wZg\nJDATWAxMzyfxi6TJ5s1w441h7PL++8PSpXDWWUr8kk4a7CaSg2eegYsugu7d4bbb4MAD445Iql3Z\nlnqWknr+klTLlsHFF8Pvfgd33QUnnBB3RFLtYun5l4Iqf0mijRvhuuvCss3Ro0PV37593FGJbBf3\nUk+RiuIODz8cBq+tXg3z54fkr8QvlSYRbR+RJHjnHRg1Cj77DKZP1x66UtkSUflrG0eJ01//Cuef\nHzZXOfvssJ2iEr8klbZxFClSY2PYQ7euLkzevOYa2HPPuKMSyU1FrPYRKZcNG2DRotDLv/fekOxf\nfBF69ow7MpHyUuUvFWnLlrCRyoIFYeha088ffRTW6vfsCd/7Hpxyih7SknQqtvJX8pdU27YNVqz4\nYpJ///0wYK1nT+jRY/vPBxwA7drFHbVI8ZT8pSK4h/EJn33W9o/168OGKQsXhp/33POLSb5bN/jy\nl+P+XYmUjpK/JNq2bSFJz5oVNi7/4x9DAt+06YtJvX17+MpX2v7RqVNo3fToEX7svnvcv0uR8lPy\nl0RpbAzr5V9+eXvC32sv6NcPjjkmtGJaSugdOqgdI5IPJX+J1dat8NZbIdHPmgVz5sC++4Zk35Tw\n99037ihFKk9FLPXUYLf02LwZ5s3bnux/+9tQzffrF7YxnDQJOneOO0qRyqXBblI2v/sdTJ0akv3c\nuXDwwdsr+759Q1tHRMpLbR8pCffQt7/1VnjttTD24DvfgT59dINVJAkqou0jybF1K/zP/4Sk/8kn\ncMklMG1auCkrIpWjpJW/mXUDLgT2Bma6+4QWPkeVfwJ88gmMHw933hl6+JdeGp6A3SkRo/9EpLlU\ntH3MbCdgmruf0cLHlPxj9OGHYYeqSZPg+OND0j/yyLijEpG2lGUzFzObaGZrzWxBs/MDzWypmS0z\nszGtvPb7wDPAtEKDlOi9+WbYfPzww0N//623QntHiV+kOuRU+ZtZX2Aj8JC798ycawe8BxwHrAbm\nAoOBI4BvAje7+5qsazzl7oNauLYq/zLZti1sRH7rrfD738MFF8CIEbqBK5JGZbnh6+6zzaxrs9O9\ngOXuviITyDRgkLvfCEzJnOsHnAJ8GXiptetnb0yg9f7R27QJHnoIbr8dOnYMrZ3TT4cvfSnuyEQk\nV1Gt72+Sc88/k/xnZFX+pwED3H1E5ngI0NvdR+UVgCr/ktm4EW67De65B3r1Ckm/Xz+NMBapBHFu\n4K6MnVDbtsEvfwn/8i+wbFl4OGvGDKitVeIXkaCYdf6rgZqs4xpgVSEX0niH6MyeDRddBLvsAk88\nESp+EakcZR/v0ELbZ2fCDd/+wBrgDWCwuy/JKwC1fSLxwQcwZgy8/jqMHQtnnqkqX6SSlWup51Tg\nVeBgM1tpZkPdvQEYCcwEFgPT8038UrwNG+DnPw9LNA89FJYuDZuRK/GLyI7kutpncCvnfw38utgg\n1PbJX2MjTJ4MV14J//qv8O670KVL3FGJSKlpqmcVe/nl0Nfv0AHuuEMPZolUIw12qyK//z2MHh3G\nKt90E5xxhto7IlIYje1KgU8+gcsuCxX+4YeHvr5u6IpIMZT8E6yxESZMgG7dYO1aWLAArrgitHtE\nRIqRiLaPbvh+0ezZYfZOx47w9NNwxBFxRyQiSaAbvhVq3brQ13/uuTCATX19EWlJnOMdJELu8N//\nDV//emjrLF6svr6IlE4i2j7V7v334Sc/gT/9CZ58Enr3jjsiEal0qvxjtHUr3HBDSPbHHQfz5inx\ni0h5qPKPyW9/C+edB/vtF9btf+1rcUckItVEyb/MPv4YLr8cnnoqbK6iG7oiEge1fcrEHR59NNzQ\ndYdFi3RDV0Tio8q/DP7wB/jpT8PY5UcegaOPjjsiEal2iaj86+rqIt2bMikaGsJa/W99C779bXj7\nbSV+ESlOfX395/Y9L5Qe8iqRuXPDDd2994b77oMDD4w7IhGpJHrIK2E2bIALL4Tvfx8uuQReeEGJ\nX0SSp+TJ38w6mtlcMzux1O8Vt+efh5494W9/Czd0zz5bN3RFJJnKccN3NDC9DO8Tm/XrQ5X/0ktw\n//0wYEDcEYmI7Fiue/hONLO1Zrag2fmBZrbUzJaZ2ZgWXnc8YX/fP0cTbvI8/jj06AG77hpGLivx\ni0ga5HTD18z6AhuBh9y9Z+ZcO+A94DhgNTAXGAwcAXwTuBn4CdAROATYBPxb87u7ab3hu3YtjBwJ\n8+eHmft9+sQdkYhUk7Lc8HX32cD6Zqd7AcvdfYW7bwWmAYPcfYq7X+zua9z9Sne/GHgYeCCVWb4Z\nd5gyBQ49FP75n+Gdd5T4RSR9iun5dwFWZh2vAlocS+buk3d0oew1q0ne1OXDD+HHP4bVq+HZZ8P6\nfRGRcohqE5cmOa/zN7OuwIysts+pwEB3H5E5HgL0dvdReQWQgrbPtm3hRu5VV4VlnKNHQ/v2cUcl\nItWs2LZPMZX/aqAm67iGUP1XlGXLYPhw2LwZZs2CQw6JOyIRkeIVs85/HnCQmXU1s/bAmcDThVwo\nieMdGhrgllvCWIaTT4Y5c5T4RSR+ZR3vYGZTgX7A3sBHwFXuPsnMTgDuANoBE9z9hrwDSGDbZ8EC\nOPdc2G03ePBBOOCAuCMSEfm8Yts+mu2TZcsWuP56uOee8PPw4XpCV0SSKc6ef2Tq6upiX+Uzbx4M\nHQpdu4bpm/vtF1soIiKtimrVT9VX/n//O1xzDUycCLfdBmedpWpfRJKvIir/uLz2GgwbBt27w7vv\nwj77xB2RiEh5VGXy37QJ/uM/4L/+C+66C04/XdW+iFSXqpvnP2cOHHYYrFwZVvVoA3URqUaJqPzL\nccP300/hiivCHrrjxsEpp5TsrURESkY3fPMwa1ZYt3/UUXDnnWFrRRGRNNMN3x3YuBEuuwyefBLu\nvRdOOinuiEREkqFie/4vvhi2VPz009DbV+IXEdmu4ir/Tz6Bn/0sjFx+4AE44YS4IxIRSZ6Kqvxn\nzgzV/rZtsHChEr+ISGsqovL/+GO49NLQ6hk/Ho4/Pu6IRESSLfWV/7PPhmq/ffvQ21fiFxFpW2or\n//Xr4eKLwzLOyZPhO9+JOyIRkfRIROWf72YuzzwTqv2OHUO1r8QvItWirJu5lFI+D3mtXw8XXQSv\nvBJ6+8ceW+LgREQSqtiHvBJR+edixgzo0QM6dQoTOJX4RUQKV9Kev5nVAtcCC4Fp7j4r32usWwcX\nXgivvgoPPwz9+kUdpYhI9Sl15b8N2ADsAqzK98VPPRV6+3vtBfPnK/GLiEQl1w3cJwInAh+5e8+s\n8wPZvoH7eHcf2+x15u5uZp2B29x9SAvX/kLP/69/hQsugNdfDztsHXNMIb81EZHKVa6e/yRgYLM3\nbgeMy5w/BBhsZt3N7Gwzu93M9s3K6h8Tqv82PflkqPY7dw7VvhK/iEj0cur5u/tsM+va7HQvYLm7\nrwAws2nAIHe/EZiSOfdvwABgD+DuHb3HX/4Sqv25c8PM/T598vltiIhIPoq54dsFWJl1vAronf0J\n7v4E8ERbFzrzzLp/rN2/++5a+vSpLSIsEZHKE9UmLk1yXuefqfxnNPX8zexUYKC7j8gcDwF6u/uo\nvAIw84MOciZNgqOPzueVIiLVK87NXFYDNVnHNRSwogfCuv0OHYqIRERE8lLMUs95wEFm1tXM2gNn\nAk8XcqGxY/Mb7yAiUq3KOt7BzKYC/YC9gY+Aq9x9kpmdwPalnhPc/Ya8AyjDHr4iIpWm2LZPqmb7\niIhIUBGzffKd6ikiUq2qcqqniIgEFVH5i4hIeSn5i4hUISV/EZEqlIjkrxu+IiK50Q1fEZEqphu+\nIiKSNyV/EZEqpOQvIlKFlPxFRKqQkr+ISBVS8hcRqUJK/iIiVSgRyV8PeYmI5EYPeYmIVDE95CUi\nInkrZgP3NpmZAf8P2A2Y5+4PlfL9REQkN6Wu/E8GugBbgFUlfq+SSss9CcUZLcUZrTTEmYYYo5BT\n8jeziWa21swWNDs/0MyWmtkyMxvTwksPBua4+/8Fzo8g3tik5RtCcUZLcUYrDXGmIcYo5Fr5TwIG\nZp8ws3bAuMz5Q4DBZtbdzM42s9vNbF9Ctf9x5iXbIoq5oD+cHb2mpY9F8Q1QqXG29fnVHGc5/swL\nfZ9iX1+pcVbT3/VsOSV/d58NrG92uhew3N1XuPtWYBowyN2nuPvF7r4GeBwYYGZ3AfVRBZ2WL3Sl\nxqnkX/h75vuaNCfVQt+n2Nen4e9QW68pR/LPeamnmXUFZrh7z8zxacAAdx+ROR4C9Hb3UXkFYKZ1\nniIiBShmqWcxq30iSdrFBC8iIoUpZrXPaqAm67iGlK/oERGpFsUk/3nAQWbW1czaA2cCT0cTloiI\nlFKuSz2nAq8CB5vZSjMb6u4NwEhgJrAYmO7uS0oXqoiIRCX22T4iIlJ+mu0jIlKFEpf8zayjmU02\nswfM7Ky442mNmX3NzMab2aNxx7IjZjYo87WcZmbHxx1Pa8ysm5n9wsweMbNz446nNZnvz7lmdmLc\nsbTGzGrNbHbm69kv7nhaY8F1ZnaXmf0w7nhaY2Z9Ml/LB81sTtzxtMbM9jOzx81sQisTFz4ncckf\nOAV4xN3PA06KO5jWuPsH7j487jja4u5PZb6WPybclE8kd1/q7ucDPwAGxB3PDowGpscdRBu2ARuA\nXUj2CrxUzP5y91cy35u/An4Zczg70hN4zN3PBQ5v65PLkvzznA3UBViZ+XVjOeIrMM7YFBjnlYRx\nHGWTb5xm9n3gGcLT4omLMfM/p8XAn8sVXyFxArPd/bvAZcA1CY4zttlfBf4dOgt4uHxR5h3nq8B5\nZvYi8FybF3f3kv8A+hL+JVqQda4dsBzoCnwJeAfoDgwBTsx8ztRyxFdInFkff7ScMRbw9TRgLNA/\nyXE2e91TSYyRMJ78dsIKtyfJLJhIWpxZH29f7u/PPL+e/wc4PfM505MaZ+Zj+wMPlDPGAr6eFwF9\nM5/T5p97Sef5N3H32ZnxENn+MRsIwMymAYOAu4BxmZ5qWZ8byCdOM1sLXA8cZmZj3H1sEuMEjgP6\nA53M7EB3vz+JcZpZZ0LL78vAS0mM0d2vzByfA/zZM3/LkhanmXUjtM72AO4uV4yQ9/fmncDdZtaX\nCGd/5SLPOJcAw4CJZQwRyDvOZ4GrMvdKP2jr2mVJ/q3Ibu9A6Pn1dvfPCF/opGgtznWEPnpStBbn\nKMqcANrQWpyzgFnxhPQFLcbYdODuk8seUcta+1reCDwRT0gtai3OTUCS7pu1+ufu7nVxBNSK1r6e\n84HTcr1InDd80/KAgeKMVhriTEOMoDijVlVxxpn80zIbSHFGKw1xpiFGUJxRq6o440z+aZkNpDij\nlYY40xAjKM6oVVecZbpjPRVYA2wm9KqGZs6fALxHuHN9ebnvpCtOxZmGGBWn4ixFnJrtIyJShZL4\nhK+IiJSYkr+ISBVS8hcRqUJK/iIiVUjJX0SkCin5i4hUISV/EZEqpOQvIlKF/j+U2li2Wq0a3QAA\nAABJRU5ErkJggg==\n", "text": [ "" ] } ], "prompt_number": 24 }, { "cell_type": "code", "collapsed": false, "input": [], "language": "python", "metadata": {}, "outputs": [] } ], "metadata": {} } ] }