{ "cells": [ { "cell_type": "markdown", "id": "5cbc8640", "metadata": {}, "source": [ "### A MINIMALIST, \"get-started\", demonstration of the unimolecular elementary reaction `A <-> B`,\n", "with 1st-order kinetics in both directions, taken to equilibrium.\n", "\n", "**\"No frills!\"** Please see other experiments for advanced graphics, analysis, diagnostics, fine-tuning, etc." ] }, { "cell_type": "markdown", "id": "5a3fe1d4-ffc9-4db9-ac0f-d51d2231d32b", "metadata": {}, "source": [ "### TAGS : \"quick-start\", \"uniform compartment\"" ] }, { "cell_type": "code", "execution_count": 1, "id": "97a57e9a-039b-479a-81dc-81399e22743a", "metadata": {}, "outputs": [], "source": [ "LAST_REVISED = \"Aug. 29, 2025\"\n", "LIFE123_VERSION = \"1.0.0rc6\" # Library version this experiment is based on" ] }, { "cell_type": "code", "execution_count": 2, "id": "b5b8a8b0-d417-4432-b6a8-c196af57b105", "metadata": {}, "outputs": [], "source": [ "#import set_path # Using MyBinder? Uncomment this before running the next cell!\n", " # Importing this module will add the project's home directory to sys.path" ] }, { "cell_type": "code", "execution_count": 3, "id": "a29db1c7", "metadata": { "tags": [] }, "outputs": [], "source": [ "#import sys, os\n", "#os.getcwd()\n", "#sys.path.append(\"C:/some_path/my_env_or_install\") # CHANGE to the folder containing your venv or libraries installation!\n", "# NOTE: If any of the imports below can't find a module, uncomment the lines above, or try: import set_path\n", "\n", "import life123" ] }, { "cell_type": "code", "execution_count": 4, "id": "ccd6701b-ae96-4d20-b537-40c2c40df9aa", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "OK\n" ] } ], "source": [ "life123.check_version(LIFE123_VERSION) # To check compatibility" ] }, { "cell_type": "code", "execution_count": null, "id": "ac9eea69-174c-43e5-9eed-443cbc5e2ba7", "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "markdown", "id": "e0529a0c", "metadata": {}, "source": [ "## Initialize the System" ] }, { "cell_type": "code", "execution_count": 5, "id": "78077d8c", "metadata": { "tags": [] }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Number of reactions: 1\n", "0: A <-> B (Elementary Unimolecular reaction) (kF = 3 / kR = 2 / delta_G = -1,005.1 / K = 1.5 / Temp = 25 C)\n", "Chemicals involved in the above reactions: ['A', 'B']\n" ] } ], "source": [ "# Instantiate the simulator and specify the chemicals\n", "uc = life123.UniformCompartment() \n", "\n", "# Elementary Reaction A <-> B\n", "uc.add_reaction(reactants=\"A\", products=\"B\", \n", " kF=3., kR=2.)\n", "\n", "uc.describe_reactions()" ] }, { "cell_type": "code", "execution_count": 6, "id": "9fc3948d", "metadata": {}, "outputs": [], "source": [ "# Set the initial concentrations of all the chemicals\n", "uc.set_conc({\"A\": 80., \"B\": 10.})" ] }, { "cell_type": "markdown", "id": "987af2c5", "metadata": { "tags": [] }, "source": [ "## Run the reaction" ] }, { "cell_type": "code", "execution_count": 7, "id": "43735178-313b-48cf-a583-5181238feac3", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "47 total variable step(s) taken in 0.100 sec\n", "Number of step re-do's because of elective soft aborts: 3\n", "Norm usage: {'norm_A': 26, 'norm_B': 22, 'norm_C': 22, 'norm_D': 22}\n", "System Time is now: 1.1343\n" ] } ], "source": [ "uc.single_compartment_react(initial_step=0.1, target_end_time=1.) # Using defaults for all other parameters" ] }, { "cell_type": "code", "execution_count": 8, "id": "2d5df59c", "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
SYSTEM TIMEABstepcaption
00.00000080.00000010.000000Set concentration
10.00640078.59200011.40800011st reaction step
20.00960077.91052812.0894722
30.01344077.10584612.8941543
40.01804876.15876713.8412334
50.02357875.04845814.9515425
60.02910773.96884616.0311546
70.03463772.91908317.0809177
80.04016671.89834418.1016568
90.04569670.90582719.0941739
100.05233269.74773520.25226510
110.05896768.62806721.37193311
120.06560367.54554622.45445412
130.07223866.49894023.50106013
140.07887465.48705824.51294214
150.08550964.50874925.49125115
160.09347263.37372626.62627416
170.10143462.28389327.71610717
180.10939761.23744928.76255118
190.11736060.23266829.76733219
200.12532259.26788930.73211120
210.13487758.15624931.84375121
220.14443357.09771732.90228322
230.15398856.08975833.91024223
240.16354355.12995534.87004524
250.17500954.03321835.96678225
260.18647552.99935737.00064326
270.19794152.02476937.97523127
280.21170150.92231239.07768828
290.22546049.89570040.10430029
300.23922048.93971741.06028330
310.25573147.87145942.12854131
320.27224246.89139343.10860732
330.29205645.81240744.18759333
340.31186944.84031445.15968634
350.33564643.78936546.21063535
360.35942242.86335547.13664536
370.38795341.88424548.11575537
380.42219140.87692749.12307338
390.45642940.04205049.95795039
400.49751439.21170450.78829640
410.54681738.41997951.58002141
420.60598037.70411352.29588742
430.67697537.09919152.90080943
440.76217036.63096553.36903544
450.86440436.30843653.69156445
460.98708436.11924153.88075946
471.13430036.03147053.96853047last reaction step
\n", "
" ], "text/plain": [ " SYSTEM TIME A B step caption\n", "0 0.000000 80.000000 10.000000 Set concentration\n", "1 0.006400 78.592000 11.408000 1 1st reaction step\n", "2 0.009600 77.910528 12.089472 2 \n", "3 0.013440 77.105846 12.894154 3 \n", "4 0.018048 76.158767 13.841233 4 \n", "5 0.023578 75.048458 14.951542 5 \n", "6 0.029107 73.968846 16.031154 6 \n", "7 0.034637 72.919083 17.080917 7 \n", "8 0.040166 71.898344 18.101656 8 \n", "9 0.045696 70.905827 19.094173 9 \n", "10 0.052332 69.747735 20.252265 10 \n", "11 0.058967 68.628067 21.371933 11 \n", "12 0.065603 67.545546 22.454454 12 \n", "13 0.072238 66.498940 23.501060 13 \n", "14 0.078874 65.487058 24.512942 14 \n", "15 0.085509 64.508749 25.491251 15 \n", "16 0.093472 63.373726 26.626274 16 \n", "17 0.101434 62.283893 27.716107 17 \n", "18 0.109397 61.237449 28.762551 18 \n", "19 0.117360 60.232668 29.767332 19 \n", "20 0.125322 59.267889 30.732111 20 \n", "21 0.134877 58.156249 31.843751 21 \n", "22 0.144433 57.097717 32.902283 22 \n", "23 0.153988 56.089758 33.910242 23 \n", "24 0.163543 55.129955 34.870045 24 \n", "25 0.175009 54.033218 35.966782 25 \n", "26 0.186475 52.999357 37.000643 26 \n", "27 0.197941 52.024769 37.975231 27 \n", "28 0.211701 50.922312 39.077688 28 \n", "29 0.225460 49.895700 40.104300 29 \n", "30 0.239220 48.939717 41.060283 30 \n", "31 0.255731 47.871459 42.128541 31 \n", "32 0.272242 46.891393 43.108607 32 \n", "33 0.292056 45.812407 44.187593 33 \n", "34 0.311869 44.840314 45.159686 34 \n", "35 0.335646 43.789365 46.210635 35 \n", "36 0.359422 42.863355 47.136645 36 \n", "37 0.387953 41.884245 48.115755 37 \n", "38 0.422191 40.876927 49.123073 38 \n", "39 0.456429 40.042050 49.957950 39 \n", "40 0.497514 39.211704 50.788296 40 \n", "41 0.546817 38.419979 51.580021 41 \n", "42 0.605980 37.704113 52.295887 42 \n", "43 0.676975 37.099191 52.900809 43 \n", "44 0.762170 36.630965 53.369035 44 \n", "45 0.864404 36.308436 53.691564 45 \n", "46 0.987084 36.119241 53.880759 46 \n", "47 1.134300 36.031470 53.968530 47 last reaction step" ] }, "execution_count": 8, "metadata": {}, "output_type": "execute_result" } ], "source": [ "uc.get_history() # The system's history, saved during the run of single_compartment_react()" ] }, { "cell_type": "code", "execution_count": null, "id": "6409d11e-26a0-4785-bf66-88566aac56ea", "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "markdown", "id": "03866901", "metadata": { "tags": [] }, "source": [ "## Plots changes of concentration with time \n", "Notice that adaptive variable time steps were automatically taken" ] }, { "cell_type": "code", "execution_count": 9, "id": "6b033cc7-078f-4b94-a466-1aa91e2fca4d", "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}", "legendgroup": "A", "line": { "color": "darkturquoise", "dash": "solid", "shape": "linear" }, "marker": { "symbol": "circle" }, "mode": "lines", "name": "A", "orientation": "v", "showlegend": true, "type": "scatter", "x": [ 0, 0.006400000000000002, 0.009600000000000003, 0.013440000000000004, 0.018048000000000005, 0.023577600000000004, 0.029107200000000007, 0.03463680000000001, 0.04016640000000001, 0.045696000000000014, 0.05233152000000001, 0.05896704000000001, 0.06560256000000002, 0.07223808000000002, 0.07887360000000003, 0.08550912000000004, 0.09347174400000004, 0.10143436800000004, 0.10939699200000004, 0.11735961600000004, 0.12532224000000003, 0.13487738880000003, 0.14443253760000002, 0.15398768640000002, 0.16354283520000001, 0.17500901376000003, 0.18647519232000004, 0.19794137088000005, 0.21170078515200005, 0.22546019942400006, 0.23921961369600006, 0.25573091082240007, 0.27224220794880005, 0.29205576450048004, 0.31186932105216003, 0.335645588914176, 0.359421856776192, 0.38795337821061116, 0.42219120393191417, 0.4564290296532172, 0.4975144205187808, 0.5468168895574572, 0.6059798524038689, 0.6769754078195629, 0.7621700743183957, 0.8644036741169949, 0.9870839938753141, 1.1343003775852971 ], "xaxis": "x", "y": [ 80, 78.592, 77.910528, 77.1058458624, 76.15876717373031, 75.04845757891101, 73.96884582376929, 72.91908317443371, 71.89834436282696, 70.90582693788352, 69.7477353740692, 68.62806650892247, 67.54554556951605, 66.49894007682887, 65.48705844253587, 64.50874861235279, 63.373726382799354, 62.2838929294738, 61.237449146205506, 60.23266755485373, 59.26788945357223, 58.15624872361806, 57.09771745659816, 56.0897583084074, 55.12995515844305, 54.033217749985546, 52.999357276322066, 52.02476894664434, 50.92231177389454, 49.89570042592975, 48.93971693212988, 47.871459377140354, 46.89139341164059, 45.81240721519991, 44.84031378886753, 43.78936544572559, 42.863355249162076, 41.88424541214457, 40.87692656753266, 40.042049758158775, 39.21170378709872, 38.41997915447455, 37.704113470448206, 37.099191058819024, 36.630964980446095, 36.30843587395681, 36.11924081574702, 36.031469807322495 ], "yaxis": "y" }, { "hovertemplate": "Chemical=B
SYSTEM TIME=%{x}
Concentration=%{y}", "legendgroup": "B", "line": { "color": "green", "dash": "solid", "shape": "linear" }, "marker": { "symbol": "circle" }, "mode": "lines", "name": "B", "orientation": "v", "showlegend": true, "type": "scatter", "x": [ 0, 0.006400000000000002, 0.009600000000000003, 0.013440000000000004, 0.018048000000000005, 0.023577600000000004, 0.029107200000000007, 0.03463680000000001, 0.04016640000000001, 0.045696000000000014, 0.05233152000000001, 0.05896704000000001, 0.06560256000000002, 0.07223808000000002, 0.07887360000000003, 0.08550912000000004, 0.09347174400000004, 0.10143436800000004, 0.10939699200000004, 0.11735961600000004, 0.12532224000000003, 0.13487738880000003, 0.14443253760000002, 0.15398768640000002, 0.16354283520000001, 0.17500901376000003, 0.18647519232000004, 0.19794137088000005, 0.21170078515200005, 0.22546019942400006, 0.23921961369600006, 0.25573091082240007, 0.27224220794880005, 0.29205576450048004, 0.31186932105216003, 0.335645588914176, 0.359421856776192, 0.38795337821061116, 0.42219120393191417, 0.4564290296532172, 0.4975144205187808, 0.5468168895574572, 0.6059798524038689, 0.6769754078195629, 0.7621700743183957, 0.8644036741169949, 0.9870839938753141, 1.1343003775852971 ], "xaxis": "x", "y": [ 10, 11.408000000000001, 12.089472, 12.894154137600001, 13.841232826269698, 14.951542421088993, 16.031154176230725, 17.080916825566298, 18.101655637173042, 19.094173062116482, 20.252264625930806, 21.371933491077524, 22.45445443048395, 23.501059923171127, 24.512941557464124, 25.491251387647203, 26.62627361720064, 27.716107070526196, 28.762550853794487, 29.767332445146263, 30.732110546427762, 31.84375127638193, 32.90228254340183, 33.910241691592596, 34.87004484155694, 35.96678225001445, 37.00064272367793, 37.97523105335565, 39.07768822610545, 40.10429957407024, 41.06028306787011, 42.12854062285964, 43.1086065883594, 44.18759278480008, 45.15968621113246, 46.2106345542744, 47.13664475083792, 48.11575458785542, 49.123073432467336, 49.95795024184122, 50.78829621290127, 51.580020845525446, 52.29588652955179, 52.90080894118097, 53.3690350195539, 53.69156412604318, 53.88075918425297, 53.9685301926775 ], "yaxis": "y" } ], "layout": { "autosize": true, "legend": { "title": { "text": "Chemical" }, "tracegroupgap": 0 }, "shapes": [ { "line": { "color": "gray", "dash": "dot", "width": 1 }, "type": "line", "x0": 0, "x1": 0, "xref": "x", "y0": 0, "y1": 1, "yref": "y domain" }, { "line": { "color": "gray", "dash": "dot", "width": 1 }, "type": "line", "x0": 0.006400000000000002, "x1": 0.006400000000000002, "xref": "x", "y0": 0, "y1": 1, "yref": "y domain" }, { "line": { "color": "gray", "dash": "dot", "width": 1 }, "type": "line", "x0": 0.009600000000000003, "x1": 0.009600000000000003, "xref": "x", "y0": 0, "y1": 1, "yref": "y domain" }, { "line": { "color": "gray", "dash": "dot", "width": 1 }, "type": "line", "x0": 0.013440000000000004, "x1": 0.013440000000000004, "xref": "x", "y0": 0, "y1": 1, "yref": "y domain" }, { "line": { "color": "gray", "dash": "dot", "width": 1 }, "type": "line", "x0": 0.018048000000000005, "x1": 0.018048000000000005, "xref": "x", "y0": 0, "y1": 1, "yref": "y domain" }, { "line": { "color": "gray", "dash": "dot", "width": 1 }, "type": "line", "x0": 0.023577600000000004, "x1": 0.023577600000000004, "xref": "x", "y0": 0, "y1": 1, "yref": "y domain" }, { "line": { "color": "gray", "dash": "dot", "width": 1 }, "type": "line", "x0": 0.029107200000000007, "x1": 0.029107200000000007, "xref": "x", "y0": 0, "y1": 1, "yref": "y domain" }, { "line": { "color": "gray", "dash": "dot", "width": 1 }, "type": "line", "x0": 0.03463680000000001, "x1": 0.03463680000000001, "xref": "x", "y0": 0, "y1": 1, "yref": "y domain" }, { "line": { "color": "gray", "dash": "dot", "width": 1 }, "type": "line", "x0": 0.04016640000000001, "x1": 0.04016640000000001, "xref": "x", "y0": 0, "y1": 1, "yref": "y domain" }, { "line": { "color": "gray", "dash": "dot", "width": 1 }, "type": "line", "x0": 0.045696000000000014, "x1": 0.045696000000000014, "xref": "x", "y0": 0, "y1": 1, "yref": "y domain" }, { "line": { "color": "gray", "dash": "dot", "width": 1 }, "type": "line", "x0": 0.05233152000000001, "x1": 0.05233152000000001, "xref": "x", "y0": 0, "y1": 1, "yref": "y domain" }, { "line": { "color": "gray", "dash": "dot", "width": 1 }, "type": "line", "x0": 0.05896704000000001, "x1": 0.05896704000000001, "xref": "x", "y0": 0, "y1": 1, "yref": "y domain" }, { "line": { "color": "gray", "dash": "dot", "width": 1 }, "type": "line", "x0": 0.06560256000000002, "x1": 0.06560256000000002, "xref": "x", "y0": 0, "y1": 1, "yref": "y domain" }, { "line": { "color": "gray", "dash": "dot", "width": 1 }, "type": "line", "x0": 0.07223808000000002, "x1": 0.07223808000000002, "xref": "x", "y0": 0, "y1": 1, "yref": "y domain" }, { "line": { "color": "gray", "dash": "dot", "width": 1 }, "type": "line", "x0": 0.07887360000000003, "x1": 0.07887360000000003, "xref": "x", "y0": 0, "y1": 1, "yref": "y domain" }, { "line": { "color": "gray", "dash": "dot", "width": 1 }, "type": "line", "x0": 0.08550912000000004, "x1": 0.08550912000000004, "xref": "x", "y0": 0, "y1": 1, "yref": "y domain" }, { "line": { "color": "gray", "dash": "dot", "width": 1 }, "type": "line", "x0": 0.09347174400000004, "x1": 0.09347174400000004, "xref": "x", "y0": 0, "y1": 1, "yref": "y domain" }, { "line": { "color": "gray", "dash": "dot", "width": 1 }, "type": "line", "x0": 0.10143436800000004, "x1": 0.10143436800000004, "xref": "x", "y0": 0, "y1": 1, "yref": "y domain" }, { "line": { "color": "gray", "dash": "dot", "width": 1 }, "type": "line", "x0": 0.10939699200000004, "x1": 0.10939699200000004, "xref": "x", "y0": 0, "y1": 1, "yref": "y domain" }, { "line": { "color": "gray", "dash": "dot", "width": 1 }, "type": "line", "x0": 0.11735961600000004, "x1": 0.11735961600000004, "xref": "x", "y0": 0, "y1": 1, "yref": "y domain" }, { "line": { "color": "gray", "dash": "dot", "width": 1 }, "type": "line", "x0": 0.12532224000000003, "x1": 0.12532224000000003, "xref": "x", "y0": 0, "y1": 1, "yref": "y domain" }, { "line": { "color": "gray", "dash": "dot", "width": 1 }, "type": "line", "x0": 0.13487738880000003, "x1": 0.13487738880000003, "xref": "x", "y0": 0, "y1": 1, "yref": "y domain" }, { "line": { "color": "gray", "dash": "dot", "width": 1 }, "type": "line", "x0": 0.14443253760000002, "x1": 0.14443253760000002, "xref": "x", "y0": 0, "y1": 1, "yref": "y domain" }, { "line": { "color": "gray", "dash": "dot", "width": 1 }, "type": "line", "x0": 0.15398768640000002, "x1": 0.15398768640000002, "xref": "x", "y0": 0, "y1": 1, "yref": "y domain" }, { "line": { "color": "gray", "dash": "dot", "width": 1 }, "type": "line", "x0": 0.16354283520000001, "x1": 0.16354283520000001, "xref": "x", "y0": 0, "y1": 1, "yref": "y domain" }, { "line": { "color": "gray", "dash": "dot", "width": 1 }, "type": "line", "x0": 0.17500901376000003, "x1": 0.17500901376000003, "xref": "x", "y0": 0, "y1": 1, "yref": "y domain" }, { "line": { "color": "gray", "dash": "dot", "width": 1 }, "type": "line", "x0": 0.18647519232000004, "x1": 0.18647519232000004, "xref": "x", "y0": 0, "y1": 1, "yref": "y domain" }, { "line": { "color": "gray", "dash": "dot", "width": 1 }, "type": "line", "x0": 0.19794137088000005, "x1": 0.19794137088000005, "xref": "x", "y0": 0, "y1": 1, "yref": "y domain" }, { "line": { "color": "gray", "dash": "dot", "width": 1 }, "type": "line", "x0": 0.21170078515200005, "x1": 0.21170078515200005, "xref": "x", "y0": 0, "y1": 1, "yref": "y domain" }, { "line": { "color": "gray", "dash": "dot", "width": 1 }, "type": "line", "x0": 0.22546019942400006, "x1": 0.22546019942400006, "xref": "x", "y0": 0, "y1": 1, "yref": "y domain" }, { "line": { "color": "gray", "dash": "dot", "width": 1 }, "type": "line", "x0": 0.23921961369600006, "x1": 0.23921961369600006, "xref": "x", "y0": 0, "y1": 1, "yref": "y domain" }, { "line": { "color": "gray", "dash": "dot", "width": 1 }, "type": "line", "x0": 0.25573091082240007, "x1": 0.25573091082240007, "xref": "x", "y0": 0, "y1": 1, "yref": "y domain" }, { "line": { "color": "gray", "dash": "dot", "width": 1 }, "type": "line", "x0": 0.27224220794880005, "x1": 0.27224220794880005, "xref": "x", "y0": 0, "y1": 1, "yref": "y domain" }, { "line": { "color": "gray", "dash": "dot", "width": 1 }, "type": "line", "x0": 0.29205576450048004, "x1": 0.29205576450048004, "xref": "x", "y0": 0, "y1": 1, "yref": "y domain" }, { "line": { "color": "gray", "dash": "dot", "width": 1 }, "type": "line", "x0": 0.31186932105216003, "x1": 0.31186932105216003, "xref": "x", "y0": 0, "y1": 1, "yref": "y domain" }, { "line": { "color": "gray", "dash": "dot", "width": 1 }, "type": "line", "x0": 0.335645588914176, "x1": 0.335645588914176, "xref": "x", "y0": 0, "y1": 1, "yref": "y domain" }, { "line": { "color": "gray", "dash": "dot", "width": 1 }, "type": "line", "x0": 0.359421856776192, "x1": 0.359421856776192, "xref": "x", "y0": 0, "y1": 1, "yref": "y domain" }, { "line": { "color": "gray", "dash": "dot", "width": 1 }, "type": "line", "x0": 0.38795337821061116, "x1": 0.38795337821061116, "xref": "x", "y0": 0, "y1": 1, "yref": "y domain" }, { "line": { "color": "gray", "dash": "dot", "width": 1 }, "type": "line", "x0": 0.42219120393191417, "x1": 0.42219120393191417, "xref": "x", "y0": 0, "y1": 1, "yref": "y domain" }, { "line": { "color": "gray", "dash": "dot", "width": 1 }, "type": "line", "x0": 0.4564290296532172, "x1": 0.4564290296532172, "xref": "x", "y0": 0, "y1": 1, "yref": "y domain" }, { "line": { "color": "gray", "dash": "dot", "width": 1 }, "type": "line", "x0": 0.4975144205187808, "x1": 0.4975144205187808, "xref": "x", "y0": 0, "y1": 1, "yref": "y domain" }, { "line": { "color": "gray", "dash": "dot", "width": 1 }, "type": "line", "x0": 0.5468168895574572, "x1": 0.5468168895574572, "xref": "x", "y0": 0, "y1": 1, "yref": "y domain" }, { "line": { "color": "gray", "dash": "dot", "width": 1 }, "type": "line", "x0": 0.6059798524038689, "x1": 0.6059798524038689, "xref": "x", "y0": 0, "y1": 1, "yref": "y domain" }, { "line": { "color": "gray", "dash": "dot", "width": 1 }, "type": "line", "x0": 0.6769754078195629, "x1": 0.6769754078195629, "xref": "x", "y0": 0, "y1": 1, "yref": "y domain" }, { "line": { "color": "gray", "dash": "dot", "width": 1 }, "type": "line", "x0": 0.7621700743183957, "x1": 0.7621700743183957, "xref": "x", "y0": 0, "y1": 1, "yref": "y domain" }, { "line": { "color": "gray", "dash": "dot", "width": 1 }, "type": "line", "x0": 0.8644036741169949, "x1": 0.8644036741169949, "xref": "x", "y0": 0, "y1": 1, "yref": "y domain" }, { "line": { "color": "gray", "dash": "dot", "width": 1 }, "type": "line", "x0": 0.9870839938753141, "x1": 0.9870839938753141, "xref": "x", "y0": 0, "y1": 1, "yref": "y domain" }, { "line": { "color": "gray", "dash": "dot", "width": 1 }, "type": "line", "x0": 1.1343003775852971, "x1": 1.1343003775852971, "xref": "x", "y0": 0, "y1": 1, "yref": "y domain" } ], "template": { "data": { "bar": [ { "error_x": { "color": "#2a3f5f" }, "error_y": { "color": "#2a3f5f" }, "marker": { "line": { "color": "#E5ECF6", "width": 0.5 }, "pattern": { "fillmode": "overlay", "size": 10, "solidity": 0.2 } }, "type": "bar" } ], "barpolar": [ { "marker": { "line": { "color": "#E5ECF6", "width": 0.5 }, "pattern": { "fillmode": "overlay", "size": 10, "solidity": 0.2 } }, "type": "barpolar" } ], "carpet": [ { "aaxis": { "endlinecolor": "#2a3f5f", "gridcolor": "white", "linecolor": "white", "minorgridcolor": "white", "startlinecolor": "#2a3f5f" }, "baxis": { "endlinecolor": "#2a3f5f", "gridcolor": "white", "linecolor": "white", "minorgridcolor": "white", "startlinecolor": "#2a3f5f" }, "type": "carpet" } ], "choropleth": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "type": "choropleth" } ], "contour": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "colorscale": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ], "type": "contour" } ], "contourcarpet": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "type": "contourcarpet" } ], "heatmap": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "colorscale": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ], "type": "heatmap" } ], "heatmapgl": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "colorscale": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ], "type": "heatmapgl" } ], "histogram": [ { "marker": { "pattern": { "fillmode": "overlay", "size": 10, "solidity": 0.2 } }, "type": "histogram" } ], "histogram2d": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "colorscale": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ], "type": "histogram2d" } ], "histogram2dcontour": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "colorscale": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ], "type": "histogram2dcontour" } ], "mesh3d": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "type": "mesh3d" } ], "parcoords": [ { "line": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "parcoords" } ], "pie": [ { "automargin": true, "type": "pie" } ], "scatter": [ { "fillpattern": { "fillmode": "overlay", "size": 10, "solidity": 0.2 }, "type": "scatter" } ], "scatter3d": [ { "line": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scatter3d" } ], "scattercarpet": [ { "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scattercarpet" } ], "scattergeo": [ { "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scattergeo" } ], "scattergl": [ { "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scattergl" } ], "scattermapbox": [ { "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scattermapbox" } ], "scatterpolar": [ { "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scatterpolar" } ], "scatterpolargl": [ { "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scatterpolargl" } ], "scatterternary": [ { "marker": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "type": "scatterternary" } ], "surface": [ { "colorbar": { "outlinewidth": 0, "ticks": "" }, "colorscale": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ], "type": "surface" } ], "table": [ { "cells": { "fill": { "color": "#EBF0F8" }, "line": { "color": "white" } }, "header": { "fill": { "color": "#C8D4E3" }, "line": { "color": "white" } }, "type": "table" } ] }, "layout": { "annotationdefaults": { "arrowcolor": "#2a3f5f", "arrowhead": 0, "arrowwidth": 1 }, "autotypenumbers": "strict", "coloraxis": { "colorbar": { "outlinewidth": 0, "ticks": "" } }, "colorscale": { "diverging": [ [ 0, "#8e0152" ], [ 0.1, "#c51b7d" ], [ 0.2, "#de77ae" ], [ 0.3, "#f1b6da" ], [ 0.4, "#fde0ef" ], [ 0.5, "#f7f7f7" ], [ 0.6, "#e6f5d0" ], [ 0.7, "#b8e186" ], [ 0.8, "#7fbc41" ], [ 0.9, "#4d9221" ], [ 1, "#276419" ] ], "sequential": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ], "sequentialminus": [ [ 0, "#0d0887" ], [ 0.1111111111111111, "#46039f" ], [ 0.2222222222222222, "#7201a8" ], [ 0.3333333333333333, "#9c179e" ], [ 0.4444444444444444, "#bd3786" ], [ 0.5555555555555556, "#d8576b" ], [ 0.6666666666666666, "#ed7953" ], [ 0.7777777777777778, "#fb9f3a" ], [ 0.8888888888888888, "#fdca26" ], [ 1, "#f0f921" ] ] }, "colorway": [ "#636efa", "#EF553B", "#00cc96", "#ab63fa", "#FFA15A", "#19d3f3", "#FF6692", "#B6E880", "#FF97FF", "#FECB52" ], "font": { "color": "#2a3f5f" }, "geo": { "bgcolor": "white", "lakecolor": "white", "landcolor": "#E5ECF6", "showlakes": true, "showland": true, "subunitcolor": "white" }, "hoverlabel": { "align": "left" }, "hovermode": "closest", "mapbox": { "style": "light" }, "paper_bgcolor": "white", "plot_bgcolor": "#E5ECF6", "polar": { "angularaxis": { "gridcolor": "white", "linecolor": "white", "ticks": "" }, "bgcolor": "#E5ECF6", "radialaxis": { "gridcolor": "white", "linecolor": "white", "ticks": "" } }, "scene": { "xaxis": { "backgroundcolor": "#E5ECF6", "gridcolor": "white", "gridwidth": 2, "linecolor": "white", "showbackground": true, "ticks": "", "zerolinecolor": "white" }, "yaxis": { "backgroundcolor": "#E5ECF6", "gridcolor": "white", "gridwidth": 2, "linecolor": "white", "showbackground": true, "ticks": "", "zerolinecolor": "white" }, "zaxis": { "backgroundcolor": "#E5ECF6", "gridcolor": "white", "gridwidth": 2, "linecolor": "white", "showbackground": true, "ticks": "", "zerolinecolor": "white" } }, "shapedefaults": { "line": { "color": "#2a3f5f" } }, "ternary": { "aaxis": { "gridcolor": "white", "linecolor": "white", "ticks": "" }, "baxis": { "gridcolor": "white", "linecolor": "white", "ticks": "" }, "bgcolor": "#E5ECF6", "caxis": { "gridcolor": "white", "linecolor": "white", "ticks": "" } }, "title": { "x": 0.05 }, "xaxis": { "automargin": true, "gridcolor": "white", "linecolor": "white", "ticks": "", "title": { "standoff": 15 }, "zerolinecolor": "white", "zerolinewidth": 2 }, "yaxis": { "automargin": true, "gridcolor": "white", "linecolor": "white", "ticks": "", "title": { "standoff": 15 }, "zerolinecolor": "white", "zerolinewidth": 2 } } }, "title": { "text": "Reaction `A <-> B` . Changes in concentrations with time (time steps shown in dashed lines)" }, "xaxis": { "anchor": "y", "autorange": true, "domain": [ 0, 1 ], "range": [ -0.0005533172573586816, 1.1348536948426557 ], "title": { "text": "SYSTEM TIME" }, "type": "linear" }, "yaxis": { "anchor": "x", "autorange": true, "domain": [ 0, 1 ], "range": [ 6.111111111111111, 83.88888888888889 ], "title": { "text": "Concentration" }, "type": "linear" } } }, "image/png": "iVBORw0KGgoAAAANSUhEUgAABLkAAAFoCAYAAAC7YSngAAAgAElEQVR4Xuy9CZwU1bn//TD7BjPDLsiwKlEUBBEFMZgQRIyXEJcAIi4kUVwSTXxj8rnmc7ne1+RN4s0/5sZ7lWv8X66oQRJxR4xoIIIsomhAg6gQlUXWGWZh9pn3PNV9eqprqrtOdVVXn+7+FZ+mu6vO8pzvc6q66jfnPKdHp9gIGwiAAAiAAAiAAAiAAAiAAAiAAAiAAAiAAAikMYEeELnS2HswHQRAAARAAARAAARAAARAAARAAARAAARAwCAAkQsdAQRAAARAAARAAARAAARAAARAAARAAARAIO0JQORKexeiASAAAiAAAiAAAiAAAiAAAiAAAiAAAiAAAhC50AdAAARAAARAAARAAARAAARAAARAAARAAATSngBErrR3IRoAAiAAAiAAAiAAAiAAAiAAAiAAAiAAAiAAkQt9AARAAARAAARAAARAAARAAARAAARAAARAIO0JQORKexeiASAAAiAAAiAAAiAAAiAAAiAAAiAAAiAAAhC50AdAAARAAARAAARAAARAAARAAARAAARAAATSngBErrR3IRoAAiAAAiAAAiAAAiAAAiAAAiAAAiAAAiAAkQt9AARAAARAAARAAARAAARAAARAAARAAARAIO0JQORKexeiASAAAiAAAiAAAiAAAiAAAiAAAiAAAiAAAhC50AdAAARAAARAAARAAARAAARAAARAAARAAATSngBErrR3IRoAAiAAAiAAAiAAAiAAAiAAAiAAAiAAAiAAkQt9AARAAARAAARAAARAAARAAARAAARAAARAIO0JQORKexeiASAAAiAAAiAAAiAAAiAAAiAAAiAAAiAAAhC50AdAAARAAARAAARAAARAAARAAARAAARAAATSngBErrR3IRoAAiAAAiAAAiAAAiAAAiAAAiAAAiAAAiAAkQt9AARAAARAAARAAARAAARAAARAAARAAARAIO0JQORKexeiASAAAiAAAiAAAiAAAiAAAiAAAiAAAiAAAhC50AdAAARAAARAAARAAARAAARAAARAAARAAATSngBErrR3IRoAAiAAAiAAAiAAAiAAAiAAAiAAAiAAAiAAkQt9AARAAARAAARAAARAAARAAARAAARAAARAIO0JQORKexeiASAAAiAAAiAAAiAAAiAAAiAAAiAAAiAAAhC50AdAAARAAARAAARAAARAAARAAARAAARAAATSngBErrR3IRoAAiAAAiAAAiAAAiAAAiAAAiAAAiAAAiAAkQt9AARAAARAAARAAARAAARAAARAAARAAARAIO0JQORKexeiASAAAiAAAiAAAiAAAiAAAiAAAiAAAiAAAhC50AdAAARAAARAAARAAARAAARAAARAAARAAATSngBErrR3IRoAAiAAAiAAAiAAAiAAAiAAAiAAAiAAAiAAkQt9AARAAARAAARAAARAAARAAARAAARAAARAIO0JQORKexeiASAAAiAAAiAAAiAAAiAAAiAAAiAAAiAAAhC50AdAAARAAARAAARAAARAAARAAARAAARAAATSngBErrR3IRoAAiAAAiAAAiAAAiAAAiAAAiAAAiAAAiAAkQt9AARAAARAAARAAARAAARAAARAAARAAARAIO0JQORKexeiASAAAiAAAiAAAiAAAiAAAiAAAiAAAiAAAhC50AcSIvCDJQ/Sn9dvo0umTaTf3Ht7QmVke6bP9h+iWQt+bGB4+YlfUtXgAdmOBO0HARAAARAAARAAARAAARAAARAAgYQJQORKGB0y3v/QClr21BoINB67wpiLb0grsXDuzffSF0eO0/pVv/XYcmQHARDQgUC6XYP8ZIY/2PhJM5iycO+ROOcNW3fQzXf/mn5650KaP2d64gXFyBnk/YFqXclusxeIsi+fNXo4PbV0iZeikBcEQAAEQMBEIGtFrmlX3EFHj5/o1hl0GZkU7yZO2v7+umUp78x8k8Gb04+zvMngtMm6uUo1jFh9iu2KN1JLsln6q7to6qSzU9oM9ufOD/dG2dC3d3mUoKV6Y5nShgRU+R+efY3ue2A53TD3UvrRLfMCqjWzqkn2A0i6XEv99Kpdv4zHwavIFc+HUkTS4fpmZRzr/M3EPpNJwlAmtYX7ZJD3dMm+3gZ5f6Bal12bdbouSftwH+HnryDKAgEQyHYCWS1ysfPNo1HkD431oT4VnSTTbrLlDQWz9PsvVpKV2xsE6W+/hE2+UbX2Kf4uhSOdpyTKhz27vi9vwOVDquqNZSrOm6DrhMjlnXiyH7oy7VqqQhwilwqlkLgwsF/vbn+kycQ+k0nCUCa1BSKX2rlql0r1XkR3kYvbJu+Rdb5PTNxTyAkCIAACwROAyGWZciVvnlL9V+dMu4mTIwW4i3MsLz9+yOWDHJeZ6Ogw8wgzr2JXLJEr2Q/xXi8bMjZYvPZzf5x87hhjpJnqjaVXu5A/Owgk+/zItGtpor0iVSO5ErU32fnk74fdb30m9plMalMmtSXZ/dxafrKvt0HeH6jWlew2++VDryNq/bID5YAACIBAJhCAyKUocpmDhEvHW2+OzaKLuXPEEmDMo5tkek677+ARI9aVdZPlxPphlzd+Mp/dqByZd/nv/jkS9JzTJ3P0mvlhguvieBBuR12ZWfghbsW68eP9iYpdsUQuaW+iQlyyLzRu/4Lopg/ZTeH0o19a+zqPDuRplna+s55ndiMJ7c5dFX/Z3Ty74RPPt3Y2WduXjHNepV4VpvKGfdL4M4wpnXIzt8EsMptZyDTmh9kf/dvDkam07Jut2/9uCObWzSygW/mYr7Mcjyaoa6nd74fKdcZupJHdA761H1q/O3FQ8VWsvurGh3JxDcn9Zz/5jvF7IDd5blr7YCxWKv0wlt2qvlfpM5LfoIF9o3675e+c1c5Y1xbr9dLN76Td/YQ8F5z8z21Uucdxc21L9JrKtjjlNZ8DC7/386jQE3Z/QLPro+Z0su1W3nZ9xHr+ubXFrj/a1ePlnDTXYQ1BwL+/HKrD3Aft+g6XYccyXj/jPG76SCxf24XisPMht4U3p/igdr/TdtfRZNju5tov2eoQiiTZ950oHwRAAASSTQAil0Xksnvgt5v6YSdc8L5nX94QNfUhlsBhJ4hwWn5w49UK4/2l0u6GSN7ImH8c4+2zCg3J/AuS1d5YYpBTZ0+GuGWt03wj5eYBg8uJ1a5E2+vEw6/jbu2T/UqlD3HZ/CBrjjVmV5+bMu3O0VhTT+1inVjrtxu9qSr8xRK5WHBT4RPLh3b1880yP9DJG/pknPMq9aow5XbxNYU3O1HL/IAV76/s5odz6x8V2FYWFcyx0OzsT/W11C7eiqrwbffQI/mbr0/WNqo+1Mn+p+qrWP1VxYfmB2bZd82Cs/lh0G6/9Xqs2g9j2cxtjnWNd9tnJD9zeebfK7v95v4s227XbhUx1O5aEOs6ZydcqN7jxLtOm233ck1VyWu+Llj7lXVRFLvy7PZZBWVzfzT7yno/48aWWH0xlsilcv2Mdx8Q6z6T/+hgFblUrqUq/SyR33Iz33j34FZxWPXeRfV66Lftbq/98UaX+nW/h3JAAARAIFsIQOQyiVyxHjxiCUD8Y/zOjo8c/4rEP8QTzj7NEK94Uxlq7+YmO94DhtX2WH+9TtZfkOz+Qmp3gxnvhDPfbKqMrvHj5DWLXap1xgs8r/Kg4ofdiZTBfcRNnDSvfUieZ04PJ9wWa79009djnUOyPzmNjGQ7LzzvLJKjT+zYxhvJZf3rsuo5pjK1wg0HVX+p1KvKVIpcdv3e+jDpViBRuVaYz9lUX0vjMdv41s64K5xZH3rMYghzkAt+WH+LVB/qzCKXiq/8FLnsVmiNFSPLut9NP7Sz2XoNsKZx02di9fVYo4Ps9sc6R+2ulXbtiXWPwu2QInC8Nqne4zjZ6RS3UeWaGq8OeT2O1RY7kSBW26z1WK/PXNbDjz1v4L58xuQIR9X7OTeCRbyRXPK+Ufo91jmi2odVrvNclt05otLP/Pi9sTJW9aGb61O8kVxOv93xGJptd3vtd7ouJXJ/hzwgAAIgkK0Eslrksq6uaDeNSuUBzDq0WP5V19ypzCJCrJsAc3o3N9lu0qregPh1QtjZphIDylx/oiKXnR9UY60lKnKx3XZD59mWZE4J9eKvZItc5r90m+2M99dxmc764OHmQS3eeWYexSHrcCP0SfuSIXL5LYKrnvMq9aoyjfXgz/v9FLliTZWLN8rJ3AdjjcqwG+1iTavK1TyiJ5FpKOa+yj7a9u6H9L1vX2FM85N2MlPzg3i6ilyxmKqyl/3OaRSu5BPr98DNb2q8vh7r4dy6P1Y6JzvNwgffz8Rrd6w2ubnHieUfq3Dn5ZqqkldV5FJpmzyHrKy5rRPPGW0gfvHVTcbvumynud+o2hLvNzoZIlcsX8VionIttRtFam2X6nXR6RzjclnEj+dDlXtpLkf1eui37W6v/RC5vNzJIi8IgAAIRBPIapHLLEjIH1zrX7NjxdkyY5Q3SfLmzFqG9aFO5S9xTjcA5r+CxxshYh1ervoj7teJEm90E9fh5qEvXacrcjtVpyf5xd1NOW4FONU+ZBYnzaKB3V+4VcuM19etD4pOfc98ntrFGXF6UI5186zallg+UhnxlYxzXqVeN0xjPbj7JXJJW8yjtuKNHFURrpLB1Xz+m32uKqpyf5IPfNzmxdfNNkZ/SWFr7uyvGPEVzQ/eqg910h5VX8Xqs27EBC5DVcyS9VnTu+mHdjY7iUdufn+5fFUxy463Xcweq80qo4ntmLgRY+L9ZshzR1Xk4rISvaaq5FUVluKNprIbUS4FZXlOcbt5k+fX5wcOG/EFzfctqrbE45sMkSvWfabduap6LeU2OPUz1d+/WHHAJCf5R8F4PkyVyKVqu9trP0QuN3euSAsCIAAC8QlA5DJNV5Q/XOaHW6ebYTNe1QcFlR9mNzfZbtKq3oD4ceLEe/DxIvokQ+wy/xUz0amF8eJDyPJVhBM/2LspQzX+VKwHTrlfddRVUCKXWRxww8N8I+808i8bR3IxHzlNLh5X1euhW4GE64yVx6vIFdS11O63JhZLs038sG3+o8qBL47RnFlTuz14Z7rI5eXcZs5OD5Nu+gGX50Xkkvn9/G0wC2dSkHEayeV0rWM73Y4OMvdpKY6o1GM9F6x5VYUlt9cW2a/4nOKpinJUthSUD3xxlPicM1//VG2Jd61Mpcjl5lpqbYNdP1O9x5TcnP7QqaPIpWq7nc/jXfudrktu72OQHgRAAASymQBELkvgebsbMZWg7LHib8gH5oH9ekdujOLdQPMxjqERLxaH9SYi3o1cKmNyOYknbqfJWU9UP8QuP8QtaVc8kcuLqJfsC5RkEG9kCds/ZFB/I4C86k1srFExXkSueDe81r6uMv2O2bKd1rgn8c4psz+SIXI5PZjx9SEZ57xKvapM4z34W0cYxLuxd/sAaXcdTvW11K5/xeNjPd/N5ycfkw/YsfZzGjtfxuOgKkjGuha59aHXkVxu+mEsm83TQGP9tqiM/ovnS1XxS+UPX/HYP/X8X6IWYOC0VkaJ+N9aZyw7rXV5uaaq5HVzXXATz0mWy7+FPFXRHM+MpyzyZp4WbMdZMkt1TK5Yv7/Wa0MsO63XUv6u0s9U7w9Uf2Pj3Vurnjeqor/ftru99rvpM8m+L0T5IAACIJDuBCByWUQuecPK7/IGV/7wWP/Syj++5tXO7EQOKZpZBQS7KWLmgJVuf9jltES7FaysKy7aBfxVmarktrM7iVhOIphqffLG1O1fwuWNT6Ijt6z2xRK5VEQk1baa08k+4kesL9nH7cqyCr+qN4J2wp5ZmEwkJhe3346z7P92U4U5Vo31r8Xc9yaNPyMy7cuaRrVvJkPk4jbanc9sk3mhi2Sc8yr1yv4Qj2m8B3+7aTSxptbEepi1i+tnHlVgvhak+loq+5K5v7sVvWV8QbvVxexiMdn1y3gcvIpc8rw0/zFHXqvsfOhV5JL1OZ3b8a6r8R6Q3fYZVTFL2mNNL+uz+0MDnxvLf/fPMRfAiJXXep2M1ybVexx5fTCf+3ajlO2uEarXVJW8bkQumdZuMQrrqDLzH73M91Lma4s1jxtbYvXHZIzkcrpGSh5O6eS1VLWfqd4fMAu7ayPv5/4oVxo3pzP7ROZVuf/xW+RStd3ttT8Z9+GJ3FsiDwiAAAhkAgGIXDYil/mGxnojYHV6rGXNZTr+UWYhzO7mX94wyrTWG1xrPDCn1eDkzZYsz+7H380NiJcOrvIQp/MUvkTaHi9OjFsBTqV+P0UuWZ9dG6z90k0fsvZhFqFYXOK4JomKXOYHXGk3l8WBuO04x4qfYX5Qs1ukwGkaBdedLJGLy7aez7zPalMyznmVelWZqq7YZ41LJPPFG7FjzcPXOxYDeEqftR+k+lpqvdazL91M27ITF8x9JNbDulUUi8XBD5HLjQ/9ELnMD5rW66XKues0asJNn/Eqcrm5/tr9NpjFmVj3E7w/Vpv4mNV/shzzuST9Zl20xy5mWKLXVK7XKa9bYcmOj90oPVm3ndjodA5ay3PqX2Y/JkPksvMpXyN/9pPvGL+VdrEMzfeOdtdSlX7m5v7Ark+af9N55LjcrL85bP+zL28guz/aWs+RZIhcqra7ufarzBqxO/+xDwRAAARAoDuBrBW50BlAAAQyhwBiWWSOL9ESEAiKQKwRhEHVn271qE4PS7d2wV4QSDUB1ZGOqbYT9YMACIBAuhCAyJUunoKdIAACBgHzVEOJBA9f6BwgAAJuCWTaaGK37XebHtdZt8SQHgScCeA65MwIKUAABEDALQGIXG6JIT0IgEBKCdhNlYsXND+lxqJyEAABrQnI64lfsRm1bqxH4yByeQSI7CBgIWBe7EBl1WIABAEQAAEQUCMAkUuNE1KBAAiAAAiAAAiAAAiAAAiAAAiAAAiAAAhoTAAil8bOgWkgAAIgAAIgAAIgAAIgAAIgAAIgAAIgAAJqBCByqXFCKhAAARAAARAAARAAARAAARAAARAAARAAAY0JQOTS2DkwDQRAAARAAARAAARAAARAAARAAARAAARAQI0ARC41TkgFAiAAAiAAAiAAAiAAAiAAAiAAAiAAAiCgMQGIXBo7B6aBAAiAAAiAAAiAAAiAAAiAAAiAAAiAAAioEYDIpcYJqUAABEAABEAABEAABEAABEAABEAABEAABDQmAJFLY+fANBAAARAAARAAARAAARAAARAAARAAARAAATUCELnUOCEVCIAACIAACIAACIAACIAACIAACIAACICAxgQgcmnsHJgGAiAAAiAAAiAAAiAAAiAAAiAAAiAAAiCgRgAilxonpAIBEAABEAABEAABEAABEAABEAABEAABENCYAEQujZ0D00AABEAABEAABEAABEAABEAABEAABEAABNQIQORS44RUIAACIAACIAACIAACIAACIAACIAACIAACGhOAyKWxc2AaCIAACIAACIAACIAACIAACIAACIAACICAGgGIXGqckAoEQAAEQAAEQAAEQAAEQAAEQAAEQAAEQEBjAhC5NHYOTAMBEAABEAABEAABEAABEAABEAABEAABEFAjAJFLjRNSgQAIgAAIgAAIgAAIgAAIgAAIgAAIgAAIaEwAIpfGzoFpIAACIAACIAACIAACIAACIAACIAACIAACagQgcqlxQioQAAEQAAEQAAEQAAEQAAEQAAEQAAEQAAGNCUDk0tg5MA0EQAAEQAAEQAAEQAAEQAAEQAAEQAAEQECNAEQuNU5IBQIgAAIgAAIgAAIgAAIgAAIgAAIgAAIgoDEBiFwaOwemgQAIgAAIgAAIgAAIgAAIgAAIgAAIgAAIqBGAyKXGCalAAARAAARAAARAAARAAARAAARAAARAAAQ0JgCRS2PnwDQQAAEQAAEQAAEQAAEQAAEQAAEQAAEQAAE1AhC51DghFQiAAAiAAAiAAAiAAAiAAAiAAAiAAAiAgMYEIHJp7ByYBgIgAAIgAAIgAAIgAAIgAAIgAAIgAAIgoEYAIpcaJ6QCARAAARAAARAAARAAARAAARAAARAAARDQmABELo2dA9NAAARAAARAAARAAARAAARAAARAAARAAATUCEDkUuOEVCAAAiAAAiAAAiAAAiAAAiAAAiAAAiAAAhoTgMilsXNgGgiAAAiAAAiAAAiAAAiAAAiAAAiAAAiAgBoBiFxqnJAKBEAABEAABEAABEAABEAABEAABEAABEBAYwIQuTR2DkwDARAAARAAARAAARAAARAAARAAARAAARBQIwCRS40TUoEACIAACIAACIAACIAACIAACIAACIAACGhMACKXxs6BaSAAAiAAAiAAAiAAAiAAAiAAAiAAAiAAAmoEIHKpcUIqEAABEAABEAABEAABEAABEAABEAABEAABjQlA5NLYOTANBEAABEAABEAABEAABEAABEAABEAABEBAjQBELjVOSAUCIAACIAACIAACIAACIAACIAACIAACIKAxAYhcGjsHpoEACIAACIAACIAACIAACIAACIAACIAACKgRgMilxgmpQAAEQAAEQAAEQAAEQAAEQAAEQAAEQAAENCYAkUtj58A0EAABEAABEAABEAABEAABEAABEAABEAABNQIQudQ4IRUIgAAIgAAIgAAIgAAIgAAIgAAIgAAIgIDGBCByaewcmAYCIAACIAACIAACIAACIAACIAACIAACIKBGACKXGiekAgEQAAEQAAEQAAEQAAEQAAEQAAEQAAEQ0JgARC6NnQPTQAAEQAAEQAAEQAAEQAAEQAAEQAAEQAAE1AhA5FLjhFQgAAIgAAIgAAIgAAIgAAIgAAIgAAIgAAIaE4DIpbFzYBoIgAAIgAAIgAAIgAAIgAAIgAAIgAAIgIAaAYhcapyQCgRAAARAAARAAARAAARAAARAAARAAARAQGMCELl8cE7dyVaqa2zzoSQUAQKpJVCYn0Nlxfl0rLY5tYagdhDwiUCvknzq6OykelyjfSKKYlJNoH9FER2va6a29s5Um4L6QcAzgZweRP0ri+mL442ey0IBIKADgZLCXCrIz6Wa+hYdzIENSSQwqE9xEktH0V4IQOTyQi+cFyKXDxBRhBYEIHJp4QYY4SMBiFw+wkRRWhCAyKWFG2CETwQgcvkEEsVoQwAilzauSLohELmSjjjhCiByJYyuKyNELh8goggtCEDk0sINMMJHAhC5fISJorQgAJFLCzfACJ8IQOTyCSSK0YYARC5tXJF0QyByJR1xwhVA5EoYXSjjvffeSz/80T9HTVc8dHA/bXlzHc2+cgHt+WgX7d2zm6bPnG2kX/HYUvr6nHnUs1e58f2VF5+mM8aMo6rhoyKW+L1vzQt/ojFjJ9CQoSOiWut2/+rnVtK4CZNo8JBhUeX4tZ8L3bJxHRWXlNDY8ZO6ecbK0pyg9kQNvfz8Spq78CZbjy5b+gAtuPFWyi8o6Hb8pWdW0PhJU2jQ4CrbvPV1tfTCqidp/vWLY/aW7ds2UVtbG513wUUx06x8/BGaeflVVF5RaZvG3G+cuqVTWTL/O1s3Ek9oOXfShXGLlPyuW7RYebqiXT+NV8mBfZ/S9m2bRf+f69S8yPE/PvEozbhsDlVU9lHOIxMuf/RBuvqaRVRUXOI67yMP/oq+e/vdrvOZMzSebKCnVyyjaxfd5qmcY0cP07q1q+nKeTd4Ksea2a3/VCo/fOggbfrrWvrG1QtVkiecZud7b1NdbQ1Nvmi6Yxl+i1yq55SjYQkk2LzhdSot60VnnzMxgdzJyeLGF8mxwH2p1t9h9yWkNoedyPXBju1UffwoXThtRmqNS1Htify+pMhUbard99le2vHuNpo1++qU2pRNIteTyx6i2VddS2VlPVPKHJW7I9De3k7/+98P0KJb7lLKCJFLCVPaJ+JnhSVLlqR9OzK1ARC5hGc3bN1BN9/964iPb5h7Kf3olnmR7z9Y8iD9ef024/tZo4fTU0u7OjRELiKIXCTEG4hcThfJRB5CIHIRQeTq3rPcCCsQuZzOTG/H3fjCW03+5YbI5R9LXUpK5PdFF9tTZQdEruDJQ+QKnrkfNULk8oNi5pUBkUtvn2a9yPXZ/kM0a8GPaemv7qKpk84m6/c/PPsaPfzY87R+1W8NT869+V6aeM7oiAgGkQsiF/cLiFzOF7pEHkIgckHksutZboQViFzO56aXFG584aUeP/NC5PKTph5lJfL7ooflqbMCIlfw7CFyBc/cjxohcvlBMfPKgMilt0+zXuSSo7jeX7cs4qlpV9xBi6+bTfPnTO8mallFL86EmFx6d3JYp04AMbnUWSFlehDwW+RKj1bDykwmgJhcmezd7GtbNk1XzD7vZmeLMV0xe/yOmFz6+jrrRS52DY/O2vnhXnr5iV+KkVyH6Z5f/D4ycssseHFaO1EMIpe+HRyWuSMAkcsdL6TWnwBELv19BAvdEYDI5Y4XUutNACKX3v6Bde4JQORyzyxdc+ggct3/0Ap68dVNEe0iXVmy3WMuvoF+eudCY6CR1w0ilyAoO8fR4ycMnuaYXFbYUuRiQaxq8ABaW1dP9x86TP/Wpz+dWVjo1R/IDwIpJZCb24MK8nKpsbktpXakS+UdHZ2Uw3fo2LQlUJifKxZf6KSW1g5tbYRhIOCGQGlRHjW2tBNff7CBQLoT6CF+QkuL8qm+sTXtmoJ7gLRzWSAG5+flUK64N2wS12lsmU2gZ0l+0hsoQylZK5JaRCpFLh4MNOHs0+g3997uCweIXL5gDBXCopV55BbvY8BS6HIaycUxuf519jcpR+T7VnkF/Uu/fpRz6BD9dd1rNG/B9fThrg/o490f0tdFGt4eXfqfdNW8BVQu0vL27NMraey48TRi1GmRVvm975k/raBzJpxHw0eMjCLndv/TK5+kiZMm09Bhw6PK8Ws/F8rcSkpKRT0XdPOylaU5QU1NNa1a+QdadNOttr3jwQfup5tuvYMKbFZX/OOKx+mCKRfRkKqhtvDq4twAACAASURBVHlra0/QyieX03cWxz6Bt2zaKFZXbKULL7o4Zu/8v488RN+8ai5VVva2TXNg/75Iv3Hq4k5lyfybNr4hPnbS5Au/HLdIye+mxbdRQX4OnWx2/mG266fxKvns03/Q1s0b6aq5C5yaFzm+7NGlNHvOldS7T1/lPDLhQw/+hm5YdLOxWqfb7f/86mf0w7vvccwW7wb3ZEMDPbbsEVp8252O5cRLcOTwIVqz+gVaeMN3PJVjzezWfyqVf3HwAL2+9hW6ZuGNKskTTrP97beI++xXpl/iWAaPTuwUWkBLmz8il+o55WhYAgnWvf4q9RIr806Y2H312QSK8yWLG1/4UqEPhVh/h30oMtAiDJFL/CHCrHG9u/1tOn70CH11xqWB2qJLZYn8vuhie6rs+MfePfTOti10xdXzU2WCUS//mai0OE+IXOn3xzW3ItcjD/2O5l17PfXs2SulzFG5OwIck4ufI+646ydKGfPFH4xzc3MgcinRSt9E/KyQ7NUVWcBa9tSabqObWL/43aOrjMXwUily+e09iFw+EuWOse3dD6NWTOTpi4MG9jFUSWugeWtMLha5nhM3CNubmwyrCsWfpL7T1Exfem87XSGWCd7z0S7au2c3TZ852zhuDXj7yotP0xljxlHV8FGRVvm9b80Lf6IxYyfQkKEjosi53e/XKoqxymHjtmxcZ4gSY8d3f4izsjQ3pvZEDb38/Eqau/Am296xbOkDtODGWynfRuR66ZkVNH7SFBo0uMo2b31dLb2w6kmaf/3imD0PgeedT8pEAgMj8DwCz9v1LDfBzv2ervjO1o1CNiY6d9KFzp3e5xSbN7xOpWW96OxzJvpccuLFufFF4rX4mxOB5/3lqUNpify+6GB3Km1A4Png6SPwfPDM/agRgef9oJh5ZSQ78LwcweU0fU+KXBxP/L4HlkdAm+ON8045G00muGTaxMgILHmM6zKXwQvzfX7gcNQ+OYKMyzFrJrJcKcxZ62ENxVw2HzfPnuPvELl8PE+kU62rK8oOpbq64tNiquMva47RrtZmqjp+jGZ+8D71/vo36RIxCmP/no8gcpl8BpGr0rYHHzq4n7a8uY5mX+k80mnl44/QzMuvovIK+7JkBaoP5FIkvG7RYiorzqdjtc2OZ5mdGBsvUyIPIRC5IHLZ9Sk3wgpELsdT2VMCN77wVJGPmSFy+QhTk6IS+X3RxPSUmQGRK3j0ELmCZ+5HjRC5/KCYeWUkW+SSYpFVrLKSlOnMohWLT7zxSC/erOGWeB/PVrt8xmT60S3zIsfPGj08kucHSx6kP6/fRuZyed+BL45F0lhFLpnHbDOnYTtYU9l38IhRn9kmqcHwPohcPp8nVmXRqphKh3G1ZudLM2Tgef7r/qqGWvr1ieP0qZi6xtvA3Fz6YXkfmldWTrk+243iQMBvAgg87zdRlJdqAn6LXKluD+oHAQSeRx/IJAIIPJ9J3kRbmAACz2dPP0hm4HmroBSLqt10RbtBOnNmTY0K6G5OYyeCSX3ELFjZlStnv0mRyixaOfUEs9AGkcuJVgqOW1dXbBUBYJ4UYtcDQuw63B6KMTAqr4B+UtGHZpWUpcBCVAkCagQgcqlxQqr0IQCRK318BUvVCEDkUuOEVOlBACJXevgJVqoTgMilzirdU+oscvHUQClQ8QipWBuncSNymcs1j+SyK8Nap3UqIx+Pt+Cfl/6B1RW90AvntYpcsshGIXb9vq6aHqqtphMdoaDH4wsKaUllPzqvsNiHmlEECPhLACKXvzxRWuoJQORKvQ9ggb8EIHL5yxOlpZYARK7U8kft/hOAyOU/U11LTKbIpRpQPtZILqvIFS+2VxAiFwtiOz/cGxHe2KcYyaVrzxZ2ceD5H/7on6nOtCqMObYSB0vf/cmHtP2CqfSoELxufnUNLZsylSb27kc/FWLXnldeQOB5k38ReB4xucynO2JyISaX3eXfTRwov0Uu1Th3yfjZQuB5f6giJpc/HHUqBTG53HsDMbncM/OaAzG5vBJMTX7E5EoNd91rTXZMLnMw+PlzpnfDweIWx7dSEbnsAsSbC/RD5JKB8mNNV7SLtwWRS+NeriJyydUVj4qpi39c/t/08AVT6GhJKeWIdv3w7W301bHn0riRp0daidUVsbqiU5dH4HknQkTLH32Qrr5mERUVlzgntqTgH67v3n6363zmDI0nG+jpFcvo2kW3eSrn2NHDtG7tarpy3g2eyrFmdrtwgErlhw8dpE1/XUvfuHqhSvKE00DkwuqKCXcekREilxd6euaFyOXeLxC53DPzmgMil1eCqckPkSs13HWvNdkiF7dfxgW3jsJiwUjGCVcRuWR8LWs5LDKtX/VbX6YrSnvf2fGRUabcZOB5q6AlbcJ0RU17uhuRi5vAN9fn/tNV9Lv2VnrmZB3N27KJ3hs2nL486gy6vbySKnJyCSIXRC6n7g6Ry4kQRC4nQhC5nAjZH8dIrmgubgTHxIj7nwsil/9MU10iRC73HoDI5Z6Z1xwQubwSTE1+iFyp4a57rUGIXMxAjrIy8zCveKgicsUqRwpMfozkkvaZF+zjfdJWOdJLpuP9LIjJFR55P1ZX1KzXx4rJ5WTm7tYW+nnNUXq1scFIygLXbb0qaVHPCirq0cMpO46DgO8EEJPLd6QoMMUE/J6umOLmoHoQIMTkQifIJAKIyZVJ3kRbmABicmVPP0hmTK7soZicliLwvA9cExW5ZNXvtTTT/1tzhDY1NRq7Bgix6+d9BtClxaU+WIciQECdAEQudVZImR4EIHKlh59gpToBiFzqrJBSfwIQufT3ESx0RwAilzte6ZwaIpe+3oPI5YNvvIpc0oS/iBFd99Uco12tzcauoXn59D0xsmt+WbkPVqIIEHAmAJHLmRFSpBcBiFzp5S9Y60wAIpczI6RIHwIQudLHV7BUjQBELjVOmZAKIpe+XoTI5dE3icTk+vqcedSzV0i4souLs/y5lbRmyBBa17evkWZ4fj59Z+sWI0B91fBREYtVY3eteeFPNGbsBBoydERUa93uXy3sGjdhEg0eMiyqHL/2c6FYXRGrK5o7F1ZXxOqKdpdoN3Gg/Ba5EJMr2iNufOHx59a37IjJ5RtKbQpCTC73rkBMLvfMvOZATC6vBFOTHzG5UsNd91qDismlOwdd7YPI5dEzyRC5pHi1c8BA+j8njtPfxciua0SA+mOjTqOvnHYmLQyP7ILI1eW82hM19PLzK2nuwptsPbps6QO04MZbKb+goNvxl55ZQeMnTaFBg6ts89bX1dILq56k+dcvjtlbtm/bRG1tbXTeBRfFTOMULP7Qwf205U2IXBC5orsQVlfsfkq5EVYgcnn8kXPI7sYXybVEvXSIXOqs0iUlRC73noLI5Z6Z1xwQubwSTE1+iFyp4a57rRC59PYQRC6P/kmmyCVHba0R0xi3rH6G1g2pot0DT6G+uXl0S68KGrL+dTp7zDjH0V1uR2zFSu/XiK1Y5bArMJILIhdELohcTpdlN8IKRC4nmt6Ou/GFt5r8yw2Ryz+WupQEkcu9JyByuWfmNQdELq8EU5MfIldquOteK0QuvT0EkcsH//gVk8vJlPVNDfRftTW0oemkkbSsRw5d07MX3dyzNw3MzXXKjuMg4EgAMbkcESFBmhHwW+RKs+bD3AwkgJhcGejULG4SYnJlsfMztOmIyZWhjrVpFmJy6etriFw++CYokUuaytMXfyumMa4+WU/tYme+eH2ztBfdXl5JI/O6T8fzoYkoIksIQOTKEkdnUTMhcmWRs7OkqRC5ssTRWdJMiFxZ4ugsaiZEruxxNkQufX0NkcsH3wQtckmT94kYUP9Ze5xWNtRSU2cn9RAHvlZcSreKFRknFRb70DIUkW0EIHJlm8czv70QuTLfx9nWQohc2ebxzG4vRK7M9m82tg4iV/Z4HSKXvr6GyOXRN0HE5GIT4wWZ7zl0OD1SV0P/W3eCLtu0gbYNHUa9xL7bevWmGUL0egWrKxICz18Yt6fLwP3XLVpMZcX5dKy22fHMsOuT8TIlEjMFqytidUW7PuUmDpTfIhdWV4z2iBtfOF5UAkqAmFwBgQ6wmkR+XwI0T8uqEJMreLcgJlfwzP2oETG5/KCYeWUgJpfePoXI5dE/OohcMkB9Y0cHPf78H2nNqafS5v79jZadll9A12/eRF8ZP4mGDRsR1Vq3AekReB6rK5o7EESu+BePxpMN9PSKZXTtots8XWWwumJ3fG6EFYhcnrqfY2Y3vnAsLKAEELkCAh1gNRC53MOGyOWemdccELm8EkxNfohcqeGue63ZLnL94dnX6L4HltMNcy+lH90yTzt3QeTy6BKdRC5uCgsPp4sVF98WItfDIkg9x++6dvOb9OGo0+jS086ghaXlVJqTY7QaIlfI+S89s4LGT5pCgwZX2faG+rpaemHVkzT/eohcELnULxgQudRZuU3pRliByOWWrrv0bnzhruTkpYbIlTy2qSoZIpd78hC53DPzmgMil1eCqckPkSs13HWvNdtFrrk330sTzxlNL766idav+q127oLI5YNLUhWTS8X01xob6KG6atrU1Ggk7yVWZLyuZwXd1LOc+uTmqRSBNFlEADG5ssjZWdJUv0WuLMGGZmpMADG5NHYOTHNNADG5XCNDBs0JICaX5g7y0bxsjsk15uIb6P11y4jfl/7qLpo66WwfyXovCiKXd4aks8glm7ezpYl+W1tNa8SKjB1iZ2GPHnSVWJHxNhGkfmger8+IDQREv8jPUY7JBV4gkA4EIHKlg5dgoxsCELnc0EJa3QlA5NLdQ7DPLQGIXG6JpW/6IESuT1taaW9LS+CQhhbk0/CCAtt6eari1u1/p9/cezv9YMmDRhr+rNMGkcsHb6SDyCWb+Wlbq1iRsZr+JFZkbBYrMvLExVklZfR9IXadVVDkAw0Ukc4EIHKls/dgux0BiFzoF5lGACJXpnk0u9sDkSu7/Z+JrYfIlYletW9TECLXz744TD89+EXgUO8Z2J/uO2Wgbb08VfF7377CGL21YesOuucXv9duyiJELo9dRseYXGeImFwyGD03zy721rH2Nlrx3Ep6sWoo/S0cpH5KYTHdIsSulrUv05ixE2jI0OhA9Qg8j5hc5tMFgefjXzwQk8vjxTVOdjdxoPwWubC6YrRj3PgieT3CXcmIyeWOVzqkRkwu915CTC73zLzmQEwurwRTkx8xuVLDXfdag4rJ9fjxanr02PHAcVzbu5K+3ad3t3o/23+IFn7v51Gi1rQr7qCf/eQ7Wk1ZhMjlscukq8glxa+RZ4+n10QHfkSM7vqivd2gcdPWzTRx3Ln09dPOjKIDkQsiF0Qu9QsGRC51Vm5TuhFWIHK5pesuvRtfuCs5eakhciWPbapKhsjlnjxELvfMvOaAyOWVYGryQ+RKDXfdaw1K5NKNw/0PraBlT63pZtYl0yZqNWURIpfHnpPuIpd5xNZT9bX0H7XHaeqG9bR5xEiqPmUQzRNxu24QgeoHiSD1ELkgckHkUr9gQORSZ+U2pRthBSKXW7ru0rvxhbuSk5caIlfy2KaqZIhc7slD5HLPzGsOiFxeCaYmP0Su1HDXvdZsFbl4quKcWVNp/pzpERfxlMWb7/61EYhelw0ilw+eSKeYXCrNXSuC0y8TMbv+IlZmlNvM4lK6rqycLhbv2DKXAGJyZa5vs7Vlfotc2coR7daHAGJy6eMLWOKdAGJyeWeIEvQigJhcevkjmdYEEZMrmfa7LVuKWS8/8UuqGjwgKjtPWVx83ewo8ctt+X6mh8jlA81ME7kkks9EkPpldSfoKSF41XSEpjJOLCyirxaV0tVlvYzRXdgyiwBErszyJ1pDBJELvSDTCEDkyjSPZnd7IHJlt/8zsfUQuTLRq/ZtyjaRK508C5HLB29lqshlRvNHIXQ9JgSvd1qaIrunFpXQfDGdcU5pTx8ooggdCEDk0sELsMFPAhC5/KSJsnQgAJFLBy/ABr8IQOTyiyTK0YUARC5dPJF8OyByJZ9xojVA5EqUXDhfJsXkkijsVmPkYxyTq49YdfHlinJ6uqGOjodHd924aSPlnTWO/mnUGXSuGOklN7cxvDjflo3rqLikhMaOn9TNM3s+2kV79+ym6TNndztWe6KGXn5+Jc1deJOtR5ctfYAW3Hgr5RcUdDv+0jMraPykKTRocJVt3vq6Wnph1ZM0/3rE5DIDwuqK8S8eiMnl8eIaJ7ubOFB+i1xYXTHaMW58kbwe4a5kxORyxysdUiMml3svISaXe2ZecyAml1eCqcmPmFyp4a57rdkak0t3v0j7IHJ59FS2iVzjJkyiwUOGUWtnJ60RMbtWiGD1VX/5M20YdRrt6defRuUV0LfKetLVYoTXthdXkUxvxhxL/ILItZ+2vLmOZl+5wLFXrnz8EZp5+VVUXlEZN63qA7kUCa9btJjKivPpWG2zow0QuSByWQkcPnSQNv11LX3j6oWO/cdLAjfCCkQuL6Sd87rxhXNpwaSAyBUM5yBrgcjlnjZELvfMvOaAyOWVYGryQ+RKDXfda4XIpbeHIHJ59E+2ilxmbM88+xR9/qUz6MlePenztjbjUK543b5lM404ZyLNHvUlKujRI5IFIpe9MHXoIEQu6+n4xycepRmXzaGKyj6uz9Tljz5IV1+ziIqKS1zn5R+u795+t+t85gwYyeUJX9zMboQViFzJ8wOX7MYXybVEvXSIXOqs0iUlRC73noLI5Z6Z1xwQubwSTE1+iFyp4a57rRC59PaQdiIXR+Y/evyELTWdlqU0G5gNMblUunGnSLSh6ST9QcTvWiNWaGwWo714q8jJpTklZTRXjO4aa5rOqFIm0gRLADG5guWN2pJPwG+RK/kWowYQiE8AMbnQQzKJAGJyZZI30RYmgJhc2dMPEJNLX19rJXLNvfleGjSwD/3m3tv1JWZjGUSu7lBOiHhdHLeLBa8PWrqmvp2RX0jfEoHqrxKCV+9cHu+FTScCELl08gZs8YMARC4/KKIMnQhA5NLJG7DFKwGIXF4JIr9uBCBy6eaR5NkDkSt5bL2WrJXINebiG2jpr+6iqZPO9tquQPND5IqP+2/NTYbY9awQvWo7O4zEeeL11eJSmltWTl8T7/wdW+oJQORKvQ9ggb8EIHL5yxOlpZ4ARK7U+wAW+EcAIpd/LFGSHgQgcunhhyCsgMgVBOXE6oDIlRi3SC7E5AqtuqgSYJ6nL754so4+Wv0svTRipBGonrc+YkTXFSU9heDVi2q3bsLqigg8Hzm/EJOL6NjRw7Ru7Wq6ct4NHq9W0dndLhygUjkCz6tQSjzN5g2vU6m4Tp4tYh3qsiEmV/CesBO5PtixnaqPH6ULp80I3iANakRMLvdOQEwu98y85kBMLq8EU5MfMblSw133WhGTS28PaSVy8XTFObOm0vw50/WmZrIOIpe6yCWxsSg2cNwEerWyklbWn6Av2tsjRK/Z9Xc6t7ySpk+cTP1yo8d37floF+3ds5umz5zdrX/I1QHnLrzJtu8sW/oALbjxVsovKOh2/KVnVtD4SVNo0OAq27z1dbX0wqonaf71i2P2y+3bNlGbCLp/3gUXxUzjtCIiAs93RweRCyKX3QnlRljxeySX6oqlyfgRg8jlD1UEnveHo06lQORy7w2IXO6Zec0BkcsrwdTkh8iVGu661wqRS28PaSVybdi6g+75xe9p/arf6k0NIleUf1RHcplFLvPIr9caG2iFmM64WgSrv+SDnVRfUEhvjjqNphaV0DdEwPrLxSivXjk5BJGrC7uTYCZTqj6QS5HwukWLqaw4n47VdsVRi3Uyuh0JlMhDCEQuiFwQuboIQOTy59YAIpc/HHUqJZHfF53sT4UtELmCpw6RK3jmftQIkcsPiplXRraKXH949jW674HlUQ49a/RwemrpEq2crJXIxTG54m1YXVGrvuO7McfFiC4Wu54VUxrfNwWr54pmiLhd3xAB6y8tLqPiHj18rxsFhgggJhd6QqYR8HskV6bxQXvSjwBicqWfz2BxbAKIyYXekWkEEJMr0zwauz3ZGJNLilxmXYY1nJ/euVCr2XhaiVzpekog8Lz/ntvT2kLPCbHrmYZ6+qStJVLBsLx8GiNGes0UotcMIXjxCC9s/hGAyOUfS5SkBwGIXHr4AVb4RwAil38sUVLqCUDkSr0PYIG/BCBy+ctT59IgcoW8wyGnJp4zmn50yzxt3KWdyMVTFm+++9dRgIJYcdE89K5v7/KoKZM/WPIg/Xn9NsMmu+F4ELmS2595VBeP7npOrM64v70tUlm++DRZTGm8VExpnCkEr4EigD02bwQgcnnjh9z6EYDIpZ9PYJE3AhC5vPFDbr0IQOTSyx+wxjsBiFzeGaZLCUGIXJ+e+JT2Vu8NHMnQiqE0vGJ4t3qtI7k+23+IZi34Mb38xC+pavCAwO2MVaFWIpeEZoYkwSVzCNz9D62gZU+tIbvpkGzTw489HxG9rEolAs8nFnjebjVG7qRbNq6Lu7ri9o930ceTL6I1jfW0W4z2klvvhgZatPlNyr9ynjGlcVR+dIB5BJ6/MO5FBzG5ovHwPPvv3n63pwt148kGenrFMrp20W2eysHqit3xIfA8Vlf0clIhJpcXenrmRUwu935BTC73zLzmQEwurwRTkx8xuVLDXfdag4rJ9bM3fkY/ff2ngeO456J76L6v3hdT5LIegMgVx0XTrriDFl83u9t8TqvQ5LeXeR5pLMdYRS2rLRC5ghW5zKsr7mltNcQufu09foyu3bSR/uNrlxjdY1ReQWiEV0kpTSgoIohcELncXDcgcrmh1ZX28KGDtOmva+kbVy9MrADFXBC5IHIpdhXbZBC5vNDTMy9ELvd+gcjlnpnXHBC5vBJMTX6IXKnhrnutQYlcj//tcXp0+6OB47h27LX07fHfjilymQcHyZl4OgldWo3kYrHJbmqiBJeMwPNyRcejx09EnHjJtIn0m3tvN75bhTerLRC5Uidymc+6vdVH6bXn/0TrLptNG5tOUktnZ+QwT2O86cXnaPQ1N9LUnhWUZzldX3pmBY2fNIUGDa6yvYDU19XSC6uepPnXL455gdm+bRO1tbXReRdcFDON04qIhw7upy1vrqPZVy5wvJA5lSULwOqKi6iouMSRpzUBRC7XyIwMELkS46aaC6srqpKKnw4ilz8cdSoFIpd7b0Dkcs/Maw6IXF4JpiY/RK7UcNe91qBELt042AWeZxtj6Tipsl8rkSsVI7lirRBww9xLjeBp1tUC7JTKhqY2Otncniofol4LgbqODlp7sp5Wi6D1rzU2UL34LrdyEah+hhjhNau0jL4qgteXInB9FL2CvB5UXJhHJxpa0a8UCHQKMbUHVvtUIJW6JKVFecR+wjU6dT5Azf4S6F1WQCdOtlB710+bvxWgNBAIkADH5KrsWUjHapsDrNWfqnAP4A/HTCulKD+H8vNyqK6xK45wprUR7QkR6FdemHUo7LSTWMJXKuFoJXKlIiaX3VRIDjTPG4/mchrJxenaOzqpQ7yw6UeAR3S9Vl9Pz9fW0ou1dXRIjLaSW5EQJ6aXldE3evWiy8Wrbx4C17NgkyvuONvw9KTUmVsFp/xcrPCpBCtFibg/89UZ1+gUOQDV+k4gTzw8tbWLXm0asex7JSgQBAIkwIJAa1v6qba4Bwiwk6RRVTnivoPvp9txL51GXkvMVL52ZdtmXqzP3HadpiqyXVqJXGxQ0Ksr2k2FNItcTjG52Gasrpgepzc/6G5rbqQ1IiD4K00ijpeI6SU3lrfOKyqmmUVilJeI4zUkj9duzL4Nqytmn88zvcVYXTHTPZx97cPqitnn80xuMVZXzGTvZmfbsLpi9vg9iNUVs4emvy3VTuTyt3lqpfForctnTDamJ8rVHGVsMKyu2MVw9XMryW5VRL/2c01OqyuaA8+bvStXB5y78CZbp9sFnt/V2kyvCMGrYc0L9NJpp9E/+vaL5D2zoJBmiemMM8VKjUObmxGTy4bqKy8+TWeMGUdVw0cpnWiJxEz54xOP0ozL5lBFZR+lOsyJlj/6IF19DWJyxQLn1n8qDkBMLhVKiadBTK7E2ZlzIiaXPxx1KiWR3xed7E+FLYjJFTx1xOQKnrkfNSImlx8UM6+MbI3JlS6ehMgV9hTH3pLbT+9cGLXCI4/s+vP6bcbhs0YPp6eWLomkReB5PQLPJyJySSdy4Pmqc8+nrZWVQvSqp41itJd5+1JLC33rjfU0UQgm5xYW2Z7bCDzvfMlL5CEEIhfRsaOHad3a1XTlvK5rlDNt5xQQuZwZ2aVQXcwhsdLj54LI5Q9ViFz+cNSplER+X3SyPxW2QOQKnjpEruCZ+1EjRC4/KGZeGRC59PYpRC6P/oHIlRkil3l1xZqOdnpVBKx/WQhe68VKjUVitNcNG96gB2bMpBIRqP7iohLjdYkY6dUvN7RWI0Qu5xMpkYcQiFwQuex61s733qa62hqafNF0x47n93RFiFzRyN34wtFZASWAyBUQ6ACrSeT3JUDztKwKIlfwboHIFTxzP2qEyOUHxcwrAyKX3j7VQuSSKxje98DyuLTeX7dMS5qIyaWlW3wxqlGszLhOCF1rmxpobeNJOtoevVLKGDGtcboQu74qRK/zCot9qTOVhSAmVyrpo+5kEPBb5EqGjSgTBNwQQEwuN7SQVncCiMmlu4dgn1sCiMnlllj6pkdMLn19p4XIpS8eNcsgcqlxyoRUO1ua6DUhdr0mRK+3m5uimsRxvIaIkV1ThOA1WQheLICl2waRK908BnudCEDkciKE4+lGACJXunkM9sYjAJEL/SPTCEDkyjSPxm4PRC59fa2VyMUjumTAdzMya/B33XBC5NLNI8HYw9Ma1wnB63UxtfEvYrTXcfHdvFXm5ArBq5imFJYY76fnFwRjmIdaIHJ5gIesWhKAyKWlW2CUBwIQuTzAQ1btCEDk0s4lBN1/dQAAIABJREFUMMgjAYhcHgGmUXaIXPo6Ky1Erg1bd9DNd/+adJyuiJhcmReTy3q61tfVKq2u+IUY2bXv7PH0ZvNJ2iKC158QUx3N252vv0r7pk2nc/sNpAuF6DUyL1r0OnRwP215cx3NvnKB4xVj5eOP0MzLr6Lyisq4aVXjB8nA/dctWkxlxfl0rLbZ0Qa3gcsTiZmCmFyIyWXXEd3EgfJb5FI9pxxPoAQSIPB8AtBssiAmlz8cdSolkd8XnexPhS2IyRU8dcTkCp65HzUiJpcfFDOvDMTk0tunaSFy3f/QCnrx1U20ftVvtaMJkQsiF3dKa+B5lrd2CtHrzZZGerOpkbaK141rX6HHL5hMx0vLjH48wBjpFRrlxaO9isUqehC5ok9xiFwQuSBydRGAyOXPLQBELn846lQKRC733oDI5Z6Z1xwQubwSTE1+iFyp4a57rRC59PZQykUuOUrLCZPdNEanPEEch8gFkctO5LL2PZ7I+OTy/6bar8ygN0WsrrfESK+GzuiRXuNrauhrH7xPp1x+hSF8DcnLj9mFMZLL+exe/uiDdPU1i6iouMQ5sSUF/3B99/a7XeczZ2gUq3I+vWIZXbvoNk/lHBPi57q1q+nKeTd4Ksea2e1IPJXKDx86SJv+upa+cfVCleQJp8FIrokJs/M7oxtf+F13ouVB5EqUnL75IHK59w1ELvfMvOaAyOWVYGryQ+RKDXfda4XIpbeHUi5ymfHEismlN0IixOTS3UP62cdrNL4nRnptFLG83hSC1zbxauzsjDJ0SJ4IYh+O58UjvgaJoPbJ3hCTK9mEUX7QBPyerhi0/agPBKwEEJMLfSKTCCAmVyZ5E21hAojJlT39ADG59PW1ViKXvpjiWwaRK109p4/dLULgeles3Pgmi15iauPb4nOTRfQaJkZ2yemNF4rVG/snQfSCyKVPn4Al/hCAyOUPR5SiDwGIXPr4ApZ4JwCRyztDlKAXAYhcevkjmdZA5EomXW9lQ+Tyxs/IDZHLB4goohuBjSKA/aamkPDFgezNW1GPHjS+sIhGCeHrXCF4nS9eVXGmN6rihcilSgrp0oUARK508RTsVCUAkUuVFNKlAwGIXERNbU3UKf6R+OOm8Y/fe4TeO0RoC7lPJAgdE/+M/TK9cLQ5TdQxc5mcJ1xupJ5wmbKeKBtCVina0GU7WxPfBm6qqa0ONnSQqa3GH4BjcDC4hcp2tCFErIu3AocuxiHecgvVx1voPTdHvMR/TS08b0Pa05VHppdlRL4L35jTy/JclRG2wVoH+zVkoXyXFlvrVLBTtt1aZqzv3drVvQ7Z1phsYrSrW/qoi15027rxViyzO7tQud/60rV059Sb0uEym5U2aiVyfbb/EM1a8OOYjsDqiuOoavioCB+7mDp2+9a88CcaM3YCDRk6Ioqt2/2rn1tJ4yZMosFDhkWV49d+LnTLxnVUXFJCY8dP6tYP9ny0i/bu2U3TZ87udkyuDjh3of3FZtnSB2jBjbdSfkH0ioZc0EvPrKDxk6bQoMFVtn1PdXXFtrY2Ou+Ci2L2X6c4WvFWV+SpjG8LoWujeG0Sotfk1S9EBbHnSnk64wVhwesCEdNrVH4Bqa4Eh9UVo92GmFyJ/R4iJldi3FRzIfC8Kqn46RCTyx+OOpWCmFxq3uCHvPaOdmrvbKd9n/2D/v7eO/Tlyy4zBImO8H7+3G4IKPwu0or9xmf+Jz4baXm/8R7+HN7PaYy0prJC5YXzyf1CtAiVJdKLz6VFuVTd0GSUF11HKJ+Rn9Oz/fKztI/LsNYhUnaVZSojbIe0T/LostFSn9HukE1mHiy6GAxMPAzbwu2LMLDwuLPzTvq9+Fcr/mFLHwJ5lEc/Ef/uE/+wgYAk8K/0r7RkyRIA0ZSAViLXtCvuoMtnTKbJ546he37x+8hqinNvvpfmzJpK8+dM1w4jAs8j8Dx3SuvqinYd1YvIZS1vxfJHqP8lX6ftRUW0ORzXy5qmt1i98YqPd9OY/EIaOXEynSkC3heLEWB2G0QuiFx+XFwhcvlBMXYZELn84QuRyx+OQZRysrWBWjtaqU28WjvaQu/t4ntn+D28v+bQUTq06zMaMuX0aEHFJHZERAopuJhFiojoIgUZszASEoDMAosUbULv4bRh0UaKHVECkEUUiYg+UQKORVTicsPCSVQdLLuYhKRQu6KFJOO4+NfYejKmm0bRKLpA/Htc/MMWDIEf0g/jilxFeUXCkB7UQ9yrmf/l9MiJ7BMfQsfC75FjYnCJ3CfTmPNFyjTyifLM9ZjK5DxRNvAxUbaSDSar2QYjj209vN+opautlvZYbcix2Gzm0K2eSNld3KLbz/WGj4mhX9Zjsq0RG0Sa8rcLqO680MgseSfN9hvfw/fW8ntBfq4YzdVDjOQKLTBlPe70PVRHuGxZm6UOpzJiHe+y3mRXeCCa1zJj19m9PbHYsV+j2q/43TaPHGCnWEYsm2S7QleJLp9/8MwWiFzBXDoTqkUrkUsGnq8a3J8Wfu/nEZGLV2A0i14JtTRJmSByQeTirhW0yGUnmG0Tgew3iSmOW8X7WyKuV524yb34w13Cuk5aN/oM4ww4UwheZwuxa5yY6jiuoIjOEZ95g8gVfYHASK7ELpgQuRLjppoLIpcqqfjpMlnkmjj1y9QmRKBWIQKF3oUwxO8RoYgFIxuxSKQzxCNDNAq/G2JSeD+Xx/nMZco6jPQiHb+bBSi7MhXr5vLcbCNoBE0V/x4T/7DFJsAPcbniD2C5PXIFsZE0qeM8WpX/jPE9R+znh3vxf9dnsZ/TsyjCafgzp8nh9OGyjM+cj8sV/1g8kHUYn41XOF84jVFGuKw8cay0uICamjui6+DSTPUZn7lOkd5cB4sfXFZ0HWxjyI6IfTKfuR3cLnMdsi1xeMj2y/oidTjw4HxP/e9/0+yrrqWysp7opmlEAKsrppGzAjQVqysGCDuBqrQUuaZOOptY8JLTE1nkuvnuX0e+J9DOpGZBTK6k4kXhCRLYKYLXs+D1tni919pEe1tbbUsaK4Svc8WIsLPE+0QxVfT8ijI6VtucYK3IBgJ6EUBMLr38kW3W1LXUing7jUbMnaZ28d4qVtIV31s6mumkGGUTOiZe7eK4cUzsE585jTUf728Wr9aORqpvOSlG6XTlq2+ty0i0RXnFlJ+TT3nilZ8r3nvkhd/D33l/jphMZKTJixI0QuJMSBSxEzJY/JD7DfmkmzhjEU7Cwo5ZUIkILvGEE0OcCQs9UkiyK4vt53/dxJmwqGQjJHUJLuF8UeJMSLgqFgx13hCTS2fvwLZECCDwfCLU0jMPAs/r6zetRC6eljjxnNH0o1vmkfnz/Q+toBdf3RQZ2aUbTohcunkE9tgRqBcju95lwUuIX++1NBvv+0QcMetWIob1niVGeY0rFCO++D2/iEbk5wMqCKQlAYhcaem2pBjNU8QMEUmKQ0JI4ulchsBkEppYjGpkUUoca2RhyiRQRdJyXotAZYhYnNY4dtIY3RTkxoKOIQQZwo98DwlAIYGoSyjKE9+70oQEIuM77xfpuo6H98syjeOiTFO6UH3heuwEKQV7pN1dtucFiQ51pZAARK4UwkfVSSEAkSspWLUsFCKXlm4xjNJK5LJi4tFccnv5iV9S1eABWpKEyKWlW2CUAoFqEUPkHSF8/U0IXiyA/a21hQ63d38wK8/JEdMceXojT3MMiV+D8/AQooAYSVJMACJXih0Qp3qr6GQVoOToJhXRKZZYxaJT6Fhj4KITj/LhkUihV5EYUVMS+VyUG9onj/NoG/5svOd27TfSiLTF+Zw39PmUynJqac4VwhJ/F+Uax0IjnrCBQLoRgMiVbh6DvU4EIHI5Ecqc4xC59PWl1iKXvti6LENMLsTk4t6gQ0wuu/MlkdUVa4R2te5YbWS019/EqK8aIYZZtxu2bqGG006nqmEjaawQvc4Vcb4qxdSIWFsiq1/98YlHacZlc6iiso/ry8HyRx+kq69ZREXFJa7zIiaXa2RGBsTkSoybaq5UxOTiwN885Y6nw9W31Ic+t9RF3hs/q6P2xlY6dsqJyPQ6Q6wyT9ELj24KjXAKTcNraK1Xbbav6UryS2lx2830YvFqai/ooKL8aFGpWIpKLDYZx6LFKENUYrHKki8iRoXz9Soo99Vuc2H9K4roeF2ziH/VtYT9Bzu2U/Xxo3ThtBlJq1fnghP5fdG5PUHYtu+zvbTj3W00a/bVQVQXs45sErmeXPYQYnKltLclVjliciXGLdNzISaX3h7WSuSSgec5Jle6bBC5IHJlmshVVpzfLSYXx/MyRnuJ2F7vNTfTjtZmmrNpI20bOox2DzwlcroOzs0LB7UvDI36EsJXT2OlHqJEHkIgchEdO3qY1q1dTVfOu8HXy+IrLz5NZ4wZR1XDR/lWLkQu31DaFqQqcvHKaieaawwxKiROdYlSLFI1sFgl9psFKxador6HBS0ebRVv49XZKsS/NeKf282vkU4qI6TMI50yOfA8RK65brth1qaHyBW86yFyBc/cjxohcvlBMfPKgMilt08hcnn0D0QuiFzZIHLZnSZPv/BHahz1JfpgwADaIUZ7vSNEMLvttPwCGiOC2o87doz6frCDrvjmPOWzDiIXRC67zrLzvbeprraGJl803bEv+T1dUXV0pKNhcRKwQFXTVC2Eqmqq4Zf4XCP2NXx4nBpzm+hwxdGQgBUWseqMz7VCpKoTI6TqjLhQfm08Qqk0vyf1LOhJZQW9Qu/ie5nxvScNqOlLhW0FlDMyPMVOg5FOTm2HyOVEKP2OJ/JHlPRrpb8WQ+Tyl6dKaRC5VCjplwYil34+0cEiiFw6eCG2DVqJXBxsfs6sqTR/jvODi05YEZNLJ2/AFi8ECvNzyG4kl2qZMqA9B7Vn4et98bLbThXxvM4UwtcZIr4XT3U8PS9fBLcvUK0G6UBAmYDfIpdqxc3tYppv03FDpDrRJEQrQ6zi76HPJ1i4YiGrhQWtULoaka5WHHcaPeVkA4+QkiJUTyFMGZ/DwpT8brznlxnHzPtKxT75vbywwljpDZteBOymK+plIawBAXUC2TRdUZ0KUqYzAcTkSmfvubMdMbnc8QoytVYi12f7D9HC7/1c21UUYzkGIleQXRZ1JZOAV5HLaltTZ6cQu8KrOYrA9rtEYPsPxFTHWNsEIXh9SQhfZwjBi0Wws8QKj2V4yE6myzO+bK8iV3XTMUN8MoSpyKiqkEhVbewLiVRy9JXxWQhXXkZTsdDEAlNFUW+qLKyk8qLQ5/KCCupVWB4lSpUVsFDVNcKKBSqOO4UtcwlA5Mpc32ZjyyByZaPXM7vNELky27/m1kHk0tfXWolc5tUU7ZC9v26ZliQhcmnpFhiVAAG/Ra5YJuwRMb4+aWuh3ULw2i0+fyzePxbv9SKWkHkrEas68mqOedSDhonRXyPyCsWIr3wamVdAw8U7NhBwIsAiV31LA+2rORIRqaKEKSFUVZtEKjlFkEdfcbpO8S+RjVe6M4QpQ6wSQpUQq1iwCn0Oi1b8Ht7H+/kzp8vLwcqliTDPljwQubLF09nRTohc2eHnbGolRK7s8TZELn19rZ3ItfRXd5E18Pwfnn2NHn7seS1HeCEmF2Jy8emdSasrqk5XdBu4XCVmyoH2NvpEjPbaLV4fidfAF1fRc5Mm0+4S+xUSOd7XSDHVcaTxXmAIYKPFCLBeQhzD6orxf3jc+k/lZyyIwPP76j4jjsl14kQ15Z9WEpryZ8SskrGr5Oiq0L4jJw+pmB4zDY+MkiIVi1BnnhxtrLjXPiTH2G8IU/zOLzHSioUt3lcqRlj5vakGnve73njluYmPFqRd8epCTC5dPOGfHSq/L/7VlhklISZX8H5ETK7gmftRI2Jy+UEx88pATC69fZoWIteGrTvo5rt/TTqO5ILIBZELIpfa6nyJPITIwPMF5ZX0UVto9JcxCkyM/PpEfN8jhLBGMSXSulXm5NItL79A//j6HBpW1lOIYIWGGDZMvPJ79HC8KvMP13dvv9sxXbwEjScb6OkVy+jaRbd5KifTV1fk6YBHG4/Q0ZPiZbwfNt6PifcjjeHP4vuRhsPGSoG8uV3RrzC3MEqM4ml/hhhliFTWUVWh0VeV4dFX1phUQQSej9VhIHJ5OpUimSFy+cNRp1IS+X3Ryf5U2AKRK3jqELmCZ+5HjRC5/KCYeWVA5NLbp2khct3/0Ap68dVNGMk1ZhxVDe8SFOxGYtjtW/PCn2jM2Ak0ZOiIqN7odv/q51bSuAmTaPCQYVHl+LWfC92ycR0Vi1E7Y8dP6nbm7PloF+3ds5umz5zd7VjtiRp6+fmVNHfhTbZn3LKlD9CCG2+l/ILuwc1femYFjZ80hQYNrrLNW19XSy+sepLmX7845tmMkVzOF7pEHkKcVldkeetgePTXJ+EpkKGpkM00/8Xn6cGvTKeTIq6X3HLFhyox2mtEeNQXj/4yRoHl5tMAMR1SbhC5nP1pl0KO5Jp46cWGUMUC1qGGg8bnI0K0OiaFrLB4dbjhC9cVDSgZSJNzplAf8e/wKcdD0/3klD85qsqY9hcSr6oq+1NRXjHVN7a5rssuA0SuaCoYyeVLt3JViN10xQ92bKfq40fpwmkzXJWVKYkT+X3JlLYn2g6IXImSSzwfRK7E2aUyJ0SuVNLXt26IXPr6hi1LucglR2k5YbKbxuiUJ6jjiMkVFGnUk2wCQcXkSnY7uPwWMcKLR34Z4pcY8bUnPPKL953oiI79xek5wP3IcLwvfmchzBDAxKtIYfRXEG1KdR0cTP1A3ed0oGEfHazfTwfEy/zOn1nYcrPxdMC+Jf2pb3G/0Is/R76H95eE9rOY5XbzGnjebX1IDwLJJoCYXMkmjPKDJICYXEHSRl1BEEBMriAo61EHYnLp4Qc7K1IucpmN4sDzOotZsdwIkUvfDg7L3BHIJJErXsuPGqO/QtMfzSPAPhXfreN9eHLjKbl54bhf4fhfYRHsVDH90XnyozsfpCp1Y1sj7RfxrlioCglXISHrYMOB0Pe6fUaMK6eNg6b3Lu4rBCuTcFVq+sz7wyJW/9KBxAHak7lB5EomXZSdCgIQuVJBHXUmiwBErmSRRbmpIgCRK1Xkg68XIlfwzFVr1ErkUjVat3QQuXTzCOxJlEC2iFyx+LSLA5+KUV8c78s6CuyIEMasG4/wGm6M+BIjv0xB8DkgPo8M02U72dpA++s/j4hVUryKCFlCxOKVBJ02Hkk1sGwQDSwdRANKT4l8PoX3lZxi7OtXOkAIf/pIfxC5nLyK4+lGACJXunkM9sYjAJEL/SPTCEDkyjSPxm4PRC59fQ2Ry6NvEHgegee5CyEml/OJlEjMFKeYXPFq9Xt1xQYxxXGXMe2xhfYKEexj8flT8b6zpdnWjKq2Nlr4+qvUePW1lEudYjRYPp2an0dDxHuVGAGmujkFnq9rqe0adWWaQmgekcVprNs1dA1tE/92i38FuQWGQGUIV0LAMkQr8R4Rs4x9g0Vsq6K4ZgexuiIb4CYOlN8iF2JyRXcBN75Q7fPJTofA88kmHHz5ify+BG+lXjUiJlfw/kBMruCZ+1EjYnL5QTHzykBMLr19qp3INe2KO+jo8RO21LC6IgLPI/B86NRY+fgjNPPyq6i8In6MItUHchm4/7pFi6msOJ+O1doLN+YT026Rg3iXu0QeQnQSuWK1jYPf7xeCFge7l/G/eCTYgfpa+uZrf6Z/n3lZt6wcAH+QCHTPghdPeRzCn8X7qWEBbKD4zml4CuHfP/0b7di4lXInlESmDn4hphDuF7GxWMhqaK2P+yvDwddPYcFKClfincWs3A9aacjpI+mM08YZ0wf92CBy+UExdhlYXdEfvhC5/OGoUymJ/L7oZH8qbIHIFTx1iFzBM/ejRohcflDMvDIgcuntU61Errk330uDBvah39x7u97UTNZhJBdGcnF3wEgu51M2kYeQdBC5YrW88WQD/WnF/9BkseLn50L0+lxMd2Qx7DMxEiz0vZ2OtgiB6uRBoWYdEO/82h/+vJ9yxP6O5mM0UPybI/49LP6ZNx5VJUdesYDFn08piRazeCQWj9Ky29yKlM4eJoLIpUIp8TQQuRJnZ84JkcsfjjqVksjvi072p8IWiFzBU4fIFTxzP2qEyOUHxcwrAyKX3j7VSuRC4Hm9Owusy3wC2R6Ty28P82qE++o+pc9rxUsEdefA7p/V/kN8Fy/x+ejJwzGrzCkZTB1CtCJ+FYtXkRhtJYK2U1E/osL+VCRGaA0QAfEH8kuM/BqQm0uniKDv8rvxLl4FWb4ypN/TFf3uIygPBNwSQEwut8SQXmcCiMmls3dgWyIEEJMrEWrpmQcxufT1G0QuH3yDwPM+QEQRWhCAyOXeDZ0i3hZPHfykejd9UvNR5H2P+MwrEvJx68YrEHKMqyE9h9KpvYaK9yoaIt5PFd+reg0zphbm9silxs5OY9TXAfHaL0aCHRSjv/jzQfH5gBgVxu8NnR1xja7IEeJXWAiTwtdAsW+gEMVYHON9fcUrUzeIXJnq2extF0Su7PV9JrYcIlcmejW72wSRK3v8D5FLX19rJXLxdMU5s6bS/DnT9SVmYxlErrRyF4yNQwAiV2w4DWJq4a7j79Oemo9p74lP6OPju+hTMSJr55F3bTP1Ke5Hp/UebQhXLGCxkBUStEJCll9bnRC5DgrBS4pfhgAmxDDed7BdCGLivdZBCOPRXjwSjAUvHh0WEcXkyDDeJwSxwjQcFQaRy6+ehnJ0IQCRSxdPwA4/CEDk8oMiytCJAEQunbyRXFsgciWXr5fStRK5NmzdQff84ve0ftVvvbQp0LyIyYWYXNzhEJPL+bRLJGZK0DG52jvb6dOaPfTaE89Q+4UFxsgsHpH1SfVHdOTkoW6N5LhYw8tH0cjK02hExenG+0jxzuJWbmsOPb1iGV276DZnOHFSOK2uqFJ4o1gZcp8xEqxLDOvx+p9p74iRtLP/QEMMOyHSOG08KiwyGizPIooJkazw6FH6YONfaM7VC52K8nTczYp+fotcqos5eGpgjMyIyeUPVcTk8oejTqUk8vuik/2psAUxuYKnjphcwTP3o0bE5PKDYuaVgZhcevtUK5GLY3LF27C6IlZXxOqKoTMEqys6X1iXP/ogXX3NIioqLumWmGNhWacWspD1We1eau1opX8N/+OMPcQ/nlo4slKIWBVCxBLvIyqEsCXErFN7VRnH7TYOPK+LyGVnnzXwPE+N3B+eChkRw3g0WEdoVNgBIZDViM/xtlOrq+nSne/R8xd/TUyBzKW+QhQzXkIQC33mqZHiXYwM6yfe+4l9+QmMDoPINdH5BAgohRtfBGSSYzUQuRwRpV0CiFzuXQaRyz0zrzkgcnklmJr8ELlSw133WiFy6e0hrUQuvVHZW4eRXBjJxT0DI7mcz95EHkK8jeT6HY2bdSF9dvJTIWjtNqYZ7gnHzTrRXNPN4LL8npGRWGN3f4kGzBwRFrVGE4/Ycrulm8il0r4mIYTxtEgWvMxxweT3ziOH6cK/baffX3SxSnFGmp49ckwiWEgAM14sioWFsn5CJOsjvvNIMt7cCCsYyaXsioQSuvFFQhUkIRNEriRATXGRify+pNjklFcPkSt4F0DkCp65HzVC5PKDYuaVAZFLb59C5PLBP4jJ5QNEFKEFgXSLycVB3Q/W748K9s4jsljU4mDwHZZYVBzMvap8mDG10BiNJUdniRFaA0pP0cIHmWDEfjHy66gY9VUtxLBDIj7Y0Y42OiK+H2vjz+Jl7Gunw+K4240D5vcRIhiLXuU5OYbwVSn28edyIZjxd34Zx8T+IaWFVCxGi9U3uq/LrW1IDwJBEEBMriAoo46gCCAmV1CkUU9QBBCTKyjSqa8HMblS74NYFkDkspCZdsUdNLBfb3pq6ZLIkR8seZD+vH6b8f2s0cOjjvE+iFz6dnBY5o6AriJXY+tJ2l29KzIiS65kyPGy+Jh1613cNzKlUMbJ4vdh5SMpPyffHRSkTiqBail6hYWvo0L4YgHsmHg/Et53THw/IgSyeocA+rEM5YD5FUIAKxeiV0QAYzEsVwhhhjCWJ46FPhtphBhaKb73EmlCY8ewgYA+BCBy6eMLWOKdAEQu7wxRgl4EIHLp5Y9kWgORK5l0vZWtncjFKyzu/HCv0aqlv7qLpk46mzhW1yXTJtJv7r3dW2sdcrPAxZtZ5PrDs6/Rw489HwmGz/ZNPGc0/eiWeZHSIHIl1S0oPEACqRa5DohRWR8LMeuj4x/SR+J9jxiV9VHNh3S44QtbCqdVfsmYYsjvwytGCmHrNBrd50zqVVAeIDVUFSSBfWKU2BExMuy4EMBOCNGrht+FCFYjAucbL44dJvZzIH3+fiSB0WLm9vB0SkMAC48Qq+ARYmJfpRhN1ssYSSaORwSzsIgm0peJNNhAIBkEIHIlgyrKTBUBiFypIo96k0UAIleyyOpXLkQu/XwiLdJK5GIBadDAPoaYxYLTz37yHUPksgpNycDJdc+ZNZX2HTxC2979MDJayypqWW1BTC7E5OL+iJhczmelOWYKi1cfiymFu4//XYhZH4rphh+Kz7uoobU+qqDb6XZ6tWgt9e0zwBiZNUwIWaN7nykELfG5fETcSuMFnneylufZf/f2u52SxT2eiTG5nIAcPnSQNv11LX1Ds9UV64UQtu9kS0gMaxcCmCGOhcUw8Z1Hk0mhrEswa6daIZJ1Whp98Ye7xJ5OWjf6jLg4eLxghRDCjCmUQvRiIazSEMr4c0g0k++GgBYZSZZHeTFKxuqKTj1Q7ThicqlxSqdUiMnl3luIyeWemdcciMnllWBq8iMmV2q4614rYnLp7SGtRC4esfXyE7+kqsEDokSuDVt30M13/5qStbqiWci6/6EVUSIXi22Lr5tN8+dMNzxptQUiF0QuiFyjYl7lmtqaaOeRd+m9w+/Qx3s/oNJDhfRQ60O26Yu5+mYwAAAgAElEQVTyiulLQsA6s99YOr3yDDq73zn06V/+TjO/fgVVVPZxfSWFyBUfmXV1RdeAbTLoKnJ1iID5icbkYgGMhbGaTiGECbHswDtbqVGIX/Vjx1O1sY9HjIVGkoVEtNAIMw7Sn+hWakyhDE2drDRNszz1nbcor6wnlYwZS2XieE8hjpWJqZhl4r2nGA5RRrwvh0rEK6gNgeeDIt1Vj91Irg92bKfq40fpwmkzgjdIgxohcrl3AkQu98y85oDI5ZVgavJD5EoNd91rhcilt4e0ErlYUFr+u3/uJnIlcyQXx9viTU6FtIpcLLz99M6F3UQuKcaxyPUv/9IVv4vL+vzzz2nt2lfpxhsX0fvvv0+7dv2drrzyKqOe//iP39J1111PFRUVxvcVK/5AEyacS6effnqkp/i978knn6BJk86nUaOixQi3+x9/fDlNmXIhjRgRPYLGr/0MgLmVlJSKeqZ0O3OsLM0Jqqurie343ve+b3vG/eIX/x/98Id3UUFBQbfjjz32v/TlL0+jYcOG2eY9ceIELVv2P3THHXfGPJvfeOMNahOrzn3lK1+NmebBB39H11yzgHr37m2bxtxvnC4bTmXJ/OvXrxcfO2natIvjFhnh9/3vUw/OofB8btdPtx3YRlv3b6XtX7xDbx98m9794t1IvSNoBE0V/x4T/4aWD6WzB5xN4weOF68JQtg6k07v03UOyEz/9V//Sd/61lzq27evE5Jux//93++nW2+9TfSnEtd5/+3fup/XdoU0izhRhWLlP7utoaGBHn74Ibrrrv/Hdf3mDIcOHaJnn32Gbr55sadyrJnt/Oe1gv3799OaNS/Tt7/9Ha9Fxc2/ZcsWqqmpppkzL3WsR2hAxqbSpx0LEwlUz6lmUeFx0T+OixFj/F4txC+eYlnNn8V0S37n76FjoemXvI9FsvYYhly6cwfVFBfT5pGxhWXOyj2yZzimWC8hlPUSn1kQ43f5neONsVBmfJfvlnQcuN9JLnPjCxW+QaSx/g4HUaefdXCftvbnt956i44ePUKzZl3mZ1VpU9bevXtpw4Y3aOHC69LG5lQb+sknn9DmzZtowYJrU20K2fXplBulYEC8ewC77A888BtatOjb1KtXL4XSkUQXAixy8XPEPff8VM0kcY1WvZdWKxCpdCTAzwpLlkRrADrama02aSVyscD04qubjPhXcrpi1eD+NGvBj+mGuZdGxcHyy2Fcz9HjJ7oV17d3ecSOeCO5OGNdY2vCowT8agfKAQE/CBTm5VBpcT4dr2t2LG7Xsffp3cNv087D79G7h7bRO4fess3DsbLG9p9AY/sJQWvARGOEVmlBmWP56ZCAHzSliJIO9majjT1Ff+aRXA1N6bO6Io8MM2KKGVMrQwJZKMaYmEIpRonVi8914tUgPvM7v4x94nuDePcyiszaR4pEBy/lUWNiRFlo9FjoxaPN+J33lYo0ReJ7kZDEioQqViw+84s/G/vkZ2N/KC0La9gSI9CvvIiq65uprV3hrxGJVYFcIBAYAY7J1a+imA5VNwZWp18V4R7AL5KZVU5JQS7l5+fSiYaWzGoYWtONwCm9i0FFUwJaiVzMSE4HNPMyj6RKNkfrSC6nmFxsDwLPJ9srKD8oArECz++t+YT+JqYcvnvkbXrv0Nu0Q0xBPNna0M2sQWWn0jkDzqVx/fk1QXyeSD0L8BfLoPyHeroT6FUSErkSna6Yjkw7hNE8ddIqhpmFsDrBpE4IaPUd4l2kbeDvRp7we1hMizWqzA8uPOIsJIAJ4Uv83TskjLEIJkUzcVyIYZzGEMfEsWJDNAun4WMiH38vNtKF9/M+03fjmNiXnyGKNALP+9H7UIYuBBB4XhdPwA6/CCDwvF8k9S8Hgef19ZF2IleqUVlFLqyumGqPoP4gCbDIVd36Ba3fs9kQs3ik1o7D26mmubqbGX1L+htClhS0zh04iSqL3MfOCrJ9qCv7CGSjyOWnl0/KkWLEI8aEWMhCWEQQEyPIhBjWKL6f5Dhk4jiPImvifTyiTBjCn419fEyU0SQGH/ExHoWWio3jlUWEMbOoZoxCY6EtJLAVi9ElhgBnvMJimyGwhUaqhUS08Eg1+TmgkWoQuVLRc1BnsghA5EoWWZSbKgIQuVJFPvh6IXIFz1y1Rq1ELo6P9ef127oFmOe4WJdMmxiJm6XauETSWUUuLkPaxZ/PGj08svIif0fgeQSe536QrqsrVjcdo3e+eIt2fP4O5bzfRv9T8D90+OShbqdOeWGFmGY43hilxVMPzxEjtXau30pnjBlHVcPjxweShSUSGPiPTzxKMy6bg8Dza1fTlfNuSOSSFjMPAs8nhvOdrRuNFRfPnXRhYgV4yJVpqyvyKLFGQxwLiV+GOCZFMUMsM4lmhlAWEtSahAcM8cwsqpnSyrJCApso18jXQXLC6p1rX6FlU6ZSjYj/mMwtWSPVTulZRC2NbZQv2iZHsf1jx7tUW30MgefnzE2mSzOqbASeD96dCDwfPHM/akTgeT8oZl4ZCDyvt0+1ErmsKxlKdMkMPO/VPRC5IHKli8hV31pnxM56V4zQMqYeilFa++s+N06B3uLfteLff4h/JfkldJaIm9U1SutcGlHRXchyK5JA5ErsanPs6GFaB5ErCp6bFf38HskFkSu6H7vxRWJngH+5pKj23PJH6Lx/uopyxEqViYpqLMAZgpsU0ox3e1HNvxbELmnS3j3Ut66O1o47h/LEaLR8kTTqXezL4ymb4pUnj8l9xnunMZ0zlCf8buyn0P7wdM+uMsX+cLro+mzyxrPHXEeUPe6oJfL74q6GzEsNkSt4n0LkCp65HzVC5PKDYuaVAZFLb59qJXLxiK2lv7qLpk46O4qajNP1/rpl2tGEyAWRS1eRa6eIm7Vl83o6ULePXulYQxxXy27jgPDn95lCp346kC6+8lI6v2oCHat1DjwPkSv+5ajxZAM9vWIZXbvoNk/XLYhc3fG5EVYgcnnqfo6Z3fjCsbCAEqx4bCl9fc486tmrPOk18mIAoZFq0aPSWowpnqHRaaGRaV0j1ZrFd+v0T7Oo1i5i9te3yVFvoZFqZ+752BC5Vo8dl/Q2BVkBj1SLL7qx2EY08MhhGvX39+mdaV8jnv6WKwS4XF7hTHDNFWXwi5c64BcfMz4b+3mfeBfMe4Tz8HEO35bb2XVclpkj0hnpI/l4FbXQd6Nc3i9G2MnjRj6u29hnsiHqc8gGI62oM1KXybbIccPWrjq78nF7TDaE03EbeQqu3QaRK8ieHKoLIlfwzP2oESKXHxQzrwyIXHr7VCuRKx1HcrF7EXhe706eLdZtO7iZNh14g946uIm2HniT6lpquzX9zL5nG9MOjVFaYuohTzs0b7ECz2cLQ7Qz8wj4LXJlHiG0KN0IxIvJ1SxEnTbxahViTOidQu9CZGkT+1r5eGQff+7aZxyL5KVI2lAemc6U36jDki4qbfhYN3tknlh1h+rA5j8BlrtYJGPBLfQ5JLxFBDgW/oS/WByT+0NiYShNl0BoEvPMaY0Vh0NCXFcdoTojZZpERq6DRwKWFuZSU3N7JF9IkLQXGQ0R0ii/S0A0BMuoOsKCphQTjXZF29Ylena3zeARQ2SMtk2yY15CUDXbYKlT2ui/V1GibgQQk0s3jyTPHsTkSh5bryVrJXLxtMT7HlhOLz/xS6oaPMBo22f7D9GsBT+mIFdYdAsVIpdbYkjvlUBTW5MQs96krULQ2rz/DXpz/1+7Fdm7uC+df8oUOu+UyTR+wHl0dv/xVJwXf6lbiFxePYP8uhGAyKWbR2CPVwLZFHieFz7oEux6hAU7k+hmEvLahSjRLr7zi5c14FVVO4TQYOwT3zmWXuizPC6+82dxoCOcl/O1G99D+yPHjTyirHD6yHs4XYfYb+Qz1Slt4Lr5eKhs07uwrasOtjVUN5ctkoXt7KozUofJNk4vUou0oXZG2s11if08pRZb+hDgkYu88f8sFrK4aHwWLxbdQvJe1z5OzmKgkUb4mr/LFF15QvlkOUYao9zwvnCe0PGQuMmJI+l5ryF8yn3hGkxpDBvk97DNnLorT7g+c7mGvSa7Iu2UrTS3MzQasqsN0Rwkm8jxCAtLu60cDG7WNF1MzRy4jXLrYVxNJGvzeyhNQV4O5QqjmlpC6xPLwZSyhK73rnaY0xmfw6euudaoOsMHupcZsjLSVyzpIo0w2y/r8lhmN1tMo0gjxyztiuy3jDiN3a5w+0znSpdfZNvD7/L6Z0kbs+xI+y1+MUMLf+Yiq3Lz6fwBWEHeBo8Wu7QSuZiInJpopmM3hVELemEjIHLp5I3MtKWxrZHe/mKzELOEoLVvvRFbq7WD/0bftQ0qO5UuGDSVzh90IV0weCqNqhztGgZELtfIkEFzAhC5NHcQzHNNIJtELtdwkCEuAUPo45dJEIuIcCZRjwVDKbR1CWnRoh8La4Z4F5WvS6AzRD9Rl1nci3wW5Rs2sEgnnuZLiwuo+mRLN2FR2hYqK1polKJfu8gfEQXD6aRtEeHRaHdYBAy3n+vvNETQUNmxxE1DQDQJj04iaISXFCxlfWEW6KIgAAKZQeCH5X3o1yMGZ0ZjMrAV2olc6cYYMbkQk4v7bDJWV9wkBK2N+9fTJjFKa/P+DVGnxvfp+/RGrzfprCHj6HwhaJ078Hyq6jWs2+mjGiS79kQNvfz8Srpu0WIqK85HTK4Hf0Xfvf1uT5cjxOTyhC9uZjdxoPwWuVTPqWS0PtNWV0wGI5Uyg4zJpWKP2zR2ItcHO7ZT9fGjWF0RqysqdyfE5FJG5VvCWDG5OOaeMcAl/M6fxQAjYx+PDjS+h61gQY/3sGBo7DeEQh5bFC7DyBNKI/OF8oT3hQXMyLFwulCe0BbJGy6bRUBZV+h42IZIudLekODIW8i+0PeuutjekLjYbV8kT9gGU7t5lKORz8wh3PaIbREWpvqsHGQa2Z4oFqHybTm0t1PZiuVUu+DGCGXJSvIIgzPecsVc2tzcnMhILukZsw/N+SL7w8O3ImWbKunyjfRRl6+MssIJYqWT3o2kC49usqaX/lcrs8tXtu0JHTb6cqi8rgbZ1WtbRjhP7PZbbLCJQdgtb4RV6EPMsi2O6PJTqM6ZK5+kJUuWhFuJN90IQOTy6BGIXBC5uAt5Fbla2lto4wev04dvv0evl6+jtw9tIZ6SaN6GlY+gKadOoymDv0yNm6vp67PnUnlFZdwerPpADpErGiMHk4TI5f7iePjQQdr017X0jasXus/sIgdErokuaCU3qRtfJNcS9dIhcqmzSpeUWF3Rvacgcrln5jUHAs97JZia/Ag8nxruuteKwPN6e0g7kYuDzx89fsKWGlZXHEdVw0dF2Nitbme3b80Lf6IxYyfQkKEjori63b/6uZU0bsIkGjxkWFQ5fu3nQrdsXEfFJSU0dvykbn1gz0e7aO+e3TR95uxux6RIM3fhTbZ9Z9nSB2jBjbdSfkFBt+MvPbOCxk+aQoMGV9nmra+rpRdWPUnzr18c82x2K3K1dbTRe4ffpo1i6iHH09ompiL2a+1LXxP//q/4x9uQXkOFoDWNJg++iC489WIaVNY1JHbl44/QzMuvgsgV5/q6/NEH6eprFlFRcYnrqzBELtfIjAwQuRLjppoLI7lUScVPB5HLH446lQKRy703IHK5Z+Y1B0QurwRTkx8iV2q4614rRC69PaSVyDX35ntp0MA+9Jt7b9ebmsU6xORKK3elxNj2znb62+HtQtASota+vxorIDa01kfZwjG1eJTWlFP///buBTyK8u77+B8QMIAcxCMIEsWiDWcjFYGGFkFRixSkARFEq4JHtK1oq+/F4/Pa56318dRikSo2BeUMclBABQw2yKEoKiIICqJyNEAAEYEA78zSDdnNZHcO9+zcs/vdXF6XJHPfc8/nP7tsfsz+x/jPCLfMkCvVD3pypVqc/fktoPrjin6vl/kRSCZAT65kQvw8TAJm0/CzGmTJ9t0Hw7Rs1opApQLcXTFzTg7urqhvrbUKuXK6DhHdm8xblZKQS98TPKiVHTt+TD759sPIVVpms/gVW5fI/sP7YpZzTu1zjau0zEDLDLbyxPw4YtAPQq6gK8D+VQsQcqkWZb6gBQi5gq4A+1cpQMilUpO5dBAg5NKhCqlZAyFXapzd7IWQy41a3BhCLgWIIZ/CbCr56bcfn7j7oXG11nIj1Np7qCTmqM6sdXbko4dmsNXJuFrrwvo/0u6oCbm0KwkL8ihAyOURkOHaCRByaVcSFuRBgJDLAx5DtRQg5NKyLL4sipDLF1Ylk2oVcpkfV+zds7MM6N1NycGlYhIaz2du4/mSQ3ukcPPbsmjzm/Ldxj1y5MhhWWR8RR+nZ50hHRt1LmsW/+HcooR9tHZs2yLL3yuUXn0HJj116cmVlEjoyZXYyKp/X3LVxFvQk8urYOLx9ORS40tPLjWOOs1CTy7n1aAnl3MzryPoyeVVMJjx9OQKxl33vdKTS+8KaRVyFa1YLY/86SVZPOM5vdXKrY6QK7NCLvMjiAu/nC8LjWDr/e3Ly86En8pPpVbV2iLNTolcqWVesZVzRuuY8zhZMEXIVfFpP/XVsdL9mt5Sv0FDx68JhFyEXKaA6iu57N6x1PEJa2MAIZcNJBubEHLZQArZJoRczgtGyOXczOsIQi6vgsGMJ+QKxl33vRJy6V0hrUIusydXogd3V+Tuiqm+u+IPpT/Igs/ekK+K1smY6i9K8fc7Y07RnDPbSLfzr5acg5fI2VnnymWXd6n0FCbk2iyrVi6Ta3vn235VJOQS2VW8UwoXzJW+/RO/PtpG/c+GXMnlVOzE9oRcsW6ffPS+7N9XIh27hOcKbEIud+e+zqMIuZxXh5DLuZnXEYRcXgWDGU/IFYy77nsl5NK7QlqFXHpTVb46enKFtXLW6zY/hvj2pjdk3hezZfHXC+WH0pN3/KlVvbZ0bXql8V936ZF9rZh9ttLpQU+udKomx2IKqL6SC1UEghagJ1fQFWD/KgXoyaVSk7l0EKAnlw5VSM0a6MmVGmc3e9Eu5DI/sjh0xFMxx6L7HRcJudycenqN2frdFiPUmiVvbpwjS4zG8eUf553WVK5p3jtyxVbn87rqtXDFqyHkUgzKdIELEHIFXgIWoFiAkEsxKNMFKkDIFSg/O/dBgJDLB1RNpyTk0rQwxrK0Crkmzlwojz87Xua9+oQ0bXziCpmvtuyQngMfkkfvH6RtQ3pCLn1P8EQrW7drTSTUmrdpjqzeuSpm0x8b/bR6XvALucr4L763VjiP1t6qCbnsObFVeAQIucJTK1ZqT4CQy54TW4VDgJArHHVilfYFCLnsW4V9S0IufSuoVciV12e4DBvcq0KYZYZfL4wzPjqmYUN6Gs+Hp/H8jbfcKR/t+kDmG8HWm0awtanki8gzc4jxtdj4OrtxY7kq+xdGuNVLmtQ9v+xZ+93+fTJnxgQZcPOwSp/Jq1YuldLSUnpyJXitc9MzhZ5c9OSyOqWc9IFSHXLRkyu2Ik5qoctbIXpy6VIJdetw8/eLur2HcyZ6cqW+bvTkSr25ij3Sk0uFYvrNQU8uvWuqVchlNp63+mhi9COMNJ6n8bzTxvOHjx6Wom/ekQ1vrJKXao6VrQe3lD0jTz3lVOPjhz+TvL1dpFPHbtLiglaWz1ZCLusXMaeNy938EkLIRchFyHVSgLsrqnlDRcilxlGnWdz8/aLT+oNYCyFX6tUJuVJvrmKPhFwqFNNvDkIuvWuqVcjFlVwnThar8MDL9+bPmSY5rdtLk/MviDkbnX5/7qwp0qZ9B2ncpFnMPKq+b066fEmhZNWqJa3bdajwzNm4YZ3YCbn2H94ni7580/gY4mx558u35Lsj++UPxtf/Gl9ZNWtFemuZH0P8+fk9xGwk/8Zrk6RdhyukUeOmhFy3DpM6WdVl175DSV+5CLkSEx38/oBMn1QgN916d1LLRBtwd8WKOk6uHuJKLk+nX9LBTmqRdLIUbUDIlSLoFO6GkMs5NiGXczOvIwi5vAoGM56QKxh33fdKyKV3hbQKuejJpffJovPqdh7YLm9++caJxvHfFIp5BVf0cU7tc407IV4nPS/sJVc0zpNTqp6i86EEujZ6cgXKz859EFAdcvmwRKZEwJEAPbkccbGx5gL05NK8QCzPsQA9uRyThXYAPbn0LZ1WIZfJxN0V9T1ZdFvZxpLPjf5asyM9tj7YvkKOG1/Rx4X1fxQJtcwrttqffZluS9d2PYRc2paGhbkUIORyCccwbQUIubQtDQtzIUDI5QKNIVoLEHJpXR6liyPkUsqpdDLtQi6lR5eiybi7Yoqgjd18tNNsHH8i2Fq/e23ZjqtIFWl79qWRUMtsHN+8QYvULSqN9kTIlUbF5FAiAoRcnAjpJkDIlW4VzezjIeTK7Pqn49ETcqVjVa2PiZBL31oTcnmsDXdX9P/uimu+/Uhmbpgqn6x9X846eKZMM77MR/Wq1eXyxl0ioVbnM7vKyrcLJX/QHZYVLRjzrAy85S6pXqNGhZ/Tk0tk394SmTd7igymJ1fk/DA/Z3/7PSM8vTrQk8sTX8LBTvpAqQ65uLtibGmc1MK/M8LZzPTkcuYVhq3pyeW8SvTkcm7mdQQ9ubwKBjOenlzBuOu+V3py6V0hLUKuaC+uR+8fJAN6d4sRS/QzHWgJufwJuWpn15fZ66fJa+unyBcl6yOlbml85VTJkX3ZB42PIl4v3bOvkbo16kV+Fg1pCLlinxV2fyEn5Ip1I+Ry9+q6c8c2WfruArm+3yB3E9gc5SRYIeSyiepyMye1cLkL5cMIuZSTBj4hIZfzEhByOTfzOoKQy6tgMOMJuYJx132vhFx6V0iLkCt/6GPS6JyG8sxj91hqPTBylGzdvksmjxmpnSYhl7qQa/fBYpk6v0BWl3wkrx2YHlPr7s2ukavrXiOnH2ggPXr+ssJ5QMhl/dQg5LpVTs2q5fh1g5DLMVlkACGXOze7o5YVLZLadepKq7a5dof4vh0hl+/EFXZg9XHFT1evkj27i6VTXvfUL0iDPRJyOS8CIZdzM68jCLm8CgYznpArGHfd90rIpXeFtAi5croOkTF//q107tDKUivajH5NYYGWmvTk8laWhV/Ol0lr/ylzv5gVM1Hn834m1190g1x3UZ+yK7a87YnRyQToyZVMiJ+HTUD1lVxhO37Wm34C9ORKv5pm8hHRkyuTq5+ex05PrvSsq9VR0ZNL31oTcimoDSGXc8TNezcZwdY4mWL8t/3AtrIJcs5sI/1aDJTeP+onZ9Y62/nEjPAkQMjliY/BGgoQcmlYFJbkSYCQyxMfgzUTIOTSrCAsx7MAIZdnwtBMQMilb6m0CLny+gyXPz58W8IruR7500uyeMZzWkoSctkry6Gjh+T1z2fIpE//KUu3/EuOG1/mo2HWmdKnRb7kXzJYLmnY0t5kbOWLACGXL6xMGqAAIVeA+OzaFwFCLl9YmTQgAUKugODZrW8ChFy+0Wo3MSGXdiUpW5AWIdeToyfJyg8/q7TnVrKeXUHy0pMreU+uXTV3y0Qj2JppNJHfd3ivDDa+3qvynlzQrIX86pJBYvbbOqXqKZEyLl9SKFm1aknrdh0qlHXjhnWyaeN66XZVrwo/oyeX9bOAnlz05Krs9fHN16fLJTltpGl2c2UvofTkUkZpORE9udT40nhejaNOs9CTy3k16Mnl3MzrCHpyeRUMZjw9uYJx132v9OTSu0JahFwmkXk1l/mIv1rL/H7x7r2iaz8uQi7rkGvvoRKZOu1leefoO1K4b2HZs+DihjnS/2i+dPrJlfLji9pWeHYQchVKr74Dk75qTHnlRbnquhukXv0GCbcl5CLkIuSqLseOH5fvDpYmfV7Z2cDuc8rOXE63IeRyKma9PSGXGkedZiHkcl4NQi7nZl5HEHJ5FQxmPCFXMO6675WQS+8KaRNymUzmFV0Fk+fHiPXIy630ros60BJynQy5GjU5X4q+fkcmGk3k538xW/KNMKvI+Np9aon0vqif8XHEQdL6rPYyd9YUadO+gzRu0oyQq5zAjm1bZPl7hFzlT4qpr46V7tf0lvoNGjp+uo8fO0r63UjIRchFyOX4yeNgAHdXdIClaFPurlgRkpDL+clFyOXczOsIQi6vgsGMJ+QKxl33vRJy6V0hrUIuvakqX12m9+TadmCrTFpTIJPXjZev922OQFWrUk3yml4Z+TjiVdnXSY1qNcJa3oxaNz25MqrcGXGw9OTKiDJn1EHSkyujyp32B0tPrrQvccYdID25Mqfk9OTSt9aEXApqk4kh15FjR2T+xjmRJvLvfr3Q+DjQsYhk8wZmn62bIndIPKv2OQp0mSKVAoRcqdRmX6kQIORKhTL7SKUAIVcqtdmX3wKEXH4LM3+qBQi5Ui0e3P4IuYKzT7ZnQq5kQjZ+nkkh1/rda2XCp/+Qaesmyp4fdkV06taoJ9dfdIP0M67auvScig3jbRCyiSYChFyaFIJlKBMg5FJGyUSaCBByaVIIlqFEgJBLCSOTaCRAyKVRMXxeCiGXz8Aepifk8oBnDs2UnlzmFVsbi9bIvMNzZaPxZT66NPl5pM9WzbVVLHtsVdZ7i55c1s3i6clV8clITy6RXcU7pXDBXOnbf4jHV6vY4dxd0R0njedj3ejJ5e488jKKnlwV9ejJ5fyMoieXczOvI+jJ5VUwmPH05ArGXfe90pNL7woRcnmsTzqHXMcbVJWCj8fI5LXjpOTQHhlsfG2o/bl0aXll5KqtRnUaR/SchlmEXIRcdp92hFyEXFbnipNgRfWVXIRchFx2X7/82o6Qi5BLxblFyKVC0dkchFzOvHTZmpBLl0rotQ5CLr3qEb8aQi5DJK/PcCnevbfMZsyffyudO7Qq+/MDI0fJW4tXRv7cskW2TB4zsuxn6RhyvTJltCw/vkxmfDu17Di7N7tGuh/oJl2v6FnhroiEXCKrVi6V0tJSudEI038AACAASURBVOzyLpU+46e88qJcdd0NUq8+IZfdl0VCLkIuQq6TAsuKFkntOnWlVdtcu08h37dzEjj6vhibO5g0boxc27u/nFa3ns0Rem1GyEXIpeKMJORSoehsDkIuZ166bE3IpUsl9FoHIZde9SDkihP4assOeebvU+WZx+6J/GTizIXy+LPjZU1hQdmfXxg3WxbPeC7y5/yhj0lu2xby4J39y2ZKl55c0z6bIKNW/q9s2LMucmwNs86UAT++WQa3vF0an9ZE7zOZ1SkRoCeXEkYm0UhA9ZVcGh0aS8lQAXpyZWjh0/Sw6cmVpoXN4MOiJ1fmFJ+eXPrWmiu5LEKvngMfknmvPiFNG59dIdQyQ7DyoZc5PMwh14HD38n4NS/Jix/+VbYf2BbRuKD+RXL/ZQ9L3xYD9D1zWZkvAoRcvrAyaYAChFwB4rNrXwQIuXxhZdKABAi5AoJnt74JEHL5RqvdxIRc2pWkbEGEXHG1KVqxWoaOeKrsSi7zo4zDBveSAb27RbaM/3lYQ65dB781gq1RMu6TF2XvoZLIsXVs3EXubPeAdGt2tb5nLCvzVYCQy1deJg9AgJArAHR26asAIZevvEyeYgFCrhSDszvfBQi5fCfWZgeEXNqUosJCCLniSMxQ67ruHcs+jpjTdYg8ev+gCiFX9EovsyfXHx75P3Lk6LGymbZ8840sLlwoN950s6xb+6lsWP+Z/OL6X0Z+/vfRz0v+jQOlXr36kT/PmDZF2rRtJxc2v6hsvOrvTZsySS7NvUyyL7hQNpVslKeXPSmvrh4v/Y72k+XG149btJSHrvi9tDm7XWQN5bcvzzN10gSj51RHadYsO0ZN1ffNSQvfWSi1atWWDj+5vMLJGm9ZfoOSkj0yddJEuX3YXZbPtueeflLuvGe41KhRo8LPJ094RTp27iJNm55vOXbf3r0yccJ4GXrniY+0Wj2WvbdEjpQekS4/7VrpNi+NGS19f5UvDRqcbrlN+fMm2UtGsrmi498r+pccN746df5pwimjfsPuultOqVZVDh05mmwJluduokGbN38pyw2nXw0YmHTu6AYvvzhGru/TVxo2PMP2mOiGz//lGbnltqHG+VTL8dgn//RHefDhR5KOO3bsuFQzvKweBw4ckIKxf5e773sg6TyJNti5c4fMnTNbhvz6dk/zxA+ePnWytGnXXpqXe+3xuoNtW7fIgrffkkE33+J1qoTj31+5Qkr2lEi37j2S7qe6UR/zOVB69HjSbe1ssKToXTGmk05dEj+n7MzldJtFC96WuvXqSu5lP3E61LftndTCt0U4nHjM6FHSf8BNRn/EE38Ph+1xavVqcrj0qBgvP2WPVe+vlOJdxdK9R2b+I9XmL7+UZUuLJN+oKw97Aps2bpSV/14m/fJvtDfAp62qGPPWrFFNfjic/H2HT0twPW2p8d7/FDOls/kY/be/yk3G7wan1a1rc0Rmbqbmb2t1dmZPLvP3iN88+LCtSc1zoqrx3+HSk78b2hrIRqESMH9XGDnyZJ/uUC0+AxZLyFWuyGbA1b7VRWX9ucwfJbuSywy5Rjz8aMxfztu2fiNF774j/foPkg2frZXPP/9Mel7bO7KngrGjpc8NA4xfVE68uZ4zc6q0bN3OCKCal61E9fdmvTZF6mU3kIKvx8rrn88y3hgfk5rVasrwWvfLzzteJR1adoo51c3t27a9VM7PvjDm+zOnT5L2uT+RpufHhlyqvm/urOjdRZGQy9xP/CPesvzP9xoh18wZk+XmW4dZPm1fGPWU3HrHvZYh14ypE6TD5Z3kvCbWIdf+fXuN4O9VIzCxDtDMHf57+XtG4/kj0rFTXqUvG+NefkF6/fJXUr+SkKv8eZPstSfZXNHxy403/eZv5D/pWHlDfHPbqN+vb79TTjXebH53sDTZEizP3USDvv7qS1m5Yqn80jj/7T5eKfi7XNOrr5x+ekO7Q8q2e3H0c3LTkNslK8t5yPWXp/+f3Peb3yfd57Hjx6VqFes3uN9/f0AmjBsrtw27L+k8iTYo/nanvDV/jtw46Nee5okfbPU643UH27dvlXcXvS2/uvFmr1MlHP/hqpWyz3jO//Rn3ZPuJ6tmNTHKpOwXKLvPqaQLc7HBvwoXSB2jWXq79pe5GO3PECe18GcFzmeN/3vY+QzBjqhbu3rkNdoM2aOPjz98X3bv3iVdf548+A129f7s3c3fL/6sJDyzbv5yo6z64N/Su09+oIs2/wqtV7uGlHx3ONB1uNm5+R6gSiXvAazm+8eLz0d+N6hzGiFXIm/7saGbqjkfY4Zc5u8Rdw8fYWtwjepVI/9g/P0Pyd9L25qQjbQUMH9XIOTSsjSRRRFy/ac2VgGX+aP4RvPxPbl0v7vios1vyvK3F8n8H+bJ58ZX3Rr1ZFDL22Rou/vk3wsWS07r9tLk/AtiztD5c6ZZfl/VXRQrm8dcxPIlhZJlXHnTul2HCs+ajRvWyaaN66XbVb0q/Gzf3hKZN3uK5A+6w/LZVjDmWRl4y11S3eJKrjdemyTtOlwhjRo3tRz73f59MmfGBBlws3WAZg7i7orJX+S2frPZcFpm3NXM/htq7q7I3RWtziwnd/RT/XHFD1YsMS/kkks7xP7jQPJngPctuLuid0NzBu6uqMZRp1nc/P2i0/qDWAt3V0y9OndXTL25ij1yd0UViuk3B3dX1LumhFxGfcyPJA7JvzrmjonRssWHWmG4u2LpsVKZtWGqjP7gGVm765PIoZxT+1y5re09Mjjndqldo47eZyWrC0yAnlyB0bNjnwRUh1w+LZNpEbAtQE8u21RsGAIBenKFoEgs0ZEAPbkccYV6Y3py6Vu+jA+5zBDr8WfHV6hQj7zcso8tPjBylLy1eGVkm5YtsmXymNjP3+pyd8WDR76XV9a8LC9+9FfZsv/ryHqbN2hhNJO/X/pefKNUr1pd3zORlWkhQMilRRlYhEIBQi6FmEylhQAhlxZlYBGKBAi5FEEyjTYChFzalML3hRBy+U7segcZH3K5lis3MOiQa88Pu4xg63n55+q/S8kPuyMryz33crmr/W+kR/a1UsX44oGAHQFCLjtKbBMmAUKuMFWLtdoRIOSyo8Q2YREg5ApLpVinXQFCLrtS4d+OkEvfGhJyeaxNkD25Zs16RT6usloKtr4kP5QejIRZ99YaLpe2uUKuvPS6siOrrMeW0+/Tk4ueXOWfLm++Pl0uyWkjTbNP3jQh0dPJTc8UenLRk8vqnKInV67Hv7nUDXdSC3V79TYTPbm8+ek42s3fLzoeRyrXRE+uVGqf2Bc9uVJvrmKP9ORSoZh+c9CTS++aEnJ5rE8QIZfZZ+u5lU/IaRtqykrja1PVTdL7R/lyb+7v5POi1RWCB6dhFo3nY08KGs9bP0kIuRK/eBw07q44fVKB3HTr3Z5eZXYV75TCBXOlb/8hnuaJH+y0fnZ2vnPHNln67gK5vt8gO5u73sZJsKL6Si4az8eWzUktXBdc8UBCLsWgGkxHyOW8CIRczs28jiDk8ioYzHhCrmDcdd8rIZfeFSLk8lifVIZcpWeJjNv6shR9805k1YOqDpL6zc6SIXl3RRrLmw+rX1wJuUS4u2LiO8FF7045+NZhUieruuzadyjpM8NpSOLmlxCu5OJKLqsT0UmwQsiV9KnsaQMntfC0I4WDCbkUYmoylZu/XzRZemDLIORKPT0hV+rNVeyRkEuFYvrNQcild00JuRTUx++eXHM2TJe/vP+kfFr8cWS1p2edIUPb3ic3t7pDTqtRV8ERMAUCJwToycWZkG4CqkOudPPheMInQE+u8NWMFVcuQE8uzo50E6AnV7pVtPLjoSeXvrUm5FJQG79CrrlfzJKnV/xRzI8nmo9m9S4w7pT4gNzU8tcKVs0UCFQUIOTirEg3AUKudKsox0PIxTmQTgKEXOlUTY7FFCDkypzzgJBL31oTcimojeqQywy3nlz+37J+99rI6rLrXyj3X/Z7uaHFjQpWyxQIVC5AyMXZkW4ChFzpVlGOh5CLcyCdBAi50qmaHAshV2adA4Rc+tabkMtjbVT25Cr86m35v0v+IO13tYk0lD9S/6gMz31I+l18k2WvLaueSPTkusOyovTkoieXk6e6+Tn72+8Z4WRIhW1pPO+JL+FgJ32gVIdcNJ6PLY2TWvh3RjibmZ5czrzCsDU9uZxXiZ5czs28jqAnl1fBYMbTkysYd933Sk8uvStEyOWxPipCrprn1ZG/bXpWlmxZHFnNkGq3yI9btpWBXW4vW53dQIuQi5Cr/Clt9xdyGs/HvhAQcrl7YeTuiu7c7I5aVrRIatepK63a5tod4vt2hFy+E1fYgdWVXJ+uXiV7dhdLp7zuqV+QBnsk5HJeBEIu52ZeRxByeRUMZjwhVzDuuu+VkEvvChFyeayPl5DL/DjizJnjZd73c2W98VWvZn2559LfSfaWJpLTsp00zW5OyFWuPhs3rJNNG9dLt6t6VahaNKTJH0TIRch1QmD82FHS78Zb5dSsWo6f5YRcjskiAwi53LnZHUXIZVcq8XZcyaXGUadZCLmcV4OQy7mZ1xGEXF4FgxlPyBWMu+57JeTSu0KEXArq47Qn15b9X8sTy/5LXls/WY4dPyannnKq3NLqTrkn93dSv2YDBStiCgTcCdCTy50bo/QVUP1xRX2PlJVligA9uTKl0plxnPTkyow6Z9JR0ng+c6pNTy59a03IpaA2dkOu4u93ytP//h+ZsOYfcuTYEalaparccPFAeejy/5Jzap+rYCVMgYA3AUIub36M1k+AkEu/mrAibwKEXN78GK2XACGXXvVgNd4FCLm8G4ZlBkIufStFyKWgNslCrr2HSuRvHzwtYz/+mxw88n1kj1c26ymPdvqjXNTgYgUrYAoE1AgQcqlxZBZ9BAi59KkFK1EjQMilxpFZ9BAg5NKjDqxCnQAhlzpL3Wci5NK3QoRcHmuTqCdX9+v7yD8WjZKNn6+TiccmRvY0otoIye2eJ52a/yzyZy8N5e2OnT9nmuS0bi9Nzr8g5midfn/urCnSpn0HadykWcw8qr5vTrp8SaFk1aolrdt1qFAZenKdJJnyyoty1XU3SL36iT/eSuN5enJV9hJn9frh8eWQnlxeAZOMpyeXGmB6cqlx1GkWenI5rwY9uZybeR1BTy6vgsGMpydXMO6675WeXHpXiJDLY32sQq4tWzfLvAXT5Pmjz8tZB86Qi42vVfU/loc7Pib7lmyXa3v3l9Pq1iPkigvLCLm2yPL3CqVX34FJz0pCrqRENJ5PQkTIlfwcstrCbnDsbvbEowi51KgScqlx1GkWQi7n1SDkcm7mdQQhl1fBYMYTcgXjrvteCbn0rhAhl8f6lA+5zCbyZjP58e+NkVYHWsrLxlenmp3k6nrXys397pVqVapJ/Jtru1djednO6RVblW2v6oqtyuYh5CLkin86Tn11rHS/prfUb9DQ8TOVuysmJiPkcnxKRQYQcsW6ffLR+7J/X4l07NLNHWgAowi5AkD3eZeEXM6BCbmcm3kdQcjlVTCY8YRcwbjrvldCLr0rRMiloD5mT67pn86K3DFx3a41kRlPq1FX7mr/G7m97b2SdUqWgr0wBQL+C9CTy39j9pBaAXpypdabvfkvQE8u/43ZQ+oE6MmVOmv2lBoBenKlxlmHvdCTS4cqWK+BkMtjbYq+KpL75g6XVTs+iMxUs1pNGdzqDhmeO0IanOr86hOPy2E4Ap4ECLk88TFYQwFCLg2LwpI8CRByeeJjsGYChFyaFYTleBYg5PJMGJoJCLn0LRUhl8fadBvXTRZtWiRVq1SVX/4oP9J3q1Gd8zzOynAEghEg5ArGnb36J0DI5Z8tMwcjQMgVjDt79UeAkMsfV2YNToCQKzj7VO+ZkCvV4vb3R8hl38pyS7Mn15LspfLw5Y/LxQ1zItvs2Hayt1L8HQHpySVCTy7rOyKWP2+SnZY0nk8mJDSeT0JET67k55DVFvTkilWhJ5e788jLKKuQ69PVq2TP7mLplNfdy9ShHUtPLueloyeXczOvI+jJ5VUwmPH05ArGXfe90pNL7woRcnmsj9XdFQm5TqA6bVRvjlm+pFCyatWS1u06VKhMfGBYfoN9e0tk3uwpkj/oDsuKFox5VgbecpdUr1Gjws/feG2StOtwhTRq3NRy7Hf798mcGRNkwM3DKj1bVq1cKqWlpXLZ5V0q3SZZMEXIVZGOxvMiu4p3SuGCudK3/xCPr1axwwm53HESchFyuTtz1I0i5KpoScjl/Pwi5HJu5nUEIZdXwWDGE3IF4677Xgm59K4QIZfH+hByOQ+zuJKLK7nsPu0IuQi5rM4VJ1cPqf64IiEXIZfd1y+/tiPkIuRScW4RcqlQdDYHIZczL122JuTSpRJ6rYOQS696xK+GkEtBfcy7K+4/WKpgJqZAIFgBenIF68/e1QuoDrnUr5AZEXAmQE8uZ15srbcAPbn0rg+rcy5ATy7nZmEdQU8ufStHyKWgNoRcChCZQgsBQi4tysAiFAoQcinEZCotBAi5tCgDi1AkQMilCJJptBEg5NKmFL4vhJDLd2LXOyDkck13ciAhlwJEptBCgJBLizKwCIUChFwKMZlKCwFCLi3KwCIUCRByKYJkGm0ECLm0KYXvCyHk8p3Y9Q4IuVzTnRhITy56cpnnAY3nkz+R3DQGpicXPbmszix6cuUmf8KlaAsntUjRkpLuJv4ux0kHaLYBPbkqFsTN3y+alTXly6EnV8rJhZ5cqTdXsUd6cqlQTL856Mmld00JuTzWh5CLkIuQq7mtZ5GbX0IIuQi5CLlOCiwrWiS169SVVm0JuWy96FSyESGXFz09x7r5+0XPI0ndqgi5Umcd3RMhV+rNVeyRkEuFYvrNQcild00JuTzWh5CLkIuQi5DL6mXk4PcHZPqkArnp1rs9vcrsKt4phQvmSt/+QzzNEz/4zdenyyU5baRptr362dn5zh3bZOm7C+T6foPsbO56GydXD6n+uCJ3V4wtm5NauC644oGEXIpBNZiOkMt5EQi5nJt5HUHI5VUwmPGEXMG4675XQi69K0TIpXd9WB0CCCCAAAIIIIAAAggggAACCCCAgA0BQi4bSGyCAAIIIIAAAggggAACCCCAAAIIIKC3ACGX3vVhdQgggAACCCCAAAIIIIAAAggggAACNgQIuWwgsQkCCCCAAAIIIIAAAggggAACCCCAgN4ChFwu6/PVlh3Sc+BDZaPH/Pm30rlDK5ezMQyB1AlMnLlQHn92fNkO1xQWVLrz+PPc3DDR9qk7CvaEwEmBB0aOkrcWr4x8o2WLbJk8ZqQtnuhzgddvW1xslCIBN+8vyr+un3F6PVk847kUrZbdIJBcwMn7DnO2J0dPkoLJ83mPnZyWLTQUKFqxWoaOeIr3yxrWhiVljgAhl8ta5/UZLsMG95IBvbsJL2YuERmWcoHoL0/zXn1CmjY+O/JGcuWHn1UaCphvTM2HeZ6bDzNM2Lp9l+0QIeUHyA4zTsA8R18YN7vsl/r8oY9JbtsW8uCd/RNaRMcV794rhFwZd9pofcBO319EAwH+AULrsmbs4ty87zD/IS56PkcDMs7vjD2FQnPg8f9AwTkbmtKx0DQUIORyUVSrUKv8m1IXUzIEgZQIxIda8W8+ky0iPlBItj0/R8BvgfhQy845ar6GP/KnlyLBWE7XIYRcfheJ+W0LuHl/YZ7D0X+4sL0jNkQgRQJO33c43T5Fh8FuELAtQDBrm4oNEfBNgJDLBa3VL1F2rx5wsTuGIKBMwLwSy3w889g9ZXM6+SU/2ZVfyhbKRAjYFIj/B4ZkV9bG/9zJ+W9zSWyGgGsBp+8vooGteUVi9NEjLzfmNd71YhiIgAIBN+87zNfl6MduzfcdW7cXc04rqAVTpEaAkCs1zuwFgUQChFwuzg/zL9zX314a0/PCDLkandOQv4RdeDIkdQJW56n5ZvLR+weVfSSxstVEr/rio12pqxd7Si4Qf/5GQyyrK1usrlwk5EpuzBapE3D6/sLqlynznB6Sf3XSj+ym7qjYUyYLuHnfYY7Z/u1uiYa3vO/I5DMofMdOyBW+mrHi9BMg5HJRU6f/0upiFwxBwBcBN/+iai4kGg7YCcN8WTiTIlCJgJMrueKbH5efknObU0wHAafvL6y2t3qd1+HYWENmCjh93xF/5Vb0Hy4IujLz/AnjURNyhbFqrDndBAi5XFTUTc8MF7thCALKBdz0uiDgUl4GJlQo4KYnV/ndcyWXwmIwlWcBp+8vrLYn5PJcBiZQKOD0fYdV+w+7V5wrXDZTIeBagJDLNR0DEVAmQMjlktLp3Y9c7oZhCCgVSHaXo/i7dPEvqEr5mcwHgWR3VzR/YTIfk8eMtNw7IZcPRWFKTwLJ3l/EfxzR3P667h0jH0/kY+We6Bnsg0Cy9x3x7zPi34dEAwNuruBDcZjSFwFCLl9YmRQBRwKEXI64Tm4cf5tYLqN2CcmwlAvEf2Sr/C2O499cmlcEvLV4ZYU18tGulJeNHSYQKH+etmyRHRNoEXJx6oRNINn7C6ueW+b3og9en8NW8fRfb6L3HVb/mBb/3oOAK/3PkXQ4wvjXbvOY6I+YDpXlGMIoQMgVxqqxZgQQQAABBBBAAAEEEEAAAQQQQACBGAFCLk4IBBBAAAEEEEAAAQQQQAABBBBAAIHQCxByhb6EHAACCCCAAAIIIIAAAggggAACCCCAACEX5wACCCCAAAIIIIAAAggggAACCCCAQOgFCLlCX0IOAAEEEEAAAQQQQAABBBBAAAEEEECAkItzAAEEEEAAAQQQQAABBBBAAAEEEEAg9AKEXKEvIQeAAAIIIIAAAggggAACCCCAAAIIIEDIxTmAAAIIIIAAAggggAACCCCAAAIIIBB6AUKu0JeQA0AAAQQQQAABBBBAAAEEEEAAAQQQIOTiHEAAAQQQQAABBBBAAAEEEEAAAQQQCL0AIVfoS8gBIIAAAggggAACCCCAAAIIIIAAAggQcnEOIIAAAggggAACCCCAAAIIIIAAAgiEXoCQK/Ql5AAQQAABBBBAAAEEEEAAAQQQQAABBAi5OAcQQAABBBBAAAEEEEAAAQQQQAABBEIvQMgV+hJyAAgggAACCCCAAAIIIIAAAggggAAChFycAwgggAACCCCAAAIIIIAAAggggAACoRcg5Ap9CTkABBBAAAEEEEAAAQQQQAABBBBAAAFCLs4BBBBAAAEEEEAAAQQQQAABBBBAAIHQCxByhb6EHAACCCCAAALhFsjpOqTCAQzJv1oevLO/5PUZHvnZ4hnPxWzz1ZYd0nPgQ9IjL1eeeeyeyM8SzZM/9DH55LNNlUKdcXq9yD7M/RXv3lthu+jPzR9Etxnz599K5w6tYraNrsHqZ+aGT46eJAWT5ycs2NinHpRf//bJmG3WFBZE/vzAyFHy1uKVEvUpv1H0GMv/zM7xhPvsYfUIIIAAAggggMBJAUIuzgYEEEAAAQQQCEQgGlTFBzbm9wfd+z9lwZYZHMVvYwY65mPymJFid57oQUa3twqizFDouu4dIwFbZY9ocFQ+YDO3nThzoTz+7PjIsMpCrvg5E+3PPMZG5zQsC/HMsdGQq3zoZn4/ekzm/8eHXMmOJ5Dis1MEEEAAAQQQQMAHAUIuH1CZEgEEEEAAAQSSC0SvaopepVTZiGh4NO/VJ6Rp47PLwqTon+3OE51fRchlBkfmFVnRNZhzm4FV+1YXRa608jPkMvf1weoNMmxwLxnQu1vksMzwy3zEX+VlJ7RLXim2QAABBBBAAAEEwiFAyBWOOrFKBBBAAAEE0k4gPrxKdIDlr9wyr+x69P5BZQGPk3nMfagIucyAacWqtZElmx+XLFqxWh7500vyx4dvk6EjnvI95OrQ7hJ5YdzsmKvdzMDN/AgnV3Kl3VOFA0IAAQQQQAABmwKEXDah2AwBBBBAAAEE1AvE94yy6jVVPpgy/79li+zIxxTLP+zOYyfkstOTywy5mjQ6KxJomeHSg//9gvTu2bnse35fyWUGa9Gw75tt38rW7cWRsC3+o5305FJ/zjIjAggggAACCOgrQMilb21YGQIIIIAAAhkjUL6fVWVBVvRjieU/IhgPZGceVVdymR8VjF5htv3b3ZGrqswrulJxJZcZaJkeKz/8LNJQP2piFXLRkytjnkYcKAIIIIAAAhkvQMiV8acAAAgggAACCOglEA2q4q+GcvqxxMrmURlyRUOt6BVoqQy5zKqZoVb5K9sIufQ6l1kNAggggAACCKRWgJArtd7sDQEEEEAAAQT+I2A2SzevSIp/VBYUVRZyOZ1HZchlrt38SKB5FZf5SHXIZR5732vzpHOHVpH9E3Lx9EIAAQQQQACBTBYg5Mrk6nPsCCCAAAIIBChgBjTm3QDLN5GPhkbnnHl6hb5biUIuJ/OoDrnKE6Y65IovHyFXgCc0u0YAAQQQQACBwAUIuQIvAQtAAAEEEEAgcwWifbbKC/TIy7W8wivRxxWdzJMs5LLbeN7syRX/UBFyRddXfu41hQWRP5rBoPmwugLO/D6N5zP3ucSRI4AAAggggIAIIRdnAQIIIIAAAggggAACCCCAAAIIIIBA6AUIuUJfQg4AAQQQQAABBBBAAAEEEEAAAQQQQICQi3MAAQQQQAABBBBAAAEEEEAAAQQQQCD0AoRcoS8hB4AAAggggAACCCCAAAIIIIAAAgggQMjFOYAAAggggAACCCCAAAIIIIAAAgggEHoBQq7Ql5ADQAABBBBAAAEEEEAAAQQQQAABBBAg5OIcQAABBBBAAAEEEEAAAQQQQAABBBAIvQAhV+hLyAEggAACCCCAAAIIIIAAAggggAACCBBycQ4ggAACCCCAAAIIIIAAAggggAACCIRegJAr9CXkABBAAAEEEEAAAQQQQAABBBBAAAEECLk4BxBAAAEEEEAAAQQQQAABBBBAAAEEQi9AyBX6EnIAB+MD4wAAAuRJREFUCCCAAAIIIIAAAggggAACCCCAAAKEXJwDCCCAAAIIIIAAAggggAACCCCAAAKhFyDkCn0JOQAEEEAAAQQQQAABBBBAAAEEEEAAAUIuzgEEEEAAAQQQQAABBBBAAAEEEEAAgdALEHKFvoQcAAIIIIAAAggggAACCCCAAAIIIIAAIRfnAAIIIIAAAggggAACCCCAAAIIIIBA6AUIuUJfQg4AAQQQQAABBBBAAAEEEEAAAQQQQICQi3MAAQQQQAABBBBAAAEEEEAAAQQQQCD0AoRcoS8hB4AAAggggAACCCCAAAIIIIAAAgggQMjFOYAAAggggAACCCCAAAIIIIAAAgggEHoBQq7Ql5ADQAABBBBAAAEEEEAAAQQQQAABBBAg5OIcQAABBBBAAAEEEEAAAQQQQAABBBAIvQAhV+hLyAEggAACCCCAAAIIIIAAAggggAACCBBycQ4ggAACCCCAAAIIIIAAAggggAACCIRegJAr9CXkABBAAAEEEEAAAQQQQAABBBBAAAEECLk4BxBAAAEEEEAAAQQQQAABBBBAAAEEQi9AyBX6EnIACCCAAAIIIIAAAggggAACCCCAAAKEXJwDCCCAAAIIIIAAAggggAACCCCAAAKhFyDkCn0JOQAEEEAAAQQQQAABBBBAAAEEEEAAAUIuzgEEEEAAAQQQQAABBBBAAAEEEEAAgdALEHKFvoQcAAIIIIAAAggggAACCCCAAAIIIIAAIRfnAAIIIIAAAggggAACCCCAAAIIIIBA6AUIuUJfQg4AAQQQQAABBBBAAAEEEEAAAQQQQICQi3MAAQQQQAABBBBAAAEEEEAAAQQQQCD0AoRcoS8hB4AAAggggAACCCCAAAIIIIAAAgggQMjFOYAAAggggAACCCCAAAIIIIAAAgggEHoBQq7Ql5ADQAABBBBAAAEEEEAAAQQQQAABBBAg5OIcQAABBBBAAAEEEEAAAQQQQAABBBAIvcD/B5qAnIQTRoWlAAAAAElFTkSuQmCC", "text/html": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "uc.plot_history(show_intervals=True)" ] }, { "cell_type": "code", "execution_count": 10, "id": "23c4b3ba", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "0: A <-> B\n", "Current concentrations: [A] = 36.03 ; [B] = 53.97\n", "1. Ratio of reactant/product concentrations, adjusted for reaction orders: 1.49782\n", " Formula used: [B] / [A]\n", "2. Ratio of forward/reverse reaction rates: 1.5\n", "Discrepancy between the two values: 0.1456 %\n", "Reaction IS in equilibrium (within 1% tolerance)\n", "\n" ] }, { "data": { "text/plain": [ "True" ] }, "execution_count": 10, "metadata": {}, "output_type": "execute_result" } ], "source": [ "# Verify that the reaction has reached equilibrium\n", "uc.is_in_equilibrium()" ] }, { "cell_type": "code", "execution_count": null, "id": "7f59733f", "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3 (ipykernel)", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.9.13" } }, "nbformat": 4, "nbformat_minor": 5 }