{ "cells": [ { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import matplotlib.pyplot as plt\n", "# for inline plots in jupyter\n", "%matplotlib inline\n", "import numpy as np\n", "import math" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "# import seaborn\n", "import seaborn as sns\n", "# settings for seaborn plotting style\n", "sns.set(color_codes=True)\n", "# settings for seaborn plot sizes\n", "sns.set(rc={'figure.figsize':(5,5)})" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "from scipy.stats import beta\n", "infest_beta = beta.rvs(4, 60, size=10000)" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 4, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAVwAAAFtCAYAAABLMNIkAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/Il7ecAAAACXBIWXMAAAsTAAALEwEAmpwYAABAo0lEQVR4nO3de3wU1f3/8dfsPZtks0kIkIQ7SIhSJRCKN2wFFFTU1mqhVG21Wmtba7Vi+UoNFrUIUq0oitW2tv1S+eqvFbxQQCpKtdWCxSIFQSBcE5OQZHPd68z8/ghZCYTNJtl7Ps/Hg4dmZy9ndrLvnD3zmXMUXdd1hBBCRJ0h3g0QQoi+QgJXCCFiRAJXCCFiRAJXCCFiRAJXCCFiRAJXCCFixBTvBkRbbW0zmtax8i072059fWucWhR7sr+pTfY3seTlZZ52W5/s4ZpMxng3IaZkf1Ob7G/y6JOBK4QQ8SCBK4QQMSKBK4QQMRKzk2bf//73OXLkCAaDAbvdzv33309xcTFTpkzBYrFgtVoBuOeee5g8eTIAH330EWVlZXi9XgoLC3n00UfJzc2NVZOFECKiYha4ixcvJjOz7ezdxo0bue+++3jllVcAWLZsGaNHj+5wf03TmDt3LosWLaK0tJSnn36apUuXsmjRolg1WQghIipmQwrtYQvQ3NyMoigh779jxw6sViulpaUAzJ49m3Xr1kW1jUIIEU0xrcOdP38+7733Hrqu8/zzzwdvv+eee9B1nQkTJnD33XfjcDiorKykoKAgeJ+cnBw0TcPlcuF0OmPZ7IQXqu4v3PvU1DRFqjlCiNNQ4jEf7urVq3njjTd47rnnqKysJD8/H5/Px8MPP0xLSwtLly5l/fr1/PnPf+bXv/518HHnnHMO77zzjgSuECIpxeVKs6985SuUlZVRX19Pfn4+ABaLhTlz5nD77bcDkJ+fT0VFRfAxdXV1GAyGbodtZ1ea5eVlplSPLpwebldS7f1Ipf3piuxvYon7lWYtLS1UVlYGf37rrbfIysrCarXS1NT2xum6ztq1aykuLgZg7NixeDwetm7dCsCqVauYMWNGLJorhBBREZMertvt5s4778TtdmMwGMjKymLFihXU1tZyxx13oKoqmqYxcuRIFixYAIDBYGDJkiUsWLCgQ1mYEEIkq7iM4caSDCmEJ9Xej1Tan67I/iaWuA8pCCGEkMAVQoiYkcAVQogYkcAVQogYkcAVQogYkcAVQogYkcAVQogYkcAVQogYkcAVQogYkcAVQogYkcAVQogYicv0jKLnIjFvghAiPqSHK4QQMSI9XAGcvuecyLMyCZFspIcrhBAxIoErhBAxIoErhBAxIoErhBAxIoErhBAxIoErhBAxIoErhBAxIoErhBAxIoErhBAxIleaJZhIzpWg6zp+re3/LUYlYs8rhOgZCdwUVdUaYP3hZipaAyjAhDwbk/PtWI3ypUaIeJHATUFHW/y8+GkDFqPC5Hw7DT6VrTUeatwqXx/lwKhIb1eIeJDATTENPpX/t7+RDLOBG0Y7STe39WgHp5t541Azm462MG1QRpxbKUTfJN8vU4iu66w/3ExA0/n6yKxg2AJ8IdfG+H42ttZ4qGoNxLGVQvRdErgpZJfLx/5GPxflp5NjM56y/aJ8O2lGhbeOtqDrehxaKETfJoGbIgKazttHWxiQZmRCnq3T+9hMBi4YaOdgs58jLdLLFSLWJHBTxEfHPDT6Nb5ckI4hxEmxc/rZsBkVtlS7Y9g6IQRI4KaEgKbzz6pWhmSYGZZpDnlfs0GhpJ+NPQ0+6r1qjFoohAAJ3JTw3zovLQGd8wemoYRR8jU+z4YC/KfWE/3GCSGCJHCTnK7rfFDtZkCakaEZoXu37TLNRkY4zOyo86LJyTMhYkYCN8ntbfRR51WZ1N8eVu+23RdybDT7NQ42+aPYOiHEiWJ24cP3v/99jhw5gsFgwG63c//991NcXEx5eTnz5s3D5XLhdDpZvHgxw4YNAwi5TbT5oMqNw2xgTLalW48blWXBZlT4uM7LcEf3HiuE6JmY9XAXL17Mq6++yurVq7n55pu57777AFiwYAFz5sxh/fr1zJkzh7KysuBjQm0TUNHSVt41sX9ayMqEzpgMCkVOC3sbfAQ0GVYQIhZiFriZmZ/PgtXc3IyiKNTW1rJz505mzpwJwMyZM9m5cyd1dXUht4k2W2s8WI0KZ+dae/T40VlWfJouwwpCxEhM51KYP38+7733Hrqu8/zzz1NZWcmAAQMwGtuuijIajfTv35/Kykp0XT/ttpycnLBfMze383kDIjkNYjy0+DU+cXkZ38/W4xnAhmaasRoUdru8jMzqfFghWd+nZG13T8n+JoeYBu7DDz8MwOrVq1myZAl33nln1F+ztrYZ7aSvzHl5mdTUNEX9tXsi3F+k7bUeNB1K+nV+VVk4TAaFkVkWPm3woel62MMSifretUvk4xsNsr+JJdRnOC5VCl/5ylf44IMPGDhwIFVVVahqWwG+qqpUV1eTn59Pfn7+abf1dZqus+2Yh6EZZnJtvfubOdppwa3qHG6WYQUhoi0mgdvS0kJlZWXw57feeousrCxyc3MpLi7m9ddfB+D111+nuLiYnJyckNv6uv2Nfhr9Wq96t+1GZFowKbDb5YtAy4QQocRkSMHtdnPnnXfidrsxGAxkZWWxYsUKFEXhgQceYN68eTz99NM4HA4WL14cfFyobX3Zv4+5yTAZOMPZ+3Iui1FhhMPCngYflwzSw6rlbf/KlMhf64RIRDEJ3H79+vHSSy91um3kyJG8/PLL3d7WVzX4VPY3+jl/YFrEVm4Y7WwL3IrWAIXp4V2tJoToPlnxIQF054zrjjovAGfn9H44od0ohwUF2Nfgk8AVIork0t4kous6O+o8DMkw47SeOsF4T9lMBgrTTexvlBNnQkSTBG4SqWgNUO/VGJvTswsdQhnhsPCZO0BL+7rqQoiIk8BNIjvqvJgUKIrAybKTjTg+n8L+RqlWECJaJHCTREDT2Vnvpchp7fGVZaEMSDOSblIkcIWIIgncJLG3wYdX1aMynACgKG3lYeVNfpkjV4gokcCNo7y8zLArFHbUe8kwGxjaxRI6vTHCYcGj6lTIApNCRIUEbhLwqhrljT7GOC3dnoaxO4ZnmlGQcVwhokUCNwnsa/Sj6lDkjM5wQjspDxMiuiRwk8Ael5d0k0JhevSvU5HyMCGiRwI3wQU0nX2NPkY7rVEdTmgn5WFCRI8EboI71OzHr7VdfhsLUh4mRPRI4Ca4fY0+TAoMiWJ1womkPEyI6JHATXD7G30MzTRjNkR/OKGdlIcJER0SuAmszqNS79WC46qx0l4eVt4kwwpCRJIEbgJrD7xYB67NZKBAysOEiDgJ3AR2sMmPw2zAaYn9YRqeaaGyNYA7IOVhQkSKBG6C0nWdQ81+hmSaw1r2JtJGONpO0pU3SS9XiEiRwE1Q1W4Vj6ozNCM+KzAMtJuwGRXKpTxMiIiRJXYS1KHjy5bHqhzsZAZFYXimmf2NPnS988UlT5x4RxaUFKJrErgxFu7sYIea/TgtBrIskVtKp7uGOyzscvmodqsMsMuvihC9JUMKCUjXdY62+OO+oOPw4DiuDCsIEQkSuAmowafRGtBjMllNKJlmI3k2o5SHCREhErgJ6GhLW8DFu4cLbTXAR1r8+FS5zFeI3pLATUBHWwJYDAp5afEbv2033GFG0+FgswwrCNFbErgJ6GiLn3y7KSbTMXZlULoZswHKZVhBiF6TwE0wfk2n2q1SEOfx23Ymg8KQDLNM1yhEBEjgJpgadwCdtgsPEsUIhwWXT6Peq572Pu0LYoZb9iZEXySBm2A+a22bEjHRAhdkFQgheksCN8F81hogzaTgMCfOocm2GnFaDBK4QvRS4nyqBQCfuQMMTDPFZcKaUEY4LBxq9hPQpDxMiJ6SwE0gAU3nmFtNqOGEdsMdZvwaHGmRagUheirxPtkpKNwTSTWeABok5LwFQzMsGBTY3+hnWGZsJ0QXIlVIDzeBVLe2VQEMSEu8wLUYFYZmmNnbIOO4QvRUTD7Z9fX13HvvvRw6dAiLxcLQoUNZuHAhOTk5FBUVMXr0aAyGtuxfsmQJRUVFALz11lssWbIEVVU566yzWLRoEWlpabFoclxUewKYDcRlhYdwjMqy8OaRFuo8Kjm2+F8FJ0SyicknW1EUbrnlFtavX89rr73G4MGDWbp0aXD7qlWrWLNmDWvWrAmGbUtLC/fffz8rVqzgzTffJD09nd/85jexaG7c1LhV8myJd8Ks3aistqGETxu8cW6JEMkpJoHrdDqZNGlS8Odx48ZRUVER8jGbN29m7NixDBs2DIDZs2fz17/+NZrNjCtd16lxBxJi/oTTybK0zR62V8rDhOiRmA8WaprGiy++yJQpU4K33XDDDaiqykUXXcQdd9yBxWKhsrKSgoKC4H0KCgqorKyMdXNjpiWg41Z18myJN357olFZFt6vcuMOaKSZOv973dlJQlkRQog4BO6DDz6I3W7n+uuvB+Dtt98mPz+f5uZm5s6dy/Lly7nrrrsi9nq5uRmd3p5ol6DWuNuuMEvkHi60Be4/q9zsb/RxVo4t7MfF+v1OtOMbbbK/ySGmgbt48WIOHjzIihUrgifJ8vPzAcjIyOC6667jd7/7XfD2Dz74IPjYioqK4H27o7a2Ge2kYv28vMyY9rjC+eWoDgZuYvdwC+wm7CaFvQ3dC9xYv999qUct+5tYQn3eY3Y6/LHHHmPHjh0sX74ci6Xt5EtDQwMejweAQCDA+vXrKS4uBmDy5Ml8/PHHHDhwAGg7sXbZZZfFqrkxd8yjkm5SsJ/ma3qiUBSFUQ4L+5v8qHLVmRDdEpPu1Keffsqzzz7LsGHDmD17NgCDBg3illtuoaysDEVRCAQClJSUcOeddwJtPd6FCxdy2223oWkaxcXFzJ8/PxbNjYtjHpV+CT5+2+4Mp4XtdV4ONfsZ7pCLIIQIV0w+4WeccQa7d+/udNtrr7122sdNmzaNadOmRatZCUPXdeo8KmflWOPdlLAMz7RgMSh84vJK4ArRDcnRpUpxzQENr6aTmyQXE5gMCqOyLOxx+Zg+WA9rZYquxrETeUxOiEhJ7AHDPqLW03ZJb7IELkCR04Jb1TnULJPZCBEuCdwEkIyBO8JhwWyA3S65CEKIcEngJoBaj4rVoJCR4BUKJzIbFEY6LOx2edF0qVYQIhzJ8wlPYcc8Krk2Y8LOoXA6RU4rrQGdI82BeDdFiKQggZsA6o4HbrIZ6bBgUuATl0xmI0Q4JHDjzKfqNAc0cqzJF7gWo8IIR1u1gi7DCkJ0SQI3isJZNrx96fHsJAxcaKtWaA5oHG2RYQUhuiKBG2fJHrijsiwYZVhBiLBI4MZZe+A6rcl5KKxGA8MzLeyWYQUhupScn/IUUu9rm7TGakzeQ1HktNDk16hslWEFIUJJ3k95iqj3qkk7nNDujCwLBuATuQhCiJAkcOOs3qslfeDaTAaGZprZ7fLKsIIQIUjgxpFP1Wn2J3/gQttFEA0+jSq3Gu+mCJGwJHDjyOVL7gqFE43OsqAAe6RaQYjTksCNo2QvCTuR3WxgcIZZJrMRIgQJ3Dj6PHBT4zAUOS3UelWOuaVaQYjOpMYnPUnVe1XsSV4SdqLRzrbVH3Y3SC9XiM6kxic9SdV7k3MOhdPJNBspTDfxqQwrCNEpCdw4qveqOFMocAFGOSx85g7Q5JdqBSFOJoEbJ35Np8mfWj1cgJFZbcMK+xtk6R0hTiaBGyepVKFwojybEYfZwN5GGVYQ4mQSuHGSqoGrKG0r+h5o8hHQ5KozIU4kgRsnriSfJSyUUVkW/BocbJJhBSFOlHqf9iRR79VIMyrYUqQk7ERDMsyYDbBPhhWE6CD1Pu1JosGnkpViwwntTAaFYZkW9jbIHLlCnEgCN04afBpOS+q+/aMcFhr9GjUeKQ8Tol3qfuITmK7rbT1cS2r2cAFGZJkB2CdXnQkRFHbgbty4kUBArpGPhOaAhqpDVgr3cDPNRvqnGdnfJIErRLuwP/HLli3jwgsvZOHChfznP/+JZptSXoNXA0i5q8xONtJh4UhzAI+qxbspQiSEsAP31Vdf5YUXXsBqtXLHHXcwffp0nn76aY4cORLN9qWkhuPz4KZyDxdghMOCDhxolPIwIaCbY7hjxozhpz/9Ke+88w4LFixg3bp1XHLJJXzzm9/k1VdfRdOkJxMOl6/tfXKk8BguQGG6CatRYb+UhwkBgKm7Dzh06BCvvvoqr776Koqi8KMf/Yj8/HxWrlzJhg0beOqpp6LRzpTScHylXrNBiXdTosqgKAzPNLO/0Y+u6yhKau+vEF0JO3BXrlzJmjVrOHjwIJdddhlLlixh3Lhxwe3Tp0/n/PPPj0YbU06DV0v58dt2IxwWPnH5qHarDLB3+++7ECkl7E/A5s2buemmm5g6dSoWi+WU7WlpaTz55JMRbVyqavCpFKSb492MmBjhOD57WKNPAlf0eWF/ApYtW4bBYMBs/jwo/P62r4rtAXzhhRd2+tj6+nruvfdeDh06hMViYejQoSxcuJCcnBw++ugjysrK8Hq9FBYW8uijj5KbmwsQcluy0nSdRp9GcXZqnzBrl2E2MCDNyL5GH+cNtMe7OULEVdif+ptvvpn//ve/HW7773//y3e+850uH6soCrfccgvr16/ntddeY/DgwSxduhRN05g7dy5lZWWsX7+e0tJSli5dChByWzJr8mtokNIXPZxshMPC0ZYAnsDpT6rm5WWG9U+IZBZ24O7evZtzzjmnw21nn302n3zySZePdTqdTJo0KfjzuHHjqKioYMeOHVitVkpLSwGYPXs269atAwi5LZkFa3BTvCTsRCPby8Nk9jDRx4X9qXc4HBw7dqzDbceOHSMtLa1bL6hpGi+++CJTpkyhsrKSgoKC4LacnBw0TcPlcoXclsxc7TW4feSkGUDB8fIwmT1M9HVhj+Feeuml/OQnP+FnP/sZgwcP5tChQzzyyCNcdtll3XrBBx98ELvdzvXXX8+bb77Z7QZ3V25uRqe3x+vraftFDw5z3+nhfl4e5ut1eVi4x62vDT/I/iaHsAP3rrvu4pFHHuG6667D5/NhtVq55ppruPvuu8N+scWLF3Pw4EFWrFiBwWAgPz+fioqK4Pa6ujoMBgNOpzPktu6orW1GO2nlgby8TGpqmrr1POHq6hehwafhMBswpngN7slGHi8Pq3KrDOxFtUI4xy2axzcRyf4mllAZEPZvvtVqZcGCBZSVlVFfX092dna3eiqPPfYYO3bs4Ne//nWwqmHs2LF4PB62bt1KaWkpq1atYsaMGV1uS2YNPhVHHxq/bXdieVhvAleIZNat3/ympibKy8tpaWnpcPt5550X8nGffvopzz77LMOGDWP27NkADBo0iOXLl7NkyRIWLFjQofQLwGAwnHZbMmvwagzJ7Bs1uCdKP14etr/Rx/lSHib6qLAD9y9/+QsLFy7Ebrdjs9mCtyuKwt/+9reQjz3jjDPYvXt3p9vGjx/Pa6+91u1tyUjVdBr9WspPWnM6Ix0W/lnlxhPQsJn65nsg+rawA/fxxx/niSee4Etf+lI025PSGv1tJWF9qQb3RCMcFv5R5aa8yU9xtjXezREi5sLuZqiqetoryUR4Grx9Y1rG0ylIN2GT2cNEHxb2J//WW2/lmWeekSkYe6HB1zcmHj+dk8vDhOhrwh5SeOGFFzh27BjPP//8KaVZb7/9doSblZpcPhUFyOxDNbgnG+GwsCsC5WFCJKOwf+NToUIg3hp8Gg6LAUMfnhdWysNEXxb2b/wXv/jFaLajT0j1lXrDkW42MNBuYp+Uh4k+KOzvtj6fj8cff5ypU6cyYcIEAN59913+93//N2qNSzUur9qnJq05nTOy2mYPa/Kr8W6KEDEV9qf/F7/4BXv27GHp0qXBK8zOOOMMXnzxxag1LpX4NZ2WgN6nJq05nSJn27DCHpdUK4i+JewhhY0bN7JhwwbsdjsGQ1tODxgwgKqqqqg1LpU09pGVesPRz2Yi12Zkt8vHhLzuzTYnRDIL+9NvNptR1Y5fAevq6ro9mUxfFSwJ6+NjuO2KnBYON/tp8UuZoeg7wg7cGTNm8NOf/pTDhw8DUF1dzcKFC7niiiui1rhU4urjFz2crCjLig582iDDCqLvCPvTf9dddzFo0CCuuuoqGhsbmT59Ov379+cHP/hBNNuXMhp8GkalbY0vAf3TjGRbDex2eePdFCFiJuwxXIvFwn333cd9991HXV1dt6dn7Ovap2WU96yNoigUZVn5V7Ubd0AjTSazEX1A2IHbPpTQ7sQpGgcPHhy5FqWoBp/W52twT1aUbeH9ajefNvg4O9fW9QOESHJhB+4ll1yCoigdroFv763t2rUr8i1LMQ0+lTOyLPFuRkIZmGbCYWkbVpDAFX1B2IF78uq8NTU1PPXUU8FVdcXp+TWd1oAuPdyTKIrCGKeVrTUyR67oG3r8G56Xl8f8+fN57LHHItmelNQgNbinNcZpQdOlWkH0Db1KgP379+N2uyPVlpTV4O3bE4+Hkm834TAb2C1XnYk+IOwhhTlz5nQ4w+52u9m7d6+UhYUh2MO1Sg/3ZIqiUOS08O9jHjyqhs0o75FIXWEH7nXXXdfh57S0NMaMGcOwYcMi3aaUE6zBlTHKTo3JtrKlxsPeBh9jc+TkmUhdYQfuV7/61Wi2I6VJDW5oBXYTmceHFboK3Ly8zOD/19Q0RbtpQkRU2IH7xBNPhHW/O++8s8eNSVVSgxta+7DCtmMevKqGVYYVRIoKO3APHjzIhg0bGDt2LIWFhVRUVPDxxx9z6aWXYrXKCqyhSA1u19rKwzzsa/BzZo78PonUFHbg6rrOL3/5S6ZPnx68bcOGDaxbt45FixZFpXHJ5MSvuidqr8GVWcJCK0w3kWE28InLG3bgtr/nMrQgkkXY3902b97MtGnTOtw2ZcoU3nnnnYg3KpV8XoMrgRtK+7DC/kYfPlVW9BWpKezAHTp0KCtXruxw24svvsiQIUMi3qhUEqzBlZKwLhU5rQR02NcoNbkiNYU9pPDQQw/xwx/+kOeffz640oPJZOLJJ5+MZvuSnvRwwzco3USaSWFvg4/ibBnHFakn7MA988wzWb9+Pf/5z3+orq4mLy+PcePGYTabo9m+pNdeg5tukpKwrhgUhZEOC3sbfGi63qeXkxepqcffcydOnIjf76e1tTWS7Uk57UujSw1ueEY5LHhUnSMtgXg3RYiIC7uHu3v3bm6//XYsFgtVVVVcfvnlbNmyhVdeeYVf/epXUWxicmurwZXx23ANd5gxKLC3wceQDPn2JFJL2EnwwAMP8KMf/Yh169ZhMrXl9MSJE/nwww+j1rhU0N7DFeGxGg0MyTCzV2YPEyko7MDdu3cvV199NfD5xON2ux2vV9akOh2f2j4PrvRwu2NUloU6r0qdR+36zkIkkbCToLCwkB07dnS4bfv27VIWFkKjVCj0yChH21V5e6U8TKSYsMdw77zzTm677TZmz56N3+/n2WefZdWqVTz44IPRbF9Sa/BJDW5POK1G8mxG9jb4+GL/tHg3R4iICTsJLr74Yp5//nnq6uqYOHEiR48e5cknn+TCCy+MZvuSmtTg9tzILAuHm/14Alq8myJExITVw1VVlenTp7N27VoeeOCBHr3Q4sWLWb9+PUePHuW1115j9OjRQNvlwRaLJTgBzj333MPkyZMB+OijjygrK8Pr9VJYWMijjz5Kbm5uj14/HqQGt+dGOiy8X+XmQJOfMXIRhEgRYfVwjUYjRqOxVyfIpk6dysqVKyksLDxl27Jly1izZg1r1qwJhq2macydO5eysjLWr19PaWkpS5cu7fHrx4PU4PZcYboJq1GRy3xFSgl7SOHGG2/kxz/+Mf/61784dOgQhw8fDv4LR2lpKfn5+WE3bMeOHVit1uCqwLNnz2bdunVhPz4RSA1uzxkUhRGZZvY3+tB1mcxGpIYuhxRqamrIy8sLnhz7xz/+0eEDoCgKu3bt6lUj7rnnHnRdZ8KECdx99904HA4qKyspKCgI3icnJwdN03C5XDidzrCfOzc3o9PbTzedYiQ1+FQGZMnX4Z4amWVhl8tHlVtloP30v6qdHctYHN9EIvubHLoM3OnTp/Pvf/+bTz75BIAf/OAHLF++PGINWLlyJfn5+fh8Ph5++GEWLlwY0aGD2tpmNK1jDykvLzPic6ie/AsgNbi9NzyzrTxsX6MvZOCefCyjcXwTmexvYgn1x6DLNDj569yWLVt636ITtA8zWCwW5syZw7///e/g7RUVFcH71dXVYTAYutW7jadgDa5VKhR6Kt1sIN9uYp9cdSZSRJeBe/IJn0iOp7W2ttLU1BR83rVr11JcXAzA2LFj8Xg8bN26FYBVq1YxY8aMiL12tAVrcKWH2ysjHRYqWgO0+qU8TCS/LocUVFXl/fffDwbtyT8DnHfeeV2+0EMPPcSGDRs4duwYN910E06nkxUrVnDHHXegqiqapjFy5EgWLFgAgMFgYMmSJSxYsKBDWViykBrcyBjpMPPuZ7C/SZZQF8lP0bvosk6ZMiX0EygKf/vb3yLaqEiK1xjupqMtfFjj5ifn5EpZWC/ous6TO+oYlmnhqmGdj43JGK7sbyIJNYbbZQ/3rbfeimhj+ooGn4pDanB7TVEURsik5CJFyABjlEgNbuSMPD4peYVMSi6SXNiT14hThfrqIDW4kTM804wC7G/0MUgmJRdJTLpgUSA1uJFlMxkoTDfJdI0i6UkiRIHU4EbeGVkWqt0qLq9MSi6SlwRuFLikBjfiRjvbhmc+lYsgRBKTRIgCqcGNvOzjk5LvdsmSTiJ5SeBGQYNPwyTz4EbcaKeFIy0BWuSqM5GkJHC7KS8vM/jvdKQGNzqKZFhBJDkpC4uCBq+GU9Yxi7g8mxGnxcBul5dx/T6/zLcn0zMm8pVKInVJKkSBy6filPHbiFMUhSKnlYOy1plIUhK4EeYJaHhUHaeUhEXFaKcFTZcl1EVyksCNsPaSMKeUhEVFgd2Ew2xgV71UK4jkI6kQYe2F+dLDjQ5FUTgrx8r+Rj/NUq0gkowEboS5gjW48tZGy9gcKzqwU3q5IslIKkSYy6uRZlKwGuWtjZZcm4l8u4kddZ54N0WIbpFUiDCpUIiNsTlWqt0qVa0yZaNIHhK4EebyqnLCLAbOzLZiUJBerkgqkgwRpOk6jT5NTpjFQJrJwCiHhZ31XrQILmwqRDRJ4EZQo09DQyoUYmVsjpWWgM4+qckVSUICN4LaKxRkSCE2RmZZyDAZ+OiYDCuI5CDJEEEN3uMXPUgPNyaMisIXcttqctsnfRcikUngRpDLp2IAMs3ytsbKObk2dGB7rdTkisQnyRBB9V6VLKtBlvKOIafVyPBMM9trPXLyTCQ8CdwIcvk0qcGNg3NybTT6Ncob/fFuihAhSeBGkMuryvhtHJyRZcFuUvioVk6eicQmgRshwWkZpUIh5owGhS/k2Njb4KPJLyfPROKSdIiQ4Eq90sONi/aTZx/LyTORwCRwI+TzGlwJ3HjIsRkZmmHmP7UedDl5JhKUBG6ENATnwZW3NF7O6WejwadxoElOnonEJItIhqmrRQldPg2bUcEm0zLGzegsC2nGtpNnwx2WkPc98XjKgpIiViQdIkQqFOLPZFAYm2PlU5ePFlkNQiQgCdwIaZsHV97OeBvXz4aGTNsoEpMkRARouk6DV6ZlTATtq0HI8jsiEcUkcBcvXsyUKVMoKipiz549wdvLy8uZNWsW06dPZ9asWRw4cCCsbYmmyX98WkapUEgIZ2ZbqXKrHPPIahAiscQkcKdOncrKlSspLCzscPuCBQuYM2cO69evZ86cOZSVlYW1LdG4pEIhoRRnW1GAnXXSyxWJJSYJUVpaSn5+fofbamtr2blzJzNnzgRg5syZ7Ny5k7q6upDbElH7RQ/Sw00MGWYDQzPN7Kz3Sk2uSChx65JVVlYyYMAAjMa2kDIajfTv35/KysqQ2xKRy6uiAA45aZYwzsy24vJpVMgikyKBpHwdbm5uRqe3d1VX2x0ur4rDItMyJpIip4X1h+G/dV4K080h7xvJ34V4SYV96I5k3d+4BW5+fj5VVVWoqorRaERVVaqrq8nPz0fX9dNu667a2mY0rePXyry8zG4Xu4c6wPVejRypUEgoVqOBUVkWPnF5mTYoPeQfw2S/8KEnv8/JLNH3N1RWxO07cG5uLsXFxbz++usAvP766xQXF5OTkxNyW6LRdZ16n0q2BG7CKXZaaQ3oHG6WS31FYlD0GJxVeOihh9iwYQPHjh0jOzsbp9PJG2+8wb59+5g3bx6NjY04HA4WL17MiBEjAEJu645o93Bb/RrLdtQxtTCdif3Tut0+ET0+VWfZx7WcnWvj0sGdDy2B9HCTTaLvb6gebkwCN56iHbhHW/z8cU8D145wMCor9PX7Ivb+sr+RipYAPxibjXKaYYVE/vCGI9EDKNISfX8TckghVdR52mpws6UGNyEVOS00BzSOtki1gog/SYleqve1lYRJDW5iGpllwaDAngZfvJsihARub7m8Gg6LAaNBSsISkc1oYHimmd0uuQhCxF/K1+FGW51XKhQS3WinlX2Hmqlyqwy0n/or39mYWyKPEYrkJT3cXtB1nXqvKjW4Ce6MLAsKsNslcyuI+JLA7QW3quNVdZmWMcHZTQaGZJjZ7fKFPayQl5cZ/CdEpEjg9kK9VyoUkkWR00KdV6XWI8uoi/iRpOiF9sCVIYXEN9ppBWC3VCuIOJLA7YX647OEZUlJWMLLMBsYlG6ScVwRVxK4vVB/vCTMJCVhSWG000q1Ww1+MxEi1iRwe6FeSsKSSpGz7dJr6eWKeJHA7QUJ3OSSZTGSbzfxiUvGcUV8SOD2kDug4VF1CdwkM8Zp4bPWQHAdOiFiSQK3h6QkLDkVHa9WkLkVRDzIpb0hhCp6b6/nlJKw5OK0GhmQZmS3y8sXZf5iEWPSPeuhOq+KAeQqsyRU5LRytCVAo0+GFURsSeD2UK2n7YSZURaOTDpj2i+CkJNnIsYkcHuozquSY5PebTLKsRnpn2ZkZ72Uh4nYksDtAU3XqfOq5MpwQtIam2OjsjXAMY+sBCFiRwK3B1xeDU1HerhJ7KxsKwqwo1Z6uSJ2pEqhB+qOl4TlSuAmrXSzgREOMzvqvVxUYMcQYiw+3CkaZdJy0RXp4fZA7fGvoTKkkNy+kGOj2a9xsMkf76aIPkICtwdqvSp2k4LNJG9fMhuVZcFqVNhRJ8MKIjYkMXqgzqPKcEIKMBkUzsy2stvlxatq8W6O6AMkcHug1quSa5Xh71Rwdq6VgA7b5eSZiAEJ3G5qDWi4A7pUKKSIfLuZQekmtta40Xq5jLqsgya6IoHbTXXH51CQE2apY2L/NBp8GnvkyjMRZRK43VQrJWEp54wsCzlWI+9+1trrXq4QoUjgdlOdR8WogMMib12qMCgKF+XbOeZR5XJfEVWSGt1U61HJsRpDFsqL5FPktDAgzcg7Fa14pGJBRIkEbjfVegNywiwFKYrC9MEZNPs1Nh1tiXdzRIqSwO0Gv6ZT79XIs0lJWCoqSDfzxf5p/KfWy3/rPL16LqlWEJ2RwO2G9pml8tKkh5uqLsq3MzjDxNpDzZQ3StWCiCwJ3E6crndS426rUJAebuoyGhSuGe4gx2rk5f2N/KfWg96LygWpzRUnksDthmp3ALMBnLJwZEpLMxn45hlZDEo389dDzfy5vIlmv5xIE72XEF21KVOmYLFYsFrblj655557mDx5Mh999BFlZWV4vV4KCwt59NFHyc3NjVs7a9wq/WwmqVDoA2wmA7NHOdha42FzRQvP76rnkkHpnJltRZHjL3ooIQIXYNmyZYwePTr4s6ZpzJ07l0WLFlFaWsrTTz/N0qVLWbRoUdzaWOMJMMphidvri9gyKApf7J/GSIeZNw4289rBZj5t8HH5kEwsRgld0X0J+914x44dWK1WSktLAZg9ezbr1q2LW3ta/BqtAZ28tIT5GyViJNdm4vrRWXwp385ul48/7nHR0oMhBhnPFQkTuPfccw9XXnklDzzwAI2NjVRWVlJQUBDcnpOTg6ZpuFyuuLSvxi0VCn2ZQVE4b6Cdr4904PKp/OnThh6FrujbEqK7tnLlSvLz8/H5fDz88MMsXLiQSy65JCLPnZub0ent3e1lVB+ftKa/VCj0acMdFr4+Mov/29vAK+WNzB6VhcnQ/eGFSPdy+1qvOVn3NyHSIz8/HwCLxcKcOXO4/fbbufHGG6moqAjep66uDoPBgNPp7NZz19Y2o2kdy3ry8jJPWX+qqwNY4w6QblKwmxPmS4GIk8EZZmYOzWT1gSY2HmlhxpDO/6iHEsn1zzr7fU5lib6/obIk7unR2tpKU1Pbm6frOmvXrqW4uJixY8fi8XjYunUrAKtWrWLGjBlxa2eNW5XxWxE0JtvKpP5pfFTrYY9LJrwR4Yl7gtTW1nLHHXegqiqapjFy5EgWLFiAwWBgyZIlLFiwoENZWDxous4xT4CSfra4vL5ITJPz7Rxo8rHucDNDMsyyxp3oUtwDd/DgwaxevbrTbePHj+e1116LbYM64fJqBHSkhys6MBkULhuSye93u3inspXpg7s/tCD6FvmTHIbq4xUK/SVwxUkG2k1MyLOx7ZiHihZZbl2EJoEbhmp3AAVZ5UF0bnK+nQyTgfWHm2XFCBGSBG4YPnMH6GczYu5B+Y9IfVajgamD0qlyq2w71rtpHUVqk8Dtgq7rfNYaYKBdhhPE6Y1xWhiaYWZzZatcECFOSwK3C03HL+mVwBWhKIrCJYPT8as6b1d0vWLEiZf5yuW+fYcEbhcqW9tOmOVL4Iou9LOZmNg/jY/rvByVE2iiExK4XahqDWBASsJEeM4fmEaG2cAGOYEmOiEp0oXK1gD90uSEmQiP1WhgamE6aw408dExD+Pz0sJ+bJeXlyfw5awiPNLDDUHXdSpaAzKcILrlxBNorXICTZxAAjeEWq+KV9UpTDfHuykiiSiKwiWD0vGpOm8eaY53c0QCkcAN4WhL2wmzwnTp4Yru6Zdm4oJ8O7tcPnbVy+Q2oo0kSQhHm/3YjAo5VrnCTHTfeQPS2NfgY92hZgbaTWT38vdISseSn/RwQzjaGqAw3SSLBooeMSgKVw/PRFHgL/sb8alStdDXSeCehjugUetRZfxW9EqWxchVwzI55lFZfaBRSsX6OAnc0zjU3Fa4PjhDAlf0zgiHhemDM9jf6OfVA00Sun1Ynx7DDTUmdqjZj9kABVISJiJgXD8bXlVjU0Ur/v2NXDUsE6tR+jt9jaTJaRxs8jMo3YxRLngQETJpgB2zQeHNIy38cU8D145w4IzQCdnTdR7kYonEIn9iO9Hi1zjmURmaKcMJIrLG56Uxa6SDZr/GC7tdfCLrofUpEridOHh8/HaIjN+KKBjmsPCtIidZFgOry5tYU94YtSvSZCayxCKB24n9jT5sRkWmZBRRk201cmORk4vy7exp8PHrXfX8q9pNQJMTaqlMEuUkmq6zr9HHSIcFg9TfiigyKgrnD7QzOsvC34628NbRFj6scXPBQDtnZVvl/EEKksA9SUVLAHdAZ1SWJd5NEX1EvzQTs0ZlUd7oY1NFC2sPNfP3ylZK82yck2uLyPLrPRlWkBNukSeBe5J9jT4UYLicMBMxNtxhYVimmfImP+9XudlU0cq7n7UyNsfGhDwb/WzycU12cgRPoOs6u10+BmeYI9KrEKK7FEVhhMPCCIeFz1oDfFjjZnuth23HPAzLNDMhzxaz4S6ZnzfyJHBPUOVWqfOqfLF/+JNGCxEtA+0mrhiaycUF6Xx0PHT/vL8Jp8XA+Lw0zs6xSscgyfTZwO3sr/d/6zwYFChyyvitSBx2s4HzB9qZNCCNT10+tta4eetoC+9WtjKxv42J/dOwyVVrSaHPBm5nPnG1VSekSa9BJCCjojAm28qYbCuftQb4Z1Ur733m5sMaD5P6pzEhLw2LMXaVDTLk0H0SuCcYm2NljNMa72YI0aWBdhNfHe7gs9YAf69s4Z3KVrbUuDlvgJ2SfjZMUlKWkBRdT+2pi2prm9FOKiaXK29EqjnS7Ofvla0cbPaTaTZw/sA0zs61YUywWvJI9Hrz8jITuvccKl8kcIVIIQeafPy9spWjLQEyTAbGZFsozrZSYO9bE+nHM5AlcCVwRR+i6zr7G/18VOthf6MPVQerUaEw3USB3Uw/m5Ecm5FsqxGzDD30SmfBHipfZAxXiBSjKAojsyyMzLLgUTX2Nvg43OznaEuA/Y2tHe6baTaQZTGQZTGSZTXgtBiDPzssBrm8PcIkcIVIYTajgbE5Nsbm2ADwqTr13rZ68zqvSp1HpdGvcrjZz3/rO85YptA2yU5BuokCu4nCdDN5acY+EcKaptPY6CYz04YxgiV3ErhC9CEWo8IAu4kBncyEp2o6TX4Nl0+lwafR4FWp8ajsb/Sxo65t3l6LQWFwhokhGWaGZJgZYDclTQCratsflHAC1OPx4feruFytOJ32iIWuBK4QAgCjQcFpNZ6yCoWu6zT4NI62+DnSEuBQk599x4cmLIa2seH+aSby0ozkWI1kmg2kmxNvOMLlamtzOAHafmZL0/SIhm7CB255eTnz5s3D5XLhdDpZvHgxw4YNi3ezhOgzFOXzID4rp+22Zr/G4WY/h5r9HG3xs7XGzYmrwCtAxvHgTTcpx/97/F/7/5sV0k0GrEal2xUUmqZjOH7C78T/7+ox0Ba82dnpYT2m/XGRCt2Er1K48cYb+drXvsbVV1/NmjVr+POf/8wf/vCHsB8vVQpCRJ+q69R7VFw+jSa/SpNPo9Gv0eLXaAlotPp1WgIanYWNUQG7yYDdpBz/7+f/bzYqmA3t/8BsUNB18HkD9HdYMSng8/rJtFsxGQjZqz6xoiA93YrdfvpL+FtavLS2+jrcZjAop4Rud6sUEjpwa2trmT59Oh988AFGoxFVVZk0aRIbNmwgJycnrOeor285JXBzczOi0VwhRAi6ruNRdVoDGu6ARmtAZ2+jF78KXk3Hr+moGhgUcAc0Aj1MJpPSdhm00aBgVNr+GKgaoOvBwLcaFUzHg7M9otuzWgF8AQ2vqpNjMTA5+/OBAINBITs7PfhzbW3zKa8fKl8SekihsrKSAQMGYDS2jSkZjUb69+9PZWVl2IF74psjhIgfRVFIMykd5ioZneSX0ne38yaztAghRIwkdODm5+dTVVWFqqoAqKpKdXU1+fn5cW6ZEEJ0X0IHbm5uLsXFxbz++usAvP766xQXF4c9nCCEEIkkoU+aAezbt4958+bR2NiIw+Fg8eLFjBgxIt7NEkKIbkv4wBVCiFSR0EMKQgiRSiRwhRAiRiRwhRAiRiRwhRAiRlIqcMvLy5k1axbTp09n1qxZHDhw4JT7qKrKz3/+c6ZNm8Yll1zCyy+/HNa2RNTb/X3yySc577zzuPrqq7n66qv5+c9/HsPWd184+/vuu+9yzTXXMHbsWBYvXtxhWyoe31D7m0zHN5x9Xb58OVdccQVXXnkl11xzDX//+9+D29xuNz/+8Y+55JJLmDFjBps2bYph67tBTyE33HCDvnr1al3XdX316tX6DTfccMp9XnnlFf3mm2/WVVXVa2tr9cmTJ+uHDx/uclsi6u3+Llu2TH/kkUdi2ubeCGd/Dxw4oO/cuVN/7LHHTtm3VDy+ofY3mY5vOPu6efNmvbW1Vdd1Xd+1a5c+YcIE3e1267qu608++aQ+f/58Xdd1vby8XD///PP15ubmGLU+fCnTw62trWXnzp3MnDkTgJkzZ7Jz507q6uo63G/t2rVcd911GAwGcnJymDZtGuvWretyW6KJxP4mk3D3d+jQoRQXF2MynTpNSDK9F5HY32QR7r5OnjyZtLQ0AIqKitB1HZfLBcBf//pXZs2aBcCwYcMYO3Ysmzdvjt1OhCllAjfURDcn36+goCD4c35+Pp999lmX2xJNJPYX4I033uDKK6/k5ptvZtu2bbFpfA+Eu79dPUeqHd+uJMPx7cm+rl69miFDhjBw4EAAKioqKCwsDG5P1GObvH8WRa/Nnj2b733ve5jNZt577z2+//3vs3btWrKzs+PdNBEBqXp8//Wvf/HEE0/w29/+Nt5N6baU6eGGO9FNfn4+FRUVwZ8rKyuDfyVDbUs0kdjfvLw8zGYzABdccAH5+fl8+umnMdqD7onEREapeHxDSZbj25193bZtG3PnzmX58uUdLvEvKCjg6NGjwZ8T9dimTOCGO9HNjBkzePnll9E0jbq6OjZu3Mj06dO73JZoIrG/VVVVwfvt2rWLo0ePMnz48NjtRDdEYiKjVDy+oSTL8Q13X7dv385dd93FsmXLOOusszpsmzFjBv/3f/8HwIEDB/j444+ZPHlybHagO+J91i6S9u7dq1977bX6pZdeql977bX6vn37dF3X9VtuuUXfvn27ruu6HggE9LKyMn3q1Kn61KlT9VWrVgUfH2pbIurt/t577736FVdcoV955ZX6Nddco7/99ttx2Y9whbO/W7Zs0SdPnqyXlJTo48aN0ydPnqxv3rxZ1/XUPL6h9jeZjm84+3rNNdfokyZN0q+66qrgv08++UTXdV1vaWnR77jjDn3atGn6pZdeqr/55ptx25dQZPIaIYSIkZQZUhBCiEQngSuEEDEigSuEEDEigSuEEDEigSuEEDEigSti7sMPP+TSSy+lpKSEjRs3xrUtJSUlHD58OC6vXVZWxvLlyyPyXBUVFZSUlAQvHrjhhhsiOhvaLbfcwiuvvBKx5+uz4l2XJiLj4osv1r/whS/o48aN00tLS/Vbb71Vr6ioCOux77//vj558uQev/bhw4f10aNH636/P6z733jjjfoLL7zQ49drN3r0aP3AgQNh3//666/XX3rppV6/bjhOPB4TJkzQZ82apf/pT3/SVVXt0XO999573XpMb/Z12bJl+k9+8pMePVaEJj3cFLJixQq2bdvGu+++S25uLg8++GC8m9SpiooKzjjjjHg3I+raj8emTZu49dZbee6555g/f37EXycQCET8OUV0SOCmIKvVyowZM9i3b1/wNp/Px+LFi/nyl7/M+eefT1lZGR6Ph9bWVm699Vaqq6spKSmhpKSEqqoqtm/fzqxZsygtLeXCCy9k4cKF+Hy+sF5/3rx5/PznP+e73/0uJSUlXHfddRw6dAiAadOmcfjwYb73ve9RUlKCz+ejqamJ++67jwsvvJDJkyfz+OOPB78aHzx4kOuvv54JEyYwadIkfvzjHwPwzW9+E4Crr76akpIS1q5dS0NDA7fddhvnnnsuEydO5LbbbgvOGPX444+zdetWFi5cSElJCQsXLgTapvk7ePAgAE1NTdx7772ce+65XHzxxTz99NNomgbAX/7yF77xjW+wePFiJk6cyJQpU3jnnXfCej8yMzOZOnUqv/rVr3jllVfYs2dP8H16/PHHAairq+O2226jtLSUL37xi8yZMwdN05g7dy4VFRXB9+u5557jyJEjFBUV8fLLL/PlL3+Zb33rW8HbTgzfQ4cOce211zJ+/Hhuv/324FSGH3zwARdddFGHNk6ZMoV//OMfbN68mWeffZa//vWvlJSUcNVVVwEdhyg0TePpp5/m4osv5rzzzuPee++lqakJINiOV155hS9/+ctMmjSJZ555Jqz3qS+QwE1BbrebtWvXcs455wRvW7p0KeXl5axevZoNGzZQXV3N8uXLsdvtPPfcc/Tv359t27axbds2BgwYgMFg4H/+5394//33WbVqFf/85z/505/+FHYb1q5dyw9/+EO2bNnCkCFDgsGyceNGCgoKgr0/i8XCvHnzMJlMbNiwgdWrV/Pee+8FP9xPPPEEF1xwAVu2bGHz5s1cf/31AKxcuRKANWvWsG3bNi6//HI0TeOaa65h06ZNbNq0CavVGgzWu+66i9LSUsrKyti2bRtlZWWntPnBBx+kqamJjRs38sc//pE1a9bw5z//Obh9+/btDB8+nPfff59bbrmF+fPno3fjQs2zzz6bgQMHsnXr1lO2/e53v2PAgAH885//5L333uPuu+9GURQeffTRDu/XrbfeGnzMli1bWLt2Lb/5zW86fb3Vq1fzi1/8gnfffReTycRDDz3UZRsvuugibrvtNi677DK2bdvGq6++esp9/vKXv/DKK6/whz/8gY0bN9La2hp8n9t9+OGHrFu3jt///vcsX768wx//vkwCN4X84Ac/oLS0lNLSUt577z2+853vAKDrOi+99BL33XcfTqeTjIwMbrvtNt54443TPtfYsWMZN24cJpOJQYMGMWvWLLZs2RJ2W6ZNm8bZZ5+NyWTiqquuYteuXZ3e79ixY7zzzjvcd9992O12cnNz+fa3vx1sm8lkoqKigurqaqxWK6Wlpad9zezsbKZPn05aWhoZGRncfvvtYbdZVVXWrl3LT37yEzIyMhg0aBA33XRTh8ApKCjg61//Okajka9+9avU1NRw7NixsN8TgP79+9PQ0HDK7SaTiZqaGioqKjCbzZSWlqIoSsjnuuOOO7Db7dhstk63X3311YwePRq73c6dd97JunXrgt8ceuO1117j29/+NoMHDyY9PZ27776btWvXduhd//CHP8RmszFmzBjGjBnDJ5980uvXTQUyH24KWb58Oeeffz6qqvK3v/2NG264gTfeeAODwYDb7eaaa64J3lfX9eDX5c6Ul5fzyCOPsGPHDtxuN6qqnjJDUyj9+vUL/r/NZqO1tbXT+1VUVBAIBLjwwguDt2maFpyab+7cuTzxxBNce+21ZGVlcdNNN3Httdd2+lxut5tFixbx97//PRhqLS0tqKoanNz6dOrr6/H7/R0mKC8oKOgw49aJ+9S+8sDp9ut0qqqqyMrKOuX273znOzz11FPcfPPNAMyaNYvvfve7IZ+rq+kHT5zesKCgAL/fT319fbfa25nq6uoOk30XFhYSCASora0N3nbye9Xd9ylVSeCmIKPRyKWXXkpZWVmwBMtms/HGG28wYMCAU+7fWU/qgQce4Mwzz+SXv/wlGRkZvPDCC6xfvz7ibR04cCAWi4X333+/02Vi8vLygl+Ft27dyk033cTEiRMZOnToKff97W9/S3l5OS+99BJ5eXns2rWLr3zlK2F97c/OzsZsNlNRUcGoUaOAz1ciiJTt27dTVVXFhAkTTtmWkZHBvHnzmDdvHnv27OFb3/oWX/jCFzjvvPNO+3xd9YBPXDGhsrISs9lMdnY2aWlpeDye4DZVVTssZ9PV8/bv37/D3LMVFRWYTCZyc3MTcpWFRCJDCilI13U2btxIY2MjI0eOxGAwcN111/GLX/wi2AupqqoKrnqam5uLy+UKnviAtp5heno66enp7Nu3jxdffDEqbe3fvz8XXHABjzzyCM3NzWiaxqFDh/jXv/4FtK1V1f4hzsrKQlEUDIa2X9t+/fp1qKFtaWnBarXicDhwuVw89dRTHV7r5PufyGg0MmPGDB5//HGam5s5evQov/vd74InjXqjubmZTZs2cffdd3PVVVdRVFR0yn02bdrEwYMH0XWdzMxMjEZjMPhCtTuUV199lb179+J2u3niiSeYPn06RqOR4cOH4/V6efvtt/H7/TzzzDMdTojm5uZy9OjR034DmjlzJr///e85fPgwLS0tPP7441x22WVJva5arEjgppD2M9njx4/nV7/6FY888kiw/Gru3LkMHTqUr3/964wfP55vf/vblJeXAzBy5EiuuOIKpk2bRmlpKVVVVfz0pz/l9ddfZ/z48dx///1cfvnlUWv3kiVL8Pv9XH755UycOJEf/ehH1NTUAPDxxx9z3XXXUVJSwu233878+fMZPHgw0DZOOG/ePEpLS1m7di3f+ta38Hq9nHvuucyaNeuUCahvvPFG1q9fz8SJEzs9gXT//feTlpbGtGnTmDNnDjNnzuRrX/taj/er/Xh86UtfYsWKFdx0000sWrSo0/sePHiQm266iZKSEmbNmsU3vvENzj33XAC++93v8swzz1BaWnraE2Sdufrqq5k3bx4XXHABPp8vWJKWmZnJggUL+NnPfsZFF11EWlpah+GJGTNmADBp0iS++tWvnvK8X/va17jqqqu4/vrrmTp1KhaLhfvvvz/sdvVlMh+uEELEiPRwhRAiRiRwhRAiRiRwhRAiRiRwhRAiRiRwhRAiRiRwhRAiRiRwhRAiRiRwhRAiRiRwhRAiRv4/BUZK2M3e61gAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "ax = sns.displot(infest_beta,\n", " kde=True,\n", " bins=100,\n", " color='skyblue',\n", " linewidth=15,\n", " alpha=1)\n", "ax.set(xlabel='Beta Infestation Distribution', ylabel='Frequency')" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Infestation rate: 0.0301\n" ] } ], "source": [ "# Choose infestation rate for shipment from beta distribution\n", "infest_rate = beta.rvs(4, 60, size=1)\n", "print(\"Infestation rate:\", round(infest_rate[0],4))" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Num stems: 152656\n", "Num pests: 4590\n" ] }, { "data": { "text/plain": [ "array([ 0, 1, 2, ..., 152653, 152654, 152655])" ] }, "execution_count": 6, "metadata": {}, "output_type": "execute_result" } ], "source": [ "num_stem = np.random.randint(1, high=1000000, size=1) #QUANTITY\n", "\n", "# Calculate number of pests in shipment based on infestation rate.\n", "num_pest = int(num_stem * infest_rate)\n", "print(\"Num stems: \", int(num_stem))\n", "print(\"Num pests:\", num_pest)\n", "np.arange(num_stem)" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [], "source": [ "# Infest stems using uniform distribution\n", "infest_stem_unif = np.random.choice(int(num_stem), num_pest, replace = False)" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 8, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAUwAAAFtCAYAAABsnlPMAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/Il7ecAAAACXBIWXMAAAsTAAALEwEAmpwYAAA2e0lEQVR4nO3deXQUZb4+8Kf37ux7SAISQcAAYgIBZFMJaBACqIwDwxXm6KjjvS54XUau3gHENcCIiihzjo7z8w4jVwcURC6gMuKCIgybgGwxQCAhSychSe9d9f7+CLREAqlOuqu7k+dzjufYXd3V366uPLxV71tvaYQQAkRE1CZtqAsgIooUDEwiIoUYmERECjEwiYgUYmASESnEwCQiUkgf6gICwWptgiwrGx2VmBiFujp7kCvqmHCvMdzrA1hjIIR7fUBgakxNjVX82i7XwtTrdaEuoU3hXmO41wewxkAI9/oA9WvscoFJRNReDEwiIoUYmERECjEwiYgUYmASESnEwCQiUqhTjMMkouDzZ7zipVRXNwagktBhC5OISCEGJhGRQgxMIiKFGJhERAoxMImIFGJgEhEpxMAkIlKIgUlEpBADk4hIIQYmEZFCDEwiIoUYmERECjEwiYgUYmASESnEwCQiUoiBSUSkEAOTiEghBiYRkUIMTCIihRiYREQKMTCJiBRS7a6RLpcLL7zwAr799luYTCbk5ubi2WefRWlpKebOnYv6+nokJCSguLgY2dnZapVFRKSYaoG5ePFimEwmbNq0CRqNBjU1NQCA+fPnY+bMmZg6dSrWrl2LefPm4d1331WrLIpgbd32VcltYSP9tq+kLlUOyW02Gz766CPMmTMHGo0GAJCSkgKr1YqDBw+iqKgIAFBUVISDBw+itrZWjbKIiPyiSguzrKwMCQkJeP3117F9+3ZER0djzpw5MJvNSE9Ph06nAwDodDqkpaWhoqICSUlJitefnBzjVz2BuCF9sIV7jeFen1Kh/h6h/vy2BLq+YHxfNbehKoEpSRLKysrQv39/PPnkk9i7dy/uv/9+vPrqqwFZv9XaBFkWil6bmhob9odh4V5juNQXiD+UUH6PcNmOl/LL+sJxewdiG/rzvVQ5JM/IyIBer/cdel977bVITEyE2WxGZWUlJEkC0BysVVVVyMjIUKMsIiK/qNLCTEpKwvDhw/HNN99g9OjRKC0thdVqRXZ2NnJycrB+/XpMnToV69evR05Ojl+H46RcOLYQqGvpyD4YDvuear3kzzzzDJ566ikUFxdDr9dj0aJFiIuLw4IFCzB37ly88cYbiIuLQ3FxsVolERH5RbXA7NGjB/7nf/7noud79+6NDz74QK0yiIjajVf6EBEpxMAkIlKIgUlEpBADk4hIIQYmEZFCDEwiIoVUG1ZEwRPu1yOHs9a2XTgMkPYXL0pQB1uYREQKMTCJiBRiYBIRKcTAJCJSiIFJRBEhNTU25B2cDEwiIoUYmERECjEwiYgUYmASESnEK33CzPmT2qE+ua22zvJ9/f0eSl/Pq3DCA1uYREQKMTCJiBRiYBIRKcRzmApwJpifXXiOtbN8p0jQWc7xBsIvt4WSbROofZUtTCIihRiYREQKMTCJiBRiYBIRKcROH6Jf6Cy3rfDXhd+7K3zf9mALk4hIIQYmEZFCDEwiIoUYmERECrHTRyU8oU6RpKvOmtUWtjCJiBRiYBIRKcTAJCJSiOcwQ6CznBfqSgO8O8tv1lWd//06un+yhUlEpBADk4hIIQYmEZFCqp3DLCgogNFohMlkAgA8/vjjGDNmDPbs2YN58+bB5XIhKysLixcvRnJyslplEREppmqnz2uvvYa+ffv6HsuyjCeeeAIvvvgi8vPz8cYbb2DJkiV48cUX1SyrBZ7c7xg1bxdApLaQHpLv378fJpMJ+fn5AIAZM2Zg48aNoSyJiOiSVG1hPv744xBCYMiQIXj00UdRUVGBzMxM3/KkpCTIsoz6+nokJCQoXm9ycoxfdbAVGVrc/hQqHd33VAvMlStXIiMjA263G88//zwWLlyIm266KSDrtlqbIMtC0Wvbutsh/5iDr7Xtz+1OaujovqfaIXlGRgYAwGg0YubMmdi1axcyMjJQXl7ue01tbS20Wq1frUsiIrWoEph2ux2Njc3JLoTAhg0bkJOTg4EDB8LpdGLnzp0AgFWrVmHChAlqlERE5DdVDsmtViseeughSJIEWZbRu3dvzJ8/H1qtFosWLcL8+fNbDCsiIgpHGiGEspN/YYznMCMLz2FSqETMOUwiokjH2Yo6Oa8scKLRg9JGN6xOCXavDI1Ggyi9BilmPXrGGNAz1gC9VqNaTWxNUqRiYHZS1Q4vvqt04MhZFzwyYNACySY9YgxaCABNHhknGh34vsoBs06Da5PNuC7dAoueBx1El8LA7GRqnRL+WW7D0bNuGLTAgEQz+iYYcUXMxa1IjyxwstGDH2qd2F7lwB6rE9elWTAszQKdii1OokjBwOwkPLLAt5V2bK90QK/RYFQ3C/JTL99iNGg16B1vRO94I6ocXnxZYcfWCjt+rHdhSs9YpFi4exBdiH8RncAZuxfrjjei1iVhQKIJY7OiEWPw79A6zaLHr3rF4ehZF/7vZBP+erge47pHIy/FEqSqiSIPAzOCyULg+yoHviy3I9qgxYzecciOM3ZonX3iTci82oBPTjRiU5kN1Q4J47tHQ6vhIToRAzNCNbglrD/RhJNNHvRLMGJCj5iAddhEG7T4Ve84bC23Y3uVA7UuCbdmx8IcZh1CDW4JP9a5cLzRg2qnBKdXhl6rQbxRi+4xBuQkmJAVrYeGYa+Ks24JJxs9sLokuCWBKL0WWdF6XBFj6DTnxBmYEehQnQsby5ogCYGJV8TgmiRTwENBq9FgbFY0ks06bCxrwsqjZ/Hr3nGINeoC+jntYXV68XWFHYfq3RAAUs06ZMcaEK3XwiMLWJ0S9tY48a9qJzKj9LghMwo9YzvW8lZbg1tChd2Ls24ZHlnAeO4fgoxoPWINof8NzpOFwJGzbuyocuC0zQugeXC3QaeBS2q+mCTWoMXojCgMCsJ+qjZe6dPK8nDlkmR8dsqGH2pdyIjSY0p2LBJNwf/jOd7gxprSRph1GtzROw6pIeoMEudOQWwtt0OnBfJSLMhLMbe6DVySjIN1Lnxb6UCDW8agJBNu6hEDQxi3dNySwD6rE3utTlQ7pUu+Lt2iw1XxRlyTZEaCCr9/a4QQ+LHOja/P2FHrkpBo0uLaZDN6xxmRbNZBq9HALQkcb3Tju0oHyu1eXJ1gxKSesSH9DTp6pQ8Ds5Xl4ei0zYOPjzfirFvGiHQLRmVEQafiv9aVdi/eLzkLrwB+1SsOPWIMqn020Dxu9JMTjShtbD4FcXP3GEQr6NjyyALbztjxbaUDqWYdbu8Vp8o/Mv4QQmCP1YmvKuywewUyovTISTShR4weiUYdDFoN3LJArUtCWZMHx866fa25PvFGDE2zqPp7nLF78dmpJpyyeZFm0WFkehT6JhgveZ5bCIHtVQ58UW7HlbEGTOsVp+qFEhdiYOLygRmuAaiUVxb45owd31U6EGvUYnLPWNXD6rx6l4T3Sxpw1i1hcnYsrk4wqfK5PzW4sf5EI9ySwLju0chNNvt9aPdTgxvrjjdCALj9ytiwOURv9EhYd7wRZU1e9IjR48bMaGRFt/37Nrgl7KpxYk+NE05J4IoYA0Z1s+CKGEPQDnvtHhlfVtixx+qERa/BjRnRuCbZpLhDcK/Vif872YRBSSZM7Bmav0sGJjpvYJ5odGNTmQ21LgnXJJkwrns0zLrQdrw4vDL+8VMDTtu8GN89GvmpwRt2JMkCWyvs+L7KgRSzDlOzYzt0OqDeJeEfPzWgziXh1itj0SdencC/lNIGNz4+0QiPLHBT9/adi/bIAntqnPiu0g6bV6B7tB6ju0WhZ2zgglOSBf5V48Q3Z+xwSwJDUs0Y3S2qXZ2AX5bbsK3SgaKeMRiYZA5Iff5gYKLzBWatU8IX5TYcOetGglGLCT1iOjxcKJA8ssC64404etaN69IsuCEzKuCtmjpXc8urwu5FXooZBVnRATn35fDKeL+kAZV2L4p6xqJ/UmhCc1e1A5tP2ZBi1uHWK2ORYu7YeWGP3Hz+87tKBxo9MrKi9RjVLQpXdiA4hRAoafBgy+nmf7SvjDVgXFZ0hy5okIXA34+eRY1Twn05iYjyc7xwRzEw0XkCs8bhxffVDuy3uqDTAtelR2FYmiUsOypkIfDpKRt21zgxINGEiVfEBGzoyIFaJzaV2aDRABOviEG/AB/6u6TmVvKpJi+KesZggIotHSEEtlU68FWFHVfFGTElOxZGXeB+X68s8EOtE9+ecaDBIyMjqjk4e8cpD05ZCBypd+PbSjsqHRKSTDoUZEX7tY7LqXF48ZdD9RgYgkNzBiYiOzCFEDjR5MGOKgdKGjzQa4BByWaM7Bbl99U6ahNC4NtKB76ssKNnjAGTs2M7VLPTK+Oz0zbsr3Whe7Qek7NjER+kYUxuSeAfPzWgrMmDySq1NIUQ+Py0DTurnc1hcUVM0C4IkGSB/bUubKu046xbRqxBi74JRlwVZ0SaRX9Rh5nNI6PS4cXRs24crnfB7hVIMukwPN2CgYmmgI+j3HLahu+rHPjd1QmqjrpgYCIyA1MSAofqXPi+yoFKh4QovQaDUywYnGJW/TClo36wOrGprAkGrQYT2tkiLDnrxv+VNcHmkTGymwWjukUF/eoityTwwU9ncarJiynZschJDF5oSkJgw4kmHKhzIT/VjHFZ0aqMSTy/nx2qd+OnBjfODY2ERa+B6VwIumQBh7d5gUEL9I4zon+iCVfFX7rnu6PsXhkrDtShV5wBt14ZF5TPaA0DE5EVmB5ZYHeNEzurmg+Zkk06DEuzYECSKWRDLQLB6vTi4xNNOGP3om+8ETdmRiPJ3Hbr8Izdi68qbChp8CDVrMPEnjHIiFJvFIAaoemRBT4qbUBJgwfXZ0RhRLolJAO43ZLAaZsHNU4JNU4vvHLz83otkGLWI9WsQ1aMQbVTQFvLbfi20oF7cxKQ3MFzuEoxMBEZgek915v57bnezB4xegxP8+/cUriThMD2Sge+rbTDIzePERyQaELPWEOLyzZtHhknGj3YX+vET40emHUaDE+zYGiaJST/aLglgfdLzuK0zYup2bG4OoCh6ZRk/KOkAadsXhT24GQmF7J5ZLxxoBbXJptxc48YVT6TgYnwD8wKuwefnGhCjVPCFTEGjMmICtlYSjXYPDJ2Vjuwz+qE7dyhnkWngfHc5XLOCy6Zy00xIz/VDFOIh0u5JBkflDQPl5p6ZWDGmDa4JXxQ0gCrS8LknsE95I9UH58bbfHAwERV9gEGJsI3MIUQ+OaMA9+caZ5NaEKPGFwVHz7Dg4JNFgKnbV6ctnlQ75LhFQIGrQaJJh26R+vRLUofVrMguaTmIUflNi9u7mBrsMrhxQclDXBLArf1ikV2mAyUDzcVNg/+35GzuKl7NIYEcUzveR0NTE6+ESReWeDjE404XO/GgEQTbuoeHXaz/QSbVqNBjxhDxLSmTTotft07DuuON09tV+OUMC7L/6nt9tc6sbnMBpNOg3/rG480TsR8SRnRBqSaddhf61IlMDuKv2QQSLLAh+dO8o/NbB5L2VnOU3Z2Jp0W03rF4Z+nbdhR7fQNcFcyyYVLkvHpqeZhUT1i9JjSMzYsZncKdwOTTPhnuR1Wp1e1zp/26lpNHhUIIfB/ZU0oafCgsEc0hqcH/ioYCi6tRoNx3WNQ1DMGVQ4Jb/1Yh63lNtg8cquvd0sC/6p24M8H63Cg1oXR3aLwm6viGZYKDUgyQwPgQK0r1KW0KbzjPALtqHZi/7k/GvaIRraBSWZcEWPAF+XNsx1tr3KgZ4wBGVF6mPVauCWBKocXxxs9cMvN13GP6x2t6rCoziDGoEXPWAN+rHdhTEZ4NzAYmAF0xu7FF6dt6BtvxKhuDMvOIM6ow5TsWIzqZsFeqws/NbhxvNKB812M8UYt+ieaMDDJhO4Rcq42HPVLMPrOG4dqvlUlwreyCCOd6+SJNmhxyxUxYf2vJPkv2axHQZYeBVnRkISAW2ru8Y/kiw3CSZ94EzaV2XC43h3WgclzmAHyfZUDVqcU0HvrUHjSaTSw6LUMywCKMWjRPVqPw/XhfR6Tf9kB0OiWsK3Sjj7n7vFNRP7rE29EtVNCg/vSt+cItU4ZmKmpsb7/1LCt0gFJAOOyolX5PKLOqNe5OV9LGzwhruTSOmVgqqne1XyHwtzk0N2QiqgzSDHrEGvQoqTBHepSLomB2UHbqxzQaIAR7BUn6hCNRoNecQacaPRACtMrthmYHWD3yPjB2jwZbDjdK5ooUvWKM8IlC5SfuytmuGFgdsCuGie8AhiaxtYlUSBccW4s68mm8DyPycBsJ/ncvaR7xRo6fAMrImpm0WuRZtHhRCMDs1MpaXCjySMjN0X9W4USdWY9Yww4bfPAe4kpG0OJgdlOe2qciNFrOe6SKMB6xhohCeC0LfxamQzMdrB5ZPzU4ME1ySboeAkkUUB1j2k+xVXWFH4dP6oH5uuvv45+/frhyJEjAIA9e/ZgypQpKCwsxN133w2r1ap2SX47VO+CADCAtxwgCjizTotUsw7lXb2FeeDAAezZswdZWVkAAFmW8cQTT2DevHnYtGkT8vPzsWTJEjVLapcDtS6kWXRICeNJAogiWWa0HuV2L8LtDjqqBabb7cbChQuxYMEC33P79++HyWRCfn4+AGDGjBnYuHGjWiW1S71LQrndi/5sXRIFTVa0AU5JoNYVXteVq9ZEevXVVzFlyhR0797d91xFRQUyMzN9j5OSkiDLMurr65GQkKB43cnJ6tyiEwCOnG2+bCsQdxUkotZlRjdH02lbYG9b0dH5JVQJzN27d2P//v14/PHHg7L+X941MpiTbhypdyHVrON140RBlGzSwaTT4LTNg0HJgRu619G7Rio+JP/ss8/g9bav12rHjh0oKSnBuHHjUFBQgDNnzuB3v/sdTpw4gfLyct/ramtrodVq/WpdqsnmkXHK5kXfBA4lIgomjUaDrCh92F0iqTgwX3vtNYwePRoLFy7E3r17/fqQ++67D19//TW2bNmCLVu2oFu3bnj77bdxzz33wOl0YufOnQCAVatWYcKECf59AxUdOzeLSt94Ho4TBVtmtAHVTglOqfWbz4WC4kPydevW4dChQ1i7di0eeughWCwWTJ069aLzkv7QarVYtGgR5s+fD5fLhaysLCxevLhd61JDyVk3Yg3Nl24RUXBlnTuPecbmRXZceBzVaUQ7+u2FEPj222/x0ksv4ejRoxg8eDCmT5+OoqIiaLXqj4VX4xymJARe21eLqxONuOUKdSYmJurKnJKMV/bVYkxGFEZ1iwrIOjt6DtPvTp+TJ09i3bp1WLduHTQaDR5++GFkZGRg5cqV2Lx5M15//XV/VxkRTtu8cMnCNys0EQWXWadFilkXVpdIKg7MlStXYu3atThx4gRuueUWLFq0CLm5ub7lhYWFGDlyZDBqDAs/NbihBZAdy1upEqklM1qPI/VuCCHC4k6sigPzyy+/xF133YVx48bBaLy4lWWxWLBs2bKAFhdOTjR6kBmth0nHy++J1NLNosc+qwsNHhnxxtD3HSgOzNdeew1arRYGw88tLI/HAyGEL0BHjx4d+ArDgEuSccbu5W0oiFSWHtUcUZV2b1gEpuLm0t13340DBw60eO7AgQP43e9+F/Ciws2pJi8Efp4NmojUkXruKp8qR3hcIqk4MA8fPoxrr722xXODBg3CoUOHAl5UuDnZ5IFW03x9KxGpx6jTINmkQ6UjPAawKw7MuLg41NTUtHiupqYGFkvnP0w92eRBZpQeBm3oTzoTdTXpUXpU2iMsMG+++WY89thjOHLkCBwOBw4fPownn3wSt9xySzDrC7nz5y95OE4UGmkWHRo8Mhze0F/xozgw//M//xO9e/fGHXfc4RuofuWVV+LRRx8NZn0hd9rG85dEoZR+bt7ZcDgs9/tKHyEE6urqkJiYGBbjooDgXunzxWkbvq924JFrkmHUhcf3JepK7F4Zr/1Qi7GZURie3rErflS90qexsRGlpaWw2Wwtnh8xYoQ/q4koJ5s8yIjSMyyJQiRKr0WsQRsWPeWKA3PNmjVYuHAhoqKiYDb/PD+dRqPB559/HpTiQs0tCVTYvbguvfN3bBGFs3SLPiwOyRUH5tKlS/Hqq6/ihhtuCGY9YaXc7oEA0IPnL4lCKi1Kh5IGNzyyCOloFcWdPpIkddoreS6l4tzkpRlRvNkZUSilW/QQAKpD3MpUHJj33nsv3nzzTchy6Lv21VJu9yLRpIVFz+vHiUIpXHrKFTed/vrXv6KmpgZvvfXWRbeQ+OKLLwJcVugJIVBh86InZyciCrl4oxZGrQbVIe74URyY4TwTejA0emQ0eWXf3euIKHQ0Gg1SzDrUOCMkMIcNGxbMOsJOuZ3nL4nCSYpZ57uvVqgoPjnndruxdOlSjBs3DkOGDAEAfP311/jb3/4WtOJCqcLmhU4DpFkYmEThIMWih90rYPOErh9FcWC+8MILOHLkCJYsWeK7wqdPnz547733glZcKJXbPUi36KHnhBtEYSHV3DwfZo0zdB0/iptPn332GTZv3oyoqCjfjc7S09NRWVkZtOJCRRYCZ+zegN5Anog6JsVyPjAl9AzRfQgVtzANBgMkqeUJ19ra2ot6zDuDGqcEjwxk8vwlUdiI0Wth1oW2p1xxYE6YMAFPPvkkysrKAABVVVVYuHAhJk2aFLTiQqX83ID1TE4YTBQ2fu4pD90huV/Tu3Xv3h1TpkxBQ0MDCgsLkZaWhgceeCCY9YVEhd0Ds06DBCMHrBOFk1SLHtVOCX5OshYwfk/vBjQfinfm6d3eOVQHi16LGVfFd7Q0Igqgf1U78OkpGx4YmIhYg/83RVNterfzh+LnXTjFW48ePRR/YLiThECNU0J+6sW3Eiai0Eo531PukNoVmB2lODBvuukmaDSaFk3h8y3MH3/8MfCVhUiNQ4Ikfr69JxGFj/N3kax2SrgyTv3PV5wKv7w7ZHV1NV5//XXk5+cHvKhQqjp3cX+6JfT3QCailqIMWkTpNagJ0SQc7e7VSE1NxdNPP42XX345kPWEXKXDC4MWSDQxMInCUYpZH7JryjvUDfzTTz/B4XAEqpawUOnwItWshzZMOrSIqKXzk3CEoqdc8SH5zJkzW/SKOxwOHDt2rFMNKxJCoMouoX+SKdSlENElJJt1cMsCTV5Z9Y4fxYF5xx13tHhssVhw9dVXIzs7O9A1hcxZtwyXLHyTlRJR+Ek+d7qs1ql+T7niZLjtttuCWUdYOMMOH6Kwl3RuaJE1BNeUKw7MV199VdHr5syZ0+5iQq3K7oUGzVcTEFF4ijVoYdACVpf6HT+Kk+HEiRPYvHkzBg4ciKysLJSXl+OHH37AzTffDJOpc5zzq3R4kWLWcUo3ojCm0WiQbNKjNgQ95YoDUwiBP/3pTygsLPQ9t3nzZmzcuBEvvvhiUIpTW6VDQk/eUpco7CWZdThl86j+uYqHFX355ZcYP358i+cKCgqwdevWgBcVCjaPjCaPzCt8iCJAslmHBrcMj6zu0CLFgdmzZ0+sXLmyxXPvvfcerrjiCkXv/4//+A9MmTIFt956K2bOnOm7nLK0tBTTp09HYWEhpk+fjuPHjyuvPoDO3+84jR0+RGHvwp5yNSmerejgwYN48MEH4fV6fTOt6/V6LFu2DAMGDGjz/Y2NjYiNbe7S+uyzz7B8+XJ8+OGHmD17NqZNm4apU6di7dq1WL16Nd59912/vkQgZivaUeXA56dteGhgEqINnNaNKJxVObz4y6F6TM2ORU6i8j4U1WYr6t+/PzZt2oS9e/eiqqoKqampyM3NhcGg7Jzf+bAEgKamJmg0GlitVhw8eBDvvPMOAKCoqAjPPvssamtrkZSUpPhLBEK1w4sovYZhSRQBzl+6bFW5hdnuE3ZDhw6F3W6Hx+NBVFSUovc8/fTT+OabbyCEwFtvvYWKigqkp6dDp2v+8jqdDmlpaaioqFA/MJ2SbyYUIgpvBq0G8UYtrCrPvq44IQ4fPox///d/h9FoRGVlJSZOnIgdO3bgww8/xCuvvKJoHc8//zwA4KOPPsKiRYsCNmYzOTmmQ+8XQqDG6cW1vOkZUcRINulQ6+dYzI5OLq74+HPBggV4+OGHsXHjRuj1zTk7dOhQ/Otf//L7Q2+99VZs374d3bp1Q2Vlpe/mapIkoaqqChkZGX6tz2ptQnV1o+8/f9W7ZXhksIVJFEGSzM2B6c8kHBfmRHvyQnFgHjt2DFOnTgXw88TBUVFRcLlcbb7XZrOhoqLC93jLli2Ij49HcnIycnJysH79egDA+vXrkZOTE5LzlwCQyh5yooiRbNbBIwMNHlm1z1TcpMrKysL+/ftxzTXX+J7bt2+fomFFDocDc+bMgcPhgFarRXx8PFasWAGNRoMFCxZg7ty5eOONNxAXF4fi4uL2fZMOqD534jiFLUyiiJFsav57rXVKiDeq09hRnBBz5szB73//e8yYMQMejwd//vOfsWrVKjz77LNtvjclJQXvv/9+q8t69+6NDz74QHnFQVDt8CLeqIVRx0siiSKFbxIOl4QrVfpMxYfkY8eOxVtvvYXa2loMHToUp0+fxrJlyzB69Ohg1qeKaqfECTeIIky0XgOTTqPq4HVFKSFJEgoLC7FhwwYsWLAgyCWpyysL1Dol9IvnXSKJIknzJBw6VcdiKmph6nQ66HQ6RR08kcbqlCAApLCFSRRxksw6Vad5U3xIPnv2bDzyyCP4/vvvcfLkSZSVlfn+i2TV5wa+pprZQ04UaZJNOjR5ZLgkdXrK22xWVVdXIzU11de5s23btovuTR7J9yWvdkjQan4+gUxEkeP8322tS0JGVPAva27zE87Pf3no0CEcOnQIBQUFvv8/dOhQRIclANQ4vUg26aDjXSKJIk6yWd1rytsMzF+Oot+xY0fQigmFGqeEFLYuiSJSolEHDdSb5q3NwNT8ouUVinsBB4tbEjjrlpHMAetEEUmn1SDBpFWt46fNpJAkCd99950vKH/5GABGjBgRvAqD6PyF+2xhEkWuJJNOtRZmm4GZnJyMp556yvc4ISGhxWONRoPPP/88ONUFWc25HnIGJlHkSjLpcKLRAyHERUfEgdZmYG7ZsiWoBYSS1SlBi58nIyWiyJNk1sErgEaPjLggX1PepacXr3FKSDLroONtdYki1vkGj79zY7ZHFw9Mr29YAhFFJjVviNZlA9MrC9S7ZJ6/JIpwMQYtDFq2MIOq1nXuGnIOKSKKaBqNBokmHeoYmMFTc675zkNyosiX1I77+7RHFw5MLzRo3tBEFNmSTDrUu2RIcnAvrOmygWl1Skg06aBnDzlRxEsy6yAA1LuD28rssoFZ45R4OE7USSSpNLSoSwamJAvUcdINok4jUaWhRV0yMOtcEmTwkkiizsKi18Ki17CFGQw/95BzSBFRZ6FGT3kXD0y2MIk6iySTDnXO4N6qoksGptXpRYJRCwN7yIk6jSSTDk3e4N7fp0sek7KHnKjz6RNvxCmbJ6hDBbtcC1MWArUuiZdEEnUyKRY97ugdH9T7c3W5wKx3yZAEe8iJyH9dLjA5yzoRtVeXC8zzt+PkfciJyF9dLjBrnBLiDFqYdF3uqxNRB3W51Khxenk4TkTt0qUCUwgBK4cUEVE7danAPOuW4RXNww+IiPzVpQLz/CWRPCQnovboUoFpPTekKJmzrBNRO3SpwGz0yIgzaGHWd6mvTUQB0qVO5g1Ls2BQsjnUZRBRhOpSgRln1CEu1EUQUcRSJTDr6urwhz/8ASdPnoTRaETPnj2xcOFCJCUlYc+ePZg3bx5cLheysrKwePFiJCcnq1EWEZFfVDmZp9FocM8992DTpk34+OOP0aNHDyxZsgSyLOOJJ57AvHnzsGnTJuTn52PJkiVqlERE5DdVAjMhIQHDhw/3Pc7NzUV5eTn2798Pk8mE/Px8AMCMGTOwceNGNUoiIvKb6t3FsizjvffeQ0FBASoqKpCZmelblpSUBFmWUV9fr3ZZRERtUr3T59lnn0VUVBTuvPNOfPrppwFZZ3JyTEDWQ0SdW2pqbIfer2pgFhcX48SJE1ixYgW0Wi0yMjJQXl7uW15bWwutVouEhAS/1mu1NkGWhe9xRzcKEXVO1dWNFz3nT16odkj+8ssvY//+/Vi+fDmMRiMAYODAgXA6ndi5cycAYNWqVZgwYYJaJRER+UUjhBBtv6xjjh49iqKiImRnZ8Nsbh443r17dyxfvhy7du3C/PnzWwwrSklJ8Wv9bGESkRIdbWGqEpjBxsAkIiUi5pCciCjSMTCJiBRiYBIRKcTAJCJSiIFJRKQQA5OISCEGJhGRQgxMIiKFGJhERAoxMImIFGJgEhEpxMAkIlKIgUlEpBADk4hIIQYmEZFCDEwiIoUYmERECjEwiYgUYmASESnEwCQiUoiBSUSkEAOTiEghBiYRkUIMTCIihRiYREQKMTCJiBRiYBIRKcTAJCJSiIFJRKQQA5OISCEGJhGRQgxMIiKFGJhERAoxMImIFGJgEhEpxMAkIlKIgUlEpBADk4hIIVUCs7i4GAUFBejXrx+OHDnie760tBTTp09HYWEhpk+fjuPHj6tRDhFRu6gSmOPGjcPKlSuRlZXV4vn58+dj5syZ2LRpE2bOnIl58+apUQ4RUbuoEpj5+fnIyMho8ZzVasXBgwdRVFQEACgqKsLBgwdRW1urRklERH7Th+qDKyoqkJ6eDp1OBwDQ6XRIS0tDRUUFkpKS/FpXcnJMMEokok4mNTW2Q+8PWWAGktXaBFkWvscd3ShE1DlVVzde9Jw/eRGyXvKMjAxUVlZCkiQAgCRJqKqquujQnYgoXIQsMJOTk5GTk4P169cDANavX4+cnBy/D8eJiNSiEUKItl/WMc899xw2b96MmpoaJCYmIiEhAZ988glKSkowd+5cNDQ0IC4uDsXFxejVq5ff6+chOREp0dFDclUCM9gYmESkRMSewyQiijQMTCIihRiYREQKMTCJiBRiYBIRKcTAJCJSiIFJRKQQA5OISCEGJhGRQgxMIiKFGJhERAoxMImIFGJgEhEpxMAkIlKIgUlEpBADk4hIIQYmEZFCDEwiIoUYmERECjEwiYgUYmASESnEwCQiUoiBSUSkEAOTiEghBiYRkUIMTCIihRiYREQKMTCJiBRiYBIRKcTAJCJSiIFJRKQQA5OISCEGJhGRQgxMIiKFGJhERAoxMImIFGJgEhEpxMAkIlIoLAKztLQU06dPR2FhIaZPn47jx4+HuiQioouERWDOnz8fM2fOxKZNmzBz5kzMmzcv1CUREV0k5IFptVpx8OBBFBUVAQCKiopw8OBB1NbWKl6HVqtp8R8RUWt+mRX+5oU+SHUpVlFRgfT0dOh0OgCATqdDWloaKioqkJSUpGgdiYnRwSyRiDqJ5OSYDr0/5C1MIqJIEfLAzMjIQGVlJSRJAgBIkoSqqipkZGSEuDIiopZCHpjJycnIycnB+vXrAQDr169HTk6O4sNxIiK1aIQQItRFlJSUYO7cuWhoaEBcXByKi4vRq1evUJdFRNRCWAQmEVEkCPkhORFRpGBgEhEpxMAkIlKIgUlEpFCXCUy1Jvioq6vDvffei8LCQkyePBkPPvig7zLPfv36YfLkyZg6dSqmTp2Kw4cP+963ZcsWTJgwATfddBMeeeQROByODi+7nIKCAkyYMMFXy1dffQUA2LNnD6ZMmYLCwkLcfffdsFqtvvcEY1lrTp065atr6tSpKCgowLBhwy5btxr1FRcXo6CgAP369cORI0d8z19u31J7WWs1Xm6fBNTdLy+1DdX+Xf3dJ31EFzFr1izx0UcfCSGE+Oijj8SsWbOC8jl1dXXiu+++8z1+6aWXxH/9138JIYTo27evaGpquug9TU1NYuTIkaK0tFQIIcRTTz0lli1b1qFlbRk7dqw4fPhwi+ckSRLjx48XO3bsEEIIsXz5cjF37tygLVPqueeeE88888wl61arvh07dojy8vKLarjcvqX2stZqvNw+KYS6++WltqGav2tH9skuEZg1NTViyJAhwuv1CiGE8Hq9YsiQIcJqtQb9szdu3Ch++9vfCiEuvWNu2LBB3Hfffb7H+/btExMnTuzQsra0toPu3btXTJo0yffYarWK3NzcoC1TwuVyieHDh4v9+/dfsm6167uwhsvtW2ovu1SNv3ThPilEaPZLpYEZbvtkyCffUEMgJvhoD1mW8d5776GgoMD33KxZsyBJEq6//no89NBDMBqNqKioQGZmpu81mZmZqKio8NXenmVKPP744xBCYMiQIXj00UcvWl9SUhJkWUZ9fX1QliUkJLRZ45YtW5Ceno4BAwZcsu64uLiQ1Xe5fUsIoeoyJftya/skEB77pVq/a0d+8y5zDjMUnn32WURFReHOO+8EAHzxxRdYs2YNVq5ciWPHjmH58uUhq23lypVYt24dVq9eDSEEFi5cGLJaLmf16tWYNm2a73Gk1B2ufrlPAuGxX0bK79olAjMUE3wUFxfjxIkTeOWVV6DVan11AEBMTAzuuOMO7Nq1y/d8eXm5773l5eW+17Z3WVvOv85oNGLmzJnYtWvXReurra2FVqtFQkJCUJa1pbKyEjt27MDkyZMvW3dr20KN+s5/7qX2LbWXtaW1ffLCbRrK/VLN37Ujv3mXCEy1J/h4+eWXsX//fixfvhxGoxEAcPbsWTidTgCA1+vFpk2bkJOTAwAYM2YMfvjhB19v56pVq3DLLbd0aNnl2O12NDY2AgCEENiwYQNycnIwcOBAOJ1O7Ny507e+CRMmAEBQlrXlww8/xA033IDExMTL1h2q+oDL71tqL7uc1vZJIDz2S7V/14785l3mWnK1Jvg4evQoioqKkJ2dDbPZDADo3r077rnnHsybNw8ajQZerxd5eXl46qmnEB3dPPnxZ599hsWLF0OWZeTk5OCll15CVFRUh5ZdSllZGR566CFIkgRZltG7d2/893//N9LS0rBr1y7Mnz8fLpcLWVlZWLx4MVJSUgAgKMsup7CwEE8//TSuv/76NutWo77nnnsOmzdvRk1NDRITE5GQkIBPPvnksvuW2staq/GVV15pdZ9cvnw5du/erep+2Vp9K1asUP13be8+2WUCk4ioo7rEITkRUSAwMImIFGJgEhEpxMAkIlKIgUlEpBADM0LMmzevxRUYf//73zFy5Ejk5eWhrq4uhJW1LVxq3blzJwoLC0P2+ZMmTcL27dsDsq5169bh7rvv9j3u168fTpw4EZB1A0BeXh7KysoCtr5OQ9EV59Rhffv2FcePH2/x3GuvvSYee+wxv9fldrvFNddcI3788cdAleeX1atXixkzZih6baBq/e6778SYMWP8ek9r2zwYysrKRN++fUVubq7Izc0VI0aMEPfdd5/4+uuv270uj8fj1/s68l3vvPNO8f7777frvV0NW5gRyGq1wuVy4aqrrvL7vUIIyLIchKpa15FaI82OHTuwe/durF27FiNHjsSDDz6INWvWBPxzvF5vwNdJyjAww8T27dtx/fXX4y9/+QtGjBiB0aNHY/Xq1b7lc+fOxdKlS1FaWuq7jGvo0KGYPXs2gOYrF6ZNm4YhQ4Zg2rRpvmtxgeaZaJYuXYoZM2bg2muvRVlZGfr164eVK1fi5ptvRl5eHl555RWcPHkSM2bMwODBgzFnzhy43W5FtRcUFODtt9/G5MmTMWTIEDzyyCNwuVyXrLWkpAR33XUXhg0bhsLCQmzYsMG3rq1bt2LixInIy8vDmDFj8Pbbb8Nut+Pee+9FVVUV8vLykJeXh8rKSuzbtw/Tp09Hfn4+Ro8ejYULF/pq/rd/+zcAwNSpU5GXl4cNGzb4tvF5JSUlmDVrFvLz8zFp0iR8/vnnLbb3M888g/vuuw95eXm44447cPLkSUXbIzU1Fb/97W/x4IMPYsmSJb5/oAoKCrBt2zYAwL59+3D77bdj8ODBGDlyJF588UUA8E2KMXToUOTl5WH37t1Ys2YNZsyYgRdeeAHDhw/HsmXLsGbNGvzmN79p8blbt27FuHHjMHz4cBQXF/s+d9myZXj88cd9rzt16hT69esHr9eLpUuXYufOnVi4cCHy8vJ8k15ceIjf2NiIP/zhD7juuuswduxYvPHGG751n6+juLgYQ4cORUFBAbZu3apoO0WkUDdxu4q2Dsm/++47kZOTI1555RXhdrvFF198IQYNGiTq6+uFEEI8+eST4uWXXxZCXHzYVldXJ/Lz88WHH34oPB6P+Pjjj0V+fr6ora0VQjQfct1www3iyJEjwuPxCLfbLfr27Svuv/9+0djYKI4cOSIGDBggZs+eLU6ePCkaGhrELbfcItasWdPqd/nlIfnYsWPFtGnTxJkzZ0RdXZ2YMGGC+Pvf/95qrTabTVx//fXiH//4h/B4POLAgQNi2LBh4ujRo0IIIUaNGuWb2LW+vt43D2Zrh+Q//PCD2L17t/B4PKKsrExMmDBBvPPOO5fc5heuw+12i/Hjx4s333xTuFwusW3bNpGbmytKSkp823vYsGFi7969wuPxiEcffVQ88sgjrW6PSx1Gnzx5UvTt21ccO3bMt52++eYbIYQQv/71r8WHH34ohGiecHf37t2XXNfq1atFTk6OePfdd4XH4xEOh+Oi36Bv377izjvvFHV1deL06dPi5ptv9h1m//LUzy8/o7VD8gu33RNPPOHbV8rKylqse/Xq1aJ///7if//3f4XX6xUrV64Uo0aNErIst7qtIh1bmGFEr9fjgQcegMFgwA033ICoqCiUlpa2+b4vvvgCPXv2xK233gq9Xo+ioiL06tUL//znP32vue2229CnTx/o9XoYDAYAwD333IOYmBj06dMHffv2xahRo9CjRw/Exsbi+uuvx8GDBxXXPmvWLKSnpyMhIQFjx47Fjz/+eMlas7KyMG3aNOj1evTv3x+FhYXYuHGjbxscO3YMTU1NiI+PbzEP5i8NHDgQubm50Ov16N69O6ZPn44dO3Yoqnfv3r2w2+247777YDQaMWLECIwdOxaffPKJ7zXjx4/HoEGDoNfrMWXKlEt+p0s5fy10fX39Rcv0ej1OnjyJ2tpaREdHIzc3t811zZo1C3q93nc9+C/de++9SEhIQGZmJmbPnu2boKMjJEnChg0b8NhjjyEmJgbdu3fHXXfdhXXr1vlek5mZiV//+tfQ6XS47bbbUF1djZqamg5/djhiYKpEp9NddO7J6/X6wgsAEhISoNf/PKezxWKB3W5vc91VVVUtJkQFmnfiyspK3+PWpte6cLIBk8l00WMln31eamqqorpPnz6Nffv2IT8/3/ffxx9/jOrqagDAa6+9hq1bt2Ls2LG48847sXv37kt+ZmlpKX7/+99j1KhRGDx4MJYuXaq4F76qqgrdunVrMc3ZL7fZhdvDbDb7tT0A+NbV2rRhzz//PI4fP45bbrkF06ZNa/GPW2u6devW5udd+BtnZWWhqqrKr3pbU1dXB4/Hc9GEwJfaThaLBQD83laRokvMuB4OMjIycOrUKfTu3dv33KlTp5Cdnd3hdaelpbWY3w9onvV6zJgxvscajabDnxMIGRkZGDp0KN55551Wlw8aNAhvvvkmPB4PVq5ciUceeQRbt25ttf4FCxagf//++NOf/oSYmBj89a9/xaZNmxTVkZaWhjNnzkCWZV9oVlRUBOT3OO/TTz9FcnIyrrzyyouWZWdn4+WXX4Ysy9i8eTMefvhhbN++/ZK/k5Lfr6KiAn369AHQPAfl+RauxWLxTeEGwK/WX2JiIgwGA8rLy30dd+dnme+K2MJUycSJE/Hmm2/6/ki3bduGLVu2BGRc4A033IDjx4/j448/htfrxYYNG3Ds2DHceOONHS88wG688UYcP34cH330ETweDzweD/bt24eSkhK43W6sW7cOjY2NMBgMiI6O9oVZcnIy6uvrffMmAoDNZkN0dDSio6NRUlKC9957r8VnpaSkXHIs4aBBg2A2m/HWW2/B4/Fg+/bt2LJlCyZOnNjh71hTU4O//e1veP311/Hoo4+2aMWet3btWt/EtXFxcQAArVaLpKQkaLXado2BfPvtt3H27FlUVFTg3Xff9X2XnJwc7NixA+Xl5WhsbMSf//znFu+73HbS6XSYMGECli5diqamJpw+fRrvvPMOpkyZ4nd9nQEDUyUPPPAA8vLyMHPmTAwdOhSLFy/GkiVL0Ldv3w6vOzExEStWrMA777yD4cOH46233sKKFSuCer+i9oqJicHbb7+NDRs2YMyYMRg9ejSWLFni691eu3YtCgoKMHjwYKxatQqLFy8GAPTu3RuTJk3C+PHjkZ+fj8rKSjz55JNYv349Bg8ejD/+8Y8Xhd2DDz6IuXPnIj8/v0VPPNA8s/eKFSvw5Zdf4rrrrsMzzzyDRYsWtTgC8NfQoUORm5uLyZMnY+vWrXj11Vfxq1/9qtXXfvXVV5g0aRLy8vLw/PPPY+nSpTCbzbBYLLj//vvxm9/8Bvn5+dizZ4/izx83bhxuv/123Hrrrbjxxht9nz1q1ChMnDgRU6ZMwe23346xY8e2eN/s2bOxadMmDB06FM8999xF6/3jH/8Ii8WC8ePHY+bMmSgqKmpx25CuhPNhEhEpxBYmEZFCDEwiIoUYmERECjEwiYgUYmASESnEwCQiUoiBSUSkEAOTiEghBiYRkUL/HxOf3kxKM75MAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "ax = sns.displot(infest_stem_unif,\n", " kde=True,\n", " bins=100,\n", " color='skyblue',\n", " linewidth=15,\n", " alpha=1)\n", "ax.set(xlabel='Uniform Infestation Distribution', ylabel='Frequency')" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Size of cluster: 7842 stems\n" ] } ], "source": [ "from scipy.stats import gamma\n", "\n", "# Infest stems using gamma distribution to create cluster of pests.\n", "infest_clust = gamma.rvs(3, scale=500, size=40000) #rvs(a, loc=0, scale=1, size=1, random_state=None)\n", "infest_clust = [round(x) for x in infest_clust]\n", "\n", "# Insert infested cluster of stems into shipment at random position\n", "start_clust = np.random.randint(1, high = num_stem - max(infest_clust), size = 1)\n", "infest_stem_clust = infest_clust + start_clust\n", "\n", "print(\"Size of cluster:\", int(max(infest_stem_clust) - min(infest_stem_clust)), \"stems\")" ] }, { "cell_type": "code", "execution_count": 10, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 10, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAV0AAAFtCAYAAACk8rkaAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/Il7ecAAAACXBIWXMAAAsTAAALEwEAmpwYAABB7klEQVR4nO3deXhTVd4H8O+92ZsuaUspaYsUkWIBlUqBGRQdClJwgVd9kQ6C8+I+M6IOA8roSGeKiCwDCuIyrww+PnZgFtGCIqCiqDOvCMMmylq2Qgtt0y1p9nvP+0ea2NYuaZrc3CS/z/PMMzb3Njknab+c/u4553KMMQZCCCGS4MPdAEIIiSUUuoQQIiEKXUIIkRCFLiGESIhClxBCJEShSwghElKGuwHhZDJZIIrdz5hLTo5Dfb1VghbJE/U/dvsfy30Huu9/WlpCj5+TRrp+UCoV4W5CWFH/Y7f/sdx3IDT9p9AlhBAJUegSQoiEKHQJIURCFLqEECIhCl1CCJEQhS4hhEgopufpylEg8/7aq6kxB6ElhJBQoJEuIYRISLKR7rJly7Bjxw5cvHgRW7duRU5ODi5cuIBf//rXvnPMZjMsFgu++eYbAEBBQQHUajU0Gg0AYP78+Rg3bhwA4ODBg1i0aBEcDgcyMzOxYsUKpKamStUdQggJiGShO2HCBNx333249957fY9lZWWhrKzM9/WSJUsgCEKb71uzZg1ycnLaPCaKIhYsWIClS5ciPz8fr776KlauXImlS5eGthOEENJLkpUX8vPzYTQaOz3udDqxdetW3H333d0+15EjR6DRaJCfnw8AKCoqwvbt24PWVkIICRXZXEjbtWsX0tPTMWzYsDaPz58/H4wxjBw5EvPmzUNiYiKqqqqQkZHhOyclJQWiKKKhoQEGg8Hv10xNjff73GBc4JJKKNoaSf0PhVjufyz3HQh+/2UTuu++++6PRrmlpaUwGo1wOp1YsmQJSkpKsHLlyqC9pr+7jKWlJUg2IyDYH3Aw2i1l/+Uolvsfy30Huu9/xO4ydvnyZezduxd33HFHm8e95Qi1Wo2ZM2di//79vscrKyt959XV1YHn+R6NcgkhJBxkEbrvvfcebr75ZiQnJ/ses1qtMJs9/8IwxrBt2zbk5uYCAIYPHw673Y59+/YBADZt2oTJkydL33BCCOkhycoLzz//PHbu3Ina2lrMmTMHBoMBH374IQBP6D777LNtzjeZTJg7dy4EQYAoihg0aBCKi4sBADzPY/ny5SguLm4zZYwQQuSOY4x1X9SMUuGu6UpxgYJqur0Xy/2P5b4DUVzTJYSQWEGhSwghEqLQjSCMMYixWw0iJCrIZp4u6drpJic+vmABBw53DkxAmo4+OkIiEY10I0C1zY1/lDeBBweHIOKdE42wusRwN4sQEgAKXZljjGFHhQVaJYdZOUkouioJDpFhf6093E0jhASAQlfmTjQ6cbHZjfEZeuiUPNJ0SgxKVGF/rQ0uP6e7xfraeULkhEJX5vZW25Ck5jE8ReN7bFRfHaxuhpONzjC2jBASCApdGauyunCh2Y2RaTrwHOd7/Ip4FXRKDuUUuoREHApdGTtYa4eKB65L1bR5nOc4XJmgxukmJ00hIyTCUOjKlFtkONbgRE6SBhrFjz+mQUlq2ASGKqs7DK0jhASKQlemTjU54RBYm1pua1cmqMABKG+iEgMhkYRCV6a+q3NAr+QwIEHV4XGtkkd6nBIXLDTSJSSSUOjKkM0torzJiaHJmjYX0NrL0itR2eyCQHVdQiIGha4MHa13QGTA8BRtl+dl6VVwM6Ca6rqERAxawC8xfxYqfFfvQJpWgb46RZfnZcZ7Pr4LzW4Y9R2XIQgh8kIjXZlpcgq42OzG0GQNuC5KCwCQoFIgSc3jQrNLotYRQnqLRroS8Xcp7qmWBQ85BrVf52fpVThrdoIx1m1IE0LCj0a6MnOi0YkUjQKpWv/+PczUK9HsZmh00q5jhEQCCl0ZsQsizltcGJzk3ygXALLiPbVcKjEQEhkodGXkbJMLIgOu6kHo9tEqoOE5mq9LSISg0JWRM2YnNDyHDL3/pXae85x/kUa6hEQECl2ZYIzhjNmFKxJUUPTwglhWvAo1dgF2N9V1CZE7Cl2ZqHeIaHKKGNjJst+uZOl/mK9LCJE3mjImE2fMnqliAxP9r+d6ZehV4DnggsXVaT24oylrNTXmHr8WIaR3aKQrExUWFxJVPJI1Xa9C64iK52CMU6KC6rqEyB6FrgwwxnDB4vZN/wpEll6FKqvbr/umEULCh0JXBhqdIixuEf3jA6/2ZMUrITKgiuq6hMgaha4MVFg8ZYGsXmxa4/1eKjEQIm8UujJwodkFrYJDH23P67leOiWPNK0CFywUuoTIGYWuDFQ2u5GhV/Z6w5qseBUuNrvpZpWEyBiFbpi5RIZau4B+cb2fvddfr4JTZKi2CUFoGSEkFCh0w+yy1Q0GoJ+u96Gb1XIhroJKDITIFi2OCCF/9tC9ZPPMNjAGYaSbqPZsan7e4sKovrpePx8hJPgkG+kuW7YMBQUFGDJkCE6cOOF7vKCgAJMnT8a0adMwbdo0fPnll75jBw8exNSpU1FYWIj7778fJpPJr2OR5JLVDb2SQ7wqOB9FdoIK5y0uqusSIlOShe6ECRNQWlqKzMzMHx1bs2YNysrKUFZWhnHjxgEARFHEggULsGjRIuzYsQP5+flYuXJlt8cizSWrG/3ien8RzSs7QQ2HwHCJblZJiCxJFrr5+fkwGo1+n3/kyBFoNBrk5+cDAIqKirB9+/Zuj0USt8hgsgtID0I912tAy6q2s2aq6xIiR7Ko6c6fPx+MMYwcORLz5s1DYmIiqqqqkJGR4TsnJSUFoiiioaGhy2MGg8Hv101Njff7XH/vcdYTJrsABiAtiKEbp+LRV6fAObMLY/t1fW5P+hSK/keSWO5/LPcdCH7/wx66paWlMBqNcDqdWLJkCUpKSiQrFZhMFoh+7FWQlpYQ0I5c3X1YtXZPCaA3iyI6kp2gxn9qbHCJDCq+87KFv30KtP/RIpb7H8t9B7rvfyCBHPYpY96Sg1qtxsyZM7F//37f45WVlb7z6urqwPM8DAZDl8ciSa1dAA8gJYCdxboyIF4FgYFWpxEiQ2ENXavVCrPZ868IYwzbtm1Dbm4uAGD48OGw2+3Yt28fAGDTpk2YPHlyt8ciSa1dQLJWAUUXo9FA9I/37K97rpu6blpawo/+RwgJLcnKC88//zx27tyJ2tpazJkzBwaDAa+//jrmzp0LQRAgiiIGDRqE4uJiAADP81i+fDmKi4vhcDiQmZmJFStWdHssktTa3egbxHqul1rBISNOSRfTCJEhjrHYndAZzpquS2RYdciEsf10GGfU9/i5u/NVlRVfXbLiiWtSoFP6/wdNR/2kul7s9j+W+w5EaU03VtW1zFzoow3NHxvZLfdaO0d1XUJkhUI3TOocnk1pgn0RzcuoV0LNc93WdQkh0gr7lLFo5M+fHA0toRvIPdH8oeA49I9X4mzLDS8JIfJAI90wqXcI0Cs5qBXBnbnQWnaCGvUOEY1O2uqRELmg0A2TeqcAQ4hGuV4DWuq656nEQIhsUOiGSYNDDFlpwStNq4CG53CRblZJiGxQ6IaBS2Qwu0QY1KENXY7jkKFX4iLdrJIQ2aDQDYNG30W00L/9mXoVauwCHIIY8tcihHSPQjcM6p2hnbnQWqbeM0GlikoMhMgChW4Y1Ds8o85QX0gDPPN1AeAChS4hskChGwaNTgFqnoMuhNPFvLQKHqlaBd1JghCZoNANgyaniEQ1H7Rb9HSnn06JyzYKXULkgEI3DJqcAhLV0r31fXUKmF0irC66mEZIuFHohkGTS0SiKvT1XK9+Lbd3p9EuIeFHoSsxl8hgczNJR7reG19SXZeQ8KPQlVhTy3QxKUNXq+SRpOZppEuIDFDoSqzJ6amrJoZ4NVp76XQxjRBZoNCVWFPLxaxElbRvfZpOgQaHCLcfd8oghIQOha7EvOWFBAnLC4DnDhUMgMlO2zwSEk4UuhJrcopIUPFQSDRH16uP1lPOqLVTiYGQcKLQlZg3dKWWrFGAA410CQk3Cl2JmV2i5KUFAFDyHFI0CtRQ6BISVhS6ErO4RMSHYaQLAKlaBY10CQkzCl0JOQQRTpGFpbwAAH10CtQ7BJrBQEgYUehKqNnlCbtwjXS9Mxi8t38nhEiPQldCZpcn7MIXut4ZDBS6hIQLha6ELC0LI8JVXkhpmcFQSyvTCAkbCl0JeUM3XCNdJc8hWaOgkS4hYUShKyGzS4Sa56BRhO9t76Ol0CUknJThbkAsCed0Ma8+WgVONjrhFhmU/I9XxaWlJfj+u6bGLGXTCIkJNNKVkDxCl2YwEBJOFLoSkkXo6mgGAyHhRKErEcYYLK7w7LvQWkrLbd9NtPENIWFBoSsRh8DgZoA+zKGr5DkY1DzqaKRLSFhIdiFt2bJl2LFjBy5evIitW7ciJycH9fX1eOqpp3D+/Hmo1WoMGDAAJSUlSElJAQAMGTIEOTk54HlPUC1fvhxDhgwBAOzatQvLly+HIAgYNmwYli5dCp1OJ1V3eqzZ3TJdTBn+f+dStQqYqKZLSFhIlgATJkxAaWkpMjMzfY9xHIcHH3wQO3bswNatW9G/f3+sXLmyzfdt2rQJZWVlKCsr8wVuc3MznnvuObz++uv4+OOPodfrsX79eqm6EhDvHF29Stp9dDuSqlWizi5AZLQHAyFSkyx08/PzYTQa2zxmMBgwZswY39cjRoxAZWVlt8/1xRdfYPjw4cjOzgYAFBUV4aOPPgpqe4Ot2e0JOL0cRroaBdzsh/u1EUKkI5t5uqIoYuPGjSgoKGjz+OzZsyEIAm666SbMnTsXarUaVVVVyMjI8J2TkZGBqqqqHr9mamq83+e2nr8aiGbfSFcGoav1XkwTYNB0foPM1n3ubf8jXSz3P5b7DgS//7IJ3cWLFyMuLg6zZs3yPfb555/DaDTCYrFgwYIFWLduHX7zm98E7TVNJgtEP7Y5TEtL6NFCgY4+pGaXCJ4DtAo5lBdaQtchYFAX53n73NP+R5tY7n8s9x3ovv+BBHL4h13wXGQ7d+4cXnrpJd9FMwC+ckR8fDymT5+O/fv3+x5vXYaorKz8UelCbprdIuKVPDiJ743WEZ2Sh07J0bQxQsIg7KG7atUqHDlyBOvWrYNarfY93tjYCLvdDgBwu93YsWMHcnNzAQDjxo3Dt99+i7NnzwLwXGybMmWK5G3viWaXKIvSgleqhu4iQUg4SFZeeP7557Fz507U1tZizpw5MBgMeOmll/DGG28gOzsbRUVFAICsrCysW7cOp0+fxqJFi8BxHNxuN/Ly8vDEE08A8Ix8S0pK8Mgjj0AUReTm5uLZZ5+VqisBsbhFJKo6r59KLVWrwIlGZ5fnxHotj5BQ4BiL3XlDUtZ0135rwlVJaky5Qh5B9k21DbsuNuPxa1IQ58eMCqrrxWb/Y7nvQBTXdKOdyBisbiaL6WJeqZofZjAQQqQjnxSIYjY3A4M8pot5eWcw0HJgQqQlnxSIYt4lwHIK3SQ1DyUH1NIMBkIkJZ8UiGK+hREyKi9wHIcUrYL21SVEYvJJgSjmG+nKKHQBT12X9tUlRFrySoEoZW3ZdyFOGf6FEa2lapVodIpw+TGDgxASHBS6ErC6PUuANTJYAtwaXUwjRHqy2Xsh0nU1X8/qEqGXyRLg1nyh6xCQHkc/CoRIgUa6Emh2i7IrLQBAcg/m6qalJdAKNUKCgEJXAlY382vVl9RUPIckNU8b3xAiIfklQRSyukVZhi4A9KFb9xAiKXkmQZSxyrS8AHjuDky37iFEOhS6IeYUGFyivFajtdZHq6Rb9xAiIXkmQRSxtiyM0Mm1vKDzXEyrobouIZKQZxJEEZtMV6N59WmZNlZro7ouIVKQZxJEkWaZrkbz0ih4JKp41NACCUIkQaEbYt7yglxnLwBAmk6BGhuVFwiRgnyTIEpEQuj20SpR5xAg0AwGQkJOvkkQJaxuBhUPqGW270JraToFBAbU03xdQkKOQjfEml3yXRjh1Ufr2XeBLqYREnq0y0kvdbcfgZxXo3mlahXg4Jk2djU04W4OIVFN3mkQBeS8Gs1LxXNI1ihQQyNdQkKOQjfE5LrZTXt9tHQXCUKkIP80iGCMMVjdomwXRrSWplOg3iHQXSQICTG/0+CTTz6B201zOXvCITIIDIiT6b4LraVplWDwb29dQkjg/E6DNWvW4MYbb0RJSQkOHToUyjZFDZvMV6O15t2DgW7JTkho+R26W7ZswVtvvQWNRoO5c+eisLAQr776Ki5cuBDK9kU0763XI6Gmm6xRQMGBLqYREmI9SoOrr74aTz/9NHbv3o3i4mJs374dt9xyC+69915s2bIFokjbA7YWCavRvBQch1StAtW0HJiQkOrxPN3z589jy5Yt2LJlCziOw+OPPw6j0YjS0lLs3LkTr7zySijaGZG8t17XR0B5AQD6xSlxssEJxpjsbqJJSLTwO3RLS0tRVlaGc+fOYcqUKVi+fDlGjBjhO15YWIixY8eGoo0RS+576bbXT6fEYZMDTS4RSWpFuJtDSFTyO3S/+OILzJkzBxMmTIBarf7RcZ1Oh7Vr1wa1cZGu2S1Co+Cg5CNj1Niv5Tbsl6xuCl1CQsTv0F2zZg14nodKpfI95nK5wBjzhfCNN94Y/BZGMKtL/qvRWkvTKcEBuGx1Y4ih4+XA3mXPNTVmCVtGSPTw++/e+++/H999912bx7777js88MADQW9UtLC6WUQsjPBS8Rz6aBW4ZKWLaYSEit+JcPz4cVx33XVtHrv22mtx7NixoDcqWljdYsTUc736xSlxyeYGo711CQkJvxMhMTERtbW1bR6rra2FTqfr9nuXLVuGgoICDBkyBCdOnPA9fubMGcyYMQOFhYWYMWMGzp492+tjcmITGHQRVF4APKFrdTOYXTT9j5BQ8Dt0J02ahN/+9rc4ceIEbDYbjh8/jqeffhpTpkzp9nsnTJiA0tJSZGZmtnm8uLgYM2fOxI4dOzBz5kwsWrSo18fkgjEGm1tEnCLyRroAqMRASIj4nQi/+c1vMGjQIEyfPh3XX389ZsyYgYEDB2LevHndfm9+fj6MRmObx0wmE77//nvcfvvtAIDbb78d33//Perq6gI+JicuERAYIm6k29d7MY0WSRASEn7PXtBoNCguLsaiRYtQX1+P5OTkXk2gr6qqQnp6OhQKz9QkhUKBvn37oqqqCoyxgI6lpKT0qA2pqfF+n9vdZuXtRdocXS8V71mZ1t1It6fvR6SLtf62Fst9B4Lf/x6tSDObzThz5gyam5vbPP7Tn/40qI2SislkgejHVoZpaQmdTpHq7AOxC57njbSRLuApMZxpcnZ5TixNGevq8492sdx3oPv+BxLIfofu5s2bUVJSgri4OGi1Wt/jHMfh008/7fELG41GXL58GYIgQKFQQBAEVFdXw2g0gjEW0DE58Y10I6ymCwDpOiWO1DlgdgpIoEUShASV36G7evVqvPzyy7j55puD8sKpqanIzc3FBx98gGnTpuGDDz5Abm6ur0QQ6DG5sEXQZjftZeg9PxaVVjeGdBK6rf+Fj+WRECE9xTE/J2SOHTsWX375pa+W2hPPP/88du7cidraWiQnJ8NgMODDDz9EeXk5Fi5ciKamJiQmJmLZsmW48sorASDgYz0RyvLCvmobPrnYjCeuSYm4uq5bZFh12ITRaTr8LFPf7fnRHrqx/Cd2LPcdCE15we/Q3bBhA5qbm/GrX/0KPB9ZIdKZUIbuF1XN+L9LNjw1IjUid+x663gD1DyHmYOTuj032n8pYzl4YrnvQJhrum+99RZqa2vx5ptvwmAwtDn2+eef9/iFo53NzaBVchEZuACQEeep64qMgY/QPhAiR36H7ooVK0LZjqhjc4sReRHNK0OvxP5aO2rtAvrqerztMiGkE37/No0ePTqU7Yg6NjeLqB3G2suI8+wmV9Xs7jZ0aecxQvzn91DM6XRi9erVmDBhAkaOHAkA+Oqrr/DOO++ErHGRLBI3u2ktWcNDq+BQaXWFuymERBW/U+GFF17AiRMnsHLlSl+dcvDgwdi4cWPIGhfJ7AKDThG5I12O45ARp0RlMy0HJiSY/C4vfPLJJ9i5cyfi4uJ8sxfS09Nx+fLlkDUuUjHGIn6kCwBGvRKnL9ngEERoIrg+TYic+P2bpFKpIAhtb89dV1f3o5kM5IfNbiK5pgv8UNelHccICR6/Q3fy5Ml4+umnUVFRAQCorq5GSUkJbrvttpA1LlLZBM9qNG2Ej3R9K9OoxEBI0PRoa8esrCxMnToVTU1NKCwsRN++ffHrX/86lO2LSLaWW69H+khXp+SRrOFRSSNdQoLG75quWq3GM888g2eeeQZ1dXW93toxmkXyZjftZcSpcM7suQEpfd6E9J7foestK3i13t6xf//+wWtRFLD59tKN/JDK0CvxXb0DZpeIRNpxjJBe8zt0b7nlFnAc1+aGhd6Rz9GjR4Pfsgj2Q3khGka6P9R1KXQJ6T2/Q7f9XX9ramrwyiuvID8/P+iNkrvuNrmwCiI4AJoInqfr1VenhILzbPN4dbIm3M0hJOIFPBRLS0vDs88+i1WrVgWzPVHB7mbQKrio2ChGwXNI1ylR2Uwr0wgJhl79/Xv69GnYbLZgtSVqRMPCiNYy9Epcsroh+LcLKCGkC36XF2bOnNnm6rXNZsOpU6doylgHIn2zm/Yy4lTYx+yosQm+W7QTQgLj92/Q9OnT23yt0+lw9dVXIzs7O9htinhWt4gkTfRcdPIukqiyuroMXbqFDyHd8zt077zzzlC2I6rYBQZjFFxE80pS84hTcqhsdiOvT7hbQ0hk8zt0X375Zb/Oe+KJJwJuTDSIls1uWuM4DsY4Ja1MIyQI/A7dc+fOYefOnRg+fDgyMzNRWVmJb7/9FpMmTYJGQ1OJvLyb3UTDwojWMvQqlDdZYRdEaKNgpR0h4eJ36DLG8Kc//QmFhYW+x3bu3Int27dj6dKlIWlcJPJudhNNI13gh0USl5rdyE5Uh7k1hEQuv5Phiy++wMSJE9s8VlBQgN27dwe9UZHMuxotkjcw74jRuzKNSgyE9IrfoTtgwACUlpa2eWzjxo244oorgt6oSObb7CbKRrpaJY9UjQIXaZEEIb3id3nh+eefx2OPPYY333zTd8cIpVKJtWvXhrJ9EcfuHelGWU0X8EwdK29y0o5jhPSC36E7dOhQ7NixA4cOHUJ1dTXS0tIwYsQIqFSqULYv4tiF6NnWsb0MvRLf1jnQ6BRhiKJ5yIRIKeBkGDVqFFwuF6xWazDbE/Fsgmekq42ymi4AGFtu30N1XUIC5/dI9/jx4/jlL38JtVqNy5cv49Zbb8XevXvx3nvv4aWXXgphEyOL3S1CzXNQ8NEXumk6BZQcUNXswlDacYyQgPg90v3DH/6Axx9/HNu3b4dS6cnqUaNG4T//+U/IGheJbAKLylEuACg4Dv3ilKiikS4hAfM7dE+dOoVp06YB+GHz8ri4ODgcjtC0LELZ3QzaKLyI5mWMox3HCOkNv0M3MzMTR44cafPY4cOHacpYO3ZBjMqLaF5GvQpuBtTahHA3hZCI5HdN94knnsAjjzyCoqIiuFwuvPHGG9i0aRMWL14cyvZFHJvA0EcbvaHru32P1YV02uaRkB7zOx3Gjx+PN998E3V1dRg1ahQuXryItWvX4sYbbwxl+yKO3S1GbU0X8Ow4plNwqGruuq6blpbQ7W2NCIlFfg1VBEFAYWEhtm3bhj/84Q8hblLkYozBLrCoLi9wHAejni6mERIov9JBoVBAoVDQRbNueHcYi+YLaYDnYlqNXYCjZSEIIcR/fhfl7rvvPjz55JN45JFH0K9fvzbLQPv37x9wAy5cuNDmlj9msxkWiwXffPMNCgoKoFarfVtHzp8/H+PGjQMAHDx4EIsWLYLD4UBmZiZWrFiB1NTUgNsRDNG8Gq21jDgVABsuWwVckRDdfSUk2LoN3ZqaGqSlpfkumP373/8GazVdiOM4HD16NOAGZGVloayszPf1kiVLIAg/XBlfs2YNcnJy2nyPKIpYsGABli5divz8fLz66qtYuXJl2LeY9O4wFvUj3Va377kigZaBE9IT3Q5TvPvnHjt2DMeOHUNBQYHvv48dO9arwG3P6XRi69atuPvuu7s878iRI9BoNMjPzwcAFBUVYfv27UFrR6C8I91ovpAGAHFKHgY1T8uBCQlAtyNd1m4S/N69e0PWmF27diE9PR3Dhg3zPTZ//nwwxjBy5EjMmzcPiYmJqKqqQkZGhu+clJQUiKKIhoYGGAwGv18vNTXe73P9uRLv3Xch2rZ17IgxTokL3cxgAPx73yJBtPQjELHcdyD4/e82dNtv4dc+hIPp3XffbTPKLS0thdFohNPpxJIlS1BSUoKVK1cG7fVMJgtEsfv+pKUltLm7bWcfgsMdvZvdtJcZr8LRBicanQKS1J3vOBYNdwVu//nHkljuO9B9/wMJ5G5DVxAEfP31176wbf81APz0pz/t8Qu3d/nyZezduxfLly/3PWY0GgEAarUaM2fOxC9/+Uvf45WVlb7z6urqwPN8j0a5oWDzlReif6SbpffUci9YXEhKoW0eCfFXt6GbmpqKZ555xve1wWBo8zXHcfj000973ZD33nsPN998M5KTkwEAVqsVgiAgISEBjDFs27YNubm5AIDhw4fDbrdj3759yM/Px6ZNmzB58uRet6G37G4GBQeooj9z0VengJrncKHZjWEp4W4NIZGj29DdtWuXFO3Ae++9h2effdb3tclkwty5cyEIAkRRxKBBg1BcXAwA4Hkey5cvR3FxcZspY+FmEzyr0WLhrgo8xyFTr0SFhW7fQ0hPyGbx/I4dO9p83b9/f7z//vudnn/99ddj69atIW5Vz9gFFhMX0bz6x6vwRZUVNrcYU/0mpDfoNyWIbO7o3Uu3I766Lt2skhC/UegGkV0QoY2hEZ9RrwTPARcsnU8d8258E+vTjgjxip2EkIA9xka6Kp5rma9LI11C/EWhG0SeHcZiJ3QBT4mhyuqGy4/5zoQQCt2gEUQGp8hiqrwAAFnxSogM3e6vSwjxiK2ECCG7dwlwDI50AeC8H1PHqLZLCIVu0PhWo8XYSFen5NFPp8RZszPcTSEkIsRWQoSQ3R2bI10AyE5QobLZTZuaE+IHCt0e6OrP4x9GujEYuokqiAAqupg6RgjxoNANErtvh7HYe0uz9CooOVCJgRA/xF5ChEisXkgDACXPIStehbNmmq9LSHcodIPEW17QxGDoAp66bq1dgMVFdV1CukKhGyTe1WixsMNYR7IT1ACoxEBIdyh0g8Szw1hsBi4ApOsU0Ck4v0oMtB8DiWUUukFic4sxeRHNi+M4DEjw1HVDeUsnQiJd7KZEkMXivgvtZSeoYXGJMDmEcDeFENmi0A2SWNvWsSPZCZ4lwWebaBYDIZ2J7ZQIoljbwLwjBo0CBjVPU8cI6QKFbhAwxmAXWEyuRmsvO0GN8xYXBKrrEtIhCt0gcPgWRtDbmZ2oglNktNUjIZ2glAgCm+BdAkwj3QHxLXVdKjEQ0iEK3SCwuz2rsOiOuC1bPcbRVo+EdIZSIghopNvWQNrqkZBOUegGgW+HMbqQBgAYkEBbPRLSGQrdILC3jOjoQppHll4FFQ+cbqISAyHtUUoEAZUX2lLyHK6IV+F0k5OWBBPSDoVuENjdItQ8BwVPoes1KFGNBqeIegfVdQlpjUI3CGwCrUZr78pEz1aP5VRiIKQNCt0gsLtpNVp7Bo0CKRpFt3Vd2uKRxBoK3SCwCyJdROvAlYkqnLe44BKprkuIFyVFENho34UODUpUQ2DAOVqdRogPhW4Q2N0i1XQ70D+epo4R0p4y3A2IdN4dxqi88GPeqWPlLVPHurp/XOu6bk2NWYrmERIWlBS95BIBgdFqtM7kJGnQ6BRx2UZ3kyAEkEnoFhQUYPLkyZg2bRqmTZuGL7/8EgBw8OBBTJ06FYWFhbj//vthMpl839PVMSnZaDValwYb1OAAHG9whLsphMiCbJJizZo1KCsrQ1lZGcaNGwdRFLFgwQIsWrQIO3bsQH5+PlauXAkAXR6Tmo32XehSnJLHgAQVjjU4aHUaIZBR6LZ35MgRaDQa5OfnAwCKioqwffv2bo9JzbfvAm3r2KmrDRrUO0RUU4mBEPlcSJs/fz4YYxg5ciTmzZuHqqoqZGRk+I6npKRAFEU0NDR0ecxgMPj9mqmp8X6f29kEfu9IN9bvBNyVwUlq7KjwlBjS47r/kZPjYgk5tkkqsdx3IPj9l0XolpaWwmg0wul0YsmSJSgpKcEtt9wS8tc1mSwQ/Zi4n5aWgJoac4dvvo02MO+WXsXjingVjjU4Mc4Y1+UsBkB+sxe8n38siuW+A933P5BAlkVSGI1GAIBarcbMmTOxf/9+GI1GVFZW+s6pq6sDz/MwGAxdHpOaTaCRrj+uTlajziGgxt59iYGWBpNoFvbQtVqtMJs9/5IwxrBt2zbk5uZi+PDhsNvt2LdvHwBg06ZNmDx5MgB0eUxqtMOYf3KSNOAAHKunWQwktoW9vGAymTB37lwIggBRFDFo0CAUFxeD53ksX74cxcXFcDgcyMzMxIoVKwCgy2NSox3G/KNX8chOUOG7eodfJQZCohXHYngeTzBquv8ob4TFJWLO1cmhaGJUOVJnxwfnLJg1OAlZLXcN7opcaomxXNeM5b4DUVzTjWR2gdFFND/lJGmg4oHvqMRAYhilRS/Z3FRe8JdawWFwkgZH6x0QaLtHEqModHvJ5hZppNsDw5I1sAuM7ihBYhalRS/8sMMYjXT9NTBRhTglh++pxEBiFIVuLzgEBgZASyNdv/Ech9xkDU42On1LqAmJJZQWvUALIwIzPFkDgQHHG6jEQGIPhW43vFNCaAlw8PSLUyJZw+O7OioxkNhDadELvs1uaFvHHuE4DsOStThvcaHJ2fmyYO9yYFoSTKIJhW4v0AbmgRueogEAHDLZw9wSQqRFadELNNINnEGjwJWJKhwyOSDG7qJIEoModHvBO9LV0IW0gOT10cLiEnGqsfsLalRmINGCQrcX7C2r0XjavCUggxLVSFDxOFjrf4mB6rwk0lHo9oJnNRoFbqB4jsN1qVqcNrvQ4KBb+ZDYQKHbCzaB0UW0Xrou1bPPbk9Gu4REMkqMXrC7GY10eylBrcDgJDUOmuxwCnRBjUQ/Ct1esAkitDTS7bUx6TrYBYbDdTTaJdGPEqMXbDTSDYpMvQpZeiX2Vtto+hiJehS6ARJEBqdIG5gHy+i+OjQ6RdqPgUQ9SowA2Vvqj7SBeXAMTlIjRaPAnss2+HsHKZo6RiIRhW6AaLOb4OI4DqP76nDJ5sZ5iyvczSEkZCgxAkTbOgbf8BQN9EoOe6pt4W4KISFDoRsgGukGn5LnMDJNh9NNLtTY3OFuDiEhQYkRIBvVdEMir48WKh402iVRi0I3QHbfSJdCN5h0Sh7XpmrxfZ2jy712CYlUFLoBsrkZeA5Q8xS6wTYqTQcGYF8NLZYg0YdCN0BWQYROwYGjHcaCzqBRYIhBjcMmO9wiLZYg0YVCN0BWF4NeRW9fqIzoo4VdYDjeQPdRI9GFUiNAVreIOJq5EDID4lUwqHkcpNv5kChDqREgq1uEnkI3ZLiWvXYrLG6Y7DR9jEQPSo0AWWmzm5C7JlULHsAhU9clBrqbBIkkFLoBcLVsdkMj3dCKV/EYbFDjW7qgRqIIpUYArC1zdOPoQlrIjUjVwiYwnPTj5pWERAJKjQBYXS2hS+WFkMtOUCFJ3bObVxIiZxS6AWh2e/7UpfJC6HkvqJ2zuFBnpxVqJPIpw92A+vp6PPXUUzh//jzUajUGDBiAkpISpKSkYMiQIcjJyQHPe8Jt+fLlGDJkCABg165dWL58OQRBwLBhw7B06VLodDpJ2uwrL1DoSuKaVA2+rLJif60NE7Piuzy3/cW0mhpzKJtGSI+FPTU4jsODDz6IHTt2YOvWrejfvz9WrlzpO75p0yaUlZWhrKzMF7jNzc147rnn8Prrr+Pjjz+GXq/H+vXrJWszha60ElQKDE3W4JDJ7nvvCYlUYU8Ng8GAMWPG+L4eMWIEKisru/yeL774AsOHD0d2djYAoKioCB999FEom9mG1c2g5AC6jiadn6Tr4BKB/9TQ7mMksoW9vNCaKIrYuHEjCgoKfI/Nnj0bgiDgpptuwty5c6FWq1FVVYWMjAzfORkZGaiqqurx66Wmdv2nameaXSLiVDztuyChNJ0SOUlq7K22I6+PDvF+/osXrLm7sTwHOJb7DgS//7IK3cWLFyMuLg6zZs0CAHz++ecwGo2wWCxYsGAB1q1bh9/85jdBez2TyQKxm/mfHb3htBotPH6Wocebx+rxZVUzplzR81+EQOu7aWkJMVsbjuW+A933P5BAlk1yLFu2DOfOncNLL73ku3BmNBoBAPHx8Zg+fTr279/ve7x1CaKystJ3rhQ8+y7QKFdqKVoFRvbR4pDJQfdRIxFLFqG7atUqHDlyBOvWrYNarQYANDY2wm73zM10u93YsWMHcnNzAQDjxo3Dt99+i7NnzwLwXGybMmWKZO21uES//7wlwTXOqIdBzePDc2Y4BLqoRiJP2MsLJ0+exBtvvIHs7GwUFRUBALKysvDggw9i0aJF4DgObrcbeXl5eOKJJwB4Rr4lJSV45JFHIIoicnNz8eyzz0rSXpExNLsZhW6YqBUc7shOwDsnGvHphWbcOiC2640k8oQ9dAcPHozjx493eGzr1q2dft/EiRMxceLEUDWr01pNc8tqNArd8MnUq/CTdB3+77INVyWpkWPQhLtJhPiNkqOHLBS6snBjvzik6xT46LwFZj/vpUY7kRE5oOToIXNL6CaoFGFuSWxT8BymZSfCzRi2nDNDZLQLGYkMFLo9RCNd+UjRKjApKx4VFjf+77L/iyZa779Lo18iNUqOHrK4RHCgHcbkYniKBsOSNfiqyooKmkZGIgCFbg9ZXCL0Kh48rUaTBY7jMKm/HklqHlvPmmGjvRmIzFHo9hDN0ZUfjYLHtIEJsLhFfHTeAkb1XSJjlB49ZKbQlSVjnAo3G+NwotGJA7ThOZExSo8esrhFJFDoytLovjpcmaDCpxebUW2jOwgTeaL06AFBZLDRajTZ4jgOtw1IgFbB4b0zTb6FLITICaVHDzT55ujS2yZXehWPu65MhNkp4u/ljbQ/A5EdSo8eaHR4Vj4lqeltk7NMvQp3DkxEtU3A5tNmv27fTnN3iVQoPXqgwekZNRk0tBpN7gYlqXHbgHics7jw3pkmv4KXEClQ6PZAo1MADyovRIrhKVpM7h+P8qbAgpdGvCQUKD16oNEpIlFNCyMiyYg+WhT216O8yYX3z5gh0IiXhBmFbg80OAQkqam0EGny+ugwKUuPU01OvH82sOClUS8JFgrdHmh0CjBo6C2LRNen6XBLlh4nG1uCl1atkTChBPGTS/TcMYJGupFrZJoOE1uCt+wMBS8JDwpdP9F0seiQn6bDhEw9TjQ6sfl0E1wBXlyjUgMJFCWIn2i6WPQY1ddT4y1vcuFvpxppZzIiKQpdP5nsnrX8qRS6UeH6NB3+KzsBVVY3Sk82osnPW/4Q0lthvzFlpKi2CUhQ8dAq6d+paHF1sgY6JYd3T5vxzolG3HNVIvpo/f+V6KjEUFNjDmYTSRSiBPFTjd2NNC2NcqPNgAQ1Zg5OgsAY3jnRiIvNdPcJEloUuq10dpFEZAwmu4A+OvrDIBr1i1Nido4BWgWHv51q6tVtf2gPB9IdCl0/1DsECAw00o1iBo0C9+YkIV7F4+/ljThvphEvCQ0auvmhxua5yJJGI92olqBSYObgJGw81Yi/lzdianYCcgyaXj9vV6NdqgHHHhrp+qHG7gYHIJVGulEvXsVj5lVJSNMpsfmMGV9VWSHSIgoSRDR080N2ghpKjoOKp41uYoFexWPm4CRsP2/BV5esONnoQEGmHlfEq8AFebMjf2u+NCKOHhS6fugfr0L/eFW4m0EkpOI53JGdgKuS1Nh1sRkbTzWhX5wSo9K0GGLQQCnxP8DdhTOFcuSg0CWkC7nJGlyVpMaROjv2Vtux9ZwFH19oxtUGDYalaJClVwZ99BsIGjFHDo6x2C1YmUwWiK3W3tP0HtIVxhjOml04UufAiUYHXKJnL46rktS4KlGN/vEqyUfA4RBLwZ2WltBlfwPJDBrpEuInjuMwMFGNgYlqOAQ9TjQ4cazBgUO1dvynxg41zyE7QYVBSWoMSlRH7V2jk5J0UKuVcDrdUKs9EeL970h/rLnZAavVGdL3j0KXkABoFDyuSdXimlQtXCLDObML5U1OnGp04kSj55fWGKfEoEQ1rkpSI12nkEUZgoQfhS4hvaTiOU+JIUmNSVkM1TbBF8BfXbLiq0tWxCt5XJmkQh+tEokqHvEqHjolB52Sh1bB0S2gYgiFLiFBxHEc0uOUSI9TYmy/OFhdIsqbnChvcuJ4gxOHBUeH36dTcIhT8UhQ8UhU8UhQ80hUK5Co4pGo5mFQK6CIgXpxLIjo0D1z5gwWLlyIhoYGGAwGLFu2DNnZ2eFuFiE+caofyhCMMdgFBrNLhMUlwuYWYROY5//dDBaXCLNLxGmbC5Z2e/zy8CxVTtUqkKJRIFHtCeNElee/tQouoPIFY4zKHhKL6NAtLi7GzJkzMW3aNJSVlWHRokV4++23w90sQjrEcZyvpNBX1/W5gugJZ7NLRKNTQJ1dQK1dgMnhKV20v+EFzwEankP/eBXSdAqkaZVI1iigUXBQKziAAW7GIDDALXr+3+EW4HAKUKmVUPAcVByg4DkoOQ5KHi3/7/lvKn8ET8ROGTOZTCgsLMSePXugUCggCALGjBmDnTt3IiUlxa/nqK9vbjNlLDU1PlTNJSRoGGN477QZbsZgcghIVPFocnlGxklqBRpDtCG7kgMUHAcFBygVHHh4wthb9eDBgQODkufAcQDHPOdxjEGh4FvaLkLB8RAZ802vY4yB53kwUQTP8+A4QBREKBU8RNHzvRwAURShUPAQBdH3fK0f472Ptfx3Tx9TcBwGaDhwLrevz6mp8TCZLJ2+J4FkRsSOdKuqqpCeng6FwrMfgkKhQN++fVFVVeV36CYn60PZREJCguM43DUoMdzNiBnBHoxF50RCQgiRqYgNXaPRiMuXL0MQPH9KCYKA6upqGI3GMLeMEEI6F7Ghm5qaitzcXHzwwQcAgA8++AC5ubl+lxYIISQcIvZCGgCUl5dj4cKFaGpqQmJiIpYtW4Yrr7wy3M0ihJBORXToEkJIpInY8gIhhEQiCl1CCJEQhS4hhEiIQpcQQiREoduNM2fOYMaMGSgsLMSMGTNw9uzZcDepV+rr6/HQQw+hsLAQd9xxBx577DHU1dUBAA4ePIipU6eisLAQ999/P0wmk+/7Aj0mV6+88gqGDBmCEydOAIidvjscDhQXF2PSpEm444478NxzzwHo+uc80GNy9Nlnn+G//uu/MG3aNEydOhU7d+4EIHH/GenS7Nmz2fvvv88YY+z9999ns2fPDnOLeqe+vp59/fXXvq9ffPFF9rvf/Y4JgsAmTpzI9u7dyxhjbN26dWzhwoWMMRbwMbk6cuQIe+CBB9j48ePZ8ePHY6rvixcvZkuWLGGiKDLGGKupqWGMdf1zHugxuRFFkeXn57Pjx48zxhg7evQoGzFiBBMEQdL+U+h2oba2lo0cOZK53W7GGGNut5uNHDmSmUymMLcseLZv385+8YtfsEOHDrHbbrvN97jJZGIjRoxgjLGAj8mRw+Fg99xzD6uoqPCFbqz03WKxsJEjRzKLxdLm8a5+zgM9JkeiKLLRo0ezffv2McYY++abb9ikSZMk73/EbngjhWBsqiNnoihi48aNKCgoQFVVFTIyMnzHUlJSIIoiGhoaAj5mMBik7I5fXn75ZUydOhVZWVm+x2Kl7xUVFTAYDHjllVewZ88e6PV6PPHEE9BqtZ3+nDPGAjomx98PjuPw0ksv4Ve/+hXi4uLQ3NyMP//5z13+noei/1TTjWGLFy9GXFwcZs2aFe6mSOLAgQM4cuQIZs6cGe6mhIUgCKioqMDQoUOxefNmzJ8/H3PnzoXVag130yThdrvxxhtv4NVXX8Vnn32G1157DU8++aTk/aeRbhdab6rj3bM3WjbVWbZsGc6dO4fXX38dPM/DaDSisrLSd7yurg48z8NgMAR8TG727t2L8vJyTJgwAQBw6dIlPPDAA5g9e3bU9x3w/DwrlUrcfvvtAIDrrrsOycnJ0Gq1nf6cM8YCOiZHR48eRXV1NUaOHAkAGDlyJHQ6HTQajaT9p5FuF6J1U51Vq1bhyJEjWLduHdRqNQBg+PDhsNvt2LdvHwBg06ZNmDx5cq+Oyc3DDz+Mr776Crt27cKuXbvQr18/rF+/Hg8++GDU9x3wlD/GjBmDf/3rXwA8V95NJhOys7M7/Tnv6ncg0n4/+vXrh0uXLuH06dMAPHu3mEwmDBgwQNL+094L3Yi2TXVOnjyJ22+/HdnZ2dBqtQCArKwsrFu3Dvv370dxcTEcDgcyMzOxYsUK9OnTBwACPiZnBQUFeP3115GTkxMzfa+oqMAzzzyDhoYGKJVKPPnkk7j55pu7/DkP9JgcbdmyBf/7v//ruy/c448/jokTJ0rafwpdQgiREJUXCCFEQhS6hBAiIQpdQgiREIUuIYRIiEKXEEIkRKEbYdauXYv58+eHuxk9tnDhQqxevTooz/XXv/4VY8eORV5eHurr64PynIHYt28fCgsLw/b6t912G/bs2ROU59qyZQvuv/9+39dDhgzBuXPngvLcAJCXl4eKioqgPV8ko9CVoa1bt+Kuu+5CXl4ebrzxRjz44IO+CfjBcOHCBQwZMgRutztoz9kbmzdvxs9//nO/znW5XHjxxRfxl7/8BQcOHEBycnJAr7lnzx7cdNNNPfqe9kGUn5+PHTt2BPT6XfF+Pnl5ecjLy8PYsWPxyCOP+BY1eH344YcYM2aMX8/V3Wc9depU/OUvf+l12wFg9uzZ+Mc//tHmsQMHDqB///5Bef5IR6ErMxs2bMALL7yARx99FP/617/w2WefYebMmfj000/D3TSfcIa1yWSCw+HAVVddFbY2SGXv3r04cOAAysrKMHbsWDz22GPYvHlz0F9HLv/4xgoKXRkxm81Ys2YNFi1ahEmTJiEuLg4qlQoFBQV4+umnf3R+R6O1goIC/Pvf/wYAHD58GHfddReuv/56jB07FkuXLgUA3wY3o0aNQl5eHg4cOAAA+Oc//4kpU6Zg1KhReOCBB3Dx4kXf8w4ZMgSlpaWYNGkSJk2aBMCzIfS0adOQn5+PoqIiHDt2zHf+999/jzvvvBN5eXl48skn4XA4/H4fCgoKsH79etxxxx0YOXKk7/vPnDnjW2Y7atQo3HfffQA8q4LmzJmD0aNHo7CwENu2bfM91+7du3HrrbciLy8P48aNw/r162G1WvHQQw+hurraN5q8fPkyDh8+jBkzZiA/Px833ngjSkpK4HQ6AQD33nsvAGDatGnIy8vDtm3bfvT+l5eXY/bs2cjPz8dtt93W5h/KhQsX4o9//CMefvhh5OXlYfr06Th//rxf70daWhp+8Ytf4LHHHsPKlSshiqLvfQrks968eTOKiorwwgsvYMyYMVi7dm2Hf23s3r0bEyZMwJgxY7Bs2TLf67YvcbUeTa9evRr79u1DSUkJ8vLyUFJSAqDtXwlmsxlPPfUUfvKTn2D8+PF49dVXfc/tbceyZcswatQoFBQUYPfu3X69TxEjyFtWkl7YvXs3y83NZS6Xq9Nz1qxZw377298yxhj7+uuv2bhx49ocHz9+PPvXv/7FGGPsnnvuYe+99x5jzLOX6oEDBxhjjFVUVLCcnJw2r/Pxxx+ziRMnslOnTjGXy8XWrVvHZsyY4Tuek5PD/ud//ofV19czm83GvvvuO/aTn/yEHTx4kLndbrZ582Y2fvx45nA4mMPhYD/72c/Yhg0bmNPpZB999BEbOnQoW7VqVYd9evfdd1lRUVGbPtx9993s0qVLrL6+nk2ePJn99a9/7bDtzc3N7KabbmL//Oc/mcvlYt999x0bPXo0O3nyJGOMsRtuuMG3yXhDQwM7cuRIp+/dt99+yw4cOMBcLherqKhgkydPZhs2bGjzHpw9e9b3devncDqdbOLEiey1115jDoeD/fvf/2YjRoxg5eXljDHGnn76aTZ69Gh26NAh5nK52Lx589iTTz7Z4fvR0efDGGPnz59nOTk57NSpU773KZDP+t1332W5ubns7bffZi6Xi9lsth99Bjk5OWzWrFmsvr6eXbx4kU2aNIn9/e9/Z4y1/Rns6DVmzZrlO7ej927BggXs0UcfZWazmVVUVLR57nfffZcNHTqU/e1vf2Nut5uVlpayG264wbfpejSgka6MNDQ0IDk5GUplcDZ/UyqVOH/+POrq6qDX6zFixIhOz920aRMefvhhDBo0CEqlEo8++iiOHj3aZrT78MMPw2AwQKvV4m9/+xtmzJiB6667DgqFAnfeeSdUKhUOHjyIQ4cOweVy4Re/+AVUKhUmT56Ma665pkdtnz17NtLT02EwGDB+/HgcPXq0w/M+//xzZGZm4u6774ZSqcTQoUNRWFiI7du3+96DU6dOwWKxICkpCcOGDev0NYcPH44RI0ZAqVQiKysLM2bMwN69e/1q76FDh2C1WvHwww9DrVbjpz/9KcaPH48PP/zQd87EiRNx7bXXQqlUYurUqZ32qTN9+/YF4Pk5aa8nn7X3uWbPng2lUunbg6O9hx56CAaDARkZGbjvvvt8G7v0hiAI2LZtG377298iPj4eWVlZmDNnDrZs2eI7JyMjA/fcc4/v56qmpga1tbW9fm25oNCVEYPBgPr6+qDV2JYsWYKzZ89iypQpuPvuu/HZZ591em5lZSVeeOEF5OfnIz8/H6NHj/ZtXefVesu6yspKbNiwwXd+fn4+Ll26hOrqalRXVyM9Pd23qQiANpt9+yMtLc333zqdrtM9Ty9evIjDhw+3acfWrVtRU1MDAFizZg12796N8ePHY9asWb5SSkfOnDmDRx55BDfccAOuv/56rF692u/ZEdXV1ejXrx94/odfqYyMjDbvX+uNcLRabY/3cfU+V0dbR/bkswY8O251p/XnnZmZierq6h61tyP19fVwuVxtfh66ep90Oh0ARNWev7Sfrozk5eVBrVbjk08+8WuLQJ1OB7vd7vtaEATfTSYBIDs7G6tWrYIoiti5cycef/xx7Nmzp00YehmNRjz66KOYOnVqp6/X+vu85//yl7/80XnffPMNLl++DMaY73sqKytDcvXaaDRi1KhR2LBhQ4fHr732Wrz22mtwuVwoLS3Fk08+id27d3f4HvzhD3/A0KFD8ac//Qnx8fF46623/J6d0LdvX1y6dAmiKPqCt6qqCtnZ2QH3rb2PP/4YqampGDhw4I+O9eSzBtDp461VVVVh8ODBADyfn3ek3f7nriej0OTkZKhUKlRWVvouhnrv3BAraKQrIwkJCXj88cdRUlKCTz75BDabDS6XC7t378by5ct/dP7AgQPhcDjw+eefw+Vy4bXXXvNd+AGAsrIy38baiYmJAACe55GSkgKe59vMmywqKsKf//xnnDx5EoDnYsdHH33UaVunT5+OTZs24dChQ2CMwWq14vPPP4fFYvH9if7222/D5XJh586d+Pbbb4P1NrXxs5/9DGfPnsX7778Pl8sFl8uFw4cPo7y8HE6nE1u2bIHZbIZKpYJer/cFYmpqKhoaGmA2m33P1dzcDL1eD71ej/LycmzcuLHNa/Xp06fTuabXXnsttFot3nzzTbhcLuzZswe7du3Crbfe2us+1tbW4p133sErr7yCefPmtRlNe/Xks/bX+vXr0djYiKqqKrz99tu+vuTm5mLv3r2orKyE2WzGG2+80eb7unqfFAoFJk+ejNWrV8NiseDixYvYsGFDl//YRxsa6crM/fffjz59+uDVV1/F/PnzodfrMWzYMDz66KM/OjchIQHFxcX4/e9/D0EQ8OCDD7b5s/HLL7/Eiy++CLvdjoyMDKxevdpXv3v00Ufx85//HG63G2+++SZuueUWNDc3Y968ebh48SISEhIwduxYTJkypcN2XnPNNVi8eDFKSkpw7tw5aLVaXH/99cjPz4darcbatWvx3HPP4aWXXsLNN9+MW265JSTvV3x8PNavX48XX3wRL774IhhjGDJkCH73u98B8ITR4sWLIQgCBg4ciBUrVgAABg0ahNtuuw0TJ06EIAj48MMP8fTTT+O5557D+vXrkZubi1tvvRVff/2177Uee+wxLFy4EHa7HSUlJUhNTfUdU6vVeP311/HHP/4Rb7zxBtLT07F8+XIMGjQo4L6NGjUKjDHodDoMHz4cL7/8cqdzi3vyWftrwoQJuOuuu2CxWHDnnXfiv//7vwEAN9xwA2699VZMnToVycnJeOihh7Br1y7f9913331YuHAhNm7ciGnTpuH3v/99m+d97rnnsHjxYkycOBEajQbTp0/H3Xff3dO3J2LRfrqEECIhKi8QQoiEKHQJIURCFLqEECIhCl1CCJEQhS4hhEiIQpcQQiREoUsIIRKi0CWEEAlR6BJCiIT+Hw9nTG4QNsdhAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "ax = sns.displot(infest_clust,\n", " kde=True,\n", " bins=100,\n", " color='skyblue',\n", " linewidth=15,\n", " alpha=1)\n", "ax.set(xlabel='Clustered Infestation Distribution', ylabel='Frequency')" ] }, { "cell_type": "code", "execution_count": 11, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "764\n" ] } ], "source": [ "num_boxes = math.ceil(num_stem / 200)\n", "print(num_boxes)" ] } ], "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.8.5" } }, "nbformat": 4, "nbformat_minor": 4 }