{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# QuTiP example: Calculate the quasi-steadystate of a time-dependent (period) quantum system" ] }, { "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": "markdown", "metadata": {}, "source": [ "Find the steady state of a driven qubit, by finding the eigenstates of the propagator for one driving period" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "%matplotlib inline" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "import matplotlib.pyplot as plt" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "import numpy as np" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "from qutip import *" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [], "source": [ "def hamiltonian_t(t, args):\n", " #\n", " # evaluate the hamiltonian at time t. \n", " #\n", " H0 = args['H0']\n", " H1 = args['H1']\n", " w = args['w']\n", "\n", " return H0 + H1 * np.sin(w * t)" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [], "source": [ "def sd_qubit_integrate(delta, eps0, A, w, gamma1, gamma2, psi0, tlist):\n", "\n", " # Hamiltonian\n", " sx = sigmax()\n", " sz = sigmaz()\n", " sm = destroy(2)\n", "\n", " H0 = - delta/2.0 * sx - eps0/2.0 * sz\n", " H1 = - A * sx\n", " \n", " H_args = {'H0': H0, 'H1': H1, 'w': w}\n", " # collapse operators\n", " c_op_list = []\n", "\n", " n_th = 0.5 # zero temperature\n", "\n", " # relaxation\n", " rate = gamma1 * (1 + n_th)\n", " if rate > 0.0:\n", " c_op_list.append(np.sqrt(rate) * sm)\n", "\n", " # excitation\n", " rate = gamma1 * n_th\n", " if rate > 0.0:\n", " c_op_list.append(np.sqrt(rate) * sm.dag())\n", "\n", " # dephasing \n", " rate = gamma2\n", " if rate > 0.0:\n", " c_op_list.append(np.sqrt(rate) * sz)\n", "\n", "\n", " # evolve and calculate expectation values\n", " output = mesolve(hamiltonian_t, psi0, tlist, c_op_list, [sm.dag() * sm], H_args) \n", "\n", " T = 2 * np.pi / w\n", "\n", " U = propagator(hamiltonian_t, T, c_op_list, H_args)\n", "\n", " rho_ss = propagator_steadystate(U)\n", "\n", " return output.expect[0], expect(sm.dag() * sm, rho_ss) * np.ones(shape(tlist))" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [], "source": [ "delta = 0.3 * 2 * np.pi # qubit sigma_x coefficient\n", "eps0 = 1.0 * 2 * np.pi # qubit sigma_z coefficient\n", "A = 0.05 * 2 * np.pi # driving amplitude (sigma_x coupled)\n", "w = 1.0 * 2 * np.pi # driving frequency\n", "\n", "gamma1 = 0.15 # relaxation rate\n", "gamma2 = 0.05 # dephasing rate\n", "\n", "# intial state\n", "psi0 = basis(2,0)\n", "tlist = np.linspace(0,50,500)" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [], "source": [ "p_ex, p_ex_ss = sd_qubit_integrate(delta, eps0, A, w, gamma1, gamma2, psi0, tlist)" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAtUAAAGDCAYAAADkjOwcAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzs3Xl4JGd9L/rvr1t7q1v7Ls2MZpVmsY092IY4wMEEMEtI\n8pAFsgBJLsc5BHJzksPJyU1ycpKQHHITIGASH0IIYAOGa8LuFbwbLzPj2Wc02jXa1a2l1Yt6f+8f\nVdWq7q6q7qFHeGb8/TwPDx7pVauq1F31rV+9iyilQEREREREPznXy70BRERERERXO4ZqIiIiIqIy\nMVQTEREREZWJoZqIiIiIqEwM1UREREREZWKoJiIiIiIqE0M1EV3zRORBEXnfZXy9XxeRRy7X620V\nEXm/iDzzE/7sG0RkxuH7d4vIn5fS9nISkX0ickJEQiLykcv82n8pIvc6fP+siLzhcv5OIrp2VLzc\nG0BEZCYikwA6AKRNX/6iUur3f9LXVErdYXr99wP4XaXUbSVuzw4AEwAqlVIp/fW+AuArP+n2XAuU\nUnfafU//G/6uUuqHW/CrPwrgcaXUDVvw2o6UUgeM/xaRvwSwWyn1Gz/t7SCiKxNDNRFdid65RYHs\nmiIibqVUunjLa8p2APe93BtBRJSP3T+I6KohIv8iIt80/fvjIvIjERH93+/Suwasi8iYiLxV//oT\nIvK7IjII4G4ArxGRsIis6d9/u4gc139uWq9CGp7S/39N/5nX5HerEJHXisgREQnq//9a0/eeEJG/\nFpFn9S4Lj4hIq83+vUFEZkTkT0UkICKTIvLrpu9/UT8GD4hIBMB/EpEGEfmyiPhFZEpE/kxEXLkv\nK3fp2zYkIrebvvEBETmvb9e4iPxni21y2pa/sWh/D4BtAL6nH6+PisgPROTDee1Oicgv2hyHn9e7\nWqzpx29Q//pjAP4TgLv0195r8bP9IvKkvk+P6vt+r/n45rWfFJE3mb5UIyJf13/+JRG5Pr+t/r76\nUwC/qm/HSav9IKJXFoZqIrqa/BGAQ3qo/VkAvwPgfUopJSI3A/gygP8GoBHA6wBMmn9YKXUewJ0A\nnlNK1SulGvVvRQD8lv5zbwfweyLyC/r3Xqf/f6P+M8+ZX1NEmgH8AMCnAbQA+ASAH4hIi6nZewF8\nAEA7gCoAf+ywj50AWgH0AHgfgM+JyL681/oYAC+AZwB8BkADgJ0AXq/vxwdM7W8BMKa/5v8E8B/6\nNgPAEoB3APDpP/NJEbnxEralgFLqNwFchPa0oV4p9fcAvgQg201CD6o90I5bDj0ofw3A/w2gDcAD\n0AJ6lVLqjQCeBvD7+msPW2zCVwEc07f7r/XtvhTvAvD/AWjWX+vbIlKZt48PAfhbAF/Xt+P6wpch\nolcahmoiuhJ9W69SGv/7vwBAKRUF8JvQguu9AD6slDIqj78D4AtKqUeVUhml1KxSaqiUX6aUekIp\ndVr/uVPQQt3rS9zWtwMYUUrdo5RKKaW+BmAIwDtNbf5dKTWslNoA8A0AxfoD/7lSKq6UehJa8PwV\n0/e+o5R6VimVAZAE8GsA/odSKqSUmgTwj9COkWEJwKeUUkml1NcBXNC3GUqpHyilxpTmSQCPAPjZ\nS9iWUn0XwF4R2aP/+zehBdKERdtfBfAD/e+YBPAPAGoBvNaibQ4R2Qbg1aZtfgrA9y5xW48ppe7X\nf/cnANQAuPUSX4OIXoEYqonoSvQLSqlG0//+1fiGUuoFAOMABFpANfRBq8heMhG5RUQe17tQBKFV\nsy27aFjoBjCV97UpaJVYw4Lpv6MA6h1eb1UpFcl7rW7Tv6dN/90KoDLv9+f/7lmllLJ6PRG5Q0Se\nF5EVvSvM25C738W2pSRKqRiArwP4Db1rynsA3GPTPOd46jcP03n7ZKfbZpsvRfb46r97Bj/BPhPR\nKw9DNRFdVUTkQwCqAcxBmwnCMA1gVwkvoSy+9lVo1dQ+pVQDtH7X4tDebA7a4DmzbQBmS9gWK00i\n4sl7rTnTv83bE4BWrd6e1978u3tERPK+Pyci1QC+Ca0S3KF3hXkAm/tdyrbYsTpmXwLw6wBuBxDN\n70ZjknM89W3vQ2nHc95mmw0RAHWm13ZD62Ji1mf6vgtAL6z3udj7goheYRiqieiqofe3/Rto/XN/\nE8BHRcToSvFvAD4gIreLiEtEekRkwOJlFgH0ikiV6WteACtKqZjeN/u9pu/5AWSg9Vm28gC0rg3v\nFZEKEflVAPsBfP8n3U8A/0tEqvR+4++A1se3gD7zxzcAfExEvCKyHcB/hdY1xtAO4CMiUikivwxg\nUN/mKmg3J34AKRG5A8Cbf9JtybOIvOOlh+gMtO4pdlVq6Pvzdv3vWAmtH30cwI+L/VKl1BSAo6Zt\nvg253XCGoQ1EfLv+2n8G7RiY3SQivyQiFdD6dccBPG+zjzvyBoUS0SsYTwZEdCUyZo4w/vctPeTc\nC+DjSqmTSqkRaDMw3CMi1UqpF6EPtgMQBPAkCivIAPAYgLMAFkQkoH/tvwD4KxEJAfgLmLqV6P24\nPwbgWb1/d07/WqXUMrSw+UcAlqFVz9+hlArgJ7MAYBVadfQrAO4s0jf8w9AqsOPQBi5+FcAXTN9/\nAcAeaFXtjwF4t1JqWSkVAvARfV9Xod1IfLfMbTH8HYA/04+XeVDmlwEcQm7oz6GUugDtpukz+ja/\nE9qgR6v+11beC21w5gq0gZlfNr12ENrf+vPQKt8RaN07zL4DrV/3KrQbt1/S+1fnM24ulkXkpRK3\njYiuYZLb1Y6IiF4uoq3Wd69Sqvfl3patICK/BeCDpS68c5l+51+Ci7QQ0U8BK9VERLTlRKQOWpX4\ncy/3thARbYUtDdUi8gURWRKRMzbfFxH5tIiMirYQwI1W7YiI6OolIm+B1nd7EVr3FCKia86Wdv8Q\nkdcBCAP4slLqoMX33watP+DboPWB+yel1C1btkFERERERFtgSyvV+sT7Kw5N3gUtcCul1PMAGkWk\nayu3iYiIiIjocnu5+1T3IHchgxmUNsE/EREREdEVo+Ll3oBSicgHAXwQADwez00DA1bTzxIRERER\nXT7Hjh0LKKXyF4oq8HKH6lmYVq+CtnKV5apZSqnPQR81fvjwYXX06NGt3zoiIiIiekUTkalS2r3c\n3T++C+C39FlAbgUQVErNv8zbRERERER0Sba0Ui0iXwPwBgCtIjIDbXWrSgBQSt0NbanctwEYBRCF\nthoaEREREdFVZUtDtVLqPUW+rwB8aCu3gYiIiIhoq73c3T+IiIiIiK56DNVERERERGViqCYiIiIi\nKhNDNRERERFRmRiqiYiIiIjKxFBNRERERFQmhmoiIiIiojIxVBMRERERlYmhmoiIiIioTAzVRERE\nRERlYqgmIiIiIioTQzURERERUZkYqomIiIiIysRQTURERERUJoZqIiIiIqIyMVQTEREREZWJoZqI\niIiIqEwM1UREREREZWKoJiIiIiIqE0M1EREREVGZGKqJiIiIiMrEUE1EREREVCaGaiIiIiKiMjFU\nExERERGViaGaiIiIiKhMDNVERERERGViqCYiIiIiKhNDNRERERFRmRiqiYiIiIjKxFBNRERERFQm\nhmoiIiIiojIxVBMRERERlYmhmoiIiIioTAzVRERERERlYqgmIiIiIioTQzURERERUZkYqomIiIiI\nysRQTURERERUJoZqIiIiIqIyMVQTEREREZWJoZqIiIiIqEwM1UREREREZWKoJiIiIiIqE0M1ERER\nEVGZGKqJiIiIiMrEUE1EREREVCaGaiIiIiKiMjFUExERERGViaGaiIiIiKhMDNVERERERGViqCYi\nIiIiKhNDNRERERFRmRiqiYiIiIjKxFBNRERERFQmhmoiIiIiojIxVBMRERERlWnLQ7WIvFVELojI\nqIj8icX3G0TkeyJyUkTOisgHtnqbiIiIiIgupy0N1SLiBvBZAHcA2A/gPSKyP6/ZhwCcU0pdD+AN\nAP5RRKq2cruIiIiIiC6nra5U3wxgVCk1rpRKALgPwLvy2igAXhERAPUAVgCktni7iIiIiIgum60O\n1T0Apk3/ntG/ZnYXgEEAcwBOA/gDpVRmi7eLiIiIiOiyuRIGKr4FwAkA3QBuAHCXiPjyG4nIB0Xk\nqIgc9fv9P+1tJCIiIiKytdWhehZAn+nfvfrXzD4A4D+UZhTABICB/BdSSn1OKXVYKXW4ra1tyzaY\niIiIiOhSbXWoPgJgj4j064MPfw3Ad/PaXARwOwCISAeAfQDGt3i7iIiIiIgum4qtfHGlVEpEfh/A\nwwDcAL6glDorInfq378bwF8D+KKInAYgAP67UiqwldtFRERERHQ5bWmoBgCl1AMAHsj72t2m/54D\n8Oat3g4iIiIioq1yJQxUJCIiIiK6qjFUExERERGViaGaiIiIiKhMDNVERERERGViqCYiIiIiKhND\nNRERERFRmRiqiYiIiIjKxFBNRERERFQmhmoiIiIiojIxVBMRERERlYmhmoiIiIioTAzVRERERERl\nYqgmIiIiIioTQzURERERUZkYqomIiIiIysRQTURERERUJoZqIiIiIqIyMVQTEREREZWJoZqIiIiI\nqEwM1UREREREZWKoJiIiIiIqE0M1EREREVGZGKqJiIiIiMrEUE1EREREVCaGaiIiIiKiMjFUExER\nERGViaGaiIiIiKhMDNVERERERGViqCYiIiIiKhNDNRERERFRmRiqiYiIiIjKxFBNRERERFQmhmoi\nIiIiojIxVBMRERERlYmhmoiIiIioTAzVRERERERlYqgmIiIiIioTQzURERERUZkYqomIiIiIysRQ\nTURERERUJoZqIiIiIqIyMVQTEREREZWJoZqIiIiIqEwM1UREREREZWKoJiIiIiIqE0M1EREREVGZ\nGKqJiIiIiMrEUE1EREREVCaGaiIiIiKiMjFUExERERGViaGaiIiIiKhMDNVERERERGViqCYiIiIi\nKhNDNRERERFRmRiqiYiIiIjKxFBNRERERFQmhmoiIiIiojJteagWkbeKyAURGRWRP7Fp8wYROSEi\nZ0Xkya3eJiIiIiKiy6liK19cRNwAPgvg5wDMADgiIt9VSp0ztWkE8M8A3qqUuigi7Vu5TURERERE\nl9tWV6pvBjCqlBpXSiUA3AfgXXlt3gvgP5RSFwFAKbW0xdtERERERHRZbXWo7gEwbfr3jP41s70A\nmkTkCRE5JiK/ZfVCIvJBETkqIkf9fv8WbS4RERER0aW7EgYqVgC4CcDbAbwFwJ+LyN78Rkqpzyml\nDiulDre1tf20t5GIiIiIyNaW9qkGMAugz/TvXv1rZjMAlpVSEQAREXkKwPUAhrd424iIiIiILout\nrlQfAbBHRPpFpArArwH4bl6b7wC4TUQqRKQOwC0Azm/xdhERERERXTZbWqlWSqVE5PcBPAzADeAL\nSqmzInKn/v27lVLnReQhAKcAZAB8Xil1Ziu3i4iIiIjochKl1Mu9DZfs8OHD6ujRoy/3ZhARERHR\nNU5EjimlDhdrdyUMVCQiIiIiuqoxVBMRERERlYmhmoiIiIioTAzVRERERERlYqgmIiIiIioTQzUR\nERERUZkYqomIiIiIysRQTURERERUJoZqIiIiIqIyMVQTEREREZWJoZqIiIiIqEwM1UREREREZWKo\nJiIiIiIqE0M1EREREVGZGKqJiIiIiMrEUE1EREREVCaGaiIiIiKiMjFUExERERGVqeRQLSJvsvja\n+y7v5hARERERXX0upVL9FyLyLyLiEZEOEfkegHdu1YYREREREV0tLiVUvx7AGIATAJ4B8FWl1Lu3\nZKuIiIiIiK4ilxKqmwDcDC1YxwFsFxHZkq0iIiIiIrqKXEqofh7AQ0qptwJ4NYBuAM9uyVYRERER\nEV1FKi6h7ZuUUhcBQCm1AeAjIvK6rdksIiIiIqKrx6VUqqdF5DdE5C8AQES2AYhtzWYREREREV09\nLiVU/zOA1wB4j/7vEIDPXvYtIiIiIiK6ylxK949blFI3ishxAFBKrYpI1RZtFxERERHRVeNSKtVJ\nEXEDUAAgIm0AMluyVUREREREV5FLCdWfBvAtAO0i8jFoc1X/7ZZsFRERERHRVaTk7h9Kqa+IyDEA\ntwMQAL+glDpvfF9EmpRSq1uwjUREREREV7RL6VMNpdQQgCGbb/8IwI1lbxERERER0VXmUrp/FMPV\nFYmIiIjoFelyhmp1GV+LiIiIiOiqcTlDNRERERHRKxK7fxARERERlanoQEURqQFwJ4DdAE4D+Del\nVMqi6e2XeduIiIiIiK4KpVSqvwTgMLRAfQeAf7RqpJRauYzbRURERER01ShlSr39SqlDACAi/wbg\nxa3dJCIiIiKiq0spleqk8R823T6IiIiIiF7RSqlUXy8i6/p/C4Ba/d8CQCmlfFu2dUREREREV4Gi\noVop5f5pbAgRERER0dWK81QTEREREZWJoZqIiIiIqEwM1UREREREZWKoJiIiIiIqE0M1EREREVGZ\nGKqJiIiIiMrEUE1EREREVCaGaiIiIiKiMjFUU8nC8RQmApGXezOIiIiIrjgM1VSyT/9oBG/55FM4\nP79evDGAYDSJf3tmAql0Zou3jIiIiOjlxVBNJRtaCCGRzuAP7juORKp4UL7/pRn89ffP4UdDSz+F\nrSMiIiJ6+TBUU8kmAxF4qyswvBjG2blg0fYnp9cAAPcfm9nqTSMiIiJ6WTFUv4Itrsfw/n9/EW/9\n1FNQSjm2TaQymFmN4rW7WwAA06sbRV//5IwWqh8fWkIgHC9pm07PBJHOOG8LERER0ZWGofoV7H/8\nx2k8ccGPoYUQAuGEY9vp1SgyCrhtT5v275WoY/uVSAJTy1G864ZupDIKT1zwF92e0zNBvPOuZ3D/\nsenSd4KIiIjoCsBQ/Qo27g/DV1MBAJhadp7VY1Kf9WN/lxet9dVFQ7VRpX73Tb1wCXCxSHsAePjs\nAgDggdMLRdsSERERXUm2PFSLyFtF5IKIjIrInzi0e7WIpETk3Vu9TQQopbCwHsPN/Vp3jmJT5Rnf\n39HiQV9zbdGQfHJ6DS4BbtzWhE5fDWZWi4fqH55fBAD8eCyA9ViylN0gIiIiuiJsaagWETeAzwK4\nA8B+AO8Rkf027T4O4JGt3B7atB5LIZbM4KbtTXC7BFPLzqF3cjkCb00Fmj1V6Guqw3SRkDwZiKCr\noRae6gr0NtVhZsW5D/b0ShRDCyHccbATybTC4yXOGDK8GGIAJyIiopfdVleqbwYwqpQaV0olANwH\n4F0W7T4M4JsAOPfaT8niegwA0NtUi96mWkwW6f4xtRxFf6sHIoK+5lrMrcUc559eXI+js6FG+x3N\ntUUr1c+NLQMA/vDn9qKhthLPj68U3YfVSALv/Mwz+N8PDhVtS0RERLSVtjpU9wAwjzqb0b+WJSI9\nAH4RwL84vZCIfFBEjorIUb+/+KA3cmaE6g5fDba3eIqG6oVgDN0NtQCAbc11SGcU5oMxx9fv9Omh\nuqkOC+sxx7mtL65E4RKgv9WDnW2eon28AeCbL80gnsrg4TMLnDGEiIiIXlZXwkDFTwH470opx9VE\nlFKfU0odVkodbmtr+ylt2tUnlc6UFDAX9EDc6avBjpY6TAWijtPqLUcSaKmvAgD0NdUBsB98aPTX\n7siG6lpk1ObvtDK7toFOXw0q3S5sb64r2h1FKYX7jkyjptKF5UgCx6ZWHdsTERERbaWtDtWzAPpM\n/+7Vv2Z2GMB9IjIJ4N0A/llEfmGLt+ua9Qf3ncCd9x4r2m4ppM0b3e6rxo4WD0LxFJYj1tPqpTMK\nq9EEWjx6qG7WQrVdl45QPIVoIo3OhmoAWqh2ag8As6sb6NXD+rYWD+aDG4in0rbthxfDGF0K44/f\nvA9VFS48dKa0GUO++OwE7nl+qqS2RERERKXa6lB9BMAeEekXkSoAvwbgu+YGSql+pdQOpdQOAPcD\n+C9KqW9v8XZdk6KJFB49t4gXxpeLLuayEIyhobYSNZVubNNDst00eavRBJQCWuq1kGxUoO26fywG\nYzntjMq20+DG2bUN9Ojhe3tzHTIKmHFYYGbcHwYA3LqzBbf0N+P58WXbtoaVSAJ/++AQ/v6hIcSS\n9oGdtl6piwEZvvLCFM7Pr5fc/q7HRvDZx0dLbn/fixfxvi+8WHI3omNTK3jHZ54ueT+Ww3H8yt3P\n4dxcafuglMJ/v/9UyQN2AeCe56fwjaOlz/H+/Pgy/s+TYyW3n1vbwD88fAFJh7EUZrFkGp945AKW\nL+Fv/YVnJi7p7/zg6Xk8fqH0Y3T84iq+9uLFktvPBzdw12MjjuNHzGLJND7zoxGsRZ3n/TcopXDP\n81MYXQqVvE2PDS2WdL4zDC2s48HT8yW3Xw7H8Y0j00WvIYZkOoMHT8+X/L4AtBmeghulDzAfXQqX\nNIOUYTkcx5h+jShFLJnGyGLpfwOlFMb84ZKPEQAsrccu6RhF4qlLuk6l0hnHQpQVdpu8/LY0VCul\nUgB+H8DDAM4D+IZS6qyI3Ckid27l734lem5sGYl0BuuxlGN/ZyC3z7MxoHBx3frit6wvDGN0/6iq\ncKG1vtq2O4fxOkao7mqogdsltiE5lc5gYT2GnkYtVO9o1buXOHQBmdJvALa31GFPuxcTgQgyRU4Q\nX3vxIhKpDEKxFJ4o8UKcSmc4u0gR/lAcn3181LHPvNnoUgg3f+yH+MoLpT0xmFvbwP/zrTP4w6+f\nKOkisJFI47OPj+ETjw5n51d3opTCvz49jieH/XjwTGnh49+fncSZ2XV86ceTJbX/5kszeHFyBXc9\nPlJS+xcnVvD1o9P43w8OlXThDkaT+Jvvn8Nff+8cIvFU0fZKKfzld8/i7x4cwuhSaeHjk48O467H\nR/HoucWS2n/9yDQ+/dgovljiMToyuYK/+v45/MPDF0pq7w/F8YffOIE//Y/TRT/7gBb+/us3TuJP\nv3UaS+vO50fD//zOWfzDI8N4dqy0EPvPj4/iHx8dxtdeLO3m5qEzC/jzb5/Bp35Y2vviwkIId97z\nEv70W6dLar8aSeD9XziCD3/teEnnsWQ6g/98zzF89JuncHx6rWh7pRT+4jtn8HtfeQlfP1LaPv/7\nsxN477++gL9/qLQB5o8PLeFtn34af3DfiZLan5kN4q3/9DR+8bPPlhRiZ1aj+KV//jHe8qmnSjpf\nBDeS+O0vHsHt//hkSZ+FVDqDP/nmKdzydz/CvzxR2k3sP/1wBIf/5of4o2+cLKn9N45M49a/+xHe\n9k9Pl3S+eOjMPH727x/DDX/1CFZsnlCbnZhew5s+8SQO/MVDODVT/H2xuB7Dmz/5JF71V4/gOyfy\nOycUSqUz+Pm7nsHhv3kUn/lRaZ+FK9WW96lWSj2glNqrlNqllPqY/rW7lVJ3W7R9v1Lq/q3epmuV\nedXCoQXnas/iegztPq3ybPz/Usj6QrMc0UJyi6c6+7Wuhhrb4L6wvtlfGwAq3C59rmrrUL2wHkM6\no7KV6m3NHgDOC9JMLUfQ4qmCt6YSO9s82Eims7/XilIK9z4/hdfsbEFrfTW+fXzOtq3Z3z4whDf+\nw5PYSFy9lW2l1CVVJL53cg4f+spLJYfkr7wwhf/34QslV4YfPruIjAI+/uBQSZVe48I1tBDCN4/N\nFG3/5PASNpJppDMK/1TCCfr8fAhj/gjcLsFdj40WDWihWBKPnluES4AvPzdVNMQqpXC/vt0PnVko\nunASgGxAubAYKmkmHGPQbiiewndPFn9vvzixgqEFrTJ3XwmV20A4ju+cmMvZNiepdAaff2Zc27Zj\nMyW9/4z3zxPD/pLeF597agyxZAbzwRhemCh+jO4/NoOJQARKAT8ooXL7wvgyHtHfez84VfyYTq9E\ncfdT2j6XcnMWiafwl987C0ALjsWqkkopfPT+k0ikMxj3R0q6Gfqr75/DwnoMqYzCU8PFB/h/7qlx\nHJ1ahQjw2PnihYeHzy7iay9OwyUo6anK6FII/+t751DhEjxxwV80AMaSaXzoqy8BSnvKUMoTgI/e\nfwrBjSTWYykcv1g8AH7ikWGM+cPIKOCpkeLH6L4XL+LxC35UV7jweAkrBT87tqyN/6lw47ESjtHS\negyf/OEwAOCJC0tFbwyUUvj4Q0OIxNMY80cwWWQ8EgDc+/xFBEIJhGKpkp56/ODUHC4uR7GRTOOH\nJdxIPDMSwPBiGLFkBt87WfyzMLwYxqmZICLxNL5+CU/brkRXwkBFukyeHPbj1p3NALSg4GRxPZ6t\nJLd6quF2iW3lOb9SDWih2r5SrYdqvQIOaP2q7R7fGWHbqFS31lehrsqdrUZbmQxEsb1Fq2jvbNNC\n+LjfPoT7w3HMB2N484EOvP1QJx6/sFT0kW44nsLXj1zUA0Xxu21Au7D+2bdPI5ooXi00zK5tXNJj\nxKeG/ZgPOs/7bfbBe47hNz7/QsmPsL/440n84PQ8PqWf2IsxLiyffXy0pMfYPzy/iN6mWkQTafzb\nMxNF2z96bhE7Wz040O0rqbr94JkFNNVV4tdv2Ybvn5or+kj0e6fm4HYJ/vjN+zC0EMK5It0PHjqz\ngHgqgz+5YwDBjWTRatXp2SCGF8P4yBt3Q0SKdtFYjyXxg9PzePdNvWisqyxpn7/64kXc0NeIgU4v\n7i1hzMA9z0+hobYSbxxox/0vzRQNdPe9eBGJdAbvuK4LT434Mbvm/P57bGgJ0ysb+PnruzEXjGWn\nzLQzEYjgiQt+vOuGbqQzCt894Rxik+kM7n3+Iu442AlPlbukz+c9z03hut4GDHb5SrrxuO/INJo9\nVbjjYCcePrtY9CbzgdPzSKQyeN9rtuPUTLDozdORyRUsrsfxgZ/ZgUgijWdHA47tl0JxnJwJ4ndu\n6wcAPHLOeRyJUgpPXFjCL93Yg2ZPVUlh6NnRAK7rbcDNO5qzi3E5eXFiBTWVLrzn5m14dixQ9H30\nkh5y3//aHZhd2yh6YzC6FEY0kcZvvWY7Mgp4psgxiiXTuLAYwntv3ga3S/B0CSH59GwQP7unDX3N\ntXh6xPn1AeDM3Dp6Gmvxur1tRf9mgFY5B4BffXUfTs2sFX1icGZus30kkcZp/eftzAdjWI4k8Kuv\n1oavFQv/lcfoAAAgAElEQVTJSimcng3iHdd1obbSjRdKCNUnZ4LY3+3DwZ6Gkm5gT82soa7Kjbcd\n6sLRqZWihQpjBeZfPtyLmdWNS7q+XWkYqq8RqXQG06tR3Nzfgt6mWsd+iUop+MNxtHu1yrPLJWj3\nVjt0/zAq1bmh2u6Nb+6vbehtqsO0zQIws0ao1ivVIoJtRWYAubgSxfYWLUzvaqsHAIwH7E/QRuDe\n1VaPQ72NiKcyjqEdAL59fBaRRBotnip86bmpkoLv558ex73PX8SXflxa14ahhXXc9vHHSgqXgPY4\n9wNfPILf/uLRkirJ67EkHh9awnPjy/h0CVXbtWgCxy+uoqG2Ev/y5FjRfoyBcBynZtbwm7duRyqj\n8Og550pMIBzHiek1/PJNfTjQ04ATRSpJwY0knh9fxpsPdOJndrfi/ELIsXKTzig8dn4JP7e/A6/Z\n1YJkWmFk0fnC/cQF7Wb0joOdAICzc84XsWdHA+jwVeO3f6YfVRWuou2NQPn+n+nHvg4vTs04tz85\nvYZ4KoNffFUPXr+3DS8VmdlmJZLA6FIYbz/UhXdc14Wzc+sIFblwH5taxRsH2vErh/uwFk3ibJG+\n3scvrmFfhxd/9OZ9UApFu0+dmgnC7RL87S8dgqfKXTQAGsHjztfvwmCXDw+ddW4/GYhgI5nGmw90\n4M0HOvHQ2QXHz2cqncHoUhiv2dmCn7++G8cvrjnORgRoXS0O9TTg3Tf1IriRxHNFwseFxRA6fNX4\nndt2AgAeLrIPw3of3t97wy54qyuKDra+oD9ZeNNgB67rbcDDZ51Drz8cx2o0iUM9DXjjQDseG3Ku\neiqlMLQQwmCnD7cPtmNoIVT05mloYR37Orx40/4OxJKZooFuaD6EmkoX3vfaHQByn65aMW5wf+3m\nPjTUVuLJIu1HFsNIZxRu6W/GDX2NRavzsWQaY/4w9nf7cNvuNjw3tly0Mnx2LogD3T7ctrsVF1ei\njt0UjfbbmuvwlgOdyCjgxSJPns7MrkME2ZunYjekRuh+5/XdaPdWF20/vbKB4EYSr9rWhMM7moqG\n5HRG4exsENfrN1vHp9eK3jydnAniYE8Dbt3ZjLVoEiNFbp5OzayhobYS776pFwBwZPLqnc2Lofoa\nEdxIaoMJPVUY7PJlH+1aCcVTSGcUmk0hucNXY9v9YyWSgEuAxrrN9p0NtViPpSwffZv7axv6mmux\nGIpZVg2NE7dRqQa05dDtun/EkmnMBTeylep2bzU8VW7HSrUxaGVnmwd7O7QQXmxgyrePz2Kg04s/\n/Lm9OD+/XvTEkEhlshWw//PUWNFgA2hdLZQCPvHoMOaKXMAA4InhJaQzCufn13HXY8VD8jMjAaQy\nCvs6vLj7yfGiVdunRwLIKOC/vUULTy8VCb1PDfuhFPArh/vQ11ybDUdO26MUcPtgOw52+3BmLugY\nhs7MBpHKKNy2uxUHun1IpDKO1a25tQ2E4inc0NeEA90NAJxDciajMBEIY7DTh23NdfBWV+DMrHPA\nHPNHsLfDiwq3C/s6vI6fNa19GK31VWj2VGGgy1u0a5bxPt7TXo+BTh/mgjHHx97GoN3dHfXY3+0D\nAMdtCse1MRe72+txsEdrX2xw4Jg/jN3t9djRoh2joSJPwsb8YWxvrkN9dQUGu3xFX39kKZydp/6m\n7Y04P7/u+L4wPot72r24aXsT1qJJxwA4tRJFIp3Bng4vbu7XnuY5VQDTGW0g2t6OetyyswUAcKpI\nH+PhxRD2dnixraUOfc21RbseDC+G0eatRru3BrfuasGxi85Bwgjh+zq9eP3eNpyecQ43Rgjf16G1\nX4+lHP9u/nAcK5EE9nV68YZ97QCcA50Rwgc6fXjNzhZUVbiKVm6NEN7XXIc97fV4ukj78/PrqK10\no7+1Hrftbi1aqT43r/1NB7t8+Nk9rTg1G3T87FxYCCGjgP16+3A8hZMOf+dIPIWJQAT7u334md3a\n+6LYNp2dW8eBbh9eta0R1RUuPDvm3P7MbBD9rR70NddhoNNbNCSfnQ3CJdo+3LqzBc8Xmajg1Ky2\nf9f1NuCW/mYMLYQc+1VPBMKIJNI41NuIm/ubkUhlHAsDiVQG5+bXcX1vA27p147Ri5POwf3EdBDX\n9zVif5cPnio3jpRQDb9SMVRfI4wPRZOnCvs6tIF7do/71yJa2Guorcx+rcNXne22kS8QSaDZUwW3\nS7Jf69K7dlj1Yzb31zb0NtVBKWB+rbD97OoGWuurcyrb21u0yrZVX8yZ1SiU0oI3oFW2d7bVO472\nHvdHUFPpQndDLXa3a6F62KGCaVwwbu5vzl6Ei1Ukn7iwhNVoEh+5fQ/Wokk8eLr449kHTi9gsMuH\nVEbh808Xr1b/8NwS2r3VeMO+Nny7yCNyQOvn6KupwIfeuBuJdKZoGHr8whIa6yrxy4d7UV3hKhok\nXhhfQbOnCge6fTjU01D0UeWFxRAqXIJ9nV4c7GlAKJayne8c2AyMu9q17h8AHKuq4/pAo51tnmyo\nc2o/F9xALJnBzrZ6uFyC/XrQt6OUwrg/jJ2t2ntvoNNbNDCO+yPY2aq95/Z3+bC4Hne8iI35w/BW\nV6DNW43BLi8A55BsvO93tdZjsKt4SJ7IPrXxoKexFr6aCsf28VQaF1ei2NWmrag60FV8n8f8YezU\nnyANdHkxNB9yvNCPLoWwvcWDmko3Brt8CMVSjrP/jCyGIaI9edrcZ/tjZNxA72mvx75O/Zg67MPF\nlSjiKS2E11dXYFtznePfIJ1RGF0KY2+H9tqDnT6cL3LzNLwYwj6jfZdPq747jN24sBBCm7cazZ4q\n7O/yIaM2g7Zde0AL4cZnx2mbjPYDnV7saqtHTaXL8e9shPCBLi9qKt3Y1+F1/BsopRUDjL/X9X2N\nRWfDOT+/jn2dXrhdgldta8R8MObY3/78fAieKm1Gq8Pbm6GU8/nCqIQf6Pbh1Tu08/xJh8A4tBCC\nUsCB7gbsaqtHi6cKJ6btb4bWY0lMLUdxoNuHmko3XrWtsWih4uzcOg7qBYFb+pvx0sVVx+4Tp2eD\n2NPuRW2VG7fsbMZSKO54Tj09G0SV24W9Hd7sDeNRh9BrBOjrehuy18IXJ+yD/vBiCIlUBtf1NqKv\nuRYdvmq86BCSo4kUhhdDuKG3ARVuF27c3oQjRUL4lYyh+hphXKRbPFXoaqxBOqN18bCytqEH8Lrc\nSrV9n+p4TlUb2OwvbfUzCxaV6s25qgsvlObp9AzbWuqQ0GcFyWd0CzEq1YAWopwq1eP+MPpbteBU\nV1WBvuZaxwvS4noc4XgKe9rr0d/qQZXbVTSQ/ui8Fkg//Mbd8NZUZPuJ2RlaCGEiEMFv3LoNN/Q2\nOp6cAS3cPDnsx+2D7bilvwUXV6JYdQhnSik8MezH6/a24cZtjQBQdOT2iYtruKW/GdUVbuzv9uFU\nkZA85g9jT7t2XA/2NODiShTBqH2FfsIfwbaWOlS6XdkLh1NleDwQQV2VG52+GvS31qO20u1YDR83\nPZFwuQSDXV7nEO7fDOEAcLCnAefn121vSJdCcUQS6WxgHOzyIRBOwB+yv9CPByLY1W6E8OKhd9wf\nwU49wO4vISSP+SOoqnChp6kWnb4aNNZVOoYbo5vUzrZ6PSQ7V5KnlqPIKOTs89BCyPZCn84oTAai\n2X0e7PIhFC8eko2bXeMYOYXY4aUQ+prqUFvlxkAJIdnoArS7vR711RXY3lLnGDCHTSFc2wfnG4np\nlShiyUz2KVixkJzJqGxlGwD2d3mLh2RTCB/Q3xdO56ShhRBa66vRUl+N7S0e1Fa6HffBHMLdLtGf\nwti3N363cZNiHCO7m6elkNYdZaBz80YiEI7bfna0EB7K3lgOlrDP5+bWMdDly372AefPzrm5dXir\nK9DbVIs2bzVa66uLtNfOPfu7fRAR/SmM/facnzNCu3auG+j0YWTR/rOzEklgdm0j+wRpsMuHaCJt\nOx2t1j96HQd7tNc3zhdOn50zs0EMdHlRVeHK3mw5tT81E0RdlRu72urRWFelfXYc9tkI4df3NkJE\ncH1vY/a4WTk/v450RuG6Xu0adUNfI0aWwlfttLcM1deI1ehmUDYCrV1IXtVDT2OduVJdg/VYyvIi\nsBJJ5Mz8AWxWqvNnAEmlM/CH4jmDFAHnBWBm1zbQ25gbqrc7zAAymQ3VnuzXdrbWY3Ztw/YiNuaP\nYFfbZvu97V7HvrYj+oC7Xe31qHS7sKejHueLPOa/sKj1R6x0u3Bdb0PRvrM/1h/rGX0kz86tO/bn\nOzu3jnA8hdfvbcf1fdpJ1Cm4G2Hvpu1N6GmsRYunyrEKk0xncHElij3t2sXoup4GnJ0NOs7cMB6I\nZAPpdT3aSdGpWj0eCGertns761HpFsfK8Lg/gv5WLWC69QulU3VrIhDRqrz6nOr79cBotw/mEA4A\nB3t8iCUz2Yp3vrG89gNFLtyrkQRWIonsPpdyoR/zh7PjBNq81WjxVBUJ4WH0t3jgdokWkotUz8f0\nrhbGTen+IiE5+7TAFKrDDiF5ZlXramFu77TPyXQGE4FINsAOdHoh4nyMRhfD2faeUkLyUhg9jbXw\nVFdo29Tpcwxn2cq2EWI7fZhYtg/Jw3ntB4uE5OlVLYTv68y/kbDeh/wQvr25DrWVbsdBtcOLoWyA\ndetPh5z2+fy8Vgk31iMwAqNdSN6sbG8GwOWI/Q2m8fc0bgiKfRYW1mMIbiSz759i7yOlFM4vrGdf\nt6W+Gm3easdjdG5+HYN6QDa2yTFUz6+jsa4S3fr1bbDLiwuLIdubcHMlHNDe29FE2vazY/z993dp\n53fjhuWCzbVnJZJAIBzPdvsy3h/DDteqCwvh7Puirkp7CnPB4WZudCmMPR3e7JPqvR1ex/Zj/jBq\nKl3oa67N7sPkctQ2JI/phQ3jpnpfpzfb/epqxFB9jTBWQ2z2VJnmnbYO1UYfs8a8SjVgPa3ecjiR\nM/OHuf1C3mDFQDiBjNr8vqHTp81VnX/Hnckoy0q1ccG3GgQytRyBt6YCTaabAiPkTFiEoVgyjZnV\naLbSBmgXv/FA2DbEGoHbCJgDnc7VPKWUfvLRfsehnkYMLaw79mEeXQqhqa4S7d5qHOptQDyVcaxU\njerbNNjlxaGeBojAMbgbfY93t2sVSS3o24fwqeUoUhmVPZaHehsRSaQxYTMAdC2aGxiN6opdqE5n\nFCaXo9nXr65wY2+H17HyPBHQQrXhQHcDzs4FHQJgBP16lddoH02kMWnTP388EEG9KYQfyFbPrbfJ\nPOAVQNFK8mZVWNuHlvpqtHurbSs90YTW39lob3S3cO7+sVkJB7TwcWEhZHsjMeaPYFtzHaortO5W\nxoXerho2llfNNy7IdiF2LC+E7+swQrL1Phjvu93mkNxcZ3tMU+kMxgNh7O7Y/DwPdjpXDEcWQ9kq\nMqDdDE0sR2xn6Rle1EJ4vRHCu3xQCrZhYrOPd2k3EkZIMkL4tuY61FW5bffBCOHGsXcZIblICDdC\nmbZNXpxfsK8km0O4sQ8rkUR29d18QwshtOvdUQDTUxib96pxbstWqovcSBjd84zqfLOnCh0++0ry\nUiiOUCyVDZbGPjjdSIwubd6cGe1HFu2vC2P+CHbrT3iM9olUxvK6A5hu8vVJAbJdj2z2eSKQ+1kz\n3h92odo4rxnd0TzV2lPYIZv3aSiWRCAcR3/r5j7v7fDavr6xTf2mp8IDnVr3Urtr22Qggh0tm+fg\nYiF5MhBBhUuyhbd9Rfb5SsdQfY1YzfaprixaqV6zqFQbP2M1A0ggHEdrfW6luqbSjWZPFebyfocR\n5PNDdYXbhe7GwrmqA5E4EqlM9gNl6G6sRaVbLGfomFqO5nxoAdO0ehYB0Hh8nVOp7qhHMq1sB0OO\n+sNorKtEq34zMdjlhT8Ut+3Pt7Ae07qL6CeE63obkEwrxxPDyGIYe9q92UdkgHNIHlkKoarChd6m\nOnhrKrGrrd5xUE1+uLmutxGjS2HbeZU3q7ZG+wbHbcoPW411VehprLWtDM2tbSCRymQvAMBmVwIr\n8VThzdBglw+RRNp2UJq5v7PRHrCv3EwENrtaANrFqabSZdslZdwfQW2lO/t5aayrQldDje0+jOWF\ncACO3S02u6PkBsYLC9bVsERKe7qwszX3GG0k07bvbXN/Z6M9YB8Ax5bC6PTVZKu8+4pUkseWNvts\nA8VDsjENo3EDCzjfxE6tRJFMq5z2g10+TNqE5HRGYdwfyX42jfZK2Y+rGFna7I6itXeuqo4shtDd\nUANvjXZO7Wuqg6fKvrvFqD83hBsh2f4YGWMLcv9udpVkY6zA7rz2a9GkZZc6Y6xAzvu003mfJwLh\n7GcfKH6DORGIoqmuMlvMafJoT1XtbiSMhVj6Tb9joNNnG9qNQNqf8/n3YnTJOiSvRRMIbiQL2ifS\n9iF5ajmCHRbnF7tz3kRAa2+cX4qG5EAE1RWu7PmlvkhInggUdoXc1+GzfX2j6+SOEkPy5qQA5mun\nHpKXbG4kliPZ8U7a9ujVc5t9mFyOoK+5DhVuLY7uaPWg0i2O1fArGUP1NWIlkkR9dQWqK7SwW+V2\nYd62Uq2H6ryBikDhwMN4Ko31WKqgTzUA9DXVFszFmr/wi1lvY11B+9nVwpk/AO1xZW9TnWUwmFrW\n+uWaGaHCql91/uNrwPSYzOaiOroYLqhIAPb9+YYXcy+Sh3qcA6lSSrtw69Wz7S118NVUOFaSR5e0\nwGg8hruutwEnZ+xnzxjzh1FX5c521bmutwEZZX8ByA/JO1s9qKpw2QbG/BAOaCfoC0UqmOaL2L4O\n7WbFauDeRYubIeNxudVFI5pIYS4Yy9kerUpv32dQG0S4+foVbhcGu3y2g1LHA2H0t2r9tc377BSS\nK92Sc9PodKHPvxHS2vsQT2Usq+0XVyJIZ1RupdqhT7I220kk75h64RLgnM17eyyQWwmvq6pAf4vH\nsTrf4qnKeRKm3TyV9r4z2k+tRC1vAMeMgGkOW11erZJssc+zqxt6d5TCY2S1D5mMwmQgkvM3MEKy\nXb9t7eZss302JNsFQH8ErfXV2RBu7LNdn2SrwLi/y4vgRtJyES6jvTncZLuYWPyd/WFtrEBO2Coy\nAHRyOZqzPQ16twi798VkIDeQAs7dLSYCEXiq3NmnSFp7H0aXQpbTiU5a7PNgpw+JdMaySmocI3Ng\ndLrBjCZSWFyP5+zzrjatC5vTUxjzPhsh2S4wTujrL5jPL/s6fLZFgallbdGqvmZTqO6stw3JxjnE\nvE17O+1DsjEpQM452+iSslh4jNIZhemV3H02QrLdOXgyEM1531W6XdjVVs9KNb28ViJxNHm0E7SI\noN1XjUW7SvVGAt7qiuydIbA58HA+rwJo9I8z5rQ2297iKbjQZyvVDYXt+9s8GPNHci4aM3lzVOe+\nfuFc1cl0BjOrGzkfQgCorXKjp7E2G/TMxi0uSLvatLBld/c8shTKduUANqs2dsFgJG9gU29TLZo9\nVbYh2R+OI7iRzLbXumc0OvZHHvXnVs8Odjc4DvTR+pFb3Bg4hOQ2bzV8+oW+wu3CnvZ6+/YBLTD2\nmf52+zq9GPdHLC96m48263Paa9tUeFyNsGX+uxk3Q1YXpUm9amMOZ7VVbuxo8VieoDf0ird5ewDt\nuJ6bW7fsYmIMIjTTLvRhy4vYmD+M7S2enM+acaG3vgGMQCS38rTZb7twH7KB1FSp3tNRD7dLLIPB\n7NoG4qlMzj7XVLrR32odkpVSGF/KrWAa+2wXJMaWIpbt7ULyuD+SUwnX2nttu1tk30emfXYaoGU8\nvTI/8u5t0rp2WIXkxVAMG8l0ToXU5TIGdBa+vlJKr0jmnpOcQvLkcu7NnNF+PZayDMmTyxE01Fbm\ndHnLDla0+OxMWpzzjPeR1U218dkxh6GG2kr0NFqveRDcSGIlksgJsMY22RUeJpcj6G+x/uxYhuTl\nCLbnPZEc7PIimbbuSjCxHEGV24XuRvMNrH1INq5d/aa/2662elS5XY7HyPzZrKrQAqDV6ydSGcys\nRguuVc6V5EjBMd3XWY9xm5A8EYigt6kWlabzy75On21ItrrxMK5tVtfCCYv3Rb9RSV4o/BvMrW0g\nmVY5x9QIyVY3BkopTC4X3mzt6/Q69gu/kjFUXyNWokk0mwYTdvpqbJftXosm0eipzPmat6YS3pqK\ngrmSs6HaZxWq67QqkOmEuBCMocIlaPUUtt/TXo/gRhKB8GZV0mqO6uzrN9fh4nI056I0t7aBVEbl\nVBcMO9s8lgPM8h9fA1rY2tZcZzlYcVlfNGG36fFysb6wo0tadc4Y5CMiONRjP1hxNK/PNqCNKB9e\nDFs+5tf6hW8UPM4F7Pswji2Fc6pzXQ018NVYBwnAGCBXeHKzaz/uD2Ob6bGd0T6VUZbdcMb9Wl/4\nVlP//OwJ3TEMbW6Tt0a70FtdlKzaA1o13PqCUVghBbS+4aF44VR/Vt1RAC1IpGwuYvndUQDnC/2Y\nP4y+prqc6SV3t9ejwiYk5w+cBLSQvNMmJFtVwo19sApn/lAcoXiqsH2nFxdXoghbVZL94ZzKttFe\nKevQO6FX/82KhaFmTxUaTAHTCMlW7a2qvC6XMaDT/n2RHwAHOq37JK9EEliPpXJCO6Ad01AsVdBF\nTvsd0cIQ7tDdIr8bAbB5Q2q9D1HUVrqzTyABwKd/dqz+BlYhHNBCrFNot6o8j/kLbzA3EmnMB2MW\n7bXPjtXc85N54ymM9oDN+yIQQV9zbc7Urzvb7GdumghE4RLkVHkr3S7sbq+3PKbZKm/e+2K/zWdn\nelV70pbffrDLi/FApGDgXiajMLUSLdjnAT0kWx4jixDuVACaCETR6atBbdXm+cUIyVZjJDZDeGFI\ntn79wuo/oL1Xrc7Z/lAc0UTa4kbCi7mgNlD1asNQfY1YjSTQbO4j3VBju0LiWjSBxtrC7hw9jbUF\nfVWXspXqwu4c21s8yCjk/MzCegzt3uqcx1cGI0COmJaynl3dgK+mIucxqGFXez1C8dzKTXY6vea6\ngvY7W7Vp9fIvevmPr83bYxW2RkwD/Myc+sIOL4YK2l/X24CRpbDljAHZgU3mwVZdXiRS1jNPjPnD\nUAol9/OMJlKYXdvICUPG9GlWF1WlFMb8kcLA2OnVpsKy6J4xbtleu+jZhd6decGgzVuNxrpK64qk\nP4I2b3XBe8PuBD1uUdkGtMebk8uFF7HsIMK8MJQdrJjXBcSqbz6gPYYHCv8OKX02lV157wvjQm91\nEbOqhFdXuPULvXX7dotjZFdJ3uzjXRgMplc2ChYssuqaYbw+gIKuPmvRBJYjCctKNWB3oY/kVIUB\nLSTbLTJjzAhj5jTryaQ+GLU1b7D1gM3AvWyottjnUCxVcI60qngCpvdF3mw1xmAxq+ocYH3jMRmI\n5gwWA7SQbLd6rlblrcv5rBn7YBnal7XBYvnFjcEuH8b8hZ+dzX22Dsn5xYqpFfsQDhR+dpLpDKZX\nNwpuPJy6pE0GCgOpMXOTVeV5ajmC7sba7IBd8z44Vbatbgys5p63u/E40K2F5Pxz2FxQK1DlB1K7\n+fmVUgVdJwDtGFVXuCynEjXeF2aV+pzVVrMqTSxH0GjqB2/Y3+WzbO/0vpgLxgquIxN2N2clTD16\npWKovsLFU2nHKc0MK5EEmkz9njt92jLiVo8eV6PJnEGKht6m2oKBhEsO3T+MD7O5C8jSehztFv2p\ngc0Aab7j1mb+KAzIgPXJZMrmxAZo3QrC8VROd4jNxTrqC9rv7dD6nuU/eswfyW+w6wtr9I82B2RA\n61edzijLE/rIUgjemoqc4+pUhRldKqxsZwfJ2YQtAAWBblAPpPldG1Yi2qCdwoqk9WP1dEZhyjST\nh8GoelhWJC1CuIjYrko4Hih8RA5o4WPMX/jIeCIQQXdDDeqqKnK+PtCpTW+Wf6G3DeEdXm2qv7zB\nikZf3vz30o4W7SKWHxinV7VHofn7YFcNM/o7W71X7UNyYdcMQAuMs2sbBZWecX8YDbWVBWMk7BaZ\nsatsD3YbA7Ty21uH8N6mWngtFplZjSSwGk0WHCOnRWbyZ4TZ3Aef5SIz44HNaRnz21uGZH2wWFfe\necxuXEX2EXlBtc36RiLbTSlvH7w1lehrLhzoawwWszrn2QZAm2O0v8uLcX/hHMCTAW2civmpE2Bf\nJZ0IaN2UtjUXdnnR9jn3GE3aVP/tPjszq9rCX/nHtMLtwt6OwhtMrcpbWLU19sGy8mz7PrIelD4Z\n0G7y66sr8tpbn7etnpAA9jftm11wco/pjhYPPFVunM3rGhgIJxCOpwreF8a4EKsZjOz2+WB3A87M\nFo7PMWbyKGjf04ClUBxLeU/DJ/R1BfLzgjHGKL97YzaE5/0OY97tYiv0XokYqq9wb/nkU/jEoxeK\ntluJJNCcs4x4DWLJDNY3Ch/PBjeSBXeegHWl2r8egwgsBypus5j2zmrhF4NWUavICTdTy5GcPrlm\nA50+iOSuZDi5HEVNpcsy5Bsh2Nwdwh/WplnKr8wBWnhKZVRBv/CxpTA8pgF+Bru+sIvrhVM5AdqK\nYYD1gisj+jy75gt9tj+fRQVgVJ9b2LrfpnXYMl7TbECfYzj/72xeiTCnvc2jxOxcxHkBsKrChZ2t\nhYNMsoMILU7oA3r/ufygnz/4y9w+ZTFF03jerBaGfTaP1cf92rRp5kehxj7s6/QWDFYct6lgVhiV\nnnmbEG6xTVZhaH5d68tr9VRloNOLhfXcSo92w2j9FGYzAOZtk97FJz9g2i1KYwx2zf9Md+tdiaza\nA4XvOxGxnPZuwqayZWxT/vzZ4XgKS6G4dfsur+UiM3YhfHOf80OyFiTyn7bZvY8mAuGCwWIAsisx\n2u2zZUi2mPVkarlwsFi2fZcPE3ldCYwnJFavP6CvxJh/g6lNm2YdMK32eTIQQXdDbU43JUALgFYr\nMU7YBMYKt0uf9cQmhJd4jBZDMcSSGZsbD2/B2JNsP3jLfba+edKeFlzCMVqOwJc39Sug3WA21FYW\nVEgwlq4AACAASURBVJLtPgubK70Wvj5g/T462OMrGBcS3EhiOZKwbt/bgNVosqCrktVTIe31rW8M\nxvyRgpm5AGwu9JXXfnQpjOoKbWYwszZvNTp9NUVX6L0SMVRfwVYiCUwuR/Gtl2Ydl/jdSKSxkUyj\nub5w3mmrftWr0UTBBx3QBguGYimsmx4BL4XiaPFUF1QwAKCtvhp1Ve6cULoYjBUs/GIQEexpr892\n/wjHUxgPRLJ37vk81RXob/XknHzOza1np6HLd6hXm7v5hGkZWKspygxGZTk/AI4shbJzO5sNdFkH\nTGN/8rt/dPhq0O6txmmLftXmJY0NRgXTqrI9uqQNeCt8VGndh3HMH7EM4UZItguA+SHZWBY5P+hn\nq7wWNysHun04nVf1sHukDmgn6EginROSN+fALmxvVD1OmKYTNAKm1QWgv8UDb00FjudNPzhuE7YA\n68rNmD+MDl9hpQrQ3nunpnMXyjG6l1jd0B3q8cEfiufc3GRnU7GoVBsXMfNiP8v60wWr9ge6rOcM\nH/cXDiIEtP72jXWVBZUh45jmB0yjK1F+MBjzh1Hl1qZ9zKctV567EI/d0wJgc5EZ8/zZjmHLomIY\nT2mDUS0Dps3UgFaDDgFkV2IsrMJG0Zc3WGxzm7wF3XyMfTAWuMrfh/yVGK1m8si27yxcZGZmVRt3\nYh0AC5crV0oVzFJh2J4Nyfk3BoV9wgFkV2K0CuGt9VWW3fyMbjvmz5pxA2t3Y5C/iqnT+2i/Rdej\nQNjoB1/a+0h7GhkqKDoAm+Nt8s+pw4th7LK4jhgrpeZXni8saKs7dlh0tTzQ3VCwiJVxw7wv7zoC\naOevUN5nx2hvnot8s71+vjBdq1YjCSysx7I3DWbaipLA6ZncfTYvQ2/WUFeJbc11BeeX8/PagkZW\n+eJgTwNDNV1exkV2LhhzfHOt6Iu5mCvVRpU1P1RnMkqrVNcWntyMUdPmwYpLobhlVRjQTg7b9MGE\ngHYTEIqnLAcdGva0e7OPErXQAlzXZx2qAe1kYgS6VDqDkzNr2SW383lrKrGnvT5nue+RpcKBXIbd\n7fWoqnAVVJK1JZMLTyTG9En5J0+j6pMfkgGtX3X+Ut/L4TiWI4mCEA7YV55HLWZgMNqnMgrDeSOx\nx5bCOQt8GAY6fXC7pCDojwc2l7o20+bQbsgJsIBpgJzFRelV2xoLAqPVjA2b7ZsAAMfNN0M2lXNA\nu3A21VXipanNv7M/rA2os2rvcgleta0Jxy9utjdCuFV7ADjQU1i50abfK9x+ADi8vQmheCon3Iz7\nIwVTy2Xb72gGABydXMl+LXtjY1F5vqGvEW6X4JhpnzfbF25Tu68Gfc21OGJ6/fVYEkuhuOUNpojg\npm1NODq5mvN1u+4lAHDT9iacnQ3mBEAjhLstxlTcuK0JkUQ6J9zYVXmNfQaQs892j9QBLTxVugXH\nTH/ni3qV1+p96qmuwK62+pz3hVHlzR90aDjY3YDjF1cLAqBV+AO089dEIIJgdLNQMeYPo7uhpuAJ\nCaCFFW3ay2BOe8A6YBor6ZmvDyMO7yNjkRlzoLu4EsVGMm15PtJWYsydYjKVzmB4IZTTFc3MataT\ns/PB7JMBq/b5KzGenQuiXV82vLB9YWXY2P/9FoFuwCIkG+HO6GJoZrXIzFwwhtVoEgd6rK9V+eft\nTEbh3Nx6tkKb72CPNt+2uSvh2TltdUer8UgHe7RFrMzzZ5+d01Z3zH+iarTX9tNUkNL3x+oYDXZp\n1wXz39nYf6uQXK8XvMyVZ+NpgFUIN/Y5P8cMmVbAzHeoR/vsWA2GvpIxVF/BzN0MHj67YNvOeCRs\ntUJi/oqH67EklIJt9w9gc+5oQFth0WrmD8PONk92kNlphxOVYX+3VmWYCESyYdZY+MTKgW4fZtc2\nsBpJ4MJiCNFEGjdub7Jtf0NfI05Mr2VP6McmV9BaX20Z9Ksr3LihtxFHTEEiuKEFj/z+0YBWSd7f\n5cPxqbwQrq+M2GLRReZQTyPG/OGcE0O2f7RFCD/U40MgnBtIU2ltjmKri54RPI5PlxaGaqvcGOzy\n5gQVQAtoxlLX+W7c1oSRpXBO/9zxQMSyby6wGZJfsnhiYFXd2tnqQUNtJV4yhRunJwwiekg2Bf0J\nh/baPjTiwmIo+xTGH44jHE9Zhi1gs3JjXHyzi2NYBBUAOLxdD8lTuftgF9oHOr3wVLlzQuyYP3eJ\ndTNPdQUOdPtyQvLmdHrWv+PVO5pxdHIzAGYXELHZplf3N2M8EMmGm0i8cLCr2c07mpHKqJxQemEh\nZHuMbu7XjtEL4ys57be31FlWefd1euGrqchpf35+HRUusQzVNZVuXN/bmNPeuOhbfZ4B4Jb+ZhyZ\nXM3OuDO0EEIyrWwv9LfubMZcMIbpFe3zGU1oN1J24emW/mYoBbwwsZz92snptWzoyfdq/WbrubHN\n9scvrmU/I/m2Ndehw1eN5037fGJ6FW6XZJe6NnO7BDdtb8Jz48um9s7n4cPbtc+a0cVkZCmMjWQ6\ne+7Jd+O2JqxGk9lwH0umMTQfwqFe6302XueFic19OD0TzC4+le9gTwNckntDeno2iN6m2pxxRYZm\nj7Yolfmcd3o2CBHYhuSD3Q055xfjPHDQ5tp2oNuHkcVQdspIY2YcY5VZq31IpDLZp6TpjMLQfMj2\n2mm8jnmxr3Pz6zhgWmLdbE+HVgAyF4zOza2jtb4qu7qjWU2lG7vb6nNmqzrnEKoBLfSeNq2TcN4h\ntBv7PL2ykV3RWeu3nrC92bqutwFKXX39qhmqr2BjgTAq3YLD25vwzEjAtp0RFBpyFnMxQnXuYAur\n1RQNRpXSHOj8objlRd5weHszZlY3MLu2sXn3b3OiAoA37e8AADx0ZgEnZ4L/f3v3HSfXVd99/HOm\nbe+rlVarsqpWsbos9yIXsI3BTowJDgFDXCihPQ8lJBCSB0LLExJCYgIESCA8mBehk9gYbINs4yJL\ntootCVlWXXl7mdkydec8f9yZ0Za5s7tabdX3/Xr5Je3MsfbMPXfu/d3f/d1zMvM5u0kfWJ96pT0T\npG1e5B5Ub0od0E+kpuJ79lgHFy+tzHrgAdhSX8GLA7Jt6VtkQx9STNu2pJI9DV2DahgHrow4rP8L\nnQPDS1kySdl+R7YMZnoFuWxBdV15AXNL8wadMPqTlqNtvVkzmABbFlWwt6Fr0NR9R11mSIEzQfLA\nA7pTvzy8dg6cYCjf7xkUbLk9RAhOJnnjwvLBmerWnkFL1w61eZGzOmQ6A3gkR+YcnKyqtWc+Qzqr\nlO2OBAzI3KTG7VRHmFAkkfVWK8DCygLmlOSxOzVusUSS/aeDricYn9fDpkUVg4Lw3Sc6WVuX/SQJ\nzndtz6muzAOau453UFUUcN1G2+orae+NZYLv9D610eVOTzroTQfuu050Yi1sXpy9/Zb6Cow5Eww1\nhyKc7Ohz/X7OLy9gQUUBO1Ptk0nLrhOdbHFp7/UYti2pHBSQ7jreydq6smG1vGmXLK1i/+lgJrjZ\ndbyT4jyf64n7kqVV9EQTmTKW9IXdFpcL90uWVgHwTCoo3dfglPy4baONi8rJ83kyQWx7T5Tj7X2u\niYHKogCra0t5KhVUW2vZc6rLNYA1xnDZsmqefqUtE9zsOdXFqnklWTPhAJctq+Zwcw8t3c5dmL2n\nguT7PYOWcR/o8uVVxBLJzDEmHai5Bb2Xr6gGyJyzDjaGSCQtG1zar19QTkm+j98dcdr3RhMcae1h\nXV32z1ya72fDwnKeOHLmnJgrCAe4ckU1Tx1pzxzz9p8OsqS6KGspl/OZqznW1ptZrOyl00G8HuMa\nYF6+vJpE0mb2i/TFnFtp46Wp/WjH4VbAuWMTjve7tl9ZU0JVUYAnXnbax/uTHGrqdj2+5Pm8bF5U\nweMD4oYDqdIMt+PLtiWV7DzWkTm3HWgMud4tAOdc1RSKZO4qpoPqVS59Sn/P0xeMZ9pnP6ZuWFiO\nMWe+azOFgupp7Giq6H/F3GLXZZkBQuHhQXXA56GqKDCs/CNdKpLtir66KI+Az5PJVPcnLW09sZyZ\n6vRJ5tmj7exr6KK+qjBrRiWtrryADQvKeOjFRvae6sp5IAS4eEkVdeUFPLDzJC+c7KS62D2IgDNZ\nj+eOd9DQGaYxGOHiVLCQzUX1FSSSNpOt+d0r7XjMmczj8PaVxBLJzBW928wfadlWVjzi8iAkOBnM\n4jzfoIzkkRxBuDGGrYsrB2U803OHu2UkNy+uoC/Wn7nDkG2p64E2LHRq1Ydmkt3a+70e1i8YHCQf\naenJWk+d6dOiCg63dGemdNt/Opgqt8l+iEoHbukM/TNH3e9IgLNfGENmO/3uSBt+r3ENhvL9XpbN\nKco8HPTEEedkdtny6qztnXE4EyS/cLKTcLzftT3A1voKDjWFCEXidPbGONAY4vJl7u0vqq8gEk/y\n4qtOdujpo+1csrTKPQgfcoH2zNEOls0pyjo9JjjZuXy/JxP0Pnu0HV8qs5lNab6fNbVnsufp4Pri\nJVWun2Hbkkp2Hu/IfG+6+uKZYD6bi5dUcby9j+ZQhGiinz0NXVyU407VxUsr6U8F6+BcqGxaVJ71\nDky6PZw5ce8+0cncUvf9aHlNMdXFgUHtATYtzN6nPJ+XrfUVmUAifZzZ5BIkA1y2rIrdJzozc9O3\n9UTZ5HIhBHDpsiraemK83NJDMmnZdyroGoSDEyQDA/rUybq6sqx1rQDbllTh85hM0LvnVJDSfF/W\nGm9wjvFLqosy7dPHvvUumXCvx3DZsiqeeNm5MHjp1ZBTFpgzSJ7D3lNdBPviBPvinOzoc83+p9t3\nRxOZZxJePB3MHJuzuWrlHAAeTwWxL74aYvmcYteLua31FRT4vZkg+cXTIfxe43peqCnNZ01taaZ9\n+qLOLVPt8RiuXjmHHYdbnaRJaoEttyAcYPuqGg42hmgKRoglkrzc3JMpF8rm2lU1hOP9me/xwcZu\n14uIdHuAxw62AM6DnXNL81yTZFsWV1BW4OeRdPtUGdhqlwveyqIAmxaW82iqfV8sMWzKz+nI5HoA\nbrraunWr3bVr1+T/4oc+Bk37J+3X7WnoosDvpSjPS0NnmG31lXiynEBbup2rxU0LywfV0O473UXA\n6xmUpenoi3G4uZt1dWUUZcka7m3oIt/v5YK5JcT6kzx/spP6qiLXGT0szgmssjBAMBynJN/nWmuX\ndrornHmAYklVUSar7qahq4+GzjAe45StrMzx71ucANnvdWYIOdrWy/q6sqwZUoBEMsmuE50sqChg\nQXmhE7AA61wOVvGkk7FZWFFIXXkBffEE+xqCLK12D1aeP9k5aLscaAzRb63r7zjYFCKeSGZOQunt\ndVF9Jd4s498YDHMilSEMeD109jmlMmvnl1KSN/wCJ5LoZ8+prsy4pj/DsjnFrncl9jZ0EfB5MjOg\nPH+yk0WVhcwvyx58nOzoozEYdjLEOMFHXXkBC12mTwyG4xxsCrGippjywgC7jncwrzQ/6yI/AP3W\nsvuEE0gvqS5i94lOygsDLHcpVQAne9SftKxf4DwA4/N4XDM94GTXW3sibF5UwdG2XnoiCTYtKseQ\nPUBrDkU41t7LhfPL6OyLcborzNbFFfg82YOVUCTOgcZQZun5l1t6XMcMyGz3+WUF1JTksaehK+f3\nx2J5/kQXxfk+Vs4tZtfxTqqKA64XQ+Dsm7FEkg0LyzInerfSBnBmIGjpdrbRqY4+2npibK2vcN1G\n6WPVuroyuiMJjrf3snFhOfm+7MFKTzTBi686+2a+35l7d+XckkHPjwzUb62z75TlU1de4Hy3ywuy\nPjiZtqehi7zUvv38yU5ne+U4xhxu6c7sC4ebuwnHk2zMUcJ2uquPU51htiyqoDEUobErzFaX7zKQ\n+f6uri0l3p/kSEuP6/EanIcxXzjVxeJUQmOk49HAY/aSOUU8d7yDeSXu3zVwapyT1kkS7Dvdhd/j\nyRlwHWvroTW1Lxxt7SXYF2PzYvf9oikU4Xh7LxsWlNPVFxt0PMumOxLnpUbneOHzOHO+r5pXknX9\nBThznK8rL2BeaT67Rzh+WSwvnOyiKM/HBXNL2H2yk7ICf87jy6GmEOF4P5sWVjjH8P4k612y7eAc\nI18NOseI011hmkIRLqqvxOOyjdp6ohxpdY4Rkbiz9Pr6BWUU+rPvF32xBPtOO/tCUZ6P/aeDLJ9T\n7Jp5TlrLrhMd1JTks6iykOeOd1BbVjBs2sSB9jV04fM6Gfx9DUHyfB7Xu0LglEoGw3G2LHZKCnsi\niZx3ntPnvs2LKmjrifJo51xe++F/d923J5IxZre1dutI7ZSpnqYslmi8nwK/N3NgiWVZaQ8gkXoi\neGg2JuD1DPt/0g9GuGUACwPezK3TM22zf8kBDIbSfD+dfTFi/clBqxa6mVOSR2UqAMqVBU+rKcnH\nYwzFee7ZkYH9qS0roCc1tZbPY1xvgwL4PB6KAj7aemIkkkl6oomsD3Gm+T0eCvzeTMnNmXIa9xKW\n4jxfZpsmraUnmqA4R59K8nz0xftJJJ3t3xtNEPB6XE/C6Sfq031K/64Cl6xKns+D3+vJ3OFIT7vo\ndisUnLsgoXCcRDKZqYnLdUeisiiAJTX/dWobVeTYRqUFPgJeD63dUULhOJbc29RrDJVFebSnnuJP\nJG3OcQNnSsdwvJ/O3jh9sf6c/QdnFdGkJdOnsgK/a1AAUFUcwGMMTaEIoXCcooDPNaAGKMn3UeD3\nZtp7jMn5/Ql4PVQWBmjujtCZGoPSHJ/BYJhTmkdnX4z2nhj91uZsD842iiT66eiN0RtNZJasdzOn\nxNlGzaEooUiCknxfzm1UUeRso8ZghO5IHL/XQ57PfRsV5XnJ93lpDIYz+2lJjm3kNYbywgAt3VG6\nUvt3thknBqoqchICXWHnGJbr3wdn3471J2nvidEdSYzYPr3fN4UidIfjFAZ8rt9lp7/O+03B9H5B\nzmNYns9Lgd9Lcyia+a6V5Lv3yWAoL/DT0RujszeGtVCUoz045YK9sQSdfTHCsf4Rj/NlhQGS1ma+\nO0V5ufeLdDliS7fzQGDA63ENqMHpr9djaO2OZu6+5jp++TweivN8tPfEMvNP52pvMJQX+lPT0EWJ\n9ydztnc+Q4BoIkl7b5TuSHwU7Z3P3NbjfD+LAz7XgBrIrCDa3hOjpTuC32Ncj/Hg7DMBr4f23hgt\noQiG3PuFxzjn8o7eGM2hCJbcx3iA8qIAoUiC9t4Y4Xh/zmM2ON//RNLS3uPse7nOCU779GeO0hSM\nUF0cmJKAeixGjoDkjJs+P2m/6mR7L3/0f3/L3129nvllBfzJN5/l+9ddkim3GOh7D/+er/z2CEfe\ncTMMCKy/85P9/PLFJp5/xw2Z177/yMv84yOHefkdN0GWg9ZjO17hcw8d4vk33cBvDrXwoWN7+dUd\nV1HlUksKsHfXKT76w30A/Oz2y5mf49YjQABYOdIGGNJ+VSTuHHRcbuMOVBbr571/9xiJpOWLb9qA\nWT03Z/sT+17lvd97gdcG5vJwrJnv33YJC7Ns57Rf/vIQX9vxCjveuZ0P/ddeeqoSPHjPla7tdzx+\nlM88eJBHb7ma42293H14F99+yzaWpG4xDvXqsQ7e/LWn+eLGDbxufS23ffrX3Lapjs1/sC5r+7z+\nJO/93KNsLq3g62/byn3/9ASFtV5+ePdlWdsb4H9+cYDvPnOCZ993HR9+4AUaE2Eeve8a18/Q+2qQ\nP/ryk3x6zVp2HG7jYF+IJ+/bDi7BQZG1fOIfH6fM7yyN/FSwjZ33XT9o/xzap0cfPsS//vYVttfW\n8DTtvHDfDeCSwQToONbBm772NEt7izgW72XXO66HHPX/RdEE93z2UXwdhq5YnF/ccQV1OW4xFwFf\n+MrvONgYIhJP8qVbN7JsU51rex/wi5++yPd2niRpLe++ehnrblzl2t4A+585wSd++iKBsIcrV1Sz\n7e0XubYHCB7v4I6vPk2g1UNVcYCn7r3WdQwA8kMR3vaF32CanYvvp++6FnKclMr6k/zvv/sNHc0x\nIvEk333DxSxa4V6SUgTc/81n2Xmsg2giycdvXs2qq5a6tvcDv/rvA3zrd8fwegyvXz+fLX+00bW9\nAV7edYqP/HAfgYiH5XOKc37XAGINQe74lycJNDoB+zN3Xwc5ApySvhj3fOE3RBr6SVrLL++8itoc\nx7uKpOWv/+kJjjT20J+0fPWNW1h24TzX9oXAd7/3PA/ubyRp4SOvvYB125e7tvcBzz72Mn//q8MY\nA7dtrGNbjm0E0HKgmXu/swt6YcOCMn5yz+Wu3zUA09LDnV96nESDpbYsn1/+6VWQ5TmbtNJInPd9\ncQetp6IU+L384q1XgMszGwDlScvn//WpTLnLd+7cxiqX4x1APvDD/9rLD3c3APDXr1/D5suXuLb3\nAC88/gqfffAQAG+5eBGXuBwf01oPtfCO/3gO+pwHTh+495Kc351kSzdv/fKTRBuS1JUX8Mi7r4Yc\nFzf5fTE++KUnaGqIkO/38Mg9V0OOOyRFQ7bRA/deAsvczzt+4L9+sJcfPe9so0/dupYtl9a7tjfA\nzh2v8PmHnG30J5cs4uLbcm+jhgPN3POdXRB2ys1+8M5Lc26jaHM3d/2zs42qi/N48j3bIUegnxeJ\n896/30Hb6Shej+G3d18DOTLhBdby2ft/l5qiFf79zbmPj9OBMtXTVGcq41BdHKC2PP3Q4fA5p8HJ\nUJYW+IcFnPNK8+nojQ2aw7i1J0JFod81U52uS3vp1SB7G7ooCnhdn/5Pe9PWhez55A288Fc3ZBY8\nOddK84d/PjcFAS8/f+8VPPaha7huhIAa4KYLa1lSXcTDLzWzfkFZzttRAG+/rB6vx/C5hw6y+0Qn\n21e5nywAbttUR8Dr4dtPHefXB5opzvNxydLcdd5Lq4v43s6TPPlyG32xfl671v2k7fd6uH3zAh47\n1MLuE50caAzlbA9wx9YFxPqT/OczJ3jmaDs3rMndfu38MlbXlvKdp0/w5JFWrl9d41rLC06N8Ru3\nLGD3iU5++WITV6+sGXH87tiyEAs8eqiFK1dUD5sOcKiL6itYV1fG8fZeXrtmHlU5AmpwZtB419VL\nqS7O485ti3LOUpP2wetXsnZ+Ge+5Zhk3r6sdsf29Vy5ly6IK7r1yKe++ZtmI7f9wcx2XL6/i9i0L\n+NRtF47YfuviCm5ZX8tVK6v5t7dtzTkG4NRu/ukVS1g5t4R/e9uWEbM8fq+HD96wkvKCAJ++dW2m\n/jaXD1y3goDPw7uuXsbbLls8Yvt7r1pKeWGA16ydxydfv2bE9rdtquOCuSWsnV/K1966ZcT26xaU\nccv6WuaU5PHAfZeMmFUtLwzwnu3LyPN5+MZdW7NOjTmQ12P42E2r8Bgn+LsxR0Cd9uHXXECB38tb\nLl7Eu68eeb+4+4qlLKos5OIllXzuD3MHQgDXr67hmgvmsLCygH9729YRv2vLa4q596qlFAa8fP2t\nWzNZUDcl+X7+6pY1eD2Gz9++LutD0wN5PIbP/MGFBHwe3nPNskyNci4fv3k1NSV5XHPBHO7KESym\n3X3FUrYsrmBxVSF/efPqEdtvX1XDndsWURjw8rk/XD/id2d5TQmfeN1qPMYJYHPdLQBnP/rSmzfi\n8xjed+2KnCVH4OxH979lM5VFAV63rpZLcwTUaX9724WsqS1leU0xd25bNGL7d161lFvW11KS5+N9\n164Ysf31a+bykddegNdj+NBrLhhxG62YW8K//PHm1Gde7lpznlaS7+c/795GRaGf2zfXZZ1KcyBj\nDN+86yJW1pSwrq6Ma0axH0011VRPUzsOt3LXt3byo3dfyuraUtZ88mE+euMFvOea4RmOD3z/Bfac\n6mLHR7YPev0Hz53ioz/axxMf3Z7Zed/1n7s52tbDr/7X1Vl/b1dfjI2f+jUfu2kVD73YRIHfw/fv\nu/Tcf8Bp5sCrIY619XLjhfNcH2oa6KM/3MsPdjWQ5/Pws/denrOODOBDP9jLg/sb8XsNV66cw/1/\nvDln+288cZS//Z+DmaBx9yduIJDjNvmRlh6u/4cdmVUxH//I9syKl25u+ecnMg8F/fg9l414MfHA\nzpP8xY+dZwoeuPeSEU8Cwb44n3voIB29Md5/3YqcDxKl7W8I0tYTZePC8qwP0w6VTFosw0ufZHJZ\na0c8AQ+UTNpRXySDM62k24N02fS7lMTlEk30j3ghN1Ak3j9iEDHe9nk+z6i3a7w/SdLaUX8Gay19\noyjlGKg7Eh+xnGagUCQ+YgnRQD3RBAV+76jHLZZIkkgmXZ+ZGcpaSyicGPEiYqCuvtiIZQ1D25cV\n+Ec9bt2ROAV+76j371giSXyUpZaQ+syRxIilHAMFU+VuozXWcQ7H+gn4PKMe5/6kJd6fHNP351wb\nbU21yj+mqfScwKX5fgoDPsoK/DR2uWSqw9l36PTKho3BSCaobu2Juj6oAM7Vdl15AS+c7OTgqyHe\ncXn9OD/JzLBmfmnOJ6OH+subV3PlijlctqxqxAwpwDuvXsqzx9opyfePKgvzxi0LuP83R9h/Osid\n2xbmDKjByTy9/bJ6fvLCaa5YXj1iQA3w5Tdv4oGdJ50p7XI8aJV257ZF3Lh2HsFw3HWxi4HKCv18\n/vb1I7YbyG0uWzdjCcxk4owloIaxj9tYAmo4u4ussQTUwJhP8BPd3u3uoxszQg1/NmMJqIExBVqQ\nu845m4DPQ2AMN9yNMWMKqCH3sx3nov1Yt2nA5xnxfDCQMWZMATKMXEs91FjHeaSs/1Bej8HrmbqA\neiwUVE9TQ6fJqy3Lp9Gl/MPtqnJellUV23qiORdbAWf6sYdfanLmFp2gco6ZrrwwwOs3zB91+5Vz\nS3jyz68d07+/8+PX0xftz/lwyUB/84a1/M0b1o76dyydU8zHXzfy7feBKooCo8ogi4iInG9UUz1N\nZTLVg4Lq7HNVhyIJSguGB17peZAbUtPXQWoxF5dlx9M+dtOqzFzQCqqnjt/roaxw9LXkIiIiDL0R\nPAAAEwVJREFUMnWUqZ6mQuE4AZ8ncwuwtrxg0AIiA7llqkvy/c5czakV1XqjCfpi/TnLPwAWVhby\n4/dczqHGkOsiCCIiIiJyhjLV01QoMjhQri3Np703Nmh57Exbl5pqgGVzijMr8qXn5xwpUw3OXKy5\nVoMTERERkTMUVE9TQ7PPtamM8dBp9SLxfqKJpOuiDstrinmltQdrbSaori5WTayIiIjIuaSgepoK\nhuOUDnhAbf6AmTwGSq+i5/a07rI5RXRHErT2RGntHn2mWkRERERGT0H1NBUKD55X8sz0eOEh7QY/\n0DjU8hpnIYMjLT28mpqSb7ov8ykiIiIy0yionqaGlX+UOeUfQzPVwXACyJGprnHmE36ltZfdJzuZ\nX5av8g8RERGRc0xB9TQVDMcHZZ8LAl7KC/3umWqXuYznleZTFPDycnM3O491sG1J5ZgXahARERGR\n3DSl3jSUTFq6I8OnyastKxi2quJINdXGGC5eWsUPdzfQF+tn25LcS0uLiIiIyNgpUz0N9cQSJO3w\nQDnbqoodvTEg99Ko77xqKX0xZyq+bUsqz3FvRURERERB9TQU7EuXdGQLqgeXf7R2R/F5DOUumWpw\nAumtiyuYU5LHsjlF577DIiIiIuc5lX9MQ+mSjqEzeswvL6CzL0441k9BwFlpsa0nSnVxXs6lrI0x\nfOUtm+kKx1VPLSIiIjIBlKmehoLh7HXS80qdqfCaQmdKQFq7o6Oad7qmNJ+Vc0vOYS9FREREJE1B\n9TR0Zu7pwTcS5qdWVWzo7Mu81tYT0xR5IiIiIlNMQfUk23OqiwOvhkgmrWubkMvc0/XVhQCcaD8T\nVI82Uy0iIiIiE0c11ZOoobOP2+7/HQDvvmYZf37jqqzt3Mo/5pbkk+fzcKK9F3Cm3kvXVIuIiIjI\n1FGmehKdHJBh3tfQ5douGI7jMVCcN/iax+MxLK4qzGSqg+E4iaRVplpERERkiimonkTpOaY3LSrn\nlZZe13ahiLOaYraZOhZXFWWC6taeKIAy1SIiIiJTTEH1JErPMX3F8mqaQhF6ooms7YLh4asppi2u\nLORERy/JpKW12wmqlakWERERmVoKqifRq8EIFYV+1s4vA+Boa0/WdjmD6uoiIvEkLd1R2noUVIuI\niIhMBwqqJ1FTMEJtWQHLa5xVDY+0ZA+qQ+H4sNUU0+qrnBlAjrf3ZjLVKv8QERERmVoKqifRq11h\nasvyWVRZhNdjeOUsMtX1VU5Afrytl9aeKAGfh9J8TeIiIiIiMpUUVE+ixmCE2vJ8Aj4PiysLXR9W\nDIYTw5YoT5tfXkBFoZ9nj3VwuKmbBeUFWnpcREREZIopxTlJ+mIJguE4tWXOqohL5xRztM2l/CMS\nH7aaYprXY9h+QQ2PHGwmHO/n7ZfVT1SXRURERGSUlKmeJOnp9OaX5wOwoKKA051hrB28smIk3k8s\nkXQt/wC4fs1cQpEE8X7LzetqJ67TIiIiIjIqCqonSWOXE1TPK3Uy1QsqCuiN9WdWT0xzW01xoCtX\nVOP3GuaX5bNxYfkE9VhERERERkvlH5OkpTsVVJedyVQDNHSGKS8MZNqFUkG12+wfACX5ft5/7Qrm\nleWrnlpERERkGlBQPUm6+pxguTyVgV5Q4UyN19AZ5sK6sky70WSqAd533YqJ6KaIiIiInAWVf0yS\nrnQGOhUs15U7merTXeFB7UYbVIuIiIjI9KGgepKEwnFK8n14PU65Rnmhn6KAl4bOvsHtIoODbxER\nERGZ/iY8qDbG3GiM+b0x5ogx5mNZ3n+LMWafMWa/MeYpY8yGie7TVBi6oIsxhrrUDCCD2vUpUy0i\nIiIy00xoUG2M8QL3AzcBa4A7jTFrhjQ7BlxtrV0HfBr4+kT2aapkWyVxQUUhDUOD6nACQKskioiI\niMwgE52p3gYcsdYetdbGgO8Dtw5sYK19ylrbmfrxGWDBBPdpSnT1xSgvHBxU15UXDKupDkXiFAW8\n+LyqzBERERGZKSY6cqsDTg34uSH1mpu7gYcmtEdTJHumuoBgOJ6powbo6I1RURQY+r+LiIiIyDQ2\nbdKhxpjtOEH1n7u8f58xZpcxZldra+vkdu4cCIYTw4LqxVVFAJxsP/OwYnMowrzS/Entm4iIiIiM\nz0QH1aeBhQN+XpB6bRBjzHrgG8Ct1tr2bP+Qtfbr1tqt1tqtc+bMmZDOThRrLaFwnLKCwRno+mpn\nrupjbb2Z15pCEeYqqBYRERGZUSY6qH4OWGGMWWKMCQBvBn4+sIExZhHwY+Ct1trDE9yfKRGO9xPr\nTw7LVC+qdILqE+1nguqWUJSa0rxJ7Z+IiIiIjM+ETjFhrU0YY94LPAx4gW9Za18yxrwr9f5XgU8C\nVcBXUktuJ6y1WyeyX5PNbUGXwoCPuaV5HE+Vf/REE/REE8pUi4iIiMwwEz5vm7X2QeDBIa99dcDf\n7wHumeh+TKV0UD109g9w6qrTmeqWUASAucpUi4iIiMwo0+ZBxdmsK8eCLvVVhZlMdXMoCsDcEmWq\nRURERGYSBdWTwK38A6C+uojW7ii90QTN6Ux1mYJqERERkZlEQfUkyBlUp6bVO97eeyaoVk21iIiI\nyIyioHoShNJBdZaa6lXzSgDY3xCkORSlKOClOE9LlIuIiIjMJIreJkFXXxyPgeLA8M29pLqI6uI8\ndh7rINqfVJZaREREZAZSUD0JguE4pQV+PB4z7D1jDBcvqeTZYx3UlOZpjmoRERGRGUjlH5MgGI5T\nnqWeOm3bkkpOd4V54WQXFy+pmsSeiYiIiMi5oKB6EgTD8awPKaZtW1IJQF15Ae+8eulkdUtERERE\nzhGVf0yCrlT5h5sL5pZw++YFvHHLAgqz1F2LiIiIyPSmCG4ShMJxFlYUuL7v8Ri++KYNk9gjERER\nETmXVP4xCYLheNYlykVERERkdlBQPU6t3VEu/uwjPHaoOev71toRa6pFREREZGZTUD1OP33hNM2h\nKN979mTW93uiCfqTVkG1iIiIyCymoHocrLX86PkGAHYcbiXYFx/WJr1EeXlBYFL7JiIiIiKTR0H1\nOBxs7OZQUzd3bFlAvN/y8EtNw9qkg+pcs3+IiIiIyMymoHoc9p/uAuD9161gflk+Ow63DmuTzl6r\n/ENERERk9lJQPQ6NwQgAc0vzubCujIONoWFtMuUfmv1DREREZNZSUD0OzaEI1cV5BHweVteWcqy9\nl75YYlCbdFCtTLWIiIjI7KWgehwagxHmleUBsLq2FGvh903dg9ooqBYRERGZ/RRUj0NTMMK8Umel\nxDW1pYDz8OJAXeE4Po+hMOCd9P6JiIiIyORQUD0OTaEzmeoFFQUU5/k41DS4rjq9mqIxZiq6KCIi\nIiKTQEH1WYrE++nqi1Nb5mSqPR7Dqnklwx5WDIbjmk5PREREZJZTUH2WmgbM/JG2uraUQ43dWGsz\nrwX7tES5iIiIyGynoPosNYWcoLq2bHBQ3R1N0NAZzrzW0RujolCrKYqIiIjMZgqqz1L2THUJAAcG\nlIA0dPaxoKJgcjsnIiIiIpNKQfVZSmeq5w3IVF8wrwRjyNRVB/vihCIJFlYUTkkfRURERGRyKKg+\nS82hCMV5PorzfJnXCgM+6quKMkH1qc4+ABZWKlMtIiIiMpspqD5LHb0xKouG10qvri3hUGoBmFMd\n6aBamWoRERGR2UxB9VlyC6rX1JZyor2PYDjOSQXVIiIiIucFBdVnyS2ovqi+EoCdxzo41dlHeaGf\n0nxNqSciIiIymymoPkudLkH1xkXl5Ps9PPVKGyc7wnpIUUREROQ84Bu5iQxlraXdJajO83m5qL6S\np19pJ5ZIsio1zZ6IiIiIzF7KVJ+FcLyfaCKZNagGuHRZFYeaujnR0ad6ahEREZHzgILqs9DeEwOg\n0mWlxOtWzcXnMVy2rIq7Lq2fxJ6JiIiIyFRQ+cdZ6OxzguoKl0z1BfNKOPTpG/F5dc0iIiIicj5Q\n1HcW2ntTmWqXoBpQQC0iIiJyHlHkdxY6RxFUi4iIiMj5Q0H1WehQUC0iIiIiAyioPgsdvTF8HkNp\nvkrSRURERERB9Vnp6I1RURTAGDPVXRERERGRaUBB9Vno6I25TqcnIiIiIucfBdVnwclU+6e6GyIi\nIiIyTSioPgst3VHmlORPdTdEREREZJpQUD1G1lpauiPMLcmb6q6IiIiIyDShoHqMuqMJIvEkNaUK\nqkVERETEoaB6jFpCUQBqVP4hIiIiIikKqseopTsCQI3KP0REREQkRUH1GLV2pzLVKv8QERERkRQF\n1WOULv/Q7B8iIiIikqageoxauiPk+z1aolxEREREMhRUj1FLd5SaknwtUS4iIiIiGQqqx6g5FNFD\niiIiIiIyiILqMWrpjuohRREREREZREH1GLWGopqjWkREREQGUVA9Bs2hCN3RBIsqC6e6KyIiIiIy\njSioHoM9p7oA2LCwfIp7IiIiIiLTyYQH1caYG40xvzfGHDHGfCzL+8YY8+XU+/uMMZsnuk9na++p\nLnwew9r5pVPdFRERERGZRiY0qDbGeIH7gZuANcCdxpg1Q5rdBKxI/Xcf8K8T2afx2HOqi9W1peT7\nvVPdFRERERGZRiY6U70NOGKtPWqtjQHfB24d0uZW4DvW8QxQboypneB+jVkyadnXEGSjSj9ERERE\nZIiJDqrrgFMDfm5IvTbWNlPuaFsPPdGE6qlFREREZJgZs9a2MeY+nPIQgB5jzO+noh93fIFqoG0q\nfrdMGo3x+UHjfH7QOJ8fNM6z31SO8eLRNJrooPo0sHDAzwtSr421DdbarwNfP9cdHCtjzC5r7dap\n7odMHI3x+UHjfH7QOJ8fNM6z30wY44ku/3gOWGGMWWKMCQBvBn4+pM3PgbelZgG5BAhaaxsnuF8i\nIiIiIufMhGaqrbUJY8x7gYcBL/Ata+1Lxph3pd7/KvAgcDNwBOgD3jGRfRIREREROdcmvKbaWvsg\nTuA88LWvDvi7Bf5sovtxDk15CYpMOI3x+UHjfH7QOJ8fNM6z37QfY+PEtCIiIiIicra0TLmIiIiI\nyDgpqB6lkZZbl5nJGPMtY0yLMebFAa9VGmN+bYx5OfVnxVT2UcbHGLPQGPMbY8wBY8xLxpgPpF7X\nOM8ixph8Y8xOY8ze1Dj/n9TrGudZyBjjNca8YIz579TPGudZxhhz3Biz3xizxxizK/XatB5nBdWj\nMMrl1mVm+g/gxiGvfQx41Fq7Ang09bPMXAngQ9baNcAlwJ+lvr8a59klClxrrd0AbARuTM0opXGe\nnT4AHBzws8Z5dtpurd04YCq9aT3OCqpHZzTLrcsMZK19HOgY8vKtwLdTf/82cNukdkrOKWtto7X2\n+dTfu3FOxHVonGcV6+hJ/ehP/WfROM86xpgFwOuAbwx4WeN8fpjW46ygenRmxFLqcs7MHTBXehMw\ndyo7I+eOMaYe2AQ8i8Z51kmVBOwBWoBfW2s1zrPTl4CPAskBr2mcZx8LPGKM2Z1aVRum+TjPmGXK\nRaaCtdYaYzRFzixgjCkGfgR80FobMsZk3tM4zw7W2n5gozGmHPiJMebCIe9rnGc4Y8wtQIu1drcx\n5ppsbTTOs8YV1trTxpga4NfGmEMD35yO46xM9eiMail1mTWajTG1AKk/W6a4PzJOxhg/TkD9/6y1\nP069rHGepay1XcBvcJ6X0DjPLpcDbzDGHMcpxbzWGPNdNM6zjrX2dOrPFuAnOKW403qcFVSPzmiW\nW5fZ4+fAXam/3wX8bAr7IuNknJT0N4GD1tp/GPCWxnkWMcbMSWWoMcYUADcAh9A4zyrW2r+w1i6w\n1tbjnIsfs9b+CRrnWcUYU2SMKUn/HXgN8CLTfJy1+MsoGWNuxqnjSi+3/pkp7pKcA8aYB4BrgGqg\nGfhr4KfAD4BFwAngTdbaoQ8zygxhjLkCeALYz5kazL/EqavWOM8Sxpj1OA8ueXESRj+w1n7KGFOF\nxnlWSpV/fNhae4vGeXYxxizFyU6DU6r8PWvtZ6b7OCuoFhEREREZJ5V/iIiIiIiMk4JqEREREZFx\nUlAtIiIiIjJOCqpFRERERMZJQbWIiIiIyDhpRUURkRkqNb3Uo6kf5wH9QGvq5z5r7WVT0jERkfOQ\nptQTEZkFjDF/A/RYa/9+qvsiInI+UvmHiMgsZIzpSf15jTFmhzHmZ8aYo8aYzxtj3mKM2WmM2W+M\nWZZqN8cY8yNjzHOp/y6f2k8gIjKzKKgWEZn9NgDvAlYDbwVWWmu3Ad8A3pdq80/AP1prLwJuT70n\nIiKjpJpqEZHZ7zlrbSOAMeYV4Fep1/cD21N/vx5YY4xJ/z+lxphia23PpPZURGSGUlAtIjL7RQf8\nPTng5yRnzgMe4BJrbWQyOyYiMluo/ENERMDJXqdLQTDGbJzCvoiIzDgKqkVEBOD9wFZjzD5jzAGc\nGmwRERklTaknIiIiIjJOylSLiIiIiIyTgmoRERERkXFSUC0iIiIiMk4KqkVERERExklBtYiIiIjI\nOCmoFhEREREZJwXVIiIiIiLjpKBaRERERGSc/j/Kek4npQf4ugAAAABJRU5ErkJggg==\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fig, ax = plt.subplots(figsize=(12,6))\n", "\n", "ax.plot(tlist, np.real(p_ex))\n", "ax.plot(tlist, np.real(p_ex_ss))\n", "ax.set_xlabel('Time')\n", "ax.set_ylabel('P_ex')\n", "ax.set_ylim(0,1)\n", "ax.set_title('Excitation probabilty of qubit');" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Software version:" ] }, { "cell_type": "code", "execution_count": 10, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
SoftwareVersion
QuTiP4.3.0.dev0+6e5b1d43
Numpy1.13.1
SciPy0.19.1
matplotlib2.0.2
Cython0.25.2
Number of CPUs2
BLAS InfoINTEL MKL
IPython6.1.0
Python3.6.2 |Anaconda custom (x86_64)| (default, Jul 20 2017, 13:14:59) \n", "[GCC 4.2.1 Compatible Apple LLVM 6.0 (clang-600.0.57)]
OSposix [darwin]
Thu Jul 20 22:25:24 2017 MDT
" ], "text/plain": [ "" ] }, "execution_count": 10, "metadata": {}, "output_type": "execute_result" } ], "source": [ "from qutip.ipynbtools import version_table\n", "\n", "version_table()" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.6.2" } }, "nbformat": 4, "nbformat_minor": 1 }