{ "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 assetlife.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 assetlife.lifetime_models import Gompertz, Weibull" ] }, { "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": [ "[None None]\n" ] } ], "source": [ "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": [ "Gompertz(shape=0.003907810339258283, rate=0.0757954646405569)" ] }, "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": [ "[3.7267452 0.01232326]\n", "[0.00390781 0.07579546]\n" ] } ], "source": [ "print(weibull.get_params())\n", "print(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": "iVBORw0KGgoAAAANSUhEUgAAAjEAAAGdCAYAAADjWSL8AAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjExLjEsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvctoD+AAAAAlwSFlzAAAPYQAAD2EBqD+naQAAjPpJREFUeJzs3Xd83HX9wPHX9/ZKLjtp0jTpnnTRQSnQllmQUZmioPwAcQAKDhRRlCFDUUEFVFBkioCooMhsC5RC6YAW6EyaZu9xye31/f3xTa5Nm3GXnfT9fDzukbvvfdelI+98Pu/P+62oqqoihBBCCDHK6Ib7BoQQQggh+kKCGCGEEEKMShLECCGEEGJUkiBGCCGEEKOSBDFCCCGEGJUkiBFCCCHEqCRBjBBCCCFGJQlihBBCCDEqGYb7BuIRjUapqqoiKSkJRVGG+3aEEEIIEQdVVWlrayM3NxedbuDHTUZFEFNVVUV+fv5w34YQQggh+qC8vJzx48cP+HlHRRCTlJQEaN+E5OTkYb4bIYQQQsSjtbWV/Pz82M/xgTYqgpiOKaTk5GQJYoQQQohRZrBSQSSxVwghhBCjkgQxQgghhBiVJIgRQgghxKg0KnJihBBCjA2RSIRQKDTctyEGiNFoRK/XD9v1JYgRQggx6FRVpaamhpaWluG+FTHAUlJSyMnJGZY6bhLECCGEGHQdAUxWVhY2m00Kl44Bqqri9Xqpq6sDYNy4cUN+DxLECCGEGFSRSCQWwKSnpw/37YgBZLVaAairqyMrK2vIp5YksVcIIcSg6siBsdlsw3wnYjB0/LkOR66TBDFCCCGGhEwhjU3D+efapyCmqKiIF154gbKysriP2bFjBy+//DK7d+/uyyWFEEKIMWHr1q187nOfi73esGEDF1544YCf9/DXY1FCQcy2bds444wzWL16NRdddBHvvPNOr8cEg0HWrFnDqlWruP/++1myZAlXXXUVqqr2+aaFEEKIwbR7925Wr159xC/rP/rRjzj77LOPmDq5/PLLeeKJJ+I6d319Pa+99lrsdU1NDW+++Wa/7/nw8x7+eixKKIhpbGzkxhtvZN++fXEfc//99/Pee++xfft23nrrLd5//32eeeYZnnrqqYRvVgghhBgKhYWFrF+/vlNwEY1Geeihh1i3bh2bN2+Oba+srOSpp56KO2l50aJFvPLKKwN+z0ejhIKY0047jdWrVyc0//Xkk09yySWXxFpwz549mzPPPJMnn3wysTsVQgghhojFYmHZsmWsW7cutm3btm0YDAYuuuiiTtvXrVuHXq/nxBNPBMDtdnPfffdxwQUXcMUVV/C3v/2t07lLS0t54IEHjrjmJ598wg033MDnP/957r77bvx+f+y9V155hSuvvLLT/jt27GD16tVEo9EB+cyj0aAusQ6Hw+zatYvrr7++0/a5c+fyhz/8odvjAoEAgUAg9rq1tRWA9x/8KmarnZDOREgxEVZMBHVmwoqJkGIm1P46pFgI6CyEdBaCioWQ3kxQsRLUmUHRln8pioKigE5RUNC+6nQACrqO7e1fDToFvb79q07X/lXp/FWvw2w49KHXvhp1mPR6zMaD2+1mPXaTAZ1OktyEEEcfVVXxhSLDcm2rUR/3L+KrVq3iT3/6U+z1unXrWLFiBStXruSpp57illtuiW1ftGgRycnJuN1uli9fTmFhIZdeeiler5dbb72VzZs38+tf/xroeprH4/Fw1lln8YMf/ID58+dzxx13sGXLFv7xj38AUFZWdkQKR1NTE6+99hrRaBSd7uhcpzOoQYzb7SYSiZCamtppe3p6eo9VG++++25uu+22I7Yva/0fyYH+/eAPqAZ8mPFixqNa8WDBrVrwYMWN5ZBt7V+x0qraaMXe6WsbVtR+Lu6ymfTYTHrsZgN2kwGH2YDNrL1OMhtItZtIt5tId5jJcJjITDKTZjeRajNh1B+df2GFEKOfLxRh1q3Dk6ux8/YzsJni+9G3atUqfvrTn7Jv3z6mTp3KunXrOPPMM1mxYgXf/OY3CQaDmEwm1q1bxyWXXALAAw88gMVi4V//+lcsWFq0aBELFizghz/8IVlZWV1eKxwO87vf/Y41a9bEjpk7dy6bN29m8eLF/f/gY9SgBjFmsxkAr9fbabvb7cZisXR73M0338x3vvOd2OvW1lby8/P5eOLVOE0q+mgQfTSAIRpAH9Ge66MB9JEAhqgfQ8SHIeJvf/gwRP0oaInEZiWMmTApeKAf8VAUhYDOhk+x49XZ8egctCpJuJRkWkiiSdUejWoS9REHdREHdVEHzSEjEVW7sDcYwRuM0OAOJnx9u1mP02IkxaYFNzlOC7lOK7kpFsanWBmfZsNpM2LU6TC0jyLJ8kYhhIjf0qVLsdlsrF+/nkmTJrFhwwbuvfdeJk6cSEZGBps2baKgoICSkhJOPvlkAN544w3q6uo455xzUFU19ohGo+zevbvbIEZRFM4888zY6zlz5lBQUMCmTZskiOnBoAYxVquVnJycI7K7y8rKmDRpUrfHmc3mWAB0qPlf+BnJycmJ34iqQtgPIR8EPRDyal+Dbgi427+2Hfwae7SCv/0RcLW/boNIAB0q1qgHKx7S4h0V1YFqM6FaUohYUgma0/Cb0vEaU3EbMnAZ0mnWp1FPKjVqCvU+hUZvhCZfhCZvmBZfmDZ/GBXwBCJ4AhGqXH6o7vpydrOeDLuZjCQzmQ4z49OsTMqwMy3bQVaStX2qS9c+1aVHL9NbQoghYDXq2Xn7GcN27XiZTCaOP/541q1bx9y5czGZTMyZMweAFStWsG7dOgoKCjCZTCxfvhwAl8vFsccey9VXX93pXN/+9reZNWtWj9c6/OdecnIyzc3Ncd/v0WjAg5idO3dSVFTEueeeC8CZZ57Jiy++yC233IJOp8Pv9/Pyyy/zla98ZaAv3T1FAaNVe9jS+n++kL89oHGBrwW8DeCpB0+D9vA2gLcJfE3a+36X9oiGUCJBFE8dOk8dRsAOdJnPrujAmgb2TEjOgnGZYM8kYs+izZxNjZJNbcROrRfK21QqWqNUu8PUuCPUu0N4gpH2QMdLaZP3iNMnWwyMT7UxPtXK+FQb+alWCjJsZDjMWI16rCY9NpMBm0mPJYF/9EII0RtFUeKe0hluJ598Mr/97W855phjWLFiRWxEuyMvpqCgIDZiAzBx4kSam5tZvXp1QtcJBAJUVlaSl5cHaNVvS0tLKSwsBLRBgUMTfUHLrTnaJfS3qK6urlNi0ebNm7FYLEycOJFjjz0WgOeee477778/lvNy6623snjxYi644ALOOussnn32WQwGQ6fpolHHaNEejq6HBbsUDoG/BVqroLUS2qqgrRraarUAyNvYHvg0a/up0fZgqAHqd8VOowdS2h8zjHZwZEJSLiTnQk4uanIuIUcuTREHJW0KZW6FSrdCuVulrDVKeWuYOneYVn+YndWt7Kxu7XSbKVYjEzPsTMtOYlq2g0mZDuxmLZixtufw2IwGrCY9JoPk5QghxrZVq1bxox/9iEcffZQbbrghtn3FihVce+217Nu3r9OoyzXXXMOZZ57JU089xWWXXQZoAcrDDz/c6fiu3HHHHbFFLx2rl84++2wAZs2aRVVVFTt27GDu3Ln4fD4efPDBAfyko1NCQUxtbS3PPvssABdccAGVlZU8++yznHrqqbEgZtasWZx33nmxYwoLC9m2bRsPPfQQ69ev58QTT+TZZ589+pqAGYxawOHIhNx5R74fCWtTXuGANq3VVg3NpeAq0wIfd337aE+d9jXQBiEPNHug+UDsNApgAnKMdnKSc1mWnAvJeZCfRyBpPF5rDo0BA3sbI+xtgVK3gdI2KGuN0uCN0OIL8VF5Cx+VtwCgU6Aw/WBQMy07iXSHNuRpNupIthhIthhJshixmmTERggxtixatIikpCT279/PypUrY9unTJlCeno6lZWVrFq1KrZ99erVPPjgg3zzm9/kzjvvJC0tjeLiYq644ooer2O326msrGTatGnYbDb27t3LX/7yl9jCmMWLF3PFFVewfPlyFixYwIEDB1i2bNlgfORRRVFHQenc1tZWnE4nLperbzkxo100ouXxhHwHv3oboaVcC3Zaqw55VGpBDt38seoM4BwPqROJphYSSJqAxz4Bry6JBm+Y/S2wpwV2NivsaoIm/5HnSbObmJbtYOa4ZBbkp5KZpAU1JoNCksXYHtQYsJtHx3CxEGJw+f1+SkpKmDhxYo+LOkaq999/n9bWVk4//fROCyQ2b95MY2Mjp5xyCkajsdMxXq+X7du3oygKc+bMweFwxN5raGhg27ZtnH766YA2QPDpp59yyimnUFJSQmlpKXPmzCEjI+OIe9m7dy8NDQ3MmjWLaDTKhx9+yBlnnIGiKEec9/DXg6WnP9/B/vktQcxo1pGw3JGkHPRoj0AbtFYfnLZqrQJXhTayE/J0fS5LCqROJJJSiC95Iq1JU3ApTmq8Crua0R5NsL8Voof9jclPs7EgP4WFE1KZmuWI1b8x6hUcFgMpVhNpdpNMPwlxlBrtQYzo2XAGMfKr8mh2aMKy/ZCIPRLqHNQE3droTTSqjdI0l0BTiTYN1VyiBTn+Fqj+CH31RzgABzDOmsbUjBksSplK28QpuGZPwBsxsrcFdjbBtnotsClv8lLe5OWl7VUkmfXMzU9l4YQU5o5PIRRRafaEKGnwkGQxkGbXAhpJFhZCCNFfEsSMRXojWFO1R4doVFtRFWiDrJna80h7A7OQH1pKDwY3DXugaT+KrwlD+UaSyjeSBIzTGYmkTmZa6lROSp1CW+E0WnCytQ4+rIOtddAWiPBeUQPvFTWgU2B6lo1jCzNYPkULstr8YUobvTjMBlLtRtLtZsmlEUII0ScSxBwtdDqwpmiPDiFfew2cNrCnQdYMbYoKtMCmcR/U79YedbtQAq0YGneT1LibpI5TOCcyNe0Yziqci3vuNHa6DHxYC5trocwNu2q97Kot45kPy1iQa+Ok6TksKMzAHQjjDoQpb/JhM+lJs5u05d0S0AghhIiTBDFHs46pqKRs7XUkfLD+jb9FW0aec4z2nqpqScT1u6But/a1+QBGVwnprhLSS14iqrdQmDaLk9Ln4i44hjJy+LBO4e1K2N0MWyu9bK3cT5JpPycUODhpRjYF2RntlYt9VDT7cFqNZCdr7RWkwrAQQoieSGKv6F4k1F6sr0WrXxMOdH7f1wxVH0PVNqj6SNvvEEFrJu70ubRlzGe3aS5vVhlZWwGNh9RrKkyGlRMdLJ+aicOZATptJMZk0JGdbCYrySIJwUKMcpLYO7bJ6qReSBAzQgS9BwMav0tb+t1BjWr5NFXbtEfdLoiGY29HDDZaMxfRnL2M96KzeaPCwPs1EGrvIK9TYHE2nDPTydTx44haUgAtdzndbiLbaSHZ0nkJoxBidJAgZmyTIKYXEsSMQKqqBTSeRq1mTeSwJpYhH9TsgMptUPa+1oKhXdjooDVrMbUZy3jFM4s3KnXsPqQ9yIxUOH+qgYUTs4jYsogatH8UdrOe7GQLGQ6z9HkSYhSRIGZskyCmFxLEjHCqqo3MeNsDmsOnndQo1O6EA++ilr6Hcsi0U9iUTGvWUvYmH8cTDdN5o0IXG50Z74ALJsNJhXZURxYhSzroDBj1CrkpVnKSLbGaNEKIkUuCmLFNgpheSBAzyvhbDzbBDPk6vxeNQO0nqCXvopZuRBdsi70VMqdRmXMKTwRP5u/lTjzts1HpFjhvEqwuUDA50ghas4iYnZgMOsanWslKMksSsBAjmAQxY5sEMb2QIGYUC7i1fk/u+iOnnKJhqN5OZP/bKGUfoAtr3bajip7mzKW8ZDiN31dNo9GvBSh2A3yuEM6dBMkOOwHbOMKWdMwmPeNTrWQ6JJgRYiQarUHMJ598wttvv811113XaXsoFGLjxo3s2bMHm83GzJkzY/0De7Nx40Z27NhBVlYWZ555Jlartcf9g8Ega9eupbS0lAkTJnDaaadhMBxcWFxeXs4//vGPI477yle+Euu7NNgkiOmFBDFjgKpqU01tNe1dug/7axcJES7ZQHTny5ia98Y2+xwFbLCfxr31x1Pk1f5xmHSwZhJcNBUsZgsB+zhC1kysZiP5qdZYg0ohxMjQ3Q+594sbh/Q+lk1OrPHwH/7wB2644Qb8/oNLKt9++22+8pWvYDQaOemkk4hEIuzatQtVVfnHP/5Bfn5+t+e76qqr+Ne//sVZZ53F9u3b8fl8vP322+Tm5na5/4EDBzj11FOx2WwsXbqUTZs2YbFYeOONN3A6nQC8+eabnHbaaXz729/udOyPfvQjsrKyEvq8fSVtB8TYpyhaawR7hpYz466Ftlqt9xOA3ohhyiqYsopQ3T5Cn72MpeJdrO5STnM/ysmGZ9iZt4LfuE9lrWsczxXB62Vw+Qw/p00oweyuIGjPYa8/G7vVTH6qjVS7aXg/sxBiTNm5cydnnXUW3/jGN/jFL36BTnew/MO2bdvoaUzgP//5D3/961/56KOPmDt3LoFAgKVLl/KDH/yAJ598sstjvve975Gens7GjRvR6/VEIhGOP/54brvtNn7961932vf+++8fkM842kgBDjH0DGZImQD5iyFnDtgzQTn4V9GYNRXbqu8QvuCvuGZ/maA1G33YyzGN/+Mvge/ybsY9rLbtpiUIv9sB33obtteGsLSVk9TwEeGG/eypbOSzKhe+YKSHGxFCiPj9/Oc/Jysri7vuuqtTAAOwcOFCJkyY0O2xzz33HMuWLWPu3LkAmM1mrrzySl588UXC4XCXx3z88ceccsop6PVa/Sy9Xs8pp5zC3/72tyP2ffrpp/nLX/7Cpk2b+vrxRiUJYsTwsqZq7Q7yl0DaRC3AaWeyOXEuuhjO/yMNx91MW8YCVBTy3Tv4Q/R23kq9m+OMRRxogx9/ALd9CBWuCGZPNY6GjwnW7uWTsnrKm7xED2+9LYQQCVq7di2nnHIKJlPio7w7d+5kxowZnbbNmDEDr9dLaWlpl8dMmjSJrVu3dtq2detWampqaGw8OBXndDr597//zZtvvsnq1as544wz8Hq9Cd/jaCTTSWJk0BvBOR6ScrVEYFclhLR/hCaDgYzpywlMXkZNTSnm3f8kteptJvs+4Vn9J3xiW8jNbRfyYW0hW+vgrEL44rQoyWodRn8jdd7xNLrzmJjpwGmVgnlCiL5paGhg/PjxfTq2ra2NlJSUTts6Em/b2tq6OALuuOMOTjnlFM455xyWLVvG+++/T0VFReyY9PR0pk6dyr59+8jMzASgsrKS+fPnc9ttt3Hvvff26V5HEwlixMii00FSDjiytSXarnKtQSVgNugYN34ivuxvU161hqS9L5BS9S7HBLbxH9M2PjAs4SeeC3m5ZDxrK+DSqXDOxAiWtlIivjr2uAtJy8iiIN2GUS+DkEKIxKSkpFBTU9PjPl6vlz/96U+x13PnzuXkk0/GZrMdEay0trYCYLPZujzX0qVL2bt3Ly+88AJVVVVceumlqKrKZZddFkvsLSgo6HRMXl4el1xyCW+88YYEMUIMG0UBe7r28LWAq0JrdwBYjXryCybSnHED+2vWkF70As6a9zku/CGvmzezVreMO30X8OjOcayrgO8sgMJkH/bmXbi9NexoLWRCdjqZSbKKSQgRv+OPP54NGzYQjUaPyInpEIlEOHDgQOx1Xl4eANOmTaO4uLjTvsXFxRiNRgoLC7u9Zm5uLt/61rdir6+99lomT57c4/Jpk8mEy+WK4xONfvLrqBj5rClaAnDuAnBkgaKgoJBmNzGhcCqepTdSvOxeXFlLUFA5JbqRt8zf4wHTH3C3NvHtd+Dv+yASBWOgGXPddsqLd7Kzohl/SBJ/hRDxuemmm9i9eze//e1vj3ivurqauro6kpKSuP/++2OPiy66CIDzzjuPd955h7KyMgCi0ShPPfUUZ555ZizHprS0lPvvv5/mZu0Xtrq6OjweT+wapaWlPPnkk1x//fWxbbt27ep0H21tbfzzn//khBNOGNgPP0LJSIwYPcwOyJyurWxqKgFvI0a9Ql6KFbdtGjWp36GhaT+ZRS+Q3LCN83TvcJplM/eGLubJ3afxQY2OG+fDhKQoZk8VEX8Dn7nyyckrINdpkUJ5QogeLV++nL/+9a984xvf4D//+Q8rV64kHA6za9cutm7dyv/+979ua7N84Qtf4PHHH2fVqlV88YtfZPPmzezcuZP33nsvts+uXbu48cYbWb16NampqTQ2NnLGGWdw1llnEYlEeOyxxzj33HM7BTGPP/44Gzdu5KSTTkJVVf72t79ht9u56667Bv37MRJIsTsxevlaoLlEqwoMqKg0uoPUu4OYm/eRs+cJbK4iAD5VJ/GD4NXsVQq5fDqsmQz69pglbErClD2DybkZmAwyOCnEQButxe42btzIP//5T375y1922t7Y2MhLL73E3r17YxV7zznnHMzmnqeoI5EIf//739m+fTtZWVl86UtfIicnJ/b+7t27+cMf/tCpUF1ZWRnPPvssHo+HVatWsXLlyiPOu2XLFt58802CwSCzZ89mzZo1sWXZQ0Eq9vZCghjRI3cdNJfGCucFI1FqW/20eoOkVqwlu+hZ9GEvERT+Ej6T34QvZEKqhRvna00mAVSdnkjqRAomTCTFJkXyhBhIo7XtgIjPcAYx8munGP0cWZB3LKQWgs6ASa8jP9XGhHQHbQWnUXT8L3FlH4cela8aXuEt8/fJc23l+rfhn8UQVUGJRjA0FnFg90eU1bf2WHlTCCHEyCBBjBgbdDpIyYfxiyA5FxSFJIuBSZkOLM5MKuZ+i9IFNxG0ZDJOaeTPpl/xW/1v+M/ORn64EWra60IZfQ24ij9k14EKSfoVQogRToIYMbbojZA+GXIXgjUVo16hIM1GjtOMJ3M+Rcf/gvrCc1AVPav1m3nL/H0Wu17lxneibK3TTqGLBFCrP2HP7k9pbPP3fD0hhBDDRoIYMTaZbNqy7IxpoDOQbjczMcOOyWKlbuqlFC/9OV7nVOyKn58Zn+CP3MXvNrXw7F5teglVxdBaQeWuTZRU10vbAiGEGIEkiBFjW1I25C0EWxpWo55JGXZS7UYCSRMoWfxTqmZcSURvYZl+J6+Yb6Zo32f8fDN4Qtrh+pAbT8kWdhftk2aSQggxwkgQI8Y+gxmyZ0PGNHR6I7lOK/lpNgx6Pc35p7J/6Z34HflkKi6eMt7FsQ3/4jvvRClrrxCuqFHU+n3s+2wrTW6ZXhJCiJFCghhx9DhkVCbZYmBiph272UDQnsv+JbfTnLsCvaLyfeNz3Bn6Jbe928qGqoOH6731VOzaRGXj0VHOWwghRjoJYsTR5ZBRGZPRRGG6jaxkMxjMVM3+GpWzriGiM7FSv50XDT/itW17+ctOrWUBgD7opmnfZoor6yRPRgghhpkEMeLo1DEqY00l02EmP9WGTqfQkreSkiW3E7CNI1dp4u+mO0g58F9+8oGKK6AdqosE8JVtY/eBMkId0Y0QQoghJ0GMOHoZzO0rmKaSZNVGZYwGhUDSBPYvvRNX9nEYlQg/MT7Nta2/5pZ3POxr0Q5VohHUms/YtW8f3mB4WD+GEGLwRSIRKioqqK6uHnPFMBsbGyktLR3u2+gTCWKESMqBcXOxWm1MTHdgM+mJGqxUHHM9VTP+j4hi4HT9Vh6P/ojHNx7go/r241QVfWMR+3Z/SrMnOKwfQQgxOJqamrjmmmtIS0tj0aJFLFmyhKSkJC688EK2bds23Lc3IB5++GEuueSS4b6NPpEu1kIAmJNg3HyM9bso0EFViw+XL0xz/mn4kiczfscDTPDX86RyO1//8Lu0zp/NijztUENbJWV7fPgmziU31T68n0MIMWBaWlpYtmwZWVlZbNq0iRkzZgAQCAR49dVXWbt2LQsXLux0jNfrpaamhpycHGw2W6f3Ghoa8Hq9TJgwAa/XS0NDA/n5+SiKEns/EomQnZ3d43GNjY3k5+d3ec+hUIjKykqysrJ6vH5TUxONjY1kZmbS0NCA3+9n9+7dAIwbN47q6uojzq3T6Zg2bVoC38HBJyMxQnQwmCBnLrrkXMan2shM0jrS+p2T2H/cz2lLmYVD8fNn473s+PhDXtp/yKG+Jhr2baG4pmnMDTULcbS6++67qaqq4vnnn48FMABms5nzzjuP733ve7FtkUiEG2+8kbS0NE444QTS0tK48cYbiUQO1pe6//77OeecczjttNOYPHkys2bN4phjjmHr1q2ceuqpzJkzh0mTJnHyySfj9Xo7HXfeeedxzjnnMGHCBKZPn868efPYv/+Q/4SAX/3qV2RnZ3PiiSeSnZ3N5z//eRobG484z5o1a5g2bRoXXnghr732Gk8//TR79uxhzZo1rFmzhjfeeCP2vONxwgknMHfu3MH4NveLBDFCHEpRIGMKZEwlK9nK+FQrigJRo4PyhTfhylyMWQnzkPEBWne/xRO7oSNm0YfceA9sY1d5LWFJ+BWie6oKQc/wPBL4JeOf//wn55xzDjk5Ob3u+6c//YnHH3+cTZs2UVVVxQcffMBf//pXHnnkkU777dixgy984QtUV1dTVVVFIBDguOOO4//+7/+oqamhvLycvXv3HnHcxx9/zLx586ivr6epqYnc3FyuvPLK2Pt//etfeeCBB9i4cSPl5eXU1NQQCoW47rrrjjjPSSedRENDA9u3b+eSSy7h29/+NvPmzWP37t3s3r2bCy+8MPZ89+7drF+/HovFwhVXXBH3926oyHSSEF1JygGjDWfdLox6hfJmH2FMVMz7NpFdfyatch13G//ML/e38Tv/eVw7V0Gv01YuRSs/Znd4DtMm5GIyyO8JQhwh5IW7cofn2j+qAlN8076lpaVcfPHFnbZVVFTgdrsBMBgMTJkyBdDySr72ta8xb948AObPn88111zDww8/zNe//vXY8VOnTuWqq64CIDk5mVNPPZX33nuPL33pSwCkpaWxatUqPvroo07XTU1N5dZbb0VRFCwWC/feey/z5s1j7969TJs2jd/85jd88YtfxGAwsG/fPlRV5ZJLLuGaa64hGo2i02n/F+Xl5XHjjTfG/e0KBoOcf/75TJo0id/97ndxHzdUJIgRojuWZMidj61uFxP1OsqbvPhDUD3zaiKmZDJL/s33jc/xl+pW7glexveO1WHWt69cqv2UXarK9Am5WIz64f4kQog+MBqNnaZ1AH7+85+zbt06mpqaiEajNDQ0AFBUVMQxxxzTad+5c+ce8YM/N7dz8Ga327vc1nHeDlOmTMFkMsVez5o1C0VRKCoqYtq0aezcuZPGxkb+9a9/dTquoKCA5uZm0tPTAZg4cWIsByce3/jGN6isrGTLli0Yjca4jxsqEsQI0RODGXLmYmoqplBXTVmjF28wQt2USwgbkxm390muNLxKSpObn71/DbcsNeAwaoGMvu4zdqkwLX8cdrP8UxMixmjTRkSG69pxOuaYY/j00087bXv44YcBuO+++7jnnnti261WKz6fr9O+Pp/viOTavvL7O7c8CQQCqKqK1WoFtIDrRz/6Ed/85jd7PI9eH/8vVb/97W959tlnee+998jMzEz8poeAjHUL0RudDjKmos+YwoR0Gzaz9p9AU8GZVMz5JlH0nK/fwA3uX3PrhgCN7f/XKNEIxvqd7C6rptUfGsYPIMQIoyjalM5wPBIYhbj66qt544032LhxY6/7LliwgHXr1nXa9tZbb7FgwYKEvz1d2bNnDzU1NbHX69atw2g0Mnv2bACWL1/OCy+8cMRxoVDv//dYLJYj9lu7di3f+973eOyxx5g/f37/bn4Qya+HQsQrORe9omcCeylv8uIJhHGNO4GIwc74HQ9wMh+THLybm9/9Ht9f5iDPAUo0jLlhJ3tVmJyXTard1Pt1hBAjwpVXXsmGDRs47bTTuOmmmzjhhBPIyMigoqKC//73v7FREICf/exnnHzyycyYMYPTTz+d1157jX/84x9HBDZ9FY1Gufjii7nzzjtpaWnhuuuu45prriErKwvQVlKtWLGCyy+/nKuvvhqdTseGDRvYuHEjL7/8co/nnjFjBnv27GHt2rXk5uZiMpm4+OKLufjii5k7d25s6fVIXGItQYwQiUjKRq/omMBuytoDGXfmAkqP/RH5237BIvby++jtXL/xh3zv+LRYIGNp3Mk+FSbmZsWWbgshRjZFUXjssce44IILeOaZZ3jppZfQ6XQUFBRw7rnndloddMIJJ/C///2P3/zmN/z973+noKCAV199leXLl8f2yczMZMKECZ2ukZWVdUTNl5ycHMzmzv9PLF++nMsvv5w777yThoYGvvzlL3PrrbfG3l+0aBFbtmzhvvvu48YbbyQpKYmTTjqJJ554osfrA5x11lnccMMN/OhHP8LlcnHdddeRkZHBli1bWLNmTWw/s9nM9u3bE/smDjJFHQVFLVpbW3E6nbhcLpKTk4f7doQATwPR2l2UN3tx+7W2A2Z3Oflb78EcbKYkms01up/xw+VOctqnxFWdAU/qLArGZZLjtAzfvQsxxPx+PyUlJUycOBGLRf7uJ+rHP/4xGzZsYP369cN9K13q6c93sH9+S06MEH1hz0CXPYv8NDtJFm1AM+DIp3TJz/CbM5moq+VX0V9w50Yf9e25fko0jL15J6XV9ZQ3eXs4uRBCiHhIECNEX9nT0WXPYvwhgUzImkn5sT8kaExmrq6EO8K/5qcbQzTFkn3D2Jp3UVXfREmDZxhvXggxWnQ3DSQkiBGif2xp6HJmk59uJ9mqBTJB+zjKF9xEWGfhBP1nfC/0MD/eGMUV0A7RRUPYmnZS19hMcb17GG9eCDEafPvb3+6U2yIOkiBGiP6ypqJkz2F8mgNneyDjd06iYv6NRBU95+g/4P8CT3DL+ypt7c2uOwKZ+uZWiuvd0m9JCCH6QIIYIQaCNQUlZw55aQ5SbFog40k/hso530BF4QrD65zl/Tc/+QA87eUYdNEQ9ubd1LV4KK73SCAjhBAJkiBGiIFicaLkHENuWlIsR6Y153hqpn8ZgO8bn+NY91p+ugl82oImdGEftpa91Lf6JZARQogESRAjxECyJKPkHENemgNre8+kpglnUD/xPAB+bvgz+a7N3PYhtK/MxhBsxdK6n/q2gAQyQgiRAAlihBho5iT0WTPIT7NhNGglzusmX0xz7kr0isrvjL/H2rSLn2+BYEQ7xOSrx+yukEBGCCESIEGMEIPBno4xYyIT0mzodQooClUzr6I181jMSohHTb/C21DGPVshEtUOMbsrMPrqJZARQog4SRAjxGBJyceSksv4NKvWc06np+KY6/GkTCdZ8fK46R4q6+r402fQEa9YWkvQB1vbAxlZtSSEED2R3klCDKb0KTjCPnIjKpXNPlS9ibL532XiljvIdpfzuPFezjtwBzk2G5+fDIoaxdayB0/aHOrbANxMznSgJNB5V4hRo+Tdob3exBMT2v3pp5/m9ttvZ8+ePZ22P//88zzzzDPs2bMHm83GzJkz+frXv96pT1JX9u3bxy233MKOHTvIysri+uuv56KLLurxmBdffJE//elPlJaWMmHCBG666SZOOeWU2PvvvfceX/rSl4447q233mLy5MkJfNrRSUZihBhMOh1kzSIlKYmsZK2hW9TooHTBDwmZ05isq+Ye46P8eafKxmrtECUawda8ByUaor4tKCMyQgyTtrY2SktLY69VVeWyyy7jm9/8JieffDJ///vfefTRRznrrLP40Y9+xKZNm7o9V1NTEyeddBIGg4Gnn36aL3zhC3zxi1/kxRdf7PaYxx57jEsvvZQLLriAF154gQsvvJBzzz2Xd955J7aPz+ejtLSU9evXd3ocLRV+ZSRGiMGmN0L2LDKj2wlGorR4QoQtqZTP/RYTt9zB2foP2Bydzn3bzuDu42F6KugifmzNe/CkzaK+LQh4mJLlGO5PIsRR7dFHH+WZZ57hgw8+YMmSJbHt8+fP59JLLyUSiXR77B//+EeCwSCPP/44RqORY489lq1bt3L77bdz/vnnd3nMww8/zOWXX85Xv/pVAGbPns2mTZu47bbbeOuttzrtW1hY2P8POArJSIwQQ8Fkh8wZ5DqtONpryPhSplEz9YsA/Nj4FDPVIm7/EGrae0PqQ26sriIA6tsCHJBeS0IMqz//+c+cfvrpnQKYQ+n1+m6PXb9+PSeffDJGozG2bfXq1Wzfvp2WlpYuj3G73aSkpHTalpKSwrvvvksoFOq0ffHixcydO5dLL72U7du3x/eBxgAJYoQYKrY0lLRJjE+1YjFq//SaJqzGlbUEIxH+aP4tBNv42SZi7QmM/ibMbdpwdrXLL92vhRhGn376KbNnz+7TsRUVFeTk5HTa1vG6srKyy2NWr17N008/zb59+wAtp+bpp58mFApRV1cHgE6n4xvf+AYPPPAAjz76KA6Hg0WLFrF58+Y+3edoI0GMEEPJmYc+eRz5aTYM+val17OuIWDNJpsGfm95iAp3lLu2QKhj6bWnGqOvHoCKZh9VLb5h/ABCHL3C4TAWi6VPx0ajUQyGzhkcHaMy3U1D3XHHHaxevZpjjjmGrKwsVq5cyZe/rFUA78iTW7lyJQ899BDHH388S5Ys4ZFHHuH444/n9ttv79N9jjYJBzH//ve/WbJkCdnZ2ZxwwgmsX7++x/1DoRC33XYbc+fOJScnh4ULF/KrX/1KEhXF0St9CiZ7KgXpNnQ6hajRRvm8G4jqjJzAdr5lfIkdjfC77Z2XXutC2ihMaaOXulb/MH4AIY5OhYWFFBcX97hPfX09hYWFscctt9wCQEZGBo2NjZ32bWhoACAzM7PLc9ntdh577DE8Hg87duygrKyMvLw8TCYTWVlZgDYSc7jjjz+enTt3Jvz5RqOEEnvffvttLrzwQn71q19xzjnn8Nhjj7F69Wq2bNnCnDlzujzmxz/+MY899hh/+9vfmD17Nh988AGXXXYZOp2OG2+8cUA+hBCjiqJA1kwske2Mc0apbPYRSCqgesb/kbfzT3xb/zybI1N5q2I24+xw6bT2pdeufbjTZoPOwP4GDzqdQobDPNyfRoijxgUXXMD9999PVVUVubm5Xe6TlpbW6Zf7pKQkAJYsWcK///3vTvtu2LCBvLw8xo0b1+N19Xp9bOrpxRdf5NRTT8VkMnW7f0lJCWlpafF8pFEvoZGYe++9lzPOOINvfetbTJw4kdtvv52ZM2fy61//uttj3nvvPc4991xOOeUUcnJyWLNmDStWrGDjxo39vnkhRi29UVt6bbeQYteGlFvyVtKcexI6VB6x/p5MmnlqD6yt0A7RhX1YW/cD2ghNUZ2bZk9wuD6BEEedm266iQkTJnDBBRfw6aefxrZXV1fzs5/9jI8//hi9Xt9pJCY9PR2Aq6++mrKyMh588EFUVeWzzz7j0Ucf5Rvf+EbsPG+//Xan0Z5t27bx7LPPEolECIfD3HPPPWzatImf//znsWPuvfdetm7diqqqqKrK008/zfPPP8+VV145RN+V4ZVQELNhw4ZORXYATjvtNDZs2NDtMeeeey5r166lpKQEgB07dvDhhx9y3nnn9eF2hRhDTDZInci4ZEss0bd6xv/hd+Rjj7h4Nvl36InwwMfwiTbqjNHfhMmjFZRRVdhb24bLF+rmAkKIgeR0Otm4cSPz58/nxBNPJCMjg9TUVBYvXozBYGDq1KndHjtz5kyeffZZ7rzzTlJSUli4cCGXXHIJP/zhD2P7dNR86Vh5NHXqVP73v/+RmZlJamoqzz77LK+//jrz58+PHXPSSSdxww03kJqaSlJSEj/84Q958MEHOwVHY5mixpmc0traitPp5Omnn+aLX/xibPt9993Hz372M9xud7fHfuc73+E3v/kNFouFYDDI3XffzU033dTt/oFAgEAg0Ona+fn5uFwukpOT47ldIUaPmk8JuhvZX+8hElUxeaqZtOnH6CM+XrKcy7davoDDCPefCOPsgKLgSZ1JxKT9W9DrFGaOSyLJYuz5OkIME7/fT0lJCRMnTuxzYuxwaGtro6mpiYKCgiPei0aj1NfXY7PZYlNG8YhGo9TW1pKSkoLVau30ns/no7a2lry8vE5LsT0eD8FgkNTU1G7P6/F4iEQiw/Izsqc/347YYbB+fiec2Ht4EpHBYOgxSffOO+/kb3/7G6+++irFxcU8//zz3HXXXfzlL3/p9pi7774bp9MZe+Tn5yd6m0KMHhlTMRlNjHNq//iD9nFUzdaKW53rf4nLk7bhDsFdWyAQAVQVq6sIJar9thaJquypacMbDA/XJxBiTEpKSuoygAHtZ2F2dnZCAUzHcePGjTsigAGwWq0UFhZ2CmBAS/DtKYDp2Odo/CU/7iAmKSkJq9Uay6buUFdXF8uSPpyqqtx7771873vf44wzziA3N5fzzz+fr371q9x1113dXuvmm2/G5XLFHuXl5fHephCjj8EM6ZNxWo2k2bVkvdbs42jMPwOAW9WHmWGqZ38rPLhDm0bSRYJYW4piy5dCEZVd1a34Q91XDBVCiLEm7iBGURQWL17Mu+92btj19ttvs3Tp0i6PUVWVaDR6xPCSxWIhGo12ey2z2UxycnKnhxBjmiML7BlkJ5uxmbSqn7XTvoTXOQVj2MPT9gcwEeatCni1vZWLIejC7K6InSIY1gKZYLj7f1tCCDGWJDSd9O1vf5t//etfvPTSS4TDYR577DE2bdrE9ddfH9vnnnvuiQ2/6XQ6zjrrLB544IFYF9CPPvqIRx55hM997nMD+DGEGAPSp6AzmMlLtaLXKag6AxXHfIuwwU66bz8Pj/sPAH/4FHY3a4eYPZUYAi2xU/hDUfbUtBGJSh0mIcTYl1AQc/7553Pfffdx1VVXYbFY+MlPfsLjjz/eqf243+/H5XLFXv/xj39k+fLlLF26FLPZzGmnncZFF13EvffeO3CfQoixQG/UCuHpdeSlaPPlIWsGNTO+AsCqlhe5ILOSsAp3bwFXe+671VWEEjmYCO8OhNlX1yYFJYUQY17cq5MO5/P5ukxM8vv9BAIBnE7nEe8FAgHM5sSLcw12drMQI0r9XnDXUtPqp9EdBFVlwke/IKlxO23J0zjTfSsVXh3zMuCOpaDXQcTowJM2C5SDv5dkJpml87UYETpWrxQWFnb5c0OMbj6fjwMHDoyO1UkduvuLaLFYugxggD4FMEIcddImgcF8MD9GUaieeRURvYWk1r38ccKbWPSwvQGe0GZp0YfcWNobRXaobwtIw0gxInSstvF65e/jWNTx53r4qqqhkFDbASHEENAbIGMqSs2n5KVaKan3ELJmUDf1C4zb/VdmlD/LLTMX8JNPM3mhCKalwPJxYPLWEjYlE7akx05V0ezDZNCRnTx6anOIsUev15OSkhLrvGyz2VAUZZjvSvSXqqp4vV7q6upISUlBr9cP+T1IECPESGRNheRcTK1V5KVaKWvy0jT+VJJrNmJv2cua5r+wdeJN/KtE4TcfQ0ESjHeAtbUEtzEJVX+wr0pJgwejXhdbvi3EcOjo/dMRyIixIyUlJfbnO9T6nBMzlCQnRhyVohGo+ghCPmpcfho9QUyeKiZ/cDO6aIiyWd/gGwdO5LMmmOCAX58IVgOEzSl4U2d0OpVOgVm5yVLVVwy7SCQSK6svRj+j0djjCMxg//yWIEaIkczfCjU7iESj7K93EwyrZJT8i+yi5wgbHWxZ+Eu++YGTpgCclAs3LdSaZPuSJxGydS5CadQrzM51YjUN/ZCvEOLoNGITe4UQQ8CSDMl56BWF3BQrigINBWfjcxRgCLmZeeBxbl4EegXeqYJ/l7Qf1nag07JraK/qWyPF8IQQY4cEMUKMdCkFYLRhNxlIs5lAZ6Bq9jWoig5n7QcsDW/h6tnarn/eqRXCU9QoVldxrC1Bh0Aoyu6aVimGJ4QYEySIEWKk0+kgtRCArGQzJoOCP3kijQVa1etxux/jvPFeTsqFqAq/2AbeEBiCrZi81UeczhOIsLdWiuEJIUY/CWKEGA3s6WBxojtkWqlu0gUErNkYA83k7HuG6+ZCthVqvfDQJ9phZncFutCRtTlavCEONErNDiHE6CZBjBCjRdokUBRtWsluQtWbqJp1jfZW5Voy2z7jewu1f9TrKmFdRfu0UmvREdNKADUuPzUu/xB/CCGEGDgSxAgxWpgdYNdWHGUlmTEbdXjTZtI0/hQAcnc+ymxngEuna7s/+AlUe0Af8nbqdn2oA40eWrzBIbl9IYQYaBLECDGapBaATo9OURjntKAoUDv1UkLmNMy+WrKKX+CSKTArFXxhuO8jiETB7K1CH2w74nSqCvvq3HiD4WH4MEII0T8SxAgxmhjMkJwHEJtWihpsVM38PwDSS1/B5inl+wvBbtBWKj2zF1BVrK37tQJ6hwlHVHbXtMnSayHEqCNBjBCjjXM8tLcV6JhWcmceiytrCQoqOXufJsuqcu1cbffn9sGnjaAL+7C4y7o8ZSAUZW9tG1FZei2EGEUkiBFitNHptWkl0FYrHTKtFFUMOJo+xdG4nRV5cGo+RIH7tkFbUGsSqQ+4ujxtmz9Mcb17CD+IEEL0jwQxQoxGjmww2QGwtU8rhWzZNE04A4DsvU9DNMLX5sA4G9T7tURfVQVrazFEu86BaXAHKW+SpddCiNFBghghRiNFgbSJsZcd00r1E88jbHRg8VSSWvU2NgN8f6HWluDdKnizAnSRIJa2rqeVACqafdS3Bbp9XwghRgoJYoQYraypYEsDDk4rqSYH9ZM+D0Bm8fPowj6mp8Jl7cuu//AJVLrB5KtDH2zt9tT76920+qXTsBBiZJMgRojRLHWiNipD+7SSzUTz+NO0Sr5BF+kH/gPABVNgbjr4I/DLbRCKgqW1pMsieKC1L9hb04Y/dORqJiGEGCkkiBFiNDPZwJETe5mZZEZvNFI79VIAMkr/i8HfhF6B7y4AhxH2ueDpPaAP+zB5q7o9dSiisqemTZpFCiFGLAlihBjtUgtAZwBAr1PISjbTlrUYT8o0dNEgWcXPAZBhhW/N0w55oQg+aQCzuxIl0n3+izcYoahOViwJIUYmCWKEGO30Rq12TLtUqwmb2UDttMsASKl6F0vbAQCWj4MzJoAK3L9dqw9jbS3p8fRNHlmxJIQYmSSIEWIsSM4DgyX2MsdpwZ8yBVf2MhRUbcl1e/7L1bMg0wo1XnhsFxgCLRj8jT2evqLZR5NHeiwJIUYWCWKEGAt0OkgtjL20GvWk2UzUTr2kvQDeZzgaPgbAZoRvt08r/ecA7GgAS1tply0JDlUkPZaEECOMBDFCjBWOTK3TdbvMJDOqI5umCasByN73TCxQWZAJZ2pFf7VppUAQi7u8x9NHolqibzgiPZaEECODBDFCjCXO/NjTjiTfzgXw1sfev7J9WqnWC3/dBSZfLbpQz0m8/lCUvbVu1G6WZgshxFCSIEaIscSWDkZb7GWq1YTFnkz9pPMByGovgAdgM3SeVtper2pJvr0EKC5fiDJJ9BVCjAASxAgxligKOPM6bcpxWmjOP5WALQdDsJWMAy/H3jt0WumBjyHo82Dy1fZ6maoWv7QmEEIMOwlihBhr7FlgMMdeWo160hy2WAG89NJXOq1GunIWZFmh1qetVjK7y1Eiva9E2l/vxh2QRF8hxPCRIEaIsUang+TcTpsyk8z4chbjSZmhFcArej723qHTSv89ADvqItpqpV5EVdhT00YwLIm+QojhIUGMEGNR0rhYFV/QknyznVZqp30RgJTqdzEfEqjMz4SzDplWCrsbMQRaer1MMBxlb20bUWlNIIQYBhLECDEW6fSQlNNpU4rViJI5HVf2cSioZO7/Z6f3/++QaaW/7ARL6wFQex9lafOHOdDoGci7F0KIuEgQI8RYlZwHSud/4uNSLDRM1lYqOes+xOyuiL1nM8AN87Xnr5TCJ7V+zJ7uG0QeqrY1QF2bf0BuWwgh4iVBjBBjlcEEjuxOmywGPfasSbRmLQYgo+Tfnd6fl9F5Wini6rlB5KFK6j14JNFXCDGEJIgRYixz5mnLrg+RmWSmacrntbdrNmLy1nR6/9Bppcc+UzuN1vQkqsKeWqnoK4QYOhLECDGWGa1gy+i0Sa9TcObNoC1jAQoqGSUvdXr/8GmlXeX1vVby7RAIRSmqj29fIYToLwlihBjrnOOP2JRiM+Kaqo3GpFS/i9FX3+n9Q6eVfr8DlJbel1x3aPaEKJeKvkKIISBBjBBjndkB1pROmxQUnBPm4k6bjaJGyDjwnyMO+7+ZkG6BKg/847M2DP6muC9Z2eKjxdt7wTwhhOgPCWKEOBp0MRqTZDHQNv0iAFKq1mPwN3d632aEr83Rnr9QBLVVpb32VeqgqlBU58YfivTvvoUQogcSxAhxNLCmaiMyh0kpnI8nZTq6aIj00iNHY47PgaXZEFHhwW0BDIclAfckFFHZV+uWQnhCiEEjQYwQRwtn/hGbrCYDvpkXAJBW8Rb6YGun9xUFvnEMWPWwqxnWfVoO0fiXUbsDYUqkEJ4QYpBIECPE0cKWDkbbEZuTJy3FlzwJXTRIeukrR7yfaYUvz9CeP74ziqehLKHL1rUGqGuVQnhCiIEnQYwQRwtF0erGHMZk0BOYqeXGpJW/gb6L5dSfmwjTUsAThr9urkMXTiwoKWnwSMdrIcSAkyBGiKOJPQsM5iM2J009Hr9jAvqIj7Sy1454X6/AdXNBp8CGKtixtyihy0ZV2FvbRkgK4QkhBpAEMUIcTXQ6SM49YrNeryc0+0IA0speRRc+ss7LZCesmaQ9/9M2NyFPS0KXDoSiFNVJITwhxMCRIEaIo03SONAZjtjsmHYSQXsuhrCHtPI3uzz0S9Mg2wr1PvjHh8UJX7rFG6KyxZfwcUII0RUJYoQ42uj0kJR9xGZFZyDcPhqTXvoKSuTIvBeLAb45V3v+clGI0or4+iodqrzJS6s/lPBxQghxOAlihDgaOXK63GybfgohaxaGUCupFeu63GdRFpyUC1HgkY2VRCKJFbRTVdhX65b8GCFEv0kQI8TRyGQDi/PI7To90TnaaExG6X9QIl23DrhmDtiNUOxSefPjfQlfPhiW/BghRP9JECPE0Sqp69EY8/TTCFvSMQaaSal6p8t9Us1w5Uzt+d+2t9DgSjwgafGGqGiWRpFCiL6TIEaIo5UtA/TGI7frjahztCq+GQde6rZC7+kTYHYaBCLw13f3osbZV+lQFc0+XD7JjxFC9I0EMUIcrXQ6sGd2+ZZx+hlEzCmY/A2k1LzX9eHttWMMCmytDvJhcW3Ct6A1imwjGJb8GCFE4iSIEeJoljSu6+0GM8w6D9DqxnTXvXpCElw0VXv+xAdl+IKJd60OhlWK6tx9GskRQhzdJIgR4mhmsoElucu39NNXo+pNWNtKsbXs7vYUF02BHBs0+aL8c1tpn27D5ZP6MUKIxEkQI8TRrrvRGHMS6qRVAF22IojtpoevzdGev/JpXZ+TdSU/RgiRKAlihDjadZfgC+hmngNAct1mjL76bk+xJBuOy4aICn/dsL9PU0OSHyOESJQEMUIc7XpI8CW1EDVnHgoqaRVdtyLocM0cMOngsxo37+9v7NOtSH6MECIREsQIIbqfUgKUmecCkFKxtstWBB2ybXBJe5LvU+8fwBvseml2b1y+EBXNkh8jhOidBDFCiB4TfBm/CNWRgyHsIaW66+XWHc6fDLl2aPaF+ce2yj7fTmWL5McIIXqXcBDjdrv5y1/+wk9/+lOeeuopgsGuy5If7tNPP+W+++7jV7/6Ffv370/4RoUQg6y70RidHmXG2UB7gm8PUz2mQ5J8X/20mvKmviX5avkx0l9JCNGzhIKYxsZGFi1axEMPPYTf7+euu+7ihBNOwOvt+T+qn/zkJxx33HEUFRXR0tLCBRdcwIYNG/p140KIAWbLAJ2h6/emnoZqsGDxVGBv+qzH0yzKgmU5EFXhL++V9Dm/JRiOsr/e06djhRBHh4SCmJ///OdEIhHeeecd7r33Xt59912Ki4t58MEHuz3mpZde4q677uLVV1/lD3/4A3fccQfvv/8++fn5/b55IcQA0unAkdX1eyY7yuRTgPbid724ZjaYdbC7po33ivuW5AvQ5AlS4+o+D0cIcXRLKIh58cUXueiii7DZbACkp6dzzjnn8OKLL3Z7zO9+9ztOO+00TjjhhNg2i8VCQUFBH29ZCDFoumkKCUD7cuukho8wentuMZBlg0umac+f/qC0z0m+AKWNHjyBvh8vhBi74g5iAoEApaWlTJkypdP2KVOmsHfv3m6P+/DDDznppJN4++23ue2223j44YcpLe25qmcgEKC1tbXTQwgxBEz27hN8neMh71gUVNLLX+/1VOdPgjw7tPhCvLC1os+3FFVhX52bSFSWXQshOos7iOnIe0lO7vwfnNPpxOPpet5aVVVaW1t54YUXuOmmmwiHw7z55pvMnDmTV155pdtr3X333TidzthDpp6EGEKOnkZj2pdbV61HF+55GbRRD19vT/J97bMaShv7nt/iC0Y40I/jhRBjU9xBjN1uR1EUXC5Xp+0tLS04HI4uj1EUBbvdjtfr5b333uOOO+7gH//4B5deeinf/e53u73WzTffjMvlij3Ky8vjvU0hRH/ZM7tP8M1dAMnj0Yd9pFS90+upFmbB8nHaaMpj7x3oVxG7utYADe5An48XQow9cQcxJpOJSZMmsWfPnk7b9+zZw8yZM7s9btasWSxcuBCD4eB/ikuXLqW4uLjb/9DMZjPJycmdHkKIIdJTgq+ig5nacuv08tdB7X0J9Fdna/2V9tS28e6+hn7dWkmDB38o8U7ZQoixKaHE3osuuojnnnsulqNSXV3Nf/7zHy666KLYPm+++Sa33XZb7PUXvvAFNm3ahN9/cIXB+vXrmTNnDoqi9Pf+hRCDoacE38mngNGGyVuNo3FHr6fKtMKlHUm+H5b1K0k3HJG2BEKIgxIKYn74wx+SmprK0qVL+drXvsayZctYuHAhX/va12L7bNiwgd/85jex19deey0TJ05k4cKFfOtb3+KUU07h9ddf56GHHhq4TyGEGFg9JfgarTD1NAAyK7rvbn2oNZMgz6Gj1Rfi+X4k+QK0+cOUN0lbAiEEdDPx3TWn08mmTZt4+eWXKSsr45xzzuHMM89Er9fH9jn11FNJSUmJvTabzbz++uu89tpr7Nmzh+XLl3P66aeTmpo6YB9CCDEIHDng72Zl4IxzYOdL2Oq3Y/JUErTn9Xgqow6+MSfKjz+AN3bWcMqMLPLTbH2+tcoWH06rEaet6+7bQoijg6KOgnHZ1tZWnE4nLpdL8mOEGCrRKJRvgmg30z9r74DyTbQWnE75tCviOuWdWw28XxVmdm4yt5w1s19TyiaDwjF5KZgM0gJOiJFqsH9+y79+IUTXekrwhdhy66TKd9CF4lv+fPXMMEadwmdVrXx4oKlftxcMqxTXu/t1DiHE6CZBjBCie/bM7t/LmQspBShhPzl178Z1uhwbfH66GYCnPiglGO5fg8cWb4hql+THCHG0kiBGCNE9SzIYLF2/pyixVgTO0tfiWm4NcMlEP+k2Iw3uIC/vqOr3LZY1eqUtgRBHKQlihBA9c/QwGjNpJZiT0HlqyXRtj+t0Fj1cMdcKwEsfV/W7gF1UhaI6N1FpSyDEUUeCGCFEz3qaUjJYYOoZAKSXvxn3KVdmtjEzx0EwEuXpTT33UouHNxihtMnb7/MIIUYXCWKEED0z2bVHd9qDGH3Nxzgj8SXr6lC5ep4VRYEP9jexs7r/TV5rXH6aPcF+n0cIMXpIECOE6F1PozHJ42DcPEAlq7b3fkodpluaOGW6dt7HNx4YkC7V+xvc/U4WFkKMHhLECCF611MQA7HRGFPJW8Rbf06JRvjSLCN2s56yJi9rd9f28ya1Zdf7G2TZtRBHCwlihBC9M1q6b0MAMGEZmJPB20iOe2fcp82I1HPRwvEAPLelAre//6uMmj0halz+3ncUQox6EsQIIeJj76Hwnd4Ik1cBYC15E5tJ3/2+h9BFQ5xZqCM/1Yo7EOb5reUDcaeUNnrwBmXZtRBjnQQxQoj42DO02jDdmbpa+1rxIVmG+Cr4Ath81XxlWQEAb+yqpbQx/mO7E1VhX60suxZirJMgRggRH70RrD00bk3Jh6xZoEaxl62PfzQm4mdeeoSlE9NQVXji/VIGoqWbNxihTJZdCzGmSRAjhIhfnAm+7HudDHv8HabNniq+tLQAo15hZ3UrH5b0r69Sh2qXnxavLLsWYqySIEYIET9bOuh6GGEpXA5GG7hrSGrZidUY32iMPuQhxxzg3Hm5ADy1qZRAODIQd0xxvZtQRJZdCzEWSRAjhIifTg/WtO7fN1i0VgQAe18jI8kc96nNnkrOmZdLut2k9VXa3v++SiDdroUYyySIEUIkxtHDKiU4OKVU9j7JeLAY4/tvxhBwYVX9XHacluT70vb+91Xq0OwJUdsqy66FGGskiBFCJMaaqiX5did9MqRPgWgYiteS4UhkNKaKpRPTmJGTRCii8syHZQNww5rSRi/+0MBMUQkhRgYJYoQQiVEUsGX0vM8hCb5OiwFznKMxxkATumiQLy8rRAHeL25kT01b/+63XSSqUlTnHpCVT0KIkUGCGCFE4uy9BDETV4DBDK5yqN8V/2iMqmL2VDExw87K6dq01ePvHyA6QIFHmz9MRbNvQM4lhBh+EsQIIRJnTdGClO6YbFB4ovZ872s4rQmMxvjqUaIhLl40HqtRT0mDh3f31ff/nttVtvho84cG7HxCiOEjQYwQom96G43pmFI6sAEl6CHdYYrrtIoaxeSpIcVm4vML8gB49sNyfMGByWdRVSiqcw9I12whxPCSIEYI0Tc99VICyJwBKRMgEoD9b5NiNWIy9NC24BAmXy2oUVbPySE72UyLL8S/t1cOwE1r/KEoBwagvYEQYnhJECOE6BuzQyts1x1F6ZTgq6DEnRujRMMYfQ0Y9TouW6otuX7lk2rqBnCZdF1rgGaPVPMVYjSTIEYI0XeOXtoQTFoFOgM0FUNjESk2IwZ9vKMxNQAcW5DKnNxkQhGVpwdwyTXA/gY3wbBU8xVitJIgRgjRd731UrIkQ8Fy7fneV9tHY+LLjdGHvOiDrSiKwuXLClEU+LCkiZ3Vrf286YOCYZX9DVLNV4jRSoIYIUTfGa1gTup5n44ppZK3IeQj1WbCoItzNMZTDcCENBunzMgG4ImNB4gOYFKuVPMVYvSSIEYI0T+9jcbkHANJ4yDkgwMb0CkKaXGOxhiDLSgRrfXARYvGYzfpKW3ysm5vXX/vuhOp5ivE6CRBjBCif+yZWhJvdzol+L4GQJrNhD6e0RhVxeTVcmOSLUbOXzgegOc2l+MNhvt124eSar5CjE4SxAgh+sdgAouz532mnAKKHup3Q3Mpep1Cqi2+0RiTrw6i2ijJ6bOzyU2x0OoP88+PBm7JNUg1XyFGIwlihBD919uUkjUV8pdoz9tHY9IdJnRxjMYo0QhGfwMABp2Oy9u7XP/v0xqqXQMbdEg1XyFGFwlihBD9Z0vveUoJYOrp2tf96yEaxqBTSLX20A37EB1TSgDz81OZn59CJKry9KaBXXKtqlBc75FqvkKMEhLECCH6T28Ec3LP++Qu1EZkAq1QsRWANIep19gHQB/2oQ8eXFp92XEF6BWFraXN7Kho6ceNH8kXjFDW5B3QcwohBocEMUKIgdFbLyWdHiau1J4XvwWASa/DaYtzNKZ9uTVAXoqV02ZrS66f/KB0wEdOalx+WrxSzVeIkU6CGCHEwLCl977PlFO0rxUfQqANgEyHOa7RmEOXWwNcsHA8DrOBimYfb+2u7csd96i43kM4ItV8hRjJJIgRQgwMg7n3wnephZA2CaJhKHkH0EZjki2G3s9/yHJrAIfZwEWLtCXXz2+pwB0YuCXXAMFwlJIGaRIpxEgmQYwQYuDEMxozuX00pn1KCSA9zsaQhy63BjhlRjbjU624A2H+sa0ioVuNR4M7SIM70PuOQohhIUGMEGLg9JYXAzDxJFB00LAXXOUAWI16kuIYjTl0uTWAXqfElly/8VktlS0DX+elpMFDICzVfIUYiSSIEUIMHKMVTLae97GmQt6x2vPitbHN8TaGPHRKCWDu+BQWTkgloqo89UFpQrcbj3BEZX+9TCsJMRJJECOEGFi2OEZjYlNK62LTQzaTAZtZ3+uh+rAPfcDVadtlx01Ar1P4uLyFj8ubE77l3rR4Q9S4pEmkECONBDFCiIEVz5RS/hIw2cHbALWfxDZnxpsbc9hozDinldWzcwBtyXU4OvCrisqavPiCMq0kxEgiQYwQYmCZ7Nq0Uk/0Jig8SXtedDDB12E2YDX2Phpz+HJrgM8vyCPJYqCqxc+bOwd+yXUkqlJcL00ihRhJJIgRQgy8RGrGlG2E0MEKuRlJcYzGHLbcGsBuNnDxonwAXthWMSg9kNr84UFJHhZC9I0EMUKIgRdPEJMxHZJzIRyA0vdjm5MtBszG3v9rOny5NcDJ07OYkGbDE4jwwtaBX3INUNHsG/CaNEKIvpEgRggx8CzJWvG7nihKlzVjADLiyI05fLk1gE6n8OVl2pLrN3fVUj4IPZBUFYrq3ESlSaQQw06CGCHE4IhnNGbSKu1rzQ5w18U2O60GjIbeexEcPqUEMDvXyeLCVKIqPPVB6aDksEiTSCFGBglihBCDI54gxpEFOXO154fUjFFQ4hqNOby7dYcvLS3AoFPYUenio7KWeO84IdUuPy7fwOfdCCHiJ0GMEGJwWJygj6NDdceU0v612lxNuxSrEYM+jtEYX90R27KTLZw555Al14PUyLG43i1NIoUYRhLECCEGh6LENxpTcLyWP9NaBfW7Y5t1ikKavfcqvgZ/E0r0yBGRNQvycFqN1LT6ee2zgV9yDRAIRTnQKNNKQgwXCWKEEIMnniDGaIUJx2vPD0vwTbOb0Ot6Ho1R1ChGX/2RlzYZuKR9yfWLH1XQOkhTP/VtAZo8wUE5txCiZxLECCEGjyUFdL03dozVjDnwLkQOBgT6OEdjTF0EMQArpmVSkG7DG4zw/NbyeO64T/bXuwmGZVpJiKEmQYwQYvDodGBL632/nLlgz4SgB8o/7PRWmt2E0ktqjK6Lfkra5RW+sqwQgLd211HaODiNHEMRlZIGaRIpxFCTIEYIMbjimVJSdDBppfb8sCklg04h1RbPaMyRCb4AM8cls3RiGqoKj79/YNDaBjR5gtS1SZNIIYaSBDFCiMFlTdWClN50rFKq3Aq+zp2o0x29j8YYA00oka5zU760tACjXmFXdRubSpriues+KW304g9Jk0ghhooEMUKIwaXTa4FMb5zjtVYEahT2v93pLZNeh9Pay3JtVcXo7zo3JjPJzDnzcgF4elPpoOWvhCNak0ghxNCQIEYIMfjsGfHtN/lk7ethU0qgtSLobTTG5O06iAE4d14u6XYTDe4gL++oiu9++qDVF6baJU0ihRgKEsQIIQafNY1eIxCAiSdpq5maS6Bpf6e3zAYdSZaeVzrpIv4uE3y14/V8aekEAF76uIoGdyC+e++DskYvvqBMKwkx2CSIEUIMPr1BW27dG3MS5C/Rnh/ShqBDPK0ITL7uC9sdNymdGTlJBCNRntlU1vv99FG0vUnkYCURCyE0EsQIIYZGPKuUACZ1tCFYD9HOoxlWox67uefRGGOgudsEX0VR+MrxhSgKvL+/kV3VR/ZdGijuQJiKZplWEmIw9SmI8fv9lJeXEwwmVqUyHA5TVFRETc2RnWeFEGOcLT2+KaW8hWBOBn8LVG8/4u2MpF6WW6tqlxV8OxSm2zl5ehagLbmORgdvtKSyxUebX5pECjFYEg5ifvrTn5KWlsb8+fPJyMjgd7/7XdzH3nLLLUydOpXrrrsu0csKIUY7g0kLTnqjN0Lhidrz/euOeNthMmAz6Xs8hclX16mZ5OEuXpSPzaSntNHLuj1d15cZCKoKxfWeQQ2UhDiaJRTEPPnkk/zyl7/kzTffpLGxkaeeeoobb7yRN954o9dj33jjDf7973+zcuXKvt6rEGK0i3dKafIq7WvZRggdOSWT3ktujC4SwBDsOsEXINlq5MJjxwPw9y3luAPh+O6rD3zBCKVN0iRSiMGQUBDz8MMPc8EFF3D88VqztnPPPZcTTzyRhx9+uMfjamtrufLKK3nyySex2+19v1shxOgWbxCTMR2SxkE4AOUfHPF2ssWA2djzf19Gb8+dq0+blU1eipU2f5h/bKuI7776qMblx+WVaSUhBlrcQUw0GmXbtm0sW7as0/YTTjiBLVu2dHucqqp8+ctf5utf/zqLFy/u+50KIUY/owVMcfwioygH2xDsX9/lLr2tVDIGW1Ai3S+jNuh0fHlZAQCvf1ZDRfPgjpYU1bsJR6RJpBADKe4gpq2tjUAgQHp659+kMjIyaGho6Pa4X/ziF3i9Xm6++ea4byoQCNDa2trpIYQYI+JpCAkHg5iqj45oQwDgtBowGnpIFFbVbrtbd5g7PoVFBalEVXjy/dJBXRIdDEelSaQQAyzuIEan03YNhToPiQaDQfT6rpPsPv30U+644w5+9rOfsX//foqKivB6vXg8HoqKigiHu56Hvvvuu3E6nbFHfn5+vLcphBjp4p1SSs6DzPY2BCXvHPG2gkKGvZfRmF4SfAEuO64Ag05hR6WLrWVHBksDqcEdHNQie0IcbeIOYpKSknA6nUcsj66pqWH8+PFdHlNTU0NOTg5f+9rXWL16NatXr2bTpk28++67rF69mtraruesb775ZlwuV+xRXl6ewEcSQoxo5iQw9F60DoBJ7Qm+XaxSAkixGTHoux+N0UWCGIItPV4iO9nCWceMA+CpD0oJDfKUT0mDh0BYqvkKMRASSuxduXIlr732Wqdt//vf/zqtOGpqaqKkpASAU089laKiok6PVatWsXr1aoqKisjLy+vyOmazmeTk5E4PIcQYYo1zSqnwRFD00FgELUf+MqNTFNLtPdeNMfWS4AuwZn4eKTYjta0BXvmkOr5766NwRGV/vUwrCTEQEgpibrnlFt555x1uvfVWNm/ezLXXXkt5eTnf/e53Y/v89re/ZcGCBQN+o0KIMSTevBiLUyt+B92OxqTaTRh03Y/GGIKuHhN8AawmPV9covVV+udHlTR5EivkmagWb4gal39QryHE0SChIGbx4sW8/vrrbN68mSuuuILKykrefvttpkyZEtsnLS2NSZMmdXuO3Nxcxo0b1/c7FkKMfpYU0PVcsC6mo7P1/vVafsxh9IpCak+jMaqKydt7QbvlUzKYmuUgEI7y9KbS+O6tH8qapEmkEP2lqKOgQ1lraytOpxOXyyVTS0KMFXW7wNP9ysaYcACeu0wrerf6Hsiec8QukajK3jp3t5Vxo3oT7owFvbY9KGnwcMs/P0EFfvK5mczKdcbzSfosyWJgdm4ySjztGIQYhQb757c0gBRCDI94VykZzFCwXHveTc0YvU4hzdb9aEw8Cb4AEzPsnDIzG4DHNh4gHB3cJN82f5jKFmkSKURfSRAjhBge1rT4GkLCwVVKB96FSNeVb9Psph5PZ4xjSgngkkX5OMwGKpp9vP5Z70nB/VXR7BvUtgdCjGUSxAghhofeEF9DSNCmkGzpEPRAxeYudzHqFVJ7GI3RKvj2nrDrsBi4tD3J94WtFTR7BzfJV1WhqIepMCFE9ySIEUIMn3inlHT6Q9oQdL1KCSDd0cNojKpi9PdcwbfDyumZTM604wtFeGZTWXz32A++YIQyaRIpRMIkiBFCDJ94l1rDwSmlis0QaOtyF5Neh9Nq7PYUJl8cicRo9Wf+b/lEFGBDUQO7qge/9Um1NIkUImESxAghho/RCiZbfPumFmqPaBhK3+t2twyHudvRGF3Yhz4YX0AyOdPByTOyAC3JNzIE0z1F9e5BrxgsxFgiQYwQYnjFO6UEB0djitd2u4vZoCPZYuj2fWMvTSEPdcliLcm3vMnL6ztrej+gn6RJpBCJkSBGCDG84m1BADBxBaBA3U5o6z6oSHd035vJ6G/URnPikGQx8oXFWgPa57dU0DLISb4Aje4gdW1SzVeIeEgQI4QYXpZk0Pfc/yjGngHj5mrPS97udjerUU9SN6MxihrVApk4rZqexaSM9iTfDwc/yRegtNGLPyTVfIXojQQxQojh15cppf3rtPXJ3chwdB8YmXzx1YwB0OkOJvm+u6+BPTVdJxUPpHBEpajOzSgoqC7EsJIgRggx/BJZpTTheG3kxlUBTcXdn9JkwG7uejRGH/KgC8WfezIly8HK6e1Jvu+VDEmSb5s/TJU0iRSiRxLECCGGXyINIU02yD9Oe95Dgi/0NhoTf4IvwBeW5GM36ylt8vLmrsGv5AtQ3uSVar5C9ECCGCHE8NPpwJoa//6T26eUSt6BaPe5Iw6zAZup6+DI6G/osit2d5ItRi5ZpCX5PrelfEiSfKWarxA9kyBGCDEyJLJKKXeB1rLA3wLVH/e4a0Y3K5WUaBiDvyn+awKnzMhmYoYdbzDCs5vLEzq2r3zBCKVSzVeILkkQI4QYGWwJNITUGWDiSdrzHtoQACRZDFiMXf9Xl+iUkk6ncMXxhQC8vbd+SJJ8AWpc/iEZ+RFitJEgRggxMuiN8TeEhIOrlMreh1DPIxXdjcYYgi6USCD+awLTspNYOS0TgD9v2E84OjQVdoulmq8QR5AgRggxciSySiljGiTnQTgApe/3uGuy1YC5u9EYb/zLrTt8cekEkiwGypt9/GdHdcLH90UwrLK/Xqr5CnEoCWKEECNHIvViFAUmn6w9L36r511Ruh2NMfrre6w305Uki5HLjysA4MVtFdS2Ds1S6CaPVPMV4lASxAghRg6jFYxxNoSEg1NKNZ+Au+cRFWc3ozG6SBBD0JXIXQJwwpQMZucmE4qo/GVDyZAVpjvQINV8heggQYwQYmRJZErJkQU5xwAq7F/f4649jsYkUME3dj5F4aoTJmLUK+yodLGxOP5WBv0Riarsq5VqvkKABDFCiJEmkSklgEkdU0pre50WcloNmAxHroAyBppRoqHErguMc1pZMz8PgCc+KB2ywnTuQJjyJt+QXEuIkUyCGCHEyJJIQ0iAwuWgN0NrBTTs7XFXBYWMpC5GY1QVY4LLrTucMy+XvBQrrb4Qzw5Rg0iAKpcPlzfxwEuIsUSCGCHEyJPIlJLRBgXLtOe9tCEASLEaux6N6WMQY9TruOqEiQC8tbtuyGrHqCoUybJrcZSTIEYIMfL0dUrpwDsQ6Xl0orvRGH3Yhz7YtwBk5rhkVk3Xasc88u5+wkMUWATDUVl2LY5qEsQIIUaeRBpCAoybp7UtCLRB5ZZedx/o0RiALy4pINlioLJl6GrHgLbseqiWeAsx0kgQI4QYeXQ6LZCJe389TFqpPY9jSqm70Rijv6HHhpI9cVgMXL6sEIAXP6qgxjV0gcWBBg/eoHS7FkcfCWKEECNTolNKHYXvKjaDv7XX3bsajVHUKEZ/35dKL5+czpw8p1Y75r2hqx0TlW7X4iglQYwQYmRKpCEkQGohpE2GaFjLjelFd3VjTH2oGRM7p6Jw1XKtdswnQ1g7BsATiFAm3a7FUUaCGCHEyJRoQ0g4pA1B71NKACk2I8bDRmP0ITe6cN9rsOQ4LXx+wXigvXaMf+imeapdfpo90u1aHD0kiBFCjFyJTilNXAGKTqsX4yrvdXcFhcwuRmP6UsH3UOfMHRerHfO3zUNXOwa0btfBsCy7FkcHCWKEECNXokGMNQXyjtWeF6+L65CuRmNMvsSbQh7KoNdx9Yla7Zi1u+vYVd17js5ACUVUiuvdQ3Y9IYaTBDFCiJHLaAGTPbFjJp+ifd2/DtTeRyS6Go1RomEMgabErnuYGTnJnDwjC4A/vbOfQHjomja2eENUtUhbAjH2SRAjhBjZEqneC5C/BIx28NRDzadxHdL1aEz/ppQAvrR0Aml2EzWtfv6+uffprYFU3uQdsl5OQgwXCWKEECNbolNKehNMPFF7HmeCb1crlQzBVpRIILFrH8ZmMvDV9mmlVz+tYXfN0E0rRVXYW9s2ZNWDhRgOEsQIIUY2cxIYumja2JOOVUql70EovqJzqYePxqiqlhvTT/PzU1k5LRMV+OPbQzutFAhF2d8gbQnE2CVBjBBi5Et0NCZzJiSNg7APyt6P65CuRmP604bgUJcdVxCbVnpuiKeVGt3BIa0eLMRQkiBGCDHyJZoXoygJ14wBrYrvoaMxukgAfcCV2LW7YDcfnFb636c1Q9bpukNpo0fyY8SYJEGMEGLks6SAzpDYMZNWaV+rPwZPQ1yH6JQjVyoNRIIvaNNKK9qnlf7wdvGQ1nKJqrBP8mPEGCRBjBBi5FOUxEdjknIgazagQsn6uA9LsRkxGw/+12gINKNEQ4lduxuXH1dAqs2orVbaMrTTSn7JjxFjkAQxQojRwZpgEAMHp5SK1sZdvO7wujGKGsXoi28kpzfatNIkAP73STV7a4d2WknyY8RYI0GMEGJ0sKVpLQUSUXgC6IzgKoOm4rgPc1qNWI362Ov+tiE41IIJwzetBFp+jEfyY8QYIUGMEGJ00OnB4kzsGJMdJhynPS9+K6FDs5IPjsbowz50oYEr5X9Z+7RStcvPc0M8rdRRPyYS7XtbBSFGCglihBCjR6JLreGQNgRvQyT+3BaH2YDNfHA0ZqASfDvOfXX7tNIrwzCt5A9F2S/9lcQYIEGMEGL06EsQk7sArKkQaIWKDxM6NDvp4GiM0d8I0YErVLdwQionTc1oL4I39NNKDe4gta2SHyNGNwlihBCjh8GkVfBNhE4PU07Tnu99NaFDbSYDSRZtabcSjWDsZ1PIw12+rJAUm5Eql5/ntw7ttBLAgQbJjxGjmwQxQojRpS+jMVNP175WfQxtNQkdmpVsRmmvfzeQCb7QPq10gjat9N9hmFaS/kpitJMgRggxuvQliEnKgXELABX2vZHQoRaDHqdVG40xBNvQhX2JX78HxxakcuKUDFQVHlxXhDc4tCMj/lCU4nqpHyNGJwlihBCji8kGRmvix007Q/ta9AZEEwsUMpMsgzYaA/CV4wvJcJioawvw1/cODPj5e9PkCVLR7B3y6wrRXxLECCFGn0Sr9wLkL9XaF/iaoGJzQoea9DpSbEYArfBdnIXz4mU3G7hu1VQUBd4tamBD0cAU10tERbOPFm9wyK8rRH9IECOEGH36MqWkN8KU9uXWe19L+PBMhwWdTkEXDWEINCd+/V5Mz0ni/AXjAfjLhpIhXzmkqlBU58YfGrgVWEIMNglihBCjjzkZ9KbEj+tI8K3cCu7EpoWMeoU0m3bNgawZc6jPL8hjenYSvlCEB9cVEY4ObcJtKKKyt7aNqBTCE6OEBDFCiNGnLw0hAZLzIGcuoGq5MQnKcJjQ6xQMQRdKJJD49Xuh1ylcu2oKNpOefXVuXtxWOeDX6I0nEGF/gxTCE6ODBDFCiNGpL1NKcDDBd98bCRev0+sU0h0mUFVMvvq+Xb8XmUnm2LLrf31Uyc7q1kG5Tk/q26RRpBgdJIgRQoxOlhStkF2iJhyvTUd5G7RppQSl2U0Y9Iq2SmmAE3w7LJucHmsS+eC6ItzDUJDuQKOHVn/8bRqEGA4SxAghRiedTgtkEqU3HuynlGAFXwC9opDpMKOLBAclwbfDFccXkpNsockT5JF396MOUsDUHVWFfbVtQ94OQYhESBAjhBi9+jyl1JHguwU8iS9nTrEZMRoUTN7avl0/DhajnutPnoJep/BhSRPr9gzO9FVPgmEt0XeoAygh4iVBjBBi9LKlEatClwhnPmTPATXapwRfnaKQnWTBEHShCw9e7sikTAeXLMoH4In3D1DZMrDVguPR5g9zoFEK4YmRSYIYIcTopTf2bUoJDknwfb1P3amdViM2kx6TN7FeTIn63NxxzMlNJhCO8vu1+wgNQ5+jGpef+raBX40lRH9JECOEGN3smX07rmA5mBzgqYeqj/p0ihynBZO/XhvRGSQ6ReEbK6fgMBs40Ojl2c1D3+0aYH+9e1gSjIXoScJBzJ49e/jmN7/J5z73OW644QbKy3v+BxUIBHjkkUe47LLLuOSSS7jvvvvweKTZmBBigNjSQenD72N6E0w+WXu+L/EKvgBWo54Uqw6jv7FPx8crzW7i6ysmA/DKJ9V8XN4yqNfrSlSFPTWtBMJS0VeMHAn9y9+7dy9LlizB4/FwxRVXUFJSwpIlS6ip6X449fjjj2fLli2sXr2aNWvW8Mwzz3D88cfj9cocqxBiAOgNYE3p27FT26eUyjeBt6lPp8hyWLD4By/Bt8OxBamcPisb0JZd17cNfR2XYFhlT00bEanoK0YIRU0g7fyyyy5j3759bNq0CYBQKMTUqVO56KKL+OUvf9nlMU1NTaSlHaysWVNTQ25uLs8//zwXXHBBXNdtbW3F6XTicrlITk6O93aFEEcLdx3U7+nbsa98H+p3wYIvw9yL+3SKeneA/capRI2Ovt1DnILhKD97+TNKGjwUptu47dw5mAxDnxWQZjcxLduB0pekanFUGeyf3wn97X/99dc577zzYq+NRiNnn302r7/+erfHHBrAADidTvR6vUwpCSEGjjWtb1NKANNWa1/3vdbn3JZ0u4mk4OD0UzqUyaDjO6dNI8mi5cc8umHo68cANHmClDXJaLoYfnH/q/d6vdTX1zN+/PhO28ePH8+BAwfivuD999+PwWDglFNO6XafQCBAa2trp4cQQnRLbwBrat+OLVwORju4a6F6e59OoVMU8oweiA5+4muGw8y3Tp6KosC7+xp4Y+fgT2V1parFT90Qd9oW4nBxBzHBYBAAq9XaabvNZou915v//ve/3Hrrrfz+978nLy+v2/3uvvtunE5n7JGfnx/vbQohjlb2jL4dZ7DApJXa8z5U8O3gtOhJUwevgu+h5uQ5+eKSCQA88X4pu2uG5xe9/Q0eXD5pTSCGT9xBjMPhwGAw0NTUOfmtsbGR1NTefwN64403uPDCC7nnnnu46qqretz35ptvxuVyxR69rYASQog+r1KCgzVjyj4AX98DkQlGV59q7/XF544Zx7LJ6URUlQfe3EeTJ75fJgeSqsLe2jZ8QVmxJIZH3P/iDQYDc+bMYdu2bZ22b9u2jfnz5/d47Jtvvsl5553H7bffzne/+91er2U2m0lOTu70EEKIHun0WgXfvkibBBnTQY1A8Vt9vgUrQbKNQ1NVV1EUrjlxEvlpNlp8Ie5/cy/hYSiEF46o7K5pHZYifEIk9GvLFVdcwfPPP09RUREAW7Zs4fXXX+eKK66I7fPEE0+wZs2a2Ou1a9dy7rnnctttt/H9739/QG5aCCG6ZOvjlBIcHI3Z2/cEX4Dx+hYM+qEZjrEY9Xz3tGnYTXr21bl5/P0DQ3Ldw/lDUemxJIZFQkHMddddxznnnMO8efM49thjOfHEE7n22mu5+OKDyxL379/P+vXrY68vuOAC9Ho9r732Gqeeemrs8cQTTwzYhxBCCEAbidHp+3Zs4YlgtEJbNVRu633/bhiDLeQ5hm7Zc3ayhWtXTUEB3txVx7o9g79KqiutvjD7G2TVqRhaCdWJ6VBaWkpZWRmTJ08mNze303v79++nvLycFStWALBu3ToikSPnSydNmsSkSZPiup7UiRFCxK1ut9ZKoC82PwI7/w05c+GMu/p8C1FnPtvdTvyhoZtieXFbBc9vrcCgU/jZubOZnDm4NWu6U5BuIzfF2vuO4qgw2D+/+xTEDDUJYoQQcfM0QN2uvh3rroMXv6rlxnzuN5AxtW/n0ZtoSpvPnlp3347vg6iq8us39rK1tJk0u4m7Pn8MTqtxyK7fQVFgapaDdId5yK8tRp4RVexOCCFGPGs/ppQcWTDxJO35Zy/2/R4iQdKUtiENInSKwjdXTmac00KTJ8hv39o3LO0BVBWK6ty4vLL0Wgw+CWKEEGOLTqctt+6rOe3tUErfg9bqvp+ntYrCDNuQLbkGsJkMfPe06ViMOnZWt/LUB6VDd/FDRFXYU9smXa/FoJMgRggx9vRnlVJqIeQdq61Q2vmvvp/H78JGkOxkS9/P0Qd5qVa+sWIKAK9+VsMrn/QjEOuHSFRld3Wr1JARg0qCGCHE2GNNBZ2h78fPbh+NKXoD/K6+n6etmglptiFv0rhkYhqXtlf0ffKDUt4vbhjS63cIRVR2VrcSCEsgIwaHBDFCiLGnv1NKOcdA+lSIBGH3f/p+HncdeqJMyrD3/Rx9dM7ccZw+KxuAh9YXs7OqH8FYPwTDUXZVt0kxPDEoJIgRQoxNfe2lBNoSm47cmN3/gVAfGx1Gw+CuI9VuIsNh6vv99IGiKHxlWSFLCtMIR1V+9cZeyoep87QvGGF3dduwJBqLsU2CGCHE2NTfKaUJy8CRA4E2bVqpr1orQVUpSLdjHKJKvh10OoVrV01henYS3mCEe17dTaM7MKT30MEdCLOnpo2oBDJiAEkQI4QYmxSlf1NKOj3MPl97vvNfEO1jXkfIB54GTAYdE9Jtfb+fPjIZdHzv9OnkpmhLr+99dTeeYVo15PKFKKp3S3sCMWAkiBFCjF39mVICmHIKWJzgroXSDX0/j6sMgKwky7AUoHNYDPxw9UxSbEbKm3386o09w5aj0ugOSnsCMWAkiBFCjF3WVND3I2gwmGHG2drzT1/UKrn1RdAL3iYAJmXa0euGdloJIDPJzA9Wz8Bq1LOruo2H1xcTHaYRkbrWAGWNw5OfI8YWCWKEEGNXf6eUAKZ/Tgtmmoqhenvfz9OijcZYjHrGpw5Pb6HCdDs3njYNvaLw/v5Gnt5UNiz3AVDZ4qOiWQIZ0T8SxAghxrb+TilZkmHK6drzT//R9/ME2sDXAsA4pwWHuR9Jx/1wTJ6Tr63Qmu++8kn1sBXDAyhv8g3biikxNkgQI4QY2ywpoO/n8ubZa0DRQfVH0Fjc9/O4ygFt+fOkTPuQtiQ41IlTM7l0cT4wvMXwACqaJZARfSdBjBBibBuIKSVHNhSeqD3vT2NIXwv4WwGwmw3kOodnWgngnHm5sWJ4D64r5sOSpmG7l4pmH6WNkuwrEidBjBBi7OvvlBIcXG594F1ttVJfuSpiT8enWrEYh+e/4Y5ieMsnpxNRVR54ay/vFzcOy70AVLX4OSCrlkSCJIgRQox9FqeWnNsf6ZNh3IL+N4b0NkJQ+2Gt0ylMynT07776QadT+ObKKZw4JYOoCr9bt4/3ioZvaqna5adEAhmRAAlihBBjn6JoU0L91dGKYN/rsWmhPjlkNMZpNZKV3M8Aqx90OoWvr5jMimmZqCo8uL6Id/fVD9v91Lj8FEtBPBEnCWKEEEeHpBz6nUk7bh6kTYZwAPb8t+/n8dR36sdUkGbDZBimLF+0QOaakyaxanoWqgoPry/m7b11w3Y/da0Bius9EsiIXkkQI4Q4OhjM/U/wVRSY054bs+tlLZjpC1XtNBpj0OuYlDF800oAOkXh6hMncurMLFTgj2/vZ+3u4Qtk6tsCMiIjeiVBjBDi6JGU0/9zFJygTU0FWmHXS30/j7sWwsHYy1S7iXFOS//vrx90isKVyydy+qxsVOCRd/fzxs5+JDH3U31bkH11EsiI7kkQI4Q4elhTwdjPZc06Pcz/kvZ8x9+1RN2+UKPQWtFp04Q027AVweugKApXHF/ImXO0gO8v75Xw2mc1w3Y/je4ge2rbiEj3a9EFCWKEEEeXgRiNmbQSMqdD2A9b/9r387TVQCQUe6nTKUzNdmDQD19+DGiBzOXHFXD23HEA/HXjgWGt7NvsCfFZlYtAuI+dxMWYJUGMEOLo4sjWqu/2h6KDJV/Tnu9fB/W7+3aeaARaqzptshj1TMyw9+/+BoCiKHxxyQTOnZcLaJV9X95e1ctRg8cTiPBpZSueQHjY7kGMPBLECCGOLnrjwBS/y5gGU07Vnm/6ozY91BetVVowc+ipHeZhXXbdQVEUvrA4n88vyAPgmQ/LeOy9kmGb2gmGo3xW1UqLN9j7zuKoIEGMEOLoMxBTSgALv6Ll2DTug6K3+naOaBjajpyqmZhux2bS9/MG+09RFC5elM+lSyYA8PrOWn7x2m68weEZEYlEVXbXtFHb6u99ZzHmSRAjhDj6WJxgsvX/PNZUmHup9nzb4xDsYyNDVyVEO4/kdOTH6HXDmx/T4dx5udx46jRMeh07Klzc+u/Phi2QUFXYX++htFFqyRztJIgRQhydksYNzHlmngPJeeBvgR3P9u0ckSC4j1wBZDMZKEwfgGBrgCyZmMZPz5lFqs1IZYuPn/z7U/bUtA3b/VS1+Nlb65aVS0cxCWKEEEcne5a2XLq/9EZYfLX2fNdL2qhKX7SUdVqp1CEr2UKGw9SPGxxYkzId3LnmGCZm2Gnzh7nzvzuHtU1BkyfIzqpWguE+5iSJUU2CGCHE0UlvAHvmwJxr/GLIW6Tlt2x+pG/niISgqaTLtyZlOoat23VX0uwmbj17FosLUwlHVR5aX8zfN5cTHaapHXcgzKdVrmHL0xHDZ+T8qxBCiKE2UAm+AIu/CjoDVG6Bis19O4e7FvyuIzbrdQpTs5MYIekxgLYU/IZTp8WWYP/r40p++9a+YavlEghF+bSylfq2PraCEKOSBDFCiKOXOQnMA9SzyJkHM8/Vnm9+pMupobg0FmmZq4dxmA1MGEH5MaC1Kbh0yQS+vmISep3CppIm7vjPTpqHaQl0JKpSVOemqE7yZI4WEsQIIY5uA5XgCzD3C2BJ0Wq/9LWvUtDbqTnkocY5raTZR05+TIcV07K45ayZOMwGius9/Phfn7Kz6sgRpaFS3xbgk0qZXjoaSBAjhDi62TO1aaCBYLLBsVdoz3c8C77mvp2npQxCXS9fnpxpx24e/voxh5s5Lpk718whN8VCkyfInf/dxd8+LCMcGZ6EW18wwicVLuqknsyYJkGMEOLoptODI2vgzjf5ZK2ab8gHWx/v2znUKDQVd/mWQa9jek4SJsPI++87O9nCz9ccw6rpmajAS9uruPWlz6hq8Q3L/URVKK73sE8aSI5ZI+9fgRBCDLWBTPA9tK9S8ZtQv6dv5/E2gaehy7fMBj0zxyUNe6PIrliMeq45aTI3njoNh9lASYOHm1/8hDd31Q5bYboGd5AdFS24pe/SmCNBjBBCmOxaFd+BkjldG5EB+PBPfe+r1LT/iL5KHWwmA9NG2IqlQy2ZmMa9F8xlTp6TYCTKnzeU8Ks39tLq62PCcz/5Q1E+q3RR45LppbFEghghhICBTfAFWHgFGKzQsAd2/7dv5wgHoPlAt287rUYmZw3Q6qpBkGY3cfOZM7hsaQEGncLW0mZu+scOPi5vGZb7iapQ0uBhT03bsC0FFwNLghghhACwpWvVdwfsfGmw4DLt+eZHoHp7387TVg0Bd7dvZzjMI27p9aF0isLn5o7jzjVzGJ9qxeULce+ru3l844Fhq7Lb5Amyo0IblZHeS6ObBDFCCAGg04Eje2DPOfNcmLRSm056+x5oO7I/Uq9UVeuS3cMP27wUKzlOS9/vcwgUpNv5+ZpjOGO2ln/06mc13PKvT9hd3Tos9xOOqJQ0ePisqhWP5MqMWhLECCFEh4GeUlIUWHY9pE+FQBusvR1Cfeh0HXBrIzI9KEy3jcgaMocyGXRccXwhP1g9HafVSEWzj9v+s5PfvrWPRvfwVNpt84f5pNJFaaNHVjCNQhLECCFEB6MF7BkDe06DGVbdAtY0rf7Lu7/qW6JvcymEu6+EqygKU7McJFkGqObNIJqfn8ovLpzLqTOzUID39zfynee28+K2imGZYlJVrSP29ooWWoap2rDoGwlihBDiUKkTtWXSA8meoQUyOiOUb4KPnk78HNGwtlqpBzqdwvScpBHVLLI7yRYjV50wibvOP4YZOUkEI1Ge31rBd5//mE37G4clVyUQirKruo29tW3SFXuUGPl/04UQYigZLZCcO/DnzZwOx39Le/7J36Hk7cTP4anX6sf0wKjXMXNcMsYRWEOmK4Xpdm49exbfOnkK6XYTDe4g97+1jzv/u4uypj5MvQ2ARneQ7RUtkvg7CijqKPgTam1txel04nK5SE5OHu7bEUKMdZEwVG6FyCBMLWx9DD79B+hNcOYvIH1KYscbzJAzVwu2euAOhNlV3Uo4MuL/i48JhCO8tL2Kl7dXEYqoKAqcNjObi47NxzFM02QWo47xqTYyHCYUZXQEhiPJYP/8liBGCCG60lqtdZQeaNEIrL0DKrdoy7rPvh+sqYmdw2jVAhlDz4m83mCYXdWjb2qkvs3PU5vK+LBEG3VymA2cOSeH02flDFswYzPpGZ9qJd1hHpbrj1YSxCBBjBBiGKgqVG3TukoPtKAHXvmu1q06czqccU/iNWpMNsiZB/qef6j7QxF217ThC46+4m6fVbl4fOMBypu13ktmg45TZmZz1pycYQsm7GY9+ak2Ukf4SrCRQoIYJIgRQgwTbxPUfjY4526thP9+F4JumHwqLP+2tiQ7EeYkyDlGa2LZg1Akyu7qtlHZOygSVflgfyMvba+K5cjodQonTsngnHm55KZYh+W+kiwG8lNtOG0DWCBxDJIgBglihBDDqPazXpNp+6zqY3jzVm3J9aKrYfaaxM9hTYWsWVqxvh5Eoip7atpwDVPvov5SVZXtFS28tL2KXdVtACjA4sI0zp2fy+TM4Wm/kGw1MD7VhtMqwUxXJIhBghghxDAKerVppcH6r3LXS1qTSEUHJ3wXJq1I/Bz2DMic0etITjSqUlzvpsE9umuh7K1t4+XtVWwpbY5tm52bzLnzcjkmzzksCbgOs4Fsp5kMuxndSO3KOQwkiEGCGCHEMGsshtaqwTm3qsL7v4N9r2uvZ5wNi65KPEcmKQcypsa1a0mDZ0x0c65o9vLy9ireK2ok0v6jLDvZzAlTMjlxagbZyUPfisGgV8h0mMlOtmA19TzNdzSQIAYJYoQQwywSgootWsG5wRCNwMdPwSfPa6/Tp8KKH2iBSSKceZA2Ka5dK5q9lDf5ErzRkanBHeC/n1SzbncdgUNWYs3ISeLEqZkcNykNm2noVzU5rUayk82k2Y/e5dkSxCBBjBBiBHBVQFPJ4F6jYjNs+LXWZ8lohxNugAnLEjtHagGkTIhr19pWPyUNnkGbKRtq/lCEzQeaeHdfA59Wuuj4WEa9wqKCNE6alsExeSnoh3i6x2TQkZVkJivZjNlwdI3OSBCDBDFCiBEgGtVyY0KDPHrhqYe374X63drrWWvg2CtAl8BIQvrkuKsON7oDFNePveaHje4A7xU18M6+BipbDv6ZOa1Glk/JYElhGlOyHEMa0CiKdv00u4lUmwmTYewXzZcgBglihBAjhKcR6nYO/nWiYdj6OOz8p/Y6czqc9ANwZMV/jszpce/vD0UoqnPT5h99S7B7o6oqJQ0e3t3XwHvFDZ0+o82k55g8J/PyU5g3PmVIu4AripYMnO4wkWY3jdkRGglikCBGCDGCVO8Av2torlX2Abz3G604njkJTvgOjF8c//GOLK2hZS+VfUH7YV/R7KOyxTdmppcOF45G+bi8hY3FjXxS4Tqibk5+mo15453MG5/CjJwkDPqhGylxmA2kOUyk2UxjKiFYghgkiBFCjCABN1R/PHhLrg/XVqNNLzXu017PuRAWXN5rgbsYnR6c+ZCc12stGYA2f4jies+orPCbiI7l5tsrWthe4aK4zs2hf6Jmg445eU5m5iQzOcvOxAz7kI2W2Ex6UmxGki1GkiyGIQ2mBpoEMUgQI4QYYer3grt26K4XCcGWv8Dul7XXtgyYehpMPR3smfGdw2iDtIlgS+v9clGV0kYPta2Bftz06NLmD7GjwhULaloPKwqoUyA/1cbkLAeTMx1MzrSTn2ob9JowiqIFNR0BTbLViHEUBTUSxCBBjBBihAkHtQaO0SEerTiwAT54CAKt2mtFB3nHwtTVMH5RfKMztjRtGbax93L9zZ4g+xvcBMMj/sfEgIqqKqWNXj6paGFfnZuiejct3iMrHZsNOiZm2JmS5WBCmo3cFCvjnJZBX85tNelJbg9oHGYDFuPInX4akUHMjh07KC0tZerUqcyYMWPQjukgQYwQYsQJtEFzKfiae993IEWCUPo+7HsVaj45uN2apo3MTD2994ReRaetXkqZEFffpf31Hpo8o7vKb381eYIUtwc0RXVuSho8+EJdB7EpNiO5Tiu5KRbGOa3kpljJdVrIcAxONV+DXsFq1GMz6bGbDVhNeuwmw5AvJe/KiApigsEgF198Me+++y7z589n8+bNXHTRRTz66KPdFvLpyzGHkyBGCDFi+VqgpRT8rUN/bVcl7HsNit48ODqDAnkLYdpqbZRG30NSr97UPsWU0Wu+TF2rn/Jm71E3KtOdaFSlyuWjuN5NUZ2HyhYvVS3+HntTGfUKWUkW0uzaiqT09q+x1w4zdpN+wArjmY067CYDNpMeq0mP2aDDYtQP6XTUiApifvGLX/DLX/6Sjz76iPHjx/PZZ5+xaNEi/vSnP3H55ZcP2DGHkyBGCDHieZug+YC2kmioRUJQ/gHsfRWqtx/cruggaRyktBfASy3Qnifndq47o9Nrq58sTu1hSuoyqFFVFZcvRIM7QJMnNOZqywwETyBMtctPVYuPapePqhY/VS4fNS4/4Ti+X2aDjjS7iRSbkSSzEYfFQJLFgMOsfY1tMxtIshixmfXoEgx69DoFs0GH2ajDbNBjaf9qNugwG3QDmkg8ooKYY445hhUrVvD73/8+tu3888/H7Xbz+uuvD9gxh5MgRggxarjrtZGZwS6K153WKq0PU/Fb3U916QzaaqWO4MaeCSabViXYZNMCGkcWJOWCI7PLoCYSVWnyBGlwB3D5QmN2WfZAiUZV6t0B6tsCNHqCNHmCNHkCNHmCsdd9rdNjMWojLFaj/pCvnbdZjHpMBh0mvYJRr8Nk0Glf9TqMBh1GvaI912vPY8eadFgMBm26ytwR7Ogx6BUMegW9ovQY9Az2z++4s4/C4TC7du3i+uuv77R97ty5/OEPfxiwYwACgQCBwMGseJdLq8nQ2joMw7VCCJEQMyRNBXcduMq1JOChpDhh2kUw9ULwt2j34CqHlkO+RvzgOwC1B+I4nwFMdu1hsGqjNjoD6A2YdQbydEbG6XQEVT3BqI6wqtNGgOCQrtoHv6pdbIvvc8W320jWZQ1lBXBoj2gUgpEogXCUcCRKKBolHFEJd3yNRAlHVUKRKJEosaaXBNofAyB6yOm6q4akAIqitH/VNiid3gMFBUUBf/tS/cFaQxR3EON2u4lEIqSmpnbanp6eTktLy4AdA3D33Xdz2223HbE9Pz8/3tsVQggxYIY4eVmMOY2NjTidzgE/b9xBjNlsBsDr9Xba7na7sVi6bnfel2MAbr75Zr7zne/EXre0tFBQUEBZWdmgfBNGs9bWVvLz8ykvL5eptkPI96V78r3pnnxvuiffm+7J96Z7LpeLCRMmkJbWe32ivog7iLFareTk5FBWVtZpe1lZGZMmdd36vS/HgBb8dARAh3I6nfIXpBvJycnyvemCfF+6J9+b7sn3pnvyvemefG+6p4ujWnSfzpvIzmeeeSYvvvgi0WgUAL/fz8svv8yZZ54Z22fnzp289NJLCR0jhBBCCJGohIKYW2+9lYqKCi644AIeeeQRPve5z2EwGDpN/Tz33HN8+ctfTugYIYQQQohEJRTEFBYWsm3bNmbMmMH69es58cQT2bx5M+np6bF9Zs2axXnnnZfQMb0xm8389Kc/7XKK6Wgn35uuyfele/K96Z58b7on35vuyfeme4P9vRkVvZOEEEIIIQ43elphCiGEEEIcQoIYIYQQQoxKEsQIIYQQYlSKu07McCkuLuaTTz4hKyuL4447btDWmo90xcXF7Nq1i6ysLBYuXIjBcOQfndvtZsOGDUQiEZYvX05KSsrQ3+gwWrt2LXV1dZx33nlYrdZO7zU0NPD+++9jsVg44YQTjnh/rIpGo2zZsoXa2lqWLFlCdnb2EfuUl5ezbds2UlNTOf7447v8uzXWqKrKtm3bqKioIC0tjSVLlnSZeLhnzx527dpFbm4uixcvHrDuwiNJMBjkjTfeIBqNcs4553S5T0tLC++99x56vZ4TTjgBh8PRp31Gm7a2Nt544w3S09NZsWLFEe9Ho1E++eSTWO2z2bNnd3me6upqNm/eTFJSEsuXL8dk6qGz+ChRV1fH+vXrmTJlCgsXLux2P1VVefHFFzEajZx77rlHvN/vn/HqCHbLLbeoNptNPfXUU9W8vDx18eLFalNT03Df1pAqKipSV61apU6ZMkU9++yz1cmTJ6tTp05VP/300077bdiwQU1PT1cXLlyoLl26VHU6neprr702THc99NauXatarVYVUMvLyzu999xzz6l2u11dvny5Onv2bDU3N1fdvn37MN3p0NmzZ486e/ZsdeLEieqaNWvU6dOnq48//ninfe677z7VarWqJ598slpYWKjOmjVLraysHKY7Hhq1tbXqvHnz1Ly8PPW8885TZ82apebk5KibNm3qtN91112nOhwO9bTTTlOzs7PVlStXqm63e5juenD8+Mc/VvPy8tSCggK1oKCgy31ee+011el0qkuXLlUXLFigpqenqxs2bEh4n9HE4/GoX//619WcnBw1JydHPeOMM47Y56233lJnz56tzps3Tz377LPVzMxMdeXKlarL5eq035/+9CfVZrOpK1euVKdOnapOmjRJLS4uHqqPMuCqqqrUSy+9VM3NzVVTUlLUa6+9tsf977nnHtVkMql5eXlHvDcQP+NHbBCzbt06FVDfeecdVVVV1eVyqdOmTVO//vWvD/OdDa2PPvpIXbduXex1KBRSzzrrLHXx4sWdthUUFKhf+9rXYtu++93vqllZWarH4xnK2x0W9fX1akFBgXrXXXcdEcTU1dWpDodD/eUvf6mqqqpGo1H185//vDp//vzhut0h4fV61cmTJ6sXXXSRGgqFVFVVVb/fr7766quxfbZv364qiqL+85//VFVVVX0+n7po0SL1wgsvHI5bHjI33XSTmp+fHwtIotGounr1avWkk06K7fPPf/5TNRgM6rZt21RV1f6OjR8/Xv3hD384LPc8WH7961+rNTU16h133NFlEOPxeNTMzEz1Bz/4QWzbVVddpRYWFsb+XsWzz2jT1NSkPvTQQ2pra6v6pS99qcsg5pVXXlF3797d6ZjCwkL1+uuvj20rLi5WjUaj+thjj6mqqv1fvXLlSvXUU08d9M8wWPbs2aM+9dRTqt/vV5cuXdpjEPPBBx+oBQUF6g033HBEEDNQP+NHbBBz1VVXdfpBraqq+otf/EJNTk5W/7+9ew9p6n/jAP6e07lagjOLLmbZopAaUVCkFUjhpYZYrfojQwvFrtA9u/1Ry7Cwm13MCzUKiQXVLHLlBauRYEVYHTNNqMR1wdnSlc7rPr8/wvNtv82+/n7Y5ubz+u8857Px8HDO53y2Peesp6fHRVkNDjk5OUwkEjGr1coY++dgqK2t5ccYDAYmEAhYQUGBq9J0CqvVyhQKBUtLS2MlJSV2i5i8vDwmFottPkHr9XoGwO7bLE+iVquZUChknz9/7nNMamoqmzRpkk3sypUrzMfHh5nN5r+dosts27aNzZ492ya2adMmNm/ePH5bqVSyyMhImzEHDhxgQUFBTsnR2fpaxGi1WiYQCNiXL1/4WHV1NQPAHj161O8x7qyvRYwj69evZwsXLuS309PTWWBgIOvu7uZjt2/fZgKB4I/nprv40yKmubmZTZ48mRUXF7OjR4/aLWIG6ho/aBtMOI7DjBkzbGJyuRxms9nuv5iGmtLSUkyfPp3/fZ7jOIhEIkydOpUfM378eAQEBIDjOFel6RRnzpxBS0sL9u3b53A/x3EICQmBRCLhY3K5nN/nqfR6PeRyOfz9/VFcXIySkhI0NTXZjOE4jq9FL7lcjq6uLtTW1jozXafavXs3rFYrkpKScPXqVRw6dAgPHjzAqVOn+DF9zT8GgwHNzc1Ozth1OI5DYGAgxowZw8dCQ0Ph4+PDnz/9GTMUdHZ2Qq/X2xw3HMchNDQUQqGQj8nlcjDG8ObNG1ek6TQpKSmIi4tDZGSkw/0DdY0ftB18LS0tdv962fuU36E0ify3/Px83Lp1Czqdjo+1tLRAKpXajR05cqRH1+rFixc4fvw4nj17ZjNJ/M7RceTv7w+hUOjRtWlsbAQAzJkzB+PGjYPFYkFlZSUyMzORlJQE4FdtpkyZYvO6oXCOSaVShIWFobCwEEajEXV1dQgNDUVQUBA/5t/mn6HSNO+oDsCvGvYeI/0ZMxTs2rULJpPJ5gPVUL2O5eXloaamBteuXetzzEDVZtAuYnx9ffHz50+bWO+2WCx2RUoud+/ePSQlJeH8+fOIjo7m445qBfyqlyfXasOGDVi0aBEqKipQUVGBqqoqAMDdu3cRFhaGWbNmOaxNe3s7enp6PLo2YrEYL1++RFFREaKiogAA586dw+bNm7F06VKMHTt2yJ5jO3bswJMnT1BdXQ2JRAKr1Yr4+HjExcWhsrISAM0/vfoztwzV+ed3KpUKarUa9+/fR3BwMB/39fVFS0uLzVhPP47MZjO2b9+OnTt3QqvVAvj1rYvFYoFGo0F4eDiCg4MH7BwbtD8nyWQyu6+U6uvrIRQKMXHiRBdl5TqFhYVYuXIlTp48ic2bN9vsk8lkaG1thclk4mMWiwVGoxGTJ092dqpOExYWBgAoKChAQUEBysvLAQBFRUX819gymQwGg4H/F3UA+PjxIwB4dG1kMhlGjBjBL2AAQKlUorOz06Y2js4xwLNr8/jxYygUCv4nRi8vLyiVSrx8+ZL/BNhXbSQSicPb1D2VTCaD0WhEe3s7H2tqakJbWxt/jPRnjCdLS0vDiRMnUFhYiIULF9rsG4rnGGMMsbGxqKur4+fmt2/fwmKxoKCgAJ8+fQIwgNf4/7WRx1l6GzKNRiMfi4mJceuu7v+XTqdjYrGYZWZmOtxvMpmYWCxm2dnZfCw/P595e3szg8HgrDRdzlFj7+vXrxkAVlpayscOHz7MAgICWEdHhyvSdIonT54wgUBgU4ve+lRXVzPG/mnIfP/+PT9m7dq1Hn/n1uLFi5lCobCJqVQqNmLECL5ZPiMjg0mlUr7B2Wq1svDwcI+9c6uvxt6Ghgbm7e3NNBoNH7tw4QIbNmwY+/79e7/HuLM/NfYeO3aMSSQSmztIf9d700VlZSUf27p1KwsJCeGPNXf2b3cn9XLU2DtQ1/hB+3NSQkICcnNzERkZieTkZDx9+hR6vR56vd7VqTnVs2fPsGLFCixYsACjR4+GRqPh961YsQIikQhSqRQqlQo7d+5EY2MjfHx8kJ6ejtTUVIwfP96F2bueXC5HSkoK4uPjsWfPHphMJmRkZCAnJ8cjHjjVl/nz52PNmjWIjo7Gli1bYLFYcPr0aSQkJCA0NBQAEBcXh8WLF2PJkiXYsmUL3r59ixs3buDBgwcuzv7v2r9/P2JiYpCYmIiIiAjU1NTg3LlzOHLkCN8sv2nTJqjVakRFRSEhIQFlZWV48+YNLl++7OLsB1bvAyI5jkNrays/v8TGxkIikSAoKAh79+7Fhg0b8OHDB3R1dSE9PR0qlYrvC+rPGHd0584dWCwW1NfX48ePH9BoNBAKhVi1ahUAIDc3FwcPHsTGjRvx9etXvnZ+fn5QKBQAgIiICCiVSixfvhzbt29HfX09Ll26hFu3brntgxO7u7tx8+ZNAMC3b99QV1cHjUYDf39/xMTE9Pt9BuoaP6j/xbqtrQ3Z2dl49eoVRo8ejeTkZEybNs3VaTnVw4cPkZOT43BfXl4e/Pz8+G2dToc7d+7AarVi6dKlWL58ubPSHBSqqqqQlpaGrKwsm4YxxhiuX7+OsrIy+Pr6YvXq1YiIiHBdok5itVpx/fp1PHz4EMOHD+cn1N91dnYiLy8Pz58/h1Qqxbp16zBz5kwXZew87969Q35+PhoaGjBq1CgoFAq7J7KazWZkZWXxT+xNSUlBSEiIizL+O1QqFaqrq+3i58+fx6hRo/htrVYLnU4HLy8vLFu2DEuWLLF7TX/GuJONGzfaNZiKRCK+WTU3NxdlZWV2rxszZgzOnj3Lb/f09ECtVqO8vBx+fn5Yu3Yt5s6d+zdT/6s6OjqQmJhoF58wYQIyMjIcvkar1aK0tBQXL160iQ/ENX5QL2IIIYQQQvoyaBt7CSGEEEL+hBYxhBBCCHFLtIghhBBCiFuiRQwhhBBC3BItYgghhBDilmgRQwghhBC3RIsYQgghhLglWsQQQgghxC3RIoYQQgghbokWMYQQQghxS7SIIYQQQohbokUMIYQQQtzSfwCYZKBt8NFjmwAAAABJRU5ErkJggg==", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "import matplotlib.pyplot as plt\n", "import numpy as np\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 }