{ "cells": [ { "cell_type": "markdown", "id": "c04b84e5-d999-4445-a33d-64e73b19cb54", "metadata": {}, "source": [ "# Lifetime distributions" ] }, { "cell_type": "markdown", "id": "2a3ac1fd-2f2b-4b76-a087-72306e713f1e", "metadata": {}, "source": [ "A built-in dataset of circuit breaker lifetimes." ] }, { "cell_type": "code", "execution_count": 1, "id": "377eb9ae-f2a6-49ad-9f68-83641abbf85b", "metadata": { "ExecuteTime": { "end_time": "2025-09-02T12:31:37.313907Z", "start_time": "2025-09-02T12:31:36.977768Z" } }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "[34. 28. 12. ... 42. 42. 37.]\n", "[ True True True ... False False False]\n", "[33. 27. 11. ... 31. 31. 26.]\n" ] } ], "source": [ "from relife.datasets import load_circuit_breaker\n", "\n", "dataset = load_circuit_breaker()\n", "print(dataset[\"time\"])\n", "print(dataset[\"event\"])\n", "print(dataset[\"entry\"])" ] }, { "cell_type": "code", "execution_count": 2, "id": "8639e654-a7fb-40e4-b375-a264c94f0125", "metadata": { "ExecuteTime": { "end_time": "2025-09-02T12:31:44.881257Z", "start_time": "2025-09-02T12:31:44.877015Z" } }, "outputs": [], "source": [ "from relife.lifetime_models import Weibull, Gompertz" ] }, { "cell_type": "code", "execution_count": 3, "id": "6f224d6f-afaf-4fd6-96fa-70ef70ca0628", "metadata": { "ExecuteTime": { "end_time": "2025-09-02T12:31:46.110164Z", "start_time": "2025-09-02T12:31:46.107742Z" } }, "outputs": [], "source": [ "weibull = Weibull()\n", "gompertz = Gompertz()" ] }, { "cell_type": "markdown", "id": "64105fe5-1c37-4405-83e6-845771977b73", "metadata": {}, "source": [ "From now, the models parameters are unknown, thus set to `np.nan`." ] }, { "cell_type": "code", "execution_count": 4, "id": "a2d94106-c580-4e4c-a91f-2987c8b03021", "metadata": { "ExecuteTime": { "end_time": "2025-09-02T12:31:47.265727Z", "start_time": "2025-09-02T12:31:47.262003Z" } }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "('shape', 'rate')\n", "[nan nan]\n" ] } ], "source": [ "print(weibull.get_params_names())\n", "print(weibull.get_params())" ] }, { "cell_type": "code", "execution_count": 5, "id": "9f2b1fbd-a62f-44a5-b562-d80d6d43c03e", "metadata": { "ExecuteTime": { "end_time": "2025-09-02T12:31:47.945168Z", "start_time": "2025-09-02T12:31:47.927412Z" } }, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 5, "metadata": {}, "output_type": "execute_result" } ], "source": [ "weibull.fit(dataset[\"time\"], event=dataset[\"event\"], entry=dataset[\"entry\"])\n", "gompertz.fit(dataset[\"time\"], event=dataset[\"event\"], entry=dataset[\"entry\"])" ] }, { "cell_type": "code", "execution_count": 6, "id": "39aae3b5-efd1-4f2f-b75a-3174f20196d6", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "('shape', 'rate') [3.7267452 0.01232326]\n", "('shape', 'rate') [0.00390781 0.07579546]\n" ] } ], "source": [ "print(weibull.get_params_names(), weibull.get_params())\n", "print(gompertz.get_params_names(), gompertz.get_params())" ] }, { "cell_type": "markdown", "id": "a6073a91-f667-4786-8ad1-e5bacf2ad9a6", "metadata": {}, "source": [ "After fitting, parameters are setted inplace." ] }, { "cell_type": "code", "execution_count": 7, "id": "8001ef26-9a8a-4140-b0d4-29c787adc125", "metadata": { "ExecuteTime": { "end_time": "2025-09-02T12:31:48.302341Z", "start_time": "2025-09-02T12:31:48.288768Z" } }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "fitted params : [3.72675, 0.0123233]\n", "AIC : 2493.72\n", "AICc : 2493.72\n", "BIC : 2506.41\n", "fitted params : [0.00390781, 0.0757955]\n", "AIC : 2485.57\n", "AICc : 2485.57\n", "BIC : 2498.25\n" ] } ], "source": [ "print(weibull.fitting_results)\n", "print(gompertz.fitting_results)" ] }, { "cell_type": "markdown", "id": "f514593e-f892-4763-b100-69d8d7df0fb1", "metadata": {}, "source": [ "To plot the survival function, do the following." ] }, { "cell_type": "code", "execution_count": 8, "id": "197cb3d8-88dc-4169-905f-6d818c669b94", "metadata": { "ExecuteTime": { "end_time": "2025-09-02T12:31:49.260292Z", "start_time": "2025-09-02T12:31:49.178032Z" } }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjAAAAGdCAYAAAAMm0nCAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjgsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvwVt1zgAAAAlwSFlzAAAPYQAAD2EBqD+naQAAhaRJREFUeJzs3Xd8XXX9+PHXOXff3Ju907TpHnS3tJTZaqECIgoKKAqCCywIrQhURRR+LEFFAeErioKiIArIxtrBLJ2UUTrTpEnT7HVz783d5/fHSW4bmnFv1s14Px+P+8jNmZ97Csk7n/F+K5qmaQghhBBCDCNqohsghBBCCBEvCWCEEEIIMexIACOEEEKIYUcCGCGEEEIMOxLACCGEEGLYkQBGCCGEEMOOBDBCCCGEGHYkgBFCCCHEsGNMdANiEYlEOHLkCE6nE0VREt0cIYQQQsRA0zRaWlrIz89HVfu3z2RYBDBHjhyhsLAw0c0QQgghRC+Ul5czZsyYfr3msAhgnE4noD+A5OTkBLdGCCGEELFwuVwUFhZGf4/3p2ERwLQPGyUnJ0sAI4QQQgwzAzH9QybxCiGEEGLYkQBGCCGEEMOOBDBCCCGEGHaGxRwYIYQQI4OmaYRCIcLhcKKbIvqBwWDAaDQmJMWJBDBCCCEGRSAQoLKyEq/Xm+imiH5kt9vJy8vDbDYP6n0lgBFCCDHgIpEIJSUlGAwG8vPzMZvNkph0mNM0jUAgQG1tLSUlJUyePLnfk9V1RwIYIYQQAy4QCBCJRCgsLMRutye6OaKf2Gw2TCYThw4dIhAIYLVaB+3eMolXCCHEoBnMv9DF4EjUv6n8lySEEEKIYUcCGCGEEEIMO3EHMG+++SbnnXce+fn5KIrC888/3+M5GzduZP78+VgsFiZNmsRf/vKXXjRVCCGEGP42btyIoig0NTUB8Je//IXU1NQ+X3fp0qVcf/310e+Lioq4//77+3zdoSruAMbj8TBnzhweeuihmI4vKSnh3HPPZdmyZezcuZPrr7+eb3/727z++utxN1YIIYQYTI888ghOp5NQKBTd5na7MZlMLF26tMOx7YFJcXFxt9c8+eSTqaysJCUlZSCaPGrEvQrp7LPP5uyzz475+EceeYTx48fzq1/9CoDp06fz9ttv85vf/IYVK1bEe3shhBBi0Cxbtgy32822bds46aSTAHjrrbfIzc1l8+bN+Hy+6MqbDRs2MHbsWCZOnNjtNc1mM7m5uQPe9pFuwJdRb9q0ieXLl3fYtmLFig7dXJ/m9/vx+/3R710uFwDv3/cFLDY7AcVMADNBxUxAOfo1gJmAasGvWAgoFvyKDb9iwa9Yj9lmJaBYCapmFNWAooDalotAVRRURa+aqSoACgYVjKqKQVUwqor+1dD29VPbLUYDFpOKxahiMRowG9vfq1hMhuh7h8VIksWIw2okyWzEoEouBCHE6KNpGq3Bwc/IazMZYs5BM3XqVPLy8ti4cWM0gNm4cSPnn38+69ev57333ov2xGzcuJFly5YRiUS45557+MMf/kBVVRVTpkzhlltu4ctf/nKH4xobGzsMHT3//PP86Ec/ory8nDPOOIM//vGPFBYWAvDNb36TpqamDtM2rr/+enbu3MnGjRv7/EyGowEPYKqqqsjJyemwLScnB5fLRWtrKzab7bhz7rrrLn7xi18ct31ecAfJ/fTLPqIpeLHgwYpbs+HBikez4cba4b1bs+EiCZdmpxk7Li2J5rbvXSThxQL0rU1Wk0qS2YjdYtC/mg04rUZSbWbSksxkOMxkJln0rw4zWU4LGQ4LdrOk8RFCDF+twTAzfjb40wk+uW1FXD8/ly1bxoYNG7j55psBvaflxhtvJBwOs2HDBpYuXUprayubN2/myiuv5K677uJvf/sbjzzyCJMnT+bNN9/k61//OllZWZxxxhmd3sPr9XLHHXfwxBNPYDab+f73v88ll1zCO++80y+feSQakr8B16xZw+rVq6Pfu1wuCgsL2TXvZyRbVNSwH0PYhxoJ6F/D/qPbwn4M4VYModaOX8M+/X1E79lRFQ0HPhz4yFGaet3WCCo+QxKtqgOPIZkWNRm34qRZTaaJZBpJpkFzUh9xUKs5qQ45qQza8QQihCIaAL5gBF8wQL0nvnubjSopNhOZDgs5yRZyk63kpdrIT7EyJs3GmDQ7yVYTRoPea2Q2qJL5Uggh4rRs2TKuv/56QqEQra2tvP/++5xxxhkEg0EeeeQRQB9t8Pv9LF26lBkzZvC///2PJUuWADBhwgTefvtt/u///q/LACYYDPLggw+yePFiAB5//HGmT5/Oli1bWLRo0eB80GFmwAOY3NxcqqurO2yrrq4mOTm5094XAIvFgsViOW77CZ/7LsnJyX1rUCQCQS8EPBBw6y9/+9eWo9/7W8DfDK3N4GsCX7P+8rva9rkgEkYlgj3cgj3cQkawMrY2qEZITSViSyNkScNvycBnTsdtSsdtyqTJkE4t6ZRH0qltVanzhmloDdPQGqKpNUSTN0QoohEIRaht8VPb4md3F7d2WIx6j02SmfxUG0WZdiZlOZiYlUSyzXzMMJc+5CXDWUKIwWIzGfjktsGfC2kzGeI6funSpXg8HrZu3UpjYyNTpkyJ9qZcccUV+Hw+Nm7cyIQJE3C73Xi9Xs4888wO1wgEAsybN6/LexiNRk488cTo99OmTSM1NZXdu3dLANOFAQ9glixZwiuvvNJh29q1a6OR6aBTVbA49Bc5PR7eJU2DYGtbYNME7hrw1ICnDty14K0Dbz20NkBrk/7yN+uBUyQE3jpUbx1mwAw4gaxP30NRwZ4BSVn6KycbkrLQ7Fl4bblUKdlUheyUuaHcpVHREuGIO0x1S4g6bxBvIILbH8LtD1FS54FDjUcvrUCO82hPjf7VxrgMO06rCbvZgL1tWMtuMmA0SMogIUT/UhRlWAyFT5o0iTFjxrBhwwYaGxujvSj5+fkUFhby7rvvsmHDBj7zmc/gdrsBePnllykoKOhwnc7+MI+VqqpomtZhWzAY7PX1RoK4/8txu90cOHAg+n1JSQk7d+4kPT2dsWPHsmbNGioqKnjiiScAuOqqq3jwwQe58cYbufLKK1m/fj3//Oc/efnll/vvUySCooDZrr+S8yB7es/nhEMQaIHmCnBV6F9bKo8GP956aG3UX94G0MLgqdVfx94aSAImAhPNDk5x5kFKASTnw5gCwsn5BOy51AWtHGqBcrdCmVuhpEmjrCVCeXOQFn+EKpePKpePbccENqoChWl2Juc4mZLjYEqOk2ynBYvJ0BbUGLC1BTdJ5tgnwgkhxHC2bNkyNm7cSGNjIz/60Y+i208//XReffVVtmzZwtVXX82MGTOwWCyUlZV1OVzUmVAoxLZt26K9LXv37qWpqYnp0/XfLVlZWXz88ccdztm5cycmk6kfPt3wFHcAs23bNpYtWxb9vn2uyuWXX85f/vIXKisrKSsri+4fP348L7/8MqtWreK3v/0tY8aM4Y9//OPoXEJtMIItTX/lzjx+fzgIIR8EffowV/NhaCyFpnJwV7UFOm0Bjafm6DBY/X791X4bwAYUWlMpTC6IBjfhKWPxOQrxGFOpcCvsrQ9T7FIp8xg41KJR1hzGG9Q41ODlUIOX/+3Wh/5SbKZoMDMlx0lRRhJmo4qqgMNqJNlqItlqwmk1osoQlBBiBFq2bBkrV64kGAx2CEzOOOMMrrnmGgKBAMuWLcPpdHLDDTewatUqIpEIp556Ks3NzbzzzjskJydz+eWXd3p9k8nEtddey+9+9zuMRiPXXHMNJ510UjSg+cxnPsO9997LE088wZIlS/jb3/7Gxx9/3O2w1EgXdwCzdOnS47qxjtVZlt2lS5fy/vvvx3ur0cdg0l8Wp/592jgoOkV/HwroQU2w9ehXTx00l4PrSNurrWfHdaRt3k7bq2aXfnn0npsks4Ps9PHMSysilFlEa9E43LY8vCEDFW6NPQ2wuxE+aVQobtJobg2ytbSRraV6T41RVRifmcT0vGTmj01jcrYDVVVQFH3OTXsw47QaZehJCDEiLFu2jNbWVqZNm9ZhZe0ZZ5xBS0tLdLk1wO23305WVhZ33XUXBw8eJDU1lfnz5/PjH/+4y+vb7XZuuukmvva1r1FRUcFpp53Gn/70p+j+FStWcMstt3DjjTfi8/m48sorueyyy/joo48G7kMPcYrWXTQyRLhcLlJSUmhubu77JN6RJhJu64nxHO2RCXrB5+oY1DRXQFOp3qujRY6/jqJCcgFaWhGB1Im4kyfRZBuLK2jkQDPsboBP2gKb5kDHUx0WI3MLU5k/NpXZY1JJsuhxsaKA3Wwg1WYm3WHGYRn6Y91CiIHh8/koKSlh/Pjx0cRvYmTo7t92IH9/y2+U4U41gDVZf7Vrn2AcXWnVFtyEA/owVVOZPjTVeFD/2lCir6pqLkdpLsfCW1iADNVIJH0i09OmsCxjMk1jJ+IzpVHp1YOZHbWwrQbc/hBvH6jj7QN1qApMzXEyf1wa88amkZ9ixeMPU9HUisWkkm7XgxmnxSjzZ4QQQvSa9MCMJkHf0WXgvma9p0bT9FdrIzSWQMNBqN0LtXv04adPidgzCaRPxe2cRHPyVDxJY9ndpLKlGrZWQ5m74/E5DhMLxqVx+tQcxmUkRbebjQppdjMZSRaSbRLMCDHSSQ/MyCU9MGLgmaz6y5Gtfx8JHxPQpOvbCxbo+zQN3NVQuxtq9ugBTWMJqrcOq7cOK++QCUQsKYzPmMUpabNwT5zN4VAKW2v0YObDeqh2B3llVw2v7KphfKqJ06dkcvLUfJKtJqpdfqpdfkwGhVS7mSyHhRT76J1RL4QQInbSAyM6Cnj03pnWJr0HJnJMnZJgK9Tt04OZmk+gepe+auoYPuc4WjJm486YTYNjCjvqTWysgM3VEGqbemNUYEG+hTMmZzBrfB5G49GgxWY2kJNsIcthkQnAQowg0gMzciWqB0YCGNE1TdN7aFqb9CGmgFvf1i4c1HtoKnbAkfehoWMJ+YjBgidtBi2ZczmStoh1tSn8rxz2Nx89JsUMZ4wzc8akDAryctEMeqIng6qQ4TCTm2yNTgoWQgxfEsCMXBLAdEMCmCEiHGrrnWnUe2eCrR33tzZB5ftHA5pj5tBoKHjSZ+DKWcLHthN5rdrJ+sPQdLToOBNT4PwpNhZPziNiy9AnKANOq5GcZCsZSWbJMyPEMCUBzMglAUw3JIAZogJevWSCp04fejqWFtFXOFXsgLJ39aGn9l2KAXf6TJpylvCGspBXj9g7DDHl2OCLkxSWTc5AdWYTNuv/5majQrbTSnayBYsxvlomQojEkgBm5JIAphsSwAwDwVa9FIKnTp8U/GktVVD6FlrpWygNB6ObI4oRd+YcqjKX8LR7Pv86ZI3mmUk2w3nj4ZyJZmzJWQRtWUSMVlQFcpKt5KfaMBtlnowQw4EEMCOXBDDdkABmmAn59UDG2xbMfPo/sebDUPoWkZI3UZvLo5sjqoX63FN5wXgmfzg8lmqvvt1igM+NhS9OgPRUJ0FbNkFrBgaDgbwUK3kpVpnwK8QQJwHMyCUBTDckgBnGQgG9dpO7Sh9y+rTGUiIlbxE5+AZGT1V0szt1Gu8mncl9NSeyr0WfxGtQ4IwCuHAijE0zE7DnErDlYDQZ2wIZGwaZIyPEkNTdL7lNxfWD1o4lEzPiPueb3/wmTU1NPP/88wBUVVVxxx138PLLL1NRUUF2djZz587l+uuv57Of/WyX12loaODaa6/lxRdfRFVVLrzwQn7729/icDi6PKe4uJgbbriBt99+G7/fz+c+9zkeeOCBDuUMioqKOHToUIfz7rrrLm6++ea4P2tvSB4YMTIZzXoxyZQCvbyBu1oPaNqXZ6cVoaYVoc77OpGqjwjuehHzkc04mvZwVtMelplT+WTcZ7i/5TNsaEhn/WFYfxg+MybA5dPKyEiqIGDL5XAwl6pmHwVpNnKcVpnsK4QYEKWlpZxyyimkpqZy7733MmvWLILBIK+//jorV65kz549XZ576aWXUllZydq1awkGg1xxxRV897vf5e9//3unx3s8Hs466yzmzJnD+vXrAbjllls477zzeO+991DVoz3Pt912G9/5znei3zudzn76xEOXBDBi8LSXPEifoA8xuav0oAZAUVDzZmPJm03YU4dv18tYitdiCjQxp/pZHlOe50jeQv4cPIs/1k1n/WGFd47ABZPCfHliBQ5vJUFbFocC+RxpsjEmzUa20yIZfoUQ/er73/8+iqKwZcsWkpKOZhc/4YQTuPLKK7s8b/fu3bz22mts3bqVhQsXAvDAAw9wzjnncN9995Gfn3/cOe+88w6lpaW8//770d6Lxx9/nLS0NNavX8/y5cujxzqdTnJzc/vrYw4LMnFADD7VAM4cyJujZ/5NKQCDObrbkJRJ0qLLUb78GM2LVuNNm4aiRSho3MJP3f+PD9Nu4urkdwlEIvxjH3xnPfyvLILRU42zbieGur2UHqnhg8PNNLcGE/hBhRAjSUNDA6+99horV67sELy0S01N7fLcTZs2kZqaGg1eAJYvX46qqmzevLnTc/x+P4qiYLFYotusViuqqvL22293OPbuu+8mIyODefPmce+99xIKheL8dMOPBDAiscx2vUemcBFkTQPz0R8KBpOZlOmfwfz5e6n/zK9pHPNZwgYLya2HuSnwINtS1nCxbQsNfo37d8L1b8FHtRomXz2O+o9Qqnexp6yaAzVuguFOKnALIUQcDhw4gKZpTJs2Le5zq6qqyM7O7rDNaDSSnp5OVVVVp+ecdNJJJCUlcdNNN+H1evF4PNxwww2Ew2EqKyujx/3gBz/gqaeeYsOGDXzve9/jzjvv5MYbb4y7jcONBDBiaFAUcGRBwXzIOQFsqdFdRlUho3AKzmXXU/e5P1Az8cuEjXYy/OXco93PpuSfcLZpB8XNGjdvgv+3FY54wOhvIqn+I1xH9vHBoTpqW/xd318IIXoQ65qXq666CofDEX31VlZWFs888wwvvvgiDoeDlJQUmpqamD9/fof5L6tXr2bp0qXMnj2bq666il/96lc88MAD+P0j+2eezIERQ489XX/5W/Ql19560DSMqkJOZgaBJZdxZNI5WPe9QHrZa+QFSnnYcB8l1knc6vkyb1bNYmu1wufHw6VTIti1SiKtdRzyjqU2PZ8JWUlYTZIITwgRn8mTJ6MoSrcTdUGfUHvDDTd02Jabm0tNTU2HbaFQiIaGhm7nrpx11lkUFxdTV1eH0WgkNTWV3NxcJkyY0OU5ixcvJhQKUVpaytSpU2P4ZMOT9MCIocvihOzpkD8fnLmg6P+5mg0qhbnZWBd9k5Klv6Vu3OeJqGbGBw/whPluXkm6nQXs5vmD8P03YGctqJEgtuZiQod38HFxORVNrTH/NSWEEADp6emsWLGChx56CI/Hc9z+pqYmALKzs5k0aVL0BbBkyRKamprYvn179Pj169cTiURYvHhxj/fOzMwkNTWV9evXU1NTwxe+8IUuj925cyeqqh43ZDXSSAAjhj6zHTInw5gTIWUMqHrHodNqZEJBPpEFV7D/tPupH/s5IqqJGeE9PG25nWesd5DhK+Mn78FDH0JrCAxBD7b6XdQe/JCPDtXQ4pNJvkKI2D300EOEw2EWLVrEv//9b/bv38/u3bv53e9+x5IlS7o8b/r06Xzuc5/jO9/5Dlu2bOGdd97hmmuu4ZJLLomuQKqoqGDatGls2bIlet6f//xn3nvvPYqLi/nb3/7GV77yFVatWhXtWdm0aRP3338/H3zwAQcPHuTJJ59k1apVfP3rXyctLW1gH0aCyRCSGD6MZkgfrwcxTWXQUokKZDstpNjGUOm8krpx55JV8h9SKzZworaLly0/5tHQufz20JfYXmtl1RyYlQmm1jq0Iw3say4gPX884zIckjtGCNGjCRMmsGPHDu644w5++MMfUllZSVZWFgsWLODhhx/u9twnn3ySa665hs9+9rPRRHa/+93vovuDwSB79+7F6z2a9HPv3r2sWbOGhoYGioqK+MlPfsKqVaui+y0WC0899RQ///nP8fv9jB8/nlWrVrF69er+//BDjGTiFcNXwKsXjPQezeLZ3BqkyuVDcdeQs+9JUmr0v2SOkMWPA1ewMTKXL4yHy6eBtS18jxiskDWVSYV52MwyN0aIgSClBEauRGXilSEkMXyZ7ZAzA3JngUWf6Z9iMzEpy0FyZgEVc6/n0NwfErBmkk8tfzH/kgdNv2NTSSPXvgmfNOiXUcM+1OoP2bv3E2qaWxP4gYQQQsRKAhgx/NlSIX8eZE0FoxWDqpCbYmV8ZhLB/BMpXvJL6sadi6aofN7wHhusN/AZ31pueifCH3eBPwxoGiZXGZX7tnGgso5wZMh3TAohxKgmAYwYORzZembftCJQjdhMBsZnJpGc4qR6yqUcXPT/8CZPwEEr/8/0Z/5t/jl7Sg7xgzdgf5N+CWPAha90G5/sL8btH/mZLIUQYriSAEaMLKoKqYUwZiE4czEoCgUpem2kYOp4ShbdRuXUywkbbMxTD/Ci+Sd83f93bnnHz9py/RJKJIxat5fiT3ZQ2ehO7OcRQgjRKQlgxMhkMOlLr3NngtFKis3EhKwk7FYTDWNXcODke2nOXoRRiXCV8SX+ZbyFFz+o4KEPIdhWdcDoraVu72b2lVVIKQIhhBhiJIARI5stTS9P4MzDbFApyrCTnWwhbEvn8JzrOTT3BoLmVKaqh3nB/FOM5e9w87tQ1zaXVw37CB7+gD17duFqDST2swghhIiSAEaMfKoBMidB7kwUo40sh4WiDDtmo4I7az7FJ92FO/0EkhQ/vzM/xDda/sSP3gzwcfvqbE1DbSqjdNdmqhtbEvpRhBBC6CSAEaPHMb0xdrORCVkOUu0mwpYUDs1fQ834L6GhcKlxHY/xcx7eVM1/DkJ7piSD30XN3i2UVtZKGQIhhEgwCWDE6HJMb4zBZKMg1UZBmg1FVamd9BUOzbuRoMnJTLWUF8w/5tDuLdz3PvjaFiSpYR/u0h3sO3SYkMyLEUKIhJEARoxOx/TGpNpMFGXYMRoUPJlzOHjSnXhSppCstPJ/5vs5rfqv3PRWiMq22m1KJETwyMfsKS7GFwwn9nMIIcQoJQGMGL3ae2NyTsBu1RPfWU0qIWsGpQt/St24cwH4lvFVfhm4jbvfqmN7jX6qokVQavexb88umlulIKQQI11VVRXXXXcdkyZNwmq1kpOTwymnnMLDDz/coXbRcFNUVMT999+f6Gb0igQwQtjTIX8uZpuDoswknFYjqEaqp1xK2ZzVBA125qsHeEb9Mf/b+hHrDx891eAqp2T3Dmqah+8PMCFE9w4ePMi8efP473//y5133sn777/Ppk2buPHGG3nppZf43//+l+gmxi0QGP6rKiWAEQLAZIO8uRgcWRSm28hwmAFoyV5IyZI78TonkKa4+ZPpl+z/4F2eKz7m1NY6Kvdto7S6SSb3CjECff/738doNLJt2zYuuugipk+fzoQJEzj//PN5+eWXOe+88wAoKyvj/PPPx+FwkJyczEUXXUR1dXX0Oj//+c+ZO3cujz32GGPHjsXhcPD973+fcDjML3/5S3Jzc8nOzuaOO+7ocH9FUXj44Yc5++yzsdlsTJgwgX/9618djikvL+eiiy4iNTWV9PR0zj//fEpLS6P7v/nNb/LFL36RO+64g/z8fKZOncrSpUs5dOgQq1atQlEUFEUBYOnSpdHvj30de72hQAIYIdqpBsiejpJWRG6KjfxUK4oCQVs2pYtupSnnZExKmPtNDxHY+zqPfXJ0hZIx0EJLyXb2Hq6Ryb1CxErTIOAZ/Fccf2jU19fz3//+l5UrV5KUlNTpMYqiEIlEOP/882loaOCNN95g7dq1HDx4kIsvvrjDscXFxbz66qu89tpr/OMf/+BPf/oT5557LocPH+aNN97gnnvu4ac//SmbN2/ucN4tt9zChRdeyAcffMCll17KJZdcwu7duwEIBoOsWLECp9PJW2+9xTvvvIPD4eBzn/tch56WdevWsXfvXtauXctLL73Es88+y5gxY7jtttuorKyksrISgGeffTb6fWVlJRdccAFTp04lJycn5uc2GIyJboAQQ07qWDA7SKvdi8moUtHQSggTFbO+T9jkIOPwf7nN9Di/PdTCb3wX8oO5CkZVX6EUOryT3YHpTBlXgMVoSPQnEWJoC3rhzvzBv++Pj4C582Dk0w4cOICmaUydOrXD9szMTHw+HwArV65k+fLlfPTRR5SUlFBYWAjAE088wQknnMDWrVs58cQTAYhEIjz22GM4nU5mzJjBsmXL2Lt3L6+88gqqqjJ16lTuueceNmzYwOLFi6P3+8pXvsK3v/1tAG6//XbWrl3LAw88wO9//3uefvppIpEIf/zjH6O9KH/+859JTU1l48aNnHXWWQAkJSXxxz/+EbPZHL2uwWDA6XSSm5sb3Zaenh59/5vf/Ib169ezefNmbDZbbM93kEgPjBCdsadD3hwcScmMz0rCYlJBUamadjk1E74MwHXGZ1la/Wfu3BKJLrNWIiGU6l3sKSmjNSArlIQYqbZs2cLOnTs54YQT8Pv97N69m8LCwmjwAjBjxgxSU1OjPSWgT5p1Op3R73NycpgxYwaqqnbYVlNT0+F+S5YsOe779ut+8MEHHDhwAKfTicPhwOFwkJ6ejs/no7j46Hj3rFmzOgQvPXn11Ve5+eabefrpp5kyZUrM5w0W6YERoitmO+TNwVy3j/FKHeWNrXj8IWonXkDI7CRvz1/4hvF/pDW5+cWmq/nxYhNOs75CyVC7hz0RjclFY3FY5H8zITplsuu9IYm4b4wmTZqEoijs3bu3w/YJEyYAxN0rYTKZOnyvKEqn2yKR2Iei3W43CxYs4MknnzxuX1ZWVvR9V0Ngnfnkk0+45JJLuPvuu6M9OEON9MAI0R2DEXJmYEgfx7gMm75CCWgsPJPDs64lohj4vOE9fuS5l1vfbo3WUFK0CKb6vewtKafZK8usheiUouhDOYP9ahtmiUVGRgZnnnkmDz74IB6Pp8vjpk+fTnl5OeXl5dFtn3zyCU1NTcyYMaNPjwngvffeO+776dOnAzB//nz2799PdnY2kyZN6vBKSUnp9rpms5lwuGNvcV1dHeeddx4XXnghq1at6nPbB4oEMELEIm0cSvYMCjOSSLbpQYwr9yTK5t1ISLVwmuFj7gnewe1vuyhrK5ekaBEsDXvYf6icBs/wX7IoxGj1+9//nlAoxMKFC3n66afZvXs3e/fu5W9/+xt79uzBYDCwfPlyZs2axaWXXsqOHTvYsmULl112GWeccQYLFy7scxueeeYZHnvsMfbt28ett97Kli1buOaaawC49NJLyczM5Pzzz+ett96ipKSEjRs38oMf/IDDhw93e92ioiLefPNNKioqqKurA+DCCy/Ebrfz85//nKqqqujr04FOokkAI0SskjJRsmcwJj2JlLYgxpMxi0MLf0rA6GCuepD/i/yCX79Tx55G/RRFi2Bt3Etx2WFqXL4ENl4I0VsTJ07k/fffZ/ny5axZs4Y5c+awcOFCHnjgAW644QZuv/12FEXhP//5D2lpaZx++uksX76cCRMm8PTTT/dLG37xi1/w1FNPMXv2bJ544gn+8Y9/RHt27HY7b775JmPHjuWCCy5g+vTpfOtb38Ln85GcnNztdW+77TZKS0uZOHFidLjpzTff5OOPP2bcuHHk5eVFX8f2Lg0FijYMEle4XC5SUlJobm7u8R9DiAHX2ohW/QlHGt00efXZu2ZPBWO3343FX88RLZ3vhm/mypPGMDVNP0VTVLxp0yjIy6UgdWjN5BdiMPh8PkpKShg/fjxWqzXRzRlWFEXhueee44tf/GKim9Kp7v5tB/L3t/TACBEvWxpKzgnkp+nVrAECSQWULvo5rfYC8pUG/my4gz9srqa4WT9F0SLYG/dQUVnFofqux9GFEELERgIYIXrDloqSO5P8dAepSXoQE7JmcOjEn+F1jCNLaeYR5W5+vamZQ8fMibE37qG6pobiWrdk7RVCiD6QAEaI3rKmoOTMpCDNSVpbEBM2OymffyM+axZFajUPKvdwxyYvFW79lPYgpr6uln3VbiIRCWKEEN3TNG3IDh8lkgQwQvSFNRlyZ5Kfnkx6kp4gKmRJo3zBGoKmZGaqpfwy8ht+tilIZdvIUXsQ09xYx76aFglihBCiFySAEaKvLE7ImUleupOMtiAmYM+lbP6NhA1WTjHs4qfh33PLuxFq2opWR4OYpiYJYoQQohckgBGiP1gckDuL3PTkaCVrX/IEyuesbkt2t5mrQ4/zk00a9W2rqRUtgr1pL00uN3urJYgRQoh4SAAjRH8xJ7UFMSnR1UmejJlUzFyJhsJlxrVc6H+OH2+CJr9+ihoO6EGM2ydBjBBCxEECGCH6k9kOubPIT3eSZDmasbdq2uUA/ND0L05vXcdPNoGrLTmvIejF3ryPJk9AghghhIiRBDBC9DeTDSVnBmPS7XoVa6Ch8Cxqxn8JgNuNjzHVs4Vb3gN3W5kko78Za0sJTd6gBDFCCBEDCWCEGAjWFIxZUxibbsdo0AvH1U78Mg0Fn8GgaPzO9BAZrt38fDP428qLmL01mD1HJIgRQogYSAAjxEBx5mBOL2Rsuh1VVUBRqJx+Ja7sE7EoQf5ovg+aDvGr96E9VrG2lGH01UsQI0aXkrcG79UL3/zmNzvkYamqquLaa69lwoQJWCwWCgsLOe+881i3bl2312loaODSSy8lOTmZ1NRUvvWtb+F2u7s9p7i4mC996UtkZWWRnJzMRRddRHV1dYdjioqKUBSlw+vuu+/u1WcdTiSAEWIgpY3HlpJNYZoNRQEUlcMzV+JJnYZTaeXP5l+yt7KZv+w+eoqtuRhDoIUmb5A9VRLECDGUlJaWsmDBAtavX8+9997LRx99xGuvvcayZctYuXJlt+deeuml7Nq1i7Vr1/LSSy/x5ptv8t3vfrfL4z0eD2eddRaKorB+/XreeecdAoEA5513HpFIpMOxt912G5WVldHXtdde2y+fdygzJroBQoxoigJZU3GEfOSFIxxp8qEZzJTNvYHxW28l11PBb00P8o3iNeTZVc4uOrq82pM+k+ZW2FPVwrRcp96LI4RIqO9///soisKWLVtISkqKbj/hhBO48soruzxv9+7dvPbaa2zdupWFCxcC8MADD3DOOedw3333kZ+ff9w577zzDqWlpbz//vvRQoiPP/44aWlprF+/nuXLl0ePdTqd5Obm9tfHHBakB0aIgaYaIHsGaU4HWU4LABGTnfLZ1xM2WDjFsIvrjP/m9x/Dthr9FCUSwt60FyUSpLk1yL6aFqmdJESCNTQ08Nprr7Fy5coOwUu71NTULs/dtGkTqamp0eAFYPny5aiqyubNmzs9x+/3oygKFoslus1qtaKqKm+//XaHY++++24yMjKYN28e9957L6FQKM5PN/xIACPEYDBZIXs62ck2Uu16x2fAUUDl9G8DcJ3xOU5TPuDubXCwrYK1GmrF1rQftAiNniAHaqQApBCJdODAATRNY9q0aXGfW1VVRXZ2dodtRqOR9PR0qqqqOj3npJNOIikpiZtuugmv14vH4+GGG24gHA5TWVkZPe4HP/gBTz31FBs2bOB73/sed955JzfeeGPcbRxuJIARYrBYkyFzMnkpNuwWAwDNeafQMEbvBn7A/BCp4Xp+vgXqWvVTjAEXNtdBAOrcAQ7WeRLSdCEEMf8BcdVVV+FwOKKv3srKyuKZZ57hxRdfxOFwkJKSQlNTE/Pnz0dVj/76Xr16NUuXLmX27NlcddVV/OpXv+KBBx7A7/f3+t7DgQQwQgwmRzZq2ljGph3NEVM19Ru0OseTjJtHrb/F5Qvxiy3gbesBNrXWYfYcAaDG5edQvQQxQiTC5MmTURSFPXv2dHvcbbfdxs6dO6MvgNzcXGpqajocFwqFaGho6HbuyllnnUVxcTE1NTXU1dXx17/+lYqKCiZMmNDlOYsXLyYUClFaWhrzZxuOJIARYrClFWFwZuk5YlQFTTVRPuc6wkY7MznAzy1/56AL7tkO4baFBlZ3OYaAC4AjTT7KG7wJ/ABCjE7p6emsWLGChx56CI/n+D8kmpqaAMjOzmbSpEnRF8CSJUtoampi+/bt0ePXr19PJBJh8eLFPd47MzOT1NRU1q9fT01NDV/4whe6PHbnzp2oqnrckNVI06sA5qGHHqKoqAir1crixYvZsmVLt8fff//9TJ06FZvNRmFhIatWrcLn8/WqwUKMCJlTMduTKUizARC0ZVNxwtUAXKq8xhcMm9lWA498DJoGaBq2pv0oYb3+wOHGViqbWxPVeiFGrYceeohwOMyiRYv497//zf79+9m9eze/+93vWLJkSZfnTZ8+nc997nN85zvfYcuWLbzzzjtcc801XHLJJdEVSBUVFUybNq3D79Q///nPvPfeexQXF/O3v/2Nr3zlK6xatYqpU6cC+uTg+++/nw8++ICDBw/y5JNPsmrVKr7+9a+TlpY2sA8jweIOYJ5++mlWr17Nrbfeyo4dO5gzZw4rVqw4rmus3d///nduvvlmbr31Vnbv3s2f/vQnnn76aX784x/3ufFCDFuqCtkzcNht0ZVJLdkLqCs6D4D7LH9gglLJK4fguYNtp0SC2JoPtEU0UFrnpcYlfwgIMZgmTJjAjh07WLZsGT/84Q+ZOXMmZ555JuvWrePhhx/u9twnn3ySadOm8dnPfpZzzjmHU089lT/84Q/R/cFgkL179+L1Hu1h3bt3L1/84heZPn06t912Gz/5yU+47777ovstFgtPPfUUZ5xxBieccAJ33HEHq1at6nDdkUrR4lzWsHjxYk488UQefPBBACKRCIWFhVx77bXcfPPNxx1/zTXXsHv37g4ZCn/4wx+yefPm45aBdcXlcpGSkkJzc3N0LbwQI4KnDq3mE0rrvXj9YYiEKdp+B0lNe6g2j+V0120EMLNmAZzSlibCn5SP3zkW0NPMTMp2kOmwdHMTIRLP5/NRUlLC+PHjsVqtiW6O6Efd/dsO5O/vuHpgAoEA27dv75A8R1VVli9fzqZNmzo95+STT2b79u3RLrGDBw/yyiuvcM4553R5H7/fj8vl6vASYkRKykRx5DAmta1mkmrg8KxrCZmTyQmU8ee0v6ABv3ofituWV1s8RzD6GwG9M+ZAjZtGTyBxn0EIIRIgrgCmrq6OcDhMTk5Oh+05OTldrmP/2te+xm233capp56KyWRi4sSJLF26tNshpLvuuouUlJToq7CwMJ5mCjG8pE/EZLFS0FZuIGRN4/Csa9FQOLl1IzekbsQfgTu3QUtbnGJrPoAS1pdIahrsq26huTWYwA8hhBCDa8BXIW3cuJE777yT3//+9+zYsYNnn32Wl19+mdtvv73Lc9asWUNzc3P0VV5ePtDNFCJxDEbImITDbIzOh/Gkn0DNxC8D8P3AnznVeogqL9HCj0okjL1pH2j6MqWIBnurWnD7R372TSGEgDgDmMzMTAwGw3GVMKurq7tcx37LLbfwjW98g29/+9vMmjWLL33pS9x5553cddddxxWjamexWEhOTu7wEmJEs6eDM48shwWHVc/UWzf+fFoy5qBGgjxi+S0pqo+tNfD0fv0UQ9CDtaU0eolwRGNvlQtfMJyADyCEEIMrrgDGbDazYMGCDhNyI5EI69at63L5mNfr7ZAxEMBg0LOQSlp0IY6RPh5MNgpSbfp8GEWlYub3CVrScfireCznXwA8ufdozSSztwZTa230EoGQxu5KF8Fw538cCCHESBH3ENLq1at59NFHefzxx9m9ezdXX301Ho+HK664AoDLLruMNWvWRI8/77zzePjhh3nqqacoKSlh7dq13HLLLZx33nnRQEYIgV70MXMyRoPKmLb5MGGzkyPTvwXA/MZXuSr/ABpw7w6oaltpaXWVoAaPLrv0BSPsrWohHJE/EMTQI3+4jjyJ+jc1xnvCxRdfTG1tLT/72c+oqqpi7ty5vPbaa9GJvWVlZR16XH7605+iKAo//elPqaioICsri/POO4877rij/z6FECOFNQWS80lqriDbaaHa5cedNY+m3FNIrXqH6/yP8mbKHXzSbOTObXDvKWAxRLA378edPlMPgoAWX4j9NS1MzXGiKEqCP5QQYDKZAL1X3mazJbg1oj+1561p/zceLHHngUkEyQMjRpVIBCrfh4CXsgYvLb4QhoCLSe/+CGOwhZKxX+b80gtwBeDMQrhujp4PJmhNpzV1SodLZSdbmJjV+2JyQvSnyspKmpqayM7Oxm63S3A9zGmahtfrpaamhtTUVPLy8o47ZiB/f8fdAyOEGGCqCplToPIDClJtFNe5CZJM5bTLKfzoQcaVP8ed0xfxgw/GsLYcpqXB58aByddAyFtN0H40zUGNy4/ZoFKYbk/gBxJC177Yo6vM7WJ4Sk1N7bYg5UCRAEaIocjihJRCDE1lFKbZKanz4MpZgqvyXZLrdvCZqj9w+bSf8+c9Kg9/DOOTYWoaWFsOETanEDEezYZ5uLEVi1ElO1myn4rEUhSFvLw8srOzCQYlb9FIYDKZEjafVQIYIYaq1LHQ2oANN1lOCzUuP5XTryDp3d3Ymw/w3ZzX2ZN7Npuq4K5t8NvTIcUSweoqxps2Qx9XanOwzoPJoJKWZE7gBxJCZzAYZBGH6LMBT2QnhOglRdGHkhSVTIcZm8lAyJpB9ZSvAZBz4J/cNLWGgiSo9cE9OyCsgTHQgtl7pMOlNA3217hp8clfvUKIkUECGCGGMnMSpI1DQSE/zYqiQGPBMjxp01Ejfibt/yM/WahhMcAHdfC3PfppFncFatDT4VJ6orsWSXQnhBgRJIARYqhLLgCLA6vRoJcaUFSOTP82EdWEo+Fj5rjf5Po5+qH/PADbqkHRIthcxdFSA+2CYT3RXSAkie6EEMObBDBCDHWKAmnjAaJDSYGkvGitpNy9f+UzmU2cpx/Cb3ZCox8MQS8W9+HjLucLRthX3UJEEt0JIYYxCWCEGA5sqWDP6DCUVD/2HFqdRRhCXnL3/IUrp0ORE5oC8Ju2oo8WbyWGgOu4y7X4QhTXugf/cwghRD+RAEaI4SJ9PCjK0aEk1UDFCd9DUwyk1Gwho24LNy4Aswrba+HFEkDTsDUfhMjx817q3AHKG7zH30cIIYYBCWCEGC5MNnDqmS7bh5L8znHUFX0egLw9f2G81c23TtAPf2w3HGwGNezD6j7U6SUPN7ZS2+IflOYLIUR/kgBGiOEkdSyoxg5DSbXjv4Tfnocp0ETuvr9z7jhYlAOhiF700RfSq1Yb/U2dXvJgrRuXLK8WQgwzEsAIMZwYTHoQA9GhJM1gpuKE76KhkHZkI0lNe7h+DqRboMwNf/pEP9XqOogSOT5QiWiwT5ZXCyGGGQlghBhunHn6cBJHh5JaU6fSWLAUgJx9fyPFHGH1PP3wVw7BpipQwwGsrtJOLxkMa+ypaiEUluXVQojhQQIYIYYbVYW0IoCOQ0kTv0LYYMXuOkhy9XvMy4ILJuqn/G4n1PvA5KvH2FrX6WVbA2H2VbsZBgXqhRBCAhghhqWkTLDqpenbh5JCllTqis4DIGf/UyjhAJdNg4kp4ArCr9qWVttaSlDCnU/cbW4NcrDO0+k+IYQYSiSAEWK4Sp8QfZvpMGM3G6gfdw5BSzpmXx3p5a9jUuHG+URLDTxXDEokjLWl81VJADUuP0eaWgfjEwghRK9JACPEcGVxgiMbaBtKSrWB0UL1pIsAyCr5D4aAizEO+F7b0uon9sD+JjD5GjD6G7u8dFmDlwZPYKA/gRBC9JoEMEIMZ6njQNH/N7YYVTIdFprzTqXVOQ5DyEvWwecAOGssnJIHIQ1+uQNaQ+gTejtJcAd69eoDNW48/tBgfRIhhIiLBDBCDGcmKyTnR7/NdJgxmQxUT/k6AOmH/4fZU4miwLWzIdMKRzzwh12ghv1YPBVdXjoc0VcmSeFHIcRQJAGMEMNd6lgwmAFQFYW8ZBue9BNoyZyHooXJ2f8PAJxm+OE8UID/lsG7lXqtJDXYdTmBQEgv/Cgrk4QQQ40EMEIMd6ohmtwOwGk14rQaqZ78VTRFJbl2G/bG3QDMzoQvT9KPe+hDaPZpWFtKur18iy9EiaxMEkIMMRLACDESOHPBbI9+m5tiJeAcQ2PBMgBy9j0Jmj4UdOkUGNdWtfr/doEx0ILJW9Pt5atdfqpdvoFrvxBCxEkCGCFGAkWBtPHRb80GfUJvzYQvR5PbpVRtAsBkgFVzQVXgjQp9KMnqLuu0zMCxSuo8UjNJCDFkSAAjxEhhTwdrSvTbTIcZQ1IqdeO/AED2gadRwvrS6Mmp8OW2LL0PfQiu1lC3uWFAX5m0v7oFf0hqJgkhEk8CGCFGkpQx0beqopCbbKN+7Nkdktu1+9oxQ0mPfAym1joMAVe3lw+ENPZVuYlEZFKvECKxJIARYiSxp4M5Kfqt02rEkZR0THK756NByrFDSW8egXcqweYqic6V6YrbH+JgnXvAPoIQQsRCAhghRpqUwg7f5qZYceWfSquzCEOolayDz0b3HTuU9PsPocXTisVzpMdb1LYEpNyAECKhJIARYqRJygSTLfqt2aCS6bRRNeVSANIPr8N8TJDy6aEks+cIaqjnFUdlDV6avFJuQAiRGBLACDHSKAokF3TYlOkwE8qeiStzfofkdnD8UNK7FRGsru5zw0DbpN4aN76gTOoVQgw+CWCEGIkcOdHsvHB0Qm/1lPbkdtuxN+6J7p+cCl9pT3D3EXhamjG11vZ4m1BYY29VC2GZ1CuEGGQSwAgxEqlqhxpJoE/otWSMozF/KUC00GO7r07Wh5Ka24aSLC1lEOm5mKM3EOZAjUzqFUIMLglghBipnHmgGjtsyk2xUj/hfDTFgKPhI2zNB6L7Pj2UtOlwEKu7PKZbNXgCVMikXiHEIJIARoiRymAEZ06HTWaDSkrWGJryTgUg8+DzHfZPToWLjhlKam2qRg3FFpiUN3hp9kqmXiHE4JAARoiRLHkMKB3/N890mGmedD4aCsl1O7C2lHbYf8kUKGobSnr4o7ahpBhoGhyolUy9QojBIQGMECOZ0QyO7A6bVEUhLXc8zblLgON7YUwqXD9XH0p66whsKWnsMUNvu0BIY3+1ZOoVQgw8CWCEGOmSC/Sl1cdIsZlwT/mSvrtmKxb34Q77jx1KeuRjCDeU6F0sMWjxhTjU4O1zs4UQojsSwAgx0pntYM84bnNawWRc2SeioJFZ+sJx+y+eDAVJ0OCHv37YislXF/Mtq5p91Lb4+9RsIYTojgQwQowGxxR5bJdkNuKddqG+u/IdTN7qDvvNBrhmtv7+1UNQUnYIIrHPbymp8+AN9LwMWwghekMCGCFGA4sTbKnHbU4fO4OWzLkoaGSV/Oe4/bMz4cxC0IAH3w9hcPdcJ6ldOKInuQuFuy8OKYQQvSEBjBCjxafKC4C+rNo/4ysApFa+han1+GGiK2dAshkOtcBLHx5BCcde/8gXjFBc6+l9m4UQogsSwAgxWtjTweI4bnPquNl40k9A0cJkHHrxuP3JZvjOCfr7f+zVaKgsjeu2kuROCDEQJIARYjRJPn4ujFFVCM68CIC0io0Y/Y3HHbOsAOZlQjAC/7e1ASUQX+kASXInhOhvEsAIMZokZYLJdtzm5HHz8KZORY0EySh9+bj9igIrZ4NZhQ/q4L1dxXHdVq9cLUnuhBD9RwIYIUYTRel0LoyqqoRn6XNh0g+v6zRxXV4SfHWK/v6xD1vxNvdcrfpYwbCe5E6LMZ+MEEJ0RwIYIUYbRw4YzMdvHr8IX/IE1IifjLJXOz31gol6mQFXAP7+XmnMye3atfhCHKqXJHdCiL6TAEaI0UZVjyvyCKAoKtpsfS5Metl/UYPHz3MxqnDtHFCAdWVh9pQcivv2lc0+Gjyxr2QSQojOSAAjxGjkyD2uvACAbcLJBJxjMYRbySj7b6enTkuDc4r094++V0UgEH8wUlzrxheU+TBCiN6TAEaI0chkBVva8dsVFaK9MK+ihjpf/nz5NMiwwhEPvLh1X9y3D4Wl6KMQom8kgBFitHLkdrrZPOE0gkl5GEMe0g7/r9NjkkzwvZn6++d3u6moa4r79m5/iDIp+iiE6CUJYIQYrezpYLQcv101oLT1wmQeehkl3HlRxpNzYXEOhDT401sHiPRidZHMhxFC9JYEMEKMVooCjuxOdxknLSNsz8YYcJFWsbHL06+eBTYD7K4LsfGTil41Q+bDCCF6QwIYIUazLibzohpRZl4AQHr566B1XpAxywbfmKa//8fWCly++LPtynwYIURvSAAjxGhmsoI1tdNd6qTPEjHZsXircNR90OUlPl8E45PBHdR4+r2SXjXD7Q9xSObDCCHiIAGMEKOdM6/z7SYbyuSzgLZemC4YVPj+LP39+v0N7K9u6VUzqpp91Ls7n28jhBCfJgGMEKNdV5N5AWXa59FQcNZ/iNnT9RyXGemwvFB//+d3ins9HHSwziPzYYQQMZEARojRrpvJvDhzoXARABllXffCAFwxHZKMUFLvY92e6l41RebDCCFiJQGMEKLrybyAMv18AFKPvIUa9HR5iVQLXNY2offprWU0t8Y/oRdkPowQIjYSwAghup3MS+4stNRxqBF/l0uq251dBBOTwROI8I8tZb1uTpXkhxFC9EACGCGErqvJvIqCMv0LAKSX/7fLJdUABgW+P1t//8a+Wvb1ckIv6Plh/CGZDyOE6FyvApiHHnqIoqIirFYrixcvZsuWLd0e39TUxMqVK8nLy8NisTBlyhReeeWVXjVYCDFA7OlgMHe+b8IZYHFi9tXirN3e7WWmpcFZY/X3j71TQriX81na58NovcjwK4QY+eIOYJ5++mlWr17Nrbfeyo4dO5gzZw4rVqygpqam0+MDgQBnnnkmpaWl/Otf/2Lv3r08+uijFBQU9LnxQoh+pCjgzOl8n9EKk1cAkN7DZF6Ab04DhwkO1Xv53+7eTegFaPGFONzYeUFJIcToFncA8+tf/5rvfOc7XHHFFcyYMYNHHnkEu93OY4891unxjz32GA0NDTz//POccsopFBUVccYZZzBnzpw+N14I0c+6KPAIwLRzQVFxNH6CpaX7+S0pFr1iNcA/t5XT5O39fJaKptZeTwgWQoxccQUwgUCA7du3s3z58qMXUFWWL1/Opk2bOj3nhRdeYMmSJaxcuZKcnBxmzpzJnXfeSTjc9di23+/H5XJ1eAkhBoHJCra0zvclZcHYkwHI6CaxXbsV42BSCngDYf7ehwm9mgYHatwEw13PvRFCjD5xBTB1dXWEw2Fycjp2M+fk5FBVVdXpOQcPHuRf//oX4XCYV155hVtuuYVf/epX/L//9/+6vM9dd91FSkpK9FVYWBhPM4UQfeHsphdm+nkApFS+jSHQ/R8WBkXP0KsAb+2vY09l7/8QCYQiHKztegm3EGL0GfBVSJFIhOzsbP7whz+wYMECLr74Yn7yk5/wyCOPdHnOmjVraG5ujr7Ky8sHuplCiHb2jK4n82bPgPSJqJEgaRUberzU1GMn9L5b2usJvQANngBVzb5eny+EGFniCmAyMzMxGAxUV3eclFddXU1ubud/teXl5TFlyhQMBkN02/Tp06mqqiIQ6Hxc3GKxkJyc3OElhBgk3WXmVRSYoS+pzjj8P4iEerzc5dPAaYbyBi///aTzntpYHar34PH3fE8hxMgXVwBjNptZsGAB69ati26LRCKsW7eOJUuWdHrOKaecwoEDB4hEjo5f79u3j7y8PMzmLv7KE0IkVlc5YQCKTgdrKkZfPck123q8VIpFX5UE8My2wzT2YUJvRIP9Ne4+9eQIIUaGuIeQVq9ezaOPPsrjjz/O7t27ufrqq/F4PFxxxRUAXHbZZaxZsyZ6/NVXX01DQwPXXXcd+/bt4+WXX+bOO+9k5cqV/fcphBD9y2QFW2rn+wwmmHo2AFmHe57MC/ow0pQ0ldZgmL9v7v2EXoDWQJjSepkPI8RoF3cAc/HFF3Pffffxs5/9jLlz57Jz505ee+216MTesrIyKisro8cXFhby+uuvs3XrVmbPns0PfvADrrvuOm6++eb++xRCiP7XXS/MlLNBNWJt3IvVdbDHS6kKfH9mBAV4+0BdnzL0AtS4/NS5/X26hhBieFO0YZDm0uVykZKSQnNzs8yHEWKwaBqUb4FwF0M+b90HBzfiHnM6h6ZfFdMl7//QyNpDIcZnJvH/zp+JqnZeQDIWRoPCrIIUrCZDzwcLIRJiIH9/Sy0kIUTnupvMCzBNX1KddORdjIHmmC75zakhbCaVkjoPG/fV9ql5obDGgRopNSDEaCUBjBCia0lZXe/LmgpZU1EiIXKrel5SDZBqga9O1yfvP721rM8riqTUgBCjlwQwQoiuWRxgtne9f5q+pDr50FoULbZg5PxCH/kpFly+EM/uONznJlY0teLySakBIUYbCWCEEN3rrhem6BSwpaP4Gsmt774qfTujCt+ZrffCvL6rmoo+9qC0lxoISakBIUYVCWCEEN1L6mYejGqEqecAkHp4XdfHfcqiNDcLC5MJaxpPbCrt8zwWfzBCSZ0srRZiNJEARgjRPZMVLM6u908+ExQVtXY36aEYM+1qGlfONmJUFT6saGb7ocY+N7POHaCmRUoNCDFaSAAjhOhZd6uR7Bkw5kQAMivfiPmSYw1NnDtTv+5f3ztEINT3IaDSOi++YNeV7oUQI4cEMEKIntkz9WXVXZm8AgBTyQacptiGgxQtwkWTVdLsJmpa/LzyUWXPJ/UgHJGl1UKMFhLACCF6ZjSDNbXr/QUL9J4Yv4vc5vdjvmxKqIavnVgIwPM7K6jvh+y6srRaiNFBAhghRGy6W42kGmDScgDMB9fisBpjuqQSCbMsP8iUHAf+UIR/bOlbnaR2srRaiJFPAhghRGySMkHp5kfG5LMABSo/IJuGmC9raa3mm0uKUIB3iuvZW9W3OkkgS6uFGA0kgBFCxEY1gD296/2OHMifC4CtdD1Jlth6YdRwgCmOVpZN0yf0Pr6plEik73NY/MGIVK0WYgSTAEYIEbvucsJAdDIvB/5HVlLsRRYtniNcvLAQu9lASZ2HDftq+tDIo2pbAtS2SNVqIUYiCWCEELGzpenJ67pSuBisKdDaQFLN+9gtsQUxaqiVdKWFLy8YA8DTW8v7XCepXWm9R5ZWCzECSQAjhIidqupzYbpiMMHEz+rv979OlsMS86XN3krOnJFDQaqNFl+If/VDnSSQqtVCjFQSwAgh4tNdAANtk3mBiu04Qk3YzbH1whgDLVhCbi4/uQiA/+6q6nOdpHYtvhAVTbK0WoiRRAIYIUR8rKlgMHe9P2UM5MwELQIH1pLljL0XxuI5wqyCFBaMSyOiwV/fK+1zc9sdbmzF3U/DUkKIxJMARggRH0WJoRembTLv/rU4zGrsvTD+JtSgl0sXj8WgKnxwuJn3y/peJwmOLq3ujxVOQojEkwBGCBG/7mojAYw7GcxJ4KmBI++TGcdcGIv3CHkpNs6emQvA3947RCjSP/lcWgNhDjV4++VaQojEkgBGCBE/ixNMtq73Gy0wYZn+fv9/cVqN2Eyx9cKYfPUoYT9fmldAss3EkWYf/91V3Q+N1lU1+2jyBvrtekKIxJAARgjRO92VFoCjw0jlm6G1icxY58JoGhZPJXazkYsX6nWSnt1xuF9LAxTXuglKll4hhjUJYIQQvdNTAJM+HjKnQCQExetIthqxmmL7kWNqrYFIiKVTshiXYccTCPPMtv5ZVg0QCGmU1EmWXiGGMwlghBC9Y7aDxdH9MdHJvP8FTYt5LoyiRTC31qKqCpcvKQJg3Z5qyvpx/kq9O0BNi6/frieEGFwSwAgheq+nXpjxp4HRCq4KqN5Fss2IJcZeGLO3GjSN6XnJLB6fjqbBXzeV9mtCutI6r2TpFWKYkgBGCNF7SVn6suqumOww/nT9/f7XUVDIcHSTQ+YYatiHMdAEwKWLx2IyKHx8xMX2Q/2zrBogHJEsvUIMVxLACCF6z2gBS3L3x7QPIx16B/xuUm0mTMZugp5jmD1VAGQ5rZw7Kw+Av20+1K8TcFt8IY40y1CSEMONBDBCiL7pKSdM5hRIK4JwAA5uQEEhMym2uTDGQDNqSC8BcP7cAlLtJqpdfl77uKqPje6ovMErWXqFGGYkgBFC9I09A5RufpQoyjGTeV8HTSPVbsJoiLEXxlsJgNVk4JITxwLw3PsV/ZrLRbL0CjH8SAAjhOgbgwlsqd0fM2GZXj+psRTq9qEqCulJsc2FMbXW6UuxgdMmZzIxK4nWYJh/9uOyapAsvUIMNxLACCH6rqfVSBYHjDtFf7//dQDS7WYMas+9MPqS6hoAVEXhsrZl1Rv31vR7LhfJ0ivE8CEBjBCi72zp3Q8jAUw+S/9a+haE/BjU2HthzN4afZwHmJLj5OSJGWjAE/28rBqguNZDSLL0CjHkSQAjhOg7gxGsKd0fk3MCOHIg2Apl7wGQnmRGjaEX5tgl1QBfWzQWs0FlT1ULm0sa+tLy4wRCEcnSK8QwIAGMEKJ/JGV2v19RjxZ4LF4HgFFVSLOZYrp8+5JqgAyHhfPm6Muqn9x8iECof3tM6twB6tz+fr2mEKJ/SQAjhOgf9ozuk9oBTPyM/rVyJ3jrAUh3mHs8DdqWVAePTrI9b04+6Ulm6twBXvmospeN7lpJnQd/SLL0CjFUSQAjhOgfBlPPSe2S8yF7BmgROLgRALNBJcUeYy9M69FeGIvRwNcW6cuqn99ZQYOnfyffhsIaB2tlKEmIoUoCGCFE/+lpGAmO9sIUr4tOzM1yWGLqhTl2STXAyRMzmJztwB+K8PTWst60uFtN3iBVkqVXiCFJAhghRP+xZ/R8TNGpoJqgqQwaigG9FybZauzx1GOXVAMoisLlJxcB8Ob+Oopr3b1qdnfKGry0BmQoSYihRgIYIUT/MVrA4uz+GLMDxp6kv2+bzAuQ6YytvEB7lep2E7McnDZZ7/kZiGXV4YhGca0UfBRiqJEARgjRv2LphWkfRjr4RnRIyGo04IyhF0YN+zH6O1akvuTEsViMKvuq3bxbXB93k3vS4gtR0dTa79cVQvSeBDBCiP4VyzyY/PlgTQW/Cw5vj27OdMSa2K66w/fpSWa+MCcfgH9sKRuQ1UOHG1ul4KMQQ4gEMEKI/mWygdne/TGqASYs1d8fM4xkNxtJsvTcC/PpJdUAn5+dT6bDTL0nwEsf9v+yain4KMTQIgGMEKL/2WNZjfRZ/evhLeBviW6OuRfmmCXVAGajytcWjQPgxQ+OUD8AiehaA2HKpOCjEEOCBDBCiP4XyzyY9PGQPkGfA1PyZnSzw2LEbjb0ePqnl1QDnDQhnak5TvyhCP/YWh53s2NR2eyj2RsckGsLIWInAYwQov9ZHPpQUk8mHJMT5hgZjp5XJH16STXoy6ovWzIOBXjnQB37q1s6P7mPDtS6peCjEAkmAYwQYmDE0gsz4Qy9RlLdPmg+2mOSbDViMfX84+nTS6oBJmQ5OH1KFgCPbyolMgDLnwOhCKX1kqVXiESSAEYIMTBiCWBsaVCwQH9fvL7DrswYemHUsB9joPm47ZecWIjVpFJc6+GdA3UxNTdetS2BAZlnI4SIjQQwQoiBYU3WE9v1pH0yb/EGiBxd/pxiM2Iy9lxfwPSpJdUAqXYzX5xbAOjLqn3BgcmkW1Ln6fdK2EKI2EgAI4QYOLH0whQuAnMSeOug+qPoZgUlpl4YU6AJJXx8T8jZM/PIdlpo9AZ58YMjcTU7VsGwxsG6/i9fIITomQQwQoiBE0sAYzBD0en6+wMdJ/Om2kwYDT30wmga5tba4zabjSpfW6xXq37xwyPUtgzMcE+jJ0iNSwo+CjHYJIARQgwcawoYTD0f115aoOxdOCZBnaoopCf1nBfG1Fpz3GRegEVF6UzPcxIMa/xjS/9Xq25XWu8dsGEqIUTnJIARQgwcRYmtFyZrGiTnQ8gPhzZ12JVuN2NQu++FUcMBjIGmTm6vcNmSIhRg08F69lS64mh87KTgoxCDTwIYIcTAiiWAUZRjJvN2HEYyqLH1wny6PlK7oowkPjMtG4C/bCodsFIArtYQlc0ylCTEYJEARggxsKypoPZc34gJy/SvVR+Cu2OCuvQkM0oPU2GMgeZOJ/MCXLSwELvZwKF6Lxv21XR6TH8ob/DiDUjBRyEGgwQwQoiBpapgT+/5OEc25M7W3x/c0GGXUVVIs/fQC6NpmL2dByfJNhMXzh8DwNNby/EMUFXpSFvBRxlKEmLgSQAjhBh4sQwjQcdhpE8FARmOnnthTL7aTifzApx1Qg4FqTZafCH+veNwbO3pBY8/zOHG1gG7vhBCJwGMEGLg2dL0kgE9GbdET37nOgK1ezrsMhtUUmzdr2hSwwGM/sZO9xlVlcuW6NWq/7urmooBDDIqmlpp8UnBRyEGkgQwQoiBpxr0IKYnJjuMPVl//6nJvKCXF+ipF8bc2vlkXoDZY1JZMC6NsKbxxKbSARvq0dqGksIDNGFYCCEBjBBisCRlxnbcpLZhpNK3IdyxF8NiVHFau58QbPR3PZkX4OuLx2FUFT6saGZHWVNsbeoFXzBCWYO35wOFEL0iAYwQYnDY0umx+wQgZ5Z+bMANh7cetzuW8gJdLakGyE2xcs6sPAD++l4pwfDA1TKqavbR5A0M2PWFGM16FcA89NBDFBUVYbVaWbx4MVu2bInpvKeeegpFUfjiF7/Ym9sKIYYzg1FfUt0T1QATlurvD248brfNZCDJ0n0vjKm168m8AF+cW0Cq3US1y8+rH1X23KY+KK71EBrAIEmI0SruAObpp59m9erV3HrrrezYsYM5c+awYsUKamq6z61QWlrKDTfcwGmnndbrxgohhrlYh5Hac8Ic3gL+44slZjq7X1KtRoIY/Q1d7reZDXz1RL1O0nM7K2gcwF6SQChCSZ1nwK4vxGgVdwDz61//mu985ztcccUVzJgxg0ceeQS73c5jjz3W5TnhcJhLL72UX/ziF0yYMKFPDRZCDGO2GPLBAKSPh7QiiITg0NvH7XaYjdjNhm4v0VVOmHanTs5kUrYDXzDCUwNYJwmgzh2gzj0wxSSFGK3iCmACgQDbt29n+fLlRy+gqixfvpxNmzZ1ed5tt91GdnY23/rWt2K6j9/vx+VydXgJIUYAoxksztiO7WYYCSCjh7kwxkAzaqjr1P6qonB527LqN/fXcaCmJbZ29VJJnQd/SAo+CtFf4gpg6urqCIfD5OTkdNiek5NDVVVVp+e8/fbb/OlPf+LRRx+N+T533XUXKSkp0VdhYWE8zRRCDGWxJrUbvxRQoPpjcB8/KTfZasRi6v5HmKmbJdUAk7KdnD5ZH9b6y7ulRAYwg24orHGwVoaShOgvA7oKqaWlhW984xs8+uijZGbGOPYNrFmzhubm5uirvLx8AFsphBhUsZQVAH2+TO4s/f3BNzo9pKcVSebWWtC6n0B7yaKxWE0qxbUe3t5fF1vbeqnJG6RKCj4K0S/iCmAyMzMxGAxUV3f8q6a6uprc3Nzjji8uLqa0tJTzzjsPo9GI0WjkiSee4IUXXsBoNFJcXNzpfSwWC8nJyR1eQogRwpwEJltsx7ZP5j24vtNVRSk2IyZj10uzlUgIo6/rybwAaXYzX5qn10n6x5YyWgMDO8xT1uAd8HsIMRrEFcCYzWYWLFjAunVHM2RGIhHWrVvHkiVLjjt+2rRpfPTRR+zcuTP6+sIXvsCyZcvYuXOnDA0JMVrFkpUXYNwpYDBD82FoOP4PHgUlhl6YnqtPnz0zl9xkK02tQZ7fWRFb23opHNEorpWCj0L0VdxDSKtXr+bRRx/l8ccfZ/fu3Vx99dV4PB6uuOIKAC677DLWrFkDgNVqZebMmR1eqampOJ1OZs6cidncQ3VZIcTIFOs8GLMdChfr74s3dHpIqs2E0dB1L4wx4EINdV/3yGRQ+fpJ+oTeVz6qpLJ5YIsxtvhCUvBRiD6KO4C5+OKLue+++/jZz37G3Llz2blzJ6+99lp0Ym9ZWRmVlQObGEoIMcxZU8DQfWHGqPZhpJI3IHL80IuqKGQkdf/HkCmGXpj5Y1OZMyaFUETjiU2HBryHpKKpFbc/NKD3EGIkU7Rh0I/pcrlISUmhublZ5sMIMVLU7gV3z4EFkRD88zLwu2D5L6BgwXGHhDWNA9VuQl0UT9RUIy1Z83usiF3Z3MqN//qQUERj9ZlTOLEoxgnHvWQzG5hdkIKqxlBiQYhhaCB/f0stJCFEYsQ6jKQaoagtg3cXOWEMikK6o+temFgm8wLkpdj4/Gy9TtITm0oHPG9LayAsBR+F6CUJYIQQiWFL67FHJGriZ/SvZe9CsPO5I+l2M4ZuejLMrbUx3er8uQVkJJmpcwd4YeeR2NrXB5XNPpq9wZ4PFEJ0IAGMECIxVIM+FyYWmVPAmQchP5S91+khBlUhzd51L4wx0IwS7jmdv9Vk4BttGXpf/PAI1a6Bz9tyoNYtBR+FiJMEMEKIxIl1GElRjskJ0/lqJIAMh7nb+SQ91Udqt6gonVkFKQTDGo+/WxpbG/sgEIpQWi9ZeoWIhwQwQojEiTUrLxytjVS5E1obOz3EqCqk2bpe3WTy1XaaEO/TFEXhmycXYVAV3i9vYvuhzu/Xn2pbpOCjEPGQAEYIkThGC1gcsR2bnA9ZU/XSACVvdnlYusOM0kUnjBoOYAw0x3S7/FQb5846OqE3EBr4IR4p+ChE7CSAEUIkVqzDSAAT2ibzFq/v8hCzQSXF3k0vTAw5Ydp9aV4B6Ulmalr8vPDBwGboBb3gY3GNDCUJEQsJYIQQiRVPAFN0KigGvaxAU1mXh2U5LF32wpj8jSiR2Fb9WE0GvtGWofeFDwZnQm9za5AjTZKlV4ieSAAjhEgscxIYrbEda005msiui5ww0NYLYzN2vlPTMMW4pBpg8fh0ZuYnEwzrGXoHQ3mDF49k6RWiWxLACCESL57JvBPbVyNt1OfDdCGjmyKPseaEgfYJveMxKAo7yhrZUTbwE3ojGhyocRPpIrOwEEICGCHEUBBPADNmEZjs4KmB6k+6PMxqNOC0dt4Lo4ZaMQRaYr5lQZqNs2flAvD4u4MzodcbCFPeKFl6heiKBDBCiMSzpuolA2JhtMC4U/T33eSEAchydt0LE88wEsAF88aQZjdR0+LnpQ8HPkMvwJEmydIrRFckgBFCJJ6ixJkTpm0YqfRtCAe6PMxmMpBk6TwwMvnqOq1u3eW1zAa+3jah9/mdFdS2DPyEXpAsvUJ0RQIYIcTQYIsjgMmdCfZMCHrg8LZuD810dl5eQNEimHz18bSQJRMymJE3uBN6A6EIJXWytFqIT5MARggxNNjTYy/uqKhHM/P2MIzkMBuxmw2d7jPHkRMGjsnQqyhsOzQ4E3oB6twBalskS68Qx5IARggxNMRT3BGOBjCHt4LP1e2hXc2FMQTdqKH4cq4UptujE3r//E4JvuDgZM4trfcM2r2EGA4kgBFCDB3xzINJK4L0iRAJQelb3R7qsBixmTrvhYknM2+7C+ePIdNhps4d4N87Dsd9fm+EwhrFtW60GGo5CTEaSAAjhBg64snKCzCx59IC7TK76IUxt8ZW4PFYVpOBK04eD8ArH1VyaJAqSbtaQxxpHpzJw0IMdRLACCGGjniKOwKMP12fD1O3F5q77wlJtnbeC6NEQhj9DfG2lPnj0lg8Pp2IBn98u2TQks6VN3hxS5ZeISSAEUIMMfGsRrKlHS0tUNz9ZF7orhcm/mEkgMuWFGEzGThQ4+Z/u6t7dY14aW1ZesOSpVeMchLACCGGlt4OIx1c321pAdB7Yaym43/sGQMulHD8q3zSk8xcsqgQgKe2ltPg6TonTX9qDYQpHaRhKyGGKglghBBDi8URe3FHaCstkASeWqj+uMfDs5ydXFvTMLXWxdHIo5ZPy2FStoPWYJgnNpX26hq9UePyU++WpdVi9JIARggx9MSzGslogaJT9fcxTObtqhemt8NIqqrw7VPHoyqwuaSBHYcGJzcMwME6WVotRi8JYIQQQ088AQwcHUYqfQeCPa/SyeykUrUa9mPwN8d33zbjMpI4Z1YeAH9+d/Byw4TCGgdqZGm1GJ0kgBFCDD3xFHcEyJ4BjlwItUL5ph4PT7GZ+rUXBjrmhvnX9sHJDQPQ4gtxuDG+ZHxCjAQSwAghhp54izsqSlw5YaDzXhijvxEl0rvqz1aTgStP0XPDvPpx5aDWL6poasXlk6rVYnSRAEYIMTTFs5wajgYwlR+Ap+cJuSk2E5ZP9cIoWqTXk3kB5o1N46QJem6YP719cNByw7QvrZaq1WI0kQBGCDE0xVPcEcCZqw8laREoeSOmU7I66YUx+Wpjv2cn2nPDFNd6WDtIuWEA/EGpWi1GFwlghBBDU7zFHeGYYaR1MZUH6KwXxhD0ogbd8d33GGl2M19tyw3z9CDmhgG9anVNi5QaEKODBDBCiKEr3qR2RaeCaoKmMmg4GNMpnfXC9GUyL8Bnpx/NDfP4u6V9ula8Suu8srRajAoSwAghhq54AxizA8aepL8vXhfTKZ31wph89RDpfRCgKkdzw2wpbWBbafy1lnorHNHYXy1Lq8XIJwGMEGLoMprB4ozvnPZhpJI3IBJb0cNP98IokTCmXhR4PNa4jCQ+PzsfgD+9XYLbN3gFGN3+EOUNsrRajGwSwAghhrZ4e2Hy5+t5ZHzNULEjplOSbcbje2H6OIwEem6Y/FQrTa3BQS0zAHCkuZUm7+DNvxFisEkAI4QY2uLNyqsaYMIZ+vsYc8IoKMf1whgDLaihvvVimI0qV50+EUWBtw7UsX0QywxoGhTXugmEZGm1GJkkgBFCDG3mJDDZ4jtn4mf1r+WbwR/biqLOe2H6tqQaYHKOk3Pbygz88e2DuP2DN5QUCGnsr2mR+TBiRJIARggx9MXbC5M2HtKKIBKEQ2/HdEpnvTCm1tqYlmP35CsLCslPsdLkHfyhJFerlBoQI5MEMEKIoS/eeTCKAhPiKy0Ax/fCqJEgRn/fh33MRpXvnTERBXhrf92gVqwGvdRAs1dKDYiRRQIYIcTQZ0kGgzm+cyacoWfyrfkEXEdiOqWzXpi+5oRpNyXHGa1Y/eggDyVpGhyobZH5MGJEkQBGCDH0xVvcEfRem7y5+vuDG2I+LcVmwmYyRL83BppRwv747t2FixYWktc2lPTXQR5KCoQ0DtRIfhgxckgAI4QYHuIdRoKjk3mL1+s1kmKUnXxML4ym9anA47HMRpWr2oaS3txfx46ywR1Kam4NynwYMWJIACOEGB6sqfoS6XiMXayvYHJX60NJMXJYjNgtR+9l7ofVSO2m5Dg5u31V0luDO5QEMh9GjBwSwAghhgdV1YOYeBitMO5U/X0ck3kBcpxHe2HUsA9DwBXfvbtxcdtQUqM3yN/eO9Rv142FzIcRI4UEMEKI4aNXw0htq5FK34ZQ7JWa7WYjTqsx+n1/TeaFtlVJp+tDSW/sq+X9QR5KkvkwYiSQAEYIMXzY0/UJvfHIOQEcORD0Quk7cZ2anWyJ3s7oa4i5tlIspuY6OXtmLgCPvnUQzyAPJcl8GDHcSQAjhBg+DCZ9SXU8FBUmr9Df7389rlOtRgMpNr0XRtEimH39M5m33UUnFpKbrA8l/XWQh5KgbT5Mq8yHEcOTBDBCiOGlN8NIk5YfzQnTGF+gkOW0Rnth+qPA47EsRgPfO2NCwoaSNA0O1Mh8GDE8SQAjhBheehPA2NNhzCL9fZy9MGaDSppdT6JnCHpRg57479+NabnJfK5tKOmRNw8OegXpQEhjX7XUSxLDjwQwQojhxWTVCzzGa8rn9K/F6yEcX5CQ6bCgqno3TH9O5m13yYljKUyz4WoN8sgbxUQGOZho8YUorfcO6j2F6CsJYIQQw0+8WXkB8udBUhYE3HAovsm8JoNCelsvjMlXF1dSvFiYjSrXfmYyJoPCB4ebeX1XVb9ePxZVzT5qW/on47AQg0ECGCHE8NObYSTVAJPP0t/vi28YCSDTYcagKiiRMCZfffz370Fhup2vLx4HwN83l3Govn+HqmJxsNY96KuhhOgtCWCEEMOPxQlGS8/HfdqkM/XJvNUfQ/PhuE41qAoZjrZemH7MzHusM2fkMH9sGqGIxu/W78cfCg/IfboS0WBfdQuhsEzqFUOfBDBCiOGpN70wSZlQsFB/34temPQkM0aDgjHgQo0jKV6sFEXhe2dMINVu4kiTj79uGvyl1b5ghP2S5E4MAxLACCGGp97MgwGY0pYTpngdhOPLgWJQFLIces+PqbW6d/fvQbLVxPeXTgJg3Z4atpY0DMh9utPklSR3YuiTAEYIMTxZU0E19njYcQoW6r03fheUvRv36WlJJsxGRS/w2M+TedvNKkjh87P1go9/eOsg9e7Bn1x7uLGVBs/gLukWIh4SwAghhidF6V0vTB8n8yooenK7SGhAJvO2u3hhIeMzk3D7Q/x+YzGRyOAP6RTXumkNDO48HCFiJQGMEGL4svVyGGnSmYACVR+CqyLu01NsRiwmFbN3YIaRAIwGlWuXTcJiVPmk0sULHx4ZsHt1JRTWk9yFExA8CdETCWCEEMOXPV1fVRQvRzYUzNff7/9v3KcrKOQmWzEE3f2emfdYeak2vnlyEQDPbCvnQE3LgN2rK95AmOJa96DfV4ieSAAjhBi+VAPY0np3bntm3gPxT+YFcFiMOK3GAe2FAThjShYnTUgnosED6w8kZEin3h3gSJNM6hVDiwQwQojhLSmrd+eNOVEPfnxNUL65V5fITbFi9tdBZOCSvymKwrdPnUCmw0xNi58/v1syYPfqTlmDd9DrNAnRHQlghBDDmz1d74mJl2psmwtDrybzgl7oMcNuxOyr69X5sUqyGFm5bBKKAm/tr+PtAwN7v85oGuyvceMNSKZeMTT0KoB56KGHKCoqwmq1snjxYrZs2dLlsY8++iinnXYaaWlppKWlsXz58m6PF0KIuPRlGKl9NVLl+9DSu/pDWQ4LNn//F3j8tGm5yVwwrwCAP751MCGlBkJhjT1VLQQlU68YAuIOYJ5++mlWr17Nrbfeyo4dO5gzZw4rVqygpqbz/4E3btzIV7/6VTZs2MCmTZsoLCzkrLPOoqIi/pn/QgjRqd4OIzlzIW+e/r4Xk3lBLzGQa4tgCLh614Y4XDBvDDMLUvCHIvx67T7cCahb5A9G2FvVkpBl3UIcS9HizBe9ePFiTjzxRB588EEAIpEIhYWFXHvttdx88809nh8Oh0lLS+PBBx/ksssui+meLpeLlJQUmpubSU5Ojqe5QojRIBKB8vcg0osJrqVvwxt36704X/5zr5LjaWjsc1tpsE+M//5xavEF+fFzH1HnDjC3MJUfnTUVVVUG/L6fluU0MynbOej3FcPLQP7+jqsHJhAIsH37dpYvX370AqrK8uXL2bRpU0zX8Hq9BINB0tO7zt/g9/txuVwdXkII0SVV7X1OmMLFelbf1kY4vLVXl1BQKDB7UcIDP8nVaTWx+sypmAwKO8ub+PeO+IpS9pfalgCHG70JubcQEGcAU1dXRzgcJicnp8P2nJwcqqpiGz++6aabyM/P7xAEfdpdd91FSkpK9FVYWBhPM4UQo1Fvh5EMJpj0Wf39vtd6fXuH2UCW0tTr8+MxPjOJ75w2AYBn369gW+ng10sCKG9oTUiZAyFgkFch3X333Tz11FM899xzWK3WLo9bs2YNzc3N0Vd5efkgtlIIMSzZ0nq3GgmOTuat2AHu3k/IHWNsQmVw5oacNjmLFSfkAvD7jcUJy9NyoMadkLk4QsQVwGRmZmIwGKiu7pi4qbq6mtzc3G7Pve+++7j77rv573//y+zZs7s91mKxkJyc3OElhBDdUlW9SGNvJBdA7mxA6/VkXgALIQosgxdIfP2ksUzLddIaDPPrtfsSkuQuosHeKhf+kNRMEoMrrgDGbDazYMEC1q1bF90WiURYt24dS5Ys6fK8X/7yl9x+++289tprLFy4sPetFUKI7vR2GAmOycy7tneTgdvkKg2YjYMzqdaoqlz32cmk2U1UNLXyyBvFxLkuo18EQhp7q6RmkhhccQ8hrV69mkcffZTHH3+c3bt3c/XVV+PxeLjiiisAuOyyy1izZk30+HvuuYdbbrmFxx57jKKiIqqqqqiqqsLtltoaQoh+Zk3t1SoiAMYuAUsyeOuh7N1eN8EYcFHoGLxVQal2M6uWT8GgKmwpbeA/Hwx+0UcAjz/M/pqWhARQYnSKO4C5+OKLue+++/jZz37G3Llz2blzJ6+99lp0Ym9ZWRmVlZXR4x9++GECgQBf/vKXycvLi77uu+++/vsUQggBfRtGMphg2rn6+4+f1VPP9lIWjSRZejkfpxcm5zi5oq3o4z+3lvNBedOg3ftYjZ4gZQ2yMkkMjrjzwCSC5IERQsTM2wDVu3p3rq8Z/nUFhAOw4s62eTG9YDDRnDmPTyoHt6f5D28eZMPeGpIsBu744ixykrteLDGQxmcmkZuSmHuLoWXI5IERQoghz5am96b0hjXlaH2kj//d+zaEg6SEm8hwmHt/jV644pQiJmYl4fGH+c3afQmbWFtS56FOlleLASYBjBBiZFGU3g8jAcz4IigqVGyHhj5Ufm6pZGy6ncFMkmsyqKxaPoVkm4lDDV7+742DRBLUyX6gxi3Vq8WAkgBGCDHyJGX2/tzkPBh7sv5+17O9v46/BWvES36qrffX6IUMh4XrPjsZg6Kw6WA9/9hSNqj3b6dpsK/aTYsvmJD7i5FPAhghxMhjTe39MBLAzAv0ryVv9imxHS1VFKTasJkHb0IvwIy8ZL53hp6p96UPK3nt48oezhgY4YhevdobkER3ov9JACOEGHkUBex96IXJnKJP4NXCsPs/vb+OpxZVCzMhK6n31+il0yZncfGJehmWJzYdYvPB+kFvA0AorLG7sgVfUBLdif4lAYwQYmTqyzASwMwL9a/7Xgd/L1cTRcLgriLZaiIn2dK39vTC+XPyOXNGDhrw0MYD7KlMTGHcQCjC7koXgVAkIfcXI5MEMEKIkcmaAoY+rALKnw9pRRDywb5Xen8d1xGIRBibbsdsHNwfuYqi8M0lRSwcl0YwrHHff/cmrIK0LxhhT5WLUFiCGNE/JIARQoxMitK3XhhFgRPaemE+eUHPDdMbIT+4qzEaVMZnDv5QkqoqXPuZyUzJceAJhLn71T00eBKzOsjjD7O3uoWIlBwQ/UACGCHEyNWX5dQA40/T6yv5mqB4fe+v46oATSM9yTzouWEAzEaVG86aSn6KlXpPgHte25OwibWu1hD7pOSA6AcSwAghRi5rChj7MPdENep5YQB2Pdf7Io/BVvDUAlCUkYTRMIjJYdo4rSZuPnsaqTYTZQ1efr12X8KGcxo9QYprPQm5txg5JIARQoxcfU1qBzD5LDA79F6U8s29v05zOaD3hoxLt/etTb2U5bRy4+emYTWp7Dri4pE3ihOW6K62xU9JnQQxovckgBFCjGx9XY1kssHU9iKP/+59kceAFzz6UubsZCsptj7kqemD8ZlJevVqReGd4sQlugOoavZRXDu49aLEyCEBjBBiZOvrMBLA9M+DaoK6vVDTy0KREO2FAZiQlTSoZQaONXtMKt89/Wiiu1c+SkyiO4Aal58DNW6ZEyPiJgGMEGLk62svjC0NJn1Wf9+XIo/+FmhtBMBqMlCYoKEkgNOnHE1099f3DvH6rqqEtaW2xU9xrQQxIj4SwAghRr6+ZOVtN+NLgAKHt0JTH4Zdmg9H3+alWHFYjH1vWy+dPyefz8/OA+Av75YmtCemtiUgPTEiLhLACCFGPmsymPrY25FSAGOX6O/7UuSxtQl8ekZcRVGYkJWEkqChJEVR+NqisXxhTj6g98S89OGRxDQGqHMH2FctQYyIjQQwQojRwZnT92u0lxc4uBE8db2/zjG9MEkWI/kpg1ux+liKonDJiYV8aV4BAE9uLuOFnRUJa0+DJyDJ7kRMJIARQowOjhxQ+vgjL2sq5MyESAh2v9D763jrIXB0CfGYNBtWU+J+HCuKwkULC/nygjEA/GNrOc+9n7ggptETlCBG9EgCGCHE6GAw9X0yLxxT5PHVDkFI3JqOrkhSVYWJ2Y6EDSW1u3D+GC5aqE/s/ee2cv61/XAPZwycJm+QPVUthCWIEV2QAEYIMXo4c/t+jYIFkDpWz6770TO9v463Tr9Gm2SriTFpiRtKaveleQV8tW110r93HOaf28oTNieluTXIniqXBDGiUxLACCFGD2sKmPtYUFFRYf7l+vtPnterTfeGpnWYCwNQkGpLWIK7Y31hbgGXLh4LwHPvV/DU1sQFMa7WEJ8ccREISRVr0ZEEMEKI0cWZ1/drjFmk98REQrD1j72/jrtGr1bdRlEUJmU7MBsTPJYEfH52PpctGQfACx8c4e9byhIWxLj9IT4+0pywApRiaJIARggxujiyQTX07RqKAid+GxQDHN4CFdt7dx0tclwvjNmoMinLmfD5MABnz8zjipOLAD1j7xObDiWsdpI/GGHXERfN3mBC7i+GHglghBCji2qApOy+XyelEKafp7/f8iiEe/mL1V193LkpdhMFqYmfDwNw1gm5fOvU8QC8tquK+/+3D1+wl1W5+ygU1thd5aLG5UvI/cXQIgGMEGL06Y/JvABzvgrWVHAdhj0v9e4akbBe6fpTxqTZcFoTl6X3WMun57By2SSMqsLW0kZ+8eIuGjyBhLRF06C41kNZvTch9xdDhwQwQojRx+IAi7Pv1zEnwfzL9Pcf/CNa5yhurko9kDmGoihMznFgMgyBsSTg1EmZ/PTcGSRbjZTWe/np8x9xMIGVpCuaWtkvuWJGNQlghBCjU39M5gWYtBwyJkHQC+//tXfXiIQ6Xc1kMRqYkOXoYwP7z9RcJ7efP5OCVBuN3iC/ePETtpQ0JKw9de4An1S6CIZlhdJoJAGMEGJ0SsoCtR+GaBQVFn1Pf79/LdQf6N11mss7rEhql55kJi/F2ocG9q/sZCu3nX8Cc8akEAhH+M3/9vGfnRUJW6HU4gvxcUUzrYHEzMsRiSMBjBBidFJVfUVSf8ieDhOWAhps/j99oka8ImGoL+5019h0e0KrVn+a3WzkRyumseIEfS7RU1vLeeSN4oT1hPiCEXYdaaa5VVYojSYSwAghRq/+GkYCmP9NMFqhdjeUbOzdNbz14D1+SEZV9fkwBnVozIcBMKgK3zy5iCtOLkJV4M39ddz5ym5cvsQEEcGwxu5KF4cbZXLvaCEBjBBi9DLb9ey8/SEpE2ZdpL/f9ucOZQLiUl983IReAKvJwISsPmYRHgBnnZDLjSumYTMZ2FPVwi3Pf0xFYy8/ex9pGpQ3tLK7UjL3jgYSwAghRrf+7IU54YvgyIXWht7XSQr5oOlQp7syHRayky29b98AmVOYym3nn0C200JNi59b/vMx7xyoS1h7mrxBPqpokqR3I5wEMEKI0S0pEwzm/rmWwaxn6AXY9ay+PLo3XEe6rHQ9PiOJZNvQmQ/TbkyandvPn8m0XCetwTAPbjjAg+v3Jyz9fyCkJ70rb/AmbIKxGFgSwAghRjdF6b/JvACFiyFvnr40etufencNTetyNZOqKkzNcWI397EcwgBItpn46bkzuHD+GFQF3imu56Z/f8ieSldC2qNpcLixlU8qXfhDskpppJEARgghnHn0W/EhRYFF39GXV5e/B0fe7911fC5oqep0l9GgMi3Pidk49H6EG1SFLy8Yw63n6UNKde4At738CU9vLScUScy8FFdriI8ON9PkTUz2YDEwht5//UIIMdhMVr0kQH9JHQvTPq+/3/IHvTemNxpLu6yxZDEamJbrHFIrk441JcfJ3RfM5owpWWgaPL+zgp+/sIvK5sRM8NVXKbVwqN4jQ0ojhAQwQggB/TuZF2Du18CSrCeo++AfvbtGOAgNB7vcnWQxMjVnaFSu7ozNbOCqMyZy3Wcnk2QxUFzrYc2zH7F+T03CgogjTT52HXElbG6O6D8SwAghBIA9HYz9uMLH7IBF39Xff/g0lL7du+u4a6C1qcvdKXbTkFxefayTJmRwzwWzmZGXjD8U4dG3DvKb/+2jJUE5Y1p8IT483ExZvVdqKQ1jEsAIIQToc1f6q0p1uwlLYcYX9fdv/6bLTLs9qi+GbuaPZDutFKbbenftQZLhsPCTc6fztUVjMbRVtb7x3x/yzoG6hPTGaJpeEHLn4SaZGzNMSQAjhBDtHLn9N5m33YIrIH8+hP2w4f/1rmJ10Auuw90eMibNPiRzxBxLVRTOm5MfLQjZ5A3y4IYD/OLFTyip63zZ+EDzByPsrmxhf3WLJL8bZiSAEUKIdkYzOHL695qqAc64EZILwFMLG+7scmJut5rKe8zuOyEzibQkUy8bOnjGZyZx55dmcdHCQixGlb3VLfzkuY949K2DuBJUz6jOHeCDw01Uu3wJub+InwQwQghxrNRxetDRn8wO+MwtYErSayW99/v4Cz5qkR6HoBRFYXK2c0gVfuyK2ajypXkF/Oorczh5YgYasH5PDav+uZNXPqpMyJLrUFjjYK2HjyuaZZLvMCABjBBCHMtohpTC/r9uyhi9J0ZR4cBa2PNi/NdobQR3bbeHGFSFqblOLKbh8eM9w2Hh2s9M5tbzZlCUYccbCPPX9w5x878/4sPDTQlp07GTfEMJqrAteqZow2BBvMvlIiUlhebmZpKTkxPdHCHESBeJQMV2vS5Rf9v1nJ6hV1Fh+S8gf15856tGyJ0FFke3h/mCYT0DbXD4/AKORDQ27Kvh6a3ltPj0HpCF49L4+knjyEm2JqRNJoNCXqqN3GTrkM25M5QN5O9vCWCEEKIz7lqo3dP/19U0eOe3UPw/MCfBub/W58fEw2DWgxizvdvD/KEweypb8AaGVxp9jz/Ev3cc5vVdVUQ0MKoKZ0zJ4vOz88lNSUwgYzYq5KfayHFaUSWQiZkEMBLACCES4chO8Lf0/3XDQXh9jR4gpYyBc36lBzPxMFohb3aPuWuC4Qh7q1qiPRrDyeFGL09sOsRHFc2AvkBs8fh0vjCngPGZicl9YzaqjEmzke20oAzVDIJDiAQwEsAIIRLB54LKDwbm2q2N8NIq8NZBwUJ9km+8k4dNdj2IMXS/8igc0dhf00KjJzErfPpC0zT2VrXwnw+OsLO8Kbp9dkEKX5ibz4y85IQEEhaTyphUG1kSyHRLAhgJYIQQiVKzR1/+PBDqD8CrN0I4ADMv1HPGxMvigJxZYOh+5ZGmaRTXuqltGb5J2w7Ve3jxw0o2FdfRnkB3YlYS588pYEFRGmoCAgmrSWVMmp1Mh1kCmU5IACMBjBAiUYI+fUKvNkCTYUvegDfv1d8v+CaccGH8yfSsKZAzE9SeVx4dqvdwpGl45zqpcfl4+aNKNuytIRjWf4Xlp1j5/Jx8TpmYmZAq3WajSrbTQk6ydUhWCU8UCWAkgBFCJFJjqZ5IbqC8/1e9XhJA4UlwyvU9rjI6jj0dsmfEFPxUNLVSVu+Nv51DTHNrkNc+rmLtJ1V42iYq20wGTpqQwelTMtsKXQ5ur4iiQHqSmRynlRT70E8qONAkgJEARgiRSJEwHN6mD/UMBE2Dva/C1j9AJKRnAz7jZsicHN91HNmQOSWmIKamxcfBWk/c+fSGotZAmHV7qnl9VxV17qP/RtlOC6dNzuS0yVkJWYZtMxvISbaQ5bBgNIzOXhkJYCSAEUIkWksV1O0f2HvUH4CNd4O7Ss/3cuK3Yeq58Q0pOfMgc1JMhzZ4AuyvbmGkFGSOaBp7qlp4a18tm0saaA0eXT4+NcfJaVMyOWl8BkmDnKnYoCpkOMzkJFuHRZbk/iQBjAQwQohE0zQ48j4EBrjoYMCt54kp26R/P+5UOPkHPeZ86SC1ENKKYjrU5Qu2FTIc8r8K4uIPhdlW2sib+2v5qKI52tNkMigsGJfGoqJ0ZhWk4rAObkBhNxtITzKTnmQe9EAqESSAkQBGCDEUtDZB1UcDfx9Ng90vwLbHQAvrvSpL10D6hNivkVakBzIxCIYjHKz10OAZviuUutPgCfBucR1v7qulvPFoQUxFgUlZDmaPSWVuYSoTMpMGNUmdxaSSkWQmLclMsnVkzpeRAEYCGCHEUFH9CXjrB+detXvgjXv0ZdyqCRZ/DyaviH1IyZaqBz0xJsmrcfkorfcSHiljSp+iaRql9V7eLa7jg/KmDsEMgMNiZPaYFOYWpjKrIIVUu3nQ2mY2KqTZzWQkWUi2GUfMkmwJYCSAEUIMFcHWtmXVg/Sj0+eCd34Dh7fq309YCouvjj1zr6LoPTip43rMFQN6DaUDNe5hmbk3XvVuPx8ebmbn4SY+OtzcYc4MwPjMJKbnOpmY7WBilmPQsu+aDAopNhPJNhNOqxG7efgONUkAIwGMEGIoqS8G15HBu58WgV3Pwo4n9PdGK4w/Xe+NiXHVEQaTHsQ4c3s8XtM0KppaOdzYOiJWKcUiFIlwoMbNB+XNfHC4iZK64+c6OSxGJrUFM5Oyk5iQ5RiUoR+TQcFpNZFsM5JsNWE3G4ZND40EMBLACCGGknAIKrbpNY0GU/Uu2PQgNB+TkyZtPEz5nN4zE0uvjDkJMibqye964PaHOFDjpnWYFYPsD03eAB8fcbG/uoXiWjeH6r2EOhlay3ZamJjtYFy6nfwUm17wMXlgl00bDQpOqxGnta2HxmQYssu0JYCRAEYIMdSE/NBUBu7qwRtOAv1eNbtg3+tQ+jZE2oIogwXGn6b3ymRN67lXJikL0sf3WAwyEtE41OClqnl4Z+/tq2A4QlmDl+IaNwdq3BTXujnSxTNRFch2WslPtZLXFtTkp1jJT7XhtA7M/BaLScVuNpBkNmJr+2o1qQnvqRlyAcxDDz3EvffeS1VVFXPmzOGBBx5g0aJFXR7/zDPPcMstt1BaWsrkyZO55557OOecc2K+nwQwQoghK+DVA5mBqpfUHX8LHNwA+17T29AudZzeKzPmRD25ndLFX+eqQa+G7cgFY/cTVpu8AcoavHj8o683pisef4iDdR6Ka9wcbmrlSFMrlc2t+IJdl52wmY4uo05PMpNxzHv9ewtJlv4ZIlIVsLcFNHazAYtRxWrSvw5Wj82QCmCefvppLrvsMh555BEWL17M/fffzzPPPMPevXvJzs4+7vh3332X008/nbvuuovPf/7z/P3vf+eee+5hx44dzJw5M6Z7SgAjhBjy/G5oOgTehsG/t6bpK5b2vQalb3XMGGy0QMpYSBsHqWP14CZ1HNgzOvbSmOz6sFL7q4uAxhsIUdcSoNbtJxAaoPpQw5imaTR6g9Fg5kiTjyPNenBT7w4Qyy9cs0El1a4PDzks+lCRw2rEaTG2bTO1DSHp+60mA1aTAUMcS8CNBgWLUcViNGA16V8tRhWLScVkUDGqSr8EUUMqgFm8eDEnnngiDz74IACRSITCwkKuvfZabr755uOOv/jii/F4PLz00kvRbSeddBJz587lkUceiemeEsAIIYYNX7NeO8nnSsz9A244uBEOrIPGEr00QWfMSW3BzFh9OMlk17eZ7HrSPFsGOHPAma/30liOn1/T3Bqkzu2nwRMgFB7ysxESzh8KU+8O0OAJUO8J0Nj2tcEToMGjP0dXH1Z/mQ0qVrMBm0nFatSDGptZD1BsJgMmg4rZqGI2qNH3JoPS9vXodoOqYDIqx/TY6Ne0mQxYTO3XM2Ixqm3HGTC3Xe/TQc9A/v6Oa21WIBBg+/btrFmzJrpNVVWWL1/Opk2bOj1n06ZNrF69usO2FStW8Pzzz3d5H7/fj9/vj37f3NwM6A9CCCGGNgWSxoPSoE+29Q9w5t7O7l+wTH9FwnpZguZyvRhlcxk0HwZXpd5j1LILynfFdlnVpK9+Ug16mQNVRVGMZKlGMlUDYQyEUAljQEMBlLavHNPTo0TfR/eNMl0WeVAAB2hJENY0QpEI4Yh29KVpRNreRyIaoQiE0YhEQGvv1wm1vQZIa9urJwqK/q+rgLdtyHEgptvGFcDU1dURDofJycnpsD0nJ4c9e/Z0ek5VVVWnx1dVVXV5n7vuuotf/OIXx20vLIwtq6QQQgghho76+npSUnpe+RaPIZkdZ82aNR16bZqamhg3bhxlZWX9/gCGO5fLRWFhIeXl5TK8dgx5Ll2TZ9M1eTZdk2fTNXk2XWtubmbs2LGkp6f3+7XjCmAyMzMxGAxUV1d32F5dXU1ubm6n5+Tm5sZ1PIDFYsFiOX5pX0pKivzH0YXk5GR5Np2Q59I1eTZdk2fTNXk2XZNn0zVV7f9VT3Fd0Ww2s2DBAtatWxfdFolEWLduHUuWLOn0nCVLlnQ4HmDt2rVdHi+EEEII0ZO4h5BWr17N5ZdfzsKFC1m0aBH3338/Ho+HK664AoDLLruMgoIC7rrrLgCuu+46zjjjDH71q19x7rnn8tRTT7Ft2zb+8Ic/9O8nEUIIIcSoEXcAc/HFF1NbW8vPfvYzqqqqmDt3Lq+99lp0om5ZWVmHrqKTTz6Zv//97/z0pz/lxz/+MZMnT+b555+POQcM6ENKt956a6fDSqOdPJvOyXPpmjybrsmz6Zo8m67Js+naQD6bYVFKQAghhBDiWEOz+pMQQgghRDckgBFCCCHEsCMBjBBCCCGGHQlghBBCCDHsDPkA5qGHHqKoqAir1crixYvZsmVLops06O666y5OPPFEnE4n2dnZfPGLX2Tv3r0djvH5fKxcuZKMjAwcDgcXXnjhcQkER7q7774bRVG4/vrro9tG83OpqKjg61//OhkZGdhsNmbNmsW2bdui+zVN42c/+xl5eXnYbDaWL1/O/v37E9jiwREOh7nlllsYP348NpuNiRMncvvtt3eo1TJans2bb77JeeedR35+PoqiHFejLpbn0NDQwKWXXkpycjKpqal861vfwu12D+KnGBjdPZtgMMhNN93ErFmzSEpKIj8/n8suu4wjR450uMZofDafdtVVV6EoCvfff3+H7f3xbIZ0APP000+zevVqbr31Vnbs2MGcOXNYsWIFNTU1iW7aoHrjjTdYuXIl7733HmvXriUYDHLWWWfh8RwtErdq1SpefPFFnnnmGd544w2OHDnCBRdckMBWD66tW7fyf//3f8yePbvD9tH6XBobGznllFMwmUy8+uqrfPLJJ/zqV78iLS0teswvf/lLfve73/HII4+wefNmkpKSWLFiBT6fL4EtH3j33HMPDz/8MA8++CC7d+/mnnvu4Ze//CUPPPBA9JjR8mw8Hg9z5szhoYce6nR/LM/h0ksvZdeuXaxdu5aXXnqJN998k+9+97uD9REGTHfPxuv1smPHDm655RZ27NjBs88+y969e/nCF77Q4bjR+GyO9dxzz/Hee++Rn59/3L5+eTbaELZo0SJt5cqV0e/D4bCWn5+v3XXXXQlsVeLV1NRogPbGG29omqZpTU1Nmslk0p555pnoMbt379YAbdOmTYlq5qBpaWnRJk+erK1du1Y744wztOuuu07TtNH9XG666Sbt1FNP7XJ/JBLRcnNztXvvvTe6rampSbNYLNo//vGPwWhiwpx77rnalVde2WHbBRdcoF166aWapo3eZwNozz33XPT7WJ7DJ598ogHa1q1bo8e8+uqrmqIoWkVFxaC1faB9+tl0ZsuWLRqgHTp0SNM0eTaHDx/WCgoKtI8//lgbN26c9pvf/Ca6r7+ezZDtgQkEAmzfvp3ly5dHt6mqyvLly9m0aVMCW5Z4zc3NANHiWNu3bycYDHZ4VtOmTWPs2LGj4lmtXLmSc889t8Pnh9H9XF544QUWLlzIV77yFbKzs5k3bx6PPvpodH9JSQlVVVUdnk1KSgqLFy8e8c/m5JNPZt26dezbtw+ADz74gLfffpuzzz4bGN3P5lixPIdNmzaRmprKwoULo8csX74cVVXZvHnzoLc5kZqbm1EUhdTUVGB0P5tIJMI3vvENfvSjH3HCCScct7+/ns2QrEYNUFdXRzgcjmb4bZeTk8OePXsS1KrEi0QiXH/99ZxyyinRbMZVVVWYzebo/zjtcnJyqKqqSkArB89TTz3Fjh072Lp163H7RvNzOXjwIA8//DCrV6/mxz/+MVu3buUHP/gBZrOZyy+/PPr5O/v/a6Q/m5tvvhmXy8W0adMwGAyEw2HuuOMOLr30UoBR/WyOFctzqKqqIjs7u8N+o9FIenr6qHpWPp+Pm266ia9+9avRYo6j+dncc889GI1GfvCDH3S6v7+ezZANYETnVq5cyccff8zbb7+d6KYkXHl5Oddddx1r167FarUmujlDSiQSYeHChdx5550AzJs3j48//phHHnmEyy+/PMGtS6x//vOfPPnkk/z973/nhBNOYOfOnVx//fXk5+eP+mcj4hcMBrnooovQNI2HH3440c1JuO3bt/Pb3/6WHTt2oCjKgN5ryA4hZWZmYjAYjlsxUl1dTW5uboJalVjXXHMNL730Ehs2bGDMmDHR7bm5uQQCAZqamjocP9Kf1fbt26mpqWH+/PkYjUaMRiNvvPEGv/vd7zAajeTk5IzK5wKQl5fHjBkzOmybPn06ZWVlANHPPxr///rRj37EzTffzCWXXMKsWbP4xje+wapVq6IFaEfzszlWLM8hNzf3uEUVoVCIhoaGUfGs2oOXQ4cOsXbt2mjvC4zeZ/PWW29RU1PD2LFjoz+XDx06xA9/+EOKioqA/ns2QzaAMZvNLFiwgHXr1kW3RSIR1q1bx5IlSxLYssGnaRrXXHMNzz33HOvXr2f8+PEd9i9YsACTydThWe3du5eysrIR/aw++9nP8tFHH7Fz587oa+HChVx66aXR96PxuQCccsopxy2137dvH+PGjQNg/Pjx5Obmdng2LpeLzZs3j/hn4/V6OxScBTAYDEQiEWB0P5tjxfIclixZQlNTE9u3b48es379eiKRCIsXLx70Ng+m9uBl//79/O9//yMjI6PD/tH6bL7xjW/w4Ycfdvi5nJ+fz49+9CNef/11oB+fTe/nHg+8p556SrNYLNpf/vIX7ZNPPtG++93vaqmpqVpVVVWimzaorr76ai0lJUXbuHGjVllZGX15vd7oMVdddZU2duxYbf369dq2bdu0JUuWaEuWLElgqxPj2FVImjZ6n8uWLVs0o9Go3XHHHdr+/fu1J598UrPb7drf/va36DF33323lpqaqv3nP//RPvzwQ+3888/Xxo8fr7W2tiaw5QPv8ssv1woKCrSXXnpJKykp0Z599lktMzNTu/HG/9/OHaoqDAVgHN8tTkQEk0UWBMFgsfgCFpNgXBKrxXaL+AoWm8WH8AFWDJqO1WQ0WVwQDX63jatMMYjzsP8PlnY4nH3h7As7+43GpCWbMAxljJExRo7jaDKZyBgTnaR5JYd2u61Go6H1eq3lcqlqtSrf95N6pLd5ls3lclGn01G5XNZms7nZl8/nczRHGrOJc38KSXpPNl9dYCRpOp3K8zxlMhk1m02tVqukl/RxjuPEXvP5PBpzOp00GAxULBaVy+XU7Xa13++TW3RC7gtMmnNZLBaq1+tyXVe1Wk2z2ezm/vV61Xg8VqlUkuu6arVa2m63Ca32c47Ho4bDoTzPUzabVaVS0Wg0unnxpCWbIAhi95ZeryfptRwOh4N831c+n1ehUFC/31cYhgk8zXs9y2a32z3cl4MgiOZIYzZx4grMO7L5kf79fhIAAMACX/sNDAAAwCMUGAAAYB0KDAAAsA4FBgAAWIcCAwAArEOBAQAA1qHAAAAA61BgAACAdSgwAADAOhQYAABgHQoMAACwDgUGAABY5w/YOxkpyM05sAAAAABJRU5ErkJggg==", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "import numpy as np\n", "import matplotlib.pyplot as plt\n", "\n", "timeline = np.linspace(0, 140)\n", "weibull.plot(\"sf\", timeline, label=\"Weibull\")\n", "gompertz.plot(\"sf\", timeline, label=\"Gompertz\")\n", "plt.legend()\n", "plt.show()" ] } ], "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.11.0rc1" } }, "nbformat": 4, "nbformat_minor": 5 }