{ "cells": [ { "cell_type": "markdown", "id": "76d0cfef-8918-4004-ba88-d27425d05675", "metadata": {}, "source": [ "### One-bin `A <-> 3B` reaction, with 1st-order kinetics in both directions,\n", "### taken to equilibrium\n", "\n", "Diffusion not applicable (just 1 bin)\n", "\n", "LAST REVISED: June 4, 2023\n", "\n", "* [First Step](#reaction_2_sec_2_first_step)\n", "* [Numerous more steps](#reaction_2_sec_2)\n", "* [Equilibrium](#reaction_2_sec_2_equilibrium)" ] }, { "cell_type": "code", "execution_count": 1, "id": "38fd1aeb-070b-47bc-b346-471da0396a42", "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": "3b64aa9a-b253-4926-805d-61238df949ea", "metadata": {}, "outputs": [], "source": [ "from experiments.get_notebook_info import get_notebook_basename\n", "\n", "from src.modules.reactions.reaction_data import ChemData as chem\n", "from src.modules.reactions.reaction_dynamics import ReactionDynamics\n", "from src.life_1D.bio_sim_1d import BioSim1D\n", "\n", "import plotly.express as px\n", "from src.modules.html_log.html_log import HtmlLog as log\n", "from src.modules.visualization.graphic_log import GraphicLog" ] }, { "cell_type": "code", "execution_count": 3, "id": "cdc446c6-4b9f-4938-b119-8855f4c530fe", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "-> Output will be LOGGED into the file 'reaction_2.log.htm'\n" ] } ], "source": [ "# Initialize the HTML logging\n", "log_file = get_notebook_basename() + \".log.htm\" # Use the notebook base filename for the log file\n", "\n", "# Set up the use of some specified graphic (Vue) components\n", "GraphicLog.config(filename=log_file,\n", " components=[\"vue_cytoscape_1\"],\n", " extra_js=\"https://cdnjs.cloudflare.com/ajax/libs/cytoscape/3.21.2/cytoscape.umd.js\")" ] }, { "cell_type": "code", "execution_count": 4, "id": "55ae0052-2f52-4d9b-8f13-cb65adb9d04b", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "SYSTEM STATE at Time t = 0:\n", "1 bins and 2 species:\n", " Species 0 (A). Diff rate: None. Conc: [10.]\n", " Species 1 (B). Diff rate: None. Conc: [50.]\n" ] } ], "source": [ "# Initialize the system\n", "chem_data = chem(names=[\"A\", \"B\"]) # NOTE: Diffusion not applicable (just 1 bin)\n", "\n", "# Reaction A <-> 3B , with 1st-order kinetics in both directions\n", "chem_data.add_reaction(reactants=[\"A\"], products=[(3,\"B\")], forward_rate=5., reverse_rate=2.)\n", "\n", "bio = BioSim1D(n_bins=1, chem_data=chem_data)\n", "\n", "bio.set_uniform_concentration(species_index=0, conc=10.)\n", "bio.set_uniform_concentration(species_index=1, conc=50.)\n", "\n", "bio.describe_state()" ] }, { "cell_type": "code", "execution_count": 5, "id": "025dc8d4-cc77-48fe-90f6-eab7f3383310", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Number of reactions: 1 (at temp. 25 C)\n", "0: A <-> 3 B (kF = 5 / kR = 2 / Delta_G = -2,271.45 / K = 2.5) | 1st order in all reactants & products\n" ] } ], "source": [ "chem_data.describe_reactions()" ] }, { "cell_type": "code", "execution_count": 6, "id": "fd1b6f4c-1327-4f5c-8ac0-ad0a61cc728e", "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", "
SYSTEM TIMEABcaption
0010.050.0Initial state
\n", "
" ], "text/plain": [ " SYSTEM TIME A B caption\n", "0 0 10.0 50.0 Initial state" ] }, "execution_count": 6, "metadata": {}, "output_type": "execute_result" } ], "source": [ "# Save the state of the concentrations of all species at bin 0\n", "bio.add_snapshot(bio.bin_snapshot(bin_address = 0), caption=\"Initial state\")\n", "bio.get_history()" ] }, { "cell_type": "code", "execution_count": 7, "id": "6aa601cb-5f4a-4a16-bde8-abb514bb25c2", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "[GRAPHIC ELEMENT SENT TO LOG FILE `reaction_2.log.htm`]\n" ] } ], "source": [ "# Send the plot to the HTML log file\n", "graph_data = chem_data.prepare_graph_network()\n", "GraphicLog.export_plot(graph_data, \"vue_cytoscape_1\")" ] }, { "cell_type": "markdown", "id": "c505025d-9bfd-485f-8465-2204353575a7", "metadata": { "tags": [] }, "source": [ "### First step" ] }, { "cell_type": "code", "execution_count": 8, "id": "1ec91d5e-717e-45f5-a88c-f8e8df328e48", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "SYSTEM STATE at Time t = 0.1:\n", "1 bins and 2 species:\n", " Species 0 (A). Diff rate: None. Conc: [15.]\n", " Species 1 (B). Diff rate: None. Conc: [35.]\n" ] } ], "source": [ "# First step\n", "bio.react(time_step=0.1, n_steps=1, snapshots={\"sample_bin\": 0})\n", "bio.describe_state()" ] }, { "cell_type": "markdown", "id": "7a9b3abf-86e7-4ef3-9ddc-23e4a2698c7c", "metadata": {}, "source": [ "_Early in the reaction :_ \n", "[A] = 15. [B] = 35." ] }, { "cell_type": "code", "execution_count": 9, "id": "23c2fb5d-c5ea-4873-a4fa-cd28a55865b0", "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", "
SYSTEM TIMEABcaption
00.010.050.0Initial state
10.115.035.0
\n", "
" ], "text/plain": [ " SYSTEM TIME A B caption\n", "0 0.0 10.0 50.0 Initial state\n", "1 0.1 15.0 35.0 " ] }, "execution_count": 9, "metadata": {}, "output_type": "execute_result" } ], "source": [ "bio.get_history()" ] }, { "cell_type": "markdown", "id": "8eaf0873-53fd-4928-804c-2a2385fc4931", "metadata": {}, "source": [ "### Numerous more steps" ] }, { "cell_type": "code", "execution_count": 10, "id": "6c178528-e88b-4be4-ab61-a5f1f6d57c36", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "SYSTEM STATE at Time t = 1.1:\n", "1 bins and 2 species:\n", " Species 0 (A). Diff rate: None. Conc: [14.54545455]\n", " Species 1 (B). Diff rate: None. Conc: [36.36363636]\n" ] } ], "source": [ "# Numerous more steps\n", "bio.react(time_step=0.1, n_steps=10, snapshots={\"sample_bin\": 0})\n", "\n", "bio.describe_state()" ] }, { "cell_type": "markdown", "id": "260cfdf1-a424-44cb-af29-7cdc693d16f6", "metadata": { "tags": [] }, "source": [ "### Equilibrium" ] }, { "cell_type": "markdown", "id": "17457676-5007-41b3-b8c2-a83a3d72ee47", "metadata": {}, "source": [ "Consistent with the 5/2 ratio of forward/reverse rates (and the 1st order reactions),\n", "the systems settles in the equilibrium: [A] = 14.54545455 , [B] = 36.36363636" ] }, { "cell_type": "code", "execution_count": 11, "id": "bb44d88c-bb34-41d7-b8c2-ca15938495d6", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "A <-> 3 B\n", "Final concentrations: [B] = 36.36 ; [A] = 14.55\n", "1. Ratio of reactant/product concentrations, adjusted for reaction orders: 2.5\n", " Formula used: [B] / [A]\n", "2. Ratio of forward/reverse reaction rates: 2.5\n", "Discrepancy between the two values: 6.875e-10 %\n", "Reaction IS in equilibrium (within 1% tolerance)\n", "\n" ] }, { "data": { "text/plain": [ "True" ] }, "execution_count": 11, "metadata": {}, "output_type": "execute_result" } ], "source": [ "# Verify that the reaction has reached equilibrium\n", "bio.reaction_dynamics.is_in_equilibrium(conc=bio.bin_snapshot(bin_address = 0))" ] }, { "cell_type": "code", "execution_count": 12, "id": "0f7ccafe-3fa4-4934-8f6f-8f6b774d1cf6", "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", "
SYSTEM TIMEABcaption
00.010.00000050.000000Initial state
10.115.00000035.000000
20.214.50000036.500000
30.314.55000036.350000
40.414.54500036.365000
50.514.54550036.363500
60.614.54545036.363650
70.714.54545536.363635
80.814.54545436.363637
90.914.54545536.363636
101.014.54545536.363636
111.114.54545536.363636
\n", "
" ], "text/plain": [ " SYSTEM TIME A B caption\n", "0 0.0 10.000000 50.000000 Initial state\n", "1 0.1 15.000000 35.000000 \n", "2 0.2 14.500000 36.500000 \n", "3 0.3 14.550000 36.350000 \n", "4 0.4 14.545000 36.365000 \n", "5 0.5 14.545500 36.363500 \n", "6 0.6 14.545450 36.363650 \n", "7 0.7 14.545455 36.363635 \n", "8 0.8 14.545454 36.363637 \n", "9 0.9 14.545455 36.363636 \n", "10 1.0 14.545455 36.363636 \n", "11 1.1 14.545455 36.363636 " ] }, "execution_count": 12, "metadata": {}, "output_type": "execute_result" } ], "source": [ "bio.get_history()" ] }, { "cell_type": "markdown", "id": "7fbc0ddf-9903-48c2-af88-5fa2d2048371", "metadata": {}, "source": [ "Note how the simulation initially **OVERSHOT** the equilibrium values; the first step was too large!" ] }, { "cell_type": "markdown", "id": "715af433-0129-4fdf-a172-77f9d3d127a2", "metadata": { "tags": [] }, "source": [ "# Plots of changes of concentration with time" ] }, { "cell_type": "code", "execution_count": 13, "id": "9ee92b39-26ee-4fa0-a93b-a0ee5048293d", "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": "navy", "dash": "solid" }, "marker": { "symbol": "circle" }, "mode": "lines", "name": "A", "orientation": "v", "showlegend": true, "type": "scatter", "x": [ 0, 0.1, 0.2, 0.30000000000000004, 0.4, 0.5, 0.6, 0.7, 0.7999999999999999, 0.8999999999999999, 0.9999999999999999, 1.0999999999999999 ], "xaxis": "x", "y": [ 10, 15, 14.5, 14.55, 14.545000000000002, 14.5455, 14.54545, 14.545455, 14.5454545, 14.54545455, 14.545454545, 14.545454545499998 ], "yaxis": "y" }, { "hovertemplate": "Chemical=B
SYSTEM TIME=%{x}
concentration=%{y}", "legendgroup": "B", "line": { "color": "darkorange", "dash": "solid" }, "marker": { "symbol": "circle" }, "mode": "lines", "name": "B", "orientation": "v", "showlegend": true, "type": "scatter", "x": [ 0, 0.1, 0.2, 0.30000000000000004, 0.4, 0.5, 0.6, 0.7, 0.7999999999999999, 0.8999999999999999, 0.9999999999999999, 1.0999999999999999 ], "xaxis": "x", "y": [ 50, 35, 36.5, 36.35, 36.365, 36.3635, 36.36365, 36.363635, 36.363636500000005, 36.36363635, 36.363636365, 36.3636363635 ], "yaxis": "y" } ], "layout": { "autosize": true, "legend": { "title": { "text": "Chemical" }, "tracegroupgap": 0 }, "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": "Changes in concentrations with time" }, "xaxis": { "anchor": "y", "autorange": true, "domain": [ 0, 1 ], "range": [ 0, 1.0999999999999999 ], "title": { "text": "SYSTEM TIME" }, "type": "linear" }, "yaxis": { "anchor": "x", "autorange": true, "domain": [ 0, 1 ], "range": [ 7.777777777777778, 52.22222222222222 ], "title": { "text": "concentration" }, "type": "linear" } } }, "image/png": "iVBORw0KGgoAAAANSUhEUgAABTkAAAFoCAYAAACc6v8hAAAAAXNSR0IArs4c6QAAIABJREFUeF7t3Q2UVeWd5/v/qfeSKgrEANKjKBgDKIqxiXRmxeXE7js3GCdZrhUaJ65rmgxNYye5KgNNYdOG2FqkGJTbpuUSlsTuiUsG13Bv2ivpl5hhyHRCQhuNqGASTDAdBXxBqKKo93PXs6v28Zxd530/Z+/n2c/3rJWFVJ397P/z+T8nBT+evXcqnU6nhRcCCCCAAAIIIIAAAggggAACCCCAAAIIIGCpQIqQ09LOUTYCCCCAAAIIIIAAAggggAACCCCAAAIIeAKEnCwEBBBAAAEEEEAAAQQQQAABBBBAAAEEELBagJDT6vZRPAIIIIAAAggggAACCCCAAAIIIIAAAggQcrIGEEAAAQQQQAABBBBAAAEEEEAAAQQQQMBqAUJOq9tH8QgggAACCCCAAAIIIIAAAggggAACCCBAyMkaQAABBBBAAAEEEEAAAQQQQAABBBBAAAGrBQg5rW4fxSOAAAIIIIAAAggggAACCCCAAAIIIIAAISdrAAEEEEAAAQQQQAABBBBAAAEEEEAAAQSsFiDktLp9FI8AAggggAACCCCAAAIIIIAAAggggAAChJysAQQQQAABBBBAAAEEEEAAAQQQQAABBBCwWoCQ0+r2UTwCCCCAAAIIIIAAAggggAACCCCAAAIIEHKyBhBAAAEEEEAAAQQQQAABBBBAAAEEEEDAagFCTqvbR/EIIIAAAggggAACCCCAAAIIIIAAAgggQMjJGkAAAQQQQAABBBBAAAEEEEAAAQQQQAABqwUIOa1uH8UjgAACCCCAAAIIIIAAAggggAACCCCAACEnawABBBBAAAEEEEAAAQQQQAABBBBAAAEErBYg5LS6fRSPAAIIIIAAAggggAACCCCAAAIIIIAAAoScrAEEEEAAAQQQQAABBBBAAAEEEEAAAQQQsFqAkNPq9lE8AggggAACCCCAAAIIIIAAAggggAACCBBysgYQQAABBBBAAAEEEEAAAQQQQAABBBBAwGoBQk6r20fxCCCAAAIIIIAAAggggAACCCCAAAIIIEDIyRpAAAEEEEAAAQQQQAABBBBAAAEEEEAAAasFCDmtbh/FI4AAAggggAACCCCAAAIIIIAAAggggAAhJ2sAAQQQQAABBBBAAAEEEEAAAQQQQAABBKwWIOS0un0UjwACCCCAAAIIIIAAAggggAACCCCAAAKEnKwBBBBAAAEEEEAAAQQQQAABBBBAAAEEELBagJDT6vZRPAIIIIAAAggggAACCCCAAAIIIIAAAggQcrIGEEAAAQQQQAABBBBAAAEEEEAAAQQQQMBqAUJOq9tH8QgggAACCCCAAAIIIIAAAggggAACCCBAyMkaQAABBBBAAAEEEEAAAQQQQAABBBBAAAGrBQg5rW4fxSOAAAIIIIAAAggggAACCCCAAAIIIIAAISdrAAEEEEAAAQQQQAABBBBAAAEEEEAAAQSsFiDktLp9FI8AAggggAACCCCAAAIIIIAAAggggAAChJysAQQQQAABBBBAAAEEEEAAAQQQQAABBBCwWoCQ0+r2UTwCCCCAAAIIIIAAAggggAACCCCAAAIIEHKyBhBAAAEEEEAAAQQQQAABBBBAAAEEEEDAagFCzpDte3jHHvnJi0dl++Z7ZGpHe8jROBwBBOIUOH2mR1avf0Q+tmie3LtqWZylcG4EEEAAAQQQQAABBBBAAAEEEKhAgJCzAJYfdhw+8nrOO754+9Kc8IOQ8wOevfsOyMbuXfLEtvWyeNG8CpYhb/UFlOFjf/Md2dG9RubOnqUdptj4ai3v+/6Pa3Zu7ZMpc8B867LQ51ZHyHns+Juyat1WuevOz8htS2/MqZLPSJlN420IIIAAAggggAACCCCAAAIIVChAyJkHzA8igoGmH4CoQ/ydm4SchJwVfuaKvp2QU6fm2FiEnPpNGREBBBBAAAEEEEAAAQQQQAAB0wQIOQMdOfTiUfnC3ZvlgXUrJuzC8t/65N5/kqU3L/EuTyfkNG1J211PnCGn3XKVVR/XTs7KquTdCCCAAAIIIIAAAggggAACCCBQrgAhZ5bU+f5BuX/LLnnjzVNl32PTD0se6lwpG7p2in95+y03L5FNa1dIa0uTdwZ/7GefO5jTm+D71DfLHdMfyN+p5v9+4bzL5Z3TZydcLpuvhnxhrh/0ZhdaziXo+QK6SudSaOHmqylYuzrX40/tywxRaCeuut/iJ264xguz/VehUDtoq94ftAi+J9jT7MuX1fHqkn7/lT1WvnNln8+3fOSrfyqPfPNpUWvp4hnTvMvLf/bKL3PGVcf53/Mvey93/OD9ZXW7VrO+/LWr5pX9ucq35oKXiwffE5xPttWFU9oz9+Qsd41kr1n/3G+dfDdnKftrMViLvzt82a03ee/PXhv+mgzWm2+tlvvZLvcHA+9DAAEEEEAAAQQQQAABBBBAwDYBQs6sjvkBxdJP3lD2Q0f8AGLh/DmZYDTfff1UCLFl+275/G2/n7nXYqHgptwxVen57qOY756A+eaWr04/gAqGb+XcJ7JQyKmCx1I+xT44+eaoav/W7u/K6js/6x2qwunsAMyf26WzpmdCsez7rGYHoPnmXCjwVo7P/eB5+eM7bs0E18FQPLhLMDv4yg6o8nmVumemsgyGl2re6jj1yr4HZDXjZz9EK9/6DOsadn3teWZ/5nOWHexlu6pzdHbtzNxbtFj4Hgx0K1kjhdZsqXtyZn+Wss8XXBt+4FlqzZT72bbthxP1IoAAAggggAACCCCAAAIIIFCJACFnnpAz3wNDCqEWuuxVBSvZgUyh44OBjB9c5ntie3BMP9zo6lyZ86Cf7N1hfuil6jxx6r2cXXB+OJYduhSaz8uv/UpaW5qLPgynkjCpEh+147LYTtJCwWAwbCr0UBk/LJs5/cJMuF3OZeP5eqdMg+ctFHrl61OpkDPfuii0tsKOXwvXMOsruN7V7+/r2ikLPnK59Pb2ZdZ2cK1Xsi4rWSM6Q061kzM7oM7XO3W+MJ/tSn4w8F4EEEAAAQQQQAABBBBAAAEEbBMg5KxhyJlv92O+S4aDO/OKBafZYxYKoYJBSKHAJDuQ84NSv758l9GXWtyVhEnlhIh+CFsqLC4U4AYDq1IBljqffyl0oTGzDQq9Jxialgo51eXz965a5g0dJuQsdKl0qZ2A/pyC664WrmHWV3Adq7EOPv+q/B/L/r385SN/Kw92rpRZMy7ydvWWCqwrvSdnoV3X+T4T1ezkrDbkrOSzXerzy/cRQAABBBBAAAEEEEAAAQQQsFmAkDNPyFnp5eqFdl3muyz1nffOZC6jVaeudCdnObsug8FHofAre+GWujdksQcx+ePUIuQs9WCnYuFT8Hv9AwOZ+y36oaKqPfg+9bXg5e/BD3mhe6xmv8+/JD6KkNO/xUF2n8Ls5Gxpbi5oUK2rf3/afEF/OetL2frBa+dXPi9df/WkLLl+gXzqk0u8WtV//7t/e53X4zWrlmV2N1eyLisJwgv9H3+UIWeln22bf1hROwIIIIAAAggggAACCCCAAALFBAg5s3TCPHgoeG+/YLBSaIdemJCz0CXflezkLPXxyA7zSj18qJIwKck7OYOmtQ45S42fvUuwkp2itdjJWSwsLrW+1LH+mlch5tb/+795uzfVg5X8XZ2f/oPfk6898rc5/5BQybq0LeQstpOz1Geb7yOAAAIIIIAAAggggAACCCCQJAFCzkA3/QejFNtZduDgz7wH6UztaM88Cb1UyFloR2KYkDPfQ1zUdILBRyWX2v7D/p/IjUsWZZ4Kr8YrtjMtm6+SMKnckLPQHNV5/Vq/+/2Dku/WAJXek1ON6V+uXqg+ZXng4Ivy72/6WNFLy7NdSoWQ2ZerF7rPpxqv2BrKd9/SfAFYJePruidntmuY9ZW9FlWw2dE+KdOvQvfn9IPR4Poo9Q8E2T1RY1TyGSoWPAZNK7n3ZtjPdpJ+cDEXBBBAAAEEEEAAAQQQQAABBIIChJx51oR/OW32E7izAwb1336oWe79M/OFdf7Xqr0nZ77gpdATp/2g7aNXfzjn4UPBp77nm0+5gWQtQk4/3Nv3/R/n7M7L3mWo3qMuV85+ynmxp4CXE2Dl29UbvNem//ufvvyLnNr8YE39qh4mU0nIWSxQruQ+koWe2l3J+PkMwrqGWV/ZYeOzzx2U7H+IKLTuC4Wcpf6BoJw1UujHSbFAVHfImR38lvps8+MPAQQQQAABBBBAAAEEEEAAgSQLEHIW6G52SJT9lmDwWW7IqcbwgxV/PDXWJ264Rjq7dk4I8Mq5z2cw9FG/V4GpeojQ1h17JPgwk0L3kcyeU773qF2rwZ2q+dhqFXL6QdXG7l2Z0+aryb8vZbZv9r03q7kUOThmvgcy5bvHZHZ9lYSc+daJfxl3sfuTBterOv9DnStlQ9fOCesguA5Lja/TNcz68vuq6gmG3tlrJHjZe6GQPrtv/j80XDilvaz7tvr3Fy30wyHYD/8zVouQM9//DxT6DCT5hxlzQwABBBBAAAEEEEAAAQQQcFuAkDOB/fdDNf+J6QmcIlNCAAEEEEAAAQQQQAABBBBAAAEEEEAgI0DIafliULva1G7QxYvmZWZS6IExlk+V8hFAAAEEEEAAAQQQQAABBBBAAAEEEMgrQMhp+cIIXkqsphO8pN7yKVI+AggggAACCCCAAAIIIIAAAggggAACRQUIOVkgCCCAAAIIIIAAAggggAACCCCAAAIIIGC1ACGn1e2jeAQQQAABBBBAAAEEEEAAAQQQQAABBBAg5GQNIIAAAggggAACCCCAAAIIIIAAAggggIDVAoScVreP4hFAAAEEEEAAAQQQQAABBBBAAAEEEECAkJM1gAACCCCAAAIIIIAAAggggAACCCCAAAJWCxByWt0+ikcAAQQQQAABBBBAAAEEEEAAAQQQQAABQk7WAAIIIIAAAggggAACCCCAAAIIIIAAAghYLUDIaXX7KB4BBBBAAAEEEEAAAQQQQAABBBBAAAEECDlZAwgggAACCCCAAAIIIIAAAggggAACCCBgtQAhp9Xto3gEEEAAAQQQQAABBBBAAAEEEEAAAQQQIORkDSCAAAIIIIAAAggggAACCCCAAAIIIICA1QKEnFa3j+IRQAABBBBAAAEEEEAAAQQQQAABBBBAgJCTNYAAAggggAACCCCAAAIIIIAAAggggAACVgsQclrdPopHAAEEEEAAAQQQQAABBBBAAAEEEEAAAUJO1gACCCCAAAIIIIAAAggggAACCCCAAAIIWC1AyGl1+ygeAQQQQAABBBBAAAEEEEAAAQQQQAABBAg5WQMIIIAAAggggAACCCCAAAIIIIAAAgggYLUAIafV7aN4BBBAAAEEEEAAAQQQQAABBBBAAAEEECDkZA0ggAACCCCAAAIIIIAAAggggAACCCCAgNUChJxWt4/iEUAAAQQQQAABBBBAAAEEEEAAAQQQQICQkzWAAAIIIIAAAggggAACCCCAAAIIIIAAAlYLEHJa3T6KRwABBBBAAAEEEEAAAQQQQAABBBBAAAFCTtYAAggggAACCCCAAAIIIIAAAggggAACCFgtQMhpdfsoHgEEEEAAAQQQQAABBBBAAAEEEEAAAQQIOVkDCCCAAAIIIIAAAggggAACCCCAAAIIIGC1ACGn1e2jeAQQQAABBBBAAAEEEEAAAQQQQAABBBAg5GQNIIAAAggggAACCCCAAAIIIIAAAggggIDVAoScVreP4hFAAAEEEEAAAQQQQAABBBBAAAEEEECAkJM1gAACCCCAAAIIIIAAAggggAACCCCAAAJWCxByWt0+ikcAAQQQQAABBBBAAAEEEEAAAQQQQAABQk7WAAIIIIAAAggggAACCCCAAAIIIIAAAghYLUDIaXX7KB4BBBBAAAEEEEAAAQQQQAABBBBAAAEECDlZAwgggAACCCCAAAIIIIAAAggggAACCCBgtQAhp9Xto3gEEEAAAQQQQAABBBBAAAEEEEAAAQQQIORkDSCAAAIIIIAAAggggAACCCCAAAIIIICA1QKEnFa3j+IRQAABBBBAAAEEEEAAAQQQQAABBBBAgJCTNYAAAggggAACCCCAAAIIIIAAAggggAACVgsQclrdPopHAAEEEEAAAQQQQAABBBBAAAEEEEAAAUJO1gACCCCAAAIIIIAAAggggAACCCCAAAIIWC1AyGl1+ygeAQQQQAABBBBAAAEEEEAAAQQQQAABBAg5WQMIIIAAAggggAACCCCAAAIIIIAAAgggYLUAIafV7aN4BBBAAAEEEEAAAQQQQAABBBBAAAEEECDkZA0ggAACCCCAAAIIIIAAAggggAACCCCAgNUChJxWt4/iEUAAAQQQQAABBBBAAAEEEEAAAQQQQICQkzWAAAIIIIAAAggggAACCCCAAAIIIIAAAlYLEHJa3T6KRwABBBBAAAEEEEAAAQQQQAABBBBAAAFCTtYAAggggAACCCCAAAIIIIAAAggggAACCFgtQMipoX1vvntewygMgYCbAhe2N0nfwIj0D464CcCsEdAgMGNqi7xzZkBGRtMaRmMIBNwTSKVEZk5tlbfe48907nWfGesSaGqok8mTGr2fR7wQQKA6gUktDdJQn5Iz54aqG4CjIhGYNa01kvNwksoFCDkrN5twBCGnBkSGcFaAkNPZ1jNxjQKEnBoxGcpJAUJOJ9vOpDULEHJqBmU4JwUIOe1oOyGnuX0i5NTQG0JODYgM4awAIaezrWfiGgUIOTViMpSTAoScTradSWsWIOTUDMpwTgoQctrRdkJOc/tEyKmhN4ScGhAZwlkBQk5nW8/ENQoQcmrEZCgnBQg5nWw7k9YsQMipGZThnBQg5LSj7YSc5vaJkFNEHt6xRx5/al9Olx5Yt0JuW3qj97W9+w7Ixu5d3n/fcvMS2bR2hbS2NGXeT8hp7gKnMvMFCDnN7xEVmi9AyGl+j6jQbAFCTrP7Q3V2CBBy2tEnqjRbgJDT7P741RFymtsnQs7xkFO16N5VyyZ06tCLR2Xrjj2yffM9MrWj3QtEg+8l5DR3gVOZ+QKEnOb3iArNFyDkNL9HVGi2ACGn2f2hOjsECDnt6BNVmi1AyGl2fwg5ze8PIWeJkFOFmpddMjOzqzMYeqoWE3Kav9Cp0FwBQk5ze0Nl9ggQctrTKyo1U4CQ08y+UJVdAoScdvWLas0UIOQ0sy/BqtjJaW6fCDnzXK7uX6p+vn9Q7t+yS5ZcvyATch47/qbc17VTHuxcKXNnz/I6S8hp7gKnMvMFCDnN7xEVmi9AyGl+j6jQbAFCTrP7Q3V2CBBy2tEnqjRbgJDT7P741ZkUcp4+0yOr1z8ia1Ytk8WL5tkBWKBKdavIg8+/OuEWkZVMipAzoKVCzFXrtkpX50q5et4cL+T83K03ZRbLhJDzR1+Tc7M/K6NTF1TiznsRQGBcoLWpXoZGRmV4JI1JhQLpdFpS6m/mvJwXaGtpkL6BYRnlY+T8WgCgOgH1/6RtrQ3Sc364ugE4CgEEpL4uJc2NddI3MFJzDf4MVHNiThCTgPrHgrqUSP/QaEwVcNpyBNpbG8p5m7b35HuOzBdvX+rdcjHOkNM/97Jbb8psDAwzaULOMHpFjvUvUf/UJ5eU3sm5NSWSqpfBBV+UwRu+JunmC2tUFcMikEyB1uZ6GRpOy/AIP8gr7XBaUpISUq1K3ZL4/kmtDdLXPyLqL328EECgCoGUSFtLo/SeH6riYA5BAAEl4IWcTfXS11/7fyzgz0CsuaQKNKqQsy4lA4O1/8eCpBpGMa/2CxqjOE0mwLx01vSc3Y3+VcdqQ94Vl/8OOzmzusFOzjxLM/s+nCXvyXlgvcihr3ujjDZ1SO+iDdK74MuRLHhOgkASBLhcPQldZA5xC3C5etwd4Py2C3C5uu0dpH4TBLhc3YQuUIPtAlyubkcHo7pcXeVRJ069V/TybX835V13fkb+v3/6kTz73EEP0d/pmS2avSP04hnTZEf3Gu82jH5oes2CufLSq8dyxvij5Z/yQtTDR173hvJv76j+O98uUv9r/vuzzxPckbpw/pzMQ77VeOzk1LD+VQP2PXdQPn/bH3ijBS9HL+fp6qeOvyqTf7xOWn7zrDfGcMeVcmbJNhm4+CYNFTIEAskWIORMdn+ZXTQChJzROHOW5AoQcia3t8wsOgFCzuisOVNyBQg57ehtFCFnuZeC++97570zmdAy+zaM/n06VcCoXuoSd/VSWVdn107vmFkzLvKuYv7py7+YMIZ6rx+GBvOyYMiZr2Z1HvVSdXzz28/IzZ+4PvN8m2CIS8ipYf37ibWfdqshn9i2PueGrQp6Y/cu72y33LxkQoruP3io+a390nHwbmk483Pvvf3/ZqmcXbJFhtsu11ApQyCQTAFCzmT2lVlFK0DIGa03Z0ueACFn8nrKjKIXIOSM3pwzJk+AkNOOnkYRcuYLKvPp5NtNGXyIthpry2O7pWvDSpna0e4Nk/2efLdqzPcg7uDXgueuNKQM1lXp8fk8uFxdw2co5+nq6RGZ9No3pf2nD0jd4PuSrmuScwv+VHoWbZB0wyQNZ2MIBJIlQMiZrH4ym3gECDnjceesyREg5ExOL5lJfAKEnPHZc+bkCBBy2tFL20JOtZvyC3dvzourLj/XFXIGb/cYPGG+TYbZl7MTchqy/nNCzvGa6gZOS/sLm2TSa4+LpEdkpGWG9Pzu16TvijtEhKchG9I6yjBAgJDTgCZQgvUChJzWt5AJxCxAyBlzAzh9IgQIORPRRiYRswAhZ8wNKPP0UYSclV6uvmbVsswVycEdl8HbMBYKH5dcvyDzlPRqdnIWCzn9nalLP3lD5pL54OXvhJxlLsBavy1fyOmfs/H9IzL5R3dL88kfeF8amnadd7/OwQ8trnVZjI+AFQKEnFa0iSINFyDkNLxBlGe8ACGn8S2iQAsECDktaBIlGi9AyGl8i7wCowg51XmKPXjIv9dlvqer57tc/b6unfJg58rM/TCzpcsJNNX7w1yurup9+pn9Obd/JOQ0dL0XCzn9kluP/51M/sl6qT/3a28nZ9+cZdKzuEtGWmcaOivKQiAaAULOaJw5S7IFCDmT3V9mV3sBQs7aG3OG5AsQcia/x8yw9gKEnLU31nGGqEJOfzfnpbOm54SD/nNj1PNkygk5/XDyjTdPTXia+SWzpsvV8+Z4Dx4Ku5PT362pnvR+29IbPWo/jFX/7T/oKPuJ7tkPO2Inp47VqWGMckJOdZrU6IC0vfyotL3ULanhXknXT5Lea/+z9F59t6TrmjVUwhAI2CdAyGlfz6jYPAFCTvN6QkV2CRBy2tUvqjVTgJDTzL5QlV0ChJx29CuqkNPXUDs6H39qXwYn+z6W5Tx4qNA4C+fP8ULPluZmLSGnOo8fdL518l3vtMF7bvoP9VZfX7v6D+Vbu7+b2WFKyGnI+i835PTLret/WyYf6pQLjj0lImkZaZstZ3+3S85f9llDZkQZCEQnQMgZnTVnSq4AIWdye8vMohEg5IzGmbMkW4CQM9n9ZXbRCBByRuMc9ixRh5xh63XpeJ6urqHblYac/ikb331BpvzwS6J+Va/BD/2enPm335ChKfM1VMUQCNghQMhpR5+o0mwBQk6z+0N15gsQcprfIyo0X4CQ0/weUaH5AoSc5vdIVUjIaW6fCDk19KbakHPs1Gm54Nhuaf+XP5f682+JpOrl3JUrpOejX5XR5qkaqmMIBMwWIOQ0uz9UZ4cAIacdfaJKcwUIOc3tDZXZI0DIaU+vqNRcAUJOc3uTXRkhp7l9IuTU0JtwIedYAamRPml78evS9spfeffuHG2aIr3X/bn0zlvlBZ+8EEiqACFnUjvLvKIUIOSMUptzJVGAkDOJXWVOUQsQckYtzvmSKEDIaUdXCTnN7RMhp4be6Ag5/TLqe9+QyYf+TFqPf8f70nDHlXJmyTYZuPgmDZUyBALmCRBymtcTKrJPgJDTvp5RsVkChJxm9YNq7BQg5LSzb1RtlgAhp1n9KFQNIae5fSLk1NAbnSGnX07zyR/I5B/dLY3vH/G+1H/JLXL2hm4ZbrtcQ8UMgYA5AoSc5vSCSuwVIOS0t3dUboYAIacZfaAKuwUIOe3uH9WbIUDIaUYfSlVByFlKKL7vE3JqsK9FyOmVlR6VSa89Lu0vbJK6gfckXdck5xb8qfQs2iDphkkaKmcIBOIXIOSMvwdUYL8AIaf9PWQG8QoQcsbrz9mTIUDImYw+Mot4BQg54/Uv9+yEnOVKRf8+Qk4N5jULOcdrqxs8K20//Zq0vfZNkfSwjLTMkJ7FD0jf3M+ru3lqmAFDIBCfACFnfPacOTkChJzJ6SUziUeAkDMed86aLAFCzmT1k9nEI0DIGY97pWcl5KxULLr3E3JqsK51yOmX2HDm59Jx8G5pfmu/96Whadd59+sc/NBiDbNgCATiESDkjMedsyZLgJAzWf1kNtELEHJGb84ZkydAyJm8njKj6AUIOaM3r+aMhJzVqEVzDCGnBueoQk6/1Jbf7JPJP/kzaeg55u3k7Jvzh9Kz+CEZaZ2pYTYMgUC0AoSc0XpztmQKEHIms6/MKjoBQs7orDlTcgUIOZPbW2YWnQAhZ3TWYc5EyBlGr7bHEnJq8I065FQlp0aHZNKrfy3tLz4kqeFeSddPkt5r10rv1f+npOuaNcyKIRCIRoCQMxpnzpJsAULOZPeX2dVegJCz9sacIfkChJzJ7zEzrL0AIWftjXWcgZBT5Hz/oNy/ZZe88eYp2b75Hpna0a6DNvQYhJyhCUXiCDn9suv635bJ/7JRLvjlf1VPKpKRSZfJ2Y91yfnZn9EwM4ZAoPYChJy1N+YMyRcg5Ex+j5lhbQUIOWvry+huCBByutFnZllbAULO2vrqGp2QU+TY8Tdlx9/+nZzt7ZMv3r5UFi+ap4s31DiEnKH4xg6OM+T0y288fVg6/vlL0vTOIe9LAzM+IWd/b5sMTZnPVdBJAAAgAElEQVSvYYYMgUDtBAg5a2fLyO4IEHK602tmWhsBQs7auDKqWwKEnG71m9nWRoCQszauukcl5BTZu+9AhvXXvzkh965appu5qvEIOatiyz3IhJDTr6j1V3tk8qH7pL7vtyKpejl35Qrp+ehXZbR5qoaZMgQC+gUIOfWbMqJ7AoSc7vWcGesVIOTU68lobgoQcrrZd2atV4CQU69nrUaLOuT8x388JoODI7WaTsFxP/3pK/N+T12qvmX7bvn8bb/vfX/LY7ula8NKIy5ZJ+TUsExMCjnVdFIjfdL20lZpe3mbpEbOy2jTFOn56EY595E/9oJPXgiYJEDIaVI3qMVWAUJOWztH3aYIEHKa0gnqsFmAkNPm7lG7KQKEnKZ0ongdUYec06dvkbff7osc5+2318pFF10w4bzqUvUn935P1q5e7n1P3ZtzyfUL5LalN0ZeY/CEhJwaWmBayOlPqf7cb2XyT9ZL6/H/7n1puONKOfPxR71L2XkhYIoAIacpnaAOmwUIOW3uHrWbIEDIaUIXqMF2AUJO2ztI/SYIEHKa0IXSNUQdct5xx145e3agdGGa3/Htb98mkydPfLD1wzv2yGWXzMyEmurS9YPPvyqb1q6Q1pYmzVVUNhwhZ2Veed9tasjpF9t06kfS8aO7Rd23U736L/m0nL3h6zLcdrmG2TMEAuEECDnD+XE0AkqAkJN1gEA4AULOcH4cjYASIORkHSAQXoCQM7xhFCNEHXJGMadyz3H6TI+sXv+IHD7yes4hF8+YJju618jc2bPKHaom7yPk1MBqesg5NsW0TPr5E9L+/P1SN/COpOuapHfBl6R3UaekGyZpUGAIBKoTIOSszo2jEMgWIORkPSAQToCQM5wfRyNAyMkaQECPACGnHsdaj+JyyHnoxaOydcce2b75npx7cAZ3d9a6B4XGJ+TUIG9HyDk20dRQr7S/8Jcy6eh2SY0OyUjLDOn53b+Uviv+o/quBg2GQKAyAULOyrx4NwL5BAg5WRcIhBMg5Aznx9EIEHKyBhDQI0DIqcex1qO4HHKqMFO9gk9TV+Hn08/sj/2SdUJODavfppDTn27D2V/K5B+vlZbf/oP3paFp18n7H/+G9ysvBKIUIOSMUptzJVWAkDOpnWVeUQkQckYlzXmSLMDl6knuLnOLSoCQMyrpcOdxOeQMJ1f7owk5NRjbGHL6025+a790/PAr0tDzS28nZ9+cP5SexQ/JSOtMDTIMgUBpAULO0ka8A4FSAoScpYT4PgLFBQg5WSEIhBcg5AxvyAgIEHLasQYIOc3tEyGnht7YHHJ60x8dlraj26XthQelbuispOsnSe+iddJ71VckXTfxSVoayBgCgYwAISeLAYHwAoSc4Q0ZwW0BQk63+8/s9QgQcupxZBS3BQg57eg/Iae5fSLk1NAb60POcYO6gXe9BxOpBxSJjMrIpMvk7Me65Pzsz2hQYggE8gsQcrIyEAgvQMgZ3pAR3BYg5HS7/8xejwAhpx5HRnFbgJDTjv4TcprbJ0JODb1JSsjpUzS+f0Q6/vlL0vT2j7wvDcz4hJz5+KMy3HGlBi2GQCBXgJCTFYFAeAFCzvCGjOC2ACGn2/1n9noECDn1ODKK2wKEnHb0n5DT3D4RcmroTdJCTp+k9fhemXxog9T3viGSqpdzH/mi9Fx3v4w2T9WgxhAIjAkQcrISEAgvQMgZ3pAR3BYg5HS7/8xejwAhpx5HRnFbgJDTjv4Tcprbp1hDztNnemT1+kfk8JHXJwgtnD9Htm++R6Z2tJurN15ZUkNONb3USL+0vbxN2l76L5Ia6ZPRpinSc91fyLl5K73gkxcCYQUIOcMKcjwCIoScrAIEwgkQcobz42gElAAhJ+sAgfAChJzhDaMYgZAzCuXqzhFryPnwjj1e1feuWlZd9YYcleSQ0yeuP39CJv9kvbT+aqxn6tJ1dQm7upSdFwJhBAg5w+hxLAJjAoScrAQEwgkQcobz42gECDlZAwjoESDk1ONY61EIOWstXP34sYWcahdn50M7Ze1dy2Xu7FnVz0Djkef7B+X+Lbu8ETetXSGtLU3ef+/dd0A2do99/Zabl+R8T33NhZDTZ256+5B0/PBL0nj6sPel/ks+LWdv+LoMt12usRMM5ZIAIadL3WautRIg5KyVLOO6IkDI6UqnmWctBdjJWUtdxnZFgJDTjk4TcprbJ0LO8d74Aeezzx3MCTIPvXhUtu7Yk7l0Pt/uU5dCzjGutFzwy29L+7/8hdT3n5R0XbP0XvUl6b12vaQbJpm72qnMSAFCTiPbQlGWCRByWtYwyjVOgJDTuJZQkIUChJwWNo2SjRMg5DSuJXkLIuQ0t0+xhZyKRAWGl10yU25bemPsQn4tqpCDz7+a2a0ZrDEYeqr3uxdyjrUrNXxO2l98SCa9+teSGh2UkdaZ0nP9A9J3xX9U3429pxRghwAhpx19okqzBQg5ze4P1ZkvQMhpfo+o0HwBQk7ze0SF5gsQcprfI1WhyyFn9gZBv1sXz5gmO7rXGHGVdqwh57Hjb8qTe78na1cvz1waHseSzt6dqS5N90NOVYu6fH3J9QsyQayq+b6unfJg58pMA10NOf1eNfT+Sib/eJ20/OZZ70tD066T9z/+De9XXgiUEiDkLCXE9xEoLUDIWdqIdyBQTICQk/WBQHgBQs7whoyAACGnHWuAkDM3J8vO0fzbPsbVydhCzmJPVlcYUT1dXTXj1785kXn4Ub6Q83O33iSLF83zepQv5Iyrecad91//p8g//YnIe0fHdnIuuEPkxq+LTLrYuFIpCIEkCAwOj3pPMuWFAAIIIIAAAgi4JMCfgVzqNnNFAAGTBPydnNmbAfNd8RxXzbGFnHFNOHhetYvz8af2TShHPWBo/ZfvkM2PfpudnJU0Kz0ik177prT/9AGpG3xf0g1t0nvNOum9+svevTt5IRAUYCcnawKB8ALs5AxvyAhuC7CT0+3+M3s9Auzk1OPIKG4LsJPTjv5HvpPz1/8oMjoYPc6cT084Z76QM7h5MPpCPzij8yFnED+4zZZ7cla3POsGTkv7C5tk0muPi6RHvKev9yzukvOz/0N1A3JUYgUIORPbWiYWoQAhZ4TYnCqRAoSciWwrk4pYgJAzYnBOl0gBQk472hp5yPnYdJHzb0ePc9fbIq0X5Zw33z051Ru+ePvSzBXS0RdqUMiptrV+4e7NOQZPbFufuTw8apxgyMnT1cN1oOHMz6Xjh1+W5pM/8AYamPEJOfPxR2W448pwA3N0YgQIORPTSiYSowAhZ4z4nDoRAoSciWgjk4hZgJAz5gZw+kQIEHLa0cbIQ859d4gMnI0e55ZvizRNzhtyZl+unm93Z/TFjp0x1p2c+a7bV/e8XLVuq9x152dieep6vhumqq9t7N7lganL2DetXZHzoCTXHzxUzuJtPf530n6oU9RDiiRVL+c+8p+k57q/kNHmqeUcznsSLEDImeDmMrXIBAg5I6PmRAkVIORMaGOZVqQChJyRcnOyhAoQctrR2MhDToNYCgWaplyyHlvI6cNkP9TH75sKP59+Zv+EMNGgvuaUQshZXmdSowPS9vKj0vZSt6SGe72As2fRRjk3b6UXfPJyU4CQ082+M2u9AoScej0ZzT0BQk73es6M9QsQcuo3ZUT3BAg57eg5IWfu09XZySki6unqnQ/tlLV3LZe5s2flrGS1m3PLY7ula8NKmdrRbvwqJ+SsrEV1/W/L5EOdcsGxp0Qk7V26ri5hV5ey83JPgJDTvZ4zY/0ChJz6TRnRLQFCTrf6zWxrI0DIWRtXRnVLgJDTjn4Tcu6SZ587mNOsB9atiOVq7OCKYSenhs8QIWd1iI3vviBTfvglUb+qV/+lt8qZj22RkbZLqxuQo6wUIOS0sm0UbZgAIadhDaEc6wQIOa1rGQUbKEDIaWBTKMk6AUJOO1rmcshpeodiCzkVjLpmf88z+2X75nsyOzbjvidnNQ0j5KxGzT8mLRcc2y3t//LnUn/+LUnXNUvvVV+W3mv/TNINk8IMzLGWCBByWtIoyjRagJDT6PZQnAUChJwWNIkSjRcg5DS+RRRogQAhpwVNEhFCTnP7FGvIqVhMe7p6Na0i5KxGLfeY1EiftL34dWl75VFJjfbLSOtM6bn+L6XvitvV87HCn4ARjBUg5DS2NRRmkQAhp0XNolQjBQg5jWwLRVkmQMhpWcMo10gBQk4j2zKhKEJOc/sUe8hpLk35lRFylm9V6p31vW/I5EN/Jq3Hv+O9dWjadfL+x7/h/cormQKEnMnsK7OKVoCQM1pvzpY8AULO5PWUGUUvQMgZvTlnTJ4AIacdPSXkNLdPhJwaekPIqQExMETzyR/I5B/dLY3vH/F2cvbNvV3OLu6S0ZYP6T8ZI8YqQMgZKz8nT4gAIWdCGsk0YhMg5IyNnhMnSICQM0HNZCqxCRByxkZf0YkJOSviivTNhJwauAk5NSDmGyI9KpNee1zaX9gkdQPvSbqhTXqvWSe9V3/Zu3cnr2QIEHImo4/MIl4BQs54/Tm7/QKEnPb3kBnEL0DIGX8PqMB+AUJOO3pIyGlunyIPOU+f6ZHV6x+RP/rD/12+9d/+Xg4feT2vzsL5c3IeSGQuoQghZ227Uzd4Rtp/+jWZdHS7d6KRSZfJmRu6pf/ST9f2xIxeMwHV09TgGakbOisdjedlYGBQBodHx883dg/WtLoXq/pbp/ca/1X9Pp313/73st7nHecdknVsvnGC71Pjjr8vM0awhsztYbPPkXu+nLr9WgvNJaeG8XkXqzurxhyT8XHS3rF+PX778jhKStL1LTXrLwNHL0DIGb05Z0yWACFnsvrJbOIRIOSMx52zJkuAkNOOfhJymtunyENOn0KFnZ0P7ZS1dy2XubNn5QiphxE9/cx+2bR2hbS2NJmrN14ZIWc0LWp4/zXp+PG90vzW//BOODDzRjmzZJsMT5kXTQGcRVIj56Vu8Kykhs6O/Tp8TuoGTo/9fuispNTXBs5kvl83dEZSw+elbuBdSQ31jB/Ti6RDAum68TA1J2gOBtEfhNK5IXWewDrzILLsQLdEMF5O8J1OyVhI7CXk+QP2TG488X05wXh2jeUE7Cp0Dr4vXw3euBMDdnVsU0O9DI2MSjrt0OJiqghoFFCfrKbGOhkY8v/BTePgDIWAIwJ1KZH6+pQMDfPDyJGWM80aCNTXpbx9F8MjfI5qwKttyOY7Dmgbi4H0ChgZch47/qZseWy3dG1YKVM72vXOuAajEXLWALXIkC3/+vcy+eB/lobesV3A5z7yn6Tno5tktHlqtIVYdrbM7snhnrFdlFlhpfpvGd9ZmQkwvVBS7bhUAeYZ75YBvBBAAAEEEEAAAQQQQAABBBBwWmANIbSp/Tcy5Ny774AcfP5VdnKaumoMqavtlf9L2l94UFLDvV7A2bNoo5yb/yeGVKevDLV7cmwXpAoce7wdk95/D3/wNRVa+sGlt6NyuK9muydHmy+S0aYOSTe2y2jTZEk3Th7/tWPs697Xxr/X1CGjje2Sbhz/Vf2+aXIODvfk1LdWTBkpNdI/Xor64T/+B4DMFsN05kv+f6S87+V73/jXUvm/P3acemUfW2Qc9T5vLHXI2K8p9bWi4/iqE9+Xc6wat8J6vOP9+vMdm1Nrnroz8xaZ2t4kZ84Nyiib0Ez5GFCHZQJq14z6efTu2UHLKqdcBMwRaKhPibrU9sy5IXOKohIELBNoaaqX+jqRc/0jllXuVrnTrvrf3JqwRbONPORUuzRXrdsqb518tyDTxTOmyY7uNRMuYzfVlZ2c8XVGXQY9+dB9csEvvy0iozLccaWc+fijMjDjE/EVlTlzWuoGez64tFvthhzK3kXZI6mB98e+pr6XfRm4t3ty7LLw1KiOv3DVyWhjm6SzQ8fsgLKpXdJNU2S0MSu0bJo89nv1q/e/qSKpeu2uhJzaSRnQQQHuyelg05myVgHuyamVk8EcFeCenI42nmlrFeCenFo5azYY9+SsGW3ogSMPOf2Ki92TM/SsIh6AkDNi8Dynazx9WDr++UvS9M4h77vnL/0PcvZj3TLSdmlVxaVGB8Yu01Y7JP17SY7fh9K7z6QKJgfeH7+Uezy4HA8lVVg5FlCqe0+G38aerm8d3xn5we7JzE7KxrGvpZvVDsqx3ZLe9xrGf/V/39hWlUMUBxFyRqHMOZIuQMiZ9A4zv1oLEHLWWpjxXRAg5HShy8yx1gKEnLUW1jM+Iacex1qMElvIWYvJxDUmIWdc8hPP2/qrp2XyoQ1S3/db75s9C9fJ0MzfG7+Ue2xXpXfJ9+D747sn/Z2W/qXg6t6ThXcZVzrTsZ2QKnxUv6qdkm0fXN6tdlL63/Mv+W7ukNEGtatyfPdk80WVntK69xNyWtcyCjZQgJDTwKZQklUChJxWtYtiDRUg5DS0MZRllQAhpx3tIuQ0t0+EnBp6Q8ipAVHjEOp+gG0v/Rdp/9lDVY9acvekd+9JFUaO755snibqGP+elGOXfZu7e7JqmBocSMhZA1SGdE6AkNO5ljNhzQKEnJpBGc5JAUJOJ9vOpDULEHJqBq3RcIScNYLVMGysIWex+3MunD9Htm++h6era2iyq0PUn/tXaf/pV6Wu/70PdkY2ThZRuyVz7j3pPxRn7B6Uo83TXCWLZd6EnLGwc9KECRByJqyhTCdyAULOyMk5YQIFCDkT2FSmFLkAIWfk5FWdkJCzKrZIDoot5DzfPyj3b9klS65fINdedYU8ufd7snb1cmltaZKHd+yRT9xwjSxeNC8ShLAnYSdnWEGOd1mAkNPl7jN3XQKEnLokGcdVAUJOVzvPvHUKEHLq1GQsVwUIOe3oPCGnuX2KLeTMfvCQ4tny2G7p2rDS27l56MWj8vQz+2XT2hVe6Gn6i5DT9A5Rn8kChJwmd4fabBEg5LSlU9RpqgAhp6mdoS6bBAg5beoWtZoqQMhpamdy6yLkNLdPRoScF05pl66/elI6v/J5L+RUl7Fnh57m8o1VRshpeoeoz2QBQk6Tu0NttggQctrSKeo0VYCQ09TOUJdNAoScNnWLWk0VIOQ0tTOEnHZ0RiS2kDP7cvXblt7oXaJ+2SUzRf333n0H5ODzr7KT05ZVRJ0IhBAg5AyBx6EIjAsQcrIUEAgnQMgZzo+jEVAChJysAwTCCxByhjeMYgR2ckahXN05Ygs5g+Wqy9dXr39EDh95XS6eMU12dK+RubNnVTeriI9iJ2fE4JwuUQKEnIlqJ5OJSYCQMyZ4TpsYAULOxLSSicQoQMgZIz6nTowAIacdrSTkNLdPxoSc5hKVroyQs7QR70CgkAAhJ2sDgfAChJzhDRnBbQFCTrf7z+z1CBBy6nFkFLcFCDnt6D8hp7l9ii3kzH7wkC07Ngu1kZDT3AVOZeYLEHKa3yMqNF+AkNP8HlGh2QKEnGb3h+rsECDktKNPVGm2ACGn2f3xqyPkNLdPhJwaekPIqQGRIZwVIOR0tvVMXKMAIadGTIZyUoCQ08m2M2nNAoScmkEZzkkBQk472k7IaW6fYgs5FYl62NAnbrhGFi+aZ65QGZURcpaBxFsQKCBAyMnSQCC8ACFneENGcFuAkNPt/jN7PQKEnHocGcVtAUJOO/pPyGlun2INOY8df1Oe3Ps9Wbt6ubS2NJmrVKIyQk5rW0fhBggQchrQBEqwXoCQ0/oWMoGYBQg5Y24Ap0+EACFnItrIJGIWIOSMuQFlnp6Qs0yoGN4WW8iZ/TT1fPNeOH+ObN98j0ztaI+BpbJTEnJW5sW7EcgWIORkPSAQXoCQM7whI7gtQMjpdv+ZvR4BQk49jozitgAhpx39J+Q0t0+xhZzmklReGSFn5WYcgYAvQMjJWkAgvAAhZ3hDRnBbgJDT7f4zez0ChJx6HBnFbQFCTjv6T8hpbp9iCzmLPV390ItH5eln9sumtSusuIydkNPcBU5l5gsQcprfIyo0X4CQ0/weUaHZAoScZveH6uwQIOS0o09UabYAIafZ/fGrI+Q0t09GhpzqXp1bHtstXRtWcrm6uWuHyhDQIkDIqYWRQRwXIOR0fAEw/dAChJyhCRkAASHkZBEgEF6AkDO8YRQjEHJGoVzdOYwMOffuOyAHn381sp2c6invjz+1LyP4xLb1OU98V/Vs7N7lff+Wm5dMqIudnNUtPo5CQAkQcrIOEAgvQMgZ3pAR3BYg5HS7/8xejwAhpx5HRnFbgJDTjv4Tcprbp8hDTrVLc9W6rfLWyXcLqlw8Y5rs6F4jc2fPqrmcumz+W7u/K6vv/Kx3abyq776unfJg50rv/OrS+a079mQegqQCUfW6d9WyTG2EnDVvEydIsAAhZ4Kby9QiEyDkjIyaEyVUgJAzoY1lWpEKEHJGys3JEipAyGlHYwk5ze1T5CGnT1HsnpxxcvlPfV+zapm3m1OFmpddMlNuW3qjV1Yw9FRfI+SMs2Oc23YBQk7bO0j9JggQcprQBWqwWYCQ0+buUbspAoScpnSCOmwWIOS0o3uEnOb2KbaQ01QSFWJ2du30dpLOmnGR3L9llyy5fkEm5Azu9CTkNLWT1GWLACGnLZ2iTpMFCDlN7g612SBAyGlDl6jRdAFCTtM7RH02CBBy2tAlEUJOc/tEyDnem+zL6P17cp7vH/RCzs/delPmHp35Qs6e88PmdpjKEDBcoLWpXoZGRmV4JG14peaVl06nJaX+Zs7LeYG2lgbpGxiWUT5Gzq8FAKoTUP9P2tbaIPyZrjo/jkJACdTXpaS5sU76BkZqDsKfgWpOzAliElD/WFCXEukfGo2pAk5bjkB7a0M5b+M9MQjEGnL6l4YfPvL6hKkvnD8ncx/MKF2yL1e/et6csnZy9vQNRVki50IgUQKtzfUyNJyW4RF+kFfa2LSkJCWkWpW6JfH9k1obpK9/RNRf+nghgEAVAimRtpZG6T3Pn+mq0OMQBDwBL+Rsqpe+/tpvAOHPQCy6pAo0qpCzLiUDg7X/x4KkGkYxr/YLGqM4DeeoQiDWkDPfQ3yqmIP2Q7Lvw8k9ObXzMiACOQJcrs6CQCC8AJerhzdkBLcFuFzd7f4zez0CXK6ux5FR3BbgcnU7+s/l6ub2KbaQ05QHD6nLz5/7wfPyx3fc6nXJv2y9q3Old4k6T1c3d/FSWTIECDmT0UdmEa8AIWe8/pzdfgFCTvt7yAziFyDkjL8HVGC/ACGnHT0k5DS3T86HnP59N5997mCmS/49Of0v7N13QDZ27/J+e8vNS2TT2hXS2tKUeT9PVzd3gVOZ+QKEnOb3iArNFyDkNL9HVGi2ACGn2f2hOjsECDnt6BNVmi1AyGl2f/zqCDnN7VNsIaciCV4Kbi5T8coIOW3tHHWbIEDIaUIXqMF2AUJO2ztI/XELEHLG3QHOnwQBQs4kdJE5xC1AyBl3B8o7PyFneU5xvCvWkFNdGv7k3u/J2tXLc3ZGxgER5pyEnGH0ONZ1AUJO11cA89chQMipQ5ExXBYg5HS5+8xdlwAhpy5JxnFZgJDTju4Tcprbp9hCzmJPVldccT1dvZpWEXJWo8YxCIwJEHKyEhAIL0DIGd6QEdwWIOR0u//MXo8AIaceR0ZxW4CQ047+E3Ka26fYQk5zSSqvjJCzcjOOQMAXIORkLSAQXoCQM7whI7gtQMjpdv+ZvR4BQk49jozitgAhpx39J+Q0t0+EnBp6Q8ipAZEhnBUg5HS29UxcowAhp0ZMhnJSgJDTybYzac0ChJyaQRnOSQFCTjvaTshpbp9iDTmzn2x+8YxpsqN7jcyacZHcv2WXLLl+gdy29EZz5bIqI+S0ok0UaagAIaehjaEsqwQIOa1qF8UaKEDIaWBTKMk6AUJO61pGwQYKEHIa2JQ8JRFymtunWENO/+nqn/rkEtmyfbd8/rbfl7mzZ8mhF4/K08/sl01rV1jxQCJCTnMXOJWZL0DIaX6PqNB8AUJO83tEhWYLEHKa3R+qs0OAkNOOPlGl2QKEnGb3x6+OkNPcPsUWcqoHD3U+tFPW3rXc272ZHXKqp65veWy3dG1YKVM72s3VG6+MkNP4FlGgwQKEnAY3h9KsESDktKZVFGqoACGnoY2hLKsECDmtahfFGipAyGloYwJlEXKa2ycjQ052cpq7YKgMAd0ChJy6RRnPRQFCThe7zpx1ChBy6tRkLFcFCDld7Tzz1ilAyKlTs3ZjEXLWzjbsyLGFnKrwvfsOyMHnX5XOr3xeHt31/3iXq184pV1Wr39Elt16E/fkDNtdjkfAAgFCTguaRInGCxByGt8iCjRcgJDT8AZRnhUChJxWtIkiDRcg5DS8QePlEXKa26dYQ07FonZtfuHuzTlCT2xbL4sXzTNXLVAZl6tb0yoKNVCAkNPAplCSdQKEnNa1jIINEyDkNKwhlGOlACGnlW2jaMMECDkNa0iBcgg5ze1T7CGnuTTlV0bIWb4V70QgKEDIyZpAILwAIWd4Q0ZwW4CQ0+3+M3s9AoScehwZxW0BQk47+k/IaW6fYg051dPVT5x6L+cp6uf7B+X+LbtkyfULuFzd3HVDZQhoEyDk1EbJQA4LEHI63HymrkWAkFMLI4M4LkDI6fgCYPpaBAg5tTDWfBBCzpoTV32C2EJOP8z83K03Tbg0nQcPVd1PDkTAOgFCTutaRsEGChByGtgUSrJKgJDTqnZRrKEChJyGNoayrBIg5LSjXYSc5vYptpDz9Jke6Xxop6y9a7nMnT0rR+jY8Tdly2O7pWvDSpna0W6u3nhlXK5ufIso0GABQk6Dm0Np1ggQclrTKgo1VICQ09DGUJZVAoScVrWLYg0VIOQ0tDGBsgg5ze1TbCEnOznNXRRUhkCUAoScUWpzrqQKEHImtbPMKyoBQs6opDlPkgUIOZPcXeYWlQAhZ1TS4c5DyBnOr5ZHxxZyqkmpy9I7u3bKju41md2cahfnqnVb5a47P8M9OWvZecZGwBABQk5DGkEZVmwc23AAACAASURBVAsQclrdPoo3QICQ04AmUIL1AoSc1reQCRggQMhpQBPKKIGQswykmN4Sa8ip5uyHmm+dfDdD8MS29RPu0xmTT1mn5XL1sph4EwJ5BQg5WRgIhBcg5AxvyAhuCxByut1/Zq9HgJBTjyOjuC1AyGlH/wk5ze1T7CGnuTTlV0bIWb4V70QgKEDIyZpAILwAIWd4Q0ZwW4CQ0+3+M3s9AoScehwZxW0BQk47+k/IaW6fCDk19IaQUwMiQzgrQMjpbOuZuEYBQk6NmAzlpAAhp5NtZ9KaBQg5NYMynJMChJx2tJ2Q09w+xRpyqiesr17/iBw+8voEoYXz58j2zffwdHVz1w6VIaBFgJBTCyODOC5AyOn4AmD6oQUIOUMTMgACQsjJIkAgvAAhZ3jDKEYg5IxCubpzxBpyPrxjj1f1vauWVVe9IUexk9OQRlCGlQKEnFa2jaINEyDkNKwhlGOdACGndS2jYAMFCDkNbAolWSdAyGlHywg5ze1TbCGn2sXZ+dBOWXvX8syT1c1lKl4ZIaetnaNuEwQIOU3oAjXYLkDIaXsHqT9uAULOuDvA+ZMgQMiZhC4yh7gFCDnj7kB55yfkLM8pjncRcmpQJ+TUgMgQzgoQcjrbeiauUYCQUyMmQzkpQMjpZNuZtGYBQk7NoAznpAAhpx1tJ+Q0t0+xhZyKRF2uftklM+W2pTeaK1RGZYScZSDxFgQKCBBysjQQCC9AyBnekBHcFiDkdLv/zF6PACGnHkdGcVuAkNOO/hNymtunWEPOY8fflCf3fk/Wrl4urS1N5iqVqIyQ09rWUbgBAoScBjSBEqwXIOS0voVMIGYBQs6YG8DpEyFAyJmINjKJmAUIOWNuQJmnJ+QsEyqGt8UWchZ7srpy4OnqMawGTolADAKEnDGgc8rECRByJq6lTChiAULOiME5XSIFCDkT2VYmFbEAIWfE4FWejpCzSrgIDost5IxgbpGdgp2ckVFzogQKEHImsKlMKXIBQs7IyTlhwgQIORPWUKYTiwAhZyzsnDRhAoScdjSUkNPcPhFyaugNIacGRIZwVoCQ09nWM3GNAoScGjEZykkBQk4n286kNQsQcmoGZTgnBQg57Wg7Iae5fYo95Dz04lH5wt2bc4Se2LZeFi+aZ65aoDJCTmtaRaEGChByGtgUSrJOgJDTupZRsGEChJyGNYRyrBQg5LSybRRtmAAhp2ENKVAOIae5fYo15FQB59Yde2T75ntkake7p6QeRrRq3Va5687PWPPUdUJOcxc4lZkvQMhpfo+o0HwBQk7ze0SFZgsQcprdH6qzQ4CQ044+UaXZAoScZvfHr46Q09w+xRZynu8flPu37JLP3XrThF2bKvx8+pn9smntCiueuk7Iae4CpzLzBQg5ze8RFZovQMhpfo+o0GwBQk6z+0N1dggQctrRJ6o0W4CQ0+z+EHKa35/YQk71dPXOh3bK2ruWy9zZs3Kk1G7OLY/tlq4NKzM7PGtF6Yetzz53MHOK4OXye/cdkI3du7zv33LzkgnhKyFnrbrDuC4IEHK60GXmWGsBQs5aCzN+0gUIOZPeYeYXhQAhZxTKnCPpAoScdnSYnZzm9im2kNOUnZwqbP3W7u/K6js/6+0aVbtIO7t2yo7uNV74Gryk/uEde7xu3rtqWaarhJzmLnAqM1+AkNP8HlGh+QKEnOb3iArNFiDkNLs/VGeHACGnHX2iSrMFCDnN7o9fHSGnuX2KLeRUJGqH5J5n9ht1T04Veq5e/4isWbXMu4xehZqXXTIzc3/QfPcRJeQ0d4FTmfkChJzm94gKzRcg5DS/R1RotgAhp9n9oTo7BAg57egTVZotQMhpdn8IOc3vT6whp+Ix7enq6lL5+7p2yoOdK2XWjIu8+4YuuX5BJuTM/r5/mT0hp/kLnQrNFSDkNLc3VGaPACGnPb2iUjMFCDnN7AtV2SVAyGlXv6jWTAFCTjP7EqyKnZzm9in2kNMkGv8Sej/UzHdJfb6Q06Q5UAsCCLgjMDg8KuovFLwQQAABBBBAAAGXBPgzkEvdZq4IIIBA+QKxhpzqUvATp97LeZBPMGgsfyrh3umfd+b0CzP328xXCzs5wzlzNAJBAXZysiYQCC/ATs7whozgtgA7Od3uP7PXI8BOTj2OjOK2ADs57eg/OznN7VNsIacpDx5SrckXcPot456c5i5eKkuGACFnMvrILOIVIOSM15+z2y9AyGl/D5lB/AKEnPH3gArsFyDktKOHhJzm9im2kFM94KfzoZ2y9q7l3lPMs19qt+SWx3ZL14aVMrWjvaZ6pXaO8nT1mvIzOAJCyMkiQCC8ACFneENGcFuAkNPt/jN7PQKEnHocGcVtAUJOO/pPyGlun2ILOU3ZyakC1VXrtspbJ9/N6dIXb1+auWxdPQV+Y/cu7/u33Lwk5/J69TUePGTuAqcy8wUIOc3vERWaL0DIaX6PqNBsAUJOs/tDdXYIEHLa0SeqNFuAkNPs/vjVEXKa26fYQk5FonZJdnbtlB3dazK7Of3Q8a47P5N5orm5fGOVEXKa3iHqM1mAkNPk7lCbLQKEnLZ0ijpNFSDkNLUz1GWTACGnTd2iVlMFCDlN7UxuXYSc5vYp1pBTseTbSfnEtvWyeNE8c9UClRFyWtMqCjVQgJDTwKZQknUChJzWtYyCDRMg5DSsIZRjpQAhp5Vto2jDBAg5DWtIgXIIOc3tU+whp7k05VdGyFm+Fe9EIChAyMmaQCC8ACFneENGcFuAkNPt/jN7PQKEnHocGcVtAUJOO/pPyGlunwg5NfSGkFMDIkM4K0DI6WzrmbhGAUJOjZgM5aQAIaeTbWfSmgUIOTWDMpyTAoScdrSdkNPcPhFyaugNIacGRIZwVoCQ09nWM3GNAoScGjEZykkBQk4n286kNQsQcmoGZTgnBQg57Wg7Iae5fSLk1NAbQk4NiAzhrAAhp7OtZ+IaBQg5NWIylJMChJxOtp1JaxYg5NQMynBOChBy2tF2Qk5z+0TIqaE3hJwaEBnCWQFCTmdbz8Q1ChByasRkKCcFCDmdbDuT1ixAyKkZlOGcFCDktKPthJzm9omQU0NvCDk1IDKEswKEnM62nolrFCDk1IjJUE4KEHI62XYmrVmAkFMzKMM5KUDIaUfbCTnN7RMhp4beEHJqQGQIZwUIOZ1tPRPXKEDIqRGToZwUIOR0su1MWrMAIadmUIZzUoCQ0462E3Ka2ydCTg29IeTUgMgQzgoQcjrbeiauUYCQUyMmQzkpQMjpZNuZtGYBQk7NoAznpAAhpx1tJ+Q0t0+EnBp6Q8ipAZEhnBUg5HS29UxcowAhp0ZMhnJSgJDTybYzac0ChJyaQRnOSQFCTjvaTshpbp8IOTX0hpBTAyJDOCtAyOls65m4RgFCTo2YDOWkACGnk21n0poFCDk1gzKckwKEnHa0nZDT3D4RcmroDSGnBkSGcFaAkNPZ1jNxjQKEnBoxGcpJAUJOJ9vOpDULEHJqBmU4JwUIOe1oOyGnuX0i5NTQG0JODYgM4awAIaezrWfiGgUIOTViMpSTAoScTradSWsWIOTUDMpwTgoQctrRdkJOc/tEyKmhN4ScGhBrNMTLh0/JuXOD0thQJw2N9dKgfm2ok8ZG9fs6aWyoz/o15X2vubmhRtUwbD4BQk7WBQLhBQg5wxsygtsChJxu95/Z6xEg5NTjyChuCxBy2tF/Qk5z+0TIqaE3hJwaEEMO0d8/LK8cPiWHXzopL790Ql55+ZS89LOTIUcVueCCxvEQNE9I6oWlE0PShob6sRDV/9ULUwsc39SQ9d6U1NdPHDMTxDbWe99vaqqTenXu8TEzvzbmP761tTG0Qy0HIOSspS5juyJAyOlKp5lnrQQIOWsly7guCRByutRt5lorAULOWsnqHZeQU6+nztEIOTVoEnJqQKxgiJ6zg/LSz05kAs2XD5+UX/z8vQpGcPetapdqY2PK29UaDEkb1NezQtaxkFbtbh0PbVWgW6/C2tRYuOuHvFk7Y5tUEKve7wWw4ztms0JfbyetOn58B62q46IpzdI/OCqDQ6MTGpNOp/M2q8CXvfcWOmbse/l7X/SYQsulSBFVnadAbcXmVI1DNceI7tqKDKjTrrq1UOT/H3Q7VLO+C9hNbm2UnvNDMpqnxmocCn5Yqv0cVfPZ03hM0c9RsR8J1fRI4zHF/3+rcOHVfI4KfSyrWT/pKj7jkf5/3ejEnzeSEmlraZTe80Pu/iGBmSMQUqC+LiXNTfXS1z8cciQOR8BdAfX3qLq6lAwMjriLYMHMt3b/gQVVulkiIaeGvhNyakAsMMTbp86NhZmH1c7ME/LKSyfljTfO5H33lfOmydULZ8g11870fr36mhnS3t5UdXGDgyMyNDQiw0NpGRoekZHhURkaGpXhYfX19PivozI8NCJDw6MyMpwee7/3Pv9X9f2xY7yvD4+93x9zSJ1Dfc97z8Qx1df9Y9S5R0bGa/C+Pnbc2K/psXH97/vf88Yd+x8vBBBAAAEEEEAAAQQQQAABBBAIJ5BO3x9uAI6umQAhpwZaQk4NiCLyxvH3vUvOVZipLjV/5fBJefvtvryDX7topiy4evp4oDldrlo4nXtpltGGvr6hsTA2E4JmhbHj4awX5nphbHao64e742FugeNVmDo4MPxBGOsHvllhbDAwllGR4dFRybuxRl0/WOCVShXeVpcqcFyR4UTrMcV6obE2dZpCcyo0nyiPUTujCvcv/zeL9qjAgJH1teh8ijS9wHE6e9TaVC8DQ6OSbwdd0fMUK7uKtVpoQUbVoyjXd8HlUPT/t/KDV7MWxuZaxeeo4Hqs/PNazWe8aI+KDFjz/69LibS3spOzjD9K8BYECgqwk5PFgUB4AXZyhjeMYgR2ckahXN05CDmrc8s5ipCzMsTR0bT88hfvZQLNlw+f8HZqnj0zMGGgtrZGuerqGXKVv0Pzmuny4SuneZdK80qGAPfkTEYfmUW8AtyTM15/zm6/APfktL+HzCB+Ae7JGX8PqMB+Ae7JaUcPuSenuX0i5NTQG0LOwojqcuwjR97xdmWqHZoqzDzyyttyPs89r6ZNa/UuMVe7Mv1Lzi+7fGrBnWoaWscQBggQchrQBEqwXoCQ0/oWMoGYBQg5Y24Ap0+EACFnItrIJGIWIOSMuQFlnp6Qs0yoGN5GyKkBnZBzDFFdCq2eau7dQ3P8wUA///m73mXPwdes32mXhdeo+2bOHPt14XS5eFa7hm4whG0ChJy2dYx6TRQg5DSxK9RkkwAhp03dolZTBQg5Te0MddkkQMhpR7cIOc3tEyGnht64GHK+f7rfCzO9QNN7MNAJ+dXr74u6FD37pZ4Mp3ZjjgWaM7xf1f+mTG3RIM8QSRAg5ExCF5lD3AKEnHF3gPPbLkDIaXsHqd8EAUJOE7pADbYLEHLa0UFCTnP7RMipoTdJDzlPvNU7/oTzk3J4/JLz3/7r2QlyDY11cuWV08YCTe8J59O9p5xfcEGjBmWGSKoAIWdSO8u8ohQg5IxSm3MlUYCQM4ldZU5RCxByRi3O+ZIoQMhpR1cJOc3tEyGnht4kJeRMp0WO//p0ZofmK+M7Nd999/wEpdbWRpl/1YcyOzPVg4Hmz79IGpvqNYgyhEsChJwudZu51kqAkLNWsozrigAhpyudZp61FCDkrKUuY7siQMhpR6cJOc3tEyGnht7YGHKOjKTlFz9/17vU3H8g0Ksvn5SensEJIpM7mr1dmQu9+2fO9B4MdMWHLxR1KTovBMIKEHKGFeR4BEQIOVkFCIQTIOQM58fRCCgBQk7WAQLhBQg5wxtGMQIhZxTK1Z2DkLM6t5yjTA85BwaGvSeaqyebjwWaJ+Xoq++I+nrwNX3GJO8Sc//emVddM0MuvbRDgxJDIJBfgJCTlYFAeAFCzvCGjOC2ACGn2/1n9noECDn1ODKK2wKEnHb0n5DT3D4RcmrojUkhZ2/vkLxy+ETmknP1YKBjv3hPhocnPuH80tlT5Oprpss13v0zZ8jCa2fIRRddoEGEIRAoX4CQs3wr3olAIQFCTtYGAuEECDnD+XE0AkqAkJN1gEB4AULO8IZRjEDIGYVydecg5KzOLeeouEJOda9M7+nmP1NPOVc7NE9599RU99bMftXXp2Tuhy/M2qE503vSeXt7k4bZMwQC4QQIOcP5cTQCSoCQk3WAQDgBQs5wfhyNACEnawABPQKEnHocaz0KIWethasfn5CzervMkVGEnOpp5l6gOf4wIHXJuXrqefDV1FQv8xd8yLtv5tg9NGd4DwhqaWnQMFOGQEC/ACGnflNGdE+AkNO9njNjvQKEnHo9Gc1NAXZyutl3Zq1XgJBTr2etRiPkrJVs+HEJOcMbis6QU+3CfP2YesL5CXl5fHemCjffP90/odK2tka56uoZop5sri41V4Hmh6+cJg0NdRpmxRAIRCNAyBmNM2dJtgAhZ7L7y+xqL0DIWXtjzpB8AULO5PeYGdZegJCz9sY6zkDIqUOxNmMQcma5Prxjj1x2yUy5bemNOdp79x2Qjd27vK/dcvMS2bR2hbS2fHCpd7Uh5/DQqLz22jveZeaH1QOBXjopr77ytpw7N/EJ5xde2OoFmdk7NC+7fKqoP5TzQsBmAUJOm7tH7aYIEHKa0gnqsFWAkNPWzlG3SQKEnCZ1g1psFSDktKNzhJzm9omQU0SyQ8wH1q3ICTkPvXhUtu7YI9s33yNTO9pFBaHqde+qZZmulhNy9vcPyyuHT3lPNvcCzcOn5OjRd2RocGTC6pj1O+3jDwJSDwRSl53PkItntZu7iqgMgRAChJwh8DgUgXEBQk6WAgLhBAg5w/lxNAJKgJCTdYBAeAFCzvCGUYxAyBmFcnXnIOTMcsu3kzP4tWDoqQ4Phpw9ZwfHn26unnI+tkPz2C9Py+ho7hOB1B+oL5+jHgg0XRaOP+H8mmtnyJSpLdV1k6MQsFCAkNPCplGycQKEnMa1hIIsEyDktKxhlGukACGnkW2hKMsECDntaBghp7l9IuQsEnKe7x+U+7fskiXXL8js7jx2/E25r2unPNi5UubOniUnTvTK9/7n8fFQ86S88tJJeeONM3k7vuCqD3lPNR8LNKfLVVdPl0mTeMK5uR8PKotCgJAzCmXOkXQBQs6kd5j51VqAkLPWwozvggAhpwtdZo61FiDkrLWwnvEJOfU41mIUQs4yQs7P3XqTLF40z3tnMORMpTZN6EvrBY1y9dXTRe3KXHTdxXLttTNFBZzqyee8EEAgV6C1qV6GRkZleCR3pzNOpQXS6bSkuDFvaSgH3tHW0iB9A8MSuGDAgZkzRQT0CKhbnLe1NkjP+WE9AzIKAg4K1NelpLmxTvoGJt6OSzcHfwbSLcp4pgiofyyoS4n0D42aUhJ15BFob23AxVABQs4yQs5iOzk/8pFvyIyZk7wg89pFY/+78sppUqf+n4kXAgiUFGhtrpeh4bQMj/CDvCRW4A1pSUlKCIcrdUvi+ye1Nkhf/4iov/TxQgCBKgRSIm0tjdJ7fqiKgzkEAQSUgBdyNtVLX3/t/7GAPwOx5pIq0KhCzrqUDOR5dkdS52zjvNovaLSxbCdqJuQsEnKqb1VzT04nVg6TRECTAJera4JkGKcFuFzd6fYzeQ0CXK6uAZEhnBfgcnXnlwAAGgS4XF0DYgRDcLl6BMhVnoKQs0TIqevp6lX2h8MQSLwAIWfiW8wEIxAg5IwAmVMkWoCQM9HtZXIRCRByRgTNaRItQMhpR3sJOc3tEyGniOzdd0A2du/KdOniGdNkR/ca78FC6pX9/VtuXiKb1q6Q1pYPHhgUfLq6ue2mMgTMEyDkNK8nVGSfACGnfT2jYrMECDnN6gfV2ClAyGln36jaLAFCTrP6UagaQk5z+0TIqaE3hJwaEBnCWQFCTmdbz8Q1ChByasRkKCcFCDmdbDuT1ixAyKkZlOGcFCDktKPthJzm9omQU0NvCDk1IDKEswKEnM62nolrFCDk1IjJUE4KEHI62XYmrVmAkFMzKMM5KUDIaUfbCTnN7RMhp4beEHJqQGQIZwUIOZ1tPRPXKEDIqRGToZwUIOR0su1MWrMAIadmUIZzUoCQ0462E3Ka2ydCTg29IeTUgMgQzgoQcjrbeiauUYCQUyMmQzkpQMjpZNuZtGYBQk7NoAznpAAhpx1tJ+Q0t0+EnBp6Q8ipAZEhnBUg5HS29UxcowAhp0ZMhnJSgJDTybYzac0ChJyaQRnOSQFCTjvaTshpbp8IOc3tDZUhgAACCCCAAAIIIIAAAggggAACCCCAQBkChJxlIPEWBBBAAAEEEEAAAQQQQAABBBBAAAEEEDBXgJDT3N5QGQIIIIAAAggggAACCCCAAAIIIIAAAgiUIUDIWQYSb0EAAQQQQAABBBBAAAEEEEAAAQQQQAABcwUIOavszd59B2Rj9y7v6FtuXiKb1q6Q1pamKkfjMASSK3C+f1Du37JLnn3uoDfJB9atkNuW3lhwwg/v2COPP7Uv8/1S70+uHDND4AOB02d6ZPX6R+Twkde9Lz6xbb0sXjSvJJH/+VNv5OdUSS7ekHCBaj5Hh148Kl+4e7Mns3D+HNm++R6Z2tGecCmmh0BhgUo/E9l/Z7p4xjTZ0b1G5s6eBTECCBQQOHb8Tdny2G7p2rCSnzesEgSqECDkrAJN/XDfumNP5g+6KpRRr3tXLatiNA5BINkC2Z8P/y+Ya1YtyxvQqEBm+9/8v/JHyz/l/VBXP+RXrdsqXZ0rywp0ki3J7FwV8IPKJdcv8P6BQH0u7uvaKQ92riz6F8Xsf2DgH+NcXT3M2xeo5nMU/PMemgi4LhD8+aMCzIPPv1rwH9GCnyE+U66vIOZfTCD7H+L4RzXWCgLVCxByVmGnQpvLLpmZ2Y3GD+wqEDnECQH1w7rzoZ2y9q7lmTCmkn8UCP6l1Ak0JolAQCD4L/rlfi78n1VquGJ/CQUcARcEKv0c5fv55YITc0SgmIAKNX/9mxOZjR2l/tEtGIKWej/6CCAg3j9ms5OTlYBA9QKEnBXa5fvLJT+wK0Tk7c4I5PtslPpX/2ycUjs/nYFkok4L5PuHtFL/WJD9/Uo+c05DM/lEC1T6OfKvJHjr5LsZly/evpSrdhK9SphcKYHgz55Sf07zv3/prOnebs/vfv9gTkha6nx8HwEXBQg5Xew6c9YpQMhZoaYfcn7u1psyl88SclaIyNudEcj3Q7qSwKVUkOMMJBN1WkCFM08/sz/ncsBin43gTptKPnNOQzP5RAtU+jkKvt8Pa5bdelPR+0onGpHJOS8QvJqtVMipwNQxrx37jfyvnxwW7snp/BICoAwBQs4ykHgLAkUECDkrXB7s5KwQjLc7LRBmJ6f6Q/GJU+/xsBSnVxCTVwKV7kALPrzLV+S+nKwnlwUq/RzlC0X5BwOXVxBz9wNL9av/HIJSIWfwH93U56qzaycPH2I5IVBEgJCT5YFAOAFCzir8uCdnFWgc4qRAtffkJOB0crkw6QICld5LMDgMwQxLC4GJ9zgrdW/bQlciZN+PEFcEXBOo9J6c1ez8dM2U+SIQFCDkZE0gEE6AkLMKP56uXgUahzgrUOrp6sFAk0vUnV0qTLyAQKmnQpe6jJaQk6WFgEipz5F/D86uzpXe7Yj898+cfqG3a63UjjWMEXBBoNTT1dXPmz3P7Jftm++RqR3tEvw9OzldWCXMMawAIWdYQY53XYCQs8oVoH5ob+ze5R3NJYBVInKYEwL+XxSffe6gN98H1q3IuZ9ZdsjZPzAgq9c/IoePvJ5jw2fMiaXCJIsI+AGL/9l4Ytv6zH2hCTlZOgiUJ1DscxQMOdWIwfcHf36Vd1behUCyBFRQ+YW7N3uTWjh/TibQVL8Phprqa9m3UOGenMlaC8xGr0DwZ44anQfe6TVmNDcECDnd6DOzRAABBBBAAAEEEEAAAQQQQAABBBBAILEChJyJbS0TQwABBBBAAAEEEEAAAQQQQAABBBBAwA0BQk43+swsEUAAAQQQQAABBBBAAAEEEEAAAQQQSKwAIWdiW8vEEEAAAQQQQAABBBBAAAEEEEAAAQQQcEOAkNONPjNLBBBAAAEEEEAAAQQQQAABBBBAAAEEEitAyJnY1jIxBBBAAAEEEEAAAQQQQAABBBBAAAEE3BAg5HSjz8wSAQQQQAABBBBAAAEEEEAAAQQQQACBxAoQcia2tUwMAQQQQAABBBBAAAEEEEAAAQQQQAABNwQIOd3oM7NEAAEEEEAAAQQQQAABBBBAAAEEEEAgsQKEnIltLRNDAAEEEEAAAQQQQAABBBBAAAEEEEDADQFCTjf6zCwRQAABBBBAAAEEEEAAAQQQQAABBBBIrAAhZ2Jby8QQQAABBBBAAAEEEEAAAQQQQAABBBBwQ4CQ040+M0sEEEAAAQQQQAABBBBAAAEEEEAAAQQSK0DImdjWMjEEEEAAAQQQQAABBBBAAAEEEEAAAQTcECDkdKPPzBIBBBBAAAEEEEAAAQQQQAABBBBAAIHEChByJra1TAwBBBBAAAEEEEAAAQQQQAABBBBAAAE3BAg53egzs0QAAQQQQAABBBBAAAEEEEAAAQQQQCCxAoSciW0tE0MAAQQQQAABBBBAAAEEEEAAAQQQQMANAUJON/rMLBFAAAEEEEDAQIHTZ3pk9fpH5PCR13Oqe2DdCvnUJ5fI/Vt2eV/ftHaFtLY0Zd5z7PibsmrdVrnrzs/IbUtvlGLjqO8/vGOPPP7UvoICC+fPkYe/+qey7ZtPy7PPHZzwvltuXuLVoF6qJvWeJ7atl8WL5mXee75/sOD3/Dft3XdANnaPzSnf6+IZ06R7459I918/lTFRtW3ffI9M7WjPzEP5qHllv/w5+t/Lrid4Ln8+2aYGLg9KQgABBBBAAAEEEKhAgJCzAizeigACCCCAAAII6BIIBpX+uOrrT+79nqxdvVz6Bwa8EHTZrTflhHoq0FOvoUsfLgAACB5JREFUe1ctk3LGyQ7z/EB0zapleUPKmdMv9MbN98oODr94+9Kc9x168ah84e7N3mHBALTYWEuuXzAhsPTPE6zFDzKDIaVv8NbJdyUYchabj65eMg4CCCCAAAIIIIBA/AKEnPH3gAoQQAABBBBAwEEBtatxzzP7M7sUCxGo8LCza6fs6F4jc2fPEvX7rTv2ZI4rdxx/fB0h5xWX/4789PAvZO1dy72a/FDymgVz5Yk9fy9dnStzAlSdIWdvX7/09vbJ5269KXMOFX62TWqV7//zC5lAuFBQ6uBSY8oIIIAAAggggIATAoScTrSZSSKAAAIIIICAaQLB8LJYfSrEO3HqPbnnjz8n93z1r3N2dlYyjjqHjpBT7b789W9OeCX7u0m3PLZb1O5OFcjWMuRU57zskply8PlXvUvo1W7Xzod2eudW4a+/65WQ07QVTz0IIIAAAggggEBtBQg5a+vL6AgggAACCCCAQF6BfPeMzHevSXVw9uXYwUu1KxmnnJCznHtyqpDz2quukPu6dsqDnSvlO3//v7zgUX1N3Su01iHnHy3/lHcZv7rk/jdvnvICV/9rwZCz2Hy4JycfTgQQQAABBBBAIDkChJzJ6SUzQQABBBBAAAFLBbLvZ6mmELzfpfqauiz9sb/5Tuay9XxTLWccXTs5/Qca/eSFIzKlo126NqyU997viSTkVLtHvcv0/+5/eAwqaL1wSnvO/UvZyWnph4GyEUAAAQQQQACBKgUIOauE4zAEEEAAAQQQQKAWAoUuPw/ei7PUuQuNozPkDD70yP99rXdyqpDTn8fHFs3zLpn3f8/l6qVWBt9HAAEEEEAAAQSSKUDImcy+MisEEEAAAQQQMFzgwMGfycL5c2RqR3tOpSoo9C8DVw/18V+FQs5Kx9EZcqrantz7T7L05iXePKIMOdW5/2H/T+SKy/+N9/AjQk7DFzzlIYAAAggggAACNRYg5KwxMMMjgAACCCCAAAL5BNTl1hu7d8kT29ZnnhLuX2Kt3q8eqpN9z8hCIWel4+gOObPnFnXImX1uQk4+ZwgggAACCCCAgNsChJxu95/ZI4AAAggggECMAn5AmV1Cvvtxqu8Xu1y9knFKhZzlPnhI3ZMz+NIRcvr1HT7yuje82u26ffM93k5R9ZR59VKXpwdfhUJOHjwU4wLn1AgggAACCCCAQIQChJwRYnMqBBBAAAEEEEAAAQQQQAABBBBAAAEEENAvQMip35QREUAAAQQQQAABBBBAAAEEEEAAAQQQQCBCAULOCLE5FQIIIIAAAggggAACCCCAAAIIIIAAAgjoFyDk1G/KiAgggAACCCCAAAIIIIAAAggggAACCCAQoQAhZ4TYnAoBBBBAAAEEEEAAAQQQQAABBBBAAAEE9AsQcuo3ZUQEEEAAAQQQQAABBBBAAAEEEEAAAQQQiFCAkDNCbE6FAAIIIIAAAggggAACCCCAAAIIIIAAAvoFCDn1mzIiAggggAACCCCAAAIIIIAAAggggAACCEQoQMgZITanQgABBBBAAAEEEEAAAQQQQAABBBBAAAH9AoSc+k0ZEQEEEEAAAQQQQAABBBBAAAEEEEAAAQQiFCDkjBCbUyGAAAIIIIAAAggggAACCCCAAAIIIICAfgFCTv2mjIgAAggggAACCCCAAAIIIIAAAggggAACEQoQckaIzakQQAABBBBAAAEEEEAAAQQQQAABBBBAQL8AIad+U0ZEAAEEEEAAAQQQQAABBBBAAAEEEEAAgQgFCDkjxOZUCCCAAAIIIIAAAggggAACCCCAAAIIIKBfgJBTvykjIoAAAggggAACCCCAAAIIIIAAAggggECEAoScEWJzKgQQQAABBBBAAAEEEEAAAQQQQAABBBDQL0DIqd+UERFAAAEEEEAAAQQQQAABBBBAAAEEEEAgQgFCzgixORUCCCCAAAIIIIAAAggggAACCCCAAAII6Bcg5NRvyogIIIAAAggggAACCCCAAAIIIIAAAgggEKEAIWeE2JwKAQQQQAABBBBAAAEEEEAAAQQQQAABBPQLEHLqN2VEBBBAAAEEEEAAAQQQQAABBBBAAAEEEIhQgJAzQmxOhQACCCCAAAIIIIAAAggggAACCCCAAAL6BQg59ZsyIgIIIIAAAggggAACCCCAAAIIIIAAAghEKEDIGSE2p0IAAQQQQAABBBBAAAEEEEAAAQQQQAAB/QKEnPpNGREBBBBAAAEEEEAAAQQQQAABBBBAAAEEIhQg5IwQm1MhgAACCCCAAAIIIIAAAggggAACCCCAgH4BQk79poyIAAIIIIAAAggggAACCCCAAAIIIIAAAhEKEHJGiM2pEEAAAQQQQAABBBBAAAEEEEAAAQQQQEC/ACGnflNGRAABBBBAAAEEEEAAAQQQQAABBBBAAIEIBQg5I8TmVAgggAACCCCAAAIIIIAAAggggAACCCCgX4CQU78pIyKAAAIIIIAAAggggAACCCCAAAIIIIBAhAKEnBFicyoEEEAAAQQQQAABBBBAAAEEEEAAAQQQ0C9AyKnflBERQAABBBBAAAEEEEAAAQQQQAABBBBAIEIBQs4IsTkVAggggAACCCCAAAIIIIAAAggggAACCOgXIOTUb8qICCCAAAIIIIAAAggggAACCCCAAAIIIBChACFnhNicCgEEEEAAAQQQQAABBBBAAAEEEEAAAQT0CxBy6jdlRAQQQAABBBBAAAEEEEAAAQQQQAABBBCIUICQM0JsToUAAggggAACCCCAAAIIIIAAAggggAAC+gX+f7Im4QLvnyXvAAAAAElFTkSuQmCC", "text/html": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fig = px.line(data_frame=bio.get_history(), x=\"SYSTEM TIME\", y=[\"A\", \"B\"], \n", " title=\"Changes in concentrations with time\",\n", " color_discrete_sequence = ['navy', 'darkorange'],\n", " labels={\"value\":\"concentration\", \"variable\":\"Chemical\"})\n", "fig.show()" ] }, { "cell_type": "code", "execution_count": 14, "id": "6b496fc2-c001-40a8-83e1-13ae119b0aaa", "metadata": {}, "outputs": [ { "data": { "application/vnd.plotly.v1+json": { "config": { "plotlyServerURL": "https://plot.ly" }, "data": [ { "hovertemplate": "Chemical=A
SYSTEM TIME=%{x}
concentration=%{y}", "legendgroup": "A", "line": { "color": "navy", "dash": "solid", "shape": "spline" }, "marker": { "symbol": "circle" }, "mode": "lines", "name": "A", "orientation": "v", "showlegend": true, "type": "scatter", "x": [ 0, 0.1, 0.2, 0.30000000000000004, 0.4, 0.5, 0.6, 0.7, 0.7999999999999999, 0.8999999999999999, 0.9999999999999999, 1.0999999999999999 ], "xaxis": "x", "y": [ 10, 15, 14.5, 14.55, 14.545000000000002, 14.5455, 14.54545, 14.545455, 14.5454545, 14.54545455, 14.545454545, 14.545454545499998 ], "yaxis": "y" }, { "hovertemplate": "Chemical=B
SYSTEM TIME=%{x}
concentration=%{y}", "legendgroup": "B", "line": { "color": "darkorange", "dash": "solid", "shape": "spline" }, "marker": { "symbol": "circle" }, "mode": "lines", "name": "B", "orientation": "v", "showlegend": true, "type": "scatter", "x": [ 0, 0.1, 0.2, 0.30000000000000004, 0.4, 0.5, 0.6, 0.7, 0.7999999999999999, 0.8999999999999999, 0.9999999999999999, 1.0999999999999999 ], "xaxis": "x", "y": [ 50, 35, 36.5, 36.35, 36.365, 36.3635, 36.36365, 36.363635, 36.363636500000005, 36.36363635, 36.363636365, 36.3636363635 ], "yaxis": "y" } ], "layout": { "autosize": true, "legend": { "title": { "text": "Chemical" }, "tracegroupgap": 0 }, "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": "Changes in concentrations with time (smoothed)" }, "xaxis": { "anchor": "y", "autorange": true, "domain": [ 0, 1 ], "range": [ 0, 1.0999999999999999 ], "title": { "text": "SYSTEM TIME" }, "type": "linear" }, "yaxis": { "anchor": "x", "autorange": true, "domain": [ 0, 1 ], "range": [ 7.777777777777778, 52.22222222222222 ], "title": { "text": "concentration" }, "type": "linear" } } }, "image/png": "iVBORw0KGgoAAAANSUhEUgAABTkAAAFoCAYAAACc6v8hAAAAAXNSR0IArs4c6QAAIABJREFUeF7s3QuUVdWd7/v/rl1PqgooUF4KKBADRJSEoJWXx479CrYjGZ4rTTo50ZDQHMw9PTQcbMpc2/bmxKJhoPbp29KElpiceGLjvdzu9kqSTsyxyQs1GhQfJAYUH4Ao76Leu/YdcxVru2uxd+3X3GvNueZ3j8EAqvae8z8//7Up6ldzrZVIp9Np4YEAAggggAACCCCAAAIIIIAAAggggAACCFgqkCDktLRzlI0AAggggAACCCCAAAIIIIAAAggggAACngAhJwcCAggggAACCCCAAAIIIIAAAggggAACCFgtQMhpdfsoHgEEEEAAAQQQQAABBBBAAAEEEEAAAQQIOTkGEEAAAQQQQAABBBBAAAEEEEAAAQQQQMBqAUJOq9tH8QgggAACCCCAAAIIIIAAAggggAACCCBAyMkxgAACCCCAAAIIIIAAAggggAACCCCAAAJWCxByWt0+ikcAAQQQQAABBBBAAAEEEEAAAQQQQAABQk6OAQQQQAABBBBAAAEEEEAAAQQQQAABBBCwWoCQ0+r2UTwCCCCAAAIIIIAAAggggAACCCCAAAIIEHJyDCCAAAIIIIAAAggggAACCCCAAAIIIICA1QKEnFa3j+IRQAABBBBAAAEEEEAAAQQQQAABBBBAgJCTYwABBBBAAAEEEEAAAQQQQAABBBBAAAEErBYg5LS6fRSPAAIIIIAAAggggAACCCCAAAIIIIAAAoScHAMIIIAAAggggAACCCCAAAIIIIAAAgggYLUAIafV7aN4BBBAAAEEEEAAAQQQQAABBBBAAAEEECDk5BhAAAEEEEAAAQQQQAABBBBAAAEEEEAAAasFCDmtbh/FI4AAAggggAACCCCAAAIIIIAAAggggAAhJ8cAAggggAACCCCAAAIIIIAAAggggAACCFgtQMhpdfsoHgEEEEAAAQQQQAABBBBAAAEEEEAAAQQIOTkGEEAAAQQQQAABBBBAAAEEEEAAAQQQQMBqAUJOq9tH8QgggAACCCCAAAIIIIAAAggggAACCCBAyMkxgAACCCCAAAIIIIAAAggggAACCCCAAAJWCxByWt0+ikcAAQQQQAABBBBAAAEEEEAAAQQQQAABQk6OAQQQQAABBBBAAAEEEEAAAQQQQAABBBCwWoCQ0+r2UTwCCCCAAAIIIIAAAggggAACCCCAAAIIEHJyDCCAAAIIIIAAAggggAACCCCAAAIIIICA1QKEnFa3j+IRQAABBBBAAAEEEEAAAQQQQAABBBBAgJCTYwABBBBAAAEEEEAAAQQQQAABBBBAAAEErBYg5LS6fRSPAAIIIIAAAggggAACCCCAAAIIIIAAAoScHAMIIIAAAggggAACCCCAAAIIIIAAAgggYLUAIafV7aN4BBBAAAEEEEAAAQQQQAABBBBAAAEEECDk5BhAAAEEEEAAAQQQQAABBBBAAAEEEEAAAasFCDmtbh/FI4AAAggggAACCCCAAAIIIIAAAggggAAhJ8cAAggggAACCCCAAAIIIIAAAggggAACCFgtQMhpdfsoHgEEEEAAAQQQQAABBBBAAAEEEEAAAQQIOTkGEEAAAQQQQAABBBBAAAEEEEAAAQQQQMBqAULOCtt3z+Zt8tTuvbJp3a3SNq61wtF4OQIIRClw/ORpWbX2Xrli4Vz56sqlUZYSydzbd+yUO9ZvlS99dkms1x/Wv9v+8bT0uqvl+iVXeT19evdeuemWdXLtNe1y15rl0tRYH0mvmRQBBBBAAAEEEEAAAQQQiJsAIWeejvrfnO55ef+IZwS/+Q/rm2UbDjw/IHnwvrWyeOFcG0o2rkZleP+3/0U2r18ts2dO017faOOrY3nHT56s2tzaF1PkgLmOy3zvWx0h574DB2XlbRvl5hs/nQm2/FJNfo+MVneR1EY9rdCxHsYPp3KFnArJt17yyStjHSYbdUBQDAIIIIAAAggggAACCMRegJAzR4vz7Wbyv2FVL/F3bhJyvgdocoBjyzuZkFN/pwg5C5v29PbLnRu2ek+My+5Ck0NO5Vzt93rhrvMMBBBAAAEEEEAAAQQQQCBeAoScgX76pxJ+/bbl5+zC8p/60PYfyZJr2r3T0wk54/WGiHo11Q4+qj1+1H7Fzh/VTs5i6wv7ef6/e3HahW16yKlj13DYxwnzIYAAAggggAACCCCAAAImCxByZnXH3830+sEjRV9j0w9L7u5YIbd3bhH/9Pbg9db8sR97fNeI4yHXddmKHdMfyN+p5v99wdyL5d3jp845XTZXDbnCXD/wyC60mPAjV6hQ6lryvVly1RSsXc31wPd2ZIYIXlogO1T4xJWXedfF8x/5Qu2grXp+0CL4nGBPs08DVq9X1zz0H9lj5Zorez7f8t6//orc+81HRB1LUydP9E4vf+7F340YV73O/5x/2nux4wevL6vbtZzjK99Ow1zHXPC06+BzguvJtpowvjVzTc5ij5HsY9af+9DbR0ccyv6xGKwl+3Tm4LHhH5PBenMdq8W+t/O9vwr9sCZYw4J5s0b8G6k+f/jIMfnC0j+SW/7q/xJ//f574YW9+0e83/Jdj7LQsebXH3Qu91gv9O+2P1/wmA2uv9R/i7Ofv+3RJ4r+emPyfyaoDQEEEEAAAQQQQAABBBCIWoCQM6sD5Vwnzf+mPPub3lw7dFQIsWHTw/K5638/c63FfMFNsWOq0nNdRzHXtfVyrS1Xnbl2dBW7+y9fyKmCx0I+o70Rcq1R1f6th78vq278jPfS4Km2/tpmTJuUOf02+zqr2QForjXnC7yV4+M/fUb+/PPXSb7nBAOj7EAmO6DK5VVo95myDAY6av3qderh39zE/1jw+p6Fxs++TmGu47NS10qPr+xAKDvYy3ZVc3R0bslcW3S08D0Y6JZyjOQ7ZgtdkzO7J9nzBY8NPwwvdMwU+97OV2+hHYXB4znXceH/m5UdXmYf97k+nn3N0mKPNbWGfMeQ8gr+0CDf9W0r/Tc23zEVvKZtoeucxnEHbdT/qWF+BBBAAAEEEEAAAQQQcFeAkDNHyJnrhiH5DpF8O6DUN8HF7NAJBjJ+cJnrphjBMf1voDs7Voy40U+um134O62C19vLtcst19wv/OZVaWpsGPVmOKWESaX4qB2Xo+0kzRfcBQOGfGGOH7BMmTQhcxOQYoLdXL1T/QvOmy/oyNWnUkLIQv9sVTp+NVzzvV+KOb6Cx7v6+9c6t8j8918sXV3dmTA7eKyXclyWcozoDDmz776txs13w5pK3tvl1DtaAPrDJ56Sq9oXencHz/fvS7Efr/RYy/XvZjnvpVL/jb1i4Vzv34xS/i3O7kO+Phd6b/N5BBBAAAEEEEAAAQQQQACBcwUIOasYcubaRZTrlOHgzrzRgtPsMfN9Ex/8xnm0b6SD35z79eU7nXS0N1EpYVIxIaKaq5gwNF+QEgxoCgVYaj4/BM43Zvb68z0nGJoWCjn9oMRf72i7z0a7I3S+U6UL7QT01xQ87qrhWsnxFTyO1Vi7nnnJO0X6v937HflGxwqZNvk8b1dvocA633uslGOknNBwtNPVs3fhFhtylvLeLlRv8Icl6vn5dssGxyo2zMw+1tTp7YXeb8F+5AsTc713yg05g//G5vthVfaav/+TXZLrfVsoxCz0ef7TggACCCCAAAIIIIAAAgggULwAIWeOkHPJJ6/M7OgrRFlsIOl/M/vusZOZ02jV2KXu5Mz+RrpQUOPvDssXfmWvrdC1IUe7EZM/TjVCzkLXChztrtDBz/X29WWut6h2X/mP4PPUxwvdaTrfNVazTf1T4sMIOf3Tb7P7VMlOzsaGhrwG5bqqHX9+GJV9XVL1sWKOL/U8P1jq+IvPSed/f0jaF82XT32y3atV/fn3PvZBr8erVy7N7G4u5bi0LeQs9b2d69+z0YLDYNDpvz7XD2ayQ8t8YWauj4/2fgsea/61PXPt7A6e+q0j5AxeIzTo5/8waNO3/1ly/QCiUIhZ6POFvv7weQQQQAABBBBAAAEEEEAAgfcECDmzjoZKbjwUvLZf8BvsfN9wVxJy5tvlWMpOzkJvhuwwr9DNh0oJk+K8kzNoWu2Qs9D42adClxL8VGMnZ9CmlOPLD0jVzjoVYm78h3/ydm+qGyv5uzr/5A8+Iv/nvd8Z8YOEUo5L20JOHSFZoetG5juez5swLnPDHBd3cma7FPtvcdBSR/8K/RvO5xFAAAEEEEAAAQQQQAABVwQIOQOd9ncDjbazbOeu57wb6bSNa/V2luW7fmYxuy4rCTnz3bQi+I3zaLsdgwd69nX2/M8VG4KUEiYVG3KOdmMOv9Z8p4qWek1OtV7/9Nl89SnLnbt2yx9dfYUXrOU7tTzbtVAImX26er7rfKrx8h1rxR4HaoxSxq/0Oom5jrtKji9Vv2+pgs1xrc2ZfuW7PqcfjOa6AVOu05B1hJyjBVe6T1cv5b2d74vaaNfdVJ/b8/J+uar98hEvD66j0pCz0mMt1/ujlGPdX1ywjmJvDFTKezAbstjxXfkPCetEAAEEEEAAAQQQQAABBCoRIOTMoedfNzD7DtzqaX4YoP7s79wsNuTM9c2s/7Fyr8mZK+DIdw09Pxz60KXvywRDah3Bu77nWk+xYV41Qk4/vAjetTg7VFHPUacrv37wSKYvfq9y3V09O1T0DYKnp+fa1Ru81qb/92dfeGXEzkE/WFO/q+sslhJyjhYoF7o8Qfa6/PWrgCo7sC9l/FwGlbpWcnxl9+qxx3eNWNdo147MdVwWCqWKOUby/cM7WvCoO+TMDn4LvbdH+0JRyrGV7+7q5Z6uri5jUOyxptaQ6wdR/r/Z2bvNSznW84Wc/r8/D3xvxzk3P1N1/PTJ571Lm5Tyb3F2H4q55nAlX+B5LQIIIIAAAggggAACCCDgkgAhZ55uZ4dE2U8JBp/FhpzZ35z746mxPnHlZdLRuWVESFbKmMFrQ6rAVN1AZOPmbRK8Y3O+60hmrynXc9Su1eDp+LnYqhVy+qFh9nUcc9UUvH5evpC6lAArOGauGzL5AUu2SXZ9pYScuY4TP7gZ7fqkweNVzX93xwq5vXPLOceBHxL59RYaX6drJceXX6+qJxh6Zx8jwcsq5Avps/vm/6BhwvjWoq7b6l9fNN8Xi2A//GOxGiGnqqGY9/ZoX9hG21GY67qfuf4drCTkzO6tChSz/43Mvn6u//FgTcEfFPnPK/VYH+2yIjfdsm4EYXDOUv4tVgONtoPWpf+EsFYEEEAAAQQQQAABBBBAQJcAIacuSYPGKXQjEYNKpRQEEDBAQMdp7wYsw6oSit0hb9WiKBYBBBBAAAEEEEAAAQQQiFCAkDNCfB1Tq11tajfo4oVzM8Pluz6ejvkYAwEE4ilQ7LV347n6cFeFdbjezIYAAggggAACCCCAAAJuCBByWt7n4KnEajnBU0ktXyLlI4BASAL5rkcc0vROTOOfQp/r0hdOALBIBBBAAAEEEEAAAQQQQKBKAoScVYJlWAQQQAABBBBAAAEEEEAAAQQQQAABBBAIR4CQMxxnZkEAAQQQQAABBBBAAAEEEEAAAQQQQACBKgkQclYJlmERQAABBBBAAAEEEEAAAQQQQAABBBBAIBwBQs5wnJkFAQQQQAABBBBAAAEEEEAAAQQQQAABBKokQMhZJViGRQABBBBAAAEEEEAAAQQQQAABBBBAAIFwBAg5w3FmFgQQQAABBBBAAAEEEEAAAQQQQAABBBCokgAhZ5VgGRYBBBBAAAEEEEAAAQQQQAABBBBAAAEEwhEg5AzHmVkQQAABBBBAAAEEEEAAAQQQQAABBBBAoEoChJxVgmVYBBBAAAEEEEAAAQQQQAABBBBAAAEEEAhHgJAzHGdmQQABBBBAAAEEEEAAAQQQQAABBBBAAIEqCRByVgmWYRFAAAEEEEAAAQQQQAABBBBAAAEEEEAgHAFCznCcmQUBBBBAAAEEEEAAAQQQQAABBBBAAAEEqiRAyFklWIZFAAEEEEAAAQQQQAABBBBAAAEEEEAAgXAECDnDcWYWBBBAAAEEEEAAAQQQQAABBBBAAAEEEKiSACFnlWAZFgEEEEAAAQQQQAABBBBAAAEEEEAAAQTCESDkDMeZWRBAAAEEEEAAAQQQQAABBBBAAAEEEECgSgKEnFWCZVgEEEAAAQQQQAABBBBAAAEEEEAAAQQQCEeAkDMcZ2ZBAAEEEEAAAQQQQAABBBBAAAEEEEAAgSoJEHJWCZZhEUAAAQQQQAABBBBAAAEEEEAAAQQQQCAcAULOcJyZBQEEEEAAAQQQQAABBBBAAAEEEEAAAQSqJEDIWSVYhkUAAQQQQAABBBBAAAEEEEAAAQQQQACBcAQIOcNxZhYEEEAAAQQQQAABBBBAAAEEEEAAAQQQqJIAIWeVYBkWAQQQQAABBBBAAAEEEEAAAQQQQAABBMIRIOQMx5lZEEAAAQQQQAABBBBAAAEEEEAAAQQQQKBKAoScVYJlWAQQQAABBBBAAAEEEEAAAQQQQAABBBAIR4CQMxxnZkEAAQQQQAABBBBAAAEEEEAAAQQQQACBKgkQclYJlmERQAABBBBAAAEEEEAAAQQQQAABBBBAIBwBQs5wnJkFAQQQQAABBBBAAAEEEEAAAQQQQAABBKokQMhZJViGRQABBBBAAAEEEEAAAQQQQAABBBBAAIFwBAg5w3FmFgQQQAABBBBAAAEEEEAAAQQQQAABBBCokgAhZ5VgGRYBBBBAAAEEEEAAAQQQQAABBBBAAAEEwhEg5AzHmVkQQAABBBBAAAEEEEAAAQQQQAABBBBAoEoChJxVgmVYBBBAAAEEEEAAAQQQQAABBBBAAAEEEAhHgJAzHGdmQQABBBBAAAEEEEAAAQQQQAABBBBAAIEqCRByVgmWYRFAAAEEEEAAAQQQQAABBBBAAAEEEEAgHAFCznCcmQUBBBBAAAEEEEAAAQQQQAABBBBAAAEEqiRAyFklWIZFAAEEEEAAAQQQQAABBBBAAAEEEEAAgXAECDnDcWYWBBBAAAEEEEAAAQQQQAABBBBAAAEEEKiSACFnlWAZFgEEEEAAAQQQQAABBBBAAAEEEEAAAQTCESDkDMeZWRBAAAEEEEAAAQQQQAABBBBAAAEEEECgSgKEnFWCZVgEEEAAAQQQQAABBBBAAAEEEEAAAQQQCEeAkDMcZ2ZBAAEEEEAAAQQQQAABBBBAAAEEEEAAgSoJEHJWCZZhEUAAAQQQQAABBBBAAAEEEEAAAQQQQCAcAULOcJyZBQEEEEAAAQQQQAABBBBAAAEEEEAAAQSqJEDIWSVYhkUAAQQQQAABBBBAAAEEEEAAAQQQQACBcAQIOcNxZhYEEEAAAQQQQAABBBBAAAEEEEAAAQQQqJIAIWeVYBkWAQQQQAABBBBAAAEEEEAAAQQQQAABBMIRIOTU4HzwaI+GURgCATcFJrTWS3dfSnr7U24CsGoENAhMbmuUd0/2SWoorWE0hkDAPYFEQmRKW5McOsb/6dzrPivWJVBfWyNjm+u8r0c8EECgPIHmxlqpTSbk5JmB8gbgVaEITJvYFMo8TFK6ACFn6WbnvIKQUwMiQzgrQMjpbOtZuEYBQk6NmAzlpAAhp5NtZ9GaBQg5NYMynJMChJx2tJ2Q09w+EXJq6A0hpwZEhnBWgJDT2dazcI0ChJwaMRnKSQFCTifbzqI1CxByagZlOCcFCDntaDshp7l9IuTU0BtCTg2IDOGsACGns61n4RoFCDk1YjKUkwKEnE62nUVrFiDk1AzKcE4KEHLa0XZCTnP7RMgpIvds3iYPfG/HiC59/bblcv2Sq7yPbd+xU+5Yv9X787XXtMtda5ZLU2N95vmEnOYe4FRmvgAhp/k9okLzBQg5ze8RFZotQMhpdn+ozg4BQk47+kSVZgsQcprdH786Qk5z+0TIeTbkVC366sql53Tq6d17ZePmbbJp3a3SNq7VC0SDzyXkNPcApzLzBQg5ze8RFZovQMhpfo+o0GwBQk6z+0N1dggQctrRJ6o0W4CQ0+z+EHKa3x9CzgIhpwo1L5o+JbOrMxh6qhYTcpp/oFOhuQKEnOb2hsrsESDktKdXVGqmACGnmX2hKrsECDnt6hfVmilAyGlmX4JVsZPT3D4RcuY4Xd0/Vb2nt1/u3LBV2hfNz4Sc+w4clK91bpFvdKyQ2TOneZ0l5DT3AKcy8wUIOc3vERWaL0DIaX6PqNBsAUJOs/tDdXYIEHLa0SeqNFuAkNPs/vjVmRRyHj95WlatvVdWr1wqixfOtQMwT5XqUpG7nnnpnEtElrIoQs6AlgoxV962UTo7Vsilc2d5IecN112dOVjOCTmfWidnLrpBhlpnluLOcxFA4KxAU31SBlJDMphKY1KiQDqdloT6zpyH8wItjbXS3TcoQ7yNnD8WAChPQP1L2tJUK6d7BssbgFchgIAkaxLSUFcj3X2pqmvwf6CqEzNBRALqhwU1CZHegaGIKmDaYgRam2qLeZq25+S6j8yXPrvEu+RilCGnP/fS667ObAysZNGEnJXojfJa/xT1T32yvfBOzo0JSdeNld5r/lEGL/5MlSpiWATiK9DUkJSBwbQMpvhCXmqX05KQhJBqleoWx+c3N9VKd29K1Dd9PBBAoAyBhEhLY5109QyU8WJeggACSsALOeuT0t1b/R8W8H8gjrm4CtSpkLMmIX391f9hQVwNw1hX65i6MKbJBJgzpk0asbvRP+tYbcibc/EF7OTM6gY7OXMcmtnX4Sx4Tc4fflnkhQe8Uc6874ty8mN/H8rBziQIxEWA09Xj0knWEaUAp6tHqc/ccRDgdPU4dJE1RC3A6epRd4D54yDA6ep2dDGs09VVHnX4yLFRT9/2d1PefOOn5f/70S/lscd3eYj+Ts9s0ewdoVMnT5TN61d7l2H0Q9PL5s+W51/aN2KMLy77lBei7nl5vzeUf3lH9edcu0j9j/nPz54nuCN1wbxZmZt8q/HYyanh+FcN2PH4Lvnc9X/gjRY8Hb2Yu6sff+a7Mv6nKyWR6paB8fPk+DWPyGDrLA3VMQQC8Rcg5Ix/j1lh9QUIOatvzAzxFiDkjHd/WV04AoSc4TgzS7wFCDnt6G8YIWexp4L7z3v32MlMaJl9GUb/Op0qYFQPdYq7eqisq6Nzi/eaaZPP885ifvaFV84ZQz3XD0ODeVkw5MxVs5pHPVQd3/zuo3LNJxZl7m8TDHEJOTUc/35i7afdasgH71s74oKtCvqO9Vu92a69pv2cFF3deKj29H5pe/wGqTvxsqSTzXLyI38r3XP+TEOFDIFAvAUIOePdX1YXjgAhZzjOzBJfAULO+PaWlYUnQMgZnjUzxVeAkNOO3oYRcuYKKnPp5NpNGbyJthprw/0PS+ftK6RtXKs3TPZzcl2qMdeNuIMfC85dakgZrKvU1+fy4HR1De8h/+7qiaE+GftUhzTv/Yfhg2bmf5QTV22WdHKMhlkYAoF4ChByxrOvrCpcAULOcL2ZLX4ChJzx6ykrCl+AkDN8c2aMnwAhpx09tS3kVLspb7plXU5cdfq5rpAzeLnH4IS5Nhlmn85OyGnI8e+HnH45jW/+QMb/+01SM3BKBltmyfHff8Q7jZ0HAgicK0DIyVGBQOUChJyVGzKC2wKEnG73n9XrESDk1OPIKG4LEHLa0f8wQs5ST1dfvXJp5ozk4I7L4GUY84WP7YvmZ+6SXs5OztFCTn9n6pJPXpk5ZT54+jshpyHHfzDkVGUlz7wpE37yp1J39NeSrmmUU1dukDPv/5IhFVMGAuYIEHKa0wsqsVeAkNPe3lG5GQKEnGb0gSrsFiDktLt/VG+GACGnGX0oVEUYIaeqYbQbD/nXusx1d/Vcp6t/rXOLfKNjReZ6mNlrLCbQVM+v5HR1Ve8jjz4x4vKPhJyFjrSIPp8r5PRKGRqUsc/eKS0v3CciaemdcZ0cv+oBSde2RFQp0yJgngAhp3k9oSL7BAg57esZFZslQMhpVj+oxk4BQk47+0bVZgkQcprVj3zVhBVy+rs5Z0ybNCIc9O8bo+4nU0zI6YeTrx88cs7dzKdPmySXzp3l3Xio0p2c/m5Ndaf365dc5fH5Yaz6s3+jo+w7umff7IidnIYc/3lDzrP1Nbz1uLTtvElq+o5KqmWmHLtmmwy0LTCkespAIFoBQs5o/Zk9HgKEnPHoI6uIToCQMzp7Zo6PACFnfHrJSqITIOSMzr6UmcMKOf2a1I7OB763I1Ni9nUsi7nxUL5xFsyb5YWejQ0NWkJONY8fdB56+6g3bfCam/5NvdXH16z6U/nWw9/P7DAl5CzlKKzicwuFnGrqmt4jMuEnfyb1R34h6Zp6Ob34G9I17ytVrIqhEbBDgJDTjj5RpdkChJxm94fqzBcg5DS/R1RovgAhp/k9okLzBQg5ze+RqjDskNMOFTOq5O7qGvpQTMjpTZMektbn/0Zad98tkk5J7wV/LCf+w4MyVD9WQxUMgYCdAoScdvaNqs0SIOQ0qx9UY58AIad9PaNi8wQIOc3rCRXZJ0DIaUfPCDnN7RMhp4beFB1ynp1L7eZUuzrV7s5U84Vy/OrvSv/5V2iohCEQsE+AkNO+nlGxeQKEnOb1hIrsEiDktKtfVGumACGnmX2hKrsECDnt6Bchp7l9IuTU0JtSQ041pbo+Z9u/3yQNBx8XSdTK6YUdcvqy20QSSQ0VMQQC9ggQctrTKyo1V4CQ09zeUJkdAoScdvSJKs0WIOQ0uz9UZ4cAIacdfSLkNLdPhJwaelNOyOlP2/LC30rrs38liaEB6T9vkZy4+jsy2HKxhqoYAgE7BAg57egTVZotQMhpdn+oznwBQk7ze0SF5gsQcprfIyo0X4CQ0/weqQoJOc3tEyGnht5UEnKq6euOPSdt/+s/Se3p30k62Syn2jfImffdpKEyhkDAfAFCTvN7RIXmCxBymt8jKjRbgJDT7P5QnR0ChJx29IkqzRYg5DS7P351hJzm9omQU0NvKg0252b+AAAgAElEQVQ5VQmJVK+Mfeovpfk3W7yKei9cIic+/k0ZapygoUKGQMBcAUJOc3tDZfYIEHLa0ysqNVOAkNPMvlCVXQKEnHb1i2rNFCDkNLMvwaoIOc3tEyGnht7oCDn9Mhrf+pGM3/llqel7R4YazpcTV/2j9F7wBxqqZAgEzBQg5DSzL1RllwAhp139olrzBAg5zesJFdknQMhpX8+o2DwBQk7zepKrIkJOc/tEyKmhNzpDTlVOTe8xGf+zP5fGN3d41Z15/wo5dcXfSDrZqKFahkDALAFCTrP6QTV2ChBy2tk3qjZHgJDTnF5Qib0ChJz29o7KzREg5DSnF6NVQshpbp8IOTX0RnfI6Zc05pVvy7gn10hisEsGW+fI8d/7HzIw4XINFTMEAuYIEHKa0wsqsVeAkNPe3lG5GQKEnGb0gSrsFiDktLt/VG+GACGnGX0oVAUhZyGh6D5PyKnBvlohpyqttutVGf/EF6T+3WdEErVy+oNfk9ML/qtIIqmhcoZAIHoBQs7oe0AF9gsQctrfQ1YQrQAhZ7T+zB4PAULOePSRVUQrQMgZrX+xsxNyivT09sudG7bK6wePyKZ1t0rbuNZi+ar6PEJODbzVDDm98tIpaX1+g7TuvlskPSj95y2S41c/JKmWGRqqZwgEohUg5IzWn9njIUDIGY8+soroBAg5o7Nn5vgIEHLGp5esJDoBQs7o7EuZmZBTZN+Bg7L5O/8qp7q65UufXSKLF84thbBqzyXk1EBb9ZDzbI11x56Ttp/8mbe7M13bIifbN0r3nP+kYQUMgUB0AoSc0dkzc3wECDnj00tWEo0AIWc07swaLwFCznj1k9VEI0DIGY17qbMScops37Ezw/baG4flqyuXlspYlecTcmpgDSvkVKUmUt0y9qm10vybf/Qq771wiZz4+DdlqHGChpUwBALhCxByhm/OjPETIOSMX09ZUbgChJzhejNbPAUIOePZV1YVrgAhZ7je5c4Wdsj5b/+2T/r7U+WWW/br/uRPLsn5WnWq+oZND8vnrv997/Mb7n9YOm9fYcQp64ScZbf7vReGGXL6sza+9SMZv/PLUtP3jgw1nCcnPr5Jeqdfq2E1DIFAuAKEnOF6M1s8BQg549lXVhWeACFneNbMFF8BQs749paVhSdAyBmedSUzhR1yTpq0Qd55p7uSkst67TvvrJHzzhtzzmvVqeoPbf+xrFm1zPucujZn+6L5cv2Sq8qaR+eLCDk1aEYRcqqya/qOy7hf3iJNrz3iraJn5vVy8qP3eaEnDwRsESDktKVT1GmyACGnyd2hNhsECDlt6BI1mi5AyGl6h6jPBgFCThu6JBJ2yPn5z2+XU6f6Qsf57nevl7FjG86Z957N2+Si6VMyoaY6dX3XMy/JXWuWS1Njfeh1Zk9IyKmBP6qQ0y+98Y3HZPzPvyI1vUdkqGGCnLzyHumZZcb1EDTwMkTMBQg5Y95glheKACFnKMxMEmMBQs4YN5elhSZAyBkaNRPFWICQ047mhh1ymqRy/ORpWbX2Xtnz8v4RZU2dPFE2r18ts2dOi7RcQk4N/FGHnGoJNf0nZeyTq2XMvv/praj3gj+Ukx/bJKkxUzWskCEQqJ4AIWf1bBnZHQFCTnd6zUqrI0DIWR1XRnVLgJDTrX6z2uoIEHJWx1X3qC6HnE/v3isbN2+TTetuHXENzuDuTt3mxY5HyFms1CjPMyHk9MtrfOvfZNzPVkmy55AM1Y2VU1eul+45X9CwSoZAoDoChJzVcWVUtwQIOd3qN6vVL0DIqd+UEd0TIOR0r+esWL8AIad+02qM6HLIqcJM9QjeTV2Fn488+kTkp6wTcmo44k0KOdVyEgNdMvapv5TmV77lra5v6u/JiY9vllTzhRpWyxAI6BUg5NTryWhuChByutl3Vq1PgJBTnyUjuStAyOlu71m5PgFCTn2W1RzJ5ZCzmq46xibk1KBoWsjpL6nh8E9l/E+/JMkzb0q6tkVOffi/yZm5K1QMqmHVDIGAHgFCTj2OjOK2ACGn2/1n9ZULEHJWbsgICBBycgwgULkAIWflhmGMQMgZhnJ5cxByluc24lWmhpyqyMTgGRn7q/9Dmvd+U0TS0n/+R+TEVd+UwdbZGlbOEAhULkDIWbkhIyBAyMkxgEBlAoSclfnxagSUACEnxwEClQsQclZuGMYIhJxhKJc3ByFneW7WhJx+ofVHfiltO5dLsuuApJNNcvpDd0rX/P9dJFGjQYAhEChfgJCzfDteiYAvQMjJsYBAZQKEnJX58WoECDk5BhDQI0DIqcex2qMQclZbuPzxCTnLt8u80uSdnNnLS6R6pPWZv5aWl/5eRIZkYPylcuIT35SBiQs1KDAEAuUJEHKW58arEMgWIOTkeECgMgFCzsr8eDUChJwcAwjoESDk1ONY7VEIOastXP74hJzl21kXcvoF17/ztIz/6Zel9tQrIlIjZ96/XE4v+roM1Y/ToMEQCJQmQMhZmhfPRiCXACEnxwUClQkQclbmx6sRIOTkGEBAjwAhpx7Hao9CyFlt4fLHjzTkPH7ytKxae6/seXn/OStYMG+WbFp3q7SNay1/dSG90padnNkciVS/tLx4r7Q8t17UDs+hhvPk1OJO6Z7zuZDUmAaBYQFCTo4EBCoXIOSs3JAR3BYg5HS7/6xejwDX5NTjyChuCxBy2tF/Qk5z+xRpyHnP5m2ezFdXLjVXqIjKbAw5/WUlzxyQ8T//L9Jw8Mfeh/rPb5cTH7tfBsfPLWLlPAWBygUIOSs3ZAQECDk5BhCoTICQszI/Xo2AEiDk5DhAoHIBQs7KDcMYgZAzDOXy5ogs5FS7ODvu3iJrbl4ms2dOK696za/q6e2XOzds9Ua9a81yaWqs9/68fcdOuWP98MevvaZ9xOfUx2wOOX3Cxjcek3FPrpZk1+siiVrpmv8VOf2hOySdHKNZmeEQGClAyMkRgUDlAoSclRsygtsChJxu95/V6xEg5NTjyChuCxBy2tF/Qk5z+0TIebY3fsD52OO7RgSZT+/eKxs3b8ucOp9r92kcQk7FkEj1SsvuTml58W8lMdQvqTEXyKkrNkjPRZ8x9wimMusFCDmtbyELMECAkNOAJlCC1QKEnFa3j+INESDkNKQRlGG1ACGnHe0j5DS3T5GFnIpEBYYXTZ8i1y+5KnIhvxZVyK5nXsrs1gzWGAw91fPjEnL6Tag9vV/G/ew/S8PbP/M+1DflP8iJT/yDpJpnRt4nCoifACFn/HrKisIXIOQM35wZ4yVAyBmvfrKaaAQIOaNxZ9Z4CRBy2tFPl0PO7A2CfremTp4om9evNuIs7UhDzn0HDspD238sa1Yty5waHsUhnb07U52a7oecqhZ1+nr7ovmZIFbV/LXOLfKNjhWZBsYt5PR70PTa/y1jn/xLSfYcknRNg3QtWC1dl6/x/swDAV0ChJy6JBnHZQFCTpe7z9p1CBBy6lBkDNcFCDldPwJYvw4BQk4ditUfg5BzZE6WnaP5l32sfhdyzxBZyDnandVVqWHdXV0147U3DmdufpQr5Lzhuqtl8cLhG/HkCjmjal4o8w50ifz8r0Se/TuR9KDI2ItEPnG3yNxl6gT3UEpgEgQQyC3QPzjkXeSfBwIIIIAAAggg4JIA/wdyqdusFQEETBLwd3JmbwbMdcZzVDVHFnJGteDgvGoX5wPf23FOOeoGQ2v/y+dl3d9919mdnNkotSf2yvif3yz17+zyPjzQdpmc/Mi90j/pI6a0kjosFWAnp6WNo2yjBNjJaVQ7KMZCAXZyWtg0SjZOgJ2cxrWEgiwUYCenHU0LfSfna/8mMtQfPs6sPzlnzlwhZ3DzYPiFvjej8yFnED+4zdbFa3KOdkA2vfb/SOuv/kpqu171ntZ74RI5dUWnDI59X5THMXNbLEDIaXHzKN0YAUJOY1pBIZYKEHJa2jjKNkqAkNOodlCMpQKEnHY0LvSQ8/5JIj3vhI9z8zsiTeeNmDfXNTnVE7702SWZM6TDL9SgkFNta73plnUjDB68b23m9PCwcYIhp0t3Vy/WOjE0IM2/+aZ3J/aavmMiiaSced9NcnrRnTLUMPINUOyYPM9dAUJOd3vPyvUJEHLqs2QkNwUIOd3sO6vWK0DIqdeT0dwUIOS0o++hh5w7Pi/Sdyp8nGu/K1I/NmfImX26eq7dneEXOzxjpDs5c523r655ufK2jXLzjZ+O5K7ruS6Yqj52x/qtHpg6jf2uNctH3CgprjceKnRQJgZOS+tz66T5pU2SGOqVdG2LdC24VbouvUXSyaZCL+fzCHgChJwcCAhULkDIWbkhI7gtQMjpdv9ZvR4BQk49jozitgAhpx39Dz3kNIglX6BpyinrkYWcPkz2TX38vqnw85FHnzgnTDSoryNKcTXk9BGSZ96U1mfvlDH7HhaRtKSapkjXh/5Kzsz5gkiCm6KYetyaUhchpymdoA6bBQg5be4etZsgQMhpQheowXYBQk7bO0j9JggQcprQhcI1EHKOvLs6OzlFRN1dvePuLbLm5mUye+a0EUeR2s254f6HpfP2FdI2rrXwERbxM1wPOX3+uuMvydhdt0rD2z/1PjQwfp6cXny39F7wRxF3iOlNFiDkNLk71GaLACGnLZ2iTlMFCDlN7Qx12SRAyGlTt6jVVAFCTlM7M7IuQs6t8tjjwzel9h9fv215JGdjB48YdnJqeA8Rco5EbDj4Yxn7VIfUnXjR+0T/+e3StXCt9F7whxq0GSJuAoScceso64lCgJAzCnXmjJMAIWecuslaohIg5IxKnnnjJEDIaUc3XQ45Te9QZCGnglHn7G979AnZtO7WzI7NqK/JWU7DCDlzqKWHZMy+h6T12bsk2X3Qe8LAhIXSdfla6Zl5nbocbDnUvCaGAoScMWwqSwpdgJAzdHImjJkAIWfMGspyIhEg5IyEnUljJkDIaUdDCTnN7VOkIadiMe3u6uW0ipAzv1oi1StjXvmOtLxwjyS7Xh8OO8fPk67L/lJ6Lv7fuGZnOQdczF5DyBmzhrKcSAQIOSNhZ9IYCRByxqiZLCUyAULOyOiZOEYChJx2NJOQ09w+RR5ymktTfGWEnEVYDQ3KmP3/JC3Pr5faU694LxhsnS1dl/1X6Z79OZGa2iIG4SlxFCDkjGNXWVPYAoScYYszX9wECDnj1lHWE4UAIWcU6swZNwFCTjs6Sshpbp8IOTX0hpCzBMT0kDS9/s/Ssnu91B1/3nthqvlC6VqwWrrfd6Okk40lDMZT4yBAyBmHLrKGqAUIOaPuAPPbLkDIaXsHqd8EAUJOE7pADbYLEHLa0UFCTnP7RMipoTeEnOUhNr71Q2nZ/TdS/87wXblSjZPlzKW3yJm5X5Z0bXN5g/Iq6wQIOa1rGQUbKEDIaWBTKMkqAUJOq9pFsYYKEHIa2hjKskqAkNOOdhFymtun0EPO4ydPy6q198oX//SP5Vv/9APZ8/L+nDoL5s0acUMicwlFCDkr607D4Z3S8vwGaTj4uDdQurZFuuf8mXRfslwGJlxW2eC82liBZM/bkug7LuNru6Svr0/6B4ZG1poeEhkalISkRNKp4T8PpdSFDkSGhiShPpYeHP59aNB7znt/fu/z3mu9R83Za8AO/55OqN8Twx9Tf5dk5s+Zj539nKRzv0YkeXacs2P74/pzJZIl+6drkiI19ZKuqRNJ1Ek6WZf15/r3PsYPAkq2jfMLCDnj3F3WFoYAIWcYyswRdwFCzrh3mPWFIUDIGYZy5XMQclZuWK0RQg85/YWosLPj7i2y5uZlMnvmtBHrUzcjeuTRJ+SuNculqbG+WmvXNi4hpx7KuneflZY9G6TpwL9kBlR3ZO++5Ebvup3puhY9EzGKdoGavqNS0/OOJPvekZred6Sm510Z/tjb3se9P/celZr+E1LTd1wSqW7tNbg+YLqm0bu2rRfeete4TcpwYForkkgOh7je3/0/n/145rlcF9fmY0h9YzmQGpJ02uZVUDsC0QkkRKS+rkb6gj9wi64kZkbAOoGahEgymZCBQb4YWdc8CjZGIFmT8PZhDKZ4HxnTlByFNHx+p8nlOV2bkSHnvgMHZcP9D0vn7SukbVyr8Q0i5NTbotquV6Vp7wMyZt//lGTPYW9wFeD0XHyDdL//Rumf9FG9EzLaqAKJVK8kTx+QZPdbUtv1mtR0vSnJrtckeeZNqT3zhiS7DpQlONRwngw1tEnNmPMklahVmzNHPjK7LFVwNxzWeb97v2olXZMV5p39vBfmJWqHwzxvV6X6c52I2hWqfsnw7wnv93Tm4wnvc2pn6NnnqedkPd97jfe5s6/xXq92jr73fDV28DXlJE6JoQFJDPWLZH5Xu1UHRFL9oj43/PEBSaR6ynLnRQgggAACCCCAAAIIIIAAAhUIrCaErkCvqi81MuTcvmOn7HrmJXZyVrX1dgze+PqjMua3D0rjm9/PFDw47hLpvuSL0j3nC15IxqNygWT3Iak9/TtJntovyVP7pPbUK8Mh5un93s7LQo90XaukGifJUON5MtR4vqQaz5ehpsmSbpo4/OfGyTLUMF6G6tu8nqnn+w+uyVlIt7jPJwbPDAe4Z0/p9373T/NXAa76+9DZU/rTWaf3qz97Hw+mzMXNy7PMEGhrrZeTZ/rP/WGBGeVRBQLGC6ifq6mvR0dP9RtfKwUiYKpAbTIh6lTbk2cGTC2RuhAwXqCxPinJGpEzvf4lt4wv2ckCJ37gD51ctw2LDj3kVLs0V962UQ69fTSvz9TJE2Xz+tXnnMZuKig7OavfGXX9xjG/+440/eZBUTs9vUeiVvomfUT6pn9K+i78IxkYP6/6hVg8gzJMnvqdF1wmT74itV6Yuc/7e2KwK+/K1LUhU2MulFTLdEk1Xyip5hky1DpDBptnSKp5uqRaZ0i6pqFsGULOsul4IQIZAa7JycGAQGUCXJOzMj9ejYAS4JqcHAcIVC7ANTkrNwxjBK7JGYZyeXOEHnL6ZY52Tc7ylhLdqwg5w7VvOPTv0vTKgzJm/z+NmFgFbr0X/rH0XfiH0jft9ySdHBNuYQbMVtN7RGpVkHlqnyRPqhDzd14oXHviN6NeB1PtrBwcO0cGW+fI4Lg5klJ/HnuxF2imxoy8Zq7uZRJy6hZlPBcFCDld7Dpr1ilAyKlTk7FcFSDkdLXzrFunACGnTs3qjUXIWT3bSkeOLOSstHCTXk/IGU031O7Dxrd+JA1vfF8a39ghNX3HRhTSN+0ab4dn3+SPy8DEhdEUWYVZ1anlSRVcqiDz1H6pPal2Zw4Hm6PuyEyOkcHx75fBllnDQea42TLYOtsLN9Up5lE9CDmjkmfeOAkQcsapm6wlCgFCzijUmTNuAoScceso64lCgJAzCvXS5yTkLN0srFcQcmqQJuTUgKhhiPojTw4Hnm/9UOqOPXfOiAPjL5XB8XNlcMKlMjDhMhlomyep5pkaZtY7hLoGprrhUo26TqYXZA6fUp48vV9qT7866g1n1O7V4SDzIhkcO1tS49SOzDmSap0tqabJegvVNBohpyZIhnFagJDT6fazeA0ChJwaEBnCeQFCTucPAQA0CBByakAMYQhCzhCQy5wi0pBztOtzLpg3Szatu5W7q5fZWNdfpkLCxgOPSsObP5D6I7+Qmv6TOUnStS0y0PYBGWi7VFJt82Rg/AckXdcs6dpmGaptknTtGO+093RdS0WkNX3vSk3PUanpPyrJ7iNS03d0+Ff3IanpOSzJ7sOS7FE7NF8vOI86tXxg3CWSavaDzFnejkx1irmpQeZoiyLkLNhynoBAQQFCzoJEPAGBUQUIOTlAEKhcgJCzckNGQICQ045jgJDT3D5FFnL29PbLnRu2Svui+XL5B+bIQ9t/LGtWLZOmxnq5Z/M2+cSVl8nihXPNlcuqjJ2c5rep9sReqX/3Kal7e5fUv/O01J14sayih+rHDweftc2STjZ54Wc6kRwxlrpLdaL/mNT0n5KagVOSGDhd0lzpZLOkxkzxroWZap0pqbHqtPKLvd2Ygy0Xy1DjhJLGM/3JhJymd4j6bBAg5LShS9RosgAhp8ndoTZbBAg5bekUdZosQMhpcnfeq42Q09w+RRZyZt94SPFsuP9h6bx9hbdz8+nde+WRR5+Qu9Ys90JP0x+EnKZ3KHd9tSd/K7VdByR5+jVJdr0mydMHRF3vUt28p6b3HS+g1PUYqh8nQw3ne9e+HGpok6HG8yTVNEnS6vcxU2VIhZpN6tfUineN6qo5rHEIOcOSZp44CxByxrm7rC0MAULOMJSZI+4ChJxx7zDrC0OAkDMM5crnIOSs3LBaIxgRck4Y3yqd//0h6fiLz3khpzqNPTv0rNbidY1LyKlL0rxxvGtjeqeWHxNJF1dfuiYp6bqxMlTfKulaFW6OL+6Fjj6LkNPRxrNsrQKEnFo5GcxBAUJOB5vOkrULEHJqJ2VABwUIOe1oOiGnuX2KLOTMPl39+iVXeaeoXzR9iqg/b9+xU3Y98xI7Oc09bqgMAW0ChJzaKBnIYQFCToebz9K1CBByamFkEMcFCDkdPwBYvhYBQk4tjFUfhJCz6sRlTxBZyBmsWJ2+vmrtvbLn5f0ydfJE2bx+tcyeOa3shYX5QnZyhqnNXHETIOSMW0dZTxQChJxRqDNnnAQIOePUTdYSlQAhZ1TyzBsnAUJOO7pJyGlun4wJOc0lKlwZIWdhI56BQD4BQk6ODQQqFyDkrNyQEdwWIOR0u/+sXo8AIaceR0ZxW4CQ047+E3Ka26fIQs7sGw/ZsmMzXxsJOc09wKnMfAFCTvN7RIXmCxBymt8jKjRbgJDT7P5QnR0ChJx29IkqzRYg5DS7P351hJzm9omQU0NvCDk1IDKEswKEnM62noVrFCDk1IjJUE4KEHI62XYWrVmAkFMzKMM5KUDIaUfbCTnN7VNkIaciUTcb+sSVl8nihXPNFSqiMkLOIpB4CgJ5BAg5OTQQqFyAkLNyQ0ZwW4CQ0+3+s3o9AoScehwZxW0BQk47+k/IaW6fIg059x04KA9t/7GsWbVMmhrrzVUqUBkhp7Wto3ADBAg5DWgCJVgvQMhpfQtZQMQChJwRN4DpYyFAyBmLNrKIiAUIOSNuQJHTE3IWCRXB0yILObPvpp5r3QvmzZJN626VtnGtEbCUNiUhZ2lePBuBbAFCTo4HBCoXIOSs3JAR3BYg5HS7/6xejwAhpx5HRnFbgJDTjv4Tcprbp8hCTnNJSq+MkLN0M16BgC9AyMmxgEDlAoSclRsygtsChJxu95/V6xEg5NTjyChuCxBy2tF/Qk5z+xRZyDna3dWf3r1XHnn0CblrzXIrTmMn5DT3AKcy8wUIOc3vERWaL0DIaX6PqNBsAUJOs/tDdXYIEHLa0SeqNFuAkNPs/vjVEXKa2ycjQ051rc4N9z8snbev4HR1c48dKkNAiwAhpxZGBnFcgJDT8QOA5VcsQMhZMSEDICCEnBwECFQuQMhZuWEYIxByhqFc3hxGhpzbd+yUXc+8FNpOTnWX9we+tyMj+OB9a0fc8V3Vc8f6rd7nr72m/Zy62MlZ3sHHqxBQAoScHAcIVC5AyFm5ISO4LUDI6Xb/Wb0eAUJOPY6M4rYAIacd/SfkNLdPoYecapfmyts2yqG3j+ZVmTp5omxev1pmz5xWdTl12vy3Hv6+rLrxM96p8aq+r3VukW90rPDmV6fOb9y8LXMTJBWIqsdXVy7N1EbIWfU2MUGMBQg5Y9xclhaaACFnaNRMFFMBQs6YNpZlhSpAyBkqN5PFVICQ047GEnKa26fQQ06fYrRrckbJ5d/1ffXKpd5uThVqXjR9ily/5CqvrGDoqT5GyBllx5jbdgFCTts7SP0mCBBymtAFarBZgJDT5u5RuykChJymdII6bBYg5LSje4Sc5vYpspDTVBIVYnZ0bvF2kk6bfJ7cuWGrtC+anwk5gzs9CTlN7SR12SJAyGlLp6jTZAFCTpO7Q202CBBy2tAlajRdgJDT9A5Rnw0ChJw2dEmEkNPcPhFynu1N9mn0/jU5e3r7vZDzhuuuzlyjM1fIebpn0NwOUxkChgs01SdlIDUkg6m04ZWaV146nZaE+s6ch/MCLY210t03KEO8jZw/FgAoT0D9S9rSVCv8n648P16FgBJI1iSkoa5GuvtSVQfh/0BVJ2aCiATUDwtqEiK9A0MRVcC0xQi0NtUW8zSeE4FApCGnf2r4npf3n7P0BfNmZa6DGaZL9unql86dVdROztPdA2GWyFwIxEqgqSEpA4NpGUzxhbzUxqYlIQkh1SrVLY7Pb26qle7elKhv+ngggEAZAgmRlsY66erh/3Rl6PESBDwBL+SsT0p3b/U3gPB/IA66uArUqZCzJiF9/dX/YUFcDcNYV+uYujCmYY4yBCINOXPdxKeMNWh/SfZ1OLkmp3ZeBkRghACnq3NAIFC5AKerV27ICG4LcLq62/1n9XoEOF1djyOjuC3A6ep29J/T1c3tU2Qhpyk3HlKnnz/+02fkzz9/ndcl/7T1zo4V3inq3F3d3IOXyuIhQMgZjz6yimgFCDmj9Wd2+wUIOe3vISuIXoCQM/oeUIH9AoScdvSQkNPcPjkfcvrX3Xzs8V2ZLvnX5PQ/sH3HTrlj/Vbvr9de0y53rVkuTY31medzd3VzD3AqM1+AkNP8HlGh+QKEnOb3iArNFiDkNLs/VGeHACGnHX2iSrMFCDnN7o9fHSGnuX2KLORUJMFTwc1lGr0yQk5bO0fdJggQcprQBWqwXYCQ0/YOUn/UAoScUXeA+eMgQMgZhy6yhqgFCDmj7kBx8xNyFucUxbMiDTnVqeEPbf+xrFm1bMTOyCggKpmTkLMSPV7rugAhp+tHAOvXIUDIqUORMVwWIOR0ufusXZcAIacuScZxWYCQ047uE3Ka26fIQs7R7qyuuKK6u3o5rSLkLEeN1yAwLEDIyZGAQOUChJyVGzKC2wKEnG73n9XrESpUVNAAACAASURBVCDk1OPIKG4LEHLa0X9CTnP7FFnIaS5J6ZURcpZuxisQ8AUIOTkWEKhcgJCzckNGcFuAkNPt/rN6PQKEnHocGcVtAUJOO/pPyGlunwg5NfSGkFMDIkM4K0DI6WzrWbhGAUJOjZgM5aQAIaeTbWfRmgUIOTWDMpyTAoScdrSdkNPcPkUacmbf2Xzq5Imyef1qmTb5PLlzw1ZpXzRfrl9ylblyWZURclrRJoo0VICQ09DGUJZVAoScVrWLYg0UIOQ0sCmUZJ0AIad1LaNgAwUIOQ1sSo6SCDnN7VOkIad/d/VPfbJdNmx6WD53/e/L7JnT5Onde+WRR5+Qu9Yst+KGRISc5h7gVGa+ACGn+T2iQvMFCDnN7xEVmi1AyGl2f6jODgFCTjv6RJVmCxBymt0fvzpCTnP7FFnIqW481HH3Fllz8zJv92Z2yKnuur7h/oel8/YV0jau1Vy9s5URchrfIgo0WICQ0+DmUJo1AoSc1rSKQg0VIOQ0tDGUZZUAIadV7aJYQwUIOQ1tTKAsQk5z+2RkyMlOTnMPGCpDQLcAIaduUcZzUYCQ08Wus2adAoScOjUZy1UBQk5XO8+6dQoQcurUrN5YhJzVs6105MhCTlX49h07ZdczL0nHX3xO/m7r/+udrj5hfKusWnuvLL3uaq7JWWl3eT0CFggQclrQJEo0XoCQ0/gWUaDhAoSchjeI8qwQIOS0ok0UabgAIafhDTpbHiGnuX2KNORULGrX5k23rBsh9OB9a2XxwrnmqgUq43R1a1pFoQYKEHIa2BRKsk6AkNO6llGwYQKEnIY1hHKsFCDktLJtFG2YACGnYQ3JUw4hp7l9ijzkNJem+MoIOYu34pkIBAUIOTkmEKhcgJCzckNGcFuAkNPt/rN6PQKEnHocGcVtAUJOO/pPyGlunyINOdXd1Q8fOTbiLuo9vf1y54at0r5oPqerm3vcUBkC2gQIObVRMpDDAoScDjefpWsRIOTUwsggjgsQcjp+ALB8LQKEnFoYqz4IIWfVicueILKQ0w8zb7ju6nNOTefGQ2X3kxciYJ0AIad1LaNgAwUIOQ1sCiVZJUDIaVW7KNZQAUJOQxtDWVYJEHLa0S5CTnP7FFnIefzkaem4e4usuXmZzJ45bYTQvgMHZcP9D0vn7SukbVyruXpnK+N0deNbRIEGCxByGtwcSrNGgJDTmlZRqKEChJyGNoayrBIg5LSqXRRrqAAhp6GNCZRFyGlunyILOdnJae5BQWUIhClAyBmmNnPFVYCQM66dZV1hCRByhiXNPHEWIOSMc3dZW1gChJxhSVc2DyFnZX7VfHVkIadalDotvaNzi2xevzqzm1Pt4lx520a5+cZPc03OanaesREwRICQ05BGUIbVAoScVreP4g0QIOQ0oAmUYL0AIaf1LWQBBggQchrQhCJKIOQsAimip0Qacqo1+6HmobePZggevG/tOdfpjMinqGk5Xb0oJp6EQE4BQk4ODAQqFyDkrNyQEdwWIOR0u/+sXo8AIaceR0ZxW4CQ047+E3Ka26fIQ05zaYqvjJCzeCueiUBQgJCTYwKBygUIOSs3ZAS3BQg53e4/q9cjQMipx5FR3BYg5LSj/4Sc5vaJkFNDbwg5NSAyhLMChJzOtp6FaxQg5NSIyVBOChByOtl2Fq1ZgJBTMyjDOSlAyGlH2wk5ze1TpCGnusP6qrX3yp6X958jtGDeLNm07lburm7usUNlCGgRIOTUwsggjgsQcjp+ALD8igUIOSsmZAAEhJCTgwCBygUIOSs3DGMEQs4wlMubI9KQ857N27yqv7pyaXnVG/IqdnIa0gjKsFKAkNPKtlG0YQKEnIY1hHKsEyDktK5lFGygACGngU2hJOsECDntaBkhp7l9iizkVLs4O+7eImtuXpa5s7q5TKNXRshpa+eo2wQBQk4TukANtgsQctreQeqPWoCQM+oOMH8cBAg549BF1hC1ACFn1B0obn5CzuKcongWIacGdUJODYgM4awAIaezrWfhGgUIOTViMpSTAoScTradRWsWIOTUDMpwTgoQctrRdkJOc/sUWcipSNTp6hdNnyLXL7nKXKEiKiPkLAKJpyCQR4CQk0MDgcoFCDkrN2QEtwUIOd3uP6vXI0DIqceRUdwWIOS0o/+EnOb2KdKQc9+Bg/LQ9h/LmlXLpKmx3lylApURclrbOgo3QICQ04AmUIL1AoSc1reQBUQsQMgZcQOYPhYChJyxaCOLiFiAkDPiBhQ5PSFnkVARPC2ykHO0O6srB+6uHsHRwJQIRCBAyBkBOlPGToCQM3YtZUEhCxByhgzOdLEUIOSMZVtZVMgChJwhg5c5HSFnmXAhvCyykDOEtYU2BTs5Q6NmohgKEHLGsKksKXQBQs7QyZkwZgKEnDFrKMuJRICQMxJ2Jo2ZACGnHQ0l5DS3T4ScGnpDyKkBkSGcFSDkdLb1LFyjACGnRkyGclKAkNPJtrNozQKEnJpBGc5JAUJOO9pOyGlunyIPOZ/evVduumXdCKEH71srixfONVctUBkhpzWtolADBQg5DWwKJVknQMhpXcso2DABQk7DGkI5VgoQclrZNoo2TICQ07CG5CmHkNPcPkUacqqAc+PmbbJp3a3SNq7VU1I3I1p520a5+cZPW3PXdUJOcw9wKjNfgJDT/B5RofkChJzm94gKzRYg5DS7P1RnhwAhpx19okqzBQg5ze6PXx0hp7l9iizk7Ontlzs3bJUbrrv6nF2bKvx85NEn5K41y6246zohp7kHOJWZL0DIaX6PqNB8AUJO83tEhWYLEHKa3R+qs0OAkNOOPlGl2QKEnGb3h5DT/P5EFnKqu6t33L1F1ty8TGbPnDZCSu3m3HD/w9J5+4rMDs9qUfph62OP78pMETxdfvuOnXLH+q3e56+9pv2c8JWQs1rdYVwXBAg5Xegya6y2ACFntYUZP+4ChJxx7zDrC0OAkDMMZeaIuwAhpx0dZienuX2KLOQ0ZSenClu/9fD3ZdWNn/F2japdpB2dW2Tz+tVe+Bo8pf6ezdu8bn515dJMVwk5zT3Aqcx8AUJO83tEheYLEHKa3yMqNFuAkNPs/lCdHQKEnHb0iSrNFiDkNLs/fnWEnOb2KbKQU5GoHZLbHn3CqGtyqtBz1dp7ZfXKpd5p9CrUvGj6lMz1QXNdR5SQ09wDnMrMFyDkNL9HVGi+ACGn+T2iQrMFCDnN7g/V2SFAyGlHn6jSbAFCTrP7Q8hpfn8iDTkVj2l3V1enyn+tc4t8o2OFTJt8nnfd0PZF8zMhZ/bn/dPsCTnNP9Cp0FwBQk5ze0Nl9ggQctrTKyo1U4CQ08y+UJVdAoScdvWLas0UIOQ0sy/BqtjJaW6fIg85TaLxT6H3Q81cp9TnCjlNWgO1IICAOwL9g0OivqHggQACCCCAAAIIuCTA/4Fc6jZrRQABBIoXiDTkVKeCHz5ybMSNfIJBY/FLqeyZ/rxTJk3IXG8zVy3s5KzMmVcjEBRgJyfHBAKVC7CTs3JDRnBbgJ2cbvef1esRYCenHkdGcVuAnZx29J+dnOb2KbKQ05QbD6nW5Ao4/ZZxTU5zD14qi4cAIWc8+sgqohUg5IzWn9ntFyDktL+HrCB6AULO6HtABfYLEHLa0UNCTnP7FFnIqW7w03H3Fllz8zLvLubZD7VbcsP9D0vn7SukbVxrVfUK7Rzl7upV5WdwBISQk4MAgcoFCDkrN2QEtwUIOd3uP6vXI0DIqceRUdwWIOS0o/+EnOb2KbKQ05SdnCpQXXnbRjn09tERXfrSZ5dkTltXd4G/Y/1W7/PXXtM+4vR69TFuPGTuAU5l5gsQcprfIyo0X4CQ0/weUaHZAoScZveH6uwQIOS0o09UabYAIafZ/fGrI+Q0t0+RhZyKRO2S7OjcIpvXr87s5vRDx5tv/HTmjubm8g1XRshpeoeoz2QBQk6Tu0NttggQctrSKeo0VYCQ09TOUJdNAoScNnWLWk0VIOQ0tTMj6yLkNLdPkYaciiXXTsoH71srixfONVctUBkhpzWtolADBQg5DWwKJVknQMhpXcso2DABQk7DGkI5VgoQclrZNoo2TICQ07CG5CmHkNPcPkUecppLU3xlhJzFW/FMBIIChJwcEwhULkDIWbkhI7gtQMjpdv9ZvR4BQk49jozitgAhpx39J+Q0t0+EnBp6Q8ipAZEhnBUg5HS29SxcowAhp0ZMhnJSgJDTybazaM0ChJyaQRnOSQFCTjvaTshpbp8IOTX0hpBTAyJDOCtAyOls61m4RgFCTo2YDOWkACGnk21n0ZoFCDk1gzKckwKEnHa0nZDT3D4RcmroDSGnBkSGcFaAkNPZ1rNwjQKEnBoxGcpJAUJOJ9vOojULEHJqBmU4JwUIOe1oOyGnuX0i5NTQG0JODYgM4awAIaezrWfhGgUIOTViMpSTAoScTradRWsWIOTUDMpwTgoQctrRdkJOc/tEyKmhN4ScGhAZwlkBQk5nW8/CNQoQcmrEZCgnBQg5nWw7i9YsQMipGZThnBQg5LSj7YSc5vaJkFNDbwg5NSAyhLMChJzOtp6FaxQg5NSIyVBOChByOtl2Fq1ZgJBTMyjDOSlAyGlH2wk5ze0TIaeG3hByakBkCGcFCDmdbT0L1yhAyKkRk6GcFCDkdLLtLFqzACGnZlCGc1KAkNOOthNymtsnQk4NvSHk1IDIEM4KEHI623oWrlGAkFMjJkM5KUDI6WTbWbRmAUJOzaAM56QAIacdbSfkNLdPhJwaekPIqQGRIZwVIOR0tvUsXKMAIadGTIZyUoCQ08m2s2jNAoScmkEZzkkBQk472k7IaW6fCDk19IaQUwMiQzgrQMjpbOtZuEYBQk6NmAzlpAAhp5NtZ9GaBQg5NYMynJMChJx2tJ2Q09w+EXJq6A0hpwZEhnBWgJDT2dazcI0ChJwaMRnKSQFCTifbzqI1CxByagZlOCcFCDntaDshp7l9IuTU0BtCTg2IVRji0MHTcvhQl5w40StnzvRLV1e/dHu/D0hXV5/0nBmQ7u7hX+rzPT0D0t83JAODKRkcGJLB1JAMDqRkYCAtg4Pq9yHvd/U59efe3sERVTc01EoymfB+1SQTUpuskZpkzdmPqd/9P7/3efUx9Vz1u3p+bV1CxoyplzFj6qS5uU6axtRLa2u9NI2p8z42pll9rlaam+u9X01NtdLUVCdNzcOfV39XH7fpQchpU7eo1VQBQk5TO0NdtggQctrSKeo0WYCQ0+TuUJstAoScdnSKkNPcPhFyaugNIacGxBKHUAHmgddOyMG3Tov686FDw7+OvH1G3nrzlLx9uKvEEeP1dBW4jmmu9QJTFXx6Aaj/50wgejY49QJUFaSq1wyHqCo4VWGpF7a2qEB1+O9jxzVohyLk1E7KgA4KEHI62HSWrFWAkFMrJ4M5KkDI6WjjWbZWAUJOrZxVG4yQs2q0FQ9MyFkxoQghpwbEHEOowPLV/cdk//4T8tqrx2Tf747L66+dkN+9ckz6+kbuosxVwcSJTXLh9HEybnyDF961tAz/GtNcJy0tDV5wNxzuDYd96nO1dTVSV6t2VCalrq5GamvV74mzv5/9uPf5Gu812Q+1EzSVEhlSO0BTQ5JKpb0/q9/V34f/LJLyfh/+uPp9KPB5tcv0zJkB6Tm7w1TtNPV2nnYPejtOu8/4O08HvZ2pat7u7kHp7u73nhPcYaq7OypAVXa+5dixjVLfUCP19bXS0FgrDfVJqW9ISn19UtRzGxqS0qg+3lAnjU3qz3Xe3xubar3fJ01skoFUWvoHhnSXauV46rhIp9MyNDT8ezot3p/9v3u/e58TGcp63tCQvPdx9XzvtVnPDY5xztjD8+QcO0VvTD+Ymptqpbs35fWcBwIIlCGQEGlprJOunoEyXsxLEEBACSRrEt7/A7sDZzuhgwACxQuo70VrahLS158q/kU8M3SBjev/IPQ5mbA4AULO4pxGfRYhZ/mIKqDb+9K78tqrKrx81wsyX3v1uLzy26OjhnVtExpl9uyJMmVqi0yd1jri1+QpLTJjxrjyi4rBK4dPv1fB54D0+gHo2b+rMFSFoipE9U7XV8Fp98DZsHT4Y+pz6vVnuoZP8x8OV4dP9eeBAAIIIIAAAggggAACCCCAgKsC6fSdri7d+HUTcmpoESFnYUR1avmr+4/L3pffHd6due+YtyNT7dbM91DXopxzyUS56KLxMmvORLn44ja56OLx3sfU53hEIzB8DdMBOXN2x6kKQHt7Ut7u2v7+Qe+6pr3qz97fU9Lfl/J+EtnTM/y83t4BL8D2fvUMSkpd91TtaGWzoNfQRCIhNTXi/QRXnT5Zk0hIwvuz+vjwL/Vx9THvzxL8uP/asx9Xr89+fvY4ar7MmO/NK2pe9XHvtcPXmeVhtgA7Oc3uD9VZIMBOTguaRImmC7CT0/QOUZ8NAuzktKFLIuzkNLdPhJwaekPI+R6i2un3wvOH5aUXj8iLe47Iyy+/I8/9+vCoype8f6LMmjNBLrlkosyaPVFmz2mTi2a1yYQJTRq6wxCmC3BNTtM7RH02CHBNThu6RI0mC3BNTpO7Q222CHBNTls6RZ0mC3BNTpO7815tXJPT3D4RcmrojYshp7rs2+sHTnhB5osvvC0vvfiOF2y++capnKLqFHIVXs6eM1HmvG+C97v6+/QZ471dZjzcFSDkdLf3rFyfACGnPktGclOAkNPNvrNqvQKEnHo9Gc1NAUJOO/pOyGlunwg5NfQm7iGnCjT37zsuzz93WJ7ffdj7/cU9b8vp0/3n6Kkbybx/3nky/wOT5AOXTpL5HzhfPrBgirS0jLxJjwZ2hoiJACFnTBrJMiIVIOSMlJ/JYyBAyBmDJrKEyAUIOSNvAQXEQICQ044mEnKa2ydCTg29iVPIqQLNfb87Jnuee1ue231Inn/ubXlxz+GcN5xRN/xRIeY8FWReOtkLNmfNbvOu5ccDgWIFCDmLleJ5COQXIOTk6ECgMgFCzsr8eDUCSoCQk+MAgcoFCDkrNwxjBELOMJTLm4OQszy3Ea+yNeQsJdC8cPpYWXD5FLl84RTv94ULp8j4tkYNegzhugAhp+tHAOvXIUDIqUORMVwWIOR0ufusXZcAIacuScZxWYCQ047uE3Ka2ydCTg29sSXkPH6sV55+6k159lcH5dfPHpTdzx7KuUMzO9C87GywSaCp4UBhiJwChJwcGAhULkDIWbkhI7gtQMjpdv9ZvR4BQk49jozitgAhpx39J+Q0t0+EnBp6Y2LIOdCf8k41V2HmM0+rUPOQvPH6yXNWO33GOFlw+WS5/PKpw7+zQ1PDEcEQpQgQcpaixXMRyC1AyMmRgUBlAoSclfnxagSUACEnxwEClQsQclZuGMYIhJxhKJc3ByFneW4jXmVCyPnaq8fl2V8dkmefOej9evGFIzI4MDSiztax9bLwg1PlQ4umyYc+PE0WX3GhjBvfoEGAIRAoX4CQs3w7XomAL0DIybGAQGUChJyV+fFqBAg5OQYQ0CNAyKnHsdqjEHJWW7j88Qk5y7fLvDLskLO/PyW7f31IfvXkQe/081899ZYcO9YzYiXJZEIuef958sFFU2XR4gu8YPN9l0wU9Z94HgiYJEDIaVI3qMVWAUJOWztH3aYIEHKa0gnqsFmAnZw2d4/aTREg5DSlE6PXQchpbp8IOTX0ptoh59GjPfL0rjflqSfflF89/ZZ3Gro6HT37cf6kZvmQCjQ/fIF8cNE0+eCiKdLUVKdhdQyBQHUFCDmr68vobggQcrrRZ1ZZPQFCzurZMrI7AoSc7vSalVZPgJCzerY6Rybk1KmpdyxCTg2eOkNOdcfz3+59V55+6i1vl+bTT74lB147MaLKmpqEvH/uebL4ygvkiiuny6IrpsmMGeM0rIQhEAhfgJAzfHNmjJ8AIWf8esqKwhUg5AzXm9niKUDIGc++sqpwBQg5w/UudzZCznLlqv86Qk4NxpWEnL29g961NIcDzTflmV8dlFMn+0ZU1dxc7+3MVNfQXHzlhd7p5y0t7NLU0DqGMECAkNOAJlCC9QKEnNa3kAVELEDIGXEDmD4WAoScsWgji4hYgJAz4gYUOT0hZ5FQETyNkDML/Z7N2+Si6VPk+iVXjWjF9h075Y71W72PXXtNu9y1Zrk0NdZnnlNKyNnVNSBP/vIN+eXPD8iuX74pe557WwYHR94gaNoFrWcDzQu8UHPe/PNFXWOTBwJxFCDkjGNXWVPYAoScYYszX9wECDnj1lHWE4UAIWcU6swZNwFCTjs6Sshpbp8IOUUkO8T8+m3LR4ScT+/eKxs3b5NN626VtnGtooJQ9fjqyqVFhZwnT/TJL352QH75ize8cPOlF9+RoaF05rUqvJz/gUmy+IoLZHH7hXJl+4UyeUqLuUcMlSGgWYCQUzMowzkpQMjpZNtZtEYBQk6NmAzlrAAhp7OtZ+EaBQg5NWJWcShCziriVjg0IWcWYK6dnMGPBUNP9fLsnZzvvtstP//pAdn1ize8X6/89qio62z6j7r6pHzwg1Plyo9cKO0fmyFXtl/ADYIqPIh5ud0ChJx294/qzRAg5DSjD1RhrwAhp729o3JzBAg5zekFldgrQMhpR+8IOc3tEyHnKCFnT2+/3Llhq7Qvmp/Z3bnvwEH5WucW+UbHCpk9c5q8+eYp2f6vv5Vf/vx17/Tz1149PqLb6g7nH/rwVGn/6Az5yEeny0c+Nt3co4HKEIhAgJAzAnSmjJ0AIWfsWsqCQhYg5AwZnOliKUDIGcu2sqiQBQg5QwYvczpCzjLhQngZIWcRIecN110tixfO9Z4ZDDkTibtGtKm1tcE75fwTV82Qj35shiz68DSpra0JoZVMgYCdAk31SRlIDclgKmvLs51LCb3qdDotCfWdOQ/nBVoaa6W7b1CyrobivAkACJQioP4lbWmqldM9g6W8jOcigECWQLImIQ11NdLdl6q6C/8HqjoxE0QkoH5YUJMQ6R0Yed+OiMph2jwCrU212BgqQMhZRMg52k7Oyy//B5k+Y6x87OMz5WMfmy6XL5wiNepfJR4IIFCUQFNDUgYG0zKY4gt5UWBZT0pLQhJCOFyqWxyf39xUK929KVHf9PFAAIEyBBIiLY110tUzUMaLeQkCCCgBL+SsT0p3b/V/WMD/gTjm4ipQp0LOmoT09Vf/hwVxNQxjXa1j6sKYhjnKECDkHCXkVJ8q9ZqcZfSAlyDgtACnqzvdfhavSYDT1TVBMoyzApyu7mzrWbhGAU5X14jJUM4KcLq6Ha3ndHVz+0TIWSDkrPTu6ua2nsoQMEOAkNOMPlCF3QKEnHb3j+qjFyDkjL4HVGC/ACGn/T1kBdELEHJG34NiKiDkLEYpmucQcorI9h075Y71WzMdmDp5omxev9q7sZB6ZH/+2mva5a41y6WpsT7z/Oy7q0fTRmZFwF4BQk57e0fl5ggQcprTCyqxU4CQ086+UbVZAoScZvWDauwUIOS0o2+EnOb2iZBTQ28IOTUgMoSzAoSczraehWsUIOTUiMlQTgoQcjrZdhatWYCQUzMowzkpQMhpR9sJOc3tEyGnht4QcmpAZAhnBQg5nW09C9coQMipEZOhnBQg5HSy7SxaswAhp2ZQhnNSgJDTjrYTcprbJ0JODb0h5NSAyBDOChByOtt6Fq5RgJBTIyZDOSlAyOlk21m0ZgFCTs2gDOekACGnHW0n5DS3T4ScGnpDyKkBkSGcFSDkdLb1LFyjACGnRkyGclKAkNPJtrNozQKEnJpBGc5JAUJOO9pOyGlunwg5NfSGkFMDIkM4K0DI6WzrWbhGAUJOjZgM5aQAIaeTbWfRmgUIOTWDMpyTAoScdrSdkNPcPhFymtsbKkMAAQQQQAABBBBAAAEEEEAAAQQQQACBIgQIOYtA4ikIIIAAAggggAACCCCAAAIIIIAAAgggYK4AIae5vaEyBBBAAAEEEEAAAQQQQAABBBBAAAEEEChCgJCzCCSeggACCCCAAAIIIIAAAggggAACCCCAAALmChByltmb7Tt2yh3rt3qvvvaadrlrzXJpaqwvczRehkB8BXp6++XODVvlscd3eYv8+m3L5folV+Vd8D2bt8kD39uR+Xyh58dXjpUh8J7A8ZOnZdXae2XPy/u9Dz5431pZvHBuQSL//aeeyNepglw8IeYC5byPnt69V266ZZ0ns2DeLNm07lZpG9cacymWh0B+gVLfE9nfM02dPFE2r18ts2dOgxgBBPII7DtwUDbc/7B03r6CrzccJQiUIUDIWQaa+uK+cfO2zH90VSijHl9dubSM0XgJAvEWyH5/+N9grl65NGdAowKZTd/+Z/nisk95X9TVF/mVt22Uzo4VRQU68ZZkda4K+EFl+6L53g8I1Pvia51b5BsdK0b9RjH7Bwz8MM7Vo4d1+wLlvI+C/99DEwHXBYJff1SAueuZl/L+EC34HuI95foRxPpHE8j+QRw/VONYQaB8AULOMuxUaHPR9CmZ3Wh8wS4DkZc4IaC+WHfcvUXW3LwsE8aU8kOB4DelTqCxSAQCAsGf6Bf7vvC/VqnhRvsmFHAEXBAo9X2U6+uXC06sEYHRBFSo+dobhzMbOwr90C0YghZ6PvoIICDeD7PZycmRgED5AoScJdrl+uaSL9glIvJ0ZwRyvTcK/dQ/G6fQzk9nIFmo0wK5fpBW6IcF2Z8v5T3nNDSLj7VAqe8j/0yCQ28fzbh86bNLOGsn1kcJiyskEPzaU+j/af7nZ0yb5O32/P5Pdo0ISQvNx+cRcFGAkNPFrrNmnQKEnCVq+iHnDdddnTl9lpCzRESe7oxAri/SpQQuhYIcZyBZqNMCKpx55NEnRpwOONp7I7jTppT3nNPQLD7WAqW+j4LP98OapdddPep1pWONyOKcFwiezVYo5FRg6jW/2feG/OypPcI1OZ0/hAAoQoCQswgknoLAKAKEnCUeHuzkLBGMX7U/dQAADJNJREFUpzstUMlOTvWf4sNHjnGzFKePIBavBErdgRa8eZevyHU5OZ5cFij1fZQrFOUHBi4fQazdDyzV7/59CAqFnMEfuqn3VUfnFm4+xOGEwCgChJwcHghUJkDIWYYf1+QsA42XOClQ7jU5CTidPFxYdB6BUq8lGByGYIZDC4Fzr3FW6Nq2+c5EyL4eIa4IuCZQ6jU5y9n56Zop60UgKEDIyTGBQGUChJxl+HF39TLQeImzAoXurh4MNDlF3dlDhYXnESh0V+hCp9EScnJoISBS6H3kX4Ozs2OFdzki//lTJk3wdq0V2rGGMQIuCBS6u7r6erPt0Sdk07pbpW1cqwT/zk5OF44S1lipACFnpYK83nUBQs4yjwD1RfuO9Vu9V3MKYJmIvMwJAf8bxcce3+Wt9+u3LR9xPbPskLO3r09Wrb1X9ry8f4QN7zEnDhUWOYqAH7D4740H71ubuS40ISeHDgLFCYz2PgqGnGrE4PODX7+Km5VnIRAvARVU3nTLOm9RC+bNygSa6u/BUFN9LPsSKlyTM17HAqvRKxD8mqNG54Z3eo0ZzQ0BQk43+swqEUAAAQQQQAABBBBAAAEEEEAAAQQQiK0AIWdsW8vCEEAAAQQQQAABBBBAAAEEEEAAAQQQcEOAkNONPrNKBBBAAAEEEEAAAQQQQAABBBBAAAEEYitAyBnb1rIwBBBAAAEEEEAAAQQQQAABBBBAAAEE3BAg5HSjz6wSAQQQQAABBBBAAAEEEEAAAQQQQACB2AoQcsa2tSwMAQQQQAABBBBAAAEEEEAAAQQQQAABNwQIOd3oM6tEAAEEEEAAAQQQQAABBBBAAAEEEEAgtgKEnLFtLQtDAAEEEEAAAQQQQAABBBBAAAEEEEDADQFCTjf6zCoRQAABBBBAAAEEEEAAAQQQQAABBBCIrQAhZ2xby8IQQAABBBBAAAEEEEAAAQQQQAABBBBwQ4CQ040+s0oEEEAAAQQQQAABBBBAAAEEEEAAAQRiK0DIGdvWsjAEEEAAAQQQQAABBBBAAAEEEEAAAQTcECDkdKPPrBIBBBBAAAEEEEAAAQQQQAABBBBAAIHYChByxra1LAwBBBBAAAEEEEAAAQQQQAABBBBAAAE3BAg53egzq0QAAQQQQAABBBBAAAEEEEAAAQQQQCC2AoScsW0tC0MAAQQQQAABBBBAAAEEEEAAAQQQQMANAUJON/rMKhFAAAEEEEAAAQQQQAABBBBAAAEEEIitACFnbFvLwhBAAAEEEEAAAQQQQAABBBBAAAEEEHBDgJDTjT6zSgQQQAABBBAwUOD4ydOyau29sufl/SOq+/pty+VTn2yXOzds9T5+15rl0tRYn3nOvgMHZeVtG+XmGz8t1y+5SkYbR33+ns3b5IHv7cgrsGDeLLnnr78i933zEXns8V3nPO/aa9q9GtRD1aSe8+B9a2XxwrmZ5/b09uf9nP+k7Tt2yh3rh9eU6zF18kRZf8d/lvV//72Miapt07pbpW1ca2YdyketK/vhr9H/XHY9wbn89WSbGnh4UBICCCCAAAIIIIBACQKEnCVg8VQEEEAAAQQQQECXQDCo9MdVH39o+49lzapl0tvX54WgS6+7ekSopwI99fjqyqVSzDjZYZ4fiK5euTRnSDll0gRv3FyP7ODwS59dMuJ5T+/eKzfdss57WTAAHW2s9kXzzwks/XmCtfhBZjCk9A0OvX1UgiHnaOvR1UvGQQABBBBAAAEEEIhegJAz+h5QAQIIIIAAAgg4KKB2NW579InMLsV8BCo87OjcIpvXr5bZM6eJ+vvGzdsyryt2HH98HSHnnIsvkGf3vCJrbl7m1eSHkpfNny0PbvuBdHasGBGg6gw5u7p7paurW2647urMHCr8bGlukp/8/NeZQDhfUOrgocaSEUAAAQQQQAABJwQIOZ1oM4tEAAEEEEAAAdMEguHlaPWpEO/wkWNy65/fILf+9d+P2NlZyjhqDh0hp9p9+dobh72S/d2kG+5/WNTuThXIVjPkVHNeNH2K7HrmJe8UerXbtePuLd7cKvz1d70Scpp2xFMPAggggAACCCBQXQFCzur6MjoCCCCAAAIIIJBTINc1I3Nda1K9OPt07OCp2qWMU0zIWcw1OVXIefkH5sjXOrfINzpWyL/84Gde8Kg+pq4VWu2Q84vLPuWdxq9OuX/j4BEvcPU/Fgw5R1sP1+TkzYkAAggggAACCMRHgJAzPr1kJQgggAACCCBgqUD29SzVEoLXu1QfU6el3//tf8mctp5rqcWMo2snp39Do6d+/bKMH9cqnbevkGMnTocScqrdo95p+v/6vzwGFbROGN864vql7OS09M1A2QgggAACCCCAQJkChJxlwvEyBBBAAAEEEECgGgL5Tj8PXouz0Nz5xtEZcgZveuT/vdo7OVXI6a/jioVzvVPm/b9zunqhI4PPI4AAAggggAAC8RQg5IxnX1kVAggggAACCBgusHPXc7Jg3ixpG9c6olIVFPqngaub+viPfCFnqePoDDlVbQ9t/5EsuabdW0eYIaea+4dPPCVzLr7Qu/kRIafhBzzlIYAAAggggAACVRYg5KwyMMMjgAACCCCAAAK5BNTp1nes3yoP3rc2c5dw/xRr9Xx1U53sa0bmCzlLHUd3yJm9trBDzuy5CTl5nyGAAAIIIIAAAm4LEHK63X9WjwACCCCAAAIRCvgBZXYJua7HqT4/2unqpYxTKOQs9sZD6pqcwYeOkNOvb8/L+73h1W7XTetu9XaKqrvMq4c6PT34yBdycuOhCA9wpkYAAQQQQAABBEIUIOQMEZupEEAAAQQQQAABBBBAAAEEEEAAAQQQQEC/ACGnflNGRAABBBBAAAEEEEAAAQQQQAABBBBAAIEQBQg5Q8RmKgQQQAABBBBAAAEEEEAAAQQQQAABBBDQL0DIqd+UERFAAAEEEEAAAQQQQAABBBBAAAEEEEAgRAFCzhCxmQoBBBBAAAEEEEAAAQQQQAABBBBAAAEE9AsQcuo3ZUQEEEAAAQQQQAABBBBAAAEEEEAAAQQQCFGAkDNEbKZCAAEEEEAAAQQQQAABBBBAAAEEEEAAAf0ChJz6TRkRAQQQQAABBBBAAAEEEEAAAQQQQAABBEIUIOQMEZupEEAAAQQQQAABBBBAAAEEEEAAAQQQQEC/ACGnflNGRAABBBBAAAEEEEAAAQQQQAABBBBAAIEQBQg5Q8RmKgQQQAABBBBAAAEEEEAAAQQQQAABBBDQL0DIqd+UERFAAAEEEEAAAQQQQAABBBBAAAEEEEAgRAFCzhCxmQoBBBBAAAEEEEAAAQQQQAABBBBAAAEE9AsQcuo3ZUQEEEAAAQQQQAABBBBAAAEEEEAAAQQQCFGAkDNEbKZCAAEEEEAAAQQQQAABBBBAAAEEEEAAAf0ChJz6TRkRAQQQQAABBBBAAAEEEEAAAQQQQAABBEIUIOQMEZupEPj/27FjGgAAAIRh/l1jYuGqAUh6jgABAgQIECBAgAABAgQIECBAgEAvIHL2phYJECBAgAABAgQIECBAgAABAgQIEDgKiJxHbFcECBAgQIAAAQIECBAgQIAAAQIECPQCImdvapEAAQIECBAgQIAAAQIECBAgQIAAgaOAyHnEdkWAAAECBAgQIECAAAECBAgQIECAQC8gcvamFgkQIECAAAECBAgQIECAAAECBAgQOAqInEdsVwQIECBAgAABAgQIECBAgAABAgQI9AIiZ29qkQABAgQIECBAgAABAgQIECBAgACBo4DIecR2RYAAAQIECBAgQIAAAQIECBAgQIBALyBy9qYWCRAgQIAAAQIECBAgQIAAAQIECBA4CoicR2xXBAgQIECAAAECBAgQIECAAAECBAj0AiJnb2qRAAECBAgQIECAAAECBAgQIECAAIGjgMh5xHZFgAABAgQIECBAgAABAgQIECBAgEAvIHL2phYJECBAgAABAgQIECBAgAABAgQIEDgKiJxHbFcECBAgQIAAAQIECBAgQIAAAQIECPQCImdvapEAAQIECBAgQIAAAQIECBAgQIAAgaOAyHnEdkWAAAECBAgQIECAAAECBAgQIECAQC8gcvamFgkQIECAAAECBAgQIECAAAECBAgQOAqInEdsVwQIECBAgAABAgQIECBAgAABAgQI9AIiZ29qkQABAgQIECBAgAABAgQIECBAgACBo4DIecR2RYAAAQIECBAgQIAAAQIECBAgQIBALyBy9qYWCRAgQIAAAQIECBAgQIAAAQIECBA4CoicR2xXBAgQIECAAAECBAgQIECAAAECBAj0AgMpMjhNrtTcTQAAAABJRU5ErkJggg==", "text/html": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# Same plot, but with smooth line\n", "fig = px.line(data_frame=bio.get_history(), x=\"SYSTEM TIME\", y=[\"A\", \"B\"], \n", " title=\"Changes in concentrations with time (smoothed)\",\n", " color_discrete_sequence = ['navy', 'darkorange'],\n", " labels={\"value\":\"concentration\", \"variable\":\"Chemical\"},\n", " line_shape=\"spline\")\n", "fig.show()" ] }, { "cell_type": "markdown", "id": "8d1fb634-3575-44fb-aa0c-0a9110fe54bc", "metadata": {}, "source": [ "The early **OVERSHOOTING** of the equilibrium values shows prominently in the last plot!" ] }, { "cell_type": "code", "execution_count": null, "id": "6e80c184-9e85-4f2a-8426-35aaba9f2628", "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 }