{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# QuTiP Example: Parallel Steady State Jaynes-Cummings Calculation" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "J.R. Johansson and P.D. Nation\n", "\n", "For more information about QuTiP see [http://qutip.org](http://qutip.org)" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "%matplotlib inline" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "import matplotlib.pyplot as plt" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "import numpy as np" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "from qutip import *" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## System Setup" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [], "source": [ "kappa = 2\n", "gamma = 0.2\n", "g = 1\n", "wc = 0\n", "w0 = 0\n", "N = 5\n", "E = 0.5\n", "nloop = 101\n", "wlist = np.linspace(-5, 5, nloop)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Define Function to Run in Parallel" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [], "source": [ "def probss(wl, E, kappa, gamma, g, wc, w0, N):\n", " # construct composite operators\n", " ida = qeye(N)\n", " idatom = qeye(2)\n", " a = tensor(destroy(N), idatom)\n", " sm = tensor(ida, sigmam())\n", " \n", " # Hamiltonian\n", " H = (w0 - wl) * sm.dag() * sm + (wc - wl) * a.dag() * a + \\\n", " 1j * g * (a.dag() * sm - sm.dag() * a) + E * (a.dag() + a)\n", "\n", " # Collapse operators\n", " C1 = np.sqrt(2 * kappa) * a\n", " C2 = np.sqrt(gamma) * sm\n", " C1dC1 = C1.dag() * C1\n", " C2dC2 = C2.dag() * C2\n", "\n", " # find steady state\n", " rhoss = steadystate(H, [C1, C2])\n", "\n", " # calculate expectation values\n", " count1 = expect(C1dC1, rhoss)\n", " count2 = expect(C2dC2, rhoss)\n", " infield = expect(a, rhoss)\n", " return count1, count2, infield" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Run Simulation" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [], "source": [ "count1, count2, infield = parfor(probss, wlist, E=E, kappa=kappa, gamma=gamma, g=g, wc=wc, w0=w0, N=N)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Plot Results" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYwAAAEKCAYAAAAB0GKPAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xd8VfX9+PHX+96bTQgjCTsJG9lLULCKG0fFveuodWuH\nba3dX9va/r7Vtmq1teLXWuveIg5cCA4Qwt7DJEAgkEAgk4x77+f3xzk3xBiSG7j3nntv3s/H4z5y\n77lnvG8g930+W4wxKKWUUu1xOR2AUkqp2KAJQymlVFA0YSillAqKJgyllFJB0YShlFIqKJowlFJK\nBUUThlJKqaBowlBKKRUUTRhKKaWC4nE6gFDKzMw0eXl5ToehlFIxY9myZXuNMVnB7BtXCSMvL4/8\n/Hynw1BKqZghItuC3VerpJRSSgVFE4ZSSqmgaMJQSikVFE0YSimlgqIJQymlVFA0YSillAqKJgyl\nlFJBiatxGEo199mWvTy3ZBv9uqVwbF4Pjs3rQfe0RKfDiqh91fUsLSpnSeF+9lTWce20PKYM7OF0\nWCpGacJQccfnNzz04Wb+Pn8rgSXrZ39aiNsl/OmCMVx67ABnA4yQpz4v5N6565t+BwDvrdvNXacP\n49aTBuNyiXPBqZikCUPFlX3V9dz+3HIWF5TjErhlxmA8LuHLwnKWFJbzmzlrmZzXnUFZXZwONazW\n7argvnc2YAwcP6gnUwb2oLrey/99Vsj98zaxpLCchy+fQEZqgtOhqhiiCUPFDWMMP3hhJYsLysns\nksTDV4xn2uDMpvfvenElr63YyY9fXsXLNx+Pxx2fTXj1Xh8/fmkVjT7DNcfn8rtZo5veO2FIJj96\naSULNpdx96ureOzqSYhoSUMFJz7/YlSn9MLSHXy2dS/dUhN4687pX0sWAL89bxS9uyazYvsB/rWw\nwKEow++hD7ewcXcVeT1TueesEV977+QR2bx5+3S6JHmYt24Pc1eXOBSlikWaMFRc2HngIPe9vQGA\ne88bRZ+MlG/sk5GSwJ8vHgvAgx9uZv2uyojGGAnLt+/nsQVf4RL4y6XjSE38ZiVCbs80fnH2MQD8\n5s217K2uj3SYKkZpwlAxzxjDz19bQ3W9lzNG9uK8cX0Pu++Jw7K4amoOjT7Do/O3RjDKyHjowy34\nDdz4rUFMyj18b6grpgzghCGZ7K9t5LdvrotghCqWacJQMe/V5TtZuLmMbqkJ/OGC0e3Wyd928hAA\nPtywh+p6byRCjIh91fV8tnUvHpdw80mD29xXRPjThWNIS3Tz9poS3lurVVOqfZowVExr8Pr52web\nAfj1OSPJTk9u95h+3VKYnNudeq+fD9fvCXeIEfPO2t34/IYThmbSI4jxJgN6pHL3TKuN44H3N+P3\nm3aOUJ2dJgwV015fUczOAwcZnJXG+RP6BX3ceeOtaqs5q3aFK7SIe2ul9VnaqpJr6YopOfTrlsLW\n0mreXbs7XKGpOKEJQ8Usr8/Po/O/AuDOU4bi7sBAtLPH9MHtEhZuLmN/TUO4QoyYXQcOsqSonCSP\nizNG9Q76uESPi1tmWNVXf/94i5YyVJs0YaiYNWfVLraX15LXM5Vzx/bp0LGZXZKYNrgnXr+Jizvr\nuaut0sWpx2TTJaljw6sumdSfXl2T2Li7ig83xE8VnQo9TRgqJvn8hkfsXk63nTzkiAbhBapu5qza\nGdLYnBCoWutIdVRAcoKbm08MlDK2YoyWMlTrNGGomPTOmhIKymro3z2FCzrQdtHcmaN7k+hx8WVh\nObsr6kIcYeQUlFWzdmcl6UkeZgzPPqJzXDElh8wuSazZWcEnm8tCHKGKF5owVEx64rNCAG45aTAJ\nRzjFR9fkBE4enoUxh6p0YtFbq6wusWeM6k1ygvuIzpGS6OamEwcC8H+fFoYsNhVfNGGomLO6+ACr\ndhwgIyWBiyb2P6pznT3GavtYuGVvKEJzxILNpQCcPSb4xu7WXHZsDikJbj7bupevyqpDEZqKM5ow\nVMx5etE2wGqsTUk8sjvqgMDaECu27ccXgz2E6hp9rN1ZiQhMzju6dS4yUhI4f4LVBvJf+3esVHOa\nMFRM2V/TwFt2A+/Vx+Ue9fn6ZKTQr1sKVfVetpRWHfX5Im3tzgoafH6GZaeTkXL0U5UHfqevLium\ntiF+RsGr0NCEoWLKy8t2UO/1c+KwLPIy00Jyzkm53QHIL9ofkvNF0rJtVsyT8rqH5Hyj+mYwKbc7\nVfVe3lgRu+06Kjw0YaiY4fcbnlm8HYBrQlC6CAgkjOXbYi9h5AcSRk5oEgbAd+zf7dOLirSLrfoa\nTRgqZizYUsb28lr6dUvh5BFH1n20NU0ljBhLGMaYpiQ3OUQlDICzxvSmZ1oiG3dXNZVglIIwJwwR\nmSkim0Rkq4jc08r7V4nIahFZIyJfiMi4YI9Vnc+zi62G2KuOy+nQNCDtGdE7ndREN9vLaymtip3x\nGEX7atlX00Bml0RyeqSG7LxJHjeX2eueP7NYG7/VIWFLGCLiBh4FzgJGAleIyMgWuxUCJxljxgC/\nBx7vwLGqEymtrGP+pjI8LuGSSQNCem6P28X4Ad2A2KqWCtz9T8zpHvJlVq+YkgPAu2t3U3GwMaTn\nVrErnCWMKcBWY0yBMaYBeAGY1XwHY8wXxpjAX+hioH+wx6rO5bUVO/H5DaeMyCYrPSnk559sV0vF\nUhXMsjBURwUM6JHKtME9qff6m3qlKRXOhNEP2NHsdbG97XBuAN49wmNVHDPG8FK+9d/h0smhLV0E\nTIzBdoxl28qBQ20woRb4Xb+cv6OdPVVnERWN3iJyMlbC+NkRHHuTiOSLSH5Zmc6BE4+Wb99PQVkN\nWelJzBieFZZrTMjpjog1rqGu0ReWa4RSxcFGNu+pJtHtYnS/jLBcY+bo3qQne1hVXMHG3fG3/rnq\nuHAmjJ1A89vB/va2rxGRscATwCxjzL6OHAtgjHncGDPZGDM5Kys8XybKWS8tLQbgwon9jmhW2mBk\npCQwLDudRp9hzc6KsFwjlJZvt0pCY/pnkOQ5utHuh5Oc4G6a/fbl/OKwXEPFlnAmjKXAUBEZKCKJ\nwOXAnOY7iEgO8BrwHWPM5o4cqzqHmnpv08SAoW7sbikw+C0W2jGautOGqToqIFAt9fqKnTR4/WG9\nlop+YUsYxhgvcAcwD9gAvGSMWScit4jILfZuvwF6Av8QkZUikt/WseGKVUWvd9aUUNPgY1Jud4Zk\ndwnrtQKD31Zsj/6EsXLHAcCqSgunsf0zGN4rnfKaBj7eqIsrdXYdW5qrg4wx7wDvtNj2WLPn3wO+\nF+yxqvN5eZlVFXLp5KOblTYYI/t2BWDznuifqXXjbmveq1F2zOEiIlx67AB+P3c9L+cXM3N0x1Y2\nVPElKhq9lWrNjvJalhSWk5zg4pyxHV9JrqMGZaXhdglF+2qiuuG7vKaBsqp60hLd9OuWEvbrzRrf\nF49L+GRzGWVV9WG/nopemjBU1HptudXPYeao3h1ep/pIJHncDMxMwxjYWhq9pYxNduliaK90XCEc\n8X44mV2s3mk+v2laClZ1TpowVFQyxvDaikDvqPBXRwUM75UOHKryiUab91ixBWKNhMC/wavLtLdU\nZ6YJQ0WlZdv2s21fLb26JjF9SGbErju8t/UlHPhSjkabAgmjd+QSxqnHZNM12cP6kko2lOiYjM5K\nE4aKSq8ut+5kz5/QL6QTDbZnmH3XvimaSxi7I58wkjxuvm2PyXhtuZYyOitNGCrq1DX6mLu6BOCo\n1+zuqGgvYRhjmkoYwyJYJQVw0STr3+KNlbvw+nRMRmekCUNFnQ/W76GqzsuYfhkR/1LM6ZFKcoKL\nkoq6qJyltaSijqo6Lz3SEsnskhjRa08Y0I2BmWmUVdXz6da9Eb22ig6aMFTUCVR5XDgx8vNNul3C\n0GwrSW2JwlLGpmYN3qGe0rw9IsKFE6x/E2387pw0YaioUlpVx8Ite/G4pGkeo0gbFsU9pZxov2ju\nAjuJf7B+D5V10VcCU+GlCUNFlTkrd+HzG2YMz6Znl9CvexGM4b2tKUiisR3DqfaLgP7dUzluUA/q\nvX7eXVPiSAzKOZowVFR5fYU1WM+J6qiAaO4ptamphBHeebXaEhiTERhYqToPTRgqamzcXcm6XZV0\nTfZwyohsx+IY0Tswp1QVxhjH4mjJ5zdssUegO1XCADhrdG+SE1x8WVjOjvJax+JQkacJQ0WN1+07\n1nPH9SU5ITxrPASjV9ckuiZ72F/bSFl19MydtG1fDQ1eP/26pZCenOBYHOnJCZw5qjcAb6zQUkZn\noglDRQWf3zRVR13kYHUUWL2BmsZj7I6eOaU2N7VfOFcdFdBULbViZ1SVwlR4acJQUeHzrXspraon\nt2cqE8O8xkMwDvWUip5pMDbZyWuYQz2kmps+uCfZ6UkU7q1hhb02h4p/mjBUVGgaezGhf8THF7Rm\nRO/oa/jetMdKXiOiIGF43C7O1zEZnY4mDOW4qrpG3lu3G4ALJjhbHRUwxB6891VZ9FRJBaZcH5Ll\nfMKAQz3Z5q4uier1Q1ToaMJQjntnTQl1jX6mDOxBTs9Up8MBYGBmGgBF+6KjF5Dfb9hmxzIwK83h\naCwjendlVN+uVBxs5KMNpU6HoyJAE4Zy3KvLrMbuiyM80WBbenVNIiXBTXlNQ1TMKVVSWUe9109W\nelJEFpMKVmByyFd1BttOQROGctS2fTUsKSonJcHN2WOjZ71oESEvUMrYW+NwNIdiGNgzOkoXAYHl\nWxdsLqO0qs7pcFSYacJQjno1sAzr6Mgsw9oRAzOt6rGifc4njEI7YeRlRkeVXUDPLkmcPCIbn9/w\n5gpdvjXeacJQjvH7TVPvqIsnRU91VECefTdfUBZNCSO6Shjw9WopHZMR3zRhKMcsKSqneP9B+mYk\nc/ygnk6H8w1NVVJRUMKI1iopgFNGZNM9NYGNu6tYtyt6xq2o0NOEoRzzyrLAuhf9cUVwGdZgDYyi\nNozCfdFbwkj0uJg13upi+4qOyYhrmjCUI6rrvbxjT4/t5My0bQlUSRXurXG0qsXr8zdN8pcXhSUM\nOFQt9ebKndR7dUxGvNKEoRzxzuoSaht8HJvXnUFZzs+N1JrMLomkJ3morPOyv9a5rrW7DtTR6DP0\nyUgmJdG5SRnbMrpfV0b0Tmd/rY7JiGeaMJQjXszfAcClkwc4HMnhNe9aW7jXuRHfBfa1o7V0Adbv\n6rJjrX/LF5fucDgaFS6aMFTEbS2tZtm2/aQlujl7TPSMvWjNoYTh3IjvoijuIdXc+eP7keh2sXBL\nGbsOHHQ6HBUGmjBUxL1sly6+Pa4vaVE29qKlgfZUJU42fAemJxkYZWMwWuqelsjpo3phjE5IGK80\nYaiIavT5mwbrXRLF1VEBTSUMB7vWNo3BiOIqqYBAFePLy4rx+3VMRrzRhKEiav7GUvZW1zMkuwsT\nc7o5HU67omF6kMA4kIFRXiUFcMKQTPpmJLO9vJbFhfucDkeFmCYMFVEv5VtVFZdOjo51L9ozKNPZ\nrrUNXqtLrQhRM5NvW9wuaRq1/5I2fscdTRgqYnZX1DF/Uykel3DBhOibCqQ13VIT6ZaaQG2Dj7Kq\nyK/vvWN/LX4D/bqlkOSJzi61LQWqGt9Zu5sDtQ0OR6NCSROGipiX8nfg8xvOGNWLrPQkp8MJWvMB\nfJHWNCVIDFRHBQzokcq3hmbS4PXzmt1epeJDWBOGiMwUkU0islVE7mnl/REiskhE6kXkJy3eKxKR\nNSKyUkTywxmnCj+f3/DCku0AXDkl1+FoOmagg3NKxVKDd3NXTc0B4Lkl23VCwjgStoQhIm7gUeAs\nYCRwhYiMbLFbOfB94IHDnOZkY8x4Y8zkcMWpImPB5lJ2VdSR2zOVaYOjb6LBthwqYUR+LEZRFM8h\n1ZZTj7FKkVtLq1latN/pcFSIhLOEMQXYaowpMMY0AC8As5rvYIwpNcYsBZxf0kyF1XNfWqWLy4/N\nicqJBtsSWIPCidHeRXtjYwxGSwluF5dOttqpnvtym8PRqFAJZ8LoBzTvJlFsbwuWAT4UkWUiclNI\nI1MRtevAQT7eWEqCW7hkcmw0djcXqJLa5sD63k0ljBirkgLr5kDEavzeX6ON3/Egmhu9TzDGjMeq\n0rpdRE5sbScRuUlE8kUkv6ysLLIRqqC8lL8Dv4EzRvUms0vsNHYH5Pawvqy3l9dGtD6+wetn14GD\niED/7rFVwgCr8fvEoVk0eP265necCGfC2Ak0H8rb394WFGPMTvtnKfA6VhVXa/s9boyZbIyZnJWV\ndRThqnDw+vxNk9FdOSXH4WiOTEZqAl2TPdQ2+NgXwTvlXQcO4jfQNyOFRE8039sd3hVTtPE7noTz\nf+FSYKiIDBSRROByYE4wB4pImoikB54DZwBrwxapCpsPN+yhpKKOgZlpUbmqXrBye0a+WmqbvQZG\nTo/YK10EnHpMNr27JlNQVsPnW3Xkd6wLW8IwxniBO4B5wAbgJWPMOhG5RURuARCR3iJSDNwF/EpE\nikWkK9AL+ExEVgFLgLeNMe+FK1YVPv/5wmrwvOb43Jhr7G4u8KW9vTxyXWu3x0HCSHC7mrrY/mdR\nkaOxqKPX7lShIjIdWGmMqRGRq4GJwEPGmHa7Phhj3gHeabHtsWbPd2NVVbVUCYxr7/wqum3eU8Wi\ngn2kJrq5aFLsNXY3F5iWY/u+yE3bvd1u8I6FKUHacvmUHB7+eAsfbdjDjvJaBsRwAuzsgilh/BOo\nFZFxwI+Br4CnwxqVigtPLyoCrCVYuyYnOBrL0Qrc5W/TEkaHZaUncc6YPvgNPGt3r1axKZiE4TVW\na9Us4BFjzKNAenjDUrGusq6xaVqIa47PczaYEMi1v7QDa2tHQqC9JDfGSxgA107LA+DFpdupa9Q1\nv2NVMAmjSkR+DnwHeFtEXEBs3y6qsHt1WTG1DT6OH9STYb1i//4iUI0SqUZvY0zclDAAxg/oxtj+\nGeyvbeStVbucDkcdoWASxmVAPfDdZm0O94c1KhXT/H7DfxdZTVzXTouteaMOp2+3FBLcQmlVPQcb\nwn+HvK+mgdoGH12TPXRLTQz79cJNRJpKmv9ZVKRdbGNUuwnDThKvAoERV3uxxkUo1ar5m0op2FtD\n34xkTjuml9PhhITbJU2D53bsD38p41B1VOyN8D6cc8f2oWdaImt3VrKksNzpcNQRaDdhiMiNwCvA\nv+xN/YA3whmUim2zPy0A4PrpA/G4Y3PAWWsiWS21I46qowKSE9xcfZxV4pz9aaHD0agjEcxf8+3A\ndKyurhhjtgDZ4QxKxa61OytYXFBOlyQPl02J/jW7OyK3aSxG5EoYsd6ltqXvHJ9LosfFRxv3UFAW\n+ckc1dEJJmHU27PNAiAiHqyJAZX6hifs0sVlxw6I+a60LTUN3ovAuhjx1ODdXGaXJC6c0A9j4MnP\ntZQRa4JJGAtE5BdAioicDrwMvBXesFQsKqk4yNzVJbgErp+e53Q4IRe4298WgRJGYER5bpwlDIAb\nThgIwCvLinUW2xgTTMK4BygD1gA3A+8YY34Z1qhUTHrqiyK8fsNZY/rE5Oyq7cmJYJVU4BrxOCp6\naK90ZgzPoq7Rz7O6VkZMCSZh3GmMmW2MucQYc7ExZraI/CDskamYUl3v5Xl7FO+N3xrkcDThEUgY\nxeUH8fnDVytb1+hjT2U9HpfQt1tK2K7jpO+dYP0feeqLbTqQL4YEkzCubWXbdSGOQ8W4Zxdvo7LO\ny5S8Howf0M3pcMIiLclDZpdEGnx+9lTWhe06gR5S/bun4I7hCRvbMn1IT0b26cre6npdKyOGHDZh\niMgVIvIWMFBE5jR7zMdai1spwLojfuIzqwHztpMHOxxNeOVEoGvtoR5S8TMGoyURafq/8tiCr/D6\n/A5HpILR1my1XwAlQCbwl2bbq4DV4QxKxZZXlhVTVlXPqL5dOWlYfC9ildszjeXbD7CjvJbjB4dn\nfY9Ao3o8Nng3d9boPgzM3Ezh3hrmri7h/AkdWcFZOeGwJQxjzDZjzCfGmOONMQuaPZbba10ohdfn\n57EFXwFw24whiMRnFUrAgAjMWhuPg/Za43YJt5xktWX845Ot+MPYLqRCI5iR3seJyFIRqRaRBhHx\niUhlJIJT0e+t1bso3n+QQZlpzBzd2+lwwi43IlVS8bEORjAumNCfPhnJbN5TzUcbS50OR7UjmEbv\nR4ArgC1ACvA94NFwBqVig99v+OcnVunilpMGx20DbXNNCymFsWttvA7aa02ix8X37F51j8zfqpMS\nRrmgJvoxxmwF3MYYnzHm38DM8IalYsF763azeU81fTKSO039c2BtinCVMHx+w47yg1+7Vry7YsoA\neqQlsmrHARZu2et0OKoNwSSMWhFJBFaKyJ9F5EdBHqfimM9v+NsHmwG4/eQhJHo6x3+JrC5JpCa6\nqTjYyIHa0I9SLqk4SIPPT3Z6EqmJ7a6gHBdSEz3cdKJVyvjrB5u1lBHFgvkr/4693x1ADTAAuCic\nQano9/aaEraUVtOvWwqXTo6vSQbbIiJNU46Ho5QROGdeHHepbc01x+fS0y5lzN+kbRnRqs2EISJu\n4I/GmDpjTKUx5l5jzF12FZXqpHx+w0MfWqWLO0/pPKWLgEDDd1EYJiGM11lq25Oa6OHWGda4jL99\nsEVLGVGqzb90Y4wPyLWrpJQC4K1Vu/iqrIYBPVK4aFJ/p8OJuNzM8LVjBHpI5XWyhAFw1dRcMrsk\nsWZnBR9u0FJGNArm1rAA+FxEfi0idwUe4Q5MRSevz89DH20B4M5ThpIQRwskBStQXRSOEkbgnPG0\n0l6wUhLd3NZUytis4zKiUDB/7V8Bc+1905s9VCf08rJiCvfWkNszlQs7Sc+olsLZU6qztmEEXDk1\nh15dk1hfUsncNSVOh6NaaLcbhjHm3kgEoqJfbYO3qWfUT84YHlfLr3ZEXpgavY0xTSWMztaGEZCc\n4OaHpw3j56+t4YF5m5g5qnenayOLZvovoYL278+LKK2qZ2z/DM4Z08fpcBzTu2syiR4Xe6vrqa4P\n3Sw5pVX11DX66ZGWSEZKfK1W2BGXTOrP4Kw0tpfX8pyulxFVNGGooJTXNPCYPar7npkjcHWCUd2H\n43JJs1lrQ9eO0dRDqhOM8G6Lx+3iZzNHAPDwx1upqmt0OCIVEMxcUtOD2abi298/3kJVvZeThmUx\nbUim0+E4Li8M7RhFnbiHVEunj+zFpNzulNc0MHthgdPhKFswJYy/B7lNxalt+2p4ZvE2RGi68+vs\ncsPQU2pbJ+4h1ZKI8Iuzrf9rsz8tDOuCVSp4h230FpHjgWlAVotutF0Bd7gDU9HjD29voNFnuGhi\nf0b27ep0OFEhUArYHtISht1DKlNLGACTcnswc1Rv3lu3m/99dyN/vWy80yF1em2VMBKBLlhJpXl3\n2krg4vCHpqLBp1vK+GD9HtIS3fxs5nCnw4kaWsKIjF+cfQyJHhevrdjJ8u37nQ6n0ztsCcMYswBY\nICJPGWO0q0In1Ojz87u31gNw+ylDyO6a7HBE0SPUXWuNMWzb27nHYLQmp2cqN35rII/O/4p756zj\n9dumd+oOF04Lpg0jSUQeF5H3ReTjwCPskSnHPbt4G1tKq8npkcp3pw90Opyo0rdbMh6XUFJRR12j\n76jPV17TQFW9l/RkD91TO2+X2tbcNmMIvbomsaq4gtdW7HQ6nE4tmITxMrAC+BXw02YPFcfKaxr4\nqz1I71fnHENygjZbNedxu+jfPQUIzWJKTet490yN+2VuOyotycM9Z1kN4P/73kbtZuugYBKG1xjz\nT2PMEmPMssAj7JEpR9339gYq67x8a2gmp4/s5XQ4USmU05xr+0XbZo3rx8ScbpRV1fOX9zc7HU6n\nFUzCeEtEbhORPiLSI/AI5uQiMlNENonIVhG5p5X3R4jIIhGpF5GfdORYFT6LvtrHq8uLSfS4+N2s\n0XrHexiHxmIcfcN3UVP7hfaQao3LJfzh/DG4XcJ/FhWxuviA0yF1SsEkjGuxqqC+AJbZj/z2DrLX\n0ngUOAsYCVwhIiNb7FYOfB944AiOVWFQ7/XxyzfWAHD7jCEMzNQ73sMJZU8pLWG0b2Tfrnx3eh7G\nwC9eX4PX53c6pE6n3YRhjBnYymNQEOeeAmw1xhQYYxqAF4BZLc5daoxZCrSslGz3WBUe/1pQQEFZ\nDYOy0rhlRjD/zJ1XXgjXxSjq5LPUBuuHpw2jb0Yya3dW8vQi7bwZacFMDXJNa48gzt0P2NHsdbG9\nLRhHc6w6QgVl1Twy31pM8Q/njybJow3dbQmUBgrKjq6EYYyhcK9OCxKMtCQP984aDcBf3t/ErgMH\nHY6ocwmmSurYZo9vAf8DnBfGmDpERG4SkXwRyS8rK3M6nJjl8xvufmU1DV4/F03sz7TBOl9Ue3J7\npJLgFnYeOEjNUcxaW1ZdT8XBRrome8hKTwphhPHp9JG9OHNUL2oafNzz2hpdzjWCgqmSurPZ40Zg\nItYI8PbsBAY0e93f3haMoI81xjxujJlsjJmclZUV5OlVS099UUT+tv1kpyfxm3O1uSgYHreLQZnW\nn8LW0uojPs+WPdaxQ3ulaweDIP3+/NFkpCSwcHMZL+cXOx1Op3Ek05vXAMGM4loKDBWRgfaa4JcD\nc4K8xtEcqzqocG8N98/bCMAfLxhDhg4cC9rQXlbC2Lyn6ojPETh2WK9g7sMUQHZ6MveeNwqA389d\nT0mFVk1FQjBtGG+JyBz78TawCXi9veOMMV7gDmAesAF4yRizTkRuEZFb7HP3FpFi4C7gVyJSLCJd\nD3fskX5IdXh+v+HuV1ZR1+jnwgn9OE3HXHTIsF7WasVbjqKEsTlQwsjWlY87Ytb4vpw+shdV9V7u\neVWrpiKh3SVa+XqXVy+wzRgTVBnQGPMO8E6LbY81e74bq7opqGNV6M3+tIClRfvJSk/iN9/WqqiO\nGhaCEsaWphKGJoyOEBHuO380SwrLWbC5jOeWbOeqqblOhxXXgmnDWABsxJqptjvQEO6gVGSs3VnB\nA+9vAuDPF42lW2qiwxHFniF2qSDQDtFRxpim0slQrZLqsOyuydx3gdVr6vdz17O19MgTt2pfMFVS\nlwJLgEspU2PhAAAe7UlEQVSAS4EvRUSnN49xtQ1evv/CChp9huum5XHyiGynQ4pJeT1TSXS7jrin\nVFnVoR5S2dpD6oicO7YvF03sT12jn+8/v5J679FPBqlaF0yj9y+BY40x1xpjrsEaVPfr8Ialwu33\nczdQUFbD8F7pTRO7qY7zuF0MyrLGYxxJO0ag/WKY9pA6KvfOGkVOj1TWl1TywLxNTocTt4JJGC5j\nTGmz1/uCPE5FqXfXlPD8ku0kelw8dMV4nYn2KA212x6OpB0jcMxQbb84Kl2SPDx4+XjcLmH2p4V8\nsqm0/YNUhwXzxf+eiMwTketE5DrgbeDd8IalwqVwbw0/fWU1AL84awQjeuuSq0draPaRj8Voar/I\n1vaLozUxpzs/Om0oAD96cSU7dRR4yAXT6P1T4F/AWPvxuDHm7nAHpkLvYIOPW59ZRnW9l3PG9OHa\naXlOhxQXjqanlPaQCq3bZgxhxvAs9tc2cvuzy2nw6gSFoXTYhCEiQ0RkOoAx5jVjzF3GmLuAMhEZ\nHLEIVcj85s21bNxdxaDMNP7fRWO0zjxEAtVJHe0pZYzRQXsh5nIJf7t0PH0zklm54wB/fGeD0yHF\nlbZKGA8Cla1sr7DfUzHkuS+38/KyYpITXPzj6omkJ+to7lDJ7XGop1R1B3pKlVbVU1nn1TmkQqx7\nWiKPXjWRBLfw1BdFvLlSl3UNlbYSRi9jzJqWG+1teWGLSIXclwX7+M2bawFr6g9ttwit5j2lOtKO\nsUV7SIXNhJzuTXOi3f3Kal1wKUTaShjd2ngvJdSBqPDYUV7Lrc8ux+s33PitgVw4sdWB9eooHUlP\nKe0hFV5XH5fLFVMGUO/1c9PTyyitrHM6pJjXVsLIF5EbW24Uke9hrbqnolxNvZcbn86nvKaBk4Zl\ncc9ZxzgdUtwaZvdy2tKBhLGlVNsvwklEuPe80Ryb153dlXXc9N9l1DXqoL6j0VbC+CFwvYh8IiJ/\nsR8LgBuAH0QmPHWkvD4/dz6/oqmR++ErJuB2abVHuByatTb4KimddDD8Ej0u/nn1JPp1S2HljgP8\n5OVV+P06SeGROmzCMMbsMcZMA+4FiuzHvcaY4+1JA1WUMsbwmznr+HhjKd1SE3ji2slkpGgjdzgF\nusVu3F0Z1KypPr9h824tYURCZpcknrh2Ml2SPMxdXcL/vrfR6ZBiVjDjMOYbY/5uPz6ORFDq6Pzj\nk6947ktrJPcT10xmUJZ+IYVbXs80uqcmsKeynu3l7a/xvaGkkqp6L/27p5DdNTkCEXZux/Tpyj+v\nnojHJfxrYQFPLypyOqSYpFN8xJlXlxVz/7xNiMBDl41ncl4Pp0PqFFwuYerAngAsLtjX7v6BfY4b\n1DOscalDvjU0iz9dOAaA385Zx7trShyOKPZowogj763dzU9fWQXAr84ZyVlj+jgcUedy3CArOS8u\nKG9338A+mjAi65LJA7jr9GEYA99/YQULN5c5HVJM0YQRJz7dUsb3n1+B38D3TxnCDScEs4quCqXj\nBh8qYbTVjuHzG5YUWiWMqQO1BBhpd54yhOun59HoM9z033zyi9pP8MqiCSMO5BeVc9PTy2jw+blu\nWh4/On2Y0yF1SsOy0+memkBJRV2b7RgbSiqprLPaLwb0SI1ghAqs7ra/PmckF0+y1tC4/qmlrN1Z\n4XRYMUETRoxbWlTOtU8u4WCjj4sm9uc3547UUcMOCbYdQ9svnOdyCf/vwjGcNbo3VXVernriS00a\nQdCEEcOWFFrJoqbBx/nj+/Lni8fi0rEWjgq0Yyz6ShNGtPO4XTx0+QTOGNmLioONXDl7MWuKNWm0\nRRNGjFpcsI/r/r2E2gYfF0zox18uHa8D86LA8YMzAatRu7V2DJ/f8GVhoMFb2y+cluhx8ciVEzlz\nVC8q67xc9cRiVu7QeacORxNGDPpowx6ufdJKFhdO6McDl4zTZBElhmZ3oUdaIrsr69i275vtGBtK\nKqmq8zKgRwr9u2v7RTQIJI2Zo3pbSWP2Yr7YutfpsKKSJowY88aKndz032XUe/1cOTWH+zVZRBWr\nHSPQvfab1VJN1VEDtToqmiS4Xfz9ygmcP74vNQ0+rvv3Uuat0wktWtKEEUOe/KyQH764Ep/fcOuM\nwdx3/mhNFlEo0DbRZsLQ9ouok+B28ddLx3PN8bk0+Pzc+swyXliy3emwooomjBjg8xv+Z846fjd3\nPQD3nDWCn80cob2holQgGSzYXMb+moam7aWVdU2N4VO1/SIquVzCveeN4vunDMFv4J7X1vDAvE1B\nzQ/WGWjCiHKBdbif+qKIBLfwt8vGcctJukJuNBvWqwuTc7uzv7aRX7y+BmMMfr/hxy+voqbBx4nD\nsrT9IoqJCHedMZz7LhiNS+CR+Vv54Ysrqffq1OiaMKJYScVBLv3XIt5fv4euyR6e/u5ULpigCyBF\nOxHhb5eNp0uSh3fX7ual/B38+4siPt2yl+6pCdx/8VinQ1RBuGpqLv937bGkJrp5c+Uurpr9JWVV\n9U6H5SiJp6LW5MmTTX5+vtNhhMSybeXc/N/l7K2uZ0CPFJ689lhdmS3GvL6imB+9uIqUBDc+v6HB\n5+fx70zijFG9nQ5NdcDanRV87z/57K6so29GMo9fM5nR/TKcDitkRGSZMWZyMPtqCSPKGGN4fsl2\nLn98MXur65k2uCdzbj9Bk0UMOn98P84b15eDjT4afFavNk0WsWd0vwzm3DmdiTnd2FVRx8WPfcEb\nK3Y6HZYjNGFEkYMNPn788ip+/toaGn2G66bl8Z/vTqF7WqLToakjICL84YLRjO2fwZS8Hvz6nJFO\nh6SOUHZ6Ms/fdByXTrbmn/rhiyv51RtrOl27hlZJRYmvyqq57ZnlbNpTRXKCi/vOH8NFk7S9Ih4E\n/sa0V1vsM8bw7Jfb+d1b62nw+RnbP4NHr5wY05NIapVUDDHG8NLSHZz78Gds2lPFoKw03rz9BE0W\ncURENFnECRHh6uNyefXWafTvnsLq4grOfuhT5qza5XRoEaEJw0EVtY3c8dwK7n51NQcbfcwa35c5\nd5zA8N7aXqFUNBvTP4O37/wWZ47qRVW9l+8/v4Ifv7SK6nqv06GFlSYMh3yyqZQzH1zI22tKSEt0\n87fLxvHQ5RPokuRxOjSlVBAyUhN47OpJ/PGCMSQnuHh1eTEzH1zY5kzFsU7bMCKsut7LfW+v5/kl\nOwCYkNONBy8bT27PNIcjU0odqa2lVfzghZWs21UJwHXT8vjZzBGkJLodjqx9UdOGISIzRWSTiGwV\nkXtaeV9E5GH7/dUiMrHZe0UiskZEVopIdGeBIH20YQ9n/HUBzy/ZQaLbxT1njeCVW6ZpslAqxg3J\nTueN26fzg1OH4nEJT31RxJkPLuSzLfE1623YShgi4gY2A6cDxcBS4ApjzPpm+5wN3AmcDUwFHjLG\nTLXfKwImG2OC/o1HawmjtKqOe+es5+01JQCM7Z/BA5eMY5iOrVAq7qzdWcFPXl7Fxt1VAFw4sR+/\nOmckPaK0e3y0lDCmAFuNMQXGmAbgBWBWi31mAU8by2Kgm4j0CWNMEeX1+Xnys0JOfWABb68pITXR\nza/PHcnrt03XZKFUnBrdL4O37jyBn545nESPi9eW7+SUv3zCs19uw+eP7SaAcCaMfsCOZq+L7W3B\n7mOAD0VkmYjcdLiLiMhNIpIvIvllZWUhCDs0Fhfs49y/f8bv5q6nqt7LqSOyef9HJ3LDCQN1SnKl\n4lyC28XtJw9h3g9P5IQhmRyobeSXr6/l/Ec/Z/n2/U6Hd8SiuUvOCcaYnSKSDXwgIhuNMQtb7mSM\neRx4HKwqqUgH2VLR3hr+9O4G5q3bA8CAHin89txRnDayl8ORKaUibWBmGv+9YQrvrNnNH95ez5qd\nFVz4jy84b1xf7p45POZmLQ5nwtgJDGj2ur+9Lah9jDGBn6Ui8jpWFdc3Eka02FddzyPzt/LM4m00\n+gwpCW5uOWkwN580iOSE6O8poZQKDxHhnLF9mDE8i0fnb+WJzwqZs2oX763bzXenD+TWkwaTkZrg\ndJhBCWejtwer0ftUrCSwFLjSGLOu2T7nAHdwqNH7YWPMFBFJA1zGmCr7+QfA74wx77V1TScavavr\nvTzxaQGzFxZQ0+BDBC6e2J+fnDmcXl2TIxqLUir6Fe+v5c/vbWoaHd412cOtM4Zw3bQ8R7rhdqTR\nO6zjMOxeUA8CbuBJY8x9InILgDHmMbHmS3gEmAnUAtcbY/JFZBDwun0aD/CcMea+9q4XyYRRU+/l\n6UXbeHzhV+yvbQRgxvAsfnrmcEb1jZ+pj5VS4bFyxwH+992NLLKX7c1KT+KWkwZz1dSciNZKRE3C\niLRIJIyqukaeWbyd2Z8WUG4vvzkptzt3nzmcqbpOs1KqA4wxfLZ1L/fP28Tq4grAShw3nziIK6fm\nkJoY/mZmTRhhsK+6nn9/XsTTi4qorLPmi5mQ040fnTaMbw3N1MnllFJHzBjDRxtK+duHm5tGi3dP\nTeC6aQO5dlou3VLDN4ZDE0YIbS2t5snPC3l1WTH1Xj8AU/J6cNvJgzlpWJYmCqVUyBhj+HBDKY/O\n38rKHQcASE10c+nkAVw/PS8ss0JowjhKfr9VTHzqiyI+3ljatP2UEdncNmMwk/N6HPU1lFLqcIwx\nLCrYxz8/+YpP7elFROCMkb24dloexw/qGbKbVU0YR6jiYCOvLy/m6UXbKNhbA0Cix8VFE/txwwkD\nGZKto7OVUpG1oaSSJz4tZM6qnTT6rO/r4b3S+c7xuZw/od9Rz3CtCaMDjDGs2HGA577cztzVu6hr\ntKqd+mQkc/VxuVx27AAyuySFI1yllApaaWUdz365neeWbKesqh6wqqtmje/LlVNyGdP/yHpnasLo\ngLmrd3HHcyuaXk8f0pOrpuZyxsheeNy6XIhSKro0eP28u7aEZ7/czpLC8qbtT11/LDOGZ3f4fB1J\nGNE8NUhEnDIim7yeqZw5ujdXHJtDXqZONa6Uil6JHhezxvdj1vh+bC2t4rkvd7BwSxnTBmeG/dqd\nvoQBVrWU9nZSSsWqo/kOi5bpzWOGJgulVCyL1HeYJgyllFJB0YShlFIqKJowlFJKBUUThlJKqaB0\n+m61SoVU40HY9jl8NR/KC6B6D1SXAgJdsqFLL8gcCoNPhpzjwaODQlXs0ISh1NHy++Grj2HpbCj4\nBLx1re9Xsd36uQn4/EFISIXBp8CUG2HgSdZkQUpFMU0YSh0pbwOseBoWPwb7thza3mccDD4V+k6A\n9N6QlmVtry6FqhLYucxKMHvWwsa51iPrGDj+Nhh3Jbj1z1JFJx24p1RHGQPrXoeP7oX9Rda2rv1h\nyvdg/FVW1VMwKnfBimdg6f9B9W5rW+YwOO1eGH6WljhUROhcUkqFS8kqePvHULzUep05DGb8HI45\n78hLBt4GWP8GzP8j7C+0tuWeAOf8BbJHhCZupQ5DE4ZSoVZfDZ/8CRb/A4wf0rLh5J/DhGtCV4Xk\nbYD8J2HB/8LBcnAlwPQfwIk/gYSU0FxDqRZ0ahClQmnLh/CP42DRI9br426D7y+Hyd8NbXuDJxGO\nu8U696Trwd8Inz4A/zgeCj8N3XWUOkKaMJQ6nIP74Y3b4NmLoGKH1Zh948cw80+QFMbFtFK6w7cf\nhO++D9kjrWqq/5wLc++C+qrwXVepdmjCUKo1G+bCo1Nh5bPgToLTfwff+9jq+RQpOVPhpgUw4xdW\n9VT+/8Gjx1klHqUcoAlDqeaqy+Dl6+DFq6xBdwOOg1s/t9oSnOju6kmEGT+DmxdYyaqy2CrxvH4r\n1Ja3f7xSIaQJQymwusqueBYenWJ1mU1Ig7P+DNe/a43MdlqvUXDDh1aXW08yrHrOKgGtecWKXakI\n0IShVOlGeOocePM2q3fSoBlw2yKYejO4ouhPxO2BE34It3wOOdOgphRevQGeuRD2feV0dKoTiKK/\nBqUi7OABmPdLeOwEa/6ntCy4cDZ85w3onut0dIeXOQSuexu+/TAkd7NGjf/jePjwXm0UV2GlCUN1\nPr5GWDIbHp5gdZX1e2HSdXDHUhh7aWyMsHa5YNK1cOcyGHcF+Orhs7/CwxNh+dPg8zodoYpDOnBP\ndR6+Rlj1Aiy8Hw5ss7blngAz/2h1mY1lxfnw3j2HRqD3HAIn/QxGXwQut7OxqaimI72Vaq6hBlY9\nD188cmjqjcxhcOpvYMS5sVGiCIYxsPZV+Pj3h+a4yhwG0+6EMZdCQrKj4anopAlDKbC+NJc9Bfn/\nhroD1rYeg2HGPfF9591UkvozHLCnVE/NhGNvgInXQEZ/Z+NTUUUThuq86qtg4zuw8hkoXHhoe7/J\n1vThx8zqPNOHexusLsKLHoHdq+2NYq3BMeEqGHYWJKY6GqJyniYM1blUl8HWD2DDW7D1I6sBGKzx\nCsd8G6bcBAOmOBujk4yxeoEtmQ2b3gFfg7XdkwJDT7Nm2h1yGqT2cDZO5YiOJIxOcqul4kptudXI\nu+0zq0vp7jXN3hRrdPa4y2DUhZDSzbEwo4YI5J1gPWrLrcF+q1+EnflWkt3wFiDQdzwMOhnyplsl\nMv3dqRa0hKGiW80+KF1vVamUrIZdK2Dvpq/v40mG3Gkw4hyrETu9tzOxxpqKYtj4trXi3/bFh0oe\nAAhkH2NNR9J7DPQea73WUkjciZoqKRGZCTwEuIEnjDH/r8X7Yr9/NlALXGeMWR7Msa3RhBGDvA3W\nsqWVu6wZYfcXWY/yAijbZI28bsmdBP0mWtVMg06GnOO1B9DRaqiFbV9AwXzYscRKzP7Gb+6XmglZ\nw6HHQOieB93yoNsA6NoX0vuAOyHSkaujFBVVUiLiBh4FTgeKgaUiMscYs77ZbmcBQ+3HVOCfwNQg\nj1XRwhjr7rShxmp0bqi2ftZVWI+DB6ypwmv3WQmgpsxqd6jeY22jjZuWxHTIGnboLrfPOOu5Jyli\nH69TSEy12jOGnma9bqyzVhfcvdp+rIGyzVC7F7bttdpEvkEgLRO69LJGzadlWSWS1J7WlO3J3SA5\nA5K7WtPDJ6VDYhfr4UmM6MdVRyacbRhTgK3GmAIAEXkBmAU0/9KfBTxtrGLOYhHpJiJ9gLwgjg0N\nX6P1hRb40vpaiav5tsM9t/dr73nguNZ+Gn8rz/3WZb72+nAPA36f/dxnP/eBv9lrv/frz/1e67P7\nfdadpK/R+tL3e62f3gbrp6/eeu6ts143HgRvPXgPWs8baqGxxjruSIgL0vtad6hd+1l3rYFH1nDr\nrjVexknEkoRka3r1nKmHtvn9ULXLKvkFSoL7i6Byp1VCrNpt3QzUlHX8eq4EK2klpFqrCyakWjcF\nnmTrpzvJSipNP+2Hy2M/T7DO4fZY25oebhB3i+cu+6f9Wlz2c7GeH+5B4H1avJZmr6WNnxx63aHn\nfP0cNPt7aL4tpXvYewCG8+z9gB3NXhdjlSLa26dfkMeGxu7VMPuUsJy6U3F5IDHNultMsu8gk7ta\nd5Up3ayfqT0PPbpkW3eiqT07TzfXWOdyWWM4DjeOw+e1SiDVe6C61Co9Bh4HD1hjYQ4esEqf9VVQ\nX2mXRqutm5ZAiVQdmduXWqXxMIr5v1QRuQm4CSAnJ6fjJ3B5rHpZ62SBsza/wKFtbWX7r90VtLK9\nrZ/N706+9tp1aD+X+9D2wN1Q8+3iana31OIOyuU5tM2dYL0Wt303Zr92Jxy6S2t+B+dOPHSHl5B8\n6I4vcCfoSdYqBWVxe6wOBx3tdGCMVWptrLVLsHXWc2/9oRKtr94u2dbbJWK7FBx47vNaz/1e67nx\nHSpFN5Wsfa2UwO1SeKC0Hiipg/U8ULpvet6yFqB5bYG/7ZqEb9Q+0Oz5YWooAr+fdmtAiMhA1HAm\njJ3AgGav+9vbgtknIYhjATDGPA48Dlajd4ej7DMO7tapoZVyjIh1M6IdF6JeOGerXQoMFZGBIpII\nXA7MabHPHOAasRwHVBhjSoI8VimlVASFrYRhjPGKyB3APKyusU8aY9aJyC32+48B72B1qd2K1a32\n+raODVesSiml2qcD95RSqhPryDgMXUBJKaVUUDRhKKWUCoomDKWUUkHRhKGUUioomjCUUkoFJa56\nSYlIGbDN6Tg6KBPY63QQEaafuXPQzxwbco0xWcHsGFcJIxaJSH6wXdrihX7mzkE/c/zRKimllFJB\n0YShlFIqKJownPe40wE4QD9z56CfOc5oG4ZSSqmgaAlDKaVUUDRhRBER+bGIGBHJdDqWcBOR+0Vk\no4isFpHXRaSb0zGFg4jMFJFNIrJVRO5xOp5wE5EBIjJfRNaLyDoR+YHTMUWKiLhFZIWIzHU6lnDR\nhBElRGQAcAaw3elYIuQDYLQxZiywGfi5w/GEnIi4gUeBs4CRwBUiMtLZqMLOC/zYGDMSOA64vRN8\n5oAfABucDiKcNGFEj78Bd9Ns4cZ4Zox53xjjtV8uxlpVMd5MAbYaYwqMMQ3AC8Ash2MKK2NMiTFm\nuf28CusLtJ+zUYWfiPQHzgGecDqWcNKEEQVEZBaw0xizyulYHPJd4F2ngwiDfsCOZq+L6QRfngEi\nkgdMAL50NpKIeBDrhs/vdCDhFM41vVUzIvIh0LuVt34J/AKrOiqutPWZjTFv2vv8Eqsa49lIxqbC\nS0S6AK8CPzTGVDodTziJyLlAqTFmmYjMcDqecNKEESHGmNNa2y4iY4CBwCoRAatqZrmITDHG7I5g\niCF3uM8cICLXAecCp5r47N+9ExjQ7HV/e1tcE5EErGTxrDHmNafjiYDpwHkicjaQDHQVkWeMMVc7\nHFfI6TiMKCMiRcBkY0ysTWDWISIyE/grcJIxpszpeMJBRDxYDfqnYiWKpcCV8bw+vVh3Pf8Byo0x\nP3Q6nkizSxg/Mcac63Qs4aBtGMopjwDpwAcislJEHnM6oFCzG/XvAOZhNf6+FM/JwjYd+A5wiv3v\nutK+81ZxQEsYSimlgqIlDKWUUkHRhKGUUioomjCUUkoFRROGUkqpoGjCUEopFRRNGCrsRMRnd69c\nJyKr7Fl5W/2/JyJ9ReSVEFwzT0QONuvauVJEEo/2vNFCRK4TkTJ7dtQtIjJPRKYFcdz5RzsZoIi8\nE6+zC6u26UhvFQkHjTHjAUQkG3gO6Ar8tvlOIuIxxuwCLg7Rdb8KXLc19vW8h3s/BrxojLkDQERO\nBl4TkZONMW3NmHo+MBdYf6QXNcbouIpOSksYKqKMMaXATcAdYrlOROaIyMfAR3bJYC2AiCwWkVGB\nY0XkExGZLCJpIvKkiCyx77CDngFWRP5HRP4rIp8D/7XXMLhfRJbaa3PcbO8nIvKIvZbFh/Zd9cX2\ne0WBNUvseD6xn7cal/0ZXxOR9+zSwJ+bxTNTRJbbJa+PRMRl75Nlv++y19LIauf3Oh9redCb7OMG\n29dbJiKfisgIuwRyHnC/XeIaHPid2sdk2jMNtBdzkb1vnohsEJHZdunxfRFJsfc51v59rrR/v2uD\n/TdS0UsThoo4Y0wB4Aay7U0TgYuNMSe12PVF4FIAEekD9DHG5GNN2PixMWYKcDLWF2BaK5ca3Kw6\n6tFm20cCpxljrgBuACqMMccCxwI3ishA4AJguL3vNUC71T3txDUeuAwYA1wm1kJDWcBs4CJjzDjg\nEmOMH3gGuMo+7jRgVZDTpywHRtjPHwfuNMZMAn4C/MMY8wUwB/ipMWa8Meards73jZhb2Wco8Kgx\nZhRwALjI3v5v4Ga7hOcLInYVA7RKSkWDD4wx5a1sfwl4H6vq6lIg0LZxBtZkbz+xXycDOXxz8ZrD\nVUnNMcYcbHausYHSA5CB9SV4IvC8McYH7LJLQO05XFwAHxljKgBEZD2QC3QHFhpjCgGa/Q6eBN7E\nmjL7u1hfvsEQ+/xdsBLcy2JNaAmQFOQ5mmst5h0t9ik0xqy0ny8D8uz2jXRjzCJ7+3NYk0yqGKcJ\nQ0WciAzCuusstTfVtLafMWaniOwTkbFYd7q3BE6BdVe+6QhDaH49wboTn9cixrbq6b0cKp0ntzjX\nN+ISkalAfbNNPtr42zPG7BCRPSJyCtYiTFcdbt8WJmAlTRdwoK32m2YO91kIMuaW+6QEF6qKRVol\npSLKroZ5DHgkyCnNX8RamCbDGLPa3jYPuFPs22cRmXAUIc0DbhVrSm5EZJhdjbQQqxrGbVeHndzs\nmCJgkv38ohbn6khci4ET7SowRKRHs/eewKqaetku5bRJRE7Car+Yba8/USgil9jviYiMs3etwpr0\nsbXPEpLOBsaYA0CVnSgBLg/FeZXzNGGoSEix2xHWAR9iVTPdG+Sxr2B94bzUbNvvgQRgtX3O3x9F\nbE9g9RhabjfM/gvrTvp1YIv93tPAombH3As8JCL5fL1+vkNx2e0SN2H1blqFlRwD5gBdaLs66jL7\n97oZaxGui5r1kLoKuME+7zoOLQ37AvBTu1F+MPAAVsJcAWS2FW8H3QDMFpGVQBpQEcJzK4fobLVK\nBUFEngLmGmOOeoxIkNebDPzNGPOtSFwv1ESkizGm2n5+D1aHhR84HJY6StqGoVSUsb9gbyX4toto\ndI6I/BzrO2YbcJ2z4ahQ0BKGUkqpoGgbhlJKqaBowlBKKRUUTRhKKaWCoglDKaVUUDRhKKWUCoom\nDKWUUkH5/+ANnM7K3u2JAAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZEAAAEKCAYAAADTgGjXAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xd8VfX9+PHXOxNImAlhhSRMGYICYThx0IpixVm1jmqt\n2Kq1ttrW0e+3an9++7X2a+tqLVpLHdRtQXGCG2SEvWcCJIyEJISEkHXv+/fHOdGIIblJ7s25yX0/\nH488yD3n3HPeQcn7ftb7I6qKMcYY0xxRXgdgjDGm7bIkYowxptksiRhjjGk2SyLGGGOazZKIMcaY\nZrMkYowxptksiRhjjGk2SyLGGGOazZKIMcaYZovxOoBQS05O1oyMDK/DMMaYNmX58uUHVLVnY9e1\n+ySSkZFBVlaW12EYY0ybIiI7A7nOurOMMcY0W9glERF5WEQ2icgaEXlTRLrVOXe3iGwTkc0ico6X\ncRpjjAnDJAJ8CByvqqOBLcDdACIyArgCGAlMBf4qItGeRWmMMSb8koiqfqCqNe7LxUCq+/104CVV\nrVTVbGAbMMGLGI0xxjjCLokc5UfAu+73/YDddc7luse+RURmiEiWiGQVFBSEOERjjIlcnszOEpH5\nQO96Tt2rqnPca+4FaoAXm3p/VZ0JzATIzMy0XbeMMSZEPEkiqjqlofMich1wPnC2fr31Yh7Qv85l\nqe4xY4wxHgm7dSIiMhX4NTBZVcvrnJoLzBaRR4C+wBBgqQchGmNM2PH5la35pazcdZCVu4pZtfsg\nM6/JJCM5IaTPDbskAjwBxAMfigjAYlX9iaquF5FXgA043Vy3qKrPwziNMcYzZZU1rNxVTFZOMSt2\nFbNy10HKKmu+cc3K3cWRl0RUdXAD5x4EHmzFcIwxJiwUH65iSXYRS7ILWZZTxIY9h/AfNeLbr1tH\nxqZ358T+3RiT1o0RfbqEPK6wSyLGGGOg5Eg1S3YUsmh7IYt3FLJpX+k3zsdECaNSu5KZ3p3M9O6M\nTe9Ory4dWj1OSyLGGBMGqn1+Vuws5vOtB/hi2wHW5B78RksjLiaKMf27MXFgEpMG9ODEtG50ivP+\nV7j3ERhjTITac/AIH2/O55PNBXy5vfAbYxoxUcK49G6cNCiZkwYmMSatGx1iw69IhyURY4xpJX6/\nsir3IPM37OejTfnf6qIanJLIaUOSOW1IMhMHJJEQH/6/osM/QmOMacMqa3ws2lbIe+v2sWDTfg6U\nVX11rlNcNKcOTubMYSmcPrQn/bp19DDS5rEkYowxQVZR7eOTzQXMW7uXjzflf6ObKrV7R6YM78WU\n4b0YP6A78THh10XVFJZEjDEmCKp9fj7fWsCcVXuYv2E/h6u+XsY2vE8Xpo7szXdH9mJY7864a+Da\nBUsixhjTTKrKil0HeWNFLu+s3UtxefVX50b168q00X049/jepCeFdsGflyyJGGNME+0tOcJrWbm8\nsTKP7AOHvzo+JCWRC8f04/zRfdp14qjLkogxxgSg2udnwcZ8Xl62i0+3FHy1hqNn53guGtOPC0/s\nx/A+7aurKhCWRIwxpgF7Dh7hpaW7eGnZbvJLKwGIjRbOHdmbS8elctrgZGKiw31rptCxJGKMMUdR\nVZZmFzFrUQ7vr9/3VatjcEoiV4zvz8VjU+mREOdtkGHCkogxxriqfX7eWr2HZz7PZsPeQ4Czcvy8\nUb25elI6Ewf0iLjuqsZYEjHGRLyyyhpmL9nJPxfmsLekAoCkhDiumpjGVZPSPSls2FZYEjHGRKyD\n5VX8c2EOsxblUHLEmZ47OCWRGacN5IIT+4ZlrapwY0nEGBNxDpZX8fTnO5i1MOerRYGZ6d25+cxB\nnDE0hago67IKlCURY0zEKK2o5pnPs3n2i2xK3VIkpw1J5tYzBzNxYJLH0bVNlkSMMe1eZY2PFxfv\n4vGPtn61qvy0IcncPmUo49K7exxd22ZJxBjTbqkqb6/Zy0PvbSK3+AgA4zO6c+d3j7OWR5BYEjHG\ntEtrcg/ywFsbyNpZDMDQXon8ZuowzhqWYtN0gyhsk4iI3AH8CeipqgfcY3cDNwA+4DZVfd/DEI0x\nYaj4cBUPvbeJl7N2o+pM1b3znOP4fmZ/om3APOjCMomISH/gu8CuOsdGAFcAI4G+wHwRGaqqvvrv\nYoyJJH6/8ury3fzvu5soLq8mNlq4/pQB3HrWYLp0iPU6vHYrLJMI8Gfg18CcOsemAy+paiWQLSLb\ngAnAlx7EZ4wJI9sLyrj79bUszSkC4ORBSTww/XgGpyR6HFn7F3ZJRESmA3mquvqofst+wOI6r3Pd\nY8aYCFXj8zPz8x38Zf5Wqmr8JCfG81/nD+eCE/rauEcr8SSJiMh8oHc9p+4F7sHpymrJ/WcAMwDS\n0tJacitjTJjaXlDGL19exercEgAuHZfKb6cNp1snK4zYmjxJIqo6pb7jIjIKGADUtkJSgRUiMgHI\nA/rXuTzVPVbf/WcCMwEyMzM1eJEbY7zm9yvPL97JH97dSEW1n37dOvKHi0dx+tCeXocWkcKqO0tV\n1wIpta9FJAfIVNUDIjIXmC0ij+AMrA8BlnoSqDHGEwfKKrnjldV8uqUAgEvGpvK7C0bYwLmHwiqJ\nNERV14vIK8AGoAa4xWZmGRM5Fm0/wO0vrSK/tJLunWL5w8WjmHp8H6/DinhhnURUNeOo1w8CD3oT\njTHGC36/8vhH2/jLgi2owoQBPXj0ihPp07Wj16EZwjyJGGMiW8mRan758ioWbMpHBG47ewi3nTU4\norejDTeWRIwxYWnzvlJuej6LnMJyunaM5bErxzDZBs/DjiURY0zY+XhTPrfOXsHhKh/D+3Rh5jXj\n6N+jk9dhmXpYEjHGhJVZC7N54O0N+BXOH92Hhy89gY5xtsNguLIkYowJCz6/8vu3NzBrUQ7gjH/8\nYsoQW3ke5iyJGGM8V1nj45evrGbemr3ERUfx0KWjuGhMqtdhmQBYEjHGeKqssoabns9i4bZCEuNj\nmHntOE4elOx1WCZAlkSMMZ45WF7FD59dyurcEpIT45l1/XiO79fV67BME1gSMcZ4ovhwFVc9s4QN\new/Rv0dHXrhhIulJCV6HZZrIkogxptUVllVy1TNL2LSvlAHJCcy+caKtQG+jLIkYY1pV8eEqfvD0\nEjbvL2VQzwRm3ziJXl06eB2WaSZLIsaYVlNaUc21zy5l8/5SBqckMvvGiaR0tgTSllkBGmNMqzhS\n5eOGWVmszSshPakTL/7YEkh70GgSEZGHAjlmjDHHUu3zc9MLy1maU0TvLh144YaJ1oXVTgTSEvlO\nPcfODXYgxpj2SVW56/W1fLalgKSEOF748USrg9WOHHNMRER+CtwMDBKRNXVOdQYWhjowY0z78Jf5\nW3l9RS4dY6P55/XjGZyS6HVIJogaGlhfDLwL/AG4q87xUlUtCmlUxph24ZVlu3l0wVaiBJ74wRhG\np3bzOiQTZA0lkWdUdZyIpKjqzlaLyBjTLizeUcg9b64F4IHpx3P28F4eR2RCoaEkEiUi9wBDReSX\nR59U1UdCF5Yxpi3bXVTOzS+uoMavzDh9IFdPSvc6JBMiDQ2sXwH4cBJN53q+jDHmW8qrapjx/HKK\nDlcxeWhPfjN1mNchmRA6ZktEVTcDD4nIGlV9txVjQkR+BtyCk8Tmqeqv3eN3Aze4x29T1fdbMy5j\nTMNUlV+9toaNew8xIDmBx64cQ3SU7QfSnjU0O+tqVX0BGCEiw48+H6ruLBE5E5gOnKCqlSKS4h4f\ngdM6Ggn0BeaLyFBV9YUiDmNM0/1zYQ7z1uylc3wMT1+bSdeOsV6HZEKsoTGR2nKarT0f76fA/6pq\nJYCq5rvHpwMvucezRWQbMAH4spXjM8bUY9Xug/zh3Y0APHzZaJvKGyEa6s76u/vn/a0XDgBDgdNE\n5EGgArhTVZcB/XCmHdfKdY8ZYzx2sLyKW15cQbVPuf6UDKYe38frkEwrabQAo4j0BG4EMuper6o/\nau5DRWQ+0LueU/e6z+gBTALGA6+IyMAm3n8GMAMgLS2tuWEaYwKgqtz56hryDh7hhNSu3H3ut3q/\nTTsWSBXfOcDnwHycAe0WU9UpxzrnrpR/Q1UVWCoifiAZyAP617k01T1W3/1nAjMBMjMzNRgxG2Pq\n9+KSXczfuJ8uHWJ44gdjiYuxuq6RJJAk0klVfxPySL72H+BM4GMRGQrEAQeAucBsEXkEZ2B9CLC0\nFeMyxhxlR0EZD85zxkH+5+JRVhMrAgWSRN4WkfNU9Z2QR+N4FnhWRNYBVcAP3VbJehF5BdgA1AC3\n2MwsY7xT4/Pzi1dWc6Tax0Vj+nH+6L5eh2Q80NAU31JAAQHuEZFKoNp9raraJRQBqWoVcPUxzj0I\nPBiK5xpjmubJj7ezevdB+nbtwH0XjPQ6HOORhmZn2ap0Y5pAVZm/MZ//fXcjPRLimHX9BBLi2+fm\noevySnjso60A/On7J9h6kAgWyKZUp4hIgvv91SLyiIjYlCdj6sg+cJjrZy3jxuey2F5wmGU5xfzi\n5VX4/e1vXke1z8+vX1uDz+9M5z15ULLXIRkPBTKN4m9AuYicANwBbAeeD2lUxrQhucXlXPD4F3yy\nuYDOHWL45XeG0rlDDB9s2M+f52/xOryge+bzbDbsPURq94786pzjvA7HeCyQJFLjDmxPB55Q1Sex\nAozGAE4X1n/PWU9pZQ2nDk7mozvO4Lazh/DkD8YSJfD4R9t4a/Uer8MMmuwDh/mLmxj/56JRdIpr\nn911JnCBJJFSt/Dh1cA8EYkCrAPUGOCdtfv4aFM+nTvE8Mj3T6Bn53gATh/ak99OGwHA3W+spaK6\n7U8k9PuVu15fQ2WNn4vH9uP0oT29DsmEgUCSyOVAJXCDqu7DWeT3cEijMqYNKDlSzX1vrQfgN1OH\nkdKlwzfOX39KBkN7JVJWWcOa3BIvQgyq11fksiS7iKSEOP7LTZDGNJpEVHWfqj6iqp+7r3ep6nOh\nD82Y8PbQe5soKK1kXHp3fjDh23NNRIRJA5MAWJpd2NrhBVXJkWoeem8TAPdOG073hDiPIzLhwuoT\nGNMMm/eVMnvJLmKjhT9cPIqoY+yZMWFADwCWZBe1ZnhB95f5WzhQVsX4jO5cNMbqnpqvWRIxphn+\n9sk2AK6ckMbQXseeZzIhw0kiy3cWU+Pzt0pswbZp3yGe+3InUQL3X3A8IrbJlPlaQElERDqKiM3l\nMwbYVVjO3NV7iIkSZpzecIHplC4dGJCcQHmVj/V7DrVShMFTO/vM51eumZTOiL4hKVRh2rBAFht+\nD1gFvOe+PlFE5oY6MGPC1VOfbcevcOGYfqR2b7zgYG1rZGkb7NKat3YvS93B9F9+xz5Hmm8LpCVy\nH84OggcBVHUVMCCEMRkTtvIPVfBaVi4i8JPJgwJ6T1sdF6ms8fHH9zYD8MvvDqVrJ5vZb74tkCRS\nrapHz09sf7UcjAnAM19kU+XzM3Vk74C3f61NIstyitpUGZTnv9zJrqJyBqckcnlm/8bfYCJSIElk\nvYj8AIgWkSEi8jiwKMRxGRN2SsqreWHxTgBuPmNwwO9L7d6Rvl07UHKkmi35paEKL6hKyqt5/CNn\n8sDd5w4jJtrm4Jj6BfJ/xs+AkTgLDv8NHAJuD2VQxoSjl7N2UV7l45TBSYxK7Rrw+0Tk69ZIG+nS\nevKTbZQcqeakgUmcNSzF63BMGAtksWG5qt6rquOBicBDqloR+tCMCR8+v/Lcl04r5PqTmz4kOGGA\ns+iwLYyL7C4qZ9bCHADuOW+4Tek1DQpkdtZsEeniloNfC2wQkV+FPjRjwsdHm/LJLT5CWo9OnNmM\nT+a1LZGl2UU49UzD12MLtlLl8zP9xL5NanGZyBRId9YIVT0EXAi8izMz65qQRmVMmPnXohwArj0p\nnehjrE5vyKCeCXTtGEt+aSX5pZVBji54dhSU8fqKXKKjhF9MGep1OKYNCCSJxIpILE4Smauq1djs\nLBNBtu4v5YttB+gYG81lzZylJCIMSE4AYGdheTDDC6pHF2zFr3DZuFQy3HiNaUggSeTvQA6QAHwm\nIuk4g+vGRIR/fZkDwEVj+7VoG9iMJGdhYk7h4SBEFXyb95Uyd/UeYqOFW88KfPaZiWyN7iijqo8B\nj9U5tFNEzgxdSMaEj0MV1byxIg+AH56U0aJ7pSU5n+x3hWlL5C/zt6Dq1AMLZCW+MRBAEgEQkWk4\n03zrbpjwQCgCEpETgafcZ9UAN6vqUvfc3cANgA+4TVXfD0UMxtR6Y3ku5VU+ThqYxHG9W7ahZzi3\nRNbllfDuun3Ex0Rxy5nWCjGBazSJiMhTQCfgTOAZ4FJgaQhj+iNwv6q+KyLnua/PEJERwBU4yawv\nMF9Ehqpq298yzoQlVeWFJbsAZ0C9pdKTwndM5MmPnYWFV01Mp9dRm2sZ05BAxkROVtVrgWJVvR84\nCQjltA0FakuFdgVqN6ieDrykqpWqmg1sw6npZUxILN5RxLb8MlI6xzNlRK8W3y+9TksknKb5bssv\n5b31+4iLjuKmyQ1XJTbmaIF0Zx1x/ywXkb5AIdAndCFxO/C+iPwJJ8md7B7vByyuc12ue+xbRGQG\nMAMgLe3bO84ZE4jaEidXTkgjNghlP5IS4kiMj6G0ooaD5dVhszvgXz/ejipcmplqrRDTZIH8y3hb\nRLrh7Ku+Amem1r9b8lARmS8i6+r5mg78FPiFqvYHfgH8o6n3V9WZqpqpqpk9e/ZsSagmQuUfquD9\n9fuIjhKurGfr2+YQkW+0RsLB7qJy5qzeQ3SU8NMAqxIbU1cgs7N+7377uoi8DXSop6pvk6jqlGOd\nE5HngJ+7L1/FGYcByAPqTtJPdY8ZE3QvLdtNjV+ZOrI3vbsG79N5elIn1u85xM7CcsakdQ/afZvr\nqU+34/MrF4/pR/8eNiPLNF2gs7NOBjJqrxcRVPW5EMW0B5gMfAKcBWx1j88FZovIIzgD60MI7QC/\niVA1Pj//XuoMqF89qeUD6nWF0+D6/kMVvOrujXLzmdYKMc0TyOys54FBOLsb1s6EUiBUSeRG4FER\niQEqcMc2VHW9iLwCbMCZ+nuLzcwyobBgUz57SyoYkJzAyYOSgnrv2mm+O8OgO+vZhXX3RmnZ9GUT\nuQJpiWTi1M9qlekkqvoFMO4Y5x4EHmyNOEzket6t1nv1pHSimlEnqyG1LRGvx0TKKmuY7U5fthlZ\npiUCGVhfB/QOdSDGhIPtBWV8se0AHWKjuHRsatDvXzuwvqvI2+6sl5ftprSihvEZ3cNibMa0Xcds\niYjIWzjdVp1xyr8vxdmYCgBVvSD04RnTumqn9U4/oV9I9hTv1bkD8TFRHCirorSims4dWn/f8hqf\nn2e/yAbgxtOsFWJapqHurD+1WhTGhIHyqhpeW54LwDVBWKFen6goZ5rvlv1l7Cws5/h+rb9fx7vr\n9pF38AgDkhOYMrzliyhNZDtmd5aqflr7BWzGWT3eBdjsHjOmXZmzag+lFTWMSesW0l/uaT3cQowe\ndGmpKjM/2wHADacOCPqYj4k8gexs+GOcqbQX49TNWiwiPwp1YMa0JlX9akA9GHWyGuJlIcYl2UWs\nzSuhR0Icl4RgzMdEnkBmZ/0KGKOqhQAikgQsAp4NZWDGtKYVu4rZsPcQPRLiOPf4UFb1gfTazakO\ntH5L5J8LnbGQqyel0zEuutWfb9qfQGZnFQKldV6XuseMaTf+uTAHgMvH96dDbGh/uab38KYlsruo\nnA837Cc2Wrh6ktWUM8ERSEtkG7BERObgzNaaDqwRkV8CqOojIYzPmJDbV1LBe+ucOlnXBHmFen0y\nkrwZE3lh8U78CtNH9yWlsxVaNMERSBLZ7n7VmuP+aUtcTbvw4pKd1PiVaaP60Ldbx5A/r2+3DsRE\nCXtLKqio9oW85QPOzLPaUi7XnZwR8ueZyBFIAcb7WyMQY7xQUe37auX2dadktMozY6KjSO3ekZzC\ncnYVlTO0V+g/j/1n5R4OuTPPTujfLeTPM5Gj5ZskGNOGvb1mL4WHqxjZtwuZ6a23crs1CzGqKrMW\nOQPq1goxwWZJxEQsVf1qttJ1J2cg0nprJlqzEOOi7YVs2e/s0BjqmWcm8gSyTiS4ZUyNCRNZO4tZ\nv8eZ1vu9E/q26rPTWrEQ478W5QDO/ulxMfa50QRXIP9HLRaRV0XkPGnNj2rGhNgznzsrt6+amNYq\ng9t1fd0SCW131p6DR5i/cT8xUcKVE/s3/gZjmiiQJDIUmAlcA2wVkf8RkaGhDcuY0Mo+cJgPNuwn\nLjqKa0/KaPXnt9aYyOwlu/ArTD2+t03rNSHRaBJRx4eqeiXOhlE/BJaKyKciclLIIzQmBP7xxQ5U\n4aIx/ejZOb7Vn9+/R0dEILe4nKoaf0ieUVXj56VlzswzLxKliQwBjYmIyM9FJAu4E/gZkAzcAcwO\ncXzGBF3R4aqvqvX++LQBnsQQHxNN364d8SvkHTwSkme8u24vB8qqOK5XZ8Zn2J4hJjQC6c76Eqd6\n74WqOk1V31DVGlXNAp4KbXjGBN8Li3dSUe3nzON6MqQV1mgcS3qIZ2jV7o1yzUnprTrzzESWQJLI\nb1X196qaW3tARC4DUNWHQhaZMSFQUe3juS9zALjxdG83ZArluMjGvYdYllNMYnwMF47pF/T7G1Mr\nkCRyVz3H7g52IMa0hjdW5HGgzFlceNJAb2evp4ewJHxtK+Tisf1IjA+kupExzXPMJCIi54rI40A/\nEXmsztcsoKYlDxWRy0RkvYj4RSTzqHN3i8g2EdksIufUOT5ORNa65x6z6camqWp8fp761CkDd9Pk\nQZ538YRqmm9pRTVvrswDaJWCkiayNdQS2QNkARXA8jpfc4FzGnhfINbhbHL1Wd2DIjICuAIYCUwF\n/ioitRP4/4YzO2yI+zW1hTGYCPPWmj3sKipnQHIC00Z5v3L76+6s4LZE3lyZR3mVj4kDeng65mMi\nwzHbuaq6GlgtIi+qaotaHvXceyNQ3yfB6cBLqloJZIvINmCCiOQAXVR1sfu+54ALgXeDGZdpv/x+\n5a8fO62Qn04eRHQYbAub5u4rsrvoCD6/BiUmVf3GgLoxoXbMJCIir6jq94GVIqJHn1fV0SGIpx+w\nuM7rXPdYtfv90ceNCcgHG/axNb+Mvl07hM1Ac0J8DD07x1NQWsnekiOkdu/U4nsuzS5iy/4ykhPj\n+e6I3kGI0piGNTTi9nP3z/Obc2MRmQ/U93/xvao6p57jQSMiM4AZAGlptoNbpFNVnvh4G+CMhYRT\n/aiMpE4UlFays7A8KEnkBbes/ZUT+ofVz2nar4a6s/a6316C08W0pyk3VtUpzYgnD6hb4CfVPZbn\nfn/08WM9eyZOqRYyMzO/1YoykeXTLQWsyztEcmIcl48Pr/pRaT0SWJZTzM7Cck4Z3LJ7FZRW8t66\nvUQJXDnBPjyZ1hHIR5XOwIci8rmI3CoivUIYz1zgChGJF5EBOAPoS92EdkhEJrmzsq7l6x0WjTkm\nVeXPH24B4MenDWz1QouNCWZJ+FeydlPtU84e3qtVdmg0BgKrnXW/qo4EbgH6AJ+6XVXNJiIXiUgu\ncBIwT0Ted5+1HngF2AC8B9yiqj73bTcDz+Ds+b4dG1Q3AViwMZ/VuSUkJ8ZxbRgONKcnB6ckfI3P\n//WAuk3rNa2oKauQ8oF9QCGQ0pKHquqbwJvHOPcg8GA9x7OA41vyXBNZ/H7l/9xWyE/PGEynuPBb\ndBestSIfbNjP3pIKBvZM4NTBycEIzZiABFKA8WYR+QRYACQBN4ZoZpYxQfXe+n1s3HuI3l06cNXE\n8BwjSO/xdekT1eYP381yN5764UkZRIXB9GUTOQL5aNYfuF1VV4U6GGOCxedXHnFbIbecNTjsxkJq\nde0US7dOsRwsr2b/oUp6d236nh8b9hxiaXYRifExXDIutfE3GBNEgYyJ3K2qq0QkRUTSar9aIzhj\nmmvOqjy25ZfRr1tHLs8MrxlZRxveuwsAq3YfbNb7a7e/vXRcqtXJMq0ukO6s74nIViAb+BTIwQa1\nTRirqPbxfx84rZDbpwwJ+/USEwb0AJyFgk1VfLiK/6xyZruH48QB0/4F8q/r/wGTgC2qOgA4m2+u\nKjcmrMxalEPewSMM692Zi8eGf/fOxNokklPY5Pe+tGw3lTV+Jg/tycCeicEOzZhGBZJEqlW1EIgS\nkShV/RjIbOxNxnih+HAVT7qr0+8+b3hY1MhqzJi07sRECRv2HOJQRXXA76uo9vHPhdkAXHdKRoii\nM6ZhgSSRgyKSiFNx90UReRQIzVZsxrTQYx9tpbSihtOGJDN5aE+vwwlIx7hoRqV2xa+wfGdxwO97\nfUUu+aWVDO/ThTPayM9q2p9Aksh0oBz4Bc4CwO3A90IZlDHNsbPwMC8s3okI3HPecK/DaZKmjovU\n+Pz8/dMdANx8hvd7o5jIFUgSSQHi3H3V/wU8jVMKxZiw8v/mbaTap1wyNpXhfbp4HU6TTGxiEpm3\ndi+7isrJSOrEeWGwN4qJXIEkkVcBf53XPveYMWHj4835fLhhPwlx0fzqnOO8DqfJxqX3QATW5B6k\notrX4LV190b5SZjsjWIiVyBJJEZVq2pfuN/HhS4kY5qmssbH/XPXA3D7lKH06tL0BXte69oxluG9\nu1DtU1buani9yEeb8tm8v5TeXTpw0djw2BvFRK5AkkiBiFxQ+0JEpgMHQheSMU3zzOfZ5BSWMzgl\nsU3PUgpkXKSqxs8f398EwI9PG0B8THiuxDeRI5Ak8hPgHhHZJSK7gd8AN4U2LGMCk3fwCI9/tBWA\n+y8YSWx0eC8sbEgg60We/nwHW/aXkZHUiautWq8JA43WSFDV7cAkd5ovqloW8qiMCYCqct/c9VRU\n+5k2qg+ntPHqtePdJLJ8ZzFVNf5vrbTPPnCYRxc4CfPBi0aFbT0wE1kCKrQjItOAkUCH2qmEqvpA\nCOMyplHvrN3Hhxv2kxgfw2/Pb1tTeuuTnBjPkJREtuaX8eTH2/jFd4Z+dU5VuffNtVTV+Ll4bL82\nnzBN+xFI7ayngMuBnwECXAZYO9p46mB5Fb+buw6A35w7jD5d28dOfvdOG06UwKMLtjJvjbNDtary\n7MIcFm15cY4xAAAUfklEQVQvpHunWH47bYTHURrztUBaIier6mgRWaOq94vI/2EFGI3HHpy3kQNl\nVUzI6MFV7Wg/8TOOS+Huc4fz4DsbuePVVfhVeXHJThbvcAbbfzttBD0SbHKkCR+BJJEK989yEemL\ns7OhrW4ynvli6wFeXZ5LXHQUf7hkVLvbhOnHpw1g8/5SXluey8/+vRKAHglx3HXuMC62Kb0mzASS\nRN4SkW7Aw8AKQHFWrRvT6kqOVPPr11YDcNvZgxnUDivXiggPXnQ8u4rKycop4upJ6fzyO0Pp1sla\nICb8NJhERCQKWKCqB4HXReRtoIOqlrRKdMYc5b6569lTUsEJqV25afIgr8MJmfiYaGb/eCKlFTV0\nt+4rE8YaHFhXVT/wZJ3XlcFIICJymYisFxG/iGTWOf4dEVkuImvdP8+qc26ce3ybiDwmVnEu4ry9\nZg9vrsyjQ2wUf778xDa9JiQQMdFRlkBM2AvkX+ECEbkkyL+01wEX45SXr+sA8D1VHQX8EHi+zrm/\nATcCQ9yvqUGMx4S5fSUV3PumMxvr3mkjbAMmY8JEIEnkJpyCi5UickhESkXkUEseqqobVXVzPcdX\nquoe9+V6oKOIxItIH6CLqi5WVQWeAy5sSQym7fD5lTteXUXJkWomD+3J1RPbz2wsY9q6QFase1X2\n/RJghapWikg/ILfOuVzApqlEiMc/2srCbYUkJcTx8KWjbe8MY8JIIIsNFwRyrJ5r5ovIunq+pgfw\n3pHAQzSzRpeIzBCRLBHJKigoaM4tTJhYtO0Ajy7Yigg8esUYUtpghV5j2rNjtkREpAPQCUgWke44\nq9UBuhBAK0BVpzQnIBFJBd4ErnXrdgHkAal1Lkt1jx3r2TOBmQCZmZnanDiM9/JLK7jtpVWowm1n\nD+HUIVbqw5hw01B31k3A7UBfYDlfJ5FDwBOhCMZdjzIPuEtVF9YeV9W97njMJGAJcC3weChiMOGh\nqsbPrbNXcqCskpMHJfHzs4d4HZIxph7H7M5S1UdVdQBwp6oOVNUB7tcJqtqiJCIiF4lILnASME9E\n3ndP3QoMBv5bRFa5XynuuZuBZ4BtOPu8W+mVduyBt9ezNLuIXl3i+csVJ9rufcaEKXEmOzVykcjJ\nQAZ1Wi6q+lzowgqezMxMzcrK8joM0wQvLtnJvW+uIy4milduOokT+3fzOiRjIo6ILFfVzMaua3R2\nlog8DwwCVuHsrw5O6ZM2kURM27I0u4jfzXG2uv3DRaMsgRgT5gKpnZUJjNBAmizGtMCOgjJuej6L\nGr/y41MHcMm41MbfZIzxVCCLDdcBvUMdiIlsB8oque6fyygur+asYSncde4wr0MyxgQgkJZIMrBB\nRJYClbUHVfWCkEVlIkp5VQ03zFrGrqJyRvXryhM/GENMO6+LZUx7EUgSuS/UQZjIVe1zpvKuzi0h\ntXtH/nFdJp3iAtq12RgTBgIpe/JpawRiIo/Pr/zi5VV8tCmfbp1imXX9BFI624p0Y9qShlasl+LM\nwvrWKUBVtUvIojLtnt+v3PPGWt5es5fE+Bie+9EEBqdYZV5j2ppjJhEPCy+adk5VeeDtDbyctZsO\nsVE8e914RqfaVF5j2iLrfDatyu9X7ntrPc99uZPYaOHv12QyYUAPr8MyxjSTJRHTavx+5d7/rOPf\nS3cRFx3FU9eMZfLQnl6HZYxpAUsiplXU+Pzc9cZaXlueS3xMFDOvzbQEYkw7YEnEhFxFtY9bZ69k\n/sb9dIiN4h8/HM8pg62suzHtgSURE1IlR6q58V9ZLM0pomvHWJ69bjzj0rt7HZYxJkgsiZiQyS0u\n54ZZWWzeX0rvLh147oYJDO1lk/6MaU8siZiQWLGrmBnPZXGgrIpBPRN47oaJ9OvW0euwjDFBZknE\nBN1bq/dw56urqazxc8rgJP561Ti6doz1OixjTAhYEjFBU+Pz8/AHm/n7pzsAuHJCGg9MH0msFVM0\npt2yJGKCouhwFbf9eyVfbDtAdJRw73nDuf6UDERsW1tj2jNLIqbFlu8s5rZ/ryTv4BGSE+N44gdj\nmTQwyeuwjDGtwJKIaTa/X/n7Zzv40web8fmVE/t3429Xj6VPVxtANyZSWBIxzbL/UAW/em0Nn20p\nAGDG6QO587vHERdj4x/GRBJP/sWLyGUisl5E/CKSWc/5NBEpE5E76xwbJyJrRWSbiDwm1tnumXlr\n9nLOXz7jsy0FdO8Uy7PXZXLPecMtgRgTgbxqiawDLgb+fozzjwDvHnXsb8CNwBLgHWBqPdeYECo+\nXMV9b61nzqo9AEwe2pM/XjqaXl1sIyljIpUnSURVNwL1ztwRkQuBbOBwnWN9gC6quth9/RxwIZZE\nWoWqMm/tXn43Zz2Fh6voGBvNvdOGc9XENJt9ZUyEC6sxERFJBH4DfAe4s86pfkBunde57rFj3WcG\nMAMgLS0t+IFGkLyDR7hv7no+3LAfgIkDevDQJaPJSE7wODJjTDgIWRIRkflA73pO3auqc47xtvuA\nP6tqWUs+4arqTGAmQGZmZn1b/JpGVNX4eXZhNo/O38qRah+J8THcc95wrhjfn6goa30YYxwhSyKq\nOqUZb5sIXCoifwS6AX4RqQBeB1LrXJcK5LU8SlOfz7YU8MDbG9iWXwbAtFF9+K/zR9C7q419GGO+\nKay6s1T1tNrvReQ+oExVn3BfHxKRSTgD69cCj3sSZDu2o6CMB+dtZMGmfAAykjpx//TjbfMoY8wx\neZJEROQinCTQE5gnIqtU9ZxG3nYzMAvoiDOgboPqQXKgrJLHFmxl9pJd1PiVxPgYbj1rMNefkkF8\nTLTX4RljwphXs7PeBN5s5Jr7jnqdBRwfwrAiTmlFNf/4IpunP9vB4SofUQKXZ/bnjnOGktLZuq6M\nMY0Lq+4s0zoOV9Ywa1EOT3++g4Pl1QBMGZ7Cr84ZxnG9bdMoY0zgLIlEkEMV1Tz/5U7+8UU2RYer\nAJiQ0YM7zzmOCQN6eBydMaYtsiQSAQrLKpm1KIdZi3IoragBYExaN+74znGcMjjJFgwaY5rNkkg7\nln3gMM98voPXludSWeMHnMWCt541mFMHJ1vyMMa0mCWRdkZVWbitkFmLslmwKR91l1qePSyFn5wx\niPEZ1m1ljAkeSyLtxKGKav6zMo/nv9zJVneRYFx0FBeO6cuNpw1kSC8bMDfGBJ8lkTZMVVmbV8K/\nl+5izqo9lFf5AEjpHM81k9K5cmIayYnxHkdpjGnPLIm0QcWHq5izKo+Xlu1m077Sr46fNDCJqyal\n8d0RvW1vD2NMq7Ak0kZU1vj4ZHMBry/P5ePN+VT7nMGOHglxXDSmH1dOSGNwSqLHURpjIo0lkTDm\n8ytLs4uYuzqPd9buo+SIszAwSpwNob6f2Z8pI1KsNIkxxjOWRMKM369k7Sxm3po9vLtuH/mllV+d\nG9a7MxeP7ceFJ/YjxXYTNMaEAUsiYaCqxs/iHYW8v34fH2zYT0GdxNG/R0cuOKEv00/sx1CbYWWM\nCTOWRDxSUl7NJ1vymb8xn08253+1khwgtXtHpo3qw7TRfRjVr6stCjTGhC1LIq1EVdm4t5RPtuTz\nyaYClu8qxuf/etPFISmJTD2+N+eM7M3Ivl0scRhj2gRLIiFUUFrJou0H+HzrAT7bUvCN8Y2YKOGk\ngUlMGdGLKcNTSE+yPcuNMW2PJZEgOlRRzbLsIhZtL2ThtgPfWMMB0KtLPKcP6clZw1I4dUgynTvE\nehSpMcYEhyWRFig5Uk1WThFLsotYsqOQtXkl1OmhIj4migkDenDq4GQmH9eT43p1tm4qY0y7Ykkk\nQKpK3sEjLN9ZTFZOMctyiti8v/SrAofgdFGd2L8rJw9K5qRBSYxL706HWFvDYYxpvyyJHMORKh/r\n9pSwclcxK3cdZMWuYvYfqvzGNXHRUYxO7crEgT2YOCCJsendSYy3v1JjTOSw33j12Lq/lKmPfv6N\n2VMAXTvGMi69O+PSuzM+owejU7taS8MYE9E8SSIichlwHzAcmKCqWXXOjQb+DnQB/MB4Va0QkXHA\nLKAj8A7wc1VVQiAjOYEOMVGkJSUwJq0bY/p3Y0xaNwYmJxIVZWMaxhhTy6uWyDrgYpxk8RURiQFe\nAK5R1dUikgRUu6f/BtwILMFJIlOBd0MRXGx0FCv++ztWk8oYYxrhSb1wVd2oqpvrOfVdYI2qrnav\nK1RVn4j0Abqo6mK39fEccGEoY7QEYowxjQu3TSeGAioi74vIChH5tXu8H5Bb57pc95gxxhgPhaw7\nS0TmA73rOXWvqs5pIJ5TgfFAObBARJYDJU189gxgBkBaWlpT3mqMMaYJQpZEVHVKM96WC3ymqgcA\nROQdYCzOOElqnetSgbwGnj0TmAmQmZkZksF3Y4wx4ded9T4wSkQ6uYPsk4ENqroXOCQik8RZ8n0t\ncKzWjDHGmFbiSRIRkYtEJBc4CZgnIu8DqGox8AiwDFgFrFDVee7bbgaeAbYB2wnRzCxjjDGBkxAt\ntQgbmZmZmpWV1fiFxhhjviIiy1U1s7Hrwq07yxhjTBvS7lsiIlIA7PQ6jiZKBg54HUQrs585MtjP\n3Hakq2rPxi5q90mkLRKRrECake2J/cyRwX7m9se6s4wxxjSbJRFjjDHNZkkkPM30OgAP2M8cGexn\nbmdsTMQYY0yzWUvEGGNMs1kSCXMicoeIqIgkex1LqInIwyKySUTWiMibItLN65hCRUSmishmEdkm\nInd5HU+oiUh/EflYRDaIyHoR+bnXMbUGEYkWkZUi8rbXsYSKJZEwJiL9cfZY2eV1LK3kQ+B4VR0N\nbAHu9jiekBCRaOBJ4FxgBHCliIzwNqqQqwHuUNURwCTglgj4mQF+Dmz0OohQsiQS3v4M/BqIiIEr\nVf1AVWvcl4v5ZuXm9mQCsE1Vd6hqFfASMN3jmEJKVfeq6gr3+1KcX6ztek8gEUkFpuHU/Gu3LImE\nKRGZDuTV7vIYgX5E+y2y2Q/YXed1RG2yJiIZwBicra7bs7/gfAj0ex1IKHm1x7qh4Y27gHtwurLa\nlUA2KxORe3G6P15szdhM6IlIIvA6cLuqHvI6nlARkfOBfFVdLiJneB1PKFkS8dCxNu4SkVHAAGC1\ns30KqcAKEZmgqvtaMcSga2yzMhG5DjgfOFvb7/zzPKB/ndcNbrLWXohILE4CeVFV3/A6nhA7BbhA\nRM4DOgBdROQFVb3a47iCztaJtAEikgNk1u742F6JyFSc/WQmq2qB1/GEirvh2hbgbJzksQz4gaqu\n9zSwEHI3k/sXUKSqt3sdT2tyWyJ3qur5XscSCjYmYsLJE0Bn4EMRWSUiT3kdUCi4kwduxdnJcyPw\nSntOIK5TgGuAs9z/tqvcT+mmjbOWiDHGmGazlogxxphmsyRijDGm2SyJGGOMaTZLIsYYY5rNkogx\nxphmsyRiPCEiPnea53oRWe1WK673/0cR6SsirwXhmRkicqTOFNNVIhLX0vuGCxG5TkQK3KqxW0Xk\nfRE5OYD3XdjSYogi8k57rrpsjs1WrBuvHFHVEwFEJAWYDXQBflf3IhGJUdU9wKVBeu722ufWx31e\nzbHOtwEvq+qtACJyJvCGiJypqg1Vkr0QeBvY0NyHqqqt+YhQ1hIxnlPVfGAGcKs4rhORuSLyEbDA\nbUGsAxCRxSIysva9IvKJiGSKSIKIPCsiS91P4gFXxRWR+0TkeRFZCDzv7gHxsIgsc/c2ucm9TkTk\nCXcfkPnup+9L3XM5tXu+uPF84n5fb1zuz/iGiLznthr+WCeeqSKywm2hLRCRKPeanu75KHcfkp6N\n/L1+jLM16wz3fYPc5y0Xkc9FZJjbUrkAeNhtmQ2q/Tt135PsVkxoLOYc99oMEdkoIk+7rcwPRKSj\ne8149+9zlfv3uy7Q/0YmfFkSMWFBVXcA0UCKe2gscKmqTj7q0peB7wOISB+gj6pm4RSt/EhVJwBn\n4vxSTKjnUYPqdGU9Wef4CGCKql4J3ACUqOp4YDxwo4gMAC4CjnOvvRZotKuokbhOBC4HRgGXi7Nx\nU0/gaeASVT0BuExV/cALwFXu+6YAqwMsDbMCGOZ+PxP4maqOA+4E/qqqi4C5wK9U9URV3d7I/b4V\ncz3XDAGeVNWRwEHgEvf4P4Gb3JagL4DYTRtg3VkmXH2oqkX1HH8F+ACn2+v7QO1YyXdxCt7d6b7u\nAKTx7Q2BjtWdNVdVj9S51+jaVgbQFecX4+nAv1XVB+xxW0qNOVZcAAtUtQRARDYA6UB34DNVzQao\n83fwLDAHp7z4j3B+IQdC3Psn4iS9V8Up6gkQH+A96qov5t1HXZOtqqvc75cDGe54SWdV/dI9Phun\n0KZp4yyJmLAgIgNxPp3mu4cO13edquaJSKGIjMb5RPyT2lvgfHrf3MwQ6j5PcD6xv39UjA31+9fw\ndcu+w1H3+lZcIjIRqKxzyEcD/x5VdbeI7BeRs3A2tbrqWNceZQxOIo0CDjY0HlTHsX4WAoz56Gs6\nBhaqaYusO8t4zu3CeQp4IsDy7y/jbPbTVVXXuMfeB34m7sdsERnTgpDeB34qTulyRGSo2wX1GU4X\nTrTblXZmnffkAOPc7y856l5NiWsxcLrbfYaI9Khz7hmcbq1X3dZQg0RkMs54yNPu3h3ZInKZe05E\n5AT30lKcwpf1/SxBmdCgqgeBUjd5AlwRjPsa71kSMV7p6I5LrAfm43RR3R/ge1/D+SX0Sp1jvwdi\ngTXuPX/fgtiewZmptMId/P07zifuN4Gt7rnngC/rvOd+4FERyeKb/f1Nissd55iBM6tqNU7CrDUX\nSKThrqzL3b/XLTgbm11SZ2bWVcAN7n3X8/WWvC8Bv3IH/gcBf8JJoiuB5IbibaIbgKdFZBWQAJQE\n8d7GI1bF15hmEpFZwNuq2uI1LAE+LxP4s6qe1hrPCzYRSVTVMvf7u3AmRfzc47BMC9mYiDFtgPtL\n96cEPhYSjqaJyN04v3d2Atd5G44JBmuJGGOMaTYbEzHGGNNslkSMMcY0myURY4wxzWZJxBhjTLNZ\nEjHGGNNslkSMMcY02/8HGDmcuYRZxtQAAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fig = plt.figure(1)\n", "ax = fig.add_subplot(111)\n", "ax.plot(wlist, count1, wlist, count2, lw=2)\n", "ax.set_xlabel('Drive Frequency Detuning')\n", "ax.set_ylabel('Count rates')\n", "\n", "# plot phase shift of cavity light\n", "fig2 = plt.figure(2)\n", "ax2 = fig2.add_subplot(111)\n", "ax2.plot(wlist, 180.0 * np.angle(infield) / np.pi, lw=2)\n", "ax2.set_xlabel('Drive Frequency Detuning')\n", "ax2.set_ylabel('Intracavity phase shift');" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
SoftwareVersion
QuTiP4.3.0.dev0+6e5b1d43
Numpy1.13.1
SciPy0.19.1
matplotlib2.0.2
Cython0.25.2
Number of CPUs2
BLAS InfoINTEL MKL
IPython6.1.0
Python3.6.2 |Anaconda custom (x86_64)| (default, Jul 20 2017, 13:14:59) \n", "[GCC 4.2.1 Compatible Apple LLVM 6.0 (clang-600.0.57)]
OSposix [darwin]
Thu Jul 20 22:07:48 2017 MDT
" ], "text/plain": [ "" ] }, "execution_count": 9, "metadata": {}, "output_type": "execute_result" } ], "source": [ "from qutip.ipynbtools import version_table\n", "\n", "version_table()" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.6.2" } }, "nbformat": 4, "nbformat_minor": 1 }