Skip to content

Instantly share code, notes, and snippets.

@asafpm
Created June 8, 2021 13:57
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/b38e09fc0b026116ac90477129f341f6 to your computer and use it in GitHub Desktop.
Save asafpm/b38e09fc0b026116ac90477129f341f6 to your computer and use it in GitHub Desktop.
Eliminación adiabática (Adiabatic Elimination)
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"id": "07777a5a",
"metadata": {},
"outputs": [],
"source": [
"import qutip as qu\n",
"import numpy as np\n",
"import pylab as pl\n",
"\n",
"π = np.pi"
]
},
{
"cell_type": "code",
"execution_count": 2,
"id": "8c917f5d",
"metadata": {},
"outputs": [],
"source": [
"k1, k2, k3 = [qu.basis(3,j) for j in range(3)]\n",
"t = np.linspace(0, 60, 1000)"
]
},
{
"cell_type": "markdown",
"id": "374e3bf7",
"metadata": {},
"source": [
"## Dinámica Completa"
]
},
{
"cell_type": "code",
"execution_count": 3,
"id": "e5a1a446",
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.legend.Legend at 0x7f5e7afbffd0>"
]
},
"execution_count": 3,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXQAAAEGCAYAAAB1iW6ZAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/Z1A+gAAAACXBIWXMAAAsTAAALEwEAmpwYAABkI0lEQVR4nO2dd3hcxfW/31n13m314t6rbLlI7oABYwi9E1IoISH5JuRHEhKSkAYhJPRuOjbFmGJjMK4qltXcVVxk9d572zK/P6688rpg2Za0Wmne59GjvUd3935mdO/ZuefOnCOklCgUCoXC9tFZW4BCoVAo+gbl0BUKhWKIoBy6QqFQDBGUQ1coFIohgnLoCoVCMUSwt9aB/f39ZWRkpLUOr1AoFDbJ3r17a6SUAWf7m9UcemRkJBkZGdY6vEKhUNgkQojCc/1NhVwUCoViiKAcukKhUAwRlENXKBSKIYJy6AqFQjFEUA5doVAohgjndehCiLeEEFVCiMxz/F0IIZ4XQuQKIQ4JIWb1vUyFQqFQnI/ejNDfAVZ+z9+vBMZ2/9wHvHLpshQKhUJxoZzXoUspE4C679nlWuA9qZECeAshgvpK4OnUtHTyt03ZtHQa+usQCsW5qcqBxpIz7fWFUJc38HoUilPoixh6CFB8ynZJt+0MhBD3CSEyhBAZ1dXVF3Ww5BO1rEnK589fZlnYDxY38KN30kk8fnGfq1Ccl+I0eHk+vLIA2k4Z4zQUafaXYqDqiPX0KYY9feHQxVlsZ62aIaV8XUoZLaWMDgg468rV87J6ejB3xISz6VAZraeM0h/74jA7jlTx6PpDmEyqaIfiEqjJhZfmwTePWtrT3gAkdDTCoY977HvfBX0rGLsgY02PXUrY9GtYcwU0VwyIdMXwpi8cegkQdsp2KFDWB597TlZPD6bTYGLn0SoAiuvayCxtYlqoF2WNHewtqjfvW9nUwSOfHmRrdmV/SlIMJeKfguocSH0VKrM1m8kEJ7bDtFshaDoc+qRn/2PfQkQsjL8ajn6jOXKAwt2agy9Ogd3PDXw7FMOOvnDoXwF3d892mQc0SinL++Bzz0l0pC++bo7syNEc+nfdzvrJ66fhaKfju6ye0dA/vs5h/d4SfvXRftq6VNxdcRqNJWA85bwwGSF3K4y9AnQOcOBDzV5+ANpqYcxymLAKyvZBS5X2/spMGHc5jF8JjcVQ1f0lcHg9OLjB6GWQs7HH0YP2BdHZMmDNVAwPejNtcR2wBxgvhCgRQvxYCPGAEOKB7l02A3lALvAG8LN+U9uNnU6weFwAu45VYzJJvsuqYPxIDyYFezI7woek3FoAWjsNbOn+W2uXke3dXwAKBQAZb8H/JsP6H/bYyg9Cez1MvUlz3tlfao74xHZAaM557OXavse3aj+gfQGMXqa9PrFTc9jHvtU+Y/L1mqOv7J75KyV8fAc8GaaN6BWKPqI3s1xuk1IGSSkdpJShUso1UspXpZSvdv9dSikfklKOllJOlVIOSArFJeMDqGvtYufRKjIK67ls0kgAYsf6k1PeRHVzJzuPVtFpMPHn1ZPwd3cyj+RPcqSiidqWzoGQqxiM7HlJ+52zEaqPaq9P7NB+j1oCE1drjrj8ABzfpoVa3Py13+6BcHyL5pC9wiFgPHiFgt9YyNsJ5fuhuRzGX6U5dYC8Xdrvyiw4uhmkCRL/O4ANVgx1bHal6OJxAegE/GLdfowmyTXTgwGIHeMPQPKJGjYfLsff3YmYKD/ixvqTnFtjfmC6PaeSlc8msvrF3XTojVZrh8JKNJZAbS7E/h/o7OHAWs2etwsCp4J7AIy/EoSd9jC0OFVzzgBCaCGW7C81pz71Bs0G2ii9YDfse08L2Yy7AjyDwX98j0PPXK997tz7oSRdC90oFH2AzTp0b1dHlk0YQVuXkZnh3owP9ABgSogXXi4OfHO4gu05VVw9NRA7nWDhGH9qW7s4WtkMwKvxJwAobWhn5xF1QQ15Opost/Pitd9TbtRG41mfa/sUpfSETlx9ISquO44uYfJ1Pe+P/jEgNKc9654e+9QbwdAOe9/RQjOuvpp91BIoTAZ9B2R+BqOXwsw7tc/N3dYPDVYMR2zWoQM8ecM0fn/lBF68vSfbgJ1OcPmkkXybVUGnwcS1M7Up8QvH+AGwO7eG0oZ20gvq+dWKsfi7O7HpUL8+w1VYm6T/afHqnf/qseXHg6s/jJgEk66DhkLY/SyY9D0xcoDFvwNnL5hxpxZWOUnwDLg/Hh5IBN+oHnvYXG22i6M7LHqkxz5qCejbIPkFbd76lBtg5BRw8dFG9CeJ/zf8PRAOnjItUqHoJVarWNQX+Ls7cf/i0WfYH14+luzyJuZE+jIr3AeAIC8XRgW4kZRbg7E77HLdjBBK6tvZllOJ0SSx051tSr3CptG3w+7ntde7n4N5D2oOOi8eohaBTgcTroavfw2Jz4CLL4TF9Lw/Yj78rujsnx00/ez229ZqDz7FKedTZCzYOcHOv4OTF0y6Vjt2xEIoSNT2aauDhKe1+ew7/w7Tbrb8DIXiPNj0CP1chPm68vXDcfxl9WQLe+wYf1LyanlvTyEzw72J9Hcjdow/DW16ssu0W3IpJWtTi9h0qF+n0isGioIkaK+DZX/SQiGHP9UegLZUaKNm0MIisb8GoYMlvwc7h0s/7umO2NkT5j+kvV74MDi6aa8jY7W7g4ZiLRRj7NJi6w1FUHH40nUorIqUkpd25vJq/Amk7P8FjzY9Qr9QVk4J5L09hZQ2tPOby8cBsLD7IWpSbg1TQ73YllPFHz7XLqRQH1dmhHlbS66iLyjcrT30nPczLU5+cJ02agctjn2Spb+3dLT9wYo/a07dzb/HFhnbo/PAhzByKiz6LaS9pj1wDZrWf3oU/c63mRU8vUWbQTUmwJ0V3bPx+oshOUI/F/NH+fH4qkk8vHws183QYusBHk5MCPRgd24NAB+mFuLhZI+jnY4v9pdaU66iLyhMhuCZ4OgKM26H0r3aStDAaeAdbrlvfzrzk5zqzAFGTAa3EbDtL1C2X3tQ6h6gxfZPja0rbJLNmRX4ujkS6OnMm0la8rYT1S3msG9fM6wcuhCCH8VG8evLxqE7JV6+cIw/aQV1nKhuIf5YNfcujGThGD92HVWzX2yKzM/ggxt7luvr26F0H0Qs0Lan3gz2LtDVAtNvs57OU9HpIPpH2px1BzeYfotmj4zVpkoa9dp2TS589QsoGZBlHoo+QEpJal4tC8f4c9f8CFLy6sitaua6F3fzxMas83/ARTCsHPq5iB3jT5fBxK8+OgDAzXPCWDQugILaNkrq26wrTtE7Opvhq4e1Zfub/k+zlWRos1bCux26ewDcsxFWvwAx91tP6+ksfBiW/hHu+FSb9QLaw1J9mzZqB/j8fm1u+2c/tkwhoBi0FNa2UdXcSUyULz+YGYIQ8PO1+2nuNLBsYv+EXpRDB+ZG+eJor+NwaSOrpwcT6uPK3Cht/nBGQf153q0YFORu00beE1ZpybAqsyE/QXvQGX7KrJWwOTDrbtDZWU/r6Ti6weLfQuTCHltE9+uCJCg7AKUZEDQD6gu0uw7FoCc1X0tBMm+UL8HeLiyfMIIjFc34uzuxYLRfvxxTOXTAzcme/9w0nVXTgnh05QQAJgR64u5kT3pBT97r9XtLeCMhT6XnHYwU7gEHV7j6GbBzhL1va04+dE7PqNeWcA+AgAmaQ9/7Dtg7w83vag94j2yytjrFabR2GvjdZ4d4M7GnyElqXh1+bo6MDnAH4O/XTeXu+RG8cucsHOz6x/UOq1ku38fq6cGs7k4fANoCpVkRPuYR+t7Ceh759CAAPm6O3Dg71Co6FeegKBlCo8EjUFv9mfa6Zl/2R+vquhSiFvW0Y+Zd4BMJwbM0J68YVKxJyuejdK3Oz6wIH2aF+5CaX8fcKF9E9xTWQC9nnrh2Sr/qUCP072FupA9HK5tpbNPz1u58PJztGeHhxMaDao76oKKjUUt4dTJWHvdrbbTu4gsz77autkth9r1aagE7Jy3nDGgPS8v2qdS7g4ztOZVMCvLE29WB1+JPUFDTSmlDO/NG9U9o5VyoEfr3MCdSi6N/nFHEt5kV/GhhJCYJ76cU0tZlwNVRdZ9V6GiE5koI0NYSUJyuZS6MmK9t+4+FXx7S4uQnc6nYIiMnwYPJWjv8uldER8ZC0n+1GTAnszh2NmszZHRqfGYNWjoNZJY18bMlozGaJK/GnyDIywWAZRNGDKgWdQZ8D9GRvgR5OfPPzUcQwL0Lo1g0LoAug4l9hQ3Wljc8MerhzcvgpTmQuUGzFSVr2QtDonv2cw+wbWd+koBxPc4ctLQEOvuesEvWF/BkOHxyl1XkKbRwrNEkmRvly13zIxBC8E5yARMCPQjzdR1QLcqhfw92OsHPl41BJ+Cni0YR7O3CrHBvdALSTnlYqhhAjn8HNd25y7f9Ras2VJis5VVxcreqtAHByb0njm4ywtY/aXcnRzapAtVWIi2/FnudYHaED0FeLvxq+VgCPJz4/VUTB1yLcujn4Y6YCLKfWGme/eLh7MCkYE8ylEO3DvkJWnz8pne0HCj73rEMPwwHTsbRD3+q5Xy56j/a9MzMz6ytbFiSll/HlBAvcwj2F8vHkv7YChaPCxhwLcqh9wJnB8s5y9ERvuwvakBvNAGwv6iej9OLMHRvK/qRoj3abJaJq7VZH1//RhuhTlhlbWUDR2QsmAzaYiOfSO3haegcrVKSol9paOvib5uyzTUUOvRGDhY3EhM1OMJ7yqFfBHOjfGnXG8kqa6KqqYNbX0/h0c8O89bufGtLG9p0NmsZCMPnaw8Klz8OCBiz4typbIcikXHgN0Z7veQPYGev2Ur3aX2k6Df+891R1iTlc//7e6lo7GB/UQNdRpN5IaK1UdM0LoLoSG2hSnp+HTWtneiNJgI8nNiwr5T7Fp2Zn13RR5R0z2YJn6dtT7kBRi0FJ8/hlTfc3hF+sg2aymBkd4roqDhI/I+2wGrc5d//fsVFIaVkW3YV08O8ySxt5PWEPFwcdeiENoFiMKAc+kUwwsOZKH83Pt9fSlFdG1dNDWJGmDd//zqH4rq2AX+yPWwoStFixaFzemxDYSbLxeDiY7kCNixGWyGbH68cej9RUt9ORVMHDy0bwyh/N9alFWHfXd7Sy6UPcuj3ASrkcpFcNyOE7PImWrsM/N9l41gyXnsAknyixsrKhhAH1nZX8DFo20V7tALOTh7W1TUYcXCB0Lk91Y9MJjj0qcr70oek5msTIeZG+vLgktG06400dxoG1apxNUK/SO5fPAqjycSkYE9GB7gjpcTPzZHU/DpumRN+/g9QfD+le+GLB7XXDq5aitnidJj9Q6vKGtRELYJd/9JK2WVt0B4Y2znCrw5rKREUl0Rafi3erg6MHeGOTif45P75lDW0W6QMsTbKoV8kzg52/PrynqLBQgjmRPpaJPNSXAKZGzRnFDxTK5zsGayVkBtO0xMvlKg42PVPOLED4p8GtwBorYZDn2gpehWXRFp+HXMifc21FAbLg9BTUSGXPmROlC/Fde2UN7ZbW4rtU5QCIbPhyn9DRwN8+kOtOEXEwvO9c/gSMlvro89+rNVMveldLWPjie3WVmbzVDZ1UFDbxtxB8vDzXCiH3oecnIua1h1ra2jrYnduTb+VmxqydLVB+QFtemLwDJh6k2aPuV8rJac4O/ZOPQ9EJ6zS8quPXq7NfNF3WFebjdHYpuetpHyKarUCNyev6cE4Kj8V5dD7kIlBWg71tPw69EYTN7+2hzveTOW57cetLc22KN2rLZwJ7062dd2r8LMUWPEXq8qyCa59GW75AK7vTrsbsQCMnVB+0Lq6bIy/fZ3NE5uyuWNNCh16I3vyanFztGNysKe1pX0vyqH3IXbd+RzS8uv4JKOYY5Ut2OsEa1OL1Cj9QihKAYRWXQi0hTMjJg6vueYXi5M7TLymp+B1WHe1puIU62myMepau/jyQCkTAj0ormvntfg8duRUETc2APt+KkzRVwxudTbI3Chfjle18NjnmcyO8OGZm6dT09LJviJVyq7XFO3Rqt7bYqWhwYZ7APiOguI0ayuxGdIL6tAbJX+/bgpXTw3if9uOUdHUwYpJ/VMHtC9RDr2PuWFWKJ7O9jjYCf549USWjB+BEJCcW2ttabaByag5n5OrQRWXTliMlsBMFZfuFfsK63G00zElxIs/XD0RXzdHJgR6sGpakLWlnRc1bbGPCfRyZscjS2jvMppXjI4f6UFGoZrOeFZaa+G9a8HOAe76XMug2NXcEz9XXDphMXBwHdTlacm8Nv0fIOHq/2r9rrBgb2E9U0I8cXawI8TbheTfLcPRTmeerjiY6dUIXQixUghxVAiRK4T43Vn+7iWE2CiEOCiEyBJC3Nv3Um0Hf3cni+X/syN82F/UoOLoZyP9Dag8rKWD3fWklu8ctIyCir7BHEdPhfQ3Yd+7sO89OLzeuroGIZ0GI4dKG5kd0RPuc3awswlnDr1w6EIIO+Al4EpgEnCbEGLSabs9BGRLKacDS4BnhBCOfazVZomO9KGl08CRiiZrSxl85O3SKg3NugfSXoOEZ7Ql7J6D//bWZgiYAG4jIPkF2P4EjF4GXuGQ85W1lQ06ssqa6DKYLBy6LdGbEfpcIFdKmSel7AI+Aq49bR8JeAitvLU7UAcY+lSpDRMdoc1d3VuoHoxaoO/QpihGzIelj2lZEw3t2jJ/Rd+h08Hc+6AqW3tGsep/2nz1vF1g6LS2ukHFvu5rdFb40HXoIUDxKdsl3bZTeRGYCJQBh4FfSinPqPYghLhPCJEhhMiorq6+SMm2R6iPCyM9ncgo0E4Wg9FETnmTCsGUHwBjlxYv9xgJ9yfADzfD9FutrWzoseAX2qrbH27S4uijl4G+Dcr2W1uZValt6eTz/SU0d+gBbQFRmK8LIzydrazs4uiNQz9b8Oh0T3QFcAAIBmYALwohzpiBL6V8XUoZLaWMDggY+PJM1kIIQXSEr7ls3cMf7efK5xL581eZVlZmZYr2aL9Pxnh9IrTVjWq+ed/j4KyttA3tLqR9ss9P/g+GKQ9+sI//+/ggv1i3H73RxJ4TtcSO8be2rIumNw69BAg7ZTsUbSR+KvcCG6RGLpAPTOgbiUOD6Egfyho7eG9PAZsPV2CvE3ySUUJT98hgWFK4B/zHgZvtXkA2i5u/1vdFw3fBUU1LJ2kFdfi5ObLraDV//DyT5k4DsWNsd7DZG4eeDowVQkR1P+i8FTj9aUoRsBxACDESGA/k9aVQW2fFxJEIAY9/mcWYEe58+JMYugwmtmVXWluadTCZtNWLar659Qifpzl00/CshXsyP8vrd89mVrg3H2cU4+XiwNIJQ9ihSykNwM+BLUAO8ImUMksI8YAQ4oHu3f4GLBBCHAa2A49KKVWlh1MI83XlsasmMjPcm+duncGcSF+8XBzMJ9Wwo/oIdDSq+ebWJHy+lsmy5qi1lViFtPw6XBzsmBbqzbO3zOTm6FBevmMWro62uzynV8qllJuBzafZXj3ldRmg6l6dh5/EjeIncaPM2zPCvNlf1GA9QQOJyQSHPwGfKAiPgcLdml2N0K3Hyb4vStFy5TSWaJkuA8ZZV9cAkZpfx+wIHxzsdIT7ufLvG22/0Lha+m9FZoX7cKyq2fyEfUhz4AP4/H5463Kt8tDxrZpz94mytrLhi08UuI/UHHpNLrwUAy/NhfwEayvrdxrb9BypaBr06XAvFNu9txgCzAz3Rko4WNxI7Ngh/mAw9TXwCtNmsLx/HXS1aHOj1YwW6yGENkrPj9dSLpiMWpWo1Ne0cnZDmPSCOqTsqWEwVFAjdCsyI9wbIYbBgqOmMqjMhJgH4Ia3tPnPoOqDDgam3AjN5dr0xauehpl3aiXs9EO76lZqfi2O9jqmh3lbW0qfokboVsTT2YEJgZ5DP3HXyalxEfO1Mmk/3aFl/hs52bq6FFru9FX/A4TmzD0CIWON5uBHL7O2un4jLb+OGWHeODvYWVtKn6JG6FZmTqQP+wrrMRhNSClJOFZNY/sQi6kXpYCDKwRO07aDZ0LILOtqUmgIoaVaiL63JwQjdFCYbG1lfUpmaSOvxp+gqUNPY7uezLKmIRduATVCtzrRkb68t6eQzLImko5X85/vjjElxJONP49FDJX4ctEebYWiStU6+HHygKDpQ8qht3QauHNNKg1tetLz67hmejBGk2TphBHWltbnqBG6lVk01h9nBx0PfbiPZ7YeAyCztGnozE/vatXi56Fzra1E0VsiFkJJxpApLL0tu5KGNj2rpgWx/UgVv/n0IAEeTswI9ba2tD5HOXQr4+3qyC3RYZQ2tDMl2Iv0x1ZgrxPsOjZEkpeV7QdpgjDl0G2Gk4Wly/ZZW0mfkF5Qh4eTPc/eMoMVE0dgNEn+b8U4m8lxfiGokMsg4I+rJnHF5ECmh3nj5mTPrAgfEo5V8+jKIZAOpyRd+x06x7o6FL3n5Ordwt2ac7dx9hbWMzPCB3s7HW/cHU1TuwEv16EZ/lMj9EGAg52OBWP8cXPSvl/nj/Ijp7zJdhcc1eT2THsrTgff0eA69B5ADVlcfWHEZMs4uo3WI21s13O0spno7oIVQogh68xBOfRBSXSkDyaJbaYFSH8TXpwNay6DlmqtiEJUnLVVKS6UiAVQlKqlAvjgBngqUtu2MfYV1SMlZoc+1FEOfRAyM9wHncCcP92mSH0dHN2hMgtemQ/6Vpj8A2urUlwoEQu0/93bV0LuNi2J186/W1vVBZNRUIedTjAj3NvaUgYEFUMfhLg72TMxyJMMW1tB2lKtZe5b/mcwGWDnPyBoBkSqEbrNMXopOHtrVaXm3g/OnpD4DLRUgbvtTPfLKKhncrCnTWdQvBCGRyttkDmRvnycXozeaMLBzkZupE5Wv4mM1Wa1TLoOvEJBN7RW4w0LXHzgx1u1KaeTroWKQ5DwNOTFw7SbrK2uV3QZTBwsaeD2uRHWljJg2IinGH7MjvChXW8kp7zJ2lJ6T2Ey2Ltoo3LQ0rA6ulpVkuISCBgHU67XvpBHTgUHt55ZSzbAgeIGOvQm5kYNj/g5KIc+aImO1E7C9IJ6imrbuP/9DN5NLrCuqPNRuBvC5oC9o7WVKPoaO3stXUNJmrWVnBOjSXL/+xnEPrWDvYX1xB+rwk4nmD96iGcyPQUVchmkBHm5EObrwsaDZXyUVsTxqha2ZFUyN8qXiUFn1N+2Pu0NUHEYlvzO2koU/UVoNCS/oE1JdXCxtpozSDhWzZYsraTjXWtSaesyEhOlVQYbLqgR+iDm1jnhHChu4ER1Cy/fMQs7nWDTodPrcw8SilMBOSQWoijOQehc7WF32X5rKzkrn+4txt/dkV2PLMG724k/sHi0lVUNLGqEPoi5f9EofN0cifJ3Y94oP95NLmB7ThW/vWIQriAt3A06BwiJtrYSRX8RFqP9LkoZdF/cUkpS8upYOn4Ekf5ufPHQQmpaupgUPAjvZvsRNUIfxNjb6bhtbjjzRvkBsGhcAEcqmqlr7bKysrNQsFuLsaqHoEMXNz/wH9eT334QkVfTSl1rF3O6nz2N8HQeds4clEO3KU6udts3GOan69sh+ytoq4PWWi2RU9Ria6tS9Dfh86A4RSv6LSU0V1hbEQB7C7RrIjpyeKeYUA7dhpge5o29TgyOBUdf/Aw+uQveugIOrtUyKo6/0tqqFP1N5CLoaISCRNj4MDwzHva8ZG1VZBTW4ePqwOgAN2tLsSrKodsQzg52TA7xYq+1S9Y1V0L2F9pDsvpC+O6P4De2Z/65Yugy4Wpw8oL3VsO+9zRb4jNgtG4iub2F9cwK9xk6RWEuEuXQbYzoCB8OljTSaTBaT0TOV9qI/NoX4da1MOUGuHEN6NTpNORxdIUrnwTvCFj5FNz2MbTVaoWlrURjm54T1a3MGiYJuL4PNcvFxoiO8GFNUj5ZZU3MCrfSCVy4G7zCIGC89jN2hXV0KKzDjNu1HwBDJ9g5aiGYcVdYRc7+Yi0EOXOYJOD6PtSQysaY3f0U/+RDIKtQnKYqECk07J20ot/F1ltBmlFQj07A9CFYUu5CUSN0G2OEhzPhvq6k5tdx69wwnt12nFAfF364IHJg4oeNJdBU2jMnWaEImwupr2k1SB2c+/1w7V1GfvnRfpo69Dxz8wy2ZFUwJ9IXR50kPz+fjo6hUQvV2dmZ0NBQHBx6v9JVOXQb5PJJI3kzKZ/rX07meFULAFH+biwZPwBpTU/OQVYl5RQnCYvRUgKUH9CmNfYza9OK+C67Ep2AZf/ZRafBxN+unUxJSQkeHh5ERg7Q4KYfkVJSW1tLSUkJUVFRvX6fCrnYID+KjcLdyZ78mlaeuWk6vm6ObNhXOjAHL04DB1cInDowx1MMfk7erQ1Q2CXhWDVjR7jzxUML8XC2Z3qoFzfODqOjowM/Pz+bd+aglcrz8/O74LsNNUK3QYK9XdjxyGI69SbCfF3JKKxj48FyugwmHO37+Tu6OAVCZoPd8El4pDgP7iPAJ6o7n0//YjCa2FtYz7UzgpkW6k3qH1agE5id+FBw5ie5mLb06uoXQqwUQhwVQuQKIc6aTk8IsUQIcUAIkSWEiL9gJYoLYoSHM2G+2jL72DEBtHQa+j93emcLVGSq+LniTMJiNIfez8Wks8ubaOk0ENOdDsNOJ4aUE79UzuvQhRB2wEvAlcAk4DYhxKTT9vEGXgZWSyknA7ZR0mSIMCvCG9AK4vYL+u7bvuJUkMYBiZMqbIzwGGithtoT/XqYtHxtUV1M1PBe4n8uejNCnwvkSinzpJRdwEfAtaftczuwQUpZBCClrOpbmYrvI8jLhSAvZ/YVNfT9h3/9G/jHSNj4KzjwITh5aiXmFIpTGdO9FiH7c9j/AbwaC5mf9flhUvLqiPRzZaRn/8+msUV649BDgOJTtku6bacyDvARQuwSQuwVQtx9tg8SQtwnhMgQQmRUV1dfnGLFWZkV7tP3Sbvq8iD9TfAbA3vf1i7QKTcMyuIGCivjHQ4RC2HH3+HLh7RiJ189rBU+6SNMJkl6QR0xUX599pn9xcsvv8yUKVOIiIjghRdeGLDj9uah6NkCVKcHyuyB2cBywAXYI4RIkVIes3iTlK8DrwNER0f3b7BtmDE7woevD5dT1tBOsHcfOdyDH4HQwT0boewA5G6F5Y/3zWcrhh4r/6UlbQubq33xv3M1HN/aZ0Wlj1Y209iuJ2bU+cMtf92YRXZZ3z5TmhTsyZ+vmXze/T777DO2bt3K/v37qampYerUqTz44IPY2/f/HJTeHKEECDtlOxQ4vWxOCVAjpWwFWoUQCcB04BiKAWFud0wxvaCOa2ecfgN1keQnaqsAPYO1nwlX9c3nKoYmQdPhwd3aa5MJ3ALg2Dd95tBT82qBnnN9sPL888/zxhtv4ODgQFBQEA4ODuTl5fHUU0/R2NjI+vXrycvL4x//+Id5u6/ojUNPB8YKIaKAUuBWtJj5qXwJvCiEsAccgRjgf32mUnFeJgZ54u5k33cO3dAJpXth7k8v/bMUww+dTourH9+qzXzpg5kou0/UEuLtQqjP+Yuo9GYk3R/o9XoOHTrEuHHjACgvL8ff359x48axZs0abrzxRgBGjRplsd1XnDeGLqU0AD8HtgA5wCdSyiwhxANCiAe698kBvgUOAWnAm1LKzD5Vqvhe7HSCWRE+JOfWIqWkrrWLmpbOi//A0n1g7ITw+X0nUjG8iFgIbTVQffSi3i6lZHtOJXsL62jvMpJ4vJrlEwdgNfQlkJ2dTWNjI3l5eZhMJn7/+9/zi1/8YsCO36ugjpRyM7D5NNurp20/DTzdd9IUF8rVUwN59LPD/HNzDmtTizBKyfoHFjAlxOvCP6xoj/ZbOXTFxXJyNlRBIoy48Dq46/eW8Nv1hwAtzNKhN3HF5MC+VNjn7N+/nzvuuIPbbruN1tZWrr/+eu67774BO75a+j+EuGZ6MAEeTryRmE+Yrys6IXg1/iLnBRenaUUr3Ab/jALFIMUnEjxDtHTLF8G7ewoY5e/GNdODScuvY8FoPxaMHtzn44EDB1i1ahWpqalkZmbyxBNPAFBbW8sDDzzA/v37+de//nXGdl+hlv4PIVwd7dnw4AL2FdVzxeRA/vF1Dp9kFNPSacDd6QL+1VJqi4jGq4egiktACG2UfmLnBcfR82taySxt4o9XT+THsVH85rJxhPi4DPpVoQcOHODBBx88w+7n58err1oENc7Y7gvUCH2IEebryrUzQnB2sOOySSPpNJgufH56bS6016mc54pLJ2IhtFZBzfELeltK94yWZRNGIIQg0t8NB7vB76527drF+PHjrXb8wd9DiotmVoQPOgEZBRdYg7QgUfutcrYoLpWTcfTCpAt6W0ZBPX5ujkT5D++izxeKcuhDGHcneyYFe5LRmxF62QHY8zK010P2V1r2vADrjTQUQwTfUeARBAVJYDJCcbqW5O08ZBTWMTtCFX2+UFQMfYgTHeHLx+nF6I2mc9+yttZqq/q6WiD+KehogLhH+mTusGKYIwRExsHRzfDOKihKhpFT4L5d50zBXNXcQWFtG3fEhA+s1iGAGqEPceZE+tKuN37/Mui9b2vO/NqXtdFU0AyY/9CAaVQMcRb8XBudl+2HqTdDZSYc+fqcu2d018udEzm4V4QORtQIfYgT3V1UOr2gjulh3mff6di3EDoXZt6h/SgUfUnQdPjVIRB24OINeTsh+wuYfN1Zd08vqMPZQcfk4ItYPzHMUSP0Ic5IT62o9MlRzxl0tWkjp8iFAytMMbxwH6GtadDZwbiVkLtDy/dyFtIL6pgZ5tP/1beGIKrHhgHRkT6kFdRhMJ7lAirNAJMBwhcMvDDF8CRiAXQ2QvWRM/5U29JJdlkTc7rvLBUXhgq5DAMumziSDftK2X2ilvYuA6UNHdw1L0IbARUmA0LNOVcMHOai0ikwchIl9W38bVM240d64Opkj0nCVdOCrKvRRlEOfRiwbOIIfN0c+em7GXR1j9JrWjp5dOUEzaGPnKLFNhWKgcB3lJZatygVon/En77IZOfRarZkVQIwK9ybCYGeVhZ5abz88su8/PLLNDc388gjjwxYgi7l0IcBTvZ2/O+WGTy37RgrpwSyr7CBdWlF/HJJJM4l6TDzTmtLVAwnhDAXlS5taGfXsWoeXj6WOZE+xB+t5kexUZd+jG9+p1VN6ksCp8KVT553t8Fe4EIxBFg8LoDF4wIAiD9WzbdZFWRmxBOtb9NimgrFQBI+D45sYs+BTKSE62YEMyrAnbixAdZWdsmcrcDFhg0b2Lp1K1VVVTz00EO0tbXx9ddfm7cvv/zyPjm2cujDkJgoXxztdTRnb9cMkXHWFaQYfnSnBGg+spNAz9l9v8S/FyPp/uBcBS5uvvlmbr75Zurr63nkkUdYs2YN1113nXm7rxy6muUyDHF2sGN2uA8BVYla/NzN39qSFMONwGlIZy/8KpOJGeU7ZJb4n6/Axd///nceeuihc25fKsqhDydMJtjyGKy7nR+6JjHFkE3HhB9YW5ViOKKzoznycpaaUlga2AUJ/4Fj31lb1SVzaoGLadOmER4ezn333YeUkkcffZQrr7ySWbNmnbHdV6iQy3Di0Mew50Wwc+QK49c0SRcy/VahIugKa5Diex2Xi0+5btfJcIPQcrwEz7CiqkvjZIGLW265xcL+wgsvsG3bNhobG8nNzaWrq8ti+4EHHuiT4yuHPpzY/z74jYGf7qBj/yfcssnEigrBgmnWFqYYjmxpCGOPeJDHpzchZt4FH1wPGW/B6uetLe2iOVeBi4cffpiHH374DFtfo0Iuw4XWGq0U2NSbwdkL5/k/xT5wMukXmitdoegjMgrrKI26EXHdyxAxHyauhqwvwGiwtrSLRhW4UAwMJ+s6jl5qNkVH+nCguIFOg9FKohTDlbKGdgpr25gbdUpGxbGXaSkBKvt4/vgwQjn04UJhMti7aKlxu1k8LoAOvYmEYzXW06UYVuiNJqSUbD5cDsDyiSN7/nhyPUTBxRWVVqgY+vChcLeWr8Xe0WxaOMYfH1cH3kjII+l4NVllTfxp1aRzp9lVKC6B91MK+ctXWUwO9qSsoYNpoV6W8889g7VKWYW7tRzqigtGjdCHA+31UJGpFew9BQc7HQ8tHUNaQR3vpxRypKKZRz49iJTSSkIVQ5WmDj1PfXOE8SM9qGvtorXTwGNXTTxzx8iF2t3kOVLrKr4fNUIfDhTuAWRPwd5T+HFsFDPDfRjp6URafh2//uQgySdqWThGLTZS9B1fHyqnpdPAkzdMZXKwF3qjCWcHuzN3jIiF/R9AVZaWO0VxQagR+nCgcDfYOUHI7DP+JIRgdoQPoT6uXDU1CCd7HTuOVFlBpGIoszu3hkBPZ6aGeGGnE2d35tBTaEXF0S8K5dCHA/nxWvzcwfl7d3N2sCNmlB+7jiqHrug7pJSk5NUxrzdL/L3DtZ+CxIERN8RQDn2o0tUGXa1QdURLIzpuZa/eFjvGjxPVrVQ1dfSzQMVw4UR1KzUtncwb5de7N4xZAbnboeN7Cpsrzopy6EORymz470R4egx8dDvoHGDaLed/H5gvupR8teBI0TfsyasFYP7oXjr06beDoR2+fAheWwQf3Ghzzv3ll19mypQpRERE8MILLwzYcdVD0aHIzn9ovydcrcUir/4PuPcuz/SkIE88nOxJyatl9fTgfhSpGC6k5NUS5KUVK+8VodEw+XrI2qClqijfBrufg+V/6vUxn0p7iiN1Z9YsvRQm+E7g0bmPnnc/VeBC0Xe0VMPRb7R5vJc9ccFvt7fTMSfKl5TuUZVCcSmYTJLUvFrixgb0PkWuEHDjW3D537W56etug33vwpLfgZ1D/wruA85W4CIrK4tXXnmFmpoali9fzpIlS3juuefM22fL/3Ix9MqhCyFWAs8BdsCbUsqzZo8XQswBUoBbpJTr+0Sh4sLIjwdphEnXXvRHzBvly44jVVQ1dTDC8/sfpCoU30d2eRM1LV3EXug0WCHAK0R7Pf0WOPYNlO3vdTHz3oyk+4NzFbiYPn06r776KiaTiZ/+9Kc8+OCDFtt9xXlj6EIIO+Al4EpgEnCbEGLSOfZ7CtjSZ+oUF05BIjh5QuD0i/6Ik3H0PWqUrrhEdnZPgY0bdwnrGqIWAwLydvWJpv7k+wpcfPXVV8TGxrJ8+fKzbvcFvXkoOhfIlVLmSSm7gI+Asw3/fgF8Bqg5b9YkP1HLiWF38dG0ycFeBHs5s35vCZ0GI5mljRhNavWooneUN7bz/p4Cssua+Ci9mHmjfBnhcQl3eq6+2iKjvPi+E9lPnKvABcDq1atJTk7mww8/POt2X9Cbqz4EKD5luwSIOXUHIUQI8ANgGTDnXB8khLgPuA8gPDz8QrUqzkdTGdSdgOgfXdLH2OkEd86P4N/fHmXWE1tp7TKycnIgr9w5a8iUClP0Dx16Ize/tofiunaz7e8/mHLpHzxqMaS8qk3Hdezlw1UrcK4CF7t27WLDhg10dnZy1VVXnbHdV/TGoZ/tCj59uPYs8KiU0vh9F7yU8nXgdYDo6Gg15Otr8rsXY0RdetHnn8SOorali7rWLux1gk/3lnCguIGZ4T6X/NmKocvXh8oprmvnvzdPp7FdT6iPK0vHj7j0D45cBMkvQEkajFpy6Z/XT5yrwMWSJUtYsmTJGba+pjcOvQQIO2U7FCg7bZ9o4KNuZ+4PXCWEMEgpv+gLkYpekp8Azt4w8tJzYDja6/jTKu1RSXOHno2Hyvh0b4ly6Irv5fP9pUT4ufKDmSF9ezcXMR+EnTZoGcQOfdeuXVY9fm8cejowVggRBZQCtwK3n7qDlDLq5GshxDvAJuXMBxiTCXK3agUsdH27XszD2YGl40ew60gVUkoVdlGclfYuI2n5ddyzIKLvzxEnDwieqVICnIfzXvlSSgPwc7TZKznAJ1LKLCHEA0KIvqlsqrh0SjOgpbLXS/wvlAVj/Clr7KCgtq1fPl9h++wtrKfLaGJBf2XqHLMCitO0Z0WKs9KrqRBSys3A5tNsr55j3x9euixFr2hvgISnwcEVivaAsxeM77sHLKeysHvZ9u7cGsuiBApFN7tP1GCvE8yN9D3/zhfD1Jsg/klt1WhkLEgJE6/R5qx3M5TuIC+mLoFaKWrLbP4tHP5Eiy1KI1zxT3D27JdDRfm7EeTlTPKJGu6cF9Evx1DYNsm5NcwI88bNqZ/civ8YLSdR6qvaD2jn/PyHAHB2dqa2thY/Pz+bd+pSSmpra3F2vrDpnsqh2yoNxZC5Hub/HBY8DO11MOIsFWD6CCEEC8f4sy2nEqNJYqez7QtG0bc0tHVxqLSRXy4f278HuuY5iIzTVpEmPQuJ/9Wm6Tq4EBoaSklJCdXV1f2rYYBwdnYmNDT0gt6jHLqtcuxbkCbtZPYYqf30M7Fj/Fm/t4Tssiamhnr1+/EUtsOeE7VIyYUv8b9QHFxg1l3aa2mCD26AEzthwlU4ODgQFRX1/e8f4qj0ubbKiZ1aIQDfUQN2yJNl6RKOD40RkKLv2Hm0Cncn+4EtMB65SEtzcXTz+fcdJiiHbosY9dr0rdHLLB4I9TcBHk5ER/jw2d4SVUhagcFooqyhnaYOPZsOlXPV1EAc7AbQpdg7atN0c7drD0gVKuRik5Tuhc4mzaEPMHfOi+BXHx/gBy8nU1LfznUzgnns6ok2/xBKcWF0Gozc+MoeDpc24mSvo8to4u75kQMvZPQyyP4Sao5BwPiBP/4gQzl0W+TEDhA6iFo04IdePT2YfUX17DlRy+gAN95Mymfx+ADixvaugIZiaPD+nkIOlzby07goalu6WD5xJFNCrPBcZdRS7feJncqhoxy6bZK7DYJngcvAL8PX6QRPXKslW+o0GFn45E7e2V2gHPowY+OhcqaHefPY1Wdk0h5YfCLAd7Q2yJmn1jmqGLqt0VIFpftg3BXWVoKTvR3XzQgm8XgNbV0Ga8tRDBA1LZ0cKmlg+YQ+SLrVF4xeCgVJYOiythKroxy6rZH9JSD7bYn/hbJk/Ai6jCb2nFDFMIYLCceqkZK+yaLYF4xeBvpWlecFFXKxDVprIH2Nlug/+XkIidYS/g8C5kT54Opox66j1Syf2P9z4RXWZ+fRavzdnZgc3D+rki+Y0cvB1R92P6uFITsaIGpJnyepswWUQx/sSAkf3wVFydq2nSNc+/KATlf8Ppzs7Vgw2o/4Y2pu+nDAYDSRcKyayyaNRDdYVgs7OEPcb2DL7+GN7oek8x6Clf+0ri4roBz6YCdvp+bMV/1PywPt4AoegdZWZcGicQFsy6misLaVCD+VuGsoc6C4gcZ2/eAJt5xk3oPgPw66muHIZi3Xy9yfDOjCu8HA8LsnsTVyNoGDG8y4Qzs5B5kzB8wzXBLUKH3IszW7Egc7QezYfl7if6EIAWNXwOQfwGV/1dICHPrU2qoGHOXQBzNSwrEt2lN8eydrqzknkX6uRPi58k1mhbWlKPoRk0myObOcBaP98XJxsLacc+MZrBVKz/rc2koGHOXQBzNV2dBUMiimKH4fQghumBVK8olajlY009Shx2RSS7GHCuWN7aTl1/FxRjHFde1cPyvE2pLOz4RVUJ0DdfnWVjKgqBj6YObYt9rvsZdbV0cvuG1uOO8kF3DV84kYTZIJgR58dN88vF0drS1NcQlklzVxwyvJtOuNAEwK8mTVtGArq+oF467QHpIe2zKsFhypEfpg5tgWCJoxKOPmpxPg4cS7987l9rnh/GzJaI5VNvPKrhPWlqW4RJ7ecgRnBx2v3jmbf10/lbU/jbGNXPh+o7WHpCcHRcMENUIfrLTWaPUTF/8/ayvpNVNDvcx50ksb2vkwtYj/u2wczg52VlamuBgqmzrYebSaXy4fy8opg39QcQbjroCUV6Gjqd8qeQ021Ah9sJL9BdBdM9EGuWFWKC2dBjXzxYbZnlMFwFVTg6ys5CIZdyWY9NrU32GCcuiDEZMRMt6GgIkwcoq11VwU80f74eXiwJasSmtLUVwk23MqCfVxYdxId2tLuTjCYsDZGzI/s7aSAUOFXAYTic9A8gtaaty2Wrjx7UGzIvRCcbDTsXR8ADuPVqkapDZIe5eRpNwabpsbbru57u3sYc5PIPE/sPYW0LdB7K+1acBDFDVCHywUp8P2J7QR+aglsPJJbZGEDbN84kjqWrvYX1RvbSmKC2R3bg2dBhPLJw6yFaEXStxvYMoNUHNc+/nkHmgqt7aqfkON0AcL+94FRw+47SNwstFb3NNYPD4Ae51gW04V0ZG+1pajuAC2ZFXg4WRPTJSftaVcGo6ucONb2uuaXHhpDqS/Ccv/ZF1d/YQaoQ8GjAY48jWMXzlknDmAp7MDMaN82Zaj4ui2RKfByJasCi6bNBJH+yHkIvzHwJgVcGAtmEzWVtMvDKH/lg1TmATtdTDpWmsr6XOunBJEblULKXkqX/pgR0qJlJJPMkpo6jBw/axQa0vqe6beBM1lWl3eIYgKuQwGsr/SsiiOXm5tJX3OjbNDeX77ce5/fy+jA9wI8nLhnz+YipfrIM4FMgw5WtHMgx/spaS+HaOUzI3yZeEYGw+3nI2xl4POAXK+grA51lbT56gRurUxGSFno3aiObpaW02f4+xgx7s/msuMMG/s7XR8m1XB377OtrYsxSmYTJLfrj9IfVsXd82P4CexUbxyxyzbnd3yfbh4w6jF2jUnh16+ITVCtzZFe6C1CiattraSfmNikCfv/mguAP/cnMMbiXn8asVYQn2G3heYLbK/uIFDJY08ef1Ubp0bbm05/c/Ea2DjL6EyCwJtc53HuVAjdGtz8CNwdB80NUL7m7vnRwDw2d5SKytRnOTrQ+U42uu4epqNrgi9UMZfDQhtlD7E6JVDF0KsFEIcFULkCiF+d5a/3yGEONT9kyyEmN73UocYnS1QewIOf6rNN3ccHpV+Qn1ciYny5auDpcgheMtra5hMks2Hy1k8LgAP52HyXMM9ACJjYf/70NVmbTV9ynkduhDCDngJuBKYBNwmhJh02m75wGIp5TTgb8DrfS10yGDogrW3wr9C4IVZoLOHJWd8Rw5pVk0L5kR1K0cqmq0tZdizr6ieiqYOrrbVfC0Xy5LfQVMpPD8TXpgNSf+ztqI+oTcj9LlArpQyT0rZBXwEWMyvk1ImSylPLgdMAYbgfKc+Ys8LcOwbmPczWPoY/GgLeA2v7rpySiB2OsHGg2XWljLs2dQdbrH5FaEXSmQs3PQuhEaDiy9s+wtkf2ltVZdMbxx6CFB8ynZJt+1c/Bj45mx/EELcJ4TIEEJkVFcPwyx8UsK+9yEyDlb+S0uNO8QeyvQGP3cnFoz2Y9OhchV2sSJ6o4mvh1u45VQmXwe3fgj3fqPlTk942uZnvvTGoZ9t7tJZWy2EWIrm0B8929+llK9LKaOllNEBAQG9VzlUKE6D+nyYcbu1lVida6YHU1TXxsGSRmtLGbZsy66kurmTm6PDrC3FutjZw7wHoeIwlO23tppLojcOvQQ49T8eCpxxryyEmAa8CVwrpVTLAs/GwXVg72KzOc77kismB+LmaMffN2WzJauCg8UN1pY0LOjQG1mXVsQz3x3lz19lMcrfjaXjh+Hg6nQmXw92TtqsMxumNw49HRgrhIgSQjgCtwJfnbqDECIc2ADcJaU81vcyhwCGTsjaoDlzJw9rq7E6Xi4O/PXaKewvbuD+9/dy7Uu7eWH7cWvLGtJIKfnFuv38fsNhXtiRiwT+e8sM7O3U7GVcvGHCVZC5Xpu4YKOcd2GRlNIghPg5sAWwA96SUmYJIR7o/vurwOOAH/By9+oyg5Qyuv9k2yDHvoWORph+i7WVDBpunB3KvFG+1LZ08XpiHv/bdoxV04OJ8h8eUzgHmj15tWzNruS3V4znoaVjrC1n8DHtVsj6HHK3ac7dBunVV7OUcrOUcpyUcrSU8h/dtle7nTlSyp9IKX2klDO6f5QzP52974J7IEQtsbaSQUWojyvTw7z58zWTsNfpeHt3vrUlDVnWZ5Tg4WzPj2OjrC1lcDJmObj6a6FRG0Xdaw0ERalwYjvM/Yn2AEZxBiM8nLlyaiBf7C+lQ2+0tpwhR3OHns2Z5VwzPVgV7T4Xdg5aNsZj39psEQzl0PuLtjp4/3r4WwC8dTl4R8Ccn1pb1aDmptlhNHUY2Jqt8qf3NV8fKqdDb+Km2cNrzcMFM+cn2mK/N5bChzfBzn9qCfRsBOXQ+4uND0NBIkT/CBb/Du7ZqD14UZyT+aP9CPZyZv3eEmtLGXJ8ureE0QFuzAjztraUwY3/GLj9YwicCs3lEP8U7HrS2qp6jbr/7w/KD2mJf5Y+pi0eUvQKO53g+lmhvLwrl4rGDgK9nK0taUiQU97E3sJ6/nDVhKGZEreviVqk/QCs/7FWuH3ufVoOmEGOGqH3B/ve1ea0zlUhlgvlhtmhmCR8kFJobSlDAiklz207jrODTi0guhgW/RYM7XDoY2sr6RVqhN7X6Nu1DIqTVoOLj7XV2BxR/m5cOyOYl3flcrSyGT83R34SN4oxI4ZOrdWBYGt2JS/uzKWmuZPShnZ+e8V4vF0drS3L9hgxAULnaJkZ5z8Eg/wOR43Q+5rsL7X55rPutrYSm+Vv103h+lmh5FW38MWBUm5+bQ81LZ3WlmUzZJY28uAHe2lu1zMj3JvHV03iwcWjrS3Ldpl5J1QfgZIMays5L2qE3tdkvA2+o7UEXIqLwtPZgf/cpKXUP1rRzKoXEnlh+3H+eu3wS2R2Mfx36zG8XBzY8LMFalTeF0y5Ab77E6S+OujrkKoRel+Sux2KU2DOjwf9rZmtMD7Qg2umB/Pp3hIa2/XWljPoKaxtZefRKu6YF6GceV/h5KHdcWd9rhWlGcQoh36pNJXB+h/Bi3Pho9u10Xn0j62takjxo4VRtHUZ+SS9+Pw7D3Pe31OInRDcETMMaoMOJPMfAid3eO86+Pb32p24vsPaqs5AOfRLwWjQFh8c2QwB47Rv8R9uAgc13a4vmRLixdwoX97dU4DRZNv5qvuTlk4DH2cUs3JKICM91TnYp3gGw20fg0cg7H0HNv0KPr/f2qrOQDn0SyH7C6jMhB+8Ard8AFc9rf3jFX3OvQsiKalvZ1uOWkV6Lj5OL6a5w8BP4kZZW8rQJGI+/GQr/KEMlv1Ru/5zt1lblQXKoV8sUkLSs+A/HiZee97dFZfGZZNGEurjwlPfHqGisQOD0WRtSYMCk0mSW9XMjiOVPLftGPNG+arVoP2NELDgYfAKg/inra3GAjXL5WLJ3Q6Vh+Hal0Gnvhf7G3s7HU/dMI0fvp3GvH9tRyfg+lmh/Ov6qTgM03zeda1d/PDtNA51V33yd3fiyeunWVnVMMHeCeb/HL59FIpSIHyetRUByqFfHFJC4jPgGaJlZ1MMCAvH+PP1w3HEH60mr6aVdWlFjPBw4v+tnGBtaVbhme+OklPexN+unUy4n5anxctlGNYGtRaz7tJyvSQ9C7cPjkpHyqFfDEc2QVEyXPUfsFdTwwaScSM9GDdSq/jUoTeyJimfu+ZHEOTlYmVlA0thbSsfpxdze0w4d82PtLac4Ymjm5bjJf5JLX9TkPXvjobnverFUHMcvvw5vP8D+OwnWja2WfdYW9Ww5teXjcMkJS/syLW2lAHnv1uPYW8n+LmqPGRdYu4HtwBYezNs/BUk/Q/aG6wmRzn03lBfAGsuh8wN2j9r8g/g9k/V6NzKhPm6cnN0GJ9mFFPa0G5tOQNGTnkTXx0s496FUYxQ0xOti6sv3PGpFn7N+Qq2/QXeWgkdTVaRoxx6b9jyGBi74IFEuG8n/OBV8AyytioF8LPuEeqzW4dHbXKjSfLXjVl4ONlz/yI1PXFQEDwTfrod/l8e3LlBy/sS/5RVpKgY+vko3KPFzJf+EfxUgqPBRoi3Cz+KjeK1+DwMJsm4kR5ER/owJ9LX2tL6DKNJ8tm+EtLy6zhR3cL+ogaeumGqWto/GBmzXHtYmvqqVtxmgH2Gcujfh8kEW/4AHkEw/2fWVqM4B49cPp72LiOf7S3h8/2lAPxp1aQhUwz5j18cZl1aMQEeTvi5OfK3aydzyxy1tH/QsvSPWnh26+Nw64cDemjl0L+PjDVQtg9+8Jr2RFsxKHGw0/HEtVP46+rJNHUY+H/rD/LPzTksGO3HxCBPa8u7JLZmV7IurZj7F4/idytVxSGbwGMkxP4f7Pibtl5lzPIBO7SQ0jq5MaKjo2VGxiDML5z1BRz8SKsnWH4ARi+HO9arxUM2RENbF8ueiWeUvxuf3D8fnc42nWB7l5EV/43HzcmOrx+OG7YLqGwSfTu8thgai2H0Mi0HzIw7IGTWJX+0EGKvlDL6bH9TZ8ipHPoEPr0HqrK1p9dxv4Fb3lfO3MbwdnXkdysnkFFYz4dpRdaWc9E8v+M4pQ3t/O3aKcqZ2xoOLnD3FzD+Km3K88GPtJlyx7b062FVyOUktSdg0/9B+AK45yuwUyvubJmbokPZeKiMf36dg5eLAwtH++Hr5jjoQxbFdW0U1bWxt7CeV+NPcNPsUGJG+VlbluJi8AyGG9dor9sb4N1VsOE+eHA3eIX2yyFt/mu/Td92TvvZwklt+jaMJqOl0dBJ2/ofotfZww1vWDjzDkMHncYzy5/pjXo6DGfmQzZJE+2GM+dESykvSqtJnpmEqk3fhsFkOMPebminy9h1hr3L2HXWNhhMhovSeqFtOKO/u+1645kFKzqNnWfvb9O5+/tcmv7+g7GE+7ry8Lr9zP77Nmb9bStrkvIvqr/1pjO1fl9/n02r0WQ8Z3+3dLby+JeZxP17J3e8mcp/tx5j6fgR/P7qs09NHGz93dfn99n6u8PQcdb+Pte1+H393Zfn97muRYv+dvGGm96l06in8/P7tQkX/YBNO/RXDrxCzNoY1hxeY2H/6sRXxKyN4YmUJyzs6RXpxH4Uy4PbHtT+OfoOqMkl/7N7WOFQx82jxtLpHmDev6GjgWu+uIaVn62kuq3abNcb9dyx+Q4WfbyIo3VHLY7xaMKjzF87n8SSRAv7s/ueJWZtDJ8c/cTC/mHOh8SsjeF/+/5nYd9ZtJP56+bz2/jfWtiza7NZ8skS7v7mbosTqaK1gpWfrWT1F6tp7mo229v0bdzw1Q0s/WQpxU09BSJM0sR9W+9j4bqF7KvcZ3GMv6f8nXlr57E5b7OF/dWDrxKzNoY3Dr1hYd94YiMxa2P4656/WthP9vf92+63uBgKGgtYsX4FN228yeJCbOxsZPXnq7li/RVUtVWZ7XqTnrs238WijxdxpO6IxTF+l/g7FqxbQEJJgoX9+f3Ps2rjYu6+opQPfhzDn1ZNYkqIF0/uXkPM2hieyXjGYv/44ngWrFvAb3b9xsKeU5vDkk+WcNfmuyz6u6qtiis/u5JrPr+Gpq6eRSTthnZu2ngTyz5ZRmFTodlukibu33o/C9YtIKPC8tnRP1P/yfyP5rE26wt+uCCSdT+dR8JvlzJv1n4Wf7qA1w6+ZrH/13lfM2/tPB5PftzCvrdyL7EfxXLf1vss+ruoqYgV61dww8YbLBxcY2cj135xLZevv5zK1p60xHqTnrs3382ijxeRXZttcYzHkh5j/rr57CreZWF/8cCLxKyNYW3OWgv7x0c+JmZtDE9nWGYlTChJYMG6Bfxq568s7EfrjrLkkyXcuflOC6de3VbNlRuuZNXnq2jsbDTb2w3t3LzpZpZ9soyCxgKzXUrJA9seYMG6BaRXpFsc48m0J4lZG8NXJ76ysL95+E1i1sbwysFXLOzf5n/LvLXz+OPuP1rY91ftJ+6jOH763U8tvpiKm4pZ8ekKrv/qeosviCaPAK6NGsXlpnwqEvtnnrrNOvROYydvZ70NwJrMNeZvQyml+QLYcHyDhSN+49Ab6E169pTvIfOtZfCPQHhxNu9WpdJspyO3vZIdRTvM+68/vp6K1gpq2mv4+OjHZvv2ou3k1OXQbmjnnax3zPYjdUf4tuBbjNLI64deN9vrOur4MEebvvTaodfM/3yTNJmd44fZH1o44jcOv4FJmviu8DvyGvPM9rcy36Ld0M7hmsOklKeY7euOrKOuo47SllK+zvvabP86/2sKmgpo7mrmvez3zPbU8lTSK9LRm/S8efhNs72kuYT1x9cjkbx2qMeR6I16c1vfznzbYqR0sq0bjm+wcMRrDq9Bb9KTWp7KweqDZvv72e/T3NXMicYTbCvqySf92fHPKGsto7ajlo+O9CQ72lm0k6zaLNoN7byV+ZbZfrTuKN/kf4NRGi201nfU80H2B5qGzDdYMMaXH8dGseae2biP3KX1d85aC0f8+uHXMUoj24q2caKhp8zY25lv025oJ6s2i+SyZLP9oyMfUdtRS1lrGZtObDLbv8n/hrzGPJr1zbyf/b7Znl6RTmpFKgaTwaK/y1rKzF/y/qEJPL5qIvNH+xHk7cDbmdr5/XbW2xaj6NcPvY5E8kXuF1S0VpzR32kVaRyoPmC2v5f9Hs1dzeQ35rOtsKe/v8j9gtKWUuo66lh3ZJ3ZHl8cT2ZtJu2GdrMGgOP1x9mUtwmTNFl8yTR2Nprb+sbhN8x3CFJKXj+snRvrjqyzcMRvHHoDozSys3gnx+p7FoW9naX1d3ZtNrtLd/f099GPqGmvoby1nI0nNprt3+Z/S25DLs36Zt7Nftdsz6jMIKU8BYPJYDEAqWitMPf3awdfM3/x6U1687n1TuY7FgON1w69hkTy1YmvKG8pt+jvLlMXGZUZFoOiD3I+oKmriYKmAr4r/K6nv49/Qam+kTo7O97u7J/qWzbr0DMqMmg3tHPr+Ftp7mo2X2wFTQUUNRdx+4TbMUkTWwq0hxCt+lbSK9O5KXgxTlKyqasSFv0Wufol4v1DuSziMoLcgixOll3Fu5jiN4XokdFsL9putu8o3oGvsy+rR68mvjje/GXyXcF32Ak77p50NweqD5hHPTuLdtJp7OTeKfdS1VZldm5ZNVnUdtRy24Tb6DJ1mb9MatprOFxzmFvH34pO6MwjZb1Jz+7S3awevRpvJ28LrfHF8cwLmscE3wkW9p1FOwnzCGN5+HJ2FO8wf5nsKNqBi70Lt46/lT3le8xfJtuLtmOSJu6dfC95jXnk1mt5UtIr02nVt2r9rW82X2yFTYUUNBVwx8Q7kEi+zf8W0O4M0irSuHHcjTjbObMpT3N6Ukp2lexiefhyQtxDzPaTbZjkN4mYwBiLL9YdxTvwcfLhujHXkVCSYP4y2Vq4FZ3QcfekuzlUfcjs3HYV76LD2MG9k++lur3afLEdbzyKQTTi0r4Ig9Tz/J4NFNe1Ud5SzeFqrb/thJ35C9FgMpBUmsQ1o67Bx8nHwnHvLN5JTGAME30nntHfIe4hXBZxGTuKevp7Z/FOnOycuH3C7aSWp5JXW01WWSP/TvwUEyZ89JfTKsvJbdT6e2/VXlr0Ldwy/hZa9a0klSQB2ugvrzGP2yfcDmDR36nlqdww9gZc7F3MWqWUxJfEsyxsGaHuoRZadxTtYILvBOYFzbM4v3cW78TLyYvrx15PQkmC+cvkZH/fM+keMmszzc4tviSedkM79065l5r2GvZVaf19pO4IVW1V3DbhNgwmg/kY9R31HKw+yC3jb8FO2JnPb6PJSGJJIlePuhpfZ1+Lc2NX8S6iR0Yz2W8yG/M2WmgNdgvm8ojL2VG0w/xlsqNoB446R26fcDvpFenmL5PtRdsxSAP3TrmXouYijtZrd9gHqg7Q3NXMreNvpc3QRlKp1t+lLaXkNuSa+/ubgm8A7c4gpTyF68dej6u9q8X5HV8Sz5LQJUR4RlicM7tKdjHOZxxPL36an8Ra3s32FTbr0HcV78LZzplfzvol7g7u7CzeCWC+9b5n8j2M8R5jtu8p24PBZODKQ18z36Aj3j8YufQPHImYTU1XA4tDF7M8fDlpFWm0G9pp6GjgUPUh4kLjWBa+jNyGXIqbizGajCSXJRMbEstlEZfRrG8mo1K7hU4qTWJ6wHSuH3s9oJ3oALvLdjPCdQT3Tb0Pe5292VnFl8SjEzp+Nv1njHQdab6NPRmuuWHcDcwImGG2H6g6QIu+hWVhy1gcupik0iSMJiOlLaWcaDxBXIim9XDNYeo66ugwdJBWkWa2V7VVkVObg5SSpNIk5gTO4apRV2EwGcwOOqk0iTHeY7hr0l1aP5fsMverk50Tv5z1SzwcPc7o7zsn3sk4n3Fme0p5CnqTnpWRK1kQrIVEpJQcqz9GVVuVub/Ty9Np07fR2NnIgeoDxIbEsjR8KScaT1DUVIRJmkguTWZByAIui7iMVn2rOWSRVJrENP9p3DDuBvM5cbK/A1wCuH/6/TjqHNlRrPV3QnECAsGbq3+PzujDB4e/Ie7fO1ny4qtIJMtDVzFr5Czz/+1A1QGa9c0sC1/GkrAlJJUmYTAZKG8pJ7chl7jQOJaHLyezNpOa9ho6jZ2kVqQSFxLH0rClVLdXk1WTZf6fzgmcQ8yI5Rikgctfe52rn09iS148jqYg1vzgNxZtiC+Ox1HnyC9n/RJPR8+eNpT29PcE3wnm/k6rSKPL1MUVkVewMHgh8SXxSCnJbcilorWCRaGLWBGxgvRKrb+bupo4WH3QrLWgqYD8xnxM0kRSaRILgrX+bjO0mUMWu0t3M8V/CjeOuxHAfOyk0iT8Xfx5YNoDONk5WZzfAsH90+4nxD2EnUU9+0sk1425juiR0eb+PlRziKauJpaGLWVp2FJ2l+5Gb9JT2VrJsfpjxIXGsSJiBdm12dS019Bl7CKlPMX8f6jrqONwzWHzMeYEzmHVqFUYpMF8niaVJhHpGckPJ/8QgbDobwedAw/PehhvJ29zG06+7/aJtzPRd6J5//SKdDqNnVwRcQWxIbHm8zu/MZ/SllKz39hbuZeWrhZaulrYX7mfuJA4VkauJMC1J7Tbl9ikQ5dSklCSwLzgebg7upsdhkmaSChJYIz3GILdg4kbEc2+inSav3yQxIQn8DCZmKFzIy76IUrbKslvzCexVHOeC0MWEhcSR6exk/SKdJLLkpFIYkNiWRS6CNBOhszaTBo7G4kNiWVu4FwcdA7sLt1NTXsNOXU5xIbEMsprFEFuQSSXJWMwGUgpSyE2JBZ3R3fmjJxjPoETShKYETADb2dvFoUuIrksmS5jF4mliYxwHcF4n/EsCl3E0fqjVLZWkliSiL3OnnnB84gLjaOpq4nDNYfNo7e40DgWhSxCIkkuSyajMoNOYydxoXHEhsQiECSWJlLUXERJSwmxIbFM85+Gp6MnSaVJtOnb2Fu5l4XBCwlwDWCC7wR2l+7WRtXFu4gJisHd0Z3Y4FiL/h7tNZpQj1AWhS5if9V+mrqaSChJwM3BjVkjZrEodBHlreWcaDhh2d+hcXSZukivSGdP2R5M0kRciNYGgMTSRLJqsqjvrCc2JJY5gXNw1DmSWJpIXUcdWbVZLAxZSJRnFCHuIewu243BZGBP2R4WBC/AzcGNOUFzzBdlfEk80wOmMy04lOsnXIabVx7/vH4CoyKKkAYPfru2lgmeMRyrP0ZFawWJpYnYC3vmBWn93axv5mD1QXMbYkNiiQuNAyC5LJm9FXtpN7Rb9HdSaRJFTUUUNRcR7jKL337YgDS6Mn1cBc/fPgkXz0JunryCsX7BTPSdaO7v+JJ45gbNxcPRg9iQWJJKk8z9HekZSZhnGHEhcRysPkhjZyMJJQm42rsSPTKauNA4Kts0J2ihNSROOx/LU0gpS8Eojdo5c8r5nVOXQ11HHXEhccwJnIOznTNJpUk0dDRwuOYwsSGxRHpFEuYRRnJZMkaT0dzfrg6uzA2cy67iXeZrdGrAVPxc/IgLiSO1IpUOQwcJJQn4OfsxyW8ScaFx5DbkUtZSRmJJInbCjvnB84kLjaNF38KBqgPm0XJciNavJ7Xuq9pHu6Gd2JBYFoYsRCd0JJUmUdJcQkFTAbEhsUz2n4yPkw+7y3bTYeggoyKDhSEL8XX2ZbLfZPOdfXxJPHMDe/o7sTQRo8lIfEk8EZ4RRHhGsCh0EQerD9LQ0UBCSQIu9i5EB2r9Xd1ezdH6o+b+jguJ0/pbGkgtT9XCP9Jg1t9f9MqhCyFWCiGOCiFyhRC/O8vfhRDi+e6/HxJCXPrs+e/haP1RylrLWBy6GIDFYYupbq8mpTyFfZX7WBQSC3teZnHyGxiQ7C7cQYKhngVOI3H46Q7ixqwGNIcaXxLPZL/J+Lv4MztwNi72LiSWJJJQmoCPkw+T/SYT4RlBmEcYiSWJ7CrehZ2wM5/A0SOjSSpNMjuNuNA4hBDEhsSSUp5CWkUazfpm8z9ycdhi8hvzSS1PJacux3wxLQpdRJuhjT1le0guSyYuRPuck39PLE0kviSe2SNn4+bgxvzg+dgJOxJKEkgoTSDUPZRIz0gm+k3E19nXrNXF3oXokdHmEzipNMk8yogLicNOp7Vld9lu9pTtQW/SExuqaV0YvJADVQc4WH2Q0pZSc3/HhcZR21FLSnkKGZUZFm0wSiO7S3eTWJLIguAFONg5mNueUJpAQkkCE30nMsJ1BLNGzMLV3pXE0kQSShLwcvJiqv9UwjzDiPSMJLE0kZ3FO9EJHQuDF+Ji78KcwDkklSYRXxxv0d8LgxeSVp5GRmUGTV1N5jYsDl1MYVMhqeWpZNVmmbUuDV9Ml6mDoMAiGmQmi8MW0aE38fZWLXvhs7u/5OvcHUzynY67ozvzguZhL+xJLEkksSSRYLdgRnmNYoLvBPxd/Eks0bQ62zkzJ3AO3k7ejPGaxJfHtvNUwgYA3tzihIeTA3GhC6k1HsbNOw+9qYu4bq2xIbEcrD7IweqDFDcXsyR0ibkNdR117CnbQ3pFuvn/cLK/T55/84LmWfR3Ymki8cXxjPcZz0i3kcwcMRNXe1fz/p6Onkz1n0qoRyhRXlHmc0MndCwMWYiTnZO5vxNKE5BI4kLizOdGWoXW3w2dDT3nd+hiSlpKSK1I5XDNYfOX8+KwxbQb2kkuS2Z32W5iQ2LRCZ35C/HktThjxAw8HT21/tbZm8+NQLdAxniPYbzPeAJcAsznt5OdE3MD5+Ll5MU0/2mW53doHDqhY0HIApJLk0kpT6HD2NHThpCFHKw+yIGqAxQ0FZjPjcWhi2nobGBP+R7Sy9Mtzm+TNJk1xQTF4Gjn2HN+d7dhjPcYgtyDmDFiBu4O7ub9PRw8mD5i+vkd3CVwXocuhLADXgKuBCYBtwkhJp2225XA2O6f+4BX6Ee2FGzBTtixLGwpNBQT29aOAP60/WEM0sAVyW/Dlt8zPWAmng7uPB0USo0OLlvwKLiPIMg9iDHeY3g/+30OVR9iRcQKAJzsnIgJjOGbgm/YXridpeFLsdPZAZrzS6tI46vcr4gOjMbLyQvQLsITjSdYc3gNIe4hjPcZb7a36lt5fPfjONs5szB4IYD55PhD4h8AuDzicgBigmJwsnPiiZQnaNW3cnmkZh/jPYYgtyBePfgqeY155v09HT2ZMWIGn+d+TnJpMisiViCEQCd02gijJJHvCr5jYfBCnO01JxUbGsvhmsOsO7KOCb4TCPUINWutaa/hyfQn8XbyZvbI2VqbQ7URxm/if4NO6FgWvsy8v07o+NPuP2EwGbgi8goApvlPw8vJi/9k/Ieq9iqz1pFuIxnvM54Psz9kf9V+c3872jkyL2ge3+R/w7aibSwN6+nv2JBY0svT+fLEl8weORsfZx+zvaCpgDcOv0GwWzCTfCeZtbYZ2vhj0h9xsnPqcSTdzs/c3939OjdwLs52zvwj9R+06Fu4ffI1fP2LWJaOnoKpy4eNRe9T2VFAWlYoVz+fyPr0Gib7TWfD8Q0klSaxNGwZnQYTnXpJTOACEkuT+LZgCzP85/JJWiWX/S+BzNwgStqOsav8c2RnELfMnMnGX8Ry1eil1HbU8mTak3g6ejIncI65DUZp5JH4R9AJHcsjtCXjC0MWYifseHz34+hNenMbpvpPxcfJh//u/S+VbZVm+wjXEUz0nci6nHXsq9pn7m8HOwfmB89nS8EWviv8jiVhS7DX2ZvP7/SKdD4//jkzR8zE19nX3N+FTYW8dvA1At0Cmew32ay13dDOY0mP4ahz7HGGYVp/P5b4mEV/zwmcg4u9C/9I+QfNXc3mcybKM4pQ91DePPwmx+qPmc8ZNwc3Zo+czZe5X5JQmsCKcO38PjlYSipN4pv8b5gfNB9XB1ez1qzaLD7M+ZAx3mOI8Iwwt62+s55/pP4DDwcP5gbONbfBJE38Jv43CIS5nxaELMBO2PHn3X+my9Rl1jTFfwq+zr48u+9ZylvLzXZ/F38m+03moyMfkVGRwWURl2n9rdP6+7vC7/i24FsWhy3GQde/61t6M0KfC+RKKfOklF3AR8DpVZGvBd6TGimAtxCiX/LLJqQ9z4eH3iS204jvczPh2Sn4fnYfMe3tVJk6GWMUTPQZB3esx/6uDVwWqU2B83D0MF/cAFePupqqdm1GxqpRq3rso6+msbORDmOHORYOsHr0ajqNnVS1V3F11NVm+5VRV2In7ChqLmLVqFXmhStxoXEEuwVT2VbJ8ojl5pMuzCOMyX6TqWqvYkbADMI8wwBwsXdhcehiqtqqGOk6kpjAGACEEFwZdSWVbZU46BzMF8LJNtS012CQBq4bc53ZvmrUKpr1zdR31rNqdE/brhl1DSZporSl1KLNl0dejq+zrzb1MXKl+aSbNWIW43zGUdVWxYLgBfi7+APg4+zDvKB5VLVVEeUVxSQ/zana6ey4POJyrb8dPMwXN8BVo64y9/fq0ast2tDU1US7oZ0bxt7Qo3X0NXSZuqhqq+KaUdeY7SujVmIv7CluLubqUVeb+zs2JJZQ91Aq2ypZFr7M3N/B7sFM9Z9KVXsV0wKmmS9yZ3tnloYtpbKtkgCXAGKCYhjh6czLd8zmxzOvR+fQiJ2w41fzbkInBH/dmE3q4QjqO+sxSAOvb/Znwp++ZeLj3/LprpG06lto6Kxn594w/vxVFu5O9vwy5laEkOgc6/j1/Fv55w+m4uHswGURl521v2cEzGCi70Qq2yqZHzTf3N9eTl7MC55HVXsVkZ6RTPWf2tPfkVp/uzm4sTRsqUW/nuzva0b39N81o64x9/ep5/c1o6/R4tVtlRbnxhWRV+Cgczjj/F4QvIAwjzBzf7s5aLmOAt0CmR4wnar2Kqb4TSHKS0uQ5mTnxNKwpVS1V+Hv4s+8YK0GpxCCq0ZdRWVbJXbCjiujruxpQ9TV1HXUYTBZnt/XjL6GNkMbdR11Fuf31aO067KkpcT8GmBFxAr8Xfy1/o5aiUP3OpNp/tOY5DeJqrYq5gXNY4TrCEAbLC0IXkBVexXhHuFMD9BG1Tqh44puf+Jq78ry8J4cLVdFXUV1ezUSaXF+rxq1iuau5jP6u9+QUn7vD3Aj8OYp23cBL562zyYg9pTt7UD0WT7rPiADyAgPD5cXQ+GRjfKB9xbI8vU/lHLTr6VMfV3KwhRZXHNUPpn6pDxSe8Ri/4aOBvl02tMyrTzNwt5h6JDP73teflfwnYXdZDLJNYfXyI+PfHzGsT879pl849Ab0mA0WNh3FO6Qz+19TrZ2tVrY91ful0+nPS2r26ot7CfqT8gnU5+U+Q35Fvaq1ir577R/ywNVByzsrV2t8r8Z/5XxxfEWdr1RL1858Ir8MvfLM7R+mP2hfDfzXWkymSzsm/M2y5f2vyQ7DZ0W9pSyFPlMxjOyoaPBwp5TmyOfSntKljaXWthLmkvO2d//Sf+PTC1LtbB3Gjrl8/uel9/mf2thN5lM8q3Db8mPcj46ow0bjm2Qrx98/Yz+3lm0Uz6799kz+vtA1YGz93eD1t95DXkW9uq2avnvtH/L/ZX7LeytXa3yfxn/s+jvnPJG+VZSrrxnwz/kz798Tb6447h8aedx+cquXPlafK78+cYX5IMbn5Hv78mX2WWN5vd9k/+NfHH/i2f0d1p5mnwm/cz+PlJ7RD6V9pQsaS6xsJc2l8qn0p6SObU5FvaGjgb5TPozck/ZHgt7p6FTvrj/RflN/jcWdpPJJN/JfEeuy1l3xrnxxfEv5GsHX5N6o97CHl8cf9b+PlR1SD6d9rSsaq2ysOc15MknU5+UJ+pPWNi/r7+f3fus3Fm008JuMBrkawdfkxuObZCn8/GRj+Xbh98+ow1b8rfIF/a9IDsMHRb2k/1d315vYT9ad1Q+lfaULG4qtrCXNZfJp9Keklk1WRb2xs5G+Uz6MzK5NNnC3mXoki/tf0luzttsYT/Z32tz1p6h9WIBMuQ5/PV5k3MJIW4CrpBS/qR7+y5grpTyF6fs8zXwLyllUvf2duD/SSn3nutzB21yLoVCoRjEXGpyrhIg7JTtUKDsIvZRKBQKRT/SG4eeDowVQkQJIRyBW4GvTtvnK+Du7tku84BGKWX56R+kUCgUiv7jvNkWpZQGIcTPgS2AHfCWlDJLCPFA999fBTYDVwG5QBtwb/9JVigUCsXZ6FX6XCnlZjSnfart1VNeS+ChvpWmUCgUigvBJleKKhQKheJMlENXKBSKIYJy6AqFQjFEUA5doVAohgjnXVjUbwcWohooPO+OZ8cfqOlDOdZEtWVwMlTaMlTaAaotJ4mQUp41/67VHPqlIITIONdKKVtDtWVwMlTaMlTaAaotvUGFXBQKhWKIoBy6QqFQDBFs1aG/fv5dbAbVlsHJUGnLUGkHqLacF5uMoSsUCoXiTGx1hK5QKBSK01AOXaFQKIYINufQz1ewejAjhHhLCFElhMg8xeYrhNgqhDje/dvHmhp7gxAiTAixUwiRI4TIEkL8sttui21xFkKkCSEOdrflr912m2sLaDWAhRD7hRCburdttR0FQojDQogDQoiMbputtsVbCLFeCHGk+5qZ319tsSmH3suC1YOZd4CVp9l+B2yXUo5FK91nC19SBuA3UsqJwDzgoe7/gy22pRNYJqWcDswAVnbn9LfFtgD8Esg5ZdtW2wGwVEo545T52rbalueAb6WUE4DpaP+f/mnLuWrTDcYfYD6w5ZTt3wO/t7auC2xDJJB5yvZRIKj7dRBw1NoaL6JNXwKX2XpbAFdgHxBji21BqxS2HVgGbOq22Vw7urUWAP6n2WyuLYAnkE/3BJT+botNjdCBEKD4lO2SbpstM1J2V3fq/j3CynouCCFEJDATSMVG29IdpjgAVAFbpZS22pZngf8HmE6x2WI7ACTwnRBirxDivm6bLbZlFFANvN0dCntTCOFGP7XF1hy6OItNzbu0EkIId+Az4FdSyiZr67lYpJRGKeUMtBHuXCHEFCtLumCEEKuAKvk9hdltjIVSyllo4dWHhBCLrC3oIrEHZgGvSClnAq30Y6jI1hz6UCxGXSmECALo/l1lZT29QgjhgObMP5RSbug222RbTiKlbAB2oT3nsLW2LARWCyEKgI+AZUKID7C9dgAgpSzr/l0FfA7MxTbbUgKUdN/1AaxHc/D90hZbc+i9KVhta3wF3NP9+h60ePSgRgghgDVAjpTyv6f8yRbbEiCE8O5+7QKsAI5gY22RUv5eShkqpYxEuy52SCnvxMbaASCEcBNCeJx8DVwOZGKDbZFSVgDFQojx3ablQDb91RZrPzS4iIcMVwHHgBPAY9bWc4Ha1wHlgB7tm/vHgB/ag6zj3b99ra2zF+2IRQt1HQIOdP9cZaNtmQbs725LJvB4t93m2nJKm5bQ81DU5tqBFnc+2P2TdfI6t8W2dOueAWR0n2NfAD791Ra19F+hUCiGCLYWclEoFArFOVAOXaFQKIYIyqErFArFEEE5dIVCoRgiKIeuUCgUQwTl0BVDBiGEX3d2vgNCiAohRGn36xYhxMvW1qdQ9Ddq2qJiSCKE+AvQIqX8j7W1KBQDhRqhK4Y8Qoglp+QHd+vOS5/enSzp2m77D4UQXwghNgoh8oUQPxdC/Lp7nxQhhG/3fruEEM8KIZKFEJlCiLnddt/u9x/q3n+a9VqsGK4oh64YbjyGtix+DrAUeLp7eTnAFOB2tLwh/wDapJZQaQ9w9ymf4SalXAD8DHir2/ZXYL+UchrwB+C9fm+JQnEa9tYWoFAMMJejJbF6pHvbGQjvfr1TStkMNAshGoGN3fbDaCkCTrIOQEqZIITw7M4FEwvc0G3f0R3P95JSNvZvcxSKHpRDVww3BHCDlPKohVGIGLTqRScxnbJtwvJaOf3Bk0SldlYMAlTIRTHc2AL8ojtjJEKImRfxGbd0vzcWaOwehScAd3TblwA10oZzxCtsEzVCVww3/oZW2edQt1MvAFZd4GfUCyGS0cqL/ajb9he0qjSHgDZ6UqMqFAOGmraoUFwAQohdwCNSygxra1EoTkeFXBQKhWKIoEboCoVCMURQI3SFQqEYIiiHrlAoFEME5dAVCoViiKAcukKhUAwRlENXKBSKIcL/B95+PKKYQldMAAAAAElFTkSuQmCC\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"Ω1 = 1\n",
"Ω2 = 1\n",
"Δ1 = 5\n",
"Δ2 = 5\n",
"Γ1 = 1\n",
"Γ2 = 1\n",
"H = Ω1/2*(k3*k1.dag()+k1*k3.dag())+\\\n",
" Ω2/2*(k3*k2.dag()+k2*k3.dag())+\\\n",
" Δ1*k1*k1.dag()+\\\n",
" Δ2*k2*k2.dag()\n",
"\n",
"rho0 = k1*k1.dag()\n",
"rslt = qu.mesolve(H,rho0,t,[],[k1*k1.dag(),k2*k2.dag(),k3*k3.dag()])\n",
"\n",
"\n",
"pl.plot(t,rslt.expect[0],label=r\"$\\rho_{11}$\")\n",
"pl.plot(t,rslt.expect[1],label=r\"$\\rho_{22}$\")\n",
"pl.plot(t,rslt.expect[2],label=r\"$\\rho_{33}$\")\n",
"pl.xlabel(\"Tiempo\")\n",
"pl.legend(loc=0)"
]
},
{
"cell_type": "markdown",
"id": "b11d169f",
"metadata": {},
"source": [
"## Eliminación adiabática"
]
},
{
"cell_type": "code",
"execution_count": 20,
"id": "7ae660ad",
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.legend.Legend at 0x7fdc082268e0>"
]
},
"execution_count": 20,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXQAAAEGCAYAAAB1iW6ZAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/Z1A+gAAAACXBIWXMAAAsTAAALEwEAmpwYAABDX0lEQVR4nO3dd1yV5f/H8dfFHjIEwQUooKCIG/feW7OstGyYZdoetndpZcOGlWa2d5lZWe69Ffdii4KKLEVkj+v3x03fH7k4wDncZ1zPx8NHnMM59/2+H8GH+9z3dV0fIaVEURRFsXx2egdQFEVRjEMVdEVRFCuhCrqiKIqVUAVdURTFSqiCriiKYiUc9NpxgwYNZPPmzfXavaIoikXas2dPppTS70rf062gN2/enOjoaL12ryiKYpGEECeu9j11yUVRFMVKqIKuKIpiJVRBVxRFsRKqoCuKolgJVdAVRVGsRJUFXQjxhRAiXQhx+CrfF0KID4UQCUKIg0KITsaPqSiKolTFkDP0r4Dh1/j+CKBlxb9pwPzax1IURVGqq8px6FLKTUKI5td4yTjgG6mtw7tDCOEthGgspTxjrJCVxZ3NZdnBM7g72ePm7IC/hzPNfd1p5uuGi6O9KXap2KryMjiXDJnxkJcOBeegpADsHcHeGTwagVcA+LYA9wZ6p1UUo0wsagqkVHqcWvHcZQVdCDEN7SyeoKCgGu0s7mwuH66Nv+x5eztBZBNPopr70C/Mjx6hvjjaq1sESjVICanRkLgWjm+CU3ugtNCw93oHQdMoaDEYwoapAq/owhgFXVzhuSt2zZBSLgQWAkRFRdWos8bodk0YGdmYwtIyLhaVkpZTSHJWPrFpF4hOPsd3O07w+ZbjeLk6MiKyEZO7NyOyqVdNdqXYinPJsOdrOLwYzp8EBDRuD1F3gX8E+IVrZ+Ou9cHRDcpKoLQAcs9CTgqkH4VTe+HENjiyRHt/cB/odAe0HgMOzjofoGIrjFHQU4HASo8DgNNG2O5V2dkJ3JwccHNywN/DhXYB3tC+CQCFJWVsistg+eE0/th/mp92p9AxyJv7+rdgcGt/hLjS3x/FJqXshq3vQ+w/gIDQATDgOWg5FNx8rv4+O3twdNEKvH8raDlEe15KOHNA296BH+G3qeDmC93vg673gIs6sVBMSxjSgq7iGvoyKWXkFb43CngAGAl0Az6UUnataptRUVHS1Gu55BSU8NueVL7ensyJrHzaBXjxxLBw+rS84ro2iq1Ij4G1r0Ls31pR7jxFK7ieTYy3j/JySFoPO+ZDwmpw9oLej2jF3dHFePtRLMrFolI+WZ9A7xYN6NmiZpflhBB7pJRRV/xeVQVdCPEj0B9oAJwFXgIcAaSUC4R2yvsR2kiYfGCKlLLKSl0XBf1fpWXl/L7vFB+uiyclu4ARkY14YXQETbxd62T/ipkoytUK+e5F4OgOvR6G7jPAuZ5p93t6H2yYA3HLoX5zGPY6hI8E9WnRZkgp+fvQGV5bdpSzF4p4ZHBLHhkcVqNt1aqgm0pdFvR/FZWWsWjzceati0cgeGJYOHf2bI6dnfrFsnqxy+Hvx+HCaehyN/R/Btx96zZD4jpY/jRkxkLEOBg1V908tQFncgp4cvFBNsdn0qaJJ7Oui6RjUP0ab08V9EuknsvnxT+OsC4mnT4tG/DOje1p6Kk+Blul4nxY/iTs+1a7wTnmQwjsol+eshLYNg/Wvw6u3lqeViP1y6OY1LKDp3l2ySFKyyVPDW/F5O7NsK/lCaQq6FcgpeTHXSm8tuwozo52vHdzBwaE++uWRzGBs0dh8RTIiIXej2pn5Q5OeqfSpB2G36fD2UPQ4wEY/ArY69aeQDGywpIyXlh6mF/3pNIh0Jv3b+5A8wbuRtm2KujXkJRxkft/2EdM2gVmDg3nvv6haiSMNTj6Byy5F5w94PpPIXSg3okuV1oMK5+F3Z9Bs14w4UvwaKh3KqWWTp8vYPp3eziYmsMDA1rw8OCWRp0Tc62CbvMzb0L86rFkRk9Gt2vC2ytjuf+HvRSWlOkdS6kpKWHT2/DL7dAoEqZvMc9iDtqnhVHvwPiF2jj2zwZonyoUi7XnRDZjP9pCUkYeC2/rzMxh4XU6wdHmCzqAq5M9H07swLMjW7H8cBqTF+3kfH6x3rGU6iotht/vhXWzoO1NcMcyyzjjbX8zTF2lLTXwxXBI2qh3IqUGVh89yy2f7aSeswNL7+/J0DaN6jyDKugVhBBM6xvKvEkdOZiaw40LtnP6fIHesRRDFefDT7fAwZ9hwPNw/ULLGu/duB3cvQY8G8N3N8DBX/VOpFTDT7tOcu+30bRq5MFvM3rSwt9DlxyqoF9idLsmfHVXF87kFHLD/G0kZ+bpHUmpSuEF+H4CJKyBMR9Avycsc4y3dyDctRICu8GSe2Dvt3onUgwwf0MiTy85RJ+WfvxwT3d86+m31IMq6FfQM7QBP9/bncKSMiYu3KGKujkrOA/fjIOUnXDDIuh8p96JasfVGyYv1pYh+PMBiP5S70TKNczfkMicFTGMbd+ERXdE4e6s70glVdCvok0TL364pztFpWVM+kwVdbNUlKudmacdgpu/g7YT9E5kHI6uMPFHbU2ZZY/Ars/0TqRcwYKN/1/M597U3ixWd9U/gRlr3diT7+/WztQnfbaD1HP5ekdS/lWcDz9M1EaH3PgVhI/QO5FxObpof6TCRsA/M+HAT3onUipZtDmJN5fHMKaimDuYQTEHVdCrFNFEK+oXi0q5/YtdZOep0S+6Ky2GnyfDia3azc/Wo/VOZBoOztofq+C+sPQ+bfkCRXdL9qYy6+9jjGzbiPfMqJiDKugGiWjiyed3dOHUuQKmfLWbvKJSvSPZLinhzwe1JhRjP7SeyyxX4+gCE3/QRsH8eickb9U7kU3bEJvOk4sP0jPUl/du7mBWxRxUQTdY12AfPrqlE4dSzzPj+72UlJXrHck2rX8dDv6kDU3sdLveaeqGswfc+pvWFenHSdpSBkqd259ynvu+30tYQw8+va0zzg7m1/JSFfRqGBLRkDeub8umuAxe/vMIei2bYLP2fgub3oKOt0HfmXqnqVvuvjD5N+0yzA83QV6W3olsSkp2PlO/2o1vPSe+uqsLHi6Oeke6IlXQq+nmLkFM7xfK9ztP8s32E3rHsR2J6+CvhyF0EIx+zzLHmdeWdxBM+hEunIGfb4XSIr0T2YS8olLu+SaakrJyvprSFX8P852wpgp6DTw5LJzBrRvy6rKjbI7P0DuO9ctOgl+ngF8ruOlrsDfPs6M6ERAF4+fDye3w50PaPQXFZMrLJY/8vJ+4s7l8dEsnQv1M3AylllRBrwE7O8H7EzvQ0r8e932/l8SMi3pHsl7FefDTZO3rid9r15NtXeQNWu/Tgz/B9o/0TmPV5q6OY/XRszw/KoK+YebfulIV9Bqq5+zAojuicLK3Y/q3e9TIF1OQEv64HzKOwY1fgk+w3onMR98noPVYWP0SJG/RO41V+ufQGT5an8DELoFM6dVc7zgGUQW9FgLqu/HhpI4kZlzkud8PqZukxrb1AzjyOwx6yXyXwNWLEDDuY/AJ0YYzXjitdyKrcjwzjycXH6RDoDevjou0mB4JqqDXUq8WDXhsSBhL95/mu50n9Y5jPY5vgrWvQJvrtWbOyuVcPLXZpMX5WlEvVZPejKGwpIwZ3+3BwV7w8a2dcHKwnDJpOUnN2H39WzAg3I/X/jrKgZTzesexfBcz4Ld7wCcUxs6zzREthvJvBeM+0hYnW/2i3mmswkt/HCEmLZf3bu5AU29XveNUiyroRmBnJ3jv5g74eThz3/d7VXOM2igvh6UzoOCcNu3d2bxHFZiFyOuh672wcz7ErdQ7jUVbvCeVn6NTeGBAC4vsMawKupF4uznxya2dOHuhkGfV9fSa2/4RJKyG4a9rLeQUwwx5FRq21f4Y5qbpncYiJWZc5IWlh+kR4sujQ8L0jlMjqqAbUftAbx4bGsY/h9JYvCdV7ziWJzVau27eeixETdU7jWVxdIEJX0BJASyZpn3SUQxWXFrOIz/tx8XRjvcndsDezjIv86mCbmT39g2le4gPL/95RK2hXh1FufDbVPBooq6b15RfGAx/E45vhG0f6J3GonywNo5Dp3J44/p2NPQ035mgVVEF3cjs7QRzb9L+wj/88361iJehVj0P505oy+G6euudxnJ1uh0irtMaZZ/aq3cai7DreDafbEjk5qhAhkfWfWNnY1IF3QSaeLvyxvXtOJByng/Xxusdx/zFrYI9X0Gvh6BZD73TWDYhtL6q9Rpq19NLCvVOZNYuFJbw6M/7CfJx48UxEXrHqTVV0E1kVLvGTOgcwMfrE9RQxmvJz9Z6Z/pHaNPZldpz9dbWis+IgfWz9U5j1l764whpFwp5/+YOuvcDNQZV0E3ohdER+Hu4MPPXAxSVlukdx/xICcse1Yr6+E+1pWEV42gxWGuYvW0enNypdxqztOpIGr/vO8UDA1rQMai+3nGMQhV0E/JydeSNG9oSn36RD9aoSy+XOfwbHF0KA57ROvIoxjV0FngHapdeilU/3Mpy8kt4fulhWjf25IGBLfSOYzSqoJvYgHB/buwcwIKNierSS2UXM7TmxwFdoKea2m8Szh7aei/ZidpwUOV/Xvv7KFl5xbw9oR2OZtZGrjYMOhIhxHAhRKwQIkEI8fQVvu8lhPhLCHFACHFECDHF+FEt1/Pq0svllj+pLY077mOwt/xrl2YruG/FLNIFqh9phQ2x6Szek8r0fiFENvXSO45RVVnQhRD2wMfACCACmCSEuPR28P3AUSlle6A/8K4QwsnIWS2WuvRyidjlcGSJtgSsX7jeaazf4Je0bkd/PWzzXY5yC0t4ZskhWvjX46FBLfWOY3SGnKF3BRKklElSymLgJ2DcJa+RgIfQ1pisB2QDaoHwSgaE+zOhcwALNyURk3ZB7zj6KbwAfz+ujWrp9YjeaWyDk7vWti8rHja/q3caXb2xPIazFwp5e0I7s2zyXFuGFPSmQEqlx6kVz1X2EdAaOA0cAh6WUl42o0YIMU0IES2EiM7IsL3Wbc+NbI2nqyPPLDlEebmNrvWy9hVt7e6x88BBfYirMy0GQ9ubYPNcSI/RO40udiZl8cPOk0ztHWw1o1ouZUhBv9Ic7Eur0TBgP9AE6AB8JITwvOxNUi6UUkZJKaP8/My/nZOx1Xd34vlRrdl38jw/7LLBtdNPbIfdi6DbdK03plK3hr2urV7510M2t9ZLcWk5zy09TEB9Vx4bYr2X+Qwp6KlAYKXHAWhn4pVNAZZITQJwHGhlnIjWZXzHpvRq4cucFTGkX7ChWXylRVoh8QqCgc/rncY21fPTinrKTtjzpd5p6tRnm5NISL/Ia+MicXWyvkst/zKkoO8GWgohgitudE4E/rzkNSeBQQBCiIZAOJBkzKDWQgjBrOvaUlRazivLjuodp+5smweZcTB6rlrjXE/tJ0FwP1jzss20rTuZlc+Ha+MZEdmIAa0sb43z6qiyoEspS4EHgJXAMeAXKeURIcR0IcT0ipe9BvQUQhwC1gJPSSkzTRXa0gU3cOfBAS34++AZ1sek6x3H9M6dgE3vQOsx0HKI3mlsmxDaDdKyYlj5rN5pTE5KyYt/HsbBTljFWi1VMWgcupTyHyllmJQyVEo5u+K5BVLKBRVfn5ZSDpVStpVSRkopvzNlaGswrV8ILfzr8fzSwxQUW/nY9BXPaIVk+Jt6J1EAfEOh92NaA+6kDXqnMakVh9PYEJvBY0PDaexlWe3kasJ6pkhZGGcHe2ZfF8mp8wXM35iodxzTiVsJsX9DvyfBK0DvNMq/ej0M9ZvDP09YbXPpi0WlvPzXESIae3JHj2Z6x6kTqqDrqFuIL2PbN2HBxkROZlnhWhslBVrBaBAO3e/XO41SmaMLjHhLu6+xc77eaUxi7qo40nOLmD0+Egcrmt5/LbZxlGbs2ZGtcbATvGqNN0i3vAfnT8Cod9SYc3MUNgzCRsCGOVZ3gzTubC5fb09mUtcgqx1zfiWqoOuskZcLDw1qyZpjZ1kfa0U3SLMSYcv7EDlBW09EMU/D34DyUq1jlJWQUvLqX0dxd7Jn5lDrHXN+Jaqgm4G7egUT0sCdV/86aj2Ld618FuydtCVcFfPlEwx9HtOWMj6+Se80RrHq6Fm2JGTy+NBwfNxt65OhKuhmwMnBjpfHtuF4Zh6fbzmud5zaS1gDcSug3xPg2VjvNEpVej0M3s3g75lQVqJ3mlopLClj1t9HCW/owa3dgvSOU+dUQTcTfcP8GBrRkHlrEziTU6B3nJorK4EVz0L9YG2Kv2L+HF21IaWZsRD9hd5pamXR5iRSsgt4aUyEzdwIrcz2jtiMvTA6gnIpeeMfC148KfpLrTAMm61aylmS8BEQ0h82vKG1BLRAZ3IK+Hh9IsPbNKJniwZ6x9GFKuhmJNDHjWl9Q/jzwGn2njynd5zqy8+GDa9rU8vDR+qdRqkOIbR1XgpzYONbeqepkTf+iaFMSp4b1VrvKLpRBd3MTO8Xip+HM7OWHUVKC1tid8ObWkEY/oZWIBTL0rANdLoDdn8GmZbViGV3cjZ/HjjN9L4hBPq46R1HN6qgmxl3ZweeGBrO3pPnWXbwjN5xDJceoy2N2/lOrTAolmnAc+DoZlHDGMvLtWGKjb1cmN4/VO84ulIF3Qzd0DmAiMaevLk8hsISCxjGKKU2TNGpnlYQFMtVzw/6ztRGKSWs1TuNQf48cJpDp3J4Ylg4bk623Z9WFXQzZG8neH5Ua06dL+CLrRYwjDF+FSSuhf5Pgbtt3oyyKt2ma+u8rHwOysy7k2RhSRlvr4ylTRNPrutwaSM126MKupnq2aIBg1s35JP1iWTkmnFj37IS7RfftwV0uUfvNIoxODhrE8IyjsHer/VOc01fbk3m1PkCnhvVGjs7dd9GFXQz9uzIVhSWlDF3dZzeUa5u7zda8+Ehr6n1WqxJq9HQrDesnw0F5/VOc0VZF4v4ZH0Cg1v70zNUfTIEVdDNWohfPW7r0Yyfd58kJu2C3nEuV5SrjVsO6qmNY1ashxAw/HVtKOrW9/VOc0Ufro0nv6SMp0eobpf/UgXdzD08qCUeLo7M/vuY3lEut+0jyMuAoa+pYYrWqHF7aHcT7JgPOaf0TvMfiRkX+X7nSSZ1DaSFv4feccyGKuhmztvNiQcHtmBzfCZb4s2oq19umtYnNOI6CIjSO41iKgOeA1kO61/XO8l/zFkeg4ujPY8MDtM7illRBd0C3NajGU29XXlzxTHKy81kstGGN6GsCAa9qHcSxZTqN4Ou0+DAD3DWPNbs35GUxaqjZ5nRP5QG9dTyEpWpgm4BnB3seXxoGIdPXWDZITOYbJQRq90MjZqq9adUrFufx8HZA9a8rHcSysslr/9zjMZeLtzVK1jvOGZHFXQLMa5DU1o18uCdlbEUl5brG2bNK9pswn5P6ptDqRtuPlpT6fiVcHyzrlH+Oniag6naJCJXJ3tds5gjVdAthL2d4OkRrTiZnc8PO0/oF+TENq3pc+9H1CQiW9LtXvBsCqtf1GYG66CkrJz3VsfRqpGHmkR0FaqgW5B+YX70CPHlw3UJ5Bbq0IhASlj1Ang0hu731f3+Ff04umo3SE/vhSO/6xLhtz2pJGflM3NouJpEdBWqoFsQIbSz9Oy8Yj7blFT3AY7+AaeitV9sJ9td0c5mtZ8I/m1g7StQWlynuy4sKePDtfF0CPRmUGv/Ot23JVEF3cK0D/RmVLvGfLb5OOm5hXW347ISWPsq+LWGDrfU3X4V82FnD0NegXPJsOfLOt31DztPcjqnkCeGhSPUnIerUgXdAs0cGk5JWTkfrq3DNav3/wDZidowRTt1M8pmtRgMzfvAxjlQWDezl/OLS/lkQwI9QnzpZaOdiAylCroFCm7gzqSuQfy4K4WkjIum32FpkdbFpmlnNcXf1gmhnaXnZ2kzSOvAl1uTybxYzMxh4XWyP0umCrqFemhQS5wd7Hh7Zazpdxb9JVxIhYEvqCn+ivaHvdVobaawifuP5hSU8OnGRAa18qdzs/om3Zc1UAXdQvl5OHN3nxCWH07jUGqO6XZUnAeb39U+Zof0N91+FMsy8Hkovghb5pp0N59tSuJCYSmPDVVT/A2hCroFu6dPMN5ujryzyoRn6bsWQl669guszs6Vf/m31ka97PoMLpw2yS4yLxbxxdbjjGrXmDZNvEyyD2tjUEEXQgwXQsQKIRKEEE9f5TX9hRD7hRBHhBAbjRtTuRIPF0dm9AtlY1wGO5OyjL+DwhzY8j60GAJB3Y2/fcWy9X8aystg09sm2fz8DYkUlpTx2BB1dm6oKgu6EMIe+BgYAUQAk4QQEZe8xhv4BBgrpWwD3Gj8qMqV3N6jOf4ezryzKhZp7Bl82z+GwvPa2bmiXKp+c60p+N5vINu48yLO5BTw7Y4T3NApgFC/ekbdtjUz5Ay9K5AgpUySUhYDPwHjLnnNLcASKeVJACllunFjKlfj6mTPgwNbsDv5HBvjMoy34bwsraC3HgtNOhhvu4p16TsT7Bxh/RtG3eyHaxOQUvLQoJZG3a61M6SgNwVSKj1OrXiusjCgvhBigxBijxDi9ittSAgxTQgRLYSIzsgwYvGxcTd3CSKgvqtxz9K3vq/dEB3wnHG2p1gnj0baOi+HfoWzR4yyyRNZefwancItXYMI9FEzkqvDkIJ+pTthl1YNB6AzMAoYBrwghLjswpeUcqGUMkpKGeXn51ftsMqVOTnY8ehgbXndFYfTar/B3DTtZle7m8FftfdSqtDrYXD2hHWzjLK599fE42AvuH9gC6Nsz5Y4GPCaVCCw0uMA4NLb2qlAppQyD8gTQmwC2gNm3N3YulzXsSnzNyby7uo4hrZphH1tFi/a9A6Ul0D/p4wXULFebj7Q60GtoKfshsAuNd5U3Nlclu4/xbS+Ifh7uFTrvSUlJaSmplJYWIdLYpiQi4sLAQEBODo6GvweQwr6bqClECIYOAVMRLtmXtkfwEdCCAfACegGvGdwCqXW7O0Ejw8JY8b3e1m67xQ3dA6o2YbOnYA9X0HHyeATYtSMihXrNgN2fqot3HXHXzUe4vruqljqOTkwvW/1G6ekpqbi4eFB8+bNLX69FyklWVlZpKamEhxseCOPKi+5SClLgQeAlcAx4Bcp5REhxHQhxPSK1xwDVgAHgV3AIinl4Roch1ILwyMb0bapF++tiat5E4yNb4Gwg76qeYVSDc71oM9MSN4MSRtqtImDqedZeeQsd/cJob67U7XfX1hYiK+vr8UXc9BWVvX19a32pw2DxqFLKf+RUoZJKUOllLMrnlsgpVxQ6TVvSykjpJSRUsr3q5VCMQohBI8PDSP1XAE/7z5Z/Q1kxmu9I7tMBS/VQECppqgp4BWorcpZg5vz76yKo76bI3f1bl7jCNZQzP9Vk2NRM0WtTL8wP7o292HeugQKisuq9+YNb4CDq9ZuTFGqy8FZm2x0ei/ELKvWW3cmZbEpLoMZ/UPxcDH8mrHyX6qgWxkhBDOHhZOeW8Q325MNf2PaYTj8G3SfDvXUCCSlhtpNBN+WsG62NovUAFJK3lkVi7+HM7f3aG7afFZOFXQr1DXYh35hfszfmGh4q7r1s8HZC3o+aNpwinWzd4ABz0LGMe0EwQAb4zLYnXyOBwe1xMXROtba/+STT4iMjKRZs2bMmzevzvZryCgXxQLNHBrOmI+2sGjzcR6tai2M1GiI/Ueb4u+qlihVainiOmg4F9a/Dm3Gg/3VL6H8e3YeUN+Vm6MCr/q66nrlryMcPW3cBhwRTTx5aUybKl/322+/sXr1avbt20dmZiZt27ZlxowZODiYvtyqM3Qr1TbAixGRjfh8y3HO5VXR/3Hda+DmC92m1004xbrZ2WknB+eOw77vrvnSFYfTOHzqAo8MDsPJwTrK0YcffsicOXNwdHSkcePGODo6kpSUxNSpU5kwYQLAZY+NRZ2hW7HHhoSx4kgaCzYm8szI1ld+0fGKYWZDZ4OzR53mU6xY2DAI6KqtxNh+EjhePkmorFzy7uo4Qv3cGd/RuKOqDDmTNoWSkhIOHjxIWJj2qfjMmTM0aNCAsLAwPv/88/8V8JCQkP88Nhbr+JOoXFHLhh6M79CUr7cnk37hCuNZpdRm93k01oYqKoqxCAGDXoALpyD6iyu+5I/9p0hIv8jjQ8NrN7PZjBw9epScnBySkpIoLy/nmWee4cEH6+6+lCroVu6RwWGUlkk+Wp9w+TcT1kDKDuj7BDi61n04xboF94XgflrHq6L/9r4tLi3n/TXxtGniyfA2jXQKaHz79u3j1ltvZdKkSbRr146goCCmTZtWZ/tXBd3KBfm6cVOXQH7cdZKU7Pz//4aU2rVz72bQ8Tb9AirWbdCLkJ8JO//bUPqX6BROZuczc2g4dlZydg6wf/9+Ro8ezc6dOzl8+DCvvvoqAFlZWUyfPp19+/bxxhtvXPbYWNQ1dBvw4MAWLN6Tyrx18bw1ob325LE/4cwBuG4BOFR/mrWiGCQgCsJHwtZ50OVucK1PYUkZ89bF07lZffqHW9ech/379zNjxozLnvf19WXBggX/ee7Sx8agztBtQGMvVyZ3a8Zve0+RlHFRm/CxbjY0CIN2N+kdT7F2A56DoguwTRuP/d2OE5y9UMQTw8Ktaqo+wIYNGwgPD9dt/6qg24j7BoTi7GDHe2vi4dBiyIzVJoDYWcdEDsWMNYqEyOthxwLysk/zyYZE+rRsQPcQX72TWR1V0G1Eg3rOTOnVnOUHTlK8dhY0agutL+0kqCgm0v9ZKC0kbvGrZOcV8/hQ/c5irZkq6DZkWp9QbnPZgtOFkzDwBW0CiKLUhQYtKIqcSJtTv3JTmKBDoLfeiayS+o22IV6OZTzm9Ad7yluyz7nmXWUUpSY+d7gJCTzr/pfeUayWKui2JPoLPIrP8qnDrby7Ol7vNIoNSc8tZF50ITt8xuId8zNkJeodySqpgm4rii7ClrkQ3I+uA8axJSGT7YlZeqdSbMQn6xMpLisn5LoXwd5JW3tfMTpV0G3Frk8hLwMGvsDk7s1o6OnMu6tikTXoLKMo1XHqfAE/7DzJjZ0DCGwWrK25f2gxnD2qdzSrowq6LSg4D1s/gLDhENgFF0d7HhzYkugT59gQl6F3OsXKfbhGu7z30KCW2hM9H9IWgls/W8dU1kkVdFuw/SMozNEmeFS4KSqQQB9XdZaumFRSxkUW703l1u5BNPGuWC/IzUdrpBKzDFL36BvQRFSDC8U08jJhx3yt0UDjdv972snBjkcGhfH4rwdYcTiNEW0b6xhSsVbvrYnHyd6O+/q3+O83us+AnQu09YRuX2r8HS9/GtIOGXebjdrCiDerfJlqcKGYzpb3oCRfm9hxies6NiXUz513V8dRVq7O0hXjOnbmAn8dOM2UXs3x83D+7zedPbRm5EnrIXmLPgFN5EoNLpYsWcI999zDuHHjWLVqFUuXLv3PY6ORUuryr3PnzlIxsZxTUr7mL+XvM676kmUHTstmTy2TS/am1GEwxRZM/Wq3jHxphTyfV3zlFxTnS/lOuJSLhkpZXl7r/R09erTW26it4uJi6e3t/b/Hp0+flpGRkf97nJ2dLe+6666rPr7UlY4JiJZXqavqDN2abXpHW4ir31NXfcmIyEZENPbkvdXxlJSV12E4xZrtO3mONcfOcm/fELzcrtJT1NFVW4s/ZQfEr67bgCZSVYOLWbNmcf/991/1cW2pgm6tso/D3q+h0+1Qv9lVX2ZnJ5g5LIyT2fn8Gp1ahwEVa/bOqlh83Z2Y0iv42i/seJu2Jv+616Dc8k8ortbgQkrJU089xYgRI+jUqdNlj41F3RS1VhvngJ2DdgZUhQHh/nQK8mbeuniu79QUF0e1AqNSc9sSM9makMXzo1rj7lxFiXFw0lb9/P1ebY3+NtfVSUZT+bfBxc033/yf5+fNm8eaNWvIyckhISGB4uLi/zyePt04DdpVQbdG6cfgwE/a0DDPqkevCCGYOSycWz7byfc7TzK1dxVnVYpyFVJK3lkZS2MvFyZ3v/onw/9oe6N28379bGg9xqKXdL5ag4uHHnqIhx566LLnjE1dcrFG62eDUz3o/ajBb+kZ2oCeob58sj6BvKJSE4ZTrNm6mHT2njzPgwNbGv5Jz85eO0vPjIODP5s2oImpBheKcZ3aC8f+gp4PaBM4qmHmsHCy8or5aluyabIpVq28XPL2ylia+7pxY1RA9d7ceiw0bq+t8VJabJqANkAVdGuzbha4+kD3+6r91k5B9RnUyp9PNyaSU1BignCKNVt26Awxabk8OiQMR/tqlhYhYOCLcP4k7PvGNAFtgCro1iR5KySu1S61uHjWaBOPDQ3jQmEpizYnGTmcYs1KysqZuyqWVo08GNOuSc020mIQBPWAjW9DSYFxA9oIgwq6EGK4ECJWCJEghHj6Gq/rIoQoE0JMMF5ExSBSakO/PBpD13tqvJk2TbwY1a4xX2w5TtbFIiMGVKzZ4j2pJGflM3NoOHZ2NWz8LITWSetiGuz6rEabkFa0LlFNjqXKgi6EsAc+BkYAEcAkIUTEVV43B1hZ7RRK7SWsgZPboe9MbcJGLTw6OIyCkjLmb1BNCJSqFZaU8cGaeDoFeTOotX/tNta8F4QO1Ea9FF6o1ltdXFzIysqyiqIupSQrKwsXF5dqvc+QYYtdgQQpZRKAEOInYBxw6WLGDwK/Aaq3WV0rL4e1r2oTNDreXuvNtfCvx/iOAXyz4wR39wmhkVf1fqgU2/LdjhOkXSjkvZs7IEQNz84rG/gCfDZAW1Su/9VnOV8qICCA1NRUMjKsY0loFxcXAgKqd3PZkILeFEip9DgV6Fb5BUKIpsB4YCDXKOhCiGnANICgoKBqBVWu4difkHYQrlugTdQwgkcGt+TPA6eYty6e2ePbGmWbivXJLSzh4/UJ9GnZgB6hvsbZaNNO0Gq0tuxz13sMHq3l6OhIcLBtz6Ew5Br6lf7kXvqZ5n3gKSll2bU2JKVcKKWMklJG+fn5GRhRuaayUm3ceYNwaHeT0TYb6OPGzV0C+Xl3CinZ+UbbrmJdPt9ynHP5JTwxzMhjrwc+D0W5sPV9427XyhlS0FOBwEqPA4DTl7wmCvhJCJEMTAA+EUJcZ4yAShUO/qxNyBj4vNFn2D04sCX2doL316iG0srlsvOKWbT5OMPbNKJdgLdxN+7fWptBunMh5KYZd9tWzJCCvhtoKYQIFkI4AROBPyu/QEoZLKVsLqVsDiwG7pNSLjV2WOUSpUWw4U1o3EGbMm1kDT1duL1HM37fl0pCeq7Rt69YtvkbEsgvLuXxoWGm2UH/p6G8BDa/a5rtW6EqC7qUshR4AG30yjHgFynlESHEdCGEcVaUUWpm7zeQcxIGvaAN+TKB6f1CcXW0573V6ixd+X9ncgr4evsJxncMoGVDD9PsxDcUOk6G6C+1CUdKlQwahy6l/EdKGSalDJVSzq54boGUcsEVXnunlHKxsYMqlyi6qK2o2KwXhA4y2W586zlzV+9g/j50hsOncky2H8WyzFuXgJSSRwa3NO2O+j4Jwg42zDHtfqyEmilqqbZ/BHkZMPgVk52d/+vuPiF4ujgwd3WcSfejWIbkzDx+2Z3CpK5BBPq4mXZnXk2hy1Q48ANkqk+JVVEF3RJdTIdt87Tr5oGmH/bv5erIvf1CWReTzp4T50y+P8W8vbMqFgd7wQMDW1T9YmPo/Rg4uML61+tmfxZMFXRLtKlirYtBL9XZLqf0ak6Dek7MWRFjFTPxlJo5kHKeZQfPcE+fEPw96mjCWT0/6D4DjiyBtEN1s08LpQq6pclKhOgvtNZyDUx8/bISNycHHh7Ukl3Hs1kXk15n+1XMh5SSN5fH4OPuxLS+IXW7854PgouXtpqoclWqoFuadbPA3kkb0lXHJnYNIriBO28uj6FUNZS2ORviMtielMVDA1vg4XKVxs+m4uoNvR6GuBWQsrtu921BVEG3JKf2ah87e9wPHo3qfPeO9nY8NTyc+PSLLN6jGkrbkrJyyZv/xNDM141buhnYWs7Yuk0Hdz9Y96o++7cAqqBbCilhzUvg5gs9jd+L0FDD2jSiU5A3c1fHkV+sWtXZiiV7U4k9m8sTw8JxctCpbDi5Q5+ZcHyTtrqochlV0C1F4lrtB7nvkzVuXmEMQgieHdma9NwiPt98XLccSt0pLClj7uo42gd4Mapt1U3HTSrqLqgfDKtegPJrLh1lk1RBtwTl5bD6ZW153KgpeqchqrkPw9o0ZMHGRDJVEwyr99W2ZM7kFPL0iNbGWR63NhycYPDLkH4U9n+vbxYzpAq6JTi8GM4e0taJdnDWOw0ATw5vRWFpOR+uVZM9rNm5vGI+Xp/AwFb+xlset7YixkFAV1g3W5sxrfyPKujmrqRAa17RqB1E3qB3mv8J9avHpK6B/LDzJEkZ6pfKWn28PoG8olKeGt5K7yj/TwgYNltrVbf9I73TmBVV0M3djk8gJ0X7AbYzr/9dDw8Kw8nBjrdXxuodRTGBlOx8vtl+ggmdAwhvZKIFuGoqsCtEXAdbP4ALZ/ROYzbMq0Io/5V7FjbPhfBRENxX7zSX8fNw5t6+oSw/nKaWBLBCc1bEYGcHjw4x0fK4tTX4JSgr0Rq8KIAq6OZt/WwoLYShr+md5Kru6RuMn4czs/8+qpYEsCLRydksO3iGe/uG0tirdk3HTcYnBLpOg33fQdphvdOYBVXQzVXaYdj3rfYD6xuqd5qrcnNy4Imh4ew9eZ4/D1zayEqxROXlkleXHaWRpwv39qvjKf7V1XemtiTA6hf1TmIWVEE3R1LCque0H9R+T+qdpkoTOgcQ2dSTN5fHUFCsxgZbuqX7T3EwNYenRoTj5mRIH3kduflovyOJa9VkI1RBN0/xqyBpA/R7Glzr652mSnZ2ghdHt+FMTiGfbkrUO45SC/nFpcxZEUP7AC/GtW+qdxzDdLkb6jfXJhuV2fbsZVXQzU1ZCax6HnxbaAv7W4iuwT6MateYBRsTOX2+QO84Sg0t2JjE2QtFvDgmAjs7nScRGcrBGYa8qk022vOl3ml0pQq6uYn+EjLjYOgssK/jFe1q6enhrSiX2ugIxfKcPl/Awk2JjGnfhM7NfPSOUz2tx0LzPtpqpPnZeqfRjSro5iQ/Gza8rg1RDBuud5pqC/RxY1qfEP7Yf1oNY7RAb62IQUp4ani43lGqTwgYMQeKcm16GKMq6OZk3SwovADD3zR5n1BTmdE/FH8PZ1796wjl5WoYo6XYd/IcS/ef5p4+IQTUN3GfUFNp2Ea7TBn9hc0OY1QF3Vyc3q/9IHadpv1gWih3ZweeGt6KA6k5/L7vlN5xFAOUlUte+vMI/h7OzOhvvkNkDdL/GXDxhhVPa6PFbIwq6OagvBz+eQLcG+jSicjYxndsSvtAb+asiCG3sETvOEoVftp9koOpOTw3qjXuzmY+TLEqbj4w8HlI3gxHl+qdps6pgm4ODv4Mqbu0ZUFdvfVOU2t2doJXxrYh42IR769RqzGas+y8Yt5eGUu3YB/Gtm+idxzj6HwnNGyrDWMsztc7TZ1SBV1vhRe0WW5No6D9LXqnMZoOgd5M7BLEV9uSiUm7oHcc5SreXhlDbmEpr46L1H+tc2Oxs9dukOakaIt32RBV0PW2cQ7kZcDIt81uNcXaenJYOJ4uDry49Iha58UMHUg5z0+7U5jSs7n5raZYW817actNb3kPsmxnspt1VRBLkx4DOxdAp9uhaSe90xhdfXcnnhreil3J2SzZq26QmpOycskLfxymQT1nHh7cUu84pjF0tjbp6J+ZNnODVBV0vZSXw7JHwakeDHpJ7zQmc1NUIB2DvHlj+TFyCtQNUnPx8+4U7UboyNZ4uFjWBDaDeTbWunwlroMjS/ROUydUQdfL/u/g5DZtaVx3M2ntZQJ2doLXxkWSnVfMu6tUIwxzkHWxiLdWxtA12IdxHazkRujVdJkKjTvAimegMEfvNCanCroeLmZod+CDekKHyXqnMbnIpl7c1r0Z3+04weFT1v9LZe5m/X2MvKJSZl1nRTdCr8bOHsa8r92nWjdL7zQmZ1BBF0IMF0LECiEShBCXDZQWQtwqhDhY8W+bEKK98aNakVXPQXGe9oNmZTdCr+axoeH4uDvz9JKDlJaV6x3HZm2Ky+D3faeY0S+UsIZWdiP0app0hC73wK7P4NRevdOYVJXVRAhhD3wMjAAigElCiIhLXnYc6CelbAe8Biw0dlCrkbheG3fe+1Hws8A1M2rIy9WRV8e14fCpC3y+5bjecWxSfnEpzy09RIifO/cNaKF3nLo18Dmo1xCWPQLl1rtmvyGnh12BBCllkpSyGPgJGFf5BVLKbVLKf1dj2gEEGDemlSgpgL8f01pn9Xlc7zR1bkRkI4ZGNGTu6jiSM/P0jmNzPlgTT0p2Aa+Pb4uLo73eceqWixcMfwPOHIAd8/VOYzKGFPSmQEqlx6kVz13NVGD5lb4hhJgmhIgWQkRnZGQYntJabHoHspNg9Hvg6KJ3mjonhODVcZE42dvxzJJDamx6HTp8KodFW44zsUsg3UOs9yb8NbUZr61ium6W1Y5NN6SgX+muyRV/E4UQA9AK+lNX+r6UcqGUMkpKGeXn52d4Smtwer82yaH9JAjpr3ca3TTycuGZka3ZnpTFL9EpVb9BqbXSsnKeWXKI+m5OPDOitd5x9COEdjJl7wR/PqQNHbYyhhT0VCCw0uMA4LJuwEKIdsAiYJyUMss48axEaTEsvQ/c/bSPfTZuYpdAugb7MOvvY6RfKNQ7jtVbtOU4h07l8PLYCLzcrHTMuaE8m8Cw2XBiC0R/rncaozOkoO8GWgohgoUQTsBE4M/KLxBCBAFLgNuklHHGj2nhNr0N6UdgzAcW0SPU1OzsBG9e35ai0nKeX3pYXXoxobizucxdFcfwNo0Y1bax3nHMQ8fJEDIAVr8E507oncaoqizoUspS4AFgJXAM+EVKeUQIMV0IMb3iZS8CvsAnQoj9QohokyW2NKf3w+Z3tUst4ZbXhchUQvzq8fiQMFYdPavWTTeRkrJyHv/lAPVcHJg13gbGnBtKCBj7ofbfvx62qmUBDBoELaX8R0oZJqUMlVLOrnhugZRyQcXXd0sp60spO1T8izJlaIuhLrVc0919QujSvD4v/XFENZY2gU/WJ3LoVA6zr4ukQT1nveOYF+8gGPIKJK2Hvd/oncZobGNWi142vaUutVyDvZ3g3Rs7UCYlTyw+oFrWGdHhUznMWxfPuA5NGKEutVxZ57u0xtIrn9VGn1kBVdBN5cT2ikstt6hLLdcQ5OvGC6Mj2JqQxbc7rOt6pl6KSsuY+esBfNydeGWs5bYzNDk7Oxi/QFseYMk0KCvVO1GtqYJuCoU52g+IdxCMfEvvNGZvYpdABoT78cbyYyRmXNQ7jsV7Z2UsMWm5vHlDW7zdnPSOY968ArShjKm7YfM7eqepNVXQTeHvmXDhFFy/CJxtZL2MWhBCMOeGdrg42vPIT/spLrW+8cF1ZUNsOp9tPs7tPZoxsFVDveNYhsgboN1E2PgWpOzWO02tqIJubAd/gUO/aM2eA7voncZi+Hu6MOeGdhw6lcOcFTF6x7FI6bmFzPz1AOENPXh2pA1PIKqJkW+DV1NYcjcU5eqdpsZUQTemc8nw9+MQ2B16P6Z3GoszrE0j7ujRjM+3HGftsbN6x7Eo5eWSx385QG5hKfNu6Wh7a7XUlosnjF8I509qn7AtdCijKujGUlIIv9wBCLj+U7B30DuRRXpmZGsiGnvy+K8HOJOjhjIa6vMtx9kcn8mLYyJsZ1lcY2vWA/o9DQd/gr1f652mRlRBN5aVz8CZ/dpd8/rN9U5jsVwc7fnolo4Ul5bz8I/71drpBohOzmbOihiGt2nELV2D9I5j2fo+AaED4Z8ntZUZLYwq6MZw4GeI/gJ6PQytRuqdxuKF+NVj9vhIdiVn87ZqW3dN6bmF3Pf9XgLquzJnQjs1G7S27Ozg+s/AzVf7xG1hbetUQa+ts0e1RfOb9YaBL+qdxmqM7xjArd2C+HRjEssOXrYWnII2tf+BH/ZxobCEBbd1xsvVxhfeMhb3BnDjl9r19KX3WdT1dFXQa6MwB365XRuaOOELdd3cyF4a04aoZvV54teDxKRd0DuO2ZmzPIZdx7N58/p2tGrkqXcc6xLUXVsaIGYZbJmrdxqDqYJeU2Wl8OsUbWTLhC/BQ435NTYnBzs+mdwJT1cHpn2zh/P5xXpHMhtL951i0Zbj3NGjGdd1vFa/GaXGejwAkRNg7asQ87feaQyiCnpNrXoOEtfC6LnQvJfeaayWv4cL8yd3Ji2nkAd/3KdukqLdBH1y8UG6Bfvw3KhL2/sqRiMEjPsImnSC3+6BtMN6J6qSKug1Ef0l7FwA3e+HTrfrncbqdQqqz6zrItkcn8kLfxyx6fXTT2blM+3bPTSt78qnt3XGyUH9CpuUoytM/EEbp/7jJMjL1DvRNamfhupKXA//zISWQ2Hoa3qnsRk3dQnkvv6h/LjrJPM3Wmc/yKrkFJRw19e7KSuXfH5HlFqnpa54NoaJ30NeOvw4EYrNt8G5KujVcWov/DwZGoTDDZ9rq7QpdWbm0HDGdWjCWyti+WO/bTXFKCwpY/q3e0jOzGP+5E6E+NXTO5JtadpZG854ao9276ysRO9EV6QKuqGyEuH7G8HNByb/pn0EU+qUnZ3grQnt6Bbsw8xfD7A1wbw//hpLaVk5D/24j+1JWbw1oR09QxvoHck2RYyFUe9C/Eqz7XSkCrohLpyBb68DJEz+XfsIpujC2cGehbdFEdKgHnd/HU10crbekUxKSskzSw6x6uhZXhoTwfWdAvSOZNui7tKWB9j/Pax52eyKuiroVclNg69HQ14W3PorNGihdyKb5+XmyLd3d6WxlwtTvtzNwdTzekcyCSklry47yq97UnloUEum9ArWO5IC2kqqUXfB1vdhg3m1llQF/Vpy0+CrUdoZ+uTftOtoilnw93Dh+3u64eXmyO1f7OLIacuaol0VKSWv/HWUL7cmM6VXcx4d3FLvSMq/hICR70LHybBxDqw3n6KuCvrV5KbBV6P/v5g366F3IuUSjb1c+fGe7rg52jNp4Q72nDindySjKC+XPL/0MF9tS+bu3sG8ODpCrdFibuzsYMw86DAZNr6pFXUzuPyiCvqVZCbA50PgwmmYvFgVczMW6OPGL9N74OPuxG2f77T4G6UlZeU89dtBvt95knv7hfDcqNaqmJsrOzsYW6moL38Sysv0jaTr3s1RajR8MRSK8+HOv6BZT70TKVUIqK8V9SAfN6Z8uZu/D57RO1KNXCwqZerX0fy6J5WHB7Xk6eGtVDE3d/8W9R4PwK6FsPguKC3SL45uezZHscvh6zHgVA+mrlLXzC2Iv4cLP03rTtsAL+7/YS/z1sZb1IzS9AuF3PzpdrYmZPLm9W15dEiYKuaWws4Ohs2GobPg6FL47gbI12f0lSroAOXlsGGONgusQRhMXQ2+oXqnUqrJ282J7+/uxviOTXl3dRyP/LyfwhJ9PwIbYs+Jc4z9aCvHM/NYdEcUE1WTCsvU80Ft8lHKTljYD9IO1XkEVdALzmmzPze8rnX+vmuFWjnRgrk42jP3pvY8MSycP/afZvwn20jMuKh3rCuSUvL1tmQmLtyOk4Mdi6f3ZEC4v96xlNpodxNMWaGtxrpoiNY0vg7ZdkE/vgnm99Zmfg1/U2sf5+iqdyqlloQQ3D+gBV/cGUVaTgFj5m3htz2pZnUJJutiETO+28tLfx6hb0s//nqgNxFN1OxjqxDQGaZtgCYdYMk92kqNddT5SOj1Qx4VFSWjo6N12TfFebD+ddj+sXZp5fqF6nq5lTqTU8DDP+1n1/FshkQ05LVxkTTyctE108ojaTy75BC5haU8NjSMaX1CsLNT18utTlkpbH5XG6vu2US7eRo6oNabFULskVJGXfF7NlXQpYRjf8GKZ+BCKnSeot3McHKv2xxKnSorlyzanMR7a+JwtLPjieHh3NI1CAf7uv2AmpyZx6y/j7Hm2FnaNPFk7k0dCG/kUacZFB2kRsOSaZCdCG2u12qOZ5Mab04VdICUXbBuFhzfCA0jYeQ7any5jTmRlcdzvx9mS0ImIQ3ceWJYOMMjG5l8NEnWxSIWbkriy63JONoLHhjYkqm9g9Va5rakpBC2faidsds5aEtvR91Vo03ZbkEvL4ekdbBjPiSsAXc/6PsERE1V/T9tlJSSNcfSeWtFDPHpF4lo7MnU3sGMad/E6AU2JTufr7cl8/3OkxSWljG+Y1OeHt4Kf099L/koOso+Dqueh4jroN2NNdpErQu6EGI48AFgDyySUr55yfdFxfdHAvnAnVLKvdfapkkLemY8HP0D9v+gfcxx99MG/ne9R11eUQBtSdrf951i4aYk4tMv0qCeM2PbN2Fchya0C/Cq8Vl7TkEJG2LTWbwnlS0JmdgJwbj2TbhvQAta+Ks1zJXaq1VBF0LYA3HAECAV2A1MklIerfSakcCDaAW9G/CBlLLbtbZrtIJelKv91Us/Cid3wImtkBmnfS+oh3Y2HjEWHJxrvy/F6kgp2RyfyXc7TrAhNoPisnL8PZzpGepL12BfwhvVo4WfB15ujpe9t6SsnLScQmLTcjmYep7dyefYnZxNabmkiZcLN0YFclOXQJp6q5FTivFcq6Abct2hK5AgpUyq2NhPwDjgaKXXjAO+kdpfhx1CCG8hRGMppfHnYMevhpXPaiNVii5CUaXhQM6eENhVK+Ktx4CX6oauXJsQgr5hfvQN8yOnoISVR9LYHJ/JloRMlu4//b/XOTvY4e3miKujPSVlkqLScrLyiv63HpOdgLCGHtzdJ4TBrf3pGFQfezVyRaljhhT0pkBKpcepaGfhVb2mKfCfgi6EmAZMAwgKquFsOBdvaNhGu3Ti6A4ejcAnBBq0BL9Wqi2cUmNero7cFBXITVGBSClJyS4gPj2XhPSLZOUVcz6/mIKSchztBc4Odvh7uNDE24XgBvWIbOqJm5O6L6Poy5CfwCudZlx6ncaQ1yClXAgsBO2SiwH7vlxgFwj8qkZvVRRDCSEI8nUjyNeNQa3VzGHFMhhyWz8VCKz0OAA4XYPXKIqiKCZkSEHfDbQUQgQLIZyAicCfl7zmT+B2oekO5Jjk+rmiKIpyVVVecpFSlgohHgBWog1b/EJKeUQIMb3i+wuAf9BGuCSgDVucYrrIiqIoypUYdBdHSvkPWtGu/NyCSl9L4H7jRlMURVGqQ809VhRFsRKqoCuKolgJVdAVRVGshCroiqIoVkK31RaFEBnAiRq+vQGQacQ4elLHYp6s5Vis5ThAHcu/mkkp/a70Dd0Kem0IIaKvtjiNpVHHYp6s5Vis5ThAHYsh1CUXRVEUK6EKuqIoipWw1IK+UO8ARqSOxTxZy7FYy3GAOpYqWeQ1dEVRFOVylnqGriiKolxCFXRFURQrYXEFXQgxXAgRK4RIEEI8rXee6hBCfCGESBdCHK70nI8QYrUQIr7iv/X1zGgIIUSgEGK9EOKYEOKIEOLhiuct8VhchBC7hBAHKo7llYrnLe5YQOsBLITYJ4RYVvHYUo8jWQhxSAixXwgRXfGcpR6LtxBisRAipuJ3poepjsWiCnpFw+qPgRFABDBJCBGhb6pq+QoYfslzTwNrpZQtgbUVj81dKfC4lLI10B24v+L/gyUeSxEwUErZHugADK9Y098SjwXgYeBYpceWehwAA6SUHSqN17bUY/kAWCGlbAW0R/v/Y5pjkVJazD+gB7Cy0uNngGf0zlXNY2gOHK70OBZoXPF1YyBW74w1OKY/gCGWfiyAG7AXrWeuxR0LWqewtcBAYFnFcxZ3HBVZk4EGlzxncccCeALHqRiAYupjsagzdK7ejNqSNZQV3Z0q/uuvc55qEUI0BzoCO7HQY6m4TLEfSAdWSykt9VjeB54Eyis9Z4nHAVpP4lVCiD0VzeXBMo8lBMgAvqy4FLZICOGOiY7F0gq6Qc2olbohhKgH/AY8IqW8oHeempJSlkkpO6Cd4XYVQkTqHKnahBCjgXQp5R69sxhJLyllJ7TLq/cLIfrqHaiGHIBOwHwpZUcgDxNeKrK0gm6NzajPCiEaA1T8N13nPAYRQjiiFfPvpZRLKp62yGP5l5TyPLAB7T6HpR1LL2CsECIZ+AkYKIT4Dss7DgCklKcr/psO/A50xTKPJRVIrfjUB7AYrcCb5FgsraAb0rDa0vwJ3FHx9R1o16PNmhBCAJ8Dx6SUcyt9yxKPxU8I4V3xtSswGIjBwo5FSvmMlDJAStkc7fdinZRyMhZ2HABCCHchhMe/XwNDgcNY4LFIKdOAFCFEeMVTg4CjmOpY9L5pUIObDCOBOCAReE7vPNXM/iNwBihB+8s9FfBFu5EVX/FfH71zGnAcvdEudR0E9lf8G2mhx9IO2FdxLIeBFyuet7hjqXRM/fn/m6IWdxxo150PVPw78u/vuSUeS0XuDkB0xc/YUqC+qY5FTf1XFEWxEpZ2yUVRFEW5ClXQFUVRrIQq6IqiKFZCFXRFURQroQq6oiiKlVAFXbEaQgjfitX59gsh0oQQpyq+viiE+ETvfIpiamrYomKVhBAvAxellO/onUVR6oo6Q1esnhCif6X1wd0r1qXfXbFY0riK5+8UQiwVQvwlhDguhHhACPFYxWt2CCF8Kl63QQjxvhBimxDisBCia8XzPhXvP1jx+nb6HbFiq1RBV2zNc2jT4rsAA4C3K6aXA0QCt6CtGzIbyJfagkrbgdsrbcNdStkTuA/4ouK5V4B9Usp2wLPANyY/EkW5hIPeARSljg1FW8RqZsVjFyCo4uv1UspcIFcIkQP8VfH8IbQlAv71I4CUcpMQwrNiLZjewA0Vz6+ruJ7vJaXMMe3hKMr/UwVdsTUCuEFKGfufJ4Xohta96F/llR6X89/flUtvPEnU0s6KGVCXXBRbsxJ4sGLFSIQQHWuwjZsr3tsbyKk4C98E3FrxfH8gU1rwGvGKZVJn6IqteQ2ts8/BiqKeDIyu5jbOCSG2obUXu6viuZfRutIcBPL5/6VRFaXOqGGLilINQogNwEwpZbTeWRTlUuqSi6IoipVQZ+iKoihWQp2hK4qiWAlV0BVFUayEKuiKoihWQhV0RVEUK6EKuqIoipX4P3QwkONSiUZOAAAAAElFTkSuQmCC\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"Ω1 = 1\n",
"Ω2 = 1\n",
"Δ1 = 5\n",
"Δ2 = 5\n",
"\n",
"Δ = (Δ1+Δ2)/2\n",
"\n",
"H = Ω1*Ω2/(4*Δ)*(k2*k1.dag()+k1*k2.dag())+\\\n",
" (Δ1-Δ+Ω1**2/(4*Δ))*k1*k1.dag()+\\\n",
" (Δ2-Δ+Ω2**2/(4*Δ))*k2*k2.dag()\n",
"\n",
"rho0 = k1*k1.dag()\n",
"rslt = qu.mesolve(H,rho0,t,[],[k1*k1.dag(),k2*k2.dag()])\n",
"\n",
"\n",
"pl.plot(t,rslt.expect[0],label=r\"$\\rho_{11}$\")\n",
"pl.plot(t,rslt.expect[1],label=r\"$\\rho_{22}$\")\n",
"pl.xlabel(\"Tiempo\")\n",
"pl.legend(loc=0)"
]
},
{
"cell_type": "markdown",
"id": "30d28118",
"metadata": {},
"source": [
"## Comparación"
]
},
{
"cell_type": "code",
"execution_count": 21,
"id": "534e6f17",
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.legend.Legend at 0x7fdc080548b0>"
]
},
"execution_count": 21,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXQAAAEGCAYAAAB1iW6ZAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/Z1A+gAAAACXBIWXMAAAsTAAALEwEAmpwYAACEwUlEQVR4nOyddXgU1/eH39mNu7sHYjgEgiQEL05bKNSot1D3b11+VeotNaSUOqUtUpziEUhIkBCBQNzdfWV+fwwspECBNmGTMO/z8MCcmZ35zOXu2Tt3zj1HEEURGRkZGZnuj0LfAmRkZGRkOgbZocvIyMj0EGSHLiMjI9NDkB26jIyMTA9BdugyMjIyPQQDfV3YwcFB9PHx0dflZWRkZLolhw4dqhBF0fFC+/Tm0H18fEhMTNTX5WVkZGS6JYIg5F5snzzlIiMjI9NDkB26jIyMTA9BdugyMjIyPQTZocvIyMj0EGSHLiMjI9NDuKRDFwThW0EQygRBSLnIfkEQhMWCIGQIgnBMEITBHS9TRkZGRuZSXM4I/Ttg8j/snwL0Pv3nAeDr/y5LRkZGRuZKuWQcuiiKUYIg+PzDIbOAH0QpD2+cIAg2giC4iqJY3FEiz6UkKZXYb1Yzoq8nplYWGDg6YujlhZGXFwoTk864pMw1ilYrUlfRTE1JE031bbQ0qlBXl6I0NkFpYYO5tTEWtsbYOJthqi4GUQN2fvqWLXMN0xELi9yB/HO2C07bznPogiA8gDSKx8vL619dLO1AEiE7fqd2B9SeYxcVSkqcvbEZFkrAjEmYhw1DMDT8V9eQuTYRRZHS7Dry0qooTK+mNKcOjUp7gSMbgIp2FktlGc5GGXhNnorPyL6YWhpdFc0yMuciXE6Bi9Mj9E2iKPa9wL7NwLuiKMac3t4F/E8UxUP/dM7Q0FDx364UfWlNEhsTsol5dDiGVRW05eXy4y+7cMxJJ6gmHyONCoW1NVaTJmJz882Y9unzr64jc21QV9FMakwRpxJKqa9sAQEcTQpwc1VjN2oydq7mmFkbYbLzKQzSfkWLEvW4t2nqfTv1VS1URv1J2clCituCadLagQDuAbb0GeWKX/GHKCuSYe73YOmi71uV6QEIgnBIFMXQC+3riBF6AeB5zrYHUNQB570oMwd58HNCAVGVWqb370uFux8fRkP/ETdyIreCXwcLOB+JpXbTZmp+/4Myz94obruL8DtvQBCEzpQm040oyarlyF95ZCeVgyDgGWxL2AxfvHPewuTET9AE9BoOzm6g1ULudhgwF0X5cQzSf8UkcgF2buZ4R6+AATaIxsepyK0iq89XpMeX8Ne3aZgIoxhgXkW/vV9iPONNfd+yTA+nI8IWNwB3nI52GQ7Udtb8+RlCfeywMzdi9/EyAP5KKwVg0Y39wciYrVa9cX//fXrv20vMdfNRV1XhsOglsubcRENsbGdKk+kGVBU1suXrY6x5/xCFJysZNNGTO94eyYxHBxI4zAmT3M3Q+zpQGMLRn6UPFR+FpkroNR6CpkPRYWgog9oCKE2BgEkIQZNxbD1AWFgz898cwYxhB3Ayzia+4XZ+2BrGoW05qFUa6XxaLbQ26K0NZPSDpqGRso8/oTEurlPOf8kRuiAIq4AxgIMgCAXAa4AhgCiKS4AtwFQgA2lMc3enKD0HpUIgMsCRvSfL0WpF/kotIdDZkhA3K4Z42xKTUQlAi7EZH1oOwu/WEbgn7uOxwijy770Py0mTcH7heQxdXTtbqkwXoq1FTdz6LFL2FWBgrCRsSDX98x/CSDURbH+SDipOguZq6HcTCAKk/QmT3oLMXYAA/uMkJ77nbTi1AzRt0ud6XwfGFtK/M/cgOAbjVf0jXiOGUmbjRsKWHOLWm5EWU8SoOb3wTXsc4dQ2uPkXCJyil/aQuXqIokj9tm2UvrsIdVkZgpER5sOHd/h1LifK5ZZL7BeBhztM0WUyJtCRdUcK2ZNeRmJuNQ9G+gMQ3tuBD7anU17fSnx2Ja1qLa9c34/HmjWsmDqVVxqPULFkCQ3R0Yj3LsTxjvk4WJlebfkyV5nsYxVErUqnoaaVfqPdGTrDF9OVI0DRAsc3Qnk6OAZC5m7pA35jJGd9cps0Oj+1E1wHgLkDmNmDhQuc2g6qFrD2kj4rCGDfG7L2gPcIqC+GwKk4+Y1m2oEF5Pf9hOiUELYuScHfuA+R1jGYRn8sO/QejqqkhOKXXqYxNhbjkGA8Pl+M6YABnXKtbrtSNDLAEYUAj646gkYrMmOAGwDhvRwA2J9ZwZbkYhwsjAnztSeitwMxObXYPfAA/ps30RwyAL74hL3X30pDYafOEMnoEVWbht0/HmfLV8cwMjVg9rNDGH1LIKaaMqjMgPAnQWEAR3+RPpC1F1z6gYWj5GgFJRxcDvnxEDhVOkYQIGCSNHo/tR36zZZsII3gc2Lh8A/SlE3AdWDlBg6BeLZsZt4rwxjeN5vs1qGsqllGdgbS1I1Mj6Ru61ayZs6i6cgRnF9+Gd/ff+80Zw7d2KHbmBkxLsiJpjYNg7xsCHSxBKCvuzXWpoZsTS5h1/EypvVzQakQGNXLgcrGNtJL6zF0d+ed8PtYPHAO/qWZ5My6noaoKD3fkUxHU1nYwO/vJnJ8fzGDJ7gw98WhuPhZSzuz9kl/950jjcZT10FLHeTFSU4ZwMwOfCNOz6OL0Of6sycPvRcQJKc9+M6z9n5zQN0Mh76D3pOkc4B0jdz9KLVtDFF/yk2D1mBma8aWmheI/TEereZC4ZEy3RVtaytFL75E4ZNPYeTrg9+6tdjdfhuCUtmp1+22Dh1g0ez+vDAliC9uPZttQKkQmBTizLbUElrVWmYNcgdgVC97AGIzKiisaSYht4bA++bzyrTnqTK3IX/BQiqWLOVywjhluj6Zh8v4Y1EiLY0qZo7LZETKCJTR7509IHsfmDmAUwiEXA81uRD7KWhVkiM+Q+TzYGINA2+XplXO4DYQFuyDhdFg53vW7jkMAqeBkQWMfuas3W8MqJpg/+dQk4dD2BhuemkU/Sx3cTTJnD8/PUpjbSvsex/ecoGk1Z3TMDKdjqq4mNzbbqd27VrsFy7A56efMPL2virXvqw49M7gv8ShX4r8qiYW/nSIoT52vD7zbAz6uI/24mVnxgg/e97deoK9z4zhiz0ZRCXns7ZhD/VbtmB53XW4vbdIXnXaTRFFkUNbc4nfkIWzrxVT7uuN+TcDobkKDEzh6ROSg/4oCLxHwk0roakKPgqU5sxN7eCZk6D8j4vSRPHsNAxIo/8PeoGmFYyt4enjYGQOv95GeroReyvuxMRcyXSjx7FXZIKNFzx+rP05ZLo8TYePUPDoo4gtLbi9/x6W48d3+DX+KQ69W4/QL4annRmbH4to58xBml+Py6rkhwO5DPKywcfBnPBeDpS1CVQ/+QpOzz5L/V9/ceim29my/4Se1Mv8WzRqLTu/SyN+QxYBw5y5/qlBmFcdlJz5uFekqZDk36UXoA0l0qgZpGmR8KdAUMCYF/67M4fzHbGJFYw4HTsw6jHJmQP4hBMoruXGBz3Qtjaztvz/KPB9AWryoCT5v+uQuWrU795N3t13o7Awx+e31ViOH48oiny5J4Ml+zKvytN/j3ToF2NyXxdaVFoKa5qZP1x6BBp1+iVqTGYl9vfeQ+VTr2KceQLF4ws5eihdn3JlrgBVm4YtXydzMr6UsJm+TLg7BANDJeTGSi89hz8Ezn0haRWc+kv6kP/YsycY+wK8UABhD3SeyAmvwbOZ7adifMIBcFQlMsf/C8yNGtl4MIyTzRHSC1eZbkH1779T8MijGAcE4LNqFcb+UtTdtpQSPtiezqKtJ9h1vPNffl9TDn2Enz2vTg/hsfG9uX6gNLfuaGlMkIslsRlSbo5vDPx4Z/RCHJtr0Dx0L225F63HKtNFaGtWs+nzJPLSKhlzWyChU33PrgjO3Q9ug8DIDAbeCoWHYN974NJfmtY4lzOj5s7E3KH9tlMfMHeCna9jWbGP2TfW4eJnzY7aJ0iLr+58PTL/mYrlyyl55VXMR43C+7uVGNjZ6fZtSSnBztwIFysTvonJAiCzvAGNtnNG69eUQxcEgXvCfXlqYgAKxdlH4lG9HDiYU0VmeQP7TpYz7MaJ/Hb7S2hbWsi9407ZqXdhWptU/PnpEUoya5k0sYE+uY9DaZq0U9UMhYeluXKAfnOlefS2Bhjwj8srrh4KBYTeI8WsG5pjPOwmZjw6AE/HSvZkTiB1X550XEUGbHgUCjrnvZPMv6Ni+XLKP/oYq2nT8PzqSxTmZwcFoigSn1XJqF4OzB/hTVxWFRll9Vz/RSxvbEztFD3XlEO/GOG9HGhTa3ni16MAzB3qSVBEKM+OXICmpYXcO++SnXoXpK1FzcbPk6goaGDyvb3ofeJByNgBm56UDihIlKJWvE47dAtHuHMjzPwcwhboT/jfGfUYjH0ZbvsdTG0xMFIydY4x3saJ7F2VQfLeAli3QIptX3Ov9MJVRu9UfvONzpm7vbfovOyuuZVNlNW3EuZrxw2D3BEEeOSXI9S3qhkX7NwpmmSHDgzztcPIQEFyYS0zB7jhYWvGMF87cqzdyH7hfcTTTl1VWKhvqTKnUbVp2PzlMcpy67nu/r74miRKI++g6ZAfJ43Ss6OkF51eYWc/6DkUBt8Bis6NB74ijMwh8lnwGaUzGfiPYorNe/i41xL160nSM8zBdSBU50hPHTJ6pfK77yj78COspk6VnLnB+Yvu47OlFCTD/exwszFlfJATJ0rqcbAwZqS/fafokh06YG5swIc3DWB6f1eemxwEQJCLFRbGBsSItnh9txJtYyOpt9/Nys1H0HbS/JfM5aFRa9m2JJmijBom3B2M30BHyD0AhmYw7SNQGsGhlZCxEzyGgqmtviVfORaOKJ16cZ3XL7jbl7Or9lGy+30lveA9sUnf6q5pav/8k7JF70khzu+/h2BgQGOrmufXHOOb6CzdcfFZVdibG+HvKOX4eev6ftwxwpuvbx+MobJzXK/s0E8zc4AbX9w6GDcbKa+LUiEw2NuWxJxqTIKCaHrtPSgtxuat51h/IEPPaq9dRFFkz48nyEurYuztQQQMPZ1jPG8/eIRKOcf7zoGDy6SMiL0n6lfwf8F3NAbZfzFV+RiO1nVsX1VOkeUsyInRt7JrloboaIpeehmz4cNx++B93ch8RUw2vybk89bm4xzOk15mx2dXMczXTveC3sXahDdm9WWoj91Fz/9fkR36PzDMx5b00npqm1R8U2fNp6PuIqCmAOO3XkJUqfQt75rk4MZs0uNLCJvpS8goKX8PLbVQmnp2rjziKWm0bmoHg+7Qn9j/ypC7QWGIkaHI9If7YWVvwuasm6nKLZVT7+qB5mPHKHj8CYx798bji89RGJ2tSrXreCkhrlbYmBmydF8mORWNFNY0M9yvc6ZWLobs0P+BM7+kqxPz2JZSQu/rp3B4zkL8spMpeOMtOU3AVSYttojELTkED7dnSGjb2R35CSBqpQyHAA69pVWWjx4Cy855+XRVcA6BB/fDQwcw9Qpg+qMDUBoq2Vz1PM3p8WePa62X8qvLdBptBYXkL3wQAzs7vJYtRWlhodvX0KompaiO8cFO3DrMix1ppXy3PweAcUFOV1Wn7ND/gVAfO1ytTXhnywkE4O5RvnjNv4Xfeo+l4fffqP75F31LvGbIT6ti78/peAXbElm/AOGrYZCyVtqZt1/Kiuh+zmpoC8ezibG6M44BYC8tUrGyN2Xqgj40auzYuvp0vdPU9bDIC36br1+dPRhtYyMFDz2EqFbjuWwZBo6O7fYfyq1GoxUZ5mvH/BHeCILAd/tzCHKxxNPO7KpqlR36P6BUCDwyrhcKAe4f7YebjSmDvWz4oc8USvsNo/Tdd+UKSFeB2vImtn+Tgp2rGdeNLUJZeVzasfN10KilxUOuA84WmOjBuAS6MN5vC8WV1uz58TjiX69ITycnNkGZnK6ioxG1Wgr/9xytGRm4f/wxxn6+5x1zMLsSA4XAEG9bXK1NeWJ8bxwtjXlhavBV1ys79EtwW5g3aW9M1kW/WJoYEuxuw3dj7sLY35/CJ56kNStbzyp7LqpWDVuXSDlNpizsj1FhlDQ/ftN3UobEw99Jucp7dXwSpK5K78H2DLP8lfSDpRwtGghTP5TCM1PW6Ftaj6N88WIadu3C+fnnsAgfdcFjDmZX0dfdGjMj6QXpo+N7k/DSBCIDHC94fGciO/TLwMSwfcxyqLcd8SWtuHz5JYKhIRkPPszv0emo5ZzWHYooiuz+4ThVRY1cd19frB1NIe+AFM0SPBNsfWDz09IINWi6vuVePXzCCTVbjb/xfg7U30Gh1Q1SeGbWHn0r61HUbdtO5ZKl2Nw0B9v50pRWTVMbb25KY88JKS9Li0pDUn4tYb5dY3pPduj/gmG+djSrNJzEAvM33kHMzaHs/97g25isS39Y5rI58lceGYfKGH69P54hdtLLv5Jk8BohLQwa/yogQK8J0pTLtYJPBIJDL8ZZf4G1jcj2FcdpcB4vLThqrde3uh5BW04OxS+9hMmA/ri88oou9PDDv9JZEZPNgh8PUVLbwpG8Gto0WobJDr37EuojLVRJyK7iB5UTP4dcx7iCw5T8tErPynoOBenVxK3PpFeoE4MmnU6iVXA6msXrdHHdvrPhf1lwy+prK2+4gRHctxOjh3cx5bFRqNq0bD8yFI1WkBZYyfwntC0tFDz+BIKBAR6ffIJwOjxRFEV2ppUxwNMGjSiyLCqLmIxyFIIUQNEVkB36v8DJ0gRfB3PWHSnk57g86mffRlW/UK6PXU1OTIK+5XV7mura2PFtKtZOZoy9Pehs5sS8OGmu2GPo2YPN7EB5yVrnPQ9TW3Dug52bOePmB1FSCPsb7pYqMcn8J0reeovW9HTcPngfQzc3nb2gupmSuhbmDPFg5gA3Vh3M44f9uYzq5YC1aQfk0O8AZIf+L7l+oDtpxXU0tql5clIQTu++S5WxFbXPPYOmpkbf8rotolZk1/fHaW1Uc93oPIwOfiJFsoA0f+7SD4wt9Suyi9E71Jl+Yz041jiVnKQSyajVwrHf5bwvV0jNuvXU/rEG+4ULsBg9ut2++OwqAIb52PHgGH+aVRrqW9XMGeKhD6kX5Boc2nQMCyL90Gi1hLhZ4e9ogehgzouR9/D6X59S/NrruH/6ydmRpcxlc3RnPnmplUROMcEh+n7JaGgmpZjNT4Ahd+lVX1dl5I3+FB05ya7s6dxcUop5/kbphbHSCJ5IllIiyPwjrVnZlLzxBmZhYTg++uh5+w9mV2JjZkhvJwsUCoHfFoygqKaZmQPcLnA2/SCP0P8lJoZKnpoUyOS+roCUa91+yCDWD5pO/fbt1K5br1+B3ZCS7Fri1mfiP8iRPsabJWfkGSYVTj65TSohdw2FJ14JBoZKrptrgxpjdn5zFHHvh2DuKNVJPfabvuV1ecS2NoqefRaFsTFu77+PoDw/G+fB7CqG+tjpaikM87Xj+kHuXWrgJjv0DmSorx0rPEZhMDiU0rfeknOoXwFtLWp2rEjF3MaYsfODEPLjwH0ITHkfWmrg97uk4hTeF44FlgHb/kOJsPmBggJDjpQOh5u+B8cgyNylb2ldnvIvv6IlNRWXN9/A0Pn85fqldS3kVDYxrIu8/LwYskPvQMJ87dAKCk7e8xQYGJD39DPEnijptHJTPYnYNRnUVbYw4Z4QjA1UUHxUCk90Gwj9bpIOClsglZKTuTAGxgQPNMLfJJb4xtspE/qB/3gp8kXVom91XZamxEQqly3Des5srCZK2Tlrm1R8G5NNXmUTII3OgS4TnngxZIfegQS7SjnU99cb4PT666hSUtj2v7f5bNcpfUvr0uQkV5AWXcSgiV649bKR6n5q1ZJDB7h+CTwUBxNe16fMboFw/VeMvWsgZtYm7PzuOGr3kaBpheIkfUvrkmjq6yn633MYenri8sILOvubm9N4Y1Mat62Io0Wl4UBWJeZGSvq4WelR7aWRHXoHojydz+FgdhVb7EP4yyuUeSd3cWBztDxKvwgtDSr2/HgCOzdzwmb4Sca8OECQqguBFJboFHxtxZr/W4wtMB44g7Hzg6kubuRg6ukY/vw4/erqopS8+Saq0lLcP3hfVw+0qrGNP48WEuRiSX5VM0v3ZbH7eBkRvR0x6KTCFB1F11bXDRnma8epsgZeWpfCwal3orWz597YHzmcUapvaV0OURTZ+0s6LY0qJtwdgtLwdHfMOwBOId2z0lAXwauPPSERbhzZV0GxyTjIP6hvSV2O+l27qNuwEYeFCzEdcHalcUJOFSqNyFvX92VaP1c+2XmSkroWJoR0/VTMskPvYGYP9sDKxABDpcAzNw7B+f/+D+/6Usq++FLf0rocpxJLyTxcxrAZvjh6no4t12ok53NmNajMv2bU7F5Y2pmwq+wuVLlH5OLS56CpraXk9f/DOCgIh4XtC4Yfzq3GSKmgr7s1L04Lxs7ciCAXS6b3d9WT2stHjkPvYFysTdj9zBia2zRSLmSvcWwNHMnQHWtoTp6Lab9++pbYJWiqayPq15M4e5sxKOsOWKaE+eukDIpt9Wfnz2X+NUYmBoy7I5g/P2khrnQyEVVZUkKzTU8CIkz7GJRdY4Xj1aZ00Xuoq6rwWPI1gmH7NjiUW01fdytMDJW425iy//lxGCkVunDFrsxljdAFQZgsCEK6IAgZgiA8f4H91oIgbBQEIUkQhFRBEO7ueKndBwcL43aJ7XPm3U+1iRVFz7+Atq3tHz557RD920lUrRrGhcShKEuS6n/uXQSn/pIO8AnXr8AegkegLf3CzDjWNJ2i+MOQ8A0c/h4O/wDJf+hbnl5oiI6mdt067O+7D9M+fdrta1VrOFZYyxDvs9N9JobKbuHM4TIcuiAISuBLYAoQAtwiCELI3w57GEgTRXEAMAb4SBAEI2QAGBDswacD59CWmUmFPPVC9rEKMhLLCJ3ig13lNqnS0OA74eBSiPoIPIaBVdd/vO0ujLh5CJYGFezZrkWz8x3wHwfWXnB8g76lXXU0DQ0Uv/IqRv7+ODz80Hn7U4vqaFNr2zn07sTljNCHARmiKGaJotgG/ArM+tsxImApSEumLIAqQN2hSrsxod52JDoHUx4+kcpvv6Ul/aS+JemNtmY1UavSsXMzZ/A4ZylE0XsEjH0JjK2k1aCh9+hbZo/C0NSQMeHV1LQ6klgzHaZ/AgGTIGsvqFv1Le+qUvbBh6jLynB75+12RZ7PcDi3GoDBXj3XobsD+edsF5y2ncsXQDBQBCQDj4uieF61B0EQHhAEIVEQhMTy8vJ/Kbn74WFrirOVMRtHzEFpaUnxq69yvLDmmgxlPLA+k4aaVsbOD0JZfkxamu41QirmvCAK7toCA27Wt8weh9fsOwnwq+Nw42yqmh2lUbqqCYqO6FvaVaMpIYGa1auxu/NOXVRLZUMr644UUN+iAqQFRJ52pjhZmehT6r/mchz6hSaP/u6JrgOOAm7AQOALQRDOi8AXRXGZKIqhoiiGOjpe/fJM+kIQBEK97YgtU+H8/HO0JCXx+VMf8NqGFH1Lu6oUZdSQsq+Q/mM9cPG1lsITQcrXAmDrDT6j5HjzzsDQhFELp2JoasCen04gug+T7HnXRv50sa2N4tf/D0N3dxwfO5t468GfDvPk6iQeXXUElUbLgcxKwns56FHpf+NyHHoB4HnOtgfSSPxc7gbWihIZQDYQ1DESewahPrYU1bbwp2N/jjj25u60LeyITqPu9Migp6NRadn70wks7UwIm3l6AVHuAXAIAPPu+wXqTphZGRE+pzclWbWkHmmT2j7v2lhwVLnyO9oyM3F59RUUpqYAVDS0cjCnCntzI/aml/PyuhTqW9WE9+q+g83LcegJQG9BEHxPv+i8Gfj725Q8YDyAIAjOQCAg12M7hwnBzggCvLohjY3jbsccLXcfXcfOtGtjwdGRnXlUlzQReWsgRiYGUr7u/Dg53vwqEzjcBY8gWw6szaDBaYzk0LU9uxZuW34+FV99heWkSVhERursZ/KzLLtjCIO9bFidmI+1qSFjg3qwQxdFUQ08AmwHjgO/iaKYKgjCQkEQFp4+7E1gpCAIycAu4DlRFCs6S3R3xNPOjJemBjPIy4aXF1yH44MLiSxMonBbz8+EV1fRzKEtOfgNcsS7r71kLD8BLbVyvPlVRhAEIm8NRKMRic0dI2WyrEjXs6rOQxRFSt58E0GpxPnFF9rtO5hdhamhkv4eNnw6bxBzQz346rbBmBl13+U5l6VcFMUtwJa/2Zac8+8iYFLHSut53Bfhx30R0nSD9r57OfHzHwxeuxztk3N1j4E9kZjfT4EA4X2PQ14TeIVBbqy0Ux6hX3VsnMwYMtmbgxuzCbHth2denJQrp7YA2prAMUDfEjuM+r920BgVjfMLz2Po0r7IR3x2FUO8bTFUKvCyN+P9Od2/0Li89F9PKIyMyJ7/CHYNlRR9vVTfcjqNnOQKspMqGDqwGssd98O3k6TKQ6d2gK2v9EfmqjNokhdWDiZENzyIJucgVGTAl2Hw5TDIjtK3vA5B09BI6dtvYxwcjO1tt7XbV9uk4kRJXZdPh3ulyA5dj/iMC2ePxyDqVq6kLT//0h/oZqjbNESvPomtixkDmj8Ba0+w8YIfr4dT26H3RDmiRU8YGCqJmBdAtcqVpKNK2PCIlEdHaQTxPWOAUfH5YtTl5bi+/hqCQfvJiIScKkRRqmHQk5Aduh4Z6GXDt32noxEUlL67SN9yOpxD23Opq2hh9Ex7lOVJELYQZn8rxT+DXB9Uz/j0c8DHV01C5RQastNh6gcw6HbI3A2qZn3L+0+0njpF1U8/YzN3brtMimeIz67EyEDBAE+bqy+uE5Eduh6xMjHE0deTqOEzadi9m4aonvGoC1BT1sSR7Xn0HuqMh1GyZPQeIeU4v3833L8HnPv880lkOp3wu8MRFSbst1wsOfPAKdIPbjeOTxdFkdJ330Vhbo7jE49f8JiD2VUM9LTBxPD82qHdGdmh65mhPrYscxyKoY8PJW+/Q1RqIbXN3T82Pfb3UygMBEbN7iWFxhmagUt/aafbIHAfrF+BMgBYO5kxeIofp7LMKThZI72kFhSQu1/f0v41Dbt20bj/AI6PPYaBrbSEP6WwliX7MqlrUVHbrCKlqK7HTbeA7ND1TqiPHbUagcYFT6DKzWXjKx9z2zdxiN04d3VeaiU5yZWETvXB3MZYGu15hF6zqVq7OoNPvyCNWpWOxsAcXAd0W4eubW2ldNF7GPfuje3N8wBoaFVz+4p4Fm09wZO/HmXPiTI0WpGxQecXg+7uyA5dz4zu7YCJoYKH0g3Z79qXW9N3UJKRp1v00N3QaLTE/H4KK0dTBoz1hLZGKE2RMijKdEkMjJSE39Sb6pImUqMKwXsUFCR2y8LSVStXoioowPmlF3UvQnemlVLTpGJ6f1d2nSjj6d+TcLQ0ZqCHjX7FdgKyQ9czNmZGzAv1pLCmmX0TbsXEQMF9qZvZe7J7Ji9LjSqiuqSJUbN7SSXlio6AqAVP2aF3ZXz6O+ARZMvBTdm0OJ8uLF10WN+yrghVSQkVS5dhOXEi5sPPrm9IyKnC0tiAT+cNZEKwExqtyJMTArpNjvMrofsuiepBvDw9hOv6uDDA04ZGizwiv17C4qh4mNy90uG0NKo4uCkLjyBbfAeczs9SkCD97TFUf8JkLokgCITf1JvVbx0kIdWTCJAWf3mP1Le0y6bsgw9Bo8Hpuf+1sx/KrWaQty0GSgXL7wilrlmNtVnPnP7rUg5dpVJRUFBAS0v3e9T7r7iamGCksMb0vvso/Hk1E/euoq75FqxMu0+dkIObsmlrUhM+wRhB3QKGptIiIjt/MOt5L6B6GvbuFoSEu5ESW0zfXpHYnjuPLopdes1A06FD1G3ejMNDD2Lk4aGz1zarSC+tZ2o/qWCKIAg91plDF3PoBQUFWFpa4uPjg9CFO09HI4oilZWVFBQU4Ovri+quBwhe/B7Hf1lL2L3dIzd4VVEjKfsKCQmsw/6PG8ClH9y+Tiqi0P8mfcuTuUyGzfDjVEIp+2tvYVrz01IqgN/mS3Pqt/4mpW3oYohaLaXvvIuBiwv2993Xbt/hvGpEEUK7aQWiK6VLzaG3tLRgb29/TTlzkEYN9vb2uieT4LtuIdPaDYMVX6Ft7foVZURRJPaPUxgaKwnjMzCygNJU+HoEqBqhzw36lihzmZhZGTFkqg85pc7kNfSClVMgY6eUxGvPW/qWd0HqNm+mJTUVpyefQGFm1m5fYk4VSoXAQC8b/Yi7ynQphw5cc878DOfet6WZMTsib8asqpyq73/Qo6rLIzelkry0KoZOdMS05jBEPA1jXoDGcnAdCD4R+pYocwUMGOuJlb0RsQ33oS06BsMWwOhnIScGGsr0La8d2pYWyj75BOOQYKxmzDhvf2JONX3crLp1BsUrocs5dBkJ2/CRHHTrS8XSpagrum4mYo1GS+wfGdg4m9HP85Rk9AmHyP/Bwwlw91ZQ9KzVeD0dpaGCUXMCqVJ5kBb8I0x+F4KmSdFKWfv0La8dVT/+iLqoGOf/PYegaO/O2tRakgpqCPW+dt7fyA69izLE25ZlIdPQtrRQvvhzfcu5KMdji6kpbWLkjf4oC/aDgak0KgcpDauR2T9+XqZr4jvQAbfeNsQnWtPaogXnfmBofjZqqQugrqqicukyLMaOxXz4+XP7R/NraFFpGeZ7bcyfg+zQuyyhPrYUWjhSPHY6NX/8wYsfreP7/Tn6ltWOthY1Bzdm4drLGp/+DlKYm+dQMOg+kTkyF+ZMGGNLo4rD2/NAaSClayg4qG9pOiq+/AptczNOzzwNgEYrsuDHRMLf282h3Gr2nSxDqRAY4X/tlDiUHfpF+Oqrr+jbty/e3t58/vnVHyG7WpviaWfKN75jaTQ0IeTP73htQyrHi+uuupaLcXRHHs31KkbO7oXQUgslydIqQ5kegaOXJQHDnEnanU9DdYuUvqEkuUtkYmzNyqZ69Wps5t6Esb8/AFEny9meWkpBdTPzV8Tz5Z5MQr1tsTbtuWGKf6fLvin4v42ppBV1rPMKcbPitRmXzvC3Zs0aduzYwZEjR6ioqKBfv348+OCDGBhc3ea6eagXH2xP55fe47k/ZSODK06x6Zg/wa5WV1XHhWisbeXIznz8Bzvh4msNJ7cDYrdaiCJzacJm+JFxqIz4jdmMHzYMtGpp9a+e/5/LPv4IhbExjo88orP9figfBwsj/lg4kluXx9HUpmFhpL8eVV595BH6BVi8eDHvvfcehoaGuLq6YmhoSFZWFvfeey9z5swBOG+7M1gw2o93b+zH3LefxNDNjYdPbmN3WkmnXe9KSNiUjValZfj1Ukk9cmNBYQjuofoVJtOhWDmY0m+MB+kHiqk0OJ0tMy9Or5oaDx6kYecu7O+/HwN7qUatKIrEZVURGeCEj4M56x8exZbHInpkAq5/osuO0C9nJN0ZqFQqjh07RkCAVFexuLgYBwcHAgICWLFihc6B+/n5tdvuDAyUCm4Z5gVA7eOPoXrueZwOxVDVOBI7c/3NU1cVN5IWW0zfSHdsnE6/9MyJleZY5ZegPY7QKT6c2F/Mge2VTHcI0KtDF7Vayt7/AAMXF+zuvENnz6popKqxjaE+0gtQJysTnKxM9CVTb8gj9L+RlpZGbW0tWVlZaLVaXnjhBR599FF9y8Jq+nQ0vv7cmbaNwxn6jQWOW5+JgZGCof4noKkKGiulRE6+kXrVJdM5mJgbMvg6b3KTKyk0nw75caDVSukA6q/uE2Pdlq20pKRIi4jOKax+KKcakNJRX8vIDv1vHDlyhNtuu41bbrmF/v374+XlxQMPPKBvWQhKJW7/exbXpkoqf12tNx1Fp2rITqpgsHsippvmw7fXQdIvUoxy4BS96ZLpXPqP9cDC1pj92SMQm2shJxo2PgYfBcKBL6+KBlGlovzzxRgHBp63iCgxtwpbM0P8Hc2vipauiuzQ/8bRo0eZPn068fHxpKSk8MYbbwBQWVnJwoULOXLkCO++++5521cD2zGjyfQIwn/br2gaGq7KNc9FFEX2r83A3MqAAY3vSznOq3Phr5fBvvfZ+HOZHoeBkZJhM/woKzMgQzMBfpgJh0+vYo7+CDSdX2WrZv16VLl5OD7++HmLiA7lVjPYy/aaXWl+hi47h64vjh49yoMPPnie3d7eniVLlrSz/X27sxEEgZwb78Z/8XOUfbMC14vUS+wsMg+XU5pdx9jwMgwzWmDWF1CTL43QRz0OCnl80JMJHO5C0q484uoX4GefhXLEArD1gVXzpMLSAdd12rW1ra1UfPU1JgP6YzF2TLt9tU0qMssbuXGwxwU/ey0hfwP/xt69ewkMDNS3jIviHzGUfe4DqF75Heryq1cEQ6PREvdnJnZu5gQZbQdrT3AMhN4TYM63UtkymR6NQiEw4oZe1NUZkDpwAwxfCP5jQWkkTcF0IjWrV6MuLsbpiSfOG4UfyZfmzwddIwm4/gnZoXczhvjY8kPwZGk+8auvrtp10w+UUFvWzPBZfigK4uUKRNcoXn3scA+0IWFzDm3NajAwlop+53feClJtUxMVS5dhFhaG+YgR5+1PzKlGIcCAHlhS7kqRHXo3w8nSBAMvb44NHEPNb7/zycqdrIzN7tSi0hqVloTN2Tj5WOHj3QJ1heDZ9fJiy3Q+giCN0lsaVCTtzpeMnsOkxUadVIO06sef0FRW4nh6irG5TcMDPyRy87IDFNY0sz21hKE+dpgbyzPIskPvhkwKceY9p3BaBQMsflzG/21MY18n1iBNiS6kobqV4bP8EPLjJaNcUu6axdnHCr+BjhzZkUdLg0r6cde0QfHRDr+Wpq6OyhUrsBgzBrNBgwD45WAef6WVcjC7inEf7uVUWQPT+7t2+LW7I7JD74bcE+6L2tqWtb1GE16UzOCWEtYeLuyUa6laNRzalot7gA0eQbbSo7WhmVSRSOaaZdhMX6lvbM89+7TWCdMuld9+i7auDsfHH9PZok6W09vJgvUPj8LSxIABHtbMGeLZ4dfujsjPKN0QNxtTdj8TSUvNIFrmHOShrJ08Ye1Bm1qLkUHH/kYn7y2gua6NsAf6Si+j8uPAfQgor52ERzLnY+9mQWCYC8l7CxgwbgQWtr5w5umtg1BXVlL1w49YTpmMSXCwZNNoOZRbzayBbvT3sCH+xQkohGu3MM7fuaxvvyAIkwVBSBcEIUMQhOcvcswYQRCOCoKQKghC18qC3wNxsjTBy9MJ+wfux/3kUXwK0zs8E2Nrs5rD23Px6mOPay8baG2AkhR5/lwGgGHTfRG1Iolbc6Q+kR8vrR7tICqXLUdsacHx0bOj87TiOhpa1YT5STlclApBdubncEmHLgiCEvgSmAKEALcIghDyt2NsgK+AmaIo9gHkqsBXCdtbb0VwcODO49s4nFvVoec+ujOP1iY1w6e5S4b8eBA14DW8Q68j0z2xcjClT7gbx2OKqLUZJZUcrMzskHOrSkqoXrUK6+uvx9jPV2c/mC318TDfa3uJ/8W4nBH6MCBDFMUsURTbgF+BWX875lZgrSiKeQCiKHatwoM9GIWpKU4PPUjfymzK90R12HmbG9pI2pmPv2sxjj/4w8Yn4OjPYGwllZiTkQGGTPVBoRQ4eDJIMqStgyM/wZJwSFnzr89b8dXXiKKIw0MPtbPHZVXhY2+G8zWYeOtyuByH7g7kn7NdcNp2LgGArSAIewVBOCQIwh1cAEEQHhAEIVEQhMTyq7go5t+g7wIXV4LtnDnU2TjSf9svHRa+eGR7Hqo2DcNU74B9Lzi0UvqC9p0NhqaXPoHMNYG5tTH9x3lwMqmBSsdZsPst+PNhqRDGhsegueaKz9mWl0fN2rXYzp2LkcdZV6PViiTkVBHma9+Bd9CzuJyXoheaoPq71zAAhgDjAVPggCAIcaIonmz3IVFcBiwDCA0N/WfPs/V5qVN0JC79YMqiSx7WVQpcXC6CkRGlN95B728/Im/dZrxvnP6fztdY20ry3gICvcqxayuCO1Oh6Chk7IDxr3aMaJkew6BJ3qREFRHXdj/TnDOluPS+s+G7aXBqB/S/shnY8i++QDAwwGHhgnb29NJ6aptVhPnJ0y0X43I8VAFwbkyQB1B0gWMqRFFsBBoFQYgCBgAn6YYsXryY5cuXtytwsXbtWnbs2EFZWRkPP/wwTU1NbN68Wbc9adIkvWr2u/lG8n7/EecvP8dr1hQEpfJfn+vQlhy0GpGh1uvBcBBYuUl/gqZ2nGCZHoOJuSGDJnoRvyGLkv9twcXPWkqva+4IJ7dekUNvPXWKuo2bsL/3HgwcHdvti8+qBGCYPH9+cURR/Mc/SE4/C/AFjIAkoM/fjgkGdp0+1gxIAfr+03mHDBki/p20tLTzbFebtrY20cbGRrddVFQk9u3bV7ddVVUl3nPPPRfd/i/8l/tXa7TiPXe/J6YFBonV69b96/PUljeJXz20W9z9Q4oovuEoitte/Nfnkrl2aG1WiSueiRLXfXxI1Gq1knHtAlF8z08Uz2xfBvmPPCKeGBIqqqqqztt33/cJ4sh3d3WU5G4LkChexK9ecg5dFEU18AiwHTgO/CaKYqogCAsFQVh4+pjjwDbgGHAQ+EYUxZQO/N25alyqwMVbb73Fww8/fNFtfaFUCKhGRZLr4EXF519QWdNIRUPrFZ8nYUsOgiAwdGAtaFrB6/zcGTIyf8fIxIAhU3woTK+h4ISULAvvUdBUAeXpl3WO5uQU6nfsxO7uuzCwtUUURXYdL+VQbhXNbRqiT5UzPvjaKil3pVxWHLooiltEUQwQRdFfFMW3T9uWiKK45JxjPhBFMUQUxb6iKH7aSXo7nYsVuBBFkeeee44pU6YwePDg87a7AtP6u7K89yRUhYW8/dC7hL+3m5TC2sv+fHVJI+kHiuk72h2LmtNlxmSHLnOZ9I1wx8LOmLj1mdKT+5loqMvMxFj+2WcobWywu/NOAP44VMC93ycy++sD3LnyIC0qLdf1ceks+T2CrvmWT4+cKXAxb968dvbPP/+cnTt3UltbS0ZGBm1tbe22Fy5cqCfFZ5kxwI0P/fuTfNKXuSd2sttnGEv2ZfLFrZf3g5OwKRulkZLBk71h00GpaIW5HFEgc3koDRUMm+7L7h9OkH20Ar+BPmDlLhUQH3b/P362KSGBxpgYnJ59FqWFBQDfH8jBz8GcPu7WbEwqYqS/PSP95f74T8gO/W9crMDFY489xmOPPXaerSthZmTA2odGcdxfg/X/Pcn/mpP5vzRDGlrVWFwiE11FQQOnEssYMtkbM0tDaRFRoPwSVObKCAxz4fD2POI3ZuEzwAGFTzhk7pFWkF5kRacoipR9+hkGjo7Y3nYrANkVjaQU1vHytGDuDffl6YkBuNuayqtCL4GcnOtvdPUCF5fC086MSbdMxjwiggExG1A2N3E4t/qSn4vfkIWRqQEDJ3pBZQY0V8k5z2WuGIVSwbAZvlQVNXIqoVSaR28sg4pTF/1MY0wMzYcO4fDQgyhMpAVDcacjWsYFOSEIAj4O5hgqZXd1KeQW6qE4Pv44ivo6bszcR2LOP6cEKMmuJedYBYMmemFibnh2zlPO2SLzL+g12Al7DwsObspG4zlKMubGXPBYURQp/+RTDN3dsZk9W2dPzKnG3twIX4dru+jzlSI79B6Kad8+WE6axOzMaFJP5P/jsfF/ZmFiJtDfYhs0V0PaBrD1lUrMychcIYJCYPhMP+rKmzlxwhQsXSEnBrQayE+Qkrydpv6vHbSkpeHwyCMIRkY6e2JuFUO85aLPV4rs0Hswjo89irG6Ff9da1FptBc8pjC9moIT1Qwx+gmj3c/BZwMha4+00k/+Msn8S7z72ePiZ0Xi1hzUnpFwcjt8Nx1WTIBvrwONClGjoXzxYoz8/LCeOUP32bL6FnIrmwj1sdXjHXRPZIfegzHu1YuGiAlMyYghLTnrvP2iKBK/IQtz0zb6Gm+AWV9JoynXgTBC/7H1Mt0XQRAIm+VPQ3UrqYZ3S6PzoiPQby6UpsCJzdRt2kRbZiaOjz3WbmVzYo70zmeoj7wi9EqRo1x6OB5PPkZl9C7Kly6Frz9oty8vtYrizFoivXZh4DQQBt0m/ZGR6QA8Am3xCLLl0P4Ggp8/gpGJEZjaQNYexGNrKf+2BOOQYCwnTWz3uYScKkwMFfRxs9aP8G6MPELv4bgF9yI2YCRO+7bSVnC2TN2Z0bmVvTHBqh/AZ5QeVcr0VMJm+dFcr+JYfKu0pkGhhIDJ1GyPQZWfj9PjjyMo2ruhhJwqBnnadnj1rWsBucWuAQqn34wGgfKvvtLZso6UU55Xz9DhKpRiC3iN1KNCmZ6Ki681Pv0dpILSjSoAtK5DqUgywLRvIOajR7c7vrKhlbSiOobK8+f/CnnK5Rpg5PAQNvuO4Pr1f3J8zCxyTBww3V2GrYsZAeZ7AUGOOZfpNMJm+rH67YMc3ZHH8Ov9qU4oR92sxG3WIARBoKC6iTc3pRHobImZsQFaEab2d9W37G6JPEK/CN2pwMWlGBfsxF8DJ9MiKDn13iesXX+SmpImhs3wQ5G/H5z7SnObMjKdgIOHBb2HOJG0p4D64hoqf1qLuZsWc3MpC/cr61PYnlrK4t0ZLNp6gsFeNgS5WOlZdfeky47Q3zv4HieqTnToOYPsgnhu2HOXPK67Fbi4FMYGSv7vzgjicycw+tA22oKVVBiIuAdbwrYEGHS7viXK9HCGzfAj43A5B77ciWd1NY6zgyE/nsKaZvaeLOex8b0Z6mPLvvRy7gn3vfQJZS5I9/RQncyFClykpqby9ddfU1FRwfjx4xkzZgyfffaZbvtC+V+6EpEBjoR/+Sp7b27FQK1kn3kr0w5HEapqAm95/lymc7FxNiNwsB3pB1X4jJuO6Ugv+OtlDhxNQRTh+oFu+DlaENHb8dInk7koXdahX85IujNQqVQcO3aMgIAAAIqLi3FwcGDAgAEsWbIErVbL/fffz4MPPthuuzsgmlmQ22sGVhVZmKGmPu10bLpPhH6FyVwT+FdEkU4AeSFzCPCRVoXWn9iDi9UQeYl/ByHPof+NfypwsWHDBsLDwxk/fvwFt7s6qVFFNKsM6V2+hwezduBYFi3Nn5s76FuaTA9HXV5Oy+qV+BoXcjKtiRpFL0QTa+xL9xPmZycv8e8gZIf+Ny5W4AJg5syZ7N+/n59//vmC212ZthY1h7bl4OFYRe+QUvzy0vAtzKAl6AZ9S5O5BqhYugyxrY0RCyNRKgUObs6j3mcSY7VxjHVpg6gP4eRf+pbZ7emyUy764mIFLvbu3cvatWtpbW1l6tSp5213dZL3FtBcryLM8D1sbXOoMrWlJNmG0oXTkGfQZToTVVERNatXY3PjDdj08af/ODj8Vy6t425ghvA71+89U2BdgAf2gttAPart3sgO/W9crMDFmDFjGDNmzHm27kBrk4ojf+XhY5OJi6sW7s/A2vglKn/YRcaeI4wcEKxviTI9mDML2hweegiAQZO8SNlXwPEkKw4bPsirA+oQBs2Hn26ExG9h5mJ9yu3WyFMuf6O7F7i4EEd35tPapGaY8gspOZKJNY7PfEyllQPua3+Q6j/KyHQCrdnZ1K5bj80tN2PoKi0WMjE3ZOBEL8wrVFQ53IBw/VfgPQKCZ0LqetCo9Su6GyM79B5Oc30bSbvy6dVbhaNhDviPBUAwMiJjys24leVStW27fkXK9FgqPv8CwcgIh9Pvoc7gFOpIkyASWK45a+w9EVproTT5KqvsOcgOvYdzaHsu6jYNwzzjwcBUSo17Gp95N5Jn4UThR58iajQXP4mMzL+gJT2dui1bsJs/HwMHKZJKpdEiiiI7TpYRb6xGLGmh8OTpEoln1kPkxOpJcfdHdug9mIbqVlL2FRI43AXbqr+kfC0GZ6vCjAp0Zt3AqRgW5LLyjaXM+Xo/Sfk1+hMs06Mo/2wxCktL7O+9B4Af43IJemUbs76MZcm+LNR+ZphbGxH/Z5Y07WflJlXKypUd+r9Fdug9mENbcxC1IkMn2EFJilSw9xwMlQpC588m09oNn02/cKqohmd+T5Ln1GX+M81JSTTs3o39vfegtLamrkXFe1tPEOhsSVVjG42tal6YHkLoNF+KM2vJTZGKQuMzCnL3g/bCFbZk/hnZofdQasubSYspImSUG1b1hwERfMLPO+7e0f64P/0kro2VfGKVz6myBvZnVl59wTI9irJPP0VpZ4fd/PkAbD5WTEOrmkWz+7Hv2bEceXUiYX72BI90xcrBhPgNWYhaEbzDoaUGylL1ewPdFNmh91ASNmcjKAVCp/pIj7BKY3Afct5xgiAw8KZpmA4ciMemVVgIGnafKLv6gmV6DI1x8TQdiMP+gftRmEtL+mMzKnCxMqGfuzVKhYCJoVRyTmmgYNh0XyryG8g8Un620Io8j/6vkB16D6SyqIH0+BL6j/HA3MYYsvdJ8+eGJhc8XhAEHJ94Ak1pKQtqk9ibLjt0mX+HKIqUf/opBi4u2N5yi84Wl1XF8Iss8e89zAVbV3MObsxCa+UJNl6QE321pfcIZIfeAzm4MRtDIyWDxzlA2QkoSYaAyf/4GfPhYZgNH87ohM0UFldRVtdyldTK9CQa9u6l+ehRHB58EIWxMQCZ5Y1UNLQy3M/+gp9RKATCZvhSXdLEyfgS6DUBMnZBS93VlN4jkB36ReiuBS7KcuvIOlLOQNM1mHwVDL/eCgpD6D/vkp91euJxjOprmZUZQ1x21VVQK9OTELVayj/9DENvL2xuPJsj6ECW9E5mhP+FHTqA3yBHHL0sObgpG03fW0DdDH8+DEtHw09zZOd+mXTZpf8l77xD6/GOLXBhHByEy4svXvK47lzgIv7PLEwMWhhotR0Cp0lzkdM+BItL55k2HTgQ88hI5uzfy6bUucwc4HYVFMv0FOq2bqU1PR23Dz9EMDTU2eOyKnG1NsHLzuyinxUEgbBZfmz6PInjOQH07XMjpK4F+15QvBNiP4Pxr1yN2+jWyCP0C7B48WLee++9dgUusrKyuPfee5kzZw7AedsX49133+WZZ57hxIkT7f7dGRSdqiYvrYrBpr9hNHQezP4Gnj4OQ+667HM4PfE4FqpmrDf+1ikaZXomokpF+eLFGAcEYDV1is6u1YrEZ1Uy3M/+kilyvULscO1lTcLWHNQzl8OTafBIojRdePh70Kg6+za6PZc15BQEYTLwGaAEvhFFcdFFjhsKxAHzRFH8478Iu5yRdGdwsQIXAQEBrFixQufA/fz82m2f+exrr71GU1MTbW1t3HnnnaxatYoFCxZQW1ur+3dQUFCH6xZFkbg/szA319LPbDOEbPtX5zEJDqY8NILIo7spyS3CxVsepctcmpp161Dl5uHx1VcIirPjxLTiOioa2gjvdemc+4IgMHyWH+s+OkLyviIGTfKSdgyYBye3QtERuZj5JbjkCF0QBCXwJTAFCAFuEQQh5CLHvQd068Qg/1Tg4lIsW7aM5uZmbGxsaGhoICAggDFjxvDwww+3+3dnkJdaRXFGLaG+yRiYmIDLgH99LpuHHsJY3Ubm4q87UKFMT0Xb0kLFl19hOnAgFmPHtNu353QIbETA5RVRcetti2eIHYe359LWfDpJl28kIEDW3g7T3FO5nCmXYUCGKIpZoii2Ab8Csy5w3KPAGqBbx7z9U4GLy/nsokWLeP311/nhhx84evQoAwZIjvXcf3c0olYk7s9MrBxMCNb8IuXEUP77+f4+wwcQ5z8Mq23raSwsIqWwFo1WXj0qc2GqV/2KurQUxyefRBAEimub+fFADmlFdfyakM9wPzucLC8cMnshhs/yo6VRRdLufMlgZgcu/SBrXyfdQc/hchy6O5B/znbBaZsOQRDcgRuAJf90IkEQHhAEIVEQhMTy8vIr1XpVOFPgIj4+npSUFN544w0AKisrWbhwIUeOHOHdd989bxtg1qxZ3HXXXfzvf/9j27ZtJCUlMXDgQIB2/+5oMo+UU5HfwNDxtiirT/7nGqFKhYDi7vsQtSLLH3iF6Z/H8PDPh+WUADLnoWlooHLpUsxHjsQ8bBgtKg1zlx7glT9Tmbo4msKaZhZE+l/ROZ28rfAb6MjRHXm0NJyeN/eLhPx4aGvqhLvoOVzOMO5CbzL+/s3+FHhOFEXNP734EEVxGbAMIDQ0tEt6h4sVuLC3t2fJkva/V3/fnjFjBjNmzNBtT558Nvb7iSee6Fihp9FqtBzcmIWtixkBNscko+9/L/p8+/Uj+WvzBCbG76R21s2sTC3haH4Ng7xs//O5ZXoOVd99j6amBscnnwSkJf75Vc18PHcAtc0qPGzNGBvodMXnHTbTl6ykcg7/lcvIG3uBz2jY/zkUHAS/MR18Fz2HyxmhFwCe52x7AEV/OyYU+FUQhBxgDvCVIAjXd4TAq013K3CRHl9KdUkTYbP8UORGgYkNOPf7z+c1MlBw3aIXURoacNepnZgYKvj9UMF/FyzTY1BXV1O1ciWWEydi2q8vAOuOFOJtb8YNg9y5e5QvE0Oc/9W57d0sCBjqTPKeAhprW6UCGIISsuUVpP/E5Tj0BKC3IAi+giAYATcDG849QBRFX1EUfURR9AH+AB4SRXF9R4uVaY9GpSVhUzaOXpb49beHjB1SAQtFx0SjGjo7YXvrrTRt3sSN9ir2niiTp11kdFQuW462uRnHxx8DoLlNw8HsKiaFOF8yRPFyGDrdF61G5NDWXDC2BLdBckqAS3DJb74oimrgEaTolePAb6IopgqCsFAQhIWdLVDm4qTFFlFf1cLwWX4IRYegofSSS/yvFPv770NhYsK0Q5soqm0hp1Kew5QBVUkJ1T//jPXMmRj36gXAodxq2jRaRl5GiOLlYONkRtAoV1KjC6mrbJZSAuQfhLq/TxDInOGyhnKiKG4RRTFAFEV/URTfPm1bIorieS9BRVG867/GoMtcmrYWNQmbsnCzr8azZAnsegNMrCFwaodex8DODts778AuIRr/mkJiMyo69Pwy3ZOKr5cgiiIOj5wNw43NrMBAITDMx67DrjN0qg+CIJCwOQf63QSI0qrR4xshbQPIT4ztkFeKdlOO7synuUHNCHERQsxH0qNo5HNgYtXh17K/+24UVlbcn/EX+zNlh36t05abS82aNdjedBNGHh46+/6MCgZ62mBu3HEpMixsTeg72p30A8VUq12lnETxS2D17fDbfIj7qsOu1ROQHXo3pKmujaN/5eBnEodL5CR46jg8FAcjOmfRktLKCvt772VAQSrlcYlyTPo1TvlnnyEYGODw4NkZ15qmNo4V1hLeu2OmW85l8GRvlEZKDm7KhhmfwcwvYP46acFR9Megau7wa3ZXZIfeDUnckoNapWW4xY8Qeg9YOoNTcKde027+7aitbLjp8AbSCms79VoyXZfm5GTqtmzF7u67MHA8m/DtQGYloshlLfG/UsysjBgwzoOMxDIqStUweD74j4NRj0FTBWTu6fBrdldkh97NqClrIjWqkBCXE9g6GICd31W5rsLMDOuFD9KvMovUdVuuyjVluhaiKFL24UcobW2xv/fedvv2pJdhYWzAAE+bTrn2oIleGJsZEP9n1lmjz2gwtoJ0uT+eQXbo3Yz4DVkoDASGKr6WRikdEB52uXjOv4UKG2ecfvkGrUrOfHet0RgdTVN8PA4PPYTSwgK1RktRTTN1LSo2HStmaj8XDJWd41KMzQwZNMmLnORKSrJOPyEaGElhuhm75Jejp5Ed+kXoigUuynLryEgsY+BQBeaafMmhX0UEQ0Na716Ia00x7zz+IaFv7eStTWlybPo1gKjRUPbBhxh6eWE7by6tag03fLWfkYt2M/StnTSrNNwxwqdTNfQf64mppSFx547S/cdBfRFUnOzUa3cXumzFhujfTlKR39Ch53TwtCBibsAlj+uKBS5EUWT/2kxMLAwZ5LwfshTgO/qq6xh731xi16xiQtyfnLhnBN/EZBMZ6EhE70sX0JDpvtT+uYHWU6dw/+RjBCMjfozOIrmwlvsjfKlsaGN8sDN93a07VYOhsZIhU3yI+e0UeamVePWxB7+x0s7MPeDYfVZ4dxZd1qHrk8WLF7N8+fJ2BS7Wrl3Ljh07KCsr4+GHH6apqYnNmzfrtidNmnTBc51J5HXfffexbt063b+vNCd6floVhenVhM/tjVHma+A2GEyvfl4VpVJB6KLXyL31Nj4U0phuEcB3sTmyQ+/BaFtaKF+8GJN+/bA8nZ9o47FiBnja8NK08zJpdyp9R7tzbE8BsWsy8Ai2Q2HrDXb+kLkbhsvrHLusQ7+ckXRncLECF3PnzmXu3LlUV1fzzDPPsGLFCq6//nrd9qRJkzqtwIWoFdm/TkqP23ewIUQdhrH6KQACYDZ4MJYTJ1C7YgXznv2C5SkVNLWpMTPqst1J5j9Q/dNPqEtKcHvvPQRBoKKhlWMFNTw54ep/R5UGCkZc78/25SmcOFBMyCg3aR796CpQt0nz6tcw8hz637hUgYu33nqrXZGKc7c7q8DFyYRSKgsaCJvph/LkRkDs8CX+V4rjk0+hbW1lQuIm2jRaDmRW6lWPTOegrq6mYukyLCIjMQ+TqgVFnSxHFPlXWRQ7Av/Bjrj4WRG/IYu2FrU0j65qlPO80IVH6Pri3AIXjY2N3HjjjTzwwAOIosjzzz/PlClTGDx48HnbZz775ZdfYmxsDMCePXv+c4ELdZuGuPWncLBtpLe4EQ4sBvdQKeG/HjH288Vm7k3U/P4H/hMD2Ztezvjgf5dZT6brUrl0GdrGRhyffkpn25NejoOFMX3cOn5V8uUgCAKj5vRmzfuHOLozn2HXjQczB4j9VJqGbKkB3zEdlqSuOyE79L9xpsDFvHnz2tk///xzdu7cSW1tLRkZGbS1tbXbXrhwoa7AhaenJ+PGjePEiRNEREi5yZOSknT/vhKSduXTUK1ivO0ihK0poDSCWV9d1XDFi+H48MPU/rmBx7J28pGr56U/INOtaCsolBJw3XA9JqenINUaLVEny5kY4oxCob8+6OJnjf9gJ478lUufcDfMI56G7S/A8tMvSYc/DJPf0Zs+vSGKol7+DBkyRPw7aWlp59muNpGRkeKJEyf0cu2/339DTYu49JFd4ubH3hHFhBWiWJkpinXFetF2Mcq++EJMCwwSr7v/CzGnokHfcmQ6kIInnxSPDxgothWf7XMJ2ZWi93ObxE1JRXpUJlFT1ih+9dBucdcPaaKo1YriyR2imLJWFP+4TxRft5W+Lz0QIFG8iF+99p5JLkFXKnBxcGM2GrWWkba/wcDbpFWhli76ltUO+7vvBnt77k/ZSFR6ty4nK3MOTYcPU7dlK/b33ouhy9k+tyOtFEOl0Ck5W64Ua0cz+o3x4Pj+YioKG6H3BOhzA0z8PxC1cOx3fUu86sgOvYtSUdDA8dgi+lnvwyYoBAyM9S3pgijMzHB58glCqnLJ++NPfcuR6QBErZbSd97FwNkZ+3vv0dm1WpEtKcWM9HfA2tRQjwrPEjrVB2NTAw6szThrtHKTCqWnrtOfMD0hO/QuiCiKxP5xCiMTgVCjFRBwnb4l/SM2N95Inac/EXt+5UROOXUtKrRyRsZuS93GjbSkpOD09FMozMworm3mYHYVqxPzya9q5sbB7pc+yVXCxNyQ0Kk+5KVVkZd6TqRV0HQoPw5V2foTpwe6nEMXr9Fl5Ofed25KJQUnqhkako+JohF6X3jRUldBUCjwfOUlHJtr+enJt+j/+l9MXRxNTVObvqXJXCHapibKPvoYk379sJo+nbSiOsZ9uI+5Sw/wwtpkQlytmN7fTd8y29Ev0gMrBxNi12Sg1Wgl45lB0Mnt+hOmB7qUQzcxMaGysvKac+qiKFJZWYmJiQkajZb9azKwcTajL7+A68AuN29+IdxGj0AbOZ55GXt5aoA1J0vr+Xpvpr5lyVwhld+sQF1WhvMLLyAoFHyw/QQmhgqW3D6Ed2/sxy/3h6HUY3TLhVAaKhh5Yy+qihpJjT5dns7eHxwC4OQ2/Yq7ynSpsEUPDw8KCgooLy/Xt5SrjomJCR4eHqRGFVFd0sTUe3xQbomDyP/pW9pl0/uVF8iaGs3sQ3+SGXorP8fn8eTEAEwMlfqWJnMZqIqLqfz2W6ymTsVs8CBK61rYk17O4+N7M7lv1x5U+A1yxD3QhvgNWfQOdcbEwlAapcctgZa6Tqnk1RXpUg7d0NAQX19ffcvQGy0NKg5uysI90BYfdgMiBM/Qt6zLxsjDHbt77qZyyVLmjZ7Kn61qok6WM6lP13YGMhJlH38CoojT6UVEu45LUUtT+7nqU9ZlIQgCEXMDWP12AvEbs4i8JRACpsD+zyFrD4TM0rfEq0KXmnK51onfkEVbs4aIOX4Ih1aCYzA499W3rCvC4f77MXB0xOWHr7AxUbI9tVTfkmQug+akJOo2bsTu7rswdJdeeu46XoqHrSkBzhZ6Vnd52Ltb0He0O6lRhVQUNIBnGJjYQMoafUu7asgOvYtQnldPSlQB/Sz+wn7VIChNkaZbusCK0CtBYW6O49NP0ZqSwv0t6exJL5NrkHZxRI2GkjffwsDREYf77weguU1DTEYFE4KdEbpRHxw2wxdjM0NifjuJqFDC0Psg7U/4ZR58P6PHl6uTHXoXQNSKRH1/CFNFLcMCToLfGJi8SFok0Q2xnjkTk/79Gb3vd1pq6jiSV61vSTL/QM3vf0hhis89h8LcHIDYjApa1VrGB+snAde/xcTckLBZfhSerCHzcDlEPA19Z0PFKenPb3dCXbG+ZXYaskPvAqQfLKGkUMsIm98xnv89zPkWhj/Y7UbnZxAUClxefgmDmirmp//FzuPyCtKuirq6mvJPPsFs6FCspk3V2benlmBpbECYr70e1f07QsLdsPewIHbNKVQYS9+nxw7DnZugrR4SvtG3xE5Dduh6pq1Zzf41GTgbZxI02ByMu8d85aUw7d8fm5tuYmZmDKmxh/UtR+YilH/8CZqGBpxfeVk3tdKq1rA9tYSJIc4YGXQ/F6FQCIye15uGqlaObM89u8OhF/SaAEd/Aa1WfwI7ke73v9XDOLg5m+Z6FaMtvkbo07PexDs++QRacwtm7PmJuMwKfcuR+RvNycnU/PEHdvPnYxIQoEvw9FtiAXUtam4c7KFvif8at9629A514vD2PGrKms7u6HeTVIO08JD+xHUiskPXI1VFjSTvLiDEKxcns2LwH69vSR2Kga0tLs88Rb/KbH56cwk3fhXLwz8fprZJpW9p1zyiRkPJ/72B0sEeh0ceJr2knvEf7SPw5W28viGVYb52jOrV/aZbzmXUnN4oDQSifj15drFi70mgMITjG/QrrpOQHbqeELUie385gaGJkuGKz6SOZmSmb1kdjsPcmyCkL3cd24iFqoVtqSW8uTlN37KueWr+WENLSgrO//sfgpk5z/6RRHVTG/NHeHNfuC9f3za4W0W3XAhzG2PCZvmTn1ZFxqHT73FMbcAvEo5vhB64Il126Hri+IFiijNqGRkpYtqSDSEz9S2pUxAUCnzf+j/MWxp4v3Y/94b7suZwAQXVTZf+sEynoK6qovzjjzELDcVq+nSO5NdwrKCW5yYH8cr0EF6YGoy9RdfM7nml9I10x9HLkpjfTtHarJaMwTOgOhtKU/UrrhOQHboeaKprY/+aDFx7WRMs/g5GFnqvEdqZmISEYHvLLVT/+iu32UmOfM2hQj2runYpXbQITVMTLq+9iiAIbD5WjJGBgmn9u/6K0CtFoRAYc1sgzfVtxP+ZJRkDpwGCNErvYVyWQxcEYbIgCOmCIGQIgvD8BfbfJgjCsdN/9guCcOXFM68hYtecQtWqYcw0c4SU36V4cyNzfcvqVBwffwylnR3a999huLc1G5IKr7kkbF2BhphY6jZsxOH++zDu3VvKcZ5cTGSAI5YmXSPHeUfj5G1F3zEeJO8roCy3DiwcwSccjvwIbT3rSfGSDl0QBCXwJTAFCAFuEQQh5G+HZQORoij2B94ElnW00J5CfkoZJ+NLGWyyGrvVYaAwgDHn/Ub2OJRWVri8/DItaWncUxRHZnkjJ0rq9S3rmkLb1ETJ669j5OuL/YIFABzOq6akroVp3SBfy38hbKYfZlZG7P05XcrVP+Z5qCuExYPg8yEQ84m+JXYIlzNCHwZkiKKYJYpiG/Ar0C6+ThTF/aIonlkOGAd033inTkTdpmHf94lYK4sYMs4exr4E92wH62ujuSyvm4TFhPF4/PkTHk2VbEwq0reka4ryL79EVVCAy/+9jsJYmiPfdHq6pbutCL1SjE0NCL+pN+V59RzbnS+N0G/6HjxCwdQOdr4upQjo5lyOQ3cH8s/ZLjhtuxj3AlsvtEMQhAcEQUgUBCHxWkyRm7glh9p6EyIDYjGY9raUq8WleyXf+i8IgoDLK6+gMDTkpRPr2ZRUJE+7XCVa0tKo+u57bG6ag/mwYQCoNFo29/DplnPpNcQJn372xP+ZJcWm97kebv4Z7t4q5U6P+qDbR75cjkO/UOzSBe9aEISxSA79uQvtF0VxmSiKoaIohjo6Ol6+yh5AeV49h7fnEmiyB8/ICH3L0RuGzs44PfssPnnHCTqyl6SCWn1L6vGIajXFr7yK0tYWp2ee0dl3ppVSXt/K3FBPPaq7egiCQOStQSgMFOz58QTimaRxSgMp1UZJMhQd0a/I/8jlOPQC4Nz/cQ/gvGdlQRD6A98As0RRrPz7/msZjVrLru/TMDVqIdzul26V47wzsLlpDkaDh/BA6kY+Xb2f7aklJOXX6FtWj6Xqu+9oSU3F5aUXUZlZsOpgHh/9lc5rG1LxczBnbOC1M7iysDVm1JxeFJ2qISXqnEirPjeC0hiSftWfuA7gchx6AtBbEARfQRCMgJuBdsusBEHwAtYC80VRPNnxMrs3iVtyqCxsZKz115j0GQfGlvqWpFcEhQKPt9/EFA2jN69kwQ+JzPoyls93ndK3tB5H66lTlH+2GMuJE7G47joeXXWEF9Ym8/nuDETg43kDMVBeW9HLwSNd8Qy2Zf+6TOoqmiWjqQ0ETYWUP0DdfWvhXvJ/UhRFNfAIsB04DvwmimKqIAgLBUFYePqwVwF74CtBEI4KgpDYaYq7GeV59RzalktgkAofRRQMmKdvSV0CY19fXJ54nBHFKWwMrGdaf1c+2XmS7IpGfUvrMYgqFUXPv4DCwgKX118jLruKHWmlPHtdIDmLppHw0gQGetroW+ZVRxAExtwehADs/fnE2fc4/W+GpkrI2KlXff+Fy/ppFkVxiyiKAaIo+oui+PZp2xJRFJec/vd9oijaiqI48PSf0M4U3V3QTbVYGhJu/R1YuIDvGD2r6jrY3XUXpkOGYPTVx7w81A4DhYKVsdn6ltVjqFi2TJpqef11DOzt+SOxAEsTA+4Nv3bLPJ7Byt6UkTf6k3+8muOxp/Oj9xoPZg6QtEq/4v4D19az1lVGN9UyScAkdwsMu096ASMDgKBU4rboXUStlrZ3/o8pfZxYf6SQFpVG39K6PS1paVR8vQSr6dOxum4S9S0qtqQUM2OAm1y0+zR9ItxxD7Qh5vdT1JY3gdJQysZ4clu3LYIhO/ROoiglj0Nbsggy24tP9FSw8Yah9+tbVpfDyNMT5+efo+lAHPNLE6lrUbMjTa5D+l/QtrVR9PwLUrbLl18CYPOxYlpUWm4acm2sebgcBIXA+DtDUCgFdnybhlajlUrWKQxg+Vj4+SbY8w5ou88AQ3bonUBrs5qd3xzBUllGxBg1RD4Pd26UXrzInIfNTTdhHjkai++XMIha/jhUoG9J3ZryTz6l9eRJXN58A6WNDQC/HyrA39H8mpwz/ycs7UyIvDWQ0uw6ErfmSkUwbl0NLv2gvhj2vQd7F+lb5mUjO/ROIOrbeBpazJg4rhqjme/A2BfA1lvfsrosgiDg+uabKIyN+d+hXziQXkxJbYu+ZXVLGqKjqVq5Ettbb8VyzBgAjhfXcSi3mnlDPbt9StzOoHeoM4FhLiRuyaEkqxZ8R8Ntv8PCGOg7B/Z/Dg3dYyGk7NA7mPT4Ek4mtzLUcg0uU2/Vt5xug6GTE65vv4VNfiZ3pmzhp7jcS39Iph3q8nKKnn8B4969cfrfswCIoshnO09hYqi4ZhYQ/RtG3xyAha0xO75Npa1Ffc6OZ0HdDMdW60/cFSA79A6krqKZqFXpuBqfZMgwDZja6ltSt8JywgRsb7uNGzOjSFi1gft/SOT5NcfIKGvQt7Quj6jVUvT8C2gbGnD/+CN2ZdUy68tYwt/bw7bUEh4d1xsbMyN9y+yyGJkaMOHuEOorW4hadU6FI6cg8BgqZWbsBmkBZIfeQahVGrYtSwGtigmWH6EIna9vSd0Sp/89i2FgEM8n/UZVdj7rjxYyd+kBKhpa9S2tS1P13fc0xsbi/MILnDJz4sGfDlHfrGKglw2vTg/hwUh/fUvs8rj1smHodF/S40tIizlnMfyg26H8BBR0/eU1skPvIGJ+z6A8r54JXuuwcrIEn2s3X8t/QWFsjOenH2MqaPno1BrWLxhOfYtKXkX6DzQdPkzZxx9jOXEiNvPm8vGOk1ibGrL2oZF8eetg7gn3RaGQ584vh9ApPniG2BG9+hTleafTO/edDcZWEL9Ev+IuA9mhdwDp8SWkRhUyaKiIb8NPMPRekF8+/WuMfX1xff01mhMPYffrCmYMcOP3QwXUNsvFpf+OurycwsefwNDdDde33yKvqok96WXcNtxbnmL5FwgKgYl3h2BiYci25SlS2TpjSxh8B6Sug8pMfUv8R2SH/h+pPJnN3u+P4WaWxfCi28HOH0Lv1besbo/1zJnY3DyPym9WcK8qk6Y2Db8l5F/6g9cQokpF4ZNPoamvx2Px5yitrPjxQC5KQeC2MC99y+u2mFoacd19faivbGH398el+fQRD4OxBfxwPWx7ARJXgqrrRWLJDv0/0NrQwrYvD2JEPZMGxqMYcjvctQkMTfQtrUfg8uKLmA4ejPLDt5lu2cT3B3LQaLv+i6mrRdlHH9OUmIjrm29gEhhAQ6ua1Yn5TO7rgrOV3Af/C669bBhxgz9ZR8s5vD0XrNzgltVg6QKHvoNNT8C6BfqWeR6yQ/+XaDVa/lq8l7pWG6673hDzO5bD1A+k/3iZDkEwMsLjs09RWlrywF9fU1Nayc7j8ipSgNqNG6n67jtsb7sN6xlSOubVCfnUt6i5L8JPz+p6BgMneNJ7qDNx67PIOloO3iPgvh3wYhGMexnS1ne5RF6yQ/+XxP6RQV6eEZFuf+I2Ybq+5fRYDBwd8Vj8GUZVFbx+dBUfbEmlpLYFtUarb2l6o+nwYYpffAmzoUNxfPZZMsrq2X2ilM92nmS4n528GrSDEASBcfODcPK2ZMfKNCoKGs7sgJGPgbUn7PtAvyL/huzQ/wWp0YUc21PAALMNhMwYBQq5GTsT04EDcXn9NfoUHWfa7p8Y/s5OAl7eyjO/J6G6xhx7W34+BQ8/gqGbG+bvfcQN3yQw4eMo7vkuESMDJYtu7K9viT0KAyMlUx/sj7GJki1fHaO5/nSudANjGPEI5MdBXpx+RZ6D7ImukPzjVUStOom3dSYj3XdI2dlkOh2b2bOxf+ABJmfH8bUihXlDvfjjUAGf7Lh26qlo6urIX/ggolaLx5Kv+SS+hOPFdbw5qw/f3zOMXU9H4uNgrm+ZPQ5zG2OmPNifpvo2Nn91DFXr6WRdg+dLBaZjPtWrvnORHfoVUJZbx9Ylydjaaphk/AqKiCfAQA4Nu1o4PvE4VtOn47Pue543zOWGQe6siMmmuLZZ39I6HW1rKwWPPkZbbi4en31GiZUTqxPyuWWYF/NH+BAZ4Ii1ac8v9KwvnH2smHhPCGU5dWz/JgWNRgtG5jDsATi5FYqP6VsiIDv0y6bmRBqbPorBhCpmGD2CkVsvGHynvmVdUwgKBa7vvI3Z0KEUvfgij9nWoBVFPt+doW9pnYqoVlP49NM0xcfj9vZbmA8P4+MdJzFQCjwytpe+5V0z+A9yYvQtgeQmV7L3p9OVjsIWgLkj/DIXNj4BMZ9Ac43eNMoO/TJozM1gwxfJiOo2Zvr/hHn/MXDr7/LoXA8ojIzw+OJzjH18aH3uKR52aOT3xHwKa3rmKF0URYpffY2GnbtwfvFFrGfN4nhxHRuSirh7lC9OcnjiVaXvaHeGTvPhxIES4tZnIpraSpkZrdzh+AbY+Tp8Oxla6vSiT3bol6CxtpX1nx6jWWPJ9AUB2Dy2Bm5YAlau+pZ2zaK0tsbr2xUYOjtz3U/v06s6n0974Fy6KIqUvvsutWvX4vDQQ9jdMR+NVuT/NqZiaWzAgtFyeKI+GDrdlz6j3Tm8PY+Dm7LBbRDcvwv+lwW3r5Xyvux7Ty/a5Hpo/0BjbSvr39tPQ4sxM8bn4Dxgqr4lyZzGwNERr+9Wknvb7SyK+4YntQJPakUCnC0J9bFlqI+dviX+J0RRpPTtd6j+6Sds5s9nz4hZHPw9iczyBo7k1fDe7H7y0n49IQgCkTcHoFVrSdycgwAMm3H6x7XXeOllafwSCL0H7K9uUjR5hH4RGmtbWf/xYRpqWpnh/hVuM+X58q6GoYsLXt9/h5mVBZ/ELyNnXxzvbTvBTUsOsCKm+xabFrVaSv7v/6j+6Sfs7rqLxUHT+N+aZPadLKe5TcObs/owb6i8tF+fCAqBsbcHETTSlYTNORzcmHU25e7Yl8HABHa8etV1yQ79AtSUNrHm/UM0VDYyw+Z13GbMl95oy3Q5jDw88Pn5J8yd7Hk3bhkHJ1txXR9n3tlynOPF+pnH/C+IKhXFL79Cza+rsb/vXo7NuJNVCQUsiPTj4Ivj2fbEaOaP8NG3TBkkpz7uHKcevfoUWq0Ils4Q/iSc2AQZu66uJlFPSdtDQ0PFxMSul1+4ZNdGNq9XIIhaplm/jnOQN9z2h7x4qIujLi8n7977aMvOxvqtd5h2zBg/B3N+WzCi26SO1TQ0UvjEEzTGxODw8MOYP7CQiZ9EYW6sZPNjERgq5T7YFRG1IvvXZnB0Zz7+g52YeHcISlphaSTU5oP/OCkHzMDbwH3wf76eIAiHRFEMvdA+uYecQ/aG9fz5hxIjGrmx/xqcx82CeT/KzrwbYODoiPcP32PSty81zz3LRy2HSMyp4ueDefqWdlmoysrIvWM+jQcO4PLmGzg++gif78mgsKaZN2f1lZ15F0ZQCIya05uRs3uRebiMjZ8fpaXVAO5YD4FToeIUJP0KKybBye2dq0UeoUu/sIlrDnNwVy2OZiVMe3kW5nbyFEt3RNvaSvErr1C3YSMpwcN5K/gG3pwXyih/e+zMjbpkkeSmI0cofOJJNPX1GL+5iKLAQRzKreaTnSeZM9iDD24aoG+JMpdJenwJu388jrm1MVMf7IeDh6W0o7kGvp8ONfnwYCxYe/zra/zTCL3bO/QmVRNmhmYXtJsamJ73BW5SNWGsNEapUALQ0qhi9/epZB+rord5LONeuBMDh7MvnFrULQiCgLHSuN15VBoVGlGDiUH7OGCtqKVV04qpgWk7uyiKNKubr1iriYEJCkFxnt1IaYSBon2QUrO6GaWgxEjZPvqhTdOGiHjePai1alRa1RVr/S/tfa7dUGGIobL96sZWjVRq7rz21qrQaC/c3i3qlnaaRFGkYulSKj79jEJ7D14beCuFlk7YmhnyyLje3BLmfMXtbag0xFDRXus/tbdW1J6nVaPV0KZt07W3KIpU//wLpe8tQunszKbZT7A496ymcUFOfDg3CDszS/5OV2rvM/fSkf37Qu3dom5BISjOa++LfRf/3t6Xq/VK7+Hv38XS7Dq2LD1Ga6OasbcHERjmIu2ozKR1SQS4D8L4jo3/+sm/x065fH30a8J+CWNF8op29g2ZGwj7JYw34t5oZ08oSSD813Ae3PkgoihSkFrK6jf2k5NcQaLXH3w9bB8aW2fd8TUtNcxYP4PJayZT3lSus6s0Km7bchujV48mvSq93TWei3qOEb+MILogup3908OfEvZLGL+l/9bO/vPxnwn7JYxPDn/Szr4nbw8jVo3g2X3PtrOnVaYx5rcx3LH1DtTas9XJSxpLmLxmMjPXz6S+rV5nb1I1MXvDbMb+Npb8urMFIrSilgd2PMCoVaM4XHq43TXeinuL4b8MZ0vWlnb2JUlLCPsljOXHlrezb8zcSNgvYfzfgf9rZz/T3gt2LuDcgUNObQ4T/pjATRtvokV9tkhAbWstM9fN5Lo/rqOsqUxnV2lVzN8yn9GrR3Oi6kS7azwf/TwjV40kqiBKZxMEgV+GtfLuTQqc2ypZHvs5i+2K6etmxaLYFYT9EsZHiR+1O8++/H2MXDWSp/c+3c5+vPI4Y34bw/wt89u1d1lTGVPWTGHGuhnUtZ19+dqsbuamjTcx7rdx5NbltmvvBTsWMHLVSBJLElFXVVH42OOUvvUWxSHO3DGnhK+qj3LXSB9W3T+cqGfHMnzwESJ/H8nSpKXtNG3O2szwX4bz6v72URSHSg8R/ms4D+x4oF1759XlMeGPCczeOJtm9dkFWLWttcxaP4tJf0yitPFsWmKVVsUdW+5g9OrRpFWmtbvGSzEvMWLVCPbm721n/+LoF4T9EsYvx39pZ199YjVhv4TxQWL7rIRRBVGMXDWSJ/Y80c6eXpXOmN/GcPuW21Fpz1aoKm8qZ8raKUxfN53a1tp27T1301zG/TaOnNocnV0URRbuXMjIVSNJKElod41FBxcR9ksYGzI3tLN/k/wNYb+E8XXS1+3s27K3MfyX4bwc+3I7+5GyI0T8GsH9f92PVjybJK7Nvpafgt+k3DyfnSvT2PFtKq3NauosHZnl68ckbTYl0Z0Tp95tHXqrppWVqSsBWJGyApVG+s8XRVH3BVh7am07R7z82HJUWhUJBYfY8Po7/Pl5MgYNeVT5vk2iezQZLaXsztutO/6PU39Q0lhCRXMFq9NX6+y78nZxvOo4zepmvkv9Tmc/UXWCbTnb0Igalh1bprNXtVTx8/GfAVh6bKnuP18ranXO8ee0n9s54uXJy9GKWv7K/Yus2iyd/duUb2lWN5NckUxc8dksb6tOrKKqpYrChkI2Z23W2TdnbyanLof6tnp+SPtBZ48vjiehJAGVVsU3yd/o7AX1Bfxx6g9ERJYeO+tIVBqV7l5XpqykTdOm23fmXteeWtvOEa9IXoFKqyK+OJ6k8iSd/ce0H6lvqyezNpOdeWfzSa85tYaixiIqWyr59cSvOvuevD2kVqbSrG7m25Rvdfb0qnS2Zm9FI2raaa1uqeantJ840kvBaw9YYta3D72//Yi3Dv+Ap5kUdfDz8V/aOeJlycvQiBp25u0ks+ZsmbGVKStpVjeTWpnK/qL9OvuvJ36lsqWSosYiNmVu0tm3Zm8lqzaLelU9P6b9qLMnlCQQXxKPWqtm90+LyJo+g4a9ezF59H6emlJGo6mAg0cUr04PZoS/Pa42hqxMkfr3ytSVupH0mfYWEVmfsZ6SxpLz2vtgyUGOlh/V2X9I+4H6tnqya7PZmXu2vddnrKewoZCqlipWnVils+/L30dKZQrN6madBoBT1afYlLUJraht9yNT21qru9flycvRaKXkVaIosixZ6hurTqxq54iXH1uORtSwJ38PJ6vPLgpbmSq1d1plGrGFsWfbO/1XKporKG4sZmPmRp19W/Y2MmoyqFfV833a9zp7YmkiccVxqLXqdgOQksYS3aBqadJS3Q+fSqvS9a3vUr5rN9BYemwpIiIbMjdQ3FDcrr3btG0klia2GxT9dPwnyiji96APsRzZwqnEMn59M54/dm2lUFVLlVLJytbOqb7VbR16Ykkizepmbg68mfq2et2XLacuh7z6PG4NuhWtqGV7jvQSolHVSEJJAreJt3Dz0RcoKB1BH78ybrpDwS5PNRO9J+Jq7tqus+zN30tf+76EOoeyK+9s+NHu/N3Ymdgx038m+/L36X5M/sr5C6Wg5I6QOzhaflQ36tmTt4dWTSt3972bsqYynXNLrUilsqWSW4JuoU3bpvsxqWiuILkimZsDb0YhKHQjZZVWRWxhLDP9Z2JjbNNO6778fQx3HU6QXVA7+568PXhaejLeazy783frfkx25+3G1MCUmwNv5kDxAd2Pya68XWhFLXf3uZus2iwyqqU8KQmlCTSqGqX2VtXrvmy5dbnk1OVwW/BtiIhsy94GSE8GB0sOMidgDiZKEzZlSU5PFEX2FuxlvNd43C3cdfYz9xBiH0KYS1i7H9bd+buxNbbl+l7XE1UQpfsx2ZG7A4Wg4I6QOzhWfkzn3Pbm76VF08Ldfe7mlFEVFe89itOzz9AYG8sH39QwPc4PraaNxQfWkl/VRHFDOcnlUnsrBaXuB1GtVRNTGMMMvxnYGtu2c9x78vcQ5hJGsF3wee3tbuHORO+J7M4729578vfgWWvIZ9tcmPFNKloHezRLvuMzPyUaQcRWNYlGsZiMWqm9D5UdokHVwLzAeTSqGokpiAEgvy6frNosbg26FaBde8cXxzO792xMDUx1WkVRZF/BPsZ5jsPDwqOd1t15uwmyC2K46/B2/XtP/h6sja25sfeNRBVE6X5MzrT3nSF3klKZonNu+wr20axu5u6+d1PRXMHhMsm5nag6QVlTGbcE3YJaq9Zdo7qlmqTyJOYFzkMpKHX9W6PVEF0QzTS/adiZ2LXrG3vz9xLqHEof+z5szNrYTqubuRuTvCexO2+37sdkd95ujBRG3Bp0KwklCbofk115u1CLau7uezd59XmkV0tP2EfLjlLfVs/NgTfTpG4iplBq78KGQjJqMnTtvTVnKyA9GcQVx3Fj7xsxMzBr17/3FexjjMcYvKw9iXJey43PDkZpoKD5T0duyHmEdwd/wH3h7Z9mO4pu69D35u/FRGnC44Mfx8LQgj35ewB0j9539rmTXja9dPbdBw8wOfUBLOOGY6BsJnbID0Q+ewuZvQZT0VZDpEck473Gc7DkIM3qZmpaajhWfowIjwjGeY0joyaD/Pp8NFoN+4v2E+4ezkTvidSr6kksld4FxBTGMMBxADf2vhGQOjpAbFEsTmZOPNDvAQwUBjpnta9gHwpBwUMDHsLZzFn3GHtmumZ2wGwGOg7U2Y+WHaVB1cA4z3FEekQSUxiDRquhsKGQzNpMItwlrckVyVS1VNGibuFgyUGdvaypjOOVUo3EmMIYhroMZarfVNRatc5BxxTG0MumF/ND5kvtXLBX167GSmMeH/w4lkaW57X37cG3E2AboLPHFceh0qqY7DOZkW7SlIgoipysPklZU5muvROKE2hSNVHbWsvR8qOEu4cz1mssmbWZ5NXloRW17C/cz0j3kUz0nkijqpHEkrPt3d+hP7MDZuv6xJn2djR1ZMGABRgpjNhdtBf7e+8l/r2byXAVuGPPST5eBhnrfiXivd2M+WIJIiLjPaYz2Hmw7v/taNlR6lX1jPMaxxjPMcQUxqDWqiluKCajJoMIjwjGe40npTKFiuYKWjWtxJfEE+EewVjPsZQ3l5NakYq6qgqr5et5f1kLzicq+HmMguuHhjJtYzHbs/ZhpHVlxQ1Pt7uHffn7MFIY8fjgx7EysmJ3vtRnogrPtneQXZCuvQ+WHKRN28Z1Ptcxym0U+wr2IYoiGTUZlDSWMNpjNBO8J5BQKrV3XVsdSeVJOq05dTlk12ajFbXEFMYw0k1q7yZ1k27KIrYwlr4OfZkTMAdAd+2YwhgcTB1Y2H8hxkrjdv1bQGBB/wW4W7izJ+/s8SIi1/e6nlDnUF17H6s4Rl1bHWM9xzLWcyyxhbGotCpKG0s5WX2SCI8IJnhPIK0yjYrmCto0bcQVx+n+H6paqkiuSNZdY6jLUKb7TUctqnX9NKYwBh8rH+7qcxcCQrv2NlQY8tjgx7AxttHdw5nP3Rp8K8F2wbrjE0oSaNW0cp33dYS7h+v6d3ZtNoUNhTq/caj0EBbuSmY814dEz604lflT9I0ppYlnn3A7km7p0EVRJKogiuFuw7EwstA5DK2oJaogil42vXCzcCPCMZTSE/Wse/ln8n5UYN/oykinP/CY30Sy0SGya7OJLpSc5yj3UUS4R9CqaSWhJIH9RfsREQl3D2e0x2hA6gwplSnUttYS7h7OMJdhGCoMiS2MpaK5guNVxwl3D8fP2g9Xc1f2F+1HrVUTVxRHuHs4FkYWDHUequvAUQVRDHQciI2JDaM9RrO/aD9tmjaiC6NxMnMi0DaQ0R6jSa9Op7SxlOiCaAwUBgx3G06ERwR1bXUkVyTrRm8RHhGMdh+NiMj+ov0klibSqmklwiOCcPdwBASiC6PJq8+joKGAcPdw+jv0x8rIipjCGJpUTRwqPcQot1E4mjkSZBdEbGGsNKrO30uYaxgWRhaEu4W3a29/a388LD0Y7TGaI2VHqGurI6ogCnNDcwY7DWa0x2iKG4vJrMls394eEbRp20goSeBA0QG0opYId+keAKILo0mtSKW6tZpw93CGugzFSGFEdGE0VS1VpFamMsp9FL5WvrhbuBNbFItaq+ZA0QFGuo3E3NCcoa5DdV/K7eokNj06GI+vvsTc1ILntuWy5sjXzK1JQNlqwbO/VBJkFcbJ6pOUNJYQXRiNgWDAcFepvetV9SSVJ+nuIdw9nAiPCAD2F+3nUMkhmtXNuvZ2qoHiRe9yatw4IqNrKRwWzGNTXmT9MAv6BJez+NYQTK1ymdtnAr3t3Qi2C9a1976CfQxzHYalkSXh7uHEFMbo2tvHygdPK08i3CNIKk+itrWWqIIozAzMCHUOJcIjgtImyQm20+oeIfXH4jjiiuLQiBqpz5zTv49XHaeqpYoI9wiGugzFRGlCTGEMNS01JFckE+4ejo+1D56Wnuwv2o9Gq9G1t5mhGcNchrE3f6/uO9rPsR/2pvZEuEcQXxJPi7qFqIIo7E3sCbEPIcIjgoyaDIoaioguiEYpKBnhNoIIjwgaVA0cLTuqGy1HuEvtekbr4bLDNKubCXcPZ5T7KBSCgpjCGArqC8ipyyHcPZw+Dn2wNbYltiiWFnULiSWJjHIfhZ2JHX3s++ie7PcV7GOYy9n2ji6MRqPVsK9gH95W3nhbeTPaYzRJ5UnUtNQQVRCFqYEpoS5Se5c3l5Nena5r7wj3CKm9RbU0vVlxkESPbQx82ArPYDuMTNq/tO4oLsuhC4IwWRCEdEEQMgRBeP4C+wVBEBaf3n9MEIT/Hj3/D6RXp1PUWESkRyQAkZ6RlDeXE1ccx+HSw4w2H0/i0tVYr/JlyvGFlFabc8zzT0rDvmbQ/15ndMBMQHKo+wr20ce+Dw6mDgxxGYKpgSnRBdFEFUZha2xLH/s+eFt542npSXRBNHvz96IUlLoOHOocSkxhjM5pRHhEIAgC4e7hxBXHcbDkIPWqel1HjPSMJLs2m/jieI5XHdd9mUZ7jKZJ3cSBogPsL9pPhLt0njP7owuj2VewjyHOQzA3NGeE2wiUgpKogiiiCqPwsPDAx8qHYPtg7EzsdFpNDUwJdQ7VdeCYwhjdKCPCPQKlQrqX2KJYDhQdQKVVEe4haR3lNoqjZUdJKk+isKFQ194RHhFUtlQSVxxHYmliu3vQiBpiC2OJLohmpNtIDJWGunuPKowiqiCKYLtgnMycGOw0GDMDM6ILo4kqiMLa2Jp+Dv3wtPLEx8qH6MJo9uTvQSEoGOU2ClMDU4a6DCWmMIZ9+fvatfcot1EcLD5IYmkidW11unuI9Igkty6X+OJ4UitTGe0ZieW4cbStfJcvpykwFJq5df0pvlnezI37fyf612oQRT6N/ZPNGbsJsRuAhZEFw12HYyAYEF0QTXRBNG7mbvhZ+xFkF4SDqQPRBZJWO5UxwYcrqX/4WRYvUeO29Qgn+njx1ANKHvKYhcbWgQiPUVRqkjG3yUKlbSPitNZw93CSypNIKk8ivz6fMR5jdPdQ1VLFgaIDJJQk6P4fzrT3mf433HV4u/aOLoxmX/4+Am0DcTZ3ZpDTIMwMzHTHWxlZ0c+hHx6WHvha++r6hkJQMMp9FMZKY117RxVGISIS4R6h6xsHS6T2rmmtOdu/PSIpaCggviSe5Ipk3Y9zpGckzepm9hftJ7YolnD3cBSCQveDeOa7ONBpIFZGVlJ7Kwx0fcPF3IVeNr0ItA3E0dRR17+NlcYMcxmGtbE1/R36t+/fHhEoBAUj3Ueyv3A/ccVxtGhazt6D+yiSypM4WnaUnLocXT+O9IikprWGA8UHSChOaNe/taJWpynMNQwjpdHZ/n36HnrZ9MLVwpWBTgOxMLTQHW9paElY0CCmLOxHwDCXK/B4l88lHbogCErgS2AKEALcIghCyN8OmwL0Pv3nAeBrOpHtOdtRCkrGeY6lrTSXPrltBJSHsmbpXmYfeRajP0KIP+KIjYmW2KDVrBn9Ffs9djMh4imwcMLVwpVeNr34Me1HjpUfY4L3BEAK3QpzCWNrzlZ25e5irNdYXfhXhHsEB0sOsiFjA6EuoVgbWwPSlzCzNpMVyStwt3An0DZQZ29UNfJq7KuYKE0Y5TYKQNc5Xox+EYBJ3pMACHMNw1hpzBtxb9CoamSSj2TvZdMLV3NXliQtIas2S3e8lZEVA50Gsi5jHfsL9zPBewKCIKAQFNIIoyCav3L+YpTbKF04V7hHOMkVyaw6sYoguyA8LD10WiuaK1iUsAgbYxuGOA+R7tlDGmE8ve9pFIKCcV7jdMcrBAWvxL6CWqvmOp/rAOjv0B9rY2s+TPyQsuYynVZnc2cCbQP5Oe1njpQd0bW3kdKI4a7D2Zq9lZ15Oxnreba9w93DSShO4M/MPxniPARbE1udPacuh+XJy3EzdyPELkSntUndxMsxL2OsND7rSE47P117n27XYe4jiB9kxlP3Crw1T4Fh/xAmZ8by/vbvWfK5iO/XX9H3QBZ1u02Z8/42/kiooI/9ANaeWktMYQxjPcfR0txKQ3Y+s0t9sf9pByGvr+bLj5uo+N+L5B9NY/VAfx5ZYMDLY2optHBj3qBBbHw0nKn+Y6lsqWTRwUVYGVkx1GWo7h40ooZn9j2DQlAw3ns8IDkepaDk1dhXUWlVunvo59APW2NbPj70MaVNpTq7k5kTwXbBrDq+isNlh3Xtbag0ZITbCLbnbOev3L8Y4zlGF24X4R5BQkkC606tY5DTIOxM7HTtnVuXy9KkpbiYu9DHvo9Oa7O6mZdiXsJIYXTWGXpK7f1S9Evt2nuoy1BMDUx5O+5t6tvqdX3G18oXDwsPvkn+hpPVJ3V9xtzQnCHOQ/gz40+iCqOY4CX17zODpZjCGLZmb2WE6whdmGG4eziplan8fPxnetn0wtvKW3dv1a3VvB3/NpaGlgxzGaa7B62o5el9TyMg6NpppPtIlIKS12Jfo03bptPU16EvdiZ2fHr4U4obi3V2B1MH+tj34dcTv5JYkshE74lSeyuk9v4r9y+25Wwj0jPyvFDMjuZyRujDgAxRFLNEUWwDfgVm/e2YWcAPokQcYCMIQqfkl9285lsafnDkroRXWPPUAZa/lsm23+wYlzEf94oBaIyqCfeN4s4HNMx+/zYChrlS0lKMpZGl7ssNMM1vGmXNUkTGdL+zRZ6n+U+jtrWWFk2Lbi4cYKb/TFo1rZQ1lzHNd5rOPsV3CkpBSV59HtP9puviVCM8InAzd6O0qZTx3uN1nc7T0pM+9n0oay5joONAPK08ATA1MCXSI5KypjKczZwJcwkDpBC8Kb5TKG0qxVBhqPsinLmHiuYK1KKa63tdr7NP95tOvaqe6tZqpvufvbcZfjPQiloKGwrb3fMkn0nYmdhJoY8+k3WdbrDTYAJsAyhrKmOk20gcTB0AsDWxZbjrcMqayvC19iXEXnKqSoWSSd6TKGsqw9LQUvflBpjqN1XX3jP9Z7a7h7q2OprVzczuPfusVv8ZtGnbKGsqY4bfDJ19su9kDAQD8uvzmeY3Tdfe4e7heFh4UNpUyjivcbr2drNwo59DP8qay+jv2F/3JTcxMGGs51hKm8so7uNMn+U/EBAbg+vbb9PSz5f+uW0s2Kbl/S17efPbJxl0/0wWvnKM1z8r59MvWxi94DeyBg+mcOpkxn95gBkxrZg0qVgX3IenIx7mw7vew/mhV6i0AYVRFU+NuJl3buiHpYkhE70nXrC9BzoOJNgumNKmUka4jtC1t7WxNcPdhlPWXIaPlQ/9HPqdbW8fqb3NDc0Z6zm2Xbueae8Z/mfbb4bfDF17n9u/Z/jPkOarm0rb9Y3rfK7DUGF4Xv8e6TYST0tPXXubG0oL8VzMXRjgOICy5jL62vfF19oXkAZLYz3HUtZchoOpA8PdhgNS/57qN5XSplKUgpIpvlPO3oPvNKpaqlBr2/fvGf4zaFI3UdVS1a5/T/OTvpcFDQW6fwNM8J6Ag6mD1N6+k3Wx+P0d+hNiH0JZUxnDXYfjZOYESIOlkW4jKWsuw8vSiwGO0sIuhaDgOh8ppNbMwIzxXuN115jqO5Xy5nJExHb9e7rfdOrb6s9r787ikguLBEGYA0wWRfG+09vzgTBRFB8555hNwCJRFGNOb+8CnhNFMfFv53oAaQSPl5fXkNzcXK6UpL1b2LzjFAOMTbE0VmJuZ4W1lxttHjZsrN/E9QHXE2gXqDu+trWW5ceWE+kZqRsJgRT2uOzYMoLsgnS/qCDNz69MXYmFoQVzA+e2u/baU2upaqni7j53t1u4sSdvD8kVydzX7752ixKOlh1lZ+5O7up7l+7LCZBVk8XvJ39nXuA8fKx9dPbypnK+S/2OST6TdJ0IpAiGpceWMsR5iG4kBFIUxjfJ3+Bm4dauEwH8cvwX1Fo180Pmt1sMsTV7K9m12dzX7752CzTii+OJLYrl3r736p4+QIpU2JC5gduDb8fNwk1nL2wo5Ke0n7i+1/nt/U3yN0S4RzDMdZjO3qZpY+mxpQTYBrT7URJFke9Sv8PMwIx5QfPa3cO6U+uoaK7gnr73tGvvvfl7SSpP4v5+97dr76TyJHbk7Di/vWuz+D39d+YGztU5GJCiiVamrGSi90QGOg1s197LkpYSqvViQKMdbVlZlOUVU1ZYRkXdKRSGpliaetFmbUebnSMtLm7E22TRaiwy0uFGhnjbEexqBcC2nG1k1mRyf7/727V3QkkC0QXR3NuvfXunV6XzZ+af3BZ8G+4W7jp7UUMRPx3/iZn+MwmyC2rX3iuSVzDSfSTDXYe3a+/lycvxt/Fnss/kdu39Q9oPGCuNmRc4r13f+DPjT0qbSrmn7z3tFspEFURxpOzIee2dXJ7M9pzt3NnnThzNHHX27Npsfkv/jZsCbsLP5mzO9n9q7+XJyxngOIAxnmN0do1Ww4qUFTiaOnJD7xs4l9/Sf/v/9u4txKoqjuP490cmll3MbthVg8hITIOsKQmzKyL1EJFYaBlBZFFURBZEEdFDEfVSUWbRhXroYnYhDUt66Gpp45g6FQ1kWVNgmg2U5b+HvUbPnDRnRo/7rD2/DxzO3uvsOazfzJz/7Fnn7LXo2tzFzJNm9siwqGMR7evbuWbsNT0ulur+fs8aM4thQ4ZtbW9f3878b+YzffT0rf+xAqzbtI7nVj3H1OOmbj1hAdj410bmts6l5YgWWo5o2dq++Z/NPLniSUYdOKrHH6Xu7/fgvQYz7YRpu+VK5V26UlTSpcAFdQV9QkTcUHPMW8D9dQX9toj4fEfP20yX/puZ5WJXrxRdCxxds38U8GM/jjEzswbqTUH/DDhe0ihJg4FpwIK6YxYAM9KnXU4HNkTEuvonMjOzxtnpEnQR8bek64GFwF7AvIhYKena9PjjwNvAFOAboAu4qnFdNjOz7enVmqIR8TZF0a5te7xmO4DZu7drZmbWF1leKWpmZv/lgm5mVhEu6GZmFeGCbmZWEaUtQSfpF6Dvl4oWDgF+3Y3dKZOzNKeqZKlKDnCWbsdGxKHbe6C0gr4rJC3d0ZVSuXGW5lSVLFXJAc7SGx5yMTOrCBd0M7OKyLWgP7HzQ7LhLM2pKlmqkgOcZaeyHEM3M7P/yvUM3czM6rigm5lVRHYFfWcLVjczSfMkdUpqq2kbLuldSV+n+4PK7GNvSDpa0vuSVklaKenG1J5jliGSPpX0ZcpyT2rPLgsUawBLWpZWEcs5R4ekFZKWS1qa2nLNMkzSy5JWp9dMS6OyZFXQe7lgdTN7Briwru12YHFEHA8sTvvN7m/glog4ETgdmJ1+Djlm+ROYHBEnA+OAC9Oc/jlmAbgRWFWzn2sOgLMjYlzN57VzzfII8E5EjAZOpvj5NCZLRGRzA1qAhTX7c4A5ZferjxlGAm01+2uAEWl7BLCm7D72I9PrwHm5ZwH2Bb4ATssxC8VKYYuBycCbqS27HKmvHcAhdW3ZZQEOAL4jfQCl0VmyOkMHjgS+r9lfm9pydnik1Z3S/WEl96dPJI0ExgOfkGmWNEyxHOgE3o2IXLM8DNwGbKlpyzEHQACLJH2eFpeHPLMcB/wCPJ2GwuZKGkqDsuRW0Le3ZLY/d1kSSfsBrwA3RcTGsvvTXxHxT0SMozjDnSBpTMld6jNJU4HO+J+F2TNzZkScQjG8OlvSWWV3qJ8GAacAj0XEeOAPGjhUlFtBr+Ji1D9LGgGQ7jtL7k+vSNqbopi/EBGvpuYss3SLiN+AJRTvc+SW5UzgIkkdwEvAZEnPk18OACLix3TfCbwGTCDPLGuBtem/PoCXKQp8Q7LkVtB7s2B1bhYAM9P2TIrx6KYmScBTwKqIeKjmoRyzHCppWNreBzgXWE1mWSJiTkQcFREjKV4X70XEFWSWA0DSUEn7d28D5wNtZJglIn4Cvpd0Qmo6B/iKRmUp+02DfrzJMAVoB74F7iy7P33s+4vAOmAzxV/uq4GDKd7I+jrdDy+7n73IMZFiqKsVWJ5uUzLNMhZYlrK0AXel9uyy1GSaxLY3RbPLQTHu/GW6rex+neeYJfV7HLA0/Y7NBw5qVBZf+m9mVhG5DbmYmdkOuKCbmVWEC7qZWUW4oJuZVYQLuplZRbigW2VIOjjNzrdc0k+SfkjbmyQ9Wnb/zBrNH1u0SpJ0N7ApIh4suy9me4rP0K3yJE2qmR98aJqX/rM0WdLFqf1KSfMlvSHpO0nXS7o5HfOxpOHpuCWSHpb0oaQ2SRNS+/D09a3p+LHlJbaBygXdBpo7KS6LPxU4G3ggXV4OMAaYTjFvyH1AVxQTKn0EzKh5jqERcQZwHTAvtd0DLIuIscAdwLMNT2JWZ1DZHTDbw86nmMTq1rQ/BDgmbb8fEb8Dv0vaALyR2ldQTBHQ7UWAiPhA0gFpLpiJwCWp/b00nn9gRGxobByzbVzQbaARcElErOnRKJ1GsXpRty01+1vo+Vqpf+Mp8NTO1gQ85GIDzULghjRjJJLG9+M5LktfOxHYkM7CPwAuT+2TgF8j4zniLU8+Q7eB5l6KlX1aU1HvAKb28TnWS/qQYnmxWantbopVaVqBLrZNjWq2x/hji2Z9IGkJcGtELC27L2b1PORiZlYRPkM3M6sIn6GbmVWEC7qZWUW4oJuZVYQLuplZRbigm5lVxL8mLFrdw7Go1gAAAABJRU5ErkJggg==\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"Ω1 = 1\n",
"Ω2 = 1\n",
"Δ1 = 5\n",
"Δ2 = 5\n",
"\n",
"## Dinámica completa\n",
"H = Ω1/2*(k3*k1.dag()+k1*k3.dag())+\\\n",
" Ω2/2*(k3*k2.dag()+k2*k3.dag())+\\\n",
" Δ1*k1*k1.dag()+\\\n",
" Δ2*k2*k2.dag()\n",
"\n",
"rho0 = k1*k1.dag()\n",
"rslt = qu.mesolve(H,rho0,t,[],[k1*k1.dag(),k2*k2.dag(),k3*k3.dag()])\n",
"\n",
"\n",
"pl.plot(t,rslt.expect[0],label=r\"$\\rho_{11}$\")\n",
"pl.plot(t,rslt.expect[1],label=r\"$\\rho_{22}$\")\n",
"pl.plot(t,rslt.expect[2],label=r\"$\\rho_{33}$\")\n",
"\n",
"\n",
"## Eliminación adiabática\n",
"Δ = (Δ1+Δ2)/2\n",
"\n",
"H = Ω1*Ω2/(4*Δ)*(k2*k1.dag()+k1*k2.dag())+\\\n",
" (Δ1-Δ+Ω1**2/(4*Δ))*k1*k1.dag()+\\\n",
" (Δ2-Δ+Ω2**2/(4*Δ))*k2*k2.dag()\n",
"\n",
"rho0 = k1*k1.dag()\n",
"rslt = qu.mesolve(H,rho0,t,[],[k1*k1.dag(),k2*k2.dag()])\n",
"\n",
"\n",
"pl.plot(t,rslt.expect[0],label=r\"$\\rho_{11eff}$\")\n",
"pl.plot(t,rslt.expect[1],label=r\"$\\rho_{22eff}$\")\n",
"pl.xlabel(\"Tiempo\")\n",
"\n",
"pl.legend(loc=0)"
]
},
{
"cell_type": "code",
"execution_count": null,
"id": "b7e8bdfe",
"metadata": {},
"outputs": [],
"source": []
}
],
"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.9.5"
}
},
"nbformat": 4,
"nbformat_minor": 5
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment