Skip to content

Instantly share code, notes, and snippets.

@asafpm
Created November 21, 2016 20:23
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save asafpm/a86cd9f33f9ec7e5e360c348c1b0506e to your computer and use it in GitHub Desktop.
Save asafpm/a86cd9f33f9ec7e5e360c348c1b0506e to your computer and use it in GitHub Desktop.
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 2,
"metadata": {
"collapsed": false
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Using matplotlib backend: Qt5Agg\n",
"Populating the interactive namespace from numpy and matplotlib\n"
]
}
],
"source": [
"%pylab\n",
"%matplotlib inline\n",
"import qutip\n",
"from scipy.special import erf"
]
},
{
"cell_type": "code",
"execution_count": 67,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"b = [qutip.basis(2,m) for m in range(2)] #Define basis\n",
"tx = b[0]*b[1].dag()+b[1]*b[0].dag()"
]
},
{
"cell_type": "code",
"execution_count": 68,
"metadata": {
"collapsed": false
},
"outputs": [],
"source": [
"def f(t,args):\n",
" a = args['a']\n",
" b = args['b']\n",
" c = args['c']\n",
" return (erf((t-c+a/2)/b)-erf((t-c-a/2)/b))/2"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"# Time-dependent hamiltominan given in list callback format"
]
},
{
"cell_type": "code",
"execution_count": 105,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"(0, 1.2)"
]
},
"execution_count": 105,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAgkAAAFkCAYAAACq4KjhAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAAPYQAAD2EBqD+naQAAIABJREFUeJzt3Xt4XPV97/v3byRZ19HVtiTLNyRzMRc7tkMSSAJpKJCE\nJKT0tNSlZT9lt91J6VPqps1JaNqcXNukG7yTtuykZ+/TQChOs5sL0N2GENJdCoRAsDEJGBtfZUmW\nLEuyNCPJus3v/LG0ZFuakTUza81aM/N5PY8e2aN1+bGQZ33m+7ssY61FREREZL5I0A0QERGRcFJI\nEBERkaQUEkRERCQphQQRERFJSiFBREREklJIEBERkaQUEkRERCQphQQRERFJSiFBREREklJIEBER\nkaQUEkRERCQphQQRERFJSiFBREREklJIEBERkaRKg25AMsaYJuBm4ChwJtjWiIiI5JUKYD3whLV2\nIJsDhTIk4ASEfwi6ESIiInnsDuCRbA4Q1pBwFODhhx9m48aNATclWDt27GDnzp1BNyNwug5n6Vo4\ndB3O0rVw6Do49u3bx2/8xm/A7L00G2ENCWcANm7cyNatW4NuS6Dq6uqK/hqArsO5dC0cug5n6Vo4\ndB0WyLq7XgMXRUREJCmFBBEREUlKIUFERESSUkgIue3btwfdhFDQdThL18Kh63CWroVD18F7xlob\ndBsWMMZsBV566aWXNAhFREQkDbt372bbtm0A26y1u7M5lioJIiIikpRCgoiIiCSlkCAiIiJJKSSI\niIhIUgoJIiIikpRCgoiIiCSlkCAiIiJJKSSIiIhIUgoJIiIikpRCgoiIiCSlkCAiIiJJKSSIiIhI\nUgoJIiIikpRCgoiIiCSlkCAiIiJJKSSIiIhIUgoJIiIikpRCgoiIiCSlkCAiIiJJKSSIiIhIUgoJ\nIiIikpRCgoiIiCSlkCAiIiJJKSSIiIhIUgoJIiIikpRCgoiIiCSlkCAiIiJJKSSIiIhIUgoJIiIi\nkpRCgoiIiCSlkCAiIiJJKSSIiIhIUgoJIiIikpRCgoiIiCSlkCAiIiJJKSSIiIhIUgoJIiIikpRC\ngoiIiCSlkCAiIiJJpR0SjDHvNMY8ZozpNsYkjDEfXMI+7zLGvGSMOWOMOWCM+U+ZNVdERERyJZNK\nQjXwMnA3YC+0sTFmPfDPwFPAZuDLwP8wxtyYwblFREQkR0rT3cFa+33g+wDGGLOEXT4CHLbWfmz2\n7/uNMe8AdgBPpnt+ERERyY1cjEl4G/DDea89AVyTg3OLiIhIhtKuJGSgBeib91ofUGuMKbfWTuSg\nDZJjiQQ89xzs3QsDAzA87LwGYK3ztdifUx6XaabNOMtsNOU2S6pvebhfWM/56rK/p3HmMlpnzubx\nSARqaqC2Fq66Cq67DiorM2+DiBS2XISEZNy3t0VvCTt27KCuru6817Zv38727dv9apdkKZGAT38a\n/uf/hO5uKCuDxkaoq4OSkrM3NmPO/5r/WipHt97F0JpvUHqmmeqht7HupYcpmam5YLiYL93tM9kn\nyDYNbPgyvVv+kMh0Nev+/SmqBt8KwMwMxONw+rTzvbwcbrvN+f+lsCCSf3bt2sWuXbvOe214eNiz\n4+ciJPQCzfNeWwmMWGsnF9tx586dbN261beGifc+9zn47GfhIx+B7dvh2mudT69eeP3U61z+tw9z\n99V301jZyBf+4wv87l0Pcvdb7vbmBAXimz//Jtu//Yf84Vv/kBd7XuS1mvfy2G89zZUrr5zbxlrY\ntw/+5V/gz/8cRkbgO9+BZcsCbLiIpC3ZB+fdu3ezbds2T46fi5DwY+C98167afZ1KSCPPQaf+pQT\nEj75Se+P/7mnP0dbbRv33XQf5aXl7Du1j79+4a/5yNUfIWK05AfA4aHD3PndO7lz853cd/N9jEyM\n8K6vv4vb/vE29v/+ftyxxsbA5Zc7X5s2wQc+AHfcAd/8plPxERGBzNZJqDbGbDbGvGn2pfbZv6+Z\n/flfGGMePGeXrwIdxpgvGmMuNcb8HvB/Afdn3XoJjSNH4Dd+A37pl+Dee70//oGBA+z6+S4+8Y5P\nUF5aDsA9b72H/QP7+cGhH3h/wjz17de+TWmklAfe9wARE6G+op77brqPNwbf4OXel5Puc9NN8I//\nCP/0T/D44zlusIiEWiYfv94M7AFewhlTcB+wG/j07M9bgDXuxtbao8AtwC/irK+wA/jP1tr5Mx4k\njz3yiPP9wQe961441+f/4/O01LRw15a75l57+5q3s6VlC1/5yVe8P2GeenT/o9zYcSPVy6rnXrtu\n3XXUV9Tzvde/l3K/D30IrrzS6XIQEXGl/XZurf13a23EWlsy7+uu2Z//lrX23Un22WatrbTWXmyt\n/YZX/wESDt//PtxwA0RTTzrI2NjUGLt+tosdb9tBRWnF3OvGGO556z3868F/5cDAAe9PnGdOjp7k\nuePPceult573ellJGbdcfAvf2586JIAzgPHxx2Fy0ZFCIlJM1JErWRsehh//GG6+2Z/jP9/1PFOJ\nKd6z4T0Lfnb7lbfTUNHAIz97xJ+T55F/PvDPALz/kvcv+Nmtl97KK32vcGToSMr9b7vNmfXwf/6P\nXy0UkXyjkCBZe+opZ2qdXyHh34/+O02VTVy+4vIFP6soreAda9/Bc8ef8+fkeeTR/Y9y7ZprWVm9\ncsHP3rPhPSwrWcaj+x9Nuf+mTdDeri4HETlLIUGy9sQTcMklcNFF/hz/6c6neee6d6acwXDtmmt5\nvut5ZhIz/jQgD4xNjfHkoSf54KXJn7cWLY/yi+2/uOi4BGOcasL3vueEPhERhQTJirXOeIT3LOwJ\n8MTE9ATPdz3PdWuvS7nNtWuuJTYZ4+cnf+5PI/LAk4eeZHx6fMF4hHPdeumt/Efnf3Bq7FTKbW67\nDfr6nO4jERGFBMnK/v3Q2elfSPhpz085M32G69alDglvXvVmSiOlRd3l8NSRp2hvaOfS5Zem3Oa9\nG95LwiYWvU5vfSu0tGgqpIg4FBIkK9//vrO07/XX+3P8p489TXRZlM0tm1NuU1VWxZaWLTx7/Fl/\nGpEH9vTu4epVVy+6zera1SyvWs6eE3tSbhOJwJYtzmqMIiIKCZKVZ5+Ft70Nqqr8Of7TnU/zjrXv\noDSy+OKgb1/z9qKtJCRsgr29e3lTy5sW3c4Yw5aWLezpTR0SADo64NAhL1soIvlKIUGycuSIM2jR\nD9OJaZ7tfHbRrgbXtWuu5cjpI5yInfCnMSF2eOgwsckYW1q2XHDbpYaEw4fPPrVTRIqXQoJk5ehR\nWL/en2Pv7d1LbDK2pJBwzRrnccg/7iq+EXdu98GW1iWEhNYtdA53MjA2kHKbjg44cwZOFF/eEpF5\nFBIkY/E4DAzAunX+HH9P7x4iJsLW1gs/CXR17WrW1q3l2c7iG5ewp3cPq6Krkq6PMJ9bbUj1HAdw\nQgKoy0FEFBIkC8eOOd/9qiS81v8aHQ0d5y3FvJhr11zL893P+9OYEHu59+ULjkdwXdx0MdVl1Yt2\nObS3O2smKCSIiEKCZOzoUee7nyEh2SqLqWxauYlXT76KtdafBoXUnt49SxqPABAxETa3bF40JFRU\nQFubQoKIKCRIFo4ehbIyaG315/iv9r+aVki4bPllDE8M0zfa50+DQqg33ktvvHfJIQFmBy8uMg0S\nNMNBRBwKCZKxo0dh7Vp/Hg09MjFC10gXV6y4Ysn7bFyxEYB9/cUzyT+dQYuuLS1b2D+wn7GpsZTb\ndHTAwYNZN09E8pxCgmTs2DH/uhrcG306lYSOhg5KI6XsO1U8IeHl3pepLa9lff36Je+zpXULCZvg\nlb5XUm6jSoKIgEKCZMHP6Y+v9r+KwSy6zPB8ZSVlbGjcwOunXvenUSG0p3cPb2p5U8qHXyVzxYor\nKI2ULtrl0NEBQ0POl4gUL4UEydjRo/5Nf3yt/zXaG9qpKktvKceNyzcWVSVhb99e3tS8tJkNrvLS\ncjYu33jBSgKomiBS7BQSJCNjY9DfH56ZDa7Lll9WNJWE6cQ0h4cOc9nyy9Le95KmSzg4lHrQgUKC\niIBCgmTI7zUS0p3Z4Nq4fCNdI13EJmI+tCpcOoc7mU5Ms6FxQ9r7Xtx4MW8MvJHy5w0NzpdCgkhx\nU0iQjLhrJPjR3RCbiNE53JlZSJid4VAM1QT3Jp9JSNjQuIHO4U4mpidSbqPBiyKikCAZOXoUSkth\n1Srvj+3e4NOZ/ui6tOnS845RyA4OHqQsUsaaujVp73tx08VYLIeHDqfcRiFBRBQSJCPHjsGaNU5Q\n8Nqr/a8CZNTXHi2Psrp2dVEMXjw4eJD2hvYLPkY7Gbf6cHBw8XEJCgkixU0hQTLi5/TH1/pfY339\neqqXVWe0/8blG4ujkjB0MKOuBoDWmlaqyqoWDQktLXDyJBTZKtcicg6FBMmIn9MfDwwcmOs2yMRl\nyy8rikrCGwNvZBwSjDFsaNzAG4OpBy82NcHkJIyOZtpCEcl3CgmSET9XWzxy+ggX1V+U8f4bl2/k\n4OBBpmamPGxVuMwkZjg8dDjjkABOl8NilYSmJuf7wEDGpxCRPKeQIGmbnobeXli92vtjW2s5MnSE\nixoyDwmXLb+M6cQ0R04f8bBl4XJ85DhTiSkubrw442Nc3HjxBSsJoJAgUswUEiRtp087392biJcG\nxweJTcayqiS4AePo6aMetSp8spn+6LrQNMjly53vCgkixUshQdI2OOh8b2z0/tjup//2hvaMj7G6\ndjUlpqSgQ8LBwYOURkpZV5/5wJCLGy8mYRMpKy6qJIiIQoKkzdeQMOTcsLLpbiiNlLK6dvXcsQrR\nwcGDrK9fn9H0R9eFpkHW1EBZmUKCSDFTSJC0+V1JqC2vpaGiIavjrK9fz9Hho940KoQODh3MajwC\nwKroKipLK1Muz2yMU01QSBApXgoJkjY3JDRkdx9P6siQM7PBGJPVcdbXry/o7oZspj+63GmQF5rh\noJAgUrwUEiRtg4NQWel8ee3I6exmNrguqr+oYEPCTGKGQ0OHsg4J4CzPfKEZDqdOZX0aEclTCgmS\ntsFBf6oIAIeHDmc1s8G1vn49vfFexqfGPWhVuHTHupmcmfQkJGxoUCVBRFJTSJC0DQ76Mx4hYRMc\nGz7mWUgAODZ8LOtjhY07IDObGSCu9fXrOT5ynJnETNKfL1+ukCBSzBQSJG1+hYSeWA+TM5OedDe4\nIaEQuxw6hzsBWFOb/tMf51tXv47pxDS98d6kP1clQaS4KSRI2vwKCXPTHz2oJLTVtlEaKS3IkHB8\n5DhNlU0ZPwDrXGvr1gJng8d8CgkixU0hQdLmW0iYXdTHrQJkozRSypraNQUZEjqHO+du7tlyj5Oq\nW6apCUZGYKpwH4MhIotQSJC0DQ35V0lorm725BMyFO40SC9DQm15LXXldYtWEuDstFcRKS4KCZI2\nPysJXoxHcK2vX1+QD3nyMiSAU024UEhQl4NIcVJIkLQkEj6HBA/GI7gKsZJgreXY8DFPBi26FgsJ\n7kOetFaCSHHKKCQYY+42xhwxxowbY543xlx9ge3/0BjzujFmzBjTaYy53xhTnlmTJUixmBMU/Opu\n8DIkXFR/ESdHTzI2NebZMYM2PDFMfDKuSoKI5ETaIcEYcztwH/ApYAuwF3jCGLM8xfa/DvzF7PaX\nAXcBtwOfz7DNEiC/ntswnZimO9ad1VMN55tbK+F04ayVcHz4OICnIWFd3bqUIaGhwXmGg0KCSHHK\npJKwA/iatfYha+3rwIeBMZybfzLXAM9Ya//RWttprf0hsAt4S0YtlkD5FRJOxE6QsAlPy+huSCik\ncQnuzdzrSsLQmSFiE7EFPyspgfp6hQSRYpVWSDDGlAHbgKfc16y1FvghThhI5jlgm9slYYxpB94H\n/O9MGizB8uvhTl0jXQCsrl3t2TFXRVdRFikrqEdGdw53UhoppaWmxbNjaq0EEUkl3UrCcqAE6Jv3\neh+Q9F3LWrsLp6vhGWPMJPAG8G/W2i+meW4JAb8qCX6EhJJICW21bXPHLgSdw52srl1NSaTEs2Mq\nJIhIKl7NbjCATfoDY94F3IvTLbEFuA14vzHmkx6dW3JocNApQUej3h63a6SLqrIq6ivqPT3u6trV\ndMUKKCSMeDv9EaA12kqJKVFIEJEFStPc/hQwAzTPe30lC6sLrs8AD1lr/372768aY2qArwGfW+xk\nO3bsoK6u7rzXtm/fzvbt29NstnjFnf5ojLfH7RrpYnXtaozHB26LttE90u3pMYN0fPi45yGhNFJK\nW23botMgDx/29JQi4pFdu3axa9eu814bHh727PhphQRr7ZQx5iXgBuAxAOO8q98AfCXFblVAYt5r\nidldzeyYhqR27tzJ1q1b02mi+MyvNRK6Yl2edjW4VteuZveJ3Z4fNyidw528Y+07PD/uurp1dI6k\nriS88ILnpxQRDyT74Lx79262bdvmyfHTrSQA3A88OBsWXsCZ7VAFfB3AGPMQ0GWtvXd2+8eBHcaY\nl4GfABfjVBceXSwgSDj5FhJGutjQuMHz47ZF2+iOdWOt9bxKkWsziRm6Rro8ryTAhddKUHeDSHFK\nOyRYa781uybCZ3C6HV4GbrbW9s9ushqYPmeXz+JUDj4LtAH9OFUIjUnIQ349t6FrpIt3rXuX58dd\nXbuasakxTp85TUOlx1MycuxE/AQzdsa3kPBM5zNJf9bU5IRDa73vZhKRcMukkoC19gHggRQ/e/e8\nv7sB4bOZnEvCZXAQ1q/39pgziRl6Yj2sqfNujQRXW20bAN2x7rwPCX6skeBaW7eWrpEuZhIzC2ZO\nNDXBzAwMDztrJohI8dCzGyQtfnQ3nBw9yXRi2rcxCUBBTIP0Y7VF19q6tczYGU7ETyz4mZZmFile\nCgmSFj9Cgh9rJLhaa1oxmIKY4dA53ElteS215bWeH9sNHsmWsHYf8qSQIFJ8FBJkyaz1JyQcH3E+\nIfsREspKymiuaS6ISoI7TdQPbdGz3TLzuV0Mp0/7cmoRCTGFBFmy8XGYmPCnklBeUk5TZZO3B561\nunZ10ptfvumJ98zdzL1WX1FPZWll0oqLu3BWbOGjHUSkwCkkyJL5+dwGPxZScrVFC2Np5u6R7rmB\nmF4zxtBW25Y0TCkkiBQvhQRZMj+f2+BXGR0KqJIQ62FVzSrfju+uKTFfaSlUVsLIiG+nFpGQUkiQ\nJcvXkFAIlQRrrRMSoj6GhNrUS1hHo6okiBQjhQRZMr+7G/yyunY1g+ODjE+N+3YOvw2MDzCVmPI3\nJKSoJIBCgkixUkiQJXPLzbUezsBL2ATdsW5/Kwm1qUfu5wv3E77vIWHEWcJ6vtpadTeIFCOFBFmy\neByWLXO+vHJq7BSTM5O+VxKAvF4roSfWA+DbwEX32BMzEwyODy74mSoJIsVJIUGWLBY7O9LdK34u\npORypw3m87iEnlgPBkNz9fyntHtnLkylmOGgkCBSfBQSZMn8DAl+zf8HqF5WTX1Ffd6HhJXVKykr\nKfPtHHMLKiWpuKi7QaQ4KSTIkvkREnrjvURMhJXVK7098Dz5Pg2yO9bt63gEgJaaFmcJa1USRGSW\nQoIsmV8hYUXVigVPHvRavk+D9Hv6I5xdwjrVqosKCSLFRyFBlsyvkNBS0+LtQZPI90pCT8y/JZnP\nlWoapLobRIqTQoIsWT6HhNaaVk7EFj4GOV/kopIALLo0syoJIsVHIUGWLK9DQrSV3nhv0jUAwm46\nMU3faF9uQkI0+aqL0agzBTaR8L0JIhIiCgmyZPF4HoeEmlamElMMjA/4fi6v9cX7SNhE7kLCIg95\nisd9b4KIhIhCgixZLAY1Nd4dz1qb00oCkJddDrlYSMnVVtvGqbFTTExPnPe6u8qmuhxEiotCgiyZ\n190NsckY49PjOaskAJyI529IyFUl4dxzuvS4aJHipJAgS+Z1SOiN9wLkJCS458jHSkJ3rJvSSCnL\nq5b7fq5Uz7lw/79rhoNIcVFIkCWZmICpqfwNCeWl5TRWNuZtJaG1ppWI8f+fa6pVF9XdIFKcFBJk\nSdybQ76GBMjfaZC5mv4IUFteS3VZdcpKgkKCSHFRSJAl8SsklJeUU1de591BF9FS05K3lYRcDFoE\nMMY4ayWMKCSIiEKCLJFfIaGlpgVjjHcHXURrtDUvQ0J3rHtu4GUuJJsGWV7uPCJcYxJEiotCgiyJ\nOz/ej5CQK/na3dAb781ZdwNo1UUROUshQZbEvTl4uU5CICEhfiKvVl2cmpni1NipnF6nxVZdVEgQ\nKS4KCbIkfnY35EprtJWxqTFik/lzpzs5ehKA5urmnJ2zLdpGT6xnQZjSQ55Eio9CgixJQYSEmvxb\ndTHXM0DA6W6YmJlYsIS1KgkixUchQZYkFoOKCigt9eZ4M4kZTo6ezHklAfJr1cW+0T4gxyEhxVoJ\nCgkixUchQZbE69UWB8YHmLEzqiRcgFtJWFm9MmfnXGzVRXU3iBQXhQRZknxektkVLY9SXVadV5WE\n3ngvy6uWU1ZSlrNzttS0EDGRpKsuqpIgUlwUEmRJCiEkwOxaCXlWScjloEWA0kgpzdXNSSsJCgki\nxUUhQZbEr5CQ6xugOw0yX+R6cKcr1aqL6m4QKS4KCbIk8bj3IaGuvI7KskrvDroE+bbqYt9oXzAh\nIcmqi+puECk+CgmyJLFYfi+k5Mq3VRcDqyQkCQlud0MerUUlIllSSJAl8aO7IbCQkEeVhCDGJEDq\n7oZEAsbGct4cEQmIQoIsScGEhGgrp8+cZnxqPOfnTtf41DgjEyOBVRIGxgc4M31m7rXaWue7uhxE\niodCgixJwYSE2bUS3IGTYRbEQkoud62EnljP3Gt6XLRI8VFIkCUplJDgntO9AYdZUNNEIfmqi+7/\nf81wECkeCglyQdZ6GxImpicYOjMUyM2vucbp38+HSkKgISHJqouqJIgUn4xCgjHmbmPMEWPMuDHm\neWPM1RfYvs4Y87fGmJ7ZfV43xrwnsyZLrp05AzMz3oWEIMvoTZVNlJgS+uL5UUkoMSU0VTXl/Ny1\n5bXULKs5r5KgMQkixSftx/UYY24H7gN+F3gB2AE8YYy5xFp7Ksn2ZcAPgV7gNqAHWAeczqLdkkPx\nuPPdq5AQ5CfkkkgJK6tX5k0lYWX1SiImmILf/GmQ6m4QKT6ZPNNvB/A1a+1DAMaYDwO3AHcBX0qy\n/X8G6oG3WWtnZl/rzOC8EhCvHxMdZEgAp8shH0JCXzyYhZRcbbXnh4TKSohEVEkQKSZpfUSZrQps\nA55yX7PWWpxKwTUpdvsA8GPgAWNMrzHmZ8aYTxgT0McjSZt7U/BqMaXeeC8RE2FF1QpvDpimlpqW\n/Bi4OBrM4E5XW/T8tRKM0aqLIsUm3Rv1cqAEmP8O2wekejdrB35l9lzvBT4LfBS4N81zS0D8qCSs\nqFpBSaTEmwOmqbk6PyoJQc0AcaVadVHdDSLFI5PuhmQMkGqx1ghOiPjd2arDHmNMG/DHwOcWO+iO\nHTuoq6s777Xt27ezffv27FssS+ZHSAjy5tdS08LTx54O7PxL1Rvv5V3r3hXY+dtq2+iJ9WCtxRgD\n6EmQImGza9cudu3add5rw8PDnh0/3ZBwCpgB5q8Tu5KF1QXXCWByNiC49gEtxphSa+10qpPt3LmT\nrVu3ptlE8VohhoS+0b7zbn5hY60N/Dq1RduYnJnk1NgpVlQ7XUPR6NmBrCISvGQfnHfv3s22bds8\nOX5a3Q3W2ingJeAG9zXjvMveADyXYrdngQ3zXrsUOLFYQJDw8GNMQpA3v+bqZsamxohPhvduF5uM\ncWb6TOADF+H8tRJqahQSRIpJJoMH7wd+1xhzpzHmMuCrQBXwdQBjzEPGmC+cs/1/B5qMMV82xlxs\njLkF+ATwN9k1XXIlFoOqKijxaAhB0CEhH1ZdDHoGCCRfdVEhQaS4pD0mwVr7LWPMcuAzON0OLwM3\nW2v7ZzdZDUyfs32XMeYmYCewF+ie/XOy6ZISQvG4d10NYSiju+fujfeyoXF+kSscwhASmmuaiZjI\neZWE6mro6VlkJxEpKBkNXLTWPgA8kOJn707y2k+AazM5lwTPyyWZY5MxxqfHA7/5QbiXZnbb5rY1\nCKWRUlpqWuga6Zp7TZUEkeKitQrkgmIxb8cjQLCfkBsqGiiLlIV6aea+eB/lJeXUldddeGMfzV8r\nQSFBpLgoJMgFeVlJCENIMMbQUtMS+kpCS01L4LMv5q+6qJAgUlwUEuSC4vHCqiRA+JdmDnrchmv+\ngkoKCSLFRSFBLsjLgYu98d5QlNHDvjRz0Esyu+Z3N1RXw+hogA0SkZxSSJAL8rqSEIYyetiXZu6N\n99JcHdygRVdbbRtDZ4YYnxoHnN+DqSmYnAy4YSKSEwoJckF+hISghb2SEPQTIF1zayXMdjm4vwfq\nchApDgoJckGjo06Z2QthCgm98V7OXy08HBI2Qd9oSEJC7fkLKikkiBQXhQS5oEKsJDRXNzM5M8np\nM6eDbsoCg+ODTCemQ3Gd5lcS3LCocQkixUEhQS6oEENCmJdmDsNCSq5oeZTosqgqCSJFSiFBFjUz\nA+Pj3oSEmcQMJ0dPhiokhHHwYlimibrOXStBIUGkuCgkyKLcsrIXIWFgfIAZOxOKm1+Yl2Z2V4IM\nw+wGOH+tBIUEkeKikCCLckOCFwMXw/QJObosSmVpZSiXZu6N9xJdFqV6mUejRbPUVtum7gaRIqWQ\nIItybwZeVBLCFBKMMaFddTEs4zZc51YSKiud1zRwUaQ4KCTIovwICWEpo7fUtNA7GsKQMNobikGL\nrrZoGz2xHhI2QSTiVJVUSRApDgoJsiivQ0JdeR2VZZXZH8wDLTUtoe1uCFUlobaN6cQ0/aP9gJ7f\nIFJMFBJkUV6HhDDd/MK6NHNfvI+W6vBcp2SrLiokiBQHhQRZVCGHhLAuzRy267Sufh0AncOdgB7y\nJFJMFBJkUV7PbgjTzc/tbkjYRNBNmTM1M8WpsVOhuk4rqlZQXVbN4aHDgCoJIsVEIUEWFY9DeTmU\nlmZ/rLCFhObqZmbsDANjA0E3ZU7/WD8WG6qBi8YYLmq4SCFBpAgpJMiiCnFJZlcYl2Z2B1KG6ToB\ntDe0KyTWX5pGAAAgAElEQVSIFCGFBFmUVyFhYnqCoTNDobr5hXHVxTCtJXGu9nqFBJFipJAgi/Iq\nJJwcPQmE6+bnrtcQxpCwsnplwC05X3tDO0dOHyFhExq4KFJEFBJkUaOjhbcks6t6WTXRZdFQrZXQ\nG++lqbKJZSXLgm7Kedob2pmcmaQn1qNKgkgRUUiQRXlVSQhjSABCtzRzbzxcqy262hvaATg8dFgh\nQaSIKCTIorwMCRETYUXViuwP5qGwrZXQO9pLa01r0M1YYH39egCODB1RSBApIgoJsigvQ8KKqhWU\nREqyP5iHWmpaQldJCFu1BaCyrJLWmtbzKgnWBt0qEfGbQoIsysuQEMabX9iWZg7rdYLZaZCnD1Nd\n7QSEM2eCbpGI+E0hQRbl2cDF0XDe/ELX3RD2kDBbSQB1OYgUA4UEWVShVxJaalroH+1nOjEddFMY\nmxpjZGIklNcJFBJEipFCgiyq0ENCc3UzFjv3GOQghXW1RVd7Qzu98V5KK8cAhQSRYqCQIIvyIiRY\na0MbEsK0NHNYp4m6Lqq/CIAhjgBaUEmkGCgkSEqTkzA1lX1IiE/GGZsaC+XNL0xLM4c9JLhrJZya\ncZZnViVBpPApJEhK7k0g25AQ5ptfmJZm7o33UhoppbGyMeimJNUabaW8pJyTk04lQSFBpPApJEhK\nbjk529kNYQ4J5aXlNFQ0hGJp5t54L83VzURMOP9ZRkyEixouomv8IKCQIFIMwvluJKFQDJUECM/S\nzGEdt3GuK1Zcwb5TP6OkRCFBpBgoJEhKXoaE8pJy6srrsm+UD8KyVkJY15I41+bmzezt20tVtdXA\nRZEioJAgKXkZElpqWjDGZN8oH4RlaeZ8qCRsbtnM0Jkhqpq7VEkQKQIKCZKS1yEhrMKyNHPYrxM4\nlQSAkra9CgkiRUAhQVLybOBiyMvoYehuCPNaEudaW7eW+op6bLNCgkgxUEiQlNybQFVVdscJ+82v\nubqZwfFBJqYnAmvD6TOnmZyZDPV1AjDGsLl5M1ONCgkixUAhQVKKx50qQiTL35KwhwS3bSdHTwbW\nhrDPADnX5ubNjNfu1cBFkSKQ0du/MeZuY8wRY8y4MeZ5Y8zVS9zv14wxCWPMdzI5r+SWF0syJ2yC\nvnhfqG9+YViaOa9CQstmRiveYHhMKUGk0KUdEowxtwP3AZ8CtgB7gSeMMcsvsN864K+ApzNopwTA\ni5AwMDbAjJ0J9c0vDEsz51VIaN4MxnLS/CzopoiIzzKpJOwAvmatfcha+zrwYWAMuCvVDsaYCPAw\n8Ocw+3QYCT0vQkI+3PxWVK3AYAJddbE33kt1WTU1yzx45KbPrlh5BcaWMLhsb9BNERGfpRUSjDFl\nwDbgKfc1a60Ffghcs8iunwJOWmv/PpNGSjBGRwt7SWZXWUkZy6uWB15JCPM1OldFaQVN9lLi1QoJ\nIoWuNM3tlwMlwPyPXH3Apcl2MMa8HfgtYHParZNAeVlJcB+kFFZBL80c9mmi860q2cxrtQoJIoUu\n3ZCQigHsgheNqQG+AfyOtXYo3YPu2LGDurrzl/Ldvn0727dvz7SdkgavQkJdeR2VZZXeNMonQa+V\nkE+VBICLKrbxyvLvMT41Hvr/tyKFbNeuXezateu814aHhz07froh4RQwA8z/WLiShdUFgA5gHfC4\nObsmbwTAGDMJXGqtTTlGYefOnWzdujXNJopX4nFYuTK7Y+TLza+5upnO4c7Azt8b7+XitRcHdv50\nvbn+vTw6/sf88NC/8YHL3hd0c0SKVrIPzrt372bbtm2eHD+tMQnW2ingJeAG97XZm/8NwHNJdtkH\nXAW8Cae7YTPwGPCj2T8fz6jVkhOeVBLypIwe9PMb8iVMuS5r3AiD7Xz3tceDboqI+CiT7ob7gQeN\nMS8BL+DMdqgCvg5gjHkI6LLW3mutnQReO3dnY8xpnPGO+7JpuPjPq4GL+XDzC7K7YToxTf9of15c\nJ1dtrYEDH+D7q7+NtQ+E9uFdIpKdtKdAWmu/BXwU+AywB9gE3Gyt7Z/dZDWQP+92kpJXYxLy4ebX\nXN3MyMQI41PjOT93/2g/FpsX18kVjQIH3s+JsS5e7n056OaIiE8yGrhorX0AeCDFz959gX1/K5Nz\nSu4VU0g4d9XF9fXrc3rufJgmOl9tLXDsOqpLa3n8wONsad0SdJNExAd6doMkZe3ZZzdkamJ6gsHx\nwby4+bltDGJcQj6GhGgUmFnG1rqbefyAxiWIFCqFBElqbMwJCrW1mR/DfWBSPtz8glya2T3nyuos\np5LkUDTqfN9c8QF+2vNTemI9wTZIRHyhkCBJjYw437MJCfn0CbmpsokSUxLI0sy98V6aKptYVrIs\n5+fOlBsSLuEWosui/OUzfxlsg0TEFwoJklSxhYSSSAkrq1cGVknIh2t0rtJSqKwEO9bIn133Zzzw\n4gP8/OTPg26WiHhMIUGS8iokREyEFVUrvGmUz4Jamjlf1pKYLxp1fk/ueds9tDe0c8/378F5lIuI\nFAqFBEnKq5CwomoFJZESbxrls6DWSsjHSgI4ISEWg2Uly9h5805+dORHPPzKw0E3S0Q8pJAgScVi\nzne37zkT+Xbza64OqJKQZ9fJVVt79vfklktuYfuV27nze3fy8R9+nKmZqWAbJyKe8OoBT1Jg3EpC\nNiHhRPxEXt38WmpaePrY0zk/b76GBLe7wfXwbQ+zpWUL9/7oXh7b/xi3Xnor16+/njW1a2isbKSy\nrBKDwRhDxESS/tn9LiKZm05Me3YshQRJamQEysudr0z1xHq4fMXl3jXKZ253g7U2ZzeqsakxRiZG\n8jYkuJUEgIiJ8Cdv/xOuX389O5/fyd+//Pf85bOa9SCScx7OSFZIkKRGRrIbjwBOJeGGi2648IYh\n0VzdzNjUGPHJONHyLEooaXCnXOZjSKithd4kvTNvaXsLu355F9ZaDg8dpm+0j4GxAc5Mn8FisdZi\nsSRsIumfRSQ7R/cd5dN/92lPjqWQIEllGxKstZyInWBVdJV3jfLZuasu5iokuGMgWmtac3I+L0Wj\n8MYbqX9ujKGjsYOOxo7cNUpE2J3YzafxJiRo4KIkFYtlNx5hYHyAqcQUrdH8ufm5bT0RP5Gzc+bT\nWhLzzR+TICKFRyFBksq2knAi5txo8+kTslv1yOUSw73xXsoiZTRUNuTsnF6ZPyZBRAqPQoIklW1I\ncG+0+dTdEF0Wpbqsei7g5EJvvJfmmmYiJv/+KZ47BVJEClP+vTNJTmRdSZgt2edTGd0Yw6roqpxX\nEvLpGp0rGnWeFJpIBN0SEfGLQoIk5UUloamyifLSLOZQBqA12kpPPIchIU+XZIazY1bi8WDbISL+\nUUiQpGKx7Mck5NOgRVcglYTq/AwJ7u+HuhxECpdCgiQ1MpL9aov5NB7BtapmVc7HJOR7JUEzHEQK\nl0KCJOVFd0M+zWxwtUZbc1ZJsNYWREhQJUGkcCkkyAJTUzA+nv3AxbysJERXEZuMEZ/0v6P99JnT\nTM5M5m1IUHeDSOFTSJAF3Df9TEOCtTZvKwlusMlFl0M+L6QE6m4QKQYKCbKA+6afaUgYOjPE5Mxk\n3g5chNwsqOSeI99DgioJIoVLIUEWyLaS4H4Kz8fuBrf6kYuQ0B3rBvLzOoHzhNCyMoUEkUKmkCAL\nuJWETGc3uDfYfOxuiJZHqVlWk5PnN3SPdNNY2UhlWaXv5/KLVl0UKWwKCbJAtt0N7g02H7sbIHdr\nJXTHummLtvl+Hj/pIU8ihU0hQRbINiT0xHpoqGigorTCu0blUGtNbqZBdse6aavN/5CgSoJI4VJI\nkAVGRsAYqK7ObP8Tsfyc/uhaFV2Vk+6GnlhP3lcS1N0gUtgUEmSBWMz5hBjJ8LfjRDw/l2R25ay7\nYUTdDSISbgoJsoAXqy3mcyUhF90NM4kZeuO9eX2dQN0NIoVOIUEW8OK5Dfk4s8G1KrqK+GSc2IR/\nd7++0T5m7Ezej0lQd4NIYVNIkAWyqSTk82qLrrlVF30cl9A94qyRoO4GEQkzhQRZIJuQMDg+yJnp\nM6yuXe1to3IoF6suugsp5XslQd0NIoVNIUEWiMUyDwldI10ArKlb42GLcssddOnn8xu6R7opi5Sx\nvGq5b+fIBYUEkcKmkCALZFNJOD5yHCCvKwk1y2qILovOfdr3gzu4M2Ly+5+gOybB2qBbIiJ+yO93\nKPFFNgMXjw8fpzRSSnN1s7eNyrG22ra5cQN+KISFlMD5PUkkYGws6JaIiB8UEmSBbCoJXSNdrIqu\noiRS4m2jcmxN7Zq5qogfumPdeT/9EfQkSJFCp5AgC2Tb3bCmNn/HI7h8DwkFsJASnP09UUgQKUwK\nCXIea7MfuJjP4xFca+rWcHzY30pCIYQEt5KgaZAihUkhQc4zPg4zM6okrKldQ2+8l8mZSc+PHZ+M\nMzIxUhBjEtzfE4UEkcKkkCDnyeYJkNbagqokWKwvayUUykJKAA0NzvfBwWDbISL+yCgkGGPuNsYc\nMcaMG2OeN8Zcvci2v22MedoYMzj79eRi20uw3JCQyeyGgfEBzkyfyes1Elxu0PGjy6FQFlICqK93\nHgQ2MBB0S0TED2mHBGPM7cB9wKeALcBe4AljTKpVYa4HHgHeBbwNOA78wBiTv+v2FrBsKgnuQkoF\nUUmY7TLxY/CiW50ohNkNJSVONUEhQaQwZVJJ2AF8zVr7kLX2deDDwBhwV7KNrbW/aa39qrX2FWvt\nAeC3Z897Q6aNFv9kExLcT92FMCYhWh6lrrzOn0rCSDf1FfVUlVV5fuwgNDUpJIgUqrRCgjGmDNgG\nPOW+Zq21wA+Ba5Z4mGqgDFAvZgi5U9kyrSSURkpZWb3S20YFZE2dP9MgC2XchqupCU6dCroVIuKH\ndCsJy4ESoG/e631AyxKP8UWgGydYSMi4A9Dq6tLf9/jIcdqibXm/kJJrTe2auS4ULx0bPsa6unWe\nHzcoqiSIFC6vZjcY4IKrtxtjPg78KvAha633c8ska/39zmC0srL09z0+crygPiH7taCSQoKI5IvS\nNLc/BcwA8xfmX8nC6sJ5jDF/DHwMuMFa++pSTrZjxw7q5n2k3b59O9u3b19ygyU9p07BihWZ7ds1\n0lUQMxtca+rW8N3Xv+v5cY+dPsYdV93h+XGDopAgEpxdu3axa9eu814bHh727PhphQRr7ZQx5iWc\nQYePARhjzOzfv5JqP2PMnwD3AjdZa/cs9Xw7d+5k69at6TRRstTfD8szfHrx8eHjvLn1zd42KEBr\natfQP9bPmekzVJRWeHLM4TPDDE8MF1QlYflyhQSRoCT74Lx79262bdvmyfEz6W64H/hdY8ydxpjL\ngK8CVcDXAYwxDxljvuBubIz5GPBZnNkPncaY5tmv6qxbL57LtJLgLqRUSJUEt+vEy3EJx4aPAbCu\nvnBCQlMTDA05K3WKSGFJOyRYa78FfBT4DLAH2ATcbK3tn91kNecPYvwIzmyGfwJ6zvn6aObNFr9k\nWkk4NXaKiZmJwhqTMBt4vJwGeez0bEgooEpCU5PzuOjTp4NuiYh4Ld0xCQBYax8AHkjxs3fP+/tF\nmZxDgtHfn1klwf20XQhrJLjmVl30cPDiseFjLCtZRnPN/GE9+aupyfk+MHD2zyJSGPTsBjnPqVOZ\nVRLcG2khVRKqyqpoqmzyvJKwtm4tEVM4//TODQkiUlgK551KsjYx4ay4mEkl4fDQYSpKK2ipWepy\nGfnB6wWVCm36IygkiBQyhQSZ477JZ1JJODx0mPaGdpzJLoXD67USFBJEJJ8oJMic/tmhp5lUEg4N\nHaK9od3bBoXAmto1dA53ena8Y6ePFdTMBoDycqipUUgQKUQKCTLHXX8/0+6GjoYObxsUAu0N7RwZ\nOoLziJLsnJk+Q99oX8FVEkALKokUKoUEmeNWEtLtbkjYBEeGjhRkJaG9oZ3RqVFOjp7M+lhuRaLQ\nKgmghzyJFCqFBJnT33+2dJyOnlgPEzMTBRkSOhqd6sjhocNZH6sQ10hwqZIgUpgUEmSOO/0x3bGH\n7g20ELsbLqp3lvk4NHQo62MdGz5GxEQKapqoSyFBpDApJMicTBdSOjTo3EDX16/3tkEhEC2PsrJ6\npWeVhFXRVZSVZPCIzZBTSBApTAoJMifThZQODx2mLdpGZVml940KgfaGds8qCYXY1QB6yJNIoVJI\nkDkZVxIKdPqjq6Ohw5tKwnDhTX90uZUEDyaBiEiIKCTInEwf7uQupFSoOho65rpUsnH09NGCrSQ0\nNTkrdo6OBt0SEfGSQoLMyfQx0YW6RoKrvaGdE/ETjE2NZXyM8alxjg8f5+LGiz1sWXho1UWRwqSQ\nIIDzqN9MxiTEJmL0j/UXdiVhdhrkkaEjGR/j4OBBLJZLmi7xqlmhopAgUpgUEgSA4WGYmUm/kjA3\n/bGxsCsJkN00yAMDBwC4dPmlnrQpbBQSRAqTQoIAma+26N44C7mS0FrTSmVpZVaDF/cP7KehooGm\nyiYPWxYeCgkihUkhQYDMn9tweOgw1WXVrKjKYDBDnjDGONMgsxi8eGDgAJcuv7TgnpLpikahrEwh\nQaTQKCQIkEUlYfAQHY0dBXvzc7U3tHP4dOaVhAMDBwp2PAI4q3RqQSWRwqOQIMDZkNCUZjX89YHX\nC/rm58p2GuT+gf1c0ljY10kPeRIpPAoJAjhv7o2NUFq69H2stbzS9wqbVm7yr2Eh0d7QzpHTR0jY\nRNr7DowNMDg+WLCDFl2rV0NnZ9CtEBEvKSQIkNlCSj2xHgbHB9nUXPghoaOxg8mZSbpGutLe153Z\nUOgVl44OOJT9mlMiEiIKCQJktpDSK32vABRFSLhs+WUAvNb/Wtr77h/YD8CGxg2etils3JCgpZlF\nCodCggBOmXjVqvT2eaXvFaLLogX7PIJzra9fT82yGn7W97O09z0wcIC1dWupKqvyoWXh0dEB4+Nw\n4kTQLRERrygkCAD798OlaXaZv3LyFa5qvoqIKfxfo4iJcOXKK/nZyfRDwv6B/QXf1QCwYbZQoi4H\nkcJR+O/uckEjI86nv7RDQpEMWnRdtfKquS6WdBwYOFDwMxsA2mfX01JIECkcCgnCAWdcXVohYWJ6\ngtdPvV4U4xFcV628in2n9jE1M7XkfRI2wRsDbxT8zAaAykqny0ohQaRwKCQI+51xdWmFhNdPvc50\nYrqoQsKm5k1MzkzyxuAbS96nc7iTiZmJouhuAM1wECk0CgnC/v3Q2gq1tUvfx+2bv3LllT61Knyu\nar4KIK0uB3eg48blG31pU9goJIgUFoUEyWzQYt8rrK9fT11FnT+NCqHGykZWRVelNcPhxZ4XWVG1\ngrV1a31sWXgoJIgUFoUEyTgkXLXyKn8aFGJXrbwqrRkOL/a8yFva3lLwz7ZwdXQ4z28YHg66JSLi\nBYWEIpdIOAMXM5rZUETjEVybmjctubvBWssL3S9w9aqrfW5VeHR0ON9VTRApDAoJRa6ry1kAJ52Q\ncGjwECfiJ4rq5ue6auVVHBs+xsjEyAW3PXL6CIPjg1zdVjzXSSFBpLAoJBS5TGY2PHHoCUojpfzC\nRb/gT6NCzB28+POTP7/gti92vwhQVGGqsRHq6hQSRAqFQkKR278fli2D9euXvs8Th57g2jXXUlue\nxnSIArFx+UZKTMmSuhxe7HmRdXXrWFGd5kMx8pgxzsqLCgkihUEhocjt3++8qZeULG37yZlJfnTk\nR9zccbO/DQup8tJyNjVv4tnjz15w2xd7XiyqrgaXZjiIFA6FhCKX7syG544/R3wyXrQhAeDG9ht5\n8tCT2EUedziTmOGlnpd4y6q35LBl4dDRAW8sfb0pEQkxhYQil25IeOLgE6yoWsGW1i3+NSrkbuq4\nib7RvkWnQu47tY/RqdGirCRcc40zIHbfvqBbIiLZUkgoYj09ziOir0xj0cQnDj3BjR03FsWTH1N5\n+9q3U1FawQ8O/SDlNi92v4jBsK11Ww5bFg433gjV1fDd7wbdEhHJVvG+0wv/8A9QUQHvf//Stu+L\n97Gnd09RdzUAVJRWcP2663ny8JMpt/m3o//GFSuvIFoezWHLwqGiAm65Bb7znaBbIiLZUkgoUtbC\ngw/Chz7kTFlbim/+/JtETISbOm7yt3F54KaOm3j62NOcmT6z4Gejk6N89/Xv8quX/2oALQuH226D\nl16CY8eCbomIZEMhIeR27drly3H37IFXX4U771za9vHJOF945gvcuflOWmpafGnTYvy6Dpm6sf1G\nzkyf4ZnOZxb87NH9jxKfjHPHpjt8OXfYrkUy73ufM7XWzy6HfLgOuaJr4dB18F5GIcEYc7cx5ogx\nZtwY87wxZtHRWcaYXzHG7Jvdfq8x5r2ZNbf4+PVL/+CD0NLi9B8vxVd+8hWGxof41PWf8qU9FxK2\nf/xXrrySlpqWpOMSHn7lYa5dcy3tDe2+nDts1yKZaBRuusnfLod8uA65omvh0HXwXtohwRhzO3Af\n8ClgC7AXeMIYszzF9tcAjwD/L/Am4HvA94wxl2faaMnO1BQ88gjccQeUll54+8HxQb707Jf48Js/\nzPr69b63Lx8YY7ip4ya+s+87TExPzL3eF+/jB4d+wG9u+s0AWxcOt90GzzzjDJAVkfyUSSVhB/A1\na+1D1trXgQ8DY8BdKba/B/hXa+391tr91tpPAbuB38+oxZK1//pf4dSppXU1TCem+egPPspUYoo/\nfeef+t+4PPLRaz7KseFjfPHZL8695o7b+JXLfyXAloXDrbdCQwP88i9DLBZ0a0QkE2mFBGNMGbAN\neMp9zToryvwQuCbFbtfM/vxcTyyyvfgkkYAdO+Dee+HjH4dNF3iIY1+8jxu/cSPf2PsNvvyeL9Nc\n05ybhuaJTc2b+Ni1H+Pz//F59vXv4/VTr/M3L/4Nt1xyC01VTUE3L3CNjfDEE87Ylw9+0HmQmIjk\nlyUUm8+zHCgB+ua93gekWpKnJcX2i41+qwC476Fv0/yjF9NsosMm+VNa+yfZLa0jLbrx0o5kLfz8\n4FH+4Iv/fZEjpT5WwsLkBIyNQ9dxZxW8EyfgV/8I1rwNHnj0/O0nE5OMTo7SP9bP3r69/LTnpyyL\nLOOrN36VrWxl9+7dS2q3H4aHhwM9fyrvj76fh2MPc8MXb6B/tJ+WmhZuv/R2X9sa1muRTCQCO3fC\n7/8+rFrlLLS0ebMzo6a62vk5OM98cLl/Pve1ZD87dmyYv/u7/LgOftO1cOg6ODo751Yyq8j2WGax\npWUXbGxMK9ANXGOt/ck5r38JeIe19tok+0wAd1pr//Gc134P+KS1dlWK8/w68A9LbpiIiIjMd4e1\n9pFsDpBuJeEUMAPMrzuvZGG1wNWb5vbgdEfcARwFFk5EFxERkVQqgPU499KspFVJADDGPA/8xFp7\nz+zfDdAJfMVa+1dJtv8mUGmtvfWc154F9lprfy+bxouIiIh/0q0kANwPPGiMeQl4AWe2QxXwdQBj\nzENAl7X23tntvwz8uzHmj4D/DWzHGfz4O9k1XURERPyUdkiw1n5rdk2Ez+B0I7wM3Gyt7Z/dZDUw\nfc72PzbGbAc+P/v1BnCrtfa1bBsvIiIi/km7u0FERESKg57dICIiIkkpJIiIiEhSoQsJ6T48qhAZ\nYz5hjHnBGDNijOkzxnzXGHNJ0O0K2ux1SRhj7g+6LblmjFlljPmGMeaUMWZs9kFpW4NuV64ZYyLG\nmM8aYw7PXoeDxphPBt0uvxlj3mmMecwY0z37b+CDSbb5jDGmZ/a6PGmM2RBEW/222LUwxpQaY75o\njHnFGBOf3ebB2TV+CspSfifO2fZrs9v8QbrnCVVISPfhUQXsncBfA28FfhEoA35gjKkMtFUBmg2L\nv4PzO1FUjDH1wLPABHAzsBH4KDAUZLsC8nHgvwC/B1wGfAz4mDGm0J8FU40zSPxukizZaoz5v3Ge\nh/NfgLcAozjvncty2cgcWexaVOE8SPDTOPeQX8JZDXje+rIFYdHfCZcx5kM4vxPdmZwkVAMXU6zB\ncBxnDYYvBdq4AM2GpJPAddbaZ4JuT64ZY2qAl4CPAH8G7LHW/lGwrcodY8xf4qxyen3QbQmaMeZx\noNda+zvnvPZPwJi1dgmPLMt/xpgE8CFr7WPnvNYD/JW1dufs32txFqz7T9babwXTUv8luxZJtnkz\n8BNgnbW2K2eNy6FU18EY0wb8GOfDxb8AO621X0nn2KGpJGT48KhiUY+TFAeDbkhA/hZ43Fr7o6Ab\nEpAPAD81xnxrtvtptzHmt4NuVECeA24wxlwMYIzZDLwd5w2wKBljLsJ5Fs65750jODfGYn/vhLPv\nn6eDbkguzX7Ifgj4krV234W2TyWTxZT8ksnDowre7P/o/wY8U4xrSxhjfg2nfPjmoNsSoHacKsp9\nOGuNvBX4ijHmjLX24UBblnt/CdQCrxtjZnA+6PyptfabwTYrUC04N8F0H6RX8Iwx5Ti/M49Ya+NB\ntyfHPg5MWmv/JpuDhCkkpGLI9FGOheEB4HKcT0tFxRizGicg3WitnQq6PQGKAC9Ya/9s9u97jTFX\n4ASHYgsJtwO/Dvwa8BpOgPyyMabHWvuNQFsWPkX93mmMKQX+F841KKpHABhjtgF/gDMuIyuh6W4g\ns4dHFTRjzN8A7wPeZa09EXR7ArANWAG8ZIyZMsZMAdcD9xhjJmerLMXgBDC/XLgPWBtAW4L2JeAv\nrLX/y1r7qrX2H4CdwCcCbleQenECgd47Z50TENYANxVhFeEdOO+dx89571wH3G+MOZzOgUITEmY/\nKb4E3OC+NnsTuAGnH7KozAaEW4FfsNZ2Bt2egPwQuArn0+Lm2a+f4nx63mzDNOrWX8+ysMvtUuBY\nAG0JWhULPx0nCNF7Wa5Za4/gBIVz3ztrcbqlivG90w0I7cAN1tpinAX0ELCJs++bm4EenJB9czoH\nClt3w6IPjyoWxpgHcB6E9UFg1BjjfkIYttYWzaOzrbWjOCXlOcaYUWAgm4E4eWgn8Kwx5hPAt3De\n/HYbTO4AAAE1SURBVH+b4nxI2uPAnxpjjgOvAltx3if+R6Ct8pkxphrYgFMxAGifHbQ5aK09jtMt\n90ljzEHgKPBZoIsCnPq32LXAuRF+G+eDxfuBsnPePwcLqdtyCb8TQ/O2n8KZGfRGWiey1obqC6fv\n6CgwjjN1481BtymAa5DA6XqZ/3Vn0G0L+gv4EXB/0O0I4L/7fcArwBjOzfGuoNsU0HWoxvkwcQRn\nLYA3cObElwbdNp//u69P8b7w/52zzf+Dc5McA54ANgTd7lxfC5yS+vyfuX+/Lui25/p3Yt72h4E/\nSPc8oVonQURERMKjaPvxREREZHEKCSIiIpKUQoKIiIgkpZAgIiIiSSkkiIiISFIKCSIiIpKUQoKI\niIgkpZAgIiIiSSkkiIiISFIKCSIiIpKUQoKIiIgk9f8DxUc9xN3IiooAAAAASUVORK5CYII=\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x7fe5f1e1e590>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"\n",
"\n",
"k=2\n",
"\n",
"#Fails for t0 < -5 or t0 > 5\n",
"t0 = 0\n",
"t = linspace(0,14,200)+t0\n",
"\n",
"args = {'a':5,'b':0.2,'c':5+t0}\n",
"\n",
"\n",
"H = [[sqrt(k)*tx,f]]\n",
"rho0 = b[0]*b[0].dag()\n",
"c_ops = []\n",
"obs_ops = [b[0]*b[0].dag(), b[1]*b[1].dag()]\n",
"output = qutip.mesolve(H, rho0, t, c_ops, obs_ops ,args=args)\n",
"\n",
"plot(t,f(t,args))\n",
"plot(t,output.expect[1])\n",
"ylim(0,1.2)"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"# Time-dependent hamiltominan given in list array format"
]
},
{
"cell_type": "code",
"execution_count": 110,
"metadata": {
"collapsed": false
},
"outputs": [
{
"data": {
"text/plain": [
"(0, 1.2)"
]
},
"execution_count": 110,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAgkAAAFkCAYAAACq4KjhAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAAPYQAAD2EBqD+naQAAIABJREFUeJzt3Xt8XGd97/vPo6t1tyRbki0ptiU7cUJiJzYhJFzCJiSh\nkJY0Pd3UQOHV7F6A9JS6cDgNpeWUS7tpSVzYbVpa9i4JFBd2yyWhLSFAN4GEQIgdBxLHju+yLMm6\n2NKM7tI854+lJcvSjKyZWWvWmjXf9+ull+LRujxekWd957n8lrHWIiIiIrJYUdANEBERkXBSSBAR\nEZGkFBJEREQkKYUEERERSUohQURERJJSSBAREZGkFBJEREQkKYUEERERSUohQURERJJSSBAREZGk\nFBJEREQkKYUEERERSUohQURERJJSSBAREZGkSoJuQDLGmEbgduAEMBFsa0RERPLKKmAj8Ki1djCb\nA4UyJOAEhH8KuhEiIiJ57O3Al7I5QFhDwgmAL37xi1x55ZUBNyVYu3fvZs+ePUE3I3C6DhfoWjh0\nHS7QtXDoOjgOHjzIO97xDpi7l2YjrCFhAuDKK69kx44dQbclUHV1dQV/DUDXYSFdC4euwwW6Fg5d\nhyWyHq7XxEURERFJSiFBREREklJIEBERkaQUEkJu165dQTchFHQdLtC1cOg6XKBr4dB18J6x1gbd\nhiWMMTuAZ5555hlNQhEREUnDvn372LlzJ8BOa+2+bI6lngQRERFJSiFBREREklJIEBERkaQUEkRE\nRCQphQQRERFJSiFBREREklJIEBERkaQUEkRERCQphQQRERFJSiFBREREklJIEBERkaQUEkRERCQp\nhQQRERFJSiFBREREklJIEBERkaQUEkRERCQphQQRERFJSiFBREREklJIEBERkaQUEkRERCQphQQR\nERFJSiFBREREklJIEBERkaQUEkRERCQphQQRERFJSiFBREREklJIEBERkaQUEkRERCQphQQRERFJ\nSiFBREREklJIEBERkaQUEkRERCQphQQRERFJSiFBREREklJIEBERkaQUEkRERCQphQQRERFJSiFB\nREREklJIEBERkaTSDgnGmNcYYx42xnQbYxLGmF9awT6vM8Y8Y4yZMMYcNsa8K7PmioiISK5k0pNQ\nBTwL3APYS21sjNkIfBP4LrAd+DTwOWPMrRmcW0RERHKkJN0drLXfAr4FYIwxK9jlPcAxa+0H5/58\nyBjzamA38Fi65xcREZHcyMWchFcC31n02qPAjTk4t4iIiGQo7Z6EDLQAfYte6wNqjTHl1trJHLRB\nciyRgCefhAMHYHAQhoed1wCsdb6W+++Ux2WGGTNOma1Juc2K+rc83C+s53y+7H9RP7uV9bM3zb9W\nVATV1VBbC9dcA699LVRUZN4GEYm2XISEZNy3t2VvCbt376auru6i13bt2sWuXbv8apdkKZGAP/1T\n+J//E7q7obQUGhqgrg6Kiy/c2Iy5+Gvxa6mc2HE359q/QMlEM1XnXsmGZ75I8Wz1JcPFYulun8k+\nQbZpcMtf0Xvtbopmqtjw/e9QOfRKAGZnIR6H8+ed7+XlcNddzv8vhQWR/LN371727t170WvDw8Oe\nHT8XIaEXaF70WhMwYq2dWm7HPXv2sGPHDt8aJt77+MfhYx+D97wHdu2Cm25yPr164cWBF7nqb77I\nPdffQ0NFA3/2gz/jt+9+kHtecY83J4iIvT/by9u+upvfv+H3+WnPT3m++k08/BuPc3XT1fPbWAsH\nD8K//zv8yZ/AyAh89atQVhZgw0Ukbck+OO/bt4+dO3d6cvxchIQfAb+w6LXb5l6XCHn4YfjIR5yQ\n8OEPe3/8jz/+cVprW7nvtvsoLynn4MBBPvOTz/Ce699DkVHJD4Bj547xrq+/i1/f9uvcd/t9jEyO\n8LrPv45f/vIvc/h3D+PONTYGrrrK+dq2DX7xF+Htb4d//menx0dEBDKrk1BljNlujLl27qWOuT+3\nz/38z40xDy7Y5e+ATmPMJ40xVxhj3gv8X8D9WbdeQuP4cXjHO+CXfxk+9CHvj3948DB7f76Xe199\nL+Ul5QC874b3cXjwMN8++m3vT5in/vWFf6W4qJi/ffPfUmSKWL1qNffddh9Hho7wbO+zSfe57Tb4\n8pfhX/4FHnkkxw0WkVDL5OPXy4H9wDM4cwruA/YBfzr38xag3d3YWnsCeDPwBpz6CruB/2atXbzi\nQfLYl77kfH/wQe+GFxb6xA8+QUt1C3dfd/f8a69qfxXXtVzHp3/8ae9PmKe+cegb3NpxK1VlVfOv\nvXbDa1m9ajVff/HrKfe78064+mpnyEFExJX227m19vvW2iJrbfGir7vnfv4b1trXJ9lnp7W2wlq7\nxVr7Ba/+AhIO3/oW3HIL1KRedJCxsekx9v5sL7tfuZtVJavmXzfG8L4b3se3jnyLQwOHvD9xnjk7\nepYnu57kLVe85aLXS4tLefOWN/P1Q6lDAjgTGB95BKaWnSkkIoVEA7mSteFh+NGP4Pbb/Tn+U6ef\nYjoxzRs3v3HJz9569VupX1XP3p/vTbJnYfnm4W8CcMfldyz52Z1b7+S5vuc4fu54yv3vustZ9fB/\n/o9fLRSRfKOQIFn77nedpXV+hYTvn/g+jRWNXLX2qiU/W1Wyildf9mqe6HrCn5PnkW8c+gY3tt9I\nc/XixURwe+ftlBWX8Y1D30i5/7Zt0NGhIQcRuUAhQbL26KNw+eWwaZM/x3/81OO8ZsNrUq5guKn9\nJp46/RSziVl/GpAHxqbHeOzoY0uGGlw15TW8oeMNy85LMMbpTfj6153QJyKikCBZsdaZj/DGpSMB\nnpicmeSp00/x2stem3Kbm9pvIj4V5+dnf+5PI/LAY0cfY3xmPGVIALjzijv5wakfMDA2kHKbu+6C\nvj5n+EhERCFBsnLoEJw65V9I+OmZnzIxM8FrN6QOCS9f/3JKikp4sutJfxqRB757/Lt01HdwxZor\nUm7zxs1vJGETPHEq9dDMDTdAS4uWQoqIQyFBsvKtbzmlfW++2Z/jP37ycWrKatjesj3lNpWllexY\nt6Og5yXs793Py9e/fNlt2mrbWFO5hv29+1NuU1QE113nVGMUEVFIkKw88QS88pVQWenP8R8/9Tiv\nuuxVlBQtXxz0prabCrYnIWETHOg9wHUt1y27nTGG61quWzYkAHR2wtGjXrZQRPKVQoJk5fhxZ9Ki\nH2YSMzxx6glu3nDpboqb2m/i+Pnj9MR6/GlMiB07d4zYVOySIQFwQkLPpUPCsWMXntopIoVLIUGy\ncuIEbNzoz7EP9B4gNhVbdj6C66Z253HIhdib4N70r2259hJbwnXrrqNrpIvBscGU23R2wsQE9BRe\n3hKRRRQSJGPxOAwOwoYN/hx/f+9+ikwRO9Zd+kmgrbWtbKjbUJghoXc/66rXJa2PsJjb27DckENn\np/NdQw4iopAgGTt50vnuV0/CC/0v0FnfeVEp5uXc2H4jPzpdeGv3nu19luvWXXqoAWBL4xaqSquW\nHXLo6HBqJigkiIhCgmTsxAnnu58hIVmVxVS2NW3jhf4XsNb606CQ2t+7f0XzEQCKTBHbW7Yv25Ow\nahW0tiokiIhCgmThxAkoLYV16/w5/vP9z6cVErau2crw5DC98V5/GhRCvfFeeuO9Kw4JgFY4iMiK\nKSRIxk6cgMsu8+fR0COTI5weOc3L1r5sxftcufZKAF4ceNH7BoXUs73PAiubtOi6ruU6Dg0cYnRq\nNOU2nZ1w5EjWzRORPKeQIBk7edK/oYaD/U41n3R6EjrrOykpKuHgQOFUAtrfs5/a8lo21a/8wRnX\nrbsOi+W5vudSbqOeBBEBhQTJgp/LH5/vfx6DWbbM8GKlxaVsbthcUD0J+3v3c23LtSkffpXMy9a+\njJKikkuucDh3zvkSkcKlkCAZO3HCv+WPL/S/wKb6TVSWplfK8co1VxZUT8KBvgNc27zyoQaA8pJy\nrlxzJT/r+1nKbbQMUkRAIUEyNDYG/f3+rmxIZz6C68o1V84PVUTdTGKGY+eOpdXb4rq88XKOnEs9\n6UAhQURAIUEy5HeNhHRXNri2rtlKd6yb2GTMh1aFy6nhU8wkZtjSsCXtfbc0bOGlwZdS/ry+3vlS\nSBApbAoJkhG3RoIfww2xyRinhk9lFBIKaYXDkSGnJ2Bzw+a0993csJlTw6eYnJlMuY0mL4qIQoJk\n5MQJKCmB9eu9P7Z7g89kuGHrmq0XHSPKXhp8idKiUtrr2tPed0vjFiyWY+eOpdxGIUFEFBIkIydP\nQnu7ExS89nz/88CFG346qsuqaattK4jJi0eGjrCpftMlH6OdjNv74PZGJKOQICIKCZIRP5c/vtD/\nAhtXb6SqrCqj/QtlhcORc0cymo8AsK56HZWllbw0lHpeQksLnD0LBVblWkQWUEiQjPi5/PHw4GGu\naEx/xr7ryjVXFsxwQybzEQCMMWxu2LxsT0JjI0xNwWjqwowiEnEKCZIRP6stHj9/nE2rV15BcLGt\na7ZyZOgI07PTHrYqXGYTsxw7dyzjkABzKxyW6UlobHS+Dw5mfAoRyXMKCZK2mRno7YW2Nu+Pba3l\n+LnjaZUZXmzrmq3zNQSiqmuki+nEdFYhYSU9CaCQIFLIFBIkbefPO9/dm4iXhsaHiE3FsupJcAPG\nyeGTXjUrdNybe6ZzEtx9l1sGuWaN810hQaRwKSRI2oaGnO8NDd4f+/j54wBZ9SS01bZRbIo5fu64\nV80KnZcGX6KkqIQNqzOfGLK5YTMJm5i/5oupJ0FEFBIkbb6GhLkbe0d9R8bHKCkqoa22jRPnT3jU\nqvA5MnSEjas3ZrT80bWl0emFSFV5sboaSksVEkQKmUKCpM3vnoTa8lrqV9VndZxN9Zs4MXzCm0aF\nUDbLH13uMshU8xKMcXoTFBJECpdCgqTNDQn12d3Hkzp+zlnZYIzJ6jgbV2+MdE9CNssfXe4yyEut\ncFBIEClcCgmStqEhqKhwvrx2/Hx2KxtcG+uiGxJmE7McPXc065AAK1vhMDCQ9WlEJE8pJEjahob8\n6UWA7GskuDau3khvvJfx6XEPWhUu3bFupmanvAkJ9epJEJHUFBIkbUND/sxHSNgEJ86f8CwkQDSX\nQbqTO724TpvqN3F65DSzidmkP1+zRiFBpJApJEja/AoJZ2JnmJqd8mS4wT1GFIccTg2fAuCyusuy\nPtZldZcxk5ihJ96T9OfqSRApbAoJkja/QoKXn5DX16ynpKgkkiGha6SLxorGjB+AtZAbNNzgsZhC\ngkhhU0iQtPkWEuaK+rhDBdkoKSqhvbY9kiHh1PApT3oRYGUhYWQEpqP7GAwRWYZCgqTt3Dn/ehKa\nqpo8+YQMTthIVU0wn3kZEmrLa1m9avWyIQEuLHsVkcKikCBp87MnIZtKi4tFtVbCqeFTtNe2e3a8\ny+ouu2RI0JCDSGFSSJC0JBL+hgQv5iO4Nq3eFNmQ4FVPAjghIdUqEPchT6qVIFKYMgoJxph7jDHH\njTHjxpinjDHXX2L73zfGvGiMGTPGnDLG3G+MKc+syRKkWMwJCn4NN3gZEjau3sjZ0bOMTY95dsyg\nDU8ME5uKeRsSatWTICLJpR0SjDFvBe4DPgJcBxwAHjXGrEmx/duAP5/bfitwN/BW4BMZtlkC5Ndz\nG2YSM3THurN6quFi87USzkenVoKXyx9dyw031Nc7z3BQSBApTJn0JOwGPmutfcha+yLwbmAM5+af\nzI3AD621X7bWnrLWfgfYC7wioxZLoPwKCT2xHhI24elYuxsSojR50a+QcH7iPCOTI0t+VlwMq1cr\nJIgUqrRCgjGmFNgJfNd9zVprge/ghIFkngR2ukMSxpgO4E3Av2XSYAmWXw93Oj1yGoC22jbPjrm+\nZj2lRaXz9Rei4NTwKUqKSmipbvHsmG7vTddwV9Kfq1aCSOFKtydhDVAM9C16vQ9I+q5lrd2LM9Tw\nQ2PMFPAS8J/W2k+meW4JAb96EvwICcVFxbTWts4fOwpODZ+irbaN4qJiz47p9kqkmryokCBSuLxa\n3WAAm/QHxrwO+BDOsMR1wF3AHcaYD3t0bsmhoSGnC7qmxtvjnh45TWVpJatXrfb0uG21bZyORSck\ndI10eTokA7Cueh3FplhVF0VkiZI0tx8AZoHmRa83sbR3wfVR4CFr7T/O/fl5Y0w18Fng48udbPfu\n3dTV1V302q5du9i1a1eazRavuMsfjfH2uKdHTtNW24bx+MCtNa10j3R7eswgeb38EZwel7batpQh\nYc0aOHbM01OKiEf27t3L3r17L3pteHjYs+OnFRKstdPGmGeAW4CHAYzzrn4L8JkUu1UCiUWvJeZ2\nNXNzGpLas2cPO3bsSKeJ4jO/aiScjp32dKjB1Vbbxr6efZ4fNyinhk/x6ste7flxL1VQ6Sc/8fyU\nIuKBZB+c9+3bx86dOz05fro9CQD3Aw/OhYWf4Kx2qAQ+D2CMeQg4ba390Nz2jwC7jTHPAj8GtuD0\nLnxjuYAg4eRbSBg5zeaGzZ4ft622je5YN9Zaz3spcm02McvpkdOe9ySAM3kx1VJRDTeIFK60Q4K1\n9itzNRE+ijPs8Cxwu7W2f26TNmBmwS4fw+k5+BjQCvTj9EJoTkIe8uu5DadHTvO6Da/z/LitNa2M\nTY9xfuI89RUeL8nIsZ54D7N21peQcFntZfzg5A+S/qyx0QmH1no/zCQi4ZZJTwLW2geAB1L87PWL\n/uwGhI9lci4Jl6Eh2LjR22POJmY5Ezvj23ADQHesO+9DgrtE0ZeQUHcZp0dOM5OYoaTo4reFxkaY\nnYXhYadmgogUDj27QdLix3DD2dGzzCRmfAkJrbWtAJFYBunOGfB6dQM4IWHWztIT61nyM5VmFilc\nCgmSFj9CgnsDb6/z/ua3rnodBhOJFQ6nhk9RW15L3aq6S2+cJrd3ItnkRfchTwoJIoVHIUFWzFp/\nQ4IfPQmlxaU0VzdHoifBXSbqB7fHpTu2NEy5Qwznz/tyahEJMYUEWbHxcZic9D4kdI10UV5cTmNF\no7cHnuOucMh3Z+JnaK1p9eXYdeV1VJVWJQ1TbuGsWMyXU4tIiCkkyIr5+dwGPwopuVprolGauXuk\nm/U16305tjGG1trkhacUEkQKl0KCrJifz23wqxsd5kozRyAknIn515MAc9Upk/S4lJRARQWMLH1I\npIhEnEKCrFg+h4R8H26w1nImdsa3ngRw5iWkuk41NepJEClECgmyYn4PN/iltaaVofEhxqfHfTuH\n3wbHB5lOTPsbEpZ5zoVCgkhhUkiQFXO7m2trvTtmwibojnX73pMAyWfu5wv35u17SJgrYb1Yba2G\nG0QKkUKCrFg8DmVlzpdXBsYGmJqd8rcnIQIFlc7EzgA+h4TaVqZmpxgYG1jyM/UkiBQmhQRZsVjs\nwkx3r/hZI8HlTvbL54JKZ2JnMBhaqlt8O8dyPS4KCSKFSSFBVszPkODnrP2qsipWr1qd9z0JTVVN\nlBaX+naO5cKUhhtECpNCgqyYHyGhN95LkSmiqarJ2wMvku8rHLpj/tVIcDVXN1NkitSTICLzFBJk\nxfwKCWsr11JcVOztgRfJ94JKfi9/BCgpKqGluiVlQSWFBJHCo5AgK+ZXSPBznN2V7wWV/C6k5EpV\nUEnDDSKFSSFBViyfQ8K66nX0xJc+Bjlf5KInAVIXVFJPgkhhUkiQFcvrkFCzjt54Lwmb8P1cXptJ\nzNA32pebkJBiWKamxlkCm8i/yyciWVBIkBWLx/M4JFSvYyYxw+DYoO/n8lpfvI+ETeQkJLTVti37\nkKd43PcmiEiIKCTIisViUF3t3fGstTntSQDycsghF4WUXK01rZybOLekhLVbZVNDDiKFRSFBVszr\n4YbYVIzxmfGc9SQA9MTyNyS4lSP95J5j8bwEPS5apDApJMiKeR0SeuO9ADkJCe458rUnoaSohDWV\na3w/V6qCSu7/d61wECksCgmyIpOTMD2dvyGhvKSchoqG+XPmk+5YN+uq11Fk/P/nmqonQcMNIoVJ\nIUFWxL055GtIgLllkHk63JCL+QgA1WXV1JbXpuxJUEgQKSwKCbIifoWE8uJy6srrvDvoMtbV5Get\nhFyGBEi+DFIhQaQwKSTIivgVElqqWzDGeHfQZeRrQaVch4Rkz7koL3ceEa45CSKFRSFBVsRdH+9H\nSMiVfB1u6In35LYnQVUXRWSOQoKsiHtz8LJOQs5Dwtxwg7U2Z+fM1vTsNANjAzRXNefsnK01rXrI\nk4gACgmyQn4ON+RKS3ULY9NjxKby5053dvQskLvJneCEhJ54z5IS1nrIk0jhUUiQFYlCSMjHgkp9\no31AjkNCbSsziZn5gOJST4JI4VFIkBWJxWDVKigp8eZ4s4lZzo6ezflwA+RXQaVcLxOF5QsqKSSI\nFBaFBFkRr6stDo4PMmtn1ZNwCW5IaKpqytk53YJKyZZBarhBpLAoJMiK5HNJZldNeQ1VpVV515PQ\nWNFIaXFpzs7ZVNVESVFJ0qqL6kkQKSwKCbIiUQgJMLfCIc96EnJ9jYpMEetr1mu4QUQUEmRl/AoJ\nuVzaB/lXUKlvtC/nIQHmlkEmeRKkhhtECotCgqxIPO59SKgrr6OitMK7g65AvpVmDqInAZIXVNJw\ng0jhUUiQFYnF8ruQkivfqi4GFhKSFFRyhxvyqBaViGRJIUFWxI/hhsBCQp71JOR6SAZSDzckEjA2\nlvPmiEhAFBJkRSITEmrWcX7iPOPT4zk/d7rGp8cZmRwJbLhhZHKE2OSF8YXaWue7hhxECodCgqxI\nZELCXK0Ed+JkmAVRbdE1X1BpQW+CHhctUngUEmRFohIS3HO6N+AwC2qZKDiPi4aLqy66//+1wkGk\ncCgkyCVZ621ImJyZ5NzEuUBDQj70JAQZEtxHU6snQaSwZRQSjDH3GGOOG2PGjTFPGWOuv8T2dcaY\nvzHGnJnb50VjzBsza7Lk2sQEzM56FxKC7EZvrGyk2BTnTUgoNsU0Vjbm/NwVpRU0VDRc1JOgOQki\nhSftx/UYY94K3Af8NvATYDfwqDHmcmvtQJLtS4HvAL3AXcAZYANwPot2Sw7F4853r0JCkJ+Qi0wR\nTVVN9MXDP9zQF++jqaqJIhNMh9/iFQ4abhApPJk802838Flr7UMAxph3A28G7gb+Isn2/w1YDbzS\nWjs799qpDM4rAfH6MdFBhgSA5urmvOlJCOoawdKCShUVUFSkngSRQpLWR5S5XoGdwHfd16y1Fqen\n4MYUu/0i8CPgAWNMrzHmZ8aYe40J6OORpM29KXhVTKk33kuRKWJt5VpvDpimluqW/Ji4OBpwSKhp\nvehJkMao6qJIoUn3Rr0GKAYWv8P2AanezTqAX5071y8AHwPeD3wozXNLQPzoSVhbuZbiomJvDpim\nluqWvOlJaK7OfSElV6qqixpuECkcmQw3JGOAVMVai3BCxG/P9TrsN8a0Ah8APr7cQXfv3k1dXd1F\nr+3atYtdu3Zl32JZMT9CQpCfkJurmvn+ie8Hdv6V6o338roNrwvs/G21bfSN9jGTmKGkyHmr0JMg\nRcJl79697N2796LXhoeHPTt+uiFhAJgFFn+8aWJp74KrB5iaCwiug0CLMabEWjuT6mR79uxhx44d\naTZRvBa1kOAON1hrMcYE1o7lWGvpiwfzBEhXa20rCZugN947XzehpubCRFYRCV6yD8779u1j586d\nnhw/reEGa+008Axwi/uacd5lbwGeTLHbE8DmRa9dAfQsFxAkPPyYkxB0SBibHiM+Fd67XWwqxvjM\neOBzEuDigkrV1QoJIoUkk8mD9wO/bYx5pzFmK/B3QCXweQBjzEPGmD9bsP3fAo3GmE8bY7YYY94M\n3Av8dXZNl1yJxaCyEoo9mkIQdEhwH5gU5nkJQa8AAacnAS4uqKSQIFJY0p6TYK39ijFmDfBRnGGH\nZ4HbrbX9c5u0ATMLtj9tjLkN2AMcALrn/jvZckkJoXjcu6EGay19o8F2oy8szbylcUtg7ViOGxKC\nnLjYWNFIeXH5RT0JVVVw5kxgTRKRHMto4qK19gHggRQ/e32S134M3JTJuSR4XpZkjk/FGZseC7Yn\noVo9CSthjGF9zfqLlkGqJ0GksKhWgVxSLObtfAQI9uZXv6qe0qLSUIeEvngf5cXl1JXXXXpjHy0u\nqKSQIFJYFBLkkrzsSQhDSDDGOCscQlya2Z23EfTqi7baNoUEkQKmkCCXFI9HqycBwl+aOejJna7F\nBZUUEkQKi0KCXJKXExd7472h6EYPe2nm3tFgqy263Ic8uWVOqqpgdDTgRolIzigkyCV53ZMQhm70\nlqpwl2bui/fRUhWCnoTaVsamxxiedCq4VVfD9DRMTQXcMBHJCYUEuSQ/QkLQNNywMosLKrm/Bxpy\nECkMCglySaOjTjezF4J+sqFrYWnmsEnYROC1JFxuQSV3GaRCgkhhUUiQS4pkT0JVM1OzU5yfOB90\nU5YYGh9iJjETiuu0vmY9cKHqohsWNS9BpDAoJMglRTEkuG0I45BDGKotusqKy2iqatJwg0iBUkiQ\nZc3Owvi4NyEhYROBP9nQtbA0c9i49RvCcJ3gwgoHUEgQKTQKCbIst1vZi5AwODbIrJ0Nxc0vzKWZ\n53sSqoLvSYCLqy4qJIgUFoUEWZYbEryYuBiWQkoANWU1VJRUhLLqYm+8l5qyGqrKPJotmqWFBZUU\nEkQKi0KCLMu9GXjRkxCmkOCWZg5rT0IY5iO4Fg43VFQ4r2niokhhUEiQZfkREsLSjd5c3UzvaAhD\nQkiWibpaa1s5O3qWyZlJioqcXiX1JIgUBoUEWZbXIaGuvI6K0orsD+aBsD7kKSyTO11uQaWeeA+g\n5zeIFBKFBFmW1yEhTDe/5qpwVl3sjfeGoiSzq622Dbi4oJJCgkhhUEiQZXkaEkLWjR7mOQlhuk4b\nVm8A4MT5E4Ae8iRSSBQSZFler24I082vpbqFs6NnSdhE0E2ZN5OYYWBsIFQTF6vLqmmqauLYuWPO\nn9WTIFIwFBJkWfE4lJdDSUn2xwpbSGiuambWzjI4Nhh0U+b1j/ZjsaG6TgAd9R0KCSIFSCFBlhXF\nksyuMFbge4vaAAAgAElEQVRdDNMy0YUUEkQKk0KCLMurkDA5M8nQ+FCobn5hfH5DaEPCaoUEkUKk\nkCDL8ioknB09C4Tr5hfG0sxuW5qqmgJuycU66jvojnUzMTOhiYsiBUQhQZY1Ohq9ksyuytJKaspq\nQlUroTfeS0NFA2XFZUE35SId9R2As8JBPQkihUMhQZblVU9CGEMChG8ZZN9oX2gqUi7khoRj544p\nJIgUEIUEWZaXIaHIFLG2cm32B/NQ2Eoz98Z7WVezLuhmLLG+Zj1lxWUKCSIFRiFBluVlSFhbuZbi\nouLsD+ahsJVmDtsKEFdxUTEbV2+8KCRYG3SrRMRvCgmyLC9DQhhvfmErzRy2kswLddR3cPTcUaqq\nnIAwMRF0i0TEbwoJsizPJi6GrCSzK2xzEsIapuDCMkg3NGrIQST6FBJkWVHvSWipbmFgbICZxEzQ\nTWF8epzhyeFQXie4UFCpqsoZZ1BIEIk+hQRZVtRDQnNVMxZL/2h/0E2Zr/wYxusETkgYmx5jqsyp\neaGQIBJ9CgmyLC9CgrU2tCEhTKWZw7pM1OUugxyyTuVFFVQSiT6FBElpagqmp7MPCfGpOGPTY6G8\n+YWpNLPbhjAugQTYVL8JgIFZJySoJ0Ek+hQSJCX3JpBtSAjzJ2S3/HEYQkJPrIeSohIaKhqCbkpS\nteW1rKlcQ9+kQoJIoVBIkJTc7uRsVzeEOSSUl5RTv6o+FLUSeuO9NFc1U2TC+89yS8MWTsQPAQoJ\nIoUgvO9GErhC6EmAuaqLIehJCOu8jYW2NW/j5wMHKC5WSBApBAoJkpKXIaG8uJy68rrsG+WDluqW\nUJRmDmstiYW2N2/nxYEXqayd1MRFkQKgkCApeRkSWqpbMMZk3ygfhKU0cz70JGxv2c5MYoZVbS+o\nJ0GkACgkSEpeh4SwCktp5rBfJ4Brmq4BoGj9AYUEkQKgkCApeTZxMeTd6C3VLYHXSQhzLYmFaspr\n6KzvxDYpJIgUAoUEScm9CVRWZnecsN/8WqpbGBofYnJmMrA2nJ84z9TsVKivk2t7y3YmGxQSRAqB\nQoKkFI87vQhFWf6WhD0kNFc1A3B29GxgbQj7CpCFrm2+lrHqA8RH9axokajL6O3fGHOPMea4MWbc\nGPOUMeb6Fe73a8aYhDHmq5mcV3LLi5LMCZugL94X6ptfGEoz51NI2N6ynenSIQanuoNuioj4LO2Q\nYIx5K3Af8BHgOuAA8KgxZs0l9tsA/CXweAbtlAB4ERIGxwaZtbOhvvk1Vzs9CUFOXsyrkNC8HYD+\nogMBt0RE/JZJT8Ju4LPW2oestS8C7wbGgLtT7WCMKQK+CPwJcDyThkrueRES8uHmt7ZyLQYTeEio\nKq2iusyDR2767LK6yyhPrOZcuUKCSNSlFRKMMaXATuC77mvWWgt8B7hxmV0/Apy11v5jJo2UYIyO\nRrsks6u0uNR5JkGAtRLCPm9jIWMMTWwjXv1s0E0REZ+l25OwBigGFr+b9gFJ3+GMMa8CfgP4zbRb\nJ4HysifBnRwYVkGXZg77MtHFWou3M7FaIUEk6ko8Oo4Blkx1NsZUA18Afstaey7dg+7evZu6uotL\n+e7atYtdu3Zl2k5Jg1choa68jorSCm8a5ZOgSzPnU08CwBWVr+Ip+z/oGu6iva496OaIFKy9e/ey\nd+/ei14bHh727PjphoQBYBZY/LGwiaW9CwCdwAbgEXOhJm8RgDFmCrjCWptyjsKePXvYsWNHmk0U\nr8Tj0NSU3THy5ebXUt3CyfMnAzt/b7yXLZdtCez86bq+/nYejBXz8Iv/xj03vDvo5ogUrGQfnPft\n28fOnTs9OX5aww3W2mngGeAW97W5m/8twJNJdjkIXANcC2yf+3oY+N7cf3dl1GrJCU96EvKkGz3o\n0sz5EqZcLatXw8nX8PCL3wy6KSLio0yGG+4HHjTGPAP8BGe1QyXweQBjzEPAaWvth6y1U8ALC3c2\nxpzHme94MJuGi/+8mriYDze/IEszzyRm6B/tz4vr5KqpAQ7fwfc3f5ix6TEqS7MsyykioZT2Ekhr\n7VeA9wMfBfYD24DbrbX9c5u0kWISo+QXr+Yk5MPNr7mqmZHJEcamx3J+7v7Rfiw2L66Tyw0Jk7MT\nfO/494Jujoj4JKOKi9baB6y1G621FdbaG621P13ws9dba1PWTLDW/oa19q5Mziu5VUghYb7qYgDL\nIPNhmehitbXA4OW0VW7mm4c15CASVXp2gyRl7YVnN2RqcmaSofGhvLj5BVmaOR9DQk0NgOH6ujv4\n5uFv4pRLEZGoUUiQpMbGnKBQW5v5MdwHJuXDzS/I0szuOZuqslxKkkNOSIBryu6gO9bNj7t/HGyD\nRMQXCgmS1MiI8z2bkJBPn5AbKxopNsWBhYTGikbKistyfu5MuSGhPXEzW9ds5Y++90fqTRCJIIUE\nSarQQkJxUTFNVU2BzUnIh2u0UEkJVFTAWLyET936Kb53/Hs8cviRoJslIh5TSJCkvAoJRaaItZVr\nvWmUz4IqzZwvtSQWq6lxfk/etOVN3NpxKx/49geYmp0Kulki4iGFBEnKq5CwtnItxUXF3jTKZ0HV\nSsjHngRwQkIs5jzw6f7b7+fouaP8yX/+iYYdRCJEIUGSisWc7+7Ycyby7eYXVNXFfLtOrtraC78n\nVzddzcf/y8f55BOf5Fe+8isMT3hXO15EguPVA54kYtyehGxCQk+8J69ufi3VLTx+8vGcnzdfQ4I7\n3OC69zX3ctXaq3jX19/F1r/Zyh1b7uDWzlvZtHoTa6vWUlVaRZEpWvJVXFSMwaQ+kYikZWrGu2E/\nhQRJamQEysudr0ydiZ3hqrVXedcon7nDDdZaLjyPzF9j02OMTI7kbUhwexJcb9n6Fvb9zj4+8+PP\n8OjRR/nc/s8F0ziRQnbGu0MpJEhSIyPZzUcApyfhlk23XHrDkGiuamZseoz4VJya8iy6UNLgrqbI\nx5BQWwu9SUZnOuo7+Ks3/hXg/P3OxM7QP9bP+PQ4CZtY8jVrZ3PccpFoO3HwBH/893/sybEUEiSp\nbEOCtZaeWA/ra9Z71yifuTfq3nhvzkJCPi0TXaymBl56afltmqub5wtViUhu7JvZxx/jTUjQxEVJ\nKhbLbj7C4Pgg04lp1tWs865RPnMDTU+8J2fnzPeQsHBOgohEj0KCJJVtT0JPzLnRrqvOn5DgBpoz\nMQ8H9C6hN95LSVEJDRUNOTunV5LNSRCRaFFIkKSyDQnujTafhhtqymqoKq3KeUhormqmyOTfP8WF\nSyBFJJry751JciLrnoS5Lvt86kY3xrC+Zv18L0gu5OvyR3B6EuJxSCSCbomI+EUhQZLyYrihsaKR\n8pIs1lAGYF3NOs7Ec9iTkKclmeHCnJV4PNh2iIh/FBIkqVgs++GGfJq06Fpfsz7nww35GhLc3w8N\nOYhEl0KCJDUykn21xXyaj+BaX5374YZ8mty5kPv7oRUOItGlkCBJeTFxMR9vfrnsSbDW5nVPghsS\n1JMgEl0KCbLE9DSMj2c/cTEfexLW1awjNhUjPuX/QPv5ifNMzU7lbUjQcINI9CkkyBLum36mIcFa\nm9c9CUBOhhzyuZASaLhBpBAoJMgS7pt+piHh3MQ5pman8rInwW1zLoYc3HPke0hQT4JIdCkkyBLZ\n9iTMV1vMw9UNbu9HLkJCd6wbyK+CUwuVl0NpqUKCSJQpJMgSbk9Cpqsb3BtsPg431JTXUF1WnZPn\nN3SPdNNQ0UBFaYXv5/KLqi6KRJtCgiyR7XCDe4PNx54EyN0Kh+5YN601rb6fx096yJNItCkkyBLZ\nhoQzsTPUr6pnVckq7xqVQ+uq1+VsTkJrbf6HBPUkiESXQoIsMTICxkBVVWb798Tyc/mjSz0JK6fh\nBpFoU0iQJWIx5xNiUYa/HWfi+VmS2bW+Zn3O5iTkc5gCDTeIRJ1CgizhxcOd8vnml4uehNnELL3x\n3rzvSdBwg0i0KSTIEl48tyEfVza41lWvIz4VJzbp392vb7SPWTub93MSNNwgEm0KCbJENj0J+Vxt\n0TVfddHHIYfuEadGQhR6EjTcIBJdCgmyRDYhYWh8iImZCdpq27xtVA7louqiW0gp33sSNNwgEm0K\nCbJELJZ5SOga6QKgva7dwxblljvp0s+QcCZ2htKiUtZUrvHtHLmgkCASbQoJskQ2PQmnR04D5HVP\nQnVZNTVlNf72JMytbCgy+f1P0J2TYG3QLRERP+T3O5T4IpuJi13DXZQUldBc1exto3KstbZ1PvD4\noTuW/8sfwfk9SSRgbCzoloiIHxQSZIlsexLW16ynuKjY20blWHtt+/zQiR+6Y915Px8B9CRIkahT\nSJAlsgkJXSNdtNfm73wEV3ttu789CSP5X20RLvyeKCSIRJNCglzE2uwnLubzpEVXe107XcM+9yRE\nICS4PQlaBikSTQoJcpHxcZidzW64oa0mfyctutpr2+mN9zI1O+X5seNTcUYmRyIx3OD+nigkiEST\nQoJcJJsnQFprOT1yOjI9CRbrywqHqBRSAqivd74PDQXbDhHxR0YhwRhzjzHmuDFm3BjzlDHm+mW2\n/U1jzOPGmKG5r8eW216C5YaETFY3DI4P5n0hJZc7r8KPIQc3eERhdcPq1c6DwAYHg26JiPgh7ZBg\njHkrcB/wEeA64ADwqDEmVVWYm4EvAa8DXgl0Ad82xuRv3d4Iy6Ynwb2hRmHioht0/FjhEJVqiwDF\nxU5vgkKCSDRl0pOwG/istfYha+2LwLuBMeDuZBtba3/dWvt31trnrLWHgd+cO+8tmTZa/JNNSIhC\nISVXTXkNdeV1vvQkdI90s3rVaipLKz0/dhAaGxUSRKIqrZBgjCkFdgLfdV+z1lrgO8CNKzxMFVAK\naBQzhNylbBn1JIzMFVKqzu9CSq72On9qJZweOR2J+QiuxkYYGAi6FSLih3R7EtYAxUDfotf7gJYV\nHuOTQDdOsJCQcSeg1dWlv2/XcBetNa15X2rY5VdBpZPDJ9mweoPnxw2KehJEosurd3MDXLJ6uzHm\nD4H/CtxprfV+bZlkrb/fmYxWWpr+vqdj0VjZ4PKroNLJ4ZNsqFNIEJHwK0lz+wFgFljcn9zE0t6F\nixhjPgB8ELjFWvv8Sk62e/du6hZ9pN21axe7du1acYMlPQMDsHZtZvt2DXdFYj6Cq72una+9+DXP\nj3vy/EnedvXbPD9uUBQSRIKzd+9e9u7de9Frw8PDnh0/rZBgrZ02xjyDM+nwYQBjjJn782dS7WeM\n+X+ADwG3WWv3r/R8e/bsYceOHek0UbLU3w9rMnx68emR07yi9RXeNihA7bXt9I/1MzEzwaqSVZ4c\nc3himOHJ4UgNN6xZo5AgEpRkH5z37dvHzp07PTl+JsMN9wO/bYx5pzFmK/B3QCXweQBjzEPGmD9z\nNzbGfBD4GM7qh1PGmOa5r6qsWy+ey7QnYb6QUgSWP7rcoRMvhxxODp8EiNxww7lzTqVOEYmWtEOC\ntfYrwPuBjwL7gW3A7dba/rlN2rh4EuN7cFYz/AtwZsHX+zNvtvgl056E/rF+JmcnIzXcMF8rwcNl\nkCfPz4WECPUkNDY6j4s+fz7oloiI19KdkwCAtfYB4IEUP3v9oj9vyuQcEoz+/sx6EtxP21GauOhH\nQaWTwycpKy6jpXqli4HCr7HR+T44eOG/RSQaorFWTTwzMJBZT4L7aTtKPQmVpZU0VjR63pPQXtse\nmWWicHFIEJFoic47lWRtctKpuJhJT8Kxc8eoKKmguSoahZRcXhdUilqNBFBIEIkyhQSZ577JZ9KT\ncOzcMTrqO3AWu0SH1wWVolYjARQSRKJMIUHm9c9NPc2kJ+HouaN01Hd426AQaK9t59TwKc+Od/J8\n9EJCeTlUVyskiESRQoLMc+vvZzrc0Fnf6W2DQqCjvoPj547jPKIkOxMzE/SN9kVuuAFUUEkkqhQS\nZJ7bk5DucMNsYpbj549Hsieho76D0elRzo6ezfpYbo9E1HoSQA95EokqhQSZ199/oes4HWdiZ5ia\nnaKzIXo9Ce7f6di5Y1kfK4o1ElzqSRCJJoUEmecuf0x37uHRc0cBItmTsGm1U+bD/Ttm4+TwSQwm\nUstEXQoJItGkkCDzMi2kdOzcMQyGjas3et6moNWU19BU1eRZT8L6mvWUFZd50LJwUUgQiSaFBJmX\naSGlo0NHaa1t9ewhSGHTWd/pWU9CFIcaQA95EokqhQSZl3FPwvlormxwddR3eNOTEMEaCS63J8GD\nRSAiEiIKCTIv04c7HR2KZo0EV2d9J0eHsu9JOHH+RKRDwuQkjI4G3RIR8ZJCgszL9DHRUa2R4Oqo\n76An3sPY9FjGxxifHqdruIvLGy/3sGXhoaqLItGkkCCA86jfTOYkDE8MMzg+GO2ehLllkMfPHc/4\nGEeGjmCxCgkiklcUEgSA4WGYnU2/J8Edq49ijQSXG4Cymbx4ePAwAFesucKTNoWNQoJINCkkCJB5\ntcUo10hwrateR0VJRVaTFw8PHqZ+VT2NFY0etiw8FBJEokkhQYDMn9tw7NwxaspqInvzAzDG0FHf\nkdXkxUODh7i88fLIPSXTVVMDpaUKCSJRo5AgQBY9CUNH6WzojOzNz9VR38Gx89n1JER1qAGcKp0q\nqCQSPQoJAlwICY1pdggcHDgY2cl4C2W7DPLQ4CEub4j2ddJDnkSiRyFBAOfNvaEBSkpWvo+1luf6\nnmN783b/GhYSHfUdHD9/nIRNpL3v4NggQ+NDkQ9TbW1w6lTQrRARLykkCJBZIaWukS6GJ4fZ1rzN\nn0aFSGdDJ1OzU5weOZ32vlFf2eDq7ISj2decEpEQUUgQILNCSs/1PQdQED0JW9dsBeCF/hfS3vfQ\n4CEANjds9rRNYeOGBJVmFokOhQQBnG7i9evT2+e5vudYvWp1JB99vNjG1RupLqvmZ30/S3vfw4OH\naa9tp7K00oeWhUdnJ4yPQ09P0C0REa8oJAgAhw7BFWn2hj/X9xzbmrdFfmUDQJEp4uqmq/nZ2cxC\nQtSHGgA2z3WUaMhBJDoUEoSREefTX7oh4UDfAbY1RX8+guuapmvmh1jSUQgrGwA65uppKSSIRIdC\ngnDYmVeXVkgYnx7n8ODhgpi06NrWvI2DAweZnp1e8T4Jm+ClwZcKoiehosIZslJIEIkOhQThkDOv\nLq2Q8EL/CyRsoqBCwjVN1zA1O8VLQy+teJ9Tw6eYnJ2M/PJHl1Y4iESLQoJw6BCsWwe1tSvf57m+\n5zAYXtb0Mv8aFjLXNF8DkNaQgzvR8co1V/rSprBRSBCJFoUEyXjSYmdDJ9Vl1f40KoQaKhpYX7M+\nrRUOT595mrWVa7ms7jIfWxYeCgki0aKQIJmFhLPPFdRQg2tb87a0Vjg8feZprm+9viBWgIATEgYH\nnUePi0j+U0gocImEM3ExnZDglmMupJUNrnRWOFhrebr7aV6x/hU+tyo8Ojud7+pNEIkGhYQCd/q0\nUwAnnZBwcOAgA2MD3NB2g38NC6lrmq7h5PBJRiZHLrnt8fPHGRwf5PrW63PQsnBQSBCJFoWEApfJ\nyoZHjzxKeXE5r93wWn8aFWLu5MWfn/35Jbd9uvtpAK5fXzghoaEB6uoUEkSiQiGhwB06BGVlsHHj\nyvf51tFvcfPGmyNfZjiZK9dcSbEpXtGQw9NnnmZD3QbWVqX5UIw8ZoxTeVEhQSQaFBIK3KFDzpt6\ncfHKth+fHufxk49ze+ft/jYspMpLytnWvI0nup645LbupMVCoxUOItGhkFDg0l3Z8PjJx5mYmSjY\nkABwa8etPHb0MewyjzucTczyzJlnCmqowdXZCS+tvN6UiISYQkKBSzckPHr0UVprWrlq7VX+NSrk\nbuu8jb7RvmWXQh4cOMjo9GhBhoQbb3QmxB48GHRLRCRbCgkF7MwZ5xHRV1+98n0ePfoot3feXjDr\n/pN51WWvYlXJKr599Nspt3m6+2kMhp3rd+awZeFw661QVQVf+1rQLRGRbCkkFLB/+idYtQruuGNl\n23cNd/FC/wvcvrlwhxoAVpWs4uYNN/PYscdSbvOfJ/6TlzW9jNryNGpdR8SqVfDmN8NXvxp0S0Qk\nWwoJBcpaePBBuPNOZ8naSjx44EHKist4Q8cb/G1cHrit87b5+RmLjU6N8rUXv8avXvWrAbQsHO66\nC555Bk6eDLolIpINhYSQ27t3ry/H3b8fnn8e3vnOlW0/ND7Ep578FO95+XtoqGjwpU3L8es6ZOrW\njluZmJngh6d+uORnDx96mPhUnHdse4cv5w7btUjmTW9yltb6OeSQD9chV3QtHLoO3ssoJBhj7jHG\nHDfGjBtjnjLGLDs7yxjzq8aYg3PbHzDG/EJmzS08fv3SP/ggtLQ448cr8ZdP/CXTiWnuffW9vrTn\nUsL2j//qpqtpqW5JOi/hC899gZvab6KjvsOXc4ftWiRTUwO33ebvkEM+XIdc0bVw6Dp4L+2QYIx5\nK3Af8BHgOuAA8KgxZk2K7W8EvgT8A3At8HXg68aYwp0eH7DpafjSl+Dtb4eSkktv3xvv5dM//jS/\nf8Pv01zd7H8D84Axhts6b+OrB7/K5Mzk/Ot98T6+ffTb/Pq2Xw+wdeFw113wwx86E2RFJD9l0pOw\nG/istfYha+2LwLuBMeDuFNu/D/gPa+391tpD1tqPAPuA382oxZK1T30KBgZWNtQwMTPBe//tvZSX\nlPOBmz7gf+PyyPtvfD8nh0/yySc+Of/al5//MkWmqKDnI7je8haor4df+RWIxYJujYhkIq2QYIwp\nBXYC33Vfs05Fme8AN6bY7ca5ny/06DLbi08SCdi9Gz70IfjDP4Rtl3iIY9dwF6/5x9fwH0f+g3/4\nxX+gvqI+Nw3NE9uat/HBmz7IJ37wCQ72H+TQwCH++id/zZu2vInGysagmxe4hgZ49FFn7ssv/ZLz\nIDERyS8r6Gy+yBqgGOhb9HofkKokT0uK7VuWOc8qgPse+leav/d0mk0ES+pKeCva35LFEexF39Lf\n+8KO1sLPjpzg//7kA2nvC2ATMDkJ4xPQ1eVUwevtgf/6B9D2Svibb1y8/9TsFKPTo/TF+9jXs499\nPftYvWo1n7vtc3RMdLBv377M/lIeGB4eDvT8qdxRcwdfjH2RWz55C/2j/bRUt/BrV/yar20N67VI\npqgI9uyB3/1dWL/eKbS0fbuzoqaqyvk5OM98cLn/nawUx8KfnTw5zN//fX5cB7/pWjh0HRynTs1X\nMluV7bHMcqVll2xszDqgG7jRWvvjBa//BfBqa+1NSfaZBN5prf3ygtfeC3zYWrs+xXneBvzTihsm\nIiIii73dWvulbA6Qbk/CADALLJ691sTS3gJXb5rbgzMc8XbgBLB0IbqIiIiksgrYiHMvzUpaPQkA\nxpingB9ba98392cDnAI+Y639yyTb/zNQYa19y4LXngAOWGvfm03jRURExD/p9iQA3A88aIx5BvgJ\nzmqHSuDzAMaYh4DT1toPzW3/aeD7xpg/AP4N2IUz+fG3smu6iIiI+CntkGCt/cpcTYSP4gwjPAvc\nbq3tn9ukDZhZsP2PjDG7gE/Mfb0EvMVa+0K2jRcRERH/pD3cICIiIoVBz24QERGRpBQSREREJKnQ\nhYR0Hx4VRcaYe40xPzHGjBhj+owxXzPGXB50u4I2d10Sxpj7g25Lrhlj1htjvmCMGTDGjM09KG1H\n0O3KNWNMkTHmY8aYY3PX4Ygx5sNBt8tvxpjXGGMeNsZ0z/0b+KUk23zUGHNm7ro8ZozZHERb/bbc\ntTDGlBhjPmmMec4YE5/b5sG5Gj+RspLfiQXbfnZum99L9zyhCgnpPjwqwl4D/A/gBuANQCnwbWNM\nRaCtCtBcWPwtnN+JgmKMWQ08AUwCtwNXAu8HzgXZroD8IfA7wHuBrcAHgQ8aY6L+LJgqnEni95Ck\nnqsx5v/FeR7O7wCvAEZx3jvLctnIHFnuWlTiPEjwT3HuIb+MUw14UX3ZSFj2d8JljLkT53eiO5OT\nhGriYooaDF04NRj+ItDGBWguJJ0FXmut/WHQ7ck1Y0w18AzwHuCPgf3W2j8ItlW5Y4z57zhVTm8O\nui1BM8Y8AvRaa39rwWv/AoxZa1fwyLL8Z4xJAHdaax9e8NoZ4C+ttXvm/lyLU7DuXdbarwTTUv8l\nuxZJtnk58GNgg7X2dM4al0OproMxphX4Ec6Hi38H9lhrP5POsUPTk5Dhw6MKxWqcpDgUdEMC8jfA\nI9ba7wXdkID8IvBTY8xX5oaf9hljfjPoRgXkSeAWY8wWAGPMduBVOG+ABckYswnnWTgL3ztHcG6M\nhf7eCRfeP88H3ZBcmvuQ/RDwF9bag5faPpVMiin5JZOHR0Xe3P/ovwJ+WIi1JYwxv4bTffjyoNsS\noA6cXpT7cGqN3AB8xhgzYa39YqAty73/DtQCLxpjZnE+6PyRtfafg21WoFpwboLpPkgv8owx5Ti/\nM1+y1saDbk+O/SEwZa3962wOEqaQkIohm4cy5r8HgKtwPi0VFGNMG05AutVaOx10ewJUBPzEWvvH\nc38+YIx5GU5wKLSQ8FbgbcCvAS/gBMhPG2POWGu/EGjLwqeg3zuNMSXA/8a5BgX1CABjzE7g93Dm\nZWQlNMMNZPbwqEgzxvw18CbgddbanqDbE4CdwFrgGWPMtDFmGrgZeJ8xZmqul6UQ9ACLuwsPApcF\n0Jag/QXw59ba/22tfd5a+0/AHuDegNsVpF6cQKD3zjkLAkI7cFsB9iK8Gue9s2vBe+cG4H5jzLF0\nDhSakDD3SfEZ4Bb3tbmbwC0445AFZS4gvAX4L9baU0G3JyDfAa7B+bS4fe7rpzifnrfbMM269dcT\nLB1yuwI4GUBbglbJ0k/HCUL0XpZr1trjOEFh4XtnLc6wVCG+d7oBoQO4xVpbiKuAHgK2ceF9cztw\nBidk357OgcI23LDsw6MKhTHmAZwHYf0SMGqMcT8hDFtrC+bR2dbaUZwu5XnGmFFgMJuJOHloD/CE\nMelftvYAAAE7SURBVOZe4Cs4b/6/SWE+JO0R4I+MMV3A88AOnPeJzwXaKp8ZY6qAzTg9BgAdc5M2\nh6y1XTjDch82xhwBTgAfA04TwaV/y10LnBvhv+J8sLgDKF3w/jkUpWHLFfxOnFu0/TTOyqCX0jqR\ntTZUXzhjRyeAcZylGy8Puk0BXIMEztDL4q93Bt22oL+A7wH3B92OAP7ebwKeA8Zwbo53B92mgK5D\nFc6HieM4tQBewlkTXxJ023z+e9+c4n3hfy3Y5v/DuUmOAY8Cm4Nud66vBU6X+uKfuX9+bdBtz/Xv\nxKLtjwG/l+55QlUnQURERMKjYMfxREREZHkKCSIiIpKUQoKIiIgkpZAgIiIiSSkkiIiISFIKCSIi\nIpKUQoKIiIgkpZAgIiIiSSkkiIiISFIKCSIiIpKUQoKIiIgk9f8D8PFPNuWzX/wAAAAASUVORK5C\nYII=\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x7fe5f54e6850>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"k=2\n",
"\n",
"#Fails for t0 != 0\n",
"t0 = 0\n",
"t = linspace(0,14,200)+t0\n",
"\n",
"args = {'a':5,'b':0.2,'c':5+t0}\n",
"\n",
"H = [[sqrt(k)*tx,f(t,args)]]\n",
"rho0 = b[0]*b[0].dag()\n",
"c_ops = []\n",
"obs_ops = [b[0]*b[0].dag(), b[1]*b[1].dag()]\n",
"output = qutip.mesolve(H, rho0, t, c_ops, obs_ops)\n",
"\n",
"plot(t,f(t,args))\n",
"plot(t,output.expect[1])\n",
"ylim(0,1.2)"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": []
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 2",
"language": "python2",
"name": "python2"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 2
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython2",
"version": "2.7.12"
}
},
"nbformat": 4,
"nbformat_minor": 1
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment