Files
bio-plots/bionet/3-transcription-regulation.ipynb
2024-01-05 11:55:33 +01:00

170 lines
70 KiB
Plaintext

{
"cells": [
{
"cell_type": "code",
"execution_count": 20,
"metadata": {},
"outputs": [],
"source": [
"import matplotlib.pyplot as plt\n",
"import numpy as np\n",
"from math import sqrt\n",
"\n",
"def pq (a,b,c):\n",
" return (-b + sqrt(b**2 - 4*a*c)) / 2*a"
]
},
{
"cell_type": "code",
"execution_count": 51,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"-3.792751279232142\n"
]
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAhYAAAGzCAYAAABzfl4TAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/OQEPoAAAACXBIWXMAAA9hAAAPYQGoP6dpAACKKElEQVR4nOzdd3gUVRfA4d/sJtn0SgmBQEILvfeOgoCKoChVIYigFAEpIipdQOlNAfWTItKbKIgUKdI70iSUFHooSUgvu/P9EbKyJIEENtmU8z7PPtmdnblzZibJnr33zr2KqqoqQgghhBBmoLF0AEIIIYTIOySxEEIIIYTZSGIhhBBCCLORxEIIIYQQZiOJhRBCCCHMRhILIYQQQpiNJBZCCCGEMBtJLIQQQghhNpJYCCGEEMJsJLEQ+c7YsWNRFMXSYeQZKefz3r17Wb6vqVOnUrJkSbRaLdWqVcvy/QkhMk8SC5GrLF68GEVRjA8rKyuKFi2Kv78/N27csHR4Igtt27aNTz/9lIYNG7Jo0SImTZpk9n0sX76cWbNmmb3cp9m2bRu9evWiUqVKaLVafHx80l3XYDAwZcoUfH19sbW1pUqVKqxYsSLNdS9cuEDr1q1xdHTE3d2d9957j7t372bRUQjxHytLByDE8xg/fjy+vr7ExcVx6NAhFi9ezL59+zh79iy2trZP3fbLL7/ks88+y6ZIhbn89ddfaDQa/ve//2FjY5Ml+1i+fDlnz55l8ODBWVJ+evtctWoVNWrUwMvL66nrfvHFF3z99df07t2b2rVr8+uvv9K1a1cURaFz587G9a5fv06TJk1wcXFh0qRJREVFMW3aNM6cOcORI0ey7PwJAYAqRC6yaNEiFVCPHj1qsnzEiBEqoK5atcpCkT1dVFSUpUPIMmPGjFEB9e7du1m6n549e6oODg5Zuo/XXntNLVGihFnLTExMVOPj49N9/8aNG2pCQsIz93/9+nXV2tpa7d+/v3GZwWBQGzdurBYrVkxNSkoyLu/bt69qZ2enBgcHG5dt375dBdSFCxe+4BEJ8XTSFCLyhMaNGwNw5cqVZ66bVh8LRVEYMGAAGzdupFKlSuh0OipWrMjWrVtTbX/jxg169eqFl5cXOp0OX19f+vbtS0JCAvBfc82ePXvo168fhQoVolixYsbt//jjDxo3boyDgwNOTk689tprnDt3zmQf//zzD/7+/pQsWRJbW1s8PT15//33uX//vsl6kZGRDB48GB8fH3Q6HYUKFaJly5acOHHCZL3Dhw/TunVrXFxcsLe3p2nTpuzfv/+5ykrPvXv36NixI87Oznh4eDBo0CDi4uKM7zdt2pSqVaumua2fnx+tWrVKt2xFUVi0aBHR0dHGZrDFixcDsGjRIl566SUKFSqETqejQoUKzJ8/P81y/vjjD5o2bYqTkxPOzs7Url2b5cuXA9CsWTM2b95McHCwcR+PN0uEhobSq1cvChcujK2tLVWrVmXJkiUm5QcFBaEoCtOmTWPWrFmUKlUKnU7H+fPn0z02Ly8vrK2t030/xa+//kpiYiL9+vUzOS99+/bl+vXrHDx40Lh83bp1vP766xQvXty4rEWLFpQtW5bVq1c/c19CvAhpChF5QlBQEABubm7PXca+fftYv349/fr1w8nJiTlz5tChQwdCQkLw8PAA4ObNm9SpU4fw8HD69OlDuXLluHHjBmvXriUmJsakirlfv34ULFiQ0aNHEx0dDcDPP/9Mjx49aNWqFd988w0xMTHMnz+fRo0acfLkSeMH2fbt27l69So9e/bE09OTc+fO8f3333Pu3DkOHTpkTIw++ugj1q5dy4ABA6hQoQL3799n3759XLhwgRo1agDJTQht2rShZs2ajBkzBo1GY/ww/vvvv6lTp06Gy3qajh074uPjw+TJkzl06BBz5swhLCyMpUuXAvDee+/Ru3dvzp49S6VKlYzbHT16lICAAL788st0y/7555/5/vvvOXLkCD/++CMADRo0AGD+/PlUrFiRN954AysrK3777Tf69euHwWCgf//+xjIWL17M+++/T8WKFRk5ciSurq6cPHmSrVu30rVrV7744gsiIiK4fv06M2fOBMDR0RGA2NhYmjVrxuXLlxkwYAC+vr6sWbMGf39/wsPDGTRokEm8ixYtIi4ujj59+qDT6XB3d3/m+XuWkydP4uDgQPny5U2Wp1y/kydP0qhRI27cuEFoaCi1atVKVUadOnXYsmXLC8cixFNZuspEiMxIaQrZsWOHevfuXfXatWvq2rVr1YIFC6o6nU69du3aM8tIqbp/HKDa2Nioly9fNi47ffq0Cqhz5841Luvevbuq0WhSNcWoanK19OMxNmrUyKR6OjIyUnV1dVV79+5tst3t27dVFxcXk+UxMTGpyl+xYoUKqHv37jUuc3FxMakaTyumMmXKqK1atTLGl1K+r6+v2rJlywyXlZ6U8/nGG2+YLO/Xr58KqKdPn1ZVVVXDw8NVW1tbdcSIESbrDRw4UHVwcHhmc1GPHj3SbApJ61y1atVKLVmypPF1eHi46uTkpNatW1eNjY01Wffx85JeU8SsWbNUQF22bJlxWUJCglq/fn3V0dFRffjwoaqqqhoYGKgCqrOzsxoaGvrU40nL05pCXnvtNZNjShEdHa0C6meffaaqqqoePXpUBdSlS5emWnf48OEqoMbFxWU6NiEySppCRK7UokULChYsiLe3N2+//TYODg5s2rTJpMnhecosVaqU8XWVKlVwdnbm6tWrQHKP/I0bN9K2bds0vw0+2bzSu3dvtFqt8fX27dsJDw+nS5cu3Lt3z/jQarXUrVuXXbt2Gde1s7MzPo+Li+PevXvUq1cPwKRpwtXVlcOHD3Pz5s00j+nUqVNcunSJrl27cv/+feM+o6Ojefnll9m7dy8GgyFDZT3L47UDAB9//DGA8Ruyi4sL7dq1Y8WKFaiqCoBer2fVqlW0b98eBweH59rv4+cqIiKCe/fu0bRpU65evUpERASQfO4jIyP57LPPUnXuzcitx1u2bMHT05MuXboYl1lbWzNw4ECioqLYs2ePyfodOnSgYMGCz3U86YmNjUWn06VannI8sbGxJj8zsq4QWUGaQkSu9O2331K2bFkiIiL46aef2Lt3b5r/SDPj8fboFG5uboSFhQFw9+5dHj58aFKN/zS+vr4mry9dugTASy+9lOb6zs7OxucPHjxg3LhxrFy5ktDQUJP1Uj4sAaZMmUKPHj3w9vamZs2avPrqq3Tv3p2SJUua7LNHjx7pxhkREYGbm9szy3qWMmXKmLwuVaoUGo3G2EwF0L17d1atWsXff/9NkyZN2LFjB3fu3OG9997L0D7Ssn//fsaMGcPBgweJiYlJdWwuLi7GvjcZvXZPCg4OpkyZMmg0pt/FUpolgoODTZY/ee3Nwc7Ojvj4+FTLU/qxpCRYKT8zsq4QWUESC5Er1alTx1hr0L59exo1akTXrl25ePGisV08sx6vXXhcyrfrzHryn3dKzcDPP/+Mp6dnqvWtrP77c+zYsSMHDhxg+PDhVKtWDUdHRwwGA61btzaWk7Je48aN2bBhA9u2bWPq1Kl88803rF+/njZt2hjXnTp1aroDSqWcr2eVlVlp1QS0atWKwoULs2zZMpo0acKyZcvw9PSkRYsWmS4fkjvrvvzyy5QrV44ZM2bg7e2NjY0NW7ZsYebMmSbnKjtlxQd3kSJF2LVrF6qqmpzbW7duARhvVS1SpIjJ8sfdunULd3f3F07ChXgaSSxErqfVapk8eTLNmzdn3rx5WTZGRcGCBXF2dubs2bPPtX1KM0uhQoWe+kEaFhbGzp07GTduHKNHjzYuT6l9eFKRIkXo168f/fr1IzQ0lBo1ajBx4kTatGlj3Kezs3OGPryfVtazXLp0yeSb+uXLlzEYDCZ3Vmi1Wrp27crixYv55ptv2LhxY6omo8z47bffiI+PZ9OmTSY1To83K8F/5/7s2bOULl063fLSaxYpUaIE//zzDwaDwaTW4t9//zW+n9WqVavGjz/+yIULF6hQoYJx+eHDh43vAxQtWpSCBQty7NixVGUcOXJERiwVWU76WIg8oVmzZtSpU4dZs2aZ3OJoThqNhvbt2/Pbb7+l+U/7WTUbrVq1wtnZmUmTJpGYmJjq/ZRREVM+ZJ8s78kRIfV6vUmzCCQnLV5eXsZq8Jo1a1KqVCmmTZtGVFRUuvvMSFnP8u2335q8njt3LkCqpOS9994jLCyMDz/8kKioKN59990MlZ+WtM5VREQEixYtMlnvlVdewcnJicmTJ6f6/Xh8WwcHh1TnAeDVV1/l9u3brFq1yrgsKSmJuXPn4ujoSNOmTZ/7GDKqXbt2WFtb89133xmXqarKggULKFq0qPEuGUju4/H7779z7do147KdO3cSEBDAO++8k+WxivxNaixEnjF8+HDeeecdFi9ezEcffZQl+5g0aRLbtm2jadOm9OnTh/Lly3Pr1i3WrFnDvn37cHV1TXdbZ2dn5s+fz3vvvUeNGjXo3LkzBQsWJCQkhM2bN9OwYUPmzZuHs7MzTZo0YcqUKSQmJlK0aFG2bdtGYGCgSXmRkZEUK1aMt99+m6pVq+Lo6MiOHTs4evQo06dPB5KToR9//JE2bdpQsWJFevbsSdGiRblx4wa7du3C2dmZ3377LUNlPUtgYCBvvPEGrVu35uDBgyxbtoyuXbumGruievXqVKpUiTVr1lC+fPkM3cqanldeeQUbGxvatm1rTFR++OEHChUqZNIU4OzszMyZM/nggw+oXbs2Xbt2xc3NjdOnTxMTE2Mcj6JmzZqsWrWKIUOGULt2bRwdHWnbti19+vRh4cKF+Pv7c/z4cXx8fFi7di379+9n1qxZODk5Pfcx/PPPP2zatAlIruWJiIjgq6++AqBq1aq0bdsWgGLFijF48GCmTp1KYmIitWvXZuPGjfz999/88ssvJrU+n3/+OWvWrKF58+YMGjSIqKgopk6dSuXKlenZs+dzxypEhljuhhQhMi+9kTdVVVX1er1aqlQptVSpUia3eT4pvdtN07rVskSJEmqPHj1MlgUHB6vdu3c33uJasmRJtX///sbRFZ8Wo6qq6q5du9RWrVqpLi4uqq2trVqqVCnV399fPXbsmHGd69evq2+++abq6uqquri4qO+884568+ZNFVDHjBmjqqqqxsfHq8OHD1erVq2qOjk5qQ4ODmrVqlXV7777LtU+T548qb711luqh4eHqtPp1BIlSqgdO3ZUd+7cmemy0juf58+fV99++23VyclJdXNzUwcMGJDq1s4UU6ZMUQF10qRJzyw/RXq3m27atEmtUqWKamtrq/r4+KjffPON+tNPP6mAGhgYmGrdBg0aqHZ2dqqzs7Nap04ddcWKFcb3o6Ki1K5du6qurq4qYHLr5507d9SePXuqBQoUUG1sbNTKlSurixYtMik/5XbTqVOnZvi4Un5f0no8+bun1+vVSZMmqSVKlFBtbGzUihUrmtwC+7izZ8+qr7zyimpvb6+6urqq3bp1U2/fvp3huIR4XoqqPmfPNCGEeE6zZ8/mk08+ISgoKM27cYQQuZckFkKIbKWqKlWrVsXDwyNVJ0shRO4nfSyEENkiOjqaTZs2sWvXLs6cOcOvv/5q6ZCEEFlAaiyEENkiKCgIX19fXF1d6devHxMnTrR0SEKILCCJhRBCCCHMRsaxEEIIIYTZSGIhhBBCCLPJ9s6bBoOBmzdv4uTklKFZBYUQQghheaqqEhkZiZeXV6oJ+R6X7YnFzZs38fb2zu7dCiGEEMIMrl27RrFixdJ9P9sTi5Shb69du2YyTbQQQgghcq6HDx/i7e39zCHssz2xSGn+cHZ2lsRCCCGEyGWe1Y1BOm8KIYQQwmwksRBCCCGE2UhiIYQQQgizkblChBA5gqqqJCUlodfrLR2KEPmSVqvFysrqhYeCkMRCCGFxCQkJ3Lp1i5iYGEuHIkS+Zm9vT5EiRbCxsXnuMiSxEEJYlMFgIDAwEK1Wi5eXFzY2NjJ4nhDZTFVVEhISuHv3LoGBgZQpU+apg2A9jSQWQgiLSkhIwGAw4O3tjb29vaXDESLfsrOzw9ramuDgYBISErC1tX2ucqTzphAiR3jeb0dCCPMxx9+h/CULIYQQwmwksRBCCCGE2UhiIYQQ+ZS/vz/t27d/4XIWL16Mq6vrC5eTWeaKPzfYvXs3iqIQHh5u6VCeSRILIYR4Ds2aNWPw4MGWDuOFzJ49m8WLF2dqGx8fH2bNmmWyrFOnTgQEBJgvMJFKgwYNuHXrFi4uLpYO5ZnyzF0hM7ZdJDI+ydJhCJGvWWs1vFOzGGUKP332w/xCVVX0ej1WVjnrX61er0dRFLN9SNnZ2WFnZ2eWsrJSynFnRUfhhISEFxr74VlsbGzw9PTMsvLNKWf9tr+AlUevERoZb+kwhMj31hy7xpZBjSni8vwfNKqqEpuY/SNw2llrMzSGhr+/P3v27GHPnj3Mnj0bgMDAQIKCgmjevDlbtmzhyy+/5MyZM2zbtg1vb2+GDBnCoUOHiI6Opnz58kyePJkWLVoYy/Tx8aFPnz5cvnyZNWvW4ObmxpdffkmfPn2A5A+uIUOGsG7dOsLCwihcuDAfffQRI0eOBCA8PJwRI0awceNGIiIiKF26NF9//TWvv/46ixcvZvDgwSxdupTPPvuMgIAALl++zNixYwkPD2fjxo1Aci1MpUqVAPj555+xtramb9++jB8/HkVRaNasGcHBwXzyySd88sknQPK1Sin/8Wr6+fPnM23aNK5du4avry9ffvkl7733nvF9RVH44Ycf2Lx5M3/++SdFixZl+vTpvPHGG0ByEtCnTx/++usvbt++TfHixenXrx+DBg3K8PVM77iLFCnCF198wYoVKwgPD6dSpUp88803NGvWzGS7xYsXM3z4cK5du0bTpk358ccf8fb2BmDs2LFs3LiRAQMGMHHiRIKDgzEYDISHhzNs2DB+/fVX4uPjqVWrFjNnzqRq1aoAnD59msGDB3Ps2DEURaFMmTIsXLiQWrVqERwczIABA9i3bx8JCQn4+PgwdepUXn31VXbv3k3z5s0JCwszNjutW7eO0aNHG4/p448/ZujQoRn+ncoqeSax8G/oQ7TUWAhhUTvOh3LxTiQfLz/Jyj71sNI+3zfD2EQ9FUb/aebonu38+FbY2zz73+Ls2bMJCAigUqVKjB8/HoCCBQsSFBQEwGeffca0adMoWbIkbm5uXLt2jVdffZWJEyei0+lYunQpbdu25eLFixQvXtxY7vTp05kwYQKff/45a9eupW/fvjRt2hQ/Pz/mzJnDpk2bWL16NcWLF+fatWtcu3YNSB5krE2bNkRGRrJs2TJKlSrF+fPn0Wq1xrJjYmL45ptv+PHHH/Hw8KBQoUJpHtuSJUvo1asXR44c4dixY/Tp04fixYvTu3dv1q9fT9WqVenTpw+9e/dO9/xs2LCBQYMGMWvWLFq0aMHvv/9Oz549KVasGM2bNzeuN27cOKZMmcLUqVOZO3cu3bp1Izg4GHd3dwwGA8WKFWPNmjV4eHhw4MAB+vTpQ5EiRejYseMzr9HTjnvAgAGcP3+elStX4uXlxYYNG2jdujVnzpyhTJkyxu0mTpzI0qVLsbGxoV+/fnTu3Jn9+/cby758+TLr1q1j/fr1xnP9zjvvYGdnxx9//IGLiwsLFy7k5ZdfJiAgAHd3d7p160b16tWZP38+Wq2WU6dOYW1tDUD//v1JSEhg7969ODg4cP78eRwdHdM8ruPHj9OxY0fGjh1Lp06dOHDgAP369cPDwwN/f3/jek/7ncoqeSax6NestKVDECLf61jLm9fn7ONYcBjTtwcwonU5S4eUJVxcXLCxscHe3j7N6unx48fTsmVL42t3d3fjN1aACRMmsGHDBjZt2sSAAQOMy1999VX69esHwIgRI5g5cya7du3Cz8+PkJAQypQpQ6NGjVAUhRIlShi327FjB0eOHOHChQuULVsWgJIlS5rElJiYyHfffWcSR1q8vb2ZOXMmiqLg5+fHmTNnmDlzJr1798bd3R2tVouTk9NTq+WnTZuGv7+/8VhSamumTZtmklj4+/vTpUsXACZNmsScOXM4cuQIrVu3xtramnHjxhnX9fX15eDBg6xevTpTicWTxx0SEsKiRYsICQnBy8sLgGHDhrF161YWLVrEpEmTjNvNmzePunXrAskJV/ny5Tly5Ah16tQBkmuRli5dSsGCBQHYt28fR44cITQ0FJ1OZzwXGzduZO3atfTp04eQkBCGDx9OuXLJfxspiUxKbB06dKBy5cpA6mv4uBkzZvDyyy8zatQoAMqWLcv58+eZOnWqSWLxtN+prJJnEgshhOWV8HDg6w5V6L/8BPN3X6GOrzvN/dL+Zvw0dtZazo9vlQURPnu/5lCrVi2T11FRUYwdO5bNmzdz69YtkpKSiI2NJSQkxGS9KlWqGJ8rioKnpyehoaFA8odwy5Yt8fPzo3Xr1rz++uu88sorAJw6dYpixYoZk4q02NjYmJSfnnr16pk0B9WvX5/p06ej1+tNakCe5sKFC6mq2xs2bGhsNkrxeDwODg44Ozsbjxfg22+/5aeffiIkJITY2FgSEhKoVq1ahmJI8eRxnzlzBr1en+pcxcfH4+HhYXxtZWVF7dq1ja/LlSuHq6srFy5cMCYWJUqUMCYVkNzMERUVZVIOQGxsLFeuXAGSk6wPPviAn3/+mRYtWvDOO+9QqlQpAAYOHEjfvn3Ztm0bLVq0oEOHDuleswsXLtCuXTuTZQ0bNmTWrFkm1+ppv1NZRRILIYRZvValCIcDS7D0YDBDVp16rv4WiqJkqEkip3JwcDB5PWzYMLZv3860adMoXbo0dnZ2vP322yQkJJisl1IlnkJRFAwGAwA1atQgMDCQP/74gx07dtCxY0datGjB2rVrM9Rx0s7OLsfNwfK04125ciXDhg1j+vTp1K9fHycnJ6ZOncrhw4cztY8njzsqKgqtVsvx48dTJUrpNTuk58nrHBUVRZEiRdi9e3eqdVP6RYwdO5auXbuyefNm/vjjD8aMGcPKlSt58803+eCDD2jVqhWbN29m27ZtTJ48menTp/Pxxx9nKq7HPe0cZxW53VQIYXZfvFaeSkWdCYtJ5OPlJ0nUZ+0/MkuwsbHJ8BTv+/fvx9/fnzfffJPKlSvj6elp7I+RGc7OznTq1IkffviBVatWsW7dOh48eECVKlW4fv26WW75fPKD+9ChQ5QpU8b4IZyR4y5fvrxJXwRIPgcVKlTIcBz79++nQYMG9OvXj+rVq1O6dGnjt/4XUb16dfR6PaGhoZQuXdrk8XjzTlJSEseOHTO+vnjxIuHh4ZQvXz7dsmvUqMHt27exsrJKVXaBAgWM65UtW5ZPPvmEbdu28dZbb7Fo0SLje97e3nz00UesX7+eoUOH8sMPP6S5r/TOcdmyZTNcs5RVJLEQQpidzkrLt11r4KSzSu5vsS3vjXHg4+PD4cOHCQoK4t69e0/9FlimTBnWr1/PqVOnOH36NF27ds30t8YZM2awYsUK/v33XwICAlizZg2enp64urrStGlTmjRpQocOHdi+fbuxZmPr1q2ZPq6QkBCGDBnCxYsXWbFiBXPnzjW5E8PHx4e9e/dy48YN7t27l2YZw4cPZ/HixcyfP59Lly4xY8YM1q9fz7BhwzIcR5kyZTh27Bh//vknAQEBjBo1iqNHj2b6eJ5UtmxZunXrRvfu3Vm/fj2BgYEcOXKEyZMns3nzZuN61tbWfPzxxxw+fJjjx4/j7+9PvXr1jM0gaWnRogX169enffv2bNu2jaCgIA4cOMAXX3zBsWPHiI2NZcCAAezevZvg4GD279/P0aNHjcnK4MGD+fPPPwkMDOTEiRPs2rUr3URm6NCh7Ny5kwkTJhAQEMCSJUuYN29eps5xVpHEQgiRJUp4OPDN28ntuwv2XGHXv1nbrpvdhg0bhlarpUKFChQsWDBVf4nHzZgxAzc3Nxo0aEDbtm1p1aoVNWrUyNT+nJycmDJlCrVq1aJ27doEBQWxZcsW45gM69ato3bt2nTp0oUKFSrw6aefZrhG5XHdu3cnNjaWOnXq0L9/fwYNGmTSX2L8+PEEBQVRqlQpk/4Fj2vfvj2zZ89m2rRpVKxYkYULF7Jo0SLj7ZwZ8eGHH/LWW2/RqVMn6taty/37942dEF/UokWL6N69O0OHDsXPz4/27dtz9OhRkzt07O3tGTFiBF27dqVhw4Y4OjqyatWqp5arKApbtmyhSZMm9OzZk7Jly9K5c2eCg4MpXLgwWq2W+/fv0717d8qWLUvHjh1p06aNsZOqXq+nf//+lC9fntatW1O2bFm+++67NPdVo0YNVq9ezcqVK6lUqRKjR49m/PjxJh03LUVRVVXNzh0+fPgQFxcXIiIicHZ2zs5dCyEsYMyvZ1lyMBg3e2s2D2yMl6tpf4C4uDgCAwPx9fV97mmahXk0a9aMatWqpRpZM79Ja1yO/OJpf48Z/fyWGgshRJb6/PH+FivyZn8LIcR/JLEQQmSpx/tbHA8OY9qfFy0dkhAiC0liIYTIciU8HJjyqL/Fwr1X2Xr2toUjEmnZvXt3vm8GgeQxQ/JjM4i5SGIhhMgWbSoX4YNGvgAMW3Oaq3ejLByRECIrSGIhhMg2I9qUo46PO1HxSfRddoKYBJnfR4i8RhILIUS2sdZqmNe1OgWddFy8E8nI9WfI5hvThBBZTBILIUS2KuRsy7dda6DVKPx66iYbT920dEhCCDOSxEIIke3q+Lozsk3y7I4Ldl8mIUluQRUir5DEQghhEb0a+fJalSIkGVTuRyfI+BZC5BGSWAghLEJRFL7pUAVvd3v0BpXbEbHS3+IF+fj4ZMntollVrrns3r0bRVFe+BZRc5WT30liIYSwGEedFePeqIhGgZgEPbcfxlk6pBzheT/Ijx49ajKvh0hfs2bNGDx4sMmyBg0acOvWLVxcXCwTVB4hiYUQwqJKeDjgZm8DwN3IeCJiEy0cUdZJSEjI0vILFiyIvb19lu4jq2T1uckIGxsbPD09URTF0qHkapJYCCEszs5Ga0wurj2IIS4hCRKis/+RiaaYZs2aMWDAAAYMGICLiwsFChRg1KhRJs05Pj4+TJgwge7du+Ps7GysTVi3bh0VK1ZEp9Ph4+PD9OnTTcoNDg7mk08+QVEUkw+5ffv20bhxY+zs7PD29mbgwIFER0eb7O/xmg5FUfjxxx958803sbe3p0yZMmzatOmpxxUaGkrbtm2xs7PD19eXX375JdU64eHhfPDBBxQsWBBnZ2deeuklTp8+bbLOV199RaFChXBycuKDDz7gs88+o1q1asb3/f39ad++PRMnTsTLyws/Pz8Afv75Z2rVqoWTkxOenp507dqV0FDTmXG3bNlC2bJlsbOzo3nz5gQFBZm8f//+fbp06ULRokWxt7encuXKrFixwmTfe/bsYfbs2cZzHBQUlGZTyNOuVco5nzRpEu+//z5OTk4UL16c77///qnnOK+zsnQAQggBUMBJR1KUnqj4JELu3KPs//yyP4jPb4KNQ4ZXX7JkCb169eLIkSMcO3aMPn36ULx4cXr37m1cZ9q0aYwePZoxY8YAcPz4cTp27MjYsWPp1KkTBw4coF+/fnh4eODv78/69eupWrUqffr0MSnnypUrtG7dmq+++oqffvqJu3fvGhObRYsWpRvjuHHjmDJlClOnTmXu3Ll069aN4OBg3N3d01zf39+fmzdvsmvXLqytrRk4cGCqD/Z33nkHOzs7/vjjD1xcXFi4cCEvv/wyAQEBuLu788svvzBx4kS+++47GjZsyMqVK5k+fTq+vr4m5ezcuRNnZ2e2b99uXJaYmMiECRPw8/MjNDSUIUOG4O/vz5YtWwC4du0ab731Fv3796dPnz4cO3aMoUOHmpQbFxdHzZo1GTFiBM7OzmzevJn33nuPUqVKUadOHWbPnk1AQACVKlVi/PjxQHJtz5MJyrOuVYrp06czYcIEPv/8c9auXUvfvn1p2rSpMVnKd9RsFhERoQJqREREdu9aCJEDxcbGqufPn1djY2PVxCS9euFmhPrP1RuqOsY5+x/xURmOu2nTpmr58uVVg8FgXDZixAi1fPnyxtclSpRQ27dvb7Jd165d1ZYtW5osGz58uFqhQgWT7WbOnGmyTq9evdQ+ffqYLPv7779VjUajxsbGprkdoH755ZfG11FRUSqg/vHHH2ke08WLF1VAPXLkiHHZhQsXVMBY7t9//606OzurcXFxJtuWKlVKXbhwoaqqqlq3bl21f//+Ju83bNhQrVq1qvF1jx491MKFC6vx8fFpxpLi6NGjKqBGRkaqqqqqI0eONDlXqpp83gE1LCws3XJee+01dejQocbXTZs2VQcNGmSyzq5du0zKyei1evfdd42vDQaDWqhQIXX+/PlPPa6c6vG/xydl9PNbaiyEEDmGlVZDCQ97roQaONPjAgWdbPB0tsu+AKwz1z+hXr16Jk0V9evXZ/r06ej1erRaLQC1atUy2ebChQu0a9fOZFnDhg2ZNWuWyXZPOn36NP/8849J04SqqhgMBgIDAylfvnya21WpUsX43MHBAWdn51Q1EI/HZmVlRc2aNY3LypUrh6urq0kcUVFReHh4mGwbGxvLlStXALh48SL9+vUzeb9OnTr89ddfJssqV66MjY2NybLjx48zduxYTp8+TVhYGAZD8m3IISEhVKhQgQsXLlC3bl2TberXr2/yWq/XM2nSJFavXs2NGzdISEggPj4+0/1PMnqtHj/HiqLg6emZ7jnODySxEELkKHY2VhRztyfkAYTGga29Na72Ns/eMIdycMh408rTREVF8eGHHzJw4MBU7xUvXjzd7aytrU1eK4pi/LB+3jiKFCnC7t27U733eAKSEU+em+joaFq1akWrVq345ZdfKFiwICEhIbRq1SpTnTunTp3K7NmzmTVrFpUrV8bBwYHBgwdnWQdRc5/j3E4SCyFEjuNqb0Nsop67kfFcD4tFZ6XFzibtb/KWdPjwYZPXhw4dokyZMunWOgCUL1+e/fv3myzbv38/ZcuWNW5nY2ODXq83WadGjRqcP3+e0qVLmyn61MqVK0dSUhLHjx+ndu3aQHLtw+OdGWvUqMHt27exsrLCx8cnzXL8/Pw4evQo3bt3Ny47evToM/f/77//cv/+fb7++mu8vb0BOHbsmMk65cuXT9UB9dChQyav9+/fT7t27Xj33XcBMBgMBAQEUKFCBeM6aZ3jJ2XkWonU5K4QIUSO5Olsi6POCoOqEnw/mqQcODJnSEgIQ4YM4eLFi6xYsYK5c+cyaNCgp24zdOhQdu7cyYQJEwgICGDJkiXMmzePYcOGGdfx8fFh79693Lhxg3v37gEwYsQIDhw4wIABAzh16hSXLl3i119/ZcCAAWY7Hj8/P1q3bs2HH37I4cOHOX78OB988AF2dv81R7Vo0YL69evTvn17tm3bRlBQEAcOHOCLL74wJgEff/wx//vf/1iyZAmXLl3iq6++4p9//nnmbZzFixfHxsaGuXPncvXqVTZt2sSECRNM1vnoo4+4dOkSw4cP5+LFiyxfvpzFixebrFOmTBm2b9/OgQMHuHDhAh9++CF37twxWcfHx4fDhw8TFBTEvXv30qxhyMi1EqlJYiGEyJEURaG4uz02VhoS9AZCHsTkuJE5u3fvTmxsLHXq1KF///4MGjTomQNU1ahRg9WrV7Ny5UoqVarE6NGjGT9+vMldBuPHjycoKIhSpUpRsGBBILkdf8+ePQQEBNC4cWOqV6/O6NGj8fLyMusxLVq0CC8vL5o2bcpbb71Fnz59KFSokPF9RVHYsmULTZo0oWfPnpQtW5bOnTsTHBxM4cKFAejWrRsjR45k2LBh1KhRg8DAQPz9/bG1tX3qvgsWLMjixYtZs2YNFSpU4Ouvv2batGkm6xQvXpx169axceNGqlatyoIFC5g0aZLJOl9++SU1atSgVatWNGvWDE9PT9q3b2+yzrBhw9BqtVSoUMHY5PKkjFwrkZqiZvNf6sOHD3FxcSEiIgJnZ+fs3LUQIgeKi4sjMDAQX1/fND944hL1XA6NwqCqFHDU4eWajZ05n6JZs2ZUq1YtRw91nZO0bNkST09Pfv75Z0uHIp7iaX+PGf38lj4WQogczdZai7e7HcH3Y7gXFY+dtRY3h9zbmTM/iImJYcGCBbRq1QqtVsuKFSvYsWOHyXgVIu+SxEIIkeO52NlQyMlAaGQc18NjsbHS4KCTf185VUpzycSJE4mLi8PPz49169bRokULS4cmsoH8ZQohcoXCzjriEvU8jEsk+H4MpQs5YGNluZ75ad1uKZLZ2dmxY8cOS4chLEQ6bwohcgVFUfB2t8fOWkuSwUDQ/Rj0hpzVmVMIIYmFECIX0WoUSng4YKXREJeo51oOvFNEiPxOEgshRK5iY5U87LeiKDyMS+T2wzhLhySEeIwkFkKIXMdBZ4W3W/Jtp3cj43kQnTVDNQshMi9TiYVer2fUqFH4+vpiZ2dHqVKlmDBhglRFCiGynau9DYWcku+zvxEeS3R8koUjEkJAJu8K+eabb5g/fz5LliyhYsWKHDt2jJ49e+Li4pLmxDhCCJGVCjvriE/SExGbM+4UEUJkssbiwIEDtGvXjtdeew0fHx/efvttXnnlFY4cOZJV8QkhRLoURaGYm9wp8rz8/f1TDXX9PBYvXpzpmU3NwVzx59T9Wuq8vqhMJRYNGjRg586dBAQEAHD69Gn27dtHmzZt0t0mPj6ehw8fmjyEEMJcjHeKaLP3TpFmzZoxePDgLN9PVpo9e3aqCbyexcfHJ9Uw5p06dTJ+LuQlQUFBKIrCqVOnTJY/z3l7lrx0XjPVFPLZZ5/x8OFDypUrh1arRa/XM3HiRLp165buNpMnT2bcuHEvHKgQQqTHxkqDj4c9V+9G8zAukZsRcXi52D5zNs2spqoqer0eK6ucNRahXq9HURRcXFzMUp6dnZ3JDKh5nbnO27Pk1vOaqRqL1atX88svv7B8+XJOnDjBkiVLmDZtGkuWLEl3m5EjRxIREWF8XLt27YWDFkLkbaqqEpMYk6kHSgIFnCBOH8uNiHCuhUdkuoyM1nT4+/uzZ88eZs+ejaIoKIpCUFAQu3fvRlEU/vjjD2rWrIlOp2Pfvn1cuXKFdu3aUbhwYRwdHaldu3aqkSl9fHyYNGkS77//Pk5OThQvXpzvv//e+H5CQgIDBgygSJEi2NraUqJECSZPnmx8Pzw8nA8//JDChQtja2tLpUqV+P3334H/qtQ3bdpEhQoV0Ol0hISEpKrSb9asGQMGDGDAgAG4uLhQoEABRo0aZTwvzZo1Izg4mE8++cR43I+X/7j58+dTqlQpbGxs8PPzSzX5mKIo/Pjjj7z55pvY29tTpkwZNm3aZHxfr9fTq1cv480Cfn5+zJ49O0PXJ8X9+/fp0qULRYsWxd7ensqVK7NixQqTdQwGA1OmTKF06dLodDqKFy/OxIkTAfD19QWgevXqKIpCs2bNANOmkO+//x4vL69U0663a9eO999/H+CZ1z87z2t2yFQaPXz4cD777DM6d+4MQOXKlQkODmby5Mn06NEjzW10Oh06ne7FIxVC5BuxSbHUXV432/d7uOth7K3tn7ne7NmzCQgIoFKlSowfPx5InvI7KCgISK7dnTZtGiVLlsTNzY1r167x6quvMnHiRHQ6HUuXLqVt27ZcvHiR4sWLG8udPn06EyZM4PPPP2ft2rX07duXpk2b4ufnx5w5c9i0aROrV6+mePHiXLt2zfhFzWAw0KZNGyIjI1m2bBmlSpXi/PnzaLX/dWSNiYnhm2++4ccff8TDw8NkKvTHLVmyhF69enHkyBGOHTtGnz59KF68OL1792b9+vVUrVqVPn360Lt373TPz4YNGxg0aBCzZs2iRYsW/P777/Ts2ZNixYrRvHlz43rjxo1jypQpTJ06lblz59KtWzeCg4Nxd3fHYDBQrFgx1qxZg4eHBwcOHKBPnz4UKVKEjh07PvMaQfJMnTVr1mTEiBE4OzuzefNm3nvvPUqVKkWdOnWA5C+/P/zwAzNnzqRRo0bcunWLf//9F4AjR45Qp04dduzYQcWKFbGxST353TvvvMPHH3/Mrl27ePnllwF48OABW7duZcuWLQBERUU99fpn53nNDplKLGJiYtBoTCs5tFptqkxNCCHyMhcXF2xsbLC3t8fT0zPV++PHj6dly5bG1+7u7lStWtX4esKECWzYsIFNmzYxYMAA4/JXX32Vfv36ATBixAhmzpzJrl278PPzIyQkhDJlytCoUSMURaFEiRLG7Xbs2MGRI0e4cOECZcuWBaBkyZImMSUmJvLdd9+ZxJEWb29vZs6ciaIo+Pn5cebMGWbOnEnv3r1xd3dHq9Xi5OSU5nGnmDZtGv7+/sZjGTJkCIcOHWLatGkmH4D+/v506dIFgEmTJjFnzhyOHDlC69atsba2NmlG9/X15eDBg6xevTrDiUXRokUZNmyY8fXHH3/Mn3/+yerVq6lTpw6RkZHMnj2befPmGb8clypVikaNGgHJySKAh4dHusfr5uZGmzZtWL58uTGxWLt2LQUKFDAea9WqVZ96/bPzvGaHTCUWbdu2ZeLEiRQvXpyKFSty8uRJZsyYYazuEUIIc7CzsuNw18PPvb2qqoQ8iCUyPhErjYZSBR2w1j675dfOyjzt2bVq1TJ5HRUVxdixY9m8eTO3bt0iKSmJ2NhYQkJCTNarUqWK8bmiKHh6ehIaGgokf1i0bNkSPz8/Wrduzeuvv84rr7wCwKlTpyhWrJgxqUiLjY2NSfnpqVevnknflPr16zN9+nT0er1JDcjTXLhwgT59+pgsa9iwYaqmjMfjcXBwwNnZ2Xi8AN9++y0//fQTISEhxMbGkpCQQLVq1TIUAyQ3p0yaNInVq1dz48YNEhISiI+Px97e3hhnfHy8MSF4Xt26daN3795899136HQ6fvnlFzp37mz8Ip7R6/8s5jqvWS1TicXcuXMZNWoU/fr1IzQ0FC8vLz788ENGjx6dVfEJIfIhRVEy1CTxNGUL2XP1bhSxiXpuh6uUKmiLVQaSC3NwcHAweT1s2DC2b9/OtGnTKF26NHZ2drz99tskJJiOGGptbW3yWlEUY41wjRo1CAwM5I8//mDHjh107NiRFi1asHbt2gx18LOzs7N4Z9YnPe14V65cybBhw5g+fTr169fHycmJqVOncvhwxhPOqVOnMnv2bGbNmkXlypVxcHBg8ODBxvNuro6Rbdu2RVVVNm/eTO3atfn777+ZOXOm8f2MXn9zedp5zQ6ZSiycnJyYNWtWqltihBAip9FqFHw8HLh8N4r4JD3BD2LwLeCAxkwfrjY2Nuj1+gytu3//fvz9/XnzzTeB5G+wKf0xMsPZ2ZlOnTrRqVMn3n77bVq3bs2DBw+oUqUK169fJyAg4Km1Fhnx5Af3oUOHKFOmjLG2IiPHXb58efbv32/S927//v1UqFAhw3Hs37+fBg0aGKv9IbkTZGbs37+fdu3a8e677wLJfVECAgKMcZQpUwY7Ozt27tzJBx98kGr7lD4VzzpeW1tb3nrrLX755RcuX76Mn58fNWrUMInjWdc/u85rdshZ90AJIYQZWVtp8PFw4OrdKKLjk7gRFksxN/N8c/fx8eHw4cMEBQXh6Oj41I5xZcqUYf369bRt2xZFURg1alSmv0HOmDGDIkWKUL16dTQaDWvWrMHT0xNXV1eaNm1KkyZN6NChAzNmzKB06dL8+++/KIqS6Xb1kJAQhgwZwocffsiJEyeYO3cu06dPNznuvXv30rlzZ3Q6HQUKFEhVxvDhw+nYsSPVq1enRYsW/Pbbb6xfvz7VnTBPU6ZMGZYuXcqff/6Jr68vP//8M0ePHjXeqZHRMtauXcuBAwdwc3NjxowZ3Llzx/hBbGtry4gRI/j000+xsbGhYcOG3L17l3PnztGrVy8KFSqEnZ0dW7dupVixYtja2qZ7q2m3bt14/fXXOXfunDGReTyOZ13/7Dqv2UEmIRNC5Gl2NlqKe9ijoBAWk0BoZLxZyh02bBharZYKFSpQsGDBp7aXz5gxAzc3Nxo0aEDbtm1p1aqVyTfajHBycmLKlCnUqlWL2rVrExQUxJYtW4zt+OvWraN27dp06dKFChUq8Omnn2a4RuVx3bt3JzY2ljp16tC/f38GDRpk0q4/fvx4goKCKFWqlLFz45Pat2/P7NmzmTZtGhUrVmThwoUsWrTIeLtmRnz44Ye89dZbdOrUibp163L//n2T2ouM+PLLL6lRowatWrWiWbNmeHp6phoxc9SoUQwdOpTRo0dTvnx5OnXqZOyPYGVlxZw5c1i4cCFeXl60a9cu3X299NJLuLu7c/HiRbp27WryXkauf3ad1+ygqNk8g9jDhw9xcXEhIiICZ2fn7Ny1ECIHiouLIzAwEF9fX2xtbbNsP/ej47kRFgtAUTc7PBzkNvgnNWvWjGrVqklzdz72tL/HjH5+S42FECJf8HDQGWdDvRkWy8PYRAtHJETeJImFECLfKOysw93eBhUIeRAjU60LkQWk86YQIt9QFIWibnYkGVQexiUSdD+aUgUdsbWWqdYBdu/ebekQRB4gNRZCiHxFURS83e2xt7FCb1AJuhdNYpKMHiyEuUhiIYTId5LHuLBHZ6UlQW8g8H40STI1gRBmIYmFECJfstJq8C1gj5VWQ1yinuD7MRgM2XqTnBB5kiQWQoh8y8ZKi6+HPVpFITo+iWthGZ86XQiRNkkshBD5mp2NFSU87FEUhYjYRG5GxElyIcQLkMRCCJHvOdpa4+2WPCHV/ah4s43OKUR+JImFEEIArvY2eLkmJxd3HsZxL8oyyYWiKGzcuNEi+36W3bt3oygK4eHhlg4l2/j4+MhIpJkkiYUQQjxSwFFHYedHo3OGxxIWnTXTWgOMHTuWatWqpVp+69Yt2rRpY7b95MdkwJyOHj1qMleKeDYZIEsIIR5TyEmH3qByLyqe62GxaDUKznbW2bZ/T0/PbNtXdkhMTMTaOvvOHyRPc64oinGCtheR3oRgKSxxfDmd1FgIIXIcVVUxxMRk+0NVVRRFoYiLLW72NqioBD+IISou9dDfW7dupVGjRri6uuLh4cHrr7/OlStXTNa5fv06Xbp0wd3dHQcHB2rVqsXhw4dZvHgx48aN4/Tp0yiKgqIoLF68GDBtCmnQoAEjRowwKfPu3btYW1uzd+9eAH7++Wdq1aqFk5MTnp6edO3a1Tg7Z1BQEM2bNwfAzc0NRVHw9/cHwGAwMHnyZHx9fbGzs6Nq1aqsXbvWZF9btmyhbNmy2NnZ0bx5c4KCgp557RRFYf78+bzxxhs4ODgwceJEAH799Vdq1KiBra0tJUuWZNy4cSQlJRmv99ixYylevDg6nQ4vLy8GDhxoLDM+Pp5hw4ZRtGhRHBwcqFu3rskooYsXL8bV1ZVNmzZRoUIFdDodP/74I7a2tqlqagYNGsRLL71kfL1u3ToqVqyITqfDx8fHZIp4SN0Ukt7xif9IjYUQIsdRY2O5WKNmtu/X78RxFPvkO0SKudmhfzT0d/D9aHwLOmBv89+/zOjoaIYMGUKVKlWIiopi9OjRvPnmm5w6dQqNRkNUVBRNmzalaNGibNq0CU9PT06cOIHBYKBTp06cPXuWrVu3smPHDgBcXFxSxdOtWzemTJnC119/jaIoAKxatQovLy8aN24MJH9jnjBhAn5+foSGhjJkyBD8/f3ZsmUL3t7erFu3jg4dOnDx4kWcnZ2xs0vuRzJ58mSWLVvGggULKFOmDHv37uXdd9+lYMGCNG3alGvXrvHWW2/Rv39/+vTpw7Fjxxg6dGiGzuPYsWP5+uuvmTVrFlZWVvz99990796dOXPm0LhxY65cuWJsXhgzZgzr1q1j5syZrFy5kooVK3L79m1Onz5tLG/AgAGcP3+elStX4uXlxYYNG2jdujVnzpyhTJkyAMTExPDNN9/w448/4uHhQbFixRg9ejTr1q2jV69eQHJNxqpVq4zJwPHjx+nYsSNjx46lU6dOHDhwgH79+uHh4WFMwDJyfOIJajaLiIhQATUiIiK7dy2EyIFiY2PV8+fPq7GxscZl+uho9bxfuWx/6KOjTWLT6w3qldBI9fS1MPXcjXA1NiEp3eO4e/euCqhnzpxRVVVVFy5cqDo5Oan3799Pc/0xY8aoVatWTbUcUDds2KCqqqqGhoaqVlZW6t69e43v169fXx0xYkS6cRw9elQF1MjISFVVVXXXrl0qoIaFhRnXiYuLU+3t7dUDBw6YbNurVy+1S5cuqqqq6siRI9UKFSqYvD9ixIhUZaUV/+DBg02Wvfzyy+qkSZNMlv38889qkSJFVFVV1enTp6tly5ZVExISUpUXHBysarVa9caNG6nKHDlypKqqqrpo0SIVUE+dOmWyzqBBg9SXXnrJ+PrPP/9UdTqdMf6uXbuqLVu2NNlm+PDhJsddokQJdebMmU89vrwkrb/HFBn9/JZUSwiR4yh2dvidOG6R/T5Oo1Eo4eFA4L0oYhL0BN6LplRBB2ystFy6dInRo0dz+PBh7t27h+HRkOAhISFUqlSJU6dOUb16ddzd3Z87noIFC/LKK6/wyy+/0LhxYwIDAzl48CALFy40rnP8+HHGjh3L6dOnCQsLM4mjQoUKaZZ7+fJlYmJiaNmypcnyhIQEqlevDsCFCxeoW7euyfv169fPUNy1atUyeX369Gn2799v0myg1+uJi4sjJiaGd955h1mzZlGyZElat27Nq6++Stu2bbGysuLMmTPo9XrKli1rUmZ8fDweHh7G1zY2NlSpUsVknW7dulGvXj1u3ryJl5cXv/zyC6+99hqurq7GY2zXrp3JNg0bNmTWrFno9Xq02rQnp3vy+IQpSSyEEDmOoigo9vaWDgNImVfEgSt3o4lP0hN4L4aSBR1o27YtJUqU4IcffsDLywuDwUClSpVISEi+k8TuiSTleXXr1o2BAwcyd+5cli9fTuXKlalcuTKQ3BzTqlUrWrVqxS+//ELBggUJCQmhVatWxjjSEhUVBcDmzZspWrSoyXs6ne6FY3ZwcEi1v3HjxvHWW2+lWtfW1hZvb28uXrzIjh072L59O/369WPq1Kns2bOHqKgotFotx48fT/VB7+joaHxuZ2dnbC5KUbt2bUqVKsXKlSvp27cvGzZsMPZlMefxCVOSWAghxDMkzyviwNW7UcQn6TkREMLFixf54YcfjH0d9u3bZ7JNlSpV+PHHH3nw4EGatRY2Njbo9fpn7rtdu3b06dOHrVu3snz5crp37258799//+X+/ft8/fXXeHt7A3Ds2LFU+wFM9pXSwTEkJISmTZumud/y5cuzadMmk2WHDh16ZrxpqVGjBhcvXqR06dLprmNnZ0fbtm1p27Yt/fv3p1y5cpw5c4bq1auj1+sJDQ01nuvM6NatG7/88gvFihVDo9Hw2muvGd8rX748+/fvN1l///79lC1bNt3aCvFscleIEEJkgI1VcnJhrdWgc3DG1d2dBQsXcvnyZf766y+GDBlisn6XLl3w9PSkffv27N+/n6tXr7Ju3ToOHjwIJN9tEBgYyKlTp7h37x7x8WkPyOXg4ED79u0ZNWoUFy5coEuXLsb3ihcvjo2NDXPnzuXq1ats2rSJCRMmmGxfokQJFEXh999/5+7du0RFReHk5MSwYcP45JNPWLJkCVeuXOHEiRPMnTuXJUuWAPDRRx9x6dIlhg8fzsWLF1m+fPlzf9sfPXo0S5cuZdy4cZw7d44LFy6wcuVKvvzySyD5ro7//e9/nD17lqtXr7Js2TLs7OwoUaIEZcuWpVu3bnTv3p3169cTGBjIkSNHmDx5Mps3b37mvrt168aJEyeYOHEib7/9tkmNzNChQ9m5cycTJkwgICCAJUuWMG/ePIYNG/ZcxykeyaoOIOmRzptCiMc9rbNYThSbkKSeuxGhLly+QS1V1k/V6XRqlSpV1N27d5t0vFRVVQ0KClI7dOigOjs7q/b29mqtWrXUw4cPq6qa3IGyQ4cOqqurqwqoixYtUlVVTVWGqqrqli1bVEBt0qRJqniWL1+u+vj4qDqdTq1fv766adMmFVBPnjxpXGf8+PGqp6enqiiK2qNHD1VVVdVgMKizZs1S/fz8VGtra7VgwYJqq1at1D179hi3++2339TSpUurOp1Obdy4sfrTTz9lqPPmk/Grqqpu3bpVbdCggWpnZ6c6OzurderUUb///ntVVVV1w4YNat26dVVnZ2fVwcFBrVevnrpjxw7jtgkJCero0aNVHx8f1draWi1SpIj65ptvqv/884+qqsmdN11cXNKNqU6dOiqg/vXXX6neW7t2rVqhQgXV2tpaLV68uDp16lST99PqvJnW8eUV5ui8qahq9s628/DhQ1xcXIiIiMDZ2Tk7dy2EyIHi4uIIDAzE19cXW1tbS4eTIbEJeq7ei0JvUHGwscKngANajfLsDYXI4Z7295jRz29pChFCiEyys9HiW8Ahebr1hCSC70djMMiMqEKAJBZCCPFc7B/VVGgUhaj4JEIexGCQ6daFkMRCCCGel4Puv+TiYVwi1yS5EEISCyGEeBGOOitKeCQPAx4Rm8j1B7Fkc9c1IXIUSSyEEDlCbv4wdrK1poS7PQoK4bEJXA+T5ELkTub4vZXEQghhUSlTTsfExFg4khfjbGdNcXc7FBTCYiS5ELlTyt/hi0wFLyNvCiEsSqvV4urqapzq2/7R7KK5kU4DhR0Ubkck8OBhAkkJ8Xi62Oba4xH5h6qqxMTEEBoaiqur6wuNPCqJhRDC4jw9PQGMyUVul5Sg50F0AqHATRstbvY2SG4hcgNXV1fj3+PzksRCCGFxiqJQpEgRChUqRGJioqXDMYvdF0OZuOUCBoNKy/KFGd66nAyiJXI0a2trs8yRIomFECLH0Gq1eWbyp9ZVi6NXrBm48iSLj9zkYaLC1HeqSnIh8jxJLIQQIou8VqUIigIfrzjJ+pM3ACS5EHmeJBZCCJGFXq1cBPgvuVCBaZJciDxMEgshhMhijycXG07eQFVVpr1TFSut3PEv8h75rRZCiGzwauUizOtSHa1GYeOpmwxceZKEJIOlwxLC7CSxEEKIbNKmchG+61YDa63CljO36ffLceIS9ZYOSwizksRCCCGyUauKnvzQvRY6Kw07LoTSe+kxYhMkuRB5hyQWQgiRzZr5FWKRf23sbbT8fekePRYdISo+ydJhCWEWklgIIYQFNChdgJ971cFJZ8WRwAe8++NhImLyxuBgIn+TxEIIISykZgl3lveuh6u9NaeuhdPlh0Pcj4q3dFhCvBBJLIQQwoIqF3NhRe96FHC04fyth3T+/hChD+MsHZYQz00SCyGEsLDyRZxZ2ac+hZ11XAqNotP3h7gZHmvpsIR4LpJYCCFEDlC6kCNrPmxAMTc7Au9F8/b8A1y5G2XpsITINEkshBAihyjuYc/qD+tTsqADNyPieGfBQc7eiLB0WEJkiiQWQgiRg3i52rHmw/pULurCg+gEOn9/iINX7ls6LCEyTBILIYTIYTwcdSzvXZd6Jd2Jik+ix6IjbDt329JhCZEhklgIIUQO5GRrzeKedWhZoTAJSQb6/nKCtcevWzosIZ5JEgshhMihbK21zO9Wg7drFkNvUBm25jQ//n3V0mEJ8VSSWAghRA5mpdUwpUMVPmjkC8BXmy8w7c+LqKpq4ciESJskFkIIkcNpNApfvFae4a38AJi36zJfbjyL3iDJhch5JLEQQohcQFEU+jcvzcQ3K6Eo8MvhEJl2XeRIklgIIUQu0q1uCeZ1qYGNVsOf5+7w7o+HCY9JsHRYQhhJYiGEELnMa1WKsLRXHZxsrTgWHMbbCw5yPSzG0mEJAUhiIYQQuVK9kh6s/agBRVxsuRwaxVvfHeD8zYeWDksISSyEECK38vN0Yn2/BvgVdiI0Mp6OCw+y//I9S4cl8jlJLIQQIhcr4mLH6o/qU9c3eZRO/0VH+PXUDUuHJfIxSSyEECKXc7GzZmmvOrxWpQiJepVBK0/x/d4rMtaFsAhJLIQQIg/QWWmZ27k6vR4NpDVpy7+M++28jHUhsp0kFkIIkUdoNAqjXq/Al6+VB2DxgSA+/Pk40fFJFo5M5CeSWAghRB7zQeOSzOtaHRsrDTsu3KHjwoPceRhn6bBEPiGJhRBC5EGvV/FiZZ96eDjYcO7mQ9p/u19uR80HEg2JBIQFWDSGTCcWN27c4N1338XDwwM7OzsqV67MsWPHsiI2IYQQL6BGcTc29m9I6UKO3IqI450FB9j1b6ilwxJZINGQyIZLG2i7oS0f/PkBMYmWGzAtU4lFWFgYDRs2xNramj/++IPz588zffp03Nzcsio+IYQQL8Db3Z51fRvQsLQH0Ql6ei05ypIDQZYOS5hJkiGJXy//yhsb3mD0gdHciLqBRtFwNeKqxWJS1Ezcj/TZZ5+xf/9+/v777+fe4cOHD3FxcSEiIgJnZ+fnLkcIIUTGJeoNfLnhLKuOXQPAv4EPo16vgFajWDgy8TwMqoFtQdv49tS3BD0MAsDd1p33K71PR7+O2FnZmX2fGf38zlRiUaFCBVq1asX169fZs2cPRYsWpV+/fvTu3TvdbeLj44mPjzcJzNvbWxILIYTIZqqqMn/PFaZsvQjAy+UKMadLdRx0VhaOTGSUqqrsvrabeafmGftSuOhc6FWpF538OmFvbZ9l+86SxMLW1haAIUOG8M4773D06FEGDRrEggUL6NGjR5rbjB07lnHjxqVaLomFEEJYxuZ/bjFk9SnikwyUL+LMD91rUswt6z6QxItTVZWDNw8y9+Rczt4/C4CjtSPdK3bnvfLv4WjjmOUxZEliYWNjQ61atThw4IBx2cCBAzl69CgHDx5McxupsRBCiJznREgYfZYe415UAh4ONix4rya1fdwtHZZIw7Hbx5h7ci4nQk8AYGdlR7fy3fCv6I+LziXb4shoYpGp+q8iRYpQoUIFk2Xly5dn3bp16W6j0+nQ6XSZ2Y0QQogsVqO4G78OaETvJcc4f+shXX84xFftK9GpdnFLhyYeOXvvLHNOzOHgreQv7jYaGzqV60SvSr3wsPOwcHTpy1Ri0bBhQy5evGiyLCAggBIlSpg1KCGEEFmvqKsda/vWZ/iaf9h85hYj1p3h39uRfPFqeay0MsyRpQQ/DGbOiTlsC94GgJVixVtl3qJ3ld54OnhaOLpny1Ri8cknn9CgQQMmTZpEx44dOXLkCN9//z3ff/99VsUnhBAiC9nbWDGva3XK7nRi5o4AFu0P4nJoFPO61MDF3trS4eUr92LvseD0AtYFrCNJTUJBoW2ptvSt2pdiTsUsHV6GZaqPBcDvv//OyJEjuXTpEr6+vgwZMuSpd4U8SW43FUKInOmPM7cYsvo0sYl6ShZw4IcetShVMOs7BeZ30YnRLD63mCXnlhCbFAtA46KNGVRjEH7ufhaO7j9Z0nnTHCSxEEKInOvczQj6LD3OjfBYnGytmNe1Bk3LFrR0WHlSoj6RNQFrWPjPQh7EPQCgcoHKfFLzE2p71rZwdKlJYiGEEOK53IuK56Ofj3MsOAyNAiNal6NPk5IoigymZQ4pg1vNOTmHa5HJA5aVcC7BwOoDaVmiZY49z5JYCCGEeG7xSXpGbzxnHKnztSpFmNKhigym9YJO3DnBlKNTOHf/HAAeth70rdqXt8q+hbUmZ/dpyZLbTYUQQuQPOistX3eoTKWizoz77Tyb/7nF5TtRLHyvJj4FHCwdXq5zLfIaM4/PZHvwdgDsrezpWakn3St0z9LRMi1BaiyEEEI81bGgB/T95QR3I+NxsrVidudqvFSusKXDyhUiEyL54cwPLDu/jERDIhpFw1tl3qJ/tf4UsCtg6fAyRZpChBBCmM2dh3H0++UEx4PDUBQY/HJZPn6pNBqZxCxNSYYk1l9az7envjV2zKxXpB7Daw+nrFtZC0f3fCSxEEIIYVYJSQYm/H6enw8FA9CifCFmdKqGs23O7huQ3Q7cOMDUY1O5HH4ZAB9nH4bVGkaTYk1ybMfMjJDEQgghRJZYc+waX2w8S0KSAd8CDix8ryZlCztZOiyLuxp+lanHprLvxj4gedbRvlX70tGvY47vmJkRklgIIYTIMmeuR/DRsuTxLuxttEx5uwqvV/GydFgWEZkQyfzT81lxYQVJahJWihWdy3Xmo6ofZeskYVlNEgshhBBZ6n5UPB+vOMmBK/cB6FG/BJ+/Vh6dldbCkWUPg2pg05VNzDo+i/txyeegmXczhtYcio+Lj2WDywKSWAghhMhySXoDM3cE8O2uKwBULebCvK418HbPW7dQPuncvXNMOjKJf+7+AyT3o/iszmc0LNrQwpFlHUkshBBCZJtd/4byyepThMck4mJnzYyOVXm5fN67JTUsLozZJ2az/tJ6VFTsrez5qOpHvFv+Xay1ub8fxdNIYiGEECJbXQ+LYcDyk5y6Fg7Ah01LMvwVvzwxBXuSIYk1AWuYe3IukQmRALxW8jWG1BxCIftCFo4ue0hiIYQQItslJBmYtOUCiw8EAVDHx505Xarj6WJr2cBewPE7x5l0eBIBYQEA+Ln5MbLuSGoWrmnhyLKXJBZCCCEsZsuZW3y69h+i4pPwcLBhdufqNCqTu0aaDIsLY/qx6fx65VcAnG2c+bj6x7xd9m2sNPlvRgxJLIQQQlhU4L1o+i47zr+3I1EU+PilMgx8qXSObxoxqAY2Xt7IjOMziIiPAKBDmQ4MqjEIN1s3C0dnOZJYCCGEsLi4RD1jfv1vltQ6Pu7M6lwNL1c7C0eWtkthl/jq0FecCD0BQFm3soyqN4pqhapZNrAcQBILIYQQOcbGkzf4cuNZouKTcLGzZurbVXiloqelwzKKSYxh4T8LWXpuKUlqEnZWdvSv1p9u5bvly2aPtEhiIYQQIkcJvh/NxytO8s/15OaF7vVL8Pmr5bG1tuyAWnuu7WHS4UncjL4JwEveLzGy7kg8HXJO4pMTSGIhhBAix0lIMjBt20W+33sVgHKeTszrWp3ShbJ/rpHb0bf5+sjX7AzZCUARhyKMrDOS5sWbZ3ssuYEkFkIIIXKs3RdDGbbmNPeiErCz1jL2jQp0rOWdLbN/GlQDqy+uZtaJWUQnRmOlWNG9Ync+rPIh9tZ5e8TQFyGJhRBCiBwtNDKOIatOs+/yPQBer1KESW9VztJp2K9GXGXcgXHGzplVC1ZldP3RlHUrm2X7zCsksRBCCJHjGQwqC/deZfq2iyQZVIq52TGjYzXq+LqbdT+JhkQWnV3EgtMLSDQkYm9lz6Aag+hcrjMaJWff/ppTSGIhhBAi1zgZEsbAlSe59iAWRYG+TUsxuEVZbKxe/EP/3L1zjD4w2jhyZqOijRhdbzRFHIu8cNn5iSQWQgghcpXIuETG/3aeNcevA1CpqDOzOlV77o6dsUmxfHvyW36+8DMG1YCrzpURdUbwmu9r2dKXI6+RxEIIIUSutPXsLUauP0NYTCI6Kw2fv1qe7vVLZCoZOHTrEOMOjON6VHKS8qrvq4yoMwJ3W/M2seQnklgIIYTItUIfxjFs7T/sDbgLQNOyBZn6dhUKOT99MrPoxGimHZvG2oC1AHg6eDKq3iiaFGuS5THndZJYCCGEyNVUVWXpwWAmbblAfJIBN3trJr9VmdaV0u4bcejWIcbsH2Mc6KqTXyc+qfkJDtYO2Rl2niWJhRBCiDzhcmgkg1ae4tzNhwC8XbMYo9tWMN6WGpMYw4zjM1h1cRUARR2LMqHhBGp71rZYzHmRJBZCCCHyjIQkA7N2BDB/zxVUFbxcbPnm7SrYOQUzav8oY1+KTn6dGFJziAx0lQXyX2Kx4SOIDTNfeUKIzNNYQd2PwLexpSMRedTRoAcMW3Oa4Afh6Ar9iY37fiB5OO5xDcZR36u+hSPMuzL6+Z13pmy7sguibls6CiFEQpQkFiLL1PZx5+uuTnzy1ySiDcn/861j6jOi7hfU9yph4egE5KXEotVESIy1dBRC5F+3TsPRHyAxztKRiDwqQZ/AvJPzWHxuMSoqrjYFSLj9Nnfu+NBr8Vm614/kszblsLfJOx9tuVHeOfuV37Z0BELkb/96JCcWqt7SkYg8KCAsgJF/jzSOntmuVDs+rfMpGtWeyVsu8MvhEJYeDGb3xbtMe6eq2YcEFxknA6QLIcxD8+h7iiHJsnGIPMWgGlh6bildfu9CQFgA7rbuzGk+h68afYWzjTOOOismvlmZn3vVwcvFlpAHMXT6/iDjfztPbIIkuZYgiYUQwjw0j/6dGOSfuTCP29G36bO9D1OPTSXBkEDTYk1Z98Y6mhdvnmrdxmUK8ucnTehc2xtVhZ/2B9J69l4OPJo5VWQfSSyEEOahaJN/SmIhzGBr4Fbe2vQWh28dxs7KjlH1RjH3pbkUsCuQ7jZOttZ83aEKi3rWpoiLLcH3Y+j642FGrP2HiNjEbIw+f5PEQghhHilNIdLHQryAyIRIRv49kuF7hxOZEEklj0qsfn01Hf06ZniukOZ+hdj2SRPerVccgFXHrtFixh62nr2VlaGLRySxEEKYhyalxkL6WIjnc/T2UTps6sDvV39Ho2j4qOpHLH11KT4uPpkuy8nWmq/aV2bNR/UpWdCBu5HxfLTsBB/9fJzQh3LnUlaSxEIIYR7GzptSYyEyJ8mQxNyTc+n1Zy9uRd+imGMxlrReQv9q/bHWWL9Q2bV93NkysDEDmpfGSqOw9dxtWszYw6qjIWTz+JD5hiQWQgjzkD4W4jncjLpJz609+f6f71FRebP0m6x9Yy3VClUz2z5srbUMa+XHpgGNqFLMhYdxSYxYd4ZuPx4m+H602fYjkkliIYQwj5SmEOljITJoe/B23v7tbU7dPYWjtSNTm0xlfMPxWTYbaQUvZ9b3bcAXr5bH1lrDgSv3aTVrL9/uukxCkiFL9pkfSWIhhDAP6WMhMiguKY7xB8czZPcQIhMiqVKgCmvarqG1b+ss37eVVkPvJiX5c3ATGpb2IC7RwNQ/L/LanL85fPV+lu8/P5DEQghhHtLHQmTA5bDLdNnchTUBawB4v9L7LG6zmGJOxbI1jhIeDizrVZeZnari4WDDpdAoOn1/iOFrTvMgOiFbY8lrJLEQQpiHIjUWIn2qqrL64mo6b+7M5fDLeNh6sLDlQj6p+ckLd9B8Xoqi8Gb1Yvw1tBld6ybfmrrm+HVemr6b1UevYTBI587nIYmFEMI8jH0spK1amHqY8JChe4Yy4dAE4vXxNPRqyNo31tLAq4GlQwPAxd6aSW9WZl3fBpTzdCI8JpFP1/1Dp+8PEnAn0tLh5TqSWAghzEMjd4WI1M7dP0fH3zqyPXg7VooVQ2sO5bsW3z11BE1LqVnCjd8/bsQXr5bH3kbL0aAwXp39N1//8S8xCVITl1GSWAghzEOaQsRjUpo+3tvyHjeiblDUsSg/v/oz/pX80Sg596MnpXPn9iFNeaVCYZIMKgv2XKHljL1sPXtLxr7IgJx7dYUQuYsM6S0eiUmM4bO/P2PCoQkkGhJp5t2MVa+volKBSpYOLcOKutrxffda/NC9FkVd7bgRHstHy07w3v+OcDlUmkeeRhILIYR5yO2mguS7Pjpv7syWwC1oFS1Dag5hTvM5uOhcLB3ac2lZoTA7hjRl4EulsbHSsO/yPVrP+puJm88TGScTm6VFEgshhHmk1FgAGKQDZ37025Xf6LqlK4ERgRSyK8RPrX6iZ6WeGZ48LKeys9Ey5BU/dnzSlBblk5tHfvg7kJem72Hd8ety98gTJLEQQpjH4+3mT9Ra6KOisjkYkZ3ikuIYe2Asn+/7nNikWOoVqcfqtqupUbiGpUMzq+Ie9vzYoxaLetbGt0DyxGZD15zmnYUHOXsjwtLh5RiSWAghzOPxGotH/SxUVeX6oMEE1KrNzS+/lI5vedC1h9d474/3WHdpHQoK/ar2Y0GLBXjYeVg6tCzT3K8QWwc3ZkTrctjbaDkeHEbbefv4fMMZwmRwLUkshBBmktLHAow1FtH7DxD5558ARKxdR+SOHZaITGSRv6//TafNnfj3wb+46dxY0HIBfav1Rfv470IepbPS0rdZKf4a2ow3qnqhqrD8cAjNpu3mf/sC8/XcI5JYCCHMQ3k8sUiusYjcaZpIhP2yPDsjElnEoBpYeHoh/Xf2T57ro2AVVrddnWMGvMpOni62zOlSnVV96lHO04mI2EQm/H6eVrP2su3c7XxZSyeJhRDCPEw6byYnFjFHjgJQ+MsvQVGIOXSIxBs3LBGdMJPIhEgG7xrMvFPzUFF5p+w7LGq1CE8HT0uHZlF1S3qweWBjvn6rMgUcdQTei6bPz8fp+sNhzt3MX/0vJLEQQpiH5rF/J6oeQ0ICCUFBADi1bIl9zZoARO7enf2xCbO4En6Frpu7suvaLqw11oxrMI7R9Udjo7WxdGg5glaj0LlOcXYPb0a/ZqWwsdJw8Op9Xp+7jxFr/yH0YZylQ8wWklgIIczHOMNpEgmBgaDXo3FxwapQQRybNQUgas8eCwYontf24O103dyVoIdBFLYvzNI2S3mrzFuWDitHctRZ8Wnrcvw1tCltH/W/WHXsGs2m7WbeX5eIS8zbg8hJYiGEMB/lv/lCEoKDAbDxKYGiKDg0aQJAzOEjGGJjLRWhyCS9Qc+s47MYsnsIMUkx1PasnetG0bSUYm72zO1SnXV9G1DN25WYBD3TtgXw0rTd/HrqRp4d/0ISCyGE+TxWY5F0JxQA68LJbe+6MmWwKlwYNT6e2NOnLRWhyISI+Aj67ezH/87+D4DuFbrzfcvv8/StpFmhZgk3NvRrwOzO1fByseVmRByDVp7ijW/3se/SPUuHZ3aSWAghzOexqdOTQpMTC6vChQFQFAX72rUBiDl6zCLhiYy7En6FLpu7cODmAeys7JjSZArDaw/H6vFOuiLDFEWhXbWi/DWsGcNb+eGos+LsjYe8+7/DvPe/w3lqgC1JLIQQ5vPYfCFJoXcAsCpU0Pi2fa1aAMQcPZrtoYmM23t9L922dONa5LXkWUnb/Ewb3zaWDitPsLXW0r95afZ+2pz3G/pirVX4+9I9Xp+7j0ErT3LtQYylQ3xhklgIIcznsT4WicamkMLGt+3rJNdYxJ4+jSFBRijMaVRV5aezPzFg5wCiE6OpWbgmy19bjp+7n6VDy3PcHWwY3bYCfw1tRvtqXgD8euomL03fzbjfznE/Kt7CET4/SSyEEObzeI3FnZQai/8SCxtfX7QeHqjx8cSdOWOJCEU64vXxfLHvC2Yen2kcn+KHlj/gbutu6dDyNG93e2Z1rs7vHzeicZkCJOpVFu0PounU3czdeYmYhNw3W7AkFkII80lpf1f1JN29C5g2hSiKIs0hOdDdmLv03NqT367+hlbR8nndzxlVbxTWWmtLh5ZvVCrqws+96rKsV10qFXUmKj6J6dsDaDp1N0sOBBGflHtuUX2hxOLrr79GURQGDx5spnCEELnao6YQNTEBw6MZTbVubiarpAyUFXPiRPbGJtJ07t45Om/uzJl7Z3C2cWZBywV0Kdcl1091nls1KlOATf0bMadLdYq723M3Mp4xm87RfOpuVh4JIVGf8+cgee7E4ujRoyxcuJAqVaqYMx4hRG72qClE/zDSuEjr6Giyil3N5Km0Y0+eQtXnnm9hedEfgX/QY2sPQmNCKelSkhWvraBekXqWDivf02gU3qjqxY4hTfmqfSUKO+u4GRHHZ+vP0HLGHjaevIE+B4+B8VyJRVRUFN26deOHH37A7YlvI0KIfOxRYmGITE4sNPb2KNam1em2fn4o9vYYIiOJv3w520MUyZ00vzv1HZ/u/ZR4fTyNizZm2avLKO5c3NKhicfYWGl4t14J9gxvzqjXK+DhYEPQ/RgGrzpF61l7+ePMrRw5ydlzJRb9+/fntddeo0WLFs9cNz4+nocPH5o8hBB51KM+Fik1Fhpn51SrKFZW2FerCkCsNIdkuwR9Ap/v+5z5p+cD4F/Rn7kvzcXJxsnCkYn02Fpr6dXIl72fNmd4Kz+cba24FBpF319O8Prcfez6NzRHJRiZTixWrlzJiRMnmDx5cobWnzx5Mi4uLsaHt7d3poMUQuQSj/pY6FP6V6SRWADY1XjUz+K4JBbZKTwunN7bevP71d/RKlrG1B/D0FpD0Wq0z95YWJyDzor+zUvz94iXGPhyGRxstJy7+ZCei4/SYf4B9l26lyMSjEwlFteuXWPQoEH88ssv2NraZmibkSNHEhERYXxcu3btuQIVQuQCKU0hD5MTC41z2t+C7R/1s4g5cTx74hIERQTRbUs3ToSewNHakfkt5vN22bctHZZ4Di521gxpWZa/R7zEh01KYmut4URIOO/+7zAdFx5k/2XLJhiZSiyOHz9OaGgoNWrUwMrKCisrK/bs2cOcOXOwsrJCn0ZHLJ1Oh7Ozs8lDCJFHpXTejIoGQOvskuZqdlWqgFZL0s1bJN68mW3h5VfHbh/j3T/eJSQyBC8HL35u8zP1vepbOizxgtwdbBj5ann2Dm+OfwMfbKw0HA0Ko9uPhzl09YHF4srUoO8vv/wyZ54Y1KZnz56UK1eOESNGoNVKdZoQ+dqTTSFOaddYaBwcsC1fnrizZ4k5cRIXL69sCzG/+e3Kb4w+MJokQxKVC1RmzktzKGBXwNJhCTMq5GzL2Dcq0rdZKebvvsKZGxHUK2m5gc0ylVg4OTlRqZLpVLkODg54eHikWi6EyIcedd40RCXPd6BxSb+G0r5mDeLOniX2xHFcXn8tW8LLT1RV5bvT37Hg9AIAWpZoycRGE7GzsrNwZCKrFH6UYBgMqkXHIZGRN4UQ5pPSFBL5qCnEKf3EQjpwZp0EfQIj9400JhXvV3qfaU2nSVKRT2g0lh3c7IXnv929e7cZwhBC5AkpnTdj4pJfPjE41uPsa1QHID4gAP3Dh+neQSIyJzIhksG7BnPk9hGsFCu+rPclHcp2sHRYIh+RGgshhPk86mNhiE+euVRjl/43ZKuCBbEuURxUldhTp7IjujzvTvQd/Lf6c+T2Eeyt7Pn25W8lqRDZThILIYT5POpjocY9Sizsn171bi/NIWZzJfwK7/7xLgFhARSwK8Di1otpULSBpcMS+ZAkFkII89GY1lgoT6mxgP/Gs5AROF/M8TvHee+P97gdfRsfZx9+bvMz5T3KWzoskU9JYiGEMJ+Uu0KMTSH2T13drsajxOKff1ATErI2tjxqe/B2+mzrQ2RCJFULVmVpm6UUcypm6bBEPiaJhRDCfJTkfymG+EQANHZPH6HXxtcXrasranw8cefPZ3l4ec3yC8sZunsoCYYEmnk344dXfsDNViaGFJYliYUQwnxS+lgYE4unN4UoioJdTelnkVkG1cDM4zOZfGQyKirvlH2Hmc1myu2kIkeQxEIIYT7GppDkxOJZfSwA7GukzBsiiUVGJBoS+XLfl/x09icAPq7+MaPqjcJK88KjBwhhFvKbKIQwH601AIaERzUW9k/vYwGPdeA8fhxVteyIgTldXFIcw/YMY8/1PWgVLWMbjKV96faWDksIE1JjIYQwH40Vqh7QJ8+sqMnALMi2FSqg6HTow8NJCAzM4gBzr8iESD7c/iF7ru9Bp9Uxu/lsSSpEjiSJhRDCfLTWGPT/1Tg8q48FgGJjkzzbKRBzXKZRT8u92Hu8/+f7xinPF7RYQFPvppYOS4g0SWIhhDAfjRWGpEeJhZUVio1NhjazMzaHSD+LJ92IukGPP3rw74N/cbd1Z1HrRdTyrGXpsIRIlyQWQgjzeSyxyEgzSAr7lDtDpAOniSvhV+j+R3dCIkPwcvBiaZullHMvZ+mwhHgqSSyEEOajtUZ91BSSkWaQFHbVqoGikBgSQtLdu1kUXO7yz91/6LG1B6ExoZRyKcXSNksp4VzC0mEJ8UySWAghzEdjjSEp+d+K8ox5Qh6ndXJC5+cHQMyJk1kSWm5y8OZBPtj2ARHxEVQpUIXFrRdT2KGwpcMSIkMksRBCmM/jTSHPGM77SSnTqMeeyN8dOHcG76T/zv7EJsVSv0h9fnjlB1xtXS0dlhAZJomFEMJ8tFb/NYXodJna1E5mOmXz1c0M3TOUREMiLUu0ZN7L87C3zlyCJoSlSWIhhDAfjTWqIfmpksnEImWgrLgLFzBER5s7shxvw6UNjPx7JHpVzxul3mBqk6nYaDN2V40QOYkkFkII83lsHIuM3mqawrpIEay8ioBen+/6Waz4dwWjD4xGRaVj2Y5MaDgB7aMp6IXIbSSxEEKYj+a/ppDMJhYADvXrAxB94IBZw8rJFp9dzKTDkwB4r8J7fFnvSzSK/GsWuZf89gohzEdjhWp4lFjoniOxaNAAyB+JhaqqzD89n+nHpwPQu3JvhtcaLnOliFxPEgshhPlo/+tjoXmeGosGDUBRiL94MU+PZ6GqKrNOzOK7U98ByTOUDqwxUJIKkSdIYiGEMB+N9Qs1hVi5uWFbvjwA0QcPmjW0nMKgGvj6yNfGac+H1xpOnyp9LByVEOYjiYUQwny0jzWFWD/fHQ0ODRsCEL1/v9nCyikMqoHxB8ez/N/lAIyqN4ruFbtbOCohzEsSCyGE+WisnvuukBQODZP7WUQdOICqqmYLzdIMqoGxB8ay7tI6NIqGrxp+RUe/jpYOSwizk8RCCGE+LzCORQq7GjVQbG3R371HfMAlMwZnOSlJxYbLG9AoGiY1mkS70u0sHZYQWUISCyGE+Wgf72Nh/VxFaGxssK9dG8gbzSEG1cCYA2OMScXkRpN5reRrlg5LiCwjiYUQwnwev930OZtC4L/mkOh9f5slLEtJSSo2Xt6IRtHwdeOvebXkq5YOS4gsJYmFEMJ8tNao+uSnz3O7aQrHJk0BiD56DH1UlDkiy3Z6g57R+0ebJBVtfNtYOiwhspwkFkII8zEZIOv5+lgA6Er6YlOiBCQmEr0v9zWH6A16xhwYw69XfkWjaPim8TeSVIh8QxILIYT5aB6bK+Q5bzdN4di8OQBRu3a9cFjZSW/QM/rAaH698itaRcs3jb+htW9rS4clRLaRxEIIYT5a8/SxgMcSi717UfX6Fw4tO6QkFZuubEKraPm6ydeSVIh8RxILIYT5aP7rY/GiiYV9jeponJ3Rh4URe/q0GYLLWgbVwLiD40yTCh9JKkT+I4mFEMJ8tNYvNAnZ4xRraxwbNwZyfnOIqqpMOjzJeEupJBUiP5PEQghhPhqtcRyLF7krJEVKc0hkDk4sVFVl+rHprLq4CgWFrxp+JUmFyNcksRBCmI/G2mx9LAAcmzQGa2sSLl8h/mrgC5eXFeadmseS80sAGFN/DG1LtbVwREJYliQWQgjz0ZqvjwWA1tkZh/r1AIj8c+sLl2duP/zzA9//8z0AI+uMpEPZDhaOSAjLk8RCCGE+GmsMKTUWVlZmKdK5VXKzwsOtf5qlPHNZcm4Jc07OAWBozaF0Ld/VwhEJkTNIYiGEMJ/H+lgo1lqzFOn08ktgZUX8xYvEB+aM5pCV/65k2rFpAPSv1h//Sv6WDUiIHEQSCyGE+Wgfm93Uyjz/XrSurjjUS2kOsXytxYZLG5h4eCIAH1T+gA+rfGjhiITIWSSxEEKYz+OdN83438W5dSvA8s0hW65uYcyBMQC8W/5dBlYfiKIoFo1JiJxGEgshhNmoihbUR4mFeVpCAHB8+WXQaon/91+L3R2y9/pevtj3BSoqnfw68WntTyWpECINklgIIczm8aG3zfmZa+XmZpxKPeK3TeYrOIOO3T7GkN1DSFKTeL3k63xe93NJKoRIhyQWQgjzSUw0PlUUg1mLdmnXDoCHv25CNZi37Kc5f/88H//1MfH6eJoVa8b4huPRmLOdR4g8Rv46hBBmoyYlGZ+bO7FwevllNI6OJN68ScyxY2YtOz2BEYH03dGXqMQoahWuxdSmU7HWWGfLvoXIrSSxEEKYzeOJBYp5ZyTV2Nri9KgTZ8Svv5q17LTcirpFn+19eBD3gAoeFZj70lxsrWyzfL9C5HaSWAghzMaYWCgqiiHx6Ss/B9dHzSGRW//EEBtr9vJT3I+9T5/tfbgdfRtfF1/mt5iPo41jlu1PiLxEEgshhNmoicmJhaIBkhLMXr5dzZpYFy2KITqayB07zV4+QGRCJH139CXoYRBFHIrwfcvvcbd1z5J9CZEXSWIhhDCfpORaCkVRQW/+xELRaHBp3x6A8DVrzF5+XFIcH//1MRceXMDd1p3vW36Pp4On2fcjRF4miYUQwmxSmkIUDaCPz5J9uL7dATQaYo4cIf7qVbOVm2hIZNieYRy/cxxHa0cWtFiAj4uP2coXIr+QxEIIYTbGPhYaFfTm72MBYF2kCI7NmgEQvmq1WcpUVZXxB8ez5/oedFod816eR3mP8mYpW4j8RhILIYTZGPtYKEBS1tRYALh17gRA+MaNGOLiXri8uSfnsvHyRjSKhmlNp1GzcM0XLlOI/EoSCyGE2agpfSyysMYCwKFhw+ROnBERPNy69YXKWvnvSn448wMAo+uNppl3MzNEKET+JYmFEMJ8sqGPBYCi1eLasSMAYb8sR1XV5ypnR/AOJh2eBEC/av3oULaD2WIUIr+SxEIIYTb/dd5Us7QpBMD1nbdRdDrizpwh9vjxTG9/7PYxRuwdgYrKO2Xf4aMqH2VBlELkP5JYCCHM5r8BssjSphAAK3d3462n9xctztS2l8IuMfCvgSQYEmju3Zwv6n4hk4oJYSaSWAghzOa/AbLULG0KSeHeowcAUX/9RXxgxqZTvx19m492fERkYiTVClZjSpMpaDVmnONdiHxOEgshhNn813mTLBkg60m6kr44Nm8OqsqDJUueuX5EfAQfbf+I0JhQSrqUZN7L82T+DyHMTBILIYT5pPSxUNQsGdI7LR7v9wQgYsNGku7fT3e9eH08A/8ayJWIKxSyL8SCFgtw0blkS4xC5CeSWAghzMZ05M3sSSzsatXCtkoV1Ph47v/0U5rrGFQDX+z7ghOhJ3CydmJ+i/kUcSySLfEJkd9IYiGEMJuUPhbJI29mT2KhKAoF+vUFIGz5ijRrLWafmM2fQX9ipbFiVvNZlHUrmy2xCZEfSWIhhDAbk9tNsymxAHBs2jS51iI2lvs//s/kvTUBa/jpbHJNxrgG46hTpE62xSVEfiSJhRDCbIydN7N4SO8nKYpCwQH9AQhbsYKku3cB2HdjHxMPTQSgX9V+vFHqjWyLSYj8ShILIYT5mNRYZO04Fk9yaNwY26pVUOPiuPfDD1x8cJFhe4ahV/W8UeoNPqoqA2AJkR0ksRBCmE12TJueHkVRKDhwIJDc12L06j5EJ0ZTx7MOY+uPlQGwhMgmmUosJk+eTO3atXFycqJQoUK0b9+eixcvZlVsQohcxth5MxtvN32cY8OG2DZsAElJtPnjLr4uvsxoNgNrrXW2xyJEfpWpxGLPnj3079+fQ4cOsX37dhITE3nllVeIjo7OqviEELmIJW43fVySIYkfmyViUKDuRZU57n1lrAohsplVZlbe+sT0xIsXL6ZQoUIcP36cJk2apLlNfHw88fH/VYk+fPjwOcIUQuQGptOmZ29TiKqqTD48mU2Gk3hVt+KlE0moc35CbdQKRSOtvkJklxf6a4uIiADA3d093XUmT56Mi4uL8eHt7f0iuxRC5GQmNRbZ23lz2YVlrA5YjYJChU/Ho3FwIO7sWcJXr8nWOITI7547sTAYDAwePJiGDRtSqVKldNcbOXIkERERxse1a9eed5dCiBzOOAmZokJSXLbtd+/1vUw7Ng2AobWG0rzamxQclNyRM3T6dBJDQ7MtFiHyu+dOLPr378/Zs2dZuXLlU9fT6XQ4OzubPIQQeZNx2nQNkJg9icXlsMt8uvdTDKqBt8q8RfcK3QFw69YN20qVMERGEvr119kSixDiOROLAQMG8Pvvv7Nr1y6KFStm7piEELmUycibSbFZvr+wuDAG/DWA6MRoahauyZd1vzTeVqpotRQZPw60Wh5u+YOoPXuyPB4hRCYTC1VVGTBgABs2bOCvv/7C19c3q+ISQuRCJiNvZnGNRaI+kU92f8KNqBsUcyzGzGYzU91WaluhAu7dk2swbo0ajT48PEtjEkJkMrHo378/y5YtY/ny5Tg5OXH79m1u375NbGzWfzMRQuQC2VRjoaoqEw5N4Pid4zhaOzLv5Xm42bqluW7BgR9jU7IkSaGh3BozFlVVsywuIUQmE4v58+cTERFBs2bNKFKkiPGxatWqrIpPCJGLGDtvZnEfi6Xnl7Lh8gY0ioYpTaZQyrVUuutq7OzwmjIFrKyI/PNPIn79NcviEkI8R1NIWg9/f/8sCk8IkZuoen3yE+XROBYGg9n3sff6XqYfmw7A8FrDaVys8TO3satUkYIDBgBwZ8JXxAcGmj0uIUQyGTVGCGE+huTEQkn5z2LmW04vhV3i072foqLydtm36Va+W4a39ej9Afa1amGIjubGwIEYYmLMGpsQIpkkFkIIs1GTUmosHi0wY2LxIO4BH//1MdGJ0dT2rM3ndT/P1MRiilaL14zpaAsWIP7SZW59OUr6WwiRBSSxEEKYjWqssXj0ryXRPB04Ew2JDNszjBtRN/B28mZG0xlYazI/sZh1oUIUmzULrKx4uGULD35aZJb4hBD/kcRCCGE+KTUWVo8+9M1UYzH92HSO3j6Kg7UDc1+ai6ut63OXZV+zJoVHjAAgdOpUHj4xB5IQ4sVIYiGEMBtjjYX1o8TCDDUWGy9v5JcLvwAwqdGkp94BklFu73bDrVty/4ybn44g5vjxFy5TCJFMEgshhPk8qrFQtDaPXr9YjcWZu2eYcHACAP2q9uOl4i+9UHkpFEWh8OcjcWzxMmpCAtf69ScuIMAsZQuR30liIYQwGzXl9lLrR4nFC9RY3Iu9x+Bdg0kwJNDcuzkfVv3QDBH+R9FqKTp1KnZVq2KIiCDEvyfxly6ZdR9C5EeSWAghzCdl5M0X7GORqE/kk12fEBobSkmXkkxqNAmNYv5/Vxo7O7y/X4hthQroHzwg2L8n8Zcvm30/QuQnklgIIczGWGNho0v++Zw1FpOPTObU3VM4WTsx56U5ONo4minC1LQuLhT/6X/oKpRHf/8+wd17EHvmbJbtT4i8ThILIYTZqPqUGovn72Ox+uJq1gSsQUHhmybfUMK5hDlDTJPW1ZUSP/30X81Fjx5E/f13lu9XiLxIEgshhPnoU/pYPF+NxcnQk0w+MhmAgTUGZmi4bnPRurpSfOkSHBrUR42J4VrffoStWZNt+xcir5DEQghhNsYai5TEIhM1Fnei7/DJrk9IMiTxSolX6FWpV1aE+FRaR0e8FyzAuW1bSEri9qjR3BozFkNCQrbHIkRuJYmFEMJ8HtVYKDa2ya8TMzYfR6I+kSF7hnA/7j5l3MowoeGETA3XbU6KjQ1e33xNgYEfg6IQvmoVwe+9R+Lt2xaJR4jcRhILIYTZGGc3tbFL/pmQscRiytEp/HP3H5xsnJjdfDb21vZZFGHGKBoNBfv1w3vhAjTOzsSd/oerb7Qj4vfNFo1LiNxAEgshhPk8SiwU3aPEICHqmZv8fvV3Vl5cCcDkRpPxdvLOsvAyy7FJE3zXrsG2cmUMDx9yc9gwbgwZQlJYmKVDEyLHksRCCGE2xhoLnUPyz2ckFgFhAYw7MA6APlX60NS7aVaG91xsihfHZ/kvFBgwALRaHm75g6ut2xC2evV/t9cKIYwksRBCmI+xxuJRYhGffmIRmRDJJ7s+IU4fRwOvBvSr2i87InwuirU1BQf0x2flCnRly6KPiOD26DEEde5C7KlTlg5PiBxFEgshhNmkrrGITns9VeXLfV8SEhlCEYcifN34a7QabTZF+fzsKlfGd/06Co/8DI2DA3H//ENQ5y7Jc41clLlGhABJLIQQ5pRSY2H79MTip7M/8de1v7DWWDOj2QzcbN2yK8IXplhZ4d6jByW3bMHlzTdBoyHqr78IbN+eG0OGEnvunKVDFMKiJLEQQpiNsc+BrVPyz4TIVOscvnWYOSfnAPBZnc+oVKBSdoVnVtaFC+E1eRIlf9uEU+vWoKo83LKFoA5vE9y9B5G7dkkfDJEvSWIhhDAb1Vhj8WhujydqLG5H3+bTvZ9iUA28UeoN3in7TnaHaHa6UqUoNmsmvhvW4/z666DVEnPkCNf79uNyixbcnTuPhOs3LB2mENlGEgshhPkYO28+Siwe67yZqE9k6J6hPIh7gJ+bH1/W+9Jig2BlBdvy5Sk6bSqld2zH/f330Tg5kXTzFve+/ZYrLVsS3MOfsBUrSLp719KhCpGlJLEQQpiNsfOmXeoaixnHZyQPgmXtxMxmM7GzsrNAhFnPukgRCn86nDJ/78Vr6lTs69cDVSXm8GFujxvPpSZNCer2LvcXLSYuIABVVS0dshBmZWXpAIQQeYixKcQ5+XVCFKgqO6/9xbILywCY2Ggi3s45ZxCsrKKxtcWl7eu4tH2dhOs3iPxzKw+3bSPu9D/EHj9O7PHjAFgVLIhDg/rY16+Pfc2aWBcrlqdqckT+I4mFEMIsVIMBUr592zmnLOV62GVG7RsFQI8KPWhevLllArQgm2JF8ejVC49evUi8dYvI7TuI2ruXmGPHSLp7l4hfNxHx6yYAtO7u2FWtmvyoUhmdnx9WHh4WPgIhMk5Rs7ke7uHDh7i4uBAREYGzs/OzNxBC5ApqYiL/Vq4CQNnDh9DOLEEiKt1rtORs2EWqFKzC4taLsdZYWzjSnMMQH0/syZNE7z9A9JHDxJ2/AImJqdbTFiiAbdky6MqURVe2DDYlSmDtXRyrQgWldkNkm4x+fkuNhRDCLIz9KwC0VqBzYoaDhrNhF3G2cWZqk6mSVDxBo9PhUK8eDvXqAcmJRvyFC8SePk3MqVPEnT9PYsg19PfuEX3vHtEHDppsr9jZYVOsGNbFi2NTrBhWhQtjVbgQ1oUKJT8vVAiNra0lDk3kY5JYCCHM47HEQrHSstPJhWVOyd+mJzaaiJejl6UiyzU0Oh121aphV60a7j16AGCIiSH+8mXiAwKICwgg4fJlEkKukXjzJmpsLPGXLhF/6VK6ZWpdXNAWKIDWzRWta/LDys3N+Fzr6orGyQmNgwNaR0c0Dg5oHBxQbG2lNkQ8F0kshBBm8XiNxY2YW4xyTP5Q8vdqTjPvZhaKKvfT2NtjV6UKdlWqmCxXExNJvHmThJAQEkJCSLxxk6TQUJLu3CEx9A5Jd0JR4+LQR0Sgj4h4jh1rjEmGxsEBjb09Gp0uOeHQ6dDobFBsdMnPbXXG54rOBo3OFsXGBsVKC1ZWKForFGsr0GqNzxWtFh5/bmWFYpX8XLF6tK5GA4oCGg0oGhTNo+caTXLSk/L+U99TJEHKZpJYCCHM4vHEYvi+z4hUoEpcPAMLN7RgVHmXYm2NTYkS2JQokeb7qqpiiIwk6c4dku4/QB8ehj48PPkRlvw8KTwcfVg4hshIDNHRyY+YmOROuAZD8vLI1KOn5jpK+kkHgPL4ek8+f/LnY89TbZeR7R9fhfTWSaucdLZJh8+a1Vi5uz91nawiiYUQwjweSyzOhZ3HGS3TQu9hHZ8HPphyIUVR0Do7o3V2Rlcm49upBgOGmNj/Eo3HEg41IR41Ph5DXPJPNSEeQ3w8anwCanzcY8/jMSTEQ5IeNSkJVZ/02HM9alJiuu+RmPhonSRjgqOqqvE5z3O/gaom/37q9aS1dZ4cSeTxPk/ZTBILIYRZqPrkeTGSHg27N9GhHEX0gRAbbrmgRKYpGg1aRwe0jg6WDiVNJknGE0mHalBBfbQ8JQlJ871Hz5MLfLxwk58mN02qqZ6kWv/xsozbvvB2j7+fev30aF1cnrlOVskziYX/Vn/C4sIsHYYQ+ZZLeCIjAYMC/hX9aXb/QfIbceGWDEvkMUpKE4YmOYOV3hM5T55JLEIehnA3VsbgF8JSCkY8+hal1TCwxkDYNyv5tdRYCJGv5JnEYkazGSQaUg8sI4TIHsr12zD/U2xt7JPHq7B1TX4jVmoShchP8kxiUa1QNUuHIES+Fh9zlauQfOsggJ1b8s+457jVUQiRa8nspkIIs1CTkpKfGBML1+SfUmMhRL4iiYUQwjwMyb3sjTUWDgWTf0ZL3ych8hNJLIQQZqEmPbpvPiWxcCyc/DP6Lhgsd0+9ECJ7SWIhhDCPR8mDscbCvgCgJI8XEHPfcnEJIbKVJBZCCLNIVWOhtQJ7j+TnUaGWCUoIke0ksRBCmMeTNRYAjoWSf0bdsUBAQghLkMRCCGEW/9VYPPZvJSWxkA6cQuQbklgIIczDWGPx2PA4Do8Si8jbFghICGEJklgIIcwiZRIykxoLZ6/knw9vZn9AQgiLkMRCCGEWqj55gCyTGgvX4sk/w4MtEJEQwhIksRBCmEfKAFmax/6tuJVI/hkeYoGAhBCWkGfmCrn77bcYoqItHYYQ+VZC8KNaCavHayweSyxUNXm6ayFEnpZnEovwlatIuis9z4WwNK2T038vXLyTfyZEJc8ZYu9umaCEENkmzyQWbl27YIiWGgshLEprhUu7N/57bW0LTl4QeRPuXwb7OpaLTQiRLfJMYlGgb19LhyCESEuh8smJReh58JbEQoi8TjpvCiGyVqHyyT/vnLdsHEKIbCGJhRAiaxWumPwzVBILIfIDSSyEEFnLs3Lyz5unZPp0IfIBSSyEEFmrUAXQOUNCJNw+Y+lohBBZTBILIUTW0mjBu27y8+D9lo1FCJHlJLEQQmS9ks2Sf178w6JhCCGyniQWQoisV+HR2BbB+2WmUyHyOEkshBBZz7U4eNcD1QBHf7R0NEKILCSJhRAie9Tvn/zz0AKIuGHZWIQQWUYSCyFE9ij3OhSrnXx3yOr3ID7S0hEJIbKAJBZCiOyh0UD7+WDnBjeOw4LGcGoFxIZbOjIhhBnlmblChBC5QIEy8O56WPUehAXCxo8ABVy9wbkY2LmClS1Y24HWOvk9RUn+Cf89N06//sT7QohkzT8HW2eL7Pq5Eotvv/2WqVOncvv2bapWrcrcuXOpU0cmFxJCZEDRGtB3f3InztMr4f4lCA9JfgghzKPRJ7knsVi1ahVDhgxhwYIF1K1bl1mzZtGqVSsuXrxIoUKFsiJGIUReY+cKTYYlP6JC4UEgRFyDhChIjIOkWNAngqoCavI2Kc/VR6+Nz9U0dyFEvmZjb7FdK6qqZuqvsm7dutSuXZt58+YBYDAY8Pb25uOPP+azzz575vYPHz7ExcWFiIgInJ0tk00JIYQQInMy+vmdqc6bCQkJHD9+nBYtWvxXgEZDixYtOHjwYJrbxMfH8/DhQ5OHEEIIIfKmTCUW9+7dQ6/XU7hwYZPlhQsX5vbttEfTmzx5Mi4uLsaHt7f380crhBBCiBwty283HTlyJBEREcbHtWvXsnqXQgghhLCQTHXeLFCgAFqtljt37pgsv3PnDp6enmluo9Pp0Ol0zx+hEEIIIXKNTNVY2NjYULNmTXbu3GlcZjAY2LlzJ/Xr1zd7cEIIIYTIXTJ9u+mQIUPo0aMHtWrVok6dOsyaNYvo6Gh69uyZFfEJIYQQIhfJdGLRqVMn7t69y+jRo7l9+zbVqlVj69atqTp0CiGEECL/yfQ4Fi9KxrEQQgghcp8sGcdCCCGEEOJpJLEQQgghhNlIYiGEEEIIs5HEQgghhBBmI4mFEEIIIcwm07ebvqiUm1BkMjIhhBAi90j53H7WzaTZnlhERkYCyGRkQgghRC4UGRmJi4tLuu9n+zgWBoOBmzdv4uTkhKIoZiv34cOHeHt7c+3atTw7PkZeP0Y5vtwvrx+jHF/ul9ePMSuPT1VVIiMj8fLyQqNJvydFttdYaDQaihUrlmXlOzs758lflsfl9WOU48v98voxyvHlfnn9GLPq+J5WU5FCOm8KIYQQwmwksRBCCCGE2eSZxEKn0zFmzBh0Op2lQ8kyef0Y5fhyv7x+jHJ8uV9eP8accHzZ3nlTCCGEEHlXnqmxEEIIIYTlSWIhhBBCCLORxEIIIYQQZiOJhRBCCCHMRhILIYQQQphNrkosvv32W3x8fLC1taVu3bocOXLkqeuvWbOGcuXKYWtrS+XKldmyZUs2RZp5kydPpnbt2jg5OVGoUCHat2/PxYsXn7rN4sWLURTF5GFra5tNEWfO2LFjU8Varly5p26Tm66fj49PquNTFIX+/funuX5uuHZ79+6lbdu2eHl5oSgKGzduNHlfVVVGjx5NkSJFsLOzo0WLFly6dOmZ5Wb27zirPO34EhMTGTFiBJUrV8bBwQEvLy+6d+/OzZs3n1rm8/yeZ6VnXUN/f/9U8bZu3fqZ5eaGawik+TepKApTp05Nt8ycdA0z8rkQFxdH//798fDwwNHRkQ4dOnDnzp2nlvu8f7sZlWsSi1WrVjFkyBDGjBnDiRMnqFq1Kq1atSI0NDTN9Q8cOECXLl3o1asXJ0+epH379rRv356zZ89mc+QZs2fPHvr378+hQ4fYvn07iYmJvPLKK0RHRz91O2dnZ27dumV8BAcHZ1PEmVexYkWTWPft25fuurnt+h09etTk2LZv3w7AO++8k+42Of3aRUdHU7VqVb799ts0358yZQpz5sxhwYIFHD58GAcHB1q1akVcXFy6ZWb27zgrPe34YmJiOHHiBKNGjeLEiROsX7+eixcv8sYbbzyz3Mz8nme1Z11DgNatW5vEu2LFiqeWmVuuIWByXLdu3eKnn35CURQ6dOjw1HJzyjXMyOfCJ598wm+//caaNWvYs2cPN2/e5K233npquc/zt5spai5Rp04dtX///sbXer1e9fLyUidPnpzm+h07dlRfe+01k2V169ZVP/zwwyyN01xCQ0NVQN2zZ0+66yxatEh1cXHJvqBewJgxY9SqVatmeP3cfv0GDRqklipVSjUYDGm+n5uunaqqKqBu2LDB+NpgMKienp7q1KlTjcvCw8NVnU6nrlixIt1yMvt3nF2ePL60HDlyRAXU4ODgdNfJ7O95dkrrGHv06KG2a9cuU+Xk5mvYrl079aWXXnrqOjn5Gj75uRAeHq5aW1ura9asMa5z4cIFFVAPHjyYZhnP+7ebGbmixiIhIYHjx4/TokUL4zKNRkOLFi04ePBgmtscPHjQZH2AVq1apbt+ThMREQGAu7v7U9eLioqiRIkSeHt7065dO86dO5cd4T2XS5cu4eXlRcmSJenWrRshISHprpubr19CQgLLli3j/ffff+oMvrnp2j0pMDCQ27dvm1wjFxcX6tatm+41ep6/45wkIiICRVFwdXV96nqZ+T3PCXbv3k2hQoXw8/Ojb9++3L9/P911c/M1vHPnDps3b6ZXr17PXDenXsMnPxeOHz9OYmKiyfUoV64cxYsXT/d6PM/fbmblisTi3r176PV6ChcubLK8cOHC3L59O81tbt++nan1cxKDwcDgwYNp2LAhlSpVSnc9Pz8/fvrpJ3799VeWLVuGwWCgQYMGXL9+PRujzZi6deuyePFitm7dyvz58wkMDKRx48ZERkamuX5uvn4bN24kPDwcf3//dNfJTdcuLSnXITPX6Hn+jnOKuLg4RowYQZcuXZ46Y2Rmf88trXXr1ixdupSdO3fyzTffsGfPHtq0aYNer09z/dx8DZcsWYKTk9Mzmwly6jVM63Ph9u3b2NjYpEp2n/XZmLJORrfJrGyfNl08W//+/Tl79uwz2/Xq169P/fr1ja8bNGhA+fLlWbhwIRMmTMjqMDOlTZs2xudVqlShbt26lChRgtWrV2foG0Ru8r///Y82bdrg5eWV7jq56drld4mJiXTs2BFVVZk/f/5T181tv+edO3c2Pq9cuTJVqlShVKlS7N69m5dfftmCkZnfTz/9RLdu3Z7ZSTqnXsOMfi7kBLmixqJAgQJotdpUPV3v3LmDp6dnmtt4enpmav2cYsCAAfz+++/s2rWLYsWKZWpba2trqlevzuXLl7MoOvNxdXWlbNmy6caaW69fcHAwO3bs4IMPPsjUdrnp2gHG65CZa/Q8f8eWlpJUBAcHs3379qfWVqTlWb/nOU3JkiUpUKBAuvHmxmsI8Pfff3Px4v/bt59Q9v84DuDvX9nGLiZ/tkVbO+DgoKzUTgo5UOQ0FzkRclAUF8nNQS47yEEc5UAODsv+XRYKk4UWayU1KTWmYWXP3+H326fvMNvYF596PmqH7fP+vLxfXp9X75eaQM59KcTvqGG6c0Gn04l4PC4ikUjK+kxnY3JNtvfkShaDhVKpFGazWTidTumzRCIhnE5nyl99f7JYLCnrhRBie3s77fqfBkCMjIyIjY0N4XK5hMlkyjnGy8uL8Pv9Qq/X/4Ud5tfDw4MIBoNp9yq3+iUtLy+LiooK0dHRkdN9cqqdEEKYTCah0+lSanR/fy/29vbS1ugzffyTkkPF+fm5cDgcorS0NOcYmZ7z3+bq6krc3t6m3a/capi0tLQkzGazqK+vz/nen6xhpnPBbDYLhUKRUo9AICAuLy/T1uMzvfuZjcvC6uoqVCoVVlZWcHp6ioGBAWg0GlxfXwMAent7MTk5Ka33er0oKCjA3Nwczs7OMD09DYVCAb/f/1MpfGhoaAjFxcXweDwIh8PSKxaLSWte5zgzMwO73Y5gMIiDgwP09PSgsLAQJycnP5HCh8bGxuDxeBAKheD1etHa2oqysjLc3NwAkH/9gP++HW8wGDAxMfHmmhxrF41G4fP54PP5IITA/Pw8fD6f9F8Rs7Oz0Gg02NzcxPHxMbq6umAymfD4+CjFaG5uhs1mk95n6uPfkl88HkdnZyeqqqpwdHSU0pPPz89p88v0nH+3j3KMRqMYHx/Hzs4OQqEQHA4HGhoaUF1djaenJymGXGuYdHd3B7VajYWFhXdj/OYaZnMuDA4OwmAwwOVyYX9/HxaLBRaLJSVObW0t1tfXpffZ9O5XyGawAACbzQaDwQClUonGxkbs7u5K15qamtDX15eyfm1tDTU1NVAqlairq8PW1tY37zh7Qoh3X8vLy9Ka1zmOjo5Kvw+tVov29nYcHh5+/+azYLVaodfroVQqUVlZCavViouLC+m63OsHAHa7HUIIBAKBN9fkWDu32/3uM5nMI5FIYGpqClqtFiqVCi0tLW9yNxqNmJ6eTvnsoz7+Th/lFwqF0vak2+2WYrzOL9Nz/t0+yjEWi6GtrQ3l5eVQKBQwGo3o7+9/MyDItYZJi4uLKCoqQiQSeTfGb65hNufC4+MjhoeHUVJSArVaje7uboTD4Tdx/rwnm979in/+/6FEREREXyaL71gQERGRPHCwICIiorzhYEFERER5w8GCiIiI8oaDBREREeUNBwsiIiLKGw4WRERElDccLIiIiChvOFgQERFR3nCwICIiorzhYEFERER58y+6Pbe/DPN3awAAAABJRU5ErkJggg==",
"text/plain": [
"<Figure size 640x480 with 1 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"# Transcription regulation\n",
"# Question 3.2.1, 3.2.2\n",
"\n",
"dt = 0.01\n",
"t = np.linspace(0,20,int(20 / dt))\n",
"C1 = np.zeros_like(t)\n",
"C2 = np.zeros_like(t)\n",
"C3 = np.zeros_like(t)\n",
"C4 = np.zeros_like(t)\n",
"m = np.zeros_like(t)\n",
"C5 = np.zeros_like(t)\n",
"tau = 10\n",
"\n",
"R = np.where(t <= 10, .1, 10)\n",
"\n",
"C1[0] = tau / (R[0] + 1)\n",
"C2[0] = 1 / (R[0] + 1 / tau)\n",
"C3[0] = tau * R[0] / (R[0] + 1)\n",
"C4[0] = pq(10,1,-2 * tau * R[0]/(R[0] + 1))\n",
"\n",
"print((R[0] / (R[0] + 1) * 2 / (10 * C4[0] + 1) - C4[0] / tau))\n",
"\n",
"m[0] = 1 / (R[0] + 1 / tau)\n",
"C5[0] = R[0] * m[0]\n",
"\n",
"for i in range (0, int(20 / dt) -1):\n",
" C1[i+1] = C1[i] + (1 / (1 + R[i]) - C1[i] / tau) * dt\n",
" C2[i+1] = C2[i] + (1 - (R[i] + 1 / tau) * C2[i]) * dt\n",
" C3[i+1] = C3[i] + (R[i] / (R[i] + 1) - C3[i] / tau) * dt\n",
" C4[i+1] = C4[i] + (R[i] / (R[i] + 1) * 2 / (10 * C4[i] + 1) - C4[i] / tau) * dt\n",
" m[i+1] = m[i] + (1 - m[i] / tau - R[i] * m[i]) * dt\n",
" C5[i+1] = C5[i] + (R[i] * m[i] - C5[i] ) * dt\n",
"\n",
"\n",
"plt.title(\"R increases by factor 100\")\n",
"plt.plot(t, C1, label=\"transcriptional repression\")\n",
"plt.plot(t, C2, label=\"protein degradation\")\n",
"plt.plot(t, C3, label=\"transcriptional activation\")\n",
"#plt.plot(t, C4, label=\"self repression\")\n",
"plt.plot(t, C5, label=\"activated reservoir\")\n",
"plt.legend()\n",
"plt.show()"
]
},
{
"cell_type": "code",
"execution_count": 50,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"[<matplotlib.lines.Line2D at 0x7fd56e430ee0>]"
]
},
"execution_count": 50,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAioAAAGdCAYAAAA8F1jjAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/OQEPoAAAACXBIWXMAAA9hAAAPYQGoP6dpAAA50UlEQVR4nO3de3xU9Z3/8feZyWRyv5ALIRAC4SoKSAEjqBUvK1hapXW9dF0ru6ytFrtttVboRer2Qlt87Hbr+lO7W2/rrlXXW6tWRUBrFUW5CSgIGK65EULuySQz8/39MclAJPfMzJlJXs9H5zGZme+cfA4n47z7/X7P91jGGCMAAIAo5LC7AAAAgO4QVAAAQNQiqAAAgKhFUAEAAFGLoAIAAKIWQQUAAEQtggoAAIhaBBUAABC14uwuYLD8fr9KS0uVmpoqy7LsLgcAAPSBMUb19fXKz8+Xw9F9v0nMB5XS0lIVFBTYXQYAABiAw4cPa8yYMd2+HvNBJTU1VVJgR9PS0myuBgAA9EVdXZ0KCgqC3+Pdifmg0jHck5aWRlABACDG9DZtg8m0AAAgahFUAABA1CKoAACAqEVQAQAAUYugAgAAohZBBQAARC2CCgAAiFoEFQAAELUIKgAAIGoRVAAAQNQiqAAAgKhFUAEAAFEr5i9KCEQrv9+o1edXq8+vNq9fbT6jto7HPr+8PiOf38hnjPz+U39WF88F7n1+I78x8voC9772tsYYGaPAvRT4WYHHCj42pzx/8rHa23X1mml/c1fPdzxWR7tTdPzek49P+fm0tp95/JkWZ+Sl6Zq5Bf381wcwVBBUMGwZY+Tx+lXX3Kba5jbVtbSprtnbft+mplafmlp9amkL3De3+dTc6lNTq/eUn33yeAPBo9V7MoS0tYcQhMa0/DSdNTrd7jIA2ICggiHFGKPqxlaV1rSosr5FVQ0eHav3qKqhVccaPKqq96iqwRMIJs1etfr8EastzmEpPs4hl9Mhl9NSnMMhp8OSwyE5LUsOhyWnZQWe67h3WHJa6vRcp9ctS06H5LAsWZZkyVL7/2RZVvt958ey2tsHX2t/b/uLJ9uf8ppOXor99OdP2faprFN/7PzqZ6/qbnXz2l8+qdKOo7V6eUcZQQUYpggqiDmtXr8OHm/UvsoGfVrVqCMnmnW0pllHTzSptKZFzW2+fm3PYUlpiS6lJbiUlhintASXUhPilBwfp8R4pxJdTiXFO5UYHxe4dzmVGO8M/ux2OeU+JYC4nI6Tj+Pan3M45HCc9lWOXozPPqLvPb1d75VU210KAJsQVBC1jDE6cqJZHx6p1c7SWu2taND+Yw06VN3U67BKTqpbI9PcyklxKzvFrezU9vuUeOWkuJWRFK/0JJfSEuKU4o4L9hYguswdlylJ+vBIjVrafEpwOW2uCECkEVQQNRo8Xn1woFqbSqr14ZFa7Thaq9rmti7bprjjNCE3RROykzUmM1FjMpM0OjNRozMSNSojQe44vtCGgrEjkpST6taxeo8+PFKrc8aPsLskABFGUIFtWr1+vX+gWn/Ze0zvflqtnUdrT+spcTktTc0LTKQ8Y1SqJuSkaGJuinJT3fSCDAOWZWlOYab+vLNcHxysJqgAwxBBBRFV29SmdbsrtO7jSv3lk2Oq93g7vT4mM1HnFmXpc2MzNX10uibnpdA7MszNGTdCf95ZrvdLqqUFdlcDINIIKgi7ljaf1u+u1PNbj2rDnkq1+U72mmSnxOvCybk6b2KWiouyNDoj0cZKEY3mFAbmqWw+eEJ+v2FSMjDMEFQQNnsr6vXYxoN6fuvRTj0nU0am6m+mjdQlZ+Rq5pgMvnjQo2n5aUp0OVXX4tW+Yw2aPDLV7pIARBBBBSHl9xut/bhCj75zQO/sPx58Pj89QVecPVpLZuVral6ajRUi1ricDp1dkKGNnx7X+weqCSrAMENQQUj4/EYv7yjTvev36pOKBkmB9UkuPWOkvjZvnOZPyKLnBAM2d1ymNn56XJsPnND1xYV2lwMggggqGBRjjF7ZWa57Xtuj/ccaJUmpCXG6vrhQf3/uWI3JTLK5QgwFs8cFzvZ5/yALvwHDDUEFA7bzaK3+5cWPtKl91dD0RJeWnT9eN84fp/REl83VYSj53NgMOSzpcHWzKupaNDItwe6SAEQIQQX9VtfSpl/+ebee2HRIxkgJLoe+fkGRbvp8kVITCCgIvdQEl6bkpenjsjp9cOCEFs8YZXdJACKEoIJ+WfdxhX743E6V17VIkpacna/vL5qqfE4rRpjNHZcZCCoHqwkqwDBCUEGfNHq8WvXHXfq/zUckSeOykvSrq2aouCjL5sowXMwuzNRjGw/qgwMn7C4FQAQRVNCr3eV1Wv4/W7T/WKMclnTTBUX67t9M5gJxiKi57RNqPyqrU6PHq2Q3//kChgM+6ejRM5uP6AfP7ZDH69fINLf+/bpZOpdeFNggPyNw0cmjNc3afPCEPj85x+6SAESAw+4CEJ38fqNfvbJbtz+9XR6vXwum5Ojlf76AkAJbdfz9bfz0eC8tAQwVBBWcprnVp1v+Z7Puf2O/JOlbF0/UQzfOVVaK2+bKMNzNnxAIKqeuegxgaGPoB53Ut7Rp2SMfaNOBasU7Hfr1387Qklmj7S4LkCTNaw8qO47UqK6lTWmcDg8MefSoIOhEY6v+/r/e06YD1UpNiNP/3FRMSEFUyc9I1PjsZPmNtOlTVqkFhgOCCiQFQspX//NdbT9Sq8wkl5646dzgWRZANJnH8A8wrETF0M99992nNWvWqLy8XDNnztS9996rc845x+6yho1Gj1dLH3lfu8vrlZPq1v/+U7EmcYVaRIv9G6SSN6XWJsnv1dfrGjQtrkoZO+Ili4XfgIiYuliaeKktv9r2oPLkk0/qtttu0wMPPKDi4mL95je/0cKFC7Vnzx7l5ubaXV5oGSP5fZLfe8rNJ5munvMH2sv0cq/PPPfZx921DfD4/PrXV3bLfaRWFyfG6e4vnKmC5g+lAxH6NzmN6b2JnUyU1zfU7HhK2vJYp6fGSRoXJ6lV0gc21AQMR+kFtgUVyxh7/8tbXFysuXPn6j/+4z8kSX6/XwUFBfrWt76lFStW9Pr+uro6paenq7a2VmlpaeEpsq1Zqi+T6sqkxkqpuUZqqTl531IreeqlthbJ29zFfbPkaw2EDwD9N3mRNPIsyREnWQ499u4BVdW3avGMPE0ZGabPPYCTxn9eKpwX0k329fvb1h6V1tZWbd68WStXrgw+53A4dOmll2rjxo1dvsfj8cjj8QQf19XVhae4jf9P2vrfUl1pIIyEm8MV+I+wI05yOCSrfdVXy5Jk9XKvzzz32cddtz3e2KoTja2yLGlUeqKS4qNlpVmr9yZ2sqK8vqEm9wzpy7+T4uKDTx2o/0gPvV2iCmeBfrVgho3FAQg3W4NKVVWVfD6fRo4c2en5kSNHavfu3V2+Z/Xq1br77rvDX1xLjVT50cnHriQpdZSUMlJKzJASMtrv0wM/u1MlV4IUlyi52m9xCSfvnfGS0yU5nCcDieU8GUwi7I/bS/XPT2yVJP36b2fomjkFEa8BGKgFU3L00NslevOTYzLGyCI8AkOW7XNU+mvlypW67bbbgo/r6upUUBCGL9npV0sFxVJafiCgJKQPmf8nvau0Vnc8vV2SdNMF4wkpiDnnjB+hRJdT5XUt2l1erzNGMfwDDFW2BpXs7Gw5nU5VVFR0er6iokJ5eXldvsftdsvtjsAKqdmTArchpsHj1a3/uzW4LP6Ky8+wuySg3xJcTs2bkKX1uyv1xp5jBBVgCLN1HZX4+HjNnj1b69atCz7n9/u1bt06zZsX2kk7kIwx+uFzO1RS1ahR6Qn6t2vOltMxNHqJMPwsmBK4KOEbeyptrgRAONk+9HPbbbfpxhtv1Jw5c3TOOefoN7/5jRobG/UP//APdpc25Dz9wRG9sK1UToel3351ljKT43t/ExClFkzOlbRLmw+eUH1Lm1JZTh8YkmwPKtdee62OHTumu+66S+Xl5Tr77LP1yiuvnDbBFoNz5EST7v7TLknSbX8zmVVnEfPGZiWpKDtZn1Y16u19VVp0Fou/AUNRVCyhf+utt+rgwYPyeDx67733VFxcbHdJQ4oxRiue2aHGVp/mFGbq5gsn2F0SEBIXtg//bNh9zOZKAIRLVAQVhNcTmw7rr/uqlOByaM3VM5mXgiHj4qmB1avX7a6Qz8+qwcBQRFAZ4kprmvXzlwLrwdyxcKrGZyfbXBEQOsXjs5SaEKeqhlZtPXTC7nIAhAFBZYj7+Usfq7HVp9mFmVo6f5zd5QAhFR/n0CXtvSqvfVTRS2sAsYigMoS9va9KL+0ok8OSfnrlWQz5YEi67MzAmkuv7iqXzZcuAxAGBJUhqtXr16o/Bs7yueHcQk3LZ0EsDE0XTs5RfJxDB483aU9Fvd3lAAgxgsoQ9eg7B7SvskFZyfG67bIpdpcDhE2yO04XTMyWJL22i+EfYKghqAxBtU1tunf9XknS9xdNUXoiC2FhaFt4yvAPgKGFoDIE3f/mftW1eDVlZKr+djYXHMTQd8kZuXI6LO0qrVNJVaPd5QAIIYLKEFNe26KH3y6RJN2xcAoTaDEsZKW4dV778M8ft5XaXA2AUCKoDDH/vm6vPF6/5hRm6pIzcu0uB4iYK2bmS5L+uP0oZ/8AQwhBZQg5UNWopz44LElacflUWRa9KRg+Fp45UvFxDu0/1qiPyursLgdAiBBUhpAH/7JfPr/Rgik5msNFBzHMpCa4dPGUQC/iH7cz/AMMFQSVIaKstln/t/mIJOlbF0+0uRrAHleeHRj+eXF7mfxc+wcYEggqQ8R//qVEbT6j4vEjNLuQ3hQMTxdNzVWKO05Ha5r1/oFqu8sBEAIElSHgeINH/7vpoCRp+UX0pmD4SnA5tXj6KEnSUx8csbkaAKFAUBkCHnnngFra/Jo+Ol0XTMq2uxzAVtfMDawd9PKOMtW3tNlcDYDBIqjEuJY2nx5/N9CbcsuCCZzpg2Hvc2MzNCEnWc1tPr34YZnd5QAYJIJKjPvj9lKdaGrT6IxEXTZtpN3lALazLEvXtveqPPn+YZurATBYBJUYZozRI28fkCTdMK9QcU4OJyBJX541RnEOS9sO1+gTrqgMxDS+2WLY+wdO6KOyOiW4HLpuLtf0ATrkpLp18dTAmir/+94hm6sBMBgElRj2yDuBa/p8edZoZSTF21wNEF3+/txCSdL/bT6iBo/X5moADBRBJUZV1rfo1V0VkqQb54+ztxggCp0/MVtFOclq8Hj17BZOVQZiFUElRj2z+ah8fqPPjc3Q1Lw0u8sBoo7DYenGeeMkSY++c4ALFQIxiqASg4wxwYsPXsvcFKBbV80eoxR3nPYfa9Rf91XZXQ6AASCoxKBNJdUqqWpUcrxTX5yRb3c5QNRKccfpb2ePkSQ99NcSm6sBMBAElRj0ZHtvypdm5ivZHWdzNUB0Wzp/nByWtGHPMX1cVmd3OQD6iaASY+pa2vTyjsBqm9cw7AP0alx2sr7Qfv2f+9/Yb3M1APqLoBJjXtlZrpY2vybmpmhWQYbd5QAx4ZYFEyRJL35YqgNVjTZXA6A/CCox5o/bSiUF1k7huj5A35yZn64FU3LkN9KDf/nU7nIA9ANBJYZU1rXonf2BMxeumMkkWqA/vrlgoiTpmc1HVFbbbHM1APqKoBJDXvywTH4TuDpswYgku8sBYso540fonPEj1Orz67fr9tldDoA+IqjEkBe2B4Z96E0BBuaOhVMkSU99cFglzFUBYgJBJUYcqGrU9sM1cljSYtZOAQZk7rgRumhKjnx+o39d+4nd5QDoA4JKjHip/ZTk8yZmKyfVbXM1QOz6Xnuvyp+2l2pXaa3N1QDoDUElRryys1ySgutBABiYM/PT9aX24dPVL+/mGkBAlCOoxICjNc3acbRWliX9zbSRdpcDxLzvXTZZ8XEO/XVfVfAq5ACiE0ElBry2K9CbMrdwhLJTGPYBBqswK1lfv6BIkvTTFz9Sc6vP5ooAdIegEgM6hn0uO5PeFCBUvnnRBOWnJ+hoTbMeeJOl9YFoRVCJcscbPHr/QLUkaeGZeTZXAwwdSfFx+sHiMyRJD7y5n6X1gShFUIlyr39cIb+RzsxPY5E3IMQWTx+l8yZmyeP16/vPfCi/n4m1QLQhqES5tR8FJvrRmwKEnmVZ+uVXZigp3qlNJdV6/L2DdpcE4DMIKlGspc2nt/cdlyRdckauzdUAQ1PBiCTduWiqJOmXf96tw9VNNlcE4FQElSi2qaRazW0+5aa6NW1Umt3lAEPWDecW6pxxI9TU6tPtT22X1+e3uyQA7QgqUeyNPcckSRdNyZVlWTZXAwxdDoelNVfPUHK8U5sOVOu367loIRAtCCpR7I09lZKkBVNybK4EGPoKs5L1i69MlyTdu36v3tlfZXNFACQbg8qBAwe0bNkyjR8/XomJiZowYYJWrVql1tZWu0qKKgePN+rTqkbFOSydNynb7nKAYeHKs0frmjljZIz0nT9sU2Vdi90lAcOebUFl9+7d8vv9evDBB7Vr1y7927/9mx544AH94Ac/sKukqNIx7DNnXKbSElw2VwMMHz+54kxNyk1RZb1HX//vzWppY9VawE62BZVFixbp4Ycf1mWXXaaioiJdccUV+t73vqdnn33WrpKiyslhH872ASIpKT5Ov/vaHKUnurTtcI1+8NwOLlwI2Ciq5qjU1tZqxIgRPbbxeDyqq6vrdBtqPF6fNn4aOC2Z+SlA5I3PTtZ9f/c5OR2Wnt1yVP/vDZbYB+wSNUFl3759uvfee/WNb3yjx3arV69Wenp68FZQUBChCiNn66EatbT5lZPq1pSRqXaXAwxL50/K1l1fnCZJWvPqHv1h0yGbKwKGp5AHlRUrVsiyrB5vu3fv7vSeo0ePatGiRbr66qt100039bj9lStXqra2Nng7fPhwqHfBdu/sC5xtMH9CFqclAza6cf443bJggiTpB8/t0J93lNlcETD8xIV6g7fffruWLl3aY5uioqLgz6Wlpbrooos0f/58/e53v+t1+263W263e7BlRrV39geGfeZPyLK5EgDfXzhFNU2temLTYX37D9vkdjl08VSuZA5ESsiDSk5OjnJy+jav4ujRo7rooos0e/ZsPfzww3I4omYkyjaNHq+2Ha6RJM2fwGnJgN0sy9LPlkxXXbNXL+0o0zf+e7N+e90sXT59lN2lAcOCbcng6NGjWrBggcaOHat77rlHx44dU3l5ucrLy+0qKSpsOlAtr9+oYEQiV0sGooTTYek3152tL83MV5vPaPn/btFzW4/YXRYwLIS8R6Wv1q5dq3379mnfvn0aM2ZMp9eG86mAHfNTzqM3BYgqLqdDv7n2bCXEOfT05iP67pPbdfREs5ZfNJG5ZEAY2dajsnTpUhljurwNZx3zU+YxPwWIOk6HpV9dNUPLzh8vSbrntU90+9Pb5fGyKBwQLkwKiSInGlv1UVlgXRjmpwDRyeGw9OMvTtNPl5wVXGflut+9q6M1zXaXBgxJBJUo8v6BahkjTcpNUU7q0D6zCYh1N5xbqIeWzlVqQpy2HqrRF/79La37uMLusoAhh6ASRT44eEKSNHd8z6vzAogOF07O0cv/fIFmjElXbXOblj36gX7yx11qavXaXRowZBBUosimkmpJ0txxmTZXAqCvCkYk6emb5+kfzhsnSXrknQNa+Ju/6J39VfYWBgwRBJUo0dzq086jtZKkOYX0qACxxB3n1KovnanH/vEcjc5I1OHqZv3df76n25/arsq6FrvLA2IaQSVKbDtcI6/faFR6gsZkJtpdDoAB+PzkHL3ynQt0ffFYSdIzW47oonve0H0b9qmljTODgIEgqESJDw4Ehn3mjBvBmgxADEtNcOnnX56u5745X7PGZqix1ac1r+7RhWs26JG3SwgsQD8RVKLE+x0TaZmfAgwJs8Zm6pmb5+s3156t/PQEVdR59JM/faTP/3qD/uutT1Xf0mZ3iUBMIKhEAZ/faEt7UGF+CjB0OByWlswarQ13LNDPlpyl/PQEVdZ79LOXPta5v1inu17YqX2VDXaXCUQ125bQx0kfl9WpweNVqjtOU/JS7S4HQIi545z6+3MLdfWcMXpm81E99HaJ9lU26LGNB/XYxoMqHj9CX541WpdPH6X0RJfd5QJRhaASBTa396bMKsyU08H8FGCocsc59XfFY/XVcwr09r7jeuSdA1q3u0LvlVTrvZJq3fXHXbr0jFwtOmuULpycQ2gBRFCJCtsO10iSZo9lfgowHFiWpfMnZev8SdkqrWnWC9tK9dzWI/qkokEv7yjXyzvKFeewNHfcCF1yRq4+PzlHk3JTmGiPYYmgEgU6gsrZYzNsrQNA5OVnJOqWBRN084VF+risXn/6sFSvf1ShvZUN2vjpcW389Lj00sfKSo5XcdEInVuUpbnjRmhSborinEwzxNBHULFZTVOrSqoaJUkzx6TbXA0Au1iWpWn5aZqWn6Y7F03VweONev3jSm3YXakPDlbreGNrsLdFkhJcDp2Zn67powO3M0enaXx2stxxTpv3BAgtgorNth8JrEY7LitJGUnxNlcDIFoUZiVr2fnjtez88Wr1+vXhkRq9++lxvftptbYdrlGDx6vNB08E57hJksOSxo5I0sTcFE3ITdHEnBSNHZGkMSOSlJeWwBw4xCSCis22dwz7FGTYWgeA6BUf59CccSM0Z9wI3Xqx5PcblRxv1I4jtfrwSK12HK3R7vJ61bd4deB4kw4cb9LrH1d22kacw9KojASNzkjUmMxAcMlJdSs7xd1+H6+cVLdS3HHMhUFUIajYrGN+ykyCCoA+cjgsTchJ0YScFC2ZNVqSZIzRsXqP9lY2aF/7bf+xBh2taVZpTbPafEaHq5t1uLpZUnW323bHOZSVHK+0RFfgluBSeqJLaYlxgfuEwPPJ8U4lxjuVFB+npODPTiW54pQY71R8HPNnEBoEFRsZY4I9KgQVAINhWZZy0xKUm5ag8yZmd3rN5zeqrG/R0RPNOnKiWUdONKmy3qOqBo+O1XtU1dCqY/UeNXi88nj9Kq1tUWnt4C6mGOewlBjvVILLqXinQy6npfg4h1xOR/De3X4feM0pl9MKPuewLDkdluIclhwOS07r5L3TITkdDjkdCrZzOiw5rM7tne0/W5IsS7Jktd+3P7Y6XjulTXs7Bdt18f7TttV5u/rMa52OU5fHrovnPtOyr51cn2332e10+/t6eV92SryyUtx9KyLECCo2OnKiWccbW+VyWpo2Ks3ucgAMUU6HpVHpiRqVnqg547pv19zqU1WDRyeaWlXb3Ka6Zq9qm9sCP7e0nfy5uU3NrT41tfrU3OZTU6s38HOrT16/kSR5/Ub1LV7Vt3gjs5MIq+8vmqJvLphoy+8mqNioY9jnjFFpSnAxUx+AvRLjnSoYkaSCEUkD3kar16/mUwKMx+tXq9evNl/gvtXnV5vPnPZcx+OO53zGyOeX/MbI6zPyGyOf38hnjPz+9p/bH/v8p7zul3x+v3xGwXZGRsZIRpJMYJtGgV7twH3Ha6c+NvL7FWynU57vaN/p/Z95b+C5vv2bmS4advXWrrZnumzZXdu+teuqZaKN31EEFRsxkRbAUBMfFxjaSRer6iI0mO1kow/bT02eMSbD3kIAAIhSBBWb+P1GH5XVSZLOGs38FAAAukJQscnB6iY1eLyKj3NoQk6K3eUAABCVCCo22VUaGPY5Iy9VLq7XAQBAl/iGtMmu0sCwz7R8ru8DAEB3CCo26QgqZ+YzPwUAgO4QVGxgjNFH7UM/BBUAALpHULFBZfuS1Q5LmppHUAEAoDsEFRt0TKSdmJuixHhWpAUAoDsEFRvsOtoxP4WJtAAA9ISgYoOdzE8BAKBPCCo2OHlqMkEFAICeEFQirL6lTUdONEuSzhzF0A8AAD0hqETYJxUNkqS8tASlJ3F1UQAAekJQibA95fWSpCl5qTZXAgBA9COoRNie8sD8FIIKAAC9I6hE2J6KQI/K5JEEFQAAekNQiSBjTHDoZyo9KgAA9IqgEkHHGjw60dQmhxVYlRYAAPSMoBJBHb0p47KSleBi6XwAAHpDUImgjqDC/BQAAPqGoBJBwaDC/BQAAPqEoBJBn1QwkRYAgP4gqESI32+Cq9Iy9AMAQN8QVCLk8IkmNbf5FB/n0LisJLvLAQAgJkRFUPF4PDr77LNlWZa2bdtmdzlhsbe9N6UoO1lxzqj4ZwcAIOpFxTfm97//feXn59tdRljtPxYIKqyfAgBA39keVP785z/rtdde0z333GN3KWFFUAEAoP/i7PzlFRUVuummm/T8888rKalv8zY8Ho88Hk/wcV1dXbjKC6l9lYGgMiGHoAIAQF/Z1qNijNHSpUt18803a86cOX1+3+rVq5Wenh68FRQUhLHK0DDGaP+xRkkEFQAA+iPkQWXFihWyLKvH2+7du3Xvvfeqvr5eK1eu7Nf2V65cqdra2uDt8OHDod6FkDve2Kra5jZZljQ+O9nucgAAiBkhH/q5/fbbtXTp0h7bFBUVaf369dq4caPcbnen1+bMmaPrr79ejz76aJfvdbvdp70n2u1vH/YZnZGoxHiu8QMAQF+FPKjk5OQoJyen13a//e1v9bOf/Sz4uLS0VAsXLtSTTz6p4uLiUJdlq45hHybSAgDQP7ZNph07dmynxykpgS/xCRMmaMyYMXaUFDYdZ/wwPwUAgP6x/fTk4YAzfgAAGBhbT08+1bhx42SMsbuMsDjZo8JEWgAA+oMelTBrbvXpaE2zJGkCc1QAAOgXgkqYlVQ1yhgpI8mlrOR4u8sBACCmEFTC7NSJtJZl2VwNAACxhaASZh1BpYiF3gAA6DeCSpgdPN4kSRrPRFoAAPqNoBJmJVWBxd7GZRFUAADoL4JKmB08HggqhVl9uzo0AAA4iaASRrVNbTrR1CaJHhUAAAaCoBJGB9p7U3JS3Up2R83aegAAxAyCShh1BJVxDPsAADAgBJUw6jjjh2EfAAAGhqASRgc6zvhhDRUAAAaEoBJGBzjjBwCAQSGohBFDPwAADA5BJUzqWtp0vLFVEj0qAAAMFEElTA5WBXpTslPilZrgsrkaAABiE0ElTE6emsywDwAAA0VQCZOOM34KCSoAAAwYQSVMDgQn0jI/BQCAgSKohEnwYoSsoQIAwIARVMLkyIlmSVJBZqLNlQAAELsIKmHg8fpUUd8iSSoYwdAPAAADRVAJg6MnmmWMlOhyKis53u5yAACIWQSVMDjcMewzIlGWZdlcDQAAsYugEgaHqwNn/BRkMuwDAMBgEFTCIDiRlvkpAAAMCkElDA6fCPSojOGMHwAABoWgEgZHqjuCCj0qAAAMBkElDE6dTAsAAAaOoBJijR6vqhtbJTFHBQCAwSKohFjHRNr0RJfSElw2VwMAQGwjqIRY8NRkhn0AABg0gkqIBc/4yWDYBwCAwSKohNjhaibSAgAQKgSVEDtyomPohx4VAAAGi6ASYsFTk1lDBQCAQSOohNgRJtMCABAyBJUQqmtpU73HK0kazWRaAAAGjaASQqU1gWGfzCSXEuOdNlcDAEDsI6iEUEdQyc9g2AcAgFAgqITQ0ZoWSQQVAABChaASQmUdPSrpCTZXAgDA0EBQCSGGfgAACC2CSgiVMvQDAEBIEVRCqLS2o0eFoR8AAELB9qDy0ksvqbi4WImJicrMzNSSJUvsLmlAfH6j8lp6VAAACKU4O3/5M888o5tuukm/+MUvdPHFF8vr9Wrnzp12ljRgx+o98vqNnA5Luan0qAAAEAq2BRWv16tvf/vbWrNmjZYtWxZ8ftq0aXaVNCgdwz55aQlyOiybqwEAYGiwbehny5YtOnr0qBwOh2bNmqVRo0bp8ssv77VHxePxqK6urtMtGpw844feFAAAQsW2oPLpp59Kkn7yk5/oRz/6kV588UVlZmZqwYIFqq6u7vZ9q1evVnp6evBWUFAQqZJ7xKnJAACEXsiDyooVK2RZVo+33bt3y+/3S5J++MMf6qqrrtLs2bP18MMPy7IsPf30091uf+XKlaqtrQ3eDh8+HOpdGBBOTQYAIPRCPkfl9ttv19KlS3tsU1RUpLKyMkmd56S43W4VFRXp0KFD3b7X7XbL7XaHpNZQKmVVWgAAQi7kQSUnJ0c5OTm9tps9e7bcbrf27Nmj888/X5LU1tamAwcOqLCwMNRlhd3JNVToUQEAIFRsO+snLS1NN998s1atWqWCggIVFhZqzZo1kqSrr77arrIGjKEfAABCz9Z1VNasWaO4uDjdcMMNam5uVnFxsdavX6/MzEw7y+q3ljafqhtbJUn56QQVAABCxdag4nK5dM899+iee+6xs4xB65ifkhzvVFqirf+kAAAMKbYvoT8UlLUvnT8qI1GWxWJvAACECkElBDqu8ZOXxhk/AACEEkElBMrrAkFlJEEFAICQIqiEQEV7UMlLj771XQAAiGUElRBg6AcAgPAgqIRABUM/AACEBUElBJijAgBAeBBUBsnnNzpW75Ek5XGdHwAAQoqgMkhVDR75jeR0WMpOYTItAAChRFAZpI6JtDkpbjkdLPYGAEAoEVQGKTg/hWEfAABCjqAySME1VNIY9gEAINQIKoPUMfTDGT8AAIQeQWWQKuoCZ/wQVAAACD2CyiCdHPohqAAAEGoElUEqD17nh6ACAECoEVQGqYI5KgAAhA1BZRAaPV7Ve7yS6FEBACAcCCqD0DHskxzvVIo7zuZqAAAYeggqgxAc9qE3BQCAsCCoDEJFPWf8AAAQTgSVQSivbb9qMkEFAICwIKgMQmV7j0oOy+cDABAWBJVBOFYf6FHJTaVHBQCAcCCoDEJle1DJSaVHBQCAcCCoDEJVsEeFoAIAQDgQVAbhGD0qAACEFUFlgJpbfcFVaQkqAACEB0FlgDp6U9xxDqWyKi0AAGFBUBmgYw2BU5Nz09yyLMvmagAAGJoIKgMUnJ+SwrAPAADhQlAZICbSAgAQfgSVAapksTcAAMKOoDJA9KgAABB+BJUBIqgAABB+BJUBqmQyLQAAYUdQGaDgBQm5cjIAAGFDUBkAv9+oqoGhHwAAwo2gMgA1zW3y+o0kKSuZoAIAQLgQVAagsj6wKu2I5HjFx/FPCABAuPAtOwCsSgsAQGQQVAaAU5MBAIgMgsoABM/4IagAABBWBJUBqKRHBQCAiCCoDABDPwAARIatQeWTTz7RlVdeqezsbKWlpen888/Xhg0b7CypTwgqAABEhq1B5Ytf/KK8Xq/Wr1+vzZs3a+bMmfriF7+o8vJyO8vq1fHGQFDJ5qwfAADCyragUlVVpb1792rFihWaMWOGJk2apF/+8pdqamrSzp077SqrT443tEqSslLiba4EAIChzbagkpWVpSlTpuixxx5TY2OjvF6vHnzwQeXm5mr27Nndvs/j8aiurq7TLZJ8fqPqpvagwqq0AACEVZxdv9iyLL3++utasmSJUlNT5XA4lJubq1deeUWZmZndvm/16tW6++67I1hpZyeaWmWMZFlSZpLLtjoAABgOQt6jsmLFClmW1eNt9+7dMsZo+fLlys3N1VtvvaVNmzZpyZIl+tKXvqSysrJut79y5UrV1tYGb4cPHw71LvSoY9gnMylecU5OmgIAIJwsY4wJ5QaPHTum48eP99imqKhIb731li677DKdOHFCaWlpwdcmTZqkZcuWacWKFX36fXV1dUpPT1dtbW2n7YTLO/uq9Hf/9Z4m5aZo7W0Xhv33AQAwFPX1+zvkQz85OTnKycnptV1TU5MkyeHo3CvhcDjk9/tDXVbIVDUykRYAgEixbexi3rx5yszM1I033qjt27frk08+0R133KGSkhItXrzYrrJ6dbwhcGpyFqcmAwAQdrYFlezsbL3yyitqaGjQxRdfrDlz5uivf/2rXnjhBc2cOdOusnoVPDU5mR4VAADCzbazfiRpzpw5evXVV+0sod+ON3JqMgAAkcJpK/10cuiHHhUAAMKNoNJPHT0q2QQVAADCjqDST0ymBQAgcggq/cRkWgAAIoeg0g8tbT7Ve7yS6FEBACASCCr9UN0+P8XltJSWYOsJUwAADAsElX44OezjlmVZNlcDAMDQR1Dph6pGTk0GACCSCCr9EOxRYX4KAAARQVDph+r2HpVszvgBACAiCCr90NGjMoKgAgBARBBU+qGKoR8AACKKoNIPx5lMCwBARBFU+qFj6Ifr/AAAEBkElX4IXucnmaEfAAAigaDSR8YYVTV2zFGhRwUAgEggqPRRY6tPrV6/JM76AQAgUggqfXSivTclweVQUjzX+QEAIBIIKn10oql9DZUkelMAAIgUgkofdVw5OZNhHwAAIoag0kcdPSqZ9KgAABAxBJU+qm5sk0SPCgAAkURQ6aOOybQjklw2VwIAwPBBUOmj4NAPPSoAAEQMQaWPgmf9EFQAAIgYgkofdZz1k8FkWgAAIoag0kcn2ifTso4KAACRQ1Dpo+rgHBUm0wIAECkElT4wxpw864c5KgAARAxBpQ8aPF55/UYSC74BABBJBJU+6JifkuhyKsHltLkaAACGD4JKH1RzajIAALYgqPTBiUYm0gIAYAeCSh8Er5zM/BQAACKKoNIHrEoLAIA9CCp9ELzODz0qAABEFEGlD6rbz/ohqAAAEFkElT44udgbk2kBAIgkgkofnFw+nx4VAAAiiaDSB8EeFYZ+AACIKIJKH5xoCsxRySCoAAAQUQSVXhhjOD0ZAACbEFR6Udfila/9goQZSUymBQAgkggqveiYn5IczwUJAQCINIJKLzjjBwAA+4QtqPz85z/X/PnzlZSUpIyMjC7bHDp0SIsXL1ZSUpJyc3N1xx13yOv1hqukAalpDyoM+wAAEHlx4dpwa2urrr76as2bN0+///3vT3vd5/Np8eLFysvL0zvvvKOysjJ97Wtfk8vl0i9+8YtwldVvNU2sSgsAgF3C1qNy991367vf/a6mT5/e5euvvfaaPvroIz3++OM6++yzdfnll+unP/2p7rvvPrW2toarrH7rCCrpifSoAAAQabbNUdm4caOmT5+ukSNHBp9buHCh6urqtGvXLrvKOk1Nc8caKgQVAAAiLWxDP70pLy/vFFIkBR+Xl5d3+z6PxyOPxxN8XFdXF54CO7bfEVQSGfoBACDS+tWjsmLFClmW1eNt9+7d4apVkrR69Wqlp6cHbwUFBWH9fUymBQDAPv3qUbn99tu1dOnSHtsUFRX1aVt5eXnatGlTp+cqKiqCr3Vn5cqVuu2224KP6+rqwhpWOoZ+0pijAgBAxPUrqOTk5CgnJyckv3jevHn6+c9/rsrKSuXm5kqS1q5dq7S0NE2bNq3b97ndbrnd7pDU0Bcdk2kzCCoAAERc2OaoHDp0SNXV1Tp06JB8Pp+2bdsmSZo4caJSUlJ02WWXadq0abrhhhv061//WuXl5frRj36k5cuXRzSI9Ka2mQsSAgBgl7AFlbvuukuPPvpo8PGsWbMkSRs2bNCCBQvkdDr14osv6pZbbtG8efOUnJysG2+8Uf/yL/8SrpIGhDkqAADYxzLGGLuLGIy6ujqlp6ertrZWaWlpId2232808Ycvy2+kTT+4RLlpCSHdPgAAw1Vfv7+51k8P6j1etV84mcm0AADYgKDSg9r2ibSJLq6cDACAHQgqPahpZn4KAAB2Iqj0oOOMH67zAwCAPQgqPQiuoUKPCgAAtiCo9KCGHhUAAGxFUOlBbccaKlyQEAAAWxBUesDQDwAA9iKo9CA49ENQAQDAFgSVHpy8ICFDPwAA2IGg0oNa1lEBAMBWBJUenOxRIagAAGAHgkoPmKMCAIC9CCrdMMYEV6bNSGKOCgAAdiCodKOlza9Wr18SQz8AANiFoNKNjgsSupyWkuK5cjIAAHYgqHSjYyJteqJLlmXZXA0AAMMTQaUbpwYVAABgD4JKN06uocJEWgAA7EJQ6QZrqAAAYD+CSjdYQwUAAPsRVLrBdX4AALBfnN0FRKv5E7LkdEhzCkfYXQoAAMMWQaUbn5+co89PzrG7DAAAhjWGfgAAQNQiqAAAgKhFUAEAAFGLoAIAAKIWQQUAAEQtggoAAIhaBBUAABC1CCoAACBqEVQAAEDUIqgAAICoRVABAABRi6ACAACiFkEFAABErZi/erIxRpJUV1dncyUAAKCvOr63O77HuxPzQaW+vl6SVFBQYHMlAACgv+rr65Went7t65bpLcpEOb/fr9LSUqWmpsqyrJBuu66uTgUFBTp8+LDS0tJCuu1owP7FvqG+j+xf7Bvq+8j+DZwxRvX19crPz5fD0f1MlJjvUXE4HBozZkxYf0daWtqQ/APswP7FvqG+j+xf7Bvq+8j+DUxPPSkdmEwLAACiFkEFAABELYJKD9xut1atWiW32213KWHB/sW+ob6P7F/sG+r7yP6FX8xPpgUAAEMXPSoAACBqEVQAAEDUIqgAAICoRVABAABRa1gHlfvuu0/jxo1TQkKCiouLtWnTph7bP/3005o6daoSEhI0ffp0vfzyyxGqtP9Wr16tuXPnKjU1Vbm5uVqyZIn27NnT43seeeQRWZbV6ZaQkBChivvnJz/5yWm1Tp06tcf3xNLxk6Rx48adto+WZWn58uVdto/24/eXv/xFX/rSl5Sfny/LsvT88893et0Yo7vuukujRo1SYmKiLr30Uu3du7fX7fb3cxxOPe1jW1ub7rzzTk2fPl3JycnKz8/X1772NZWWlva4zYH8rYdLb8dw6dKlp9W6aNGiXrcbLcewt/3r6vNoWZbWrFnT7Taj6fj15XuhpaVFy5cvV1ZWllJSUnTVVVepoqKix+0O9LPbV8M2qDz55JO67bbbtGrVKm3ZskUzZ87UwoULVVlZ2WX7d955R1/96le1bNkybd26VUuWLNGSJUu0c+fOCFfeN2+++aaWL1+ud999V2vXrlVbW5suu+wyNTY29vi+tLQ0lZWVBW8HDx6MUMX9d+aZZ3aq9a9//Wu3bWPt+EnS+++/32n/1q5dK0m6+uqru31PNB+/xsZGzZw5U/fdd1+Xr//617/Wb3/7Wz3wwAN67733lJycrIULF6qlpaXbbfb3cxxuPe1jU1OTtmzZoh//+MfasmWLnn32We3Zs0dXXHFFr9vtz996OPV2DCVp0aJFnWp94oknetxmNB3D3vbv1P0qKyvTQw89JMuydNVVV/W43Wg5fn35Xvjud7+rP/3pT3r66af15ptvqrS0VF/5yld63O5APrv9Yoapc845xyxfvjz42Ofzmfz8fLN69eou219zzTVm8eLFnZ4rLi423/jGN8JaZ6hUVlYaSebNN9/sts3DDz9s0tPTI1fUIKxatcrMnDmzz+1j/fgZY8y3v/1tM2HCBOP3+7t8PZaOnyTz3HPPBR/7/X6Tl5dn1qxZE3yupqbGuN1u88QTT3S7nf5+jiPps/vYlU2bNhlJ5uDBg9226e/feqR0tX833nijufLKK/u1nWg9hn05fldeeaW5+OKLe2wTrcfPmNO/F2pqaozL5TJPP/10sM3HH39sJJmNGzd2uY2Bfnb7Y1j2qLS2tmrz5s269NJLg885HA5deuml2rhxY5fv2bhxY6f2krRw4cJu20eb2tpaSdKIESN6bNfQ0KDCwkIVFBToyiuv1K5duyJR3oDs3btX+fn5Kioq0vXXX69Dhw512zbWj19ra6sef/xx/eM//mOPF9+MpeN3qpKSEpWXl3c6Runp6SouLu72GA3kcxxtamtrZVmWMjIyemzXn791u73xxhvKzc3VlClTdMstt+j48ePdto3lY1hRUaGXXnpJy5Yt67VttB6/z34vbN68WW1tbZ2Ox9SpUzV27Nhuj8dAPrv9NSyDSlVVlXw+n0aOHNnp+ZEjR6q8vLzL95SXl/erfTTx+/36zne+o/POO09nnXVWt+2mTJmihx56SC+88IIef/xx+f1+zZ8/X0eOHIlgtX1TXFysRx55RK+88oruv/9+lZSU6IILLlB9fX2X7WP5+EnS888/r5qaGi1durTbNrF0/D6r4zj05xgN5HMcTVpaWnTnnXfqq1/9ao8Xe+vv37qdFi1apMcee0zr1q3Tr371K7355pu6/PLL5fP5umwfy8fw0UcfVWpqaq/DItF6/Lr6XigvL1d8fPxpwbm378aONn19T3/F/NWT0bvly5dr586dvY6Lzps3T/PmzQs+nj9/vs444ww9+OCD+ulPfxruMvvl8ssvD/48Y8YMFRcXq7CwUE899VSf/h9OrPn973+vyy+/XPn5+d22iaXjN9y1tbXpmmuukTFG999/f49tY+lv/brrrgv+PH36dM2YMUMTJkzQG2+8oUsuucTGykLvoYce0vXXX9/rhPVoPX59/V6IBsOyRyU7O1tOp/O0mcwVFRXKy8vr8j15eXn9ah8tbr31Vr344ovasGGDxowZ06/3ulwuzZo1S/v27QtTdaGTkZGhyZMnd1trrB4/STp48KBef/11/dM//VO/3hdLx6/jOPTnGA3kcxwNOkLKwYMHtXbt2h57U7rS2996NCkqKlJ2dna3tcbqMXzrrbe0Z8+efn8mpeg4ft19L+Tl5am1tVU1NTWd2vf23djRpq/v6a9hGVTi4+M1e/ZsrVu3Lvic3+/XunXrOv0/0lPNmzevU3tJWrt2bbft7WaM0a233qrnnntO69ev1/jx4/u9DZ/Ppx07dmjUqFFhqDC0GhoatH///m5rjbXjd6qHH35Yubm5Wrx4cb/eF0vHb/z48crLy+t0jOrq6vTee+91e4wG8jm2W0dI2bt3r15//XVlZWX1exu9/a1HkyNHjuj48ePd1hqLx1AK9HDOnj1bM2fO7Pd77Tx+vX0vzJ49Wy6Xq9Px2LNnjw4dOtTt8RjIZ3cghQ9Lf/jDH4zb7TaPPPKI+eijj8zXv/51k5GRYcrLy40xxtxwww1mxYoVwfZvv/22iYuLM/fcc4/5+OOPzapVq4zL5TI7duywaxd6dMstt5j09HTzxhtvmLKysuCtqakp2Oaz+3j33XebV1991ezfv99s3rzZXHfddSYhIcHs2rXLjl3o0e23327eeOMNU1JSYt5++21z6aWXmuzsbFNZWWmMif3j18Hn85mxY8eaO++887TXYu341dfXm61bt5qtW7caSeZf//VfzdatW4NnvPzyl780GRkZ5oUXXjAffvihufLKK8348eNNc3NzcBsXX3yxuffee4OPe/scR1pP+9ja2mquuOIKM2bMGLNt27ZOn0uPxxPcxmf3sbe/9WjZv/r6evO9733PbNy40ZSUlJjXX3/dfO5znzOTJk0yLS0t3e5fNB3D3v5GjTGmtrbWJCUlmfvvv7/LbUTz8evL98LNN99sxo4da9avX28++OADM2/ePDNv3rxO25kyZYp59tlng4/78tkdjGEbVIwx5t577zVjx4418fHx5pxzzjHvvvtu8LULL7zQ3HjjjZ3aP/XUU2by5MkmPj7enHnmmeall16KcMV9J6nL28MPPxxs89l9/M53vhP89xg5cqT5whe+YLZs2RL54vvg2muvNaNGjTLx8fFm9OjR5tprrzX79u0Lvh7rx6/Dq6++aiSZPXv2nPZarB2/DRs2dPk32bEPfr/f/PjHPzYjR440brfbXHLJJaftd2FhoVm1alWn53r6HEdaT/tYUlLS7edyw4YNwW18dh97+1uPpJ72r6mpyVx22WUmJyfHuFwuU1hYaG666abTAkc0H8Pe/kaNMebBBx80iYmJpqampsttRPPx68v3QnNzs/nmN79pMjMzTVJSkvnyl79sysrKTtvOqe/py2d3MKz2XwoAABB1huUcFQAAEBsIKgAAIGoRVAAAQNQiqAAAgKhFUAEAAFGLoAIAAKIWQQUAAEQtggoAAIhaBBUAABC1CCoAACBqEVQAAEDUIqgAAICo9f8Byd5BBATCS68AAAAASUVORK5CYII=",
"text/plain": [
"<Figure size 640x480 with 1 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"# sRNA regulation\n",
"# Questions 3.3.1, 3.3.2, 3.3.3, 3.3.4\n",
"\n",
"dt = 0.01\n",
"t = np.linspace(0,20,int(20 / dt))\n",
"\n",
"m = np.zeros_like(t)\n",
"s = np.zeros_like(t)\n",
"tau = 1\n",
"\n",
"gamma = 10\n",
"alpha = np.where(t <= 10, .1, 4)\n",
"\n",
"\n",
"m[0] = pq(- gamma / tau, (gamma - alpha[0] * gamma - 1 / tau), 1)\n",
"s[0] = alpha[0] / (1 + gamma * m[0])\n",
"\n",
"\n",
"for i in range (0, int(20 / dt) -1):\n",
" m[i+1] = m[i] + (1 - gamma * alpha[i] * s[i] - m[i] / tau) * dt\n",
" s[i+1] = s[i] + (alpha[i] - gamma * alpha[i] * s[i] - s[i]) * dt\n",
"\n",
"plt.plot(t, m)\n",
"plt.plot(t, s)\n",
"\n",
"# FIXME\n"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": []
}
],
"metadata": {
"kernelspec": {
"display_name": "venv",
"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.8.10"
}
},
"nbformat": 4,
"nbformat_minor": 2
}