{ "metadata": { "name": "development-spectrum-test" }, "nbformat": 3, "nbformat_minor": 0, "worksheets": [ { "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# QuTiP development notebook for testing spectrum functions\n", "\n", "Copyright (C) 2011 and later, Paul D. Nation & Robert J. Johansson" ] }, { "cell_type": "code", "collapsed": false, "input": [ "%pylab inline" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 1 }, { "cell_type": "code", "collapsed": false, "input": [ "from qutip import *" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 2 }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Parameters" ] }, { "cell_type": "code", "collapsed": false, "input": [ "N = 4 # number of cavity fock states\n", "wc = wa = 1.0 * 2 * pi # cavity and atom frequency\n", "g = 0.1 * 2 * pi # coupling strength\n", "kappa = 0.75 # cavity dissipation rate\n", "gamma = 0.25 # atom dissipation rate\n", "\n", "# Jaynes-Cummings Hamiltonian\n", "a = tensor(destroy(N), qeye(2))\n", "sm = tensor(qeye(N), destroy(2))\n", "H = wc * a.dag() * a + wa * sm.dag() * sm + g * (a.dag() * sm + a * sm.dag())\n", "\n", "# collapse operators\n", "n_th = 0.01\n", "c_ops = [sqrt(kappa * (1 + n_th)) * a, sqrt(kappa * n_th) * a.dag(), sqrt(gamma) * sm]" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 3 }, { "cell_type": "code", "collapsed": false, "input": [ "# calculate the correlation function using the mesolve\n", "# solver, and then fft to obtain the power spectrum\n", "tlist = linspace(0, 100, 2500)\n", "corr = correlation_ss(H, tlist, c_ops, a.dag(), a)\n", "wlist1, spec1 = spectrum_correlation_fft(tlist, corr)\n", "\n", "# calculate the power spectrum using spectrum_ss, which\n", "# uses essolve to solve for the dynamics\n", "wlist2 = linspace(0.25, 1.75, 500) * 2 * pi\n", "spec2 = spectrum_ss(H, wlist2, c_ops, a.dag(), a)" ], "language": "python", "metadata": {}, "outputs": [], "prompt_number": 4 }, { "cell_type": "code", "collapsed": false, "input": [ "# plot the correlation function\n", "fig, ax = plt.subplots(1, 1, figsize=(10,4))\n", "\n", "ax.plot(tlist, real(corr), 'b', lw=1, label='real')\n", "ax.plot(tlist, imag(corr), 'r', lw=1, label='imag')\n", "\n", "ax.legend(loc=1)\n", "ax.set_xlabel('Time')\n", "ax.set_ylabel('Correlation')" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "pyout", "prompt_number": 5, "text": [ "" ] }, { "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAnkAAAEKCAYAAABuV0NTAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3Xl8VfWd//HXOXdNQkJCIAESNk2QTZAB3NpqFAOCNaJU\npJ1RfspUastYdaza6bSFWhW0dlFaR1qLuJRi61SoQipW0dYpRCutCyhBQLIQtiRAlruf3x/n5iYh\nG0sWyH0/H4/7gHtzvud+7w2278fnuxmWZVmIiIiISK9i9nQHRERERKTzKeSJiIiI9EIKeSIiIiK9\nkEKeiIiISC+kkCciIiLSCynkiYiIiPRCPRLyCgsLGTVqFLm5uSxdurTVa26//XZyc3OZMGECW7Zs\n6bBtUVER559/PhMnTmTKlCm88847Xf45RERERE5X3R7ywuEwCxcupLCwkK1bt7Jq1Sq2bdvW7Jp1\n69axY8cOiouLWb58ObfddluHbe+55x7uv/9+tmzZwg9+8APuueee7v5oIiIiIqeNbg95RUVF5OTk\nMHz4cFwuF3PnzmXNmjXNrlm7di3z5s0D4IILLqC6upqKiop22w4aNIjDhw8DUF1dTVZWVvd+MBER\nEZHTiLO737CsrIwhQ4bEnmdnZ7N58+YOrykrK6O8vLzNtkuWLOHzn/88d999N5FIhL/97W8t3tsw\njM7+OCIiIiJd5lQOJuv2St7xBq0T/VDz58/nscceY8+ePfzkJz/hlltuafO+epx5j+9///s93gc9\n9PuLx4d+d2f2Q7+/M/txqro95GVlZVFSUhJ7XlJSQnZ2drvXlJaWkp2d3W7boqIirr32WgC+9KUv\nUVRU1JUfQ0REROS01u0hb/LkyRQXF7N7924CgQCrV6+moKCg2TUFBQU888wzAGzatInU1FQyMzPb\nbZuTk8Obb74JwOuvv87IkSNbff/otD0RERGRXq3b5+Q5nU6WLVvG9OnTCYfDzJ8/n9GjR/Pkk08C\nsGDBAmbOnMm6devIyckhKSmJFStWtNsWYPny5XzjG9/A7/eTkJDA8uXLW33/I0egb9/u+azSefLy\n8nq6C3IK9Ps7c+l3d2bT7y++GVZnDPqeIQzD4JNPLNoo8omIiIicNgzDOKW5ed1eyetp9fU93QMR\nEZH40q9fP6qqqnq6G6ettLQ0KisrO/2+cRfy6up6ugciIiLxpaqqqlNWi/ZWXbXFW9ydXatKnoiI\niMQDhTwRERGRXijuQp7f39M9EBEREel6cRfywuGe7oGIiIj0Nhs3bmx29OrpQCFPREREpBeKu5AX\nCvV0D0REROR0FOplISHuQp4qeSIiItJg+PDhPPzww4wfP57k5GTefvttLr74YtLS0jjvvPNiR6YC\nrFixgjFjxpCSksLZZ5/d5ulap4u42yevl4V0EREROUW//e1vWb9+PYZhMH78eJ577jmuvPJKXnvt\nNWbPns0nn3xCeno6mZmZvPLKK4wYMYK33nqLGTNmMGXKFCZOnNjTH6FVquSJiIhIjzOMznmc+Psa\n3H777WRlZfHss88yc+ZMrrzySgCuuOIKJk+ezCuvvALAzJkzGTFiBACXXHIJ06ZN4y9/+UunfQed\nLe5Cnip5IiIipx/L6pzHyWhYFfvZZ5/xu9/9jrS0tNjj7bffpqKiAoD169dz4YUXkp6eTlpaGuvW\nrePQoUOd9RV0urgbrk0t+QA+S4Fhw3q6KyIiInIaaDhWbOjQodx4442tzrXz+/3Mnj2b5557jmuu\nuQaHw8G11157Wh/XFneVvK8sGQ8FBT3dDRERETnN/Nu//Rt//OMfefXVVwmHw/h8PjZu3EhZWRmB\nQIBAIED//v0xTZP169fz6quv9nSX2xV3IU9ERESkNdnZ2axZs4YHH3yQjIwMhg4dyqOPPoplWSQn\nJ/PYY48xZ84c+vXrx6pVq7jmmmuatTdOZlJgFzKs07nO2MkMw8ACOP982Ly5p7sjIiISFwzDOK2H\nNXtaW9/PqX5v8VnJ0+oLERER6eXiM+T5fD3dAxEREZEuFZ8hTyVjERER6eUU8kRERER6IYU8ERER\nkV5IIU9ERESkF1LIExEREemFeiTkFRYWMmrUKHJzc1m6dGmr19x+++3k5uYyYcIEtmzZclxtH3/8\ncUaPHs24ceO499572+6AQp6IiIj0ct0e8sLhMAsXLqSwsJCtW7eyatUqtm3b1uyadevWsWPHDoqL\ni1m+fDm33XZbh23feOMN1q5dy/vvv8+HH37I3Xff3XYnFPJEREQEGDduHG+99VZPd6NLOLv7DYuK\nisjJyWH48OEAzJ07lzVr1jB69OjYNWvXrmXevHkAXHDBBVRXV1NRUcGuXbvabPvEE0/w7W9/G5fL\nBcCAAQNaff9FAJWVsGgReXl55OXldc0HFRERkdPehx9+2NNdiNm4cSMbN27stPt1e8grKytjyJAh\nsefZ2dlsPuaIsdauKSsro7y8vM22xcXFvPXWW/zXf/0XXq+XH/3oR0yePLnF+y8CSE2FRYs682OJ\niIiInJJji0+LFy8+pft1+3Dt8R7ee6JntYVCIaqqqti0aROPPPIIc+bMae/mJ3RvERER6Z2GDx/O\nn//8ZxYtWsT111/PjTfeSEpKCuPHj6e4uJiHHnqIzMxMhg0bxoYNG2LtVqxYwZgxY0hJSeHss89m\n+fLlze778MMPM3jwYLKzs/nVr36FaZrs3LmzWz9bt4e8rKwsSkpKYs9LSkrIzs5u95rS0lKys7Pb\nbZudnc11110HwJQpUzBNk0OHDrXeiXC4sz6OiIiInMGaFp9efvllbrrpJqqqqpg4cSL5+fkAlJeX\n893vfpcFCxbErs3MzOSVV17hyJEjrFixgjvvvDO2ULSwsJCf/OQn/PnPf6a4uJiNGzced5GrM3V7\nyJs8eTLFxcXs3r2bQCDA6tWrKSgoaHZNQUEBzzzzDACbNm0iNTWVzMzMdtvOmjWL119/HYDt27cT\nCARIT09vvRMKeSIiIqcXw+icx0m/vcEll1xCfn4+DoeDL33pSxw6dIj77rsPh8PBDTfcwO7duzly\n5AgAM2fOZMSIEQBccsklTJs2jb/85S8AvPDCC9xyyy2MHj2ahIQEFi9efMIjlJ2h2+fkOZ1Oli1b\nxvTp0wmHw8yfP5/Ro0fz5JNPArBgwQJmzpzJunXryMnJISkpiRUrVrTbFuCWW27hlltu4dxzz8Xt\ndsdCYqsU8kRERE4vp8FUqoyMjNjfExIS6N+/f6wCl5CQAEBNTQ0pKSmsX7+exYsXU1xcTCQSoa6u\njvHjxwOwd+9ezj///Ni9jh2x7C7dHvIAZsyYwYwZM5q91rQECrBs2bLjbgvgcrl49tlnj68DCnki\nIiJykvx+P7Nnz+a5557jmmuuweFwcO2118aqdYMGDWoxvawnxN2JFxEMhTwRERFp5kSGUwOBAIFA\ngP79+2OaJuvXr+fVV1+N/XzOnDmsWLGCjz/+mLq6Ou6///6u6HKH4i7kBUyvQp6IiIjEGIYRexz7\nemvPk5OTeeyxx5gzZw79+vVj1apVXHPNNbHrrrzySm6//XYuu+wyRo4cyUUXXQSAx+Pp4k/SnGH1\nxEzAHmIYBoedaaS4/VBb29PdERERiQuGYfTIwoPTxbZt2zj33HMJBAKYZsv6Wlvfz6l+b3FXyQua\nHlXyREREpEv94Q9/wO/3U1VVxb333ktBQUGrAa8rxV3I8xsarhUREZGutXz5cjIzM8nJycHlcvHE\nE090ex/ibrh2d8I5DPNth0ikp7sjIiISF+J9uLYjGq7tJEHDbe/Fo5AnIiIivVjchbwQTnA4NGQr\nIiIivVr8hTzDBU6nQp6IiIj0aj1y4kVPUiVPRESke6WlpbXYc04apaWldcl94y7kRTAV8kRERLpR\nZWVlT3chLsXdcG0Yh0KeiIiI9HrxF/IshTwRERHp/eIu5MWGa0Ohnu6KiIiISJeJu5AX1pw8ERER\niQNxF/IilqktVERERKTXi7uQp0qeiIiIxIO4C3naQkVERETiQdyFvFDD6lotvBAREZFeLO5Cnip5\nIiIiEg/iLuSFLYU8ERER6f3iLuRpda2IiIjEgx4JeYWFhYwaNYrc3FyWLl3a6jW33347ubm5TJgw\ngS1bthx320cffRTTNNs8J0+ra0VERCQedHvIC4fDLFy4kMLCQrZu3cqqVavYtm1bs2vWrVvHjh07\nKC4uZvny5dx2223H1bakpIQNGzYwbNiwtt9fw7UiIiISB7o95BUVFZGTk8Pw4cNxuVzMnTuXNWvW\nNLtm7dq1zJs3D4ALLriA6upqKioqOmx711138fDDD7f7/hFLx5qJiIhI7+fs7jcsKytjyJAhsefZ\n2dls3ry5w2vKysooLy9vs+2aNWvIzs5m/Pjx7b7/C+EP2VZWBb/6FXnBIHl5eZ3wqUREREROzcaN\nG9m4cWOn3a/bQ55hGMd1nWVZx33P+vp6HnzwQTZs2NBh+2vMidw0YhfcdBMo4ImIiMhpIi8vr1nx\nafHixad0v24PeVlZWZSUlMSel5SUkJ2d3e41paWlZGdnEwwGW2376aefsnv3biZMmBC7ftKkSRQV\nFZGRkdHs3pqTJyIiIvGg2+fkTZ48meLiYnbv3k0gEGD16tUUFBQ0u6agoIBnnnkGgE2bNpGamkpm\nZmabbceNG8e+ffvYtWsXu3btIjs7m/fee69FwINoyNMWKiIiItLLdXslz+l0smzZMqZPn044HGb+\n/PmMHj2aJ598EoAFCxYwc+ZM1q1bR05ODklJSaxYsaLdtsdqb0g4rIUXIiIiEgcM60Qmv53hDMPg\nF3yNr32xFOOrX4WCAt5/HxITISenp3snIiIi0sgwjBNao3Csbq/k9TQLE8zGOXkXXwzBIPj9Pdwx\nERERkU4Uf8eaGQ6sYxZeBAI92CERERGRLhB3Ic8yTCzDAZEIYFfxRERERHqb+At5pollmrFKXvzM\nSBQREZF4En8hDxOryZw8t7uHOyQiIiLSBeIv5JnRkBcdrnW5erhDIiIiIl0g/kKeYdorbKOVvOM8\nZU1ERETkjBJ3IQ+z+XCtQp6IiIj0RnEX8iKGo9lwrUKeiIiI9EZxF/IswyRiNA7XRrNe7E8RERGR\n3iDuQt6xw7UN++RpQ2QRERHpTeIu5B27GXIgYK+wVcgTERGR3iTuQh6miRUdrrUsO9wlJ+vsWhER\nEeld4i7k2XPy7EpeKAROJ3i9pxDy9u7t1P6JiIiIdIa4C3lNK3mBgH3ihcdzksO1tbUweDCUl3d6\nN0VERERORdyFPMtwEIkuvGgIeUmuAMbf3z2h+xQXQ+ETu+wnpaVd0FMRERGRkxd/Ic9sHK5tCHmf\nC77BiDlTYituj8fChfD9b9XaTyoqYq/v3AmW1dm9FhERETkxxxXywuEw5eXl7NmzJ/Y4Y7UyXJtm\nHrF/Vl193LdxOKAPNfaTmprY62efDRs2dGaHRURERE6cs6MLHn/8cRYvXkxGRgYOhyP2+gcffNCl\nHesyDQsvmoS8/tZ++2eVlZCefny3MVqGvLo6+2lZWWd3WkREROTEdBjyfvrTn/LJJ5+Qfpzh57QX\nG64NxEJecjA67FpZedy3CQYh9ZiQt2+f/bS2rBoOBGHAgM7suYiIiMhx63C4dujQoaSkpHRHX7qF\nZZpYNB+uTTB99g9ra4/7PsFgy0re0aP20/wVX4EpUzqz2yIiIiInpMNK3ogRI7jsssu46qqrcLvd\nABiGwV133dXlnesKlumIDdf6/XbI89ZHN8lrGG89Dn5/k5AXDYcNU/P6H9wGRz7rzG6LiIiInJAO\nQ97QoUMZOnQogUCAQCCAZVkYhtEdfesapmlvoRKJEAzamyEnGNFKXn39cd8mELBDXiQxCTO6k3JD\nyAuY3s7utYiIiMgJ6TDkLVq0CICj0bHI5OTkU37TwsJC7rjjDsLhMP/+7//Ovffe2+Ka22+/nfXr\n15OYmMjTTz/NxIkT2237rW99i5dffhm3283ZZ5/NihUr6Nu3b8s3N0wilj1cGw7bIc/DiYe8hkpe\nODUdM7qTckPIsyLaQ0VERER6Vodz8j744AMmTpzI2LFjGTt2LJMmTeLDDz886TcMh8MsXLiQwsJC\ntm7dyqpVq9i2bVuza9atW8eOHTsoLi5m+fLl3HbbbR22nTZtGh999BH//Oc/GTlyJA899FDrHXCY\nsc2Qw2F7KxSP5SdimCdcyUs2agml9IudiVZXBy4XOMInc3yGiIiISOfpMOTdeuut/PjHP47tj/fo\no49y6623nvQbFhUVkZOTw/Dhw3G5XMydO5c1a9Y0u2bt2rXMmzcPgAsuuIDq6moqKirabZufn49p\nmrE2pW2dQmGYRLCHaxtDng9/YtoJz8lLdvsJJvaNhbxAAPr2bRLytCuyiIiI9JAOh2vr6uq47LLL\nYs/z8vKoPYFVqMcqKytjyJAhsefZ2dls3ry5w2vKysooLy/vsC3Ar3/9a7785S+3+v6vV65hb5Gf\n9PB+koZuxOHIwx0NeQknOFyb5PIT8Ka0CHneqqONF3k1P09EREQ6tnHjRjZu3Nhp9zuu1bX3338/\nN954I5Zl8fzzz3PWWWed9Bse76IN6ySrYA888ABut5uvfOUrrf788szZXHvhUUb7/sG6sXm88Qa4\nIz58ntQTHq5NcgUIeJLBb4e6YBBSUyIkhGvstFdbq5AnIiIixyUvL4+8vLzY88WLF5/S/Tocrv31\nr3/N/v37ue6665g9ezYHDhzg17/+9Um/YVZWFiUlJbHnJSUlZGdnt3tNaWkp2dnZHbZ9+umnWbdu\nHc8//3yb72+ZjtjCi1DIHq51W37qEvq1H/Kqq+Hw4dhTvx8SHAH8ruRmlbx+yQHChhOSk09o+FdE\nRESkM3VYyevXrx+PP/54p73h5MmTKS4uZvfu3QwePJjVq1ezatWqZtcUFBSwbNky5s6dy6ZNm0hN\nTSUzM5P09PQ22xYWFvLII4/w5ptv4m2nemaYjceaNczJc4V91PTJaj+UzZ4NRUWxHY8DAUjs46fe\nnQJHG0Neep8AAdy4ExMV8kRERKTHtBnyvvnNb/Kzn/2Mq6++usXPDMNg7dq1J/eGTifLli1j+vTp\nhMNh5s+fz+jRo3nyyScBWLBgATNnzmTdunXk5OSQlJTEihUr2m0L8B//8R8EAgHy8/MBuOiii/jF\nL37RsgOmSfiYhReusI86dwfDtR9/HNsjJRKxh2a9ZoB6xzGVvD4BAoYHFPJERESkB7UZ8m666SYA\n/vM//7PFz051M+QZM2YwY8aMZq8tWLCg2fNly5Ydd1uA4uLi43pvw2Ha26U0reRF/NS606B+V9sN\nnY1fVcNxaG4CHHI2D3lpSXYlj8TE2EkYO3fCO+/ADTccVxdFRERETlmbc/ImTZoEwD/+8Y/YRMCG\nx5YtW7qtg52uoZLXZDNkZ8hHjTOt/UqexxP7a0PIc0X8+FoJeX7L3ayS973vwdy5XfqpRERERJrp\ncOHFypUrW7z29NNPd0VfuoejcZ+8hoUXrrCPWmfjfnetikRif/r9duZzRgLUHTNc2zchGvKSkmKV\nPIfDbqpt80RERKS7tDlcu2rVKn7zm9+wa9euZvPyjh49Snp6erd0rkuYJhGaD9eakRC1Rh/w+dpu\n11Dlq68nEEiyK3nBAPXOZDvdYc/TS4ouvIh4vJjR+x05YjetqoJ+/bryw4mIiIjY2gx5F198MYMG\nDeLAgQPcfffdsX3rkpOTmTBhQrd1sNOZjmbDtbGQZ/aBunZCXm0tmCbU1eH3J9mVPJ+fGrP5Zshe\nM0DQcBNxeTGjrx86ZN/i4P4I/fp1WDwVEREROWVthrxhw4YxbNgwNm3a1J396XL2wovmq2vNSIg6\nktocrl22DL5RW4uRkQG1tfj9A9ocrvUYAUKmm7DTgzNayautBQ8+Ro5OgM8+g6FDu+3zioiISHzq\nsKz0t7/9jSlTptCnTx9cLhemaZKSktIdfesapknYajlce5TWh2stC+76jwCRsAWpqVBXF1t44QgH\nqDGOCXn4CZluQi5v7H4+H1zcd6t9w3fe6baPKiIiIvGrw5C3cOFCfvOb35Cbm4vP5+Opp57i61//\nenf0rUsYDpNwdDPkUMheXWuEQ9TQp9VK3oEDkEQtAVef2GKK2MKLsN8OedE5eU0reSGHJ3Y/nw9G\npBy0b7hnT7d9VhEREYlfxzVBLDc3l3A4jMPh4Oabb6awsLCr+9V1TNM+1uyY4dqaSFKrlbzqajvk\n1TuSYiGvoZJnhgLURBLtcl84TCAALitA2OEm5Gis5NXXw5CkSvuGBw9256cVERGRONXhsWZJSUn4\n/X4mTJjAPffcw8CBA2OLMM5Ix+yT53DYlbyjVuuVvMOHoyHPTIrtfef32JU8RziAL+K2E18wSCDg\nwE2AGoeboMMDfnsLFZ8PBmdW2Tesru7OTysiIiJxqsNK3jPPPEMkEmHZsmUkJiZSWlrKiy++2B19\n6xKGo/mxZk7T3v+uJpzQbiWvzmheyUtwhzGsCL6gww55gYBd4SNA2OEhaDav5A10V3I4dai9j4qI\niIhIF+uwkjd8+HAAEhISWLRoURd3pxs4HM0WXriMEJbDSW3Y22Ylb0T/Gmp9jSHP74JEZ4CI000w\nZMRCXjBoD9dGnG4Cpj0nz7LsrNffrORQ37Poq5AnIiIi3aDNkHfuuee22cgwDN5///0u6VBXa7rw\nIhwGt2mHPF/Iac+ta1iNEVVdDUP711KzJwm8dnXO74dEVzTMBQCXK1bJc1kBwk43gWglLxCwf5wS\nrmJ/0lmcVb215z68iIiIxI02Q94f//jH7uxHtzEcJmGr8Vgzr2GHukCAWIijT5/Y9dXVkJ1Wy5Hd\nSfZEvGiYS3AE7UpekGbDtS4rgOV0E8Cu5NXXQ0ICJIaPsNM7Dqre7rHPLiIiIvGjzZDXMEwLsHv3\nbnbs2MEVV1xBXV0d4XC4O/rWJew5eY3DtU7sSl4wiB3i/P5mIe/wYRiTUsuRUJId5vz+WCXPcroa\nQ14wSCAATitAxOXGb9iB0eezs6MnXM9+z2DNyRMREZFu0eHCi+XLl3P99dezYMECAEpLS5k1a1aX\nd6zLmCaRJqtrXa1V8pqoroaMPrUcDiVhue0QGAhAoiOA5YoO1zap5Dkj0dcN+1qfz67kucP17CcT\njh495Y/w9NNQ/GY5vPbaKd9LREREeqcOQ97Pf/5z/vrXv8ZOuRg5ciT79+/v8o51FcNhErIaV9c2\nhLxgEDvkHbP44vBh6Ouopd7Zh4Bhhzm/H7yOILhcLYZrneEAuNz4sANjfb19W1e4nv2RdPv+odBJ\n99+y4OabYc8tiyA/HyKRU/o+REREpHfqMOR5PB48Hk/seSgUwjCMLu1UV7Ln5JktKnmx4dpWKnl9\nzFoCriRC0RWzfr+9ujZWyWuy8MIZ9mO53fisxkqe1wvOYD2HA4n2UPCJVvNmzYL33gOgvNx+KbGy\n1P7LZ58B9i1/97tT+GJERESkV+kw5F166aU88MAD1NXVsWHDBq6//nquvvrq7uhb13A4CNHkWDMj\nhNF0uLaVSl6iVUvQlUTQbFx44TGD4G6+8CIYBEfEPg6jvkklLyHBDnnV/gRIToajRwmF4LrrYOfO\ntrtaUwMbV34Ga9bAb38LwMcfw7BhMLCmGPr3jx2T9qtfwZw5cGDlOs37ExERkY5D3tKlSxkwYADn\nnnsuTz75JDNnzuSHP/xhd/StS5jOxtW14TC4CIGr/UpeYqSWgDuJoNG48CLBYe+NEpuTF1144QgH\nMNxufJHmlTxHoJ4qX2PI27IF/vAHWLWq7b7+5jfwP//vb/aTHTsAu3B38QVhBof2YF1yCZTaFb2P\nP4Z0DjLg/10Fjz/eBd+ciIiInEna3Qw5FAoxbtw4Pv74Y2699dbu6lPXMk3CFs1W13ZUyfOGagi5\nkwgYodjCC48jWsmro9mcPEc4AB43dVbz1bWmv56qSAJk2CHv0932/T/4oO2uvvoq5LCbgxMup3+0\n5HfgAJwzoJIaIxnvwLNJKikBYO9emJr5EewD3n23s781EREROcO0W8lzOp2cc845fBad99UbmM5j\njjU7jkqeO1RLyJtEgCYLL4wAhruxkmf57eFaMxTA8HioDzffJ8/011NZnwApKXD0KLt2wdSpsG1b\n23395BOYMqiU4kGXwqefAnbIy/Ye5Ii7P5WJWbFKXkUFXHZOGTsHnB+bvyciIiLxq8Ph2srKSsaO\nHcvll1/O1VdfzdVXX01BQUF39K1LHLvwwkkIIxryrGO2UIlE7AUNLn8tYU8SfhoXXnjMIEbDnDyX\ni1B9AJcLjEAAw+umLtKkkuexoL6eg7WNw7W7d8OVV8J/ffgVIg8+1KKflmXP1xvbt5QPrbH2itza\nWvbvh8GuA9QkDOCQM9NOfdiVvPH9y/lnwkX2a9GK5L//O6xY0R3frIiIiJxOOjy79oc//CGWZTV7\n7UxeXYtpEooO14ZC9upaw+nE6QTL5cHw+9m5085ibrd9XK1RV0skIQk/gcbhWm8AmlTywnV2yCMQ\nwPS4qQ03zsnr4wmCYRDERTgpGceRI+zfD9MvC3BdZBXhR/4E374Pmnyv+/fbFcDMYCmf1A2BAQPg\nwAEOHEgiY8gB/H36U0kG7N9PJAL79sHZnjL+6hsCgwdDaSl7XGfz1FNQ8er73HzwT3D33c3eQ0RE\nRHqvdit5oVCIW2+9lby8vGaPSy+99JTetLCwkFGjRpGbm8vSpUtbveb2228nNzeXCRMmsGXLlg7b\nVlZWkp+fz8iRI5k2bRrV1dWt3td0RVfXNh2udTpxuSDitqtvX/iCXWWrrobUVKDWDnk+q3G41mME\nMDyNq2tDviBuN3bI87qpDTWurk121mMkJJCUBCGvXcmrqoLswx+xO3E0AU+yPTbbxM6dcNZZkFRV\nyraapiEP0iMHCaYOYG/IDnmVlfbOLGm+crYdHow1dCjs2cPGjTB7NtxX8U24554Tmqv30Ufw6eu7\n4cEHsT+kHTy/8x048tLr7Y8zi4iISI/rcE7eqFGjOnVOXjgcZuHChRQWFrJ161ZWrVrFtmMCw7p1\n69ixYwcucmElAAAgAElEQVTFxcUsX76c2267rcO2S5YsIT8/n+3btzN16lSWLFnSegccJuFI43Ct\no0nIC7vs6lt5ub1atVnIS+oT2/suEACXEcT0uptV8pqHvGglr94i2WlPzEtKwg500ZCXUb2dgxlj\n2HNWHmzcaCe7z30OPvqIXbsgZ2gA19FKPq7MaBby+gYOEEkfwB6fHfL27oWBA8G1v4xyM4vgQDvk\nffwxnDcuxJTIZrbOvLv9pbyWZT+Agwfh4oth95Vfw/re92DlSgAWL4Y//ugT+lyXD1ddFdvU+Ze/\nhP95tNZeKXJM1VdERER6RrfPySsqKiInJ4fhw4fjcrmYO3cua9asaXbN2rVrmTdvHgAXXHAB1dXV\nVFRUtNu2aZt58+bx0ksvtf6BG068aJiTZ9khz+2GiMuLVW/Pyaursxcz9O8P1NZCYpK9LUpDJY8A\nprvxxIuIr0nIS3DjCzrA4cBfE4yFvMRE8HtSYiEvpa4CBg3iH33z7JD36KPwf/8HK1awcydMyNiL\nlTmQigMOGDAAa/8B9u+HPvUHMDP781lNOlRVUVEWZtAgMMrLMbIGU51ih7wdO2Bi0nbq0rJYk/hl\neOUVwF7DsXmz/X2EQvDKUxWEhgyHCy+EQICVK+GWaaVM4R3+dvtv4fnnCYdh9WpYc+PvecrzDSIZ\nmbB+PZ9+CvfeC+b3/xtrxgz4wQ8Aey7jK69AzZHjOJGjSTAMh+Hvf4farZ/FFpuAvWfgvn3YyVtB\nUkREpEMdzsm7//77gcZ5eJZlndKcvLKyMoYMGRJ7np2dzeaGxNHONWVlZZSXl7fZdt++fWRmZgKQ\nmZnJvn37Wn3/V/7yY44cNeDoUSoqNuI4p2klz4v/iJ8+feBzKR+wZdNI+vf3wJZa6JOEL9K4T57b\nCGJEK3mWy0W4vjHkORoqfB4PoVo/SWa0kucBnysZjpZTWQlJRyvwDhvInyvzmPO7+fjMREqXvEju\nqh+waxIUDCjDzBoMlRBIHYBVfgDDANfhg7gHT2H/pw5IS6NqxyEGDRwAb5fj+cJg9nmHkrHnXXbs\ngFET3idy7gRWf3Ie3z5yhENFO7jwqhxu8D/DvpydPGP9G9/efjO/cP8/vpH0No7nnuO5525h9Rde\nYu+/fJFf753Bxe/ezLt/PkxWVl9GvL+G90c8xNbzxzLu2Wf51dlXs/CGA8x7ZiXfnPI3HntsBoG5\nN5E/bwRTy1aSV/YN6rMGwV13sdUxHt/BGgb0C/OB619Yv6qayz5axg21T+G4Mp8jP36Ka28dwPgP\nnuf7lbfjSTFwLvour5/7TW69voqfH/k3Lo+8hjllEsGfPcFPXp/A9hVvc4vv55w9uJ4+101na8qF\n7Pw4QEJlGQPCFXicYQzTwDIMXB4H9EujknSq9/vxVHzGEKOMPmdl4Bt2Dp8dToV3isg4vIPkoWmY\n545lq/9sjnyyl4wjxQz0HsY7fCDViYM56E+m9nAIRySIkxBOy/7TQQiHYWE2fZgWhgER00kw4iAU\ntHAEfLiNAA6Pi4jLjT/iIuIL4oz4cZkRTI+LsOkiZJlYwTCmFcIgOqXSNIDW/hs8kfDbSvt2/ru2\nNJdTRKTLbSnbxpbyzpsO1WHIy8vLo6KignfeeQfDMDj//PPJyMg46Tc83oB47GKPtq5p7X6GYbT5\nPjOnf5c9nwZZtOvX/KNfHi7WxSp5IaeHQKWPgelBCj8bz89/9wsGXHob1NZi9kmiLtJkuNayK3mm\nCZbTTdgfnZPn9+NIdNuLW71ewrU+kjzRkJcA9c5kItVH8PnAXVVB6qjP8+ZvhvPeRd/ghb8NYdfb\nM3nhk3+lpK+fEYPLICuLjH1w1DsAZ8lBBgwADh4kYfIA9v8NyMigZud+zkp1QEICGcMTKWEo4/b8\nLzt2wOCD/8T1hfHseNQkcO1M3l38CosmDOVr279HkXEVv/3wXBxfmctDR75PZkoh1/zov9lXdTO5\nH73Eoa8u5KV7k1j++S+w84k/MfeSz8PzOxj535fwi3f+hZ//6R7+11vN29f9DOe/zmHNn87nruvu\npG7GQv4teRbf4Ae89pMi/uf+A8z91hOMSnwW09uHQz6Y5nuXmUkeDl4znzkffMy17/6Cy0ZP4Udn\nXcjElHdY882N3P/TFF5edCW19X/hg37vYcy/mluq/sDg15/lWxfnc7O3D31TDf6Rdwe/ei+dkd8t\n5HznMrKSPBxJyeagcyBVESdYYFgRIqEw3voqMqxDDEpwc7TfMLaHB2Ns/pih9WsY7ajkwPAplKaM\nonZzFQNf+Q3nGZ8SSB/E/tSRFPtScb/9LpmRvZxlHrWDmOEkbLgIGU5ChosIDiIYRCwDyzIIYxC0\n7H+LphXGTQivwyDgTMBnubBCIZyRAF4CWC4XQdNDfcTECAXt4GhEsEwnEcNBGAMsCwMLw7JaDV5W\nq+HvmP8+WgmDRnv/vVlW65nyFLTVfxGReHYecB6JgP2/k78+xft1GPJeeOEFvvWtb8UWWyxcuJBH\nHnmE66+//qTeMCsri5LoBr4AJSUlZGdnt3tNaWkp2dnZBIPBFq9nZWUBdvWuoqKCgQMHsnfv3jaD\naLhPX0JWVWx1rcNqrOSFnF581T7+JbkYALP4EzJnBSASwZnopj7UuPDCZdibIbvdEHa6idQ0VvKc\nCdGQ11DJS2gcrq11JBOsPEpqKhgVFQyaPYiyR2Dy9sd47TX7mNrI2Tk4i7cy+At2yMssgyrnAPqW\n77CHjw8coM/w/vbuKUMzqN+zn7OHWzB4MEOGwI6yoYR3fobLBQnb34evfpVJk+DvI+dy7qKbuTw1\niOPll7jooovAWgaGwUPFcPGFV3KJ/5s8NPUPGG8U0f+P0xj2C/h01BdJevoV5vx7JcyYwey5Lsbe\nn8Z951zBHQd/Sf/f/Q8UFfGjfJh889381rqBr3l+hlFYSP7YMVxxO1jWpZhNJwdEQ8UQw+B3Yfjf\n/13KgZLp/IvzQ7jpf5iVmsrAK+DZP27iq66nSfj81+CKK3jGgs2b57MvPJuxaeUwahQXmSYXAfCV\nE/73eP4xz4e3cV3OCd9ZRETkFBm/OaXmx7WFyjvvvBMLTQcOHGDq1KknHfImT55McXExu3fvZvDg\nwaxevZpVxywIKCgoYNmyZcydO5dNmzaRmppKZmYm6enpbbYtKChg5cqV3HvvvaxcuZJZs2a1+v6m\niT0nL7q6tmHhhdsNIYcH/+FqchLKABgY+AxHci0kJeFNMKira1LJi9jHmrlcEDabz8lzNqnkRep8\nJNC48KLWTCZUfZS0NKCiAsfggdx7r31c2dSpMGEClLsmMPjjf5LmK4fBg8nMhINWf5L3HbBD3raD\n9M2JhrzJGYRK9jPkrABkZTFkCPzfJyMw9nxG7oQIvP8+TJjApZfCjB9dwaNn/wfzHz4HLrJjUcMQ\nXW4ufOe7Jk/97D7+e+1sWLgQEhOZNQv+88UvsuLw9+j3chE8/DCDB8N558HX35rPy5Gr4Otfh7PO\n4vqz4JxzvAwdugZnauN3bhitjAQ2ecHhAPuf0+XRh+3CC+HCC/sC32zW7MILAVKjDxEREWlNhyHP\nsiwGDBgQe56enn5cQ6ltvqHTybJly5g+fTrhcJj58+czevRonnzySQAWLFjAzJkzWbduHTk5OSQl\nJbEiuptvW20B7rvvPubMmcNTTz3F8OHDeeGFF1p9f4eDZpshN6/kJRA+Ws9gxz5qB+cwtHwP3tSj\nkJKC1wu1Rxs3Q3bRWMkLmS4sX529T15NAGeSp3G4tt5PotEY8o4aKUQOH6VfP+Aze1nsd7/b2L/8\nfPjjyglcmvZPHHv3w4RzycyEveEBnH3oAANygL8cJOWsAfh8EErPwPj7fgZZdY2VvPJE/Amp5PX/\nED4+DMOGceedsHOnwcXfuRdGt/67ueMO4Ju3wFs5cL5d4/r61+GiZ4ey97J/JT3wHsyYAcAvfgFv\nvz0T64IPMEadE7vH+PEn+y9DREREOlOHIe/KK69k+vTpfOUrX8GyLFavXs2M6P/Rn6wZM2a0uMeC\nBQuaPV+2bNlxtwXo168fr732WofvbZoQjDRZXRsJxkJewNMHamoYYO2DyVMYunYDGcOOxELe0WDj\nGbVOy67kud0QdtjHmjVU8lxJjcO1kTofXhqHa49a9hYq6dlheO8QHDOsPG0afH/RBB7PWg+f1sHw\n4WRuh9LDGbiqDjCobx0ARlIiAwZAbWIGzqr99PcfjlXySkrgQMrZzDy62i4NmiZpafDccx1+PXap\nrMk+iOnpsH07wM+aXTZ6tP2AccdxUxEREelubW6hUlxczF//+lceeeQRFixYwPvvv88HH3zAxRdf\nzK233tqdfexUDkfbw7UBdzJG7VH6h/aRdME4+jurGZJwqLGSF2qs5DkjAXC77Qqg6cYKBFuGPK8X\nq96P12qs5B2xknHUHGFowgHo1w+czXP2hRfCzPsmkFPzD9i6FcaOZeBA2FmbScLRfXa7/v0Be+u8\nw54MEo7uJ6WmPFbJKy2Ff7qmcNG7jzWMbYqIiEicaTPk3XHHHaSkpAAwe/ZsfvzjH/PjH/+YWbNm\nceedd3ZbBzubaTYO14ZC4Ig0Dtf63ck4ao+SFtwHgwZh9O8PxcXQty8JCVATbNwnzxlpOlzrxgq0\nXsnD52sW8qrDyTjqjzLUXWHvYHwMw4A7HsrEcbjK3myuXz+ysuDT/ckQDjPM2h0LeRkZUB7KoF/4\nAO5D5ZCVRWKifRTbqvJLcflq4PLLW7yHiIiI9H5tDtfu27eP8a1MsBo/fjy7du3q0k51pdhwbUMl\nr8lmyD5XH1Lrj9LXuQ8yM+0Qtn17YyUv4LKHa/0WjkjjcG3IiA7XpgKBAO4+TSp5Pj+ecONwbVVN\nMi5/DYPNva2GvJjvfhd7kp99FG35XoMqz0CGHP3ILuFh//HR/gFMcO3H2FMH0VXKQ4fC6n8U8Nwr\nhTiunNbF36iIiIicjtoMeW2d/Qrg8/m6pDPdweGASHTfskgogmmFYqtk/a5kXL4aEh3+xpBXXAzp\n6Xi9UO83weXC8gdwNKnkBbHDX0MlLxbyPB7w+3BH6iHJruQdPOgk6ExgaP329kNe9OQIgKwsKC+H\nCgYy5MCW2Dy+AQNg05YMrjT3wY59kGNv9PHEE1BSYuKYOb2rvkYRERE5zbU5XDt58mSWL1/e4vVf\n/vKXTJo0qUs71ZVMe6Q2usw2HBuutSt5ybgDR0k6Gh1KPaaSV18P9uS9AM5w4xYqQcPdLOQ1nZOH\nz4873GQLlVqo9g5k6P53Y5W3jmRm2qd57agbTNqHb8FZZwF2yCv8aAhZ9Z9CJBI9aNeehneSO9yI\niIhIL9FmJe+nP/0p1157Lc8//3ws1P3973/H7/fzhz/8ods62Nkc9kgtOBxEQpFmW6jUOZJJ8Ffj\n9dXZySozEz75BK6/Hq8X/H6wPB6ctX6McJNKnunGFQzidkYgFMKT5IxV8syAD1eoHhL62psh18JB\n1yCG7imCIZ87rj47nXZXtu4ew5d2vAhnnw3AiBFQeiiBiGFiOhxd96WJiIjIGafNkDdw4ED+7//+\njzfeeIMPP/wQwzD44he/yOVn+ET+WCXPNLFC4ehwrR3y6p3JJIWO4O87AI/TaVfywmEYOLBJJc9D\nH3cAI9A4Jy+IGyMYIMFhr7j1eI1YJc8I+qMhL4E+9g4t7DMHMbb8r/bkueN0+eXge3McfEpsWDa6\nlR3//M7vmDi1X2d/VSIiInIGa3efPMMwuPzyy8/4YNeUaTZW8qzQMcO1ETcAVlKyfXF0WJSsLBIS\nwOcDy+2mj8sPwWBsC5XAMSHP6bRP7Yq47UqeM2CHvPR0qKyEslB0Lt6QIcfd71/+EsIH8+HBb8bS\nXU4OvPsunPcv13X62aIiIiJyZmtzTl5v1XS41gpHMJtsoRIM2teYROy/TJ5s/zlpUqySF3F57JDX\npJIXsFzNQp5h2Gsuwk4vjqAfR9AOef37w8GD8F4guoHwqFHH3W/TBFdGGvz0p8321ps0qZUjw0RE\nRCTuxV3IazpcSzgcC3nR9RQcMDMITIpuIDx4MNTZx4V5vXYlL+LykOQO2BdH5+QFcGM2hDyPB7D/\nCDo8uCM+TH9jyNu7F5b7biKy7ZPYFikiIiIina3DY816m6aVvMaQlxir5J2X+hn/+IVBn4YGCQmx\nP+rrIZLkJsnZOFzrdoPfcmOEg3ijlTyI7oNseUly1mLU2yEvLQ2OHIH0dC/mqJE98vlFREQkPsRv\nJa9huDa6utZjn1jGoVovSf08Ldo1reT1cTYO17pcdsgzQwG8ZvOQVxfxkOTw2+kwISE2yqoCnoiI\niHS1uAt5sUpew+ra6HBtUpK9F10oFCveNdNQyQs73CQ4Ay0qeWYogNc4JuSFvSQ5fPaQb1JS7F7O\nuKufioiISHeLu7jRdHUtgTBmuDHkffoppKS0vpDBYx9bS9gRrc4ds/DCEQrgOSbk1YQ8JDj89uZ4\n0ZD3xBP2nnciIiIiXSkuQ15sn7xwBKOhkue2F0UkJ7fermHFbMD0knjMnDxfxI0Zbi3keUk0fXbI\nS0wE4Gtf66YPKiIiInEt7kJes4UXkSaVvASoqLAreW3xesFnRINbkzl5vogbRyTYMuQFPCQY/hbD\ntSIiIiJdLe7m5DU/uzbSbE5ee5U8iM7Li3jo4/A1q+TVh904jqnkud1wOOAlwfA1G64VERER6Q5x\nF/KOXXhhNAl5Bw50XMmri0QXUzSZk+e33DgjATz4m1Xyjvg8JFJnX+v1ds8HFBERESEOQ17TSp4R\nCWM0WXgBHVfyakPROXlNNkOuD7lwRgK4rObDtdU+LynhSns+no6lEBERkW4UdyHv2GPNjg15HVXy\naoIee05eMNg4Jy/kxGGF8RrNK3lVdR76Bg9pqFZERES6XdyFvIYtVCzTxLBaVvJSU9tu6/XCkYCX\nRMNn74zs9eJ2QzBkEDDceINHY8OyHg8crPHSJ1ipkCciIiLdLi5DXsNwrZMwRsgOeQ3DtAMHtt02\nIQEO+6OLKaKnWDSceRs03Hj8R5qdXbv/sAdvSIsuREREpPvFXchrGK61TAcuR8Q+4sLpJCPD/nnf\nvm239XrteXZJ4SP2sRUORyzkhQyXHfKaVPL2HY4utujTp+2bioiIiHSBbg15lZWV5OfnM3LkSKZN\nm0Z1dXWr1xUWFjJq1Chyc3NZunRph+03bNjA5MmTGT9+PJMnT+aNN95osw8NlTzLMHE7wrGQZ5rw\nne/Adde13f+EBKiu95AUrIqdfeZyRSt5uHEfE/LKqqMVvH79TuRrEhERETll3RrylixZQn5+Ptu3\nb2fq1KksWbKkxTXhcJiFCxdSWFjI1q1bWbVqFdu2bWu3/YABA3j55Zd5//33WblyJTfeeGObfWjc\nQsWBy2wMeQA//CH07992/71eqKzzkhiojoU8t9tegxEw3Lh8zefkVQSj4a69iX4iIiIiXaBbQ97a\ntWuZN28eAPPmzeOll15qcU1RURE5OTkMHz4cl8vF3LlzWbNmTbvtzzvvPAZGJ9ONGTOG+vp6gsFg\nq32IVfIcDpxmpFnI60hSEhwJeknwVTULeYEABCw3zvrmIS/ccKCIw3Fc9xcRERHpLN16rNm+ffvI\nzMwEIDMzk3379rW4pqysjCFDhsSeZ2dns3nz5uNu/+KLLzJp0iRcLlerfVi2bBGHDsEPfCUcsf5+\nQiEvJQUO4sFbXwV9mw/X+nHjrGs+XBvT5POIiIiItGbjxo1s3Lix0+7X6SEvPz+fioqKFq8/8MAD\nzZ4bhoHRygbBx75mWVab1x37+kcffcR9993Hhg0b2uzfnXcu4o9/hP/OepMPjo6H0HMnFPJ8eHHX\nVcPALKBpJc+Fo7Yx5DUccPHaqgNccV07m++JiIiIAHl5eeTl5cWeL168+JTu1+khr72AlZmZSUVF\nBQMHDmTv3r1kNCxpbSIrK4uSkpLY89LSUrKysjpsX1paynXXXcezzz7LiBEj2uxDbLgWs8WcvI70\n7Qt1JOI5ehAScgC7Yuf320ebmTWNIS893W7jyeoP7uO6vYiIiEin6dY5eQUFBaxcuRKAlStXMmvW\nrBbXTJ48meLiYnbv3k0gEGD16tUUFBS02766upqrrrqKpUuXctFFF7Xbh2ZbqJxgyEtPh2pSMcLh\n2Jy8xESorbVDnlHTOCdvwAC7TXv77omIiIh0lW4NeQ1DqSNHjuT111/nvvvuA6C8vJyrrroKAKfT\nybJly5g+fTpjxozhhhtuYPTo0e22X7ZsGZ9++imLFy9m4sSJTJw4kYMHD7bah4ZKXsQ88YUXgwfb\nIQ+IhbykJKiqsjdDNo40VvKiXdZ0PBEREekRhmVZVk93orsYhkFFhcX48bBzzFXcvu02nkq+E9at\ng9zcDtvX1sKLv6zkpjvTYc4cWL2a3bvhvPPg9zXTuYI/w5o1EA2s0UMxRERERE6YYRicSkyLuxMv\nGs+udeA+weHapCS4aWF0EUWTSt7hwxAxXXaJsGHFBQp4IiIi0nPiMuSFwxAxTny4Fmi81m2vpmg4\nljbgjCa6JiFPREREpKfEXciLLbwwTJxG2D6u4kRCXoPoPnwN1bp6Z7L9F4U8EREROQ3EXchrWsk7\n0dW1MVOmwBVXANCwVV+dIxrydISZiIiInAa69cSL00FjJc+BwziJ4VqAoqIWL9WZ0ZDXt28n9FJE\nRETk1MRxJc/EebKVvFbUOqPhTiFPRERETgNxF/IaKnnNhmvbOOf2hJjRr7Iz7iUiIiJyiuJuuLZh\nC5WI4cB5ssO1rfhz6nXcfUtlJ/RQRERE5NTFZchrOLvWaYTtJw7HKd3zJz+B3NwRcNUDndRLERER\nkVMTlyEPIGw48OC3A17DEtmTdMcdndAxERERkU4Ud3PyIDpka5l2yNMcOhEREemF4jLkORwQwoEX\nn0KeiIiI9EpxGfJME8I48FgKeSIiItI7xW3IC1kKeSIiItJ7xWXIczggZLjwWnUKeSIiItIrxWXI\nM00I4sIbUcgTERGR3ikuQ5698EIhT0RERHqvuAx5pgkBhTwRERHpxeIy5DkcELRceMIKeSIiItI7\nxWXIa6jkKeSJiIhIbxW3Ic9vufGGaxXyREREpFeKy5DncEAg4sITUsgTERGR3ikuQ55pgj/iwh3S\ncK2IiIj0Tt0a8iorK8nPz2fkyJFMmzaN6urqVq8rLCxk1KhR5ObmsnTp0uNuv2fPHvr06cOjjz7a\nbj8cDoU8ERER6d26NeQtWbKE/Px8tm/fztSpU1myZEmLa8LhMAsXLqSwsJCtW7eyatUqtm3bdlzt\n77rrLq666qoO+9FQyXMFFfJERESkd+rWkLd27VrmzZsHwLx583jppZdaXFNUVEROTg7Dhw/H5XIx\nd+5c1qxZ02H7l156ibPOOosxY8Z02A+HA3wRF85IUCFPREREeiVnd77Zvn37yMzMBCAzM5N9+/a1\nuKasrIwhQ4bEnmdnZ7N58+Z229fU1PDwww/z2muv8cgjj7Tbh0WLFnHoELxSt52JQJ5CnoiIiJwG\nNm7cyMaNGzvtfp0e8vLz86moqGjx+gMPPNDsuWEYGIbR4rpjX7Msq83rGl5ftGgRd955J4mJiViW\n1W7/Fi1axO9/D5/r8yfySlaB293hZxIRERHpanl5eeTl5cWeL168+JTu1+khb8OGDW3+LDMzk4qK\nCgYOHMjevXvJyMhocU1WVhYlJSWx56WlpWRlZbXbvqioiBdffJF77rmH6upqTNMkISGBr3/96632\nw+EAXzhawVMlT0RERHqhbp2TV1BQwMqVKwFYuXIls2bNanHN5MmTKS4uZvfu3QQCAVavXk1BQUG7\n7d966y127drFrl27uOOOO/jOd77TZsADe+FFfUghT0RERHqvbg159913Hxs2bGDkyJG8/vrr3Hff\nfQCUl5fHVsU6nU6WLVvG9OnTGTNmDDfccAOjR49ut/2JUiVPREREejvD6mgSWy9iGAaWZXH++TDy\n6N957uPJcNdd0MG+eiIiIiLdrSG3nKy4PfEiNlybmNiznRERERHpAnEZ8hwOqA9HV9UmJPRsZ0RE\nRES6QFyGPNOEumC0kqeQJyIiIr1QXIY8hwNqg6rkiYiISO8VlyHPNOFgsK/9RHPyREREpBeKy5Dn\ncEBVsI/9xOvt2c6IiIiIdIG4DHmmCf5g9KP369eznRERERHpAvEb8vzw93t+C5de2tPdEREREel0\nnX527ZnA4QDLggOX3QA68EJERER6obit5AE44zLiioiISDyIy5DncNh/KuSJiIhIbxWXIa8h3Cnk\niYiISG8VlyHPHd0HuaGiJyIiItLbxHXIa/hTREREpLdRyBMRERHpheI65Hk8PdsPERERka6ikCci\nIiLSCynkiYiIiPRCCnkiIiIivZBCnoiIiEgvpJAnIiIi0gvFdcjTZsgiIiLSW8VlyAuHe7oHcqI2\nbtzY012QU6Df35lLv7szm35/8a1bQ15lZSX5+fmMHDmSadOmUV1d3ep1hYWFjBo1itzcXJYuXXpc\n7d9//30uuugixo0bx/jx4/H7/W32Q5sgn3n0P1RnNv3+zlz63Z3Z9PuLb90a8pYsWUJ+fj7bt29n\n6tSpLFmypMU14XCYhQsXUlhYyNatW1m1ahXbtm1rt30oFOLGG29k+fLlfPjhh7z55pu4XK42+/HV\nr8I773TNZxQRERE5HXRryFu7di3z5s0DYN68ebz00kstrikqKiInJ4fhw4fjcrmYO3cua9asabf9\nq6++yvjx4zn33HMBSEtLwzTb/mgeD0ye3KkfTUREROS0YliWZXXXm6WlpVFVVQWAZVn069cv9rzB\n73//e/70pz/xy1/+EoDnnnuOzZs38/jjj7fZ/qc//Snvvfce+/fv58CBA8ydO5dvfetbLd7fMIwu\n/oQiIiIinedUYpqzE/sBQH5+PhUVFS1ef+CBB5o9Nwyj1dB17GuWZbV5XcProVCIv/71r7z77rsk\nJNGfcwIAAAchSURBVCQwdepUJk2axOWXX97iXiIiIiLxoNND3oYNG9r8WWZmJhUVFQwcOJC9e/eS\nkZHR4pqsrCxKSkpiz0tLS8nKymq3/ZAhQ7jkkkvo168fADNnzuS9995rEfJERERE4kW3zskrKChg\n5cqVAKxcuZJZs2a1uGby5MkUFxeze/duAoEAq1evpqCgoN3206ZN44MPPqC+vp5QKMSbb77J2LFj\nu+lTiYiIiJx+unVOXmVlJXPmzGHPnj0MHz6cF154gdTUVMrLy/nqV7/KK6+8AsD69eu54447CIfD\nzJ8/n29/+9vttgd4/vnneeihhzAMg6uuuqrVlbsiIiIiccOKE+vXr7fOOeccKycnx1qyZElPd0fa\nsWfPHisvL88aM2aMNXbsWOtnP/uZZVmWdejQIeuKK66wcnNzrfz8fKuqqqqHeyrtCYVC1nnnnWd9\n8YtftCxLv78zRVVVlTV79mxr1KhR1ujRo61Nmzbpd3cGefDBB60xY8ZY48aNs7785S9bPp9Pv7/T\n1M0332z9//buJ6TJP44D+HvmgohBFJnKYySGj8w/7Qkz8BCRyTJaUXYoyGREFyn6Cx261KFJeCgP\nnYIudVhQHaQ/q8SEYNvBVhhOMOIZDC2JhQcjcNNPhx+NVr/NJdXzPNv7dZnPdw97PvCG8fH7+Hws\nKyuThoaG9FqurHw+n2zcuFFUVZWnT5/mdY2i+I8XuWbvkfnY7XZcu3YNY2NjCIfDuHHjBsbHx/Oa\ns0jm0d/fD6fTmX5AivlZw6lTp7B7926Mj49jdHQUdXV1zM4iYrEYbt68iUgkgrdv32J+fh5+v5/5\nmZTX60UgEMhYy5ZVNBrF3bt3EY1GEQgE0NPTg4WFhcUv8sdbUxMKBoPidrvTx729vdLb22tgRfQ7\n9u3bJ8+fPxdVVeXjx48iIvLhwwdRVdXgyiibeDwubW1tMjQ0lN7JY37mNzMzI9XV1b+sMztrSCQS\nUltbK58/f5ZkMil79uyRZ8+eMT8T03U9YycvW1Y+ny/jLqTb7ZZQKLTo5xfFTt7k5CSqqqrSx4qi\nYHJy0sCKKF+xWAyvX7/G1q1bMT09jXXr1gH470nr6elpg6ujbM6cOYO+vr6MoeTMz/x0XcfatWvh\n9XqxefNmHD9+HF++fGF2FrF69WqcO3cO69evR2VlJVatWoX29nbmZyHZspqamoKiKOnz8u1jiqLJ\n4xBka5qdnUVnZyf6+/vhcDgy3ss2Z5GM9/DhQ5SVlUHTtKyzKZmfOaVSKUQiEfT09CASiWDlypW/\n3Npjdub1/v17XL9+HbFYDFNTU5idncWdO3cyzmF+1rFYVvnkWBRN3s+z9+LxeEZHTOaTTCbR2dmJ\nrq6u9Kic73MSAWSds0jGCwaDGBgYQHV1NQ4fPoyhoSF0dXUxPwtQFAWKomDLli0AgIMHDyISiaC8\nvJzZWcDIyAhaW1uxZs0alJaW4sCBAwiFQszPQrJ9T+aaIZxLUTR5uWbvkfmICI4dOwan04nTp0+n\n1/OZs0jG8/l8iMfj0HUdfr8fO3bswO3bt5mfBZSXl6OqqgoTExMAgMHBQdTX18Pj8TA7C6irq0M4\nHMbXr18hIhgcHITT6WR+FpLte3Lv3r3w+/2Ym5uDrut49+4dWlpaFv/AP/kHhGb2+PFjqa2tlZqa\nGvH5fEaXQzm8fPlSbDabbNq0SVwul7hcLnny5IkkEglpa2vjGAALGR4eFo/HIyLC/CzizZs30tzc\nLE1NTbJ//36ZmZlhdhZy9erV9AiVo0ePytzcHPMzqUOHDklFRYXY7XZRFEVu3bqVM6srV65ITU2N\nqKoqgUAgr2v802HIRERERPRvFMXtWiIiIqJiwyaPiIiIqACxySMiIiIqQGzyiIiIiAoQmzwiop8k\nEglomgZN01BRUQFFUaBpGhwOB06cOGF0eUREeeHTtUREOVy+fBkOhwNnz541uhQiot/CnTwiokV8\n/114eHgYHo8HAHDp0iV0d3dj27Zt2LBhAx48eIDz58+jqakJHR0dSKVSAIBXr15h+/btaG5uxq5d\nu9LT7ImI/jY2eURES6TrOl68eIGBgQEcOXIE7e3tGB0dxYoVK/Do0SMkk0mcPHkS9+/fx8jICLxe\nLy5evGh02URUJEqNLoCIyIpsNhs6OjqwbNkyNDQ0YGFhAW63GwDQ2NiIWCyGiYkJjI2NYefOnQCA\n+fl5VFZWGlk2ERURNnlEREu0fPlyAEBJSQnsdnt6vaSkBKlUCiKC+vp6BINBo0okoiLG27VEREuQ\nzzNrqqri06dPCIfDAIBkMoloNPq3SyMiAsAmj4hoUTabLf36fz//eM6Px3a7Hffu3cOFCxfgcrmg\naRpCodC/K5yIihpHqBAREREVIO7kERERERUgNnlEREREBYhNHhEREVEBYpNHREREVIDY5BEREREV\nIDZ5RERERAXoG1aJN+VzkeOwAAAAAElFTkSuQmCC\n" } ], "prompt_number": 5 }, { "cell_type": "code", "collapsed": false, "input": [ "# plot the spectra given by the two methods\n", "fig, ax = plt.subplots(1, 1, figsize=(10,4))\n", "\n", "ax.plot(wlist1/(2*pi), spec1, 'b', lw=2, label='eseries method')\n", "ax.plot(wlist2/(2*pi), spec2, 'r--', lw=2, label='me+fft method')\n", "\n", "ax.legend(loc=1)\n", "ax.set_xlabel('Frequency')\n", "ax.set_ylabel('Power spectrum')\n", "ax.set_title('Vacuum Rabi splitting')\n", "ax.set_xlim(wlist2[0]/(2*pi), wlist2[-1]/(2*pi));" ], "language": "python", "metadata": {}, "outputs": [ { "output_type": "display_data", "png": "iVBORw0KGgoAAAANSUhEUgAAAmoAAAEXCAYAAAAdhmj9AAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3XlclWXawPHfYVV2RAHZBAUV3CBR3FLcldJKW6zJ0HRy\nfFMzm15zyrJsUqesMctJG1vMJVveXJJwS8pxBLRAzSUBZUdcAAVRlsPz/nHkJAoc0LMB1/fz4XM4\n59zPfV/PieDyXlWKoigIIYQQQgizY2HqAIQQQgghRO0kURNCCCGEMFOSqAkhhBBCmClJ1IQQQggh\nzJQkakIIIYQQZkoSNSGEEEIIMyWJmhBC3IH09HQsLCyoqqqq9f0lS5bw5z//Wa9txsXF4evrq33e\nvXt3fv7550bV4ejoSHp6ul7jEkIYjiRqQogGGTNmDK+99tptr2/dupX27dvXmbCYsylTpmBra4uj\noyNt2rRh+PDhHD9+XC91L1iwgI8//lgvddXlt99+Y/DgwQAsWrSIyZMn13g/MjKStWvX1nituLgY\nf39/g8YlhNAfSdSEEA0yZcoU1q9ff9vrX3zxBU8++SQWFk3v14lKpWL+/PkUFxeTm5uLn58fU6dO\nNXVYeqNSqUwdghDiLjW936xCCJN44IEHuHTpEvv379e+VlhYyI4dO3jqqadITEykf//+uLq64uXl\nxezZs6moqNCWPX78OCNHjsTNzQ1PT0+WLl0KaBLAhQsXasvdOrxnYWHBmTNntM9vLh8XF4ePjw9v\nv/027u7ueHl5sWXLFmJiYujcuTNubm7adnRp1aoVjzzySI0etR07dhAWFoazszN+fn68/vrrt123\ndu1avL298fLyYvny5drXa+vhqnbx4kXuv/9+XF1dcXNz0/aKAfj7+7N06VK6detGmzZtePrppykr\nK6u1Hn9/f/bu3UtsbCxLlixh8+bNODo6EhoayiuvvML+/fuZNWsWjo6OzJkz57bPc8qUKTz77LPc\nf//9ODk50a9fvxqf9a5du+jSpQsuLi48++yzDBky5LYeOiGEYUmiJoRokNatW/Poo4+ybt067Wtf\nffUVwcHB9OjRAysrK1asWMGlS5c4ePAge/fuZdWqVYBmuG3EiBFERUWRl5dHamoqw4cPBzS9Po3p\n+bm1fH5+PmVlZeTl5fHGG28wffp0NmzYQFJSEvv37+eNN94gIyOjzvqqT9G7evUqmzZtIiIiQvue\ng4MD69ev5/Lly+zYsYN//etfbN26tcb1cXFxpKamsmvXLpYtW8bevXu1cdZl+fLl+Pr6cvHiRc6f\nP8+SJUtqvL9x40Z27dpFWloap0+f5s0336z3sxgzZgx/+9vfmDRpEsXFxSQnJ/Pmm29y77338uGH\nH1JcXMz7779fax2bN29m0aJFFBYWEhgYyMsvvwxokslHHnmEZcuWUVBQQJcuXTh48KD00glhZJKo\nCSEaLDo6mm+++Yby8nIA1q1bR3R0NAD33HMPffv2xcLCgg4dOvDMM8/w008/AfD999/j5eXF888/\nj42NDQ4ODvTp00dbb2OPHL65vLW1NS+//DKWlpY89thjFBQUMHfuXOzt7QkJCSEkJITk5OQ663nn\nnXdwdXXFycmJ//73v3z11Vfa94cMGUK3bt0A6NGjB5MmTdLeU7XXXnuN1q1b0717d6ZOncqmTZt0\n3pONjQ15eXmkp6djaWnJwIEDte+pVCpmzZqFt7c3rq6uvPzyy9o6dX0mtbVZXxwqlYoJEyYQHh6O\npaUlf/rTn7SfVUxMDN27d+fBBx/EwsKCOXPm4OnpqTMOIYR+SaImhGiwgQMH0rZtW7777jvS0tI4\ndOgQTzzxBACnT5/m/vvvp3379jg7O/Pyyy9z6dIlALKysujYsaNBYnJzc9P28rRu3RoADw8P7fut\nW7fm6tWrtV6rUql48cUXKSwsJD09HVtb2xo9hgkJCQwdOhR3d3dcXFxYvXq19p6q3TxM6+fnR25u\nrs6YX3zxRQIDAxk1ahSdOnVi2bJld11nXXT1gN36WZWUlACQm5uLj49PjbK3PhdCGJ4kakKIRnnq\nqadYt24d69evZ8yYMbRr1w6AmTNnEhISQmpqKpcvX+bvf/+7diWon59fjblPN7O3t6e0tFT7/Ny5\nczXet7Ozq/F+Xl6eXoffqnucfH19ef/991m8eDHFxcUAPPHEEzz44INkZ2dTVFTEX/7yl9tWt2Zm\nZtb43tvbW2ebDg4OvPPOO6SlpbFt2zbeffdd9u3bV2edXl5eOuus7TO5m8/Jy8uL7Oxs7XNFUWo8\nF0IYhyRqQohGeeqpp9i9ezf//ve/tcOeACUlJTg6OmJnZ8epU6f417/+pX3vvvvuIy8vjxUrVlBW\nVkZxcTGJiYkAhIaGEhMTQ2FhIefOneOf//xnjfZCQ0PZsGEDarWa2NjYRu8bVp9bhwVHjBhBYGCg\ndm5dSUkJrq6u2NjYkJiYyMaNG29Lft58802uXbvG8ePH+eyzz3jsscd0trtjxw5SU1NRFAUnJycs\nLS21q2YVRWHVqlXk5ORQUFDA3//+dyZNmqSzTk9PT9LT02vck4eHB2lpaQ2+/5tFRUVx7Ngxtm7d\nSmVlJR9++OFtSbQQwvAkURNCNEqHDh0YOHAgpaWljB8/Xvv6O++8w8aNG3FycuKZZ55h0qRJ2qTG\n0dGR3bt3s337dtq3b0/nzp2Ji4sDYPLkyfTq1Qt/f3/GjBlT4zqAFStWsH37dlxdXdm4cSMPPfRQ\njXhuTZzuZmECaIYl33//fSoqKli1ahWvvvoqTk5OLF68+LYkTKVSMWTIEAIDAxkxYgQvvvgiI0aM\nqLPuaikpKYwcORJHR0cGDBigXVFZfd0TTzyhHRYNCgrilVde0Xl/jzzyCKAZCg4PDwfgueee45tv\nvqFNmzbMnTu3Qfdf/bxt27Z8/fXX/O///i9t27bl5MmThIeHY2trW/uHKYQwCJXS2Fm8QgghDCYg\nIIC1a9cybNgwU4dSQ1VVFb6+vmzcuFGbVAohDM+gPWqxsbF07dqVoKCg2ybLVpszZw5BQUH06tWL\npKQkAK5fv05ERAShoaGEhISwYMECbflFixbh4+NDWFgYYWFhxMbGGvIWhBCixdq1axdFRUWUlZXx\n1ltvAdCvXz8TRyVEy2JlqIrVajWzZs1iz549eHt706dPH8aPH09wcLC2TExMDKmpqaSkpJCQkMDM\nmTOJj4+nVatW7Nu3Dzs7OyorKxk0aBAHDhxg4MCBqFQq5s2bx7x58wwVuhBCCODgwYM88cQTlJeX\n061bN7Zs2SJDn0IYmcEStcTERAIDA7Vnyk2aNImtW7fWSNS2bdumnYwcERFBUVER+fn5eHh4YGdn\nB0B5eTlqtRpXV1ftdTJaK4Rors6ePWvqELRee+21Ws93FUIYj8EStZycnBp7Afn4+JCQkKCzTHZ2\nNh4eHqjVanr37k1aWpp22X+1lStXsm7dOsLDw1m+fDkuLi416pWds4UQQgjRlNTVCWWwOWoNTZZu\nDaz6OktLS5KTk8nOzubnn3/WrhCbOXMmZ8+eJTk5mfbt2/PCCy/UWW9z/HrttddMHkNT+ZLPSj4n\n+azkczLnL/mc5LOq/qqPwRI1b29vsrKytM+zsrJu29X61jLZ2dm3bRbp7OzMfffdx+HDhwFwd3fX\nLimfPn26di8mIYQQQojmxmCJWnh4OCkpKaSnp1NeXs7mzZtr7LkEMH78eO1xLfHx8bi4uODh4cHF\nixcpKioC4Nq1a+zevZuwsDBAsyt5te+++44ePXoY6haEEEIIIUzKYHPUrKys+OCDDxg9ejRqtZpp\n06YRHBzM6tWrAZgxYwZRUVHExMQQGBiIvb09n376KaBJxqKjo6mqqqKqqorJkyczfPhwAObPn09y\ncjIqlYqAgABtfS1FZGSkqUNoMuSzahj5nBpOPquGkc+pYeRzariW/Fk1yw1vVSqVzjFfIYQQQghz\nUF/eYrAeNSGEEKKladOmDYWFhaYOQ5gpV1dXCgoKGnWN9KgJIYQQeiJ/f0R96vr5qO/nRg5lF0II\nIYQwU5KoCSGEEEKYKUnUhBBCCCHMlCRqQgghhBBmShI1IYQQQuhVVFQUX3zxhanDaBR/f3/27t2r\nl7oiIyNZu3atXuqS7TmEEEIIoVcxMTGmDqFeU6ZMwdfXl8WLF2tfqz6eUh/0WZf0qAkhhBBCLxpy\nyLhoHEnUhBBCiBYiNzeXiRMn4u7uTseOHVm5cqX2vcTERMLDw3F2dsbT05MXXnhB+158fDwDBgzA\n1dWV0NBQfvrpJ+17kZGRvPLKKwwcOBAHBwfOnDlz29DfJ598QkhICG3atGHMmDFkZmZq33v++efx\n8PDA2dmZnj17cvz48Vpjj4yMZOHChQwcOBBHR0fGjx/PxYsX+dOf/oSzszN9+/YlIyNDW/7UqVOM\nHDkSNzc3unbtytdffw3AmjVr2LhxI//4xz9wdHTkgQce0F6TlJREr169cHFxYdKkSZSVlWnf+/jj\njwkKCsLNzY0HHnigxtnju3fvpmvXrri4uDB79mz9JqxKM9RMb0sIIYSZ0/X3B/T31VhqtVq55557\nlMWLFysVFRXKmTNnlI4dOyo7d+5UFEVR+vXrp6xfv15RFEW5evWqEh8fryiKomRnZytubm7KDz/8\noCiKouzevVtxc3NTLl68qCiKogwZMkTp0KGDcuLECUWtVisVFRVKZGSksnbtWkVRFGXLli1KYGCg\ncurUKUWtVitvvvmmMmDAAEVRFCU2Nlbp3bu3cvnyZUVRFOXUqVNKXl5erfEPGTJECQoKUs6cOaNc\nvnxZCQkJUQIDA5W9e/cqlZWVylNPPaVMnTpVURRFKSkpUXx8fJTPPvtMUavVSlJSktK2bVvlxIkT\niqIoypQpU5SFCxfWqL9Dhw5KRESEkpeXpxQUFCjBwcHKRx99pCiKouzdu1dp27atkpSUpJSVlSmz\nZ89WBg8erCiKoly4cEFxdHRUvv32W6WyslJ57733FCsrK+3936yun4/6fm6kR00IIYRoAQ4dOsTF\nixd55ZVXsLKyIiAggOnTp/Pll18CYGNjQ0pKChcvXsTOzo6IiAgA1q9fT1RUFGPGjAFgxIgRhIeH\ns2PHDkAzH2vKlCkEBwdjYWGBlVXN6e8fffQRCxYsoEuXLlhYWLBgwQKSk5PJzMzExsaG4uJiTp48\nSVVVFV26dMHT07PW+FUqFVOnTiUgIAAnJyfGjh1L586dGTZsGJaWljzyyCMkJSUB8P333xMQEEB0\ndDQWFhaEhoYyYcIEba+aUkuPl0qlYs6cOXh6euLq6sq4ceNITk4GYMOGDUybNo3Q0FBsbGxYsmQJ\nBw8eJCMjg5iYGLp3786ECROwtLRk7ty5dd7DnZBETQghhDASffapNVZGRga5ubm4urpqv5YsWcL5\n8+cBWLt2LadPnyY4OJi+fftqE7GMjAy+/vrrGtcdOHCAc+fOaev29fWtt93nnntOe62bmxugGYYd\nOnQos2bN4tlnn8XDw4MZM2ZQXFxcZ10eHh7a71u1aoW7u3uN5yUlJdo2ExISasS8ceNG8vPzAeqc\n6H9zgtW6dWuuXr0KQF5eHh06dNC+Z29vj5ubGzk5OeTl5eHj41Ojnvo+j8aSVZ9CCCFEC+Dn50dA\nQACnT5+u9f3AwEA2btwIwLfffsvDDz/MpUuX8PPzY/LkyaxZs6bOuutb4ejn58fChQt5/PHHa31/\n9uzZzJ49mwsXLvDoo4/y9ttv88Ybb+i8H11tDhkyhF27djX62tp4eXmRnp6ufX716lUuXbqEj48P\n7du3JysrS/ueoig1nt8t6VETQgghWoC+ffvi6OjIP/7xD65du4Zarea3337j8OHDgGaI88KFCwA4\nOzujUqmwtLTkySefZPv27ezatQu1Ws3169eJi4sjJydHW/etw4g3+8tf/sJbb73FiRMnALh8+bJ2\nCPLw4cMkJCRQUVGBnZ0drVq1wtLSss66bm6nvjbvu+8+Tp8+zfr166moqKCiooJDhw5x6tQpQNMz\nd+bMGV0fmbaNxx9/nE8//ZQjR45QVlbG3/72N/r164efnx9RUVEcP36c7777jsrKSt5///0avY13\nSxI1IYQQogWwsLDg+++/Jzk5mY4dO9KuXTueeeYZrly5AsDOnTvp3r07jo6OPP/883z55ZfY2tri\n4+PD1q1beeutt3B3d8fPz4/ly5fXSJTq66F68MEHmT9/PpMmTcLZ2ZkePXqwc+dOAK5cucIzzzxD\nmzZt8Pf3p23btrz44ot11nVzO7XtVVb93NHRkV27dvHll1/i7e1N+/btWbBgAeXl5QBMmzaNEydO\n4OrqyoQJE+psq7q+4cOHs3jxYiZOnIiXlxdnz57Vzu1r27YtX3/9NS+99BJt27YlNTWVQYMG1XkP\njaVS6ktJmyiVSiX7uAghhDA6+fsj6lPXz0d9PzfSoyaEEEIIYaYkURNCCCGEMFOSqAkhhBBCmClJ\n1IQQQgghzJQkakII81JcDNevmzoKIYQwC5KoCSHMQlkZJCRA7Mwt5LftxpSHSyRfa6D8fPj8c3js\nMRg2DG5sNC+EaAYMmqjFxsbStWtXgoKCWLZsWa1l5syZQ1BQEL169dKe0XX9+nUiIiIIDQ0lJCSE\nBQsWaMsXFBQwcuRIOnfuzKhRoygqKjLkLQghjGDfPnBzg3794IcNBXhcPUPHb//B7t2mjsy8ZWZC\n//7g6QlTpsBXX2k+y82bTR2ZEEJfDJaoqdVqZs2aRWxsLCdOnGDTpk2cPHmyRpmYmBhSU1NJSUlh\nzZo1zJw5E9Cc17Vv3z6Sk5M5evQo+/bt48CBAwAsXbqUkSNHcvr0aYYPH87SpUsNdQtCCCP56iu4\nehX8/cHzvt4A/JV3SDtcaNrAzNzXX0N8PLRqBVFRsKzfdxymN2f3pZs6NCGEnhgsUUtMTCQwMBB/\nf3+sra2ZNGkSW7durVFm27ZtREdHAxAREUFRUZH2wFQ7OzsAysvLUavVuLq63nZNdHQ0W7ZsqbX9\nRYsWab/i4uIMcYtCCD359VfN4yefwILvB5EbNBg7rmER96NpAzNz537NJYJ43l5SyY4dMM35G3rz\nKx32rzd1aEJo5efnM3jwYJycnLSnDkydOpU2bdrQr18/E0enYWFh0aAjpRrC39+fvXv31lsmLi6u\nRp5SH4Mdyp6Tk1Pj9HgfHx8SEhJ0lsnOzsbDwwO1Wk3v3r1JS0tj5syZhISEAJr/4B4eHoDmrK7q\nxO5Wum5cCGEeKirgyBHN92FhN16LHAUpP+N1Yg8w0WSxmbsOB7/kbV4ge9efYe4anJ5+GHZu5J6L\nO7ly5RWcnEwdoRCwZs0a3N3dtUdV7d+/nz179pCbm0urVq1YtGgRaWlpfPHFF0aJJzIyksmTJzNt\n2jSD1F/b0Va1xRAZGal9/vrrr9dZ1mA9ag09mf7WIxOqr7O0tCQ5OZns7Gx+/vnnWnvFGvJhCCHM\n28kTCsPKYugRUIKLi+Y110dGAND9UhxVVSYMzsx1zdkDQKvRgwGwHj2MSpUV/TnIkZ9k/q4wns8+\n+4ypU6fW+l5GRgbBwcE1nvv7+9OqVStjhVdDU8sbDJaoeXt7k5WVpX2elZWFj49PvWWys7Px9vau\nUcbZ2Zn77ruPX375BdD0olWfSp+Xl4e7u7uhbkEIYQRntx0jhvvYk98DbvzDzWlob6Jdt9FXSSAj\nw8QBmqnC/HL6l/8EQJuHh2tedHbmbPsBWKHm0lf1D72Ilsnf35933nmHnj174ujoyLRp08jPz2fs\n2LE4OzszcuTIGov04uPjGTBgAK6uroSGhvLTTz/VWm9dyc+UKVNYt24d//jHP3B0dGTNmjX8+c9/\n5uDBgzg6OtKvXz+WLFnC5s2bcXR0JKy6W91Acb/88svs37+fWbNm4ejoyJw5c7TX7N69m86dO+Pq\n6sqsWbO0ryuKwptvvom/vz8eHh5ER0drewcBvvjiCzp06EDbtm156623GvBfoZEUA6moqFA6duyo\nnD17VikrK1N69eqlnDhxokaZHTt2KGPHjlUURVEOHjyoREREKIqiKBcuXFAKCwsVRVGU0tJS5d57\n71X27NmjKIqivPjii8rSpUsVRVGUJUuWKPPnz7+tbQPelhBCz74evEJRQDl6z1M1Xh8+XFFAUb7/\n3kSBmbkjXxxRFFAybAJrvP7LhDcVBZQdnWaZKLKWrUF/fzT/JLn9qzHl75C/v7/Sv39/5fz580pO\nTo7i7u6uhIWFKcnJycr169eVYcOGKa+//rqiKIqSnZ2tuLm5KT/88IOiKIqye/duxc3NTblw4cJt\n9X722WfKlClTam1zypQpysKFC2uUHTRokPb5okWLlMmTJxs87osXLyqKoiiRkZHK2rVra9SvUqmU\ncePGKZcvX1YyMzOVdu3aKbGxsYqiKMratWuVwMBA5ezZs0pJSYkyYcIEbbzHjx9XHBwclP379ytl\nZWXKvHnzFCsrK2Xv3r213kddPx/1/dwYrEfNysqKDz74gNGjRxMSEsJjjz1GcHAwq1evZvXq1QBE\nRUXRsWNHAgMDmTFjBqtWrQI0PWXDhg0jNDSUiIgIxo0bx/Dhmn8xvvTSS9qs98cff+Sll14y1C0I\nIYzA4eQhACwH1pxUfGNaKidOGDuipuHy/qMA5Ln3rPG6w8zJRBDP3PK3TRGWaAJmz55Nu3bt8PLy\n4t5776V///706tULW1tbHnroIe1WWevXrycqKooxY8YAMGLECMLDw4mJibmtTuWWaUz1vX9rWUVR\ndF6vj7h37NhRb7wvvfQSTk5O+Pr6MnToUI7cmDy7YcMGXnjhBfz9/bG3t2fJkiV8+eWXqNVqvvnm\nG8aNG8egQYOwsbFh8eLFWFjoN7Uy2GICgLFjxzJ27Ngar82YMaPG8w8++OC263r06MGv1cvAbtGm\nTRv27NmjvyCFECajVkPARU2i5jW+T433unXTPB4/buyomoa0wjYUMxa6Danxeqehfhy39+Nqlmbj\nW5kdYoYakJTcVXkdqhfkAbRu3brG81atWlFSUgJo5pJ9/fXXbN++Xft+ZWUlw4YNA+B//ud/2LRp\nE6DZoaGyslK7E0OHDh1ITk42y7ih9qFaT09P7fd2dnba+vLy8ujQoYP2PT8/PyorK8nPzycvL6/G\ntC47Ozvc3Nzu5jZvY9BETQgh6pP6y2W6KL9TjjUu9/ao8Z70qNVvW2UU3xHFxuiar1taQu/e8PPP\ncOgQ3HefaeITTUddvVl+fn5MnjyZNWvW1Pr+qlWrtCNhn3/+OT/99BOffPJJo9u/0x6oO427sYsJ\nvLy8SE9P1z7PzMzEysoKT09P2rdvX2OP2NLSUi5dutSo+nWRI6SEECZzPKGYzTzKL+3vB1vbGu9V\nJ2qpx8uoKq80QXTm7dQpzWOXLre/1+dG5+ShQ8aLRzQ/Tz75JNu3b2fXrl2o1WquX79OXFwcOTk5\nt5Wtb/hS17Cmh4cH6enpDRr+1EfcHh4epKWl1VvHzffz+OOP895775Genk5JSQl/+9vfmDRpEhYW\nFkycOJHvv/+eAwcOUF5ezquvvkqVnpeqS6ImhDCZA+k+TGIze2b+323vubnB97YTOF9qT37MLyaI\nznxVVkJqqub7zp1vf79vX81jYqLxYhJN1809TDdve+Xj48PWrVt56623cHd3x8/Pj+XLl9eaiNS3\nXdat7936/JFHHgHAzc2N8PBwg8VdnXg999xzfPPNN7Rp04a5c+fqjPnpp59m8uTJDB48mI4dO2Jn\nZ8fKlSsB6NatGx9++CFPPPEEXl5etGnTpsb+sPqgUvSVwpoRlUqlt8xcCGE4Q4dCXBx8/33tQ3R7\nPZ5g+PlNHH3+E3q+W/seTS1RaioEBYGPD9y0w5HW2bPQsSP4upWSccGOJrZtVJMmf39Efer6+ajv\n50Z61IQQJlFV9cfRUb17116m3C8QgKtH6x+maGnqG/YE8G9bQqZFB3675EnGWdkxWIimTBI1IYRJ\nnDkDV65A+/Zw02KrGmxCOgFgcUYStZuVfv8jU/iUgV5na31f5eiAvVUZThRz/IdMI0cnhNAnWfUp\nhDCJ6t60e+6pu4xbeEdYBw7n9XNYcnPht/cTPmUDe9RrgYBayxR496DN2Xzyf/wNnvU3anwtmaur\na5M7okgYj6ura6OvkURNCGES52J+5Vn+S7D3YKBnrWV8IzU9aralhSgKMtfqBod8TQ+jc1inOstY\n9OgOZ/egTjoG3G+kyERBQYGpQxDNjAx9CiFMwmX/dj5gNiMvbKizjFv39gS6FRKknCY724jBmTmP\nEk2i5jWoY51l3IZ0B8Ap6zfUaqOEJYQwAEnUhBBGpyjgmK3Zydbt3m51F1Sp8OnuAsgJBdUKM4tp\np1zgOra0D/eus5xz/xCqUOFYWVjrylAhRNMgiZoQwugyM6FjuWbpYpsBXestKycU1JT1k2a+Xq5t\nABZW9fwK79OH0QOvch8xpKQYKTghhN5JoiaEMLqTx6vozGkAVF3r2GPihuozPyVR00jJd+Jt/kpi\n5yfrL2hlRYeurQE4fdoIgQkhDEIWEwghjO5cYiatuc7l1p44OzvXW7a6R02GPjUOXwpgKW/z6kO6\ny1afWiA9akI0XZKoCSGMLiWrFa/zKiPutWCgjrIhIWDLdS4dLwLq2HCtBfn9d81j1/pHjAHN6QUg\nPWpCNGWSqAkhjO7XXE9ieZ3Q/9Fd1j37V67Tm6TiUEpKknBwMHx85kzXqQQ3kx41IZo+maMmhDC6\n6sShusenPipvLwD8yCQvz4BBNQG6DmO/VadO4MxlLM6kUlFh2NiEEIYhiZoQwqgqKiA9XbN5bce6\ntwH7g7s75Sob3Cgg/2ypocMza9nZms/Py4sG9Sy2OnuSIlzYXhXF2dpPmxJCmDlJ1IQQRpWeDmo1\n+PpCq1YNuMDCggI7HwAu/9ayNwTLy6zgLRYwy3p1wy4ICKAKFR05Q8rxcsMGJ4QwCEnUhBBG1Zhh\nz2rFLr4AXD/dsg8YLzyeywKW8kz+Gw27oFUrChw7YIWaC4nSpSZEUySJmhDCqPLjz7KSWfxJ+aLB\n11zz8CdlWz8+AAAgAElEQVSX9lzJv2bAyMxf6e+aRPWyk2+Drynx1qw6KE2WpZ9CNEWSqAkhjKri\ncDKz+JAheV82+Jqk2Z/iTS577ccbMDLzV3lGM/R7rZ1fg69RgjSrDixTfzdITEIIw5JETQhhVBZp\nN5YtBgY2+Jr2XioAcnMNEVHTocrVnExf5eXT4Gvsw0NIoyN5F60NFZYQwoAMmqjFxsbStWtXgoKC\nWLZsWa1l5syZQ1BQEL169SIpKQmArKwshg4dSrdu3ejevTvvv/++tvyiRYvw8fEhLCyMsLAwYmNj\nDXkLQgg9s8/TTFKz69nwSWpemh06Wvz2HNYXNB+Apa9Xg69xXfAXulim8XrRc1xr2SPHQjRJBtvw\nVq1WM2vWLPbs2YO3tzd9+vRh/PjxBAcHa8vExMSQmppKSkoKCQkJzJw5k/j4eKytrXnvvfcIDQ2l\npKSE3r17M2rUKLp27YpKpWLevHnMmzfPUKELIQykvBw8izWJWpt+DU/U2rfXPLb0RO0HiyiO48Cj\ngwc1+Bpra802KCkpkJYG3bsbMEAhhN4ZrEctMTGRwMBA/P39sba2ZtKkSWzdurVGmW3bthEdHQ1A\nREQERUVF5Ofn4+npSWhoKAAODg4EBweTk5OjvU5RFEOFLYQwoLNnIQhNomYT3PChzzZtwMYGLl+G\n0ha8ldp3xSN5lcU4j+7XqOvkKCkhmi6D9ajl5OTg6/vHyiQfHx8SEhJ0lsnOzsbDw0P7Wnp6OklJ\nSURERGhfW7lyJevWrSM8PJzly5fj4uJyW/uLFi3Sfh8ZGUlkZKQe7koIcTdSUmAd7xLV6TRTOnRo\n8HUqFXR1L0CdnUtedjCdOlsaMErzVFYGly6BhQW0a9e4azt3hpgYOUpKCHMRFxdHXFxcg8oaLFFT\nqVQNKndr79jN15WUlPDwww+zYsUKHG5swz1z5kxeffVVABYuXMgLL7zA2rVrb6v35kRNCGEeUlPh\nax6l3WiY0sjfPnvzu9GWcyT8lk2nzt6GCdCMnTunefT0BMtG5qnSoyaEebm1A+n111+vs6zBhj69\nvb3JyvpjF/GsrCx8fHzqLZOdnY23t+YXcEVFBRMnTuTJJ5/kwQcf1JZxd3dHpVKhUqmYPn06iYmJ\nhroFIYSe3clmt9Wu2Gsmql35vWUu/ayen1c9X68xgr2v0I+DqJOP6TcoIYTBGSxRCw8PJyUlhfT0\ndMrLy9m8eTPjx9fcA2n8+PGsW7cOgPj4eFxcXPDw8EBRFKZNm0ZISAhz586tcU3eTbOJv/vuO3r0\n6GGoWxBC6Fl1otaInTm0Sp01Gcq1tJa5oqB6axKvhi/41Op5fBMHGcDYE8v1G5QQwuAMNvRpZWXF\nBx98wOjRo1Gr1UybNo3g4GBWr9acUTdjxgyioqKIiYkhMDAQe3t7Pv30UwAOHDjA+vXr6dmzJ2Fh\nYQAsWbKEMWPGMH/+fJKTk1GpVAQEBGjrE0KYv7vpUato5wUZUJnVMnvUried5D1W41jSB/hTo651\n7av5wH2up3DlCjg5GSBAIYRBqJRmuIRSpVLJylAhzExZGdjZab6/dk2zirMxkse/Suj2xWzpsZAH\njzbwrMtm5MsJXzHpu8c4ETyBkBPfNu7irCzw8+M87cg6fJ7evQ0ToxDiztSXt8jJBEIIozhzBr6s\neoT19s9gU1bc6OstgwJIpRMXSu0NEJ35q8rRDPlWedzBJDVvb8osWuHOBdKPXNZzZEIIQ5JETQhh\nFGd/u8ojfMPDJZ9B69aNvl6JnkoQqfzTdr7+g2sCrKpPJfC5g0TNwoKCNpqJgQWJqfoMSwhhYJKo\nCSGM4mKCJkEocO4IVo2fHls9ib6lnvfZqkiTqLUOuINEDbgUOpxtjCMjS37tC9GUyP+xQgijuHZM\nk6hd9b6DlQSAm5vmOKSiIlrkmZVOJZpEzbHznSVqBQv/yQNsY8+lMH2GJYQwMJ3/rK2srGTHjh2k\np6dTWVkJoD1vUwghGsoyTbPkU3UnSz7RnE7g6amZF5+Xpzm/sqWoqIB3K2axk2EsHnhnh3XKprdC\nNE06E7Vx48bRunVrevTogYWFdMAJIe6MwzlNomYfemeJGmiGP1tionbuHGxnPIc8x7Mk4M7q8PQE\nBwcoLNQcReXmpt8YhRCGoTNRy8nJ4ejRo8aIRQjRTF27BguvLuBbi/vY+PidD70Fu56jikwupgTB\nQFc9Rmje7uZUgmoqlaZXLSlJ06vWv79+YhNCGJbOLrJRo0axc+dOY8QihGim0tIglUCOBk7AuvMd\ndgkBL5ycRiIRWBzYr8fozF91onYnpxLcrHNnzaMczi5E06GzR23AgAE89NBDVFVVYW1tDWjmqF25\ncsXgwQkhmofff9c8VicKd6qibfXpBC3rGKnqla5306MGMNjhV+xJIv/QSHjK7+4DE0IYnM4etXnz\n5hEfH09paSnFxcUUFxdLkiaEaJTqCex3m6jhpclUVHkta48OfQx9Atx35O+sZTo2hw/cfVBCCKPQ\nmaj5+fnRrVs3WUgghLhj1Ylaly53V4+1nyZTsbnUshK1Ngk/sIY/0//CtruqxypYs5DDOl3GPoVo\nKnQOfQYEBDB06FDGjh2LzY3D+WR7DiFEY+irR80uUDNJy/5yyxr6dD9zkCf4N6cKvYDxd1yPc3gQ\nfAGuF1JQFM0CAyGEedPZTRYQEMCwYcMoLy+npKREO/wphBAN9eSv89jPILpfjLurelx6duAIPUmp\n8NdLXE1Fq0JNYtrK3/Ou6nG4sTVKgDpFO5wqhDBvOnvUFi1aZIQwhBDNVUEB9LyewED+i9JGuau6\nXCJD8bA6QmUZTC4DW1s9BWnmHIrPaR6D7nKSWqDmvM8gUvgt5e5XkQohDE9nojZ06NDbXlOpVPz4\n448GCUgI0bycPq1JDABUne98s1sACwvNxq3Z2ZoJ9v7+egjQzFVWQptyTfeXS9e7TNTatye+w2P8\nlNEBt5OVDBnS+DNXhRDGpfP/0rffflv7/fXr1/n222+xuoMDlYUQLdPZ5Mv04wJllq2x1UMXTvv2\nLStRy8+H9mgSNSvfu0zUVCr2TP+ShQvhxTN6CE4IYXA6M67w8PAazwcNGkSfPn0MFpAQonkpPKTp\nTSt0C8RTD6vHq3O93Bay8DMvD95iJRFe2cz3vLs5avDHmZ+y6a0QTYPORK2goED7fVVVFYcPH5Z9\n1IQQDVb+myYjKO9wd8Oe1ar3Emspk+Hz8uA7JnC9F8y3ufv65HB2IZoWnYnaPffcg+rGGm4rKyv8\n/f1Zu3atwQMTQjQPX5RO5CNO8vWL+qmvs30OkZzm2vEOQPM/mV1fpxJUq07U0tJArQZLS/3UK4Qw\njHoTtaqqKjZs2MDAgQONFY8QohmpqoKTaTZcoyu+I/VT57DfVvA8b/NN4lvAAv1Uasb0dSpBNUdH\nzYKMc+cgK6tlzPMToimrd8KIhYUFzz77rLFiEUI0Mzk5cO0auLuDi4t+6rS+MaHe6mLLGPvU14Hs\nWpWVvOjwL95igcxTE6IJ0Dmzd8SIEXzzzTcoyt3tfySEaHn0dsbnTew6aibU2xW1jERN30OfWFry\nPxn/ywKWkpF0SU+VCiEMRWei9tFHH/Hoo49iY2ODo6Mjjo6OODk5GSM2IUQT9/vvmse7PePzZk5d\nNBmLU+k5/VVqxu49spL1/Iku5/frp0KVisvtNBvfFv8iXWpCmDudiVpJSQlVVVVUVFRoj49q6KrP\n2NhYunbtSlBQEMuWLau1zJw5cwgKCqJXr14kJSUBkJWVxdChQ+nWrRvdu3fn/fff15YvKChg5MiR\ndO7cmVGjRlFUVNSgWIQQxpfyexWg3x415xubvratzKO8XH/1mqtuF+L4Extxr9TffiQV/poVBepT\nkqgJYe50JmrDhw9v0Gu3UqvVzJo1i9jYWE6cOMGmTZs4efJkjTIxMTGkpqaSkpLCmjVrmDlzJgDW\n1ta89957HD9+nPj4eD788ENOnToFwNKlSxk5ciSnT59m+PDhLF26tEE3KoQwPu+4DVyiDQ8k/k1v\ndVr6tOcX637E049zzbxTTa0G1+uaIV7nuz2V4Ca23TWZc+us3/VWpxDCMOpM1K5du8alS5e4cOEC\nBQUF2q/09HRycnJ0VpyYmEhgYCD+/v5YW1szadIktm7dWqPMtm3biI6OBiAiIoKioiLy8/Px9PQk\nNDQUAAcHB4KDg7Vt3nxNdHQ0W7ZsubM7F0IYnH3mSdpQSBt3a/1V6uTEjJ4Hmcz6Zr/p7YULf5xK\nYO2nv0TNuV9XADwKf6eiQm/VCiEMoM7tOVavXs2KFSvIzc2ld+/e2tcdHR2ZNWuWzopzcnLw9fXV\nPvfx8SEhIUFnmezsbDw8PLSvpaenk5SUREREBAD5+fna9z08PMjPz6+1/ZsPk4+MjCQyMlJnzEII\n/SkvB88iTS+6y43EQF+8vOCXX5r/pre5OQoh6Hl/DsDm3n6sdH6Z2Mv96XlWv0PTQgjd4uLiiIuL\na1DZOhO1uXPnMnfuXFauXMns2bMbHUT1Jrm63Lqa9ObrSkpKePjhh1mxYgUODg61tlFXOzcnakII\n40tLg65opixY9wzWa93VOUuz71FLLaIVZZRaOmBXy+/AOxYYyLY+b7JnD8w8LYmaEMZ2awfS66+/\nXmdZnXPUVCoVhYWF2ueFhYWsWrVKZxDe3t5kZWVpn2dlZeHj41NvmezsbLy9vQGoqKhg4sSJPPnk\nkzz44IPaMh4eHpy7MTElLy8Pd3d3nbEIIYwv5UQFgaRShUrvmUD1nmLNvUct51JrotjBZ4P0fxpM\n9X8SOUpKCPOmM1H7+OOPcXV11T53dXVlzZo1OisODw8nJSWF9PR0ysvL2bx5M+PHj69RZvz48axb\ntw6A+Ph4XFxc8PDwQFEUpk2bRkhICHPnzr3tms8//xyAzz//vEYSJ4QwH+cOZaGgotCpA9jZ6bXu\nlnLeZ/bFVvxAFNkDHtV73V1vjEafOKH3qoUQeqQzUauqqqKqqkr7XK1WU9GA2adWVlZ88MEHjB49\nmpCQEB577DGCg4NZvXo1q1evBiAqKoqOHTsSGBjIjBkztD11Bw4cYP369ezbt4+wsDDCwsKIjY0F\n4KWXXmL37t107tyZH3/8kZdeeumOblwIYViJFztiRylb//eA3usOsM0lih3Ynzyk97rNid5PJbhJ\njx6ax99+03/dQgj9USk6jhz461//SmZmJjNmzEBRFFavXo2fnx/Lly83VoyNplKp5CQFIUysXz9I\nSIAff4ShQ/Vbd/rCf+P/5p/Z5hrN+ILP9Fu5GXngAdi2Db75BiZO1G/dFy9Cu3Zgbw9XroCFzn+2\nCyEMpb68pd5D2QGWLVvGmjVr+Ne//gXAyJEjmT59un4jFEI0K1VVf/TUVPfc6JNjkGbs06GkeW+k\nZsgetbYVeXxst4SrV1VkZKwgIED/bQgh7p7ORM3S0pLo6GiGDh1K1676XWIvhGie0tPh6lXw9IS2\nbfVff/Xmr24VeVRUgLUet2kzJ3o/5/NmlpZML11JEc78fPSfBAQ0bKW+EMK4dHZ2b9u2jbCwMMaM\nGQNAUlLSbYsChBDiZseOaR4N0ZsGYOWjOZi9PXnN9nSCqir4W+4s/o+HaH/5lP4baNeOUltXXLjM\n2fhm+iEK0QzoTNQWLVpEQkKCduVnWFgYZ86cMXhgQoim6+Sv1/Akjx7dDTRX1N2dKlS05SLnsprn\n1voXL8IwZQ8PsQVbK7X+G1CpKPbV7G93JcEAiaAQQi90JmrW1ta4uLjUvEhmnQoh6lER9x/y8GLe\n7rGGacDKioOeD7GRJzifed0wbZhYbu4fx0cZZuwTLEI001lUp07qKCmEMBWdGVe3bt3YsGEDlZWV\npKSkMHv2bAYMGGCM2IQQTZTFjT/8Np07GKyNz8d9y2TWk1noaLA2TCn/bCnOXKFcZQM37WWpT04R\nmh4113MnKS83SBNCiLukM1FbuXIlx48fx9bWlscffxwnJyf++c9/GiM2IUQTdP06tL2g2UXVuZ9+\nj466WXPf9PbKac28sSt2ntDAI/kay/ah+3i+3XreU56TEwqEMFM6V33a29vz1ltvMX/+fFQqFU5O\nTsaISwjRRJ08CT2UowDY3GOg1QQ0//M+r6ZqMtBSZ8MMewIQHMzZAcGkbdVsp9K9u+GaEkLcGZ09\naocOHaJHjx707NmTHj160KtXLw4fPmyM2IQQTdCxI1V058Ymaj17Gqyd5n7e5zF6EMk+Eh9+26Dt\nVCdn1St1hRDmRWei9vTTT7Nq1SoyMjLIyMjgww8/5OmnnzZGbEKIJijtcCFn6EiRs59m63sDae49\namcuOvETkSiD7jVoO3KUlBDmTWeiZmVlxb33/vGLYtCgQVhZ6RwxFUK0UAmpboSRzL61Zw3ajo/t\nBR7hK0LSdxi0HVMx5KkEN5MeNSHMm86zPufOncu1a9d4/PHHAdi8eTOtWrVi8uTJANxzzz2Gj7KR\n5KxPIUzH21vTy5WaCp06Ga6dyrj9WA0dzH/pT9+K/9Lc/v3YoQNkZkJaGnTsaLh2Kio0531WVEBx\nMTg4GK4tIUTt6stbdCZqkZGRqG5acaQoSo3n+/bt01OY+iOJmhCmcemS5sgooxz0nZICnTtzhgBs\ns8/g7W3AtoxMUcDWVpM8lZZC69YGbCwjg/hu07hy1QKXhF307WvAtoQQtbqrQ9nj4uL0HY8Qopmq\nHj7r1s3ASRpoJ6m1J4/juQre3s3nrMpLlzRJmouLgZM0AFdX+l3dSxk2bEiqoG/fZnpwqhBNlM5f\npStWrODKlSsoisK0adO455572LlzpzFiE0I0MYY+47MGBwdKLR1ozXXOp14xQoPGk5sL+4hke8Vo\nTdekITk5UeTqjy3lnD+QYti2hBCNpjNRW7t2LU5OTuzatYuCggLWrVvHSy+9ZIzYhBBNTGZCHiPZ\nRV//fKO0V2yvOZy9+HTz2qPjXFYFg/mZ/lf3aMaRDex6kCazrvz1qMHbEkI0js5ErXrMdMeOHUye\nPJnusiOiEKIOTvE72cVoHtg7xyjt/d59Iqt5htxLtkZpz1gun8rDAoUrrT3A0tLg7bXqq9nvziFd\nln4KYW50Jmq9e/dm1KhRxMTEMGbMGK5cuSKHsgshblNVBc6Zmj/0Dv0Nt9HtzY4/uZS/sJqT1wOM\n0p6xlKbkAFDs7GOU9pwGanrUOl49xoULRmlSCNFAOhcTrF27luTkZDp16oSdnR2XLl3i008/NUZs\nQogmJCMDgis0Q2f2/YwxSa35nvdZmaFJ1MrbGWcpq8XoUUT3TOLLo8FsT4JRo4zSrBCiAXR2jVla\nWtK7d29cXFwAcHNzo6cBj4URQjRNR48o3MOvmiehoUZps3oz2OZ2OoEqT5OoKe2NtOeIqyuuQ0Mp\nx5aEBOM0KYRoGBnDFELoReqedNwooMSuHfj6GqXN5tqjttnmKfqSwOVo48z1A4iI0DxKoiaEeWlm\ne3kLIUzl6K+VfMGTDLrXAQeVcfY089Qs+iQ/H9Rqo8y7N4rTF1xJpy8uRtx89uZETVHASP8JhRA6\n1HsyQWVlJd26deP33383Zkx3TU4mEMK41GpwddUcQZSX90cCZXDl5cxquwmr4kLm587V9rA1ZYoC\nrVpBeTmUlBhldw5tu+7ucPGi4Y+tEkLUVF/eUu/Qp5WVFV27diUjI+OOGo6NjaVr164EBQWxbNmy\nWsvMmTOHoKAgevXqRVJSkvb1p59+Gg8PD3rcsnPmokWL8PHxISwsjLCwMGJjY+8oNiGE/pw4oUnS\n/P2NmKQBqFS8XzyVd5lHbmaFERs2nMJCTZLm5GS8JA00PWgREdCO8zL8KYQZ0TlHraCggG7dujFs\n2DDGjRvHuHHjGD9+vM6K1Wo1s2bNIjY2lhMnTrBp0yZOnjxZo0xMTAypqamkpKSwZs0aZs6cqX1v\n6tSptSZhKpWKefPmkZSURFJSEmPGjGnIfQohDCg+XvPYr5+RG7a2psjWAwsUCk+cM3LjhlE9387o\nvYOKwqeJIZzHg+P7jLNhsRBCN51z1BYvXnzba6oGTF5ITEwkMDAQf39/ACZNmsTWrVsJDg7Wltm2\nbRvR0dEAREREUFRUxLlz5/D09OTee+8lPT291rplWFMI82KyRA247OhDm7JzFJ/KAYyziMGQqlew\nVq9oNRqVCkvPdnDhJKU//wJEGTkAIURtdCZqkZGRpKenk5qayogRIygtLaWyslJnxTk5OfjetPLL\nx8eHhFv602srk5OTg6eOsZOVK1eybt06wsPDWb58uXbrkJstWrSoxj1ERkbqjFkIcWdMmaiVuXnD\nxcNcS8k2fuMGUHroOGmMIyd9ELDOqG3bDQ6HYz/jknqYsrIobJvXgQ9CmI24uDji4uIaVFZnorZm\nzRo+/vhjCgoKSEtLIzs7m5kzZ7J37956r2tIrxvc3jum67qZM2fy6quvArBw4UJeeOEF1q5de1u5\nmxM1IYThFBVBrxMbGWuZT6jrQ4C/UdtXvL3hd1Bn5hi1XUMpS8mkI2epVHcyetutBobDh9BL/QtH\njkBfI646FaIlubUD6fXXX6+zrM45ah9++CH/+c9/cHJyAqBz586cP39eZxDe3t5kZWVpn2dlZeHj\n41NvmezsbLy969/g0d3dHZVKhUqlYvr06SQmJuqMRQhhOIcOwTOs4R31PGxTjhu9/crIkbzHXBLK\njbPJrqFpTyXwMM7xUTWEh2seOCwLCoQwEzoTNVtbW2xv6v+urKxsUG9ZeHg4KSkppKenU15ezubN\nm29bhDB+/HjWrdN07cfHx+Pi4oKHh0e99ebdtLPld999d9uqUCGEccUfUP9xIkHv3kZv3/5PDzKP\n99hWNMTobRuCKkczhGvlZ6RTCW7WqRPX7VwpxJUj/yk2fvtCiNvoHPocMmQIf//73yktLWX37t2s\nWrWKcePG6a7YyooPPviA0aNHo1armTZtGsHBwaxevRqAGTNmEBUVRUxMDIGBgdjb29c4Q/Txxx/n\np59+4tKlS/j6+vLGG28wdepU5s+fT3JyMiqVioCAAG19QgjTyNtzHCeKudrWD3uj7s2h4eOj2Voi\nJwcqKsDa2ugh6JXNBU2Pml2QCRI1CwtO7s3jnv62BP5q/OaFELerd8NbgKqqKv7973+za9cuAEaP\nHs306dMbPAfNFGTDWyGMQ1HgRftVvHPtWa4++Cfsv1tvkjh8fDSJ2pkzEBBgkhD0QlFgt3UUo9Q/\nULxhG45P6P5Hsb5VVGj2cLt+XbP5rZub0UMQosWpL2/R2aO2b98+Jk+ezDPPPKP3wIQQTVtKCoRd\nOwCA3ciBJoujQwdNopaR0bQTtcJCeES9iS522STcZ5qtRqytNSPYBw5AYiKMHWuSMIQQN+ico/b5\n55/Tq1cvIiIiePHFF9m+fTuFhYXGiE0IYebi42Els9kYvBjVqJEmi6NDB83jHR6iYjbS0+EKzlzv\n1A2Vs5PJ4pAD2oUwHzp71Kon++fm5vLNN9/w7LPPkpub26C91IQQzdvBg5BAP7Kn9INA08Vxf+lX\nhHCUwiPPAH6mC+QuVe/xXZ14mookakKYD52J2hdffMF//vMfjh49Srt27Zg1axaDBg0yRmxCCDNX\nvdFt//6mjWPw7x/zBHt47/ggmnKiVt0jeONAF5Pp36uUSBJQ/ccStXowlpamjUeIlkxnojZ37lw6\nderEzJkziYyMJKApTwARQuhNSQkcOwaWlibZlaMGlY83nIKqzKZ9OkF1j5qpEzXfozvYx6PsLBnF\nL78Mlo1vhTAhnXPULl68yCeffML169d5+eWX6du3L08++aQxYhNCmLF9+0Cthj59wM7OtLHYdNaM\nFdrkZ5o2kLuUflaz6svUiRoDNQtD+nOQ3bFqEwcjRMumM1ErLi4mMzOTjIwM0tPTKSoqwsJC52VC\niGZu507N4+jRpo0DwKm7JlFzvZxOVZWJg7kLIxP/Ti7t6ZP4oWkD8fLiqkcAThSTtuWYaWMRooXT\nmXENGjSI7du307NnT7766itOnz6tXWAghGi5rv7fTo7QkykFy00dCrZd/AHwq0rn3DnTxnKnFAWc\nCs7SnnO08dA5K8XgrCM1c5Gdj/xEsRxSIITJ6PxtcPToUUDTs2bOm9wKIYwnLQ265+2iJ8eosi8w\ndTgQEsKnni+x+1wPZmeAl5epA2q8oiLwqtCsJrAPNvGyT8BmzDDY/AVDq/YSF/ccDTiQRghhADp7\n1I4dO0ZYWBjdunUjJCSE3r1789tvvxkjNiGEmdq5E0ayGwCLUSNMHA3g6cn3A5awiSea7F5q6eng\nTzoAqgB/0wVSbcQITnV5gC08yI2DaYQQJqAzUXvmmWd49913yczMJDMzk+XLl8spBUK0cAlbz9GT\nY1TY2MGAAaYOB/hjAn6TTdTS1PhxYzGEqTdSA/DxoeCTLXzK05KoCWFCOhO10tJShg4dqn0eGRnJ\n1atXDRqUEMJ8lZeD1U97AVAPGAy2tiaOSKOpn05w4bd8bKjgip0HtG5t6nAA6NtXc+7n6dNN93MV\noqnTmagFBASwePFi0tPTOXv2LG+++SYdO3Y0RmxCCDN08CB0KTsCQKv7zWDY84bqRK16L7Km5kSR\nF7ZcZ9Pzh0wdipaVFQwbpvl+927TxiJES6UzUfv00085f/48EyZMYOLEiVy4cIFPPvnEGLEJIczQ\nzp0wn3/wxrQMiI42dThaTb1HLT0dyrHFLdQ0h7HXZdQozaMMfwphGnWu+rx27RofffQRqamp9OzZ\nk3fffRdra2tjxiaEMEPV+6dFPOIHbU0by80Crh1nOWvJSQ1AUWbT1Bapm8upBLeqTtT27NFscCzH\nSQlhXHX2qEVHR/PLL7/Qo0cPfvjhB/76178aMy4hhBnKz4dff4VWrWDwYFNHU5NTSR7zeI9x5d9S\nYAY7hjSWuSZqnUqO8H/2k3mu8DV++cXU0QjR8tTZo3by5EmOHdPsSD19+nT69OljtKCEEOapep7S\n4HgNU14AACAASURBVMFmM99dSxXgD2i2uMjIADc3U0bTOEVFcPky2NubYdzXr/PQ1fX0xpd1OxfR\nt28T66oUoomrs0fNysqq1u+FEC1X9bDnmDGmjaNWvr5UocKHbDLPVJo6mkZJTwcHivH3x/yGbPv0\n4Zpre/zI4uSmZFNHI0SLU2eidvToURwdHbVfx44d037v5ORkzBiFEGagrAwytiQRQTxjR5vhgZq2\ntly298IKNQVHs00dTaOkn6niEm7E/+4C16+bOpyaLCywfkhzLEHQya2kpJg4HiFamDoTNbVaTXFx\nsfarsrJS+/2VK1eMGaMQwgzExsKskiXE05+ucR+ZOpxalbbVLP0sPZFu2kAa6cLRPGyoQLG20UwA\nNDNWDz8IwANs5csvTRyMEC2Mzu05hBAC4KsvyhjLD5onUVGmDaYO6RNfIJrP+PVqF1OH0ihXT2j2\nFClp62/aQOoybBiVrR0I4QSxX1xAUUwdkBAthyRqQgidiouhZNuPOFJCeUgv81uaeIPVoxNYRzTJ\n+e1NHUqjVJ1JB6DSx9+kcdTJ1hZid9LZ9QL/TWmHHPcshPFIoiaE0GnLFhhbsRUAm4cfMHE0dWuq\nm95a56YDYNXJDM74rIPV4AGMeUQzP1mGP4UwHoMmarGxsXTt2pWgoCCWLVtWa5k5c+YQFBREr169\nSEpK0r7+9NNP4+HhQY8ePWqULygoYOTIkXTu3JlRo0ZRVFRkyFsQQgCb11cwkW81Tx56yLTB1MPd\nXdP5U1Cg6QVsKioLiinHGoeenUwdSr0mTdI8fvklMvwphJEYLFFTq9XMmjWL2NhYTpw4waZNmzh5\n8mSNMjExMaSmppKSksKaNWuYOXOm9r2pU6cSGxt7W71Lly5l5MiRnD59muHDh7N06VJD3YIQAjh/\nHuL2VPKWxSuUPTQJevUydUh1srCA6qOIm8rqxKIimFe2hLatS3GYPdXU4dRr8GBo3x7OnIFD5nMk\nqRDNmsEStcTERAIDA/H398fa2ppJkyaxdevWGmW2bdtG9I2zAiMiIigqKuLcuXMA3Hvvvbi6ut5W\n783XREdHs2XLFkPdghAC+PpruFrVmpSxz2H7f5vMcKOvmrp10zweP27aOBqqepjWN8AKVStb0waj\ng6UlPPqo5nsZ/hTCOAy2k21OTg6+vn8cLuzj40NCQoLOMjk5OXh6etZZb35+Ph4eHgB4eHiQn59f\na7lFixZpv4+MjCQyMvIO7kIIsXGj5vGJJ0wbR0NNv7ycaezm9I9L4P/bu/O4qMr28eOfYdiRwBUX\nUNzFEtRwK038qpnZgy2aS6m5lFlplk9lm6mlj61atlhPmlvl1lP2yyXMRE1Dctc0VzBAcUNF9mHm\n/v1xB7I7FnAGuN6v13nNnMM5c645wJlr7nVYO6PDua6ckr/GjY2Nw15D7ksn4f0f2PZlB2zvBOIk\nLZ2FuGGRkZFERkbatW+ZJWomO791qwINHew9Lmff4vbPm6gJIf6emBjYvh08PSE83Oho7BOUtpOG\n/MjhPUMAx0/U/pqpj1tuMTYOe3X8+hlW8inTz73K5s3T6dHD6IiEqHgKFiBNmzat2H3L7LtQgwYN\niIuLy12Pi4vD39+/xH3i4+Np0KBBia/r5+eXWz165swZ6tSpU4pRCyHyWrpUP/bvD9WqGRuLvTyC\nWwDgGnvU4EjskzPURYF+Uw7LNHAgAA/xJR9/JD0KhChrZZaohYaGcuzYMWJjY8nKymL58uWEF/hK\nHh4ezuLFiwGIiorC19c3t1qzOOHh4SxatAiARYsWce+995bNGxCiisvKgkUfpeBCFqNGGR2N/ap3\n0YPd+l05QmqqwcHY4fyeeOpyhlturiBJT1gY1rr1acpJLn67hT//NDogISq3MkvUnJ2d+fDDD+nT\npw+tW7dm0KBBBAUF8emnn/Lpp58CcPfdd9OkSROaNWvG2LFj+fjjj3OPHzJkCLfddhtHjx4lICCA\nL774AoDJkyezYcMGWrRowc8//8zkyZPL6i0IUaUtXw4jzr7JaeeG9ExaaXQ4dnMO0iVqLThKgY7m\nDic9HUbFvMoZ6tN6+3+NDsc+ZjPmR0cD8JjtEz76yOB4hKjkTKpgI7FKwGQyFWr7JoSwn1LQqb2F\n1XsbUo9EiIyE7t2NDss+ycng40M67qxYkMqIkY7b2n33bsi8tQtdiIJNm6CidHqKi0MFBmKxmWnj\nE8ee0354ehodlBAVV0l5i+PewYQQhtmyBRruXU09ErEFtdYDaFUUN93EygdX0o2tHDpkdDAlO3hA\nEcRfxX6tWhkbzI0ICMD02mu80mw5J67UzG3LKIQofZKoCSEKmT0bHmceAE7jHnf4sdMKcnpwALsI\n5eAhx77FxUYl4ssV0t194Trtcx3OlCm0nXYfVpz54AOZqUCIsuLYdzEhRLk7fhyOrP6DXmxEeXjC\nsGFGh3TDcoa6cPRBb1N3/QFAWqOgCpcMAwwYoGcq+P132LjR6GiEqJwkURNC5PPBB+BKJofrhmEa\nPgx8fY0O6YY1bQqurnrU/5QUo6MpXvzJLA5wC+b2jjstV0lcXeGJJ/Tz9983NhYhKivpTCCEyHX5\nMvj7Q2oq7N0LIa0t4OJidFh/S0gI7N8PO3ZAx45GR1PYpUtQowZ4eOgJ5M1moyP6e86dg4AAsFj0\n30xwsNERCVHxSGcCIYRdZs3SSVqvXn/NvV5BkzS4NudnzoCyjianWrZ164qbpAHUqQPPjLzMQLWc\nF180OhohKh9J1IQQAPz5J8yZo5/PmGFsLKVh6p7+JFCfE7uvGB1KkSra1FHFysxk5jctWc5gTq09\nyObNRgckROUiiZoQAoApUyAzEwYNcsyqwhtVJyuO+pwh/TfHLFKraFNHFcvNDadBelqpl5nBCy9I\nD1AhSpMkakII9u2DHxZdJMD5DDNnGh1N6XAK0Y2l3I4eMDiSouUkahW+RA3g+edRbm4MYRmZO/bw\n7bdGByRE5SGJmhCC55+HKUzjOE1p8uuXRodTKqp10UVV/pcPcMXBaj+VAtOe3fwfGwkOuGR0OP9c\nw4aYnnwSgJm8xIsvQna2wTEJUUlIoiZEFRcRAUcjYnicebhYMypNtz2nEJ2oBbPf4WYoOHMGhl/9\nkI30ou7GypEY8+KLKG9v+rIel6MHmT/f6ICEqBwkUROiCrNY4LnnYAYv44oF08MPV4JGU3/5K+Fs\nxR/8ftCxGk0dPAht0FWyppDKkRhTqxamefPYMCOa37mFV17RQ3cIIf4ZSdSEqML+8x+osX8TQ/ka\n5e4O06cbHVLp8fNj4ZO/EUAcvx9yrFH/f99ryU3UKksJJgBDh9LrxQ707AkXLsBTTxkdkBAVnyRq\nQlRRe/fC69MV7/EsAKaXX4bAQGODKk0mEzX7hJKJu8NNJZW07TDuZHKldtMKOfNDSUwm+Pxz8PKC\nlSth1SqjIxKiYpNETYgqKCsLRoyAbKuJ//fQchg9WteBVjI5hVU7d4LNZmwseTnv2wVA1i3tDY6k\nbAQGwltv6edPPKFL14QQf48kakJUQW+8oadXatoUJn3aQheBuLkZHVapa9hQT2906ZLjzFBgsUD0\nmQCWMQjPe/sYHU6ZefxxCAuDzPNXGD/e6GiEqLgkUROiitm1C2bO1FVUX3yhq6gqK5MJunfXzx1l\nxPzffoO1Wb2Y1moZXhNGGx1OmXHKSOO7eo+z3xTCmmXJrFxpdERCVEySqAlRhVy8CAMHgtUKEydC\nt25GR1T2wsKgNufYsTHF6FAA2LRJP/boYWwcZc7FBZ+jO2mkTvEezzJqFA7XVlCIikASNSGqCIsF\nHnwQMmJOc2t7VSnm87THgB/HcA4/PH/+wSGmNqpKiRqLF6Pc3BjDfHqn/I9779XV0EII+0miJkQV\nMWkS7P75ElFOt7O11r14WJKNDqlc3NQmEICbr/7K4cPGxpKZCdu26edhYYaGUj5at8b05psALHIa\nCcePMXSoLtEVQthHEjUhqoDPP4cP59pYahpOQ1ssHkmnK2XngaKYbusCQBd+Nbyd2o4dkJGh5/es\nXdvYWMrNhAnwwAN425J5z3Uy69fDq68aHZQQFYckakJUchs36iESpjKVfuoHqF5dD3BVRRI1OnRA\nmUy0Yw/bN6YbGsqmTTCbiUyp+SGkGxtLuTGZYMECePppblr1BWazHmj500+NDkyIisGklCO02ihd\nJpOJSvi2hLhhmzdD374wIH0xixkBTk7www96YxWS0TIY96MHuLfGFr690A2TQRMV3H9bIv/7tR7Z\n7l44X70Mzs7GBGKgefNg3Dj9fOFCPZ6fEFVdSXlLmZaorV+/nlatWtG8eXPe/KudQkETJkygefPm\nhISEsGfPnuseO3XqVPz9/WnXrh3t2rVj/fr1ZfkWhKiwfvkF+vWD9HTFC/WX6I1z51a5JA3Arfcd\nHDO3xJKUzLFjxsSQng7u0VsAUJ1vq5JJGujx1d5+Wz8fNQqWLTM2HiEcnioj2dnZqmnTpiomJkZl\nZWWpkJAQdejQoXz7rFmzRvXt21cppVRUVJTq1KnTdY+dOnWqevfdd0s8dxm+LSEqhF9/VapaNaVA\nqWHDlMpOzVDqq6+MDss4FosaOFBfj88+MyaEjRuV+ohxOogZM4wJwoFMn64vhdms1DffGB2NEMYq\nKW8psxK16OhomjVrRmBgIC4uLgwePJjVq1fn2+f7779nxF/l3p06deLy5cskJiZe91gl1ZpCFGvD\nBrjzTkhJgcGD9aC2Zk83GDLE6NCM4+xs+MC3P/8MYUTqlUo/Nsf1vfJcJvtbDeQu6w88+CB89pnR\nEQnhmMqs7D0hIYGAgIDcdX9/f3bs2HHdfRISEjh9+nSJx86dO5fFixcTGhrKu+++i28RkxpPnTo1\n93lYWBhhVaIvvKjq/vtf3f7HaoVBg2DJEjCbjY7KMeRN1JSi3Nup7f3xLG9wmGw3T5xDQ8v35A7I\ntHQJbf5YxXfm7xloXcbYsfcREwMzZuimlEJUZpGRkURGRtq1b5klaiY774I3Wjo2btw4pkyZAsCr\nr77KpEmTmD9/fqH98iZqQlR2Nhu89BK8+SZ0Iorez7Rh2jte8oGXR+vWULMmxMdDTAw0aVJ+505J\ngZ/3VKePKYLvZp3E2cWl/E7uqEaPhoMHcX7/fb5xGsgotYBZs4YTG6s7GVSVTsmiaipYgDRt2rRi\n9y2z23iDBg2Ii4vLXY+Li8Pf37/EfeLj4/H39y/x2Dp16mAymTCZTIwZM4bo6OiyegtCVAhJSXD/\n/TpJe9j0Jducu/P6wftwsmQaHZpDcXKCO+7Qz8u7+nPbNki3unIptDceE8eW78kdlckEs2fDlCk4\n2awsVCN43fV1li1TdO0KJ04YHaAQjqHMErXQ0FCOHTtGbGwsWVlZLF++nPDw8Hz7hIeHs3jxYgCi\noqLw9fXFz8+vxGPPnDmTe/y3335LmzZtyuotCOHwtm6Ftm3h+9U2/uP2GkvUw5izs6BVKz2Fj8hn\nUPPdvMLr/Lz8fLmet8pMG3WjTCaYNg3mzAGTifFdfqNxQxs7d0L79rBihdEBCuEAyrIXw9q1a1WL\nFi1U06ZN1cyZM5VSSs2bN0/Nmzcvd58nn3xSNW3aVAUHB6tdu3aVeKxSSg0bNky1adNGBQcHq/79\n+6vExMRC5y3jtyWE4bKzlZo6VSknJ6VqcEFt9e2nu9A5OSn1wQdGh+ewMrr3VgrUcKcl6ty58jmn\nzaZUixb61xMRUT7nrJA2bFAqOVklJSl1//36eoFSjz2mVEqK0cEJUbZKyltkwFshKphdu2DsWP1o\nMsH6bm9w55ZX9YwDX35ZJcdJs9ucOfDMM6xgIGfmrODpp8v+lFu36irX+vXh1KkqO3zaDVEKPvkE\nnn1Wz4/aqJFelz9tUVkZNuCtEKL0XL0KzzwDHTvqJC0gACIi4M6fXoCnnoK9e+WT7Hruuw+Au1nL\n1wvKZwqnzz+H2pzjkUckSbOXyaSnPduxA+66OY4Lp1K4+249wszZs0ZHJ0T5kkRNCAdntcKiRbrX\n4pw5etuzz8KhQ9CrF7ot2ty50LChoXFWCI0aYWt/K9VIpc7+DezbV7anu3IF9i8/TCJ1eWlH/7I9\nWSUUEpTFWrd7ia8Zwj2uESxbBi1b6o4zVWWqVCEkURPCQSkFa9bozgKPPAKX468yPOg3du6Ed9+F\natWMjrBichrwAAAPsoJFi8r2XMuWQXjmCpxQePlXL9uTVUbnzmHKzsb34kn+X1YfNtUbivuVRCZP\nhhYt9DAeVqvRQQpRtiRRE8LBKAVr10K3bnDPPXDooJV/11jAWZ8WLEz6F+2aJhsdYsU2eDAJj01j\nCtNZuhQslrI71fzPFcP4a57VoUPL7kSVlb8/7NwJs2aBhwdhZ74mzqsVr/nPJz4eRo6ENm1g8eKy\n/T0KYSRJ1IRwENnZ8NVXEBKiJ1Pftk0xsNo6TjfowNtJo/G8koipcWO4cMHoUCu2xo2pP28KHq2b\ncP48rFtXNqfZvx/MO6NoxglU3XrQs2fZnKiyc3GBF16A33+Hvn1xSb3ClOlmFi/Wtf2HD8OIEdC8\nOXz0EaSmGh2wEKVLEjUhDJaQoIeSCgyEhx6CAwd0D8Hdt49nRcrd+CXs0T0HvvoKtm8v3yH1KymT\nSX+4g64+Kwvz58NIvtDne2iozOX1TzVurNsC/PwzTsOHMWwYHDum57Jt2VL3qH3qKWjQACZOhKNH\njQ5YiNIhw3MIYYCsLPjxR50krF59rZ1Ny5bw73/DsGHg9muknrDz+ef1BJ6enkaGXOmcPq3zX7NZ\nJ8u1a5fea2dm6mR7XNIbvFZjLi7bNutBiEWZsFrh++XpHH1lMW/EDCUFb0AXYj7yiO7s6+VlbIxC\nlKSkvEUSNSHKic2mhxv48kvdyPziRb3dx5xC7/uqMW6cHrk+3zS56eng4WFIvFXBv/4FP/wATz4J\nH35Yeq+7cKFuP9W2LeyJtsgsEeVh7lyYMAGrlzebGw1n0okn2ZsZBOgk7YEHdIl1jx7y6xCORxI1\nIQySnQ1btsD//gfffqtLcQBcyeTxgDU8WW0RzU5txOnkCfDzMzbYKujAfsWkdhvZq9oSsbsWbdv+\n89e8dEkXnp07pxO2nCpWUcbWrYP//EePMPyXhBY9mOY2g/8e6JK7zdcXwsN14tarlxRUC8cgiZoQ\n5ejMGVi/Xn9uRETosbRy3O23kxfrLqRz7Nc4X0nSG81mWLoUBg82JuCqbNIkeO89pvMqEbdPZ+vW\nAiWaf8PYsfDZZ3o2gsjIf/564gbt3w8ff6z/p1JTYdMmjvuHsXQprFypxx/M4eYGXbvCnXdC7966\nI4+TtNwWBpBETYgyFBenS81ylj/+yP/zli3h/vv1cuuSpzF98IH+QUiILm4ZOlRK04yybRt07coV\nkw9N1AlmL6rJ8OF//+W2b4fbb9dVa/v2QVBQ6YUqbtCVK/DNN7qRWp7s648/dAl3/KKfWHK0Mylc\nG5Cwdm2dsPXurdu3BQQYELeokiRRE6KUWK1w5Ij+QM5JzE6dyr9PK/dY+tx6geZDQunbt0Anzeho\nWL5c9xYojXo28c/16QMREczhaWb5zeHIEfDxufGXsVigfXvIOHiMIS83ZfobUjTjsE6fhgYNUG5u\nnGndk0jXPnwe24tNZ4OAa0Wg/v7QubNeunTRv193d+PCFpWXJGpC/A3Z2bqaZPduvezapafTTEvL\nv19z70Qeab6Nu9w2ERQfgUfcMejUCaKijAlc3Jj9+1Ft22LFTFu1h54TbuH992/8Zd58E96ZfJ4j\nTkH4BjfEaUME1KpV+vGKf27vXj2ZaFSUHmH6L2kduvP5w5FERMAvv+RvtgC6pLRtW524tWsHwcF6\najfp7yP+KUnUhCiBzQaxsXrgzEOH9OPBg3o8s4yMwvs3bAgdOug2SL2axdK6X+P8O/j46FKar7+W\nBi8VxbhxMG8eSxjGcBYzYwa8+KL97cvWrIEBA+DzjId4iK90vdmGDdJAzdElJurGpBs3wk8/6V/i\nX91/bTZdTRoVBUc3nOLUb+f45kRbLOTvMurkpKezCg6+ttx8MzRqJEPnCftJoiaqPKXg7Fk4eRJO\nnNCPR4/qxOzIkaIneDaTTW//w9xTbzcdXfdy6eV3aBdqzj/ellL6Lh0YqOd86t1bZ3HOzuX11kRp\nSE6GTz5hvvdEHn3KDaXgmWfgnXeun2t/8okeaPUh22IWM0IXrxw8KAMTVzRK6eLyogZce/11mDIF\n5e5OcmAwJ25qx25bO1Zd7sVPMU2LnG/U1RWaNtW3h4KLn5/k8CI/SdREpZedrXtbxsfrxv05jzEx\n1xKzglWWedWvrxt+t24NYw5OJPD0NrxPHcSUt0jt0KGiW4crJXfdSmTFCnj4Yd3mbPhw+Pzzosfd\nstn0zEbvvAMdiOYXlx64WtL0AaNHl3/gouzMmQPz5ulvdXl98AEZj47n8GHd2TRnSfz9IkfOeGPB\ntciXc3fX3+2KWgICoE4d+a5X1UiiJiosqxWSknRp2Llz+jE+vvCSmKg/OItjwkaQ7xlur3uStjed\npKVbDBfDR9Gwa0OCggo0Hu/WTTdQAV0q0r69XkaM0BmdqPQiIvRo9mlp+oMzPBz699dtk3bt0jVl\na9fqtovOzrAv7Gla//QBjBqlEzVJ3CunpCTdvm3PHr08+6y+NxQ0YABq9Woy/ZtysVYrYj1acSi7\nJd9m9GXHqbokJZV8GicnXepWv37+pUGD/Os1a0rrispCEjXhEJTSwxpduqSXpKRrjzlJ2Llz+Z9f\nuFByApbDmWzq17XhF+CKv7/+cG3QQE8P2GPpKGr++BWmzMz8B61apUe9LGjLFv1Be8stUL166bx5\nUeHs2AHjHzjNzoS6qGKmRfb11WNz9fo/m550cvhwGfZeQN++eo64gp9DW7ZAt24kJ+ve4rGxevH5\ncTkXzin2JQey/0IDDl6sSzbX/zsym3V/lVq19NAitWvnf553vVYt/ffq4SHfIxyRJGqiVNhsOtFK\nToarV/WS8zw5OX/yVdxzi+XGz1ujhq4KyFnuyvyOdpd+pnZWAr4p8XgkJWA+dwbTwoV62IuCHnsM\n/vtffbdq0uTaMmCADJEhipeaimrfnjSnavzQ8t+8f+Qufv2jOiG3WOnRy0zPntC9O3h7Gx2ocEjp\n6XrW+MOHda+Eo0fh3Xehbt3C+7Zqla9aVZlM2Gr7cfCDTZx0bcXp0+QuCQngFHuCmLNenEiuhZUb\nqyN1ddUJW/Xq+tGe5z4++u+8WjXdhE9K8UqfJGpVTFaWrrLJWVJT868XtS1nvagELGfb1av/PDYP\nD/3PX6PGtcdOmVton7qFOuoc1bPP4Z1+Ds+Uc2SOfw6vJ0YULqCYMEHP61fQnDnw9NOFt1+4oIcg\nl09UcSP27oV+/a7N+wUod3dMd92l5wMTorS88opO6P78U2diiYm6NC4xsejBsBs1gj//1Amdbw0s\n1euQXq02kU+sJD6rDufPw/nz+tZ3/jx4xx/m9BUvYq9UJ8lSjbxjxf0d1arpJSd58/Yu/rmXl56m\nK2fx8Ch+3c2t6pb2SaJWDpTS7akslmtLZqZeMjIKP7/e443sm56eP+nKzi6b9+iElbqeV6lX7Sp+\nHslk+9TAUrNe7j9mTvIVGv8drf/4hmpZl/DISMIt7RLOyUnYnnsBl+efLfzCL72k5+graOpUeO21\nwts3bdLtQ/z9df2mvz/Uq6e/KgpRmtLTYcECWLZMN05LT9elH4cPGx2ZqMyys3WSVr9+4eIrpfQg\nbgkJcPFi/urV5OSiv5BWrw6XL+vDnZ2x+VQn27s6+z7dwYVsXy5f1j++dEk/ttj5FVfSXTmfeRPn\n0r05m+bN2fSbOJrmT0pa2RWnmUz5E7eCz93ddTLn5pb/ecH1kn5W1Lqrq15cXPTi7Fz+CWOVTNRm\nz1b5kqaCS1ZW8T/7uz93DApPp0xqeaXh6uGM8r4p37cXLy9omXWAoNTfuMkpFS9TGp6mNDxVGmeD\nwkjq0o+bbrr2reimm6Del2/jM3sqpoLdJl99FaZPLxzCG2/onxX0wgswa1bh7Rs26OQrb/1m7dp6\nwDJpIyYchc2mvx25ukqXPOEYrFadrJ07p4vOwsIKZxhK6SYeFy/qTCzvfTwrq3CbSqV0JpOVVfh8\naWlYXT3y1b6kpEDjx+/EgisZLt6kmb1JNXuTihdr2k8hJcO5UI1Os3PbuZLlzqVML5IyPUlK9yQp\nw5MrWe7809K+0pKTtLm45E/i8j6/3vr19s1JCp2dYeLE4hO1Snu3eeaZ8j6jwmxSOLs65fuluLtD\nPefzNHSKo5pzJl7OmXg6Z+JlzuSSbyDn/ILzZfju7tA8cStBJ37AjUxcbRm4qkxcVCbnQ/uysV4A\nHTqE5TvGb9VH1PxoGk4Z+r/AZFNwFRj3vB4uvaC31umkqaCOwIh+hbdXd9L/XSZT/gzO17foS/Gv\nf+nW/AXrOGvUKHr/nMn1SllkZCRhYWGl/rqVjVwnOzk5ERkdLdfKDvI3ZZ9/fJ3M5mtfbotjMumJ\nZ3NkZV0rOiuq44vVCkOGFG77kpYG7u6Y83wMALr0b+eGIk/d/eepFOqHoxQ43V7k/tYsK+mZJtLS\n8tQUpSqajerG5rRUOvnUJdvsjsXJjSwnNzYMXUhmlilfTVRmhqLTro/IsLmRptxJt7mRanUnzerG\ndq87C+2fka7wyTpPmsWFVIsraRYXMmwuWCwmhymAKdNEbf369UycOBGr1cqYMWN4oYjkYMKECaxb\ntw5PT08WLlxIu3btSjw2KSmJQYMGcerUKQIDA1mxYgW+RSQM/+s2GxcsuJgsuGLBGQuJgV34M/ie\nQllto4NraLZ5PmabBbPNgtNfj0n/N4CLQycUyoB9V3yGz9svg8WCKVsXp5ksFpjwtG4nVdB7S2DS\npMLbJ06E2bMLb5+9E754q9DmwNDarMk4Qc+eYfl/8HM2XDx/bd3dXRefFdf7LCREDyNQsKFAbYf/\nlQAAC51JREFUp05F7//EEzB2rC6Os6c8OCRELwaTDwv7yHWyn1wr+8h1so8h18nVVbd7K6rtG+ji\nnYUL7X89k0nXiORt4JySorOgonodWCxw222FG07bbJhdnKjmotu35cqywIFt7AbC876Oiwstti8q\n/PpZFnAbX3i72Vx0uyBLNrgWvhbKzY3UCxmFa9HSs2n4r2CU2QWrsys2sws2swvZzu7smPZjof2z\nM610/nI8FpMLFlzIxgWLyYVM3NjY6WWys/V+RTW7zlFmiZrVauWpp57ip59+okGDBnTo0IHw8HCC\n8gwYunbtWo4fP86xY8fYsWMH48aNIyoqqsRjZ82aRe/evXn++ed58803mTVrFrOKqE67b2sRbaGC\nn4JJ9xTePvckRBduHOx7RwhNOhTx5jwz4eKFwtuLaxxWr54ues5bKe7mpkdXLUrXrrrNVsHK9KAg\n+P77wvuPHg1Dh+qEy939+vOW9OmjF3vJRHZCCCGKYjbrKld7ubrCtm039vpbt+qe+4MGXWucXVJj\n7CeeKNyguzgWix67JG+GlZ2NyWTKnzDmSMuCk0W0UXV3p18RFVKkZ8HjnxTe7ubGnZEv564akqhF\nR0fTrFkzAgMDARg8eDCrV6/Ol6h9//33jBgxAoBOnTpx+fJlEhMTiYmJKfbY77//ns2bNwMwYsQI\nwsLCikzUePrpwpXBHYrKutBj3qxaVXh/f/+i9x81Sv/BFNy/uD7LQ4boxV4dOhQfa1GJWk4XHCGE\nEKIyMZt14cVPP8Hdd19/f1dX+Ogj+1/f01O378tLqeITQXd3PUVcwaKz4gb8dHbW88cWSARvqLeC\nKiMrV65UY8aMyV1fsmSJeuqpp/Ltc88996ht27blrvfs2VPt3LlTrVq1qthjfX19c7fbbLZ86zkA\nWWSRRRZZZJFFlgqzFKfMStRMdmaLyo5Op0qpIl/PZDIVud2e1xRCCCGEcHRlNiBKgwYNiIuLy12P\ni4vDv0BVYsF94uPj8ff3L3J7gwYNAPDz8yMxMRGAM2fOUKek3i5CCCGEEBVYmSVqoaGhHDt2jNjY\nWLKysli+fDnh4fn6bBAeHs7ixYsBiIqKwtfXFz8/vxKPDQ8PZ9Ei3dNj0aJF3HvvvWX1FoQQQggh\nDFVmVZ/Ozs58+OGH9OnTB6vVyujRowkKCuLTTz8FYOzYsdx9992sXbuWZs2a4eXlxRdffFHisQCT\nJ0/mwQcfZP78+bnDcwghhBBCVEo32EdAlKN169apli1bqmbNmqlZs2YVu190dLQym83qm2++Kcfo\nHIc912nTpk2qbdu26uabb1bdu3cv3wAdxPWu0/nz51WfPn1USEiIuvnmm9UXX3xR/kE6gJEjR6o6\ndeqoW265pdh9xo8fr5o1a6aCg4PV7t27yzE6x3K9a7V06VIVHBys2rRpo2677Ta1b9++co7QMdjz\nN6WU3MvtuU5V8V4uiZqDys7OVk2bNlUxMTEqKytLhYSEqEOHDhW5X48ePVS/fv3UqlWrDIjUWPZc\np0uXLqnWrVuruLg4pZROSKoae67Ta6+9piZPnqyU0teoRo0aymKxGBGuobZs2aJ2795d7IfFmjVr\nVN++fZVSSkVFRalOnTqVZ3gO5XrXavv27ery5ctKKf1Foapeq+tdJ6XkXq7U9a9TVb2Xl93squIf\nyTsOnYuLS+5YcgXNnTuXAQMGULt2bQOiNJ491+mrr77igQceyO3MUqtWLSNCNZQ916levXokJycD\nkJycTM2aNXGugnNaduvWjeolzDFb1PiPZ8+eLa/wHMr1rlWXLl3w8fEB9LWKj48vr9AcyvWuE8i9\nHK5/narqvVwSNQeVkJBAQEBA7rq/vz8JCQmF9lm9ejXjxo0D7B8SpTKx5zodO3aMpKQkevToQWho\nKEuWLCnvMA1nz3V69NFH+f3336lfvz4hISG8//775R1mhVDUtayqCciNmD9/PnfbM2BpFST3cvtU\n1Xt51fu6XEHY8486ceJEZs2ahclkQulq7HKIzLHYc50sFgu7d+9m48aNpKWl0aVLFzp37kzz5s3L\nIULHYM91mjlzJm3btiUyMpITJ07Qu3dv9u3bh3fu7MsiR8H/NflgLdmmTZtYsGAB225k6qAqRO7l\n9qmq93JJ1ByUPePQ7dq1i8GDBwNw4cIF1q1bh4uLS6FhUCoze65TQEAAtWrVwsPDAw8PD+644w72\n7dtX6f+587LnOm3fvp2XX9ZzzzVt2pTGjRtz5MgRQkNDyzVWR1fSOI+isP379/Poo4+yfv3661b/\nVVVyL7dPVb2XS9Wng7JnHLqTJ08SExNDTEwMAwYM4JNPPqly/9j2XKf+/fvzyy+/YLVaSUtLY8eO\nHbRu3dqgiI1hz3Vq1aoVP/30EwBnz57lyJEjNGnSxIhwHVpx4z+Kwv7880/uv/9+li5dSrNmzYwO\nx2HJvdw+VfVeLiVqDsqeceiEfdepVatW3HXXXQQHB+Pk5MSjjz5aJf6587LnOr300kuMHDmSkJAQ\nbDYbb731FjVq1DA48vI3ZMgQNm/ezIULFwgICGDatGlYLBag5PEfq6LrXavp06dz6dKl3LZXLi4u\nREdHGxmyIa53nYR2vetUVe/lJiWV4UIIIYQQDkmqPoUQQgghHJQkakIIIYQQDkoSNSGEEEIIByWJ\nmhBCCCGEg5Jen0KISsFsNhMcHJy7vnr1aho2bGhgREII8c9Jr08hRKXg7e3N1atXi/xZzm1OZhAQ\nQlQ0UvUphKiUYmNjadmyJSNGjKBNmzbExcXx9ttv07FjR0JCQpg6dWruvjNmzKBly5Z069aNoUOH\n8u677wIQFhbGrl27AD1ifOPGjQGwWq0899xzua/12WefARAZGUlYWBgDBw4kKCiIhx9+OPccv/32\nG7fffjtt27alc+fOpKSk0L17d/bt25e7T9euXTlw4EBZXxohRAUiVZ9CiEohPT2ddu3aAdCkSRPe\ne+89jh8/zpIlS+jYsSMREREcP36c6OhobDYb/fv3Z+vWrXh6erJ8+XL27duHxWKhffv2udNmmUym\nIkvh5s+fj6+vL9HR0WRmZtK1a1fuvPNOAPbu3cuhQ4eoV68et99+O9u3byc0NJTBgwezYsUKbr31\nVlJSUvDw8GD06NEsXLiQ2bNnc/ToUTIzM2nTpk35XTQhhMOTRE0IUSl4eHiwZ8+e3PXY2FgaNWpE\nx44dAYiIiCAiIiI3mUtNTeXYsWNcvXqV+++/H3d3d9zd3e2auiciIoIDBw6watUqAJKTkzl+/Dgu\nLi507NiR+vXrA9C2bVtiYmLw9vamXr163HrrrQBUq1YNgAEDBvD666/z9ttvs2DBAkaOHFl6F0QI\nUSlIoiaEqLS8vLzyrb/44os89thj+ba9//775G2qm/e5s7MzNpsNgIyMjHzHffjhh/Tu3TvftsjI\nSNzc3HLXzWYz2dnZxbaN8/T0pHfv3nz33XesXLmS3bt338C7E0JUBdJGTQhRJfTp04cFCxaQmpoK\nQEJCAufPn+eOO+7gu+++IyMjg6tXr/LDDz/kHhMYGMjOnTsBckvPcl7r448/Jjs7G4CjR4+SlpZW\n5HlNJhMtW7bkzJkzua919epVrFYrAGPGjGHChAl07NgRHx+f0n/jQogKTUrUhBCVQlGlVnm39e7d\nm8OHD9OlSxdA9xJdunQp7dq1Y9CgQYSEhFCnTh06dOiQW6r273//mwcffJDPPvuMfv365b7emDFj\niI2NpX379iilqFOnDt9++22xbdpcXFxYvnw548ePJz09HU9PTzZs2ICXlxft27fHx8dHqj2FEEWS\n4TmEECKPadOmUa1aNSZNmlQu5zt9+jQ9evTgyJEj5XI+IUTFIlWfQghRQHmNt7Z48WI6d+7MzJkz\ny+V8QoiKR0rUhBBCCCEclJSoCSGEEEI4KEnUhBBCCCEclCRqQgghhBAOShI1IYQQQggHJYmaEEII\nIYSDkkRNCCGEEMJB/X/T37+rt20/ygAAAABJRU5ErkJggg==\n" } ], "prompt_number": 6 }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Software versions" ] }, { "cell_type": "code", "collapsed": false, "input": [ "from qutip.ipynbtools import version_table\n", "version_table()" ], "language": "python", "metadata": {}, "outputs": [ { "html": [ "
SoftwareVersion
Cython0.16
SciPy0.10.1
QuTiP2.2.0.dev-d83b39f
Python2.7.3 (default, Sep 26 2012, 21:51:14) \n", "[GCC 4.7.2]
IPython0.13
OSposix [linux2]
Numpy1.7.0.dev-3f45eaa
matplotlib1.3.x
Thu Feb 14 12:41:03 2013
" ], "output_type": "pyout", "prompt_number": 7, "text": [ "" ] } ], "prompt_number": 7 } ], "metadata": {} } ] }