{ "cells": [ { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "%matplotlib inline\n" ] }, { "cell_type": "code", "execution_count": 41, "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "from scipy.integrate import odeint\n", "import dnaplotlib\n", "import matplotlib.pyplot as plt\n", "from matplotlib import gridspec" ] }, { "cell_type": "code", "execution_count": 45, "metadata": {}, "outputs": [], "source": [ "# Initialize Parts\n", "# Initialize Simulation\n", "# Initial concentration of mRNA and Protein for each repressor\n", "mtet, mlac, mgamma, tet, lac, gamma = initial = [1, 1, 1, 2, 1, 1]\n", "# Non-dimensionalized production rate\n", "alpha = 15\n", "# Degradation Rate\n", "beta = 2000\n", "# Repressor/Promoter Leak\n", "leak = 1\n", "# Hill Coefficient\n", "n = 8\n", "# tetr is orange [1.00, 0.75, 0.17]\n", "# lacI is green [0.38, 0.82, 0.32]\n", "# gamma is blue [0.38, 0.65, 0.87]\n", "plac = {'name':'P_lac', 'start':1, 'end':10, 'type':'Promoter', 'opts': {'color':[0.38, 0.82, 0.32]}}\n", "rbs1 = {'name':'RBS', 'start':11, 'end':20, 'type':'RBS', 'opts':{'linewidth': 0, 'color':[0.0, 0.0, 0.0]}}\n", "tetr = {'name':'tetR', 'start':21, 'end':40, 'type':'CDS', 'opts':{'label': 'tetR', 'fontsize': 8, 'label_y_offset': 0, 'label_x_offset': -2, 'label_style':'italic', 'color':[1.00, 0.75, 0.17]}}\n", "term1 = {'name':'Term', 'start':41, 'end':55, 'type':'Terminator'}\n", "pgamma = {'name':'P_gamma', 'start':56, 'end':65, 'type':'Promoter', 'opts': {'color':[0.38, 0.65, 0.87]}}\n", "rbs2 = {'name':'RBS', 'start':66, 'end':75, 'type':'RBS', 'opts':{'linewidth': 0, 'color':[0.0, 0.0, 0.0]}}\n", "laci = {'name':'lacI', 'start':76, 'end':95, 'type':'CDS', 'opts':{'label': 'lacI', 'fontsize': 8, 'label_y_offset': 0, 'label_x_offset': -2, 'label_style':'italic', 'color':[0.38, 0.82, 0.32]}}\n", "term2 = {'name':'Term', 'start':96, 'end':110, 'type':'Terminator'}\n", "ptet = {'name':'P_tet', 'start':111, 'end':120, 'type':'Promoter', 'opts': {'color':[1.00, 0.75, 0.17]}}\n", "rbs3 = {'name':'RBS', 'start':121, 'end':130, 'type':'RBS', 'opts':{'linewidth': 0, 'color':[0.0, 0.0, 0.0]}}\n", "gammaimg = {'name':'gamma', 'start':131, 'end':150, 'type':'CDS', 'opts':{'label': 'gamma', 'fontsize': 8, 'label_y_offset': 0, 'label_x_offset': -1, 'label_style':'italic', 'color':[0.38, 0.65, 0.87]}}\n", "term3 = {'name':'Term', 'start':151, 'end':165, 'type':'Terminator'}\n", "\n", "lac_repress = {'from_part':laci, 'to_part':plac, 'type':'Repression', 'opts':{'linewidth':1, 'color':[0.38, 0.82, 0.32]}}\n", "gamma_repress = {'from_part':gammaimg, 'to_part':pgamma, 'type':'Repression', 'opts':{'linewidth':1, 'color':[0.38, 0.65, 0.87]}}\n", "tet_repress = {'from_part':tetr, 'to_part':ptet, 'type':'Repression', 'opts':{'linewidth':1, 'color':[1.00, 0.75, 0.17]}}\n", "\n", "def repressilator(y, t):\n", " mtet, mlac, mgamma, tet, lac, gamma = y\n", " \n", " dmtet = -mtet + (alpha / (1 + lac**n)) + leak\n", " dtet = -beta * (tet - mtet)\n", " \n", " dmlac = -mlac + (alpha / (1 + gamma**n)) + leak\n", " dlac = -beta * (lac - mlac)\n", " \n", " dmgamma = -mgamma + (alpha / (1 + tet**n)) + leak\n", " dgamma = -beta * (gamma - mgamma)\n", " return [dmtet, dmlac, dmgamma, dtet, dlac, dgamma]" ] }, { "cell_type": "code", "execution_count": 61, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAsIAAAIcCAYAAAD8JsiGAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAAewgAAHsIBbtB1PgAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzs3Xd4VFX+x/H3SSGFBEIPkkDoTSEUFUQFFetaQUFlbbCr67qoq2vZXVHQdUUX14LruhbEig3F7s8KiggK0nsnAUILSYCEkGTO7487CQlJIGVKMvfzep55wtx759wz4RI+OXPu+RprLSIiIiIibhMW7A6IiIiIiASDgrCIiIiIuJKCsIiIiIi4koKwiIiIiLiSgrCIiIiIuJKCsIiIiIi4koKwiIiIiLiSgrCIiIiIuJKCsIiIiIi4koKwiIiIiLiSgrCIiIiIuJKCsIiIiIi4koKwiIiIiLiSgrCIiIiIuJKCsIiIiIi4koKwiIiIiLiSgrCIiIiIuJKCsIiIiIi4koKwiIiIiLiSgrCIiIiIuJKCsIiIiIi4koKwiIiIiLiSgrCIiIiIuFJEsDsgIhIIo6dtPABEAR5gZ5C7IyLu0RJn4DF/ylXtGwa7M1KWgrCIuEUUEO59tAlyX0TEfaKC3QEpT0FYRNzCgxOCi4CMIPdFRNwjEednjyfYHZHyFIRFxC124owEZ0y5qn1SsDsjIu4wetrGdJyfPZqSVQfpZjkRERERcSUFYRERERFxJQVhEREREXElBWERERERcSUFYRERERFxJQVhEREREXElBWERERERcSUFYRERERFxJQVhEREREXElBWERERERcSWVWHah0UuGHwCicOqeq+SjuEPE/YkUNoKInMTRS4anB7s7LtYSZxAmf0qv6Q2D3RkRcTcFYXeKAsK9jzZB7otIgHiKv+q6rxuigt0BEREFYXfy4ITgIiAjyH0RCZCwRCAcwnTdB5f376HkNxMRkaBREHannTgjYhlTek1PCnZnRAJh9PKN6UAbChvpug8i77SUNmhalojUAbpZTkRERERcSUFYRERERFxJQVhEREREXElBWERERERcSUFYRERERFxJQVhEREREXElBWERERERcSUFYRERERFxJQVhEREREXEmV5USCZebxB4AonFKzqrLlZ43DXknM9jSlcVhmIjMvSg92f9yqcULHxOywCBp7ChOZebz+HvyvJc6gVz5DljUMdmdE6hoFYZHgiQLCvY82Qe5LyAvDU/xV3+8gCsMWf9XfQ2BFBbsDInWRgrBI8HhwQnARkBHkvoQ8D2GJQLiHMH2/g8iD8f49GP09BEYizs8ZT7A7IlIXKQiLBM9OnBGxDIYsSwp2Z0Jd9rSN6UCbbE9Tfb+DKHvJcOfvISxCfw+B4Ew/aYOmX4lUSDfLiYiIiIgrKQiLiIiIiCspCIuIiIiIKykIi4iIiIgrKQiLiIiIiCspCIuIiIiIKykIi4iIiIgrKQiLiIiIiCupoIY7ZRzxVcQNdN2LG+m6FzkKBeEQNHrJ8Pk4ZTWPJXG0U+WpMhlTek3v76NuifjV6Gkbq37dO1XmKpMx5ar2uu6lfph5fJWve2+VucpkMGSZrntxHQXh0JSIU1JTxE103Ysb6boXqQUF4dDkq4/A9FGa1Ce67sWNdN2L1IKCcAjSdAZxI01nEFfSdAaRWtGqESIiIiLiSgrCIiIiIuJKCsIuNnrJ8J6jlwzX9BhxldHTNp4wetpGE+x+iATUzON7BbsLInWRgrC7/RlYOXrJ8GsViMVFngLmj5628SIFYnGRD5l5/NfMPP60YHdEpC5REJZOwCsoEIu79AU+QoFY3OUs4HsFYpHDFISlmAKxuJECsbiRArGIl4KwHEmBWNxIgVjcSIFYXE9BWCqjQCxupEAsbqRALK6lICzHokAsbqRALG6kQCyuo1BTz41eMvwAEAV4gJ3VfHmTahxbHIjHjV4y/CHgzSm9phdW83wSIowx84HEKh7eEueX7qpeoxnW2qNWyxo9bWNtrvvm1Ti2OBD/OnraxvHAJ1Ouam+reT4JEcG+7pl5fG2u+6r2G5xAfBYzj/8GmMCQZT9U81wSQqpx3Vf3moeqXPd+piBc/0UB4d5HmwCcLxnoB3wO7ArA+aRuSqT615svr9FAX/cdca77b4EDATif1E1uu+57An2ZefxPDFmmgQ/3qu51H6jr0ycUhOs/D85FVwRkVPO1TYDYKh6bD/wPeHRKr+nbqnkeqVjGEV9rzBgTAURba/fXtq0qqk6fE6neNVqVY2pz3TfHCRRVkQ08CTw55ar2WdU8j1TMl9d9DGCttQdr21YV1efrvrg/VZEBTASeZ8iyvGqeR8rz2TUPYIxpBOy31np80V4VVLXf1b3mq9O23xhr9SlffTZ6yfB0nN+8tk7pNT2pmq99ERhzjMMUgOswY0wT4F7gWmAzcKa1Nje4vSrLGFNyjVprq3WNVmb0tI2Hr/ur2lfvup+28VvgjGMcpgBchxlj2uJc91cB3wFXWGuLgtursvxx3TPz+JI2GbKsem3OPH4jkHKMoxSA6zBjzPEYcy/WDgOes9beEew+leaXaz4ANCIslVEArh+6AHd7/5wIvGWMGWat1ceYNaMAXD8MAm72/vky4CljzFirkZ2aUgCuHy7G2lHeP//ZGJNurf13UHsUAhSE5UgKwPXbRcB/jDF/UCioFgXg+u0WIA14NNgdqWcUgOu3x40x26y1bwW7I/WZgrAUUwAOHTfihIJ/BLsj9YACcOiY6A0FrwW7I/WAAnDoeMUYs8Na+12wO1JfKQiLAnCI6JrcgNVph4qfPmSM2WqtfTmYfarDFIBDROMO8WRv2Ff8dIoxJsNa+1Uw+1SHKQCHiMZtOpG9dR1AA2CGMeZUa+3SIHerXlIQdrf3gfsVgEPDeSfF8bvfRHLXczuKN73gDQWfB7NfddB/gG8UgEND18vbk7NpH6ve2QjO/2nTjTGnW2sXBblrdc3fgQ8UgENDvxF3smbmO6Qv/A6gEfC5MWagtTYtyF2rdxSEXWxKr+mfBbsP4lt3jmhG2s4Cnn4/E5xlbN41xgyx1s4PctfqjClXtZ8e7D6I7xjg5HtTyd11kC3fbQeI53Ao2BTUztUlQ5a9GewuiO+YsAiGjJ3MFw+PYvf6xeCs1vC5d2RYv+RXg0osi4QQYwz//mMilw9uVLypIfCpMaZDELsl4ldh4YbBj5xIi95NizclAl8YY5oFsVsifhUZHcvQv7xAfKt2xZt64kyTiA5it+odBWGREBMebnjtb204rVdJrZSWOKGgRRC7JeJXETERDH16II3axRVv6gp85C26IRKSYho355x7pxLdqOR3vsE4N9Ap31WRpkZISDLGNMCpIBbqKgy30Q3C+PAfbTl17EZWbM4H6Ax8bIypcwU3xHeMMbFAQrD7EQBNKtoYnRDFOc8O4tNrZ5K3Jx/gFOANY0ydK7ghvmOMiceZEhPqGlW4sVU7ht71Il/842oK8/MARgBbgTpVcKOuUhCWUHUiMDvYnQimJvHhfP5oWwb+aSPbdhcCnIwKboS64cCrwe5EMMUnNWTof07h89E/UJhbCCq44QZ/BiYEuxPB1KJjb4bc+gzfPH4j1lMEKrhRZQrCEhTGGAOchBNY+3J4ZHMvsAiYD8xRYKudtq0a8PnEdpx220ZyDnhABTeCyhgTCQwE+gOpHB7Z3An8CvwC/KK/m9pp3r0JZz5+Ml+NnYMttKCCG0Hl/aTiVKAf0Atn9LYI2IZz3c/V0l+1l9znDE4Z8w9+fOGvxZtUcKMKFIQloIwxjXAKPvwB6FjJYdd4v24zxrwIPGut3VHJscfUqU0DenWIqunL641eHSq+P6JXx2g+eDCZ8+7ZQoETClRwI8CMMa1wwtjvgNaVHDba+3WdMeY54AVrbU5Nz9m0a2PikxrW9OX1Rnxyxe+xzSmtOHV8X364b0HxJhXcCDDvTbq3AtdxjCk7xpiFwH+BV621+TU9Z8su/YhpHPqz4mISKn6PXc4YyYHMDBZNf6p4kwpuHIOCsASMMeZ84AWcZV6q4jjgfmCsMea2mv4HdsmgeCbdnFiTl4aMM/vGMfWe4xj18NbiTSq4ESDGmOtwindUde5uJ2AScLsx5vfW2i9qct4uw1PoPrKy3zXdodNF7TiwI49fJ68o3qSCGwFgjAkH7sSZrlDVFQz6AM8Dtxljrq/pko+9L72FpNQhNXlpyEgddiu5mRms+e5tUMGNY9JdheJ3xvFv4DOqHoJLawK8aoyZ5r0JTmrg6qEJPHZTq9KbXvD+ciJ+YIxpYIx5G5hKzW5gS8JZF3SSdyqR1ECvMV3pNqJ98dPighupQexSSDPGJADf4kxDqckyXj2BucaYW33aMRcxxjBw9EMk9TmjeFNxwY3kIHarzlIQlkB4Fudmhtq6EqdAhD7JqKG/jGzGrcNK1lotLrjRP4hdCkneucDv49y9XVt3ApN90I4rGWM4+d5U2p5RMiOluOBGStA6FaK8U9++Bk6vZVPhODc43ln7XrlTWLhTcKN5h17Fm4oLbrhhVZlqURAWvzLG3IIzH9hXLgb+6cP2XKW44Mbw01Vww88eAX7jw/ZuMcbc7MP2qu2XJ5bxcu/32bU0M5jdqJFKCm58boxpepSXSfW9gnNDnK9MMsac58P2XCUyOpahd72oghvHoCAsfuMNV/64S/tOY8wAP7Rbobufy8CcsZyfV9Z++d1Ww1bR4/p1Fe679entmDOWl3m0uHQVp926kRmza3zPVDnh4YbX/66CG/5ijDkF33wCcqTHgjmKmbkqCxNuaNK5sd/PNePyr3m59/vs3+67Ja+LC240TikpuNENFdzwGWPM1cClfmj6Be9Is9RATOPmnHPPy0cW3HhVBTcO0zdC/OlxnBFHXwsD/uOHdiv069qDhIfBCZWsylBV6bsK2Lm3iD6dKm7n17UHAbjjimY8cF0L7rumOYN7N2TO8lwuG5fGlM/21ur8pUU3CGPGQ8l0b1eymkZxwY3Yo7xMquZZ/POzNQ7nJrqgyFydTeOUeCKiw/16nsL8IrI27iOqSQPiWvv2coxOiOLs/wwiplnJdT8Ip+CGf99UiPOOMD7pp+aTgPv81LYrNEpMYehdLxIRVfI73xU4/z8LCsLiJ8aYdjjTGPylr3fkze9+XXuQbm2jiImq3T+XBavzAOjTuXwQ9ngsi9YdJC4mjH/9oRXjr2/JQ6Nb8d6EZJ74o7PixX9m+PYj6aaNIvji0bYc17xkynVxwQ3Nwa4hY8xpQG8/nuJSY0ySH9uv0IEduRzcm0+zbv4fDd67OhtbaGnWzT9TGYsLbkTEllzmxQU3dENizY2kkiqXPjJGI/e14xTcmIwJK/md73ZjjCrPoeXTxH9G4/9ftG4E5vir8bFPb+eZD5zwuXdfEeaM5SX7Nk7rTEpiA7btLuDxd/bwydx9bM4ooGFMGGf3a8hjN7WibStngQtrLY0vXMW+XA8Adz23g7uec5ZFHtw7lplPtmd12iEOHPQw6PhYwsLK/n98zonOR7m7sn1fIbZtqwZ89ohTcMPbPxXcqJ0b/dx+ODCGAFfR2rMyG4CmpcLpriWZrP90Cxm/7uFARi6eAg9xbRrS5dJ2dL+qI2ER5f/571mVxfLX15Exfxd5u/OJahRJ066N6X51J5JPSyw5BqB5D//d06OCGz73ez+33xS4HKjyEpoFeftZ8tFzbPzpE3L3ZtAosT2pw28julFTPn/wSlKH30af4bcBsGvdItbPnkHGqp85sGc7nsIC4lom0WXwFXQ/9zrCwg9HpYM5mUz7Q3/anXQefa+4k4XTn2T7sjl4igo47vhBnDLmYaIbNWXXusUsnvEMO1bPx4SFk3LSuZx0zTgiGkT7pS2ArUt/YONPn7Jz7a/kZmYAzmhw93OupfPgy0nuc6YKblRAQVj8ZXAAzlHbO5OPamCPGHZkNuLdWTlccHIcJ3ZzBiTCDLRrFcnPK3O54N4t5OR6uGhgHJecEs+qtEO8/V0OPy7LZfGLHWnaKIK8fMsdVzTjzW+yWZt+iL+Nak5khBN2i6dJLFjjjBb3rWC0eMO2QwB0TfbPynG9OzkFN86/VwU3fMCv16TXaQE4RxmZ3nBaepR2/lPLyN11kOY9m5A0qBUFBwpI+yGDnyctJT/7EH3/1LNMGyveXM/Pk5YQHhVO2yGtaZgYw4GMPLb+uIPdyzJLgvDuFc4UoGbdm+BPlRTc2Gqtfd2vJw4x3mkRJwfgVKdTxSB8cN9ePn/wSrK2rqV1z4GknHw+WdvW891Tt9B+gHMPa7OUw9fn/LceI3fvDpp36EVS78EU5O0nbdF3/Pz6w+Tvz6bviMMDp3s2OwMiBbn7+HjcpbTuMYDOgy9n65Lv2fzL/2HCwmnbbyhzpowjKXUInYdcweafv2D1N9OIadycPpf/2S9tAcx58e80iG1Ei469iek3lPx9e9k8/ytm/+9urKeILmeMVMGNCigIi895P2LsE4BTtTfGNLXW+uU29quHJrB5RwHvzsrhtuHNSkZmwZnve+HfthAbHcacZ9rTJflw5bpnPtjD2Kcz+N/He/nrqBbERocx/vqWTPs2m9bNInj4d63KnevXNc784COnTezJLuTeF5zR47+M8F+1pLP6qeBGbRljmgNtA3CqvgE4RxnFo7RNS02NGDS+L42S48oc1z+3kBnDv2btR1vKBOENn6cx79HFtExtyuCJJ5WZ+1twoID8fQWHz7XSG7q7+3+VpwoKbrzsLbjxtd9PHjp6EZgsUeXVKH547i9kbVvHqTc9RufBl5dsXzj9qZIA2DSlR8n2Qb9/hEaHV1YAoP/BA8y4+zzWfj+9TBDO3LQSgJ1rF3L2PVNI7HYSAKnDxvLO2FNJ+/Vrdq1byIUPTqdJUhcAepx7He/edjpbl8wuE1592RbAuX97rdz76HP5n5l+x5ms+346Xc4Y6W1fBTdKUxCuR4wx83GW/SkR0zw60YSB9ZD48m6THqSuHSkMZwHvQFhpjCmoYLtPhk8XrXMCauoRN7jd+/wOdmUVMe/ZtmVCMMAN5zXh1skZzPfOCQY4kOdh3dZDnHti2eBQrHhE+IcluWzeUYDHA1t2FvDRnH0UFlleuus4zj2p4tf6ytVDE9i6u5C7/1dSzXqKMeYxoMblTr2Kr9lEY3xzjcYktEg0YeFYT1Hiy1fvqivXfaCKvTQxxmwDPBXs88vNjpmrs4k7LpaoRoff4pEhGCAyNoJGbePY/vPOkm35OYeYO3ExsS2jGTr5lDJtAEQ2jCSyYSQARQUestbl0KBRZMDKQ/ca05XcHXmsemcjOP8nfmmM2QVU9HOlOnx+3bduGpEYFgYeD4nbM+vMz/tAzd1NPcr3Mb74D3s2LiN94Xd0GzqqTAgG6HT6cBZNf4qouATimh1Xsv3I8AgQGd2QRq3bs31Z2dl3md5R3BOvvrckuAJExsQR1yKJzM0rGHjDgyXBFSCueRsaxMRTkH/Ab21V9j5im7QktklLDmSW/FwvKbiRu3cH6YtmgvP/9ULvdV+bOXi+uuYzrLUBW99eQbh+SeSIymx5uw8W/zH8yH0u0dKfjS9cd5DjmkfQssnhfyqZOYW89W02LZuE89m8/Xw2b3+Fry09wXbx+oN4PFS4YoS1loXewP3yF1ll9sVGG96f0NbvIbjYX0Y2442vs1m8vuS68uUwtM+u0bysXT5vs55pfexDfCM/5xD7t+XS9szDwcF6LOs/3cL6T9PIXJ1Nfs6h4rm2AMQddziPr/8kjfysQwz4a+9yIfhIe9dk4ym0ARkNLmaM4cQ7e7Hpq20c3JsPYPDtzxWfXaPbMwt93mY9YqjCe/YGO064qPzy9VENnU80mrY7PBpsPR7Wz57B+h8/JHPzCvIPZGOLSr7PxDUve8o9m1cS3iCajqddVq79A3u2Eds0kaQ+Z5bZXnDwAIdyc2jRuewHpb5sq6ggn9Xfvs2meZ+Rlb6WQ7k5WM/hTJvYvezslbDwCAb9/hHevfV0PEUF4FxTZQbaaqFeXZ8KwvVLxpEboptGtTRhJsx6rOdgZv7Oil7kUg2o5V3M+/OKWLf1EBecXDaEfrcwlyIP7NxbxIRXdlXyakhqHlny54VrK18xYm36IfblehjYM4Y5zzh1LbL2F/HGV1mMnZzBqIfT2fRWZ+Ji/L/C0wMv7yodggF2AYdq2WxLnE8JPIBPrtHoRs1amrCwMOvxeA7m7NF1f1gsTklynymZquCdFlGYX8TXY+ewfd4u4pMaknRqK2JbRBMWGcah/YWseH0dTbsenkKx7SdnJKrtGceVb7yyc/l5fnBp1lrmPLSwOASDc53uAgorf1WV+Py6b9kkvGW4MWFF1np27i3SdX9YPN5PIfdsXkHj1h2Ia1E+h+VlOz+vi+cHFx7K5+t/jWH78jnEt2xLUuoQYhNaEhYRyaG8faz4/OUyobnw0EFytm+kVbeTiIwu+4nF/l1byd+fRcdTL+PIBUgyt6wCoGm77n5p62BOJv/3yDVkbl5JQnIX2p14LtGNmhIWHsH+3VtZ9/30MseDE5xn/efPxSEYnJHgnVT8SVNV+eqaL5d1/ElBuB4J5EcF9Z0xZhAwuzZtLF6Xj7Xlw2vaLucHxwcPJXPpqVWbAbLIGy77dCr/KeIC7/zgvp0P70uIC+eWy5oxe1kub32bw9vf5jDmN/4NB89/nMlDr5UEewtcaa19x68nFZ8yxlwDvOrLNjNXOStGFN8ot+yVtWyft4vUm7qRenP3Mv9Rr/t4C1B2dYkDO/IwEYaGrY79CfrhG+UCNyL86+TlrP9kS/HTg8BQa+2PAeuA1Jox5n68K6kUHTpIw2YVf2CSsWIecDhELvv0BbYvn0PqsLGkDr+97LX8wwfOsaXmEu/dsgrrKaJZqW3F9mxaBkCz9j3L7cvc5EyBaJZyvF/amv/WY2RuXsmgGyfSZUjZqu6L3p/sfR+H27IeDz88dxcZK+YWb9oNnGKtXVvuhC6gdYRFKrFkgxNQe3csG4SLVxXbXY3lzJaszyc+NowOx0WW2/fr2spXjLjhPCf8vjMzu8rnqomP5+zj5ie3l950h0KwQOkb5ZxwuvVHZ7DmhDFdywQHT5Fl1dvrAcqsN2wt2EJLfs6xP1jIDMDSaaWtfHs9S15aU/zUAlcrBNd/+QfK/7z0FBWy4oupwOER4a2LZwFwwsV/LHste4pY9dVrZY4FZ7T5yG0l+zYV7zv+KPsOh15ftrV18SxiElqUC8EFeftZO+vdcsf/Mm0iG3/6pPhpHnChW0MwKAiLVGpPjhN042LK/jM5ubsz/3Hy+3vYn1c+DKfvKihZ8qx0Ww2jw8p9zAWlRoS7lA/CZ/RpSKOGYXy3KJes/b5fRxhg3opcRj6YhufwB2KPW2v9VSVK6pnMVVlEN4kqGdENi3Sm6ORsOTw3vuhQET89vJBdS50R3dIjwsWhdsWb68u1nbNlP54i5xdLT4GHvWtziIyLID7Z/zfKbf5mK3MfWVx601hr7Qd+P7H4VXTj5mRuXknm5pUl2zxFhcx56T6ytq4lvEE0jY5zpqCFRThz1nMyNpYcW1SQz09T7mfXukVA2RHhzE1HC6/LwJhyUxCKXxcZE0d8qZvZfNlWWEQDDh3I5mDO4QWUDuXmMOs/f2b/rnTCIxuQcFwnAJZ/PoXln75Y8q0BRlpr55U7kYtoaoRIJYpvbBv7dAaXD84lKtJwVt+GnNarIZedFs8HP+yjw9VrufiUeFokRJCRWcj6bYf4cVkus55MocNxh28M6tM5mvdm5XDpfVvo3TEaY+DOEc2Ijw1n4do8oiINPVPKB+HICMO5/eN4d1YOn/60j1Fn+3akbG16Phf+bQt5+SU3Or0F3O3Tk0i9VXiwiOzN+2l90uHp9u3PaUPGL7v48qbZpJybhOeQh21zdxLbIpqY5lF4Ci1xiYdvljthdFc2f7ONRf9dyfafd9G8ZxMK8wrJXJ3N/q0HuPJbZ13XvetzKDrkoUWvphX+wuhLOxbuYdZffyl9R+tEa23AyraL/6ScdB6rvnqdzx4cSfuBFxEZHUv6opmEN4gmIiqWJsldCfNWV2t/8vlkrPiJLx+5lpQBF+IpzGfb0h+JbdKSmIQWeAoLy6wukbl5RZkgXdqejctp1KodDWLjy2z3FBawN30NLTqllrmufdlWyoDzWfbx83xy/zDanngOh/Znkb74e+e4sHASkroQFhHJxrmf8fPrD5du8mZr7cfV+f6GIo0Ii1TiNwPjeeT3LSnyWCa9vZsJr+wiw3vX9jsPJDPp5lYc1yySN7/J5sn39vDdwgPExYTx1J8SObFr2fmQT/0pkeGnN2L20lwefHUXj7+zh4bRYazfeois/R6Obx9VUmTjSJcMcn4Yzvhxn0/f347MQs67Z3PpKR4zgeuttbW5WUJCyN412diisuWOu17Rnv63H094VDhrpm9kx6+76TIshaGTT+FgZn6ZG+UAGreL48LXh9D+vCSyN+5jxZvr2PzNNsKjwun/58Mf/e5ZEZj1g7M27uPrW+dQlF9ymb8G/M2vJ5WASUo9g34j/0JkTDzrvp9O2sJvaT/gAk4Z8w8K83Np1f3EkmO7Dh1F/6vuIbxBNGu+ncaOVfPpcuZIht71Igez95QZkbUeD3vTVtO0bbeSIF0sL3sXeVk7K5zKkLV1HZ7CQ2WmJviyLYC+l/+Z4y+8kaKCfFZ99TqZm1eROvxW+o28C+spomm7HmSsnMf3z97hzFVyPGStfb6q39dQZlRFVUJR6Zvl7hzRjEk3+2pVmNCwP6+IIbdvKpmWASwDTrPWZh3lZVLHlb5ZbsDfetN9ZMcg96huyd2Vx6fXzmL/ttziTV8Dv7HW1nZlFAmi0jfLnX33FJJSh5Q7Zu7U8az88lUuevgjmrcvHzJD2d70NXw2fgSHcnOKN70MjLEKgICmRoi4TkGhZcSE9NIhOB04XyFYQtmh/QV8dcuc0iF4ETBcITi0HMoru667tZbV305j5Zev0rbf2a4LwQf2bOfLiTeUDsFfADcpBB+mICziItZabnp8G58fLgKSBZxnra0rVapEfK6owMN3d84jc3XJagKbgQustTlHeZlY53aWAAAgAElEQVTUQz/89y9smvsp8S3bUnBwPzvXLmTvllU0aduNQTdODHb3Air/QA5fPnYDuZklKwItAK6w1ta2amJIURAWcZEHXt5VunrdIeBSa+3yIHZJxK+stfz4wAK2zS1Z3z8T55e/7Ud5mdRTzVJ6sHPtr6Qt/I6wiEgat25P/6vuofu51xHRoPwNyaGqqCCfb5/4A1lpJcsDbsCZBlRxKVQXUxAWcYkKCmZcY62dFcQuifjdgqeXs/7TtOKnB4GLrbWrgtgl8aPUYbdWOEfYTSopmHGetXZHELtVZykIS8h7d1bOkWWDQ9Jlpzbij5c2rXCfCma4z8o317Pl29Af9Ow2oj3tzipfThdg5VvrWTpFBTPcZMHbk1j+2UvB7obf9bn8z7Ts0rfCfSqYUT0KwhLytuwoYMuO0J8S1TMlqsLtKpjhTtmb9pO9KfQ/BU0+veIVYTZ/s5W5E1Uww20yvRXbQl33c6+vcLsKZlSf1hEWCWEqmCFupIIZ4kYqmFEzGhGWUDUXiD/mUfXficC3Fe1QwQxXegtww6jnFcCUinZkbcg5smDG66hgRqibCPw72J0IgLuA+yvaoYIZNacgLCHJWlsEhPznwsaY3Iq2788r4jd/3cyGbSVTQpYCl1lr8wPWOQk477JIIT8PyBhT4aT/3F15fHXLHA7llHwLvkaFA0Kedy3okF8P2hhT4c/vvelr+Obxm/AUlnwLXgYeCFjH6jlNjRAJMZUUzLhABTMklKlghriRCmbUnoKwSAhRwQxxo0oKZvxGBTMklKlghm8oCIuEEBXMELexUFnBjG3B65WIf3kKVTDDVzRHWCREfPjjPjZllAwEqGCGuMKKN9axf2vJdAgVzBBX+OXNR9m/q6RQjApm1IJGhEVCRKkQDCqYIS5RKgSrYIa4RqkQrIIZtaQgLBJ6VDBD3EgFM8RtVDDDBxSERUKLCmaIG6lghriRCmb4gNEKGyL1lzGmEfAn4AYgDThfawVLqDPGtAZuBa7FKShzrZaLklBnjOmMMbdj7VXAM9baCotrSPUoCIuEAGOMARooBIubGGMicP4f03JR4hrGmCjgkH758w0FYRERERFxJc0RFhERERFXUhAWEREREVdSEBYRERERV1IQFhERERFXUhAWEREREVdSEBYRERERV1IQFhERERFXUhAWEREREVdSEBYRERERV1IQFhERERFXUhAWEREREVdSEBYRERERV1IQFhERERFXUhAWEREREVdSEBYRERERV1IQFhERERFXUhAWEREREVdSEBYRERERV1IQFhERERFXUhAWEREREVdSEBYRERERV1IQFhERERFXUhAWEREREVdSEBYRERERV1IQFhERERFXUhAWEREREVdSEBYRERERV1IQFhERERFXUhAWEREREVdSEBYRERERV1IQFhERERFXUhAWEREREVdSEBYRERERV1IQFhERERFXUhAWEREREVdSEBYRERERV1IQFhERERFXUhAWEREREVdSEBYRERERV1IQFhERERFXUhAWEREREVcK+SBsjPmHMcZ6H/fVsI2BxpgPjTG7jDF5xpgVxphxxphoX/dXRERERAIjpIOwMaY7cFct2xgF/ABcDOQDK4FOwIPA98aY2Nr2U0REREQCL2SDsDHGAP8DCoBva9hGCvASEA7cDSRba/sCnYHVwInAYz7oroiIiIgEWMgGYWAMcBrOyG1aDdu4C4gCvrTW/staawGstZuB0d5jbjTGtKptZ0VEREQksEIyCBtjWgCPAiuAJ2rYhgEu8z596cj91to5wCogErikZj0VERERkWCJCHYH/OQJoCkwzFpb4GTaamsLtPb++cdKjvkR6AacDDxfncaNMUmV7GrgbXMnsAsoqk67IiIiIiEiHGjh/fNSa22+r08QckHYGHMWMAp43Vo7qxZNdfZ+zQe2VXLMhiOOrY6aTtcQERERcZsTgfm+bjSkpkZ4lzN7DsgG/lLL5pp4v2YVzw2uwN4jjhURERGReiLURoTvw1na7E/W2h21bKt4jeBDRzmmeIg+pgbtJ1eyvQ0wF+Dnn3+mdevWlRwmIiIiErq2b9/OSSedVPx0lz/OETJBuNSawb8C//VBkwe9Xxsc5Zgo79e86jZurU2vaHvp+cytW7cmKamyqcQiIiIiruGXe6ZCaWrEszjB/mZrrccH7RVPe0gwld9t1+SIY0VERESkngiZEWGgD2CBjyrIrY29X+8xxvwJSLPWnniM9tZ6v0YBxwFbKzimwxHHioiIiEg9EUpBGJxlNo5W3CLO+zh4lGOKbQEygERgEPBOBccM8n6dV40+ioiIiEgdEDJTI6y1CdZaU9EDeMV72DjvtpQqtGeBD7xPxxy53xhzCs56vwXAR755FyIiIiISKCEThGvKGHO7MWaTMeatCnb/C2fViHOMMXcVzxU2xrQDpniPedFamxGg7oqIiIiIj7g+CAMJQDucKRBlWGs3Ar8HPMBjQJox5lecOcFdgQU4K1WIiIiISD2jIHwM1tpXgdOAT3DWC+6BU1FuPHCqtfZA8HonIiIiIjVlKi+aJsFgjEnCW345LS1N6wiLiIiIK6Wnp5OcXFJ/LLmyGgy1oRFhEREREXElBWERERERcSUFYREfefLJJxk/fjxZWVk+a3PmzJmMHz+emTNnVnrMkCFDMMaUeTRo0ICkpCSGDRvG119/7bP+iIiIhBIFYREfefLJJ5kwYYLPg/CECROOGoSLJScnM2jQIAYNGkTPnj3Jzs7mgw8+4Oyzz+af//ynz/okIiISKhSERULE6NGjmT17NrNnz2bhwoXs3LmTP/3pTwCMGzeO1atXB7mHIiIidYuCsEiIiomJ4YknniAlJQWPx8OMGTOC3SUREZE6RUFYpJamTp2KMYbNmzcD0L59+zLzdUtPa8jMzOTvf/87xx9/PA0bNiQ+Pp4BAwbwwgsv4PF4yrRrjGHChAkATJgwoUyb119/fZX6FhERQZ8+fQDYtGlTrd+riIhIKIkI5MmMMY2BDkAcYCo7zlr7fcA6JVJLrVq1YtCgQcyfP5/8/Hz69+9PVFRUyf7GjRsDsHz5cs4991y2bt1KgwYN6NSpE/n5+fz888/MmzePL7/8knfeeQdvJW8GDRrEli1bSEtLIzk5mbZt25a02aVLlyr3Lzc3F4DY2FhfvF0REZGQEZAgbIw5HZgInFyFwy0BDuhSQ9YDBb67MSwoIhPA1O6DkfPPP5/zzz+flJQUNm/ezLvvvktKSkqZYw4cOMAll1zC1q1bufXWW3nooYdo1KgRACtWrGDEiBG89957PPvss9xyyy0AzJ49m/HjxzNhwgRGjx7N+PHjq923zMxMfv75ZwBSU1Nr9T5FRERCjd8DpzHmDOALIBLIBzYBOwHPUV4m9UFBFsw5Pdi9qJ1TvocGTf1+milTprB+/Xouu+wynnrqqTL7evTowZtvvklqair//ve/S4JwbRw4cIAlS5Zw1113sXfvXjp06MCIESNq3a6IiEgoCcTI6wScEPwGcLu1dk8AzilSp7z//vsA/O53v6twf69evUhJSWHDhg2kp6fXqLT2hAkTSuYUl3bOOefw/PPPl5muISIiIoEJwn2BLOB6a21RAM4nUucsXboUgPvvv7/SNX13794NwNatW2sUhEvPI969ezdr1qzBGMMJJ5xQZn6xiIiIOAIRhA8B6xSCxc2ys7MBWLBgwTGPzcvLq9E5jpxHvGLFCi6++GIef/xxEhISuO+++2rUroiISKgKRBCeB/QzxhhrrQ3A+SRQIhOcObb1WWRCQE4TFxdHVlYWa9eupVOnTgE5Z48ePXj77bc56aSTePDBBxk1ahTt27cPyLlFRETqg0DNEZ4F/AX4VwDOJ4FiwgJyo1l9UbzsWUV69OjBnDlzWLZsWbWC8NHarIp+/foxYsQI3nrrLSZMmMDUqVNr1Z6IiEgoCURBjW3AHcBDxpgPjDEXGmO6G2PaVvYIQJ9EfC4mJgaoeGrDsGHDAHj66aepzgcjR2uzqu69914A3njjjZKiHyIiIhKYILwJeBpoAFwMfAgsAzZW8tgQgD6J+FyHDh0AmDVrVrl9N910Ex06dOC7775j1KhRbN++vcz+/fv3884773DHHXdU2OacOXMoLCysUb969+7NueeeS2FhIY899liN2hAREQlFgQjCW7yPzaX+fLRHWgD6JOJzI0eOBODmm2/mhBNOYMiQIQwZMoRFixYRFxfHp59+Svv27Zk2bRpJSUn06NGDAQMG0LVrVxISEhg5ciRz5swp0+Y555xDkyZNmD17Nm3btuXUU09lyJAhTJw4sVp9u/vuuwFnPeOMjAzfvGEREZF6zu9zhK21Kf4+h0hdcM0117B3715eeukl1q5dy7JlywDIynKq73Xr1o3Fixfz7LPP8sEHH7By5Uo2bNhA69atGTx4MBdccAHDhw8v02ajRo348ssvuf/++5k3bx4//fQTHo+nXOW6YznzzDPp378/8+fP59///rdGhkVERACjhRzqFmNMEt5R8bS0tBqtJysiIiJS36Wnp5OcnFz8NNlam+7rcwRiaoSIiIiISJ0TiOXTADDGRAFXAecAXYB4YB+wBvg/4C1rbX6g+iMiIiIi7haQIGyM6Qu8C6QARy6M2hcYCYwzxoyw1v4aiD6JiIiIiLv5PQh757x+BTQBdgMvAMuBHUAroCfwO6AD8H/GmFRr7VZ/90tERERE3C0QI8J/xQnB7wPXWGvLVQYwxjwEvAYM9x7/pwD0S0RERERcLBA3y50PHACurygEA1hrDwI3eI+7IAB9EhERERGXC0QQPg5Yaa3df7SDvPtXAq0D0CcRERERcblABOF9OHOBq6IVzqhwjRhjLjXG/M8Ys8AYs90Yc8gYk2WMmWOMuc0Y06Ca7Y03xthjPLrVtL8SOlJSUjDGsGnTpmB3RURERKooEEF4AZBkjLnyaAcZY64CkoH5tTjXX4AbcW7AywMWA/uBgcCTwBxjTEIN2k0DfqzkkVuL/orU2pAhQzDGMH78+GB3RUREpF4JxM1yk3HWDn7FGHMi8Iy1dmPxTmNMe2As8EfAAk/X4lwvAvcBP1prC0qdYwDO8m39gIeBW6rZ7hRr7fha9EtERERE6hi/jwhbaz8FHgUigduBdcaYXGPMZmNMLrAOuA1oAEy01n5Wi3NNtdbOLB2CvdvnAnd4n15a0/ZFREREJHQEpMSytfavwMXATzijvtE40yCivc9/BC6y1v7dj91Y5f0a68dziIiIiEg9EZAgDGCt/cRaeyrQGOgDnOb92thae7p35NifBnq/1qRy3RnGmHeNMd8aY94zxtxtjEn0Zeck9GRkZDB58mTOPfdcUlJSiI6OpkmTJgwePJjXXnvtqK/Nzc1l0qRJDBgwgISEBGJjY+ncuTPXXHMNs2bNCtA7EBERCW0BKbFcmrX2AM5NbH5njAnHWY7tYmAizooUf61BU6cf8Xw4MN4Y80dr7dQa9i2pkl0K2CHixRdfZNy4ccTExHDcccdxwgknsHPnTr7//nu+//575syZw3//+99yr9uyZQvnnXceK1euBKBz587Ex8ezadMmXn/9ddLS0pg5c2aA342IiEjoCXgQDgRjzO3AE0dsngGMs9Yuq0ZT24F/Ah8AG3BWouiDc0Pe+cAUY8wea+3HNehmWg1eU6d4rIf9RfuC3Y1aiQuPJ8z454ORIUOG8O2333L66acTHh5esn3JkiVceeWVPPfcc1x55ZUMHjy4ZF9RURHDhg1j5cqV9O/fn1dffZXu3buX7F+0aBE//fSTX/orIiLiNj4NwsaYa71/zLbWfnjEtiqz1r5ay65sxZl3HAm0w1mf+AzgKmPM/dbaoir2438VbJ5jjPkNMB24DHjCGPOJtdbWss/1zv6ifdy+YnSwu1ErT/aYQqOIxn5p+9RTT61we69evZg8eTJDhw7ljTfeKBOE33//fRYsWEDLli354osvaNasWZnXpqamkpqa6pf+ioiIuI2vR4Sn4tz8thr48Iht1VGrIGytfRdnuTQAjDEnA/8D/gY0BW6uZfvWGHMvThDuCPSi+tM9kivZngj8UovuSR2yb98+3nrrLWbPns327dvJy8vDWkt+fj4AixeXvWw+/ND5ZzN69OhyIVhERER8y9dB+FWc0Lu9gm1BY62dZ4y5AGd6w43GmInW2s21bHONMSYTJ1h3oppB2FqbXtF2Y0xtuiV1yMKFC7nwwgvZtm1bpcdkZmaWeV48L3jAgAF+7ZuIiIj4OAhba6+vyrZgsNZuM8YsAk4GegO1CsJexesVh+Rca6m5oqIiRowYwbZt27jgggu455576NmzJwkJCYSHh7Nu3To6d+5MQUGZJa/JyckBICGhJgUQRUREpDrcFuAijvhaY8aY5kBL79MKR3dDXVx4PE/2mBLsbtRKXHi8X9r9+eefWbduHe3ateP9998nKiqqzP60tIrvlYyPd/qTlZXll36JiIjIYX4PwsaYb4El1trbq3DsE0Ava+1ZfuhHCs5IMPhm+bY7AANk49I5vWEmzG83mtV3mzZtAqBfv37lQjCUnxtcrGfPnixcuJC5c+dyySWX+LOLIiIirheIghpDgL5VPDbVe3y1GWP6GWMmGGM6VLDvPOBznOD/mbV2fal9txtjNhlj3jriNT2NMc8aY3oesT3aGPM34B7vpkettYdq0mcJXTExMQDs2LGj3L6CggKefPLJCl936aVOBfApU6aUmz8sIiIivhWwynJV1ACo0tJmFYgH7gfWG2O2G2N+McYsNsbsxQnB3XBGbq874nUJOEusHVnIIhJndYllxpidxpj5xpj5wB7gYZzv3Us4hTpEyhgwYAARERH8+OOPvPrq4UVQsrOzGTVqVIUBGZwg3L9/f3bu3MkFF1zA6tWry+xfvHhxhUU4REREpPrqTBA2xsQAXXCCZk0sBm4DPsKpINfN+8jDCcI3AKdYa3dXsb1NwDjva/cDXYETgEzgPeA8a+3v3Lh+sBxbYmIit9/uzAa67rrraNeuHf3796d169bMmDGDJ544st6LIzw8nOnTp9O1a1fmzZtHt27d6Nq1K/3796d58+akpqby9ttvB/KtiIiIhCyfzxE2xlwCHDm5sbMx5mh3VcUAJ+IsRfZeTc5rrd0LPO19VOd144HxFWzPAv5Rk76IADz22GMkJSXx3HPPsWHDBnJzcxk6dCh///vfadWqVaWva9u2LQsWLGDy5Mm8++67rFmzBo/HQ5s2bbjwwgsZM2ZMAN+FiIhI6DK+HtA0xjwAPFBqk8W5qawq1gIXlJ7D6zbGmCS85ZfT0tJISkoKco9EREREAi89PZ3k5JL6Y8mV1WCoDX+sGjEDZ1oBOAF4CrAGeKSS4y3O9IUNwK+aaiAiIiIigeDzIGytXUyp5cmMMeOBxdbaV3x9LhERERGRmvL7OsLW2hR/n0NEREREpLrqzKoRIiIiIiKBFNASy8aYATjV3ZrirNNbEWutfShwvRIRERERNwpIEDbGnI5TfKJc1bcjD8W5eU5BWERERET8yu9B2BjTA6coRSTwBjAYSAL+CSTjjBD3xlk54r/APn/3SUREREQkECPC9wLRwO+stS8bY34Akqy144oPMMacgzNifC4wMAB9EhERERGXC8TNckOAbKDS5dOstV8Cw4CewP0B6JOIiIiIuFwggnBLYJO11uN9XghgjIkpfZC19hdgNU4gFhERERHxq0AE4WwgvNTzTO/XdhUcewho4/ceiYiIiIjrBSIIbwFal3q+1Pv1otIHGWNSgK5ATgD6JCIiIiIuF4gg/B3QzBt0AabhLJH2sDHmH8aY3xhjRgNf4qws8VkA+iQiIiIiLheIIDwd2ACcCmCtXQ2Mw1mx4q/AR8ALQCdgI84qEyL1zvr167nqqqto0aIFsbGxpKam8txzzwGQkpKCMYZNmzaVHL9hwwYeffRRhgwZQnJyMlFRUbRo0YLzzjuPTz/9tMJzzJw5E2MMQ4YMoaioiEcffZTu3bsTExNDSkoK48ePp7CwEIC8vDzGjRtHp06diI6OpmPHjjz22GNYa8u1e/3112OMYerUqWzevJnf/va3tGrViri4OAYOHMhXX31VcuzSpUsZPnw4LVu2JDY2ltNPP525c+dW2N9ly5bxwAMPMHDgQFq3bk2DBg1o3bo1w4YNY86cOTX9VouIiPiGtTYoD+A04DngC+AD4E4gPlj9qSsPnDWWLWDT0tJsXVbk8djsvMJ6/SjyeHzyvVi8eLFNSEiwgI2JibH9+vWz7dq1s4C99dZbS/68cePGkteMGTPGAjYuLs526dLF9u/f37Zu3doW//1PnDix3Hm+++47C9jBgwfb4cOHW8B2797ddu3a1RpjLGBvuOEGm5eXZ08++WQbHh5ue/XqZVNSUkravf/++8u1e91115Xsa968uW3YsKHt16+fbd68uQVsRESE/eabb+wPP/xgGzZsaBMSEmy/fv1s48aNLWBjY2PtsmXLyrV71llnWcAmJCTY7t272759+5a0GR4ebt944w2ffP9FRCT0pKWllfzfhbP0rs9zl7EVjA5J8BhjkoA0gLS0NJKSkoLco8rlHCzi9g+2BLsbtfLkZW1pFB1+7AOPwuPxkJqaytKlSzn//PN54403aNKkCQDTp09n1KhReDweCgoK2LhxIykpKQB8/vnnNG3alJNOOgljTEl7P/zwAyNHjmTnzp2sXr2ajh07luybOXMmZ5xxBpGRkbRq1YqPP/6Y1NRUAGbNmsU555xDQUEBF110ERs3bmTGjBl06OAUdHzzzTcZNWoUUVFRbN++vaSP4IwIv/LKK0RGRjJ8+HCef/554uPj8Xg8jB07lmeffZbevXuTnZ3NsGHDeOSRR2jQoAH5+flcccUVfPzxx4wYMYK33367zPfmvffeo2vXrpxwwgkl26y1fPTRR/z2t78lLCyM9PR04uPja/V3ICIioSc9PZ3k5OTip8nW2nRfn8PvUyOMMd8aYz4zxjTw97lEguGrr75i6dKlNGvWjGnTppUJmMOHD+fee++loKCg3OvOP/98Tj755DIhGOC0007joYceoqioqFywLFZQUMDkyZNLQjDA4MGDGT58ONZaPv74Y1555ZWSEAxw9dVXM2DAAPLz85k5c2aF7TZv3pyXXnqpJJiGhYXxz3/+k+joaBYvXkyTJk2YNGkSDRo4/5yjoqL417/+BcAXX3xRrr3LL7+8TAgGMMZwySWXcPvtt5OTk8PHH39cYV9ERET8LRCV5QYCy621hwJwLpGAK54/O2zYMBo3blxu/w033MCECRMqfO2uXbt48803mTdvHjt37uTgwYMAZGdnA7B48eIKX9e0aVMuvfTScttTU1OZNm0affr0oU+fPuX29+nTh7lz57Jhw4YK273qqquIjY0ts61x48a0b9+elStXcsMNN5QL7l27diUmJoacnBz27NlDs2bNyuzfsmULb775Jr/++iu7d+/m0CHnR8HOnTtL3uPVV19dYX9ERET8KRBBeAtOiWWRkLR27VoAevXqVeH+du3a0ahRI3Jyyq4M+OWXXzJixIiS0FuRzMzMCreXni5RWosWLaq0f//+/dVud+XKlUfdv2XLFvbv318mCL/yyiv84Q9/KAn4FansPYqIiPhbIILwdOBuY0wXa+2aAJxPAiQuKownL2sb7G7USlxU7WcHHThwAOCo81zj4+PLBOGsrCyuvPJKsrOzufbaa/njH/9I165dadSoEWFhYXz99decffbZFU6pAMqN2hYrHq091v7K7g3wZbvr16/n97//PQUFBdx555389re/pWPHjsTFxWGM4cUXXyzZLyIiEgyBCML/AM4BPjTGXGmtrfizXql3woyp9Y1moaBhw4ZA5aOsAPv27Svz/PPPP2fv3r0MHDiQqVOnlptukJaW5vuOBtg777xDQUEBV155JZMmTSq3PxTeo4iI1G+BCMLPAGuBy4FfjTHLgZXAgUqOt9baMQHol4hPdOnSBYAlS5ZUuH/Lli3lpkUUryc8cODAciEYKp8bXJ8Uv8dTTjmlwv2h8B5FRKR+C0QQvh5n/bfi/+2P9z4qYwEFYak3zj77bCZNmsT777/PpEmTyk2RmDp1arnXxMTEALBjx45y+/bs2cNLL73kl74G0tHe46pVq7RahIiIBF0gKsvdAIz2fq3KY3QA+iTiM0OHDqVXr17s3r2bq6++mqysrJJ9M2bM4JFHHiEyMrLMa0477TTAmT7w9ddfl2zfvn07w4cPL6kOV5+deuqpADz77LMsWrSoZPuaNWu44oorSpZgExERCRa/jwhba1/x9zlEgiksLIzXXnuNwYMH88knn9CmTRt69OjB7t272bRpE2PHjuWjjz5i8+bNhIc7c6r79evH5ZdfznvvvcfZZ59Np06diIuLY9myZcTExDBx4kRuv/32IL+z2rn00ksZMGAAc+fOpX///nTp0oXw8HCWL19OYmIi9913H/fdd1+wuykiIi4WiIIabY0xLat4bEtjTP1ehkBcqVevXsyfP58rr7ySmJgYli1bRnx8PM888wxPP/10hStLvPHGG4wbN46UlBQ2b95MRkYGl19+Ob/88gu9e/cO1lvxmYiICP7v//6PsWPH0qpVK9atW0dWVhZjxoxhwYIFtGnTJthdFBERl/N7iWVjjAf4wVo7uArHfgecZq0NxNzlOqk+lViWqtmzZw/NmzcnISGBvXv3Brs7IiIi9UJIlFj2Kn9bvG+OLftCYy41xvzPGLPAGLPdGHPIGJNljJljjLmtpmWejTEDjTEfGmN2GWPyjDErjDHjjDEqFCLH9PLLLwOVr54gIiIiwRGoIFxVjYD8Wrz+L8CNQE8gD1gM7Mcp8/wkMMcYk1CdBo0xo4AfgIu9fVsJdAIeBL43xlRcYUBcZenSpTz//PNl1hK21vL6668zbtw4AP7wh0pjKVoAACAASURBVD8Eq3siIiJSgToRhI0xUcaYc4BewKZaNPUicAYQb63tYK090VqbhBOE04F+wMPV6FcK8BIQDtyNMyzfF+gMrAZOBB6rRX8lROzZs4ebbrqJhIQEOnbsyEknnUSLFi245pprOHjwIDfddBMXXXRRsLspIiIipfg8CBtjHjDGFBU/vJsHld525APIBT7HCZxv1fTc1tqp1tqZ1tqCI7bPBe7wPr20Gk3eBUQBX1pr/2W9E6qttZs5vMzbjcaYVjXts4SGHj16cPfdd3PCCSeQnZ3NwoULsdZy1lln8dZbb/Hcc88Fu4siIiJyBH/clGYoO8/Xcux5v3nABuBtYKIf+gSwyvu1SlMZjFPu6zLv03LVDay1c4wxq4BuwCXA877opNRPLVu25NFHH+XRRx8NdldERESkinwehK2144Hxxc+9q0bMttae7utzVdNA79dfq3h8W6C1988/VnLMjzhB+GSqGYS9q0NUJLE67YiIiIhIzQRimbIJwP+z997xcV1l/v/7TFEZ9WLJRXKvcYkdOwkJ6QnZJBBIgITedjdLWRbYZYEFlv3BLkvZ7JfelhIWNoQUSAGSQHp3HJfYjrslW5Zk9T4z0vTz++OOpDl3imZkaxLbz/v1mtfonHvuuUeW5Pnc537O87Tm4TpJKKWcWGL2zViRZj/w+SxPXxZ/DwIdacYcsY3NhbZpnCMIgiAIgiCcJPJRWe4rM30NO0qpTwHftnXfD3xJa70ny2mq4u9DOn2y5UHbWEEQBEEQBOEU4XQtXHEcy7bgBhYA9VjZJN6llPo3rXU008lxxnMEhzKMGU/1VjyNNTam6Z8NbJ3GfIIgCIIgCEIO5EUIK6WqsXL8XgssBkozDNcnWllOa30PcE/C9c8H/gf4AlANfDSLaQLx90xFOArj72PTWGPK6ijWHj1BEARBEARhppnxPMJKqUXAbuBzwNlAGZOZJVK9TvqatNZbgOuwIrh/p5RakMVp47aHSpVenVbZxgqCIAiCIAinCPkoqPFNYC7wCtamtTmAU2vtSPeaiUVorTuAnVjf89lZnHI4/l4YX38qFtvGCoIgCIIgCKcI+RDCV2LZDK7WWv9Ja92dYfPZTOOyvWeiFeiKf/36NGPG+7ecyKIEQRAEQRCE/JMPIVwEHNBa9+ThWmmJl0sejwTvmmp8XKzfF2/+TYr5LsTKIRwG/nBSFikIgiAIgiDkjXwI4f1A5UxfRCm1USn1FaXU4hTHrsEq4ewCHtJaNycc+5RSqkUplaq0861YWSOuVkp9ZtwrHPcY3xYf83OtdVeKcwVBEARBEITXMPkQwt8BFiqlrp7h65QB/wY0K6U6lVJblVK7lFKDWCJ4JVZasg/YzqvESrGWVNFNa30UuAWIAf8FtCmldmB5glcA24HPzND3IwiCIAiCIMwgMy6Etda3Y1V1u0sp9UmlVNkMXWoX8Eksm4IfS/iuxEpt9jDwIeBCrXVfLpNqrX8NXAz8CStf8FlYFeW+DFyktfafpPULgiAIgiAIeUTN9L41pdR4GeIGwBn/ug9LrKZCa62XzOiiXsMopRqIl19ua2ujoaHhVV6RIAiCIAhC/mlvb6excaL+WGO6GgwnQj4KaixM0Tcr/krFq5VRQhAEQRAEQTiDyIcQXpSHawiCIAiCIAhCTsy4ENZaH5vpawiCIAiCIAhCruQja4QgCIIgCIIgvObIhzUCAKVUAXATcCkwDyjSWl+ZcPwCrBRoj2uto/lalyAIgiAIgnBmkhchrJR6HXAXVuYIFe+2b4p7C1ZO3uuAv+RjXYIgCIIgCMKZy4xbI+KV3v4MNAL3YhW02Jti6O1YIvltM70mQRAEQRAEQciHR/hfgXLgi1rrm7TW/wcM2QdprfcAA8C5eViTIAiCIAiCcIaTDyH8BmAYq7rcVLRg2ScEQRAEQRAEYUbJhxCeBTTr7ErYRYHSGV6PIAiCIAiCIORFCA9hZYnIhiVA9wyuRRAEQRAEQRCA/Ajhl4A6pdTFmQYppW4AqoFn87AmQRAEQRAE4QwnH0L4h1jZIG5TSq1LNUApdQnwU6yUaj/Mw5oEQRAEQRCEM5wZF8Ja678A38OyPWxTSr0ILAdQSv1aKbUDeBKoBb6htX5xptckCIIgCIIgCHkpqKG1/pRSaj/wZeC8hEPvjb/3Af+mtf5JPtYjCIIgCIIgCHkrsay1/h+l1G3ABcBaoALwAfuAZ7XWwXytRRAEQRAEQRDyJoQBtNZh4Jn4SxAEQRAEQRBeNfKxWU4QBEEQBEEQXnPMuBBWSr1ZKXVEKfXpKcZ9Oj7uuplekyAIgiAIgiDkIyL8fmABcN8U4x4AFsbHC4IgCIIgCMKMkg8hvAHo0VofyTRIa92EVVVuUx7WJAiCIAiCIJzh5EMIzwVasxzbBsyZwbUIgiAIgiAIApAfIewHZmU5thaQNGqCIAiCIAjCjJMPIfwKsEApldHyED++ENiThzUJgiAIgiAIZzj5EMJ3AAr4jVJqcaoBSqlFwG8AHR8vCIIgCIIgCDNKPgpq3AZ8ALgQ2KOUuhfYAgwBlcDrgBuAYuAF4Gd5WJMgCIIgCIJwhjPjQlhrHVVKvQn4JfAW4N3AuxKGqPj7fcDfaK2jM70mQRAEQRAEQchLiWWt9RBwY9wH/BZgFVAOeIG9wP1a6x35WIsgCIIgCIIgQJ6E8Dha623AtpmYWymlgNdjCe2LgZWAB+gDNgM/0Fo/meOcXwb+vymGrdJaH8h5wYIgCIIgCMKrSl6F8AxzBfBY/OsY0ISVum0Z8FbgrUqpr2qtvzSNudtInwt5dBrzCYIgCIIgvOaJ6ij+qI9yV8WrvZQZIe9CWClVBZQy6Q1OQmudbQEOY2os8fst4E6t9WD8egXAl4HPA/+qlNqitf5TjnPfprX+8jTWJAiCIAiCcErii3j5cev/o2W0mbfNfg+X1bwBh3K+2ss6qeRFCCullmOJ0WuAqW4pNNNb10tYNoWIMZnWIeALSqn1wLXALUCuQlgQBEEQBOGMoXWshR+0fJO+cA8At3f8jOcHn+QDDR9hfvGiV3l1J48ZzyMcF6BbgXdgpUsLAu1YVoNUr7bpXEdrPWIXwTYejb8vn878giAIgiAIZwLHxo7wtaYvTIjgcY6ONfGto18lGAu8Sis7+eQjIvw1oAx4HPhHrfWrVTmuKP4+No1zL1dKrQZqgAGs6POvtdZd012MUqohzaHZ051TEARBEAThRHmo5z5COpjy2EhkiMf7Hua6uhvzvKqZIR9C+ELAB9ygtfbn4XpJxDNK3BRvPj+NKS6xtd8GfFkp9TGt9f9Oc1nTinwLgiAIgiDMFL2hbrYNv5hxzMO993NZzdV4nCV5WtXMkY8Syw7g4KslguPcAmwAQsB3cjivEyuifS5WNNiDlaLtYaxKeLcppa4/uUsVBEEQBEF4dXi0909oYhPtQkcRn170b8YYf9THI72nx3arfAjhncCcPFwnJUqpc4Dvxpv/qrVuzvZcrfX/aK2/qLXeprUe0FqPaa1fAN6IVQlPAd+OR5xzpTHN69xpzCUIgiAIgnBC+CM+nh18wui7tPoqVpedzcby843+pwYeIXYaFAPOhxD+OjBHKfW+PFzLQCm1CCtDRBFwB/DfJ2NerbUG/iXeXAKsm8Yc7alewLR9x4IgCIIgCNNl58hWYyOcAwdX1b4RgLfUv8MYOxIZ4pB/f17XNxPMuBDWWj8MfAz4kVLq20qpNUqp4pm+rlJqNlamiDnAg8AH4wL2pKC1PoS1cQ5g6cmaVxAEQRAEYUp0DPb+Ezx/yeTrlY/DCURpd3t3GO115RupLagDoKF4AY1FC43jW4c3T/tarxXykT4tCvwIy1/7CWAX4FNKRdO8MqVAy/aa1VgieAnwNHCT1jp8ovOmYHzO06lCnyAIgiAIr3Xa/hd6H4HwwOSr/yk4+oNpTRfVUfZ4dxp968s2Ge1NFRcY7R3DL57y9oh8WCNUjq8TWpNSqhR4CFiDlb/4eq31dFKmTXWdWqAu3mw/2fMLgiAIgiCkZPhlOPLd1MdafwZ9T6Q+loHm0YOMxUaNvrXl5xjtcysvNJcRGeKw/0DO13otkQ9rhCPX13SvpZQqBB4Azgf2Atdorb0n63ux8U9Ywn0YS3ALgiAIgiDMLNEg7PsMkCESu/8LEB7KadrdI6YtYn7RIqrc1Ubf7MK5NBTNN/q2neL2iHxEhPOCUsoJ3AlcATQDb9BaD2Q+C5RSn1JKtSil7rT1r1ZK/SheSCOxv0gp9QXgc/Gub8bLOAuCIAiCIMwsfY9DcIp99VEfdN6X07R2f/Da8g0px22qMKPC24ZfJKZjKceeCpw2Qhi4Gbgh/nUMuEcp9VyK1z228yqBBSRXdHMDHwX2KKV6lFLblFLbgH7gP7H+7X4BfGOmviFBEARBEASDrnvNdtk6uHQX1Fxu9nfcZW2oy4KRyDDtgWNG37qyjSnHnpvgEy50FLGi5CzGoqMpx54K5HWTl1LqMuBqYDlW2WUvcAj4i9b66ROcvjDh62XxVyqOpem30wJ8Casy3kpgBVAA9GB5kH+utf7LtFYqCIIgCIKQK2PHYXCL2dfwblBOmP8h6H9ysj/QDgPPQ83FU07b7D9otAsdRSz2pJZRc4oauKb2zSwpWcnasvUUOApTjjtVyIsQVkotxMrjO56NObEAhQY+p5TaDLxXa90ynWvESx3/7zTO+zLw5RT9Q8BXp7MWQRAEQRAEAF8wyoHuAI1VBcwqdeGYVg2uOF33Y8mmOM4yqL3K+rp8A5QsB/+hyePHf5uVEG4aNYXwYs8ynMqZdvzNcz+Qy6pf08y4EFZKVQFPYtkPQsDvsTaydQP1wGrgbViR1yeUUhu11oMzvS5BEARBEISZZDQU41tPddEyYG0lKnQpltUW8fcX11HoytGdqjV0/9Hsq38jOIusr5WCee+CQ1+ZPD7wvLVpzl2Zcerm0UNGe4lneW5rO4XJh0f4c1gi+Dlgsdb6PVrrr2mtfxF/fw+wOH58AfDZPKxJEARBEARhxgiEY3zn6UkRDBCMaPZ0jfGb7f25Tzh2zLI7JDL7BrNd90ZwFCV0RKH/mYzTRnSEltEmo2+JZ0Xu6ztFyYcQfgsQBN6ute5INSDefxNWgYob87AmQRAEQRCEGeMnL/TQ1BdMeey5Iz42H/XlNuHAC2a7YBaUrTb7XB6oMote0P9UxmnbxloI2ZJfSUT45LIA2KO17sk0SGvdDewB5mcaJwiCIAiC8FqmuS/A7o7Mtbx+va2PAX8OxXQHbUK46kLLDmGn9jKzPfAcxNJnebXbIuoL5lLqKst+Xac4+RDCQawUZdlQHh8vCIIgCIJwSvKXA8NGu6TAwaVLTHEZjGieODyS3YSxMAy9ZPZVX5h6bM2lGDkJoqMwlL7uV5OtMtzSkjMnGgz5EcK7gcVKqSsyDYofXwrsysOaBEEQBEEQTjo9vjDb2828utedVcEHzqvlgoUlRv8zzV5CkSxy/Y7stgRtIlWvSz22oBbK15l9GUouHx07c/3BkB8h/DOsW5N7lVL/oJQqTjyolPIopT6BlU1Cx8cLgiAIgiCccjx2cASdkOGs0KUmosFvPMt8QO4LxdjS6p96UrstonQlFNSkH28vrmH3F8cZjfrpDXUbfYs8S6dez2nEjAthrfXtwG+xbA/fAfqUUvuVUk8rpfYDvcC3gQrgDq31b2Z6TYIgCIIgCOMEYwEO+w+wZeg5Hu65n82DmTMtpCMW02w5Zm6Cu3RJGZ4CKyfv3IoCzqovMo4/fmgEnaicU2G3NlSlsUWMU32R2Q60QSA5X0HbWIvRdioX8wobM899mpGXghpa6/fEC2Z8BmjEqtKWGHtvBW7VWv8wH+sRBEEQBEEAiOooP2j5L/b6TGfmYLif6+pyS2TV1B/EGzStDpcvKzfaVy4vZ193YKLdOhiiYyTMvIqC1JPGwuDda/ZVbsq8kNLl4KqASIJXefAlmGOmWzs2dsRozytsxOVwZ577NCMf1ggAtNY/0FovwCqgcSPwvvj7aq31QhHBgiAIgiDkm3s6f50kggF+13U793fdNXW0NoGdNm/wvAo39WWmsDx7rocqj1m1bXtbBnuE7wDEbHkEys/OvBDlgMrzzL6hLUnDjo0dNdrzixdlnvc0JG9CeByt9X6t9QNa69/E3/fnew2CIAiCIAjPDzzJI31/Snv8Dz138/LIS2mP23n5uCmENzR4ksY4HIqNDeamuR1to0njJhjZbbaLF4G7YurFVJ1vtgdfApuobw2IEJ4RIayUqldKrVNKLcly/JL4+PqZWI8gCIIgCEIioViQuzp/NeW4+7rvJKajU47rHAnR7Q0bfRvmlaQcu7HRFMitQyF6fOGUYxmxRasrpogGj2OPCIe6rep0481YkE5bpToRwicBpZQLeAbYAazJ8rQ1wMvA40qpvEepBUEQBEE4s3hh8Gl8Ua/R96GGv+edcz5k9B0PtPLSUOqsC4nYbRGVxU4WVKf2/S6rLaKs0JQ7aaPCdiFsT42WDs8iq/pcIoOT9ojjgTZiTPqZFYrGooXZzX0aMROi823AMuAerfUD2ZwQH3cXsAopsSwIgiAIwgwS0zEetVki1pSu5+LqK7iq9joaiswit/d330l0iqjwvm6zktz6eR4cqSq/YdkjNtjtEe0pfMLBPggcN/um8gePo1RyVHh428SX9o1ydQVzKHYWc6YxE0L4Jqx8wLfmeN43sfIN33zSVyQIgiAIwumBjsLQDohk8NVOwR7vTjqDpsC8etb1ADiUgxvr320c6wl18Yr35bTzRWKapj5zQ9vq2ZlF5Tk2/3Bzf5DRkE1s26PBTg+U5JDnt3Kj2R7aMeETbrVtlFtwBtoiYGaE8EagT2u9I5eTtNa7gB7g3BlYkyAIgiAIpzrBPtj1Ydj5fnjuAtj+bjj2P0mbwKbi6YFHjPa8wkZWl05GWteXb2JBsbnN6fnBJ9POd2wgSDBirmF5XVGa0RYr64pwOycjxlpjpFUDwPuK2S5bC8rMOJGRCpsQDnVP5BNuD7QahxqLF2Y/72nETAjheqy8wNOhDZh9EtciCIIgCMLpwOBLsO1tMPRivCMK3t1w9PvQ8qOspwlEx9jjNSOtV9Zeh0qwMSiluLT6KmPMzpFt+CKmp3icgz2mgJ1X4aasMLNgLXA5WD7LFMt7O017Bb4DZrss261XcTyLwWVWs2N4O1prjtuE8DybHeRMYSaEcBSYbjZmV/x8QRAEQRAEi/AQ7P1HCPenPn7sx9D9YFZT7fHuJKxDE20nTs6tSK7Udl7l63GpSTkT1RG2DD2Xcs4DNiE8VTR4HLt9Yk/X2GTeYq3Bu888oWxlVvNOoBRUnGP2De9gINzHWMy0lth90WcKMyGEu4GFSuUSu5/INrE4fr4gCIIgCILF0e+bVdJSceBLMDr1A+kdI2ZhiRWlqylxlSaN8zhLOKfc3GyWyh4RjWkO95pCeEWWQnjNHFMI9/sjdHsjViPUC+EB84TSVVnNa1BpF8Lbk2wRxQ4PNW5bhokzhJkQwi8AZcC1OZ53bfy8qXOUCIIgCIJwZuA9AB33mH1la6DujWafDll+4QxEYmF2jWw3+s4pPz/NaHh99eVGu2Wsmb5Qj9HXOhhK8gevmJWdEJ5X4aai2Iwb7u2K2yN8tnpjTg8UTyNqa/cJjx7luN+0XMwrajSsIWcSMyGEf4uV/eFWpVRVNifEx/03VraJO2dgTYIgCIIgnIoc/Q4k5LvFUQSrvw1nfRPmvM0c2/1HGD1GOvb79yRZAjZUnJdmNJxVuo5SZ5nR9/LIVqPd1GdGg2eXuakodqWdMxGlFGts9ogJIey1+YNLV1qlk3OldCU4zGu02zbhNRQtyH3e04STLoS11g8DTwErgKeVUhmzQCilzsMqwLEUeEZr/dDJXpMgCIIgCPlnJBDl2WYvrYNBorHcMjsAMNYOAzZf7vxboGiO9fWiT9hEXgyO/STtdK+MmAmtlniWU+WuTjveqZycXW5GVHfahPDRfjNt2pLawrTzpeKselOkHuwJEIvp5IhwaY7+4HEc7qTcw+221HFn6kY5sDanzQQ3A89jVYx7USn1MvA0cBTwAyXAIuBSYANWBLkZySEsCIIgCKc8Wmu2HPNzx/Z+fCErmlvoUpzbWMIHzqvF6cjyMXzX/WbbVQ6N759sF9TAvHdB222Tfd0PwuJ/hMK6pOn2+8xI6LqyjUlj7KwvP4/nB5+aaB/07cUf8U34io8OmEJ4cU1uQnhlvWmjGAvHODYYYlGSEJ6GP3icynMmsm1EgM7YqKW84pypG+VghoSw1rovHgn+EfBO4BwswWtHYT3vuBP4e6310EysRxAEQRCE/KC15ucv9rK5xayUFoxonjvqIxLT/O3rZuGYSgzrKHTeZ/bVvwns1c8aPwTHfwux8dRjMej6Ayz4W2PYcHiQ48E2o29V6dopv581ZWfjVgUTmSZixNjt3c4FVZfiD0UnN7fFWVidmxCu8riYXeamyxue6DvQOcgie0W5ExHCCT7hbmcBUds/vVgjZgCt9YjW+r1YFokvAg8Bu4Aj8feH4v0rtdbvEREsCIIgCKc+zx/1JYngRF485uf27WnSoCUy8LxVACKROW9NHldQBXW2/fld9yUV2djv22O0ixzFLPJMXaWt0FHE6tJ1Rt+4T7ilP2T0uxzQWFkw5Zx2Vtmiwvs7bRkylBtKzAIfOVG+DpQV+2x3mkK9yl2dMmvGmcKMCeFxtNbNWuuva62v11qfo7VeFn+/Pt7fNNNrEARBEARh5hkcjfDbHQNTjnuqycu+rrHMg+x5gctWp/fJzrnRbI8dg2GzJLLdFrGi5CycWWZ6XV9hbnfa59tNVEeTbBGNVYW4nLlnX1hl8wkfHnAS0QkP7UuWWV7f6eIshtKzgGQhPO8MjgZDHoSwIAiCIAhnBrdv72csHDP63rOxho9dVIfTpjju3NFvbQpLRSwM/c+YfbNvSH/h8vVQvMjs67rXaNqF8CpblDcTa8tMd+do1E/LaHPSRrnFOdoixllhiwgHY06OhpdNdkx3o1wi8XzCx21C+Ez2B8NpJISVxUVKqVuVUi8qpYaUUiGlVIdS6vdKqcunniXt3BcopR5QSvUqpcaUUvuUUl9SSmWXKFAQBEEQXuNorQnFgoxEhukL9UxWOMuS48MhXm43U5NdsLCEK5eXs6mxhHefU2Mcax8O88yR1CWLGdoGUdux2ivTX1yp5Khw72MQs6wLvaFu+sJm/t+zsvAHj1PlrmFeYaPRt8e3MykivLAmd1sEQFmhkwabpWJ/MEGol52AP3icuE+43WleR4Tw6cMVwLPAPwPnYlWo24NVpOOtwBNKqf/IdVKl1Hvi874ZCAL7sVK9/TvwjFLKc1JWLwiCIAivEgOhPv696XN8ZM+7+dS+v+azBz7K91q+znA4++07jx40fa2lBQ7elSB+L11SliT2HtgzRCSaQnD3PWG2y9akzAJhUP8mjFQIUR8MbgbgkM8sVVzmLM85ZdjqsvVGe1f/IYbGokZfrhkjErH7hA+EEoTwiWyUG6diA2Mo+mxC+ExOnQanlxBWQBPwMaBWa71Ca30OUAN8PT7mX5VSb8p6QqUWAr8AnMBngcb4nMuAg1iC+79O1jcgCIIgCPnmkH8fXzn8GY6NNRv9u7zb+dKhf2SP9+U0Z04yEojywlFzg9xly8opLZz04DocindtMHP2Do9F2d5u21inNfTbhHBtFg91C+ugwpagqvdRAA6PmsUplpWsyrmS2hqbED42YGaLKHYr6sum7+O1+4SbQisJ6QJAQenyac87gbuSDtuGOweKuYUNJz73KczpJIRfAlZprX+stR4c79Rah7TWXwAejnfdksOcnwEKgUe01rfq+HMirfUx4K/jY/5OKVV/4ssXBEEQhPzSEWjnv4/8O97oSMrjvugIPzz233QHOzPO81TTCJEEv6/LAVcuK0sat2p2MctnmVHTxw7aru3bB0FbtohMtohEZr3BbPc9AbEwzf6DRvfSktw9t8tLVuFWk9FUPWYKyIXVhThOoEzx8llFJJ4eoYDm0ErwLLTKK58E2kvMjXH1OHA7pmfnOF04bYRwPF1bJMOQR+PvWd1WKetWcdxw9IsU13sBOAC4gbfksFRBEARBeNWJ6Si/bP8hER3OOC4YC/CT1m8RjqUep7XmuSM+o+/8BaVpywxftbzCaDf3BzmSuOls4AXzhKJG8GSZOqz2KrMdGcE/8ExS/uBlntyFcIGjkOUlCRaFMdMznGv+YDueAgcLq2w+4dC6iWwPJ4N2d7nRnhf2gY6lGX1mcNoI4SwYN99Mka9lgvlAvIYjz6cZM95//nQX9aoxvAN6/gxdD0DHPeA7OPU5giAIwmnDY30P0Tx6yOhbWbKGf1nyVVaXmiV5j40d4f7uO1POc6Q/SJ/fjENdsaw85ViADQ0eqjxm2rInDiVEhYe2mCfUXAzZRlqL5iSVE27uMdOwuVUBC+wZJrJkwiesVZIQPhF/8DgrbfaI/cF1JydjRJzjmDczDWE/+M/sLLZnhBCOR3dvijfTiVo743lLgkBHmjFHbGNzWVNDqhcwO9e5pkXLj2DfP8OBL8Khr8Dgi3m5rCAIgnASiIWg+0+w472w9UY49B/Q/ZBVjS0LfBEv99mEbX3BHD616AssL1nFxxd+ljmF84zjj/Y9yEAouRDGlmOmx7e+zM3C6vSP250OlSSUt7f7CYRjEA0m5f+l8rxsvqVJbPaIw6NmqeJFniW4ppmTd8InHKqBmClaF51gRBhgVZ25rqPh5YwVn5yIsNaa9lCX0dcQDVqBsTOYM0IIY/mCNwAh4DtZnlMVfx/S6XPIDNrG5kJbmtfWacyVOw7bH2wslHqcIAiC8NrCux82vwH2/wuM7AT/Yei4C/Z/FvZ8AiKjU07xeP/DBGMBo+9DjR+jIP7ZSauKnAAAIABJREFUUOgo4iPz/wmnmrQ3RHSYP/TcbZwTjWleajWF8OsWlEy5Ee2iRaUkVlgORjQ72v0wsgtiiSnJFFSem3R+RmrMjXXNtqUs8azIbb4E5hU2UumqhjEz00JFkTMpyj0dlnk6cSZEbWM4OTy2+ITnBRiODOGzpaSzhPD2kzL/qUpehLBSyq2U+jul1H1KqV1KqWal1JE0r+apZ8zp2ucA3403/1Vrne3841aKTApx/K+1OMOY1ybKdrceC6YeJwiCILx28DfD7r+DcJoSxf1Pw873Q7Av7RTBWIDH+ky7wOsqL2F5iRl5bCxeyCXVpuf2uYEn6A5OPiQ90B1gJGBGoc9fMHW53opiF2vmmB+dLxz1wdBL5sDSVeA2PcVT4lkAHks8RoAjLjMt2bJpbJQbRynF6rKzk2wRi2oKc85CkYrCwH6WFJhWxQP9Jy6wAdoDLea1dIzaWBiGdiSVoz6TmHEhrJSqBbYDP8baVLYWWAQszPA6WddeBPwJS9TeAfx3DqeP3ypn2k45HlbN1necSGOaV463vtPEvks0zSYIQRAE4TVCsBd23QLhwczjfAfgwL+k3QT1zMBj+KPm5rbrZqWu2nZ93dspSAicxIjxUM/9E2176rOF1QXMLs/OdnDhQlMw7+8OMNBj2hiomuYWnHhUuM1ZREiZUmfpNDbKJbKmbH2SEJ5flXpjYM74DrCywKyAt797OhIjmbaxY0a7IRq0RGCoGwLHT8o1TkXyERH+BrAGOA58ArgAWIIlhlO9TsozAKXUbKxMEXOAB4EPZrA4pGL8f5pKlf42r8o2Nmu01u2pXkDXlCefDJKsERIRFgRBmGkiUU1zX4ChsUxJjtLQfCuEzOpoFDVA3RvBWWL2D74I7f+XNEVMx3is7yGj7+yyjTQUL0gaC1DpruKK2muNvheHnsUbGUFrze4O04axqdG2jgysn+eh2D358aqBLX2zbAuYphCuvQyAwy4z6jyncB6lruS0brmwwrMWAqZ/uqCkO83oHPEdYFXBbqOrdTCEL5id9zsTbbaIcGMk4XP/DLZHnKRbmIy8CQgDV2it87I1USlVjSWClwBPAzdpPUV+mGQOx98LgblYQt7OYtvYUwd7RFiLR1gQBGGmaOoLcO+uQY70BwlFNQpYUlvI6xeVcsmSsqkfqw9shh5TwFK2Ds7+GbhKYPQY7PyQKZSPfAeqXw8lSye6Dvn30RsyRdu1damjweNcXXs9j/Y9SDSeoTSsQzwz8BhnF7yJgVFToK2bm32+2wKXg02NJTybkHpt29gFXFvyO6uhXFBxTtbzGZSvA3c1TW7TFnGi0WCAYV8RaDPqPeTaDZzgpjatwXeQxQWjuAkSjj901sDBngAbc7jJSEVbIDkiPMHwDph9ZmaCzUdEuAI4mEcRXAo8hBWF3gpcr7WeznOFViajs69PM2a8f0ua469dZLOcIAjCjKO15snDI3zz8U4O9AQIxcsJa6CpL8ivtvZz25a+1GWGx4mF4fB/mn2uclj7fUsEg+WLXfV1jBLDOgzN3zJOe3bArNg2r7CRZZ7M5Xsr3VWcW3Gh0fdE/5/Zedy0V9R4nMyryC0bw3nzTXvE0fBy+iLxUspla8E1zUISyomuvoQmmz94acn0N8qNc3TA9gS1oI/DwZMQUQ0ch8gIbhVhWYFZEvpAdyDNSdkRjoXpstkfGhOF8NCZmzkiH0K4icw+25OGUqoQeAArr+9e4BqttTfzWamJ2yjuizf/JsW1LgRWYkW7/zCtBb+a2FPHiBAWBEE46fxx7xD/t62faIaaBc8f9fGtp7oIpxPD3X+CsRazb/GnoKDG7Ks6Hxo/YPYNPANDVjKi0aif7cObjcOvr74iq01eV9VeZ7QHw/1sbus1+tbN9eS8YWxFfRElBaYU2RaIi+7p+oPj9FdtYsj2Wbe0YF6a0dljFP8AKG7l2NgRRiLDJzaxb9IfvarQtEcc6Dkxn3BnsJ0oZvS+IdEaMXYUUqTGOxPIhxD+ObBMKbVxJi+ilHICdwJXAM3AG7TWA1mc9ymlVItSKlWm8FuxskZcrZT6zLhXWCm1ALgtPubnWuv8+HpPJuIRFgRBmJKYjtET7GKfdzdj0dzEyO6OUR54ZSirsQd6Atz1cgohoqPQ+nOzr2wNzHl76okWfBTc1WZf87dAa14aep5Qgg3OiZMLKi/Jan2LPctY4kkozBoppnPQdFfmYosYx+VQnNNgnrctEH/Ymmv+YBuH3aaVoDQWYbbvxB9Ot9iFcJFVtW6fd3eK0TngOzDxpX3D3PHhMMNj0/cJt9lupGrddRQ7zWj5mZpPOB8e4e9hZUK4Xyn1ca31AzN0nZuBcaNTDLgnzZ1pp9b6poR2JbAAaLEP1FofVUrdAvwS+C/gk0qpHizbhRsrG8ZnTtY3kFeSskZIRFgQBGEcf8TH3Z2/5qXh5yfy7TqVi5Ulq7mu7kZWla7NeH6/P8LPNvdij/Fetbyca1dV8ErnGHds75+wSgA8cdjL0toiXpeYTaH3UbDt9mfhx0GliWO5SmDhx+DwVyf7vK/AwHNsGXrOGLqufCMV7sqM30cil1RfNVmJzr+CxFia26lYWV+U+sQpOHe+6RM+El5Jf2weNRXrpzXfOE2Bo0Z7aWQM1f8UzH7TtOcMRmIcH7FtOfJYQnivbyevq7p42nMnCuGF7sMUOcIEYpMR7QM9Y1mlpkuF3R/cWLwAyiMwmPCEYHhHUjGSM4F8COHH4+91wL1KqUGsiK0/zXittb5yGtdJDHEuI321t2Np+tMt5tdKqSbg88CFWG74I8BvgW9qrU/MuPNqIUJYEAQhJa94X+a2th8wHDGjuVEdYa9vF/t8u7m+7u28uf4mHCp1jtc7dvTjD5l+iJvWV3HtKkt4XrKkjMbKAr7xeKdhifj11j5W1RdRUeyyNk+1/sKcuHSVtQEuE3PeZmWMSBDQw60/45DTzPBwYdWlmeexsaniQn5z/BeEdBC8pq94VX0Rha7pPWReWV9MiTOIPzr5Mb6dt3G1/XMqRw77DxjtpZEADDxnea6nWVnu2EDIlnI3CkWW93aPdxda6+nnE/ZOrtepYiyvCrC7f3Kd+7sDJyCEW4x2Y9FCoNIUwkNnZuaIfFgjLou/3Fgu/mqsCPFlGV45o7X+X621yuK10Hbel+P9aa+rtX5Ba3291rpGa12ktV6ptf7KKSuCQQpqCIIgpGCP92W+1/KNJBGciEbzh557+Gnrd4mlyNW7t2uMl9tN0bmxwcM1K83CEItqCnnfJtPnG4ho7t4Zz8jp3WP4RgGYfwtMJbQcbphvbm3ZHjiCTohPFzqKWFu2IfM8NoqdxWyseB1oBT4z+8J0bBHjuByKDZ6dRt/W0RNLqT8a9XM80Gr0LY2MQdQ34ZmeDkn+4KIucFgR4uHIIMeDbdObODRg5fNNYNVs09pxYJr5hLXWyTmEixckZ+TwHYBIuhjl6Us+IsKXTz1EyDuSNUIQhNONaBA677Z2wPsOABrKz7ZK9M5+y5RRwGb/IX7QcutEmrCpeGn4eaq7arl5zvsn+iIxzW+3m17f0kIHHzy/NmWk8KLFZezrGuPFY5MCZHOLj0uXlrF84F5zcNE8mJXlA9P6N8HRH0ykU9taYObOXV++aaKcci5cVH05mztaIGqKtHVzT6DAarCPTa6HeY7JzXHN3koGRiNUe6YnU46MHjKEv0vHWBSJx676n4LqC1OfOAUttowRhZ4+Env2eHfSUGSWX84K+w2Po4hVDXNh7+QWpB5fhH5/hJqS3P5NrNLKI0ZfY9FCcFVZKeomft9jVsnuqZ44nGbMeERYa/10rq+ZXpOAWCMEQTi96HsStr4Fmr4JfY9CoA0C7dDzIBz6Mmy7CYa2pT3dH/Hxo9Zbrcf+Cawt28B/LP8OP1z9f9xQ/w4cto/NP/c+wLMDj0+0Nx/10WHzkL51XRUlBenL5L5jQ41RWALgzu296O6HzYGzb4Q0VowkHAXQ8F4AhpWTQ7bCEvZ0aNmyomQ1xaPm3vcSj4/akulZDQAY2sJZhbsoVmY6tu1t049O2m0RCyNB3OPCuO/JaZcUtkeEF1SbonSv14xsZ43PXC8ly2moSs6oMZ2ocFJpZUcRswrqwVkEZavNwWfghrl8WCOE1yJSUEMQhNcYvmCUSGwaAqX1NtjzD5bwTcdok1VwouOelIdv7/gZg2Ez0dD5lRfxyYVfYF5RI8VOD2+uv5lPLfoCTkwxekfHbfSGuonGNH/aZ1oq5lcVcMnizJXMKoqd3LC2yuhrGYywczRRpKjcCx7MvQmcJWwrKEMnRKMLHUVWmeBp4FAO3P6zjb5QyU7CsVxrViUwuAWXirChyEzJv+0EhHDT6EGjvTSSICCDneA/SK54g1H6/ObTgvWz6432If9+QtOxGtqFcNlKHEqxss7cgLh/GvmEk2wRRfNxjG+2rLAl9DoDfcIihM9UJH2aIAivAQ71BvjRcz18+oFWPnFvKx++u4V/fqCVX2/to9ubhbg6+gM48q2pxwGg4dBXoO1XRu/WoReSMiqsLFnD3zR8fFIwxFlTtoEPNHzE6AvGAvyi7QdsbvHS6zOF0tvWVeFwTL156opl5dSXmVHV+7zvIabj51a/HormTDmPgasM6q9nm80WsaFs47RsEQCDoxFGfKYtIlyym93eExBQQy8BsKnoeaO7qTfI4GjupagjOkKz/5DRt8xh8zD3PZnzvEdt0eACp+LC+pWoBCkV1qGkaHRWeG3nlFoe7JX1ZiR/f88YOsdotn2jXENRQjltu0/Y+8oZ94T4pHqElVLjuXU7tdZftPVli9ZaJxWwEE4ySdaIE7ibFwRByJHB0Qi/2d7PDtumMq1hYDTKU01enm7ycunSMt51TjVuZ4q4Tefv4dhPkvvLz4b6660cvJ2/B5soovlWKKyDumsJxgLc2flL47DHWcIt8z+BK42n+KLqK2gaPcgzA49N9B3y7aerpZPE+lGLagpZMyc776zToXjLmkp+unmySEV7ZBHbAq/nvOLnYPZbs5rHzlD9dRwaMx93b2J6ac7Ayo1s4BgFTwvPDTxhbaTLlbH2iUj+6sKXKVZ+xrQltDWwvd3PVcsrMkyQTNtYS5LFZUnF68CfEBntexIWfjSnee3+4PlVBZQXlLHYs3QyrRyWPWJ12dn209MTGU0umFJqZeVYZUtJNzgapccXSbppykRS6rSihZONig1YeQzi4joWtDZpTre09SnIyd4s98H4+wHgi7a+bNGkqOQmnGSShLBEhAVByI5gLMAB3172+16hN9SNW7kpc5Wzrnwjq0rWpBWQ4xwbCPLdZ7oZmqJAgAaeavJypC/I319cx6zShHmHd8GhryaftPjT0PjBycwK894JR38IrT81xx34EngW8aBvR5Il4r1zb6HKbavaZuMdcz7IPu9u+sLWZjS8ZzEyav6/+ubVlTml0jpvQQkP7hvi+PBkYOJB302cW7YXVTu9fefbw92mLULHWDOwGRZ8bFrz7e6weVRLD4KKsce7k5HIMOWu3ETreDQYwK0irC/eyebRyc1a21pHcxbCTaNmdHV24VzKZ10NHXdNdvr2QaALimZnPa/dH7yoxoqqry5dbwjh3d4d3Iytwl8m/AfByDjthBIrA+yccjcVRU6GA5N/K/u7x7IWwqlKKzcUJ0SE3RVQshT8hyf7BreIED4BPhR/H07RJ7yWkM1ygiDkSDgW5umBR/lj9z14bbvQAR7vf5gSZyk31L+Ty2quxpliY9eB7jG+83S3UUhiKlqHQnzj8U7+5co5lhgOj8C+fwJte5K1/N9g7s1mn3LC4k9YhSaOfHuyPxagd+8/8mePKbJWl57N+ZUXTbmmYmcxH2r8GLce+bLV0W9WaJtfWZBzJgWHUrxldRk/emFSmLdFFrO3+BbWTDPv7TZbSeUNIR8F/sPgOwilK3KaKxzV7LNv1iqzsh1EibJl6DneUPvG3BY4aPqCN80aYnNCAPNwb4ChsQiVxdnLlaT8wZ6VlrBzlUMk4fe2/2mY946s5tRaJ1WUW1RtCeE1Zev5Q8/dE/0dwXZ6gl3UFWYpsu3+YM8iayMboJRVqGRLQlaRA90BLltantXUKUsr27NaVF1gE8Kbc46Wn8qcVCGstf5VNn3CawC7P0yHQcfSVysSBOHUJjpmiY7BFyDUD9FRyx5QugqqL4bieRlP7w528MNjt9Juy81qxx/18ZuOn/PMwGN8bMGnqS+cO3Hs2ECQ7z2TLIJLCxxcu6qC1XOKicRgW6ufJw6PGOMGR6P81xNdfP7KOVS3/RcEzZyrNHwgWQQnMv9vrMfwnZOb5R5whIkkiGkHDt4590NZR3FXla5lU8UFbOs6BqNLjGNXryyfVmGFczzbqXe66Y42TPQ9NHARa3KeCYbDgxzy7zP6NoW81hcdd1k3DjlwsGeMYCTxZ6etiHCczYNP5yaEtTYiwgBrGusoPK4mrqOBHe2jXLEsO+GntU4SwstKVlqp86ovtrKIjNP/VNZCuH80wkjQzBk9HhFe7FlKmbPcuDnc5d3GGwqzrGDn3Wu2y8wczavqiw0hvL97LOvCHalKK3ucpsebqgug/deT7eFdEPGBa3rFO041RPWcqdgLaoD4hAXhNUxMa8bCsZw3yhAegaPfhxcuhT0fh+N3QO9fYOBZyz97+Kuw5Vp45RPJH8hxdo9s598Pf25KEZxIW6CFrzZ9nr3eXYBVcvjbT3URiJjrX1FXxNfe1MC1Z1Uyv6qQxTWF3Lyhmi/91VzmlJtR0H5/hO892USw05ZWrPJcWPyPUy9q2ect/zDQ4Shgc4Epri6vuYZ5RY1Zf48AN895P46By4w+t3uM8+ZPT0Q4uu7l2tL7jL4DfSppo1Y2bBt+0SyioWOsCccFVdcfLbGTA3ZbxJzKKLgmBVrLWDMdmTJ32Bk9AqFeo8tdcz7r55kb27a2Zp89oi/cw3Bk0Ohb6okLy9rLzMGDL2ZdQKKl33xqWlLgoK7UiiU6lJN15Wb2hZ0j6VP1JTGyx2yXmbc99tLV3mDMsM9kIqmiXKItYpzKjaAS/9aiSTcopzN5FcJKqUKl1IVKqZuUUu+f+gxhxkhVulJ8woLwmiES1exo9/Pzzb189g9tfPjuFv7+d8f48N3H+Pyf2rl9Wx97O6fYQd77qCVyj/2PFQFOSwz6n4Dt74KmW63ocZxtQ5v5fss3GYsln7/Us4JrZr2Fa2fdwBLP8qTj/qiP7xz9T7YObOWHz3UnRdTOafDw6ctmU1qYbKGYV1HA56+aQ2Ol+X9Vq7eQXw79w2QaWGcprPwaOLJ4wOkogFXfAEcxDxTX2FKKFXJ9/dunnsOGM1KNHjZTkYWrnqIt2JzzXAR7of9ZLih+kgqH6Vt+eH/6Snfp2Dr8gtHeEPJTMLEpagy6/5j1XFrrpI1y5zXWUOGqNPo2D+ZQCsAutgrroXg+mxrNiOWhnkDW2SOabNHgUmc5s8efSlRfZBWQGEeHrSck2czbZ6YtW1hdaERk15eblfAO+fYxGs1CZEf8MGr7XbEJ4VklrqQiGvu6sssnfGzsqNE2NsqN4/TEN80lMLA5edxpSj4qy6GUKgT+HfgIkHib/OuEMb8ArgUu11rnnuBPyI1UqXPEJywIJ0xMxzg61kTLaDOdwXaCsQAKB1XuahqLF7KiZDVlrvSPeSMxzXNHvDy4d4j+0eTNZJGYptsbptsb5onDXhorC3jzmkrOafBMfjDHwnD4P6Hzd7muHtp/BUNbYM332RZo4Set3yKGKWDnFTbyvoYPs7xkldF/0LeX24//zCgzGyXKT15qRg/WGmNX1hXx4Qtn4XKmf7xbWujk05fP5uuPdRqp1F4KXMqy0X1cWfIQLP1sbmnFihtpXfi3bB16xOh+g7Mu941ewBOHR9A6IaakwlC1mXu7+vj04tysB3Q9AERxqyhvKPkDv/N+cOLQ9rZRur3hrDdJDYUHOew3q5WdW7QQ/B2THcfvgrnvnLpkM9DlDdNjSw23fm4JgcAl/KXvDxN9m4ee4cbZ70pKO5cSmz+YyvNBKdbNLabIpSaeHmisqPDVK6f++STbIlZM/l24yqBykxUJHqfvSZj1hinnbe4zA0VLa83P0NWl63ApF5F4lbYoUfZ4d3Je5RRV2nz7MDbKKddE6rSJLqVYVV/Ec0cmI/i7O0an/PeI6RgtY6bIXuhZknpw1QXmjUmWNwinAzMuhJVSBcAjwEWAH3gKWAPU2obei7Wx7u3Af870us54UkWEpaiGcCYTDVgbRvzNEPFa0SJ3lSWyytZam60yMBge4LG+B3lx6JmkLASJOHGyumw9V9Rcw9qyDUZU6dhAkF9u6aN1KPu/xbahED98roez5xbz/nNrqXL7Ye+nYGhritEOa9NQxXpwFFuRqP6nwR658h3gwM738lNPTZIIPr/yIj7Y8FEKHckpuFaUruaLS7/OL9q+z/aRuMgZXocePN8YV1fq4uMX16VOiWajvMjJJy6p46t/aWcsMvlvdffIX3NWrWLO7BunnMPO/VHz5+OJRbmm9xlY0AKehVnPE4zEeKrJa3ZWbgPXKHt9u9jve4VVpWuzm0xr6JosqXyZ52Ee9L+bsZj1f7UG/rx/mA+cZ//oTM324c2GLaLIUcyaOe+F/gSBM9pkVRKr3JhiBhO7LaKy2Mn8qgJU4FJDCA+E+zjo3zv1961jyRHhKuv3xO10cE5DCS+0TAq/LdkK4dEUG+USqbncFML9z0AskvGJQjgao2XQLoTN3/8iZzErS9eyx/vyRN/OkW1TC+GRV8x2yfKJjXKJnD3XYwjhg70BxsIxit3p/4a6gx0EY2Yke0FxGiFcfSEc/e5ke+wYjLVBcW5WoVORfESEPwFcDDwLvENr3aWUepZkIfwoEAKuRoTwzCPWCOE0IBLVtA2F6PaGGQ5EicY0Loei0uNkdpmbuRUFuDIVM4iOQc9foPfP1oejTvf41Qnla6zctHXXgXsyouuP+Li/+y6eHnhkIhqUiShRdnu3s9u7ncWeZbxzzodY4lnOn/cP8/vdg0ynsBrAro4x/u2hdj4868es0SlEcM0VsPiTUGL7IIz4of3/LPtEfOPYcWcB3y8qI2LbbX55zV/xnrl/mzHaV+Qs5qML/plftf+EZ3u2QqdpN3A74eMX1+PJUHLYzpwyN7fMvovvtb9zoi9MIT/t/zBfjIEr+6loHj3ETq/p3/yrwACeWMiKoq/7aVYRUoDNLT78IfNGgZpnJ768t+sOvrDka9ltmhveBmOTHmyPY5TLFsZ4+MjkkOeP+rhhbSUVWWRQSLJFlJ+Lu/oiKJ5vXIeOu7IUwqYtYt3cYpRSzC9eSEPRfMM/vnnwmamFsO+gmcEBoPK8iS/PX2AK4aP9QXp8YepK00fER6N+OgJtRt+yErsQvgyavj7ZjgzByK6M/wYtAyGiCT9mBSyuTX6qur5skyGEX/HuIKqjKbOnTOC1+4NXpxx21uxinA4m1hGNWfaIjY3pb9Dt0eBKVzVV7urUg0tXgbsaEm/i+5+FhnenX/tpQj6E8HuAMPAurXVXukFa65BSqglI4eQWTjrKhZFEG8QaIcw4Wmv6wj0cHW2iJ9TFYLifUCyERlPqLKXSXUNj0QIWFC+m1JW6LK0vGGVbm5/tbaMc6g0QzpCGq9itWFlfzLmNJWxs9ExGIMPD1i7p43dCZDjt+ZNErQ/LkV3Q/P+sD4eGD/Kifw93dNyGL0UqsWw4MnqYrx36CrW9n6Svvy7lmCU1hZy3oIQltUWUFDgYCUQ51BPghRYfnSPmhhl/WPPtjr/jprICrhnfcOUohuX/CvVvTi3wXCWw8CPW4+E9n8QbaOO7pfMYc5gf3pdVX817596SlahzKAfvm/dhdr+yieGoufGprPEx6spz3CLS8zDro7/hSk8Jj49eP9F9bFjxh72DvHVdmg/3FNzX9VujXRqLcFUgvrlqcLPlq667esp5YlrzlwPmz31ujZ+Owp6JdvPoIXaObGNDxbn205PpvNdsexZx1ZqlPNrSTiQufiIxzWOHRnjb2Zm/38HwQJJFYFPFhVZWoLk3Q/N/Tx7ofQRCn4OC9HmTx8IxDvWYkcV1cyd/rhdWXcbdnZNZB7YNb+Y98/6WwkzV6wZtHtTiBYbFZdXsYkoLHPgSbjS2HvPzxtWmJzmRJv8BIwruUu7kCGjxPCvqmlhkpe+JjELY7g+eV1mQMhJ7dvlGbu/42UTbH/XR5D/AitLU4hZI3ihXnvoGotjtYPmsIqPE8q6O0YxC+Ohok9FOa4sA63ej+iLonozu0/+0COGTxHLgsNa6Y8qR4AUy/KSEk4ZSlk848bGJCOHTjpjW9PsjDI1F8QWjE9FGT4GDskIndaUuClwZHk/rmBW5GdkFvv0wegxCPdZucx2zfocKqqBorvXhUr4WKjYZNoKYjrLft4ftwy+yy7sto21gHIWDxZ5lbCg/lwuqLqHKXcPx4RAP7xvmpVY/kSzDpmNhzcvto7zcPkrJdgeXLy3hqrK/UN7xw+RoVLbExgi03sZvBh7h+QJP2mGzC+eyqHgpFa5KwjpMR7CdJv9BwokWpHA5tP41fYFkEdxQ4ead59Swqr7IEJ/1ZW6WzSri2lUVvNTq586XBxhJSLavcXC396/xxsp5e81DqHU/hLKzpv6+SpYQ2fBrfrz/o/QpM7K9MRzgvVVX5ZQO7IlDPoaHZ5md5TsZ8PyZ24+H+OvGv89uoojPqgQHvL38V+wLrqczOvm49sF9w6yb42HprKmrpe33vcI+326j743hMMWJAYHmb1qCwJX+Zwuw6/hoUgnom9cu5I7h2fSEJmM+93XfwdnlGzN7ZiNeS4AnMvutVHncXLCwlGcTHok/edjLdWdVZnwkvt2WLaLY4WHNeKWz2TfAke9NWuF0BLrut9LLpWFv1xiJ95suB5yVUPr3/MqLuafzdnTcRhMSlxSnAAAgAElEQVSIjbFjeAsXVF1in2qSAbOc8rgtYvIaik3zSwzryZZjvoxCeL/PFJWLPUtxp8q/XHuFKYT7n4Ql/5z2SUDTFP7gcWoKZjG/aBGtgckNartGtqUXwqE+CNqkUVn6RHlnz/UYQviVjjFiWuNIs+4kf3A6W8TEN3CJKYSHtlqbbJ2Z/xZOdfIhhCNAtpnAa7B8xEI+cBSIELYRjQssZ6bH6amIhSHYBYGOuEgMWX06ZlXucVdBQTUU1KX0f2mtGYkM0RY4Rnugle5gB/6oD3/Ux2jUT1RHcasCylzl1BbMoragnlkFdSz2LDMqYI2GouzqGKOpN0DbUIi2oZAt72cydaUultcVsaKuiJV1xdR4lGUT6H4QBp4zH5WlItRtJYTve8JqKxdUbGJ41pU86Yjw7BSe2VRoYjSPHqR59CC/a/0zVYPvZLBvAdZTjOnhD8X40z4vj6hzeGPJNfxV6f0UKPvvvILihVZ+XeW0PqhGjxiWiWPOQn5SOpduZ7K9yKVcXFR1BVfWXss8e9J6YCw6xtah53mw9156RxS0/g2Eq5K+++tXV3H96sqMG8kcDsXrFpayZlaEXz3+DNv9ZiTpYf/b8dfezPtL5medHuiOnns4YBPBiyJj3OJtw7H7I3DO7dbO/iloGwzxu122n7lrCOb8HhQ8N/gES0tWcEn1VVMv6ugPJlJsFaogf1f1//hq/7eJauvfRmv42Yu9/Pu18yjMcFOnteb3Xb8x+ipd1VxefwPs/9xkZ7AbWn9m2Ugy8OcD5pOExsoC1s4u4Yaid/LTtu9M9LcHWtk6/ELmIh3dD5n/FysXzH4zANesqjCE8Gg4xrPN3ox+2W02W8T68nNxj9vh3JVQ91dmxoiOe6DxQ2nzyO+0lcFeUVdEUYIQr3JXc1bpWvb6dk30PTvwWHohHBmF4e1mX3Xyv8/5C0wh3D4c5vhwiHkVKax9WDc6iawsSWPPqLnMLM091gqjR6FkcdJQrTXNvWZEOJ0QBisqnCiEd3q3pa8yZ/cHO4qTrUuJc8/zcOfLk39Xw4EoLQMhFtckryeqo7TaMkYs8ixNOzcAVRcCThi3ROmQtaFxmlUNTxXyIYQPAeuUUrO01r3pBimllgBLgRxyrwgnRD7LLMfCEB6E0ACE+61onNago0AMHEVW1R93ObgqLSGScCcf0zEGwn10Bo/TGWhnKDzAWGyMsegogVgABwqnw4VLuSh0FFHhqqTSbfmhZhXUW0n9tZPO4TDHBi2vWZ8vQp8/wuBohEBEE4zEJh5Bup2KIpeiyO2g2uOirtRFXZmb2WVuFlVCVeQAyrsThndaEdNgN5C5XKyFAzwLoHQlQ8UL2OsuZW90kP3+PQxHck+PBFDNEioDFxMYWkLXQHHOHtMeX4Qen29iI0aDu53zix7n/KLN1LpyE7AArQ4nj4Zb2NJ3N5ETKdASLofeq2DwfAbJbAKdVeqixuOiwKUIRTS9vnDKjAsAIV3Efb738ezY1dxcdhsbi15A1VxkWQdqLklOIh8dg6GX0B338ahvB/d4aommiMCcU9DAOxd9gdoMQrHYWcwlNVdRFXod39/dRSRq+76cfmi4ncis1Tgc74Ypvm8iXkoPfJSPlb3CQ46383uv+YH7zLEYo7EebrmgDncGUQ3wRN+feWrAzKRQGQvzD97jVsqtYCfs/gis/5XhkbYTjMT4nxd6Jv6WwLp9cTXeS9g1ueHq9uM/Z37RosyPa737rLzHCSyon8tb5lRx7yuTfy+9vgj37BzgvZvSbyTbObKNI6OHjb43199EQfUboOs+cwNV2//CnBstP20KmvsCHO41/7+8ZmUFSinOq3w9D/bey/EEz+z9XXexqeKC9F7RLpstoubSCavCnPIC1s/zsPP4pBh95OAwVywvT+l/Hwz3J9kizq28wBw09x2mEA60w8ALUJMsRiMxza4kf3ByhPDi6isNIXzAv5fuYIdRUGWCoZdMP75yWRkjbCybVURlsdMoxf3SMT83rksWwr6INyln7qrSNNHVsrOsoERo0sZC36NQ8uGkoT2+5EIa9o1yiawvP5c/9kxma+kKdtAV7JhM4ZbI8Mtmu2y1dQOehvoyN/VlbuNJxI42f0oh3BloJ2TbAL+gOFnoG7jLrTRqwwke+v5nRAifBH4HbAS+rZR6v9Y6Zh8QzyzxYyzD6p15WJMAyUU1phkRjkQ13mB0YtOIOzqMa3QPBaMHKQntwzF6xNp9moVQjAI9DjedziI6CqvpdJXS4XTRqcOEshKa49+LC4KzYWweBBrirzmgs/uVD0c14ajGG4zR64twsMc8XuGARW7P/8/eeYfHUV19+L3bm7TqvbrJveJucAEDppfQAwkhgSQQQvIlIZBGypeQhBRK8oUk1IQaeuim2Bg3sI1x77as3rtWW+/3x6yknV2tihs2vu/zzCPNzN27s6PV7pkz5/x+DLekUmJJpMhci2mAZKUPwS6Tja2yna1d2ykP7Ov/Af3hTYfW8dA2nkZPIUMPV+NT7s+j3P8lnm/7EiPNW5nnWMp020pshq64j5HAFrODN20pbDf3r64AmrZnob2YVEs6doMDSYi2QCvV3irK2usI1M2Dxrkg49xMEkFGZ8PC4mzGZNr61KFt6fSzedc6Vu9vYHtX7K3J+mAmf22+g7FpkmuL8+PLUhnttLon83Drh2yS6TG7bTLIdR01zPLtBH8QSn6hZf/jsHJfG49+VE9QRh2zpQ4KHgJrPW/U7abKW8GN+d/GZoxj0xto0wLTts0IAee6niPJ0MgjLbcSigig15V14gvWcPO8+CoNm9s+4cnKh3TbzMLMt0JGkmTE/13HbtjyLa2hzNh3VuyZTxqpjKpdPmu0m6Lis/nbwV6Xs4D089fSe/jpyN/1XQ8ug7DrFxCpWmGwwsg7OMeaxKYqj+6W9Xu725iW52RMVuz5CskgL9ToA+oMSxbzUhZpt8NH3AHrLu0NzqRf01OecH+fr/GtqGxwst3I9ELtfW8QBi7OvJIHSn/Xs7/GV8mqpuWcmrIodrK2rbFGJtmX6FaXjHHrAuHGziAfl3YwuzjWtCPaRMNucDDOpdc5JnESOEugI0KptPKZPgPhXbVdMQ2B0/Ji/8enJM7AZUygPdibwf2g8V0uy742ZmxMWYR7Sp+qLAYhmFHg5O2dvWVMa0u1hsHoMp2dHVt1r9siLAzrQ9sa0DLfaQu119xNzRtQGBsI74yqjU60GUl3xf8eKbQPw21K0iU21res4dyMS2IHt2zQrw+iaXFqnoM3tve+/9aXd3LppOSY87HPo7/oSzWnD04eMPW0qEB4GciffK5dZ49FIHwfcB1wFTBcCPEY4AYQQiwEJgA3AWOADcDDx+CYFBCrJTzIQLjDF2RTpYfNlZ3sb/RR1+7vIwtZBBRhZBFJxkaSDQ2kGOtIN9aQaqol0ViLsDTSaWmh1mSkymilymCh2miJyrYF+unkDxO0QVcOdOVqPz254M1kwGzaYdASSmGjdxYbvbMAMONluGUnIy3bGGXZwnDzTsyGLkqNNnaaHWwzO9hpsh96hlQK8ORB23gtAPYNfHs6EosREmwmTAaBlJIOr48O/8DHsts/jt0t43iq9SZmpBzgtNwOitMSEeZEEEb8wXbWtK7j7Y5NVMj+Bd4LA13M8LUxSTjJzrkYkX2Jrg7T4w/x9o4Wqva0EIhXzmHo0rryU1ayw9SONTiNrNDVuCjSj/NU4N7zK+Y1r2BeMlT683i94wus9ixERhUKbKsX/OT1Cs4f5+bsMUkxmdOtbZ/yz7L7+szYFwc83NReRUa3K2PDMi2gGnN3TM1jKCR5cXMTr22Lbc5zJTTSnnM/mHqDnY2tH/PrvT/i20V3kGqJCsD9LbDpppiO87lJm7GXmPjbOqnLyG6q9HD/ilpumZcRUxN+0HOAv5beEyOTdn3ezRQ7x8En10GkTWvLetj+Qxh3T0z2an1ZR4ycWEGShYsnJmM2zmVv5y6W1r/as6/eX8s/y+7n1qIfxtbQVv4ntqO+4Eawa6UeN8xK52dvVOhsmB/+qJ5fLMmNqZ9d2/yhLkMLcFHmlZi6zRWcwyH3Gk1DuZuG96Hhw5jgsLbNz/qoUoHFJW5ddnZK4gyK7MN1NZqv1DzLrKRTY2tWK57Rr1szIVkvuTUy3caINKsu8H9jezOzipwxAdCaphW69SnuGbHPKYRmLbzrFxGvd5lWIhCVBd9Qrq9WLEqxkOKMDR3MBjNzkhfwdn1vpvnDxve4OPNKTNHP3/ihfj05vsTYjEKXLhCubQ9woNHXY2/cTXRZxEjnmL7rg7vJWKIPhDv3QPtucI3UDdtRo/9cG51h67dW3iAMTHXP5P2Gt3q2fdy8KjYQDnpjL4ASo0wt+mBavlMXCNe0+als9ceUi+zp0NsxDFgW0U3qAtj3x951X51WwuGeNLjHn4Ac9UBYSukRQiwG/gPMBmZE7H4n/FMAa4BLpJR+FMeG6NKIAXSEGzsDvLm9hRX72gasO+0miJmGYCYNwUxNOySGEJhawdII5kYwtYPRoy0GL1qusdsJyQ7B8BJIBF+qtgQ/ez90P1Z2+CaywzcxvCWIsFUgnftA7APLARD9OXtpWGSIvKCXHL8Psy+Ftq4SGnzTqeocSVdg4GYgQBP0d+4Bx36wVYKtEp+pjUxbATOCIaY2bcTlr8cbslIVzGO3byw7vePZ5RtPu+z7lneXtPNBwxg+aNCauCYWBvEnruGjtrf7LekwyRCzva2c7m2mINj9Jd4Ie++G0r9C9mW0pF7J+2U23tvVqusQ17+mAKR8CGnv6yxdP21bz6a2DcxKOpULM68gw+iEg/+E8n/rSn1yzOV8NenPnOl6kye9P2JXq77hJhCSvLi5mTWlHVw3PZWSDDueoIdnqx5jeWNUE1OYJWnncXFrKabWKNMKXx18+lUovBEKvwEGE52+EH9fXRujxQowLd/BV2cWsKr1Gp6sfIRgxIVfeVcpv9zzQ75VdHuvc5uvCTbdqDUvRmJOhckPM9VRxHecHu5bUaP7P91S5eH+FTXccmpmTy1tk7+Rew/8OkZr9LyMS5mVfKq2MvFB+OSLeivc+qWw57daJjUcFDR0BHhkbb1uHotRcNPc9J6Li8uyr2V/5x72RGi9bmpbz+t1L3JexqW9D/TWw757dXNhL4aC63tWMxPMXDopmac29N4PaegI8OwnjTqt3UDIz0s1+mAzz1YQq+9a9HWtXCCynn3Pb7ULmoiA6u2dLUSa+dlMgtOG6zPaQgguybqKP+7/Ve+x+etY0fgui9LO7h3ob4Ha1/XHkX1Zn5q2S8a4uX9F762p8hY/W6o8TIgoU6joKmN/VCZwhjtOkJlxHuz9IwS764+l9n8z8s6eISEp2RAV9E/tIxvczWkpZ+gC4bZgKxtb13FKZGlGZylESZz1VR/cTXGKhQyXSWfm8eH+tj4CYf1F0+h4ZRHduKdqFx3emt5tta+Dq7c2XErJjqiMcEnGwJ/D091zdIHwwa79seURbVt75Ao1xKCCzaIUC8kOI00RpV8byjtjA+HOaGORKBm5eDiHgWOY1hvRTf1SFQgfLlLKSiHEPOBc4BK0LLAbaAe2oZlpvCj79QpVHHEGWSMspeTD/e08tb6hx+nnCB4EBJK0hQHql44U5kYtQDQ3hgPwpnDg7dMWgKAVQlYIOsCfAt5UjL40Ql25yNBgOmiNyK4C6CqAhgXhTe1grQNLPZjawODFLLrIDHlJ9Ztw+S2EQi4qA/l85C/ARz/SQ1EYjD6kaxsyYRO4doIx9m+5zVvKNuBfrhSGB+yMDnQyxl/JAvM+Fjv/S0gKdodmskZeybrmYXT4+856lLf4Kd8EiCmQaAL3enDuBUPvF1Wiyc0i96ks8LWTWP2Slv2IQEo40JnJ+5sdrPG0EaDvsgujgFOHJ3D2GCfbvLm8Vmunwa/PUEkkq5s/4KPmFcz3dXB+RyVu2UcZTcqpFIz8Ebfbcvm4rINnPmnUfZkAVLX6+e271UzI91Hu/j+aZFnMNImmJL6WfyvjEiZBDpC2AHb+NKqpUGq6vE0fUZV/N/evlVS3xV4JnjU6kcsmp2AQgoWpZ5NjzeeB0t/REextjmoNNPPbvT/lhvxbmOkYDZ9+TStRiCQcBOMoBjTpqf9ZkMUfl1Xr/me3Vndx7/Iavj0/k5Dwcu/+X9Pkb9BNNcM9l4sye/V6sefCxL/BJ1+KCJrQanetmVBwA4GQ5O+ra+n06y9krjkllezE3s8ZkzDxjcL/4ee7v0drhGzdi9VPM8w+krEJ4QvJPXdDUJ9ZZtSPYz6zTh+VyIbyTt3t6+V725iW72B8tvZ/+n7j29T5anSPuyTr6tgMtCkBht2m/S278eyHiicg/8sAtHYFdaYGAPNHJOCwxN5dGeeazCjnGHZFuLv9t/Y55qYs7JUVq3kltkku+1L6YlKug6wEs+599Nq2Zl0gvKppme4xblOS9j7tC5NDe67ILHjVi1B0s9bgi6bdG1mfC9qt+Xjk2PJiXvPyxqX6QLj+Xf2DLGngilPCgHZRMbvIxctbei+41x7o4MopKT2lPo2+eqq85brHDahjLAyQfrb+9de+AcW39lzc1bQFYl7/6MyBA+FRzjEkmpJojUgSfNy8Sm/hHV0W4RylvQcHwCAEU3OdvLu7N0u+vqyD8yPUNLQyM70aRYyxSH+knQEH/967XvcODPufnvMipaSlK0jSIPSsTwSO2asIB7mvhhfF8UBMaUTsF3UoJHnko3pW7m+P2XcikGksJ8+yD7elFKOtHI+9igprkDKjdchlCkHQShR8qeApgM5C6CwGbx9NEH1O4IJOl/aYMH6gPLwMFbehgam2tUyxrWG0ZTNBEWSzwcmGQAIbDU664jRdBIVgl9nBLrODV+xglJIsYSbDmk+ycyROYzXzZR0V9S5KK9NobsqAvnQHpBlapmqL8INjHwlJNczJHs6inGmk2G2a+kbh1/FWL6XhwOuUtwTY5RvHFu9UaoPxz5sAZhU5uXB8Mhnh2t0M15nMTV7Issa3ebXmOdqitHuDSN6zOPjQPIwzupo4vauZZBnQvmhH3AHpZ4IQCGBGgYuJ2Q5e2tzE0l2tRF+Cby6zQOUNkPma5hQmtAHjXZO5If9buM0RGeW0BZD4Amy/U2dLKiUsq8ri2V2teKX+y9Mo4OppqSwcqc/Al7jG8eMRd3Pfgd9Q5a3o2R6Qfh48+CeqfF4u7Dig186wZMDkh3qC4G5GpNv4n4VZ/HFZDZ6IAHVHbRd/WFYJBf/UdbeD9mV5Q/4tsUGiqwTG36dloiNLlfb9CSzpPFU2J6Z5bEaBk3l91LAmm1O4qeA73LPvFz1yW5IQDx78E3eNuofkxrWawUkkGefGlJqAFhRcPzONn71Roct+P/JRPb9ckkvA0M5L1fq2k+GOEiYlnBIzF6BJi1U+o79lvf8vkH4W2LJ5fVuzrhTDKGDxqL7rLoUQXJx5Fb/d1xtYtwSaWFr/qpb9lkEo19ctk3Y6WGPr0Ltf69lj3Dz6UW/WfVedl521Hkoy7ARlkFVN+l7z2cnz+zdzyPuilgXu7r8IeaDyWSj8GgBrSvUXnd0mNf1xWsoZukB4W/sm6n21pFnCEoHd6jLdpC4csP50brE+EO70h9hQ3snMQu39talNH1Q6ja6BG8MAMs/RB8Jd5dC6UatZBnbU6u/guG2aUc9AGISRU9yzeK+h9338cctKfSDcGtUo5x64LKKbqfkOXSB8sMmns9+OLouwCCv59qJBz096VCDcVa7Vk4etn+s7Atz+3/JexaF0GzOLXP2bFx3HHAuL5euAGinlW4MYeyaQJaV8fKCxiiNAdP1UVEY4JPsLgoNa5jFhO9gqtCyn0QNIkEatKS1k0zr/A27wu7XMqi85/DMFBpVZHRxG/OSaDlJg3tez5Jv2YzdE3Yr2aksAqDRaOWCyccBo44DJSvlggmMhwVqvLUnhD9+AHTzF0DFMC3I9eRyt+uQsYzlTbGuYaltDsXkXBhEhHg+c4m/nFH87vg7BZrOTjy0JbLS48PXzuoJCUEGACu9+8OoDI3KAdDc0nwJNM7W/XV9IM3SU0NZRwlsV8BbVCDT1jZCEQGgMWhvAwEyxruHihH+TZ3ZA5XjtNp0lFQxWzCEvi701zDNksdRTz5tmY0zA7xMGXren8qYthYnGJGZkXcH4pDm4our6bGYDV05NZVqBhX9+VEVdS9T/Q9AFlVdAwwJM6R9weclEFqWf2bcerCVNy5qWPQr772OvdxjPtN7AHn+sfq/LauDmeRmUZPTdBJdpzeZHI37D/5X+QdeFD/CKxUqlK4cvd1TjkCGwZmtBcBx1g+FpNr63UMsMRzY87anzQ/tCKNjXU2qSYcniW0W398psRZM8A0b/Grb/QLf53Q2reb9Fn31Lc5q4bnpa3FrKMa4JXJx1JS9U9waCbcFW/rr/1/ygdrVeb9OUCCO+3/cxARkuM5dPTuFf63oz202dQZ7e0Ai5z+IJ6W/tX5Z9bfwaT2GAEXfCJ9f0bgt5YNcvaBx+H+/t1mepZxW5+qyX7abENY5xrkm6v+PrtS9yWsoZJDaujS0RyLki7lwAs4tcvLKlicaIOxmvbGnm+4vsbGnbSEugSTd+bvKCfufDlq2Zh9S+0but4knI/zJBTHxcqv/sn1k4cCPsKe7ZPFn5MJ1h226J5IPGd7kk6yqt5KVV/54m7fQB50xzmRmTqTeT+HBfe9xAeJxrUv8XAN24xmpyiZE18NUv9wbCNfo7VaMz+68PjmS6e44uEC7vOki5p5Q8e6EmqxmtGOGeOqh5AUal20i0GnRqFmtL27lgvCbFGF0WMcwxsrcefjC4xoA1R69xXPdOTyDcfQemW3FoY0Unc/q46D1ROBYZ4UfR7JUHDISBO4DTABUIHwtiMsL6QPg/Gxv7DoKTPob0t8HSFLsPQAQRSFwWKwkuSYJJkmAKkWQOkWExk2FxkG5JwkEaLR5N+qihQ5My6/CF6PRpChS+gIysEMZmEjgtBhwWIwlWA+kuM+kuE+m2AOmGg5i6OqCjDbp84LWANxl8gag6LA2TMFEgrBQY0zjNng+2PAK2XCqsKRwQkgPecg549lHRVTqwba7JAwnbIGEbyeZUcszFOHzjkB3FtLQkcaAx2K/7WX8kWgKUOMoZZ1zOWMNy0kxxFQh1WJBM87czzd+Ot0OwyezkI2sim8wJxKl4iI+5BdLfhbT3oGO4FhC3TRhQgUOCLnvWHyb8zLK/z5nOl8kzh5ua2omtgw1jBy4AFgojr9tSeNeWFHMRExKCjaEWNlb+HSr/TrY1j2xrLokmN0IY6Ai0UeOrpryrlGBuEOxzoHaJdgEXiTeTQPllvN5gpLGwmWn5DgpTrDHZj3afZEvocpZ3ncrOxr7PTYF5P7eM2UNa0mXhV9E3DqOTbxffyVMH7+f9Fn1j0TpLAgeMVr7qDzJq4j81M5N+KE618r2FWfzh/Wp9DbanAPbfAgUP4XL4+E7xj0gwxZdFA7QMmq8e9mqKCCs7F/Jky1d1Q8xGwTfmZvRZLhDJOekXs7djJ5+29erJ7u06wKNWB18NtPZmvUf9RLvQ6IcFIxJYX9bBtshAaX87+OvDrdkas5NOY5RzgAsy9yStTrfqP73bGlfwatNmAqFezWejAS4cH9/coZvLsq9l2+5NPYoGXSEPL1c/w7U1UZlR11hI6t+BzmwULBmTxBPre4P+7TVd7Krr4r32N3Rji+zD+9SyjiHvS/pA2FcHta+zLbQ4RjZsVtHAwY7FYGV20mm829A754rGd7gg4wuYGpahczM1OrULrEEwrzhBFwhvq/ZQ1+4nySHZ1qY3SZmUOLD6AqDd6s+6EPZH1KPXvgkjbidksPXZKDdYRjpHk2RKoTnQWzK1qnk5l9uv03TXow19hpARNhoE0wtcuqzw2lKtPEIIESOfN2Kw9cHdCKFlhcsjQrHaN7SyGSFiXAZL0gd/gXA8cqxKIwZ7hk7cM3kiElMj3Nsst7uuKywRFPEnEQHIewISN2OUkgJTOvkJE0mzZpFsTiXFnEqiKYlEkxun0YlhEFfkbhsUJA++FjY+48NLH4QCYYMLn3YlbnRqrz3qH9eE5u9dCMwPb5NS0h5so8nfSLO/AU/IQ0iGCMkgISQOo4MEk5sEYyJucxIOY2zGJBCU1HX4qWkLUNPmp7bdT5c/RJdf0hUIEZLgMBtwWLQl1WEiL8lCXpKFRJsRGAlyAXhKNV3V9u3aT08p+Br6DPQBLXiwF2FNnMR09xSmuyfTZbCyq2Mb29u3sKNjCwc9+3WSQ/0iJBkp7YwpqGWc3UNjXT4fHfSwv9EbU1owWDKcBmYl72Yh9+EODl1OLkEGucJTx2JvE6/YUllhTULG+RSp8pbH1BH2IIDUlZC4BaouhLaJMUOaPUHe3NHCmztasBgFGS4TCTYjIak1ktZ3BMLnoe+P1fmON7ky8Z9Ya7xQ9w8tqMw4V5OyipQiC/mhdTOm6he5tuZVcixOnnRkICPer/VGC3cbYX79K1yUeaW+VKMPCpItLD6ljBfXJkEw4j3qS4d9t3HBdHvfeq99kX8deGtYsaucR1u+FaPCccPMtJhGpr4wCANfzb+Vn+/+PvX+3iaw1VY3mUE/F3Q1aOcnY8mAcwkhuH5mOj95vVzfx1B5GdjLwdKIzWDnsuxBWjsP+46mohBuEDzoL2Z5vb4E4rRhCaS5Br5NXmAvZnbyfF397vLGt1nYuY883cCvxHU20z3vcBevbmumJaJ29ZmN1ezP+ET3cT0osxKAxPGaG2SkZFbZo6zx67WHi1Ot8SUGo48x5QxdINwSaOajlpXMia4PTjk19nsoDlPzHdjXG3rKfCSwfE8bo4cfxCd7kzgCwfiEwQeVZJ4P+++jJ0APtkP9e5Sazoi5EBiTGf/iNRqDMDI7+TTeqHupZ9uaphV8Ic4dph4AACAASURBVOsaDJGa1aDd0Ymwlx4MMwr1dcJVrX7Kmn1kuw0xjnKDbpSLJONsfSDsKdVKhhLHszPKYGTUEC4QjkeOt0rnLJSz3LEjTiDsC4R4YFUpELFfBCD/MXLtGznbMpopw36IY4AszXGDwQSYDskmUghBgimRBFMiBUOpsYrAZBRkJ1p0TUOHcCDgKNKWzHN6t0up6cn6G7SAXwgwujT3qD4c7GzAxMRpTAxnTDqDHZR3HaSyq4xGfwOtgWZ8IS9BGcRisOI0Okkxp5NlzaHQXozbHOGClgpnjk6mwxdke00XW6s8HGj0Utnqj5sBd9uN5LstlGTaGJdlpzDZghCFEDpNs5ite1vTGA3F1yyOwZJBSsZZfDnrEs4yOVha/yqrmpbrvhwHjbkFZ/ELjMNM5cEJlDf3fTfAF5SUt/ihZWCRm3RTLdcmPsB4a8St0JAHqp7XFmEGe75WAhDs1DS3Q72ZqNO9zWSEfDzozKHToL+4XN64lNVNyzktdTHzkheRbyvUZWYCIT+b2zbyWt3zmplEcTqU3qh3swvZeHKtpKK+nsunpPRr3avNKXmp/Su83hJrUX1RwlPMsE8BBg5eAZwmF7cW3c7/7v4B3gid8JccaSSaElgQoWAwEKlOE1dPS+XhSOWKkB3KroPiB7g45yqSIt+//WFOhJE/gq23ISX8u+XruoDfbBSc14/VbzSXZF7Fx82reuy1Q0iecGTyg7YyLXa15WkNSoM5NKOBc8a4dWoZ++slOEogQasNdRid/dsbR5P/JV0g7GkrZ0N9G5G9AbMGURbRM529iBLnOHZ29NZav137MrObVuuzXWl96CrHwWoyMKfYxbu7et93K/a10ZWqL4sY5hg58J2NSGxZkDxbV99P9YtsNs7RDctMMPf0KwyW2cnzdYFwc6CR7e1bGBf5XKA9/xAZnmYl1WHUGQetOdDBpJFVBCKSIwLRqzgzFBImaO/LrojkQe3rNJpGU9eu/1wcjJLG8cwRD4SFEAUQLeqJWwjR33+lHS0JNwpYe6SPSRGHOIYaT2zaS1tn1L70N7jUtJIleXdgyDjzGB2gYkCE0L60+3H66g+H0cko55iBbxf3g9Ni5JR8J6fka1+UoZBmRNLSFSAY0vIsTouBBKsx/u1ygwUyz9WWkE/TrWzdCB37tDpKf6tWumOwakG+PU9r4HJP1zrOw2UR2cB1eTfxhewvsqFlLR+3rGJn+9YYh6WYp8fAGNcEprvnMDN5HlaDjdBYySflnby1o0Wn3zpYEm1Gloxxs6goCXPZGKjYRJ+mMtKvlyrqgwn+Tu7qaODvKRPZ46/W7fNJH+/Uv8Y79a/hNiWRZyvEZrDRGmilvKtUXyNrrYPi+zVr565c3TzL97bxSUUn549LYm6xS2ehC9rdkZ21XTy1oZGy5tjzebbzBc53PgnbntTcFotuju1DiHntIfKqXuKmtlLud+Xqst7/stgwtq3j1JSBa0i7mVvs4r2DezhQFRGkduWRWHMjiyYuGPQ8gHZrOO0MPjxITK33uWPdJDsG//WZYkljSfqFvFLbW26x0+xgrSWBWb42TZWiD8m0eCwYkcjbO1r07om1S8C1C4Tk1OTTsRqGEJykztfVyq7yLMQb7P37G4TWYDoUFqedqwuED3oPsstopiQQDqIMVs28YQgsGJGgC4TbvCE+KmuFiEObmDDIsohIsi7SB8JNa/jUo5eFnJgz+GxwN3m2Agpsxbqm1JWN7zKuOUox4hACYYMQzCh06TSFV5e2Y8jUl4kU2Iv7vFM5IEJon8elD/Zuq32DXcav64Y5zAbyBmigPN45Ghnh64GfRm0bD7w/wOO6PwH/3O8oxZGjjxphXyDIqr0+dG8Neyk32J9j7ohfDfmDS3HyYTAI3HYjbvshNgwaLJrD0iBcluLhMDqZl7KIeSmLCIT8HPDso9JbRp23hq6QB7/04zS6cJuSKbAXUWgfhj3qjoFBCKblO5mW76S82cfqA1pTSFVr/CywxSgYk2ljRoGL6QVOTN3mHCPvhJwrNY3jmtcYnB13GGGC3KtJK7yJ200u3q1/nRdrno7R/gXtFnRL+wBW3eZWZk3ZRKBiOOsO6udo7QryxPoGnvu0kfFZdvKTLViMBho7A2yv8VARJwN+tvMFLkt4pPfO/sF/QtNH2utOjFOy5GvQpMoaljMZuLqzliecvUYxEskj5X/FE/SwOO3cQdUgbmxdxwH3n6HhVvBl9L6uxmH8d2srF00YZEY4TE32D3lyi963Md1Yw5JiKzC0uZZkXMzKpmU0+Hvr/J92ZDDe4MYVRzItHmaj4MIJyfrsd1ceNE9DJK9nYepZQ5oPYdCywrt+jpTwXue5ut1Tch1D/n+enHgK6ZZMnXTdG7YUStrDaiipC2PtzAcg121hdIZNp+3bXjsRXL0mIpMS4yiC9EfaIu2OTLhutzXo5kCzJLLWZGL2oTV3z0lewMGq3kB4XcsariYQEbsbIGlwddLRzC7SB8ItniAbKht1N3PHumLLvAZNxjn6QNhXx87yMqD3XIzKsGE4QdUiujkagfAB4IOI9flAK7AxzngJeIB9wDNSyg/jjBsQIUQxcAaaaccMYBxa+/5PpJS/6u+xcea7C/jZAMPGSCl3DDDm+CTGUMPP8zu3Ewzorx7npjzM3MJvqiBYcUJiMpgZ4SxhhLPkkOfIS7Jw2eQULpucQmtXkPJmHw2dATy+EAYDuCxGshLN5CSaY5zbenAOgzG/1nRKa1/TmnLad0C8Gm3HMMg8D7Iu7pHUMgJnpp/PNPcsXq55llVNy2Ic4fojw5LF5dnXMdU9E1koWZbRxrMbG2MMcrwByfryzhgHtWgMAr44LZkFsgsqo3a2bYINV2oBT+Z5kDAejHboqoD6d6DiKa0UJMzp3mbaDEZesetLrp6ueoSyrv18MffGXv3dPljTtIKHyu4HYxAKHoV939JKI8K8sqUZh9nAmaMHYTML+IMhHlwXxCv1mcCrEh/EvNOv6TYPsr4VwGqwcqW5gL9EBMKtBhNPpZTwtYEy530wp8jF69ubqG6NzAqfw6RcKxnWrCHPR9ZFUPp3dramUhXQN9lFy/wNBoMwckbauTxV2WsWu8niotRopTDohazzh36MwIKRCXqTC0+xJmXpKCXLmkO+rXDokxpt2v9ZWEpti3cqMiIItprEIdfBzko+lf9U/6vHKCdAkFUWN2d6w83mCeMP+Y5eXpKF4hQL+xt779BU1+RBfu+YwwqEncM1pYj23hBnV50XXSCcfmKXRcBRCISllI8BPcJ8QogQsFlKufBIP1cffDu8HGnKgINx9g1sF3a8EvUhHgp0sWxPJ9AbCJsdO/lSkmtAWR+F4mQh0WZkbNbQb5P2YMuCghu0JdCuNT766rVslNEB5hRIGAuWOFJ1QKolna/k38x5GZeyovFdVjUvjzHGiGSEo4R5KYuYk7ygR0ZJCMHCkYmMy7bz5PqGPl3v+iPXbeb6mekMS7WCvBMcBbD3D7GW6A3va8uAGLkw72aCRj+v1b6g27OyaRk727dxRc6XmJI4XdeI2x5o47nqf/NB4zu9D7DWao29B28gMqv39CeNeIOS88a6+80wB0OSv62s40CjvgTkVPvbTLZ9rKV2dv4cRv9qUA1uAHTsY2rFc0x2prHR0mucsLprL1Nb1jLNHauT3B8Gg6CoeBvVn0Zc4AUScTRdPKR5eic0Q+FNvLNGXwaU5fAxZhAmEn1xWsrp/LfmP7RHmKO8Zk/lm94uSJ7TzyPjMzXPSYqjUSchR8N8cDzODPe8Q1cvyLm8JxDe6NVnaMdm2mOs1wdLosnNtMSZfNSysmfbMpubxd4m7Z2ZPOvQjjfMvGEJ7G+M+N9vGwsBJ5g6MAnzoTXKRZJ5YU8g3Bp0U+XV3wk50euD4dg0yy0EWgYcdWSoRzPs+Aj4GPgqMLR7Tn3zsJTyriMwz/FFVHZlaXsSfo9ezH2+6zVMJT8b/Ie9QqEYPCbXoOWj+iLDmsWl2ddwcdZV1Pmq2de5h2Z/I96QF7vRToYlmyLHcJLN8YPqDJeZ2+ZnsbPWw+vbWthS7elXBSTDZeL0UYksHJHYW/ohBORdC4mTYfvt4ImXN4iDOQXG3oNInsGlgM1g5/nqJ3RD6v21/KX09ySbUxntHE+CKYE6Xw2b2zbqmoO6mVeQQ3FGKv9apy9teHFTE9Wtfq6ZltpnzXqXP8TDa+v4pEKf48gyV3N1YoTJQM3LWqd/WFKqX4JdsONOhPRxbUcNu0wOXePjw2UPUGAvIt2S2c8keup9tWwIPg6ua6C9t4Z5zV5YPNxLYcrQ1XgOWpewwVur27bQ/jJCfkdr6hwiVoONM1OX8ELtsz3b1ptdVCTMIvcQsuAAJoPgjFFunt0Y8XdtHQ++1Fjb7KHgKITk2XgaNvJpl17GblLu4WneL0g9UxcIVxut7DLZKQl4IOUwjhlN2/npTxp7G5SlSdN9T1vOKOcYLP3cRRkUmefBPu0Cd6tvsm6XzSQoSD6x64PhGATCUsrlA486Ys+lK38QQlwZb6yCmIzwyhadmA9Gcy2X5aZpNqoKheK4xSAMZFpzBi+B1gclGXZKMuy0eAJ8WumhLOxWhQCrUVCQYmVEmpWSDBuGeIFf4gQ45UXNlrj073pL5j4RkHkBDLtV9zlzbsYlZFqzeajsgZha6CZ/A6ub+/9aWZx2HldkfwmDMNAVkPxno17zfPWBdrbXeDh/XBKzilzYzQZCIcknFZ08/2lTjB223Sz4xuxkrHsN6CpRSv+mOcRF2PLGIEOw405o2wJAsgxyVWctD7l65bI8oU7+WvoHfjj8l/2Wf/RMKSWPV/xdU0bJegX2jurR9Q5JeGhtPT89M6f3QmWQvLxV//dyiHbmmp6DiixNNu8QWBSUvBEK4gkH/lIIXjCG+NYhzaYxf0QCz2+uIRjsDqYNOJvOJ8eW1+/jBiT3Kj6psOGPsLc3CtmvrfRgKHGOI8uao7M9fseWTInXDu7J/TxyYBwWI9PyHHoXwMY5kPrB4ZVFdGNJ1kxP6t5ic5e+/np0pl1zDz3BOWbyaUIIK3AVcCaaOkQC0AbsQjPbeFrKQ9E7UhwyEYGwBKo6Rup2j3Kuwlx4wzE+KIVC8Vnitps4bXjCwAPjYbRqmrg5l0PdW1pzYNt2iLg9jr1IUyrIulBT/eiDU9yzybXm81jF33SWvf1hNdi4PPtaFqSc1XOLfMmYJEwGwdMbGnXV2M2eIP9a18C/1zWQaDPS6Q/1KftnMQq+fVoW+Rk2sP4OttyKrq774D+0DHjJXWCKOm9BD+y8S5MFjGCOMYkt7tmsbVnds63Us5cHS//IzUU/GNAVbWXT+2xpC8vxWesh7V2o622QK2/28fKWJi6dFP9OQDT7G7wxWfCznC/iMHTCgb9CxllDT4pIiaPyGRaHmnS13590bmd3x45Dvm0fFB3I5NVQ39u30tkwjrp2P+mD0HaOS+p81vr052Cccw8u6yDsmvtBCMGClDN5uurRnm0bzC5qXaeQMRgHvAFYODJRHwj7U6F9NONckw57bgCyLyFUu5TNXr373aEoaRyPDOAne2QQQkwFtgEPAVcCU4GR4Z9XAg8DW8PjjkcWCiH+I4R4TwjxnBDiB0KIQ+hG6EUIkdfXgqalfGyICIS3kknQq88mzUv1ajJVCoVCMVRMLsi+VGsqm7cKZi/TllPXwcxXNdvkOEFwN9m2PG4f9ktuzL+NQnv8YMQoTEx3z+GXo/7EwtSzY+pEF5e4ufnUDKym2OyVBFq6+nZ/tJsFt5ya0dsolbYQRv449gDq3oKPLtAstjv2QlclVL8C66/SGiN1B+tCjL+P6/JuJtOi/8zd2LaOh8oeIBCKr0yyu2MHj1c8qNuWmLWBHLc+oHptWwsbKwbXwhIISR7/uF63zSlaOcP5X20l2A67fsWQnXMaV0LHbs7yNOEK6WvHn616HHmITjwfNr1HKOV9EL3nSUrBf7cOoJgyAG0+2Nall5KcaXpZu5A7TOY5J2CTvXXNUgiWmo9MIDkizYrbpbdgMDctpMBefETmJ3k2e5lLh9Q39U04RCWN442jnhEOB3dL0bRm6oF/AFuBGiATTdnhq8Aw4C0hxGQpZcXRPq4hEi2XcClwlxDim1LKRw9xzrKBhxxlIm7BLfPr6xSFwcP0/CG48ygUCkU8hADroRnwCCGYlXwqs5JP5aBnPzs7tlLtrcIT7CDJnEKmJZtp7lm4orOxUUzNc/LLJRYe+7ierdUDG7YUpVj4+twMMqIzjLlXaJq/O+9Clxn21cHee7Ql7osxwbg/gXMEduCWou/zm70/pjPYG8Ssaf6AlkATN+bfFuMYuLN9K38p/X2M7fsX868nLSeT/327ksh4/p+r67hjcTa5A+i8vrm9hdImfWPgkqwd2IlooGx4X7PZjTT06Q8Z1GpLATshzvc08FSEPN7ezp2sbHqfeSmDN9UACMkQ7ze8BeY2SFmlNcqFWbW/nSVj3IdsXrTmQDtB2ZsftOBlinUtHDDAhPsPac5uHI0rOa2rhbftvVn6FZ3buTDQNuB7dyCEEBhSV0H74p5t/rZhVLb6B/zbD+4JDGwyXa7blGc6QKrJgRbGndgci4zwHWhB8AtAoZTyR1LKJ6WU74Z//gjN1fZ5IDU8/nihCvg1MB3t2BzAXOANNBOQh4UQh6YBczwQ0fyw06sPetPs2zGlH14Rv0KhUBxJCuzFLE47j2tzv8aNBbdxefZ1zE9dPOhAIs1l5rsLsvjO/Ewm5zr6LOstSrHwtVnp3HlGTmwQ3E32pTD+Pk17drAYnTDuz5DSa56Qayvg1qIfYopqRNvevpk7d36Ll6qfCVuib+ZfFf/g9/vu0ikwAJyeeg6nuGdTlGLl8in6UohOf4h73qvWar3jsL3Gwytb9DXUOYlmFs86U2tijGTXz6HzwOBeb9UL0LG7Z3Wht5l0g143+Nmqx2kLxDoU9seGlrW92sRp74PoDeBDEp75pDHOI/snJKXOrANgku0jbIYu7SKgdfMhzQtomfTqF1nsbcIQkQX3SS9L61899HnD1HgraXK8C0Z9Vvj1bUdOp2BT23Dd+gTrOqh4+ojN/1lyLALhJWi2yV+WUvapzyOl7EIz4ugABnm5efSRUj4YDtzXSSkbpZQeKeUq4FzgRTRdnj+JQ9NryY+zTO/vQUeUcEa4GRMdnfrbQZMTK4csdq5QKBTHO0IIJuQ4uPW0TP5yaSE/OyuHb87L4IdnZPOHC/P5yZk5zC52DdxolrYQTnl+cK5gjuEw9SlIWxCza5RzLDcXfh+L0DfJeUKdvFL7LHfv/Qm/33cX7ze8GaMZPTFhKlfmfLln/YxRicwo0OvAt3QF+c3SSrZXx3797qnv4r4PaghETCsEXD8zDbMtWbOZjiTYAVtug4A+4IrBWw/79RlUk7OEq/Nv0W1rD7bpdIYHIiiDPF/9ZMSk7SREOaltqvSwqXLoqqabKj3URlkHL3K83ruy/74hz9lD21Zo30FqKMB0n/5CZmn9a7QH2uI8cHCsb1kLhgCkrNRtX1vaTl37wDbwA1Hd6qesRW8CNMn2MVQ+M/B74QTgWATCOcB2KWW/7cPh/dvRXFKPa6RW2PTD8OpwYMitmVLK8r4WoHrABx8pwoHwJlkMQX3QuzD3xL/doVAoFP1hMxsoTLFySr6TUek2kh2moenQ2rJh0j9gyhOQfhZYet3sEBZInKKVQkx/QTNUicOkxGncPvwXuE1JccdEMyFhCjcVfFfXWCeE4PqZaYxI0wfVrd4Q9yyr5h+ra9lU2cn2ag9PrG/g7neqYsxUzipxMzwtXBOdfqam6hFJ5x7YfLPODEVHyA/bvgf+qMzsiO8zyT2dU9z6C4c1zSv4sPG9Qb3mFY3vUuPTO7dcNiEXt01fH/3UhgZ8gcEbzQAs3anPnuY72xll2dK7oWk1NKzgkKh6rufX8z0NiIiscFfIw5t1Lx/avGgKIh82hc9fyodg6NUcCEl0znOHyppSffiWYGhmuHmHpn1eeeJnhY9FINzG4ItIMtGywsc9UspdQPd/+ojP8lgOmXCz3L6g3kXIaGwhK0+VRSgUCsWgcE+CcX+AOe/BnOXactp6mPovSF8Mg1AGKHaM4K6R9zA3eUG/4wSCCzOv4NtFd2I3xjZbWU0GbpufRWGKvjZUSlh9oIM/L6/h9+9X8+6uVkJRvWrT8h18YVKEYYIQMOon4NQrCtGyDjZ9Hbw1+u0hP+z8mbY/ktRFPcYRV+fcgN2gb7L6d8U/KPMc6Pd1t/ibealGH3QV2YczJ3UWl07SmzzUtAV47lN9uUd/7K7rYnuNvm588dgchFWvq8+e30BIX0s9IIGOsK26Rk7Ix8yoJsl36l+n0RffEKc/tJr58MWBqROS1+j2r9jXRu1hZIWllKw+oA+EZ9o+wCjCFxplj2nqKCcwxyIQXg/kDaTpK4S4Cq00YF1/444zut9dx0yG7ogSDoQrgnplCIelFuy5n8URKRQKxYmNJVVbDqFizm1O5ob8b/HjEXdzeuo5ZFg0ESGLsJBqTuOc9Iv5TckDXJh5OQYR/+vbYTHwvYVZQ5K3Gptp48bZGRiidWGNdhh3L5j1wSYtG+DjS6D0H9rvtW/Ahmug5hX9OGs2lPy8ZzXJnMyX8r6uG+KTPv64/5fUeKv6PLaQDPGPsntpDeizm1/I+iIGYWBOsUtzOIzgnV2tfZaDxM4teXKDPghNsBqYWZys6UNH4jkIBx8ZcE4dFU9AKPI4DFyQdyMiIvzySS/PVf9raPOGWdagl+bLyN5JpMt7MKSZyBwqexu81EWVjMy2L+td8TdC5XOcyByLQPh+tFrax4QQfxBC6PQ8hBDFQog/Ao+gteAeRiHOsUMIkQZ03wcr/yyP5ZAJB8L1Af3Vaar18CRoFAqFQnHoDHOM5JrcG7h79F/454T/8LcJT/H7MQ/yhewvkmEdnMKm02Lk1tMy+cKkZCz91DsbDXDB+CRum58V30bYUaCVgJjc+u2BFth/L3xyHWz7PrRv0+8XZhj3R82UIYIZSXNZmHq2bltLoJnf7fsZ+zv36LaHZJAnKh9iW7u+FnhywimMTdCqEg1CcMOstJjj/9uq2gGzoav2t1MaZaV93rgkzEaDpnOdGKXFW/o3aNvR75w9+Fs1Sb1I0haQlTCRU6PUMtY0r2B3xyDnDdMaaGF961rdtkVZc1k0Ut/Euba0gwONh2bTsGq/PhucmWCmKDPqJn/ZIxA8cW0gjnogLKV8DfgtYAZuA/YIITqFEKVCiE5gD/BtwALcLaV8Pf5sxxXfRQvwW9DsnE88wjXC7X79B2uO8/CL6xUKhUJx+PSX+R34sYJzxiZxz4X5XDQhifwkCw6zNl+u28yikQn87KxcLpqQPHBzoGs0THpo8KYaBjuM/Z3mNtgHV2Z/iREOvaFGk7+BX++5k6crH2F7+2Y2tKzld+FGwUiSTCl8Of+bum3ZiRYunagPuNu8Ie5dXkObV9/o1U1Fi4+norLB2YlmFnYHksIQbhiMODfSD9t/MLhygPLHtDraSAq/AcAlWVfHlIg8Vv5/+IdQevFW3SsEI6T0zMLCnOT5nDcuCbtZ/755Yn0DoSHqNrd5gzGB8JwiJ6LoJv1Aszu2TOYE4pgYakgp7wAuAFajZX1taGUQtvD6SuD8sJTaMUUIcZsQ4oAQ4umo7eOEEH8VQoyL2m4TQtwJ3B7e9Fsp5RCLho4TDGY8CII+/QdbceLnQyRboVAoFOCyGrlgfDI/X5LLA18o5MHLi/jlOXl88ZQ08pKGoDObMBpOeUFrDOwPazZMCddHx8FssHBb8Z0U2vWyXEGCvF3/Kr/fdxcPlP6OXR36LLPAwE0Ft5EYnZ0GzihJjCkHqWr18+ullTESck2dAe77oAaPXx8cXjElBVNkeUjCWMi/Xv9Enftgx080++x4tO+OzQannwkJmkJTosnNRVlX6HZXest5vvqJ+HNG0Ohr4J16fd5wRtIcXKYEXFYj54zVn5+99V6W7RmaOsW7u1rxRQhTGwXMKU4A9xRImgEJ4zQZwVOe1+4anKAck0AYQEr5qpRyHuAGpgCnhn+6pZSnhTPHh4UQYq4Qor57QXOtA7gjcrsQIj/iYUloOsbR95vMwDeALUKIWiHEOiHEOqAB+F+0c/cQcPfhHvdnhsHKPlIgqJfbKUk/cd/QCoVCoeifuCUQg3qwG8beA1Meh6xLNC1lYQYMkDARRt0F01/UguYBcBidfLf4xzGZ4XiYhImv5d9KiWtcn/sNQnDTnAxy3Xpd5pq2AD97o4KnNzSwobyD17c189M3KmJqX2cWOpmY00ciqPhb4Bqr31b3Juz5Xd9ue0GPVioSiiwXMEDRzbphC1PPpsCmd397u/5VNrVu6PP1RfJK7bP4I3JwRoycl3Fpz/riUYlkuPTtS89tbKSxQ/+a4+ENhGJ0lWcUukh1hucc9yeY+jSkLdIy5ycwR/3ohRDXhRcrgJSyQ0r5qZRyZfjnkVSJMKMZX3Qv3dXzjqjtgzH3PgD8BM08ox0oASagKUU8B5wtpfyqPFSPyOMBYWG3LIzaGCQ7bXifwxUKhUKhQAhwT4XRv9Dss+d/Ags2wbQnIecLQ9KgTzAlcvvwX3BuxiUI4gfoLmMC3xt2F7OST+13PrvZwLfnZ5Ls0H/N+4KSt3e28sCKWp77tIkOnz6bW5Ri4foZcdwPDWYY+1swRgXJFf8Ol0lESMn5W2HLrZrMXCS5V4FT/91qEia+VvDtGEOVvx38Iwf7UdHY2vYpKxrf1W07LfUMMq29/T4Wk4HrputfT1dA8n+ragn0YScezVs7WmLO0ZIxEVlms/uQGkKPR8TRjuOEEEFgn5Ry5ICDFd2W1GUAZWVl5OXlDfCIwyDQyb3L/pdP677as8lqrub/vjAIgXiFQqFQKI4g1d5KPmh8h9VNXfLbjwAAFlVJREFUH9ARbCcog2Rbc5ifupjZSfOHZEXc1Bngz8trKGseuHIx1WHkzsU5JDsGEIBqXKlpKEdZXGNOhYyzNAvt2jc0u+1InCUw9Ukw6pUtulla/ypPVerVKNymZL5T/CMK7PqMcaOvnrt2f0/nMGgRVn47+i+4o5U9gIfW1LEyqs73jFGJXD0tNe7LrGzxcdebFTqjlYk5dm6bP7hGzSNJeXk5+fk9N/Hzw34LR5RjIftVBxy6dofi6GGwUBPQy6QlWmo/o4NRKBQKxclMljWHy7Ov4/Ls6w57rmSHiR+ekc0znzSyYl9bnxUMAKfkO7lueiou6yBuFKfMhdG/gu13oLU3hfE3QMWTfT/GYIdx98QNgkGzyd7XuZu1zR/2bGsJNPGbvT/miuwvMS9lESZhYlfHNv5+8N4Ym+2Lsq7oMwgGuGpqCjtru6iPKIl4Z1crLquBC8bHPsYflDz6UX2M2+BFE/qe//PAsQiEPwQWCyFsYStlxfGCwUSzX2/kl6mk0xQKhULxOcBuNvDlGWmcNTqRN7e3cLDZR1WLH5fVQEmGnal5DqbmOYbmJph5HpgSYNsPIThA85k5BcbfD47ifocZhIGv5N1Ck7+BXR3be7Z7Q108XvEgz1Q9hs1goyUQ+/08JXE6Z6VdELO9G4fFyM3zMvjfpVUEIhxUXtrcTKcvxKWTkjWpOMAfDPHXD2vZU6+XQjurxE1RSvxA/kTnWATCvwTOBf4IfHOAsYpjjNevV4zIc5yYAhgKhUKhUPRFdqKF62dqLnHd5aBDCn6jSZ0P056GXb+A5rV9j3GNgfH3gi2n7/1RmA1mbim8nQdKf6sLhkELiL2h2DxihiWLG/K/NeBrKUyxct30VB5eW6/b/vbOVjZVephR6MRoEKza30ZNm77sI8Nl4sIJg7f+PhE5FoGwG/g18FMhxEzgCWA7/VgpSyk/OAbHddITkiGkX/8Gz07NjzNaoVAoFIoTm8MKgCNxFMLkh6BtO1S9oNUFB9rAXtBrxDHE53KZEvif4p/xeMXfWNm0rN+xJc5xfL3gOziMzn7HdTNvWAK+oOTf6/S6ydVtfl7Z0vedYKtJ8LXZ6VhNJ7YqxEAci0B4GVoxjUCTS5s8wHjJiWpZfILhCXZBSH+7Iyl9oD+PQqFQKBQKQNMFTjhyFghmg5mv5N3CzKRTebnmWfZ27tTtdxkTmZ96BhdlXolRDEYAq5dFIxMxCHhyfYOuBrgvrCbBdxdkMTzNNtSXcMJxLALOD9BVlSuOF1p9nUQr6Lktn/83vUKhUCgUxytCCMYnTGacaxKV3nIafLU0+hvItRUw3DESwxAD4EgWjEhkRJqNf6yui6uokZ9k4brpqSdFEAzHIBCWUi442s+hODRavLEWkW6bcpVTKBQKheKzRghBri2fXNuRLVnMS7Lw07Ny2FHTxcdlHexv8GI3C5wWI7OLXUzNc2D4nGgEDwZVgnAS0+LtIvot4DKb+x6sUCgUCoXic4HRIBiXbWdctn3gwZ9zjlogLIRwAIuBbiONPcDSI+wkpzgMWn1edG8Bgw+j4eS5ClQoFAqFQnFyc1QCYSHEucAjaHbGkTQKIb4mpXzpaDyvYmi0+/y6dYNBSacpFAqFQqE4eTjimhhCiLHAc0Aa4AO2AtvCv6cCTwshJh7p51UMnQ6fPvA1mPxxRioUCoVCoVB8/jga4nD/A1iBpUCRlHKilHICUAy8C1iA7x6F51UMkQ6/XjjbZAzEGalQKBQKhULx+eNoBMLzAS/wRSllTfdGKWU1cA1aZnj+UXhexRDp9AV16yZjMM5IhUKhUCgUis8fRyMQzgF2SynrondIKWuB3UDWUXhexRDxRClqm01K7lmhUCgUCsXJw9EIhG1A3359Gs1o5RGKzxhvVEmwxTSA1YxCoVAoFArF54jPt4G0ol+8USXBNrOSTlMoFAqFQnHycLR0hDOEENfF2wcghLgW6DPyklI+fpSOSxGBz68//SoQVigUCoVCcTJxtALhkWg6wv3xaJztElCB8DHAH9TfELCb1Q0ChUKhUCgUJw9HIxA+iBbMKo5zAgH9n99hNn5GR6JQKBQKhUJx7DnigbCUsuhIz6k4OgSC+j+/02L+jI5EoVAoFAqF4tij7oWfxISC+sDXZVaBsEKhUCgUipMHFQifxMigXsXOZVWqdgqFQqFQKE4eVCB8kuIP+SFk1W1zW2yf0dEoFAqFQqFQHHtUIHyS0uLtJPrP77aqQFihUCgUCsXJgwqET1JafJ6YbW6L/TM4EoVCoVAoFIrPBhUIn6S0ertitiVaVSCsUCgUCoXi5EEFwicpLT6vfoPwYTEqHWGFQqFQKBQnDyoQPklp8+oDYWH0fUZHolAoFAqFQvHZ8LkKhIUQxUKIrwkh/iGE+FQIERBCSCHEjw9z3tlCiJeFEHVCCI8QYpsQ4idCiBO2u6zD79etG4z+OCMVCoVCoVAoPp8cDYvlz5Jvh5cjhhDiGuAxwAhUAGXAeOAXwPlCiAVSys4j+ZzHgg5fQLduVIGwQqFQKBSKk4zPVUYYqAdeBX4KLAGeP5zJhBBFwENoQfAPgHwp5VRgJLATmA787nCe47Oiw68PhE2m4Gd0JAqFQqFQKBSfDZ+rjLCU8leR60KIKw9zyu8DVuBtKeXvI56nVAjxFWAlcKMQ4pdSyprDfK5jiscX0q2bjaE4IxUKhUKhUCg+n3zeMsJHDCGEAC4Orz4UvV9KuQrYAZiBC4/hoR0RugJSt24xq0BYoVAoFArFycXnKiN8hCkAssO/r4wzZiUwGpgJ/H0okwsh8uLsyu3+paqqaihTDom6qlo6Gnr//ImmWsrLy4/a8ykUCoVCoVAMhag46KhovKpAOD4jwz+9QGWcMfuixg6FsoEGzJgx4xCmPXTuPabPplAoFAqFQjFo0oHSIz2pKo2IT3L4Z7OUUsYZ0xQ1VqFQKBQKhUJxgqAywvHp1gjuz2mi25XiULyJ8+Nst6CVW9QCdcDRknPIAj4O/z4dqD5Kz3Oyos7v0UWd36OLOr9HH3WOjy7q/B5djtX5NaJlggE2H40nUIFwfLrCPy39jLGGf3qGOrmUsr+C3H397DsiaL2APVQPcDyKIaLO79FFnd+jizq/Rx91jo8u6vweXY7x+T3i5RCRqNKI+HSXPSSJqL94BMlRYxUKhUKhUCgUJwgqEI7P7vBPK5ATZ8ywqLEKhUKhUCgUihMEFQjH5yC9NS9z44zp3r726B+OQqFQKBQKheJIogLhOISVIl4Mr94QvV8IMQetqc0PvHIMD02hUCgUCoVCcQQ46QNhIcRtQogDQoin+9j9ezTViDOFEN/vrhUWQhQCD4fH/FNKqbpRFQqFQqFQKE4wRHyJ3BMPIcRc4OWITS60Gt9O9MoOU6SUZeHH3AX8DFgupVzQx5zXAY+gXTRUoMmajUezVl4PzJdSdhzp16JQKBQKhUKhOLp83uTTzEBqH9sd4aWbQdv0SSkfF0LsAe4A5gBj0eTNngJ+K6Xs6u/xCoVCoVAoFIrjk89VRlihUCgUCoVCoRgsJ32NsEKhUCgUCoXi5EQFwgqFQqFQKBSKkxIVCCsUCoVCoVAoTkpUIKxQKBQKhUKhOClRgbBCoVAoFAqF4qREBcIKhUKhUCgUipMSFQgrFAqFQqFQKE5KVCCsUCgUCoVCoTgpUYGwQqFQKP6/vbuPtawq7zj+/cnMAJbCgG2lRnB4saEwEcW2ymDKaNEQBAUqaaBR6Qu1bUrQoq2lQzNSW0oxVgmYYgExMWKKWqkvxSBF0zDFGhRq21ARSsEBBErHKa8D8vSPvY9zOJx75tw7995zz93fT7Kzz1577TvrrHnmznP2WWttSeokE+GOSnJckq8keTjJo0m+meTMJMbETkhyZZLawbbbpNu5lCU5IMkZSf4mya1Jnm77bcMY1x6Z5JokDyZ5PMl/JDnXPt9uLv2bZOMYcX3IYr6PpSiN1yS5MMlNSbYk2Zbk3iSfSfLaHVxv/I4w1/41fseX5MQklya5Ocl9bf9uSbIpyVlJVo24dirjd8WkG6DFl+S9wPnt4Z3AI8DhwEXAMUlOqqpnJtW+ZeJ24IEZztm3o53VbrOS5FeBjwO7AJuBe4C1wHnACUnWV9Vj89nQKTWn/m3dA9w9wzn7Fl4HfKV9/QzwXeBR4KXAycDJSd5fVecOXmj8jmXO/dsyfnfs3cBRwJPAvcCtwE8DR7bbW5McU1Vb+i+a6vitKrcObTSB/AzwQ+DUvvLDgfuBAt496XZO6wZc2fbh6ZNuy7RuwAbg88C5wLHAp9s+3TDimjXAE2299wBpy18C3NaWXzzp97YUtjn278a2zsZJt38pb8AxNB+CfwfYu698FfDnbR8WcPzAdcbvwvav8Tt+H58OrAdWDpS/mia5LeCSgXNTHb8Tb4DbIv+FwxfboLx0yLnT2nMPDf4jcBu7f02EF65PRyVql7R1vjzk3Lr23DbghZN+P0ttG7N/TSTG68s9gRUjzn+p7cdrBsqN34XtX+N3fvr/lLYfNw+UT3X8Oh60Q5LsSfOJGuDyIVWuBrYCLwBGjmWTlookAU5qD58T11W1ieauxErgzYvYNHVMVW2tqqdHVLmu3f9Mr8D4Hd9c+lfz6rZ2//xewXKIXxPhbnkFzVdITwDfHDxZVU8B32gPX7WI7VqO3pLkc0n+Mcmn2omIe026UcvU/jRj2ABunKFOr9y43jmvTXJ1G9efTvIHSfaddKOmSG/S0ON9Zcbv/BnWv/2M351zZLvvzx+mPn6dLNctL233d4/4VH0n8Et9dTU3bxw4/hXgfUlOq6prJ9GgZawXq73JHcPcOVBXc/OLA8e/DGxM8rtVdeUE2jM12jtnp7SH/QmD8TsPRvRvP+N3lpLsQpPovgn4C5rJiX/UV2Xq49c7wt2yd7v/3xF1euf2HlFHM7sDOIdm8uGewI8DbwC+TtOnn0vyc5Nr3rLUi9Ut1Q5KG8K43jn30UxG+nmaoVPPp5lZ/g/A7sAVSU6YXPOmwhk038ptAz7UV278zo+Z+heM31lL8s4kBTxNM0nuEuB64NVV9S99Vac+fr0j3C29r422jajzZLvffYHbsixV1Z8OKb4uydeAfwJ+AbiA5q675odxvcCq6tIhxZuSvBH4DM0Ywb9K8oUR/xl2VpIjgA+3hxuq6o6+08bvTtpB/xq/c7OZ5s76SprVH15IM3fo1CR/UlU/bOtNffx6R7hbnmj3My6IDeza7mcaY6U5qKptNMtVAaxPsiQ/GU8p43pC2qThve3hQcDLJticJSnJAcAXaBKGTwIfGKhi/O6EMfp3RsbvzKrq6qp6TVW9qqr2pVk+7S6abzwv7qs69fFrItwt43w9Mc7wCc3NP7f75wEHTrIhy0wvVle34wSHMa4XSFV9B3i4PTx4km1ZatqJWNfRjLH8Is2yioN3HI3fORqzf0cyfsdTVV8HjqO5u/tbSV7Snpr6+DUR7pbb2/3+SWYaFnPgQF3Nn6f6Xjssaf70YnVX4EUz1DGuF1Yvto3rVpJ9aJK0g4CvAae0K/MMMn7nYBb9Ow7jdwxVdS9wC03ueHhbPPXxayLcLd+i+Qe/G3DE4MkkK2kmE0AzuUvz67C+19+bWCuWn7tpnooIzQSYYXrlxvU8S/ITwE+1h8Y1kGQPmoc7rKVZkvKEqprpa2Hjd5Zm2b87+lnG7+ysGNhPffyaCHdIVW1l+3Paf2NIlVNoVjr4H+Cri9SsLjm73d9WVZsn2pJlpP0q9O/aw+fEdZJ1wCE0HwL/fhGb1hW/DwT4AdvXIe+sJLsC19CsmfrvwLFV9X8z1Td+Z2e2/TsG43dMSdaw/U7wrbA84tdEuHv+jOZxh7+Z5NReYZLDgQ+2h3/ZTu7SLCR5fZLz28kb/eV7JbkI6PX3eYvfumXvQppZy29I8p7eWLV2HNsVbZ3Lqur+mX6AhktyWJKPJDlsoHy3JOcAf9gWXdD13xvtmqufAl5Hs5Ti66vq4dFXAcbvWObSv8bv+JK8Msn7kjxnDkuSY2mWm1sBfGlgZY6pjt+4Ukj3JPlj4P3t4Z3AIzRfMT2PZsLBm/uWRtGYkpzI9k/Gm2kWF18JHEozo7aA86pq40QaOCWSHEVzx6dnD5rxZ4/x7FnHr6iqe/quexvwMZo43gw8QBPXK4GbgaOr6tGFbf3SN9v+TfJymmFVAA/SfBUK8LNsf9Tq5cAZXV96qr258Mn28HaaGBzmvqo6pb/A+N2xufSv8Tu+JOuBG9rD+2mGiqyieXrc6rb8G8BxVfXQwLVTG78mwh2V5HjgXcAraQL1dpogvtgkeG6S7Ae8g+YxlAcDP0nzldt9NGsIf6SdeasRBn4Zj3JAVd01cO06mqcerQN+jGa5n6to7vY8MfgDumi2/ZtkNfB7NH16CE1cr6L5j+4mmjs9X16g5k6VJKfT/B7dkf+uqjVDrjd+R5hL/xq/42uX9XwrzTr3h9GsHbyKZrjkLcDfAp+Y6cm00xq/JsKSJEnqJMcIS5IkqZNMhCVJktRJJsKSJEnqJBNhSZIkdZKJsCRJkjrJRFiSJEmdZCIsSZKkTjIRlqSOSLImSSW5a9JtkaSlwERYkqZIm8jOdvvqpNstSUvRikk3QJI0KzcOKdsLWDvi/Lfb/VPAfwKbF6BdkjR1fMSyJE25JOuBGwCqKpNtjSRND4dGSJIkqZNMhCWpI0ZNluuNJ25fn5RkU5JHknw/yceT7NtX99eS3Jzk0SQPJPnrJHuN+HNfnOSiJN9J8niSLUluSPKWBXmjkjQmE2FJ0o8kORP4LLAf8F2a8cdvA65PsluSDwNXAKuB/wL2Bt4BXJPkOcMykhwN/BtwJvBi4HZgK7AeuDrJBxb6PUnSTEyEJUn9zgdOq6r9qurlwME0CfGhwFXA24FjquqgqloLHAE8DBwNHNv/g5K8iCap3hM4B9i7ql5WVfsDR9FM2js7yfGL89Yk6dlMhCVJ/S6rqqt6B1X1PeDC9vBEYGNVXd93/tvAR9vDZyXCwNnAPsCHqur8qnqy77pNwG+3h++a37cgSeMxEZYk9bt8SNktfa+vGHL+W+3+wIHyk9v9ZTP8WdcC24B1SVzOU9Ki8xePJKnfHUPKHuztq2rriPN79AqS7AGsaQ8/OmT4cL/dgBcA359VSyVpJ5kIS5J+pKoeG1bc7oed6z/fn+32ryJx1Bh/9O5j1JGkeWUiLElaCI/0vV5VVU9NrCWSNAPHCEuS5l1V/QC4tz08bJJtkaSZmAhLkhbKZ9v9OyfaCkmagYmwJGmhXECzxvDbk3wwyer+k0n2SfLrSTZMpnmSus5EWJK0INo1iN8EPESzVvADSf41yU1J7mjLLwfWTrCZkjrMyXKSpAVTVTcmORQ4CzgeOAjYheapctcCn2f7EApJWlSpqh3XkiRJkpYZh0ZIkiSpk0yEJUmS1EkmwpIkSeokE2FJkiR1komwJEmSOslEWJIkSZ1kIixJkqROMhGWJElSJ5kIS5IkqZNMhCVJktRJJsKSJEnqJBNhSZIkdZKJsCRJkjrJRFiSJEmdZCIsSZKkTjIRliRJUieZCEuSJKmTTIQlSZLUSSbCkiRJ6iQTYUmSJHWSibAkSZI6yURYkiRJnWQiLEmSpE76f+jEw/wt1CAEAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# Initial concentration of mRNA and Protein for each repressor\n", "mtet, mlac, mgamma, tet, lac, gamma = initial = [1, 1, 1, 2, 1, 1]\n", "# Non-dimensionalized production rate\n", "alpha = 15\n", "# Degradation Rate\n", "beta = 2000\n", "# Repressor/Promoter Leak\n", "leak = 1\n", "# Hill Coefficient\n", "n = 6\n", "\n", "mtet, mlac, mgamma, tet, lac, gamma = initial = [1, 1, 1, 2, 1, 1]\n", "\n", "t = np.arange(0, 30.1, 0.1)\n", "ymtet, ymlac, ymgamma, ytet, ylac, ygamma = list(zip(*odeint(repressilator, initial, t)))\n", "\n", "plt.figure(figsize=(3.5, 6.5),dpi=200)\n", "gs = gridspec.GridSpec(8, 1, height_ratios=[1, 2.5, 0.1, 1, 1, 1, 1, 1])\n", "ax = plt.subplot(gs[0])\n", "dnaplotlib.plot_sbol_designs([ax], [[plac, rbs1, tetr, term1, pgamma, rbs2, laci, term2, ptet, rbs3, gammaimg, term3]],[[lac_repress, gamma_repress, tet_repress]])\n", "ax.set_ylim([-10, 31])\n", "\n", "# Plot of repressilator dynamics\n", "ax = plt.subplot(gs[1])\n", "plt.plot(t, ytet, color=[1.00, 0.75, 0.17])\n", "plt.plot(t, ylac, color=[0.38, 0.82, 0.32])\n", "plt.plot(t, ygamma, color=[0.38, 0.65, 0.87])\n", "\n", "plt.ylim([1,4])\n", "ax.tick_params(axis='both', labelsize=8, width=0.8, length=3)\n", "ax.yaxis.tick_left()\n", "ax.xaxis.tick_bottom()\n", "ax.set_xlabel('Time', fontsize=8, labelpad=1)\n", "ax.set_ylabel('Protein Concentration', fontsize=8, labelpad=2)\n", "plt.legend(['tetR', 'lacI', 'gamma'], frameon=False, fontsize=8, labelspacing=0.15, loc=(0.06,0.65))\n", "plt.subplots_adjust(hspace=0.4, left=0.12, right=0.95, top=0.99, bottom=0.01)" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "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.7.0" } }, "nbformat": 4, "nbformat_minor": 2 }