{
"cells": [
{
"cell_type": "markdown",
"id": "49bcb5b0-f19d-4b96-a5f1-e0ae30f66d8f",
"metadata": {},
"source": [
"## Comparing the dynamics of the reaction `A <-> B` , with and without an enzyme\n",
"\n",
"LAST REVISED: June 23, 2024 (using v. 1.0 beta36)"
]
},
{
"cell_type": "code",
"execution_count": 1,
"id": "cbb1af2e-3564-460e-a4ae-41e4ec4f719f",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Added 'D:\\Docs\\- MY CODE\\BioSimulations\\life123-Win7' to sys.path\n"
]
}
],
"source": [
"import set_path # Importing this module will add the project's home directory to sys.path"
]
},
{
"cell_type": "code",
"execution_count": 2,
"id": "087c0d08",
"metadata": {
"tags": []
},
"outputs": [],
"source": [
"from life123 import ChemData\n",
"from life123 import UniformCompartment"
]
},
{
"cell_type": "markdown",
"id": "d6d3ca49-589d-49b7-8424-37c7b01bcacf",
"metadata": {},
"source": [
"# 1. WITHOUT ENZYME\n",
"### `A` <-> `B`"
]
},
{
"cell_type": "code",
"execution_count": 3,
"id": "23c15e66-52e4-495b-aa3d-ecddd8d16942",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Number of reactions: 1 (at temp. 25 C)\n",
"0: A <-> B (kF = 1 / kR = 0.2 / delta_G = -3,989.7 / K = 5) | 1st order in all reactants & products\n",
"Set of chemicals involved in the above reactions: {'A', 'B'}\n"
]
}
],
"source": [
"# Initialize the system\n",
"chem_data = ChemData(names=[\"A\", \"B\"])\n",
"\n",
"# Reaction A <-> B , with 1st-order kinetics, favorable thermodynamics in the forward direction, \n",
"# and a forward rate that is much slower than it would be with the enzyme - as seen in part 2, below\n",
"chem_data.add_reaction(reactants=\"A\", products=\"B\",\n",
" forward_rate=1., delta_G=-3989.73)\n",
"\n",
"chem_data.describe_reactions()"
]
},
{
"cell_type": "markdown",
"id": "0e771dda-1c0f-4fc0-ab21-049740643897",
"metadata": {},
"source": [
"### Set the initial concentrations of all the chemicals"
]
},
{
"cell_type": "code",
"execution_count": 4,
"id": "5563e467-a637-44fa-9ba1-d35ddd82c887",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"SYSTEM STATE at Time t = 0:\n",
"2 species:\n",
" Species 0 (A). Conc: 20.0\n",
" Species 1 (B). Conc: 0.0\n",
"Set of chemicals involved in reactions: {'A', 'B'}\n"
]
}
],
"source": [
"dynamics = UniformCompartment(chem_data=chem_data, preset=\"fast\")\n",
"dynamics.set_conc(conc={\"A\": 20.},\n",
" snapshot=True)\n",
"dynamics.describe_state()"
]
},
{
"cell_type": "markdown",
"id": "651941bb-7098-4065-a598-e50c0b641ab3",
"metadata": {
"tags": []
},
"source": [
"### Take the initial system to equilibrium"
]
},
{
"cell_type": "code",
"execution_count": 5,
"id": "76f24d9a-a788-41d8-90a4-db87386f91aa",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"16 total step(s) taken\n",
"Number of step re-do's because of negative concentrations: 0\n",
"Number of step re-do's because of elective soft aborts: 1\n",
"Norm usage: {'norm_A': 10, 'norm_B': 9, 'norm_C': 8, 'norm_D': 8}\n"
]
}
],
"source": [
"dynamics.set_diagnostics() # To save diagnostic information about the call to single_compartment_react()\n",
"\n",
"dynamics.single_compartment_react(duration=3.0,\n",
" initial_step=0.1, variable_steps=True)"
]
},
{
"cell_type": "code",
"execution_count": 6,
"id": "e58db01b-b932-4f60-91c2-a578353a3702",
"metadata": {},
"outputs": [],
"source": [
"#dynamics.explain_time_advance()"
]
},