{ "cells": [ { "cell_type": "markdown", "id": "5cbc8640", "metadata": {}, "source": [ "### A minimalist demonstration for the reaction `A <-> B`,\n", "with 1st-order kinetics in both directions, taken to equilibrium.\n", "\n", "\"No frills!\" For advanced graphics, analysis, diagnostics, fine-tuning, etc, please see other experiments.\n", "\n", "LAST REVISED: Dec. 3, 2023" ] }, { "cell_type": "code", "execution_count": 1, "id": "3792d78d-7429-4221-a263-57b07d77b5bc", "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": "a29db1c7", "metadata": { "tags": [] }, "outputs": [], "source": [ "from src.modules.reactions.reaction_dynamics import ReactionDynamics" ] }, { "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": 3, "id": "78077d8c", "metadata": { "tags": [] }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Number of reactions: 1 (at temp. 25 C)\n", "0: A <-> B (kF = 3 / kR = 2 / delta_G = -1,005.1 / K = 1.5) | 1st order in all reactants & products\n", "Set of chemicals involved in the above reactions: {'B', 'A'}\n" ] } ], "source": [ "# Instantiate the simulator and specify the chemicals\n", "dynamics = ReactionDynamics(names=[\"A\", \"B\"])\n", "\n", "# Reaction A <-> B , with 1st-order kinetics in both directions\n", "dynamics.add_reaction(reactants=\"A\", products=\"B\", \n", " forward_rate=3., reverse_rate=2.)\n", "\n", "dynamics.describe_reactions()" ] }, { "cell_type": "code", "execution_count": 4, "id": "9fc3948d", "metadata": {}, "outputs": [], "source": [ "# Set the initial concentrations of all the chemicals, in their declared index order\n", "dynamics.set_conc([80., 10.])" ] }, { "cell_type": "markdown", "id": "987af2c5", "metadata": { "tags": [] }, "source": [ "## Run the reaction" ] }, { "cell_type": "code", "execution_count": 5, "id": "43735178-313b-48cf-a583-5181238feac3", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "* INFO: the tentative time step (0.1) leads to a least one norm value > its ABORT threshold:\n", " -> will backtrack, and re-do step with a SMALLER delta time, multiplied by 0.4 (set to 0.04) [Step started at t=0, and will rewind there]\n", "* INFO: the tentative time step (0.04) leads to a least one norm value > its ABORT threshold:\n", " -> will backtrack, and re-do step with a SMALLER delta time, multiplied by 0.4 (set to 0.016) [Step started at t=0, and will rewind there]\n", "* INFO: the tentative time step (0.016) leads to a least one norm value > its ABORT threshold:\n", " -> will backtrack, and re-do step with a SMALLER delta time, multiplied by 0.4 (set to 0.0064) [Step started at t=0, and will rewind there]\n", "Some steps were backtracked and re-done, to prevent negative concentrations or excessively large concentration changes\n", "47 total step(s) taken\n" ] } ], "source": [ "dynamics.single_compartment_react(initial_step=0.1, target_end_time=1.) # Using defaults for all other parameters" ] }, { "cell_type": "code", "execution_count": 6, "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", "
SYSTEM TIMEABcaption
00.00000080.00000010.000000Initial state
10.00640078.59200011.408000
20.00960077.91052812.089472
30.01344077.10584612.894154
40.01804876.15876713.841233
50.02357875.04845814.951542
60.02910773.96884616.031154
70.03463772.91908317.080917
80.04016671.89834418.101656
90.04569670.90582719.094173
100.05233269.74773520.252265
110.05896768.62806721.371933
120.06560367.54554622.454454
130.07223866.49894023.501060
140.07887465.48705824.512942
150.08550964.50874925.491251
160.09347263.37372626.626274
170.10143462.28389327.716107
180.10939761.23744928.762551
190.11736060.23266829.767332
200.12532259.26788930.732111
210.13487758.15624931.843751
220.14443357.09771732.902283
230.15398856.08975833.910242
240.16354355.12995534.870045
250.17500954.03321835.966782
260.18647552.99935737.000643
270.19794152.02476937.975231
280.21170150.92231239.077688
290.22546049.89570040.104300
300.23922048.93971741.060283
310.25573147.87145942.128541
320.27224246.89139343.108607
330.29205645.81240744.187593
340.31186944.84031445.159686
350.33564643.78936546.210635
360.35942242.86335547.136645
370.38795341.88424548.115755
380.42219140.87692749.123073
390.45642940.04205049.957950
400.49751439.21170450.788296
410.54681738.41997951.580021
420.60598037.70411352.295887
430.67697537.09919152.900809
440.76217036.63096553.369035
450.86440436.30843653.691564
460.98708436.11924153.880759
471.13430036.03147053.968530
\n", "
" ], "text/plain": [ " SYSTEM TIME A B caption\n", "0 0.000000 80.000000 10.000000 Initial state\n", "1 0.006400 78.592000 11.408000 \n", "2 0.009600 77.910528 12.089472 \n", "3 0.013440 77.105846 12.894154 \n", "4 0.018048 76.158767 13.841233 \n", "5 0.023578 75.048458 14.951542 \n", "6 0.029107 73.968846 16.031154 \n", "7 0.034637 72.919083 17.080917 \n", "8 0.040166 71.898344 18.101656 \n", "9 0.045696 70.905827 19.094173 \n", "10 0.052332 69.747735 20.252265 \n", "11 0.058967 68.628067 21.371933 \n", "12 0.065603 67.545546 22.454454 \n", "13 0.072238 66.498940 23.501060 \n", "14 0.078874 65.487058 24.512942 \n", "15 0.085509 64.508749 25.491251 \n", "16 0.093472 63.373726 26.626274 \n", "17 0.101434 62.283893 27.716107 \n", "18 0.109397 61.237449 28.762551 \n", "19 0.117360 60.232668 29.767332 \n", "20 0.125322 59.267889 30.732111 \n", "21 0.134877 58.156249 31.843751 \n", "22 0.144433 57.097717 32.902283 \n", "23 0.153988 56.089758 33.910242 \n", "24 0.163543 55.129955 34.870045 \n", "25 0.175009 54.033218 35.966782 \n", "26 0.186475 52.999357 37.000643 \n", "27 0.197941 52.024769 37.975231 \n", "28 0.211701 50.922312 39.077688 \n", "29 0.225460 49.895700 40.104300 \n", "30 0.239220 48.939717 41.060283 \n", "31 0.255731 47.871459 42.128541 \n", "32 0.272242 46.891393 43.108607 \n", "33 0.292056 45.812407 44.187593 \n", "34 0.311869 44.840314 45.159686 \n", "35 0.335646 43.789365 46.210635 \n", "36 0.359422 42.863355 47.136645 \n", "37 0.387953 41.884245 48.115755 \n", "38 0.422191 40.876927 49.123073 \n", "39 0.456429 40.042050 49.957950 \n", "40 0.497514 39.211704 50.788296 \n", "41 0.546817 38.419979 51.580021 \n", "42 0.605980 37.704113 52.295887 \n", "43 0.676975 37.099191 52.900809 \n", "44 0.762170 36.630965 53.369035 \n", "45 0.864404 36.308436 53.691564 \n", "46 0.987084 36.119241 53.880759 \n", "47 1.134300 36.031470 53.968530 " ] }, "execution_count": 6, "metadata": {}, "output_type": "execute_result" } ], "source": [ "dynamics.get_history() # The system's history, saved during the run of single_compartment_react()" ] }, { "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": 7, "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": "red", "dash": "solid" }, "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" }, "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.00047820420640189593, 1.134778581791699 ], "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": "iVBORw0KGgoAAAANSUhEUgAABVsAAAFoCAYAAABJ+CtsAAAAAXNSR0IArs4c6QAAIABJREFUeF7snQmYFOWduH/Tc3MNlxwDwiiIcooQIlkvRFREJMZ7xWggKqLZJOriavYwZrPRx5NNjIquBqNZr4SNMd6iiP4VBQUFwajcckbua2a6e/r/fNVdM9Xd1d1VXdXV1d1v5zHDVH/n+31VXf3OV7+vLBKJRIQXBCAAAQhAAAIQgAAEIAABCEAAAhCAAAQgAAEIOCJQhmx1xI/MEIAABCAAAQhAAAIQgAAEIAABCEAAAhCAAAQ0AshWJgIEIAABCEAAAhCAAAQgAAEIQAACEIAABCAAARcIIFtdgEgREIAABCAAAQhAAAIQgAAEIAABCEAAAhCAAASQrcwBCEAAAhCAAAQgAAEIQAACEIAABCAAAQhAAAIuEEC2ugCRIiAAAQhAAAIQgAAEIAABCEAAAhCAAAQgAAEIIFuZAxCAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEXCCAbHUBIkVAAAIQgAAEIAABCEAAAhCAAAQgAAEIQAACEEC2MgcgAAEIQAACEIAABCAAAQhAAAIQgAAEIAABCLhAANnqAkSKgAAEIAABCEAAAhCAAAQgAAEIQAACEIAABCCAbGUOQAACEIAABCAAAQhAAAIQgAAEIAABCEAAAhBwgQCy1QWIFAEBCEAAAhCAAAQgAAEIQAACEIAABCAAAQhAANnKHIAABCAAAQhAAAIQgAAEIAABCEAAAhCAAAQg4AIBZKsLECkCAhCAAAQgAAEIQAACEIAABCAAAQhAAAIQgACylTkAAQhAAAIQgAAEIAABCEAAAhCAAAQgAAEIQMAFAshWFyBSBAQgAAEIQAACEIAABCAAAQhAAAIQgAAEIAABZCtzAAIQgAAEIAABCEAAAhCAAAQgAAEIQAACEICACwSQrS5ApAgIQAACEIAABCAAAQhAAAIQgAAEIAABCEAAAshW5gAEIAABCEAAAhCAAAQgAAEIQAACEIAABCAAARcIIFtdgEgREIAABCAAAQhAAAIQgAAEIAABCEAAAhCAAASQrcwBCEAAAhCAAAQgAAEIQAACEIAABCAAAQhAAAIuEEC2ugCRIiAAAQhAAAIQgAAEIAABCEAAAhCAAAQgAAEIIFuZAxCAAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEXCCAbHUBIkVAAAIQgAAEIAABCEAAAhCAAAQgAAEIQAACEEC2MgcgAAEIQAACEIAABCAAAQhAAAIQgAAEIAABCLhAANnqAkSKgAAEIAABCEAAAhCAAAQgAAEIQAACEIAABCCAbGUOQAACEIAABCAAAQhAAAIQgAAEIAABCEAAAhBwgQCy1QWIFAEBCEAAAhCAAAQgAAEIQAACEIAABCAAAQhAANnKHIAABCAAAQhAAAIQgAAEIAABCEAAAhCAAAQg4AIBZKsLECkCAhCAAAQgAAEIQAACEIAABCAAAQhAAAIQgACylTkAAQhAAAIQgAAEIAABCEAAAhCAAAQgAAEIQMAFAshWFyBSBAQgAAEIQAACEIAABCAAAQhAAAIQgAAEIAABZCtzAAIQgAAEIAABCEAAAhCAAAQgAAEIQAACEICACwSQrS5ApAgIQAACEIAABCAAAQhAAAIQgAAEIAABCEAAAshW5gAEIAABCEAAAhCAAAQgAAEIQAACEIAABCAAARcIIFtdgEgREIAABCAAAQhAAAIQgAAEIAABCEAAAhCAAASQrcyBnBCY99JC+fc7H5P/vGm6nDfp5JzUUQqF3jvnWXn0qZdk7uybZczIY0qhy/QRAhCAAAQgAAEIQAACEIAABCAAAQgULAFka8EOnf8broTrsy8skAfvuF661HX0f4N92kIlXLdu3ym3zZoutTVVPm1lW7PUuD/w+PMy584bZUD/et+3lwZCAALpCahr0IfLPi/Ja/niZZ/LD356h5x92tiCuQaX+nzmMyj7GbB6/WaZcdM9cvstV+XkD7xejo3Veye9z9de8V3fLQ7Q29a9a11JXn+zn8nkhAAEIAABCOSfALLVZAwONTbLrXc9Ji/OX5T07g//cZLcMOOivI9cuptDdYP50psf5F126Rx79eiakZneny3bdhTlath0c6p3z25px2rXnn0y8+b75KJzxuX1i0CqPiRKCC+/TOX9RMzQAH1eTxp/fMZzwO99yVf7cv1FuBCupW6zN5uX6c5bp7JVv4Z9e+QxSeeBn8+RVO0uxjlTTNftYuqLOve9vKdDtkafxNKfzvLDU0X6dahffQ/+4OP2hyHlQQACEIAABHJIANmaRraqt4yrCfWb0FHDjsr7DU+xfdnTvxx179JJ+vXp6SrfbL/M6yua3BDsuqhMnFNef5HK9lqSanWX3q+PV3zZKoyL7YtutsxUvmznnpM6iy0vstX9EUW2WmOaSjIX2+evLpeK5YmEYvsMQrZaO18TUzlZ2eon2ar6pd+D+UH+Zjca5IIABCAAAQiUHgFkqw3Zqn8hUbFI833Dk2sB4eWpYFwB23B4L9ceQTeuls1WmLpRhmKZTrbqN9F+jW9rRRi+uuBDGXhEXy1sQLF90fXyXKCuZAK5vtbluvxCGdN8rWz1K59086IY50wxXbeLqS9enx+sbPXnHgPp7iG9niPUBwEIQAACEICANQLIVpdkq9kj1onyTJdqiVWmErf6X9aN6VXarl06aTG11CP3xpcuFFOtxkkszywGnZ73V7dcJT+7/RFZvmqNVkUu49UZ/2Kv981J7Cy3BKmRrdMwB1Zka74FfqpLht1VNfoX3dm/+JH8/tlXW8NxDB98ZFLMMbM5bhZWwU6Zqh/G8VK/qzJTrZq2cl4klqfKtCLHzaSM3b6kGhez60lim6z0ze45b6XexDSJY298PPuk40doMTH1l7EPZtxVOv1aZ5Qqz7/yrraZnHqpc0m9jOXq5RvPs0zle3UtzXZ+mbXPbAVU4uPwifPS7DzUOaqN+ezOEeOc1evWP0v09/TPlM3bvtE+z4zXfH1F2uUXnSk//Y/7Wz/r9DwrPl8TN7apPp+szP9U51cqYWd3zhjDwKi61B9rE+e6vhGi2TmgH7Nyj5Huts/svNXPhUzjr8o1G8fEzyw78yTbOa/akq4v6n2719hMfUv1+W02R8zOLbVi2MrnYbrPYGPMZKvXz3TzQX8vce4NP+YI+WbX3riYrVY/p1WZSeUlfO67PTZ6P+x85idySfc5bYw9b2d+m81Ts3sgO+eBmve33P5I3kOEWZlXpIEABCAAAQhAQATZajILMj3ynbhRiNnKP7NYb+pG6Z0PPo2LWWf2xVivf8Pm7XFyStUz/52P5OrLzmmVSWZS0uwLeOKxVHXoN8rGm8J08fbcOImMbauprtbi5aqX3Q2hciFZU92U240tm2pO+X21QjaxwvQvZkZpmu7LqmJ83qS21SRmX2DtlGm2Utis/nTngPEcN/siZvWRvlRf4pRwscIn3ZfvxLjMaqx+9/TLMvOKc7WN1HJxzpuJd7N6E9uWOKZGuWFcdW7GNd0qQqMESBQ/Vq+3mR4JT7zeu83Vyfwy+/KtX8ONXBNXq1mVC4lCJtvPBSsxWxNlqxLnRolqvL6bHTfmt3pupzq/Ml2X7cwZ41w3/iHBOHcTjyc+zm/1HiNVf8zOq8RzMt1qULNrqhkDs/sHM5ZO57z6I4rxfDfri9VrbLq2JI6LccNPo/w2pks8J+18dtmVreoPGJmun5nmeOJ9Zqq5YuVzOtW10Xg/Z4eH1Xln9TM/FQur10M798dmn5WJ89TueZDuuuPGvThlQAACEIAABCDgLgFkqw3Zmupx71Rxoaw8yma2iZSVfHa+7KUSQ2ZfgNOt5DJ+0XBrGqaS0olfptLVZ/zSk224ALv9MX7xt7IiNd0GWapur9ptt5/ZiPZU89fqqgyzDcGslpluU7bE8zRVe1KtTDKucFEcVbpDjU0y7OgjUmK1+iVOFWCVjxXRm4tz3kq9iVJPB5M4j1LNK7PxyyRb7cSZtFt+4vUwF1zTraDMNL8SzxUjV6MkVnUYr9925qUaQ6efC9nI1q3bdyb9wS3VZ22253aqEzfTdc/O52+qDQ7tHHdyj5Fu/Fb8ba3U1lSnDf+S6ZpqnGep5knieeNkzqeqw0pfEq+x6aR6KjF2+y1XiVrtrebAv97+iAw5+gjZv/9g61xNHCurn13pPpsT+2zn+plOwJutkkx1DU8sJ921J3ETWRVm6OSxI7U/AlrlYXXe6X+gN9uINdV5k9gXq9dDq9dBq5+Dds+DdEzs3tuRHgIQgAAEIACB3BNAtqaRrS/OXxT3rtmjiqm+MOkyRj0eqd+cq2OpHqk0yjYrN4h2vuwlftHWO2X2RSPdzaQdqWF16lp55DVTWdnI1lTy08qj4caxVStc7chWlTdxxa4+J9R7D95xvXSp65ipy569n0k6mDUk3ZcpM4luFNfG8jKt9tIFpbHMTOeFUeCkOs8Sv9Doc9TsEcBMA2H1S5xZX1KVnep8NqbPxTnvpF7VNiPvxqYmmXnzfZK4O73VVXB6XzP9YcrK9TbTnEknLd24ljqZX2Zz9Z45z4oeBubGGRdpYihxrtuZl/rYJa7wVccz8df5eC1brZ7bqc6xdJ/rxuu/lSdL7EhV4z3CReeM01b8273HSHVNVis904UDSid+Eu9j9DoS5WWq+4fE+eZkzusrI7PpS+L9Rjq2mdKqdiz6aKWoUBe/vO/38l+3XCX1PbtrT+YYxZ/dz0Oz8bMrW83uMxLLzTRWxvvWxHses8/pVCt9E+u1yiOd9DXOO1V+YhgSvU4r99Kpzmezdlq9P073WWlskx4Oxeq9BbI10x0X70MAAhCAAAT8RQDZajIeZl/49Rv8RBmXShQZi02Mi2YUq4k3T5keX9TLtSMIUt0gmn2Jtnoz6dY01h/LMivP6g1oIhMlQHO1UtQ43lbFrGpfpnG1smrQLeZ2ysk2jICZmDfroz7+RpZ2V7YaZWs6jsYvOYqB+lKc+AcVIxuzxzNTvZ+KqR2pZXUOpDufjV8yzeSYk3Pear163FQzJrog8UK26tfsdNfbVF+0U3HM1bVUH3u780ulN36xV2Ek1EutLFNtVS8VVkLN9bGjh7SG67AzL83mjN5OP8pWu+e22Tz1k2y1c4+R7tqunw/GNHb+oJWqbGM4lFTnhxlPp3PeGPtWtc1OX/R7snT3UWbv6Z8ht/x4qtz+6z9o59RZ48e2nl+nnnCc9kck/Y8c+vlp9fMwFWO3ZWu6+xEzyWn1c9rsj9iJMditylaz+WHko5e7c9deLX6z2R+98yVb093Xqj4Y/1Bg5zxAttq5eyUtBCAAAQhAIP8EkK0mY5DqRtTqDWe6L2+pVnIZV0JYuUG0I1tzscrNjambrg9WxZNZO4xfTt2SrtlKVr19VmWrHYHrxhhYKUPNx8QYnJm+1Fv5cplq/J3I1nSrYRLPKyvnWbovv0oqZhovO1LL6px3ssLUyWp2J/Umcsz0GKxKr68AzyaMQKbyjddbP15LzT5rUs1FfZWXWommVrXqokcd139XjwsbV6rZmZeFJlv1eMVmYQisXO9Umkwr+u3MmVyubLXan8R0Rilm/GOw2XXb6mPl6eZJOl562+zMeWN/7PQl02pVY7lm12P9GqjOsXseekZbzTqgf732Bw+1ynXy6d+RX9z3+7gNjKzKxXRj6bZs1cfK7BxJHG87n9OJfdDzdu9a1/rUjlUeVuednc/8VIytXg+tLkaw8lmZzb0FsjXbKx75IAABCEAAAvkhgGw14Z5KjJnd1GeSaHrxqW5YzW6eUt2MqrQLFy2TM8d9O+2jhYk3hG7FGXQ7jEC6VVHZrKhMdaPvZKWrW+I20zyxKtrycZnQGUwaf3zc5m7GtliNl2dlFaoT2ZpOkiTKVaur8hYu+kTUKmtjeIdMMibdeW/1y2Y6uZZqJY8eG09/PDFxtY9Zu61+gUw3RzPV60S2pltlmOnR58THvM2ut+nK9+Ja6mR+GcVgXcf2smffgVapYRzrVJs6GvmkixlsdY6kmrPpJEG6FYSJIVdS/YEk23M7U3vV+2YbNdqZM05la6bPDiufCcaYmamuTanG3+q1Lp1sTSzbyZy30her19h0bM3K0OeqEqzqfDP+QcgsfqtiYrUtXsvWTNdO/Y8zme4fjeEu1IZdJ489Nq4rifVY5WF13tn5zPdKtlq9n7N7Hlj5o4WV6wFpIAABCEAAAhDwhgCy1YZsNX6x/WbnntbVC/oN0KhhR8V9MVM38nc9+LRMPW9Caywv45e3VLFG9ePGXWLthBsw+2KcuDrRrI50X5asiimr09bKF0i36rQiC83ard8wu7E61spje4nzxyrLTMLg4xVfxq20yaZcnUVinDy9X8Y6nHyZMsbYzOax0FRfbvVHaI3tN2u7zkalVy8VM9GsP1a/TFldMaPqslqmfp4mrjY2W7lrTOPGOW+2ytmsXrXqN1H0qv6988GnmrDPtPLU7DppPGYcJ7M/Apmdb5mut2ble3EtdTK/dA76qsDEa1Wq42bz0s5qzUz8za4xqYRtLmSr1XM73bUw3edPuut5qlWIupTS67QjYa3cYygBmOplxj6xf1aeNEmcX6oP6nF69Vi9+mOUWT1mn79O5ryVvlj9DFK8zNqnX4sTn1wwXkOM76U6nurzKJvrvfEPJnaun6nmhFkZxs9f/fqdKZ3OwSyd2XliZ2xS3YMlzjuzMs0+81OxsPo5beePTvq1N93noN3zIN0fxLK5ryMPBCAAAQhAAAK5JYBsNeGbSQQaVzvqN1LGm21jkWbxAvU4kSrm1Oxf/Eh+/+yrcRsqJH6B1n9PFF3GG2OVRq8r3Q2hMc6Z2QYTdm4mnUzNVF9mjGUW01/xU80Pvb+ZHkfPhnU64eCkvEwbx9n5MpU4h9UKUn1zH6OcsFOm/gXXONcV33Ubt4rZY5P6lzIjk8R4wfoXJz1NYhw6p1/i7H75NuujWYzjxL65cc4nlmlWr1kcOiMzu7Ig1bXOihCzcr3N97U02/mlzzt9TMy+2KtV0KlijSeu/E0cN708Nz4XEq+B+lzcvO2bpA1urK5gNX5WZntu2zl3jWmtzhk7UlWVnyq9lXuMVH0xy2vlvDXOp8T+mn12Jc7jdJ9v2c55K32x+3lh1rdUG1+a/cHJeD1OzGe3LWZjmErgW9lgMN1nvdnnrwqRkBh2xOrntPHeWK83UdDb5WFl3tn9zE9kkgvZavxMN9aXeO9g9TzI9L0km3s68kAAAhCAAAQgkFsCyNbc8qV0CEDABwTSPcbsg+bRBAhAwIcEUklmHzbVF02Cly+GgUYUIQE7T74UYffpEgQgAAEIQKAgCSBbC3LYaDQEIJCKgPFxdT0Nj98xXyAAAbsErMaNtFtusaZHthbryNKvfBLgOpRP+tQNAQhAAAIQyJ4AsjV7duSEAAR8SMDsEXazR2Z92HSaBAEI+IxAqnjVPmumL5qDbPXFMNCIIiKghzjo3rWudfPDIuoeXYEABCAAAQgUNQFka1EPL52DAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEvCKAbPWKNPVAAAIQgAAEIAABCEAAAhCAAAQgAAEIQAACRU0A2VrUw0vnIAABCEAAAhCAAAQgAAEIQAACEIAABCAAAa8IIFu9Ik09EIAABCAAAQhAAAIQgAAEIAABCEAAAhCAQFETQLYW9fDSOQhAAAIQgAAEIAABCEAAAhCAAAQgAAEIQMArAshWr0hTDwQgAAEIQAACEIAABCAAAQhAAAIQgAAEIFDUBJCtRT28dA4CEIAABCAAAQhAAAIQgAAEIAABCEAAAhDwigCy1SvS1AMBCEAAAhCAAAQgAAEIQAACEIAABCAAAQgUNQFka1EPL52DAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEvCKAbPWKNPVAAAIQgAAEIAABCEAAAhCAAAQgAAEIQAACRU0A2VrUw0vnIAABCEAAAhCAAAQgAAEIQAACEIAABCAAAa8IIFu9Ik09EIAABCAAAQhAAAIQgAAEIAABCEAAAhCAQFETQLYW9fDSOQhAAAIQgAAEIAABCEAAAhCAAAQgAAEIQMArAshWr0hTDwQgAAEIQAACEIAABCAAAQhAAAIQgAAEIFDUBJCtRT28dA4CEIAABCAAAQhAAAIQgAAEIAABCEAAAhDwigCy1SvS1AMBCEAAAhCAAAQgAAEIQAACEIAABCAAAQgUNQFka1EPL52DAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEvCKAbPWKNPVAAAIQgAAEIAABCEAAAhCAAAQgAAEIQAACRU0A2VrUw0vnIAABCEAAAhCAAAQgAAEIQAACEIAABCAAAa8IIFu9Ik09EIAABCAAAQhAAAIQgAAEIAABCEAAAhCAQFETQLYW9fDSOQhAAAIQgAAEIAABCEAAAhCAAAQgAAEIQMArAshWr0hTDwQgAAEIQAACEIAABCAAAQhAAAIQgAAEIFDUBJCtRT28dA4CEIAABCAAAQhAAAIQgAAEIAABCEAAAhDwigCy1SvS1AMBCEAAAhCAAAQgAAEIQAACEIAABCAAAQgUNQFka1EPL52DAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEvCKAbPWKNPVAAAIQgAAEIAABCEAAAhCAAAQgAAEIQAACRU0A2VrUw0vnIAABCEAAAhCAAAQgAAEIQAACEIAABCAAAa8IIFu9Ik09EIAABCAAAQhAAAIQgAAEIAABCEAAAhCAQFETQLYW9fDSOQhAAAIQgAAEIAABCEAAAhCAAAQgAAEIQMArAshWr0hTDwQgAAEIQAACEIAABCAAAQhAAAIQgAAEIFDUBJCtRT28dA4CEIAABCAAAQhAAAIQgAAEIAABCEAAAhDwigCy1SvS1AMBCEAAAhCAAAQgAAEIQAACEIAABCAAAQgUNQFka1EPL52DAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEvCKAbPWKNPVAAAIQgAAEIAABCEAAAhCAAAQgAAEIQAACRU0A2VrUw0vnIAABCEAAAhCAAAQgAAEIQAACEIAABCAAAa8IIFu9Ik09EIAABCAAAQhAAAIQgAAEIAABCEAAAhCAQFETQLYW9fDSOQhAAAIQgAAEIAABCEAAAhCAAAQgAAEIQMArAshWr0hTDwQgAAEIQAACEIAABCAAAQhAAAIQgAAEIFDUBJCtRT28dA4CEIAABCAAAQhAAAIQgAAEIAABCEAAAhDwigCy1SvS1AMBCEAAAhCAAAQgAAEIQAACEIAABCAAAQgUNQFka1EPL52DAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEvCKAbPWKNPVAAAIQgAAEIAABCEAAAhCAAAQgAAEIQAACRU0A2VrUw0vnIAABCEAAAhCAAAQgAAEIQAACEIAABCAAAa8IIFu9Ik09EIAABCAAAQhAAAIQgAAEIAABCEAAAhCAQFETQLYW9fDSOQhAAAIQgAAEIAABCEAAAhCAAAQgAAEIQMArAshWr0hTDwQgAAEIQAACEIAABCAAAQhAAAIQgAAEIFDUBJCtRT28dA4CEIAABCAAAQhAAAIQgAAEIAABCEAAAhDwigCy1QXS+w6FZN/BoAslUQQESo9Azy418vc9TdLSEim9ztNjCLhAoHtdtew5EJRgqMWF0igCAqVHoEvHKmlsCsuh5nDpdZ4eQ8AFAp3aVUpLJCL7D4VcKI0iIFB6BNrXVEhFeZl2P8cLAlYJ1HertZqUdHkggGx1AboSrUq48oIABOwTQLbaZ0YOCBgJIFuZDxBwRgDZ6owfuSGAbGUOQMAZAWSrM36lmhvZ6u+RR7a6MD7IVhcgUkTJEkC2luzQ03GXCCBbXQJJMSVLANlaskNPx10igGx1CSTFlCwBZGvJDr2jjiNbHeHLeWZkq0PEt912m9ww62dJK1vXr/1K/rbyUznj7PO0Gn730L3y/St/LBUVFdrvc+fMlqnTrpXKqirt9717dsvLf3lWLv7+1XEt+uqLlbJx3Ro59YzJccffW/iG1HXuKkNHjIo7/vwfn5SxJ46Xnr3q444/cv+dctWPboo7tm3LJvngvQUy5fypSRS+3rBWli9bImdNuTDpvVRtUglXLl8qu3Z+IyeccnpSvqVL3pdQKCRjxp6U9N7i9xdqLEaOHpv03rsLXpOu3XvIkGEjTUcrXVv1DIm8Ewt6e/7L0ru+rwwaPNy0jg3rVsvKFctk4uTzU86YRx+4W6Zdc4MEAoG0s8rtdGZja9YAq+lU3kgkIqqdV143y9IZEg6H5fGHZ8v0mTdaSh8KBuWJx+6XaTOuF7uy9c3XXpD+DQNlwKDBlurSE23+er0sXbJIzj73Ylv5VOL/nfugTLngMunQoaPtvHqGTHPMSsHpzj0r+c3SLHj9RenTr0GOOnpotkUk5Xvi0fvlwkunS01tO9fKNBbU3NwkT819SK64+ic5KV8V+s32rfLOW6/K9y6+Imd1mBX85ecrZNPXG2TchEmW682lbE31WWO5cS4nTPUZ53I1rhQ375nH5ZTxE6XbYT1dKS+fhTi5fuaz3VbrTiVbV332iezYvlVOPPVMq0WVZDp13zf/lb/IBZdOL8n+e9HpQwcPyJ+eniuXTb/Oi+ps14FsbUOWr/sH24NGhpwR+P0jv5aLL79aqqtrLNeBbLWMioQxAuq7/a233goPHxNAtorI6vWbZcZN98iWbTu0ofrPm6bLeZNObh22eS8tlH+/8zHt97NPGyu3zZoutTVRSYpsjZ/dyFZkq5XrHbLVCqXkNMjWKBNka/zcQLZmdz7lOheyNdeE3Ssf2eqMJbLVGT8ruZGtVij5Iw2y1R/jkM9WIFvzSb906ka2+n+sS1627tqzT2befJ/cOOMiGTPyGEn8ffGyz+WeOc/Kg3dcL13qOsq9c57VRvWGGRchWxNW2yogyFZkq5XLHrLVCiVkaypKyFZka+LTG9mdUbnNhWxFozlqAAAgAElEQVTNLV83S0e2OqOJbHXGz0puZKsVSv5Ig2z1xzjksxXI1nzSL526ka3+H+uSl61qVeu/3v6I/NctV8mA/vVyqLFZbr3rMRk7eoi2ulXJ1YbDe7WudE2Ur2qIidnq/4lOC/1LwG4YAf/2hJZBID8EcrmyNT89olYIeEuAmK3e8qa24iNAGIHiG1N65C0Bwgh4y7tYaiNmq79HsuRlqxoeJVRfevMDmXNnNN7kXQ88Lbf/7Cqpqa6OE6/qvUQ5i2z19wSndf4ngGz1/xjRQn8TQLb6e3xonf8JIFv9P0a00N8EkK3+Hh9a538CyFb/j5EfW+g32Wq2MNGP3DK1KXEBZqb0qd5HtoqINikeeka+2bVXi9uqx2zVIV94zjgtxICpbH3xRWm5b7Y0/fNNEjrl1GzHgXwQKFkCHWoq5EBTSCKRkkWQs44rpmVlOSuegn1CoF11uTQFWyTcwknkkyGhGRYJqM0Yy3xwkaqtKpdQuEWCYc4hi0NHMgjEEaiuDGj3cc2hlqIi45drVFFBpTOmBKoqAhIoE2kMFtc5xHDnlkDH2ujm6169dD/24vxFcVXq/iyfslXts/TsCwtaw386YYJsdULPkFetVNVXsqqYrHrM1ovOGSdnjR+bcWWr2iDr1p//XCsxPGq0BP/5JglOOVdWr/5SViz/VL77vQu09359351y7T/dIBUV0RPi/v++W66e+WOpqoputLV79y6Z99zTMv2qmXE9W7XyM1m79iuZdPZ3446/Nf816dylqxw36ltxx5/6w+My7tQJ0ru+T9zxe+/6ldww62dxxzZv+loWvv2mXHLp5Uk0161bIx8v/lDOu/CSpPdStUklXLb0I9m54xsZPyF5594PFv0/CQVDcsJJpySV+e7CBVJVXSXfPv4fkt6b//or0v2wHnLsyFGmo56urXqGRN6JBb368l+l7+H9ZOiwEaZ1rF3zlSxb+rF87/xorF6z1+x77pAfX3+TBAKBtLPT7XRmY2vWAKvpVF51c6naef0/35K2L/qb4XBYm9M/ueFfLKUPBoPy0G9nyz/9dJa0r62Qg41hrU4rr5f++mc5csAgOWbwECvJW9NsWL9OPvzgPbngoktt5VOJH3nofrlk6uXSsWMn23n1DJnmmJWC0517ZvkjZWVSloHryy/9RRr6HymDhw6z0gRLaR787Wz5wbSrpbZdO0vp7SZqamqS/5lzv1z34+jTCLl4bdu2Vd549SWZerm3u2uv/Gy5qLk6cdI5lrvVrqZCGpvD0pID2Zrqs8Zy41xOmOozzuVqXCnuyccflTPOmiw9evR0pbx8FuLk+pmu3REpkzKxdu3PZf9rqsslFIpowtX4+vSTpbJ921aZcMZZuay+4Mve8c038tcX5skV064u+L74tQMHDxyQ38/9H7nmup/4sonVleUSkYg0F5koUlcnu3+zztf9gy8nRok26oHf3CvTr7pWampqLBOoVLI1UCZNzWHLeUhY2gTUd/tbb73VMwj6pvKTxh/fun+Rqlz5s1t+9YjMuvYS2blrb9x+R541zuWKkK0uAVX2/bkXFshts6ZLbU1UfBo3wcoUs1XJ1p99tkoq/+9PIqGQlj905ED59Jrr5NPDussZ55yvHfvdQ/fK96/8catsnTtntkyddq1UxmTr3j275eW/PCsXfz/+RjXV7uPvLXxD6jp3laEj4gXk8398UsaeOF4SNw9RAZSv+tFNcdS2bdkkH7y3QKacPzWJ5tcb1sryZUvkrCkXJr2Xbkd0NshigywrpyYbZFmhlJwm3bmXXYkiC15/Ufr0a5Cjjh6abRFJ+Z549H658NLpUlObG9nKBlnxyHMZRiDVZ41rk8VmQak+42wW40lyNsjyBLMrlbBBljOMbJDljJ+V3GyQZYWSP9KwQZY/xiGfrWCDrHzSL526vdwgS5ePvXp0jROtibT1la1q8/lbbn9Ee2pcvebOvrn1SXH1uy5u9fd/+I+TWsvVFz9Ou3ii/O6ZV2T5qjVaGWr17LFDB8qMm+4xLVetbF300co4r5dYz/DBR2orX9Vr5s33tZatfje2Adnq0nmkD8Dtt1ylTQDjyla1QVbiUmijiFVNULJVrRg9uHq9tP/tf0v7J+dK2YH98rdjjpGPjz9ezm4YJAemXyWPPfEwsnXJ+xIKhWTM2JOSRm/x+ws18Txy9Nik995d8Jp07d5DhgwbaTrq6cSwniFRbicW9Pb8l6V3fV8ZNHi4aR0b1q2WlSuWycTJUXlu9nr0gbtl2jXIViunJrLVCqXkNMjWKBNka/zcQLZmdz7lOheyNdeE3Ssf2eqMJbLVGT8ruZGtVij5Iw2y1R/jkM9WIFvzSb906vZStiY6s1SUlTv7wU/vkLNPG9sqPRMf70+1Qb0ucnUfp+pQYlQ9fa6Xq8tS/dg9c55tTZMoW83a/OqCD2XgEX2la+eO8runX5aZV5yrLbhMdIDIVhfPI33w9CL1mBP672rg/v3Ox7RfjRNHf3/fwaDsOxRd1RrYu1faP/awtHv4t1K+fZt2rKVTJzkw7So5cO1PpKVbdxdbTlEQKHwCbJBV+GNID/JLIJeyNb89o3YIeEOADbK84UwtxUuADbKKd2zpmTcE2CDLG87FVotXG2QpX6ZWqqoN5Qf0r0+J0Sxma6JcTXxyXBVmzKd+V6tO1epYfd8kXYamO5YoWxMXSWYae5V/3cat2gpbZGsmWh6+b5SterVlzc1S+9xT0uG3s6Xii79phyPVNXLwkqmy/6ezJHx4Pw9bSFUQ8C8BZKt/x4aWFQYBZGthjBOt9C8BZKt/x4aWFQYBZGthjBOt9C8BZKt/x8bPLSs02Vrfs7u2J1LiBluKceIj/k5kqypP1WPc6D5xHBNDDKj39YWVev6xo4eIeto921dZxOquNNnWUAL5zGSrsds1r70s7X9zn1S//270cHm5HJp8ruz/55slONi9OIklgJouFiEBZGsRDipd8pQAstVT3FRWhASQrUU4qHTJUwLIVk9xU1kREkC2FuGgetAlr2SrnTACxkf7FQLjylZdtqaTmFZWsapyE9MZV7Zmkq36k+vGWLJm+ZGtHkzidFXoMVv1MAJ62vVrv5K/rfxUzjj7PO2Q2iBr+nfGSefZ90jtS3+R2//lX+TGu++WyD+cKPt/Mkv+PuJYNsgiZmvcVLMaA9Zs8zOzOWs1ncqr/gaj6r/yulmWzrBwOCyPPzxbps+0tjs8MVstYU1KRMzWKBJitsZPjVzKVjbIyu5cVbmI2Zo9O69zErPVGXFitjrjZyU3MVutUPJHGmK2+mMc8tkKYrbmk37p1O1lzNZMG2TpsVB37tor6WSrCkGQ6fF+N2SrisOarh6zUAbIVh+eO3Zk6/ev/LFUVFRI+cYN8tifn5Ib7pst1Xv3aL3aduLJ8tTEM+WiGT8VKStr7WkquZLqC3CqnZrNRNu2LZvkg/cWyJTzpyaRTbfpVDrhs3L5UlE33SeccnpSmUvZIEtjYlWiWk1nVaJaTYdsTb7Q/O/cB2XKBZdJhw4ds74KZdqEzUrByFZkq9k8QbZaOXu8T4Ns9Z55tjUiW7MlF82HbHXGz0puZKsVSv5Ig2z1xzjksxXI1nzSL526vZStiqq+unXS+OO1uKb6S4nLl978QIvnakW26vslGfdJUoJV37CqsanJlZitej3G1au6FH7+lXdl6/adrZt46X0bNewo7Zh6qTAErGzN8/mUjWxVTZ47Z7Z8/7yp0vnx/5H2jzwoeyIt8sRll8m1L74sB/7pejl48VSJVFYKslXk3QWvSdfuPWTIsJGmo51ODOsZFO+p066Vyqoq0zIyibAN61bLyhXLZOLk81POOKty1O10ViWq1XSqg6xsjR9mZKu9C+0Tj94vF146XWpq29nLaDE1K1vjQSFbLU4cj5MhWz0G7qA6ZKsDeMhWZ/As5ka2WgTlg2TIVh8MQp6bgGzN8wCUSPVey1aFVV/haoy5qsda7VLXMW6jK/W7Lmn/9fZH5L9uuap1cy2zeKm6fHVrZauqWxeu+pTQ21pTXR0XO1bFah0xZIB8unI1stVv50+mmK2Z2lvW3CTtnnpS2v/2v6VizVda8nDPXnJg5o/lwLQrJdK+Q6YieB8CBUuAmK0FO3Q03CcEcilbfdJFmgGBnBIgZmtO8VJ4CRAgZmsJDDJdzCkBYrbmFG/RFu5VzNaiBZjjjrFBlguAncpWYxNqn/+TtH/wN1K15EPtcKRDB9l3y3/Ioe+eL+FevV1oLUVAwF8EkK3+Gg9aU3gEkK2FN2a02F8EkK3+Gg9aU3gEkK2FN2a02F8EkK3+Go9CaQ2y1d8jhWx1YXzclK16c6o+eE+TrrV/fb61hQcvuFgO/vAaaR5zvAutpggI+IMAstUf40ArCpcAsrVwx46W+4MAstUf40ArCpcAsrVwx46W+4MAstUf41BorUC2+nvEkK0Ox8dJzFZjDNG9e3bLy395Vi7+/tVxLVr77luy6Z235KJ77m493nzcaHnpkkuk/ZjvyNARo+LSs0GW+YASs1WEmK3rZemSRXL2uRfbPuuJ2WoPGTFb7fEypv7y8xWy6esNMm7CJMuF5FK2ptqM0XLjXE6Y6jPO5WpcKY6Yra5g9KQQYrY6w8wGWc74WclNzFYrlPyRhpit/hiHfLaCmK35pF86decjZmvp0HWnp8hWhxxzLVv1DbImjBgt7X99r7R//FEpa2qUlydNks4VlTJswNFy8PLpEj68n9YTZCuyNdWURrYiW83mxoLXX5Q+/RrkqKOHOrwatmVHtmaPEtmanh2yNfu55STn5q+zv346qdervMhWZ6SRrc74WcmNbLVCyR9pkK3+GId8tgLZmk/6pVM3stX/Y41sdThGXsnWU8+YrLU0sGuXtH/oN/LO6s+l++bN8u0PPxQJBKRx3Gly6IofytMH98rYk8ZLz171cT0zE23btmySD95bIFPOn5pE4esNa2X5siVy1pQLk97TBbDeJmOClcuXirrpPuGU05PyLV3yvoRCIRkz9qSk9xa/v1Aqq6pk5OixSe+9u+A16dq9hwwZNtJ0tNK1Vc/AylZWtjqRBaxstXehRLba42VMjWxFtmY/e3KX08n1M3etcq9kZKszlshWZ/ys5Ea2WqHkjzTIVn+MQz5bgWzNJ/3SqRvZ6v+xRra6MEa5iNmaqVlljYek9s/zpN3vH5OqD99vTR7u2UsOXnq5HPzBlRLu0zdTMbwPgbwTIGZr3oeABhQ4gVyGEShwNDQfApYIELPVEiYSQSAlAWK2Mjkg4IwAMVud8SvV3MRs9ffII1tdGJ98yFZjsyvWrpZ2v/sfaff0kxLYuSP6Vmy1q5KujWdOEikvd6GnFAEB9wkgW91nSomlRQDZWlrjTW/dJ4BsdZ8pJZYWAWRraY03vXWfALLVfaalUCKy1d+jjGx1YXzyLVv1LpQFg1Lz8gvaatfqt98SiUS0t8K9ekdXu17xQ1a7ujDeFOEuAWSruzwprfQIIFtLb8zpsbsEkK3u8qS00iOAbC29MafH7hJAtrrLs1RKQ7b6e6SRrQ7Hx+uYrXpzU+0QrW8eUh9ukXaPPyrtnnpCyrdsltt+/nO59Re/kMZTJ2ghBhrPOEu2bd9KzNYY0Lfnvyy96/vKoMHDTWfEhnWrZeWKZTJx8vkpZ8yjD9wt0665QQKBQNpZ5XY6qxtfWU2nGh+JRES188rrZlk6Q8LhsDz+8GyZPvNGS+lDwaA88dj9Mm3G9WJXtr752gvSv2GgDBg02FJdeiInMQeJ2WoLtRCz1R4vY2pitqZnxwZZ2c8tJzmdXD+d1OtVXmK2OiNNzFZn/KzkJmarFUr+SEPMVn+MQz5bQczWfNIvnbqJ2er/sUa2Ohwjv8rW1g2ywmGpefN1+c3fPpVbf/lLkVBI67Fa7frV5dNkweF95Jx/nJZEgQ2y4pEgW9OfKMjWzBeSTEI/cwki6Tans5LfLM2C11+UPv0a5Kijh2ZbRFI+ZGv2KJGtyNbsZ0/uciJbz8wd3CIoGdma+0FEtuaesVs1IFvdIlm45SBbC3fsCqnlyNa20dq1Z5/MvPk+6VffQ26bNV1qa6p8MZTIVofD4HvZGuufOhlnXDJN2j/5uNQ+OVcq1q2Rjf36yesTJsilGzZpIQbUalc9tiuyFdnKyta2OcDKVnsXSmSrPV7G1MhWZGv2syd3OZGtyNZ0swvZmrtzTy8Z2Zp7xm7VgGx1i2ThloNsLdyxK6SWI1vbRmvxss/luRcWyN79B2XWtZfIgP71vhhKZKsLw+CXmK2WuxKJSPW7C6XdE7+Tmr8+L2XNTVrWcO/6aGzXH1yp/ZsXBLwgYDeMgBdtog4IFBIBYrYW0mjRVj8SIGarH0eFNhUSAWK2FtJo0VY/EiBmqx9Hxf9tImZrdIzunfOsnHT8CHnng0+l4fBect6kk30xeMhWF4ah4GSroc+BPXuk3TNPSu0Tv5PKVSuj7wQC0jj+9Ghs19Mntq52dQEVRUAgiQCylUkBAWcEkK3O+JEbAshW5gAEnBFAtjrjR24IIFuZA9kQ8Fy27tol8skn2TTVWZ4uXUSOPda0DBVC4PZf/0Fu+fFU+WrtJm2Fq19CCSBbnQ27lruQZaux+1VLPtRWu9b++U9SdmC/9pa22nXqFVqYAVa7ujBZKALZyhyAgMsEkK0uA6W4kiOAbC25IafDLhNAtroMlOJKjgCyteSG3JUOey5b588XmTDBlbbbKuS000TeeMM0iwohoFa03jDjItFjt9444yIZM/IYW1XkIjGy1SHVQorZetWPborr7bYtm+SD9xbIlPOnxh1XovXvz/5BPln7pXz/gQei76nVrqedoUnXFQ39ZeOGdXLqGZOT6K1cvlRU7K4TTjk96b2lS96XUCgkY8aelPTe4vcXSmVVlYwcPTbpvXcXvCZdu/eQIcNGmo5Wuviyeoa5c2bL1GnXanWYvTJtXsQGWelPFDbIynwhyTTHMpfABlk6o+bmJnlq7kNyxdU/sYItqzT5irlGzNb0w/X8H5+UsSeOl9ZNILMaXW8yzXvmcTll/ETpdlhPbyrMYS3EbCVma7rpRczWHJ58saKJ2Zp7xm7VkK/7B7faTznOCRCz1TlDSshMIC8xWz/+WOTGGzM3zu0Uo0aJ3HNPUqmHGpvl1rsekwvPGdcqV1VIAfVS8jXfL2SrwxEoRtmqkOgC85yjh0v73z0s7Z75g5Tt26fR+vSEE2XVxEky7uxzJdyvIY4gsvUGCQQCaWeV2nhq2jXupVMX2kSRbtYAq+lU3kgkImyQ1UaRDbLsXSjZIMseL2NqZGt6dsjW7OeWk5zIVmRruvmDbHVydlnLi2y1xskPqZCtfhiF/LYB2Zpf/qVSe15kq8/grl6/WWbcdI9s2bYjrmXDBx8pD95xvXSp65jXFiNbHeIvdtl61pQLWwm1++MzUvv7x+SL/Xvky4ED5bx586R5zPFy6PyL5dAFl0hL586CbHVPolqVslYlqtV0yNbkiwKy1d6FEtlqjxey1TovZKt1Vm6mRLYiW5Gtbp5R9stCttpnlq8cyNZ8kfdPvchW/4xFMbcE2Soy76WFsuijlXExWs1Wu+ZrHiBbXSBfLDFbraKoWLdGap/+g9T+8RlR/9ZfajOtQ5d+X4IjRkqo/xFWiyNdiRNgg6wSnwB03zEBYrY6RkgBJU6AmK0lPgHovmMCxGx1jJACSpwAMVtLfAJk2X3PY7Zm2c5cZNOl6tjRQ+S8SSfHVaEk7LqNW/MeSgDZ6sLIl5psNSKrWrxIap97WmrnPSuB3btb3wr3PVwOffc8aZxynjSPHuMCZYooVgLI1mIdWfrlFQFkq1ekqadYCSBbi3Vk6ZdXBJCtXpGmnmIlgGwt1pHNbb9KWbbmlqw7pSNbXeBYyrLViK/mtZel9pn/lZpXX5KyxkNt4vXwfnLonO9J43cRry5Mt6IrAtladENKhzwmgGz1GDjVFR0BZGvRDSkd8pgAstVj4FRXdASQrUU3pJ50CNnqCeasK0G2Zo0umrGUYrbqqL76YqVsXLdGTj1jchI9FbN199+3yWkHDknN//1Jal5/RcoOHdTSvXPyydLctZsc37e/HDr3AgmOHNWaf/H7C6WyqkpGjh6bVOa7C16Trt17yJBhI01HS9/MyxhfNjHh3DmzZeq0a7U6zF6ZdorfsG61rFyxTCZOPj/ljLEaY9XtdFZjsVpNpzrIBlnxw0zMVnsXSmK22uNlTM0GWenZEbM1+7nlJCcxW4nZmm7+sEGWk7PLWl5itlrj5IdUxGz1wyjktw3EbM0v/1KpnZit/h9pZKvDMUK2xgNM3CBLiVYlXJV4XXxgr4QkIuPffFPLFO7XIIe++z1NvL53cB+yNWEuWpWyViWq1XTI1uSLArLV3oUS2WqPF7LVOi9kq3VWbqZEtiJbka1unlH2y0K22meWrxzI1nyR90+9yFb/jEUxtwTZ6v/RRbY6HCNka3rZanx32fsLRb78XE59a4HUvP5qXKiB1887XwIDj5IRE6doG2wZX6xsDaSdpVYlqtV0yFZkq8PLoiBbsyfIytb07JCt2c8tJzmRrchWZKuTM8h5XmSrc4ZelYBs9Yq0f+tBtvp3bIqpZchW/48mstWFMSJmq32I2orXV1+SmufnJYnXUMOR2orXRhVqYPix9gsnR0ERIGZrQQ0XjfUhAWK2+nBQaFJBESBma0ENF431IQFitvpwUGhSQREgZmtBDZdvGkvMVt8MhWlD8ipbd+3ZJzNvvk+Wr1qT1Ljhg4+UB++4XrrUdfQ3QRFBtjobIk28vvJiVLy+8VrcitfQEQPaxOuwEc4qIrcvCSBbfTksNKqACCBbC2iwaKovCSBbfTksNKqACCBbC2iwaKovCSBbfTksvm8UstXfQ5RX2XrvnGc1OjfMuMjflDK0Dtnq3vCVHTwQL16bGlsLDx05UA5N+Z40fu8CCQ4d7l6llJRXAsjWvOKn8iIggGwtgkGkC3klgGzNK34qLwICyNYiGES6kFcCyNa84i/YypGt/h66vMlWtar1ll89IrOuvUQG9K/3N6U0rSNmazycxA2yjO8uXfK+hEIhGTP2pCSii99fmLRBli5e31m2WHp//JGMWbSoTbwOOKpVvK7t0F6WL1siZ025MOVIzZ0zW6ZOu1arw+z19vyXpXd9Xxk02Fzibli3WlauWCYTJ5+fsg6rG1q5nc5qLFar6VQHI5GIqHZeed0sS+dmOByWxx+eLdNn3mgpfSgYlCceu1+mzbhe7MrWN197Qfo3DJQBgwZbqktP5CTmIBtk2UJNzFZ7uOJSE7M1PTxitjqYXA6yOrl+OqjWs6ypZOuqzz6RHdu3yomnErM13WDs2vmNzH/lL3LBpdM9G7NSq4iYrYUz4sRsLZyxylVLidmaK7KUayRAzFb/zwdkq8MxQrbmTrbqJasNsrp1qpPRX3wpNS/8WYv1anx98Q8nyHsTz5KzTxyftLmWng7ZKoJsXS9LlyySs8+92PZZj2y1h4wNsuzxMqZGtiJbs589ucuJbEW2Iltzd35ZKRnZaoWSP9IgW/0xDvlsBbI1n/RLp25kq//HOm+yVaFRYQQaDu8l50062f+kUrQQ2eqNbO3avYcMGTZSq6zswH6pfemvUvPnP2ridfXAgfL+2LFy2ZNPSkv3w+TQxLOl6cxJ0nTqaRKpqdXyIFuRrU5kAbLV3iUa2WqPF7LVOi9Wtlpn5WZKJ9dPN9uRq7JY2eqMLCtbnfGzkhvZaoWSP9IgW/0xDvlsBbI1n/RLp25kq4jZHlB+2vspr7J19frN8od5b8ismZdIbY35492FcLoQszV/o6SJ1xf+LNVqg6235msi1vhqnHCmJl4bJ50j4Z698tdQak5JwG4YAVBCAALxBIjZyoyAgDMCxGx1xo/cECBmK3MAAs4IELPVGb9SzV3qMVt12XrjjItkzMhjtGngp32h8iZbzSy08STxk5HOdPIiWzMR8u796rfmS/Xrr0jtK3+V8g3rWyuOtO8gzceNkuDIUdI0boI0jRvvXaOoKS0BZCsTBALOCCBbnfEjNwSQrcwBCDgjgGx1xo/cEEC2MgeyIYBs3Sczb75PjLJ13ksLZdFHK+W2WdPzvqAzb7I1m8mUyzyLl30uP/jpHVoViaJXDdi/3/mY9t7Zp41NGjhkay5HJvuyKz9fpYUZqH71Jala8oFIS0ubfK2ukebv/IM0nXq6Jl6DQ803xsq+dnJaJYBstUqKdBAwJ4BsZWZAwBkBZKszfuSGALKVOQABZwSQrc74lWpur2XrrsZd8snWTzzH3aW2ixzb89ikelOtbPVLqFJkq4go0XrPnGflwTuuly51HeMGMfG9xGXJxGyNn/Mrly8VFbvrhFNOTzoZli55X0KhkIwZe1LSe4vfXyiVVVUycvTYpPfUBlnGmK2JCb7esFaWL1siZ025MOWJP3fOffLDXv2kw4I3pebtN+NWvapMf774Eqnv2EkGHfdtaTztdC32q/G1Yd1qWblimUycfH7KOh594G6Zds0NEggE0l6A3E5ndeMrq+lU4yORiKh2XnndLEsX03A4LI8/PFumz7zRUvpQMChPPHa/TJtxvdiVrW++9oL0bxgoAwYNtlSXnshJzEFittpCLcRstcfLmJoNstKzI2Zr9nPLSU4n108n9XqVl5itzkgTs9UZPyu5idlqhZI/0hCz1R/jkM9WELM1n/RLp+58xGydv3a+TPj9BM8hn3bEafLG5W8k1Zvqafn/vGm6L/aFyrtsNa4o1enNnX1za8yFXI+kGqBbfvWIzLr2EhnQvz6pusRNvBLlK7I1Hpl/ZetsmTrtWk3oqlfFmq+kesGbUv3WG1L97tvyl9NOk/7r1snIZcu094ODh2jhBppPPU2a/uFEWb91M7I1zcmIbM18pXp7/svSu76vDBqc/Srqr75YKRvXrZFTz5icue/HcPsAACAASURBVEKLKRa8/qL06dcgRx091GKOzMmQrZkZpUqBbEW2Zj97cpcT2Xpm7uAWQcnI1twPIrI194zdqgHZ6hbJwi0H2Vq4Y1dILc+HbP14y8dy42vWFle5yXJU71Fyzxn3pJStxjACZqtd3WyLnbLyKlvNVpSqTbNm3HSPXHvFdz2x0Xp9W7btaOX2w3+cJDfMuEgONTbLrXc9JmNHD2lti0r/r7c/Iv91y1WanEW2FqZsjWt1OCzvPvek9Nu0SUa98bpUfbxEJBxuTRKpqpZVZ0+WJSNGyFnjzpDgsBEiZWVJ55nbK1atlmd1xarVdKpjrGyNH15Wttr5WBFWttrDFZca2YpsdTB9cpYV2YpsTTe5kK05O/VaC0a25p6xWzUgW90iWbjlIFsLd+wKqeX5kK1+45NKrCYumMxXu/MmW3WReeE545JWsSoJ+9wLCzwJaptYlz5gF50zTs4aP1aTrcY2JspWNXAHG0NysKlNzuVrMKnXHQJl+/dJ1YK3pOLN16XqrflSvnZNXMEt3btL8JTxEhw/QZpPP1NaevR0p+ISLaVrpyrZtT8okZZIiRLIXbcV0eQ/C+SuPkrOD4HOHapk/6GQhMJtcanz0xJqhUBhEujYvlKam1ukKci9XGGOIK3ONwEVb7IlEpFDRfZ9iPuofM+s0qm/prpcKgJl2v0cLwhYJaD2bSjll5lsZWWriKR7fF8JzbseeFpu/9lVSTFU3Z5MZmJX38Hs5n+6TO74zZNpV7aq9oRbItp/vIqTQNnGDRJ45RUpe+01CSx4S2T37riORgYPkZYJEyRy+hnScsopIu3aFSeIHPWqsjygSSLOIPcBq+tSeQDd6j5Zf5VYUV6mfQZFOIn8NTC0JiMBv1yj1BdcJYq4lcs4ZCSAgCkB/V6j2L4PhcMRKS/nPoppn3sC5eqpybKoV+AFAasEqirS7xVjtZxCTZcqZquXYUnTsSv5la1mYlfJ1nUbt2qhBDLFbFVw9x0Myj7+ClWo56i9dre0aGEGqt9+U6rffEOqPvpQJNT2F8hIVZU0jxkrzadOkMZx4yU4YqRIhg2z7DWg+FLb3SCr+AjQIwg4I6D+qr3nQFCCIVa2OiNJ7lIlkGqDrFLlQb8hYJdAp3aV2h8sWJVnlxzpIRAloFaHqz+eq/s5XhCwSqC+W63VpKTLA4G8yVbVVyU1n31hgTx4x/WtK1i9jtmqhzPo1aOrJlcTlx0nxpVV8lW9VFr1ImZr/KwtlA2yEs+1TJsXbVi32nSDrLID+6X6nbelesF8ube+h/zrL38pgZY24dHStas0nXSqNI2fIE2nnSHhXr3FaixWq+msxmK1mk6xIWZr/AwhZqu9Tyc2yLLHy5iamK3p2T3/xydl7InjpWev5A0ts6eem5zznnlcThk/UbodVvihZojZSszWdGcJMVtzcw0xlkrM1twzdqsGYra6RbJwyyFma+GOXSG1nJit/h+tvMpWhUfJzB/89I44Ul4v+01cfvyfN02P25xLSeF/v/MxrY1nnzY2LpYssrW0Zaux90qOXjnlYmn3xqtS9fZbUv32fAns3RsHKDTgKLn9+1PlmuPGSnj4SGnp1CnlVQLZao7mzddekP4NA2XAoMG2rrBOZAGy1RZqNsiyhysuNbIV2epg+uQsq5PrZ84a5WLBqVa2rvrsE9mxfauceCqyFdnq4oTLoihkaxbQ8pQF2Zon8D6qFtnqo8Eo4qYgW/0/uHmXrf5HlL6FyFZkq07ATI5WLf1ICzlQ9eYbUv3eO1rS//yP/2hdARscMkyaTj5Vgt85QZpOGhcnX5GtyNZEAl99sVI2rlsjp54x2bVL54LXX5Q+/RrkqKOHulYmK1uzR4lsRbZmP3tylxPZimxFtubu/LJSMrLVCiV/pEG2+mMc8tkKZGs+6ZdO3chW/481stWFMSJmqwsQS6CIskMHpfrdhVKlwg4sfEsqV3ya1Ot08rVYERGztVhHln55RYCYrV6Rpp5iJUDM1mIdWfrlFQFitnpFmnqKlQAxW4t1ZHPbL2K25pav09KRrU4JskGWCwRLs4iyffuk+v8tlKp3F2o/Kz9bLmKI9yplZRI6apA0j/62BMccL83f+rYEjxlSdBtuIVtLc/7Ta/cIIFvdY0lJpUkA2Vqa406v3SOAbHWPJSWVJgFka2mOu9NeI1udEsxtfs9lqx4fddrFE+V3z7wiy1etMe3h8MFHxm2clVsMzkpnZaszfuSOEsgoX9XGVe07SPOo0XECtqVb94JGiGwt6OGj8T4ggGz1wSDQhIImgGwt6OGj8T4ggGz1wSDQhIImgGwt6OHLW+ORrXlDb6liz2Wr3iolXW/51SMy69pLZED/+F2F1aZZz72wIG4jKku9yUMiYrbGQ1+5fKmoXWlPOOX0pNFYuuR9CYVCMmbsSUnvLX5/oVRWVcnI0WOT3nt3wWvStXsPGTJspOkIf71hrSxftkTOmnJhyhkwd85smTrtWq0Os9fb81+W3vV9ZdDg4abvb1i3WlauWCYTJ5+fsg6rMVbtpLtmyCip+XixVC7+QKo+WiyBnTuS6r/t5z+Xm5avlOYxx0tQrYIdNkIilZVJ6VRcl6t+dJOlsyQSiYi24dd1syylD4fD8vjDs2X6zBstpQ8Fg/LEY/fLtBnXi13ZygZZxGxtbm6Sp+Y+JFdc/RNL8y2bRPmKuUbM1vSj9fwfn5SxJ46Xnr3i7xuyGeNc55n3zONyyviJ0u2wnrmuKuflE7OVmK3pJpm675v/yl/kgkun53wulmoFxGwtnJHP1/1D4RAq/pYSs7X4x9gPPSRmqx9GIX0bfClbV6/fLHc98LTc/rOrpEtdR19TRLbGDw+y9QYJBAJp56wd2TrtmvjyKtaulqoli6VyyQdSteRDLfTAbf/2b3Lrz3/eWmekukaCxx4XDTsw5tvaKthwfR9Btq6XpUsWydnnXmz7mvK/cx+UKRdcJh06ZH89yiT0rTSKDbKilJCt8bMllytb31v4htR17ipDR4yyMkVzngbZmnPEphUgW5GtyNb8nHt6rcjW/PK3Uzuy1Q6t4kyLbC3OcfVbr5CtfhuR5Pb4UrbOe2mhLPpoJStbRSSVXEn1BTjVF1Ez0bZtyyb54L0FMuX8qUkzI91q0XTCB9maW9maOFBlTY3y8CO/luvLqlsFbPmWzUnjGe5dL7+ccbWWrvlbYyQ4cpQoKZvqxcrWeDLIVnsfZk88er9ceOl0qaltZy+jxdTIVmQrK1stniwuJUO2IluRrS6dTFkWg2zNElwesiFb8wDdZ1UiW302IEXaHGSr/wfWc9mqVq3OuOke2bIt+XFoHVfvnt1kzp03JoUX8CtOYrb6dWRKs13lmzdJ1UcfSuXiD6OrXz9ZKkrKGl8qzEBo6PDo6tdvqc23xkio4ci8ALMbRiAvjaRSCPiYQC5Xtvq42zQNAq4RIGaraygpqEQJELO1RAeebrtGgJitrqEsqYKI2erv4fZctuo40sVs9Tey5NYhWwttxEqrvWXBoFR8tjwqXrXwA4ulYl3yxnRqo63m0WMkOOZ4LfRA86hvSaRDh5zDQrbmHDEVFDkBZGuRDzDdyzkBZGvOEVNBkRNAthb5ANO9nBNAtuYccVFWgGz197DmTbb6G4u91iFb7fEidf4JBHbulKolH0Q33lICdulHUrZ/f3zDAgEJHTVIk67B476l/VSrYc0233LSI2SrE3rkhYAIspVZAAFnBJCtzviRGwLIVutzoCncKBH1v4j2/7Gfov2UpOOxNIbjEpcvsZy29ObpovXo9bbWV2ZSjlk6Y/uipbSW11afVmpC/2K/l0VUA1K/H8sXxyGhnrb2K+ap+6+ns1JfJKFdsZINfWjrk6W2tbQYGEfnhuq19lMb54RjEpHK8oCUBUSaguGEdIb0sbxtZbXNu7Tlm9SpuGjtUX3X25NQvrG92ZZv1ncndSqS0XbZ67ve/rj2xMqKH5c042Q2dmmOGeu0wtt0bsSVH5s3sWOXDLlcfnLCldYvPqT0nACy1SFyNsiKB0jMVm9jtir6Vje+SpuupUUq/7Yquvp18QdS+fFiuf2iC+U/fvGLuAGOVFVLcNhwCR43WoKjxkjzcaM1IRtuaZHHH54t02feaOmMCgWD8sRj98u0GdeLXdn65msvSP+GgTJg0GBLdemJnMQcJGarLdRCzFZ7vIypv/x8hWz6eoOMmzDJciG5lK1skGV5GJISznvmcTll/ETpdljP7AvxSU4n10+fdCFtM1LJ1lWffSI7tm+VE08lZms6gLt2fiPzX/mLXHDp9EIY7oJpY0ukRVoiYWmRFjl4cL+88MxTcu73r5CIdrxFwhLWfrb+3hJNq34Pq3yRSOv70d+j+bQy9X+rHCpfLK2WLlaGVof2XqzMWH3R3yPR99T/Ii1SVSna74eaQ0l1hA1tVOm1fLFyk+uLlqf1wVCf1l5DfdE+RPuntTMSbuMQ618k1rZoH5SUi3GJHdd+j7VDK6e1vli5hvp0fq2s9Tq0MmNcY9x6tPSQyTJZHpaHC2au0VB3CdwsN8tsmS2NEh/Gzd1aKK3UCfxcfi633nprqWPwdf/zKlvTxW8dPvhIefCO66VLXfa7f3tBHtmKbNUJPPrA3TLtmgKVrQkni75B1nUjvi2Vyz6WqsWLtNiv5Zu+TjqtVKiBQ6O+JfeOHyczj/uOtHTuLMGhw9OefsjW7K5O6Tany65EkQWvvyh9+jXIUUcPzbaIpHzI1uxRIlvTs0u1CWT2xHOXE9maO7Zul4xsjSfaGD4kwXBIwi0hCUVCEmoJSqglJGHt32EJhYOx49H39u/ZK+sXfS5HnjZck1e65EolwczEly7f2iRgmxBUf9DVJJcuz5RYa4mKs3gJFmn7XZd6MRGWJO9adPGol2usLybuWuuLykUl+IzSsk1qxspoFXcxWWeQhRoXg5Rrk4NRaamYp3t1kA5yjVwjd8vdbk9/ynOZQL3UW5Kt1eU1Uqb+V6b9f+ynaD8l6XgsjeG4xOVLLCf2e6RMFZiyHr3e1voi8eUoNIlttNK2srJAXH+S+pmyXdH6jfXG1RdtTQIz9Xssh1lfDX1yxixaTxIzvU2GdtUtq5F9w5tFKkTKAoHW/uhTzdhH/Vh1ZbkEytTK1ujKSlVPa/rYv/V80eMxTjFexjz6obgyDFxTtsNYlhojVYtp+anblqnO2GC19s+0/Ax91+toy2toT+yf6dphv05j+W3zM914Jo6xWZ06C41zAu9MHFV5n/3fImSry9dvt4vLm2w91Ngst971mIwdPUSOHTpQ/jDvDZk18xKpramSe+c8KycdP0LGjDzG7f66Xh6yNR4pK1uLS7Zeed2suAHWwg+ozbc+WiyVHy+RqmUfiToWrqiQ22++Wf7tl79sTR8cfqwEjxksoaEjJDRkqAQHD5Vw73rtfWRrdpciZGuUW3Nzkzw19yG54uqfZAfSQq587SaMbE0/OMhWC5M3B0lKfWXr8SefGpOL4ah4jERFZLAlGBWSSkBGYkJS+10XlCq9EpZt70XzqPRtErO1PCU4teMhCYab2+oxHI+rPxyVodE26PWEE9rWdrytXUaRGs2r15vN9OkhPeQCuUAekAeyyU6eFAQCZQHR/pNy6VjWUaaFp8mcqodEiazo8YCUl5W3/l4eULImIGXqnTL1XvSnSq/StZanHY/9rtLG8gXKyqL5xKRMQ33RMstEr0+1o7qyQiszHJZYPcb6VLlt9Wn5WutMrK+tzXr/Wtut8sT6Ft+XeA7Rvpe3cdDqUqIsxkVinFQ6QzuiZUfTauUn1BfHOtafaJllUb4SZb3rm2/k3bdele9dfAVzu0QJ/P6RX8vFl18t1dU1lgkQs9UyKhLGCKinVlnZ6u/pkDfZatwgSyG664Gn5fafXaWtZF287HN57oUFctus6Zp89fuLmK1+HyHal0sCFevXSuXSj7SVr5UrlkvFqs+kfOsW0ypbunSR4OBhokRs+JjBEhwyTLqc8C35e3NAWloMwXdy2WDKhkCREchlGIEiQ0V3CpTAvua90hhulKZQozSFm0TFPFT/bg43SWMo9rt2rO3f+vHG0CFpDje35lF5G1UZWlnRPCFpkkPBRmkMxpd1MHSgQIk5b3a7ivZSUV4hFWUVUh6okMpAZdu/yyukvKxCKgKx98oqtTStsi5gEF1xUsoo/eIFXKvcMgixVqmoSbGoyIoTcDFRpovEqFSMiciYdIyTaUbBqMm0trRGwajJOtWHxPpiwi5eXBqkYoK4bJVzuvzUy4wJTV2EarIuUC415bXOBy5PJRCzNU/gqbZoCCBbi2YoPe0IG2R5itt2Zb6QrV07d5Tbf/0HueXHUzXZqsILGOWr7V55nAHZ6jFwqvM9gcCePVK54lOpWLVCKlYsl8pVn0nlqpVSdtDki6vaiKv/ERIcOkxCg4dqm3AFhwyV0BEDtEefeEEAAukJIFuZIbkkoB5v1iVmVHi2yc428alLyniJqUSnMU9zi3o/mrYxQZxq4lOTp1GJqqeL22Ailx1NUbZ6lK9VNgbKo/+OyUUlGyuUhIwdV2Iy+ntUUmqyUgnJsqisVHnNjutCU5XTmr+80lBP2/H4+stj5UXrbKsn/ni0XW1t0f5dbpCnsbzqOC8IZEMA2ZoNNfJAoI0AspXZkA0BZGs21LzLkzfZagwjcN6kk7XQAQ2H9xL173kvLZRFH61kZat384CaIJB7ApGIaKtgP1uhrX6t+Gy5VK78TCrWrhZpaUmqP1LbLhqGYMgwCQ0dJsEhwyU4bIQWE5YXBCDQRgDZWryzQT0yrlZhRuVj20pMTWAmrOTU5KQmM2PvKXEZPBRNF5OYeh6tvFh+sxWj+upPVX++X5XlVVJTXi1VgWqprqgRFeewujz67xr174rq2LEaqaqo0v5do6WLHq+prG39t56nqjyaTvuvoloO69hRyloqRSKVsWN6/mpNfvKCAATSE0C2MkMg4IwAstUZv1LNjWz198jnTbYmYlFhBWbefJ8sX7VGevfsJnPuvFEG9I/Gd/Tzi5it8aNDzNbijtma6lwMh8Py+MOzZfrMGy2drnExW2tFdn+wVMpVCILPVkjlyqiEDez4xrSsP069TAYEg3LUYb0kpASsWhE76BiJVKUPOeIk5uD/zn1QplxwmXTokP2GfW/Pf1l61/eVQYPTbx6WDiAxW6N0iNkaP0tyKVvfW/iG1HXuKkNHjLJ0buc6kR9itioReSB4QA7G/jsQ3K/9rn6qY9Hf94ssbZTd/Q7KwSpdeLatCtVXdaaTqJk2zMk1a718JSSPChwlx4ePl1fbvxYnLjWZGROfurhUP6vKq6Wmolp7LFoTpLoQjUnQqCStiQnUNlmq8hglaG1FO0+6yQZZzjDv2vmNzH/lL3LBpdOdFUTulAQOHTwgf3p6rlw2/TpfUkK2tg1LvmK++3JilGijiNlaogPvcbeJ2eox8Cyq841szaLtvsiCbEW26gQefeBumXYNstXKiWllg6zy7dukYuVnmnytWLkiuiL2b5/LvCnnyNGffy7DVqxoq6qiQkIDBmqrX6OrYIdFN+Q6vF9rGmSr+cgseP1F6dOvQY46eqiVobOU5olH75cLL50uNbW5ESXIVmRrz16Z/xirHj9XEvRQqwBNFKIH5EBzVJDub94XFaihqETVBerBmDg1ylS1mZCVl9o5/M/yZ9kqW60kN02jHus2ruRsFZpqladxVWeCxGzNU9m2elMTnAkrQZUEVceVHNVXguqrR5VEVatK1cvJ9TPrznuYEdnqDDay1Rk/K7mRrVYo+SMNstUf45DPViBb80m/dOpGtvp/rPMmW40bZBXCCtZUQ4lsRbYiW0UcrWztUiN/39NkeYOsBc89KQOCIRmyZm00FuzKFVK+cYPpKRpp30ELPRAcPETWHn2MfFBdLmeffb60dOtu6+rMylZbuATZao+XMfWXn6+QTV9vkHETJlkupNhWtu5s3CH7mvfI3qa9ojZG2tu0R/Y279GOhT4+ILvq98qeKiVHoytKo0K0TZzqx3O5MrRTdZ20r+wg7Svbaz/bxf07eqz3l10lNDAgVZ1iMlN/7N1kJaj+OHzbo/LVWhl+eSFbz/TLUPiyHcjW3A8LsjX3jN2qAdnqFsnCLQfZWrhjV0gtR7b6f7SQrS6MERtkuQCRIkqWQE+bstUMVNmB/Zp01WLA6rFgV60QtVGX2aulq5IgR0tw0NESHnS0hI5S/w2SUL8GkUCgZMeCjhcmgVzKVjtE1ErSfU17ZU/zbk2S7tMkqfq5V5Ol2n+Nu5OOGdO5ufu72thIPYauhKhRhmq/VxlFaTtpX9XRIE6jslSXqca8Hao6itphnFdxEUi1srW4eklvIJA7AoQRyB1bSi4NAsRsLY1xdruXxGx1m6i75eVNtqpuqE2xTjp+hIwZeYy7vfK4NGSrx8CprqgIuCFbUwEp/3qjVH6+Uiq1eLDLpWL1l1Lx5RdSduhgSoYq/EBo4FESGjhIwoOO0YSsErGRdu2LijudKR4CbslWfRXpvuBe2dO0W5Okrcea1bFdcatNtfead7ce0+KUuvBSq0Y7VdVJx6pO0Z/VnUQ/1rm6i7TT5Gl7aR8nTPUVpkquRleadvDRylAXsFBEDgkgW3MIl6JLggCytSSGmU7mkACyNYdwi7hoZKu/BzevsnX1+s3yh3lvyKyZl0htTfrNbfyMEdnq59GhbX4nkEvZmqrvgb9vl4p1a6RizRopX7taKtaulvK1a7RjgZ07k7K1dOokwWHHigTKJNyrt4QbjpRw/wbtP7UaNtz3cL9jpn1FTEDJ1s27d8nOg7r4jK0ibVKP3qvH8NV/0dWmmkQ1rDbVjylR2hJpcURJrSTtUNVBOlZFZakSpJowre4UPVZdJ3XVnbVjmkCNHdPFasfqOgSpoxEgc7YEkK3ZkiMfBKIEkK3MBAg4I4BsdcavVHMjW/098nmTrSpm68yb75Plq9aYEho++Eh58I7rpUtd9rt/e4GemK3xlFcuXyoqdtcJp5yehH/pkvclFArJmLEnJb23+P2FUllVJSNHj016790Fr0nX7j1kyLCRpkP69Ya1snzZEjlryoUph3zunNkyddq1Wh1mr0w7xW9Yt1pWrlgmEyefn7IONsiaLdNn3mjptLOyQVaqgt587QXp3zBQBgwabKkuPZHVmINl+/dL5Zd/k/L1a6Vi3VopX/2VzOnfR37wp3nS5csvUtYZOnKghPv3l1D/IySs/muIidgjBoiStZnmmJXOfPXFStm4bo2cesZkK8ktpWGDLHNMXsdcU9JzV+MO+WLVZ7Jz23apHVIXi1W6V3v0Xq021WKXagJ1X1Sexo4pgerGq11F+5gE7SRdarppq0eH7jlGIrVlEukdkE41naOC1LDaVF99qotTN9qRrozn//ikjD1xvFjZICvXbclU/rxnHpdTxk+Ubof1zJTU9+9bvX76viMpGsgGWc5GjpitzvhZyU3MViuU/JHG6/sHf/SaVhgJELOV+eAFAWK2ekHZWR15k63Omu2f3MjW+LFAtt4ggQwxP92WsupCe9WPbsp4UlhNpwqKRCKi2nnldbMylqsSeLlBVq5lq1mH9Q2yOlbXSPnG9VK+fp1UrF8n5evWSPn69VKxfq0mZ1PGiO3cWZ7/3nnSJyIypGOdJmOjUra/hA/vL5HKSkucka1RTM3NTfLU3Ifkiqt/YolbNomy/bKkVocq+ak2eVLydFfjTu2/nYd2yM7Gb2TXoZ2x99TxtvdDLUGtmSNlpDRIg7aLvZ1XTXlt9FH71lWknaRjZXQ1qTrWSa0oVT8Nx9RqU/1YXU1nUStTE1/vLXxD6jp3laEjRtlpTs7SIltzhjZtwchWNshKN0GQrbk/L5GtuWfsVg3Z3j+4VT/l5J8AsjX/Y1AKLUC2+n+U8yZb1crWW371iMy69hIZ0L8+jtTiZZ/Lcy8skNtmTfd9eAFka/wkR7YiW61c9vy8sjWdbO3QIf1K+8DevZqA1URsTMbq/35x5Ajpv2aNjFy2LL6KQEDCvetj8rVBwg3RlbEhJWKVkO3RtioO2eq9bP2HyWfIrqadsrtVmEYFqVGeKrm649A3mkzNdpWp2sipS01XOa7sOOkT7iMbe2+RTpXqcfs6USK0k/bofZ12rHNNVy1eqX7s6J59ZM+BoARDzsIAmM19ZKuVK5p5Gla2Zs/O65ysbHVGHNnqjJ+V3MhWK5T8kQbZ6o9xyGcrkK35pF86dSNb/T/WvpStKpbrXQ88Lbf/7CrfhxFQQ0zMVv9PdFroXwL5iNmaFxotLVK+ZXN0FeyG9VqM2PIN66RinRKza6V8+zbTZkVq20m4X7+ojFXxYY/QV8U2SKh/Axt3ZTGYuxt3yfaDW2X7wW2y/YDx55a4Y9mIU7U6VK0W7VLdTbrWdtMEqvqva0136dqum3Sp7ho73vZet5ruUlmefdxytzbIygIlWSBQFASI2VoUw0gn8kiAmK15hE/VRUGAmK1FMYyed4KYrZ4jt1WhL2XrvJcWyqKPVhbEylZkq635RmIIJBEoGdmaYezLGg9FQxJsiMnXtWuj/16nfq6XsgPmO723dOsuIX01bL/+Ej7iSE3KhhoaJFzfV6S8vCRmXTDcrInSvx/cJtsObo3+3B+Vp8Zj6neV1sqrvKxcOld3kS61UVnaKk41idpdO24mT1U+L1/IVi9pU1cxEkC2FuOo0icvCSBbvaRNXcVIANlajKOa+z4hW3PP2EkNnstWtWp1xk33yJZtO1K2u3fPbjLnzhuTwgs46Wgu87KyNZd0KbvYCSBbrY1wYMc30Q271GrY9eskoCSs2sRr/Xop37RRBc5NLqiiQkJ9Dm/drKul4ciohFUytv8R0tK1q7XK85hKrS6NW4F6cIu2GnXbgahQ1d47uFXUalUrr4pAhXSrPUwOq+0hh7XrKYe17yHda3tIj/Y9pXttApdPFAAAIABJREFUTzmsXQ/p0a6ndG/XQ5OrZnFMrdTjZRpkq5e0qasYCSBbi3FU6ZOXBJCtXtKmrmIkgGwtxlHNfZ+Qrbln7KQGz2Wr3th0MVuddMjrvMRsjSdOzFZitlo5B4s1Zmu6vr89/2XpXd9XBg0ebgWRaZp0MVsr1GZdGzaI+hlQQlatjNVCFqyTwM6dpuVF2neQeRddJA3NQTmmslpa+h4u4V71Eq7vE/2vb1+J1NTabu8Tj94vF146XWpq25nmVZtB6StOjSJ124Et0VWoukg9sE2aW5qSyqiRGvmp/FTukDu09yoCldJdCdR2MYFa20O6t2+Tpj3a9dKEqnpfCVQrr3zFXPvy8xWy6esNMm7CJCvN1NLkUrYSs9XyMCQlJGZr9uy8zknMVmfEidnqjJ+V3MRstULJH2nydf/gj97TCkWAmK3MAy8IELPVC8rO6sibbHXWbP/kRrbGjwWyFdlq5exEtlqhlJwm2w2yVAgCbaOu2GrY6OZd0TAFL4wYLgO++EJGfPqpaaNa6uqi4rV3VMC29Okb/dm7vu1YXV1r3jW7v5S3n/urdPpOL9kR+rsogdr6KL8WH9X6KtRe7etNBWq3qu6yc+EmOenCs20JVDvU8/VlCdmafpSe/+OTMvbE8dKzV/zGmnbG1qu0yFavSDuvB9nqjCGy1Rk/K7mRrVYo+SNNvu4f/NF7WoFsZQ54RQDZ6hXp7OtBtmbPTsuJbEW26gQefeBumXYNstXKKYVstULJPdmarrYFr78ofbt0k8HVtVK+eZP2X0D7+bW2oVf02GbZ23JAvu4ksqmjaD/1/zapf3cOyNedymRndTSUwU1yk9wv98tBORhXdVWgWrq3UytQo4/rq0f51eP62iP8+uP7sUf566o7p2x2c3OTPDX3Ibni6p9kB9JCrnx9WUK2IlstTE/Pk2z+er0sXbJIzj73Ys/r9qJCZKszyshWZ/ys5Ea2WqHkjzT5un/wR+9pBbKVOeAVAWSrV6SzryevslWFEph5832yfNWapB4MH3ykPHjH9dKlrmP2vfMoJzFbPQJNNUVJgJit+R/Wbw79Xbbs3yRbD2zWfm7Zt0k2H/hatu6P/b5/sxwMHUjZ0HbhgPTdVyZ9d4Wlz16RfntEehwQ6XlApNd+kZ7qvwMincPV0nLYYRLu0VNaDuuh/Rc+rIdEevbUfrb0UD9j73VOLVvzT8xfLchlGAF/9ZTWQCA3BIjZmhuulFo6BIjZWjpjTU9zQ4CYrbnhWuylErPV3yOcV9l675xnNTo3zLjI35QytA7ZWtDDR+PzTADZmvsBOBDcLxv2rpP1e9bKxn3Rnxv2rJMNe9fKxr0bpDF8KGUj1ArTXh3qpXf7PtK7Qx+p79hX+6n9165e+9mpOhpGoKypUcq/3ijlW7ZIYNtWKd++VQLbt0lg6xYp3xb9d7n6PUUM2cRGhHvXxwRsTMz26i2R7t1jYrZX9Kf6rwA2+srlKCNbc0mXskuBALK1FEaZPuaSALI1l3QpuxQIIFtLYZTd7yOy1X2mbpaYN9laLBtkqcFAtro5JSmr1AggW90ZcSVON+xdrwnUdXvWaDJ14771sn7PGtnVaL5BVvvKDtJQd6TUd4gJ1I59NKnap+Ph0rN9b+nTsa/UlNvfICtTj8qCQQn8fXubgI2JWE3KKhmrpKw6tm2blDUnb5CVWH6kslJauh8WL2Z79JCItlJWrZgtbjGLbM0043gfAukJIFuZIRBwRgDZ6owfuSGAbGUOZEMA2ZoNNe/yIFsdsiZmazxANsgiZquVU4qYrVYoxafZ0fiNLFv+gWzasE72NhyKClW1MnXPeu2R/1BLKKnQikCF1Hc4XPrXNcjhHRs0sXp4XX/p3+kI6dfpCOlS01VUzNY+/RrkqKOH2m9UihxPPHq/XHjpdKmpbee4zMDu3W3yVZOwWyS4fZs80qWjXP/2O23v7dplra6KCgkrMatCGKhwBj1ioQz0EAaxEAdbKwKy8KNF8r2Lr7BWrkupiNmaHiQbZLk00WwWQ8zWM20SK63kxGzN/XgTszX3jN2qgZitbpEs3HJ+/8iv5eLLr5bq6hrLnUC2WkZFwhgBYrb6fyrkTbYqNCqMQMPhveS8SSf7n1SKFiJbka06ATbImi3TZ95o6VxGtiZjUo/yb9izPvqYv1qlujv6mL96/F/9p0IBjJARMlAGyjyZ11pAt5ru0q9OydMG7b/+nY+Ufh0bpF9dg7ZitbysPO2Y+F22mjXebIOssubmtjAFWgiDaPiC1tAF+rHtarVsc8Z5uqW+Xl6YPFmu/NM8LUxBS9fu0Z9dukpLt27Jv3dRx6LvS3l65ukqR7YiWzNOzjwkQLYiW9NNO2Rr7k9KZGvuGbtVA7LVLZKFWw6ytXDHrpBajmz1/2jlVbauXr9Z/jDvDZk18xKpranyPy2TFiJbka3IVpFwOCyPP4xsTXcRe3v+S1LTrYNIj4CsV/FTd6+RDfvWxWKnrpO/H9wmEYkkFaEe4+9X119biTo4fIx0a+wq/b99jCZTlVxtV9He0bWzWGSrHQiBXbtShDCIxZXdtlXUaLw4bpxc/fDDdooWKSuTlk6dpEXJVyVllZyNydqI/rv2MyZnu6qf3SRSFf0MRLYiW+1NOG9SI1uRrchWb861VLUgW/PL307tyFY7tIozLbK1OMfVb71CtvptRJLbkzfZqmK2zrz5Plm+ao0ppeGDj5QH77heutR19D1FYrb6fohooI8JFFPM1n3Ne7XH+zcqmaptSLVGNu6J/vvr/RskGE5eURkoC0iv9vWaTO2vVqd2Vj+PkMM7RR/3P6xdTx+PXvE3TYsdu3OHBHbtlMCOHa3/LtP+/Y222Zf2nvqpft+9Oysokdp2bXK2Sxdp6VQnkY6doj+7dJFIx45xx1rq6iSipG6nOul6ZB/ZcyAowVBLVnX///beBD7K6uzf/yZACMi+i2wKoqwiiKWvqIDiAogr7tZCRYRqW7Hwiq3176+t+mpV3tZXpVTEXbGldQNUsIAWVFDAsKnsiyyySNjJJPl/zjMLs2XyzJZ5nplr/PgJmTnPOfe57vNMkisn9+EiCOQ6AWq25voKYP7JEqBma7IEuT7XCVBGINdXQGLzp2ZrYtyq6qqMydaqmmA84xw+ckwPPDbFuuTBcSMCu22nz5iv+x/1Pj/4gj4hr5nnkK3xUKYtBEIJuEm2espKvIdOFW/wSlSrbqr3UCrzp//FR/dFTW/9mg18f+Zv/ty/rdrWP8WSq2bHauu6bVU9vwbLIlsIlJVZwtUraI2QNRLW+zHPL2d9n5vXq+3erbwf9socGpbMo7ywloIFrCVi69YPPFdev4G147a8nv+jV9QGnqtf39qVywMCuUgA2ZqLWWfOqSSAbE0lTfrKRQLI1lzMevJzRrYmzzCdPSBbfXT9ovW9OZ+GCNVFS1fr8UnTArtsTZ1Z8xg76tpAXpCt6Vyi9J3tBJwmWw+UHNC6H77Vhn1rrT/xX7f3W+8u1eJ12nZga4XpaFe/vXUQlbVD1ZKpbdW6XjudXL+96hbUy/Y0Mr8kCeQXFwd2yuYdOOAVtfv2Ka+4WHn7i5W/7wflFe8LPJdv/r2/2Ptcgrtpw0Mur32CT776yiCccILKA0K2nsobNFRZXSNs/ZK2nsqMxPU/Z4QtDwi4kACy1YVJI2RHEUC2OiodBONCAshWFybNASEjWx2QhBghZFS2BgvOE5s31qRH71HL5k2s3aV9enWu0oOz/Id1GVaffrEysHs1/BCvcPlKzdbQ1bWyaInMQQnnnD8wYtktWbxQHo9HvfucG/HaooXzVaOgQD169Yl47ZO5H6hRk2bq3LVH1KW8ZdN6FS1drEuHDqtwqU+dNFE3DR9jjRHtMW/OTJ3YspU6duoW9fVNG9Zq5fKlumTI1RWOwQFZ7qnZuuvw95ZM3bhvvdb/sMb603/zv9mpuufI7ogcj9VY/bPwLTVu0NSqlWr+vN8cSmX+7N8I1RPrnKQ8xd4VWNkas/OlYs03K7V5wzr1v2iInea22uRizVY7YDJVcy3Rmq3FO/eqdI8RssWyRGyxEbXmoxGyxcrb94NXzpqPVhufqPU/d+CAVBZZhmDmoEFqvGuXzv7888qx5eWpvE6dUAFr7ab1Sdl69STfzlprR61/t63ZgesriWCuj/V46+8vq0/fAWreomXl8WS4xfQ3XtD5Ay5R46buLwVCzVZqtsa6nTggK/1vNtRsTT/jVI2Qqe8fUhU//SRPgJqtyTOkh8oJULO1ckaZbpFR2eoXmZcO6KPHnnldN111odq3bSkjNN98Z27En+unC1bwblVTMsAvW8144eLXHOr1m4cn648TRlqxIltDs4JsHav8/PyYSzXVUta80Y68c3ylt4fddqaj8vJymThv+/m4Svs1DZx6QNbW/Zu1avdyfb17pTZtWqu6O2travlUHfYcqnBeHRt1UoeGp+nUhqdZHzs0Ok1F7y3U5dfcojp1Eq8hjWy1tZRsNTp27Khem/qsbr39l7baJ9IoUz8sJSpbU1GzNULAFhdr/pYNaljiUc/DR5W3b69P5vpErSVt/YK3WHkHDySCOuIabzmE+iE7Zr3lDurrlZbN1U/V1fyEOpbYLTe7b0+oozLfx/LatX2fe1/P5APZmkn68Y1d0c7WVSuWaffO7erbH9mKbI1vTaW6NbI11UTT11+mvn9I34zoOV4CyNZ4idE+EQLI1kSoVe01GZOt5oCsCQ9N1rgx11u7WYNlqxGajz39uh6+b2TaD8gycnXD5u2BsgDRZOuwy/qpd4/TrcxEk62/+93vIs4Q//rrr7Vs6VJde9111nUPP/SQxo0fr+rVq1uf/88jj+jusWNV4NtpuXfvXr3y8su68667QlbA8qIifbtmja688sqQ52fNmqXGjRqp99lnhzz//JQpGnjRRWrVqlXI87//f/9P9//udyHPbd68WXNmz9ZPhw+PWHVr167VZ59+qhtvuinitYpiMg0XL1qk73ft0qWXXhpx3Scff6wSj0f9+/ePeO2jjz5SzYICndO3b8RrM2bMUPNmzdTrrLOi3h2xYvVfEM47vKO3335bbdu00Rk9ou+e/fbbb/XF4sW6/oYbKrxD//iHP2jCffdVKltT3S5abqMFabedudbIVhPnb++/39Y7kpGthvF9v/mNrfYlJSV6/E9/0r0TJijfnN5eXm7rOtNo+vTpOq1jR3Xp2jVwzcGSg/pq+zIVfb/c+3Hnci3fWaQfjhw/rOgUnaK+6qsX9aIa12qs05qcrk5NTtfpTTvp9Can67TGp+nkBifLHFgV/vjfiRM1fMQI1TM78xJ8VLbG7HQb696Ldr2ntFzVq8Xecfuvf/1L7U85Rd26d7cTgq02Jrejx4xR7dq1bbWPt9HRo0dlcjL+v/873kttt9+2bZvee/dd3TZypO1rUtHwq2XLtH7DBl1++eW2uzP3kLmFyiO+EtnuosKGFX2tiXqBCWLfPumHH5RnPpr/i4sDz/lfsz762pnXA23NcwcPxgz6udtu08WzZqnVli2VT87UoDVr0EhXI2Z9H3XCCYHnrNfM53Xrhj4X1DZwnXnOvAdUq1b52JIm//WvumzoULVo0cJWeyc3Wr9+vf7zySe6+ZZbUhqmx1Ou6tUzXys4z6yVKPfQl198oW3bt2vw4MEpnXe2dfb999/rH3//u+4YPTrbpuaY+Rw8eFCTnn1WY++5xzExBQfi/+uedHwdyuSE7XwfFR5fpr5/yCQnxg4l8Nijj+quX/xChYWFttFk6z1kGwAN4yZgfrZ/4IEH4r6OC6qOgCNla1XubDW7Wp97bUYEcXMQ1r133axH/vJySEmDcNlqLjxw2GP9zwMCEIifQNMGNbW7+JjKyuwJ12/2rNKq3Su0atdyrdz1lfXvLfs3RR24To26Or1JF3Vu3FVmx2rXZmdYu1UbFzaJP1AXXmF+6KmsvIELp0XIYQQa1StQ8SGPPJ7IMgBuhJVvDg3b7y+B4K9T6/vclD3Yv9/aRZt38GDIx3z/cwd8rx05nLbpW7tvze5a63/v7lrzsczacev73OyELyxUec2aKq9ZKJmP1ueF1kfrc9+/va8XeD+3nqtp9Zn1jzxjODMvWxvUqaEjx0p15Fh23ENZv26YoOMI1K1d3frF+cHDpY6LLZmA+D4qGXpcGw+B2oXVrA0SxQdxCvFwy/W2LRrZF/q5zioT88+YbDWT9e8infCLm/SXKf+0ygg0alBXo+99Utde1q9Ka7b64QfvbK1VWKDKaraa6zggKxNLlzGzhUBFB2R9f2iHV6ruLrLEqvn/mz2rdazsaMTUq+dX1ykNOuj0Rl3VuWk3nd64izo17qpWddtkCybmAYEKCTSpX1OpKCOQdYjLyy0hm38oXMoGfW4OIzO7aQ8YgXtQeYcOWp/7Za5X4B4Xu1ZbTxX+IFS9uk/ABonamgVSsJQNCNua3uf9Qtf85UytWgF5G3GNT+qqwN93uPytqfIK6pxn21rhgKxsyyjzqWoCHJBV1cQZL9sIcEBWtmW0aubDAVlVwznRUTIqW03QZhfrT3/1SEj8UyfeG/iz/UQnluh14bI1/ECs4Pqu/jGQrYnS5joISPXrSp+sX6IV3/uk6u4ird61QruP7IqKp1ntFurUxCtTOzXupk5Nuqpjw9NVo1r0w89gDIFsJ4BsrdoM5x09EthVm3/wUISM9crZA8o7dEg6fEh5R47Kusb87/u3jpjPjyrv2FHlWf/2vqbgdmncmRsPscBOW//u3IgducE7dn3CtqCmZHbnFtbyCeBgkRskhS3R63vN9B/0b+v5WukpRRI+f2RrPCuCthCIJIBsZVVAIDkCyNbk+OXq1chWZ2c+47LVaXjCZauJzzx3/6NTrFBNeYEHx42Q2fVqHhyQFZpBDsjigKyK7mnzp1ibitdbO1RXf79ChYulV+q/rg371qqsPPJPN2tVr2396b8pAdCpqVesdmnSXV/M+1ht23VQ+46d4nr7SOY07VenPqOh19zMAVk2ib/03FMaduMIFaZJlHBAVmgi0ilbF8yfrfoNGqlL9542s5/eZm/9/WX16TtAzVu0TO9AKeg9FQdk5ZWUSJVJ2SBBq2B5e/SodPiwT/QelVfw+oSvaRcseP3y13rO284a14wvaf0pp+jjvn31kxdfTAGZ+LqwdtcW+Eos+Hft+koyeHfyHpe15b5SDCHlG6zdvcdlr3eHb2jJhroN6+pofg0dqVYjqLRDTa1a+4127f6eA7IqSdnePbs0Z9bbuubGEfEll9a2CXBAlm1UGW/IAVkZT0HGA+CArIynICcC4IAs56c5o7LV7BLdvnNPiLw8fOSYHnhsSkidVCdjRLYiW/0Ennv6Txp+B7LV8Cg+uk8rTD3VPcu18vsird69Ql/vXqlDHu/hNwUq0K/1az2kh6wDqdrUO9naodqpURd1btLNqrParn77qPVGP/rgHWTrRUNS9rY498P3dFKbdjr1tC4p6xPZmjjKb1cv19Ytm9TvwkG2O0G22kZVpQ1TIVurNOAKBss7cljfbVinJV99ocv6nHd8p26QoLXErCV9w3bqmt25QTt6vTt5vW0Du3uDBK+53no+uF0cByimmtcXZ52lbS1aaPD770vVa0jVq6ncHHRavYbvY3WVBz9frbrKa5h21SXr36ZtdZVXqy7V8F0T5XmrrdWv6d93vdWPf5zjz4eOX803jvea4+OEPh8Ri2lr+vePW8N3re8Q13g5IlvjJRZ/e2Rr/MwydQWyNVPknTMustU5ucjmSJCtzs9uxmSrX6oOu6xfRMmAqjwgK9kUIVuRrbkuW81u1SXbF2vHh+s0r9V/tHL3cu08tD3qrdWgsKG1Q7VrozN00somOuuqfjq/w1k6cDDf9gFZyNZ16o9s1WtTn9Wtt/8y2bfwCq/P1A9LyNbYKc21na1pW+BxdpzMXwbEOVTU5vnFxT5p6xO6fkF7zMjZyDINXpnrl7pmd+8x307eoF27QTt9q3uOqvzwEZWH9fVl506WbB3y7rupmIar+iivfYJXFgeEbHQ5a+TvrgYN9c8zz9Btn30u5edL1aqpPC8/8G/zXLl53vq/mu9jnsqr+f7taxv43Nfu+DWmT18fpq3/Ov9HX59W+2r+cfK9MQTa5EcZ73hbE1foeL7PQ/rI87XxzyGoTzN2XvDneV4OvucDXMI4eMf0tTVlNyp4IFvdc/tk6vsH9xDK/kiRrdmfYyfMENnqhCzEjiFjsnXvvv2a8NBkjRtzvdq3Df1zwLUbv9NjT7+uh+8bqYamoKPDH9RsdXiCCC9lBPYe2aPF2z/VVzu+1Jc7FmnZzi9knov26NzYu0O1c5Pu6uQ7sKr5CSdGNK3ogKyUBU1HEMhyAunc2Zrl6JgeBCwCldVszTt2TCr1SJ5Sb2mFUo/yzEFpnhLfx1LleUqsw9PCn/c+F/5aqbcP05fpJ/DvUuVZ43gk35hWfyWewPMh45f4rrVi8Y9jYgmO7fjzEbFY4/qu9Y/LmnAOAZ+YLQ8I6iDp7Je4fqHrk9yW2A6Ibq9M9grd44LYei5IdB+X3HkBIR4Qz752AWkc+NwriC2JnJevGgVmd3Q1lZiqUNaYweOFiXVzXZAEDxHdPlkePl4gZv/crDlUIstDhL9flB9n4ZXzQXEGZLmXg8UlaLyAPA/iGSzLLdnu7zMvzznriEhcQYCara5Ik+OCpGar41ISElDGZGu27Gw1NJGtzl7kRJcYgYMlB7R0xxeWUP1i+2f6aucSfXdgS0RndWrUUbdmPdWzRW+rBIC3zmo324MiW22joiEEohJAtrIwIJAcgcpka3K9u+9qUzrCErwhItgnfIOet+StXz6XlkmlRiKXKs/UYS/z/W8+9/+7rNT3fLn3OdPe1zbP/DvQrkx5gbamnfm8LNDWus5/va+dd4xyb59lZd4Ygvr0vn58jDwTr3+88jKFjO+P2T8P07a83NcmqI9ADN7YjvdR7o0hZA5BHHyxBcdsMeeRUwTM4YMyUjbkf1mfl0c8b+RtaNuobcKvUwV9Wd0d7y/QV/AYJhtR4rAVmyXhTf/R+zAvVBq/b/yQdmEMjsenivsLuqbSMS1JHiu24/OJGZeZv//h/3ewgPf/O+i5ghrVrE3uRzzlviuDhH2U9tb4/hz5xwo85+8idh8Wv6A+An16n4zRf4x+K5lnRNzR2vvHDp5fjH5D4vaHFiuOOMcM7T91XGLyroSjyZ2nbTs173FaTr1vum2yGZOtBpQpFzDh4cma9Og9gd2tZlfrqPGPa8ytl+uqQee5giey1RVpIsgYBMwBVd/uWa0vd3yuL7Z9pi92fK41e7+OOLiqIL+mOjftph7NeunM5r11RvNe6tCwY9TaqnaBI1vtkqIdBKITQLayMiCQHAFka3L8uDpNBIywDRLSXsHsE89B8jggs31tAyLbyGdL+oaL5iCxHZDE3rZ+AR7SZ7Q+/GLbJ6RrVZPKS0t19KgnimgOitkvwYOFuYkhIOe9Ejw8Zuu5CDnvi9k3v4iYA7LczDeK4A+W5SFy3MctmL1f4of94sDL2ifu/fOwrvMLszStDbqFAARynsD+//6t6j7y+5zn4GQAGZWtBoxfrm7bsTvAaerEeyPquDoVIjVbQzOzsmiJzEEJ55w/MCJlSxYvlMfjUe8+50a8tmjhfNUoKFCPXn0iXvtk7gdq1KSZOnftEXUZbNm0XkVLF+vSocMqXCZTJ03UTcPHWGNEe8ybM1Mntmyljp2i78jctGGtVi5fqkuGXF3hGG46IMv86b/ZrWrEqhGsS3cs1oGSA4G55Stfv9Vv9XqjN9WjeS/1aHaWejQ/S12adFP1/BoRDEpLS/XCXydqxOh7bN2qnpISvTTlKQ0fdbfila3UbKVm67FjR6nZGnSnpVO2Lpg/W/UbNFKX7j1t3dvpbkTN1nQTjt5/pmu2pnvWFcnWVSuWaffO7erb/+J0h+Dq/jkgK/3po2Zr+hmnaoTKarZaNaWNjA35X97d0xHPG2kb2jZqm+DrfO2jt/OO4/8/0CZsjPD44osrdIyQvlTBHMPmHTJehbEpNbwsMR4jLrMwrCZebnZie7qwukYcPKrCYOnu/3e051SuGtXyrZ2tR0tKvUvR184az/8I7yPY6cfqP0ofx/uN0X9QHIFfINicU8y4LaC+h524g8YM9Bvv3EN+AeK7OBYXm/MM+cVKRH6OB2mLd6VjhjJ75KIL9MADD6TqrYt+0kAg47I1DXOq0i6RraG4ka1jlR/8pyNRVmOqpawpjj3yzvEVrntPWYlW7CrS59Nma1nHVfpy+yJtLF4X0b5RYWP1bHG2ep34I/Vs1lur316s234+ztb9hGytHFNlQr/yHqQ136zU5g3IVmRr6GpBttq5e6q+zfQ3XtD5Ay5R46bNq37wFI+IbEW2xlpSyNYU33BRukO2pp9xqkaoTLamahz6cS4BDshybm6yKTIOyHJ+NpGtSeYI2Yps9RNItUS121+4bN12YKuvHMDn+nL75/rq+yU6WnpE/5/vPxOv2Z3auXFXS6z2avEjq95q23qnBJJZXl4uMz6y1Yvk1anPaOg1N6tOncQP7EO2JvlmG3Q5shXZ2rxF6MGaqVtdqesJ2Zo6lunuiZ2tyRFGtibHz87VyFY7lJzRBtnqjDxkMgpkaybp587YyFbn5zqjsnXvvv0afe+TKloVucuuW6dT9Mwjd6th/cTlRlXhp2ZrVZFmnHACRqIu2/mltVvVlAUwcnX7we8iQLU4oaV316pPrJ7RrKdqVit0BNB4ywg4ImiCgICDCKRzZ6uDpkkoEEgbAWq2pg0tHecIgXq1a6isvFwHDntyZMZMEwKpJXBCYXVVr5anfQdLUtsxvWU1gZaNa2X1/Nw+uYzK1icmTbP4jR11ras5IltdnT5XBb9h39oQsbpyd5E8ZaHf2BqJ2r3pmb6SAGerZ/MeyHkbAAAgAElEQVSzdWKdkxw7T2SrY1NDYC4hgGx1SaII07EEkK2OTQ2BuYQAstUliSJMxxJAtjo2NY4ODNnq6PSYGs+ZOS7R7Gqd8NBkjRtzvdq3df6fA8ZKI7LV2YvcrdHtP7bfOrjKHGBlDrJasmOR9hw5fpCcf15t6p2sXi16q6cpCdD8RxUeYuVUDshWp2aGuNxCANnqlkwRp1MJIFudmhnicgsBZKtbMkWcTiWAbHVqZpwdF7LV2flBtiaZH2q2hgLkgKzED8gycnXh1vlasGW+9fGqXZfrD/qDylQWgFynRh2d0byXtVvVOsiqxdma/rcpMQ/I8l9c2UFawZmkZmvouqZma3xvlC8995SG3ThChbVqx3ehzdbUbA0FlU7ZumD+bNVv0Ehduve0mZ30Nnvr7y+rT98BomZrejmH984BWRyQFWvFUbM1/fcjNVvTzzhVI1CzNVUk3dsPNVvdmzs3RU7NVudnK2Oy1aAxZQTatW6hqwad53xSFUSIbEW2+gnYPdDK3+6g52CIXDUlAcrKj4vV3+l3eqPhmzqzRe+AWO3YsJPy8/JDoNuVqHbbmc6RrcjWZN6Uka2J0/t29XJt3bJJ/S4cZLsTZKttVFXakAOyqhR3UoNxQFZS+IRsTY6fnauRrXYoOaMNstUZechkFMjWTNLPnbGRrc7PdUZl69qN3+mV6bM1bvT1qlVY4HxaUSJEtiJb7crWw55DWrB1vta9+5X+3uifWr57WcSK6tS4q/q26q8+J52jrbO+0fA7Kt8pa1ei2m2HbI280dnZGt/bM7I1Pl7BrZGtsdmxszXxtZXMlexsZWdrrPWDbE3m7rJ3LbLVHicntEK2OiELmY0B2ZpZ/rkyOrLV+ZnOmGw1NVtH3/ukilati0qpW6dT9Mwjd6th/bqOp0jNVsenKGMBLtv5heZv/khzN87Wp999HFOuntOqn+oW1MtYrJkamJqtmSLPuNlCIJ07W7OFEfOAQCwC1GxlfUAgOQLUbE2OH1dDgJqtrIFECFCzNRFqVXdNxmRr1U0x/SMhW9PP2C0jbDuwVbM3zNTHm/+t+ZvnaP+x4pDQT67fQf3bDNS5bQboxyedm5NyNTyXyFa3rG7idCoBZKtTM0NcbiGAbHVLpojTqQSQrU7NDHG5hQCy1S2ZclacyFZn5SM8GmRrCvKDbE0BRJd2cchzUJ9snmuJ1Xmb5mjdD9+GzKRhYSOrLMD5bS5QvzYDdWKdk1w60/SFjWxNH1t6zg0CyNbcyDOzTB8BZGv62NJzbhBAtuZGnpll+gggW9PHNpt7RrY6O7sZl62Llq7WT3/1SAilqRPvVe8epzubnC86araGpmll0RLroIRzzh8Ykb8lixfK4/God59zI15btHC+ahQUqEevPhGvfTL3AzVq0kydu/aIuia2bFqvoqWLdenQYRWumamTJuqm4WOsMaI95s2ZqRNbtlLHTt2ivr5pw1qtXL5UFw2+UqY0wLzNczR/4xx9seMzeco81jXmQKtH8v/HOszKyNXzWl+gbk17RBxoFe9BWvn5oQdihQdotxar3Xamfw7ICqVMzdb43o6p2Rofr+DW1GyNzY6arYmvrWSupGYrNVtjrR9qtiZzd9m7lpqt9jg5oRU1W52QhczGQM3WzPLPldGp2er8TGdUthrR+vikaSG1Wc2hWaPGP64xt16uqwad53iCyNbsl61b92/W7CXvaseaLXq+dIqKj+4LmfSpDU/Tea0vVJvlzXTDbaN1QkGdmOsW2Sp5Skr00pSnNHzU3Yp3Z+tHH7yjtu06qH3HTnG9PyQjC5CtcaEWsjU+XshW+7yQrfZZpbJlMu+fqYwjXX1VtLN11Ypl2r1zu/r2R7YiW9O1+uz1i2y1x8kJrZCtTshCZmNAtmaWf66Mjmx1fqYzJlsPHzmmBx6bomGX9YvYxWok7JvvzNWD40aoVmH0nYhOQYtszT7ZeqDkgD7Z/JFVFsD8v7F4nTqqo87SWXpVr6pxYROd23qAzmt7gQa0uUhNaze3IKRaotrtz+6OVbvtzFzY2Rq6rpGt8b3jIlvj44Vstc8L2WqfVSpbIluRrcjWVN5R8feFbI2fWaauQLZmirxzxkW2OicX2RwJstX52c2YbN27b78mPDRZ48Zcr/ZtW4aQMrtbH3v6dT1830g1rF/X8RSp2er4FMUMsLS8VF/u+NwqC2DKAyzdsVjmOf+jIL+merfsY9VcPa/NBerSpLvylOfuSTso+nh3tjoodEKBgCMIULPVEWkgCBcToGari5NH6I4gQM1WR6SBIFxMgJqtLk5eBkOnZmsG4dsYOmOyNVt2thrGyFYbK81hTdbvW6P5m8zu1dlasHWe9h/bHxLhaY076/zWF1q1V/uc1FeF1Wo5bAbZEw6yNXtyyUwyQwDZmhnujJo9BJCt2ZNLZpIZAsjWzHBn1OwhgGzNnlxW5UyQrVVJO/6xMiZbTajTZ8zXtHfmurpmK7I1/kWXiSv2HyvW/M1zNHfjbH285SNtLt4YEkbTWs10bpsBlmDt13agmtRqmokwc3JMZGtOpp1Jp5AAsjWFMOkqJwkgW3My7Uw6hQSQrSmESVc5SQDZmpNpT3rSyNakEaa1g4zKVjMzU5/1p796JGSSUyfeG1HHNa0Ukuicmq2h8FYWLZE5lfac8wdGUF2yeKE8Ho969zk34rVFC+erRkGBevTqE/HaJ3M/UKMmzdS5a4+omdqyab2Kli7WpUOHhbz+zZ5Vmr1hpj5Y954Gbu+vP+lPOqZjVpua1Qr1o5b/pfOt0gAD9P3SLTqxZSt17NQt6hibNqzVyuVLdcmQqytcLXZrrKa6nd1arHbbmQlSszU0zdRsje9Nkpqt8fEKbv3t6uXaumWT+l04yHYn6ZStC+bPVv0GjdSle0/b8aSzITVb00m34r6p2UrN1lgrz3zfN2fW27rmxhGZWaA5MCo1W92TZGq2uidX6YqUmq3pIku/wQSo2er89ZBx2ep8RLEjRLY6R7YOHHKFFmz9WB+un6HZ62do0/4NgeDu0316t9EMndO2n85rc6H6nHSOTC1W/2PenJnI1qBUIluRrcm8NyNbE6eHbI3NDtma+NpK5kpkK7IV2ZrMHZT8tcjW5BlWVQ/I1qoi7dxxkK3OzU02RYZsdX42Mypbn5g0Tdt37tGD40aoVmGBRctfy7VPr866atB5jieIbM2sbD1YckBvff6Gtq5er8mev8p8HvwY0PZiDTx5kEr/c0A3D/+5tXs22gPZGkoF2YpsTebNF9maOD1kK7I18dWTviuRrchWZGv67i87PSNb7VByRhtkqzPykMkokK2ZpJ87YyNbnZ/rjMlWDshy/uJwaoSHPAf17pp/aubat/TB+vdCwjS1Vi8+5TJd0O4S63ArDrZyahaPx0XNVufniAidTSCdZQScPXOig0BqCFCzNTUc6SV3CVCzNXdzz8xTQ4CaranhmGu9ULPV2RnPmGzdu2+/Jjw0WePGXK/2bVuGUFq78Ts99vTrevi+kWpYv66zCUraf6hE+w97HB+nmwM0gnXm2rf17prpEYK1Tb2TdcWpw3TxKUPUo/lZbp5mTsaObM3JtDPpFBJAtqYQJl3lJAFka06mnUmnkACyNYUw6SonCSBbczLtSU8a2Zo0wrR2kDHZys7WtOY1Kzo3gtXsXH3n23/o3xs/1NHSI4F5nVy/gy7rcJUu63i1OjeOfqhVVkDIgUkgW3MgyUwxrQSQrWnFS+c5QADZmgNJZoppJYBsTSteOs8BAsjWHEhyGqaIbE0D1BR2mTHZauawaOlqTXh4siY9ek9gd6vZ1Tpq/OMac+vl1GyVtOabldq8YZ36XzQkJO0VnRBd0eEh0U6i37Ftqz5bMFdDr74pYklt2bReRUsX69KhwyJeqygm03Bl0RKZU2nPOX9gxHVLFi+Ux+NR7z7nRry2aOF8q55qxzO6WQdcvW0J1g8swTpEQ7Rd27Wr/l6vYD31anVp0j3QR6xY/Y2mTpqom4aPoWbrneNtvX1QszUU06tTn9HQa25WnTqJ77SvrC6wncTEuvfsXB+tzdwP39NJbdrp1NO6JNpFxHXUbE0cJTVbY7PjgKzE11YyV1KzlZqtsdaP+b5vzqy3dc2NI5JZZlwbgwA1W92zPKjZ6p5cpStSaramiyz9BhOgZqvz10NGZavB45er23bsDtCaOvFe9e5xuvPpSeKArNA0JSJbD3sOafoHL2v1Dyv06r6XdaT0cKDTtvVO0Q01b1Dntmfogj6hwtnfCNn6qEbakKjRhHtFNxmyFdmazBswsjVxeshWZGviqyd9VyJbka3I1vTdX3Z6RrbaoeSMNshWZ+Qhk1EgWzNJP3fGRrY6P9cZl63ORxQ7QmRrYrLVCNYPN8zQO99O10cb3lff0nN0REf0iT6REaxDTr1SQ0+9Rl2bnKFP5n6gRk2aqXPXHlGTgWxFttp5H0lGFrCz1Q7h422QrfHxCm6NbEW2Jr560ndlMu+f6YsqdT1XVEZg1Ypl2r1zu/r2R7YiW1O33hLpCdmaCLXMXINszQx3J42KbHVSNrI3FmSr83OLbE1Bjjggyx5EI1hnb5hplQgwgjV4B6s55GpIB69g7dY0ulS1Nwqt3EaAmq1uyxjxOo0ANVudlhHicRsBara6LWPE6zQC1Gx1WkaIx20EqNnqtow5I15qtjojDxVFgWxNQX6QrRVDNIJ1zsZZevubf+ijje/LfO5/tKnbztrBetmp16h70zNTkAm6cCMBZKsbs0bMTiKAbHVSNojFjQSQrW7MGjE7iQCy1UnZIBY3EkC2ujFrmY8Z2Zr5HMSKANmagvwgW0MhmkOtrBIB30y3RGuwYG1dr62GtPcK1jOa9UwBfbpwOwFkq9szSPyZJoBszXQGGN/tBJCtbs8g8WeaALI10xlgfLcTQLa6PYOZiR/ZmhnudkdFttolVUE7arYeB2ME6yefz9aO77fqnfJ3Ai+0OKGlrug4TL1KeqlpYXP17nNuBM1FC+erRkGBevTqE/EaNVvzY65Suwdf2W1nBuOArFDk1GyN742Smq3x8QpuTc3W2Oze+vvL6tN3gJq3aJk45Cq6cvobL+j8AZeocdPmVTRi+oahZis1W2Otrr17dmnOrLd1zY0j0rcIc7xnara6ZwFQs9U9uUpXpNRsTRdZ+g0mQM1W568HZGuSOcp12bp4+6ea/vXr+tc307Tv6A86W2eriZpoWd0iXdbhKg3pcJXObN7borxk8UJ5PB5k69N/0vA7xio/P7ZEfc5mO7sS1W47ZGvkmwKyNb43SmRrfLyQrfZ5IVvts0plS2QrshXZmso7Kv6+kK3xM8vUFcjWTJF3zrjIVufkIpsjQbY6P7s5L1sPHzmmBx6bovfmfBrI1tSJ96p3j9MDn0+fMV/3PzrF+nzwBX304LgRqlVYYH2ei7L1s2XztHjlAr1S8oo2Fq8LcLqg7SUafMIQ1TpWS0MvuSFi9SNbvUjsSlS77exKVLvtkK3I1mS/dCFbEyfIztbY7JCtia+tZK5EtiJbka3J3EHJX4tsTZ5hVfWAbK0q0s4dB9nq3NxkU2TIVudnM+dl6959+/X86zM1+tYrLIG6aOlqTXh4siY9eo/at21pff74pGl65pG71bB+XT0xaZqV1bGjrg1kN1dqtv7r22l6ZfkULdg6PzD3Xi1+pOs63WLtYK1fs4HzVzwROo4ANVsdlxICchkBara6LGGE6zgC1Gx1XEoIyGUEqNnqsoQRruMIULPVcSlxRUDUbHV2mnJetoanx8jX0fc+qXtGXWvtbjVytV3rFrpq0HlW03D5ap7LZtn6zZ7Veqlost78+hXtP1ZsMWhaq5lu7vozXdvpZrWpd7KzVzjROZ4AstXxKSJAhxNAtjo8QYTneALIVseniAAdTgDZ6vAEEZ7jCSBbHZ8iRwaIbHVkWgJBIVvD8rN243f6zcOT9ccJI9WyeROrxECfXp0DsjX4dbPzNRtl62HPIasG6ysrnteSHYusOebn5atfm4G6pettuqDdJaqWV83ZK5voXEMA2eqaVBGoQwkgWx2aGMJyDQFkq2tSRaAOJYBsdWhiCMs1BJCtrkmVowJFtjoqHRHBIFuDkPjrt/rlqv/zYZf1C9RwDZetpmbrfb+5XyWlZSFw13z7jYq+WqYrrx5mPf/kn/5Hd/3qHlWvXt36/H+feEyj7/ylCgq8tV9/+GGv3nz9NY28Y0xIPytXLNf6dWs1+LLLQ56f8+H7atiosXr2Oivk+Vdemqr+Awaq5UknhTz/2CN/1Lh7fxPy3NYtWzRv7hzdePOt1vOLty3S80v/pjdXTdOJx1qoj/roozof6Sfdh+tnPW9XyzpeuVxRTOa1JV9+od27vteFF10Ssdg+XfAflXhKdO55/SJemz/v36pZUFM/+vF/Rbz24fsz1bRZc/U4s2fUu2nD+nVa9PlnGnZdZJ1Y/wXhvMM7mvneu2rdprW6djsj6hjr1q6x5nb1sOsqvKMff/Rh3f3r/6704KtUt4uW22hB2m1nri0vL5eJ89f/fZ+td7DS0lJrTY8dd6+t9iUlJfq/Pz+pX90zXoUF1XS0pMwa087j3bf/pfYdTlWnzl3sNA+02bhxgz5b8B9de8NNcV1nGj/79F900823qm69eiHX5sXRk1ljrdq0VrcK1pidrsy9t27dWg0Jez+o6NrScqlaXuwo33v3LbVrd4q6dO1mJwRbbZ768xMacdsdql27tq328TY6evSonn36z/rl3ePivdR2++3bt+n9WTN0609/ZvuaVDRcXvSVNm3aoEGDh9rurmaNfB0rLVN56Jch29fHajj7w/fVqFEj9ezlPegw04+XX5qqAQMuVMuTWmU6lErHf+H5v+mSQUPUvHmLSts6vYF5//x0wSe67oabUxpqWXm58it5j0rpgBV0VlAjX6Wl5SotC/06tGzpl9q+fbsuvmRQVYTh2jF27fpeb/9rukbcNsq1c3B64AcPHtTUKX/Vz++625Gh1qieJ/NtnMd845FFD/OeUC32ebYRs83U9w9uwZ5dKyQ69b9MfFwj7/i5CgsLbaelerV85edJxzxp+GbOdhQ0dBMB87P9Aw884KaQcy5WZKsv5X6x2qJZo0A91nD5appGk63j7/2tjhwrDVk869Z+q5XLl2nI5ddYzz/958d0+5i7A7L12ace14jb7wrI1n0/7NW/pr+hW0fcEdLP16tWaMP6tbp4UOgP3vP+/aEaNGioM84Mla3TXntR5/a7QCeeGCpb//zEw/rF2AkhfW/7bovmzftQxacf0otFU7Ry13LrdbOL9erm16p32VkaftOYiF2sFcVkrv1q6Rfas2e3+g24KOJmWvTZAnk8JfrxOedHvLbgk7kqKKips87+ccRr/57zvpo0aapuZ0SXrRs3rNOSLxfpiqsqFqHhvMMHmf3+e2rZqrU6d+ke9U3A5MDMbeiVx2v1hjd8auL/aMwvxlUqW1PdLlpuo03CbjtzrRGfJs677rYnT41sNYx//svxtt5EjWz927P/q9F3/Vr1T6ih4kMl1jfpdh6zZrylU045VR1P72yneaDN5k0btPjzhbrymoqlfEUdTpn8fxp2/S2qWzdUttoM2erWrLGTWrVWpwrWmJ3JmHtv4/q1uijs/aDCa8ulvEp+SPhg5jtq0/Zknd65q50QbLWZ/PT/6ubhI1WrVnpk67GjR2VycsedY23Fk0ijnTu266MPZ+r6m4cncnnC16xaUaQtmzdq4CVDbPdRt3YNHT5aKk/YL/1sdxCj4byPPlSDhpFfa1LRdyJ9WF/jzr9AJ7YM/RqXSF/pvua1l6bowosHW78sdPvDev/8bKGuHBb/+2fMuZs30Xh+a5UmkGZHUYmnLOKH3OVfLdHOHTs0YGDkL5DTFIoru92ze5dmvPNP3fzTka6M3w1BHzp0UK++8JxuG/0LR4Zbq6C6ylUe8fOQI4ONIyi735sGd2m+f/j37Jm67qaq/f4hjmlltKkD3vLTPv9JTz+pn44YrZpxyFbzi/P8/Dzr+zkeELBDwPxsj2y1QypzbZCtkqKJVn9KKqvZana2jh13n/Yf9oRkceP6Nfp65Ve6aPBV1vPPP/uEbrntFwHZOnXSRN00fIxq+Ha2Fu/7QTPfnqbrbrk9pJ8136zU5g3r1P+i0B+8F8yfrfoNGqlL91ABWdFJzeEn0ZtDrt5a/Lpqbq6m5/ScNWbzE07UDZ1u1S3dbpNn91EVLV2sS4d6d+YGPyqKybRZWbREe/fs0jnnD4y4bsnihfJ4POrd59yI1xYtnG+x6NGrT8Rrn8z9QI2aNFPnrj2i3ilbNq2vMFb/BeG8wzuaN2emTmzZSh07Rd/dt2nDWq1cvlSXDLm6wrv1uaf/pOF3jK1Utqa6XXhuKwrQbjtzvZGtJs7bfm5v56CRrS/8daJGjL7H1ruZp6REL015SsNH3a14ywh89ME7atuug9p37GRrLH+jZE7TfnXqMxp6zc2qU6duXGMGN65sjdnpONa9Z+f6aG3mfvieTmrTTqeeFt9O4VjjvfTcUxp24wgVpku2Hjuq16Y+q1tv/2Wi0670ukydJvzt6uXaumWT+l1ofyddOssIVPS1plKAaWpQ0de4NA2XVLfT33hB5w+4RI2bul+2JvP+mRTEKrq4ojICq1Ys0+6d29W3/8VVFIk7hzHf982Z9bauuXGEOyfggqgPHzqof7w+VTeP+Lkjo6WMwPG0ZOr7B0cujBwN6sXJf9Z1P7ldNWva39lKGYEcXSxJTNv8bI9sTQJgFVya87I12u7VYO7hB2IZ+WoeY0cd3+HolgOydh3+Xm+sfFGvrpyqDfvWWvMwu1jPb32hJVgvbHcptVir4KZjiFAC8cpW+EEAAqEE0ilbYQ2BXCBAzdZcyDJzTCcBZGs66dJ3LhBAtuZCllM/R2q2pp5pKnvMedlqygKMGv+4tu3YHcL1ZzcMCgjV6TPm6/5Hp1ivD76gjx4cN0K1Cr21Vs3DybK1rLxM/970gV5d8bxmb5gpT5l3B67ZxXp9p5/o5q4/U8s6zq99l8pFT1/OIoBsdVY+iMZ9BJCt7ssZETuLALLVWfkgGvcRQLa6L2dE7CwCyFZn5cMt0SBbnZ2pnJetqUiPE2Xrdwe26JXlU/TG6pe07cBWa5r+XaxGsA48eRC7WFORfPpImgCyNWmEdJDjBJCtOb4AmH7SBJCtSSOkgxwngGzN8QXA9JMmgGxNGmFOdoBsdXbaka1J5sdJNVs9ZSV65dVntLDmp3pv59syu1rNo1ntFhpz6A4N/umNIbtYd2zbqs8WzNXQqyNPZ49VB5WardEXTaprsdrtz24tVrvtzOyo2RqaY2q2xvdGSc3W+HgFt6Zma2x21GxNfG0lcyU1W6nZGmv9ULM1mbvL3rXUbLXHyQmtqNnqhCxkNgZqtmaWf66MTs1W52ca2ZpkjpwgW0391ReL/qa/r35FVx65QrM0S9/lfafzWl+gW6xdrIM15f8e18g7Q0+KR7YeT35lhxdxQFbsG4UDsip/I6lsjVXeg8QBWV5KxzggK2S5pHNnKwdk2bkzo7fhgKzE2VX1lRyQlRxxZGty/OxcjWy1Q8kZbZCtzshDJqNAtmaSfu6MjWx1fq6RrUnmKFOy9eN5s7Tx2Eb9o/hNfb5tQWAWd+SPVq1T6+vGH48I2cUabVcjshXZWtHyZ2drKBl2tsb3RsnO1vh4BbdmZ2tsduxsTXxtJXMlO1vZ2Rpr/SBbk7m77F2LbLXHyQmtkK1OyEJmY0C2ZpZ/royObHV+ppGtKchRVdZs/eHIXv1t2VN6vuhZmX+bh6nFel6rAbql223WLtZqedVSMCu6gEDVEKBma9VwZpTsJZDOna3ZS42ZQeA4AWq2shogkBwBarYmx4+rIUDNVtZAIgSo2ZoItaq7BtmaAtZVIVt3Htqup794Uq+smKJDnoNW1E1rN9f1nX6in3QbGbKLNQVTogsIVBkBZGuVoWagLCWAbM3SxDKtKiOAbK0y1AyUpQSQrVmaWKZVZQSQrVWGOqsGQrY6O53I1hTkJ52ydf2+NfrLosf0xuqXApGe0ayn7uj5Kw3tcE0KoqcLCGSWALI1s/wZ3f0EkK3uzyEzyCwBZGtm+TO6+wkgW92fQ2aQWQLI1szyd+voyFZnZw7ZmmR+0lWzddH2hZq05M/auna9OqiDpmu6rux4nYZ3v0O9WvxIFR1aUlE9O2q2TtRNw8eoRkFB1IxXdngRB2TFvlE4IKvyN5LK1ljlPXBAlp8RB2SFrpZ0ylYOyLJzZ0ZvwwFZibOr6is5ICs54tRsTY6fnaup2WqHkjPaULPVGXnIZBTUbM0k/dwZm5qtzs81sjXJHKVattb+cWP9delf9OX2z63IelfrrQvqDdS1V/xMzU84MRAtsvV44rZsWq+ipYt16dBhFWZz6iRkazThXhEwDsgKJcMBWfG9UXJAVny8gltzQFZsdhyQlfjaSuZKDsjigKxY6wfZmszdZe9aZKs9Tk5ohWx1QhYyGwOyNbP8c2V0ZKvzM41sTTJHqZCtx8qO6o0vX9DOzzfribLHrYhMPdbbe9yl8wv7aeeW79T/oiEhkSJbka1+AnYlqt12pl9kK7I1mbdGZGvi9JCtyNbEV0/6rkS2IluRrem7v+z0jGy1Q8kZbZCtzshDJqNAtmaSfu6MjWx1fq6RrSnIUaI1W/cd/UFTvnpGzy97RruP7LIiObl+B43u+Stde/rNqlEt+p+8pyBkuoCAYwhQs9UxqSAQlxJIZxkBlyIhbAjERYCarXHhojEEIghQs5VFAYHkCFCzNTl+uXo1NVudnXlkawryE69s3bJ/k5758km9seolHfYcsiLo3vRM3XnWrzWo/RXKU14KoqILCLiDALLVHXkiSucSQLY6NzdE5g4CyFZ35IkonUsA2erc3BCZOwggW92RJ6dFiWx1WkZC4wa98lEAABiLSURBVEG2piA/dmXr8l3L9OdFj2rmurdUVl5mjXxu6/66s9c49W3VLwWR0AUE3EcA2eq+nBGxswggW52VD6JxHwFkq/tyRsTOIoBsdVY+iMZ9BJCt7suZEyJGtjohCxXHgGxNMj+V1WwdOPhKzd4wU2vf+0oPlz8kjzzKz8vXb/J+q3OvvlRdmp9hRVC87wfNfHuarrvl9pCI1nyzUps3rKNma5Nm6ty1R9RscUDWoxp55/hKVzI1WzdqyeJPNfiK6yplFd6AA7LiQ0bN1vh4BbemZmtsdhyQlfjaSuZKarZSszXW+uGArGTuLnvXUrPVHicntKJmqxOykNkYqNmaWf65Mjo1W52faWRrkjmqSLauXbtK8xZ9oBdLp+rbvV/rt/qtHs97Qld3vkF3nTVOs1+drpuGj1GNAm9dVmTrfItFj159IjLyydwP1AjZWuFKtStR7bYzA3FAVihuZGt8b5TI1vh4IVvt80K22meVypbIVmQrsjWVd1T8fSFb42eWqSuQrZki75xxka3OyUU2R4JsdX52ka1J5ihctu4/VqypRZP00ZczdMrRk/W6Xlfdgnoa6xmrK37yEzWr08IaceqkicjWIPaLFiJbw5fic0//ScPvGKv8/PyYq9SuRLXbDtkaiRvZGt8bJbI1Pl7IVvu8kK32WaWyJbIV2YpsTeUdFX9fyNb4mWXqCmRrpsg7Z1xkq3Nykc2RIFudn11kawpyZGq2fvP9Rj3z5US9vuoFHSw5YPXarHYLjexxl27tNlIn1KiTgpHoAgLZR4CardmXU2ZUtQSo2Vq1vBkt+whQszX7csqMqpYANVurljejZR8BarZmX06rYkbUbK0KyomPgWxNnJ11ZdHOIj347z/oX1//Q6XlpdZzJ9fvoDG97taw025SjWreMgE8IACB6ASQrawMCCRHANmaHD+uhgCylTUAgeQIIFuT48fVEEC2sgYSIYBsTYRa1V2DbE2S9aWvXKpZa2ZZvZzRrKd+3useDWp/hfKUl2TPXA6B3CCAbM2NPDPL9BFAtqaPLT3nBgFka27kmVmmjwCyNX1s6Tk3CCBbcyPPqZ4lsjXVRFPbH7I1SZ6mZuu8th/rjjPvUd9W/QK9bVy/Rl+v/EoXDb7Keu75Z5/QLbf9QtWrV7c+p2ZrKHhqtkYuRGq2Rr85P/rgHbVt10HtO3aK6+5NpuYgNVvjQi1qtsbHK7j1t6uXa+uWTep34SDbnaRTti6YP1v1GzRSl+49bceTzobUbE0n3Yr7Tub9MzMRxzdqRbJ11Ypl2r1zu/r2p2ZrLKJ79+zSnFlv65obR8QHnta2CVCz1TaqjDekZmvGU5DxAKjZmvEU5EQA1Gx1fpqRrUnmKPyALH93yNaBEWSXLF4oj8ej3n3OjXgN2YpsLSsrt3U3IlvXqf9FQ2yxstNo7ofv6aQ27XTqaV3sNLfVBtlqC1PURsjW2OyQrYmvrWSuRLYiW5GtydxByV+LbE2eYVX1gGytKtLOHQfZ6tzcZFNkyFbnZxPZmmSOkK2hAFcWLZHZ4XDO+cjWipaW3R2rdtuZN9qRd46vdCXbbWc6Ki8vlxn/tp+Pq7Rf06C0tFQv/HWiRoy+x1Z7T0mJXprylIaPulvxlhFAtiJbjx07qtemPqtbb/+lrfWWSKNM/bCEbEW2JrJe030NshXZimxN910Wu39ka2b5xzN6pr5/iCdG2qaXALI1vXzp3UsA2er8lYBsTUGO9h8q0f7DnhT0RBcQyD0C8crW3CPEjCEQm0A6ywjAHgK5QICarbmQZeaYTgLUbE0nXfrOBQLUbM2FLKd+jtRsTT3TVPaIbE0BTWRrCiDSRc4SQLbmbOqZeIoIIFtTBJJucpYAsjVnU8/EU0QA2ZoikHSTswSQrTmb+qQmjmxNCl/aL0a2pgAxsjUFEOkiZwkgW3M29Uw8RQSQrSkCSTc5SwDZmrOpZ+IpIoBsTRFIuslZAsjWnE19UhNHtiaFL+0XI1uTREzN1lCA1Gwdq/z8/Jirym4tVrvt7NZitdvOBE/N1tAUvjr1GQ295mbVqVM34XeMeXNm6sSWrdSxU7eE+1jzzUpt3kDNVmq2hi6hdMrWBfNnq36DRurSvWfC6zaVF3JAVipp2u+Lmq3UbI21Wkyt/jmz3tY1N46wv6hoGRcBarbGhSujjanZmlH8jhicmq2OSEPWB0HNVuenGNmaZI6QrchWPwG7cjTV7exKVLvtkK2RbwrI1vjeKF967ikNu3GECmvVju9Cm62RrcjW5i1a2lwtmWs2/Y0XdP6AS9S4afPMBZGikZGtyFZka4pupgS7QbYmCC4DlyFbMwDdYUMiWx2WkCwNB9nq/MQiW5PMkSVbx/9GppRA8GPj+jX6euVXumjwVdbTzz/7hG657ReqXr269fnUSRN10/AxqlFQYH1evO8HzXx7mq675faQfirayVbRbqOKdv1EE207tm3VZwvmaujVN0VQ2LJpvYqWLtalQ4dFvBZrdx07W9nZaueW8pSU6KUpT2n4qLsVbxmBjz54R23bdVD7jp3sDBVok4wsQLbGhVrI1vh4Bbf+dvVybd2ySf0uHGS7E3a22kZVpQ2RrVWKO6nBKiojsGrFMu3euV19+yNbYwFmZ2tSy8/WxchWW5gc0QjZ6og0ZDQIZGtG8efM4MhW56ca2er8HBEhBCAAAQhAAAIQgAAEIAABCEAAAhCAAAQg4AICyFYXJIkQIQABCEAAAhCAAAQgAAEIQAACEIAABCAAAecTQLY6P0dECAEIQAACEIAABCAAAQhAAAIQgAAEIAABCLiAALLVBUkiRAhAAAIQgAAEIAABCEAAAhCAAAQgAAEIQMD5BJCtCeZo7779Gn3vkypatc7qYerEe9W7x+kJ9sZlEMheAvHcK4ePHNMDj03Re3M+DQDh3sretcHM7BEIvy9+P36Erhp0nq2Ln5g0TZ8vXa1nHrlbDevXtXUNjSCQbQQSuYfWbvxOo8Y/rm07duvE5o016dF71L5ty2xDw3wgYJvA9Bnzdf+jU6z2gy/oowfHjVCtQu9Bv9Ee5uvPc6/NsF7q1ukUvg7ZJk3DXCRgvuY89vTrevi+kXy/losLgDlnJQFkawJp9X/T3qdXZ+sHXvPm+JuHJ+uPE0byjXgCPLkkewnEe68YMfv86zM1+tYrrG/gFy1drQkPT+aH3OxdIszMBgHzA6t5jB11rfy/vLhn1LWV/oLP/4MuP+TagEyTrCYQ7z3E93VZvRyYXAIEzPdjj0+aFhCmwfdUtO6MmP30i5UBIRv+eQIhcAkEspJA8KYUvl/LyhQzqRwmgGxNIPnhv3kKF0oJdMklEMhKAsneK/GIpawEyKRynoC5ByY8NFnjxlwf+GVeZT/kGmjmB9sNm7fr3B91D/kBOeeBAiDnCMR7D/m/pxt2Wb9Kf6GRczCZcM4SMF932rVuEfirinD5Gg4m/OtUZe1zFiwTh4CPADtbWQoQyD4CyNYEchrtGwY7P/wmMBSXQMDVBJK9V9hd5Or0E3wKCES7ByrbIRT8+vLV65CtKcgDXbiXQLz3UHjpGzNzO38y7V5CRA6B2ASibSqp7PszfxmOQQN+ZP1VRrishTkEIBBKANnKioBA9hFAtiaQUyOQ3nxnbkitImRrAiC5JOsJJHOvsGM865cHE7RBINo337Fka/g9x24iG5BpktUE4r2HKvqLjBbNGlnSiAcEco1AtN3elclW/zX79h/UJ58XUbM11xYN842bALI1bmRcAAHHE0C2JpCiZHfrJTAkl0DAlQQSvVf836Tzw60r007QKSQQ76684ANMgsOgDlgKk0JXriIQ7z0U7QdefmnhqpQTbIoJJLKzNXwnq/naNO2duRySleLc0F32EEC2Zk8umQkE/ASQrQmshWTrUCYwJJdAwJUEErlXEK2uTDVBp4lAvPUmw8NAEqUpMXTrGgLx3kPR2kf7Kw3XACBQCKSAQDw1WxPZCZuCEOkCAq4mgGx1dfoIHgJRCSBbE1gY8Z6wnsAQXAKBrCBQ2b3ir4137WX9rEMXKB2QFWlnEikmEOsk9fB7CNmaYvh0lxUEYt1DZoLm9e079wTKQwV/bl5/4LEp6tOrc+BwoKyAwiQgEAeB8F/chZdPi3UP1SossA5tZGdrHMBpmnMEkK05l3ImnAMEkK0JJjn8AIWpE+/l1NoEWXJZdhOIda+EiyL/gQrbduwOgfKzGwZRKy+7lwmzi0HA/0uI9+Z8arX6/fgRAemDbGXpQKByArHuoWiyNbw9X4MqZ0yL7CcQXKYm/NC4cNkafg9Ryib71wczTIxAtEMZ+ZqTGEuugoDTCCBbnZYR4oEABCAAAQhAAAIQgAAEIAABCEAAAhCAAARcSQDZ6sq0ETQEIAABCEAAAhCAAAQgAAEIQAACEIAABCDgNALIVqdlhHggAAEIQAACEIAABCAAAQhAAAIQgAAEIAABVxJAtroybQQNAQhAAAIQgAAEIAABCEAAAhCAAAQgAAEIOI0AstVpGSEeCEAAAhCAAAQgAAEIQAACEIAABCAAAQhAwJUEkK2uTBtBQwACEIAABCAAAQhAAAIQgAAEIAABCEAAAk4jgGx1WkaIBwIQgAAEIAABCEAAAhCAAAQgAAEIQAACEHAlAWSrK9NG0BCAAAQgAAEIQAACEIAABCAAAQhAAAIQgIDTCCBbnZYR4oEABCAAAQhAAAIQgAAEIAABCEAAAhCAAARcSQDZ6sq0ETQEIAABCEAAAhCAAAQgAAEIQAACEIAABCDgNALIVqdlhHggAAEIQAACEIAABCAAAQhAAAIQgAAEIAABVxJAtroybQQNAQhAAAIQgAAEIAABCEAAAhCAAAQgAAEIOI0AstVpGSEeCEAAAhCAAAQgAAEIQAACEIAABCAAAQhAwJUEkK2uTBtBQwACEIAABCAAAQhAAAIQgAAEIAABCEAAAk4jgGx1WkaIBwIQgAAEIAABCEAAAhCAAAQgAAEIQAACEHAlAWSrK9NG0BCAAAQgAAEIQAACEIAABCAAAQhAAAIQgIDTCCBbnZYR4oEABCAAAQhAAAIQgAAEIAABCEAAAhCAAARcSQDZ6sq0ETQEIAABCEAAAhCAAAQgAAEIQAACEIAABCDgNALIVqdlhHggAAEIQAACEIBAHASmz5iv+x+dEnJFt06n6JlH7taa9Vv10189oqkT71XvHqeHtHli0jR9vnS11a5h/bqK1c+eH/Zr1PjHtW3H7goj+/34EWrdspk1XrSHP4ZFS1dbbQZf0EcPjhuhWoUFgeaxXjON9u7br9H3PqmiVesqjONnNwxSu9YtQpiY2K4adJ7WbvzOmkeTRvUD8/Z3FO01fzyx5hNHqmgKAQhAAAIQgAAEIJADBJCtOZBkpggBCEAAAhCAQHYSCBem/lma58/9UXdLsBqJOu2duSFy0YjF3zw8WX+cMFLt27aUnX6CCZo+P/1iZVRZOuHhyZr06D1Wv9EefoF5YvPGIe0OHzmmBx6bovfmfBpVxFbU1+OTpkWIU9PWjBMei1+oGmkcLqANg+demyG/qDYCOlof2bmSmBUEIAABCEAAAhCAQKoIIFtTRZJ+IAABCEAAAhCAQBUS8O/yvPayftauzYoefonZolkjjR11rfyf9+nV2brObj+plK1GkA4450wdOHjYiskvR998Z67q1KmtAwcORYjcVMlWI5mHDPwvfbVybWAMI2Efe/p19ex2qj76z5KAvEW2VuGCZigIQAACEIAABCCQJQSQrVmSSKYBAQhAAAIQgEBuEQiXqLFm79/R+fCEkdr83c6Qna7x9OMfI9mdrUa23jPqWj3+7Bshu2vNn/9v2Lxd23fuSatsveeO6+SPwb/714xrxg/eBYxsza17itlCAAIQgAAEIACBVBBAtqaCIn1AAAIQgAAEIACBDBAIryka/Cfw4eH4/0zePB/+J/Tx9GOujyVb7dRs9f/p//Ovz7TCvPySvtbO0ofvGynzXLplqymfsGzFGqsUwt23D9ODT7ygcWOut54Ll62VzScDaWdICEAAAhCAAAQgAAEHE0C2Ojg5hAYBCEAAAhCAAATsEAiud2rah9dDNc/5ywWc3eP0wJ/uh/dtp5/KZKudmq1+2WoO3jJ/1m8e1w7tb5U1MFK4KmRrowZ1rcO26tc9Qae1b20xCa9vy85WO6uPNhCAAAQgAAEIQAACwQSQrawHCEAAAhCAAAQgkEUEKioLEF6rtbIpxyovkIoyAs88crfMIVThh3NVlWw1B3iZeTz9wluBg7qQrZWtCl6HAAQgAAEIQAACEKiMALK1MkK8DgEIQAACEIAABBxIwOxULVq1Tuf1OSMiOiMszcN/+JT5d0WyNd5+TF+plK2mnuya9Vt0cb+zrZirUraauc+Y86luumqgNTay1YELnZAgAAEIQAACEICAywggW12WMMKFAAQgAAEIQAAChoC/LECbls1CDpPy118Nr8saS7aaP6e320+qZWt4NqtStoaPjWzl3oIABCAAAQhAAAIQSJYAsjVZglwPAQhAAAIQgAAEMkTAL1zNDlf/I1q9VvNarDIC8fRTmWyt7EApI4P9NVtNGYF0yFYjTe9/dEqg69+PH2HVgzW7aE2NWHNAlikjYEe2VjafDKWeYSEAAQhAAAIQgAAEHEoA2erQxBAWBCAAAQhAAAIQgAAEIAABCEAAAhCAAAQg4C4CyFZ35YtoIQABCEAAAhCAAAQgAAEIQAACEIAABCAAAYcSQLY6NDGEBQEIQAACEIAABCAAAQhAAAIQgAAEIAABCLiLALLVXfkiWghAAAIQgAAEIAABCEAAAhCAAAQgAAEIQMChBJCtDk0MYUEAAhCAAAQgAAEIQAACEIAABCAAAQhAAALuIoBsdVe+iBYCEIAABCAAAQhAAAIQgAAEIAABCEAAAhBwKAFkq0MTQ1gQgAAEIAABCEAAAhCAAAQgAAEIQAACEICAuwggW92VL6KFAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEHEoA2erQxBAWBCAAAQhAAAIQgAAEIAABCEAAAhCAAAQg4C4CyFZ35YtoIQABCEAAAhCAAAQgAAEIQAACEIAABCAAAYcSQLY6NDGEBQEIQAACEIAABCAAAQhAAAIQgAAEIAABCLiLALLVXfkiWghAAAIQgAAEIAABCEAAAhCAAAQgAAEIQMChBJCtDk0MYUEAAhCAAAQgAAEIQAACEIAABCAAAQhAAALuIoBsdVe+iBYCEIAABCAAAQhAAAIQgAAEIAABCEAAAhBwKAFkq0MTQ1gQgAAEIAABCEAAAhCAAAQgAAEIQAACEICAuwggW92VL6KFAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEHEoA2erQxBAWBCAAAQhAAAIQgAAEIAABCEAAAhCAAAQg4C4CyFZ35YtoIQABCEAAAhCAAAQgAAEIQAACEIAABCAAAYcSQLY6NDGEBQEIQAACEIAABCAAAQhAAAIQgAAEIAABCLiLALLVXfkiWghAAAIQgAAEIAABCEAAAhCAAAQgAAEIQMChBJCtDk0MYUEAAhCAAAQgAAEIQAACEIAABCAAAQhAAALuIoBsdVe+iBYCEIAABCAAAQhAAAIQgAAEIAABCEAAAhBwKAFkq0MTQ1gQgAAEIAABCEAAAhCAAAQgAAEIQAACEICAuwggW92VL6KFAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEHEoA2erQxBAWBCAAAQhAAAIQgAAEIAABCEAAAhCAAAQg4C4CyFZ35YtoIQABCEAAAhCAAAQgAAEIQAACEIAABCAAAYcSQLY6NDGEBQEIQAACEIAABCAAAQhAAAIQgAAEIAABCLiLALLVXfkiWghAAAIQgAAEIAABCEAAAhCAAAQgAAEIQMChBJCtDk0MYUEAAhCAAAQgAAEIQAACEIAABCAAAQhAAALuIoBsdVe+iBYCEIAABCAAAQhAAAIQgAAEIAABCEAAAhBwKAFkq0MTQ1gQgAAEIAABCEAAAhCAAAQgAAEIQAACEICAuwggW92VL6KFAAQgAAEIQAACEIAABCAAAQhAAAIQgAAEHEoA2erQxBAWBCAAAQhAAAIQgAAEIAABCEAAAhCAAAQg4C4CyFZ35YtoIQABCEAAAhCAAAQgAAEIQAACEIAABCAAAYcSQLY6NDGEBQEIQAACEIAABCAAAQhAAAIQgAAEIAABCLiLALLVXfkiWghAAAIQgAAEIAABCEAAAhCAAAQgAAEIQMChBJCtDk0MYUEAAhCAAAQgAAEIQAACEIAABCAAAQhAAALuIoBsdVe+iBYCEIAABCAAAQhAAAIQgAAEIAABCEAAAhBwKAFkq0MTQ1gQgAAEIAABCEAAAhCAAAQgAAEIQAACEICAuwj8/5SqEqF4COLwAAAAAElFTkSuQmCC", "text/html": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "dynamics.plot_history(colors=['red', 'green'], show_intervals=True)" ] }, { "cell_type": "code", "execution_count": 8, "id": "23c4b3ba", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "0: A <-> B\n", "Final 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": 8, "metadata": {}, "output_type": "execute_result" } ], "source": [ "# Verify that the reaction has reached equilibrium\n", "dynamics.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.8.10" } }, "nbformat": 4, "nbformat_minor": 5 }