{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# QuTiP example: Dynamics of an atom-cavity system using three different solvers" ] }, { "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\n", "import matplotlib.pyplot as plt" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "import numpy as np" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "from qutip import *" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Model and parameters" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "kappa = 2\n", "gamma = 0.2\n", "g = 1 \n", "wc = 0 \n", "w0 = 0 \n", "wl = 0\n", "N = 4\n", "E = 0.5\n", "tlist = np.linspace(0, 10, 200)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### mesolve" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [], "source": [ "def solve(E,kappa,gamma,g,wc,w0,wl,N,tlist):\n", "\n", " ida = qeye(N)\n", " idatom = qeye(2)\n", "\n", " # Define cavity field and atomic operators\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 + 1j*g*(a.dag()*sm - sm.dag()*a) \\\n", " + 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", " #intial state\n", " psi0 = tensor(basis(N,0),basis(2,1))\n", " rho0 = psi0.dag() * psi0;\n", "\n", " # evolve and calculate expectation values\n", " output = mesolve(H, psi0, tlist, [C1, C2], [C1dC1, C2dC2, a]) \n", "\n", " return output.expect[0], output.expect[1], output.expect[2]" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [], "source": [ "count1, count2, infield = solve(E,kappa,gamma,g,wc,w0,wl,N,tlist)" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAtoAAAF3CAYAAACbhOyeAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3Xl4XHXZ//H3ncneLN3SNU3TlbayFCgt0gIPIAqolAcU\nEUVFsSAibs/PHX1Wt8cHBQURBBU3XABBBQFBNoHSha2lLaR7uqVtmibNvty/P2bSTtMsJ2kmJ5n5\nvK7rXHPOmXNmPilccPeb77m/5u6IiIiIiEj/Sgs7gIiIiIhIMlKhLSIiIiKSACq0RUREREQSQIW2\niIiIiEgCqNAWEREREUkAFdoiIiIiIgmgQltEREREJAFUaIuIiIiIJIAKbRERERGRBFChLSIiIiKS\nAOlhB+hPo0eP9tLS0rBjiIiIiEgSW7FixR53L+rpuqQqtEtLS1m+fHnYMUREREQkiZnZ5iDXaeqI\niIiIiEgCqNAWEREREUmAhBbaZnaema0zszIz+1In788ys+fNrNHM/q2T9yNm9pKZ/SWROUVERERE\n+lvCCm0ziwC3AOcDc4D3m9mcDpdVAtcD3+viYz4NrElURhERERGRREnkiPZ8oMzdN7h7E3APsDj+\nAnevcPdlQHPHm82sGHgn8NMEZhQRERERSYhEFtoTga1xx+Wxc0H9APgC0NafoUREREREBsKgfBjS\nzN4FVLj7igDXLjGz5Wa2fPfu3QOQTkRERESkZ4kstLcBk+KOi2PnglgIXGhmm4hOOTnbzH7V2YXu\nfru7z3P3eUVFPfYNFxEREREZEIkstJcBM8xsipllApcBDwa50d2/7O7F7l4au+8Jd/9g4qKKiIiI\niPSvhK0M6e4tZnYd8AgQAe5y99Vmdk3s/dvMbBywHCgA2szsM8Acd69OVC4RERERkYFg7h52hn4z\nb9481xLsIiIiIpJIZrbC3ef1dF3CRrQlee050MiaHdVkZ0TIyYgwpiCLMfnZYccSERERGVRUaEsg\n7s7SjZX86oXNPLJ6J82th34Tkmbw9jnjuOr0KZw8eQRmFmJSERERkcFBhbb0qL6plet+s5LH11ZQ\nkJ3OB0+dzLlzxtLa5tQ1tfLy1ip+s3QLf1u9k5Mnj+D7l86lZFRu2LFFREREQqU52tKtqromPvaL\n5azcso8vnTeLD721lJzMyBHX1TW1cO+Kcr736BsA/PD9J3LGTLVbFBERkeQTdI72oFywRgaHnfsb\nuPQnz/Na+X5uvfwkrj5zWqdFNkBuZjpXvLWUB69byPjCbD7ysxe57an1JNNf5ERERER6Q4W2dKqp\npY2P/nwZ26sa+PlHT+H848YHum/yqGHcd+1pnH/ceL798FpufXJ9gpOKiIiIDE6aoy2duvnxN3l9\nRzV3fGgep00b3at7czPT+eFlJ5KRZvzvI+soys/i0nmTer5RREREJImo0JYjvLRlH7c+WcZ7Ti7m\n3Dlj+/QZaWnGd99zAntrm/jyfa8xalgm58zu22eJiIiIDEWaOiKHqW9q5fO/f4XxhTl8/d1zjuqz\nMtPTuO2DJ/OWCQV88jcreX27FvwUERGR1KFCWw5z42Pr2LCnlu++53gKsjOO+vOGZaVz10dOoSA7\ng0/9diV1TS39kFJERERk8FOhLQftqm7gF89t5tJ5xSyc3rt52d0ZnZfF9983lw17avnPP7/eb58r\nIiIiMpip0JaDbn96A63uXHfWjH7/7IXTR/OJM6dxz7Kt/OXV7f3++SIiIiKDjQptAWDPgUZ+vXQz\nF82dmLBVHT977kxOLBnOl+97jW1V9Qn5DhEREZHBQoW2AHDnsxtpbGnj2rOmJew7MiJp3HzZibS0\nOt94YJUWsxEREZGkpkJbqKpr4u7nNvHO48YzrSgvod81aWQunzt3Jn9fU8Ejq3cm9LtEREREwqRC\nW/j5c5uobWrlurOnD8j3XbmwlDnjC/jGg6upbmgekO8UERERGWgqtFNcc2sbv3phM2+bPYZZ4woG\n5DvTI2l86+LjqKhp5P8eWTcg3ykiIiIy0FRop7in1u1mz4Em3ndKyYB+7wmThvPht5Zy9wubebW8\nakC/W0RERGQgqNBOcfeuLGfUsEz+5ZiiAf/uz799JqOGZfJff3ldD0aKiIhI0lGhncKq6pp4fE0F\nF86dQEZk4P9VyM/O4HPnHsOyTfv42yo9GCkiIiLJRYV2Cvvzqztoam3jkpOKQ8tw6bxijhmbz7f/\ntpbGltbQcoiIiIj0NxXaKezeFeXMGpfPWyYMzEOQnUmPpPHVd85m8946fvn85tByiIiIiPQ3Fdop\nav3uA7y8tYpLTirGzELNcsbMIv7lmCJuevxNKmubQs0iIiIi0l9UaKeoe1eUE0kzFp84IewoAHzl\ngtnUNrbw4yfLwo4iIiIi0i9UaKcgd+eBl7dz+ozRjMnPDjsOADPH5vOvJxZz9/Ob2VXdEHYcERER\nkaOmQjsFrdtVw7aqes57y7iwoxzm0+fMoLXN+eETb4YdRUREROSoqdBOQU+srQDgrFljQk5yuJJR\nubzvlEnc8+JWtlbWhR1HRERE5Kio0E5BT6yp4LiJhYwtGBzTRuJ96uwZRNKMmx7XqLaIiIgMbYEK\nbYsaY2YT2rdEB5PEqKxtYuWWfZw9yEaz240rzOaKUydz38pyyioOhB1HREREpM96LLTN7FpgN/AM\n8Hhs+3uCc0mCPPVGBW0O58wenIU2wCf+ZRrZGRF+pLnaIiIiMoQFGdH+HDDb3Y9x99mxbU6ig0li\nPL6mgqL8LI6dUBh2lC6Nysvig6dO5sFXtrNxT23YcURERET6JEihXQ5UJjqIJF5zaxtPvbGbs48Z\nQ1pauIvU9OSq06eQEUnjln+or7aIiIgMTekBrikDnjCzvwCN7Sfd/eaEpZKEWL5pHzUNLZw9iKeN\ntBuTn83755fwyxc28+lzZjBpZG7YkURERER6JciI9g7gaaAAKIrbZIh5Yu0uMiNpLJo+OuwogVxz\n5jQiZtz65Pqwo4iIiIj0Wo+Ftrvf4O43AN8Evhl33CMzO8/M1plZmZl9qZP3Z5nZ82bWaGb/Fnd+\nkpn9w8xeN7PVZvbp3vxQ0rkn1lZw6rRRDMsK8ouM8I0rzObSU4r544qtbK+qDzuOiIiISK8E6Toy\nx8yWAW8Cb5rZUjObHeC+CHALcD4wB3i/mXV8iLISuB74XofzLcDnYw9dngp8spN7pRcqqhtYv7uW\n04fIaHa7a86chjvc/vSGsKOIiIiI9EqQqSO3A19x92J3Lwa+CtwR4L75QJm7b3D3JuAeYHH8Be5e\n4e7LgOYO53e4+8rYfg2wBpgY4DulCy9sjD7PumDqyJCT9E7xiFwunDuB3y3byr7aprDjiIiIiAQW\npNDOd/fH2g/c/e9AfoD7JgJb447L6UOxbGalwInA0t7eK4cs3bCXvKx05owvCDtKr119xjTqm1u5\n+/nNYUcRERERCSxIob3JzL5sZsWx7UvApgTnAsDM8oB7gc+4e3UX1ywxs+Vmtnz37t0DEWtIenFj\nJfNKR5AeCbQY6KByzLh8zp41hl88v4n6ptaw44iIiIgEEqTq+igwCXgotk2KnevJtti17Ypj5wIx\nswyiRfav3f2+rq5z99vdfZ67zysqUjOUzuw90MibFQdYMGVU2FH67Jozp1FZ28QfVmzt+WIRERGR\nQaDH9hPuvhe4tg+fvQyYYWZTiBbYlwGXB7nRzAy4E1jj7jf24bslzoux+dnzpwyt+dnxTikdwYkl\nw7njmQ1cPr9kSI7Mi4iISGrpsloxs/+Lvd5vZvd13Hr6YHdvAa4DHiH6MOPv3X21mV1jZtfEPnuc\nmZUTXeb9a2ZWbmYFwELgCuBsM3s5tl1w1D9tilq6sZKcjAjHFw/eZdd7YmZcc+Y0tlbW89CqnWHH\nEREREelRdyPav4u9/qivH+7u7dNN4s/dFre/k+iUko6eBQb3GuFDyAsb9nLy5BFkDPFR4HNnj2Vq\n0TB+8tR63n38eKK/+BAREREZnLqsvNz9xdjr4+0bsBzYEduXIaCqrol1u2pYMISnjbRLSzOuPmMq\nq7dX82zZnrDjiIiIiHQryII1j5tZgZmNAF4Gfmlm/5v4aNIflm3ahzssmDp0H4SMd9GJExmTn8VP\nntICNiIiIjK4BZlLMDLWWu9i4FfufjLwjsTGkv6ydMNeMtPThvT87HhZ6RE+umgKz5bt4bXy/WHH\nEREREelSkEI73cyKgPcCf05wHulnSzdWcuKk4WRnRMKO0m8uX1BCflY6P3l6fdhRRERERLoUpND+\nH+ApYIu7v2hmU4GNiY0l/aGuqYXV2/cP6bZ+nSnIzuDyU0t46LUdbNlbF3YcERERkU71WGi7+z3u\nPsfdl8SON7j74sRHk6O1als1bQ5zJw0PO0q/++jCKaSnpXHHM5qrLSIiIoNTkIchvxV7GDLdzB4x\ns11mFmjhGQnXq+VVABxfnHyF9tiCbBbPncAfVmxlX21T2HFEREREjhBk6sj5sYch3wVsB2YDX0xo\nKukXL2+tYuLwHIrys8KOkhAfP2MqDc1t/Hrp5rCjiIiIiBwh0MOQsdcLgD+4eyXgiYsk/eXV8v2c\nMCk5uo10ZubYfM6cWcTPn9tMQ3Nr2HFEREREDhOk0H7YzFYBC4DHzGw00JjYWHK0Kmub2FJZl5TT\nRuItOWMqew408uDL28OOIiIiInKYIA9D/j/gbOBkd28GGoj21JZBrH1+9glJXmifNm0Us8cXcMcz\nG3DXL1pERERk8Oiy0DazM2OvFwKnAhfE9s8GTh6YeNJXr2zdjxkclyQL1XTFzFhyxhTerDjAk2/s\nDjuOiIiIyEHdjWifG3t9byfbexKcS47Sq+VVTC/KIy8rveeLh7h3HT+BcQXZ3PG0Wv2JiIjI4NFl\nFebuX4u9XjFwcaQ/uDuvlFdx5swxYUcZEBmRNK5cWMq3Hl7Lqm37OXZico/ii4iIyNAQpI92gZld\na2bfNbMb27eBCCd9s31/A3sONCV1x5GOLptfwrDMCD/VAjYiIiIySATpOvIQMAt4E1gdt8kg9crW\n1HgQMl5hTgaXzS/hL6/uYHtVfdhxRERERAIV2rnufr273+Hud7ZvCU8mffZKeRUZEWPW+Pywowyo\nKxeW4sDPn9sUdhQRERGRQIX2b8zsSjMrik0jKTCzgoQnkz57ZWsVc8YXkJUeCTvKgCoekcv5x47j\nt0u3UNPQHHYcERERSXFBCu0DwA+Alzg0bWRVIkNJ37W1Oau2VSf9QjVdWXLGVGoaW/jdsq1hRxER\nEZEUF6TQ/iIww92L3X1SbCtJdDDpm82VdRxobOG4FO28cXzxcOZPGcnP/rmJlta2sOOIiIhICgtS\naJcB1YkOIv1j7Y7oP6rZ41N3ds+S06eyraqeh1btDDuKiIiIpLAgq5lUAy+Z2RNAY/tJd/9cwlJJ\nn63ZUU2awYyxeWFHCc3Zs8YwtWgYdzy9gXcfPx4zCzuSiIiIpKCg7f2+B6xE7f0GvTU7a5gyehjZ\nGan1IGS8tDTjqkVTeW3bfpZurAw7joiIiKSoLke0zWyYu9d21srPzCYmNpb01dqd1SnVP7srF580\nkf97dB13PL2BU6eOCjuOiIiIpKDuRrSfad8xs0c7vPfnxMSRo1HT0MzWyvqUnp/dLjsjwhVvnczj\naysoqzgQdhwRERFJQd0V2vETW4u6eU8GiXU7awCYnWIL1XTlilMnk5Wexp3Pall2ERERGXjdFdre\nxX5nxzIIrIl1HJk1TiPaAKPysrjk5GLuXbmNPQcae75BREREpB91V2iPMbPrzezTcfvtxx1HuGUQ\nWLOzhsKcDMYXZocdZdD42KIpNLW0cffzm8OOIiIiIimmu0L7Z0QL6tFx++3HP094Mum1NTuqmTUu\nX+3s4kwryuNts8fyqxc209DcGnYcERERSSFddh1x9xsGMogcnbY2Z93OGi6dNynsKIPOx0+fwvtu\n38W9K8v5wILJYccRERGRFBGkj7YMAVv31VHX1KoHITsxf8pITigu5M5nNtLWpscLREREZGCo0E4S\na7T0epfMjKtOn8qGPbU8vrYi7DgiIiKSIlRoJ4k1O2qiS6+P0Yh2Z84/dhwTh+dwx9Nq9SciIiID\no8dC28yuM7OC2P5PzOxFMzsn8dGkN9bsqKZ09DByMlN36fXupEfS+OiiKby4qZKXt1aFHUdERERS\nQJAR7SXuXm1mbwfGAh8HvpvYWNJba3fWaNpID953yiTys9O54xmNaouIiEjiBSm0258euwD4pbu/\nEvA+zOw8M1tnZmVm9qVO3p9lZs+bWaOZ/Vtv7pVDDjS2sKWyjtnjNG2kO3lZ6Vy+oISHX9vB1sq6\nsOOIiIhIkgtSML9iZg8B7wIeNrM8AqwMaWYR4BbgfGAO8H4zm9PhskrgeuB7fbhXYsoqDgAwc6wK\n7Z585LRS0sy4658bw44iIiIiSS5IoX0l8O/AfHevA7KBjwW4bz5Q5u4b3L0JuAdYHH+Bu1e4+zKg\nubf3yiHthfb0MXkhJxn8xhfmcOEJE/j9sq3sr+/4r52IiIhI/wlSaC8gurDNLDM7DZgJ5AS4byKw\nNe64PHYuiKO5N+W8WVFDZiSNkpG5YUcZEq46fSq1Ta389sUtYUcRERGRJNblypBx4leIzAZOBl4C\nzkxIol4ysyXAEoCSkpKQ04RjfcUBpoweRnpE3RqDmDOhgEXTR3PXsxu5cmEpWenq1CIiIiL9r8fK\nzN3Pj9vOAo4Hgqz6sQ2IXw+8OHYuiMD3uvvt7j7P3ecVFRUF/PjkUlZxQNNGeumaM6dRUdPIfSuD\n/ispIiIi0ju9HgJ1903AWwJcugyYYWZTzCwTuAx4MODXHM29KaWhuZUtlXVMU6HdKwunj+KE4kJu\ne2o9La1tYccRERGRJNTj1BEz+z6HuoykAScCr/R0n7u3mNl1wCNABLjL3Veb2TWx928zs3HAcqAA\naDOzzwBzYn27j7i39z9e8tu4p5Y214OQvWVmXHvWdK7+5Qr++toOFs/VIwAiIiLSv4LM0V4Vt98C\n3O/uTwX5cHd/CHiow7nb4vZ3Ep0WEuheOVJ7x5EZKrR77dzZY5kxJo9b/7Gedx8/gbQ0CzuSiIiI\nJJEgc7TvBH4B/BNYGnuVQaKs4gBpBlNGDws7ypCTlmZce9Y01u2q4Ym1QR47EBEREQmux0LbzE4H\nyoA7gbuAN8xsYaKDSTBlFQeYNDKX7Ax1zuiLdx8/geIROfzoH2W497gOk4iIiEhgQR6G/D5wgbsv\ndPfTgHcCNyU2lgRVVnGA6UWaNtJX6ZE0rjlzGi9vreL5DXvDjiMiIiJJJEihnenur7cfuPsaIDNx\nkSSoltY2Nu6pZfpYFdpH4z0nF1OUn8Wt/1gfdhQRERFJIkEK7ZVmdpuZLYptPya6YI2EbEtlHU2t\nbRrRPkrZGRGuWjSFZ8v28MrWqrDjiIiISJIIUmhfA2wAvhDbNgBXJzKUBNPecUSt/Y7eB06dTEF2\nOrc+WRZ2FBEREUkSPbb3c/cGM/sBcJ+7qwoZRMp2RwttLVZz9PKy0vnIwinc/PibvLmrhhlj88OO\nJCIiIkNckK4j7wJeAx6LHc81s/sTHUx6VlZxgLEFWRRkZ4QdJSlceVopuZkRfvyk5mqLiIjI0Qsy\ndeQ/gAVAFYC7vwxMT2QoCaas4gAzxmjktb+MGJbJ5fNLeOCV7WzaUxt2HBERERnighTaze7e8Qkx\nNRwOmbuzvuKA5mf3syVnTiUjYtz8xJthRxEREZEhLkihvcbMLgXSzGyKmX0feCHBuaQHO/Y3UNvU\nqvnZ/WxMfjYfemspf3pp28GHTUVERET6IkihfR1wMtAG3A80AZ9JZCjp2cbY1IZpWnq93119xlSy\nMyLc/LhGtUVERKTveiy03b3W3b/o7ie6+9zYft1AhJOubYgV2lOKVGj3t1F5WXz4tFL+/Op23thV\nE3YcERERGaKCdB2Zbma3mtlDZvZo+zYQ4aRrm/bUkp2Rxtj87LCjJKUlp08lNyPCTX/XqLaIiIj0\nTY99tIE/AncCvwJaExtHgtq0p5bSUcNIS7OwoySlEcMy+eiiKfzwiTKu21HN7PEFYUcSERGRISbI\nHO02d/+huz/n7kvbt4Qnk25t3FPLFM3PTqirFk0lPzudH/z9jbCjiIiIyBAUpNB+wMyWmFmRmRW0\nbwlPJl1qaW1jS2UdpSq0E6owN4OrFk3lkdW7WLVtf9hxREREZIgJUmhfBdwArARWx7ZViQwl3dtW\nVU9Lm2tEewBcuaiUwpwMjWqLiIhIrwXpOjKpk61kIMJJ5w52HFGhnXAF2RksOWMqf19TwStbO67b\nJCIiItK1ICPamNksM7vYzC5v3xIdTLrWvjx46SgV2gPhw6eVMiI3g+9rVFtERER6IUh7v68BtwO3\nAecDPwDek+Bc0o2Ne2rJz0pndF5m2FFSQl5WOlefOY0n1+3mxY2VYccRERGRISLIiPb7gLOAHe5+\nBXACoKHUEG3cU0vp6GGYqbXfQPnwW0sZk5/Ftx9eg7uHHUdERESGgCCFdr27twItZpYP7AQmJzaW\ndGfTXrX2G2g5mRE+e+5MVm6p4tHXd4UdR0RERIaAIIX2S2Y2HLgLWA68GNskBI0trWzbV6/WfiF4\n78nFTCsaxnf/tpaW1raw44iIiMggF6TryNXuXuXutwDvBK529w8lPpp0ZmtlHW0OU0bnhh0l5aRH\n0vjCebNYv7uWP6woDzuOiIiIDHJBHoZ8tH3f3cvcfWX8ORlYG/fUATBldF7ISVLT2+eM5aSS4Xz/\nsTeoa2oJO46IiIgMYl0W2maWGVsBcqyZ5cetClkMqI92SDbuOQDAFLX2C4WZ8eULZlNR08gdT28M\nO46IiIgMYt2NaH+S6CqQszi0IuRq4BGirf4kBBv31DEiN4PC3Iywo6SsU0pHcv6x47jtqfXs3N8Q\ndhwREREZpLostN39++4+Cfiiu5fErQr5Fnf/wQBmlDib9qjjyGDw5fNn09rmfO/RdWFHERERkUEq\nyMOQPzCz+WZ2qVaGDF97D20JV8moXK5cWMq9K8tZtW1/2HFERERkEAryMOTPgR8BbwNOj22LEhtL\nOlPf1MrO6gamqtAeFD559nRG5Gbyn395XYvYiIiIyBHSA1xzKjDH3dU4OGSb9tYCaER7kCjIzuCz\n587khj+t4uFVO7nguPFhRxIREZFBJMiCNauBokQHkZ5t3htt7VeqjiODxvtPmcSscfn8919ep76p\nNew4IiIiMogEKbQLgdfN7K9mdl/7luhgcqQtldER7UkjtVjNYJEeSeM/Fx/L9v0N3PpkWdhxRERE\nZBAJMnXkW339cDM7D7gJiAA/dfdvd3jfYu9fANQBH3H3lbH3PgtcBTjwGnClu6d0L7UtlXUMz82g\nMEet/QaT+VNGctHcCfzkqQ1cclKxpvaIiIgIEKzryOPAK0BGbHsldq5bZhYBbgHOB+YA7zezOR0u\nOx+YEduWAD+O3TsRuB6Y5+7HEi3ULwv4MyWtzXvrKNFo9qD05QtmkxEx/usvr4cdRURERAaJIF1H\nLgFWAlcAHwKWm9m/Bvjs+UCZu29w9ybgHmBxh2sWA3d71AvAcDNrf6IsHcgxs3QgF9ge6CdKYlsq\nVWgPVmMLsvnM22by+NoKHnt9V9hxREREZBAIMkf768Ap7v4Bd78cWAD8e4D7JgJb447LY+d6vMbd\ntwHfA7YAO4D97v5ogO9MWi2tbWzbV69CexD7yMJSjhmbzzceWEVtY0vYcURERCRkQQrtNHePH6Kr\nCHhfn5nZCKKj3VOACcAwM/tgF9cuMbPlZrZ89+7diYwVqh37G2hpcyaPUqE9WGVE0vjmxdEHI298\n7I2w44iIiEjIghTMj8Y6jnwwVuz+GXgkwH3bgElxx8Wxc0GueRuw0d13u3szcB9wWmdf4u63u/s8\nd59XVJS8XQi3VEZb+6njyOB28uSRXL6ghJ/9c6NWjBQREUlxQQrtfwN+QXTO9fzY/v8LcN8yYIaZ\nTTGzTKIPMz7Y4ZoHgQ9Z1KlEp4jsIDpl5FQzy411JjkHWBPoJ0pS7T20J6uH9qD3xXfMYuSwLL58\n32u0tmnFSBERkVQVpOuIA/8gOor9MPCkB1hv2t1bgOti960Bfu/uq83sGjO7JnbZQ8AGoAy4A7g2\ndu9S4I9EH8J8LZbz9t79aMllS2UdGRFjXEF22FGkB4W5GXz93XN4bdt+fvbPjWHHERERkZD02Efb\nzK4E/hN4CjDgNjP7urv/oqd73f0hosV0/Lnb4vYd+GQX934D+EZP35EqtlTWMmlELpE0CzuKBPDu\n48fzwEvb+N9H1nHO7LFMUW9tERGRlBNk6siXgJPc/YPu/gFgHvCVxMaSjrZU1ml+9hBiZnzz4uPI\nSk/jC398hTZNIREREUk5QQrtSqAq7rgqdk4GiLuzeW+dOo4MMWMLsvn6u9/Csk37+MXzm8KOIyIi\nIgMsSKG9DnjezL5mZl8FngPWmtn1ZnZ9YuMJwP76ZmoaWtRDewi65KSJnHVMEd/521o27akNO46I\niIgMoCCF9laiD0FmAdnA34guLFMU2yTB2juOqNAeesyMb118PBmRND73+5dpaW0LO5KIiIgMkB4f\nhnT3G9r3zSzf3WsSG0k6au+hXaKpI0PSuMJs/vuiY/n0PS9z65Pruf6cGWFHEhERkQHQ5Yi2mX3V\nzGbF9jPN7FGg3Mx2mdnZA5ZQDhXaGtEeshbPnchFcydw0+Nv8tKWfWHHERERkQHQ3dSRy4nOzwb4\nENGpI6OBs4FvJTiXxNm8t5ai/CxyM3v8BYQMYv+x+FjGFWTz2d+9TG1jS9hxREREJMG6K7Sb4ham\nOQ/4rbs3u/tqICPx0aTdlso6jWYngcKcDG689AQ2V9bx9QdWhx1HREREEqy7QrvRzGab2Siio9iP\nxr2Xk9hYEm/L3jomq9BOCgumjuJTZ8/g3pXl/GH51rDjiIiISAJ1V2h/HniQ6PLoN7v7BgAzu4Do\nsugyABpbWtlR3aDFapLIp8+ZwalTR3LDA6t4Y5eeLRYREUlWXRba7v5Pd5/h7iPc/d/jzj/k7pcO\nSDqhfF+CrMHCAAAgAElEQVQ97mixmiQSSTNuvuxE8rIyuPbXK6lr0nxtERGRZBSkj7aEqL3jiEa0\nk8uYgmxuumwu63cf4Mv3vcahxyFEREQkWajQHuTK99UDau2XjBZOH83nz53JAy9v56fPbAw7joiI\niPQzFdqDXPm+OjIjaRTlZYUdRRLgk2dN5/xjx/Gth9fw7Jt7wo4jIiIi/ajLxsxmdmF3N7r7g/0f\nRzoq31fPxBE5pKVZ2FEkAcyM7733BNbvPsB1v13Jg59cpBVARUREkkR3I9rvjW2fAH4JfCy23Q1c\nk/hoAtFCu3iEuikms2FZ6dx+xTza2pyP/WIZ++ubw44kIiIi/aC7riNXuPsVQASY4+6L3X0x8Jbu\n7pP+tW1fnQrtFFA6ehi3ffBkNu6p5dpfr6CppS3sSCIiInKUghTMxe6+Le54O1CSoDwSp76plT0H\nmigeoakEqeC06aP59iXH88+yvXztT+pEIiIiMtR1OUc7zpNm9lfgt7Hj9wFPJiyRHLStKtraTyPa\nqeM9JxezZW8tNz9RRvGIXK4/Z0bYkURERKSPghTanyQ6V/v02PHdwB8TlkgO2hpr7adCO7V89tyZ\nlFfVc+NjbzBiWCZXnDo57EgiIiLSBz0W2h79/fXvY5sMoPKDhbamjqQSM+M7lxxPdX0zX39gFQXZ\n6SyeOzHsWCIiItJLPc7RNrPFZrbGzPabWbWZ1ZhZ9UCES3XqoZ26MiJp/Ojyk5hfOpLP//4Vnli7\nK+xIIiIi0ktBHob8P+BSdy909wJ3z3f3gkQHE/XQTnXZGRF++uF5zB5fwDW/XMk/1lWEHUlERER6\nIUihvcvdX0t4EjmCemhLfnYGv/zYfGaOy+Pqu1fwj7UqtkVERIaKIIX2MjP7tZm918wubN8SnkzU\nQ1sAGJ6bya8+toBjxuVz9S9X8PgaTSMREREZCoIU2qOANuBCDq0W+Z5EhhL10JbDtRfbs8ZHi+0/\nvbSt55tEREQkVEG6jlwxEEHkcOqhLR0V5mbw66sWsOTuFXzmdy+zr66JKxdOCTuWiIiIdKHHQtvM\nsoCPEF16Pbv9vLsvSVwsUQ9t6Ux+dgY/u/IUPnPPy/zHn19nd00j//b2Y/TArIiIyCAUZOrI3UAp\n8C5gKTANaEhgJkE9tKVr2RkRbvnASVy+oIRbn1zPp377EvVNrWHHEhERkQ6CFNoz3f3LwAF3vxM4\nD5if2FiiHtrSnUia8T8XHctXL5jNQ6t28L7bn2dXtf7+KyIiMpgEKbSbY69VZjYbyAfGJC6SgHpo\nS8/MjI+fMZXbr5hHWcUBLvzRsyzfVBl2LBEREYkJUmjfaWYjgG8AjwBvAN9LaCpRD20J7Nw5Y7n3\nE6eRnRHhsttf4KfPbMDdw44lIiKS8nostN39J+6+z93/4e4l7j7a3W8diHCpTD20pTdmjy/gz59a\nxNmzxvDff13DJ361kqq6prBjiYiIpLQgI9oywNRDW/qiIDuDn1xxMl+9YDaPr93FO37wNM+8uTvs\nWCIiIikroYW2mZ1nZuvMrMzMvtTJ+2ZmN8fef9XMTop7b7iZ/dHM1prZGjN7ayKzDibqoS191T5v\n+/5rF5KfncEVd77Ivz+4mtrGlrCjiYiIpJwe+2j3lZlFgFuAc4Fyoku5P+jur8dddj4wI7YtAH4c\newW4Cfibu7/HzDKBlBneVQ9tOVrHTizkL59axLcfXsvPn9vEY6/v4r8vOpazZuk5ZpFk4+64g7fv\nQ+w4ep7YcdvB9w5dQ9x17e+1xc5x8NyR17gfOt8bRvAH/K2fewEE/TwLeGFnV3V2a8efufNrgn1B\nXz+r48/U1+ydZhrgDB2viaQZGZHBO0EjyII1S4G7gN+6e3UvPns+UObuG2Kfcw+wGIgvtBcDd3v0\nya0XYqPY44E64AyiC+Xg7k1Aykw4VQ9t6Q/ZGRH+/cK38K7jx/Ol+17jyp8v413Hj+crF8xmwnD9\nJU5ST0trG40tbTQ0t9LQ/trcSmNLG43NbbS0tdHc2kZTix/cb2716GtLGy1tTlNrGy2xc/H7B6+L\nnYu+F72nzZ3WNqetDVrb92Ov8fttTifn4q/jiPNteu5ZUtznzp3J9efMCDtGl4KMaH8YuBJ42cye\nA37m7o8HuG8isDXuuJxDo9XdXTMRaAF2Az8zsxOAFcCn3b02wPcOeeqhLf1pXulI/nr9In7y1AZu\n+UcZf1+zi6vPmMbVZ04lNzNhv9QS6bXm1jbqGls50NRCXWMLtU2t1Da2UNvYQl1TK7VNLbHjVuqa\nou83NLXS0NJKQ/OhwrmhuY3GlsNfG5pbaenHqtQMMiJpZEbSSI/YEfvpsVG2SJpFNzPS0iAjLY00\niz8XfY2kte/TybnovhlHnE+z6EihWXTkL/oazdc+gmgGaWaHzseui74Xfz763Ra7qeP1RxzHXRNE\nb5ohBb00aIelwF8d8MLORvI7i9LxVOfX9O2zOruos/gdL+vszyzQfZ1eE/DP/4jP6r8/v1NKRwTK\nEJYe/y/r7muBL5rZV4ALgbvNrInoKPcP3b0qQblOAj7l7kvN7CbgS8ANHS80syXAEoCSkpIERBl4\n6qEt/S0rPcL158zg4pMm8q2H13LT42/yu2Vb+dQ503nvyZPITB+8v3aToaGtzalpbKG6vpnqhmaq\n61tir81UN3R/vraphbrGVppa2wJ/X05GhGFZEbIz2rc0stMj5GRGGJGbSXZGhKyMtOh76bH99Nh1\nGYdes9LTyIq9RgvlNDIidth+RiQtth3aj+i/zyISQKDhLDObQ3RU+93AA8CvgUXAE0QL4s5sAybF\nHRfHzgW5xoFyd18aO/9HooX2Edz9duB2gHnz5iXFL9HUQ1sSpXhELrdcfhIfOa2Sbz20hq/ev4of\nP7me68+ewUUnTlTBLUB0ikVVfTP7apuorG1iX10TlbXNsdem6Pm6Q69Vdc0caGzpccQyPzudguwM\nCnIyKMhOp3hEDvnj88nLSic3M528rAi5mekMi71Gz0cYFnvNy0onNyudnIyICl0RGRKCzNF+keic\n6buAr7t7feytf5rZwm5uXQbMMLMpRIvny4DLO1zzIHBdbP72AmC/u++Ife9WMzvG3dcB53D43O6k\ntm1fHXPmjA07hiSxU0pHcu8nTuPJN3Zz46Nv8IV7X+XGx97go4tKef/8EvKzM8KOKP2sobmV3TWN\nVNQ0sudAI7trYltsf++BRvbVNVNZ28T++uYuP2dYZoQRwzIZOSyTEbmZTC3KozAng8KcQwV09DWD\ngpxDhXVeVrqKYxFJOUFGtD/o7m/EnzCzEnff4u4XdnWTu7eY2XVEV5OMAHe5+2ozuyb2/m3AQ8AF\nQBnRYv7KuI/4FPDrWMeRDR3eS1rqoS0Dxcw465gx/MvMIp5ct5ufPL2ebz60lh8+Xsa/njSRyxeU\nMGtcQdgxpRvuTnV9Czuq69mxv+FQ8RxXQO+JHdd00uLRDEYNy2R0Xhaj87KYOCKXkbkZhxXS8a/D\nczPIzoiE8JOKiAxN1tNEdjNb6e4n9XRuMJg3b54vX7487BhHpayihrfd+DQ3XTaXxXMnhh1HUsyr\n5VXc9exGHlq1k6aWNk6ePIJLTirmguPGMTw3M+x4KaWtzdlb28TO/Q3s2F/PruoGduxviB03sLM6\nul/f3HrEvflZ6RTlZzE6P4ui/CyK8mKvccdj8rMYOSyT9EHcFktEZLAysxXuPq+n67oc0TazmcBs\noNDM4keuC4Dso48onVEPbQnT8cXD+cFlJ/KN2ibuXVnOb1/cwlfuf41vPLiKM2eO4fxjx3H2rDGM\nGKai+2i0tLZRUdN4sHCOFs31hxXSFTUNNLcePhCSnmaMLchmXGE2cyYUcM6sMYwrzGZ8YQ7jCrMY\nk5/N6LwscjI16iwiMhh0N3XkLcDFwHDgvXHna4CrExkqlamHtgwGI4ZlctXpU/nYoims3l7Nn17a\nxp9f3c7f1+wikmacUjqCM2eOYdH00bxlQoE65MRpaG49WDwfPgpdf/D87prGI/ofZ6WnMb4wWkTP\nnzIyVkBnM7Yg++D50cOy9GctIjKEdFlou/v9wP1mtsjdnx3ATClNPbRlMDEzjp1YyLETC/nKBbN5\nbdt+Hnt9F39fs4vv/G0t3wGG52ZwUskIji8u5IRJwzmheDgjk3DE292pqmtmx/5oAd0+dSO+qN5Z\n3UBV3ZEPEuZnpTMuViwfMy6fcYU5jIsroMcXZlOYkxF4RToRERkaups68nl3/z/gEjO7uOP77v65\nhCZLUeqhLYNVWppFC+lJw/m3dxxDRU0Dz5Xt5Z9le3h5axX/WFdxsL1b8YgcTigezjHj8pk8KpfJ\no4ZROip3UM7zbmxpZe+BJvYciHbj2FPTxO72/QNN0QI6Vlw3thze59kMRudlMa4gm+IRucwrHcH4\nwpzo9I6CbMYVZjGuMIe8LC0MJCKSirr7r//62OuqgQgiUeqhLUPFmPxsLjpxIhedGH1ot6ahmVXb\nqnm1vIpXy/fzSnkVf31tx2H3FOZkUDoql/GFOYzMy2RUrLvFyGGZjBqWRUFO+sEFRrIzoguJZGdE\nFxLx2PLTTnR1MCe6pHV0Oe3o6n/17asENrVS3dDM/vroVl3fcnB/f2zhlP31zeypaaS64chuHAB5\nWemMzstkTEE2cycNZ1zcNI72edJj8rPI0MOEIiLShe6mjvwp9npn+zmL/l4zN1WWQg+DemjLUJWf\nncFbp43irdNGHTzX0NzK1so6Nu2tY/PeWjbtrWXz3jo27DnA8s3RxU/6cVXsLqUZFMR6PbdvEwpz\nGDX9UGu70XmZ0U4dsWM9UCgiIkcryII1dwPXAS3Ai8AoM/tfd78x0eFSjXpoS7LJzogwY2w+M8bm\nd/p+W5uzv76ZvbWN7D3QRE1Dy8HR6YbmVhqaW2lsaaOppY00M8zAiE7ZMDPSzA5fUjv90JLcBTnp\nB4vqvKx0zX8WEZEBF2Ti4PHuXm1mlwOPAV8ElgMqtPvZtqo6QK39JHWkpRkjhmUyYlgm08eEnUZE\nRKR/BZlcmGFm6cBi4AF3bwLaerhH+kA9tEVERESSR5BC+6fAFmAE8JSZlQAHEpoqRamHtoiIiEjy\n6LHQdvfvu/sEd3+7R9dr3wqcnfhoqUc9tEVERESSR5CHITOBi4DSDtd/M0GZUpZ6aIuIiIgkjyAP\nQ94PNAArgNbExklt6qEtIiIikjyCFNqT3f3YhCcR9dAWERERSSJBHoZ8wczmJDxJilMPbREREZHk\nEmREewHwkpmVAY1E14twdz8poclSjHpoi4iIiCSXIIX2RQlPIeqhLSIiIpJkgrT3Ww8UAQtj+1VA\nc6KDpZr2HtoTh2vqiIiIiEgyCNLe72vAQmAacDeQDfwGWJTYaKmlvYf2mHz10BYRERFJBkEehnwP\ncAFQC+Du24CCRIZKReqhLSIiIpJcghTajbEVIR3AzDS3IQHUQ1tEREQkuQQptO8zs1uAQjO7EngU\n+FliY6WebfvqmDhchbaIiIhIsuhxjra7f8fMzgeagBOA/3H3hxOeLIU0NLf30FahLSIiIpIsgjwM\n+U13/wrwcCfnpB+UH2ztp1k5IiIiIskiyNSR8zo5987+DpLKyvdpsRoRERGRZNPliLaZXQ1cA8w0\ns5Vxb+UDKzu/S/pCI9oiIiIiyae7qSO/Bx4HvgV8Ke58jbtXJDRViinfV09GxNRDW0RERCSJdFlo\nu/s+YB/wXjMzoqtDpgPpZjbB3bcPUMakVx7rOKIe2iIiIiLJI8jDkJ8A/gvYC7TFTjswJ4G5Ukq0\nh7amjYiIiIgkkx4LbeDzwGx3353oMKmqfF8958waE3YMEREREelHQbqOlAOViQ6SqqI9tBvVcURE\nREQkyQQZ0S4DnjCzvwCN7Sfd/eaEpUoh26piHUdGqtAWERERSSZBCu0dsa0gwVlSklr7iYiIiCSn\nIEuw3zAQQVKVFqsRERERSU7dLVhzP9HuIp1y94sTkijFHOqhnR12FBERERHpR92NaP/oaD/czM4D\nbgIiwE/d/dsd3rfY+xcAdcBH3H1l3PsRYDmwzd3fdbR5BqPyffVMGJ5DRD20RURERJJKdwvWPH40\nHxwrkm8BziXauWSZmT3o7q/HXXY+MCO2LQB+HHtt92lgDUk8P7x9sRoRERERSS5B2vv11XygzN03\nuHsTcA+wuMM1i4G7PeoFYLiZjQcws2LgncBPE5gxdNv21Wt+toiIiEgSSmShPRHYGndcHjsX9Jof\nAF/g0GqUnTKzJWa23MyW7949tNbUaWhupaKmUR1HRERERJJQIgvtPjOzdwEV7r6ip2vd/XZ3n+fu\n84qKigYgXf/Z3t5DWyPaIiIiIkknkV1HtgGT4o6LY+eCXHMJcKGZXQBkAwVm9it3/2AP3zmkqIe2\niIiISPLqbkT7R0QfZiwnOn3jl7GtBdgS4LOXATPMbIqZZQKXAQ92uOZB4EMWdSqw3913uPuX3b3Y\n3Utj9z2RbEU2xBfaGtEWERERSTY9dh0xs++4+7z282b2J+DFnj7Y3VvM7DrgEaLt/e5y99Vmdk3s\n/duAh4i29isj2t7vyqP4WYac8n11pKcZYwvUQ1tEREQk2QRZgj3PzErdfVPsuATIC/Lh7v4Q0WI6\n/txtcfsOfLKHz3gSeDLI9w015fvqGT88Wz20RURERJJQkEL788AzZrYOMGA68ImEpkoR5fvqKB6u\n+dkiIiIiyajHQtvd/2pmM4E5sVOvu3t9YmOlhm1V9ZwxY2h1ShERERGRYHps72dmOURXaPx4rN3e\nRDM7P+HJklxjSyu7qtVDW0RERCRZBemjfVfsukWx4+3ANxOWKEVsr2oA1HFEREREJFkFKbRnuPs3\ngWYAd68jOldbjkL5vjpAhbaIiIhIsgpSaDeZWTaxxWvMbArQlNBUKeBgD+2RmjoiIiIikoyCdB35\nL+BvQLGZ/QI4E7gqoalSQPm+OiJpxtj8rLCjiIiIiEgCBOk68rCZLQdOIzpl5P+5e0XCkyW58n31\njC/MJj0S5JcKIiIiIjLUBOk68qi773b3B9z9T+5eYWaPDkS4ZLZtX73mZ4uIiIgksS4LbTPLNLMC\nYKyZ5ZtZQWwrJro6pByF8n31au0nIiIiksS6mzrySeBzwBhgNYc6jVQDt3V1k/SssaWVXTUNGtEW\nERERSWJdFtru/n3g+2b2BXf/bvx7ZhbkIUrpwo6qBtzRiLaIiIhIEgvyJN5lnZx7sb+DpJKDrf00\noi0iIiKStLocmTazMcB4IMfMjuPQ1JECQEOxR6F9sZqJw1Voi4iIiCSr7qaAvBP4KFAM3Bp3vga4\nIZGhkl35vnoiacb4wuywo4iIiIhIgnQ3R/tnwM/M7FJ3//0AZkp626rqGVegHtoiIiIiyay7qSPv\nd/ffAuPN7PqO77v7zQlNlsTK99VpfraIiIhIkutu6siI2OvogQiSSsr31XPaNP2xioiIiCSz7qaO\n3Bp71XzsftTU0sbOavXQFhEREUl2PfbDNrMS4DqgNP56d784cbGS14799bEe2iq0RURERJJZkIVn\nHgTuBh4D2hIbJ/m199CeqEJbRHqjrQ28FdpaoC326m0djlsP7R9xLu4eHNxjr21x+3TyvvdwT4dr\nO7L2zrDWybmuru3ueuv9tT2e7/iZXR13laG7495cG/C4Tzn6Oye9uD4RfwZ9ydGf+vnzlK/v0tIh\nMnjXUQySrMndb0x4khTR3kN7klaFFBl62lqhuQ6a66GpNvraXA/N7ft10FR36Jrm+P16aG2Kbi2x\n19ZGaG2Glsa492LnOr5HJ0WsiEiqO+urcOYXwk7RpSCF9g/N7GvAI0Bj+0l3fzVhqZLY1sp60gzG\nqYe2yMBrboD6SqjbG90a9kNDNTRWd3iNO99Yc2i/paH33xnJgoyc6BbJhPSs6Gv7lp4JmcOi10Uy\nYu9nxI5j70cyo6M2aRGwyKH9tHSwtA7HkUP7R5yLRK9vHxGM3yd2fHCfTs71dE8Xo1SHjXR7P5zv\nj8+Ine/4md7hniPe78Xx0dzbZZZ+yBU4J0de22+f3dvj+Cx9/Kz+0tlvbo7uA/v541Is3+TT+vfz\n+lmQQnsmcBVwPoemjjhwRqJCJbPNlXVMHJFDhnpoixy95nqo2QkHdkW39gK6Lq6Yjj/XdKD7z8vM\ng6wCyC6IvuaOghFTYsf5kJl/qGjOHBbbz43bciAzbj8jN1rciohISgpSaL8fKHX3xh6vlB5t2VtL\n6ahhYccQGdzcoXYP7N8a27bBgZ3Rorq9sK7ZER2R7kxmPuSOjBbKuaNh9DGx/fZzsf3s4dECur2w\nVlEsIiL9KEihvRrIJ27aiPTd5so63nnc+LBjiITLHQ5UQOUG2LcJ9pfD/i1QtTW2Xw4t9YffE8mE\nvHGQPxZGz4DS0yF/XHTLGwd5Y2BYUbSATs8K5ccSERGJF6TQzgfWmtlSDp+jrfZ+vbS/rpmqumYm\nj9KDkJIC2kelK9fD3vXR18oNsf2N0FRz+PXDiqBwEoydAzPfEd0fPgkKi6GgOFpA9/uT7yIiIokT\npND+n4SnSBGbK2sBmKypI5JM3KF2N1S8DhVrYfcaqFgDu9cePrXDIjC8BEZNg5K3Rl9HTo3OgS6c\nGJ3TLCIikkSCFNrPAQ3u7mY2DTgGeDSxsZLTpr3R1n4a0ZYhq7kedq2GHS9Hi+mKtdECu77y0DU5\nI6BoNhx7CYyeCSOnRYvq4SXRbhoiIiIpIkih/QxwhpkVAk8AK4HLgA8lMlgy2rI3OqJdMlKFtgwB\nTXWwaxVsfzlaWO94JVpce2v0/awCKJoFs98FY+ZE98fMhryxmuIhIiJCsEI7zd3rzOyjwI/d/dtm\n9nKigyWjTXvrGJOfRW7m4F3BSFJUW2u0iN66FMqXRYvrPetiKwIS7dwxYS7MPC/6Ou746Ai1CmoR\nEZEuBSq0zewU4APAx2Pn1AOrD7bsrVNrPxkcGmugfHm0sN66NLrfWB19b9gYmHAizLkQxs+F8SdA\nwQQV1SIiIr0UpND+HPAfwF/cfZWZTSU6nUR6adPeWs6YWRR2DElFtXtg49Ow+Z+wZSlUrI6NVhuM\nfQsc916YtABKFsDwySqqRURE+kGPhba7P0F0bnb78Qbg2kSGSkb1Ta1U1DQyWfOzZSA01sDm52DD\nU9ECe9dr0fOZeVA8D874AkyaD8WnRBdrERERkX7XY6FtZtOJjmqXxl/v7m8PcO95wE1Ep5r81N2/\n3eF9i71/AVAHfMTdV5rZJOBuYCzR5d5vd/ebAv5Mg9KWyljHkdGaOiIJ0NIYnQKy8elocb1tRfSh\nxUhWdJT67BtgypnRKSERPSMgIiIyEIL8H/ePwJ3Ar4DWoB9sZhHgFuBcoBxYZmYPuvvrcZedD8yI\nbQuAH8deW4DPx4rufGCFmT3W4d4hZVOs44hGtKXfVG+HNx+FNx6FDU9Ccy1YGkw4CRZ9JlpYT5qv\n/tQiIiIhCVJot7n7D/vw2fOBsthUE8zsHmAxEF8sLwbudncHXjCz4WY23t13ADsA3L3GzNYAEzvc\nO6RsifXQ1sOQ0mdtrbBtJbz5CLzxN9gZmw5SUAwnvA+mnwulCyG7MNycIiIiAgQrtB8wsyXA/Ry+\nBHt1D/dNBLbGHZcTHa3u6ZqJxIpsADMrBU4ElgbIOmht2ltLYU4GhblasEN6ofEAlD0G6/4Wfa3b\nGx21nrQAzvlGdKnyMXP08KKIiMggFKTQvir2ekPcOQdK+j/O4cwsD7gX+ExXhX3sLwFLAEpKEh6p\nz7ZU1lGqFSEliIb98MYj8PoDUPZ3aGmAnJEw/W3Rwnra2ZA7MuyUIiIi0oMgXUcm9fGztwHx9xbH\nzgW6xswyiBbZv3b3+7rJdztwO8C8efO8j1kTbtPeWuZOGhF2DBms6iph3UPR4nrDk9DaBPnj4aQP\nwZzFUPJWSFP7ehERkaEkUPsBM5sFzAGy28+5+296uG0ZMMPMphAtni8DLu9wzYPAdbH52wuA/e6+\nI9aN5E5gjbvfGOgnGcSaWtrYtq+ei+ZODDuKDCZ1ldHC+vU/wcZnol1ChpfA/CXR4nriPEhLCzul\niIiI9FGQ9n5fA94OzAIeAd4BPAt0W2i7e4uZXRe7JwLc5e6rzeya2Pu3AQ8Rbe1XRrS935Wx2xcC\nVwCvxS33/hV3f6h3P97gsK2qnjaHEnUckeZ6WPcwvPYHePMxaGuGUdNh4acPrcSo+dYiIiJJIciI\n9vuAucBKd7/CzMYDPw/y4bHC+KEO526L23fgk53c9yyQNNXG5vbWfuo4kpraWqP9rV/7A7z+IDTV\nRKeFLLgajr8Uxh2v4lpERCQJBSm069291cxaYj2tdwKTE5wrqbQvVqOHIVOI+/9v786D/K7rO44/\n37ubg9wJuSAHG0I4QrjiJkC1VozKZY3WqeCBlNGqM1pt7UyrPR2nnTodq3VGx4ocYrUyFuzAIEUZ\ni9Z2pmQ3CUdIICeYDSGbGEiyIdfuvvvH94dEJGEJv+9+d3/7fMzs/G72lXwhefHdz/f9ge0PF+X6\n0Tug+xkYNQHOXQ7nvQda3+Caa0mSGlx/ivbqiJgE3AJ0AHuBFaWmajBbdu3npBHNTBs/quooKtvz\nu+GR78Pqf4Uda6BpRDEp5LzfL27dPEaSpGHjuEW7dlHi5zLzOeBrEfEjYEJmrhqQdA1iY1c386eP\nJVwe0Jj6emHzA7D6O/D4D4uJIadeBFf/E5z7e47ikyRpmDpu0c7MjIj7gUW1xxsHJFWD2djVzaWn\nn1x1DNXb7i3w0L8VX3s7i1nXSz4MF74fZi6qOp0kSapYf5aOPBQRF2Xm6tLTNKC9B4+wfc9Bzpgx\nruooqofeI8VZ646biwscCThjGVz+93DWldDi8iBJklQ4ZtGOiJbM7KHY/rw9IjYB+ymmgWRmLh6g\njEPaxq5uABZMH19xEr0m+3bAqtug41bY93Qx7/qyv4IL3wsTZ1edTpIkDULHO6O9AlgMvGOAsjSk\njdoenJwAAAyiSURBVDuKon2mZ7SHnkzo7IAVN8Jj/1HMvJ6/DN7+ZVjwVqeGSJKk4zpe0Q6AzNw0\nQFka0vod+xjV0sTsyY72GzKOHIA1PygK9vaHirF8Sz5cfE09o+p0kiRpiDhe0Z4WEZ8+1ouNsDX6\nQNjQ1c38aeNobnLiyKD37FPQcQus+jYc2A3TzoGrvwTnXwOj/ImEJEl6dY5XtJuBcTTQDo1V2NjV\nzZLWyVXH0LFkwuafwopvwvr/BALOvhqWfqTYVMaRjJIk6QQdr2hvz8zPD1iSBtR9qIdtzx3gfTPm\nVh1FL3VwLzx8O7R/E3athzFT4Q2fhrYbvLhRkiTVxSuu0daJe3HiiMsOBo2dTxRnrx/+Hhzuhlmv\ng3d9A859l6P5JElSXR2vaC8bsBQNasOOfQAsmOFov0r19sD6+4qLG7f8DJpHwqJ3w5I/hNmvqzqd\nJElqUMcs2pm5eyCDNKKNXd2MbGli7hQnjlRi/y9rs69vgT1bYcJsWPY3sPh6GDu16nSSJKnB9Wdn\nSJ2g9Tv2OXGkCttWFctD1twJvYdg3hvhin+AM6+EZv+VlyRJA8PWUaINXd0snuvEkQHRcwjW3gUP\nfgO2dcCIsbD4umL29fRzqk4nSZKGIYt2SZ4/3EPnswe4pm1O1VEa255tsPJWWPkt2L8TTj4DrvxH\nuOBaGD2x6nSSJGkYs2iX5FcTR7wQsv4y4an/LS5uXHcPZB+cdWVx9vr0y6CpqeqEkiRJFu2ybNjx\nQtF2tF/dHOqGR79frL/uWgujJ8GlH4clH4LJrVWnkyRJ+jUW7ZJs6OpmZHMTpzlx5LX75SZovwlW\nfxcO7YGZ58E7vlqM6Bvp768kSRqcLNoleaTzOc6cOY6WZpcxnJDenmJL9PabYfMD0NQCC99ZbI0+\nZ6lbo0uSpEHPol2CI719rP7Fc1yzxAshX7W9T8Oqb8PK22Df0zBhFlz2l7D4gzB+ZtXpJEmS+s2i\nXYLHnt7LgSO9LGmdUnWUoaGvr9ixseNmePxeyF6Yvwyu/iIsuNzZ15IkaUiywZSgfUuxqeaSVmdo\nH9f+XfDw7cXOjbs3wUlTiosb226AKadXnU6SJOk1sWiXoP3J3Zx28himTxhddZTBp7cHNt4Pq78D\n6++Dvh6YczH8zp/DwuUwwt8zSZLUGCzadZaZdDz1LG8+e3rVUQaXnU8U5frh22F/F4yZChd/DC58\nP8xYWHU6SZKkurNo19mmnfvZvf+wy0YADu6BNXcWY/m2dUA0w5mXw0UfgAVvg+YRVSeUJEkqjUW7\nztqffGF99jC9EPLIAdhwf1Gw198HPQdh2jnwtr+D86+BcZ7plyRJw4NFu87at+xm6riRzJs6tuoo\nA6f3CGz+KTx6Bzz+Qzi8r1gactEH4ML3wamLnXstSZKGHYt2nbU/tZu206YQjV4sDz9fbCTz+L3w\nxA/hwLMweiKcu7zYsbH1jY7lkyRJw5pNqI6e2XOQrbsPcP2lrVVHKUf3zmI5yBP3wqYHoOcAjJoI\nZ76tKNfz3wwto6pOKUmSNChYtOtoRW199tJ5DbI+u+cwdLYXZ643PQDbVgIJE2bD4uvgrKug9Q1e\n1ChJkvQyLNp19ODmXzJmZDMLT5lQdZQT09cLXevgyZ8XxfrJ/4Ej+yGaYNbrilnXZ18FM893zbUk\nSdIrsGjXyZO79vPvKzu5atFMWpqbqo7TPwf3QGcHbF0BWx8s7h/eV7w25XS44FqYfxm0/jacNKna\nrJIkSUNMqUU7Iq4AvgI0Azdl5hde8nrUXr8KeB74g8xc1Z/PDiaZyV/ftYaRzU189qpzqo7zmzJh\n7zbY8RjsWFPcPrMGdq0HsjhjPf1cOP89MGcpzL0UJp9WdWpJkqQhrbSiHRHNwNeAtwKdQHtE3J2Z\na49625XAgtrXxcDXgYv7+dlB455HtvPzDbv43O8uZEZV265nwvO74dktsHvLr992rYODz7343klz\nYcai4gLGuRcXy0JGja8mtyRJUoMq84z2UmBjZm4GiIjbgeXA0WV5OfDtzEzg/yJiUkScArT247OD\nwp4DR/j8PWs5b9ZErqvXtJG+vmKjl0P74NBeOLgXDu0pHh/cA91dta8dR93ugMPdv/7PGX8qTG6F\nc98FMxcV5Xr6OcUYPkmSJJWqzKI9C9h61ONOirPWr/SeWf387KCw7uaPcuOhhzmraSzNNzUVZ5ZJ\nyD5Iavdrj1+4/6vXE7K3mO7Re6i47TkIfUde+RuPmljssjhuBpxyQXE7aS5MmQeT5xVLP0acVOqv\nXZIkScc25C+GjIiPAB8BmDt37oB+78yEEWOYMOlkxkwcD0QxjSOaXrz/q+fi5V+PZmgZCS2joXlU\nMYf6ha9R44tCPXpC7f6E4v7YaZZoSZKkQa7Mor0NmHPU49m15/rznhH9+CwAmXkjcCNAW1tbvrbI\nr05EcMlHv1oUbsfdSZIk6ShlzqFrBxZExLyIGAlcC9z9kvfcDXwwCpcAezJzez8/O2g0/HbrkiRJ\netVKO6OdmT0R8QngRxQj+m7JzMci4mO11/8FuJditN9GivF+Nxzvs2VllSRJkuotioEfjaGtrS07\nOjqqjiFJkqQGFhErM7Ptld43RLYwlCRJkoYWi7YkSZJUAou2JEmSVAKLtiRJklQCi7YkSZJUAou2\nJEmSVAKLtiRJklQCi7YkSZJUAou2JEmSVAKLtiRJklSChtqCPSJ2Ak9V8K2nArsq+L4aWB7n4cHj\n3Pg8xsODx3l4qOo4n5aZ017pTQ1VtKsSER392e9eQ5vHeXjwODc+j/Hw4HEeHgb7cXbpiCRJklQC\ni7YkSZJUAot2fdxYdQANCI/z8OBxbnwe4+HB4zw8DOrj7BptSZIkqQSe0ZYkSZJKYNF+DSLiioh4\nIiI2RsRnqs6j+ouIORHxQESsjYjHIuJTVWdSeSKiOSJWR8Q9VWdROSJiUkTcERGPR8S6iLi06kyq\nr4j4k9qf12si4nsRMbrqTHrtIuKWiOiKiDVHPTclIu6PiA2128lVZnw5Fu0TFBHNwNeAK4GFwHsj\nYmG1qVSCHuBPM3MhcAnwcY9zQ/sUsK7qECrVV4D7MvNs4AI83g0lImYBnwTaMnMR0AxcW20q1cm3\ngCte8txngJ9k5gLgJ7XHg4pF+8QtBTZm5ubMPAzcDiyvOJPqLDO3Z+aq2v19FH8pz6o2lcoQEbOB\nq4Gbqs6ickTEROCNwM0AmXk4M5+rNpVK0AKcFBEtwBjg6YrzqA4y87+B3S95ejlwW+3+bcA7BzRU\nP1i0T9wsYOtRjzuxgDW0iGgFLgIerDaJSvLPwJ8BfVUHUWnmATuBW2tLhG6KiLFVh1L9ZOY24IvA\nL4DtwJ7M/HG1qVSiGZm5vXb/GWBGlWFejkVb6oeIGAfcCfxxZu6tOo/qKyLeDnRl5sqqs6hULcBi\n4OuZeRGwn0H4o2aduNoa3eUU/1N1KjA2Ij5QbSoNhCzG6A26UXoW7RO3DZhz1OPZtefUYCJiBEXJ\n/m5m/qDqPCrF64F3RMSTFMvA3hwR36k2kkrQCXRm5gs/lbqDonircbwF2JKZOzPzCPAD4LcqzqTy\n7IiIUwBqt10V5/kNFu0T1w4siIh5ETGS4mKLuyvOpDqLiKBYz7kuM79UdR6VIzM/m5mzM7OV4r/l\n/8pMz4I1mMx8BtgaEWfVnloGrK0wkurvF8AlETGm9uf3MrzgtZHdDVxfu389cFeFWV5WS9UBhqrM\n7ImITwA/oriq+ZbMfKziWKq/1wPXAY9GxEO15/4iM++tMJOkE/dHwHdrJ0g2AzdUnEd1lJkPRsQd\nwCqKqVGrGeQ7B6p/IuJ7wJuAqRHRCfwt8AXg+xHxIeAp4D3VJXx57gwpSZIklcClI5IkSVIJLNqS\nJElSCSzakiRJUgks2pIkSVIJLNqSJElSCRzvJ0kNJCJOBn5SezgT6KXYdhzg+cx08w5JGiCO95Ok\nBhURnwO6M/OLVWeRpOHIpSOSNExERHft9k0R8bOIuCsiNkfEFyLi/RGxIiIejYj5tfdNi4g7I6K9\n9vX6an8FkjS0WLQlaXi6APgYcA7F7qdnZuZS4CaK3RMBvgJ8OTOXAO+uvSZJ6ifXaEvS8NSemdsB\nImIT8OPa848Cl9XuvwVYGBEvfGZCRIzLzO4BTSpJQ5RFW5KGp0NH3e876nEfL/7d0ARckpkHBzKY\nJDUKl45Iko7lx7y4jISIuLDCLJI05Fi0JUnH8kmgLSIeiYi1FGu6JUn95Hg/SZIkqQSe0ZYkSZJK\nYNGWJEmSSmDRliRJkkpg0ZYkSZJKYNGWJEmSSmDRliRJkkpg0ZYkSZJKYNGWJEmSSvD/CLM5zevz\nSRkAAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fig, ax = plt.subplots(figsize=(12,6))\n", "ax.plot(tlist, np.real(count1))\n", "ax.plot(tlist, np.real(count2))\n", "ax.set_xlabel('Time')\n", "ax.set_ylabel('Transmitted Intensity and Spontaneous Emission');" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### eseries" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [], "source": [ "def solve(E,kappa,gamma,g,wc,w0,wl,N,tlist):\n", "\n", " # Define cavity field and atomic operators\n", " a = tensor(destroy(N),qeye(2))\n", " sm = tensor(qeye(N),sigmam())\n", "\n", " # Hamiltonian\n", " H = (w0-wl)*sm.dag()*sm + (wc-wl)*a.dag()*a + 1j*g*(a.dag()*sm - sm.dag()*a) \\\n", " + 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", " #intial state\n", " psi0 = tensor(basis(N,0),basis(2,1))\n", " rho0 = ket2dm(psi0)\n", "\n", " # Calculate the Liouvillian\n", " L = liouvillian(H, [C1, C2])\n", "\n", " # Calculate solution as an exponential series\n", " rhoES = ode2es(L,rho0);\n", " \n", " # Calculate expectation values\n", " count1 = esval(expect(C1dC1,rhoES),tlist);\n", " count2 = esval(expect(C2dC2,rhoES),tlist);\n", " infield = esval(expect(a,rhoES),tlist);\n", "\n", " # alternative\n", " expt_list = essolve(H, psi0, tlist, [C1, C2], [C1dC1, C2dC2, a]).expect\n", "\n", " return count1, count2, infield, expt_list[0], expt_list[1], expt_list[2]" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [], "source": [ "count1, count2, infield, count1_2, count2_2, infield_2 = solve(E, kappa, gamma, g, wc, w0, wl, N, tlist)" ] }, { "cell_type": "code", "execution_count": 10, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAtoAAAF3CAYAAACbhOyeAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3Xuc1HXZ//HX9Z3ZWVAOKiIiDCwKCnhAZF1NRcvu7p9a\naXflsRStbvNOs+7sYKVZd9120tS6K7O01A6WlQplZnliFQWWgwdEk5PsIgrIgsAuOzvzvX5/zCwM\nK7v7BXZ2dmfez8djHvv5nmau7+4KFx+v7/Uxd0dERERERLpXUOwARERERERKkRJtEREREZECUKIt\nIiIiIlIASrRFRERERApAibaIiIiISAEo0RYRERERKQAl2iIiIiIiBaBEW0RERESkAJRoi4iIiIgU\ngBJtEREREZECiBc7gO60//77e1VVVbHDEBEREZESNm/evHXuPrSr80oq0a6qqqKurq7YYYiIiIhI\nCTOzV6Ocp9IREREREZECUKItIiIiIlIABU20zew0M3vZzJaY2dU7OT7ezJ42sxYz+/xOjsfMbIGZ\n/aWQcYqIiIiIdLeCJdpmFgN+DJwOTATON7OJ7U5bD1wJ3NDB23wGWFyoGEVERERECqWQM9o1wBJ3\nX+buKeAe4Kz8E9x9jbvPBVrbX2xmI4H3Ar8oYIwiIiIiIgVRyER7BFCft92Q2xfVzcAXgbA7gxIR\nERER6Qm98mFIM3sfsMbd50U491IzqzOzurVr1/ZAdCIiIiIiXStkor0KSOZtj8zti+JE4EwzW0G2\n5ORUM/v1zk5099vcvdrdq4cO7bJvuIiIiIhIjyhkoj0XGGdmY8wsAZwHTI9yobt/2d1HuntV7rpH\n3f2jhQtVRERERKR7FWxlSHdPm9kVwN+BGHCHuy8ys8tyx281swOBOmAQEJrZZ4GJ7v5WoeISERER\nEekJ5u7FjqHbVFdXu5Zg7wH1c2BFLVRNhWRNsaMRERER6VFmNs/dq7s6r2Az2lKi6ufgd74f0inC\nWILV77iOfdjEgMPepaRbREREJI8SbYnEV86mfsHDvLrsZd7RmiJuIZ5uYVjtNRhO6qkbee346xjd\nfys2RjPdIiIiIkq0pUsty56Gu8/ioLCVYRZAECPECMwwDzFC0p5ixKxrCc2xeIJg2gwl2yIiIlLW\nlGhLpzY0pfjLn+/hvLCVuIXEzLApF8LgJPQfAg9dDZkUMTM8DAkISadTvP7EHRw0WnXcIiIiUr6U\naEuH3lxcywP3/56lWyqxygSErVgsAZMu2J48D5sIK2qx/kOwh67GMykybgx55V7CJSEWr8SmTVey\nLSIiImVHibbsVOuKZ9j79x/kIk9jlQlip38Xmt98+wx1smaHpNtW1BKsX0mw4C4CQsJ0CltRq0Rb\nREREyk6vXIJdiu/pRx8g7ulsuUjYmk2yp17VecKcrIGpV1FxzEcI4pVkCGjxGI9sPbTnAhcRERHp\nJTSjLW+zYGUjtyw5gOMrK4A0xBLZmeyokjXYtOn4spl8b9EQnn98KaPf+gZja07XzLaIiIiUDSXa\nsoOWZU9Td89vGDJgIq3n3k/itad374HGZA3xZA1fHPk0wd1nEnsuTfjiTwguVjcSERERKQ9KtGW7\n+jkEd5/FJWErH4sniCVmZMtF9kD/12bhlsm2AMykSC99goQSbRERESkDqtGWbTa9/BjW1sYvbM0u\ns76nqqZisQShxWj1OL9YOWLP31NERESkD1CiLdvcu66KVuK4xXa9LrsjyRqYNp3g1K/ypyN+wj8X\nv8Hie78O9XP2/L1FREREejGVjggA6za38L1Fg9h8yM1cecjr3bvQTK4F4HmvPsPZiz5FbFGa8OWf\navVIERERKWlKtAWA259cTks65L3vPQuGDijIZ8RXPkXMMpiHZNIpfHktpkRbRERESpRKR4RNrzxF\nYtZNXDF2PYcUKMkGttdrEyPlcWb7hMJ9loiIiEiRaUa73NXPod/v/oNPWyvBa/dD/WGFK+fI1Wv7\nspl8tW4QTz2V4B/HtTKoX0VhPk9ERESkiDSjXeYyy2Zu6zQSdFenkc4ka4id8nmmnXMOaza18Ps/\n/wlqb9TDkSIiIlJyNKNd5ubbERzhcYIgQ9BdnUYimJTch68euYmPvPxp/JUMFkvAtOl6OFJERERK\nhhLtMnfHyqE0x67j9lO2wsEn92iie+HwemL/SmMe4pkUtqJWibaIiIiUDJWOlLENTSkeWbyGMZPf\nReyUz/d4kls59hQIEqQ9IAwqemw2XURERKQnKNEuYzOeW00qE/KhY0YWJ4BkDUx7gLv6fZRPV3yd\nluFTihOHiIiISAEo0S5jf5rXwPgDB3L4QYOKFkN89PGM/dDXeHDDKP7+0HQ9GCkiIiIlQzXaZWrp\n2s0srN/AV8+YgJkVNZaTDx3Kx6vW8p66L+CmByNFRESkNGhGu0w9/fjfuDz+AB8atqrYoQBw2ejV\nVHga8wxkUoVvMygiIiJSYEq0y5CvnM2HF32Kz8XvZb97z+4VpRpDj3g3YVBB2gM8pgcjRUREpO9T\nol2G1r7wCHFPEyPsPbPHyRrWf+iP3Jw5m9vG3KyyEREREenzlGiXodrW8bQSxy0GPbhITVcOPOJk\nGqd8mu8vGkz9+qZihyMiIiKyR5Rol6HfvXYg1w6+Hjv1q73uocNPnzqOWGDcP+M+dSARERGRPi1S\n1xHLtqUYmn++u79WqKCkcNZvSTF/ZSMnnnoKTD202OG8zYGD+3H1EW9x3uLP4ivUgURERET6ri5n\ntM3sU8BaoBZ4JPf6Z4HjkgJ54l9rCB3ePeGAYofSoXOGvkoF6kAiIiIifVuU0pHPARPc/TB3n5B7\nTSx0YFIYjyxew9CBlRxx0OBih9KhvQ99Jx7LdiAJ1YFERERE+qgoiXYDsL7QgUjhtWZCnvjXWk49\n7ACCoLiL1HQqWcOmc//ELeE5/Dj5A5WNiIiISJ8UpUZ7CfComf0FaGnb6e4/LFhUUhB1KxrZtDXN\nqb24bKTNfodNZdOx+3H3M6/ygfVNJPfbq9ghiYiIiOySKDPaq4GZwCCyD0S2vaSPebnun3y6Yjon\n91te7FAiueyUQ4iZ8Ze/PqAOJCIiItLndDmj7e7XAphZ/9x2c9Q3N7PTgFuAGPALd/9Ou+PjgV8C\nxwBfdfcbcvuTwF3AMMCB29z9lqifKztRP4fzFl9BRSxN7HcP9IlOHgcO7sdVEzdw0b+uxJepA4mI\niIj0LVG6jkw0s7nAK8ArZjbbzCZEuC4G/Bg4HZgInG9m7R+iXA9cCdzQbn8auCr30OXxwOU7uVZ2\nweaXH+t9q0FGcJ46kIiIiEgfFaV05DbgK+4+0t1HAl8Ffh7huhpgibsvc/cUcA9wVv4J7r7G3ecC\nre32r3b3+bnxJmAxMCLCZ0oHFgRH9MrVILsyeMKphEG2A4mrA4mIiIj0IVES7YHu/o+2DXf/JzAw\nwnUjgPq87QZ2I1k2sypgMjB7V6+V7R7aOIpP8DXCd36lb5VfJGtYfdYf+EH6bP4w8Sd9J24REREp\ne1G6jqwwsy8Dd+e2PwqsKFhEecxsAPAn4LPu/lYH51wKXAowatSongirT5qzfD0jqo4jdkrfS1RH\nH/1OXlq4N/cs2sCZ78vQPxErdkgiIiIiXYoyo/0xIAk8mHslc/u6sip3bpuRuX2RmFkF2ST7N+7+\n547Oc/fb3L3a3auHDlUzlJ15c3MLr6zZzHFjhhQ7lN122SmHsH5Linvn1Xd9soiIiEgvEKXryJvA\np3bjvecC48xsDNkE+zzggigXmpkBtwOL3f0Hu/HZkmfO8ux6QzVj9ityJLvv2Kp9mTxqH2Y9/iAf\nSTURO/hklZGIiIhIr9Zhom1mN7r7VWZ2H9kWeztw9w929sbunjazK4C/k23vd4e7LzKzy3LHbzWz\nA4E6sj26QzP7LNkOJUcBFwLPm9nC3Ft+xd0f3PVblNnL19O/IsZRI3vvsutdMTO+dMRbTHrkOuyx\nNNTe0LdqzUVERKTsdDaj/fvc1//b3TfPJcYPttt3a974dbIlJe09CfTiNcL7lmeWvcmU0ftSEYtS\nKdR71fAibmkCQjyTwlbUKtEWERGRXqvDRNvd5+S+PtK2z8wGAyPc/cUeiE26wYamFC+/sYn3Hjm8\n2KHssWDMVNJBgnQmhcUqiKnVn4iIiPRiURasecTMBpnZvsBC4G4z+37hQ5PuMHdFI+5w3MF990HI\nbZI1hBfdz89i5/PN/b6j2WwRERHp1aLUEuyXa633QeDX7j4F+H+FDUu6S8Nzj/PpiulM4uVih9It\nElXvIHbKVfyq/gCeb9hY7HBEREREOhQl0Y6b2VDgbGBGgeOR7lQ/hwteuoLPxv5A5W/+A+rnFDui\nbnHBcaMYWBnnZzOXFjsUERERkQ5FSbT/F3gCWOnuc8zsYGB5YcOS7pBa+gQxTxMjhEwKVtQWO6Ru\nMahfBRccP4rVLzzBhr9/p2T+ASEiIiKlpctE293vcfeJ7n5pbnuZu59V+NBkTy3pP5lW4oQWg1gC\nSujhwU+OWcevK65n4NPfhTvPVLItIiIivU6UhyG/nXsYMm5mfzezN8ws0sIzUlyzUgfzkdRXaDrx\n6pLrOb3f2jkkLDtb7yU0Wy8iIiKlI0rpyOm5hyHfB7wGTAC+VNCopFssrN/AG4MnMeDfvlhSSTYA\nVVOxeIK0B6QtXlKz9SIiIlIaIj0Mmft6BnCvu69nJytFSu/zXMNGJiX77mqQnUrWEEybwX37XMyl\ndh1bD5xS7IhEREREdhAl0f6bmb0AHAf8w8z2B1oKG5bsqfVbUqxc38RRI/cpdiiFk6zhoPd/lce2\nVDF94WvFjkZERERkB1EehvwCcCowxd1bga1ke2pLL/ZcwwYAJpVyog2ccMgQJgwfxJOPPYjPvFEP\nRYqIiEiv0eES7GZ2irs/YWZn5u3LP2VVIQOTPfNs/UbM4MiRJVo6kmNmXH3ERmpmXoM/msbilSX3\n4KeIiIj0TR0m2sB7yPbPPnsnxxyYXpCIpFs817CBsUMHMKCysx9xaTgp/hJuaYL8fuFKtEVERKTI\nOszC3P2a3NcLey4c6Q7uzrMNGzjl0AOKHUqPiB18Mq1PfI90ppUgVkGgDiQiIiLSC3Q53Wlmg4CP\nAlX557v75woXluyJ1zZuZd3mVOl2HGkvWUPLBffz87vuxEZO5bOazRYREZFeIEpdwYPAfOB5ICxs\nONIdnq0vjwch8w0YewKba/bhzlkrOGdDMwft07/YIYmIiEiZi9Leby93v9Ldf+7ut7e9Ch6Z7LZn\nGzZQETPGDx9Y7FB61CUnVuHAr2atKHYoIiIiIpES7d+a2SVmNjS3FPugXDmJ9FJNS2dx7eCHqFw9\nr9ih9KiR++7F6UccyIuz/0nLY99Xqz8REREpKnPvfJFHM7sM+C6wie0rQrq7jypwbLusurra6+rq\nih1GUYWvzqbljveRsDSxMmx1t2TeI4yYfh6VliYow/sXERGRwjOzee5e3dV5UWa0vwSMc/eR7p7M\nvXpdki1ZGxY/RgVpYvmt7srI2KaFJHKt/rwM719ERER6jyiJ9hLgrUIHIt3j5X5H0UoctxjEElBu\nre6qpkIsQdoDwqCi/O5fREREeo0oXUfeAhaY2aNAS9tOtffrnZ5OHcINrV/hnn/PUHHIyeVXNpGs\nwS6azh2/uZuX+03i+yOPxbq+SkRERKTbRW3v92ChA5Husfj1TWwYMpmKd76z2KEUTTD6OPZ+93D+\neN/zfHj5eo4/eEixQxIREZEy1GGibWZ7u/uWnbXyM7MRhQ1LdtdLr79VVv2zO/LBY0Zw48Mv8/OZ\ny5Roi4iISFF0VqO97SkyM3u43bEZhQlH9sSmra3Ur29mwnB1X+xXEePCd4ym8eUnefNv31arPxER\nEelxnSXa+aWtQzs5Jr3Ey69vAmBCmS1U05FLRq3hN4nr2Wf29+DOM5Vsi4iISI/qLNH2DsY725Ze\nYPHqbHOY8QdqRhtg8Buzs/3E1epPREREiqCzhyEPMLMryc5et43Jbbef4ZZeYPHrmxjcv4Lhg/sV\nO5TeoWoqFk+Qbk3hQZwKtfoTERGRHtRZov1LtifU+WOAXxUqINl9i1e/xfgDB2Kmyh4AkjUE02Zw\n35/vYcbGg7n1wCnonyAiIiLSUzpMtN392p4MRPZMGDovv76Jc6qTxQ6ld0nWMPL9h/D4bc/wp/kN\nfOS40cWOSERERMpElJUhpQ+ob2yiKZXRg5A7UTNmPyaNHMzttcsJQz1eICIiIj1DiXaJaHsQUq39\n3s7M+MTUg9nnzQUsve9/1H1EREREekSUlSGlD1i8ehOBwbgDNKO9M6fvs5L3VF5PxfNpWPxTmDa9\n/JanFxERkR7V5Yy2mV1hZoNy45+Z2Rwze3fhQ5NdsXj1W1Ttvzf9E7Fih9IrxVc+RQK1+hMREZGe\nE6V05FJ3f8vM/h0YBvwn8L3ChiW7KrZqLp+pnKGyiI5UTcXilaQJaCUOavUnIiIiBRYl0W57euwM\n4G53fzbidZjZaWb2spktMbOrd3J8vJk9bWYtZvb5XblWtmta+jQ/2Po13r/uDq2A2JFkDTZtOk+O\n/CTnt3yZ+r2PKHZEIiIiUuKiJMzPmtmDwPuAv5nZACKsDGlmMeDHwOnAROB8M5vY7rT1wJXADbtx\nreRsWPwoFaQJCEFlER1L1nDY2dfxLIdxx1PLix2NiIiIlLgoifYlwNeBGndvAvoBH49wXQ2wxN2X\nuXsKuAc4K/8Ed1/j7nOB1l29VrZ7qXISrcRxi0EsobKITgwf3J8zJx3EH+bWs7G5/a+diIiISPeJ\nkmgfR7Y7yXgzOwE4FOgf4boRQH3edkNuXxR7cm3ZmZ0+hGmZawjf+RV104jgE1MP5rDWxbx079dV\nZiMiIiIFE6W9X/4Kkf2AKcAC4JSCRLSLzOxS4FKAUaNGFTma4li6ZjMbh0wmdsrJxQ6lT5iYeYnf\n9fs2sWWt+MrbMf3jRERERAqgyxltdz897/Uu4ChgTYT3XgXkrwc+MrcvisjXuvtt7l7t7tVDhw6N\n+PalZcmazYw9YECxw+g7VtSS8FbihHhaNe0iIiJSGLu8MqS7rwAOj3DqXGCcmY0xswRwHjA94sfs\nybVlZWtrhpXrmzhEiXZ0VVMh1+ovRYz0qBOLHZGIiIiUoC5LR8zsJrZ3GQmAycCzXV3n7mkzuwL4\nOxAD7nD3RWZ2We74rWZ2IFAHDAJCM/ssMDHXt/tt1+767ZW+5eu2EDqa0d4VuVZ/y2b/javnDWJa\nY5KzRhc7KBERESk1UWq0X8gbp4H73P2JKG/u7g8CD7bbd2ve+HWyZSGRrpW3W7JmMwDjlGjvmmQN\nY0ccy6aVM/nJY0t5/1EHEQRW7KhERESkhESp0b4duBN4Cpid+yq9xJI1mwkMxuy/d7FD6XOCwPjU\nuw5h7zXzWHbf/6gDiYiIiHSrKKUjU4G7yT6MaMCBZnahuyvh7gWWrNlMcr+96FcRK3YofdL792vg\n9MrrqXg+jS/+qTqQiIiISLeJ8jDkTcAZ7n6iu58AvBe4pbBhSVRL1mxm7FCVjeyu+MqnSJAmRohr\nVU0RERHpRlES7YS7v9i24e6LgUThQpKo0pmQ5eu2MHaYEu3dVjUVy3UgaSWuVTVFRESk20RJtOeb\n2a1mdlLu9VOyC9ZIka1c30QqE2pGe0/kOpDMO/hTnL/1yzzLocWOSEREREpElET7MmAZ8MXcaxnw\nyUIGJdG0dRxRa789lKzh8HO/wZLKifzk8SXFjkZERERKRJcPQ7r7VjO7GfizuysL6UWWrM0m2lqs\nZs8NqIxz8YljePLRv7Lub/9g/yPerYciRUREZI90OaNtZu8Dngf+kds+2szuK3Rg0rUlazYzbFAl\ng/pVFDuUkvCJ0Wv5TeJ69p39PbjzTLX7ExERkT0SpXTkG8BxwAYAd18IjC1kUBLNkjWbGXfAwGKH\nUTIGvfEMlaYOJCIiItI9oiTare6+od0+3+mZ0mPcnQFr5nNx+GfNvHaXqqlYPEHa1YFERERE9lyU\nRHuxmZ0DBGY2xsxuAp4pcFzShXWLn+R2+yanvnabyhy6S7IGmzaD2uQnOb/lyyypnFjsiERERKQP\ni5JoXwFMAULgPiAFfLaQQUnXtvzrcSpIExCCyhy6T7KGo877BovjE/jhI68UOxoRERHpw7pMtN19\ni7t/yd0nu/vRuXFTTwQnHVtcOYlW4rjFIJZQmUM3GjKgkmknVDHjudf41xubih2OiIiI9FFdtvcz\ns7HA54Cq/PPd/d8LF5Z0pS4zljvDa/jtu9PYmKlqRdfNLp16MM/OepjFf3iUQz9wrr6/IiIissu6\nTLSBPwK3A78GMoUNR6JasW4LG4ZMJjj55GKHUpL2Xb+QO+Pfwta1Ev7q1wQXz1CyLSIiIrskSqId\nuvuPCh6J7JLl67Zw2IFq7VcwK2qJexqzkExbDbwSbREREdkFUR6GfMDMLjWzoWY2qO1V8MikQ+lM\nyMr1TVTtv3exQyldVVOxWIKQGCmPs3TvycWOSERERPqYKDPan8h9vTZvnwOjuj8ciWLVhmbSoTNG\niXbhJGtg2nRSS57g0if6Ufn8QH5xTLGDEhERkb6ky0Tb3ZM9EYhEt2zdFgAl2oWWrKFfsobjWcIj\nD89g9V8eZvik96iERERERCKJMqONmY0HJgL92va5+28LFZR0bkUu0a4aokS7J1wyai0fr7yeiro0\nLPwRTJuuZFtERES61GWNtpldA9wG3AqcDtwMfLjAcUknlq/bwsDKOPsPSBQ7lLKw12uzSJAmRohr\ncSARERGJKMrDkOcC7wJWu/uFwCRAU6lFtHzdFqr23xszK3Yo5aFqKhavJE1AymP46JOKHZGIiIj0\nAVES7WZ3zwBpMxsIvA6MLmxY0pkVb25RfXZPStZg06az6NBPc37LV3h4k379RUREpGtREu0FZrYP\ncAdQB8zJvaQIWtIZVjU2q7VfT0vWcPi5X2fj/pP53kMvkc6ExY5IREREerkuE213/6S7b3D3HwPv\nBT7p7hcVPjTZmfr1TYQOY/bfq9ihlJ14LOCLp41n8LoFLPr916Fe/94UERGRjkV5GPLhtrG7L3H3\n+fn7pGctX9cEwJj9BxQ5kvL07wNf5XeV13P4v36E33mmkm0RERHpUIeJtpklcitADjOzgXmrQo5E\ni9UUzfJ1mwEYo9Z+RWGvPkmCNHFCPK0OJCIiItKxzvpoXw58DjgAWAS0tbh4i2yrPymC5eua2Hev\nCgbvVVHsUMpTrgNJJt1CymM07V/DkGLHJCIiIr1Sh4m2u98E3GRmn3X3m3swJunEinXqOFJUuaXZ\nN734KJfW9mfUC4O4YUKxgxIREZHeKMoS7DebWQ1QlX++VoYsjuXrtnDCWM2hFlWyhn2SNUzOLOa2\n2mVcfEIVR4wYXOyoREREpJeJ8jDkr4D/A/4NmJp7acWOImhOZTho03Oc3fQHPYTXC1x+6lhO6b+c\n5+75Gr5ydrHDERERkV7G3L3zE8xeAia6e69vHFxdXe11dXXFDqNgVix4jGH3n0OlpQnilTBteraU\nQYqjfg7pX74fMiksniB28Qz9PERERMqAmc1z9+quzouyYM0iYOiehyR7KrV0JhWkCQgho44XRbei\nlpi3ErcQz6RoXTqz2BGJiIhIL9JljTYwGHjRzJ4BWtp2uvsHCxaV7NSLlUeSJE7MMlgsAVVTix1S\neauaisUS2SQ7jHHfm6O5oNgxiYiISK8RJdH+9u6+uZmdBtwCxIBfuPt32h233PEzgCbgYnefnzv2\n38AnAAeeBy5x9627G0spqMuM48/B17jr1FQ2yVaZQnHlOpDYilp+tmQYP5k/gBNO2UKVusKIiIgI\n0ZZgfwR4FqjIvZ7N7euUmcWAHwOnAxOB881sYrvTTgfG5V6XAj/NXTsCuBKodvcjyCbq50W8p5L1\n6ptNbBgyGaZepSS7t0jWwNSrOP9DH6YiZvz2j3+E2hv1sKqIiIhE6jryIWA+cCFwEVBnZv8R4b1r\ngCXuvszdU8A9wFntzjkLuMuzngH2MbPhuWNxoL+ZxYG9gNci3VEJW7m+iVH77VXsMGQnhg3qx/XH\nNvPfqz9P+Mi3QMuzi4iIlL0oD0N+DTjW3T/i7hcAxwFfj3DdCKA+b7sht6/Lc9x9FXADsBJYDWx0\n94cjfGbJSmdCVjU2K9Huxd43aCkJyz6s6npYVUREpOxFSbQDd38jb3tNxOt2m5ntS3a2ewxwELC3\nmX20g3MvNbM6M6tbu3ZtIcMqqtUbt5IOndFDlGj3VrGDT8ZiCdIekCauh1VFRETKXJSE+WEz+6uZ\nfTSX7M4A/h7hulVAMm97ZG5flHP+DVju7mvdvRX4M3DCzj7E3W9z92p3rx46tHS7EK5c3wRAUjPa\nvVeyhuDiGTwy/D85r+XLvBAcVuyIREREpIiiJNqfB+4kW3Ndkxt/IcJ1c4FxZjbGzBJkH2ac3u6c\n6cBFlnU82RKR1WRLRo43s71ynUneDSyOdEcl6tU3s4n26CHqaNGrJWs4/qL/5dW9juRXv/8D4Uw9\nGCkiIlKuumzv5+5uZo8BW4AQqPOulpPMXpc2syvIzn7HgDvcfZGZXZY7fivwINnWfkvItve7JHds\ntpn9kexDmGlgAXDbbtxfyVi5vomKmHHgoH7FDkW6MHivCm46oYXqJ74Cj6ZBq3iKiIiUpS4TbTO7\nBPgf4AnAgFvN7GvufmdX17r7g2ST6fx9t+aNHbi8g2uvA67r6jPKxcr1W0juuxexwIodikRwUvwl\nwrwHI21FrRJtERGRMhNlwZqrgWPcfS2AmQ0FniRbQiI9ZOX6JtVn9yE2ZipBPEE6nSLjcSpGnVTY\nJ4hFRESk14nyd/96YEPe9obcPukh7s6rbzap40hfkqzBps1g8fgrOb/lyzz4wmotZCMiIlJmosxo\nvww8bWb3k10O/QPAC2Z2JYC7/7CA8QmwsbmVTVvT6qHd1yRrOOLcYzlmw528e85/4kEGiyVUry0i\nIlImosxo1wN/AyqBfsBDZBeWGZp7SYG1dRxRot33mBmfOeQNKiyNeUYL2YiIiJSRKF1Hrm0bm9lA\nd99U2JCNErUPAAAgAElEQVSkvbYe2qNUOtInDRz/LtJP/4B0JoUHcSq0kI2IiEhZ6HBG28y+ambj\nc+OEmT0MNJjZG2Z2ao9FKNsTbc1o903JGuKXzOChAz7O+S1fYYGPK3ZEIiIi0gM6Kx25gGx9NsBF\nZEtH9gdOBb5d4Lgkz6tvbmHowEr2SkQpqZdeKVnD1I99m9UDj+Lnv72H1GPf14ORIiIiJa6zRDuV\ntzDNacDv3L3V3RcBFYUPTdqsXN+k2ewSMLh/Bbe9K8ONzV8j9sT1cOeZSrZFRERKWGeJdouZTTCz\nIWRnsR/OO9a/sGFJvpVvNjFaiXZJODz1HJWWJkZImNaDkSIiIqWss0T7KmA62eXRf+juywDM7Azg\n+R6ITYDUiqf5wJbfUxNfUuxQpDtUTcXilWQIaPEYKwdNKXZEIiIiUiC2vTqk76uurva6urpih9F9\n6ucQ3vl+wtYUxBLEL5mh/suloH4Om19+jE8/PYABlXFurHmLxCGn6GcrIiLSR5jZPHev7uo8PV3X\nm62ohXSKuIW4t2a3lYz1fckaBiRr+O+9H2bcQx8h9ngaf/JGTAvZiIiIlJQoC9ZIsVRNJQwSpD2A\nWAWo/3JJOSr9/LZ6bVe9toiISMlRot2bJWv47fgfcUt4Dn6hZjtLTrt67Vc2JaD2RnUiERERKREd\nlo6Y2ZmdXeju07s/HGlvdnosLw4+gKtGH1fsUKS7JWuwadPJLHmCW2c3ctmc/8Etg8USoDISERGR\nPq+zGu2zc1/3B04AHs9tnwLMItuRRAqsobGZkfuqm2LJStaQSNZwcct3qHgmjXmIZ1KY6vFFRET6\nvA5LR9z9Qne/EIgBE939LHc/Czi8s+uke61qbFKiXQb2PfxULJatx08RpzV5YrFDEhERkT0UJWEe\n6e6r8rZfA0YVKB7J05zKsG5zipH7arGakpesIXbxDBZPuJLzt36Zn81cis9UvbaIiEhfFqW93+Nm\n9lfgd7ntc9leRiIFtGpDE4BmtMtFsoYjz6vh7Ip7+cBzn8KXpbF4peq1RURE+qgoifblZOu123rL\n3QX8sWARyTb1jc2AEu1yc97QlYSWJsgt0x6oXltERKRP6jLR9uzSkX/IvaQHNWxLtFU6Uk5szFSC\neCWZdAspj/HSmwGTa2/M9lFXwi0iItJndJlom9lZwHeAgwDLvdzdBxU4trLX0NhEIhYwdEBlsUOR\nnpRr+xcuncld8zdy0YLrszPcKiMRERHpU6I8DHkjcI67D3b3Qe4+UEl2z2hobGbEvv0JAit2KNLT\nkjVUvPPzXDx5EIm2MpKMVo8UERHpS6Ik2m+4+/MFj0TeRj20pXLsKQTxBGkCWsIYL23U6pEiIiJ9\nRZSHIeea2W+A+4GWtp1aGbLwVjU2MXHisGKHIcWUrMGmzaD1lcf5Rd1GPjH3f1RGIiIi0kdESbSH\nACGQvyS7o5UhC0o9tGWbZA39kzX8J9+nYqa6kYiIiPQVUbqOXNgTgciO1ENb2us37hR81g1k0ilS\nHmPea85J6kYiIiLSa0XpOlIJXEx26fV+bfvd/dLChSXqoS1vkysjCZfO5I8vbOHDL36XcHF2URtT\nGYmIiEivE+VhyLuAKuB9wGzgEGBrAWMS1ENbOpDrRvKRowZs70aSbqF16cxiRyYiIiLtREm0D3X3\nLwOb3f124DRAU2cFph7a0pkgt6hNSIyUx/nZ3EY2/fO76kYiIiLSi0R5GLI193WDmU0A3gAOKFxI\nAuqhLV3ILWpjK2pZvqGCj9d9k4on02Se/gGxi2eojERERKQXiDKjfbuZ7QtcB/wd+BdwQ0GjEvXQ\nlq4la2DqVUzcp5V+QYY4IZ5O8eJDt+Iz1WtbRESk2KJ0HflZbvgYMKqw4Ugb9dCWyKqmYrEEnknh\nBBzS8ADhqvsI4glsmma3RUREiiXKjLb0sJZlT3NO871MCZYUOxTpC5I1MG06dupXiU/5KBUWEiMk\nk07x2hN3aCVJERGRIiloom1mp5nZy2a2xMyu3slxM7Mf5o4/Z2bH5B3bx8z+aGYvmdliM3tHIWPt\nNernUPHbD/C5+L184Ln/UoIk0eTKSOzoCwjilbjFyBAw5JV7CR/5Fn7nmfpdEhER6WFRHobcLWYW\nA34MvAdoILuU+3R3fzHvtNOBcbnXccBPc18BbgEecvcPm1kCKI8+dytqIZ0ibiHurdlt/a9/iapt\ndntFLbHGlQTz7yIgJJ1u4fXH72BkVa0WuBHpTP2c7J+7VVOz27s7TtZ033t1MfbRJ+EOvPok4egT\nc/uzY3cIXn2KcNSJhA7ByidJj8qeE3v1KdLJE3Gc+MqnaM2NYytnkUq+AxwqGmaRGnECqYOqia+q\nI9Ewi5aR78AdKhtm0TwiOwfWb9XTbI0wbhlxQm48i60Rxi0js+PKhj0fm20fp4ZXk1hdt/PzMBIN\ns0jl9nc07ugz2p/f2u5as+z3tW1/xU7O2TZO5sb12bFhVNTPojV5AumDjqVi9dxt2wAV9bNI58bx\nvHF2/4m5/U9tG1fUZ3/u2Tieets5hu1wfrz+qW2/P/GV2bEZxFY+Rabt9ypv3HZO2zgz6qTcOU8S\n5o0zo04iHHkssVVzCfKOBSufxLeNs7/HWO53evT2940dfHKv/XstyoI1s4E7gN+5+1u78N41wBJ3\nX5Z7n3uAs4D8RPss4C53d+CZ3Cz2cKAJOJnsQjm4ewpI7cJn911VUwmDBGEmRSxesf0PU5GokjWQ\nrCFePwee+z2eSZFxY+iSe8ksDbF4gkC129KXdFPCmg5DwmW1NI94B6m0E6x8ko3DjqM1E9Kv4Wla\nEvtwyLxvYZkUHmT/erQwTWhxDDBPE1oMgMAz28bmGUJi2QTAM2SsgnuGfIpz3/wJcW8lQ/a8GJkd\nxmkCDNthHJAh025M7pyMB7lrw9z+zseeGxshEBDk9pN3juWNg7xxRd44Tpzvt17IdRV3U0GafnnH\n+ueN99rFcf8ijRPE+UYH97Or48oI48wOn931uCLCONbuHuI7/Lx2Po51MA4ijG0PxnQxDgjJEOeb\nufsJSG/73Q3yvn+Wd80OYwuh9gbopQu3RZnRngZcAiw0s1nAL939kQjXjQDq87Yb2D5b3dk5I4A0\nsBb4pZlNAuYBn3H3LRE+t29L1vDb8T9i7fOP8N8Xfgzrhb800kd0NLvdmuL5v/6UCeOfIHHIKb3y\nDyYpAbuaHCdrSL/6DOmltWw68Hi2tmaI1T9FU2wQY+qyyW9o2b+yAk8TEsNziW2YS17bJ6bptrFn\ndkh2+hHQLzcesENiakBIYE4m04rhBAYe5o99p2NwzLPpgXsrEzc8TtxbcwlF7jzA8Fys3uE46Ghs\nO+633PsVfpzhMwctpnJddpGsnv3sAt3PiMVUru3L95PhyuEvUvlmX76HHe/n8mGLqGzcvfshk+q1\nFQBRuo68BHzJzL4CnAncZWYpsrPcP3L3DQWK6xjg0+4+28xuAa4Grm1/opldClwKMGpUaTRFmZ0e\ny4uDD+Cq0e3/XSKyi3Yyu+0EjF89neD1+0nPvAFO+zbxlkZUUiKd6iJxDkedRFNrhvTSmWyJD+LA\nWd/Ylhxvmw3OS44zeclxK3G+lb6Ia2J3UUGaQQQM4u3J7w4Jr3su4YQuE1mL9pe1E+AWw3Eslo3P\nwwzB28YGYZpY/jjYPo7HElSfdjE8dDVkUjscC/riOJZg+PHnwkMLIJPqHTHt6f0c18fvJ5bgoHec\nBw8t7Lv30O5+kieeDw89u1v3QyxBb60AiFSjbWYTyc5qvx94APgNcBLwKNmEeGdWAcm87ZG5fVHO\ncaDB3Wfn9v+RbKL9Nu5+G3AbQHV1tUe5n95OPbSl2+XNbldsbMDrfoURks604A9+PpvOxCuxXvq/\n3qQAOkmcM8tmsmn48WxqTsOKWjbaQMYv/DZBmGqXLGdT2yBX1lAB9CNk746SY5zAd5YQp/n4fs9T\n+dbbZ7OcAHLJb1uS650kuXvyl73FEthp34HmN7HuqNEeNrFHarR7ZKz76X3jUriH7ryfXvp3l2XL\nozs5wWwO2ZrpO4B73b0579h0dz+zg+viZBe3eTfZ5HkucIG7L8o7573AFcAZZMtKfujuNbljtcAn\n3P1lM/s6sLe7f6GzWKurq72urq7zO+4Dqr/1D94zcRjf/uBRxQ5FSlH9HLjzzNzstuGeIYaTJuBf\nB/0Hh4ybQOVYlZT0OV3MOLeMPIENTa2kl82kkYGMX3A9QZjaYWa5o1ILxwgIiZmTccMsmyyHWO5a\nJ8S2JcchAViQTZSDWPZImMbyZ6Daz0ad9p1tM8A7Pdb8Zs/9ZS8i0gUzm+fu1V2eFyHRPtTd/9Vu\n3yh3XxkhiDOAm4EYcIe7/6+ZXQbg7reamQH/B5xGNpm/xN3rctceDfwCSADLcscaO/u8Uki0m1MZ\nJnztIb7w/w7j8neNLXY4UqrakrL+Q/CHrsYzKdIe4O7ECAmDCl476/dUDdlbCUixdZRAJ2vwlbNp\nWTKTRgZwwFNfz5VqZBNnyy/P6Cxxxmgru+gocXYCCAJw7zhZ3tPkuLMuHfrdE5FepjsT7fnufkxX\n+3qDUki0l6zZxL/9YCa3nHc0Zx09otjhSDloS242NuDz7sQ8Q9oD7sm8i7PjT1JBKxZPYKd9d3vi\npMSne3QyCx0ur2VzbBADHr0G2iXQaSr4UeLjXN7yCypIv33GOVeqUfTEWUSkREVNtDus0TazQ4EJ\nwGAzyy8PGQT02/MQZWfqG7OVOarRlh6Te2CS+jnYwt9lH96KV3DMyH2Jr2rNdSppgb98Lls7q6Q7\nmi5qoN+yQQx+4tq3JdH5s9D9c4UYMfMdu13QyhnxOVSmsjXNYbs65rZSjc5qkfc4cc7f7mgsIlLm\nOnsY8nDgg8A+wNl5+zcBnyxkUOWsYVuiXR7r80gvkntgkhW1WNVUJgJ+51/xTCo3L5p9SC0/6W6b\n6Qy2ri+vpLuTJLp16Uw2MpB9a7+2k1KO7GN+Mc8woIMkOr9TRvsEOr+rxcRTL4SHFmWf0FfiLCLS\nK3WYaLv7fcB9ZnaSuz/ZgzGVtYbGJhKxgKEDKosdipSjttntHMsl3rG8Wm4zwHNJd7oF/+tVOE4Y\nVLB0yrWM6tfMXoe+M/sGfbGMoIMk2t3Z+soTNDKQYbO210M72dnn/AcJB5PfdWPXkuj2La86TKDb\nP6GvxFlEpNfprHTkKne/EfiQmX2w/XF3/1xBIytTDY3NjNi3P0F2BQSR4spLvC2X1MX6D4FtSbcR\neHae1jMpDplzHYbTUvtdAssmnWGQ4LV3XMfQ2Gb6j3tn9n17qqY3P2nu4GG7VPIENjalSS+byQYG\nMm7+/76tjV2aILscNSFD2yXRgbVrV7eTJHqnpRxRZ6G7SqCVOIuI9FqdlY4szX19oScCkSz10JZe\nKz+pGzYxu+JkLunOLsphEIZYrsSEXL/kMNPC8CevySbgM79L25LTITEsl8iGQYLl1dcwIHyL1uSJ\nJGIBA15/GqqmkogbFStn4aNPyvbLWPEk4eiTsouKrHiK5oPeQSoTErz6FBuHHUcqk6Ffw9NsCrb3\nfs4EFfz1oM9wxqpbiHnrDv2f3bcvjrI/lo3ffIcFUQwHa3uosAeTaBER6dO67Dqyw8nZdnx79dal\n0Euh64h6aEufk9cqsH0fZA/TuBmWS8Dzu2BksOyS1QZpzzaYM3yHNnR7Ms7vxJH2gDl2JDU8T7xd\nHLvcjWN3k2gRESkZe9x1JO+N7iK7qEwamAMMMbPvu/sP9jxMydecyrBuc0oPQkrf0m6mOz/ptBW1\nWF4Cnj/rGwti2QVzwjRBsD0Zj7JUdtTltAli2d7g8QQnnPaxncaxW904NBMtIiIRRFmC/Sh3f8vM\nLgD+AXwJqAOUaHezVRuaALX2kz6sfc1wBwl429hy446S8T1dTnuHpDnq8r5KokVEpJtESbQrcsup\nnwX81N1TZhYWOK6ypB7aUrI6SsDzx1GS4F0d7yxpVjcOERHpIVES7V8AK8k+FPmEmY0CNhc0qjKl\nHtpS1qImwbs6FhERKZKgqxPc/SZ3P8jd/92zT07WA6cWPrTyox7aIiIiIqUjysOQCeADQFW7868v\nUExlSz20RUREREpHlNKR+4CtwDwgU9hwypt6aIuIiIiUjiiJ9mh3P6LgkQirGpuYOHFYscMQERER\nkW7QZY028IyZTSx4JGVOPbRFRERESkuUGe3jgAVmtgRoAQxwdz+moJGVmTdfquVTsQc40s8CxhY7\nHBERERHZQ1ES7Q8UPIpyVz+H4Q+cy+fiKYKnHoBxM9SeTERERKSPi9LebykwFDgxN94AtBY6sLKy\nohbLpIhbiIWt2cU2RERERKRP6zLRNrNrgOuAa3K7+gG/LWRQZadqKpmggrQHEEtsX91ORERERPqs\nKA9Dfhg4A9gC4O6rgEGFDKrsJGu4ZcSN/LLyI9i06SobERERESkBURLtltyKkA5gZmqLUQC1Ww9m\n5rALlWSLiIiIlIgoifafzezHwGAzuwR4GPhlYcMqP6samxixjxarERERESkVXXYdcffvmtnpQAqY\nBPyvu/+t4JGVka2tbT20lWiLiIiIlIouE20zu97dvwL8bSf7pBs0NDYDaLEaERERkRISpXTktJ3s\ne293B1LOGhqbADSjLSIiIlJCOpzRNrNPApcBh5rZ/LxDA4H5O79KdodmtEVERERKT2elI38AHgG+\nDVydt3+Tu68paFRlpqGxmYqYccDAymKHIiIiIiLdpMNE290bgUbgbDMzsqtDxoG4mR3k7q/1UIwl\nryHXcSQIrNihiIiIiEg3ifIw5H8B3wTeBMLcbgcmFjCustLQ2KyyEREREZES02WiDVwFTHD3tYUO\nplw1NDbz7vEHFDsMEREREelGUbqONADrCx1Iucr20G5RxxERERGREhNlRnsJ8KiZ/QVoadvp7j8s\nWFRlZNWGXMeR/ZRoi4iIiJSSKIn26txrUIFjKUtq7SciIiJSmqIswX5tTwRSrrRYjYiIiEhp6mzB\nmvvIdhfZKXf/YEEiKjPbe2j3K3YoIiIiItKNOpvR/r89fXMzOw24BYgBv3D377Q7brnjZwBNwMXu\nPj/veAyoA1a5+/v2NJ7eqKGxmYP26U9MPbRFRERESkpnC9Y8sidvnEuSfwy8h2znkrlmNt3dX8w7\n7XRgXO51HPDT3Nc2nwEWU8L14W2L1YiIiIhIaYnS3m931QBL3H2Zu6eAe4Cz2p1zFnCXZz0D7GNm\nwwHMbCTwXuAXBYyx6FY1Nqs+W0RERKQEFTLRHgHU52035PZFPedm4ItsX41yp8zsUjOrM7O6tWv7\n1po6W1szrNnUoo4jIiIiIiWokIn2bjOz9wFr3H1eV+e6+23uXu3u1UOHDu2B6LrPa209tDWjLSIi\nIlJyCtl1ZBWQzNsemdsX5ZwPAWea2RlAP2CQmf3a3T/axWf2KRv/9RSfij3A+NZKsrcuIiIiIqWi\nsxnt/yP7MGMD2fKNu3OvNLAywnvPBcaZ2RgzSwDnAdPbnTMduMiyjgc2uvtqd/+yu49096rcdY+W\nWpJN/RyOfOQiPhe/l/H/+AjUzyl2RCIiIiLSjbrsOmJm33X36rb9ZnY/0GVW6O5pM7sC+DvZ9n53\nuPsiM7ssd/xW4EGyrf2WkG3vd8ke3EvfsqKWIEwRWIhnWmFFLSRrih2ViIiIiHSTKEuwDzCzKndf\nkdseBQyI8ubu/iDZZDp/3615Ywcu7+I9Hgcej/J5fUrVVNJWQeCtxGMJqJpa7IhEREREpBtFSbSv\nAmrN7GXAgLHAfxU0qnKQrOFr+1zP0ekXOO+cCzSbLSIiIlJiuky03f2vZnYoMDG360V3by5sWOXh\nsS1VZMYdy3nJScUORURERES6WZft/cysP9kVGv8z125vhJmdXvDISlxLOsMbb6mHtoiIiEipitJH\n+47ceSfltl8Dri9YRGXitQ1bAfXQFhERESlVURLtce5+PdAK4O5NZGu1ZQ80NDYBSrRFRERESlWU\nRDtlZv3ILV5jZmOAVEGjKgMNjblVIfdT6YiIiIhIKYrSdeSbwEPASDO7EzgF+ERBoyoDDY1NxAJj\n2MDKYociIiIiIgUQpevI38ysDjiBbMnIF9x9TcEjK3ENjc0MH9yPeCzK/1QQERERkb4mSteRh919\nrbs/4O73u/saM3u4J4IrZasam1WfLSIiIlLCOky0zSxhZoOAYWY20MwG5V4jya4OKXugobFZrf1E\nRERESlhnpSOXA58DDgAWsb3TyFvArR1dJF1rSWd4Y9NWzWiLiIiIlLAOE213vwm4ycy+6O7fyz9m\nZlEeopQOrN6wFXc0oy0iIiJSwqI8iXfeTvbN6e5Aysm21n6a0RYREREpWR3OTJvZAcBwoL+ZHcn2\n0pFBgKZi90DbYjUj9lGiLSIiIlKqOisBeS/wMWAk8JO8/ZuAawsZVKlraGwmFhjDB/crdigiIiIi\nUiCd1Wj/EvilmZ3j7n/owZhK3qoNzRw4SD20RUREREpZZ6Uj57v774DhZnZl++Pu/sOCRlbCGhqb\nVJ8tIiIiUuI6Kx3ZN/d1/54IpJw0NDZzwiH6toqIiIiUss5KR36S+6p67G6USoe8/pZ6aIuIiIiU\nui77YZvZKOAKoCr/fHf/YOHCKl2rNzbnemgr0RaR4mlasICmOXPZq+bY7HaRx3tNntzrYooat4hI\nR6IsPDMduAv4BxAWNpzS19ZDe4QSbZGS0FPJYf9jq8kQ0jR7NpXVU4hNmsjWBc+ydW4dsSmTCD0k\nVbcAO+YIQg/JzHuWcPJEQkJ83vNkjp5AxjMECxcTDhxAvx//BlpbIRYABplMh2OPBVgHYzAsk8GD\nAGw3xxVx1n/yA+z3s/ux1vSO79vJZ+/uOMo9RxpXxGm6/DyCt5pITzoMMyO+8KXs9/rwscQWLSX+\n7Etkjh4PQGzhS4RHT9hhbBjBwhcJJ0/M7V9MeHR2HCx8EZ98eHa8oG1suf0Ts/ezYBEckz3H5r8I\nxxye3T//BTjmiOz+BYtgcnbMghewY47Mjue3jQ0WPI8dcySGEc5/nmBK9hyf9zzBlKOIHTmRzPOL\n8fnPERxzFIaRmf8c8WMmAZCZ/yyxKZMwjPS8Z4lXH53dP+9ZKqZkx+l5z1KR299at5CK6sm5cxYS\nn3I0htE6byGJ3P7WugVUVE/GMFJ1C0gce8y2/ZXV2XFq3vZxS918Ko+dkt0/dz6Vxx6DYbTMnUe/\nY6sB2Fo3j37VU6g8+ihaFj63bTv/WJSxYWytq6Nfddv79vx4+z3s5DyDrXPn0S/3/djTcf+279/c\nuu3fyz0Z121/z+a5dfQ/tpp+R09i68Jnt23nH+ts3DJ3PgOOO67X/qM3SqKdcvcfFDySMtG8bBaf\niv2FQ7YOBE4udjgiZWt3EuQts2cTmzKJlkwLTXPmkB64F5mbfw6taYgZ2xLTIAADy4SEQXYJgrax\nmW3fbxBknDAAN9s2BghCyGTfZtsYIBZCOga//LeAS/4ZEs/seKyjMQHEcmM3CMPsdhhmsp/XydjD\nDN7BmOxd79E4TLXSMOOP7JtyAn/7eR199u6Oo9xzlHGmpZXKm+/GPJd/576/FoO7834+1snPpO1n\n7bmx5e23kG2/D+SNrZOfcyHGu/P71pvHpXA/pXAP3XU/FW6sT1Qy6pd39MpkO0qi/SMzuwb4O9DS\nttPdnytYVKWqfg6nPPMJ3hlvJXbfdBg0HZI1xY5KpM+LkjTvNXkyG+fNoXHWkzTvHcdvuR1aW7cl\nxWRCvF1SnE2W3578wvaE1bYlrLlz2DFRtNC37SfjgG8b51JzCME8tx2COZ2Og9A4f/VoEuEKzJ0g\nBHx7Uk7eeeCY7zh2DGIB7o7lMkTPhDsZG57JEOTN4u50HGYIglj2GxB5vP3aWDzOlHMvJXPTbZBO\n73Bs2zhKHBHGnsnscG+djzv6vmTHgWX/xZL9GeR/rwM+8eZEYuGLmIcd/hyKOs79vnU83vH36KI1\nY0mEr3Tx+9ZXxsZ5kf776c1j47zVo0iEr/bIPUT5c2nPxsa5q5MkwpXZPwt36Xr//+3de3CddZ3H\n8ff3eU7Skra0XEqF0jZcWirqSJlwUYRFQW4q5aKri6ul3sZBV12dcdXZnVVntOCyaFXEQURgQJAB\nd4moCwhCW0vblB5GsFVoS3rB0vv9tMm5fPeP50l7kjbJSZonz8nJ5zWTyXPOeZ7ke/pMk09//f5+\nPzyfj37WD9GgPQ34FHAlB1tHHA3H9l3rfIJSntBKUGyH1vkK2iK96ClE71q0kL0NIW233dElNBcp\n2cER40LGuP+yej76RBuZYueA7J1GKzuH4o4g3DnkcuCXTxRYDXewjtaEjkAZh/cgDKNzi8WBO66r\n44zrZrHxbzfj+TxW9lolx0FdHRO+/jWKO3ZWVa9z7u3vTL2OvhyH48aycc6h9yCoq+PE93+QjS/1\n7/5U23FQV8epM29g4/LaeT/Tj+DvTzUcR+/hxkF7D4n8HCs7DuvqePN1sw+8n75eb3V1B/5uVhtz\n955PMFsJvMXd23o8sQo0NTX50qVL0y6je+uW0Pbz95GhQJgZAbM0oi3DW7ch2mHLwufY3QDMvRvL\nFzq1YBTLRpg7jSrHXzfoehzAljNPYvzyDVjJcTMIAswdwoOhmLIf3JUcW4qB9UgnD1bjyM9Q1N09\nGKqTO7s71vupvuNaeA8D+X4G+2eamb3g7k29nldB0H4M+KS7bxmo4pJS9UEb+Pi3fsRHJ6zj8vd9\nUCFbho1Df3guYd+oOnb/11zI5ykFAeCVh2g7OMLsBgThgdBsZt0G4gOjjwMYkBVYRUSGn0qDdiWt\nI2OAv5rZYjr3aGt5vz7amcszb9+pvOvMq2DSaWmXIzLgOgXqeFR6+8gCmR/eB4UCpbgXw4reKURb\nWfsGZW0ajsU9sU4QZrpto6gkNDfMmMGIadO6Dcjlj/t6LCIicjiVBO3vJF7FMLFm214Aphw3KuVK\nROP24F0AABQYSURBVI7MoYH6WTbXtzHyxw9AvhCtokE0Kh0YeByogQN9zw5R2wYeB2c7ohBdrqdw\nrIAsIiKDpZKgvRDY7+5uZqcBZwBPJltWbWrdmgNgynENKVci0ne5bJZdixayub6NYO4v4kBtONGM\n9/ryQN11VDoTJe8kQ7SIiEi1qSRozwcuMrOxwDPAMuAjwMeTLKwWrd0ajWhPPlZBW6pT+Ui1vW06\nr8z7DVsWPsta28HZD71IWIjS84FAXbYyR38CtUK0iIjUskqCduDuOTP7BHCHu99sZi8mXVgtat2a\n44QxI2ior+SPXWRw5LJZ9i5ezJb6dkrx5iuFjHHPpQGznioyoQgnGAQefQx0oBYREalVFQVtMzsH\n+Cjw6fi5sIfzpRtrt+ZoVH+2pKR8tHp/YT+tz/2WtWzn1HufJSiUOk1OzBSd618/ifrS+mg1Dwui\nre/cFahFREQqVEnQ/jLwLeBxd3/ZzE4laieRPmrdupeLpo1PuwwZRjrC9f7Rdez83vehPU8xnqhY\nX4LTOo1Ug2VCcAjr6njL9Z9k4yuHXw5PgVpERKR3vQZtd3+GqDe74/Fq4KYki6pF+9qLbNrdxhT1\nZ0vCctksOxb9iXXBDkbd/hDki51Gq8snKvY0Ut3bcngiIiLSs16DtpmdTjSq3Vh+vrtfVsG1VwBz\niX6/3+XuN3d53eLXrwJywI3uvszMJgH3AROIBtrudPe5Fb6nqrR2W7ziyPFqHZGBUd4KUigVWP3H\nZlaxman3ziMoOA1xuO7YJdHCqOOr67a1vY1UK2CLiIj0TyWtI48APwfuB4qVfmEzC4HbgfcC64EW\nM2t29+Vlp10JTI0/zgPuiD8XgK/EoXsM8IKZPdXl2iGlNV5xRCPaciQ6wnU4bixvfHcOnm+naOA4\nmRJMO2TSYgjuhH3oqxYREZGBUUnQLrn7j/rxtc8FVsatJpjZQ8BMoDwszwTu82gf+EVmNs7MTnT3\nDcAGAHffbWYrgIldrh1S1sZraGsypPRVR7i2sUezcc4cPJ+nhB9oBQk4uAlMb60g5RSwRUREklVJ\n0H7MzD4D/A+dt2Df1ct1E4F1ZY/XE41W93bOROKQDWBmjcAMYHEFtVat1q17GXtUHWMb6tIuRYaA\njnBdHNPAllu+B+15SmV91m5EW5M7BJlMxa0gIiIiMngqCdqfij//R9lzDkwe+HI6M7PRwKPAl7oL\n9vE/Aj4DMHly4iX129ptORq1I6T0oCNc58eMZNvNt0K+c7iGOFxjZNQKIiIiUvUqWXVkUj+/9utA\n+bUnx89VdI6Z1RGF7Afc/dc91HcncCdAU1OT97PWxLVu3ctZk45JuwypMh3hum10PTtuua2bcB0d\nde2zViuIiIhIdatoi0Izmw6cCYzseM7df9nLZS3AVDM7hSg8fwS4ocs5zcDn4/7t84Cd7r4hXo3k\n58AKd7+tondSxdoLJV7fvo9rzpqYdilSBTrCdfvoEWy/5VbIF/odrkVERKR6VbK8378DlwHTgSeA\ny4EFQI9B290LZvb5+JoQuNvd/2Jmn41f/ynwO6Kl/VYSLe83O778AuBjwEtl271/w91/17e3Vx1e\n37GPksNkrTgybHWE69LRo9h88y2H9FyDwrWIiEitqWRE+8PAWcAyd/+YmZ0I3FPJF4+D8e+6PPfT\nsmMHPneY6xYQL6JQC9Z0LO2nFUeGpT3LXmDt7Nn4YcJ1oHAtIiJSsyoJ2vvcvWhmhXhN6zeAKQnX\nVVNyq5/npvC3nN52NHBh2uXIIHB3Vsz7X/76h0fZ/NoKLmjPE8YzCDRyLSIiMjxUErSzZjYOuBtY\nCuwCliRaVS1Zt4T3Lvk0l2XyhI80w6xmmHRu2lVJQjYtns+fn3qQln0ruPKxN5hWhKmBEWQyUDp0\n4xiFaxERkdrVY9COJyV+0913ALeb2RPA0e6+bFCqqwWt8wk8T2glKLZD63wF7RqSy2bZu3gxq09t\n4Ln1z3HZbQt5UxE+YBC4Ye5AwLjrrqfupJMUrkVERIaRHoO2u7uZPQW8NX68clCqqiWNF5InAxQI\nw3poVOvIUNcxsXHnyBL7/vt2yBcYEcKYs0ZQV4q2Py/fodHq6hh7zUwFbBERkWGmktaRF81shrtn\nE6+mBu0aP4Mb277BV6dv5vz3XKPR7CGqI1z72DFsmjPnkFVDglLAB069mj0v/wbP57VDo4iIiHQf\ntM0s4+4Fou3PW8xsFbCXaDUQd/ezB6nGIW3lpj0s82nsPucGmDQh7XKkH3LZLGtunI3n2yninVcN\nyYTRNuh1dRx77bUce+212qFRREREgJ5HtJcAZwNXD1ItNWnlxj0ATJswOuVKpK9yy7KseraZVX9b\nxOntbYQOgR3cBr27iY0K2CIiIgI9B20DcPdVg1RLTXpl425GZAJOPkab1QwFuWyWXYsWsrywnuPv\nfIyg4JwSEPdbQ1hXr5YQERERqUhPQXu8mX25uxdrYWv0wfDqpj2cNn40YVAz++/UnI7+611HQe7W\nH0K+yHiLJjUGDiEh467/oFYNERERkT7pKWiHwGhqaIfGNKzctIdzGo9JuwzpoiNch2PHsmHOdw7Z\ntdGxaGMZrRoiIiIi/dRT0N7g7t8etEpq0J62Aq/v2McNEyanXYqUyWWzrJ09m1J7z5Mb1SIiIiIi\nR6LXHm3pv5WboomQU0/QRMhqkMtmWT/vCV7960KmtEWTGw0IepncKCIiItIfPQXtSwatihr16sbd\nAEydMCblSoavXDbLnsWLWekbGf2ThwkKJU6OJzd6CTL1mtwoIiIiyeg2aLv7tsEspBat3LSH+kzA\n5GO14shg6ui/3j+6jh233Ab5AmM6TW4MGHf9hzS5UURERBJVyc6Q0k+vbNytFUcGWS6bpfXGG/H2\ndk1uFBERkVQpaCfo1U17OHuyVhwZDLteaOHPTz7Ia6+0MKO9ndCj54MwmuaoyY0iIiIy2BS0E5Jr\nL7B++z4+3DQp7VJqVi6bZeOCZ1jW9iqn3/Mc44rw9gAsE0IJTW4UERGRVCloJ+TAiiOaCDmgov7r\nJawLdpL54b1YocS0rpvLXKfNZURERCR9CtoJ2bpiATeFzbzVRwFvSrucmrC9ZRF//+Sn8XyBTFm4\nxgzUfy0iIiJVRkE7CeuW8K7nP8GFmTzhY80wrhkmnZt2VUPW6gW/J/vE/Wxc/TIX5QuEDg5YmDkQ\nrtUiIiIiItVGQTsJrfMJSnlCK0GxHVrnK2j30e5lS1n+1MMszq3gokdXckYRpgVGkMlAyTW5UURE\nRKqegnYCCpMvIO8Z6q1AGNZD44VplzQk5LJZNi34Iy+2reKUe55hVBHebUQ7ODpAwLjrrlf/tYiI\niAwJCtoJeDmYzrfbv8F3z97J9POv0mh2D3LZLHuXLGGNbWPEj+7HCiVO69R/HUAYqP9aREREhhwF\n7QS0vLaNZT6NYy+/BI4emXY5VeuNxfPY8unPQb7AiB4mN6pFRERERIYiBe0EtLRuY8pxDZygkH2I\n3cuWsuIPj/D0cW+wv+UFPqTJjSIiIlKjFLQHmLuzdM123jP9hLRLqRq5bJb1857gxbZVnHHfAhqK\ncFUIKz52AeHzLVAoanKjiIiI1BwF7QG2avNetu1t55zG4b31ei6bZfvzC1heWMeEnz1OUHDeXL6x\njIdcOu5cGu75HLklLQrXIiIiUnMUtAdYS+s2AM5pPDblSgZfLptl16KFrGYLY37yMFYo8aYeeq87\nwrUCtoiIiNQiBe0B1vLaNo4fXc8px49Ku5RBkctm2b14EWtsG0f9+EGsUGRMWbh2DNPERhERERmG\nFLQHWMuabTRNORYzS7uUxOSyWXYuWkgrWxn9k4exQpGjugnX6r0WERGR4UpBewC9sXM/67btY9Y7\nGtMuZcDlslm2/Ok5VvpGjr+zmaBQYrTCtYiIiEi3FLQH0KsvPM1N4WNc3DASODXtcvoll80emJxY\nKBV57dlmVvompt47j6DgjFe4FhEREamIgvZAWbeE8+bP5h2ZPOHvm+GE5iGzI2RHuLaxR7Npzhw8\nn6do4DhhCaaVT2jk4IRGhWsRERGR7iloD5Dtf3mGMaU8GStBsR1a51dd0C4frQbY/vwCNoS7GXn7\ng1i+QMnAShACAWDxR9dt0BWuRURERHqXaNA2syuAuUTZ7S53v7nL6xa/fhWQA25092WVXFtN3J0f\nrp7AV8kQWhEL66HxwtTq6Rqoc0uWsLchYO+tP4J8IRqtNgiKzggDSlGwBrAwOgrCDAZ4sahwLSIi\nItIPiQVtMwuB24H3AuuBFjNrdvflZaddCUyNP84D7gDOq/DaqvH4nzfwi7Un0HThz3jf0auikJ3Q\naPahITo+dtj6/Dx2jCwSzL0HKxQomoE5VnS8m9Fqx7BM1AoSdgnU5V9f4VpERESkb5Ic0T4XWOnu\nqwHM7CFgJlAelmcC97m7A4vMbJyZnQg0VnBtVdj6m/uw+3/FTY3n8g8nXs2Wlk00jK6DLd0E4gqO\n3Z1di5+HGW8hV8ixb8lS9r6tkf2F/Rz/bz+OQnRgUVgulijGw9FBCTDweIQ6wMsCNVFvNRCE4YHR\n6t76rBWwRURERPonyaA9EVhX9ng90ah1b+dMrPDa1OWe+CV///p3OaUIU15exWuP/4qg5JQCwAwr\nHjwODneMHTyfKCh3hOawy/GIEJ5/q3Fp3qOJiSWPrgEogXn5CHUA3nOgBo1Wi4iIiCRpyE+GNLPP\nAJ8BmDx58qB+79yzvycoQuAGRcfiEeQo+MajyRUfc8hxUAIOHAdc3vheWPEsFAqdQ3Q/A7UCtoiI\niEhykgzarwOTyh6fHD9XyTl1FVwLgLvfCdwJ0NTU5EdWct80XHwlYXMLXnICAwsz0XF3IbiPx9Yl\nQE/+8Cz48KxeW08UqEVERETSl2TQbgGmmtkpRCH5I8ANXc5pBj4f92CfB+x09w1mtrmCa1PXcPkN\nTL4tGtluuPhKOOHN/e7LruS4IzCXB+fujkVEREQkXRbNQ0zoi5tdBfyAaLGLu939O2b2WQB3/2m8\nvN+PgSuIlveb7e5Lu7u2t+/X1NTkS5cuTebNiIiIiIgAZvaCuzf1el6SQXuwKWiLiIiISNIqDdpB\nbyeIiIiIiEjfKWiLiIiIiCRAQVtEREREJAEK2iIiIiIiCVDQFhERERFJgIK2iIiIiEgCFLRFRERE\nRBKgoC0iIiIikgAFbRERERGRBChoi4iIiIgkoKa2YDezzcCaFL718cCWFL6vDC7d5+FB97n26R4P\nD7rPw0Na93mKu4/v7aSaCtppMbOllex3L0Ob7vPwoPtc+3SPhwfd5+Gh2u+zWkdERERERBKgoC0i\nIiIikgAF7YFxZ9oFyKDQfR4edJ9rn+7x8KD7PDxU9X1Wj7aIiIiISAI0oi0iIiIikgAF7SNgZleY\n2d/MbKWZfS3temTgmdkkM/ujmS03s7+Y2RfTrkmSY2ahmWXN7PG0a5FkmNk4M3vEzP5qZivM7B1p\n1yQDy8z+Nf55/bKZPWhmI9OuSY6cmd1tZpvM7OWy5441s6fM7NX48zFp1ng4Ctr9ZGYhcDtwJXAm\n8E9mdma6VUkCCsBX3P1M4Hzgc7rPNe2LwIq0i5BEzQX+z92nA29H97ummNlE4AtAk7u/FQiBj6Rb\nlQyQe4Arujz3NeBpd58KPB0/rioK2v13LrDS3Ve7ezvwEDAz5ZpkgLn7BndfFh/vJvqlPDHdqiQJ\nZnYy8D7grrRrkWSY2VjgIuDnAO7e7u470q1KEpABjjKzDNAA/D3lemQAuPs8YFuXp2cC98bH9wLX\nDGpRFVDQ7r+JwLqyx+tRAKtpZtYIzAAWp1uJJOQHwFeBUtqFSGJOATYDv4hbhO4ys1FpFyUDx91f\nB24F1gIbgJ3u/mS6VUmCJrj7hvj4DWBCmsUcjoK2SAXMbDTwKPAld9+Vdj0ysMzs/cAmd38h7Vok\nURngbOAOd58B7KUK/6tZ+i/u0Z1J9I+qk4BRZvbP6VYlg8GjZfSqbik9Be3+ex2YVPb45Pg5qTFm\nVkcUsh9w91+nXY8k4gLgajNrJWoDe4+Z3Z9uSZKA9cB6d+/4X6lHiIK31I5LgdfcfbO754FfA+9M\nuSZJzkYzOxEg/rwp5XoOoaDdfy3AVDM7xczqiSZbNKdckwwwMzOifs4V7n5b2vVIMtz96+5+srs3\nEv1dfsbdNQpWY9z9DWCdmZ0RP3UJsDzFkmTgrQXON7OG+Of3JWjCay1rBmbFx7OAx1Ks5bAyaRcw\nVLl7wcw+DzxBNKv5bnf/S8plycC7APgY8JKZvRg/9w13/12KNYlI//0L8EA8QLIamJ1yPTKA3H2x\nmT0CLCNaNSpLle8cKJUxsweBi4HjzWw98J/AzcDDZvZJYA3wj+lVeHjaGVJEREREJAFqHRERERER\nSYCCtoiIiIhIAhS0RUREREQSoKAtIiIiIpIABW0RERERkQRoeT8RkRpiZscBT8cP3wQUibYdB8i5\nuzbvEBEZJFreT0SkRpnZN4E97n5r2rWIiAxHah0RERkmzGxP/PliM3vOzB4zs9VmdrOZfdTMlpjZ\nS2Z2WnzeeDN71Mxa4o8L0n0HIiJDi4K2iMjw9Hbgs8CbiXY/nebu5wJ3Ee2eCDAX+L67nwNcH78m\nIiIVUo+2iMjw1OLuGwDMbBXwZPz8S8C74+NLgTPNrOOao81stLvvGdRKRUSGKAVtEZHhqa3suFT2\nuMTB3w0BcL677x/MwkREaoVaR0REpDtPcrCNBDM7K8VaRESGHAVtERHpzheAJjP7s5ktJ+rpFhGR\nCml5PxERERGRBGhEW0REREQkAQraIiIiIiIJUNAWEREREUmAgraIiIiISAIUtEVEREREEqCgLSIi\nIiKSAAVtEREREZEEKGiLiIiIiCTg/wF7+fpPAWxOqwAAAABJRU5ErkJggg==\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fig, ax = plt.subplots(figsize=(12,6))\n", "ax.plot(tlist, np.real(count1), tlist, np.real(count1_2), '.')\n", "ax.plot(tlist, np.real(count2), tlist, np.real(count2_2), '.')\n", "ax.set_xlabel('Time')\n", "ax.set_ylabel('Transmitted Intensity and Spontaneous Emission');" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### mcsolve" ] }, { "cell_type": "code", "execution_count": 11, "metadata": {}, "outputs": [], "source": [ "ntraj = 500 #number of Monte-Carlo trajectories\n", "# Hamiltonian\n", "ida = qeye(N)\n", "idatom = qeye(2)\n", "a = tensor(destroy(N),idatom)\n", "sm = tensor(ida,sigmam())\n", "H = (w0-wl)*sm.dag()*sm + (wc-wl)*a.dag()*a + 1j*g*(a.dag()*sm-sm.dag()*a) + E*(a.dag()+a)\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", "# intial state\n", "psi0=tensor(basis(N,0),basis(2,1))" ] }, { "cell_type": "code", "execution_count": 12, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "10.0%. Run time: 0.55s. Est. time left: 00:00:00:04\n", "20.0%. Run time: 1.05s. Est. time left: 00:00:00:04\n", "30.0%. Run time: 1.54s. Est. time left: 00:00:00:03\n", "40.0%. Run time: 2.05s. Est. time left: 00:00:00:03\n", "50.0%. Run time: 2.56s. Est. time left: 00:00:00:02\n", "60.0%. Run time: 3.24s. Est. time left: 00:00:00:02\n", "70.0%. Run time: 3.76s. Est. time left: 00:00:00:01\n", "80.0%. Run time: 4.26s. Est. time left: 00:00:00:01\n", "90.0%. Run time: 4.74s. Est. time left: 00:00:00:00\n", "100.0%. Run time: 5.22s. Est. time left: 00:00:00:00\n", "Total run time: 5.26s\n" ] } ], "source": [ "avg = mcsolve(H, psi0, tlist, [C1,C2],[C1dC1,C2dC2], ntraj)" ] }, { "cell_type": "code", "execution_count": 13, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAtoAAAF3CAYAAACbhOyeAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3Xdc1uX+x/HXxRYBB+BWwIkouPfWPJZaNmwvrU7HlnU6\nebJ9Rp3qV+d0Om2PDbOyzHZZlmW5TST3lqEoyhJkyL5+f4AcVBQsbm64eT8fDx5xf+f7viH83Nd9\nDWOtRUREREREapabswOIiIiIiLgiFdoiIiIiIg6gQltERERExAFUaIuIiIiIOIAKbRERERERB1Ch\nLSIiIiLiACq0RUREREQcQIW2iIiIiIgDqNAWEREREXEAFdoiIiIiIg7g4ewANSkoKMiGhoY6O4aI\niIiIuLANGzakWmuDqzrOpQrt0NBQoqOjnR1DRERERFyYMSahOsep64iIiIiIiAOo0BYRERERcQCH\nFtrGmPONMbuMMXuNMbMr2R9ujFljjMk3xtxXyX53Y8wvxpgvHZlTRERERKSmOayPtjHGHXgJGA8k\nAuuNMZ9ba7dXOCwdmAlcfIbL3A3sAAIclVNERERcS2FhIYmJieTl5Tk7itRzPj4+tGvXDk9Pz191\nviMHQw4E9lprYwGMMe8DU4DyQttamwwkG2MmnXqyMaYdMAl4ArjXgTlFRETEhSQmJuLv709oaCjG\nGGfHkXrKWktaWhqJiYmEhYX9qms4sutIW+BAhceJZduq69/An4GSmgwlIiIiri0vL4/AwEAV2fKb\nGGMIDAz8TZ+M1MnBkMaYyUCytXZDNY691RgTbYyJTklJqYV0IiIiUtepyJaa8Ft/jxxZaB8E2ld4\n3K5sW3UMAy4yxsQD7wNjjTHvVHagtXaOtba/tbZ/cHCV84aLiIiIONzhw4e56qqr6NSpE/369WPi\nxIns3r37nK9z6NAhpk6dCsDGjRtZvHhxTUf9zX788UdWr17t7Bh1kiML7fVAF2NMmDHGC7gK+Lw6\nJ1prH7DWtrPWhpad94O19jrHRRURERGpGdZaLrnkEkaPHs2+ffvYsGEDTz75JEeOHDnna7Vp04ZF\nixYBKrTrI4cV2tbaIuBOYAmlM4cstNZuM8bMMMbMADDGtDLGJFI62PFhY0yiMUYzjIiIiEi9tWzZ\nMjw9PZkxY0b5tl69ejFixAiys7MZN24cffv2JTIyks8++wyA2bNn89JLL5Uf/5e//IVnn32W+Ph4\nevbsSUFBAY8++igffPABvXv35oMPPqBLly6c6DZbUlJC586dObUbbXp6OhdffDFRUVEMHjyYzZs3\nn3T9E3r27El8fDzx8fGEh4dz7bXX0r17d6ZOnUpubi5QugJ3amoqANHR0YwePZr4+HheffVVnnvu\nOXr37s2KFSsc8IrWXw5dgt1auxhYfMq2Vyt8f5jSLiVnu8aPwI8OiCciIiIu7q9fbGP7oWM1es2I\nNgE8dmGPM+7funUr/fr1q3Sfj48Pn3zyCQEBAaSmpjJ48GAuuugirrzySu655x7uuOMOABYuXMiS\nJUsoLi4GwMvLi7/97W9ER0fz4osvArBz507effdd7rnnHpYuXUqvXr04tRvtY489Rp8+ffj000/5\n4YcfuOGGG9i4ceNZn9+uXbt4/fXXGTZsGDfddBMvv/wy99132nInQGnxPWPGDPz8/M54TENWJwdD\nSt2WnV/EhoR01uxLY/nuFDYnZmCtdXYsERGROs9ay4MPPkhUVBTnnXceBw8e5MiRI/Tp04fk5GQO\nHTrEpk2baNasGe3btz/rtW666SbefvttAN544w2mT59+2jErV67k+uuvB2Ds2LGkpaVx7NjZ33i0\nb9+eYcOGAXDdddexcuXKX/NUBQe3aItr2X0ki/lrEvg4JpGcguKT9vVoE8DNw8OYHNUGLw+9fxMR\nkbrhbC3PjtKjR4/yftWnevfdd0lJSWHDhg14enoSGhpaPn3c5ZdfzqJFizh8+DBXXnlllfdp3749\nLVu25IcffuDnn3/m3XffrXZGDw8PSkr+N4NyxSnsTp1p48TjiudoMaDqUUUkVSosLuHPizbxu+eW\n80H0ASb0bMXcG/qz4PeDWTRjCP+4JJL8ohLuXbiJMc/+yC/7jzo7soiIiNOMHTuW/Px85syZU75t\n8+bNrFixgszMTFq0aIGnpyfLli0jISGh/Jgrr7yS999/n0WLFnH55Zefdl1/f3+ysrJO2nbLLbdw\n3XXXcfnll+Pu7n7aOSNGjCgvwH/88UeCgoIICAggNDSUmJgYAGJiYoiLiys/Z//+/axZswaA9957\nj+HDhwOl3UQ2bCidefmjjz46ay4ppUJbzionv4ib50WzMDqRP4zqyNoHxvGvK3pzXkRLhnQKpH9o\nc64Z1IHv/jiSt6YPwM0NrnxtLR9GH6j64iIiIi7IGMMnn3zC0qVL6dSpEz169OCBBx6gVatWXHvt\ntURHRxMZGcnbb79NeHh4+Xk9evQgKyuLtm3b0rp169OuO2bMGLZv314+GBLgoosuIjs7u9JuI1A6\n6HHDhg1ERUUxe/Zs5s2bB8Bll11Geno6PXr04MUXX6Rr167l53Tr1o2XXnqJ7t27c/ToUW677Tag\ntL/33XffTf/+/U8q6i+88EI++eQTDYashHGlvrX9+/e30dHRzo7hMtKy87nprfVsOZjJk5dGcuWA\nDlWeczSngLsW/MLKvalMGxrKo5MjcHPTogEiIlJ7duzYQffu3Z0do1ZER0fzxz/+scYK3Pj4eCZP\nnszWrVtr5HquoLLfJ2PMBmtt/6rOVR9tqVRRcQm3vB3NzsNZzLm+P+dFtKzWec0ae/HW9AE8+fVO\nXl9Z+jHUYxdGaIUuERGRGvbUU0/xyiuvnFPfbKldKrSlUq8tj+WX/Rn85+o+1S6yT/Bwd+ORyREY\nYO7KOIL9vbljTGfHBBUREWmgZs+ezezZs2v0mqGhoWrNrkEqtOU02w8d499LdzMpqjUX9Wrzq6/z\n4MTupGbn88ySXQT7e3NF/7NPUyQiIiLiSlRoy0kKikq4d+FGmjTy4u9Tev6ma7m5Gf5vai/Scgp4\n4OMthAY2ZmBY8xpKKiIiIlK3adYROcmLP+xh5+Esnro0kuaNvX7z9bw83Hjlun60a9aIe97/hczc\nwhpIKSIiIlL3qdCWcqnZ+cxZEctFvdqcc7/ss/Hz9uD5q/qQnJXPg59u0SqSIiIi0iCo0JZyc1fE\nkV9Uwt3ndanxa/du35R7f9eVrzYn8eGGxBq/voiISF3yxBNP0KNHD6Kioujduzfr1q2r0evHx8fz\n3nvv1eg166pDhw4xderUcz7vlltuYfv27Q5IVH3qoy0ApOcU8PaaeC6MakOnYD+H3OMPIzuxYncq\nf/l8G4PDAukQ6OuQ+4iIiDjTmjVr+PLLL4mJicHb25vU1FQKCgpq9B4nCu1rrrmmRq9bF7Vp0+aM\nS9qfzdy5cx2Q5tyoRVsAeGNlHMcLi7lzrOOm4XN3M/zryl64GcNjn29VFxIREXFJSUlJBAUF4e3t\nDUBQUBBt2pTO4hUaGsqf//xnIiMjGThwIHv37gVKC+exY8cSFRXFuHHj2L9/PwDTpk1j5syZDB06\nlI4dO5YXnLNnz2bFihX07t2b5557jvj4eEaMGEHfvn3p27cvq1evBkqXXR89ejRTp04lPDyca6+9\ntvzf3w0bNjBq1Cj69evHhAkTSEpKAmDjxo0MHjyYqKgoLrnkEo4ePQrA6NGjObEwYGpqKqGhoQBs\n27aNgQMH0rt3b6KiotizZ89pr8m3337LkCFD6Nu3L5dffjnZ2dnlr8cDDzxA79696d+/PzExMUyY\nMIFOnTrx6quvlr82PXv2POO9cnJymDRpEr169aJnz57lq2ZWzLtgwQIiIyPp2bMn999/f3kuPz8/\nHnroIXr16sXgwYM5cuTIb/jJV8Ja6zJf/fr1s3LuMnIKbM9Hv7G3vRNdK/ebuyLWhtz/pf16y6Fa\nuZ+IiDQs27dvP3nDGxNP/1o3p3Rffk7l+2PeKd2fnXr6vipkZWXZXr162S5dutjbbrvN/vjjj+X7\nQkJC7OOPP26ttXbevHl20qRJ1lprJ0+ebN966y1rrbWvv/66nTJlirXW2htvvNFOnTrVFhcX223b\nttlOnTpZa61dtmxZ+bnWWpuTk2OPHz9urbV29+7d9kRNtGzZMhsQEGAPHDhgi4uL7eDBg+2KFSts\nQUGBHTJkiE1OTrbWWvv+++/b6dOnW2utjYyMLM/8yCOP2Lvvvttaa+2oUaPs+vXrrbXWpqSk2JCQ\nEGuttXfeead9553S1ys/P9/m5uae9HqkpKTYESNG2OzsbGuttU899ZT961//Wv56vPzyy9Zaa++5\n5x4bGRlpjx07ZpOTk22LFi2stdbGxcXZHj16nPFeixYtsrfcckv5/TIyMk7Ke/DgQdu+fXubnJxs\nCwsL7ZgxY+wnn3xirbUWsJ9//rm11tpZs2bZv//976f9PE/7fSo9L9pWozZVi7bw5uo4svKLuHNM\nzffNrsyNQ0IIb+XPX7/YTk5+Ua3cU0REpLb4+fmxYcMG5syZQ3BwMFdeeSVvvfVW+f6rr766/L9r\n1qwBSrubnOgGcv3117Ny5cry4y+++GLc3NyIiIg4Y4trYWEhv//974mMjOTyyy8/qW/ywIEDadeu\nHW5ubvTu3Zv4+Hh27drF1q1bGT9+PL179+bxxx8nMTGRzMxMMjIyGDVqFAA33ngjy5cvP+vzHTJk\nCP/4xz94+umnSUhIoFGjRiftX7t2Ldu3b2fYsGH07t2befPmkZCQUL7/oosuAiAyMpJBgwbh7+9P\ncHAw3t7eZGRkVHmvyMhIvvvuO+6//35WrFhBkyZNTjpn/fr1jB49muDgYDw8PLj22mvLn5OXlxeT\nJ08GoF+/fsTHx5/1uZ4r9dFu4AqLS3hnbQLjwlsQ0SagVu7p4e7GE5f05LJX1vCf7/fwwMTutXJf\nERFpoKZ/deZ9Xr5n39848Oz7z8Dd3Z3Ro0czevRoIiMjmTdvHtOmTQPAGFN+XMXvz+REFxTgjN0u\nn3vuOVq2bMmmTZsoKSnBx8en0vPd3d0pKirCWkuPHj3KC/0TMjMzz5jDw8ODkpISAPLy8sq3X3PN\nNQwaNIivvvqKiRMn8tprrzF27NiTMo8fP54FCxac9fm5ubmdlNXNzY2iopMb5M50r5iYGBYvXszD\nDz/MuHHjePTRR8/4PCry9PQs/xmceG1qklq0G7jlu1NIzS7gqoEdavW+/UKac2X/9ry+Mo7dR7Jq\n9d4iIiKOtGvXrpP6KW/cuJGQkJDyxyf6EH/wwQcMGTIEgKFDh/L+++8D8O677zJixIiz3sPf35+s\nrP/9+5mZmUnr1q1xc3Nj/vz5FBcXn/X8bt26kZKSUl5oFxYWsm3bNpo0aUKzZs1YsWIFAPPnzy9v\n3Q4NDWXDhg0AJw1OjI2NpWPHjsycOZMpU6awefPmk+41ePBgVq1aVd4fPScnh927d58135lUdq9D\nhw7h6+vLddddx6xZs4iJiTnpnIEDB/LTTz+RmppKcXExCxYsKH9OjqYW7Qbuo5hEAht7MbpbcK3f\n+/4Lwlm8NYl/LN7BW9MH1vr9RUREHCE7O5u77rqLjIwMPDw86Ny5M3PmzCnff/ToUaKiovD29i5v\n5X3hhReYPn06zzzzDMHBwbz55ptnvUdUVBTu7u706tWLadOmcfvtt3PZZZfx9ttvc/7559O4ceOz\nnu/l5cWiRYuYOXMmmZmZFBUVcc8999CjRw/mzZvHjBkzyM3NpWPHjuVZ7rvvPq644grmzJnDpEmT\nyq+1cOFC5s+fj6enJ61ateLBBx886V7BwcG89dZbXH311eTn5wPw+OOP07Vr1+q/qGe51/r165k1\naxZubm54enryyiuvnHRO69ateeqppxgzZgzWWiZNmsSUKVPO+d6/hjnTRxD1Uf/+/e2J0aVStYzc\nAgY+8T3XDu7AYxf2cEqGOcv38Y/FO5l/80BGdKn9Yl9ERFzPjh076N69bnZLDA0NJTo6mqCgIGdH\nkWqq7PfJGLPBWtu/qnPVdaQB+2JzEgXFJVzWt53TMtw4NJR2zRrxxFc7KC5xnTd9IiIiIiq0G7CP\nNiQS3sqfHrU0CLIy3h7u3H9+ODsPZ/GRVowUEREXFx8fr9bsBkSFdgO1LyWbjQcyuKxvu2qNeHak\nyVGt6d2+Kc9+u4vcAk33JyIiIq5BhXYD9dGGRNzdDFP6tHF2FIwxPDK5O8lZ+by1Ot7ZcURExAW4\n0hg0cZ7f+nukQrsBstby2cZDjOgSRAt/n6pPqAX9Qpozulsw/10eq0VsRETkN/Hx8SEtLU3Ftvwm\n1lrS0tJOmpP8XGl6vwZo15EsDmYcZ+a4zs6OcpK7x3XhkpdXM39tAjNGdXJ2HBERqafatWtHYmIi\nKSkpzo4i9ZyPjw/t2v36SSNUaDdA3+9IBmBMtxZOTnKyPh2aMbJraav2DUNC8PXSr6eIiJw7T09P\nwsLCnB1DRF1HGqIfdiYT1a4JLQLqRreRiu4e15m0nALeWZvg7CgiIiIiv4kK7QYmPaeAmP1HGRte\nt1qzT+gX0pzhnYOYszyW4wVnXz5WREREpC5Tod3A/LgrGWthXHhLZ0c5o7vP60JqdgHvrlOrtoiI\niNRfKrQbmO93JtPC39upi9RUZUBoc4Z2CuTVn9SqLSIiIvWXCu0GpLC4hOW7Uhgb3gI3N+cuUlOV\nu8d1ITU7nwU/73d2FBEREZFfRYV2A7I+Pp2s/KI62z+7okEdAxncsTmv/rSPvEK1aouIiEj9o0K7\nAflhRzJeHm4M6xzk7CjVcve4riRn5fO+WrVFRESkHnJooW2MOd8Ys8sYs9cYM7uS/eHGmDXGmHxj\nzH0Vtrc3xiwzxmw3xmwzxtztyJwNxQ+7khnSMZDG3vVjfuohnQIZGNacV9SqLSIiIvWQwwptY4w7\n8BJwARABXG2MiTjlsHRgJvDsKduLgD9ZayOAwcAdlZwr5+DIsTxiU3IY0aV+tGafcM+4Lhw5ls+H\n0QecHUVERETknDiyRXsgsNdaG2utLQDeB6ZUPMBam2ytXQ8UnrI9yVobU/Z9FrADaOvArC5vbWwa\nAIPCAp2c5NwM6RRInw5NmbsyjuIS6+w4IiIiItXmyEK7LVCxGTKRX1EsG2NCgT7AuhpJ1UD9HJeO\nv7cHEXV4Wr/KGGO4dURHEtJy+XbbYWfHEREREam2Oj0Y0hjjB3wE3GOtPXaGY241xkQbY6JTUlJq\nN2A9si4unf6hzXCv49P6VeZ3PVoREujLa8tjsVat2iIiIlI/OLLQPgi0r/C4Xdm2ajHGeFJaZL9r\nrf34TMdZa+dYa/tba/sHBwf/6rCuLDU7n73J2QysZ91GTnB3M9wyPIyNBzKITjjq7DgiIiIi1eLI\nQns90MUYE2aM8QKuAj6vzonGGAO8Duyw1v7LgRkbhJ/j0gEY1LG5k5P8elP7taeZryev/RTr7Cgi\nIiIi1eKwQttaWwTcCSyhdDDjQmvtNmPMDGPMDABjTCtjTCJwL/CwMSbRGBMADAOuB8YaYzaWfU10\nVFZXty42DV8vdyLbNnF2lF+tkZc71w8JZemOI+xLyXZ2HBEREZEqOXRCZWvtYmDxKdterfD9YUq7\nlJxqJVD/OhPXUevi0ukX0gxP9zrdJb9KNwwJ4bWf9jF3RRxPXhrp7DgiIiIiZ1W/Ky+pUkZuAbuO\nZDEwtP52GzkhyM+by/q146OYRFKy8p0dR0REROSsVGi7uJ/j0rEWBnWsnwMhT3XL8DAKi0uYvybe\n2VFEREREzkqFtotbF5eOt4cbvdrX3/7ZFXUM9mN895a8vTaB3IIiZ8cREREROSMV2i7u57h0+nRo\nireHu7Oj1JhbR3YkI7eQRRsSnR1FRERE5IxUaLuw3IIith3KdIn+2RX1C2lWuiz7Ci3LLiIiInWX\nCm0XtvXgMUos9O7Q1NlRapQxhj+M7Mj+9FyWaFl2ERERqaNUaLuwTQcyAIhq51qFNsD4iFaEall2\nERERqcNUaLuwTYkZtG3aiCA/b2dHqXHuboabR3Rk04EM1sdrWXYRERGpe1Rou7BNiRkuM9tIZab2\nbUfzxl7MWb7P2VFERERETqNC20Wl5xRwIP04vVyw28gJjbzcuX5wCEt3JLM3Wcuyi4iISN2iQttF\nbUp03f7ZFV0/JARvDzfmroh1dhQRERGRk6jQdlGbD2RiDES2c92uI/C/Zdk/jjlIclaes+OIiIiI\nlFOh7aI2JWbQOdgPP28PZ0dxuFuGh1FYUsL8NQnOjiIiIiJSToW2C7LWsjkxg17tXbvbyAknlmWf\nr2XZRUREpA5Roe2CDmYcJzW7gF4u3m2koj+MKl2W/cNoLcsuIiIidYMKbRe0OTETcP2BkBX1C2lO\n3w5NmbsylqLiEmfHEREREVGh7Yo2HcjAy92N8Nb+zo5Sq24d2YkD6cdZsu2Is6OIiIiIqNB2RZsS\nM+je2h9vD3dnR6lV4yNaEhroy5zl+7Qsu4iIiDidCm0XU1Ji2ZKY2WAGQlbk7ma4ZURHNiVm8nNc\nurPjiIiISAOnQtvFJKTnklNQTM82DWcgZEVT+5Uuy/7aci1gIyIiIs6lQtvF7Eg6BkD31gFOTuIc\nPp7u3DgklB92JrPz8DFnxxEREZEGTIW2i9mZdAw3A11a+jk7itPcODQEXy93XvtJrdoiIiLiPCq0\nXcz2pCw6Bvvh49mwBkJW1NTXi2sGduDzTYc4kJ7r7DgiIiLSQKnQdjE7Dx8jvFXDmtavMjePCMPN\nwH9XqFVbREREnEOFtgs5lldI4tHjDbZ/dkWtmzTikj5t+WD9AVKz850dR0RERBogFdouZNfhLAC6\nN7CFas7k1pGdKCgu4c1Vcc6OIiIiIg2QCm0XsrOBzzhyqs4t/Di/RyveXp1AZm6hs+OIiIhIA6NC\n24VsT8qiSSNPWgX4ODtKnXHX2C5k5Rfxulq1RUREpJap0HYhOw8fo3trf4wxzo5SZ0S0CeD8Hq14\nc1UcmcfVqi0iIiK1R4W2iygpsew6nEV4K3UbOdXMcV3IyivijZVq1RYREZHao0LbRexPzyW3oJgI\n9c8+zYlW7TfUqi0iIiK1SIW2izix3Hi4ZhyplFq1RUREpLap0HYR25OycDPQtaUK7cpEtAlgQo+W\nvLkqjuz8ImfHERERkQZAhbaL2Jl0jLCgxg166fWqzBjViWN5RXyw/oCzo4iIiEgDoELbRew4fIxw\n9c8+qz4dmjEorDmvr4ilsLjE2XFERETExTm00DbGnG+M2WWM2WuMmV3J/nBjzBpjTL4x5r5zOVf+\nJzu/iAPpx+neSt1GqjJjVCcOZebx5eZDzo4iIiIiLs5hhbYxxh14CbgAiACuNsZEnHJYOjATePZX\nnCtl9iZnA+qfXR2juwXTraU/r/0Ui7XW2XFERETEhTmyRXsgsNdaG2utLQDeB6ZUPMBam2ytXQ+c\nOudalefK/+w5kgVAFxXaVTLGcOvIjuw8nMVPu1OcHUdERERcmCML7bZAxVFniWXbHH1ug7M3JRsv\ndzfaN2vk7Cj1woW92tC6iQ+v/RTr7CgiIiLiwur9YEhjzK3GmGhjTHRKSsNsodyXnE1YUGM83Ov9\nj7NWeHm4cfPwMNbEprHpQIaz44iIiIiLcmRldhBoX+Fxu7JtNXqutXaOtba/tbZ/cHDwrwpa3+1J\nzqZzCz9nx6hXrhrYAX8fD+YsV6u2iIiIOIYjC+31QBdjTJgxxgu4Cvi8Fs5tUPIKizmQnksnFdrn\nxM/bg+sHh/D11iTiU3OcHUdERERckMMKbWttEXAnsATYASy01m4zxswwxswAMMa0MsYkAvcCDxtj\nEo0xAWc611FZ67O41BxKLHRRoX3Opg0LxcPNjbkr1aotIiIiNc/DkRe31i4GFp+y7dUK3x+mtFtI\ntc6V052Y2k9dR85dC38fLuvXlg+jE7nnvK4E+Xk7O5KIiIi4EI2eq+f2JGfjZiAsqLGzo9RLt4zo\nSEFxCW+vjnd2FBEREXExKrTruX3J2bRv7ouPp7uzo9RLnYL9GN+9JfPWJJCTX+TsOCIiIuJCVGjX\nc3uTs9U/+zf6w6hOZB4vZGH0gaoPFhEREakmFdr1WFFxCXGpOZpx5DfqF9KMAaHNmLsijsLiEmfH\nERERERehQrse25+eS0FxCZ2DVWj/Vn8Y2YmDGcdZvCXJ2VFERETERajQrsc040jNGRvegs4t/Hj1\np1istc6OIyIiIi5AhXY9tjeltNBW15Hfzs3NcOvIjuxIOsaKPanOjiMiIiIuQIV2Pbb3SDatAnwI\n8PF0dhSXMKV3G1oGePPqT/ucHUVERERcgArtemxvSra6jdQgbw93bhnekdX70lgbm+bsOCIiIlLP\nqdCup6y17EtWoV3Trh8SQqsAH576eqf6aouIiMhvokK7nkrKzCOnoFj9s2uYj6c7947vysYDGSzZ\ndtjZcURERKQeU6FdT8Wl5gDQSUuv17hL+7alSws//u+bXRRpXm0RERH5lVRo11OxZYV2WLAK7Zrm\n4e7GrAndiE3NYWF0orPjiIiISD2lQrueik/NwcfTjZb+Ps6O4pLGR7SkX0gznlu6m+MFxc6OIyIi\nIvWQCu16Ki41h9DAxri5GWdHcUnGGP48oRspWfm8uy7B2XFERESkHlKhXU/Fp+YQpv7ZDjWoYyBD\nOwXy6k/7yC0ocnYcERERqWdUaNdDRcUl7E/PJVSFtsP9cXxXUrMLeHftfmdHERERkXpGhXY9lHj0\nOEUlVi3atWBAaHOGdw5Sq7aIiIicMxXa9VBcWtmMIyq0a8U953UhLaeAd9aqr7aIiIhUnwrteii+\nbGq/0EAV2rWhf2hzRnQJ4rWfYtWqLSIiItWmQrseikvNwd/bgyA/L2dHaTBOtGq/vUat2iIiIlI9\nVRbaxphOxhjvsu9HG2NmGmOaOj6anElcag6hQY0xRlP71ZZ+IaWt2nOWx5KTr1ZtERERqVp1WrQ/\nAoqNMZ2BOUB74D2HppKzik/T1H7O8MfxXUlXq7aIiIhUU3UK7RJrbRFwCfCCtXYW0NqxseRM8ouK\nOXj0uKZTRzk0AAAgAElEQVT2c4K+HZoxqmswc5bvU6u2iIiIVKk6hXahMeZq4Ebgy7Jtno6LJGdz\nID2XEgthQb7OjtIg3XNeF47mFjJvTbyzo4iIiEgdV51CezowBHjCWhtnjAkD5js2lpxJXGouAGFB\nfk5O0jD16dCM0d2CmbM8lqy8QmfHERERkTqsykLbWrsduB+IKXscZ6192tHBpHJxqdkAhGlqP6e5\nd3xXMnILeWnZPmdHERERkTqsOrOOXAhsBL4pe9zbGPO5o4NJ5eJSc2nm60kTX/XecZaodk25rG87\n3lgZVz6nuYiIiMipqtN15C/AQCADwFq7EejowExyFvGpmnGkLrj//G54uhse/2qHs6OIiIhIHVWt\nwZDW2sxTtpU4IoxU7cQc2uJcLQJ8uHNsF5buOMLy3SnOjiMiIiJ1UHUK7W3GmGsAd2NMF2PMC8Bq\nB+eSShwvKObwsTw6qtCuE24aHkpIoC9/+3I7hcV67ykiIiInq06hfRfQA8indKGaTOBuR4aSysWn\nlfYHVot23eDt4c5DE7uzNzmbBT/vd3YcERERqWOqU2hPstY+ZK0dUPb1MHCRo4PJ6RLSSqf2C2mu\nQruuGB/RkkFhzfn30j0c03R/IiIiUkF1Cu0HqrlNHGx/emmLdodALVZTVxhjeGhSd9JzCnj1R033\nJyIiIv9zxkLbGHNBWX/stsaY/1T4eguo1vrTxpjzjTG7jDF7jTGzK9lvyq651xiz2RjTt8K+Pxpj\nthljthpjFhhjfH7F83Mp+9NzaerrSZNGmtqvLolq15SLe7fh9ZVxHMo47uw4IiIiUkecrUX7EBAN\n5AEbKnx9Dkyo6sLGGHfgJeACIAK42hgTccphFwBdyr5uBV4pO7ctMBPob63tCbgDV1X7WbmohLRc\nOjRXa3ZddN+Ebljg2W93OTuKiIiI1BEeZ9phrd0EbDLGvGet/TWdTwcCe621sQDGmPeBKcD2CsdM\nAd621lpgrTGmqTGmdYVsjYwxhYAvpYV/g3YgPZeebZs4O4ZUol0zX6YPC2XO8limDQ0lql1TZ0cS\nERERJ6tOH+1QY8wiY8x2Y0zsia9qnNcWOFDhcWLZtiqPsdYeBJ4F9gNJQKa19ttq3NNlFRWXkHj0\nuFq067A7xnQmyM+bBz7eQpGm+xMREWnwqlNov0lpl44iYAzwNvCOI0MZY5pR2todBrQBGhtjrjvD\nsbcaY6KNMdEpKa67cEhSZh5FJZYQDYSsswJ8PPnrRT3YdugYb62Od3YcERERcbLqFNqNrLXfA8Za\nm2Ct/QswqRrnHQTaV3jcrmxbdY45D4iz1qaUdVv5GBha2U2stXOstf2ttf2Dg4OrEat+2p9eOrVf\ne7Vo12kX9GzFed1b8M9vd3Og7GcmIiIiDVN1Cu18Y4wbsMcYc6cx5hLArxrnrQe6GGPCjDFelA5m\n/PyUYz4HbiibfWQwpV1EkijtMjLYGONrjDHAOGBHdZ+UKyqfQztQc2jXZcYY/jalJ24GHv50K6XD\nD0RERKQhqk6hfTelgxFnAv2A64AbqzrJWlsE3AksobRIXmit3WaMmWGMmVF22GIgFtgL/Be4vezc\ndcAiIAbYUpZzTvWfluvZn56Lp7uhVUCDn+WwzmvTtBH3TejGT7tT+OSXUz/EERERkYbijLOOQPkU\nfVdaa+8DsoHp53Jxa+1iSovpitterfC9Be44w7mPAY+dy/1c2f70HNo388XdzTg7ilTDDUNC+Wpz\nEo99vo3BHQNp07SRsyOJiIhILTtri7a1thgYXktZ5Cz2p+eqf3Y94u5m+OcVvSgusdz34SZKStSF\nREREpKGpTteRX4wxnxtjrjfGXHriy+HJpJy1loS0XM04Us+EBDbmkckRrN6XpllIREREGqCzdh0p\n4wOkAWMrbLOUzgQitSDzeCFZeUWaQ7seumpAe5ZuP8LT3+xkZNcgOrfwd3YkERERqSVVFtrW2nPq\nly0178SMIyq06x9jDE9eFsmE55Zz34eb+ei2oepnLyIi0kBUp+uIONmJObQ7qOtIvdTC34e/XNSD\njQcyeHNVnLPjiIiISC1RoV0PlBfaatGuty7q1YbzurfkmSW7iEvNcXYcERERqQVVFtrGmLDqbBPH\n2Z+WS5CfN75e1elSL3WRMYYnLumJl4cb9y/arFlIREREGoDqtGh/VMm2RTUdRM4sIT1HM464gJYB\nPjwyOYKf49OZtybe2XFERETEwc7YRGqMCQd6AE1Omc4vgNKZSKSW7E/LZVDHQGfHkBpweb92fLP1\nME9+vZNhnYPo2lKzkIiIiLiqs7VodwMmA02BCyt89QV+7/hoApBfVEzSsTz1z3YRxhieviwKf28P\nZi74hfyiYmdHEhEREQc5Y4u2tfYz4DNjzBBr7ZpazCQVJB49jrUaCOlKgv29eebyKG56K5pnvtnF\nw5MjnB1JREREHKA6o+v2GmMeBEIrHm+tvclRoeR/NLWfaxob3pIbhoQwd2UcI7sGM7JrsLMjiYiI\nSA2rTqH9GbACWAroc+5alnj0OADtm6nQdjUPTuzOuth07lrwC5/eMYywoMbOjiQiIi4oKfM481Yn\n8HFMIp2C/Zg2LJTzurfUAmq1oDqzjvhaa++31i601n504svhyQSAxKO5eLobWvh7OzuK1DAfT3f+\ne0N/3N0MN7+1nszcQmdHEhGROux4QTGLtyRRVFxSrePzCou5d+FGhj+9jDnL9xHZtgkJaTn8Yf4G\nRj2zjK+3JDk4sVSn0P7SGDPR4UmkUolHj9O2aSPc9K7TJXUI9OXV6/px4Ggud7wXU+0/niIi0vA8\n9vlWbn83hpeW7avyWGstsxZt5uOYg9w4JJSfZo3h9WkDWP7nMbx6XV+a+npy27sxPPfdbq3t4EDV\nKbTvprTYPm6MOWaMyTLGHHN0MCmVePQ47dRtxKUNDGvOExdHsnJvKv9YvNPZcUREpA76fscRFkYn\nEuzvzQs/7GHrwcyzHv/c0j18sekQfz6/G49eGEH7skkVPNzdOL9naz66bShT+7Xj+e/3cOeCGHIL\nimrjaTQ4VRba1lp/a62btbaRtTag7HFAbYQTOHg0l3bNGjk7hjjYFQPaM21oKG+simPJtsPOjiMi\nInXI0ZwCZn+8hfBW/nw1cziBfl7cu3AjeYWVD5379JeD/Of7PVzerx23jepU6THeHu48MzWKhyZ2\n55uth7nytbUkZ+U58mk0SFUOhjTGjKxsu7V2ec3HkYqOFxSTml2gQruBeGBiODH7jzLrw01EtA4o\nb30QEZGG7eHPtpKRW8C86QNp4e/DU5dFMf3N9fzru91MGxrK/vRc4lJz2HQgg5j9R9mTnM3gjs15\n4pJIjDlz11NjDL8f2ZGOwY25871fuPTl1bw1fSCdW/jV4rNzbdWZdWRWhe99gIHABmCsQxJJuYMZ\npVP7qetIw+Dt4c6LV/dl0gsruPO9GD6cMRRPd8PxwmIaebqf9Y+liIhU3w87j/DIp9s4XliMh5vB\n18udf1/Vh97tmzr83tZaCopL8PZwr9bxC9cf4KvNScya0I2INqUdCsZ0a8HVAzswZ3ksc5bHlh/b\npJEnfTo05cKoNtwwNBQvj+r0EIZx3Vvy/q2DuXneei57ZTVPXhrJ+IiWeLpX73w5syoLbWvthRUf\nG2PaA/92WCIpd6Bsaj+1aDccHQJ9eWZqL2a8s4G+f/+O44XFFJdYBoY2543pA/Dzrs57YxER1/dz\nXDor96bi7+1BE19P2jfzZXDH5mdtlLDW8saqeJ74ajtdW/ozNrwFhcUlfLHpEAvW7a+VQvvpb3bx\n6k/7CPDxoGWADxFtApg1oVuljWqr96Xy4CdbGN45iD+M7HjSvkcmd6dds0Y08/WiQ3NfOjT3pX3z\nRr+6UaZX+6Z8fNswpr/1M7e/G0OQnzeX9W3LdYND9Anrb/Br/tVOBLrXdBA5XWJ5oa1f8Ibk/J6t\nePbyXmw6kEFAIw9KLMxZHstNb61n3vSBNPKqXiuIiIgrWhebxr+X7mFNbNpp+0Z0CeLxi3sSEnj6\nugSFxSU89vk23lu3nwk9WvLclb3x9Sotg/IKi1my/TCPF/d0aCtuSYnl45hEerYNoF+HZhw+lsd3\n24/w7bYj/Ol3XZk2NBSPsvvvS8lmxvwNhAU15qVr+5ZvP8HXy4M7xnSu0XwdAn355p6R/LgrhYXR\nB3h9ZRzvrE3g6alRTI5qU6P3aiiq00f7BeDEvC9uQG8gxpGhpJTm0G64pvZrx9R+7cofR7QO4O73\nf+H3b0cz98b++Hiq2BaRhuVoTgEPfLyFb7YdJtjfm0cmR3D1wPYUlVgycwtZuuMI//x2N797bjkz\nx3Xh9yM6lnedyCss5s73Yli6I5nbRndi1u+6nTRt7gWRrfl04yHWxaYzvEuQw57DxsQMkrPyeXBi\ndy7u0xaAgxnHefTTrTz+1Q7eWZtAj7ZNCA305cvNSXi6u/HGtAE0aeTpsEyn8nR3Y3xES8ZHtORg\nxnFmLviFO9/7hfVx6cw6P5ysvELSsgtoGeBDsOqTKlWnRTu6wvdFwAJr7SoH5ZEKNIe2nHBhrzbk\nF5Vw34ebuO2dDbx2ff9q970TEanvVu5J5U8fbiQ9p4BZE7px8/CwkxocAnw8mT4sjAt6tuavX2zj\nmSW7+PSXg/zj0kjCW/nz+7ejWRubzt+m9OCGIaGnXX9U12B8vdxZvDXJoYX2km2H8XAzjAlvUb6t\nbdNGzL2xP19vPcwH6w+w7WAm32w9TCNPd+bdNNCp3TbaNm3E+7cO5umvdzJ3ZRzz1iSU7/Pz9uDl\na/sysmuw0/LVB8baqicpN8Z4AV3LHu6y1tbJJez69+9vo6Ojqz6wnpjy0ir8vT1455ZBzo4idcR7\n6/bz4CdbmNCjJS9e01cDVURcUH5RMXe+9wsjugRVWhQ2NK/8uI+nv9lJp+DGPH9VH3q2bVLlOd/v\nOMKjn23jYMZxWjfxITkrn39e3qu8Fbkyd74Xw5p9afz80HkOWZrcWsuYZ3+kfXNf5t989n/XC4tL\nKC6xderTy+W7U9h0IINAP2+a+nryn+/3sDc5m39cGskV/ds7O16tM8ZssNb2r+q4Kv+VNsaMBvYA\nLwEvA7vPNOWf1CzNoS2numZQBx67MIIl245w78JNFGs1LxGX8+ySXXy3/QhPfLWDxKO5zo7jVO//\nvJ+nv9nJhb3a8OVdI6pVZEPpLBrf3TuSW0d2pLDYMuf6fmctsgEmRrYmLaeAdXGn9/2G0jdAhb9h\n9d7dR7KJT8vl/J6tqjzW092tThXZACO7BnPXuC5cM6gDEyNb8+GMIQzpFMifF23mn9/u0uqSZ1Cd\n5rB/Ar+z1o6y1o4EJgDPOTaWaA5tOZPpw8KYfUE4X2w6xB8/OPOCBSJS/6zYk8J/V8QxKbI1xpTO\nUNFQLd1+hAc/2cKorsH864pe5zwQ3NfLgwcndif64fMY171llceP6daCRp7ufL3l5EXDrLUsjD5A\n779+x/Cnf+DFH/aQlp1/TlkAvtl6GGNgfETVWeoDfx9P3pg2gCv6t+OFH/Yy/a31pOcUVOvc3/KG\npb6pTqHtaa0t/z/dWrsbqL1e+Q2U5tCWs5kxqhN/Pr8bn286xKUvryYhLcfZkUTkN0rLzufehZvo\n0sKPZy/vxa0jO/HFpkNsSEh3drRa98v+o9y5IIaebZvw8rW1002ukZc7Y8KD+Wbb4fJPC7Pzi7jn\ng438edFmoto1oWtLf579djdDnvqBv32xnez8ypctzyssZvGWJFbsSSnftmTbYfp2aEYLfx+HP5fa\n4unuxtOXRfH4xT1Zsy+Nyf9ZQcz+o2c957/LY+n+yDfcMi+apduPUOTiRXe1BkMaY+YC75Q9vpaT\nB0iKA2gObanK7aM7071VAPd8sJHJL6zkX1f0dpmWEpGGxlrL/R9tJjO3kLdvKp3Gc8aojnywfj9/\n+3IHn9w29KwD499eE8/e5GxmTeiGv8/Z28JKSizvrz/AN9sOExboS3jrAMKCGmOAYlvaL7hP+6ZO\nXSTr719uJ7CxN29MG0DjWlw/4IKerVm85TA3vbWe7Pwi4lJzyMgt4N7xXbljTGfc3Qx7k7P47/I4\n3lwdxzdbk/jblJ6M6hZM4tHjxKVms3RHMl9sOkRWXmkRPjmqNb8f0ZHtScd4cGJ4rT2X2mKM4brB\nIfRq15Tb39vAla+t4cGJ3Zk2NPS036EvNh3iicU76NOhKZsSM1i64witAnx4emoUo1x0UGWVgyGN\nMd7AHcDwsk0rgJettef+uYmDudJgyPlrE3jk062se3AcLQNc592v1LwD6bnc/m4MWw9l8tDE7tw8\nPEyrSIrUM/PXxPPIZ9t4dHIENw0PK9/+0YZE/vThJh64IJybh4edNpcywKINidz34SYAQgJ9eeHq\nPkS1q3zhlR1Jx3jwky38sj+DkEBfUrPyySk4vftZnw5NeWRyBH07NKuZJ3gO4lJzGPPsjzxwQTh/\nGNWpVu+dk1/EpS+vpqC4hFYBPrRu4sPVgzowILT5acduSDjKgx9vYdeRLNzdTHkruI+nGxf0bM1l\nfdux8cBRnv9+DyUWikssP80aXekc364iM7eQP324iaU7jjApsjVPXRZZ/sbv57h0rpu7jl7tmzD/\n5kF4uBmW7Urhn9/uYteRLGaO7cLMcV0cMhDVEao7GLI6hXZjIM9aW1z22B3wttbWuREarlRoP/n1\nDt5cGc/Ov5+v6f2kSnmFxfzxg418vfUw04aG8sjkiHrzx6ou2JKYSZeWfnVu8JE0DLuPZHHhCysZ\n3DGQt6YPOOmNckmJ5dq561gTm0abJj7cODSUK/q3p1ljLwBW7U3lxjd+ZlDH5twxpjOzPtxMclYe\nt43uzKV92hIaVFrUbTuUyZur4vnkl4M0beTJQ5O6c0mftlhbOo9zQlouxoC7m2FfSjb/XrqHlKx8\nLuzVhjvGdCK8VcBJmWJTcwgN9K208P+t/vXdbl74YQ9rZo+jVZO63dBUUFTCe+sSSMnOJzSwMWFB\njQlvHXDSKr47ko5x/0eb8fVy5/1bhzgxbe2w1jJneSz/t2QXTRp50jGoMS2b+LBqbyrNfb346Lah\n5b+/UDom7eFPt/JRTCLDOwfx/FW9CfSr+/Nz12ShvRY4z1qbXfbYD/jWWju0RpLWIFcqtO94L4bt\nh46x7L7Rzo4i9URJieUfi3cwd2UcY8Nb8OSlkfo0pBrWxaZx5Zy1XNCzFS9f21efBkityiss5uKX\nVpGanc/Xd4+sdAGQkhLLDzuTeX1lHGti03B3M/QLacawTkHMXRFLm6aN+PC2IQT4eJKZW8gDn2xm\ncdmAvvBW/vh5exCdcJRGnu5cOaA995zXhaa+Xqfdp6Kc/CJe+2kfc1bEkldYwtBOgVzWtx1bDmay\neEsSyVn5dAxuzP3nh/O7iJY19v+NtZaRzywjNLBxlVPg1TfW2gb192V9fDrvrk0gKTOPI8fy8PF0\nZ871/ekQePrYsxMDTh/5bBvNfb148Zo+9K/kU4Ts/CKeXLyD/em55OQXkVtQzPVDQrh2UEhtPKWT\n1GShvdFa27uqbXWBKxXaU15aRYCPh8v9oRHHm78mnr9/tQNvdzfum9CN6waHqHX7DKy1XPrKarYe\nzKSw2PLEJT2d8gdbGiZrLY99vo231yTw5rQBJy1icibbDx3j661JLN2RzI6kY7QM8OaT24fRpunJ\n43kOpOfy7fYjLNl6mPTcAq7s354rBrQ/5xUGM3ILWPDzAd5eE09SZh5eHm6M6RbMwLBA3luXwL6U\nHPqHNOOpyyLp3ML/nK5dmfXx6Vz+6hr+dUUvLu3bruoTxKVsO5TJ7e/GcPDoce4/v7S71IlP9bPz\ni5j2xs/8ciCDyLZN8PP2wNfLnSm92zIpqnWtZ63JQnsVcJe1NqbscT/gRWttnfv8w5UK7f6Pf8f4\niJY8eWmUs6NIPRSfmsMjn21lxZ5UIts24cGJ3RnSKdDZseqcJdsO84f5G/jHJZF8vTWJn+PS+fzO\n4XRrdXLBUFJiWbrjCB2DG9dIMSGSeDSXBz/ZyvLdKUwfFspjF/Y452skZR7Hy92tVj5mLywuYXNi\nJt3KWsgBiopLWBidyD+/3UVBUQn/uaYPY7pV/WbhbB74eDOfbTzE+ofOq9VBkFJ3HMsrZNaHm1iy\n7QgRrQOYNaEbA8KaM/3Nn4nZn8ELV/dhYmTtF9anqslCewDwPnAIMEAr4Epr7YaaCFqTXKXQPl5Q\nTPdHv2HWhG7cMaazs+NIPWWt5fNNh3jq650kZeYxulswsy8IP6mv5QlFxSUO6WtZ06y15BeV1Ehf\n6qLiEib8ezkAS+4ZSXpuAROfX0Hzxl68f+sQmpf1IVwbm8YTX+1gy8FMvDzceGhid24YElL+EfDR\nnAK8Pd3w9VJRIFUrKbHMX5vA09/sBOD+88O5fnBIvR6LczDjOLfMi2bX4WM8NCmCm4adPttEdeQV\nFjPgiaWMj2jJv66ocx+aSy2y1vLpxoM8990e9qfn0qSRJ9n5Rfznqj5Oab2uTI0V2mUX8wS6lT2s\n9hLsxpjzgecBd2CutfapU/absv0TgVxgWoWW86bAXKAnYIGbrLVrznY/Vym09yZncd6/lvP8Vb2Z\n0vvsK1mJVCWvsJh5q+N5adlecgqKmTWhG7eO6IibmyGvsJgnF+9gwc8HmDYslJnjupw0iKeueWnZ\nXp77bjeju7Xg0r5tGRve4qxFd15hMevi0ikqLsHNGLw83OjSwo8WAT58sH4/93+0hVev68v5PUv/\ncC/fncINb/wMQFNfT4L9vNmTnE2bJj7MHNeFJdsOs2xXCuPCWzAwrDnfbT/Chv1H6dbSn49vH6pi\nW84q+Vgef/pwEyv2pDKyazD/uKSny6yVkJNfxL0LN7Jk2xEmRrbirxf1rLS/eUXWlo4riUvNZVBY\ncwpLSvi/b3bxzs2DGN4lqJaSS11WUFTCwugDvLtuP3eM6cTkqDbOjlSuJlu0PYHbgBPLrv8IvFZV\nsV02O8luYDyQCKwHrrbWbq9wzETgLkoL7UHA89baQWX75gErrLVzjTFegK+1NuNs93SVQnvZrmSm\nv7mej24bQr+Q0wcDiPwaGbkFPPTJVr7aksSYbsH8YVQnHvtsG7uOZDEorDnr4tJpGeDNQ5MiuDCq\ndZ0btJOdX8TQJ78n2N+brLwikrPyCfDxYFJUGy7r25Z+Ic3KM2fmFvLOugTeXBVPaiUruLUK8CG3\noIiOwX58cvvQk57rhoR0YhIyiE/LIfHocQaGNefm4WH4eLpjreXNVfE89fVOCopL6NEmgP4hzXh7\nbQIXRrXh+at617nX7XhBMT6ebnUuV0Pz/Y4jzFq0mdyCIh6eFMG1gzq43M+kpMTyyk/7eH7pHhp5\nufPo5Agu7dv2jM/ztZ/28eTXO2kZ4M2RY6X/n7YK8GHV7LEaVyJ1XnUL7eo0v7xC6UqQL5c9vr5s\n2y1VnDcQ2GutjS0L9D4wBdhe4ZgpwNu2tNpfa4xpaoxpTWnr9khgGoC1tgCo3rqeLiCxfLEa12jp\nkLqhadlI7sFrm/P3L3ewbFcKQX7evDV9AKO7tSBm/1Ee/WwrMxf8wscxiTxxSSRtm56+YNLe5Cxu\neyeGG4aEcP2Q0FrL/+7aBI7lFTH/5kH0bNuEVXtT+TgmkU9/OciCn/cT5OeFt4c7RSUlHM0tpKCo\nhJFdg5k2NIQgP2+KSyzHC4vZkZTFlsQM9qXk8NiFEacVAf1Cmp/xDa4xhpuGhzExsjXF1pa/Pi0C\nfHhmyS56tW/KzRXmQHamwuIS/rsiln8v3cPILsH8+6redfrTirruvXX7eXNVXPn0el4e1e9qtWxn\nMjfPi6Z76wBeuLq3y/bzd3Mz3DGmMxN6tOT+j7bwpw838Z8f9jCmWwvGhrdgcMfA8tdt2c5knvpm\nJ5MiW/PiNX1IzspnXVw6Ic19VWSLS6lOi/Yma22vqrZVct5U4Hxr7S1lj68HBllr76xwzJfAU9ba\nlWWPvwfuB4qAOZQW5b2ADcDd1tqzrjPtKi3amkNbHG3rwUy+2HyI34/oSFCFgVTFJZZ5q+N5Zsku\n3AzMntidawd2KP89PJyZx6UvryLpWB7Wwv9dFsUVA9o7PG9eYTHDn15G99b+p83Ek51fxDdbD7Nm\nXxrGgIebwc/bg8v6taN769P7ozuCtZY/zN/A9zuTeWPagFpf4Swnv4jHv9qOv48n3Vr606yxJ88s\n2c2OpGMMDGvOhoSjdApuzH9v6O/Si2U4yolFYwIbe5GWU0D75o3443ldubh325P+Rh/NKeDb7YeZ\nFNWm/E3N4cw8Lnh+Oa2aNOKT24c2mLnaS0osH8UksnhLEqv3pZFfVEKTRp5c0LMVQzsH8dDHW2jf\n3JdFtw1Rlyupl2qy60gMcLm1dl/Z447AImtt3yrO+y2FNsBaYJi1dp0x5nngmLX2kUrucytwK0CH\nDh36JSQkVPWc6zzNoS3OdiA9lwc+3sLKvalEtA7ggYnhRLVtyhWvreFgxnHm3zyQ55buYcWeFP59\npePHEpxYKXXB7wfX2dlTsvIKmfLSKmJTcohs24Sp/drRu31TUrLyOZKVR2hgY4Z1dky/08e/3M7c\nlXF4ubtRUFwCQMsAb/42pScTerRi1d5U7ngvBoDzurekqLiEohLL5KjW5f3TpXJfb0nijvdiGNIp\nkNdvHMCa2DT++e0uth48Rs+2ATwyKYIBoc1ZFJPIU1/vJD2ngC4t/Hj1+n6EBjbmmv+uZcvBTL64\nazidgv2c/XSc4nhBMav2pvLVliSWbDtMbkExgY29+Pyu4ZV+aiZSH9RkoT0OeBOIpXTWkRBKByb+\nUMV5Q4C/WGsnlD1+AMBa+2SFY14DfrTWLih7vAsYTengx7XW2tCy7SOA2dbaSWe7p6u0aGsObakL\nTsxa8uy3uziQfpzmjb3Iyitk3vSBDO0cxPGC4v9v777D46ru/I+/v+pdsmxLFnJv2MbGYAtjMNWU\nYDj0atMAACAASURBVJqTAKEFCAkhJEBIX5JfKvvbTbLZBMJCILQEAgmhJMHxQug4gGOwjcHGXa6S\nXGSrWr2d/eOMkNyFrasrjz6v55lHM3PvzHynaO7nnjn3HK773bss3lTBeccMYmj/FIZmp3DK6AEM\nye6+bk/NrW2c+d9vkJOeyLNfPrlX92utqmvmL0uKeWZxMcu3VO+1/O4rj+fiyQc/mGfjzlrmrytj\n1sRBu82gti/LiquYfe9bXDltKD+5+Bg2ltVRVF7H1OH9yEjqGDN5U1kt33r6A7ZUNhAXa9Q3tVK6\nq5FvnTuWm88cjZlRVtPIb/+5nvqmVs6flMe0Ednd/jP+nA+28O9zV3DDKSN2m1J8V0Mz89bsIDsl\ngZED08jNSDzk97qhuZW4GDvskXReX13KjY8t4tjBWTz2+WkfDTfX/r/x8xdWsaWqgSHZyRSV1zN1\nWD8uP2EIP3thFc2Rrkv/u2wrv7xsMpdM1ZjQ4EP3vDWljBqYxpjc6OxCI31Ddwbt9t+VPxp1BMA5\nt/cRRrvfLg5/MORZQAn+YMirnHPLO61zAXALHQdD3u2cmxZZ9iZwg3NutZn9GEh1zn37QI8ZLUFb\nY2hLb9LY0srjCzbz+IJNfPPcsbsd9V3T2MIP/vYh7xdVUlxRR3OrwwzOGDuQz04fxhlH5xx2UHtm\ncTHfevoDHr6ugLPG5x7u0+kxK7dWU1ReR25GEtmpCXzz6Q9Ysrniox2VfWlsaeW389Zzz+uFNLW0\nkZ4Yx42njeTzp4zY55jCLa1tzL73bXbsauTlb5z+sSYjaWxp5d+eWcrf3t/CZwoGM3xAKr95fR31\nza0kxMZQ39xKTnoik4dkEWNgGEdlJXPR5DyOG5J1SCF4/jo/XXhGUjxltU0cc1QG3zr3aN5cu5On\nFhVR09jy0bppiXF8pmAIN5856qNxolta23i/qBIzGJiWxID0BJpbHVV1zZTVNrJ4UwXz1uzgnfXl\nJMbFcMqYAZxx9EDOHp/7sceaXrC+jOseeZfROWn88YvT9/naNjS38uA/1/PKyu1cfeIwLp06mJgY\no6Syni8/vpilxVV8ekq+hqoTiULd2nVkz24i+7puP7c9H7gLP7zfI865/zCzmwCcc/dHhve7BzgP\nfwDk9c65RZHbHocf3i8B35p+vXOu4kCPFw1BW2Noy5Gqtc2xqayWvy0p4cmFRZTuaiQ/K5krp/kZ\n6XLSk2hubWPHrkZSE+LITDl4KKyqa+asX80jLzOJObfM6NWt2QdTVd/MZ+73XW8e/fw0cjMSqW1s\npby2iaKKOjaX1fHCh1tZt6OWC47N45rpw3jkrQ28tGI7aYlxpCbG4hzEx8Zw2tiBXDQ5j6XFVfzs\nhVXcd/UUZh3CBA7OOe58eQ13v1YIwNnjc7h91niOykri1ZWlzF26hU1ldTgHbc6xqbyOppY2hvVP\n4aSR/SPvh8M5/AlHcnwsY3LTGZ+XztjcdNIjreqrt+3i0vvnMygjiWduOpn563bywznL2bGrkbgY\n48Jj87hy2lBa2hzrd9by3qYKnnu/hOT4WK49eTg7dzXy8srtVNYdeHTZ0TlpnD52ILWNLbyxegfb\nqhtIio/hymlDufG0keRlHryrwvtFlVz94ALyspL5843TD2lCmIbmVl5cvo1zJwwiOaFv9MsW6UsO\nO2ib2SAgH3gcuArfbQQgA7jfOTeum2rtNtEQtDWGtkSD5tY2Xlmxncff2cTbhWXExRhZKfHsrOkY\nPGjkgFQmD8nirPE5nD8xb58H/t7+7FKeXlzMczfPYGJ+Zk8+hUBsrarnkt/MZ0tVw17LYmOMMTlp\n3D5rHGd0ml3vvc0VPLO4mNZWR0yMD+xvrN5BXVMr4PtcP3jt1MPaCXlt1XZSE+I4ceSB+79XNzTz\n4ofbeO79LazaVg0YZn7jYJFW75rGlt1apgf3S2bcoHQ+LKmmzTn+evOMj/rlVtU38/qqUk4a1Z/c\njKS9Hq+wtIZfvbya55dtIz0xjrPG53DuMYNISYhlx65GdtQ0khAbQ1ZKAlnJ8YzLS99ttCbnHCu2\nVvO7tzfy1yUlxJoxY3R/8vslk5eZzAnDs5k2YvcRZj4sqeLqh94hIzmOp790MoMy965LRKQ7gvZ1\n+OH1CoDO6XUX8Hvn3F+6oc5uFQ1BW2NoS7RZv6OGpxYVU93QTG56EjkZiZTXNvFBUSVLiirZsauR\nifkZ3H7e+N0mqViwvowrHljAjaeN5Hvnjw/xGXSvksp6Xl6+jeSEWFIT48hKTmBodgp5WUnEd7FP\ncX1TK6+u2s4768u5ZebofYbUsDjnKKmsZ/W2XayKnFZvq6appY17r57CMUd9/B2m7dUN9EtJ+FhD\n6u2pqLyOB99cz8KNFWytqv+oZfya6cP43vnjSYqP4ffzN/LTF1bRPzWBp750UrceayAi0aU7u45c\n4px7ttsqC1A0BO320RXe+d5ZvWrjKRKEtjZ/UNkvXlxNSWU9BcP6cf6kPM44eiA3PLqIljbHi187\nTT+9S7eraWzhzpfX8PBbGxg5MJUh/VKYt8bP+vlflx57SN1FRKTv6M4Ja141s1/RMTPkPOAO51zV\n4RQo+1ZcUUdCbAwD9SUvfUBMjPHJ4/OZNWkQf/jXJp5aVMQdc1dwx1y//A9fmKaQLYFIS4zjBxdO\nYOa4HL719Af8a30ZP7n4GK49adgRfSyAiPQuXQnaDwMfAp+JXL4GP9zfp4Mqqi8rrqgnv1+yJqqR\nPiUxLpYbTh3JDaeOZMPOWl5cvo2UhFhOHdOzE79I3zNj9ABe/sbp7Gpo7tKBkiIiH0dXgvYo59wl\nnS7/xMzeD6qgvq64op7B/fRlL33XiAGp3HT6qLDLkD4kLTFO09OLSCC6cmRJvZmd0n7BzGYA9cGV\n1LeVVNQpaIuIiIhEga7swt8EPGZmmfhRnMrxo5FIN6tvamVnTdNuw1OJiIiIyJHpoEHbOfcBMNnM\nMiKX955XWLpFSWUdgFq0RURERKLAQYN2ZAr2S4DhQFz70djOuTsCrawPKqrwPXIUtEVERESOfF3p\nOvIcUAUsBhqDLadvK/4oaKvriIiIiMiRritBe7Bz7rzAKxGNoS0iIiISRboy6sh8M5sUeCWiMbRF\nREREosh+W7TNbBngIutcb2br8V1HDHDOuWN7psS+Q2Noi4iIiESPA3UdubDHqhDAj6E9YUJu2GWI\niIiISDc4UNDejh9DezSwDHjYOdfSI1X1QRpDW0RERCS6HKiP9qNAAT5kzwJ+2SMV9VEaQ1tEREQk\nuhyoRXuCc24SgJk9DLzbMyX1TRpDW0RERCS6HKhFu7n9jLqMBK99DO38LHUdEREREYkGB2rRnmxm\n7dOtG5Acudw+6khG4NX1Ie1jaOekawxtERERkWiw36DtnIvtyUL6uuKKeo7KStIY2iIiIiJRoisT\n1kgPKKmo14gjIiIiIlFEQbuX0GQ1IiIiItFFQbsXaGhuZWdNo4K2iIiISBRR0O4Fij8a2k9dR0RE\nRESihYJ2L1BcoclqRERERKKNgnYvoBZtERERkeijoN0LFFfUEx9rGkNbREREJIooaPcCJZX1HJWV\nrDG0RURERKKIgnYvUFxRp/7ZIiIiIlFGQbsXKK6oZ3CW+meLiIiIRBMF7ZA1NLeyY5fG0BYRERGJ\nNgraISupjIw4kq2gLSIiIhJNFLRDpqH9RERERKKTgnbI2ieryc9Si7aIiIhINFHQDllxRT1xMUZu\nRlLYpYiIiIhINwo0aJvZeWa22swKzez2fSw3M7s7snypmU3ZY3msmS0xs7lB1hmmkgo/hnasxtAW\nERERiSqBBW0ziwXuBWYBE4ArzWzCHqvNAsZETjcC9+2x/DZgZVA19gYaQ1tEREQkOgXZoj0NKHTO\nrXfONQFPArP3WGc28JjzFgBZZpYHYGaDgQuAhwKsMXTFFfUK2iIiIiJRKMignQ8UdbpcHLmuq+vc\nBXwHaDvQg5jZjWa2yMwW7dix4/Aq7mENza2U7mrUiCMiIiIiUahXHgxpZhcCpc65xQdb1zn3gHOu\nwDlXMHDgwB6orvtsaR9DWy3aIiIiIlEnyKBdAgzpdHlw5LqurDMDuNjMNuK7nMw0s8eDKzUcGkNb\nREREJHoFGbQXAmPMbISZJQBXAHP2WGcOcG1k9JHpQJVzbqtz7rvOucHOueGR273mnPtsgLWGoj1o\n56tFW0RERCTqxAV1x865FjO7BXgRiAUecc4tN7ObIsvvB54HzgcKgTrg+qDq6Y1KKuv8GNrpiWGX\nIiIiIiLdLLCgDeCcex4fpjtfd3+n8w64+SD38QbwRgDlha64op68rCTiYntlV3kREREROQxKeCEq\nrqhncJb6Z4uIiIhEIwXtEGmyGhEREZHopaAdksaWVrZXawxtERERkWiloB2SLZUNgMbQFhEREYlW\nCtohKa6oAzS0n4iIiEi0UtAOScdkNQraIiIiItFIQTskJRX1xMYYgzKSwi5FRERERAKgoB2S4oo6\n8jI1hraIiIhItFLKC0lxRb26jYiIiIhEMQXtkPigraH9RERERKKVgnYIGlta2b6rQS3aIiIiIlFM\nQTsEWysbcA7ysxS0RURERKKVgnYIOob2U9cRERERkWiloB2Ckko/WY26joiIiIhELwXtEBRHxtDO\ny9QY2iIiIiLRSkE7BMUV9QzK0BjaIiIiItFMSS8ExRV16jYiIiIiEuUUtEOgMbRFREREop+Cdg9r\namljW3UD+WrRFhER6R5NdVC+HlpbPv5tnYO6cmis8edFulFc2AX0NVur6nFOI46IiIgctuLFMPc2\n2L4CXCvEJcHAcXDhnZA/BXauhfVv+CBdWwq1O6GtBS74JaQPggX3wat3QLMfDQyLhcR0uPU9SO0P\nSx6HTfNh4NH+fvsNh8zBkJC6dy3OQc12aKiGgWP9dduXQ2wCZOT7+kpXQlsrDDvJL1/6FLQ0Qnyy\nrz0hBTKHwoDRPfHqSQ9Q0O5hHWNoK2iLBKq+wm8gty2DpCyY+jmIT4LaMmhrhsQMv3EzO/h9Odex\n3jsPwM7VUFXccRp5OnzmMb/8gTPAtUFaLiRlQky836hOudYvX/5X3+q2Y6Xf6FYVwZTrYNoX/Qb6\nyasgOcvXnJwFyf1g+KkwZFpHa1tXapae1dII6+dB8buQX+A/E/FR+j3fuAtqd0D2SH95y/v+M5+Y\n7gNoQpr/GxN78PtyDrZ+AI3V/vPfuMufz5kAI06FplqYc6t/fc3AYqB6i/+fmXINpPSDlP5wytch\nayjsXOPDbXKWv/8N/4Tnv+XPJ/eDlAE++LY0+usGjvPfDZmDobW5o47ENL+8oRoKX4H3n+ioOSYO\nvl/qn99bd0LxIv99U7rC/+0/Bm5d5Ned+w0oWrD7cx56Enz+H/78Kz+B6uLdl4+/GC7/gz//6MW+\ntvwpkD8Vco6BuISDv669iXP+fdnyPqQOgNFn+euXPePfu+wRkJztPz9R+N2moN3DSiJBe4j6aIt8\nPO1hd+PbsOgRv8G1mI6N75nf8xukFc/BvP/yG1vafwY2OOEL/uwb/wkLH/LnY+J8GE7Ohq/8C2Lj\noWghbF8GdWU+lG/9wLcyffZZf5v3HvXhOnOIPw2dDoOnddR51PFQVeJbtsoKfahOyfbLWlvg6et9\nXTFx0H+0byFLyvTLWxp9a9vOQr/BbqiElgaY9QsftMvWwW9Pg/6jIrcd5oNAwedh0ERY9zr89Uv+\nJ/CkDB/20wfB2T+GnPH+uS18CFqb/Cm5n38OU66FjDx/+w+f8bWk5vjbJ6TCpEt9aCxdCRWbIC3H\nbxRbGnzN+VP3vYF0zrcUtrf+VRb54NRcB8ULYdPb0NIEVz3ply96xP/NHOrfy6RM/zjtoaeztlZo\nqPLLY+M/zifp0Lz+U/jgjxCf0hEmk/vBZx71y39/gX9O7eKS/es+/SZoroe7j/fvRf/R0G8ENNXA\nqJkw5hwf8D74Ewyb4V/nujLfAjvwaH+b5gb/eWh/3JhD7PVZu9N/Jis3Q+UmOOUbPixWbob6St/i\n2trs1ytfByff6m/39PVQsth/jqu3+M/2UVPgxtf98jm3+B3azoadAtf/rz8/9+uR/98Y/5mur/Rh\n8/Rv++UPzvSP3dm0G33Qjk2ALUv864mLBPqMjs9U9ki49rn9P+fjrobxF/kd130F1FFn+tP+nPQV\nf6qvgB1rIq9VRcdOREO1/79MyoAJs/0OQv9OrdGf+E//WlZu9vXnTIDciR3LvzTP/z80N0BLvf//\nSIh83psb/OdhzYsdQT82Ec75CUz/sv9cffis31kYOG7v/5OydbD4d7D8b/5+Y+L86/mFlyAz3++E\nbHjTXxcb7/+Xkvv51ys2HmpK/e3a/7+c85/bnPH+8ss/hMLX/OehpdH/v+aMh6uf8svn3wNb3oON\nb/l1AI4+3wfttjZ47hb/nNtZrG9wmPVz/1iPX+Jf16SsyPd0lv/cDJ3u61r6FKQOhPEX7v/96wUU\ntHtYUUUdMQaDNIa2iN9Q1JVDa6P/Mk3J9hv69g3/9g/9aePbcN5PfSthfbnf8Lq2yMkBDmbc5u/T\nYn04GX8xDD/Fh8Dmuo6NxbGX+41dQ1VH61VMbMfyeT+Hwpf9+fhUyBkHg6Z31PzF1yAucf/P6cI7\n97/MYuCWRT7k9h+19/2kDexo6froNWroCLFxCT4UlxX6DdiK5/yGaPTZPminD4Kx5/kNdUMV1Gzz\nob+lwd9+1xbYPN9vWGPifGCo2Q7HfArIg4oNsPYVX2dtqQ/9AGM/4Tf4S//sW/D29IOd/vV781ew\n+gX/2PXl/v6T+8G3C/16/7gdVs3tuF3mEB+E2r39a6jYuPt9j53VEcTvmeY39E21/jFwMPFSuPRh\nvxPz99s6bpeQ4lsvh8/wn4PWZtj4pl+vdoc/xcb7sNseHNq1NMGaF/yG/JO/8Rv5zHwYfIJ/TRoj\nNdTu6LjNKV/3r+vQ6VD0Lqz5B9Tt9Mvik2HMuf7Xi6J3fEtefIp/v8ac43fm5ty69+t62e/9e7N5\nPvzhUx3XJ2X5naCLfu1/LVn7CrzyI98a3FTj6zODzz0Pg6fCe3/woai+fPf7n3yVf17v/Bb+dc/u\ny+JT/Q5cQqoPNzFx/vnkHAP9R0Le5I51L7zLvxZNtZEaav1za7dzbWTHF/95SM7qCNZmcOWf/GuU\nmOHDXlKmPw/+Pfrqkr1fm66KT/Knw5XcD4ae6E+dnf0jf9qfwVP9aX9SB+x/WXwSXPVn/x1XVeR3\ndkoW+51zgPIN8NzNHetnDYWEdF/P2E/478kF9/nPXka+/yWvtbnjMde+BPP/Z+/H/UGZ//vGTzt2\nftslpMHtRX5nLzaho7U9Lsl/nyZlday7YZ5/30ec5k9DTuz4XJjBrYv9TkjFpo6Ghbzj/PKWBv8/\nXrkp8n1S6es/7dv+f6y5HuZ+zf/a18uDtrko6vhfUFDgFi1aFHYZB/TVPy1hSVEFb35nZtiliATP\nOd8qUrHRB8Kc8f4L9Y9X+OtqtnWsO/P7/ku0cjPcNanj+oQ0/+V7+nd80A5aVYkPminZBw7UvUHn\nLi2HqqXRd2/Zs5W0rc0Hs+Z6yDjK74zU7vQbxZrtPkzFJfqANGqmX/72r/3P7O07TcnZvsXppK/4\n+9y8AKpLfGg7agpkDdnjMVth1zYfKqqK/YY7I9+HBoAX/5///MSn+J+ck7P8Z2rkGT5A/7o9/EVa\n3hqq/Gdq5vd97b8YtffzP+9nvnWwdBU8cZkPeru2+ueengeXP3HgoHQoWlv869X+3jnnw+jm+X7n\nMaW/f+1yj/GvY1WxD0WNNf55tfc3PuO7/vlvXgBv3+1bNBPS/F/XBifdCum5vuVy2TO+hbz/mE79\njCO/rO5Y47syxcT7nbmEdB+kuyOgSnBaW3wQLV3pu63sXOP/X6d/xe9gNjf4HZ+0gfu/D+f8zmNL\no1+3obJjx7N4MexYtfuvDcnZfmc+tofbadt/HQO/89fW5rcfCWl+2xICM1vsnCs46HoK2j1r9j1v\nkZ4Uz+M3nHjwlUWOJO2hzzn482f9z5aVmzq+HKdeDxfd5cPUY7Mha5jf4KcN9K0hg46F3Ak+TKya\n679Mc4+BrOGH/lO59G2tzT5ExCf7VuqSRT5Mpg7wOwAtDT7wJvfzrYPzfu7DRkIaTLrMdynoSj9j\nEelzuhq01XWkh20qr+OCSXlhlyHSPZrq/MFfHzzpu0Nc+ogP2y0NvmvE6LN8mO433PchBB9cPjd3\n//eZmAaTr+iJ6iXaxcZ3dAmKS4BhJ+++vHOf1uwR8Kn7e642EekTFLR7UFVdM5V1zQzrrwMh5Qg3\n/x5Y9hRs+9D/rJiQDpMv72jVbj9wUEREpA9T0O5Bm8prARjWfx/jb4r0VrU7/QGJJYvhxJt8l464\nRN8P95Svw+ACf6DLvsaVFRER6cMUtHvQxjLfV1Ut2tLrtLZA2VrfbzV1AJS8By98xx/4VlvasV7u\nRH9g2rQv+pOIiIjsl4J2D9pc5lu0h2YraEsv4JyfaGHRI7Dib/6gxQvvgoLr/agO8Skw9lwYON4P\nHZc78cBDUYmIiMhuFLR70KayOnLSE0lJ0MsuIWjc5ScMyZ3gLz98jp9gIyHdj7Aw7GQ/3jD4saOv\nmxNerSIiIlFAia8HbSqrU7cR6VltbbD+dT8qyMq/+4kFrn/eL+s/xk8WMulSP3awiIiIdCsF7R60\nqbyWU8ccYOB4ke609mV46ft+woGkTDjuSji207B5n7ovvNpERET6AAXtHlLf1Mr26kaGqX+2BKmq\nODJrXrafYQ/g0w/B+Is0y5uIiEgPC3S6NTM7z8xWm1mhmd2+j+VmZndHli81symR64eY2etmtsLM\nlpvZbUHW2RM2l0dGHBmgIdCkGzkHW5bA6z+F+0+FO4+BpU/5ZcddBV+eD8deppAtIiISgsBatM0s\nFrgXOAcoBhaa2Rzn3IpOq80CxkROJwL3Rf62AN90zr1nZunAYjN7eY/bHlE2RkYcUYu2HLb2SWGa\nG+B/pkJ1MWAwZBqc/RM4+jy/nqaOFhERCVWQXUemAYXOufUAZvYkMBvoHJZnA4855xywwMyyzCzP\nObcV2ArgnNtlZiuB/D1ue0TZHBlDe7gmq5FDUbYO1r0Gha+AxcKVf/St1JMvh/6jYcy5GnpPRESk\nlwkyaOcDRZ0uF+Nbqw+2Tj6RkA1gZsOB44F3giiyp2wsqyUzOZ7MlPiwS5EjyaLfwb/ugbJCf7nf\ncDj6go5W7bN+GGp5IiIisn+9+mBIM0sDngW+5pyr3s86NwI3AgwdOrQHq/t4NpfXMVxD+8nBNO7y\nQ/Ed/1mIT/aTyGQOhmlfgtFnQf9RYVcoIiIiXRRk0C4BhnS6PDhyXZfWMbN4fMh+wjn3l/09iHPu\nAeABgIKCAnf4ZQdjY1ktxw3pF3YZ0hs5ByWL/UGMH/wJGqshLQcmzIaTbvYnEREROeIEGbQXAmPM\nbAQ+PF8BXLXHOnOAWyL9t08EqpxzW83MgIeBlc65XwVYY49obm1jS2UDsyerRVv2UFcOD5wBlZsg\nNgHGXwzTvwyDC8KuTERERA5TYEHbOddiZrcALwKxwCPOueVmdlNk+f3A88D5QCFQB1wfufkM4Bpg\nmZm9H7nue86554OqN0glFfW0tjnNCilQvh6WPQstDXDWD/x416NmwuATYPyFfmIZERERiQqB9tGO\nBOPn97ju/k7nHbDX7+LOubcAC7K2nvTR0H4acaRvqt4Ky/8Cy56BLe/560af3XFA40V3hVufiIiI\nBKJXHwwZLdonq9HBkH1IQzUkpEFMDCz4Dcy/GwYdC+fcAcd8GrKGHPw+RERE5IimoN0DNuysJTk+\nloHpiWGXIkEreQ8W/853D7nicd8t5MSbYMq1MGBM2NWJiIhID1LQ7gGFpTWMyknFH+MpUae12XcL\nefcB3zUkPgUmfhoy8v3yzPxw6xMREZFQKGj3gMLSGqaP7B92GdLdmuv9WNcAr/7EdxWZ9V8w+Qod\n1CgiIiIK2kGrbmhma1UDY3LTwi5FuoNzsPHNSOv1+/DVJRAbD194CTKH+IMbRURERFDQDlxhaQ0A\nY3LSQ65EDkv7jI0LH4IdqyA52/e7bmnwQTur985KKiIiIuFQ0A5Y4XYftMeqRfvI1NbmRw4pehee\n/xbkHQefvM+PHBKfFHZ1IiIi0ospaAdszfZdJMbFMLifhvY7YpSvh5V/hxVzYNjJcO6/w8gz4Yuv\nw1HHq3uIiIiIdImCdsDWltYwamAasTEKZ73epvnwxs9gwzx/Oe+4ji4hMTGQPyW82kREROSIo6Ad\nsMLSGk4Y3i/sMmR/2mdnBHjvMShdCWf9CCZeAv2GhVubiIiIHNFiwi4gmtU0tlBSWc+YXB0I2es0\n1cLCh+HeabBlib/u3P8PX1sKp35DIVtEREQOm1q0A9Q+4sjoHB0I2WuUrfMjhyx5AhqrfPeQlka/\nLHVAuLWJiIhIVFHQDtDa7bsAGKsW7d6hpREePNO3Zk+YDdNuhCEn6uBGERERCYSCdoAKS2tIiIth\naLZGHAlF7U7f73rT23D1MxCXCJc+ArmTID037OpEREQkyiloB2jN9l0acaSnOQcli+HdB2H5X6G1\nEYafCvUVkJINo88Ou0IRERHpIxS0A7S2tIYpQzXiSI9a+zL88TJISIMp18AJN0DO+LCrEhERkT5I\nQTsgdU0tFFfUc3nBkLBLiW4VG/3oIZlD4MQbYeQZcOGdMPFSSMoIuTgRERHpyxS0A9I+4oiG9gtA\nWxuse9V3D1n7EliMb7kGiEuAgs+HW5+IiIgICtqBWbu9PWhraL9uN/c2f5Bjag6c9m2Y+jnIzA+7\nKhEREZHdKGgHZG1pDQmxMQzTiCOHb+sHvvX6tG/7iWSOvwZGnA7jL/Yt2CIiIiK9kIJ2QJYWVzIm\nN424WE2+eUhaGmHFcz5gF78L8Skw5lwftIdM8ycRERGRXkxBOwDNrW0s2VzJ5SfoQMhD0tIIN56H\n5wAACnlJREFUdx8P1SWQPQo+8VM47ipIzgq7MhEREZEuU9AOwPIt1dQ3t3LC8OywSzkytLXCmheh\naAGcc4efWGb6V/ywfCPPhBj9KiAiIiJHHgXtACzcUA7ACcM1hvYBVW+FJX+AxY9CdTGk58EpX4fk\nfnDyLWFXJyIiInJYFLQDsHBjOcP6p5CTkRR2Kb3XmpfgT1eAa/Wt1uf9FI6eBbHxYVcmIiIi0i0U\ntLuZc45FmyqYOS4n7FJ6l/L18MGTMGAsTLoUhp4IJ93sh+brPyrs6kRERES6nYJ2N1u3o5by2iZ1\nGwFo3AXL/wbv/xE2zwfM972edCkkZcK5/x52hSIiIiKBUdDuZgs3tvfP7qMHQra1QkysP//k1bBh\nHvQfDWf9EI69QhPLiIiISJ+hoN3NFm4oZ0BaAiMGpIZdSs9pbfGB+sO/wJp/wC0LISUbTv8OzPw+\nDD4BzMKuUkRERKRHKWh3s4WbyikYlo31hWBZvh7eugtW/S/U7YTEDBh3ITTXAdkw/JSwKxQREREJ\njYJ2N9pW1UBReT3XnTQ87FKCUbMD1r4I/UbA8BngnG/FHnM2TLwERp8D8RppRURERAQUtLvVu5H+\n2dNGRFH/7I1vw7rXYP3rUPIe4GDq9T5o9x8F31kPcQlhVykiIiLS6yhod6N31peRkhDLhLyMsEs5\nNG2tsH05VBXDuPP9dXO/DmWFkD8Vzrjdj3U96NiO2yhki4iIiOyTgnY32bizlqcXF3P+xEHExR5B\nU4YXL/bdQYregeJF0FQDiZnwbxv86CGX/d6PFJKUGXalIiIiIkeUQBOhmZ1nZqvNrNDMbt/HcjOz\nuyPLl5rZlK7etjdxzvGD5z4kITaG754/Puxy9uYcVJXAmhfhzV/CM5+Hxhq/bNXf4Z+/gLoymHwF\nfPpB+PLbHUP05U5QyBYRERE5BIG1aJtZLHAvcA5QDCw0sznOuRWdVpsFjImcTgTuA07s4m17jblL\nt/Lm2p38+KIJ5IY17bpzUFcOFRugfAOMOA3Sc/2IIH/7CjRUdqybORR2bYXEMXDSrXDKNyDpCO3u\nIiIiItJLBdl1ZBpQ6JxbD2BmTwKzgc5heTbwmHPOAQvMLMvM8oDhXbhtr1BV38wdc1cwKT+Ta7pr\ntJG2NsD5VuXmeti+Ahqr/EyLDVVQUwpjzoW8Y/0Bis/eADXbfbePdpc/DuMvgvQ8OOaTkDsxctqj\nhTq1f/fULCIiIiK7CTJo5wNFnS4X41utD7ZOfhdv2yv88qXV/Lbx3xgXk0rsQzG+Zdm1waTLYMZX\nfVB+cKa/Htex/IQvwPQv+yHzHjgdWhqgpcn/bWuGc/8DTr7FH5j40My9Hzgl2wft5H6QNxnSciBr\nqB96L3uE/wuQP8WfRERERKRHHfEHQ5rZjcCNAEOHDu3Rx3bOkZIQR0bWQFIyUwEDi/GzILZ3xbAY\nyB4ZmRmx0/K0HL88PglGnQlxSRCbCHGRU/5UvzwjH656yk8Gk5ju7zd1IMQn++XZI+Cy3/Xo8xYR\nERGRgwsyaJcAQzpdHhy5rivrxHfhtgA45x4AHgAoKChwh1fyx2Nm3D5rHO68F/Y/xXhcIlzxxP7v\nJDEdZt+7/+UJKTD2E4dXqIiIiIj0uCBHHVkIjDGzEWaWAFwBzNljnTnAtZHRR6YDVc65rV28ba/R\nJ6ZbFxEREZGPJbAWbedci5ndArwIxAKPOOeWm9lNkeX3A88D5wOFQB1w/YFuG1StIiIiIiLdzfyA\nH9GhoKDALVq0KOwyRERERCSKmdli51zBwdY7gqYwFBERERE5cihoi4iIiIgEQEFbRERERCQACtoi\nIiIiIgFQ0BYRERERCYCCtoiIiIhIABS0RUREREQCoKAtIiIiIhIABW0RERERkQAoaIuIiIiIBCCq\npmA3sx3AphAeegCwM4THlZ6l97lv0Psc/fQe9w16n/uGsN7nYc65gQdbKaqCdljMbFFX5ruXI5ve\n575B73P003vcN+h97ht6+/usriMiIiIiIgFQ0BYRERERCYCCdvd4IOwCpEfofe4b9D5HP73HfYPe\n576hV7/P6qMtIiIiIhIAtWiLiIiIiARAQfswmNl5ZrbazArN7Paw65HuZ2ZDzOx1M1thZsvN7Law\na5LgmFmsmS0xs7lh1yLBMLMsM3vGzFaZ2UozOynsmqR7mdnXI9/XH5rZn8wsKeya5PCZ2SNmVmpm\nH3a6LtvMXjaztZG//cKscV8UtA+RmcUC9wKzgAnAlWY2IdyqJAAtwDedcxOA6cDNep+j2m3AyrCL\nkED9GviHc24cMBm931HFzPKBrwIFzrmJQCxwRbhVSTf5PXDeHtfdDrzqnBsDvBq53KsoaB+6aUCh\nc269c64JeBKYHXJN0s2cc1udc+9Fzu/Cb5Tzw61KgmBmg4ELgIfCrkWCYWaZwGnAwwDOuSbnXGW4\nVUkA4oBkM4sDUoAtIdcj3cA590+gfI+rZwOPRs4/CnyyR4vqAgXtQ5cPFHW6XIwCWFQzs+HA8cA7\n4VYiAbkL+A7QFnYhEpgRwA7gd5EuQg+ZWWrYRUn3cc6VAP8NbAa2AlXOuZfCrUoClOuc2xo5vw3I\nDbOYfVHQFukCM0sDngW+5pyrDrse6V5mdiFQ6pxbHHYtEqg4YApwn3PueKCWXvhTsxy6SB/d2fid\nqqOAVDP7bLhVSU9wfhi9XjeUnoL2oSsBhnS6PDhynUQZM4vHh+wnnHN/CbseCcQM4GIz24jvBjbT\nzB4PtyQJQDFQ7Jxr/1XqGXzwluhxNrDBObfDOdcM/AU4OeSaJDjbzSwPIPK3NOR69qKgfegWAmPM\nbISZJeAPtpgTck3SzczM8P05VzrnfhV2PRIM59x3nXODnXPD8f/Lrznn1AoWZZxz24AiMzs6ctVZ\nwIoQS5LutxmYbmYpke/vs9ABr9FsDnBd5Px1wHMh1rJPcWEXcKRyzrWY2S3Ai/ijmh9xzi0PuSzp\nfjOAa4BlZvZ+5LrvOeeeD7EmETl0twJPRBpI1gPXh1yPdCPn3Dtm9gzwHn7UqCX08pkDpWvM7E/A\nGcAAMysGfgT8DHjKzL4AbAI+E16F+6aZIUVEREREAqCuIyIiIiIiAVDQFhEREREJgIK2iIiIiEgA\nFLRFRERERAKgoC0iIiIiEgAN7yciEkXMrD/wauTiIKAVP+04QJ1zTpN3iIj0EA3vJyISpczsx0CN\nc+6/w65FRKQvUtcREZE+wsxqIn/PMLN5Zvacma03s5+Z2dVm9q6ZLTOzUZH1BprZs2a2MHKaEe4z\nEBE5sihoi4j0TZOBm4Dx+NlPxzrnpgEP4WdPBPg1cKdz7gTgksgyERHpIvXRFhHpmxY657YCmNk6\n4KXI9cuAMyPnzwYmmFn7bTLMLM05V9OjlYqIHKEUtEVE+qbGTufbOl1uo2PbEANMd8419GRhIiLR\nQl1HRERkf16ioxsJZnZciLWIiBxxFLRFRGR/vgoUmNlSM1uB79MtIiJdpOH9REREREQCoBZtERER\nEZEAKGiLiIiIiARAQVtEREREJAAK2iIiIiIiAVDQFhEREREJgIK2iIiIiEgAFLRFRERERAKgoC0i\nIiIiEoD/A7Q+XCnND+UPAAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fig, ax = plt.subplots(figsize=(12, 6))\n", "ax.plot(tlist, avg.expect[0],\n", " tlist, avg.expect[1],'--')\n", "ax.set_xlabel('Time')\n", "ax.set_ylabel('Photocount rates')\n", "ax.legend(('Cavity ouput', 'Spontaneous emission'));" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Versions" ] }, { "cell_type": "code", "execution_count": 14, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
SoftwareVersion
QuTiP4.2.0
Numpy1.13.1
SciPy0.19.1
matplotlib2.0.2
Cython0.25.2
Number of CPUs2
BLAS InfoINTEL MKL
IPython6.1.0
Python3.6.1 |Anaconda custom (x86_64)| (default, May 11 2017, 13:04:09) \n", "[GCC 4.2.1 Compatible Apple LLVM 6.0 (clang-600.0.57)]
OSposix [darwin]
Wed Jul 19 21:56:31 2017 MDT
" ], "text/plain": [ "" ] }, "execution_count": 14, "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.1" } }, "nbformat": 4, "nbformat_minor": 1 }