{
"cells": [
{
"cell_type": "markdown",
"id": "49bcb5b0-f19d-4b96-a5f1-e0ae30f66d8f",
"metadata": {},
"source": [
"## Coupled pair of reactions: `A <-> B` , and `A` + `E` <-> `B` + `E`\n",
"A direct reaction and same reaction, catalyzed\n",
"### Enzyme `E` initially zero, and then suddenly added mid-reaction\n",
"\n",
"LAST REVISED: Nov. 4, 2023"
]
},
{
"cell_type": "code",
"execution_count": 1,
"id": "cbb1af2e-3564-460e-a4ae-41e4ec4f719f",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Added 'D:\\Docs\\- MY CODE\\BioSimulations\\life123-Win7' to sys.path\n"
]
}
],
"source": [
"import set_path # Importing this module will add the project's home directory to sys.path"
]
},
{
"cell_type": "code",
"execution_count": 2,
"id": "087c0d08",
"metadata": {
"tags": []
},
"outputs": [],
"source": [
"from src.modules.chemicals.chem_data import ChemData\n",
"from src.modules.reactions.reaction_dynamics import ReactionDynamics"
]
},
{
"cell_type": "code",
"execution_count": 3,
"id": "23c15e66-52e4-495b-aa3d-ecddd8d16942",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Number of reactions: 2 (at temp. 25 C)\n",
"0: A <-> B (kF = 1 / kR = 0.2 / delta_G = -3,989.7 / K = 5) | 1st order in all reactants & products\n",
"1: A + E <-> B + E (kF = 10 / kR = 2 / delta_G = -3,989.7 / K = 5) | Enzyme: E | 1st order in all reactants & products\n",
"Set of chemicals involved in the above reactions (not counting enzymes): {'B', 'A'}\n",
"Set of enzymes involved in the above reactions: {'E'}\n"
]
}
],
"source": [
"# Initialize the system\n",
"chem_data = ChemData(names=[\"A\", \"B\", \"E\"])\n",
"\n",
"# Reaction A <-> B , with 1st-order kinetics, favorable thermodynamics in the forward direction, \n",
"# and a forward rate that is much slower than it would be with the enzyme - as seen in the next reaction, below\n",
"chem_data.add_reaction(reactants=\"A\", products=\"B\",\n",
" forward_rate=1., delta_G=-3989.73)\n",
"\n",
"# Reaction A + E <-> B + E , with 1st-order kinetics, and a forward rate that is faster than it was without the enzyme\n",
"# Thermodynamically, there's no change from the reaction without the enzyme\n",
"chem_data.add_reaction(reactants=[\"A\", \"E\"], products=[\"B\", \"E\"],\n",
" forward_rate=10., delta_G=-3989.73)\n",
"\n",
"chem_data.describe_reactions() # Notice how the enzyme `E` is noted in the printout below"
]
},
{
"cell_type": "markdown",
"id": "0e771dda-1c0f-4fc0-ab21-049740643897",
"metadata": {},
"source": [
"### Set the initial concentrations of all the chemicals - starting with no enzyme"
]
},
{
"cell_type": "code",
"execution_count": 4,
"id": "5563e467-a637-44fa-9ba1-d35ddd82c887",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"SYSTEM STATE at Time t = 0:\n",
"3 species:\n",
" Species 0 (A). Conc: 20.0\n",
" Species 1 (B). Conc: 0.0\n",
" Species 2 (E). Conc: 0.0\n",
"Set of chemicals involved in reactions (not counting enzymes): {'B', 'A'}\n",
"Set of enzymes involved in reactions: {'E'}\n"
]
}
],
"source": [
"dynamics = ReactionDynamics(chem_data=chem_data)\n",
"dynamics.set_conc(conc={\"A\": 20., \"B\": 0., \"E\": 0.},\n",
" snapshot=True) # Initially, no enzyme `E`\n",
"dynamics.describe_state()"
]
},
{
"cell_type": "markdown",
"id": "651941bb-7098-4065-a598-e50c0b641ab3",
"metadata": {
"tags": []
},
"source": [
"### Advance the reactions (for now without enzyme), but stopping well before equilibrium"
]
},
{
"cell_type": "code",
"execution_count": 5,
"id": "76f24d9a-a788-41d8-90a4-db87386f91aa",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Some steps were backtracked and re-done, to prevent negative concentrations or excessively large concentration changes\n",
"9 total step(s) taken\n"
]
}
],
"source": [
"dynamics.set_diagnostics() # To save diagnostic information about the call to single_compartment_react()\n",
"\n",
"# All of these settings are currently close to the default values... but subject to change; set for repeatability\n",
"dynamics.set_thresholds(norm=\"norm_A\", low=0.5, high=0.8, abort=1.44)\n",
"dynamics.set_thresholds(norm=\"norm_B\", low=0.08, high=0.5, abort=1.5)\n",
"dynamics.set_step_factors(upshift=1.2, downshift=0.5, abort=0.4)\n",
"dynamics.set_error_step_factor(0.25)\n",
"\n",
"# Perform the reactions\n",
"dynamics.single_compartment_react(reaction_duration=0.25,\n",
" initial_step=0.05, variable_steps=True, explain_variable_steps=False)"
]
},
{
"cell_type": "code",
"execution_count": 6,
"id": "e58db01b-b932-4f60-91c2-a578353a3702",
"metadata": {},
"outputs": [],
"source": [
"#dynamics.explain_time_advance()"
]
},