
{
"cell_type": "code",
"execution_count": 7,
"id": "4a19ad2a-fbd2-420a-b958-2daf88bcc841",
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
" \n",
" "
]
},
"metadata": {},
"output_type": "display_data"
},
{
"data": {
"application/vnd.plotly.v1+json": {
"config": {
"plotlyServerURL": "https://plot.ly"
},
"data": [
{
"hovertemplate": "Chemical=A
SYSTEM TIME=%{x}
Concentration=%{y}
SYSTEM TIME=%{x}
Concentration=%{y}
Reaction `A <-> B` . Changes in concentrations with time (time steps shown in dashed lines)"
},
"xaxis": {
"anchor": "y",
"autorange": true,
"domain": [
0,
1
],
"range": [
-0.002338945466491458,
3.5622139454664903
],
"title": {
"text": "SYSTEM TIME"
},
"type": "linear"
},
"yaxis": {
"anchor": "x",
"autorange": true,
"domain": [
0,
1
],
"range": [
-1.1111111111111112,
21.11111111111111
],
"title": {
"text": "Concentration"
},
"type": "linear"
}
}
},
"image/png": "iVBORw0KGgoAAAANSUhEUgAAA6sAAAFoCAYAAACxAW22AAAgAElEQVR4Xuy9CZgdRbn//84+k2QmyWTfF8ISwo5EQFZBRNkCVyBcEaIsgope8aci3D+oP1G5XvFBRBaXC4KXAP4gbIrIvkoUWUIWIIQshGxkmUyS2Wf+9daZOunpdJ9T1V29nfn2PPOcmXNq/bzVffrbVe9bZT3iIBwgAAIgAAIgAAIgAAIgAAIgAAIgkCICZRCrKbIGmgICIAACIAACIAACIAACIAACICAJQKxiIIAACIAACIAACIAACIAACIAACKSOAMRq6kyCBoEACIAACIAACIAACIAACIAACECsYgyAAAiAAAiAAAiAAAiAAAiAAAikjgDEaupMggaBAAiAAAiAAAiAAAiAAAiAAAhArGIMgAAIgAAIgAAIgAAIgAAIgAAIpI4AxGrqTIIGgQAIgAAIgAAIgAAIgAAIgAAIQKxiDIAACIAACIAACIAACIAACIAACKSOAMRq6kyCBoEACIAACIAACIAACIAACIAACECsYgyAAAiAAAiAAAiAAAiAAAiAAAikjgDEaupMggaBAAiAAAiAAAiAAAiAAAiAAAhArGIMgAAIgAAIgAAIgAAIgAAIgAAIpI4AxGrqTIIGgQAIgAAIgAAIgAAIgAAIgAAIQKxiDIAACIAACIAACIAACIAACIAACKSOAMRq6kyCBoEACIAACIAACIAACIAACIAACECsYgyAAAiAAAiAAAiAAAiAAAiAAAikjgDEaupMggaBAAiAAAiAAAiAAAiAAAiAAAhArGIMgAAIgAAIgAAIgAAIgAAIgAAIpI4AxGrqTIIGgQAIgAAIgAAIgAAIgAAIgAAIQKxiDIAACIAACIAACIAACIAACIAACKSOAMRq6kyCBoEACIAACIAACIAACIAACIAACECsYgyAAAiAAAiAAAiAAAiAAAiAAAikjgDEaupMggaBAAiAAAiAAAiAAAiAAAiAAAhArGIMgAAIgAAIgAAIgAAIgAAIgAAIpI4AxGrqTIIGgQAIgAAIgAAIgAAIgAAIgAAIQKxiDIAACIAACIAACIAACIAACIAACKSOAMRq6kyCBoEACIAACIAACIAACIAACIAACECsYgyAAAiAAAiAAAiAAAiAAAiAAAikjgDEaupMggaBAAiAAAiAAAiAAAiAAAiAAAhArGIMgAAIgAAIgAAIgAAIgAAIgAAIpI4AxGrqTIIGgQAIgAAIgAAIgAAIgAAIgAAIQKxiDIAACIAACIAACIAACIAACIAACKSOAMRq6kyCBoEACIAACIAACIAACIAACIAACECsYgyAAAiAAAiAAAiAAAiAAAiAAAikjgDEaupMggaBAAiAAAiAAAiAAAiAAAiAAAhArGIMgAAIgAAIgAAIgAAIgAAIgAAIpI4AxGrqTIIGgQAIgAAIgAAIgAAIgAAIgAAIZFqsLnxnOZ118ffpM5/8OP331Zf2sebPbp5Lt9/zGH33q+fQeWd+us9nn//aj+j1t5bSwmduJ1WGM8HokY30yx99XZatc6g6Zhwzhw7YZxr98Vf/uUs2Vc+cs0+kb186e5fPjzvrclq7flOf973azgm4Hq8+82ePP/sP+uY1N5Gqh9PqHH7t0smLNCAAAiAAAiAAAiAAAiAAAiBgm0CmxSrDYJHHx5P3Xt+HjRKkXqLOS+wpcfuLH3yVTjj6EE/OTpHrlSCIWP3DfX+l6266exeR6xTRLKqdh4lYdbfz//zwZvrLU6/Qvbd9n2bsMdn2eEJ5IAACIAACIAACIAACIAACIGCFQObFqp/4UsKRZyudQlbNPLpnLZMQq0qQFpuNdX8OsWpl7KMQEAABEAABEAABEAABEACBFBPIvFj1Ep/qPZ4l5SWxzllEJUrdM4tJiFWdOr3EeFJiVc0CO8eze9aXZ7p5GfX5Yuk1s1eHW3CrWWqvc4OXJH/2uEN9l3hzHucstnOJ9boNm+XMsTrUQwl3fX4z6O5l035LsVN8TqNpIAACIAACIAACIAACIFASBDIvVpVwcS73ZRH4+sKl0neUxYdTcLBocc+2chk6wtH2MmC/tjhHlpcYT0KsKtHsFHleQlr53jrFaTF/XdVflVc9SPCbNXc/cHAumXbaWqXj8r3edz6wUJyd48jt/1sSZzw6AQIgAAIgAAIgAAIgAAIZIVASYtUt+vj/A2ZMk4GM+G8+VNAjP6FnS6wWs7szkJGfv62XWHXmi1us+i2d5nZyHw7cZ/d8gCs1s+oOMlVMmCth6hSVfiLXXadfOr92e73v1241LtwzyMXsjM9BAARAAARAAARAAARAAATCESgJseqeaWMxp2YAnWKjkOiyJVaL+Z9mUaz6LZ3moed+GBBErKrlxV4Rid0i17nEWwXC8hOruu8XmvlVbUNAqnAXGuQGARAAARAAARAAARAAAVMCJSFWnSKUAXB0XTUTpoQIi9c3Fr0nt7PxmiVLQqwWm23kvgRdBuznaxkkGrDK4ze42EdVBbEyFavFgky5++8Wx9wmXVGq2u9Or+oodPIUihJtetIhPQiAAAiAAAiAAAiAAAiAQHECJSFWuZtqaSj/vWb9xj57nfJnJx47U/qx8uG1D2oSYlWnTj+/UBaIXv0oNHvMfQ8iVk2WwpqKVRXQqNAyWyVQr/z6uTLokluIhxWrKj+CKRW/YCAFCIAACIAACIAACIAACMRFoGTEKouw1956V3JjYcr+qurgz1jAvv7WUvJaasrpdISj7QBLxWYV/T5XgrPQDLGf+AsiVosJYOdgNRGrimexWUu1FJeDH7GN3XvqhhWr6mGH3wOAuE5G1AMCIAACIAACIAACIAACILCTQMmIVee2Km7xU+gzhSIJscp1KyHoXEpb6H3VXufWLeq9Qr6fTuHOW7uY+mB6RQPmMrnONxcvMw6wpHjrzmaqSMFe6W2IVa9owE7WCLCEyyYIgAAIgAAIgAAIgAAIxEugZMSqc/sSt7DQ/UyhdwtH9b7tmVWnqb32HS0m5LzyFJulDDKz6hbD7iHqrFN3ZlWJT6/h7tw+xi2yvUSjDbHqFKbuNnm1J97TFLWBAAiAAAiAAAiAAAiAQP8jUDJitf+Zrv/0WHef1v5DBD0FARAAARAAARAAARAAgdInALFa+jbOfA/DzAZnvvPoAAiAAAiAAAiAAAiAAAj0UwIQq/3U8FnptppVxVLcrFgM7QQBEAABEAABEAABEAABOwQgVu1wRCkgAAIgAAIgAAIgAAIgAAIgAAIWCUCsWoSJokAABEAABEAABEAABEAABEAABOwQgFi1wxGlgAAIgAAIgAAIgAAIgAAIgAAIWCQAsWoRJooCARAAARAAARAAARAAARAAARCwQwBi1Q5HlAICIAACIAACIAACIAACIAACIGCRAMSqRZhZLEpF2+W233vb92nGHpOz2I3E26y210HU4sRNgQaAAAiAAAiAAAiAAAiUCAGI1RIxZNhuzDhmDkFohaP4s5vn0u33PJYZ0a8eVMw5+0T69qWzw3UeuUEABBIn8If7/krX3XQ3/eIHX6UTjj4k8fbE3QD+HuOjv/Y/bt426sO9R3CKx511OY0e2Uh//NV/Bi/EJ2ec9wePP/sP+uY1N9F3v3oOnXfmpwv25fNf+xGtXb+Jnrz3eut9Dlsg24PbptOPsHX1t/yZF6vqy9nLcGmZKfS7GKu2p0EsqIuFzpc8Xyxef2upvEim8YIR9iRWotOrnGK2YjZ8RPHlYdIvZU93HudFNM4vI5O2J5WWz9NSHdNxMY3y5on7kIVrqW3W6gZo4TO3y6ILnbc2xKpaJaLqU/1R9R6wz7TEr29ejNX3kvt7329Mpun712TMlNp1u5TEatxjKsrrbZzjLIxYTdt1SV0/06I/TK4taU5bMmLV/SRDfXGl4QlHKd1gOZcN88DWEbcmJ4C6MTM90f1uVEzqVmmVWHX3Le4voiBt5zx+T/dU+9XNZpxfRkH7Emc+iNXwtKO8eYJYhVj1G6GFbnYhVsOf11GWALEanG6U19s47w9KSayq76m0PtQLPtqSzVmyYjVNA6aULsZOIcfLNmwtHVZiU84c9M4gmJ4azjJMxa6zLj+xqoRgVMtuTPvrlV7nqR5z4pnfOL+MbPQNZaSfQJQ3T4XEavrJ2GthUjOr9npgv6RCSwOjHpP2e1O4xFK7bpfS/VHcYyHKsR3nOAsjVuNmrlOfjRUuOvX0pzT9Uqy6l3l6Lf1Ts1POweC3RNA928h5VFrlQ+NVjt8J6lWee5bPmffNxcvoL0+9kq8iytlk58VRicOg4pIbbEOkuk/YsKK1kFhle6b1iZkaN7oPEJxfRsyQ/W3V4V7u7Les2J3OpExVl/tcY768zNxvzDvt7fVQwl2e7tJe9xe/zXPM63ribHtU57yNep3LQt3XM+e571UX24rTOMfFqOFDpV8lH+pc8rpOus+ztFxLFY9i49D5ud+13utm070M1ynEvMaJk6PzJukO4b/K55E6dFbBePWN86u8bmHh7Ndfhd+ZV33uceF1zuqMf78bs0I31YXGpJdNnA/7zrr4+/kq1TXEfR30+y7wck8y+Z50j3X3ahg3C3c7nN+BnNb9nWB6bXOXp85rP5s43y92PVZjar/pU/PXBa82qzKL9c0rdoTfGHGff6Ztcfe/0Jhi+xe6fhZj6TWmeFw6H577XR+87gm80jrTBfkuL2YbPxtyH3R9Pb0eTLmvS0Gug27beN1Hm1z78RCm2Ig2+7xkxarfjbvXclEv0cUXMR6sziAV/B4fTj9NdVK4Bzan/eWPvi6j6/oNWq8Lm1d5Xu85vzSddUcZ5MfdXr++6wzBKESqu96gdfiJ1UIiVqfPUacxtYfzy8opOr2eCrLtWWA4x74aD868JmUyDy/xr74QnDfXir3XWHemc5dnIuD9xCq3M+g55lc/95FvzjiYRBTnvEm9XrZ39tf5Be0Wp+5VBn5P+p3jwuvGie3mFjFeYyPpa6nfd4X7e8HrXPcTepzW2Xc3Q/fNmc7MKpfpPC/8/FC92lkobaE+eNWn0w4v9wqTa5lK67eaxm9MFhIW7gdc6kbW6333ePa6fumseJEC8J3lxCLZWab7OltsxsvLncF9LhW7f3Da0uu+R/ehrc71WLH16rP7vsrvmuAU6152Vd8fTvt5cTRpi9e5U2hMuQW+yayo372HH1/ndVpnTHHb3OdckO9yr/PD/SDF5D7ai7GJWNW5/iibFTrnuBzTa3+aA0FFfS8aRfklK1YLnRBeT5j5pC8WPMfrS1Hn6YnJDZbfBcx9A1FsVrZYX4IMJq+bGN0vLVWfc1bb5ElzkPaqPKYBoQoFWHJ/4YRpl+28pmK60E2PzrhWF3B+VQGlTMr0uyl2i+VCZToDWhVaSsRsikU8LjSz6oxQWOxm0WlXnS+sKM55nXr9bMx2ee2td/MPJvzs5PVgrJhYNbkuqfKd14mkr6V+9euML7dNOM/rC5fK2UjnTbn7uyiIWC20EqdYtM0gYtUtKvyWwXm9rzsO/a6XXuPEmTaIWHV/N/mdT+73i/nOHrjP7vTfV1/qe+kvtHxQjbFC1x8/27kDKBZrp3oIVaiuYmNe93rsZ3+33fwexHvV4y5TfU/wuaYeariZsFF02+JnwGIzq858uhMLxVYO6Lglua/pfnVzXa++8bZ8iGrju9zN2MSGfoxNxKrOddDv+uC+rphe+4tdl2zf/5V6eSUjVt2GKjSF7yWSeMC6v0j8RIs6AXTX2eveYOk8MS9Wt8mNtOng9uqH7hPjMGLVyw4mS3GDilWvhxpeT81NOUaVPg6xWmx5vMkXnN+XhPuGrdANnPNL2PkkWGfJo9sOUYjVYqI/qnO+WL3FbiR5Sbi6TsYlVv2WoDpnzJK+lqrxbyK61ThzP+zkss4XN4XsxrFm/Ub5wEfZxdnnLIpVv/Hlft9kHPpdN9034sXOa/W5ibDwOwfc7xcSIM4Ha8XEDn/uN1McRLyoPOq+qBB3d5/UbKOpe5Hu9VhXIBZ6AOcuw+t7gXl+/T9/SSceO1M+uGRbPfb0/D6rhXTbUsx+XitT3PedumK10Pef33eo17Jt54N255Jiv0kDG9/l7nHnZ0Pde2n1gNy9dY3bbn7M3PXofP+q89D02g+xavdOt2TEqvvJtNfF3u8EVki9/EKcXxgmN9FOM+neYBU6Yd1LNOKeWfXyl3D20fTmLegSXZPhH7SOYqLPdDbZpM1h0notlS1UnsmXkSrbyzfK+cVhUqbfeeE+z4rNdDu/hJ03SM6+68zi2xarOg+Oojjndeotdj47r59Ri1XnkkSnnbxu5tJwLfV6YKPzcMR548azIRygjvvrFKh/fvLvu9xAl7JYNRmHWRCrfg9cVNt1/Oed54NXvmLX2ELXfPU9bSJWuTwlWJ1l6wQx1Lke6wrEQstm3d/JzgdDfE7xCgZ+GOQUqHxejRk5rM9Mt25b4hSrhUStm4nqt3uceV3Dvc495/d7sXHmXDbrNT687g9NVjn4MbY5s+p1rrnrdS+J5/sd5+F37YdYDXM3u2vekhSrzguk8+an2BNYhcfv4uC+idZ9GqR7g6XzlCepmdVCTzW9llzrDtOggrJQ+c4ydb5Q3WUVE6th+qvLJUg6NX6CBFhyL5F1j9lCT3CjFqvFfNIKsVJ5dWbibYtVdZNXyB5RnfO6M6s6IitqsWpSftqupX7fNX5jUq3gGTViaP4GWo0TfuDKQYrcN9ClLFa9lmKaXvuK3RSa3CDrzqCqNvrNrOo8HNPtp+qfuo4Um1ktttSY6y1071LMhUDl1RHe7j56XY91BaLJzKrzwRCfUwfMmCZnU50Phtg32Msn1ut6retfajJbrzuzaiJW/Rjp+KyrNOqBholY9VqZ6DW+Tc7FOMSqe+ZX95zkdMWu/cWuSyZ1IS1RSYpVNqzXRbHQcgrnYND1pSt0M8p1Hbz/njLAkt+J7HVhi8J/LexALyaCiok7nfptiNawItX9sMLvRj6tM6vcfsWgkAjh8c1+UyZfRn43FCY30l5+RO7lPM5zV/VBZ6ZQ5eMbKGdQNMXEqx73uIxCrPrdPHCflG9QFOe8Tr3FBK3fDbn7PHEvWeXPlQ+zSlvIhn5t9bpJS/Jayn3gWRr3gx3dm05moR6Y8t+8BFj5j/L7vBTY7b/qN37Zdl6rWXSXvxW6Jhe6yXKPGd3lvqo+Hf9Cne8LZ5piD7P8lt+aCAtdEav7ANuvj8z+s8cdKu8bnIf7GuFnfx1RwuUWaqfTxpyOr51uP2edetguOtdjXbFq6u+omDn9VLnvfu8Xup9LUqwWspW7XbrfJeoewD0OnbYwuT/QGQ/q+sc7V7gf5picNzZnVp3jwf2d5WQT5Npf7KGP6XWuv6cvWbHKhnU/kVRf/O6Ll0rLr3wj4l5yy+87l0wUi1Lq/vIsFvTA+XRP5S0WoTPOZcDFxGgxMWtykqkldqYzokqkmebzaluh/gZtXzEGqk7T5dRe5ao2+gU9cS9396rTS1jyeeP8kvHyBTb5glNj2C+acLHzTN108RYd/EXjde6YjM0oxKrqo9eWAMo+UZzzQetlpsyMfbtU5GeTmU+/G8pC48LrfHMu/Xae00leS1UfdKJb+p3vfr5izuVo7muY102P342QDbGqK2QKiR4TEes1/r3GoR/TYg+zTESOrihVbfFKr95zPzDkfrJ/cqEAS16uHF7XyWKzjF6znpyHH5DwwzyvMp33R+o673UdUaKu2GoV3euxrlhV9epEnHUKI3c7FWMvRiZt8RqPJg9ATB5yqe9ar+9fZ//U2PO6ZspzqncPe68x6r7vNfku97ONGlNq3Hl9HztnKHX8om2LVb8xrvrk3HrN5Nqv+zC42L0hPs8RKGmx6rxgOS9Mfj5wXv6paqDwzabaB0xn/0d3Gqefk3vPNvcJ6uXnoRPVTH3B8/IWG6JH9V1n2avXxTSrJ1khH8kgS590ONgUq1yfny+YcxyZfhm5fb65LBaKQZcBO292FSPmy19svE2O+6bdz7/Ea1sZJ3PdcyEKsara4fbp8duf1tluG+d8kHq5Dc72mYhV5w2v6ot7n1WvqMzuc47r54MDPbnHQZLXUq8xWOym3WnTQg9P/K6zfsLEaVvVBhtildvrtofOPqvO2TcTser83nJfK3XP3WKzGO5rl9NXWCcYjqmI1bn++n0vePm9et3Ee9lflenuL7/v/O5yLuV1+uB5fb959UXX1cTru9RtU1OB6O6bX1u8Jh2cY81vC620LQP2symfk/w96Y4G7B4/PHb4IYl7RtNrjIS5P/C69rvHnde5zu3nNrIPfxJi1e87mt/32s7Ged4W22dZx83G7zqA9/sSyLxYhUFBAARKiwB8PUrLnugNCERNwGQZYdRtyUL54JUFK6GNWSWQZlexrDKFWM2q5dBuECgBAnxR99rAXHdGpQQQoAsgAAIWCHgtgbRQbEkWAbFakmZFp1JAANehaIwAsRoNV5QKAiCgQcAr5L3OUiCNopEEBECgnxFQ1xMsvytseIjVfnZioLuxEPCLFxJL5SVeCcRqiRsY3QMBEAABEAABEAABEAABEACBLBKAWM2i1dBmEAABEAABEAABEAABEAABEChxAhCrJW5gdA8EQAAEQAAEQAAEQAAEQAAEskgAYjWLVkObQQAEQAAEQAAEQAAEQAAEQKDECUCslriB0T0QAAEQAAEQAAEQAAEQAAEQyCIBiNUsWg1tBgEQAAEQAAEQAAEQAAEQAIESJwCxWuIGRvdAAARAAARAAARAAARAAARAIIsEIFazaDW0GQRAAARAAARAAARAAARAAARKnADEaokbGN0DARAAARAAARAAARAAARAAgSwSgFjNotXQZhAAARAAARAAARAAARAAARAocQIQqyVuYHQPBEAABEAABEAABEAABEAABLJIAGI1i1ZDm0EABEAABEAABEAABEAABECgxAlArJa4gdE9EAABEAABEAABEAABEAABEMgiAYjVLFoNbQYBEAABEAABEAABEAABEACBEicAsVriBkb3QAAEQAAEQAAEQAAEQAAEQCCLBCBWs2g1tBkEQAAEQAAEQAAEQAAEQAAESpwAxGqJGxjdAwEQAAEQAAEQAAEQAAEQAIEsEoBYzaLV0GYQAAEQAAEQAAEQAAEQAAEQKHECEKslbmB0DwRAAARAAARAAARAAARAAASySABiNYtWQ5tBAARAAARAAARAAARAAARAoMQJQKyWuIHRPRAAARAAARAAARAAARAAARDIIgGI1SxaDW0GARAAARAAARAAARAAARAAgRInALFa4gZG90AABEAABEAABEAABEAABEAgiwQgVrNoNbQZBEAABEAABEAABEAABEAABEqcAMRqiRsY3QMBEAABEAABEAABEAABEACBLBKAWLVgteYdHdTc0mmhJBQRFYHysjIaObSG1m5qjaoKlGuJQFVlOQ0ZWEUbmtoslYhioiJQW11BA2oqaFNze1RVoFxLBAbWVlJlRRk1be+wVCKKiYrAYHH96+zqoe2tuK+IirGtchvrq2lHWxe1tnfZKhLlRERgWEMNbWvpoLaObqMaxg6rM0qPxPYJQKxaYAqxagFixEVArEYM2GLxEKsWYUZcFMRqxIAtFg+xahFmxEVBrEYM2GLxEKsWYUZcFMRqxIAjLB5i1QJciFULECMuAmI1YsAWi4dYtQgz4qIgViMGbLF4iFWLMCMuCmI1YsAWi4dYtQgz4qIgViMGHGHxEKsh4f7gBz+gy799pfYy4OatTfTovLk0+7wvh6w5l33uH26lk2bNpvqGwcbl/Wv+izLPQTM/YZw3bD8eeWAuHTzzcBozbqJx3cuXvUvvLnmLPvXZ07Xz6orVNatX0qvzX6KTT5+tXbZuwo0fradnn/gznTF7jm6WyNKFsb3tRrnHcNbF6h2/uYHOOe8Sqq6psY0qkfJ+86v/oou+9h3PuiFWvU3S3tZGd//hFjr/om8kYjOvSrMgVtN0jUzScKUkVoN8XyfJ3rRuiFVTYsmlDyJW+fvvmmuuSa7RqFkS6Pdi9UuXX0ev/Gtxn+Gw8Jnb+/x/2pyraOny1fK9aZPH0YO3X5v/HGI1mOiGWE32CgSxGh1/iNXo2GalZIjVYJaCWM1xg1gNNn6SyAWxmgT1YHVCrAbjloZc/V6sHjnrMnp+3o15W3zvx7fRC/MX5N9jMbtx09a8QGXhOqyxgX5//XdlHohViFXdEzlNN2IQq7pWM08HsWrOrNRyQKwGs2iarpHBemAnF8SqHY5xlAKxGgdlO3VArNrhmEQp/V6suqEvWLyMZl/6Q5p789W07/SpxGL2W5ecTbNOPEImnffYC/TzW+7pI3Dhs5rE0DWrU3cZsFmpSB0FgawvA46CSVrLxDLgtFpm13ZlYRlwdmhG29JSEqvRkkq+dIjV5G2g24IgYpXLRjRgXcLRpYNYdbG98Xf3070PPy3FqFu4clKv9yBWoxugtkqGWLVFMvpyIFajZ2yrBohVWySjLwdiNXrGtmqAWLVFMvpyIFajZ2yrhqyIVfcKTlv9j6McL41ko16IVQdFBfnaKy6UM6laYrVpOXUuvpda9vumDXugjIgIlIl9VgfWVog9trBvXUSIrRVbXl5GLIJ2YI9Ba0yjKqiyopyqKsuoRewziCPdBPghEJ9bbdgPMt2GEq2rEde/7u4e6ug02w8y9R0rwQbWiX2mOzp7xL64sFXazVtXU0ntnV3UJfYwNjnqB1SZJC+a1itWTuOQ+vyKzSTEKq8aveqnvyWlf4p2wicBxGpQcpr5FOBLvnAqXXbBGTKXjlhln9VrBn2fOid+mlqP+x/qqR1esMatTU107z130YUXf1WzZYWT/fa2m+iss8+lhsHm0YBfful5Wfhhhx9p3Jaw/bhv7h/p0E8cQRMmTDKue+m779CihW/SqbM+p51XV6yuWrWC/v7iC3Tm7M9rl62bcP36dfT4Xx6hc8+/QDdLZOnC2N52o9xjOOti9aZfXk8XfvsM6GwAACAASURBVPkrVFNTaxtVIuVd/7Mfy4jnXgfEqrdJ2tpa6be3/pq++vXLE7GZV6VZEKtpukYmabhSEqtBvq+TZG9aN8SqKbHk0gcRq/z9ZzMa8Ixj5pBTmCoaLGBHDR9KP7nyYkpCrNqyCsSqLZIe5agnCspP1ZnEy2eVnz6oiMEsVv9ztz9Txbr51D7qCNo+/VJqmey/pUrYLV/czcfWNXoDQ3cZMLau0eNpMxW2rrFJ035Z2LrGnCkCLJkz4xwIsJTjVkrLgLF1TbBzAbnsEwiyDNjm1jUsSN9d9kGfmDdevVRilT9Tu5X4CVznbiZODcPa5YiZ+8qAsZu2NMtqeDJuwriRcgZVHSqPl8h0zwCrybxCu6hArNoft7JEr4BJzqp0ogF/+7zjqey1X1Ld+/dS14AxQrB+hbbv/RXqqajbpdUQqzkk2LomogGtWSyiAWuCCpAM0YADQCuxLBCrwQwKsQqxGmzkJJcLPqvJsTetOWmxyrOqp55wuJw9LXSo7TKdKz1ZfO4+dXx+JxK3NuF4O7fc+VB+Io3Ts0hVYlR97l5uzO3g7TjdItMtrPnzX/zmPlk/f/bNi86UQWj54Pb6lWNqI7/0/dpnVRnHC45z3XahfVY5LwdY2r55LQ1c/GsatOjXVNaxlXbs8UUpWjuGzrBlK5QTgoDuzGqIKpDVEgEEWLIEMoZiEGApBsiWqkCAJUsgYyimlGZWY8CVaBUQq4niN6o8iFjlCmxEA3bHxCkmVp1bZHJa3lZz0TsrPIWlKosF6lmnHCtdGdXMqhLGXjOezq06nZ9zebwriq7/arHAtEZG8kncr8WqDYBKrDb3Bu4ZsPROGrjwJqra/Ca1jT5azrC2TjzFVlUoJyABiNWA4BLIBrGaAPSAVUKsBgSXQDaI1QSgB6wSYjUguASyQawmAD1glaUiVpXrohcGNRvrJ1adAtRPZL634kO5VFi5O3rVo2ZunZ9xeiwDDjg448jm3rqmet0Lcoa1dsU86ho0IbcseK9LxbLg6jiagzo8CECsZmdYQKxmx1YQq9mxFcRqdmwFsZodW0GsZsdWSYpVpmSyDLjQzKoSq8XEJPusumdWbYhV7sfHD5qeX5LsXIIMsZrS84EDLHGkTDWzqppZseNDGigE68BFN1FZdxtt3/NCOcu6uWwMPTpvLs0+78tWeoQAS3oYdcUqAizp8bSZCgGWbNK0XxYCLJkzhc+qOTPOAZ/VHLdSEqsIsBTsXEAu+wSCiNU4Ayzxsly/aMBey4ALLdMNM7PK5P2WAXsJZYhV+2PVeol+YlVVNPDt30lf1soti6lt7HG0ZtIFdP8r6yBWH5hLB888nMaMm2hskyBffhCrfTEjwJLxsNPOgABL2qhKNiHEajDTQqxCrAYbOcnlwsxqcuxNa05arHJ7vbauUQJQBV/y2rrGKVa5HBWR1zm7ymk+ftDeNOvEI3x9VnVmVjlwErdh05at+cjFKsASB1ZyC1nuEx9YBmw6ImNMX0ysclMqm9+joU9/nqo2vSlmVkfR7Z1fo7Mv+JaVVmJmVQ8jxCrEqt5ICZ8KYjU8w6yXALEazIIQqxCrwUZOcrkgVpNjb1pzGsSqU2g62+8O6lpoGbDKV2gLmTAzq84ov0uXr843U7WRRfFDj7+Uf5/9ZFUkYiwDNh2VMaZ3+6x6VV3W1UqDX/4GcQAmPngv1qbDf03d1YNjbGn/rUpXrPZfQunpOXxW02OLYi2Bz2oxQun5HD6r6bFFsZaU0jLgYn3N+ucQq9mxYBCxyr2zEQ04O5TS2VJEA7ZgFx2xqqqpe+9uGvzK5VTe3kRdA8fRlk/cKpYHf9JCK1BEIQIQq9kZHxCr2bEVxGp2bAWxmh1bQaxmx1YQq9mxFcRqdmzlbinEqgXbmYhVrq563Us0cMltVPf+vdRTVS+CL10kf7vqJ1loDYrwIgCxmp1xAbGaHVtBrGbHVhCr2bEVxGp2bAWxmh1bQaxmx1YQq5ZtpeOz6qyyeWuTjAb877Nn04Alt9LAt39DFdtXU9v4E4RgvZhaJ3zWqIXwWdXDpStWEQ1Yj6fNVIgGbJOm/bIQDdicKXxWzZlxDvis5riVklgNEhAx2OhJJhfEajLcg9QaRKzajAYcpM3IkyOAmdWQIyGoWFVb19SuelTMst5KNaufoK4B48SerBfTDiFadX1ZIVb1DAix2pcTogHrjZsgqRBgKQi10soDsRrMnhCrEKvBRk5yuSBWk2NvWjPEqimx9KSHWA1pi7Bilauv2Pa+2N7mNjnLWta5g1qmnCVE60XUPvITRVsHsVoUkUwAsQqxqjdSwqeCWA3PMOslQKwGsyDEKsRqsJGTXC6I1eTYm9YMsWpKLD3pIVYt2MLUZ9WvSo4UzL6sVR+9Sp1D95bLgtmXlcrKLLSyfxehK1b7N6V09B4+q+mwg04r4LOqQykdaeCzmg476LSilJYB6/Q3y2kgVrNjvSBilXuHaMDJ2xhi1YINbIlVbkrVxn9JwTrg3T/Ilm3f80K5LLijcR8LLe2/RUCsZsf2EKvZsRXEanZsBbGaHVtBrGbHVhCr2bEVxGp2bOVuKcSqBdvZFKvcnLLuVrkseIAQrZXNy8Ry4MOkL2vLlLMttLZ/FgGxmh27Q6xmx1YQq9mxFcRqdmwFsZodW0GsZsdWEKvZsRXEqmVb2fBZ9WtSzYdPyFnW2pWPiIBLDXKGVW5xM2hCPgt8VvUMqitWEQ1Yj6fNVIgGbJOm/bIQDdicKXxWzZlxDvis5riVklhFNOBg5wJy2ScQRKwiGrB9OwQpETOrQag58kQpVrmaipY1NICDL4mIweXtW6h1/Im0Yy+xxY145QNiVc+AEKt9OSEasN64CZIKAZaCUCutPBCrwewJsQqxGmzkJJcLM6vJsTetGWK1MLEZx8yhaZPH0YO3X2uKNvL0EKshEUctVlXz6t6/Ty4Nrl7/InUNnEDb97pIzrT+79y76aRZs6m+YbBxT8IIFrVfrNqCx7TyRx6YSwfPPJzGjJtompWCPKmFWIVYNR5oATNArAYEV0LZIFaDGRNiFWI12MhJLhfEanLsTWuGWPUnduPv7qcnnn+VNm3ZSr/+yTdp3+lTTfFGmh5i1QJe2z6rfk2q3LxIzrDyFjd8tEydLUUr+7TiKExAV6yCY/IE4LOavA10WwCfVV1SyaeDz2ryNtBtQSktA9btc1bTQaxmx3JBxCr3rj9EAz5tzlV0/JEH02sL36VRw4fST668OFWGhVi1YI64xKpqat17d9PgVy4Xy4KbhC/rYGo67Jci+NKZFnpSukVArGbHthCr2bEVxGp2bAWxmh1bQaxmx1YQq9mxVWrE6vrXidq2xA9u5AFENUN2qXfB4mU0+9If0tybr6b3VnxIP7/lHnp+3o3xt69AjRCrFswRt1jlJlc0r6ChL1xA1etekj1oG/tJavr4z6lz8J4WelR6RUCsZsemEKvZsRXEanZsBbGaHVtBrGbHVhCr2bFVasTqvccSrXomfnBnPU004Zhd6lVLgJWvKvuusnBN01JgiNWQwyUun1XPZnZ30j3/cwOdX/t7aux4n6i8krbN+A9q3v871FM5qGjP4LO6KyJEAy46bKwnQDRg60itFohowOY44bNqzoxzwGc1x62UxGqQGBPBRk8yuSBWk+EepNYgYjWSaMDPfJOIZ1fjPo75BRHPrroOtQT4sgvOkJ986fLrUrcUGGI15GBJVKyKtvON/qxjD6KRa/9EA5beReVtm4UP66G0Y9q51LLbudRTUe3bQ4hViNWDZn4i5BkQPjvEaniGUZYAsWpOF2LVnBnE6k5mEKvBxk8SuSBWk6AerM7UiNVgzY8kl1oC7C68cUh9qpYCQ6yGNH8axKqKBly76lEa8O6dYl/Wh2SvWqacJUVr27jjPXsJsQqxCrEa8gLgkR3RgO0zzVqJEKvBLIaZ1Rw3iNVg4yeJXFkVq2u2raZOsTqv2NHW1Uobdqwrlkx+vqllI23v2KaV9oPmlVrptndsp02tG7XSbtixllo723zTVlWWUVdXD3WIfnP/dY4Ltn6RrrnmGp2kmUzjXgKsOsFLga+94kKadeIRqegXxKoFMyThs+rX7LKOrUKw3iVmWe+kqk1vUFfdaGrZ/QtStHY27G6ht9ksAj6r2bEbfFazYyv4rGbHVvBZzY6tSkmsmlLv6umi7u4ukj/itbunW77u/D/3vkwnfvk1n673f5m/933OJ/+X73W78vTmLcvVIcuTdYl0jjLy9eTrypXFdVRXipghnR3U1tEpcu1smyxP/p9Lly+vt92yX6o8v3TONjjS7Cy7l4MXL8VC1tdtagak7yVwyYH/QTefKpbPluhx5KzL6KxTjiW1BFh1k5cC8/H767+bip5DrFowQ5rEqupO1aa3qO69O+VMa3n7FrE0+ONCsH6BWqZ9nnrKayz0OltFQKxmx14Qq9mxFcRqdmwFsZodW8UtVre2NVGTIzopzzp19eRm3Ximimes1MGzXDvEbJc6VjlmyNo6W2m9YxZuc+tHYqZtZ1qnBbg+rhdHMgTGDBpHlSLOSbGjpqKWRgwYVSyZ/LyxbhgNrCoeL4XTTqifqFXmgKqB1Fg7TCvtiAGjqbbS//62YUA1tbSLhwpd5cT91zkGi+i5e40drZMUaSIkALFqAW4axarqlufS4N3F0uCx3kuDLeBIZREQq6k0i2ejIFazYyuI1ezYCmI1G7ZisVddLWbMxHLF5tbWPssVm4S427qLqOySHesUM2jOpY2cjtOrwylA3aIy7WTSLKzqaiqoo7OHOru6yaawUjapKKs0ElYNNYPTbs7E2hfEZ5Ub2x/2WU3MKJoVQ6xqgvJLliafVb82lrdvFbOsYmmwmGV1Lg1+bvvB1F0zjIL4LTZvbaJH582l2ed9ORDBRx6YSwfPPJzGjNN7uuasJEh0QV2ximjAgcwZKhMCLIXCF3lmBFgyRwyfVXNmnCNtPqurtq7Id4RnCHmmUB1u/zynD57bz87pS8czlk5RyT5+XFaSBwscnkFSB4tDFkl88OybcxZqsEjbIPZ3V8fogTxDVyH/ramspZGOWbihtcNp25ot9OGy9+nAY/r6vnF9pSCssuqzmuR4S6ruIGI1kmjASQHIcL0QqyGNlwWxqrpYtXkB1YmIwezPWi6euj5VdS51NO5H+336IuOlwRCr5gMnTTdiYYJrmfe8cA6IVdtE7ZYHsWrOE2LVnJmfWG3tbKG2rjaxHLWFWkWwF54VzL+qzzrEZ93iffEZf85BYfrkU3lU/j5pWqmlozdPl6hLLHvl/HH7+VVWVFFteS3VVNTQgOoB8rVaLMHkZZi1lXVyeaP8u6JOisJa8Stfxf+1Im1tlXif0ws3H/W3zMefidd8nnyZnLeWqgrsGBDMin1zBXm4bKPeuMqAWI2LdPh6IFbDM0yqBIjVkOSzJFZVV2tXPkIDxEzrK+9ulW/N3HNkwajBXoggVs0HDsSqNzOIVfOxFGcOiFVz2hCrIuiMEHzrt6+Tvo8q8qiaUVSzkMrfUX3e2dxOh26fSfMGPpjPY07fXg7n8lO3757bP8/pg+deDtpYN1z48g3MN2x8/aT83+zjx2U5j7h9Vu0R27UkiNUo6aJsEwIQqya00pU2EbHK0ac2bWn2JLHwmdvTRUijNWn2WfVrvlwaLGZYeW/W/NJgETFYRg0evIdGr7OVRHcZcLZ6VZqthc9qduwKn9Xs2MrUZ1UJSN62QgpPfhUzkirAztZ2DsqT86Hk1x28nFUE34liWat7KSovQ21wLFsdy8tWe4PF8IykcylqvUjLS1fVYSJAk7JuKYnVpBjGVS9mVuMiHb6eIGKVa4XPanj2YUuIXayeNucqGtbYkJpwyGEBcv4silXVb14azL6s7NPKS4PbR8zMb3VTSlGDIVZtjPR4yoBYjYezjVogVm1QtF+Giu6q/Cz5dXvnZmrvbqNVWz7MB+PpEvsNfti73+AHzTn/TKefZpiWOQWmmkmc0JCbUZSCsayClL+jiuLJPo5qBnLkwFFyWWt/PCBWs2N1iNXs2ApiNTu2crc0drGato1mbZguy2JV9V8uDRYzrbUrH5ZvtUw5U2510zauNKIGQ6zaGOnxlAGxGg9nG7VArNqguLMM9pPc1t4st/TY1uF4FbOY6v1m8Xmz+F++inTO/7eJFTNb27ZKURrmGFA5kOprGqi+Wvzya1XvK//f+9sgXgdV18tgO/zqTqu7hUWYdpZqXojV7FgWYjU7toJYzY6tIFZ9bLVg8TKafekPae7NV9O+06fmU8177AW66qe/3SWXWq6cRZ9V1Rl3kB330uDuulFSsO4QW910NvRdGgyfVfOTHj6r3szgs2o+luLMAZ9Vc9rNO5rovjt/R9NP/lifpbIsMlWkWLVclpfY8lLbTjHL6YwSa17rzhzKD1L5WXJwnjH1Y0TAnkoaVjNGJhzfu8+hc7ZTbZOhs/9imPb55U3TNTKK/umWWUpiFT6rulZHuqgJBBGriAYctVX0yo99ZpWXAR9/5MF02QVn6LUwhlROH1ovsfrzW+6h5+fd6NmSUhKrqoO+S4N3+zz19C7Lglg1H5hpuhFDNGBz++nmuOM3N9A5511C1TX+m5PrlpWGdBCrfa3A24yw/yYLS/5du1288m/v//zZjpbt9B/i56fiJ8jhFpDsg6m2BWGfzKG1w+RyWRW4h/03eYbTK1iPqt/UZzVIu8PmSdM1MmxfwuSHWA1DL968mFmNl3eY2iBWw9BLNm/sYpVnKguJv6RwFJpZ7W9iVdnAc2nwbudS2/hPEcSq+UhN040YxKq5/XRzQKzqkkpfOqfgzAnRD6UIZb9O9RnPhBY7xtSOpS+1fZGeGvt8H1HZKERmdW8AICU+1Wwml6lEarHyg3wOsRqEWjJ5IFaT4R6kVojVINSSyQOxuit3pX3cn1x7xYU068S+eyMnY7VcrbGLVfZZLXQkFQ3YZBmwu42l4LPqZ5NysWxN7s0qgjBVbX6T8kuDp31eRA3eM8mxa1Q3fFaNcCWaGD6rieI3qjwrPqubWzfRuu1rxXLbtbRO/K4Xf6+Xr+vk3+o99hMtdPDM5QgxyzlqwGgaOWh07nXgzteRA8bQKPH/0NpGI45xJM6CWI2DQxbqKCWxmgXeYdoIsRqGXrx5g4hVbmEpRwP20j7f+/Ft9ML8Bb4rSuO1WkJiNYlO6tTpJ1bdeb90+XW0cdNWevD2a3WKLZ00IlIwPf89ogXCf1f4VlH9eKKjfka01+zS6SN6AgIgkCkCHJBoTfMaMQO6Vr6u2SZ++f9t4n/1d+9n7V3tBfs2tHao9OscM0j8FnitF8GEcIAACIAACIBA1gl4aR8VqyepyUMvprHPrKbVsLpiVaVzGrGUZ1bd9qr66F80+JXLqXrDfPlRR+N+tHXmf1Hb6KPSalrZLsyspto8fRqHmdXs2CrKmdXWrhY5GypnQeVMqJgF5ZlR+fea/Kzohpb1RYGNqBuZmwEVvyMHjhH7cPb+zTOjYqZU/i9mSkt5qxTMrBYdJqlJgJnV1JiiaEMws1oUUWoSpGVm9fW1r9OWVjEBFPNxwOgDaEjtkD61emkfnpTj4/fXfzfmFvpXl4hY9Yqwm/T6aF2x6n7iUIoBlnRGZ+cbv6MHXlpO36j7uUzOYrX5wKupfdThOtnpkQfm0sEzD6cx4yZqpXcmChJdUFesrlm9kl6d/xKdfLr9GWP4rHqbGtGAjU+BWDNEEWCpSazUWLd9jYyCq5bg8usGsSw3J0ZZlK4R27JsLdhXjlrLIlSJzxEigrkSpPLVIUrLy8pj49be1kZ3/+EWOv+ib8RWZ7GKsiBW03SNLMYzys9LSawG+b6Okq3tsiFWbRONrrwgYjWKaMDH3nEsPbP8meg66lPy0+c/TcdMPsZTrLqzXPKFU1MVCDd2sXrj7+6nW+58qM8WMUooJgnHT6xypGBnJGCOZjyssSH/xKG/ilUZYOn+u+jCj7VR3bL7qGrTGyJScDW1TD5T7tHaNv6EgicixGrs16k+FSLAUnT8+3uAJTkL2tLrCyoEp3tWdIP4bJ1YplssUFFtRV2v8OzrF5rzExW+ob1+o8Nqh0dnzIAlQ6wGAwexmuMGsRps/CSRC2I1CerB6kyLWP3mX79JPLsa9/GLT/+CeHbVeWAZsI8VWPyddcqxuyh2FrH3Pvx0Ig69zq1ruNmNQ+rz7WBxunT56nxvPn7Q9D5T4/1arM6bS7PP+zJVNr1Dde//Sf5WNi2hbrFZfOvkz1HLVCFaRx/tORIgVuO+TPWtD2I1Ov6lKlbbu9ukyFRLcXlWdGOb+BXic+WW1blZ0t5Z0WJ068U1YpRrOS4vw3XOhPLnnC6LB8RqMKtBrEKsBhs5yeWCWE2OvWnNaRGrpu2OMr3fRB0Hw3Vv5RllO4qVHfvMKgPwWvKbRofeYvDU5/3JZ7UQE96fNSda76OK5uUicvBIMdMqRKuYaW0f+XFdnJGk010GHEnlKNSIAHxWjXBZTdwpgqct3fwOLW96T7y+Te9teZdWbV1BHzSvkFu38OfFDvb7HDNonPzlrVjYJ3T0wLH59/h//gxHvASysAw4XiLpra2UZlbTS9lOyyBW7XCMo5QgYpXb1d+iAasVsP06wFIaZ1bDniQQq30JVm/4hxSstUK4VrSspa5BE6VgbZnyORGQaf+wuAPlh1gNhC2RTBCr0WPnYEXvCTHKwvS9Le/Q+1tYnOZEaqGjoWYwjRVCc4RYgsuCc8qQyTRpyHhqqMz9z0K0sW5Y9B1ADcYEIFaNkSWWAWI1MfTGFUOsGiNLLAPE6q7o/fZZTZNQ5VbHPrOaVp/VMGcPxKo3vep1z0t/1rrlf6JyEVClc/AeOdE6+d+oc8heYZAb54VYNUaWWAaIVTvo2SdUClApRHOzpPz6vhCkW9uafCvh2dBpQ/egyYN3o92G5F4nNEyUYpT3GXUeUUYDtkMBpSgCEKvZGQsQq9mxFcRqdmwFsZodW7lbGrtY5QakMRpwUBPCZ/XLRdFVbXqT6l+/lmpXPizTdlcPpv/p/j90wNGn0ZgJ04rmdycIEl1QV6wiGrCxOUJnQDTg4AhZkK7culz+rmh6X7y+n3sVy/BXNi2nHR3bPQsfXz9RCFAxKzp4Ck0Sr/LvhinydcSAkX3yRBENOHiPs5ETPqvB7ASf1Ry3UhKrQb6vg42eZHJBrCbDPUitQcRqFNGAg7S9v+dJRKyWEnSI1eJiVdmb92ht+OeVVLP2ObqjZQ4dVf8GNR54Du3Y60IRSbhWe1gE+fKDWO2LFwGWtIebccIoAizxNi5KjEphKgVqTpjyZ17H4JohOQE6OCdEWZhOqJ/U+zqZKsortPoGsaqFqU8iiFVzZpwDYhViNdjISS4XxGpy7E1rhlg1JZae9BCrIW0BsaovVhVqFqsPP/IoHVP2ME2uEIGYaofR9r0uoe17f03OuhY7IFaLESr+OcRqcUZBUwQVqzwLulIEM8rPjvaK0VUiuBGLUq+tXnif0Qn1PEPaOzsqBOlEIUwn9s6Uso9p2ANi1ZwgxKo5M4jVncwwsxps/CSRC2I1CerB6oRYDcYtDbliE6scBZj3UeU9VgsdaXPq1TESfFZ1KPVNU9bRTHUrHqRa/l31qPyQfVpbJ54mfFpnUcewA80LLZBDd2bVaqUoLBCBUvdZ/bD5AzEz+r6Ists7QyqEKP+/UrxuaFnvyYwDF+WX6sqZ0slSkLIwHVc/IRBnG5ngs2qDYjxlwGc1Hs42aiklsWqDR5rLgFhNs3X6ti2IWOUSSjkacFasF5tYzQqQIO2EWA1CrTeP2ApDitaVOeFa1t2Rix486TQpXNtHHR6i8J1ZIVatYIylkFIQq9vEwxi3H+kq4Ue6QviRsjDt7OrYhWVd5YDcDKmcKc35j/JrbqZ0EvHnaTsgVtNmEf/2QKxmx1YQq9mxFcRqdmwFsZodW7lbGrtY9dtnlaME3/vw0/T8vBszRxNi1Y7Jalf9WQpWFq9lHVvF8uAROdE6aRa1jf1kqEogVkPhizVz1sQqL89dsnERzf/wJVqw4TV6Y/2/5N6kXst2GSQvzZ0iIuxOG7qniLa7u3ydPIT/34N4j9IsHRCr2bEWxGp2bAWxmh1bQaxmx1YQq9mxVWrFqooQnLVlwPBZNfdZ5UH4yANz6eCZh9OYcRN3OXtqPnwqL1rLW9dTT1VDr2gVwnXCZwk+q+EvOPBZDc6Q9yJ9de18KUxZoLJQdQrTK+gK+lX5r2j04HFi25epOTEqXvds3FuKUl7Sm6UDPqvm1oLPqjkzzoEASzlupSRWg3xfBxs9yeSCWE2Ge5Bag4hVRAMOQtp+ntTMrH7vx7fRC/MXZG5mFWLVvlhVw7x63ctyeTDPtFZsW0E95VVylnVhzfG0cNNA+tRnz9A+I3RnVrF1jTZSawnTvHUN70c6f81L9KaYLeUZU/7ba4/SGcP3o/1HHUQHjpxJbS9uprO/cDENqBtojVGSBUGsmtOHWDVnBrG6kxnEarDxk0QuiNUkqAerE2I1GLc05IpFrHrtq+rV+WuvuJBmnXhEGrhotwFiNTqxqoxQ3rKOBi65lQYuvpnK25toSede9FrFsXTisUdQy8RTiERE1GIHxGpfQphZ3XXEfNSyQcySvkWLP1pISzaJ34/E3xsX7rKcd8/G6bTXsH1o+nDxK173GjaDeN9SdQSNBlxsDCf1OcSqOXmIVXNmEKsQq8FGTbK5IFaT5W9SO8SqCa10pY1FrDq77Oezmi4sZq2Bz6oZr6CpWagOePs3NGjRTcQClo+ugeNo+/Sv0o49vlhw2xtdsRq0bchnj0AcPqvtXW1SiEpxKl9z4tQdjXf0wLFClM6gvRp7xakQqHs1zqDysnJ7Hc5wSfBZzY7x4LOaHVuV0sxqdqgHaynEajBuSeQKIla5nYgGnIS1+tYZu1hNvsv2WwCxap9pvXrJBQAAIABJREFUoRLLOreJmdbfytlWXh7MR48ITLNjjzlir9avUmf9brtkh1iN10ZhaotCrL7ftDQ3YyrEKQtTFqjLtrzbp5kDqgbKWdLpQozKmdMR4leI1CG1Q8N0p6TzQqxmx7wQq9mxFcRqdmwFsZodW0GsZsdW7pZCrFqwHcSqBYgBiihv3US1Hzwq9mn9M9WIX972hoMxcRCm1gknia1vTpIilg+I1QCAE8oSVqzqLuflwEdSnIpftax3ktgmBoc+AYhVfVZJp4RYTdoC+vVDrOqzSjolxGrSFtCvH2JVn1XaUsYuVhcsXkazL/2hLwdEAzYbIu7gNCa5w/gtNm9tokfnzaXZ50Xvs+ruk1d0wZoPHheiVQjXD/5MFdtXyyxt446n1vE54dojfApHDq2htZtaCyJCgCWTEWQnbZgAS+3dYjlv74xpoeW8HIGXBeleYkkvz5xOH76vFKmVImiX7QM+q7aJZq88+KwGsxmiAee4lZJYRTTgYOcCctknEESsIhqwfTsEKTF2sXrkrMvoiJn70scP2pt+fss9+ei/p825io4/8mC67AL9CK9BOmw7DwIspUOsKrtWbXw9J1rFTGvVxtfk2x2N+1HbxJNpwIzTaW3VDIhVQSDMgwrb55CJWNVZzst7lUo/UxanLEx59nTEvjSsdrjtpnuWB7EaC+ZUVwKxGsw8EKsQq8FGTnK5MLOaHHvTmiFWTYmlJ33sYlUFWNpt0lj6yvd+kRerHDHYKV7Tg6hwSyBW0yVWlbUqdnyYWx68UgjX1X/NvT1IBGMa9xkx03oytY0/wdOwmFmN/8zzE6uL137QNzpvb6Re556m3NopQ6bll/KqZb1Th+wef0d6a4RYTQx9aiqGWA1mCohViNVgIye5XBCrybE3rRli1ZRYetInJlZ5ixoWrmrZr9reJmvLgNmU8FlNz4DepSXdnVK01onlwfxLLR/J/VrbhD9rbonwydRdMyTFHegfTXMu531n0yJ6d8tCenPdAtqwY30fAI11w4QwzS3hzQnT3NYxtZV1/QNUynoJn9WUGaRAc+Czmh1bldIy4OxQD9ZSiNVg3JLIFUSscjsRDTgJa/WtM3axyst9995jEv3kyovJ+ff3fnwbvTB/QX6mNXk0+i2AWNVnlVRKGWCp5VXasXCeDMpUuXmxbEr7yE9Q6yQhXIVfa2dDcrNxSXFJst7lTe/R86uepmdW/o1eWv0cbW1r2qU5Exom0UGjZtKBow6hg0fPpBkj9iNe5osjHQQgVtNhB51WQKzqUEpHGojVdNhBpxUQqzqU0pEGYjUddgjSitjFqruRPLuqjrk3X037Tp8apB+J5oFYTRS/VuXOaMCVzctk9GC5THjNMzJ/5+A981GE20cdrlUmEpkRWL9jHb0sROkzK5+Qr6u25rYdUkdDzWDaf+TBdMjYj9PhEw6hPQcfQjyTiiO9BCBW02sbd8sgVrNjK4jV7NgKYjU7toJYzY6t3C1NXKxmF12u5fBZTafPqntc+W1dU/PB32jQkl8Rv/KxvHs3eppOpVnHHiC2vjklv/WNjXGaJn+sOAIsbe/YRvM/fJle+OBpMYP6FC386M0+GHnW9LBxR9G+y/ekgz51FB0w8WPy87Bb19iwVZgy4LMahl5p5IXPajA7pukaGawHdnKVklhFNGA7YwKlhCcQRKwiGnB47jZKiF2sqgBL7LNaCgfEarbFqhqDlU1v08Alt9L6RU/Rc62H0/l1t1N39WBqmXIm7dh9DnUMPyj0cE3TjVgUYnVzyyZ6ff2r9Pq6f8jXN8Sv0+d08uDdxHLej9EB4le9VpRVkEk04NBGiKEAiNUYIKe8CojVYAZK0zUyWA/s5IJYtcMxjlIwsxoHZTt1QKza4ZhEKRCrIalDrJaGWFXDYO3yxfTai3+h8+rvpapNO2cCO4fuTdunzaGW3c6h7tpgS1PTdCNmQ6x2dLXTa+v+KUUpv76+7lVasXVZ/owaMWCkXNabE6aHyNfBHsGsIFZDXoQizs5Pli/62nc8a8EyYG/4EKvBBmWarpHBemAnF8SqHY5xlAKxGgdlO3VArNrhmEQpsYvVrO6nWsg48FlNYuia1em3DLhQKTWrn6CaD58Ufq1PCOG6UCbtrhtFreM+RW1jj6M28dpd02jWkIynXrRxQU6grv2nfOX/1cEReVmQ7hSoH6Px9RONe5z1ZcDGHc5wBojV7BgPPqvZsVUpidXsUA/WUojVYNySyBVErHI7EQ04CWv1rTN2sbpg8bI++6smjyB8CyBWwzOMuoQgYlW1qWLHGiFahXD9QPwK4Vretll+1DHsQClYW4VwbR99ZNRdSKR8ninlGVPnDCrPqKpjnxH70wFi9jS3tPcQuY1M2ANiNSzB+PJDrMbHOmxNEKthCcaXH2I1PtZha4JYDUswvvwQq/Gxtl1T7GLVGf3XqzPYZ9W2iVEeEwgjVp0Eqzb+S2x98zc541q97gX5UY/YSkXNtLaNO44663fLLHT2MVXLet9YL5b2iiW+7IuqjkmDp9KBI8WyXhaoo4VAFX9XVlRZ7S/EqlWckRYGsRopXquFQ6xaxRlpYRCrkeK1WjjEqlWckRYGsRop3kgLj12sRtqbBAqHz2pp+ayuWb2SXp3/Ep18+uyCo6msuzU308ozrmK5MG+Hw0dn/dSccB2fWyrcU1GXLydN/ljKZ3X6QQeKmdOcKJWzp2IW9YPmlfk2D6sdng+ItL8QprzMd2id3aXP8FlN4MJlUCV8Vg1g9SaFz6o5M86RpmtksB7YyVVKYhXRgO2MCZQSnkAQsYpowOG52yghdrHqFw34xt/dT/c+LLa4mHejjX7FVgbEav8Uq84BJvdtZd/W3hnXsq4W+XH7qCOkYG0VwrVj2EGpuRHb1LKRHn5qLq3cupzub/kT8f6nzmPa0D3p2EmfoiPGHyu2ljmSBlYNivR8gliNFG/owiFWzRFCrJozg1jdyQxiNdj4SSIXZlaToB6sTojVYNzSkCs1YnXeYy/QVT/9LWVtGTDEKsSq80SuXv+inHGtFTOuVR+9Kj/qrhkqROun6IP6o+jxt7vp9HMujP3cZ4H62LKH6aGlf6KXVz9HR3Tnto56RvyMHDCKjpl0Ah029gj5yv/HeUCsxknbvC6IVXNmEKvmzCBWIVaDjZpkc0GsJsvfpHaIVRNa6UqbGrH6vR/fRi/MX5C5mVU2JwIspWtQe7XGls+qdk+7O8Vs69M0YOntVLvqL1TW1ZoTrmLbm9aJp1DL5M9R2+ijhDNtpXaRpgndArVTtEkdM8ceTp+ecjIdP/kzxDOpaTrgs5omaxRuC3xWs2Mr+Kxmx1alNLOaHerBWgqxGoxbErmCiFVuJ6IBJ2GtvnXGIlbVrGmx7l57xYU068TcjE/cB0cpnn3pD2nuzVfTvtOn9qmet9tZuny1fG/a5HH04O3X9vkcYjVua5nXF7tYdTSxcvMiOdOa2wrnifwnHY37iKXCRwrReqR8Dbp/q5PG25sWy5nTl1c/Ty9/8BxtbP1Iflxf3SCX9B4+7ig6VLzuO+IAc4gx5YBYjQm0hWogVi1AjKkIiNWYQFuoBmLVAsSYioBYjQm0hWogVi1ATKiIWMSqs29+PqsJ9V9We+Ssy2jTlmb5t1usfuny62jjpq15gcrCdVhjA/3++u/mmwyxmqT19OpOUqw6W1iz9jmqXvu8/K0RvyKWsPy4c+je1CYEK2+Bw+K1u3aEXsdEqiUbF+UE6odCoIpXnlHlY3DNEClQDxMClV9nDN9Pu8wkE0KsJknfrG6IVTNeSaaGWE2SvlndEKtmvJJMDbGaJH2zuiFWzXilKXXsYjVNnXe2xW9mlYXsty45Oz/jy7PEP7/lnvxyZfiswmdVd0y7I13y1jc1a4R45VchYqmnOydch0wXwvWIncK1blcf0iUbF9JLQpj+Xcyg8uvm1tz2MkNrG3cK1LFH0vTh+3g2T0UDPmjmJ3SbH1k6+KxGhtZKwfBZNccIn1VzZpwD0YBz3EpJrCIacLBzAbnsEwgiVhEN2L4dgpQIsdpLzUus6rwHsQqxqnviFboRk4GZWLiuVcK1KydcG/agtjG5GdfXq0bQCxsX09/FDCoL1C2tm2Ua3l7m0PG9S3yFQN1r2N5FmwSxWhRR4AR3/OYGOue8S6i6piZwGWnKCLFqbg2IVXNmEKs7mUGsBhs/SeTCzGoS1IPVCbEajFsaciUiVp3Lbt0QkooGrCNMua3udCxWr7jyP6m9IzcrVuzY2tRE995zF1148VeLJdX6/Le33URnnX0uNQwerJXemejll3gZKtFhhx9pnDdsP+6b+0c69BNH0IQJk4zrXvruO7Ro4Zt06qzPaectKyujgbUVtK1lZ5Ahr8yrVq2gv7/4Ap05+/PaZesmXL9+HT3+l0fo3PMvKJilYs2LVPHhs1QpROlbH4h4vR1d9Kxo9jPd5bS5d/Z1ZN1wOnLysfSJ8UfREROOEgJ1um4zZLowtjeqSCOxewyXl5cRLy/d0VrYVhpFJ5Lkpl9eTxd++StUU1ObSP22K73+Zz+my799pWexlRXlVFVZRi1tuYcrOHIE2tpa6be3/pq++vXLU4OEl9fzudXWnl5b6V4jUwM1oobUiOtfd3cPdXTq3VdE1AwrxQb5vrZScUyF1NVUCDv1UGdX9m0VE7LEqqmrqaT2zi7q6sq5X+kc/P13zTXX6CRFmggJxC5WvXw+I+yfdtFBxaq8MRFiQlesajcICa0S0BWrVisNWNib69+gF1Y+Sy+IAEnPr3iKmtq3yZLGlBEdI4IHH10u/KyHTKDdJxxHXeOOpq4xR1F3/YSAtaUvW9bFavqIRtciiNXo2NouOQti1Xafs1peKYnVrNpAt90Qq7qkkk8XRKxyq+sHVCXf+H7egtjFahoDLPEYMPFZde8HiwBL6T+L0hJgyY/UWx+9IaP3cpCkl8Rrc/tWmXT0wLG5CL5ime9RNfU0ffsyqhE+rtXCx7Wsq12m6aqf3Cc4U9cg89nqNFkQAZbSZI3CbUGApezYCgGWsmOrUloGnB3qwVqKZcDBuCWRK8gyYG4ntq5Jwlp964RY7eXhJ1YRDTj5QWqjBWkUqws2vJ4XqLzVjBKoYweN791i5ggZLGny4N12QVC9/hUpWGvWicjCa56nsu62nHAdNHGncBW+rl2DptjAF2sZEKux4g5VGcRqKHyxZoZYjRV3qMogVkPhizUzxGqsuENVBrEaCl+imWMXq7wM+PgjD6bLLjgj0Y47K3f70DYOqc9H++V0hfZZRYAlBFjSHcgcYOmxx/4fbdunvXcf1OdpW0duy6QJDZPoUBEc6bBxR8h9UCc16IvM6g3zdwpXMTNb1t2aE64Dx+f2ce0N0NRZv1P0IsCSrtXM0yHAkjmzUsuBAEvBLIpowDlupSRWEQ042LmAXPYJBBGriAZs3w5BSoxdrLq3fgnS6DTlgViFWC02HldtXUHPrnqC/vHuizRydSP9uufXMgsLUhamUqAKocqCNexRveEffWdcu1pywnXA2D7Cdf6S9fJ9bF0Tlviu+SFW7TPNWokQq8EsBrEKsRps5CSXCzOrybE3rRli1ZRYetLHLlbZZ7XQkVQ04KAmgViFWPUbO/9aO5/uXXIX3f/2XNresY1Gi58zys6gdXtsps/uNouOnfQpqqmILmJsxfbVVLv8AapbOU/s5fpSvpk9os6nqs+nzsF70P5Hn0HdHvu4Bj0fguTDPqtBqMWXB1vXmLOGWDVnxjkgViFWg42c5HJBrCbH3rRmiFVTYulJH7tYTU/X7bUEAZbssYyqpLh8VtkP9dmVT9Az4vdlsRcqH9UVNXS0iNx7jBCn/DplyLSouulbrp9w5Qydg/ektvGfEsuFj6U2sZ9rT+Wg2NvnrBA+q4niN6ocPqtGuBJNDJ/VRPEbVV5Ky4CNOp7BxBCr2TFaELHKvUOApeRtDLFqwQYQqxYgRlxElGJ1Q8t6enbF34RIfZKeWfU32tSyUfZm3xEH0DEThUCdeJxY6ntUxD3UL14J19o1fxNLhl+mss7c1jjyKBf7kI2YmROu4pf/5vfiPCBW46Qdri6I1XD84swNsRon7XB1QayG4xdnbojVOGmHqwtiNRy/JHMnIladAYuuveJCmnXiEcTLgz9+0HT6/fXfTZJHoLohVgNhizVTFGL1hQ+epmdWPCFnUhdtXCD7M6JuJB3NM6hCoB4z4VPUWDcs1n6aVlbeup6qP3qVqjb8k6o3/lO+lrdtzhfTVT+F2ocfTB3DPyZePyZfeyqqTasxSg+xaoQr0cQQq4niN6ocYtUIV6KJIVYTxW9UOcSqEa5EE0OsJoo/VOWxi1UWqsMaG6Qo5Si837rkbClWb/zd/XTvw0/3icIbqmcxZYbPav/yWX130xIRLOlJIVBzM6ldPV1ypB0x4dj8Ut+9h+3rOfrS5I/lFw24snkZVX0kRKsQsNVCuFZtfDW/nyt3qmPovjnhOuIQ+drRuE/oMw0+q6ERRloAfFbN8cJn1ZwZ50jTNTJYD+zkKiWximjAdsYESglPIIhYRTTg8NxtlBC7WOUZ1Lk3X037Tp/aR6xylOCrfvpbQoAlM7O6b/RNcofZvqR5axM9Om8uzT6v9MUq73/Ks6cc0ZcF6urmVRLztKF75mZQxVLfo4QvamWR5bJpuhHTtX3VpjdzwvWjf4hZVzEDuzk3g8xHT3lNr3DNzbi2j/iY2NfVPKIxxKrJWRt/WohVc+YQq+bMIFZ3MoNYDTZ+ksiFmdUkqAerE2I1GLc05IpdrPJs6q9/8s1dxCpmVoMNB4hVPW66y4DXrF5Jr85/iU4+fTb9c83fhQ9qbpkvR/blY1B1vZxBlUt9xx9H4xsm6jVApMqiWHV2rqy7TQpWuWy4d/aVZ2LV0V0zXArW/MzrsIOpu7axKB+I1aKIEk0AsWqOH2LVnBnEKsRqsFGTbC6I1WT5m9QOsWpCK11pYxer3/vxbfTC/AVyua9aBrzbpLE0+9If0qknHE4/ufLidBHSaA18VjUgJZxEV6y2dbXSI0sfoLsX3ZGP5stNHzNoHP373nPozL3OtbIfasI4rFRf3rpJLhPmvV1zs6/C37V1Q77szvqpO2dehXDtEEK2p7yqaN3wWS2KKDUJ4LOaGlMUbQh8VosiSk2CUppZTQ3UiBoCsRoR2AiKDSJWuRmIBhyBMQyLjF2scvvUkl9nWy/5wql02QVnGDY/HckhVtNhh0KtKCZW12xbTbcvuJX+d9H/5KP58h6oJ049hc7Z+3w6csIn09/JhFtYsW2lFKx9/F07d+Rb1dG4f37mVfq7Dp3h2WKI1YQNaVA9xKoBrISTQqwmbACD6iFWDWAlnBRiNWEDGFQPsWoAK2VJExGrKWMQujkQq6ERRl6An1jl5b2/feMmevS9B6izu1O2Y+bYw+nUaZ+jf9vzHGqoGRx520q1gvKWdVSz9jmqWf2EfK3YtqJPVzuH7k1tYz+Z2yJn2EHUXTdKfg6xmp0RAbGaHVtBrGbHVhCr2bEVxGp2bAWxmh1buVsau1j90uXX0Sv/WrxLIKWsbl2DaMDZC7C0fsc6+tvyP9Pf3n9U/P5ZnhO8zPeEKSfRkfVHU9eKNumzavvIus9qWB6Vze9RzYdPU/Wap8XrU1Te3iSLvGHHf9D5tbdTff1Aua9r98iZVDPpcNpQM516KgeFrTb2/Hf85gY657xLqLqmJva6o6gQPqvmVOGzas6Mc6TpGhmsB3ZylZJYRTRgO2MCpYQnEESsIhpweO42SohdrLKf6lmnHLvLkl8EWApmTgRY0uPGM6tr2pbQ3DceoMeFUH193T9lxn1HHCBF6glTT6Z9hu9PzgBLeiXrp0rTjZhuNGD93pmnrN4wn6rFjOvvX9xO59fdTkN71vUtRERX7hiytxSwHcMPkfu98mxs2g+I1bRbKPr2QawGY5yma2SwHtjJBbFqh2McpWBmNQ7KduqAWLXDMYlSYherPIN67RUXyr1VnQe2rglmfojV4tw4mu8Ty/9CT678M63YkluKeqyI5vupyUKkCqHKs6rqgFgtztN2Ch7Dp3zm09TYuZwqNy+kmi3it2mRmGZ5i8ocPq9cb1f9FOnr2jF0H+ocwq8zxOtetpsUqjyI1VD4SiIzxGowM0Ks5rhBrAYbP0nkglhNgnqwOiFWg3FLQ67YxWqpzayyEeGzmoah3LcNvDfqX99/RC7z5V+O8jtY+J8qgcoitaqiOn0NR4vyPqubP1gsxWvV5reoUgjYqk0LqbLp7T6EeiqqekWrEK8sYoWA7RRCtmvAaJCMgQB8VmOAbKkK+KxaAhlDMaUkVmPAlWgVEKuJ4jeqPIhY5QoQDdgIcySJYxervNz3ljsfork3Xy33WuVjweJlcuuarEYEhliNZGwGKnTZlnfp8V5f1L9/+IIsY8rgafTpqSfR2fudTnvUzwxULjLFR8AvwFJZd0evaBXilUWsELCVm96iipa1fRrXXTs8N/PqEK8dw/YR2+bUxteJflITxGp2DA2xmh1bQaxmx1YQq9mxFcRqdmzlbmnsYpUb4LV1jdfS4KxghVhN3lIsTOUsqvBHfW/zO7JBh4w9TMykfpZOEMt99xw2nUYOraG1m1qTbyxaUJCASTTg8taNQrT2ileegZUzsQt3WT7c2TAtJ16FiJWzr0K8dtZPgyVCEoBYDQkwxuwQqzHCDlkVxGpIgDFmh1iNEXbIqiBWQwJMMHsiYjXB/lqvGtGAk4sGfPSJn6XH3xMRfYVAfXz5o7S1rYnKy8pzAZPk78k0tLZR2rzYPqtqYMBn1fopUrRAt9+1iVj1KryyednO5cO94nXX5cO1fcVroxCx4re7ZnjR9hZLAJ/VYoRK/3P4rAazMXxWc9xKSawiGnCwcwG57BMIIlYRDdi+HYKUCLEahJojD8Rq/GL13aUL6YVXn6C7uu+kJRtFIB5xVIrIsaftfiZdsP9Xaf+RB+1iVYjVvkjSEA1Ytci2WHUbX3v5sNjnVc68ctAmfhXilYM49ZRXGV0lIFaNcJVkYojVYGaFWIVYDTZyksuFmdXk2JvWDLFqSiw96RMRqxxkadOWZk8KC5+5PT10NFoCsRqfWOUgSf+78HZ6/J8P0dgdo2mu+GmsG0b/vvcXpUgdOWCUr8UgVvuvWPUaFPnlwyJo085lxB7Lhwfv2RtxWCwhHrav/Ltr0JSCVwaIVY0LZ4kngVgNZmCIVYjVYCMnuVwQq8mxN60ZYtWUWHrSxy5WT5tzFQ1rbKDfX//d9FAI2RL4rIYEWCQ7L+/930X/Q7e+dgOt35Hbi3OvYXvTnH0vobOmn0s1FcUD5+iK1Wh7gtJ1CIRdBqxTh1careXDlQN3itdGIV4beescXj48NGi1mc4Hn9XsmA8+q9mxVSktA84O9WAthVgNxi2JXEHEKrcT0YCTsFbfOmMXq377rCaPIngLIFaDsyuUc+FHb9Kf33uQ/vLePHp702KZlAMmfXa30+gz4re+ukG7YohVbVSJJ0xKrHoL2Peo6qN/UfVH/6Cq9fPF9jlvUpmY4XcePeJhScfwg6h9xEzxeoh87Rq4c+/exIFG2ACI1QjhWi4aYtUy0AiLg1iNEK7loiFWLQONsDiI1QjhRlw0xKoFwBCrFiA6inh17Su9IvVBWrH1faooq5DilEXqSbvNEv6pZj6EXDTEql0bRVlamsTqLv3s7pSCtXrD/JyAFULWHbyJ83RXD6auht2ok5cRD5lOXfXib/n/HsTitlQOiNXsWBJiNTu2gljNjq0gVrNjK4jV7NjK3dLYxSovAz7+yIPpsgvOyC41R8vhs2rPZ/XFD56RIpV/1+9YSwOrBuVF6qdFZF91BIkuqCtWEQ04/tMy6gBLUfeovL1JzrpWb3yVqjb8g254+wD6xoBfUG2Z9zZJXYMmiaXDe4itc1jA9r4O3Tu1s7EcDfGir33HEyPEqvfogs9qsLMOPqs5bqUkVoN8XwcbPcnkglhNhnuQWoOIVUQDDkLafp7YxSrvsfrzW+6h5+fdaL83CZQIsRperD654rFekTpPbj/DQZM+O1Us9Z06i46ZdPwuVg3y5Qex2hdjf4oGHPdlgQMsff7M02hA+2qq5JUB25ZTZfNy8Sr+5v9b1no2qVtEI+4cNFnMwk4RInaKEK9C1MpX8d6gCXF3I18fxKo5eohVc2acA2IVYjXYyEkuF8RqcuxNa4ZYNSWWnvSxi1X2WS10IBqw2eBwz0qZ5A4jWJq3NtGj8+bS7POCitW7qX1cDz25+XH687IHqb2rjcYMGieX+fKS30PHHuHbFYhVEyt7pw1j+/C19y0h6zOrbh6FogGXdbUIwSrEKwtXKWLFq/iVr9tWUFnnjl3xim2ZOlmwNkwWs7FCvApBq15Z2PKS4ygPiFVzuhCr5swgVncyw8xqsPGTRC6I1SSoB6sTYjUYtzTkil2spqHTttsAn1Uzohw46Zrnv0Mvr35OZuTlvv8+Yw595aBvFdx+xqyWvql1Z1bD1IG8dgik2mfVThc9SylvWddXvG5f0TszK4Ts9tWeebqrhzhmYlnQillZMSObe50snLUrImwxEZYBR4rXauHwWbWKM9LCSkmsRgoqBYVDrKbACJpNCCJWuWhEA9YEHGEyiFULcCFW9SAu3fw2/eTlq+mxZQ/LDA01g+mi/b9GF4pf/jvKA2I1Srp2y+6vYtWPYllXu1xCnFtKnFtOnFtezK/LxWzsNs+sXYMmimXFYia2npcWi9lYtcSYZ2NrR1gxGsSqFYyxFAKxGgtmK5VArFrBGEshEKuxYLZSCcSqFYyJFJKIWGW/1at++ts+Hb72igtp1on+Sz8ToaNZKcRqYVDLm96jG/75X3T/23dTp4imyvuifu3gb8UiUlXLIFY1B3MKkkGs6huhvO0jh3DNLS3OLy/ettKzoB6xV2wnC1gWr70zsPJVitnJIlpxnXYTOaH3AAAgAElEQVQDIFa1USWeEGI1cRNoNwBiVRtV4gkhVhM3gXYDIFa1UaUuYexi9cbf3U+33PkQzb35atp3+lQJZMHiZTT70h/SJV84NXNRghFgyd9n9fV1r9JD7/6JHnz3Plq7/UMaWz+eTp32OTptjzNp5Qtv08EzD6cx4yYanxTwWTVGtksG+KyGZ+hXQiGf1ehq3bXk8taNcludqi2LhF/sMqrc8raYnRWv4j2/o6dykAjqNLV3m509RdTi6fTLR9+nL19wPnEAKPcBsepNEj6rwUY6AizluJWSWA3yfR1s9CSTC2I1Ge5Bag0iVhENOAhp+3liF6tHzrqMzjrl2F1EKYvYex9+OnNRgiFWdxWrSzYukgL1wXfuE/ukLqNhdcPp1N2FSN39TDpkzGFyFD/ywFyIVfvns3aJEKvaqIwTpkWs+jZcrG6obHpHCNf3xKws/4q/Ny+S/7PAdR8/2PZ9umbQ94X/ayV1DRgnt9hRr+XiAVRVwxhqKh+Ve69+kjGvUswAsRrMqhCrEKvBRk5yuSBWk2NvWjPEqimx9KSPXaxyNGCvJb9qaTCiAZsNjjRFA+blvixQWai+vWmxDJx02h45kXrE+GP7dAxi1czOtlNDrNomurO81IvVAl2Xs7EsYlm8sogVYva6RQfT/9f4C+L9ZHUOFrM8C7tT2I6Vf3ezyO0Vuj3CFaCUD4jVYNaFWIVYDTZykssFsZoce9OaIVZNiaUnfexiNWszq17+tWw+p6ju7z6ra7atzs+kvrnhNaoUMzAsUHm573GTTkzFaIfPairMoNUI+KxqYYo9UXnrBrlHbPmOdeJ1jXhdS9Xt66mqdR11NX8o/68QEY3LOrcXbFtPeRV1DxhNXbWjhKgVr+Lv3OuY3Pvib/U+UVns/SzVCuGzmh3LltIy4OxQD9ZSiNVg3JLIFUSscjsRDTgJa/WtM3axmjWfVRarP7/lnoLLk/urWN3UslGK1IeW/onmf/iSHFknTTtdClXeLzVNB8RqmqxRuC0Qq9mxldtntaxjK1VI0SpErRCxFdvFq/ib35P/7xAilz9r21y0k3J2VolZIWSliBUCNy9u+f+Bo6mnvKZoWUggtgirraTKijJq2t4BHCknALGacgM5mgexmh1bQaxmx1bulsYuVrkBWYoGDLHqPbhZpP7X339IvPSXjxnD96PvHHoNHT/5M6k8GyBWU2kWz0ZBrGbHVkECLJV1t/XOzvIsrUPYspDl//PCdq0A0VMQRnfN0J2zsixs+8zUCjHbO4PbU1WfHagRtRRiNSKwERQLsRoB1IiKhFiNCGwExUKsRgA1piITEasx9c1KNV7C2rkEuL8FWOK9Uv/zuW/RglX/ojniZ97gh+nKw34oZ1RNDvismtCynxY+q/aZqhKz7LPqRYWjIV70te94AgsiVnXJl4slxc5lx3kRK5ci94pcsSS5rKulYJEc4XjnUuOdS4x3LjcWy495llYIX1sHfFaDkYTPao5bKYlVRAMOdi4gl30CQcQqogHbt0OQEiFWDal96fLraOOmrfTg7dfKnCxWr7nmGu1StmzZQnfccQd94xvf0M5TKOENN9xA559/Pg0ZMsS4vGeffVbmOfroo4vm3da+ja59/lr675f+W+6VOrFmIl1UfhFd8X+ukD6qpgcz4HonT55smpWWLFlCb7zxBp199tnGeYtlWL58OTEXZmr7WLt2LT344IP05S/7b/dju06/8kxsH3WbwozhqNsWpPzrrrtOnt+1taURRMj0GheEWeA8bVuItq0h2i5+hYDN/83v8f/8Pv/N6QodFWIp8aAxRGImlgaKV/6bX91/i2XHOkdrayvxuP7ud7+rkxxpegmk6RoJo9ghEOX3tZ0WohQQ8CeQ6u+/fmS42MSq8lX12ku10Gdps4XaE1bNrpb6zGpndwfds+QuunfxnfTPNX+n+uoGOmv6uXTqhH+jpc+9SbPPCya8MLOa7MjGzGp0/DGzGh3bICXzzGsuKJTbf5ZnZ8XSY+FXWyH8aXkmt/BR1tdflpcYc1Aoufy4N1gUz9CK99s6uunuP9xC519k56FkkH6782RhGTBmVnNWw8yqjREfTxlYBhwPZxu1YGbVBsVkyohNrJ425yoa1thAv7/e+0mze8YyGRzFa/XaYqdUAyw9t+pJIVLvogfeuUeCOWHKSUKofoE+M/XU4qBSlgI+qykzSIHmwGc1O7aKchlw3BR4ax4ZDGr76t7flVLAyv935H699qF1t5O35VHb9vRUD6bu6iHUUzVQLDMeRj1iBpeDRPVU1srtfXrKcnvXyj1sxbY+UR5ZEKtR9j9LZZeSWM0S9yBthVgNQi2ZPEHEKrcU0YCTsZez1tjEqt/+qqoxad1nlbfaeX7ejXlmXqK71MTq+01LpUi9R8ymrhPL6KYP34fO3usLdOZe59KQWnt+XXEOf4jVOGmHqwtiNRy/OHOXkljV4iZcIFi8Vm5fIWZkhYCVfrVi2x71N79uW6FVlF+irkGT5EcsXlnM8v608pXFLYtdfhURkLvrhgvRK0SwEMPdVYOpp0a8CnHsd0CshjJLrJkhVmPFHaoyiNVQ+GLNDLEaK26rlUGsFsHJ4nTp8tX5VB8/aPous8OlIlb9lvzybOo+w/e3OvDiLgxiNW7iweuDWA3OLu6c/U6sagLOz8YK8VrWsY3KeNa2o0nuQcuzs2Wd23KvXa25vWnFa/FlyHqVc0Cp7lqexRUzvELYUmWN2NN2NFVWVRHVjxdLlIUQ5tlcFsT1DmFMvbO7AWIQ6LUMqXQJQKzqkko+HcRq8jbQbQHEqi6p9KWLTazyDOW3LjmbZp14hCcFnS1i0ocvF2Dp8m9fSc0tnVrNa97aRI/OmxvY19Ndydw/3EonzZpN9Q3+T9T9Gub0WzRd8hu2H/BZ1RoukSWCz2pkaAk+q9GxzUrJgaIB93RTeftWIW637nwVf/f53/V5mfi/nNPwq/ysmXhroDCHFLvV9SSXMIstf3Z5rWoQy5obRBrXqzNt1aBATYDPag5bKYlVRAMOdCogUwQEgohVRAOOwBABioxNrH7vx7fRondW5KPouttazKc1QN9iyZJ1sbq5dRP9s+pVuvftu2jttg9p7+H70lliuS/Ppg4Wy8r8DohV8+GVphsxiFVz++nmgFjVJVW66QKJVUs42Pe2TEQ/5lnciraPcrO5LR9RTXmH+H89tbfskLO4ZT2dcvkyH5W9S5fDLmHOd4F9cMXsrfTJ7fXFdS5vprKK3Owuv4qgVHL2V8wCb9jSSk+//CqdMftLlmhksxiI1ezYDTOr2bEVxGp2bOVuaWxilSvm2VU+nD6g6v1NW5rJuX9pVpBmVax2dnXQXY+JBwgfvUl/bL6T6mtElN9ekaqz5Bdi1XyEQqx6M3OvDsj6MmCIVfNzo9RyJClW/Via+KzyMmZi/1wRVEqKWha3nWKpchsvYRbil187xGu7EMW8xFm8SpHcu7w5qD3Xdo+mB1tn0ZcH3JIvQolcfoMDVXXVDM9/Jpc7i5lgdXQNmpj/W/rziqBW+c96RbH6v3Ngbgl0rtzc0um0HBCrabFE8XZArBZnlJYUEKtpsYR5O2IVq9w8nmF96PGX+rTUyw/UvCvJ5ciaz+pzK5+k+5b8ke5/Z66E9umpJ8vgSVmM8qtrdfis6pJKPl3WxWryBONrAXxW42MdtiYTsRq2LumDK7YLKuvO+eRK4ds7i1uxfaUsvmJb76sQxOrzvDAW+ZM8csGscvskO2eH5f+9EZ5V+1TgK/m/mjHu/bBbrE7i9OqQQbOEb7AspzcitFc/S0msJmnHOOqGWI2Dsp06gohVrhnRgO3wD1NK7GI1TGPTmjcrYvX9LSLKrxCpvGfq2u1iye8wseRX7JlabMlvWrmbtAti1YRWsmkhVpPlb1I7xKoJrWTTxilWbfSUfW95NpeEcC3rbJFBqHLv9f7Ngjb/mXhf7Kcr08j3xGdCKOfSFvmM08i0ubzU02Wj+fplCIHbU1knIz33lNfKv8urcq9dJN4TopaFc+63N534jGTa3vc5P//PZfBn+fS5cjh6dK4OIcD5M84n3sMRngDEaniGcZUAsRoXafv1QKxaYJp2sdopnmrzVjT3LrmT/rHmZaoXgTHOFj6pPJu6z4hsR/nVNR/Eqi6p5NNBrCZvA90WQKzqkko+XdbEapLE1DJobkMZ7ZwVlv/zsmfhE6wOXiqtBC4LY2dUZ7VMOp+Wo0OLZdWyHIsRoIOyUpGjVX5eCt3dZ1m1c6l032XVagsllVdFlub/uQznEuzuATtnqoO2Na35IFbTapld2wWxmh1buVsKsRrSdmn3WX1+1VNSqOaX/E7pXfK726kUJsgOfFbNBw58Vr2ZwWfVfCzFmYOjIV70te94Vgmx6m2JrPusxjm+nHWl6RoZNwO5lZHY5oiP+gFV1NP0PrW2987ydraJ5dRr801SvsN5Edy7pFqKYJffsAyyJfyLd6YNtw+wKZclnXvRG50H0Nm1ObejgkdvYK5iyfhz3lOY9xjWOdRexcXSqn2Mi6WT9Qt/aPaL5odA7Z3d1CF+/Q6nkC9UtlvkF0xbwg8AdPgHSRNErCIacBDS9vNArIZkmlaxurzpPbqHZ1P7LPn9ggyiNLg2d4GHWN3V+GtWr6RX579EJ58+O+TI2DV7mm7EwtjeNhiIVdtE7ZYHsWrOE2LVnBnnSNM1MlgP7OSKy2dVCtuWjflGV7Tm9vyVh8PPmP9VgbVUYhl0S8wky88ckaX5fxVwi/9+u3UiLdg+ns6uudsOHJRiRiCDDwB0OhjkAcDQ+mra3tpJ7R3eDxa8VgBArOpYI/o0EKshGadNrPot+WW/1BnD9+vT2zCCBTOr5gMnTTdiYWxv3vPCOSBWbRO1Wx7EqjlPiFVzZpwjTdfIYD2wkysusWqntYVL8d9ntUcI4i6xzFrMHrOfMP8t9hqWf4vfst73ckusxf/i89xn3bm/OZ/Mo953/C3TinTqM+rN01tWPp8spzddb/m7lNebxq8uMbFKnZ2d1CV++7RblNun3fm+ivcL1qX6mmOR72uekbO/jrJ6eYhMcZi139Txg23fp2uuuabf9DetHYVYtWCZtPisei35ZZF64tRTLPQy20XAZzU79oPPanZshWXA2bEVfFazY6tSEqvZoR6spWn3WVWRuXV6x8vLy8Vsu85R0ay3nNw9I1+obOdsfaF07pn8gmX2bqvFafg+sKenx1fOS/9zsaLAeTQfcBXVH/cjHSRIEyEBiFULcJMWq2rJ731L7qI121b3Rvn9ggyi1FCzM2y+ha5mtgiI1eyYDmI1O7aCWM2OrSBWs2MriNXs2CrtYjU7JKNvaRCfVW4Vtq6J3jbFaoBYLUZI4/MkxeqD795H//fFK6VI5ePfZ8yh7x36f6mxLj0bnGsgjDwJxGrkiK1VALFqDWXkBUGsRo7YWgUQq9ZQRl4QxGrkiK1VALFqDWXkBUGsRo44sgogVkOiTcpndZMIivDNJy+mfZbvSbeLn0kjp9KPjrqeDho9U7tHYfwW4bOqjTmfME3+WGFsb97zwjngs2qbqN3y4LNqzhM+q+bMOEearpHBemAnVymJVX+fVTuski4FYjVpC+jXH0SsIsCSPt8oU0KshqSbhFh9YvlfpFCVgrXscmo8ZCyd97GLqVJEfTM5wggWiFUT0rm0aboRC2N7855DrNpmFmd5EKvmtCFWzZml7RoZrAd2ckGs2uEYRykQq3FQtlMHxKodjkmUArEaknqcYnXF1mV0+4Lb6A8LfkOtnS30mamn0hFrD6V/+9wXqb7B3Dc1jGCBWDUfOBCr3swws2o+luLMAbFqThti1ZwZxOpOZhCrwcZPErkgVpOgHqxOiNVg3NKQC2LVghXi8Fn96/uP0B1CqD678gkaVjucztvvYjpvn4topNgYGkdxAvBZLc4oLSngs5oWSxRvB3xWizNKSwr4rKbFEsXbUUpitXhvs50CYjU79gsiVrl3CLCUvI0hVi3YIEqx2tS2RYjUW+mOt35Da7d9SJ8YfwydL0TqSdNOt9Dy/lMExGp2bA2xmh1bQaxmx1YQq9mxFcRqdmwFsZodW0GsZsdW7pZCrFqwXVRi9ZU1L9Idb95GHPG3sryK5vTOpu42ZA8Lre5fRUCsZsfeEKvZsRXEanZsBbGaHVtBrGbHVhCr2bEVxGp2bAWxatlWUfms/kHMpN7x1m205KOFtO+IA8SS34vltjTuw+3vZ9I9+KzuSmvN6pX06vyX6OTTZ5ug1EoLn1VvTPBZ1Ro+iSWCz6o5evismjPjHGm6RgbrgZ1cpSRWEQ3YzphAKeEJBBGriAYcnruNEjCzGpKibbH69qbFctkvB1HqET9n7nUunb/vRXTgqEM8WwqxqmdA3ZlViFU9njZTQazapGm/LIhVc6YQq+bMIFZ3MoNYDTZ+ksiFmdUkqAerE2I1GLc05IJYDWkFm2L1oXf/RDyj+vLq52nS4KlSpHIQpbrKAb6thFjVMyDEal9OYWbV9Yjrp4JY1WeVREqIVXPqEKvmzCBWIVaDjZpkc0GsJsvfpHaIVRNa6UoLsWrBHmF9VtfvWNu7Jc1ttLl1E316yslCqF5MR0883kLrUAQT0BWroJU8AfisJm8D3RbAZ1WXVPLp4LOavA10W1BKM6u6fc5qOojV7FguiFjl3iEacPI2hli1YIMwYvW5VU/KJb9/WfYQDa1tpPOESOVov6MGjrHQMhShCECsZmcsQKxmx1YQq9mxFcRqdmwFsZodW0GsZsdWEKvZsZW7pRCrFmwXRKy2d7XJSL+8Jc37TUvpsHFHSZF6yu7/ZqFFKMJNAGI1O2MCYjU7toJYzY6tIFazYyuI1ezYCmI1O7aCWM2OrSBWLdsqiM/qA//vD/TmpMV096LbxZY0ldIv9Tzhn7r70L2MWwefVT1kumIVAZb0eNpMBZ9VmzTtlwWfVXOm8Fk1Z8Y5EA04x62UxCqiAQc7F5DLPoEgYhXRgO3bIUiJmFkNQs2Rx1SsPrLwfnrnmdfo5z0/p8a6YfTTo39JJ007PXArIFb10EGs9uWEAEt64yZIqjt+cwOdc94lVF1TEyR76vJArJqbBGLVnBnE6k5mEKvBxk8SuTCzmgT1YHVCrAbjloZcEKshraArVju7O+knL19Nc1+7neaInzcmLaSfHftrGjNoXKgWQKzq4YNYhVjVGynhU0GshmeY9RIgVoNZEDOrOW4Qq/9/e3ccq1dZ3wH8yWIi/kGtReMmxDUEddXUKSyQkDL8a+EvaFDXzszRtWZQtCHBLdU1oQlZHf1DDGkmzEFTYpS6OWm7LSHEZJA2qGwysCGtgTRdJrPOgbUss8n4Yz1vc27ufbn3nud5znnec973fprwB73P83ue9/M757332/Oec/OOnz5mCat9qOetKazmuQ1hlrDaQRea7ln96f+8Eu588rbw7H8+M/rY759dd0/43DV/2sHKSsQKxIbV2HrGlRNwz2o5264ru2e1a9Fy9dyzWs6268qzFFa7thlaPWF1aB1Zej85YbWq5mnA/fdYWO2gB8uF1W+f/EZ45IWvhh/9/N/Cte+5Pmxdv91DlDowTy0hrKaK9TdeWO3PPnVlYTVVrL/xwmp/9qkrC6upYv2NF1b7s09dWVhNFRvOeGG1g14sFlbPnv9FeORHX73w31+FX54/Gz75W38Ytv32nWH9uz7SwYpKpAoIq6li/Y0XVvuzT11ZWE0V62+8sNqfferKwmqqWH/jhdX+7FNXFlZTxYYzXlht2YvF7ll9/mc/DPsvBNW///Fj4bK3vTNs/fD2sO3Dnw2XvnVVeP3cL8M/HToYNv/R7S1XvjjdPatxjLFh1dOA4zy7HOVpwF1qdl/LA5bSTd2zmm5WzXDP6kW3WQqrngacdy6Y1b1ATlj1NODu+5BTUViNULtly67w8ulXRiOvWnt5OHxgz9ys8bB65KVvj66o/utPvx+u/vVrR1dTN77v9+fGC6sXKf7x8YPhmmuvD79x+XsjOrBwSM43P2F1oaGnAScfdtETPGApmmpmBwqrea0VVoXVvCOnv1murPZnn7qysJoqNpzxwmpDL7bevTe8+tq5uYBaBdfL1qwK++/fOZpZh9Uz586O7k3df+G/n//qv8LHP/AHoyuqH3n37yxYQVgVVm/dvKX3dwBhtVwLhNVyttNSWVjN65SwKqzmHTn9zRJW+7NPXVlYTRUbznhhtaEXN2zcET5/x6aw8aYNo5GHnjgWvvzQt8LRQ/vmZn7v358L+/5lX3jsxQPh7W9dfeFq6mcvfOz3zrD6kncMp9MrfCexV1ZXONMgXr57VgfRhqhNuGc1imkQg9yzOog2RG1ilj4GHPWCp3iQsDo9zcsJq9Wr8zTg/nssrC7Tg+MnToXN2+8NBx+8J6xfd+Vo5PjfPX7y8fCVZx4IR//j6QtXUa8Jf3zhauonPvCp/jtrBwsEhNXpOSCE1enplbA6Pb0SVqenV8Lq9PRKWJ2eXgmr09Or8Z0Kqy3D6kf/+qPh+TPPhw3v3RAe+/hj4YpVV0zv0WDnBAgQIECAAAECBAgQGIiAsNoyrFb3rL7xu2+E3TfuDm/5tbc0tvXs2bPh0UcfDXfddVfj2JgBDzzwQLjtttvC6tWrY4YvGPP000+P/v/GG29Mntv2dVQG1bpr165NXvvkyZPhhRdeCJs2bUqe2zTh9OnToXKpTLv+c+bMmXD48OFw++3dPAm6zf7a9L7NuovNbXMMd72XLurt3bt3dH5fcsklXZTrvUb1Hrd79+7e9zFNGzh//nyojuudOy8+28CfOIEhvUfG7dioJoGS36+b1vZ1Am0FfP9rK9jNfGG1wXGxe1Z33fdwePGpA6OZi/3qmuVKesDSRR1PA+7mBM6t4gFLuXLN8zxgqdlo1kd4wFJehz1g6aLbLH0MOOfp/XlHTz+zfAy4H/ecVXM+BuxX1+RIdz9HWG0wjX0a8Ou/eiOqO8KqsOppwAtPFb9nNeqto7dBfs9qOr2wmm5WzRBWhdW8I6e/WcJqf/apKwurqWLDGS+sRvRiud+zWk1//X//L8SG1YjlDCkg4AFLBVALlfSApUKwBcp6wFIB1EIlPWCpEGyBsrN0ZbUAz6BKCquDaseym8kJq1VBTwPuv8fCagc9EFY7QCxcQlgtDNxheWG1Q8zCpYTVwsAdlhdWO8QsXEpYLQzcYXlhtUPMwqWE1cLABcsLqx3gCqsdIBYuIawWBu6wvLDaIWbhUsJqYeAOywurHWIWLiWsFgbusLyw2iFm4VLCamHgguWF1Za4HrCU92RbD1hqeeC1nO4BSy0Bl5nuAUvlbKelsntW8zrlntWLbrMUVj1gKe9cMKt7gZyw6gFL3fchp6KwmqM2b46wKqzGHkJD+kFMWI3tWvo4YTXdbNZmCKt5HR3Se2TeK+hmlrDajeMkqriyOgnlbtYQVrtx7KOKsNpSXVgVVmMPoSH9ICasxnYtfZywmm42azOE1byODuk9Mu8VdDNLWO3GcRJVhNVJKHezhrDajWMfVYTVPtStSYAAAQIECBAgQIAAAQLLCgirDhACBAgQIECAAAECBAgQGJyAsDq4ltgQAQIECBAgQIAAAQIECAirjgECBAgQIECAAAECBAgQGJyAsNqiJbds2RVePv3KqMJVay8Phw/saVHN1LYCKf049MSxsOu+h9+05ItPHWi7DfNbChw/cSps3n5vOPjgPWH9uitbVjO9C4HYnjivutDutsbWu/eGHzx3YkFR73PdGqdWS+2J8ypVuPz4L37pa+HIk884r8pTR6+Q2hPnVTRt7wOF1cwWVN9sXn3t3FxArYLSZWtWhf3378ysaFobgdR+VG9SX37oW+HooX1tljW3Y4EbNu4Ir519fVRVWO0YN7NcSk+cV5nIBadV/Zv/Plf9QHfs2ePe+wqaN5VO7Ynzqkl08l+vfub7i53b5v5Bdd8j3wl/+w//7LyafCvmVkztifOqx2YlLi2sJoLVw6tvNp+/Y1PYeNOG0V856DMhO5qW2g/96gi+QJnYq3gFllZyCYHYnjivhn8IxfZy+K9kdnbY1BPn1fB73dTD4b+C2dthU0+cV9PTc2E1o1eLnQBNJ0XGMqZECuT0Y7GPf/hoXCR44WHOpcLAGeVje+K8ysCd8BRXgCYMHrFcU0+cVxGIPQ+pPt310qmfuLLacx/mL9/UE+fVgJrVsBVhNaNXOeEoYxlTIgW66Mf4x4gjlzasgEBsMCqwtJJLCOT2xHk1rEOq7uOeL3xm7lNBw9rhyttNTk+cV8M5TubfJuEfvIfRl9yeOK+G0b/FdiGsZvSmi3CUsawpCT9Ip/5wXY/3zab/wyy1d/3vePZ3kNsT59Vwjo26F3d8+uawY9utw9nYCt5Jbk+cV8M7aKqr4w99/UjwM8RwepPaE+fVcHo3vhNhNbM3i90jWT1d1htVJmjLaW37UX8cRP9aNqKD6bnBqIOllUj4B6EYLOdVjFL5MXUfPLSsvHXsCm164ryKVZ7suA99bIsHA06WvHG1lJ44rxo5exsgrGbSpz59NnMZ0yIFmvpRPSWu+lP/eqHxpzF6mnMk9ASGCasTQE5cYqmeOK8SIXsY7iEiPaA3LNnUE+fV8Ho2vqPUJzoP/xVN/w6beuK8mt4eC6stepfyez1bLGNqpMBy/Rh/k5o/tip/3dXr/NqhSOeSw+bfa1Kts2b1pR5YURI8ovZyPXFeRQD2OKT+R4bFtuC+1X4aE9MT51U/vUlZdfxniGquT2alCHY/tqknzqvuzSdVUVidlLR1CBAgQIAAAQIECBAgQCBaQFiNpjKQAAECBAgQIECAAAECBCYlIKxOSto6BAgQIECAAAECBAgQIBAtIKxGUxlIgAABAgQIECBAgAABApMSEFYnJW0dAgQIECBAgAABAgQIEIgWEFajqQwkQIAAAQIECBAgQIAAgUkJCHmeX0QAAAd9SURBVKuTkrYOAQIECBAgQIAAAQIECEQLCKvRVAYSIECAAAECBAgQIECAwKQEhNVJSVuHAAECBAgQIECAAAECBKIFhNVoKgMJECBAgAABAgQIECBAYFICwuqkpK1DgAABAgQIECBAgAABAtECwmo0lYEECBAgQIAAAQIECBAgMCkBYXVS0tYhQIAAAQIECBAgQIAAgWgBYTWaykACBAgQIECAAAECBAgQmJSAsDopaesQIECAAAECBAgQIECAQLSAsBpNZSABAgQIECBAgAABAgQITEpAWJ2UtHUIECBAgAABAgQIECBAIFpAWI2mMpAAAQIECBAgQIAAAQIEJiUgrE5K2joECBAgQIAAAQIECBAgEC0grEZTGUiAAAECXQnse+Q74aGvH3lTuTs+fXPYse3WcMPGHaOvHT20701jqq+tWb0qHD6wZ/S1plof+tiWZbe9ZvWlo3W23r03/OC5E4uO3fOFz4SNN20It2zZFV4+/Uqo/78efOiJY2HXfQ+Hq9ZePrev8UIx+9hw7fpw5Mln5qbe/HvXh7/88z9JWjfmdXTVR3UIECBAgEBJAWG1pK7aBAgQIPAmgTpMHXzwnrB+3ZVzX69C53eP/nAu7FXh7rqr14X99++cG/PFL30tHHv2+FyIja01HirHw2b19arWq6+dWzJsVmPqsDq+r/rvlwur8yHqcLvYPhb7Wsq6Ma/DYUmAAAECBKZBQFidhi7ZIwECBGZIoAqh9RXD5V7WeGg7fuJU2Lz93gVXNWNrdRlWL1uzanQFtg7b9b6qANsUdmP2sVRYjV1XWJ2hk8VLIUCAwAoXEFZX+AHg5RMgQGDSAtXHeN935RULrpgutYcqeL106iejK6nV1cUqsM2/0ppSq1pjuSuaMSGv2sMH3/+b4Wf//Yvw7ne+Y/QR3epqb/Wn+ruSYTV23ZjXMemeW48AAQIECOQICKs5auYQIECAQLZAHRjrAvU9o0sVnH+v54tPHVgwLLVWU1iNuWe1Co3XXf3B0T2q1X6q/VVXWb/yN39XPKzGrOue1exD00QCBAgQGJiAsDqwhtgOAQIEVpJA/RHa+jUv9vHgOmDWD19ayielVpt7VquwWj/0qNpLfbU35Ypmzj2rseum7GMlHWteKwECBAhMn4CwOn09s2MCBAjMpED1cdrqSbjjV08Xu1e1CWCpWk1XVps+xlt/DLgKq/VTiOvgmxIS24TVpnVT9tHk6OsECBAgQKBPAWG1T31rEyBAYIUJVMHzm49/d3RlcvxPHcLGnxK8VFjNqdVlWK32X90zW/96nZSQ2CasNq2bso8Vdvh5uQQIECAwZQLC6pQ1zHYJECAwzQLzP6o7/wrq/Cfqzn+AUvValwur1dOBqz+xtboOq/N7kRIS24bV5dZN2cc0H0v2ToAAAQKzLyCszn6PvUICBAgMTmD+Q5PqzS11T2rTx4BTajWF1dgHLC12ZTglJC61j/rjy7VJfQ/v/I8fjzdzfF0PWBrc4W5DBAgQIJApIKxmwplGgAABAgQIECBAgAABAuUEhNVytioTIECAAAECBAgQIECAQKaAsJoJZxoBAgQIECBAgAABAgQIlBMQVsvZqkyAAAECBAgQIECAAAECmQLCaiacaQQIECBAgAABAgQIECBQTkBYLWerMgECBAgQIECAAAECBAhkCgirmXCmESBAgAABAgQIECBAgEA5AWG1nK3KBAgQIECAAAECBAgQIJApIKxmwplGgAABAgQIECBAgAABAuUEhNVytioTIECAAAECBAgQIECAQKaAsJoJZxoBAgQIECBAgAABAgQIlBMQVsvZqkyAAAECBAgQIECAAAECmQLCaiacaQQIECBAgAABAgQIECBQTkBYLWerMgECBAgQIECAAAECBAhkCgirmXCmESBAgAABAgQIECBAgEA5AWG1nK3KBAgQIECAAAECBAgQIJApIKxmwplGgAABAgQIECBAgAABAuUEhNVytioTIECAAAECBAgQIECAQKaAsJoJZxoBAgQIECBAgAABAgQIlBMQVsvZqkyAAAECBAgQIECAAAECmQLCaiacaQQIECBAgAABAgQIECBQTkBYLWerMgECBAgQIECAAAECBAhkCgirmXCmESBAgAABAgQIECBAgEA5AWG1nK3KBAgQIECAAAECBAgQIJApIKxmwplGgAABAgQIECBAgAABAuUEhNVytioTIECAAAECBAgQIECAQKaAsJoJZxoBAgQIECBAgAABAgQIlBMQVsvZqkyAAAECBAgQIECAAAECmQLCaiacaQQIECBAgAABAgQIECBQTkBYLWerMgECBAgQIECAAAECBAhkCgirmXCmESBAgAABAgQIECBAgEA5AWG1nK3KBAgQIECAAAECBAgQIJApIKxmwplGgAABAgQIECBAgAABAuUEhNVytioTIECAAAECBAgQIECAQKaAsJoJZxoBAgQIECBAgAABAgQIlBMQVsvZqkyAAAECBAgQIECAAAECmQLCaiacaQQIECBAgAABAgQIECBQTkBYLWerMgECBAgQIECAAAECBAhkCvw/ZVvJd7BP4YoAAAAASUVORK5CYII=",
"text/html": [
"
\n", " | SYSTEM TIME | \n", "A | \n", "B | \n", "E | \n", "caption | \n", "
---|---|---|---|---|---|
0 | \n", "0.000000 | \n", "20.000000 | \n", "0.000000 | \n", "30.0 | \n", "Initialized state | \n", "
1 | \n", "0.000250 | \n", "18.500000 | \n", "1.500000 | \n", "30.0 | \n", "\n", " |
2 | \n", "0.000375 | \n", "17.817500 | \n", "2.182500 | \n", "30.0 | \n", "\n", " |
3 | \n", "0.000500 | \n", "17.165712 | \n", "2.834288 | \n", "30.0 | \n", "\n", " |
4 | \n", "0.000625 | \n", "16.543255 | \n", "3.456745 | \n", "30.0 | \n", "\n", " |
5 | \n", "0.000750 | \n", "15.948809 | \n", "4.051191 | \n", "30.0 | \n", "\n", " |
6 | \n", "0.000875 | \n", "15.381112 | \n", "4.618888 | \n", "30.0 | \n", "\n", " |
7 | \n", "0.001000 | \n", "14.838962 | \n", "5.161038 | \n", "30.0 | \n", "\n", " |
8 | \n", "0.001125 | \n", "14.321209 | \n", "5.678791 | \n", "30.0 | \n", "\n", " |
9 | \n", "0.001250 | \n", "13.826755 | \n", "6.173245 | \n", "30.0 | \n", "\n", " |
10 | \n", "0.001375 | \n", "13.354551 | \n", "6.645449 | \n", "30.0 | \n", "\n", " |
11 | \n", "0.001525 | \n", "12.813405 | \n", "7.186595 | \n", "30.0 | \n", "\n", " |
12 | \n", "0.001675 | \n", "12.301481 | \n", "7.698519 | \n", "30.0 | \n", "\n", " |
13 | \n", "0.001855 | \n", "11.720345 | \n", "8.279655 | \n", "30.0 | \n", "\n", " |
14 | \n", "0.002071 | \n", "11.068171 | \n", "8.931829 | \n", "30.0 | \n", "\n", " |
15 | \n", "0.002330 | \n", "10.346417 | \n", "9.653583 | \n", "30.0 | \n", "\n", " |
16 | \n", "0.002589 | \n", "9.692012 | \n", "10.307988 | \n", "30.0 | \n", "\n", " |
17 | \n", "0.002900 | \n", "8.980003 | \n", "11.019997 | \n", "30.0 | \n", "\n", " |
18 | \n", "0.003274 | \n", "8.221264 | \n", "11.778736 | \n", "30.0 | \n", "\n", " |
19 | \n", "0.003647 | \n", "7.564476 | \n", "12.435524 | \n", "30.0 | \n", "\n", " |
20 | \n", "0.004095 | \n", "6.882233 | \n", "13.117767 | \n", "30.0 | \n", "\n", " |
21 | \n", "0.004543 | \n", "6.309997 | \n", "13.690003 | \n", "30.0 | \n", "\n", " |
22 | \n", "0.004991 | \n", "5.830030 | \n", "14.169970 | \n", "30.0 | \n", "\n", " |
23 | \n", "0.005528 | \n", "5.346939 | \n", "14.653061 | \n", "30.0 | \n", "\n", " |
24 | \n", "0.006066 | \n", "4.957322 | \n", "15.042678 | \n", "30.0 | \n", "\n", " |
25 | \n", "0.006711 | \n", "4.580248 | \n", "15.419752 | \n", "30.0 | \n", "\n", " |
26 | \n", "0.007485 | \n", "4.232821 | \n", "15.767179 | \n", "30.0 | \n", "\n", " |
27 | \n", "0.008413 | \n", "3.932073 | \n", "16.067927 | \n", "30.0 | \n", "\n", " |
28 | \n", "0.009528 | \n", "3.691843 | \n", "16.308157 | \n", "30.0 | \n", "\n", " |
29 | \n", "0.010865 | \n", "3.519230 | \n", "16.480770 | \n", "30.0 | \n", "\n", " |
30 | \n", "0.012470 | \n", "3.411824 | \n", "16.588176 | \n", "30.0 | \n", "\n", " |
31 | \n", "0.014396 | \n", "3.357403 | \n", "16.642597 | \n", "30.0 | \n", "\n", " |
32 | \n", "0.016707 | \n", "3.337375 | \n", "16.662625 | \n", "30.0 | \n", "\n", " |
33 | \n", "0.019480 | \n", "3.333337 | \n", "16.666663 | \n", "30.0 | \n", "\n", " |
34 | \n", "0.022808 | \n", "3.333329 | \n", "16.666671 | \n", "30.0 | \n", "\n", " |
35 | \n", "0.026802 | \n", "3.333331 | \n", "16.666669 | \n", "30.0 | \n", "\n", " |
36 | \n", "0.031594 | \n", "3.333330 | \n", "16.666670 | \n", "30.0 | \n", "\n", " |
37 | \n", "0.037345 | \n", "3.333331 | \n", "16.666669 | \n", "30.0 | \n", "\n", " |
38 | \n", "0.044245 | \n", "3.333330 | \n", "16.666670 | \n", "30.0 | \n", "\n", " |
39 | \n", "0.052526 | \n", "3.333332 | \n", "16.666668 | \n", "30.0 | \n", "\n", " |
40 | \n", "0.062463 | \n", "3.333327 | \n", "16.666673 | \n", "30.0 | \n", "\n", " |
41 | \n", "0.074388 | \n", "3.333341 | \n", "16.666659 | \n", "30.0 | \n", "\n", " |
42 | \n", "0.088697 | \n", "3.333285 | \n", "16.666715 | \n", "30.0 | \n", "\n", " |
43 | \n", "0.105869 | \n", "3.333568 | \n", "16.666432 | \n", "30.0 | \n", "\n", " |