{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Oxygen comsumption in a lift\n", "\n", "According to [this page](https://health.howstuffworks.com/human-body/systems/respiratory/question98.htm), the average adult breathes around 7–8 L of air per minute. Starting with air of 20% oxygen, this produces air with 15% oxygen.\n", "\n", "Let's interpret that as a per-minute oxygen consumption of\n", "\n", "$$\\frac{\\mathrm{d}V_{\\mathrm{ox}}}{\\mathrm{d}t} = - V'_0 \\times C_{\\mathrm{ox}} \\times \\eta \\qquad \\mathrm{where}~~\\eta = \\frac{1}{4} ~~\\mathrm{and}~~ V'_0 = 7.5\\;\\mathrm{L\\;min}^{-1}$$\n", "\n", "Now given that\n", "\n", "$$C_{\\mathrm{ox}} \\equiv \\frac{V_{\\mathrm{ox}}}{V_\\Omega}$$\n", "\n", "this results in the straightforward differential equation\n", "\n", "$$ \\frac{\\mathrm{d}C_{\\mathrm{ox}}}{\\mathrm{d}t} = - \\frac{ N_\\mathrm{people} V'_\\mathrm{ox} \\eta } { V_\\Omega } \\times C_{\\mathrm{ox}}$$\n", "\n", "\n", "$$C_\\mathrm{ox}(t) = C_\\mathrm{ox}(0) \\exp(-kt) \\qquad \\mathrm{with~~} k = \\frac{ N_\\mathrm{people} V'_\\mathrm{ox} \\eta } { V_\\Omega } $$" ] }, { "cell_type": "code", "execution_count": 1, "metadata": { "ExecuteTime": { "end_time": "2018-08-24T11:00:35.052841Z", "start_time": "2018-08-24T11:00:34.856480Z" } }, "outputs": [], "source": [ "%matplotlib inline\n", "import numpy as np\n", "from matplotlib import pyplot as plt" ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "ExecuteTime": { "end_time": "2018-08-24T11:00:35.297906Z", "start_time": "2018-08-24T11:00:35.056480Z" } }, "outputs": [], "source": [ "import pint\n", "ureg = pint.UnitRegistry()" ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "ExecuteTime": { "end_time": "2018-08-24T11:01:13.124812Z", "start_time": "2018-08-24T11:01:13.118199Z" } }, "outputs": [ { "data": { "text/html": [ "0.00997340425531915 1/minute" ], "text/latex": [ "$0.00997340425531915\\ \\frac{1}{\\mathrm{minute}}$" ], "text/plain": [ "" ] }, "execution_count": 3, "metadata": {}, "output_type": "execute_result" } ], "source": [ "N = 11\n", "V_prime = 7.5 * ureg.L / ureg.min\n", "η = 1/4\n", "V_person = 0.1 * ureg.metre**3\n", "V_room = 1.2 * 1.2 * 2.2 * ureg.metre**3\n", "V_air = V_room - N * V_person\n", "\n", "k = N * V_prime * η / V_air\n", "\n", "k.to('min^-1')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "> Oxygen deprivation is one form of asphyxiation. While it is desirable to maintain the atmospheric oxygen level at 21% by volume, the body can tolerate deviation from this ideal. **When the oxygen level falls to 17%, the first sign of hypoxia is a deterioration to night vision which is not noticeable until a normal oxygen concentration is restored. Physiologic effects are increased breathing volume and accelerated heartbeat.** Between 14-16% physiologic effects are increased breathing volume, accelerated heartbeat, very poor muscular coordination, rapid fatigue, and intermittent respiration. Between 6-10% the effects are nausea, vomiting, inability to perform, and unconsciousness. **Less than 6%, spasmatic breathing, convulsive movements, and death in minutes.**\n", "\n", "https://web.archive.org/web/20040405214417/http://www.osha.gov/SLTC/smallbusiness/sec12.html\n" ] }, { "cell_type": "code", "execution_count": 4, "metadata": { "ExecuteTime": { "end_time": "2018-08-24T11:03:34.597308Z", "start_time": "2018-08-24T11:03:34.393861Z" } }, "outputs": [ { "data": { "text/plain": [ "Text(0.5,1,'Number of people: 11')" ] }, "execution_count": 4, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZYAAAEWCAYAAABFSLFOAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzs3Xd4FdX28PHvIlTpIlx676EEUoBLCUWQIogCGkQEFbhRUfyhSPGl4xUpFsSGKFIDSFdARSA0aQHpgtKrgCAlUiSw3j/OkHsIgeRAkpOyPs8zT87M7D2zJuhZ2bNn9hZVxRhjjEko6bwdgDHGmNTFEosxxpgEZYnFGGNMgrLEYowxJkFZYjHGGJOgLLEYY4xJUJZYTJolIl+LyDAvnVtEZIKI/CUiG7wRQ2xEZJCITPF2HCZls8Rikg0ROSgiJ0Ukq9u2LiIS7sWwEksdoDFQWFWDvB1MYhCR7iISISJXReTrGPsyisgs599cRaS+d6I0icESi0lu0gM9vB2Ep0TEx8MqxYCDqvp3YsSTTBwHhgFf3WH/auAZ4I8ki8gkCUssJrkZCbwhIrli7hCR4s5ft+ndtoWLSBfnc2cRWSMi74vIORHZLyL/drYfEZFTItIpxmEfEpElInJRRFaISDG3Y5d39p0VkT0i8qTbvq9F5FMRWSQifwMNYom3oIgscOrvFZGuzvYXgPFALRGJFJHBsdS9eS0fich5EdktIo3c9ucUkS9F5ISIHBORYTeTm4ikE5H/JyKHnGueJCI5Y/wOu4nIcaf+63f6xxCRmiLys/P73OpJy0JV56jqPOBMLPv+UdUPVHU1cD2+xzQpgyUWk9xEAOHAG/dYvwawDcgDTAOmA4FAaVx/HY8VkWxu5TsAQ4GHgC3AVADndtwS5xj5gPbAJyLi61b3aeBtIDuuv75jCgOOAgWBtsB/RaSRqn4JhAJrVTWbqg68y7Xsd2IbCMwRkQedfROBKOe6qgFNgC7Ovs7O0gAoCWQDxsY4dgOgjFOvj4g8HPPkIlIIWIir1fEgrn+T2SKS19nfR0S+u0PsJg2zxGKSowHAKze/wDx0QFUnqOp1YAZQBBiiqldV9UfgH1xfxjctVNWVqnoVeAtXK6II8CiuW1UTVDVKVTcDs3EliJvmq+oaVb2hqlfcg3COUQforapXVHULrlZKRw+u5RTwgapeU9UZwB6ghYj8C2gGvKaqf6vqKeB9IMSp1wF4T1X3q2ok0BcIcW/pAYOdutuBCbgSZ0zPAItUdZFzjUtwJf7mAKo6XFUf9eB6TBqRPu4ixiQtVd3h/CXcB/jVw+on3T5fdo4Xc5t7i+WI23kjReQsrhZGMaCGiJxzK5semBxb3VgUBM6q6kW3bYeAgPhchOOY3jpK7CG32DIAJ0Tk5r50bvEUdMq610sP/OsOsR8CKsdy/mJAOxFp6bYtA7Dcg2swaZAlFpNcDQQ2A6Pdtt3s6H4AuOB8zn+f5yly84Nzi+xBXJ3OR4AVqtr4LnXvNjT4ceBBEcnullyKAsc8iK2QiIhbcikKLHBiuwo8pKpRdzh3Mbf1orhum50ECjvbigC73fYfj+U4R4DJqtrVg5iNsVthJnlS1b24bmW96rbtNK4v5mdExEdEngdK3eepmotIHRHJiKuvZb2qHgG+A8qKSEcRyeAsgSJSIZ7xHwF+Bt4RkcwiUgV4AacPJ57yAa86524HVMB1a+oE8CMwWkRyOJ31pUQk2KkXBvyfiJRwkuV/gRkxklB/EXnA6TN6DtfvOqYpQEsRecT5fWcWkfoiUjiWsrcRkfQikhnwAW7Wd3/wIpOzHyCjs19iPZhJUSyxmORsCJA1xrauQC9cTxr54vryvh/TcLWOzgL+uPoncFoZTXD1WxzH9Ujsu0AmD47dHiju1J8LDHT6KeJrPa4O9j9xPSTQVlVvPmH1LJAR2AX8BcwCCjj7vsJ1y24lcAC4ArwS49grgL3AUmCU0/90Cyc5Pgb0A07jasH0wvneEJF+IrL4LvH/P1y3Hvvg6q+57Gy7aY+zrRDwg/O5GCbFE5voy5jkR0Q6A11UtU4CH7c4rmST4Q630Yy5b9ZiMcYYk6AssRhjjElQdivMGGNMgrIWizHGmASVJt9jeeihh7R48eLeDsMYY1KMTZs2/amq8RoNI00mluLFixMREeHtMIwxJsUQkUNxl3KxW2HGGGMSlCUWY4wxCcoSizHGmARlicUYY0yCssRijDEmQSVpYhGRps4Ur3tFpE8s+3uKyC4R2SYiS2NME9tJRH53lk5u2/1FZLtzzDE2OqoxxnhXkiUWZz7uj3HNfFcRaC8iFWMU+wUIUNUquEZrHeHUfRDXCLQ1gCBgoIjkdup8CnTDNQpsGaBpIl+KMcaYu0jKFksQsNeZLvUfXHORP+ZeQFWXq+olZ3Ud/5uU6BFgiaqeVdW/cM1F3lRECgA5VHWtMxnSJKB1YgR/48YNPvzwQ5Yvt8nzjDHmbpLyBclC3Dod6lFcLZA7eQG4OddDbHULOcvRWLbfRkS64WrZULRoUU/iBuDq1auMHTuW69evs23bNrJlyxZ3JWOMSYOSssUSW99HrCNgisgzuOYGHxlH3XgfU1XHqWqAqgbkzRuvUQlukSVLFr766isOHjxIv379PK5vjDFpRVImlqO4zS+O6zbXbfNsi8jDwFtAK1W9Gkfdo/zvdtkdj5lQ6tatS/fu3fnoo49YtWpVYp3GGGNStKRMLBuBMs483BlxTfm6wL2AiFQDPseVVE657foBaCIiuZ1O+ybAD87c3xdFpKbzNNizwPzEvIh33nmHEiVK8Pzzz3Pp0qW4KxhjTBqTZInFmQa1O64k8SswU1V3isgQEWnlFBsJZAO+EZEtIrLAqXsWGIorOW0EhjjbAF4ExuOav3sf/+uXSRRZs2Zl/Pjx7N27lwEDBiTmqYwxJkVKkxN9BQQE6P2ObhwaGsoXX3zBmjVrqFmzZgJFZowxyZOIbFLVgPiUtTfv79GIESMoVKgQzz33HFeuXPF2OMYYk2xYYrlHOXLk4IsvvmD37t0MHjzY2+EYY0yyYYnlPjzyyCM8//zzjBw50iYOM8YYhyWW+zR69Gj+9a9/0alTJ7slZowxWGK5b7ly5eLLL79k165d9O/f39vhGGOM11liSQBNmzYlNDSU0aNHs3LlSm+HY4wxXmWJJYGMHDmSkiVL0qlTJy5cuODtcIwxxmsssSSQbNmyMXHiRA4fPkzPnj29HY4xxniNJZYEVLt2bd58802+/PJLvv32W2+HY4wxXmFv3iewq1evEhQUxB9//MGOHTu4l5GUjTEmubE3770oU6ZMTJ48mXPnzhEaGkpaTNzGmLTNEksiqFKlCkOHDmXOnDlMmTLF2+EYY0ySssSSSF5//XVq165N9+7dOXDggLfDMcaYJGOJJZH4+PgwefJkAJ555hmioqK8HJExxiQNSyyJqESJEnz66af8/PPPDBs2zNvhGGNMkrDEksiefvppOnbsyNChQ1mzZo23wzHGmERniSUJjB07luLFi9OhQwfOnTvn7XCMMSZRWWJJAjly5GDq1KkcPXqUF1980R5BNsakapZYkkjNmjUZPHgw06dPj+7UN8aY1ChJE4uINBWRPSKyV0T6xLK/nohsFpEoEWnrtr2BiGxxW66ISGtn39cicsBtn19SXpMn+vTpQ7169Xj55ZfZt2+ft8MxxmvGjBlDhQoV6NChAwsWLGD48OHxrnvw4EGmTZvm0fmaN2/uldvQ33zzDRUqVKBBgwa3bA8PD+fRRx9N8ngAunTpwq5duxL3JKqaJAvgA+wDSgIZga1AxRhligNVgElA2zsc50HgLPCAs/71ncreafH391dvOXTokObKlUuDgoL06tWrXovDGG8qV66c7t+//65lrl27Fuv25cuXa4sWLRIjrAT3yCOP6LJly27bnpKu4SYgQuP5HZuULZYgYK+q7lfVf4DpwGPuBVT1oKpuA27c5ThtgcWqeinxQk08RYsWZfz48WzYsIG+fft6OxxjklxoaCj79++nVatWvP/++3z99dd0794dgM6dO9OzZ08aNGhA7969WbFiBX5+fvj5+VGtWjUuXrxInz59WLVqFX5+frz//vu3HPvEiRPUq1cPPz8/KlWqxKpVqwAoXrw4f/75JwBDhw6lfPnyNG7cmPbt2zNq1CgA6tevT+/evQkKCqJs2bLRdXfu3ElQUBB+fn5UqVKF33///bZrCgsLo3LlylSqVInevXsDMGTIEFavXk1oaCi9evW6rU5kZCRt27alfPnydOjQAVVl6dKlPP7449FllixZwhNPPAG4RlB//fXXqV69Oo0aNeL06dMAbNmyhZo1a1KlShUef/xx/vrrL6KioggMDCQ8PByAvn378tZbb0Vf582xEl988UUCAgLw9fVl4MCBnv5T3ll8M9D9LrgSwni39Y7A2DuU/Zo7t1iWAY/GKLsH2Aa8D2S6Q71uQAQQUbRo0QTJ4Pfj5ZdfVkDnz5/v7VCMSXLFihXT06dPq6rqhAkT9OWXX1ZV1U6dOmmLFi00KipKVVUfffRRXb16taqqXrx4Ua9du3bXv/ZHjRqlw4YNU1XVqKgovXDhwi3n27hxo1atWlUvXbqkFy5c0NKlS+vIkSNVVTU4OFh79uypqqoLFy7URo0aqapq9+7ddcqUKaqqevXqVb106dIt5zx27JgWKVJET506pdeuXdMGDRro3Llzo4+5cePG2+Jcvny55siRQ48cOaLXr1/XmjVr6qpVq/TGjRtarlw5PXXqlKqqtm/fXhcsWKCqqkB0HIMHD47+nVWuXFnDw8NVVbV///7ao0cPVVXdsWOHli9fXn/88Uf18/OLvkPiHtOZM2eif1fBwcG6devWO/6bkUxbLBLLNo8ejxKRAkBl4Ae3zX2B8kAgrttkvWOrq6rjVDVAVQOSw4jDo0aNolq1anTu3JlDhw55Oxxjko127drh4+MDuKai6NmzJ2PGjOHcuXOkT5/+rnUDAwOZMGECgwYNYvv27WTPnv2W/atXr+axxx4jS5YsZM+enZYtW96y/2brwN/fn4MHDwJQq1Yt/vvf//Luu+9y6NAhsmTJckudjRs3Ur9+ffLmzUv69Onp0KFDvGaSDQoKonDhwqRLlw4/Pz8OHjyIiNCxY0emTJnCuXPnWLt2Lc2aNQMgXbp0PPXUU4BrNI/Vq1dz/vx5zp07R3BwMACdOnWKPrevry8dO3akZcuWfPXVV2TMmPG2GGbOnEn16tWpVq0aO3fuTLC+l6RMLEeBIm7rhYHjHh7jSWCuql67uUFVTzgJ9SowAdctt2Qvc+bMzJw5k6ioKEJCQrh27VrclYxJA7JmzRr9uU+fPowfP57Lly9Ts2ZNdu/efde69erVY+XKlRQqVIiOHTsyadKkW/ZrHI/6Z8qUCXANyXRzGKann36aBQsWkCVLFh555BGWLVvm0THjOlfM8z333HNMmTKFsLAw2rVrd8dkKhLb3+q32r59O7ly5eLkyZO37Ttw4ACjRo1i6dKlbNu2jRYtWnDlypV7upaYkjKxbATKiEgJEckIhAALPDxGeyDMfYPTikFcv+XWwI4EiDVJlC5dmvHjx7Nu3Tr69evn7XCMSXb27dtH5cqV6d27NwEBAezevZvs2bNz8eLFWMsfOnSIfPny0bVrV1544QU2b958y/46derw7bffcuXKFSIjI1m4cGGcMezfv5+SJUvy6quv0qpVK7Zt23bL/ho1arBixQr+/PNPrl+/TlhYWHQL4l4ULFiQggULMmzYMDp37hy9/caNG8yaNQuAadOmUadOHXLmzEnu3Lmj+4MmT54cfe45c+Zw5swZVq5cyauvvnrbU3EXLlwga9as5MyZk5MnT7J48eJ7jjmmu7crE5CqRolId1y3sXyAr1R1p4gMwXXvboGIBAJzgdxASxEZrKq+ACJSHFeLZ0WMQ08Vkby4brVtAUKT5IISyJNPPkl4eDijRo0iODjYa48gGpMcffDBByxfvhwfHx8qVqxIs2bNSJcuHenTp6dq1ap07tyZ//u//4suHx4ezsiRI8mQIQPZsmW7rcUSGBhIq1atqFq1KsWKFSMgIICcOXPeNYYZM2YwZcoUMmTIQP78+RkwYMAt+wsUKMA777xDgwYNUFWaN2/OY489doejxU+HDh04ffo0FStWjN6WNWtWdu7cib+/Pzlz5mTGjBkATJw4kdDQUC5dukTJkiWZMGECf/75J3369GHp0qUUKVKE7t2706NHDyZOnBh9vKpVq1KtWjV8fX0pWbIktWvXvq+Y3dkMksnAlStXqFWrFocPH+aXX36haNGi3g7JmFQrMjKSbNmycenSJerVq8e4ceOoXr26t8O6Rffu3alWrRovvPBC9LZs2bIRGRnptZg8mUEyTSaWglJQ/8N/vB3GLc5whs/5nHzkozOdSZ90jUlj0pRZzOI0p4kiCj/8qEtdb4d0i8/5nAxk4FmeveV74G3e5i3e8lpcgxhkieVukluL5aYZM2YQEhLCq6++yocffujtcIwxJprNeZ9CPfXUU/To0YMxY8YQFhYWdwVjjEmGLLEkMyNHjqR27dp06dKFHTtSzANuxhgTzRJLMpMhQwa++eYbcuTIQZs2bTh//ry3QzLGGI9YYkmGChQowMyZM9m3bx/PPfeczd9ijElRLLEkU3Xr1mXkyJHMnTuXkSNHejscY4yJN0ssydhrr73Gk08+Sd++fW8bRsIYY5IrSyzJmIjw5ZdfUq5cOUJCQjh8+LC3QzLGmDhZYknmsmXLxty5c7l69SqPP/44ly6lyGlojDFpiCWWFKBcuXJMmzaNX375hS5dulhnvjEmWbPEkkK0aNGC//73v4SFhTFixAhvh2OMMXdkiSUF6d27N0899RR9+/Zl0aJF3g7HGGNiZYklBRERvvrqK/z8/Gjfvj179uzxdkjGGHMbSywpzAMPPMC8efPIlCkTrVq1um3yHmOM8TZLLClQ0aJFmT17Nvv376dDhw5cv37d2yEZY0w0SywpVN26dRk7diyLFi3izTff9HY4xhgTzWaTSsH+85//sGvXLt577z3KlClDaGiKmpXZGJNKWWJJ4d577z327dtH9+7dKVmyJE2aNPF2SMaYNC5Jb4WJSFMR2SMie0WkTyz764nIZhGJEpG2MfZdF5EtzrLAbXsJEVkvIr+LyAwRyZgU15Jc+Pj4EBYWhq+vL+3atWPnzp3eDskYk8bFO7GISCYReVpE+onIgJuLB/V9gI+BZkBFoL2IVIxR7DDQGZgWyyEuq6qfs7Ry2/4u8L6qlgH+Al6Ib0ypRfbs2fn222954IEHePTRRzl16pS3QzLGpGGetFjmA48BUcDfbkt8BQF7VXW/qv4DTHeOF01VD6rqNuBGfA4oIgI0BGY5myYCrT2IKdUoWrQoCxYs4OTJkzz22GNcvnzZ2yEZY9IoT/pYCqtq0/s4VyHgiNv6UaCGB/Uzi0gErsQ2XFXnAXmAc6oa5XbMQvcRY4oWGBjIlClTaNOmDc899xzTpk0jXTp78M8Yk7Q8+db5WUQq38e5JJZtnoymWFRVA4CngQ9EpJQnxxSRbiISISIRp0+f9uC0KcsTTzzB8OHDmTFjBoMGDfJ2OMaYNMiTxFIH2OR0vm8Tke0iss2D+keBIm7rhYHj8a2sqsedn/uBcKAa8CeQS0RutrzueExVHaeqAaoakDdvXg/CTnnefPNNnn/+eYYOHcpXX33l7XCMMWmMJ7fCmt3nuTYCZUSkBHAMCMHV+oiTiOQGLqnqVRF5CKgNjFBVFZHlQFtcfTadcPUFpWkiwqeffsqRI0fo1q0b+fPnp3nz5t4OyxiTRsS7xaKqh4BcQEtnyeVsi2/9KKA78APwKzBTVXeKyBARaQUgIoEichRoB3wuIjefna0ARIjIVmA5rj6WXc6+3kBPEdmLq8/ly/jGlJplzJiR2bNnU6VKFdq1a8fGjRu9HZIxJo2Q+E4aJSI9gK7AHGfT48A4Vf0okWJLNAEBARoREeHtMJLEH3/8Qa1atfj777/5+eefKV26tLdDMsakQCKyyennjpMnfSwvADVUdYCqDgBq4ko0JhnLnz8/33//PTdu3KBp06b2josxJtF5klgEcB9G9zqxP5Vlkply5crx3Xffcfz4cVq0aEFkZKS3QzLGpGKeJJYJwHoRGSQig4B1WH9GilGzZk1mzJjB5s2befLJJ7l27Zq3QzLGpFKedN6/BzwPnMU1dMpzqvpBYgVmEl7Lli357LPPWLx4Md26dSO+/WvGGOMJj0Y3VtVNwKZEisUkga5du3Ls2DEGDx5M3rx5effdd3GNjGOMMQkjzsQiIqtVtY6IXOTWt9oFUFXNkWjRmUQxcOBATp8+zciRI8mdOzd9+/b1dkjGmFQkzsSiqnWcn9kTPxyTFESEjz76iPPnz9OvXz9y5crFiy++6O2wjDGphCfD5r8bn20mZUiXLh0TJkygZcuWvPzyy0ydOtXbIRljUglPngprHMu2+x3mxXhRhgwZmDlzJvXr16dTp058++233g7JGJMKxJlYRORFEdkOlHMGn7y5HAA8GYTSJEOZM2dm/vz5VK9enXbt2hEeHu7tkIwxKVx8WizTcI0NtoD/jRPWEvBX1WcSMTaTRLJnz87ixYspVaoULVu2tHHFjDH3Jc7EoqrnnZkd26vqIbflbFIEaJJGnjx5WLJkCXnz5qVp06Zs377d2yEZY1IoTzrvJ4pILrf13CJik32kIgULFuSnn34iS5YsNGrUiF27dsVdyRhjYvCk876Kqp67uaKqf+GabMukIiVLlmT58uWkT5+ehg0bsnv3bm+HZIxJYTxJLOmcCbcAEJEH8fDNfZMylClThmXLlgHQsGFDfv/9dy9HZIxJSTxJLKNxzXs/VESGAj8DIxInLONt5cuXZ+nSpURFRdGgQQP27dvn7ZCMMSmEJ4NQTgLaACeBU8ATqjo5sQIz3ufr68vSpUu5cuUKDRo04MCBA94OyRiTAnjSYgE4AWwAtgIPiUi9hA/JJCeVK1fmp59+IjIykoYNG3L48GFvh2SMSeY8eSqsC7AS15z1g52fgxInLJOc+Pn5sWTJEv766y8aNGhgycUYc1eetFh6AIHAIVVtgOuJsNOJEpVJdvz9/fnxxx85c+YM9erVY//+/d4OyRiTTHmSWK6o6hUAEcmkqruBcp6cTESaisgeEdkrIn1i2V9PRDaLSJSItHXb7icia0VkpzOczFNu+74WkQMissVZ/DyJycRfUFAQS5cu5eLFi9SrV4/ffvvN2yEZY5IhTxLLUecFyXnAEhGZDxyPb2UR8QE+xjVwZUWgvYhUjFHsMNAZ1zAy7i4Bz6qqL9AU+MD9ZU2gl6r6OcsWD67JeMjf35/ly5fzzz//EBwcbC9RGmNuE6/EIq4pBl9V1XOqOgjoj2u++9YenCsI2Kuq+1X1H2A68Jh7AWfomG3AjRjbf1PV353Px3E9lZbXg3ObBFSlSpXowSqDg4PZunWrdwMyxiQr8Uos6pocfZ7b+gpVXeAkiPgqBBxxWz/qbPOIiAQBGQH3Fyvedm6RvS8ime5Qr5uIRIhIxOnT1jV0vypWrMjKlSvJnDkzDRo0YNMmm7HaGOPiya2wdSISeB/nim1idY1l250PIFIAmAw8p6o3WzV9gfK4Hix4EOgdW11VHaeqAaoakDevNXYSQpkyZVi5ciU5c+akUaNGrFu3ztshGWOSAU8SSwNgrYjsc1oH20XEk/lYjgJF3NYL41kfTQ5gIfD/VDX6G0xVT6jLVWACrltuJomUKFGCFStW8NBDD9G4cWOWL1/u7ZCMMV4Wn4m+br5dPw4oBTTENR/Lo87P+NoIlBGREiKSEQjBNcdLnJzyc4FJqvpNjH0FnJ+Cq89nhwcxmQRQtGhRVq5cSbFixWjatCnz5s2Lu5IxJtWKT4vFX0SKAc8BF4CLMZZ4UdUooDuuFyt/BWaq6k4RGSIirQBEJFBEjgLtgM9FZKdT/UmgHtA5lseKpzozXG4HHgKGxTcmk3AKFizIypUrqVatGm3atGHChAneDskY4yXi6pe/SwGRV4EXgZLAMW7tK1FVLZl44SWOgIAAjYiI8HYYqVJkZCRPPPEES5YsYfTo0fTs2dPbIRljEoCIbFLVgPiUjc8MkmNUtQLwlaqWVNUSbkuKSyomcWXLlo1vv/2Wdu3a8frrr9OvXz/i+uPFGJO6xHs+FVV9MTEDMalHpkyZCAsLI1euXLzzzjucPXuWjz/+GB8fH2+HZoxJAjZRl0kUPj4+fP755+TJk4fhw4dz9uxZJk+eTKZMsb5mZIxJRSyxmEQjIrzzzjvkyZOHXr16cerUKebNm0euXLnirmyMSbE8nY/FGI+98cYbTJkyhZ9//pk6derYsPvGpHKezMciIvKMiAxw1os6w6sYE6cOHTrw/fffc+TIEWrVqmXjixmTinnSYvkEqAW0d9Yv4hqt2Jh4adiwIatXryZdunTUrVuXJUuWeDskY0wi8CSx1FDVl4ErAKr6F67BII2Jt8qVK7N27VpKlChB8+bNmTRpkrdDMsYkME8SyzVnThUFEJG8xBje3pj4KFy4MCtXriQ4OJhOnToxbNgwe9fFmFTEk8QyBtd4Xf8SkbeBNcA7iRKVSfVy5szJokWL6NixI/3796dz585cvXrV22EZYxKAJy9IThWRTUAjZ1MrZ3piY+5JxowZmThxIqVLl2bgwIHs37+fOXPmYNMaGJOyefJUWACuFkoXIBSY6eGw+cbcRkQYMGAAM2bMICIigho1ath0x8akcJ7cCpuKa76TJ/jfkPmeDJtvzB09+eSTrFixgkuXLlGrVi1++OEHb4dkjLlHniSW0850xAdU9dDNJdEiM2lOUFAQGzdujH5ibOzYsd4OyaQgn332WfRThl9//TXHj/9vHsEuXbokm5ZwtmzZvB1Cootz2PzogiKNcL3DshSI7mVV1TmJE1risWHzk7fIyEg6dOjAggULeOmll/jggw/IkCGDt8MyKUj9+vUZNWoUAQHxGuU9SWXLlo3IyEhvh+GxBB02381zgB/QlP/dBnvU8/CMubts2bIxZ84cevXqxSeffELjxo05vsjSAAAgAElEQVQ5deqUt8MyCaR379588skn0euDBg1i9OjRqCq9evWiUqVKVK5cmRkzZgAQHh5OcHAwTz75JGXLlqVPnz5MnTqVoKAgKleuzL59+6KPM2rUKGbNmkVERAQdOnTAz8+Py5cvU79+fW7+MZktWzbeeustqlatSs2aNTl58iQA+/bto2bNmgQGBjJgwIA7tixat26Nv78/vr6+jBs3Lnr7nY574MABatWqRWBgIP3794/1mAcPHqRChQp07doVX19fmjRpwuXLlwH44osvCAwMpGrVqrRp04ZLly4B0LlzZ2bNmnXL+QFOnDhBvXr18PPzo1KlSqxatQqAH3/8kVq1alG9enXatWsXndyGDBlCYGAglSpVolu3bgnz6L+qxmsBtse3bHJf/P391aQMU6ZM0cyZM2uRIkV048aN3g7HJIDNmzdrvXr1otcrVKighw4d0lmzZunDDz+sUVFR+scff2iRIkX0+PHjunz5cs2ZM6ceP35cr1y5ogULFtQBAwaoquoHH3ygPXr0UFXVgQMH6siRI1VVNTg4+Jb/XtzXAV2wYIGqqvbq1UuHDh2qqqotWrTQadOmqarqp59+qlmzZo01/jNnzqiq6qVLl9TX11f//PPPux63ZcuWOnHiRFVVHTt2bKzHPXDggPr4+Ogvv/yiqqrt2rXTyZMnq6pGH19V9a233tIxY8aoqmqnTp30m2++id5387ijRo3SYcOGqapqVFSUXrhwQU+fPq1169bVyMhIVVUdPny4Dh48+JbrUVV95plnoq8hJiBC4/kd60mLZZ2IVLz/VGZM/HXo0IE1a9YgItSpU8fe1E8FqlWrxqlTpzh+/Dhbt24ld+7cFC1alNWrV9O+fXt8fHz417/+RXBwMBs3bgQgMDCQAgUKkClTJkqVKkWTJk0A10gOBw8e9Oj8GTNm5NFHXTdb/P39o+uvXbuWdu3aAfD000/fsf6YMWOiWyVHjhzh999/v+tx16xZQ/v2rpGwOnbseMfjlihRAj8/v9vq79ixg7p161K5cmWmTp3Kzp0773gMcP2uJkyYwKBBg9i+fTvZs2dn3bp17Nq1i9q1a+Pn58fEiRM5dMjVRb58+XJq1KhB5cqVWbZsWZzHjw9Phs2vA3QSkQO4+lgE19TEVe47CmPuonr16kRERPDUU0/RqVMnNm3axKhRo6zfJQVr27Yts2bN4o8//iAkJATgrrdg3OfxSZcuXfR6unTpiIqK8ujcGTJkQMQ1w7qPj49H9cPDw/npp59Yu3YtDzzwAPXr1+fKlStxHvfm9rtxv0YfH5/oW2GdO3dm3rx5VK1ala+//prw8HAA0qdPz40brsFPVJV//vkHgHr16rFy5UoWLlxIx44d6dWrF7lz56Zx48aEhYXdcs4rV67w0ksvERERQZEiRRg0aFD09dwPT1osTYEyQBP+17/i0ePGItJURPaIyF4R6RPL/noisllEokSkbYx9nUTkd2fp5LbdX0S2O8ccI/H5FzQpTt68efnxxx957bXXGDNmjPW7pHAhISFMnz6dWbNm0bat63/1evXqMWPGDK5fv87p06dZuXIlQUH3NoB69uzZuXjxokd1atasyezZswGYPn16rGXOnz9P7ty5eeCBB9i9ezfr1q2L87i1a9eOPt7UqVM9igng4sWLFChQgGvXrt1Sv3jx4mzatAmA+fPnc+3aNQAOHTpEvnz56Nq1Ky+88AKbN2+mZs2arFmzhr179wJw6dIlfvvtt+gk8tBDDxEZGXlLn839iHdiUbdHjPUeHjd2xhn7GGgGVATax3Jr7TDQGZgWo+6DwECgBhAEDBSR3M7uT4FuuJJeGVwJ0KRC6dOn5/3332fSpEmsX7+egIAA1q9f7+2wzD3w9fXl4sWLFCpUiAIFCgDw+OOPU6VKFapWrUrDhg0ZMWIE+fPnv6fjd+7cmdDQ0OjO+/j44IMPeO+99wgKCuLEiRPkzJnztjJNmzYlKiqKKlWq0L9/f2rWrBnncT/88EM+/vhjAgMDOX/+vMfXMnToUGrUqEHjxo0pX7589PauXbuyYsUKgoKCWL9+PVmzZgVcrSo/Pz+qVavG7Nmz6dGjB3nz5uXrr7+mffv2VKlShZo1a7J7925y5cpF165dqVy5Mq1btyYwMNDj+GLjyePGA2LbrqpD4lm/FjBIVR9x1vs69W8bb0xEvga+U9VZznp7oL6q/sdZ/xwId5blqlo+tnJ3Yo8bp3ybNm2ibdu2HDt2jNGjR9O9e/d43W4w5k4uXbpElixZEBGmT59OWFgY8+fP93ZYyYYnjxt70sfyt9vnzLhuhf3qQf1CwBG39aO4WiD3WreQsxyNZfttRKQbrpYN/iWAafYllJL5A5v7wbOfwquvvsrqqa8yvitkz+LtyExKtWk3dJ8IqpDrAfiqG/Y9cY88GYRytPu6iIwCFnhwrtj+heL7wPSd6sb7mKo6DhgHrhYLT1uLJaXLDcx/4QYjR46kX79+bD1XhlmzZlGpUiVvh2ZSoLrA1njdf0mjOsQ/yd7PnPcPACU9KH8UKOK2Xhg4foey8a171Pl8L8c0qUC6dOno3bs3y5Yt4/z58wQFBdkjycZ4mSejG28XkW3OshPYA3zowbk2AmVEpISIZARCiH+L5wegiYjkdjrtmwA/qOoJ4KKI1HSeBnsWsJuiaVBwcDC//PILNWrUoFOnTnTt2jX6DWVjTNLypMXiPqJxE6CgqsZ7lEBVjQK640oSvwIzVXWniAwRkVYAIhIoIkeBdsDnTgJDVc8CQ3Elp43AEGcbwIvAeGAvsA9Y7ME1mVQkf/78LFmyhL59+zJ+/HgCAwPZvn27t8MyJs2J91NhqYk9FZb6LVmyhI4dO3Lu3DlGjx7NSy+9ZE+NGXMfEnQQShG5KCIX3JaL7j/vP1xjEl7jxo3Ztm0bDRs2pHv37rRu3ZozZ854Oyxj0oQ4E4uqZlfVHG5LdvefSRGkMfciX758fPfdd7z//vssXryYqlWrRg+HYYxJPB49FSYiVUWku7PYGGEm2UuXLh2vvfYa69atI2vWrDRs2JD+/ft7PL6UMSb+PHkqrAeu6YnzOctUEXklsQIzJiFVr16dTZs20blzZ4YNG0bdunWjR6U1xiQsT1osLwA1VHWAqg4AagJdEycsYxJetmzZ+OqrrwgLC2P37t34+fnx6aefJszERsaYaJ4kFgGuu61fJ/Y3341J1kJCQtixYwd16tThpZdeolmzZhw7dszbYRmTaniSWCYA60VkkIgMAtYBXyZKVMYkskKFCvH999/z8ccfs3LlSipVqnTbXBXGmHvjybD57+Ga9/4s8BfwnKp+kFiBGZPYRISXXnqJrVu3Ur58eZ5++mlCQkI4e/Zs3JWNMXfkSef9/wGnVHWMqn6oqr8kYlzGJJkyZcqwatUq3n77bWbPnk2lSpVYtGiRt8MyJsXy5FZYDuAHEVklIi+LyL8SKyhjklr69Onp168fGzZsIE+ePLRo0YJnn33WXqo05h54citssKr6Ai8DBYEVIvJTokVmjBdUq1aNiIgI+vfvT1hYGBUrVkyw6VqNSSvuZdj8U8AfwBlc77MYk6pkypSJIUOGEBERQZEiRWjXrh1t2rThjz/+8HZoxqQInvSxvCgi4cBS4CGgq6ra2/cm1apatSrr1q1j+PDhLFy4kIoVKzJx4kR778WYOHjSYikGvKaqvqo6UFV3JVZQxiQX6dOnp3fv3mzduhVfX186d+5Ms2bNOHTokLdDMybZ8qSPpY+qbknMYIxJrsqVK8eKFSv46KOPWL16NRUrVmTEiBFcu3bN26EZk+zcz9TExqQp6dKlo3v37uzatYsmTZrQu3dvqlevzpo1a7wdmjHJiiUWYzxUtGhR5s6dy/z58zl//jx16tSha9eu9miyMQ5LLMbco1atWrFr1y569erFhAkTKF++vHXuG8M9JBYRaSwiX4iIn7PeLeHDMiZlyJYtGyNGjGDz5s2ULVuWzp0706BBA3bu3Ont0IzxmntpsbwE9AKeEZGGgF98K4pIUxHZIyJ7RaRPLPszicgMZ/96ESnubO8gIlvclhtuiS3cOebNffZujUlyVapUYdWqVXzxxRds27aNqlWr0qNHD/766y9vh2ZMkruXxHJaVc+p6htAEyAwPpVExAf4GGgGVATai0jFGMVeAP5S1dLA+8C7AKo6VVX9VNUP6AgcjPGEWoeb+1X11D1ckzH3LV26dHTp0oXff/+dbt26MXbsWMqWLcsXX3zB9evX4z6AMamEJy9IZhaRSsAhEckMrkeQgUnxPEQQsFdV96vqP8B04LEYZR4DJjqfZwGNRCTmnC/tARvf3CRbefLk4ZNPPmHTpk1UqFCBbt26ERQUZE+PmTQjzsQiIulFZARwFNeXflvgiIiMEJEMqvpRPM9VCDjitn7U2RZrGVWNAs4DeWKUeYrbE8sE5zZY/1gS0c3r6CYiESIScfr06XiGbMy98/PzY8WKFYSFhXHy5Enq1KlDx44dOX78uLdDMyZRxafFMhJ4ECihqv6qWg0oBeQCRnlwrti+8GM+PnPXMiJSA7ikqjvc9ndQ1cpAXWfpGNvJVXWcqgaoakDevHk9CNuYeycihISEsGfPHt566y1mzpxJ2bJlGTp0KH///be3wzMmUcQnsTyKa1ywizc3qOoF4EWguQfnOgoUcVsvDMT80y26jIikB3LimljsphBitFZU9Zjz8yIwDdctN2OSlaxZszJs2DB27drFI488woABAyhbtiwTJkyw/heT6sQnsajG8mC+ql7n9hbH3WwEyohICRHJiCtJLIhRZgHQyfncFlh289wikg5oh6tvBmdbehF5yPmcAVcS3IExyVSpUqWYPXs2q1atonDhwjz//PP4+/vz0082A4VJPeKTWHaJyLMxN4rIM8Du+J7I6TPpDvwA/ArMVNWdIjJERFo5xb4E8ojIXqAn4P5Icj3gqKrud9uWCdfkY9uALcAx4Iv4xmSMt9SpU4d169Yxffp0zp8/T+PGjWnevDk7dtjfRSblk7jeEhaRQsAc4DKwCVcrJRDIAjx+81ZUShIQEKARERHeDsMYAK5evcrYsWMZNmwYFy5c4Pnnn2fgwIEULlzY26EZE01ENqlqQHzKxtliUdVjqloDGAIcBA4DQ1Q1KCUmFWOSm0yZMvH666+zd+9eXnnlFSZOnEjp0qV54403+PPPP70dnjEei7PFkhpZi8UkZwcPHmTw4MFMmjSJrFmz0rNnT3r27EmOHDm8HZpJwxK0xWKMSVrFixdnwoQJ7NixgyZNmjB48GBKlizJqFGjuHz5srfDM4nIx8cHPz8/fH19qVq1Ku+99x43btwAIDw8nJw5c+Ln5xe9uD/0MXfuXESE3btdXd/bt2+PLvfggw9SokQJ/Pz8ePjhhzl48CCVKlW65dyDBg1i1ChP3iC5M0ssxiRTFSpUYNasWWzcuJGAgAB69epF6dKl+fzzz22CsVQqS5YsbNmyhZ07d7JkyRIWLVrE4MGDo/fXrVuXLVu2RC8PP/xw9L6wsDDq1KnD9OmuB2crV64cXa5Vq1aMHDmSLVu2JMkTiJZYjEnmAgIC+P777wkPD6d48eKEhoZSrlw5xo8fzz///OPt8EwiyZcvH+PGjWPs2LFxTsUQGRnJmjVr+PLLL6MTizdZYjEmhQgODmb16tV899135MmTh65du1K2bFk+//xzSzCpVMmSJblx4wanTrnG1l21atUtt8L27dsHwLx582jatClly5blwQcfZPPmzXEee9++fbcc67PPPkuwuC2xGJOCiAgtWrRgw4YNLFy4kPz58xMaGkrp0qX55JNPuHr1qrdDNAnMvbUS81ZYqVKlANdtsJCQEABCQkIIC4t7nN5SpUrdcqzQ0NAEi9kSizEpkIjQvHlz1q5dyw8//ECRIkV4+eWXKVWqFGPHjuXKlSveDtEkgP379+Pj40O+fHeeZurMmTMsW7aMLl26ULx4cUaOHMmMGTO8OpOpJRZjUjARoUmTJqxevZqffvqJkiVL8sorr1CyZElGjx7NxYsX4z6ISZZOnz5NaGgo3bt35w6DtgMwa9Ysnn32WQ4dOsTBgwc5cuQIJUqUYPXq1UkY7a0ssRiTCogIjRo1YuXKlSxfvpzy5cvzxhtvULRoUd566y1Onjzp7RBNPFy+fDn6ceOHH36YJk2aMHDgwOj9MftYZs2aRVhYGI8//vgtx2nTpg3Tpk1L6vCj2QuSxqRSGzZsYMSIEcyZM4eMGTPy3HPP8cYbb0TflzfGE568IJk2E4uIWloxacVvuCZOmghE4Ro2vDdQ3ZtBmRRHwN68vyt/f1C1xZY0sZRVZZwqB48fp1fv3nyfIwf+QOOHH2bxokXcuH7d6zHakgIWD6TNxGJMGlSgQAGGDx/O4cOHeffdd9m5cyfNmzfH19eXzz77zGa0NAnGEosxaUzOnDl58803OXjwIFOmTCFr1qy8+OKLFClShD59+nDkyBFvh2hSOEssxqRRGTNmpEOHDmzcuJHVq1fTqFEjRo4cSYkSJQgJCWHdunXeDtGkUJZYjEnjRITatWvzzTffsG/fPl577TUWL15MrVq1qFGjBpMmTbIXLo1HLLEYY6IVL16cUaNGcfToUcaMGcP58+fp1KkThQoV4o033mDv3r3eDtGkAJZYjDG3yZ49O6+88gq//vorS5cupWHDhnz44YeUKVOGRx55hPnz5xMVFeXtME0ylaSJRUSaisgeEdkrIn1i2Z9JRGY4+9eLSHFne3ERuSwiW5zlM7c6/iKy3akzRu429oExxiMiQsOGDfnmm284dOgQgwcPZufOnbRu3ZoSJUowdOhQTpw44e0wTTKTZIlFRHyAj4FmQEWgvYhUjFHsBeAvVS0NvA+867Zvn6r6OYv7MJyfAt2AMs7SNLGuwZi0rGDBggwYMICDBw8yd+5cKlSowIABAyhatCht2rRh0aJF1ooxQNK2WIKAvaq6X1X/AaYDj8Uo8xiuF4QBZgGN7tYCEZECQA5VXauuIQQmAa0TPnRjzE3p06endevW/Pjjj/z222/06NGDlStX0qJFC4oVK8Zbb71lfTFpXFImlkKA+wPyR51tsZZR1SjgPJDH2VdCRH4RkRUiUtet/NE4jgmAiHQTkQgRiTh9+vT9XYkxBoAyZcowatQojh07xuzZs/Hz82P48OGUKVOG+vXrM2nSJC5duuTtME0SS8rEElvLI+Y4AXcqcwIoqqrVgJ7ANBHJEc9jujaqjlPVAFUNyJs3rwdhG2PikjFjRp544gkWLlzI4cOHefvttzl69CidOnWiQIEChIaGsmHDBtLi2IRpUVImlqNAEbf1wsDxO5URkfRATuCsql5V1TMAqroJ2AeUdcoXjuOYxpgkVKhQIfr168fvv/9OeHg4rVu3ZtKkSdSoUYMKFSowbNgwDhw44O0wTSJKysSyESgjIiVEJCMQAiyIUWYB0Mn53BZYpqoqInmdzn9EpCSuTvr9qnoCuCgiNZ2+mGeB+UlxMcaYuxMRgoODmThxIidOnOCLL74gf/789O/fn5IlS1K3bl0+//xz/vrrL2+HahJYkiUWp8+kO/AD8CswU1V3isgQEWnlFPsSyCMie3Hd8rr5SHI9YJuIbMXVqR+qqmedfS8C44G9uFoyi5Pkgowx8ZYzZ066dOlCeHg4Bw8e5L///S9nzpwhNDSU/Pnz88QTTzBnzhyuXr3q7VBNAkib87HYRF/GeJ2qsmXLFiZPnsy0adM4efIkuXLlom3btjz11FPUr1+f9OnTeztM47CJvuJgicWY5CUqKoply5YxefJk5s2bR2RkJA899BBt2rThySefJDg4GB8fH2+HmaZZYomDJRZjkq/Lly/z/fffM3PmTL799lv+/vtv8uXLR5s2bXjqqaeoU6eOJRkvsMQSB0ssxqQMly5dYvHixcycOZPvvvuOS5cukT9/ftq2bUu7du2oXbu2JZkkYoklDpZYjEl5/v77bxYuXMjMmTNZuHAhV65cIW/evLRq1YrHH3+cRo0akTlzZm+HmWpZYomDJRZjUrbIyEgWLlzIvHnzWLhwIRcvXiRr1qw0a9aM1q1b06JFC3LlyuXtMFMVSyxxsMRiTOpx9epVwsPDmTdvHvPnz+fEiROkT5+eBg0a0Lp1a1q1akXhwoXjPpC5K0sscbDEYkzqdOPGDTZs2MC8efOYO3cuv/32GwD+/v60aNGC5s2bExgYSLp0NhWVpyyxxMESizFpw6+//sq8efP47rvvWLduHTdu3CBv3rw0a9aMFi1a0KRJE7tlFk+WWOJgicWYtOfMmTP88MMPLFy4kO+//56zZ8/i4+ND7dq1adGiBS1atKBixYrYXIGxs8QSB0ssxqRt169fZ/369SxcuJBFixaxZcsWAIoWLUqTJk1o3LgxjRo1Ik+ePHEcKe2wxBIHSyzGGHfHjh1j0aJFLF68mGXLlnH+/HlEBH9/fxo3bkzjxo3597//TaZMmbwdqtdYYomDJRZjzJ1ERUURERHBjz/+yJIlS1i3bh1RUVE88MADBAcHRycaX1/fNHXbzBJLHCyxGGPi68KFC4SHh7NkyRKWLFnCnj17AChQoAANGzakfv361K9fn1KlSqXqRGOJJQ6WWIwx9+rw4cPRSSY8PJyTJ08CrgnObiaZBg0aULJkyVSVaCyxxMESizEmIagqe/bsITw8PHq5mWgKFy4cnWjq16+f4hONJZY4WGIxxiSGuyWaggULUrt2berUqUPt2rWpWrVqippvxhJLHCyxGGOSgqqye/duwsPDWbVqFWvWrOHw4cMAZM2alRo1alC7dm1q165NrVq1yJEjh5cjvjNLLHGwxGKM8ZajR4+yZs0aVq9ezZo1a9i6dSs3btxARKhSpUp0ovn3v/9NsWLFks3tM0sscbDEYoxJLi5evMj69etZs2YNa9asYe3atURGRgKQL18+goKCqFGjBkFBQQQGBpI7d26vxGmJJQ6WWIwxyVVUVBTbt29n7dq1bNiwgQ0bNvDrr79G7y9btuwtyaZq1apJ8uJmsk0sItIU+BDwAcar6vAY+zMBkwB/4AzwlKoeFJHGwHAgI/AP0EtVlzl1woECwGXnME1U9dTd4rDEYoxJSc6fP09ERAQbNmxg/fr1rF+/nj/++AOAjBkz4ufnR1BQEP7+/vj7+1OhQoUEfzAgWSYWEfEBfgMaA0eBjUB7Vd3lVuYloIqqhopICPC4qj4lItWAk6p6XEQqAT+oaiGnTjjwhqrGO1NYYjHGpGSqyrFjx1i/fn10qyYiIiL6FlrmzJmpWrUq1atXx9/fn+rVq+Pr60vGjBnv+ZzJNbHUAgap6iPOel8AVX3HrcwPTpm1IpIe+APIq25Biqsn60+goKpetcRijDGugTV///13Nm/ezKZNm9i0aRObN2/m4sWLgKtlU6VKFdasWXNPCcaTxJKUD1EXAo64rR8FatypjKpGich5IA+uRHJTG+AXVb3qtm2CiFwHZgPDNJZsKSLdgG7gGsHUGGNSEx8fH8qXL0/58uV5+umnAdfEZ/v27YtOMidOnLivVkt8JWViie2ZuZgJ4K5lRMQXeBdo4ra/g6oeE5HsuBJLR1z9NLceRHUcMA5cLRbPQjfGmJQnXbp0lClThjJlyhASEpJ0502yM7laKEXc1gsDx+9UxrkVlhM466wXBuYCz6rqvpsVVPWY8/MiMA0ISqT4jTHGxENSJpaNQBkRKSEiGYEQYEGMMguATs7ntsAyVVURyQUsBPqq6pqbhUUkvYg85HzOADwK7Ejk6zDGGHMXSZZYVDUK6A78APwKzFTVnSIyRERaOcW+BPKIyF6gJ9DH2d4dKA30F5EtzpIPyAT8ICLbgC3AMeCLpLomY4wxt7MXJI0xxsTJk6fCkvJWmDHGmDTAEosxxpgEZYnFGGNMgrLEYowxJkGlyc57ETkNHLrH6g9x60gAyVlKihUs3sSUkmKFlBVvSooV7j3eYqqaNz4F02RiuR8iEhHfJyO8LSXFChZvYkpJsULKijclxQpJE6/dCjPGGJOgLLEYY4xJUJZYPDfO2wF4ICXFChZvYkpJsULKijclxQpJEK/1sRhjjElQ1mIxxhiToCyxGGOMSVCWWO5ARJqKyB4R2SsifWLZn0lEZjj714tI8aSPMjqWuGLtLCKn3UaG7uKNOJ1YvhKRUyIS6/QG4jLGuZZtIlI9qWOMEU9c8dYXkfNuv9sBSR2jWyxFRGS5iPwqIjtFpEcsZZLN7zee8SaL36+IZBaRDSKy1Yl1cCxlktN3QnziTbzvBVW1JcYC+AD7gJJARmArUDFGmZeAz5zPIcCMZBxrZ2Cst3+vTiz1gOrAjjvsbw4sxjWbaE1gfTKPtz7wnbd/r04sBYDqzufswG+x/LeQbH6/8Yw3Wfx+nd9XNudzBmA9UDNGmWTxneBBvIn2vWAtltgFAXtVdb+q/gNMBx6LUeYxYKLzeRbQSERim1o5scUn1mRD/397dxcqVRWGcfz/pPZF0KE81FELw5Ig8QNBKisIuugiNMhAkqzoqrCIqIhugqAIBCGCiKLogwjK4nAsM5IKoxIyK9TswvIiUdAkP0IRzLeLtY6N455zdmf2nL2N5wfC7Jl19jyzcNY7e82wVsQG8q6gHSwG3opkI9AnaWB80p2uRN7GiIg9EbE53z5M2vdoaluzxvRvybyNkPvrr3w4Kf9r/+VTU8aEsnl7xoWl2FTg95bjXZz+H/5km0ibmB0ELh6XdB1yZEVZAe7IUx+rJV1W8HhTlH09TXJdnnL4RNI1dYcByNMw80ifVFs1sn9HyAsN6V9JEyT9COwFPouIjn1b85gAlMoLPRoXXFiKFX3KaK/2ZdqMhzI51gDTI2I2sJ5/P1U1UVP6tazNpDWU5gAvAoM150HSBcAHwCMRcaj94YI/qbV/R8nbmP6NiL8jYi4wDTA0RVIAAANlSURBVFggaVZbk0b1bYm8PRsXXFiK7QJaq/c0YHenNpImAhdSz5TJqFkjYn9EHMuHrwLzxynbWJTp+8aIiEPDUw4RsRaYJGlyXXkkTSIN0u9ExIcFTRrVv6PlbVr/5hwHgC+BW9seasqYcIpOeXs5LriwFPsOuErSFZLOJn0RN9TWZgi4J99eAnwe+RuxcTZq1rY59EWkueymGgKW518vXQscjIg9dYfqRNKlw/PokhaQ3lP7a8oi4DVge0Ss6tCsMf1bJm9T+ldSv6S+fPs84Bbgl7ZmTRkTSuXt5bgwsaoT/Z9ExHFJK4BPSb+6ej0itkl6BtgUEUOkN8TbknaQPpUsbXDWhyUtAo7nrPfWkRVA0rukX/pMlrQLeJr0xSIR8TKwlvTLpR3AEeC+epImJfIuAR6QdBw4CiytazABFgJ3A1vy3DrAU8Dl0Mj+LZO3Kf07ALwpaQKpuL0XER81cUzIyuTt2bjgJV3MzKxSngozM7NKubCYmVmlXFjMzKxSLixmZlYpFxYzM6uUC4tZFyT1SXqw5XiKpNU9eq7bh1f3lfSGpCUFbfolrevF85uV5cJi1p0+0qq2AETE7og4bcCvyBPASyM1iIh9wB5JC3uUwWxULixm3XkemJH3s1gpabry3i15v4tBSWsk7ZS0QtKjkn6QtFHSRbndDEnrJH0v6StJV7c/iaSZwLGI+KPl7pskfSPpt7arl0FgWQ9fs9mIXFjMuvMk8GtEzI2IxwsenwXcRdre4FngSETMA74Fluc2rwAPRcR84DGKr0oWkhZkbDUA3ADcRipwwzYBN47t5Zh1z0u6mPXWF3mvkcOSDpJWlAXYAszOK/teD7zfsnXHOQXnGQD2td03GBEngJ8lXdJy/15gSlUvwOy/cmEx661jLbdPtByfIL3/zgIO5OXNR3KUtFpup3O3Ltl+bm5vVgtPhZl15zBpW90xyfuP7JR0J5zck35OQdPtwJUlTzsT2DrWTGbdcmEx60JE7Ae+lrRV0soxnmYZcL+kn4BtFG8tvQGYJ5Xa6vZm4OMxZjHrmlc3NjtDSHoBWBMR60dptwFYHBF/jk8ys1P5isXszPEccP5IDST1A6tcVKxOvmIxM7NK+YrFzMwq5cJiZmaVcmExM7NKubCYmVmlXFjMzKxS/wD4QGParLMfEAAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "C_0 = 0.2\n", "\n", "hypoxia = 0.17\n", "vomiting = 0.1\n", "death = 0.06\n", "\n", "t = np.arange(0, 3.5, 0.01) * ureg.hour\n", "concentration = C_0 * np.exp(-k * t)\n", "\n", "plt.axhline(hypoxia, color='purple', lw=1)\n", "plt.axhline(vomiting, color='orange', lw=1)\n", "plt.axhline(death, color='r', lw=1)\n", "\n", "plt.plot(t, concentration, 'k')\n", "\n", "plt.xlabel('time (h)')\n", "plt.ylabel('O$_2$ volume fraction')\n", "plt.text(3.2, death+.002, \"DEATH\")\n", "plt.text(2.4, vomiting+.002, \"vomiting and nausea\")\n", "plt.text(2.4, hypoxia+.002, \"first signs of hypoxia\")\n", "\n", "plt.title(f'Number of people: {N}')\n", "\n", "#plt.savefig('death.pdf')" ] }, { "cell_type": "code", "execution_count": 5, "metadata": { "ExecuteTime": { "end_time": "2018-08-24T11:01:24.222951Z", "start_time": "2018-08-24T11:01:24.220295Z" } }, "outputs": [], "source": [ "def time_of_death(N):\n", " k = N * V_prime * η / V_air\n", " return (- np.log(death / C_0) / k).to('min')\n", "\n", "def time_of_suffocation(N):\n", " k = N * V_prime * η / V_air\n", " return (- np.log(vomiting / C_0) / k).to('min')\n", "\n" ] }, { "cell_type": "code", "execution_count": 6, "metadata": { "ExecuteTime": { "end_time": "2018-08-24T11:01:24.757010Z", "start_time": "2018-08-24T11:01:24.457740Z" } }, "outputs": [ { "data": { "text/plain": [ "Text(0,0.5,'time (h)')" ] }, "execution_count": 6, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYIAAAEKCAYAAAAfGVI8AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvIxREBQAAIABJREFUeJzt3Xl8VOX1+PHPSQgJgSTs+xJEFpElCCKLCgJVquJWLQrutVrX2latSlttvz++pVbbqnVDcUetWkGr9evG4soqiyAoiwERZF8NgSzn98czEybJTDKZubMkOe/X675m5s7cc59hycm9z73niKpijDGm/kpJ9ACMMcYkliUCY4yp5ywRGGNMPWeJwBhj6jlLBMYYU89ZIjDGmHouZolARDqJyGwRWSUiK0XklxXev0VEVERaxmoMxhhjqtcghrGLgd+o6ucikgUsFpH3VPVLEekE/AjYGMP9G2OMCUPMjghUdYuqfu57vh9YBXTwvf134DbA7mYzxpgEi+URQRkRyQUGAPNF5CzgO1VdJiJVbXM1cDVARkbGwM6dO3s2ntLSUlJSvMuBFs/iWTyLl4zxvv766x2q2qraD6pqTBegCbAYOA/IBOYDOb738oGW1cXo0aOHemn27NkWz+JZPItX5+MBizSMn9MxvWpIRNKAfwPTVfU1oBvQFVgmIvlAR+BzEWkby3EYY4wJLWanhsSd95kGrFLVvwGo6hdA64DP5AODVHVHrMZhjDGmarE8IhgOXAKMEpGlvuX0GO7PGGNMBGJ2RKCqHwOhZ4PdZ3JjtX9jvFRUVMSmTZsoLCyMOlZOTg6rVq3yYFQWz+I5GRkZdOzYkbS0tIj2EZerhoyp7TZt2kRWVha5ublUdbVbOPbv309WVpZHI7N49T2eqrJz5042bdpE165dI9qHlZgwJgyFhYW0aNEi6iRgjNdEhBYtWkR1tGqJwJgwWRIwySraf5uWCIwxpp6zRGBMLbBnzx4efvjhstebN2/m/PPP93w/c+bMYf78+Z7HjcTpp5/Onj17Kn33LVu2ePrd3377bXr37k2fPn2YNGmSZ3FrE0sExsTCN9NhZi68kOIev5keVbiKPwzbt2/Pq6++Gt0Yg0imRPDf//6Xpk2bVvru7dq18/S7//a3v+Wtt95ixYoVXHXVVZ7FrU0sERjjtW+mw4KroWADoO5xwdVRJYPbb7+ddevWkZeXx6233kp+fj59+vQBYPr06ZxzzjmMGzeOrl278s9//pO//e1vDBgwgCFDhrBr1y4A1q1bx9ixYxk4cCAnnXQSq1evLreP/Px8Hn30UR566CHy8vL46KOP2LBhA6NHj6Zfv36MHj2ajRvLFwwuLS0lNzeXPXv2lK07+uij2bp1a9m2Q4cOLbft5ZdfzrXXXsspp5zCUUcdxdy5c7nyyis55phjuPzyy8vi5ObmsmPHjkrffcOGDWXf/emnn+a8885j7NixdO/endtuu61s+2nTptGjRw9GjhzJz3/+c2644Yagf7YNGzZk06ZNAFVedXP33Xdz5ZVXMnLkSI466igeeOCBsvfOOeccBg4cyLHHHstTTz1Vtr5JkyZlz1999dWy7/fKK6/Qp08f+vfvz8knnwxASUkJt956K8cffzz9+vXjscceA+DAgQOMHj2a4447jr59+/L666+HHGOk7PJRY2pq8c2we2no93fMg9JD5deVFMD8n8G6x2lUUgKpqeXfb5YHA/8RMuSUKVNYsWIFS5e6/ebn55d7f8WKFSxZsoTCwkKOPvpo/vKXv7BkyRJ+9atf8eyzz3LzzTdz9dVX8+ijj9K9e3fmz5/Pddddx6xZs8pi5Obm8otf/IK0tLSyUyTjxo3j0ksv5bLLLuPJJ5/kpptuYubMmWXbpKSkcPbZZzNjxgyuuOIK5s+fT25uLm3atCnb9rzzzuOVV14pt+3u3buZNWsWb7zxBuPGjeOTTz7hiSee4Pjjj2fp0qXk5eWF/O4rVqwo992XLl3KkiVLSE9Pp2fPntx4442kpqbyP//zP3z++edkZWUxatQo+vfvX+nPtbS0lJ49e3LllVfy7rvvVnv55erVq5k9ezb79++nZ8+eXHvttaSlpfHkk0/SvHlzDh48yMCBA5k4cSItWrQIGedPf/oT77zzDh06dChLotOmTSMnJ4eFCxdy6NAhhg8fzqmnnkrTpk2ZMWMG2dnZ7NixgyFDhnDWWWd5evFC/ToimD4dcnMZMWoU5Oa618Z4rWISqG69B0455RSysrJo1aoVOTk5jBs3DoC+ffuSn5/PgQMH+PTTT7ngggvIy8vjmmuuYcuWLdXG/eyzz5gwYQIAl1xyCR9//HGlz4wfP55//etfALz00kuMHz++2m3HjRuHiNC3b1/atGlD3759SUlJ4dhjj62U5KozevRocnJyyMjIoHfv3mzYsIEFCxYwYsQImjdvTlpaGhdccEHQbR988EF69erFI488wrhx49i+fTsLFiwI+fkzzjiD9PR0WrZsSevWrdm6dSsADzzwAP3792fIkCF89913rFmzpsoxDx8+nMsvv5zHH3+ckpISAN59912effZZ8vLyOOGEE9i5cydr1qxBVbnzzjvp168fY8aM4bvvvivbr1diWWuoE/As0BYoBaaq6v0i8ldgHHAYWAdcoap7QkfyyPTpcPXVUFDgbnfesMG9Bpg4Mea7N3VIFb+5A25OoGBD5fWZXWDMHA56fMMRQHp6etnzlJSUstcpKSkUFxdTWlpK06ZNy36rjlSw30KHDh3K2rVr2b59OzNnzuR3v/tdtdsGjq/i2IuLi2s0psDtU1NTKS4u9lc+rtY777zDDTfcwJgxY/jDH/7AGWecweDBg8uSWTj7mjNnDu+//z6fffYZmZmZnHTSSWXX9Ad+58Dr/B999FHmz5/PW2+9RV5eHkuXLkVVefDBBznttNPK7fPRRx9l+/btLF68mLS0NHJzcz25wz1QLI8I/B3KjgGGANeLSG/gPaCPqvYDvgbuiOEYjpg0CQoKyq8rKHDrjfFS/8mQmll+XWqmWx+hrKws9u/fH/H22dnZdO3alVdeeQVwd6MuW7as2v0MGzaMl156CXBzESeeeGKlbUSEc889l1//+tccc8wxZadEwtk2HJF898GDBzN37lx2795NcXEx//73v4N+bsCAAfzrX/+itLSUn/70p3Tv3p0XXniBM844I+x97d27l2bNmpGZmcnq1atZuHBh2Xtt2rRh1apVlJaWMmPGjLL169at44QTTuBPf/oTLVu25Ntvv+W0007jkUceoaioCICvv/6aH374gb1799K6dWvS0tKYPXs2GzYE+SUjSnHvUKaq76qqP+XPw5Wijr0Kk1zVrjcmUl0nwuCp7ggAcY+Dp7r1EWrRogXDhw+nT58+3HrrrRHFmD59OtOmTaN///4ce+yxQScdx40bx5tvvlk2WfzAAw/w1FNP0a9fP5577jnuv//+oLHHjx/P888/X+43af+2Q4cOrXLb6kTy3Tt06MCdd97JCSecwJgxY+jduzc5OTmVPjdp0iRUlT59+jBw4EDatGnDNddcw4QJEygtLQ1rX2PHjqW4uJh+/frx+9//nuOPP77svSlTpnDmmWcyatQo2rVrV7b+1ltvpW/fvvTp04eTTz6Z/v37c9VVV9G7d2+OO+44+vTpwzXXXENxcTHjx49n0aJFDBo0iOnTp9OrV6+wxlUTEu4hVFQ7cR3KPsQdCewLWP8f4F+q+nyQbco6lLVq1Wrgyy+/HNUYhlx4IRlBzqsVtmnDPN9vLZE6cOBAuasDomXxki9ehw4dOProoz2JV1JSQmrFyWKL53k8/7+D4uJiJkyYwCWXXFI2d5IM4/M63tq1a9m7d2+5daeccspiVR1U7cbhdK+JZiGgQ1mF9ZOAGfiSUVWLJx3Knn9eNTNTFY4smZlufZSSqSORxYtNvC+//NKzePv27fMslsUL7Te/+Y32799fe/bsqTfeeKOWlpZGFS9ciYoX7N8oYXYoi+nlo0E6lPnXXwacCYz2DTb2/BPCt90GmzdD8+bwwAM2UWxMHXXvvfcmegi1RszmCIJ1KPOtHwv8FjhLVQtCbR8TEyfCt99S3Lgx/PSnlgSMMYbY3lDm71D2hYj4r1m7E3gASAfe811aNU9VfxHDcZSXksL+Hj1otnhx3HZpjDHJLBEdyv4bq32Ga3+PHjSbMQMOH4aGDRM9HGOMSaj6dWexz/6ePV0SWLky0UMxxpiEq5+JoEcP92TRosQOxJgwxasMdSiHDh1izJgx5OXllZWTCCawGF5NVfyOc+bM4cwzz4wolqmZepkICtu3h6ZNLRGY2PHVtSIlxZO6VvEqQx3KkiVLKCoqYunSpSHLL0Sr4nc08VMvEwEiMHAg2ISxiQV/XasNG9wdK/66VlEkg3iUoQbYtWsXF110Ef369WPIkCEsX76cbdu2cfHFF5dVBV23bl25bRYvXkz//v0ZOnQoDz30UNl6f1nlESNGhFVWueJ39H/2/PPPp1evXkycODHsGkKmZupvGepBg+Bvf4NDhyCgkJQx1br5ZqiqeNu8ee7fVaCCAvjZz+DxEGWo8/LgH4ktQw1w11130a9fP958801mzZrFpZdeytKlS3niiSe49957efPNNyuN7YorruDBBx9kxIgR5UpA+Msqz507l4YNG5aVVe7UqVPQssoVv+OcOXNYsmQJK1eupH379gwfPpxPPvkkaDlpE536nQiKimDFCnd0YIxXKiaB6tZ7wF+GOisrq1IZ6uXLl5crQ31kOJXH8/HHH/PMM88AMGrUKHbu3FmpbEGgvXv3smfPHkaMGAG4ctNvv/024MoqL1++nJdffpmUlBT27t3LmjVr6NixI3feeScffvghKSkpVZZVHjx4MB07unJkeXl55OfnWyKIgfqbCPw//BctskRgaqaK39wBNycQrEJkly4wJ7nLUAc79VJVAxRVDfm++soqDxs2rNz3ffrpp8Muqxys7LPxXv2cIwD3n7V5c5swNt6bPBkyK5Shzsx06yMUrzLUJ598Mv4Cj3PmzKFly5ZkZ2eHjNu0aVNycnLKms5MD5gHqWlZ5Wi/o4lc/U0EIu70kE0YG69NnAhTp7ojABH3OHVqVCVN4lWG+u6772bJkiX069eP22+/vew0UVWeeuoprr/+eoYOHUqjRo3K1vvLKp900knlyipPnDgxaFllL76jiVA4lekiWYBOwGxcH4KVwC9965vjmtOs8T02qy6WJ9VHA5RVp7zjDtUGDVQPHvQmnkcsXvLFs+qjFi/Z40VTfTQRHcpuBz5Q1e7AB77XiTFoEBQXw/LlCRuCMcYkWtw7lAFnA/7jzWeAc2I1hmr5J4nt9JAxph6Le4cyYKOqNg14b7eqNguyjacdygKVdbBSZdi557Jz2DC+uu226ON5PT6LlzTxOnToQLdu3aq8giZcdaUjlsVLnniqyrp162pPhzJgT4X3d1cXI2ZzBKqqY8eq9u/vXTwPWLzki7d+/Xrdvn17yC5XNVFXzklbvOSIV1paqtu3b9f169dXeo8k7lC2VUTaqeoWEWkHbIvlGKo1cCBMmQIHD0LAFQ/GBOrYsSObNm1i+/btUccqLCwkIyPDg1FZPIvnZGRklN14F4mYJYJQHcqAN4DLgCm+x8rXsMXToEFQUgLLlsGQIQkdikleaWlpdO3a1ZNYc+bMYcCAAZ7EsngWzwuxvGrI36FslIgs9S2n4xLAj0RkDfAj3+vEsQljY0w9l4gOZQCjY7XfGuvYEVq3tjuMjTH1Vv29s9jPf4exJQJjTD1liQDc6aEvv3Slgo0xpp6xRADuiKC0tOoa88YYU0dZIgCXCMAmjI0x9ZIlAoD27aFtW5snMMbUS5YI/GzC2BhTT1ki8Bs4EFavhgMHEj0SY4yJK0sEfjZhbIyppywR+AX2MDbGmHokZolARJ4UkW0isiJgXZ6IzPOVm1gkIoNjtf8aa9fOTRrblUPGmHomlkcETwNjK6y7B/ijquYBf/C9Th42YWyMqYdi2aHsQ2BXxdVAtu95DrA5VvuPyMCB8NVXsH9/okdijDFxE9MOZb7OZG+qah/f62OAd3DF6FKAYaq6IcS2se9QVkHzefPod8cdLPnHP9jbv3/U8bwen8WzeBbP4tVEUnQoA3KBFQGvHwB+4nv+U+D9cOLEtENZoO+/VwXV++7zJl6ELJ7Fs3gWz4t4hNmhLN5XDV0G+DuVvQIkz2QxQJs2riy1TRgbY+qReCeCzcAI3/NRwJo47796NmFsjKlnYtmq8kVgJNBSRDYBdwE/B+4XkQZAIb45gKQyaBDMnAn79kF2dvWfN8aYWi6WHcouCvHWwFjt0xP+G8s+/xxGjkzoUIwxJh7szuKK7A5jY0w9Y4mgolatoHNnmzA2xtQblgiCsQljY0w9YokgmEGDYO1a2LMn0SMxxpiYs0QQTOCEsTHG1HGWCIKxCWNjTD1iiSCYFi0gN9cmjI0x9YIlglBswtgYU09YIghl0CBYvx52VaykbYwxdUtcO5T51t8oIl+JyEoRSa7GNIFswtgYU0/EtUOZiJwCnA30U9VjgXtjuP/o2ISxMaaeiHeHsmuBKap6yPeZbbHaf9SaNYNu3SwRGGPqvHh3KFsKvI47UigEblHVhSG2jXuHsop6//GPZK1ezfwXX/QkXrgsnsWzeBbPi3jJ2qFsBa5LmeCa0nyDLxlVtcStQ1lF99zjOpbt2OFNvDBZPItn8SyeF/FI0g5lm4DXfGNcAJQCLeM8hvD55wnsfgJjTB0W70QwE9eZDBHpATQEdsR5DOE77jj3aPMExpg6LN4dyp4EnvRdUnoYuMx3+JKcmjaF7t0tERhj6rREdCi7OFb7jImBA+HTTxM9CmOMiRm7s7g6gwbBxo2wLXmvdDXGmGhYIqiOTRgbY+o4SwTV8U8YWyIwxtRRlgiqk50NPXvahLExps6yRBCOgQPtiMAYU2dZIgjHoEGwaRN8/32iR2KMMZ6zRBAOmzA2xtRhlgjCMWAAiFgiMMbUSZYIwpGVBb162YSxMaZOinuHMt97t4iIikjyFpyraOBASwTGmDoprh3KAESkE/AjYGMM9+29QYNgyxbYvDnRIzHGGE/Fu0MZwN+B24DkLTYXzCBfbwebJzDG1DHx7lB2FjBaVX8pIvnAIFUNWoY6GTqUBUo5eJCTzjyTDZdcQv7ll0cdz+vxWTyLZ/EsXkVJ16EMyATmAzm+1/lAy3DiJKxDWUXHHqt6xhnexQvB4lk8i2fxvIiHVx3KRGSoiDwkIstFZLuIbBSR/4rI9SKSU4Pk1A3oCizzHQ10BD4XkbY1iJFY/gnjJG6hYIwxNVVlIhCRt4GrgHdwE7/tgN7A74AM4HXf6Z5qqeoXqtpaVXNVNRfXtvI4Va09t+sOGgRbt9qEsTGmTqmuMc0lWvkc/gHgc99yX6hLQIN1KFPVaVGON7H8E8aLFkGHDokdizHGeKTKRBAkCYT9GQ3docz/fm51sZNO//6QkuKuHDr77ESPxhhjPBHW5aMicp6IrBGRvSKyT0T2i8i+WA8u6WRmwrHH2o1lxpg6Jdz7CO4BzlLVHFXNVtUsVc2O5cCSlk0YG2PqmHATwVZVXRXTkcTDN9NhZi4jNo+CmbnudU0NGgTbt7uy1MYYUwdUOUcgIuf5ni4SkX8BM4FD/vdV9bUYjs1b30yHBVdDSQECULDBvQboOjH8OIETxp06eT1KY4yJu+qOCMb5lmygADg1YN2ZsR2ax5ZNgpKC8utKCtz6mujXD1JTbZ7AGFNnVHfV0BXxGkjMFYSocRdqfSiNGkGfPlZzyBhTZ1R3Q9nvRKRZFe+PEpHacWSQ2blm66syaJBNGBtj6ozqTg19AbwpIh+IyF9F5DYR+YOIPCciX+BOEc2P/TA90H8ypGaWX5ea6dbX1MCBsHMnbKxdlbSNMSaY6k4NvY4rI9EdGI4rMbEPeB64WlUPxn6IHvFPCC+bhBZsQCQVBj9Ws4liv8AJ4y5dvBujMcYkQFiXj6rqGlV9WlX/rKr/UNV3qksCwTqU+Y4qVvsK2M0QkabRfoEa6ToRzslnddM7QEugSdfI4vTtCw0a2ISxMaZOiHeHsveAPqraD/gauCOG+w9pR8aJkJoB+S9GFiAjwyUDmzA2xtQBce1Qpqrvqmqx7+U8XCnquCtJyYQO42Djy1BaFFkQmzA2xtQRce1QVuG9/wD/UtXnQ2wb0w5lualL6bP79yxv/hd2ZQyucYx2b7xBz7//nXnTp7MjOztpOhJZPItn8Syen6cdyoAewAcc6TbWD/hdGNvl+repsH4SMANfIqpuiUmHsuJC1ZdzVD+9NLIgixapgurLLydVRyKLZ/EsnsXzw6sOZT6P487nF/mSx3LgwrDTUgARuQx3V/JE30ATIzUdOv0Evn0NiiO4+KlPH0hLswljY0ytF24iyFTVBRXWFQf9ZBVEZCzwW1wl04LqPh9zuROg+ABsfqvm26anu3ITNmFsjKnlwk0EO0SkG6AAInI+sKWqDXwdyj4DeorIJhH5GfBPIAt4T0SWisijkQ/dA61HQkZbyH8hsu1twtgYUwdU16rS73pgKtBLRL4DvgEurmoDDd6hLLlaVaakQpfxsOYROLwHGtbwtobiYti7lxGjRrkbyyZPhokR3KBmjDEJFO4NZetVdQzQCuilqieqan5MRxYvXSZA6WH4dkbNtps+3S3gylpv2ABXX122zhhjaouwjgh8dwBfirsKqIGIAKCqN8VsZPHS4nho0g02vADdalBsddIkKCwsv66gwK23owJjTC0S7qmh/+JuAPsCKI3dcBJABLpcBF/+Lxz8Hhq1DW+7UAXnrBCdMaaWCTcRZKjqr2M6kkTKnQAr/5+707hnmAc5nTu700HB1htjTC0S7lVDz4nIz0WknYg09y8xHVk85RwDTfvX7OqhyZMhs0JZ64wMt94YY2qRcBPBYeCvuMtBF/uWunUnVe4E2Dkf9q8L7/MTJ8LUqdClCyoCKSnQsSNcFOxiKWOMSV7hJoJfA0eraq6qdvUtR8VyYHHXxXej9IaXwt9m4kTIz2furFkwbRqsXQtPPx2T4RljTKyEmwhW4prX112NO0OrEyF/emQ3iF16KQwfDr/9LezaVf3njTEmSYSbCEqApSLymIg84F9iObCEyJ0A+1bBnuU13zYlBR5+GHbvdpeQGmNMLRFuIpgJTAY+5cgcQZVFdkJ0KGsuIu+JyBrfY7NIBx4TnS4AaQAbImxY068f3HgjPPYYLFzo7diMMSZGwr2z+JlgSzWbPU3lDmW3Ax+oandcWevbazziWMpoCW1/5DqXaYS3S9x9N7RpA9ddByUlng7PGGNiocpEICIv+x6/8PUZLrdUta0G6VAGnA34E8gzwDkRjjt2cidAwUbY/mlk2+fkwH33uWJ0jz/u7diMMSYGquxQJiLtVHWLiHQJ9r6qBrmjqtz2uQR0KBORParaNOD93aoa9PRQrDuUher4k1pawLCt5/F9o7GsaXpzZPFU6f+b39Bk7VoWPPssRU1rVswumTocWTyLZ/FqbzyvO5T9JZx1QT6TS0CHMmBPhfd3h7P/mHQoq8pHP1V9taVqyeHI461cqdqggeqVV3o/Potn8SyexQsDHnco+1GQdT8Oc9tAW0WkHbijDWBbBDFiL3cCHNoB378feYzeveHXv4Ynn4RPIzzNZIwxcVDdHMG1IvIFrrlM4PzAN0AE11jyBnCZ7/llwOsRxIi9dmMhrambNI7G73/v7ja+7jrXu8AYY5JQdUcELwDjcD/AxwUsA1W1ysY0ITqUTQF+JCJrcEcZU6Icf2ykpkPnn8CmGVAcxX10TZrA3/8Oy5a5ewyMMSYJVVl9VFX3AnuBGhfQ0eAdygBG1zRWQnSZAOumuX7GnS+IPM5PfgKnnuqODn76U2gbZplrY4yJk3DnCOqf1iOgUbvI+xn7icA//+ma2Nx6qzdjM8YYD1kiCCUlFTqPh83/df2Mo9G9O9x2Gzz/PMyd6834jDHGI5YIqpLr72f8WvSx7rgDcnPh+uuhqCj6eMYY4xFLBFVpPgiaHB396SFwTWweeABWroT7748+njHGeMQSQVVEIPci2DoLDm6JPt64cW65+27YtCn6eMYY4wFLBNXpchGgsMGjEhf33++K0f3mN97EM8aYKFkiqE7OMdBsAGzw4PQQQNeucOed8PLL8H4Udy4bY4xHLBGEo8tFsHMB7F/rTbxbb4Wjj3YTx4cOeRPTGGMiZIkgHJH0M65KRoa7t+Drr13JamOMSaCEJAIR+ZWIrBSRFSLyoohkJGIcYWvcCVqfHHk/42BOO83ddfz//h/k53sT0xhjIhD3RCAiHYCbgEHq+hSkAhfGexw11uUi2Lca9izzLubf/+6uTPrVr7yLaYwxNZSoU0MNgEYi0gDIBDYnaBzh63S+62ccbUXScjE7wR/+ADNnwn//611cY4ypgSo7lMVspyK/BCYDB4F3VXVikM8kpENZVfruvIPGxeuZ1/pFkCM5NJoOQlJUxKCf/5yUoiIWPvkkpenpSdXhyOJZPItXe+N52qHMywVoBswCWgFpwEzg4qq2iXuHslDWP686HdWtH3oTz2/WLFVQvesub+JVYPEsnsWrn/HwuEOZl8YA36jqdlUtAl4DhiVgHDXX8WxIbQQbPDw9BHDKKXDRRW7iuEMHRowa5eoSTZ/u7X6MMSaIRCSCjcAQEckUEcH1J1iVgHHUXFoTlww2vgylHheOGz7c3XG8eTOiChs2wNVXWzIwxsRc3BOBqs4HXgU+B77wjWFqvMcRsS4XwaGdsOU9b+P+9a+V1xUUwKRJ3u7HGGMqSMhVQ6p6l6r2UtU+qnqJqtae22vbjYWGzbw/PbRxY83WG2OMR+zO4ppKbQidPOhnXFHnzsHXd+rk3T6MMSYISwSRyJ0AxT/Ad296F3PyZNezoKKmTWHvXu/2Y4wxFVgiiESrk6FRe+8qkgJMnAhTp0KXLqgIdOkCV14JX34JQ4bAmjXe7csYYwJYIohEuX7Gu72LO3Ei5Oczd9YsV39o2jRXqnr7djjhBPjgA+/2ZYwxPpYIIpU7wV1C6kU/46qMGAELF0L79q4JZ7SJAAAb4ElEQVRQ3cMPx3Z/xph6xxJBpJoPhPQ2sPB6RmweBTNz4ZsYXfPftSt8+imcfrrrYXDttVDk8X0Mxph6yxJBpPJfgMO7oPQQgkLBBlhwdeySQXY2zJgBt98Ojz4Kp54KO3fGZl/GmHrFEkGklk0CrfBbeUmBWx8rqanw5z/Dc8/BZ5/B4MFuMtkYY6KQqMY0TUXkVRFZLSKrRGRoIsYRlYIQN3qFWu+liy+GuXPdncdDhsBbb8V+n8aYOitRRwT3A/+nqr2A/tSWWkOBMkPcABZqvddOOMFNInfvDuPGuRIVCSgpboyp/RLRoSwbOBmYBqCqh1V1T7zHEbX+kyE1yA1gHc6M3xg6doSPPoILLoDbboPLL4fCwvjt3xhTJyTiiOAoYDvwlIgsEZEnRKRxAsYRna4TYfBUyOyCIpDZCbJ6wpqHYd20+I0jMxNeegn+9Cd49llX0vr77+O3f2NMrRf3DmUiMgiYBwxX1fkicj+wT1V/X+FzSdehrLp4KaWFHLv7LlocWsCa7Ov5rsn5cR1fy7lzOWbKFIqyslgxeTIHunePKp7X47N4Fs/ixTdeMncoawvkB7w+CXirqm2SpkNZOPGKC1Xnnus6ma2YHH28mlqyRLVTJ9VGjVRvukm1SxctFVHt0kX1+ecjj+vV+CyexbN4cYtHsnYoU9XvgW9FpKdv1Wig7lwDmZoOJ74MuRe7S0mX3hHfSdy8PDeJ3KEDPPAAbNhgjW6MMVVqkKD93ghMF5GGwHrgigSNIzZSGsDQZ6BBY/hyChQfgIH3l2t4H1Nt2sDhw5XX+xvdTJwYn3EYY2qFhCQCVV0KVH/eqjaTFDj+EWjQBFbf55LB4Cdcwbp4+Pbb4Out0Y0xpgK7sziWRGDAX6Hv3bD+afh0ApQE+U09FkI1ulGFkSNh9my778AYA1giiD0R6HsXDLjXNb3/6CdQEodr/YM1umnUCC65BL7+GkaNcpVNP/jAEoIx9Zwlgng55jfuVNHmt2DOGVB0ILb7C9bo5vHH3b0G69fDgw+6xzFj4KST4L33LCEYU09ZIoin7r9wk8jb5sDsU+FwjG+ortjoxj9JnJEBN9wAa9fCQw+5K4pOPRWGD4d33rGEYEw9Y4kg3rpeAie+ArsWwQejoHB74saSkQHXXecSwiOPwKZNMHYsDB0Kb79tCcGYesISQSJ0Og9OfgP2rYL3R0DB5sSOJz0dfvELlxAee8yVqDj9dFfY7q23LCEYU8dZIkiU9mNh5P9Bwbfw/klwID/RI4KGDd1NZ19/7eYTtm+HM890fQ/+8x93M1puLiNGjYLcXLs5zZg6whJBIrUZAaPeh0O7XDL48j6YmRv71pfVadgQrrrKJYRp02DXLjjrLHfFkd2pbEydY4kg0VqeAGPmwKG9sPQWKNgQn9aX4UhLgyuvhNWroUWLyqeICgrgzjsTMzZjjGcSlghEJNVXhvrNRI0haTTrDw2zKq+PdevLcKWluaOCYDZuhFtugfnzbS7BmFoqkUcEv6Q2diaLlYNbgq+PR+vLcIS6U7lRI1fcbsgQd6/Cr3/t+imXlsZ3fMaYiCWqZ3FH4AzgiUTsPymFanGZ0gA2/hs0wT9Yg92pnJnpJpW3bYNnnnGVTx96CIYNc0nh5pvhk08sKRiT5BJ1RPAP4DbAfkL4BWt9mdIQ0prBx+fDW8fC+megtCgx4wt2p/LUqW5906Zw6aXwxhsuKTz3HAwcCI8+CieeCJ06wU03ubaagUnBrkIyJikkokPZmcDpqnqdiIwEblHVSo1+a2OHsmi1Lnifo/Y/QXrJNg6ltmZ91lVsa3QKrQrn0mX/CzQpXkdhahs2NhnP95mnUyrpcR1fTeOl/vADLT77jFYffkjz+fNJPXyYQy1asOOkkyhq3JhOr75K6qFDZZ8vSU/nq1tuYduYMXEZn8WzeHU9XjJ3KPszsAnIB74HCoDnq9qmVnUoi1W80lLVTW+pvjPMdT/7d2vVFX9WPbQnOcZXnX37VF98UfUnP3Hd09zUcuWlS5fEjM/iWbw6GI8k7lB2h6p2VNVc4EJglqpeHO9x1Doi0OF0+NHHMGYuNBsAy+6A17vAst8ltlRFOLKy4MIL4dVX3ekjkeCf27ABXnghdD8FY4zn7D6C2kYEWp8Mp/wfjF0EbX8EK//XJYTFN8MPteAHaJMmoa9CEnHzDp07u3mIiy92cw0rV9qkszExktBEoKpzNMj8gAlT84Fw0itwxpfQZTx8/RD8pxvMvwr2rXE3oyXDncrBhLoK6Zln4PPP4f77Xa2jDz6Aa6+FPn2gZUsYNw7uuQc+/bRyO06bfDYmIonqWWy8lNMLhjzlOqGtuhfWPQHrpoGkgpYgcOROZYCuSdCz2F8Se9IkdONGpHNnlxz86wcMcFcaqbq+CR9/7K46+vhjeNN3D2JGhquDdNJJLin8859w8KD7vv4SGIH7MsYEZaeG6pLGXWDQg3BWPjTIBi0p/35JASxLopIQofolBBKBbt3gssvgiSdcuYutW+G119yRwsGDMGUK/PWv7nmgggK4/fa4fBVjajM7IqiLGrWB4v3B3yvYCJ9dBu1+DO1OhfTm8R2bF1q3hnPPdQvAgQOQnR28xMWmTdCxI/Tt65Z+/dxjr16u/LYxxhJBnZXZ2Z0Oqig1E757E755FiQFWpzgkkL7H0Pz49y62sY/+bwhyPdt2tT1Z/7iC5g168i8QoMG0LNn5QTRufORK5qmT4dJkxixcaNbH3jqypg6xBJBXdV/spsTKCk4si41EwZPhS4Xwq6FsPltt3xxF3zxB8hoDW1Pc0mh3amQ3iJx46+pyZPdnEBBwPfNzHTzBv4f3kVFsGaNSwrLl7vHefPgpZeObJOd7RJCRgZ8+CEUFdmcg6nzLBHUVf4J4WWT0IKNSGZnlxz861sOcUu/P7p7ELa840sMb0H+c8GPFvJfhGWTGFGwEWZWiJdo1U0+g6ui2ru3W8aPP7J+3z5YsaJ8gpg1K3jZ7auvdp/t1s0tRx8NHTpAShhHUnaEYZKUJYK6rOtE6DqRuXPmMHLkyNCfy2gFXS92S2lJ8KOFBlnu6CJZr0IC90N1Yhjft6LsbFcob9iwI+tC/WAvKID77nNHF37p6dC165HEEJgkcnNdo5/p08uOWOwIwyQbSwSmvJTU4EcLC64JfhXSwutchdRmeZDVvXbOMQQTas6hSxdYt87d+bx2rXvuX9auhTlz4Icfjnw+JcUV3fv+ewioqwQcaexjicAkmCUCUzX/0cJnlwZ/v3gffHKhe96gMTTt78pfNMuD5gMgpw+khrg655vpyXuqKdScw+TJkJrqftPPzYWKBfJUXQkNf2LwJ4lQN7dt3Ajt27tk0bHjkSXwdfv27qiiIjvVZDwS90QgIp2AZ4G2uDLUU1X1/niPw9RQqKuQMjvDiNdh1xLY7Vu+eRbWPOTelwaQ09slhmYDjiSJ794sm8xOylNN4cw5BCMCbdq4JfBU08cfBz/CyMmBH//YXea6ejW89x7sr3Dprz9mYKLYvh1mzIDDh7071WSJpd5KxBFBMfAbVf1cRLKAxSLynqp+mYCxmHCFugqp///6fsjnAVe49VoKB9b7EsNSlyS2vOsShJ/vrudy/K05kyERQORzDsGEOsJ46KHKP2z37XOJYdMmdwrK/3zTJneUMXs27N1beR8FBe7GuwceOJKM/EvbtuVf5+SUL/xncxj1WtwTgapuAbb4nu8XkVVAB8ASQTKr7iqkQJICWUe7pfMFR9Yf/N4lht1LQt/hXLAB3hkKWd2gyVHQpJtbsrpBRtvQVUshuU811eQIIzv7yNVNoaSkBL+BrqTE3TuxcSMsXOhOUwUr1tewYfnEMGdO+SQF7vVvfwtnneXu1ajqzz4YO8KoNeLemKbczkVygQ+BPqq6r8J79a4xTX2KN2TrhWSUbK20vlgasT+tF41KNpNesh0JaGJXIukUprbjYIP2HExtz8EG7Sn0PWYdWknPff8gVQ+V+/xXObewLbPuNboZcuGFZGyt/OdX2KYN8wLviygpIW3fPhru3k3D3btJ27Wr7HnD3btJ8z02WbOGqn7MlzZoQFF2NkU5ORRnZ7vnvtfB1uUsWUL3hx/2tPFQ6/ff56gnniB92zYOtW7N+quuirqJESTn369X8cJtTJOwRCAiTYC5wGRVfa2qz/bs2VO/+uorz/Y9x4tDfYsXXZBvpoe+4c3/W3zJYfgh351mOrAO9q9zjwfWuXUlB4OGLie9FYz4DzRqD43aQkpazcYYzhFQDXny5xdwKqdMZuaR9qE1lZsbfA6jRQtXr2nnTtixwz36F//r4uLw95OZeaS9adOm0KzZkecVX2dkxO77+mPWdA4oDMn0/01EwkoECblqSETSgH8D06tLAqaOCudUU2pDyO7hlopUofD7I8lh3uXB93NoO7w7xPdC3FVQjdpDow6+x/aQ2f7I80btXfLY8FLdnMwOJdQcxv33Vx1T1U1uV0wUF4foNVVQAP/5D+zZA4WFVY8pI+NIUli/vnLZ8YICuPFGl4iys13zo+zs8kvjxsFPacViTqQWnwpLxFVDAkwDVqnq3+K9f5NEwr3hLRgRaNTOLa1PhOV3Bb+qKaMtnPA4HNwMBZvdo3/ZtQgKtwEVjool1T0Gm8xe/EtXeiOjFaS3dEuDxuGN2es5DC8ns6O5Ssr/Q7dr1yPrJ00KfR9Gfr57XljoEsKePbB795HnFV/v3u2uqApm9264/PLQ40tJCZ4gPvoo+JzIL3/pEmCTJkeWrKwjz9PT62RiScQRwXDgEuALEVnqW3enqv43AWMxdUWoq5oG3Asdquh9VFoEhVuh4LvySWLl/wb//OGdMOfH5delNjqSFNIDEkRgsti1FFbfB6WF3h1hJHNiqeo+DL+MDHc1U9u21ccLdeqqY0eYO9ddaRVs2b+/8rrduysnAb+dO+G880KPIzW1fJLwLwsWBC+DfuON7gbDzEx3dNK48ZHnFR8bNixf8DCOV3El4qqhj6HKeSljaq4mVzUFSkmDzI5uCfTN9OBHGI3aw4mvuFNOh3a4pTDg+aEdcGCteyzaV3n7QCUFMO8KWPsIpDWDhk3dktYUGlbxukF2bE5deZlY4nXqasoUOOqomscLlVjat4e33nKlzWuyVEwCfrt3wzXXhDem1NQjSWHHjspzLwUF7kirLiQCY2ImmlNNFYU6wsi7B1oNC71doJLDR5LD23lUOgUFoEWQku6OQvauhKI9cHhP8M+W8f8eVeEzJQWw8Bew+3NIy3b1odKyKjyv8DrFd6ojYPLes8QyDPgHUABkAv0jCwPEL7Hccw/k5dU8XlVHLPPnu6OCH35w+wt8DLXu8ceD72fjxpqPLQyWCIwJJtIjjECpDd1EdGb7Ku7M7gKjPyi/Tkuh+AAc3u2SwuE9vgQR8HrFH4Pvs/gArHm0fAKrijRwiaFob4haUr+AXYvdPEhaE0ht7J6HWvzvf/ta8ieW7Z/A/05Ft5cgrVLhzstic8TSvn3N4737bvDE0rlzZOOrhiUCY0KJxxFG/8mVPysp7odzWrZrPxrM+qdDJ5Zz8l0V2eIDrlNd0T4o2l/+edF+VyfK/9xfEqSi4gOwbioU/xD8/ZooKYD5P3NlzlMbue/fwPeY2gga+B7LPfc9bv/Y9eMuN8fycyg9DF0vdcUSa+Kb6dDuGfiHr5ouJZD6DHwzPLJEFYvE8vMr4WDAlVKNGpafY/GQJQJj4sGLI4xA1SWWlFRomOOWcHz3ZtWJRdXdt1H8Q4XlwJHnJQHrQ905XnrIHdkUf+filRyE4gL3WHoo+DahlByE+Ve6JSUNUjIgtcISat2mGZWPmkoKYPGN7julprvTZv7HlIaV15W9lw6bXvc2sQwDfqbwErADaAlcqG59DFgiMCZevDzCiHdiEXG/pTfIBFpVH2/NY6ETy2nzg29TWgKlhVB80I0jMEm8dyIh5036/hFKCt1SWnjkeeC64gNursa/rvhA8FiHd7vTYV4oKYB5l8LySS6R+BNKSkN32jCl4pJ+ZP03z8PQIhgaGLAoZrW4LBEYU1vV5sQSTEoqpDQOfl9GVXMsff9Q8/HNzA1xVVhHl6hKD0HJofKPldYdPvJ88U3B96Ol0HqE77OH3QUE/ufFBVC6O8R7+4PHK7DJYmNMLNW1xBJJvLwpbnK/plbdFzpRDX2m5vFCJarM2EwW15F2UsaYpNN1IpyTz9z2s9w8QzSnNLpOdHWoMrugiPsBG1iXKtHx+k92iSRQtInKy3jVsERgjKkdvEwsXsdL9kRVjYQkAhEZKyJfichaEbk9EWMwxhhPJXOiqkbcE4GIpAIPAT8GegMXiUgVHTiMMcbEUiKOCAYDa1V1vaoexl0pe3YCxmGMMYYENKYRkfOBsap6le/1JcAJqnpDhc9ZhzKLZ/EsnsWLQ4cyVDWuC3AB8ETA60uAB6vapkePHuql2bNnWzyLZ/EsXp2PByzSMH4uJ+LU0CagU8DrjsDmBIzDGGMMiTk11AD4GhgNfAcsBCao6soqttkPeNe02FXu2GHxLJ7Fs3h1PF4XVa22JkgiGtMUi8gNwDtAKvBkVUnA5ysN5zxXmERkkcWzeBbP4tXHeMEkpMSEuraU1prSGGOSgN1ZbIwx9VxtSQRTLZ7Fs3gWz+LFRtwni40xxiSX2nJEYIwxJkYsERhjTD2X1InA6yqlIvKkiGwTkRUexOokIrNFZJWIrBSRX3oQM0NEFojIMl/MP3oQM1VElojIm9HG8sXLF5EvRGSpiCyKMlZTEXlVRFb7/hyHVr9VyFg9fWPyL/tE5OYox/cr39/DChF5UUQyooz3S1+slZGOLdi/YRFpLiLvicga32OzKONd4BtjqYjU6LLFEPH+6vs7Xi4iM0SkaZTx/scXa6mIvCsiYXeSqepngIjcIiIqIi2jHN/dIvJdwL/F06Mdn4jc6PtZuFJE7gk3XtjCuf04EQvuHoN1wFFAQ2AZ0DvKmCcDxwErPBhfO+A43/Ms3E1y0Y5PgCa+52nAfGBIlDF/DbwAvOnR30s+0NKjWM8AV/meNwSaevhv53vczTSRxugAfAM08r1+Gbg8inh9gBVAJu6y7feB7hHEqfRvGLgHuN33/HbgL1HGOwboCcwBBnkwvlOBBr7nf/FgfNkBz28CHo0mnm99J9y9TRtq8u87xPjuBm6J8N9JsHin+P69pPtet47032GoJZmPCDyvUqqqHwK7vBicqm5R1c99z/cDq3A/PKKJqarq76qd5lsins0XkY7AGcAT0YwrFkQkG/ePfhqAqh5W1T0ehR8NrFPVIL3+aqQB0Mh3N3wm0ZVCOQaYp6oFqloMzAXOrWmQEP+Gz8YlVXyP50QTT1VXqWpEd/KHiPeu7zsDzMOVlYkm3r6Al42pwf+RKn4G/B24rSaxqokXkRDxrgWmqOoh32e2ebU/v2ROBB2AbwNebyLKH7SxIiK5wADcb/DRxkoVkaXANuA9VY0m5j9w/7hLox1XAAXeFZHFvgqxkToK2A485Tt19YSIBOlaHpELgRejCaCq3wH3AhuBLcBeVX03ipArgJNFpIWIZAKnU77mVjTaqOoWcL+gAK09ihsLVwJvRxtERCaLyLfARCCC7vXlYp0FfKeqy6IdV4AbfKevnqzJqboQegAnich8EZkrIsd7McBAyZwIJMi6pLvWVUSaAP8Gbq7wm0pEVLVEVfNwvzUNFpE+EY7rTGCbqi6OdkwVDFfV43CNha4XkZMjjNMAdwj8iKoOAH7AndaIiog0BM4CXokyTjPcb9pdgfZAYxG5ONJ4qroKd1rkPeD/cKc6i6vcqI4RkUm47zw92liqOklVO/li3VDd56sYUyYwiSiTSQWPAN2APNwvEfdFGa8B0AwYAtwKvCwiwX4+RiyZE0HSVykVkTRcEpiuqq95Gdt3mmQOMDbCEMOBs0QkH3dabZSIPO/BuDb7HrcBM3Cn8CKxCdgUcMTzKi4xROvHwOequjXKOGOAb1R1u6oWAa8Bw6IJqKrTVPU4VT0Zd/i/Jsox+m0VkXYAvkfPTx1ES0QuA84EJqrvRLdHXgB+EsX23XDJfpnv/0pH4HMRaRtpQFXd6vuFrhR4nMj/j/htAl7znTpegDvCD3tCOxzJnAgWAt1FpKvvt7wLgTcSPKYyvow8DVilqn/zKGYr/xUVItII98NodSSxVPUOVe2oqrm4P7tZqhrxb7S+MTUWkSz/c9wkYERXYKnq98C3ItLTt2o08GU04/O5iChPC/lsBIaISKbv73o0bh4oYiLS2vfYGTgPb8YJ7v/FZb7nlwGvexTXEyIyFvgtcJaqFngQr3vAy7OI8P8IgKp+oaqtVTXX939lE+4ikO+jGF+7gJfnEuH/kQAzgVG+2D1wF1Z4Wd00ea8a8v3ScDruapx1wCQP4r2IO1Qrwv2F/yyKWCfiTlUtB5b6ltOjHF8/YIkv5grgDx79OY7Eg6uGcOf1l/mWldH+neAOnRf5vu9MoFmU8TKBnUCOR39uf8T9kFkBPIfvqo0o4n2ES3bLgNERxqj0bxhoAXyAO8L4AGgeZbxzfc8PAVuBd6KMtxY33+f/f1KTq3yCxfu37+9kOfAfoEM08Sq8n0/NrhoKNr7ngC9843sDaBdlvIbA877v/Dkwyot/34GLlZgwxph6LplPDRljjIkDSwTGGFPPWSIwxph6zhKBMcbUc5YIjDGmnrNEYOodEfmziIwUkXPEg6q2UYwjvyaVLo2JFUsEpj46AVcXagTu2n5j6jVLBKbe8NXFXw4cD3wGXAU8IiKV6syIyNMi8qiIfCQiX/tqN/mLAv5VRBb6iopd41svvvUrxPVrGO9bP1JEPvTV4f/SF7PS/zsRuVhcL4qlIvKYiKTG8s/CmEANEj0AY+JFVW8VkVeAS3B9Guao6vAqNsnFHTV0A2aLyNHApbhKpMeLSDrwiYi8i6uTlAf0x9WBWSgiH/riDAZ642rd/x+uvMSr/p2IyDHAeFxBvyIReRhXVfNZb765MVWzRGDqmwG4Mge9qL620cvqCoetEZH1vm1OBfqJyPm+z+QA3XElR15U1RJcEbi5uCOPfcACVV0PICIv+j77asB+RgMDcckDoBFJWDjO1F2WCEy9ICJ5wNO46pI7cHWJxNf7YaiqHgyyWcX6K4orj36jqr5TIX5V7QiDxSm3OfCMqt5R5ZcwJkZsjsDUC6q6VF2fh69xp2lmAaepal6IJABwgYikiEg3XMG9r3DtDK/1lSBHRHr4KrF+CIz3zSG0wnVfW+CLM9hXRTcFdwro4wr7+QA4P6A6aXMR6eLVdzemOnZEYOoN3w/o3apaKiK9VLW6U0Nf4VpKtgF+oaqFIvIEbu7gc1956u241pAzgKG4yqIK3Kaq34tIL9zE9BSgLy5hzAjciap+KSK/w3V+S8FVnrweN6dgTMxZ9VFjghCRp3Glu1+t7rPVxBmJa2R+phfjMiYW7NSQMcbUc3ZEYIwx9ZwdERhjTD1nicAYY+o5SwTGGFPPWSIwxph6zhKBMcbUc/8fsI/H6mLU4ZIAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "NN = np.arange(1, 17)\n", "suffocation_times_h = [time_of_suffocation(n).m_as('hour') for n in NN]\n", "death_times_h = [time_of_death(n).m_as('hour') for n in NN]\n", "plt.plot(NN, suffocation_times_h, 'o-', color='orange', label='time to vomiting & nausea')\n", "plt.plot(NN, death_times_h, 'o-', color='r', label='time of death')\n", "plt.legend()\n", "plt.yticks(np.arange(0, 25, 2))\n", "plt.xticks(np.arange(0, 17, 1))\n", "plt.grid()\n", "plt.xlabel('# people')\n", "plt.ylabel('time (h)')\n", "#plt.savefig('death2.pdf')" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "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.7.0" }, "toc": { "nav_menu": {}, "number_sections": true, "sideBar": true, "skip_h1_title": false, "toc_cell": false, "toc_position": {}, "toc_section_display": "block", "toc_window_display": false }, "varInspector": { "cols": { "lenName": 16, "lenType": 16, "lenVar": 40 }, "kernels_config": { "python": { "delete_cmd_postfix": "", "delete_cmd_prefix": "del ", "library": "var_list.py", "varRefreshCmd": "print(var_dic_list())" }, "r": { "delete_cmd_postfix": ") ", "delete_cmd_prefix": "rm(", "library": "var_list.r", "varRefreshCmd": "cat(var_dic_list()) " } }, "types_to_exclude": [ "module", "function", "builtin_function_or_method", "instance", "_Feature" ], "window_display": false } }, "nbformat": 4, "nbformat_minor": 2 }