diff --git a/miTrabajo/Tabla3.ipynb b/miTrabajo/Tabla3.ipynb new file mode 100644 index 0000000000000000000000000000000000000000..e491894b571a073b775be5fd629719dc6c9d6b2b --- /dev/null +++ b/miTrabajo/Tabla3.ipynb @@ -0,0 +1,131 @@ +{ + "cells": [ + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "#### Con los datos de la tabla 3 del artÃculo, verifique que la frecuencia de dicéntricos más anillos por célula aumenta con el cuadrado de la dosis de radiación inmediatamente después del accidente y 15 meses después del accidente." + ] + }, + { + "cell_type": "code", + "execution_count": 8, + "metadata": {}, + "outputs": [], + "source": [ + "import numpy \n", + "import matplotlib.pyplot as plt \n", + "import scipy.stats as stats" + ] + }, + { + "cell_type": "code", + "execution_count": 36, + "metadata": {}, + "outputs": [ + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYIAAAEGCAYAAABo25JHAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8rg+JYAAAACXBIWXMAAAsTAAALEwEAmpwYAAAwUElEQVR4nO3dd3hUZfbA8e8hBAg1AqEkEAgIUQQkEIplFzsoFuyKFUUsq+uuu6yyuurq+tNd7CsrgqIiKCoiooLYwUYJhA6B0JMgoQUCSUg7vz/mBoeQcgOZzGTmfJ4nT3LL3HuiYc7c9z3v+4qqYowxJnTV8XcAxhhj/MsSgTHGhDhLBMYYE+IsERhjTIizRGCMMSGurr8DqKqWLVtqx44d/R2GMcbUKosXL96lqlFlHat1iaBjx44kJSX5OwxjjKlVRGRLecesacgYY0KcJQJjjAlxlgiMMSbEWSIwxpgQZ4nAGGNCnCUCY4wJcZYIjDEmxFkiMMaYAKeqvPT1etZs3++T69e6AWXGGBNKVJX/zEnh1e83kFdYxMltm1b7PeyJwBhjAtgLX63j1e83MKx/LKMuiPfJPVwlAhHpICLnOT9HiEgTn0RjjDHmsJe+Xs/L36YyIK4536/NpPPfZ3HGM98yIzm9Wu9TaSIQkTuAacBrzq52wIxqjcIYY8wRXvl2PS98vY5+HZuzdFsWGfvyUCA9K5fR01dUazJw80TwB+AMYD+Aqq4HWlVbBMYYY47w6vcbePbLdVyeEEPa3hzyCouPOJ5bUMSYOSnVdj83ieCQquaXbIhIXcBWvDfGGB+YMG8j//5iLZeeGs2zV5/K9n15ZZ6XkZVbbfd0kwjmisjfgQgROR/4EPi02iIwxhgDwMQfN/HUrDUM6dmW5685lbA6QnRkRJnnlrf/WLhJBA8BO4EVwJ3ALFV9uNoiMMYYw6RfNvPEZ6u5sHsbXry2F3XDPG/PowbFExEedsS5EeFhjBpUfRVEbsYR3KeqLwETSnaIyP3OPmOMMcdp8vwtPPrJKi7o1pqXr08gPOy3z+hDE2IAGDMnhYysXKIjIxg1KP7w/uogqhU394vIElXtXWpfsqomVFsUVZCYmKi2QpkxJlhMXbiVh6av4LyTW/G/G/pQr65vhneJyGJVTSzrWLlPBCJyPTAMiBORmV6HmgB7XNx0InAxkKmq3cs4fgPwICBANnC3qi6r7LrGGBMsPkjaxuiPV3B2fBRjb+jtsyRQmYqahn4GtgMtgee89mcDy11c+y3gFWBSOcc3AQNVda+IXAiMB/q7uK4xxtR605ek8eBHyznzxJa8emMf6tcNq/xFPlJuIlDVLcAW4LRjubCqzhORjhUc/9lrcz6egWrGGBP0Plmazl8/XMbpnVsw4eZEGoT7LwmAu5HFV4jIehHZJyL7RSRbRKp7CrzbgdkVxDBSRJJEJGnnzp3VfGtjjKk5ny7L4M/vL6VfXHNev7mv35MAuCsf/Q9wqao2U9WmqtpEVatt+jsRORtPIniwvHNUdbyqJqpqYlRUVHXd2hhjatTsFdv50/tLSezQnIm39iWinv+TALgrH92hqmt8cXMR6Qm8Dlyoqrt9cQ9jjAkEc1b9yn3vJZPQPpKJw/vSsF7grALgJpIkEXkfz0Rzh0p2qur047mxiMQC04GbVHXd8VzLGGMC2derd3Dvu0vo0a4Zbw7vS+P6gZMEwF0iaArkABd47VM8b+LlEpH3gLOAliKSBjwGhAOo6jjgUaAF8D8RASgsr8bVGGNqq+/WZnLPlCV0a9uUt2/rR5MG4f4O6SiVJgJVHX4sF1bV6ys5PgIYcSzXNsaY2mDuup3cOXkx8W2aMOn2/jQNwCQA7qqGuorINyKy0tnuKSKP+D40Y4ypvX5cv4uRk5I4Maox79zej2YRgZkEwF3V0ARgNFAAoKrLget8GZQxxtRmP2/YxYhJi4hr2YgpI/oT2bCev0OqkJtE0FBVF5baV+iLYIwxprZbsHE3t7+VRGzzhkwZ0Z8TGgV2EgB3iWCXiHTGWYxGRK7CM/WEMcYYL0mb9zD8rUXEnBDBlBEDaNG4vr9DcsVN1dAf8MwDdJKIpOOZI+hGn0ZljDG1TNLmPdwycSFtmjXg3Tv6E9WkdiQBcFc1tBE4T0QaAXVUNdv3YRljTO3xU+ouRrydRNtmDXjvjgG0atLA3yFVSUXTUD9Qzn4AVPV5H8VkjDG1xjdrdnD3lCXEtWjE5BG160mgREVPBE2c7/FAX6BkTYJLgNKdx8YYE3I+X76d+6cm0y26KW8P71crOobLUtE01P8EEJF5QO+SJiEReRz4vEaiM8aYADVtcRp/m7aMPh1O4I1b+wbsYDE33HQWtwbyvbbznX3GGBOS3vllM//4ZBVnntiS8Tf3CagJ5I6Fm+gnAQtF5GNneyie1ceMMSbkjJ+3gf+btZbzTm7NK8MSAmI9gePlpmroKRGZDfzO2TVcVZN9G5YxxgQWVeXFr9fz0jfrubhnW164thfhYf5ZY7i6VVQ11FRV94tIc2Cz81VyrLmqVrqAvTHGBANV5f9mrWHCD5u4uk87nrmyJ2F1xN9hVZuKngjeBS4GFuOMKnaIs93Jh3EZY0xAKC5W/vHJSqYs2Motp3XgsUtOoU4QJQGouGroYud7XM2FY4wxgaOwqJi/TVvO9OR07hrYmQcHxx8eSxVM3ExDfbmINPPajhSRoT6Nyhhj/Cy/sJg/Tk1menI6fzm/a9AmAXA36dxjqrqvZENVs/CsNmaMMUEpr6CIO99JYtaKX3lkyMncd26XoE0C4K58tKxkUbuLZo0xphwHDxUydOxPrM88AMCbP22mZeP6DE2I8XNkvuPmiSBJRJ4Xkc7O1/N4OpCNMSao7MstYMjLPxxOAgDpWbmMnr6CGcnpfozMt9wkgvvwjCZ+H5gK5OGZmtoYY4LGnoP5DJswn827c446lltQxJg5KX6Iqma4GVB2EHioBmIxxhi/yNyfxw2vL2DrnqOTQImMrNwajKhmuaka+kpEIr22TxCROS5eN1FEMksWvS/juIjIyyKSKiLLRaR3lSI3xphqkLY3h6tf+4X0rFzeGt6PmMiIMs+LLmd/MHDTNNTSqRQCQFX3Aq1cvO4tYHAFxy8EujhfI4FXXVzTGGOqzaZdB7lm3C/sPZjP5BH9Oa1zC0YNiiei1PxBEeFhjBoU76cofc9NIigWkdiSDRHpwJEjjcukqvOAiqahuAyYpB7zgUgRaesiHmOMOW7rdmRzzWu/kFdYzHsjB9A79gQAhibE8PQVPYiJjECAmMgInr6iR1BXDbkpA30Y+FFE5uKZXuJ3eD7BH68YYJvXdpqzb3vpE0VkZMk9Y2NjSx82xpgqWZm+j5veWEB4WB0+uHMAJ7ZqcsTxoQkxQf3GX5qbzuIvnPb7Ac6uP6nqLt+GdVQM44HxAImJiZU+jRhjTHkWb9nDrRMX0TQinHfv6E+HFo38HZLfuR0YVgRkAg2AbiJS0vRzPNKB9l7b7Zx9xhjjEz+n7mLEpCRaN23AlBH9g7oDuCrcVA2NAOYBc4B/Ot8fr4Z7zwRudqqHBgD7VPWoZiFjjKkOny3P4NY3F9H+hIa8f+cASwJe3DwR3I9n8fr5qnq2iJwE/F9lLxKR94CzgJYikoZnfqJwAFUdB8wCLgJSgRxg+LH8AsYYU5mJP27iyc9X07dDcybcnEizhrV3fWFfcJMI8lQ1T0QQkfqqulZEKq2jUtXrKzmu2AhlY4wPFRcrz3yxlvHzNjL4lDa8eF2voFhasrq5SQRpzoCyGcBXIrIX2OLLoIwx5njlFxbzt2nLmLE0g5udBWWCaVWx6uSmauhy58fHReQ7oBnwhU+jMsaY45CdV8Ddk5fwY+ouRg2K556zOgf1NNLHq0rTSavqXF8FYowx1SFzfx63vrmIdTuyefbqU7mqTzt/hxTwbF0BY0zQ2LDzALdMXMieg/m8fksiZ8W7mQ3HWCIwxgSFJVv3cvtbiwirI0wdOYCe7SL9HVKtYYnAGFPrfb16B/e+t4TWTRsw6bZ+Nlq4ispNBCKSzW+Ty5X0sqjzs6pqUx/HZowxlZq6cCt//3gF3WOaMfHWvrRsXN/fIdU65SYCVW1S3jFjjPE3VeWlb9bz4tfrOSs+irHDetOovjVyHAtX/9VE5Eygi6q+KSItgSaqusm3oRljTNkKi4r5xycreW/hNq7q046nr+hBeJibWfVNWSpNBCLyGJAIxANvAvWAycAZvg3NGGOOlptfxH3vLeHrNZnce/aJ/OWCrjZG4Di5eSK4HEgAlgCoaoaIWLORMabG7TmYz+1vL2LptiyeHNqdmwZ08HdIQcFNIshXVRURBRAR6443xtS4bXtyuGXiQtKycnn1hj4M7t7G3yEFDTeJ4AMReQ3PUpJ3ALcBE3wbljHG/GZl+j6Gv7WI/MJipozoT9+Ozf0dUlBxM9fQsyJyPrAfTz/Bo6r6lc8jM8YY4Mf1u7hr8mKaNqjLu3edRpfW1jJd3VxVDTlv/Pbmb4ypUTOS0/nrh8s4sVVj3hrejzbNGvg7pKDkdkDZUWxAmTHGV1SV8fM28vTstQzo1JzXbkqkWYQtJuMrlQ4oE5Enge3AO3hGFd8AtK2R6IwxIae4WHny89W8+dNmhvRsy/PXnEr9uraYjC+5aRq6VFVP9dp+VUSWAY/6KCZjTIjKKyjiLx8u4/Pl27ntjDgeGXIydWwxGZ9zkwgOisgNwFQ8TUXXAwd9GpUx5pjNSE5nzJwUMrJyiY6MYNSgeIYmxPg7rEplZudx5zuLSd6axegLT2Lk7zvZQLEa4iYRDANecr4U+MnZZ4wJMDOS0xk9fQW5BUUApGflMnr6CoCATgarMvZxx9tJ7M0p4NUbenNhD2t9rkluykc3A5f5PhRjzPEaMyflcBIokVtQxJg5KQGbCL5YuZ0/v7+MyIbhfHjXaXSPaebvkEKOTdVnTBDJyMqt0n5/UlXGfpfKs1+uo1f7SMbf1IdWTa081B98Ol2fiAwWkRQRSRWRh8o4Hisi34lIsogsF5GLfBmPMcEuOjKiSvv9Ja+giD+/v5Rnv1zHZb2imTpygCUBP/JZIhCRMGAscCHQDbheRLqVOu0R4ANVTQCuA/7nq3iMCQWjBsUTEX5kqWVEeBijBsX7KaKjZWbncd34+cxYmsGoQfG8eG0vGoRbeag/VZoIROR+EWkqHm+IyBIRucDFtfsBqaq6UVXz8VQdle5rUKBkYFozIKMqwRtjjjQ0IYanr+hBTGQEAsRERvD0FT0Cpn9gZfo+LnvlJ1J+zWbcjb35w9knWmVQAHDTR3Cbqr4kIoOAE4Cb8Awu+7KS18UA27y204D+pc55HPhSRO4DGgHnlXUhERkJjASIjY11EbIxoWtoQkzAvPF7s07hwOWmaagkXV8EvKOqq7z2Ha/rgbdUtV3J9UXkqJhUdbyqJqpqYlRUVDXd2hhTE1SVV75dz12TlxDfpgmf3HuGJYEA4+aJYLGIfAnEAaOdRWmKXbwuHWjvtd3O2eftdmAwgKr+IiINgJZApovrG2MCXF5BEX+btpyZyzIY2iuaZ67saf0BAchNIrgd6AVsVNUcEWkBDHfxukVAFxGJw5MAruPogWhbgXOBt0TkZKABsNNl7MaYAJa5P4873lnMsm1ZjBoUzz1ndbb+gADlJhEonqqfi4En8LTlV1rnpaqFInIvMAcIAyaq6ioReQJIUtWZwF+ACSLyZ+c+t6pquTOeGmNqh5Xp+7hjUhJZOQWMu9FWEwt0Utn7roi8iqcp6BxVPVlETgC+VNW+NRFgaYmJiZqUlOSPWxtjXJi9Yjt//mApzRvWY8ItiZwSbf0BgUBEFqtqYlnH3DwR9FfV3iKSDKCqe0WkXrVGaIyp9Tydwqk899U6EmIjee2mPrRqYoPEagM3iaDAGRxWsnh9FO46i40xIcK7U/hyZyyDdQrXHm4SwcvAx0ArEXkKuArPiGBjjPF0Ck9KYnn6Pv42OJ67B1qncG3jZvbRKSKyGE91jwBDVXWNzyMzxgS8len7GPF2EvvzPJ3Cg06xTuHaqNJEICKxQA7wqfc+Vd3qy8CMMYFt1ortPOB0Ck+763S6Rdsy5rWVm6ahz/H0DwiestE4IAU4xYdxGWMClKry8jepvPD1OnrHRvLaTYlENanv77DMcXDTNNTDe1tEegP3+CwiY0zAOnCokAenLefzFdu5IiGG/7NO4aBQ5YVpVHWJiJSePM4YE+TW7cjmrsmL2bzroK0pHGTc9BE84LVZB+iNTRdtTEgpWQu5Uf26TBkxgNM6t/B3SKYauXkiaOL1cyGePoOPfBOOMSaQHCos4l+freGd+Vvo17E5rwxLsJXEgpCbPoJ/AohIU8+mZvs8KmOM36Vn5XLPlCUs25bFyN93YtSgeMLDfLq6rfETN01DicCbOE8GIrIPuF1VbcIfY4LU3HU7+dPUZAqKlHE39mZw97b+Dsn4kJumoYnAPar6A4CInOns6+nLwIwxNa+4WPnvt6m8+M06urZqwqs39qZTVGN/h2V8zE0iKCpJAgCq+qOIFPowJmOMH+w9mM+f3l/K3HU7uSIhhn9d3p2G9apcWGhqITf/l+eKyGvAe3gGll0LfO+MJ0BVl/gwPmNMDVi2LYt7pixhZ/Yhnrq8O8P6xVppaAhxkwhOdb4/Vmp/Ap7EcE61RmSMqTGqyuQFW3ny09VENanPtLtPo2e7SH+HZWqYm6qhs2siEGNMzcrJL+Thj1fycXI6Z8VH8cI1vTihkS01EooqrQUTkftFpKl4vC4iS0TkgpoIzhjjGxt3HuDysT8zY2k6D5zflYm39LUkEMLcFAXfpqr7gQuAFsBNwDM+jcoY4zOzV2zn0ld+IjM7j0m39eOP53ahTh3rDwhlbvoISv5CLgImOQvQ21+NMbVMQVEx/569ltd/3ESv9pGMvaE3MZER/g7LBAA3iWCxiHyJZ/rp0SLSBFuq0phaZcf+PO59dwmLNu/lltM68PCQbtSra6OEjYebRHA70AvYqKo5ItICGO7m4iIyGHgJCANeV9WjmpRE5BrgcTwVSMtUdZi70I0xbvyyYTf3vbeEg4eKeOm6XlzWK8bfIZkA4yYRKNANuBh4AmiEZ4GaCjkL3o8FzgfSgEUiMlNVV3ud0wUYDZyhqntFpFXVfwVjTFlUlXFzNzJmzlriWjbivTsG0KV1k8pfaEKOm0TwPzxNQefgSQTZeGYf7VvJ6/oBqaq6EUBEpgKXAau9zrkDGKuqewFUNbNK0RtjyrQvt4C/friMr1bvYEjPtvz7yp40rm+jhE3Z3Pxl9FfV3iKSDOB8cndTZxYDbPPaTgNKL2jTFUBEfsLTfPS4qn5R+kIiMhIYCRAbG+vi1saEruSte7l/6lIysnJ59OJuDD+jo40SNhVykwgKnGYeBRCRKKqvs7gu0AU4C2gHzBORHqqa5X2Sqo4HxgMkJiZqNd3bmKBSVKz877tUXvxmPW2aNuD9OwfQp0Nzf4dlagE3ieBl4GOglYg8BVwFPOLidelAe6/tds4+b2nAAlUtADaJyDo8iWGRi+sbYxzb9uTw5/eXkrRlL5f1iuaJy7rTLCLc32GZWsLNFBNTRGQxcC6eMQVDVXWNi2svArqISByeBHAdULoiaAZwPfCmiLTE01S00X34xpgZyen8Y8ZKAF68thdDE6wqyFRNuYlARLyfKTPxzD56+Jiq7qnowqpaKCL3AnPwtP9PdAajPQEkqepM59gFIrIaKAJGqeruY/91jAkd+/MK+MeMlXyyNIPEDifwwrW9aN+8ob/DMrWQqJbd5C4im/D0CwgQC+x1fo4EtqpqXA3FeITExERNSrLF0UxoW7hpD39+fym/7s/j/nO7cM9Znalry0iaCojIYlVNLOtYuU8EJW/0IjIB+FhVZznbFwJDfRCnMaYSBUXFvPzNesZ+l0r75g2ZdtdpbNmdw8Ax35ORlUt0ZASjBsVb85CpEjcfIQaUJAEAVZ0NnO67kIwxZdm86yBXjfuF/36bypW92/H5H3/Hlt05jJ6+gvSsXBTPgvOjp69gRnLpugxjyuemaihDRB4BJjvbNwAZvgvJGONNVfkwKY3HP11F3TrC2GG9GdLTs5j8mDkp5BYUHXF+bkERY+ak2FOBcc1NIrgez+pkH+PpM5jn7DPG+FhWTj6jp69g9spfGdCpOc9f04torxlDM7Jyy3xdefuNKYub8tE9wP01EIsxxsvPqbt44INl7D54iIcuPIk7fteJsFLrBkRHRpBexpt+tE0vbarAygyMCTCHCot4etYabnhjAQ3rhzH97jO4a2Dno5IAwKhB8USEhx2xLyI8jFGD4msqXBMEbBYqYwJIamY2909dyqqM/dzQP5aHh5xMw3rl/zMt6QcYMyfFqobMMbNEYEwAUFWmLNjKvz5fTcN6dZlwcyLnd2vt6rVDE2Lsjd8cl0oTgYh0BV4FWqtqdxHpCVyqqv/yeXTGhIBdBw7x0EfL+XpNJr/r0pLnrj6VVk0rXfLDmGrjpo9gAp7FYwoAVHU5nnmDjDHH6fuUTAa/+APz1u/i0Yu78fbwfpYETI1z0zTUUFUXlprPvNBH8RgTEvIKinhm9lre+nkz8a2bMHlEP05q09TfYZkQ5SYR7BKRzvy2HsFVwHafRmVMEFu4aQ8PfbScjbsOcuvpHXnowpNoUKryx5ia5CYR/AHPojAniUg6sAnP6GJjTBUcOFTIv2ev5Z35W2h3QgSTb+/PmV1a+jssY9wtXq+q54lII6COqmY7awwYY1z6LiWTh6evYPv+PG47I46/XNCVRraGsAkQbv4SPwJ6q+pBr33TgD6+CcmY4LH3YD5Pfraa6cnpnNiqMdPuOp0+HU7wd1jGHKGihWlOAk4BmonIFV6HmgJW1mBMBVSVWSt+5bGZK8nKKeCP55zIH845kfp1rS/ABJ6KngjigYvxLERzidf+bOAOH8ZkTK2WuT+PR2as5MvVO+gR04xJt/WnW7RVBJnAVdHCNJ8An4jIaar6Sw3GZEytVDJd9JOfrya/sJjRF57E7WfG2cphJuC56SNIFZG/Ax29z1fV23wVlDG1zbY9ngVifkzdRb+45jxzRQ86RTX2d1jGuOImEXwC/AB8jWeBeWOMo6hYefvnzYyZk0JYHeFfQ7szrF8sdcqYKdSYQOV2ZPGDPo/EmFpm/Y5sHvxoOUu2ZnF2fBRPXd7D1gEwtZKbxsvPROSiY7m4iAwWkRQRSRWRhyo470oRURFJPJb7GFOTCoqK+e836xny8o9s2nWQF6/txcRb+1oSMLWWmyeC+4G/i0g+kA8InkFmFZZBiEgYMBY4H0gDFonITFVdXeq8Js49FhxD/MbUqBVp+xg1bRlrf83mklOjeeySbrRsXN/fYRlzXNwsVdnkGK/dD0hV1Y0AIjIVuAxYXeq8J4F/A6OO8T7G+FxeQREvfL2OCfM2EtWkfpXWCzAm0LlZj0DwzC0Up6pPikh7oK2qLqzkpTHANq/tNKB/qWv3Btqr6uciUm4iEJGRwEiA2NjYykI2plot2Libh6avYNOug1zfrz0PXXgyzSLC/R2WMdXGTdPQ/4Bi4Bw8n94P4Gny6Xs8NxaROsDzwK2Vnauq4/FMfEdiYqIez32NcSs7r4B/f7GWyfO3Etu8Ie+O6M/pJ9okcSb4uEkE/VW1t4gkA6jqXhGp5+J16UB7r+12zr4STYDuwPfOWgdtgJkicqmqJrmK3hgfUFVmLsvg6VlryczOY8SZcTxwQdcK1w42pjZz85dd4HT8lqxHEIXnCaEyi4Auzkyl6XhWNRtWclBV9wGHP16JyPfAXy0JGH9anbGfx2euYuHmPfSIaca4m/rQq32kv8MyxqfcJIKXgY+BViLyFHAV8EhlL1LVQhG5F5gDhAETVXWViDwBJKnqzOOI25hqtfdgPs99lcK7C7YS2bAez1zRg2sS29vAMBMSRLXyJndnJtJz8ZSOfqOqa3wdWHkSExM1KckeGkz1KCpW3lu4lWe/TCE7r5CbBnTgz+d1pVlD6ww2wUVEFqtqmWO13DZ67sAzzURdIEJEeqvqkuoK0Bh/WLR5D499sorV2/czoFNzHr/0FFs32IQkN+WjT+Kp7NmA00/gfD/Hd2EZ4zu/7svjmdlrmLE0g+hmDRg7rDcX9WiDU7RgTMhx80RwDdBZVfN9HYwxvnSosIiJP27mv9+up7BY+eM5J3LXWZ2tGsiEPDf/AlbiWZwm07ehGOM7363N5InPVrNp10HO79aafwzpRmyLhv4Oy5iA4CYRPA0ki8hK4FDJTlW91GdRGVNNNu86yJOfreabtZl0imrE27f1Y2DXKH+HZUxAcZMI3sYzF9AK3I0fMMbvDh4qZOx3qbz+wybCw4S/X3QSt54eR726tlqYMaW5SQQ5qvqyzyMxphp4jwr+dX8eV/SO4aHBJ9GqaQN/h2ZMwHKTCH4QkaeBmRzZNGTloyageI8K7h7TlLE3JNCnQ3N/h2VMwHOTCBKc7wO89ln5qAkYWTn5PPflOqYs2EJkw3o87YwKDrNRwca44mY9grNrIhBjqqqwqJipi7YdHhV882kdbVSwMcfAVQG1iAwBTgEON7Sq6hO+CsqYiqgqs1f+ynNfprBh50EbFWzMcXIzsngc0BA4G3gdz6RzlS1KY0y1U1V+WL+LMXNSWJG+jy6tGjPuxj4MOqV1maOCZySnM2ZOChlZuURHRjBqUDxDE2L8ELkxgc3NE8HpqtpTRJar6j9F5Dlgtq8DM8bbkq17+c8Xa5m/cQ8xkRE8e/WpXJ4Qc7gfoPSb/tknRfHR4nRyC4oASM/KZfT0FQCWDIwpxU0iyHW+54hINLAbaOu7kIz5Tcqv2YyZk8LXa3bQsnE9Hr+kG9f3j6V+3bDD58xITmf09BVHvOlPmb+V0vPq5hYUMWZOiiUCY0pxkwg+E5FIYAywBE/F0Ou+DMqYrbtzeOHrdcxYmk7j+nX56wVdGX5GHI3qH/0nO2ZOyuEkUKK8ydUzsnLLOWJM6HKTCP6jqoeAj0TkMzwdxnm+DcuEqsz9efz321TeW7iVsDrCyN934u6BnYlsWP7qqFV5c4+OjKiOMI0JKm4SwS9AbwAnIRwSkSUl+4ypDvtyChg3bwNv/rSJwiLl2r7t+eO5XWjtYkRwdGQE6WUkA+HIJ4OI8DBGDYqvvqCNCRLlJgIRaQPE4FmIJgHPvyuApniqiIw5bjn5hbz502Zem7uB7EOFXHpqNA+c35UOLRq5vsaoQfFH9BGA503/yj4xfLd2p1UNGVOJip4IBuFZkKYd8LzX/mzg7z6MyYSA/MJipi7aysvfpLLrwCHOO7kVf7kgnpPbHj0WoLIy0JKfrVTUmGNT6ZrFInKlqn5UQ/FUytYsrt2KipVPlqbzwtfr2LYnl35xzXlwcHy5cwKVrggCz6f9p6/oYW/0xlTBca1ZrKof2chic7xUla9W7+DZL1NYt+MAp0Q35a3h3RnYNarCJSLLqgiyMlBjqpeNLDY+93PqLv4zJ4Wl27Lo1LIRY4f15sLubajjYlK48iqCrAzUmOrj05HFIjIYeAkIA15X1WdKHX8AGAEUAjuB21R1S5V+AxOQVJV563cx7vsN/LJxN22bNeDfV/bgyt7tqBvmfnGY8iqCrAzUmOrjs5HFIhIGjAXOB9KARSIyU1VXe52WDCSqao6I3A38B7i2Kr+ACSyFRcXMWvkr477fwOrt+2ndtD6PDDmZGwd0oEF4WOUXKKW8iiArAzWm+vhyZHE/IFVVNwKIyFTgMuBwIlDV77zOnw/c6C5sE2jyCor4cHEaE+ZtZOueHDpFNeI/V/bksoToI6aDqCqrCDLG99x0Fj/p/Hh4ZLGq7nNx7Rhgm9d2GtC/gvNvp5wmJxEZCYwEiI2NdXFrU1P25RTwzvzNvPXzZnYdyKdX+0geHnIy55/c2lUfgBtDE2Lsjd8YH6poQNk5qvqtiFxRxjFUdXp1BSEiNwKJwMCyjqvqeGA8eMpHq+u+5tj9ui+PN37cyLsLtnIwv4iBXaO4+6zO9I9rXmEVkDEm8FT0RDAQ+Ba4pIxjClSWCNKB9l7b7Zx9RxCR84CHgYHOFBYmgKVmHmD8vA18nJxOUbFyyanR3Pn7znSLtkVhjKmtyk0EqvqY8334MV57EdBFROLwJIDrgGHeJzhTV7wGDFbVzGO8j6kByVv3Mm7uBr5cvYN6YXXoH9eC9Tuymbk0g6TNe63d3pharKKmoQcqeqGqPl/J8UIRuReYg6d8dKKqrhKRJ4AkVZ2JpwO6MfCh05ywVVUvreLvYHxEVZm7bievfr+BBZv20CwinHvPPpFWTerzf7PW2qIvxgSJipqGmjjf44G+wExn+xJcDihT1VnArFL7HvX6+TzXkZoaU1hUzOcrtjNu7kbWbN9Pm6YNeGTIyVzXL5bG9etyxjPf2mhfY4JIRU1D/wQQkXlAb1XNdrYfBz6vkehMjcrNL+LDxduY8MNGtu3JpXNUI8Zc1ZPLesVQr+5vg8BstK8xwcXNOILWQL7Xdr6zzwSJfTkFTPrFUwK6+2A+CbGR/GNIN84rpwTURvsaE1zcJIJJwEIR+djZHgq87bOITI1ZkbaPdxdu4ZOlGeTkF3F2fBR3DexMv0pKQG20rzHBxc2AsqdEZDbwO2fXcFVN9m1Yxldy8gv5dFkGUxZsZXnaPhqE1+GSntHcdmZcmWsBlMVG+xoTXNw8EYBn9tH9qvqmiESJSJyqbvJlYKZ6pfyazbsLtjB9STrZhwrp2rox/7z0FIYmxNAsIrzK17PRvsYEDzfTUD+GZ9RvPPAmEA5MBs7wbWjmeOUVFDF75XamzN9K0pa91Aurw0U92nDDgA4kdjjBRgAbYwB3TwSXAwl4JpxDVTNEpEnFLzH+tHHnAd5dsJVpS9LIyikgrmUjHr7oZK7s047mjer5OzxjTIBxkwjyVVVFRAFExP2q4qbG5BcW89XqHUxZsIWfN+ymbh1h0CltGNY/ltM6tai2CeCMMcHHTSL4QEReAyJF5A7gNmCCb8Mybm3bk8N7C7fyQdI2dh3IJ8bpuL06sR2tmjSo/ALGmJDnpmroWRE5H9iPp5/gUVX9yueRmXIVFhXz7dpMpizYyrz1OxHgnJNac8OAWH7fJYow+/RvjKkCV1VDqvqViCwoOV9EmqvqHp9GZo7y1k+beO6rdWTnFQLQLCKc+87pwnV929tgLmPMMXNTNXQn8E8gDygGBM801J18G5oBOFRYxNyUnbz87XpWpu8/8lhBEZ1aNrIkYIw5Lm6eCP4KdFfVXb4OJpTNSE4/PECrbbMGXHJqNLsP5jNn1a9k5xVSVmtPXmGxTfRmjDlubhLBBiDH14GEEu83/ejICM4+KYppSWnkFRYDkLEvj9fmbaRBeB0u7hnNJadGc8vEsid8tYnejDHHy00iGA387PQRHF5BTFX/6LOogtiM5HRGfbiMgmLPipvpWblMnr+1zHObN6zHs1efCkCMTfRmjPGROpWfwmt4lqycDyz2+jLH4LFPVh5OApXZvi/v8M+jBsUTER52xHGb6M0YUx3cPBGEq2qFq5WZym3ceYBPl21nn1Px44b3p32b6M0Y4ytuEsFsERkJfMqRTUNWPlqJtL05fLZ8O58uy2BVxn6qMrVPWZ/2baI3Y4wvuEkE1zvfR3vtq/Xlo6U7bKvr03Xm/jw+X+F581+yNQuAXu0j+cfF3RjSoy0XvjSPvTkFR72uUb0wIhvWs0/7xpga52ZkcVxNBFKTZiSnH7GwyvEuvp62N4fvU3Yya8V25m/cTbHCSW2aMGpQPJf0jCa2RcPD5z52ySmMmraMgqLf+gnCw4SnLu9hb/zGGL9wNbJYRE4HOnqfr6qTfBRTtSv96f/gocLjWnw9r6CIRZv38H3KTuau20lq5gEA4lo24t5zunBJz7Z0aV32BK3W1m+MCTRuRha/A3QGlgIl756KZwnLyl47GHgJCANeV9VnSh2v71ynD7AbuFZVN7sPv3JlffovT0U1+Vt2Hzz8xv/Lht3kFhRRL6wO/Ts157q+7TkrPorOUY1dzfFvbf3GmEDi5okgEeimqu5qHh0iEgaMBc4H0oBFIjJTVVd7nXY7sFdVTxSR64B/A9dW5T6VGTMn5ahP/+XxrtLJzS9i/sbdzF23k+9TMtm82zOmrkOLhlyT2I6B8VEM6NSChvXcLvJmjDGByc272EqgDbC9itfuB6Sq6kYAEZkKXAZ4J4LLgMedn6cBr4iIVDXpVMTtyNsGdetw82kdeOPHTXyfksmCTXvILyymQXgdTuvUgltP78hZ8a3o2NKWYzDGBBc3iaAlsFpEFnJk+eillbwuBtjmtZ0G9C/vHFUtFJF9QAvgiHmNnPLVkQCxsbEuQv5NdDkjck9oGE6D8DC278ujUb0w6oeH8fTstQB0jmrEjf07cFZ8FP3imtOg1EAuY4wJJm4SweO+DqIyqjoeGA+QmJhYpaeFUYPij+gjAKhbR2jZuD6bdx/0XB/oHXsCZ8VHMbBrFO2bNyznasYYE3zclI/OPcZrpwPtvbbbOfvKOidNROoCzfB0Glebkk7ZJz9bze6D+QAUFit1RLjtjDgGdo0isWNz6tV1M9uGMcYEn3ITgYhk4/mwfNQhQFW1aSXXXgR0EZE4PG/41wHDSp0zE7gF+AW4Cvi2OvsHSgxNiKFr6yb899v1DOwaxcD4KNo2s8najDEGKkgEqlp2IbxLTpv/vcAcPOWjE1V1lYg8ASSp6kzgDeAdEUkF9uBJFj7RLbopr97Yx1eXN8aYWsuntY+qOguYVWrfo14/5wFX+zIGY4wxFbOGcWOMCXGWCIwxJsRZIjDGmBBnicAYY0KcJQJjjAlxlgiMMSbEWSIwxpgQJz4YyOtTIrIT2FLGoZaUmqyuFrCYa4bF7Hu1LV4IvZg7qGpUWQdqXSIoj4gkqWqiv+OoCou5ZljMvlfb4gWL2Zs1DRljTIizRGCMMSEumBLBeH8HcAws5pphMftebYsXLObDgqaPwBhjzLEJpicCY4wxx8ASgTHGhLigSAQiMlhEUkQkVUQe8nc8lRGRiSKSKSIr/R2LWyLSXkS+E5HVIrJKRO73d0wVEZEGIrJQRJY58f7T3zG5JSJhIpIsIp/5OxY3RGSziKwQkaUikuTveNwQkUgRmSYia0VkjYic5u+YyiMi8c5/25Kv/SLyp2q9R23vIxCRMGAdcD6QhmeJzOtVdbVfA6uAiPweOABMUtXu/o7HDRFpC7RV1SUi0gRYDAwN1P/OIiJAI1U9ICLhwI/A/ao638+hVUpEHgASgaaqerG/46mMiGwGElW11gzOEpG3gR9U9XURqQc0VNUsP4dVKef9Lh3or6plDaw9JsHwRNAPSFXVjaqaD0wFLvNzTBVS1Xl4luasNVR1u6oucX7OBtYAMf6NqnzqccDZDHe+Av5Tj4i0A4YAr/s7lmAlIs2A3+NZKhdVza8NScBxLrChOpMABEciiAG2eW2nEcBvUMFARDoCCcACP4dSIaeJZSmQCXylqgEdr+NF4G9AsZ/jqAoFvhSRxSIy0t/BuBAH7ATedJrgXheRRv4OyqXrgPeq+6LBkAhMDRKRxsBHwJ9Udb+/46mIqhapai+gHdBPRAK6GU5ELgYyVXWxv2OpojNVtTdwIfAHp+kzkNUFegOvqmoCcBCoDX2L9YBLgQ+r+9rBkAjSgfZe2+2cfaaaOW3tHwFTVHW6v+Nxy3ns/w4Y7OdQKnMGcKnT5j4VOEdEJvs3pMqparrzPRP4GE9zbSBLA9K8nhCn4UkMge5CYImq7qjuCwdDIlgEdBGROCdjXgfM9HNMQcfpfH0DWKOqz/s7nsqISJSIRDo/R+ApJljr16AqoaqjVbWdqnbE83f8rare6OewKiQijZziAZzmlQuAgK6GU9VfgW0iEu/sOhcIyKKHUq7HB81C4HlEqtVUtVBE7gXmAGHARFVd5eewKiQi7wFnAS1FJA14TFXf8G9UlToDuAlY4bS7A/xdVWf5L6QKtQXedqos6gAfqGqtKMesZVoDH3s+J1AXeFdVv/BvSK7cB0xxPjxuBIb7OZ4KOUn2fOBOn1y/tpePGmOMOT7B0DRkjDHmOFgiMMaYEGeJwBhjQpwlAmOMCXGWCIwxJsRZIjDGJRF5XET+WsXXXFobZsQ1oa3WjyMwJpCp6kxsgKMJcPZEYEwFRORhEVknIj8C8c6+XiIyX0SWi8jHInKCs/+PznoNy0VkqrPvVhF5xfn5ahFZ6ayRMM9vv5QxpdgTgTHlEJE+eKZ66IXn38oSPOswTALuU9W5IvIE8BjwJzwTl8Wp6qGS6S1KeRQYpKrp5Rw3xi/sicCY8v0O+FhVc5yZVmcCjYBIVZ3rnPM2nrntAZbjmbbgRqCwjOv9BLwlInfgmQ7FmIBgicCY6jMEGItnJstFInLEE7eq3gU8gme23MUi0qLmQzTmaJYIjCnfPGCoiEQ4M2xegmfu+r0i8jvnnJuAuSJSB2ivqt8BDwLNgMbeFxORzqq6QFUfxbMwivf06cb4jfURGFMOZ33m94FleFY5W+QcugUYJyIN+W3myjBgsrMMogAvq2qWMytniTEi0sU5/o1zXWP8zmYfNcaYEGdNQ8YYE+IsERhjTIizRGCMMSHOEoExxoQ4SwTGGBPiLBEYY0yIs0RgjDEh7v8BOBO/cOfV83cAAAAASUVORK5CYII=\n", + "text/plain": [ + "<Figure size 432x288 with 1 Axes>" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + }, + { + "name": "stdout", + "output_type": "stream", + "text": [ + " 2\n", + "0.01723 x + 0.06768 x - 0.03289\n" + ] + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYgAAAEGCAYAAAB/+QKOAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8rg+JYAAAACXBIWXMAAAsTAAALEwEAmpwYAAAsv0lEQVR4nO3dd3zV9dn/8ddFWGFIAJG9wVQUGUZQcVStonVA7RCquFD0p1att/TWu1Zbu2y5u0UtilbFbZUbW5RaUVvbMsKQAIIEZCRMCTuQef3+ON9giIfkG8jJGXk/H4/zOOc7z5W05uIzvtfH3B0REZGqGsU7ABERSUxKECIiEpUShIiIRKUEISIiUSlBiIhIVI3jHUBdOvbYY71Xr17xDkNEJKksWLDgM3fvUHV/SiWIXr16kZ2dHe8wRESSipmti7ZfXUwiIhKVEoSIiESlBCEiIlEpQYiISFRKECIiElVKzWISEWlopi/KZ9KslWzcuZ8uGelMHJnJ6CFd6+TeShAiIklq+qJ87ns9h/0lZQDk79zPfa/nANRJklAXk4hIkpo0a+XB5FBhf0kZk2atrJP7K0GIiCSpjTv312p/bSlBiIgkqS4Z6bXaX1tKECIiSWriyEzSm6Qdsi+9SRoTR2bWyf01SC0ikqQqBqI1i0lERL5g9JCudZYQqoppF5OZXWRmK80s18zujXL8bjNbbmZLzOxdM+tZ6ViZmS0OXjNiGaeIiHxRzFoQZpYGTAYuAPKA+WY2w92XVzptEZDl7oVm9v+AXwJXBsf2u/vgWMUnIiLVi2ULYhiQ6+5r3L0YeAkYVfkEd3/P3QuDzTlAtxjGIyIitRDLBNEV2FBpOy/YdzjjgbcqbTc3s2wzm2Nmow93kZlNCM7L3rZt21EFLCIin0uIQWozuxrIAs6ptLunu+ebWR9gtpnluPvqqte6+xRgCkBWVpbXS8AiIgnE3TGzOr9vLFsQ+UD3Stvdgn2HMLOvAN8HLnf3oor97p4fvK8B3geGxDBWEZGktKGgkEv/8CHLN+6u83vHMkHMB/qbWW8zawqMAQ6ZjWRmQ4A/EkkOWyvtb2tmzYLPxwIjgMqD2yIiDd62PUVcPXUueTv20zit7lsQMeticvdSM7sdmAWkAU+5+zIzewjIdvcZwCSgFfBq0Dxa7+6XAycAfzSzciJJ7OEqs59ERBq0XftLuOapeWzdXcTzNw3n+I6t6/w7YjoG4e4zgZlV9j1Q6fNXDnPdv4GBsYxNRCRZHSgp46Znssnduocnrz2VoT3axuR7EmKQWkREwiktK+f2FxYyf10Bvx8zhHOO7xCz71KxPhGRJFFe7vz3n3P4+8dbeWjUSVw2qEtMv08JQkQkCbg7P535MX9emMfdFxzPuNN61nzRUVKCEBFJAo++v5qpH37KdWf04jvn9auX71SCEBFJcC/MXc+kWSsZPbgLD1w6ICYPxUVTqwRhZi1iFYiIiHzRzJxN3D89h/O+dByTvjmIRo3qJzlAyARhZmeY2XJgRbA9yMwejWlkIiIN3IerPuOulxYztEdbJn97KE3S6rfTJ+y3/QYYCWwHcPePgLNjFZSISEO3eMNOJjyXTZ8OLZl67amkN02r+aI6FjodufuGKrvK6jgWEREBcrfu4fqn59G+VVOevWEYbVo0iUscYR+U22BmZwBuZk2AO4GPYxeWiEjDlL9zP+OmziOtUSOmjR/Occc0j1ssYVsQtwC3EVnPIR8YDNwao5hERBqk7XuLGDd1LnuLSnn2hmH0bN8yrvGEbUFkuvtVlXeY2QjgX3UfkohIw7O3qJTr/zSf/B37eW78cAZ0OSbeIYVuQfwh5D4REamlotIyJjybzbKNu3n0qqEM690u3iEBNbQgzOx04Aygg5ndXenQMURKeIuIyFEoK3fufHEx/169nV9/axDnn9Ax3iEdVFMXU1Mi6zU0BioXG98NfCNWQYmINATuzv3Tc3h72WZ+cOkArhjaLd4hHaLaBOHuHwAfmNmf3H1dPcUkItIgTJq1khfnbeD2c/sx/sze8Q7nC8IOUjczsylAr8rXuPt5sQhKRCTVPfGPNTz6/mq+PbwH/3Xh8fEOJ6qwCeJV4HHgSfSAnIjIUXltQR4/nfkxlwzszI9HnVRvxfdqK2yCKHX3x2IaiYhIA/DO8i3895+XcGa/Y/n1lYNIq8fie7UVdprrm2Z2q5l1NrN2Fa+YRiYikmLmrNnObS8s5KSubfjjuFNo1jixJ4OGbUFcG7xPrLTPgT51G46ISGpamr+Lm57Jpke7Fjx93am0bBb2z2/8hIrQ3RNveF1EJEl8+tk+rnt6HsekN+G58cNo17JpvEMKJex6EC3M7P5gJhNm1t/MLo1taCIiyW/d9n1c9cQcyh2eHT+Mzm3S4x1SaGHHIJ4Giok8VQ2Rgn0/iUlEIiIpYt32fYydMofCkjKmjR9O3w6t4h1SrYRNEH3d/ZdACYC7FwKJO/QuIhJn67cXHkwOL9x4WkIU36utsAmi2MzSiQxMY2Z9gaKYRSUiksTWby9kzJT/JHVygPCzmB4E3ga6m9nzwAjgulgFJSKSrConh+dvTIyy3Ucq7Cymd8xsIXAaka6lO939s5hGJiKSZNZvL2TsE3MOJocTu7SJd0hHpaZy30Or7NoUvPcwsx7uvjA2YYmIJJcNBZHksK+4lGnjkz85QM0tiF8F782BLOAjIi2Ik4Fs4PTqLjazi4DfEVk74kl3f7jK8buBG4FSYBtwQ0XVWDO7Frg/OPUn7v5MyJ9JRKRebSgoZMyUOewtKuX5G4dzUtfkTw5QwyC1u5/r7ucSaTkMdfcsdz8FGEJkquthmVkaMBm4GBgAjDWzAVVOWwRkufvJwGvAL4Nr2xEZ9xgODAMeNLO2tf3hRERiLVWTA4SfxZTp7jkVG+6+FDihhmuGAbnuvsbdi4GXgFGVT3D394IpswBzgIrVMkYC77h7gbvvAN4BLgoZq4hIvUjl5ADhZzEtMbMngWnB9lXAkhqu6QpsqLSdR6RFcDjjgbequbZrtIvMbAIwAaBHjx41hCQiUjdSPTlA+BbE9cAy4M7gtTzYVyfM7GoiYxyTanutu08Jur6yOnToUFchiYgcVkNIDhB+musB4DfBK6x8oHul7W5EGbcws68A3wfOcfeiStd+ucq179fiu0VEYqKhJAeooQVhZq8E7zlmtqTqq4Z7zwf6m1lvM2sKjAFmVLn/EOCPwOXuvrXSoVnAhWbWNhicvjDYJyISNw0pOUDNLYg7g/daV25191Izu53IH/Y04Cl3X2ZmDwHZ7j6DSJdSK+DVYMm99e5+ubsXmNmPiSQZgIfcvaC2MYiI1JWK5xwaSnIAMHev+SSz3sCmoKuJoC5TR3dfG9vwaicrK8uzs7PjHYaIpJiK5LB7fwkv3HRayiUHM1vg7llV94cdpH4VKK+0XRbsExFJaXk7Pk8Oz9+YesmhOmETROPgWQYAgs/JsSSSiMgRytsRGXOoSA4DuzWc5ADhE8Q2M7u8YsPMRgEq1iciKatycph24/AGlxwg/INytwDPm9kjwXYecE1sQhIRia/K3UrTbhzOyd0y4h1SXIR9DmI1cJqZtQq298Y0KhGROMnfuZ+xT8xhZ2EJzzfg5AAhu5jM7GdmluHue919b/B8gtakFpGUkr9zP2Om/EfJIRB2DOJid99ZsREU0PtqTCISEYmDsMlh+qJ8Rjw8m973/pURD89m+qJqC1sntbBjEGlm1qyiFEbwHESz2IUlIlJ/1m3fx9VT57KzsIRp46tPDve9nsP+kjIgklTuez1S6Hr0kKj1RJNa2BbE88C7ZjbezMYTKb+tBXxEJOkt37ibrz/2H/YciKwEN6h7xmHPnTRr5cHkUGF/SRmTZq2McZTxEXaQ+hdB7aXzg10/dnfVRhKRpDbv0wLGPzOfVs0a89KE0+l3XOtqz9+4c3+t9ie7sF1MuPtbfL5eg4hIUvv78i3c9sJCurZN57nxw+makV7jNV0y0smPkgy6hLg2GYWdxXSamc03s71mVmxmZWa2O9bBiYjEwmsL8rh52gK+1Kk1r91yRqjkADBxZCbpTdIO2ZfeJI2JIzNjEWbchW1BPEKkXPerRBb2uQY4PlZBiYjEyhP/WMNPZ37Mmf2O5fFxp9CqWeiOlIMD0ZNmrWTjzv10yUhn4sjMlByghtp1MeWaWZq7lwFPm9ki4L7YhSYiUnfcnV+8vZLHP1jNJQM78+srB9GscVrNF1YxekjXlE0IVYVNEIXBoj+LzeyXwCbCz4ASEYmr0rJyvv/GUl7O3sBVw3vw0KiTSGtk8Q4r4YX9Iz8uOPd2YB+RpUS/HqugRETqyoGSMm59fiEvZ2/gjvP68ZPRSg5hhZ3mui74eAD4UezCERGpO3sOlHDTs9nMWVPAg5cN4PoRveMdUlIJPzojIpJEtu0p4rqn57Fy8x5+N2YwowY3jHGDuqQEISIpZ0NBIeOmzmXz7gM8cW0W52YeF++QkpIShIiklBWbd3PN1HkUlZbz/I2ncUrPtvEOKWlVmyDM7E3AD3fc3S8/3DERkfqWvbaAG/40n/Smabx6y+kc37H60hlSvZpaEP9bL1GIiByl2Su2cOvzC+ncJp1nbxhG93Yt4h1S0qs2Qbj7BxWfgxLfPdw9NcsWikjSemNRHve8uoQTOrfmT9cP49hWWo2gLoStxXQZsBh4O9gebGYzYhiXiEgoT334Kd99+SOG927HizedpuRQh8I+KPdDYBiwE8DdFwOaUCwicePu/O+slTz0l+VcdGInnrruVFo3bxLvsFJK2FlMJe6+y+yQpw8PO3gtIhJLZeXOD/5vKS/MXc/YYd35yeiBejo6BsImiGVm9m0iS4/2B+4A/h27sEREoisqLeO7Ly9mZs5mbv1yXyaOzKTKP16ljoTtYvoOcCJQBLwI7AbuilFMIiJR7S0q5fqn5zMzZzP3X3IC37voS0oOMRS2FlMh8P3gJSJS7z7bW8QNf5rPso27+dU3B/H1U7rFO6SUF9MH5czsIuB3QBrwpLs/XOX42cBvgZOBMe7+WqVjZUBOsLleD+WJNFwrN+9h/DPz2baniCnjTuH8EzrGO6QGIeyDclcAnYBpwfZYYEt1F5pZGjAZuADIA+ab2Qx3X17ptPXAdcA9UW6x390H1xCfiKS4dz/ewh0vLqJls8a8cvPpDOqeEe+QGoxQD8qZ2a/cPavSoTfNLLuGew8Dct19TXCPl4BRwMEE4e5rg2PltQ9dRFKZu/PEP9fw87dWcGKXY3jymlPp1KZ5vMNqUMIOUrc0sz4VG2bWG2hZwzVdgQ2VtvOCfWE1N7NsM5tjZqMPd5KZTQjOy962bVstbi8iiaqotIzvvbaEn81cwcUndeLVm89QcoiDsNNcvwu8b2ZrAAN6AhNiFlVET3fPDxLTbDPLcffVVU9y9ynAFICsrCw9myGS5LbvLeKWaQuYv3YHd5zfn7vO708jPeMQF2FnMb0dPP/wpWDXCncvquGyfCJLk1boFuwLxd3zg/c1ZvY+MAT4QoIQkdRReTD692OHcPmgLvEOqUEL28WEuxe5+0fBq6bkADAf6G9mvc2sKTAGCFW/yczamlmz4POxwAgqjV2ISOqZvWILVzz6L4pLy3n55tOVHBJAzBYMcvdSM7sdmEVkmutT7r7MzB4Cst19hpmdCrwBtAUuM7MfufuJwAnAH4PB60bAw1VmP4lICNMX5TNp1ko27txPl4x0Jo7MZPSQxFp609158p+f8rO3PubELsfwxDVZdG6THu+wBDD31Om2z8rK8uzsmiZXiTQM0xflc9/rOewvKTu4L71JGj+/YmDCJIni0nK+/0YOry7I46sDO/Grbw4mvWlavMNqcMxsQZWZqkD4ct8jzKxl8PlqM/u1mfWs6yBFpO5MmrXykOQAsL+kjEmzEmNJl+17i7j6ybm8uiCPO87vzyNjhyo5JJiwYxCPAYVmNgj4LyKDxc/GLCoROWobd+6v1f76tHLzHkZN/hcf5e3k92OHcPcFx2umUgIKmyBKPdIXNQp4xN0nA1rsVSSBdcmI3o9/uP31ZfaKLXz9sX9rMDoJhE0Qe8zsPmAc8FczawRoZQ6RBDZxZCbpTQ7tsklvksbEkZlxiScyGL2G8c9k0+vYFvzf7SMYrLIZCS3sLKYrgW8DN7j7ZjPrAUyKXVgicrQqBqITYRZTcWk590/P4ZXsyGD0/35zEC2axmwSpdSRsA/KbTazPwP9g12fEZmeKiIJbPSQrnGfsVSwr5hbpi1g3qcF3HFeP+76isYbksVhE4SZdXP3vODzTURKa7QD+hKpqfQ4cH59BCkiyemTLZEno7fsLuJ3YwYzanBiTK+VcKobgzjLzO4IPt9G5Gnm3QDuvgo4LsaxiUgSe2/FVq549N8cKCnnlZtPV3JIQodNEO7+IrA32Cx29+KKY2bWmGoWEhKRhuvzwej59GzfghkajE5aNa0H8VTw8X0z+x8g3cwuAG4F3ox1cCKSXIpLy/nB9KW8nL2Bi0/qxK++pcHoZBb2f7l7gfFElgC9GZgJPBmroEQk+WzedYDvvLiQ+Wt38J3z+vFdDUYnvbCzmMqBJ4KXiMghPvhkG999eTEHSspUpjuFhEoQZvYpUcYc3L1PlNNFpIEoLSvnN3//hMnvrSazY2smXzWUfse1indYUkfCdjFVrvLXHPgmkSmvItJAbd51gDteXMS8tQWMObU7D152oortpZiwXUzbq+z6rZktAB6o+5BEJNFV7lL6zZWD+NqQbvEOSWIgbBfT0EqbjYi0KDQ1QaSBUZdSwxL2j/yvKn0uBdYC36rzaEQkYalLqeEJ28V0bqwDEZHE9f7Krdz9ykccKCnjt1cOjnt9J6kfYbuY2gMPAmcSmc30IfBQlLEJEUkhpWXl/PqdT3j0/dV8qVNrHvm2upQakrBdTC8B/wC+HmxfBbwMfCUWQYlI/FV+8G3ssEiXUvMm6lJqSMImiM7u/uNK2z8xsytjEZCIxJ+6lATCJ4i/mdkY4JVg+xvArNiEJCLxUrVLafJVQ+nbQV1KDVXYBHETcBfwHGBEprruM7ObAXf3Y2ITnojUl0279nPHi4vUpSQHhZ3F1DrWgYhI/KhLSaIJO4tpBLDY3feZ2dXAUOC37r4+ptGJSEyVlpXzq3c+4TF1KUkU1a0oV9ljQKGZDQL+C1hNpLtJRJLUpl37GfvEHB57fzVjh3Vn+m0jlBzkEGHHIErd3c1sFPCIu081s/GxDExEYue9lVu5++XFFJWWq0tJDitsgthjZvcB44isVd0IaBK7sEQkFkqCWUrqUpIwwiaIK4FvAze4+2Yz6wFMil1YIlLXlm/czcTXPmLZxt2MHdaDBy8boFlKUq1QYxDuvhn4M9As2PUZ8EZN15nZRWa20sxyzezeKMfPNrOFZlZqZt+ocuxaM1sVvK4NE6eIfFFJWTm/+/sqLn/kQ7bsPsDjVw/l51cMVHKQGoWdxXQTMIHIIkF9ga7A48D51VyTBkwGLgDygPlmNsPdl1c6bT1wHXBPlWvbEan9lEWk9tOC4Nod4X4sEYFIq+GeVz9i+abdXD6oCz+6/ETatmwa77AkSYTtYroNGAbMBXD3VWZ2XA3XDANy3X0NgJm9BIwCDiYId18bHCuvcu1I4B13LwiOvwNcBLwYMl6RBq24tJxH38/lkdm5ZLRoyuNXn8JFJ3WKd1iSZMImiCJ3LzYzAMysMVHWqK6iK7Ch0nYeMDzk90W7Nuo0CzObQKR1Q48ePULeXiR1Ldu4i3teXcLHm3YzanAXfniZWg1yZMImiA/M7H+AdDO7ALgVeDN2YYXn7lOAKQBZWVk1JS2RlFVcWs7k93KZ/F6k1fDHcacw8kS1GuTIhU0Q9wLjgRzgZmAm8GQN1+QD3Sttdwv2hZEPfLnKte+HvFakwVGrQWIhbC2mcuCJ4BXWfKC/mfUm8gd/DJGpsmHMAn5mZm2D7QuB+2rx3SINgloNEkvVJggzy6GasQZ3P7maY6VmdjuRP/ZpwFPuvszMHgKy3X2GmZ1KZLpsW+AyM/uRu5/o7gVm9mMiSQYiq9cV1O5HE0ltS/N3cc+rH7Fi8x6+NqQrD142gIwWajVI3TH3w3fbm1nP4ONtwXtF/aWriZT5/sKzDfGUlZXl2dnZ8Q5DJKaKS8t55L1cHn0vl7Ytm/Kzrw3kggEd4x2WJDEzW+DuWVX3V9uCcPd1wcUXuPuQSof+28wWEhmbEJF6olaD1Kewg9RmZiPc/V/BxhmErwQrIkepuLScR2avYvL7q2nXsilPXJOlVoPEXNgEMR54yszaBNs7gRtiEpGIHKJyq+GKIV15QK0GqSdhZzEtAAZVJAh33xXTqESEotIyHpmdy6NBq+HJa7L4iloNUo/CtiAAJQaR+pKTF2k1rNyyhyuGduWBS9VqkPpXqwQhIrG1t6iUR2bn8sQ/19C+ZVOmXpvF+Seo1SDxoQQhKW/6onwmzVrJxp376ZKRzsSRmQm3glp5uTN9cT4Pv7WCrXuK+MYp3fjBJQNo00Lrckn8hC33/U3gbXffY2b3A0OBn7j7wphGJ3KUpi/K577Xc9hfUgZA/s793Pd6DkDCJImPNuzkh28uY9H6nQzq1obHx53C0B5ta75QJMbCTlX9QZAczgS+AkwFHotdWCJ1Y9KslQeTQ4X9JWVMmrUyThF9buueA9zz6keMmvwvNhTs55ffOJk3bh2h5CAJI2wXU8V/YZcAU9z9r2b2kxjFJFJnNu7cX6v99aG4tJyn//Upf5idS1FpGTef3Yfbz+tH6+bqTpLEEjZB5JvZH4msDvcLM2uGHpSTJNAlI538KMmgS0Z6HKKB2Su28OO/fMynn+3jvC8dx/2XnECfDq3iEotITcL+kf8WkaJ7I919J5GlRyfGKiiRujJxZCbpVdZeTm+SxsSRmfUax+pte7nu6Xnc8KdsDHj6+lN56rpTlRwkoYV9UK7QzLYCZwKrgNLgXSShVQxEx2sW0+4DJfzh3VU8/a+1pDdJ4/tfPYFrz+hF08ZqgEviCzuL6UEgC8gEngaaANOAEbELTaRujB7Std5nLJWXO68tyOOXs1awfV8x3zqlO/eMzKRD62b1GofI0Qg7BvE1YAiwEMDdN5pZ65hFJZLEFqzbwY/eXMaSvF0M7ZHBU9edysndMuIdlkithU0Qxe7uZuYAZtYyhjGJJKUtuw/w8FsreGNRPh2PacZvrxzMqMFdMLN4hyZyRMImiFeCWUwZZnYTkUqutVl+VCRlHSgpY+qHnzL5vVxKy53bzu3LrV/uR8tmKlQgyS3sIPX/mtkFwG4i4xAPuPs7MY1MJMG5O39bvoWf/vVj1hcUcuGAjtx/yQB6tG8R79BE6kTYQeqWwGx3f8fMMoFMM2vi7iWxDU8kMa3asoeH/rKcf676jOM7tmLa+OGc2f/YeIclUqfCtoH/AZxlZm2Bt4Fs4ErgqlgFJpKI1mzbyyOzc5m+OJ9WzRrzw8sGcPVpPWmcpmmrknpCLzkaPAsxHnjM3X9pZotjGJdIQlkdJIb/W5xP08aNuGFEb249tx/tWmqNBkldtVmT+nQiLYbxwb60as4XSQm5W/fyyOxVzPhoI80ap3HjWX246aw+Sfs8QzKUPpfEETZB3AXcB7zh7svMrA/wXsyiEomz3K17+UOQGJo3TuOms/pw09l9OLZVciYGSI7S55JYws5i+gD4wMxaBNtrgDtiGZhIPORu3cPv383lzSUbSW+SxoSz+zDhrD60T+LEUKG60udKEBJN2FlMpxNZA6IV0MPMBgE3u/utsQxOpL58smUPv393FX/N2UR6kzRuOacvN53VJ6XGGBKx9LkktrBdTL8FRgIzANz9IzM7O1ZBidSXT7bs4XfvrmJmziZaNEnj/53TlxtTLDFUSLTS55L4Qj/q6e4bqpQMKDvcuSKJbuXmSIth5tJIYrj1y3258cw+tE3BxFBh4sjMQ8YgID6lzyV5hE0QG8zsDMDNrAlwJ/Bx7MISiY0Vm3dHEkPOZlo1a8xtX+7H+DN7p3RiqBDv0ueSfMImiFuA3wFdgXzgb8BtsQpKpK59vCmSGN5aupnWzRpzx3n9uOHM3mS0SP3EUFk8Sp9L8go7i+kzjuCpaTO7iEhiSQOedPeHqxxvBjwLnAJsB65097Vm1otIC6ViZfk57n5Lbb9fZNnGXfz+3VXMWrYlkhjO78/4Eb1p00LrP4vUJOwspt7Ad4Bela9x98uruSYNmExkHes8YL6ZzXD35ZVOGw/scPd+ZjYG+AWREh4Aq919cPgfRSTC3Vm0YSd//GB1JDE0b8yd5/fnBiUGkVoJ28U0ncg01zeB8pDXDANyg2cmMLOXgFFA5QQxCvhh8Pk14BFT8Xw5QoXFpcxYvJHn5qxj2cbdtG7emLu+0p/rR/SmTboSg0hthU0QB9z997W8d1dgQ6XtPGD44c5x91Iz2wW0D471NrNFREqM3+/u/4z2JWY2AZgA0KNHj1qGKKkgd+teps1Zx58X5rHnQClf6tSan37tJEYP7qo1GUSOQtj/en4XrEv9N6CoYqe7L4xJVLAJ6OHu283sFGC6mZ3o7rurnujuU4ApAFlZWR6jeCTBlJSV8/flW3huzjr+vXo7TdMacfHATow7rSen9GyrVdxE6kDYBDEQGAecx+ddTB5sH04+0L3SdrdgX7Rz8sysMdAG2O7uTpCI3H2Bma0GjidSZlwasM27DvDivPW8NH89W3YX0TUjne9dlMm3srondZ0kkUQUNkF8E+jj7sW1uPd8oH8wwJ0PjAG+XeWcGcC1wH+AbxBZlMjNrANQ4O5lQWHA/sCaWny3pBB35z+rt/PcnHX8bfkWyt055/gO/OxrPfly5nGkNVJrQSQWwiaIpUAGsDXsjYMxhduBWUSmuT4VVIJ9CMh29xlEBr6fM7NcoIBIEgE4G3jIzEqItFhucfeCsN8tqWHX/hJeX5jHtDnrWL1tHxktmnDjmb359vAe9GzfMt7hiaQ8i/Tm1HCS2fvAyURaBZXHIA47zTUesrKyPDu74fVCJVqN/6ONZ2n+Lp6fu47pizayv6SMwd0zGHdaTy45uTPNmxy6DMmRfFei/b5E4s3MFrh7VtX9YVsQD9ZxPFJHEq3G/5HGc6CkjJk5m3huzjoWrd9J8yaNGD24K1ef1pOTuraps+9KtN+XSCKrzXoQkoASrcZ/beNZv72Q5+et45X5G9hRWEKfY1vywKUD+PrQbjU+1HYkP3ui/b5EEpkmiSe5RKvxHyaeHfuKeWf5Ft5cspEPcz+jkRkXnNCRcaf35Iy+7UNPUT2Snz3Rfl8iiUwJIsklWo3/w8XT8ZjmvDhvPTNzNvHv1dspK3e6tU3nO+f2Y+zwHnRuU/t4j+RnT7Tfl0giaxTvAOToTByZSXqVgdt41viPFk8jgy17DnDf6zmsLyhkwtl9ePP2M/nn987l7gszjyg5HO67avrZE+33JZLIqm1BmNlF7v528LkN8GvgVCLTXr/r7ltiH6JUJ9Fq/J/Rtz1fHdiJN5dsorg08kxl+1bNuDKrOxcP7MSAzsfU2VPOR/KzJ9rvSySRVTvN1cwWuvvQ4POTwGbgCeAK4Bx3H10fQYbVUKe5xtumXft5e+lm3srZzPx1BbhDv+Na8dWBnfnqwE5kdmyt0hciCexop7kCZFUqv/0bM7u2TiKTpLRx535m5mziraWbWbBuBwCZHVtz5/n9uWRgZ/p3bB3nCEXkaNWUII4zs7sBA44xM/PPmxwav2hgNhQU8tbSTczM2cziDTsBOKHzMfzXBcdz8cDO9DuuVXwDFJE6VVOCeAKo+KfgM8CxwDYz6wQsjmFckiDWbd/HzJzNvLV0E0vydgFwUtdjmDgyk68O7EzvY1XyQiRVVZsg3P1Hh9m/GbgmJhFJXG0oKGTupwXMXbOduZ8WsL6gEIBB3dpw78Vf4uKTOqkOkkgDccTPQZjZ9e7+dF0GI/XL3Vm7vfBgMpj3acHBZwTapDdhWO92XHtGLy4c0JHu7VrEOVoRqW9H86DcjwAliCTi7qzetpc5awoOthK27onUXmzfsinD+7Rjwtl9GNa7HZkdW9NIZbRFGrSanoNYcrhDQMe6D0fqUnm5s3LLHuau2c68tZEWwmd7I0t6HNe6GcP7tGd473ac1qcdfTu00lRUETlETS2IjsBIYEeV/Qb8OyYRJYB4loO+f3oOL87dQJk7aWaMHd6dn4weGOrasnLn4027mRN0Gc1fW8DOwhIAumakc3b/Dgzv047hvdvTs30LJQQRqVZNCeIvQCt3X1z1QLBGRMqJZzno+6fnMG3O+oPbZe4Ht6MliYJ9xazasofFG3YeTAh7DpQC0LN9Cy4c0JFhvSOtBI0hiEht1TSLaXw1x6ouH5oS4lkO+sW5G6Luf2HOekae2IlVW/aSu20vucF7wb7PV4Dt06Ell57chdP6tGNY73ZHXN9IRKSCqrlWEa9y0GXlTtlhyp6UA+OmzgMgo0UT+nVoxcgTO9K3Qyv6d2zNgM7H0KF1s5jGJyINjxJEFbEuB11UWsbazwrJ3bqX3K17WbV1D7lb97Lms32HvcaAF246jf4dW9G+ZVONHYhIvVCCqGLiyMxDxiAgXDlod2dPUSk79hWzo7CEHfuKKdhXzI7CYrbtLWLNtn2s3rqXdQWFlJVHWgpm0L1tC/od14pzju/Akryd/GdNwRfufdVpPTi9b/u6/UFFRGqgBFHF6CFdcXd+OWslm3YdoEOrZlwxtCuOM/XDTyN/+AuL2VkYJIB9JRQUFrNjXzGl5dG7iJqkGb3atySzU2suOTlSs6jfca3o26EVzausTXA0s5hEROpSteW+k82Rlvt+7j9reWvp5oP/4t+xr4TisvKo5zYyaNeyKRktmtKuRVPatmxCu5ZNadui6aHvLZvStkUT2rZsSutmjdUtJCIJqy7KfaesfcVlFJeW071dCwZ1yyCjZZPgj3/Tz9+Dz62bN9YTxiLSIChBALec05dbzukb7zBERBKK1nQQEZGolCBERCQqJQgREYlKCUJERKJSghARkahiniDM7CIzW2lmuWZ2b5Tjzczs5eD4XDPrVenYfcH+lWY2MhbxTV+Uz4iHZ9P73r8y4uHZTF+UH4uvOWKJHp+IpK6YTnM1szRgMnABkAfMN7MZ7r680mnjgR3u3s/MxgC/AK40swHAGOBEoAvwdzM73t0PLbV6FOJZ2juMRI9PRFJbrFsQw4Bcd1/j7sXAS8CoKueMAp4JPr8GnG+Rx45HAS+5e5G7fwrkBverM9WV9k4EiR6fiKS2WCeIrkDlRQ7ygn1Rz3H3UmAX0D7ktZjZBDPLNrPsbdu21Sq4eJX2DivR4xOR1Jb0g9TuPsXds9w9q0OHDrW69nAlvOuqtPfRSvT4RCS1xTpB5APdK213C/ZFPcfMGgNtgO0hrz0qE0dmkl6lmmqY0t71JdHjE5HUFusEMR/ob2a9zawpkUHnGVXOmQFcG3z+BjDbIyVmZwBjgllOvYH+wLy6DG70kK78/IqBdM1Ix4CuGen8/IqBCTMAnOjxiUhqi3m5bzP7KvBbIA14yt1/amYPAdnuPsPMmgPPAUOAAmCMu68Jrv0+cANQCtzl7m9V911HWu5bRKQhO1y5b60HISLSwB0uQST9ILWIiMSGEoSIiESlBCEiIlEpQYiISFQpNUhtZtuAdVEOHQt8Vs/hHC3FHHvJFi8o5vrS0GLu6e5feNI4pRLE4ZhZdrQR+kSmmGMv2eIFxVxfFHOEuphERCQqJQgREYmqoSSIKfEO4Ago5thLtnhBMdcXxUwDGYMQEZHaaygtCBERqSUlCBERiSrlE4SZXWRmK80s18zujXc8NTGzp8xsq5ktjXcsYZhZdzN7z8yWm9kyM7sz3jHVxMyam9k8M/soiPlH8Y4pLDNLM7NFZvaXeMcShpmtNbMcM1tsZglfSdPMMszsNTNbYWYfm9np8Y6pOmaWGfxuK167zeyuOrt/Ko9BmFka8AlwAZElS+cDY919eVwDq4aZnQ3sBZ5195PiHU9NzKwz0NndF5pZa2ABMDrBf8cGtHT3vWbWBPgQuNPd58Q5tBqZ2d1AFnCMu18a73hqYmZrgSx3T4qHzszsGeCf7v5ksIZNC3ffGeewQgn+3uUDw9092gPDtZbqLYhhQK67r3H3YuAlYFScY6qWu/+DyLoYScHdN7n7wuDzHuBjoqwdnkg8Ym+w2SR4Jfy/lMysG3AJ8GS8Y0lFZtYGOBuYCuDuxcmSHALnA6vrKjlA6ieIrsCGStt5JPgfr2RmZr2ILPw0N86h1CjoqlkMbAXecfeEj5nIwlvfA8rjHEdtOPA3M1tgZhPiHUwNegPbgKeDbrwnzaxlvIOqhTHAi3V5w1RPEFJPzKwV8GciK//tjnc8NXH3MncfTGSt82FmltDdeWZ2KbDV3RfEO5ZaOtPdhwIXA7cFXaiJqjEwFHjM3YcA+4CEH7cECLrDLgdercv7pnqCyAe6V9ruFuyTOhT04/8ZeN7dX493PLURdCG8B1wU51BqMgK4POjTfwk4z8ymxTekmrl7fvC+FXiDSLdvosoD8iq1Jl8jkjCSwcXAQnffUpc3TfUEMR/ob2a9gww7BpgR55hSSjDgOxX42N1/He94wjCzDmaWEXxOJzKJYUVcg6qBu9/n7t3cvReR/x/Pdver4xxWtcysZTBxgaCr5kIgYWfnuftmYIOZZQa7zgcSdrJFFWOp4+4liDSpUpa7l5rZ7cAsIA14yt2XxTmsapnZi8CXgWPNLA940N2nxjeqao0AxgE5QZ8+wP+4+8z4hVSjzsAzwayPRsAr7p4U00aTTEfgjci/IWgMvODub8c3pBp9B3g++AflGuD6OMdToyD5XgDcXOf3TuVpriIicuRSvYtJRESOkBKEiIhEpQQhIiJRKUGIiEhUShAiIhKVEoRIHTCzH5rZPbW85vJkqDAsDVdKPwchksjcfQZ6cFMSmFoQIkfIzL5vZp+Y2YdAZrBvsJnNMbMlZvaGmbUN9t8RrJmxxMxeCvZdZ2aPBJ+/aWZLgzUq/hG3H0qkErUgRI6AmZ1CpOTFYCL/HS0kshbGs8B33P0DM3sIeBC4i0jRt97uXlRR5qOKB4CR7p5/mOMi9U4tCJEjcxbwhrsXBtVrZwAtgQx3/yA45xki6wsALCFSwuFqoDTK/f4F/MnMbiJSFkYk7pQgROrHJcBkItVB55vZIa13d78FuJ9I9eEFZta+/kMUOZQShMiR+Qcw2szSg4qllxFZP2CHmZ0VnDMO+MDMGgHd3f094L+BNkCryjczs77uPtfdHyCyaE3lMvUicaExCJEjEKzB/TLwEZFV6eYHh64FHjezFnxeDTQNmBYsaWnA7919Z1DltMIkM+sfHH83uK9IXKmaq4iIRKUuJhERiUoJQkREolKCEBGRqJQgREQkKiUIERGJSglCRESiUoIQEZGo/j+SchBagLYHeQAAAABJRU5ErkJggg==\n", + "text/plain": [ + "<Figure size 432x288 with 1 Axes>" + ] + }, + "metadata": { + "needs_background": "light" + }, + "output_type": "display_data" + }, + { + "name": "stdout", + "output_type": "stream", + "text": [ + " 2\n", + "0.004728 x + 0.0002015 x + 0.01276\n" + ] + } + ], + "source": [ + "x = [7,6.2,4.4,4.3,3,2.9,2.7,2.1,1.3,1.1,1,1,0.3,0.2]\n", + "y= [1.253,1.030,0.824,0.540,0.270,0.262,0.228,0.160,0.070,0.060,0.053,0.050,0.02,0.005]\n", + "z= [0.26,0.18,0.06,0.14,0.04,0.08,0.04,0.04,0.04,0,0.02,0.02,0.02,0]\n", + "\n", + "poly = numpy.polyfit(x, y, 2)\n", + "polynomial = numpy.poly1d(poly)\n", + "plt.scatter(x, y)\n", + "\n", + "xs = numpy.linspace(x[0], x[-1], len(x))\n", + "y_int = polynomial(xs)\n", + "\n", + "plt.xlabel(\"dosis\")\n", + "plt.ylabel(\"Inmediatamente despues del accidente\")\n", + "plt.plot(xs, y_int)\n", + "plt.show()\n", + "print(polynomial)\n", + "\n", + "\n", + "\n", + "\n", + "poly = numpy.polyfit(x, z, 2)\n", + "polynomial = numpy.poly1d(poly)\n", + "plt.scatter(x, z)\n", + "\n", + "xs = numpy.linspace(x[0], x[-1], len(x))\n", + "z_int = polynomial(xs)\n", + "\n", + "plt.xlabel(\"dosis\")\n", + "plt.ylabel(\"15 meses después del accidente\")\n", + "plt.plot(xs, z_int)\n", + "plt.show()\n", + "print(polynomial)" + ] + }, + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + " Entonces el modelo indica que la frecuencia del dicéntrico aumenta efectivamente con el cuadrado de la dosis inmediatamente. " + ] + } + ], + "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.3" + } + }, + "nbformat": 4, + "nbformat_minor": 4 +}