
{
"cell_type": "code",
"execution_count": 7,
"id": "4a19ad2a-fbd2-420a-b958-2daf88bcc841",
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
" \n",
" "
]
},
"metadata": {},
"output_type": "display_data"
},
{
"data": {
"application/vnd.plotly.v1+json": {
"config": {
"plotlyServerURL": "https://plot.ly"
},
"data": [
{
"hovertemplate": "Chemical=A
SYSTEM TIME=%{x}
concentration=%{y}
SYSTEM TIME=%{x}
concentration=%{y}
SYSTEM TIME=%{x}
concentration=%{y}
Changes in concentration for `A <-> B` and `A + E <-> B + E` (time steps shown in dashed lines)"
},
"xaxis": {
"anchor": "y",
"autorange": true,
"domain": [
0,
1
],
"range": [
-0.00011732081911262797,
0.2751173208191126
],
"title": {
"text": "SYSTEM TIME"
},
"type": "linear"
},
"yaxis": {
"anchor": "x",
"autorange": true,
"domain": [
0,
1
],
"range": [
-1.1111111111111112,
21.11111111111111
],
"title": {
"text": "concentration"
},
"type": "linear"
}
}
},
"image/png": "iVBORw0KGgoAAAANSUhEUgAABUwAAAFoCAYAAACMirH9AAAAAXNSR0IArs4c6QAAIABJREFUeF7svQuUHdV5qPl3q59CDws8IMvhYeQ4PMzDELAcXxMS4rm2bEIuBgVj35goUQR4lsdAxJLwMArjZcSShsfyJWCNbhTy0DWRHZwMQY5fWQrOxGBsLjaOITeBGDsXGxmQQUj97p61d9U+Xae66pyqU7uq9j7n67W01H1O1X58/z6n63z917/75ubm5oQvCEAAAhCAAAQgAAEIQAACEIAABCAAAQhAAAIQkD6EKasAAhCAAAQgAAEIQAACEIAABCAAAQhAAAIQgEBAAGHKSoAABCAAAQhAAAIQgAAEIAABCEAAAhCAAAQgEBJAmLIUIAABCEAAAhCAAAQgAAEIQAACEIAABCAAAQggTFkDEIAABCAAAQhAAAIQgAAEIAABCEAAAhCAAASaCZBhyoqAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEIBASQJiyFCAAAQhAAAIQgAAEIAABCEAAAhCAAAQgAAEIIExZAxCAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEmgmQYcqKgAAEIAABCEAAAhCAAAQgAAEIQAACEIAABCAQEkCYshQgAAEIQAACEIAABCAAAQhAAAIQgAAEIAABCCBMWQMQgAAEIAABCEAAAhCAAAQgAAEIQAACEIAABJoJkGHKioAABCAAAQhAAAIQgAAEIAABCEAAAhCAAAQgEBJAmLIUIAABCEAAAhCAAAQgAAEIQAACEIAABCAAAQggTFkDEIAABCAAAQhAAAIQgAAEIAABCEAAAhCAAASaCZBhyoqAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEIBASQJiyFCAAAQhAAAIQgAAEIAABCEAAAhCAAAQgAAEIIExZAxCAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEmgmQYcqKgAAEIAABCEAAAhCAAAQgAAEIQAACEIAABCAQEkCYshQgAAEIQAACEIAABCAAAQhAAAIQgAAEIAABCCBMWQMQgAAEIAABCEAAAhCAAAQgAAEIQAACEIAABJoJkGHKioAABCAAAQhAAAIQgAAEIAABCEAAAhCAAAQgEBJAmLIUIAABCEAAAhCAAAQgAAEIQAACEIAABCAAAQggTFkDEIAABCAAAQhAAAIQgAAEIAABCEAAAhCAAASaCZBhyoqAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEIBASQJiyFCAAAQhAAAIQgAAEIAABCEAAAhCAAAQgAAEIIExZAxCAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEmgmQYcqKgAAEIAABCEAAAhCAAAQgAAEIQAACEIAABCAQEkCYshQgAAEIQAACEIAABCAAAQhAAAIQgAAEIAABCCBMWQMQgAAEIAABCEAAAhCAAAQgAAEIQAACEIAABJoJkGHKioAABCAAAQhAAAIQgAAEIAABCEAAAhCAAAQgEBJAmLIUIAABCEAAAhCAAAQgAAEIQAACEIAABCAAAQggTFkDEIAABCAAAQhAAAIQgAAEIAABCEAAAhCAAASaCZBhyoqAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEIBASQJiyFCAAAQhAAAIQgAAEIAABCEAAAhCAAAQgAAEIIExZAxCAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEmgmQYcqKgAAEIAABCEAAAhCAAAQgAAEIQAACEIAABCAQEkCYshQgAAEIQAACEIAABCAAAQhAAAIQgAAEIAABCCBMWQMQgAAEIAABCEAAAhCAAAQgAAEIQAACEIAABJoJkGEaWxGPPfG0XPXx2+S+uzbLeWef0vTsHTv3yr6/e1R2br9BVp+4qvHc2PikbN2xW/98y6b18r2nn9VtRL8+eeN6Oev0N8vGG2+XH7/wUtt1qPp/85veKNdsvlPOP/sUuX7jugXnPPPc87q9az9yiVy69oK2bXIABCAAAQhAAAIQgAAEIAABCEAAAhCAAAQg0JoAwjTG5+ArhxIlpXn8yaeeFSU/o4IySVwaifrD5w/IvbddJyuWL02MhJKwPznwshatoyNDTcekjcUchDDl5Q0BCEAAAhCAAAQgAAEIQAACEIAABCAAAbsEEKYJPJXE/OYTTzeJTiUnP7Ftlz76hDce1yQ4VVbqlm27mjJPEaZ2FyqtQQACEIAABCAAAQhAAAIQgAAEIAABCECgCgII0wTKD+x7WO75k79uEqDqsUe+/X15/7vfoZ+LZo0mCVaXhKnJRE0qBfCG445pmmc0k9agiZcniPL567/9B/mjz+7Th0aPU0zM4+q53/ng2sSyAmmL3JRGMM+fcerJTcyj2bfvevuZTSUQ4hnA8bFE+1Tj+u0r3pta+iCe5WtYbtuyQb7+6Heb5qjm/9ZTTtblGR762iONbpLKO5j1ET0uPu4q3gDoAwIQgAAEIAABCEAAAhCAAAQgAAEIQKCZAMI0YUUk3equpJv6uuQ9/0HXDVXCTNU4NeJr5bFHNwlBl4Rp0qI3c1z7q29vjNtIyqi4S2KhhOnN24OarXEZGK/nqsoMGOl4wqpjE0sPxMeXVCs2LrGjYjcqY1vVoI32E+8jSZKr4+PZw1H5nCaI448nZSurNRRl3678Am9cEIAABCAAAQhAAAIQgAAEIAABCEAAAtUQQJgmcI5LUCOzbti4rpFBuObc03Qd07Q6ojaFqaqb2uorb2aimY9q02TKpolfdUw8gzZNLqpj057LWm81msEZ3XQrLhTTBGOreRiGRvhGuZn21l18YaM+bZL8TRtfnsfT6ta24lrN2wG9QAACEIAABCAAAQhAAAIQgAAEIAABCCBMU9ZAVBK+/LNDsuOe+2XbTRv05k1R4fW9p59dUL9UNWlTmJ5/9imJt7NnlZDRKZpxPf69f2m6FT9N+Klz41mWrcRemgzMmkGp2t774P7EjbKibY9PTCTeRp8kOaPzNxmoSSUC4mI4iUkeMar6jR+fJGbN+FrFgLcqCEAAAhCAAAQgAAEIQAACEIAABCAAgWoIIExTOEcl4Xf+6V91/VKzk7167vade7XU++P7v5i4y72rwtRkV8ZvpY/XDI1jidY6TROmrWRlO5Fp+mtVb1Qd876L1ug4dCJM25UGiAvLJPlbVJi2qidrGCTVPK3m7YBeIAABCEAAAhCAAAQgAAEIQAACEIAABBCmKWsgmr35gx/9RE46fmXjVm3z3P953W/J33zlGxKvX6qadFGYtsquzJPdWFeGaTRU7W7JV8cawR2NRzyzNh5+I0mv+73L5bo/+ENRZRiipQGKCtNWGaa8HUEAAhCAAAQgAAEIQAACEIAABCAAAQjUTwBhmhKDhvD8ny/IiwdfbWzyFJVvS5Yslocf+U7Tc6Y514SpEX3nvPXnEzdeynrLvJpfmTVMs27alFeYmszVdtmbJrP4gjVnyff/xw8WlAYoKkyzZtrW/9bACCAAAQhAAAIQgAAEIAABCEAAAhCAQG8SQJi2iLu5ff2MU09eIM5aPReVqj98/kBiPU7TbVrNT/V8O4mZtYZpVnmbloGqxrHt03tky8c+pGu4thKmSX21uxU+HoI0uanG9/VHv6vrueYRpmZeWTbHMuN/6GuPSNLxRYWpmmuavFZ977j3fvnQpb8mq09c1chSbpcV25tvXcwaAhCAAAQgAAEIQAACEIAABCAAAQiUQwBh2oKrEW2mbuboyFDj6CzPRZtOk3VVCNN29UmjWZdGRD751LNNZKLjz7Kbe7wWadImS62WdNKYo3VUswpT1cfWHbtFCdCkr6TYttp4yoYwjQr1+LiinNI26CrnrYBWIQABCEAAAhCAAAQgAAEIQAACEIAABBQBhCnrAAIRAu2yeoEFAQhAAAIQgAAEIAABCEAAAhCAAAQg0N0EEKbdHV9ml5NAluzZnE1yOAQgAAEIQAACEIAABCAAAQhAAAIQgIBHBBCmHgWLoZZLgOzScvnSOgQgAAEIQAACEIAABCAAAQhAAAIQ8IEAwtSHKDFGCEAAAhCAAAQgAAEIQAACEIAABCAAAQhAoBICCNNKMNMJBCAAAQhAAAIQgAAEIAABCEAAAhCAAAQg4AMBhKkPUWKMEIAABCAAAQhAAAIQgAAEIAABCEAAAhCAQCUEEKaVYKYTCEAAAhCAAAQgAAEIQAACEIAABCAAAQhAwAcCCFMfosQYIQABCEAAAhCAAAQgAAEIQAACEIAABCAAgUoIIEwrwUwnEICAywSeee552Xjj7XLOW39ebtm0XkZHhlwerrNje2Dfw3Lz9t3yyRvXy6VrL3B2nHkHxvrIS4zjOyXQra+hTnlwHgQgAAEIQAACEIAABOoigDC1QP6OnXvlm088Lffedp2sWL7UQos0AYHeJaBeT3/02X0awPsuWlOZwDRSbNuWDXLe2af0bgAKzlwJn70P7nf6/bCTWHdyTkGUnN6jBGy8hox4/Z0PrpXrN67rUZJ+TVvF7J4/+WvZuf0GWX3iKr8GX/Noy35/rvI6X/X1kwMvt732MXO+9iOXOPcHyseeeFqu+vhtlV7D1bwE6R4CEIAABLqUAMK0RWAPvnJIrtl8pzz51LNNR8U/gFR5IeX6OjQf0u67azPSqcNglf2hKa19s95PWHVs2wv1DqfW9rS6X0t5RIURu2ecerLTcrAt9JQDzGs56elWGaRj45OydcduWXns0c6Kmk4/XGddH9HfHXULK/PBNW0dVPlHiU7XYpbz0n5fq3PrfI22eh2psSWtj6KvoTSR0up3i3o/2/d3j3on68r+fZll7dk6ppvmophUuaY6fU/PGrsqr02KCFOXrsFNTNb+6tudvRbIGn+OgwAEIACB3iWAMG0jC+IfZsyHMnWaySit8kLK9aXq0sWa66zSxlf2hyZXhal5ba27+ELnsiXisTJjXb70KPmHbz4pNv9AYGTJD58/kEvEJr03FXkNpK0TI+B8vu2+7A/XitGWbbvk9SuWifT15YpjUsyyitqkc028bK7ReD+drglzng2pbNb/+Wef0vThvNPXk5ljUXlZ9vt5PBZmvOrxeHkRhGmRd8Tyz616rZQ9I4RpZ4S7RZiq2Xfbmu4sopwFAQhAAAI+E0CYJkQvy4e/PQ98RdZetEbfgo8w9fkl4N7Yy77ALLv9Ton6JEyNEFO379++c6/EJU0nDIzoeOhrj3ScEWejDTN2V8V6J2zj55QtTM0H3ve/+x36LoWistJVYRrNXu1UoJtY/PiFlxKzLbPGO02Ymg/tqrZuJ3HwTZi2EuSuvvdnjXHaHxO65Tb2boxPkdjmObeK9/SqSm8VEaZ5mFVxbKv35Sr6pw8IQAACEIBAUQII0xjBTrJRjDC9dcsGuWnbrsYt/PFbHaMyI9pt0i2RWduMyg31gdB8nXHKm+TFg69KvLZR0hiSPuhG60iqNrPe0hi/4I9eLL3r7WfqmkbmK+8H7KRbG+MfgOPjTssQVoIr63iy9Bs/Jh7T6C2Sav7RWEXnkHb7pjnGrIs7/+Cjcuf/8zlRcu0Nxx2jb6P8zj/9a1O7qh/znKmH1qr9N7/pjVruJMm/MrhGXwNRaRJ9PMomzxguec9/0Js4KQlj+5bjeAbXvX/yV4VqGNuUnIZdtM1O599OmNqQxNFYZy2B0sl7yoL3s/D9sYx6tdHX+q+8822pr6k8v7xdE6Y2RGmaxFav2by/G1Rb7YRpp1LNN2Ga9gfcrL9bzJ0z0T9gxX9nmfjEX1dJcct6zZH2emh1LdJuTqrNeEmKpGuZPNdbSSUusor4duea99y7/q//Tf5075f07/dW11/x35vx3/dJpRnSMpDj7/d5x5IUv/ha7OS9O/O6SHhPT1ofcUam/XbXvHnWSNZ1p45LiqG6M+GENx7XtjRSUnxtXIO3u55NGrd6LO19u8jvrzy/JzkWAhCAAAQgUAYBhGmMaic1d5JqGSZ9eFMXqjvuvV8+dOmvNQr6p128Zm1TDT/ptqekC6mkuSWNM36R2+oWv/iiTLtYU3Vgo/Iyz22iaRJbzedrX/+2/N6HL5akMSbV5EyrLZg0njz9xm+fjjOMXhRHLyqTpFS72ybVhkhJF/3qPPUV3Z08T/tpa1bVo1Rf5vbOolzT3sjSMkw7iW2nkjDLm2z8tZVnLUfbL0OUxsdfRJymrcOysqBU/Ld9eo9s+diHGpvnJa1JG6/hTmOWZX3E+di4A6HIB06bcy1DlMaZpr1XtmOfJkyLZjj5JEzbzbXd75ZoBl30dRb/nWX+4Nfud1nWa45WUiw6pqTfBe3mFK/LmnR80vVWUl9Jr6Ws74dZzjWiKvr7Pe36K629aCZ10to1MVHMo5tLxbMa84wla/zyvHentZl2XZbGI8s1UZZr3qLX5EnrJOlutjzX23mEaZZr8DS2adez0YSMVr9nbP4Oavd7gOchAAEIQAACtgkgTFOEaZ5dJ1tldGTZLdrc3hu/eE26/Sf+wTntNqQk+ZR2m0/0Qu7o1y1NzYj60v5vygVrzpbRkaHUddjqr9vRnXrzfAjN8oEk7Zj4BWW7OnfRjWqy9JsUOwUn3m/aJhxJcWr3ATDPbWF52k9iUwbXvMK06Bhsv2nGx5PnA44Zi/kwmjVzu+gcouI0a73IpOwcM46qxq36i7/n5XkNt3t92s4wTVoLNnYytiFM09ZQlj8upMmzouuy1flRcZolgy8tQ9n00UnWqjo3z++qpPm0eh2p47O+HrOwbrfW8vxuSfsDVp7Hs1xzpO0G30r+Rq9F2v1+iL/Gs/yR2LCOi560a73v/fO/yejIcMud7bOcmzaX+PtYKzbxfuLvHernL//9t/QU/9df/kX9B9Y81wlp76lJ67NVhmmn14NF39Pjc826zrJe57e7Jjd3ZrR6X0l73cQZ5xGmabWdo9e87diaz0WtXnNj4xPy1l94U9NQ094zsryncQwEIAABCECgbgII05KFadJtgEkfoOIZg60uzqJtpl24pF0UJm2oE73Ae+spJ+sdrtWtYJ18wMwrTBX++MYU8RdFlovHtGPiF8PtZEt0PEX6jV8MtxOm0YvZPB9q0y6g1W2t0a92mUDq2LQPkj858PKCGBXhmvaml3ZRXTS2Nt9ks97O2K7PToRp0vtGFtkVlT7qdZ1V0LRbhyrTOYvIasci6fWu2o5+RQVt3tdw0h8Yyqp3l+cW2HavA5UR1O4ry3u0jeyeToRpmsDMumY6FaZJpSJMW+e89efb/r5Jum06KQ5pt/bGj83yh7d2cc76fLt13e41nZRhGr9myCpMW0mSduOMv2e1Wudpc2r1R4b475S0663469m8B2d9343GLcu5rYSpKmdkXjut+MXbiEqwVce9Xl/frTn3ND20R779ff16eP6FF3UJm6hczjqWVmszrzDNej2Y5z09fsu7Ga9ZU9E/JrZaZ3muydMSJaLrzjBPSs7Icu2p5mFDmGa55o1fz5r3yKx/OEWYZn0H5zgIQAACEHCRAMI0FpVOb8lPywaNyk1z0fDiy6803QqVN8M02mbaRVz8AiXtojE6fXMxnlRzrNMPh3nkRtILJEv2Xqtj4s+NT0wkZtDGj1Njid+KHh9fEqf4MUZOVSFMzS1j0Yv+PJkj8ViVwbVVdnLSWG2MweYbb7vXURaBZcbj6y35avytsnI65W3Yvv7o5U27yufNMDUfAFu9hrMIm07m0SqbMOt7aFK/NjJMs4rKVvP28ZZ8M59OxaWNDNNO66fmXYPt1nWVwrTde6WaW7s1meVaJG1O5vdhGsOo9Mx6HaXaSnqNZ33fb3duVkmZ5/bn6O/Vs05/s3xi2y751JYNGov5/uWDr+rNC039WjPPpHWb5w8wtoVpq+uBpLWf9ZooyzrLKkyzrrvvPf2srumf9BqoQ5ia35embm7S6yaprFb0uLQ/xiJM876TczwEIAABCLhEAGEai0ZaDZ9WQct6IZX1divVV542k/6anSfDtN2CTBMZaR/skyRx2u1ARm60EmlZLh6LZiEmXYgX6TfOpmxh2q79aJZQu6zkaKzK4Jq23nzIME17XWYR+2nzLkOcRtvsJBuq1Qdm9ZyJ1Qmrjm2bsdfu/cU8n/WWxzx/hElbv+3EUtYxR49rJZHb3Srdrj9XhKkZZxniNCrYsgqorPyjoqudpIvHwkdhmlZSqEphWoYkSboW6STDNB7jrBmmaWtDSaZO11X03KzCNE+GaXT9nnT8ykZWqRFkKtv0Bz/6icTvJMk6ljzXxnneu9PazfqenueaKO16LfrHu6LX5Gl9JJWFyXLtqdqrKsO03e8n85lF3RmS9N5dxntBljFxDAQgAAEIQMAGAYRpAsWkQuzxwx5+5Dt65/gVy5dmlptpF1xFMkzT/tofv0DJKnTUeepW0AvWnNU05azZOWXckp/Wt5rTw488If/xwvN11kdSNkT8gjLPBXuRftMujuMfZJPG06pGWKs1lJStkHShmtZ+0ljK4Jr2xpV2UV10DDbeKFUb7bIqTfZQ3g/OZnw2xKmNNsx4Wr3m27HIy7zV+1OnGabqjzDt1o7NGqatMq86+UNclKFrwtSMzYY4jYrSrOUiktZXuzWZ9XeYz8K0HYM8v1uy3npveHV6zdHq90GWa5G0ORXJhIyu7y3bdjXuCEqq4571jyFZzs0qKVvFOekawbSrdl9f9+u/0tgYUj0er2fa7v2/CNc8119p6yLre3rWa+Os17xZhWlWPu1imFQKKct1ZZFr8KzvkdHPP/HXf1JJlKxM8l43cDwEIAABCECgCgII0xTKRn7EP8CZixx1mrl9qciFlLmQ6LSGaZJoSKvJlFbHTR2/49775UOX/pokbfqUVbYqJkUu1tIWfJJsiGf9JB2TlAWX54I9T7+Pf+9fmsosGBbqf7WpQrtsh+hFZqsPYO1uHYy2Y+aqPnRG/+qfZyxlcM0rTIuOwdYbabuL/nayIus4OpVrSe9NWftMOq7VHwxUuYqkNW+7P/M+3GkN06SYRF8XncrtpHmmvTbjcrGTPl0VpvG55c0MNa+pIqI0zwf2TvrxKcNUsWi1DvP8bikqTNVYslxz5Nn0KelapN2ckmotq3X39Ue/K2bjoSRmSeWZWsnI6KadWd8f4u+xWYWpaj/pD/tpf7Qzc1HnRceZ9njStVwn72O2b8lXY8j6nt7uOPNelXRc2h1HWUpvmddglnWXFG8Twyx3htjOMDXzbnc9mzTuVtdHRX5/Fbmm4FwIQAACEICADQII0xYUox+so4fFP3RlFabRi1zTnmrrXW8/U6JZDK0+9CRdqMTrLyn5qjKnVE2q+IYNSbWaVH/ROSXVHsv6QbMMYWpYxWtDJV1Qxo9JE955SgRk6ddc5EbXSVTy5JGUSevESJZWH4bj61X1f+uWDXLTtl0L1kF8YxPV/pvf9MbE+q7RDwDRdRvd5TaPiE57yaV9QE+Lf9bY2nijzCox20kzG2Opqo2kNW36zvJhrpNxxvtUH2rVV7TsSN61lvS6uGHjOv2eayvDNEmuxOdfRKgX+cDZbhOjsmLZSfyLnJP2+9q02YmoVufaEKY3b9+dOrWsv1uzsmn3h52k9/7zzj5lgWi1IUyj/OK1EbPMO+u1SNqckn6XqseS/kAd32hOHRf/A0DS9VPWjW+ynJtHmKrxxfmk1UlO+uNxNDY/fP5AU/1S9VzesSStzzKEqeon63t61muiLOssz3V+1nVnOEffH9SaSyqRkMTXtjA1fbS7nk26Jmy39pIyT7O+p3EcBCAAAQhAoE4CCNM66ZfYdxk1+kocLk1DAAIQgAAEIFCQQJ47Qgp21TWnd9MfuromKEykKwhkvc2/KybLJCAAAQhAoCsJIEy7IKzqYl9lqaosEfOVtWh8F0yfKUAAAhCAAAQgEBLIWlcTYAEBhCkrAQL2CfA+ZJ8pLUIAAhCAQPUEEKbVM7feY/yWcdVBltvdrA+EBiEAAQhAAAIQqJ1AWh322gfm4AAQpg4GhSF5TcCUyeiWki9eB4PBQwACEIBAIQII00L4OBkCEIAABCAAAQhAAAIQgAAEIAABCEAAAhDoJgII026KJnOBAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEChFAmBbCx8kQgAAEIAABCEAAAhCAAAQgAAEIQAACEIBANxFAmHZTNJkLBCAAAQhAAAIQgAAEIAABCEAAAhCAAAQgUIgAwrQQPk6GAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEuokAwrSboslcIAABCEAAAhCAAAQgAAEIQAACEIAABCAAgUIEEKaF8HEyBCAAAQhAAAIQgAAEIAABCEAAAhCAAAQg0E0EEKbdFE3mAgEIQAACEIAABCAAAQhAAAIQgAAEIAABCBQigDAthI+TIQABCEAAAhCAAAQgAAEIQAACEIAABCAAgW4igDDtpmgyFwhAAAIQgAAEIAABCEAAAhCAAAQgAAEIQKAQAYRpIXycDAEIQAACEIAABCAAAQhAAAIQgAAEIAABCHQTAYRpN0WTuUAAAhCAAAQgAAEIQAACEIAABCAAAQhAAAKFCCBMC+HjZAhAAAIQgAAEIAABCEAAAhCAAAQgAAEIQKCbCCBMuymazAUCEIAABCAAAQhAAAIQgAAEIAABCEAAAhAoRABhWggfJ0MAAhCAAAQgAAEIQAACEIAABCAAAQhAAALdRABh2k3RZC4QgAAEIAABCEAAAhCAAAQgAAEIQAACEIBAIQII00L4OBkCEIAABCAAAQhAAAIQgAAEIAABCEAAAhDoJgII026KJnOBAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEChFAmBbCx8kQgAAEIAABCEAAAhCAAAQgAAEIQAACEIBANxFAmHZTNJkLBCAAAQhAAAIQgAAEIAABCEAAAhCAAAQgUIgAwrQQPk6GAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEuokAwrSboslcIAABCEAAAhCAAAQgAAEIQAACEIAABCAAgUIEEKaF8HEyBCAAAQhAAAIQgAAWXt4SAAAgAElEQVQEIAABCEAAAhCAAAQg0E0EEKbdFE3mAgEIQAACEIAABCAAAQhAAAIQgAAEIAABCBQigDAthI+TIQABCEAAAhCAAAQgAAEIQAACEIAABCAAgW4igDC1EM1DY9Ny6MiUhZZoAgIQWHn0qBw4OCazc7CAAARsEDj2dSPy8qEJmZ7hRWWDJ21A4Jhlw/La2JRMTM0CAwIQsEBgxZIhGZ+akbGJGQut0QQEILD8qEF93Xd4fBoYHhBYdcyoB6PszSEiTC3EXclSJU35ggAEihNAmBZnSAsQiBJAmLIeIGCXAMLULk9agwDClDUAAbsEEKZ2eZbdGsK0bMKdt48w7Zxd40yEqQWINAGBkADClKUAAbsEEKZ2edIaBBCmrAEI2CWAMLXLk9YggDD1aw0gTN2NF8K0YGxuueUWuX7TTWSYFuTI6RAwBPbsvlsuu3K9DI8sBgoEIGCBwF/vvU8ufPf7ZfmK11tojSbKIvDvP/w3efKJb8l7f/3ysrqgXUsEvvLQ5+Wst50vx646wVKLNFMGgZdePCD7v7pPPnDFVWU0T5sWCfzj/i/KG48/SU5cfarFVmnKNoGxI4flL++/Tz68/qO2m6Y9ywQef3S/LFmyTN5y+jmWW6Y52wR23b1dtm7dartZ2rNEoOeF6dj4pGzdsVse+tojDaT33bVZzjv7lMbPD+x7WG7evlv//L6L1sgtm9bL6MiQ/hlhamkl0gwEQgIIU5YCBOwSQJja5VlWawjTssjabxdhap9pGS0iTMugWk6bCNNyuNpuFWFqm2h57SFMy2Nru2WEqW2idtvreWF68JVD8sf3f1Gu+chvaAn62BNPy5Ztu2Tn9htk9Ymr9M+379wr9952naxYvlTu2LlXR+D6jesQpnbXIq1BQBNAmLIQIGCXAMLULs+yWkOYlkXWfrsIU/tMy2gRYVoG1XLaRJiWw9V2qwhT20TLaw9hWh5b2y0jTG0TtdtezwvTOE4lUK/ZfKfcsHGdzjJVgvSk41fKpWsv0IfGBap6jBqmdhclrfU2AWqY9nb8mb19AtQwtc+UFnubADVMezv+zN4+AWqY2mdKi71NgBqmfsWfGqbuxgthGovNM889L5/Ytks+tWWDrDru9fp2/TXnntYQptHnVQYqwtTdxc3I/CSAMPUzbozaXQIIU3djw8j8JIAw9TNujNpdAghTd2PDyPwkgDD1K26uCVNVkvKRb3+/qRSlX0SD0SYlO+adB8I0QszUMzWC1Px8+cUXNmqaLhCmj3xK5NFPydyiUZkbGBUZGNXfy+Di4P+BxTI3OCKyaHHwfPh4n3le/a8eV+cPmvMXB+0MmP/Dx9Vx0pc3xhwPAa8ILBkdkMNj0zLn1ajdGezc3Jz09fE+4U5E6h/JUSMDMjYxLbO8qOoPRg+OoBvfkxYPL5LJqVmZ5kXVgyuaKZdBYHRokUzPzMrUTPm/qObmRLhMKiOKtOkSgZHBfn3dNzk969KwGEsKgaWjA5WyUU5r4423y49feKnR7xuOO6ZRlrIuYWr828pjj26UwCwCBmFahF7s3KTgxAWqOiUuTNWmT1uX/IHFkbRpatFIKFJHQkkbl6oR2aqOHVwsfUraavmqzg0lbkPqjgRS1zw+EPkZOVtdXOmpQeAzf3iXXLX+92RkdDFUOiCgPmqgSzsA18Wn/Pmf/FdZ+/5L5Ohj/pcunqX/U/vBD56Vxx/7plx6+RX+TyYygznpk74u+xPYFz5/v/zi29fI8cef1FWx6rbJ/PSnB+RvH/p/5T9f9bvdNrWum8+X//ZBOfGkk+UXTjm99LnN9fVJn7KmfOUmcOTwYfnT+/6rXP3R/z33uZxQLYGv//3XZOmyZXL2286rtmN6y03gjh23ytatW3Of1+kJZkPz+EbnSi5+7sH9Oqv0i3/3CBmmIWAyTEWklcluV8NUCdPrr7tODh86JH0zY9I3c0T6Zsalbyr8X/88Jn1T4+FzY9I3PSain1fHq5+PSN/0ePB/4zH1XOSx6eBnqfBDx1y/EqzzYnYulLVKrpqM2sb3OoN2RGTwKNHH6YxZlWWrvg+ybVUGrf5f/xzIXNMmiqfTt7zuO49Nn7ovpsyoXgJs+lQv/6y9s+lTVlL1H8emT/XHIMsI2PQpCyU3jmHTJzfi0G4UbPrUjpA7z7PpkzuxaDeSKjd9Mpml27ZsaNxBnTQ+k2H6/ne/Q+/vo76iGajmHCNfzc9RCWuyO3/7N98j1//BPU1tfOef/lVu3r5bP3bGqSc3NlhPSlhUx8T7+Z0PrtUZqEmZsp+8cX3L/YfaxSP+fM8L07SgGFDxNF4lUNWXCpD60sJ0001yaGw6L/uOjg8Ea4Jc1UI1/Ddlvg+krUweDuRrQ95GRK2WspFzjbydGatYzg6HcjaUq7okQVCWoCFqjbDV5Q9UJqySs+a48H8jZHU2bdhWQ9YG8lb6+jtiz0nVEECYVsOZXnqHAMLUj1gjTP2IkxolwtSPWCFM/YiTGiXC1I9YIUz9iJMaJcLUn1hVKUyVeNz74P6GoEyjZASlEZPqOOXBfnLg5UZd0/ht+3EZqzzaVR+/TeJt/NFn9y14zPi1JDcXH7M65i8f+nv5wPt+WZ5/4UX52te/Lb/34Yv1VJLGcPvOvW3n22q19LwwTbLSClg0sFGj/b6L1iwofnvoyFRlwrTKl36Q4RqXqyZzNjkzVmfONgTskVhmbbrorTZzdiiQsUly1Yha/Vz4T9eWVfI1yKJtlrgRUduQu+Yx5Gwn65VNnzqhxjkQSCfApk+sDgjYJcCmT3Z50hoE2PSJNQABuwTY9Mkuz7Jbq2rTp7j0bCVM45s+RRMJ1Xlbbt0lm669QsxG6EaqGvmZVD+03WMjw8NNm64ffOWQznC9YeO6lhmx0XlE7xCnhmnZKzdj+90qTDNO38phgZxtLlsQiNfwsUbW7Lyolekwc9aUMmg6JkH06vZV5mx1xa/n+ociGbKRbNh4WYKIvA3KGkSPNSUQQmHbKHEQKXegatD2LbISi7obQZjWHQH67zYCCNNuiyjzqZsAwrTuCNB/txFAmHZbRJlP3QQQpnVHIF//vgnTl392aMGmUWbGJvGwnRxdsXypPiV6XFyYquTGHffcL9tu2iDm+DhZk8kafbzVGPJFRlTNa6pe54UWPx5hWpRgtef3zU5IX1Su6jIH4T9dbzZef3ZMJFpfNkHeNsodmHbCLNtq5ezgfFas3sgrVpZASdVFIzKrNgBbNCqz6mddRzZWU3ZBzdlY/dmS5SzCtNrXA711PwGEaffHmBlWSwBhWi1veut+AgjT7o8xM6yWAMK0Wt5Fe6tKmOa5Jb9VhqkSpp/Ytks+tWVDU4ZplEMVwlRlk+77u0dl5/YbGuOIltAkw7ToyrRwftU1TC0MmSYqJKDlbJjZ2tjcy9SSjclVLV2nx0WmVOZsmF2bIG+TNxcbE5mbqWxmc/2DgXxVNWGVfFVlC7SEVT8rURspaWDkrcmcXZA9azJng5q1f/KXX5Z1l/2GDC1eoQWw9A9UNi86gkA3EqCGqR9RpYapH3FSo6SGqR+xooapH3FSo6SGqR+xooapH3FSo6SGqT+xqrKGaatNn6K1Qb/4d49Iu1vy290qb0OYtrol39Q7vfziC5tu10eYOrb2EaaOBaSHh9M3M5m8uVcoYrWw1aUP1P9BRm0jczZF3jaO1Zm385uDlSlndxy+UT66+G5Z3HckiGbfQHPN2cFoKYKFWbQNWasErcqyDTNmteAdDISuFr0DYZZtVPQqEcwXBLqMAMLUj4AiTP2IE8LUnzghTP2JFcLUj1ghTP2IE8LUnzipkVYpTFV/Zn+e6I72RkyesOpYvV9PO2Gqbo9Pyu5UkvRHzx/Qu9TbEKZqvKqfbz7xdGPjJiN21170Drntv/y5rDz26MaG7PGNpsgwdeC1gDB1IAgMoXICfbNTzZt7qSzaiIid3/grWnNWidogizZJ3gYS94jc+eP/JNcs3yNHzR3Uglfmpqubn5KzOkNWZb2qLNmIkDVZs43NveZryqpNwUyWbWOzsOgGYibL1pQ7UM8hZ6uLa4/3hDD1YwEgTP2IE8LUnzghTP2JFcLUj1ghTP2IE8LUnzjVIUxVn0kbn8c3PW+VYWrqiUY3R1ftvuG4Yxq3x9sSpkaa/tFn9zUCa8ZqRO+TTz2rn1OPm6/rN65LlLZ5Vwc1TPMSSzieGqYWINIEBEIC8RqmSs7KzLj0Tx8RmRmTfp0ZOyb9M+rncemfCv/XPyvpOi59pqxBtDZtKGSTNxcLM24rlbOL5mvODiaUMIjK1Yh8lQGTORuUMJgXvJHM20Z74fOLhlhfPUyAGqY9HHymXgoBapiWgpVGe5gANUx7OPhMvRQC1DAtBWtpjVZVw7S0CXRxwwhTC8FFmFqASBMQSBGmlYKZnQ7rx85nxjaVJQg382raJEyXKzgcyZwNSxdEZK0WvVNjgfANRa8WwVV99Rk5OyJzeqOvFnLVPK/q0SZlxkYzZ5XUjW8uhpytKqqZ+0GYZkbFgRDIRABhmgkTB0EgMwGEaWZUHAiBTAQQppkwOXMQwtSZUCwYCMLUQmwQphYg0gQEXBCmVUZhbiYsaxDWlDU1YlVpA7NRmN4gTGXMhv8r2Tp1JMiy1dm2YXZt+L8pazBf9iCUt1XKWekPMmeNWG1sAJYgVyOZtVrOxjNjtbyNtKW+1z+HbfUPVxkxb/tCmHobOgbuKAGEqaOBYVjeEkCYehs6Bu4oAYSpo4FJGRbC1N14IUwLxoYapgUBcjoEYgT27L5bLrtyvQyPLIaNLQINOasE6hHpmxoP/jclDEz92XAzMFODVkzmbOPYoHSBlri6Hq35OczIVY/PTtoadYZ2jJw1tWZTas42ZcGOSKOsgS5pEEpYLWvD7NtGFm74mN48bCTDeNw8hBqmbsYlPipqmPoRJzXKrzz0eTnrbefLsatO8GfQPThSapj6E3RqmPoRK2qY+hEnNcrHH90vS5Ysk7ecfo4/g+7RkVa96VOPYu542gjTjtEFJyJMCwLkdAggTLtrDczNNstUU8bAyFhTS1bJ1UY27RER9b3JsjWZtXEhqzNvjegdq1jO9rXOdjU1Zwcj5Q4GRkVMlmw0MzYqahM3F7MrZxGmfrzEEKZ+xAlh6k+cEKb+xAph6kesEKZ+xAlh6k+c1EgRpm7HC2FaMD4I04IAOR0CCFPWQAECfVOvhXVnTeZsmEWr5GrWzFlT9mB2QvqmVT3a5szZ/qlXC4yws1OD0gPN8lWXOWiUJQhLFRgRqzJgB4+a3wRMPR6WOPiLv39W3r3mVFm+4hiZ61dZuMMiup0RmR1a3tkAOcs6AYSpdaSlNUiGaWlorTaMMLWKs9TGEKal4rXWOMLUGsrSGyLDtHTE1jpAmFpDWUpDCFMLWKlhagEiTUAgJLDy6FE5cHBMZudAAgGHCOjM2bAmrM6EjWS7xjNjdc3ZQLrqzNl4ZmzT5mFhfVojd9VzsxMVT1xlzyqROhIIVS1hlVQNv180HIraQLbq57V4HdHHiK5X23xusLGYkrLB8bqtsLSBPldl4vaPiPT1VzxXuoNAcQLUMC3OkBYgECVADVPWAwTsEqCGqV2eZbdGDdOyCXfePsK0c3aNMxGmFiDSBAQQpqwBCDQI6ExXJV0bMtbUiQ2zYKP1ZkMBG8jZUOpqARtk2Q73Tcr0+GERLXnDf2ZjsenXaqc+N7AkFLZB1mvjnxKxamMvLVqD73XWrJa6w0G5g1Dc6hqz5tzo8eb7UPRKv8qqXVb7nBmA3wQQpn7Hj9G7RwBh6l5MGJHfBBCmfsUPYepuvBCmFmKDMLUAkSYggDBlDUCgFALHvm5EXj40IdMzKWnbKnt2djzYDEz9r0XreChYlYCdCEWrkrbqe5UFq46Z0Bm0WsQqARs5NxCz8+eKFrzN56pjRGZLmXP7RtWGYZGs1zAbtpFV25RJG82wDTJqlXhtZMyGGbnR7NymrFoteOePF+lrPzyOcJoAwtTp8DA4DwkgTD0MGkN2mgDC1OnwLBgcwtTdeCFMC8aGGqYFAXI6BGIE9uy+Wy67cr0MjyyGDQQgYIGA65s+9U++ooVsQ6pq8apk6+T8xmCqvuyMKnEwnyUrOpvW/BwRsrqtUM6a47XgNTLYgazawaVBdqzOmA2E7DNTx8s3Xz1Z1v3c00FWbVhvVstWU3u2IWuHw0zbsExCo0RC+HijRMJwUDpB96H+2d1QzMLy9LIJapj6ETZqmPoRJzVKapj6EStqmPoRJzVKapj6EytqmLodK4RpwfggTAsC5HQIIExZAxAolYDrwrTUyac1bmrSRrNeTaasFrMTYcZss4hV2ba69EE0CzcqY00mbSh9dTta/ipZG2TqiiRn+j4z82b5xuQa+fDon5eIpC8QtCqDVtWQ1SULIpJVy9mwDEIuYWtq3obtNcojDAfZuPrn7hG2CNMSl6jFphGmFmGW3BTCtGTAlppHmFoCWUEzCNMKIFvqAmE6D/KxJ56Wqz5+m3zyxvVy6doLLBEu1gzCtBg/QZgWBMjpEECYsgYgUCoBhGmpeHM3bmRrUMIglLMz4/Kj//m8fOef/10ueccJYYmDiea6s0rW6qxbVd82fE5l3mopG5ZEUGUSGgI3lL4OZdYaWHMmw7axEViQCRuUMhgKBGtM2MqAenyxiJa9gXw1m42JkrHRrNxQ0gbHmpII6hw7dy4gTHMv+1pOQJjWgr2jThGmHWGr/CSEaeXIO+4QYdoxuspPRJjOI79j5179w08OvCy3bFovoyNDlccj3iHC1EIIqGFqASJNQCAksPLoUTlwcExmU8otAgoCEMhHoG0N03zNcbTnBHTGa5hBG9SfjYhVXfogrGGrha762dSwVdm1Y6HEbZa0pjRCUAohsrmYFrqmHq7qZ6J2enP9Q6IFa7QOrc66DUsXxCWrysTVcjbMyh0clcWLj5Lx2UGZ6gsfyyJslQgeOKr2+TMACLhIgBqmLkaFMflMgBqmfkWPGqYiB185JNs+vUeu/sglsv0PPyubrr1CVp+4qvZAIkwthABhagEiTUAAYcoagEApBBCmpWCl0Q4JNJUn0Nmwk+EmYkGWrJauZhOxhOxZXbu2UTJhXvjq+rYR4Rtk8KpM3Hnh26okQofTyX9a30BTiQJTriDYcCxW0sCUMYiWNBgYbZQ30DVpwzIKQYZtpNxCRPIG2buh3O0fyD9mzoBAyQQQpiUDpvmeI4Aw9SvklQvTnz0jcuhH1UN63WqRpccn9qtux//6o9+V6zeuE5VpetLxK524LR9hamGZIEwtQKQJCCBMWQMQKIUAwrQUrDTqKQG1mZgkiNVGlmyCjA02JQtr0E6PyeiiaZmeOCKzk2NB1myYSRuVwVr+xoRvIGxn6yXXtyh5Q7FoGYRG9m1YgzZSOiGoextk3Dbq00azcheF9WyNwI2XTugfrHf+9O4kAYSpk2FhUB4TQJj6FbzKhenDN4o8tqN6SBdsFzlvU2K/SpK+6+1nynlnnyJKnt6+c6/ce9t1smL50urHGekRYVoQPzVMCwLkdAjECOzZfbdcduV6GR6xU2sOwBDodQLUMPVjBfz7D/9NnnziW/LeX7/cjwH38CiL1DDtm50K5Gt0s7B46YJITVotWSOZuDIVSNoFpROiWbk6UzdSbiHyvczN1By5/oZwNfVlm2rN6qzacIOwaPZsmGUrOsNWZdaGdWwX1LoNH180Ii++Oi5f+/8el8su/U/BObr2bf310GoOgJPdU8PUybAsGBQ1TP2IkxolNUz9iVUtNUy/8xmRf/6L6iGddbXIL/zmgn6fee552XHP/bLtpg1akKrb86/ZfKfcsHGdFqh1fiFMC9JHmBYEyOkQQJiyBiBQKgGEaal4rTWOMLWGsvSGigjT0gfXroPZ6UDWxkofaMGqRey8jNXZs9GatKrEgSqJoB+LPhetXRtm4obHzNe0DTNxlTCu6OuF2ZXyhfHfkKsXfybSY19Y+mBIZ9oGGbMLa9rqx2MyVmfVDowGmbUJtW51HVzd1ry0bcrENSK4ovn71A3C1I9oIUz9iBPC1J84qZHWIkwdQ/TAvofl5u27F4zqdz64Vt+iX+cXwrQgfYRpQYCcDgGEKWsAAqUSQJiWitda4whTayhLb8hrYVo6nfYd9M0cCerQhpm0C0oXhNm3SZm4QYZtkD0bSN4wkzYUvo2atdPjcmBisTx48O3ye8vvn6+POzvZfoAlHxEI1YTSBWEd20b2bLQOrcmwVZmyOss2EL6mPm0gccNyCeEmZkb4BnVuQ8m7aLTk2XXWPMK0M25Vn4UwrZp45/2RYdo5u6rP7HVhOjY+KVt37JY1557WVLM0nnVadVxMfwhTC+SpYWoBIk1AICSw8uhROXBwTGbnQAIBCNggQA1TGxRpAwLzBI5ZNiyvjU3JxFTN9UgJSkcE+mbmNw5rKl0Q2zhMl0JIyMQNsm7nhW9D/kZq3TYya+PlFlTWbs1fc/1DYSZsgmRtKneQvImYDI7OZ9GGmbha8kazcpWg1XI3zN418ndRcrklapjWvCjovusIUMPUr5BWXsPUITxKjH5i2y751JYNsvrEVY2RGZF6+cUX1npbPsLUwmJBmFqASBMQQJiyBiBQCgGEaSlYabSHCSBMezj4FqYe3RwsqE/bXNKgsYlYU/ZsmFWrNxObF75BLdugJq7EhK8pvRA8bo5Twrbev0jPqY2/GmULgvIF/YMjMtuv/s3XsNWbi0Xr2ZrvdYbtfNmDIIM2FLbRc/Rj0XILwffSP2AhijQBAbcJIEzdjk98dL0sTF2PFMLUQoQQphYg0gQEEKasAQiUQgBhWgpWGu1hAgjTHg5+F0y9b2ZStJTVsjWUqSZjVmXJhgI2KmP15mPhMUEd21a1bifTyy2odqTmzOy+Rcn1aaOlC1qUOwgybAOZ2yifEM2kTSq3EM3EVcKYLwiUTABhWjJgy80jTC0DtdgcwrQgTGqYFgTI6RCIEdiz+2657Mr1MjySfNsWwCAAgXwEqGGaj1ddR1PDtC7y+fulhml+ZnWc8dKLB2T/V/fJB664qo7u6TOFQJ/a+EsJ17BWrZKxX//GY3L8ccvl5DcsTyx3EM3E1cJ2NpSy0+OSVDpBZ9xGsm+bMnHnZmqOTX+w8VhEsjaVLjDZs2oTMZUlG8ueDWrYhhuMNdWnDWrV6rbCerjxDc2CzN6hjudPDdOO0VV+IjVMK0fecYe9XsO0Y3AVnYgwLQgaYVoQIKdDAGHKGoBAqQQQpqXitdY4wtQaytIbQpiWjthKBwhTKxgraaTqTZ/6pl5rErZ9CbVmF5Q7CCWvzExJ3/ThRoZuLmE7c7gSnu06aWwCpgWrEq1DEckakbExsXt4dkT+9LtL5OpfGgjKIjRk7rA0Nh6LlkJo1LE1pRFUmYSj2g2P5y0QQJhagFhREwjTikB32A3CtENw5jSEaUGAnA4BhClrAAKlEkCYlorXWuMIU2soS28IYVo6YisdIEytYKykkaqFaSWTatGJFq669IGpPxvbRCzMvo1n4qrjZUpl2IY1bVUWbUqt20b7kdILQXuTHU//8NwSuffI1fL7R/3fHbehTpyvY9ucDTtfj7a59qzOjDWbiqmSB2Ed20ZJBF1CwbQVlkowZRX0/0OB0I3WuC00A/dPRpi6HyMzQoSp27FCmFqIDzVMLUCkCQiEBFYePSoHDo7JbL17EhAPCHQNAWqYdk0omYgjBKhh6kggGEbXEFixZEjGp2ZkbKLu2+W7BmnLifTNHJnfbExvGJZD2OpjxxJq3QZtmFq3OkNX18qdCDcki2xQ5gjmuYElYXmEoJxBU13ayMZiuoxBo66tErBH6c3DTGmFeRkbZurGSyOYzcx0mYURmR1aVjoBapiWjthqB9QwtYrTamO1CtODrxySazbfKU8+9eyCSZ1x6sly723XyYrlS61OuIzGEKZlUKXNXiWAMO3VyDPvsgggTMsiS7u9SgBh2quRZ95lEUCYlkXW3Xb7pg5FNh4Ls21VTVolWs2mYrpUgsnCDR9XWbZ647HgZ328OS7MuNVlE2bGwk3NYhucTR1yCsrs4DLRmbU6gzZa5mBUZNFQ8FikNEJwbHCcybQ1NWuDTNzg+MVLlsvMXL+MzQwEJRfU8YsGwzIKwyJ9Su4Ggpev+gkgTOuPQdoIahWmd+zcq8d1/cZ17hLKMDKEaQZIHAKBjAQQphlBcRgEMhJAmGYExWEQyEgAYZoRFIdBICMBhGlGUBxmjYCRrI0NxVRpAyNqdaZsuLFYpORBkD0bClxdGiGeVRvKWS1zJ4MNyaaN6I1k2BYsjWANQtjQfImEQKg2RG3/UETaqixbJViVrA1kq9lELCpjTU1c3YaSsgOqvq2SuGHbpo3web1JmT4ubM98Hwpd23N1tT2EqauREalNmKrs0i237pJN114hq09c5QQhJXBPOn6lXLr2gsZ4nnnuedl44+3y4xdeajwWzX6lhqkToWMQXURgz+675bIr18vwyOIumhVTgUB9BKhhWh/7PD1TwzQPrXqPpYZpvfyz9k4N06yk6j+u12qY1k+8sxGMHTksf3n/ffLh9R/trAHOaiIQlEaY1HVodaZsmE2rJevcpJatQcZtJIvW1K/VpRRUpm0odtU5Kvs2zLzdf+B4Wdb/mpy35H8EZRJmVT+TgQjWole1/5oXEQlKJ5gsWfV/8E/XsI3K2EZGbrPcNQJYy9lopq2qhxvK30ACx+VuRPRG2ra9cRk1TN1ehghTEXlg38Ny8/bdOlKfvHH9AmH6iW275FNbNiSKXYSp2wuc0flHAGHqX8wYsdsEEKZux8eMDmHqR5zUKBGmfsQKYepHnNQoEaZ+xAph6kec1CjzbvrUN/VaQ6SKFqpT8zVolWwNReu8cJ2S4LhQxuqM2smwXq0St+b8QNCq5xptzKqSCWF9W/P4nDo3LM2gSyoc9gL23ODSUNtw2QUAACAASURBVNrGyh7EZGwruXvHN5bK1q1bvZhvLw6yNmGqYCdldNYZhLQMU4RpnVGh714jgDDttYgz37IJIEzLJmynfYSpHY5VtIIwrYJy8T4QpsUZVtUCwrQq0sX6QZgW41fl2XmFaZVjy9NXkDUbZMcGcrZZxgbS1mTPzsvb+fMC+auOE52FG4hbfV7j+/Zt63q5Uek7Z2+Dulte+4OeF6aPPfG0XPXx25qWxu98cK0TpTtrFabqdvc9D3xVNl1zhYyODOV57ZRybJZb8pM2o6KGaSnhoNEeJUAN0x4NPNMujQA1TEtDS8M9SoAapj0aeKZdGgFqmJaGloZ7lMDyowZlemZODo9P9yiB8qetSyqEMlaXU1BlFMJMW5N5Gy2B0Cxp1YZmk/qciePeKa8/9aLyB+xwD0qY3r5zb2PTd7M5/A0b18l5Z59S68hrE6YGwpNPPZsIIElMlk0qS8arOuYnB16WWzatd0Lyls2E9iEAAb8ITE3PyuBAv1+DZrQQgEDXEpicnpUh3pO6Nr5MDAK+EZianpPBgT7fhs14IQABCHQtgbgwHRuflK07dsuac09rKpdZB4DahGkdk23XZxZhqrJid9xzv2y7aYOsWL5UN0mGaTuyPA+B7ATIMM3OiiMhkIUAGaZZKHEMBLITIMM0OyuOhEAWAmSYZqHEMRDIToAM0+ysXDhy1TGjlQ7jmYPPyI9e+VGlfarOVh+9Wo5fdvyCfpMyTF3ZIB5hGglXJ8KUTZ8qf53RYZcToIZplweY6VVOgBqmlSPvqENqmHaErZaTqGFaC/bcnVLDNDey2k6ghmlt6HN1TA3TXLhqPbhbapjWCrGiznfdvb3yGqY3fuVG2fGPOyqa4Xw329+9XTb90qZEYRqvYfqG446RndtvSNx4vcqB1y5Mkwq83nfX5lpqFSQJ0y/t/6a8+U0/1wiUOkZ9Xb9xnf4fYVrlcqWvXiCAMO2FKDPHKgkgTKuk3XlfCNPO2VV9JsK0auKd9Ycw7YxbHWchTOugnr9PhGl+ZnWdgTCti3z+fusQpp/51mfkL/7pL/IPtuAZV//i1fKbp/9mojCN1jBVB8SzTgt23fHptQrTJAjqlveNN94u137kksrqFTyw72G5efvuBsSozY4L3fddtKapfinCtOO1x4kQSCSAMGVhQMAuAYSpXZ5ltYYwLYus/XYRpvaZltEiwrQMquW0iTAth6vtVhGmtomW1x7CtDy2tluuQ5jankPR9pK8oNrzyIXb8msTpqaQ6+UXX7ggm1QB+9yD+73ZWIkapkVfIpwPgXkC1DBlNUDALgFqmNrlSWsQoIYpawACdglQw9QuT1qDADVM/VoDVdcwdY1OkjDt+QzTVsY4aWMl14IaHQ/C1OXoMDbfCCBMfYsY43WdAMLU9QgxPt8IIEx9ixjjdZ0AwtT1CDE+3wggTP2KGML0aYnXMD3j1JPl3tuua2y0XldEyTC1QB5hagEiTUAgJIAwZSlAwC4BhKldnrQGAYQpawACdgkgTO3ypDUIIEz9WgO9LkxdjlZtwlRBUbVD9z64v8kc11HDtEiAqGFahB7nQmAhAWqYsiogYJcANUzt8iyrNWqYlkXWfrvUMLXPtIwWqWFaBtVy2qSGaTlcbbdKDVPbRMtrjxqm5bG13TI1TG0TtdtercJUTSW+qZJ67L67Ni+oa2p32vZaQ5jaY0lLEFAEEKasAwjYJYAwtcuzrNYQpmWRtd8uwtQ+0zJaRJiWQbWcNhGm5XC13SrC1DbR8tpDmJbH1nbLCFPbRO22V7swtTud6ltDmFbPnB67mwDCtLvjy+yqJ4AwrZ55Jz0iTDuhVs85CNN6uOftFWGal1h9xyNM62Ofp2eEaR5a9R6LMK2Xf57eEaZ5aFV/LMLUAnNqmFqASBMQCAlQw5SlAAG7BKhhapcnrUGAGqasAQjYJUANU7s8aQ0C1DD1aw1Qw9TdeCFMLcQGYWoBIk1AAGHKGoBAKQQQpqVgpdEeJoAw7eHgM/VSCCBMS8FKoz1MAGHqV/ARpu7Gq3JhevCVQ3LN5jvlt3/zPfLHf/G38uRTzybSOePUk5s2g3IXoQjC1OXoMDbfCJBh6lvEGK/rBBCmrkeI8flGAGHqW8QYr+sEEKauR4jx+UYAYepXxBCm7sarcmFqUChxuuXWXbLp2itk9YmrmgipjaA+9+B+uWXTehkdGXKXnohQw9Tp8DA4DwlQw9TDoDFkpwlQw9Tp8DQGRw1TP+KkRkkNUz9iRQ1TP+KkRkkNUz9iRQ1TP+KkRkkNU39iRQ1Tt2PlpDB95rnnZcc998u2mzbIiuVLnSaIMHU6PAzOQwIIUw+DxpCdJoAwdTo8CFM/wtM0SoSpH0FDmPoRJ4SpP3FCmPoTK4SpP7FCmLodKyeF6QP7HpZHvv19MkzdXjuMDgKlEECYloKVRnuYAMLUj+CTYepHnNQoEaZ+xAph6kecEKb+xAlh6k+sEKb+xAph6nasKhemKnt04423y49feCmVzBuOO0Z2br9hwa36rqKkhqmrkWFcPhKghqmPUWPMLhOghqnL0WFsPhKghqmPUWPMLhOghqnL0WFsPhKghqlfUev1GqaqJOdVH79tQdA+eeN6uXTtBbUGs3JhambbqoZprUQ66Bxh2gE0ToFACgGEKUsDAnYJIEzt8qQ1CCBMWQMQsEsAYWqXJ61BAGHq1xpAmD4tt+/c6+Sm77UJU7+WcOvRIky7KZrMpW4CCNO6I0D/3UYAYdptEWU+dRNAmNYdAfrvNgII026LKPOpmwDCtO4I5OsfYYowzbdiPDqaTZ88ChZD9YIANUy9CBOD9IgANUz9CBY1TP2IkxolNUz9iBU1TP2IkxrlP+7/orzx+JPkxNWn+jPoHhwpNUz9CTo1TP2JVR01TGcPzsrsq7OVQ+pf0S/9y/oX9KtuySfDNCEcreqZnnHqyU6m5MangTCt/HVGh11OAGHa5QFmepUTQJhWjryjDhGmHWGr5SSEaS3Yc3eKMM2NrLYTEKa1oc/VMcI0F65aD0aY1oo/V+d1CNOxr47J+DfGc43TxsGjF43KyC+NJApTapjGsIyNT8rWHbtlzbmnyVmnv1n2PPBV2XTNFTI6MiR37Nwr73r7mXLe2afYiEupbSBMS8VL4z1IAGHag0FnyqUSQJiWitda4whTayhLbwhhWjpiKx0gTK1grKQRhGklmAt3gjAtjLCyBhCmlaEu3FEdwnTi2xMy+U+Thceet4Hhc4dl6PShRGFKhmkMS3TTJ/XUjnvul203bZAVy5eKSsn93IP75ZZN67VAdf2LGqauR4jx+USAGqY+RYux+kCAGqY+RIkx+kSAGqY+RYux+kCAGqY+RIkx+kSAGqY+RUuEGqbckr9gxUaF6dGvWyrbPr1HtnzsQ1qYqlv1owLV9eWOMHU9QozPJwIIU5+ixVh9IIAw9SFKjNEnAghTn6LFWH0ggDD1IUqM0ScCCFOfooUwpYZpwnqN3pJ/6doL9G34Jx2/UtT3D+x7WB759vfJMPXrdc5oIWCFAMLUCkYagUCDAMKUxQABuwQQpnZ50hoEEKasAQjYJYAwtcuz7NbIMH1aqGHaZpWpjNNrNt8pTz71rLzhuGNk5/YbZPWJq8pem4Xbp4ZpYYQ0AIEmAtQwZUFAwC4Bapja5VlWa9QwLYus/XapYWqfaRktUsO0DKrltEkN03K42m6VGqa2iZbXHjVMy2Nru+U6apjankM3t9c3Nzc3180TLHtuCNOyCdN+rxFAmPZaxJlv2QQQpmUTttM+wtQOxypaQZhWQbl4HwjT4gyragFhWhXpYv0gTIvxq/JshGmVtIv1hTAtxq/ss2sTptEapj5kkqYFAmFa9hKl/V4jgDDttYgz37IJIEzLJmynfYSpHY5VtIIwrYJy8T4QpsUZVtUCwrQq0sX6QZgW41fl2QjTKmkX6wthWoxf2WcjTC0QZtMnCxBpAgIhAWqYshQgYJcANUzt8qQ1CFDDlDUAAbsEqGFqlyetQYAapn6tgV6vYepytGoTpgqK2ujpXW8/U847+xSXGbUdG8K0LSIOgEBmAgjTzKg4EAKZCCBMM2HiIAhkJoAwzYyKAyGQiQDCNBMmDoJAZgII08yonDgQYepEGBIHUaswfea552XPA1+VTddcIaMjQ+5SajMyhKm3oWPgDhJAmDoYFIbkNQGEqdfhY/AOEkCYOhgUhuQ1AYSp1+Fj8A4SQJg6GJQWQ0KYuhuv2oSpqmF6zeY75cmnnk2kc8apJ8u9t10nK5YvdZeeiFDD1OnwMDgPCVDD1MOgMWSnCVDD1OnwNAZHDVM/4qRGSQ1TP2JFDVM/4qRGSQ1TP2JFDVM/4qRGSQ1Tf2JFDVO3Y1WbMHUbS/bRIUyzs+JICGQhgDDNQoljIJCdAMI0O6s6j0SY1kk/X98I03y86joaYVoX+fz9IkzzM6vjDIRpHdQ76xNh2hm3Os5CmNZBPXuftQlTlWG65dZdsunaK2T1iauaRvzYE0/L5x7cL7dsWu/8rfoI0+yLjSMhkIUAwjQLJY6BQHYCCNPsrOo8EmFaJ/18fSNM8/Gq62iEaV3k8/eLMM3PrI4zEKZ1UO+sT4RpZ9zqOAthWgf17H06KUxVbdMd99wv227a4Pwt+Qo1NUyzLziOhEA7AtQwbUeI5yGQjwA1TPPx4mgItCNADdN2hHgeAvkIUMM0Hy+OhkA7AtQwbUfIreepYepWPKKjcVKYPrDvYXnk29+vPMP0jp175aTjV8qlay9oipgaz83bd+vH3nfRmgXjQpi6u8AZmX8EEKb+xYwRu00AYep2fBidfwQQpv7FjBG7TQBh6nZ8GJ1/BBCmfsUMYepuvCoXpip7dOONt8uPX3gplcobjjtGdm6/YcGt+mVhjArRT964vkmYqvIAt+/c29iASklV9XX9xnWN4SBMy4oM7fYiAYRpL0adOZdJAGFaJl3a7kUCCNNejDpzLpMAwrRMurTdiwQQpn5FHWHqbrwqF6YGRasapnXhSsowjT8WF6jUMK0rWvTbrQSoYdqtkWVedRGghmld5PP1Sw3TfLzqPJoapnXSz943NUyzs6r7SGqY1h2BbP1TwzQbJxeOooapC1HINgZqmGbjVNdRtQnTuibcqt+4HB0bn5StO3bLmnNPa2SdqgzZT2zbJZ/askFnwCJMXYwkY/KZAMLU5+gxdhcJIExdjMrCMSFM/YiTGiXC1I9YIUz9iJMaJcLUj1ghTP2IkxolwtSfWCFM3Y4VwjQSnzRhevnFF8p5Z5+ij0wSpptvulkmp2fdjjSjg4AnBD5z951y1fqNMrJ4sScjdmuYc3Nz0tfX59agGE2tBP78vl2y9v2XyNGvP7bWcdB5awI/+Ldn5fFvPSqXXv7BrkLVje9JX/j8Z+W8898hP3fCSV0Vq26bzE8PvCB/u+9B+c9X/W63Ta3r5vOVLz4oJ570JnnLqW8tfW5zcyJcJnWG+cjhw/Kn9+2Sqz/68c4a4KzKCPzD/q/K0mXL5axzzqusTzrqjMAd2z8lW7du7exkziqdQK3CVN2Wf83mO+XJp55dMNEzTj25UTe0dAphB51kmKpTJ6ZmZHIKYVpVnOinuwksGR2Uw+NToi5o+cpPQGFDl+bn1s1nHDUyIGOTMzI7y4uqm+Ps6tzmpE/6pLvW3ujwgExOz8jMTHfNy9U1xLi6n8DI0CKZnp2T6QoSUOb6+qSPi8zuX1Q9PsPhoUX6um+qgtdUj6O2Mv2liwettEMj9gnUKkyTNlCyP8XsLXZSw1S1zqZP2RlzJATaEWDTp3aEeB4C+Qiw6VM+XhwNgXYE2PSpHSGeh0A+Amz6lI8XR0OgHQE2fWpHyK3n2fTJrXhER1ObMPVl06f4Jk9Jkhdh6u4CZ2T+EUCY+hczRuw2AYSp2/FhdP4RQJj6FzNG7DYBhKnb8WF0/hFAmPoVM4Spu/FCmIrIA/selpu3725E6Q3HHSM7t9+gN3VSX9Hn33fRGrll03oZHRnSz7Hpk7uLm5H5SYBNn/yMG6N2lwCbPrkbm+jI2PTJjzipUbLpkx+xYtMnP+KkRsmmT37Eik2f/IiTGiWbPvkTKzZ9cjtWtQlThSXpFni3cS0cHcLUt4gxXtcJIExdjxDj840AwtSPiCFM/YgTwtSfOCFM/YkVwtSPWCFM/YgTwtSfOKmRIkzdjletwlTtOL/nga/KpmuuaGRsuo0LYepbfBivfwQQpv7FjBG7TQBh6nZ8zOgQpn7ECWHqT5wQpv7ECmHqR6wQpn7ECWHqT5wQpu7HqjZhqmqYXrP5TnnyqWcTKZ1x6sly723XyYrlS52nSA1T50PEAD0iQA1Tj4LFUL0gQA1TL8LEID0iQA1Tj4LFUL0gQA1TL8LEID0iQA1Tj4IlItQwdTdetQlTd5HkHxnCND8zzoBAGgGEKWsDAnYJIEzt8qQ1CCBMWQMQsEsAYWqXJ61BAGHq1xpAmLobL4SphdggTC1ApAkIhAQQpiwFCNglgDC1y5PWIIAwZQ1AwC4BhKldnrQGAYSpX2sAYepuvGoVpmPjk7J1x2556GuPiNmZftVxr9ePrTn3NLl07QXukgtHxqZPzoeIAXpGgBqmngWM4TpPgBqmzodID5Aapn7ESY3yKw99Xs562/ly7KoT/Bl0D46UGqb+BJ0apn7EihqmfsRJjfLxR/fLkiXL5C2nn+PPoHt0pGz65HbgaxWmd+zcKycdv1Le+6trZMe998uHLv01WX3iKnnsiaflcw/ul1s2rXd+MyiEqdsLnNH5RwBh6l/MGLHbBBCmbsfHjA5h6kecEKb+xAlh6k+sEKZ+xAph6kecEKb+xEmNFGHqdrxqE6Zq06ctt+6STddeISqrNCpMn3nuedlxz/2y7aYNzm/6hDB1e4EzOv8IIEz9ixkjdpsAwtTt+CBM/YhPdJRkmPoRM4SpH3FSo0SY+hErhKkfcUKY+hMnhKn7sXJSmPqUYapCTA1T9xc6I/SHADVM/YkVI/WDADVM/YgTo/SHADVM/YkVI/WDADVM/YgTo/SHADVM/YmVGik1TN2NV23CVCF5YN/D8si3vy9bPvYh+S+7v6BvyT/6dUvlms13yrqLL/SihinC1N3Fzcj8JIAw9TNujNpdAghTd2PDyPwkgDD1M26M2l0CCFN3Y8PI/CSAMPUrbghTd+NVqzBVWFQ26VUfv62J0H13bZbzzj7FXWqxkZFh6k2oGKgHBBCmHgSJIXpFAGHqVbgYrAcEEKYeBIkhekUAYepVuBisBwQQph4EKTJEhKm78apdmLqLJtvIqGGajRNHQSArAWqYZiXFcRDIRoAaptk41X0Umz7VHYHs/VPDNDurOo+khmmd9PP1TQ3TfLzqOpoapnWRz9/v44/ulyVLlslbTj8n/8mcUSkBNn2qFHfuzmoVpnfs3Cs/OfCy3LJpvYyODOnBj41PytYdu2XNuad5cUs+wjT3muMECLQkgDBlgUDALgGEqV2eZbWGMC2LrP12Eab2mZbRIsK0DKrltIkwLYer7VYRpraJltcewrQ8trZbRpjaJmq3vdqEqRGjl1984YLb733a9AlhandB0hoEEKasAQjYJYAwtcuzrNYQpmWRtd8uwtQ+0zJaRJiWQbWcNhGm5XC13SrC1DbR8tpDmJbH1nbLCFPbRO22V5swPfjKIdly6y7ZdO0VsvrEVU2zeua552XHPffLtps2yIrlS+3OuITWqGFaAlSa7FkC1DDt2dAz8ZIIUMO0JLA027MEqGHas6Fn4iURoIZpSWBptmcJUMPUr9BTw9TdeNUmTLslw1SFFmHq7gJnZP4RQJj6FzNG7DYBhKnb8WF0/hFAmPoXM0bsNgGEqdvxYXT+EUCY+hUzhKm78apNmCok6tb7Ldt2yc7tNzSyTFV26cYbb5drP3KJFzVMEabuLm5G5icBhKmfcWPU7hJAmLobG0bmJwGEqZ9xY9TuEkCYuhsbRuYnAYSpX3FDmLobr1qFqcJiBOmPX3ipQem+uzYvqGvqKkJqmLoaGcblKwFqmPoaOcbtKgFqmLoameZxUcPUjzipUVLD1I9YUcPUjzipUVLD1I9YUcPUjzipUVLD1J9YUcPU7VjVLkzdxtN+dAjT9ow4AgJ5CCBM89DiWAi0J4Awbc/IhSMQpi5EIdsYEKbZONV9FMK07ghk7x9hmp1VnUciTOukn69vhGk+XnUejTCtk377vhGm7Rm1PAJhWhAgp0MgRgBhypKAgF0CCFO7PMtqDWFaFln77SJM7TMto0WEaRlUy2kTYVoOV9utIkxtEy2vPYRpeWxtt4wwtU3Ubnu1CtODrxySazbfKU8+9eyCWZ1x6sly723XyYrlS+3OuITW2PSpBKg02bMEqGHas6Fn4iURoIZpSWBptmcJUMO0Z0PPxEsiQA3TksDSbM8SoIapX6Gnhqm78apVmN6xc68mc/3Gde4SyjAyhGkGSBwCgYwEEKYZQXEYBDISQJhmBMVhEMhIAGGaERSHQSAjAYRpRlAcBoGMBBCmGUE5chjC1JFAJAyjNmGqsku33LpLNl17haw+cZW7hDKMDGGaARKHQCAjAYRpRlAcBoGMBBCmGUFxGAQyEkCYZgTFYRDISABhmhEUh0EgIwGEaUZQjhyGMHUkEAhT+4Gghql9prTY2wSoYdrb8Wf29glQw9Q+0zJapIZpGVTLaZMapuVwtd0qNUxtEy2vPWqYlsfWZsvUMLVJs9y2qGFaLl+brVPD1CZN+23VlmGqpqJuyT/p+JVy6doL7M+sohYRphWBppueIYAw7ZlQM9GKCCBMKwJdsBuEaUGAFZ6OMK0QdoGuEKYF4FV8KsK0YuAddocw7RBcDachTGuA3mGXCNMOwVV0Wq3C9Jnnnpc9D3xVNl1zhYyODFU0ZbvdIEzt8qQ1CCBMWQMQsEsAYWqXZ1mtIUzLImu/XYSpfaZltIgwLYNqOW0iTMvhartVhKltouW1hzAtj63tlhGmtonaba82YapqmF6z+U558qlnE2d0xqkny723XScrli+1O+MSWqOGaQlQabJnCVDDtGdDz8RLIkAN05LA0mzPEqCGac+GnomXRIAapiWBpdmeJUANU79CTw1Td+NVmzB1F0n+kSFM8zPjDAikEUCYsjYgYJcAwtQuT1qDAMKUNQABuwQQpnZ50hoEEKZ+rQGEqbvxQphaiA3C1AJEmoBASABhylKAgF0CCFO7PGkNAghT1gAE7BJAmNrlSWsQQJj6tQYQpu7Gq3Zh+tgTT8tVH7+tidB9d22W884+xV1qkZFRw9SLMDFIjwhQw9SjYDFULwhQw9SLMAk1TP2IkxolNUz9iBU1TP2IkxolNUz9iBU1TP2IkxolNUz9iRU1TN2OVa3CVMnS23fubapVqjaC2njj7XLtRy6RS9de4DY9EUGYOh8iBugZAYSpZwFjuM4TQJg6HyI9QISpH3FCmPoTJ4SpP7FCmPoRK4SpH3FCmPoTJzVShKnb8apNmI6NT8rWHbvl8osvXJBNqkTq5x7cL7dsWi+jI0NOE0SYOh0eBuchAYSph0FjyE4TQJg6HZ7G4BCmfsQJYepPnBCm/sQKYepHrBCmfsQJYepPnBCm7seqNmF68JVDsuXWXbLp2itk9YmrmkipLNMd99wv227aICuWL62Vosl4/fELLzXGccapJzdlxVLDtNYQ0XmXEaCGaZcFlOnUToAaprWHgAF0GQFqmHZZQJlO7QSoYVp7CBhAlxGghqlfAaWGqbvxqk2Y+pJhqoTpJ7btkk9t2bBA7JqwIkzdXeCMzD8CCFP/YsaI3SaAMHU7PozOPwIIU/9ixojdJoAwdTs+jM4/AghTv2KGMHU3XrUJU4XkgX0Py94H9ztdwxRh6u7iZWTdSQBh2p1xZVb1EUCY1seenruTAMK0O+PKrOojgDCtjz09dycBhKlfcUWYuhuvWoWpwqLqlV718duaCN131+YFdU3rQhi/JT9+Oz41TOuKDP12KwFqmHZrZJlXXQSoYVoX+Xz9UsM0H686j/7KQ5+Xs952vhy76oQ6h0HfbQhQw9SfJUINUz9iRQ1TP+KkRvn4o/tlyZJl8pbTz/Fn0D06UjZ9cjvwtQtTt/EsHN0dO/fKTw683NiQSgnTrVu3+jYNxgsBZwns2LFDPvrRj8rixYudHaPLA5uanpXBgX6Xh8jYKiZw7733ygc+8AE59thjK+6Z7vIQeOaZZ+Qb3/iGfPjDH85zmvPHTk7PylCXvSf92Z/9mbzzne+Uk08+2Xn+vTzAF154Qb7whS/I1Vdf3csYvJi7itPq1avlzDPPLH28U9NzMjjQV3o/3djB4cOHRV1T/P7v/343Tq+r5vSlL31Jli9fLmvWrOmqeXXjZPBJbke1VmEal48Klaltuubc0+TStRc4Ry++IRUZps6FiAF5ToAMU88DyPCdI0CGqXMhSRwQGaZ+xEmNkgxTP2JFhqkfcVKjJMPUj1iRYepHnNQoyTAtP1Zj00dkamZKJmYnZHJ6QqZmJ2VyZlIm1f/TEzKpHp+ZksmZCZmamQyOm5lsfK8eUz+PfGOOBLzyw9VxD7UJU182fYqTjQtT9TybPnW8/jgRAgsIUMOURQEBuwSoYWqXJ61BgBqmrAEI2CVADVO7PGkNAt1Sw1RJx4npyVBGKuEYSsjZUELqn4PnJ6LSMnxcna+l5oxqZzxoZzqUmlp0ht/r49PbDgRoKERnJmRmbsbaIvv98/8P2fHeT1prj4bsEqhNmB585ZBsuXWXbLr2igW7zydJSbvTzt7al/Z/U978pp9rjFFlxaqv6zeuazSCMM3OkyMh0I4AwrQdIZ6HQD4CCNN8vDgaAu0IIEzbEeJ5COQjgDDNx4ujIdCOQF5h+trUa82ZkFFJ2ZCMgZwMMiaTMysnpsYXZFZq6RlmVkazLrWgTGn7yPThdlN04vmlQ0tlqH9YBhcNydCiYRlaNKh/AkZwbQAAIABJREFUHhoYlqF+9diQDPYPybD+WR03GB43FB43JBee8G657Mz/6MR8GMRCArUJU18yTOObUr3vojWN+qUGJ8KUlxYE7BFAmNpjSUsQUAQQpqwDCNglgDC1y5PWIIAwZQ30IgF1S7cSkBMz4zI+M64zJNX3KptRZ0Sqx6fH9f/qZyUY9ffh8/r86eB8lT2pjlHtqONmZFLGpyaC803mZSRDUvV7eOo1L7AfNbhEi8dASCrR2F5CBrJyOJSVQzKyaFQGGrJyodTUxy4akuGI1Ix+H/SrROiQLB44yjq3VceMWm+TBu0QqE2YquErGbll2y7Zuf2GRgan2ZX+2o9c4mQN0zh2apjaWYi0AgFDgBqmrAUI2CVADVO7PMtqjRqmZZG13y41TO0zLaNFapiWQbWcNqlhWg5X2612Uw3T4FbvUDCGotLcth0IylBERkSmEZVaZE6NNaSmkZtGeAbSczKUoIHI1I8ZCRrWtLQdn2h775H3yM/kZ/KIPNK2m4H+QRnWwjAQiSobUmVFmgzJqKhUmZP6OSUXQ8moRWIoI4cHRhqSMjhuuKntqPg0mZdN7YVtaym6aKjt2LvhgF13b6eGqcOBrFWYKi5GkP74hZcamO67a7Ocd/YpDmObHxrC1IswMUiPCCBMPQoWQ/WCAMLUizAJwtSPOKlRIkz9iBXC1I84qVEiTP2IlU1hemjyUCMzMrhlO6gx2ZRpqR+LZloGG+lMTU+KumVbZ1ROhyKyIT3Vz1MyPjPWyLqcz8AcF9WvS19Lh5bJyKIRUaLRyEX1/fCiERkJb+NW3yuJqR9bNNLIchwcCKRl8PyIzo40x439y0sydNSIHLf6pJa3i6tz+aqXAMK0Xv7teq9dmLYboOvPI0xdjxDj840AwtS3iDFe1wkgTF2PUDA+hKkfcUKY+hMnhKk/sUKYVhsrdSt3NBtyPCIbG49reWluAQ9uC58cn5DBJ2blxTNflfg5Onsydo4RoY3bysN+qp1tem9LBpfMy0YlHZV8DMVk438lKBepW7pDKbloRI4aPEoWLRoIzx0OnguP0/IylJ1GZEZFqDpWSdKyvx5/dL8sWbJM3nL6OWV3RfsFCSBMCwIs+XSEqQXA1DC1AJEmIBASoIYpSwECdglQw9QuT1qDADVMWQMQsEug12qYdiosdUalvhU8qEupszGjt4xHbyUPa1pGsyuVAJ2ambQbvA5aU7eAj+h6lEpQBoJRCUqdPRnJkgyeC27pbhy3aFhGBkeDc/Rt4fPna2mp24lkXUaFZyRTs4Nhe3VK3k2fvJpcFw6WGqbuBhVhaiE2CFMLEGkCAghT1gAESiGAMC0FK432MAGEaQ8Hn6mXQqBqYap2BI+KxPnalfO3gEdFZJAhGWRbTs9MyeEpdTt4sNFOIC2D7Mqmc+Iic3pC30buwlf0Fu+GkIxkUGrhGBGZUSk5Ori4cev3vPAMsyzDc4L6ls2ZlsGxw6VsmOMCU9fGgDB1LSKtx4MwdTdeCFMLsUGYWoBIExBAmLIGIFAKAYRpKVhptIcJIEx7OPhMvTCB6VlVX7L5Vu/hoVl5deKwvDLWfKt4o35lKB/Vbd9Nm+3kFJZKaM7MzRSeQ5EG+vv6GxmVC27fjmRANmpSxrMwB4LsSlXfMp5dGb8VvKm+pa5/GewEzlf3E0CY+hVjhKm78UKYFowNNUwLAuR0CMQIUMOUJQEBuwSoYWqXZ1mtUcO0LLL222XTJ/tMy2iRGqbJVNUt2XFhGc2QTKxnqXYFb2y+ozbTiQrPcKfxhgQNn4ucE828nJ2bXTCwS+VS+Vf5V/mufLeMpdDU5qK+RQuEpbkdvCnzMrwdPFq7UgnI+dvBk+tdNtpq3E4+X/tSCUt1O7qvXzY3ffKVgS/jpoapL5ESoYap27FCmBaMD8K0IEBOhwDClDUAgVIJIExLxWutcYSpNZSlN4QwLR2xlQ5cFaZpt23rTXOUZEy4vTtas3J8OiIsI7d9R28bb9weHu4irjfsmZ7U7c/JnBW+nTZi6leaepPq/1+ZuEB+OvSSHDjqxflNd6I7g0fkYyAs52tUBhvuJNxCnlgLc0QG+gc6HXrPn4cw9WcJIEz9iRXC1O1YIUwLxgdhWhAgp0MAYcoagECpBBCmpeK11jjC1BrK0htCmJaO2EoHacJ0fGZeOOr6k7FsybTalEo4KlnZdE7s/LTMTXOO2jW87i91S3Z8U5yocDS3ekc33wmyLANJGReW+rGUmpfx28ZHBxYnTv8f939R3nj8SXLi6lPrxkP/LQggTP1ZHghTf2KFMHU7VghTC/GhhqkFiDQBgZDAyqNH5cDBMZmtNwGCeECgawhQw7RrQslEHCFADdNigShjh3AlInV2Zor8NJmbLuwQbkTk/K3ezTt6t9pwR50zMjAqQwPhpjopt43Hdx6f34V8tFjwSjq76k2fSpoGzULAGQLUMHUmFJkGQg3TTJhqOQhhagE7wtQCRJqAAMKUNQCBUgggTEvBSqM9TMBnYTo9Ox3s5B2/7TsUjkm1KZtu754aa5wf3AIeydCM3h4e+X7S1L6cmRC14U+dX33Sp7MkjXAcSagzqetUpuwQrmRlNNsyXvPSZG7qHcHDrMzglvH5DM065+9q3whTVyPDuHwlgDD1K3IIU3fjhTC1EBuEqQWINAEBhClrAAKlEECYloKVRnuYQBFhmrRDuJGX6bd6q+zJMINSC8v5DXf0OQk1L9Pa8mWHcCUck2pTKmHJDuHd9+JDmHZfTJlRvQQQpvXyz9s7wjQvseqOR5gWZE0N04IAOR0CMQJ7dt8tl125XoZHkutcAQwCEMhHgBqm+XjVdTQ1TKshb2OH8OX/MiIvH/MzeXn4lUBWquzMhGzLpN3Gk3YIr2bmQS82dgjX2ZMDRmjOZ042ZVQ6sEO4q5s+VRlvX/qihqkfkaKGqR9xUqOkhqk/saKGqduxQpgWjA/CtCBATocAwpQ1AIFSCSBMS8VrrfFeEaZJO4Tr+pOhcLS9Q/j8LeZBVqaNHcJ/S35L/kH+QZ6VZ3PHX+0QHuzw3byrd/y27egt3NHbxtvtED60KKht2XTbeERe9tIO4QjT3MuzthMQprWhz9UxwjQXrloPRpjWij9X5wjTXLgqPxhhWhA5wrQgQE6HAMKUNQCBUgkgTEvFa63xqoRp3h3CTf1JJR57YYfw4VA4NgnL2A7hR/1zvwycMCpDxywJ5WSz/NTn6s15wjqW4flpO4RbW0Q01EQAYerPgkCY+hErhKkfcVKjRJj6EyuEqduxQphaiA81TC1ApAkIhARWHj0qBw6OyewcSCAAARsEqGFqg6K9NoKdvIPNclSdyfEp8/2YTM9My5Hpw2G25YR+vtWGO5O6dmVYz3JmXCZnpiTYgTysaxk5/7Wp1+xNomBLS4eWhpvghJvhhNmWI4tGZdjs/h2RlE3yMeVWb7UruBKSg4tUBme5wrJIDdOC6DgdAl1JgBqmXRlWJlUjAWqY1gi/g66pYdoBtIpOQZhaAI0wtQCRJiCAMGUNQKAUAgjThVhVHUmTaamyJo28VKIxLjGNkNTHTU8EQjIUncGGOxEBGh4TtDcmqoalad+0Y+OW8E4Xitoh3NSfDCSkEotmN+95ednpDuHD/e1vNe907C6dhzB1KRqMpRsIIEy7IYrMwSUCCFOXotF+LAjT9ozqOgJhaoE8wtQCRJqAAMKUNQCBUgi4LkxfnVAb54zL2EwgJAMBOSaTM5NN2ZL6dvBQaAbCMvjZyMnGuWoX8dn5c80xRoYediDT0mRYBjt+D4v+fyCoO6l3B9ciM2EznYHR5tqU0WP0+Qvlp2l36dCyUtZXLzaKMO3FqDPnMgkgTMukS9u9SABh6lfUEabuxgthWjA21DAtCJDTIRAjsGf33XLZletleGQxbCAAAQsE2tUwVdmW0SzJaPZkU73LMHtSP6Zv9Z6QsakjQUam2SE8dou5kZhjCeeqNuraMby/rz+sMWlu3VbCMtj1O7gtPJSYoZDUj6nnBlVWZiAtjdQMZGcgKqPnRo8x56r/VZZn0ldVNUwtLKmeb+IrD31eznrb+XLsqhN6noXLAKhh6nJ0msdGDVM/YkUNUz/ipEZJDVN/YkUNU7djhTAtGB+EaUGAnA4BhClrAAILCKg6lkoyqixK/f+M+v+Irnep/9c/B8+b75W8HGv8fERLTPXYL/74LPnmsm/JT+XFoCZmWNdSZWy6UNdyyWC4cc5gKCW1pAzEZJB9GYpI9X14zGiYkdm4rTzMzgx2Bp8/N8iunBeg6nlXMy0Rpv68ESBM/YgVwtSPOKlRIkz9iBXC1I84IUz9iZMaKcLU7XghTAvGB2FaECCnQwBhyhrwgICqO6lE5bykDMXl9FiipIxKTC0uleRUx+rbyiP/hxmajeenxmRydsIqkWvlWvm8fF4OyIHEdlXGYzQDsiEYdRalybwMZGY08zIQmWHG5WCYqRnKyUBkzp+rjwuPma+dOSIq05OvgADC1J+VgDD1I1YIUz/ihDD1J04IU39iRYapP7FCmLodK4SphfhQw9QCRJqAQEhg5dGjcuDgmMzOgQQC7QmoDMlAQAYZldHv0yTlkanDzceq+plTaiOfCdHPqexNXR8zaPfQ5KvtB2L5CCUtlWBUu36r71VGpcquVOJRPzYwKqP6Z/XcYlk8uDj1+VUrVsjERJ8M9s/XxRwNbzFfNrzc8shpDgLdT4Aapt0fY2ZYLQFqmFbLm966nwA1TP2KMTVM3Y0XwtRCbBCmFiDSBAQQpl2xBlRNykYmZiSTMi4h1a3huv5lmH2pz4lIymCn8ta3o0/NTFbGTGVhakmpxGXT/0peBtmUgciMCE0tNRdHzokdq+SnEp3mHCU/w1vObU7M9U2fbM6VtiBQBQGEaRWU6aOXCCBMeynazLUKAgjTKijb6wNhao+l7ZYQphaIIkwtQKQJCCBMS1sDM3MzkVvB5yWlzqAMa14m3S5uMjHjmZqNW8xjmZhKclYpMdXt3DoTUwlHJR4Tsi/jmZjq2MWDRy0Qm1HhmSRE1W3kvn4hTH2NHON2lQDC1NXIMC5fCSBMfY0c43aVAMLU1cgkjwth6m68EKYFY0MN04IAOR0CMQJ7dt8tl125XoZHFnc1GyMxo7Urkzb2MRv36IxMnYWpNv4Zk7EZdcu4qqmpMjGPyJjeqXz+1vSo5JyenaqMpZGY5hZy/X9423jyLeZB1qbOxIzcYj4vQdNvRx/qH65sXj539Nd775ML3/1+Wb7i9T5Po+vHTg1Tf0JMDVM/YkUNUz/ipEbJpk9+xIoapn7ESY2SGqb+xIoapm7HCmFaMD4I04IAOR0CDgnT6dnptjuSmx3Lm24XV3UvE3Yvn9/dPNztPCI5q5SYi/oW6VvGdQZmwi3gjezKyPNKcupMzIQ6mnEJGhWiSEz3XtIIU/dikjQihKkfcVKjRJj6ESuEqR9xQpj6EyeEqT+xQpj6EyuEqduxQpgWjA/CtCBATodAG2GqxKLOolTZlEo4mrqW+ucwu7LpeVUbM9y4R+1I3qYOpnre7GiuhGlVX0pimnqYSbeANzb2SaiZqTMxTc3MFhLUbBY0uGioqmnRj4MEEKYOBiVhSAhTP+KEMPUnTghTf2JFhqkfsUKY+hEnNUqEqT+xQpi6HSuEqYX4UMPUAkSa8IqAqlOpN+VJuQW8uR5mWDNT3TaudjKPbfYT3wxoam5CXptQu5grkTkuVUrMgf6B+TqYLXYkj9fMbMrabLExUFSCDvQPehVzBusvAWqY+hs7Ru4mAWqYuhkXRuUvAWqY+hs7Ru4mAWqYuhmXtFFRw9TdeCFMLcQGYWoBIk0UJqAkZuN28chmPkEty3BH8rAOptnsJ9jNPHg+UXLq58KdyiP/q/qbVX0pidnqFnAtIXWNzFEZGQz/17uNx28nH5ERVSdzkdrRPNy5XB2vNw1SP48IErOqqNJPlQQQplXSpq9eIIAw7YUoM8cqCSBMq6RNX71AAGHqV5QRpu7GC2FqITYIUwsQu7SJydmJxNqW5hbwpNvF9c7kavfx8HZys9lPdEfzuARVz1UrMQe1YNSiUQnHqKiMi0slNFUdTCUwE8Rm87mL5fjXv06OHOmXof5AfiphyhcEINA5AYRp5+w4EwJJBBCmrAsI2CWAMLXLk9YggDD1aw0gTN2NF8K0YGyoYVoQYA2nK4k5pupWmtvJI7eIB7eHhxmZ+vZx9XPw2PwGPmqX8vFQagbfB7eaz0tOc+zs3GxlM1TZkaYmZnDr94gsrIMZk5wq41JnaC7cDKhps59IzUzVh6q/WdbXnt13y2VXrpfhkcVldUG7EOgpAtQw9SPc1DD1I05qlGz65EesqGHqR5zUKKlh6kesqGHqR5zUKKlh6k+sqGHqdqwQphni88C+h+Xm7bv1ke+7aI3csmm9jI4Em6ggTDMAzHDIhK6HGYrHhFvAk24X15mY5pxQbAZ1LwPJmSZEq5SYarMdLSO1wFS3iaudylXmZCg01Y7l+hbx+cf0zuQt6mBGj40K0TIlZoYQWjsEYWoNJQ1BQBNAmPqxEBCmfsQJYepPnBCm/sQKYepHrBCmfsQJYepPnNRIEaZuxwth2iY+jz3xtNy+c6/ce9t1smL5Urlj5159xvUb13W9MFUSU4tHnTmZXAczLjkbmZgtNgNKqok5J3OVvVKUxIzXxEyrg2kkp3o+LjnjEjQuOdXz/X39lc2rWzpCmHZLJJmHKwQQpq5EovU4EKZ+xAlh6k+cEKb+xAph6kesEKZ+xAlh6k+cEKbuxwph2iZGSpCedPxKuXTtBfrIuEBVj1VZw1TJy1a3gGuxGblFXAnNI5OqJmYkezNFgkZ3PFffVykxh/qHgzqXCzbwSb9dXGdiJm32Y2pmxtoymZ5ITLffmFYePSoHDo7JbHUO3W0gjA4CBQlQw7QgQE6HQIwANUxZEhCwS4AapnZ50hoEqGHq1xqghqm78UKYtojN2PikbN2xW9ace1pDmD7z3PPyiW275FNbNsjqE1dpqfjTVw/JTw8dSrwFvFHbMn6L+OTCmphpdTDNZj8q47NKiTmsbxWP3C7eqHMZ7kTedNt4cIt5cDv5SCAy1e3mpo6m/lntRh7eeh6rmdknfe6+ShhZpQQQppXiprMeIIAw7YEgM8VKCSBMK8VNZz1AAGHaA0FmipUSQJhWirtwZwjTwghLawBhmkGYXn7xhXLe2afoI+PC9G/u+xt554/eWVqAaBgCEIAABCAAAQhAAAIQgAAEIAABCECguwiMXDAio7882l2T6qLZIEwzCNNWGaZq06eP9X+si5YEU4FAvQR2ze6SD/d9WEb7+MVRbyTovVsI7JndI+/pe48c03dMt0ypK+fx3Nxz8t/n/rv8Rv9vdOX8umlSfzX3V3KOnCMn9J3QTdPqurm8KC/Kl2e/LFf2X9l1c+u2CX157stygpwgp/QFCSp8uUngiByR/zb73+R3+3/XzQEyqgaBr899XZbIEnlb39ug4jiBT89+WrZu3er4KHt3eAjTNrFvV8NUCdPrN90kh8ame3cVMXMIWCTApk8WYdIUBESETZ/8WAZs+uRHnNQov/LQ5+Wst50vx65CmLocNTZ9cjk6zWNj0yc/YsWmT37ESY3y8Uf3y5Ily+Qtp5/jz6B7dKS77t6OMHU49gjTNsGJb/KkBKr6un7jOv0/wtTh1c3QvCSAMPUybAzaYQIIU4eDExkawtSPOCFM/YkTwtSfWCFM/YgVwtSPOCFM/YmTGinC1O14IUwzxOeBfQ/Lzdt36yPfd9EauWXTehkdGWqceejIFBmmGThyCASyEGDTpyyUOAYC2Qmw6VN2VhwJgSwE2PQpCyWOgUB2Amz6lJ0VR0IgCwE2fcpCyZ1j2PTJnVjER4IwtRAbhKkFiDQBgZAAwpSlAAG7BBCmdnnSGgQQpqwBCNglgDC1y5PWIIAw9WsNIEzdjRfC1EJsEKYWINIEBBCmrAEIlEIAYVoKVhrtYQII0x4OPlMvhQDCtBSsNNrDBBCmfgUfYepuvBCmBWNDDdOCADkdAjEC1DBlSUDALgFqmNrlWVZr1DAti6z9dtn0yT7TMlqkhmkZVMtpkxqm5XC13So1TG0TLa89Nn0qj63tlqlhapuo3fYQpgV5IkwLAuR0CCBMWQMQKJUAwrRUvNYaR5haQ1l6QwjT0hFb6QBhagVjJY0gTCvBXLgThGlhhJU1gDCtDHXhjhCmhRGW2gDCtCBeLUxv/ISo2/L5ggAEihMgw7Q4Q1qAQJQAwtSP9YAw9SNOapQIUz9ihTD1I05qlAhTP2KFMPUjTmqUCFN/YoUwdTtWCFO348PoIAABCEAAAhCAAAQgAAEIQAAC/397dxcrR1XAAfw8EVCaSjGlNKkQiBFIqigBGz+wSojaiiY1XCG8YEmt8ERp2vRCCCEEb9OmwBO1aajwBJZINE2vTzWVxORGQkQbaR9KgaBQSAAJKB9P5kycdXe6d3Z37s6dj/3dJ+3OnI/fOcPM/nfODAECBAgQWEQBgekiYquKAAECBAgQIECAAAECBAgQIECAAIF6CwhM6z0+WkeAAAECBAgQIECAAAECBAgQIECAwCIKCEwXEVtVBAgQIECAAAECBAgQIECAAAECBAjUW0BgWnB83n3v/XD7jofDseOnkhIef2RHuPrKywqWZjcC7RX48KNPwn27D4TDR+aSTj6wfWPYsO7aeTuct332uIuFXHjB+WHfrq3h0otWthdRzwjkCBQ9Hz2072C4eNWK3OMRPIFJEXhm9tlw764DSXfXX7cm3L9tYzjn7LNyu//Sq6+H3Y8+FWbu3hTOW7qks208th57crZn30Hnvklx1s/JFBj1WjB7DDl+JnPe6HWvwCjnqedeOBFuvXNnp4Dsea27rHSj225eF+7aPIWdAIEuAYFpgemQnvTXXHVF8kUzXjDfM7M/PDi9SWhTwNMu7RaIF73xL56A02Bn6+apeX9gyNt+mP3bral3BHoFipyPui+SfQk1owiEEL9Y7tl3MOzduSUJPrvPQ/18un+kWH35JZ390m0H7c+cwKQJjHItGM9re5/4bfjpTd9Pjsf4PWvz9j1hZnqTm1MmbeLob0dg1PNUvNZbtXJ5csyk14orli/rBKLx87nnXxzqx0HDQGCSBQSmBUY/e0dB9gtrgSLtQqCVAvFL5fQv9odtd9zU+TEh74vkoO0Fpq2cJjq1AIGFnI/cYboAeLu2SiB7LGS/mM7X2bw7TNMfClsFpTMECggMurYbVKTvWYOEfD4JAkXPU6lNNiAVmE7CrNHHcQgITAso9ruQdjdBAUi7tF6g393XeSfoQdt/9PHHPY/CsBy/9VNIBwcILOR8JDA1vQiEzp036aqhaDLsyqFhl+S7k9tMm2SBQdd2gx594cfySZ49+h4F+v1oMOx5KhXMZhXZJfmW45trBPoLCEwLzIz4BfXpQ0d7bmEXmBaAtEvrBfp9mRwUmGafB5e3ffzs4KGjZyyHbD2sDhL4n8BCzkcCU9OIwP+/iN54w9rOct9hv4jOF5h2u1pObJZNusCo14JZL9+xJn0G6X8amBY5T0W9Qasm0h8lpm5Y67n2phuBjIDAtMCUWMgdPQWqswuBxgqMelfBqNv3W+bVWCwNJ1BAYCHnI4FpAXC7tE5gIXfuDBOYRjDHWuumjQ6NIDDqtV130fHYOf3WO56zOIK3TdsnsJDzVLxOnJ7ZP/AFufEmlFdeO+2lT+2bPnq0QAGBaQHAhTwzrkB1diHQWIFRn1s1ju0bi6XhBAoILOR8JMQpAG6XVgoUfTacwLSV00Gnxiww6rVdWr2wdMwDobhGCxQ5Tw0blkYYgWmjp4fGlyggMC2AW+StxAWqsQuBVggM89b77iUgedvHE3/8i298TE/u3vDYimmiEwUFBp2P8pZZCUwLotutdQKD3j48X3DTLzCNx9zskblwy4brE6dhl/e3DlWHCHQJ5F3bxc2yx5hl+KYPgV6BQeep7GPK8pbhx2vH3xz+Y/jx+m+F+Axhzwk22wjMLyAwLTg70v+wHDt+Kinh8Ud2dEKcgkXajUArBdJA5/CRuaR/3S+/6Bfm5G2fPgvujTffTspaffklnl/aylmjU6MI5J2P+h1j2Qf9e3naKNq2batA93Gx/ro1ZzynvntZcPaYiybpCzOy5zDXiG2dMfo1ikDetV02MM2+4DOtJ3tcjlK/bQm0QSDvPJUNTOOPDo89OdvT7e7rveznXk7YhhmiD2UICEzLUFUmAQIECBAgQIAAAQIECBAgQIAAAQKNFBCYNnLYNJoAAQIECBAgQIAAAQIECBAgQIAAgTIEBKZlqCqTAAECBAgQIECAAAECBAgQIECAAIFGCghMGzlsGk2AAAECBAgQIECAAAECBAgQIECAQBkCAtMyVJVJgAABAgQIECBAgAABAgQIECBAgEAjBQSmjRw2jSZAgAABAgQIECBAgAABAgQIECBAoAwBgWkZqsokQIAAAQIECBAgQIAAAQIECBAgQKCRAgLTRg6bRhMgQIAAAQIECBAgQIAAAQIECBAgUIaAwLQMVWUSIECAAAECBAgQIECAAAECBAgQINBIAYFpI4dNowkQIECAAAECBAgQIECAAAECBAgQKENAYFqGqjIJECBAgAABAgQIECBAgAABAgQIEGikgMC0kcOm0QQIECBAgAABAgQIECBAgAABAgQIlCEgMC1DVZkECBAgQIAAAQIECBAgQIAAAQIECDRSQGDayGHTaAIECBAgQIAAAQIECBAgQIAAAQIEyhAQmJahqkwCBAgQIECAAAECBAgQIECAAAECBBopIDBt5LBpNAECBAgQIECAAAECBAgQIECAAAECZQgITMtQVSYBAgQIECBAgAABAgQIECBAgAABAo0UEJg2ctg0mgABAgQIECBAgAABAgQIECBAgACBMgQEpmWoKpMAAQIECBAgMEaBZ2aLTMKNAAAIpElEQVSfDffuOtBT4urLLwl7d24JJ1/+Z7j1zp3h8Ud2hKuvvKxnm4f2HQx/fuFEst15S5eEvHLe+df7YfP2PeGNN9+et+UPbN8YVq1cntTX7y9tw3MvnEi2WX/dmnD/to3hnLPP6mye91nc6N333g+373g4HDt+at523HbzunDxqhU9JrFtG9ZdG1569fWkH59dtrTT77Sgfp+l7cnrzxiHUlEECBAgQIAAAQINEBCYNmCQNJEAAQIECBCYXIFs6JlKxH//5le/mISkMQg9eOhoT0AYw8F7ZvaHB6c3hUsvWhmGKadbOZY59/yLfQPP6Zn9Yd+urUm5/f7SEPLCC87v2e7Djz4J9+0+EA4fmesbps5X1p59B88IP+O2sZ5sW9JQNAa/2RA5Gjz25GxIw+YYIvcrY3Jnm54TIECAAAECBAhEAYGpeUCAAAECBAgQqKlAerfl1A1rk7sn5/tLg8gVy5eFuzZPhfT/r7nqimS/YcsZZ2AaQ87vfP3L4YN/f5i0KQ04nz50NJx77qfCBx/854wwdlyBaQyKf3D918LfXnypU0cMUnc/+lT4yurPhz/86S+dAFZgWtPJr1kECBAgQIAAgQoFBKYV4quaAAECBAgQIJAnkA1C87ZN76ycmd4UXnv9rZ47TkcpJ61joXeYxsB06+apsOeXv+65yzUupX/ltdPh9FvvlBqYbv35T0LahvQu3FhvrL/7blyBqWOQAAECBAgQIEAgKyAwNScIECBAgAABAjUWyD5js3s5ebbZ6ZLz+O/Z5eijlBP3zwtMh3mGabqM/ldP/T5p5o++943kDs+ZuzeF+G9lB6bxUQR//fvJ5LECW352Y7j/oSfCtjtuSv4tG5gO6k+Np4emESBAgAABAgQIlCAgMC0BVZEECBAgQIAAgXELdD//M5adfT5o/Ld06f01V17WWQafbccw5QwKTId5hmkamMaXScUl8vFv6offTh4REIPdxQhMl31mSfICqaVLPh2+cOmqxCT7vFd3mI57piqPAAECBAgQINB8AYFp88dQDwgQIECAAIEJE5hviX322aWDWPKW6o9jSf7enVtCfLFS9oVTixWYxpdSxX48+sTvOi+fEpgOmhU+J0CAAAECBAgQEJiaAwQIECBAgACBmgrEO0aPHT8Vrl3zpTNaGEPH+Je+UCn+7/kC01HLiWWNMzCNz1c9+fI/wnfXXpO0eTED09j32SNz4ZYN1yd1C0xrOtk1iwABAgQIECBQIwGBaY0GQ1MIECBAgAABAt0C6RL7z61c3vOCpPR5pNnnlOYFpnFp+rDljDswzY7qYgam2boFpo4xAgQIECBAgACBQQIC00FCPidAgAABAgQIVCiQhqbxTtP0r9/zS+NneUvyRylnUGA66CVJMdBNn2Eal+SXEZjG4PPeXQc6RT+wfWPyfNR4N2t8Zmp86VNckj9MYDqoPxUOv6oJECBAgAABAgQqEBCYVoCuSgIECBAgQIAAAQIECBAgQIAAAQIE6ikgMK3nuGgVAQIECBAgQIAAAQIECBAgQIAAAQIVCAhMK0BXJQECBAgQIECAAAECBAgQIECAAAEC9RQQmNZzXLSKAAECBAgQIECAAAECBAgQIECAAIEKBASmFaCrkgABAgQIECBAgAABAgQIECBAgACBegoITOs5LlpFgAABAgQIECBAgAABAgQIECBAgEAFAgLTCtBVSYAAAQIECBAgQIAAAQIECBAgQIBAPQUEpvUcF60iQIAAAQIECBAgQIAAAQIECBAgQKACAYFpBeiqJECAAAECBAgQIECAAAECBAgQIECgngIC03qOi1YRIECAAAECBAgQIECAAAECBAgQIFCBgMC0AnRVEiBAgAABAgQIECBAgAABAgQIECBQTwGBaT3HRasIECBAgAABAgQIECBAgAABAgQIEKhAQGBaAboqCRAgQIAAAQIECBAgQIAAAQIECBCop4DAtJ7jolUECBAgQIAAAQIECBAgQIAAAQIECFQgIDCtAF2VBAgQIECAAAECBAgQIECAAAECBAjUU0BgWs9x0SoCBAgQIECAAAECBAgQIECAAAECBCoQEJhWgK5KAgQIECBAgAABAgQIECBAgAABAgTqKSAwree4aBUBAgQIECBAgAABAgQIECBAgAABAhUICEwrQFclAQIECBAgQIAAAQIECBAgQIAAAQL1FBCY1nNctIoAAQIECBAgQIAAAQIECBAgQIAAgQoEBKYVoKuSAAECBAgQIECAAAECBAgQIECAAIF6CghM6zkuWkWAAAECBAgQIECAAAECBAgQIECAQAUCAtMK0FVJgAABAgQIECBAgAABAgQIECBAgEA9BQSm9RwXrSJAgAABAgQIECBAgAABAgQIECBAoAIBgWkF6KokQIAAAQIECBAgQIAAAQIECBAgQKCeAgLTeo6LVhEgQIAAAQIECBAgQIAAAQIECBAgUIGAwLQCdFUSIECAAAECBAgQIECAAAECBAgQIFBPAYFpPcdFqwgQIECAAAECBAgQIECAAAECBAgQqEBAYFoBuioJECBAgAABAgQIECBAgAABAgQIEKingMC0nuOiVQQIECBAgAABAgQIECBAgAABAgQIVCAgMK0AXZUECBAgQIAAAQIECBAgQIAAAQIECNRTQGBaz3HRKgIECBAgQIAAAQIECBAgQIAAAQIEKhAQmFaArkoCBAgQIECAAAECBAgQIECAAAECBOopIDCt57hoFQECBAgQIECAAAECBAgQIECAAAECFQgITCtAVyUBAgQIECBAgAABAgQIECBAgAABAvUUEJjWc1y0igABAgQIECBAgAABAgQIECBAgACBCgQEphWgq5IAAQIECBAgQIAAAQIECBAgQIAAgXoKCEzrOS5aRYAAAQIECBAgQIAAAQIECBAgQIBABQIC0wrQVUmAAAECBAgQIECAAAECBAgQIECAQD0F/guvyvrw4xnYfAAAAABJRU5ErkJggg==",
"text/html": [
"
| \n", " | SYSTEM TIME | \n", "A | \n", "B | \n", "E | \n", "caption | \n", "
|---|---|---|---|---|---|
| 0 | \n", "0.000000 | \n", "20.000000 | \n", "0.000000 | \n", "0.0 | \n", "Initial state | \n", "
| 1 | \n", "0.050000 | \n", "19.000000 | \n", "1.000000 | \n", "0.0 | \n", "\n", " |
| 2 | \n", "0.100000 | \n", "18.060000 | \n", "1.940000 | \n", "0.0 | \n", "\n", " |
| 3 | \n", "0.125000 | \n", "17.618200 | \n", "2.381800 | \n", "0.0 | \n", "\n", " |
| 4 | \n", "0.150000 | \n", "17.189654 | \n", "2.810346 | \n", "0.0 | \n", "\n", " |
| 5 | \n", "0.175000 | \n", "16.773964 | \n", "3.226036 | \n", "0.0 | \n", "\n", " |
| 6 | \n", "0.200000 | \n", "16.370745 | \n", "3.629255 | \n", "0.0 | \n", "\n", " |
| 7 | \n", "0.225000 | \n", "15.979623 | \n", "4.020377 | \n", "0.0 | \n", "\n", " |
| 8 | \n", "0.250000 | \n", "15.600234 | \n", "4.399766 | \n", "0.0 | \n", "\n", " |
| 9 | \n", "0.275000 | \n", "15.232227 | \n", "4.767773 | \n", "0.0 | \n", "\n", " |
| 10 | \n", "0.275000 | \n", "15.232227 | \n", "4.767773 | \n", "30.0 | \n", "Set concentration of `E` | \n", "
| 11 | \n", "0.275200 | \n", "14.372651 | \n", "5.627349 | \n", "30.0 | \n", "\n", " |
| 12 | \n", "0.275400 | \n", "13.575171 | \n", "6.424829 | \n", "30.0 | \n", "\n", " |
| 13 | \n", "0.275600 | \n", "12.835300 | \n", "7.164700 | \n", "30.0 | \n", "\n", " |
| 14 | \n", "0.275800 | \n", "12.148878 | \n", "7.851122 | \n", "30.0 | \n", "\n", " |
| 15 | \n", "0.276000 | \n", "11.512043 | \n", "8.487957 | \n", "30.0 | \n", "\n", " |
| 16 | \n", "0.276200 | \n", "10.921213 | \n", "9.078787 | \n", "30.0 | \n", "\n", " |
| 17 | \n", "0.276440 | \n", "10.263435 | \n", "9.736565 | \n", "30.0 | \n", "\n", " |
| 18 | \n", "0.276728 | \n", "9.542527 | \n", "10.457473 | \n", "30.0 | \n", "\n", " |
| 19 | \n", "0.277074 | \n", "8.767428 | \n", "11.232572 | \n", "30.0 | \n", "\n", " |
| 20 | \n", "0.277419 | \n", "8.089086 | \n", "11.910914 | \n", "30.0 | \n", "\n", " |
| 21 | \n", "0.277834 | \n", "7.376689 | \n", "12.623311 | \n", "30.0 | \n", "\n", " |
| 22 | \n", "0.278249 | \n", "6.771007 | \n", "13.228993 | \n", "30.0 | \n", "\n", " |
| 23 | \n", "0.278663 | \n", "6.256054 | \n", "13.743946 | \n", "30.0 | \n", "\n", " |
| 24 | \n", "0.279161 | \n", "5.730676 | \n", "14.269324 | \n", "30.0 | \n", "\n", " |
| 25 | \n", "0.279659 | \n", "5.299738 | \n", "14.700262 | \n", "30.0 | \n", "\n", " |
| 26 | \n", "0.280256 | \n", "4.875569 | \n", "15.124431 | \n", "30.0 | \n", "\n", " |
| 27 | \n", "0.280853 | \n", "4.542897 | \n", "15.457103 | \n", "30.0 | \n", "\n", " |
| 28 | \n", "0.281570 | \n", "4.229802 | \n", "15.770198 | \n", "30.0 | \n", "\n", " |
| 29 | \n", "0.282430 | \n", "3.951341 | \n", "16.048659 | \n", "30.0 | \n", "\n", " |
| 30 | \n", "0.283462 | \n", "3.720982 | \n", "16.279018 | \n", "30.0 | \n", "\n", " |
| 31 | \n", "0.284700 | \n", "3.547589 | \n", "16.452411 | \n", "30.0 | \n", "\n", " |
| 32 | \n", "0.286186 | \n", "3.432586 | \n", "16.567414 | \n", "30.0 | \n", "\n", " |
| 33 | \n", "0.287969 | \n", "3.368656 | \n", "16.631344 | \n", "30.0 | \n", "\n", " |
| 34 | \n", "0.290109 | \n", "3.341352 | \n", "16.658648 | \n", "30.0 | \n", "\n", " |
| 35 | \n", "0.292677 | \n", "3.333912 | \n", "16.666088 | \n", "30.0 | \n", "\n", " |
| 36 | \n", "0.295758 | \n", "3.333265 | \n", "16.666735 | \n", "30.0 | \n", "\n", " |
| 37 | \n", "0.299456 | \n", "3.333353 | \n", "16.666647 | \n", "30.0 | \n", "\n", " |
| 38 | \n", "0.303893 | \n", "3.333317 | \n", "16.666683 | \n", "30.0 | \n", "\n", " |
| 39 | \n", "0.309218 | \n", "3.333343 | \n", "16.666657 | \n", "30.0 | \n", "\n", " |
| 40 | \n", "0.315608 | \n", "3.333314 | \n", "16.666686 | \n", "30.0 | \n", "\n", " |